Soft matter physics of the ground beneath our feet.

The soft part of the Earth's surface - the ground beneath our feet - constitutes the basis for life and natural resources, yet a general physical understanding of the ground is still lacking. In this critical time of climate change, cross-pollination of scientific approaches is urgently needed to better understand the behavior of our planet's surface. The major topics in current research in this area cross different disciplines, spanning geosciences, and various aspects of engineering, material sciences, physics, chemistry, and biology. Among these, soft matter physics has emerged as a fundamental nexus connecting and underpinning many research questions. This perspective article is a multi-voice effort to bring together different views and approaches, questions and insights, from researchers that work in this emerging area, the soft matter physics of the ground beneath our feet. In particular, we identify four major challenges concerned with the dynamics in and of the ground: (I) modeling from the grain scale, (II) near-criticality, (III) bridging scales, and (IV) life. For each challenge, we present a selection of topics by individual authors, providing specific context, recent advances, and open questions. Through this, we seek to provide an overview of the opportunities for the broad Soft Matter community to contribute to the fundamental understanding of the physics of the ground, strive towards a common language, and encourage new collaborations across the broad spectrum of scientists interested in the matter of the Earth's surface.

Development of broad modulus profile upon polymer-polymer interface formation between immiscible glassy-rubbery domains.

Interfaces of glassy materials such as thin films, blends, and composites create strong unidirectional gradients to the local heterogeneous dynamics that can be used to elucidate the length scales and mechanisms associated with the dynamic heterogeneity of glasses. We focus on bilayer films of two different polymers with very different glass transition temperatures ([Formula: see text]) where previous work has demonstrated a long-range (∼200 nm) profile in local [Formula: see text] is established between immiscible glassy and rubbery polymer domains when the polymer-polymer interface is formed to equilibrium. Here, we demonstrate that an equally long-ranged gradient in local modulus [Formula: see text] is established when the polymer-polymer interface ([Formula: see text]5 nm) is formed between domains of glassy polystyrene (PS) and rubbery poly(butadiene) (PB), consistent with previous reports of a broad [Formula: see text] profile in this system. A continuum physics model for the shear wave propagation caused by a quartz crystal microbalance across a PB/PS bilayer film is used to measure the viscoelastic properties of the bilayer during the evolution of the PB/PS interface showing the development of a broad gradient in local modulus [Formula: see text] spanning [Formula: see text]180 nm between the glassy and rubbery domains of PS and PB. We suggest these broad profiles in [Formula: see text] and [Formula: see text] arise from a coupling of the spectrum of vibrational modes across the polymer-polymer interface as a result of acoustic impedance matching of sound waves with [Formula: see text] nm during interface broadening that can then trigger density fluctuations in the neighboring domain.

Ticks' attraction to electrically charged hosts.

Ticks are blood-feeding parasites with limited locomotion, known for transmitting multiple pathogens to vertebrates. England et al. suggest that ticks can be easily pulled, via electrostatic induction, toward charged hosts with fluffy coats that are prone to accumulate higher electrostatic potentials. Thus, static electricity may influence ticks' ecology and management.

Extracting forces from noisy dynamics in dusty plasmas.

Extracting environmental forces from noisy data is a common yet challenging task in complex physical systems. Machine learning (ML) represents a robust approach to this problem, yet is mostly tested on simulated data with known parameters. Here we use supervised ML to extract the electrostatic, dissipative, and stochastic forces acting on micron-sized charged particles levitated in an argon plasma (dusty plasma). By tracking the subpixel motion of particles in subsequent images, we successfully estimated these forces from their random motion. The experiments contained important sources of non-Gaussian noise, such as drift and pixel locking, representing a data mismatch from methods used to analyze simulated data with purely Gaussian noise. Our model was trained on simulated particle trajectories that included all of these artifacts, and used more than 100 dynamical and statistical features, resulting in a prediction with 50% better accuracy than conventional methods. Finally, in systems with two interacting particles, the model provided noncontact measurements of the particle charge and Debye length in the plasma environment.

The lifetime of charged dust in the atmosphere.

Wind-blown dust plays a critical role in numerous geophysical and biological systems, yet current models fail to explain the transport of coarse-mode particles (>5 µm) to great distances from their sources. For particles larger than a few microns, electrostatic effects have been invoked to account for longer-than-predicted atmospheric residence times. Although much effort has focused on elucidating the charging processes, comparatively little effort has been expended understanding the stability of charge on particles once electrified. Overall, electrostatic-driven transport requires that charge remain present on particles for days to weeks. Here, we present a set of experiments designed to explore the longevity of electrostatic charge on levitated airborne particles after a single charging event. Using an acoustic levitator, we measured the charge on particles of different material compositions suspended in atmospheric conditions for long periods of time. In dry environments, the total charge on particles decayed in over 1 week. The decay timescale decreased to days in humid environments. These results were independent of particle material and charge polarity. However, exposure to UV radiation could both increase and decrease the decay time depending on polarity. Our work suggests that the rate of charge decay on airborne particles is solely determined by ion capture from the air. Furthermore, using a one-dimensional sedimentation model, we predict that atmospheric dust of order 10 µm will experience the largest change in residence time due to electrostatic forces.

