Source author record

Atsushi Ikeda

Atsushi Ikeda 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

39works
5topics
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

39 published item(s)

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

Phonon transport properties of particulate physical gels

Particulate physical gels are sparse, low-density amorphous materials in which clusters of glasses are connected to form a heterogeneous network structure. This structure is characterized by two length scales, $ξ_s$ and $ξ_G$: $ξ_s$ measures the length of heterogeneities in the network structure, and $ξ_G$ is the size of glassy clusters. Accordingly, the vibrational states of such a material also exhibit a multiscale nature with two characteristic frequencies, $ω_\ast$ and $ω_G$, which are associated with $ξ_s$ and $ξ_G$, respectively: (i) phonon-like vibrations in the homogeneous medium at $ω< ω_\ast$, (ii) phonon-like vibrations in the heterogeneous medium at $ω_\ast < ω< ω_G$, and (iii) disordered vibrations in the glassy clusters at $ω> ω_G$. Here, we demonstrate that the multiscale characteristics seen in the static structures and vibrational states also extend to the phonon transport properties. Phonon transport exhibits two distinct crossovers at the frequencies $ω_\ast$ and $ω_G$~(or at wavenumbers of $\sim ξ_s^{-1}$ and $\sim ξ_G^{-1}$). In particular, both transverse and longitudinal phonons cross over between Rayleigh scattering at $ω< ω_\ast$ and diffusive damping at $ω>ω_\ast$. Remarkably, the Ioffe--Regel limit is located at the very low frequency of $ω_\ast$. Thus, phonon transport is localized above $ω_\ast$, even where phonon-like vibrational states persist. This markedly strong scattering behavior is caused by the sparse, porous structure of the gel.

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.

preprint2021arXiv

Computational simulations of the vibrational properties of glasses

Glasses show vibrational properties that are markedly different to those of crystals which are known as phonons. For example, excess low-frequency modes (the so-called boson peak), vibrational localization, and strong scattering of phonons have been the most discussed topics, and a theoretical understanding of these phenomena is challenging. To address this problem, computational simulations are a powerful tool, which have been employed by many previous works. In this chapter, we describe simulation methods for studying the vibrational properties of glasses (and any solid-state materials). We first present a method for studying vibrational eigenmodes. Since vibrational motions of particles are excited along eigenmodes, the eigenmodes are fundamental to descriptions of vibrational properties. The eigenmodes in glasses are non-phonon modes in general, and some of them are even localized in space. We next present a method of analysing phonon transport, which is also crucial for understanding vibrational properties. Since phonons are not eigenmodes in glasses, they are decomposed into several different, non-phonon eigenmodes. As a result, phonons in glasses are strongly scattered. In addition, we describe how to analyse the elastic response. The elastic response of glasses is also anomalous with respect to that of crystals. Finally, we briefly introduce recent advances that have been achieved by means of large-scale computational simulations.

preprint2021arXiv

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 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 new glass phase called the 2RSB double glass where both components vitrify simultaneously but with an energy landscape topography distinct from the 1RSB double glass.

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.

preprint2021arXiv

Structural, mechanical, and vibrational properties of particulate physical gels

Our lives are surrounded by a rich assortment of disordered materials. In particular, glasses are well known as dense, amorphous materials, whereas gels exist in low-density, disordered states. Recent progress has provided a significant step forward in understanding the material properties of glasses, such as mechanical, vibrational, and transport properties. In contrast, our understanding of particulate physical gels is still highly limited. Here, using molecular dynamics simulations, we study a simple model of particulate physical gels, the Lennard-Jones (LJ) gels, and provide a comprehensive understanding of their structural, mechanical, and vibrational properties, all of which are markedly different from those of glasses. First, the LJ gels show sparse, heterogeneous structures, and the length scale $ξ_s$ of the structures grows as the density is lowered. Second, the gels are extremely soft, with both shear $G$ and bulk $K$ moduli being orders of magnitude smaller than those of glasses. Third, many low-frequency vibrational modes are excited, which form a characteristic plateau with the onset frequency $ω_\ast$ in the vibrational density of states. Structural, mechanical, and vibrational properties, characterized by $ξ_s$, $G$, $K$, and $ω_\ast$, respectively, show power-law scaling behaviors with the density, which establishes a close relationship between them. Throughout the present work, we reveal that gels are multiscale, solid-state materials: (i) homogeneous elastic bodies at long lengths, (ii) heterogeneous elastic bodies with fractal structures at intermediate lengths, and (iii) amorphous structural bodies at short lengths.

