## Brittleness of metallic glasses dictated by their state at the fragile-to-strong transition temperature

Achraf Atila, 1, ∗ Sergey V. Sukhomlinov, 1 Marc J. Honecker, 1 and Martin H. M¨ser u 1

1

Department of Material Science and Engineering, Saarland University, Saarbr¨cken, u 66123, Germany

(Dated: August 2, 2024)

The effect of cooling on the brittleness of glasses in general, and bulk metallic glasses (BMGs) in particular, is usually studied with continuously varying cooling rates; slower cooling rates lead to stiffer, harder, and more brittle glasses than higher cooling rates. These protocols obscure any potential discontinuity that a glass might experience depending on whether its microstructure resembles that of a fragile or a strong glass-forming liquid. Here, we use large-scale molecular dynamics to simulate the nanoindentation behavior of model BMGs (Zr 0 6 . Cu 0 3 . Al 0 1 . ) obtained by rapidly quenching equilibrium melts from temperatures above and below the fragile-to-strong transition temperature T fst , leading to fragile and strong glasses, respectively. While the contact modulus deduced from the indentation simulation evolves smoothly with the temperature T q from which the melt is quenched, the plastic response changes quasi-discontinuously as T q passes through T fst . In particular, strong glasses develop highly asymmetric flow profiles with mature shear bands, in contrast to fragile glasses. Quantitative differences reveal themselves not only through a formal von Mises localization parameter analysis but also through image analysis of flow patterns using pretrained artificial intelligence models. Moreover, seemingly erratic flow profiles for our indentation geometry produced surprisingly reproducible and, thus, deterministic features. It remains to be determined to what extent other classes of glass formers follow our observation that the degree of brittleness is significantly influenced by whether the melt is fragile or strong when it falls out of equilibrium at the glass transition temperature.

Keywords: Bulk metallic glasses, Deformation behavior, Shear bands, fragile-to-strong transition, Atomistic simulations

## I. INTRODUCTION

Bulk metallic glasses (BMGs) have elastic strain limits of up to 5% [1] and high tensile [2] as well as yield [3] strengths of order 1 GPa exceeding that of crystalline and poly-crystalline metals by an order of magnitude. At the same time, BMGs have defied the rule that materials cannot be simultaneously strong and tough [3]. In fact, BMGs are known to have the highest known damage tolerance [3], which is the product of fracture yield strength and toughness. On the downside, BMGs generally strain to soften, making them prone to localized deformations in the form of shear bands [4-6], whereby undesirable surface markings can be produced, which ultimately evolve into cracks. This contrasts conventional metals, which usually strain hardens as the dislocations are the main carriers of plasticity [7, 8]. This behavior limits the moldability of most BMGs and their widespread technological applicability.

quickly. In contrast, silica-based oxide glasses are always brittle. A distinguishing feature of silica-based glasses is the ratio of the fragile-to-strong transition (FST) temperature T fst and the glass transition temperature T g , i.e., typical values of T fst /T g clearly exceed unity for silicates but are of order unity for BMGs. This prompts us to hypothesize that the brittleness of a glass is significantly influenced by the degree to which its atomic structure resembles that of a strong or fragile melt. In particular, the question arises whether the quasi-discontinuity of both dynamic and static properties of the melt at the FST translates into discontinuities in the mechanical properties of the corresponding glass.

One subject of ongoing research is how the thermal history of a BMG affects its atomic structure and, thereby, its deformation mechanisms, most notably shear bands [4, 9, 10]. It is established both experimentally [11] and from simulations [12, 13] that (bulk-metallic) glasses produced with a small cooling rate are stronger but also more prone to shear banding than those quenched more

∗ achraf.atila@uni-saarland.de

A few notes on the FST, which has been argued to occur in any glass-forming melt [14], and the glass transition itself may be in place. The (equilibrium) viscosity and atomic diffusion coefficients cross over from a nonArrhenius to an Arrhenius dependence upon cooling in a rather narrow temperature window near T fst , which is accompanied by a breakdown of the Stokes-Einstein relation below T fst [15, 16]. The specific heat of melts c p at T &gt; T fst is substantially above that associated with the rule of Dulong-Petit. This means that structural rearrangements go beyond (quasi-) harmonic vibrations so that the melt is called (thermodynamically) fragile . After passing through a maximum near T fst upon cooling, c p quickly decreases to values satisfying Dulong-Petit. This is indicative of vibrations around barely changing reference structures. The melt is called strong . The potentially most localized abrupt change at the FST is the

jump of the correlation length describing the asymptotic decay of density oscillation [17, 18]. In contrast to the FST, the glass transition refers to a more continuous process in which a melt falls out of equilibrium during cooling, however, with a strong dependence on the cooling rate. Similar to the FST, the glass transition leads to a reduction of the specific heat toward the rule of DulongPetit.

Ascertaining to what degree FST discontinuities translate into corresponding features in the glass requires equilibrium melts to be quenched instantly from both above and below T fst to well below T g in the first step. This way, the glass represents the frozen-in structure of a fragile and strong melt, respectively, and hence will be called either a fragile glass or a strong glass in the following. Whenever T g and T fst are not clearly separated, conventional cooling protocols using constant cooling rates plus potentially holding periods below T g generally produce a mix of fragile and strong glass, since the melt would not be fully equilibrated at T &gt; T fst , while some lowtemperature relaxations would occur at T &lt; T fst . The lower the cooling rate ˙, q the more the structure leans toward a strong glass, but changes are continuous in ˙. q Thus, previous studies [13, 19, 20] revealing that the failure mechanisms evolves smoothly from being homogeneous in a high- ˙ sample to localized shear bands using q low ˙ q do not contradict our hypothesis that the brittleness of a (bulk metallic) glass is dictated by its state at the FST.

