NVE ensemble: Difference between revisions

From VASP Wiki
No edit summary
No edit summary
Line 1: Line 1:
The NVE 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 an internal energy E fluctuating around an equilibrium value E. This page describes how to sample the NVE ensemble from a [[Molecular dynamics calculations|molecular dynamics]] run.
The NVE 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 an internal energy E fluctuating around an equilibrium value E. This page describes how to sample the NVE ensemble from a [[Molecular dynamics calculations|molecular-dynamics]] run.
   
   
''' Instructions for setting up a NVE ensemble '''
''' Instructions for setting up a NVE ensemble '''


For the NVE ensemble there are two equivalent choices of thermostats to set up a [[Molecular dynamics calculations|molecular dynamics]] run.
For the NVE ensemble there are two equivalent choices of thermostats to set up a [[Molecular dynamics calculations|molecular-dynamics]] run.
The user can either set the [[Andersen thermostat]] or the [[Nose-Hoover thermostat]] to give a NVE ensemble. The [[Andersen thermostat]] thermostat can be used by setting the collision probability ({{TAG|ANDERSEN_PROB}}) with the fictitious heat bath to zero and [[Nose-Hoover thermostat]] by setting the mass of the virtual degree to minus three. Both settings will switch the thermostat of such that the velocities
The user can either set the [[Andersen thermostat]] or the [[Nose-Hoover thermostat]] to give a NVE ensemble. The [[Andersen thermostat]] thermostat can be used by setting the collision probability ({{TAG|ANDERSEN_PROB}}) with the fictitious heat bath to zero and [[Nose-Hoover thermostat]] by setting the mass of the virtual degree to minus three. Both settings will switch the thermostat of such that the velocities
are determined by the [[Hellmann-Feynman forces]] or [[Machine-learned force fields]] only.
are determined by the [[Hellmann-Feynman forces]] or [[Machine-learned force fields]] only.
Line 16: Line 16:
The additional tags in the column for every thermostat have to be set to the given values. Otherwise the NVE ensemble will not be realized.
The additional tags in the column for every thermostat have to be set to the given values. Otherwise the NVE ensemble will not be realized.
There are two implementations of the [[Nose-Hoover thermostat]] in VASP which will give the same results. The {{TAG|MDALGO}}=0 version can be used even if the code was compiled without the precompiler option [[Precompiler_options#-Dtbdyn|-Dtbdyn]].
There are two implementations of the [[Nose-Hoover thermostat]] in VASP which will give the same results. The {{TAG|MDALGO}}=0 version can be used even if the code was compiled without the precompiler option [[Precompiler_options#-Dtbdyn|-Dtbdyn]].
Other tags related to molecular dynamics simulations can be found [[Molecular dynamics calculations|here]].
Other tags related to molecular-dynamics simulations can be found [[Molecular dynamics calculations|here]].


''An example INCAR file for the [[Andersen thermostat]] could look like''  
''An example INCAR file for the [[Andersen thermostat]] could look like''  


   #INCAR molecular dynamics tags NVE ensemble  
   #INCAR molecular-dynamics tags NVE ensemble  
   IBRION = 0                  # choose molecular dynamics  
   IBRION = 0                  # choose molecular-dynamics  
   MDALGO = 1                  # using Andersen thermostat
   MDALGO = 1                  # using Andersen thermostat
   ISIF = 2                    # compute stress tensor but do not change box volume/shape  
   ISIF = 2                    # compute stress tensor but do not change box volume/shape  
Line 29: Line 29:
   ANDERSEN_PROB = 0.0          # setting friction coefficient in inverse time units for two atom types
   ANDERSEN_PROB = 0.0          # setting friction coefficient in inverse time units for two atom types


{{NB|mind| 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.}}
{{NB|mind| 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.}}





Revision as of 13:38, 20 July 2022

The NVE 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 an internal energy E fluctuating around an equilibrium value E. This page describes how to sample the NVE ensemble from a molecular-dynamics run.

Instructions for setting up a NVE ensemble

For the NVE ensemble there are two equivalent choices of thermostats to set up a molecular-dynamics run. The user can either set the Andersen thermostat or the Nose-Hoover thermostat to give a NVE ensemble. The Andersen thermostat thermostat can be used by setting the collision probability (ANDERSEN_PROB) with the fictitious heat bath to zero and Nose-Hoover thermostat by setting the mass of the virtual degree to minus three. Both settings will switch the thermostat of such that the velocities are determined by the Hellmann-Feynman forces or Machine-learned force fields only.

NVE ensemble Andersen Nose-Hoover
MDALGO 1 0 or 2
additional tags to set ANDERSEN_PROB=0.0 SMASS=-3

The additional tags in the column for every thermostat have to be set to the given values. Otherwise the NVE ensemble will not be realized. There are two implementations of the Nose-Hoover thermostat in VASP which will give the same results. The MDALGO=0 version can be used even if the code was compiled without the precompiler option -Dtbdyn. Other tags related to molecular-dynamics simulations can be found here.

An example INCAR file for the Andersen thermostat could look like

 #INCAR molecular-dynamics tags NVE ensemble 
 IBRION = 0                   # choose molecular-dynamics 
 MDALGO = 1                   # using Andersen 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 
 ANDERSEN_PROB = 0.0          # setting friction coefficient in inverse time units for two atom types


Mind: This INCAR file only contains the parameters for the molecular-dynamics part. The electronic minimization or the machine learning tags have to be added.


Related tags and articles

REPORT, Molecular dynamics calculations