Evolution of force networks during stick-slip motion of an intruder in a granular material: Topological measures extracted from experimental data.

In quasi-two-dimensional experiments with photoelastic particles confined to an annular region, an intruder constrained to move in a circular path halfway between the annular walls experiences stick-slip dynamics. We discuss the response of the granular medium to the driven intruder, focusing on the evolution of the force network during sticking periods. Because the available experimental data do not include precise information about individual contact forces, we use an approach developed in our previous work [Basak et al., J. Eng. Mech. 147, 04021100 (2021)0733-939910.1061/(ASCE)EM.1943-7889.0002003] based on networks constructed from measurements of the integrated strain magnitude on each particle. These networks are analyzed using topological measures based on persistence diagrams, revealing that force networks evolve smoothly but in a nontrivial manner throughout each sticking period, even though the intruder and granular particles are stationary. Characteristic features of persistence diagrams show identifiable slip precursors. In particular, the number of generators describing the structure and complexity of force networks increases consistently before slips. Key features of the dynamics are similar for granular materials composed of disks or pentagons, but some details are consistently different. In particular, we find significantly larger fluctuations of the measures computed based on persistence diagrams and, therefore, of the underlying networks, for systems of pentagonal particles.

Universal features of the stick-slip dynamics of an intruder moving through a confined granular medium.

Experiments and simulations of an intruder dragged by a spring through a two-dimensional annulus of granular material exhibit robust force fluctuations. At low packing fractions (ϕ<ϕ_{0}), the intruder clears an open channel. Above ϕ_{0}, stick-slip dynamics develop, with an average energy release that is independent of the particle-particle and particle-base friction coefficients but does depend on the width W of the annulus and the diameter D of the intruder. A simple model predicts the dependence of ϕ_{0} on W and D, allowing for a data collapse for the average energy release as a function of ϕ/ϕ_{0}. These results pose challenges for theories of mechanical failure in amorphous materials.

Stress propagation in locally loaded packings of disks and pentagons.

The mechanical strength and flow of granular materials can depend strongly on the shapes of individual grains. We report quantitative results obtained from photoelasticimetry experiments on locally loaded, quasi-two-dimensional granular packings of either disks or pentagons exhibiting stick-slip dynamics. Packings of pentagons resist the intruder at significantly lower packing fractions than packings of disks, transmitting stresses from the intruder to the boundaries over a smaller spatial extent. Moreover, packings of pentagons feature significantly fewer back-bending force chains than packings of disks. Data obtained on the forward spatial extent of stresses and back-bending force chains collapse when the packing fraction is rescaled according to the packing fraction of steady state open channel formation, though data on intruder forces and dynamics do not collapse. We comment on the influence of system size on these findings and highlight connections with the dynamics of the disks and pentagons during slip events.

Static and dynamic features of granular material failure due to upward pulling of a buried sphere by a slowly increasing force.

A spherical intruder embedded in a confined granular column is extracted by pulling it upward by an attached string. As the tension of the string gradually increases, a failure event occurs at a certain pulling force, leading to rapid upward acceleration of the intruder. The threshold force and the dynamics of the failure event are experimentally investigated for different filling heights and column diameters, using Ottawa sand and glass beads. For the Ottawa sand, we find that the failure force can be fit by a model describing the weight of the granular material in a cone with the vertex at the bottom of the intruder and a vertex angle of 72°. The agreement between the model and experiments is good for heights less than the column (tube) diameter, but measured values deviate from the model for larger heights. We also report on experiments with glass beads that reveal unexpected effects for relatively small ratios of tube diameters to grain size. The dynamics of the intruder during the failure event is studied using high-speed video analysis. The granular drag force monotonically decays during the pullout for sufficiently large tube diameters. In narrow columns, a monotonic decay of drag force after failure is observed for low heights, whereas a secondary peak can be seen in sufficiently deep and narrow columns, indicating the existence of different mechanisms of failure. The normalized drag force declines with intruder displacement closely for all tube diameters within small displacements.

