Journal of Nanomaterials

Volume 2014 (2014), Article ID 318757, 6 pages

http://dx.doi.org/10.1155/2014/318757

Research Article

Molecular Dynamics Simulation of Structural Characterization of Elastic and Inelastic Deformation in ZrCu Metallic Glasses

Shidong Feng,¹ Li Qi,¹ Gong Li,1,2 and Riping Liu¹

¹State Key Laboratory of Metastable Materials Science and Technology, Yanshan University, Qinhuangdao 066004, China

²Department of Materials Science and Engineering, The University of Tennessee, Knoxville, TN 37996, USA

Received 4 March 2014; Accepted 10 May 2014; Published 27 May 2014

Academic Editor: Alireza Khataee

Copyright © 2014 Shidong Feng et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The nanoscopic deformation behaviors in a ZrCu metallic glass model during loading-unloading process under uniaxial compression have been analyzed on the basis of the molecular dynamics (MD). The reversible degree of shear origin zones (SOZs) is used as the structural indicator to distinguish the elastic deformation and inelastic deformation of ZrCu metallic glass at the atomic level. We find that the formation of SOZs is reversible at the elastic stage but irreversible at the inelastic stage during the loading and unloading processes. At the inelastic stage, the full-icosahedra fraction in SOZs is quickly reduced with increased strain and the decreasing process is also irreversible during the unloading processes.

1. Introduction

Bulk metallic glasses (BMGs) have attracted much attention as novel potential engineering materials because of their specific mechanical properties [1–3]. However, compared with their crystalline counterparts, many extraordinary behaviors in glasses are still poorly understood [4–6]. One of the most well-known features is the very high elastic strain limit of BMGs, which is about one order of magnitude higher than that in conventional crystalline metals and alloys [7–9]. However, it is known that the global plasticity of BMGs is very low [10, 11], which is caused by avalanche of shear transformation zones (STZs). The shear transformation zone (STZ) was hypothesized as the small cluster containing tens of atoms confined within elastic medium [11]. High-purity metallic-glass materials only show the elastic deformation or very limited plastic deformation under unconstrained conditions [12]. It is of great significance to understand the microcosmic differences between elastic strain and plastic strain. Therefore, studying structural characterization of elastic and inelastic deformation in BMGs at the atomic level is particularly important.

Recently, some researchers provide some details on the elastic deformation of amorphous solids [13–16]. For instance, Hufnagel reports the use of high-energy X-ray scattering to measure strain in a Zr₅₇Ti₅Cu₂₀Ni₈Al₁₀ bulk metallic glass in situ during uniaxial compression in the elastic regime [15]. In addition, much work has been devoted to the inhomogeneous plastic responses (shear transformation events) in MGs [17–21]. However, the exact correlation of the degree of elastic and inelastic deformation with the internal structure in MGs at the atomic level is still unresolved. This condition can be attributed to the very limited existing experimental method for detecting atomic motion. Considering the experimental difficulty of directly measuring the movement of an atom or determining which local structures are prone to form, we have focused on computer simulation of simple alloy systems. This technique eliminates some of the problems encountered in experiments and studies the deformation processes at the atomic level.

Given the lack of a structural indicator for making a distinction between elastic and inelastic deformation, understanding the microscopic deformation mechanism in MGs on the whole is difficult [22]. In this letter, the nanoscopic deformation behaviors in a ZrCu metallic-glass model during loading-unloading process under uniaxial compression have been analyzed on the basis of the molecular dynamics (MD). And we employ the reversible degree of shear origin zones (SOZs) as structural indicator to distinguish the elastic and inelastic deformation processes in MGs.

2. Theoretical Model and Calculation Details

2.1. Preparation of ZrCu MG

In our molecular dynamics (MD) simulations, the model system of ZrCu MG was supplied in a large-scale atomic/molecular massively parallel simulator (LAMMPS) [23]. The force field used in this work is a realistic embedded atom method (EAM) potential [24], which is suitable for computing pairwise interactions for metals and metal alloys. The dimensions of the model structures used in the calculations are 17.1 nm × 25.6 nm × 4.3 nm in the -, -, and -directions, respectively. The structures contained approximately 100,000 atoms. The model was heated from 300 K to 2500 K at a constant rate of 10 K ps⁻¹. To make the state of the system as natural as possible, we relaxed the liquid system for 50 ps at 2500 K within the NPT ensemble (namely, constant number, constant pressure, and constant temperature, and the temperature was controlled by Nose-Hoover thermostat method) under periodic boundary conditions. Next, we cooled the system from 2500 K down to 100 K in 100 K decrements at a constant quenching rate of 5 K ps⁻¹ at zero external pressure in the NPT ensemble.

