Emergent solidity of amorphous materials as a consequence of mechanical self-organisation

Amorphous GeSnSe nanoparticles as a Li-Ion battery anode with High-Capacity and long cycle performance

A new ternary amorphous GeSnSe (GSS) nanopowder was effectively synthesized by using ball milling under inert atmosphere. Its topographical, microstructural and elemental characterizations revealed the formation of nanoparticles with undefined shape, short-range order and the tailored stoichiometry. Remarkably, this novel amorphous material demonstrates its competences as a promising Li-ion host anode, exhibiting a high cycle performance with a specific charge capacity of 963 mAh g-1 at an applied C-rate of 0.2C with a coulombic efficiency > 99.4 % after 300 cycles. Its high specific capacity, large rate capability, acceptable capacity retention and long cycle life could be attributed to a dual Li-ion storage mechanism that consists mostly of multiple reversible electrochemical processes as conversion and alloying reactions and capacitive processes. Moreover, its stable volume expansion (34 %), moderate electrode polarization (248.9 mV), reasonable charge transfer resistance (83 Ω) and apparent Li-ion diffusion coefficients between 10-9 - 10-14 cm2 s-1 could be promoted by a synergistic effect between Ge (capacity), Sn (conductivity) and Se (stability), which plays an important role on the stability and high cycle performance of the promising GSS-based anode.

Long-range correlations in elastic moduli and local stresses at the unjamming transition

We explore the behaviour of spatially heterogeneous elastic moduli as well as the correlations between local moduli in model solids with short-range repulsive potentials. We show through numerical simulations that local elastic moduli exhibit long-range correlations, similar to correlations in the local stresses. Specifically, the correlations in local shear moduli exhibit anisotropic behavior at large lengthscales characterized by pinch-point singularities in Fourier space, displaying a structural pattern akin to shear stress correlations. Focussing on two-dimensional jammed solids approaching the unjamming transition, we show that stress correlations exhibit universal properties, characterized by a quadratic p2 dependence of the correlations as the pressure p approaches zero, independent of the details of the model. In contrast, the modulus correlations exhibit a power-law dependence with different exponents depending on the specific interaction potential. Furthermore, we illustrate that while affine responses lack long-range correlations, the total modulus, which encompasses non-affine behavior, exhibits long-range correlations.

Short-range correlation of stress chains near solid-to-liquid transition in active monolayers

Using a three-dimensional model of cell monolayers, we study the spatial organization of active stress chains as the monolayer transitions from a solid to a liquid state. The critical exponents that characterize this transition map the isotropic stress percolation onto the two-dimensional random percolation universality class, suggesting short-range stress correlations near this transition. This mapping is achieved via two distinct, independent pathways: (i) cell-cell adhesion and (ii) active traction forces. We unify our findings by linking the nature of this transition to high-stress fluctuations, distinctly linked to each pathway. The results elevate the importance of the transmission of mechanical information in dense active matter and provide a new context for understanding the non-equilibrium statistical physics of phase transition in active systems.

Universal stress correlations in crystalline and amorphous packings

We present a universal characterization of stress correlations in athermal systems, across crystalline to amorphous packings. Via numerical analysis of static configurations of particles interacting through harmonic as well as Lennard-Jones potentials, for a variety of preparation protocols and ranges of microscopic disorder, we show that the properties of the stress correlations at large lengthscales are surprisingly universal across all situations, independent of structural correlations, or the correlations in orientational order. In the near-crystalline limit, we present exact results for the stress correlations for both models, which work surprisingly well at large lengthscales, even in the amorphous phase. Finally, we study the differences in stress fluctuations across the amorphization transition, where stress correlations reveal the loss of periodicity in the structure at short lengthscales with increasing disorder.

Viscoelastic relaxation and topological fluctuations in glass-forming liquids

A method for characterizing the topological fluctuations in liquids is proposed. This approach exploits the concept of the weighted gyration tensor of a collection of particles and permits the definition of a local configurational unit (LCU). The first principal axis of the gyration tensor serves as the director of the LCU, which can be tracked and analyzed by molecular dynamics simulations. Analysis of moderately supercooled Kob-Andersen mixtures suggests that orientational relaxation of the LCU closely follows viscoelastic relaxation and exhibits a two-stage behavior. The slow relaxing component of the LCU corresponds to the structural, Maxwellian mechanical relaxation. Additionally, it is found that the mean curvature of the LCUs is approximately zero at the Maxwell relaxation time with the Gaussian curvature being negative. This observation implies that structural relaxation occurs when the configurationally stable and destabilized regions interpenetrate each other in a bicontinuous manner. Finally, the mean and Gaussian curvatures of the LCUs can serve as reduced variables for the shear stress correlation, providing a compelling proof of the close connection between viscoelastic relaxation and topological fluctuations in glass-forming liquids.

Solid-that-Flows Picture of Glass-Forming Liquids

This perspective article reviews arguments that glass-forming liquids are different from those of standard liquid-state theory, which typically have a viscosity in the mPa·s range and relaxation times on the order of picoseconds. These numbers grow dramatically and become 1012 - 1015 times larger for liquids cooled toward the glass transition. This translates into a qualitative difference, and below the "solidity length" which is roughly one micron at the glass transition, a glass-forming liquid behaves much like a solid. Recent numerical evidence for the solidity of ultraviscous liquids is reviewed, and experimental consequences are discussed in relation to dynamic heterogeneity, frequency-dependent linear-response functions, and the temperature dependence of the average relaxation time.

Control of Mechanical and Fracture Properties in Two-Phase Materials Reinforced by Continuous, Irregular Networks

Composites with high strength and high fracture resistance are desirable for structural and protective applications. Most composites, however, suffer from poor damage tolerance and are prone to unpredictable fractures. Understanding the behavior of materials with an irregular reinforcement phase offers fundamental guidelines for tailoring their performance. Here, the fracture nucleation and propagation in two phase composites, as a function of the topology of their irregular microstructures is studied. A stochastic algorithm is used to design the polymeric reinforcing network, achieving independent control of topology and geometry of the microstructure. By tuning the local connectivity of isodense tiles and their assembly into larger structures, the mechanical and fracture properties of the architected composites are tailored at the local and global scale. Finally, combining different reinforcing networks into a spatially determined meso-scale assembly, it is demonstrated how the spatial propagation of fracture in architected composite materials can be designed and controlled a priori.

Emergence of rigidity percolation in flowing granular systems

Jammed granular media and glasses exhibit spatial long-range correlations as a result of mechanical equilibrium. However, the existence of such correlations in the flowing matter, where the mechanical equilibrium is unattainable, has remained elusive. Here, we investigate this problem in the context of the percolation of interparticle forces in flowing granular media. We find that the flow rate introduces an effective long-range correlation, which plays the role of a relevant perturbation giving rise to a spectrum of varying exponents on a critical line as a function of the flow rate. Our numerical simulations along with analytical arguments predict a crossover flow rate [Formula: see text] below which the effect of induced disorder is weak and the universality of the force chain structure is shown to be given by the standard rigidity percolation. We also find a power-law behavior for the critical exponents with the flow rate [Formula: see text].

Local and global measures of the shear moduli of jammed disk packings

