<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">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">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">si.scf.efield2.in</a> after running <a href="http://si.scf.efield.in">si.scf.efield.in</a>)<br></div><div>-## INPUT1: <a href="http://si.scf.efield.in">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">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">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>