2.2. Compression of ZrCu MG

Afterwards, compression at a constant strain rate of 1 × 10⁸ s⁻¹ along the -direction was loaded at a constantly low temperature of 50 K. Periodic boundary conditions were imposed in the - and -directions and free surfaces in the -direction to allow the occurrence of shear offset on the free surfaces. To simulate plane strain state, we maintained the pressure of the -direction at zero.

3. Results and Discussion

3.1. The Compressive Stress-Strain Curve of ZrCu MG

The compressive stress-strain curve of ZrCu MG is shown in Figure 1. The yield stress is linear when the strain is <4%, suggesting an elastic deformation stage. The simulation result is in good agreement with the result obtained by Cao [25]. With increased strain, the rise of stress values becomes slow, showing a platform corresponding to the inelastic deformation. The stress reaches a maximum of 2.3 GPa at a strain of about 6%. To clarify the microstructure evolution at different stages of the system, 2% and 5% compressive strain of the system at the imposed -direction, respectively, are analyzed.

Figure 1: Typical stress-strain curve of ZrCu MG as a function of overall sample strain.

3.2. Change in -Displacement under Relaxation Process of Elastic and Inelastic Deformation

To display the relaxation process of elastic and inelastic deformation, we used the degree of -displacement of each atom as the structural indicator, as shown in Figure 2. It is important to note that there are also -displacement and -displacement, and we just put the -displacement as a representative to characterize the relaxation process. Figures 2(a)–2(d) show the relaxation process by monitoring -displacements after unloading 2% stress. Figures 2(a′)–2(d′) reveal the relaxation process under 5% strain. In Figure 2(a), some -displacements exist after compressing the ZrCu system. However, after 150 ps relaxation (Figures 2(b)–2(d)), -displacements gradually disappear, suggesting that the -displacements of each atom are reversible at the elastic deformation stage. On the other hand, the -displacements of the system with 5% strain are strongly evident, as shown in Figure 2(a′). After 150 ps relaxation, the -displacements gradually recover but numerous -displacements still exist in the system, as shown in Figure 2(d′). Figures 2(a′)–2(d′) suggest that these -displacements are irreversible at a strain of 5%, which correspond to the inelastic deformation stage. This observation also indicates that some atoms cannot go back to the initial location at the inelastic deformation stage. The author Zhang et al. also reported that the instantaneous elastic strain is observed to vanish immediately as the applied stress is removed (unloading process), while inelastic strain that decays beyond the simulation timescale is still remained [26]. This validates our results and implies that, with the increase in strain, some atoms deviate from their original positions, causing the local structure to be destroyed. Even after a long period of relaxation, these atoms cannot be restored to their original locations. To eliminate system size effects, we also structured a bigger model whose dimensions were 27.5 nm × 5.5 nm × 55 nm in the -, -, and -directions, respectively, containing approximately 500,000 atoms. And we found that the main results presented above do not change.

Figure 2: Relaxation process observed by monitoring -displacement after unloading the stress of ZrCu MG. (a)–(d) compressed by 2% strain. (a′)–(d′) compressed by 5% strain. The red region represents the atoms moving along the positive -direction, and the blue region represents the atoms moving along the negative -direction.

3.3. Change in under Relaxation Process of Elastic and Inelastic Deformation

To further explain the phenomenon of the -displacement in the relaxation process, we use the atomic local shear strain as a structural indicator to monitor deformation processes [27]. We use relaxed glass prior to performing uniaxial compression deformation in the -direction as a reference. Atoms in Figure 3 are colored according to their atomic shear strain . Regions that have undergone larger shear strain (light blue) have a relatively large .

Figure 3: Relaxation process observed by monitoring after unloading the stress of ZrCu MG. (a)–(d) compressed by 2% strain. (a′)–(d′) compressed by 5% strain. Regions that have undergone larger shear strain (light blue) have relative large .