preprint2020arXiv

Mechanical and Vibrational Properties of Three-Dimensional Dimer Packings Near the Jamming Transition

We comprehensively study mechanical and vibrational properties of dimer packings in three-dimensional space with particular attention on critical scaling behaviors near the jamming transition. First, we confirm the dependence of the packing fraction at the transition on the aspect ratio, the isostatic contact number at the transition, and the scaling dependence of the excess contact number on the excess density. Second, we study the elastic moduli, bulk and shear moduli, and establish power-law scaling of them. Finally, we study the vibrational density of states and its characteristic frequency. The vibrational density of states shows two plateaus in the lower- and higher-frequency regions, which are characterized by rotational and translational vibrational modes, respectively. The onset frequency of the lower-frequency plateau scales linearly to the square root of excess density. The scaling laws in the mechanical and vibrational properties are consistent between two- and three-dimensional dimers, and they are identical to those in spheres.

preprint2020arXiv

Novel elastic instability of amorphous solids in finite spatial dimensions

Recently, progress has been made in the understanding of anomalous vibrational excitations in amorphous solids. In the lowest-frequency region, the vibrational spectrum follows a non-Debye quartic law, which persists up to zero frequency without any frequency gap. This gapless vibrational density of states (vDOS) suggests that glasses are on the verge of instability. This feature of marginal stability is now highlighted as a key concept in the theories of glasses. In particular, the elasticity theory based on marginal stability predicts the gapless vDOS. However, this theory yields a quadratic law and \textit{not} the quartic law. To address this inconsistency, we presented a new type of instability, which is different from the conventional one, and proposed that amorphous solids are marginally stable considering the new instability in the preceding study~[M. Shimada, H. Mizuno, and A. Ikeda, Soft Matter, {\bf 16}, 7279, 2020]. In this study, we further extend and detail the results for these instabilities. By analyzing various examples of disorder, we demonstrate that real glasses in finite spatial dimensions can be marginally stable by the proposed novel instability.

preprint2020arXiv

Phase separation kinetics of a symmetric binary mixture of glass-forming liquids

Mixtures of glass-forming fluids sometimes exhibit glass-glass phase separation at low temperatures. Here, we use a molecular dynamics simulation to study one of the simplest examples of the glass-glass phase separation. We consider a mixture composed of type A and B particles, in which the A-A and B-B interactions are the identical Lennard-Jones interactions and the A-B interaction is repulsive only. To avoid crystallization, we also introduce the polydispersity in the particle sizes for each component. We study the phase separation kinetics of this model at a 50:50 concentration at various temperatures. We find that hydrodynamic coarsening takes place when the temperature is higher than the onset temperature of the glassy dynamics. At lower temperatures, diffusive coarsening is observed over a long duration, and a further slower coarsening appears within a shorter time. Below the glass transition temperature, the domain growth does not stop but becomes logarithmically slow or even slower than logarithmic. By analyzing two-time correlation functions, we show that these slow coarsening processes are accompanied by a slowing down of the microscopic dynamics, which has qualitative similarities with the aging dynamics without phase separation. Based on the results, we discuss a possible link between the slow coarsening and the aging-like microscopic slowing down in the glass-glass phase separation.

preprint2020arXiv

Spatial structure of unstable normal modes in a glass-forming liquid

