Nose-Hoover thermostat: Difference between revisions

Revision as of 08:45, 31 May 2019

In the approach by Nosé and Hoover^[1]^[2]^[3] an extra degree of freedom is introduced in the Hamiltonian. The heat bath is considered as an integral part of the system and has a fictious coordinate $s$ which is introduced into the Lagrangian of the system. This Lagrangian for an $N$ is written as

${\mathcal {L}}=\sum \limits _{i=1}^{N}{\frac {m_{i}}{2}}s^{2}{\dot {\mathbf {r}}}_{i}^{2}-U({\mathbf {r}})+{\frac {Q}{2}}{\dot {s}}^{2}-gk_{B}T\mathrm {ln} \,s$

where $m_{i}$ and $k_{B}$ are the mass of ion $i$ and the Boltzmann constant, respectively. The first two terms are the kinetic and potential energy of the system. The third and fourth term represent the kinetic and potential energy of the fictitious coordinate $s$ . The parameter $g$ is usually equal to the number of degrees of freedom of the system $g=3N-N_{\mathrm {constraint} }$ , where $N_{\mathrm {constraint} }$ is equal to the number of constraint set (fixed coordinates in the POSCAR file). The parameter $Q$ is an effective "mass" of $s$ , which controls the coupling of the system to the heat bath. It is set by the INCAR tag SMASS.

References

[nose:jcp:1984-1] S. Nosé, J. Chem. Phys. 81, 511 (1984).

[nose:ptp:1991-2] S. Nosé, Prog. Theor. Phys. Suppl. 103, 1 (1991).

[hoover:pra:1985-3] W. G. Hoover, Phys. Rev. A 31, 1695 (1985).

[1]

[2]

[3]

@@ Line 2: / Line 2: @@
 <math>
-\mathcal{L} = \sum\limits_{i=1}^{N} \frac{m_{i}}{2} s^{2} \dot{\bold{r}}_{i}^{2} - U(\bold{r}) + \frac{Q}{2} \dot{s}^{2}-\frac{g}{\beta}\mathrm{ln} \, s.
+\mathcal{L} = \sum\limits_{i=1}^{N} \frac{m_{i}}{2} s^{2} \dot{\bold{r}}_{i}^{2} - U(\bold{r}) + \frac{Q}{2} \dot{s}^{2}-g k_{B} T \mathrm{ln} \, s
 </math>
+where <math>m_{i}</math> and <math>k_{B}</math> are the mass of ion <math>i</math> and the Boltzmann constant, respectively. The first two terms are the kinetic and potential energy of the system. The third and fourth term represent the kinetic and potential energy of the fictitious coordinate <math>s</math>. The parameter <math>g</math> is usually equal to the number of degrees of freedom of the system <math>g=3N - N_{\mathrm{constraint}}</math>, where <math>N_{\mathrm{constraint}}</math> is equal to the number of constraint set (fixed coordinates in the {{TAG|POSCAR}} file). The parameter <math>Q</math> is an effective "mass" of <math>s</math>, which controls the coupling of the system to the heat bath. It is set by the {{TAG|INCAR}} tag {{TAG|SMASS}}.