Experimental measures of affine and nonaffine deformation in granular shear

An explicit granular-mechanics approach to marine sediment acoustics

Here, we theoretically and computationally study the frequency dependence of phase speed and attenuation for marine sediments from the perspective of granular mechanics. We leverage recent theoretical insights from the granular physics community as well as discrete-element method simulations, where the granular material is treated as a packing of discrete objects that interact via pairwise forces. These pairwise forces include both repulsive contact forces as well as dissipative terms, which may include losses from the fluid as well as losses from inelasticity at grain-grain contacts. We show that the structure of disordered granular packings leads to anomalous scaling laws for frequency-dependent phase speed and attenuation that do not follow from a continuum treatment. Our results demonstrate that granular packing structure, which is not explicitly considered in existing models, may play a crucial role in a complete theory of sediment acoustics. While this simple approach does not explicitly treat sound propagation or inertial effects in the interstitial fluid, it provides a starting point for future models that include these and other more complex features.

Effects of polydispersity on the plastic behaviors of dense two-dimensional granular systems under shear

We study particle-scale motion in sheared highly polydisperse amorphous materials, in which the largest particles are as much as ten times the size of the smallest. We find strikingly different behavior from the more commonly studied amorphous systems with low polydispersity. In particular, an analysis of the nonaffine motion of particles reveals qualitative differences between large and small particles: The smaller particles have dramatically more nonaffine motion, which is induced by the presence of the large particles. We characterize how the nonaffine motion changes from the low- to high-polydispersity regimes. We further demonstrate a quantitative way to distinguish between "large" and "small" particles in systems with broad distributions of particle sizes. A macroscopic consequence of the nonaffine motion is a decrease in the energy dissipation rate for highly polydisperse samples, which is due both to a geometric consequence of the changing jamming conditions for higher polydispersity and to the changing character of nonaffine motion.

The role of non-affine deformations in the elastic behavior of the cellular vertex model

The vertex model of epithelia describes the apical surface of a tissue as a tiling of polygonal cells, with a mechanical energy governed by deviations in cell shape from preferred, or target, area, A₀, and perimeter, P₀. The model exhibits a rigidity transition driven by geometric incompatibility as tuned by the target shape index, . For with p_*(6) the perimeter of a regular hexagon of unit area, a cell can simultaneously attain both the preferred area and preferred perimeter. As a result, the tissue is in a mechanically soft compatible state, with zero shear and Young's moduli. For p₀ < p_*(6), it is geometrically impossible for any cell to realize the preferred area and perimeter simultaneously, and the tissue is in an incompatible rigid solid state. Using a mean-field approach, we present a complete analytical calculation of the linear elastic moduli of an ordered vertex model. We analyze a relaxation step that includes non-affine deformations, leading to a softer response than previously reported. The origin of the vanishing shear and Young's moduli in the compatible state is the presence of zero-energy deformations of cell shape. The bulk modulus exhibits a jump discontinuity at the transition and can be lower in the rigid state than in the fluid-like state. The Poisson's ratio can become negative which lowers the bulk and Young's moduli. Our work provides a unified treatment of linear elasticity for the vertex model and demonstrates that this linear response is protocol-dependent.

Tensor electromagnetism and emergent elasticity in jammed solids

The theory of mechanical response and stress transmission in disordered, jammed solids poses several open questions of how nonperiodic networks-apparently indistinguishable from a snapshot of a fluid-sustain shear. We present a stress-only theory of emergent elasticity for a nonthermal amorphous assembly of grains in a jammed solid, where each grain is subjected to mechanical constraints of force and torque balance. These grain-level constraints lead to the Gauss's law of an emergent U(1) tensor electromagnetism, which then accounts for the mechanical response of such solids. This formulation of amorphous elasticity has several immediate consequences. The mechanical response maps exactly to the static, dielectric response of this tensorial electromagnetism with the polarizability of the medium mapping to emergent elastic moduli. External forces act as vector electric charges, whereas the tensorial magnetic fields are sourced by momentum density. The dynamics in the electric and magnetic sectors naturally translate into the dynamics of the rigid jammed network and ballistic particle motion, respectively. The theoretical predictions for both stress-stress correlations and responses are borne out by the results of numerical simulations of frictionless granular packings in the static limit of the theory in both 2D and 3D.