Random logic networks: From classical Boolean to quantum dynamics.

We investigate dynamical properties of a quantum generalization of classical reversible Boolean networks. The state of each node is encoded as a single qubit, and classical Boolean logic operations are supplemented by controlled bit-flip and Hadamard operations. We consider synchronous updating schemes in which each qubit is updated at each step based on stored values of the qubits from the previous step. We investigate the periodic or quasiperiodic behavior of quantum networks, and we analyze the propagation of single site perturbations through the quantum networks with input degree one. A nonclassical mechanism for perturbation propagation leads to substantially different evolution of the Hamming distance between the original and perturbed states.

Yielding, rigidity, and tensile stress in sheared columns of hexapod granules.

Granular packings of nonconvex or elongated particles can form freestanding structures like walls or arches. For some particle shapes, such as staples, the rigidity arises from interlocking of pairs of particles, but the origins of rigidity for noninterlocking particles remains unclear. We report on experiments and numerical simulations of sheared columns of "hexapods," particles consisting of three mutually orthogonal sphero-cylinders whose centers coincide. We vary the length-to-diameter aspect ratio, α, of the sphero-cylinders and subject the packings to quasistatic direct shear. For small α, we observe a finite yield stress. For large α, however, the column becomes rigid when sheared, supporting stresses that increase sharply with increasing strain. Analysis of x-ray microcomputed tomography (micro-CT) data collected during the shear reveals that the stiffening is associated with a tilted, oblate cluster of hexapods near the nominal shear plane in which particle deformation and average contact number both increase. Simulation results show that the particles are collectively under tension along one direction, even though they do not interlock pairwise. These tensions comes from contact forces carrying large torques, and they are perpendicular to the compressive stresses in the packing. They counteract the tendency to dilate, thus stabilizing the particle cluster.

Intruder in a two-dimensional granular system: Effects of dynamic and static basal friction on stick-slip and clogging dynamics.

We present simulation results for an intruder pulled through a two-dimensional granular system by a spring using a model designed to mimic the experiments described by Kozlowski et al. [Phys. Rev. E 100, 032905 (2019)2470-004510.1103/PhysRevE.100.032905]. In that previous study the presence of basal friction between the grains and the base was observed to change the intruder dynamics from clogging to stick-slip. Here we first show that our simulation results are in excellent agreement with the experimental data for a variety of experimentally accessible friction coefficients governing interactions of particles with each other and with boundaries. We then use simulations to explore a broader range of parameter space, focusing on the friction between the particles and the base. We consider both static and dynamic basal friction coefficients, which are difficult to vary smoothly in experiments. The simulations show that dynamic friction strongly affects the stick-slip behavior when the coefficient is decreased below 0.1, while static friction plays only a marginal role.

Dispersion tuning and route reconfiguration of acoustic waves in valley topological phononic crystals.

The valley degree of freedom in crystals offers great potential for manipulating classical waves, however, few studies have investigated valley states with complex wavenumbers, valley states in graded systems, or dispersion tuning for valley states. Here, we present tunable valley phononic crystals (PCs) composed of hybrid channel-cavity cells with three tunable parameters. Our PCs support valley states and Dirac cones with complex wavenumbers. They can be configured to form chirped valley PCs in which edge modes are slowed to zero group velocity states, where the energy at different frequencies accumulates at different designated locations. They enable multiple functionalities, including tuning of dispersion relations for valley states, robust routing of surface acoustic waves, and spatial modulation of group velocities. This work may spark future investigations of topological states with complex wavenumbers in other classical systems, further study of topological states in graded materials, and the development of acoustic devices.

Shear-Jammed, Fragile, and Steady States in Homogeneously Strained Granular Materials.

