<div dir="ltr"><div dir="ltr"><div dir="ltr">Dear all,<div><br></div><div>thank you for your help so far. Unfortunately I am still not able to reach self consistency during the ph.x calculation. I have tested all smearing methods (and tetrahedra_opt), several smearing widths between 0.001 and 0.02, in-plane k-point grids up to 21x21x1 (121 k points) and relaxed the structure down to forces of 1e-6 a.u.. <br clear="all"><div><br></div><div>whilst I can reach arbitrary scf accuracy</div><div><br></div><div><div>!    total energy              =   -3634.01568752 Ry</div><div>     Harris-Foulkes estimate   =   -3634.01568752 Ry</div><div>     estimated scf accuracy    <          5.9E-17 Ry</div></div><div><br></div><div>there seems to be no way to avoid the "oscillating" behavior in ph.x </div><div><div>     Pert. #  1: Fermi energy shift (Ry) =    -9.6471E-02     0.0000E+00</div><div><br></div><div>      iter # 100 total cpu time : 23216.2 secs   <a href="http://av.it">av.it</a>.:   7.2</div><div>      thresh= 1.733E-04 alpha_mix =  0.100 |ddv_scf|^2 =  1.665E-09</div><div>     End of self-consistent calculation<br></div><div>     No convergence has been achieved<br></div></div><div><br></div><div>This is very odd and very frustrating - I am curious if anyone else has encountered this for thin (7-11 ML) metal layers and has found a strategy to overcome this?</div><div><br></div><div>Thanks in advance for your help, it is very much appreciated!</div><div><br></div><div>Best,</div><div>Chris </div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr">Postdoctoral Researcher<br>Center for Quantum Nanoscience, Institute for Basic Science<br>Ewha Womans University, Seoul, South Korea<blockquote type="cite" style="font-size:12.8px"><div dir="ltr"><div><div dir="ltr"></div></div></div></blockquote></div></div></div></div></div></div>