<div dir="ltr">Thanks a lot Matt for your suggestions and additional clarification.<br><br>On Monday, 10 July 2017 19:35:44 UTC+2, Matt W  wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr"><blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><br></div><div>Thanks for your suggestions. After playing around with the parameters, I figured the Extrapolation issue and cut-off issue (so fixed that). </div><div>I already see improvement in time after some trials. I am wondering though which method would be best for extrapolation: ASPC etc.</div></div></blockquote><div><br></div><div>I'd stick with ASPC. PS is maybe a bit faster for geo_opt, but ASPC should be good and robust.</div><div> </div><blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>I haven't switched to GAPW though. I will try that. </div></div></blockquote><div><br></div><div>Only needed really when you have some elements that need stupidly large cutoff (F, Na, Mg) with normal pseudos. It does tend to make forces a bit more accurate, so can speed up final convergence a lot too.</div><div> </div><blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>Can you please elaborate on "Mixing and broyden optimizer will be moderately inefficient for insulators"</div><div>and a possible alternative.</div><div><br></div></div></blockquote><div>The Orbital Transform (OT) method is CP2Ks flagship optimizer for insulators. For a smallish system</div><div><br></div><div><div>    &SCF</div><div>      &OUTER_SCF</div><div>         EPS_SCF 5.0E-7</div><div>         MAX_SCF 1</div><div>      &END OUTER_SCF</div><div>      SCF_GUESS RESTART</div><div>      EPS_SCF 5.0E-7</div><div>      MAX_SCF  15</div><div>      &OT</div><div>        MINIMIZER DIIS</div><div>        PRECONDITIONER FULL_ALL</div><div>        ENERGY_GAP 0.001</div><div>      &END</div><div>    &END SCF</div></div><div><br></div><div>would be my default setup. </div><div><br></div><div><div>        PRECONDITIONER FULL_SINGLE_INVERSE</div><div>        ENERGY_GAP 0.1</div></div><div><br></div><div>for large systems (>500 atoms or so)</div><div><br></div><div>Matt</div><div><br></div><div> </div><blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div></div><div>Thanks</div><div><div><blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><br></div></div></blockquote></div></div></blockquote></div></div></div></blockquote></div></blockquote></div>