Next Article in Journal
Immediate Post-Extraction Short Implant Placement with Immediate Loading and without Extraction of an Impacted Maxillary Canine: Two Case Reports
Next Article in Special Issue
Rejuvenation in Deep Thermal Cycling of a Generic Model Glass: A Study of Per-Particle Energy Distribution
Previous Article in Journal
Electron-Phonon Coupling Parameter of Ferromagnetic Metal Fe and Co
Previous Article in Special Issue
Enhanced Mechanical Properties of Metallic Glass Thin Films via Modification of Structural Heterogeneity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Molecular Dynamics Study of the Nanoindentation Behavior of Cu64Zr36/Cu Amorphous/Crystalline Nanolaminate Composites

1
Department of Engineering Mechanics, School of Civil Engineering, Wuhan University, Wuhan 430072, China
2
State Key Laboratory of Water Resources & Hydropower Engineering Science, Wuhan University, Wuhan 430072, China
3
Erich Schmid Institute of Materials Science, Austrian Academy of Sciences, Jahnstraße 12, A-8700 Leoben, Austria
4
Institut für Materialwissenschaft, Technische Universität Darmstadt, Otto-Berndt-Straße 3, D-64287 Darmstadt, Germany
5
Department of Materials Science, Chair of Materials Physics, Montanuniversität Leoben, Jahnstraße 12, A-8700 Leoben, Austria
*
Author to whom correspondence should be addressed.
Materials 2021, 14(11), 2756; https://doi.org/10.3390/ma14112756
Submission received: 29 March 2021 / Revised: 16 May 2021 / Accepted: 21 May 2021 / Published: 23 May 2021
(This article belongs to the Special Issue Advances in Metallic Glass Matrix Composites)

Abstract

:
Amorphous/crystalline nanolaminate composites have aroused extensive research interest because of their high strength and good plasticity. In this paper, the nanoindentation behavior of Cu64Zr36/Cu amorphous/crystalline nanolaminates (ACNLs) is investigated by molecular dynamics (MD) simulation while giving special attention to the plastic processes occurring at the interface. The load–displacement curves of ACNLs reveal small fluctuations associated with shear transformation zone (STZ) activation in the amorphous layer, whereas larger fluctuations associated with dislocations emission occur in the crystalline layer. During loading, local STZ activation occurs and the number of STZs increases as the indentation depth in the amorphous layer increases. These STZs are mostly located around the indenter, which correlates to the high stresses concentrated around the indenter. When the indenter penetrates the crystalline layer, dislocations emit from the interface of amorphous/crystalline, and their number increases with increasing indentation depth. During unloading, the overall number of STZs and dislocations decreases, while other new STZs and dislocations become activated. These results are discussed in terms of stress distribution, residual stresses, indentation rate and indenter radius.

1. Introduction

Metallic glasses (MGs) exhibit excellent mechanical, chemical and physical performances owing to the lack of grain boundaries and crystal defects [1,2,3]. However, their limited ductility with almost no global plasticity has impeded the wide application of MGs as structural materials. To overcome this limitation, amorphous/crystalline nanolaminates consisting of nanoscale amorphous and ductile crystalline layers (ACNLs) have received attention due to their high strength and good plasticity [4,5,6].
At the nanoscale, the co-deformation of glassy and crystalline layers can effectively prevent catastrophic and localized shear banding [7,8,9]. Due to the extreme difficulties of in situ experimental observation and the microstructural complexity of amorphous/crystalline interfaces, molecular dynamics (MD) simulation is suitable for probing atomic-level structural changes and the complex deformation mechanism. The STZ activation process in the amorphous layers and dislocations nucleation and dynamics in crystalline layers have recently been revealed by many MD simulations of ACNLs [10,11,12].
Previous studies [8,13,14] have demonstrated that ACNLs have both high-strength and large plastic deformability because of the special mechanical behavior at the amorphous/crystalline interfaces. Ass a popular experimental method, nanoindentation can effectively investigate the mechanical properties of materials at the nanometer scale [15,16]. Atomistic models as a useful supplement can provide some important information for the understanding complex microstructure evolution characteristics during nanoindentation; especially they have the indisputable advantage of characterizing the dislocation movement in crystalline materials [17,18,19,20] and shear bands (SB) in single-phase MGs [21,22,23,24,25,26,27], respectively. Incorporating crystalline phases mitigates the poor failure resistance in MGs, which accommodates plastic strain by dislocation slip, effectively preventing shear band propagation [4,5,6,27,28,29]. More recent simulation studies have confirmed that the crystalline phase inserted into the amorphous phase can impede the propagation of SBs, and plays a vital role in the plastic deformation mechanism of ACNLs [30,31,32,33,34,35,36,37]. Furthermore, these studies also provide insight into how the microstructure, layer thickness and individual phase properties can be tailored to optimize the mechanical properties of the composites [35,36,37,38]. Although atomic simulation of the nanoindentation behavior and the influence of different factors on mechanical properties have been considered comprehensively, the change of the local stress distribution from the amorphous layer to the crystalline layer and the deformation behavior during the unloading process were not discussed yet. Moreover, the activation of STZs, dislocation nucleation and the change of the stress distribution during nanoindentation are seldom observed in experiments, especially in the unloading process. Therefore, the details of the activation of STZs, dislocation nucleation and stress distribution characteristics during nanoindentation loading and unloading processes still remain largely speculative, and the question of how STZs and dislocations contribute to the serrated flow of ACNLs during nanoindentation is not yet fully understood. It is also important to investigate the microstructure evolution and the change of the local stress distribution during nanoindentation, which is helpful for further understanding the plastic deformation behavior and the deformation mechanism of ACNLs during nanoindentation.
In this paper, nanoindentation tests on ACNLs were performed by using MD simulations to examine the evolution of dislocation nucleation, STZ activation and stress distribution from amorphous to crystalline layers. On an atomistic level, the activation and distribution of STZs and dislocation nucleation at maximum indentation depth and after complete unloading are discussed in detail. The aim was to investigate the serrated flow behavior of ACNLs under nanoindentation and reveal the activation and evolution processes of STZs and dislocations and the deformation mechanisms during loading and unloading. The present work provides important information for understanding the nanoindentation behavior and plastic deformation mechanisms of ACNLs from an atomistic perspective.

