Comparing zero-parameter theories for the WCA and harmonic-repulsive melting lines.

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 are 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.

Comparing four hard-sphere approximations for the low-temperature WCA melting line.

By combining interface-pinning simulations with numerical integration of the Clausius-Clapeyron equation, we accurately determine the melting-line coexistence pressure and fluid/crystal densities of the Weeks-Chandler-Andersen 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 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.

Extreme case of density scaling: The Weeks-Chandler-Andersen system at low temperatures.

This paper studies numerically the Weeks-Chandler-Andersen 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.

Isomorph Invariance of Higher-Order Structural Measures in Four Lennard-Jones Systems.

In the condensed liquid phase, both single- and multicomponent Lennard-Jones (LJ) systems obey the "hidden-scale-invariance" symmetry to a good approximation. Defining an isomorph as a line of constant excess entropy in the thermodynamic phase diagram, the consequent approximate isomorph invariance of structure and dynamics in appropriate units is well documented. However, although all measures of the structure are predicted to be isomorph invariant, with few exceptions only the radial distribution function (RDF) has been investigated. This paper studies the variation along isomorphs of the nearest-neighbor geometry quantified by the occurrence of Voronoi structures, Frank-Kasper bonds, icosahedral local order, and bond-orientational order. Data are presented for the standard LJ system and for three binary LJ mixtures (Kob-Andersen, Wahnström, NiY2). We find that, while the nearest-neighbor geometry generally varies significantly throughout the phase diagram, good invariance is observed along the isomorphs. We conclude that higher-order structural correlations are no less isomorph invariant than is the RDF.

Solid-liquid coexistence of neon, argon, krypton, and xenon studied by simulations.

The noble elements constitute the simplest group of atoms. At low temperatures or high pressures, they freeze into the face-centered cubic (fcc) crystal structure (except helium). This paper investigates neon, argon, krypton, and xenon by molecular dynamics using the simplified atomic potentials recently proposed by Deiters and Sadus [J. Chem. Phys. 150, 134504 (2019)], which are parameterized using data from accurate ab initio quantum-mechanical calculations by the coupled-cluster approach at the single-double-triple level. We compute the fcc freezing lines and find good agreement with the empirical values. At low pressures, predictions are improved by including many-body corrections. Hidden scale invariance of the potential-energy function is established by showing that mean-squared displacement and the static structure factor are invariant along the lines of constant excess entropy (isomorphs). The isomorph theory of melting [Pedersen et al., Nat. Commun. 7, 12386 (2016)] is used to predict from simulations at a single state point the freezing line's shape, the entropy of melting, and the Lindemann parameter of the crystal at melting. Finally, our results suggest that the body-centered cubic crystal is the thermodynamically stable phase at high pressures.

How a supercooled liquid borrows structure from the crystal.

Using computer simulations, we establish that the structure of a supercooled binary atomic liquid mixture consists of common neighbor structures similar to those found in the equilibrium crystal phase, a Laves structure. Despite the large accumulation of the crystal-like structure, we establish that the supercooled liquid represents a true metastable liquid and that liquid can "borrow" the crystal structure without being destabilized. We consider whether this feature might be the origin of all instances of liquids with a strongly favored local structure.

The EXP pair-potential system. IV. Isotherms, isochores, and isomorphs in the two crystalline phases.

This paper studies numerically the solid phase of a system of particles interacting by the exponentially repulsive pair potential, which is a face-centered cubic (fcc) crystal at low densities and a body-centered cubic (bcc) crystal at higher densities [U. R. Pedersen et al., J. Chem. Phys. 150, 174501 (2019)]. Structure is studied via the pair-distribution function and dynamics via the velocity autocorrelation function and the phonon density of states. These quantities are evaluated along isotherms, isochores, and three isomorphs in both crystal phases. Isomorphs are traced out by integrating the density-temperature relation characterizing configurational adiabats, starting from state points in the middle of the fcc-bcc coexistence region. Good isomorph invariance of structure and dynamics is seen in both crystal phases, which is notable in view of the large density variations studied. This is consistent with the fact that the virial potential-energy correlation coefficient is close to unity in the entire fcc phase and in most of the bcc phase (basically below the re-entrant density). Our findings confirm that the isomorph theory, developed and primarily studied for liquids, applies equally well for solids.

The EXP pair-potential system. III. Thermodynamic phase diagram.

This paper determines the thermodynamic phase diagram of the EXP system of particles interacting by the purely repulsive exponential pair potential. The solid phase is face-centered cubic (fcc) at low densities and pressures. At higher densities and pressures, the solid phase is body-centered cubic (bcc) with a re-entrant liquid phase at the highest pressures simulated. The investigation first identifies the phase diagram at zero temperature at which the following four crystal structures are considered: fcc, bcc, hexagonal close packed, and cubic diamond. There is a T = 0 phase transition at pressure 2.651 × 10-3 with the thermodynamically stable structure being fcc below and bcc above this pressure. The densities of the two crystal structures at the phase transition are 1.7469 × 10-2 (fcc) and 1.7471 × 10-2 (bcc). At finite temperatures, the fcc-bcc, fcc-liquid, and bcc-liquid coexistence lines are determined by numerical integration of the Clausius-Clapeyron equation and validated by interface-pinning simulations at selected state points. The bcc-fcc phase transition is a weak first-order transition. The liquid-fcc-bcc triple point, which is determined by the interface-pinning method, has temperature 5.9 × 10-5 and pressure 2.5 × 10-6; the triple-point densities are 1.556 × 10-3 (liquid), 1.583 × 10-3 (bcc), and 1.587 × 10-3 (fcc).