In the remainder of this work, we use molecular dynamics (MD) simulations to contrast the mechanical behavior of strong and fragile glasses during nanoindentation. To this end we recycle previously produced atomistic equilibrium configurations of a ( in-silico model for a) ternary BMG, Zr 0 6 . Cu 0 3 . Al 0 1 . , which had been shown to undergo a FST in thermal equilibrium [18]. MD is an ideal tool for this study because equilibrium melts can be quenched extremely quickly and homogeneously, preventing temperature-gradient-induced heterogeneities or residual stresses from occurring.

## II. METHODS

## A. Molecular dynamics

In this work, we simulate the ternary alloy Zr 0 6 . Cu 0 3 . Al 0 1 . using an embedded-atom potential [21] parameterized by Cheng, Ma, and Sheng [22]. This system revealed a transition from a fragile to a strong liquid at T fst ≈ 830 K [18], at which the correlation length of the long-range density fluctuations changed quasidiscontinuously by roughly 10% in a very narrow temperature range similarly as in a related experimental BMG called Vitleroy 106a [17].

Initial configurations were taken from the abovementioned previous work [18], in which a sample containing ≈ 8800 atoms was equilibrated at various tempera- tures, including a narrow temperature regime below T fst . The full equilibration times used during step-wise cooling were: 30 ns at T = 1 000 K, 64 ns at 900 K, 400 ns at , 850 K, 1.25 µ s at 812 K, 2.5 µ s at 800 K, and 4.0 µ s at 787 K. In the following, samples obtained by quenching an equilibrium sample from T &lt; T fst will be referred to as strong glasses, while those obtained from equilibrium simulations at T &gt; T fst will be called fragile glasses. All samples were duplicated n x = 3 , n y = 3 , n z = 3 times in x , y , and z -direction, respectively, equilibrated for 4 ns at the temperature from which they were obtained, and subsequently quenched to room temperature (T = 300 K with cooling rate of 1 K / ps, and annealed for another 4 ns. During these simulations, a timestep of 2 fs was used with periodic boundary conditions applied to all box directions. Such equilibrated samples were replicated again to obtain a larger sample with 230 x 16.5 x 86 nm 3 , used in the nanoindentation simulations. Using these large samples ( ≈ 18 Mio atoms), a surface was created along the z -axis, and a 1 nm thick layer of atoms at the bottom of the simulation box along the z-direction was fixed. The samples were periodic along the x -and y -direction. Next, the samples were further equilibrated for 200 ps at room temperature and 0 MPa using a timestep of 2 fs. The temperature was maintained using a Langevin thermostat [23] while the pressure was controlled using Nos´ e-Hoover chain barostat [24]. All simulations were conducted using the LAMMPS simulation package [25] and the snapshots were prepared using OVITO [26].

The indentation was performed with a rigid cylindrical indenter exerting a force of magnitude F r ( ) = -K r ( -R ) 2 on each atom, where K is the force constant which we chose to be K = 10 eV/ ˚ , A 3 r is the distance from the center of the indenter to the respective atom, and R is the radius of the indenter which was fixed to R = 100 nm. Lilleodden et al. [27] studied the response of single-crystalline gold to nanoindentation while varying the force constant of the indenter and found that the force-displacement curves are unaffected by the change of K . Displacement-controlled nanoindentation simulations were performed by displacing the indenter along the z -direction towards the sample with a fixed velocity of v i = 5m s until reaching a / depth of 20 nm and using a timestep of 2 fs.

## B. Local strain analysis

The detection of shear bands was performed by using the von Mises shear strain. The von Mises strain was calculated by getting an affine transformation matrix to the relative neighbors through a least-square fitting [28]. The undeformed glass was chosen as a reference configuration. More details of the strain calculation can be found in Refs. [26, 28]. To get the regions with a nonaffine displacement, the current atomic configuration needs to be compared with a reference configuration. The neighbors of each particle are found by defining a sampling radius

(i.e., 10 ˚ in our case). A The local strain is then calculated from the displacements of the neighboring particles to the central one and the relative displacements that they would have if they were in a region of uniform strain. The final value of the nonaffine displacement was divided by the number of atoms within that spherical region of the same radius. More details about the calculation of the nonaffine displacements are given in Ref. [29].

## C. Image analysis

Snapshots of the local von Mises shear strain taken at different depths during the nanoindentation process were image analyzed. To this end, pre-trained DINOv2 neural networks [30] were used for feature extraction because DINOv2 is known for its robustness in capturing image features through self-supervised learning. The DINOv2 architecture is based on the so-called vision transformer (ViT) [31], which is currently the de-facto standard for computer vision. In this work, we mainly use the bare backbones ViT-b/14 with registers [32] and ViT-g/14 without registers available at GitHub ( https://github.com/facebookresearch/dinov2 ), because they showed the best results for our analysis. The models had been pre-trained on 142 million images in a contrastive learning setting, which allows the model to capture robust image features, although images supposedly did not show shear bands or flow lines.

Each image was preprocessed to meet the input requirements of the DINOv2 model. The preprocessing steps included resizing the image to 224 × 224 pixels, converting it to a PyTorch [33] tensor, and applying the DINOv2-specific preprocessing function. The latter includes normalizing the images by the mean and standard deviation values of ImageNet [34] statistics: mean = [0.485, 0.456, 0.406] and standard deviation = [0.229, 0.224, 0.225]. After putting all images through the model, the cosine similarity of each pair of images was computed. It ranges from -1 (anti-parallel, minimum similarity) to 1 (parallel, maximum similarity).

## III. RESULTS

In this work, we simulate the indentation process of glasses obtained from quenching equilibrium melts from six temperatures T q to T = 300 K, where all indentation simulations were conducted. Three of the initial temperatures, that is, T 1 = 787 K, T 2 = 800 K, and T 3 = 812 K were below the fragile-to-strong crossover temperature T fst = 830 K, the remaining ones T 4 = 850 K, T 5 = 900 K, and T 6 = 1 000 K were above , T fst . To be concise, the sample quenched from T 1 may not have fully reached equilibrium at T 1 . To be specific, the relaxation time of energy fluctuations in the α -relaxation regime can be estimated to be τ = 10 ns and the Kohlrausch exponent to β = 0 73, which yields .

a 1 -exp( -( t sim /τ ) β ) = 0 995 complete . α relaxation given an simulation time of t sim = 100 ns. For the volume calculation, the relaxation is even more complete. In this case, τ = 1 ns and β = 0 79 yielding . 1 -exp( -( t sim /τ ) β ) = 1.

At first glance, all six samples show qualitatively similar force-displacement curves as revealed in Fig. 1: an initial Hertzian regime up to an indentation of d ≈ 2 nm, followed by a non-elastic regime. A maximum indentation force of F ≈ 7 ± 0 5 . µ N is reached at d ≈ 10 ± 1 nm for all T q ≤ 900 K samples. After a short drop, F starts increasing again at d ≈ 12 ± 2 nm.

FIG. 1. Force displacement curves of loading for the fragile glasses and the strong glass at T = 300K. The dashed line reflects an indentation with a contact modulus of E ∗ = 57 9 GPa. Deviations of the data from linearity at very small . d are due to an effectively finite-range repulsion between the indenter and samples.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00536v1-figure-43948de0-e748-441c-b5c3-1fb2eb8ab30d.jpg)

Closer inspection of Fig. 1 reveals qualitative differences between fragile and strong glasses. All strong samples have similar F N ( d ) relations at d &lt; d max , where d max is the depth at which the force is a local maximum. In contrast, F N ( d ) of fragile samples reveal some T q dependence at d &lt; d max . The situation reverts when d exceeds d max by a few nanometers. This time, F N ( d ) converges to a single dependence for all fragile glasses but reveals a weak T q dependence for the strong glasses. Moreover, the drop in the indentation force after reaching its maximum is noticeably more pronounced for strong than for fragile glasses.

Subtle qualitative differences between the fragile and the strong glasses can already be noticed at a relatively small indentation depth. The strong glass starts to deviate from their linear response by assuming a left-curved F N ( d ) relation, while the fragile glasses adopt a minor right-curved F N ( d ) dependence. In contrast, the contact modulus E ∗ ( T q ) has a smooth dependence on T q without detectable signs of a discontinuity of the function itself or its derivative, as is revealed in Fig. 2 These findings align

FIG. 2. Contact modulus E ∗ at T = 300 K as a function of T q as deduced from the initial (0 2 . &lt; d/ nm &lt; 2 5) linear . F N ( d ) dependence of the indentation process using F N ( d ) ≈ π 4 E Ld ∗ . The very initial d &lt; 0 2 nm domain is excluded from fit be-. cause the finite stiffness of the cylinder-BMG repulsion causes a parabolic dependence at these extremely small indentation depths.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00536v1-figure-f25666d6-afa3-4519-b380-54bd38eefc80.jpg)

