Mean-field theory, mode-coupling theory, and the onset temperature in supercooled liquids

Does the Naı̈ve Mode-Coupling Power Law Divergence Provide an Objective Determination of the Crossover Temperature in Glass Formation Behavior?

The glass formation temperature range is commonly divided into a weakly supercooled regime at higher temperatures and a deeply supercooled regime at lower temperatures, with a change in the physical mechanisms that govern dynamics often postulated to occur at the crossover between these regimes. This crossover temperature Tc is widely determined based on a fit of relaxation time vs temperature data to a power law divergence form predicted by the naı̈ve mode coupling theory (MCT). Here, we show, based on simulation data spanning polymeric, small molecule organic, metallic, and inorganic glass formers, that this approach does not yield an objective measure of a crossover temperature. Instead, the value of Tc is determined by the lowest temperature Tmin employed in the fit, and no regime of stationary or convergent Tc value is generally observed as Tmin is varied. Nor does the coefficient of determination R2 provide any robust means of selecting a fit range and thus a value of Tc. These results may require a re-evaluation of published results that have employed the fit MCT Tc value as a metric of temperature-dependent dynamics or a benchmark for depth of supercooling, and they highlight a need for the field to converge on a more objective determination of any posited crossover temperature between high and low temperature regimes of glass formation.

Freezing, melting, and the onset of glassiness in binary mixtures

We clarify the relationship between freezing, melting, and the onset of glassy dynamics in a prototypical glass-forming mixture model. Our starting point is a precise operational definition of the onset of glassiness, as expressed by the emergence of inflections in time-dependent correlation functions. By scanning the temperature-composition phase diagram of the mixture, we find a disconnect between the onset of glassiness and freezing. Surprisingly, however, the onset temperature closely tracks the melting line, along which the excess entropy is approximately constant. At fixed composition, all characteristic temperatures display nonetheless similar pressure dependencies, which are very well predicted by the isomorph theory. While our results rule out a general connection between thermodynamic metastability and glassiness, they call for a reassessment of the role of crystalline precursors in glass-forming liquids.

Jamming, relaxation, and memory in a minimally structured glass former

Structural glasses form through various out-of-equilibrium processes, including temperature quenches, rapid compression (crunches), and shear. Although each of these processes should be formally understandable within the recently formulated dynamical mean-field theory (DMFT) of glasses, the numerical tools needed to solve the DMFT equations up to the relevant physical regime do not yet exist. In this context, numerical simulations of minimally structured (and therefore mean-field-like) model glass formers can aid the search for and understanding of such solutions, thanks to their ability to disentangle structural from dimensional effects. We study here the infinite-range Mari-Kurchan model under simple out-of-equilibrium processes, and we compare results with the random Lorentz gas [J. Phys. A 55, 334001 (2022)10.1088/1751-8121/ac7f06]. Because both models are mean-field-like and formally equivalent in the limit of infinite spatial dimensions, robust features are expected to appear in the DMFT as well. The comparison provides insight into temperature and density onsets, memory, as well as anomalous relaxation. This work also further enriches the algorithmic understanding of the jamming density.

Unifying fluctuation-dissipation temperatures of slow-evolving nonequilibrium systems from the perspective of inherent structures

For nonequilibrium systems, how to define temperature is one of the key and difficult issues to solve. Although effective temperatures have been proposed and studied to this end, it still remains elusive what they actually are. Here, we focus on the fluctuation-dissipation temperatures and report that such effective temperatures of slow-evolving systems represent characteristic temperatures of their equilibrium counterparts. By calculating the fluctuation-dissipation relation of inherent structures, we obtain a temperature-like quantity T IS For monocomponent crystal-formers, T IS agrees well with the crystallization temperature T c, while it matches with the onset temperature T on for glass-formers. It also agrees with effective temperatures of typical nonequilibrium systems, such as aging glasses, quasi-static shear flows, and quasi-static self-propelled flows. From the unique perspective of inherent structures, our study reveals the nature of effective temperatures and the underlying connections between nonequilibrium and equilibrium systems and confirms the equivalence between T on and T c.

