Abstract

By using molecular dynamics method, thermal conductivity of (10, 10) carbon and boron nitride (BN) nanotubes under radial compression was investigated, and the 𝜆-𝑇 (thermal conductivity versus temperature) curves of the two nanotubes were obtained. It is found that with the increase of temperature the thermal conductivity of two nanotubes decreases; the nanotubes, under both the local compression and whole compression, have lower thermal conductivity, and the larger the compressive deformation is, the lower the thermal conductivity is; the whole compression has more remarkable effect on thermal conductivity than the local compression.

1. Introduction

As the feature sizes of MEMSs (microelectro mechanical systems) continue to be scaled down, the problem of thermal dissipation becomes even more important. It has been found that very few materials are known to exhibit high thermal conductivity at reduced dimensions, and traditional semiconductors like Si have a dramatic suppression of thermal conductivity at the nanoscale [1].

Just like carbon nanotubes can be viewed as rolling a graphite sheet into a cylinder, boron nitride (BN) tubes can also be regarded as geometrically derived by rolling a hexagonal BN sheet into a cylinder. Due to low-dimensional effects, carbon and BN nanotubes have been predicted to have unusual thermal properties [2, 3]. The thermal conductivity of an isolated carbon nanotube has been predicted to reach 1800–6000 W/mK at room temperature [4]. Similarly, the thermal conductivity (290 K) of isolated BN nanotube has been reported to be 1620 W/mK [5].

Recent Xu and Buehler [6] carried out a systematic investigation of the effects of axial tensile, compressive and torsional strain on the thermal conductivity of single-walled carbon nanotubes using MD simulation, and found that the thermal conductivity drops when the different strain is applied to carbon nanotube, respectively. However, up to now the work about the effects of radial compression on carbon nanotube as well as BN nanotube is not found. Here we perform nonequilibrium MD simulation to investigate the thermal conductivity of the (10, 10) carbon and BN nanotubes under radial compression and discuss the effects of the local and whole compressive deformation on thermal conductivity of the carbon and BN nanotubes.

2. The Investigated Objects and Method

2.1. The Investigated Objects

Figure 1 shows the investigated objects, (10, 10) zigzag carbon and BN nanotubes. Considering that within the length range of 1 nm to 25 nm nanotubes have little change in thermal conductivity [7], the present tube length is taken to be 8.85 nm to 9.81 nm. The (10, 10) carbon tube is generated by the nanotube generator [8], a program module developed by us. The (10, 10) BN tube is constructed by alternately replacing the carbon atoms of (10, 10) carbon tube with boron or nitrogen atoms. After geometrically optimizing the carbon and BN nanotubes through MD simulation at 270 K, the initial configurations are obtained. The initial diameter 𝐷 of the carbon and BN nanotubes is 1.36 nm and 1.52 nm, respectively.

2.2. Method

In classical MD, atoms are regarded as basic particles, and the temperature of the simulated system is given through the equipartition theorem of energy. When the practical temperature is lower than the Debye temperature of material, the thermal capacity will be highly correlative with the temperature of system, the equipartition theorem of energy will become invalidated [9], and the quantum correction to temperature must be performed. In the present MD calculations of thermal conductivity, the temperature of both the nanotube systems is taken above their Debye temperature so that quantum effects can be ignored. The Debye temperature of the carbon and BN nanotubes takes 420 K [10] and 400 K [11], respectively.

In this study the EMD (equilibrium molecular dynamics) method with the time step of 0.5 fs is applied to geometrically optimize the carbon and BN nanotubes, as well as to calculate the thermal conductivity of the carbon and BN nanotubes locally or wholly compressed in radial direction (see Figure 1). The process to calculate the thermal conductivity can be described as the following: firstly relaxing the tubes in the NTV ensemble for 3 × 105 time-steps, then locally or wholly compressing the carbon and BN nanotubes in the radial direction through two parallel rigid plates, where the loading rate takes 0.03 nm per 1000 time steps, and the temperature is controlled at 270 K through Nose method; fixing the two rigid plates when the carbon and BN tubes reach certain deformation, calculating the heat flux of each time step for the subsequent 4 × 106 time-steps, and obtaining the thermal conductivity of the compressed tubes through integrating the autocorrelation function of heat flux [12].

In the present MD simulations, the Velocity-Verlet algorithm is taken. Tersoff potential, one three-body potential based on the concept of bond order, is used to describe the interaction between atoms. The Tersoff potential explicitly includes an angular contribution of the force and has been widely used in various applications for III–V elements. The Tersoff potential and its parameters for the carbon, BN atomic systems see [13, 14]. In the simulations of the local compression the width of two rigid-plates takes 0.74 nm. The interaction between atoms and the rigid plates is described using the L-J potential:𝜎𝑈(𝑟)=4𝜀𝑟12𝜎𝑟6,(1) where 𝑟 is the distance between atoms and rigid-plate; the constants 𝜀 and 𝜎 take 0.440 kJ and 0.385 nm [15], respectively.

