<div dir="ltr"><div>Dear Nick,</div><div><br></div><div>I typically use something like this (for HfO2, TiO2 etc.)<br></div><div><br></div><div>      EPS_SCF 1.0E-6<br>      MAX_SCF 30<br>      &OT<br>        PRECONDITIONER FULL_SINGLE_INVERSE<br>        ALGORITHM IRAC<br>        ENERGY_GAP 0.5<br>        MINIMIZER DIIS<br>      &END OT<br>      &OUTER_SCF<br>         EPS_SCF 1.0E-6<br>         MAX_SCF 20<br>      &END OUTER_SCF</div><div><br></div><div>In my experience this is a good compromise between accuracy and speed, e.g. for MD. If you need more accurate forces, e.g. CELL_OPT comment the following two lines out<br></div><div>!        ALGORITHM IRAC<br>!        ENERGY_GAP 0.5</div><div><br></div><div>You should compare this to GC to make sure this works for your system.</div><div>      EPS_SCF 1.0E-6<br>      MAX_SCF 30<br>      &OT<br>        PRECONDITIONER FULL_SINGLE_INVERSE<br>        MINIMIZER GC</div><div>        LINESEARCH 3PNT<br>      &END OT<br>      &OUTER_SCF<br>         EPS_SCF 1.0E-6<br>         MAX_SCF 20<br>      &END OUTER_SCF</div><div><br></div><div>I hope this helps,</div><div>Fabian<br></div></div>