with the existing literature, which also reports more significant effects of cooling rate on plastic than on elastic properties [13].

Despite subtle differences between fragile and strong plasticity at the early loading stages, effects are more pronounced at large indentation. To analyze plasticity, we investigated the local von Mises local shear strain, as described in Sec. II. Exemplary results are shown in Fig. 3 with representatives corresponding to the three different plasticity regimes of the F d ( ) curve, i.e., before, near, and after the maximum indentation force. Fig. S1 depicts more snapshots of the von Mises strain, i.e., for all six samples at seven depths, including one snapshot taken in the elastic regime. We also note that visualization of the shear bands using non-affine displacement squared, defined similarly as by Falk and Langer [29] but normalized to the number of neighbors, produced snapshots, which closely resembled those based on the von Mises strain, as is evidenced in Fig. S2.

Shear bands, which are the main carrier of plasticity in metallic glasses [6, 35-37] start to occur once F d ( ) deviates from the elastic response. They mature, i.e., the shear bands reach a stage where their growth slows down or stops after F d ( ) reaches its maximum, at which point they have hit the surface. Contrary to previous simulations [4, 13, 37-39] and our fragile glasses, the strong samples develop visible, radially shaped shear bands with an average width of w ≈ 8 nm. This is on par with w ≈ 10 nm [36, 40, 41], which is the intrinsic thickness of mature shear bands in ex-silico metallic glasses. The shear bands of strong glasses are more localized and more asymmetric than in fragile glasses, to the extent that they could be said to differ qualitatively.

Going beyond a visual inspection of the shear bands requires a quantitative analysis to be conducted. To this end, we study the standard deviation of the atomic von

Mises strain divided by the average atomic von Mises strain, Ψ = 〈 ∆ η 2 n 〉 1 / 2 / η 〈 n 〉 , where η n is the von Mises strain projected to atom n , 〈 ... 〉 an average over all atoms, and ∆ η n = η n -〈 η n 〉 . Ψ is similar to a previously defined localization parameter [42] but normalized so that it does not change when atoms with zero von Mises strain are added. Like the previous parameter, Ψ is sensitive to localization because atoms with large variance tend to cluster in shear bands. Results for Ψ are shown in Fig. 4. At very small d , the value for Ψ is essentially identical for all samples. This can be understood from the fact that elastic deformation is scale-invariant and, thereby, the value for Ψ at d = 0 5 nm is merely 2% between the . most fragile and the strongest sample compared to a 25% difference in E ∗ .

