Comparison of compression versus shearing near jamming, for a simple model of athermal frictionless disks in suspension.

Using a simplified model for a non-Brownian suspension, we numerically study the response of athermal, overdamped, frictionless disks in two dimensions to isotropic and uniaxial compression, as well as to pure and simple shearing, all at finite constant strain rates ε[over ̇]. We show that isotropic and uniaxial compression result in the same jamming packing fraction ϕ_{J}, while pure-shear- and simple-shear-induced jamming occurs at a slightly higher ϕ_{J}^{*}, consistent with that found previously for simple shearing. A critical scaling analysis of pure shearing gives critical exponents consistent with those previously found for both isotropic compression and simple shearing.

Universality of stress-anisotropic and stress-isotropic jamming of frictionless spheres in three dimensions: Uniaxial versus isotropic compression.

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.

Critical scaling of compression-driven jamming of athermal frictionless spheres in suspension.

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.

Dynamic length scales in athermal, shear-driven jamming of frictionless disks in two dimensions.

We carry out numerical simulations of athermally sheared, bidisperse, frictionless disks in two dimensions. From an appropriately defined velocity correlation function, we determine that there are two diverging length scales, ξ and ℓ, as the jamming transition is approached. We analyze our results using a critical scaling ansatz for the correlation function and argue that the more divergent length ℓ is a consequence of a dangerous irrelevant scaling variable and that it is ξ, which is the correlation length that determines the divergence of the system viscosity as jamming is approached from below in the liquid phase.

Depletion forces in athermally sheared mixtures of frictionless disks and rods in two dimensions.

We carry out numerical simulations to study the behavior of an athermal mixture of frictionless circular disks and elongated rods in two dimensions, under three different types of global linear deformation at a finite strain rate: (i) simple shearing, (ii) pure shearing, and (iii) isotropic compression. We find that the fluctuations induced by such deformations lead to depletion forces that cause rods to group in parallel oriented clusters for the cases of simple and pure shear, but not for isotropic compression. For simple shearing, we find that as the fraction of rods increases, this clustering increases, leading to an increase in the average rate of rotation of the rods, and a decrease in the magnitude of their nematic ordering.

Shear-driven flow of athermal, frictionless, spherocylinder suspensions in two dimensions: Particle rotations and orientational ordering.

We use numerical simulations to study the flow of a bidisperse mixture of athermal, frictionless, soft-core two-dimensional spherocylinders driven by a uniform steady-state simple shear applied at a fixed volume and a fixed finite strain rate γ[over ̇]. Energy dissipation is via a viscous drag with respect to a uniformly sheared host fluid, giving a simple model for flow in a non-Brownian suspension with Newtonian rheology. Considering a range of packing fractions ϕ and particle asphericities α at small γ[over ̇], we study the angular rotation θ[over ̇]_{i} and the nematic orientational ordering S_{2} of the particles induced by the shear flow, finding a nonmonotonic behavior as the packing ϕ is varied.

Shear-driven flow of athermal, frictionless, spherocylinder suspensions in two dimensions: Spatial structure and correlations.

We use numerical simulations to study the flow of athermal, frictionless, soft-core two-dimensional spherocylinders driven by a uniform steady-state simple shear applied at a fixed volume and a fixed finite strain rate γ[over ̇]. Energy dissipation is via a viscous drag with respect to a uniformly sheared host fluid, giving a simple model for flow in a non-Brownian suspension with Newtonian rheology. We study the resulting spatial structure of the sheared system, and compute correlation functions of the velocity, the particle density, the nematic order parameter, and the particle angular velocity.

Shear-driven flow of athermal, frictionless, spherocylinder suspensions in two dimensions: Stress, jamming, and contacts.

We use numerical simulations to study the flow of a bidisperse mixture of athermal, frictionless, soft-core two-dimensional spherocylinders driven by a uniform steady-state shear strain applied at a fixed finite rate. Energy dissipation occurs via a viscous drag with respect to a uniformly sheared host fluid, giving a simple model for flow in a non-Brownian suspension and resulting in a Newtonian rheology. We study the resulting pressure p and deviatoric shear stress σ of the interacting spherocylinders as a function of packing fraction ϕ, strain rate γ[over ̇], and a parameter α that measures the asphericity of the particles; α is varied to consider the range from nearly circular disks to elongated rods.

Orientational Ordering in Athermally Sheared, Aspherical, Frictionless Particles.

