Negative Thermophoresis in Concentric Carbon Nanotube

Sep 14, 2016 - Nanosurfactants create droplet-sized reaction flasks. By adding a suit of nanoparticle surfactant armor to otherwise dull droplets, che...
0 downloads 15 Views 1MB Size
Subscriber access provided by Northern Illinois University

Communication

Negative thermophoresis in concentric carbon nanotube nanodevices Jiantao Leng, Zhengrong Guo, Hongwei Zhang, Tienchong Chang, Xingming Guo, and Huajian Gao Nano Lett., Just Accepted Manuscript • DOI: 10.1021/acs.nanolett.6b02815 • Publication Date (Web): 14 Sep 2016 Downloaded from http://pubs.acs.org on September 16, 2016

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Nano Letters is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 24

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

Negative thermophoresis in concentric carbon nanotube nanodevices Jiantao Leng1, Zhengrong Guo1, Hongwei Zhang1, Tienchong Chang1,* , Xingming Guo1, Huajian Gao2 1

Shanghai Institute of Applied Mathematics and Mechanics, Shanghai Key Laboratory of Mechanics in Energy Engineering, Shanghai University, Shanghai 200072, People’s Republic of China 2

School of Engineering, Brown University, Providence, RI 02912, USA

ABSTRACT Positive and negative thermophoresis in fluids has found widespread applications from mass transport to molecule manipulation. In solids, although positive thermophoresis has been recently discovered in both theoretical and experimental studies, negative thermophoresis has never been reported. Here we reveal via molecular dynamics simulations that negative thermophoresis does exist in solids. We consider the motion of a single walled carbon nanotube nested inside of two separate outer tubes held at different temperatures. It is found that a sufficiently long inner tube will undergo negative thermophoresis, whereas positive thermophoresis is favorable for a relatively short inner tube. Mechanisms for the observed positive thermophoresis and negative thermophoresis are shown to be totally different. In positive thermophoresis, the driving force comes mainly from the thermally induced edge force and the interlayer attraction force, while the driving force for negative thermophoresis is mainly due to the thermal gradient force. These findings have enriched our knowledge of the fundamental driving mechanisms for thermophoresis in solids and may stimulate its further applications in nanotechnology. *

Email address: [email protected].

ACS Paragon Plus Environment

Nano Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Keywords: Negative thermophoresis, thermal gradient force, edge force, interlayer attraction force, molecular dynamics simulations

ACS Paragon Plus Environment

Page 2 of 24

Page 3 of 24

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

The discovery1 and understanding2 of thermophoresis have provided ample opportunities for applications including thermal precipitation,3, energy conversion,8,

9

4

particle separation,5,

6

drug delivery,7

biomolecule manipulation,10-14 and interaction detection,15 etc.

Thermophoresis was first discovered in liquids in 1856 by Ludwig1 who observed salt diffusing in water against a temperature gradient. The phenomenon was confirmed by Soret16 who believed that the temperature gradient leads to a homogeneous salt concentration gradient in liquid. Thermophoresis in aerosols was first reported in 1870 by Tyndall17 who found that dust particles can be pushed away from a heated surface in a dust-filled room. However, in spite of the wide applications3-7, 10-14, 18 of thermophoresis in fluids, it was not until more than a century later that thermophoresis at solid-solid interfaces was first predicted. In 2006, Schoen et al.19, 20 demonstrated via molecular dynamics (MD) simulations that an axial temperature gradient along a single-walled carbon nanotube (SWCNT) would drive a confined nanoparticle to move toward the direction of decreasing temperature. The phenomenon was immediately confirmed in an elegant experimental study by Barreiro et al.9 Such a temperature gradient can also induce the motion of many kinds of nanoscale particles, including fullerene clusters,21 capsule-like nanotubes,22 graphene sheet23 and water droplets,24-26 indicating that thermophoresis in solids has the potential of opening up new opportunities in nanoscale mass transport and energy conversion technologies.9 In contrast, there also exists negative thermophoresis, i.e., induced motion towards a direction of increasing temperature. Negative thermophoresis in fluids was first noticed in 1967 by Dwyer27 in a theoretical solution, and the name was coined by Sone.28 The phenomenon was soon confirmed by a series of numerical calculations.29-31 Very recently,

ACS Paragon Plus Environment

Nano Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

