Beyond linear elasticity: Jammed solids at finite shear strain and rate

The shear response of soft solids can be modeled with linear elasticity, provided the forcing is slow and weak. Both of these approximations must break down when the material loses rigidity, such as in foams and emulsions at their (un)jamming point -- suggesting that the window of linear elastic response near jamming is exceedingly narrow. Yet precisely when and how this breakdown occurs remains unclear. To answer these questions, we perform computer simulations of stress relaxation and shear startup experiments in athermal soft sphere packings, the canonical model for jamming. By systematically varying the strain amplitude, strain rate, distance to jamming, and system size, we identify characteristic strain and time scales that quantify how and when the window of linear elasticity closes, and relate these scales to changes in the microscopic contact network. Our findings indicate that the mechanical response of jammed solids are generically nonlinear and rate-dependent on experimentally accessible strain and time scales.Comment: 10 pages, 9 figure

Softening and Yielding of Soft Glassy Materials

Solids deform and fluids flow, but soft glassy materials, such as emulsions, foams, suspensions, and pastes, exhibit an intricate mix of solid and liquid-like behavior. While much progress has been made to understand their elastic (small strain) and flow (infinite strain) properties, such understanding is lacking for the softening and yielding phenomena that connect these asymptotic regimes. Here we present a comprehensive framework for softening and yielding of soft glassy materials, based on extensive numerical simulations of oscillatory rheological tests, and show that two distinct scenarios unfold depending on the material's packing density. For dense systems, there is a single, pressure-independent strain where the elastic modulus drops and the particle motion becomes diffusive. In contrast, for weakly jammed systems, a two-step process arises: at an intermediate softening strain, the elastic and loss moduli both drop down and then reach a new plateau value, whereas the particle motion becomes diffusive at the distinctly larger yield strain. We show that softening is associated with an extensive number of microscopic contact changes leading to a non-analytic rheological signature. Moreover, the scaling of the softening strain with pressure suggest the existence of a novel pressure scale above which softening and yielding coincide, and we verify the existence of this crossover scale numerically. Our findings thus evidence the existence of two distinct classes of soft glassy materials -- jamming dominated and dense -- and show how these can be distinguished by their rheological fingerprint.Comment: 9 pages, 11 figures, to appear in Soft Matte

Aggregation in a mixture of Brownian and ballistic wandering particles

In this paper, we analyze the scaling properties of a model that has as limiting cases the diffusion-limited aggregation (DLA) and the ballistic aggregation (BA) models. This model allows us to control the radial and angular scaling of the patterns, as well as, their gap distributions. The particles added to the cluster can follow either ballistic trajectories, with probability $P_{ba}$, or random ones, with probability $P_{rw}=1-P_{ba}$. The patterns were characterized through several quantities, including those related to the radial and angular scaling. The fractal dimension as a function of $P_{ba}$ continuously increases from $d_f\approx 1.72$ (DLA dimensionality) for $P_{ba}=0$ to $d_f\approx 2$ (BA dimensionality) for $P_{ba}=1$. However, the lacunarity and the active zone width exhibt a distinct behavior: they are convex functions of $P_{ba}$ with a maximum at $P_{ba}\approx1/2$. Through the analysis of the angular correlation function, we found that the difference between the radial and angular exponents decreases continuously with increasing $P_{ba}$ and rapidly vanishes for $P_{ba}>1/2$, in agreement with recent results concerning the asymptotic scaling of DLA clusters.Comment: 7 pages, 6 figures. accepted for publication on PR

Conserved Growth on Vicinal Surfaces

A crystal surface which is miscut with respect to a high symmetry plane exhibits steps with a characteristic distance. It is argued that the continuum description of growth on such a surface, when desorption can be neglected, is given by the anisotropic version of the conserved KPZ equation (T. Sun, H. Guo, and M. Grant, Phys. Rev. A 40, 6763 (1989)) with non-conserved noise. A one--loop dynamical renormalization group calculation yields the values of the dynamical exponent and the roughness exponent which are shown to be the same as in the isotropic case. The results presented here should apply in particular to growth under conditions which are typical for molecular beam epitaxy.Comment: 10 pages, uses revte

A Hybrid Monte Carlo Method for Surface Growth Simulations

