# Molecular Dynamics Study on the Effect of Temperature on the Tensile Properties of Single-Walled Carbon Nanotubes with a Ni-Coating

Authors: Fulong Zhu; Hengyou Liao; Kai Tang; Youkai Chen; Sheng Liu
Journal: Journal of Nanomaterials
 (2015)

---

## Abstract

The effect of temperature on the tensile behavior of the armchair (6, 6) single-walled carbon nanotubes with a Ni-coating (SWCNT-Ni) was investigated using molecular dynamics (MD) methods. The mechanical properties of SWCNT-Ni and SWCNT were calculated and analyzed at different temperatures in the range from 220 K to 1200 K. From the MD results, temperature was determined to be the crucial factor affecting the mechanical properties of SWCNT-Ni and SWCNT. After coating nickel atoms onto the surface of a SWCNT, the Young’s modulus, tensile strength, and tensile failure strain of SWCNT were greatly reduced with temperature rising, indicating that the nickel atoms on the surface of SWCNT degrade its mechanical properties. However, at high temperature, the Young’s modulus of both the SWCNT and the SWCNT-Ni exhibited significantly greater temperature sensitivity than at low temperatures, as the mechanical properties of SWCNT-Ni were primarily dominated by temperature and C-Ni interactions. During these stretching processes at different temperatures, the nickel atoms on the surface of SWCNT-Ni could obtain the amount of energy sufficient to break the C-C bonds as the temperature increases.


---

## Body

## 1. Introduction

Since discovery in 1991 [1] and due to the excellent properties such as high thermal and electrical conductivity, ultra mechanical strength, low density, and high aspect ratio, carbon nanotubes (CNTs) have gained increasing attention from both academia and industry. The thermal conductivity of multiwalled carbon nanotubes (MWCNT) can reach up to 3000 W/mK [2], and the thermal conductivity of single-walled nanotubes (SWCNT) can be even higher [3]. Some types of CNTs are considered to be ideal candidates as reinforcing materials to enhance the properties of composites and electrical and thermal devices [4, 5].Because of the remarkable mechanical properties of SWCNT, it is believed that the mechanical properties of a metal matrix can be greatly improved by mixing the metal with SWCNT coated by a metal film, such as nickel or copper. For some types of metals, during mixing, the metal could agglomerate into discrete deposits on SWCNT bundles [6, 7], which could decrease the adhesive force between SWCNT and the metal matrix, thus preventing the improvement of the mechanical properties of the metal by mixing the metal with SWCNTs. To increase the adhesive force between SWCNT and some metal atoms, the surface of each SWCNT must be treated by coating some metal films, such as nickel or copper. After coating certain metal films onto the surfaces of SWCNTs, the SWCNTs with a coating of metal film are suitable to mix with the metal matrix, thereby improving the efficiency of load transfer from the metal matrix to the carbon nanotubes. Coating nickel or copper is one conventional method of decorating the surface of SWCNT. Nickel atoms can be coated onto suspended single-walled carbon nanotubes using electron-beam evaporation, as the coating of nickel atoms on the suspended carbon tubes has been demonstrated to be continuous and quasi-continuous [6]. Some metals were also deposited onto a vertically aligned SWCNT (VA-SWCNT) film by evaporative deposition [7]; and the continuous spreading of nickel atoms onto carbon nanotubes was investigated using molecular dynamics methods [8]. Some experiments also indicated that the mechanical and tribological properties of the copper matrix composites could be significantly enhanced by the incorporation of SWCNT reinforcements after coating SWCNT with nickel metal [9]. Meanwhile, the decorated CNTs could also be functionalized to satisfy interconnecting processing or for the application of nanoscale sensors [10].To determine the effect of temperature on the mechanical properties of armchair (6, 6) single-walled carbon nanotubes with a Ni-coating (SWCNT-Ni), a perfect SWCNT decorated with a coating of nickel atoms on its surfaces is considered in this work. Although the influence of nickel atoms on SWCNT was investigated at 0 K and 300 K [11], a wide temperature range will be defined in this paper by considering temperatures of 220 K, 300 K, 500 K, 700 K, 1000 K, and 1200 K. After modeling and analyzing the stretching process of SWCNT and SWCNT-Ni at different temperatures of 220 K, 300 K, 500 K, 700 K, 1000 K, and 1200 K using molecular dynamics methods, the mechanical behaviors of the SWCNT and the SWCNT-Ni can then be determined. In the following molecular dynamics simulation for stretching SWCNT and SWCNT-Ni, a more complicated second-generation Brenner reactive empirical bond-order potential will be introduced for analyzing the behaviors of SWCNT-Ni and SWCNT. In the following section, the stress-strain curves of SWCNT and SWCNT-Ni are determined with the data from these molecular dynamics computational results, their fracture surfaces are also determined, and the influence of temperature on the coating of nickel onto the surface of SWCNT is revealed.


