Isomorphs, hidden scale invariance, and quasiuniversality

NVU view on energy polydisperse Lennard-Jones systems

When energy polydispersity is introduced into the Lennard-Jones (LJ) system, there is little effect on structure and dynamics [T. S. Ingebrigtsen and J. C. Dyre, J. Phys. Chem. B 127, 2837 (2023)10.1021/acs.jpcb.3c00346]. For instance, at a given state point both the radial distribution function and the mean-square displacement as a function of time are virtually unaffected by even large energy polydispersity, which is in stark contrast to what happens when size polydispersity is introduced. We here argue-and validate by simulations of up to 30% polydispersity-that this almost invariance of structure and dynamics reflects an approximate invariance of the constant-potential-energy surface. Because NVU dynamics defined as geodesic motion at constant potential energy is equivalent to Newtonian dynamics in the thermodynamic limit, the approximate invariance of the constant-potential-energy surface implies virtually the same physics of energy polydisperse LJ systems as of the standard single-component version. In contrast, the constant-potential-energy surface is significantly affected by introducing size polydispersity.

Entropy of strongly coupled Yukawa fluids

The entropy of strongly coupled Yukawa fluids is discussed from several perspectives. First, it is demonstrated that a vibrational paradigm of atomic dynamics in dense fluids can be used to obtain a simple and accurate estimate of the entropy without any adjustable parameters. Second, it is explained why a quasiuniversal value of the excess entropy of simple fluids at the freezing point should be expected, and it is demonstrated that a remaining very weak dependence of the freezing point entropy on the screening parameter in the Yukawa fluid can be described by a simple linear function. Third, a scaling of the excess entropy with the freezing temperature is examined, a modified form of the Rosenfeld-Tarazona scaling is put forward, and some consequences are briefly discussed. Fourth, the location of the Frenkel line on the phase diagram of Yukawa systems is discussed in terms of the excess entropy and compared with some predictions made in the literature. Fifth, the excess entropy scaling of the transport coefficients (self-diffusion, viscosity, and thermal conductivity) is reexamined using the contemporary datasets for the transport properties of Yukawa fluids. The results could be of particular interest in the context of complex (dusty) plasmas, colloidal suspensions, electrolytes, and other related systems with soft pairwise interactions.

Melting curve of two-dimensional Yukawa systems predicted by isomorph theory

The analytical expression for the conditions of the solid-fluid phase transition, i.e., the melting curve, for two-dimensional (2D) Yukawa systems is derived theoretically from the isomorph theory. To demonstrate that the isomorph theory is applicable to 2D Yukawa systems, molecular dynamical simulations are performed under various conditions. Based on the isomorph theory, the analytical isomorphic curves of 2D Yukawa systems are derived using the local effective power-law exponent of the Yukawa potential. From the obtained analytical isomorphic curves, the melting curve of 2D Yukawa systems is directly determined using only two known melting points. The determined melting curve of 2D Yukawa systems well agrees with the previous obtained melting results using completely different approaches.

Universal scaling of transverse sound speed and its isomorphic property in Yukawa fluids

Molecular dynamical simulations are performed to investigate the scaling of the transverse sound speed in two-dimensional (2D) and 3D Yukawa fluids. From the calculated diagnostics of the radial distribution function, the mean-squared displacement, and the Pearson correlation coefficient, the approximate isomorphic curves for 2D and 3D liquidlike Yukawa systems are obtained. It is found that the structure and dynamics of 2D and 3D liquidlike Yukawa systems exhibit the isomorphic property under the conditions of the same relative coupling parameter Γ/Γ_{m}=const. It is demonstrated that the reduced transverse sound speed, i.e., the ratio of the transverse sound speed to the thermal speed, is an isomorph invariant, which is a quasiuniversal function of Γ/Γ_{m}. The obtained isomorph invariant of the reduced transverse sound speed can be useful to estimate the transverse sound speed, or determine the coupling strength, with applications to dusty (complex) plasma or colloidal systems.

