<div dir="ltr"><div dir="ltr"><div>Dear Iurii, <br></div><div><br></div><div>I learned a lot with your guidance, thank you so much! Unfortunately, I am afraid that none of your propositions (use SSSP pseudopotentials, change density cutoff) did solve the matter (1). To further explore this question, I also did DFPT+U phonon calculation on FeO (2), with an input almost identical to your example02 on NiO, and there was a divergence for the second representation, very similar to what I have for Fe2SiO4. At last (3), I tried to compute the +U value based on your example02 as well. For fayalite, I succeeded to compute U for 1x1x1 q-point grid, but 2x1x2 and 2x2x2 grids failed. For FeO, the situation is more complex, depending on the k-point grid and on the q-point grid, sometimes it works, sometimes it does not converge, sometimes it crashes. But generally for large grids it is not possible to reach the end of the calculation. Below are more detailed descriptions of these problems:<br></div><div><br></div><div><div>The procedures were tested for pseudopotentials <br></div><div>O.pz-n-rrkjus_psl.0.1.UPF + Fe.pz-spn-rrkjus_psl.0.2.1.UPF + Si.pz-n-rrkjus_psl.0.1.UPF <br>Fe.pbesol-spn-kjpaw_psl.0.2.1.UPF + Si.pbesol-n-rrkjus_psl.1.0.0.UPF + O.pbesol-n-kjpaw_psl.0.1.UPF<br></div><div>and the formerly described pseudopotentials<br></div><div><br></div></div><div>(1) For DFPT+U phonon calculations on Fayalite, I followed your suggestion to adopt the approach described in example 02, with a first scf calculation in the metallic state: <br></div><div>&system<br> (...)<br> nspin=2,<br> occupations='smearing', <br> degauss=0.001, <br> smearing= 'gauss',<br> starting_magnetization(1)=0.5, <br> starting_magnetization(2)=-0.5,<br> lda_plus_u = .true., <br> lda_plus_u_kind = 0,<br> U_projection_type = 'ortho-atomic',<br> Hubbard_U(1)=2.5, <br> Hubbard_U(2)=2.5, <br>/&end <br></div><div><br></div><div>&electrons<br> electron_maxstep = 200,<br> conv_thr = 1.d-10,<br> mixing_beta = 0.4,<br>/&end</div><div>(...)<br></div><div>K_POINTS {automatic}<br>2 2 2 1 1 1 <br></div><div><br></div><div>and a second calculation as insulating material: <br></div><div>&system<br> ibrav =8, <br> celldm(1)=9.107535, <br> celldm(2)=2.17425,<br> celldm(3)=1.263056,<br> nat =28, <br> ntyp = 4, <br> ecutwfc =80.0, <br> ecutrho = 320.0,<br> nbnd=${nstates},<br> nspin=2,<br> occupations='fixed',<br> tot_magnetization= 0, <br> lda_plus_u = .true., <br> lda_plus_u_kind = 0,<br> U_projection_type = 'atomic',<br> Hubbard_U(1)=2.5, <br> Hubbard_U(2)=2.5, <br></div><div>/&end <br> &electrons<br> electron_maxstep = 200,<br> conv_thr = 1.d-15,<br> mixing_mode = 'plain',<br> mixing_beta = 0.4,<br> startingpot = 'file',<br> startingwfc = 'file',<br>/&end<br></div><div><br></div><div>followed by the phonon calculation: <br></div><div>Liz m1<br> &inputph<br> alpha_mix(1) = 0.3,<br> tr2_ph = 1.0D-17,<br> prefix='FAYA2',<br> fildyn='mat.$PREFIX',<br> lraman=.false.,<br> epsil =.false.,<br> trans =.true.,<br> zue = .false.,<br> zeu= .false.,<br> start_irr=23,<br> last_irr=23,<br> niter_ph=20, <br> outdir ='./',<br>/&end<br>0.0 0.0 0.0<br></div><div><br></div><div>The procedure was tested for pseudopotentials <br></div><div>O.pz-n-rrkjus_psl.0.1.UPF + Fe.pz-spn-rrkjus_psl.0.2.1.UPF + Si.pz-n-rrkjus_psl.0.1.UPF <br>Fe.pbesol-spn-kjpaw_psl.0.2.1.UPF + Si.pbesol-n-rrkjus_psl.1.0.0.UPF + O.pbesol-n-kjpaw_psl.0.1.UPF<br></div><div>and the formerly described<br></div><div><br></div><div>I also changed cutoffs as you suggested. There is no difference that I can see in the result. <br></div><div><br></div><div>(2) For DFPT+U phonon calculations on FeO, I adopted the same approach, with two scf calculations very similar to the one shown in HP example02, except for U_projection_type = 'atomic', instead of 'ortho-atomic' for the second scf . The following phonon input looks like this: <br></div><div><br>Liz m1<br> &inputph<br> alpha_mix(1) = 0.3,<br> tr2_ph = 1.0D-15,<br> (...)<br> lraman=.false.,<br> epsil =.false.,<br> trans =.true., <br> zue = .false.,<br> zeu= .false.,<br> start_irr=2,<br> last_irr=2,<br> niter_ph=20,<br>/&end<br>0.0 0.0 0.0<br></div><div>Here also, we obtain a very clear divergence, whatever the Fe and O pseudopotentials considered above. <br></div><div><br></div><div>(3) Then I tried to compute the value of the +U parameter, based on your HP example02. I did it for fayalite and for FeO (in which case the input is essentially the same than example02, except for the pseudopotentials). Following your article (Timrov et al 2018, PRB 98, 085127), I tried different values for the k-point and q-point grids . Some work, some do not converge, some diverge. For fayalite, it was only possible to converge for 1x1x1 q-point grid. For FeO as well, large grids tend to be just impossible to complete, as some particular q-points do not converge. Once again, whatever the pseudo or the cutoff. However, these are not exactly the same configurations that fail depending on the pseudopotentials. <br></div><div><br></div><div>If you have any hints that may permit to progress on those issues that would be awesome! <br></div><div><br></div><div>Best regards,<br></div><div><br></div><div>Merlin <br></div><div><br></div><div> -- <br><div dir="ltr"><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: (+33) 5 61 33 26 17<br></div></div></div></div><div><a href="http://www3.obs-mip.fr/get/profils/Meheut_Merlin">http://www3.obs-mip.fr/get/profils/Meheut_Merlin</a></div><div><br></div><div><br></div><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Le jeu. 20 févr. 2020 à 10:05, Timrov Iurii <<a href="mailto:iurii.timrov@epfl.ch" target="_blank">iurii.timrov@epfl.ch</a>> a écrit :<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div>
<div id="gmail-m_5681557970449129721gmail-m_8275223278710308024gmail-m_4073407159693916210divtagdefaultwrapper" style="font-size:12pt;color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif" dir="ltr">
<p>Dear Merlin,</p>
<p><br>
</p>
<p>By the way, since you are using DFPT+U to compute phonons including the Hubbard U correction, I would like to let you (and others) know that the paper describing DFPT+U has been just published in Physical Review B:</p>
<p></p>
<h5><span style="font-family:Arial,Helvetica,sans-serif;font-size:12pt">A. Floris, I. Timrov, B. Himmetoglu, N. Marzari, S. de Gironcoli, and M. Cococcioni, "</span><span style="font-family:Arial,Helvetica,sans-serif;font-size:12pt">Hubbard-corrected
density functional perturbation theory with ultrasoft pseudopotentials</span><span style="font-family:Arial,Helvetica,sans-serif;font-size:12pt">", </span><span style="font-family:Arial,Helvetica,sans-serif;font-size:12pt">Phys. Rev. B
</span><span style="font-family:Arial,Helvetica,sans-serif;font-size:12pt">101</span><span style="font-family:Arial,Helvetica,sans-serif;font-size:12pt">, 064305 (2020).</span></h5>
Best regards,
<p></p>
<p>Iurii<br>
</p>
<p><br>
</p>
<div id="gmail-m_5681557970449129721gmail-m_8275223278710308024gmail-m_4073407159693916210Signature">
<div id="gmail-m_5681557970449129721gmail-m_8275223278710308024gmail-m_4073407159693916210divtagdefaultwrapper" dir="ltr" style="font-size:12pt;color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols">
<div name="divtagdefaultwrapper">
<div name="divtagdefaultwrapper">
<font size="3" face="'Times New Roman', Times, serif" color="808080">--<br>
Dr. Iurii Timrov<br>
Postdoctoral Researcher<br>
<font color="808080"><font face="'Times New Roman', Times, serif"></font></font></font></div>
<font color="808080"></font>
<div name="divtagdefaultwrapper">
<font size="3" face="'Times New Roman', Times, serif" color="808080">STI - IMX <font color="808080">
<font face="'Times New Roman', Times, serif">- THEOS</font></font></font><font size="3" face="'Times New Roman', Times, serif" color="808080"> and NCCR - MARVEL<br>
</font></div>
<div name="divtagdefaultwrapper">
<font size="3" face="'Times New Roman', Times, serif" color="808080"><font size="3" face="'Times New Roman', Times, serif" color="808080">Swiss Federal Institute of Technology Lausanne (EPFL<font color="808080"><font face="'Times New Roman', Times, serif">)</font></font></font><br>
</font></div>
<font color="808080"></font>
<div name="divtagdefaultwrapper">
<font size="3" face="'Times New Roman', Times, serif" color="808080">CH-1015 Lausanne, Switzerland<br>
+41 21 69 34 881</font></div>
<div name="divtagdefaultwrapper">
<a href="http://people.epfl.ch/265334" id="gmail-m_5681557970449129721gmail-m_8275223278710308024gmail-m_4073407159693916210LPNoLP" target="_blank">http://people.epfl.ch/265334</a><br>
</div>
</div>
</div>
</div>
</div>
<hr style="display:inline-block;width:98%">
<div id="gmail-m_5681557970449129721gmail-m_8275223278710308024gmail-m_4073407159693916210divRplyFwdMsg" dir="ltr"><font style="font-size:11pt" face="Calibri, sans-serif" color="#000000"><b>From:</b> users <<a href="mailto:users-bounces@lists.quantum-espresso.org" target="_blank">users-bounces@lists.quantum-espresso.org</a>> on behalf of Timrov Iurii <<a href="mailto:iurii.timrov@epfl.ch" target="_blank">iurii.timrov@epfl.ch</a>><br>
<b>Sent:</b> Wednesday, February 19, 2020 10:05:01 AM<br>
<b>To:</b> Quantum ESPRESSO users Forum<br>
<b>Subject:</b> Re: [QE-users] non-convergency of DFPT+U calculation</font>
<div> </div>
</div>
<div>
<div id="gmail-m_5681557970449129721gmail-m_8275223278710308024gmail-m_4073407159693916210divtagdefaultwrapper" style="font-size:12pt;color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif" dir="ltr">
<p>Dear Merlin,</p>
<p><br>
</p>
<p>Total magnetization must be zero for an antiferromagnetic insulator. Hence, some deviation from zero that you obtain is probably a numerical noise (I would check the convergence of the results: reduce
<span>conv_thr</span>, increase the size of the k mesh, etc., and relax well the structure).
<br>
</p>
<p><br>
</p>
<p><span>> If I wanted to be closer to the smearing calculation, I would have to consider nocolin=.true., is that right?
</span><br>
</p>
<p><br>
</p>
<p>DFPT+U with noncolin=.true. is not implemented.</p>
<p><br>
</p>
<p>Greetings,</p>
<p>Iurii<br>
</p>
<p><br>
</p>
<div id="gmail-m_5681557970449129721gmail-m_8275223278710308024gmail-m_4073407159693916210Signature">
<div id="gmail-m_5681557970449129721gmail-m_8275223278710308024gmail-m_4073407159693916210divtagdefaultwrapper" dir="ltr" style="font-size:12pt;color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols">
<div name="divtagdefaultwrapper">
<div name="divtagdefaultwrapper">
<font size="3" face="'Times New Roman', Times, serif" color="808080">--<br>
Dr. Iurii Timrov<br>
Postdoctoral Researcher<br>
<font color="808080"><font face="'Times New Roman', Times, serif"></font></font></font></div>
<font color="808080"></font>
<div name="divtagdefaultwrapper">
<font size="3" face="'Times New Roman', Times, serif" color="808080">STI - IMX <font color="808080">
<font face="'Times New Roman', Times, serif">- THEOS</font></font></font><font size="3" face="'Times New Roman', Times, serif" color="808080"> and NCCR - MARVEL<br>
</font></div>
<div name="divtagdefaultwrapper">
<font size="3" face="'Times New Roman', Times, serif" color="808080"><font size="3" face="'Times New Roman', Times, serif" color="808080">Swiss Federal Institute of Technology Lausanne (EPFL<font color="808080"><font face="'Times New Roman', Times, serif">)</font></font></font><br>
</font></div>
<font color="808080"></font>
<div name="divtagdefaultwrapper">
<font size="3" face="'Times New Roman', Times, serif" color="808080">CH-1015 Lausanne, Switzerland<br>
+41 21 69 34 881</font></div>
<div name="divtagdefaultwrapper">
<a href="http://people.epfl.ch/265334" id="gmail-m_5681557970449129721gmail-m_8275223278710308024gmail-m_4073407159693916210LPNoLP" target="_blank">http://people.epfl.ch/265334</a><br>
</div>
</div>
</div>
</div>
</div>
<hr style="display:inline-block;width:98%">
<div id="gmail-m_5681557970449129721gmail-m_8275223278710308024gmail-m_4073407159693916210divRplyFwdMsg" dir="ltr"><font style="font-size:11pt" face="Calibri, sans-serif" color="#000000"><b>From:</b> users <<a href="mailto:users-bounces@lists.quantum-espresso.org" target="_blank">users-bounces@lists.quantum-espresso.org</a>> on behalf of Merlin Meheut <<a href="mailto:merlin.meheut@gmail.com" target="_blank">merlin.meheut@gmail.com</a>><br>
<b>Sent:</b> Wednesday, February 19, 2020 12:38:08 AM<br>
<b>To:</b> Quantum ESPRESSO users Forum<br>
<b>Subject:</b> Re: [QE-users] non-convergency of DFPT+U calculation</font>
<div> </div>
</div>
<div>
<div dir="ltr">
<div dir="ltr">
<div>Dear Iurii, <br>
</div>
<div><br>
</div>
<div>Thank you again for your help! So far, the Example02 seems quite clear to me, although it will take me some time to apply it to my system and check that everything is working. One quick question, though: in my case, the "total magnetization" obtained from
the first calculation (with smearing) is non-zero and non-integer: "total magnetization = -0.63 Bohr mag/cell" (which makes sense since I have two different sites). As I understand it, for the second calculation, the variable tot_magnetization should
be integer, so I can only take it equal to 0 or 1 here (I tried with 0, I don't know which is better). If I wanted to be closer to the smearing calculation, I would have to consider nocolin=.true., is that right?
<br>
</div>
<div><br>
</div>
<div>Greetings, <br>
</div>
<div><br>
</div>
<div>Merlin <br>
</div>
</div>
<br>
<div class="gmail_quote">
<div dir="ltr" class="gmail_attr">Le mar. 18 févr. 2020 à 17:32, Timrov Iurii <<a href="mailto:iurii.timrov@epfl.ch" target="_blank">iurii.timrov@epfl.ch</a>> a écrit :<br>
</div>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div>
<div id="gmail-m_5681557970449129721gmail-m_8275223278710308024gmail-m_4073407159693916210gmail-m_8546846971254182096divtagdefaultwrapper" style="font-size:12pt;color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif" dir="ltr">
<p>Dear Merlin,</p>
<p><br>
</p>
<p><span>> Does that mean that magnetic insulators cannot be treated by Quantum Espresso?
</span><br>
</p>
<p><br>
</p>
<p>It is possible to model magnetic insulators in Quantum ESPRESSO. I am using the procedure described in
<span>q-e-qe-6.5/HP/examples/example02</span>/<span>README</span>. If the description in README file is not clear enough, please let me know.</p>
<p><br>
</p>
<p>Greetings,</p>
<p>Iurii<br>
</p>
<p><br>
</p>
<div id="gmail-m_5681557970449129721gmail-m_8275223278710308024gmail-m_4073407159693916210gmail-m_8546846971254182096Signature">
<div id="gmail-m_5681557970449129721gmail-m_8275223278710308024gmail-m_4073407159693916210gmail-m_8546846971254182096divtagdefaultwrapper" dir="ltr" style="font-size:12pt;color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif,"EmojiFont","Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols">
<div name="divtagdefaultwrapper">
<div name="divtagdefaultwrapper"><font size="3" face="'Times New Roman', Times, serif" color="808080">--<br>
Dr. Iurii Timrov<br>
Postdoctoral Researcher<br>
<font color="808080"><font face="'Times New Roman', Times, serif"></font></font></font></div>
<font color="808080"></font>
<div name="divtagdefaultwrapper"><font size="3" face="'Times New Roman', Times, serif" color="808080">STI - IMX
<font color="808080"><font face="'Times New Roman', Times, serif">- THEOS</font></font></font><font size="3" face="'Times New Roman', Times, serif" color="808080"> and NCCR - MARVEL<br>
</font></div>
<div name="divtagdefaultwrapper"><font size="3" face="'Times New Roman', Times, serif" color="808080"><font size="3" face="'Times New Roman', Times, serif" color="808080">Swiss Federal Institute of Technology Lausanne (EPFL<font color="808080"><font face="'Times New Roman', Times, serif">)</font></font></font><br>
</font></div>
<font color="808080"></font>
<div name="divtagdefaultwrapper"><font size="3" face="'Times New Roman', Times, serif" color="808080">CH-1015 Lausanne, Switzerland<br>
+41 21 69 34 881</font></div>
<div name="divtagdefaultwrapper"><a href="http://people.epfl.ch/265334" id="gmail-m_5681557970449129721gmail-m_8275223278710308024gmail-m_4073407159693916210gmail-m_8546846971254182096LPNoLP" target="_blank">http://people.epfl.ch/265334</a><br>
</div>
</div>
</div>
</div>
</div>
<hr style="display:inline-block;width:98%">
<div id="gmail-m_5681557970449129721gmail-m_8275223278710308024gmail-m_4073407159693916210gmail-m_8546846971254182096divRplyFwdMsg" dir="ltr"><font style="font-size:11pt" face="Calibri, sans-serif" color="#000000"><b>From:</b> users <<a href="mailto:users-bounces@lists.quantum-espresso.org" target="_blank">users-bounces@lists.quantum-espresso.org</a>>
on behalf of Merlin Meheut <<a href="mailto:merlin.meheut@gmail.com" target="_blank">merlin.meheut@gmail.com</a>><br>
<b>Sent:</b> Tuesday, February 18, 2020 5:21:46 PM<br>
<b>To:</b> Quantum ESPRESSO users Forum<br>
<b>Subject:</b> Re: [QE-users] non-convergency of DFPT+U calculation</font>
<div> </div>
</div>
<div>
<div dir="ltr">
<div dir="ltr">
<div>Hi Iurii, <br>
</div>
<div><br>
</div>
<div>First thank you very much for this remark, I did not know that one could not use tot_magnetization and starting_magnetization together. Indeed, this is reported in the Documentation INPUT_PW.def.
<br>
</div>
<div>At this point, however, I don't see how to make an input for an antiferromagnetic insulating material (which is the goal here). Because with this input, the state is non magnetic:
<br>
</div>
<div><br>
</div>
<div>
<div>
<div>&system<br>
ibrav =8, celldm(1)=9.107535, celldm(2)=2.17425,<br>
celldm(3)=1.263056,<br>
nat =28, ntyp = 4, ecutwfc =80.0, ecutrho = 320.0,<br>
nspin=2, tot_magnetization=0,<br>
lda_plus_u = .true., U_projection_type = 'atomic',<br>
Hubbard_U(1)=2.4, Hubbard_U(2)=2.6, ! valeurs Yu et al 2013 =2.4/2.6<br>
/&end </div>
</div>
</div>
<div><br>
</div>
<div>And if you use this input: <br>
</div>
<div><br>
</div>
<div>&system<br>
ibrav =8, celldm(1)=9.107535, celldm(2)=2.17425,<br>
celldm(3)=1.263056,<br>
nat =28, ntyp = 4, ecutwfc =${a}.0, ecutrho = ${b}.0,<br>
nspin=2,<br>
starting_magnetization(1)=0.5, starting_magnetization(2)=-0.5,<br>
lda_plus_u = .true., U_projection_type = 'atomic',<br>
Hubbard_U(1)=2.5, Hubbard_U(2)=2.5,<br>
/&end </div>
<div><br>
</div>
<div>,you have an error message:</div>
<div><br>
</div>
<div>%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br>
Error in routine iosys (1):<br>
fixed occupations and lsda need tot_magnetization<br>
<br>
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</div>
<div><br>
</div>
<div>So as you say, using starting_magnetization imposes using smearing. But if I am correct, smearing means that your material is metallic, and therefore you cannot compute dielectric quantities or effective charges with phonon. Actually, in the INPUT_PW.def
file, the word "antiferromagnetic" is only present once, in the "starting_magnetization" definitions. Does that mean that magnetic insulators cannot be treated by Quantum Espresso?
<br>
</div>
<div><br>
</div>
<div>Would you have any hint as to how to extract myself from such a dead end? <br>
</div>
<div><br>
</div>
<div>Thanks again for your help, <br>
</div>
<div><br>
</div>
<div>Greetings, <br>
</div>
<div><br>
</div>
<div>Merlin<br>
</div>
</div>
<br>
<div class="gmail_quote">
<div dir="ltr" class="gmail_attr">Le dim. 16 févr. 2020 à 23:11, Timrov Iurii <<a href="mailto:iurii.timrov@epfl.ch" target="_blank">iurii.timrov@epfl.ch</a>> a écrit :<br>
</div>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div>
<div id="gmail-m_5681557970449129721gmail-m_8275223278710308024gmail-m_4073407159693916210gmail-m_8546846971254182096gmail-m_2454870331033590078divtagdefaultwrapper" style="font-size:12pt;color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif" dir="ltr">
<p>Dear Merlin,</p>
<p><br>
</p>
<p>In your input for the PW calculation you have:</p>
<p><br>
</p>
<p><span style="color:rgb(33,33,33);font-family:wf_segoe-ui_normal,"Segoe UI","Segoe WP",Tahoma,Arial,sans-serif,serif,EmojiFont;font-size:15px">> tot_magnetization=0,</span><br style="color:rgb(33,33,33);font-family:wf_segoe-ui_normal,"Segoe UI","Segoe WP",Tahoma,Arial,sans-serif,serif,EmojiFont;font-size:15px">
<span style="color:rgb(33,33,33);font-family:wf_segoe-ui_normal,"Segoe UI","Segoe WP",Tahoma,Arial,sans-serif,serif,EmojiFont;font-size:15px">> starting_magnetization(1)=0.5, starting_magnetization(2)=-0.5</span><br>
</p>
<p><br>
</p>
<p>tot_magnetization and starting_magnetization must not be used together (please check the documentation). Moreover, starting_magnetization must be used with occupations='smearing', but I do not see smearing in your input. Is the ground state metallic or insulating? </p>
<p><br>
</p>
<p>How many representations do you have in total? And only #23 diverges?</p>
<p><br>
</p>
<p>Not sure this will help, but you may try to use dual = ecutrho/ecutwfc = 12, with ecutwfc=80 Ry. </p>
<p><br>
</p>
<p><span style="color:rgb(33,33,33);font-family:wf_segoe-ui_normal,"Segoe UI","Segoe WP",Tahoma,Arial,sans-serif,serif,EmojiFont;font-size:15px">> Si 27.9769 Si.pz2.UPF</span><br style="color:rgb(33,33,33);font-family:wf_segoe-ui_normal,"Segoe UI","Segoe WP",Tahoma,Arial,sans-serif,serif,EmojiFont;font-size:15px">
<span style="color:rgb(33,33,33);font-family:wf_segoe-ui_normal,"Segoe UI","Segoe WP",Tahoma,Arial,sans-serif,serif,EmojiFont;font-size:15px">> O 15.9949 O.pz.UPF</span><br>
</p>
<p><br>
</p>
<p>I do not know what are these pseudos, and how good/bad are they. You can try to check the SSSP library.</p>
<p><br>
</p>
<p><span style="color:rgb(33,33,33);font-family:wf_segoe-ui_normal,"Segoe UI","Segoe WP",Tahoma,Arial,sans-serif,serif,EmojiFont;font-size:15px">> Fe1 55.9349 Fe.pz-n-rrkjus_psl.1.0.0.UPF </span><br style="color:rgb(33,33,33);font-family:wf_segoe-ui_normal,"Segoe UI","Segoe WP",Tahoma,Arial,sans-serif,serif,EmojiFont;font-size:15px">
<span style="color:rgb(33,33,33);font-family:wf_segoe-ui_normal,"Segoe UI","Segoe WP",Tahoma,Arial,sans-serif,serif,EmojiFont;font-size:15px">> Fe2 55.9349 Fe.pz-n-rrkjus_psl.1.0.0.UPF</span><br>
</p>
<p><br>
</p>
<p>Also for these I would check the <span style="font-size:12pt">SSSP library.</span></p>
<p><span style="font-size:12pt"><br>
</span></p>
<p><span style="font-size:12pt">Greetings,</span></p>
<p><span style="font-size:12pt">Iurii</span></p>
<p><br>
</p>
<div id="gmail-m_5681557970449129721gmail-m_8275223278710308024gmail-m_4073407159693916210gmail-m_8546846971254182096gmail-m_2454870331033590078Signature">
<div id="gmail-m_5681557970449129721gmail-m_8275223278710308024gmail-m_4073407159693916210gmail-m_8546846971254182096gmail-m_2454870331033590078divtagdefaultwrapper" dir="ltr" style="font-size:12pt;color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif,Helvetica,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols">
<div name="divtagdefaultwrapper">
<div name="divtagdefaultwrapper"><font size="3" face="'Times New Roman', Times, serif" color="808080">--<br>
Dr. Iurii Timrov<br>
Postdoctoral Researcher<br>
<font color="808080"><font face="'Times New Roman', Times, serif"></font></font></font></div>
<font color="808080"></font>
<div name="divtagdefaultwrapper"><font size="3" face="'Times New Roman', Times, serif" color="808080">STI - IMX
<font color="808080"><font face="'Times New Roman', Times, serif">- THEOS</font></font></font><font size="3" face="'Times New Roman', Times, serif" color="808080"> and NCCR - MARVEL<br>
</font></div>
<div name="divtagdefaultwrapper"><font size="3" face="'Times New Roman', Times, serif" color="808080"><font size="3" face="'Times New Roman', Times, serif" color="808080">Swiss Federal Institute of Technology Lausanne (EPFL<font color="808080"><font face="'Times New Roman', Times, serif">)</font></font></font><br>
</font></div>
<font color="808080"></font>
<div name="divtagdefaultwrapper"><font size="3" face="'Times New Roman', Times, serif" color="808080">CH-1015 Lausanne, Switzerland<br>
+41 21 69 34 881</font></div>
<div name="divtagdefaultwrapper"><a href="http://people.epfl.ch/265334" id="gmail-m_5681557970449129721gmail-m_8275223278710308024gmail-m_4073407159693916210gmail-m_8546846971254182096gmail-m_2454870331033590078LPNoLP" target="_blank">http://people.epfl.ch/265334</a><br>
</div>
</div>
</div>
</div>
</div>
<hr style="display:inline-block;width:98%">
<div id="gmail-m_5681557970449129721gmail-m_8275223278710308024gmail-m_4073407159693916210gmail-m_8546846971254182096gmail-m_2454870331033590078divRplyFwdMsg" dir="ltr">
<font style="font-size:11pt" face="Calibri, sans-serif" color="#000000"><b>From:</b> users <<a href="mailto:users-bounces@lists.quantum-espresso.org" target="_blank">users-bounces@lists.quantum-espresso.org</a>> on behalf of Merlin Meheut <<a href="mailto:merlin.meheut@gmail.com" target="_blank">merlin.meheut@gmail.com</a>><br>
<b>Sent:</b> Sunday, February 16, 2020 9:22:57 PM<br>
<b>To:</b> Quantum Espresso users Forum<br>
<b>Subject:</b> [QE-users] non-convergency of DFPT+U calculation</font>
<div> </div>
</div>
<div>
<div dir="ltr">
<div>Dear PWSCF users, <br>
</div>
<div><br>
</div>
<div>I am trying to realize a phonon calculation (at gamma) with +U on fayalite (olivine structure Fe2SiO4) with anti-ferromagnetic configuration , and whether the calculation goes fine without +U, I have one representation (#23) that systematically fails to
converge when I add a +U. I tried several values for the +U (0.5, 2.5, 4), several functionals (PZ, PBE), but the result is always the same. I have used QE versions 6.4 and 6.4.1 on two different computing centers . Would you have any idea about the reasons
for such a failure? I am specifically worried by one output message: "Atomic wfc used for the DFT+U projector are NOT orthogonalized", but I don't know how to fix that.
<br>
</div>
<div><br>
</div>
<div>Here are examples of my scf and ph inputs, and the failing ph output: <br>
</div>
<div><br>
</div>
<div>scf:<br>
&control<br>
calculation = 'scf',<br>
(...)<br>
tprnfor = .true.,<br>
tstress = .true.,<br>
/&end<br>
&system<br>
ibrav =0, celldm(1)=9.107535, <br>
nat =28, ntyp = 4, ecutwfc =80.0, ecutrho = 320.0,<br>
nspin=2, tot_magnetization=0,<br>
starting_magnetization(1)=0.5, starting_magnetization(2)=-0.5,<br>
lda_plus_u = .true., U_projection_type = 'atomic',<br>
Hubbard_U(1)=2.4, Hubbard_U(2)=2.6<br>
/&end <br>
&electrons<br>
electron_maxstep = 150,<br>
conv_thr = 1.d-11,<br>
mixing_mode = 'plain',<br>
startingwfc = 'atomic',<br>
mixing_beta = 0.3,<br>
/&end<br>
ATOMIC_SPECIES<br>
Fe1 55.9349 Fe.pz-n-rrkjus_psl.1.0.0.UPF <br>
Fe2 55.9349 Fe.pz-n-rrkjus_psl.1.0.0.UPF<br>
Si 27.9769 Si.pz2.UPF<br>
O 15.9949 O.pz.UPF<br>
(...)</div>
<div><br>
</div>
<div>ph input: <br>
</div>
<div> &inputph<br>
amass(1)=55.9349,<br>
amass(2)=55.9349,<br>
amass(3)=27.9769,<br>
amass(4)=15.9949,<br>
alpha_mix(1) = 0.5,<br>
! ldisp=.true., nq1=2, nq2=2, nq3=2,<br>
tr2_ph = 1.0D-17,<br>
prefix='FAYA2',<br>
fildyn='mat.$PREFIX',<br>
lraman=.false.,<br>
epsil =.false.,<br>
trans =.true.,<br>
zue = .false.,<br>
zeu= .false.,<br>
start_irr=23,<br>
last_irr=23, <br>
outdir ='./',<br>
/&end<br>
0.0 0.0 0.0<br>
</div>
<div><br>
</div>
<div>ph output: <br>
</div>
<div> Representation # 23 mode # 23<br>
<br>
Self-consistent Calculation<br>
<br>
iter # 1 total cpu time : 58.7 secs <a href="http://av.it" target="_blank">
av.it</a>.: 8.5<br>
thresh= 1.000E-02 alpha_mix = 0.500 |ddv_scf|^2 = 1.045E-06<br>
<br>
iter # 2 total cpu time : 62.9 secs <a href="http://av.it" target="_blank">
av.it</a>.: 23.0<br>
thresh= 1.022E-04 alpha_mix = 0.500 |ddv_scf|^2 = 2.668E-04<br>
<br>
iter # 3 total cpu time : 66.4 secs <a href="http://av.it" target="_blank">
av.it</a>.: 17.5<br>
thresh= 1.633E-03 alpha_mix = 0.500 |ddv_scf|^2 = 3.033E-03<br>
<br>
iter # 4 total cpu time : 69.0 secs <a href="http://av.it" target="_blank">
av.it</a>.: 11.0<br>
thresh= 5.508E-03 alpha_mix = 0.500 |ddv_scf|^2 = 5.602E-02<br>
<br>
iter # 5 total cpu time : 73.0 secs <a href="http://av.it" target="_blank">
av.it</a>.: 15.5<br>
thresh= 1.000E-02 alpha_mix = 0.500 |ddv_scf|^2 = 1.238E+00<br>
<br>
iter # 6 total cpu time : 76.7 secs <a href="http://av.it" target="_blank">
av.it</a>.: 19.5<br>
thresh= 1.000E-02 alpha_mix = 0.500 |ddv_scf|^2 = 2.841E+01<br>
<br>
iter # 7 total cpu time : 80.9 secs <a href="http://av.it" target="_blank">
av.it</a>.: 23.0<br>
thresh= 1.000E-02 alpha_mix = 0.500 |ddv_scf|^2 = 6.534E+02<br>
<br>
iter # 8 total cpu time : 85.3 secs <a href="http://av.it" target="_blank">
av.it</a>.: 24.5<br>
thresh= 1.000E-02 alpha_mix = 0.500 |ddv_scf|^2 = 1.446E+04<br>
<br>
iter # 9 total cpu time : 90.1 secs <a href="http://av.it" target="_blank">
av.it</a>.: 27.5<br>
thresh= 1.000E-02 alpha_mix = 0.500 |ddv_scf|^2 = 3.214E+05<br>
<br>
iter # 10 total cpu time : 95.5 secs <a href="http://av.it" target="_blank">
av.it</a>.: 30.0<br>
thresh= 1.000E-02 alpha_mix = 0.500 |ddv_scf|^2 = 7.200E+06</div>
<div><br>
</div>
<div>(.....)<br>
</div>
<div><br>
</div>
<div>Thank you in advance for your help! Feel free to ask for any further calculation details.<br>
</div>
<div><br>
</div>
<div>Regards, <br>
</div>
<div>-- <br>
<div dir="ltr">
<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: (+33) 5 61 33 26 17<br>
<br>
</div>
</div>
</div>
</div>
</div>
</div>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso" rel="noreferrer" target="_blank">www.max-centre.eu/quantum-espresso</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">
users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote>
</div>
<br clear="all">
<br>
-- <br>
<div dir="ltr">
<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: (+33) 5 61 33 26 17<br>
<br>
</div>
</div>
</div>
</div>
</div>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso" rel="noreferrer" target="_blank">www.max-centre.eu/quantum-espresso</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">
users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote>
</div>
<br clear="all">
<br>
-- <br>
<div dir="ltr">
<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: (+33) 5 61 33 26 17<br>
<br>
</div>
</div>
</div>
</div>
</div>
</div>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu/quantum-espresso" rel="noreferrer" target="_blank">www.max-centre.eu/quantum-espresso</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div><br clear="all"><br><div dir="ltr"><div dir="ltr"><br></div></div>
</div>