negative thermophoresis in fluid was observed in experiments.11, 32, 33 However, to this date, the existence of negative thermophoresis in solids still has not been verified in either theoretical or experimental studies, despite its potential applications in nanotechnology. Here, we demonstrate via MD simulations that negative thermophoresis can occur at solid-solid interfaces. By investigating the motion of a SWCNT nested inside of two separate outer tubes held at different temperatures, we find that the direction of motion of the inner tube is dependent on the geometry of the device. A sufficiently long inner tube with both ends extending out of the outer tubes will undergo negative thermophoresis, while positive thermophoresis is favorable for a relatively short inner tube. Mechanisms for the observed geometry dependent thermophoresis will be discussed in detail.

Simulation models and methods Two models (i.e., Model I and Model II, as shown in Figure 1) are considered in our simulations. In both models, a narrow inner SWCNT is commensurately nested in two open-ended coaxial CNTs. The only difference between the structures of the two models is that the inner tube in Model II is substantially longer and extends out of the outer tubes, while that in Mode I is relatively short with both ends confined inside the outer tubes. Unless specifically stated, the inner tube is a (5, 5) tube, while the two outer tubes are (10, 10) tubes held at 300 and 500 K. All simulations are carried out using the LAMMPS package.34 A time step of 0.5 fs is used in all simulations. A reactive empirical bond order (REBO) potential35 is used to describe the in-plane C-C bond interactions and the interlayer vdW interactions are described by a

ACS Paragon Plus Environment

Page 4 of 24

Page 5 of 24

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

Lennard-Jones (L-J) 6-12 potential (with a well-depth of ɛ = 2.968 meV and an equilibrium distance of σ = 0.3407 nm). We note that the qualitative behavior of the inner tube in the present system should be insensitive to the choice of atomic potentials because the relative sliding behavior of concentric CNTs is governed by the structure as well as the relative strengths of the intertube and intratube interactions, as Tangney et al.36, 37 showed in studying the interlayer friction between concentric CNTs.

Figure 1. Models for molecular dynamics simulations. An inner SWCNT (whose length is different in the two models) is nested inside two concentric outer SWCNTs (whose configurations are exactly the same in the two models). In Model I, the inner nanotube is 36.83 nm long, while that in Model II is 99.28 nm long. The two outer nanotubes, each being 19.35 nm long, are symmetrically placed relative to the inner nanotube, with their neighboring ends spaced at 15 nm apart. In the simulation, the temperature of the left outer tube is lower than that of the right tube.

The shrink-wrapped boundary condition34 is imposed in the longitudinal direction of the inner tube to allow its free motion in the outer tubes, while fixed boundary conditions are used in the other two perpendicular directions. Under these non-periodical boundary conditions, atoms do not interact across the boundary or move from one side of the box to the other. The two outer tubes are constrained against rotation and motion of their centers of mass

ACS Paragon Plus Environment

Nano Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

in any direction, and maintained at different temperatures via a Berendsen thermostat to form a stable temperature difference.38, 39 We first fix the longitudinal motion of the inner tube by removing the velocity of its center of mass in all directions, and relax the system for 1.7 ns in a NVT (canonical) ensemble (in which the system temperature is set to the mean temperature of the two outer tubes) to allow a stable temperature gradient to be developed along the inner tube. We then start the simulation in an NVE (microcanonical) ensemble by releasing the artificial constraint on the longitudinal motion of the inner tube.

Results Geometry dependent thermophoretic motion. When the simulations start, both of the inner tubes in the two models start to move spontaneously along the outer tubes at an accelerated speed (see Figures 2a and 2b for snapshots of the motion of (5,5)(10,10) devices). However, the initial directions of motion of the inner tubes in the two models are opposite to each other, one toward the cooler side whereas the other toward the hotter side. In Model I, the inner tube drifts toward the low temperature CNT (LTCNT), corresponding to positive thermophoresis, while in Model II, the inner tube drifts toward the high temperature CNT (HTCNT), i.e. negative thermophoresis. For convenience, we shall regard the initial direction of motion of the inner tube as its intended direction of motion. In both cases, the initially resting inner tube accelerates until its rear end hits the end of one of the outer CNTs and then bounces back and decelerates while moving in the opposite direction, much like a ball bouncing on the ground under gravity (Figures 2c and 2d). As soon as the motion of the inner tube in the opposite direction decelerates to zero speed, it starts to accelerate toward the

ACS Paragon Plus Environment

Page 6 of 24

Page 7 of 24

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

intended direction again. It is also seen from Figures 2c and 2d that the speed of the inner tube in both models depends on the diameter and the temperature difference. However, the initial moving direction is only determined by the length of the inner tube. In this letter, we will only focus our attention on the driving mechanism of the positive and negative thermophoretic motion, and will leave more detailed analyses of various system parameters (including length, diameter and temperature difference) to future work.