The values of Ψ( d ) bifurcate at a depth of d /greaterorsimilar 4 nm, where Ψ( d ) keeps increasing for the strong samples while they saturate otherwise. An even stronger qualitative difference of strong and fragile samples becomes apparent at d /greaterorsimilar 9 nm, where the slope of Ψ( d ) suddenly increases near the point of the local F N ( d ) minima for the strong samples, while Ψ( d ) has no such discontinuity in the fragile samples.

From the results shown so far, a qualitative difference in the deformation behavior of glasses reflecting the structure of a fragile and a strong melt has become obvious. To ensure that these differences are not merely a kinetic effect resulting from the ratio of the indentation rate and the relaxation time at T q surpassing some critical value, we performed another indentation simulation of the T q = 850 K sample, however, using a ten-fold increased indentation velocity. This exceeds the ratio of relaxations times between the 'hottest strong' and the 'coldest fragile' melt, for both energy and volume relaxation.

The ten-fold increased indentation velocity produced a noticeable increase in the F N ( d ) relationship. However, the higher indentation velocity produces noticeably more symmetric and less localized shear bands or rather flow patterns, as is revealed in Fig. 5 both visually but also in terms of the strain-localization parameter. The higher symmetry of the flow images at the higher velocity could be argued to arise from the reduced time for an initial local symmetry breaking to migrate from small scales to large scales. In our simulations, the symmetry breaking arises due to structural heterogeneities, potentially also, to some degree, by thermal fluctuations. In case of instability, which is unavoidable in strain-softening materials, numerical round-off errors, even in a perfect symmetric finite-element simulation, might suffice to cause similar asymmetries in the shear bands as in the MD simulations when the loading is slow.

To explore the degree to which our findings can be tested experimentally if images of flow profiles are available but there is no spatially resolved von Mises strain, we quantified the similarity between the shear-band signature using two different models, as described in the Methods section. The result produced with ViT g/14 is

FIG. 3. Local shear strain maps for the strong and fragile glasses loaded to an indentation depth of (a) d = 10 nm, (b) d = 12 nm, and (c) d = 15 nm.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00536v1-figure-ca8496eb-7c91-410a-875d-158603146eca.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00536v1-figure-798f2ef5-a3f2-4796-a5c7-4c205e65d5c6.jpg)

FIG. 4. Strain localization parameter Ψ as a function of the indentation depth for the strong and fragile glasses. Different symbols are used for their respective temperatures T q : circles for 787 K, diamonds for 800 K, stars for 812 K, triangles down for 850 K, pluses for 900 K, and crosses for 1000 K.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00536v1-figure-1931fe8e-128e-4dbb-bbb9-d02464902c5f.jpg)

shown in Fig. 6. At each indentation depth, inspection of both colors and numbers included in the figure generally gives the largest difference in similarity between the hottest strong and the coolest fragile samples. This claim holds for both models. Moreover, the ViT g/14 judges the similarity between the upper left and upper right panel in Fig. 3 to be more than that between the top

FIG. 5. Contrasting the indentation process of the T q = 850 K fragile glass using two different indentation velocities. Panel (a) shows the normal force F N and (b) the strain-localization parameter Ψ as a function of displacement d for both velocities. Panel (c) and (d) depict the von Mises shear strain for v 0 = 5 m/s and v 0 = 50 m/s, respectively.

right and the bottom right, i.e., 0.94 and 0.88 cosine similarity, respectively. This is an interesting finding given that the covered colored area is similarly large for the two samples indented to 15 nm.

Figuratively speaking, the human eye can clearly identify the top left panel as the child and the upper right as

its adult, rugose version, while the upper and lower right panels represent two different adults.

Since neither the color maps nor the numbers included in Fig. 6 allow us to conclude whether the two models reveal discontinuities between fragile and strong glasses, a normalized similarity S n ( T q ) of a sample quenched from arbitrary T q and indented to d = 15 nm with that quenched from T q = 1 000 K is shown in Fig. 7. , It is defined as the normalized cosine similarity between the sample quenched from T q = 1000 K and all other samples at an indentation depth of 15 nm, as defined by

$$S _ { n } ( T ) = \frac { S _ { \cos } ^ { 1 0 0 0 K } ( T ) - \min ( S _ { \cos } ^ { 1 0 0 0 K } ( T _ { q } ) ) } { 1 - \min ( S _ { \cos } ^ { 1 0 0 0 K } ( T _ { q } ) ) }, \quad ( 1 ) \quad \ t r a n \colon \, \ s i m u$$

where S T 1 cos ( T 2 ) is the cosine similarity of the flow profiles of quench temperature T 1 and T 2 . This measure makes the two studied models operate on a similar domain, because it is normalized such that the similarity between the least similar and the T q = 1000 K sample is zero.

Although the similarity indices show a clustering for the strong samples and a jump between strong and fragile samples, the judgment appears less clear than that by humans or the von Mises localization parameter. It is nevertheless remarkable that both models produce similar results, although the ViT-b/14 has less than 8% of the parameters compared to ViT-g/14.

A final question to be addressed here is to what extent the shape of the flow profiles and shear bands, particularly those obtained for the strong glasses, is erratic or deterministic. All strong samples develop a long shear band to the right with spider-web-like features between the shear band and the surface, plus a shorter but thicker shear band to the left, which is very close to the surface. The probability for breaking the symmetry when quenching from three independent samples each time in the same direction is 1/4, which is small compared to unity but not very small so that it is unclear whether some hidden memory in the strong melt causes the symmetry to be broken each time in the same direction. We remind the reader that all final configurations share the same ancestor, e.g., the sample obtained after equilibrating the melt at T = 812 K was quenched to yield the T q = 812 K glass at 300 K, but also copied and further equilibrated at 800 K for times much exceeding the energy and volume relaxation time.