Effect of nonaffine displacement on the mechanical performance of degraded PCL and its graphene composites: an atomistic investigation

Evaluating the mechanical properties of biodegradable implants can be challenging for in situ experiments and time-consuming for materials with a slow degradation rate, such as polycaprolactone (PCL). In this work, the effects of chain scission and water erosion on the mechanical properties of degraded PCL are investigated by molecular dynamics simulation. The decrease of the mechanical performance is correlated with the increase of the nonaffine displacement during the degradation. The nonaffine squared displacements (NSD) during the tensile deformation are calculated by subtracting the affine squared displacements from the mean squared displacements. After chain scission, short polymer chains increase the NSD of the system and weaken the modulus of the polymer matrix. The effect of the NSD is also observed in a water erosion model. When the bond break ratio is less than 5%, PCL still maintains a well-entangled network, which constrains the diffusion of the water molecules, resulting in a higher modulus of the erosion model than the chain scission model at a low degradation rate. The effect of NSD is also found in the PCL/graphene composites. For the degraded polymer chains, the diffusion of PCL is constrained by the graphene network, and such an effect increases during the degradation. As a result, the addition of graphene nanosheets slows down the decreasing trend of Young's modulus. Such findings can also explain the size effect of the graphene reinforcement on the mechanical properties of the polymer composites. This work provides atomistic insights into the mechanical property evolution during polymer degradation, revealing the possibility of tuning the mechanical performance by controlling the diffusion, which could be beneficial for the design and lifetime prediction of degradable implants.

Label-free cell tracking enables collective motion phenotyping in epithelial monolayers

Collective cell migration is an umbrella term for a rich variety of cell behaviors, whose distinct character is important for biological function, notably for cancer metastasis. One essential feature of collective behavior is the motion of cells relative to their immediate neighbors. We introduce an AI-based pipeline to segment and track cell nuclei from phase-contrast images. Nuclei segmentation is based on a U-Net convolutional neural network trained on images with nucleus staining. Tracking, based on the Crocker-Grier algorithm, quantifies nuclei movement and allows for robust downstream analysis of collective motion. Because the AI algorithm required no new training data, our approach promises to be applicable to and yield new insights for vast libraries of existing collective motion images. In a systematic analysis of a cell line panel with oncogenic mutations, we find that the collective rearrangement metric, D²_min, which reflects non-affine motion, shows promise as an indicator of metastatic potential.

•

Versatile AI algorithm identifies individual cell tracks in phase contrast images

•

Motion of cells relative to nearby neighbors may indicate cancer progression

Characterising shear-induced dynamics in flowing complex fluids using differential dynamic microscopy

Microscopic dynamics reveal the origin of the bulk rheological response in complex fluids. In model systems particle motion can be tracked, but for industrially relevant samples this is often impossible. Here we adapt differential dynamic microscopy (DDM) to study flowing highly-concentrated samples without particle resolution. By combining an investigation of oscillatory flow, using a novel "echo-DDM" analysis, and steady shear, through flow-DDM, we characterise the yielding of a silicone oil emulsion on both the microscopic and bulk level. Through measuring the rate of shear-induced droplet rearrangements and the flow velocity, the transition from a solid-like to liquid-like state is shown to occur in two steps: with droplet mobilisation marking the limit of linear visco-elasticity, followed by the development of shear localisation and macroscopic yielding. Using this suite of techniques, such insight could be developed for a wide variety of challenging complex fluids.

Friction and Adhesion Govern Yielding of Disordered Nanoparticle Packings: A Multiscale Adhesive Discrete Element Method Study