Shear origin zones (SOZs) are the regions in a metallic glass where strain becomes inhomogeneous [28]. A typical experimentally observed shear band width in a metallic glass is on the order of 10 nm to 20 nm [29], but the size of SOZs is about 1 nm to 2 nm. Therefore, SOZs are smaller than shear bands. In Figure 3(a), some SOZs can be clearly observed after 2% strain compression of the ZrCu systems. However, after 150 ps relaxation (Figure 3(d)), fewer SOZs are observed, suggesting that the formation of SOZs is reversible at the elastic deformation stage. In Figure 3(a′), many SOZs can be observed after 5% strain compression of the ZrCu systems. After 150 ps relaxation (Figure 3(d′)), many SOZs are still observed, suggesting that the formation of SOZs is irreversible at the inelastic deformation stage. The quantity and size of the SOZs in Figure 3(a′) are larger than those of the SOZs in Figure 3(a). This observation indicates that new SOZs form with increased strain and the origins increase in size. With increased strain, these SOZs interact with one another and develop into shear bands. Smaller amounts of SOZs can be seen in Figure 3(d), whereas numerous SOZs can be observed in Figure 3(d′). This observation implies that, with increased SOZ size to a certain value, the formation of SOZs becomes irreversible. Hence, we propose that elastic deformation is related to these SOZs in MGs. The degree of SOZs in MGs can be controlled, and the plasticity of metallic glass can be improved [28].

3.4. Change in the Fraction of FI in SOZs under Relaxation Process

To further explain the structural features of SOZs at different deformation stages of ZrCu metallic glass, the Voronoi tessellation method was used [30]. In this method, the simulation cell is divided into Voronoi polyhedra (VPs) around each atom. A VP is defined as the minimal polyhedron whose planar faces bisect at right angles the lines joining an atom to its neighboring atoms. The polyhedra can be characterized by the Voronoi index [, , , ], where denotes the number of -edged faces of the VP. An -edged face reflects the local symmetry of the central atom with some nearest neighbor atoms in a certain direction.

In previous studies, the Cu-centered Voronoi polyhedron with index (full icosahedron (FI)) has been found to be a key structural motif in amorphous Cu-Zr alloys [31]. FI units with high packing density and high shear resistance reportedly determine the plasticity of Cu-Zr glasses.

For our model ZrCu MG, we observed a clear difference between the responses of Cu atoms in Cu-centered full icosahedra [Cu-FI; Voronoi index (0, 0, 12, 0) clusters)] and those that do not have full-icosahedra surroundings (Cu-non-FI). The Voronoi tessellation analysis shows that 21.5% of Cu atoms are centers of full icosahedra in this sample. In our Voronoi analysis, the cutoff is chosen large enough (i.e., 5 Å) so that the Voronoi indices distributions are converged. Meanwhile, we also performed the analyses using a cutoff of 3.8 Å, which corresponds to the first minimum of , and found that the main results do not change.

Figure 4(a) shows the fraction of Cu-centered full icosahedra in SOZs during deformation before the strain reaches ~5%. Evidently, the fraction of Cu-centered full icosahedra is maintained at 21% before the strain reaches 4%. Lekka et al. reported that the numbers of icosahedral clusters (ICO) destroyed and created within the elastic regime are almost the same, resulting in a dynamically constant number of ICO clusters [32]. However, a sudden and rapid drop in the fraction of Cu-centered full icosahedra is observed at , corresponding to the nonlinear section of the stress-strain curve (Figure 1). This observation indicates that elastic deformation translates to inelastic deformation associated with a drop of the numbers of the full icosahedra in SOZs. Meanwhile, along with the drop of the numbers of the full icosahedra in SOZs, the efficient packing for this composition is destroyed to some degree and the energy of these zones increases. Hence, the stability of SOZs is weakened and the structures become inhomogeneous.

(a)

(b)

Figure 4: (a) Fraction of Cu-centered full icosahedra at different strains during compression deformation. (b) Fractions of Cu-centered full icosahedra as a function of relaxation time. To show whether the drop of the numbers of the Cu-centered full icosahedra can recover to the initial fraction, relaxation of the sample compressed by 0%, 2%, and 5% compressive strain, respectively, is analyzed.

To check whether the drop of the numbers of the Cu-centered full icosahedra can recover to the initial fraction, the ZrCu systems compressed by 2% and 5% strain were analyzed, as shown in Figure 4(b). With increased relaxation time, the fraction of Cu-centered polyhedra of the ZrCu systems compressed by 5% returns to ~19.3% from ~18.33%, whereas the fraction of Cu-centered polyhedra of the ZrCu systems compressed by 2% returns to ~21.5% from ~21.1%. The fraction of Cu-centered polyhedra of the ZrCu systems before being compressed is ~21.5%, as shown by the line of 0% strain in Figure 4(b). This fraction testifies the deformation along with decreasing the fractions of Cu-centered full icosahedra. Before = 4%, the decreasing process is reversible; after = 4%, the decreasing process becomes irreversible. This finding indicates that FI in SOZs is partially destroyed at the inelastic deformation stage. Interestingly, both fractions of Cu-centered polyhedra of the ZrCu systems compressed by 2% strain and 5% strain increase at the initial stage of relaxation. This phenomenon is facilitated by the additional energy introduced by the large stress overshoot caused by the high simulated strain rate. Therefore, the phenomenon leads to a recovery of the atomic structure, that is, a decrease in free volume and an increase in FI density.

