<div dir="ltr">Dear all,<div><br></div><div>I am curious if this is a common problem: when calculating the phonons of an insulator and a metal I face no problems in slab geometry. But when both are put together there will always be a q-point where convergence cannot be reached, for example:</div><div><br></div><div><div><i>     Computing dynamical matrix for </i></div><div><i>                    q = (   0.0000000   0.1666667   0.0000000 )</i></div><div><i>       2 Sym.Ops. (no q -> -q+G )</i></div><div><i>     G cutoff =  609.4453  (  20928 G-vectors)     FFT grid: ( 50, 50,400)<br></i></div><div><i>     number of k points=   144  gaussian smearing, width (Ry)=  0.0100<br></i></div><div><i><br></i></div><div><i><br></i></div><div><div><i>      iter #  15 total cpu time : 24287.3 secs   <a href="http://av.it">av.it</a>.:  13.0</i></div><div><i>      thresh= 4.091E-07 alpha_mix =  0.700 |ddv_scf|^2 =  8.763E-11</i></div></div><div><i>...</i></div><div><i>      iter #  59 total cpu time : 27386.7 secs   <a href="http://av.it">av.it</a>.:  11.2<br></i></div><div><div><i>      thresh= 5.160E-05 alpha_mix =  0.700 |ddv_scf|^2 =  1.219E-06</i></div></div><div><i>..<br></i></div><div><div><i>      iter # 100 total cpu time : 30214.5 secs   <a href="http://av.it">av.it</a>.:  10.8</i></div><div><i>      thresh= 1.000E-02 alpha_mix =  0.700 |ddv_scf|^2 =  4.154E-02</i></div><div><i><br></i></div><div><i>     End of self-consistent calculation</i></div><div><i><br></i></div><div><i>     No convergence has been achieved</i></div></div><div><br></div><div>Due to the oscillating behavior of ddv_scf I doubt that increasing the iterations will in the end work. </div><div><br></div><div>I tried increasing k points and cutoff well beyond SCF convergence (SCF shows convergence at 12x12x1 and ecutwfc~80 Ry but I went up to 36x36x1 and 120 Ry.)</div><div><br></div><div><div>&inputph</div><div>  prefix   = 'mgo',</div><div>  trans    = .true.</div><div>  epsil    = .false.,</div><div>  fildyn   = 'mgo.dyn',</div><div>  ldisp    = .true.</div><div>  fildvscf = 'dvscf'</div><div>  amass(1) = 107.86820</div><div>  amass(2) = 24.30500</div><div>  amass(3) = 15.99990</div><div>  outdir='./',</div><div>  nq1=6,  nq2=6,  nq3=1,</div><div>  tr2_ph   =  1.0d-12,<br></div><div> /</div></div><div><br></div><div>The system is MgO stacked on top of Ag. The pseudos are NC from the ONCV library and work very well for the "pure" Ag and MgO separately.</div><div><br></div><div>Any hint is very welcome!</div><div><br></div><div>Best,</div><div>Chris </div>-- <br><div 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>