## 2. Computational Methods and Models

To investigate the effect of nickel atoms and temperature on the mechanical behavior of SWCNT-Ni, the surface of an armchair (6, 6) SWCNT was coated with nickel using metal cluster deposition [8]. After coating nickel atoms onto the surface of SWCNT, nickel clusters were demonstrated to spread and form a continuous smooth surface of a VA-SWCNT, in which the equilibrium position for nickel atoms was the center of the hexagons of the carbon atoms of SWCNT [12]. In this work, molecular dynamics simulations of the stretching process at different temperatures for SWCNT-Ni and SWCNT were performed, based on the large-scale atomic/molecular massively parallel simulator (LAMMPS) code [13]. In the molecular dynamics (MD) modeling of SWCNT and SWCNT-Ni, to simplify the cluster deposition process, nickel atoms were simply placed in the vicinity of the center of a carbon hexagon, with every computing step of the molecular dynamics simulation following a period of energy minimization; MD simulations for the stretching process of SWCNT-Ni and SWCNT were conducted at different temperatures of 220 K, 300 K, 500 K, 700 K, 1000 K, and 1200 K, and Figures 1 and 2 show the stretching models. In the stretching diagrams of Figures 1 and 2, the molecular model of the armchair (6, 6) SWCNT consisted of 576 carbon atoms, and the molecular model of the armchair (6, 6) SWCNT-Ni had 576 carbon atoms and 264 nickel atoms. In the two models shown in Figures 1 and 2, the first two rings of atoms located at each model’s right end were fixed, the last two rings of atoms located at the left end were also restricted in their x- and y-directions, and the stretching process occurred along the z-direction. In the MD simulations for SWCNT and SWCNT-Ni, the entire stretching process occurred in an isothermal environment; the C-C interactions were modeled using the second-generation Brenner reactive empirical bond-order potential [14], which has been successfully used to investigate various carbon and hydrocarbon systems; the C-Ni and Ni-Ni interaction parameters were also determined using the Brenner-type potential, and its potential parameters were used from these results [14]. This type of potential could be described by the following formulae:
(1)EijC-Ni=DeS-1exp⁡-β2Srij-Re-B*DeSS-1exp⁡-β2/Srij-Re,B*=1+bNC-1δ,EijNi-Ni=DeS-1exp⁡-β2Srij-Re,-DeSS-1exp⁡-β2/Srij-Re,De=De1+De2exp⁡-CDNM-1,Re=Re1-Re2exp⁡-CRNM-1,NiC=1+∑carbonfrik,NiM=1+∑metalfrik.Figure 1
Stretching model diagram of a SWCNT.Figure 2
Stretching model diagram of a SWCNT-Ni.Here,rij was the distance between the two atoms i and j. De and Re were the binding energy and equilibrium bond length, respectively. S, β, b, and δ determined the shape of the potential functions. For the C-Ni interactions, the many-body effects were defined by using an additional term B*, which was expressed as a function of the carbon coordination number of a metal atom NC. For the Ni-Ni interactions, the binding energy De and the equilibrium bond length Re were expressed as direct functions of the metal coordination number NM. Consider
(2)frij=1r<R1,121+cos⁡πrij-R1R2-R1R1≤r≤R2,0r>R2.The cutoff function equation in the second-generation Brenner reactive empirical bond-order potential introduced a dramatic increase in the interatomic force atr=R1 (similar to a camelback on the force curve) [15], which could rise sharply with a peak near the bond breaking length. Therefore, R1=R2=2.0 Å [16] was used to avoid overestimating the force in the molecular dynamics simulation of this paper.In this work, the structure of SWCNT-Ni was two concentric cylinders, and the cross-sectional of SWCNT-Ni was two concentric circles; the cross-sectional area of SWCNT-Ni was defined by the area between SWCNT and nickel layer, and the thickness of SWCNT-Ni was the distance between SWCNT and nickel layer. The distance between SWCNT and nickel layer was 1.5 Å and the diameter of (6, 6) SWCNT was 8.14 Å. So the cross-sectional area of SWCNT-Ni was defined as 93.74 Å × Å. In the atomistic simulation, the stress ofσmn on the m and n directions is given by the following formula [17]:
(3)σmn=1NS∑imivimvinVi-12Vi∑j∂Urij∂rijrijmrijnrij,
where NS is the number of particles contained in the region S, and S is defined as the region of atomic interaction, vi is the velocity of the particle i, and rijm and rijn are two components of the vector from the atoms i to j. Vi is the volume assigned to atom i.For each atomic model from this work, the constant pressure and constant temperature (NPT) integration was performed using a Nose/Hoover temperature thermostat and a Nose/Hoover pressure barostat, and the bulk pressure was defined as zero. Subsequently, under the conditions of constant volume and constant temperature (NVT), these models were equilibrated to the desired temperatures using a Nose/Hoover temperature thermostat. The velocity Verlet method was used to integrate the classical equations of motion with a time step of 0.5 fs.During the entire stretching process shown in Figures1 and 2, the 47th and 48th rings of carbon atoms at the right end of both models were fixed, the first and second rings of carbon atoms at the left end were also restricted in their x- and y-directions, and the models of SWCNT and SWCNT-Ni were gradually stretched along the z-direction; the tensile load was applied along the z-direction. After every step of stretching, it was necessary that the simulation system provide sufficient time for the model to relax. In the modeling of the stretching process, the relaxation time was defined to be 10 ps. The external applied force was obtained by summing the internal forces acting on the atoms of the stretching section. All the data used to evaluate the mechanical properties of SWCNT and SWCNT-Ni was extracted from the back 1000 steps after the relaxation time of every instance of stretching. Considering unintentional sudden variations, the average of 10 sets of data at an interval of 100 steps during the above 1000 steps was used to determine the mechanical properties of SWCNT and SWCNT-Ni.