Figure 2. Snapshots of thermophoretic motion in Model I (a) and Model II (b) of (5,5)(10,10) devices. The displacements and velocities of the inner tubes in Model I and Model II for devices with different configurations under different temperature differences are shown in (c) and (d), respectively. The elargement of the circled area of the curve for the (5,5)(10,10) device under 300-500 K temperature difference in Figure 2d will be presented in Figure 6d to clearly show the trap behavior of the inner tube.

ACS Paragon Plus Environment

Nano Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The behavior of the inner tube clearly shows that there exists a persistent driving force for the inner tube to move toward the cooler region in Model I but toward the hotter region in Model II. This means that the direction of the driving force in the present device can be designed by tuning the length of the inner tube (or the distance between the outer tubes, as will be shown shortly). In particular, our results show that the directional motion toward the high temperature region (i.e., negative thermophoresis) can be achieved via the design of the geometry of a nanodevice, which is extremely important for the nanoscale actuation and energy conversion.

Mechanics of geometry dependent thermophoresis. Thermophoretic motion of a shorter carbon nanotube nested inside or outside a longer concentric tube with a temperature gradient has been extensively studied,9, 21, 22, 40, 41 and it is found that the motion is always toward cooler regions. For the present devices, different temperatures are imposed on the two outer tubes which provide driving forces for the inner tube, and the motion is dependent on the device geometry. This immediately raises the question about the driving mechanism for the observed geometry dependence of the thermophoretic motion. To understand where the driving force comes from, we show in Figures 3a-3e the longitudinal distribution of the lateral interlayer van der Waals force (the inerlayer van der Waals force projected onto the tube axis) exerted on the inner tube by the outer tubes when the inner tube is accelarated toward its intended direction of motion. In order to eliminate the effect of the potential corrugation, the value of the force is calculated by averaging the force on each ring over the time period in which the inner tube slides a lattice length (~0.246 nm) in

ACS Paragon Plus Environment

Page 8 of 24

Page 9 of 24

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

the outer tubes. Figure 3f defines the ring number on the inner tube. It is clearly seen that the inner tubes in the two devices experience different force distributions. To gain a deeper insight into the driving mechanism, we split the force distribution into six regions: from Region I to Region VI (Figures 3a and 3b). Despite the similar force distribution profiles of the two models in Regions II-V, the force distributions of the two models are quite different in Regions I and VI where the interlayer contacts in the two models are different. In Model I, the ends of the inner tube are in contact (near the equilibrium distance of van der Waals interactions) with the walls of the outer tubes, while in Model II, it is the ends of the outer tubes that are in contact with the wall of the inner tube. In what follows, we will discuss in detail the forces exerted on the inner tube by the outer tubes in these regions.

Figure 3. The distributions of lateral force on the inner tube in two models. (a) and (b) The overview of force distribution for Model I and Model II, respectively. The colored regions indicate the positions of the outer tubes. (c) and (d) Typical force distritubions in Regions II and V for Model I and Model II. (e) Force distribution near the ends of the inner tube in

ACS Paragon Plus Environment

Nano Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Model I. (f) Definition of the carbon ring number on the inner tube.

The difference in temperature between the two outer CNTs can generate a temperature gradient on the inner tube (typical temperature distributions along the inner tubes of the two models are shown in Figure 4). Such a temperature gradient may induce thermal gradient forces

9, 41-43

on the inner tube atoms (i.e., atoms in Regions II and V in Figures 3a and 3b)

whenever they interact with the outer tubes. In addition, these atoms in contact with the outer tubes would experience frictional forces when the inner tube is sliding against the outer tube.36, 40 The total force exerted on Region II or Region V can thus be expressed as Fg + f, with Fg and f being the gradient and frictional forces, respectively. In fact, there are also random forces exerted on the inner tube atoms from their collision with the outer tubes. However, since the random forces cannot contribute a net driving force for a persistent directional motion, here we exclude the effect of the random force on the driving force.

Figure 4. Temperature gradients developed on the inner tubes in Model I (a) and Model II (b). Results for the inner tubes on three positions where they are accelerated to the intended

ACS Paragon Plus Environment

Page 10 of 24

Page 11 of 24

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

direcition are presented.

Previous studies have shown that the thermal gradient force on a carbon nanotube nested inside or outside another one with a thermal gradient is toward the direction of decreasing temperature,9,

41-43

