Abstract
Nanoindentation simulations on a binary metallic glass were performed under various strain rates by using molecular dynamics. The ratedependent serrated plastic flow was clearly observed, and the spatiotemporal behavior of its underlying irreversible atomic rearrangement was probed. Our findings clearly validate that the serration is a temporally inhomogeneous characteristic of such rearrangements and not directly dependent on the resultant shearbanding spatiality. The unique spatiotemporal distribution of shear banding during nanoindentation is highlighted in terms of the potential energy landscape (PEL) theory.
Keywords:
Metallic glasses; Serrated plastic flow; Strain rate; Molecular dynamicsNano Express
Bulk metallic glasses (BMGs), because of their longrange atomic disorder, deform uniquely [14]: the plastic deformation is highly localized into narrow shear bands at room temperature [58]. Under deformationconstrained loading modes such as compression [911] and nanoindentation [1219], serrated plastic flow phenomena have been widely observed and found to be ratedependent: as strain rate decreases, the flow serrations become more distinct. Moreover, ratedependent shearband patterns were observed on the surfaces or inside the postdeformed specimens under compression [10] and nanoindentation [1517]. Therefore, currently, it is well accepted that the macroscopic serrated plastic flow behavior is associated with the shearbanding operations on a nanoscale within BMGs. Considerable efforts have been made to uncover the relationship between them. Schuh et al. [13], Schuh and Nieh [14] and Zhang et al. [18] suggested that the simultaneous operations of multiple shear bands at high strain rate result in the smooth plastic flow; with increasing strain rate, the deformation mode transits from inhomogeneous to homogeneous. However, Jiang and Atzmon [15], Greer et al. [17], and Jang et al. [19] considered that the disappearance of flow serrations at high strain rate arises from the limitation of data acquisition and the instrumental response; even at extremely high strain rate, the plastic deformation is still inhomogeneous. Since these hypotheses were mainly gained from ex situ experimental observations on shearband patterns, the physical origin of ratedependent serrated plastic flow, even if well accepted as a temporal event or process behavior [10,11,15], is still not well grounded.
Very recently, infrared camera technique has been used for in situ observing dynamic shear banding operations in the study of the serrated plastic flow during compression by Jiang et al. [10]. Based on the information from experimental observations, they conjectured a spatiotemporal picture of shearbanding during serrated flow, and further explained the inhomogeneous deformation during the indentation [10]. However, as they pointed out the propagation of a shear band is very fast; hence only mature shear bands can be captured [10], not to mention the fine shearbanding events on atomic scale. In addition, during indentation, since the shear bands develop underneath indenter, direct observation in situ on them is very difficult, even impossible. The molecular dynamics (MD) simulation is generally believed to be an effective way in modeling various indentation processes [2022], providing an in situ observation on atomic motions. In this aspect, remarkable progress has been made by Falk, Langer, Shi et al. [2023] by developing Argon’s shear transformation zone (STZ) concept [2]. They found that the loaddrop events in a loaddisplacement curve, i.e., serrated plastic flow, correspond closely to the bursts in deformation activity (irreversible atomic rearrangement) associated with shear bands [21]. Spatially, the suppression of winglike shear bands of postdeformed specimen at higher strain rate leads to milder serrations [22]. However, how do the strain rates affect this inherent correlation between serrations and atomic rearrangements? How does the temporal behavior of such rearrangements produce the final shearbanding patterns after loading? These questions are still not well understood and deserve further investigation. In this letter, we rely on MD computer simulations to roundly probe how spatiotemporal distribution of shearbanding events on an atomic scale is related in situ to macroscopic ratedependent serrated plastic flows in BMGs undergoing nanoindentation. Strain rate effect on this relationship and its underlying physics are discussed.
A binary amorphous alloy system, Cu_{46}Zr_{54}, was used in our MD simulations. In this system, atoms interact via a modified LennardJones 4–8 potential of the form [24,25]:
where r_{ij} is the distance between the atoms i and j, A, B, C, and D are constants whose values are available in Ref. [25], and r_{t} is the truncation distance with the values of 5.08, 5.58, and 6.00 Å for Cu–Cu, Cu–Zr, and Zr–Zr pairs, respectively. The motion of each atom was evaluated by integrating the Newtonian equations of motion using velocityVerlet method with a time step of 1 fs. To form an amorphous sample, an initial structure for the sample was built by placing all atoms into a facecentered cubic (fcc) crystal lattice in a random order, and the initial velocities of all atoms were set to be zero [26]. The initial structure was gradually heated to 2400 K for sufficiently melting, and then cooled to 1 K with the cooling rate of 25 K/ps. In this process, the NPT ensemble was used, and the pressure was kept at zero; periodic boundary conditions were used in all three directions. Then, for the subsequent indentation simulations, we set the top boundary free and fixed a layer of 6.0 Å in thickness at the bottom. Another 100 ps was carried out to a new equilibration. Finally, a threedimensional sample (sample I) which contains 432,000 atoms with the size of 250 × 250 × 125 Å^{3} and a twodimensional sample (sample II) which includes 250,000 atoms with the size of 1950 × 1050 Å^{2} were formed in this way.
A spherical indenter whose atomic nature is ignored was used in the nanoindentation simulations. The indenter was modeled by a purely repulsive potential with the form [27,28]:
where r is the distance from the indenter center to a sample atom, R is the radius of the indenter which is chosen as 20 Å for sample I and 400 Å for sample II, θ(R − r) is the standard step function, and E is a constant which is equal to 3.0 and 3.9 nN/Å^{2} for a Cu atom and a Zr atom, respectively [28]. The indenter was displaced toward the top surface of the sample at a constant strain rate by keeping an invariable displacement interval of 0.1 Å and adjusting the relaxation time for each displacement interval. In this process, the control to the temperature (1 K) was only allowed in a layer of 10.0 Å in thickness which is just above the fixed layer at the bottom. For sample I, the total indentation depth was 15 Å, and three strain rates, 10^{11}, 10^{10}, and 10^{9} s^{−1}, were executed; for sample II, the total indentation depth was 50 Å, and three strain rates, 4 × 10^{10}, 4 × 10^{9}, and 4 × 10^{8} s^{−1}, were performed. Parallel computing was used in all the simulation processes.
The load–displacement (p–h) curves for the indentation simulations are shown in Fig. 1. Obviously, the ratedependent serrated plastic flow phenomena can be observed: when strain rate decreases from 1a–c, flow serrations become more prominent. The result is consistent with a series of experimental observations for real metallic glasses under indentation [1118]. It has been recognized that the serrated plastic flow, relating to shearbanding operations, occurs as a result of a number of structural rearrangements at atomic scale. The parameter D_{min}^{2}(−Δtt), therefore, is adopted to identify such irreversible rearrangement with the form [23,29]:
where the subscript n runs over the atoms within the interaction range of the reference atom (n = 0) and r_{n}(t) is the position vector of the nth atom at time t. The parameter D_{min}^{2}(t − Δtt) then denotes the local deviation between the true deformation denoted by r_{n}(t) − r_{0}(t)] and the affine deformation indicated by (X Y^{−1}) r_{n}(t − Δt)−r_{0}(t − Δt)] during the time interval t − Δtt[29]. We calculated D_{min} values of all atoms during each displacement interval (0.1 Å) to get information of in situ deformation. We selected 1.5 Å, which is about half of the average distance between a Cu atom and a Zr atom in the samples, as a cutoff of D_{min} to characterize the rearrangements that make up a plastic event at all strain rates [23,29]. It is important to point out that the method of choosing the cutoff may be a little rough, considering its value may be affected by strain rates; nevertheless, it is efficient to judge the plastic deformation, and should not significantly change the trend. In addition, we find that choosing different cutoffs cannot change the trend of plastic flow under various strain rates. Any atom whose D_{min} value is greater than the cutoff is considered to be rearranged, and the numbers of the rearranged atoms at all intervals are displayed in Fig. 1 as histograms under p–h curves. Note that the numbers of the rearranged atoms get larger, but their distribution becomes more inhomogeneous when strain rate decreases. Moreover, when comparing the histograms with the p–h curves, a strong correlation between them was surprisingly discovered: the loaddrop events (i.e., flow serrations) in the p–h curves correspond to the peak values (i.e., large numbers of the rearranged atoms) in the histograms; the more obvious the flow serration is, the larger the number of the rearranged atoms in that interval is. The phenomenon is consistent with the simulation results presented by Shi and Falk [21]. In fact, the number of the rearranged atoms can be regarded as an indication of the degree of plastic deformation in the displacement interval. Thus, we can conclude that the serrated plastic flow strongly depends on the temporal characteristic of the atomic rearrangement underpinning plastic deformation: successive low degree of plastic deformation at high strain rate leads to less pronounced serrated flow, and intermittent high degree of plastic deformation at low strain rate results in more distinct serrated flow.
Figure 1. The load–displacement curves and the temporal distribution of the number of rearranged atoms at various strain rates for sample I. The strain rate decreases from (a–c)
To ferret out how shear bands operate during the indentation processes, we examined the spatial distributions of the rearranged atoms when flow serrations occur. As shown in Fig. 2, we found that, during a displacement interval, at high strain rate, few rearranged atoms form many small atomic clusters, and at low strain rate, many rearranged atoms form few large atomic clusters. The figure indicates that the degree of atomic rearrangement underlying flow in an individual loading step varies with strain rate or loading timescale. This kind of discrete flow event (i.e., rearranged atomic cluster) finally leads to distinguishing shearband patterns at the maximum indenting depth, which are displayed in Fig. 3: more and thinner shear bands formed at high strain rate, while fewer and coarser shear bands nucleated at low strain rate. The patterns are drawn by coloring the atoms according to their D_{min}(0, t) values; here, the darker the color the larger the D_{min}(0, t) value. They are quite similar to those observed from instrumented [1519] and simulated [20,22] indentations. The results in Figs. 2 and 3 were taken from sample II for the reason that larger planar size is available to display the shearband patterns. Furthermore, the potential energy versus displacement curves, as shown in Fig. 4a, were also found to be ratedependent. With decreasing strain rate, the potential energy grows more slowly, but fluctuates more prominently; the energydrop events essentially correspond to the loaddrop events. Since our MD processes proceed at very low temperature, the potential energy which is large compared to the thermal energy must dominate the flow [30]. Thus, the unique spatiotemporal characteristic of deformation can be understood in terms of potential energy landscape (PEL) theory [3033].
Figure 2. The ratedependent spatial distributions of the rearranged atoms in the displacement intervals where flow serrations occur for sample II. The strain rate decreases from (a–c);astrain rate 4 × 10^{10} s^{−1},bstrain rate 4 × 10^{9} s^{−1},cstrain rate 4 × 10^{8} s^{−1}
Figure 3. The ratedependent shearband patterns of the maximum indenting depth for sample II. The strain rate decreases from (a–c);astrain rate 4 × 10^{10} s^{−1},bstrain rate 4 × 10^{9} s^{−1},cstrain rate 4 × 10^{8} s^{−1}
Figure 4. (a) The potential energy versus displacement curves at different strain rates for sample I; (b) the schematic of the disappearance of a local energy minimum induced by loading; (c) the schematic of the ratedependent systematic energy states when the system transits to a new energy basin
Before loading, the system stays at a local minimum of the potential energy surface (PES), and the configuration of the system is metastable. Loading will tilt the PES, and induce the disappearance of some local minima. As a result, the system will became unstable, and move to a new energy local basin. At the same time, atomic rearrangement occurs. The disappearance of a potential energy basin induced by loading is schematically shown in Fig. 4b. When transiting from one energy basin to another, the system at high strain rate will stay at a higher energy state in the new basin in a displacement interval than at low strain rate because of the shorter relaxation time it has. The comparison of the systematic states at different strain rates is schematically shown in Fig. 4c. Thus a whole transition of the system at high strain rate usually costs more displacement intervals, and the degree of plastic deformation at a particular displacement is relatively lower. On the contrary, at low strain rate, the system experiences less displacement intervals to reach the new basin, and the degree of plastic deformation at a particular displacement is generally higher. Hence the ratedependent temporal distribution of deformation appears. It is this temporal characteristic of atomic rearrangement that dominates the macroscopic serrated plastic flow. On the other hand, since the potential energy of the system at the same indenting depth at high strain rate is higher, there will be more atoms at high energy state in the plastic zone. These atoms do not have enough time to rearrange. Thus, the rearrangement occurs at multiple regions, but the number of the rearranged atoms is small (see Fig. 2a). When next loading step is applied, atoms preferentially rearrange at the same positions owing to their relatively higher local energy. In other words, the shear banding preferentially operates at the same locations at higher strain rate [10]. So it is hard to conjecture from the shearband patterns how many shear bands operate in a moment. On the other hand, the system has lower potential energy at low strain rate, so the number of the atoms with high energy state will be smaller. Thus the rearrangement happens at fewer regions. Nevertheless, the rearrangement in these few regions can develop more sufficiently, finally producing coarse shear bands (see Fig. 3c). In addition, since local atoms rearrange sufficiently, leading to a lower local energy, the atomic rearrangement at the next step must occur at other regions having high energy level. In other words, the shear bands nucleate and develop sufficiently at different positions under low strain rate. If the plastic deformation zone or spatial distribution is not related to its time, such as the uniaxial compression case [10], more and finer shear bands can finally be produced at low strain rates. However, in nanoindentation, the unique spatial patterning, fewer and coarser shear bands (see Fig. 3c), can be observed owing to space–time relationship of plastic deformation.
In summary, ratedependent serrated plastic flow was observed in our MD simulations. The temporally inhomogeneous characteristic of the plastic deformation was revealed as the main determining factor of serrated flow behavior. It is not proper to directly relate the numbers of shear bands with the flow serrations. For example, although there is no serration in some Febased or Cebased BMGs during nanoindentation, a number of fine shear bands are observed under the indents [34]. The unique ratedependent spatiotemporal distributions of shear banding can be understood in terms of PEL theory. We believe that these findings can shed light on the relationship between microstructure and inhomogeneous plastic flow in BMGs.
Acknowledgements
The work was supported by the Natural Science Foundation of China (Grants Nos. 10725211, 10721202, 10534030, 10674163), the Ministry of Science and Technology of China (2006CB921300, 2007CB925000), and the Knowledge Innovation Project & Key Project of Chinese Academy of Sciences (Nos. KJCX2YWM04 and KJCXSWL08). All computation of this work were carried out by Supercomputer DeepComp 6800, and we thank Dr. Yangde Feng of Super Computing Center of Chinese Academy of Science for his help in the computations.
References