Statistics of small length scale density fluctuations in supercooled viscous liquids.

Many successful theories of liquids near the melting temperature assume that small length scale density fluctuations follow Gaussian statistics. This paper presents a numerical investigation of density fluctuations in the supercooled viscous regime using an enhanced sampling method. Five model systems are investigated: the single component Lennard-Jones liquid, the Kob-Andersen binary mixture, the Wahnström binary mixture, the Lewis-Wahnström model of ortho-terphenyl, and the TIP4P/Ice model of water. The results show that the Gaussian approximation persists to a good degree into the supercooled viscous regime; however, it is less accurate at low temperatures. The analysis suggests that non-Gaussian fluctuations are related to crystalline configurations. Implications for theories of the glass transition are discussed.

Experimental Evidence for a State-Point-Dependent Density-Scaling Exponent of Liquid Dynamics.

A large class of liquids obey density scaling characterized by an exponent, which quantifies the relative roles of temperature and density for the dynamics. We present experimental evidence that the density-scaling exponent γ is state-point dependent for the glass formers tetramethyl-tetraphenyl-trisiloxane (DC704) and 5-polyphenyl ether (5PPE). A method is proposed that from dynamic and thermodynamic properties at equilibrium estimates the value of γ. The method applies at any state point of the pressure-temperature plane, both in the supercooled and the normal liquid regimes. We find that γ is generally state-point dependent, which is confirmed by reanalyzing data for 20 metallic liquids and two model liquids.

Phase Diagram of Kob-Andersen-Type Binary Lennard-Jones Mixtures.

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 a 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 cutoff in a narrow interval of B particle concentrations around 26% at which the bipyramidal orthorhombic PuBr_{3} structure is the thermodynamically stable phase. The melting temperature's variation with B particle concentration at two constant pressures, as well as at the constant density 1.2, is estimated from 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.

Communication: Simple liquids' high-density viscosity.

This paper argues that the viscosity of simple fluids at densities above that of the triple point is a specific function of temperature relative to the freezing temperature at the density in question. The proposed viscosity expression, which is arrived at in part by reference to the isomorph theory of systems with hidden scale invariance, describes computer simulations of the Lennard-Jones system as well as argon and methane experimental data and simulation results for an effective-pair-potential model of liquid sodium.

Thermodynamics of freezing and melting.

Although the freezing of liquids and melting of crystals are fundamental for many areas of the sciences, even simple properties like the temperature-pressure relation along the melting line cannot be predicted today. Here we present a theory in which properties of the coexisting crystal and liquid phases at a single thermodynamic state point provide the basis for calculating the pressure, density and entropy of fusion as functions of temperature along the melting line, as well as the variation along this line of the reduced crystalline vibrational mean-square displacement (the Lindemann ratio), and the liquid's diffusion constant and viscosity. The framework developed, which applies for the sizable class of systems characterized by hidden scale invariance, is validated by computer simulations of the standard 12-6 Lennard-Jones system.

Computing the crystal growth rate by the interface pinning method.

An essential parameter for crystal growth is the kinetic coefficient given by the proportionality between supercooling 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 a 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 and Si from first principles. A generalized version of the method may be used for computing the rates of crystal nucleation or other rare events.

Direct calculation of the solid-liquid Gibbs free energy difference in a single equilibrium simulation.

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.

Factors contributing to the glass-forming ability of a simulated molecular liquid.

The time scales of crystallization of two model liquids, an atomic liquid and a molecular liquid of bent trimers originally introduced as a model of the glass the former o-terphenyl, are determined using molecular dynamics simulations. The molecular liquid is found to have a minimum crystallization time, on supercooling, that is 10(4) times larger than that of the atomic liquid. We present evidence that this enhanced glass-forming ability is due, in equal parts, to the slower dynamics and the larger crystal-liquid interfacial free energy in the molecular liquid.

Pressure-energy correlations in liquids. V. Isomorphs in generalized Lennard-Jones systems.

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. 131, 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 a single curve. Two applications are given. One tests the prediction that the solid-liquid coexistence curve follows an isomorph by comparing to recent simulations by Ahmed and Sadus [J. Chem. Phys. 131, 174504 (2009)]. Excellent agreement is found on the liquid side of the coexistence curve, whereas the agreement is less convincing 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 the potential energy on isochores. It is shown that the new equation of state agrees well with simulations.

Crystallization of the Lewis-Wahnström ortho-terphenyl model.

Crystallization is observed during microsecond long molecular dynamics simulations of bent trimers, a molecular model proposed by Lewis and Wahnströ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 is almost optimal for this structure. The crystal exhibits orientational disorder with the molecules aligned randomly along the three Cartesian axis, i.e., 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 trimer over the monatomic liquid.

Geometry of slow structural fluctuations in a supercooled binary alloy.

The liquid structure of a glass-forming binary alloy is studied using molecular dynamics simulations. The analysis employs the geometrical approach of Frank and Kasper to establish that the supercooled liquid contains extended clusters characterized by the same short range order as the crystal. The steep increase in the heat capacity on cooling is 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.