Since equilibration and indentation simulations require several million CPU hours, averaging over independently produced samples is impractical. To nevertheless ascertain how deterministic features are, we repeated the simulations for a larger cell, in which the quenched samples were not only repeated 14 × 1 × 5 times in x , y , and z -direction respectively but 15 × 1 × 6 times ( ≈ 21.3 Mio). Moreover, the original simulation cell kept its orientation one time but was rotated by 90 ◦ so that the original yz surfaces became the new to-be-indented xy surface. Von Mises strains of these three different T q = 812 K samples are contrasted in Fig. 8. The two new configurations show the same characteristics as the old one, i.e., one prominent long shear band, which breaks the symmetry and has the spider-web feature underneath it. Thus, not only the gross features of the shear profile appear to be predetermined.

## IV. DISCUSSION

The FST is conjectured to exist in all liquids [43, 44]. MGs prepared experimentally are quenched with cooling rates that are fast enough to avoid crystallization and, at the same time, slow enough to allow for this FTS transition [45]. However, the glasses prepared in MD simulations are usually cooled with cooling rates several decades faster than those used in experiments [13, 20, 4648], preventing them from experiencing the FTS liquid transition [18]. The results of this work show that this renders a one-to-one comparison between simulation and experiments problematic. However, suppose glasses are prepared in a way that makes them undergo FTS liquid transition. In that case, the differences between the simulated and experimentally prepared glasses will reduce, and a one-to-one comparison of the mechanisms will be much trusted. This was shown to be the case recently when authors combined MD and Monte Carlo swaps to equilibrate the system and produce stable glasses with properties closer to the experimental ones [13].

Regarding the mechanical behavior of MGs, Packard and Schuh [49] showed that for isotropic materials such as metallic glasses, the formation of narrow shear bands in the area with concentrated shear stresses that exceed the yield stress. Nevertheless, from our results presented above, the first shear transformation zones in the bulk sample appeared at an indentation depth of about 2.5 nm and ≈ 15 nm under the indenter, but not all of them are exactly along the axis of indentation, suggesting that the local shear stresses do not uniquely determine the STZ formation. Similar observations were made by Wang et al. [50]. This can be explained by the fact that the local variation of the atomic configuration results in BMG having different zones with low or high shear stresses distributed in the sample and with STZ preferring to nucleate at the sites with low shear resistance [4, 51], which dictate the location where STZ form and also the propagation of the shear bands afterward. This is in line with the observations we presented, where asymmetric shear bands are observed in the strong glass, although the cylindrical indenter is symmetric. In addition to the strain field applied by the indenter, the shear band will naturally propagate towards the surface because it provides no constraints on their propagation [50].

The differences between strong and fragile glasses can be explained by the degree of stability of the glasses and the concept of the potential energy landscape (PEL) [52]. The strong glass was cooled using a much lower effective cooling rate ≈ 10 7 K/s, giving the atoms more time to relax and move towards deeper energy minima. On the other hand, the fragile glasses were cooled from temper-

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00536v1-figure-98b34daf-d4a5-472f-90d4-33dbf21a4103.jpg)

FIG. 6. Cosine similarity matrix using ViTg14 for the images shown in Fig. S1. The first number in each row and column gives T q in Kelvin, while the second gives d in nm.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00536v1-figure-5d38ba34-83f4-42b4-95ab-4e8f4bd0e476.jpg)

ble glasses [13, 52]. Although their focus was not related to the FST, experimentally, it has been proven that the processing history of the glass dictates the yielding and the shear band dynamics in the sample [53]. Shi and Falk [12, 54] made closely similar conclusions using binary Lennard-Jones glasses, where they observed a stronger localization of the local shear strain in the samples cooled slowly.

FIG. 7. Normalized similarities measure of a flow profile obtained from a glass quenched from various T q and from 1 000 K at d = 15 nm for the ViTb/14 reg (plus symbols) and the ViTg/14 (star symbols) models.

atures above the FTS transition temperature, where the samples did not explore deeper energy minima. Thus, those glasses are less stable than the strong ones. This implies that, when comparing the mechanical properties and the mechanical behavior of the strong and the fragile glasses, the strong glass will have larger mechanical properties and be more brittle, which is consistent with experimental and other simulation results that show that well-relaxed glasses have a rougher PEL, which results in lower atomic mobility, that makes the glasses more brittle and has a higher strength than the less sta-

The presence of ≈ 8 nm thick shear bands in the strong glass, which are well compared to those found experimentally [41], is reported here for the first time. The radial shape of the shear bands is consistent with the slip-lines that develop during an indentation using a cylindrical indenter and are also observed experimentally [55-57]. Previous MD simulations [37, 39] showed the presence of radial shear bands as well, but with a much lower thickness. These interesting differences in the deformation behavior of fragile and strong glasses shed light on the overall deformation mechanisms of MGs, where both localized and homogeneous deformation were observed for the same glass composition.

## V. CONCLUSIONS

In this paper, we used large-scale molecular dynamics simulations to explore how the fragility of the melt affects the brittleness of metallic glasses. We showed that melts

FIG. 8. Von Mises strain of (a) the default set up for T q = 812 K and a depth of d = 15 nm and (b) and (c) for a slightly larger sample at a depth of d = ?? nm. The starting configurations in (b) and (c) differed, i.e., the original glass sample in (c) was rotated by 90 ◦ before duplication so that the original yz surface rather than the xy surface was indented.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00536v1-figure-a3d346bb-059e-42d3-a2f6-748d49496040.jpg)

quenched from just below and just above the fragile-tostrong transition temperature show qualitatively different plasticity. The strong glasses exhibit asymmetric and localized shear bands. In contrast, the fragile glasses form broader and less defined shear bands. These differences between the strong and the fragile glasses were quantitatively observed using a modified strain localization parameter as well as by using pre-trained artificial intelligence models for image similarities. The distinctness between the deformation behavior and plastic flow profile of the fragile and strong glasses arises because strong glasses can relax into deeper energy minima, resulting in more stable atomic configurations enhancing their mechanical strength.