Memory Formation in Jammed Hard Spheres

Liquids equilibrated below an onset condition share similar inherent states, while those above that onset have inherent states that markedly differ. Although this type of materials memory was first reported in simulations over 20 years ago, its physical origin remains controversial. Its absence from mean-field descriptions, in particular, has long cast doubt on its thermodynamic relevance. Motivated by a recent theoretical proposal, we reassess the onset phenomenology in simulations using a fast hard sphere jamming algorithm and find it to be both thermodynamically and dimensionally robust. Remarkably, we also uncover a second type of memory associated with a Gardner-like regime of the jamming algorithm.

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.

An energy-landscape-based crossover temperature in glass-forming liquids

The systematic identification of temperature scales in supercooled liquids that are key to understanding those liquids' underlying glass properties, and their formation-history dependence, is a challenging task. Here, we study the statistics of particles' squared displacements δr² between equilibrium liquid configurations at temperature T and their underlying inherent states, using computer simulations of 11 different computer glass formers. We show that the relative fluctuations of δr² are nonmonotonic in T, exhibiting a maximum whose location defines the crossover temperature T_X. Therefore, T_X marks the point of maximal heterogeneity during the process of tumbling down the energy landscape, starting from an equilibrium liquid state at temperature T down to its underlying inherent state. We extract T_X for the 11 employed computer glasses, ranging from tetrahedral glasses to packs of soft elastic spheres, and demonstrate its usefulness in putting the elastic properties of different glasses on the same footing. Interestingly, we further show that T_X marks the crossover between two distinct regimes of the mean ⟨δr²⟩: a high temperature regime in which ⟨δr²⟩ scales approximately as T^0.5 and a deeply supercooled regime in which ⟨δr²⟩ scales approximately as T^1.3. Further research directions are discussed.

The polarizability response of a glass-forming liquid reveals intrabasin motion and interbasin transitions on a potential energy landscape

Potential energy landscape (PEL) concepts have been useful in conceptualizing the effects of intermolecular interactions on dynamic and thermodynamic properties of liquids and glasses. "Basins", or regions of reduced potential energy associated with locally preferred molecular packing are important PEL features. The molecular configurations at the bottom of these basins are referred to as inherent structures (ISs). Experimental methods for directly characterizing PEL features such as these are rare, largely relegating PEL concepts to theory and simulation studies, and impeding their exploration in real systems. Recently, we showed that quasielastic neutron scattering (QENS) data from propylene carbonate (PC) exhibit signatures of picosecond timescale motion that are consistent with intrabasin motion and interbasin transitions [Cicerone et al., J. Chem. Phys., 2017, 146, 054502]. Here we present optically-heterodyne-detected optical Kerr effect (OHD-OKE) spectroscopy studies on PC. The data exhibit signatures of motion within and transitions between basins that agree quantitatively with and extend the QENS results. We show that the librational component of the OKE response corresponds to intrabasin dynamics, and the enigmatic intermediate OKE response corresponds to interbasin transition events. The OKE data extend the measurement range of these parameters and reveal their utility in characterizing PEL features of real systems.

Connecting glass-forming ability of binary mixtures of soft particles to equilibrium melting temperatures

The glass-forming ability is an important material property for manufacturing glasses and understanding the long-standing glass transition problem. Because of the nonequilibrium nature, it is difficult to develop the theory for it. Here we report that the glass-forming ability of binary mixtures of soft particles is related to the equilibrium melting temperatures. Due to the distinction in particle size or stiffness, the two components in a mixture effectively feel different melting temperatures, leading to a melting temperature gap. By varying the particle size, stiffness, and composition over a wide range of pressures, we establish a comprehensive picture for the glass-forming ability, based on our finding of the direct link between the glass-forming ability and the melting temperature gap. Our study reveals and explains the pressure and interaction dependence of the glass-forming ability of model glass-formers, and suggests strategies to optimize the glass-forming ability via the manipulation of particle interactions.

Glass-forming ability is an important parameter for manufacturing glassy materials, but it remains challenging to be characterized due to its nonequilibrium nature. Nie et al. provide a solution by linking it to the pressure dependence of melting temperature of constituent components in binary mixtures.

