Cut-off problrm
simin pahlavi
simin.... at gmail.com
Tue Apr 24 11:16:54 UTC 2018
Hello Matteu,
I would like to thank you very much for your Helpful comments again.
Have a nice day
Simin
On Tuesday, April 24, 2018 at 12:42:37 PM UTC+2, Matt W wrote:
>
> 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.
>
> I would look at forces on atoms and see how they converge and see how
> errors in forces propagate to energy conservation.
>
> Matt
>
> On Tuesday, April 24, 2018 at 10:24:34 AM UTC+1, simin pahlavi wrote:
>>
>> Hello Matteu
>> Thank you very much fir your helpful comment.
>> As I am supposed to use this DFT calculations for every step of an ab
>> initio molecular dynamics, I wonder if I should
>> go for the heavy smoothing and get reasonable forces or I should be worry
>> about convergence and cutoff dependence of minimization?
>> I appreciate any reply in advance
>> Simin
>>
>> On Monday, April 23, 2018 at 10:41:00 PM UTC+2, Matt W wrote:
>>>
>>> Sorry, I posted too quickly without really looking at your setup.
>>>
>>> The big problem for you is the XC_GRID settings.
>>>
>>> &XC_GRID
>>> XC_SMOOTH_RHO NN50
>>> XC_DERIV NN50_SMOOTH
>>> &END XC_GRID
>>>
>>> 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.
>>>
>>> Try changing to the defaults
>>>
>>> &XC_GRID
>>> XC_SMOOTH_RHO NONE
>>> XC_DERIV PW
>>> &END XC_GRID
>>>
>>> you should see at least a significant improvement.
>>>
>>> Matt
>>>
>>>>
>>>>>>>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20180424/08c5e19c/attachment.htm>
More information about the CP2K-user
mailing list