[CP2K-user] [CP2K:13012] AIMD of ACN doesn't converge after 100 ps
Thomas Kühne
tku... at gmail.com
Tue Mar 24 14:02:18 UTC 2020
Dear Abdullah Bin Faheem,
there is nothing wrong. What you are observing is not due to incomplete equilibration, but
rather a consequence, of incomplete SCF convergence. Beside time-step errors, which in
your case may be non-negligible due to a discretized time-step of 1.0 fs, the so-called
„conserved quantity“ is in a numerical calculation only conserved at perfect SCF convergence.
Hence, reducing the integration time-step and more importantly tightening the EPS_SCF
should lead to a better conservation of the conserved quantity. However, in a numerical
calculation the latter will never be perfectly conserved and it is up to you to decided what
level gives you the best efficiency / accuracy compromise.
I personally, would suggest 0.5-0.75fs for the time-step and 1e-7 - 1e-6 for EPS_SCF when
using OT. I have commented on the differences between TD and OT wrt to EPS_SCF in an
earlier posting.
Cheers,
Thomas
> Am 24.03.2020 um 14:03 schrieb Abdullah Bin Faheem <abdullahb... at gmail.com>:
>
> Hello everyone,
>
> I performed AIMD of 60 Acetonitrile molecules using cp2k 5.1.0. However, I noticed that "Cons Qty" keeps increasing and does not converge even after 100 ps. While on the other hand, the kinetic and potential energies reach equilibrium within 2-3 ps. Does this mean my system has not converged? Or is there something wrong with my input file?
>
> If so, any suggestions would be greatly appreciated for possible improvement/fixes.
>
> I have attached the input and the energy output files, where the initial geometry for ACN was taking from classical MD.
>
> Thank you for your time
>
> Best Regards,
> Abdullah Bin Faheem
>
> <Energies.png> <about:invalid#zClosurez>
>
>
> --
> You received this message because you are subscribed to the Google Groups "cp2k" group.
> To unsubscribe from this group and stop receiving emails from it, send an email to cp... at googlegroups.com <mailto:cp... at googlegroups.com>.
> To view this discussion on the web visit https://groups.google.com/d/msgid/cp2k/71cf94c9-77b2-4561-8af7-cbf9d428f828%40googlegroups.com <https://groups.google.com/d/msgid/cp2k/71cf94c9-77b2-4561-8af7-cbf9d428f828%40googlegroups.com?utm_medium=email&utm_source=footer>.
> <files.rar><Energies.png>
==============================
Thomas D. Kühne
Dynamics of Condensed Matter
Chair of Theoretical Chemistry
University of Paderborn
Warburger Str. 100
D-33098 Paderborn
Germany
tdku... at mail.upb.de
+49/(0)5251/60-5726
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20200324/87fc9cc2/attachment.htm>
More information about the CP2K-user
mailing list