Dynamic Heterogeneities in Colloidal Supercooled Liquids: Experimental Tests of Inhomogeneous Mode Coupling Theory

The dynamics in supercooled liquids slow enormously upon approaching the glass transition, albeit without significant change of liquid structure. This empirical observation has stimulated development of many theoretical models which attempt to elucidate microscopic mechanisms in glasses and glass precursors. Here, quasi-two-dimensional colloidal supercooled liquids and glasses are employed to experimentally test predictions of widely used models: mode coupling theory (MCT) and its important extension, inhomogeneous MCT (IMCT). We measure two-point dynamic correlation functions in the glass forming liquids to determine structural relaxation times, τ_α, and mode coupling exponents, a, b, and γ; these parameters are then used to extract the mode coupling dynamic crossover packing area-fraction, ϕ _c. This information, along with our measurements of supercooled liquid spatiotemporal dynamics, permits characterization of dynamic heterogeneities in the samples and facilitates direct experimental tests of the scaling predictions of IMCT. The time scales at which dynamic heterogeneities are largest, and their spatial sizes, exhibit power law growth on approaching ϕ _c. Within experimental error, the exponents of the measured power laws are close to the predictions of IMCT.

Determination of onset temperature from the entropy for fragile to strong liquids

In this paper, we establish a connection between the onset temperature of glassy dynamics with the change in the entropy for a wide range of model systems. We identify the crossing temperature of pair and excess entropies as the onset temperature. Below the onset temperature, the residual multiparticle entropy, the difference between excess and pair entropies, becomes positive. The positive entropy can be viewed as equivalent to the larger phase space exploration of the system. The new method of onset temperature prediction from entropy is less ambiguous, as it does not depend on any fitting parameter like the existing methods.

Role of the Pair Correlation Function in the Dynamical Transition Predicted by Mode Coupling Theory

In a recent study, we have found that for a large number of systems the configurational entropy at the pair level S_{c2}, which is primarily determined by the pair correlation function, vanishes at the dynamical transition temperature T_{c}. Thus, it appears that the information of the transition temperature is embedded in the structure of the liquid. In order to investigate this, we describe the dynamics of the system at the mean field level and, using the concepts of the dynamical density functional theory, show that the dynamical transition temperature depends only on the pair correlation function. Thus, this theory is similar in spirit to the microscopic mode coupling theory (MCT). However, unlike microscopic MCT, which predicts a very high transition temperature, the present theory predicts a transition temperature that is similar to T_{c}. This implies that the information of the dynamical transition temperature is embedded in the pair correlation function.

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.

Phase Transition for Quenched Coupled Replicas in a Plaquette Spin Model of Glasses

We study a three-dimensional plaquette spin model whose low temperature dynamics is glassy, due to localized defects and effective kinetic constraints. The thermodynamics of this system is smooth at all temperatures. We show that coupling it to a second system with a fixed (quenched) configuration leads to a phase transition, at finite coupling. The order parameter is the overlap between the copies, and the transition is between phases of low and high overlap. We find critical points whose properties are consistent with random-field Ising universality. We analyze the interfacial free energy cost between the high- and low-overlap states that coexist at (and below) the critical point, and we use this cost as the basis for a finite-size scaling analysis. We discuss these results in the context of mean-field and dynamical facilitation theories of the glass transition.

Equilibrium phase diagram of a randomly pinned glass-former

We use computer simulations to study the thermodynamic properties of a glass-former in which a fraction c of the particles has been permanently frozen. By thermodynamic integration, we determine the Kauzmann, or ideal glass transition, temperature [Formula: see text] at which the configurational entropy vanishes. This is done without resorting to any kind of extrapolation, i.e., [Formula: see text] is indeed an equilibrium property of the system. We also measure the distribution function of the overlap, i.e., the order parameter that signals the glass state. We find that the transition line obtained from the overlap coincides with that obtained from the thermodynamic integration, thus showing that the two approaches give the same transition line. Finally, we determine the geometrical properties of the potential energy landscape, notably the T- and c dependence of the saddle index, and use these properties to obtain the dynamic transition temperature [Formula: see text]. The two temperatures [Formula: see text] and [Formula: see text] cross at a finite value of c and indicate the point at which the glass transition line ends. These findings are qualitatively consistent with the scenario proposed by the random first-order transition theory.

