[CP2K-user] [CP2K:18524] PW cutoff determination for aqueous system
Rajorshi Chattopadhyay
rajorshichat at gmail.com
Wed Mar 8 11:08:35 UTC 2023
Hello all,
I am trying to run AIMD simulations with a box containing water and La+3
ion (net charge +3). I have done an energy convergence test as mentioned in
a previous discussion. The results are as follows:
Cutoff
Rel cutoff
Energy
Charge
200
200
-1479.12144663111371
3.0001076208
300
300
-1479.04319059475324
3.0000001878
400
400
-1479.02723340271655
2.9999999983
500
500
-1479.02438781296655
2.9999999495
600
600
-1479.02417065297800
2.9999999885
700
700
-1479.02386062640312
2.9999999880
800
800
-1479.02394157188951
2.9999999934
900
900
-1479.02387722893354
2.9999999825
1000
1000
-1479.02390064546512
2.9999999922
Cutoff
Rel cutoff
Energy
Charge
1000
30
-1479.02482154760651
2.9999999769
40
-1479.02395654562406
2.9999999781
60
-1479.02384588278346
2.9999999882
80
-1479.02392425354356
2.9999999783
100
-1479.02388590916007
2.9999999797
500
-1479.02378029727311
2.9999999871
1000
-1479.023810336304450
2.9999999969
>From the data, it seems that 700,60 should be a good choice for production
calculations.
However, I have two questions:
1. The difference between the charge of the system and the total charge on
the grids is more than 10**-8. Is this ok ?
2. I see energy convergence only till 3 decimal places. Is this tight
enough ?
I have attached the input script and initial coordinates with this email. I
thank you in advance for your suggestions/help.
--
With Regards,
Rajorshi Chattopadhyay,
PhD Researcher in Mineralogy/Crystallography,
Institut für Geologie und Mineralogie,
Universität zu Köln, Germany
--
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 cp2k+unsubscribe at googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/cp2k/CACTJSRcNGVQm1FSpa7mEE28%2BfsnHesTEaC6p7aw9uy9%3DrCwemw%40mail.gmail.com.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20230308/ea8b1588/attachment-0001.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: main.inp
Type: application/octet-stream
Size: 1613 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20230308/ea8b1588/attachment-0002.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: coord_eql
Type: application/octet-stream
Size: 16445 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20230308/ea8b1588/attachment-0003.obj>
More information about the CP2K-user
mailing list