Minimum Leidenfrost Temperature on Smooth Surfaces.

During the Leidenfrost effect, a thin insulating vapor layer separates an evaporating liquid from a hot solid. Here we demonstrate that Leidenfrost vapor layers can be sustained at much lower temperatures than those required for formation. Using a high-speed electrical technique to measure the thickness of water vapor layers over smooth, metallic surfaces, we find that the explosive failure point is nearly independent of material and fluid properties, suggesting a purely hydrodynamic mechanism determines this threshold. For water vapor layers of several millimeters in size, the minimum temperature for stability is ≈140 °C, corresponding to an average vapor layer thickness of 10-20 µm.

Metabolic excretion associated with nutrient-growth dysregulation promotes the rapid evolution of an overt metabolic defect.

In eukaryotes, conserved mechanisms ensure that cell growth is coordinated with nutrient availability. Overactive growth during nutrient limitation ("nutrient-growth dysregulation") can lead to rapid cell death. Here, we demonstrate that cells can adapt to nutrient-growth dysregulation by evolving major metabolic defects. Specifically, when yeast lysine-auxotrophic mutant lys- encountered lysine limitation, an evolutionarily novel stress, cells suffered nutrient-growth dysregulation. A subpopulation repeatedly evolved to lose the ability to synthesize organosulfurs (lys-orgS-). Organosulfurs, mainly reduced glutathione (GSH) and GSH conjugates, were released by lys- cells during lysine limitation when growth was dysregulated, but not during glucose limitation when growth was regulated. Limiting organosulfurs conferred a frequency-dependent fitness advantage to lys-orgS- by eliciting a proper slow growth program, including autophagy. Thus, nutrient-growth dysregulation is associated with rapid organosulfur release, which enables the selection of organosulfur auxotrophy to better tune cell growth to the metabolic environment. We speculate that evolutionarily novel stresses can trigger atypical release of certain metabolites, setting the stage for the evolution of new ecological interactions.

Pore-size dependence and slow relaxation of hydrogel friction on smooth surfaces.

Hydrogels consist of a cross-linked polymer matrix imbibed with a solvent such as water at volume fractions that can exceed 90%. They are important in many scientific and engineering applications due to their tunable physiochemical properties, biocompatibility, and ultralow friction. Their multiphase structure leads to a complex interfacial rheology, yet a detailed, microscopic understanding of hydrogel friction is still emerging. Using a custom-built tribometer, here we identify three distinct regimes of frictional behavior for polyacrylic acid (PAA), polyacrylamide (PAAm), and agarose hydrogel spheres on smooth surfaces. We find that at low velocities, friction is controlled by hydrodynamic flow through the porous hydrogel network and is inversely proportional to the characteristic pore size. At high velocities, a mesoscopic, lubricating liquid film forms between the gel and surface that obeys elastohydrodynamic theory. Between these regimes, the frictional force decreases by an order of magnitude and displays slow relaxation over several minutes. Our results can be interpreted as an interfacial shear thinning of the polymers with an increasing relaxation time due to the confinement of entanglements. This transition can be tuned by varying the solvent salt concentration, solvent viscosity, and sliding geometry at the interface.

Tuning Contact Line Dynamics and Deposition Patterns in Volatile Liquid Mixtures.

The spreading of a pure, volatile liquid on a wettable substrate has been studied in extensive detail. Here we show that the addition of a miscible, nonvolatile liquid can strongly alter the contact line dynamics and the final liquid deposition pattern. We observe two distinct regimes of behavior depending on the relative strength of solutal Marangoni forces and surface wetting. Fingerlike instabilities precede the deposition of a submicron thick film for large Marangoni forces and small solute contact angles, whereas isolated pearl-like drops emerge and are deposited in quasicrystalline patterns for small Marangoni forces and large solute contact angles. This behavior can be tuned by directly varying the contact angle of the solute liquid on the solid substrate.

High-throughput quantification of microbial birth and death dynamics using fluorescence microscopy.

