Evidence of a one-dimensional thermodynamic phase diagram for simple glass-formers

Invariant dynamics in a united-atom model of an ionic liquid

We study a united-atom model of the ionic liquid 1-butyl-1-methylpyrrolidinium bis(trifluoromethyl)sulfonylamide to determine to what extent there exist curves in the phase diagram along which the microscopic dynamics are invariant when expressed in dimensionless, or reduced, form. The initial identification of these curves, termed isodynes, is made by noting that contours of reduced shear viscosity and reduced self-diffusion coefficient coincide to a good approximation. Choosing specifically the contours of reduced viscosity as nominal isodynes, further simulations were carried out for state points on these, and other aspects of dynamics were investigated to study their degree of invariance. These include the mean-squared displacement, shear-stress autocorrelation function, and various rotational correlation functions. These were invariant to a good approximation, with the main exception being rotations of the anion about its long axis. The dynamical features that are invariant have in common that they are aspects that would be relevant for a coarse-grained description of the system; specifically, removing the most microscopic degrees of freedom in principle leads to a simplification of the potential energy landscape, which allows for the existence of isodynes.

Thermobaric history as a tool to govern properties of glasses: case of dipropylene glycol

The elastic properties of high- and low-pressure glasses of dipropylene glycol were determined for the first time under conditions of isothermal compression up to 1 GPa at 77 K and isobaric heating of 77-300 K at 0.05 GPa and 1 GPa. A strong dependence of the elastic properties of glasses on their thermobaric history has been revealed: glasses obtained at high pressure have not only higher densities (3.9%), but also noticeably higher elastic moduli. This effect is especially pronounced in the shear modulus: high-pressure glass has a 30% higher shear modulus than low-pressure glass. The behavior of elastic moduli during the glass-to-liquid transition also depends on the thermobaric history. Glass produced at low pressure but heated at high pressure has anomalous temperature dependences of the elastic moduli. Heating dipropylene glycol glasses at different pressures allowed us to refine the Tg(P) dependence.

Scaling Properties of Liquid Dynamics Predicted from a Single Configuration: Small Rigid Molecules

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. Schrøder, Phys. Rev. Lett. 2022, 129, 245501]. 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öm 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.

Isomorphs in sheared binary Lennard-Jones glass: Transient response

We have studied shear deformation of binary Lennard-Jones glasses to investigate the extent to which the transient part of the stress strain curves is invariant when the thermodynamic state point is varied along an isomorph. Shear deformations were carried out on glass samples of varying stability, determined by cooling rate, and at varying strain rates, at state points deep in the glass. Density changes up to and exceeding a factor of two were made. We investigated several different methods for generating isomorphs but none of the previously developed methods could generate sufficiently precise isomorphs given the large density changes and nonequilibrium situation. Instead, the temperatures for these higher densities were chosen to give state points isomorphic to the starting state point by requiring the steady-state flow stress for isomorphic state points to be invariant in reduced units. In contrast to the steady-state flow stress, we find that the peak stress on the stress strain curve is not invariant. The peak stress decreases by a few percent for each ten percent increase in density, although the differences decrease with increasing density. Analysis of strain profiles and nonaffine motion during the transient phase suggests that the root of the changes in peak stress is a varying tendency to form shear bands, with the largest tendency occurring at the lowest densities. We suggest that this reflects the effective steepness of the potential; a higher effective steepness gives a greater tendency to form shear bands.

Connecting Entropy Scaling and Density Scaling

It is shown that the residual entropy (entropy minus that of the ideal gas at the same temperature and density) is mostly synonymous with the independent variable of density scaling, identifying a direct link between these two approaches. The residual entropy and the effective hardness of interaction (itself a derivative at constant residual entropy) are studied for the Lennard-Jones monomer and dimer as well as a range of rigid molecular models for carbon dioxide. It is observed that the density scaling exponent appears to be related to the two-body interactions in the dilute-gas limit.

A density scaling conjecture for aging glasses

The aging rate of glasses has traditionally been modeled as a function of temperature, T , andfictive temperature, while density, ρ, is not explicitly included as a parameter. However, this de-scription does not naturally connect to the modern understanding of what governs the relaxationrate in equilibrium. In equilibrium it is well known that the relaxation rate, γeq , depends on tem-perature and density. In addition a large class of systems obey density scaling which means therate specifically depends on the scaling parameter, Γ = e(ρ)/T , where e(ρ) is a system specificfunction. This paper present a generalization of the fictive temperature concept in terms of a fic-tive scaling paramter, Γfic , and a density scaling conjecture for aging glasses in which the agingrate depends on Γ and Γfic .

A sample holder for simultaneous neutron and dielectric spectroscopy – dielectric tests with glycerol, glycerol-water, water and phosphoric acid

