Trond S. Ingebrigtsen

Scaling of viscous dynamics in simple liquids: theory, simulation and experiment

Supercooled liquids are characterized by relaxation times that increase dramatically by cooling or compression. From a single assumption follows a scaling law according to which the relaxation time is a function of h() over temperature, where is the density and the function h() depends on the liquid in question. This scaling is demonstrated to work well for simulations of the Kob-Andersen binary Lennard-Jones mixture and two molecular models, as well as for the experimental results for two van der Waals liquids, dibutyl phthalate and decahydroisoquinoline. The often used power- law density scaling, h() / , is an approximation to the more general form of scaling discussed here. A thermodynamic derivation was previously given for an explicit expression for h() for liquids of particles interacting via the generalized Lennard-Jones potential. Here a statistical mechanics derivation is given, and the prediction is shown to agree very well with simulations over large density changes. Our findings effectively reduce the problem of understanding the viscous slowing down from being a quest for a function of two variables to a search for a single-variable function.Supercooled liquids are characterized by relaxation times that increase dramatically by cooling or compression. Many liquids have been shown to obey power-law density scaling, according to which the relaxation time is a function of density to some power over temperature. We show that power-law density scaling breaks down for larger density variations than usually studied. This is demonstrated by simulations of the Kob-Andersen binary Lennard-Jones mixture and two molecular models, as well as by experimental results for two van der Waals liquids. A more general form of density scaling is derived, which is consistent with results for all the systems studied. An analytical expression for the scaling function for liquids of particles interacting via generalized Lennard-Jones potentials is derived and shown to agree very well with simulations. This effectively reduces the problem of understanding the viscous slowing down from being a quest for a function of two variables to a search for a single-variable function.

Communication: Thermodynamics of condensed matter with strong pressure-energy correlations

We show that for any liquid or solid with strong correlation between its NVT virial and potential-energy equilibrium fluctuations, the temperature is a product of a function of excess entropy per particle and a function of density, T = f(s)h(ρ). This implies that (1) the systems isomorphs (curves in the phase diagram of invariant structure and dynamics) are described by h(ρ)/T = Const., (2) the density-scaling exponent is a function of density only, and (3) a Grüneisen-type equation of state applies for the configurational degrees of freedom. For strongly correlating atomic systems one has h(ρ) = ∑(n)C(n)ρ(n/3) in which the only non-zero terms are those appearing in the pair potential expanded as ν(r) = ∑(n)ν(n)r(-n). Molecular dynamics simulations of Lennard-Jones type systems confirm the theory.

Predicting how nanoconfinement changes the relaxation time of a supercooled liquid.

The properties of nanoconfined fluids can be strikingly different from those of bulk liquids. A basic unanswered question is whether the equilibrium and dynamic consequences of confinement are related to each other in a simple way. We study this question by simulation of a liquid comprising asymmetric dumbbell-shaped molecules, which can be deeply supercooled without crystallizing. We find that the dimensionless structural relaxation times-spanning six decades as a function of temperature, density, and degree of confinement-collapse when plotted versus excess entropy. The data also collapse when plotted versus excess isochoric heat capacity, a behavior consistent with the existence of isomorphs in the bulk and confined states.

Isomorphs in Model Molecular Liquids

Isomorphs are curves in the phase diagram along which a number of static and dynamic quantities are invariant in reduced units (Gnan, N.; et al. J. Chem. Phys.2009, 131, 234504). A liquid has good isomorphs if and only if it is strongly correlating, i.e., if the equilibrium virial/potential energy fluctuations are more than 90% correlated in the NVT ensemble. Isomorphs were previously discussed with a focus on atomic systems. This paper generalizes isomorphs to liquids composed of rigid molecules and study the isomorphs of systems of small rigid molecules: the asymmetric dumbbell model, a symmetric inverse power-law dumbbell, and the Lewis–Wahnström o-terphenyl (OTP) model. For all model systems, the following quantities are found to a good approximation to be invariant along an isomorph: the isochoric heat capacity, the excess entropy, the reduced molecular center-of-mass self-part of the intermediate scattering function, and the reduced molecular center-of-mass radial distribution function. In agreement with theory, we also find that an instantaneous change of temperature and density from an equilibrated state point to an isomorphic state point leads to no relaxation. The isomorphs of the Lewis–Wahnström OTP model were found to be more approximative than those of the asymmetric dumbbell model; this is consistent with the OTP model being less strongly correlating. The asymmetric dumbbell and Lewis–Wahnström OTP models each have a “master isomorph”; i.e., the isomorphs have identical shape in the virial/potential energy phase diagram.

