Journal of Nanomaterials

Volume 2013 (2013), Article ID 386068, 4 pages

http://dx.doi.org/10.1155/2013/386068

Research Article

Mass Transport Induced by Heat Current in Carbon Nanotubes

Wei-Rong Zhong,¹ Zhi-Cheng Xu,¹ Ming-Ming Yang,¹ and Bao-Quan Ai²

¹Department of Physics and Siyuan Laboratory, College of Science and Engineering, Jinan University, Guangzhou 510632, China

²Laboratory of Quantum Engineering and Quantum Materials, ICMP and SPTE, South China Normal University, Guangzhou 510006, China

Received 26 September 2013; Revised 15 November 2013; Accepted 22 November 2013

Academic Editor: Miguel A. Correa-Duarte

Copyright © 2013 Wei-Rong Zhong 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.

Transport of helium atoms in the carbon nanotubes is investigated in the presence of temperature gradients. The heat current flowing along the carbon nanotubes can induce a stable directed transport of helium; it is demonstrated that the heat current density rather than the temperature gradient performs as a fundamental physical factor to the mass transport. We provide an alternative route to control the mass transport by using heat. Our results reported here are also relevant for understanding the transition from thermal energy to mechanical energy.

1. Introduction

The study of transport is always a challenge issue with incalculable value to us. Thermal and mass transports are the two important transports and have attracted substantial concerns [1]. In thermal transport, thermal conduction is the transfer of heat energy by microscopic diffusion and collisions of particles or quasiparticles within a body due to a temperature gradient. The mass transport refers to the orientated flow of particles. Even if thermal transport is a deviation from mass transport in physical principle, researchers are always trying to find out their inherent relationship [2, 3].

Recently, an experimental result shows that thermal gradients can drive the subnanometer motion of cargoes along carbon nanotubes (CNTs) [4]. Another numerical result shows that the fullerene encapsulated in carbon nanotubes can also be driven by temperature gradients [5]. More and more studies focus on nanoparticle manipulation using heat [6–8]. A model of carbon nanotubes filled with fullerene (C60) reported that the temperature can pumping large particles [9]. The Langevin model of a thermally -driven double-walled nanotubes motor has described the dynamics of nanoelectromechanical device activated by heat [10]. It was also reported that the velocity of C60 induced by temperature gradients had a linear relationship with heat flux [6]. These few previous efforts suggest that the heat originated from temperature gradients can induce the motion of particles. However, the underline physics of mass transport driven by heat is still not so clear.

In this paper, we will study the carbon nanotubes filled with helium (He) atoms and investigate the transport of He atoms induced by the heat flowing through the CNTs. By constructing the relationship between the mass flux and the heat current, we try to reveal some fundamental physics about the energy transport and its transformation. The obtained results are of significance for understanding the thermally induced mass transport in quasi-one-dimensional materials and may be useful for controlling the mass transport.

Figure 1 shows the carbon nanotubes filled with He atoms. Here we use He atoms as the fillings just because of their stable chemical properties. Two heat baths (with temperature and , resp.) are connected to the two sides of the CNTs. The two ends of the CNTs are fixed to avoid their spurious global rotation. The fixed region and the heat baths occupy one layer and six layers of atoms, respectively. In order to keep a stable number of helium atoms, the periodic boundary is applied for He atoms along the axial direction of the CNTs. That means a helium atom will be inset in one end of the CNTs when another Helium atom goes out of another end of the CNTs. In our simulations, we have used classical molecular dynamics method based on the Tersoff-Brenner potential [11] of carbon-carbon (C–C) bonding interactions. The Lennard-Jones (LJ) potentials form [12, 13]

is used for van der Waals interactions of He–He and He–C couplings. The parameters of LJ are and for He–He interaction and and for He–C interaction, respectively. The equations of motion for atoms in either the left or right Nosé-Hoover thermostat are [14, 15]

where is the momentum and is the force applied on the th atom; , where is the relaxation time, which is set to . is the dynamic parameter of the thermostat; is the instant temperature of the heat baths at time , which is defined as , where is the time-dependent velocity; ( or ) is the set temperature of the heat baths. is the number of the atoms in the heat baths, is the Boltzmann constants and is the mass of the carbon atom. The set temperatures of the heat baths, and , are placed in the two ends of the CNTs and the temperature gradient is denoted by , where is the length of the CNTs. For the convenience of comparison, we use the average temperature 300 K, which is near the room temperature. The atoms between the two thermostats follow Newton’s law of motion:

