<div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div>Dear all,</div><div><br></div><div>I have observed that the convergence of phonon calculations sometimes strongly depends on the choice of q-point without being able to understand which points should better be avoided.</div><div><br></div><div>Take for example graphene. I sample it with 33x33x1 k points for the scf step. This choice of grid contains, among others</div><div><br></div><div>     .....<br>        k(  106) = (   0.3030303   0.5248639   0.0000000), wk =   0.0110193<br>        k(  107) = (   0.3030303   0.5598548   0.0000000), wk =   0.0220386<br>   **   k(  108) = (   0.3333333   0.5773503   0.0000000), wk =   0.0036731 **</div><div><br></div><div>The latest k point (in cartesian units) is the K-point where the diract point lies. Subsequently I find that a phonon calculation that contains q=k(108) such as a 15x15x1 grid</div><div><br></div>     Dynamical matrices for (15,15, 1)  uniform grid of q-points<br>     (  27 q-points):<br>       N         xq(1)         xq(2)         xq(3)<br>       1   0.000000000   0.000000000   0.000000000<br>       2   0.000000000   0.076980036   0.000000000<br>....<br>      25   0.266666667   0.461880215   0.000000000<br>      26   0.266666667   0.538860251   0.000000000<br>  *** 27   0.333333333   0.577350269   0.000000000 ***<div><br></div><div>will fail as follows:</div><div></div><div><br>     Representation #  3 modes #   3  4<br>     Self-consistent Calculation<br>      iter #   1 total cpu time : 64587.1 secs   <a href="http://av.it">av.it</a>.:   7.1<br>      thresh= 1.000E-02 alpha_mix =  0.700 |ddv_scf|^2 =  2.917E+03<br>      ......<br>      iter #  14 total cpu time : 65207.4 secs   <a href="http://av.it">av.it</a>.:  45.2<br>      thresh= 1.000E-02 alpha_mix =  0.700 |ddv_scf|^2 =  5.244E+40<br></div><div><br></div><div>(playing with alpha_mix(niter) does not seem to help much here).</div><div><br></div><div>Now the question I have been asking myself is if there is a "good" and a "bad" choice of SCF k points and phonon q-points? If someone has any experience related to the proper choice of points I would be very glad to hear!</div><div><br></div><div>Thanks in advance and with best regards,</div><div>Chris <br></div><div><br>-- <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></div>