## The effect of thermal history on the atomic structure and mechanical properties of amorphous alloys

Nikolai V. Priezjev 1 2 ,

1

Department of Mechanical and Materials Engineering,

Wright State University,

Dayton, OH 45435 and

2 National Research University Higher School of Economics, Moscow 101000, Russia

(Dated: August 2, 2019)

## Abstract

The influence of thermal processing on the potential energy, atomic structure, and mechanical properties of metallic glasses is examined using molecular dynamics simulations. We study the three-dimensional binary mixture, which was first relaxed near the glass transition temperature, and then rapidly cooled deep into the glass phase. It was found that glasses prepared at higher annealing temperatures are relocated to higher energy states and their average glass structure remains more disordered, as reflected in the shape of the pair correlation function. The results of mechanical testing demonstrate that both the shear modulus and yielding peak increase significantly when the annealing temperature approaches T g from above. Moreover, the shear modulus becomes a strong function of strain rate only for samples equilibrated at sufficiently high temperatures. Based on the spatial distribution of nonaffine displacements, we show that the deformation mode changes from brittle to ductile upon increasing annealing temperature. These results can be useful for the design and optimization of the fabrication processes of bulk glassy alloys with improved plasticity.

Keywords: metallic glasses, deformation, thermal treatment, yield stress, molecular dynamics simulations

## I. INTRODUCTION

Recent progress in the thermal treatment and deformation processing of amorphous alloys has resulted in a wider range of accessible energy states and optimized physical and mechanical properties [1]. It is well known that metallic glasses are much stronger than crystalline materials due to the absence of topological defects, or dislocations, but their deformation often proceeds via the formation of the so-called shear bands. Generally, aged samples become more brittle, while rejuvenated, or higher energy glasses, are more ductile. Common methods to rejuvenate metallic glasses and improve their plasticity include cold rolling, high pressure torsion, ion irradiation, and shot peening [2]. Interestingly, an application of elastostatic loading over extended time intervals can significantly rejuvenate metallic glasses due to plastic deformation of relatively soft domains [3-9]. A novel method to control rejuvenation consists of heating up metallic glasses slightly above the glass transition temperature, followed by rapid quenching back to the glass phase [10-13]. More recently, it was found that cryogenic thermal cycling might lead to higher energy states via local plastic events that are triggered because of the heterogeneity in the local thermal expansion [14-18]. In addition, it was demonstrated that the atomic structure and internal energy of amorphous materials can be tuned by cyclic deformation [19-31]. Finally, recent experimental studies showed that the ductile-to-brittle transition at room temperature is controlled by the socalled fictive temperature, which characterizes the average glass structure during a rapid quench below the glass transition temperature [32, 33]. However, despite extensive efforts, the effect of thermo-mechanical processing on microstructure and mechanical properties of amorphous alloys remains not fully understood.

In this paper, we investigate the dependence of the potential energy, atomic structure and mechanical properties of binary glasses on the preparation protocol that consists of equilibration at an annealing temperature near the glass transition point and a subsequent quench to the test temperature. It will be shown that glasses prepared at lower annealing temperatures settle down at deeper potential energy minima and their short-range structure becomes more ordered. As a result, the shear modulus and stress overshoot depend strongly on the annealing temperature. We find that the rate dependence of the shear modulus becomes significant when the annealing temperature is greater than the glass transition temperature. The analysis of nonaffine displacements indicates that upon increasing an-

nealing temperature, deformation changes from being localized within a shear band to more homogeneous.

This paper is organized as follows. The description of the molecular dynamics simulation model as well as the preparation and deformation protocols is provided in the next section. The dependence of the potential energy, radial distribution function, mechanical properties, and distribution of nonaffine displacements on the annealing temperature is presented in section III. A brief summary is given in the last section.

## II. MOLECULAR DYNAMICS SIMULATION MODEL

In this work, the structure and mechanical properties of an amorphous alloy are investigated by means of molecular dynamics simulations. We consider a well-studied binary (80:20) mixture originally introduced by Kob and Andersen (KA) about two decades ago [34]. The interaction parameters of the KA mixture are similar to the parametrization used by Weber and Stillinger to study the amorphous metal alloy Ni 80 P 20 [35]. More specifically, the atoms of types α, β = A, B in the KA model interact via the Lennard-Jones (LJ) potential:

