Source author record

Ludovic Berthier

Ludovic Berthier appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

85works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

85 published item(s)

preprint2023arXiv

Microscopic observation of two-level systems in a metallic glass model

The low-temperature quasi-universal behavior of amorphous solids has been attributed to the existence of spatially-localized tunneling defects found in the low-energy regions of the potential energy landscape. Computational models of glasses can be studied to elucidate the microscopic nature of these defects. Recent simulation work has demonstrated the means of generating stable glassy configurations for models that mimic metallic glasses using the swap Monte Carlo algorithm. Building on these studies, we present an extensive exploration of the glassy metabasins of the potential energy landscape of a variant of the most widely used model of metallic glasses. We carefully identify tunneling defects and reveal their depletion with increased glass stability. The density of tunneling defects near the experimental glass transition temperature appears to be in good agreement with experimental measurements.

preprint2022arXiv

Collective dynamics in a glass-former with Mari-Kurchan interactions

We numerically study the equilibrium relaxation dynamics of a two-dimensional Mari-Kurchan glass model. The tree-like structure of particle interactions forbids both non-trivial structural motifs and the emergence of a complex free-energy landscape leading to a thermodynamic glass transition, while the finite-dimensional nature of the model prevents the existence of a mode-coupling singularity. Nevertheless, the equilibrium relaxation dynamics is shown to be in excellent agreement with simulations performed in conventional glass-formers. Averaged time-correlation functions display a phenomenology typical of supercooled liquids, including the emergence of an excess signal in relaxation spectra at intermediate frequencies. We show that this evolution is accompanied by strong signatures of collective and heterogeneous dynamics which cannot be interpreted in terms of single particle hopping and emerge from dynamic facilitation. Our study demonstrates that an off-lattice interacting particle model with extremely simple structural correlations displays quantitatively realistic glassy dynamics.

preprint2022arXiv

Disordered collective motion in dense assemblies of persistent particles

We explore the emergence of nonequilibrium collective motion in disordered non-thermal active matter when persistent motion and crowding effects compete, using simulations of a two-dimensional model of size polydisperse self-propelled particles. In stark contrast with monodisperse systems, we find that polydispersity stabilizes a homogeneous active liquid at arbitrary large persistence times, characterized by remarkable velocity correlations and irregular turbulent flows. For all persistence values, the active fluid undergoes a nonequilibrium glass transition at large density. This is accompanied by collective motion, whose nature evolves from near-equilibrium spatially heterogeneous dynamics at small persistence, to a qualitatively different intermittent dynamics when persistence is large. This latter regime involves a complex time evolution of the correlated displacement field

preprint2022arXiv

Microscopic origin of excess wings in relaxation spectra of supercooled liquids

Glass formation is encountered in diverse materials. Experiments have revealed that dynamic relaxation spectra of supercooled liquids generically become asymmetric near the glass transition temperature, $T_g$, where an extended power law emerges at high frequencies. The microscopic origin of this "wing" remains unknown, and was so far inaccessible to simulations. Here, we develop a novel computational approach and study the equilibrium dynamics of model supercooled liquids near $T_g$. We demonstrate the emergence of a power law wing in numerical spectra, which originates from relaxation at rare, localised regions over broadly-distributed timescales. We rationalise the asymmetric shape of relaxation spectra by constructing an empirical model associating heterogeneous activated dynamics with dynamic facilitation, which are the two minimal physical ingredients revealed by our simulations. Our work offers a glimpse of the molecular motion responsible for glass formation at relevant experimental conditions.

preprint2022arXiv

Rare events and disorder control the brittle yielding of well-annealed amorphous solids

We use atomistic computer simulations to provide a microscopic description of the brittle failure of amorphous materials, and we assess the role of rare events and quenched disorder. We argue that brittle yielding originates at rare soft regions, similarly to Griffiths effects in disordered systems. We numerically demonstrate how localized plastic events in such soft regions trigger macroscopic failure via the propagation of a shear band. This physical picture, which no longer holds in poorly annealed ductile materials, allows us to discuss the role of finite size effects in brittle yielding and reinforces the similarities between yielding and other disorder-controlled nonequilibrium phase transitions.

preprint2022arXiv

Relaxation dynamics in the energy landscape of glass-forming liquids

We numerically study the zero-temperature relaxation dynamics of several glass-forming models to their inherent structures, following quenches from equilibrium configurations sampled across a wide range of initial temperatures. In a mean-field Mari-Kurchan model, we find that relaxation changes from a power-law to an exponential decay below a well-defined temperature, consistent with recent findings in mean-field $p$-spin models. By contrast, for finite-dimensional systems, the relaxation is always algebraic, with a non-trivial universal exponent at high temperatures crossing over to a harmonic value at low temperatures. We demonstrate that this apparent evolution is controlled by a temperature-dependent population of localised glassy excitations. Our work unifies several recent lines of studies aiming at a detailed characterisation of the complex potential energy landscape of glass-formers, and challenges both mean-field and real space descriptions of glasses.

preprint2022arXiv

Static self-induced heterogeneity in glass-forming liquids: Overlap as a microscope

We propose and numerically implement a local probe of the static self-induced heterogeneity characterizing glass-forming liquids. The method relies on the equilibrium statistics of the overlap between pairs of configurations measured in mesoscopic cavities with unconstrained boundaries. By systematically changing the location of the probed cavity, we directly detect spatial variations of the overlap fluctuations. We provide a detailed analysis of the statistics of a local estimate of the configurational entropy and we infer an estimate of the surface tension between amorphous states, ingredients that are both at the basis of the random first-order transition theory of glass formation. Our results represent the first direct attempt to visualize and quantify the self-induced heterogeneity underpinning the thermodynamics of glass formation. They pave the way for the development of coarse-grained effective theories and for a direct assessment of the role of thermodynamics in the activated dynamics of deeply supercooled liquids.

preprint2022arXiv

Statistical mechanics of coupled supercooled liquids in finite dimensions

We study the statistical mechanics of supercooled liquids when the system evolves at a temperature $T$ with a field $ε$ linearly coupled to its overlap with a reference configuration of the same liquid sampled at a temperature $T_0$. We use mean-field theory to fully characterize the influence of the reference temperature $T_0$, and we mainly study the case of a fixed, low-$T_0$ value in computer simulations. We numerically investigate the extended phase diagram in the $(ε,T)$ plane of model glass-forming liquids in spatial dimensions $d=2$ and $d=3$, relying on umbrella sampling and reweighting techniques. For both $2d$ and $3d$ cases, a similar phenomenology with nontrivial thermodynamic fluctuations of the overlap is observed at low temperatures, but a detailed finite-size analysis reveals qualitatively distinct behaviors. We establish the existence of a first-order transition line for nonzero $ε$ ending in a critical point in the universality class of the random-field Ising model (RFIM) in $d=3$. In $d=2$ instead, no phase transition is found in large enough systems at least down to temperatures below the extrapolated calorimetric glass transition temperature $T_g$. Our results confirm that glass-forming liquid samples of limited size display the thermodynamic fluctuations expected for finite systems undergoing a random first-order transition. They also support the relevance of the physics of the RFIM for supercooled liquids, which may then explain the qualitative difference between $2d$ and $3d$ glass-formers.

preprint2022arXiv

Violation of the fluctuation-dissipation theorem and effective temperatures in spin ice

We present numerical tests of the fluctuation-dissipation theorem (FDT) in the dumbbell model of spin ice with parameters suitable for dysprosium titanate. The tests are made for local spin variables, magnetic monopole density, and energy. We are able to achieve local equilibrium in which the FDT is satisfied down to $T=0.4$ K below which the system completely freezes. Non-equilibrium dynamics, together with violation of the FDT, are nonetheless observed following a thermal quench into the non-contractable monopole pair regime. Despite FDT violation, an approximate linear response regime allows for the identification of effective non-equilibrium temperatures which are different for each variable. The spin variable appears hotter than the heat reservoir, the monopole concentration responds with a lower effective temperature while the energy has a negative effective temperature. Results are discussed in the context of the monopole picture of spin ice and compared to the structure of FDT violations in other glassy materials. Prospectives for future experiments are reviewed.

preprint2021arXiv

Relaxation dynamics of non-Brownian spheres below jamming

