Isochronal temperature–pressure superpositioning of the α-relaxation in type-A glass formers

Scaling Properties of Liquid Dynamics Predicted from a Single Configuration: Pseudoisomorphs for Harmonic-Bonded Molecules

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

A study of OH···O hydrogen bonds along various isolines in 2-ethyl-1-hexanol. Temperature or pressure - which parameter controls their behavior?

The nature of H-bonding interactions is still far from being understood despite intense experimental and theoretical studies on this subject carried out by the leading research centers. In this paper, by a combination of unique high-pressure infrared, dielectric and volumetric data, the intramolecular dynamics of hydroxyl moieties (which provides direct information about H-bonds) was studied along various isolines, i.e., isotherms, isobars, isochrones, and isochores, in a simple monohydroxy alcohol (2-ethyl-1-hexanol). This allowed us to discover that the temperature controls the intermolecular hydrogen bonds, which then affect the intramolecular dynamics of OH units. Although the role of density fluctuations gets stronger as temperature rises. We also demonstrated a clear connection between the intra- and intermolecular dynamics of the associating liquid at high pressure. The data reported herein open a new perspective to explore this important aspect of the glass transition phenomenon and understand H-bonding interactions at varying thermodynamic conditions.

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 .

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.

Dynamic Properties of Supercooled Chlorinated Biphenyls

A study of the dynamics of a series of biphenyl compounds having varying chlorine levels was carried out. Increasing the chlorine content increases the glass transition temperature and makes the dynamics substantially more sensitive to density changes. Nonetheless, in the vicinity of their respective glass transitions, the different liquids display very similar extents of dynamic correlation and dynamic heterogeneity. The slight narrowing of the relaxation peak with increasing chlorine follows the general trend of the effect of increasing molecular polarity. This relationship between the peak breadth and dipole moment was reproduced in molecular dynamics simulations of a simplified model of the Aroclor molecule.

Studies on the internal medium-range ordering and high pressure dynamics in modified ibuprofens

Broadband dielectric spectroscopy (BDS), combined with the X-ray diffraction (XRD) and Fourier transform infrared (FTIR) techniques, was used to study the dynamics of the primary (α) relaxation process and slow mode (SM), as well as structural properties and intermolecular interactions, in the methyl-, isopropyl-, hexyl-, and benzyl derivative of a well-known pharmaceutical, ibuprofen (IBU). Unexpectedly, the XRD and FTIR methods revealed the formation of medium-range ordering together with some molecular organization, which probably leads to the creation of small aggregates at the scale of several microns at lower temperatures. Moreover, high pressure dielectric experiments revealed that the SM (observed in the ambient pressure data) is not detected in the loss spectra of compressed IBU esters, which is consistent with the results reported previously for propylene carbonate and dioxolane derivatives. This finding can be interpreted as connected to either the comparable time scale of the structural dynamics and slow mode or suppression of the motions responsible for the latter process at elevated pressure. Additionally, it was found that the pressure coefficient of the glass transition temperature (dTg/dp) and activation volume (ΔV) change with molecular weight (Mw) in a non-monotonic way. It might be related to various chemical structures, conformations, and intermolecular interactions, as well as different architecture of supramolecular aggregates in the investigated compounds.

Chain Flexibility and the Segmental Dynamics of Polymers

Using molecular dynamics simulations, we examine the dynamics of a family of model polymers with varying chain length and torsional potential barriers. We focus on features of the dynamics of polymers that are seen experimentally but absent in simulations of freely rotating and freely jointed chains. The reduced effect of volume on the segmental dynamics with increasing chain length, a capacity for pressure densification, and the deviation from constant Johari-Goldstein relaxation time at a constant segmental relaxation time all have a common origin, torsional rigidity, and these effects become increasingly apparent for more rigid chains.

Elastic properties of liquid and glassy propane-based alcohols under high pressure: the increasing role of hydrogen bonds in a homologous family

