<div style="line-height:1.7;color:#000000;font-size:14px;font-family:Arial"><div>Dear All,</div><div><br></div><div>                I am trying the quantum espresso 5.1 to perform the phonon calculation of a compound.</div><div>firstly, i use VASP to get a better geometry calculation , then the opted model was used to carried out the self consistent calculation and phonon calculation on quantum espresso 5.1 .</div><div>my input file for the self consistent calculation  is given below:</div><div><div><span style="color: rgb(0, 128, 0); ">&control</span></div><div><span style="color: rgb(0, 128, 0); ">   title='sbcrse3'</span></div><div><span style="color: rgb(0, 128, 0); ">   calculation = 'scf' ,</span></div><div><span style="color: rgb(0, 128, 0); ">   restart_mode='from_scratch',</span></div><div><span style="color: rgb(0, 128, 0); ">   wf_collect=.true.,</span></div><div><span style="color: rgb(0, 128, 0); ">   prefix = 'sbcrse3',</span></div><div><span style="color: rgb(0, 128, 0); ">   outdir = './',</span></div><div><span style="color: rgb(0, 128, 0); ">   pseudo_dir = './',</span></div><div><span style="color: rgb(0, 128, 0); ">   tstress=.true.</span></div><div><span style="color: rgb(0, 128, 0); ">   tprnfor=.true.</span></div><div><span style="color: rgb(0, 128, 0); ">   etot_conv_thr=1.0D-6,</span></div><div><span style="color: rgb(0, 128, 0); ">   forc_conv_thr=1.0D-6,</span></div><div><span style="color: rgb(0, 128, 0); ">/</span></div><div><br></div><div><span style="color: rgb(0, 128, 0); ">&system</span></div><div><span style="color: rgb(0, 128, 0); ">   ibrav = 8,</span></div><div><span style="color: rgb(0, 128, 0); ">   celldm(1) = 17.721473993588537072269931508765,</span></div><div><span style="color: rgb(0, 128, 0); ">   celldm(2) = 0.40987224478027270642574239356715,</span></div><div><span style="color: rgb(0, 128, 0); ">   celldm(3) = 1.4354965598661794909645534975973,</span></div><div><span style="color: rgb(0, 128, 0); ">   nat = 20 ,</span></div><div><span style="color: rgb(0, 128, 0); ">   ntyp = 3 ,</span></div><div><span style="color: rgb(0, 128, 0); ">   ecutwfc = 40,</span></div><div><span style="color: rgb(0, 128, 0); ">   ecutrho = 300,</span></div><div><span style="color: rgb(0, 128, 0); ">   occupations = 'fixed' ,</span></div><div><span style="color: rgb(0, 128, 0); ">   nbnd = 150,</span></div><div><span style="color: rgb(0, 128, 0); ">/</span></div><div><br></div><div><span style="color: rgb(0, 128, 0); ">&electrons</span></div><div><span style="color: rgb(0, 128, 0); ">   mixing_beta=0.7</span></div><div><span style="color: rgb(0, 128, 0); ">   conv_thr = 1.0d-10 ,</span></div><div><span style="color: rgb(0, 128, 0); ">/</span></div><div><br></div><div><span style="color: rgb(0, 128, 0); ">ATOMIC_SPECIES</span></div><div><span style="color: rgb(0, 128, 0); ">  Sb    121.7       Sb.pbe-mt_fhi.UPF</span></div><div><span style="color: rgb(0, 128, 0); ">  Se    78.96       Se.pbe-mt_fhi.UPF</span></div><div><span style="color: rgb(0, 128, 0); ">  Cr    51.996      Cr.pbe-mt_fhi.UPF</span></div><div><span style="color: rgb(0, 128, 0); ">ATOMIC_POSITIONS crystal</span></div><div><span style="color: rgb(0, 128, 0); ">Sb     0.471599996         0.250000000         0.655629992</span></div><div><span style="color: rgb(0, 128, 0); ">Sb     0.528400004         0.750000000         0.344370008</span></div><div><span style="color: rgb(0, 128, 0); ">Sb     0.028400004         0.750000000         0.155629992</span></div><div><span style="color: rgb(0, 128, 0); ">Sb     0.971599996         0.250000000         0.844370008</span></div><div><span style="color: rgb(0, 128, 0); ">.....</span></div><div><br></div><div><span style="color: rgb(0, 128, 0); ">K_POINTS automatic</span></div><div><span style="color: rgb(0, 128, 0); ">  13 13 13   0  0  0</span></div><div><br></div><div>the above calculation was running successfully. the details shows as follows:</div><div><br></div><div><br></div><div><span style="color: rgb(0, 128, 0); ">convergence has been achieved in  76 iterations</span></div><div><br></div><div><span style="color: rgb(0, 128, 0); ">     Forces acting on atoms (Ry/au):</span></div><div><br></div><div><span style="color: rgb(0, 128, 0); ">     atom    1 type  1   force =     0.00330300   -0.00000173    0.01711914</span></div><div><span style="color: rgb(0, 128, 0); ">     atom    2 type  1   force =    -0.00330300    0.00000173   -0.01711914</span></div><div><span style="color: rgb(0, 128, 0); ">     atom    3 type  1   force =    -0.00330300    0.00000173    0.01711914</span></div><div><span style="color: rgb(0, 128, 0); ">     atom    4 type  1   force =     0.00330300   -0.00000173   -0.01711914</span></div><div><span style="color: rgb(0, 128, 0); ">     atom    5 type  2   force =    -0.00590819   -0.00028964    0.01684058</span></div><div><span style="color: rgb(0, 128, 0); ">    ........</span></div><div><br></div><div><span style="color: rgb(0, 128, 0); ">     Total force =     0.064539     Total SCF correction =     0.000006</span></div><div><br></div><div><br></div><div><span style="color: rgb(0, 128, 0); ">     entering subroutine stress ...</span></div><div><br></div><div><span style="color: rgb(0, 128, 0); ">          total   stress  (Ry/bohr**3)                   (kbar)     P=  -77.36</span></div><div><span style="color: rgb(0, 128, 0); ">  -0.00044232   0.00000132   0.00000000        -65.07      0.19      0.00</span></div><div><span style="color: rgb(0, 128, 0); ">   0.00000132  -0.00068747   0.00000000          0.19   -101.13      0.00</span></div><div><span style="color: rgb(0, 128, 0); ">   0.00000000   0.00000000  -0.00044787          0.00      0.00    -65.88 </span></div><div><br></div><div><span style="font-size: 18px; ">But when i begin the phonon calculation there is a problem, some tips will be show in the "out " file:</span></div></div><div><div><span style="font-size: 18px; color: rgb(255, 0, 0); ">Electric Fields Calculation</span></div><div><span style="font-size: 18px; color: rgb(255, 0, 0); ">     ik   1 ibnd  59 linter: root not converged  5.993E-09</span></div><div><span style="font-size: 18px; color: rgb(255, 0, 0); ">     ik   1 ibnd  59 linter: root not converged  3.084E-09</span></div><div><span style="font-size: 18px; color: rgb(255, 0, 0); ">     ik   2 ibnd  59 linter: root not converged  2.722E-07</span></div><div><span style="font-size: 18px; color: rgb(255, 0, 0); ">     ik   2 ibnd  59 linter: root not converged  3.810E-09</span></div><div><span style="font-size: 18px; color: rgb(255, 0, 0); ">     ik   2 ibnd  59 linter: root not converged  1.129E-09</span></div><div><span style="font-size: 18px; color: rgb(255, 0, 0); ">     ik   3 ibnd  59 linter: root not converged  1.634E-09</span></div><div><span style="font-size: 18px; color: rgb(255, 0, 0); ">     ik   3 ibnd  59 linter: root not converged  9.282E-09</span></div><div><span style="font-size: 18px; color: rgb(255, 0, 0); ">     ik   3 ibnd  59 linter: root not converged  1.749E-09</span></div><div><span style="font-size: 18px; color: rgb(255, 0, 0); ">     ik   4 ibnd  59 linter: root not converged  2.191E-09</span></div><div><span style="font-size: 18px; color: rgb(255, 0, 0); ">     ik   4 ibnd  59 linter: root not converged  2.338E-09</span></div><div><span style="font-size: 18px; color: rgb(255, 0, 0); ">     ik   4 ibnd  59 linter: root not converged  1.895E-09</span></div><div><span style="font-size: 18px; color: rgb(255, 0, 0); ">     ik   5 ibnd  59 linter: root not converged  1.241E-08</span></div><div><span style="font-size: 18px; color: rgb(255, 0, 0); ">     ik   5 ibnd  59 linter: root not converged  1.440E-09</span></div><div><span style="font-size: 18px; color: rgb(255, 0, 0); ">     ik   5 ibnd  59 linter: root not converged  3.946E-09</span></div><div><span style="font-size: 18px; color: rgb(255, 0, 0); ">     ik   6 ibnd  59 linter: root not converged  5.932E-08</span></div><div><span style="font-size: 18px; color: rgb(255, 0, 0); ">     ik   6 ibnd  59 linter: root not converged  6.026E-09</span></div><div><span style="font-size: 18px; color: rgb(255, 0, 0); ">     ik   6 ibnd  59 linter: root not converged  5.804E-08</span></div><div><span style="font-size: 18px; color: rgb(255, 0, 0); ">     ik   7 ibnd  59 linter: root not converged  6.012E-09</span></div><div><span style="font-size: 18px; color: rgb(255, 0, 0); ">     ik   7 ibnd  59 linter: root not converged  1.088E-09</span></div><div><span style="font-size: 18px; color: rgb(255, 0, 0); ">     ik   7 ibnd  59 linter: root not converged  1.407E-08</span></div><div><span style="font-size: 18px; color: rgb(255, 0, 0); ">     ik   8 ibnd  59 linter: root not converged  1.962E-09</span></div></div><div><span style="font-size: 18px; color: rgb(255, 0, 0); ">................</span></div><div><span style="font-size: 18px; ">the calculation is still on progress, i really don`t know how to conquer the problem as a new beginner, please help me, Thanks a lot!</span></div></div><br><br><span title="neteasefooter"><span id="netease_mail_footer"></span></span>