A Nosé-Hoover Thermostat Adapted to a Slab Geometry

Abstract

A Nosé-Hoover (NH) type thermostat is considered for Molecular Dynamics (MD) simulations of confined systems. This is based on a generalised velocity of the same generic form as the NH thermostat of Allen and Schmid, [Mol. Sim. 33, 21 (2007)]. An unthermostatted confined region is sandwiched between two walls which are thermostatted. No external shearing is imposed. Non-equilibrium Molecular Dynamics (NEMD) simulations were carried out of the time evolution of the wall and confined region temperature after a jump in temperature of the walls. Relaxation of the confined region temperature to the target value was found to be typically slower than that of the wall. An analysis of the system parameter dependence of the lag time, τ , and departures from what would be expected from Fourier’s law suggest that a boundary transmission heat flux bottleneck is a significant factor in the time delay. This delayed thermal equilibration would therefore become an important factor when a time-dependent (e.g., oscillatory) temperature or shearing of the walls is implemented using NEMD. Adjustments between the response time of the wall thermostat should be made compatible with that of the rest of the system, to minimise its effects on the observed behaviour.