We report on dielectric test measurements of a rectangular flat sample holder which serves as capacitor and which is aimed for simultaneous neutron and dielectric (n-DE) spectroscopy of acidic liquid samples. We describe technical details of the sample holder assembly and the dielectric and neutron equipment as well as the sample preparation procedure of the air sensitive acidic samples. The sample holder was characterised off-line from the neutron spectrometer by dielectric spectroscopy, but using the standard IN16B cryofurnace with a dielectric sample stick with 4-wire connection and a Novocontrol equipment, previously setup by a collaborative effort between ILL and Roskilde University. Temperature-dependent dielectric scans on standard samples (glycerol, glycerol-water, and Milli-Q water) were measured in the frequency range between 0.27 Hz and 1 MHz. Step-like temperature changes allowed to probe the temperature equilibrium conditions and continuous temperature changes were made to mimic typical IN16B backscattering neutron fixed window scans. Both type of scans were carried out in cooling and in heating. The standard samples show that our dielectric setup with flat sample holder is well suited for simultaneous n-DE-experiments. On the other hand, the dielectric scan on phosphoric acid reveals the limitations of our setup in case of high sample conductivities, but also shows that the DC-conductivity can still be accessed in a sufficiently wide low temperature range where the onset of conductivity can be simultaneously probed with the change in proton dynamics as seen by neutron spectroscopy.

Isochronal superpositioning of the caged dynamics, the α, and the Johari-Goldstein β relaxations in metallic glasses

The superposition of the frequency dispersions of the structural α relaxation determined at different combinations of temperature T and pressure P while maintaining its relaxation time τ_α(T, P) constant (i.e., isochronal superpositioning) has been well established in molecular and polymeric glass-formers. Not known is whether the frequency dispersion or time dependence of the faster processes including the caged molecule dynamics and the Johari-Goldstein (JG) β relaxation possesses the same property. Experimental investigation of this issue is hindered by the lack of an instrument that can cover all three processes. Herein, we report the results from the study of the problem utilizing molecular dynamics simulations of two different glass-forming metallic alloys. The mean square displacement 〈Δr²t〉, the non-Gaussian parameter α₂t, and the self-intermediate scattering function F_sq,t at various combinations of T and P were obtained over broad time range covering the three processes. Isochronal superpositioning of 〈Δr²t〉, α₂t, and F_sq,t was observed over the entire time range, verifying that the property holds not only for the α relaxation but also for the caged dynamics and the JG β relaxation. Moreover, we successfully performed density ρ scaling of the time τ_α2,maxT,P at the peak of α₂t and the diffusion coefficient D(T, P) to show both are functions of ρ^γ/T with the same γ. It follows that the JG β relaxation time τ_β(T, P) is also a function of ρ^γ/T since τ_α2,maxT,P corresponds to τ_β(T, P).

Long Wavelength Thermal Density Fluctuations in Molecular and Polymer Glass-Forming Liquids: Experimental and Theoretical Analysis under Isobaric Conditions

We establish via an in-depth analysis of experimental data that the dimensionless compressibility (proportional to the dimensionless amplitude of long wavelength thermal density fluctuations) of one-component normal and supercooled liquids of chemically complex nonpolar and weakly polar molecules and polymers follows extremely well a surprisingly simple and general temperature dependence over an exceptionally wide range of pressures and temperatures. A theoretical basis for this behavior is shown to exist in the venerable van der Waals model and its more modern interpretations. Although associated hydrogen-bonding (and to a lesser degree strongly polar) liquids display modestly more complex behavior, rather simple temperature and pressure dependences are also discovered. A new approach to collapse the temperature- and pressure-dependent dimensionless compressibility data onto a master curve is formulated that differs from the empirical thermodynamic scaling approach. As a practical matter, we also find that the dimensionless compressibility scales well as an inverse power law with temperature with an exponent that is system dependent and decreases with pressure. At very high pressures and low temperatures, the thermal liquid behavior appears to approach (but not reach) a repulsion-dominated random close packing limit. All these findings are relevant to our recent theoretical work on the problem of activated relaxation and vitrification of supercooled molecular and polymeric liquids.

Theoretical Estimate of the Glass Transition Line of Yukawa One-Component Plasmas

The mode coupling theory of supercooled liquids is combined with advanced closures to the integral equation theory of liquids in order to estimate the glass transition line of Yukawa one-component plasmas from the unscreened Coulomb limit up to the strong screening regime. The present predictions constitute a major improvement over the current literature predictions. The calculations confirm the validity of an existing analytical parameterization of the glass transition line. It is verified that the glass transition line is an approximate isomorphic curve and the value of the corresponding reduced excess entropy is estimated. Capitalizing on the isomorphic nature of the glass transition line, two structural vitrification indicators are identified that allow a rough estimate of the glass transition point only through simple curve metrics of the static properties of supercooled liquids. The vitrification indicators are demonstrated to be quasi-universal by an investigation of hard sphere and inverse power law supercooled liquids. The straightforward extension of the present results to bi-Yukawa systems is also discussed.

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.

