Search

The hypersphere model is a simple one-parameter model of the potential energy landscape of viscous liquids, which consists of a percolating system of hyperspheres of equal sizes randomly distributed in $R^{3N}$ where $N$ is the number of particles. We study random walks within overlapping hyperspheres in 12 to 45 dimensions, utilizing an algorithm for on-the-fly placement of the hyperspheres in conjunction with the kinetic Monte Carlo method. We find behavior typical of viscous liquids; thus decreasing the hypersphere density (corresponding to decreasing the temperature) leads to a slowing down of the dynamics by many orders of magnitude. The shape of the mean-square displacement as a function of time is found to be very similar to that of the Kob-Andersen binary Lennard-Jones mixture and the Random Barrier Model, which predicts well the frequency-dependent fluidity of nine glass-forming liquids of different chemistry [Bierwirth et al., Phys. Rev. Lett. 119, 248001 (2017)].

Isomorphs are curves in the thermodynamic phase diagram of invariant excess entropy, structure, and dynamics, while pseudoisomorphs are curves of invariant structure and dynamics, but not of the excess entropy. The latter curves have been shown to exist in molecular models with flexible bonds [A. E. Olsen et al., J. Chem. Phys. 145, 241103 (2016)]. We here present three force-based methods to trace out pseudoisomorphs based on a single configuration and test them on the asymmetric dumbbell and 10-bead Lennard-Jones chain models with bonds modeled as harmonic springs. The three methods are based on requiring that particle forces, center-of-mass forces, and torques, respectively, are invariant in reduced units. For each of the two investigated models we identify a method that works well for tracing out pseudoisomorphs, but these methods are not the same. Overall, it appears that the more internal degrees of freedom there are in the molecule studied, the less they affect the gross dynamical behavior. Moreover, the "internal" degrees of freedom (including rotation) do not appear to significantly affect the scaling behavior of the dynamical/transport coefficients provided some "quenching" is performed.

A simple model for viscous liquid dynamics is introduced. Consider the surface of the union of hyper-spheres centered at random positions inside a hypercube with periodic boundary conditions. It is argued and demonstrated by numerical simulations that at high dimensions geodetic flows on this surface is a good model for viscous liquid dynamics. It is shown that this simple model exhibits viscous dynamics for densities above the percolation threshold in $8$, $12$ and $16$ dimensions. Thus the slowing down of the dynamics, measured by the mean-squared displacement, extends to several orders of magnitude similarly to what is observed in other models for viscous dynamics. Furthermore, the shape of the mean-squared displacement is to a very good approximation the same as for the standard model in simulations of viscous liquids: the Kob-Andersen binary Lennard Jones mixture., Comment: 4 pages, 5 Figures

Isomorphs are curves in the thermodynamic phase diagram along which structure and dynamics are invariant to a good approximation. There are two main ways to trace out isomorphs, the configurational-adiabat method and the direct-isomorph-check method. Recently a new method based on the scaling properties of forces was introduced and shown to work very well for atomic systems [T. B. Schroder, Phys. Rev. Lett. 129, 245501 (2022)]. A unique feature of this method is that it only requires a single equilibrium configuration for tracing out an isomorph. We here test generalizations of this method to molecular systems and compare to simulations of three simple molecular models: the asymmetric dumbbell model of two Lennard-Jones spheres, the symmetric inverse-power-law dumbbell model, and the Lewis-Wahnstr\"om o-terphenyl model. We introduce and test two force-based and one torque-based methods, all of which require just a single configuration for tracing out an isomorph. Overall, the method based on requiring invariant center-of-mass reduced forces works best., Comment: 20 pages, 12 figures

Time-dependent dynamical properties of a fluid can not be estimated from a single configuration without performing a simulation. Here we show, however, that the scaling properties of both structure and dynamics can be predicted from a single configuration. The new method is demonstrated to work very well for equilibrium dynamics of the Kob-Andersen Binary Lennard-Jones mixture. Furthermore, the method is applied to isobaric cooling where the liquid falls out of equilibrium and forms a glass, demonstrating that the method requires neither equilibrium nor constant volume conditions to work, in contrast to existing methods., Comment: 11 pages, 5 figures

The dramatic slowing down of relaxation dynamics of liquids approaching the glass transition remains a highly debated problem, where the crux of the puzzle resides in the elusive increase of the activation barrier $\Delta E(T)$ with decreasing temperature $T$. A class of theoretical frameworks -- known as elastic models -- attribute this temperature dependence to the variations of the liquid's macroscopic elasticity, quantified by the high-frequency shear modulus $G_\infty(T)$. While elastic models find some support in a number of experimental studies, these models do not take into account the spatial structures, length scales, and heterogeneity associated with structural relaxation in supercooled liquids. Here, we propose that viscous slowing down is controlled by a mesoscopic elastic stiffness $\kappa(T)$, defined as the characteristic stiffness of response fields to local dipole forces in the liquid's underlying inherent states. First, we show that $\kappa(T)$ -- which is intimately related to the energy and length scales characterizing quasilocalized, nonphononic excitations in glasses -- increases more strongly with decreasing $T$ than the macroscopic inherent state shear modulus $G(T)$ in several computer liquids. Second, we show that the simple relation $\Delta E(T)\propto\kappa(T)$ holds remarkably well for some computer liquids, implying a direct connection between the liquid's underlying mesoscopic elasticity and enthalpic energy barriers. On the other hand, we show that for other computer liquids, the above relation fails. Finally, we provide strong evidence that what distinguishes computer liquids in which the $\Delta E(T) \propto \kappa(T)$ relation holds, from those in which it does not, is that the latter feature highly granular potential energy landscapes, where many sub-basins separated by low activation barriers exist. [Rest of abstract abridged], Comment: 15 pages, 12 figures. Accepted manuscript