Recent studies have demonstrated that amorphous materials, from granular packings to atomic glasses, share multiple striking similarities, including a universal onset strain level for yield. This is despite vast differences in length scales and in the constituent particles' interactions. However, the nature of localized particle rearrangements is not well understood, and how local interactions affect overall performance remains unknown. Here, we introduce a multiscale adhesive discrete element method to simulate recent novel experiments of disordered nanoparticle packings indented and imaged with single nanoparticle resolution. The simulations exhibit multiple behaviors matching the experiments. By directly monitoring spatial rearrangements and interparticle bonding/debonding under the packing's surface, we uncover the mechanisms of the yielding and hardening phenomena observed in experiments. Interparticle friction and adhesion synergistically toughen the packings and retard plastic deformation. Moreover, plasticity can result from bond switching without particle rearrangements. These results furnish insights for understanding yielding in amorphous materials generally.

Using delaunay triangularization to characterize non-affine displacement fields during athermal, quasistatic deformation of amorphous solids

We investigate the non-affine displacement fields that occur in two-dimensional Lennard-Jones models of metallic glasses subjected to athermal, quasistatic simple shear (AQS). During AQS, the shear stress versus strain displays continuous quasi-elastic segments punctuated by rapid drops in shear stress, which correspond to atomic rearrangement events. We capture all information concerning the atomic motion during the quasi-elastic segments and shear stress drops by performing Delaunay triangularizations and tracking the deformation gradient tensor F_α associated with each triangle α. To understand the spatio-temporal evolution of the displacement fields during shear stress drops, we calculate F_α along minimal energy paths from the mechanically stable configuration immediately before to that after the stress drop. We find that quadrupolar displacement fields form and dissipate both during the quasi-elastic segments and shear stress drops. We then perform local perturbations (rotation, dilation, simple and pure shear) to single triangles and measure the resulting displacement fields. We find that local pure shear deformations of single triangles give rise to mostly quadrupolar displacement fields, and thus pure shear strain is the primary type of local strain that is activated by bulk, athermal quasistatic simple shear. Other local perturbations, e.g. rotations, dilations, and simple shear of single triangles, give rise to vortex-like and dipolar displacement fields that are not frequently activated by bulk AQS. These results provide fundamental insights into the non-affine atomic motion that occurs in driven, glassy materials.

Depletion attraction impairs the plasticity of emulsions flowing in a constriction

We study the elasto-plastic behavior of dense attractive emulsions under a mechanical perturbation. The attraction is introduced through non-specific depletion interactions between the droplets and is controlled by changing the concentration of surfactant micelles in the continuous phase. We find that such attractive forces are not sufficient to induce any measurable modification on the scalings between the local packing fraction and the deformation of the droplets. However, when the emulsions are flowed through 2D microfluidic constrictions, we uncover a measurable effect of attraction on their elasto-plastic response. Indeed, we measure higher levels of deformation inside the constriction for attractive droplets. In addition, we show that these measurements correlate with droplet rearrangements that are spatially delayed in the constriction for higher attraction forces.

Dynamics phenotyping across length and time scales in collective cell migration

Processes in collective migration span many length and time scales. In this review, we focus on length scales ranging from tens of microns (single cells) to a few millimeters (cell clusters) and the motion of these cells and cell groups on time scales of minutes to hours. We focus on epithelial cell sheets and metrics of motion developed to measure migration phenotypes in this system. Comparisons between cell motion and fluid flows, facilitated by the popular image analysis technique particle image velocimetry, yield metrics that can be used to study migration across a range of length and time scales. Measuring collective cell migration across these scales provides a complex, quantitative phenotype useful for migration models, in particular those that compare and contrast collective cell migration to movement of particles near a transition to jamming. Contrasting the motion of epithelial cells and the jamming transition illustrates aspects of collective motion that can be attributed to the jammed character of cell clusters, and highlights aspects of collective behavior that likely involve active motility and cell-cell guidance. The application of multiple migration metrics, which span multiple scales of the system, thus allows us to link cell-scale signals and mechanics to collective behavior.

Stable small bubble clusters in two-dimensional foams

