Density scaling and decoupling in o-terphenyl, salol, and dibutyphthalate

Combined description of pressure-volume-temperature and dielectric relaxation of several polymeric and low-molecular-weight organic glass-formers using SL-TS2 approach

We apply our recently-developed mean-field "SL-TS2" (two-state Sanchez-Lacombe) model to simultaneously describe dielectric α-relaxation time, τ_α, and pressure-volume-temperature (PVT) data in four polymers (polystyrene, poly(methylmethacrylate), poly(vinyl acetate) and poly(cyclohexane methyl acrylate)) and four organic molecular glass formers (ortho-terphenyl, glycerol, PCB-62, and PDE). Previously, it has been shown that for all eight materials, the Casalini-Roland thermodynamical scaling, τ_α = f(Tvγsp) (where T is temperature and v_sp is specific volume) is satisfied (R. Casalini and C. M. Roland, Phys. Rev. E, 2004, 69(6), 62501). It has also been previously shown that the same scaling emerges naturally (for sufficiently low pressures) within the "SL-TS2" framework (V. V. Ginzburg, Soft Matter, 2021, 17, 9094-9106). Here, we fit the ambient pressure curves for the relaxation time and the specific volume as functions of temperature for the eight materials and observe a good agreement between theory and experiment. We then use the Casalini-Roland scaling to convert those results into "master curves", thus enabling predictions of relaxation times and specific volumes at elevated pressures. The proposed approach can be used to describe other glass-forming materials, both low-molecular-weight and polymeric.

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.

Time-scale ordering in hydrogen- and van der Waals-bonded liquids

The time scales of structural relaxation are investigated on the basis of five different response functions for 1,2, 6-hexanetriol, a hydrogen-bonded liquid with a minor secondary contribution, and 2,6,10,15,19,23-hexamethyl-tetracosane (squalane), a van der Waals-bonded liquid with a prominent secondary relaxation process. Time scales of structural relaxation are derived as inverse peak frequencies for each investigated response function. For 1,2,6-hexanetriol, the ratios of the time scales are temperature-independent, while a decoupling of time scales is observed for squalane in accordance with the literature. An alternative evaluation approach is made on the squalane data, extracting time scales from the terminal relaxation mode instead of the peak position, and in this case, temperature-independent time-scale ratios are also found for squalane, despite its strong secondary relaxation contribution. Interestingly, the very same ordering of response-function-specific time scales is observed for these two liquids, which is also consistent with the observation made for simple van der Waals-bonded liquids reported previously [Jakobsen et al., J. Chem. Phys. 136, 081102 (2012)]. This time-scale ordering is based on the following response functions, from fast to slow dynamics: shear modulus, bulk modulus, dielectric permittivity, longitudinal thermal expansivity coefficient, and longitudinal specific heat. These findings indicate a general relation between the time scales of different response functions and, as inter-molecular interactions apparently play a subordinate role, suggest a rather generic nature of the process of structural relaxation.

Thermodynamic or density scaling of the thermal conductivity of liquids

Thermodynamic or density scaling is applied to thermal conductivity (λ) data from the literature for the model Lennard-Jones (12-6) fluid; the noble gases neon to xenon; nitrogen, ethene, and carbon dioxide as examples of linear molecules; the quasi-spherical molecules methane and carbon tetrachloride; the flexible chain molecules n-hexane and n-octane; the planar toluene and m-xylene; the cyclic methylcyclohexane; the polar R132a and chlorobenzene; and ammonia and methanol as H-bonded fluids. Only data expressed as Rosenfeld reduced properties could be scaled successfully. Two different methods were used to obtain the scaling parameter γ, one based on polynomial fits to the group (TV^γ) and the other based on the Avramov equation. The two methods agree well, except for λ of CCl₄. γ for the thermal conductivity is similar to those for the viscosity and self-diffusion coefficient for the smaller molecules. It is significantly larger for the Lennard-Jones fluid, possibly due to a different dependence on packing fraction, and much larger for polyatomic molecules where heat transfer through internal modes may have an additional effect. Methanol and ammonia, where energy can be transmitted through intermolecular hydrogen bonding, could not be scaled. This work is intended as a practical attempt to examine thermodynamic scaling of the thermal conductivity of real fluids. The divergence of the scaling parameters for different properties is unexpected, suggesting that refinement of theory is required to rationalize this result. For the Lennard-Jones fluid, the Ohtori-Iishi version of the Stokes-Einstein-Sutherland relation applies at high densities in the liquid and supercritical region.

