Comment on "Density scaling of the diffusion coefficient at various pressures in viscous liquids"

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.

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.

Activation volume of selected liquid crystals in the density scaling regime

In this paper, we demonstrate and thoroughly analyze the activation volumetric properties of selected liquid crystals in the nematic and crystalline E phases in comparison with those reported for glass-forming liquids. In the analysis, we have employed and evaluated two entropic models (based on either total or configurational entropies) to describe the longitudinal relaxation times of the liquid crystals in the density scaling regime. In this study, we have also exploited two equations of state: volumetric and activation volumetric ones. As a result, we have established that the activation volumetric properties of the selected liquid crystals are quite opposite to such typical properties of glass-forming materials, i.e., the activation volume decreases and the isothermal bulk modulus increases when a liquid crystal is isothermally compressed. Using the model based on the configurational entropy, we suggest that the increasing pressure dependences of the activation volume in isothermal conditions and the negative curvature of the pressure dependences of isothermal longitudinal relaxation times can be related to the formation of antiparallel doublets in the examined liquid crystals. A similar pressure effect on relaxation dynamics may be also observed for other material groups in case of systems, the molecules of which form some supramolecular structures.

The complex, non-monotonic thermal response of the volumetric space of simple liquids

We show that a non-monotonic solution of the equation ∂α_p(p,T)/∂T = 0 divides the phase diagram of simple liquids into two parts.

On the scaling behavior of electric conductivity in [C4mim][NTf2]

In this work we examine, for the first time, the molar conductivity behavior of the deeply supercooled room temperature ionic liquid [C₄mim][NTf₂] in the temperature, pressure and volume thermodynamic space in terms of density scaling regime (TV^γ)⁻¹ combined with the equation of state (EOS).

Scaling of volumetric data in model systems based on the Lennard-Jones potential

The crucial problem for better understanding the nature of glass transition and related relaxation phenomena is to find proper interrelations between the molecular dynamics and thermodynamics of viscous systems. To make progress towards this goal the recently observed density scaling of viscous liquid dynamics has been very intensively and successfully studied in the past few years. However, previous attempts at related scaling of volumetric data yielded results inconsistent with those found from the density scaling of molecular dynamics. In this paper, we show that volumetric data obtained from simulations in simple molecular models based on the Lennard-Jones (LJ) potential, such as the Kob-Andersen binary LJ liquid, its repulsive inverse power-law version, and the Lewis-Wahnström o-terphenyl model, can be scaled by using the same value of the exponent, which scales dynamic quantities and is directly related to the exponent of the repulsive inverse power law that underlies short-range approximations of the LJ potential.

The role of the isothermal bulk modulus in the molecular dynamics of super-cooled liquids

Elastic models imply that the energy expended for a flow event in ultra-viscous matter coincides with the elastic work required for deforming and re-arranging the environment of the moving entity. This is quite promising for explaining the strong non-Arrhenius behavior of dynamic quantities of fragile super-cooled liquids. We argue that the activation volume obtained from dielectric relaxation and light-scattering experiments for super-cooled liquids should scale with the Gibbs free energy of activation, with a proportionality constant determined by the isothermal bulk modulus and its pressure derivative, as described by an earlier thermodynamic elastic model. For certain super-cooled liquids the bulk compression transpiring in the local environment, as governed by the isothermal bulk modulus, play a significant role in the reorientational dynamics, with far-field density fluctuations and volume changes avoided by shear deformation.

Density scaling in viscous systems near the glass transition

In this paper, a general equation of state (EOS) valid for fluids in the vicinity of the glass transition is derived on the basis of its isothermal precursor. This EOS is able to predict the density scaling of both isobaric and isothermal PVT data and it explicitly involves the scaling exponent γ(EOS), which is most likely straightforwardly related to the exponent of the inverse power law of some effective potential valid for viscous systems. This EOS and the density scaling are very successfully tested for representatives of several material classes (van der Waals liquids, polymer melts, ionic liquids, and even strongly hydrogen-bonded systems). Additionally, if the thermodynamic scaling of primary relaxation times can be achieved with the scaling exponent γ for a given material, then the value γ(EOS) found from fitting its PVT data to the EOS enables us to evaluate the value γ, which is always considerably smaller than γ(EOS).