Search

Understanding the structure and dynamics of liquids is pivotal for the study of larger spatiotemporal processes, especially in glass-forming materials at low temperatures. Density scaling, observed in many molecular systems through experiments, offers an efficient means for exploring a vast range of time scales along a one-dimensional phase diagram. However, the theoretical foundation provided by isomorph theory is of limited use for molecular systems, since currently no first-principles theory exists that can explain the origins of density scaling or make predictions based on it. In this work, we propose a first-principles framework employing coarse-graining in space and time. Spatial coarse-graining reduces a molecule to a center-of-mass-level description by eliminating fast degrees of freedom, while temporal coarse-graining involves averaging fluctuations or correlation functions over characteristic time scales. We show that both approaches enable ab initio estimation of the density scaling coefficient for ortho-terphenyl, consistent with experimental values. Building on these findings, we employ excess entropy scaling to derive a microscopic theory that underpins density scaling from fully atomistic simulations. Our results illuminate the role of coarse-graining in assessing slow fluctuations in molecules and unravel the microscopic nature of density scaling. Ultimately, our proposed framework enables systematic bottom-up approaches for predicting transport coefficients that are otherwise experimentally inaccessible and computationally prohibitive., Comment: 34 pages (26 main, 8 supplemental) and 24 figures (18 main, 6 supplemental)

The melting line of the Weeks-Chandler-Andersen (WCA) system was recently determined accurately and compared to the predictions of four analytical hard-sphere approximations [Attia et al., J. Chem. Phys. 157, 034502 (2022)]. Here, we study an alternative zero-parameter prediction based on the isomorph theory, the input of which relate to properties at a single reference state point on the melting line. The two central assumptions made are that the harmonic-repulsive potential approximates the WCA potential and that pair collisions are uncorrelated. The new approach gives excellent predictions at high temperatures, while the hard-sphere-theory based predictions are better at lower temperatures. Supplementing the WCA investigation, the face-centered-crystal to fluid coexistence line is determined for a system of harmonic-repulsive particles and compared to the zero-parameter theories. The results indicate that the excellent isomorph-theory predictions for the WCA potential at higher temperatures may be partly due to a cancellation of errors between the two above-mentioned assumptions., Comment: 14 pages, 11 figures, data available at doi.org/10.5281/zenodo.7646884

We present a numerical study of the asymmetric dumbbell model consisting of ``molecules'' constructed as two different-sized Lennard-Jones spheres connected by a rigid bond. In terms of the largest (A) particle radius, we report data for the structure and dynamics of the liquid phase for the bond lengths 0.05, 0.1, 0.2, and 0.5, and analogous data for the plastic-crystal phase for the bond lengths 0.05, 0.1, 0.2, and 0.3. Structure is probed by means of the AA, AB, and BB radial distribution functions. Dynamics is probed via the A and B particle mean-square displacement as functions of time and via the rotational time-autocorrelation function. Consistent with the systems' strong virial potential-energy correlations, the structure and dynamics are found to be isomorph invariant to a good approximation in reduced units, while they generally vary considerably along isotherms of the same (20%) density variation. Even the rotational time-autocorrelation function, which due to the constant bond length is not predicted to be isomorph invariant, varies more along isotherms than along isomorphs. Our findings provide the first validation of isomorph-theory predictions for plastic crystals for which isomorph invariance, in fact, is found to apply better than in the liquid phase of asymmetric-dumbbell models.

By combining interface-pinning simulations with numerical integration of the Clausius-Clapeyron equation we determine accurately the melting-line coexistence pressure and fluid/crystal densities of the Weeks-Chandler-Andersen (WCA) system covering four decades of temperature. The data are used for comparing the melting-line predictions of the Boltzmann, Andersen-Weeks-Chandler, Barker-Henderson, and Stillinger hard-sphere approximations. The Andersen-Weeks-Chandler and the Barker-Henderson theories give the most accurate predictions, and they both work excellently in the zero-temperature limit for which analytical expressions are derived here., Comment: 12 pages, 9 figures

This paper studies numerically the Weeks-Chandler-Andersen (WCA) system, which is shown to obey hidden scale invariance with a density-scaling exponent that varies from below 5 to above 500. This unprecedented variation makes it advantageous to use the fourth-order Runge-Kutta algorithm for tracing out isomorphs. Good isomorph invariance of structure and dynamics is observed over more than three orders of magnitude temperature variation. For all state points studied, the virial potential-energy correlation coefficient and the density-scaling exponent are controlled mainly by the temperature. Based on the assumption of statistically independent pair interactions, a mean-field theory is developed that rationalizes this finding and provides an excellent fit to data at low temperatures and densities.

