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.

On De Gennes narrowing of fluids confined at the molecular scale in nanoporous materials.

Beyond well-documented confinement and surface effects arising from the large internal surface and severely confining porosity of nanoporous hosts, the transport of nanoconfined fluids remains puzzling in many aspects. With striking examples such as memory, i.e., non-viscous effects, intermittent dynamics, and surface barriers, the dynamics of fluids in nanoconfinement challenge classical formalisms (e.g., random walk, viscous/advective transport)-especially for molecular pore sizes. In this context, while molecular frameworks such as intermittent Brownian motion, free volume theory, and surface diffusion are available to describe the self-diffusion of a molecularly confined fluid, a microscopic theory for collective diffusion (i.e., permeability), which characterizes the flow induced by a thermodynamic gradient, is lacking. Here, to fill this knowledge gap, we invoke the concept of "De Gennes narrowing," which relates the wavevector-dependent collective diffusivity D0(q) to the fluid structure factor S(q). First, using molecular simulation for a simple yet representative fluid confined in a prototypical solid (zeolite), we unravel an essential coupling between the wavevector-dependent collective diffusivity and the structural ordering imposed on the fluid by the crystalline nanoporous host. Second, despite this complex interplay with marked Bragg peaks in the fluid structure, the fluid collective dynamics is shown to be accurately described through De Gennes narrowing. Moreover, in contrast to the bulk fluid, the departure from De Gennes narrowing for the confined fluid in the macroscopic limit remains small as the fluid/solid interactions in severe confinement screen collective effects and, hence, weaken the wavevector dependence of collective transport.

Structural and Dynamical Properties of Elastin-Like Peptides near Their Lower Critical Solution Temperature.

Elastin-like peptides (ELPs) are artificially derived intrinsically disordered proteins (IDPs) mimicking the hydrophobic repeat unit in the protein elastin. ELPs are characterized by a lower critical solution temperature (LCST) in aqueous media. Here, we investigate the sequence GVG(VPGVG)3 over a wide range of temperatures (below, around, and above the LCST) and peptide concentrations employing all-atom molecular dynamics simulations, where we focus on the role of intra- and interpeptide interactions. We begin by investigating the structural properties of a single peptide that demonstrates a hydrophobic collapse with temperature, albeit moderate, because the sequence length is short. We observe a change in the interaction between two peptides from repulsive to attractive with temperature by evaluating the potential of mean force, indicating an LCST-like behavior. Next, we explore dynamical and structural properties of peptides in multichain systems. We report the formation of dynamical aggregates with coil-like conformation, in which valine central residues play an important role. Moreover, the lifetime of contacts between chains strongly depends on the temperature and can be described by a power-law decay that is consistent with the LCST-like behavior. Finally, the peptide translational and internal motion are slowed by an increase in the peptide concentration and temperature.

Permanent shear localization in dense disordered materials due to microscopic inertia.

In this work using computer simulations of 3D model of dense disordered solids we show, for the first time, the appearance of shear localization in the stationary flow under homogeneous driving conditions. To rationalize our simulation results we develop a continuum model, that couples the dynamics of the local flow to the evolution of a kinetic temperature field related to the local inertial dynamics. Our model predicts that the coupling of the flow field to this additional destabilizing field appears only as a necessary condition for shear localization, a minimum system size is necessary to accommodate the flow instability. Moreover we show that this size criterion resulting from our continuum description is in quantitative agreement with our particle-based simulation results.

Elastic avalanches reveal marginal behavior in amorphous solids.

Mechanical deformation of amorphous solids can be described as consisting of an elastic part in which the stress increases linearly with strain, up to a yield point at which the solid either fractures or starts deforming plastically. It is well established, however, that the apparent linearity of stress with strain is actually a proxy for a much more complex behavior, with a microscopic plasticity that is reflected in diverging nonlinear elastic coefficients. Very generally, the complex structure of the energy landscape is expected to induce a singular response to small perturbations. In the athermal quasistatic regime, this response manifests itself in the form of a scale-free plastic activity. The distribution of the corresponding avalanches should reflect, according to theoretical mean-field calculations [S. Franz and S. Spigler, Phys. Rev. E 95, 022139 (2017)], the geometry of phase space in the vicinity of a typical local minimum. In this work, we characterize this distribution for simple models of glass-forming systems, and we find that its scaling is compatible with the mean-field predictions for systems above the jamming transition. These systems exhibit marginal stability, and scaling relations that hold in the stationary state are examined and confirmed in the elastic regime. By studying the respective influence of system size and age, we suggest that marginal stability is systematic in the thermodynamic limit.

