Journal of Nanomaterials
Volume 2016 (2016), Article ID 9265948, 8 pages
http://dx.doi.org/10.1155/2016/9265948

Research Article

Molecular Dynamics Simulation of Nanoindentation of Cu/Au Thin Films at Different Temperatures

Qibin Li,1,2 Cheng Huang,² Yunpei Liang,¹ Tao Fu,² and Tiefeng Peng¹

¹State Key Laboratory of Coal Mine Disaster Dynamics and Control, Chongqing University, Chongqing 400044, China
²College of Aerospace Engineering, Chongqing Key Laboratory of Heterogeneous Material Mechanics, Chongqing University, Chongqing 400044, China

Received 4 September 2016; Accepted 25 October 2016

Academic Editor: Fengyuan Lu

Copyright © 2016 Qibin Li 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.

Two methods, deposition method and ideal modeling based on lattice constant, are used to prepare three modulation periods’ (1.8 nm Cu/3.6 nm Au, 2.7 nm Cu/2.7 nm Au, and 3.6 nm Cu/1.8 nm Au) thin films for nanoindentation at different temperatures. The results show that the temperature will weaken the hardness of thin films. The deposition method and the formation of coherent interface will result in a lot of defects in thin films. These defects can reduce the residual stress in the thin films which is caused by the external force. The proposed system will provide potential benefits in designing the microstructures for thin films.

1. Introduction

Coatings play important roles in the modern manufacturing and processing industry [1]. The tools and engines obtain certain characteristics (e.g., super-hardness, excellent antiwear capability, and high heat-resistance) by coating specific materials on the surface. Recently, nanomultilayer coatings (NMLCs) [2] have drawn significant attentions for their multifunction applications in comparison to monolayer coatings, as the NMLCs are promising candidatures to use in the field of aerospace, manufacture, automobile, and so on.

Several novel experimental methods (e.g., nanoindentation, nanoscratch) were proposed to examine the properties of NMLCs. Li et al. [3] found two different deformation mechanisms, that is, elastic buckling-assisted grain boundary sliding and dislocation plasticity-dominated shearing, of Au/Cu multilayer via nanoindentation tests. However, the micromechanisms of NMLCs are difficult to reveal by these experiments. Therefore, analytical models [4] and numerical methods [5–7] were used to cover the shortage of experimental measurement.

Molecular dynamics (MD) simulations [8] are developed based on the classical Newton and statistical mechanics, and have become a powerful tool to investigate a serial of complex systems [9–12] or conditions [13–17], which can be widely applicable in the area of nanomaterials. Especially, MD can be used to examine the materials in extreme conditions. Landman et al. [18] used a nickel tip to perform the indentation of a gold surface by MD. The gold surface experienced onset of plastic deformation. Zientarski and Chocyk [19] examined the structure and stress in Cu/Au and Fe/Au thin films. The compressive stress was observed in Cu/Au thin films while tensile stress was found in Fe/Au thin films.

Many interesting mechanisms of NMLCs had been revealed based on MD [20, 21]. However, the studies of coatings are still intensive because of the importance of their application. Meanwhile, there are certain thermal conditions (e.g., high temperature) that exist during the real working process of thin films. Thus, in the present paper, MD is employed to perform the nanoindentation of Cu/Au thin films at different temperatures, to investigate the detail molecular mechanisms regarding the preparation methods.

2. Model and Computational Method

2.1. Preparation of Cu/Au Thin Films

The preparation of Cu/Au thin films is based on the deposit system (PVD method) and is shown in Figure 1. The simulation box is 7.2 × 7.2 × 30 nm³ in dimensions. The Au atoms are randomly generated in the region of 25–27.5 nm in -axis. Then, Au atoms deposit on Cu-substrate with the incident kinetic energy of 0.15 eV.

Figure 1: Simulation model system of deposition.

The substrate in Figure 1 consists of 16,000 Cu atoms which is 20 × 20 × 10 unit cells. This system is adapted from solid-liquid interface model [20, 22, 23]. The atoms at bottom two layers of substrate are fixed. Meanwhile, the other Cu atoms of substrate equilibrate in NVT ensemble at 300 K for 1 ns. Next, the atoms at top six layers of substrate are simulated in NVE ensemble for 1 ns, while the other atoms’ interplaying remains unchanged. After that, Au atoms incident to the substrate in 1 ns with the deposition rate of 7.5 × 10²⁴ atom/s·cm². The Au particles are simulated in NVE ensemble. While the Cu atoms in three circumstances, that is, NVE ensemble, NVT ensemble and fixed condition, interact with each other, as the Figure 1 presented.