2. Simulation Procedure and Stress Calculation

In this work, MD simulations were carried out using an open-source LAMMPS code [38] to investigate dislocation movement, STZ activation and stress distributions for ACNLs during nanoindentation loading and unloading. The geometries of the model were shown in Figure 1. The model had a cubic orientation (i.e., X-[100], Y-[010], and Z-[001]), and the size of the box (X × Y × Z) was 250 Å × 100 Å × 250 Å. The thicknesses of the Cu64Zr36 MG and Cu crystalline layers were both 50 Å. Free boundary conditions were used on the top surface, while periodic boundary conditions (PBC) were applied in the x and z directions. The atomic positions in the last three layers (10.83 Å) at the bottom were frozen in, resembling a hard substrate. The system sizes in the directions perpendicular to the indentation direction were chosen to be large enough to avoid spurious effects of the PBC. In the present MD simulations, the modified Finnis–Sinclair-type potential for Cu–Zr binary alloys proposed by Mendelev et al. [39] were used, which had previously been successfully applied to simulate the deformation behaviors of Cu–Zr MGs and their composites [25,35,36,37,40]. The atoms in the indenter were kept fixed (i.e., the indenter was assumed to be an infinitely rigid body). At the start of the simulation, the indenter was relaxed with the conjugate gradient method at a time step of 2 fs for 20 ps to reach a minimum energy level. After relaxation, the indenter was inserted into the free surface under displacement control at an average indentation speed of v = 2.5 m/s to a maximum penetration depth of 85 Å. The whole indentation process was completed when the indenter returns to its initial position at the same velocity. To better highlight the material response upon mechanical loading, the initial temperature of the system was maintained at 50 K to avoid excess thermal activation [6,36,40]. Finally, the atomic configurations and the microstructural evolution were analyzed using the visualization tool OVITO [41], which provided details of the microstructural evolution of samples during nanoindentation.
The repulsive force exerted by a rigid indenter on the system in MD simulations is given by [42]:
F ( r ) = K ( r R ) 2
where r is the distance of an atom to the center of the indenter and R is the radius of the rigid indenter. Similar to previous studies, for the Cu–Zr MG system, the stiffness constant of the indenter is set to K = 10 eV/Å3 [24,43,44], and the force can be calculated by the above Equation (1).

3. Simulation Results and Discussion

3.1. Force–Displacement Curves

Figure 2 shows the force–displacement curve for the sample with a loading rate of 2.5 m/s and an indenter radius of 25 Å during the nanoindentation loading and unloading process. The loading force increases linearly with indentation depth after the indenter tip touches the surface of the sample. Furthermore, small fluctuations were observed in the Cu64Zr36 MG layer that were related to the activation of shear transformation zones (STZs). The first STZ activation occurs when the sample enters the plastic deformation stage. At this time, the indentation depth was about 10 Å, and the corresponding atomic shear strain was about 0.2. With increasing indentation depth, the indenter reached the interface, the loading force had a relatively larger change compared with the Cu64Zr36 MG layer due to the presence of interface, and then, the indenter penetrated the Cu crystalline layer. At the same time, the observed larger fluctuations were connected with dislocation slip. The first dislocation nucleation occurred at the “pop-in” region, as shown in Figure 2. With further increasing indentation depth, the indenter penetrated the Cu crystalline layer, while the observed larger fluctuations were connected with dislocation slip. Further penetration of the indenter overcame the resistance of dislocation movement, leading to more dislocation slip. This implies that larger stresses induced plastic rearrangements and caused work hardening [32,33]. In the loading stage, the loading curves increased until reaching a maximum depth. During the unloading process, the loading force largely decreased due to the significantly large adhesion between the indenter tip and the contact region. Moreover, in the retraction part, the adhesive phenomena occur before the force returned to zero when the indenter had no impact on the substrate. These results were consistent with the results in previous studies [34,36].
Figure 3a,b show force–displacement curves at different loading/unloading rates and different tip radii, respectively. It can be seen that a higher force value was achieved for performing simulation at a larger loading/unloading rate because the atoms did not have enough time to release energy, resulting in higher stresses for increasing loading rates [34,36]. Moreover, for different tip radii, this was explained by the stronger structural recovery with increasing tip radius caused by the slip mechanism and dislocations [36,45].

