On the pressure dependence of the thermodynamical scaling exponent γ

Insight into properties of sizable glass former from volumetric measurements

Sizable glass formers feature numerous unique properties and potential applications, but many questions regarding their glass transition dynamics have not been resolved yet. Here, we have analyzed structural relaxation times measured as a function of temperature and pressure in combination with the equation of state obtained from pressure-volume-temperature measurements. Despite evidence from previous dielectric studies indicating a remarkable sensitivity of supercooled dynamics to compression, and contrary to intuition, our results demonstrated the proof for the almost equivalent importance of thermal energy and free volume fluctuations in controlling reorientation dynamics of sizable molecules. The found scaling exponent γ = 3.0 and Ev/Ep ratio of 0.6 were typical for glass-forming materials with relaxation dynamics determined by both effects with a minor advantage of thermal fluctuations involvement. It shows that the high values of key parameters characterizing the sensitivity of the glass transition dynamics to pressure changes, i.e., activation volume ΔV and dTg/dP, are not a valid premise for a remarkable contribution of volume to glass transition dynamics.

Role of anisotropy in understanding the molecular grounds for density scaling in dynamics of glass-forming liquids

Molecular Dynamics (MD) simulations of glass-forming liquids play a pivotal role in uncovering the molecular nature of the liquid vitrification process. In particular, much focus was given to elucidating the interplay between the character of intermolecular potential and molecular dynamics behaviour. This has been tried to achieve by simulating the spherical particles interacting via isotropic potential. However, when simulation and experimental data are analysed in the same way by using the density scaling approaches, serious inconsistency is revealed between them. Similar scaling exponent values are determined by analysing the relaxation times and pVT data obtained from computer simulations. In contrast, these values differ significantly when the same analysis is carried out in the case of experimental data. As discussed thoroughly herein, the coherence between results of simulation and experiment can be achieved if anisotropy of intermolecular interactions is introduced to MD simulations. In practice, it has been realized in two different ways: (1) by using the anisotropic potential of the Gay-Berne type or (2) by replacing the spherical particles with quasi-real polyatomic anisotropic molecules interacting through isotropic Lenard-Jones potential. In particular, the last strategy has the potential to be used to explore the relationship between molecular architecture and molecular dynamics behaviour. Finally, we hope that the results presented in this review will also encourage others to explore how 'anisotropy' affects remaining aspects related to liquid-glass transition, like heterogeneity, glass transition temperature, glass forming ability, etc.

The origin of the density scaling exponent for polyatomic molecules and the estimation of its value from the liquid structure

In this article, we unravel the problem of interpreting the density scaling exponent for the polyatomic molecules representing the real van der Waals liquids. Our studies show that the density scaling exponent is a weighted average of the exponents of the repulsive terms of all interatomic interactions that occur between molecules, where the potential energy of a given interaction represents its weight. It implies that potential energy is a key quantity required to calculate the density scaling exponent value for real molecules. Finally, we use the well-known method for potential energy estimation and show that the density scaling exponent could be successfully predicted from the liquid structure for fair representatives of the real systems.

Density Scaling of Translational and Rotational Molecular Dynamics in a Simple Ellipsoidal Model near the Glass Transition

In this paper, we show that a simple anisotropic model of supercooled liquid properly reflects some density scaling properties observed for experimental data, contrary to many previous results obtained from isotropic models. We employ a well-known Gay–Berne model earlier parametrized to achieve a supercooling and glass transition at zero pressure to find the point of glass transition and explore volumetric and dynamic properties in the supercooled liquid state at elevated pressure. We focus on dynamic scaling properties of the anisotropic model of supercooled liquid to gain a better insight into the grounds for the density scaling idea that bears hallmarks of universality, as follows from plenty of experimental data collected near the glass transition for different dynamic quantities. As a result, the most appropriate values of the scaling exponent γ are established as invariants for a given anisotropy aspect ratio to successfully scale both the translational and rotational relaxation times considered as single variable functions of density^γ/temperature. These scaling exponent values are determined based on the density scaling criterion and differ from those obtained in other ways, such as the virial–potential energy correlation and the equation of state derived from the effective short-range intermolecular potential, which is qualitatively in accordance with the results yielded from experimental data analyses. Our findings strongly suggest that there is a deep need to employ anisotropic models in the study of glass transition and supercooled liquids instead of the isotropic ones very commonly exploited in molecular dynamics simulations of supercooled liquids over the last decades.

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

This paper studies numerically the Weeks-Chandler-Andersen system, which is shown to obey hidden scale invariance with a density-scaling exponent that varies from below 5 to above 500. This unprecedented variation makes it advantageous to use the fourth-order Runge-Kutta algorithm for tracing out isomorphs. Good isomorph invariance of structure and dynamics is observed over more than three orders of magnitude temperature variation. For all state points studied, the virial potential-energy correlation coefficient and the density-scaling exponent are controlled mainly by the temperature. Based on the assumption of statistically independent pair interactions, a mean-field theory is developed that rationalizes this finding and provides an excellent fit to data at low temperatures.

A Simple New Way To Account for Free Volume in Glassy Dynamics: Model-Free Estimation of the Close-Packed Volume from PVT Data

In this article we focus on the important role of well-defined free volume (V_free) in dictating the structural relaxation times, τ, of glass-forming liquids and polymer melts. Our definition of V_free = V - V_hc, where V is the total system volume, means the use of V_free depends on determination of V_hc, the system's volume in the limiting closely packed state. Rejecting the historically compromised use of V_free as a dynamics-dependent fitting function, we have successfully applied a clear thermodynamics-based route to V_hc using the locally correlated lattice (LCL) model equation of state (EOS). However, in this work we go further and show that V_hc can be defined without the use of an equation of state by direct linear extrapolation of a V(T) high-pressure isobar down to zero temperature (T). The results from this route, tested on a dozen experimental systems, yield ln τ vs 1/V_free isotherms that are linear with T-dependent slopes, consistent with the general ln τ ∼ f(T) × (1/V_free) form of behavior we have previously described. This functional form also results by implementing a simple mechanistic explanation via the cooperative free volume (CFV) rate model, which assumes that dynamic relaxation is both thermally activated and that it requires molecular segmental cooperativity. With the degree of the latter, and thus the activation energy, being determined by the availability of free volume, the new route we demonstrate here for determination of V_free expands the potential for understanding and predicting local dynamic relaxation in glass-forming materials.

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.