It was recently shown that the real part of the frequency-dependent fluidity for several glass-forming liquids of different chemistry conforms to the prediction of the random barrier model (RBM) devised for ac electrical conduction in disordered solids [S. P. Bierwirth \textit{et al.}, Phys. Rev. Lett. {\bf 119}, 248001 (2017)]. Inspired by these results we introduce a crystallization-resistant modification of the Kob-Andersen binary Lennard-Jones mixture for which the results of extensive graphics-processing unit (GPU)-based molecular-dynamics simulations are presented. We find that the low-temperature mean-square displacement is fitted well by the RBM prediction, which involves no shape parameters. This finding highlights the challenge of explaining why a simple model based on hopping of non-interacting particles in a fixed random energy landscape can reproduce the complex and highly cooperative dynamics of glass-forming liquids.

The exponentially repulsive EXP pair potential defines a system of particles in terms of which simple liquids' quasiuniversality may be explained [A. K. Bacher et al., Nat. Commun. 5, 5424 (2014); J. C. Dyre, J. Phys. Condens. Matter 28, 323001 (2016)]. This paper and its companion present a detailed simulation study of the EXP system. Here we study how structure monitored via the radial distribution function and dynamics monitored via the mean-square displacement as a function of time evolve along the system's isotherms and isochores. The focus is on the gas and liquid phases, which are distinguished pragmatically by the absence or presence of a minimum in the radial distribution function above its first maximum. An NVU-based proof of quasiuniversality is presented, and quasiuniversality is illustrated by showing that the structure of the Lennard-Jones system at four selected state points is well approximated by those of EXP pair-potential systems with the same reduced diffusion constant. The companion paper studies the EXP system's isomorphs, focusing also on the gas and liquid phases.

This paper continues the investigation of the exponentially repulsive EXP pair-potential system of Paper I with a focus on isomorphs in the low-temperature gas and liquid phases. As expected from the EXP system's strong virial potential-energy correlations, the system's reduced-unit structure and dynamics are isomorph invariant to a good approximation. Three methods for generating isomorphs are compared: the small-step method that is exact in the limit of small density changes and two versions of the direct-isomorph-check method that allows for much larger density changes. Results from approximate methods are compared to those of the small-step method for each of three isomorphs generated by 230 one percent density increases, covering in all one decade of density variation. Both approximate methods work well.

Crucial to gaining control over crystallisation in multicomponent materials or accurately modelling rheological behaviour of magma flows is to understand the mechanisms by which crystal nuclei form. The microscopic nature of such nuclei, however, makes this extremely hard in experiments, while computer simulations have hitherto been hampered by their short timescales and small system sizes due to limited computational power. Here we use highly-efficient GPU simulation techniques to access system sizes around 100 times larger than previous studies. This makes it possible to elucidate the nucleation mechanism in a well-studied binary glassformer. We discover that the supercooled liquid is inherently unstable for system sizes of 10,000 particles and larger. This effect is due to compositional fluctuations leading to regions comprised of large particles only which rapidly nucleate. We argue that this mechanism provides a minimum rate of crystallisation in mixtures in general, and use our results to stabilise a model binary mixture and predict glassforming ability for the CuZr metallic glassformer., Comment: 11 pages

The binary Kob-Andersen (KA) Lennard-Jones mixture is the standard model for computational studies of viscous liquids and the glass transition. For very long simulations the viscous KA system crystallizes, however, by phase separating into a pure A particle phase forming an FCC crystal. We present the thermodynamic phase diagram for KA-type mixtures consisting of up to 50\% small (B) particles showing, in particular, that the melting temperature of the standard KA system at liquid density $1.2$ is $1.028(3)$ in A particle Lennard-Jones units. At large B particle concentrations the system crystallizes into the CsCl crystal structure. The eutectic corresponding to the FCC and CsCl structures is cut-off in a narrow interval of B particle concentrations around 26\% at which the bipyramidal orthorhombic ${\rm PuBr_3}$ structure is the thermodynamically stable phase. The melting temperature's variation with B particle concentration at two other pressures, as well as at the constant density $1.2$, is estimated from the simulations at pressure $10.19$ using isomorph theory. Our data demonstrate approximate identity between the melting temperature and the onset temperature below which viscous dynamics appears. Finally, the nature of the solid-liquid interface is briefly discussed.

Computer simulations show that liquids of molecules with harmonic intramolecular bonds may have "pseudoisomorphic" lines of approximately invariant dynamics in the thermodynamic phase diagram. We demonstrate that these lines can be identified by requiring scale invariance of the inherent-structure reduced-unit low-frequency vibrational spectrum evaluated for a single equilibrium configuration. This rationalizes why excess-entropy scaling, density scaling, and isochronal superposition apply for many liquids with internal degrees of freedom.

The recent theoretical prediction by Maimbourg and Kurchan [arXiv:1603.05023] that for regular pair-potential systems the virial potential-energy correlation coefficient increases towards unity as the dimension $d$ goes to infinity is investigated for the standard 12-6 Lennard-Jones fluid. This is done by computer simulations for $d=2,3,4$ going from the critical point along the critical isotherm/isochore to higher density/temperature. In all cases the virial potential-energy correlation coefficient increases significantly. For a given density and temperature relative to the critical point, with increasing number of dimension the Lennard-Jones system conforms better to the hidden-scale-invariance property characterized by high virial potential-energy correlations (a property that leads to the existence of isomorphs in the thermodynamic phase diagram, implying that it becomes effectively one-dimensional in regard to structure and dynamics). The present paper also gives the first numerical demonstration of isomorph invariance of structure and dynamics in four dimensions. Our findings emphasize the need for a universally applicable $1/d$ expansion in liquid-state theory; we conjecture that the systems known to obey hidden scale invariance in three dimensions are those for which the yet-to-be-developed $1/d$ expansion converges rapidly.

