Correlation between particle motion and Voronoi-cell-shape fluctuations during the compaction of granular matter

Evident structural anisotropies arising from near-zero particle asphericity in granular spherocylinder packings

With magnetic resonance imaging experiments, we study packings of granular spherocylinders with merely 2% asphericity. Evident structural anisotropies across all length scales are identified. Most interestingly, the global nematic order decreases with increasing packing fraction, while the local contact anisotropy shows an opposing trend. We attribute this counterintuitive phenomenon to a competition between gravity-driven ordering aided by frictional contacts and a geometric frustration effect at the marginally jammed state. It is also surprising to notice that such slight particle asphericity can trigger non-negligible correlations between contact-level and mesoscale structures, manifested in drastically different nonaffine structural rearrangements upon compaction from that of granular spheres. These observations can help improve statistical mechanical models for the orientational order transformation of nonspherical granular particle packings, which involves complex interplays between particle shape, frictional contacts, and external force field.

Thermal shakedown in granular materials with irregular particle shapes

Granular materials with irregular particle shapes undergo a myriad of temperature variations in natural and engineered systems. However, the impacts of cyclic temperature variations on the mechanics of granular materials remain poorly understood. Specifically, little is known about the response of granular materials to cyclic temperature variations as a function of the following central variables: particle shape, applied stress level, relative density, and temperature amplitude. This paper presents advanced laboratory experiments to explore the impacts of cyclic temperature variations on the mechanics of granular materials, with a focus on sands. The results show that cyclic temperature variations applied to sands induce thermal shakedown: the accumulation of irreversible bulk deformations due to microstructural rearrangements caused by thermal expansions and contractions of the constituting particles. The deformation of sands caused by thermal shakedown strongly depends on particle shape, stress level, relative density, and temperature amplitude. This deformation is limited for individual thermal cycles but accumulates and becomes significant for multiple thermal cycles, leading to substantial compaction in sands and other granular materials, which can affect various natural and engineered systems.

Rearrangements during slow compression of a jammed two-dimensional emulsion

As amorphous materials get jammed, both geometric and dynamic heterogeneity are observed. We investigate the correlation between the local geometric heterogeneity and local rearrangements in a slowly compressed bidisperse quasi-two-dimensional emulsion system. The compression is driven by evaporation of the continuous phase and causes the area packing fraction to increase from 0.88 to 0.99. We quantify the structural heterogeneity of the system using the radical Voronoi tessellation following the method of Rieser et al. [Phys. Rev. Lett. 116, 088001 (2016)]0031-900710.1103/PhysRevLett.116.088001. We define two structural quantities characterizing local structure, the first of which considers nearest neighbors and the second of which includes information from second-nearest neighbors. We find that droplets in heterogeneous local regions are more likely to have local rearrangements. These rearrangements are generally T1 events where two droplets converge toward a void, and two droplets move away from the void to make room for the converging droplets. Thus, the presence of the voids tends to orient the T1 events. The presence of a correlation between the structural quantities and the rearrangement dynamics remains qualitatively unchanged over the entire range of packing fractions observed.

Soft matter dynamics: A versatile microgravity platform to study dynamics in soft matter

We describe an experiment container with light scattering and imaging diagnostics for experiments on soft matter aboard the International Space Station (ISS). The suite of measurement capabilities can be used to study different materials in exchangeable sample cell units. The currently available sample cell units and future possibilities for foams, granular media, and emulsions are presented in addition to an overview of the design and the diagnostics of the experiment container. First results from measurements performed on ground and during the commissioning aboard the ISS highlight the capabilities of the experiment container to study the different materials.

Phase behaviour and dynamical features of a two-dimensional binary mixture of active/passive spherical particles

