Soft modes and nonaffine rearrangements in the inherent structures of supercooled liquids

Local Plastic Response and Slow Heterogeneous Dynamics of Supercooled Liquids

We demonstrate, via numerical simulations, that the relaxation dynamics of supercooled liquids correlates well with a plastic length scale measuring a particle's response to impulsive localized perturbations and weakly to measures of local elasticity. We find that the particle averaged plastic length scale vanishes linearly in temperature and controls the super-Arrhenius temperature dependence of the relaxation time. Furthermore, we show that the plastic length scale of individual particles correlates with their typical displacement at the relaxation time. In contrast, the local elastic response only correlates with the dynamics on the vibrational timescale.

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.

Atomic vibration as an indicator of the propensity for configurational rearrangements in metallic glasses

In a metallic glass (MG), the propensity for atomic rearrangements varies spatially from location to location in the amorphous solid, making the prediction of their likelihood a major challenge. One can attack this problem from the "structure controls properties" standpoint. But all the current structure-centric parameters are mostly based on local atomic packing information limited to short-range order, hence falling short in reliably forecasting how the local region would respond to external stimuli (e.g., temperature and/or stress). Alternatively, one can use indicators informed by physical properties to bridge the static structure on the one hand, and the response of the local configuration on the other. A sub-group of such physics-informed quantities consists of atomic vibration parameters, which will be singled out as the focus of this article. Here we use the Cu₆₄Zr₃₆ alloy to systematically demonstrate the following two points, all using a single model MG. First, we show in a comprehensive manner the interrelation among common vibrational parameters characterizing the atomic vibrational amplitude and frequency, including the atomic mean square displacement, flexibility volume, participation fraction in the low-frequency vibrational modes and boson peak intensity. Second, we demonstrate that these vibrational parameters fare much better than purely static structural parameters based on local geometrical packing in providing correlation with the propensity for local configurational transitions. These vibrational parameters also share a correlation length similar to that in structural rearrangements induced by external stimuli. This success, however, also poses a challenge, as it remains to be elucidated as to why short-time dynamical (vibrational) behavior at the bottom of the energy basin can be exploited to project the height of the energy barrier for cross-basin activities and in turn the propensity for locally collective atomic rearrangements.

Structural and short-time vibrational properties of colloidal glasses and supercooled liquids in the vicinity of the re-entrant glass transition

We investigate the short-time vibrational properties and structure of two-dimensional, bidisperse, colloidal glasses and supercooled liquids in the vicinity of the re-entrant glass transition, as a function of interparticle depletion attraction strength. The long-time spatiotemporal dynamics of the samples are measured to be non-monotonic, confirming that the suspensions evolve from repulsive glass to supercooled liquid to attractive glass with increasing depletion attraction. Here, we search for vibrational signatures of the re-entrant behavior in the short-time spatiotemporal dynamics, i.e., dynamics associated with particle motion inside its nearest-neighbor cage. Interestingly, we observe that the anharmonicity of these in-cage vibrations varies non-monotonically with increasing attraction strength, consistent with the non-monotonic long-time structural relaxation dynamics of the re-entrant glass. We also extract effective spring constants between neighboring particles; we find that spring stiffness involving small particles also varies non-monotonically with increasing attraction strength, while stiffness between large particles increases monotonically. Last, from study of depletion-dependent local structure and vibration participation fractions, we gain microscopic insight into the particle-size-dependent contributions to short-time vibrational modes in the glass and supercooled liquid states.

Residual stress distributions in amorphous solids from atomistic simulations

The distribution of local residual stresses (threshold to instability) that controls the statistical properties of plastic flow in athermal amorphous solids is examined with an atomistic simulation technique. For quiescent configurations, the distribution has a pseudogap (power-law) form with an exponent that agrees well with global yielding statistics. As soon as deformation sets in, the pseudogap region gives way to a system size dependent plateau at small residual stresses that can be understood from the statistics of local residual stress differences between plastic events. Results further suggest that the local yield stress in amorphous solids changes even if the given region does not participate in plastic activity.

Computational study of transient shear banding in soft jammed solids

We have designed three-dimensional numerical simulations of a soft spheres model, with size polidispersity and in athermal conditions, to study the transient shear banding that occurs during yielding of jammed soft solids. We analyze the effects of different types of drag coefficients used in the simulations and compare the results obtained using Lees-Edwards periodic boundary conditions with the case in which the same model solid is confined between two walls. The specific damping mechanism and the different boundary conditions indeed modify the load curves and the velocity profiles in the transient regime. Nevertheless, we find that the presence of a stress overshoot and of a related transient banding phenomenon, for large enough samples, is a robust feature for overdamped systems, where their presence do not depend on the specific drag used and on the different boundary conditions.