$$V _ { \alpha \beta } ( r ) = 4 \, \varepsilon _ { \alpha \beta } \left [ \left ( \frac { \sigma _ { \alpha \beta } } { r } \right ) ^ { 1 2 } - \left ( \frac { \sigma _ { \alpha \beta } } { r } \right ) ^ { 6 } \right ],$$

with the following LJ parameters: ε AA = 1 0, . ε AB = 1 5, . ε BB = 0 5, . σ AA = 1 0, . σ AB = 0 8, . σ BB = 0 88, . and m A = m B [34]. Note that this parametrization represents a system with strongly non-additive interaction between different types of atoms, which suppresses crystallization near the glass transition temperature [34]. The total number of atoms is N = 60000. To alleviate the computational load, the LJ potential was truncated at the cutoff radius r c, αβ = 2 5 . σ αβ . In what follows, all physical quantities are expressed in the LJ units of length, mass, energy, and time; namely, σ = σ AA , m = m A , ε = ε AA , and τ = σ √ m/ε . The MD simulations were performed using the LAMMPS open-source, parallel code [36] with the integration time step glyph[triangle] t MD = 0 005 . τ .

The first step in the preparation procedure was to arrange all atoms at the sites of the fcc lattice and then equilibrate the binary mixture at the temperature T LJ = 1 0 . ε/k B and zero pressure. Here, k B and T LJ indicate the Boltzmann constant and temperature, respectively. In all simulations, the system temperature was regulated via the Nos´-Hoover e thermostat [36, 37]. Periodic boundary conditions were imposed along all three dimensions.

Next, the system was linearly cooled down to an annealing temperature T a (in the vicinity of the glass transition temperature) during 10 4 τ at zero pressure. Once at T a , the binary mixture was relaxed at constant temperature and pressure during the time interval 2 × 10 5 τ (i.e., 4 × 10 7 MD steps). Following the relaxation period, the system was instantaneously quenched to the very low temperature T LJ = 0 01 . ε/k B and then further relaxed during the additional time interval 10 4 τ at zero pressure.

After the thermal treatment, the mechanical properties were obtained by imposing shear strain deformation at constant volume and temperature T LJ = 0 01 . ε/k B . The shear modulus and the peak value of the stress overshoot were computed from the stress-strain curves, and the data were averaged over 15 independent samples. The deformation procedure was repeated for different values of the strain rate. This required substantial computational resources (about 1000 processors). In all simulations, the system temperature, potential energy, stress components, and atomic configurations were periodically saved for the postprocessing analysis.

## III. RESULTS

It has been long realized that physical and mechanical properties of metallic glasses depend crucially on the annealing history [2]. An auxiliary concept to characterize the structural degrees of freedom is the so-called fictive temperature defined as the temperature at which the glass would be in equilibrium if suddenly brought to it from its current state [38, 39]. In the limiting case of an infinitely fast quenching rate, used in the present study, the fictive temperature essentially coincides with the annealing temperature at which the system is equilibrated before the quench. Interestingly, it was demonstrated experimentally for different metallic glasses that an abrupt toughening transition between brittle and ductile regimes occurs as functions of the fictive temperature and strain rate [33]. It was argued that the mechanical transition originates from the competition between two time scales; the loading time, which is inversely proportional to the applied strain rate, and the plastic deformation time scale, determined by the density of plasticity carriers or shear transformation zones [33].

In the recent MD study, it was found that the glass transition temperature of the KA mixture of N = 60000 particles at zero pressure is T g ≈ 0 35 . ε/k B [13]. This value was

obtained by linearly extrapolating the low and high temperature dependence of the potential energy during gradual cooling with the computationally slow rate of 10 -5 ε/k B τ . Note that higher cooling rates generally result in an increase of the glass transition temperature, which can be understood from the assumption that the relaxation time of the system follows a Vogel-Fulcher dependence on temperature [41]. For reference, the mode-coupling critical temperature of the KA model at the fixed density ρ = 1 2 . σ -3 is T c = 0 435 . ε/k B , which was determined from the fit of the diffusion coefficients of either A or B types of particles to the power-law functions of temperature at constant volume [34].