We introduce an algorithm for treating growth on surfaces which combines important features of continuum methods (such as the level-set method) and Kinetic Monte Carlo (KMC) simulations. We treat the motion of adatoms in continuum theory, but attach them to islands one atom at a time. The technique is borrowed from the Dielectric Breakdown Model. Our method allows us to give a realistic account of fluctuations in island shape, which is lacking in deterministic continuum treatments and which is an important physical effect. Our method should be most important for problems close to equilibrium where KMC becomes impractically slow.Comment: 4 pages, 5 figure

Numerical Study of the Stress Response of Two-Dimensional Dense Granular Packings

We investigate the Green function of two-dimensional dense random packings of grains in order to discriminate between the different theories of stress transmission in granular materials. Our computer simulations allow for a detailed quantitative investigation of the dynamics which is difficult to obtain experimentally. We show that both hyperbolic and parabolic models of stress transmission fail to predict the correct stress distribution in the studied region of the parameters space. We demonstrate that the compressional and shear components of the stress compare very well with the predictions of isotropic elasticity for a wide range of pressures and porosities and for both frictional and frictionless packings. However, the states used in this study do not include the critical isostatic point for frictional particles, so that our results do not preclude the fact that corrections to elasticity may appear at the critical point of jamming, or for other sample preparation protocols, as discussed in the main text. We show that the agreement holds in the bulk of the packings as well as at the boundaries and we validate the linear dependence of the stress profile width with depth.Comment: 7 pages, 5 figure

Laplacian growth with separately controlled noise and anisotropy

Conformal mapping models are used to study competition of noise and anisotropy in Laplacian growth. For that, a new family of models is introduced with the noise level and directional anisotropy controlled independently. Fractalization is observed in both anisotropic growth and the growth with varying noise. Fractal dimension is determined from cluster size scaling with its area. For isotropic growth we find d = 1.7, both at high and low noise. For anisotropic growth with reduced noise the dimension can be as low as d = 1.5 and apparently is not universal. Also, we study fluctuations of particle areas and observe, in agreement with previous studies, that exceptionally large particles may appear during the growth, leading to pathologically irregular clusters. This difficulty is circumvented by using an acceptance window for particle areas.Comment: 13 pages, 15 figure

Effect of boundaries on the force distributions in granular media

The effect of boundaries on the force distributions in granular media is illustrated by simulations of 2D packings of frictionless, Hertzian spheres. To elucidate discrepancies between experimental observations and theoretical predictions, we distinguish between the weight distribution {\cal P} (w) measured in experiments and analyzed in the q-model, and the distribution of interparticle forces P(f). The latter one is robust, while {\cal P}(w) can be obtained once the local packing geometry and P(f) are known. By manipulating the (boundary) geometry, we show that {\cal P}(w) can be varied drastically.Comment: 4 pages, 4 figure

Conformal mapping methods for interfacial dynamics

The article provides a pedagogical review aimed at graduate students in materials science, physics, and applied mathematics, focusing on recent developments in the subject. Following a brief summary of concepts from complex analysis, the article begins with an overview of continuous conformal-map dynamics. This includes problems of interfacial motion driven by harmonic fields (such as viscous fingering and void electromigration), bi-harmonic fields (such as viscous sintering and elastic pore evolution), and non-harmonic, conformally invariant fields (such as growth by advection-diffusion and electro-deposition). The second part of the article is devoted to iterated conformal maps for analogous problems in stochastic interfacial dynamics (such as diffusion-limited aggregation, dielectric breakdown, brittle fracture, and advection-diffusion-limited aggregation). The third part notes that all of these models can be extended to curved surfaces by an auxilliary conformal mapping from the complex plane, such as stereographic projection to a sphere. The article concludes with an outlook for further research.Comment: 37 pages, 12 (mostly color) figure

Diffusion limited aggregation as a Markovian process: site-sticking conditions

Cylindrical lattice diffusion limited aggregation (DLA), with a narrow width N, is solved for site-sticking conditions using a Markovian matrix method (which was previously developed for the bond-sticking case). This matrix contains the probabilities that the front moves from one configuration to another at each growth step, calculated exactly by solving the Laplace equation and using the proper normalization. The method is applied for a series of approximations, which include only a finite number of rows near the front. The fractal dimensionality of the aggregate is extrapolated to a value near 1.68.Comment: 27 Revtex pages, 16 figure