Configurational thermostats for molecular systems

L. Lue, O. G. Jepps, Jerome Delhommelle*, D. J. Evans

*Corresponding author for this work

    Research output: Contribution to journalArticlepeer-review

    42 Citations (Scopus)


    A new method of thermostatting non-equilibrium molecular dynamics (NEMD) simulations is described. The thermostat is based on a recently developed, entirely configurational expression for the temperature. To demonstrate this method, thermostatted NEMD simulations are performed on WCA atoms, linear, freely jointed Lennard-Jones 8-mer chains and a united-atom model of n-decane under a constant applied strain rate. The results of simulations thermostatted kinetically (the standard method) and configurationally are compared. As expected, both types of thermostat yield identical system properties for low strain rates. For higher strain rates, both thermostats yield the same qualitative dependence of system properties on applied strain rate. The great advantage of the configurational thermostat is that no a priori knowledge of the streaming velocity is required. For molecular systems and atomic systems in most flow geometries, the analytical form for the atomic streaming velocity is not known. This makes the implementation of standard kinetic thermostats highly problematic.

    Original languageEnglish
    Pages (from-to)2387-2395
    Number of pages9
    JournalMolecular Physics
    Issue number14
    Publication statusPublished - 20 Jul 2002


    Dive into the research topics of 'Configurational thermostats for molecular systems'. Together they form a unique fingerprint.

    Cite this