>Dear all,<br>>I have a problem in calculating phonons of a material.I think the ecutwfc = 44.0982(600eV) <br>>and the K_POINTS is enough(in vc-relax I set it 10 10 8). The forces of atoms is also list <br>>below. But the frequencies in *.dyn1 is negative(about -1.67895THz). <br>&CONTROL<br>   calculation =   'scf' ,<br>   restart_mode = 'from_scratch' ,<br>   outdir='/home/user/ecnuphyk/SUN/tmp/tmp-p4/' ,<br>   pseudo_dir = '/home/user/ecnuphyk/SUN/pseudo/' ,<br>   prefix='...' ,<br>   tstress = .true. ,<br>   tprnfor = .true. ,<br>   etot_conv_thr = 1.0D-5 ,<br>   forc_conv_thr = 1.0D-4 ,<br> /<br> &SYSTEM<br>   ibrav = 0 ,<br>   celldm(1) = ...  ,<br>                         nat  =  16
 ,<br>                        ntyp  =  2 ,<br>                       nspin  =  1 ,<br>                     ecutwfc  = 44.0982 ,#about 600eV<br> /<br> &ELECTRONS<br>                    conv_thr  =  1.0d-7  ,<br>                  mixing_beta = 0.7 ,<br> /<br>CELL_PARAMETERS (alat)<br>...<br>ATOMIC_SPECIES<br>...<br>ATOMIC_POSITIONS (crystal)<br>...<br>K_POINTS automatic<br> 8 8 6 0 0 0<br>*************The forces of
 the atoms in the self-consistent calculation:<br>     atom   1 type  1   force =     0.00000085   -0.00000187   -0.00000180<br>     atom   2 type  1   force =    -0.00000200    0.00000092    0.00000220<br>     atom   3 type  1   force =     0.00000200   -0.00000092   -0.00000220<br>     atom   4 type  1   force =    -0.00000085    0.00000187    0.00000180<br>     atom   5 type  2   force =    -0.00002835    0.00003010   -0.00000078<br>     atom   6 type  2   force =   
 -0.00002239   -0.00001414    0.00000083<br>     atom   7 type  2   force =     0.00001800    0.00001850    0.00000131<br>     atom   8 type  2   force =     0.00000012   -0.00000010   -0.00000065<br>     atom   9 type  2   force =     0.00002757   -0.00002933   -0.00000124<br>     atom  10 type  2   force =    -0.00000019   -0.00000006    0.00000066<br>     atom  11 type  2   force =     0.00000019    0.00000006   -0.00000066<br>     atom  12 type  2   force =   
 -0.00002757    0.00002933    0.00000124<br>     atom  13 type  2   force =    -0.00000012    0.00000010    0.00000065<br>     atom  14 type  2   force =    -0.00001800   -0.00001850   -0.00000131<br>     atom  15 type  2   force =     0.00002239    0.00001414   -0.00000083<br>     atom  16 type  2   force =     0.00002835   -0.00003010    0.00000078<br>********************ph.in<br>phonons of ...<br>&inputph<br>  tr2_ph=1.0d-12,<br>  prefix='...',<br>  ldisp=.true.,<br>  nq1=4, nq2=4, nq3=3<br>  epsil=.true. ,<br>  ...<br>  outdir= '/home/user/ecnuphyk/SUN/tmp/tmp-p4/' ,<br> 
 fildyn='....dyn',<br> /<br>>*********************************<BR><BR><p>



      <hr size=1><a href="https://member.cn.yahoo.com/cnreg/reginfo.html?id=89034" target=blank>Ż䣬飡</a>