Bridge functions of classical one-component plasmas

In a recent paper, Lucco Castello et al. [arXiv:2107.03537] performed systematic extractions of classical one-component plasma bridge functions from molecular dynamics simulations and provided an accurate parametrization that was incorporated in their isomorph-based empirically modified hypernetted chain approach for Yukawa one-component plasmas. Here the extraction technique and parametrization strategy are described in detail, while the deficiencies of earlier efforts are discussed. The structural and thermodynamic predictions of the updated version of the integral equation theory approach are compared with extensive available simulation results revealing a truly unprecedented level of accuracy in the entire dense liquid region of the Yukawa phase diagram.

Generalized hydrodynamics of the Lennard-Jones liquid in view of hidden scale invariance

In recent years lines along which structure and dynamics are invariant to a good approximation, so-called isomorphs, have been identified in the thermodynamic phase diagrams of several model liquids and solids. This paper reports computer simulation data of the transverse and longitudinal collective dynamics at different length scales along an isomorph of the Lennard-Jones system. Our findings are compared to corresponding results along an isotherm and an isochore. Confirming the theoretical prediction, the reduced-unit dynamics of the transverse momentum density is invariant to a good approximation along the isomorph on all time and length scales. Likewise, the wave-vector dependent shear-stress autocorrelation function is found to be isomorph invariant (with minor deviations at very short times). A similar invariance is not seen along the isotherm or the isochore. Using a spatially nonlocal hydrodynamic model for the transverse momentum-density time-autocorrelation function, the macroscopic shear viscosity and its wave dependence are determined, demonstrating that the shear viscosity is isomorphic invariant on all length scales studied. This analysis implies the existence of a length scale that is isomorph invariant in reduced units, i.e., which characterizes each isomorph. The transverse sound-wave velocity, the Maxwell relaxation time, and the rigidity shear modulus are also isomorph invariant. In contrast to the isomorph invariance of all aspects of the transverse dynamics, the reduced-unit dynamics of the mass density is not invariant on length scales longer than the interparticle distance. By fitting to a generalized hydrodynamic model, we extract values for the wave-vector-dependent thermal diffusion coefficient, sound attenuation coefficient, and adiabatic sound velocity. The isomorph variation of these quantities in reduced units on long length scales can be eliminated by scaling with the density-scaling exponent, a fundamental quantity in the isomorph theory framework; this is an empirical observation that remains to be explained theoretically.

Excess entropy scaling law: A potential energy landscape view

The relationship between excess entropy and diffusion is revisited by means of large-scale computer simulation combined to supervised learning approach to determine the excess entropy for the Lennard-Jones potential. Results reveal a strong correlation with the properties of the potential energy landscape (PEL). In particular the exponential law holding in the liquid is seen to be linked with the landscape-influenced regime of the PEL whereas the fluidlike power-law corresponds to the free diffusion regime.

Effect of attractions on the hard-sphere dynamic universality class

The fundamental understanding of the dynamic and transport properties of liquids is crucial for the better processing of most materials. The usefulness of this understanding increases when it involves general scaling rules, such as the concept of the hard-sphere dynamic universality class, which provides a unifying scaling of the dynamics of soft-sphere repulsive systems. A relevant question is how far this concept extends to systems that also involve attractive interactions. To answer this question, in this work we performed systematic molecular and Brownian dynamics simulations with the Lennard-Jones system in a wide range of temperatures and densities and verify the extent to which its static and dynamic properties map onto those of the hard-sphere system. We determine that in most of the fluid regime, the Lennard-Jones liquid exhibits the same dynamic equivalence with the hard-sphere system as most purely repulsive fluids, thus establishing the degree of its inclusion in the hard-sphere dynamic universality class.

Understanding simple liquids through statistical and deep learning approaches

