Ductile-to-brittle transition and yielding in soft amorphous materials: perspectives and open questions.

Soft amorphous materials are viscoelastic solids ubiquitously found around us, from clays and cementitious pastes to emulsions and physical gels encountered in food or biomedical engineering. Under an external deformation, these materials undergo a noteworthy transition from a solid to a liquid state that reshapes the material microstructure. This yielding transition was the main theme of a workshop held from January 9 to 13, 2023 at the Lorentz Center in Leiden. The manuscript presented here offers a critical perspective on the subject, synthesizing insights from the various brainstorming sessions and informal discussions that unfolded during this week of vibrant exchange of ideas. The result of these exchanges takes the form of a series of open questions that represent outstanding experimental, numerical, and theoretical challenges to be tackled in the near future.

Self-Consistent Stochastic Dynamics for Finite-Size Networks of Spiking Neurons.

Despite the huge number of neurons composing a brain network, ongoing activity of local cell assemblies is intrinsically stochastic. Fluctuations in their instantaneous rate of spike firing ν(t) scale with the size of the assembly and persist in isolated networks, i.e., in the absence of external sources of noise. Although deterministic chaos due to the quenched disorder of the synaptic couplings underlies this seemingly stochastic dynamics, an effective theory for the network dynamics of a finite assembly of spiking neurons is lacking. Here, we fill this gap by extending the so-called population density approach including an activity- and size-dependent stochastic source in the Fokker-Planck equation for the membrane potential density. The finite-size noise embedded in this stochastic partial derivative equation is analytically characterized leading to a self-consistent and nonperturbative description of ν(t) valid for a wide class of spiking neuron networks. Power spectra of ν(t) are found in excellent agreement with those from detailed simulations both in the linear regime and across a synchronization phase transition, when a size-dependent smearing of the critical dynamics emerges.

Spatial population genetics with fluid flow.

The growth and evolution of microbial populations is often subjected to advection by fluid flows in spatially extended environments, with immediate consequences for questions of spatial population genetics in marine ecology, planktonic diversity and origin of life scenarios. Here, we review recent progress made in understanding this rich problem in the simplified setting of two competing genetic microbial strains subjected to fluid flows. As a pedagogical example we focus on antagonsim, i.e., two killer microorganism strains, each secreting toxins that impede the growth of their competitors (competitive exclusion), in the presence of stationary fluid flows. By solving two coupled reaction-diffusion equations that include advection by simple steady cellular flows composed of characteristic flow motifs in two dimensions (2D), we show how local flow shear and compressibility effects can interact with selective advantage to have a dramatic influence on genetic competition and fixation in spatially distributed populations. We analyze several 1D and 2D flow geometries including sources, sinks, vortices and saddles, and show how simple analytical models of the dynamics of the genetic interface can be used to shed light on the nucleation, coexistence and flow-driven instabilities of genetic drops. By exploiting an analogy with phase separation with nonconserved order parameters, we uncover how thesegeneticdrops harness fluid flows for novel evolutionary strategies, even in the presence of number fluctuations, as confirmed by agent-based simulations as well.

Self-similar properties of avalanche statistics in a simple turbulent model.

In this paper, we consider a simplified model of turbulence for large Reynolds numbers driven by a constant power energy input on large scales. In the statistical stationary regime, the behaviour of the kinetic energy is characterized by two well-defined phases: a laminar phase where the kinetic energy grows linearly for a (random) time [Formula: see text] followed by abrupt avalanche-like energy drops of sizes [Formula: see text] due to strong intermittent fluctuations of energy dissipation. We study the probability distribution [Formula: see text] and [Formula: see text] which both exhibit a quite well-defined scaling behaviour. Although [Formula: see text] and [Formula: see text] are not statistically correlated, we suggest and numerically checked that their scaling properties are related based on a simple, but non-trivial, scaling argument. We propose that the same approach can be used for other systems showing avalanche-like behaviour such as amorphous solids and seismic events. This article is part of the theme issue 'Scaling the turbulence edifice (part 1)'.

Fixation probabilities in weakly compressible fluid flows.

Competition between biological species in marine environments is affected by the motion of the surrounding fluid. An effective 2D compressibility can arise, for example, from the convergence and divergence of water masses at the depth at which passively traveling photosynthetic organisms are restricted to live. In this report, we seek to quantitatively study genetics under flow. To this end, we couple an off-lattice agent-based simulation of two populations in 1D to a weakly compressible velocity field-first a sine wave and then a shell model of turbulence. We find for both cases that even in a regime where the overall population structure is approximately unaltered, the flow can significantly diminish the effect of a selective advantage on fixation probabilities. We understand this effect in terms of the enhanced survival of organisms born at sources in the flow and the influence of Fisher genetic waves.

Stress Overshoots in Simple Yield Stress Fluids.