In our EMD simulations, based on the Green-Kubo linear response theory, the thermal conductivity has the following formula [16]1𝜆=𝐾𝐵𝑇2𝑉0+𝑞(𝑡)𝑞(0)𝑑𝑡,(2) where 𝐾𝐵 is the Boltzmann constant, 𝑇 the simulative temperature, 𝑉 the volume of simulated system, that is, the cross-section area (=𝜋𝐷2/4) multiplying the tube length, 𝑡 time, and 𝑞 the effective heat flux.

The effective heat flux 𝑞 has the following expression [11] of:𝑑𝑞(𝑡)=𝑑𝑡𝑖𝑟𝑖(𝑡)𝑒𝑖(𝑡),thatis,𝑞(𝑡)=𝑖𝑣𝑖𝑒𝑖+𝑖𝑟𝑖𝑑𝑒𝑖,𝑑𝑡(3) where 𝑣𝑖 is the velocity of the 𝑖th atom, 𝑒𝑖 the total energy of the 𝑖th atom, and 𝑟𝑖 the coordinate.

3. Results and Discussion

It is proved that the contribution to thermal conductivity of nanotube mainly comes from two parts, that is, (1) the interaction between electrons and phonons, that is dependent on electronic structures (molecular configuration) and phonon scattering, and (2) the interaction between phonons, that is, correlative to atomic vibration.

Figure 2 gives the 𝜆-𝑇 (thermal conductivity versus temperature) curves of the carbon nanotube under different local compressive strain and whole compressive strain 𝜀, and Figure 3 gives those of the BN nanotube. For convenience to discuss, Figure 4 presents the 𝜆-𝑇 curves of the carbon nanotube under different compression modes, that is, the local compression and whole compression, and Figure 5 presents those of the BN nanotube. In these figures the compressive strain 𝜀 is defined as the ratio of the radial compressive deformation to the initial tube diameter 𝐷.

From Figures 2 and 3, it can be seen that,(1)with the increase of temperature, the thermal conductivity of both the carbon and BN nanotubes decreases, and when temperature >1,800 K, their thermal conductivity trends to be constant. What is the reason? It can be explained as the following: due to the increase of temperature, the U scattering process of phonons increases, and the average free range of phonons decreases, which causes the decreases in thermal conductivity of both the carbon and BN nanotubes;(2)under the local compression, both the carbon and BN nanotubes have lower thermal conductivity; the larger the compressive strain (deformation) is, the lower the thermal conductivity is. The phenomenon can be explained as the following: firstly, the local compression causes local cave-in on the tube wall, that is, the change in local molecular structure, and further the change of the local atoms in vibration mode; secondly, the local cave-in on the tube wall changes the heat exchange path and increases the local phonon U scattering. In a word, the local radial compression makes a bottleneck of heat exchange on the nanotube wall; the larger the compressive deformation is, the worse the thermal conductivity becomes;(3)at the same temperature, the whole compression also makes the thermal conductivity of both the nanotubes decrease; the larger the compressive deformation is, the lower the thermal conductivity is. The reason lies in that under whole compression the chemical-bonds are tightened up, and the larger the compressive deformation is, the tighter the chemical-bonds become, which results in the change of atoms in vibration mode, that is, vibration amplitude and frequency, as well as the decrease in thermal conductivity;(4)at relatively low temperature, the effect of both the local and whole compression on the thermal conductivity of the carbon and BN nanotubes is remarkable. At high temperature, the effect becomes slight, which is because that the effect of compressive deformation has been inundated by that of high temperature.

According to Figures 4 and 5, it can be observed that, under the same temperature and compressive-deformation, the wholly compressed carbon (or BN) nanotube has lower thermal conductivity than the locally compressed one. The phenomenon can be explained through the different heat-transfer mechanism of the nanotubes under different compression modes, that is, the tightened up chemical-bonds in the wholly compressed nanotubes results in the change of atom vibration mode as well as the change of thermal conductivity; the bottleneck of heat exchange on the tube wall decreases the thermal conductivity of the locally compressed nanotubes.

The present (10, 10) carbon nanotube without compressive deformation is predicted to have the thermal conductivity of 180 W/mK and 94 W/mK, respectively, at 800 K and 1200 K. The corresponding MD results obtained in [17] are 190 W/mK and 120 W/mK, respectively. Our results are close to those of [17], which implies the effectivity of the MD method in the present paper. In [18, 19], the measurement of the thermal conductivity of single-layer graphenes was reported, and the room temperature values of the thermal conductivity were found to be in the range 3000 to 5300 W/mK. Conceptually, single-wall carbon nanotubes can be considered to be formed by the rolling of a single-layer of graphite/graphene. By the cubic polynomial interpolation of the curves in Figure 2, it is found that the thermal conductivity of the (10, 10) carbon nanotube without deformation is about 430 W/mK at 300 K, lower by one order of magnitude than that of graphenes in [18, 19].

4. Conclusions

As the important components of future NEMSs (nanomechanical systems), one-dimensional nano-materials like carbon and BN nanotubes have broad application prospect, and their thermal conductivity and thermal dispersion will directly affect the system operation. In this study the EMD method is used to investigate the effect of local and whole radial compression on the thermal conductivity of carbon and BN nanotubes. It is found that both the local and whole compression can decrease the thermal conductivity of nanotubes. The conclusion will be very helpful for the design of carbon and BN nanotube thermal-conductivity-devices.