Abstract
A systematic algorithm to design multiple thermostat systems in the framework of the Nosé–Hoover type non-Hamiltonian formulation is presented. Using ‘non uniform’ time transformations in a generalised Hamiltonian equation, we develop the non-Hamiltonian equations of motion for multiple thermostat systems having an arbitrary number of thermostats and arbitrary connections between a physical system and thermostats (‘Nosé–Hoover network’). We then present the algorithm to construct the Nosé–Hoover network equations based on a simple diagram only. On the basis of this algorithm, recursively attached Nosé–Hoover thermostats are introduced as an example of the Nosé–Hoover network and its high efficiency in sampling the canonical distribution for an one-dimensional double-well system is illustrated by numerical calculations.
Acknowledgements
The author thanks Peter Daivis, Ian Snook, and Masuhiro Mikami for helpful discussion and careful reading of the manuscript. This work is partly supported by the Next Generation Super Computing Project, the Nanoscience Program, Japan.
Notes
Notes
1. In d (>1) dimensional systems of N particles, i denotes both the particle number and a spatial component (i.e. i = 1, …, dN).
2. We do not incorporate the auxiliary function f aux in Citation6 for clarity. In fact, the generalised Nosé–Hoover dynamics can be constructed on the basis of H gMT with f aux, which results in extra harmonic oscillator-like terms in the dynamical equations of the thermostat variables [p η l : Equation (Equation22)]. These extra terms might obviate the resonance problem in multiple time step integrators (see Citation14 for details).
3. For simplicity, we assume that all particles are attached to the same set of thermostats (e.g. ), but the discussion here can be generalised to ‘individual’ or ‘massive’ thermostatting Citation13,Citation14. In this case, the product of j and j ′ in Equations (Equation12) and (Equation13) depends on i as in Equation (Equation14).
4. h i (s) and f i (s) are, in this case, given as h i (s) = s 1 s 2 and f i (s) = s 3. u l (s), on the other hand, depends on l; u 1 = s 2 s 4, u 2 = s 5, u 3 = s 5, and u 4 = u 5 = 1 (i.e. u l consists of the thermostat variables that control the ‘lth’ thermostat). Note that s 2 is included both in s j [Equation (Equation12): h i ] and s j ′′(1) [Equation (Equation14): u l=1].