We have measured the elastic moduli of liquid and glassy n-propanol and propylene glycol (PG) under pressure by ultrasonic techniques and have recalculated similar characteristics for glycerol from the previous experiment. All three substances form a ternary homologous family with the common formula C3H8-n(OH)n (n = 1, 2, 3), where the number of hydrogen bonds per molecule increases with the number of oxygen atoms approximately as ≈2n. In turn, the enhancement of hydrogen bonding results in an increase in elastic moduli (bulk modulus for liquids or bulk and shear moduli for glasses) from n-propanol to glycerol at all pressures, while the volume per molecule Vm shows the opposite trend at atmospheric pressure in spite of an increase in the molecular size. Nevertheless, the ratios between the Vm values at pressure P > 0.05 GPa are inverted in liquids and tend to the ratios of molecule volumes which indicates a decrease of the relative contribution of hydrogen bonds to the repulsive intermolecular forces with increasing pressure regardless of increase or decrease in the number of hydrogen bonds and their strength. A similar volume behavior is observed for glasses at T = 77 K. We have also established that the relative difference between corresponding moduli of liquid or glassy n-propanol and PG is remarkably less than that between corresponding values for PG and glycerol. We explain this property by the formation of a three-dimensional network of hydrogen bonds in glycerol, where the number of hydrogen bonds per molecule is close to six.

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.

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 dielectric studies on 1,6-anhydro-β-D-mannopyranose (plastic crystal) and 2,3,4-tri-O-acetyl-1,6-anhydro-β-D-glucopyranose (canonical glass)

Broadband Dielectric Spectroscopy was applied to investigate molecular dynamics of two anhydrosaccharides, i.e., 1,6-anhydro-β-D-mannopyranose, anhMAN (hydrogen-bonded system) and 2,3,4-tri-O-acetyl-1,6-anhydro-β-D-glucopyranose, ac-anhGLU (van der Waals material), at different thermodynamic conditions. Moreover, the reported data were compared with those recently published for two other H-bonded systems, i.e., 1,6-anhydro-β-D-glucopyranose (anhGLU) and D-glucose (D-GLU). A direct comparison of the dynamical behavior of the materials with a similar chemical structure but significantly differing by the degrees of freedom, complexity, and intermolecular interactions made it possible to probe the impact of compression on the fragility, Temperature-Pressure Superpositioning and pressure coefficient of the glassy crystal/glass transition temperatures (dT_gc/dp ; dT_g/dp). Moreover, the correlation between dT_gc/dp determined experimentally from the high-pressure dielectric data and the Ehrenfest equation has been tested for the plastic crystals (anhGLU and anhMAN) for the first time. Interestingly, a satisfactory agreement was found between both approaches. It is a quite intriguing finding which can be rationalized by the fact that the studied materials are characterized by the low complexity (lower degrees of freedom with respect to the molecular mobility) as well as ordered internal structure. Therefore, one can speculate that in contrast to the ordinary glasses the dynamics of the plastic crystals might be described with the use of a single order parameter. However, to confirm this thesis further, pressure-volume-temperature (PVT) experiments enabling calculations of the Prigogine Defay ratio are required.

Communication: Effect of density on the physical aging of pressure-densified polymethylmethacrylate

The rate of physical aging of glassy polymethylmethacrylate (PMMA), followed from the change in the secondary relaxation with aging, is found to be independent of the density, the latter controlled by the pressure during glass formation. Thus, the aging behavior of the secondary relaxation is the same whether the glass is more compacted or less dense than the corresponding equilibrium liquid. This equivalence in aging of glasses formed under different pressures indicates that local packing is the dominant variable governing the glassy dynamics. The fact that pressure densification yields different glass structures is at odds with a model for non-associated materials having dynamic properties exhibited by PMMA, such as density scaling of the relaxation time and isochronal superposition of the relaxation dispersion.

A test for the existence of isomorphs in glass-forming materials

We describe a method to determine whether a material has isomorphs in its thermodynamic phase diagram. Isomorphs are state points for which various properties are invariant in reduced units. Such materials are commonly identified from strong correlation between thermal fluctuations of the potential energy, U, and the virial W, but this identification is not generally applicable to real materials. We show from molecular dynamic simulations of atomic, molecular, and polymeric materials that systems with strong U-W correlation cannot be pressure densified, that is, the density obtained on cooling to the glassy state and releasing the pressure is independent of the pressure applied during cooling.