We numerically study the relaxation dynamics and associated criticality of non-Brownian frictionless spheres below jamming in spatial dimensions $d=2$, $3$, $4$, and $8$, and in the mean-field Mari-Kurchan model. We discover non-trivial finite-size and volume fraction dependences of the relaxation time associated to the relaxation of unjammed packings. In particular, the relaxation time is shown to diverge logarithmically with system size at any density below jamming, and no critical exponent can characterise its behaviour approaching jamming. In mean-field, the relaxation time is instead well-defined: it diverges at jamming with a critical exponent that we determine numerically and differs from an earlier mean-field prediction. We rationalise the finite $d$ logarithmic divergence using an extreme-value statistics argument in which the relaxation time is dominated by the most connected region of the system. The same argument shows that the earlier proposition that relaxation dynamics and shear viscosity are directly related breaks down in large systems. The shear viscosity of non-Brownian packings is well-defined in all $d$ in the thermodynamic limit, but large finite-size effects plague its measurement close to jamming.

preprint2020arXiv

Analogies between growing dense active matter and soft driven glasses

We develop a minimal model to describe growing dense active matter such as biological tissues, bacterial colonies and biofilms, that are driven by a competition between particle division and steric repulsion. We provide a detailed numerical analysis of collective and single particle dynamics. We show that the microscopic dynamics can be understood as the superposition of an affine radial component due to the global growth, and of a more complex non-affine component which displays features typical of driven soft glassy materials, such as aging, compressed exponential decay of time correlation functions, and a crossover from superdiffusive behaviour at short scales to subdiffusive behaviour at larger scales. This analogy emerges because particle division at the microscale leads to a global expansion which then plays a role analogous to shear flow in soft driven glasses. We conclude that growing dense active matter and sheared dense suspensions can generically be described by the same underlying physics.

preprint2020arXiv

Depletion of two-level systems in ultrastable computer-generated glasses

Amorphous solids exhibit quasi-universal low-temperature anomalies whose origin has been ascribed to localized tunneling defects. Using an advanced Monte Carlo procedure, we create {\it in silico} glasses spanning from hyperquenched to ultrastable glasses. Using a multidimensional path-finding protocol, we locate tunneling defects with energy splittings smaller than $k_{B}T_Q$, with $T_Q$ the temperature below which quantum effects are relevant ($T_Q \approx 1 \,$K in most experiments). We find that as the stability of a glass increases, its energy landscape as well as the manner in which it is probed tend to deplete the density of tunneling defects, as observed in recent experiments. We explore the real-space nature of tunneling defects, and find that they are mostly localized to a few atoms, but are occasionally dramatically delocalized.

preprint2020arXiv

Glass stability changes the nature of yielding under oscillatory shear

We perform molecular dynamics simulations to investigate the effect of a glass preparation on its yielding transition under oscillatory shear. We use swap Monte Carlo to investigate a broad range of glass stabilities from poorly annealed to highly stable systems. We observe a qualitative change in the nature of yielding, which evolves from ductile to brittle as glass stability increases. Our results disentangle the relative role of mechanical and thermal annealing on the mechanical properties of amorphous solids, which is relevant for various experimental situations from the rheology of soft materials to fatigue failure in metallic glasses.

preprint2020arXiv

How to "measure" a structural relaxation time that is too long to be measured?

It has recently become possible to prepare ultrastable glassy materials characterised by structural relaxation times which vastly exceed the duration of any feasible experiment. Similarly, new algorithms have led to the production of ultrastable computer glasses. Is it possible to obtain a reliable estimate of a structural relaxation time that is too long to be measured? We review, organise, and critically discuss various methods to estimate very long relaxation times. We also perform computer simulations of three dimensional ultrastable hard spheres glasses to test and quantitatively compare some of these methods for a single model system. The various estimation methods disagree significantly and it is not yet clear how to accurately estimate extremely long relaxation times.

preprint2020arXiv

On the overlap between configurations in glassy liquids

