Identification and characterization of potential shear transformation zones in metallic glasses

Topological characterization of rearrangements in amorphous solids

In amorphous materials, plasticity is localized and occurs as shear transformations. It was recently shown by Wu et al. that these shear transformations can be predicted by applying topological defect concepts developed for liquid crystals to an analysis of vibrational eigenmodes [Z. W. Wu et al., Nat. Commun. 14, 2955 (2023)10.1038/s41467-023-38547-w]. This study relates the -1 topological defects to the displacement fields expected of an Eshelby inclusion, which are characterized by an orientation and the magnitude of the eigenstrain. A corresponding orientation and magnitude can be defined for each defect using the local displacement field around each defect. These parameters characterize the plastic stress relaxation associated with the local structural rearrangement and can be extracted using the fit to either the global displacement field or the local field. Both methods provide a reasonable estimation of the molecular-dynamics-measured stress drop, confirming the localized nature of the displacements that control both long-range deformation and stress relaxation.

Using delaunay triangularization to characterize non-affine displacement fields during athermal, quasistatic deformation of amorphous solids

We investigate the non-affine displacement fields that occur in two-dimensional Lennard-Jones models of metallic glasses subjected to athermal, quasistatic simple shear (AQS). During AQS, the shear stress versus strain displays continuous quasi-elastic segments punctuated by rapid drops in shear stress, which correspond to atomic rearrangement events. We capture all information concerning the atomic motion during the quasi-elastic segments and shear stress drops by performing Delaunay triangularizations and tracking the deformation gradient tensor F_α associated with each triangle α. To understand the spatio-temporal evolution of the displacement fields during shear stress drops, we calculate F_α along minimal energy paths from the mechanically stable configuration immediately before to that after the stress drop. We find that quadrupolar displacement fields form and dissipate both during the quasi-elastic segments and shear stress drops. We then perform local perturbations (rotation, dilation, simple and pure shear) to single triangles and measure the resulting displacement fields. We find that local pure shear deformations of single triangles give rise to mostly quadrupolar displacement fields, and thus pure shear strain is the primary type of local strain that is activated by bulk, athermal quasistatic simple shear. Other local perturbations, e.g. rotations, dilations, and simple shear of single triangles, give rise to vortex-like and dipolar displacement fields that are not frequently activated by bulk AQS. These results provide fundamental insights into the non-affine atomic motion that occurs in driven, glassy materials.

The role of configurational disorder on plastic and dynamic deformation in Cu64Zr36 metallic glasses: A molecular dynamics analysis

The varying degrees of configurational disorder in metallic glasses are investigated quantitatively by molecular dynamics studies. A parameter, the quasi-nearest atom, is used to characterize the configurational disorder in metallic glasses. Our observations suggest configurational disorder play a role in structural heterogeneity, plasticity and dynamic relaxations in metallic glasses. The broad configurational disorder regions distribution is the indicator of abundant potential deformation units and relaxations. Plastic flow, as well as relaxation, is believed to start at configurational disorder regions. The width of the shear bands and dynamic relaxations can then be regulated by the degree of configurational disorder regions in metallic glasses.

Unusual fast secondary relaxation in metallic glass

The relaxation spectrum of glassy solids has long been used to probe their dynamic structural features and the fundamental deformation mechanisms. Structurally complicated glasses, such as molecular glasses, often exhibit multiple relaxation processes. By comparison, metallic glasses have a simple atomic structure with dense atomic packing, and their relaxation spectra were commonly found to be simpler than those of molecular glasses. Here we show the compelling evidence obtained across a wide range of temperatures and frequencies from a La-based metallic glass, which clearly shows two peaks of secondary relaxations (fast versus slow) in addition to the primary relaxation peak. The discovery of the unusual fast secondary relaxation unveils the complicated relaxation dynamics in metallic glasses and, more importantly, provides us the clues which help decode the structural features serving as the ‘trigger' of inelasticity on mechanical agitations.