Three modulation periods’ thin films are obtained based on the model of deposition, which are 1.8 nm Cu/3.6 nm Au, 2.7 nm Cu/2.7 nm Au, and 3.6 nm Cu/1.8 nm Au. Also, the same modulation periods are used to generate ideal thin films without defect, based on the lattice constant of Cu and Au. The thin films studied in this work are denoted as Deposit 1.8, Deposit 2.7, Deposit 3.6, Ideal 1.8, Ideal 2.7, and Ideal 3.6. For example, Deposit 1.8 represents the 1.8 nm Cu/3.6 nm Au thin film prepared by the deposition model. The microstructures and thermophysical properties of these thin films had been discussed in our previous paper [20, 24].

2.2. Nanoindentation of Thin Films

The six obtained thin films are annealed (heated to 1,000 K and equilibrated for a short period, then cooled to 300 K) firstly. Then some atoms of the samples were removed to flatten the sample surface, for that the surface of the deposited sample is not flat. Next, the thin films interacted with an ideal spherical indenter with radius of 2 nm at a constant velocity of 50 m/s in -axis to obtain the load-depth (loading and unloading) curves. The parameter of radius and the impact velocity had been studied in another paper [6]. The load-depth is 1.6~2.2 nm for each thin film. Besides, the nanoindentation of these thin films is tested at temperatures of 0, 300, 700, and 1000 K, respectively. It should be noted that residual stress remains in thin film during the preparation. The goal of this paper is to reveal the atomic mechanisms in thin films regarding the preparation methods. We assumed there is no residual stress in the film before indentation.

2.3. Forcefield and Computational Parameters

MD simulations are performed by using LAMMPS (large-scale atomic/molecular massively parallel simulator) [25]. The EAM forcefield [26] is chosen as the interaction potential. Periodic boundary conditions are applied in and coordinates. Two ideal walls are created perpendicular to -axis at bottom and up surfaces of the simulation box as the boundary. The time-step is set as 1 femtosecond (fs) in the simulations. This time-step is tested for energy conservation. The Langevin method is employed to conserve the temperature [24].

3. Results and Discussion

The different microstructures will result in different properties of materials. The discussion of this paper will focus on the mechanical properties and deformation mechanisms of these thin films at different temperatures, since the microstructures of these thin films had been analyzed previously [20, 24].

3.1. Indentation Force-Depth Curves

The loading and unloading process of indentation can be illustrated by the force-depth (P-h) curves, which is the relationship between the force calculated from the indenter and depth of the indenter. The P-h curves of the investigated thin films are presented in Figure 2. Typically, in Figure 2(a), OA curve represents the loading process whose indenter moves into the thin films. And AB curve represents the unloading process whose indenter moves out of the thin films. The changes in these curves are as a result of the deformation in thin films.

(a) Deposit 1.8

(b) Ideal 1.8

(c) Deposit 2.7

(d) Ideal 2.7

(e) Deposit 3.6

(f) Ideal 3.6

Figure 2: Indentation P-h curves for the deposit and ideal system.

Usually, P-h curves at 0 K are used to demonstrate the mechanism of indentation tested by MD. As shown in Figure 2, the force at 0 K is larger than that of other higher temperatures, that is, 300 K, 700 K, and 1000 K, even if the indenter reaches deeper depth. There is little difference of P-h curves among 300 K, 700 K, and 1000 K. However, there is a trend among the forces at different higher temperatures, which is . In other words, all of the P-h curves agree with the previous reports [27, 28] that the hardness decreases with the temperature increasing. The atoms vibrate harder as the temperature increase. The atoms are easier moving by the external force at higher temperature. Thus, the force on the indenter is smaller at higher temperature.

3.2. Deformation Mechanism

The deformation of thin films during the nanoindentation can be visualized and illustrated by the lattice evolution. The face centered cubic (FCC) structure is the lattice structure for the pure Cu/Au metal. The indentation will move the atoms and form the defects in thin films. By visualization, this makes the lattice structure transform from FCC to other structures (such as amorphous, HCP, and BCC). Therefore, the other structures (except the FCC) are analyzed to study the deformation. And it should be noted that the lattice structure is complex and hard to distinguish by naked eye when the atoms vibrate at high temperature. Thus, the lattice evolution (without FCC) of thin films during nanoindentation at 0 K is presented as Figure 3. The lattice structure is visualized and analyzed based on the open software OVITO [29].

(a) Deposit 1.8

(b) Ideal 1.8

(c) Deposit 2.7

(d) Ideal 2.7

(e) Deposit 3.6

(f) Ideal 3.6

Figure 3: The lattice evolution of thin films during nanoindentation at 0 K.