where runs over all of the atoms between the two thermostats.

Figure 1: Diagram of the carbon nanotubes (filled with helium atoms, the green balls) coupling to two heat baths (the red and blue balls). The temperatures of the red part and blue part are and , respectively. As , the heat current will flow from the left side to the right side. Accordingly, a mass flux induced by the temperature gradient will occur along the axial direction of the CNTs. The length of CNTs is 12 nm.

We integrate these equations of motion by Verlet method [16]. The time step is 0.55 fs. Generally, as shown in Figure 1, the heat current can stabilize around the set value when the temperature gradients are applied to the CNTs after 2 ns. At the same time, the mass flux driven by heat current occurs along the axial direction of CNTs. The heat bath acts on the particle with a force ; thus the power of heat bath is , which can also be regarded as the heat flux that comes out of the high temperature heat bath and is injected into the low temperature heat bath. The total heat flux injected from the heat bath to the system can be obtained by , where the subscript runs over all the particles in the thermostat [17]. The average velocity of He is obtained by , where is the velocity of the th helium atom, means a time average, and is the total number of He atoms. The mass flux is defined as

where is the mass of single He atom and is the mass density of He in CNTs.

When applying the temperature gradients on the CNTs, a stable flowing of He atoms along CNTs will occur. Figure 2 shows the temperature gradient dependence of average velocity of He atoms for different diameters of CNTs. To consider the diameter dependence of CNTs on the mass transport, here we study armchair chirality of the CNTs (5, 5), (7, 7), (10, 10), and (12, 12). Obviously, the average velocity increases with temperature gradient. However, the average velocity of He atoms depends on the diameter of CNTs with slight irregularities.

Figure 2: Temperature gradient dependence of average velocity of He atoms in armchair CNTs for different diameters of the CNTs. The mass densities of He in (5, 5), (7, 7), (10, 10), and (12, 12) are 66.39, 33.89, 16.6, and 16.6 amu/nm³, respectively.

From our model, when applying temperature gradient on CNTs, the motion of He atoms occurs with a heat current flow from the hot heat baths to the cold heat baths along the surface of CNTs. If one uses the heat current, which is corresponds to the temperature gradient, to plot the -axis, one would find an organized result comparing to Figure 2. As shown in Figure 3, if the heat current flowing through the CNTs is the same, the CNTs with chirality of (5, 5) can drive higher average velocity of He atoms than that of (7, 7). When the diameter of armchair CNTs increases, the ability to drive the motion of He atom reduces. Then it is not difficult to find the conflicts between temperature gradients and heat current. In order to solve the conflicts, we plot the heat current density as the transverse axis and the average velocity as the longitudinal axis.

Figure 3: Heat current dependence of average velocity of He atoms in armchair CNTs for different diameters of CNTs. The CNTs are the same as in Figure 2.

As illustrated in Figure 4, the curves for different diameters as shown in Figure 3 all run together and the average velocity in the CNTs is proportional to the heat current density via the CNTs. This relationship is independent of the diameter as well as the length of CNTs. In Figure 4, we also provide the results of two zigzag CNTs, (20, 0) and (16, 0). When the driving heat current density ranges from about −2 to 2 uW/nm², the linear relationship between average velocity and heat current density is still valid. This indicates that the heat current density determines the average velocity, which may be independent of the structure of the CNTs. Since the mass flux is obtained from (3), which means that the mass flux does not depend on the cross section area, then the mass flux is proportional to the average velocity as well as the heat current density.

Figure 4: Average velocity versus heat current density for different chiralities, diameter, and lengths of CNTs. The heat current density is obtained by , where is the total heat flux injected from the heat bath to the system and is the area of cross section of CNTs, respectively. The mass density of He in (16,0) and (20, 0) is 16.6 amu/nm³.

In summary, we have studied the temperature gradient driven mass flux in the CNTs using nonequilibrium molecular dynamics. It is reported that relationship between mass flux and heat current density rather than average velocity versus temperature gradients can describe the fundamental physics of temperature gradient induced mass transport better. Our results try to connect energy transport and mass transport from a new point of view, which would help researchers understand the physical mechanism of transport. Although we demonstrate a universal relationship between heat and motion, more physical properties such as the transfer efficiency need to be addressed in detail. Even so, our results are an illumination to control mass flux by using heat conduction.