In this work we have characterized the phase behaviour and the dynamics of bidimensional mixtures of active and passive Brownian particles. We have evaluated state diagrams at several concentrations of the passive components finding that, while passive agents tend to hinder phase separation, active agents force crystal-like structures on passive colloids. In order to study how passive particles affect the dynamics of the mixture, we have computed the long-time diffusion coefficient of each species, concluding that active particles induce activity and super-diffusive behaviour on passive ones. Interestingly, at the density at which the system enters a MIPS state the active particles' diffusivity shows an inflection point and the passive particles' one goes through a maximum, due to the change in the dynamics of the active components, as shown in the displacement's probability distribution function.

Machine learning characterization of structural defects in amorphous packings of dimers and ellipses

Structural defects within amorphous packings of symmetric particles can be characterized using a machine learning approach that incorporates structure functions of radial distances and angular arrangement. This yields a scalar field, softness, that correlates with the probability that a particle is about to be rearranged. However, when particle shapes are elongated, as in the case of dimers and ellipses, we find that the standard structure functions produce imprecise softness measurements. Moreover, ellipses exhibit deformation profiles in stark contrast to circular particles. In order to account for the effects of orientation and alignment, we introduce structure functions to recover the predictive performance of softness, as well as provide physical insight into local and extended dynamics. We study a model disordered solid, a bidisperse two-dimensional granular pillar, driven by uniaxial compression and composed entirely of monomers, dimers, or ellipses. We demonstrate how the computation of softness via a support vector machine extends to dimers and ellipses with the introduction of orientational structure functions. Then we highlight the spatial extent of rearrangements and defects, as well as their cross correlation, for each particle shape. Finally, we demonstrate how an additional machine learning algorithm, recursive feature elimination, provides an avenue to better understand how softness arises from particular structural aspects. We identify the most crucial structure functions in determining softness and discuss their physical implications.

Scale-free crystallization of two-dimensional complex plasmas: Domain analysis using Minkowski tensors

Experiments of the recrystallization processes in two-dimensional complex plasmas are analyzed to rigorously test a recently developed scale-free phase transition theory. The "fractal-domain-structure" (FDS) theory is based on the kinetic theory of Frenkel. It assumes the formation of homogeneous domains, separated by defect lines, during crystallization and a fractal relationship between domain area and boundary length. For the defect number fraction and system energy a scale-free power-law relation is predicted. The long-range scaling behavior of the bond-order correlation function shows clearly that the complex plasma phase transitions are not of the Kosterlitz, Thouless, Halperin, Nelson, and Young type. Previous preliminary results obtained by counting the number of dislocations and applying a bond-order metric for structural analysis are reproduced. These findings are supplemented by extending the use of the bond-order metric to measure the defect number fraction and furthermore applying state-of-the-art analysis methods, allowing a systematic testing of the FDS theory with unprecedented scrutiny: A morphological analysis of lattice structure is performed via Minkowski tensor methods. Minkowski tensors form a complete family of additive, motion covariant and continuous morphological measures that are sensitive to nonlinear properties. The FDS theory is rigorously confirmed and predictions of the theory are reproduced extremely well. The predicted scale-free power-law relation between defect fraction number and system energy is verified for one more order of magnitude at high energies compared to the inherently discontinuous bond-order metric. It is found that the fractal relation between crystalline domain area and circumference is independent of the experiment, the particular Minkowski tensor method, and the particular choice of parameters. Thus, the fractal relationship seems to be inherent to two-dimensional phase transitions in complex plasmas. Minkowski tensor analysis turns out to be a powerful tool for investigations of crystallization processes. It is capable of revealing nonlinear local topological properties, however, still provides easily interpretable results founded on a solid mathematical framework.

Note: Eliminating stripe artifacts in light-sheet fluorescence imaging

We report two techniques to mitigate stripe artifacts in light-sheet fluorescence imaging. The first uses an image processing algorithm called the multidirectional stripe remover method to filter stripes from an existing image. The second uses an elliptical holographic diffuser with strong scattering anisotropy to prevent stripe formation during image acquisition. These techniques facilitate accurate interpretation of image data, especially in denser samples. They are also facile and cost-effective.

Nucleation in Sheared Granular Matter