Do the repulsive and attractive pair forces play separate roles for the physics of liquids

According to standard liquid-state theory repulsive and attractive pair forces play distinct roles for the physics of liquids. This paradigm is put into perspective here by demonstrating a continuous series of pair potentials that have virtually the same structure and dynamics, although only some of them have attractive forces of significance. Our findings reflect the fact that the motion of a given particle is determined by the total force on it, whereas the quantity usually discussed in liquid-state theory is the individual pair force.

RUMD: A general purpose molecular dynamics package optimized to utilize GPU hardware down to a few thousand particles

RUMD is a general purpose, high-performance molecular dynamics (MD) simulation package running on graphical processing units (GPUs). RUMD addresses the challenge of utilizing the many-core nature of modern GPU hardware when simulating small to medium system sizes (roughly from a few thousand up to hundred thousand particles). It has a performance that is comparable to other GPU-MD codes at large system sizes and substantially better at smaller sizes.RUMD is open-source and consists of a library written in C++ and the CUDA extension to C, an easy-to-use Python interface, and a set of tools for set-up and post-simulation data analysis. The paper describes RUMDs main features, optimizations and performance benchmarks.

Effect of Energy Polydispersity on the Nature of Lennard-Jones Liquids

In the companion paper [ Ingebrigtsen , T. S. ; Tanaka , H. J. Phys. Chem. B 2015 , 119 , 11052 ] the effect of size polydispersity on the nature of Lennard-Jones (LJ) liquids, which represent most molecular liquids without hydrogen bonds, was studied. More specifically, it was shown that even highly size polydisperse LJ liquids are Roskilde-simple (RS) liquids. RS liquids are liquids with strong correlation between constant volume equilibrium fluctuations of virial and potential energy and are simpler than other types of liquids. Moreover, it was shown that size polydisperse LJ liquids have isomorphs to a good approximation. Isomorphs are curves in the phase diagram of RS liquids along which structure, dynamics, and some thermodynamic quantities are invariant in dimensionless (reduced) units. In this paper, we study the effect of energy polydispersity on the nature of LJ liquids. We show that energy polydisperse LJ liquids are RS liquids. However, a tendency of particle segregation, which increases with the degree of polydispersity, leads to a loss of strong virial-potential energy correlation but is mitigated by increasing temperature and/or density. Isomorphs are a good approximation also for energy polydisperse LJ liquids, although particle-resolved quantities display a somewhat poorer scaling compared to the mean quantities along the isomorph.

NVU dynamics. III. Simulating molecules at constant potential energy

This is the final paper in a series that introduces geodesic molecular dynamics at constant potential energy. This dynamics is entitled NVU dynamics in analogy to standard energy-conserving Newtonian NVE dynamics. In the first two papers [T. S. Ingebrigtsen, S. Toxvaerd, O. J. Heilmann, T. B. Schrøder, and J. C. Dyre, J. Chem. Phys. 135, 104101 (2011); T. S. Ingebrigtsen, S. Toxvaerd, T. B. Schrøder, and J. C. Dyre, ibid. 135, 104102 (2011)], a numerical algorithm for simulating geodesic motion of atomic systems was developed and tested against standard algorithms. The conclusion was that the NVU algorithm has the same desirable properties as the Verlet algorithm for Newtonian NVE dynamics, i.e., it is time-reversible and symplectic. Additionally, it was concluded that NVU dynamics becomes equivalent to NVE dynamics in the thermodynamic limit. In this paper, the NVU algorithm for atomic systems is extended to be able to simulate the geodesic motion of molecules at constant potential energy. We derive an algorithm for simulating rigid bonds and test this algorithm on three different systems: an asymmetric dumbbell model, Lewis-Wahnström o-terphenyl (OTP) and rigid SPC/E water. The rigid bonds introduce additional constraints beyond that of constant potential energy for atomic systems. The rigid-bond NVU algorithm conserves potential energy, bond lengths, and step length for indefinitely long runs. The quantities probed in simulations give results identical to those of Nosé-Hoover NVT dynamics. Since Nosé-Hoover NVT dynamics is known to give results equivalent to those of NVE dynamics, the latter results show that NVU dynamics becomes equivalent to NVE dynamics in the thermodynamic limit also for molecular systems.