Static triplet correlations in glass-forming liquids: A molecular dynamics study

A comparative study of the correlation between the structure and the dynamics for systems interacting via attractive and repulsive potentials

We present the study of the structure-dynamics correlation for systems interacting via attractive Lennard-Jones (LJ) and its repulsive counterpart, the Weeks-Chandler-Andersen (WCA) potentials. The structural order parameter (SOP) is related to the microscopic mean-field caging potential. At a particle level, the SOP shows a distribution. Although the two systems have similar pair structures, their average SOP differs. However, this difference alone is insufficient to explain the well known slowing down of the dynamics in the LJ system at low temperatures. The slowing down can be explained in terms of a stronger coupling between the SOP and the dynamics. To understand the origin of this system specific coupling, we study the difference in the microscopic structure between the hard and soft particles. We find that for the LJ system, the structural differences of the hard and soft particles are more significant and have a much stronger temperature dependence compared to the WCA system. Thus, the study suggests that attractive interaction creates more structurally different communities. This broader difference in the structural communities is probably responsible for stronger coupling between the structure and dynamics. Thus, the system specific structure-dynamics correlation, which also leads to a faster slowing down in the dynamics, appears to have a structural origin. A comparison of the predictive power of our SOP with the local energy and two body excess entropy in determining the dynamics shows that in the LJ system, the enthalpy plays a dominant role and in the WCA system, the entropy plays a dominant role, and our SOP can capture both these contributions.

Emergent structural correlations in dense liquids

The complete quantitative description of the structure of dense and supercooled liquids remains a notoriously difficult problem in statistical physics. Most studies to date focus solely on two-body structural correlations, and only a handful of papers have sought to consider additional three-body correlations. Here, we go beyond the state of the art by extracting many-body static structure factors from molecular dynamics simulations and by deriving accurate approximations up to the six-body structure factor via density functional theory. We find that supercooling manifestly increases four-body correlations, akin to the two- and three-body case. However, at small wave numbers, we observe that the four-point structure of a liquid drastically changes upon supercooling, both qualitatively and quantitatively, which is not the case in two-point structural correlations. This indicates that theories of the structure or dynamics of dense liquids should incorporate many-body correlations beyond the two-particle level to fully capture their intricate behavior.

Dimensionality reduction of local structure in glassy binary mixtures

We consider unsupervised learning methods for characterizing the disordered microscopic structure of supercooled liquids and glasses. Specifically, we perform dimensionality reduction of smooth structural descriptors that describe radial and bond-orientational correlations and assess the ability of the method to grasp the essential structural features of glassy binary mixtures. In several cases, a few collective variables account for the bulk of the structural fluctuations within the first coordination shell and also display a clear connection with the fluctuations of particle mobility. Fine-grained descriptors that characterize the radial dependence of bond-orientational order better capture the structural fluctuations relevant for particle mobility but are also more difficult to parameterize and to interpret. We also find that principal component analysis of bond-orientational order parameters provides identical results to neural network autoencoders while having the advantage of being easily interpretable. Overall, our results indicate that glassy binary mixtures have a broad spectrum of structural features. In the temperature range we investigate, some mixtures display well-defined locally favored structures, which are reflected in bimodal distributions of the structural variables identified by dimensionality reduction.

Many-Body Correlations Are Non-negligible in Both Fragile and Strong Glassformers

It is widely believed that the emergence of slow glassy dynamics is encoded in a material's microstructure. First-principles theory [mode-coupling theory (MCT)] is able to predict the dramatic slowdown of the dynamics from only static two-point correlations as input, yet it cannot capture all of the observed dynamical behavior. Here we go beyond two-point spatial correlation functions by extending MCT systematically to include higher-order static and dynamic correlations. We demonstrate that only adding the static triplet direct correlations already qualitatively changes the predicted glass-transition diagram of binary hard spheres and silica. Moreover, we find a nontrivial competition between static triplet correlations that work to stabilize the glass state and dynamic higher-order correlations that destabilize it for both materials. We conclude that the conventionally neglected static triplet direct correlations as well as higher-order dynamic correlations are, in fact, non-negligible in both fragile and strong glassformers.

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.

Multi-component generalized mode-coupling theory: predicting dynamics from structure in glassy mixtures

