<div dir="ltr"><div><div><div><div><div>Hi Jibiao,<br></div>Tuning of "alpha_mix(niter)" may help you in convergence, so you can try that.<br><br><br></div>Thanks,<br></div>Rita<br></div>NISER, Bhubaneswar<br></div>INDIA<br><div><div><div><div><span id="gmail-:1zg" tabindex="0" class="gmail-ams gmail-bkH gmail-amt"></span><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Oct 26, 2017 at 4:05 PM, jibiaoli <span dir="ltr"><<a href="mailto:jibiaoli@foxmail.com" target="_blank">jibiaoli@foxmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div><div>Dear all,</div><div><br></div><div>I try to calculate phonons at gamma for a fully relaxed structure(CO/Fe(100)), but the calculation failed to converge in 100 interations. Any suggestions to get better results? Thank you very much!</div><div><br></div><div><div>     Calculation of q =    0.0000000   0.0000000   0.0000000</div></div><div>...</div><div><div>     There are  58 irreducible representations</div></div><div>...</div><div><div>     Representation #  1 mode #   1</div><div><br></div><div>     Self-consistent Calculation</div><div> </div><div>     Pert. #  1: Fermi energy shift (Ry) =     1.5373E-02    -3.9705E-23</div><div><br></div><div>      iter #   1 total cpu time :  2186.1 secs   <a href="http://av.it" target="_blank">av.it</a>.:  10.8</div><div>      thresh= 1.000E-02 alpha_mix =  0.700 |ddv_scf|^2 =  1.220E-05</div><div> ...</div></div><div><br></div><div><div>     Pert. #  1: Fermi energy shift (Ry) =    -4.1145E+08    -1.3878E-17</div><div><br></div><div>      iter #  99 total cpu time :114535.6 secs   <a href="http://av.it" target="_blank">av.it</a>.:  76.5</div><div>      thresh= 1.000E-02 alpha_mix =  0.700 |ddv_scf|^2 =  1.192E+16</div><div>     kpoint   3 ibnd 211 solve_linter: root not converged  1.297E-02</div><div>     kpoint   6 ibnd 210 solve_linter: root not converged  2.155E-02</div><div>     kpoint  12 ibnd 166 solve_linter: root not converged  3.906E-02</div><div>     kpoint  13 ibnd 167 solve_linter: root not converged  4.197E-02</div><div>     kpoint  15 ibnd 167 solve_linter: root not converged  1.025E-02</div><div>     kpoint  16 ibnd 167 solve_linter: root not converged  3.094E-02</div><div>     kpoint  17 ibnd 169 solve_linter: root not converged  1.276E-02</div><div> </div><div>     Pert. #  1: Fermi energy shift (Ry) =     1.7443E+09    -5.5511E-17</div><div><br></div><div>      iter # 100 total cpu time :116795.0 secs   <a href="http://av.it" target="_blank">av.it</a>.:  75.2</div><div>      thresh= 1.000E-02 alpha_mix =  0.700 |ddv_scf|^2 =  1.063E+17</div><div><br></div><div>     End of self-consistent calculation</div><div><br></div><div>     No convergence has been achieved </div></div><div><br></div><div>Sincerely</div><div><br></div><div>Jibiao Li</div><div><br></div><div>Yangtze Normal University, China</div><div><br></div><div>My input files</div><div>phonons at Gamma</div><div> &inputph</div><div>  tr2_ph=1.0d-12,</div><div>  prefix='T.1ml',</div><div>  amass(1)=15.999,</div><div>  amass(2)=12.001,</div><div>  amass(3)=55.850,</div><div>  outdir='./',</div><div>  fildyn='T.1ml.dynG',</div><div> /</div><div>0.0 0.0 0.0</div></div><div><br></div><div><br></div><div> &CONTROL</div><div>                 calculation = 'scf' ,</div><div>                restart_mode = 'from_scratch' ,</div><div>                      outdir = './' ,</div><div>                  pseudo_dir = '/home/bmllzr/codes/pseudo/' ,</div><div>                      prefix = 'T.1ml' ,</div><div>                     tstress = .true. ,</div><div>                     tprnfor = .true. ,</div><div> /</div><div> &SYSTEM</div><div>                       ibrav = 6,</div><div>                   celldm(1) = 10.847032278,</div><div>                   celldm(3) = 3.5,</div><div>                         nat = 28,</div><div>                        ntyp = 3,</div><div>                     ecutwfc = 35,</div><div>                     ecutrho = 280 ,</div><div>                 occupations = 'smearing' ,</div><div>                     degauss = 0.02D0 ,</div><div>                    smearing = 'methfessel-paxton' ,</div><div>                       nspin = 2 ,</div><div>   starting_magnetization(1) = -0.1,</div><div>   starting_magnetization(2) = -0.1,</div><div>   starting_magnetization(3) = 2.9,</div><div>                    vdw_corr = 'DFT-D',</div><div> /</div><div> &ELECTRONS</div><div>                 mixing_beta = 0.2D0 ,</div><div>             diagonalization = 'david' ,</div><div> /</div><div>&IONS</div><div>/</div><div>ATOMIC_SPECIES</div><div>    O   15.99900  O.pbe-van_ak.UPF </div><div>    C   12.001    C.pbe-van_ak.UPF </div><div>   Fe   55.85000  Fe.pbe-sp-van_ak.UPF </div><div>ATOMIC_POSITIONS angstrom </div><div>C        2.870000022   0.000000000  10.963681264</div><div>O        2.870000002   0.000000000   9.820542707</div><div>C        0.000000000   2.870000002  10.963681264</div><div>O        0.000000000   2.870000002   9.820542707</div><div>O        2.870000022   2.869999989   8.484765130</div><div>C        2.870000002   2.869999976   7.310891865</div><div>O        0.000000000   0.000000000   8.485269386</div><div>C        0.000000000   0.000000000   7.311308054</div><div>Fe      -0.000000001   0.000000000   5.545607710</div><div>Fe       2.870000008   0.000000007   5.773720938</div><div>Fe       0.000000006   2.870000008   5.773720988</div><div>Fe       2.870000000   2.870000009   5.545277888</div><div>Fe       1.435015580   1.435015326   4.303533738</div><div>Fe       4.304984355   1.435015405   4.303533575</div><div>Fe       1.435015640   4.304984586   4.303533577</div><div>Fe       4.304984433   4.304984688   4.303533741</div><div>Fe       0.000000000   0.000000000   2.848495086</div><div>Fe       2.870000000   0.000000000   2.879568308</div><div>Fe       0.000000000   2.870000000   2.879568308</div><div>Fe       2.870000000   2.870000000   2.848497766</div><div>Fe       4.304935556   4.304935556   1.421163132</div><div>Fe       1.435064444   1.435064444   1.421163132</div><div>Fe       4.304935556   1.435064444   1.421163132</div><div>Fe       1.435064444   4.304935556   1.421163132</div><div>Fe       0.000000000   0.000000000   0.066873495</div><div>Fe       2.870000000   0.000000000   0.071056483</div><div>Fe       0.000000000   2.870000000   0.071056483</div><div>Fe       2.870000000   2.870000000   0.066732928</div><div>K_POINTS automatic </div><div>  6 6 1   0 0 0 </div></div><div><br></div><br>______________________________<wbr>_________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" rel="noreferrer" target="_blank">http://pwscf.org/mailman/<wbr>listinfo/pw_forum</a><br></blockquote></div><br></div></div></div></div></div></div>