Potential energy landscape and inherent dynamics of a hard-sphere fluid

Hard-sphere models exhibit many of the same kinds of supercooled-liquid behavior as more realistic models of liquids, but the highly nonanalytic character of their potentials makes it a challenge to think of that behavior in potential energy landscape terms. We show here that it is possible to calculate an important topological property of hard-sphere landscapes, the geodesic pathways through those landscapes, and to do so without artificially coarse-graining or softening the potential. We show, moreover, that the rapid growth of the lengths of those pathways with increasing packing fraction quantitatively predicts the precipitous decline in diffusion constants in a glass-forming hard-sphere mixture model. The geodesic paths themselves can be considered as defining the intrinsic dynamics of hard spheres, so it is also revealing to find that they (and therefore the features of the underlying potential energy landscape) correctly predict the occurrence of dynamic heterogeneity and nonzero values of the non-Gaussian parameter. The success of these landscape predictions for the dynamics of such a singular model emphasizes that there is more to potential energy landscapes than is revealed by looking at the minima and saddle points.

Coherent neutron scattering and collective dynamics on mesoscale

By combining, and modestly extending, a variety of theoretical concepts for the dynamics of liquids in the supercooled regime, we formulate a simple analytic model for the temperature and wavevector dependent collective density fluctuation relaxation time that is measurable using coherent dynamic neutron scattering. Comparison with experiments on the ionic glass-forming liquid Ca-K-NO3 in the lightly supercooled regime suggests the model captures the key physics in both the local cage and mesoscopic regimes, including the unusual wavevector dependence of the collective structural relaxation time. The model is consistent with the idea that the decoupling between diffusion and viscosity is reflected in a different temperature dependence of the collective relaxation time at intermediate wavevectors and near the main (cage) peak of the static structure factor. More generally, our analysis provides support for the ideas that decoupling information and growing dynamic length scales can be at least qualitatively deduced by analyzing the collective relaxation time as a function of temperature and wavevector, and that there is a strong link between dynamic heterogeneity phenomena at the single and many particle level. Though very simple, the model can be applied to other systems, such as molecular liquids.

Jamming transition and inherent structures of hard spheres and disks

Recent studies show that volume fractions φ(J) at the jamming transition of frictionless hard spheres and disks are not uniquely determined but exist over a continuous range. Motivated by this observation, we numerically investigate the dependence of φ(J) on the initial configurations of the parent fluid equilibrated at a volume fraction φ(eq), before compressing to generate a jammed packing. We find that φ(J) remains constant when φ(eq) is small but sharply increases as φ(eq) exceeds the dynamic transition point which the mode-coupling theory predicts. We carefully analyze configurational properties of both jammed packings and parent fluids and find that, while all jammed packings remain isostatic, the increase of φ(J) is accompanied with subtle but distinct changes of local orders, a static length scale, and an exponent of the finite-size scaling. These results are consistent with the scenario of the random first-order transition theory of the glass transition.

The role of attractive forces in viscous liquids

We present evidence from computer simulation that the slowdown of relaxation of a standard Lennard-Jones glass-forming liquid and that of its reduction to a model with truncated pair potentials without attractive tails are quantitatively and qualitatively different in the viscous regime. The pair structure of the two models is however very similar. This finding, which appears to contradict the common view that the physics of dense liquids is dominated by the steep repulsive forces between atoms, is characterized in detail, and its consequences are explored. Beyond the role of attractive forces themselves, a key aspect in explaining the differences in the dynamical behavior of the two models is the truncation of the interaction potentials beyond a cutoff at typical interatomic distance. This leads us to question the ability of the jamming scenario to describe the physics of glass-forming liquids and polymers.

Locally preferred structures and many-body static correlations in viscous liquids