Statistical and deep learning-based methods are employed to obtain insights into the quasi-universal properties of simple liquids. In the first part, a statistical model is employed to provide a probabilistic explanation for the similarity in the structure of simple liquids interacting with different pair potential forms, collectively known as simple liquids. The methodology works by sampling the radial distribution function and the number of interacting particles within the cutoff distance, and it produces the probability density function of the net force. We show that matching the probability distribution of the net force can be a direct route to parameterize simple liquid pair potentials with a similar structure, as the net force is the main component of the Newtonian equations of motion. The statistical model is assessed and validated against various cases. In the second part, we exploit DeepILST [A. Moradzadeh and N. R. Aluru, J. Phys. Chem. Lett. 10, 1242-1250 (2019)], a data-driven and deep-learning assisted framework to parameterize the standard 12-6 Lennard-Jones (LJ) pair potential, to find structurally equivalent/isomorphic LJ liquids that identify constant order parameter [τ=∫₀ ^ξ_cf gξ-1ξ²dξ, where gξ and ξ(=rρ¹³) are the reduced radial distribution function and radial distance, respectively] systems in the space of non-dimensional temperature and density of the LJ liquids. We also investigate the consistency of DeepILST in reproducibility of radial distribution functions of various quasi-universal potentials, e.g., exponential, inverse-power-law, and Yukawa pair potentials, quantified based on the radial distribution functions and Kullback-Leibler errors. Our results provide insights into the quasi-universality of simple liquids using the statistical and deep learning methods.

Testing the isomorph invariance of the bridge functions of Yukawa one-component plasmas

It has been recently conjectured that bridge functions remain nearly invariant along phase diagram lines of constant excess entropy for the broad class of R-simple liquids. To test this hypothesis, the bridge functions of Yukawa systems are computed outside the correlation void with the Ornstein-Zernike inversion method employing structural input from ultra-accurate molecular dynamics simulations and inside the correlation void with the cavity distribution method employing structural input from ultra-long specially designed molecular dynamics simulations featuring a tagged particle pair. Yukawa bridge functions are revealed to be isomorph invariant to a very high degree. The observed invariance is not exact, however, since isomorphic deviations exceed the overall uncertainties.

Isomorph theory beyond thermal equilibrium

This paper generalizes isomorph theory to systems that are not in thermal equilibrium. The systems are assumed to be R-simple, i.e., to have a potential energy that as a function of all particle coordinates R obeys the hidden-scale-invariance condition U(R_a) < U(R_b) ⇒ U(λR_a) < U(λR_b). "Systemic isomorphs" are introduced as lines of constant excess entropy in the phase diagram defined by density and systemic temperature, which is the temperature of the equilibrium state point with the average potential energy equal to U(R). The dynamics is invariant along a systemic isomorph if there is a constant ratio between the systemic and the bath temperature. In thermal equilibrium, the systemic temperature is equal to the bath temperature and the original isomorph formalism is recovered. The new approach rationalizes within a consistent framework previously published observations of isomorph invariance in simulations involving nonlinear steady-state shear flows, zero-temperature plastic flows, and glass-state isomorphs. This paper relates briefly to granular media, physical aging, and active matter. Finally, we discuss the possibility that the energy unit defining the reduced quantities should be based on the systemic rather than the bath temperature.

Advancing predictions of protein stability in the solid state

The β-relaxation associated with the sub-glass transition temperature (T_g,β) is attributed to fast, localised molecular motions which can occur below the primary glass transition temperature (T_g,α). Consistent with T_g,β being observed well-below storage temperatures, the β-relaxation associated motions have been hypothesised to influence protein stability in the solid state and could thus impact the quality of e.g. protein powders for inhalation or reconstitution and injection. Why then do distinct solid state protein formulations with similar aggregation profiles after drying and immediate reconstitution, display different profiles when reconstituted following prolonged storage? Is the value of T_g,β, associated with the β-relaxation process of the system, a reliable parameter for characterising the behaviour of proteins in the solid state? Bearing this in mind, in this work we further explore the different relaxation dynamics of glassy solid state monoclonal antibody formulations using terahertz time-domain spectroscopy and dynamical mechanical analysis. By conducting a 52 week stability study on a series of multi-component spray-dried formulations, an approach for characterising and analysing the solid state dynamics and how these relate to protein stability is outlined.