In our study, the system was relaxed at zero pressure in the range of annealing temperatures, 0 31 . ε/k B glyph[lessorequalslant] T a glyph[lessorequalslant] 0 50 . ε/k B , near the glass transition temperature. The variation of the potential energy per particle during the relaxation time interval of 2 × 10 5 τ and a subsequent quench to T LJ = 0 01 . ε/k B is shown in Fig. 1 for selected values of the annealing temperature. It can be observed that the potential energy remains nearly constant at T a glyph[greaterorsimilar] 0 38 . ε/k B , while it decreases gradually during aging at lower temperatures. To estimate the degree of structural relaxation during 2 × 10 5 τ , we computed the mean square displacement (MSD) of atoms at low annealing temperatures. For example, the averaged MSD is r 2 ≈ 1 64 . σ 2 at T a = 0 32 . ε/k B , indicating that each atom undergoes a number of cage jumps during the relaxation time interval. This estimate is based on two atomic configurations separated by the time interval 2 × 10 5 τ .

After the system was equilibrated at a certain temperature near T g , it was instantaneously quenched to the test temperature T LJ = 0 01 . ε/k B and then relaxed during 10 4 τ at zero pressure. This procedure aims at freezing out the thermal vibrations but keeping intact the average glass structure. The inset in Fig. 1 shows the dependence of the potential energy after 10 4 τ at T LJ = 0 01 . ε/k B as a function of the annealing temperature. It can be seen that the potential energy, which is mostly determined by the glass structure, decreases towards a minimum when T a approaches the glass transition temperature from above. These results confirm that the potential energy of inherent structures becomes lower when the temperature is reduced towards the glass transition temperature [20, 42]. Test simulations at lower annealing temperatures, 0 28 . ε/k B glyph[lessorequalslant] T a glyph[lessorequalslant] 0 30 . ε/k B , have shown that U , measured at T LJ = 0 01 . ε/k B , continues to increase upon further reducing T a (not shown). This trend can be explained by realizing that the effective cooling rate increases when the system

temperature is reduced from T LJ = 1 0 . ε/k B to lower values of T a during the fixed time interval of 10 4 τ . This way the system is annealed to a glass phase with higher energy states.

We next analyze the changes in the glass structure for systems prepared at different annealing temperatures. It was previously realized that the most sensitive measure of the atomic structure of the KA model glass is the shape of the pair correlation function of the smaller atoms of type B [40, 41]. This conclusion stems from the fact that the interaction energy, ε BB in Eq. (1), is the lowest among the three energies, and, therefore, the number of nearest-neighbor contacts between smaller, more mobile, atoms is reduced at lower temperatures [40]. In our study, the averaged pair distribution function is plotted in Fig. 2 for systems prepared at T a = 0 32 . ε/k B , 0 38 . ε/k B , and 0 50 . ε/k B . It can be readily observed that the height of the first peak is significantly reduced when the annealing temperature is varied from 0 50 . ε/k B to 0 32 . ε/k B , while the opposite trend is seen for the second peak. We comment that the variation of the first peak height in Fig. 2 is much more pronounced than the change reported during the heat treatment process, where the glass structure was reset by heating above T g and then the system was cooled back to the glass phase with different rates, resulting in slightly rejuvenated states [13].

The examples of stress-strain curves during shear deformation with the strain rate ˙ ε xz = 10 -5 τ -1 are presented in Fig. 3 for selected values of T a . The deformation was imposed at constant volume and temperature T LJ = 0 01 . ε/k B up to 20% shear strain. As is evident, the stress response gradually changes from a monotonic increase and saturation at T a = 0 50 . ε/k B to a pronounced overshoot followed by an abrupt stress drop at T a = 0 32 . ε/k B . Such behavior typically reflects either a uniform plastic flow of poorly annealed glasses or the formation of a shear band in well annealed samples (see discussion below).