Pair Potential That Reproduces the Shape of Isochrones in Molecular Liquids

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

Isochronal superposition and density scaling of the intermolecular dynamics in glass-forming liquids with varying hydrogen bonding propensity

We have tested the idea of thermodynamic scaling T⁻¹ρ^γ and isochronal superposition in glass-forming liquids with varying propensity to form hydrogen bonds.

Communication: High pressure specific heat spectroscopy reveals simple relaxation behavior of glass forming molecular liquid

The frequency dependent specific heat has been measured under pressure for the molecular glass forming liquid 5-polyphenyl-4-ether in the viscous regime close to the glass transition. The temperature and pressure dependences of the characteristic time scale associated with the specific heat is compared to the equivalent time scale from dielectric spectroscopy performed under identical conditions. It is shown that the ratio between the two time scales is independent of both temperature and pressure. This observation is non-trivial and demonstrates the existence of specially simple molecular liquids in which different physical relaxation processes are both as function of temperature and pressure/density governed by the same underlying "inner clock." Furthermore, the results are discussed in terms of the recent conjecture that van der Waals liquids, like the measured liquid, comply to the isomorph theory.

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.

Isomorphs in model molecular liquids

Isomorphs are curves in the phase diagram along which a number of static and dynamic quantities are invariant in reduced units (Gnan, N.; et al. J. Chem. Phys.2009, 131, 234504). A liquid has good isomorphs if and only if it is strongly correlating, i.e., if the equilibrium virial/potential energy fluctuations are more than 90% correlated in the NVT ensemble. Isomorphs were previously discussed with a focus on atomic systems. This paper generalizes isomorphs to liquids composed of rigid molecules and study the isomorphs of systems of small rigid molecules: the asymmetric dumbbell model, a symmetric inverse power-law dumbbell, and the Lewis-Wahnström o-terphenyl (OTP) model. For all model systems, the following quantities are found to a good approximation to be invariant along an isomorph: the isochoric heat capacity, the excess entropy, the reduced molecular center-of-mass self-part of the intermediate scattering function, and the reduced molecular center-of-mass radial distribution function. In agreement with theory, we also find that an instantaneous change of temperature and density from an equilibrated state point to an isomorphic state point leads to no relaxation. The isomorphs of the Lewis-Wahnström OTP model were found to be more approximative than those of the asymmetric dumbbell model; this is consistent with the OTP model being less strongly correlating. The asymmetric dumbbell and Lewis-Wahnström OTP models each have a "master isomorph"; i.e., the isomorphs have identical shape in the virial/potential energy phase diagram.

Dielectric and mechanical relaxation in isooctylcyanobiphenyl (8*OCB)

The dynamics of isooctylcyanobiphenyl (8*OCB) was characterized using dielectric and mechanical spectroscopies. This isomer of the liquid crystalline octylcyanobiphenyl (8OCB) vitrifies during cooling or on application of pressure, exhibiting the typical features of glass-forming liquids: non-Debye relaxation function, non-Arrhenius temperature dependence of the relaxation times, τ(α), a dynamic crossover at T ∼ 1.6T(g). This crossover is evidenced by changes in the behavior of both the peak shape and the temperature dependence of τ(α). The primary relaxation time at the crossover, 2 ns at ambient pressure, is the smallest value reported to date for any molecular liquid or polymer. Interestingly, at all temperatures below this crossover, τ(α) and the dc conductivity remain coupled (i.e., conform to the Debye-Stokes-Einstein relation). Two secondary relaxations are observed in the glassy state, one of which is identified as the Johari-Goldstein process. Unlike the case for 8OCB, no liquid crystalline phase could be attained for 8*OCB, demonstrating that relatively small differences in chemical structure can effect substantial changes in the intermolecular potential.

Pressure-energy correlations in liquids. IV. "Isomorphs" in liquid phase diagrams

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 phase 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 co-workers. 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, and it explains isochronal superposition for strongly correlating viscous liquids.