Density Scaling Based Detection of Thermodynamic Regions of Complex Intermolecular Interactions Characterizing Supramolecular Structures

In this paper, applying the density scaling idea to an associated liquid 4-methyl-2-pentanol used as an example, we identify different pressure-volume-temperature ranges within which molecular dynamics is dominated by either complex H-bonded networks most probably leading to supramolecular structures or non-specific intermolecular interactions like van der Waals forces. In this way, we show that the density scaling law for molecular dynamics near the glass transition provides a sensitive tool to detect thermodynamic regions characterized by intermolecular interactions of different type and complexity for a given material in the wide pressure-volume-temperature domain even if its typical form with constant scaling exponent is not obeyed. Moreover, we quantify the observed decoupling between dielectric and mechanical relaxations of the material in the density scaling regime. The suggested methods of analyses and their interpretations open new prospects for formulating models based on proper effective intermolecular potentials describing physicochemical phenomena near the glass transition.

On the pressure dependence of the thermodynamical scaling exponent γ

Since its initial discovery more than fifteen years ago, the thermodynamical scaling of the dynamics of supercooled liquids has been used to provide many new important insights in the physics of liquids, particularly on the link between dynamics and intermolecular potential. A question that has long been discussed is whether the scaling exponent γ_S is a constant or does it depends on pressure. An alternative definition of the scaling parameter, γ_I = ∂ ln T/∂ ln ρ|_X has been presented in the literature, and has been erroneously considered equivalent to γ_S. Here we offer a simple method to determine the pressure dependence of γ_I using only the pressure dependence of the glass transition and the equation of state. Using this new method we find that for the six nonassociated liquids investigated, γ_I always decreases with increasing pressure. Importantly in all cases the value of γ_I remains always larger than 4. Liquids having γ_I closer to 4 at low pressure show a smaller change in γ_I with pressure. We argue that this result has very important consequences for the experimental determination of the functional form of the repulsive part of the potential in liquids. Comparing the pressure and temperature dependence of γ_S and γ_I we find, contrary to what has been assumed in the literature to date, that these two parameters are not equivalent and have very different pressure and temperature dependences.

Temperature and pressure dependence of the alpha relaxation in ortho-terphenyl

Molecular dynamics (MD) simulations of ortho-terphenyl using an all-atom model with the optimized potentials for liquid simulations (OPLS) force field were performed both in the high temperature Arrhenian region and at lower temperatures that include the onset of the super-Arrhenian region. From the MD simulations, the internal energy of both the equilibrium liquid and crystal was determined from 300 K to 600 K and at pressures from 0.1 MPa to 1 GPa. The translational and rotational diffusivities were also determined at these temperatures and pressures for the equilibrium liquid. It is shown that within a small offset, the excess internal energy Ū_x from the MD simulations is consistent with the experimentally determined excess internal energy reported earlier [Caruthers and Medvedev, Phys. Rev. Mater. 2, 055604, (2018)]. The MD mobility data {including extremely long-time 1 atm simulations from the study by Eastwood et al. [J. Phys. Chem. B 117, 12898, (2013)]} were combined with experimental data to form a unified dataset, where it was shown that in both the high temperature Arrhenian region and the lower temperature super-Arrhenian region, the mobility is a linear function of 1/Ū_x(T,p), albeit with different proportionality constants. The transition between the Arrhenian and super-Arrhenian regions is relatively sharp at a critical internal energy Ū_x ^α. The 1/Ū_x(T,p) model is able to describe the mobility data over nearly 16 orders-of-magnitude. Other excess thermodynamic properties such as excess enthalpy and excess entropy (i.e., the Adam-Gibbs model) are unable to unify the pressure dependence of the mobility.

Exploring the connection between the density-scaling exponent and the intermolecular potential for liquids on the basis of computer simulations of quasireal model systems

In this paper, based on the molecular dynamics simulations of quasireal model systems, we propose a method for determination of the effective intermolecular potential for real materials. We show that in contrast to the simple liquids, the effective intermolecular potential for the studied systems depends on the thermodynamic conditions. Nevertheless, the previously established relationship for simple liquids between the exponent of the inverse power law approximation of intermolecular potential and the density-scaling exponent is still preserved when small enough intermolecular distances are considered. However, our studies show that molecules approach each other at these very short distances relatively rarely. Consequently, only sparse interactions between extremely close molecules determine the value of the scaling exponent and then strongly influence the connection between dynamics and thermodynamics of the whole system.