## 3. Results and Discussions

After the MD computations of the stretching process for SWCNT and SWCNT-Ni at the different temperatures of 220 K, 300 K, 500 K, 700 K, 1000 K, and 1200 K, the internal forces acting on the atoms and the stretching distance were obtained, which could be transformed into the corresponding stress and strain values. Due to the difficulty of defining the thickness of SWCNT-Ni, the definition of stress in the stretching process was from the atomistic stress of the MD modeling [17]. Figures 3 and 4 show all the simulated stress-strain curves over the prescribed temperature range for SWCNT and SWCNT-Ni, which clearly indicates that temperature was the crucial factor that determines the diversity of the stress-strain curves. From the curves in Figure 4, the stress-strain curves of SWCNT-Ni over the temperature range underwent some phase of linear elastic deformation, nonlinear elastic deformation, yield deformation, strengthening, and rupture, for which the length of any phase exhibited the temperature effect. Compared with SWCNT-Ni’s stress-strain curves in Figure 4, SWCNT’s stress-strain curves in Figure 3 were different. During the stretching process over all the prescribed temperatures, the stress-strain curves of the SWCNT did not exhibit any yielding phase. Unlike SWCNT’s stress-strain curves, the stress-strain curve of the SWCNT-Ni at 220 K clearly exhibited a yielding section after elastic deformation in Figure 4, which was also found in the other curves of SWCNT-Ni. In Figure 4, SWCNT-Ni’s stress-strain curves were different from those of SWCNT, which revealed that the nickel atoms coated onto the surface of SWCNT play an important role in affecting the mechanical behavior of SWCNT-Ni. In SWCNT-Ni’s stretching process at different temperatures, nickel atoms’ rearrangement can occur, which would change some mechanical properties of SWCNT-Ni. The strain-stress curve of the armchair (5, 5) carbon nanotubes with nickel-coating exhibited a flat phase at the beginning of the stretching at 300 K [11]; this phenomenon did not occur in the stress-strain curves shown in Figure 4, which might be due to the use of a different initial equilibrium method and the change of the cutoff function. In this work, NPT integration was performed using a Nose/Hoover temperature thermostat and Nose/Hoover pressure barostat, the bulk pressure was defined as zero, and the internal stress could be released fully.Figure 3
Stress-strain curves of SWCNT at different temperatures.Figure 4
Stress-strain curves of SWCNT-Ni at different temperatures.According to the force and deformation data from the MD computations for the stretching at different temperatures, the various curves, such as the failure force and temperature, the failure stress and temperature, the failure strain and temperature, and Young’s modulus and temperature, could be analyzed, thus enabling the mechanical properties of SWCNT and SWCNT-Ni to be determined at the different temperatures. Figures5, 6, 7, and 8 show the failure force, failure stress, failure strain, and Young’s modulus, respectively, at temperatures of 220 K, 300 K, 500 K, 700 K, 1000 K, and 1200 K. In the failure force and temperature curves in Figure 5, the failure force of SWCNT and SWCNT-Ni decreases as the temperature increases. In the same way, this tendency of decreasing as the temperature increases was observed in the curves of failure stress, failure strain, and Young’s modulus in Figures 6, 7, and 8, respectively.Figure 5
Failure force versus temperature curves of SWCNT and SWCNT-Ni.Figure 6
Failure stress versus temperature curves of SWCNT and SWCNT-Ni.Figure 7
Failure strain versus temperature curves of SWCNT and SWCNT-Ni.Figure 8
Young’s modulus versus temperature curves of SWCNT and SWCNT-Ni.As shown in Figure6, the failure stress of SWCNT decreases more rapidly than that of SWCNT-Ni as the temperature increases. In the MD computation in this work, the attenuation distance R1 of the cutoff function increased and the failure stress of SWCNT was still higher than the experimental values that range from 11 GPa to 63 GPa at room temperature [18], which resulted from the perfect SWCNT used in the molecular dynamics modeling in this work. At the same time, the failure stress at different temperatures in Figure 6 was still larger than that of the modeling for a (5, 5) pristine SWCNT using a modified Tersoff-Brenner potential [19], which could be caused by differences in the SWCNT modeled. In Figure 6, the failure stress of SWCNT-Ni was nearly half that of SWCNT over the same temperature range, which was also found in a previous study [11], and it could be attributed partly to the increase of the cross-sectional diameter and the decrease of the failure force of the SWCNT-Ni. Compared with SWCNT, the nickel atoms on the surface of SWCNT-Ni could weaken the strength of SWCNT. When being stretched, the six-member ring of SWCNT-Ni suffered local distortion, the C-C bond length became longer along the z-axis (tensile direction), and the C-C interaction became equal to the C-Ni interaction, which enabled thermal fluctuation to be a key factor for determining the rupture of SWCNT-Ni. In this work, the MD simulation of SWCNT-Ni was different from those using the pristine Brenner potential and the dissimilar attenuation distance [11].For the failure strains shown in Figure7, the failure strain of SWCNT was still larger than the experimental results of 10–13% failure strain [18] and 25% smaller than those results [20]. Figure 7 clearly shows that the failure strain of SWCNT and SWCNT-Ni would decrease almost linearly with temperature, thus indicating that the effect of temperature is dominant when the temperature increases. During the stretching processes, carbon atoms could prefer to exhibit a sp2 structure and tend to maintain this structure at the Brenner potential [21]. For the stretching of SWCNT-Ni, once a particular C-C bond is broken, the nearby nickel atoms would terminate the dangling bonds of the carbon atoms and form bonds with each other [11]. Unlike the failure stress and failure strain shown in Figures 6 and 7, Young’s modulus of SWCNT and SWCNT-Ni was not severely sensitive to temperature. In Figure 8, the variation of Young’s modulus of SWCNT was approximately 6% from 220 K to 1000 K, which was consistent with the results of [22]; however, when the temperature increased to 1200 K, Young’s modulus of SWCNT decreased by approximately 10.6%. In Figure 8, all Young’s modulus values of SWCNT at temperatures from 220 K to 1200 K were larger than the results of [11]. However, compared to those simulations using the Brenner potential function [22] and the modified Morse potential function [15, 23], in this work, Young’s modulus of SWCNT was in the range of the computed and experimental values of [18]. After coating nickel onto the surface of SWCNT, Young’s modulus of SWCNT-Ni decreased over the temperature range from 220 to 1200 K by approximately 27.6%. Apparently, the sectional diameter of SWCNT-Ni was larger than that of SWCNT, and the C-Ni interactions decreased the C-C bond strength, thereby affecting the mechanical behavior of SWCNT-Ni.When SWCNT or SWCNT-Ni was stretched to some extent, rupturing failure was inevitable. As shown in Figure9, all the final rupture surfaces at different temperatures were found to be completely brittle, in agreement with the past report [15]. From the stretching rupture surfaces shown in Figure 9, the rupture surfaces of SWCNT-Ni at different temperatures were completely separated, while the rupture surfaces of SWCNT were slightly connected by a single chain of carbon atoms after breaking failure; similar rupture surfaces were also noted in the past [24]. During the stretching simulation for SWCNT and SWCNT-Ni, when breakage of the C-C bonds occurred, a significant amount of elastic potential energy was released. As the temperature increased, nickel atoms would obtain enough energy to overcome the energy barrier, thereby breaking the nearby C-C bonds, followed by the occurrence of rupture, with the rupture surfaces 45° from the z-direction, as shown Figure 9.Figure 9
Rupture surfaces of SWCNT and SWCNT-Ni at different temperatures.In the MD simulations of stretching SWCNT and SWCNT-Ni at different temperature ranging from 220 to 1200 K, the initial working length, that is, the initial equilibrium total length minus the length of the fixed layer of carbon atoms, could also indicate the effect of temperature on thermal expansion behaviors. The initial working lengths of SWCNT and SWCNT-Ni shown in Figure10 indicated that the working length of SWCNT-Ni was longer than that of SWCNT at the same temperature, the equilibrium bond length of the nickel atoms increased with temperature, and the SWCNT-Ni had a larger thermal expansion coefficient; therefore, SWCNT-Ni was more temperature sensitive than SWCNT. At higher temperature during stretching, the C-C, C-Ni and Ni-Ni interactions were weakened, leading to the diversification of the mechanical properties, while the nickel atoms on the surface of the SWCNT-Ni played a role in interfering with the C-C bonds at a variety of temperatures.Figure 10
Initial working lengths of SWCNT and SWCNT-Ni at different temperatures.


## 4. Conclusions

In this work, using the second-generation Brenner reactive empirical bond-order potential with a modified attenuation distance of the attenuation function of the Brenner potential, MD simulations of stretching the armchair (6, 6) SWCNT and the SWCNT-Ni were conducted at different temperatures of 220 K, 300 K, 500 K, 700 K, 1000 K, and 1200 K. The simulation results revealed that the mechanical properties of both SWCNT and SWCNT-Ni were influenced by temperature. In particular, the failure stress and failure strain of SWCNT were obviously reduced at higher temperature, but Young’s modulus of SWCNT was not very sensitive to temperature; however, when the temperature was above 1000 K, Young’s modulus became greatly reduced. Compared with SWCNT, Young’s modulus of SWCNT-Ni was lower at the same temperature. Meanwhile, the failure force and the failure stress of SWCNT and SWCNT-Ni were greatly weakened at higher temperatures. When the surface of SWCNT was coated with nickel atoms, the nickel atoms tended to bind to the carbon atoms. With increasing temperature, nickel atoms on the surface of SWCNT-Ni became more active, and the influence of the nickel interference with the C-C bond became more significant, thereby leading to increased temperature sensitivity of the mechanical properties of SWCNT-Ni.



---

*Source: 767182-2015-04-27.xml*