Communication: Relationships between Intermolecular potential, thermodynamics, and dynamic scaling in viscous systems

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.

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.

The effect of molecular architecture on the physical properties of supercooled liquids studied by MD simulations: Density scaling and its relation to the equation of state

Theoretical concepts in condensed matter physics are typically verified and also developed by exploiting computer simulations mostly in simple models. Predictions based on these usually isotropic models are often at odds with measurement results obtained for real materials. One of the examples is an intriguing problem within the density scaling idea that has attracted attention in recent decades due to its hallmarks of universality, i.e., the fact that the difference between the density scaling exponent and the exponent of the equation of state is observed for real materials, whereas it has not been reported for the model system. In this paper, we use new model molecules of simple but anisotropic architecture to study the effect of molecular anisotropy on the dynamic and thermodynamic properties of the system. We identify the applicable range of intermolecular interactions for a given physical process, and then we explain the reason for observed differences between the behavior of the model and real systems. It demonstrates that the new model systems open broad perspectives for simulation and theoretical research, for example, into unifying concepts in the glass transition physics.

Determination of derived volumetric properties and heat capacities at high pressures using two density scaling based equations of state. Application to dipentaerythritol hexa(3,5,5-trimethylhexanoate)

Scaling based EoSs describe the complex behavior of derived properties for broad temperature and pressure ranges from diPEiC₉ experimental densities.

In search of invariants for viscous liquids in the density scaling regime: investigations of dynamic and thermodynamic moduli

In this paper, we report the nontrivial results of our investigations of dynamic and thermodynamic moduli in search of invariants for viscous liquids in the density scaling regime by using selected supercooled van der Waals liquids as representative materials. Previously, the dynamic modulus M_p-T (defined in the pressure-temperature representation by the ratio of isobaric activation energy and activation volume) as well as the ratio B_T/M_p-T (where B_T is the thermodynamic modulus defined as the inverse isothermal compressibility) have been suggested as some kinds of material constants. We have established that they are not valid in the explored wide range of temperatures T over a dozen decades of structural relaxation times τ. The temperature dependences of M_p-T and B_T/M_p-T have been elucidated by comparison with the well-known measure of the relative contribution of temperature and density fluctuations to molecular dynamics near the glass transition, i.e., the ratio of isochoric and isobaric activation energies. Then, we have implemented an idea to transform the definition of the dynamic modulus M_p-T from the p-T representation to the V-T one. This idea relied on the disentanglement of combined temperature and density fluctuations involved in isobaric parameters and has resulted in finding an invariant for viscous liquids in the density scaling regime, which is the ratio of thermodynamic and dynamic moduli, B_T/M_V-T. In this way, we have constituted a characteristic of thermodynamics and molecular dynamics, which remains unchanged in the supercooled liquid state for a given material, the molecular dynamics of which obeys the power density scaling law.

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 γ.

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.

Recent developments in the experimental investigations of relaxations in pharmaceuticals by dielectric techniques at ambient and elevated pressure

In recent years, there is a growing interest in improving the physicochemical stability of amorphous pharmaceutical solids due to their very promising applications to manufacture medicines characterized by a better water solubility, and consequently by a higher dissolution rate than those of their crystalline counterparts. In this review article, we show that the molecular mobility investigated both in the supercooled liquid and glassy states is the crucial factor required to understand molecular mechanisms that govern the physical stability of amorphous drugs. We demonstrate that pharmaceuticals can be thoroughly examined by means of the broadband dielectric spectroscopy, which is a very useful experimental technique to explore different relaxation processes and crystallization kinetics as well. Such studies conducted in the wide temperature and pressure ranges provide data needed in searching correlations between properties of molecular dynamics and crystallization process, which are aimed at developing effective and efficient methods for stabilizing amorphous drugs.

Adam-Gibbs model in the density scaling regime and its implications for the configurational entropy scaling

To solve a long-standing problem of condensed matter physics with determining a proper description of the thermodynamic evolution of the time scale of molecular dynamics near the glass transition, we have extended the well-known Adam-Gibbs model to describe the temperature-volume dependence of structural relaxation times, τ_α(T, V). We also employ the thermodynamic scaling idea reflected in the density scaling power law, τ_α = f(T⁻¹V^−γ), recently acknowledged as a valid unifying concept in the glass transition physics, to differentiate between physically relevant and irrelevant attempts at formulating the temperature-volume representations of the Adam-Gibbs model. As a consequence, we determine a straightforward relation between the structural relaxation time τ_α and the configurational entropy S_C, giving evidence that also S_C(T, V) = g(T⁻¹V^−γ) with the exponent γ that enables to scale τ_α(T, V). This important findings have meaningful implications for the connection between thermodynamics and molecular dynamics near the glass transition, because it implies that τ_α can be scaled with S_C.

Role of entropy in the thermodynamic evolution of the time scale of molecular dynamics near the glass transition

In this paper, we investigate how changes in the system entropy influence the characteristic time scale of the system molecular dynamics near the glass transition. Independently of any model of thermodynamic evolution of the time scale, against some previous suppositions, we show that the system entropy S is not sufficient to govern the time scale defined by structural relaxation time τ. In the density scaling regime, we argue that the decoupling between τ and S is a consequence of different values of the scaling exponents γ and γ(S) in the density scaling laws, τ=f(ρ(γ)/T) and S=h(ρ(γ(S))/T), where ρ and T denote density and temperature, respectively. It implies that the proper relation between τ and S requires supplementing with a density factor, u(ρ), i.e., τ=g(u(ρ)w(S)). This meaningful finding additionally demonstrates that the density scaling idea can be successfully used to separate physically relevant contributions to the time scale of molecular dynamics near the glass transition. The relation reported by us between τ and S constitutes a general pattern based on nonconfigurational quantities for describing the thermodynamic evolution of the characteristic time scale of molecular dynamics near the glass transition in the density scaling regime, which is a promising alternative to the approaches based as the Adam-Gibbs model on the configurational entropy that is difficult to evaluate in the entire thermodynamic space. As an example, we revise the Avramov entropic model of the dependence τ(T,ρ), giving evidence that its entropic basis has to be extended by the density dependence of the maximal energy barrier for structural relaxation. We also discuss the excess entropy S(ex), the density scaling of which is found to mimic the density scaling of the total system entropy S.