Degeneracy in molecular scale organization of biological membranes

The scale-rich spatiotemporal organization in biological membranes has its origin in the differential inter- and intra-molecular interactions among their constituents. In this work, we explore the molecular-origin behind that variety and possible degeneracy in lateral organization in membranes. For our study, we post-process microsecond long all-atom molecular dynamics trajectories for three systems that exhibit fluid phase coexistence: (i) PSM/POPC/Chol (0.47/0.32/0.21), (ii) PSM/DOPC/Chol (0.43/0.38/0.19) and (iii) DPPC/DOPC/Chol (0.37/0.36/0.27). To distinguish the liquid ordered and disordered regions at molecular scales, we calculate the degree of non-affineness of individual lipids in their neighbourhood and track their topological rearrangements. Disconnectivity graph analysis with respect to membrane organization shows that the DPPC/DOPC/Chol and PSM/DOPC/Chol systems exhibit funnel-like energy landscapes as opposed to a highly frustrated energy landscape for the more biomimetic PSM/POPC/Chol system. We use these measurements to develop a continuous lattice Hamiltonian and evolve that using Monte Carlo simulated annealing to explore the possibility of structural degeneracy in membrane organization. Our data show that model membranes with lipid constituents that are biomimetic (PSM/POPC/Chol) have the ability to access a large range of membrane sub-structure space (higher degeneracy) as compared to the other two systems, which form only one kind of substructure even with changing composition. Since the spatiotemporal organization in biological membranes dictates the "molecular encounters" and in turn larger scale biological processes such as molecular transport, trafficking and cellular signalling, we posit that this structural degeneracy could enable access to a larger repository to functionally important molecular organization in systems with physiologically relevant compositions.

Vibrational Modes and Dynamic Heterogeneity in a Near-Equilibrium 2D Glass of Colloidal Kites

Using video microscopy and particle-tracking techniques developed for dense Brownian systems of polygons, we study the structure-dynamics relationship in a near-equilibrium 2D glass consisting of anisotropic Penrose kite-shaped colloids. Detailed vibrational properties of kite glasses, both translational and rotational, are obtained using covariance matrix techniques. Different from other colloidal glasses of spheres and ellipsoids, the vibrational modes of kite glasses at low frequencies show a strong translational character with spatially localized rotational modes and extended translational modes. Low-frequency quasilocalized soft modes commonly found in sphere glasses are absent in the translational phonon modes of kite glasses. Soft modes are observed predominantly in the rotational vibrations and correlate well with the spatial distribution of Debye-Waller factors. The local structural entropy field shows a strong correlation with the observed dynamic heterogeneity.

Structural and topological nature of plasticity in sheared granular materials

Upon mechanical loading, granular materials yield and undergo plastic deformation. The nature of plastic deformation is essential for the development of the macroscopic constitutive models and the understanding of shear band formation. However, we still do not fully understand the microscopic nature of plastic deformation in disordered granular materials. Here we used synchrotron X-ray tomography technique to track the structural evolutions of three-dimensional granular materials under shear. We establish that highly distorted coplanar tetrahedra are the structural defects responsible for microscopic plasticity in disordered granular packings. The elementary plastic events occur through flip events which correspond to a neighbor switching process among these coplanar tetrahedra (or equivalently as the rotation motion of 4-ring disclinations). These events are discrete in space and possess specific orientations with the principal stress direction.

It is a general consensus that the structural defects are the plasticity carriers in amorphous solids, but its microscopic view remains largely unknown. Cao et a. show that highly distorted coplanar tetrahedra act as defects in granular packings, which flip under shear to carry local plasticity.

Yielding of a model glass former: An interpretation with an effective system of icosahedra

We consider the yielding under simple shear of a binary Lennard-Jones glass former whose super-Arrhenius dynamics are correlated with the formation of icosahedral structures. We recast this glass former as an effective system of icosahedra [Pinney et al., J. Chem. Phys. 143, 244507 (2015)JCPSA60021-960610.1063/1.4938424]. Looking at the small-strain region of sheared simulations, we observe that shear rates affect the shear localization behavior particularly at temperatures below the glass transition as defined with a fit to the Vogel-Fulcher-Tamman equation. At higher temperature, shear localization starts immediately on shearing for all shear rates. At lower temperatures, faster shear rates can result in a delayed start in shear localization, which begins close to the yield stress. Building from a previous work which considered steady-state shear [Pinney et al., J. Chem. Phys. 143, 244507 (2015)JCPSA60021-960610.1063/1.4938424], we interpret the response to shear and the shear localization in terms of a local effective temperature with our system of icosahedra. We find that the effective temperatures of the regions undergoing shear localization increase significantly with increasing strain (before reaching a steady-state plateau).