4. Conclusion

In conclusion, elastic deformation is related to the reversible degree of SOZs. We find that the formation of SOZs is reversible at the elastic stage but is irreversible at the inelastic stage. At the elastic deformation stage, the fraction of FI in SOZs remains at a high grade. However, the fraction of FI in SOZs quickly decreases with increased strain at the inelastic deformation stage, and the decreasing process is also irreversible. This finding implies that the FI structure is destroyed and that the local structural environments of the SOZs are changed. Therefore, we can change the degree of SOZ localization to adjust plastic deformation and accordingly provide foundations for basic research on and future applications of MG.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors thank Shaopeng Pan from Peking University for helpful discussion. This work was supported by the National Basic Research Program of China (2013CB733000) and the National Natural Science Foundation of China (51271161 and 51271162).

References

C. A. Pampillo, “Flow and fracture in amorphous alloys,” Journal of Materials Science, vol. 10, no. 7, pp. 1194–1227, 1975.

H. S. Chen and D. E. Polk, “Mechanical properties of NiFe based alloy glasses,” Journal of Non-Crystalline Solids, vol. 15, no. 2, pp. 174–178, 1974.

K. Kondoh, K. Kawabata, T. Serikawa, and H. Kimura, “Structural characteristics and crystallization of metallic glass sputtered films by using Zr system target,” Advances in Materials Science and Engineering, vol. 2008, Article ID 312057, 4 pages, 2008.

D. Srolovitz, V. Vitek, and T. Egami, “An atomistic study of deformation of amorphous metals,” Acta Metallurgica, vol. 31, no. 2, pp. 335–352, 1983.

J. Lu, G. Ravichandran, and W. L. Johnson, “Deformation behavior of the Zr41.2Ti13.8Cu12.5Ni₁₀Be22.5 bulk metallic glass over a wide range of strain-rates and temperatures,” Acta Materialia, vol. 51, no. 12, pp. 3429–3443, 2003.

C. D. Lorenz and M. J. Stevens, “Fracture behavior of Lennard-Jones glasses,” Physical Review E—Statistical, Nonlinear, and Soft Matter Physics, vol. 68, no. 2, Article ID 021802, 6 pages, 2003.

C. A. Schuh, T. C. Hufnagel, and U. Ramamurty, “Mechanical behavior of amorphous alloys,” Acta Materialia, vol. 55, no. 12, pp. 4067–4109, 2007.

W. H. Wang, “The elastic properties, elastic models and elastic perspectives of metallic glasses,” Progress in Materials Science, vol. 57, no. 3, pp. 487–656, 2012.

G. Li, Y. Y. Wang, P. K. Liaw, Y. C. Li, and R. P. Liu, “Electronic structure inheritance and pressure-induced polyamorphism in lanthanide-based metallic glasses,” Physical Review Letters, vol. 109, no. 12, Article ID 125501, 2012.

W. L. Johnson and K. Samwer, “A universal criterion for plastic yielding of metallic glasses with a (T/Tg)2/3 temperature dependence,” Physical Review Letters, vol. 95, no. 19, Article ID 195501, 2005.

M. Chen, “Mechanical behavior of metallic glasses: microscopic understanding of strength and ductility,” Annual Review of Materials Research, vol. 38, pp. 445–469, 2008.

D. Qiao, L. Huang, G. Wang, and P. K. Liaw, “Strain rate effects on the deformation behavior of (Zr₅₅Al₁₀Ni₅Cu₃₀)₉₉Y₁ and Cu₄₅Zr47.5Ag7.5 bulk metallic glasses,” Materials Science and Engineering A, vol. 526, no. 1-2, pp. 69–73, 2009.

M. Tsamados, A. Tanguy, C. Goldenberg, and J.-L. Barrat, “Local elasticity map and plasticity in a model Lennard-Jones glass,” Physical Review E—Statistical, Nonlinear, and Soft Matter Physics, vol. 80, no. 2, Article ID 026112, 17 pages, 2009.

M. J. Demkowicz and A. S. Argon, “Autocatalytic avalanches of unit inelastic shearing events are the mechanism of plastic deformation in amorphous silicon,” Physical Review B—Condensed Matter and Materials Physics, vol. 72, no. 24, Article ID 245206, 2005.