Mechanical relaxation processes in glasses can provide information on the structural and mechanical properties of glasses. Here, the authors observe a fast secondary relaxation process in La-based metallic glasses, providing information on the inelasticity of metallic glasses.

Crossover from Localized to Cascade Relaxations in Metallic Glasses

Thermally activated deformation is investigated in two metallic glass systems with different cooling histories. By probing the atomic displacements and stress changes on the potential energy landscape, two deformation modes, a localized process and cascade process, have observed. The localized deformation involves fewer than 30 atoms and appears in both systems, and its size is invariant with cooling history. However, the cascade deformation is more frequently observed in the fast quenched system than in the slowly quenched system. The origin of the cascade process in the fast quenched system is attributed to the higher density of local minima on the underlying potential energy landscape.

Rheological properties vs. local dynamics in model disordered materials at low temperature

We study the rheological response at low temperature of a sheared model disordered material as a function of the bond rigidity. We find that the flow curves follow a Herschel-Bulkley law, whatever is the bond rigidity, with an exponent close to 0.5. Interestingly, the apparent viscosity can be related to a single relevant time scale t rel, suggesting a strong connection between the local dynamics and the global mechanical behaviour. We propose a model based on the competition between the nucleation and the avalanche-like propagation of spatial strain heterogeneities. This model can explain the Herschel-Bulkley exponent on the basis of the size dependence of the heterogeneities on the shear rate.

Heterogeneous relaxation dynamics in amorphous materials under cyclic loading

Molecular dynamics simulations are performed to investigate heterogeneous dynamics in amorphous glassy materials under oscillatory shear strain. We consider three-dimensional binary Lennard-Jones mixture well below the glass transition temperature. The structural relaxation and dynamical heterogeneity are quantified by means of the self-overlap order parameter and the dynamic susceptibility. We found that at sufficiently small strain amplitudes, the mean square displacement exhibits a broad subdiffusive plateau and the system undergoes nearly reversible deformation over about 10(4) cycles. Upon increasing strain amplitude, the transition to the diffusive regime occurs at shorter time intervals and the relaxation process involves intermittent bursts of large particle displacements. The detailed analysis of particle hopping dynamics and the dynamic susceptibility indicates that mobile particles aggregate into clusters whose sizes increase at larger strain amplitudes. Finally, the correlation between particle mobilities in consecutive time intervals demonstrates that dynamic facilitation becomes increasingly pronounced at larger strain amplitudes.

Shear-induced phase transition of nanocrystalline hexagonal boron nitride to wurtzitic structure at room temperature and lower pressure

Disordered structures of boron nitride (BN), graphite, boron carbide (BC), and boron carbon nitride (BCN) systems are considered important precursor materials for synthesis of superhard phases in these systems. However, phase transformation of such materials can be achieved only at extreme pressure-temperature conditions, which is irrelevant to industrial applications. Here, the phase transition from disordered nanocrystalline hexagonal (h)BN to superhard wurtzitic (w)BN was found at room temperature under a pressure of 6.7 GPa after applying large plastic shear in a rotational diamond anvil cell (RDAC) monitored by in situ synchrotron X-ray diffraction (XRD) measurements. However, under hydrostatic compression to 52.8 GPa, the same hBN sample did not transform to wBN but probably underwent a reversible transformation to a high-pressure disordered phase with closed-packed buckled layers. The current phase-transition pressure is the lowest among all reported direct-phase transitions from hBN to wBN at room temperature. Usually, large plastic straining leads to disordering and amorphization; here, in contrast, highly disordered hBN transformed to crystalline wBN. The mechanisms of strain-induced phase transformation and the reasons for such a low transformation pressure are discussed. Our results demonstrate a potential of low pressure-room temperature synthesis of superhard materials under plastic shear from disordered or amorphous precursors. They also open a pathway of phase transformation of nanocrystalline materials and materials with disordered and amorphous structures under extensive shear.

Plasticity and dynamical heterogeneity in driven glassy materials

