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.

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.

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.

Short-Time Transport Properties of Bidisperse Suspensions of Immunoglobulins and Serum Albumins Consistent with a Colloid Physics Picture.

The crowded environment of biological systems such as the interior of living cells is occupied by macromolecules with a broad size distribution. This situation of polydispersity might influence the dependence of the diffusive dynamics of a given tracer macromolecule in a monodisperse solution on its hydrodynamic size and on the volume fraction. The resulting size dependence of diffusive transport crucially influences the function of a living cell. Here, we investigate a simplified model system consisting of two constituents in aqueous solution, namely, of the proteins bovine serum albumin (BSA) and bovine polyclonal gamma-globulin (Ig), systematically depending on the total volume fraction and ratio of these constituents. From high-resolution quasi-elastic neutron spectroscopy, the separate apparent short-time diffusion coefficients for BSA and Ig in the mixture are extracted, which show substantial deviations from the diffusion coefficients measured in monodisperse solutions at the same total volume fraction. These deviations can be modeled quantitatively using results from the short-time rotational and translational diffusion in a two-component hard sphere system with two distinct, effective hydrodynamic radii. Thus, we find that a simple colloid picture well describes short-time diffusion in binary mixtures as a function of the mixing ratio and the total volume fraction. Notably, the self-diffusion of the smaller protein BSA in the mixture is faster than the diffusion in a pure BSA solution, whereas the self-diffusion of Ig in the mixture is slower than in the pure Ig solution.

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.

Adsorption and Desorption of Polymers on Bioinspired Chemically Structured Substrates.

Natural biological surfaces exhibit interesting properties due to their inhomogeneous chemical and physical structure at the micro- and nanoscale. In the case of hair or skin, this also influences how waterborne macromolecules ingredients will adsorb and form cosmetically performing deposits (i.e., shampoos, cleansers, etc.). Here, we study the adsorption of hydrophilic flexible homopolymers on heterogeneous, chemically patterned substrates that represent the surface of the hair by employing coarse-grained molecular dynamics simulations. We develop a method in which the experimental images of the substrate are used to obtain information about the surface properties. We investigate the polymer adsorption as a function of polymer chain length and polymer concentration spanning both dilute and semidilute regimes. Adsorbed structures are quantified in terms of trains, loops, and tails. We show that upon increasing polymer concentration, the length of tails and loops increases at the cost of monomers belonging to trains. Furthermore, using an effective description, we probe the stability of the resulting adsorbed structures under a linear shear flow. Our work is a first step toward developing models of complex macromolecules interacting with realistic biological surfaces, as needed for the development of more ecofriendly industrial products.

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.

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.

Protein Short-Time Diffusion in a Naturally Crowded Environment.

The interior of living cells is a dense and polydisperse suspension of macromolecules. Such a complex system challenges an understanding in terms of colloidal suspensions. As a fundamental test we employ neutron spectroscopy to measure the diffusion of tracer proteins (immunoglobulins) in a cell-like environment (cell lysate) with explicit control over crowding conditions. In combination with Stokesian dynamics simulation, we address protein diffusion on nanosecond time scales where hydrodynamic interactions dominate over negligible protein collisions. We successfully link the experimental results on these complex, flexible molecules with coarse-grained simulations providing a consistent understanding by colloid theories. Both experiments and simulations show that tracers in polydisperse solutions close to the effective particle radius Reff = ⟨ Ri3⟩1/3 diffuse approximately as if the suspension was monodisperse. The simulations further show that macromolecules of sizes R > Reff ( R < Reff) are slowed more (less) effectively even at nanosecond time scales, which is highly relevant for a quantitative understanding of cellular processes.

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.

Correlation and shear bands in a plastically deformed granular medium.

Recent experiments (Le Bouil et al., Phys. Rev. Lett., 2014, 112, 246001) have analyzed the statistics of local deformation in a granular solid undergoing plastic deformation. Experiments report strongly anisotropic correlation between events, with a characteristic angle that was interpreted using elasticity theory and the concept of Eshelby transformations with dilation; interestingly, the shear bands that characterize macroscopic failure occur at an angle that is different from the one observed in microscopic correlations. Here, we interpret this behavior using a mesoscale elastoplastic model of solid flow that incorporates a local Mohr-Coulomb failure criterion. This differs from the interpretation of Le Bouil et al., which is based on purely elastic considerations ignoring the potential role of local friction on deformation patterns. We show that the angle observed in the microscopic correlations can be understood by combining the elastic interactions associated with Eshelby transformation with the local failure criterion. At large strains, we also induce permanent shear bands at an angle that is different from the one observed in the correlation pattern. We interpret this angle as the one that leads to the maximal instability of slip lines.

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.

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.

Inertia and universality of avalanche statistics: The case of slowly deformed amorphous solids.

By means of a finite elements technique we solve numerically the dynamics of an amorphous solid under deformation in the quasistatic driving limit. We study the noise statistics of the stress-strain signal in the steady-state plastic flow, focusing on systems with low internal dissipation. We analyze the distributions of avalanche sizes and durations and the density of shear transformations when varying the damping strength. In contrast to avalanches in the overdamped case, dominated by the yielding point universal exponents, inertial avalanches are controlled by a nonuniversal damping-dependent feedback mechanism, eventually turning negligible the role of correlations. Still, some general properties of avalanches persist and new scaling relations can be proposed.

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.

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.