Recent progress on dielectric properties of protic ionic liquids

Protic ionic liquids (PILs) are key materials for a wide range of emerging technologies. In particular, these systems have long been envisioned as promising candidates for fuel cells. Therefore, in recent years special attention has been devoted to thorough studies of these compounds. Amongst others, dielectric properties of PILs at ambient and elevated pressure have become the subject of intense research. The reason for this lies in the role of broadband dielectric spectroscopy in recognizing the conductivity mechanism in protic ionic systems. In this paper, we summarize the dielectric results of various PILs reflecting recent advances in this field.

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).

Effect of temperature and density fluctuations on the spatially heterogeneous dynamics of glass-forming Van der Waals liquids under high pressure

In this Letter, we show how temperature and density fluctuations affect the spatially heterogeneous dynamics at ambient and elevated pressures. By using high-pressure experimental data for van der Waals liquids, we examine contributions of the temperature and density fluctuations to the dynamics heterogeneity. We show that the dynamic heterogeneity decreases significantly with increasing pressure at a constant structural relaxation time (isochronal condition), while the broadening of the relaxation spectrum remains constant. This observation questions the relationship between spectral broadening and dynamic heterogeneity.

Systematic coarse-graining of potential energy landscapes and dynamics in liquids

Recent efforts have shown that the dynamic properties of a wide class of liquids can be mapped onto semi-universal scaling laws and constitutive relations that are motivated by thermodynamic analyses of much simpler models. In particular, it has been found that many systems exhibit dynamics whose behavior in state space closely follows that of soft-sphere particles interacting through an inverse power repulsion. In the present work, we show that a recently developed coarse-graining theory provides a natural way to understand how arbitrary liquids can be mapped onto effective soft-sphere models and hence how one might potentially be able to extract underlying dynamical scaling laws. The theory is based on the relative entropy, an information metric that quantifies how well a soft-sphere approximation to a liquid's multidimensional potential energy landscape performs. We show that optimization of the relative entropy not only enables one to extract effective soft-sphere potentials that suggest an inherent scaling of thermodynamic and dynamic properties in temperature-density space, but that also has rather interesting connections to excess entropy based theories of liquid dynamics. We apply the approach to a binary mixture of Lennard-Jones particles, and show that it gives effective soft-sphere scaling laws that well-describe the behavior of the diffusion constants. Our results suggest that the relative entropy formalism may be useful for "perturbative" type theories of dynamics, offering a general strategy for systematically connecting complex energy landscapes to simpler reference ones with better understood dynamic behavior.

Pressure coefficient of the glass transition temperature in the thermodynamic scaling regime

We report that the pressure coefficient of the glass transition temperature, dT(g)/dp, which is commonly used to determine the pressure sensitivity of the glass transition temperature T(g), can be predicted in the thermodynamic scaling regime. We show that the equation derived from the isochronal condition combined with the well-known scaling, TV(γ) = const, predicts successfully values of dT(g)/dp for a variety of glass-forming systems, including van der Waals liquids, polymers, and ionic liquids.

Temperature-Volume Entropic Model for Viscosities and Structural Relaxation Times of Glass Formers

In this Letter, an entropic model recently formulated by Mauro et al. for the temperature dependence of viscosity in glass-forming materials is generalized to describe the temperature-volume dependences of viscosities and structural relaxation times near the glass transition. It is found that the generalization shows limitations of its temperature precursor. The extended model describes well the structural dielectric relaxation times τα(T,V) of supercooled van der Waals liquids. The obtained results are discussed in the context of the thermodynamic scaling law for molecular dynamics of viscous systems.

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.

Pair force distributions in simple fluids

Analytic expressions are derived for the frequency distribution, P(f), of pair forces, f, and those of their α-Cartesian component, f(α), or P(f(α)), for some typical model simple fluids, expressed in terms of the radial distribution function and known constants. For strongly repulsive inverse power (IP), exponential and Yukawa purely repulsive potentials, P(f) diverges at the origin approximately as ∼f(-1), but with different limiting analytic forms. P(f(α)) is also shown to diverge as ∼f(-1) as f → 0 for the IP fluid. For the Lennard-Jones potential fluid, P(f) is finite for all f ≥ 0 but has two singularities for negative f, corresponding to the zero force limit (i.e., f → 0(-)) and the point of inflection in the potential. The corresponding component force distribution is singular as f(α) → 0 from both positive and negative force sides. The large force limit of P(f), which originates from the close neighbor interactions, is nearly exponential for the IP and LJ fluids, as is also found for granular materials. A more complete picture of force distributions in off-lattice particulate systems as a function of force law and state point (particularly the extent of "thermalization" of the particles) is provided.

Communication: Thermodynamic scaling of the Debye process in primary alcohols

The molecular dynamics of hydrogen-bonded liquids usually does not satisfy the thermodynamic scaling. However, very recently, two opposite conclusions about validity of thermodynamical scaling in monohydroxy alcohol, 2-ethyl-1-hexanol, were presented by Reiser et al. [J. Chem. Phys. 132, 181101 (2010)] and Fragiadakis et al. [J. Chem. Phys. 132, 144505 (2010)]. In this communication we present new experimental results that can explain this ostensible contradiction.

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).