These findings highlight a potentially overlooked link between the thermodynamic state of the liquid at the FST and the deformation behavior of the resulting glass. This link can be expected to generalize to other classes of glass-forming liquids than BMGs. In fact, recent work

- [1] O. Glushko, R. Pippan, D. Sopu, C. Mitterer, and J. Eck-¸ ert, How to catch a shear band and explain plasticity of metallic glasses with continuum mechanics, Nature Communications 15 , 10.1038/s41467-024-49829-2 (2024).
- [2] A. Inoue, B. Shen, H. Koshiba, H. Kato, and A. R. Yavari, Cobalt-based bulk glassy alloy with ultrahigh strength and soft magnetic properties, Nature Materials 2 , 661 (2003).
- [3] R. O. Ritchie, The conflicts between strength and toughness, Nature Materials 10 , 817 (2011).
- [4] D. S ¸opu, STZ-vortex model: The key to understand STZ percolation and shear banding in metallic glasses, Journal of Alloys and Compounds 960 , 170585 (2023).
- [5] X. Mu, M. R. Chellali, E. Boltynjuk, D. Gunderov, R. Z. Valiev, H. Hahn, C. K¨ ubel, Y. Ivanisenko, and L. Velasco, Unveiling the local atomic arrangements in the shear band regions of metallic glass, Advanced Materials 33 , 2007267 (2021).
- [6] Z. Liu, Y. Yang, and C. Liu, Yielding and shear banding of metallic glasses, Acta Materialia 61 , 5928 (2013).
- [7] Z. Xie, D. Chauraud, A. Atila, E. Bitzek, S. KorteKerzel, and J. Gu´ enol´, e Thermally activated nature of synchro-shockley dislocations in laves phases,

suggested that the propensity for ductility is controlled by the topography of the energy landscape based on the simulation of various glass-forming systems [52]. It changes from being relatively smooth above T fst (high c p ) to deep quasi-harmonic minima below T fst ( c p in accordance with Dulong-Petit). However, it remains to be determined to what extent the brittleness of other glassforming systems changes abruptly depending on whether they are cooled from above or below T fst .

## ACKNOWLEDGMENTS

The authors gratefully acknowledge the Gauss Centre for Supercomputing e.V. (www.gauss-centre.eu) for funding this project by providing computing time through the John von Neumann Institute for Computing (NIC) on the GCS Supercomputer JUWELS at J¨ ulich Supercomputing Centre (JSC).

Scripta Materialia 235 , 115588 (2023).

- [8] Z. Xie, D. Chauraud, A. Atila, E. Bitzek, S. KorteKerzel, and J. Gu´ enol´, e Unveiling the mechanisms of motion of synchro-shockley dislocations in laves phases, Phys. Rev. Mater. 7 , 053605 (2023).
- [9] D. S ¸ opu, S. Scudino, X. Bian, C. Gammer, and J. Eckert, Atomic-scale origin of shear band multiplication in heterogeneous metallic glasses, Scripta Materialia 178 , 57-61 (2020).
- [10] S. Walley, Shear localization: A historical overview, Metallurgical and Materials Transactions A 38 , 2629-2654 (2007).
- [11] H. Yang, L. Zhang, J. Zhang, H. Tang, and S. Chen, Effect of sample size and cooling rate on the plastic deformation behavior of bulk metallic glasses: A comparative study, Journal of Non-Crystalline Solids 589 , 121643 (2022).
- [12] Y. Shi and M. L. Falk, Stress-induced structural transformation and shear banding during simulated nanoindentation of a metallic glass, Acta Materialia 55 , 4317 (2007).
- [13] Z. Zhang, J. Ding, and E. Ma, Shear transformations in metallic glasses without excessive and predefinable defects, Proceedings of the National Academy of Sciences 119 , 10.1073/pnas.2213941119 (2022).