We present an experiment on crystallization of packings of macroscopic granular spheres. This system is often considered to be a model for thermally driven atomic or colloidal systems. Cyclically shearing a packing of frictional spheres, we observe a first order phase transition from a disordered to an ordered state. The ordered state consists of crystallites of mixed fcc and hcp symmetry that coexist with the amorphous bulk. The transition, initiated by homogeneous nucleation, overcomes a barrier at 64.5% volume fraction. Nucleation consists predominantly of the dissolving of small nuclei and the growth of nuclei that have reached a critical size of about ten spheres.

Can exotic disordered "stealthy" particle configurations tolerate arbitrarily large holes?

The probability of finding a spherical cavity or "hole" of arbitrarily large size in typical disordered many-particle systems in the infinite-system-size limit (e.g., equilibrium liquid states) is non-zero. Such "hole" statistics are intimately linked to the thermodynamic and nonequilibrium physical properties of the system. Disordered "stealthy" many-particle configurations in d-dimensional Euclidean space [Doublestruck R]^d are exotic amorphous states of matter that lie between a liquid and crystal that prohibit single-scattering events for a range of wave vectors and possess no Bragg peaks [Torquato et al., Phys. Rev. X, 2015, 5, 021020]. In this paper, we provide strong numerical evidence that disordered stealthy configurations across the first three space dimensions cannot tolerate arbitrarily large holes in the infinite-system-size limit, i.e., the hole probability has compact support. This structural "rigidity" property apparently endows disordered stealthy systems with novel thermodynamic and physical properties, including desirable band-gap, optical and transport characteristics. We also determine the maximum hole size that any stealthy system can possess across the first three space dimensions.

Kinetics of fluid demixing in complex plasmas: Domain growth analysis using Minkowski tensors

A molecular dynamics simulation of the demixing process of a binary complex plasma is analyzed and the role of distinct interaction potentials is discussed by using morphological Minkowski tensor analysis of the minority phase domain growth in a demixing simulated binary complex plasma. These Minkowski tensor methods are compared with previous results that utilized a power spectrum method based on the time-dependent average structure factor. It is shown that the Minkowski tensor methods are superior to the previously used power-spectrum method in the sense of higher sensitivity to changes in domain size. By analysis of the slope of the temporal evolution of Minkowski tensor measures, qualitative differences between the case of particle interaction with a single length scale compared to particle interactions with two different length scales (dominating long-range interaction) are revealed. After proper scaling the graphs for the two length scale scenarios coincide, pointing toward universal behavior. The qualitative difference in demixing scenarios is evidenced by distinct demixing behavior: in the long-range dominated cases demixing occurs in two stages. At first, neighboring particles agglomerate, then domains start to merge in cascades. However, in the case of only one interaction length scale only agglomeration but no merging of domains can be observed. Thus, Minkowski tensor analysis is likely to become a useful tool for further investigation of this (and other) demixing processes. It is capable to reveal (nonlinear) local topological properties, probing deeper than (linear) global power-spectrum analysis, however, still providing easily interpretable results founded on a solid mathematical framework.

Divergence of Voronoi Cell Anisotropy Vector: A Threshold-Free Characterization of Local Structure in Amorphous Materials

Characterizing structural inhomogeneity is an essential step in understanding the mechanical response of amorphous materials. We introduce a threshold-free measure based on the field of vectors pointing from the center of each particle to the centroid of the Voronoi cell in which the particle resides. These vectors tend to point in toward regions of high free volume and away from regions of low free volume, reminiscent of sinks and sources in a vector field. We compute the local divergence of these vectors, where positive values correspond to overpacked regions and negative values identify underpacked regions within the material. Distributions of this divergence are nearly Gaussian with zero mean, allowing for structural characterization using only the moments of the distribution. We explore how the standard deviation and skewness vary with the packing fraction for simulations of bidisperse systems and find a kink in these moments that coincides with the jamming transition.

Universality of slip avalanches in flowing granular matter