Effective hardness of interaction from thermodynamics and viscosity in dilute gases

The hardness of the effective inverse power law (IPL) potential, which can be obtained from thermodynamics or collision integrals, can be used to demonstrate similarities between thermodynamic and transport properties. This link is investigated for systems of increasing complexity (i.e., the EXP, square-well, Lennard-Jones, and Stockmayer potentials; ab initio results for small molecules; and rigid linear chains of Lennard-Jones sites). These results show that while the two approaches do not yield precisely the same values of effective IPL exponent, their qualitative behavior is intriguingly similar, offering a new way of understanding the effective interactions between molecules, especially at high temperatures. In both approaches, the effective hardness is obtained from a double-logarithmic temperature derivative of an effective area.

Role of Attractive Forces in the Relaxation Dynamics of Supercooled Liquids

The attractive tail of the intermolecular interaction affects very weakly the structural properties of liquids, while it affects dramatically their dynamical ones. Via the numerical simulations of model systems not prone to crystallization, both in three and in two spatial dimensions, here we demonstrate that the nonperturbative dynamical effects of the attractive forces are tantamount to a rescaling of the activation energy by the glass transition temperature T_{g}: systems only differing in their attractive interaction have the same structural and dynamical properties if compared at the same value of T/T_{g}.

Entropy Scaling of Viscosity - I: A Case Study of Propane

In this work, a broadly-applicable and simple approach for building high accuracy viscosity correlations is demonstrated for propane. The approach is based on the combination of a number of recent insights related to the use of residual entropy scaling, especially a new way of scaling the viscosity for consistency with the dilute-gas limit. With three adjustable parameters in the dense phase, the primary viscosity data for propane are predicted with a mean absolute relative deviation of 1.38%, and 95% of the primary data are predicted within a relative error band of less than 5%. The dimensionality of the dense-phase contribution is reduced from the conventional two dimensional approach (temperature and density) to a one-dimensional correlation with residual entropy as the independent variable. The simplicity of the model formulation ensures smooth extrapolation behavior (barring errors in the equation of state itself). The approach proposed here should be applicable to a wide range of chemical species. The supporting information includes the relevant data in tabular form and a Python implementation of the model.

Thermodynamics, dynamics, and structure of supercritical water at extreme conditions

Molecular dynamics (MD) simulations to understand the thermodynamic, dynamic, and structural changes in supercritical water across the Frenkel line and the melting line have been performed.

Modified Entropy Scaling of the Transport Properties of the Lennard-Jones Fluid

Rosenfeld proposed two different scaling approaches to model the transport properties of fluids, separated by 22 years, one valid in the dilute gas, and another in the liquid phase. In this work, we demonstrate that these two limiting cases can be connected through the use of a novel approach to scaling transport properties and a bridging function. This approach, which is empirical and not derived from theory, is used to generate reference correlations for the transport properties of the Lennard-Jones 12-6 fluid of viscosity, thermal conductivity, and self-diffusion. This approach, with a very simple functional form, allows for the reproduction of the most accurate simulation data to within nearly their statistical uncertainty. The correlations are used to confirm that for the Lennard-Jones fluid the appropriately scaled transport properties are nearly monovariate functions of the excess entropy from low-density gases into the supercooled phase and up to extreme temperatures. This study represents the most comprehensive metastudy of the transport properties of the Lennard-Jones fluid to date.

Topological extension of the isomorph theory based on the Shannon entropy

