CP2K Open Source Molecular Dynamics

User Tools

Site Tools

Trace:
exercises:common:ensemble

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revisionPrevious revision
Next revision		Previous revision
exercises:common:ensemble [2022/10/31 16:15] – jglanexercises:common:ensemble [2022/10/31 17:28] (current) – [NVT Ensemble] jglan
Line 217: Line 217:
 ===== NVT Ensemble  ===== ===== NVT Ensemble  =====
  
-In the previous sections, you have already run NVE ensemble molecular dynamics calculations for liquid Ar. In this section, we will focus on the NVT and NPT ensembles.+In the previous sections, you have already run NVE ensemble molecular dynamics calculations for liquid Ar. In this section, we will focus on the NVT ensembles.  
 + 
 +Although the most popular Nose-Hoover thermostat is commonly-used in other MD codes, the original Nose-Hoover thermostat has an ergodic issue. This has been solved by Mark Tuckerman et al. See [[https://doi.org/10.1063/1.463940 |Nose-Hoover Chain Thermostat]]  
 +In CP2K, the default length of the Nose-Hoover chain is set to 3. (See Manual [[https://manual.cp2k.org/cp2k-2022_1-branch/CP2K_INPUT/MOTION/MD/THERMOSTAT/NOSE.html| %MOTION%MD%THERMOSTAT%NOSE]]) 
  
 To set up an NVT calculation, change the settings in the &MD section as shown below: To set up an NVT calculation, change the settings in the &MD section as shown below:
Line 230: Line 234:
       REGION MASSIVE       REGION MASSIVE
       &NOSE                    #Uses the Nose-Hoover thermostat       &NOSE                    #Uses the Nose-Hoover thermostat
-        TIMECON 1000           #timeconstant of the thermostat chain, how often does thermostat adjust your system +        TIMECON 100           #timeconstant of the thermostat chain, how often does the thermostat adjust your system 
       &END NOSE       &END NOSE
     &END     &END
Line 236: Line 240:
      
  
 +Alternatively, one can also use [[https://manual.cp2k.org/cp2k-2022_1-branch/CP2K_INPUT/MOTION/MD/THERMOSTAT/CSVR.html| Canonical sampling through velocity rescaling
 +(CSVR)]] as developed by Giovanni [[https://doi.org/10.1063/1.2408420 |Bussi et al. ]]
 +  &MD
 +    ENSEMBLE NVT
 +    STEPS 3000
 +    TIMESTEP 5
 +    TEMPERATURE 298
 +    &THERMOSTAT
 +      &CSVR                    
 +        TIMECON 100           #timeconstant of the CSVR, how often does the thermostat adjust your system 
 +      &END CSVR
 +    &END
 +  &END MD
 +  
 ===== NPT Ensemble  ===== ===== NPT Ensemble  =====
  
exercises/common/ensemble.1667232938.txt.gz · Last modified: 2022/10/31 16:15 by jglan