Many liquids have curves (isomorphs) in their phase diagrams along which structure, dynamics, and some thermodynamic quantities are invariant in reduced units. A substantial part of their phase diagrams is thus effectively one dimensional. The shape of these isomorphs is described by a material-dependent function of density $h(\rho)$, which for real liquids is well approximated by a power law $\rho^\gamma$. However, in simulations, a power law is not adequate when density changes are large; typical models such as the Lennard-Jones liquid show that $\gamma(\rho) \equiv \mathrm{d} \ln h(\rho)/\mathrm{d} \ln \rho$ is a decreasing function of density. This paper presents results from computer simulations using a new pair potential that diverges at a nonzero distance and can be tuned to give a more realistic shape of $\gamma(\rho)$. Our results indicate that the finite size of molecules is an important factor to take into account when modeling liquids over a large density range., Comment: 5 pages, 4 figures

The invariance of several structural and dynamical properties of the Lennard-Jones (LJ) system along the freezing and melting lines is interpreted in terms of the isomorph theory. First the freezing/melting lines for LJ system are shown to be approximated by isomorphs. Then we show that the invariants observed along the freezing and melting isomorphs are also observed on other isomorphs in the liquid and crystalline phase. Structure is probed by the radial distribution function and the structure factor and dynamics is probed by the mean-square displacement, the intermediate scattering function, and the shear viscosity. Studying these properties by reference to the isomorph theory explains why known single-phase melting criteria holds, e.g., the Hansen-Verlet and the Lindemann criterion, and why the Andrade equation for the viscosity at freezing applies, e.g., for most liquid metals. Our conclusion is that these empirical rules and invariants can all be understood from the isomorph theory and that the invariants are not peculiar to the freezing and melting lines, but hold along all isomorphs., Comment: 21 pg, 12 figures Accepted from PCCP (Physical Chemistry Chemical Physics)

The previous paper [Veldhorst et al., J. Chem. Phys. 141, 054904 (2014)] demonstrated that the isomorph theory explains the scaling properties of a liquid of flexible chains consisting of ten Lennard-Jones particles connected by rigid bonds. We here investigate the same model with harmonic bonds. The introduction of harmonic bonds almost completely destroys the correlations in the equilibrium fluctuations of the potential energy and the virial. According to the isomorph theory, if these correlations are strong a system has isomorphs, curves in the phase diagram along which structure, dynamics and the excess entropy are invariant. The Lennard-Jones chain liquid with harmonic bonds does have curves in the phase diagram along which the structure and dynamics are invariant. The excess entropy is not invariant on these curves, which we refer to as "pseudoisomorphs". In particular this means that Rosenfeld's excess-entropy scaling (the dynamics being a function of excess entropy only) does not apply for the Lennard-Jones chain with harmonic bonds., Comment: 11 pages, 11 figures