Isomorph theory is one of the promising theories for understanding the quasiuniversal relationship between thermodynamic, dynamic, and structural characteristics. Based on the hidden scale invariance of the inverse power law potentials, it rationalizes the excess entropy scaling law of dynamic properties. This work aims to show that this basic idea of isomorph theory can be extended by examining the microstructural features of the system. Using the topological framework in conjunction with the entropy calculation algorithm, we demonstrate that Voronoi entropy, a measure of the topological diversity of single atoms, provides a scaling law for the transport properties of soft-sphere fluids, which is comparable to the frequently used excess entropy scaling. By examining the relationship between the Voronoi entropy and the solidlike fraction of simple fluids, we suggest that the Frenkel line, a rigid-nonrigid crossover line, be a topological isomorphic line at which the scaling relation qualitatively changes.

Probing the link between residual entropy and viscosity of molecular fluids and model potentials

This work investigates the link between residual entropy and viscosity based on wide-ranging, highly accurate experimental and simulation data. This link was originally postulated by Rosenfeld in 1977 [Rosenfeld Y (1977) Phys Rev A 15:2545-2549], and it is shown that this scaling results in an approximately monovariate relationship between residual entropy and reduced viscosity for a wide range of molecular fluids [argon, methane, [Formula: see text], [Formula: see text], refrigerant R-134a (1,1,1,2-tetrafluoroethane), refrigerant R-125 (pentafluoroethane), methanol, and water] and a range of model potentials (hard sphere, inverse power, Lennard-Jones, and Weeks-Chandler-Andersen). While the proposed "universal" correlation of Rosenfeld is shown to be far from universal, when used with the appropriate density scaling for molecular fluids, the viscosity of nonassociating molecular fluids can be mapped onto the model potentials. This mapping results in a length scale that is proportional to the cube root of experimentally measurable liquid volume values.

Determination of derived volumetric properties and heat capacities at high pressures using two density scaling based equations of state. Application to dipentaerythritol hexa(3,5,5-trimethylhexanoate)

Scaling based EoSs describe the complex behavior of derived properties for broad temperature and pressure ranges from diPEiC₉ experimental densities.

Reference Correlation for the Viscosity of Carbon Dioxide

A comprehensive database of experimental and computed data for the viscosity of carbon dioxide (CO2) was compiled and a new reference correlation was developed. Literature results based on an ab initio potential energy surface were the foundation of the correlation of the viscosity in the limit of zero density in the temperature range from 100 K to 2000 K. Guided symbolic regression was employed to obtain a new functional form that extrapolates correctly to T → 0 K and to 10 000 K. Coordinated measurements at low density made it possible to implement the temperature dependence of the Rainwater-Friend theory in the linear-in-density viscosity term. The residual viscosity could be formulated with a scaling term ργ /T the significance of which was confirmed by symbolic regression. The final viscosity correlation covers temperatures from 100 K to 2000 K for gaseous CO2, and from 220 K to 700 K with pressures along the melting line up to 8000 MPa for compressed and supercritical liquid states. The data representation is more accurate than with the previous correlations, and the covered pressure and temperature range is significantly extended. The critical enhancement of the viscosity of CO2 is included in the new correlation.

Simple liquids' quasiuniversality and the hard-sphere paradigm

This topical review discusses the quasiuniversality of simple liquids' structure and dynamics and two possible justifications of it. The traditional one is based on the van der Waals picture of liquids in which the hard-sphere system reflects the basic physics. An alternative explanation argues that all quasiuniversal liquids to a good approximation conform to the same equation of motion, referring to the exponentially repulsive pair-potential system as the basic reference system. The paper, which is aimed at non-experts, ends by listing a number of open problems in the field.

Freezing and melting line invariants of the Lennard-Jones system

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 isomorph theory. First the freezing/melting lines of the 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 phases. The structure is probed by the radial distribution function and the structure factor and dynamics are probed by the mean-square displacement, the intermediate scattering function, and the shear viscosity. Studying these properties with reference to isomorph theory explains why the known single-phase melting criteria hold, e.g., the Hansen-Verlet and the Lindemann criteria, 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 isomorph theory and that the invariants are not peculiar to the freezing and melting lines, but hold along all isomorphs.