3.2. Analysis of Strain/Stress Localization and Distribution

In MGs, clusters of atoms in regions with high free volume overcome the energy barrier under the action of external shear stress and undergo a synergistic shearing motion relative to the substrate, thus forming STZs [46,47,48]. As stated previously, the plastic deformation of MGs is mainly due to STZ activation and percolation processes [49], while in metallic crystals, it is due to the slip of dislocations [19]. Upon indentation, the accumulation of localization plastic events caused by shear deformation leads to STZ proliferation under the indenter and eventually forms the shear deformation region (SDR) [36,37,50,51,52,53].
Here, the local atomic shear strain for each atom was calculated for quantifying plastic deformation at the atomic level, and atoms with a shear strain larger than 0.2 were considered to visualize the local plastic deformation within the sample [36,43,54]. For visualization, a color code of the shear strain between 0 and 0.5 was generally used [36,44,54]. Figure 4 shows the shear strain distribution of Cu64Zr36/Cu amorphous/crystalline nanolaminates (ACNLs) at different indentation depths. In the simulation, the radius of the indenter was 25 Å, and the loading rate was 2.5 m/s. All the atoms are colored according to the shear strain value during the nanoindentation loading, where the red and blue atoms represent higher and lower shear strain, respectively. The atomic shear strain distribution is shown at different indentation depths in the range of 0 to 85 Å. The zones of high shear strain or STZs were concentrated around the indenter, and the number of atoms with high shear strain value increased with the increase of the indentation depth. When the indenter was pressed into the Cu64Zr36 amorphous layer, STZs activate at “weak spots” (defined as regions with large free volume) around the contact zones, and the SDR emerged around the indenter [55]. The shear strain was further increased at an indentation depth of 35 Å, and STZs were generated at the periphery of the previous SDR, resulting in a further expansion of the SDR. As the indentation depth increased to 45 Å, the indenter was still in the amorphous layer. The SDR extended further outwards.
Similar to the distribution of the shear strain in Figure 4, the von Mises stress distribution also showed obvious localization; higher von Mises stress regions represent clusters of atoms with the plastic flow, and the STZs gather from the local areas of these plastic flow clusters. Figure 5 shows the von Mises stress distribution of the Cu64Zr36/Cu nanolaminates during indentation loading. The high von Mises stress was also concentrated around the indenter more than at other locations of the Cu64Zr36 amorphous layer. Moreover, the local stress increases by the pressure of the indenter with increasing indentation depth. When the indenter penetrated the Cu crystal layer, the high von Mises stress was gradually dispersed in the Cu64Zr36 amorphous layer, but the local stress was still concentrated around the indenter in the crystalline layer. The reason was that the surface of the substrate was directly affected by the pressure of the indenter; the local stress of the atomic zone around the indenter was higher than in other areas. In brief, the high von Mises stress concentration zones were concentrated around the indenter during nanoindentation loading. The local stress in the amorphous layer was larger than in the crystalline layer, which was consistent with findings by Goryaeva et al. [56]. This also indicated that the ACNLs with incorporated crystalline layers, which dissipated energy and accommodated plastic strain by dislocation slip, could avoid catastrophic localization, effectively impeding SB propagation.

3.3. Dislocation Analysis

In the amorphous layer, the high shear strain first appeared along the sides of the indenter because the pressure of the indenter induced deformation (see Figure 6). With further penetration of the amorphous layer, larger plastic zones formed around the indenter, but there was no nucleation of dislocations or stacking faults up to an indentation depth of 45 Å. Hence, none of these plastic zones were a potential origin of dislocation nucleation in the crystalline layer, as previously reported in the literature [32,33,34,35]. One reason for that could be that in our case, no dominant shear band forms; thus, the local stress fluctuations were too small to generate dislocation slip. When the indentation depth increased to 50 Å, the indenter penetrated the crystalline layer; dislocation emitted from the interface between the amorphous and crystalline layers, expressed through the appearance of dislocation in the crystalline layer, as shown in Figure 7. As the indentation depth increased from 55 to 85 Å, the number of slip dislocations increased. These dislocations mainly included 1 / 2 110 perfect dislocations, 1 / 6 112 Shockley partial dislocations, 1 / 6 110 stair–rod partial dislocations, and 1 / 3 100 Hirth dislocations. The plastic deformation of the crystalline layers was much smaller than that of the amorphous layers. This was consistent with the local stress distribution in the amorphous layers being much larger than in the crystalline layers. As the indentation depth increased, the plastic deformation slightly increased in the crystalline layers. The plastic deformation was significantly localized and mainly concentrated under the indenter.