Acknowledgments

The authors would like to thank Siyuan clusters for running a part of our programs. Wei-Rong Zhong thanks Y. Xiao at City University of Hong Kong for his contribution in plotting. This work was supported in part by the National Natural Science Foundation of China (Grant nos. 11004082, 11175067, and 11291240477); the Natural Science Foundation of Guangdong Province, China (Grant nos. 10451063201005249 and S201101000332); the PCSIRT (Grant no. IRT1243); and the Fundamental Research Funds for the Central Universities, JNU (Grant no. 21611437).

References

K. F. Rinne, S. Gekle, D. J. Bonthuis, and R. R. Netz, “Nanoscale pumping of water by AC electric fields,” Nano Letters, vol. 12, no. 4, pp. 1780–1783, 2012.

M. J. Longhurst and N. Quirke, “Temperature-driven pumping of fluid through single-walled carbon nanotubes,” Nano Letters, vol. 7, no. 11, pp. 3324–3328, 2007.

H. Somada, K. Hirahara, S. Akita, and Y. Nakayama, “A molecular linear motor consisting of carbon nanotubes,” Nano Letters, vol. 9, no. 1, pp. 62–65, 2009.

A. Barreiro, R. Rurali, E. R. Hernández et al., “Subnanometer motion of cargoes driven by thermal gradients along carbon nanotubes,” Science, vol. 320, no. 5877, pp. 775–778, 2008.

R. Rurali and E. R. Hernández, “Thermally induced directed motion of fullerene clusters encapsulated in carbon nanotubes,” Chemical Physics Letters, vol. 497, no. 1–3, pp. 62–65, 2010.

N. Wei, H.-Q. Wang, and J.-C. Zheng, “Nanoparticle manipulation by thermal gradient,” Nanoscale Research Letters, vol. 7, article 154, 2012.

H. A. Zambrano, J. H. Walther, and R. L. Jaffe, “Thermally driven molecular linear motors: a molecular dynamics study,” Journal of Chemical Physics, vol. 131, Article ID 241104, 2009.

I. Hanasaki, A. Nakatani, and H. Kitagawa, “Molecular dynamics study of Ar flow and He flow inside carbon nanotube junction as a molecular nozzle and diffuser,” Science and Technology of Advanced Materials, vol. 5, no. 1-2, pp. 107–113, 2004.

X. wang and J. D. Lee, “Electronic transport in a new type nano-junction: carbon atomic chain inserted into a carbon nanotube,” Journal of Computational and Theoretical Nanoscience, vol. 9, pp. 1–5, 2012.

S. Holek, D. Reguera, and J. M. Rubi, “Carbon-nanotube-based motor driven by a thermal gradient,” The Journal of Physical Chemistry C, vol. 117, pp. 3109–3113, 2013.

D. W. Brenner, “Empirical potential for hydrocarbons for use in simulating the chemical vapor deposition of diamond films,” Physical Review B, vol. 42, no. 15, pp. 9458–9471, 1990.

R. E. Tuzun, D. W. Noid, B. G. Sumpter, and R. C. Merkle, “Dynamics of fluid flow inside carbon nanotubes,” Nanotechnology, vol. 7, no. 3, pp. 241–246, 1996.

R. E. Tuzun, D. W. Noid, B. G. Sumpter, and R. C. Merkle, “Dynamics of He/C60 flow inside carbon nanotubes,” Nanotechnology, vol. 8, no. 3, pp. 112–118, 1997.

J. Hu, X. Ruan, and Y. P. Chen, “Thermal conductivity and thermal rectification in graphene nanoribbons: a molecular dynamics study,” Nano Letters, vol. 9, no. 7, pp. 2730–2735, 2009.

G. Wu and B. Li, “Thermal rectification in carbon nanotube intramolecular junctions: molecular dynamics calculations,” Physical Review B, vol. 76, no. 8, Article ID 085424, 8 pages, 2007.

H. Rai-Tabar, Computational Physics of Carbon Nanotubes, Cambridge University Press, New York, NY, USA, 2008.

W. R. Zhong, M. P. Zhang, B. Q. Ai, and D. Q. Zheng, “Chirality and thickness-dependent thermal conductivity of few-layer graphene: a molecular dynamics study,” Applied Physics Letters, vol. 98, Article ID 113107, 2011.