Many amorphous glassy materials exhibit complex spatio-temporal mechanical response and rheology, characterized by an intermittent stress strain response and a fluctuating velocity profile. Under quasistatic and athermal deformation protocols this heterogeneous plastic flow was shown to be composed of plastic events of various sizes, ranging from local quadrupolar plastic rearrangements to system spanning shear bands. In this paper, through numerical study of a 2D Lennard-Jones amorphous solid, we generalize the study of the heterogeneous dynamics of glassy materials to the finite shear rate (gamma not equal to 0) and temperature case (T not equal to 0). In practice, we choose an effectively athermal limit (T approximately 0) and focus on the influence of shear rate on the rheology of the glass. In line with previous works we find that the model Lennard-Jones glass follows the rheological behavior of a yield stress fluid with a Herschel-Bulkley response of the form, sigma = sigmaY + c1gamma(beta). The global mechanical response obtained through the use of Molecular Dynamics is shown to converge in the limit gamma --> 0 to the quasistatic limit obtained with an energy minimization protocol. The detailed analysis of the plastic deformation at different shear rates shows that the glass follows different flow regimes. At sufficiently low shear rates the mechanical response reaches a shear-rate-independent regime that exhibits all the characteristics of the quasistatic response (finite-size effects, cascades of plastic rearrangements, yield stress, ...). At intermediate shear rates the rheological properties are determined by the externally applied shear rate and the response deviates from the quasistatic limit. Finally at higher shear the system reaches a shear-rate-independent homogeneous regime. The existence of these three regimes is also confirmed by the detailed analysis of the atomic motion. The computation of the four-point correlation function shows that the transition from the shear-rate-dominated to the quasistatic regime is accompanied by the growth of a dynamical cooperativity length scale xi that is shown to diverge with shear rate as xi is proportional to gamma(-nu), with nu approximately 0.2 -0.3. This scaling is compared with the prediction of a simple model that assumes the diffusive propagation of plastic events.

Local elasticity map and plasticity in a model Lennard-Jones glass

In this work we calculate the local elastic moduli in a weakly polydispersed two-dimensional Lennard-Jones glass undergoing a quasistatic shear deformation at zero temperature. The numerical method uses coarse-grained microscopic expressions for the strain, displacement, and stress fields. This method allows us to calculate the local elasticity tensor and to quantify the deviation from linear elasticity (local Hooke's law) at different coarse-graining scales. From the results a clear picture emerges of an amorphous material with strongly spatially heterogeneous elastic moduli that simultaneously satisfies Hooke's law at scales larger than a characteristic length scale of the order of five interatomic distances. At this scale, the glass appears as a composite material composed of a rigid scaffolding and of soft zones. Only recently calculated in nonhomogeneous materials, the local elastic structure plays a crucial role in the elastoplastic response of the amorphous material. For a small macroscopic shear strain, the structures associated with the nonaffine displacement field appear directly related to the spatial structure of the elastic moduli. Moreover, for a larger macroscopic shear strain we show that zones of low shear modulus concentrate most of the strain in the form of plastic rearrangements. The spatiotemporal evolution of this local elasticity map and its connection with long term dynamical heterogeneity as well as with the plasticity in the material is quantified. The possibility to use this local parameter as a predictor of subsequent local plastic activity is also discussed.

Distribution of thermally activated plastic events in a flowing glass

The potential energy landscape of a flowing metallic glass is revealed using the activation-relaxation technique. For a two-dimensional Lennard-Jones system initially deformed into a steady-state condition through quasistatic shear, the distribution of activation energies is shown to contain a large fraction of low-energy barriers, consistent with a highly nonequilibrium flow state. The distribution of plastic strains has a fundamentally different shape than that obtained during quasistatic simulations, exhibiting a peak at finite strain and, after elastic unloading, a nonzero mean plastic strain that evidences a polarization of the flow state. No significant correlation is found between the activation energy of a plastic event and its associated plastic strain.