T. C. Hufnagel, R. T. Ott, and J. Almer, “Structural aspects of elastic deformation of a metallic glass,” Physical Review B—Condensed Matter and Materials Physics, vol. 73, no. 6, Article ID 064204, 8 pages, 2006.

J. Ding, Y. Q. Cheng, and E. Ma, “Correlating local structure with inhomogeneous elastic deformation in a metallic glass,” Applied Physics Letters, vol. 101, no. 12, Article ID 121917, 2012.

M. L. Falk, J. S. Langer, and L. Pechenik, “Thermal effects in the shear-transformation-zone theory of amorphous plasticity: comparisons to metallic glass data,” Physical Review E—Statistical, Nonlinear, and Soft Matter Physics, vol. 70, no. 1, Article ID 011507, 11 pages, 2004.

V. V. Bulatov and A. S. Argon, “A stochastic model for continuum elasto-plastic behavior. III. Plasticity in ordered versus disordered solids,” Modelling and Simulation in Materials Science and Engineering, vol. 2, no. 2, pp. 203–222, 1994.

L. Qi, H. F. Zhang, Z. Q. Hu, and P. K. Liaw, “Molecular dynamic simulation studies of glass formation and atomic-level structures in Pd-Ni alloy,” Physics Letters, A: General, Atomic and Solid State Physics, vol. 327, no. 5-6, pp. 506–511, 2004.

C. A. Schuh and A. C. Lund, “Atomistic basis for the plastic yield criterion of metallic glass,” Nature Materials, vol. 2, no. 7, pp. 449–452, 2003.

H. L. Peng, M. Z. Li, and W. H. Wang, “Structural signature of plastic deformation in metallic glasses,” Physical Review Letters, vol. 106, no. 13, Article ID 135503, 4 pages, 2011.

D. Pan, A. Inoue, T. Sakurai, and M. W. Chen, “Experimental characterization of shear transformation zones for plastic flow of bulk metallic glasses,” Proceedings of the National Academy of Sciences of the United States of America, vol. 105, no. 39, pp. 14769–14772, 2008.

S. Plimpton, “Fast parallel algorithms for short-range molecular dynamics,” Journal of Computational Physics, vol. 117, no. 1, pp. 1–19, 1995.

M. I. Mendelev, D. J. Sordelet, and M. J. Kramer, “Using atomistic computer simulations to analyze x-ray diffraction data from metallic glasses,” Journal of Applied Physics, vol. 102, no. 4, Article ID 043501, 7 pages, 2007.

A. J. Cao, Y. Q. Cheng, and E. Ma, “Structural processes that initiate shear localization in metallic glass,” Acta Materialia, vol. 57, no. 17, pp. 5146–5155, 2009.

Y. Zhang, N. Mattern, and J. Eckert, “Effect of uniaxial loading on the structural anisotropy and the dynamics of atoms of Cu₅₀Zr₅₀ metallic glasses within the elastic regime studied by molecular dynamics simulation,” Acta Materialia, vol. 59, no. 11, pp. 4303–4313, 2011.

F. Shimizu, S. Ogata, and J. Li, “Theory of shear banding in metallic glasses and molecular dynamics calculations,” Materials Transactions, vol. 48, no. 11, pp. 2923–2927, 2007.

S. D. Feng, L. Qi, G. Li, W. Zhao, and R. P. Liu, “Effects of pre-introduced shear origin zones on mechanical property of ZrCu metallic glass,” Journal of Non-Crystalline Solids, vol. 373-374, pp. 1–4, 2013.

Y. F. Shi and M. L. Falk, “Atomic-scale simulations of strain localization in three-dimensional model amorphous solids,” Physical Review B—Condensed Matter and Materials Physics, vol. 73, no. 21, Article ID 214201, 10 pages, 2006.

W. Brostow, M. Chybicki, R. Laskowski, and J. Rybicki, “Voronoi polyhedra and Delaunay simplexes in the structural analysis of molecular-dynamics-simulated materials,” Physical Review B—Condensed Matter and Materials Physics, vol. 57, no. 21, pp. 13448–13458, 1998.

Y. Q. Cheng, A. J. Cao, H. W. Sheng, and E. Ma, “Local order influences initiation of plastic flow in metallic glass: effects of alloy composition and sample cooling history,” Acta Materialia, vol. 56, no. 18, pp. 5263–5275, 2008.

C. E. Lekka, A. Ibenskas, A. R. Yavari, and G. A. Evangelakis, “Tensile deformation accommodation in microscopic metallic glasses via subnanocluster reconstructions,” Applied Physics Letters, vol. 91, no. 21, Article ID 214103, 3 pages, 2007.