We numerically simulate the uniform athermal shearing of bidisperse, frictionless, two-dimensional spherocylinders and three-dimensional prolate ellipsoids. We focus on the orientational ordering of particles as an asphericity parameter α→0 and particles approach spherical. We find that the nematic order parameter S_{2} is nonmonotonic in the packing fraction ϕ and that, as α→0, S_{2} stays finite at jamming and above.

Compression-driven jamming of athermal frictionless spherocylinders in two dimensions.

We simulate numerically the compression-driven jamming of athermal, frictionless, soft-core spherocylinders in two dimensions, for a range of particle aspect ratios α. We find the critical packing fraction ϕ_{J}(α) for the jamming transition and the average number of contacts per particle z_{J}(α) at jamming. We find that both are nonmonotonic, with a peak at α≈1.

Anomalous stress fluctuations in athermal two-dimensional amorphous solids.

We numerically study the local stress distribution within athermal, isotropically stressed, mechanically stable, packings of bidisperse frictionless disks above the jamming transition in two dimensions. Considering the Fourier transform of the local stress, we find evidence for algebraically increasing fluctuations in both isotropic and anisotropic components of the stress tensor at small wave numbers, contrary to recent theoretical predictions. Such increasing fluctuations imply a lack of self-averaging of the stress on large length scales.

Shear banding, discontinuous shear thickening, and rheological phase transitions in athermally sheared frictionless disks.

We report on numerical simulations of simple models of athermal, bidisperse, soft-core, massive disks in two dimensions, as a function of packing fraction ϕ, inelasticity of collisions as measured by a parameter Q, and applied uniform shear strain rate γ[over ̇]. Our particles have contact interactions consisting of normally directed elastic repulsion and viscous dissipation, as well as tangentially directed viscous dissipation, but no interparticle Coulombic friction. Mapping the phase diagram in the (ϕ,Q) plane for small γ[over ̇], we find a sharp first-order rheological phase transition from a region with Bagnoldian rheology to a region with Newtonian rheology, and show that the system is always Newtonian at jamming.

Effect of collisional elasticity on the Bagnold rheology of sheared frictionless two-dimensional disks.

We carry out constant volume simulations of steady-state, shear-driven flow in a simple model of athermal, bidisperse, soft-core, frictionless disks in two dimensions, using a dissipation law that gives rise to Bagnoldian rheology. Focusing on the small strain rate limit, we map out the rheological behavior as a function of particle packing fraction ϕ and a parameter Q that measures the elasticity of binary particle collisions. We find a Q^{*}(ϕ) that marks the clear crossover from a region characteristic of strongly inelastic collisions, QQ^{*}, and give evidence that Q^{*}(ϕ) diverges as ϕ→ϕ_{J}, the shear-driven jamming transition.

Critical scaling of Bagnold rheology at the jamming transition of frictionless two-dimensional disks.

We carry out constant volume simulations of steady-state shear-driven rheology in a simple model of bidisperse soft-core frictionless disks in two dimensions, using a dissipation law that gives rise to Bagnoldian rheology. We discuss in detail the critical scaling ansatz for the shear-driven jamming transition and carry out a detailed scaling analysis of our resulting data for pressure p and shear stress σ. Our analysis determines the critical exponent β that describes the algebraic divergence of the Bagnold transport coefficients lim_{γ[over ̇]→0}p/γ[over ̇]^{2},σ/γ[over ̇]^{2}∼(ϕ_{J}-ϕ)^{-β} as the jamming transition ϕ_{J} is approached from below.

Search for hyperuniformity in mechanically stable packings of frictionless disks above jamming.

We numerically simulate mechanically stable packings of soft-core, frictionless, bidisperse disks in two dimensions, above the jamming packing fraction ϕ(J). For configurations with a fixed isotropic global stress tensor, we investigate the fluctuations of the local packing fraction ϕ(r) to test whether such configurations display the hyperuniformity that has been claimed to exist exactly at ϕ(J). For our configurations, generated by a rapid quench protocol, we find that hyperuniformity persists only out to a finite length scale and that this length scale appears to remain finite as the system stress decreases towards zero, i.

Maximum entropy and the stress distribution in soft disk packings above jamming.

We show that the maximum entropy hypothesis can successfully explain the distribution of stresses on compact clusters of particles within disordered mechanically stable packings of soft, isotropically stressed, frictionless disks above the jamming transition. We show that, in our two-dimensional case, it becomes necessary to consider not only the stress but also the Maxwell-Cremona force-tile area as a constraining variable that determines the stress distribution. The importance of the force-tile area had been suggested by earlier computations on an idealized force-network ensemble.

Statistics of conserved quantities in mechanically stable packings of frictionless disks above jamming.

