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

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.

Structurally Related Scaling Behavior in Ionic Systems

We examine the density scaling properties of two ionic materials, a classic aprotic low molecular weight ionic liquid, 1-butyl-3-methylimidazolium bis(perfluoroethylsulfonyl)imide ([BMIm][BETI]), and a polymeric ionic liquid, poly(3-methyl-1,2,3-triazolium bis(trifluoromethylsulfonyl)imide) (TPIL). Density scaling is known to apply rigorously to simple liquids lacking specific intermolecular associations such as hydrogen bonds. Previous work has found that ionic liquids conform to density scaling over limited ranges of temperature and pressure. In this work, we find that the dc-conductivity of [BMIm][BETI] accurately scales for density changes of 17%; however, there is a departure from scaling for TPIL for even more modest variations of temperature and pressure. The entropy of both ionic samples conforms to density scaling only if the scaling exponent is allowed to vary linearly with the magnitude of the entropy.

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.

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.

Thermodynamics of freezing and melting

Although the freezing of liquids and melting of crystals are fundamental for many areas of the sciences, even simple properties like the temperature–pressure relation along the melting line cannot be predicted today. Here we present a theory in which properties of the coexisting crystal and liquid phases at a single thermodynamic state point provide the basis for calculating the pressure, density and entropy of fusion as functions of temperature along the melting line, as well as the variation along this line of the reduced crystalline vibrational mean-square displacement (the Lindemann ratio), and the liquid's diffusion constant and viscosity. The framework developed, which applies for the sizable class of systems characterized by hidden scale invariance, is validated by computer simulations of the standard 12-6 Lennard-Jones system.

Melting is a classic first-order phase transition, but an accurate thermodynamic description is still lacking. Here, Pedersen et al. develop a theory, validated by simulations of the Lennard-Jones system, for the melting thermodynamics applicable to all systems characterized by hidden scale invariance.

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.

Thermodynamic consequences of the kinetic nature of the glass transition

In this paper, we consider the glass transition as a kinetic process and establish one universal equation for the pressure coefficient of the glass transition temperature, dTg/dp, which is a thermodynamic characteristic of this process. Our findings challenge the common previous expectations concerning key characteristics of the transformation from the liquid to the glassy state, because it suggests that without employing an additional condition, met in the glass transition, derivation of the two independent equations for dTg/dp is not possible. Hence, the relation among the thermodynamic coefficients, which could be equivalent to the well-known Prigogine-Defay ratio for the process under consideration, cannot be obtained. Besides, by comparing the predictions of our universal equation for dTg/dp and Ehrenfest equations, we find the aforementioned supplementary restriction, which must be met to use the Prigogine-Defay ratio for the glass transition.

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.