RUMD is a general purpose, high-performance molecular dynamics (MD) simulation package running on graphical processing units (GPU's). 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 RUMD's main features, optimizations and performance benchmarks., Comment: 22 pages. Resubmission to SciPost Physics

This paper shows that several known properties of the Yukawa system can be derived from the isomorph theory, which applies to any system that has strong correlations between its virial and potential-energy equilibrium fluctuations. Such "Roskilde-simple" systems have a simplified thermodynamic phase diagram deriving from the fact that they have curves (isomorphs) along which structure and dynamics in reduced units are invariant to a good approximation. We show that the Yukawa system has strong virial potential-energy correlations and identify its isomorphs by two different methods. One method, the so-called direct isomorph check, identifies isomorphs numerically from jumps of relatively small density changes (here 10%). The second method identifies isomorphs analytically from the pair potential. The curves obtained by the two methods are close to each other; these curves are confirmed to be isomorphs by demonstrating the invariance of the radial distribution function, the static structure factor, the mean-square displacement as a function of time, and the incoherent intermediate scattering function. Since the melting line is predicted to be an isomorph, the theory provides a derivation of a known approximate analytical expression for this line in the temperature-density phase diagram. The paper's results give the first demonstration that the isomorph theory can be applied to systems like dense colloidal suspensions and strongly coupled dusty plasmas., Comment: 12 pages, 12 figures

It has been known for a long time that many simple liquids have surprisingly similar structure as quantified, e.g., by the radial distribution function. A much more recent realization is that the dynamics are also very similar for a number of systems with quite different pair potentials. Systems with such non-trivial similarities are generally referred to as "quasi-universal". From the fact that the exponentially repulsive pair potential has strong virial potential-energy correlations in the low-temperature part of its thermodynamic phase diagram, we here show that a liquid is quasi-universal if its pair potential can be written approximately as a sum of exponential terms with numerically large prefactors. Based on evidence from the literature we moreover conjecture the converse, i.e., that quasi-universality only applies for systems with this property.

Asphaltene and smaller aromatic molecules tend to form linear nanoaggregates in bitumen.Over the years bitumen undergoes chemical aging and during this process, the size of the nanoaggregate increases. This increase is associated with an increase in viscosity and brittleness of the bitumen, eventually leading to road deterioration. This paper focuses on understanding the mechanisms behind nanoaggregate size and stability. We used molecular dynamics simulations to quantify the probability of having a nanoaggregate of a given size in the stationary regime. To model this complicated behavior, we chose first to consider the simple case where only asphaltene molecules are counted in a nanoaggregate. We used a master equation approach and a related statistical mechanics model. The linear asphaltene nanoaggregates behave as a rigid linear chain. The most complicated case where all aromatic molecules are counted in a nanoaggregate is then discussed. The linear aggregates where all aromatic molecules are counted seem to behave as a flexible linear chain., Comment: 31 pages, 9 figures

The theory of isomorphs is reformulated by defining Roskilde-simple systems (those with isomorphs) by the property that the order of the potential energies of configurations at one density is maintained when these are scaled uniformly to a different density. Isomorphs remain curves in the thermodynamic phase diagram along which structure, dynamics, and excess entropy are invariant, implying that the phase diagram is effectively one-dimensional with respect to many reduced-unit properties. In contrast to the original formulation of the isomorph theory, however, the density-scaling exponent is not exclusively a function of density and the isochoric heat capacity is not an exact isomorph invariant. A prediction is given for the latter quantity's variation along the isomorphs. Molecular dynamics simulations of the Lennard-Jones and Lennard-Jones Gaussian systems validate the new approach.

This paper shows by computer simulations that some crystalline systems have curves in their thermodynamic phase diagrams, so-called isomorphs, along which structure and dynamics in reduced units are invariant to a good approximation. The crystals are studied in a classical-mechanical framework, which is generally a good description except significantly below melting. The existence of isomorphs for crystals is validated by simulations of particles interacting via the Lennard-Jones pair potential arranged into a face-centered cubic (FCC) crystalline structure; the slow vacancy-jump dynamics of a defective FCC crystal is also shown to be isomorph invariant. In contrast, a NaCl crystal model does not exhibit isomorph invariances. Other systems simulated, though in less detail, are the Wahnstrom binary Lennard-Jones crystal with the ${\rm MgZn_2}$ Laves crystal structure, monatomic FCC crystals of particles interacting via the Buckingham pair potential and via a novel purely repulsive pair potential diverging at a finite separation, an ortho-terphenyl molecular model, and SPC/E hexagonal ice. Except for NaCl and ice, the crystals simulated all have isomorphs. Based on these findings and previous simulations of liquid models, we conjecture that crystalline solids with isomorphs include most or all formed by atoms or molecules interacting via metallic or van der Waals forces, whereas covalently- or hydrogen-bonded crystals are not expected to have isomorphs. Crystals of ions or dipolar molecules constitute a limiting case for which isomorphs are only expected when the Coulomb interactions are relatively weak. We briefly discuss the consequences of the findings for theories of melting and crystallization.

Koperwas {\it et al.} showed in a recent paper, Phys. Rev. Lett. {\bf 111}, 125701 (2013), that the dynamic susceptibility $\chi_4$ as estimated by dielectric measurements for certain glass-forming liquids decreases substantially with increasing pressure along a curve of constant relaxation time. This observation is at odds with other measures of dynamics being invariant and seems to pose a problem for theories of glass formation. We show that this variation is in fact consistent with predictions for liquids with hidden scale invariance: measures of dynamics at constant volume are invariant along isochrones, called isomorphs in such liquids, but contributions to fluctuations from long-wavelength fluctuations can vary. This is related to the known non-invariance of the isothermal bulk modulus. Considering the version of $\chi_4$ defined for the NVT ensemble, data from simulations of a binary Lennard-Jones liquid show in fact a slight increase with increasing density. This is a true departure from the formal invariance expected for this quantity., Comment: Expanded to include simulation data, and accepted for publication in Phys. Rev. E

This paper investigates two conjectures for calculating the density dependence of the density-scaling exponent of a single-component, pair-potential liquid with strong virial potential-energy correlations. The first conjecture gives an analytical expression for the density-scaling exponent directly in terms of the pair potential. The second conjecture is a refined version of this, which involves the most likely nearest-neighbor distance determined from the pair-correlation function. The two conjectures for the density-scaling exponent are tested by simulations of three systems, one of which is the standard Lennard-Jones liquid. While both expressions give qualitatively correct results, the second expression is more accurate.

We derive exact results for the rate of change of thermodynamic quantities, in particular the configurational specific heat at constant volume, $C_V$, along configurational adiabats (curves of constant excess entropy $S_{\textrm{ex}}$). Such curves are designated isomorphs for so-called Roskilde liquids, in view of the invariance of various structural and dynamical quantities along them. Their slope in a double logarithmic representation of the density-temperature phase diagram, $\gamma$ can be interpreted as one third of an effective inverse power-law potential exponent. We show that in liquids where $\gamma$ increases (decreases) with density, the contours of $C_V$ have smaller (larger) slope than configurational adiabats. We clarify also the connection between $\gamma$ and the pair potential. A fluctuation formula for the slope of the $C_V$-contours is derived. The theoretical results are supported with data from computer simulations of two systems, the Lennard-Jones fluid and the Girifalco fluid. The sign of $d\gamma/d\rho$ is thus a third key parameter in characterizing Roskilde liquids, after $\gamma$ and the virial-potential energy correlation coefficient $R$. To go beyond isomorph theory we compare invariance of a dynamical quantity, the self-diffusion coefficient along adiabats and $C_V$-contours, finding it more invariant along adiabats.

The isomorph theory provides an explanation for the so-called power law density scaling which has been observed in many molecular and polymeric glass formers, both experimentally and in simulations. Power law density scaling (relaxation times and transport coefficients being functions of $\rho^{\gamma_S}/T$, where $\rho$ is density, $T$ is temperature, and $\gamma_S$ is a material specific scaling exponent) is an approximation to a more general scaling predicted by the isomorph theory. Furthermore, the isomorph theory provides an explanation for Rosenfeld scaling (relaxation times and transport coefficients being functions of excess entropy) which has been observed in simulations of both molecular and polymeric systems. Doing molecular dynamics simulations of flexible Lennard-Jones chains (LJC) with rigid bonds, we here provide the first detailed test of the isomorph theory applied to flexible chain molecules. We confirm the existence of isomorphs, which are curves in the phase diagram along which the dynamics is invariant in the appropriate reduced units. This holds not only for the relaxation times but also for the full time dependence of the dynamics, including chain specific dynamics such as the end-to-end vector autocorrelation function and the relaxation of the Rouse modes. As predicted by the isomorph theory, jumps between different state points on the same isomorph happen instantaneously without any slow relaxation. Since the LJC is a simple coarse-grained model for alkanes and polymers, our results provide a possible explanation for why power-law density scaling is observed experimentally in alkanes and many polymeric systems. The theory provides an independent method of determining the scaling exponent, which is usually treated as a empirical scaling parameter., Comment: 9 pages, 10 figures

We propose a four-component molecular model of bitumen. The model includes realistic chemical constituents and introduces a coarse-graining level that suppresses the highest frequency modes. Molecular dynamics simulations of the model are being carried out using Graphic-Processor-Units based software in time spans in order of microseconds, and this enables the study of slow relaxation processes characterizing bitumen. This paper focuses on the high-temperature dynamics as expressed through the mean-square displacement, the stress autocorrelation function, and rotational relaxation. The diffusivity of the individual molecules changes little as a function of temperature and reveals distinct dynamical time scales as a result of the different constituents in the system. Different time scales are also observed for the rotational relaxation. The stress autocorrelation function features a slow non-exponential decay for all temperatures studied. From the stress autocorrelation function, the shear viscosity and shear modulus are evaluated at the highest temperature, showing a viscous response at frequencies below 100 MHz. The model predictions of viscosity and diffusivities are compared to experimental data, giving reasonable agreement. The model shows that the asphaltene, resin and resinous oil tend to form nano-aggregates. The characteristic dynamical relaxation time of these aggregates is different from the homogeneously distributed parts of the system, leading to strong dynamical heterogeneity., Comment: 10 pages, 12 figures

We study chemical aging in "COOEE bitumen" using molecular dynamic simulations. This model bitumen is composed of four realistic molecule types: saturated hydrocarbon, resinous oil, resin, and asphaltene. The aging reaction is modelled by the chemical reaction: "2 resins $\rightarrow$ 1 asphaltene". Molecular dynamic simulations of four bitumen compositions, obtained by a repeated application of the aging reaction, are performed. The stress autocorrelation function, the fluid structure, the rotational dynamics of the plane aromatic molecules, and the diffusivity of each molecule, are determined for the four different compositions. The aging reaction causes a significant dynamics slowdown, which is correlated to the aggregation of asphaltene molecules in larger and dynamically slower nanoaggregates. Finally, a detailed description of the role of each molecule types in the aggregation and aging processes is given., Comment: 36 pages, 15 figures

We investigate the accuracy of the expression of Rosenfeld and Tarazona (RT) for the excess isochoric heat capacity, C_V^{ex} \propto T^{-2/5}, for eighteen model liquids. Previous investigations have reported no unifying features of breakdown for the RT expression. Here liquids with different stoichiometric composition, molecular topology, chemical interactions, degree of undercooling, and environment are investigated. We find that the RT expression is a better approximation for liquids with strong correlations between equilibrium fluctuations of virial and potential energy, i.e., Roskilde simple liquids [Ingebrigtsen et al., Phys. Rev. X 2, 011011 (2012)]. This observation holds even for molecular liquids under severe nanoscale confinement, the physics of which is completely different from the original RT bulk hard-sphere fluid arguments. The density dependence of the specific heat is predicted from the isomorph theory for Roskilde simple liquids, which in combination with the RT expression provides a complete description of the specific heat's density and temperature dependence., Comment: 6 pages

Non-equilibrium molecular dynamics simulations were performed to study the thermodynamic, structural, and dynamical properties of the single-component Lennard-Jones and the Kob-Andersen binary Lennard-Jones liquids. Both systems are known to be strongly correlating, i.e., have strong correlations between equilibrium thermal fluctuations of virial and potential energy. Such systems have good isomorphs, i.e., curves in the thermodynamic phase diagram along which structural, dynamical, and some thermodynamic quantities are invariant when expressed in reduced units. The SLLOD equations of motion were used to simulate Couette shear flows of the two systems. We show analytically that these equations are isomorph invariant provided the reduced strain rate is fixed along the isomorph. Since isomorph invariance is generally only approximate, a range of shear rates were simulated to test for the predicted invariance, covering both the linear and non-linear regimes. For both systems, when represented in reduced units the radial distribution function and the intermediate scattering function collapse for state points that are isomorphic. The strain-rate dependence of the viscosity, which exhibits shear thinning, is also invariant along an isomorph. Our results extend the isomorph concept to the state-state non-equilibrium situation of a shear flow, in which the phase diagram is three dimensional because the shear rate defines the third dimension.

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.

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., Comment: 7 pages, 5 figures

Simple liquids are traditionally defined as many-body systems of classical particles interacting via radially symmetric pair potentials. We suggest that a simple liquid should be defined instead by the property of having strong correlation between virial and potential energy equilibrium fluctuations in the NVT ensemble. There is considerable overlap between the two definitions, but also some notable differences. For instance, in the new definition simplicity is not a property of the intermolecular potential only because a liquid is usually only strongly correlating in part of its phase diagram. Moreover, according to the new definition not all simple liquids are atomic (i.e., with radially symmetric pair potentials) and not all atomic liquids are simple. The main part of the paper motivates the new definition of liquid simplicity by presenting evidence that a liquid is strongly correlating if and only if its intermolecular interactions may be ignored beyond the first coordination shell (FCS). This is demonstrated by NVT simulations of structure and dynamics of 15 atomic and molecular model liquids with a shifted-forces cutoff placed at the first minimum of the radial distribution function. No proof is given that the chemical characterization follows from the strong correlation property, but it is shown to be consistent with the existence of isomorphs in strongly correlating liquids' phase diagram. Finally, we note that the FCS characterization of simple liquids calls into question the basis for standard perturbation theory, according to which the repulsive and attractive forces play fundamentally different roles for the physics of liquids.

This paper compares the Wolf method to the shifted forces (SF) method for efficient computer simulation of isotropic systems interacting via Coulomb forces, taking results from the Ewald summation method as representing the true behavior. We find that for the Hansen-McDonald molten salt model the SF approximation overall reproduces the structural and dynamical properties as accurately as does the Wolf method. It is shown that the optimal Wolf damping parameter depends on the property in focus, and that neither the potential energy nor the radial distribution function are useful measures for the convergence of theWolf method to the Ewald summation method. The SF approximation is also tested for the SPC/Fw model of liquid water at room temperature, showing good agreement with both the Wolf and the particle mesh Ewald methods; this confirms previous findings [Fennell & Gezelter, J. Chem. Phys. {\bf 124}, 234104 (2006)]. Beside its conceptualsimplicity the SF approximation implies a speed-up of a factor 2 to 3 compared to the Wolf method (which is in turn much faster than the Ewald method).

Isomorphs are curves in the phase diagram along which a number of static and dynamic quantities are invariant in reduced units. A liquid has good isomorphs if and only if it is strongly correlating, i.e., the equilibrium virial/potential energy fluctuations are more than 90% correlated in the NVT ensemble. This paper generalizes isomorphs to liquids composed of rigid molecules and study the isomorphs of two systems of small rigid molecules, the asymmetric dumbbell model and the Lewis-Wahnstrom OTP model. In particular, for both systems we find that the isochoric heat capacity, the excess entropy, the reduced molecular center-of-mass self part of the intermediate scattering function, the reduced molecular center-of-mass radial distribution function to a good approximation are invariant along an isomorph. In agreement with theory, we also find that an instantaneous change of temperature and density from an equilibrated state point to another isomorphic state point leads to no relaxation. The isomorphs of the Lewis-Wahnstrom OTP model were found to be more approximative than those of the asymmetric dumbbell model, which is consistent with the OTP model being less strongly correlating. For both models we find "master isomorphs", i.e., isomorphs have identical shape in the virial/potential energy phase diagram., Comment: 20 pages

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(\rho)$. This implies that 1) the system's isomorphs (curves in the phase diagram of invariant structure and dynamics) are described by $h(\rho)/T={\rm Const.}$, 2) the density-scaling exponent is a function of density only, 3) a Gr{\"u}neisen-type equation of state applies for the configurational degrees of freedom. For strongly correlating atomic systems one has $h(\rho)=\sum_nC_n\rho^{n/3}$ in which the only non-zero terms are those appearing in the pair potential expanded as $v(r)=\sum_n v_n r^{-n}$. Molecular dynamics simulations of Lennard-Jones type systems confirm the theory.

It is demonstrated by molecular dynamics simulations that liquids interacting via the Buckingham potential are strongly correlating, i.e., have regions of their phase diagram where constant-volume equilibrium fluctuations in the virial and potential energy are strongly correlated. A binary Buckingham liquid is cooled to a viscous phase and shown to have isomorphs, which are curves in the phase diagram along which structure and dynamics in appropriate units are invariant to a good approximation. To test this, the radial distribution function, and both the incoherent and coherent intermediate scattering function are calculated. The results are shown to reflect a hidden scale invariance; despite its exponential repulsion the Buckingham potential is well approximated by an inverse power-law plus a linear term in the region of the first peak of the radial distribution function. As a consequence the dynamics of the viscous Buckingham liquid is mimicked by a corresponding model with purely repulsive inverse-power-law interactions. The results presented here closely resemble earlier results for Lennard-Jones type liquids, demonstrating that the existence of strong correlations and isomorphs does not depend critically on the mathematical form of the repulsion being an inverse power law., Comment: 7 pages, 9 figures

We present molecular dynamics simulations studying the influence of pressure on the correlation between the constant-volume thermal equilibrium fluctuations of virial W and potential energy U, focusing on liquids that are not strongly correlating at low pressure, i.e., do not have a WU correlation coefficient above 0.9. The systems studied are the two hydrogen-bonded liquids GROMOS methanol and TIT5P water, the ionic liquid defined by a united-atom model of the 1-butyl-3-methyl-imidazolium nitrate and, for reference, the standard single-component Lennard-Jones liquid. The simulations were performed for pressures varying from 0 GPa to 10 GPa. For all systems studied we find that the virial / potential energy correlation increases with increasing pressure. This suggests that if crystallization is avoided, all liquids become strongly correlating at sufficiently high pressure., Comment: Brief report