Soft glassy materials such as mayonnaise, wet clays, or dense microgels display a solid-to-liquid transition under external shear. Such a shear-induced transition is often associated with a nonmonotonic stress response in the form of a stress maximum referred to as "stress overshoot." This ubiquitous phenomenon is characterized by the coordinates of the maximum in terms of stress σ_{M} and strain γ_{M} that both increase as weak power laws of the applied shear rate. Here we rationalize such power-law scalings using a continuum model that predicts two different regimes in the limit of low and high applied shear rates. The corresponding exponents are directly linked to the steady-state rheology and are both associated with the nucleation and growth dynamics of a fluidized region. Our work offers a consistent framework for predicting the transient response of soft glassy materials upon startup of shear from the local flow behavior to the global rheological observables.

Rayleigh-Bénard convection of a model emulsion: anomalous heat-flux fluctuations and finite-size droplet effects.

We present mesoscale numerical simulations of Rayleigh-Bénard (RB) convection in a two-dimensional model emulsion. The systems under study are constituted of finite-size droplets, whose concentration Φ0 is systematically varied from small (Newtonian emulsions) to large values (non-Newtonian emulsions). We focus on the characterisation of the heat transfer properties close to the transition from conductive to convective states, where it is well known that a homogeneous Newtonian system exhibits a steady flow and a time-independent heat flux. In marked contrast, emulsions exhibit non-steady dynamics with fluctuations in the heat flux. In this paper, we aim at the characterisation of such non-steady dynamics via detailed studies on the time-averaged heat flux and its fluctuations. To quantitatively understand the time-averaged heat flux, we propose a side-by-side comparison between the emulsion system and a single-phase (SP) system, whose viscosity is suitably constructed from the shear rheology of the emulsion. We show that such local closure works well only when a suitable degree of coarse-graining (at the droplet scale) is introduced in the local viscosity. To delve deeper into the fluctuations in the heat flux, we furthermore propose a side-by-side comparison between a Newtonian emulsion (i.e., with a small droplet concentration) and a non-Newtonian emulsion (i.e., with a large droplet concentration), at fixed time-averaged heat flux. This comparison elucidates that finite-size droplets and the non-Newtonian rheology cooperate to trigger enhanced heat-flux fluctuations at the droplet scales. These enhanced fluctuations are rooted in the emergence of space correlations among distant droplets, which we highlight via direct measurements of the droplets displacement and the characterisation of the associated correlation function. The observed findings offer insights on heat transfer properties for confined systems possessing finite-size constituents.

A multi-component lattice Boltzmann approach to study the causality of plastic events.

Using a multi-component lattice Boltzmann (LB) model, we perform fluid kinetic simulations of confined and concentrated emulsions. The system presents the phenomenology of soft-glassy materials, including a Herschel-Bulkley rheology, yield stress, ageing and long relaxation time scales. Shearing the emulsion in a Couette cell below the yield stress results in plastic topological re-arrangement events which follow established empirical seismic statistical scaling laws, making this system a good candidate to study the physics of earthquakes. One characteristic of this model is the tendency for events to occur in avalanche clusters, with larger events, triggering subsequent re-arrangements. While seismologists have developed statistical tools to study correlations between events, a process to confirm causality remains elusive. We present here, a modification to our LB model, involving small, fast vibrations applied to individual droplets, effectively a macroscopic forcing, which results in the arrest of the topological plastic re-arrangements. This technique provides an excellent tool for identifying causality in plastic event clusters by examining the evolution of the dynamics after 'stopping' an event, and then checking which subsequent events disappear. This article is part of the theme issue 'Fluid dynamics, soft matter and complex systems: recent results and new methods'.

Unified Theoretical and Experimental View on Transient Shear Banding.

Dense emulsions, colloidal gels, microgels, and foams all display a solidlike behavior at rest characterized by a yield stress, above which the material flows like a liquid. Such a fluidization transition often consists of long-lasting transient flows that involve shear-banded velocity profiles. The characteristic time for full fluidization τ_{f} has been reported to decay as a power law of the shear rate γ[over ˙] and of the shear stress σ with respective exponents α and ß. Strikingly, the ratio of these exponents was empirically observed to coincide with the exponent of the Herschel-Bulkley law that describes the steady-state flow behavior of these complex fluids. Here we introduce a continuum model, based on the minimization of a "free energy," that captures quantitatively all the salient features associated with such transient shear banding. More generally, our results provide a unified theoretical framework for describing the yielding transition and the steady-state flow properties of yield stress fluids.

Avalanche statistics during coarsening dynamics.

