Tittaya Thaiyanurak

Mean atomic bond strain (MABS), based on the globally averaged bond length, has recently emerged as a new strain metric that retains clear physical meaning even as severe atomic neighborhood reconstruction occurs. It has been shown to exhibit a nearly perfect linear relationship with global stress throughout the elastic and plastic deformation in single-crystal face-centered cubic (FCC) metals, contradicting conventional expectations based on nonlinear dislocation activity. Whether this near-linear relationship holds in other materials stands out as an important and intriguing question. In this study, we examine the MABS–stress relationship in representative unary, binary, and ternary metallic glasses (MGs), where neither a crystal structure nor dislocations are present. Large-scale molecular dynamics simulations of uniaxial tensile tests and statistical analysis of millions of atomic bonds are performed. Irrespective of their differing compositions, all the MGs exhibit a persistent near-linear relationship between total MABS (all bonds included) and global stress up to fracture, even in the presence of significant local nonaffine displacements (shear transformation zones and shear bands), with the Pearson correlation coefficient consistently exceeding 0.99. Unlike the nonaffine displacements, the spatial distribution of individual atomic bond strain does not localize under the uniaxial loading. In the MGs containing more than one element, MABS computed for a single bond type may not correlate as linearly with global stress as total MABS. The results demonstrate that the persistent near-linear total MABS–stress relationship over the entire deformation process, recently discovered in single-crystal FCC metals, also applies to MGs despite their vastly different atomic structures. This strengthens the candidacy of total MABS as a universal stress descriptor across materials classes and deformation regimes. With further development and implementation in atomistic simulations and constitutive modeling, the MABS concept has the potential to reshape our understanding of materials mechanics and generate new insights into the design of stronger, tougher, and more thermally and chemically stable materials.

Atomic-scale strain is the basis of a material’s macroscopic deformation behavior. The current measure of atomic-scale strain in the form of the Green–Lagrange tensor loses its physical meaning beyond the yield point, as atomic neighborhoods undergo significant reconstructions. We have recently introduced a new atomic-scale strain measure, namely, atomic bond strain, through our study of bond behavior in multicomponent metallic glasses. Here, we apply this new strain measure to uniaxial tensile tests (simulated using molecular dynamics) of several representative single-crystal FCC (face-centered cubic) metals under varied strain rates. We show that this new strain measure displays remarkable near-linear correlation with stress, not only in the elastic regime, but also in the plastic regime where complex dislocation dynamics (nucleation, bursting, motion, annihilation, regeneration) and stress fluctuations take place. This suggests that the overall stress of the materials even in the plastic regime is predominantly determined by the degree of bond stretching among all atoms. This appears to contradict the common conceptions that the plastic flow stress of a crystalline material is governed by dislocation events involving only a small fraction of atoms around dislocations, and that the stress–strain relationship is highly non-linear for plastic deformation. The contradictions can be reconciled by considering the causal sequence: dislocation events alter bond stretching, and bond stretching directly determines the stress. This brings a novel insight into the nature of plastic deformation, owing to the newly introduced atomic bond strain. How well the near-linear correlation between the stress and the atomic bond strain holds in other materials (e.g., non-FCC single crystals, polycrystals, quasicrystals, elements, alloys, and compounds) is an intriguing and important topic for future investigation, following the example of this work.

The mechanical properties of materials are fundamentally determined by the behavior of atomic bonds under stress. Probing bond behavior during deformation, however, is highly challenging, particularly for materials with complex chemical compositions and/or atomic structures, such as metallic glasses (MGs). As a result, a significant gap exists in the current understanding of the mechanical properties of MGs in relation to the atomic bond behavior and how this relationship is influenced by metallurgical factors (e.g., alloy composition, processing conditions). Here, we present our study of the compositional effects on the tensile behavior of atomic bonds in Cu_93−xZr_xAl₇ (x = 40, 50, 60 at.%) MGs using large-scale molecular dynamics (MD) simulations and statistical analysis. Specifically, we examine the populations (fractions), mean bond lengths, mean bond z-lengths, and mean bond z-strains of the different bond types before and during tensile loading (in the z-direction), and we compare these quantities across the different alloy compositions. Among our key findings, we show that increasing the Zr content in the alloy composition leads to shortened Zr-Zr, Al-Cu, Al-Zr, and Cu-Zr bonds and elongated Cu-Cu bonds, as evidenced by their mean bond lengths. During deformation, the shorter Zr-Zr bonds and longer Cu-Cu bonds in the higher-Zr-content alloys, compared with those in the x = 40 alloy, appear stronger (more elastic stretching in the z-direction) and weaker (less z-stretching), respectively, consistent with general expectations. In contrast, the Al-Cu, Al-Zr, and Cu-Zr bonds in the higher-Zr-content alloys appear weaker in the elastic regime, despite their shortened mean bond lengths. This apparent paradox can be reconciled by considering the fractions of these bonds associated with icosahedral clusters, which are known to be more resistant to deformation than the rest of the glassy structure. We also discuss how the compositional effects on the bond behavior relate to variations in the overall stress–strain behavior of the different alloys.