The influence of static correlations beyond the pair level on the dynamics of selected model glass formers is investigated. The pair structure, angular distribution functions, and statistics of Voronoi polyhedra of two well-known Lennard-Jones mixtures as well as of the corresponding Weeks-Chandler-Andersen variants, in which the attractive part of the potential is truncated, are compared. By means of the Voronoi construction, the atomic arrangements corresponding to the locally preferred structures of the models are identified. It is found that the growth of domains formed by interconnected locally preferred structures signals the onset of the slow-dynamics regime and allows the rationalization of the different dynamic behaviors of the models. At low temperature, the spatial extension of the structurally correlated domains, evaluated at fixed relaxation time, increases with the fragility of the models and is systematically reduced by truncating the attractions. In view of these results, proper inclusion of many-body static correlations in theories of the glass transition appears crucial for the description of the dynamics of fragile glass formers.

Free-energy functional for freezing transitions: hard-sphere systems freezing into crystalline and amorphous structures

A free-energy functional that contains both the symmetry-conserved and symmetry-broken parts of the direct pair correlation function has been used to investigate the freezing of a system of hard spheres into crystalline and amorphous structures. The freezing parameters for fluid-crystal transition have been found to be in very good agreement with the results found from simulations. We considered amorphous structures found from molecular dynamics simulations at packing fractions η lower than the glass close packing fraction η(J) and investigated their stability compared to that of a homogeneous fluid. The existence of a free-energy minimum corresponding to a density distribution of overlapping Gaussians centered around an amorphous lattice depicts a deeply supercooled state with a heterogeneous density profile.

Activated dynamics in dense fluids of attractive nonspherical particles. I. Kinetic crossover, dynamic free energies, and the physical nature of glasses and gels

We apply the center-of-mass versions of naïve mode coupling theory and nonlinear Langevin equation theory to study how short-range attractive interactions modify the onset of localization, activated single-particle dynamics, and the physical nature of the transiently arrested state of a variety of dense nonspherical particle fluids (and the spherical analog) as a function of volume fraction and attraction strength. The form of the dynamic crossover boundary depends on particle shape, but the reentrant glass-fluid-gel phenomenon and the repulsive glass-to-attractive glass crossover always occur. Diverse functional forms of the dynamic free energy are found for all shapes including glasslike, gel-like, a glass-gel form defined by the coexistence of two localization minima and two activation barriers, and a "mixed" attractive glass characterized by a single, very short localization length but an activation barrier located at a large displacement as in repulsive-force caged glasses. For the latter state, particle trajectories are expected to be of a two-step activated form and can be accessed at high attraction strength by increasing volume fraction, or by increasing attraction strength at fixed high enough volume fraction. A new classification scheme for slow dynamics of fluids of dense attractive particles is proposed based on specification of both the nature of the localized state and the particle displacements required to restore ergodicity via activated barrier hopping. The proposed physical picture appears to be in qualitative agreement with recent computer simulations and colloid experiments.

Glass formation and shear elasticity in dense suspensions of repulsive anisotropic particles

Kinetic vitrification, shear elasticity, and the approach to jamming are investigated for repulsive nonspherical colloids and contrasted with their spherical analog. Particle anisotropy dramatically increases the volume fraction for kinetic arrest. The shear modulus of all systems increases roughly exponentially with volume fraction, and a universal collapse is achieved based on either the dynamic crossover or random close packing volume fraction as the key nondimensionalizing quantity. Quantitative comparisons with recent microscopic theories are performed and good agreement demonstrated.

Jamming transitions in amorphous packings of frictionless spheres occur over a continuous range of volume fractions

We numerically produce fully amorphous assemblies of frictionless spheres in three dimensions and study the jamming transition these packings undergo at large volume fractions. We specify four protocols yielding a critical value for the jamming volume fraction which is sharply defined in the limit of large system size, but is different for each protocol. Thus, we directly establish the existence of a continuous range of volume fractions where nonequilibrium jamming transitions occur. However, these jamming transitions share the same critical behavior. Our results suggest that, even in the absence of partial crystalline ordering, a unique location of a random close packing does not exist.

Numerical and theoretical study of a monodisperse hard-sphere glass former

