<div dir="ltr">Hello Matteu,<br>I would like to thank you very much for your Helpful comments again.<div>Have a nice day</div><div>Simin</div><div><br>On Tuesday, April 24, 2018 at 12:42:37 PM 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">Personal opinion is that converging total energy is not very important in this case - it is normally just a rigid shift based on inaccuracy of the core part of the wavefunctions. I am still (again personally) not really a big fan of the smoothing, I'd rather have a slightly bigger cutoff but a not messed with functional.<div><br></div><div>I would look at forces on atoms and see how they converge and see how errors in forces propagate to energy conservation.</div><div><br></div><div>Matt<br><br>On Tuesday, April 24, 2018 at 10:24:34 AM UTC+1, simin pahlavi wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hello Matteu <div>Thank you very much fir your helpful comment.</div><div>As I am supposed to use this DFT calculations for every step of an ab initio molecular dynamics, I wonder if I should</div><div>go for the heavy smoothing and get reasonable forces or I should be worry about convergence and cutoff dependence of minimization?</div><div>I appreciate any reply in advance</div><div>Simin<br><br>On Monday, April 23, 2018 at 10:41:00 PM 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">Sorry, I posted too quickly without really looking at your setup.<div><br></div><div>The big problem for you is the XC_GRID settings. </div><div><br></div><div><div>   &XC_GRID</div><div>    XC_SMOOTH_RHO NN50</div><div>    XC_DERIV NN50_SMOOTH</div><div>   &END XC_GRID</div></div><div><br></div><div>You really do not want to be using heavy smoothing, especially for this purpose. These settings (see original paper) help get reasonable forces with a rather low cutoff, but as you increase the cutoff lead to you continuously redefining the XC functional you use (you are smoothing over different areas of real space), so they do not converge.</div><div><br></div><div>Try changing to the defaults</div><div><br></div><div><div><div>   &XC_GRID</div><div>    XC_SMOOTH_RHO NONE</div><div>    XC_DERIV PW</div><div>   &END XC_GRID</div></div></div><div><br></div><div>you should see at least a significant improvement.</div><div><br></div><div>Matt<br><blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><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"><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></blockquote></div></blockquote></div></div></div></blockquote></div></div></blockquote></div></div></blockquote></div></div></blockquote></div></div>