A model-free analysis of configurational properties to reduce the temperature- and pressure-dependent segmental relaxation times of polymers

The segmental relaxation time data for poly(vinyl acetate), poly(vinyl chloride), and linear and star polystyrene are analyzed using a model-free method to determine how the temperature- and pressure-dependent relaxation times, τ, scale with the relative configurational thermodynamic properties. The model-free method assumes no specific mathematical form, such as reciprocal linearity, and the configurational properties are referred to an isochronal state to eliminate the bias associated with the definition of the ideal glassy state. The scaling ability of a given configurational property is strongly material-dependent with the logarithm of τ scaling better with TS_c and H_c for poly(vinyl acetate), with TS_c, H_c, and U_c for poly(vinyl chloride), and with TS_c, H_c, and V_c for linear and star polystyrene. The choice of the isochronal reference state does not qualitatively affect the results.

On the experimental determination of the repulsive component of the potential from high pressure measurements: What is special about twelve?

In this paper, we present an overview of results in the literature regarding the thermodynamical scaling of the dynamics of liquids and polymers as measured from high-pressure measurements. Specifically, we look at the scaling exponent γ and argue that it exhibits the limiting behavior γ → 4 in regimes for which molecular interactions are dominated by the repulsive part of the intermolecular potential. For repulsive potentials of the form U(r) ∝ r^-n, γ has been found to be related to the exponent n via the relation γ = n/3. Therefore, this limiting behavior for γ would suggest that a large number of molecular systems may be described by a common repulsive potential U(r) ∝ r^-n with n ≈ 12.

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.

Glass Transition Temperature and Density Scaling in Cumene at Very High Pressure

We present a new method that allows direct measurements of the glass transition temperature T_{g} at pressures up to 4.55 GPa in the glass-forming liquid cumene (isopropylbenzene). This new method uses a diamond anvil cell and can measure T_{g} at pressures of 10 GPa or greater. Measuring T_{g} at the glass→liquid transition involves monitoring the disappearance of pressure gradients initially present in the glass, but also takes advantage of the large increase in the volume expansion coefficient α_{p} at T_{g} as the supercooled or superpressed liquid is entered. Accurate T_{g}(P) values in cumene allow us to show that density scaling holds along this isochronous line up to pressures much higher than any previous study, corresponding to a density increase of 29%. Our results for cumene over this huge compression range yield ρ^{γ}/T=C, where C is a constant and where γ=4.77±0.02 for this nonassociated glass-forming system. Finally, high-pressure cumene viscosity data from the literature taken at much lower pressures and at several different temperatures, corresponding to a large dynamic range of nearly 13 orders of magnitude, are shown to superimpose on a plot of η vs ρ^{γ}/T for the same value of γ.

Structural relaxation and thermal conductivity of high-pressure formed, high-density di-n-butyl phthalate glass and pressure induced departures from equilibrium state

We report a study of structural relaxation of high-density glasses of di-n-butyl phthalate (DBP) by measuring thermal conductivity, κ, under conditions of pressure and temperature (p,T) designed to modify both the vibrational and configurational states of a glass. Various high-density glassy states of DBP were formed by (i) cooling the liquid under a fixed high p and partially depressurizing the glass, (ii) isothermal annealing of the depressurized glass, and (iii) pressurizing the glass formed by cooling the liquid under low p. At a given low p, κ of the glass formed by cooling under high p is higher than that of the glass formed by cooling under low p, and the difference increases as glass formation p is increased. κ of the glass formed under 1 GPa is ∼20% higher at ambient p than κ of the glass formed at ambient p. On heating at low p, κ decreases until the glass to liquid transition range is reached. This is the opposite of the increase in κ observed when a glass formed under a certain p is heated under the same p. At a given high p, κ of the low-density glass formed by cooling at low p is lower than that of the high-density glass formed by cooling at that high p. On heating at high p, κ increases until the glass to liquid transition range is reached. The effects observed are due to a thermally assisted approach toward equilibrium at p different from the glass formation p. In all cases, the density, enthalpy, and entropy would change until the glasses become metastable liquids at a fixed p, thus qualitatively relating κ to variation in these properties.

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.