implying that the latter feels a reaction force toward the direction of

increasing temperature. This can also be considered as a result of the net current of phononic excitations 9, 44 induced by the tempearture gradient on the inner tube that hits the cooler outer tube and pulls the hotter outer tube. That is to say, the thermal gradient force in Model I provides a resistance force, while in Model II, it is a driving force. This is confirmed by our calculations. In Model I, the line densities of FII and FV are about -0.10 pN/nm, while in Model II, they are about 0.063 pN/nm. In other words, FII + FV provides a resistance force to the inner tube of about 1.81 pN in Model I, but a driving force of about 2.25 pN in Model II. Since the thermal gradient force is opposite to the sliding direction of the inner tube in Model I, there should be some other forces that drive the tube to move. Now we turn our attention to the interlayer attraction forces exerted on Regions III and IV in Model I (see Figure 3a) and those on Regions I, III, IV, VI in Model II (see Figure 3b), where the ends of the outer tubes contact with the wall of the inner tube. When an end of the inner tube extends out of the outer tube (see Figure 5), an interlayer attraction force (Fa) is induced between the end of the outer tube and the protruding part of the inner tube.36, 45-48 This means that the forces on Regions III and IV in Model I and those on Regions I, III, IV, VI in Model II, i.e., FIII and FIV in Figure 3a and FI, FIII, FIV, FVI in Figure 3b can be expressed as Fa + f.36, 40 Here we take eight rings of the inner tube atoms in contact

ACS Paragon Plus Environment

Nano Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 24

with each outer tube edge to account for the interlayer attraction forces exerted on the inner tube, as the edge attraction occurs generally within a range of three to four rings of atoms. In Model I, FIII + FIV = (-866.822 + 860.843) pN = -5.98 pN, which pushes the inner tube toward the cool region and is thus a driving force.

Outer tube

Fe

Inner tube

Fa

Figure 5. Illustration of the origin of the interlayer attraction force Fa on an inner tube when it extends out of an outer tube, and the origin of the edge force Fe when the end of the inner tube contacts with the wall of the outer tube.

Unlike in Model I, in Model II the interlayer attraction forces Fa from an outer tube appear in pair, approximately equal in magnitube but opposite in direction, and thus almost cancel each other and cannot generate an significant driving force. In the present calculation, FI + FIII + FIV + FVI = (868.09 – 867.67 + 862.01 – 861.55) pN = 0.88 pN (where the non-zero summation of the interlayer attraction forces may have resulted in part from the interlayer friction,36, 40 and in part from the slight difference in temperatures between the regions of the inner tube contacting with the outer tube ends49). Besides the forces mentioned above, there is another force called the thermally induced edge force (Fe)40, 41, 50 acting on each end of the inner tube in Model I (see Regions I and VI in Figures 3a and 3e) where the inner tube end is in contact with the outer tube wall. This edge force is induced by the differential atomic vibrations of the outer tube atoms in contact with and those not in contact with the inner tube (see Figure 5), and is linearly dependent on the

ACS Paragon Plus Environment

Page 13 of 24

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

temperature of the outer tube.40, 50, 51 Here we take five rings of atoms on the inner tube to calculate the edge force. It is found that FI is smaller than FVI, and the unbalanced edge force,

FI + FVI, is (7.191-11.976) = -4.79 pN, providing a driving force comparable to FIII + FIV (= -5.98 pN). It is thus clear that the driving mechanism in the two models are totally different. In Model I, the net driving force is a result of competitions among the thermal gradient force, the frictional force, the interlayer attraction force and the thermally induced edge force. The interlayer attraction force and the thermally induced edge force provide the driving force, whereas the thermal gradient force acts as a resisting force similar to the frictional force. This is why the inner tube in Model I tends to travel toward the cooler region (i.e., experiences positive thermophoresis). In contrast, in Model II, the net driving force results from the thermal gradient force, the frictional force, and the interlayer attraction force. It is mainly the thermal gradient force that drives the inner tube to move towards the hotter region (i.e., undergoes negative thermophoresis), whereas the interlayer attraction forces almost cancel each other out.

Discussion Thermal gradient. Both numerical40 and analytical41 studies have clearly shown that the thermal gradient force increases linearly with the temperature gradient, with a coeficient of 5.33 µeV/atom/K. However, in the present system, the gradient force is relatively small. The reason is that the real temperature gradient in the inner tube is very small (about 0.3 K/nm, see Figure 4) compared with the apparent temperature gradient between the two outer tubes