As discussed in our previous paper [20, 24], there are lots of defects in the thin films which prepared by deposition method. As shown in Figures 3(a), 3(c), and 3(e), there are lots of visualized lattice structures in the Deposit 1.8, Deposit 2.7, and Deposit 3.6 systems. And there are defects located at the coherent interface of Ideal 1.8, Ideal 2.7, and Ideal 3.6 systems. There are stacking faults in Deposit 1.8 system.

During the loading process, the indenter contacts with the thin film and causes the dislocation and slipping in the thin film. It should be noted that there are lots of defects in the thin film, and the dislocation and slipping is found under the stacking faults. The indenter did not pass through the Cu/Au coherent interface; however, the defects caused by the indenter transferred to the coherent interface. In the unloading process of Deposit 1.8 system, the defects reduced as the indenter removed from the thin film. However, there is residual stress in the thin film after the indenter reached its initial position. The residual stress will hold the atoms and hinder the thin film recover to its original configuration. The deformation process is clearer in the ideal systems.

In the Ideal 1.8 system, the indenter caused large defects in the thin film. The dislocation and slipping is more pronounced than that of the Deposit 1.8 system. And stronger residual stress remains in the thin films. It can be concluded that the defects caused by the deposition method could reduce the residual stress in a certain degree. The similar phenomena are also found in Deposit 2.7, Deposit 3.6, Ideal 2.7, and Ideal 3.6 systems. Interestingly, in the Deposit 3.6 and Ideal 3.6 systems, the indenter passed through the coherent interface and caused the defects in the Cu region. While these defects were eliminated after the indenter was removed from the thin films, as discussed earlier, there are defects in the coherent interface, which could also reduce or even eliminate the residual stress of the region under the coherent interface.

4. Conclusions

The Cu/Au thin films, which are obtained via deposition method and modeled based on lattice constant, with three modulation periods (1.8 nm Cu/3.6 nm Au, 2.7 nm Cu/2.7 nm Au, and 3.6 nm Cu/1.8 nm Au) are employed to perform the nanoindentation test at different temperatures. The results could draw the following conclusions.

The hardness of the thin films is weaker as the temperature increase. The lattice evolution shows that thin films that obtained based on deposition method own more defects. And there are defects in the coherent interface of thin films. These defects could reduce the residual stress in the thin films after nanoindentation.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

This work is financially supported by the State Key Research Development Program of China (Grant no. 2016YFC0801404), National Natural Science Foundation of China (nos. 11332013 and 11272364), and Fundamental Research Funds for the Central Universities (nos. 104167693 and 106112015CDJZR245520).

References

L. B. Freund and S. Suresh, Thin Film Materials: Stress, Defect Formation and Surface Evolution, Cambridge University Press, 2003.

J. S. Koehler, “Attempt to design a strong solid,” Physical Review B, vol. 2, no. 2, pp. 547–551, 1970.

Y. P. Li, X. F. Zhu, G. P. Zhang, J. Tan, W. Wang, and B. Wu, “Investigation of deformation instability of Au/Cu multilayers by indentation,” Philosophical Magazine, vol. 90, no. 22, pp. 3049–3067, 2010.

Y. Song, X. Zhuan, T. J. Wang, and X. Chen, “Evolution of thermal stress in a coating/substrate system during the cooling process of fabrication,” Mechanics of Materials, vol. 74, pp. 26–40, 2014.

D. Yin, X. Peng, Y. Qin, J. Feng, and Z. Wang, “Quantifying adhesion energy of mechanical coatings at atomistic scale,” Applied Surface Science, vol. 258, no. 4, pp. 1451–1455, 2011.

T. Fu, X. Peng, C. Huang, D. Yin, Q. Li, and Z. Wang, “Molecular dynamics simulation of VN thin films under indentation,” Applied Surface Science, vol. 357, pp. 643–650, 2015.

T. H. Wang, T.-H. Fang, and Y.-C. Lin, “Finite-element analysis of the mechanical behavior of Au/Cu and Cu/Au multilayers on silicon substrate under nanoindentation,” Applied Physics A, vol. 90, no. 3, pp. 457–463, 2008.

D. Frenkel and B. Smit, Understanding Molecular Simulation: from Algorithms to Applications, Academic Press, New York, NY, USA, 2nd edition, 2002.

Y. Qiao, L. Liu, and X. Chen, “Pressurized liquid in nanopores: a modified Laplace-Young equation,” Nano Letters, vol. 9, no. 3, pp. 984–988, 2009.

Y. Liu and X. Chen, “High permeability and salt rejection reverse osmosis by a zeolite nano-membrane,” Physical Chemistry Chemical Physics, vol. 15, no. 18, pp. 6817–6824, 2013.