We study the coarsening dynamics of a two-dimensional system via numerical simulations. The system under consideration is a biphasic system consisting of domains of a dispersed phase closely packed together in a continuous phase and separated by thin interfaces. Such a system is elastic and typically out of equilibrium. The equilibrium state is attained via the coarsening dynamics, wherein the dispersed phase slowly diffuses through the interfaces, causing the domains to change in size and eventually rearrange abruptly. The effect of rearrangements is propagated throughout the system via the intrinsic elastic interactions and may cause rearrangements elsewhere, resulting in intermittent bursts of activity and avalanche behaviour. Here we aim at quantitatively characterizing the corresponding avalanche statistics (i.e. size, duration, and inter-avalanche time). Despite the coarsening dynamics is triggered by an internal driving mechanism, we find quantitative indications that such avalanche statistics displays scaling-laws very similar to those observed in the response of disordered materials to external loads.

Economic Complexity: Correlations between Gross Domestic Product and Fitness.

In this paper we study the causal relation between country Economic Fitness F c and its Gross Domestic Product per capita ( G D P ). Using the Takens' theorem, as first suggested in (Sugihara, G. et al. 2012), we show that there exists a reasonable evidence of causal correlation between G D P and F c for relatively rich countries. This is not the case for relatively poor countries where F c and G D P do not show any significant causal relation. We also present some preliminary results to understand whether G D P or F c are driving factor for economic growth.

Competition between fast- and slow-diffusing species in non-homogeneous environments.

We study an individual-based model in which two spatially distributed species, characterized by different diffusivities, compete for resources. We consider three different ecological settings. In the first, diffusing faster has a cost in terms of reproduction rate. In the second case, resources are not uniformly distributed in space. In the third case, the two species are transported by a fluid flow. In all these cases, at varying the parameters, we observe a transition from a regime in which diffusing faster confers an effective selective advantage to one in which it constitutes a disadvantage. We analytically estimate the magnitude of this advantage (or disadvantage) and test it by measuring fixation probabilities in simulations of the individual-based model. Our results provide a framework to quantify evolutionary pressure for increased or decreased dispersal in a given environment.

Turbulence on a Fractal Fourier Set.

A novel investigation of the nature of intermittency in incompressible, homogeneous, and isotropic turbulence is performed by a numerical study of the Navier-Stokes equations constrained on a fractal Fourier set. The robustness of the energy transfer and of the vortex stretching mechanisms is tested by changing the fractal dimension D from the original three dimensional case to a strongly decimated system with D=2.5, where only about 3% of the Fourier modes interact. This is a unique methodology to probe the statistical properties of the turbulent energy cascade, without breaking any of the original symmetries of the equations. While the direct energy cascade persists, deviations from the Kolmogorov scaling are observed in the kinetic energy spectra. A model in terms of a correction with a linear dependency on the codimension of the fractal set E(k)∼k(-5/3+3-D) explains the results. At small scales, the intermittency of the vorticity field is observed to be quasisingular as a function of the fractal mode reduction, leading to an almost Gaussian statistics already at D∼2.98. These effects must be connected to a genuine modification in the triad-to-triad nonlinear energy transfer mechanism.

Selective advantage of diffusing faster.

We study a stochastic spatial model of biological competition in which two species have the same birth and death rates, but different diffusion constants. In the absence of this difference, the model can be considered as an off-lattice version of the voter model and presents similar coarsening properties. We show that even a relative difference in diffusivity on the order of a few percent may lead to a strong bias in the coarsening process favoring the more agile species. We theoretically quantify this selective advantage and present analytical formulas for the average growth of the fastest species and its fixation probability.

Spinodal decomposition in homogeneous and isotropic turbulence.

We study the competition between domain coarsening in a symmetric binary mixture below critical temperature and turbulent fluctuations. We find that the coarsening process is arrested in the presence of turbulence. The physics of the process shares remarkable similarities with the behavior of diluted turbulent emulsions and the arrest length scale can be estimated with an argument similar to the one proposed by Kolmogorov and Hinze for the maximal stability diameter of droplets in turbulence. Although, in the absence of flow, the microscopic diffusion constant is negative, turbulence does effectively arrest the inverse cascade of concentration fluctuations by making the low wavelength diffusion constant positive for scales above the Hinze length.

Population genetics in compressible flows.

We study competition between two biological species advected by a compressible velocity field. Individuals are treated as discrete Lagrangian particles that reproduce or die in a density-dependent fashion. In the absence of a velocity field and fitness advantage, number fluctuations lead to a coarsening dynamics typical of the stochastic Fisher equation. We investigate three examples of compressible advecting fields: a shell model of turbulence, a sinusoidal velocity field and a linear velocity sink. In all cases, advection leads to a striking drop in the fixation time, as well as a large reduction in the global carrying capacity. We find localization on convergence zones, and very rapid extinction compared to well-mixed populations. For a linear velocity sink, one finds a bimodal distribution of fixation times. The long-lived states in this case are demixed configurations with a single interface, whose location depends on the fitness advantage.