We study the jamming phase diagram of sheared granular material using a novel Couette shear setup with a multiring bottom. The setup uses small basal friction forces to apply a volume-conserving linear shear with no shear band to a granular system composed of frictional photoelastic discs. The setup can generate arbitrarily large shear strain due to its circular geometry, and the shear direction can be reversed, allowing us to measure a feature that distinguishes shear-jammed from fragile states. We report systematic measurements of the stress, strain, and contact network structure at phase boundaries that have been difficult to access by traditional experimental techniques, including the yield stress curve and the jamming curve close to ϕ_{SJ}≈0.75, the smallest packing fraction supporting a shear-jammed state. We observe fragile states created under large shear strain over a range of ϕ<ϕ_{SJ}. We also find a transition in the character of the quasistatic steady flow centered around ϕ_{SJ} on the yield curve as a function of packing fraction. Near ϕ_{SJ}, the average contact number, fabric anisotropy, and nonrattler fraction all show a change of slope. Above ϕ_{F}≈0.7 the steady flow shows measurable deviations from the basal linear shear profile, and above ϕ_{b}≈0.78 the flow is localized in a shear band.

Dynamics of a grain-scale intruder in a two-dimensional granular medium with and without basal friction.

We report on a series of experiments in which a grain-sized intruder is pushed by a spring through a two-dimensional granular material composed of photoelastic disks in a Couette geometry. We study the intruder dynamics as a function of packing fraction for two types of supporting substrates: A frictional glass plate and a layer of water for which basal friction forces are negligible. We observe two dynamical regimes: Intermittent flow, in which the intruder moves freely most of the time but occasionally gets stuck, and stick-slip dynamics, in which the intruder advances via a sequence of distinct, rapid events. When basal friction is present, we observe a smooth crossover between the two regimes as a function of packing fraction, and we find that reducing the interparticle friction coefficient causes the stick-slip regime to shift to higher packing fractions. When basal friction is eliminated, we observe intermittent flow at all accessible packing fractions. For all cases, we present results for the statistics of stick events, the intruder velocity, and the force exerted on the intruder by the grains. Our results indicate the qualitative importance of basal friction at high packing fractions and suggest a possible connection between intruder dynamics in a static material and clogging dynamics in granular flows.

Seismicity in sheared granular matter.

We report on experiments investigating the dynamics of a slider that is pulled by a spring across a granular medium consisting of a vertical layer of photoelastic disks. The motion proceeds through a sequence of discrete events, analogous to seismic shocks, in which elastic energy stored in the spring is rapidly released. We measure the statistics of several properties of the individual events: the energy loss in the spring, the duration of the movement, and the temporal profile of the slider motion. We also study certain conditional probabilities and the statistics of mainshock-aftershock sequences. At low driving rates, we observe crackling with Omori-Utsu, Båth, and waiting time laws similar to those observed in seismic dynamics. At higher driving rates, where the sequence of events shows strong periodicity, we observe scaling laws and asymmetrical event shapes that are clearly distinguishable from those in the crackling regime.

Hyperuniformity and anti-hyperuniformity in one-dimensional substitution tilings.

This work considers the scaling properties characterizing the hyperuniformity (or anti-hyperuniformity) of long-wavelength fluctuations in a broad class of one-dimensional substitution tilings. A simple argument is presented which predicts the exponent α governing the scaling of Fourier intensities at small wavenumbers, tilings with α > 0 being hyperuniform, and numerical computations confirm that the predictions are accurate for quasiperiodic tilings, tilings with singular continuous spectra and limit-periodic tilings. Quasiperiodic or singular continuous cases can be constructed with α arbitrarily close to any given value between -1 and 3. Limit-periodic tilings can be constructed with α between -1 and 1 or with Fourier intensities that approach zero faster than any power law.

Phase diagram and aggregation dynamics of a monolayer of paramagnetic colloids.