3.4. Analysis of STZs and Dislocations during Loading and Unloading

Atomic-scale deformation is often characterized by the local shear invariant [27] or the non-affine squared displacement (D2min) as a structural parameter providing information about non-elastic deformations [55]. To elucidate the activation of STZs and dislocation movement during nanoindentation loading and unloading, we calculated the atomic shear strain, D2min, and the von Mises stress at maximum indentation depth and after complete unloading stages. Figure 8 shows cross-section snapshots of atomic shear strain, D2min, and von Mises stress of nanolaminates at maximum depth and complete unloading stages. The distributions of atomic shear strain and D2min reveal that high shear strain and high D2min concentration zones developed around the indenter at the maximum indentation. The high shear strain and high D2min concentration zones were still concentrated around the indenter after complete unloading. A high von Mises stress was also concentrated around the indenter at the maximum indentation stage. Moreover, the von Mises stress was dispersed in the Cu64Zr36 amorphous layer, but high local stress was still concentrated around the indenter even after complete unloading. Figure 8 quantitatively presents the variation of the percentage of STZ atoms with indentation depth during the loading and unloading stages. To quantitatively demonstrate the changes in the STZ activation during the loading and unloading process, as mentioned above, atoms with a shear strain larger than 0.2 were considered to visualize the local plastic deformation within the sample [36,42,53], which was also consistent with our previous study on the value of atomic shear strain of STZ activation during creep [57]. The percentage of STZ atoms in the total number of atoms increases slowly from 0% to 4.5% when the indentation depth increases from 0 to 30 Å, and the percentage of STZs shows a very rapid increase from 4.5% to 27.0% at indentation depths from 30 to 85 Å (maximum depth). During the unloading stage, the percentage of STZ atoms very rapidly decreases from 27.0% to 7.0%. Subsequently, the percentage of STZ shows a small increase and basically maintains at 9.0% until complete unloading, as shown in Figure 9. This indicated that some new STZs still activate in the unloading stage, and these activated STZs were mainly concentrated around the indenter both at the loading and unloading stages.
Figure 10 shows the dislocation evolution in the nanolaminates at maximum depth and complete unloading stated. It could be clearly seen that a larger number of dislocations emits from the interface in the crystalline layer; these dislocations were mainly composed of 1 / 6 112 Shockley partial dislocations and 1 / 6 110 stair–rod partial dislocations. The number of dislocations decreases, but some dislocations survive the unloading process. Figure 11 shows the variation of the dislocation density with indentation depth during loading and unloading. During loading, the dislocation density was close to zero in the amorphous layer, indicating no dislocations nucleated when the indenter was pressed into the amorphous layer. Once the indenter penetrated the crystalline layer, dislocations were emitted from the interface, and their number increases rapidly with increasing indentation depth. The dislocation density decreases rapidly during unloading, and when the indenter returns to the initial position, the dislocation density remains unchanged. This indicated that some dislocations still existed in the crystalline layer even after complete unloading. These results were consistent with those in Figure 9.

4. Conclusions

MD simulations were performed to investigate the nanoindentation behavior of Cu64Zr36/Cu amorphous/crystalline nanolaminate composites. The characteristics of plastic deformation (STZ activation and dislocations emission) of the nanolaminate composites are observed during nanoindentation loading and unloading. The main findings are as follows:
(1) The load–displacement curves of the nanolaminate composites showed small fluctuations in the amorphous layer, which are related to STZ activation, whereas larger fluctuations occur in the crystalline layer, which is related to dislocation emission. The force value increases with the increasing loading rate and the tip radius of the indenter.
(2) The local STZ activation occurred during loading, and the number of STZs increased with an increase of the indentation depth in the amorphous layer. These activated STZs were mostly concentrated around the indenter. When the indenter penetrated the crystalline layer, dislocations emitted from the amorphous/crystalline interface. They were localized around the indenter, and their number increased with increasing indentation depth. During unloading, the overall number of STZs and dislocations decreased, but also new STZs activated, and dislocations nucleated.
(3) During loading, high stresses were rather concentrated around the indenter than at other locations in the amorphous layer. When the indenter penetrated the crystalline layer, the high stresses were gradually dispersed in the amorphous layer, but the local stresses were still concentrated around the indenter in the crystalline layer. During unloading, the high stresses were dispersed in the amorphous layer, and the stressed in the amorphous layer were larger than in the crystalline layer.

Author Contributions

Conceptualization, D.Ş. and J.E.; methodology, W.-P.W.; software, W.-P.W.; validation, W.-P.W., D.Ş. and J.E.; formal analysis, W.-P.W.; investigation, W.-P.W.; resources, D.Ş.; data curation, W.-P.W.; writing—original draft preparation, W.-P.W.; writing—review and editing, J.E.; visualization, W.-P.W.; supervision, J.E.; project administration, J.E.; funding acquisition, J.E. All authors have read and agreed to the published version of the manuscript.

Funding