ACS Paragon Plus Environment

Nano Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

((500-300)/35 ≈ 8.6 K/nm). This is due to the fact that the intralayer thermal conductivity of a CNT is orders higher than the interlayer thermal conductivety,52-56 which leads to a very small temperature gradient along the inner tube but a quite high temperature difference between the outer and inner tubes. Temperature dependent interlayer attraction force and edge force. Owing to an entropy difference,49 the interlayer attraction force on the inner tube decreases with the temperature. In spite of the existence of the frictional force, we can see that the interlayer attraction force from the cooler outer tube is higher than that from the hotter one in both models (see Figure 6a, open squares). In trying to eliminate the effect of friction on the calculation of the interlayer attraction force, we imposed a constant small velocity of 2.5 nm/ns at the center of the inner tube in the z direction (yet maintain the temperatures of two outer tubes at 300 and 500 K, respectively) to obtain an approximately quasi-static mean value of the interlayer attraction force. The difference between the results (Figure 6a, open squares) obtained for the thermophoretic motion and the corresponding ones (Figure 6a , solid squares) for the imposed motion should be the friction force (which is estimated about 0.5 pN from each end of the outer tubes, on the same order of the previously calculated result40). The calculated interlayer attraction force in four regions exhibits a clear decreasing trend with temperature. To more clearly show how the temperatures of both the inner tube and outer tube influence the attraction force, by imposing a constant velocity of 2.5 nm/ns on the inner tube of the Model I system, and separately varying the temperatures of the inner tube and two outer tubes, we obtain an approximately linear temperature dependence of the attraction force on the

ACS Paragon Plus Environment

Page 14 of 24

Page 15 of 24

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

temperatures of both the inner and outer tubes (see Figure 6b).

Figure 6. (a) Amplitude of the interlayer attraction forces (open squares) in Regions I, III, IV and VI in Figure 4(f) and the corresponding forces (solid squares) under uniform motion of the inner tube. (b) Dependence of the interlayer attraction force on the temperatures of the inner and outer tubes. (c) Dependence of the edge force on the temperature. (d) The random force can help the inner tube to escape out of the trap (the overall motion behavior of the inner tube can be seen in Figure 2d).

Similar procedure to the calculation of the interlayer attraction force can be used to obtain the edge force acting on the inner tube in Model I. The results show that the edge force is approximately linearly dependent on the outer tube temperature and less sensitive to the inner tube temperature (shown in Figure 6c). This means that the linear temperature dependence of the edge force on a caped CNT is the same as that of the edge force on an open-ended CNT40

ACS Paragon Plus Environment

Nano Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

or on a graphene layer50 from an adjacent layer. Role of random force. As has been mentioned, the long term average of the random forces vanishes and thus it has no net contribution to the driving force. However, at the initial motion of the inner tube at rest, the random force may help it get out the possible trap of potential energy corrugation,57-59 if the driving force is small compared with the potential barrier. As shown in Figure 6d, when its velocity decelerates to zero, the inner tube in Model II may fall into a potential trap (possibly due to its larger length which may lead to a larger structural distortion (see Figure 2b) and thus a larger interlayer resistance), and oscillate in a unit lattice with a vibration amplitude of about half of the lattice length la.60 With the aid of the random force, the inner tube can escape from the trap and start to move forward at some point. Geometry dependence of driving force. To give a clearer picture of the geometry dependent thermophoresis, we investigate the effect of a change in distance between the two outer tubes on the driving force of the device under consideration. The inner and outer tubes used in the simulations are the same as those in Model II, with distance between the outer tubes varying from 10 to 80 nm, and the temperatures of the two outer tubes still set to 300 and 500 K, respectively. The two ends of the inner tube is exposed outside of the outer tubes when the distance is smaller than 60 nm, i.e., the configuration of the device is similar to Model II, while when the distance is larger than 60 nm, the two ends of the inner tube is covered inside of the outer tubes, i.e., the configuration of the device is similar to Model I. It is seen from Figure 7 that the driving force first decreases with the distance between the outer tubes, toward the hotter side, in Model II like configurations. The decrease of the driving force is due to the decrease of the thermal gradient in the inner tube caused by the

ACS Paragon Plus Environment

Page 16 of 24

Page 17 of 24

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