Predicting Shear Transformation Events in Metallic Glasses

Shear transformation is the elementary process for plastic deformation of metallic glasses, the prediction of the occurrence of the shear transformation events is therefore of vital importance to understand the mechanical behavior of metallic glasses. In this Letter, from the view of the potential energy landscape, we find that the protocol-dependent behavior of shear transformation is governed by the stress gradient along its minimum energy path and we propose a framework as well as an atomistic approach to predict the triggering strains, locations, and structural transformations of the shear transformation events under different shear protocols in metallic glasses. Verification with a model Cu_{64}Zr_{36} metallic glass reveals that the prediction agrees well with athermal quasistatic shear simulations. The proposed framework is believed to provide an important tool for developing a quantitative understanding of the deformation processes that control mechanical behavior of metallic glasses.

Local yield stress statistics in model amorphous solids

We develop and extend a method presented by Patinet, Vandembroucq, and Falk [Phys. Rev. Lett. 117, 045501 (2016)PRLTAO0031-900710.1103/PhysRevLett.117.045501] to compute the local yield stresses at the atomic scale in model two-dimensional Lennard-Jones glasses produced via differing quench protocols. This technique allows us to sample the plastic rearrangements in a nonperturbative manner for different loading directions on a well-controlled length scale. Plastic activity upon shearing correlates strongly with the locations of low yield stresses in the quenched states. This correlation is higher in more structurally relaxed systems. The distribution of local yield stresses is also shown to strongly depend on the quench protocol: the more relaxed the glass, the higher the local plastic thresholds. Analysis of the magnitude of local plastic relaxations reveals that stress drops follow exponential distributions, justifying the hypothesis of an average characteristic amplitude often conjectured in mesoscopic or continuum models. The amplitude of the local plastic rearrangements increases on average with the yield stress, regardless of the system preparation. The local yield stress varies with the shear orientation tested and strongly correlates with the plastic rearrangement locations when the system is sheared correspondingly. It is thus argued that plastic rearrangements are the consequence of shear transformation zones encoded in the glass structure that possess weak slip planes along different orientations. Finally, we justify the length scale employed in this work and extract the yield threshold statistics as a function of the size of the probing zones. This method makes it possible to derive physically grounded models of plasticity for amorphous materials by directly revealing the relevant details of the shear transformation zones that mediate this process.

Local thermal energy as a structural indicator in glasses

Identifying heterogeneous structures in glasses-such as localized soft spots-and understanding structure-dynamics relations in these systems remain major scientific challenges. Here, we derive an exact expression for the local thermal energy of interacting particles (the mean local potential energy change caused by thermal fluctuations) in glassy systems by a systematic low-temperature expansion. We show that the local thermal energy can attain anomalously large values, inversely related to the degree of softness of localized structures in a glass, determined by a coupling between internal stresses-an intrinsic signature of glassy frustration-anharmonicity and low-frequency vibrational modes. These anomalously large values follow a fat-tailed distribution, with a universal exponent related to the recently observed universal [Formula: see text] density of states of quasilocalized low-frequency vibrational modes. When the spatial thermal energy field-a "softness field"-is considered, this power law tail manifests itself by highly localized spots, which are significantly softer than their surroundings. These soft spots are shown to be susceptible to plastic rearrangements under external driving forces, having predictive powers that surpass those of the normal modes-based approach. These results offer a general, system/model-independent, physical/observable-based approach to identify structural properties of quiescent glasses and relate them to glassy dynamics.

Cage Size and Jump Precursors in Glass-Forming Liquids: Experiment and Simulations

Glassy dynamics is intermittent, as particles suddenly jump out of the cage formed by their neighbors, and heterogeneous, as these jumps are not uniformly distributed across the system. Relating these features of the dynamics to the diverse local environments explored by the particles is essential to rationalize the relaxation process. Here we investigate this issue characterizing the local environment of a particle with the amplitude of its short time vibrational motion, as determined by segmenting in cages and jumps the particle trajectories. Both simulations of supercooled liquids and experiments on colloidal suspensions show that particles in large cages are likely to jump after a small time-lag, and that, on average, the cage enlarges shortly before the particle jumps. At large time-lags, the cage has essentially a constant size, which is smaller for longer-lasting cages. Finally, we clarify how this coupling between cage size and duration controls the average behavior and opens the way to a better understanding of the relaxation process in glass-forming liquids.

Structure in sheared supercooled liquids: Dynamical rearrangements of an effective system of icosahedra