We numerically simulate mechanically stable packings of soft-core, frictionless, bidisperse disks in two dimensions, above the jamming packing fraction ϕ(J). For configurations with a fixed isotropic global stress tensor, we compute the averages, variances, and correlations of conserved quantities (stress Γ(C), force-tile area A(C), Voronoi volume V(C), number of particles N(C), and number of small particles N(sC)) on compact subclusters of particles C, as a function of the cluster size and the global system stress. We find several significant differences depending on whether the cluster C is defined by a fixed radius R or a fixed number of particles M.

Pressure distribution and critical exponent in statically jammed and shear-driven frictionless disks.

We numerically study the distributions of global pressure that are found in ensembles of statically jammed and quasistatically sheared systems of bidisperse, frictionless disks at fixed packing fraction ϕ in two dimensions. We use these distributions to address the question of how pressure increases as ϕ increases above the jamming point ϕ(J), p ∼ |ϕ-ϕ(J)(y). For statically jammed ensembles, our results are consistent with the exponent y being simply related to the power law of the interparticle soft-core interaction.

Universality of jamming criticality in overdamped shear-driven frictionless disks.

We investigate the criticality of the jamming transition for overdamped shear-driven frictionless disks in two dimensions for two different models of energy dissipation: (i) Durian's bubble model with dissipation proportional to the velocity difference of particles in contact, and (ii) Durian's "mean-field" approximation to (i), with dissipation due to the velocity difference between the particle and the average uniform shear flow velocity. By considering the finite-size behavior of pressure, the pressure analog of viscosity, and the macroscopic friction σ/p, we argue that these two models share the same critical behavior.

Athermal jamming versus thermalized glassiness in sheared frictionless particles.

Numerical simulations of soft-core frictionless disks in two dimensions are carried out to study the behavior of a simple liquid as a function of temperature T, packing fraction φ, and uniform applied shear strain rate γ[over ·]. Inferring the hard-core limit from our soft-core results, we find that it depends on the two parameters φ and T/γ[over ·]. Here T/γ[over ·]→0 defines the athermal limit in which a shear-driven jamming transition occurs at a well defined φ(J) and T/γ[over ·]→∞ defines the thermalized limit where an equilibrium glass transition may take place at φ(G).

Herschel-Bulkley shearing rheology near the athermal jamming transition.

We consider the rheology of soft-core frictionless disks in two dimensions in the neighborhood of the athermal jamming transition. From numerical simulations of bidisperse, overdamped particles, we argue that the divergence of the viscosity below jamming is characteristic of the hard-core limit, independent of the particular soft-core interaction. We develop a mapping from soft-core to hard-core particles that recovers all the critical behavior found in earlier scaling analyses.

Combined topical and intracameral anesthesia for Descemet's stripping automated endothelial keratoplasty.

To evaluate the use of combined topical and intracameral anesthesia for Descemet's stripping automated keratoplasty (DSAEK). This was a retrospective comparative cohort analysis consisting of 10 eyes in 10 consecutive patients undergoing DSAEK surgery with combined topical and intracameral anesthesia. These cases were compared with 21 randomly selected controls during the same time period undergoing DSAEK surgery performed under retrobulbar anesthesia.

Glassiness, rigidity, and jamming of frictionless soft core disks.

The jamming of bidisperse soft core disks is considered, using a variety of different protocols to produce the jammed state. In agreement with other works, we find that cooling and compression can lead to a broad range of jamming packing fractions ϕ{J}, depending on cooling rate and initial configuration; the larger the degree of big particle clustering in the initial configuration, the larger will be the value of ϕ{J}. In contrast, we find that shearing disrupts particle clustering, leading to a much narrower range of ϕ{J} as the shear strain rate varies.

Finite-size scaling at the jamming transition: corrections to scaling and the correlation-length critical exponent.

We carry out a finite-size scaling analysis of the jamming transition in frictionless bidisperse soft core disks in two dimensions. We consider two different jamming protocols: (i) quench from random initial positions and (ii) quasistatic shearing. By considering the fraction of jammed states as a function of packing fraction for systems with different numbers of particles, we determine the spatial correlation length critical exponent ν ≈ 1 and show that corrections to scaling are crucial for analyzing the data.

Critical scaling of shearing rheology at the jamming transition of soft-core frictionless disks.

We perform numerical simulations to determine the shear stress and pressure of steady-state shear flow in a soft-disk model in two dimensions at zero temperature in the vicinity of the jamming transition ϕ{J}. We use critical point scaling analyses to determine the critical behavior at jamming, and we find that it is crucial to include corrections to scaling for a reliable analysis. We find that the relative size of these corrections are much smaller for pressure than for shear stress.