Are strongly confined colloids good models for two dimensional liquids?

Quasi-two-dimensional (quasi-2D) colloidal hard-sphere suspensions confined in a slit geometry are widely used as two-dimensional (2D) model systems in experiments that probe the glassy relaxation dynamics of 2D systems. However, the question to what extent these quasi-2D systems indeed represent 2D systems is rarely brought up. Here, we use computer simulations that take into account hydrodynamic interactions to show that dense quasi-2D colloidal bi-disperse hard-sphere suspensions exhibit much more rapid diffusion and relaxation than their 2D counterparts at the same area fraction. This difference is induced by the additional vertical space in the quasi-2D samples in which the small colloids can move out of the 2D plane, therefore allowing overlap between particles in the projected trajectories. Surprisingly, this difference in the dynamics can be accounted for if, instead of using the surface density, one characterizes the systems by means of a suitable structural quantity related to the radial distribution function. This implies that in the two geometries, the relevant physics for glass formation is essentially identical. Our results provide not only practical implications on 2D colloidal experiments but also interesting insights into the 3D-to-2D crossover in glass-forming systems.

Temperature dependence of thermodynamic, dynamical, and dielectric properties of water models.

We investigate the temperature dependence of thermodynamic (density and isobaric heat capacity), dynamical (self-diffusion coefficient and shear viscosity), and dielectric properties of several water models, such as the commonly employed TIP3P water model, the well-established four-point water model TIP4P-2005, and the recently developed four-point water model TIP4P-D. We focus on the temperature range of interest for the field of computational biophysics and soft matter (280-350 K). The four-point water models lead to a spectacularly improved agreement with experimental data, strongly suggesting that the use of more modern parameterizations should be favored compared to the more traditional TIP3P for modeling temperature-dependent phenomena in biomolecular systems.

Time correlation functions for quantum systems: Validating Bayesian approaches for harmonic oscillators and beyond.

The quantum harmonic oscillator is the fundamental building block to compute thermal properties of virtually any dielectric crystal at low temperatures in terms of phonons, extended further to cases with anharmonic couplings, or even disordered solids. In general, Path Integral Monte Carlo or Path Integral Molecular Dynamics methods are powerful tools to determine stochastically thermodynamic quantities without systematic bias, not relying on perturbative schemes. Addressing transport properties, for instance calculating thermal conductivity from PIMC, however, is substantially more difficult. Although correlation functions of current operators can be determined by PIMC from analytic continuation on the imaginary time axis, Bayesian methods are usually employed for the numerical inversion back to real time response functions. This task not only strongly relies on the accuracy of the PIMC data but also introduces noticeable dependence on the model used for the inversion. Here, we address both difficulties with care. In particular, we first devise improved estimators for current correlations, which substantially reduce the variance of the PIMC data. Next, we provide a neat statistical approach to the inversion problem, blending into a fresh workflow the classical stochastic maximum entropy method together with recent notions borrowed from statistical learning theory. We test our ideas on a single harmonic oscillator and a collection of oscillators with a continuous distribution of frequencies and provide indications of the performance of our method in the case of a particle in a double well potential. This work establishes solid grounds for an unbiased, fully quantum mechanical calculation of transport properties in solids.

Local versus Global Stretched Mechanical Response in a Supercooled Liquid near the Glass Transition.