Key features of the mechanical response of amorphous particulate materials, such as foams, emulsions, and granular media, to applied stress are determined by the frequency and size of particle rearrangements that occur as the system transitions from one mechanically stable state to another. This work describes coordinated experimental and computational studies of bubble rafts, which are quasi-two dimensional systems of bubbles confined to the air-water interface. We focus on small mechanically stable clusters of four, five, six, and seven bubbles with two different sizes with diameter ratio σ_L/σ_S ≃ 1.4. Focusing on small bubble clusters, which can be viewed as subsystems of a larger system, allows us to investigate the full ensemble of clusters that form, measure the respective frequencies with which the clusters occur, and determine the form of the bubble-bubble interactions. We emphasize several important results. First, for clusters with N > 5 bubbles, we find using discrete element simulations that short-range attractive interactions between bubbles give rise to a larger ensemble of distinct mechanically stable clusters compared to that generated by long-range attractive interactions. The additional clusters in systems with short-range attractions possess larger gaps between pairs of neighboring bubbles on the periphery of the clusters. The ensemble of bubble clusters observed in experiments is similar to the ensemble of clusters with long-range attractive interactions. We also compare the frequency with which each cluster occurs in simulations and experiments. We find that the cluster frequencies are extremely sensitive to the protocol used to generate them and only weakly correlated to the energy of the clusters.

Microscopic Description of the Granular Fluidity Field in Nonlocal Flow Modeling

A recent granular rheology based on an implicit "granular fluidity" field has been shown to quantitatively predict many nonlocal phenomena. However, the physical nature of the field has not been identified. Here, the granular fluidity is found to be a kinematic variable given by the velocity fluctuation and packing fraction. This is verified with many discrete element simulations, which show that the operational fluidity definition, solutions of the fluidity model, and the proposed microscopic formula all agree. Kinetic theoretical and Eyring-like explanations shed insight into the obtained form.

Measurement of Stress Redistribution in Flowing Emulsions

We study how local rearrangements alter droplet stresses within flowing dense quasi-two-dimensional emulsions at area fractions ϕ≥0.88. Using microscopy, we measure droplet positions while simultaneously using their deformed shape to measure droplet stresses. We find that rearrangements alter nearby stresses in a quadrupolar pattern: stresses on neighboring droplets tend to either decrease or increase depending on location. The stress redistribution is more anisotropic with increasing ϕ. The spatial character of the stress redistribution influences where subsequent rearrangements occur. Our results provide direct quantitative support for rheological theories of dense amorphous materials that connect local rearrangements to changes in nearby stress.

Experimental observation of local rearrangements in dense quasi-two-dimensional emulsion flow

We experimentally study rearranging regions in slow athermal flow by observing the flow of a concentrated oil-in-water emulsion in a thin chamber with a constricting hopper shape. The gap of the chamber is smaller than the droplet diameters, so that the droplets are compressed into quasi-two-dimensional pancakes. We focus on localized rearrangements known as "T1 events" where four droplets exchange neighbors. Flowing droplets are deformed due to forces from neighboring droplets, and these deformations are decreased by nearby T1 events, with a spatial dependence related to the local structure. We see a tendency of the T1 events to occur in small clusters.

Shear-induced rigidity in athermal materials: A unified statistical framework

Recent studies of athermal systems such as dry grains and dense, non-Brownian suspensions have shown that shear can lead to solidification through the process of shear jamming in grains and discontinuous shear thickening in suspensions. The similarities observed between these two distinct phenomena suggest that the physical processes leading to shear-induced rigidity in athermal materials are universal. We present a nonequilibrium statistical mechanics model, which exhibits the phenomenology of these shear-driven transitions, shear jamming and discontinuous shear thickening, in different regions of the predicted phase diagram. Our analysis identifies the crucial physical processes underlying shear-driven rigidity transitions, and clarifies the distinct roles played by shearing forces and the packing fraction of grains.

Affine and nonaffine motions in sheared polydisperse emulsions

We study dense and highly polydisperse emulsions at droplet volume fractions ϕ≥0.65. We apply oscillatory shear and observe droplet motion using confocal microscopy. The presence of droplets with sizes several times the mean size dramatically changes the motion of smaller droplets. Both affine and nonaffine droplet motions are observed, with the more nonaffine motion exhibited by the smaller droplets which are pushed around by the larger droplets. Droplet motions are correlated over length scales from one to four times the mean droplet diameter, with larger length scales corresponding to higher strain amplitudes (up to strains of about 6%).