We have developed a tunable colloidal system and a corresponding theoretical model for studying the phase behavior of particles assembling under the influence of long-range magnetic interactions. A monolayer of paramagnetic particles is subjected to a spatially uniform magnetic field with a static perpendicular component and a rapidly rotating in-plane component. The sign and strength of the interactions vary with the tilt angle θ of the rotating magnetic field. For a purely in-plane field, θ=90^{∘}, interactions are attractive and the experimental results agree well with both equilibrium and out-of-equilibrium predictions based on a two-body interaction model. For tilt angles 50^{∘}≲θ≲55^{∘}, the two-body interaction gives a short-range attractive and long-range repulsive interaction, which predicts the formation of equilibrium microphases. In experiments, however, a different type of assembly is observed. Inclusion of three-body (and higher-order) terms in the model does not resolve the discrepancy. We further characterize the anomalous regime by measuring the time-dependent cluster size distribution.

Formation of limit-periodic structures by quadrupole particles confined to a triangular lattice.

We have performed Monte Carlo (MC) simulations on two-dimensional systems of quadrupole particles confined to a triangular lattice in order to determine the conditions that permit the formation of a limit-periodic phase. We have found that limit-periodic structures form only when the rotations of the particles are confined to a set of six orientations aligned with the lattice directions. Related structures including striped and unidirectional rattler phases form when π/π66 rotations or continuous rotations are allowed. Order parameters signaling the formation of the limit-periodic structure and related structures are measured as a function of temperature. Our findings on the formation of the limit-periodic structure elucidate features relevant to the experimental creation of such a structure, which is expected to have interesting vibrational and electromagnetic modes.

Hard sphere packings within cylinders.

Arrangements of identical hard spheres confined to a cylinder with hard walls have been used to model experimental systems, such as fullerenes in nanotubes and colloidal wire assembly. Finding the densest configurations, called close packings, of hard spheres of diameter σ in a cylinder of diameter D is a purely geometric problem that grows increasingly complex as D/σ increases, and little is thus known about the regime for D > 2.873σ. In this work, we extend the identification of close packings up to D = 4.00σ by adapting Torquato-Jiao's adaptive-shrinking-cell formulation and sequential-linear-programming (SLP) technique. We identify 17 new structures, almost all of them chiral. Beyond D ≈ 2.85σ, most of the structures consist of an outer shell and an inner core that compete for being close packed. In some cases, the shell adopts its own maximum density configuration, and the stacking of core spheres within it is quasiperiodic. In other cases, an interplay between the two components is observed, which may result in simple periodic structures. In yet other cases, the very distinction between the core and shell vanishes, resulting in more exotic packing geometries, including some that are three-dimensional extensions of structures obtained from packing hard disks in a circle.

Phase transformations in binary colloidal monolayers.

Phase transformations can be difficult to characterize at the microscopic level due to the inability to directly observe individual atomic motions. Model colloidal systems, by contrast, permit the direct observation of individual particle dynamics and of collective rearrangements, which allows for real-space characterization of phase transitions. Here, we study a quasi-two-dimensional, binary colloidal alloy that exhibits liquid-solid and solid-solid phase transitions, focusing on the kinetics of a diffusionless transformation between two crystal phases. Experiments are conducted on a monolayer of magnetic and nonmagnetic spheres suspended in a thin layer of ferrofluid and exposed to a tunable magnetic field. A theoretical model of hard spheres with point dipoles at their centers is used to guide the choice of experimental parameters and characterize the underlying materials physics. When the applied field is normal to the fluid layer, a checkerboard crystal forms; when the angle between the field and the normal is sufficiently large, a striped crystal assembles. As the field is slowly tilted away from the normal, we find that the transformation pathway between the two phases depends strongly on crystal orientation, field strength, and degree of confinement of the monolayer. In some cases, the pathway occurs by smooth magnetostrictive shear, while in others it involves the sudden formation of martensitic plates.

Regulatory logic and pattern formation in the early sea urchin embryo.