Amorphous materials have a rich relaxation spectrum, which is usually described in terms of a hierarchy of relaxation mechanisms. In this work, we investigate the local dynamic modulus spectra in a model glass just above the glass transition temperature by performing a mechanical spectroscopy analysis with molecular dynamics simulations. We find that the spectra, at the local as well as on the global scale, can be well described by the Cole-Davidson formula in the frequency range explored with simulations. Surprisingly, the Cole-Davidson stretching exponent does not change with the size of the local region that is probed. The local relaxation time displays a broad distribution, as expected based on dynamic heterogeneity concepts, but the stretching is obtained independently of this distribution. We find that the size dependence of the local relaxation time and moduli can be well explained by the elastic shoving model.

Mean-Field Scenario for the Athermal Creep Dynamics of Yield-Stress Fluids.

We develop a theoretical description based on an existent mean-field model for the transient dynamics prior to the steady flow of yielding materials. The mean-field model not only reproduces the experimentally observed nonlinear time dependence of the shear-rate response to an external stress, but also allows for the determination of the different physical processes involved in the onset of the reacceleration phase after the initial slowing down and a distinct fluidization phase. The fluidization time displays a power-law dependence on the distance of the applied stress to an age-dependent yield stress, which is not universal but strongly dependent on initial conditions.

Creep dynamics of athermal amorphous materials: a mesoscopic approach.

Yield stress fluids display complex dynamics, in particular when driven into the transient regime between the solid and the flowing state. Inspired by creep experiments on dense amorphous materials, we implement mesoscale elasto-plastic descriptions to analyze such transient dynamics in athermal systems. Both our mean-field and space-dependent approaches consistently reproduce the typical experimental strain rate responses to different applied steps in stress. Moreover, they allow us to understand basic processes involved in the strain rate slowing down (creep) and the strain rate acceleration (fluidization) phases. The fluidization time increases in a power-law fashion as the applied external stress approaches a static yield stress. This stress value is related to the stress over-shoot in shear start-up experiments, and it is known to depend on sample preparation and age. By calculating correlations of the accumulated plasticity in the spatially resolved model, we reveal different modes of cooperative motion during the creep dynamics.

Role of the Intercrystalline Tie Chains Network in the Mechanical Response of Semicrystalline Polymers.

We examine the microscopic origin of the tensile response in semicrystalline polymers by performing large-scale molecular dynamics simulations of various chain lengths. We investigate the microscopic rearrangements of the polymers during tensile deformation and show that the intercrystalline chain connections known as tie chains contribute significantly to the elastic and plastic response. These results suggest that the mechanical behavior of semicrystalline polymers is controlled by two interpenetrated networks of entanglements and tie chains.

Nonlinear Rheology in a Model Biological Tissue.

The rheological response of dense active matter is a topic of fundamental importance for many processes in nature such as the mechanics of biological tissues. One prominent way to probe mechanical properties of tissues is to study their response to externally applied forces. Using a particle-based model featuring random apoptosis and environment-dependent division rates, we evidence a crossover from linear flow to a shear-thinning regime with an increasing shear rate. To rationalize this nonlinear flow we derive a theoretical mean-field scenario that accounts for the interplay of mechanical and active noise in local stresses. These noises are, respectively, generated by the elastic response of the cell matrix to cell rearrangements and by the internal activity.

Acoustic excitations and elastic heterogeneities in disordered solids.

In the recent years, much attention has been devoted to the inhomogeneous nature of the mechanical response at the nanoscale in disordered solids. Clearly, the elastic heterogeneities that have been characterized in this context are expected to strongly affect the nature of the sound waves which, in contrast to the case of perfect crystals, cannot be completely rationalized in terms of phonons. Building on previous work on a toy model showing an amorphization transition, we investigate the relationship between sound waves and elastic heterogeneities in a unified framework by continuously interpolating from the perfect crystal, through increasingly defective phases, to fully developed glasses. We provide strong evidence of a direct correlation between sound wave features and the extent of the heterogeneous mechanical response at the nanoscale.

Effects of Inertia on the Steady-Shear Rheology of Disordered Solids.

