Experimental Evidence for a State-Point-Dependent Density-Scaling Exponent of Liquid 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.

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.

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.

Virial-potential-energy correlation and its relation to density scaling for quasireal model systems

In this paper, we examine the virial- and the potential-energy correlation for quasireal model systems. This correlation constitutes the framework of the theory of the isomorph in the liquid phase diagram commonly examined using simple liquids. Interestingly, our results show that for the systems characterized by structural anisotropy and flexible bonds, the instantaneous values of total virial and total potential energy are entirely uncorrelated. It is due to the presence of the intramolecular interactions because the contributions to the virial and potential energy resulting from the intermolecular interactions still exhibit strong linear dependence. Interestingly, in contrast to the results reported for simple liquids, the slope of the mentioned linear dependence is different than the values of the density scaling exponent. However, our findings show that for quasireal materials, the slope of dependence between the virial and potential energy (resulting from the intermolecular interactions) strongly depends on the interval of intermolecular distances that are taken into account. Consequently, the value of the slope of the discussed relationship, which enables satisfactory density scaling, can be obtained. Interestingly, this conclusion is supported by the results obtained for analogous systems without intermolecular attraction, for which the value the slope of the virial-potential-energy correlation is independent of considered intermolecular distances, directly corresponds to the exponent of the intermolecular repulsion, and finally leads to accurate density scaling.

Hidden Scale Invariance in Polydisperse Mixtures of Exponential Repulsive Particles

Polydisperse systems of particles interacting by the purely repulsive exponential (EXP) pair potential are studied in regard to how structure and dynamics vary along isotherms, isochores, and isomorphs. The sizable size polydispersities of 23%, 29%, 35%, and 40%, as well as energy polydispersity 35%, were considered. For each system an isomorph was traced out covering about one decade in density. For all systems studied, the structure and dynamics vary significantly along the isotherms and isochores but are invariant to a good approximation along the isomorphs. We conclude that the single-component EXP system's hidden scale invariance (implying isomorph invariance of structure and dynamics) is maintained even when a sizable polydispersity is introduced into the system.

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.

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.

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.

Experimental Evidence for a State-Point-Independent Density-Scaling Exponent in Ionic Liquids

This Letter addresses a fundamental issue of condensed-matter physics, which is the validity of the density-scaling concept. For this purpose, the ambient and high-pressure conductivity measurements of two selected ionic liquids (ILs), with the different contribution of H-bonding interactions, were performed in the dynamic range of 13 orders of magnitude and corresponding to the density changes as large as 20%. All experimental data obtained within one compound are shown to superimpose each other when plotted as a function of ρ^{γ}/T. These results clearly show that for studied ILs the scaling exponent is a state-point-independent parameter that is in odds with the recent findings for van der Waals liquid [Sanz et al., Phys. Rev. Lett. 122, 055501 (2019)PRLTAO0031-900710.1103/PhysRevLett.122.055501].