The stress-strain curves were used to estimate the peak value of the stress overshoot as a function of T a . In turn, the shear modulus was computed from the slope of the best linear fit to the data at ε xz glyph[lessorequalslant] 0 01. . Both σ Y and G , averaged over 15 samples, are plotted in Fig. 4 for the indicated strain rates. It can be clearly seen that the mechanical properties depend strongly on the annealing temperature and the transition occurs when T a ≈ T g . Furthermore, it was recently shown experimentally that a toughening transition as a function of strain rate becomes more pronounced when the fictive temperature is slightly higher than T g [33]. In our MD setup, the yielding peak becomes higher with increasing strain rate,

as expected; but the increase in σ Y is rather insensitive to T a (see Fig. 4). On the other hand, as shown in the inset to Fig. 4, the shear modulus is a strong function of strain rate only when T a &gt; T g . In this sense, the results of MD simulations are consistent with the conclusions of the experimental study [33]. It should also be commented that it is difficult to access strain rates outside the range 10 -5 glyph[lessorequalslant] ˙ ε xz τ glyph[lessorequalslant] 10 -3 reported in Fig. 4. Specifically, the averaging becomes computationally expensive at lower strain rates, while temperature oscillations are detected at higher strain rates.

The microscopic details of the deformation process can be captured via the analysis of the so-called nonaffine displacements of atoms [43]. It is well realized that the local displacement of an atom with respect to its neighbors can be computed via the transformation matrix J i , which linearly transforms a group of atoms during the time interval ∆ t and at the same time minimizes the following quantity:

$$D ^ { 2 } ( t, \Delta t ) = \frac { 1 } { N _ { i } } \sum _ { j = 1 } ^ { N _ { i } } \left \{ { \mathbf r } _ { j } ( t + \Delta t ) - { \mathbf r } _ { i } ( t + \Delta t ) - { \mathbf J } _ { i } \left [ { \mathbf r } _ { j } ( t ) - { \mathbf r } _ { i } ( t ) \right ] \right \} ^ { 2 },$$

where the summation is carried over the neighbors within a sphere of radius 1 5 . σ with the center at the position vector r i ( ). t It was demonstrated that the nonaffine measure is an excellent diagnostic of localized shear transformations in deformed amorphous materials [43]. Thus, recent studies on periodic deformation of disordered solids have shown that the amplitude of the nonaffine measure is approximately power-law distributed, while the majority of atoms rearrange reversibly during one or several subyield cycles [24, 25]. Above the yielding point, the temporal evolution of spatial distributions of atoms with large nonaffine displacements illustrated the existence of a transient regime of periodic deformation, followed by the formation of shear bands in both poorly [29] and well [27] annealed binary glasses. More recently, the analysis of nonaffine displacements was used to study aging and rejuvenation during elastostatic loading [9] as well as the mechanical [28, 29, 31] and thermal [16] annealing processes in amorphous alloys.

The spatial distribution of nonaffine displacements during steady shear deformation with the strain rate ˙ ε xz = 10 -5 τ -1 is plotted in Fig. 5 for T a = 0 32 . ε/k B , in Fig. 6 for 0 38 . ε/k B , and in Fig. 7 for 0 50 . ε/k B . The sequence of instantaneous snapshots are taken at shear strains ε xz = 0 05, 0 10, 0 15, and 0 20, and the nonaddine measure, Eq. (2), was computed . . . . with respect to the atomic configuration at zero strain. It can be clearly seen in Fig. 5

that the shear deformation is accompanied by relatively small clusters of atoms with large nonaffine displacements at ε xz glyph[lessorequalslant] 0 10, . and the shear band is formed and become fully developed along the xy plane at higher strains. This behavior is consistent with the stress response shown in Fig. 3 for T a = 0 32 . ε/k B , where the stress drop at ε xz ≈ 0 11 is associated . with the formation of the shear band.