Strong noise limit for population dynamics in incompressible advection.

Genetic diversity is at the basis of the evolution process of populations and it is responsible for the populations' degree of fitness to a particular ecosystem. In marine environments many factors play a role in determining the dynamics of a population, including the amount of nutrients, the temperature, and many other stressing factors. An important and yet rather unexplored challenge is to figure out the role of individuals' dispersion, due to flow advection, on population genetics. In this paper we focus on two populations, one of which has a slight selective advantage, advanced by an incompressible two-dimensional flow. In particular, we want to understand how this advective flow can modify the dynamics of the advantageous allele. We generalize, through a theoretical analysis, previous evidence according to which the fixation probability is independent of diffusivity, showing that this is also independent of fluid advection. These findings may have important implications in the understanding of the dynamics of a population of microorganism, such as plankton or bacteria, in marine environments under the influence of (turbulent) currents.

Deep learning velocity signals allow quantifying turbulence intensity.

Turbulence, the ubiquitous and chaotic state of fluid motions, is characterized by strong and statistically nontrivial fluctuations of the velocity field, and it can be quantitatively described only in terms of statistical averages. Strong nonstationarities impede statistical convergence, precluding quantifying turbulence, for example, in terms of turbulence intensity or Reynolds number. Here, we show that by using deep neural networks, we can accurately estimate the Reynolds number within 15% accuracy, from a statistical sample as small as two large-scale eddy turnover times. In contrast, physics-based statistical estimators are limited by the convergence rate of the central limit theorem and provide, for the same statistical sample, at least a hundredfold larger error. Our findings open up previously unexplored perspectives and the possibility to quantitatively define and, therefore, study highly nonstationary turbulent flows as ordinarily found in nature and in industrial processes.

The collective effect of finite-sized inhomogeneities on the spatial spread of populations in two dimensions.

The dynamics of a population expanding into unoccupied habitat has been primarily studied for situations in which growth and dispersal parameters are uniform in space or vary in one dimension. Here, we study the influence of finite-sized individual inhomogeneities and their collective effect on front speed if randomly placed in a two-dimensional habitat. We use an individual-based model to investigate the front dynamics for a region in which dispersal or growth of individuals is reduced to zero (obstacles) or increased above the background (hotspots), respectively. In a regime where front dynamics is determined by a local front speed only, a principle of least time can be employed to predict front speed and shape. The resulting analytical solutions motivate an event-based algorithm illustrating the effects of several obstacles or hotspots. We finally apply the principle of least time to large heterogeneous environments by solving the Eikonal equation numerically. Obstacles lead to a slow-down that is dominated by the number density and width of obstacles, but not by their precise shape. Hotspots result in a speed-up, which we characterize as function of hotspot strength and density. Our findings emphasize the importance of taking the dimensionality of the environment into account.

Continuum modeling of shear startup in soft glassy materials.

Yield stress fluids (YSFs) display a dual nature highlighted by the existence of a critical stress σ_{y} such that YSFs are solid for stresses σ imposed below σ_{y}, whereas they flow like liquids for σ>σ_{y}. Under an applied shear rate γ[over ̇], the solid-to-liquid transition is associated with a complex spatiotemporal scenario that depends on the microscopic details of the system, on the boundary conditions, and on the system size. Still, the general phenomenology reported in the literature boils down to a simple sequence that can be divided into a short-time response characterized by the so-called "stress overshoot," followed by stress relaxation towards a steady state. Such relaxation can be either (1) long-lasting, which usually involves the growth of a shear band that can be only transient or that may persist at steady state or (2) abrupt, in which case the solid-to-liquid transition resembles the failure of a brittle material, involving avalanches. In the present paper, we use a continuum model based on a spatially resolved fluidity approach to rationalize the complete scenario associated with the shear-induced yielding of YSFs. A key feature of our model is to provide a scaling for the coordinates of the stress overshoot, i.e., stress σ_{M} and strain γ_{M} as a function of γ[over ̇], which shows good agreement with experimental and numerical data extracted from the literature. Moreover, our approach shows that the power-law scaling σ_{M}(γ[over ̇]) is intimately linked to the growth dynamics of a fluidized boundary layer in the vicinity of the moving boundary. Yet such scaling is independent of the fate of that layer, and of the long-term behavior of the YSF, i.e., whether the steady-state flow profile is homogeneous or shear-banded. Finally, when including the presence of "long-range" correlations, we show that our model displays a ductile to brittle transition, i.e., the stress overshoot reduces into a sharp stress drop associated with avalanches, which impacts the scaling σ_{M}(γ[over ̇]). This generalized model nicely captures subtle avalanche-like features of the transient shear banding dynamics reported in experiments. Our work offers a unified picture of shear-induced yielding in YSFs, whose complex spatiotemporal dynamics are deeply connected to nonlocal effects.