NVT ensemble: Difference between revisions

Revision as of 13:19, 19 July 2022

NVT ensemble

The NVT ensemble is a statistical ensemble that is used to study material properties under the conditions of a constant particle number N, constant volume V and a temperature fluctuating around an equilibrium value T. This page describes how to sample the NVT ensemble from a molecular dynamics run.

Instructions for setting up a NVT ensemble There are 4 choices of thermostats which can either be stochastic or deterministic to simulate the NVT ensemble. The stochastic Andersen thermostat or Langevin thermostat, the deterministic Nose-Hoover thermostat or Multiple Anderson thermostats can be used. See table for the corresponding MDALGO tags.

NVT ensemble	Andersen	Langevin	Nose-Hoover	Multiple Andersen
MDALGO	1	3	2	13
additional tags	ANDERSEN_PROB	LANGEVIN_GAMMA	SMASS	NSUBSYS, TSUBSYS, PSUBSYS

The user is not allowed to set ISF>2 because then the volume of the simulation would be allowed to change and the simulation would not result in a NVT ensemble. Other flags related to molecular dynamics simulations can be found here.

A simple example INCAR file for the Langevin thermostat could look like'

 #INCAR molecular dynamics tags NVT ensemble 
 IBRION = 0                   # choose molecular dynamics 
 MDALGO = 3                   # using Langevin thermostat 
 ISIF = 2                     # compute stress tensor but do not change box volume/shape 
 TEBEG = 300                  # set temperature 
 NSW = 10000                  # number of time steps 
 POTIM = 1.0                  # time step in femto seconds 
 LANGEVIN_GAMMA = 10.0   10.0 # setting friction coefficient in inverse time units

Note this INCAR file only contains the parameters for the molecular dynamics part. The electronic minimization or the machine learning tags have to be added.

@@ Line 1: / Line 1: @@
- [[:Category: NVT ensemble]]
+[[:Category: NVT ensemble]]
- == NVT ensemble ==
+== NVT ensemble ==
- The NVT ensemble is a statistical ensemble that is used to study material properties under the conditions of a
+The NVT ensemble is a statistical ensemble that is used to study material properties under the conditions of a
- constant particle number N, constant volume V and a temperature fluctuating around an equilibrium value T.
+constant particle number N, constant volume V and a temperature fluctuating around an equilibrium value T.
- This page describes how to sample the NVT ensemble from a [[Molecular dynamics calculations|molecular dynamics]] run.
+This page describes how to sample the NVT ensemble from a [[Molecular dynamics calculations|molecular dynamics]] run.
- ''' Instructions for setting up a NVT ensemble '''
+''' Instructions for setting up a NVT ensemble '''
- There are 4 choices of thermostats which can either be stochastic or deterministic to simulate the NVT ensemble.
+There are 4 choices of thermostats which can either be stochastic or deterministic to simulate the NVT ensemble.
- The stochastic [[Andersen thermostat]] or [[Langevin thermostat]],
+The stochastic [[Andersen thermostat]] or [[Langevin thermostat]],
- the deterministic [[Nose-Hoover thermostat]] or [[MDALGO#MDALGO.3D13:_Multiple_Anderson_thermostats|Multiple Anderson thermostats]] can be used.
+the deterministic [[Nose-Hoover thermostat]] or [[MDALGO#MDALGO.3D13:_Multiple_Anderson_thermostats|Multiple Anderson thermostats]] can be used.
- See table for the corresponding {{TAG|MDALGO}} tags.
+See table for the corresponding {{TAG|MDALGO}} tags.
- {|class="wikitable" style="margin:aut
+{|class="wikitable" style="margin:aut
- ! NVT ensemble !! Andersen !! Langevin !! Nose-Hoover !! Multiple Andersen
+! NVT ensemble !! Andersen !! Langevin !! Nose-Hoover !! Multiple Andersen
- |-
+|-
- |MDALGO           ||     1        ||     3            ||      2      ||        13
+|MDALGO           ||     1        ||     3            ||      2      ||        13
- |- {{TAG|PSUBSYS}
+|- {{TAG|PSUBSYS}
- | additional tags || {{TAG|ANDERSEN_PROB}} ||   {{TAG|LANGEVIN_GAMMA}} ||  {{TAG|SMASS}}      || {{TAG|NSUBSYS}}, {{TAG|TSUBSYS}}, {{TAG|PSUBSYS}}
+| additional tags || {{TAG|ANDERSEN_PROB}} ||   {{TAG|LANGEVIN_GAMMA}} ||  {{TAG|SMASS}}      || {{TAG|NSUBSYS}}, {{TAG|TSUBSYS}}, {{TAG|PSUBSYS}}
- |}
+|}
- The user is not allowed to set ISF>2 because then the volume of the simulation would be allowed to change
+The user is not allowed to set ISF>2 because then the volume of the simulation would be allowed to change
- and the simulation would not result in a NVT ensemble. Setting ISIF=0 will speed up the calculation because the stress
+and the simulation would not result in a NVT ensemble.
- tensor is not computed. The tag {{TAG|IBRION}} has to be set to zero.
+Other flags related to molecular dynamics simulations can be found [[Molecular dynamics calculations|here]].
+''A simple example INCAR file for the [[Langevin thermostat]] could look like'''
+  #INCAR molecular dynamics tags NVT ensemble
+  IBRION = 0                   # choose molecular dynamics
+  MDALGO = 3                   # using Langevin thermostat
+  ISIF = 2                     # compute stress tensor but do not change box volume/shape
+  TEBEG = 300                  # set temperature
+  NSW = 10000                  # number of time steps
+  POTIM = 1.0                  # time step in femto seconds
+  LANGEVIN_GAMMA = 10.0   10.0 # setting friction coefficient in inverse time units
+Note this {{FILE|INCAR}} file only contains the parameters for the molecular dynamics part. The
+[[Electronic minimization|electronic minimization]] or the [[Machine-learned force fields|machine learning]] tags
+have to be added.