Configurational temperature thermostat for fluids undergoing shear flow: application to liquid chlorine

Non-equilibrium 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 centre 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. In this work, we propose a thermostat based on a configurational expression for the temperature and apply it to NEMD simulations of chlorine undergoing Couette flow. The results so obtained are compared with those obtained using a COM kinetic thermostat. At equilibrium the properties of systems thermostatted in the two different ways are of course equivalent. We show that the two responses only differ far from equilibrium. In particular, we show that the formation of a string phase for extremely high shear rates is an artefact of the COM thermostat. At the largest shear rates studied with the configurational thermostat, no string phase is observed.

[1]  Antoinette J. Taylor,et al.  Thermodynamic and structural properties of liquids modelled by ‘2-Lennard-Jones centres’ pair potentials , 1977 .

[2]  Hoover,et al.  Canonical dynamics: Equilibrium phase-space distributions. , 1985, Physical review. A, General physics.

[3]  Evans,et al.  Shear thickening and turbulence in simple fluids. , 1986, Physical review letters.

[4]  Denis J. Evans,et al.  On the nonlinear Born effect , 1987 .

[5]  J. Clarke,et al.  Investigation of the homogeneous-shear nonequilibrium-molecular-dynamics method. , 1992, Physical review. A, Atomic, molecular, and optical physics.

[6]  Evans,et al.  Conditions for the existence of a reentrant solid phase in a sheared atomic fluid. , 1992, Physical review. A, Atomic, molecular, and optical physics.

[7]  P. Daivis,et al.  Comparison of constant pressure and constant volume nonequilibrium simulations of sheared model decane , 1994 .

[8]  P. Daivis,et al.  Thermostats for molecular fluids undergoing shear flow: Application to liquid chlorine , 1995 .

[9]  H. D. Cochran,et al.  Molecular dynamics simulations of the rheology of normal decane, hexadecane, and tetracosane , 1996 .

[10]  H. H. Rugh,et al.  Dynamical Approach to Temperature , 1997, chao-dyn/9701026.

[11]  J. Pablo,et al.  Rheological, thermodynamic, and structural studies of linear and branched alkanes under shear , 1997 .

[12]  R. L. Rowley,et al.  Predicting the viscosity of alkanes using nonequilibrium molecular dynamics: Evaluation of intermolecular potential models , 1997 .

[13]  Denis J. Evans,et al.  Configurational temperature: Verification of Monte Carlo simulations , 1998 .

[14]  Denis J. Evans,et al.  ON THE VALIDITY OF FOURIER'S LAW IN SYSTEMS WITH SPATIALLY VARYING STRAIN RATES , 1999 .

[15]  Evans,et al.  Configurational temperature for systems with constraints , 2000, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[16]  Microscopic expressions for the thermodynamic temperature , 1999, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[17]  Jerome Delhommelle,et al.  Configurational temperature profile in confined fluids. I. Atomic fluid , 2001 .