
&gt;Dear all,<br>&gt;I have a problem in calculating phonons of a material.I think the ecutwfc = 44.0982(600eV) <br>&gt;and the K_POINTS is enough(in vc-relax I set it 10 10 8). The forces of atoms is also list <br>&gt;below. But the frequencies in *.dyn1 is negative(about -1.67895THz). <br>&amp;CONTROL<br>&nbsp;&nbsp; calculation =&nbsp;&nbsp; 'scf' ,<br>&nbsp;&nbsp; restart_mode = 'from_scratch' ,<br>&nbsp;&nbsp; outdir='/home/user/ecnuphyk/SUN/tmp/tmp-p4/' ,<br>&nbsp;&nbsp; pseudo_dir = '/home/user/ecnuphyk/SUN/pseudo/' ,<br>&nbsp;&nbsp; prefix='...' ,<br>&nbsp;&nbsp; tstress = .true. ,<br>&nbsp;&nbsp; tprnfor = .true. ,<br>&nbsp;&nbsp; etot_conv_thr = 1.0D-5 ,<br>&nbsp;&nbsp; forc_conv_thr = 1.0D-4 ,<br>&nbsp;/<br>&nbsp;&amp;SYSTEM<br>&nbsp;&nbsp; ibrav = 0 ,<br>&nbsp;&nbsp; celldm(1) = ...&nbsp; ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nat&nbsp; =&nbsp; 16

 ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ntyp&nbsp; =&nbsp; 2 ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nspin&nbsp; =&nbsp; 1 ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ecutwfc&nbsp; = 44.0982 ,#about 600eV<br>&nbsp;/<br>&nbsp;&amp;ELECTRONS<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; conv_thr&nbsp; =&nbsp; 1.0d-7&nbsp; ,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; mixing_beta = 0.7 ,<br>&nbsp;/<br>CELL_PARAMETERS (alat)<br>...<br>ATOMIC_SPECIES<br>...<br>ATOMIC_POSITIONS (crystal)<br>...<br>K_POINTS automatic<br>&nbsp;8 8 6 0 0 0<br>*************The forces of

 the atoms in the self-consistent calculation:<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp; 1 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00000085&nbsp;&nbsp; -0.00000187&nbsp;&nbsp; -0.00000180<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp; 2 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp; -0.00000200&nbsp;&nbsp;&nbsp; 0.00000092&nbsp;&nbsp;&nbsp; 0.00000220<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp; 3 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00000200&nbsp;&nbsp; -0.00000092&nbsp;&nbsp; -0.00000220<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp; 4 type&nbsp; 1&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp; -0.00000085&nbsp;&nbsp;&nbsp; 0.00000187&nbsp;&nbsp;&nbsp; 0.00000180<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp; 5 type&nbsp; 2&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp; -0.00002835&nbsp;&nbsp;&nbsp; 0.00003010&nbsp;&nbsp; -0.00000078<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp; 6 type&nbsp; 2&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;

 -0.00002239&nbsp;&nbsp; -0.00001414&nbsp;&nbsp;&nbsp; 0.00000083<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp; 7 type&nbsp; 2&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00001800&nbsp;&nbsp;&nbsp; 0.00001850&nbsp;&nbsp;&nbsp; 0.00000131<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp; 8 type&nbsp; 2&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00000012&nbsp;&nbsp; -0.00000010&nbsp;&nbsp; -0.00000065<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp;&nbsp; 9 type&nbsp; 2&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00002757&nbsp;&nbsp; -0.00002933&nbsp;&nbsp; -0.00000124<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 10 type&nbsp; 2&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp; -0.00000019&nbsp;&nbsp; -0.00000006&nbsp;&nbsp;&nbsp; 0.00000066<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 11 type&nbsp; 2&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00000019&nbsp;&nbsp;&nbsp; 0.00000006&nbsp;&nbsp; -0.00000066<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 12 type&nbsp; 2&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;

 -0.00002757&nbsp;&nbsp;&nbsp; 0.00002933&nbsp;&nbsp;&nbsp; 0.00000124<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 13 type&nbsp; 2&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp; -0.00000012&nbsp;&nbsp;&nbsp; 0.00000010&nbsp;&nbsp;&nbsp; 0.00000065<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 14 type&nbsp; 2&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp; -0.00001800&nbsp;&nbsp; -0.00001850&nbsp;&nbsp; -0.00000131<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 15 type&nbsp; 2&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00002239&nbsp;&nbsp;&nbsp; 0.00001414&nbsp;&nbsp; -0.00000083<br>&nbsp;&nbsp;&nbsp;&nbsp; atom&nbsp; 16 type&nbsp; 2&nbsp;&nbsp; force =&nbsp;&nbsp;&nbsp;&nbsp; 0.00002835&nbsp;&nbsp; -0.00003010&nbsp;&nbsp;&nbsp; 0.00000078<br>********************ph.in<br>phonons of ...<br>&amp;inputph<br>&nbsp; tr2_ph=1.0d-12,<br>&nbsp; prefix='...',<br>&nbsp; ldisp=.true.,<br>&nbsp; nq1=4, nq2=4, nq3=3<br>&nbsp; epsil=.true. ,<br>&nbsp; ...<br>&nbsp; outdir= '/home/user/ecnuphyk/SUN/tmp/tmp-p4/' ,<br>&nbsp;

 fildyn='....dyn',<br>&nbsp;/<br>&gt;*********************************<BR><BR><p>&#32;


      <hr size=1><a href="https://member.cn.yahoo.com/cnreg/reginfo.html?id=89034" target=blank>雅虎邮箱，终生伙伴！</a>