We study the finite-shear-rate rheology of disordered solids by means of molecular dynamics simulations in two dimensions. By systematically varying the damping strength ζ in the low-temperature limit, we identify two well-defined flow regimes, separated by a thin (temperature-dependent) crossover region. In the overdamped regime, the athermal rheology is governed by the competition between elastic forces and viscous forces, whose ratio gives the Weissenberg number Wi∝ζγ[over ˙]; the macroscopic stress Σ follows the frequently encountered Herschel-Bulkley law Σ=Σ_{0}+ksqrt[Wi], with yield stress Σ_{0}>0. In the underdamped (inertial) regime, dramatic changes in the rheology are observed for low damping: the flow curve becomes nonmonotonic. This change is not caused by longer-lived correlations in the particle dynamics at lower damping; instead, for weak dissipation, the sample heats up considerably due to, and in proportion to, the driving. By thermostating more or less underdamped systems, we are able to link quantitatively the rheology to the kinetic temperature and the shear rate, rescaled with Einstein's vibration frequency.

Driving Rate Dependence of Avalanche Statistics and Shapes at the Yielding Transition.

We study stress time series caused by plastic avalanches in athermally sheared disordered materials. Using particle-based simulations and a mesoscopic elastoplastic model, we analyze system size and shear-rate dependence of the stress-drop duration and size distributions together with their average temporal shape. We find critical exponents different from mean-field predictions, and a clear asymmetry for individual avalanches. We probe scaling relations for the rate dependency of the dynamics and we report a crossover towards mean-field results for strong driving.

Analytical correlation functions for motion through diffusivity landscapes.

Diffusion of a particle through an energy and diffusivity landscape is a very general phenomenon in numerous systems of soft and condensed matter. On the one hand, theoretical frameworks such as Langevin and Fokker-Planck equations present valuable accounts to understand these motions in great detail, and numerous studies have exploited these approaches. On the other hand, analytical solutions for correlation functions, as, e.g., desired by experimentalists for data fitting, are only available for special cases. We explore the possibility to use different theoretical methods in the specific picture of time-dependent switching between diffusive states to derive analytical functions that allow to link experimental and simulation results to theoretical calculations. In particular, we present a closed formula for diffusion switching between two states, as well as a general recipe of how to generalize the formula to multiple states.

Simulation of entangled polymer solutions.

We present a computer simulation of entangled polymer solutions at equilibrium. The chains repel each other via a soft Gaussian potential, appropriate for semi-dilute solutions at the scale of a correlation blob. The key innovation to suppress chain crossings is to use a pseudo-continuous model of a backbone which effectively leaves no gaps between consecutive points on the chain, unlike the usual bead-and-spring model. Our algorithm is sufficiently fast to observe the entangled regime using a standard desktop computer. The simulated structural and mechanical correlations are in fair agreement with the expected predictions for a semi-dilute solution of entangled chains.

Edwards Thermodynamics for a Driven Athermal System with Dry Friction.

We obtain, using semianalytical transfer operator techniques, the Edwards thermodynamics of a one-dimensional model of blocks connected by harmonic springs and subjected to dry friction. The theory is able to reproduce the linear divergence of the correlation length as a function of energy density observed in direct numerical simulations of the model under tapping dynamics. We further characterize analytically this divergence using a Gaussian approximation for the distribution of mechanically stable configurations, and show that it is related to the existence of a peculiar infinite temperature critical point.

On the relevance of disorder in athermal amorphous materials under shear.

We show that, at least at a mean-field level, the effect of structural disorder in sheared amorphous media is very dissimilar depending on the thermal or athermal nature of their underlying dynamics. We first introduce a toy model, including explicitly two types of noise (thermal versus athermal). Within this interpretation framework, we argue that mean-field athermal dynamics can be accounted for by the so-called Hébraud-Lequeux (HL) model, in which the mechanical noise stems explicitly from the plastic activity in the sheared medium. Then, we show that the inclusion of structural disorder, by means of a distribution of yield energy barriers, has no qualitative effect in the HL model, while such a disorder is known to be one of the key ingredients leading kinematically to a finite macroscopic yield stress in other mean-field descriptions, such as the Soft-Glassy-Rheology model. We conclude that the statistical mechanisms at play in the emergence of a macroscopic yield stress, and a complex stationary dynamics at low shear rate, are different in thermal and athermal amorphous systems.