MDALGO: Difference between revisions

From VASP Wiki
VisualWikitext
No edit summary
 
(3 intermediate revisions by one other user not shown)
Line 1: Line 1:
{{TAGDEF|MDALGO|0 {{!}} 1 {{!}} 2 {{!}} 3 {{!}} 4 {{!}} 5 {{!}} 11 {{!}} 21 {{!}} 13  |0}}
{{TAGDEF|MDALGO|0 {{!}} 1 {{!}} 2 {{!}} 3 {{!}} 4 {{!}} 5 {{!}} 11 {{!}} 21 {{!}} 13  |0}}


Description: {{TAG|MDALGO}} specifies the molecular-dynamics-simulation protocol (in case {{TAGO|IBRION|0}}). {{NB|mind|All options except {{TAGO|MDALGO|0}} require a {{VASP}} executable which was compiled with the [[Precompiler options#-Dtbdyn|<code>-Dtbdyn</code>]] precompiler option enabled. This is usually the case because the option is present by default in all shipped [[makefile.include]] templates since {{VASP}} 5.4.4.}}
Description: Specifies the thermostat for [[MD calculations]] (in case {{TAGO|IBRION|0}}).
----
----
{{NB|mind|All options except {{TAGO|MDALGO|0}} require a {{VASP}} executable which was compiled with the [[Precompiler options#-Dtbdyn|<code>-Dtbdyn</code>]] precompiler option enabled. This is usually the case because the option is present by default in all shipped [[makefile.include]] templates since {{VASP}} 5.4.4.}}
{{TAG|MDALGO}}=1,2,3,4,5 can be applied in the context of
* standard [[molecular-dynamics calculations]]
* [[constrained molecular dynamics]]
* [[metadynamics calculations]]
* the [[slow-growth approach]]
* monitoring geometric parameters using the {{FILE|ICONST}} file
* [[Biased molecular dynamics]]


== {{TAG|MDALGO}}=0: Standard molecular dynamics ==
The main output file is the {{FILE|REPORT}} file.
Selects a [[Nosé-Hoover thermostat]] which allows to sample the [[NVT ensemble]]. {{NB|deprecated|If possible, we recommend to use one of the newer Nosé–Hoover thermostat implementations {{VASP}} provides ({{TAGO|MDALGO|2 or 4}}). While the results (ensemble averages) should be identical this variant comes with some drawbacks regarding post-processing: the atom coordinates in output files will always be wrapped back into the box if atoms cross the periodic boundaries. This makes it impossible to carry out certain analysis, e.g., computing the mean squared displacement (MSD).}}


=== Standard molecular dynamics ===
== {{TAGDEF|MDALGO|0}}: Standard molecular dynamics ==
Selects a [[Nosé-Hoover thermostat]] which allows sampling the [[NVT ensemble]] at temperature {{TAG|TEBEG}}. The [[Nosé-Hoover thermostat]] requires an appropriate setting for {{TAG|SMASS}}. To sample an [[NVE ensemble]] set {{TAGDEF|SMASS|-3}}. {{NB|deprecated|If possible, we recommend using one of the newer Nosé-Hoover thermostat implementations {{VASP}} provides ({{TAGO|MDALGO|2 or 4}}). While the results (ensemble averages) should be identical ,this variant comes with some drawbacks regarding post-processing: the atom coordinates in output files will always be wrapped back into the box if atoms cross the periodic boundaries. This makes it impossible to carry out certain analysis, e.g., computing the mean squared displacement (MSD).}}


* For a standard molecular dynamics run with Nose-Hoover thermostat, one has to:
== {{TAGDEF|MDALGO|1}}: [[Andersen thermostat]] ==
#Set the standard MD-related tags: {{TAG|IBRION}}=0, {{TAG|TEBEG}}, {{TAG|POTIM}}, and {{TAG|NSW}}.
#Set {{TAG|MDALGO}}=0, and choose an appropriate setting for {{TAG|SMASS}}.


=== Special case: NVE ensemble ===
The [[Andersen thermostat]] can be used to sample an [[NVT ensemble]], which requires setting an appropriate value for {{TAG|ANDERSEN_PROB}}.


*To perform a calculation in the [[NVE ensemble]]:
For {{TAGDEF|ANDERSEN_PROB|0.0}}, the thermostat is inactive, such that the [[NVE ensemble]] is sampled. This is usually done after thermalization to a certain target temperature. {{NB|tip|Leave the value for {{TAG|TEBEG}} that was set in the thermalization. For {{TAGDEF|TEBEG|<0.1}}, some part of the code assumes it is used for [[structure optimization]] and not an [[molecular-dynamics calculations|MD run]].}}  
#Set the standard MD-related tags: {{TAG|IBRION}}=0, {{TAG|TEBEG}}, {{TAG|POTIM}}, and {{TAG|NSW}}.
#Set {{TAG|MDALGO}}=0 and set {{TAG|SMASS}}=-3.


== {{TAG|MDALGO}}=1: Andersen thermostat ==
== {{TAGDEF|MDALGO|2}}: [[Nosé-Hoover thermostat]] ==


For the description of the Andersen thermostat see: {{TAG|Andersen thermostat}}.
The [[Nosé-Hoover thermostat]] is currently only available for the [[NVT ensemble]]. It requires setting an appropriate value for {{TAG|SMASS}}.
{{NB|tip|The [[Nosé-Hoover thermostat]] is a special case of the [[Nosé-Hoover chain thermostat]] ({{TAGDEF|MDALGO|4}} with {{TAGDEF|NHC_NCHAINS|1}}). The control tags for {{TAGDEF|MDALGO|4}} may be more convenient to use than the older implementation ({{TAGDEF|MDALGO|2}}).}}


The Andersen thermostat is only available for the NVT ensemble.
== {{TAGDEF|MDALGO|3}}: [[Langevin thermostat]] ==


=== Standard molecular dynamics in ===
The [[Langevin thermostat]] is available for sampling an [[NVT ensemble]] as well as sampling an [[NpT ensemble]].
* For an [[NVT ensemble]], fix the cell shape and volume with {{TAGDEF|ISIF|2}} and set an appropriate value for the friction coefficients for all species in the {{FILE|POSCAR}} file by means of the {{TAG|LANGEVIN_GAMMA}} tag.
* The Langevin dynamics in the [[NpT ensemble]] is calculated by the method of Parrinello and Rahman{{cite|parrinello:prl:1980}}{{cite|parrinello:jap:1981}} combined with a [[Langevin thermostat]].
#Set {{TAGDEF|ISIF|3}} to allow for relaxation of the cell volume and shape. At the moment, dynamics with ''fixed volume+variable shape'' ({{TAG|ISIF}}=4) or ''fixed shape+variable volume'' ({{TAG|ISIF}}=7) are not available.
#Specify friction coefficients for all species in the {{FILE|POSCAR}} file by means of the {{TAG|LANGEVIN_GAMMA}} tag.
#Specify a separate set of friction coefficient for the lattice degrees-of-freedom using the {{TAG|LANGEVIN_GAMMA_L}} tag.
#Set a mass for the lattice degrees-of-freedom using the {{TAG|PMASS}} tag.
#Optionally, one may define an external pressure (in kB) by means of the {{TAG|PSTRESS}} tag.


*For a standard molecular dynamics run with Anderson thermostat, one has to:
Also see [[stochastic boundary conditions]].
#Set the standard MD-related tags: {{TAG|IBRION}}=0, {{TAG|TEBEG}}, {{TAG|POTIM}}, and {{TAG|NSW}}.
#Set {{TAG|MDALGO}}=1, and choose an appropriate setting for {{TAG|ANDERSEN_PROB}}


=== Constrained molecular dynamics ===
== {{TAGDEF|MDALGO|4}}: [[Nosé-Hoover chain thermostat]] ==
For a description of constrained molecular dynamics see {{TAG|Constrained molecular dynamics}}.
The [[Nosé-Hoover chain thermostat]] can be used to sample an [[NVT ensemble]] and requires selecting the number of thermostats in the chain via {{TAG|NHC_NCHAINS}} as well as choosing an appropriate setting for the thermostat parameter {{TAG|NHC_PERIOD}}.  


* For a constrained molecular dynamics run with Andersen thermostat, one has to:
== {{TAG|MDALGO}}=5: [[CSVR thermostat|Canonical sampling through velocity-rescaling (CSVR thermostat)]] ==
#Set the standard MD-related tags: {{TAG|IBRION}}=0, {{TAG|TEBEG}}, {{TAG|POTIM}}, and {{TAG|NSW}}
#Set {{TAG|MDALGO}}=1, and choose an appropriate setting for {{TAG|ANDERSEN_PROB}}
#Define geometric constraints in the {{FILE|ICONST}}-file, and set the STATUS parameter for the constrained coordinates to 0
#When the free-energy gradient is to be computed, set {{TAG|LBLUEOUT}}=.TRUE.
 
=== Slow-growth approach ===
For a description of slow-growth approach see {{TAG|Slow-growth approach}}.
* For a slow-growth simulation, one has to perform a calcualtion very similar to {{TAG|Constrained molecular dynamics}} but additionally the transformation velocity-related {{TAG|INCREM}}-tag for each geometric parameter with <tt>STATUS=0</tt> has to be specified:
#Set the standard MD-related tags: {{TAG|IBRION}}=0, {{TAG|TEBEG}}, {{TAG|POTIM}}, and {{TAG|NSW}}
#Set {{TAG|MDALGO}}=1, and choose an appropriate setting for {{TAG|ANDERSEN_PROB}}
#Define geometric constraints in the {{FILE|ICONST}}-file, and set the STATUS parameter for the constrained coordinates to 0
#When the free-energy gradient is to be computed, set {{TAG|LBLUEOUT}}=.TRUE.
 
<ol start="5">
<li>Specify the transformation velocity-related {{TAG|INCREM}}-tag for each geometric parameter with <tt>STATUS=0</tt>.</li>
</ol>
 
=== Monitoring geometric parameters ===
Geometric parameters with <tt>STATUS = 7</tt> in the {{FILE|ICONST}}-file are monitored during the MD simulation.
The corresponding values are written onto the {{FILE|REPORT}}-file, for each MD step, after the lines following the string <tt>Monit_coord</tt>.
 
Sometimes it is desirable to terminate the simulation if all values of monitored parameters get larger that some predefined upper and/or lower limits. These limits can be set by the user by means of the {{TAG|VALUE_MAX}} and {{TAG|VALUE_MIN}}-tags.
 
*To monitor geometric parameters during an MD run:
#Set the standard MD-related tags: {{TAG|IBRION}}=0, {{TAG|TEBEG}}, {{TAG|POTIM}}, and {{TAG|NSW}}.
#Set {{TAG|MDALGO}}=1, and choose an appropriate setting for {{TAG|ANDERSEN_PROB}}.
#Define geometric constraints in the {{FILE|ICONST}}-file, and set the <tt>STATUS</tt> parameter for the constrained coordinates to 7.
#Optionally, set the upper and/or lower limits for the coordinates, by means of the {{TAG|VALUE_MAX}} and {{TAG|VALUE_MIN}}-tags, respectively.
 
=== Metadynamics ===
For a description of metadynamics see {{TAG|Metadynamics}}.
 
* For a metadynamics run with Andersen thermostat, one has to:
#Set the standard MD-related tags: {{TAG|IBRION}}=0, {{TAG|TEBEG}}, {{TAG|POTIM}}, and {{TAG|NSW}}.
#Set {{TAG|MDALGO}}=1 (or {{TAG|MDALGO}}=11 in VASP 5.x), and choose an appropriate setting for {{TAG|ANDERSEN_PROB}}.
#Set the parameters {{TAG|HILLS_H}}, {{TAG|HILLS_W}}, and {{TAG|HILLS_BIN}}.
#Define collective variables in the {{FILE|ICONST}}-file, and set the {{TAG|STATUS}} parameter for the collective variables to 5.
#If needed, define the bias potential in the {{FILE|PENALTYPOT}}-file.
 
The actual time-dependent bias potential is written to the {{FILE|HILLSPOT}}-file, which is updated after adding a new Gaussian. At the beginning of the simulation, VASP attempts to read the initial bias potential from the {{FILE|PENALTYPOT}}-file. For the continuation of a metadynamics run, copy {{FILE|HILLSPOT}} to {{FILE|PENALTYPOT}}. The values of all collective variables for each MD step are listed in {{FILE|REPORT}}-file, check the lines after the string <tt>Metadynamics</tt>.
 
=== Biased molecular dynamics ===
 
For a description of biased molecular dynamics see {{TAG|Biased molecular dynamics}}.
 
* For a biased molecular dynamics run with Andersen thermostat, one has to:
#Set the standard MD-related tags: {{TAG|IBRION}}=0, {{TAG|TEBEG}}, {{TAG|POTIM}}, and {{TAG|NSW}}.
#Set {{TAG|MDALGO}}=1 (or {{TAG|MDALGO}}=11 in VASP 5.x), and choose an appropriate setting for {{TAG|ANDERSEN_PROB}}.
#In order to avoid updating of the bias potential, set {{TAG|HILLS_BIN}}={{TAG|NSW}}.
#Define collective variables in the {{FILE|ICONST}}-file, and set the <tt>STATUS</tt> parameter for the collective variables to 5.
#Define the bias potential in the {{FILE|PENALTYPOT}}-file if necessary.
 
The values of all collective variables for each MD step are listed in the {{FILE|REPORT}}-file, check the lines after the string <tt>Metadynamics</tt>.
 
=== Special case: NVE ensemble ===
 
NVE ensemble calculations can be also run by selecting the Anderson thermostat and setting the update collision probability ({{TAG|ANDERSEN_PROB}}) to zero.
*To run an NVE ensemble:
#Set the standard MD-related tags: {{TAG|IBRION}}=0, {{TAG|TEBEG}}, {{TAG|POTIM}}, and {{TAG|NSW}}.
#Set {{TAG|MDALGO}}=1 and {{TAG|ANDERSEN_PROB}}=0.0.
 
== {{TAG|MDALGO}}=2: Nose-Hoover thermostat==
 
For the description of the Nose-Hoover thermostat see: {{TAG|Nose-Hoover thermostat}}.
 
The Nose-Hoover thermostat is currently only available for the NVT ensemble.
 
=== Standard molecular dynamics ===
 
* For a standard molecular dynamics run with Nose-Hoover thermostat, one has to:
#Set the standard MD-related tags: {{TAG|IBRION}}=0, {{TAG|TEBEG}}, {{TAG|POTIM}}, and {{TAG|NSW}}.
#Set {{TAG|MDALGO}}=2, and choose an appropriate setting for {{TAG|SMASS}}.
 
=== Constrained molecular dynamics ===
For a description of constrained molecular dynamics see {{TAG|Constrained molecular dynamics}}.
 
* For a constrained molecular dynamics run with Nose-Hoover thermostat, one has to:
#Set the standard MD-related tags: {{TAG|IBRION}}=0, {{TAG|TEBEG}}, {{TAG|POTIM}}, and {{TAG|NSW}}.
#Set {{TAG|MDALGO}}=2, and choose an appropriate setting for {{TAG|SMASS}}.
#Define geometric constraints in the {{FILE|ICONST}}-file, and set the STATUS parameter for the constrained coordinates to 0.
#When the free-energy gradient is to be computed, set {{TAG|LBLUEOUT}}=.TRUE.
 
=== Slow-growth approach ===
 
For a description of slow-growth approach see {{TAG|Slow-growth approach}}.
 
* For a slow-growth approach run with Nose-Hoover thermostat, one has to:
#Set the standard MD-related tags: {{TAG|IBRION}}=0, {{TAG|TEBEG}}, {{TAG|POTIM}}, and {{TAG|NSW}}
#Set {{TAG|MDALGO}}=2, and choose an appropriate setting for {{TAG|SMASS}}
#Define geometric constraints in the {{FILE|ICONST}}-file, and set the <tt>STATUS</tt> parameter for the constrained coordinates to 0
#When the free-energy gradient is to be computed, set {{TAG|LBLUEOUT}}=.TRUE.
<ol start="5">
<li>Specify the transformation velocity-related {{TAG|INCREM}}-tag for each geometric parameter with <tt>STATUS=0</tt></li>
</ol>
 
=== Monitoring geometric parameters ===
Geometric parameters with <tt>STATUS = 7</tt> in the {{FILE|ICONST}}-file are monitored during the MD simulation.
The corresponding values are written onto the {{FILE|REPORT}}-file, for each MD step, after the lines following the string <tt>Monit_coord</tt>.
 
Sometimes it is desirable to terminate the simulation if all values of monitored parameters get larger that some predefined upper and/or lower limits. These limits can be set by the user by means of the {{TAG|VALUE_MAX}} and {{TAG|VALUE_MIN}}-tags.
 
To monitor geometric parameters during an MD run:
#Set the standard MD-related tags: {{TAG|IBRION}}=0, {{TAG|TEBEG}}, {{TAG|POTIM}}, and {{TAG|NSW}}
#Set {{TAG|MDALGO}}=2, and choose an appropriate setting for {{TAG|SMASS}}
#Define geometric constraints in the {{FILE|ICONST}}-file, and set the <tt>STATUS</tt> parameter for the constrained coordinates to 7
#Optionally, set the upper and/or lower limits for the coordinates, by means of the {{TAG|VALUE_MAX}} and {{TAG|VALUE_MIN}}-tags, respectively.
<div id="Langevin"></div>
 
=== Metadynamics ===
For a description of metadynamics see {{TAG|Metadynamics}}.
 
* For a metadynamics run with Nose-Hoover thermostat, one has to:
#Set the standard MD-related tags: {{TAG|IBRION}}=0, {{TAG|TEBEG}}, {{TAG|POTIM}}, and {{TAG|NSW}}.
#Set {{TAG|MDALGO}}=2 (or {{TAG|MDALGO}}=21 in VASP 5.x), and choose an appropriate setting for {{TAG|SMASS}}.
#Set the parameters {{TAG|HILLS_H}}, {{TAG|HILLS_W}}, and {{TAG|HILLS_BIN}}.
#Define collective variables in the {{FILE|ICONST}}-file, and set the <tt>STATUS</tt> parameter for the collective variables to 5.
#If needed, define the bias potential in the {{FILE|PENALTYPOT}}-file.
 
The actual time-dependent bias potential is written to the {{FILE|HILLSPOT}}-file, which is updated after adding a new Gaussian. At the beginning of the simulation, VASP attempts to read the initial bias potential from the {{FILE|PENALTYPOT}}-file. For the continuation of a metadynamics run, copy {{FILE|HILLSPOT}} to {{FILE|PENALTYPOT}}. The values of all collective variables for each MD step are listed in {{FILE|REPORT}}-file, check the lines after the string <tt>Metadynamics</tt>.
 
=== Biased molecular dynamics ===
 
For a description of biased molecular dynamics see {{TAG|Biased molecular dynamics}}.
 
* For a biased molecular dynamics run with Nose-Hoover thermostat, one has to:
#Set the standard MD-related tags: {{TAG|IBRION}}=0, {{TAG|TEBEG}}, {{TAG|POTIM}}, and {{TAG|NSW}}.
#Set {{TAG|MDALGO}}=2 (or {{TAG|MDALGO}}=21 in VASP 5.x), and choose an appropriate setting for {{TAG|SMASS}}.
#In order to avoid updating of the bias potential, set {{TAG|HILLS_BIN}}={{TAG|NSW}}.
#Define collective variables in the {{FILE|ICONST}}-file, and set the <tt>STATUS</tt> parameter for the collective variables to 5.
#Define the bias potential in the {{FILE|PENALTYPOT}}-file.
 
The values of all collective variables for each MD step are listed in the {{FILE|REPORT}}-file, check the lines after the string <tt>Metadynamics</tt>.
 
== {{TAG|MDALGO}}=3: Langevin thermostat ==
 
For the description of the Langevin thermostat see: {{TAG|Langevin thermostat}}.
 
The Langevin thermostat is available for [[NVT ensemble|NVT]] and [[NpT ensemble|NpT]] ensembles.
 
=== ''NVT''-simulation with Langevin thermostat ===
 
*To run an ''NVT''-simulation with a Langevin thermostat, one has to:
{{Template:NVT_Langevin_thermostat_recipe}}
 
=== ''NpT''-simulation with Langevin thermostat ===
The Langevin dynamics in the isobaric-isothermal ensemble is calculated by the method of Parrinello and Rahman{{cite|parrinello:prl:1980}}{{cite|parrinello:jap:1981}} (see {{TAG|NpT ensemble}} for more descriptions) combined with a {{TAG|Langevin thermostat}}.
 
*To run an NpT-simulation (Parinello-Rahman dynamics) with a Langevin thermostat, one has to:
{{Template:NpT_Langevin_thermostat_recipe}}
 
Note that the advanced simulation protocols, such metadynamics, or geometric constraints and restraints are available within the {{TAG|MDALGO}}=3 setting.
 
=== Stochastic boundary conditions ===
In some cases it is desirable to study approach of initially non-equilibrium system to equilibrium. Examples of such simulations include the impact problems when a particle with large kinetic energy hits a surface or calculation of friction force between two surfaces sliding with respect to each other. As shown by Toton ''et al.''<ref name="Toton10"/>, this type of problems can be studied using the stochastic boundary conditions (SBC) derived from the generalized Langevin equation by Kantorovich and Rompotis.<ref name="Kantorovich08"/> In this approach, the system of interest is divided into three regions: (a) fixed atoms, (b) the internal (Newtonian) atoms moving according to Newtonian dynamics, and (c) a buffer region of Langevin atoms (''i.e.'', atoms governed by [[#LangevinEOM|Langevin equations of motion]]) located between (a) and (b).
 
The role of the Langevin atoms is to dissipate heat, while the fixed atoms are needed solely to create the correct potential well for the Langevin atoms to move in. The Newtonian region should include all atoms relevant to the process under study: in the case of the impact problem, for instance, the Newtonian region should contain atoms of the molecule hitting the surface and several uppermost layers of the material forming the surface. Performing molecular dynamics with such a setup guarantees that the system (possibly out of equilibrium initially) arrives at the appropriate canonical distribution.
 
*To run a simulation with stochastic boundary conditions, one has to:
#Set the standard MD-related tags: {{TAG|IBRION}}=0, {{TAG|TEBEG}}, {{TAG|POTIM}}, and {{TAG|NSW}}
#Set {{TAG|ISIF}}=2
#Set {{TAG|MDALGO}}=3 to invoke the Langevin thermostat
#Prepare the {{FILE|POSCAR}} file in such a way that the Newtonian and Langevin atoms are treated as different species (even if they are chemically identical). In your {{FILE|POSCAR}}, use "selective dynamics" and the corresponding logical flags to define the frozen and moveable atoms.
#Specify friction coefficients &gamma;, for all species in the {{FILE|POSCAR}} file, by means of the {{TAG|LANGEVIN_GAMMA}}-tag: set the friction coefficients to 0 for all fixed and Newtonian atoms, and choose a proper &gamma; for the Langevin atoms.
 
==== Practical example ====
 
Consider a system consisting of 16 hydrogen and 48 silicon atoms. Suppose that eight silicon atoms are considered to be Langevin atoms and the remaining 32 Si atoms are either fixed or Newtonian atoms. The Langevin and Newtonian (or fixed) atoms should be considered as different species, ''i.e.'', the {{FILE|POSCAR}}-file should contain the line like this:
 
Si H Si
40 16 8
 
As only the final eight Si atoms are considered to be Langevin atoms, the {{FILE|INCAR}}-file should contain the following line defining the friction coefficients:
 
LANGEVIN_GAMMA = 0.0  0.0  10.0
 
''i.e.'', for all non-Langevin atoms, &gamma; should be set to zero.
 
<div id="multiAnderson"></div>
 
== {{TAG|MDALGO}}=4: Nose Hoover chains thermostat ==
For a description of Nose Hoover chains theromastat, see {{TAG|NHC thermostat}}.
*To run an ''NVT''-simulation with a Nose Hoover chains thermostat, one has to:
#Set the standard MD-related tags: {{TAG|IBRION}}=0, {{TAG|TEBEG}}, {{TAG|POTIM}}, and {{TAG|NSW}}
#Set {{TAG|ISIF}}=2
#Set {{TAG|MDALGO}}=4 to invoke the thermostat, select the number of thermostats in chain via {{TAG|NHC_NCHAINS}} and choose an appropriate setting for the thermostat parameter {{TAG|NHC_PERIOD}}.
Note that the advanced simulation protocols, such metadynamics, or geometric constraints and restraints are available within the {{TAG|MDALGO}}=4 setting.
 
== {{TAG|MDALGO}}=5: Canonical sampling through velocity-rescaling (CSVR thermostat) ==
{{NB|mind|This option is available as of VASP 6.4.3.}}
{{NB|mind|This option is available as of VASP 6.4.3.}}


For a description of CSVR theromastat, see {{TAG|CSVR thermostat}}.
The [[CSVR thermostat]] can be used to sample an [[NVT ensemble]]. It requires setting {{TAG|CSVR_PERIOD}}.
*To run an ''NVT''-simulation with a CSVR thermostat, one has to:
#Set the standard MD-related tags: {{TAG|IBRION}}=0, {{TAG|TEBEG}}, {{TAG|POTIM}}, and {{TAG|NSW}}
#Set {{TAG|ISIF}}=2
#Set {{TAG|MDALGO}}=5 to invoke the thermostat and choose an appropriate setting for the thermostat parameter {{TAG|CSVR_PERIOD}}.
Note that the advanced simulation protocols, such metadynamics, or geometric constraints and restraints are available within the {{TAG|MDALGO}}=5 setting.


== {{TAG|MDALGO}}=13: Multiple Andersen thermostats ==
== {{TAG|MDALGO}}=13: Multiple Andersen thermostats ==
Up to three user-defined atomic subsystems may be coupled with independent Andersen thermostats<ref name="Andersen80"/> (see remarks under {{TAG|MDALGO}}=1 as well).
<div id="multiAnderson"></div>
Up to three user-defined atomic subsystems may be coupled with independent [[Andersen thermostat|Andersen thermostats]]<ref name="Andersen80"/> (see remarks under {{TAG|MDALGO}}=1 as well).
The {{FILE|POSCAR}} file must be organized such that the positions of atoms of subsystem ''i+1'' are defined after those for the subsystem ''i'', and the following flags must be set by the user:  
The {{FILE|POSCAR}} file must be organized such that the positions of atoms of subsystem ''i+1'' are defined after those for the subsystem ''i'', and the following flags must be set by the user:  
*{{TAG|NSUBSYS}}=[int array]  
*{{TAG|NSUBSYS}}=[int array]  
Line 247: Line 65:
:Collision probability for atoms in each subsystem. Only the values 0&le;{{TAG|PSUBSYS}}&le;1 are allowed.
:Collision probability for atoms in each subsystem. Only the values 0&le;{{TAG|PSUBSYS}}&le;1 are allowed.


== Related Tags and Sections ==
== Related tags and articles ==
{{TAG|IBRION}},
{{TAG|IBRION}},
{{TAG|ISIF}},
{{TAG|ISIF}},
Line 278: Line 96:
<references>
<references>
<ref name="Andersen80">[http://dx.doi.org/10.1063/1.439486 H. C. Andersen, J. Chem. Phys. 72, 2384 (1980).]</ref>
<ref name="Andersen80">[http://dx.doi.org/10.1063/1.439486 H. C. Andersen, J. Chem. Phys. 72, 2384 (1980).]</ref>
<ref name="Toton10">[http://dx.doi.org/10.1088/0953-8984/22/7/074205 D. Toton, C. D. Lorenz, N. Rompotis, N. Martsinovich, and L. Kantorovich, J. Phys.: Condens. Matter 22, 074205 (2010).]</ref>
<ref name="Kantorovich08">[http://dx.doi.org/10.1103/PhysRevB.78.094305 L. Kantorovich and N. Rompotis, Phys. Rev. B 78, 094305 (2008).]</ref>
</references>
</references>
----
----


[[Category:INCAR tag]][[Category:Molecular dynamics]][[Category:Advanced molecular-dynamics sampling]][[Category:Howto]]
[[Category:INCAR tag]][[Category:Molecular dynamics]][[Category:Advanced molecular-dynamics sampling]][[Category:Howto]]

Latest revision as of 13:46, 21 October 2025

MDALGO = 0 | 1 | 2 | 3 | 4 | 5 | 11 | 21 | 13
Default: MDALGO = 0 

Description: Specifies the thermostat for MD calculations (in case IBRION = 0).

Mind: All options except MDALGO = 0 require a VASP executable which was compiled with the -Dtbdyn precompiler option enabled. This is usually the case because the option is present by default in all shipped makefile.include templates since VASP 5.4.4.

MDALGO=1,2,3,4,5 can be applied in the context of

The main output file is the REPORT file.

MDALGO = 0 : Standard molecular dynamics

Selects a Nosé-Hoover thermostat which allows sampling the NVT ensemble at temperature TEBEG. The Nosé-Hoover thermostat requires an appropriate setting for SMASS. To sample an NVE ensemble set SMASS = -3 .

Deprecated: If possible, we recommend using one of the newer Nosé-Hoover thermostat implementations VASP provides (MDALGO = 2 or 4). While the results (ensemble averages) should be identical ,this variant comes with some drawbacks regarding post-processing: the atom coordinates in output files will always be wrapped back into the box if atoms cross the periodic boundaries. This makes it impossible to carry out certain analysis, e.g., computing the mean squared displacement (MSD).

MDALGO = 1 : Andersen thermostat

The Andersen thermostat can be used to sample an NVT ensemble, which requires setting an appropriate value for ANDERSEN_PROB.

For ANDERSEN_PROB = 0.0 , the thermostat is inactive, such that the NVE ensemble is sampled. This is usually done after thermalization to a certain target temperature.

Tip: Leave the value for TEBEG that was set in the thermalization. For TEBEG = <0.1 , some part of the code assumes it is used for structure optimization and not an MD run.

MDALGO = 2 : Nosé-Hoover thermostat

The Nosé-Hoover thermostat is currently only available for the NVT ensemble. It requires setting an appropriate value for SMASS.

Tip: The Nosé-Hoover thermostat is a special case of the Nosé-Hoover chain thermostat (MDALGO = 4  with NHC_NCHAINS = 1 ). The control tags for MDALGO = 4  may be more convenient to use than the older implementation (MDALGO = 2 ).

MDALGO = 3 : Langevin thermostat

The Langevin thermostat is available for sampling an NVT ensemble as well as sampling an NpT ensemble.

  1. Set ISIF = 3  to allow for relaxation of the cell volume and shape. At the moment, dynamics with fixed volume+variable shape (ISIF=4) or fixed shape+variable volume (ISIF=7) are not available.
  2. Specify friction coefficients for all species in the POSCAR file by means of the LANGEVIN_GAMMA tag.
  3. Specify a separate set of friction coefficient for the lattice degrees-of-freedom using the LANGEVIN_GAMMA_L tag.
  4. Set a mass for the lattice degrees-of-freedom using the PMASS tag.
  5. Optionally, one may define an external pressure (in kB) by means of the PSTRESS tag.

Also see stochastic boundary conditions.

MDALGO = 4 : Nosé-Hoover chain thermostat

The Nosé-Hoover chain thermostat can be used to sample an NVT ensemble and requires selecting the number of thermostats in the chain via NHC_NCHAINS as well as choosing an appropriate setting for the thermostat parameter NHC_PERIOD.

MDALGO=5: Canonical sampling through velocity-rescaling (CSVR thermostat)

Mind: This option is available as of VASP 6.4.3.

The CSVR thermostat can be used to sample an NVT ensemble. It requires setting CSVR_PERIOD.

MDALGO=13: Multiple Andersen thermostats

Up to three user-defined atomic subsystems may be coupled with independent Andersen thermostats[3] (see remarks under MDALGO=1 as well). The POSCAR file must be organized such that the positions of atoms of subsystem i+1 are defined after those for the subsystem i, and the following flags must be set by the user:

Define the last atom for each subsystem (two or three values must be supplied). For instance, if total of 20 atoms is defined in the POSCAR file, and the initial 10 atoms belong to the subsystem 1, the next 7 atoms to the subsystem 2, and the last 3 atoms to the subsystem 3, NSUBSYS should be defined as follows:
NSUBSYS= 10 17 20
Note that the last number in the previous example is actually redundant (clearly the last three atoms belong to the last subsystem) and does not have to be user-supplied.
Simulation temperature for each subsystem
Collision probability for atoms in each subsystem. Only the values 0≤PSUBSYS≤1 are allowed.

Related tags and articles

IBRION, ISIF, SMASS, ANDERSEN_PROB, RANDOM_SEED, LBLUEOUT, SHAKETOL, SHAKEMAXITER, HILLS_H, HILLS_W, HILLS_BIN, INCREM, VALUE_MIN, VALUE_MAX, LANGEVIN_GAMMA, LANGEVIN_GAMMA_L, PMASS, NSUBSYS, TSUBSYS, PSUBSYS, ICONST, PENALTYPOT, HILLSPOT, REPORT

Examples that use this tag

References

  1. ↑ M. Parrinello and A. Rahman, Phys. Rev. Lett. 45, 1196 (1980).
  2. ↑ M. Parrinello and A. Rahman, J. Appl. Phys. 52, 7182 (1981).
  3. ↑ H. C. Andersen, J. Chem. Phys. 72, 2384 (1980).
Retrieved from "https://vasp.at/wiki/index.php?title=MDALGO&oldid=32646"