Role of Attractive Forces in the Relaxation Dynamics of Supercooled 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.

Heterogeneous attenuation of sound waves in three-dimensional amorphous solids

Sound waves are attenuated as they propagate in amorphous materials. We investigate the mechanism driving sound attenuation in the Rayleigh scattering regime by resolving the dynamics of an excited phonon in time and space via numerical simulations. We find sound attenuation is spatiotemporally heterogeneous. It starts in localized regions, which identify soft regions within the material and correlate with low-frequency vibrational modes. As time progresses, the regions where sound is primarily attenuated invade the system via an apparent diffusive process.

Comparing four hard-sphere approximations for the low-temperature WCA melting line

By combining interface-pinning simulations with numerical integration of the Clausius–Clapeyron equation, we accurately determine the melting-line coexistence pressure and fluid/crystal densities of the Weeks–Chandler–Andersen system, covering four decades of temperature. The data are used for comparing the melting-line predictions of the Boltzmann, Andersen–Weeks–Chandler, Barker–Henderson, and Stillinger hard-sphere approximations. The Andersen–Weeks–Chandler and Barker–Henderson theories give the most accurate predictions, and they both work excellently in the zero-temperature limit for which analytical expressions are derived here.

Intrinsic viscuit probability distribution functions for transport coefficients of liquids and solids

A reformulation of the Green–Kubo expressions for the transport coefficients of liquids in terms of a probability distribution function (PDF) of short trajectory contributions, which were named “viscuits,” has been explored in a number of recent publications. The viscuit PDF, P, is asymmetric on the two sides of the distribution. It is shown here using equilibrium 3D and 2D molecular dynamics simulations that the viscuit PDF of a range of simple molecular single component and mixture liquid and solid systems can be expressed in terms of the same intrinsic PDF ( P0), which is derived from P with the viscuit normalized by the standard deviation separately on each side of the distribution. P0 is symmetric between the two sides and can be represented for not very small viscuit values by the same gamma distribution formulated in terms of a single disposable parameter. P0 tends to an exponential in the large viscuit wings. Scattergrams of the viscuits and their associated single trajectory correlation functions are shown to distinguish effectively between liquids, solids, and glassy systems. The so-called viscuit square root method for obtaining the transport coefficients is shown to be a useful probe of small and statistically zero self-diffusion coefficients of molecules in the liquid and solid states, respectively. The results of this work suggest that the transport coefficients have a common underlying physical origin, reflecting at a coarse-grained level the traversal statistics of the system through its high-dimensioned potential energy landscape.

Finite-size scaling and thermodynamics of model supercooled liquids: long-range concentration fluctuations and the role of attractive interactions

We compute partial structure factors, Kirkwood-Buff integrals (KBIs) and chemical potentials of model supercooled liquids with and without attractive interactions. We aim at investigating whether relatively small differences in the tail of the radial distribution functions result in contrasting thermodynamic properties. Our results suggest that the attractive potential favours the nucleation of long-range structures. Indeed, upon decreasing temperature, Bathia-Thornton structure factors display anomalous behaviour in the k→0 limit. KBIs extrapolated to the thermodynamic limit confirm this picture, and excess coordination numbers identify the anomaly with long-range concentration fluctuations. By contrast, the purely repulsive system remains perfectly miscible for the same temperature interval and only reveals qualitatively similar concentration fluctuations in the crystalline state. Furthermore, differences in both isothermal compressibilities and chemical potentials show that thermodynamics is not entirely governed by the short-range repulsive part of the interaction potential, emphasising the nonperturbative role of attractive interactions. Finally, at higher density, where both systems display nearly identical dynamical properties and repulsive interactions become dominant, the anomaly disappears, and both systems also exhibit similar thermodynamic properties.

Unifying Description of the Vibrational Anomalies of Amorphous Materials

The vibrational density of states D(ω) of solids controls their thermal and transport properties. In crystals, the low-frequency modes are extended phonons distributed in frequency according to Debye's law, D(ω)∝ω^{2}. In amorphous solids, phonons are damped, and at low frequency D(ω) comprises extended modes in excess over Debye's prediction, leading to the so-called boson peak in D(ω)/ω^{2} at ω_{bp}, and quasilocalized ones. Here we show that boson peak and phonon attenuation in the Rayleigh scattering regime are related, as suggested by correlated fluctuating elasticity theory, and that amorphous materials can be described as homogeneous isotropic elastic media punctuated by quasilocalized modes acting as elastic heterogeneities. Our numerical results resolve the conflict between theoretical approaches attributing amorphous solids' vibrational anomalies to elastic disorder and localized defects.

How attractive and repulsive interactions affect structure ordering and dynamics of glass-forming liquids