BACKGROUND: Microbes live in dynamic environments where nutrient concentrations fluctuate. Quantifying fitness in terms of birth rate and death rate in a wide range of environments is critical for understanding microbial evolution and ecology. METHODS: Here, using high-throughput time-lapse microscopy, we have quantified how Saccharomyces cerevisiae mutants incapable of synthesizing an essential metabolite (auxotrophs) grow or die in various concentrations of the required metabolite. We establish that cells normally expressing fluorescent proteins lose fluorescence upon death and that the total fluorescence in an imaging frame is proportional to the number of live cells even when cells form multiple layers. We validate our microscopy approach of measuring birth and death rates using flow cytometry, cell counting, and chemostat culturing. RESULTS: For lysine-requiring cells, very low concentrations of lysine are not detectably consumed and do not support cell birth, but delay the onset of death phase and reduce the death rate compared to no lysine. In contrast, in low hypoxanthine, hypoxanthine-requiring cells can produce new cells, yet also die faster than in the absence of hypoxanthine. For both strains, birth rates under various metabolite concentrations are better described by the sigmoidal-shaped Moser model than the well-known Monod model, while death rates can vary with metabolite concentration and time. CONCLUSIONS: Our work reveals how time-lapse microscopy can be used to discover non-intuitive microbial birth and death dynamics and to quantify growth rates in many environments.

Universal scaling of polygonal desiccation crack patterns.

Polygonal desiccation crack patterns are commonly observed in natural systems. Despite their quotidian nature, it is unclear whether similar crack patterns which span orders of magnitude in length scales share the same underlying physics. In thin films, the characteristic length of polygonal cracks is known to monotonically increase with the film thickness; however, existing theories that consider the mechanical, thermodynamic, hydrodynamic, and statistical properties of cracking often lead to contradictory predictions. Here we experimentally investigate polygonal cracks in drying suspensions of micron-sized particles by varying film thickness, boundary adhesion, packing fraction, and solvent. Although polygonal cracks were observed in most systems above a critical film thickness, in cornstarch-water mixtures, multiscale crack patterns were observed due to two distinct desiccation mechanisms. Large-scale, primary polygons initially form due to capillary-induced film shrinkage, whereas small-scale, secondary polygons appear later due to the deswelling of the hygroscopic particles. In addition, we find that the characteristic area of the polygonal cracks, A_{p}, obeys a universal power law, A_{p}=αh^{4/3}, where h is the film thickness. By quantitatively linking α with the material properties during crack formation, we provide a robust framework for understanding multiscale polygonal crack patterns from microscopic to geologic scales.

Precursors to Molecular Slip on Smooth Hydrophobic Surfaces.

Experiments and simulations suggest that simple liquids may experience slip while flowing near a smooth, hydrophobic surface. Here we show how precursors to molecular slip can be observed in the complex response of a liquid to oscillatory shear. We measure both the change in frequency and bandwidth of a quartz crystal microbalance during the growth of a single drop of water immersed in an ambient liquid. By varying the hydrophobicity of the surface using self-assembled monolayers, our results show little or no slip for water on all surfaces. However, we observe excess transverse motion near hydrophobic surfaces due to weak binding in the corrugated surface potential, an essential precursor to slip. We also show how this effect can be easily missed in simulations utilizing finite-ranged interaction potentials.

Quantifying flow and stress in ice mélange, the world's largest granular material.

Tidewater glacier fjords are often filled with a collection of calved icebergs, brash ice, and sea ice. For glaciers with high calving rates, this "mélange" of ice can be jam-packed, so that the flow of ice fragments is mostly determined by granular interactions. In the jammed state, ice mélange has been hypothesized to influence iceberg calving and capsize, dispersion and attenuation of ocean waves, injection of freshwater into fjords, and fjord circulation. However, detailed measurements of ice mélange are lacking due to difficulties in instrumenting remote, ice-choked fjords. Here we characterize the flow and associated stress in ice mélange, using a combination of terrestrial radar data, laboratory experiments, and numerical simulations. We find that, during periods of terminus quiescence, ice mélange experiences laminar flow over timescales of hours to days. The uniform flow fields are bounded by shear margins along fjord walls where force chains between granular icebergs terminate. In addition, the average force per unit width that is transmitted to the glacier terminus, which can exceed 107 N/m, increases exponentially with the mélange length-to-width ratio. These "buttressing" forces are sufficiently high to inhibit the initiation of large-scale calving events, supporting the notion that ice mélange can be viewed as a weak granular ice shelf that transmits stresses from fjord walls back to glacier termini.

Emergent Bistability and Switching in a Nonequilibrium Crystal.

Multistability is an inseparable feature of many physical, chemical, and biological systems which are driven far from equilibrium. In these nonequilibrium systems, stochastic dynamics often induces switching between distinct states on emergent time scales; for example, bistable switching is a natural feature of noisy, spatially extended systems that consist of bistable elements. Nevertheless, here we present experimental evidence that bistable elements are not required for the global bistability of a system. We observe temporal switching between a crystalline, condensed state and a gaslike, excited state in a spatially extended, quasi-two-dimensional system of charged microparticles. Accompanying numerical simulations show that conservative forces, damping, and stochastic noise are sufficient to prevent steady-state equilibrium, leading to switching between the two states over a range of time scales, from seconds to hours.