We model the endomesoderm tissue specification process in the vegetal half of the early sea urchin embryo using Boolean models with continuous-time updating to represent the regulatory network that controls gene expression. Our models assume that the network interaction rules remain constant over time and the dynamics plays out on a predetermined program of cell divisions. An exhaustive search of two-node models, in which each node may represent a module of several genes in the real regulatory network, yields a unique network architecture that can accomplish the pattern formation task at hand--the formation of three latitudinal tissue bands from an initial state with only two distinct cell types. Analysis of an eight-gene model constructed from available experimental data reveals that it has a modular structure equivalent to the successful two-node case. Our results support the hypothesis that the gene regulatory network provides sufficient instructions for producing the correct pattern of tissue specification at this stage of development (between the fourth and tenth cleavages in the urchin embryo).

Emergence of limit-periodic order in tiling models.

A two-dimensional (2D) lattice model defined on a triangular lattice with nearest- and next-nearest-neighbor interactions based on the Taylor-Socolar monotile is known to have a limit-periodic ground state. The system reaches that state during a slow quench through an infinite sequence of phase transitions. We study the model as a function of the strength of the next-nearest-neighbor interactions and introduce closely related 3D models with only nearest-neighbor interactions that exhibit limit-periodic phases. For models with no next-nearest-neighbor interactions of the Taylor-Socolar type, there is a large degenerate class of ground states, including crystalline patterns and limit-periodic ones, but a slow quench still yields the limit-periodic state. For the Taylor-Socolar lattic model, we present calculations of the diffraction pattern for a particular decoration of the tile that permits exact expressions for the amplitudes and identify domain walls that slow the relaxation times in the ordered phases. For one of the 3D models, we show that the phase transitions are first order, with equilibrium structures that can be more complex than in the 2D case, and we include a proof of aperiodicity for a geometrically simple tile with only nearest-neighbor matching rules.

Delayed transition to new cell fates during cellular reprogramming.

In many embryos specification toward one cell fate can be diverted to a different cell fate through a reprogramming process. Understanding how that process works will reveal insights into the developmental regulatory logic that emerged from evolution. In the sea urchin embryo, cells at gastrulation were found to reprogram and replace missing cell types after surgical dissections of the embryo. Non-skeletogenic mesoderm (NSM) cells reprogrammed to replace missing skeletogenic mesoderm cells and animal caps reprogrammed to replace all endomesoderm. In both cases evidence of reprogramming onset was first observed at the early gastrula stage, even if the cells to be replaced were removed earlier in development. Once started however, the reprogramming occurred with compressed gene expression dynamics. The NSM did not require early contact with the skeletogenic cells to reprogram, but the animal cap cells gained the ability to reprogram early in gastrulation only after extended contact with the vegetal halves prior to that time. If the entire vegetal half was removed at early gastrula, the animal caps reprogrammed and replaced the vegetal half endomesoderm. If the animal caps carried morpholinos to either hox11/13b or foxA (endomesoderm specification genes), the isolated animal caps failed to reprogram. Together these data reveal that the emergence of a reprogramming capability occurs at early gastrulation in the sea urchin embryo and requires activation of early specification components of the target tissues.

Causal structure of oscillations in gene regulatory networks: Boolean analysis of ordinary differential equation attractors.

A common approach to the modeling of gene regulatory networks is to represent activating or repressing interactions using ordinary differential equations for target gene concentrations that include Hill function dependences on regulator gene concentrations. An alternative formulation represents the same interactions using Boolean logic with time delays associated with each network link. We consider the attractors that emerge from the two types of models in the case of a simple but nontrivial network: a figure-8 network with one positive and one negative feedback loop. We show that the different modeling approaches give rise to the same qualitative set of attractors with the exception of a possible fixed point in the ordinary differential equation model in which concentrations sit at intermediate values. The properties of the attractors are most easily understood from the Boolean perspective, suggesting that time-delay Boolean modeling is a useful tool for understanding the logic of regulatory networks.