The noble elements constitute the simplest group in the periodic table. At low temperatures or high pressures, the liquid phase solidifies into a face-centered cubic crystal structure (except helium). In the companion paper, we investigated the fcc solid-liquid coexistence of argon in the light of hidden scale invariance. Here we extend the investigation to neon, argon, krypton, and xenon. Computations are done using the SAAP potential, suggested by Deiters and Sadus [J. Chem. Phys 150, 134504 (2019)], derived from accurate ab initio calculations. The systems exhibit hidden scale invariance in the investigated part of the phase diagram, which makes it possible to predict the shape and property variations along the solid-liquid coexistence lines.

Melting temperatures at ambient pressure of systems of isosceles Lennard-Jones trimers with angles ranging from 70 degrees to 100 degrees are determined. Two crystal structures are considered: a distorted body centered cubic structure and a distorted face centered cubic structure with preferred angles of 77 and 96 degrees, respectively. Liquid dynamics are slowed down when the angle is increased. A trimer angle of 83 degrees yields the largest distance between isochrones and the melting temperature, suggesting that this value gives the optimal glass-forming ability. It is conjectured that better glass-formers may be found at angles larger than the ones considered in this study., Comment: 5 pages, 6 figures

Many successful theories of liquids near the melting temperature assume that small length scale density fluctuations follow Gaussian statistics. In this paper I present numerical investigations of fluctuations in the supercooled viscous regime using an enhance sampling method. I present results for the single component Lennard-Jones liquid, the Kob-Andersen binary mixture, the Wahnstr\"om binary mixture, the Lewis-Wahnstr\"om model of ortho-terphenyl and for the TIP4P/Ice model of water. Results show that the Gaussian approximation persist to a good degree into the supercooled viscous regime, however, the approximation is less accurate at low temperatures. I relate the non-Gaussian fluctuations to crystalline configurations. Implications to theories of the glass transition are discussed., Comment: 10 pages, 10 figures

A large class of liquids have hidden scale invariance characterized by a scaling exponent. In this letter we present experimental evidence that the scaling exponent of liquid dynamics is state-point dependent for the glass-forming silicone oil tetramethyl-tetraphenyl-trisiloxane (DC704) and 5-polyphenyl ether (5PPE). From dynamic and thermodynamic properties at equilibrium, we use a method to estimate the value of $\gamma$ at any state point of the pressure-temperature plane, both in the supercooled and normal liquid regimes. We find agreement between the average exponents and the value obtained by superposition of relaxation times over a large range of state-points. We confirm the state-point dependence of $\gamma$ by reanalyzing data of 20 metallic liquids and two model liquids.

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.

We investigate the variation of the driving force for crystallization of a supercooled liquid along isomorphs, curves along which structure and dynamics are invariant. The variation is weak, and can be predicted accurately for the Lennard-Jones fluid using a recently developed formalism and data at a reference temperature. More general analysis allows interpretation of experimental data for molecular liquids such as dimethyl phthalate and indomethacin, and suggests that the isomorph scaling exponent $\gamma$ in these cases is an increasing function of density, although this cannot be seen in measurements of viscosity or relaxation time., Comment: Corresponds to published version, i.e., including changes made during proof

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