There exists a variety of theories of the glass transition and many more numerical models. But because the models need built-in complexity to prevent crystallization, comparisons with theory can be difficult. We study the dynamics of a deeply supersaturated monodisperse four-dimensional (4D) hard-sphere fluid, which has no such complexity, but whose strong intrinsic geometrical frustration inhibits crystallization, even when deeply supersaturated. As an application, we compare its behavior to the mode-coupling theory (MCT) of glass formation. We find MCT to describe this system better than any other structural glass formers in lower dimensions. The reduction in dynamical heterogeneity in 4D suggested by a milder violation of the Stokes-Einstein relation could explain the agreement. These results are consistent with a mean-field scenario of the glass transition.

Soft colloids make strong glasses

Glass formation in colloidal suspensions has many of the hallmarks of glass formation in molecular materials. For hard-sphere colloids, which interact only as a result of excluded volume, phase behaviour is controlled by volume fraction, phi; an increase in phi drives the system towards its glassy state, analogously to a decrease in temperature, T, in molecular systems. When phi increases above phi* approximately 0.53, the viscosity starts to increase significantly, and the system eventually moves out of equilibrium at the glass transition, phi(g) approximately 0.58, where particle crowding greatly restricts structural relaxation. The large particle size makes it possible to study both structure and dynamics with light scattering and imaging; colloidal suspensions have therefore provided considerable insight into the glass transition. However, hard-sphere colloidal suspensions do not exhibit the same diversity of behaviour as molecular glasses. This is highlighted by the wide variation in behaviour observed for the viscosity or structural relaxation time, tau(alpha), when the glassy state is approached in supercooled molecular liquids. This variation is characterized by the unifying concept of fragility, which has spurred the search for a 'universal' description of dynamic arrest in glass-forming liquids. For 'fragile' liquids, tau(alpha) is highly sensitive to changes in T, whereas non-fragile, or 'strong', liquids show a much lower T sensitivity. In contrast, hard-sphere colloidal suspensions are restricted to fragile behaviour, as determined by their phi dependence, ultimately limiting their utility in the study of the glass transition. Here we show that deformable colloidal particles, when studied through their concentration dependence at fixed temperature, do exhibit the same variation in fragility as that observed in the T dependence of molecular liquids at fixed volume. Their fragility is dictated by elastic properties on the scale of individual colloidal particles. Furthermore, we find an equivalent effect in molecular systems, where elasticity directly reflects fragility. Colloidal suspensions may thus provide new insight into glass formation in molecular systems.

Glass transition of dense fluids of hard and compressible spheres

We use computer simulations to study the glass transition of dense fluids made of polydisperse repulsive spheres. For hard particles, we vary the volume fraction, phi , and use compressible particles to explore finite temperatures, T>0 . In the hard sphere limit, our dynamic data show evidence of an avoided mode-coupling singularity near phi(MCT) is approximately 0.592; they are consistent with a divergence of equilibrium relaxation times occurring at phi(0) is approximately 0.635, but they leave open the existence of a finite temperature singularity for compressible spheres at volume fraction phi>phi(0). Using direct measurements and a scaling procedure, we estimate the equilibrium equation of state for the hard sphere metastable fluid up to phi(0), where pressure remains finite, suggesting that phi(0) corresponds to an ideal glass transition. We use nonequilibrium protocols to explore glassy states above phi(0) and establish the existence of multiple equations of state for the unequilibrated glass of hard spheres, all diverging at different densities in the range phi in [0.642, 0.664]. Glassiness thus results in the existence of a continuum of densities where jamming transitions can occur.

Theory of coupled translational-rotational glassy dynamics in dense fluids of uniaxial particles

