1
|
Kawasaki T, Miyazaki K. Unified Understanding of Nonlinear Rheology near the Jamming Transition Point. PHYSICAL REVIEW LETTERS 2024; 132:268201. [PMID: 38996305 DOI: 10.1103/physrevlett.132.268201] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/01/2023] [Accepted: 05/09/2024] [Indexed: 07/14/2024]
Abstract
When slowly sheared, jammed packings respond elastically before yielding. This linear elastic regime becomes progressively narrower as the jamming transition point is approached, and rich nonlinear rheologies such as shear softening and hardening or melting emerge. However, the physical mechanism of these nonlinear rheologies remains elusive. To clarify this, we numerically study jammed packings of athermal frictionless soft particles under quasistatic shear γ. We find the universal scaling behavior for the ratio of the shear stress σ and the pressure P, independent of the preparation protocol of the initial configurations. In particular, we reveal shear softening σ/P∼γ^{1/2} over an unprecedentedly wide range of strain up to the yielding point, which a simple scaling argument can rationalize.
Collapse
|
2
|
Interiano-Alberto KA, Morse PK, Hoy RS. Critical-like slowdown in thermal soft-sphere glasses via energy minimization. Phys Rev E 2024; 109:L062603. [PMID: 39020966 DOI: 10.1103/physreve.109.l062603] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2023] [Accepted: 06/03/2024] [Indexed: 07/20/2024]
Abstract
Using hybrid molecular dynamics/SWAP Monte Carlo (MD/SMC) simulations, we show that while the terminal relaxation times τ(ϕ) for FIRE energy minimization of soft-sphere glasses can decrease by orders of magnitude as sample equilibration proceeds and the jamming density ϕ_{J} increases, they always scale as τ(ϕ)∼(ϕ_{J}-ϕ)^{-2}∼[Z_{iso}-Z_{ms}(τ)]^{-2}, where Z_{iso}=2d and Z_{ms}(τ) is the average coordination number of particles satisfying a minimal local mechanical stability criterion (Z≥d+1) at the top of the final potential-energy-landscape (PEL) sub-basin the system encounters. This scaling allows us to collapse τ datasets that look very different when plotted as a function of ϕ, and to address a closely related question: how does the character of the PEL basins that dense thermal glasses most typically occupy evolve as the glasses age at constant ϕ and T?
Collapse
|
3
|
Hedström L, Olsson P. Considerations on the relaxation time in shear-driven jamming. Phys Rev E 2024; 109:064904. [PMID: 39020915 DOI: 10.1103/physreve.109.064904] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2024] [Accepted: 05/08/2024] [Indexed: 07/20/2024]
Abstract
We study the jamming transition in a model of elastic particles under shear at zero temperature, with a focus on the relaxation time τ_{1}. This relaxation time is from two-step simulations where the first step is the ordinary shearing simulation and the second step is the relaxation of the energy after stopping the shearing. τ_{1} is determined from the final exponential decay of the energy. Such relaxations are done with many different starting configurations generated by a long shearing simulation in which the shear variable γ slowly increases. We study the correlations of both τ_{1}, determined from the decay, and the pressure, p_{1}, from the starting configurations as a function of the difference in γ. We find that the correlations of p_{1} are longer lived than the ones of τ_{1} and find that the reason for this is that the individual τ_{1} is controlled both by p_{1} of the starting configuration and a random contribution which depends on the relaxation path length-the average distance moved by the particles during the relaxation. We further conclude that it is γ_{τ}, determined from the correlations of τ_{1}, which is the relevant one when the aim is to generate data that may be used for determining the critical exponent that characterizes the jamming transition.
Collapse
|
4
|
Shivers JL, Sharma A, MacKintosh FC. Strain-Controlled Critical Slowing Down in the Rheology of Disordered Networks. PHYSICAL REVIEW LETTERS 2023; 131:178201. [PMID: 37955486 DOI: 10.1103/physrevlett.131.178201] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/24/2022] [Revised: 06/19/2023] [Accepted: 09/25/2023] [Indexed: 11/14/2023]
Abstract
Networks and dense suspensions frequently reside near a boundary between soft (or fluidlike) and rigid (or solidlike) regimes. Transitions between these regimes can be driven by changes in structure, density, or applied stress or strain. In general, near the onset or loss of rigidity in these systems, dissipation-limiting heterogeneous nonaffine rearrangements dominate the macroscopic viscoelastic response, giving rise to diverging relaxation times and power-law rheology. Here, we describe a simple quantitative relationship between nonaffinity and the excess viscosity. We test this nonaffinity-viscosity relationship computationally and demonstrate its rheological consequences in simulations of strained filament networks and dense suspensions. We also predict critical signatures in the rheology of semiflexible and stiff biopolymer networks near the strain stiffening transition.
Collapse
Affiliation(s)
- Jordan L Shivers
- Department of Chemical and Biomolecular Engineering, Rice University, Houston, Texas 77005, USA
- Center for Theoretical Biological Physics, Rice University, Houston, Texas 77005, USA
| | - Abhinav Sharma
- Institute of Physics, University of Augsburg, 86159 Augsburg, Germany
- Leibniz-Institut für Polymerforschung Dresden, Institut Theorie der Polymere, 01069 Dresden, Germany
| | - Fred C MacKintosh
- Department of Chemical and Biomolecular Engineering, Rice University, Houston, Texas 77005, USA
- Center for Theoretical Biological Physics, Rice University, Houston, Texas 77005, USA
- Department of Chemistry, Rice University, Houston, Texas 77005, USA
- Department of Physics and Astronomy, Rice University, Houston, Texas 77005, USA
| |
Collapse
|
5
|
Taghizadeh K, Ruf M, Luding S, Steeb H. X-ray 3D imaging-based microunderstanding of granular mixtures: Stiffness enhancement by adding small fractions of soft particles. Proc Natl Acad Sci U S A 2023; 120:e2219999120. [PMID: 37339218 PMCID: PMC10293847 DOI: 10.1073/pnas.2219999120] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2022] [Accepted: 05/21/2023] [Indexed: 06/22/2023] Open
Abstract
This research focuses on performing ultrasound propagation measurements and micro-X-ray computed tomography (µXRCT) imaging on prestressed granular packings prepared with biphasic mixtures of monodisperse glass and rubber particles at different compositions/fractions. Ultrasound experiments employing piezoelectric transducers, mounted in an oedometric cell (complementing earlier triaxial cell experiments), are used to excite and detect longitudinal ultrasound waves through randomly prepared mixtures of monodisperse stiff/soft particles. While the fraction of the soft particles is increasing linearly from zero, the effective macroscopic stiffness of the granular packings transits nonlinearly and nonmonotonically toward the soft limit, remarkably via an interesting stiffer regime for small rubber fractions between 0.1 ≲ ν ≲ 0.2. The contact network of dense packings, as accessed from µXRCT, plays a key role in understanding this phenomenon, considering the structure of the network, the chain length, the grain contacts, and the particle coordination. While the maximum stiffness is due to surprisingly shortened chains, the sudden drop in elastic stiffness of the mixture packings, at ν ≈ 0.4, is associated with chains of particles that include both glass and rubber particles (soft chains); for ν ≲ 0.3, the dominant chains include only glass particles (hard chains). At the drop, ν ≈ 0.4, the coordination number of glass and rubber networks is approximately four and three, respectively, i.e., neither of the networks are jammed, and the chains need to include particles from another species to propagate information.
Collapse
Affiliation(s)
- K. Taghizadeh
- Institute of Applied Mechanics, University of Stuttgart, Stuttgart70569, Germany
- Multi-Scale Mechanics, University of Twente, AE Enschede7500, The Netherlands
| | - M. Ruf
- Institute of Applied Mechanics, University of Stuttgart, Stuttgart70569, Germany
| | - S. Luding
- Multi-Scale Mechanics, University of Twente, AE Enschede7500, The Netherlands
| | - H. Steeb
- Institute of Applied Mechanics, University of Stuttgart, Stuttgart70569, Germany
- Stuttgart Center for Simulation Science, University of Stuttgart, Stuttgart70569, Germany
| |
Collapse
|
6
|
Barik S, Majumdar S. Origin of Two Distinct Stress Relaxation Regimes in Shear Jammed Dense Suspensions. PHYSICAL REVIEW LETTERS 2022; 128:258002. [PMID: 35802438 DOI: 10.1103/physrevlett.128.258002] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/01/2021] [Revised: 05/05/2022] [Accepted: 06/08/2022] [Indexed: 06/15/2023]
Abstract
Many dense particulate suspensions show a stress induced transformation from a liquidlike state to a solidlike shear jammed (SJ) state. However, the underlying particle-scale dynamics leading to such striking, reversible transition of the bulk remains unknown. Here, we study transient stress relaxation behaviour of SJ states formed by a well-characterized dense suspension under a step strain perturbation. We observe a strongly nonexponential relaxation that develops a sharp discontinuous stress drop at short time for high enough peak-stress values. High resolution boundary imaging and normal stress measurements confirm that such stress discontinuity originates from the localized plastic events, whereas system spanning dilation controls the slower relaxation process. We also find an intriguing correlation between the nature of transient relaxation and the steady-state shear jamming phase diagram obtained from the Wyart-Cates model.
Collapse
Affiliation(s)
- Sachidananda Barik
- Soft Condensed Matter Group, Raman Research Institute, Bangalore 560080, Karnataka, India
| | - Sayantan Majumdar
- Soft Condensed Matter Group, Raman Research Institute, Bangalore 560080, Karnataka, India
| |
Collapse
|
7
|
Olsson P. Relaxation times, rheology, and finite size effects for non-Brownian disks in two dimensions. Phys Rev E 2022; 105:034902. [PMID: 35428108 DOI: 10.1103/physreve.105.034902] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2021] [Accepted: 03/01/2022] [Indexed: 06/14/2023]
Abstract
We carry out overdamped simulations in a simple model of jamming-a collection of bidisperse soft core frictionless disks in two dimensions-with the aim to explore the finite size dependence of different quantities, both the relaxation time obtained from the relaxation of the energy and the pressure equivalent of the shear viscosity. The motivation for the paper is the observation [Nishikawa et al., J. Stat. Phys. 182, 37 (2021)0022-471510.1007/s10955-021-02710-8] that there are finite size effects in the relaxation time, τ, that give problems in the determination of the critical divergence, and the claim that this is due to a finite size dependence, τ∼lnN, which makes τ an ill-defined quantity. Beside analyses to determine the relaxation time for the whole system we determine particle relaxation times which allow us to determine both histograms of particle relaxation times and the average particle relaxation times-two quantities that are very useful for the analyses. The starting configurations for the relaxation simulations are of two different kinds-completely random or taken from steady shearing simulations-and we find that the difference between these two cases are bigger than previously noted and that the observed problems in the determination of the critical divergence obtained when starting from random configurations are not present when instead starting the relaxations from shearing configurations. We also argue that the the effect that causes the lnN dependence is not as problematic as asserted. When it comes to the finite size dependence of the pressure equivalent of the shear viscosity we find that our data don't give support for the claimed strong finite size dependence, but also that the finite size dependence is at odds with what one would normally expect for a system with a diverging correlation length, and that this calls for an alternative understanding of the phenomenon of shear-driven jamming.
Collapse
Affiliation(s)
- Peter Olsson
- Department of Physics, Umeå University, 901 87 Umeå, Sweden
| |
Collapse
|
8
|
Peshkov A, Teitel S. Universality of stress-anisotropic and stress-isotropic jamming of frictionless spheres in three dimensions: Uniaxial versus isotropic compression. Phys Rev E 2022; 105:024902. [PMID: 35291159 DOI: 10.1103/physreve.105.024902] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/08/2021] [Accepted: 01/28/2022] [Indexed: 06/14/2023]
Abstract
We numerically study a three-dimensional system of athermal, overdamped, frictionless spheres, using a simplified model for a non-Brownian suspension. We compute the bulk viscosity under both uniaxial and isotropic compression as a means to address the question of whether stress-anisotropic and stress-isotropic jamming are in the same critical universality class. Carrying out a critical scaling analysis of the system pressure p, shear stress σ, and macroscopic friction μ=σ/p, as functions of particle packing fraction ϕ and compression rate ε[over ̇], we find good agreement for all critical parameters comparing the isotropic and anisotropic cases. In particular, we determine that the bulk viscosity diverges as p/ε[over ̇]∼(ϕ_{J}-ϕ)^{-β}, with β=3.36±0.09, as jamming is approached from below. We further demonstrate that the average contact number per particle Z can also be written in a scaling form as a function of ϕ and ε[over ̇]. Once again, we find good agreement between the uniaxial and isotropic cases. We compare our results to prior simulations and theoretical predictions.
Collapse
Affiliation(s)
- Anton Peshkov
- Department of Physics and Astronomy, University of Rochester, Rochester, New York 14627, USA
| | - S Teitel
- Department of Physics and Astronomy, University of Rochester, Rochester, New York 14627, USA
| |
Collapse
|
9
|
Matsuyama H, Toyoda M, Kurahashi T, Ikeda A, Kawasaki T, Miyazaki K. Geometrical properties of mechanically annealed systems near the jamming transition. THE EUROPEAN PHYSICAL JOURNAL. E, SOFT MATTER 2021; 44:133. [PMID: 34718887 DOI: 10.1140/epje/s10189-021-00142-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/04/2021] [Accepted: 10/19/2021] [Indexed: 06/13/2023]
Abstract
Geometrical properties of two-dimensional mixtures near the jamming transition point are numerically investigated using harmonic particles under mechanical training. The configurations generated by the quasi-static compression and oscillatory shear deformations exhibit anomalous suppression of the density fluctuations, known as hyperuniformity, below and above the jamming transition. For the jammed system trained by compression above the transition point, the hyperuniformity exponent increases. For the system below the transition point under oscillatory shear, the hyperuniformity exponent also increases until the shear amplitude reaches the threshold value. The threshold value matches with the transition point from the point-reversible phase where the particles experience no collision to the loop-reversible phase where the particles' displacements are non-affine during a shear cycle before coming back to an original position. The results demonstrated in this paper are explained in terms of neither of universal criticality of the jamming transition nor the nonequilibrium phase transitions.
Collapse
Affiliation(s)
| | - Mari Toyoda
- Department of Physics, Nagoya University, Nagoya, 464-8602, Japan
| | - Takumi Kurahashi
- Department of Physics, Nagoya University, Nagoya, 464-8602, Japan
| | - Atsushi Ikeda
- Graduate School of Arts and Sciences, University of Tokyo, Tokyo, 153-8902, Japan
| | - Takeshi Kawasaki
- Department of Physics, Nagoya University, Nagoya, 464-8602, Japan
| | | |
Collapse
|
10
|
Saitoh K, Mizuno H. Sound damping in soft particle packings: the interplay between configurational disorder and inelasticity. SOFT MATTER 2021; 17:4204-4212. [PMID: 33881038 DOI: 10.1039/d0sm02018d] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
We numerically investigate sound damping in disordered two-dimensional soft particle packings. We simulate evolution of standing waves of particle displacements and analyze time correlation functions of particle velocities and power spectra. We control the strength of inelastic interactions between the particles in contact to show how the inelasticity affects anomalous sound characteristics of disordered systems: Increasing the strength of inelastic interactions, we find that (i) sound softening vanishes and (ii) attenuation coefficients exhibit a transition from the Rayleigh law to quadratic growth. We also report (iii) how the Ioffe-Regel limit frequencies depend on the strength of inelasticity as useful information for experiments and applications of the sound in disordered media. Our findings suggest that sound damping in soft particle packings is determined by the interplay between elastic heterogeneities and inelasticity.
Collapse
Affiliation(s)
- Kuniyasu Saitoh
- Department of Physics, Faculty of Science, Kyoto Sangyo University, Motoyama, Kamigamo, Kita-ku, Kyoto 603-8555, Japan.
| | - Hideyuki Mizuno
- Graduate School of Arts and Sciences, The University of Tokyo, Tokyo 153-8902, Japan
| |
Collapse
|
11
|
Ikeda H, Hukushima K. Nonaffine displacements below jamming under athermal quasistatic compression. Phys Rev E 2021; 103:032902. [PMID: 33862705 DOI: 10.1103/physreve.103.032902] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/28/2020] [Accepted: 02/23/2021] [Indexed: 11/07/2022]
Abstract
Critical properties of frictionless spherical particles below jamming are studied using extensive numerical simulations, paying particular attention to the nonaffine part of the displacements during the athermal quasistatic compression. It is shown that the squared norm of the nonaffine displacement exhibits a power-law divergence toward the jamming transition point. A possible connection between this critical exponent and that of the shear viscosity is discussed. The participation ratio of the displacements vanishes in the thermodynamic limit at the transition point, meaning that the nonaffine displacements are localized marginally with a fractal dimension. Furthermore, the distribution of the displacement is shown to have a power-law tail, the exponent of which is related to the fractal dimension.
Collapse
Affiliation(s)
- Harukuni Ikeda
- Graduate School of Arts and Sciences, The University of Tokyo 153-8902, Japan
| | - Koji Hukushima
- Graduate School of Arts and Sciences, The University of Tokyo 153-8902, Japan.,Komaba Institute for Science, The University of Tokyo, 3-8-1 Komaba, Meguro, Tokyo 153-8902, Japan
| |
Collapse
|
12
|
Peshkov A, Teitel S. Critical scaling of compression-driven jamming of athermal frictionless spheres in suspension. Phys Rev E 2021; 103:L040901. [PMID: 34006006 DOI: 10.1103/physreve.103.l040901] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2020] [Accepted: 03/08/2021] [Indexed: 11/07/2022]
Abstract
We study numerically a system of athermal, overdamped, frictionless spheres, as in a non-Brownian suspension, in two and three dimensions. Compressing the system isotropically at a fixed rate ε[over ̇], we investigate the critical behavior at the jamming transition. The finite compression rate introduces a control timescale, which allows one to probe the critical timescale associated with jamming. As was found previously for steady-state shear-driven jamming, we find for compression-driven jamming that pressure obeys a critical scaling relation as a function of packing fraction ϕ and compression rate ε[over ̇], and that the bulk viscosity p/ε[over ̇] diverges upon jamming. A scaling analysis determines the critical exponents associated with the compression-driven jamming transition. Our results suggest that stress-isotropic, compression-driven jamming may be in the same universality class as stress-anisotropic, shear-driven jamming.
Collapse
Affiliation(s)
- Anton Peshkov
- Department of Physics and Astronomy, University of Rochester, Rochester, New York 14627, USA
| | - S Teitel
- Department of Physics and Astronomy, University of Rochester, Rochester, New York 14627, USA
| |
Collapse
|
13
|
Ishikawa H, Takada S, Matsumoto Y. Rheology of two-dimensional crushable granular materials. EPJ WEB OF CONFERENCES 2021. [DOI: 10.1051/epjconf/202124907007] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The rheology of two-dimensional crushable granular materials under shear is numerically studied using the discrete element method. We find that the mean fragment size changes as the shear strain increases while the shear stress is almost independent of this mean size. The fragment size distribution is found to follow a power law. In particular, the exponent in the intermediate fragment size regime becomes approximately – 11/6, which is almost independent of the shear rate.
Collapse
|
14
|
Affiliation(s)
- Harukuni Ikeda
- Graduate School of Arts and Sciences, The University of Tokyo, Tokyo 153-8902, Japan
| |
Collapse
|
15
|
Hossain T, Rognon P. Mobility in immersed granular materials upon cyclic loading. Phys Rev E 2020; 102:022904. [PMID: 32942415 DOI: 10.1103/physreve.102.022904] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2020] [Accepted: 07/24/2020] [Indexed: 11/07/2022]
Abstract
We study the mobility of objects embedded in an immersed granular packing and subjected to cyclic loadings. With this aim, we conducted uplift experiments whereby a horizontal plate is embedded in the packing and subjected to a vertical cyclic force oscillating between zero and a maximum amplitude. Tests performed at different cyclic force frequencies and amplitudes evidence the development of three mobility regimes whereby the plate stays virtually immobile, moves up steadily, or slowly creeps upwards. Results show that steady plate uplift can occur at lower force magnitudes when the frequency is increased. We propose an interpretation of this frequency-weakening behavior based on force relaxation experiments and on the analysis of the mobility response of theoretical viscoelastoplastic mechanical analog. These results and analysis point out inherent differences in mobility response between steady and cyclic loadings in immersed granular materials.
Collapse
Affiliation(s)
- Tanvir Hossain
- Particles and Grains Laboratory, School of Civil Engineering, University of Sydney, Sydney, NSW 2006, Australia
| | - Pierre Rognon
- Particles and Grains Laboratory, School of Civil Engineering, University of Sydney, Sydney, NSW 2006, Australia
| |
Collapse
|
16
|
Hwang S, Ikeda H. Force balance controls the relaxation time of the gradient descent algorithm in the satisfiable phase. Phys Rev E 2020; 101:052308. [PMID: 32575263 DOI: 10.1103/physreve.101.052308] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2019] [Accepted: 04/15/2020] [Indexed: 06/11/2023]
Abstract
We numerically study the relaxation dynamics of the single-layer perceptron with the spherical constraint. This is the simplest model of neural networks and serves a prototypical mean-field model of both convex and nonconvex optimization problems. The relaxation time of the gradient descent algorithm rapidly increases near the SAT-UNSAT (satisfiable-unsatisfiable) transition point. We numerically confirm that the first nonzero eigenvalue of the Hessian controls the relaxation time. This first eigenvalue vanishes much faster upon approaching the SAT-UNSAT transition point than the prediction of the Marchenko-Pastur law in random matrix theory derived under the assumption that the set of unsatisfied constraints are uncorrelated. This leads to a nontrivial critical exponent of the relaxation time in the SAT phase. Using a simple scaling analysis, we show that the isolation of this first eigenvalue from the bulk of spectrum is attributed to the force balance at the SAT-UNSAT transition point. Finally, we show that the estimated critical exponent of the relaxation time in the nonconvex region agrees very well with that of frictionless spherical particles, which have been studied in the context of the jamming transition of granular materials.
Collapse
Affiliation(s)
- Sungmin Hwang
- LPTMS, Université Paris-Sud 11, UMR 8626 CNRS, Bât. 100, 91405 Orsay Cedex, France
| | - Harukuni Ikeda
- Laboratoire de Physique de l'École Normale Supérieure, Université PSL, CNRS, Sorbonne Université, Université de Paris, 75005 Paris, France
| |
Collapse
|
17
|
Saitoh K, Hatano T, Ikeda A, Tighe BP. Stress Relaxation above and below the Jamming Transition. PHYSICAL REVIEW LETTERS 2020; 124:118001. [PMID: 32242697 DOI: 10.1103/physrevlett.124.118001] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/22/2019] [Accepted: 02/26/2020] [Indexed: 06/11/2023]
Abstract
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, which is a key ingredient of linear viscoelastic theory. The vibrational density of states exhibits an extra slow mode that emerges below jamming, which we utilize to demonstrate the anomalous exponent below jamming.
Collapse
Affiliation(s)
- Kuniyasu Saitoh
- Research Alliance Center for Mathematical Sciences, Tohoku University, 2-1-1 Katahira, Aoba-ku, Sendai 980-8577, Japan
| | - Takahiro Hatano
- Department of Earth and Space Science, Osaka University, 560-0043 Osaka, Japan
| | - Atsushi Ikeda
- Graduate School of Arts and Sciences, University of Tokyo, Tokyo 3-8-1, Japan
- Research Center for Complex Systems Biology, Universal Biology Institute, University of Tokyo, Komaba, Tokyo 153-8902, Japan
| | - Brian P Tighe
- Delft University of Technology, Process & Energy Laboratory, Leeghwaterstraat 39, 2628 CB Delft, Netherlands
| |
Collapse
|