<div dir="ltr"><div>Dear pwscf users and developpers, <br><br></div><div>Best wishes for the New Year! <br></div><div><br></div>I have trouble with phonon calculations on magnetic insulators: the calculations works well at gamma (k=0) point, and I am able to compute dielectric tensor and effective charges, but as soon as I try out of the zone center, the calculation never works, with errors such as (in the ph output): <br><br>------------------------------<wbr>------------------------------<wbr>---------------------------<br>    Representation #  1 mode #   1<br><br>     Self-consistent Calculation<br>     kpoint   1 ibnd  75 solve_linter: root not converged  1.071E-01<br>     kpoint   2 ibnd  75 solve_linter: root not converged  1.924E-01<br>     kpoint   3 ibnd  75 solve_linter: root not converged  1.745E-02<br>     kpoint   4 ibnd  75 solve_linter: root not converged  1.722E-01<br>------------------------------<wbr>------------------------------<wbr>------------------------------<wbr>-<br clear="all"><div><div><br></div><div>On the other hand, when I do the very same calculation, but just changing the insulator state to a metallic state (in the preceeding scf step, that is), then the calculation goes perfectly well. Is there a particular reason why that shouldn't work? The electronic bands show that this system is indeed an insulator with a significant gap.  <br><br></div><div>Thank you very much for your help,<br><br></div><div>Best Regards, <br><br></div><div>Merlin Méheut<br></div><div><br></div><div>1/ scf input for insulator state:<br>&system<br>    ibrav =0,<br>    celldm(1) =  12.459185,<br>    nat =20, ntyp = 4,<br>    ecutwfc =80.0, ecutrho = 320.0,<br>    nspin=2, starting_magnetization(1)=0.5,<br>    tot_magnetization=4,<br>/&end<br><br></div><div>2/ for metallic state<br></div><div>&system<br>    ibrav =0,<br>    celldm(1) =  12.459185,<br>    nat =20, ntyp = 4,<br>     ecutwfc =80.0, ecutrho = 320.0, <br>    nspin=2, starting_magnetization(1)=0.5,<br>    occupations='smearing', smearing='mp', degauss=0.04,<br>/&end<br><br></div><div>3/ phonon input: <br>&inputph<br>   amass(1)=57.9353,<br>   amass(2)=39.9626,<br>   amass(3)=27.9769,<br>   amass(4)=15.9949,<br>   alpha_mix(1) = 0.5,<br>   ! ldisp=.true., nq1=2, nq2=3, nq3=2,<br>   tr2_ph =  1.0D-17,<br>   prefix='NiCPX',<br>   fildyn='mat.$PREFIX',<br>   lraman=.false.,<br>   epsil =.false.,<br>   trans =.true.,<br>   zue = .false.,<br>   zeu= .false.,<br>   outdir='/tmpdir/$LOGNAME/'<br>/&end<br>  -0.0006214  -0.0004751  -0.5041257<br><br></div><div><br>-- <br><div class="gmail-m_2827493750489094115gmail-m_-2626443323123650728gmail_signature"><div dir="ltr">Merlin Méheut<br>adresse labo:<br>GET - OMP   - Université Paul Sabatier <br>14 avenue Edouard Belin<br>31400 Toulouse<br>FRANCE <br>tel: <a href="tel:+33%205%2061%2033%2026%2017" value="+33561332617" target="_blank">(+33) 5 61 33 26 17</a><br><br></div></div>
</div></div></div>