An algorithm is derived for computer simulation of geodesics on the constant potential-energy hypersurface of a system of N classical particles. First, a basic time-reversible geodesic algorithm is derived by discretizing the geodesic stationarity condition and implementing the constant potential energy constraint via standard Lagrangian multipliers. The basic NVU algorithm is tested by single-precision computer simulations of the Lennard-Jones liquid. Excellent numerical stability is obtained if the force cutoff is smoothed and the two initial configurations have identical potential energy within machine precision. Nevertheless, just as for NVE algorithms, stabilizers are needed for very long runs in order to compensate for the accumulation of numerical errors that eventually lead to "entropic drift" of the potential energy towards higher values. A modification of the basic NVU algorithm is introduced that ensures potential-energy and step-length conservation; center-of-mass drift is also eliminated. Analytical arguments confirmed by simulations demonstrate that the modified NVU algorithm is absolutely stable. Finally, simulations show that the NVU algorithm and the standard leap-frog NVE algorithm have identical radial distribution functions for the Lennard-Jones liquid.

In the companion paper [Ingebrigtsen et al., arXiv:1012.3447] an algorithm was developed for tracing out a geodesic curve on the constant-potential-energy hypersurface. Here simulations of this NVU dynamics are compared to results for four other dynamics, both deterministic and stochastic. First, NVU dynamics is compared to the standard energy-conserving Newtonian NVE dynamics by simulations of the Kob-Andersen binary Lennard-Jones liquid, its WCA version (i.e., with cut-off's at the pair potential minima), and the Gaussian Lennard-Jones liquid. We find identical results for all quantities probed: radial distribution functions, incoherent intermediate scattering functions, and mean-square displacement as function of time. Arguments are then presented for the equivalence of NVU and NVE dynamics in the thermodynamic limit; in particular to leading order in 1/N these two dynamics give identical time-autocorrelation functions. In the final section NVU dynamics is compared to Monte Carlo dynamics, to a diffusive dynamics of small-step random walks on the constant-potential-energy hypersurface, and to Nose-Hoover NVT dynamics. If time is scaled for the two stochastic dynamics to make their single-particle diffusion constants identical to those of NVE dynamics, the simulations show that all five dynamics are equivalent at low temperatures except at short times.

A well-known result of liquid state theory is that the structure of dense fluids is mainly determined by repulsive forces. The WCA potential, which cuts intermolecular potentials at their minima, is therefore often used as a reference. However, this reference gives quite wrong results for the viscous dynamics of the Kob-Andersen binary Lennard-Jones liquid [Berthier and Tarjus, Phys. Rev. Lett. 103, 170601 (2009)]. We show that repulsive inverse-power law potentials provide a useful reference for this liquid by reproducing its structure, dynamics, and isochoric heat capacity.

This series of papers is devoted to identifying and explaining the properties of strongly correlating liquids, i.e., liquids with more than 90% correlation between their virial W and potential energy U fluctuations in the NVT ensemble. Paper IV [N. Gnan et al., J. Chem. Phys. v131, 234504 (2009)] showed that strongly correlating liquids have "isomorphs", which are curves in the phase diagram along which structure, dynamics, and some thermodynamic properties are invariant in reduced units. In the present paper, using the fact that reduced-unit radial distribution functions are isomorph invariant, we derive an expression for the shapes of isomorphs in the WU phase diagram of generalized Lennard-Jones systems of one or more types of particles. The isomorph shape depends only on the Lennard-Jones exponents; thus all isomorphs of standard Lennard-Jones systems (with exponents 12 and 6) can be scaled onto to a single curve. Two applications are given. One is testing the prediction that the solid-liquid coexistence curve follows an isomorph by comparing to recent simulations by Ahmed and Sadus [J. Chem. Phys. v131, 174504 (2009)]. Excellent agreement is found on the liquid side of the coexistence, whereas the agreement is worse on the solid side. A second application is the derivation of an approximate equation of state for generalized Lennard-Jones systems by combining the isomorph theory with the Rosenfeld-Tarazona expression for the temperature dependence of potential energy on isochores. It is shown that the new equation of state agrees well with simulations., Comment: 12 pages, 14 figures, Section on solid-liquid coexistence expanded

The liquid structure of a glass-forming binary alloy is studied using molecular dynamics simulations. The analysis combines common neighbour analysis with the geometrical approach of Frank and Kasper to establish that the supercooled liquid contains extended clusters characterised by the same short range order as the crystal. Fluctuations in these clusters exhibit strong correlations with fluctuations in the inherent structure energy. The steep increase in the heat capacity on cooling is, thus, directly coupled to the growing fluctuations of the Frank-Kasper clusters. The relaxation of particles in the clusters dominates the slow tail of the self-intermediate scattering function.

In this third paper of the series, which started with [N. P. Bailey et al., J. Chem. Phys. 129, 184507 and 184508 (2008)], we continue the development of the theoretical understanding of strongly correlating liquids - those whose instantaneous potential energy and virial are strongly correlated in their thermal equilibrium fluctuations at constant volume. The existence of such liquids was detailed in previous work which identified them, based on computer simulations, as a large class of liquids, including van der Waals liquids but not, e.g., hydrogen-bonded liquids. We here discuss the following: (1) The scaling properties of inverse power-law and extended inverse power-law potentials (the latter include a linear term which "hides" the approximate scale invariance); (2) results from computer simulations of molecular models concerning out-of-equilibrium conditions; (3) ensemble dependence of the virial / potential energy correlation coefficient; (4) connection to the Gruneisen parameter; (5) interpretation of strong correlations in terms of the energy-bond formalism.

This paper is the fourth in a series devoted to identifying and explaining the properties of strongly correlating liquids, i.e., liquids where virial and potential energy correlate better than 90% in their thermal equilibrium fluctuations in the NVT ensemble. For such liquids we here introduce the concept of "isomorphic" curves in the state diagram. A number of thermodynamic, static, and dynamic isomorph invariants are identified. These include the excess entropy, the isochoric specific heat, reduced-unit static and dynamic correlation functions, as well as reduced-unit transport coefficients. The dynamic invariants apply for both Newtonian and Brownian dynamics. It is shown that after a jump between isomorphic state points the system is instantaneously in thermal equilibrium; consequences of this for generic aging experiments are discussed. Selected isomorph predictions are validated by computer simulations of the Kob-Andersen binary Lennard-Jones mixture, which is a strongly correlating liquid. The final section of the paper relates the isomorph concept to phenomenological melting rules, Rosenfeld's excess entropy scaling, Young and Andersen's approximate scaling principle, and the two-order parameter maps of Debenedetti and coworkers. This section also shows how the existence of isomorphs implies an "isomorph filter" for theories for the non-Arrhenius temperature dependence of viscous liquids' relaxation time, as well as explains isochronal superposition for strongly correlating viscous liquids.

A liquid is termed strongly correlating if its virial and potential energy thermal equilibrium fluctuations in the NVT ensemble are more than 90% correlated [Phys. Rev. Lett. 100, 015701 (2008)]. The fluctuations of a strongly correlating liquid are well approximated by those of an inverse power-law intermolecular potential. Building on this fact we here define "isomorphic lines" in the state diagram of a strongly correlating liquid. It is shown from computer simulations of the Kob-Andersen binary Lennard-Jones liquid that no aging is associated with jumps between two isomorphic points. Isomorphic state points have the same excess entropy, the same reduced average relaxation time, the same (reduced) dynamics, and the same scaled radial distribution functions. Finally we calculate the equation for isomorphs in the virial / potential energy diagram for Lennard-Jones type liquids and show that all such 12-6 liquids have the same isomorphs; these may be scaled into a master isomorph.

It was recently shown by computer simulations that a large class of liquids exhibits strong correlations in their thermal fluctuations of virial and potential energy [Pedersen et al., Phys. Rev. Lett. 100, 015701 (2008)]. Among organic liquids the class of strongly correlating liquids includes van der Waals liquids, but excludes ionic and hydrogen-bonding liquids. The present note focuses on the density scaling of strongly correlating liquids, i.e., the fact their relaxation time tau at different densities rho and temperatures T collapses to a master curve according to the expression tau propto F(rho^gamma/T) [Schroder et al., arXiv:0803.2199]. We here show how to calculate the exponent gamma from bulk modulus and specific heat data, either measured as functions of frequency in the metastable liquid or extrapolated from the glass and liquid phases to a common temperature (close to the glass transition temperature). Thus an exponent defined from the response to highly nonlinear parameter changes may be determined from linear response measurements.

Computer simulations recently revealed that several liquids exhibit strong correlations between virial and potential energy equilibrium fluctuations in the NVT ensemble [U. R. Pedersen {\it et al.}, Phys. Rev. Lett. {\bf 100}, 015701 (2008)]. In order to investigate whether these correlations are present also far from equilibrium constant-volume aging following a temperature down jump from equilibrium was simulated for two strongly correlating liquids, an asymmetric dumbbell model and Lewis-Wahnstr{\"o}m OTP, as well as for SPC water that is not strongly correlating. For the two strongly correlating liquids virial and potential energy follow each other closely during the aging towards equilibrium. For SPC water, on the other hand, virial and potential energy vary with little correlation as the system ages towards equilibrium. Further proof that strong pressure-energy correlations express a configuration space property comes from monitoring pressure and energy during the crystallization (reported here for the first time) of supercooled Lewis-Wahnstr{\"o}m OTP at constant temperature.

Results from molecular dynamics simulations of two viscous molecular model liquids -- the Lewis-Wahnstrom model of ortho-terphenyl and an asymmetric dumbbell model -- are reported. We demonstrate that the liquids have a ``hidden'' approximate scale invariance: Equilibrium potential energy fluctuations are accurately described by inverse power law (IPL) potentials, the radial distribution functions are accurately reproduced by the IPL's, and the radial distribution functions obey the IPL predicted scaling properties to a good approximation. IPL scaling of the dynamics also applies -- with the scaling exponent predicted by the equilibrium fluctuations. In contrast, the equation of state does not obey the IPL scaling. We argue that our results are general for van der Waals liquids, but do not apply, e.g., for hydrogen-bonded liquids., Comment: 7 pages, 9 figures

This paper reports all-atom computer simulations of five phospholipid membranes (DMPC, DPPC, DMPG, DMPS, and DMPSH) with focus on the thermal equilibrium fluctuations of volume, energy, area, thickness, and chain order. At constant temperature and pressure, volume and energy exhibit strong correlations of their slow fluctuations (defined by averaging over 0.5 nanosecond). These quantities, on the other hand, do not correlate significantly with area, thickness, or chain order. The correlations are mainly reported for the fluid phase, but we also give some results for the ordered (gel) phase of two membranes, showing a similar picture. The cause of the observed strong correlations is identified by splitting volume and energy into contributions from tails, heads, and water, and showing that the slow volume-energy fluctuations derive from van der Waals interactions of the tail region; they are thus analogous to the similar strong correlations recently observed in computer simulations of the Lennard-Jones and other simple van der Waals type liquids [U. R. Pedersen et al., Phys. Rev. Lett. 2008, 100, 015701]. The strong correlations reported here confirm one crucial assumption of a recent theory for nerve signal propagation proposed by Heimburg and Jackson [T. Heimburg and A. D. Jackson, Proc. Natl. Acad. Sci. 2005, 102, 9790-9795]., Comment: Accepted by J. Phys. Chem. B