Glass transition of hard spheres in high dimensions

Simple Fluctuations in Simple Glass Formers

Critical single-particle fluctuations associated with particle displacements are inherent to simple glass-forming liquids in the limit of large dimensions and leave a pseudocritical trace across all finite dimensions. This characteristic could serve as a crucial test for distinguishing between theories of glass formation. We here examine these critical fluctuations, as captured by the well-established non-Gaussian parameter, within both mode-coupling theory (MCT) and dynamical mean-field theory (DMFT) across dimensions for hard sphere fluids and for the minimally structured Mari-Kurchan model. We establish general scaling laws relevant to any liquid dynamics theory in large dimensions and show that the dimensional scalings predicted by MCT are inconsistent with those from DMFT. Simulation results for hard sphere fluids in moderately high dimensions align with the DMFT scenario, reinforcing the relevance of mean-field theory for capturing glass physics in finite dimensions. We identify potential adjustments to MCT to account for certain mean-field physics. Our findings also highlight that local structure and spatial dimensionality can affect single-particle critical fluctuations in nontrivial ways.

An alternative, dynamic density functional-like theory for time-dependent density fluctuations in glass-forming fluids

We propose an alternative theory for the relaxation of density fluctuations in glass-forming fluids. We derive an equation of motion for the density correlation function that is local in time and is similar in spirit to the equation of motion for the average non-uniform density profile derived within the dynamic density functional theory. We identify the Franz-Parisi free energy functional as the non-equilibrium free energy for the evolution of the density correlation function. An appearance of a local minimum of this functional leads to a dynamic arrest. Thus, the ergodicity breaking transition predicted by our theory coincides with the dynamic transition of the static approach based on the same non-equilibrium free energy functional.

Efficient d-dimensional molecular dynamics simulations for studies of the glass-jamming transition

We develop an algorithm suitable for parallel molecular dynamics simulations in d spatial dimensions and describe its implementation in C++. All routines work in arbitrary d; the maximum simulated d is limited only by available computing resources. These routines include several that are particularly useful for studies of the glass-jamming transition, such as SWAP Monte Carlo and FIRE energy minimization. The scalings of simulation runtimes with the number of particles N and number of simulation threads n_{threads} are comparable to popular molecular dynamics codes such as LAMMPS. The efficient parallel implementation allows simulation of systems that are much larger than those employed in previous high-dimensional glass-transition studies. As a demonstration of the code's capabilities, we show that supercooled d=6 liquids can possess dynamics that are substantially more heterogeneous and experience a breakdown of the Stokes-Einstein relation that is substantially stronger than previously reported, owing at least in part to the much smaller system sizes employed in earlier simulations.

Dynamics of liquids in the large-dimensional limit

In this paper we analytically derive the exact closed dynamical equations for a liquid with short-ranged interactions in large spatial dimensions using the same statistical mechanics tools employed to analyze Brownian motion. Our derivation greatly simplifies the original path-integral-based route to these equations and provides insight into the physical features associated with high-dimensional liquids and glass formation. Most importantly, our construction provides a route to the exact dynamical analysis of important related dynamical problems, as well as a means to devise cluster generalizations of the exact solution in infinite dimensions. This latter fact opens the door to the construction of increasingly accurate theories of vitrification in three-dimensional liquids.

Mean-Field Caging in a Random Lorentz Gas

The random Lorentz gas (RLG) is a minimal model of both percolation and glassiness, which leads to a paradox in the infinite-dimensional, d → ∞ limit: the localization transition is then expected to be continuous for the former and discontinuous for the latter. As a putative resolution, we have recently suggested that, as d increases, the behavior of the RLG converges to the glassy description and that percolation physics is recovered thanks to finite-d perturbative and nonperturbative (instantonic) corrections [Biroli et al. Phys. Rev. E 2021, 103, L030104]. Here, we expand on the d → ∞ physics by considering a simpler static solution as well as the dynamical solution of the RLG. Comparing the 1/d correction of this solution with numerical results reveals that even perturbative corrections fall out of reach of existing theoretical descriptions. Comparing the dynamical solution with the mode-coupling theory (MCT) results further reveals that, although key quantitative features of MCT are far off the mark, it does properly capture the discontinuous nature of the d → ∞ RLG. These insights help chart a path toward a complete description of finite-dimensional glasses.

