<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 data-smartmail="gmail_signature" dir="auto">--<br>Professor Laurence Marks<br>Department of Materials Science and Engineering<br>Northwestern University<br><a href="http://www.numis.northwestern.edu" target="_blank" rel="noreferrer">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>