<div dir="ltr"><div><span style="font-family:georgia,serif">   Dear QE users,</span></div><div><span style="font-family:georgia,serif">I have massive problems with converging the phonon scf for a layered TMD (bulk calculation).</span></div><div><span style="font-family:georgia,serif">Here are the last lines of the PHONON output</span></div><div><br> **************************************************************************</div><div>      iter #  95 total cpu time :  3975.1 secs   <a href="http://av.it" target="_blank">av.it</a>.:   8.3<br>      thresh= 4.694E-07 alpha_mix =  0.300 |ddv_scf|^2 =  6.601E-12<br><br>      iter #  96 total cpu time :  3987.5 secs   <a href="http://av.it" target="_blank">av.it</a>.:   8.3<br>      thresh= 2.569E-07 alpha_mix =  0.300 |ddv_scf|^2 =  7.102E-12<br><br>      iter #  97 total cpu time :  3999.9 secs   <a href="http://av.it" target="_blank">av.it</a>.:   9.1<br>      thresh= 2.665E-07 alpha_mix =  0.300 |ddv_scf|^2 =  2.096E-11<br><br>      iter #  98 total cpu time :  4011.2 secs   <a href="http://av.it" target="_blank">av.it</a>.:   8.4<br>      thresh= 4.579E-07 alpha_mix =  0.300 |ddv_scf|^2 =  1.120E-12<br><br>      iter #  99 total cpu time :  4023.8 secs   <a href="http://av.it" target="_blank">av.it</a>.:  10.8<br>      thresh= 1.058E-07 alpha_mix =  0.300 |ddv_scf|^2 =  8.827E-12<br><br>      iter # 100 total cpu time :  4035.6 secs   <a href="http://av.it" target="_blank">av.it</a>.:   8.9<br>      thresh= 2.971E-07 alpha_mix =  0.300 |ddv_scf|^2 =  1.223E-11<br><br>     End of self-consistent calculation<br><br>     No convergence has been achieved <br> **************************************************************************</div><div><br></div><div><span style="font-family:georgia,serif">The scf hits</span> |ddv_scf|^2 =  7.292E-11 <span style="font-family:georgia,serif">at iteration # 30.</span></div><div><span style="font-family:georgia,serif">Single layer calculations never had this problem.</span></div><div><span style="font-family:georgia,serif">The structure has been relaxed up to    "forc_conv_thr                 = 2.0D-5".</span></div><div><span style="font-family:georgia,serif">This problem never comes out at Gamma, and trying coarser meshes usually gives no problems with the convergence.<br></span></div><div><span style="font-family:georgia,serif">I attach the input and output of the PHONON for reference.</span></div><div><span style="font-family:georgia,serif">Any help will be highly appreciated.</span></div><div><span style="font-family:georgia,serif"><br></span></div><div><span style="font-family:georgia,serif">Thank you in advance,</span></div><div><span style="font-family:georgia,serif">Fabrizio</span></div><br clear="all"><br>-- <br><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><b>Fabrizio Cossu<br>postdoctoral fellow at APCTP (Asia Pacific Center for Theoretical Physics)</b>,</div><div dir="ltr"><pre>                       |<br>  .----. .----. .----. |----  .----.   ===========================<br>  ,----| |    | |      |      |    |   <a href="http://www.apctp.org/?JrgId=16" style="font-family:"comic sans ms",sans-serif;white-space:normal;color:rgb(17,85,204)" target="_blank">http://www.apctp.org/?JrgId=16</a><br>  `----^ |----' `----' `----' |----'   ===========================<br>         |                    |</pre></div><div dir="ltr"><font face="comic sans ms, sans-serif"><span><span></span></span></font></div></div></div></div></div></div></div></div></div></div></div></div></div></div>