Evolution of atomic rearrangements in deformation in metallic glasses

Atomic rearrangements induced by shear stress are fundamental for understanding deformation mechanisms in metallic glasses (MGs). Using molecular dynamic simulation, the atomic rearrangements characterized by nonaffine displacements (NADs) and their spatial distribution and evolution with tensile stress in Cu50Zr50 MG were investigated. It was found that in the elastic regime the atomic rearrangements with the largest NADs are relatively homogeneous in space, but exhibit strong spatial correlation, become localized and inhomogeneous, and form large clusters as strain increases, which may facilitate the so-called shear transformation zones. Furthermore, initially they prefer to take place around Cu atoms which have more nonicosahedral configurations. As strain increases, the preference decays and disappears in the plastic regime. The atomic rearrangements with the smallest NADs are preferentially located around Cu atoms, too, but with more icosahedral or icosahedral-like atomic configurations. The preference is maintained in the whole deformation process. In contrast, the atomic rearrangements with moderate NADs distribute homogeneously, and do not show explicit preference or spatial correlation, acting as matrix during deformation. Among the atomic rearrangements with different NADs, those with largest and smallest NADs are nearest neighbors initially, but separating with increasing strain, while those with largest and moderate NADs always avoid to each other. The correlations in the fluctuations of the NADs confirm the long-range strain correlation and the scale-free characteristic of NADs in both elastic and plastic deformation, which suggests a universality of the scaling in the plastic flow in MGs.

Effect of inertia on sheared disordered solids: critical scaling of avalanches in two and three dimensions

Molecular dynamics simulations with varying damping are used to examine the effects of inertia and spatial dimension on sheared disordered solids in the athermal quasistatic limit. In all cases the distribution of avalanche sizes follows a power law over at least three orders of magnitude in dissipated energy or stress drop. Scaling exponents are determined using finite-size scaling for systems with 10(3)-10(6) particles. Three distinct universality classes are identified corresponding to overdamped and underdamped limits, as well as a crossover damping that separates the two regimes. For each universality class, the exponent describing the avalanche distributions is the same in two and three dimensions. The spatial extent of plastic deformation is proportional to the energy dissipated in an avalanche. Both rise much more rapidly with system size in the underdamped limit where inertia is important. Inertia also lowers the mean energy of configurations sampled by the system and leads to an excess of large events like that seen in earthquake distributions for individual faults. The distribution of stress values during shear narrows to zero with increasing system size and may provide useful information about the size of elemental events in experimental systems. For overdamped and crossover systems the stress variation scales inversely with the square root of the system size. For underdamped systems the variation is determined by the size of the largest events.

Localization of folds and cracks in thin metal films coated on flexible elastomer foams

Thin metal films coated on soft elastomeric foam substrates exhibit enhanced electromechanical performance. The open-cell foam structure conveys highly anisotropic mechanical properties within the top, thin capping elastomer at the surface of the foam. Upon stretching, large strain fields inducing cracks and folds localize above the foam cells, while the surrounding cell ligaments remain almost strain-free, enabling stable electrical conduction in the metallic coating.

Quantifying stretching and rearrangement in epithelial sheet migration

Although understanding the collective migration of cells, such as that seen in epithelial sheets, is essential for understanding diseases such as metastatic cancer, this motion is not yet as well characterized as individual cell migration. Here we adapt quantitative metrics used to characterize the flow and deformation of soft matter to contrast different types of motion within a migrating sheet of cells. Using a Finite-Time Lyapunov Exponent (FTLE) analysis, we find that - in spite of large fluctuations - the flow field of an epithelial cell sheet is not chaotic. Stretching of a sheet of cells (i.e., positive FTLE) is localized at the leading edge of migration and increases when the cells are more highly stimulated. By decomposing the motion of the cells into affine and non-affine components using the metric D2min , we quantify local plastic rearrangements and describe the motion of a group of cells in a novel way. We find an increase in plastic rearrangements with increasing cell densities, whereas inanimate systems tend to exhibit less non-affine rearrangements with increasing density.

