<div dir="ltr">Dear Nicola,<div>I am looking forward to reading the "short thing on this"! </div><div>Can you comment on the widespread (?) practice of using much lower than recommended (say, 2/3 or 1/2) ecutwfc for geometry optimization until a minimum is found, and then increasing cutoffs until convergence is achieved? Is it generally a good idea?</div><div><br></div><div>Best regards,</div><div><br></div><div>Michal Krompiec</div><div><br></div><div>Associate Director at Merck KGaA and Adjunct Prof at University of Southampton</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, 19 May 2019 at 13:57, Nicola Marzari <<a href="mailto:nicola.marzari@epfl.ch">nicola.marzari@epfl.ch</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"><br>
<br>
Hi Robert,<br>
<br>
this issue of optimal parameters is a bit of a neverending story, we are<br>
preparing a short thing on this. But in a nutshell - yes, the electronic <br>
convergence should be tighter, especially for a follow-up phonon <br>
calculation. that number actually depends on the unit cell used, but for <br>
5 atoms, I would indeed try d-12 and d-14. For ecutwfc and ecutrho, <br>
well, my advise is to use the SSSP suggestions (google it); and say if <br>
you have 30 120, 40 160, 25 300, well, take the max of the ecutwfc <br>
suggested, for the gloval ecutwfc, and same for ecutrho.<br>
<br>
                        nic<br>
<br>
<br>
On 17/05/2019 21:24, Appleton, Robert J wrote:<br>
> Hello,<br>
> <br>
> I am trying to study the phonons of the material BaZrS3 to determine if <br>
> it is stable for a given structure. I first use phonopy and vasp and <br>
> found negative frequency values for phoning band structure. My advisor <br>
> is concerned in the optimization of the structure and so we now want to <br>
> use quantum espresso. I have only used it a few times and mostly used <br>
> default or example inputs. The input I want to make sure is sufficient <br>
> is my criteria for convergence.<br>
> <br>
> Under the &control:<br>
> force convergence criteria right now is<br>
> forc_conv_thr = .0001 <— this value seems to big but how do I choose a <br>
> value that is best for optimization but does not take to much computer time?<br>
> <br>
> Under the &electrons:<br>
> conv_thr = 1.d-9 <— I am considering changing to 1.d-12. Is this a good <br>
> idea?<br>
> <br>
> Under &system:<br>
> I am using ecutrho = 12 • ecutwfc<br>
> <br>
> If you can give me some advice on how best to choose convergence <br>
> criteria for studying phonons that would be very good thanks.<br>
> <br>
> Robert Appleton<br>
> Cal State LA - MS Physics Student<br>
> <br>
> Get Outlook for iOS <<a href="https://aka.ms/o0ukef" rel="noreferrer" target="_blank">https://aka.ms/o0ukef</a>><br>
> <br>
> _______________________________________________<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><br>
> <br>
<br>
<br>
-- <br>
----------------------------------------------------------------------<br>
Prof Nicola Marzari, Chair of Theory and Simulation of Materials, EPFL<br>
Director, National Centre for Competence in Research NCCR MARVEL, EPFL<br>
<a href="http://theossrv1.epfl.ch/Main/Contact" rel="noreferrer" target="_blank">http://theossrv1.epfl.ch/Main/Contact</a> <a href="http://nccr-marvel.ch/en/project" rel="noreferrer" target="_blank">http://nccr-marvel.ch/en/project</a><br>
_______________________________________________<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><br>
</blockquote></div>