The search for scale-bridging relations in the deformation of amorphous materials presents a current challenge with tremendous applications in material science, engineering and geology. While generic features in the flow and microscopic dynamics support the idea of a universal scaling theory of deformation, direct microscopic evidence remains poor. Here, we provide the first measurement of internal scaling relations in the deformation of granular matter. By combining macroscopic force fluctuation measurements with internal strain imaging, we demonstrate the existence of robust scaling relations from particle-scale to macroscopic flow. We identify consistent power-law relations truncated by systematic pressure-dependent cutoff, in agreement with recent mean-field theory of slip avalanches in elasto-plastic materials, revealing the existence of a mechanical critical point. These results experimentally establish scale-bridging relations in the flow of matter, paving the way to a new universal theory of deformation.

Whether the deformation of amorphous materials is governed by universal scaling, a characteristic feature of the critical phenomena, is currently under debate. Here, Denisov et al. provide experimental evidence by linking the internal strains at microscales to the fluctuations in the applied force.

Geometric order parameters derived from the Voronoi tessellation show signatures of the jamming transition

A jammed packing of frictionless spheres at zero temperature is perfectly specified by the network of contact forces from which mechanical properties can be derived. However, we can alternatively consider a packing as a geometric structure, characterized by a Voronoi tessellation which encodes the local environment around each particle. We find that this local environment characterizes systems both above and below jamming and changes markedly at the transition. A variety of order parameters derived from this tessellation carry signatures of the jamming transition, complete with scaling exponents. Furthermore, we define a real space geometric correlation function which also displays a signature of jamming. Taken together, these results demonstrate the validity and usefulness of a purely geometric approach to jamming.

Structural characterization of submerged granular packings

We consider the impact of the effective gravitational acceleration on microstructural properties of granular packings through experimental studies of spherical granular materials saturated within fluids of varying density. We characterize the local organization of spheres in terms of contact connectivity, distribution of the Delaunay free volumes, and the shape factor (parameter of nonsphericity) of the Voronoï polygons. The shape factor gives a clear physical picture of the competition between less and more ordered domains of particles in experimentally obtained packings. As the effective gravity increases, the probability distribution of the shape factor becomes narrower and more localized around the lowest values of the shape factor corresponding to regular hexagon. It is found that curves of the pore distributions are asymmetric with a long tail on the right-hand side, which progressively reduces while the effective gravity gets stronger for lower densities of interstitial fluid. We show that the distribution of local areas (Voronoï cells) broadens with decreasing value of the effective gravity due to the formation of lose structures such as large pores and chainlike structures (arches or bridges). Our results should be particularly helpful in testing the newly developed simulation techniques involving liquid-related forces associated with immersed granular particles.

Effects of grain shape on packing and dilatancy of sheared granular materials

A granular material exposed to shear shows a variety of unique phenomena: Reynolds dilatancy, positional order and orientational order effects may compete in the shear zone. We study granular packing consisting of macroscopic prolate, oblate and spherical grains and compare their behaviour. X-ray tomography is used to determine the particle positions and orientations in a cylindrical split bottom shear cell. Packing densities and the arrangements of individual particles in the shear zone are evaluated. For anisometric particles, we observe the competition of two opposite effects. On the one hand, the sheared granules are dilated, on the other hand the particles reorient and align with respect to the streamlines. Even though aligned cylinders in principle may achieve higher packing densities, this alignment compensates for the effect of dilatancy only partially. The complex rearrangements lead to a depression of the surface above the well oriented region while neighbouring parts still show the effect of dilation in the form of heaps. For grains with isotropic shapes, the surface remains rather flat. Perfect monodisperse spheres crystallize in the shear zone, whereby positional order partially overcompensates dilatancy effects. However, even slight deviations from the ideal monodisperse sphere shape inhibit crystallization.

Granular dynamics during impact