increase of the distance between the outer tubes. However, when the two ends of the inner tube are covered inside of the outer tubes to form a Model I like configuration, the driving force reverses its direction, i.e., toward the cooler side, and maintains almost a constant magnitude. This is because the driving force mainly comes from the interlayer attraction force and the edge force, both of which are strongely dependent on the outer tube temperatures (which are independent of the distance between outer tubes) but less dependent on the inner tube temperature (which may change with the distance between outer tubes). The results in Figure 7 clearly show that the direction of the driving force on the inner tube can be reversed by a change in the distance betweent the outer tubes, which means that such a device can generate repeated mechanical operations61 through the geometry dependent thermophoresis of the inner tube by varying the relative axial position of the outer tubes.

Figure 7. Dependence of the driving force on the distance between the outer tubes. Exellent agreement between theoretical prediction and MD calculations confirms the reported driving mechanisms.

Theoretical model for the driving force. From the above discussion, we know that the

ACS Paragon Plus Environment

Nano Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 24

driving force in Model I device can be expressed as F I = ∑ Fe + ∑ Fa + ∑ Fg ,

(1)

while in Model II, F II = ∑ Fa + ∑ Fg .

(2)

The gradient force on the inner tube, i.e., the reaction force on the outer tube caused by the gradient on the inner tube, can be calculated as Fg = ηin ∇TπdmLcont/A0 (see Supporting Information), where ηin is a potential coefficient, Lcont the contact length between tubes, dm the mean diameter of the inner and outer tubes, ∇T the temperature gradient along the tube, and

A0 the area per atom. The edge force on an inner tube end is calculated as Fe = (ηedge−ηin) πdmTouter /A0 (see Supporting Information), where Touter is the outer tube temperature. The attraction force is linearly dependent on both of the inner and outer tube temperatures, i.e., Fa = F0 − αinTinner – αoutTouter, where F0 is the attraction force at zero temperature, and αin and αout are numerically determined coefficients for the device composed of (5, 5) and (10, 10) tubes. By establishing an equivalent thermal circuit of the device under consideration, we can determine the temperature gradient along the inner tube and thus determine all above force components (details are given in Supporting Information), and the initial driving force for the device considered in Figure 7 can be obtained as F II =

71.14L2 pN, for δ < l−2L = 60 nm, (189nm + L)( L + δ )

F I = −15.53 −

(l − δ ) ( 29.36δ − 41.78(l − δ ) ) pN, for δ > l−2L = 60 nm. (378nm + l − δ )(l + δ )

(3)

(4)

It is seen from Figure 7 that the theoretical predictions from Equations. (3) and (4) agree well with the MD simulations, further confirming the above revealed driving mechanism.

ACS Paragon Plus Environment

Page 19 of 24

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

Conclusion By investigating the thermophoretic motion of a single walled carbon nanotube nested inside of two separate outer tubes held at different temperatures, we have shown that the direction of motion of the inner tube is dependent on the system geometry. When both ends of the inner tube extended out of the outer tubes, the inner tube undergoes negative thermophoresis (i.e., moves towards higher temperature), whereas positive thermophoresis is favorable when the ends of the inner tube stay inside of the outer tubes. To the best of our knowledge, this is the first time that negative thermophoresis is observed at solid-solid interfaces. The mechanisms for the positive and negative thermophoretic motions in the devices under consideration are quite different. In positive thermophoresis, the interlayer attraction force and a thermally induced edge force serve as the driving force, whereas the gradient force produces a resisting force. By contrast, in negative thermophoresis, the driving force comes mainly from the thermal gradient force. A proposed theoretical model further comfirms the simulation revealed mechanisms. Since its discovery in liquids more than a century ago, thermophoresis has long been studied for a variety of applications ranging from mass transport to molecule manipulation. The newly discovered example of negative thermophoresis in solids fills an important knowledge gap and is expected to stimulate further interests on this intriguing phenomenon along with potential applications in nanoscale actuation and energy conversion systems.

ACS Paragon Plus Environment

Nano Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Associated Content Supporting Information The Supporting Information is available free of charge on the ACS Publications website at DOI: [A theoretical model for predicting the positive and negative thermophoretic forces.]

Author Information Corresponding author *E-mail: [email protected]

Author Contributions T.C. supervised the project and initiated the study. J.L. performed simulations. Z.G. and H.Z. provided assistance in the simulations. J.L., T.C. and X.G analyzed the data. J.L, T.C., and H.G. wrote the manuscript. All authors contributed to the discussion of the results.

Notes The authors declare no competing financial interest.

Acknowledgements The authors acknowledge the financial supports from the NSF of China (Grant No. 11425209) and Shanghai Pujiang Program (Grant No. 13PJD016). The MD simulations were carried out on the computing platform of the International Center for Applied Mechanics in Energy