The theory developed in our previous papers [Phys. Rev. E 99, 030101(R) (2019)10.1103/PhysRevE.99.030101; Phys. Rev. E 103, 032611 (2021)10.1103/PhysRevE.103.032611] is applied in this paper to investigate the dependence of slowing down of dynamics of glass-forming liquids on the attractive and repulsive parts of intermolecular interactions. Through an extensive comparison of the behavior of a Lennard-Jones glass-forming liquid and that of its WCA reduction to a model with truncated pair potential without attractive tail, we demonstrate why the two systems exhibit very different dynamics despite having nearly identical pair correlation functions. In particular, we show that local structures characterized by the number of mobile and immobile particles around a central particle markedly differ in the two systems at densities and temperatures where their dynamics show large difference and nearly identical where dynamics nearly overlap. We also show how the parameter ψ(T) that measures the role of fluctuations embedded in the system on size of the cooperatively reorganizing cluster (CRC) and the crossover temperature T_{a} depend on the intermolecular interactions. These parameters stemming from the intermolecular interactions characterize the temperature and density dependence of structural relaxation time τ_{α}. The quantitative and qualitative agreements found with simulation results for the two systems suggest that our theory brings out the underlying features that determine the dynamics of glass-forming liquids.

Mechanical disorder of sticky-sphere glasses. II. Thermomechanical inannealability

Many structural glasses feature static and dynamic mechanical properties that can depend strongly on glass formation history. The degree of universality of this history dependence and what it is possibly affected by are largely unexplored. Here we show that the variability of elastic properties of simple computer glasses under thermal annealing depends strongly on the strength of attractive interactions between the glasses' constituent particles-referred to here as glass "stickiness." We find that in stickier glasses the stiffening of the shear modulus with thermal annealing is strongly suppressed, while the thermal-annealing-induced softening of the bulk modulus is enhanced. Our key finding is that the characteristic frequency and density per frequency of soft quasilocalized modes becomes effectively invariant to annealing in very sticky glasses; the latter are therefore deemed "thermomechanically inannealable." The implications of our findings and future research directions are discussed.

Mechanical disorder of sticky-sphere glasses. I. Effect of attractive interactions

Recent literature indicates that attractive interactions between particles of a dense liquid play a secondary role in determining its bulk mechanical properties. Here we show that, in contrast with their apparent unimportance to the bulk mechanics of dense liquids, attractive interactions can have a major effect on macro- and microscopic elastic properties of glassy solids. We study several broadly applicable dimensionless measures of stability and mechanical disorder in simple computer glasses, in which the relative strength of attractive interactions-referred to as "glass stickiness"-can be readily tuned. We show that increasing glass stickiness can result in the decrease of various quantifiers of mechanical disorder, on both macro- and microscopic scales, with a pair of intriguing exceptions to this rule. Interestingly, in some cases strong attractions can lead to a reduction of the number density of soft, quasilocalized modes, by up to an order of magnitude, and to a substantial decrease in their core size, similar to the effects of thermal annealing on elasticity observed in recent works. Contrary to the behavior of canonical glass models, we provide compelling evidence indicating that the stabilization mechanism in our sticky-sphere glasses stems predominantly from the self-organized depletion of interactions featuring large, negative stiffnesses. Finally, we establish a fundamental link between macroscopic and microscopic quantifiers of mechanical disorder, which we motivate via scaling arguments. Future research directions are discussed.

Excess-entropy scaling in supercooled binary mixtures

Transport coefficients, such as viscosity or diffusion coefficient, show significant dependence on density or temperature near the glass transition. Although several theories have been proposed for explaining this dynamical slowdown, the origin remains to date elusive. We apply here an excess-entropy scaling strategy using molecular dynamics computer simulations and find a quasiuniversal, almost composition-independent, relation for binary mixtures, extending eight orders of magnitude in viscosity or diffusion coefficient. Metallic alloys are also well captured by this relation. The excess-entropy scaling predicts a quasiuniversal breakdown of the Stokes-Einstein relation between viscosity and diffusion coefficient in the supercooled regime. Additionally, we find evidence that quasiuniversality extends beyond binary mixtures, and that the origin is difficult to explain using existing arguments for single-component quasiuniversality.

Role of Attractive Interactions in Structure Ordering and Dynamics of Glass-Forming Liquids

A key question in glass physics is what the origin of slow glassy dynamics is. The liquid structure is a natural candidate; however, an apparently severe counterexample has been known. Two model glass-forming liquids, with the standard Lennard-Jones interaction potential and its Weeks-Chandler-Andersen variation without the attractive tail, exhibit very similar structures at the two-body level but drastically different dynamical behaviors in the supercooled states. Here we look at the liquid structure through a (many-body) structural order parameter Θ characterizing the packing capability of local particle arrangements. We show that the structures of these two systems seen by Θ are actually very different at a many-body level, but, quite surprisingly, the macroscopic structure (Θ)-dynamics (τ_{α}) relationships commonly follow a Vogel-Fulcher-Tammann-like function. Furthermore, the mutual information analysis reveals strong local structure-dynamics correlations. Therefore, we conclude that attractive interactions affect the liquid structure in a nonperturbative manner, but a general structural origin of slow dynamics holds for these systems.