Density functional theory (DFT) calculations of 58 liquid elements at their triple point show that most metals exhibit near proportionality between thermal fluctuations between virial and potential-energy in the isochoric ensemble. This demonstrates a general "hidden" scale invariance of metals making the dense part of the thermodynamic phase diagram effectively one dimensional with respect to structure and dynamics. DFT computed density scaling exponents, related to the Gr{\"u}neisen parameter, are in good agreement with experimental values for 16 elements where reliable data were available. Hidden scale invariance is demonstrated in detail for magnesium by showing invariance of structure and dynamics. Computed melting curves of period three metals follow curves with invariance (isomorphs). The experimental structure factor of magnesium is predicted by assuming scale invariant inverse power-law (IPL) pair interactions. However, crystal packings of several transition metals (V, Cr, Mn, Fe, Nb, Mo, Ta, W and Hg), most post-transition metals (Ga, In, Sn, and Tl) and the metalloids Si and Ge cannot be explained by the IPL assumption. Thus, hidden scale invariance can be present even when the IPL-approximation is inadequate. The virial-energy correlation coefficient of iron and phosphorous is shown to increase at elevated pressures. Finally, we discuss how scale invariance explains the Gr{\"u}neisen equation of state and a number of well-known empirical melting and freezing rules., Comment: 12 pages, 11 figures

An essential parameter for crystal growth is the kinetic coefficient given by the proportionality between super-cooling and average growth velocity. Here we show that this coefficient can be computed in a single equilibrium simulation using the interface pinning method where two-phase configurations are stabilized by adding an spring-like bias field coupling to an order-parameter that discriminates between the two phases. Crystal growth is a Smoluchowski process and the crystal growth rate can therefore be computed from the terminal exponential relaxation of the order parameter. The approach is investigated in detail for the Lennard-Jones model. We find that the kinetic coefficient scales as the inverse square-root of temperature along the high temperature part of the melting line. The practical usability of the method is demonstrated by computing the kinetic coefficient of the elements Na, Mg, Al and Si from first principles. It is briefly discussed how a generalized version of the method is an alternative to forward flux sampling methods for computing rates along trajectories of rare events.

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.

Computing phase diagrams of model systems is an essential part of computational condensed matter physics. In this paper we discuss in detail the interface pinning (IP) method for calculation of the Gibbs free energy difference between a solid and a liquid. This is done in a single equilibrium simulation by applying a harmonic field that biases the system towards two-phase configurations. The Gibbs free energy difference between the phases is determined from the average force that the applied field exerts on the system. As a test system we study the Lennard-Jones model. It is shown that the coexistence line can be computed efficiently to a high precision when the IP method is combined with the Newton-Raphson method for finding roots. Statistical and systematic errors are investigated. Advantages and drawbacks of the IP method are discussed. The high pressure part of the temperature-density coexistence region is outlined by isomorphs., Comment: 11 pages, 10 figures

We propose an approach for computing the Gibbs free energy difference between phases of a material. The method is based on the determination of the average force acting on interfaces that separate the two phases of interest. This force, which depends on the Gibbs free energy difference between the phases, is computed by applying an external harmonic field that couples to a parameter which specifies the two phases. Validated first for the Lennard-Jones model, we demonstrate the flexibility, efficiency and practical applicability of this approach by computing the melting temperatures of sodium, magnesium, aluminum and silicon at ambient pressure using density functional theory. Excellent agreement with experiment is found for all four elements, except for silicon, for which the melting temperature is, in agreement with previous simulations, seriously underestimated.

Crystallization is observed during long molecular dynamics simulations of bent trimers, a molecular model proposed by Lewis and Wahnstr\"{o}m for ortho-terphenyl. In the crystal, the three spheres that make up the rigid molecule sit near sites of a body centered cubic lattice, the trimer bond angle being almost optimal for this structure. The crystal exhibits orientational disorder with the molecules aligned randomly along the three Cartesian axis (an example of cubatic orientational order). The rotational and translational mobilities exhibit only modest decreases on crystallization, by factors of 10 and 3 respectively. The rotational relaxation does change from Debye-like in the liquid to large angle jumps in the crystal. We consider the origin of the superior glass forming ability of the the trimer over the monatomic liquid., Comment: 6 pages

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

This paper summarizes the properties of strongly correlating liquids, i.e., liquids with strong correlations between virial and potential energy equilibrium fluctuations at constant volume. We proceed to focus on the experimental predictions for strongly correlating glass-forming liquids. These predictions include i) density scaling, ii) isochronal superposition, iii) that there is a single function from which all frequency-dependent viscoelastic response functions may be calculated, iv) that strongly correlating liquids are approximately single-parameter liquids with close to unity Prigogine-Defay ratio, and v) that the fictive temperature initially decreases for an isobaric temperature up jump. The "isomorph filter", which allows one to test for universality of theories for the non-Arrhenius temperature dependence of the relaxation time, is also briefly discussed.

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