S. Lee, “Molecular Dynamics Study of the Separation Behavior at the Interface between PVDF binder and copper current collector,” Journal of Nanomaterials, vol. 2016, Article ID 4253986, 12 pages, 2016.

F. L. Zhu, H. Y. Liao, K. Tang, Y. K. Chen, and S. Liu, “Molecular dynamics study on the effect of temperature on the tensile properties of single-walled carbon nanotubes with a Ni-coating,” Journal of Nanomaterials, vol. 2015, Article ID 767182, 7 pages, 2015.

H. Du, J. C. Rasaiah, and J. D. Miller, “Structural and dynamic properties of concentrated alkali halide solutions: a molecular dynamics simulation study,” The Journal of Physical Chemistry B, vol. 111, no. 1, pp. 209–217, 2007.

H. Hu, C. R. Weinberger, and Y. Sun, “Model of meniscus shape and disjoining pressure of thin liquid films on nanostructured surfaces with electrostatic interactions,” Journal of Physical Chemistry C, vol. 119, no. 21, pp. 11777–11785, 2015.

H. Hu, C. R. Weinberger, and Y. Sun, “Effect of nanostructures on the meniscus shape and disjoining pressure of ultrathin liquid film,” Nano Letters, vol. 14, no. 12, pp. 7131–7137, 2014.

H. F. Zhan, G. Zhang, J. M. Bell, and Y. T. Gu, “Thermal conductivity of configurable two-dimensional carbon nanotube architecture and strain modulation,” Applied Physics Letters, vol. 105, no. 15, Article ID 153105, 2014.

H. Zhan, G. Zhang, Y. Zhang, V. B. C. Tan, J. M. Bell, and Y. Gu, “Thermal conductivity of a new carbon nanotube analog: the diamond nanothread,” Carbon, vol. 98, pp. 232–237, 2016.

U. Landman, W. D. Luedtke, N. A. Burnham, and R. J. Colton, “Atomistic mechanisms and dynamics of adhesion, nanoindentation, and fracture,” Science, vol. 248, no. 4954, pp. 454–461, 1990.

T. Zientarski and D. Chocyk, “Structure and stress in Cu/Au and Fe/Au systems: a molecular dynamics study,” Thin Solid Films, vol. 562, pp. 347–352, 2014.

Q. Li, X. Peng, T. Peng, Q. Tang, X. Zhang, and C. Huang, “Molecular dynamics simulation of Cu/Au thin films under temperature gradient,” Applied Surface Science, vol. 357, pp. 1823–1829, 2015.

T. Fu, X. Peng, Y. Zhao et al., “MD simulation of effect of crystal orientations and substrate temperature on growth of Cu/Ni bilayer films,” Applied Physics A: Materials Science and Processing, vol. 122, no. 2, article 67, pp. 1–9, 2016.

Q. Li and C. Liu, “Molecular dynamics simulation of heat transfer with effects of fluid-lattice interactions,” International Journal of Heat and Mass Transfer, vol. 55, no. 25-26, pp. 8088–8092, 2012.

Q. Li, Q. Tang, T. Peng, X. Zhang, C. Liu, and X. Shi, “Molecular characteristics of H₂O in hydrate/ice/liquid water mixture,” International Journal of Modern Physics B, vol. 29, no. 27, Article ID 1550185, 2015.

Q. Li, X. Peng, T. Peng, Q. Tang, C. Liu, and X. Shi, “Molecular dynamics simulations of coating process: influences of thermostat methods,” Journal of Computational and Theoretical Nanoscience, vol. 13, no. 7, pp. 4629–4633, 2016.

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

X. W. Zhou, R. A. Johnson, and H. N. G. Wadley, “Misfit-energy-increasing dislocations in vapor-deposited CoFe/NiFe multilayers,” Physical Review B, vol. 69, no. 14, Article ID 144113, 2004.

H. Wang, M. Hu, M. Xia, F. Ke, and Y. Bai, “Molecular/cluster statistical thermodynamics methods to simulate quasi-static deformations at finite temperature,” International Journal of Solids and Structures, vol. 45, no. 13, pp. 3918–3933, 2008.

S. Amaya-Roncancio, D. F. Arias-Mateus, M. M. Gómez-Hermida, J. C. Riaño-Rojas, and E. Restrepo-Parra, “Molecular dynamics simulations of the temperature effect in the hardness on Cr and CrN films,” Applied Surface Science, vol. 258, no. 10, pp. 4473–4477, 2012.

A. Stukowski, “Visualization and analysis of atomistic simulation data with OVITO–the Open Visualization Tool,” Modelling and Simulation in Materials Science and Engineering, vol. 18, no. 1, Article ID 015012, 2010.