Recent Advances in Nanocarbon MaterialsView this Special Issue
Thermal Resistance across Interfaces Comprising Dimensionally Mismatched Carbon Nanotube-Graphene Junctions in 3D Carbon Nanomaterials
In the present study, reverse nonequilibrium molecular dynamics is employed to study thermal resistance across interfaces comprising dimensionally mismatched junctions of single layer graphene floors with (6,6) single-walled carbon nanotube (SWCNT) pillars in 3D carbon nanomaterials. Results obtained from unit cell analysis indicate the presence of notable interfacial thermal resistance in the out-of-plane direction (along the longitudinal axis of the SWCNTs) but negligible resistance in the in-plane direction along the graphene floor. The interfacial thermal resistance in the out-of-plane direction is understood to be due to the change in dimensionality as well as phonon spectra mismatch as the phonons propagate from SWCNTs to the graphene sheet and then back again to the SWCNTs. The thermal conductivity of the unit cells was observed to increase nearly linearly with an increase in cell size, that is, pillar height as well as interpillar distance, and approaches a plateau as the pillar height and the interpillar distance approach the critical lengths for ballistic thermal transport in SWCNT and single layer graphene. The results indicate that the thermal transport characteristics of these SWCNT-graphene hybrid structures can be tuned by controlling the SWCNT-graphene junction characteristics as well as the unit cell dimensions.
Based on the exceptional thermal transport characteristics as well as the electronic mobility of one- and two-dimensional carbon nanomaterials (e.g., carbon nanotube and graphene), there have been several attempts by researchers to synthesize three-dimensional (3D) hybrid nanostructures that synergistically combine two or more low-dimensional carbon nanomaterials [1–7]. Amongst the several potential 3D carbon nanostructured materials, 3D pillared single-walled carbon nanotube- (SWCNT-) graphene architectures (PGS) that incorporate SWCNTs as pillars and single layer graphene as floors have attracted particular attention [1, 4, 6, 8, 9]. These materials have the potential to exploit the extraordinary thermal conductivity of their base nanomaterials to enable tunable thermal transport in both out-of-plane and in-plane directions. Varshney et al. , using molecular dynamics, were the first to report thermal transport properties of PGS. More recently, Loh et al.  have shown that long wavelength out-of-plane phonon modes contribute significantly to thermal transport in such structures.
In the present paper we focus attention on thermal resistance across dimensionally mismatched SWCNT-graphene junctions in the 3D PGS. Using classical molecular dynamics with adaptive intermolecular reactive empirical bond order (AIREBO) interatomic potential interfacial thermal resistances are computed and the thermal conductivities of 3D PGS are compared with those of pristine SWCNT and single layer graphene of comparable size. The presence of the interfacial thermal resistance at the junctions is hypothesized to be due to the change in dimensionality (from nearly 1D to 2D to 1D) as well as phonon spectra along the thermal transport path as the phonons propagate from the SWCNT to the graphene sheet and then back again to the SWCNT. In order to investigate this we compute the local density of states at various locations along the phonon transport path from SWCNT to graphene to SWCNT.
2. Simulation Method
In the present study, SWCNT-graphene junction structures were constructed by translating carbon atoms of (6,6) armchair SWCNTs to match with corresponding holes in free standing single layer graphene nanostructures that were built using the visual molecular dynamics (VMD)  software. The procedure resulted in the creation of SWCNT-graphene junctions with regularly distributed septagon/hexagon carbon rings [12, 13] with sp2 type hybridization at the SWCNT-graphene junctions. Figure 1(a) illustrates the creation of such junction. It is well known that, in graphene, the three outer-shell electrons of each carbon atom occupy the planar sp2 hybrid orbital to form three in-plane σ bonds with an out-of-plane π orbital (bond). This makes graphene a planar hexagonal network. The length and energy of the σ bonds in the sp2 hybridization are 0.14 nm and 420 kcal/mol, respectively, and 0.15 nm and 360 kcal/mol in the sp3 configuration. The planar hexagonal structure along with the strong covalent (sp2) bonding is understood to be responsible for the high in-plane phonon group velocities and thus high in-plane thermal conductivity of graphene . When a graphene sheet is rolled to form a nanotube, the circular curvature of the nanotubes promotes quantum confinement and σ-π rehybridization (i.e., σ-π bond mixing); the three planar σ bonds are forced slightly out of plane and for compensation the π orbital is more delocalized outside the tube. This rehybridization together with the π electrons confinement gives nanotubes their unique, extraordinary mechanical, thermal, and electronic properties. Moreover, they allow for topological defects, such as pentagons and septagons, to be incorporated into the hexagonal networks to form bent, toroidal, and helical nanotubes. For the purposes of the present study, the formation of the SWCNT-graphene junctions with sp2 hybridized covalent bonds between SWCNT and graphene involves creating six septagonal and six distorted hexagonal carbon rings, as shown in Figure 1(a); the existence of these seven member rings and distorted hexagon rings in a SWCNT-graphene junction has been confirmed by experiments .
Figures 1(b) and 1(c) illustrate the detailed configuration of a hexagonal carbon ring in a (6,6) SWCNT and a septagon/hexagon carbon ring in a SWCNT-graphene junction obtained from the molecular dynamics simulations. The average bond length in a hexagonal carbon ring of a (6,6) SWCNT is 0.1398 nm while the average bond length in a typical distorted hexagonal carbon ring from the SWCNT-graphene junction region is 0.1425 nm. Moreover, the average bond length in a typical septagonal carbon ring from the SWCNT-graphene junction region is found to be 0.1434 nm. Charlier et al.  using tight-binding models to investigate topological defects in SWCNTs estimated the bond length in a hexagon ring of a single 5/6/7 defect to be 0.1421 nm and the bond length in a septagon ring of a single 5/7 defect to be 0.1433 nm. This confirms that the distorted hexagonal and septagonal carbon rings in SWCNT-graphene junctions have similar molecular structures to the topological defects observed in pristine SWCNTs and are expected to play the role of defects in SWCNTs [16, 17] giving rise to phonon scattering and thus an interfacial thermal resistance at the SWCNT-graphene junctions.
In the present research reverse nonequilibrium molecular dynamics (RNEMD) [16, 18–20] is used to estimate the thermal conductivities in the hybrid structures as well as the interfacial thermal resistance at the SWCNT-graphene junctions. Two types of simulation structures (unit cells), depicted in Figure 2, are constructed according to the direction in which thermal conductivity is to be evaluated. Temperature profiles are obtained along the SWCNT axis, that is, along the direction, as shown in Figure 2(a), to investigate the effect of SWCNT-graphene junction on the thermal transport characteristics in the cross-plane direction. We will denote this as “out-of-plane” thermal conductivity. On the other hand, the “in-plane” thermal conductivity is defined as the thermal conductivity along the graphene layer, that is, along the direction in Figure 2(b), which is useful for understanding the effects of SWCNT-graphene junction on thermal transport (phonon scattering) along the graphene layers. Two control parameters used in the construction of SWCNT-graphene junction structures are the pillar height, , and the interpillar distance, . Pillar height is the length of the SWCNT including the junction region that connects the two adjacent graphene floors. It also determines the sample length for thermal conductivity calculations in the out-of-plane direction (Figure 2(a)). Three different pillar heights, that is, = 50 nm, 100 nm, and 200 nm, are chosen for the study. Interpillar distance, , is defined as the distance between two SWCNT-graphene junctions and determines the sample length for the in-plane thermal conductivity calculations—the sample size is three times the interpillar distance (Figure 2(b)). For the four interpillar distances ( = 5 nm, 10 nm, 15 nm, and 20 nm) chosen in the present study, the in-plane thermal conductivity is calculated for four different sample sizes, that is, four different lengths of graphene layers, which are 15 nm, 30 nm, 45 nm, and 60 nm.
All simulations are performed using the large-scale atomic/molecular massively parallel simulator (LAMMPS)  code using the AIREBO  interatomic potential. The simulation structure is energy-minimized to reach the equilibrium configuration by iteratively adjusting atomic coordinates. After elevating the temperature (500 K) using isenthalpic (NPH) ensemble with Langevin thermostat for faster equilibration, the system is cooled to room temperature with the isothermal-isobaric (NPT) ensemble. RNEMD is carried out on the equilibrated structures for 2 ns with the microcanonical (NVE) ensemble with a time step of 1 fs. Figure 3 shows the total (cumulative) energy transported and the evolution of the temperature profile as a function of time for the SWCNT-graphene junction structure with 50 nm pillar height and 5 nm interpillar distance. As seen from the figure the system quickly approaches a steady state—the quick stabilization of heat flow and the temperature distribution is due to the well-equilibrated molecular structure. In addition, during RNEMD with NVE ensemble the total energy is found to be conserved adequately with only small energy fluctuations of ~10 eV. These observations are true for all SWCNT-graphene junction structures utilized in this work. Also, in all simulated structures thermal conductivities and interfacial thermal resistances are obtained from the temperature profiles after 2 ns of RNEMD simulations to ensure achievement of steady state transport conditions.
Using Fourier’s law, the thermal conductivity, k, is obtained as where is the heat flux and is temperature gradient averaged over time and space. From the temperature jumps in the temperature profiles, interfacial thermal resistance, , can be calculated by using In the calculation of the thermal conductivity, only areas that directly contribute to thermal transport are used—annular area of one (6,6) SWNCT (0.335 nm thickness) is used in the out-of-plane thermal conductivity calculation while sectional areas of two graphene layers (each with 0.335 nm thickness) are used for the in-plane thermal conductivity calculation, as shown in Figure 4.
3. Results and Discussion
Table 1 lists the out-of-plane and in-plane thermal conductivities and the interfacial thermal resistances obtained from all the simulations conducted in the present study. The interpillar distance was varied from 5 nm to 20 nm while the pillar height was changed from 25 to 200 nm. To save computational time the in-plane thermal conductivities for unit cells with 20 nm interpillar distance and pillar heights of 50 nm and 100 nm were omitted. The validation of the simulation methods was provided in our previous work  in which the theoretical predictions of the thermal conductivity of pristine (6,6) SWCNT and single layer graphene were compared with experimental measurements reported by other researchers.
Figure 5(a) shows the temperature profile for out-of-plane thermal conductivity in a unit cell with a 50 nm pillar height obtained after imposing a heat flux for a total duration of 2 ns. The temperature gradient is determined from the slope of the linear region (shown as a solid black line in Figure 5(a)) excluding the hot and cold bath regions and also excluding the region of the sudden temperature drop. The temperature slips at the specimen boundary are a consequence of quasi-ballistic thermal transport due to the relatively short length of the sample when compared to the phonon mean free path. Fourier theory of thermal transport considers heat transport as a diffusive process where energy flow is driven by a temperature gradient. However, this is not valid at length scales smaller than the mean free path for the energy carriers in a material, which can be tens of nanometers in carbon-based nanomaterials at room temperature. In this case, the heat flow can become ballistic, as seen by the temperature slips near the ends (hot and cold baths) of the simulation box. Distinct temperature jumps can also be observed at the SWCNT-graphene junction locations in the temperature profile. The distinct temperature jumps indicate the existence of an interfacial thermal resistance as heat flows from one SWCNT to the other via the two SWCNT-graphene junctions and the adjoining graphene region.
Based on the observed temperature drops in the temperature profiles, interfacial thermal resistance is calculated for the nine different SWCNT-graphene junction structures and the results are shown in Table 1. The interfacial thermal resistance for out-of-plane thermal transport is observed to decrease slightly with an increase in pillar height. Also, it is observed that the out-of-plane thermal conductivity depends weakly on the interpillar distance. For the case of the in-plane thermal conductivity the temperature gradient is determined from the slope of the linear region (shown as a solid black line in Figure 5(d)) excluding the hot and cold bath regions. It is interesting to note that unlike the case with out-of-plane thermal transport the sudden temperature jumps are not evident at the locations of the SWCNT-graphene junctions in the in-plane temperature profiles. An important implication of this observation is that phonon scattering during in-plane thermal transport is weaker (or nearly absent) when compared to the out-of-plane thermal transport in these SWCNT-graphene hybrid structures. From the linear regression of the slope, equivalent in-plane thermal conductivities for 10 different SWCNT-graphene junction structures are calculated, and the results are shown in Table 1. The in-plane thermal conductivity shows a small dependency on pillar height and the effect becomes smaller as the interpillar distance increases.
The slight decrease in the interfacial thermal resistance with an increase in SWCNT pillar height during out-of-plane thermal transport can be attributed to scattering of low frequency (i.e., large wavelength) phonons [23, 24]. On the other hand, the reasons for the slight increase in in-plane thermal conductivity with an increase in SWCNT pillar height when the interpillar distance is small (5 nm) are much more complex. Further studies, including the wave packet analysis, will be undertaken in the future to better understand these effects. However, the small dependency of in-plane thermal conductivity on pillar height when the interpillar distance is large can be attributed to the near diffusive thermal transport in graphene in the region between the junctions. Along the same vein, Bagri et al.  studied the effect of twin grain boundary on the thermal transport of graphene and showed a similar trend in thermal conductivity with different graphene lengths in between the defected regions. They showed that the thermal conductivities of graphene with 13.2-degree grain boundary orientation were higher than the ones with 5.5-degree grain boundary orientation for structures with relatively short graphene lengths between defects (25 nm and 50 nm) but could not observe a clear difference between the thermal conductivities for structures with 13.2-degree grain boundary orientation and the 5.5-degree grain boundary orientation when the distance between defects was relatively larger (125 nm).
With regard to the distinct temperature jumps in the temperature profiles (Figure 5) during out-of-plane thermal transport, we hypothesize that the jumps correspond to phonon dispersion due to the sudden change in dimensionality which the phonons experience as they flow from the quasi-1D structure (SWCNT) to 2D (graphene) and then back again to quasi-1D (SWCNT) via the SWCNT-graphene junctions. As a consequence, the low frequency longitudinal acoustic phonon modes, which carry thermal energy efficiently in a SWCNT, must transform into phonon modes which are efficient energy carriers of thermal energy in graphene, for example, the in-plane modes, and then back to the preferred phonon modes in SWCNTs. In order to better understand the contributions of SWCNT-graphene junctions we investigate the local phonon density of states (LDOS) using velocity autocorrelation along the phonon transport path. It is expected that the subtle differences in the LDOS in various regions along the path of the phonon transport, that is, the SWCNT-graphene junctions, graphene floor, and the adjoining SWCNTs, would help to provide insights into the origin of thermal interfacial resistance at the SWCNT-graphene interface in these systems. For example, a mismatch in the LDOS in the SWCNT and the graphene regions would represent a mismatch in the phonon spectra between the two regions and hence to a higher thermal resistance across the interface.
Figure 6(a) shows the LDOS in the SWCNT and the graphene floor from regions away from the SWCNT-graphene junction. It can be observed that the LDOS distribution is similar except in the low energy (frequency) regime between 60 meV and 110 meV. Since these low energy (frequency) phonon modes are expected to contribute more to the thermal transport characteristics when compared to the high energy (frequency) phonon modes, the change in dimensionality from 1D to 2D to 1D is expected to contribute to the observed thermal interfacial resistance during out-of-plane thermal transport in these structures. On the other hand, the LDOS for in-plane thermal transport in regions close to and away from the SWCNT-graphene junctions is shown in Figure 6(b). It can be observed that the LDOS in the two regions is very similar except for the higher energy (frequency) phonon modes with energy greater than 190 meV. Since most of the higher energy phonons reside in the higher frequency optical modes which have relatively low group velocities they are expected to contribute little to the overall thermal conductivity and hence to a negligible interfacial thermal resistance.
Figure 7 shows a comparison of the thermal conductivity of the 3D SWCNT-graphene hybrid structure with the thermal conductivity of pristine (6,6) SWCNT and single layer graphene of corresponding size. In Figure 7(a) thermal conductivities of individual (6,6) SWCNT are compared with out-of-plane thermal conductivities of SWCNT-graphene hybrid structures. The out-of-plane thermal conductivity of the 3D hybrid structures is observed to increase with increasing size of the unit cell, that is, with an increase in pillar height (length of the SWCNT). It is interesting to note that the increase in out-of-plane thermal conductivity in SWCNT-graphene hybrid structure is similar to that observed in pristine (6,6) SWCNT. In our previous study of thermal transport in (6.6) SWCNTs  we observed that thermal conductivity of SWCNTs is linearly dependent on specimen length as long as the length of the SWCNTs is smaller than the mean free path (ballistic thermal transport) and then transitions to a regime where it becomes nearly length-independent (diffusive and governed by Fourier’s law). Consequently, the out-of-plane thermal conductivity of the 3D hybrid structures is expected to occur nearly linearly with pillar height as long as the height is smaller than the critical length for ballistic thermal transport in (6,6) SWCNT and then approaches a plateau as the pillar height becomes greater than the ballistic length ~75 nm. In addition, due to the presence of interfacial thermal resistance at the SWCNT-graphene junction, the thermal conductivity of the hybrid structure is consistently lower than that of pristine (6,6) SWCNT. As seen in the temperature profiles describing the out-of-plane thermal transport (Figures 5(a), 5(b), and 5(c)), thermal vibrations must travel from SWCNT via SWCNT-graphene junctions to reach the next SWCNT; this leads to an approximately 5% decrease in thermal conductivity in the out-of-plane direction.
Figure 7(b) compares thermal conductivity of single layer graphene with in-plane thermal conductivity obtained from simulations using various sized SWCNT-graphene hybrid structures. Like in the case of out-of-plane thermal conductivity, the thermal conductivities of the hybrid structures are observed to increase linearly with an increase in cell size (interpillar distance). It is postulated that the increase in thermal conductivity with interpillar distance is primarily due to the size of the structures evaluated in the present paper, which are smaller than the critical ballistic length in single layer graphene for a given graphene floor width. Moreover, since the hybrid structures investigated herein are created such that the width of graphene layer () is the same as the interpillar distance (), it is expected that the ratio of defected (junction) area to total graphene area becomes smaller as the interpillar distance becomes larger. Use of larger sized structures (greater than the critical ballistic length in graphene) in the present analysis becomes computationally prohibitive. As the interpillar distance increases, it is expected that the difference between the in-plane thermal conductivities of the 3D hybrids and the single layer graphene (of corresponding size) converges.
In summary, in the present paper we present results of a reverse nonequilibrium molecular dynamics study with AIREBO interatomic potential to investigate thermal transport in 3D SWCNT-graphene hybrid structures comprising SWCNT-graphene junctions. The results of the simulations reveal the presence of interfacial thermal resistance for thermal transport in the out-of-plane direction (along the axis of the SWCNT pillars) while the thermal resistance along the graphene floor (in-plane direction) is found to be negligible. The origin of thermal resistance in the out-of-plane direction is understood to be a combination of phonon scattering due to the change in dimensionality of the medium along the phonon path (i.e., from quasi-1D in SWCNT to 2D in graphene) and a mismatch of phonon spectra along the phonon transport path. In addition, the interfacial thermal resistance is found to be weakly dependent on the pillar (SWCNT) height as well as the interpillar distance. For out-of-plane thermal transport the highest interfacial thermal resistance across the SWCNT-graphene junction was computed to be 1.93 × 10−10 K-m2/W for a pillar height of 50 nm and an interpillar distance of 15 nm. Moreover, the highest out-of-plane thermal conductivity of a unit cell representing the SWCNT-graphene hybrid structure was obtained to be W/m-K for a cell with pillar height of 200 nm and an interpillar distance of 5 nm. In comparison the thermal conductivity of a pristine (6,6) SWCNT with 200 nm length is computed to be W/m-K. The highest in-plane thermal conductivity of a unit cell with two SWCNT-graphene junctions and a cell size of 60 nm with 25 nm pillar height was estimated to be W/m-K; the thermal conductivity of pristine single layer graphene with the same dimensions was computed to be W/m-K. Moreover, the thermal conductivity of the hybrid unit cells was observed to increase nearly linearly with an increase in dimensions of the cell, that is, pillar height as well as interpillar distance. The thermal conductivity of the cell is observed to approach a plateau as the pillar height and the interpillar distance approach the critical lengths for ballistic thermal transport in SWCNT and single layer graphene of comparable sizes. These results indicate that the out-of-plane and in-plane thermal transport characteristics of these 3D SWCNT-graphene hybrid structures can be tuned by controlling the unit cell size.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
This work made use of the High Performance Computing Resource in the Core Facility for Advanced Research Computing at Case Western Reserve University. The authors also would like to acknowledge the support of the Air Force Office of Scientific Research (AFOSR) MURI Grant no. FA9550-12-1-0037 (Program Manager: Dr. Joycelyn Harrison) for conducting this research.
J.-C. Charlier, T. W. Ebbesen, and P. Lambin, “Structural and electronic properties of pentagon-heptagon pair defects in carbon nanotubes,” Physical Review B—Condensed Matter and Materials Physics, vol. 53, no. 16, pp. 11108–11113, 1996.View at: Google Scholar
F. Müller-Plathe, “A simple nonequilibrium molecular dynamics method for calculating the thermal conductivity,” Journal of Chemical Physics, vol. 106, no. 14, pp. 6082–6085, 1997.View at: Google Scholar
S. Plimpton, P. Crozier, and A. Thompson, LAMMPS-Large-Scale Atomic/Molecular Massively Parallel Simulator, Sandia National Laboratories, 2007.
S. J. Stuart, A. B. Tutein, and J. A. Harrison, “A reactive potential for hydrocarbons with intermolecular interactions,” Journal of Chemical Physics, vol. 112, no. 14, pp. 6472–6486, 2000.View at: Google Scholar