Multicomponent metallic glasses (MGs) are a fascinating class of advanced alloys known for their exceptional properties such as limit-approaching strength, high hardness and corrosion resistance, and near-net-shape castability. One important question regarding these materials that remains unanswered is how the different elements and atomic bonds within them control their strength and deformability. Here, we present a detailed visual and statistical analysis of the behaviors of various elements and atomic bonds in the Zr₄₇Cu₄₆Al₇ (at%) MG during a uniaxial tensile test (in the z-direction) simulated using molecular dynamics. Specifically, we investigate the identities of atoms undergoing significant shear strain, and the averaged bond lengths, projected z-lengths, and z-angles (angles with respect to the z-direction) of all the atomic bonds as functions of increasing strain. We show that, prior to yielding, the Zr element and the intermediate (Zr-Zr, Cu-Al) and stronger (Zr-Al, Zr-Cu) bonds dominate the elastic deformation and strength, while the Cu and Al elements and the weaker Al-Al and Cu-Cu bonds contribute more to the highly localized shear transformation. The significant reconstruction, as signified by the cessation of bond-length increment and bond-angle decrement, of the intermediate and the stronger bonds triggers yielding of the material. After yielding, all the elements and bonds participate in the plastic deformation while the stronger bonds contribute more to the residual strength and the ultimate (fracture) strain. The results provide new insights into the atomic mechanisms underlying the mechanical behavior of multicomponent MGs, and may assist in the future design of MG compositions towards better combination of strength and deformability.

It is generally known that the incorporation of crystals in the glass matrix can enhance the ductility of metallic glasses (MGs), at the expense of reduced strength, and that the deformation of MGs, particularly during shear banding, can induce crystal formation/growth. Here, we show that these known trends for the interplay between crystals and deformation of MGs may hold true or become inverted depending on the size of the crystals relative to the shear bands. We performed molecular dynamics simulations of tensile tests on nanocrystal-bearing MGs. When the crystals are relatively small, they bolster the strength rather than the ductility of MGs, and the crystals within a shear band undergo redissolution as the shear band propagates. In contrast, larger crystals tend to enhance ductility at the cost of strength, and the crystal volume fraction increases during deformation. These insights offer a more comprehensive understanding of the intricate relationship between deformation and crystals/crystallization in MGs, useful for fine-tuning the structure and mechanical properties of both MGs and MG–crystal composites.

The recently discovered Cu₄₆Zr_33.5Hf_13.5Al₇ (at.%) bulk metallic glass (BMG) presents the highest glass-forming ability (GFA) among all known copper-based alloys, with a record-breaking critical casting thickness (or diameter) of 28.5 mm. At present, much remains to be explored about this new BMG that holds exceptional promise for engineering applications. Here, we report our study on the crystallization behavior of this new BMG, using isochronal and isothermal differential scanning calorimetry (DSC), X-ray diffraction (XRD), and transmission electron microscopy (TEM). With the calorimetric data, we determine the apparent activation energy of crystallization, the Avrami exponent, and the lower branch of the isothermal time–temperature–transformation (TTT) diagram. With XRD and TEM, we identify primary and secondary crystal phases utilizing samples crystallized to different degrees within the calorimeter. We also estimate the number density, nucleation rate, and growth rate of the primary crystals through TEM image analysis. Our results reveal that the crystallization in this BMG has a high activation energy of ≈360 kJ/mole and that the primary crystallization of this BMG produces a high number density (≈10²¹ m⁻³ at 475 °C) of slowly growing (growth rate < 0.5 nm/s at 475 °C) Cu₁₀(Zr,Hf)₇ nanocrystals dispersed in the glassy matrix, while the second crystallization event further produces a new phase, Cu(Zr,Hf)₂. The results help us to understand the GFA and thermal stability of this new BMG and provide important guidance for its future engineering applications, including its usage as a precursor to glass–crystal composite or bulk nanocrystalline structures.

Understanding crystallization mechanisms in nano-sized metallic glasses (MGs) is important to the manufacturing and application of these new nanomaterials that possess a unique combination of structural and functional properties. Due to the two-dimensional projections and limited spatial and/or temporal resolutions in experiments, significant questions (e.g., whether nucleation takes place on the free surface or in a near-surface layer) regarding this subject remain under debate. Here, we address these outstanding questions using molecular dynamics simulations of crystallization in MG nanorods together with atomistic visualization and data analysis. We show that nucleation in the nano-sized MGs predominantly takes place on the surface by converting the high-energy liquid surface to a lower-energy crystal surface (the most close-packed atomic plane). This is true for all the nanorods with different diameters studied. On the other hand, the apparent growth mode (inward/radial, lateral or longitudinal) and the resulting grain structure are more dependent on the nanorod diameter. For a relatively big diameter of the nanorod, the overall growth rate does not differ much among the three directions and the resulting grains are approximately semispherical. For small diameters, grains appear to grow more in longitudinal direction and some grains may form relatively long single-crystal segments along the length of the nanorod. The reasons for the difference are discussed. The study provides direct atomistic insights into the crystallization mechanisms in nano-sized MGs, which can facilitate the manufacturing and application of these new advanced materials.