Acta Met.. 1977, 25:407.
COI number [1:CAS:528:DyaE2sXltVGntrg%3D]
Publisher Full Text 
Acta Metal.. 1979, 27:47.
COI number [1:CAS:528:DyaE1MXhvFSqtrg%3D]
Publisher Full Text 
Schun CA, Hufnagel TC, Ramamurty U:
Acta Mater.. 2007, 55:4067. Publisher Full Text

Science. 2007, 318:1880.
COI number [1:CAS:528:DC%2BD2sXhsVygtb%2FO]
Publisher Full Text 
Hufnagel TC, EIDeiry P, Vinci RP:
Scr. Mater.. 2000, 43:1071.
COI number [1:CAS:528:DC%2BD3cXoslSls7w%3D]
Publisher Full Text 
Philos. Mag.. 2002, 82:2623.
COI number [1:CAS:528:DC%2BD38XntFGlsL0%3D]

Liu LF, Dai LH, Bai YL, Wei BC:
J. NonCryst. Solids. 2005, 351:3259.
COI number [1:CAS:528:DC%2BD2MXhtVKjtrrN]
Publisher Full Text 
Appl. Phys. Lett.. 2006, 89:071907. Publisher Full Text

Liu LF, Dai LH, Bai YL, Wei BC, Eckert J:
Mater. Chem. Phys.. 2005, 93:174.
COI number [1:CAS:528:DC%2BD2MXlt1Gisb4%3D]
Publisher Full Text 
Jiang WH, Fan GJ, Liu FX, Wang GY, Choo H, Liaw PK:
Int. J. Plast.. 2007, 24:1.
COI number [1:CAS:528:DC%2BD1cXitlOltrs%3D]
Publisher Full Text 
Mater. Trans JIM. 2001, 42:642.
COI number [1:CAS:528:DC%2BD3MXktlKjsLc%3D]
Publisher Full Text 
Acta Mater.. 2007, 55:295.
COI number [1:CAS:528:DC%2BD28Xht1KhsrnE]
Publisher Full Text 
Schuh CA, Argon AS, Nieh TG, Wadsworth J:
Philos. Mag.. 2003, 83:2585.
COI number [1:CAS:528:DC%2BD3sXmtVGktbk%3D]
Publisher Full Text 
Acta Mater.. 2003, 51:87.
COI number [1:CAS:528:DC%2BD38Xps1OgtLk%3D]
Publisher Full Text 
J. Mater. Res.. 2003, 18:755.
COI number [1:CAS:528:DC%2BD3sXivVeltr0%3D]
Publisher Full Text 
Dai LH, Liu LF, Yan M, Wei BC, Eckert J:
Chin. Phys. Lett.. 2004, 21:1593.
COI number [1:CAS:528:DC%2BD2MXltlKlsbw%3D]
Publisher Full Text 
Greer AL, Castellero A, Madge SV, Walker IT, Wilde JR:
Mater. Sci. Eng. A. 2004, 375:1182. Publisher Full Text

Zhang GP, Wang W, Zhang B, Tan J, Liu CS:
Scr. Mater.. 2005, 52:1147.
COI number [1:CAS:528:DC%2BD2MXjtV2nsrs%3D]
Publisher Full Text 
Appl. Phys. Lett.. 2007, 90:211906. Publisher Full Text

Appl. Phys. Lett.. 2005, 86:011914. Publisher Full Text

Acta Mater.. 2007, 55:4317.
COI number [1:CAS:528:DC%2BD2sXntFensbg%3D]
Publisher Full Text 
Phys. Rev. E Stat. Phys. Plasmas Fluids Relat. Interdiscip. Topics. 1998, 57:7192.
COI number [1:CAS:528:DyaK1cXjvVKmsL4%3D]

Kobayashi S, Maeda K, Takeuchi S:
J. Phys. Soc. Jpn.. 1980, 48:1147.
COI number [1:CAS:528:DyaL3cXit1CntLg%3D]
Publisher Full Text 
Acta Mater.. 2003, 51:5399.
COI number [1:CAS:528:DC%2BD3sXnsl2rt7s%3D]
Publisher Full Text 
Phys. Rev. B. 2004, 69:144205. Publisher Full Text

Kelchner CL, Plimpton SJ, Hamilton JC:
Phys. Rev. B. 1998, 58:11085.
COI number [1:CAS:528:DyaK1cXntVaqsbg%3D]
Publisher Full Text 
Nanotechnology. 2003, 14:1208.
COI number [1:CAS:528:DC%2BD2cXhvFGktQ%3D%3D]
Publisher Full Text 
Gannepalli A, Mallapragada SK:
Nanotechnology. 2001, 12:250.
COI number [1:CAS:528:DC%2BD3MXotlGls7Y%3D]
Publisher Full Text 
J. Chem. Phys.. 1969, 51:3728.
COI number [1:CAS:528:DyaE3cXmsVSl]
Publisher Full Text 
Science. 1984, 225:983.
COI number [1:CAS:528:DyaL2cXlsFyisbs%3D]
Publisher Full Text 
J. Phys. Chem. B. 2004, 108:19619.
COI number [1:CAS:528:DC%2BD2cXmtlOmtrg%3D]
Publisher Full Text 
Phys. Rev. Lett.. 2006, 97:195501.
COI number [1:STN:280:DC%2BD28jitlCnsw%3D%3D]
Publisher Full Text 
Li WH, Wei BC, Zhang TH, Xing DM, Zhang LC, Wang YR:
Intermetallics. 2007, 15:706. Publisher Full Text