We study the impact of a projectile onto a bed of 3 mm grains immersed in an index-matched fluid. We vary the amount of prestrain on the sample, strengthening the force chains within the system. We find this affects only the prefactor of the linear depth-dependent term in the stopping force. We propose a simple model to account for the strain dependence of this term, owing to increased pressure in the pile. Interestingly, we find that the presence of the fluid does not affect the impact dynamics, suggesting that dynamic friction is not a factor. Using a laser sheet scanning technique to visualize internal grain motion, we measure the trajectory of each grain throughout an impact. Microscopically, our results indicate that weaker initial force chains result in more irreversible, plastic rearrangements, suggesting static friction between grains does play a substantial role in the energy dissipation.

Characterization of void space in polydisperse sphere packings: Applications to hard-sphere packings and to protein structure analysis

The implementation of a method for the exact evaluation of the volume and surface area of cavities and free volumes in polydisperse sphere packings is described. The generalization of an algorithm for Voronoi tessellation by Tanemura et al. is presented, employing the radical plane construction, as a part of the method. We employ this method to calculate the equation of state for monodisperse and polydisperse hard-sphere fluids, crystals, and for the metastable amorphous branch up to random close packing or jamming densities. We compute the distribution of free volumes, and compare with previous results employing a heuristic definition of free volume. We show the efficacy of the method for analyzing protein structure, by computing various quantities such as the distribution of sizes of buried cavities and pockets, the scaling of solvent accessible area to the corresponding occupied volume, the composition of residues lining cavities, etc.

Nucleation and crystal growth in sheared granular sphere packings

We investigate the nucleation of ordered phases, their symmetries, and distributions in dense frictional hard sphere packings as a function of particle volume fraction ϕ, by imposing cyclic shear and constant applied pressure conditions. We show, with internal imaging, that the nucleating crystallites in the bulk consist of 10-60 spheres with hexagonal close packed (hcp) order and nonspherical shape, that are oriented preferentially along the shear axis. Above ϕ=0.62±0.005, crystallites with face centered cubic (fcc) order are observed with increasing probability, and ordered domains grow rapidly. A polycrystalline phase with domains of fcc and hcp order is observed after hundreds of thousands of shear cycles.

Onset of irreversibility in cyclic shear of granular packings

We investigate the onset of irreversibility in a dense granular medium subjected to cyclic shear in a split-bottom geometry. To probe the micro- and mesoscale, we image bead trajectories in three dimensions throughout a series of shear strain oscillations. Although beads lose and regain contact with neighbors during a cycle, the global topology of the contact network exhibits reversible properties for low oscillation amplitudes. With increasing reversal amplitude, a transition to an irreversible diffusive regime occurs.

Invited Article: Refractive index matched scanning of dense granular materials

We review an experimental method that allows to probe the time-dependent structure of fully three-dimensional densely packed granular materials and suspensions by means of particle recognition. The method relies on submersing a granular medium in a refractive index matched fluid. This makes the resulting suspension transparent. The granular medium is then visualized by exciting, layer by layer, the fluorescent dye in the fluid phase. We collect references and unreported experimental know-how to provide a solid background for future development of the technique, both for new and experienced users.

Packings of deformable spheres

We present an experimental study of disordered packings of deformable spheres. Fluorescent hydrogel spheres immersed in water together with a tomography technique enabled the imaging of the three-dimensional arrangement. The mechanical behavior of single spheres subjected to compression is first examined. Then the properties of packings of a randomized collection of deformable spheres in a box with a moving lid are tested. The transition to a state where the packing withstands finite stresses before yielding is observed. Starting from random packed states, the power law dependence of the normal force versus packing fraction or strain at different velocities is quantified. Furthermore, a compression-decompression sequence at low velocities resulted in rearrangements of the spheres. At larger packing fractions, a saturation of the mean coordination number took place, indicating the deformation and faceting of the spheres.

Path to fracture in granular flows: dynamics of contact networks

