Andersen thermostat: Difference between revisions
No edit summary |
No edit summary |
||
| (8 intermediate revisions by one other user not shown) | |||
| Line 1: | Line 1: | ||
In the approach proposed by Andersen | In the approach proposed by Andersen {{cite|andersen:jcp:1980}} the system is thermally coupled to a fictitious heat bath with the desired temperature. The coupling is represented by stochastic collisions that act occasionally on randomly selected particles. In particular the momentum of the ''lucky'' particle at every collision step is instantaneously chosen at random from the Boltzmann distribution at the selected temperature. The collision probability is defined as an average number of collisions per atom and time-step and the collision frequency occurs with the following distribution | ||
<math> | ::<math> | ||
P(t)=\nu e^{-\nu t}. | P(t)=\nu e^{-\nu t}. | ||
</math> | </math> | ||
The exponent of this the distribution (<math>\nu t</math>) is controlled by the flag {{TAG|ANDERSEN_PROB}}. Since <math>t</math> is the time step in the calculation {{TAG|ANDERSEN_PROB}} has to be scaled if the time step changes. The total number of collisions with the heat-bath is written out to the file {{FILE|REPORT}} for each MD step. | |||
A very good implementation of the Andersen thermostat can be found in chapter 6.1.1 of reference {{cite|frenkel:book:1996}}. | |||
[[Category:Molecular | The Andersen thermostat is selected by setting {{TAG|MDALGO}}=1. | ||
== References == | |||
<references/> | |||
<noinclude> | |||
---- | |||
[[Category:Molecular dynamics]][[Category:Thermostats]][[Category:Theory]][[Category:Howto]] | |||
Latest revision as of 13:50, 7 April 2022
In the approach proposed by Andersen [1] the system is thermally coupled to a fictitious heat bath with the desired temperature. The coupling is represented by stochastic collisions that act occasionally on randomly selected particles. In particular the momentum of the lucky particle at every collision step is instantaneously chosen at random from the Boltzmann distribution at the selected temperature. The collision probability is defined as an average number of collisions per atom and time-step and the collision frequency occurs with the following distribution
The exponent of this the distribution () is controlled by the flag ANDERSEN_PROB. Since is the time step in the calculation ANDERSEN_PROB has to be scaled if the time step changes. The total number of collisions with the heat-bath is written out to the file REPORT for each MD step.
A very good implementation of the Andersen thermostat can be found in chapter 6.1.1 of reference [2].
The Andersen thermostat is selected by setting MDALGO=1.
References