This research was founded by the NSFC, grant Nos. 11772236, 11472195, and the ERC Advanced Grant INTELHYB, grant No. ERC-2013-ADG-340025.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing is not applicable to this article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Schuh, C.A.; Hufnagel, T.C.; Ramamurty, U. Mechanical behavior of amorphous alloys. Acta Mater. 2007, 55, 4067–4109. [Google Scholar] [CrossRef]
  2. Eckert, J.; Das, J.; Pauly, S.; Duhamel, C. Mechanical properties of bulk metallic glasses and composites. J. Mater. Res. 2007, 22, 285–301. [Google Scholar] [CrossRef]
  3. Trexler, M.M.; Thadhani, N.N. Mechanical properties of bulk metallic glasses. Prog. Mater. Sci. 2010, 55, 759–839. [Google Scholar] [CrossRef]
  4. Wang, Y.; Li, J.; Hamza, A.V.; Barbee, T.W. Ductile crystalline-amorphous nanolaminates. Proc. Natl. Acad. Sci. USA 2007, 104, 11155–11160. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Jian, W.R.; Wang, L.; Li, B.; Yao, X.H.; Luo, S.N. Improved ductility of Cu64Zr36 metallic glass/Cu nanocomposites via phase and grain boundaries. Nanotechnology 2016, 27, 175701. [Google Scholar] [CrossRef] [PubMed]
  6. Sha, Z.-D.; Branicio, P.S.; Lee, H.P.; Tay, T.E. Strong and ductile nanolaminate composites combining metallic glasses and nanoglasses. Int. J. Plast. 2017, 90, 231–241. [Google Scholar] [CrossRef]
  7. Guo, W.; Jägle, E.A.; Choi, P.P.; Yao, J.; Kostka, A.; Schneider, J.M.; Raabe, D. Shear-induced mixing governs code formation of crystalline-amorphous nanolaminates. Phys. Rev. Lett. 2014, 113, 035501. [Google Scholar] [CrossRef] [PubMed]
  8. Cui, Y.; Abad, O.T.; Wang, F.; Huang, P.; Lu, T.-J.; Xu, K.-W.; Wang, J. Plastic Deformation Modes of CuZr/Cu Multilayers. Sci. Rep. 2016, 6, 23306. [Google Scholar] [CrossRef] [Green Version]
  9. Şopu, D.; Albe, K.; Eckert, J. Metallic glass nanolaminates with shape memory alloys. Acta Mater. 2018, 159, 344–351. [Google Scholar] [CrossRef]
  10. Cui, Y.; Shibutani, Y.; Li, S.; Huang, P.; Wang, F. Plastic deformation behaviors of amorphous-Cu50Zr50/crystalline-Cu nanolaminated structures by molecular dynamics simulations. J. Alloy. Compd. 2017, 693, 285–290. [Google Scholar] [CrossRef]
  11. Luan, Y.W.; Li, C.H.; Zhang, D.; Li, J.; Han, X.J.; Li, J.G. Plastic deformation mechanisms and size effect of Cu50Zr50/Cu amorphous/crystalline nanolaminate: A molecular dynamics study. Comp. Mater. Sci. 2017, 129, 137–146. [Google Scholar] [CrossRef]
  12. Zhou, H.; Qu, S.; Yang, W. An atomistic investigation of structural evolution in metallic glass matrix composites. Int. J. Plast. 2013, 44, 147–160. [Google Scholar] [CrossRef]
  13. Brandl, C.; Germann, T.; Misra, A. Structure and shear deformation of metallic crystalline–amorphous interfaces. Acta Mater. 2013, 61, 3600–3611. [Google Scholar] [CrossRef]
  14. Wang, Y.; Zhang, J.; Liang, X.; Wu, K.; Liu, G.; Sun, J. Size- and constituent-dependent deformation mechanisms and strain rate sensitivity in nanolaminated crystalline Cu/amorphous Cu–Zr films. Acta Mater. 2015, 95, 132–144. [Google Scholar] [CrossRef]
  15. Fischer-Cripps, A.C. Nanoindentation, 3rd ed.; Springer: New York, NY, USA, 2011. [Google Scholar]
  16. Burgess, T.; Ferry, M. Nanoindentaion of Metallic Glasses. Mater. Today 2009, 12, 24–32. [Google Scholar] [CrossRef]
  17. Li, J.; Van Vliet, K.J.; Zhu, T.; Yip, S.; Suresh, S. Atomistic mechanisms governing elastic limit and incipient plasticity in crystals. Nat. Cell Biol. 2002, 418, 307–310. [Google Scholar] [CrossRef] [PubMed]
  18. Liang, H.; Woo, C.H.; Huang, H.; Ngan, A.H.W.; Yu, T.X. Dislocation nucleation in the initial stage during nanoindentation. Philos. Mag. 2003, 83, 3609–3622. [Google Scholar] [CrossRef]
  19. Lilleodden, E.; Zimmerman, J.; Foiles, S.; Nix, W. Atomistic simulations of elastic deformation and dislocation nucleation during nanoindentation. J. Mech. Phys. Solids 2003, 51, 901–920. [Google Scholar] [CrossRef]
  20. Fu, T.; Peng, X.; Chen, X.; Weng, S.; Hu, N.; Li, Q.; Wang, Z. Molecular dynamics simulation of nanoindentation on Cu/Ni nanotwinned multilayer films using a spherical indenter. Sci. Rep. 2016, 6, 35665. [Google Scholar] [CrossRef] [Green Version]
  21. Shi, Y.; Falk, M.L. Stress-induced structural transformation and shear banding during simulated nanoindentation of a metallic glass. Acta Mater. 2007, 55, 4317–4324. [Google Scholar] [CrossRef]
  22. Jiang, S.Y.; Jiang, M.Q.; Dai, L.H.; Yao, Y.G. Atomistic Origin of Rate-Dependent Serrated Plastic Flow in Metallic Glasses. Nanoscale Res. Lett. 2008, 3, 524–529. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Qiu, C.; Zhu, P.; Fang, F.; Yuan, D.; Shen, X. Study of nanoindentation behavior of amorphous alloy using molecular dynamics. Appl. Surf. Sci. 2014, 305, 101–110. [Google Scholar] [CrossRef]
  24. Avila, K.E.; Küchemann, S.; Alhafez, I.A.; Urbassek, H.M. Shear-Transformation Zone Activation during Loading and Unloading in Nanoindentation of Metallic Glasses. Materials 2019, 12, 1477. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Kumar, A.; Gupta, P.; Yedla, N. Nanoindentation studies of Zr50Cu50 metallicglass thin film nanocomposites via moleculardynamics simulations. Metall. Res. Technol. 2016, 113, 602. [Google Scholar] [CrossRef]
  26. Wu, C.-D. Atomistic simulation of nanoformed metallic glass. Appl. Surf. Sci. 2015, 343, 153–159. [Google Scholar] [CrossRef]
  27. Shimizu, F.; Ogata, S.; Li, J. Theory of Shear Banding in Metallic Glasses and Molecular Dynamics Calculations. Mater. Trans. 2007, 48, 2923–2927. [Google Scholar] [CrossRef] [Green Version]
  28. Guo, W.; Jägle, E.; Yao, J.; Maier-Kiener, V.; Korte-Kerzel, S.; Schneider, J.M.; Raabe, D. Intrinsic and extrinsic size effects in the deformation of amorphous CuZr/nanocrystalline Cu nanolaminates. Acta Mater. 2014, 80, 94–106. [Google Scholar] [CrossRef]
  29. Şopu, D.; Yuan, X.; Moitzi, F.; Stoica, M.; Eckert, J. Structure-Property Relationships in Shape Memory Metallic Glass Composites. Materials 2019, 12, 1419. [Google Scholar] [CrossRef] [Green Version]
  30. Wang, C.-H.; Chao, K.-C.; Fang, T.-H.; Stachiv, I.; Hsieh, S.-F. Investigations of the mechanical properties of nanoimprinted amorphous Ni–Zr alloys utilizing the molecular dynamics simulation. J. Alloy. Compd. 2016, 659, 224–231. [Google Scholar] [CrossRef]
  31. Narayan, R.L.; Boopathy, K.; Sen, I.; Hofmann, D.; Ramamurty, U. On the hardness and elastic modulus of bulk metallic glass matrix composites. Scr. Mater. 2010, 63, 768–771. [Google Scholar] [CrossRef]
  32. Doan, D.-Q.; Fang, T.-H.; Chen, T.-H. Nanotribological characteristics and strain hardening of amorphous Cu64Zr36/ crystalline Cu nanolaminates. Tribol. Int. 2020, 147, 106275. [Google Scholar] [CrossRef]
  33. Song, H.; Yin, P.; Zuo, X.; An, M.; Li, Y. Atomic simulations of plastic deformation mechanism of MgAl/Mg nanoscale amorphous/crystalline multilayers. J. Non-Cryst. Solids 2018, 500, 121–128. [Google Scholar] [CrossRef]
  34. Jian, W.; Wang, L.; Yao, X.; Luo, S. Tensile and nanoindentation deformation of amorphous/crystalline nanolaminates: Effects of layer thickness and interface type. Comput. Mater. Sci. 2018, 154, 225–233. [Google Scholar] [CrossRef]
  35. Song, H.; Xu, J.; Zhang, Y.; Li, S.; Wang, D.; Li, Y. Molecular dynamics study of deformation behavior of crystalline Cu/amorphous Cu50Zr50 nanolaminates. Mater. Des. 2017, 127, 173–182. [Google Scholar] [CrossRef]
  36. Hua, D.; Ye, W.; Jia, Q.; Zhou, Q.; Xia, Q.; Shi, J.; Deng, Y.; Wang, H. Molecular dynamics simulation of nanoindentation on amorphous/amorphous nanolaminates. Appl. Surf. Sci. 2020, 511, 145545. [Google Scholar] [CrossRef]
  37. Avila, K.E.; Vardanyan, V.H.; Küchemann, S.; Urbassek, H.M. Response of an amorphous/crystalline interface to nanoindentation: An atomistic study. Appl. Surf. Sci. 2021, 551, 149285. [Google Scholar] [CrossRef]
  38. Plimpton, S. Fast Parallel Algorithms for Short-Range Molecular Dynamics. J. Comput. Phys. 1995, 117, 1–19. Available online: http://lammps.sandia.gov/ (accessed on 10 March 2021). [CrossRef] [Green Version]
  39. Mendelev, M.; Sordelet, D.J.; Kramer, M.J. Using atomistic computer simulations to analyze x-ray diffraction data from metallic glasses. J. Appl. Phys. 2007, 102, 043501. [Google Scholar] [CrossRef]
  40. Şopu, D.; Stoica, M.; Eckert, J. Deformation behavior of metallic glass composites reinforced with shape memory nanowires studied via molecular dynamics simulations. Appl. Phys. Lett. 2015, 106, 211902. [Google Scholar] [CrossRef] [Green Version]
  41. Stukowski, A. Visualization and analysis of atomistic simulation data with OVITO—The Open Visualization Tool. Model. Simul. Mater. Sci. Eng. 2009, 18, 015012. [Google Scholar] [CrossRef]
  42. Ruestes, C.; Bringa, E.; Gao, Y.; Urbassek, H. Molecular Dynamics Modeling of Nanoindentation. In Applied Nanoindentation in Advanced Materials; Tiwari, A., Natarajan, S., Eds.; Wiley: Chichester, UK, 2017; Chapter 14; pp. 313–345. [Google Scholar] [CrossRef]
  43. Wu, W.-P.; Şopu, D.; Yuan, X.; Eckert, J. Aspect ratio-dependent nanoindentation behavior of Cu64Zr36 metallic glass nanopillars investigated by molecular dynamics simulations. J. Appl. Phys. 2020, 128, 084303. [Google Scholar] [CrossRef]
  44. Avila, K.E.; Küchemann, S.; Urbassek, H.M. Stucture and size of the plastic zone formed during nanoindentation of a metallic glass. J. Non-Cryst. Solids 2019, 523, 119593. [Google Scholar] [CrossRef] [Green Version]
  45. Doan, D.-Q.; Fang, T.-H.; Tran, A.-S.; Chen, T.-H. Residual stress and elastic recovery of imprinted Cu-Zr metallic glass films using molecular dynamic simulation. Comput. Mater. Sci. 2019, 170, 109162. [Google Scholar] [CrossRef]
  46. Greer, A.; Cheng, Y.; Ma, E. Shear bands in metallic glasses. Mater. Sci. Eng. R Rep. 2013, 74, 71–132. [Google Scholar] [CrossRef]
  47. Argon, A. Plastic deformation in metallic glasses. Acta Met. 1979, 27, 47–58. [Google Scholar] [CrossRef]
  48. Yang, B.; Liu, C.T.; Nieh, T.G. Unified equation for the strength of bulk metallic glasses. Appl. Phys. Lett. 2006, 88, 221911. [Google Scholar] [CrossRef]
  49. Li, W.; Wei, B.; Zhang, T.; Xing, D.; Zhang, L.; Wang, Y. Study of serrated flow and plastic deformation in metallic glasses through instrumented indentation. Intermetallics 2007, 15, 706–710. [Google Scholar] [CrossRef]
  50. Şopu, D.; Stukowski, A.; Stoica, M.; Scudino, S. Atomic-Level Processes of Shear Band Nucleation in Metallic Glasses. Phys. Rev. Lett. 2017, 119, 195503. [Google Scholar] [CrossRef]
  51. Cheng, L.; Jiao, Z.M.; Ma, S.G.; Qiao, J.W.; Wang, Z.H. Serrated flow behaviors of a Zr-based bulk metallic glass by nanoindentation. J. Appl. Phys. 2014, 115, 084907. [Google Scholar] [CrossRef]
  52. Limbach, R.; Kosiba, K.; Pauly, S.; Kühn, U.; Wondraczek, L. Serrated flow of CuZr-based bulk metallic glasses probed by nanoindentation: Role of the activation barrier, size and distribution of shear transformation zones. J. Non-Cryst. Solids 2017, 459, 130–141. [Google Scholar] [CrossRef]
  53. Zhong, C.; Zhang, H.; Cao, Q.; Wang, X.; Zhang, D.; Ramamurty, U.; Jiang, J. On the critical thickness for non-localized to localized plastic flow transition in metallic glasses: A molecular dynamics study. Scr. Mater. 2016, 114, 93–97. [Google Scholar] [CrossRef]
  54. Wang, Y.-C.; Wu, C.-Y.; Chu, J.P.; Liaw, P.K. Indentation Behavior of Zr-Based Metallic-Glass Films via Molecular-Dynamics Simulations. Met. Mater. Trans. A 2010, 41, 3010–3017. [Google Scholar] [CrossRef]
  55. Falk, M.L.; Langer, J.S. Dynamics of viscoplastic deformation in amorphous solids. Phys. Rev. E 1998, 57, 7192–7205. [Google Scholar] [CrossRef] [Green Version]
  56. Goryaeva, A.M.; Fusco, C.; Bugnet, M.; Amodeo, J. Influence of an amorphous surface layer on the mechanical properties of metallic nanoparticles under compression. Phys. Rev. Mater. 2019, 3, 033606. [Google Scholar] [CrossRef] [Green Version]
  57. Wu, W.-P.; Şopu, D.; Yuan, X.; Adjaoud, O.; Song, K.; Eckert, J. Atomistic understanding of creep and relaxation mechanisms of Cu64Zr36 metallic glass at different temperatures and stress levels. J. Non–Cryst. Solids 2021, 559, 120676. [Google Scholar] [CrossRef]
