Comparison of thermostatting mechanisms in NVT and NPT simulations of decane under shear

Nonequilibrium molecular dynamics (NEMD) simulations play a major role in characterizing the rheological properties of fluids undergoing shear flow. However, all previous studies of flows in molecular fluids either use an “atomic” thermostat which makes incorrect assumptions concerning the streaming velocity of atoms within their constituent molecules, or they employ a center of mass kinetic (COM) thermostat which only controls the temperature of relatively few degrees of freedom (3) in complex high molecular weight compounds. In the present paper we show how recently developed configurational expressions for the thermodynamic temperature can be used to develop thermostatting mechanisms which avoid both of these problems. We propose a thermostat based on a configurational expression for the temperature and apply it to NEMD simulations of decane undergoing Couette flow at constant volume and at constant pressure. The results so obtained are compared with those obtained using a COM kinetic thermostat. At eq...