Multiple glass transitions and higher-order replica symmetry breaking of binary mixtures

We extend the replica liquid theory in order to describe the multiple glass transitions of binary mixtures with large size disparities, by taking into account the two-step replica symmetry breaking (2RSB). We determine the glass phase diagram of the mixture of large and small particles in the large-dimension limit where the mean-field theory becomes exact. When the size ratio of particles is beyond a critical value, the theory predicts three distinct glass phases; (i) the one-step replica symmetery breaking (1RSB) double glass where both components vitrify simultaneously, (ii) the 1RSB single glass where only large particles are frozen while small particles remain mobile, and (iii) a glass phase called the 2RSB double glass where both components vitrify simultaneously but with an energy landscape topography distinct from the 1RSB double glass.

Connecting real glasses to mean-field models

We propose a novel model for a glass-forming liquid, which allows us to switch in a continuous manner from a standard three-dimensional liquid to a fully connected mean-field model. This is achieved by introducing k additional particle-particle interactions, which thus augments the effective number of neighbors of each particle. Our computer simulations of this system show that the structure of the liquid does not change with the introduction of these pseudo-neighbors and by means of analytical calculations, and we determine the structural properties related to these additional neighbors. We show that the relaxation dynamics of the system slows down very quickly with the increase in k and that the onset and the mode-coupling temperatures increase. The systems with high values of k follow the mode-coupling theory power law behavior for a larger temperature range compared to the ones with lower values of k. The dynamic susceptibility indicates that the dynamic heterogeneity decreases with the increase in k, whereas the non-Gaussian parameter is independent of it. Thus, we conclude that with the increase in the number of pseudo-neighbors, the system becomes more mean-field-like. By comparing our results with previous studies on mean-field-like systems, we come to the conclusion that the details of how the mean-field limit is approached are important since they can lead to different dynamical behavior in this limit.

Generalized mode-coupling theory of the glass transition. I. Numerical results for Percus-Yevick hard spheres

Mode-coupling theory (MCT) constitutes one of the few first-principles-based approaches to describe the physics of the glass transition, but the theory's inherent approximations compromise its accuracy in the activated glassy regime. Here, we show that microscopic generalized mode-coupling theory (GMCT), a recently proposed hierarchical framework to systematically improve upon MCT, provides a promising pathway toward a more accurate first-principles description of glassy dynamics. We present a comprehensive numerical analysis for Percus-Yevick hard spheres by performing explicitly wavenumber- and time-dependent GMCT calculations up to sixth order. Specifically, we calculate the location of the critical point, the associated non-ergodicity parameters, and the time-dependent dynamics of the density correlators at both absolute and reduced packing fractions, and we test several universal scaling relations in the α- and β-relaxation regimes. It is found that higher-order GMCT can successfully remedy some of MCT's pathologies, including an underestimation of the critical glass transition density and an overestimation of the hard-sphere fragility. Furthermore, we numerically demonstrate that the celebrated scaling laws of MCT are preserved in GMCT and that the predicted critical exponents manifestly improve as more levels are incorporated in the GMCT hierarchy. Although formally the GMCT equations should be solved up to infinite order to reach full convergence, our finite-order GMCT calculations unambiguously reveal a uniform convergence pattern for the dynamics. We thus argue that GMCT can provide a feasible and controlled means to bypass MCT's main uncontrolled approximation, offering hope for the future development of a quantitative first-principles theory of the glass transition.

Generalized mode-coupling theory of the glass transition. II. Analytical scaling laws