Capturing the dynamics of granular flows at intermediate length scales can often be difficult. We propose studying the dynamics of contact networks as a new tool to study fracture at intermediate scales. Using experimental three-dimensional flow fields with particle-scale resolution, we calculate the time evolving broken-links network and find that a giant component of this network is formed as shear is applied to this system. We implement a model of link breakages where the probability of a link breaking is proportional to the average rate of longitudinal strain (elongation) in the direction of the edge and find that the model demonstrates qualitative agreement with the data when studying the onset of the giant component. We note, however, that the broken-links network formed in the model is less clustered than our experimental observations, indicating that the model reflects less localized breakage events and does not fully capture the dynamics of the granular flow.

From frictional to viscous behavior: three-dimensional imaging and rheology of gravitational suspensions

We probe the three-dimensional flow structure and rheology of gravitational (nondensity matched) suspensions for a range of driving rates in a split-bottom geometry. We establish that for sufficiently slow flows, the suspension flows as if it were a dry granular medium, and confirm recent theoretical modeling on the rheology of split-bottom flows. For faster driving, the flow behavior is shown to be consistent with the rheological behavior predicted by the recently developed "inertial number" approaches for suspension flows.

Experimental study of random-close-packed colloidal particles

A collection of spherical particles can be packed tightly together into an amorphous packing known as "random close packing" (RCP). This structure is of interest as a model for the arrangement of molecules in simple liquids and glasses, as well as the arrangement of particles in sand piles. We use confocal microscopy to study the arrangement of colloidal particles in an experimentally realized RCP state. We image a large volume containing more than 450,000 particles with a resolution of each particle position to better than 0.02 particle diameters. While the arrangement of the particles satisfies multiple criteria for being random, we also observe a small fraction (less than 3%) of tiny crystallites (4 particles or fewer). These regions pack slightly better and are thus associated with locally higher densities. The structure factor of our sample at long length scales is nonzero, S(0)=0.049±0.008, suggesting that there are long wavelength density fluctuations in our sample. These may be due to polydispersity or tiny crystallites. Our results suggest that experimentally realizable RCP systems may be different from simulated RCP systems, in particular, with the presence of these long wavelength density fluctuations.

Spatial distribution functions of random packed granular spheres obtained by direct particle imaging

We measure the two-point density correlations and Voronoi cell distributions of cyclically sheared granular spheres obtained with a fluorescence technique and compare them with random packing of frictionless spheres. We find that the radial distribution function g(r) is captured by the Percus-Yevick equation for initial volume fraction ϕ=0.59. However, small but systematic deviations are observed because of the splitting of the second peak as ϕ is increased toward random close packing. The distribution of the Voronoi free volumes deviates from postulated Γ distributions, and the orientational order metric Q6 shows local order but no long range order. Overall, these measures show significant similarity of random packing of granular and frictionless spheres, but some systematic differences as well.

Contact network in nearly jammed disordered packings of hard-sphere chains

We present salient results of the analysis of the geometrical structure of a large fully equilibrated ensemble of nearly jammed packings of linear freely jointed chains of tangent hard spheres generated via extensive Monte Carlo simulations. In spite the expected differences due to chain connectivity, both the pair-correlation function and the contact network for chain packings are found to strongly resemble those in packings of monomeric hard spheres at the maximally random jammed (MRJ) state. A remarkable finding of the present work is the tendency of chains to form closed loops at the MRJ state as a consequence of chain collapse. Our simulations on disordered nearly jammed chain packings yield an average coordination number of 6, which fulfills the isostaticity condition and is in excellent agreement with the corresponding simulation [A. Donev, S. Torquato, and F. H. Stillinger, Phys. Rev. E 71, 011105 (2005)] and experimental [T. Aste, M. Saadatfar, and T. J. Senden, Phys. Rev. E 71, 061302 (2005)] findings for jammed packings of monatomic hard spheres. An exact correspondence between the statistical-mechanical ensembles of monomeric spheres and of hard-sphere chains offers insights regarding the structure and topology of the contact network of hard-sphere systems at the MRJ state.