- [14] C. Angell, Relaxation in liquids, polymers and plastic crystals - strong/fragile patterns and problems, Journal of Non-Crystalline Solids 131-133 , 13-31 (1991).
- [15] C. A. Angell, Formation of glasses from liquids and biopolymers, Science 267 , 1924 (1995).
- [16] C. A. Angell, K. L. Ngai, G. B. McKenna, P. F. McMillan, and S. W. Martin, Relaxation in glassforming liquids and amorphous solids, Journal of Applied Physics 88 , 3113-3157 (2000).
- [17] M. Stolpe, I. Jonas, S. Wei, Z. Evenson, W. Hembree, F. Yang, A. Meyer, and R. Busch, Structural changes during a liquid-liquid transition in the deeply undercooled Zr 58 5 . Cu 15 6 . Ni 12 8 . Al 10 3 . Nb 2 8 . bulk metallic glass forming melt, Phys. Rev. B 93 , 014201 (2016).
- [18] S. V. Sukhomlinov and M. H. M¨ser, Quasidiscontinuous u change of the density correlation length at the fragile-tostrong transition in a bulk-metallic-glass forming melt, Phys. Rev. Mater. 2 , 115604 (2018).
- [19] M. Ghaemi, M. Jafary-Zadeh, K. H. Khoo, and H. Gao, Chemical affinity can govern notch-tip brittle-to-ductile transition in metallic glasses, Extreme Mechanics Letters 52 , 101651 (2022).
- [20] M. Jafary-Zadeh, R. Tavakoli, J. Koh, Z. Aitken, and Y.-W. Zhang, Effect of chemical composition and affinity on the shortand medium-range order structures and mechanical properties of Zr-Ni-Al metallic glass, Journal of Non-Crystalline Solids 456 , 68 (2017).
- [21] M. S. Daw and M. I. Baskes, Embedded-atom method: Derivation and application to impurities,surfaces,and other defects in metals, Phys. Rev. B 29 , 6443 (1984).
- [22] Y. Q. Cheng, E. Ma, and H. W. Sheng, Atomic level structure in multicomponent bulk metallic glass, Phys. Rev. Lett. 102 , 245501 (2009).
- [23] T. Schneider and E. Stoll, Molecular-dynamics study of a three-dimensional one-component model for distortive phase transitions, Physical Review B 17 , 1302 (1978).
- [24] W. Shinoda, M. Shiga, and M. Mikami, Rapid estimation of elastic constants by molecular dynamics simulation under constant stress, Phys. Rev. B 69 , 134103 (2004).
- [25] A. P. Thompson, H. M. Aktulga, R. Berger, D. S. Bolintineanu, W. M. Brown, P. S. Crozier, P. J. in 't Veld, A. Kohlmeyer, S. G. Moore, T. D. Nguyen, R. Shan, M. J. Stevens, J. Tranchida, C. Trott, and S. J. Plimpton, Lammps -a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales, Computer Physics Communications 271 , 108171 (2022).
- [26] A. Stukowski, Visualization and analysis of atomistic simulation data with OVITO-the Open Visualization Tool, Modeling and Simulation in Materials Science and Engineering 18 , 10.1088/0965-0393/18/1/015012 (2010).
- [27] E. Lilleodden, J. Zimmerman, S. Foiles, and W. Nix, Atomistic simulations of elastic deformation and dislocation nucleation during nanoindentation, Journal of the Mechanics and Physics of Solids 51 ,
- R. Howes, P.-Y. Huang, S.-W. Li, I. Misra, M. Rabbat, V. Sharma, G. Synnaeve, H. Xu, H. Jegou, J. Mairal, P. Labatut, A. Joulin, and P. Bojanowski, DINOv2: Learning robust visual features without supervision, Transactions on Machine Learning Research (2024).
- [31] A. Dosovitskiy, L. Beyer, A. Kolesnikov, D. Weissenborn, X. Zhai, T. Unterthiner, M. Dehghani, M. Minderer, G. Heigold, S. Gelly, J. Uszkoreit, and N. Houlsby, An image is worth 16x16 words: Transformers for image recognition at scale, in International Conference on Learning Representations (2021).
- [32] T. Darcet, M. Oquab, J. Mairal, and P. Bojanowski, Vision transformers need registers (2024), arXiv:2309.16588 [cs.CV].
- [33] A. Paszke, S. Gross, F. Massa, A. Lerer, J. Bradbury, G. Chanan, T. Killeen, Z. Lin, N. Gimelshein, L. Antiga, A. Desmaison, A. K¨ opf, E. Yang, Z. DeVito, M. Raison, A. Tejani, S. Chilamkurthy, B. Steiner, L. Fang, J. Bai, and S. Chintala, Pytorch: An imperative style, high-performance deep learning librar (2019), arXiv:1912.01703 [cs.LG].
- [34] O. Russakovsky, J. Deng, H. Su, J. Krause, S. Satheesh, S. Ma, Z. Huang, A. Karpathy, A. Khosla, M. Bernstein, A. C. Berg, and L. Fei-Fei, Imagenet large scale visual recognition challenge (2015), arXiv:1409.0575 [cs.CV].
- [35] M. Hassani, A. E. Lagogianni, and F. Varnik, Probing the degree of heterogeneity within a shear band of a model glass, Phys. Rev. Lett. 123 , 195502 (2019).
- [36] A. Greer, Y. Cheng, and E. Ma, Shear bands in metallic glasses, Materials Science and Engineering: R: Reports 74 , 71 (2013).
- [37] Y. Yang, J. Luo, L. Huang, G. Hu, K. D. Vargheese, Y. Shi, and J. C. Mauro, Crack initiation in metallic glasses under nanoindentation, Acta Materialia 115 , 413 (2016).
- [38] O. Adjaoud and K. Albe, Nanoindentation of nanoglasses tested by molecular dynamics simulations: Influence of structural relaxation and chemical segregation on the mechanical response, Frontiers in Materials 8 , 10.3389/fmats.2021.664220 (2021).
- [39] D. Zhao, B. Zhu, S. Wang, Y. Niu, L. Xu, and H. Zhao, Effects of pre-strain on the nanoindentation behaviors of metallic glass studied by molecular dynamics simulations, Computational Materials Science 186 , 110073 (2021).
- [40] H. Sheng, D. Sopu, S. Fellner, J. Eckert, and C. Gammer, ¸ Mapping shear bands in metallic glasses: From atomic structure to bulk dynamics, Physical Review Letters 128 , 10.1103/physrevlett.128.245501 (2022).
- [41] Y. Zhang and A. L. Greer, Thickness of shear bands in metallic glasses, Applied Physics Letters 89 , 071907 (2006).
- 901 (2003). [42] Y. Cheng, A. Cao, and E. Ma, Correlation between the elastic modulus and the intrinsic plastic behavior of metallic glasses: The roles of atomic configuration and alloy composition, Acta Materialia 57 , 3253 (2009).
- [28] F. Shimizu, S. Ogata, and J. Li, Theory of shear banding in metallic glasses and molecular dynamics calculations, Materials Transactions 48 , 2923 (2007).
- [29] M. L. Falk and J. S. Langer, Dynamics of viscoplastic deformation in amorphous solids, Physical Review E 57 , 7192 (1998).
- [43] K. Ito, C. T. Moynihan, and C. A. Angell, Thermodynamic determination of fragility in liquids and a fragile-to-strong liquid transition in water, Nature 398 , 492 (1999).
- [30] M. Oquab, T. Darcet, T. Moutakanni, H. V. Vo, M. Szafraniec, V. Khalidov, P. Fernandez, D. HAZIZA, F. Massa, A. El-Nouby, M. Assran, N. Ballas, W. Galuba,
- [44] P. Lucas, Fragile-to-strong transitions in glass forming liquids, Journal of Non-Crystalline Solids: X 4 , 100034 (2019).