The phenomenology of glass-forming liquids is often described in terms of their underlying, high-dimensional potential energy surface. In particular, the statistics of stationary points sampled as a function of temperature provides useful insight into the thermodynamics and dynamics of the system. To make contact with the real space physics, however, analysis of the spatial structure of the normal modes is required. In this work, we numerically study the potential energy surface of a glass-forming ternary mixture. Starting from liquid configurations equilibrated over a broad range of temperatures using a swap Monte Carlo method, we locate the nearby stationary points and investigate the spatial architecture and the energetics of the associated unstable modes. Through this spatially-resolved analysis, originally developed to study local minima, we corroborate recent evidence that the nature of the unstable modes changes from delocalized to localized around the mode-coupling temperature. We find that the displacement amplitudes of the delocalized modes have a slowly decaying far field, whereas the localized modes consist of a core with large displacements and a rapidly decaying far field. The fractal dimension of unstable modes around the mobility edge is equal to 1, consistent with the scaling of the participation ratio. Finally, we find that around and below the mode-coupling temperature the unstable modes are localized around structural defects, characterized by a disordered local structure markedly different from the liquid's locally favored structure. These defects are similar to those associated to quasi-localized vibrations in local minima and are good candidates to predict the emergence of localized excitations at low temperature.

preprint2019arXiv

Anharmonic properties of vibrational excitations in amorphous solids

Understanding the vibrational and thermal properties of amorphous solids is one of the most discussed and long-standing issues in condensed matter physics. Recent works have made significant steps towards understanding harmonic vibrational states. In particular, it has been established that quasi-localized vibrational modes emerge in addition to phonon-like vibrational modes. In this work, we study the anharmonic properties of these vibrational modes. We find that vibrational modes exhibit anharmonicities that induce particle rearrangements and cause transitions to different states. These anharmonicities are distinct from those in crystals, where particle rearrangements never occur. Remarkably, for both the phonon modes and quasi-localized modes, the vibrational modes exhibit strong anharmonicities, and the induced particle rearrangements are always localized in space and are composed of 1 to 1000 particles. Our findings contribute to the understanding of low-temperature thermal properties, for which anharmonic vibrations are crucial.

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

Stress relaxation above and below the jamming transition

We numerically investigate stress relaxation in soft athermal disks to reveal critical slowing down when the system approaches the jamming point. The exponents describing the divergence of the relaxation time differ dramatically depending on whether the transition is approached from the jammed or unjammed phase. This contrasts sharply with conventional dynamic critical scaling scenarios, where a single exponent characterizes both sides. We explain this surprising difference in terms of the vibrational density of states (vDOS), which is a key ingredient of linear viscoelastic theory. The vDOS exhibits an extra slow mode that emerges below jamming, which we utilize to demonstrate the anomalous exponent below jamming.

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.

preprint2019arXiv

Vibrational spectrum derived from the local mechanical response in disordered solids

The low-frequency vibrations of glasses are markedly different from those of crystals. These vibrations have recently been categorized into two types: spatially extended vibrations, whose vibrational density of states (vDOS) follows a non-Debye quadratic law, and quasilocalized vibrations (QLVs), whose vDOS follows a quartic law. The former are explained by elasticity theory with quenched disorder and microscopic replica theory as being a consequence of elastic instability, but the origin of the latter is still debated. Here, we show that the latter can also be directly derived from elasticity theory with quenched disorder. We find another elastic instability that the theory encompasses but that has been overlooked so far, namely, the instability of the system against a local dipolar force. This instability gives rise to an additional contribution to the vDOS, and the spatial structure and energetics of the mode originating from this instability are consistent with those of the QLVs. Finally, we construct a model in which the additional contribution to the vDOS follows a quartic law.

preprint2016arXiv

Mean-field dynamic criticality and geometric transition in the Gaussian core model

We use molecular dynamics simulations to investigate dynamic heterogeneities and the potential energy landscape of the Gaussian core model (GCM). Despite the nearly Gaussian statistics of particles' displacements, the GCM exhibits giant dynamic heterogeneities close to the dynamic transition temperature. The divergence of the four-point susceptibility is quantitatively well described by the inhomogeneous version of the Mode-Coupling theory. Furthermore, the potential energy landscape of the GCM is characterized by a geometric transition and large energy barriers, as expected from the lack of activated, hopping dynamics. These observations demonstrate that all major features of mean-field dynamic criticality can be observed in a physically sound, three-dimensional model.