Recently the supercooled Wahnstrom binary Lennard-Jones mixture was partially crystallized into ${\rm MgZn_2}$ phase crystals in lengthy Molecular Dynamics simulations. We present Molecular Dynamics simulations of a modified Kob-Andersen binary Lennard-Jones mixture that also crystallizes in lengthy simulations, here however by forming pure fcc crystals of the majority component. The two findings motivate this paper that gives a general thermodynamic and kinetic treatment of the stability of supercooled binary mixtures, emphasizing the importance of negative mixing enthalpy whenever present. The theory is used to estimate the crystallization time in a Kob-Andersen mixture from the crystallization time in a series of relared systems. At T=0.40 we estimate this time to be 5$\times 10^{7}$ time units ($\approx 1. ms$). A new binary Lennard-Jones mixture is proposed that is not prone to crystallization and faster to simulate than the two standard binary Lennard-Jones mixtures; this is obtained by removing the like-particle attractions by switching to Weeks-Chandler-Andersen type potentials, while maintaining the unlike-particle attraction.

We present an analysis and discuss consequences of the strong correlations of the configurational parts of pressure and energy in their equilibrium fluctuations at fixed volume reported for simulations of several liquids in the companion paper [arXiv:0807.0550]. The analysis concentrates specifically on the single-component Lennard-Jones system. We demonstrate that the potential may be replaced, at fixed volume, by an effective power-law, but not because only short distance encounters dominate the fluctuations. Indeed, contributions to the fluctuations are associated with the whole first peak of the RDF, as we demonstrate by an analysis of the spatially resolved covariance matrix. The reason the effective power-law works so well depends on going beyond single-pair effects and on the constraint of fixed volume. In particular, a better approximation to the potential includes a linear term, which contributes to the mean values of potential energy and virial, but not to their fluctuations. We also study the T=0 limit of the crystalline phase, where the correlation coefficient becomes very close, but not equal, to unity. We then consider four consequences of strong pressure-energy correlations: (1) analyzing experimental data for supercritical Ar we find 96% correlation; (2) we discuss the significance acquired by the correlations for viscous van der Waals liquids approaching the glass transition: For strongly correlating viscous liquids knowledge of just one of the eight frequency-dependent thermoviscoelastic response functions basically implies knowledge of them all; (3) we re-interpret aging simulations of ortho-terphenyl carried out by Mossa {\it et al.} in 2002, showing their conclusions follow from the strongly correlating property; and (4) we discuss correlations in model biomembranes., Comment: Some changes corresponding to those made in the proof of the accepted article

We show that a number of model liquids at fixed volume exhibit strong correlations between equilibrium fluctuations of the configurational parts of (instantaneous) pressure and energy. We present detailed results for thirteen systems, showing in which systems these correlations are significant. These include Lennard-Jones liquids (both single- and two-component) and several other simple liquids, but not hydrogen-bonding liquids like methanol and water, nor the Dzugutov liquid which has significant contributions to pressure at the second nearest neighbor distance. The pressure-energy correlations, which for the Lennard-Jones case are shown to also be present in the crystal and glass phases, reflect an effective inverse power-law potential dominating fluctuations, even at zero and slightly negative pressure. An exception to the inverse-power law explanation is a liquid with hard-sphere repulsion and a square-well attractive part, where a strong correlation is observed, but only after time-averaging. The companion paper [arXiv:0807.0551] gives a thorough analysis of the correlations, with a focus on the Lennard-Jones liquid, and a discussion of some experimental and theoretical consequences., Comment: Some changes corresponding to those made in proof of the accepted article

We address a recent conjecture according to which the relaxation time $\tau$ of a viscous liquid obeys density scaling ($\tau=F(\rho^\gamma/T)$ where $\rho$ is density) if the liquid is ``strongly correlating,'' i.e., has almost 100% correlation between equilibrium virial and potential-energy fluctuations [Pedersen {\it et al.}, PRL {\bf 100}, 011201 (2008)]. Computer simulations of two model liquids - an asymmetric dumbbell model and the Lewis-Wahnstr\"om OTP model - confirm the conjecture and demonstrate that the scaling exponent $\gamma$ can be accurately predicted from equilibrium fluctuations., Comment: 5 pages, 5 figures. Data added for Lewis-Wahnstrom OTP