- [45] C. Way, P. Wadhwa, and R. Busch, The influence of shear rate and temperature on the viscosity and fragility of the Zr 41.2 Ti 13.8 Cu 12.5 Ni 10.0 Be 22.5 metallic-glass-forming liquid, Acta Materialia 55 , 2977 (2007).
- [46] A. Atila, M. Kbirou, S. Ouaskit, and A. Hasnaoui, On the presence of nanoscale heterogeneity in Al 70 Ni 15 Co 15 metallic glass under pressure, Journal of Non-Crystalline Solids 550 , 120381 (2020).
- [47] A. Atila, E. M. Ghardi, S. Ouaskit, and A. Hasnaoui, Atomistic insights into the impact of charge balancing cations on the structure and properties of aluminosilicate glasses, Physical Review B 100 , 10.1103/physrevb.100.144109 (2019).
- [48] A. Atila, Y. Ouldhnini, S. Ouaskit, and A. Hasnaoui, Atomistic insights into the mixed-alkali effect in phosphosilicate glasses, Physical Review B 105 , 10.1103/physrevb.105.134101 (2022).
- [49] C. Packard and C. Schuh, Initiation of shear bands near a stress concentration in metallic glass, Acta Materialia 55 , 5348 (2007).
- [50] J. Wang, P. D. Hodgson, J. Zhang, W. Yan, and C. Yang, Effects of pores on shear bands in metallic glasses: A molecular dynamics study, Computational Materials Science 50 , 211 (2010).
- [51] K. Kosiba, S. Scudino, J. Bednarcik, J. Bian, G. Liu, U. K¨ uhn, and S. Pauly, Guiding shear bands in bulk metallic glasses using stress fields: A perspective from the activation of flow units, Phys. Rev. B 102 , 134113 (2020).
- [52] L. Tang, H. Liu, G. Ma, T. Du, N. Mousseau, W. Zhou, and M. Bauchy, The energy landscape governs ductility in disordered materials, Materials Horizons 8 , 1242 (2021).
- [53] Z. Liu, Y. Yang, and C. Liu, Yielding and shear banding of metallic glasses, Acta Materialia 61 , 5928 (2013).
- [54] Y. Shi and M. L. Falk, Simulations of nanoindentation in a thin amorphous metal film, Thin Solid Films 515 , 3179 (2007).
- [55] A. Antoniou, A. Bastawros, and B. Biner, Experimental observations of deformation behavior of bulk metallic glasses during wedge-like cylindrical indentation, Journal of Materials Research 22 , 514 (2007).
- [56] A. Antoniou, A. Bastawros, C. Lo, and S. Biner, Deformation behavior of a zirconium based metallic glass during cylindrical indentation: in situ observations, Materials Science and Engineering: A 394 , 96 (2005).
- [57] C. Su and L. Anand, Plane strain indentation of a zrbased metallic glass: Experiments and numerical simulation, Acta Materialia 54 , 179 (2006).

## Supplementary material to:

## Brittleness of metallic glasses dictated by their state at the fragile-to-strong transition temperature

FIG. S1. Local shear strain maps for the strong and fragile glasses loaded to different indentation depths.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00536v1-figure-bc76f2cb-1ad6-43f2-a20c-2bda9c56a15e.jpg)

FIG. S2. Local non-affine displacement squared calculated as detailed in the main text for the strong and fragile glasses loaded to an indentation depth of (a) d = 10 nm, (b) d = 12 nm, and (c) d = 15 nm.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00536v1-figure-042406d2-4334-40c1-95ac-77ab7b512805.jpg)
Achraf Atila Sergey V. Sukhomlinov Marc J. Honecker Martin H. Müser
Brittleness of metallic glasses dictated by their state at the fragile-to-strong transition temperature
The effect of cooling on the brittleness of glasses in general, and bulk
metallic glasses (BMGs) in particular, is usually studied with continuously
varying cooling rates; slower cooling rates lead to stiffer, harder, and more
brittle glasses than higher cooling rates. These protocols obscure any
potential discontinuity that a glass might experience depending on whether its
microstructure resembles that of a fragile or a strong glass-forming liquid.
Here, we use large-scale molecular dynamics to simulate the nanoindentation
behavior of model BMGs (Zr$_{0.6}$Cu$_{0.3}$Al$_{0.1}$) obtained by rapidly
quenching equilibrium melts from temperatures above and below the
fragile-to-strong transition temperature $T_\textrm{fst}$, leading to fragile
and strong glasses, respectively. While the contact modulus deduced from the
indentation simulation evolves smoothly with the temperature $T_\text{q}$ from
which the melt is quenched, the plastic response changes quasi-discontinuously
as $T_\text{q}$ passes through $T_\textrm{fst}$. In particular, strong glasses
develop highly asymmetric flow profiles with mature shear bands, in contrast to
fragile glasses. Quantitative differences reveal themselves not only through a
formal von Mises localization parameter analysis but also through image
analysis of flow patterns using pre-trained artificial intelligence models.
Moreover, seemingly erratic flow profiles for our indentation geometry produced
surprisingly reproducible and, thus, deterministic features. It remains to be
determined to what extent other classes of glass formers follow our observation
that the degree of brittleness is significantly influenced by whether the melt
is fragile or strong when it falls out of equilibrium at the glass transition
temperature.