The emergence of glassy dynamics and the glass transition in dense disordered systems is still not fully understood theoretically. Mode-coupling theory (MCT) has shown to be effective in describing some of the non-trivial features of glass formation, but it cannot explain the full glassy phenomenology due to the strong approximations on which it is based. Generalized mode-coupling theory (GMCT) is a hierarchical extension of the theory, which is able to outclass MCT by carefully describing the dynamics of higher-order correlations in its generalized framework. Unfortunately, the theory has so far only been developed for single-component systems and as a result works poorly for highly polydisperse materials. In this paper, we solve this problem by developing GMCT for multi-component systems. We use it to predict the glassy dynamics of the binary Kob-Andersen Lennard-Jones mixture, as well as its purely repulsive Weeks-Chandler-Andersen analogue. Our results show that each additional level of the GMCT hierarchy gradually improves the predictive power of GMCT beyond its previous limit. This implies that our theory is able to harvest more information from the static correlations, thus being able to better understand the role of attraction in supercooled liquids from a first-principles perspective.

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.

A structural study and its relation to dynamic heterogeneity in a polymer glass former

The relationship between structure and dynamical behavior (super-Arrhenius temperature dependence of relaxation time accompanied by heterogeneous dynamics) in glassy materials remains an open issue in the physics of condensed matter. The question of whether this dynamic phenomena have a thermodynamic origin or not still remains unanswered. In this work we analyze several dynamic and structural parameters in a polymer glass-former by means of molecular dynamics simulations. The results obtained in this work indicate that the structure does affect dynamic behavior, whereas structural conditioning becomes noticeable below the temperature at which the non-Arrhenius behavior manifests and increases as the system approaches the glass transition temperature. Moreover, we observed that the short-range order parameters are related to local dynamics at the single-particle level. These results reinforce the idea of a connection between the structure and dynamics and that could indicate the thermodynamic nature of glass transition.

Structure of the simple harmonic-repulsive system in liquid and glassy states studied by the triple correlation function

An efficient description of the structures of liquids and, in particular, the structural changes that happen with liquids on supercooling remains to be a challenge. The systems composed of soft particles are especially interesting in this context because they often demonstrate non-trivial local orders that do not allow to introduce the concept of the nearest-neighbor shell. For this reason, the use of some methods, developed for the structure analysis of atomic liquids, is questionable for the soft-particle systems. Here we report about our investigations of the structure of the simple harmonic-repulsive liquid in 3D using the triple correlation function (TCF), i.e., the method that does not rely on the nearest neighbor concept. The liquid is considered at reduced pressure (P = 1.8) at which it exhibits remarkable stability against crystallization on cooling. It is demonstrated that the TCF allows addressing the development of the orientational correlations in the structures that do not allow drawing definite conclusions from the studies of the bond-orientational order parameters. Our results demonstrate that the orientational correlations, if measured by the heights of the peaks in the TCF, significantly increase on cooling. This rise in the orientational ordering is not captured properly by the Kirkwood's superposition approximation. Detailed considerations of the peaks' shapes in the TCF suggest the existence of a link between the orientational ordering and the slowdown of the system's dynamics. Our findings support the view that the development of the orientational correlations in liquids may play a significant role in the liquids' dynamics and that the considerations of the pair distribution function may not be sufficient to understand intuitively all the structural changes that happen with liquids on supercooling. In general, our results demonstrate that the considerations of the TCF are useful in the discussions of the liquid's structures beyond the pair density function and interpreting the results obtained with the bond-orientational order parameters.

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.

Attractive versus truncated repulsive supercooled liquids: The dynamics is encoded in the pair correlation function

We compare glassy dynamics in two liquids that differ in the form of their interaction potentials. Both systems have the same repulsive interactions but one has also an attractive part in the potential. These two systems exhibit very different dynamics despite having nearly identical pair correlation functions. We demonstrate that a properly weighted integral of the pair correlation function, which amplifies the subtle differences between the two systems, correctly captures their dynamical differences. The weights are obtained from a standard machine learning algorithm.

Role of Attractive Forces in the Relaxation Dynamics of Supercooled Liquids

The attractive tail of the intermolecular interaction affects very weakly the structural properties of liquids, while it affects dramatically their dynamical ones. Via the numerical simulations of model systems not prone to crystallization, both in three and in two spatial dimensions, here we demonstrate that the nonperturbative dynamical effects of the attractive forces are tantamount to a rescaling of the activation energy by the glass transition temperature T_{g}: systems only differing in their attractive interaction have the same structural and dynamical properties if compared at the same value of T/T_{g}.