Next, as illustrated in Fig. 6, the deformation pattern remains similar for the sample prepared at the higher annealing temperature T a = 0 38 . ε/k B , except that the shear band appears to be thinner and it is oriented perpendicular to the plane of shear. Such orientation is allowed due to periodic boundary conditions, and it was observed in the previous MD simulations of binary glasses under steady and periodic shear deformation [27, 44]. In sharp contrast, a more homogeneous plastic deformation is detected in glasses equilibrated at T a = 0 50 . ε/k B , as shown in Fig. 7. The inherent structure at this annealing temperature has the highest potential energy level (see the inset in Fig. 1) and, therefore, the largest number of shear transformation zones with low instability thresholds, resulting in homogeneous deformation. Altogether, these results demonstrate the importance of controlling the annealing temperature in order to optimize the mechanical properties of bulk metallic glasses and improve their plasticity.

## IV. CONCLUSIONS

In summary, molecular dynamics simulations were carried out to investigate the influence of the thermal treatment on structure and mechanical properties of disordered solids. A model glass was represented by the three-dimensional binary mixture with highly nonadditive parametrization that prevents crystallization near the glass transition temperature. The binary mixture was first equilibrated at a temperature close to the glass transition point and then instantaneously quenched and further relaxed at a very low temperature. It was shown that the potential energy of the inherent structures is reduced and the height of the nearest-neighbor peak in the pair distribution function of smaller atoms becomes lower in samples prepared at lower annealing temperatures. The results of mechanical tests indicate that both the shear modulus and the yielding peak increase significantly when the annealing temperature is reduced towards the glass transition temperature. We also found that the strain rate dependence of the shear modulus becomes apparent only at anneal-

ing temperatures greater than the glass transition temperature. The analysis of nonaffine displacements illustrated that the shear deformation changed from brittle to ductile upon increasing annealing temperature.

## Acknowledgments

Financial support from the National Science Foundation (CNS-1531923) is gratefully acknowledged. The article was prepared within the framework of the HSE University Basic Research Program and funded in part by the Russian Academic Excellence Project '5-100'. The molecular dynamics simulations were performed using the LAMMPS software developed at Sandia National Laboratories [36]. The numerical simulations were performed at Wright State University's Computing Facility and the Ohio Supercomputer Center.

- [1] J. C. Qiao, Q. Wang, J. M. Pelletier, H. Kato, R. Casalini, D. Crespo, E. Pineda, Y. Yao, Y. Yang, Structural heterogeneities and mechanical behavior of amorphous alloys, Prog. Mater. Sci. 104 , 250 (2019).
- [2] Y. Sun, A. Concustell, and A. L. Greer, Thermomechanical processing of metallic glasses: Extending the range of the glassy state, Nat. Rev. Mater. 1 , 16039 (2016).
- [3] K.-W. Park, C.-M. Lee, M. Wakeda, Y. Shibutani, M. L. Falk, and J.-C. Lee, Elastostatically induced structural disordering in amorphous alloys, Acta Materialia 56 , 5440 (2008).
- [4] Y. Tong, W. Dmowski, Y. Yokoyama, G. Wang, P. K. Liaw, and T. Egami, Recovering compressive plasticity of bulk metallic glasses by high-temperature creep, Scripta Materialia 69 570 (2013).
- [5] A. L. Greer, and Y. H. Sun, Stored energy in metallic glasses due to strains within the elastic limit, Philos. Mag. 96 , 1643 (2016).
- [6] M. Zhang, Y. M. Wang, F. X. Li, S. Q. Jiang, M. Z. Li and L. Liu, Mechanical relaxation-torejuvenation transition in a Zr-based bulk metallic glass, Sci. Rep. 7 , 625 (2017).
- [7] J. Pan, Y. X. Wang, Q. Guo, D. Zhang, A. L. Greer, and Y. Li, Extreme rejuvenation and softening in a bulk metallic glass, Nat. Commun. 9 , 560 (2018).
- [8] M. Samavatian, R. Gholamipour, A. A. Amadeh, and S. Mirdamadi, Role of tensile elastostatic

loading on atomic structure and mechanical properties of Zr 55 Cu 30 Ni Al 5 10 bulk metallic glass, Mater. Sci. Eng. A 753 , 218 (2019).

