[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