Understanding, predicting, and tuning the fragility of vitrimeric polymers

Fragility is an empirical property that describes how abruptly a glass-forming material solidifies upon supercooling. The degree of fragility carries important implications for the functionality and processability of a material, as well as for our fundamental understanding of the glass transition. However, the microstructural properties underlying fragility still remain poorly understood. Here, we explain the microstructure-fragility link in vitrimeric networks, a novel type of high-performance polymers with unique bond-swapping functionality and unusual glass-forming behavior. Our results are gained from coarse-grained computer simulations and first-principles mode-coupling theory (MCT) of star-polymer vitrimers. We first demonstrate that the vitrimer fragility can be tuned over an unprecedentedly broad range, from fragile to strong and even superstrong behavior, by decreasing the bulk density. Remarkably, this entire phenomenology can be reproduced by microscopic MCT, thus challenging the conventional belief that existing first-principles theories cannot account for nonfragile behaviors. Our MCT analysis allows us to rationally identify the microstructural origin of the fragile-to-superstrong crossover, which is rooted in the sensitivity of the static structure factor to temperature variations. On the molecular scale, this behavior stems from a change in dominant length scales, switching from repulsive excluded-volume interactions to intrachain attractions as the vitrimer density decreases. Finally, we develop a simplified schematic MCT model which corroborates our microscopically founded conclusions and which unites our findings with earlier MCT studies. Our work sheds additional light on the elusive structure-fragility link in glass-forming matter and provides a first-principles-based platform for designing amorphous materials with an on-demand dynamic response.

Computation of static quantum triplet structure factors of liquid para-hydrogen

The instantaneous and centroid triplet structure factors, S ( 3 ) ( k 1 , k 2 ) , of liquid (one-center) para-hydrogen are computed on the crystallization line for temperatures T/K ≤ 33. The focus is on salient equilateral and isosceles features, and the methods utilized are path integral Monte Carlo (PIMC) simulations and Ornstein-Zernike (OZ) integral equations, which involve Jackson-Feenberg convolution (JF3) and other distinct closures. Long path integral simulation runs are carried out in the canonical ensemble, so as to obtain sufficiently accurate direct PI triplet results. Conclusions are drawn regarding general triplet structure features and the role and usefulness of the OZ closures. The equilateral features are studied in more detail, and one finds that (a) PIMC results point to the existence of regularity in the centroid main peak amplitudes; (b) some of the studied closures give qualitative descriptions for wave numbers below k ≈ 1 Å^-1, but they all fail to describe the main peak amplitude regions (1.75 < k/Å^-1 < 2.5); and (c) JF3 plays the role of a limit closure that is valid for increasing wave numbers (k ≥ 2.6 Å^-1). In addition, representative isosceles PI features turn out to be reasonably bounded (within Δk = 0.1 Å^-1) by those of some closures.

Local structure in deeply supercooled liquids exhibits growing lengthscales and dynamical correlations

Glasses are among the most widely used of everyday materials, yet the process by which a liquid's viscosity increases by 14 decades to become a glass remains unclear, as often contradictory theories provide equally good descriptions of the available data. Knowledge of emergent lengthscales and higher-order structure could help resolve this, but this requires time-resolved measurements of dense particle coordinates-previously only obtained over a limited time interval. Here we present an experimental study of a model colloidal system over a dynamic window significantly larger than previous measurements, revealing structural ordering more strongly linked to dynamics than previously found. Furthermore we find that immobile regions and domains of local structure grow concurrently with density, and that these regions have low configurational entropy. We thus show that local structure plays an important role at deep supercooling, consistent with a thermodynamic interpretation of the glass transition rather than a principally dynamic description.

On static triplet structures in fluids with quantum behavior

The problem of the equilibrium triplet structures in fluids with quantum behavior is discussed. Theoretical questions of interest to the real space structures are addressed by studying the three types of structures that can be determined via path integrals (instantaneous, centroid, and total thermalized-continuous linear response). The cases of liquid para-H₂ and liquid neon on their crystallization lines are examined with path-integral Monte Carlo simulations, the focus being on the instantaneous and the centroid triplet functions (equilateral and isosceles configurations). To analyze the results further, two standard closures, Kirkwood superposition and Jackson-Feenberg convolution, are utilized. In addition, some pilot calculations with path integrals and closures of the instantaneous triplet structure factor of liquid para-H₂ are also carried out for the equilateral components. Triplet structural regularities connected to the pair radial structures are identified, a remarkable usefulness of the closures employed is observed (e.g., triplet spatial functions for medium-long distances, triplet structure factors for medium k wave numbers), and physical insight into the role of pair correlations near quantum crystallization is gained.