preprint2016arXiv

The decoupling of the glass transitions in the two-component $p$-spin spherical model

Binary mixtures of large and small particles with disparate size ratio exhibit a rich phenomenology at their glass transition points. In order to gain insights on such systems, we introduce and study a two-component version of the $p$-spin spherical spin glass model. We employ the replica method to calculate the free energy and the phase diagram. We show that when the strengths of the interactions of each component are not widely separated, the model has only one glass phase characterized by the conventional one-step replica symmetry breaking. However when the strengths of the interactions are well separated, the model has three glass phases depending on temperature and component ratio. One is the "single" glass phase in which only the spins of one component are frozen while the spins of the other component remain mobile. This phase is characterized by the one-step replica symmetry breaking. The second is the "double" glass phase obtained by cooling further the single glass phase, in which the spins of the remaining mobile component are also frozen. This phase is characterized by the two-step replica symmetry breaking. The third is also the "double" glass phase, which however is formed by the simultaneous freezing of the spins of both components at the same temperatures and is characterized by the one-step replica symmetry breaking. We discuss the implications of these results for the glass transitions of binary mixtures.

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

One-dimensional Kac model of dense amorphous hard spheres

We introduce a new model of hard spheres under confinement for the study of the glass and jamming transitions. The model is an one-dimensional chain of the $d$-dimensional boxes each of which contains the same number of hard spheres, and the particles in the boxes of the ends of the chain are quenched at their equilibrium positions. We focus on the infinite dimensional limit ($d \to \infty$) of the model and analytically compute the glass transition densities using the replica liquid theory. From the chain length dependence of the transition densities, we extract the characteristic length scales at the glass transition. The divergence of the lengths are characterized by the two exponents, $-1/4$ for the dynamical transition and $-1$ for the ideal glass transition, which are consistent with those of the $p$-spin mean-field spin glass model. We also show that the model is useful for the study of the growing length scale at the jamming transition.

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

Correlation of Local Order with Particle Mobility in Supercooled Liquids is Highly System Dependent

We investigate the connection between local structure and dynamical heterogeneity in supercooled liquids. Through the study of four different models we show that the correlation between a particle's mobility and the degree of local order in nearby regions is highly system dependent. Our results suggest that the correlation between local structure and dynamics is weak or absent in systems that conform well to the mean-field picture of glassy dynamics and strong in those that deviate from this paradigm. Finally, we investigate the role of order-agnostic point-to-set correlations and reveal that they provide similar information content to local structure measures, at least in the system where local order is most pronounced.

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

Cluster and reentrant anomalies of nearly Gaussian core particles

We study through integral equation theory and numerical simulations the structure and dynamics of fluids composed of ultrasoft, nearly Gaussian particles. Namely, we explore the fluid phase diagram of a model in which particles interact via the generalized exponential potential u(r)=εexp[-(r/σ)^n], with a softness exponent n slightly larger than 2. In addition to the well-known anomaly associated to reentrant melting, the structure and dynamics of the fluid display two additional anomalies, which are visible in the isothermal variation of the structure factor and diffusivity. These features are correlated to the appearance of dimers in the fluid phase and to the subsequent modification of the cluster structure upon compression. We corroborate these results through an analysis of the local minima of the potential energy surface, in which clusters appear as much tighter conglomerates of particles. We find that reentrant melting and clustering coexist for softness exponents ranging from 2^+ up to values relevant for the description of amphiphilic dendrimers, i.e., n=3.

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

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

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 about how a fluid turns 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, we find the cage to keep a non-trivial 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 here qualitatively supported and non-mean-field corrections are found to remain small on decreasing dimension, reconsideration of its implementation is needed for it to result in a coherent description of experimental observations.

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