Generalized mode-coupling theory (GMCT) constitutes a systematically correctable, first-principles theory to study the dynamics of supercooled liquids and the glass transition. It is a hierarchical framework that, through the incorporation of increasingly many particle density correlations, can remedy some of the inherent limitations of the ideal mode-coupling theory (MCT). However, despite MCT's limitations, the ideal theory also enjoys several remarkable successes, notably including the analytical scaling laws for the α- and β-relaxation dynamics. Here, we mathematically derive similar scaling laws for arbitrary-order multi-point density correlation functions obtained from GMCT under arbitrary mean-field closure levels. More specifically, we analytically derive the asymptotic and preasymptotic solutions for the long-time limits of multi-point density correlators, the critical dynamics with two power-law decays, the factorization scaling laws in the β-relaxation regime, and the time-density superposition principle in the α-relaxation regime. The two characteristic power-law-divergent relaxation times for the two-step decay and the non-trivial relation between their exponents are also obtained. The validity ranges of the leading-order scaling laws are also provided by considering the leading preasymptotic corrections. Furthermore, we test these solutions for the Percus-Yevick hard-sphere system. We demonstrate that GMCT preserves all the celebrated scaling laws of MCT while quantitatively improving the exponents, rendering the theory a promising candidate for an ultimately quantitative first-principles theory of glassy dynamics.

Active Hard Spheres in Infinitely Many Dimensions

Few equilibrium-even less so nonequilibrium-statistical-mechanical models with continuous degrees of freedom can be solved exactly. Classical hard spheres in infinitely many space dimensions are a notable exception. We show that, even without resorting to a Boltzmann distribution, dimensionality is a powerful organizing device for exploring the stationary properties of active hard spheres evolving far from equilibrium. In infinite dimensions, we exactly compute the stationary state properties that govern and characterize the collective behavior of active hard spheres: the structure factor and the equation of state for the pressure. In turn, this allows us to account for motility-induced phase separation. Finally, we determine the crowding density at which the effective propulsion of a particle vanishes.

Active glasses

Active glassy matter has recently emerged as a novel class of non-equilibrium soft matter, combining energy-driven, active particle movement with dense and disordered glass-like behavior. Here we review the state-of-the-art in this field from an experimental, numerical, and theoretical perspective. We consider both non-living and living active glassy systems, and discuss how several hallmarks of glassy dynamics (dynamical slowdown, fragility, dynamical heterogeneity, violation of the Stokes-Einstein relation, and aging) are manifested in such materials. We start by reviewing the recent experimental evidence in this area of research, followed by an overview of the main numerical simulation studies and physical theories of active glassy matter. We conclude by outlining several open questions and possible directions for future work.

On mean-field theories of dynamics in supercooled liquids

We develop a hybrid numerical approach to extract the exact memory function K(t) of a tagged particle in three-dimensional glass-forming liquids. We compare the behavior of the exact memory function to two mean-field approaches, namely, the standard mode-coupling theory and a recently proposed ansatz for the memory function that forms the basis of a new derivation of the exact form of K(t) for a fluid with short-ranged interactions in infinite dimensions. Each of the mean-field functions qualitatively and quantitatively share traits with the exact K(t), although several important quantitative differences are manifest.

Simple Theory for the Dynamics of Mean-Field-Like Models of Glass-Forming Fluids

We propose a simple theory for the dynamics of model glass-forming fluids, which should be solvable using a mean-field-like approach. The theory is based on transparent physical assumptions, which can be tested in computer simulations. The theory predicts an ergodicity-breaking transition that is identical to the so-called dynamic transition predicted within the replica approach. Thus, it can provide the missing dynamic component of the random first order transition framework. In the large-dimensional limit the theory reproduces the result of a recent exact calculation of Maimbourg et al. [Phys. Rev. Lett. 116, 015902 (2016)PRLTAO0031-900710.1103/PhysRevLett.116.015902]. Our approach provides an alternative, physically motivated derivation of this result.

Solution of the Dynamics of Liquids in the Large-Dimensional Limit

We obtain analytic expressions for the time correlation functions of a liquid of spherical particles, exact in the limit of high dimensions d. The derivation is long but straightforward: a dynamic virial expansion for which only the first two terms survive, followed by a change to generalized spherical coordinates in the dynamic variables leading to saddle-point evaluation of integrals for large d. The problem is, thus, mapped onto a one-dimensional diffusion in a perturbed harmonic potential with colored noise. At high density, an ergodicity-breaking glass transition is found. In this regime, our results agree with thermodynamics, consistently with the general random first order transition scenario. The glass transition density is higher than the best known lower bound for hard sphere packings in large d. Because our calculation is, if not rigorous, elementary, an improvement in the bound for sphere packings in large dimensions is at hand.

Quantitative approximation schemes for glasses

By means of a systematic expansion around the infinite-dimensional solution, we obtain an approximation scheme to compute properties of glasses in low dimensions. The resulting equations take as input the thermodynamic and structural properties of the equilibrium liquid, and from this they allow one to compute properties of the glass. They are therefore similar in spirit to the Mode Coupling approximation scheme. Our scheme becomes exact, by construction, in dimension d→∞, and it can be improved systematically by adding more terms in the expansion.

Crossover from β to α Relaxation in Cooperative Facilitation Dynamics

β and α relaxation processes are dynamical scaling regimes of glassy systems occurring on two separate time scales which both diverge as the glass state is approached. We study here the crossover scaling from β to α relaxation in the cooperative facilitation scenario (CFS) and show that it is quantitatively described, with no adjustable parameter, by the leading order asymptotic formulas for scaling predicted by the mode-coupling theory (MCT). These results establish (i) the mutual universality of the MCT and CFS, and (ii) the existence of a purely dynamic realization of MCT, which is distinct from the well-established random first order transition scenario for disordered systems. Some implications of the emerging kinetic-static duality are discussed.

Validity of the Stokes-Einstein Relation in Soft Colloids up to the Glass Transition

We investigate the dynamics of kinetically frozen block copolymer micelles of different softness across a wide range of particle concentrations, from the fluid to the onset of glassy behavior, through a combination of rheology, dynamic light scattering, and pulsed field gradient NMR spectroscopy. We additionally perform Brownian dynamics simulations based on an ultrasoft coarse-grained potential, which are found to be in quantitative agreement with experiments, capturing even the very details of dynamic structure factors S(Q,t) on approaching the glass transition. We provide evidence that for these systems the Stokes-Einstein relation holds up to the glass transition; given that it is violated for dense suspensions of hard colloids, our findings suggest that its validity is an intriguing signature of ultrasoft interactions.

Renormalized dynamics of the Dean-Kawasaki model

We study the model of a supercooled liquid for which the equation of motion for the coarse-grained density ρ(x,t) is the nonlinear diffusion equation originally proposed by Dean and Kawasaki, respectively, for Brownian and Newtonian dynamics of fluid particles. Using a Martin-Siggia-Rose (MSR) field theory we study the renormalization of the dynamics in a self-consistent form in terms of the so-called self-energy matrix Σ. The appropriate model for the renormalized dynamics involves an extended set of field variables {ρ,θ}, linked through a nonlinear constraint. The latter incorporates, in a nonperturbative manner, the effects of an infinite number of density nonlinearities in the dynamics. We show that the contributing element of Σ which renormalizes the bare diffusion constant D(0) to D(R) is same as that proposed by Kawasaki and Miyazima [Z. Phys. B Condens. Matter 103, 423 (1997)]. D(R) sharply decreases with increasing density. We consider the likelihood of a ergodic-nonergodic (ENE) transition in the model beyond a critical point. The transition is characterized by the long-time limit of the density correlation freezing at a nonzero value. From our analysis we identify an element of Σ which arises from the above-mentioned nonlinear constraint and is key to the viability of the ENE transition. If this self-energy would be zero, then the model supports a sharp ENE transition with D(R)=0 as predicted by Kawasaki and Miyazima. With the full model having nonzero value for this self-energy, the density autocorrelation function decays to zero in the long-time limit. Hence the ENE transition is not supported in the model.

Glass transitions and scaling laws within an alternative mode-coupling theory

Idealized glass transitions are discussed within an alternative mode-coupling theory (TMCT) proposed by Tokuyama [Physica A 395, 31 (2014)]. This is done in order to identify common ground with and differences from the conventional mode-coupling theory (MCT). It is proven that both theories imply the same scaling laws for the transition dynamics, which are characterized by two power-law decay functions and two diverging power-law time scales. However, the values for the corresponding anomalous exponents calculated within both theories differ from each other. It is proven that the TMCT, contrary to the MCT, does not describe transitions with continuously vanishing arrested parts of the correlation functions. It is also demonstrated for a schematic model that the TMCT does not lead to the MCT scenarios either for transition-line crossings or for the appearance of higher-order glass-transition singularities.

Dimensional study of the dynamical arrest in a random Lorentz gas

The random Lorentz gas (RLG) is a minimal model for transport in heterogeneous media. Upon increasing the obstacle density, it exhibits a growing subdiffusive transport regime and then a dynamical arrest. Here, we study the dimensional dependence of the dynamical arrest, which can be mapped onto the void percolation transition for Poisson-distributed point obstacles. We numerically determine the arrest in dimensions d=2-6. Comparison of the results with standard mode-coupling theory reveals that the dynamical theory prediction grows increasingly worse with d. In an effort to clarify the origin of this discrepancy, we relate the dynamical arrest in the RLG to the dynamic glass transition of the infinite-range Mari-Kurchan-model glass former. Through a mixed static and dynamical analysis, we then extract an improved dimensional scaling form as well as a geometrical upper bound for the arrest. The results suggest that understanding the asymptotic behavior of the random Lorentz gas may be key to surmounting fundamental difficulties with the mode-coupling theory of glasses.

Equilibrium dynamics of the Dean-Kawasaki equation: mode-coupling theory and its extension

We extend a previously proposed field-theoretic self-consistent perturbation approach for the equilibrium dynamics of the Dean-Kawasaki equation presented in [Kim and Kawasaki, J. Stat. Mech. (2008) P02004]. By taking terms missing in the latter analysis into account we arrive at a set of three new equations for correlation functions of the system. These correlations involve the density and its logarithm as local observables. Our new one-loop equations, which must carefully deal with the noninteracting Brownian gas theory, are more general than the historic mode-coupling one in that a further approximation corresponding to Gaussian density fluctuations leads back to the original mode-coupling equation for the density correlations alone. However, without performing any further approximation step, our set of three equations does not feature any ergodic-nonergodic transition, as opposed to the historical mode-coupling approach.

Jammed lattice sphere packings

We generate and study an ensemble of isostatic jammed hard-sphere lattices. These lattices are obtained by compression of a periodic system with an adaptive unit cell containing a single sphere until the point of mechanical stability. We present detailed numerical data about the densities, pair correlations, force distributions, and structure factors of such lattices. We show that this model retains many of the crucial structural features of the classical hard-sphere model and propose it as a model for the jamming and glass transitions that enables exploration of much higher dimensions than are usually accessible.

Perspective: The glass transition

We provide here a brief perspective on the glass transition field. It is an assessment, written from the point of view of theory, of where the field is and where it seems to be heading. We first give an overview of the main phenomenological characteristics, or "stylised facts," of the glass transition problem, i.e., the central observations that a theory of the physics of glass formation should aim to explain in a unified manner. We describe recent developments, with a particular focus on real space properties, including dynamical heterogeneity and facilitation, the search for underlying spatial or structural correlations, and the relation between the thermal glass transition and athermal jamming. We then discuss briefly how competing theories of the glass transition have adapted and evolved to account for such real space issues. We consider in detail two conceptual and methodological approaches put forward recently, that aim to access the fundamental critical phenomenon underlying the glass transition, be it thermodynamic or dynamic in origin, by means of biasing of ensembles, of configurations in the thermodynamic case, or of trajectories in the dynamic case. We end with a short outlook.

Cooperative heterogeneous facilitation: multiple glassy states and glass-glass transition

The formal structure of glass singularities in the mode-coupling theory (MCT) of supercooled liquids dynamics is closely related to that appearing in the analysis of heterogeneous bootstrap percolation on Bethe lattices, random graphs, and complex networks. Starting from this observation one can build up microscopic on-lattice realizations of schematic MCT based on cooperative facilitated spin mixtures. I discuss a microscopic implementation of the F(13) schematic model including multiple glassy states and the glass-glass transition. Results suggest that our approach is flexible enough to bridge alternative theoretical descriptions of glassy matter based on the notions of quenched disorder and dynamic facilitation.

Dimensional study of the caging order parameter at the glass transition

The glass problem is notoriously hard and controversial. Even at the mean-field level, little is agreed upon regarding why a fluid becomes sluggish while exhibiting but unremarkable structural changes. It is clear, however, that the process involves self-caging, which provides an order parameter for the transition. It is also broadly assumed that this cage should have a gaussian shape in the mean-field limit. Here we show that this ansatz does not hold. By performing simulations as a function of spatial dimension d, we find the cage to keep a nontrivial form. Quantitative mean-field descriptions of the glass transition, such as mode-coupling theory, density functional theory, and replica theory, all miss this crucial element. Although the mean-field random first-order transition scenario of the glass transition is qualitatively supported here and non-mean-field corrections are found to remain small on decreasing d, reconsideration of its implementation is needed for it to result in a coherent description of experimental observations.

Mode-coupling theory of the glass transition for confined fluids

We present a detailed derivation of a microscopic theory for the glass transition of a liquid enclosed between two parallel walls relying on a mode-coupling approximation. This geometry lacks translational invariance perpendicular to the walls, which implies that the density profile and the density-density correlation function depends explicitly on the distances to the walls. We discuss the residual symmetry properties in slab geometry and introduce a symmetry adapted complete set of two-point correlation functions. Since the currents naturally split into components parallel and perpendicular to the walls the mathematical structure of the theory differs from the established mode-coupling equations in bulk. We prove that the equations for the nonergodicity parameters still display a covariance property similar to bulk liquids.

Time dependent stretching of aging dynamics in a generalized hydrodynamic model for supercooled liquids

The nonequilibrium dynamics and aging behavior of a supercooled liquid is investigated from an analysis of the correlation of density fluctuations at two different times. The dynamic correlation functions are computed by solving numerically the equations of nonlinear fluctuating hydrodynamics. The aging time dependence follows a modified stretched exponential form with a relaxation time which is dependent on the aging time. This is similar to the behavior seen in the aging data of dielectric response functions of a typical glass forming liquid.

Mode-coupling theory predictions for the dynamical transitions of partly pinned fluid systems

The predictions of the mode-coupling theory (MCT) for the dynamical arrest scenarios in a partly pinned (PP) fluid system are reported. The corresponding dynamical phase diagram is found to be very similar to that of a related quenched-annealed (QA) system. The only significant qualitative difference lies in the shape of the diffusion-localization lines at high matrix densities, with a reentry phenomenon for the PP system but not for the QA model, in full agreement with recent computer simulation results. This finding clearly lends support to the predictive power of the MCT for fluid-matrix systems. In addition, the predictions of the MCT are shown to be in stark contrast with those of the random first-order transition theory. The PP systems are thus confirmed as very promising models for differentiating tests of theories of the glass transition.

Glass transition and random close packing above three dimensions