Unraveling the success and failure of mode coupling theory from consideration of entropy

We analyze the dynamics of model supercooled liquids in a temperature regime where predictions of mode coupling theory (MCT) are known to be valid qualitatively. In this regime, the Adam-Gibbs (AG) relation, based on an activation picture of dynamics, also describes the dynamics satisfactorily, and we explore the mutual consistency and interrelation of these descriptions. Although entropy and dynamics are related via phenomenological theories, the connection between MCT and entropy has not been argued for. In this work, we explore this connection and provide a microscopic derivation of the phenomenological Rosenfeld theory. At low temperatures, the overlap between the MCT power law regime and AG relation implies that the AG relation predicts an avoided divergence at Tc, the origin of which can be related to the vanishing of pair configurational entropy, which we find occurring at the same temperature. We also show that the residual multiparticle entropy plays an important role in describing the relaxation time.

Microscopic Theory for the Role of Attractive Forces in the Dynamics of Supercooled Liquids

We formulate a microscopic, no adjustable parameter, theory of activated relaxation in supercooled liquids directly in terms of the repulsive and attractive forces within the framework of pair correlations. Under isochoric conditions, attractive forces can nonperturbatively modify slow dynamics, but at high enough density their influence vanishes. Under isobaric conditions, attractive forces play a minor role. High temperature apparent Arrhenius behavior and density-temperature scaling are predicted. Our results are consistent with recent isochoric simulations and isobaric experiments on a deeply supercooled molecular liquid. The approach can be generalized to treat colloidal gelation and glass melting, and other soft matter slow dynamics problems.

Brownian dynamics: from glassy to trivial

We endow a system of interacting particles with two distinct, local, Markovian, and reversible microscopic dynamics that both converge to the Boltzmann-Gibbs equilibrium of standard liquids. While the first, standard, one leads to glassy dynamics, we use field-theoretical techniques to show that the latter displays no sign of glassiness. The approximations we use, akin to the mode-coupling approximation, are famous for magnifying glassy aspects of the dynamics, supposedly through the neglect of activated events. Despite this, the modified dynamics seem to stick to standard liquid relaxation. This finding singles out as applying to a realistic system of interacting particles in low dimensions and questions the role of the dynamical rules used to explore a given static free-energy landscape. Moreover, our peculiar choice of dynamical rules offers the possibility of a direct connection with replica theory, and our findings therefore call for a clarification of the interplay between replica theory and the underlying dynamics of the system.

Role of Structure and Entropy in Determining Differences in Dynamics for Glass Formers with Different Interaction Potentials

We present a study of two model liquids with different interaction potentials, exhibiting similar structure but significantly different dynamics at low temperatures. By evaluating the configurational entropy, we show that the differences in the dynamics of these systems can be understood in terms of their thermodynamic differences. Analyzing their structure, we demonstrate that differences in pair correlation functions between the two systems, through their contribution to the entropy, dominate the differences in their dynamics, and indeed overestimate the differences. Including the contribution of higher order structural correlations to the entropy leads to smaller estimates for the relaxation times, as well as smaller differences between the two studied systems.

Structural evolution in the aging process of supercooled colloidal liquids

When a liquid is rapidly quenched to a temperature below the glass-transition point, it is driven out of equilibrium; it then slowly relaxes to a (quasi)equilibrium state. This slow relaxation process is called aging. By definition, any glasses are inevitably in the process of aging and actually slowly evolving with time. Thus the study of aging phenomena is of fundamental importance for understanding not only the nonequilibrium nature of the glass transition, but also the stability of glassy materials. Here we consider aging after a rather shallow quench, for which a system is still able to reach (metastable) equilibrium. By using polydisperse colloidal liquids as a model, we show the validity of dynamical scaling that there is only one relevant length scale not only for a quasiequilibrium supercooled state but also for a nonequilibrium process of aging, which is reminiscent of dynamical critical phenomena. Our finding indicates that the aging toward (metastable) equilibrium may be regarded as the growth process of critical-like fluctuations of static order associated with low-free-energy configurations, further suggesting that this ordering is the origin of cooperative slow dynamics in the systems studied. The generality of this statement for other glass-forming systems remains for a future study.