Langevin thermostat: Difference between revisions

From VASP Wiki
No edit summary
No edit summary
 
(4 intermediate revisions by one other user not shown)
Line 10: Line 10:
with Δ''t'' being the time-step used in the MD to integrate the equations of motion. Obviously, Langevin dynamics is identical to the classical Hamiltonian in the limit of vanishing γ.
with Δ''t'' being the time-step used in the MD to integrate the equations of motion. Obviously, Langevin dynamics is identical to the classical Hamiltonian in the limit of vanishing γ.


The friction coefficient is set by the {{TAG|LANGEVIN_GAMMA}} parameter.
*{{TAG|NVT ensemble}}:
The friction coefficient is set by the {{TAG|LANGEVIN_GAMMA}} parameter.
*{{TAG|NpT ensemble}}:
As for the {{TAG|NVT ensemble}} the {{TAG|LANGEVIN_GAMMA}} parameter has to be set. If the {{TAG|NpT ensemble}} is used (by setting {{TAG|ISIF}}=3) additionally the friction coefficient of the lattice {{TAG|LANGEVIN_GAMMA_L}} has to be provided too.


The Langevin thermostat is selected by MDALGO=3.  
The Langevin thermostat is selected by MDALGO=3.  
Line 19: Line 22:
<noinclude>
<noinclude>
----
----
[[Category:Molecular Dynamics]][[Category:Thermostats]][[Category:Theory]][[Category:Howto]]
[[Category:Molecular dynamics]][[Category:Thermostats]][[Category:Theory]][[Category:Howto]]

Latest revision as of 14:04, 7 April 2022

The Langevin thermostat[1][2][3] maintains the temperature through a modification of Newton's equations of motion

where Fi is the force acting on atom i due to the interaction potential, γi is a friction coefficient, and fi is a random force simulating the random kicks by the damping of particles between each other due to friction. The random numbers are chosen from a Gaussian distribution with the following variance

with Δt being the time-step used in the MD to integrate the equations of motion. Obviously, Langevin dynamics is identical to the classical Hamiltonian in the limit of vanishing γ.

The friction coefficient is set by the LANGEVIN_GAMMA parameter.

As for the NVT ensemble the LANGEVIN_GAMMA parameter has to be set. If the NpT ensemble is used (by setting ISIF=3) additionally the friction coefficient of the lattice LANGEVIN_GAMMA_L has to be provided too.

The Langevin thermostat is selected by MDALGO=3.

References