<div dir="ltr">Hi Lorenzo,<div><br></div><div>This is really an important piece of information. I'll give it a try.</div><div><br></div><div>Thank you,</div><div>Boyang</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Aug 30, 2024 at 11:46 AM Lorenzo Bastonero <<a href="mailto:lbastone@uni-bremen.de">lbastone@uni-bremen.de</a>> wrote:<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>Dear Boyang,<div><br></div><div>Something that I noticed in the past is that for some QE versions the electric field calculations did not converge when compiling with “-O3” flag. Using “-O2” solved the problem for me. I don’t know if this might apply here as well, but maybe it’s worth trying.</div><div><br></div><div>HTH,</div><div>Lorenzo</div><div><br><div>
<div dir="auto" style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none"><div dir="auto" style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none"><div dir="auto" style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none"><div dir="auto" style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none"><div dir="auto" style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;text-decoration:none"><div><div>********************************</div><div>Lorenzo Bastonero</div><div><br></div><div>PhD Student</div><div>U Bremen Excellence Chair,<br>Bremen Center for Computational Materials Science,<br>and MAPEX Center for Materials and Processes</div><div><br></div><div>University of Bremen</div><div>Faculty of Production Engineering</div><div>TAB-Building, Room 3.32</div><div>Am Fallturm 1</div><div>28359 Bremen, DE</div><br><a href="http://www.hmi.uni-bremen.de/" target="_blank">http://www.hmi.uni-bremen.de/</a></div></div></div></div></div></div>
</div>
<div><br><blockquote type="cite"><div>Il giorno 30 ago 2024, alle ore 17:31, Boyang Zheng <<a href="mailto:boyang@psu.edu" target="_blank">boyang@psu.edu</a>> ha scritto:</div><br><div><div dir="ltr">Dear Prof. Giannozzi,<div><br></div><div>Thanks for your reply. </div><div><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">if you use the exact same input, it works.</blockquote><div>I really hope it works, but it does not. I only deleted the "pseudo_dir" and "outdir" lines, which I think should have no effect on the calculation.</div><div><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">You may need to perform a zero-field calculation first, then do the finite-field one with initial charge density and wavefunctions from the zero-field case.</blockquote><div>I tried, still did not converge.</div><div> </div><div><br></div><div>To make things super clear, here (see the gist link below) I give my inputs and outputs of 2 tries. You can navigate by searching "-#" in there. A summary of what I did is also copied here.</div><div><a href="https://gist.github.com/boyoungzheng/c9e653b1b07766482ae46ca2826b5940" target="_blank">https://gist.github.com/boyoungzheng/c9e653b1b07766482ae46ca2826b5940</a><br></div><div><br></div><div>-# 1st try: Same input from "run_example" (i.e., restart_mode='from_scratch')<br>-## INPUT: <a href="http://si.scf.efield2.in/" target="_blank">si.scf.efield2.in</a><br></div><div>-## OUTPUT: si.scf.efield2.out  (not converged)</div><div>-# 2nd try: Restart from converged 0-field calculation (i.e., set restart_mode='restart' in <a href="http://si.scf.efield2.in/" target="_blank">si.scf.efield2.in</a> after running <a href="http://si.scf.efield.in/" target="_blank">si.scf.efield.in</a>)<br></div><div>-## INPUT1: <a href="http://si.scf.efield.in/" target="_blank">si.scf.efield.in</a><br></div><div>-## OUTPUT1: si.scf.efield.out  (converged)<br></div><div>-## INPUT2: <a href="http://si.scf.efield2.in/" target="_blank">si.scf.efield2.in</a><br></div><div>-## OUTPUT2: si.scf.efield2.out  (not converged)<br></div><div><br></div><div>Thanks,</div><div>Boyang</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Aug 30, 2024 at 6:19 AM Paolo Giannozzi <<a href="mailto:paolo.giannozzi@uniud.it" target="_blank">paolo.giannozzi@uniud.it</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">On 8/28/24 21:07, Boyang Zheng wrote:<br>
<br>
> I am running PW/examples/example10. I used the exact same input <br>
> "<a href="http://si.scf.efield2.in/" rel="noreferrer" target="_blank">si.scf.efield2.in</a>" from "run_example" (...). But the <br>
> calculation shows "convergence NOT achieved after 100 iterations".<br>
<br>
if you use the exact same input, it works.<br>
<br>
Finite-field calculations are tricky. You may need to perform a <br>
zero-field calculation first, then do the finite-field one with initial <br>
charge density and wavefunctions from the zero-field case.<br>
<br>
Paolo<br>
-- <br>
Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>
Univ. Udine, via delle Scienze 208, 33100 Udine Italy, +39-0432-558216<br>
</blockquote></div>
_______________________________________________<br>The Quantum ESPRESSO community stands by the Ukrainian<br>people and expresses its concerns about the devastating<br>effects that the Russian military offensive has on their<br>country and on the free and peaceful scientific, cultural,<br>and economic cooperation amongst peoples<br>_______________________________________________<br>Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" target="_blank">www.max-centre.eu</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" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></div></blockquote></div><br></div></div>_______________________________________________<br>
The Quantum ESPRESSO community stands by the Ukrainian<br>
people and expresses its concerns about the devastating<br>
effects that the Russian military offensive has on their<br>
country and on the free and peaceful scientific, cultural,<br>
and economic cooperation amongst peoples<br>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" rel="noreferrer" target="_blank">www.max-centre.eu</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>