The naïve mode coupling theory (NMCT) for ideal kinetic arrest and the nonlinear Langevin equation theory of activated single-particle barrier hopping dynamics are generalized to treat the coupled center-of-mass (CM) translational and rotational motions of uniaxial hard objects in the glassy fluid regime. The key dynamical variables are the time-dependent displacements of the particle center-of-mass and orientational angle. The NMCT predicts a kinetic arrest diagram with three dynamical states: ergodic fluid, plastic glass, and fully nonergodic double glass, the boundaries of which meet at a "triple point" corresponding to a most difficult to vitrify diatomic of aspect ratio approximately 1.43. The relative roles of rotation and translation in determining ideal kinetic arrest are explored by examining three limits of the theory corresponding to nonrotating, pure rotation, and rotationally ergodic models. The ideal kinetic arrest boundaries represent a crossover to activated dynamics described by two coupled stochastic nonlinear Langevin equations for translational and rotational motions. The fundamental quantity is a dynamic free-energy surface, which for small aspect ratios in the high-volume fraction regime exhibits two saddle points reflecting a two-step activated dynamics where relatively rapid rotational dynamics coexists with slower CM translational motions. For large-enough aspect ratios, the dynamic free-energy surface has one saddle point which corresponds to a system-specific coordinated translation-rotation motion. The entropic barriers as a function of the relative amount of rotation versus translation are determined.

Growing length and time scales in glass-forming liquids

The glass transition, whereby liquids transform into amorphous solids at low temperatures, is a subject of intense research despite decades of investigation. Explaining the enormous increase in relaxation times of a liquid upon supercooling is essential for understanding the glass transition. Although many theories, such as the Adam-Gibbs theory, have sought to relate growing relaxation times to length scales associated with spatial correlations in liquid structure or motion of molecules, the role of length scales in glassy dynamics is not well established. Recent studies of spatially correlated rearrangements of molecules leading to structural relaxation, termed "spatially heterogeneous dynamics," provide fresh impetus in this direction. A powerful approach to extract length scales in critical phenomena is finite-size scaling, wherein a system is studied for sizes traversing the length scales of interest. We perform finite-size scaling for a realistic glass-former, using computer simulations, to evaluate the length scale associated with spatially heterogeneous dynamics, which grows as temperature decreases. However, relaxation times that also grow with decreasing temperature do not exhibit standard finite-size scaling with this length. We show that relaxation times are instead determined, for all studied system sizes and temperatures, by configurational entropy, in accordance with the Adam-Gibbs relation, but in disagreement with theoretical expectations based on spin-glass models that configurational entropy is not relevant at temperatures substantially above the critical temperature of mode-coupling theory. Our results provide new insights into the dynamics of glass-forming liquids and pose serious challenges to existing theoretical descriptions.

Facilitation, complexity growth, mode coupling, and activated dynamics in supercooled liquids

In low-temperature-supercooled liquids, below the ideal mode-coupling theory transition temperature, hopping and continuous diffusion are seen to coexist. Here, we present a theory that shows explicitly the interplay between the two processes and shows that activated hopping facilitates continuous diffusion in the otherwise frozen liquid. Several universal features arise from nonlinear interactions between the continuous diffusive dynamics described here by the mode coupling theory (MCT)] and the activated hopping (described here by the random first-order transition theory). We apply the theory to a specific system, Salol, to show that the theory correctly predicts the temperature dependence of the nonexponential stretching parameter, beta, and the primary alpha relaxation timescale, tau. The study explains why, even below the mean field ergodic to nonergodic transition, the dynamics is well described by MCT. The nonlinear coupling between the two dynamical processes modifies the relaxation behavior of the structural relaxation from what would be predicted by a theory with a complete static Gaussian barrier distribution in a manner that may be described as a facilitation effect. Furthermore, the theory correctly predicts the observed variation of the stretching exponent beta with the fragility parameter, D. These two predictions also allow the complexity growth to be predicted, in good agreement with the results of Capaccioli et al. [Capaccioli S, Ruocco G, Zamponi F (2008) J Phys Chem B 112:10652-10658].

Collisions, caging, thermodynamics, and jamming in the barrier hopping theory of glassy hard sphere fluids

An ultralocal limit of the microscopic single particle barrier hopping theory of glassy dynamics is proposed which allows explicit analytic expressions for the characteristic length scales, energy scales, and nonequilibrium free energy to be derived. All properties are shown to be controlled by a single coupling constant determined by the fluid density and contact value of the radial distribution function. This parameter quantifies an effective mean square force exerted on a tagged particle due to collisions with its surroundings. The analysis suggests a conceptual basis for previous surprising findings of multiple inter-relationships between characteristics of the transient localized state, the early stages of cage escape, non-Gaussian or dynamic heterogeneity effects, and the barrier hopping process that defines the alpha relaxation event. The underlying physical picture is also relevant to fluids of nonspherical molecules and sticky colloidal suspensions. The possibility of a unified view of liquid dynamics is suggested spanning the range from dense gases to the zero mobility jammed state.

Understanding fragility in supercooled Lennard-Jones mixtures. I. Locally preferred structures

The existence of systematic variations of isobaric fragility in different supercooled Lennard-Jones binary mixtures is revealed by molecular dynamics simulations. The connection between fragility and local structures in the bulk is analyzed by means of a Voronoi construction. It is found that clusters of particles belonging to locally preferred structures form slow, long-lived domains, whose spatial extension increases with decreasing temperature. As a general rule, a more rapid growth, upon supercooling, of such domains is associated with a more pronounced super-Arrhenius behavior, and hence to a larger fragility.

Theory of relaxation and elasticity in polymer glasses

The recently developed activated barrier hopping theory of deeply supercooled polymer melts [K. S. Schweizer and E. J. Saltzman, J. Chem. Phys. 121, 1984 (2004)] is extended to the nonequilibrium glass state. Below the kinetic glass temperature T(g), the exact statistical mechanical relation between the dimensionless amplitude of long wavelength density fluctuations, S(0), and the thermodynamic compressibility breaks down. Proper extension of the theory requires knowledge of the nonequilibrium S(0) which x-ray scattering experiments find to consist of a material specific and temperature-independent quenched disorder contribution plus a vibrational contribution which varies roughly linearly with temperature. Motivated by these experiments and general landscape concepts, a simple model is proposed for S(0)(T). Deep in the glass state the form of the temperature dependence of the segmental relaxation time is found to depend sensitively on the magnitude of frozen in density fluctuations. At the (modest) sub-T(g) temperatures typically probed in experiment, an effective Arrhenius behavior is generically predicted which is of nonequilibrium origin. The change in apparent activation energy across the glass transition is determined by the amplitude of frozen density fluctuations. For values of the latter consistent with experiment, the theory predicts a ratio of effective activation energies in the range of 3-6, in agreement with multiple measurements. Calculations of the shear modulus for atactic polymethylmethacrylate above and below the glass transition temperature have also been performed. The present work provides a foundation for the formulation of predictive theories of physical aging, the influence of deformation on the alpha relaxation process, and rate-dependent nonlinear mechanical properties of thermoplastics.

Ideal glass transitions, shear modulus, activated dynamics, and yielding in fluids of nonspherical objects

An extension of naive ideal mode coupling theory (MCT) and its generalization to treat activated barrier hopping and glassy dynamics in fluids and suspensions composed of nonspherical hard core objects is proposed. An effective center-of-mass description is adopted. It corresponds to a specific type of pre-averaging of the dynamical consequences of orientational degrees of freedom. The simplest case of particles composed of symmetry-equivalent interaction sites is considered. The theory is implemented for a homonuclear diatomic shape of variable bond length. The naive MCT glass transition boundary is predicted to be a nonmonotonic function of the length-to-width or aspect ratio and occurs at a nearly unique value of the dimensionless compressibility. The latter quantifies the amplitude of long wavelength thermal density fluctuations, thereby (empirically) suggesting a tight connection between the onset of localization and thermodynamics. Localization lengths and elastic shear moduli for different aspect ratio and volume fraction systems approximately collapse onto master curves based on a reduced volume fraction variable that quantifies the distance from the ideal glass transition. Calculations of the entropic barrier height and hopping time, maximum restoring force, and absolute yield stress and strain as a function of diatomic aspect ratio and volume fraction have been performed. Strong correlations of these properties with the dimensionless compressibility are also found, and nearly universal dependences have been numerically identified based on property-specific nondimensionalizations. Generalization of the approach to rigid rods, disks, and variable shaped molecules is possible, including oriented liquid crystalline phases.