We consider a binary Lennard-Jones glassformer whose super-Arrhenius dynamics are correlated with the formation of particles organized into icosahedra under simple steady state shear. We recast this glassformer as an effective system of icosahedra [Pinney et al., J. Chem. Phys. 143, 244507 (2015)]. From the observed population of icosahedra in each steady state, we obtain an effective temperature which is linearly dependent on the shear rate in the range considered. Upon shear banding, the system separates into a region of high shear rate and a region of low shear rate. The effective temperatures obtained in each case show that the low shear regions correspond to a significantly lower temperature than the high shear regions. Taking a weighted average of the effective temperature of these regions (weight determined by region size) yields an estimate of the effective temperature which compares well with an effective temperature based on the global mesocluster population of the whole system.

Universal structural parameter to quantitatively predict metallic glass properties

Quantitatively correlating the amorphous structure in metallic glasses (MGs) with their physical properties has been a long-sought goal. Here we introduce ‘flexibility volume' as a universal indicator, to bridge the structural state the MG is in with its properties, on both atomic and macroscopic levels. The flexibility volume combines static atomic volume with dynamics information via atomic vibrations that probe local configurational space and interaction between neighbouring atoms. We demonstrate that flexibility volume is a physically appropriate parameter that can quantitatively predict the shear modulus, which is at the heart of many key properties of MGs. Moreover, the new parameter correlates strongly with atomic packing topology, and also with the activation energy for thermally activated relaxation and the propensity for stress-driven shear transformations. These correlations are expected to be robust across a very wide range of MG compositions, processing conditions and length scales.

Various known structural descriptors of metallic glasses have limitations in quantitatively predicting properties. Here authors define a physically-motivated measure and show it to correlate strongly with elastic properties, local structure and relaxation kinetics over a wide range of simulated compositions.

Connecting Local Yield Stresses with Plastic Activity in Amorphous Solids

In model amorphous solids produced via differing quench protocols, a strong correlation is established between local yield stress measured by direct local probing of shear stress thresholds and the plastic rearrangements observed during remote loading in shear. This purely local measure shows a higher predictive power for identifying sites of plastic activity when compared with more conventional structural properties. Most importantly, the sites of low local yield stress, thus defined, are shown to be persistent, remaining predictive of deformation events even after fifty or more such plastic rearrangements. This direct and nonperturbative approach gives access to relevant transition pathways that control the stability of amorphous solids. Our results reinforce the relevance of modeling plasticity in amorphous solids based on a gradually evolving population of discrete and local zones preexisting in the structure.

Particle jumps in structural glasses

Particles in structural glasses rattle around temporary equilibrium positions, that seldom change through a process which is much faster than the relaxation time, known as particle jump. Since the relaxation of the system is due to the accumulation of many such jumps, it could be possible to connect the single particle short time motion to the macroscopic relaxation by understanding the features of the jump dynamics. Here we review recent results in this research direction, clarifying the features of particle jumps that have been understood and those that are still under investigation, and examining the role of particle jumps in different theories of the glass transition.

Nonlinear plastic modes in disordered solids

We propose a theoretical framework within which a robust micromechanical definition of precursors to plastic instabilities, often termed soft spots, naturally emerges. They are shown to be collective displacements (modes) z[over ̂] that correspond to local minima of a barrier function b(z[over ̂]), which depends solely on inherent structure information. We demonstrate how some heuristic searches for local minima of b(z[over ̂]) can a priori detect the locus and geometry of imminent plastic instabilities with remarkable accuracy, at strains as large as γ_{c}-γ∼10^{-2} away from the instability strain γ_{c}. Our findings suggest that the a priori detection of the entire field of soft spots can be effectively carried out by a systematic investigation of the landscape of b(z[over ̂]).

Correlation between rearrangements and soft modes in polymer glasses during deformation and recovery

We explore the link between soft vibrational modes and local relaxation events in polymer glasses during physical aging, active deformation at constant strain rate, and subsequent recovery. A softness field is constructed out of the superposition of the amplitudes of the lowest energy normal modes, and found to predict up to 70% of the rearrangements. Overlap between softness and rearrangements increases logarithmically during aging and recovery phases as energy barriers rise due to physical aging, while yielding rapidly rejuvenates the overlap to that of a freshly prepared glass. In the strain hardening regime, correlations rise for uniaxial tensile deformation but not for simple shear. These trends can be explained by considering the differing degrees of localization of the soft modes in the two deformation protocols.

Spatial correlations of elementary relaxation events in glass-forming liquids