Figure 1. The geometries of the model.
Figure 1. The geometries of the model.
Materials 14 02756 g001
Figure 2. Force–displacement curve during nanoindentation loading and unloading for a loading rate of 2.5 m/s and an indenter radius of 25 Å.
Figure 2. Force–displacement curve during nanoindentation loading and unloading for a loading rate of 2.5 m/s and an indenter radius of 25 Å.
Materials 14 02756 g002
Figure 3. Force–displacement curves of ACNLs at different conditions. (a) Different loading and unloading rates, (b) different indenter radii.
Figure 3. Force–displacement curves of ACNLs at different conditions. (a) Different loading and unloading rates, (b) different indenter radii.
Materials 14 02756 g003aMaterials 14 02756 g003b
Figure 4. Cross-section snapshots of atomic shear strain distributions of ACNLs during nanoindentation loading.
Figure 4. Cross-section snapshots of atomic shear strain distributions of ACNLs during nanoindentation loading.
Materials 14 02756 g004
Figure 5. Cross-section snapshots of von Mises stress distributions of ACNLs during nanoindentation loading.
Figure 5. Cross-section snapshots of von Mises stress distributions of ACNLs during nanoindentation loading.
Materials 14 02756 g005
Figure 6. The process of dislocation motion ofACNLs during loading.
Figure 6. The process of dislocation motion ofACNLs during loading.
Materials 14 02756 g006
Figure 7. The dislocation emits from the interface of ACNLsat indentation depth of 50 Å.
Figure 7. The dislocation emits from the interface of ACNLsat indentation depth of 50 Å.
Materials 14 02756 g007
Figure 8. Cross-section snapshots of atomic shear strain, D2min, and von Mises stress of ACNLs at maximum depth and complete unloading stages. (a,c,e) are atomic shear strain, D2min, and von Mises stress of ACNLs at maximum depth stage, respectively; (b,d,f) are atomic shear strain, D2min, and von Mises stress of ACNLs at complete unloading stage, respectively.
Figure 8. Cross-section snapshots of atomic shear strain, D2min, and von Mises stress of ACNLs at maximum depth and complete unloading stages. (a,c,e) are atomic shear strain, D2min, and von Mises stress of ACNLs at maximum depth stage, respectively; (b,d,f) are atomic shear strain, D2min, and von Mises stress of ACNLs at complete unloading stage, respectively.
Materials 14 02756 g008
Figure 9. Variation of the percentage of activated STZ atoms (those atoms with shear strains > 0.2) with indentation depth during loading and unloading.
Figure 9. Variation of the percentage of activated STZ atoms (those atoms with shear strains > 0.2) with indentation depth during loading and unloading.
Materials 14 02756 g009
Figure 10. The process of dislocation motion at maximum depth and complete unloading.
Figure 10. The process of dislocation motion at maximum depth and complete unloading.
Materials 14 02756 g010
Figure 11. Variation of the dislocation density with indentation depth during loading and unloading.
Figure 11. Variation of the dislocation density with indentation depth during loading and unloading.
Materials 14 02756 g011
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wu, W.-P.; Şopu, D.; Eckert, J. Molecular Dynamics Study of the Nanoindentation Behavior of Cu64Zr36/Cu Amorphous/Crystalline Nanolaminate Composites. Materials 2021, 14, 2756. https://doi.org/10.3390/ma14112756

AMA Style

Wu W-P, Şopu D, Eckert J. Molecular Dynamics Study of the Nanoindentation Behavior of Cu64Zr36/Cu Amorphous/Crystalline Nanolaminate Composites. Materials. 2021; 14(11):2756. https://doi.org/10.3390/ma14112756

Chicago/Turabian Style

Wu, Wen-Ping, Daniel Şopu, and Jürgen Eckert. 2021. "Molecular Dynamics Study of the Nanoindentation Behavior of Cu64Zr36/Cu Amorphous/Crystalline Nanolaminate Composites" Materials 14, no. 11: 2756. https://doi.org/10.3390/ma14112756

APA Style

Wu, W. -P., Şopu, D., & Eckert, J. (2021). Molecular Dynamics Study of the Nanoindentation Behavior of Cu64Zr36/Cu Amorphous/Crystalline Nanolaminate Composites. Materials, 14(11), 2756. https://doi.org/10.3390/ma14112756

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop