Open Access Nano Express

Atomistic insight into the minimum wear depth of Cu(111) surface

Zengqiang Li1, Yanhua Huang2, Junjie Zhang1*, Yongda Yan1 and Tao Sun1

Author Affiliations

1 Center for Precision Engineering, Harbin Institute of Technology, Harbin, 150001, People’s Republic of China

2 Research Center of Laser Fusion, China Academy of Engineering Physics, Mianyang 621900, People’s Republic of China

For all author emails, please log on.

Nanoscale Research Letters 2013, 8:514  doi:10.1186/1556-276X-8-514


The electronic version of this article is the complete one and can be found online at: http://www.nanoscalereslett.com/content/8/1/514


Received:13 October 2013
Accepted:27 November 2013
Published:5 December 2013

© 2013 Li et al.; licensee Springer.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

In the present work, we investigate the minimum wear depth of single crystalline Cu(111) under single asperity friction by means of molecular dynamics simulations. The atomistic mechanisms governing the incipient plasticity are elucidated by characterizing specific defect structures and are correlated to the observed mechanical and frictional responses of the material. Furthermore, the effect of probe radius on the friction process is studied. Our simulations indicate that the formation of wear impression is closely associated with defect nucleation and the minimum wear depth is equivalent to the critical penetration depth at which plasticity initiates. It is found that the probe radius has a strong influence on the formation of defect structures and the observed mechanical responses.

Keywords:
Single asperity friction; Friction and wear; Single crystalline Cu; Incipient plasticity; Molecular dynamics

Background

With the feature size of miniaturized mechanical components shrinking down to the nanometer regime, friction and wear, as the major causes of mechanical failures and dissipative energy losses, play pronounced and even dominant role in determining the functionality of nanoelectromechanical system (NEMS) devices [1-3]. Therefore, reducing the friction and wear between contacting surfaces of components is of significant importance for the application of NEMS devices. Specifically, a fundamental understanding of the atomic scale origin of the friction-induced wear is essentially required for the rational design of the components that possess good wear resistance.

During the course of friction, wear phenomena are closely accompanied with permanent deformation and even removal of the materials under applied mechanical loads. Thus, identifying and characterizing the initiation of plasticity of the materials under friction are central to the understanding of the atomic scale origin of wear phenomena. In the past few decades, both experimental investigations and atomistic simulations have been conducted to investigate the incipient plasticity of metallic and semiconductor materials under nanoindentation [4-8]. Recently, Paul et al. performed nanoindentation experiments to study the minimum threshold of the incipient plasticity of a gold single crystal. They found that the indentation-induced elastic deformation and plastic deformation can be well identified by features observed in the force-displacement curves, and the first pop-in phenomenon reflects the onset of plasticity [9]. However, a rather limited effort has been taken to study the incipient plasticity which occurs under friction. Compared to the localized uniaxial stress state of nanoindentation, the multi-axial states of localized stress induced by friction action may lead to more complex mechanical responses at the onset of plasticity. On the other hand, it is crucial to correlate microstructure evolution that occurs within the materials with the observed features in force-displacement curves, which is of great challenge for the experimental investigations because of the involvement of nanometer length scale. As a complement to experiments, molecular dynamics (MD) simulation has been demonstrated to be one powerful tool to investigate the atomic scale phenomena of friction and wear [10-20]. Although previous MD simulations have provided valuable insights into the nanoscale friction and wear processes, our knowledge about the incipient plasticity under friction process, particularly the relationship between specific defect structures and observed wear phenomena, is still scarce.

In the present work, we perform MD simulations to investigate the incipient plasticity of single crystalline copper under single asperity friction with a spherical probe. The deformation mechanisms of the material are analyzed in detail, and the specific defect structures are particularly characterized and are correlated to the mechanical and frictional responses. Our simulations demonstrate that the minimum wear depth is determined by the formation of permanent defects such as dislocations and vacancies and is strongly probe radius-dependent. This paper is outlined as follows. In ‘Methods’ Section, we describe the simulation method. The simulation results are presented and discussed in Section ‘Results and discussion’. Finally, we summarize the results in ‘Conclusions’ Section.

Methods

As depicted in Figure 1, the MD model of single asperity friction employed in the present work consists of a substrate and a spherical probe. The substrate of single crystalline copper has a dimension of 30, 10, and 30 nm in X[2], Y [111], and Z[1-10] directions, respectively. Periodic boundary conditions are imposed in the transverse X and Z directions of the substrate. Figure 1 shows that the substrate is composed of two virtual types of atoms, as the green color stands for the fixed atoms and the red one represents the mobile atoms in which motions follow the Newton’s second law of motion. The atomic interactions within the substrate are described by an embedded atom method developed for copper [21]. The frictionless spherical probe is modeled by a strong repulsive potential [22]. To study the influence of probe radius on the friction, four probe radiuses of 6, 8, 10, and 12 nm are considered.

thumbnailFigure 1. MD model of single asperity friction of single crystalline copper. The atoms in the substrate are colored according to their virtual types, as red for mobile atoms and green for fixed atoms.

The atoms in the as-created substrate first undergo global energy minimization at 0 K, and then the substrate is relaxed to its equilibrium configuration at 30 K and 0 bar through dynamic NPT relaxation for 50 ps. After relaxation, the substrate is subjected to friction by placing the probe above the free surface of the substrate with a distance of 0.2 nm. The friction process is composed of two stages of first penetration and following scratching, as illustrated in Figure 1. In the penetration stage, the probe moves along negative Y direction with constant velocity of 20 m/s to penetrate into the substrate until a pre-determined penetration depth is reached. In the following scratching stage, the probe scratches at 12.2 nm along negative X direction with constant velocity of 20 m/s. Both the penetration and scratching velocities of 20 m/s are a few orders of magnitude higher than the typical velocities utilized in nanoscratching experiments due to the intrinsic requirement of integration timesteps to be of the order of 1 fs. All the MD simulations are completed using the IMD code with a time step of 1 fs [23]. The detailed description about the friction procedure can also be found elsewhere [24]. To identify the defects generated within the substrate, a modified bond angle distribution (BAD) method is utilized [25]. In the present work, the perfect face-centered cubic (FCC) atoms are not shown for better viewing of the defect structures, and the coloring scheme for various defects is as follows: red stands for surface atoms, blue indicates hexagonal close-packed (HCP) atoms, and the remaining atoms are categorized into defects including dislocation cores and vacancies. The software AtomEye is employed to visualize MD data and generate MD snapshots [26].

Results and discussion

Determination of minimum wear depth

In the friction process, there are three force components acting on the probe, as scratching force along X direction, penetration force along Y direction, and lateral force along Z direction, respectively. In the penetration stage, both scratching force and lateral force mainly fluctuate around constant value of 0 because the probe only applies uniaxial localized stress along Y direction. Figure 2 plots the penetration force-penetration depth curve during the penetration stage with a probe radius of 8 nm, indicating that the deformation behavior of the substrate is divided into two regimes. In the regime I, the substrate undergoes elastic deformation, accompanied with rapid increase of the penetration force. After the penetration depth reaches a critical value of 0.72 nm, the penetration force drops precipitously, indicating the occurrence of elastic deformation-plastic deformation transition. The observed phenomenon of force drop, which corresponds to the pop-in event widely observed in the load-controlled nanoindentation experiments, is caused by dislocation avalanche beneath the penetrated surface [5,7,24]. We note that the tribochemistry, e.g., the presence of cupric oxide, may significantly alter the deformation behavior of the topmost surface. In the regime II, the substrate undergoes plastic deformation dominated by dislocation activities. The action of penetration stops at a penetration depth D2 of 0.82 nm. Another penetration depth D1 of 0.65 nm in the elastic deformation regime, at which the penetration force is equal to that at D2, is also marked in Figure 2. The two insets in Figure 2 present instantaneous defect structures obtained at the two penetration depths D1 and D2, respectively. While the substrate is purely elastically deformed at D1, there is a considerable amount of defects formed beneath the penetrated surface at D2.

thumbnailFigure 2. Penetration force-penetration depth curve during the penetration with a probe radius of 8 nm. The two penetration depths D1 of 0.65 nm and D2 of 0.82 nm have the same penetration force. The two insets show instantaneous defect structures at D1 and D2, in which atoms are colored according to their BAD values and FCC atoms are not shown.

While Figure 2 shows that the defect structures at the two penetration depths are significantly different, two scratching simulations under the two scratching depths D1 and D2 are conducted with the same probe radius of 8 nm. Under the scratching depth D1, both the penetration force and scratching force remain constant values throughout the scratching stage. However, the scratching force is far smaller than the penetration force because of the absence of permanent deformation in the vicinity of the probe. We also note that the non-adhesion between the substrate and the probe in the current simulated system also contributes to the ultra-small scratching force. In contrast, under the scratching depth D2, the friction coefficient first increases rapidly and then fluctuates around constant value with large amplitudes when the scratching becomes stable. To examine the evolutions of defect structures and surface morphologies, retractions of the probe along Y direction to its initial height are conducted right after the completion of the two scratching stages. Figure 3 presents instantaneous defect structures and surface morphologies of the substrate after the completion of scratching and retraction for the two scratching depths. We note that the following observations are made based on not only the captured MD snapshots, but also the entire dynamic process provided by MD simulations: under the scratching depth D1, the substrate undergoes pure elastic deformation, and there is no defect formed beneath the surface after the completion of the scratching, as shown in Figure 3a. Accordingly, there is only one penetration impression formed on the surface shown in Figure 3e. Furthermore, Figure 3b,f demonstrates that the penetrated surface is fully recovered after the retraction, indicating that there is no permanent deformation that occurs within the substrate. Under the scratching depth D2, however, it is seen from Figure 3c that the defect zone beneath the surface extends significantly along the scratching direction. Figure 3g shows that there is one scratching-induced impression of the groove formed on the surface, and wear debris which accumulate on both sides of the groove is also observed. Although the penetrated surface undergoes tiny plastic recovery accompanied by the shrinking of the defect structures beneath the probe after the retraction, Figure 3d,h shows that both the defect structures, particularly those behind the probe, and the surface morphology are mainly unchanged. Furthermore, the height of wear debris increases slightly due to the annihilation of the dislocations at the surface [24].

thumbnailFigure 3. Defect structures and surface morphologies after scratching and retraction under D1 and D2 (a,b,c,d). Defect structures after scratching and retraction under the scratching depths D1 and D2, respectively. Atoms are colored according to their BAD values, and FCC atoms are not shown. (e,f,g,h) Surface morphologies after scratching and retraction under the scratching depths D1 and D2, respectively. Atoms are colored according to their heights in Y direction.

The above analysis indicates that the minimum wear depth is closely associated with the initiation of plasticity. To reveal the specific defect structures formed at the early stage of plastic deformation, a dynamic inspection of the defect evolution in the regime II of Figure 2 is performed. Figure 4a shows that at the critical penetration depth of 0.72 nm a dislocation loop formed on one {111} slip plane inclined to the (111) free surface, which leads to the sharp drop of the penetration force observed in Figure 2. The Burgers vector of each segment of the dislocation loop in Figure 4a is analyzed through DXA algorithm [27], and the corresponding dislocation network is presented in Figure 4b. Dislocation cores are represented by thin tubes, in which Shockley partial dislocation with 1/6 <112 > Burgers vector and perfect dislocation with 1/2 <110 > Burgers vector are colored gray and red, respectively. It is seen from Figure 4b that the dislocation loop consists of four partial dislocations and one perfect dislocation. In addition, there is one vacancy formed beneath the probe. Upon further penetration, the other three {111} slip planes are activated sequentially, and Figure 4c shows that the defect zone beneath the probe expands greatly. The glide of dislocations on adjacent slip planes leads to the formation of stair-rod dislocations with 1/6 <110 > Burgers vector highlighted by the arrows in Figure 4d. Figure 4e,f presents dislocation network after the completion of scratching and penetration, respectively. It is seen from Figure 4e that there is less dislocations but more vacancies in the wake of the probe than that in the vicinity of the probe due to the plastic recovery. In addition to the stair-rod dislocations, there are glissile prismatic dislocation loops formed by dislocation reaction and cross-slip events. In particular, the prismatic dislocation half-loops in front of the probe glide parallels to the free surface to transport the materials displaced by the probe without the formation of surface steps [24]. Although small part of the dislocations beneath the probe annihilates at the free surface during the retraction, Figure 4f shows that the defect structures are stable.

thumbnailFigure 4. Close inspections of defect structures in friction with a probe radius of 8 nm. The scratching depth is 0.82 nm. (a,c) Bottom views of defect structures at penetration depths of 0.72 and 0.82 nm, respectively. Atoms are colored according to their BAD values and FCC atoms are not shown. (b,d) Dislocation networks shown in (a) and (c), respectively. (e,f) Dislocation networks after the completion of scratching and retraction, respectively.

Effect of probe radius on minimum wear depth

To investigate the influence of probe radius on the minimum wear depth, friction simulations with another three probe radiuses of 6, 10, and 12 nm are conducted, in addition to the probe radius of 8 nm. For each probe radius, the penetration stage stops at a penetration depth that is 0.1 nm deeper than the critical penetration depth at which the phenomenon of force drop occurs. Figure 5a,b plots the contact pressure-penetration depth curves and the friction coefficient-scratching length curves during the penetration and scratching stages with the four probe radiuses, respectively. The contact pressure is defined as the ratio of the penetration force to the contact area. A detailed description about the calculation of the contact area during spherical penetration can be found elsewhere [28]. It should be noted that the maximum contact pressure shown in Figure 5a corresponds to the theoretical strength of the pristine and dislocation-free single-crystalline Cu. Furthermore, to quantitatively access the influence of probe radius on the frictional property of the substrate, the average friction coefficient is obtained by averaging more than 1,000 instantaneous points of friction coefficient in the range between 3 and 12.2 nm. Table 1 summarizes the mechanical responses of the substrate extracted during friction with the four probe radiuses. Figure 5a shows that the slope of the contact pressure-penetration depth curve in the elastic deformation regime decreases with increasing probe radius, indicating that the elastic deformation of the substrate is more compliant with the larger probe. However, the contact pressure reflecting the critical stress for initial dislocation nucleation from penetrated surface is approximately independent on the probe radius. It is seen from Table 1 that with the increase of the probe radius, both the critical force and the critical penetration depth associated with the initiation of plasticity increases, but the average friction coefficient decreases.

thumbnailFigure 5. Influence of probe radius on mechanical and frictional properties of the substrate under friction. (a) Contact pressure-penetration depth curves. (b) Friction coefficient-scratching length curves.

Table 1. Mechanical responses of the substrate under friction with different probe radiuses

Figure 6a,b,c,d presents the surface morphologies of the substrate after the completion of scratching with probe radiuses of 6, 8, 10, and 12 nm, respectively. A larger probe results in a larger volume and also wider extent of the wear debris, indicating that more atoms within the substrate are involved in the scratching action. To quantitatively characterize the scratching-induced motion of atoms, the shear strain of each atom is calculated by comparing the current atomic configuration of the substrate with the reference configuration obtained after relaxation. Figure 6e,f presents the cross-sectional views of the substrate after scratching with the four probe radiuses, respectively, in which atoms are colored according to their shear strains ranging from 0 to 1. It is seen from Figure 6 that the distributions of wear debris and shear strain are closely correlated for each probe radius. When probe radius is small, Figure 6e shows that the distribution of shear strain is compact and shallow. Furthermore, the atoms in the wear debris have significantly larger mobility than that within the material. In contrast, a lager probe leads to larger and more compliant distribution of shear strain.

thumbnailFigure 6. Influence of probe radius on the friction of the substrate. (a,b,c,d) Surface morphologies of the substrate after scratching with probe radiuses of 6, 8, 10, and 12 nm, respectively. Atoms are colored according to their height in Y direction. (e,f,g,h) Cross-sectional views of the substrate after scratching with probe radiuses of 6, 8, 10, and 12 nm, respectively. Atoms are colored according to shear strain ranging from 0 to 1.

Figure 7 presents numbers of HCP and defect atoms generated within the substrate after penetration and scratching with the four probe radiuses. For each probe, there are more HCP and defect atoms generated in the scratching stage than that in penetration stage, because of the more complex plastic deformation associated with the multi-axial localized stress states. When the probe radius is not larger than 10 nm, there are more defect atoms than HCP atoms in both penetration and scratching stages for each probe radius. However, the friction with the probe radius of 12 nm results in more HCP atoms than defect atoms generated within the material. The formation of HCP atoms is associated with the activity of partial dislocations, while defect atoms are composed of not only dislocation cores but also vacancies. Therefore, Figure 7 indicates that the dislocation activity plays more pronounced role in governing incipient plasticity for larger probe. In addition, the incipient plasticity shows strong dependence on probe radius: the larger the probe, the larger both the HCP and defect atoms.

thumbnailFigure 7. Influence of probe radius on numbers of HCP and defect atoms generated within the substrate under friction.

Conclusions

In summary, we perform MD simulations to investigate the atomic scale origin of the minimum wear depth of single crystalline Cu(111) during single asperity friction. Simulation results show that scratching impression can only be made under a scratching depth at which there are permanent defects formed. It is indicated that the minimum wear depth is equivalent to the critical penetration depth associated with the first force-drop observed in the force-depth curve. The specific permanent defects governing the wear phenomena are composed of stair-rod dislocations and prismatic dislocation loops as well as vacancies. While the contact pressure for the nucleation of initial dislocation is independent on probe radius, the minimum wear depth increases with probe radius. Further analysis of the shear strain distribution implies that a larger probe results in more compliant deformation of the material, which leads to larger volume of wear debris and wider extent of defect structures.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

JZ, YY, and TS conceived the project. LZ, YH, and JZ performed molecular dynamics simulations and analyzed data. LZ and JZ wrote the paper. All authors read and approved the final manuscript.

Acknowledgements

The authors greatly acknowledge financial supports from the NSFC (51005059 and 51222504), China Postdoctoral Science Foundation (20100471047 and 2012 M511463), and Heilongjiang Postdoctoral Foundation of China (LBH-Z11143). JZ also greatly acknowledges Dr. Alexander Hartmaier and Dr. Christoph Begau in ICAMS, Ruhr-University Bochum for helpful discussion and providing advanced analyzing technique of lattice defects.

References

  1. Spengen W, Modlinski R, Puers R, Jourdain A: Failure mechanisms in MEMS/NEMS devices. In Springer Handbook of Nanotechnology. Berlin: Springer; 2007:1663-1684. OpenURL

  2. Bhushan B: Nanotribology and nanomechanics of MEMS/NEMS and BioMEMS/BioNEMS materials and devices.

    Microelectron Eng 2007, 84:387-412. Publisher Full Text OpenURL

  3. Kim HJ, Yoo SS, Kim DE: Nano-scale wear: a review.

    Int J Precis Eng Man 2013, 13:1709-1718. OpenURL

  4. Tadmor EB, Miller R, Phillips R, Ortiz M: Nanoindentation and incipient plasticity.

    J Mater Res 1999, 14:2233-2250. Publisher Full Text OpenURL

  5. Li J, Vliet KJV, Zhu T, Yip S, Suresh S: Atomistic mechanisms governing elastic limit and incipient plasticity in crystals.

    Nature 2002, 418:307-310. PubMed Abstract | Publisher Full Text OpenURL

  6. Lund AC, Hodge AM, Schuh CA: Incipient plasticity during nanoindentation at elevated temperatures.

    Appl Phys Lett 2004, 85:1362. Publisher Full Text OpenURL

  7. Lee YM, Park JY, Kim SY, Jun S, Im SY: Atomistic simulations of incipient plasticity under Al(111) nanoindentation.

    Mech Mater 2005, 37:1035-1048. Publisher Full Text OpenURL

  8. Catoor D, Gao YF, Geng J, Prasad MJNV, Herbert EG, Kumar KS, Pharr GM, George EP: Incipient plasticity and deformation mechanisms in single-crystal Mg during spherical nanoindentation.

    Acta Mater 2013, 61:2953-2965. Publisher Full Text OpenURL

  9. Paul W, Oliver D, Miyahara Y, Grütter PH: Minimum threshold for incipient plasticity in the atomic-scale nanoindentation of Au(111).

    Phys Rev Lett 2013, 110:135506. PubMed Abstract | Publisher Full Text OpenURL

  10. Zhang LC, Tanaka H: Towards a deeper understanding of wear and friction on the atomic scale-a molecular dynamics analysis.

    Wear 1997, 211:44-53. Publisher Full Text OpenURL

  11. Fang TH, Weng CI: Three-dimensional molecular dynamics analysis of processing using a pin tool on the atomic scale.

    Nanotechnology 2000, 11:148-153. Publisher Full Text OpenURL

  12. Zhu PZ, Hu YZ, Ma TB, Wang H: Molecular dynamics study on friction due to ploughing and adhesion in nanometric scratching process.

    Tribol Lett 2011, 41:41-46. Publisher Full Text OpenURL

  13. Zhu PZ, Hu YZ, Wang H, Ma TB: Study of effect of indenter shape in nanometric scratching process using molecular dynamics.

    Mater Sci Eng A 2011, 528:4522-4527. Publisher Full Text OpenURL

  14. Khan HM, Kim SG: On the wear mechanism of thin nickel film during AFM-based scratching process using molecular dynamics.

    J Mech Sci Technol 2011, 25:2111-2120. Publisher Full Text OpenURL

  15. Liu XM, Liu ZL, Wei YG: Nanoscale friction behavior of the Ni-film/substrate system under scratching using MD simulation.

    Tribol Lett 2012, 46:167-178. Publisher Full Text OpenURL

  16. Mishra M, Egberts P, Bennewitz R, Szlufarska I: Friction model for single-asperity elastic–plastic contacts.

    Phys Rev B 2012, 86:045452. OpenURL

  17. Wu CD, Fang TH, Lin JF: Atomic-scale simulations of materials behaviors and tribology properties for FCC and BCC metal films.

    Mater Lett 2012, 80:59-62. OpenURL

  18. Kim CI, Yang SH, Kim YS: Deformation characteristics of various grain boundary angles on AFM-based nanolithography using molecular dynamics.

    J Mech Sci Technol 2012, 26:1841-1847. Publisher Full Text OpenURL

  19. Mishra M, Szlufarska I: Dislocation controlled wear in single crystal silicon carbide.

    J Mater Sci 2013, 48:1593-1603. Publisher Full Text OpenURL

  20. Zhang P, Zhao HW, Shi CL, Zhang L, Huang H, Ren LQ: Influence of double-tip scratch and single-tip scratch on nano-scratching process via molecular dynamics simulation.

    Appl Surf Sci 2013, 280:751-756. OpenURL

  21. Mishin Y, Mehl MJ, Papaconstantopoulos DA, Voter AF, Kress JD: Structural stability and lattice defects in copper: ab initio, tight-binding, and embedded-atom calculations.

    Phys Rev B 2001, 63:224106. OpenURL

  22. Vliet KJV, Li J, Zhu T, Yip S, Suresh S: Quantifying the early stages of plasticity through nanoscale experiments and simulations.

    Phys Rev B 2003, 67:104105. OpenURL

  23. Stadler J, Mikulla R, Trebin HR: IMD: a software package for molecular dynamics studies on parallel computers.

    Int J Mod Phys C 1997, 8:1131-1140. Publisher Full Text OpenURL

  24. Zhang JJ, Sun T, Hartmaier A, Yan YD: Atomistic simulation of the influence of nanomachining-induced deformation on subsequent nanoindentation.

    Comput Mater Sci 2012, 59:14-21. OpenURL

  25. Begau C, Hartmaier A, George EP, Pharr GM: Atomistic processes of dislocation generation and plastic deformation during nanoindentation.

    Acta Mater 2011, 59:934-942. Publisher Full Text OpenURL

  26. Li J: AtomEye: an efficient atomistic configuration viewer.

    Modelling Simul Mater Sci Eng 2003, 11:173-177. Publisher Full Text OpenURL

  27. Stukowski A, Albe K: Extracting dislocations and non-dislocation crystal defects from atomistic simulation data.

    Modelling Simul Mater Sci Eng 2010, 18:085001. Publisher Full Text OpenURL

  28. Ziegenhain G, Urbassek HM, Hartmaier A: Influence of crystal anisotropy on elastic deformation and onset of plasticity in nanoindentation: a simulational study.

    J Appl Phys 2010, 107:061807. Publisher Full Text OpenURL