Motivated by a recently identified severe discrepancy between a static and a dynamic theory of glasses, we numerically investigate the behavior of dense hard spheres in spatial dimensions 3 to 12. Our results are consistent with the static replica theory, but disagree with the dynamic mode-coupling theory, indicating that key ingredients of high-dimensional physics are missing from the latter. We also obtain numerical estimates of the random close packing density, which provides new insights into the mathematical problem of packing spheres in large dimensions.

Glassy dynamics and hysteresis in a linear system of orientable hard rods

We study the dynamics of a one-dimensional fluid of orientable hard rectangles with a non-coarse-grained microscopic mechanism of facilitation. The length occupied by a rectangle depends on its orientation, which is a discrete variable coupled to an external field. The equilibrium properties of our model are essentially those of the Tonks gas, but at high densities the orientational degrees of freedom become effectively frozen due to jamming. This is a simple analytically tractable model of the glassy phase. Under a cyclic variation of the pressure, hysteresis is observed. Following a pressure quench, the orientational persistence exhibits a two-stage decay characteristic of glassy systems.

Tests of mode-coupling theory in two dimensions

We analyze the glassy dynamics of binary mixtures of hard disks in two dimensions. Predictions of the mode-coupling theory (MCT) are tested with extensive Brownian dynamics simulations. Measuring the collective particle density correlation functions in the vicinity of the glass transition, we verify four predicted mixing effects. For instance, for large size disparities, adding a small amount of small particles at a fixed packing fraction leads to a speedup in the long-time dynamics, while for small size disparities it leads to a slowing-down. Qualitative features of the nonergodicity parameters and the β relaxation, which both depend in a nontrivial way on the mixing ratio, are found in the simulated correlators. Studying one system in detail, we are able to determine its ideal MCT glass transition point as φ(c)=0.7948 and test MCT predictions quantitatively.

Glass transition of the monodisperse Gaussian core model

We numerically investigate the dynamical properties of the one-component Gaussian core model in supercooled states. We find that nucleation is increasingly suppressed with increasing density. The system concomitantly exhibits glassy, slow dynamics characterized by the two-step stretched exponential relaxation of the density correlation and a drastic increase of the relaxation time. We also find a weaker violation of the Stokes-Einstein relation and a smaller non-Gaussian parameter than in typical model glass formers, implying weaker dynamic heterogeneities. Additionally, the agreement of the simulation data with the prediction of mode-coupling theory is exceptionally good, indicating that the nature of the slow dynamics of this ultrasoft particle fluid is mean-field-like. This fact may be understood as a consequence of the long-range nature of the interaction.

Application of Edwards' statistical mechanics to high-dimensional jammed sphere packings

The isostatic jamming limit of frictionless spherical particles from Edwards' statistical mechanics [Song et al., Nature (London) 453, 629 (2008)] is generalized to arbitrary dimension d using a liquid-state description. The asymptotic high-dimensional behavior of the self-consistent relation is obtained by saddle-point evaluation and checked numerically. The resulting random close packing density scaling ϕ∼d2(-d) is consistent with that of other approaches, such as replica theory and density-functional theory. The validity of various structural approximations is assessed by comparing with three- to six-dimensional isostatic packings obtained from simulations. These numerical results support a growing accuracy of the theoretical approach with dimension. The approach could thus serve as a starting point to obtain a geometrical understanding of the higher-order correlations present in jammed packings.

Mode-coupling theory as a mean-field description of the glass transition

Mode-coupling theory (MCT) is conjectured to be a mean-field description of dynamics of the structural glass transition and the replica theory to be its thermodynamic counterpart. However, the relationship between the two theories remains controversial and a quantitative comparison is lacking. In this Letter, we investigate MCT for monatomic hard-sphere fluids at arbitrary dimensions above three and compare the results with replica theory. We find grave discrepancies between the predictions of two theories. While MCT describes the nonergodic parameter quantitatively better than the replica theory in three dimensions, it predicts a completely different dimension dependence of the dynamical transition point. We find it to be due to the pathological behavior of the nonergodic parameters derived from MCT, which exhibit negative tails in real space at high dimensions.