Jamming Transition and Inherent Structures of Hard Spheres and Discs

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

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

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

preprint2011arXiv

Glass Transition of the Monodisperse Gaussian Core Model

We numerically study dynamical properties of the one-component Gaussian Core Model in the supercooled states. We find that nucleation is suppressed as density increases. Concomitantly the system exhibits glassy slow dynamics characterized by the two-step and stretched exponential relaxation of the density correlation as well as drastic increase of the relaxation time. It is found that violation of the Stokes-Einstein relation is weaker and the non-Gaussian parameter is smaller than typical model glass formers, implying weaker dynamic heterogeneities. Besides, agreement of simulation data with the prediction of mode-coupling theory is exceptionally good, indicating that the nature of slow dynamics of this ultra-soft particle fluid is mean-field-like. This fact may be understood as the consequences of multiple overlaps of the constituent particles at high densities.

preprint2011arXiv

Slow Dynamics of the High Density Gaussian Core Model

We numerically study crystal nucleation and glassy slow dynamics of the one-component Gaussian core model (GCM) at high densities. The nucleation rate at a fixed supersaturation is found to decrease as the density increases. At very high densities, the nucleation is not observed at all in the time window accessed by long molecular dynamics (MD) simulation. Concomitantly, the system exhibits typical slow dynamics of the supercooled fluids near the glass transition point. We compare the simulation results of the supercooled GCM with the predictions of mode-coupling theory (MCT) and find that the agreement between them is better than any other model glassformers studied numerically in the past. Furthermore, we find that a violation of the Stokes-Einstein relation is weaker and the non-Gaussian parameter is smaller than canonical glassformers. Analysis of the probability distribution of the particle displacement clearly reveals that the hopping effect is strongly suppressed in the high density GCM. We conclude from these observations that the GCM is more amenable to the mean-field picture of the glass transition than other models. This is attributed to the long-ranged nature of the interaction potential of the GCM in the high density regime. Finally, the intermediate scattering function at small wavevectors is found to decay much faster than its self part, indicating that dynamics of the large-scale density fluctuations decouples with the shorter-ranged caging motion.

preprint2011arXiv

Thermodynamics and Structural Properties of the High Density Gaussian Core Model

We numerically study thermodynamic and structural properties of the one-component Gaussian core model (GCM) at very high densities. The solid-fluid phase boundary is carefully determined. We find that the density dependence of both the freezing and melting temperatures obey the asymptotic relation, $\log T_f$, $\log T_m \propto -ρ^{2/3}$, where $ρ$ is the number density, which is consistent with Stillinger's conjecture. Thermodynamic quantities such as the energy and pressure and the structural functions such as the static structure factor are also investigated in the fluid phase for a wide range of temperature above the phase boundary. We compare the numerical results with the prediction of the liquid theory with the random phase approximation (RPA). At high temperatures, the results are in almost perfect agreement with RPA for a wide range of density, as it has been already shown in the previous studies. In the low temperature regime close to the phase boundary line, although RPA fails to describe the structure factors and the radial distribution functions at the length scales of the interparticle distance, it successfully predicts their behaviors at shorter length scales. RPA also predicts thermodynamic quantities such as the energy, pressure, and the temperature at which the thermal expansion coefficient becomes negative, almost perfectly. Striking ability of RPA to predict thermodynamic quantities even at high densities and low temperatures is understood in terms of the decoupling of the length scales which dictate thermodynamic quantities from the interparticle distance which dominates the peak structures of the static structure factor due to the softness of the Gaussian core potential.

preprint2010arXiv

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 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 dimension, 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.

preprint2010arXiv

The double cover of cubic surfaces branched along their Hessian

We prove the relation between the Hodge structure of the double cover of a nonsingular cubic surface branched along its Hessian and the Hodge structure of the triple cover of the ambient projective space branched along the cubic surface. And we introduce a method to study the infinitesimal variations of Hodge structure of the double cover of the cubic surface. Using these results, we compute the Néron-Severi lattices for the double cover of a generic cubic surface and the Fermat cubic surface.