The overlap, or similarity, between liquid configurations is at the core of the mean-field description of the glass transition, and remains a useful concept when studying three-dimensional glass-forming liquids. In liquids, however, the overlap involves a tolerance, typically of a fraction $a/σ$ of the inter-particle distance, associated with how precisely similar two configurations must be for belonging to the same physically relevant "state". Here, we systematically investigate the dependence of the overlap fluctuations and of the resulting phase diagram when the tolerance is varied over a large range. We show that while the location of the dynamical and thermodynamic glass transitions (if present) is independent of $a/σ$, that of the critical point associated with a transition between a low- and a high-overlap phases in the presence of an applied source nontrivially depends on the value of $a/σ$. We rationalize our findings by using liquid-state theory and the hypernetted chain (HNC) approximation for correlation functions. In addition, we confirm the theoretical trends by studying a three-dimensional glass-former by computer simulations. We show in particular that \rev{a range of $a/σ$ below what is commonly considered maximizes the temperature of the critical point, pushing it up in a liquid region where viscosity is low and computer investigations are easier due to a significantly faster equilibration.

preprint2020arXiv

Role of fluctuations in the yielding transition of two-dimensional glasses

We numerically study yielding in two-dimensional glasses which are generated with a very wide range of stabilities by swap Monte-Carlo simulations and then slowly deformed at zero temperature. We provide strong numerical evidence that stable glasses yield via a nonequilibrium discontinuous transition in the thermodynamic limit. A critical point separates this brittle yielding from the ductile one observed in less stable glasses. We find that two-dimensional glasses yield similarly to their three-dimensional counterparts but display larger sample-to-sample disorder-induced fluctuations, stronger finite-size effects, and rougher spatial wandering of the observed shear bands. These findings strongly constrain effective theories of yielding.

preprint2019arXiv

Brittle yielding of amorphous solids at finite shear rates

Amorphous solids display a ductile to brittle transition as the kinetic stability of the quiescent glass is increased, which leads to a material failure controlled by the sudden emergence of a macroscopic shear band in quasi-static protocols. We numerically study how finite deformation rates influence ductile and brittle yielding behaviors using model glasses in two and three spatial dimensions. We find that a finite shear rate systematically enhances the stress overshoot of poorly-annealed systems, without necessarily producing shear bands. For well-annealed systems, the non-equilibrium discontinuous yielding transition is smeared out by finite shear rates and it is accompanied by the emergence of multiple shear bands that have been also reported in metallic glass experiments. We show that the typical size of the bands and the distance between them increases algebraically with the inverse shear rate. We provide a dynamic scaling argument for the corresponding lengthscale, based on the competition between the deformation rate and the propagation time of the shear bands.

preprint2019arXiv

Finite-dimensional vestige of spinodal criticality above the dynamical glass transition

Finite-dimensional signatures of spinodal criticality are notoriously difficult to come by. The dynamical transition of glass-forming liquids, first described by mode-coupling theory, is a spinodal instability preempted by thermally activated processes that also limit how close the instability can be approached. We combine numerical tools to directly observe vestiges of the spinodal criticality in finite-dimensional glass formers. We use the swap Monte Carlo algorithm to efficiently thermalise configurations beyond the mode-coupling crossover, and analyze their dynamics using a scheme to screen out activated processes, in spatial dimensions ranging from $d=3$ to $d=9$. We observe a strong softening of the mean-field square-root singularity in $d=3$ that is progressively restored as $d$ increases above $d=8$, in surprisingly good agreement with perturbation theory.

preprint2019arXiv

Low-frequency vibrations of jammed packings in large spatial dimensions

Amorphous packings prepared in the vicinity of the jamming transition play a central role in theoretical studies of the vibrational spectrum of glasses. Two mean-field theories predict that the vibrational density of states $g(ω)$ obeys a characteristic power law, $g(ω)\simω^2$, called the non-Debye scaling in the low-frequency region. Numerical studies have however reported that this scaling breaks down at low frequencies, due to finite dimensional effects. In this study, we prepare amorphous packings of up to $128000$ particles in spatial dimensions from $d=3$ to $d=9$ to characterise the range of validity of the non-Debye scaling. Our numerical results suggest that the non-Debye scaling is obeyed down to a frequency that gradually decreases as $d$ increases, and possibly vanishes for large $d$, in agreement with mean-field predictions. We also show that the prestress is an efficient control parameter to quantitatively compare packings across different spatial dimensions.

preprint2019arXiv

Universal relaxation dynamics of sphere packings below jamming

We show that non-Brownian suspensions of repulsive spheres below jamming display a slow relaxational dynamics with a characteristic time scale that diverges at jamming. This slow time scale is fully encoded in the structure of the unjammed packing and can be readily measured via the vibrational density of states. We show that the corresponding dynamic critical exponent is the same for randomly generated and sheared packings. Our results show that a wide variety of physical situations, from suspension rheology to algorithmic studies of the jamming transition are controlled by a unique diverging timescale, with a universal critical exponent.

preprint2016arXiv

Criticality and correlated dynamics at the irreversibility transition in periodically driven colloidal suspensions

One possible framework to interpret the irreversibility transition observed in periodically driven colloidal suspensions is that of a non-equilibrium phase transition towards an absorbing reversible state at low amplitude of the driving force. We consider a simple numerical model for driven suspensions which allows us to characterize in great detail a large body of physical observables that can be experimentally determined to assess the existence and universality class of such a non-equilibrium phase transition. Characterizing the behaviour of static and dynamic correlation functions both in real and Fourier space we determine in particular several critical exponents for our model, which take values that are in good agreement with the universality class of direct ed percolation. We also provide a detailed analysis of single-particle and collective dynamics of the system near the phase transition, which appear intermittent and spatially correlated over diverging timescales and lengthscales, and provide clear signatures of the underlying criticality.

preprint2016arXiv

Efficient measurement of point-to-set correlations and overlap fluctuations in glass-forming liquids

Cavity point-to-set correlations are real-space tools to detect the roughening of the free-energy landscape that accompanies the dynamical slowdown of glass-forming liquids. Measuring these correlations in model glass formers remains, however, a major computational challenge. Here, we develop a general parallel-tempering method that provides orders-of-magnitude improvement for sampling and equilibrating configurations within cavities. We apply this improved scheme to the canonical Kob-Andersen binary Lennard-Jones model for temperatures down to the mode-coupling theory crossover. Most significant improvements are noted for small cavities, which have thus far been the most difficult to study. This methodological advance also enables us to study a broader range of physical observables associated with thermodynamic fluctuations. We measure the probability distribution of overlap fluctuations in cavities, which displays a non-trivial temperature evolution. The corresponding overlap susceptibility is found to provide a robust quantitative estimate of the point-to-set length scale requiring no fitting. By resolving spatial fluctuations of the overlap in the cavity, we also obtain quantitative information about the geometry of overlap fluctuations. We can thus examine in detail how the penetration length as well as its fluctuations evolve with temperature and cavity size.

preprint2016arXiv

Growing timescales and lengthscales characterizing vibrations of amorphous solids

Low-temperature properties of crystalline solids can be understood using harmonic perturbations around a perfect lattice, as in Debye's theory. Low-temperature properties of amorphous solids, however, strongly depart from such descriptions, displaying enhanced transport, activated slow dynamics across energy barriers, excess vibrational modes with respect to Debye's theory (i.e., a Boson Peak), and complex irreversible responses to small mechanical deformations. These experimental observations indirectly suggest that the dynamics of amorphous solids becomes anomalous at low temperatures. Here, we present direct numerical evidence that vibrations change nature at a well-defined location deep inside the glass phase of a simple glass former. We provide a real-space description of this transition and of the rapidly growing time and length scales that accompany it. Our results provide the seed for a universal understanding of low-temperature glass anomalies within the theoretical framework of the recently discovered Gardner phase transition.

preprint2016arXiv

Macroscopic yielding in jammed solids is accompanied by a non-equilibrium first-order transition in particle trajectories

We use computer simulations to analyse the yielding transition during large-amplitude oscillatory shear of a simple model for soft jammed solids. Simultaneous analysis of global mechanical response and particle-scale motion demonstrates that macroscopic yielding, revealed by a smooth crossover in mechanical properties, is accompanied by a sudden change in the particle dynamics, which evolves from non-diffusive motion to irreversible diffusion as the amplitude of the shear is increased. We provide numerical evidence that this sharp change corresponds to a non-equilibrium first-order dynamic phase transition, thus establishing the existence of a well-defined microscopic dynamic signature of the yielding transition in amorphous materials in oscillatory shear.

preprint2016arXiv

Point-to-set lengths, local structure, and glassiness

The growing sluggishness of glass-forming liquids is thought to be accompanied by growing structural order. The nature of such order, however, remains hotly debated. A decade ago, point-to-set (PTS) correlation lengths were proposed as measures of amorphous order in glass formers, but recent results raise doubts as to their generality. Here, we extend the definition of PTS correlations to agnostically capture any type of growing order in liquids, be it local or amorphous. This advance enables the formulation of a clear distinction between slowing down due to conventional critical ordering and that due to glassiness, and provides a unified framework to assess the relative importance of specific local order and generic amorphous order in glass formation.

preprint2016arXiv

The melting of stable glasses is governed by nucleation-and-growth dynamics

We discuss the microscopic mechanisms by which low-temperature amorphous states, such as ultrastable glasses, transform into equilibrium fluids, after a sudden temperature increase. Experiments suggest that this process is similar to the melting of crystals, thus differing from the behaviour found in ordinary glasses. We rationalize these observations using the physical idea that the transformation process takes place very close to a `hidden' equilibrium first-order phase transition, which is observed in systems of coupled replicas. We illustrate our views using simulation results for a simple two-dimensional plaquette spin model, which is known to exhibit a range of glassy behaviour. Our results suggest that nucleation-and-growth dynamics, as found near ordinary first-order transitions, is also the correct theoretical framework to analyse the melting of ultrastable glasses. Our approach provides a unified understanding of multiple experimental observations, such as propagating melting fronts, large kinetic stability ratios, and `giant' dynamic lengthscales.

preprint2016arXiv

The nonequilibrium glassy dynamics of self-propelled particles

We study the glassy dynamics taking place in dense assemblies of athermal active particles that are driven solely by a nonequilibrium self-propulsion mechanism. Active forces are modeled as an Ornstein-Uhlenbeck stochastic process, characterized by a persistence time and an effective temperature, and particles interact via a Lennard-Jones potential that yields well-studied glassy behavior in the Brownian limit, obtained as the persistence time vanishes. By increasing the persistence time, the system departs more strongly from thermal equilibrium and we provide a comprehensive numerical analysis of the structure and dynamics of the resulting active fluid. Finite persistence times profoundly affect the static structure of the fluid and give rise to nonequilibrium velocity correlations that are absent in thermal systems. Despite these nonequilibrium features, for any value of the persistence time we observe a nonequilibrium glass transition as the effective temperature is decreased. Surprisingly, increasing departure from thermal equilibrium is found to promote (rather than suppress) the glassy dynamics. Overall, our results suggest that with increasing persistence time, microscopic properties of the active fluid change quantitatively, but the broad features of the nonequilibrium glassy dynamics observed with decreasing the effective temperature remain qualitatively similar to those of thermal glass-formers.

preprint2015arXiv

Diverging viscosity and soft granular rheology in non-Brownian suspensions

We use large scale computer simulations and finite size scaling analysis to study the shear rheology of dense three-dimensional suspensions of frictionless non-Brownian particles in the vicinity of the jamming transition. We perform simulations of soft repulsive particles at constant shear rate, constant pressure, and finite system size, and study carefully the asymptotic limits of large system sizes and infinitely hard particle repulsion. Extending earlier analysis by about two orders of magnitude, we first study the asymptotic behavior of the shear viscosity in the hard particle limit. We confirm its asymptotic power law divergence at the jamming transition, but show that a precise determination of the critical density and critical exponent is difficult due to the `multiscaling' behavior of the viscosity. Additionally, finite-size scaling analysis suggests that this divergence is accompanied by a growing correlation length scale, which also diverges algebraically. We then study the effect of soft repulsion, and propose a natural extension of the standard granular rheology to account for softness effects, which we validate from simulations. Close to the jamming transition, this `soft granular rheology' offers a detailed description of the non-linear rheology of soft particles, which differs from earlier empirical scaling forms.

preprint2015arXiv

Evidence for a disordered critical point in a glass-forming liquid

Using computer simulations of an atomistic glass-forming liquid, we investigate the fluctuations of the overlap between a fluid configuration and a quenched reference system. We find that large fluctuations of the overlap develop as temperature decreases, consistent with the existence of the random critical point that is predicted by effective field theories. We discuss the scaling of fluctuations near the presumed critical point, comparing the observed behaviour with that of the random-field Ising model. We argue that this critical point directly reveals the existence of an interfacial tension between amorphous metastable states, a quantity relevant both for equilibrium relaxation and for nonequilibrium melting of stable glass configurations.

preprint2015arXiv

From single-particle to collective effective temperatures in an active fluid of self-propelled particles

We present a comprehensive analysis of effective temperatures based on fluctuation-dissipation relations in a model of an active fluid composed of self-propelled hard disks. We first investigate the relevance of effective temperatures in the dilute and moderately dense fluids. We find that a unique effective temperature does not in general characterize the non-equilibrium dynamics of the active fluid over this broad range of densities, because fluctuation-dissipation relations yield a lengthscale-dependent effective temperature. By contrast, we find that the approach to a non-equilibrium glass transition at very large densities is accompanied by the emergence of a unique effective temperature shared by fluctuations at all lengthscales. This suggests that an effective thermal dynamics generically emerges at long times in very dense suspensions of active particles due to the collective freezing occurring at non-equilibrium glass transitions.

preprint2015arXiv

Glassy dynamics of athermal self-propelled particles: Computer simulations and a nonequilibrium microscopic theory

We combine computer simulations and analytical theory to investigate the glassy dynamics in dense assemblies of athermal particles evolving under the sole influence of self-propulsion. The simulations reveal that when the persistence time of the self-propelled particles is increased, the local structure becomes more pronounced whereas the long-time dynamics first accelerates and then slows down. These seemingly contradictory evolutions are explained by constructing a nonequilibrium mode-coupling-like theory for interacting self-propelled particles. To predict the collective dynamics the theory needs the steady state structure factor and the steady state correlations of the local velocities. It yields nontrivial predictions for the glassy dynamics of self-propelled particles in qualitative agreement with the simulations.

preprint2015arXiv

Hyperuniform density fluctuations and diverging dynamic correlations in periodically driven colloidal suspensions

The emergence of particle irreversibility in periodically driven colloidal suspensions has been interpreted as resulting either from a nonequilibrium phase transition to an absorbing state or from the chaotic nature of particle trajectories. Using a simple model of a driven suspension we show that a nonequilibrium phase transition is accompanied by hyperuniform static density fluctuations in the vicinity of the transition, where we also observe strong dynamic heterogeneities reminiscent of dynamics in glassy materials. We find that single particle dynamics becomes intermittent and strongly non-Fickian, and that collective dynamics becomes spatially correlated over diverging lengthscales. Our results suggest that the two theoretical scenarii can be experimentally discriminated using particle-resolved measurements of standard static and dynamic observables.

preprint2015arXiv

Relaxation dynamics in a transient network fluid with competing gel and glass phases

We use computer simulations to study the relaxation dynamics of a model for oil-in-water microemulsion droplets linked with telechelic polymers. This system exhibits both gel and glass phases and we show that the competition between these two arrest mechanisms can result in a complex, three-step decay of the time correlation functions, controlled by two different localization lengthscales. For certain combinations of the parameters, this competition gives rise to an anomalous logarithmic decay of the correlation functions and a subdiffusive particle motion, which can be understood as a simple crossover effect between the two relaxation processes. We establish a simple criterion for this logarithmic decay to be observed. We also find a further logarithmically slow relaxation related to the relaxation of floppy clusters of particles in a crowded environment, in agreement with recent findings in other models for dense chemical gels. Finally, we characterize how the competition of gel and glass arrest mechanisms affects the dynamical heterogeneities and show that for certain combination of parameters these heterogeneities can be unusually large. By measuring the four-point dynamical susceptibility, we probe the cooperativity of the motion and find that with increasing coupling this cooperativity shows a maximum before it decreases again, indicating the change in the nature of the relaxation dynamics. Our results suggest that compressing gels to large densities produces novel arrested phases that have a new and complex dynamics.

preprint2015arXiv

Structure and dynamics of coupled viscous liquids

We perform Monte-Carlo simulations to analyse the structure and microscopic dynamics of a viscous Lennard-Jones liquid coupled to a quenched reference configuration of the same liquid. The coupling between the two replicas is introduced via a field epsilon conjugate to the overlap Q between the two particle configurations. This allows us to study the evolution of various static and dynamic correlation functions across the (epsilon, T) equilibrium phase diagram. As the temperature is decreased, we identify increasingly marked precursors of a first-order phase transition between a low-Q and a high-Q phase induced by the field epsilon. We show in particular that both static and dynamic susceptibilities have a maximum at a temperature-dependent value of the coupling field, which defines a `Widom line'. We also show that, in the high-overlap regime, diffusion and structural relaxation are strongly decoupled because single particle motion mostly occurs via discrete hopping on the sites defined by the reference configuration. These results, obtained using conventional numerical tools, provide encouraging signs that an equilibrium phase transition exists in coupled viscous liquids, but also demonstrate that important numerical challenges must be overcome to obtain more conclusive numerical evidence.

preprint2015arXiv

Thermal fluctuations, mechanical response, and hyperuniformity in jammed solids

Jamming is a geometric phase transition occurring in dense particle systems in the absence of temperature. We use computer simulations to analyse the effect of thermal fluctuations on several signatures of the transition. We show that scaling laws for bulk and shear moduli only become relevant when thermal fluctuations are extremely small, and propose their relative ratio as a quantitative signature of jamming criticality. Despite the nonequilibrium nature of the transition, we find that thermally induced fluctuations and mechanical responses obey equilibrium fluctuation-dissipation relations near jamming, provided the appropriate fluctuating component of the particle displacements is analysed. This shows that mechanical moduli can be directly measured from particle positions in mechanically unperturbed packings, and suggests that the definition of a "nonequilibrium index" is unnecessary for amorphous materials. We find that fluctuations of particle displacements are spatially correlated, and define a transverse and a longitudinal correlation lengthscales which both diverge as the jamming transition is approached. We analyse the frozen component of density fluctuations and find that it displays signatures of nearly-hyperuniform behaviour at large lengthscales. This demonstrates that hyperuniformity in jammed packings is unrelated to a vanishing compressibility and explains why it appears remarkably robust against temperature and density variations. Differently from jamming criticality, obstacles preventing the observation of hyperuniformity in colloidal systems do not originate from thermal fluctuations.

preprint2014arXiv

A novel approach to numerical measurements of the configurational entropy in supercooled liquids

The configurational entropy is among the key observables to characterize experimentally the formation of a glass. Physically, it quantifies the multiplicity of metastable states in which an amorphous material can be found at a given temperature, and its temperature dependence provides a major thermodynamic signature of the glass transition, which is experimentally accessible. Measurements of the configurational entropy require, however, some approximations which have often led to ambiguities and contradictory results. Here we implement a novel numerical scheme to measure the configurational entropy Sigma(T) in supercooled liquids, using a direct determination of the free energy cost to localize the system within a single metastable state at temperature T. For two prototypical glass-forming liquids, we find that Sigma(T) disappears discontinuously above a temperature T_c, which is slightly lower than the usual estimate of the onset temperature for glassy dynamics. This observation is in good agreement with theoretical expectations, but contrasts sharply with alternative numerical methods. While the temperature dependence of Sigma(T) correlates with the glass fragility, we show that the validity of the Adam-Gibbs relation (relating configurational entropy to structural relaxation time) established in earlier numerical studies is smaller than previously thought, potentially resolving an important conflict between experiments and simulations.

preprint2014arXiv

Clustering and heterogeneous dynamics in a kinetic Monte-Carlo model of self-propelled hard disks

We introduce a kinetic Monte-Carlo model for self-propelled hard disks to capture with minimal ingredients the interplay between thermal fluctuations, excluded volume and self-propulsion in large assemblies of active particles. We analyze in detail the resulting (density, self-propulsion) nonequilibrium phase diagram over a broad range of parameters. We find that purely repulsive hard disks spontaneously aggregate into fractal clusters as self-propulsion is increased, and rationalize the evolution of the average cluster size by developing a kinetic model of reversible aggregation. As density is increased, the nonequilibrium clusters percolate to form a ramified structure reminiscent of a physical gel. We show that the addition of a finite amount of noise is needed to trigger a nonequilibrium phase separation, showing that demixing in active Brownian particles results from a delicate balance between noise, interparticle interactions and self-propulsion. We show that self-propulsion has a profound influence on the dynamics of the active fluid. We find that the diffusion constant has a nonmonotonic behaviour as self-propulsion is increased at finite density and that activity produces strong deviations from Fickian diffusion that persist over large time scales and length scales, suggesting that systems of active particles generically behave as dynamically heterogeneous systems.

preprint2014arXiv

Crossovers in the dynamics of supercooled liquids probed by an amorphous wall

We study the relaxation dynamics of a binary Lennard-Jones liquid in the presence of an amorphous wall generated from equilibrium particle configurations. In qualitative agreement with the results presented in Nature Phys. {\bf 8}, 164 (2012) for a liquid of harmonic spheres, we find that our binary mixture shows a saturation of the dynamical length scale close to the mode-coupling temperature $T_c$. Furthermore we show that, due to the broken symmetry imposed by the wall, signatures of an additional change in dynamics become apparent at a temperature well above $T_c$. We provide evidence that this modification in the relaxation dynamics occurs at a recently proposed dynamical crossover temperature $T_s > T_c$, which is related to the breakdown of the Stokes-Einstein relation. We find that this dynamical crossover at $T_s$ is also observed for a system of harmonic spheres as well as a WCA liquid, showing that it may be a general feature of glass-forming systems.

preprint2014arXiv

Equilibrium ultrastable glasses produced by random pinning

Ultrastable glasses have risen to prominence due to their potentially useful material properties and the tantalizing possibility of a general method of preparation via vapor deposition. Despite the importance of this novel class of amorphous materials, numerical studies have been scarce because achieving ultrastability in atomistic simulations is an enormous challenge. Here we bypass this difficulty and establish that randomly pinning the position of a small fraction of particles inside an equilibrated supercooled liquid generates ultrastable configurations at essentially no numerical cost, while avoiding undesired structural changes due to the preparation protocol. Building on the analogy with vapor-deposited ultrastable glasses, we study the melting kinetics of these configurations following a sudden temperature jump into the liquid phase. In homogeneous geometries, we find that enhanced kinetic stability is accompanied by large scale dynamic heterogeneity, while a competition between homogeneous and heterogeneous melting is observed when a liquid boundary invades the glass at constant velocity. Our work demonstrates the feasibility of large-scale, atomistically resolved, and experimentally relevant simulations of the kinetics of ultrastable glasses.

preprint2014arXiv

Nonequilibrium equation of state in suspensions of active colloids

Active colloids constitute a novel class of materials composed of colloidal-scale particles locally converting chemical energy into motility, mimicking micro-organisms. Evolving far from equilibrium, these systems display structural organizations and dynamical properties distinct from thermalized colloidal assemblies. Harvesting the potential of this new class of systems requires the development of a conceptual framework to describe these intrinsically nonequilibrium systems. We use sedimentation experiments to probe the nonequilibrium equation of state of a bidimensional assembly of active Janus microspheres, and conduct computer simulations of a model of self-propelled hard disks. Self-propulsion profoundly affects the equation of state, but these changes can be rationalized using equilibrium concepts. We show that active colloids behave, in the dilute limit, as an ideal gas with an activity-dependent effective temperature. At finite density, increasing the activity is similar to increasing adhesion between equilibrium particles. We quantify this effective adhesion and obtain a unique scaling law relating activity and effective adhesion in both experiments and simulations. Our results provide a new and efficient way to understand the emergence of novel phases of matter in active colloidal suspensions.

preprint2014arXiv

Nonequilibrium glassy dynamics of self-propelled hard disks

We analyse the collective dynamics of self-propelled particles in the large density regime where passive particles undergo a kinetic arrest to an amorphous glassy state. We capture the competition between self-propulsion and crowding effects using a two-dimensional model of self-propelled hard disks, which we study using Monte-Carlo simulations. Although the activity drives the system far from equilibrium, self-propelled particles undergo a kinetic arrest, which we characterize in detail and compare with its equilibrium counterpart. In particular, the critical density for dynamic arrest continuously shifts to larger density with increasing activity, and the relaxation time is surprisingly well described by an algebraic divergence resulting from the emergence of highly collective dynamics. These results show that dense assemblies of active particles undergo a nonequilibrium glass transition which is profoundly affected by self-propulsion mechanisms.

preprint2014arXiv

Thinning or thickening? Multiple rheological regimes in dense suspensions of soft particles

The shear rheology of dense colloidal and granular suspensions is strongly nonlinear, as these materials exhibit shear-thinning and shear-thickening, depending on multiple physical parameters. We numerically study the rheology of a simple model of soft repulsive particles at large densities, and show that nonlinear flow curves reminiscent of experiments on real suspensions can be obtained. By using dimensional analysis and basic elements of kinetic theory, we rationalize these multiple rheological regimes and disentangle the relative impact of thermal fluctuations, glass and jamming transitions, inertia and particle softness on the flow curves. We characterize more specifically the shear-thickening regime and show that both particle softness and the emergence of a yield stress at the jamming transition compete with the inertial effects responsible for the observed thickening behaviour. This allows us to construct a dynamic state diagram, which can be used to analyze experiments.

preprint2013arXiv

Disentangling glass and jamming physics in the rheology of soft materials

The shear rheology of soft particles systems becomes complex at large density because crowding effects may induce a glass transition for Brownian particles, or a jamming transition for non-Brownian systems. Here we successfully explore the hypothesis that the shear stress contributions from glass and jamming physics are `additive'. We show that the experimental flow curves measured in a large variety of soft materials (colloidal hard spheres, microgel suspensions, emulsions, aqueous foams) as well as numerical flow curves obtained for soft repulsive particles in both thermal and athermal limits are well described by a simple model assuming that glass and jamming rheologies contribute linearly to the shear stress, provided that the relevant scales for time and stress are correctly identified in both sectors. Our analysis confirms that the dynamics of colloidal hard spheres is uniquely controlled by glass physics while aqueous foams are only sensitive to jamming effects. We show that for micron-sized emulsions both contributions are needed to successfully account for the flow curves, which reveal distinct signatures of both phenomena. Finally, for two systems of soft microgel particles we show that the flow curves are representative of the glass transition of colloidal systems, and deduce that microgel particles are not well suited to studying the jamming transition experimentally.

preprint2013arXiv

Intermittent dynamics and logarithmic domain growth during the spinodal decomposition of a glass-forming liquid

We use large-scale molecular dynamics simulations of a simple glass-forming system to investigate how its liquid-gas phase separation kinetics depends on temperature. A shallow quench leads to a fully demixed liquid-gas system whereas a deep quench makes the dense phase undergo a glass transition and become an amorphous solid. This glass has a gel-like bicontinuous structure that evolves very slowly with time and becomes fully arrested in the limit where thermal fluctuations become negligible. We show that the phase separation kinetics changes qualitatively with temperature, the microscopic dynamics evolving from a surface tension-driven diffusive motion at high temperature to a strongly intermittent, heterogeneous and thermally activated dynamics at low temperature, with a logarithmically slow growth of the typical domain size. These results shed light on recent experimental observations of various porous materials produced by arrested spinodal decomposition, such as nonequilibrium colloidal gels and bicontinuous polymeric structures, and they elucidate the microscopic mechanisms underlying a specific class of viscoelastic phase separation.

preprint2013arXiv

Nonequilibrium glass transitions in driven and active matter

The glass transition, extensively studied in dense fluids, polymers, or colloids, corresponds to a dramatic evolution of equilibrium transport coefficients upon a modest change of control parameter, like temperature or pressure. A similar phenomenology is found in many systems evolving far from equilibrium, such as driven granular media, active and living matter. While many theories compete to describe the glass transition at thermal equilibrium, very little is understood far from equilibrium. Here, we solve the dynamics of a specific, yet representative, class of glass models in the presence of nonthermal driving forces and energy dissipation, and show that a dynamic arrest can take place in these nonequilibrium conditions. While the location of the transition depends on the specifics of the driving mechanisms, important features of the glassy dynamics are insensitive to details, suggesting that an `effective' thermal dynamics generically emerges at long time scales in nonequilibrium systems close to dynamic arrest.

preprint2013arXiv

Overlap fluctuations in glass-forming liquids

We analyse numerically thermal fluctuations of the static overlap between equilibrium configurations in a glass-forming liquid approaching the glass transition. We find that the emergence of slow dynamics near the onset temperature correlates with the development of non-Gaussian probability distributions of overlap fluctuations, measured using both annealed and quenched definitions. Below a critical temperature, a thermodynamic field conjugate to the overlap induces a first-order phase transition, whose existence we numerically demonstrate in the annealed case. These results establish that the approach to the glass transition is accompanied by profound changes in the nature of thermodynamic fluctuations, deconstructing the view that glassy dynamics occurs with little structural evolution.

preprint2013arXiv

Probing a liquid to glass transition in equilibrium

We use computer simulations to investigate the static properties of a simple glass-forming fluid in which the positions of a finite fraction of the particles has been frozen in. By probing the equilibrium distribution of the overlap between independent configurations of the liquid, we find strong evidence that this random pinning induces a glass transition. At low temperatures, our numerical findings are consistent with the existence of a random first order phase transition rounded by finite size effects.

preprint2013arXiv

Yield stress in amorphous solids: A mode-coupling theory analysis

The yield stress is a defining feature of amorphous materials which is difficult to analyze theoretically, because it stems from the strongly non-linear response of an arrested solid to an applied deformation. Mode-coupling theory predicts the flow curves of materials undergoing a glass transition, and thus offers predictions for the yield stress of amorphous solids. We use this approach to analyse several classes of disordered solids, using simple models of hard sphere glasses, soft glasses, and metallic glasses for which the mode-coupling predictions can be directly compared to the outcome of numerical measurements. The theory correctly describes the emergence of a yield stress of entropic nature in hard sphere glasses, and its rapid growth as density approaches random close packing at qualitative level. By contrast, the emergence of solid behavior in soft and metallic glasses, which originates from direct particle interactions is not well described by the theory. We show that similar shortcomings arise in the description of the vibrational dynamics of the glass phase at rest. We discuss the range of applicability of mode-coupling theory to understand the yield stress and non-linear rheology of amorphous materials.

preprint2012arXiv

Dynamic criticality at the jamming transition

We characterize vibrational motion occurring at low temperatures in dense suspensions of soft repulsive spheres over a broad range of volume fractions encompassing the jamming transition at (T = 0, phi = phi_J). We find that characteristic time and length scales of thermal vibrations obey critical scaling in the vicinity of the jamming transition. We show in particular that the amplitude and the time scale of dynamic fluctuations diverge symmetrically on both sides of the transition, and directly reveal a diverging correlation length. The critical region near phi_J is divided in three different regimes separated by a characteristic temperature scale T*(phi) that vanishes quadratically with the distance to phi_J. While two of them, (T < T*(phi), phi > phi_J) and (T < T*(phi), phi < phi_J), are described by harmonic theories developed in the zero temperature limit, the third one for T > T*(phi) is inherently anharmonic and displays new critical properties. We find that the quadratic scaling of T*(phi) is due to nonperturbative anharmonic contributions, its amplitude being orders of magnitude smaller than the perturbative prediction based on the expansion to quartic order in the interactions. Our results show that thermal vibrations in colloidal assemblies directly reveal the critical nature of the jamming transition. The critical region, however, is very narrow and has not yet been attained experimentally, even in recent specifically-dedicated experiments.

preprint2012arXiv

Finite size effects in the dynamics of glass-forming liquids

We present a comprehensive theoretical study of finite size effects in the relaxation dynamics of glass-forming liquids. Our analysis is motivated by recent theoretical progress regarding the understanding of relevant correlation length scales in liquids approaching the glass transition. We obtain predictions both from general theoretical arguments and from a variety of specific perspectives: mode-coupling theory, kinetically constrained and defect models, and random first order transition theory. In the latter approach, we predict in particular a non-monotonic evolution of finite size effects across the mode-coupling crossover due to the competition between mode-coupling and activated relaxation. We study the role of competing relaxation mechanisms in giving rise to non-monotonic finite size effects by devising a kinetically constrained model where the proximity to the mode-coupling singularity can be continuously tuned by changing the lattice topology. We use our theoretical findings to interpret the results of extensive molecular dynamics studies of four model liquids with distinct structures and kinetic fragilities. While the less fragile model only displays modest finite size effects, we find a more significant size dependence evolving with temperature for more fragile models, such as Lennard-Jones particles and soft spheres. Finally, for a binary mixture of harmonic spheres we observe the predicted non-monotonic temperature evolution of finite size effects near the fitted mode-coupling singularity, suggesting that the crossover from mode-coupling to activated dynamics is more pronounced for this model. Finally, we discuss the close connection between our results and the recent report of a non-monotonic temperature evolution of a dynamic length scale near the mode-coupling crossover in harmonic spheres.

preprint2012arXiv

Spatial correlations in glass-forming liquids across the mode-coupling crossover

We discuss a novel approach, the point-to-set correlation functions, that allows to determine relevant static and dynamic length scales in glass-forming liquids. We find that static length scales increase monotonically when the temperature is lowered, whereas the measured dynamic length scale shows a maximum around the critical temperature of mode-coupling theory. We show that a similar non-monotonicity is found in the temperature evolution of certain finite size effects in the relaxation dynamics. These two independent sets of results demonstrate the existence of a change in the transport mechanism when the glass-former is cooled from moderately to deeply supercooled states across the mode-coupling crossover and clarify the status of the theoretical calculations done at the mean field level.

preprint2012arXiv

Unified study of glass and jamming rheology in soft particle systems

We explore numerically the shear rheology of soft repulsive particles at large volume fraction. The interplay between viscous dissipation and thermal motion results in multiple rheological regimes encompassing Newtonian, shear-thinning and yield stress regimes near the `colloidal' glass transition when thermal fluctuations are important, crossing over to qualitatively similar regimes near the `jamming' transition when dissipation dominates. In the crossover regime, glass and jamming sectors coexist and give complex flow curves. Although glass and jamming limits are characterized by similar macroscopic flow curves, we show that they occur over distinct time and stress scales and correspond to distinct microscopic dynamics. We propose a simple rheological model describing the glass to jamming crossover in the flow curves, and discuss the experimental implications of our results.

preprint2011arXiv

A microscopic mean-field theory of the jamming transition

Dense particle packings acquire rigidity through a nonequilibrium jamming transition commonly observed in materials from emulsions to sandpiles. We describe athermal packings and their observed geometric phase transitions using fully equilibrium statistical mechanics and develop a microscopic many-body mean-field theory of the jamming transition for soft repulsive spherical particles. We derive analytically some of the scaling laws and exponents characterizing the transition and obtain predictions for microscopic correlation functions of jammed states that are amenable to experimental verifications, and whose accuracy we confirm using computer simulations.

preprint2011arXiv

Highly nonlinear dynamics in a slowly sedimenting colloidal gel

We use a combination of original light scattering techniques and particles with unique optical properties to investigate the behavior of suspensions of attractive colloids under gravitational stress, following over time the concentration profile, the velocity profile, and the microscopic dynamics. During the compression regime, the sedimentation velocity grows nearly linearly with height, implying that the gel settling may be fully described by a (time-dependent) strain rate. We find that the microscopic dynamics exhibit remarkable scaling properties when time is normalized by strain rate, showing that the gel microscopic restructuring is dominated by its macroscopic deformation.

preprint2011arXiv

Influence of the glass transition on the liquid-gas spinodal decomposition

We use large-scale molecular dynamics simulations to study the kinetics of the liquid-gas phase separation if the temperature is lowered across the glass transition of the dense phase. We observe a gradual change from phase separated systems at high temperatures to nonequilibrium, gel-like structures that evolve very slowly at low temperatures. The microscopic mechanisms responsible for the coarsening strongly depend on temperature, and change from diffusive motion at high temperature to a strongly intermittent, heterogeneous and thermally activated dynamics at low temperature, leading to logarithmically slow growth of the typical domain size.

preprint2011arXiv

Inhomogeneous shear flows in soft jammed materials with tunable attractive forces

We perform molecular dynamics simulations to characterize the occurrence of inhomogeneous shear flows in soft jammed materials. We use rough walls to impose a simple shear flow and study the athermal motion of jammed assemblies of soft particles, both for purely repulsive interactions and in the presence of an additional short-range attraction of varying strength. In steady state, pronounced flow inhomogeneities emerge for all systems when the shear rate becomes small. Deviations from linear flow are stronger in magnitude and become very long-lived when the strength of the attraction increases, but differ from permanent shear-bands. Flow inhomogeneities occur in a stress window bounded by the dynamic and static yield stress values. Attractive forces enhance the flow heterogeneities because they accelerate stress relaxation, thus effectively moving the system closer to the yield stress regime where inhomogeneities are most pronounced. The present scenario for understanding the effect of particle adhesion on shear localization, which is based on detailed molecular dynamics simulations with realistic particle interactions, differs qualitatively from previous qualitative explanations and ad-hoc theoretical modelling.

preprint2011arXiv

Microscopic theory of the jamming transition of harmonic spheres

We develop a microscopic theory to analyze the phase behaviour and compute correlation functions of dense assemblies of soft repulsive particles both at finite temperature, as in colloidal materials, and at vanishing temperature, a situation relevant for granular materials and emulsions. We use a mean-field statistical mechanical approach which combines elements of liquid state theory to replica calculations to obtain quantitative predictions for the location of phase boundaries, macroscopic thermodynamic properties and microstructure of the system. We focus in particular on the derivation of scaling properties emerging in the vicinity of the jamming transition occurring at large density and zero temperature. The new predictions we obtain for pair correlation functions near contact are tested using computer simulations. Our work also clarifies the conceptual nature of the jamming transition, and its relation to the phenomenon of the glass transition observed in atomic liquids.

preprint2011arXiv

Non-monotonic temperature evolution of dynamic correlations in glass-forming liquids

The viscosity of glass-forming liquids increases by many orders of magnitude if their temperature is lowered by a mere factor of 2-3 [1,2]. Recent studies suggest that this widespread phenomenon is accompanied by spatially heterogeneous dynamics [3,4], and a growing dynamic correlation length quantifying the extent of correlated particle motion [5-7]. Here we use a novel numerical method to detect and quantify spatial correlations which reveal a surprising non-monotonic temperature evolution of spatial dynamical correlations, accompanied by a second length scale that grows monotonically and has a very different nature. Our results directly unveil a dramatic qualitative change in atomic motions near the mode-coupling crossover temperature [8] which involves no fitting or indirect theoretical interpretation. Our results impose severe new constraints on the theoretical description of the glass transition, and open several research perspectives, in particular for experiments, to confirm and quantify our observations in real materials.

preprint2011arXiv

Random pinning in glassy spin models with plaquette interactions

We use a random pinning procedure to study amorphous order in two glassy spin models. On increasing the concentration of pinned spins at constant temperature, we find a sharp crossover (but no thermodynamic phase transition) from bulk relaxation to localisation in a single state. At low temperatures, both models exhibit scaling behaviour. We discuss the growing length and time scales associated with amorphous order, and the fraction of pinned spins required to localize the system in a single state. These results, obtained for finite dimensional interacting models, provide a theoretical scenario for the effect of random pinning that differs qualitatively from previous approaches based either on mean-field, mode-coupling, or renormalization group reatments.

preprint2011arXiv

Static point-to-set correlations in glass-forming liquids

We analyze static point-to-set correlations in glass-forming liquids. The generic idea is to freeze the position of a set of particles in an equilibrium configuration and to perform sampling in the presence of this additional constraint. Qualitatively different geometries for the confining set of particles are considered and a detailed comparison of resulting static and dynamic correlation functions is performed. Our results reveal the existence of static spatial correlations not detected by conventional two-body correlators, which appear to be decoupled from, and shorter-ranged than, dynamical length scales characterizing dynamic heterogeneity. We find that the dynamics slows down dramatically under confinement, which suggests new ways to investigate the glass transition. Our results indicate that the geometry in which particles are randomly pinned is the best candidate to study static correlations.

preprint2011arXiv

Suppressed compressibility at large scale in jammed packings of size disperse spheres

We analyze the large scale structure and fluctuations of jammed packings of size disperse spheres, produced in a granular experiment as well as numerically. While the structure factor of the packings reveals no unusual behavior for small wavevectors, the compressibility displays an anomalous linear dependence at low wavectors and vanishes when q -> 0. We show that such behavior occurs because jammed packings of size disperse spheres have no bulk fluctuations of the volume fraction and are thus hyperuniform, a property not observed experimentally before. Our results apply to arbitrary particle size distributions. For continuous distributions, we derive a perturbative expression for the compressibility that is accurate for polydispersity up to about 30%.

preprint2011arXiv

Testing "microscopic" theories of glass-forming liquids

We assess the validity of "microscopic" approaches of glass-forming liquids based on the sole k nowledge of the static pair density correlations. To do so we apply them to a benchmark provided by two liquid models that share very similar static pair density correlation functions while disp laying distinct temperature evolutions of their relaxation times. We find that the approaches are unsuccessful in describing the difference in the dynamical behavior of the two models. Our study is not exhausti ve, and we have not tested the effect of adding corrections by including for instance three-body density correlations. Yet, our results appear strong enough to challenge the claim that the slowd own of relaxation in glass-forming liquids, for which it is well established that the changes of the static structure factor with temperature are small, can be explained by "microscopic" appr oaches only requiring the static pair density correlations as nontrivial input.

preprint2011arXiv

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 is 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.

preprint2011arXiv

Theoretical perspective on the glass transition and amorphous materials

We provide a theoretical perspective on the glass transition in molecular liquids at thermal equilibrium, on the spatially heterogeneous and aging dynamics of disordered materials, and on the rheology of soft glassy materials. We start with a broad introduction to the field and emphasize its connections with other subjects and its relevance. The important role played by computer simulations to study and understand the dynamics of systems close to the glass transition at the molecular level is spelled out. We review the recent progress on the subject of the spatially heterogeneous dynamics that characterizes structural relaxation in materials with slow dynamics. We then present the main theoretical approaches describing the glass transition in supercooled liquids, focusing on theories that have a microscopic, statistical mechanics basis. We describe both successes and failures, and critically assess the current status of each of these approaches. The physics of aging dynamics in disordered materials and the rheology of soft glassy materials are then discussed, and recent theoretical progress is described. For each section, we give an extensive overview of the most recent advances, but we also describe in some detail the important open problems that, we believe, will occupy a central place in this field in the coming years.

preprint2010arXiv

A critical test of the mode-coupling theory of the glass transition

The mode-coupling theory of the glass transition predicts the time evolution of the intermediate scattering functions in viscous liquids on the sole basis of the structural information encoded in two-point density correlations. We provide a critical test of this property and show that the theory fails to describe the qualitatively distinct dynamical behavior obtained in two model liquids characterized by very similar pair correlation functions. Because we use `exact' static information provided by numerical simulations, our results are a direct proof that some important information about the dynamics of viscous liquids is not captured by pair correlations, and is thus not described by the mode-coupling theory, even in the temperature regime where the theory is usually applied.

preprint2010arXiv

Can the jamming transition be described using equilibrium statistical mechanics?

When materials such as foams or emulsions are compressed, they display solid behaviour above the so-called `jamming' transition. Because compression is done out-of-equilibrium in the absence of thermal fluctuations, jamming appears as a new kind of a nonequilibrium phase transition. In this proceeding paper, we suggest that tools from equilibrium statistical mechanics can in fact be used to describe many specific features of the jamming transition. Our strategy is to introduce thermal fluctuations and use statistical mechanics to describe the complex phase behaviour of systems of soft repulsive particles, before sending temperature to zero at the end of the calculation. We show that currently available implementations of standard tools such as integral equations, mode-coupling theory, or replica calculations all break down at low temperature and large density, but we suggest that new analytical schemes can be developed to provide a fully microscopic, quantitative description of the jamming transition.

preprint2010arXiv

Equilibrium equation of state of a hard sphere binary mixture at very large densities using replica exchange Monte-Carlo simulations

We use replica exchange Monte-Carlo simulations to measure the equilibrium equation of state of the disordered fluid state for a binary hard sphere mixture up to very large densities where standard Monte-Carlo simulations do not easily reach thermal equilibrium. For the moderate system sizes we use (up to N=100), we find no sign of a pressure discontinuity near the location of dynamic glass singularities extrapolated using either algebraic or simple exponential divergences, suggesting they do not correspond to genuine thermodynamic glass transitions. Several scenarios are proposed for the fate of the fluid state in the thermodynamic limit.

preprint2010arXiv

Increasing the density melts ultrasoft colloidal glasses

We use theory and simulations to investigate the existence of amorphous glassy states in ultrasoft colloids. We combine the hyper-netted chain approximation with mode-coupling theory to study the dynamic phase diagram of soft repulsive spheres interacting with a Hertzian potential, focusing on low temperatures and large densities. At constant temperature, we find that an amorphous glassy state is entered upon compression, as in colloidal hard spheres, but the glass unexpectedly melts when density increases further. We attribute this re-entrant fluid-glass transition to particle softness, and correlate this behaviour to previously reported anomalies in soft systems, thus emphasizing its generality. The predicted fluid-glass-fluid sequence is confirmed numerically.

preprint2010arXiv

Overview of different characterisations of dynamic heterogeneity

Dynamic heterogeneity is now recognised as a central aspect of structural relaxation in disordered materials with slow dynamics, and was the focus of intense research in the last decade. Here we describe how initial, indirect observations of dynamic heterogeneity have recently evolved into well-defined, quantitative, statistical characterisations, in particular through the use of high-order correlation and response functions. We highlight both recent progress and open questions about the characterisation of dynamic heterogeneity in glassy materials. We also discuss the limits of available tools and describe a few candidates for future research in order to gain deeper understanding of the origin and nature of glassiness in disordered systems.

preprint2010arXiv

Subdiffusion and intermittent dynamic fluctuations in the aging regime of concentrated hard spheres

We study the nonequilibrium aging dynamics in a system of quasi-hard spheres at large density by means of computer simulations. We find that, after a sudden quench to large density, the relaxation time initially increases exponentially with the age of the system. After a surprisingly large crossover time, the system enters the asymptotic aging regime characterized by a linear increase of the relaxation time with age. In this aging regime, single particle motion is strongly non-Fickian, with a mean-squared displacement increasing subdiffusively, associated to broad, non-Gaussian tails in the distribution of particle displacements. We find that the system ages through temporally intermittent relaxation events, and a detailed finite size analysis of these collective dynamic fluctuations reveals that these events are not spanning the entire system, but remain spatially localized.

preprint2010arXiv

Superdiffusive, heterogeneous, and collective particle motion near the jamming transition in athermal disordered materials

We use computer simulations to study the microscopic dynamics of an athermal assembly of soft particles near the fluid-to-solid, jamming transition. Borrowing tools developed to study dynamic heterogeneity near glass transitions, we discover a number of original signatures of the jamming transition at the particle scale. We observe superdiffusive, spatially heterogeneous, and collective particle motion over a characteristic scale which displays a surprising non-monotonic behavior across the transition. In the solid phase, the dynamics is an intermittent succession of elastic deformations and plastic relaxations, which are both characterized by scale-free spatial correlations and system size dependent dynamic susceptibilities. Our results show that dynamic heterogeneities in dense athermal systems and glass-formers are very different, and shed light on recent experimental reports of `anomalous' dynamical behavior near the jamming transition of granular and colloidal assemblies.

preprint2009arXiv

Anomalous structural evolution of soft particles: Equibrium liquid state theory

We use the hyper-netted chain approximation of liquid state theory to analyze the evolution with density of the pair correlation function in a model of soft spheres with harmonic repulsion. As observed in recent experiments on jammed soft particles, theory predicts an `anomalous' (nonmonotonic) evolution of the intensity of the first peak when density is increased at constant temperature. This structural anomaly is a direct consequence of particle softness, and can be explained from purely equilibrium considerations, emphasizing the generality of the phenomenon. This anomaly is also predicted to have a non-trivial, `${\cal S}$-shaped', evolution with temperature, as a result of a competition between three distinct effects, which we describe in detail. Computer simulations support our predictions.

preprint2009arXiv

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 fraction where nonequilibrium jamming transitions occur. However, these jamming transitions share the same critical behaviour. Our results suggest that, even in the absence of partial crystalline ordering, a unique location of a random close packing does not exist, and that volume fraction alone is not sufficient to describe the properties of jammed states.

preprint2009arXiv

Nonperturbative effect of attractive forces in viscous liquids

We study the role of the attractive intermolecular forces in the viscous regime of a simple glass-forming liquid by using computer simulations. To this end, we compare the structure and the dynamics of a standard Lennard-Jones glass-forming liquid model with and without the attractive tail of the interaction potentials. The viscous slowing down of the two systems are found to be quantitatively and qualitatively different over a broad density range, whereas the static pair correlations remain close. The common assumption that the behaviour of dense nonassociated liquids is determined by the short-ranged repulsive part of the intermolecular potentials dramatically breaks down for the relaxation in the viscous liquid regime.

preprint2009arXiv

Scaling of the glassy dynamics of soft repulsive particles: a mode-coupling approach

We combine the hyper-netted chain approximation of liquid state theory with the mode-coupling theory of the glass transition to analyze the structure and dynamics of soft spheres interacting via harmonic repulsion. We determine the locus of the fluid-glass dynamic transition in a temperature -- volume fraction phase diagram. The zero-temperature (hard sphere) glass transition influences the dynamics at finite temperatures in its vicinity. This directly implies a form of dynamic scaling for both the average relaxation time and dynamic susceptibilities quantifying dynamic heterogeneity. We discuss several qualitative disagreements between theory and existing simulations at equilibrium. Our theoretical results are, however, very similar to numerical results for the driven athermal dynamics of repulsive spheres, suggesting that `mean-field' mode-coupling approaches might be good starting points to describe these nonequilibrium dynamics.

preprint2009arXiv

When gel and glass meet: A mechanism for multistep relaxation

We use computer simulations to study the dynamics of a physical gel at high densities where gelation and the glass transition interfere. We report and provide detailed physical understanding of complex relaxation patterns for time correlation functions which generically decay in a three-step process. For certain combinations of parameters we find logarithmic decays of the correlators and subdiffusive particle motion.

preprint2008arXiv

On the relaxation dynamics of glass-forming systems: Insights from computer simulations

We discuss the relaxation dynamics of a simple lattice gas model for glass-forming systems and show that with increasing density of particles this dynamics slows down very quickly. By monitoring the trajectory of tagged particles we find that their motion is very heterogeneous in space and time, leading to regions in space in which there is a fast dynamics and others in which it is slow. We determine how the geometric properties of these quickly relaxing regions depend on density and time. Motivated by this heterogeneous hopping dynamics, we use a simple model, a variant of a continuous time random walk, to characterize the relaxation dynamics. In particular we find from this model that for large displacements the self part of the van Hove function shows an exponential tail, in agreement with recent findings from experiments and simulations of glass-forming systems.

preprint2007arXiv

Heterogeneous diffusion in a reversible gel

We introduce a microscopically realistic model of a physical gel and use computer simulations to study its static and dynamic properties at thermal equilibrium. The phase diagram comprises a sol phase, a coexistence region ending at a critical point, a gelation line determined by geometric percolation, and an equilibrium gel phase unrelated to phase separation. The global structure of the gel is homogeneous, but the stress is supported by a fractal network. The gel dynamics is highly heterogeneous and we propose a theoretical model to quantitatively describe dynamic heterogeneity in gels. We elucidate several differences between the dynamics of gels and that of glass-formers.

preprint2006arXiv

The Monte-Carlo dynamics of a binary Lennard-Jones glass-forming mixture

We use a standard Monte-Carlo algorithm to study the slow dynamics of a binary Lennard-Jones glass-forming mixture at low temperature. We find that Monte-Carlo is by far the most efficient way to simulate a stochastic dynamics since relaxation is about 10 times faster than in Brownian Dynamics and about 30 times faster than in Stochastic Dynamics. Moreover, the average dynamical behaviour of the system is in quantitative agreement with the one obtained using Newtonian dynamics, apart at very short times where thermal vibrations are suppressed. We show, however, that dynamic fluctuations quantified by four-point dynamic susceptibilities do retain a dependence on the microscopic dynamics, as recently predicted theoretically.

preprint2005arXiv

Aging and ultra-slow equilibration in concentrated colloidal hard spheres

We study the dynamic behaviour of concentrated colloidal hard spheres using Time Resolved Correlation, a light scattering technique that can detect the slow evolution of the dynamics in out-of-equilibrium systems. Surprisingly, equilibrium is reached a very long time after sample initialization, the non-stationary regime lasting up to three orders of magnitude more than the relaxation time of the system. Before reaching equilibrium, the system displays unusual aging behaviour. The intermediate scattering function decays faster than exponentially and its relaxation time evolves non-monotonically with sample age.

preprint2004arXiv

What do we learn from the shape of the dynamical susceptibility of glass-formers?

We compute analytically and numerically the four-point correlation function that characterizes non-trivial cooperative dynamics in glassy systems within several models of glasses: elasto-plastic deformations, mode-coupling theory (MCT), collectively rearranging regions (CRR), diffusing defects and kinetically constrained models (KCM). Some features of the four-point susceptibility chi_4(t) are expected to be universal. at short times we expect an elastic regime characterized by a t or sqrt{t} growth. We find both in the beta, and the early alpha regime that chi_4 sim t^mu, where mu is directly related to the mechanism responsible for relaxation. This regime ends when a maximum of chi_4 is reached at a time t=t^* of the order of the relaxation time of the system. This maximum is followed by a fast decay to zero at large times. The height of the maximum also follows a power-law, chi_4(t^*) sim t^{*lambda}. The value of the exponents mu and lambda allows one to distinguish between different mechanisms. For example, freely diffusing defects in d=3 lead to mu=2 and lambda=1, whereas the CRR scenario rather predicts either mu=1 or a logarithmic behaviour depending on the nature of the nucleation events, and a logarithmic behaviour of chi_4(t^*). MCT leads to mu=b and lambda =1/gamma, where b and gamma are the standard MCT exponents. We compare our theoretical results with numerical simulations on a Lennard-Jones and a soft-sphere system. Within the limited time-scales accessible to numerical simulations, we find that the exponent mu is rather small, mu < 1, with a value in reasonable agreement with the MCT predictions.