ACS Paragon Plus Environment

Page 20 of 24

Page 21 of 24

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

Engineering (ICAMEE), Shanghai University.

References (1) Ludwig, C. Lehrbuch der Physiologie des Menschen. Sitzungsbericht. Kaiser. Akad. Wiss. 1856. (2) Duhr, S.; Braun, D. Proc. Natl. Acad. Sci. USA. 2006, 103, 19678-19682. (3) Rosner, D.; Mackowski, D.; Garcia-Ybarra, P. Combust. Sci. Technol. 1991, 80, 87-101. (4) Blums, E.; Odenbach, S.; Mezulis, A.; Maiorov, M. Phys. Fluids 1998, 10, 2155-2163. (5) Wiegand, S. J. Phys.: Condens. Matter 2004, 16, R357-R379. (6) Maeda, Y. T.; Tlusty, T.; Libchaber, A. Proc. Natl. Acad. Sci. USA. 2012, 109, 17972-17977. (7) Brown, M. B.; Martin, G. P.; Jones, S. A.; Akomeah, F. K. Drug. Deliv. 2006, 13, 175-187. (8) Maggi, C.; Saglimbeni, F.; Dipalo, M.; De Angelis, F.; Di Leonardo, R. Nat. Commun. 2015, 6, 7855. (9) Barreiro, A.; Rurali, R.; Hernandez, E. R.; Moser, J.; Pichler, T.; Forro, L.; Bachtold, A. Science 2008, 320, 775-778. (10) Thamdrup, L. H.; Larsen, N. B.; Kristensen, A. Nano Lett. 2010, 10, 826-832. (11) Vélez-Cordero, J. R.; Hernández-Cordero, J. Langmuir 2015, 31, 10066-10075. (12) Mast, C. B.; Schink, S.; Gerland, U.; Braun, D. Proc. Natl. Acad. Sci. USA. 2013, 110, 8030-8035. (13) Pedersen, J. N.; Lüscher, C. J.; Marie, R.; Thamdrup, L. H.; Kristensen, A.; Flyvbjerg, H. Phys. Rev. Lett. 2014, 113, 268301. (14) Reichl, M.; Herzog, M.; Götz, A.; Braun, D. Phys. Rev. Lett. 2014, 112, 198101. (15) Jerabek-Willemsen, M.; André, T.; Wanner, R.; Roth, H. M.; Duhr, S.; Baaske, P.; Breitsprecher, D. J. Mol. Struct. 2014, 1077, 101-113. (16) Soret, C. Sur L'etat D'équilibre, Que Prend au Point de Vue de sa Concentration une Dissolution Saline Primitivement Homogène Dont Deux Parties Sont Portées a des Températures Difféntes. Arch. Sci. Phys. Nat. 1879. (17) Tyndall, J. Nature 1870, 1, 339-342. (18) Niether, D.; Afanasenkau, D.; Dhont, J. K. G.; Wiegand, S. Proc. Natl. Acad. Sci. USA. 2016, 113,

ACS Paragon Plus Environment

Nano Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