- [9] N. V. Priezjev, Aging and rejuvenation during elastostatic loading of amorphous alloys: A molecular dynamics simulation study, Comput. Mater. Sci. 168 , 125 (2019).
- [10] M. Wakeda, J. Saida, J. Li, and S. Ogata, Controlled rejuvenation of amorphous metals with thermal processing, Sci. Rep. 5 , 10545 (2015).
- [11] S. Kuchemann, P. M. Derlet, C. Liu, D. Rosenthal, G. Sparks, W. S. Larson, and R. Maass, Energy Storage in Metallic Glasses via Flash Annealing, Adv. Funct. Mater. 28 , 1805385 (2018).
- [12] M. Wang, H. Liu, J. Mo, Y. Zhang, Z. Chen, C. Yin, and W. Yang, Thermal-pressure effects on energy state of metallic glass Cu 50 Zr 50 , Comput. Mater. Sci. 155 , 493 (2018).
- [13] N. V. Priezjev, Atomistic modeling of heat treatment processes for tuning the mechanical properties of disordered solids, J. Non-Cryst. Solids 518 , 128 (2019).
- [14] S. V. Ketov, Y. H. Sun, S. Nachum, Z. Lu, A. Checchi, A. R. Beraldin, H. Y. Bai, W. H. Wang, D. V. Louzguine-Luzgin, M. A. Carpenter, and A. L. Greer, Rejuvenation of metallic glasses by non-affine thermal strain, Nature 524 , 200 (2015).
- [15] B. Shang, P. Guan, and J.-L. Barrat, Role of thermal expansion heterogeneity in the cryogenic rejuvenation of metallic glasses, J. Phys.: Mater. 1 , 015001 (2018).
- [16] N. V. Priezjev, The effect of cryogenic thermal cycling on aging, rejuvenation, and mechanical properties of metallic glasses, J. Non-Cryst. Solids 503 , 131 (2019).
- [17] Q.-L. Liu and N. V. Priezjev, The influence of complex thermal treatment on mechanical properties of amorphous materials, Comput. Mater. Sci. 161 , 93 (2019).
- [18] W. Guo, J. Saida, M. Zhao, S. Lu, and S. Wu, Rejuvenation of Zr-based bulk metallic glass matrix composite upon deep cryogenic cycling, Materials Letters 247 , 135 (2019).
- [19] N. V. Priezjev, Heterogeneous relaxation dynamics in amorphous materials under cyclic loading, Phys. Rev. E 87 , 052302 (2013).
- [20] D. Fiocco, G. Foffi, and S. Sastry, Oscillatory athermal quasistatic deformation of a model glass, Phys. Rev. E 88 , 020301(R) (2013).
- [21] I. Regev, T. Lookman, and C. Reichhardt, Onset of irreversibility and chaos in amorphous solids under periodic shear, Phys. Rev. E 88 , 062401 (2013).
- [22] N. V. Priezjev, Dynamical heterogeneity in periodically deformed polymer glasses, Phys. Rev.

