[CP2K-user] [CP2K:21497] System exploded in NPT-F simulations
Zhuoran Long
zlong1 at umd.edu
Fri May 23 20:38:21 UTC 2025
Dear all,
I am testing NPT simulations on the tutorial water system
<https://www.cp2k.org/exercises:2018_ethz_mmm:h2o_md>. I'm including the
following changes to the input file in the tutorial example:
&FORCE_EVAL
STRESS_TENSOR ANALYTICAL
...
&END FORCE_EVAL
...
&MOTION
&MD
ENSEMBLE NPT_I
STEPS 200000
TIMESTEP 1.0
TEMPERATURE 300
COMVEL_TOL 1.0E-10
&BAROSTAT
PRESSURE 1.01325 #[bar]
TIMECON 1000
VIRIAL XYZ
&THERMOSTAT
TYPE SAME_AS_PARTICLE
&CSVR
TIMECON 200
&END CSVR
&END THERMOSTAT
&END BAROSTAT
&THERMOSTAT
TYPE CSVR
REGION MASSIVE
&CSVR
TIMECON 200
&END CSVR
&END THERMOSTAT
...
&END MD
...
&END MOTION
The Cons Qty in the .ener file converged after about 20 ps. However. when I
perform the same simulation from scratch in the NPT-F ensemble (by changing
the NPT-I to NPT-F), the Cons Qty kept increasing and exploded after about
80 ps (please see the plot below, so are the kinetic energy and potential
energy). I will need the NPT-F in future for a tetragonal covalent organic
framework system. Could someone please let me know if I had any wrong or
missing parameters?
And a relevant question, while I set the time constant to be 200 for the
CSVR thermostat of the barostat, the corresponding TIMECON in the restart
file has a value of ~9.6755. Is that relevant to my question? I'm using the
version 2024.2.
And helps and discussions are greatly appreciated.
Zhuoran
[image: t300.png]
--
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 visit https://groups.google.com/d/msgid/cp2k/542e30ec-8ae8-4b3b-bc8c-697214641c81n%40googlegroups.com.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20250523/4b99c2e8/attachment-0001.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: t300.png
Type: image/png
Size: 45051 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20250523/4b99c2e8/attachment-0001.png>
More information about the CP2K-user
mailing list