Simplicity of condensed matter at its core: generic definition of a Roskilde-simple system

The isomorph theory is reformulated by defining Roskilde-simple systems 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. If the potential energy as a function of all particle coordinates is denoted by U(R), this requirement translates into U(Ra) < U(Rb) ⇒ U(λRa) < U(λRb). 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.

Density scaling and quasiuniversality of flow-event statistics for athermal plastic flows

Athermal steady-state plastic flows were simulated for the Kob-Andersen binary Lennard-Jones system and its repulsive version in which the sign of the attractive terms is changed to a plus. Properties evaluated include the distributions of energy drops, stress drops, and strain intervals between the flow events. We show that simulations at a single density in conjunction with an equilibrium-liquid simulation at the same density allow one to predict the plastic flow-event statistics at other densities. This is done by applying the recently established "hidden scale invariance" of simple liquids to the glass phase. The resulting scaling of flow-event properties reveals quasiuniversality, i.e., that the probability distributions of energy drops, stress drops, and strain intervals in properly reduced units are virtually independent of the microscopic pair potentials.

Variation of the dynamic susceptibility along an isochrone

Koperwas et al. showed in a recent paper [Phys. Rev. Lett. 111, 125701 (2013)] that the dynamic susceptibility χ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 noninvariance of the isothermal bulk modulus. Considering the version of χ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.

Hidden scale invariance in condensed matter

Recent developments show that many liquids and solids have an approximate "hidden" scale invariance that implies the existence of lines in the thermodynamic phase diagram, so-called isomorphs, along which structure and dynamics in properly reduced units are invariant to a good approximation. This means that the phase diagram becomes effectively one-dimensional with regard to several physical properties. Liquids and solids with isomorphs include most or all van der Waals bonded systems and metals, as well as weakly ionic or dipolar systems. On the other hand, systems with directional bonding (hydrogen bonds or covalent bonds) or strong Coulomb forces generally do not exhibit hidden scale invariance. The article reviews the theory behind this picture of condensed matter and the evidence for it coming from computer simulations and experiments.

Investigating isomorphs with the topological cluster classification

Isomorphs are lines in the density-temperature plane of certain "strongly correlating" or "Roskilde simple" liquids where two-point structure and dynamics have been shown to be close to identical up to a scale transformation. Here we consider such a liquid, a Lennard-Jones glass former, and investigate the behavior along isomorphs of higher-order structural and dynamical correlations. We then consider an inverse power law reference system mapped to the Lennard-Jones system [Pedersen et al., Phys. Rev. Lett. 105, 157801 (2010)]. Using the topological cluster classification to identify higher-order structures, in both systems we find bicapped square antiprisms, which are known to be a locally favored structure in the Lennard-Jones glass former. The population of these locally favored structures is up to 80% higher in the Lennard-Jones system than the equivalent inverse power law system. The structural relaxation time of the two systems, on the other hand, is almost identical, and the four-point dynamical susceptibility is marginally higher in the inverse power law system. Upon cooling, the lifetime of the locally favored structures in the Lennard-Jones system is up to 40% higher relative to the reference system.

Universal features of dynamic heterogeneity in supercooled liquids

A few years ago it was shown that some systems that have very similar local structure, as quantified by the pair correlation function, exhibit vastly different slowing down upon supercooling. Recently, a more subtle structural quantity, the so-called "point-to-set" length, was found to reliably correlate with the average dynamics. Here we use computer simulations to examine the behavior of fluctuations around the average dynamics, i.e., dynamic heterogeneity. We study five model glass-forming liquids: three model liquids used in previous works and two additional model liquids with finite range interactions. Some of these systems have very similar local structure but vastly different dynamics. We show that for all these systems the spatial extent and the anisotropy of dynamic heterogeneity have the same correlation with the average dynamics. This result complements a recent experimental finding of a universal correlation between the number of correlated particles and the apparent activation enthalpy.