<div dir="ltr">Thank you 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 a 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>