4272-4277. (19) Schoen, P. A.; Walther, J. H.; Arcidiacono, S.; Poulikakos, D.; Koumoutsakos, P. Nano Lett. 2006, 6, 1910-1917. (20) Schoen, P. A. E.; Walther, J. H.; Poulikakos, D.; Koumoutsakos, P. Appl. Phys. Lett. 2007, 90, 253116. (21) Rurali, R.; Hernández, E. R. Chem. Phys. Lett. 2010, 497, 62-65. (22) Coluci, V. R.; Timóteo, V. S.; Galvão, D. S. Appl. Phys. Lett. 2009, 95, 253103. (23) Becton, M.; Wang, X. J. Chem. Theory Comput. 2014, 10, 722-730. (24) Zambrano, H. A.; Walther, J. H.; Koumoutsakos, P.; Sbalzarini, I. F. Nano Lett. 2009, 9, 66-71. (25) Longhurst, M.; Quirke, N. Nano Lett. 2007, 7, 3324-3328. (26) Bresme, F.; Lervik, A.; Bedeaux, D.; Kjelstrup, S. Phys. Rev. Lett. 2008, 101, 020602. (27) Dwyer, H. A. Phys. Fluids 1967, 10, 976-984. (28) Sone, Y. J. Phys. Soc. of Japan 1972, 33, 232-236. (29) Sone, Y.; Aoki, K. Progr. Astronaut. Aeronaut. 1977, 51, 417-433. (30) Loyalka, S. J. Aerosol Sci. 1992, 23, 291-300. (31) Takata, S.; Aoki, K.; Sone, Y. Rarefied Gas Dynamics- Theory and Simulations 1994, 626-639. (32) Gotsmann, B.; Dürig, U. Appl. Phys. Lett. 2005, 87, 194102. (33) Bosworth, R. W. An Experimental Study of Positive and Negative Thermophoresis. (University of Colorado Colorado Springs. Kraemer Family Library) 2016. (34) Plimpton, S. J. Comput. Phys. 1995, 117, 1-19. (35) Brenner, D.; Shenderova, O.; Harrison, J. J. Phys.: Condens. Matter 2002, 14, 7. (36) Tangney, P.; Louie, S. G.; Cohen, M. L. Phys. Rev. Lett. 2004, 93, 065503. (37) Tangney, P.; Cohen, M. L.; Louie, S. G. Phys. Rev. Lett. 2006, 97, 195901. (38) Schoen, P. A.; Walther, J. H.; Arcidiacono, S.; Poulikakos, D.; Koumoutsakos, P. Nano Lett. 2006, 6, 1910-1917. (39) Zambrano, H. A.; Walther, J. H.; Koumoutsakos, P.; Sbalzarini, I. F. Nano Lett. 2009, 9, 66-71. (40) Guo, Z.; Chang, T.; Guo, X.; Gao, H. Phys. Rev. Lett. 2011, 107, 105502. (41) Guo, Z.; Chang, T.; Guo, X.; Gao, H. J. Mech. Phys. Solids 2012, 60, 1676-1687. (42) Shenai, P. M.; Xu, Z.; Zhao, Y. Nanotechnology 2011, 22, 485702. (43) Santamaría-Holek, I.; Reguera, D.; Rubi, J. M. J. Phys. Chem. C 2013, 117, 3109-3113.

ACS Paragon Plus Environment

Page 22 of 24

Page 23 of 24

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

(44) Prasad, M. V.; Bhattacharya, B. Nano Lett. 2016, 16, 2174-2180. (45) Yu, M.-F.; Yakobson, B. I.; Ruoff, R. S. J. Phys. Chem. B 2000, 104, 8764-8767. (46) Saito, R.; Matsuo, R.; Kimura, T.; Dresselhaus, G.; Dresselhaus, M. Chem. Phys. Lett. 2001, 348, 187-193. (47) Zheng, Q.; Jiang, Q. Phys. Rev. Lett. 2002, 88, 045503. (48) Cumings, J.; Zettl, A. Science 2000, 289, 602-604. (49) Guo, Z.; Chang, T.; Guo, X.; Gao, H. arXiv preprint arXiv:1511.08980 2015. (50) Li, J.; Zhang, H.; Guo, Z.; Chang, T.; Gao, H. J.Appl. Mech. 2015, 82, 101011. (51) Chang, T.; Guo, Z. Nano Lett. 2010, 10, 3490-3493. (52) Che, J.; Cagin, T.; Goddard III, W. A. Nanotechnology 2000, 11, 65. (53) Pop, E.; Varshney, V.; Roy, A. K. MRS Bulletin 2012, 37, 1273-1281. (54) Fujii, M.; Zhang, X.; Xie, H.; Ago, H.; Takahashi, K.; Ikuta, T.; Abe, H.; Shimizu, T. Phys. Rev. Lett. 2005, 95, 065502. (55) Pop, E.; Mann, D.; Wang, Q.; Goodson, K.; Dai, H. Nano Lett. 2006, 6, 96-100. (56) Xu, Z.; Buehler, M. J. Nanotechnology 2009, 20, 185701. (57) Astumian, R. D.; Hänggi, P. Phys. Today 2002, 55, 33. (58) Tu, Z.; Hu, X. Phys. Rev. B 2005, 72, 033404. (59) Ueno, Y.; Somada, H.; Hirahara, K.; Nakayama, Y.; Akita, S. Jpn. J. Appl. Phys. 2009, 48, 06FG03. (60) Bonelli, F.; Manini, N.; Cadelano, E.; Colombo, L. Eur. Phys. J. B 2009, 70, 449-459. (61) Abendroth, J. M.; Bushuyev, O. S.; Weiss, P. S.; Barrett, C. J. ACS nano 2015, 9, 7746-7768.

ACS Paragon Plus Environment

Nano Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

TOC Graphic 252x85mm (120 x 120 DPI)

ACS Paragon Plus Environment

Page 24 of 24