The dynamical facilitation scenario, by which localized relaxation events promote nearby relaxation events in an avalanche process, has been suggested as the key mechanism connecting the microscopic and the macroscopic dynamics of structural glasses. Here we investigate the statistical features of this process via numerical simulations of a model structural glass. First we show that the relaxation dynamics of the system occurs through particle jumps that are irreversible, and that cannot be decomposed in smaller irreversible events. Then we show that each jump does actually trigger an avalanche. The characteristics of this avalanche change upon cooling, suggesting that the relaxation dynamics crossovers from a noise dominated regime, where jumps do not trigger other relaxation events, to a regime dominated by the facilitation process, where a jump triggers more relaxation events.

Memory effects in schematic models of glasses subjected to oscillatory deformation

We consider two schematic models of glasses subjected to oscillatory shear deformation, motivated by the observations, in computer simulations of a model glass, of a nonequilibrium transition from a localized to a diffusive regime as the shear amplitude is increased, and of persistent memory effects in the localized regime. The first of these schematic models is the NK model, a spin model with disordered multi-spin interactions previously studied as a model for sheared amorphous solids. The second model, a transition matrix model, is an abstract formulation of the manner in which occupancy of local energy minima evolves under oscillatory deformation cycles. In both of these models, we find a behavior similar to that of an atomic model glass studied earlier. We discuss possible further extensions of the approaches outlined.

Structural relaxation is a scale-free process

We show that in deeply supercooled liquids, structural relaxation proceeds via the accumulation of Eshelby events, i.e. local rearrangements that create long-ranged and anisotropic stresses in the surrounding medium. Such events must be characterized using tensorial observables and we construct an analytical framework to probe their correlations using local stress data. By analyzing numerical simulations, we then demonstrate that events are power-law correlated in space, with a time-dependent amplitude which peaks at the alpha relaxation time τ(α). This effect, which becomes stronger near the glass transition, results from the increasingly important role of local stress fluctuations in facilitating relaxation events. Our finding precludes the existence of any length scale beyond which the relaxation process decorrelates.

Structural relaxation in glassy polymers predicted by soft modes: a quantitative analysis

We present a quantitative analysis of the correlation between quasi-localized, low energy vibrational modes and structural relaxation events in computer simulations of a quiescent, thermal polymer glass. Our results extend previous studies on glass forming binary mixtures in 2D, and show that the soft modes identify regions that undergo irreversible rearrangements with up to 7 times the average probability. We study systems in the supercooled- and aging-regimes and discuss temperature- as well as age-dependence of the correlation. In addition to the location of rearrangements, we find that soft modes also predict their direction on the molecular level. The soft regions are long lived structural features, and the observed correlations vanish only after >50% of the system has undergone rearrangements.

Evolution of atomic rearrangements in deformation in metallic glasses

Atomic rearrangements induced by shear stress are fundamental for understanding deformation mechanisms in metallic glasses (MGs). Using molecular dynamic simulation, the atomic rearrangements characterized by nonaffine displacements (NADs) and their spatial distribution and evolution with tensile stress in Cu50Zr50 MG were investigated. It was found that in the elastic regime the atomic rearrangements with the largest NADs are relatively homogeneous in space, but exhibit strong spatial correlation, become localized and inhomogeneous, and form large clusters as strain increases, which may facilitate the so-called shear transformation zones. Furthermore, initially they prefer to take place around Cu atoms which have more nonicosahedral configurations. As strain increases, the preference decays and disappears in the plastic regime. The atomic rearrangements with the smallest NADs are preferentially located around Cu atoms, too, but with more icosahedral or icosahedral-like atomic configurations. The preference is maintained in the whole deformation process. In contrast, the atomic rearrangements with moderate NADs distribute homogeneously, and do not show explicit preference or spatial correlation, acting as matrix during deformation. Among the atomic rearrangements with different NADs, those with largest and smallest NADs are nearest neighbors initially, but separating with increasing strain, while those with largest and moderate NADs always avoid to each other. The correlations in the fluctuations of the NADs confirm the long-range strain correlation and the scale-free characteristic of NADs in both elastic and plastic deformation, which suggests a universality of the scaling in the plastic flow in MGs.

Soft spots and their structural signature in a metallic glass

In a 3D model mimicking realistic Cu64Zr36 metallic glass, we uncovered a direct link between the quasi-localized low-frequency vibrational modes and the local atomic packing structure. We also demonstrate that quasi-localized soft modes correlate strongly with fertile sites for shear transformations: geometrically unfavored motifs constitute the most flexible local environments that encourage soft modes and high propensity for shear transformations, whereas local configurations preferred in this alloy, i.e., the full icosahedra (around Cu) and Z16 Kasper polyhedra (around Zr), contribute the least.