Deformation field in indentation of a granular ensemble

An experimental study has been made of the flow field in indentation of a model granular material. A granular ensemble composed of spherical sand particles with average size of 0.4 mm is indented with a flat ended punch under plane-strain conditions. The region around the indenter is imaged in situ using a high-speed charge-coupled device (CCD) imaging system. By applying a hybrid image analysis technique to image sequences of the indentation, flow parameters such as velocity, velocity gradient, and strain rate are measured at high resolution. The measurements have enabled characterization of the main features of the flow such as dead material zones, velocity jumps, localization of deformation, and regions of highly rotational flow resembling vortices. Implications for validation of theoretical analyses and applications are discussed.

Fracture pattern formation in frictional, cohesive, granular material

Naturally, deformed rocks commonly contain crack arrays (joints) forming patterns with systematic relationships to the large-scale deformation. Kinematically, joints can be mode-1, -2 or -3 or combinations of these, but there is no overarching theory for the development of the patterns. We develop a model motivated by dislocation pattern formation in metals. The problem is formulated in one dimension in terms of coupled reaction-diffusion equations, based on computer simulations of crack development in deformed granular media with cohesion. The cracks are treated as interacting defects, and the densities of defects diffuse through the rock mass. Of particular importance is the formation of cracks at high stresses associated with force-chain buckling and variants of this configuration; these cracks play the role of 'inhibitors' in reaction-diffusion relationships. Cracks forming at lower stresses act as relatively mobile defects. Patterns of localized deformation result from (i) competition between the growth of the density of 'mobile' defects and the inhibition of these defects by crack configurations forming at high stress and (ii) the diffusion of damage arising from these two populations each characterized by a different diffusion coefficient. The extension of this work to two and three dimensions is discussed.

Microscopic structural relaxation in a sheared supercooled colloidal liquid

The rheology of dense amorphous materials under large shear strain is not fully understood, partly due to the difficulty of directly viewing the microscopic details of such materials. We use a colloidal suspension to simulate amorphous materials and study the shear-induced structural relaxation with fast confocal microscopy. We quantify the plastic rearrangements of the particles in several ways. Each of these measures of plasticity reveals spatially heterogeneous dynamics, with localized regions where many particles are strongly rearranging by these measures. We examine the shapes of these regions and find them to be essentially isotropic, with no alignment in any particular direction. Furthermore, individual particles are equally likely to move in any direction other than the overall bias imposed by the strain.

Soft random solids and their heterogeneous elasticity

Spatial heterogeneity in the elastic properties of soft random solids is examined via vulcanization theory. The spatial heterogeneity in the structure of soft random solids is a result of the fluctuations locked-in at their synthesis, which also brings heterogeneity in their elastic properties. Vulcanization theory studies semimicroscopic models of random-solid-forming systems and applies replica field theory to deal with their quenched disorder and thermal fluctuations. The elastic deformations of soft random solids are argued to be described by the Goldstone sector of fluctuations contained in vulcanization theory, associated with a subtle form of spontaneous symmetry breaking that is associated with the liquid-to-random-solid transition. The resulting free energy of this Goldstone sector can be reinterpreted as arising from a phenomenological description of an elastic medium with quenched disorder. Through this comparison, we arrive at the statistics of the quenched disorder of the elasticity of soft random solids in terms of residual stress and Lamé-coefficient fields. In particular, there are large residual stresses in the equilibrium reference state, and the disorder correlators involving the residual stress are found to be long ranged and governed by a universal parameter that also gives the mean shear modulus.

Diffusive waves in a dilating scattering medium

We consider the propagation of a diffusive wave in a scattering medium submitted to a homogeneous expansion. The light multiply scattered by a glass spheres sample is measured. We analyze the variations of the scattered light when the material and the optical wavelength are dilated. We experimentally show that an isotropic expansion of the material is equivalent to a contraction of the wavelength. Moreover, the effect of an expansion of the material on the scattered wave may be canceled by a proportional increase of the wavelength, keeping the phase of the scattered wave unchanged. Applications to the characterization of deformation of disordered materials are outlined.