Density scaling of structure and dynamics of an ionic liquid

The lines in the pressure–temperature phase diagram with constant conductivity are found to be lines where other dynamic variables as well as the molecular structure factor peak are constant, while charge ordering changes.

Testing density scaling in nanopore-confinement for hydrogen-bonded liquid dipropylene glycol

Recently, it has been demonstrated that the glassy dynamics of the molecular liquids and polymers confined at the nanoscale level might satisfy the density scaling law (ρ ^γ /T) with the same value of the scaling exponent, γ, as that determined from the high-pressure studies of the bulk material. In this work, we have tested the validity of this interesting experimental finding for strongly hydrogen-bonded molecular liquid, dipropylene glycol (DPG), which is known to violate the ρ ^γ /T scaling rule in the supercooled liquid bulk state. The results of the independent dielectric relaxation studies carried out on increased pressure and in nanopores, have led to an important finding that when the density change induced by geometrical confinement is not very large, DPG can still obey the density scaling law with the same value of the scaling exponent as that found for the bulk sample. In this way, we confirm that the information obtained from the universal density scaling approach applied to nanoscale confined systems is somehow consistent with the macroscopic ones and that in both cases the same fundamental rules governs the glass-transition dynamics.

Perspective: Searching for simplicity rather than universality in glass-forming liquids

This article gives an overview of experimental results on dynamics in bulk glass-forming molecular liquids. Rather than looking for phenomenology that is universal, in the sense that it is seen in all liquids, the focus is on identifying the basic characteristics, or "stylized facts," of the glass transition problem, i.e., the central observations that a theory of the physics of glass formation should aim to explain in a unified manner.

Perspective: Excess-entropy scaling

This article gives an overview of excess-entropy scaling, the 1977 discovery by Rosenfeld that entropy determines properties of liquids like viscosity, diffusion constant, and heat conductivity. We give examples from computer simulations confirming this intriguing connection between dynamics and thermodynamics, counterexamples, and experimental validations. Recent uses in application-related contexts are reviewed, and theories proposed for the origin of excess-entropy scaling are briefly summarized. It is shown that if two thermodynamic state points of a liquid have the same microscopic dynamics, they must have the same excess entropy. In this case, the potential-energy function exhibits a symmetry termed hidden scale invariance, stating that the ordering of the potential energies of configurations is maintained if these are scaled uniformly to a different density. This property leads to the isomorph theory, which provides a general framework for excess-entropy scaling and illuminates, in particular, why this does not apply rigorously and universally. It remains an open question whether all aspects of excess-entropy scaling and related regularities reflect hidden scale invariance in one form or other.

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.

Isomorph theory of physical aging

This paper derives and discusses the configuration-space Langevin equation describing a physically aging R-simple system and the corresponding Smoluchowski equation. Externally controlled thermodynamic variables like temperature, density, and pressure enter the description via the single parameter T_s/T, in which T is the bath temperature and T_s is the "systemic" temperature defined at any time t as the thermodynamic equilibrium temperature of the state point with density ρ(t) and potential energy U(t). In equilibrium, T_s ≅ T with fluctuations that vanish in the thermodynamic limit. In contrast to Tool's fictive temperature and other effective temperatures in glass science, the systemic temperature is defined for any configuration with a well-defined density, even if it is not close to equilibrium. Density and systemic temperature define an aging phase diagram, in which the aging system traces out a curve. Predictions are discussed for aging following various density-temperature and pressure-temperature jumps from one equilibrium state to another, as well as for a few other scenarios. The proposed theory implies that R-simple glass-forming liquids are characterized by the dynamic Prigogine-Defay ratio being equal to unity.

High-pressure cell for simultaneous dielectric and neutron spectroscopy

In this article, we report on the design, manufacture, and testing of a high-pressure cell for simultaneous dielectric and neutron spectroscopy. This cell is a unique tool for studying dynamics on different time scales, from kilo- to picoseconds, covering universal features such as the α relaxation and fast vibrations at the same time. The cell, constructed in cylindrical geometry, is made of a high-strength aluminum alloy and operates up to 500 MPa in a temperature range between roughly 2 and 320 K. In order to measure the scattered neutron intensity and the sample capacitance simultaneously, a cylindrical capacitor is positioned within the bore of the high-pressure container. The capacitor consists of two concentric electrodes separated by insulating spacers. The performance of this setup has been successfully verified by collecting simultaneous dielectric and neutron spectroscopy data on dipropylene glycol, using both backscattering and time-of-flight instruments. We have carried out the experiments at different combinations of temperature and pressure in both the supercooled liquid and glassy state.