[Pw_forum] convergence problem in phonon calculation

杨顶峰 yangxunscience at 163.com
Wed May 13 15:04:13 CEST 2015


Dear All,


                I am trying the quantum espresso 5.1 to perform the phonon calculation of a compound.
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 .
my input file for the self consistent calculation  is given below:
&control
   title='sbcrse3'
   calculation = 'scf' ,
   restart_mode='from_scratch',
   wf_collect=.true.,
   prefix = 'sbcrse3',
   outdir = './',
   pseudo_dir = './',
   tstress=.true.
   tprnfor=.true.
   etot_conv_thr=1.0D-6,
   forc_conv_thr=1.0D-6,
/


&system
   ibrav = 8,
   celldm(1) = 17.721473993588537072269931508765,
   celldm(2) = 0.40987224478027270642574239356715,
   celldm(3) = 1.4354965598661794909645534975973,
   nat = 20 ,
   ntyp = 3 ,
   ecutwfc = 40,
   ecutrho = 300,
   occupations = 'fixed' ,
   nbnd = 150,
/


&electrons
   mixing_beta=0.7
   conv_thr = 1.0d-10 ,
/


ATOMIC_SPECIES
  Sb    121.7       Sb.pbe-mt_fhi.UPF
  Se    78.96       Se.pbe-mt_fhi.UPF
  Cr    51.996      Cr.pbe-mt_fhi.UPF
ATOMIC_POSITIONS crystal
Sb     0.471599996         0.250000000         0.655629992
Sb     0.528400004         0.750000000         0.344370008
Sb     0.028400004         0.750000000         0.155629992
Sb     0.971599996         0.250000000         0.844370008
.....


K_POINTS automatic
  13 13 13   0  0  0


the above calculation was running successfully. the details shows as follows:




convergence has been achieved in  76 iterations


     Forces acting on atoms (Ry/au):


     atom    1 type  1   force =     0.00330300   -0.00000173    0.01711914
     atom    2 type  1   force =    -0.00330300    0.00000173   -0.01711914
     atom    3 type  1   force =    -0.00330300    0.00000173    0.01711914
     atom    4 type  1   force =     0.00330300   -0.00000173   -0.01711914
     atom    5 type  2   force =    -0.00590819   -0.00028964    0.01684058
    ........


     Total force =     0.064539     Total SCF correction =     0.000006




     entering subroutine stress ...


          total   stress  (Ry/bohr**3)                   (kbar)     P=  -77.36
  -0.00044232   0.00000132   0.00000000        -65.07      0.19      0.00
   0.00000132  -0.00068747   0.00000000          0.19   -101.13      0.00
   0.00000000   0.00000000  -0.00044787          0.00      0.00    -65.88 


But when i begin the phonon calculation there is a problem, some tips will be show in the "out " file:
Electric Fields Calculation
     ik   1 ibnd  59 linter: root not converged  5.993E-09
     ik   1 ibnd  59 linter: root not converged  3.084E-09
     ik   2 ibnd  59 linter: root not converged  2.722E-07
     ik   2 ibnd  59 linter: root not converged  3.810E-09
     ik   2 ibnd  59 linter: root not converged  1.129E-09
     ik   3 ibnd  59 linter: root not converged  1.634E-09
     ik   3 ibnd  59 linter: root not converged  9.282E-09
     ik   3 ibnd  59 linter: root not converged  1.749E-09
     ik   4 ibnd  59 linter: root not converged  2.191E-09
     ik   4 ibnd  59 linter: root not converged  2.338E-09
     ik   4 ibnd  59 linter: root not converged  1.895E-09
     ik   5 ibnd  59 linter: root not converged  1.241E-08
     ik   5 ibnd  59 linter: root not converged  1.440E-09
     ik   5 ibnd  59 linter: root not converged  3.946E-09
     ik   6 ibnd  59 linter: root not converged  5.932E-08
     ik   6 ibnd  59 linter: root not converged  6.026E-09
     ik   6 ibnd  59 linter: root not converged  5.804E-08
     ik   7 ibnd  59 linter: root not converged  6.012E-09
     ik   7 ibnd  59 linter: root not converged  1.088E-09
     ik   7 ibnd  59 linter: root not converged  1.407E-08
     ik   8 ibnd  59 linter: root not converged  1.962E-09
................
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!
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20150513/ef65aa01/attachment.html>


More information about the users mailing list