- E 89 , 012601 (2014).
- [23] I. Regev, J. Weber, C. Reichhardt, K. A. Dahmen, and T. Lookman, Reversibility and criticality in amorphous solids, Nat. Commun. 6 , 8805 (2015).
- [24] N. V. Priezjev, Reversible plastic events during oscillatory deformation of amorphous solids, Phys. Rev. E 93 , 013001 (2016).
- [25] N. V. Priezjev, Nonaffine rearrangements of atoms in deformed and quiescent binary glasses, Phys. Rev. E 94 , 023004 (2016).
- [26] P. Leishangthem, A. D. S. Parmar, and S. Sastry, The yielding transition in amorphous solids under oscillatory shear deformation, Nat. Commun. 8 , 14653 (2017).
- [27] N. V. Priezjev, Collective nonaffine displacements in amorphous materials during largeamplitude oscillatory shear, Phys. Rev. E 95 , 023002 (2017).
- [28] N. V. Priezjev, Molecular dynamics simulations of the mechanical annealing process in metallic glasses: Effects of strain amplitude and temperature, J. Non-Cryst. Solids 479 , 42 (2018).
- [29] N. V. Priezjev, The yielding transition in periodically sheared binary glasses at finite temperature, Comput. Mater. Sci. 150 , 162 (2018).
- [30] Y. C. Lo, H. S. Chou, Y. T. Cheng, J. C. Huang, J. R. Morris, P. K. Liaw, Structural relaxation and self-repair behavior in nano-scaled Zr-Cu metallic glass under cyclic loading: Molecular dynamics simulations, Intermetallics 18 , 954 (2010).
- [31] N. V. Priezjev, Slow relaxation dynamics in binary glasses during stress-controlled, tensioncompression cyclic loading, Comput. Mater. Sci. 153 , 235 (2018).
- [32] G. Kumar, P. Neibecker, Y. H. Liu, and J. Schroers, Critical fictive temperature for plasticity in metallic glasses, Nat. Commun. 4 , 1536 (2013).
- [33] J. Ketkaew, W. Chen, H. Wang, A. Datye, M. Fan, G. Pereira, U. D. Schwarz, Z. Liu, R. Yamada, W. Dmowski, M. D. Shattuck, C. S. O'Hern, T. Egami, E. Bouchbinder, and J. Schroers, Mechanical glass transition revealed by the fracture toughness of metallic glasses, Nat. Commun. 9 , 3271 (2018).
- [34] W. Kob and H. C. Andersen, Testing mode-coupling theory for a supercooled binary LennardJones mixture: The van Hove correlation function, Phys. Rev. E 51 , 4626 (1995).
- [35] T. A. Weber and F. H. Stillinger, Local order and structural transitions in amorphous metalmetalloid alloys, Phys. Rev. B 31 , 1954 (1985).
- [36] S. J. Plimpton, Fast parallel algorithms for short-range molecular dynamics, J. Comp. Phys.

## 117 , 1 (1995).

- [37] M. P. Allen and D. J. Tildesley, Computer Simulation of Liquids (Clarendon, Oxford, 1987).
- [38] A. Q. Tool, Relation between inelastic deformability and thermal expansion of glass in its annealing range, J. Am. Ceram. Soc. 29 , 240 (1946).
- [39] S. Richter, F. Zimmermann, A. Tunnermann, and S. Nolte, Laser welding of glasses at high repetition rates - Fundamentals and prospects, Optics and Laser Technology 83 , 59 (2016).
- [40] M. Utz, P. G. Debenedetti, and F. H. Stillinger, Atomistic simulation of aging and rejuvenation in glasses, Phys. Rev. Lett. 84 , 1471 (2000).
- [41] K. Vollmayr, W. Kob, and K. Binder, How do the properties of a glass depend on the cooling rate? A computer simulation study of a Lennard-Jones system, J. Chem. Phys. 105 , 4714 (1996).
- [42] F. H. Stillinger and T. A. Weber, Packing structures and transitions in liquids and solids, Science 225 , 983 (1984).
- [43] M. L. Falk and J. S. Langer, Dynamics of viscoplastic deformation in amorphous solids, Phys. Rev. E 57 , 7192 (1998).
- [44] G. P. Shrivastav, P. Chaudhuri, and J. Horbach, Yielding of glass under shear: A directed percolation transition precedes shear-band formation, Phys. Rev. E 94 , 042605 (2016).

## Figures

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00464v1-figure-b76b4f9a-aaed-4904-ac49-85b29ba4dd7b.jpg)

/

FIG. 1: (Color online) The time dependence of the potential energy, U/ε , at various annealing temperatures, T a , and after the quench to T LJ = 0 01 . ε/k B at P = 0. The values of the annealing temperature are T a k B /ε = 0 32, 0 34, 0 36, 0 38, 0 40, 0 42, 0 44, 0 46, 0 48, and 0 50 (from bottom . . . . . . . . . . to top). The inset shows the potential energy at T LJ = 0 01 . ε/k B after the time interval of 10 4 τ . The data in the inset are averaged over 15 samples.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00464v1-figure-7464cf8b-5f00-479f-be8b-2ee02050a8ba.jpg)