Casimir effect between pinned particles in two-dimensional jammed systems.

The Casimir effect arises when long-ranged fluctuations are geometrically confined between two surfaces, leading to a macroscopic force. Traditionally, these forces have been observed in quantum systems and near critical points in classical systems. Here we show the existence of Casimir-like forces between two pinned particles immersed in two-dimensional systems near the jamming transition. We observe two components to the total force: a short-ranged, depletion force and a long-ranged, repulsive Casimir-like force. The Casimir-like force dominates as the jamming transition is approached, and when the pinned particles are much larger than the ambient jammed particles. We show that this repulsive force arises due to a clustering of particles with strong contact forces around the perimeter of the pinned particles. As the separation between the pinned particles decreases, a region of high-pressure develops between them, leading to a net repulsive force.

Echoes from anharmonic normal modes in model glasses.

Glasses display a wide array of nonlinear acoustic phenomena at temperatures T ≲ 1 K. This behavior has traditionally been explained by an ensemble of weakly coupled, two-level tunneling states, a theory that is also used to describe the thermodynamic properties of glasses at low temperatures. One of the most striking acoustic signatures in this regime is the existence of phonon echoes, a feature that has been associated with two-level systems with the same formalism as spin echoes in NMR. Here we report the existence of a distinctly different type of acoustic echo in classical models of glassy materials. Our simulations consist of finite-ranged, repulsive spheres and also particles with attractive forces using Lennard-Jones interactions. We show that these echoes are due to anharmonic, weakly coupled vibrational modes and perhaps provide an alternative explanation for the phonon echoes observed in glasses at low temperatures.

Coalescence of bubbles and drops in an outer fluid.

When two liquid drops touch, a microscopic connecting liquid bridge forms and rapidly grows as the two drops merge into one. Whereas coalescence has been thoroughly studied when drops coalesce in vacuum or air, many important situations involve coalescence in a dense surrounding fluid, such as oil coalescence in brine. Here we study the merging of gas bubbles and liquid drops in an external fluid. Our data indicate that the flows occur over much larger length scales in the outer fluid than inside the drops themselves. Thus, we find that the asymptotic early regime is always dominated by the viscosity of the drops, independent of the external fluid. A phase diagram showing the crossovers into the different possible late-time dynamics identifies a dimensionless number that signifies when the external viscosity can be important.

Energy loss at propagating jamming fronts in granular gas clusters.

We explore the initial moments of impact between two dense granular clusters in a two-dimensional geometry. The particles are composed of solid CO(2) and are levitated on a hot surface. Upon collision, the propagation of a dynamic "jamming front" produces a distinct regime for energy dissipation in a granular gas in which the translational kinetic energy decreases by over 90%. Experiments and associated simulations show that the initial loss of kinetic energy obeys a power law in time ΔE = -Kt(3/2), a form that can be predicted from kinetic arguments.

Collision dynamics of particle clusters in a two-dimensional granular gas.

In a granular gas, inelastic collisions produce an instability in which the constituent particles cluster heterogeneously. These clusters then interact with each other, further decreasing their kinetic energy. We report experiments of the free collisions of dense clusters of particles in a two-dimensional geometry. The particles are composed of solid CO(2), which float nearly frictionlessly on a hot surface due to sublimated vapor. After two dense clusters of ≈100 particles collide, there are two distinct stages of evolution. First, the translational kinetic energy rapidly decreases by over 90% as a "jamming front" sweeps across each cluster. Subsequently, the kinetic energy decreases more slowly as the particles approach the container boundaries. In this regime, the measured velocity distributions are non-Gaussian with long tails. Finally, we compare our experiments to computer simulations of colliding, two-dimensional, granular clusters composed of circular, viscoelastic particles with friction.

The inexorable resistance of inertia determines the initial regime of drop coalescence.

Drop coalescence is central to diverse processes involving dispersions of drops in industrial, engineering, and scientific realms. During coalescence, two drops first touch and then merge as the liquid neck connecting them grows from initially microscopic scales to a size comparable to the drop diameters. The curvature of the interface is infinite at the point where the drops first make contact, and the flows that ensue as the two drops coalesce are intimately coupled to this singularity in the dynamics. Conventionally, this process has been thought to have just two dynamical regimes: a viscous and an inertial regime with a cross-over region between them. We use experiments and simulations to reveal that a third regime, one that describes the initial dynamics of coalescence for all drop viscosities, has been missed. An argument based on force balance allows the construction of a new coalescence phase diagram.