We first summarize the classical arguments that the vast majority of glass-forming liquids require more than one "order" parameter for their description. Critiques against this conventional wisdom are then presented, and it is argued that the matter deserves to be reconsidered in light of recent experimental developments. Out of the eight basic thermoviscoelastic frequency-dependent response functions, there are generally three independent functions. For stochastic dynamics we show that there are only two independent response functions; for this case it is shown how analytic continuation may be utilized to express the third response functions in terms of two others. Operational criteria are presented for the linear thermoviscoelasticity being described by a single "order" parameter, in which case there is just one independent thermoviscoelastic response function. It is shown that a single "order" parameter description applies to a good approximation whenever thermal equilibrium fluctuations of fundamental variables like energy and pressure are strongly correlated. Results from computer simulations showing that this is the case for a number of simple glass-forming liquids, as well as a few exceptions, are briefly presented. Finally, we discuss a new conjecture according to which experiments at varying temperature and pressure follow the density-scaling expression for the relaxation time, $\tau=F(\rho^x/T)$ ($\rho$ and $T$ are density and temperature), if and only if the liquid is "strongly correlating," i.e., to a good approximation is described by a single "order" parameter.

Constant-pressure molecular-dynamics simulations of phospholipid membranes in the fluid phase reveal strong correlations between equilibrium fluctuations of volume and energy on the nanosecond time-scale. The existence of strong volume-energy correlations was previously deduced indirectly by Heimburg from experiments focusing on the phase transition between the fluid and the ordered gel phases. The correlations, which are reported here for three different membranes (DMPC, DMPS-Na, and DMPSH), have volume-energy correlation coefficients ranging from 0.81 to 0.89. The DMPC membrane was studied at two temperatures showing that the correlation coefficient increases as the phase transition is approached.

We report observation of crystallization of the glass-forming binary Lennard-Jones liquid first used by Wahnstr\"om [G. Wahnstr\"om, Phys. Rev. A 44, 3752 (1991)]. Molecular dynamics simulations of the metastable liquid on a timescale of microseconds were performed. The liquid crystallized spontaneously. The crystal structure was identified as MgZn_2. Formation of transient crystallites is observed in the liquid. The crystallization is investigate at different temperatures and compositions. At high temperature the rate of crystallite formation is the limiting factor, while at low temperature the limiting factor is growth rate. The melting temperature of the crystal is estimated to be T_m=0.93 at rho=0.82. The maximum crystallization rate of the A_2B composition is T=0.60+/-0.02., Comment: 4 pages, 4 figures; corrected typos

Strong correlations between equilibrium fluctuations of the configurational parts of pressure and energy are found in the Lennard-Jones liquid and other simple liquids, but not in hydrogen-bonding liquids like methanol and water. The correlations, that are present also in the crystal and glass phases, reflect an effective inverse power-law repulsive potential dominating fluctuations, even at zero and slightly negative pressure. In experimental data for supercritical Argon, the correlations are found to be approximately 96%. Consequences for viscous liquid dynamics are discussed., Comment: Phys. Rev. Lett., in press

Molecular dynamics results for the dynamic Prigogine-Defay ratio are presented for two glass-forming liquids, thus evaluating the experimentally relevant quantity for testing whether metastable-equilibrium liquid dynamics to a good approximation are described by a single parameter. For the Kob-Andersen binary Lennard-Jones mixture as well as for an asymmetric dumbbell model liquid a single-parameter description works quite well. This is confirmed by time-domain results where it is found that energy and pressure fluctuations are strongly correlated on the alpha-time scale in the NVT ensemble; in the NpT ensemble energy and volume fluctuations similarly correlate strongly., Comment: Phys. Rev. E, in press

Thermo-viscoelastic linear-response functions are calculated from the master equation describing viscous liquid inherent dynamics. From the imaginary parts of the frequency-dependent isobaric specific heat, isothermal compressibility, and isobaric thermal expansion coefficient, we define a "linear dynamic Prigogine-Defay ratio" with the property that if this quantity is unity atone frequency, then it is unity at all frequencies. This happens if and only if there is a single-order-parameter description of the thermo-viscoelastic linear responses via an order parameter (which may be non-exponential in time). Generalizations to other cases of thermodynamic control parameters than temperature and pressure are treated in an Appendix., Comment: Replaces arXiv:cond-mat/0405709

We present a three-dimensional lattice-gas model with trivial thermodynamics, but nontrivial dynamics. The model is characterized by each particle having its own random energy landscape. The equilibrium dynamics of the model were investigated by continuous time Monte Carlo simulations at two different densities at several temperatures. At high densities and low temperatures the model captures the important characteristics of viscous liquid dynamics. We thus observe non-exponential relaxation in the self part of the density auto-correlation function, and fragility plots of the self-diffusion constant and relaxation times show non-Arrhenius behavior., Comment: 6 pages, 6 figures, Submitted for proceedings of the 5th IDMRCS