/

FIG. 2: (Color online) The pair distribution function of smaller atoms of type B g , BB ( r ), for glasses prepared at different annealing temperatures. The values of T a (in units of ε/k B ) are listed in the legend. The distribution function was computed after the samples were instantaneously quenched from T a to T LJ = 0 01 . ε/k B and then relaxed during the time interval 10 4 τ at P = 0. The data are averaged over 15 independent samples for each value of T a .

σ

FIG. 3: (Color online) The variation of shear stress, σ xz (in units of εσ -3 ), as a function of strain, ε xz , for the indicated values of the annealing temperature (see the legend). The samples were strained along the xz plane with the strain rate ˙ ε xz = 10 -5 τ -1 at T LJ = 0 01 . ε/k B and constant volume.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00464v1-figure-278b5257-a471-4fdd-bcb8-9e7c9178cd4c.jpg)

σ

FIG. 4: (Color online) The peak value of the stress overshoot σ Y (in units of εσ -3 ) as a function of the annealing temperature for the indicated values of the strain rate, ˙ ε xz (in units of τ -1 ). The inset shows the shear modulus G (in units of εσ -3 ) versus T a for the same strain rates. The data for σ Y and G are averaged over 15 samples.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00464v1-figure-1832ac67-fdf5-48d2-bb9f-0e3c119b55cb.jpg)

FIG. 5: (Color online) The sequence of snapshots for the sample prepared at T a = 0 32 . ε/k B and then quenched to T LJ = 0 01 . ε/k B . The shear strain ε xz is (a) 0 05, (b) 0 10, (c) 0 15, and (d) 0 20. . . . . The rate of strain is ˙ ε xz = 10 -5 τ -1 . The color denotes the nonaffine measure D 2 with respect to zero strain (see the legend).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00464v1-figure-0e396233-7492-4c8c-b807-a27a2101fae9.jpg)

FIG. 6: (Color online) The series of snapshots depicting atomic configurations in the strained glass prepared at the annealing temperature T a = 0 38 . ε/k B . The shear strain is (a) 0 05, (b) 0 10, (c) . . 0 15, and (d) 0 20. . . The nonaffine measure D 2 with respect to zero strain is marked according to the legend. The glass is strained at T LJ = 0 01 . ε/k B and constant volume with the rate ˙ ε xz = 10 -5 τ -1 . The atoms are not drawn to scale.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00464v1-figure-4188277c-a9b2-403c-8a67-4671529262a2.jpg)

FIG. 7: (Color online) The consecutive snapshots of the glass equilibrated at T a = 0 50 . ε/k B and then quenched to T LJ = 0 01 . ε/k B . The strain is (a) 0 05, (b) 0 10, (c) 0 15, and (d) 0 20. . . . . The nonaffine displacements of atoms with respect to their neighbors at zero strain are indicated by the color according to the figure legend. The strain rate is ˙ ε xz = 10 -5 τ -1 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00464v1-figure-3e9be3c1-bfc6-49b4-8d29-65e322ff7c6a.jpg)
Nikolai V. Priezjev
The effect of thermal history on the atomic structure and mechanical properties of amorphous alloys
The influence of thermal processing on the potential energy, atomic
structure, and mechanical properties of metallic glasses is examined using
molecular dynamics simulations. We study the three-dimensional binary mixture,
which was first relaxed near the glass transition temperature, and then rapidly
cooled deep into the glass phase. It was found that glasses prepared at higher
annealing temperatures are relocated to higher energy states and their average
glass structure remains more disordered, as reflected in the shape of the pair
correlation function. The results of mechanical testing demonstrate that both
the shear modulus and yielding peak increase significantly when the annealing
temperature approaches $T_g$ from above. Moreover, the shear modulus becomes a
strong function of strain rate only for samples equilibrated at sufficiently
high temperatures. Based on the spatial distribution of nonaffine
displacements, we show that the deformation mode changes from brittle to
ductile upon increasing annealing temperature. These results can be useful for
the design and optimization of the fabrication processes of bulk glassy alloys
with improved plasticity.