<div dir="ltr">Dear Prof. Marks,<div><br></div><div>                           Thank you for your explanations Sir. Comments are very informative.</div><div><br></div><div>with regards,</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, 6 Oct 2021 at 18:30, Laurence Marks <<a href="mailto:laurence.marks@gmail.com">laurence.marks@gmail.com</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 dir="auto">Some mixer and optimization comments, to clarify.<div dir="auto"><br></div><div dir="auto">99.999% of cases where there are convergence problems are due to:</div><div dir="auto"><br></div><div dir="auto">1) Inappropriate physics. A bad model does not have to have a fixed-point solution, that is there may be no converged solution. (Nothing says there has to be.)</div><div dir="auto"><br></div><div dir="auto">2) Inappropriate numerics. For instance, a functional that jumps around very noisily may not be smooth. I know this can happen for 4f electrons using -eece.</div><div dir="auto"><br></div><div dir="auto">3) Discontinuity at the solution. There is such a discontinuity (electronic phase transition) for systems such as isolated Fe atoms without spin, T=0.</div><div dir="auto"><br></div><div dir="auto">4) User error. I recently ran a case with 16 symmetry elements which only should have had 8. While it came close, it would not converge.</div><div dir="auto"><br></div><div dir="auto">Sometimes the pseudo-forces are quite bad with large density residues, so MSR1a can be slow. There can also be issues if the electron-phonon coupling is large.</div><div dir="auto"><br></div><div dir="auto">PORT is sometimes better, particularly with repetitive calculations such as lattice parameter optimizations. However, it is very sensitive to inconsistent forces and energies.</div><div dir="auto"><br></div><div dir="auto">In my opinion NEWT & NEWT1 should be removed. My opinion.</div><div dir="auto"><br></div><div dir="auto">Always look first at how the model is setup.</div><div dir="auto"><br><div dir="auto">--<br>Professor Laurence Marks<br>Department of Materials Science and Engineering<br>Northwestern University<br><a href="http://www.numis.northwestern.edu" rel="noreferrer" target="_blank">www.numis.northwestern.edu</a><br>"Research is to see what everybody else has seen, and to think what<br>nobody else has thought" Albert Szent-Györgyi</div></div></div>
_______________________________________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
SEARCH the MAILING-LIST at:  <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" rel="noreferrer" target="_blank">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div style="font-size:small">Dr. Shamik Chakrabarti</div><div style="font-size:small">Research Fellow </div><div style="font-size:small">Department of Physics</div><div style="font-size:small">Indian Institute of Technology Patna</div><div style="font-size:small">Bihta-801103</div><div style="font-size:small">Patna</div><div style="font-size:small">Bihar, India</div></div></div></div></div></div></div></div>