Strain-controlled isotropic compression gives rise to jammed packings of repulsive, frictionless disks with either positive or negative global shear moduli. We carry out computational studies to understand the contributions of the negative shear moduli to the mechanical response of jammed disk packings. We first decompose the ensemble-averaged, global shear modulus as 〈G〉=(1-F_{-})〈G_{+}〉+F_{-}〈G_{-}〉, where F_{-} is the fraction of jammed packings with negative shear moduli and 〈G_{+}〉 and 〈G_{-}〉 are the average values from packings with positive and negative moduli, respectively. We show that 〈G_{+}〉 and 〈|G_{-}|〉 obey different power-law scaling relations above and below pN^{2}∼1. For pN^{2}>1, both 〈G_{+}〉N and 〈|G_{-}|〉N∼(pN^{2})^{β}, where β∼0.5 for repulsive linear spring interactions. Despite this, 〈G〉N∼(pN^{2})^{β^{'}} with β^{'}≳0.5 due to the contributions from packings with negative shear moduli. We show further that the probability distribution of global shear moduli P(G) collapses at fixed pN^{2} and different values of p and N. We calculate analytically that P(G) is a Γ distribution in the pN^{2}≪1 limit. As pN^{2} increases, the skewness of P(G) decreases and P(G) becomes a skew-normal distribution with negative skewness in the pN^{2}≫1 limit. We also partition jammed disk packings into subsystems using Delaunay triangulation of the disk centers to calculate local shear moduli. We show that the local shear moduli defined from groups of adjacent triangles can be negative even when G>0. The spatial correlation function of local shear moduli C(r[over ⃗]) displays weak correlations for pn_{sub}^{2}<10^{-2}, where n_{sub} is the number of particles within each subsystem. However, C(r[over ⃗]) begins to develop long-ranged spatial correlations with fourfold angular symmetry for pn_{sub}^{2}≳10^{-2}.

Elasticity, Facilitation, and Dynamic Heterogeneity in Glass-Forming Liquids

We study the role of elasticity-induced facilitation on the dynamics of glass-forming liquids by a coarse-grained two-dimensional model in which local relaxation events, taking place by thermal activation, can trigger new relaxations by long-range elastically mediated interactions. By simulations and an analytical theory, we show that the model reproduces the main salient facts associated with dynamic heterogeneity and offers a mechanism to explain the emergence of dynamical correlations at the glass transition. We also discuss how it can be generalized and combined with current theories.

Composites Based on Polylactide Doped with Amorphous Europium(III) Complex as Perspective Thermosensitive Luminescent Materials

This work reports fabrication of polylactide (PLA) films doped with various additives of an amorphous Eu(III) complex. We study the temperature behavior of the luminescence intensity and lifetime of the PLA-Eu(III) composites in the range of 298–353 K and investigate the mechanism of luminescence temperature quenching. The peak relative sensitivity of the films reaches 20.1 %×K−1 and exceeds the respective characteristics of all known lanthanide-containing thermosensors designed for the range of physiological temperatures. The produced films can be potential novel materials for luminescent thermosensors.

Structural Measures as Guides to Ultrastable States in Overjammed Packings

Jammed, disordered packings of given sets of particles possess a multitude of equilibrium states with different mechanical properties. Identifying and constructing desired states, e.g., of superior stability, is a complex task. Here, we show that in two-dimensional particle packings the energy of all metastable states (inherent structures) is reliably classified by simple scalar measures of local steric packing. These structural measures are insensitive to the particle interaction potential and so robust that they can be used to guide a modified swap algorithm that anneals polydisperse packings toward low-energy metastable states exceptionally fast. The low-energy states are extraordinarily stable against applied shear, so that the approach also efficiently identifies ultrastable packings.

Microscopic structural origin behind slowing down of colloidal phase separation approaching gelation

The gelation of colloidal particles interacting through a short-range attraction is widely recognized as a consequence of the dynamic arrest of phase separation into colloid-rich and solvent-rich phases. However, the microscopic origin behind the slowing down and dynamic arrest of phase separation remains elusive. In order to access microscopic structural changes through the entire process of gelation in a continuous fashion, we used core–shell fluorescent colloidal particles, laser scanning confocal microscopy, and a unique experimental protocol that allows us to initiate phase separation instantaneously and gently. Combining these enables us to track the trajectories of individual particles seamlessly during the whole phase-separation process from the early stage to the late arresting stage. We reveal that the enhancement of local packing and the resulting formation of locally stable rigid structures slow down the phase-separation process and arrest it to form a gel with an average coordination number of z = 6–7. This result supports a mechanical perspective on the dynamic arrest of sticky-sphere systems based on the microstructure, replacing conventional explanations based on the macroscopic vitrification of the colloid-rich phase. Our findings illuminate the microscopic mechanisms behind the dynamic arrest of colloidal phase separation, the emergence of mechanical rigidity, and the stability of colloidal gels.

Stress hyperuniformity and transient oscillatory-exponential correlation decay as signatures of strength vs fragility in glasses

We examine and compare the local stress autocorrelation in the inherent states of a fragile and a strong glass: the Kob-Andersen (KA) binary mixture and the Beest-Kramer-Santen model of silica. For both systems, local (domain-averaged) stress fluctuations asymptotically reach the normal inverse-volume decay in the large domain limit; accordingly, the real-space stress autocorrelation presents long-range power law tails. However, in the case of silica, local stress fluctuations display a high degree of hyperuniformity, i.e., their asymptotic (normal) decay is disproportionately smaller than their bond level amplitude. This property causes the asymptotic power law tails of the real-space stress autocorrelation to be swamped, up to very large distances (several nanometers), by an intermediate oscillatory-exponential decay regime. Similar contributions exist in the KA stress autocorrelation, but they never can be considered as dominating the power law decay and fully disappear when stress is coarse-grained beyond one interatomic distance. Our observations document that the relevance of power-law stress correlation may constitute a key discriminating feature between strong and fragile glasses. Meanwhile, they highlight that the notion of local stress in atomistic systems involves by necessity a choice of observation (coarse-graining) scale, the relevant value of which depends, in principle, on both the model and the phenomenon studied.

Unifying Description of the Vibrational Anomalies of Amorphous Materials

The vibrational density of states D(ω) of solids controls their thermal and transport properties. In crystals, the low-frequency modes are extended phonons distributed in frequency according to Debye's law, D(ω)∝ω^{2}. In amorphous solids, phonons are damped, and at low frequency D(ω) comprises extended modes in excess over Debye's prediction, leading to the so-called boson peak in D(ω)/ω^{2} at ω_{bp}, and quasilocalized ones. Here we show that boson peak and phonon attenuation in the Rayleigh scattering regime are related, as suggested by correlated fluctuating elasticity theory, and that amorphous materials can be described as homogeneous isotropic elastic media punctuated by quasilocalized modes acting as elastic heterogeneities. Our numerical results resolve the conflict between theoretical approaches attributing amorphous solids' vibrational anomalies to elastic disorder and localized defects.

Towards Glasses with Permanent Stability

Unlike crystals, glasses age or devitrify over time, reflecting their nonequilibrium nature. This lack of stability is a serious issue in many industrial applications. Here, we show by numerical simulations that the devitrification of quasi-hard-sphere glasses is prevented by suppressing volume-fraction inhomogeneities. A monodisperse glass known to devitrify with "avalanchelike" intermittent dynamics is subjected to small iterative adjustments to particle sizes to make the local volume fractions spatially uniform. We find that this entirely prevents structural relaxation and devitrification over aging time scales, even in the presence of crystallites. There is a dramatic homogenization in the number of load-bearing nearest neighbors each particle has, indicating that ultrastable glasses may be formed via "mechanical homogenization." Our finding provides a physical principle for glass stabilization and opens a novel route to the formation of mechanically stabilized glasses.

Percolation transitions in compressed SiO2 glasses

Amorphous-amorphous transformations under pressure are generally explained by changes in the local structure from low- to higher-fold coordinated polyhedra^1-4. However, as the notion of scale invariance at the critical thresholds has not been addressed, it is still unclear whether these transformations behave similarly to true phase transitions in related crystals and liquids. Here we report ab initio-based calculations of compressed silica (SiO₂) glasses, showing that the structural changes from low- to high-density amorphous structures occur through a sequence of percolation transitions. When the pressure is increased to 82 GPa, a series of long-range ('infinite') percolating clusters composed of corner- or edge-shared tetrahedra, pentahedra and eventually octahedra emerge at critical pressures and replace the previous 'phase' of lower-fold coordinated polyhedra and lower connectivity. This mechanism provides a natural explanation for the well-known mechanical anomaly around 3 GPa, as well as the structural irreversibility beyond 10 GPa, among other features. Some of the amorphous structures that have been discovered mimic those of coesite IV and V crystals reported recently^5,6, highlighting the major role of SiO₅ pentahedron-based polyamorphs in the densification process of vitreous silica. Our results demonstrate that percolation theory provides a robust framework to understand the nature and pathway of amorphous-amorphous transformations and open a new avenue to predict unravelled amorphous solid states and related liquid phases^7,8.

Emergence of linear isotropic elasticity in amorphous and polycrystalline materials

We investigate the emergence of isotropic linear elasticity in amorphous and polycrystalline solids via extensive numerical simulations. We show that the elastic properties are correlated over a finite length scale ξ_{E}, so that the central limit theorem dictates the emergence of continuum linear isotropic elasticity on increasing the specimen size. The stiffness matrix of systems of finite size L>ξ_{E} is obtained, adding to that predicted by linear isotropic elasticity a random one of spectral norm (L/ξ_{E})^{-3/2} in three spatial dimensions. We further demonstrate that the elastic length scale corresponds to that of structural correlations, which in polycrystals reflect the typical size of the grain boundaries and length scales characterizing correlations in the stress field. We finally demonstrate that the elastic length scale affects the decay of the anisotropic long-range correlations of locally defined shear modulus and shear stress.

An energy-landscape-based crossover temperature in glass-forming liquids

The systematic identification of temperature scales in supercooled liquids that are key to understanding those liquids' underlying glass properties, and their formation-history dependence, is a challenging task. Here, we study the statistics of particles' squared displacements δr² between equilibrium liquid configurations at temperature T and their underlying inherent states, using computer simulations of 11 different computer glass formers. We show that the relative fluctuations of δr² are nonmonotonic in T, exhibiting a maximum whose location defines the crossover temperature T_X. Therefore, T_X marks the point of maximal heterogeneity during the process of tumbling down the energy landscape, starting from an equilibrium liquid state at temperature T down to its underlying inherent state. We extract T_X for the 11 employed computer glasses, ranging from tetrahedral glasses to packs of soft elastic spheres, and demonstrate its usefulness in putting the elastic properties of different glasses on the same footing. Interestingly, we further show that T_X marks the crossover between two distinct regimes of the mean ⟨δr²⟩: a high temperature regime in which ⟨δr²⟩ scales approximately as T^0.5 and a deeply supercooled regime in which ⟨δr²⟩ scales approximately as T^1.3. Further research directions are discussed.