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

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.

Entropy Scaling of Molecular Dynamics in a Prototypical Anisotropic Model near the Glass Transition

Dynamics and thermodynamics of molecular systems in the vicinity of the boundary between thermodynamically nonequilibrium glassy and metastable supercooled liquid states are still incompletely explored and their theoretical and simulation models are imperfect despite many previous efforts. Among them, the role of total system entropy, configurational entropy, and excess entropy in the temperature-pressure or temperature-density dependence of global molecular dynamics (MD) timescale relevant to the glass transition is an essential topic intensively studied for over half of a century. By exploiting a well-known simple ellipsoidal model recently successfully applied to simulate the supercooled liquid state and the glass transition, we gain a new insight into the different views on the relationship between entropy and relaxation dynamics of glass formers, showing the molecular grounds for the entropy scaling of global MD timescale. Our simulations in the anisotropic model of supercooled liquid, which involves only translational and rotational degrees of freedom, give evidence that the total system entropy is sufficient to scale global MD timescale. It complies with the scaling effect on relaxation dynamics exerted by the configurational entropy defined as the total entropy diminished by vibrational contributions, which was earlier discovered for measurement data collected near the glass transition. Moreover, we argue that such a scaling behavior is not possible to achieve by using the excess entropy that is in excess of the ideal gas entropy, which is contrary to the results earlier suggested within the framework of simple isotropic models of supercooled liquids. Thus, our findings also warn against an excessive reliance on isotropic models in theoretical interpretations of molecular phenomena, despite their simplicity and popularity, because they may reflect improperly various physicochemical properties of glass formers.

Breakdown of the Simple Arrhenius Law in the Normal Liquid State

It is common practice to discuss the temperature effect on molecular dynamics of glass formers above the melting temperature in terms of the Arrhenius law. Using dielectric spectroscopy measurements of dc conductivity and structural relaxation time on the example of the typical glass former propylene carbonate, we provide experimental evidence that this practice is not justified. Our conclusions are supported by employing thermodynamic density scaling and the occurrence of inflection points in isothermal dynamic data measured at elevated pressure. Additionally, we propose a more suitable approach to describe the dynamics both above and below the inflection point based on a modified MYEGA model.

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.

Thermodynamic scaling of relaxation: insights from anharmonic elasticity

Using molecular dynamics simulations of a molecular liquid, we investigate the thermodynamic scaling (TS) of the structural relaxation time [Formula: see text] in terms of the quantity [Formula: see text], where T and ρ are the temperature and density, respectively. The liquid does not exhibit strong virial-energy correlations. We propose a method for evaluating both the characteristic exponent [Formula: see text] and the TS master curve that uses experimentally accessible quantities that characterise the anharmonic elasticity and does not use details about the microscopic interactions. In particular, we express the TS characteristic exponent [Formula: see text] in terms of the lattice Grüneisen parameter [Formula: see text] and the isochoric anharmonicity [Formula: see text]. An analytic expression of the TS master curve of [Formula: see text] with [Formula: see text] as the key adjustable parameter is found. The comparison with the experimental TS master curves and the isochoric fragilities of 34 glassformers is satisfying. In a few cases, where thermodynamic data are available, we test (i) the predicted characteristic exponent [Formula: see text] and (ii) the isochoric anharmonicity [Formula: see text], as drawn by the best fit of the TS of the structural relaxation, against the available thermodynamic data. A linear relation between the isochoric fragility and the isochoric anharmonicity [Formula: see text] is found and compared favourably with the results of experiments with no adjustable parameters. A relation between the increase of the isochoric vibrational heat capacity due to anharmonicity and the isochoric fragility is derived.

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.

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.

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.

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.