This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Evidence for Phonon-Assisted Intertube Electronic Transport in an Armchair Carbon Nanotube Film

Davoud Adinehloo Department of Electrical Engineering, University at Buffalo, Buffalo, NY 14228, USA    Weilu Gao Department of Electrical and Computer Engineering, University of Utah, Salt Lake City, UT 84112, USA    Ali Mojibpour Department of Electrical and Computer Engineering, Rice University, Houston, TX 77005, USA    Junichiro Kono Department of Electrical and Computer Engineering, Rice University, Houston, TX 77005, USA Department of Physics and Astronomy, Rice University, Houston, TX 77005, USA Department of Materials Science and NanoEngineering, Rice University, Houston, TX 77005, USA The Smalley-Curl Institute, Rice University, Houston, TX 77005, USA    Vasili Perebeinos Department of Electrical Engineering, University at Buffalo, Buffalo, NY 14228, USA [email protected]
Abstract

The electrical conductivity of a macroscopic assembly of nanomaterials is determined through a complex interplay of electronic transport within and between constituent nano-objects. Phonons play dual roles in this situation: their increased populations tend to reduce the conductivity via electron scattering while they can boost the conductivity by assisting electrons to propagate through the potential-energy landscape. Here, we identify a phonon-assisted coherent electron transport process between neighboring nanotubes in temperature-dependent conductivity measurements on a macroscopic film of armchair single-wall carbon nanotubes. Through atomistic modeling of electronic states and calculations of both electronic and phonon-assisted junction conductances, we conclude that phonon-assisted conductance is the dominant mechanism for the observed high-temperature transport. The unambiguous manifestation of coherent intertube dynamics proves a single-chirality armchair nanotube film to be a unique macroscopic solid-state ensemble of nano-objects promising for the development of room-temperature coherent electronic devices.

Toward large-scale applications of nanomaterials in electronics Jariwala et al. (2013); Okimoto et al. (2010); Rao et al. (2018); Gao et al. (2020), understanding and controlling electron transport processes, not only within each nano-object but also between them, is crucial Zorn and Zaumseil (2021). The overall electrical conductivity of a macroscopic assembly of nanomaterials is determined by an array of interdependent quantities – e.g., defect density, doping level, Fermi energy, material size, the density of nano-objects, purity, homogeneity, morphology, and temperature. The role of phonons in this highly complex situation is subtle since they tend to scatter electrons to reduce the conductivity while simultaneously assisting electrons to go through the potential-energy landscape to increase the conductivity. This dual role of phonons in electronic transport across macro-objects has not been elucidated. Particularly, phonon-assisted processes, including phonon-assisted coherent electron transfer between nano-objects, have not been identified in a macroscopic sample.

Here we study macroscopic assemblies of carbon nanotubes (CNTs), which provide an ideal model system in which to address the above issues and questions. Since their discovery in the early 1990s,Iijima (1991); Iijima and Ichihashi (1993); Bethune et al. (1993) extensive studies have revealed and established the truly unique electronic and optical properties of these one-dimensional nano-objects, particularly on individual nanotube levels. Depending on the precise atomic arrangements of sp2sp^{2}-bonded carbon atoms in the honeycomb lattice, specified by a pair of integers called chirality indices, (nn,mm), both metallic and semiconducting species of CNTs exist; details of the band structure, most importantly the band gap, are determined by (nn,mm).Dresselhaus et al. (2001); Avouris et al. (2007, 2008); Nanot et al. (2013); Weisman and Kono (2019) Intertube electronic transport has also been studied theoreticallyNakanishi and Ando (2001); Buldum and Lu (2001); Maarouf and Mele (2011) and experimentally.(Fuhrer et al., 2000a, b) Based on the momentum-matching conditions on the initial and final states of electrons in a transfer between two individual CNTs with a crossing angle θ\theta (see Figure 1a), the existence of special values of θ\theta, where the intertube conductance becomes maximum and minimum has been identified Nakanishi and Ando (2001); Buldum and Lu (2001); Maarouf and Mele (2011). We show that after structural relaxation in a CNT- crossing geometry, the overall junction conductance increases by an order of magnitude with a quantitatively similar crossing angle dependence. Moreover, following Reference 19 for twisted bilayer phonon-assisted junction conductance, we demonstrate here that the intertube phonon-assisted junction conductance at room temperature is comparable to the pure electronic conductance. At some crossing angles, phonon-assisted junction conductance exceeds its electronic counterpart by an order of magnitude. These two effects – geometry relaxation and the additional phonon-assisted mechanism – enable us to achieve better quantitative agreement with prior experimentsFuhrer et al. (2000a) as well as the data reported here.

In our previous work,(Gao et al., 2021) we investigated the temperature dependence of conductivity in a set of macroscopic single-chirality CNT films with well-defined (nn,mm). Distinctly (nn,mm)-dependent and strongly temperature-dependent conductivity was observed, and the overall behaviors were explained through the Mott variable range hopping (VRH) model in a wide temperature range. However, one of the samples – an armchair (6,6) CNT film – exhibited clear deviation from the VRH behavior, especially at elevated temperatures, indicating that a different transport mechanism was at work. Through detailed quantitative analysis of the temperature-dependent conductivity, we found that the localization length was longer in (6,6) CNT film than in the other chirality samples such that intertube transport, as opposed to intratube transport, dominated in the (6,6) film at high temperatures. Here, through modeling and calculations of both electronic and phonon-assisted junction conductances, we unambiguously identify phonon-assisted conduction as the dominant mechanism for the observed high-temperature transport behaviour. The clear manifestation of this coherent dynamic process makes a single-chirality armchair CNT film a unique macroscopic object in which to study quantum transport processes at room temperature.

The device we used in this study was based on a thin film of randomly oriented, chirality-enriched (6,6) CNTs. The device had four electrodes, and the channel between the two inner electrodes, shown in the bottom panel of Figure 1a, had a length (LL), width (WW), and thickness (tt) of 8 μ\mum, 5 μ\mum, and \sim12 nm, respectively. Further experimental details are fully described in Reference 20. The experimentally measured film conductance, GfilmG_{\rm film}, which is equal to the film conductivity times Wt/LWt/L, is plotted against temperature in Figure 1b, showing that GfilmG_{\rm film} monotonically but sublinearly increases with temperature as the temperature is raised from 4.2 K to 300 K.

Refer to caption
Figure 1: (a) Schematic illustration of the device used in this study. θ\theta - crossing angle, φ\varphi - rotating angle, hh - interlayer distance, and SS - sliding shift between CNTs. (b) Experimentally measured (red circles) and theoretically calculated (blue solid) film conductance of the (6,6) CNT film as a function of temperature. (c) Theoretically obtained total resistance of the film (blue solid) and the decomposed contributions: intratube (red dashed) and a sum of the two electronic and phonon-assisted contributions to the intertube junction resistance (green dotted dashed). The calculations were performed at Fermi energy EF=10E_{\text{F}}=10 meV and applied bias across the junction V=1V=1 mV.

To understand the observed temperature dependence of the film conductance, we consider a series resistance model of the CNT conductivity and junction resistance, which contribute to the overall resistance of the film asShim et al. (2015a):

Gfilm1=C(Lcntσintra+1Ginter),G_{\text{film}}^{-1}=C\left(\frac{L_{\text{cnt}}}{\sigma_{\text{intra}}}+\frac{1}{G_{\text{inter}}}\right), (1)

where the first and second terms originate from the intratube conductivity σintra\sigma_{\text{intra}} and the intertube junction conductance GinterG_{\text{inter}}, respectively. A geometrical prefactor, C=L(d+h)/(tW)C=L(d+h)/(tW), is determined by the average number of junctions in the film, where dd is the CNT diameter, hh is the interlayer distance between CNTs, and LcntL_{\text{cnt}} is the average length of CNTs. If we assume that the average crossing angle in a randomly oriented CNT film is θ=45\theta=45^{\circ}, CC should be 0.160.16, according to our geometry, while the best fit to the data was obtained for C=0.09C=0.09, as shown in Figure 1b.

A temperature-dependent film conductance can arise from various sources. Our previous study on single-chirality CNT filmsGao et al. (2021) has shown that VRH is the dominant mechanism at low temperatures. Therefore, we model σintra\sigma_{\text{intra}} by the 1D Mott VRH model(Gao et al., 2021), i.e., σintra=σ0exp[(T0/T)1/2]\sigma_{\text{intra}}=\sigma_{0}\exp[-(T_{0}/T)^{1/2}]. In Figure 1b, the theory curve according to Equation (1) includes the electronic and phonon-assistant contributions to the intertube (junction) conductance (discussed below) as well as the intratube conductivity. Since the sample is randomly oriented and the junction conductance depends on the angle θ\theta (see Figure 1a), we averaged the intertube junction conductance over θ\theta between 2020^{\circ} and 7070^{\circ}. The relative contributions of the intratube and intertube resistances are shown in Figure 1c. At temperatures below 50 K, the intratube resistance dominates, whereas, at temperatures above 50 K, the intertube junction resistance gives a dominant contribution to the film resistance.

Refer to caption
Figure 2: (a) Electronic conductance, Ginter, elG_{\text{inter, el}}, of crossed CNTs as a function of the crossing angle, θ\theta, between armchair CNTs of different diameters. (b) Surface plot of Ginter, elG_{\text{inter, el}} versus SS and φ\varphi in CNTs. Ginter, elG_{\text{inter, el}} versus φ\varphi in (c) and versus SS in (d) at horizontal and vertical cuts depicted in (b). In (b), (c), and (d), we used θ=70\theta=70^{\circ} and (6,6) CNTs. All results were calculated at room temperature, EF=0E_{\text{F}}=0 eV, and V=1V=1 mV.

Below we discuss calculations of the intertube junction conductance. In the case of an ideally transparent contact, the maximum intertube conductance equals 4e2/h1554e^{2}/h\approx 155 μ\muS. However, an experimentally measured conductanceFuhrer et al. (2000a) is substantially smaller; therefore, we use perturbation theory (see Methods) to calculate the intertube junction conductance as a sum of two parallel channels: electronic conductance Ginter,elG_{\rm inter,el} due to the overlap of π\pi orbitals on neighboring tubes described by the tight-binding model and phonon-assistant conductance Ginter,phG_{\rm inter,ph} (see Methods). Since the distance between π\pi-orbitals in neighboring CNTs determines the hopping overlap, intertube junction conductance varies with the CNT structure registry. Atomic positions in an unrelaxed structure depend on the crossing angle (θ\theta), rotation angle (φ\varphi), and sliding shift (SS) between adjacent CNTs; see Figure 1a. In Figure 2a we show Ginter, elG_{\text{inter, el}} for an unrelaxed structure versus crossing angle θ\theta in different armchair CNTs, including (6,6) CNTs; Ginter,elG_{\rm inter,el} can vary by more than two orders of magnitude as θ\theta changes. Figures 2b–2d show that Ginter,elG_{\rm inter,el} varies with φ\varphi and SS by about 30% in (6,6) CNTs. Figure 2b is a surface map of Ginter,elG_{\rm inter,el} versus φ\varphi and SS with θ=70\theta=70^{\circ}, while Figures 2c and 2d are detailed φ\varphi and SS dependences of Ginter,elG_{\rm inter,el}, corresponding to the labeled solid lines (ø1ø4\o 1-\o 4) in Figure 2b. The variation of Ginter,elG_{\rm inter,el} with φ\varphi and SS are periodic with periods of 6060^{\circ} and 2.46 Å, respectively, and hence, Ginter,elG_{\rm inter,el} in Figure 2a was obtained as an average over φ\varphi and SS.

Refer to caption
Figure 3: (a) Electronic (dashed blue) and phonon-assisted (solid red) intertube junction-conductance of relaxed (6,6) CNTs as a function of θ\theta at room temperature. (b) Electronic (dashed) and phonon-assisted (solid) intertube junction conductance of crossed CNTs for θ=40\theta=40^{\circ} (blue) and θ=70\theta=70^{\circ} (red). (c) shows phonon-assisted intertube junction conductance as a function of Fermi energy for θ=40\theta=40^{\circ} (dashed blue) and θ=70\theta=70^{\circ} (solid red). We used a bias of V=1V=1 mV in all panels and zero Fermi energy in panels (a) and (b).

The large-diameter limit of GinterG_{\rm inter} can be understood by considering the junction conductance between twisted graphene bilayers,Maarouf and Mele (2011) where the momentum conservation law of electrons in two Brillouin zones twisted with respect to each other by angle θ\theta governs the tunneling probability. However, in smaller-diameter CNTs, finite curvature plays an important role such that the finite contact area between the CNTs breaks the translation symmetry and helps relax the momentum conservation requirement. As a result, the conductance variations with angles diminish in magnitude, and a local minimum at 9090^{\circ} vanishes in small-diameter armchair CNTs, as shown in Figure 2a.

To better simulate the experimental situation, we relaxed the atomic positions at the junction using an atomistic valence force modelPerebeinos and Tersoff (2009) for the internal distortions in each CNT and the Lennard-Jones potential for atoms on adjacent CNTs. We set the van der Waals adhesion energy to Ea=60E_{\text{a}}=60 meV per atomZacharia et al. (2004) and the interlayer graphene bilayer spacing to h=3.35h=3.35 Å. The Ginter,elG_{\rm inter,el} of the relaxed structure is increased by an order of magnitude compared to the unrelaxed case; this arises from the fact that more atoms on neighboring CNTs get closer due to CNT bending in the contact area. Figure 3a depicts the Ginter,elG_{\rm inter,el} and Ginter,phG_{\rm inter,ph} of the relaxed structure versus θ\theta, showing that overall they are comparable in magnitude at room temperature. The θ\theta dependence of Ginter,phG_{\rm inter,ph} is weaker than that of Ginter,elG_{\rm inter,el} because phonons help relax the momentum conservation law. Note also that the SS and φ\varphi dependences of the conductance do not enter into consideration for the relaxed structure.

The temperature dependence of phonon-assisted conductance Ginter, phG_{\text{inter, ph}} is nearly linear in Figure 3b, reflecting the fact that low energy phonons are responsible for the phonon-assisted conductance, while the electronic intertube conductance is independent of temperature. The crossover temperature at which the two contributions become equal depends on the crossing angle. At angles near θ30\theta\sim 30^{\circ}, the phonon-assisted contribution dominates the intertube conductance for almost all temperatures, whereas at θ90\theta\sim 90^{\circ} the electronic contribution is larger even at room temperature, as shown in Figures 3a and b.

Finally, Figure 3c show the calculated Fermi energy dependence of Ginter,phG_{\rm inter,ph}. A small variation of the Ginter,phG_{\rm inter,ph} with the Fermi energy in Figure 3c reflects the fact that electron-phonon matrix elements have a weak momentum dependence and densities of electronic states are constant. In the case of electronic intertube conductance, we find that Ginter,elG_{\rm inter,el} is almost independent of the Fermi energy (not shown).

In summary, we demonstrated that the temperature dependence of CNT film conductivity depend both on intratube and intertube (or junction) electronic transport. The dominant conduction mechanism depends on the film preparation conditions and the electronic structure of the CNTs.Gao et al. (2021) In this study, we showed that in films with a long localization length, such as a chirality-enriched (6,6) CNT film, the temperature dependence of the junction conductance can give rise to the observed overall temperature dependence of the CNT film conductance over a wide temperature range. We evaluated the temperature-independent electronic junction conductance and the temperature-dependent phonon-assisted junction conductance using microscopic theory, as shown in Figure 3. Our numerical results show that both contributions are important at room temperature, and depending on the crossing angle, the phonon-assisted junction conductance can be significantly larger than the electronic intertube conductance. After averaging over the CNT crossing angle, our model accurately reproduces the experimental data over a wide temperature range with a single parameter fit known from the film device dimensions. We found that structural relaxation increases the electronic conductance by an order of magnitude, making calculated results more consistent with prior experimental results of junction conductance at low temperatures.Fuhrer et al. (2000a) While the overall qualitative agreement between our model and the measured film conductance is evident, the nearly perfect quantitative agreement is likely accidental because of the oversimplified model used to estimate the numerical value of the parameter CC in Equation (1). At the same time, we underestimate by a factor of three the intertube conductance measured on individual CNTs,Fuhrer et al. (2000a) which could be due to a lack of angle control in the experiment. However, there can be a more fundamental reason for the observed discrepancy: the role of short-range defects in CNTs on the electronic and phonon-assistant intertube junction conductance, which should be a subject of further theoretical and experimental studies. Our work on chirality-enriched films provides guidelines for reinterpretation and analysis of other CNT network film resistances with the interplay of CNTs intratube transport and intertube junction transport.

I Methods

Conductance of 2D film:

As the thickness of the film is much smaller than the averaged length of SWCNT (LcntL_{\text{cnt}}), the film can be considered to be a 2D film. Then, the number of junctions for each percolation path can be estimated by Njunc=L/(Lcnt/2)N_{\text{junc}}=L/(L_{\text{cnt}}/\sqrt{2}) assuming a saw-type 1D chain of length LL made of CNT rodes with crossing angle of 4545^{\circ}. The number of parallel percolation paths over width WW can be estimated as Npath=W/(Lcnt/2)N_{\text{path}}=W/(L_{\text{cnt}}/\sqrt{2}) and the number of 2D layers of CNT monolayers in a film of thickness tt can be estimated as Nlayer=t/(d+h)N_{\text{layer}}=t/(d+h). If each junction contributes a conductance of GjuncG_{\text{junc}}, then the total conductance due to the junctions in the film should be:

Gfilm, junc=NlayerNpathNjuncGinter=tWL1(d+h)GinterG_{\text{film, junc}}=N_{\text{layer}}\frac{N_{\text{path}}}{N_{\text{junc}}}G_{\text{inter}}=\frac{tW}{L}\frac{1}{(d+h)}G_{\text{inter}} (2)

If we consider the conductance contribution from the nanotubes themselves (i.e., the intratube conductance), we can estimate it as follows:

Gfilm, intra=NlayerNpathσintra2L=tWL1(d+h)LcntσintraG_{\text{film, intra}}=N_{\text{layer}}N_{\text{path}}\frac{\sigma_{\text{intra}}}{\sqrt{2}L}=\frac{tW}{L}\frac{1}{(d+h)L_{\text{cnt}}}\sigma_{\text{intra}} (3)

where 2L\sqrt{2}L is the average distance that electron travels in the random film to reach from one electrode to the other and σintra\sigma_{\text{intra}} is the intratube conductivity. The net resistance of the 2D film can be estimated as:Shim et al. (2015a)

Gfilm1\displaystyle G_{\text{film}}^{-1} =\displaystyle= Gfilm, intra1+Gfilm, inter1\displaystyle G_{\text{film, intra}}^{-1}+G_{\text{film, inter}}^{-1} (4)
=\displaystyle= LtW(d+h)(1Ginter+Lcntσintra)\displaystyle\frac{L}{tW}\,(d+h)\left(\frac{1}{G_{\text{inter}}}+\frac{L_{\text{cnt}}}{\sigma_{\text{intra}}}\right)

which leads to an estimate of the C=L(d+h)/(tW)C=L(d+h)/(tW) value in Eq. (1). Note that the CC value in this simple model does not depend on the CNT volume density as opposed to more sophisticated models Komori and Makishima (1977); Shim et al. (2015b); Tripathy et al. (2020). In our model, CNT volume density is given by nCNT=NlayerNpathNjunc/(tWL)=2/(Lcnt2(d+h))n_{\rm CNT}=N_{\text{layer}}N_{\text{path}}N_{\text{junc}}/(tWL)=2/(L^{2}_{\rm cnt}(d+h)) consistent with an estimate in Ref. Gao et al. (2021).

Electronic Intertube Junction Conductance:

The electronic conductance can be evaluated by Datta (2005)

Ginter, el=4πeV𝐤ssM𝐤ss(f(E𝐤s)f(E𝐤s+eV)),G_{\text{inter, el}}=\frac{4\pi e}{\hbar V}\sum_{\begin{subarray}{c}\mathbf{k}ss^{\prime}\end{subarray}}M_{\mathbf{k}s}^{s^{\prime}}(f(E_{\mathbf{k}s})-f(E_{\mathbf{k}s^{\prime}}+eV)), (5)

where the spin degeneracy is included. Here, 𝐤\mathbf{k} is the 2D electron wave vector in 2D Brillouin zone corresponding to periodically repeated supercell containing a single junction. s(s)s(s^{\prime}) labels the band index, ee, \hbar, VV, and ff are the elementary electronic charge, the reduced Plank constant, the applied voltage bias between layers, and the Fermi-Dirac distribution function, respectively. The intertube coupling M𝐤ssM_{\mathbf{k}s}^{s^{\prime}} is given by

M𝐤ss=|ψ𝐤s|Hel|ψ𝐤s|2δ(E𝐤sE𝐤s),M_{\mathbf{k}s}^{s^{\prime}}=|\matrixelement{\psi_{\mathbf{k}s}}{H_{\text{el}}}{\psi_{\mathbf{k}s^{\prime}}}|^{2}\delta(E_{\mathbf{k}s^{\prime}}-E_{\mathbf{k}s}), (6)

where ψ𝐤s|Hel|ψ𝐤s\matrixelement{\psi_{\mathbf{k}s}}{H_{el}}{\psi_{\mathbf{k}s^{\prime}}} is the tunneling matrix element for carrier scattering from state ψ𝐤s\psi_{\mathbf{k}s} to state ψ𝐤s\psi_{\mathbf{k}s^{\prime}} on different CNTs. E𝐤sE_{\mathbf{k}s} is the carrier energy, calculated by the tight-binding method with the first nearest-neighbor hoping t=3.1t=3.1 eV. HelH_{\text{el}} is the electron intertube interaction Hamiltonian due to the hopping between atoms on different CNTs. The hopping parameter is defined as (Perebeinos et al., 2012)

tij=texp(rijh0λz)exp[(ξijλxy)α],t_{ij}=t_{\perp}\exp(-\frac{r_{ij}-h_{0}}{\lambda_{z}})\exp[-\left(\dfrac{\xi_{ij}}{\lambda_{xy}}\right)^{\alpha}], (7)

where h0=3.35h_{0}=3.35 Å is the equilibrium interlayer distance, rijr_{ij} is the distance between atoms ii and jj, λz=0.6\lambda_{z}=0.6 Å, λxy=1.7\lambda_{xy}=1.7 Å, α=2.0\alpha=2.0, and t=0.4t_{\perp}=0.4 eV. 𝝃ij\boldsymbol{\xi}_{ij} is the projection of 𝐫ij\mathbf{r}_{ij} onto the layers’ bisector, and its absolute value is given by

ξij=rijsinϕ=rij1(𝝅ij𝐫ijπijrij)2.\xi_{ij}=r_{ij}\sin\phi=r_{ij}\sqrt{1-\left(\dfrac{\boldsymbol{\pi}_{ij}\cdot\mathbf{r}_{ij}}{\pi_{ij}r_{ij}}\right)^{2}}\,. (8)

ϕ\phi is the angle between 𝐫ij\mathbf{r}_{ij} and 𝝅ij\boldsymbol{\pi}_{ij}, 𝝅ij=𝝅i𝝅j\boldsymbol{\pi}_{ij}=\boldsymbol{\pi}_{i}-\boldsymbol{\pi}_{j}, where 𝝅i\boldsymbol{\pi}_{i} is determined by (Perebeinos and Tersoff, 2009)

𝝅=3𝐫ij×𝐫ik+𝐫ik×𝐫il+𝐫il×𝐫ijrijrik+rikril+rilrij,\boldsymbol{\pi}=3\dfrac{\mathbf{r}_{ij}\times\mathbf{r}_{ik}+\mathbf{r}_{ik}\times\mathbf{r}_{il}+\mathbf{r}_{il}\times\mathbf{r}_{ij}}{r_{ij}r_{ik}+r_{ik}r_{il}+r_{il}r_{ij}}, (9)

where jj, kk, and ll are the three nearest neighbors of ii. 𝝅i\boldsymbol{\pi}_{i} is perpendicular to the surface of the CNT, and its direction points toward the outside of CNT, and 𝝃ij\boldsymbol{\xi}_{ij} is perpendicular to 𝝅ij\boldsymbol{\pi}_{ij}.

Phonon-Assisted Intertube Junction Conductance:

The phonon-assisted conductance is computed by (Perebeinos et al., 2012)

Ginter, ph=4πeV𝐤,𝐤,s,s,μM𝐤sμ𝐤s[f(E𝐤s)f(E𝐤s+eV)].G_{\text{inter, ph}}=\frac{4\pi e}{\hbar V}\sum_{\begin{subarray}{c}\mathbf{k},\mathbf{k^{\prime}},s,s^{\prime},\mu\end{subarray}}M_{\mathbf{k^{\prime}}s^{\prime}\mu}^{\mathbf{k}s}[f(E_{\mathbf{k}s})-f(E_{\mathbf{k^{\prime}}s^{\prime}}+eV)]. (10)

Here, M𝐤sμ𝐤sM_{\mathbf{k^{\prime}}s^{\prime}\mu}^{\mathbf{k}s} is given by

M𝐤sμ𝐤s=\displaystyle M_{\mathbf{k^{\prime}}s^{\prime}\mu}^{\mathbf{k}s}={} |ψ𝐤s|He-phμ|ψ𝐤s|2[n𝐪μδ(E𝐤sE𝐤s+ω𝐪μ)\displaystyle|\matrixelement{\psi_{\mathbf{k}s}}{H_{\text{e-ph}}^{\mu}}{\psi_{\mathbf{k^{\prime}}s^{\prime}}}|^{2}[n_{\mathbf{q}\mu}\delta(E_{\mathbf{k^{\prime}}s^{\prime}}-E_{\mathbf{k}s}+\hbar\omega_{\mathbf{q}\mu}) (11)
(1+n𝐪μ)δ(E𝐤sE𝐤sω𝐪μ)],\displaystyle(1+n_{-\mathbf{q}\mu})\delta(E_{\mathbf{k^{\prime}}s^{\prime}}-E_{\mathbf{k}s}-\hbar\omega_{-\mathbf{q}\mu})],

where nn is the Bose-Einstein function, ψ𝐤s|He-phμ|ψ𝐤s\matrixelement{\psi_{\mathbf{k}s}}{H_{\text{e-ph}}^{\mu}}{\psi_{\mathbf{k^{\prime}}s^{\prime}}} is the electron-phonon matrix element, and ω𝐪μ\omega_{\mathbf{q}\mu} is the phonon frequency with wave vector 𝐪=𝐤𝐤\mathbf{q}=\mathbf{k}-\mathbf{k^{\prime}} and branch number μ\mu. To obtain the phonon frequency, the structure is relaxed using the BFGS algorithm of the quasi-Newton method.Avriel (2003) The electron-phonon Hamiltonian are extracted by expanding Equation (7) in atomic displacements corresponding to the phonon normal modes.

Acknowledgements.
We gratefully acknowledge funding from the Vice President for Research and Economic Development (VPRED), SUNY Research Seed Grant Program, and computational facilities at the Center for Computational Research at the University at Buffalo (http://hdl.handle.net/10477/79221). A.M. and J.K. acknowledge support by the Basic Energy Science (BES) program of the U.S. Department of Energy through Grant No. DE-FG02-06ER46308 (for preparation of carbon nanotube films) and the Robert A. Welch Foundation through Grant No. C-1509 (for structural and electrical characterization measurements).

References

  • Jariwala et al. (2013) D. Jariwala, V. K. Sangwan, L. J. Lauhon, T. J. Marks,  and M. C. Hersam, Chem. Soc. Rev. 42, 2824 (2013).
  • Okimoto et al. (2010) H. Okimoto, T. Takenobu, K. Yanagi, Y. Miyata, H. Shimotani, H. Kataura,  and Y. Iwasa, Advanced Materials 22, 3981 (2010).
  • Rao et al. (2018) R. Rao, C. L. Pint, A. E. Islam, R. S. Weatherup, S. Hofmann, E. R. Meshot, F. Wu, C. Zhou, N. Dee, P. B. Amama, J. Carpena-Nuñez, W. Shi, D. L. Plata, E. S. Penev, B. I. Yakobson, P. B. Balbuena, C. Bichara, D. N. Futaba, S. Noda, H. Shin, K. S. Kim, B. Simard, F. Mirri, M. Pasquali, F. Fornasiero, E. I. Kauppinen, M. Arnold, B. A. Cola, P. Nikolaev, S. Arepalli, H.-M. Cheng, D. N. Zakharov, E. A. Stach, J. Zhang, F. Wei, M. Terrones, D. B. Geohegan, B. Maruyama, S. Maruyama, Y. Li, W. W. Adams,  and A. J. Hart, ACS Nano 12, 11756 (2018).
  • Gao et al. (2020) W. Gao, N. Komatsu, L. W. Taylor, G. V. Naik, K. Yanagi, M. Pasquali,  and J. Kono, J. Phys. D 53, 063001 (2020).
  • Zorn and Zaumseil (2021) N. F. Zorn and J. Zaumseil, Applied Physics Reviews 8, 041318 (2021).
  • Iijima (1991) S. Iijima, Nature 354, 56 (1991).
  • Iijima and Ichihashi (1993) S. Iijima and T. Ichihashi, Nature 363, 603 (1993).
  • Bethune et al. (1993) D. S. Bethune, C. H. Kiang, M. S. de Vries, G. Gorman, R. Savoy, J. Vazquez,  and R. Beyers, Nature 363, 605 (1993).
  • Dresselhaus et al. (2001) M. S. Dresselhaus, G. Dresselhaus,  and P. Avouris, eds., Carbon Nanotubes: Synthesis, Structure, Properties, and Applications, Topics in Applied Physics No. 18 (Springer, Berlin, 2001).
  • Avouris et al. (2007) P. Avouris, Z. Chen,  and V. Perebeinos, Nat. Nanotechnol. 2, 605 (2007).
  • Avouris et al. (2008) P. Avouris, M. Freitag,  and V. Perebeinos, Nat. Photon. 2, 341 (2008).
  • Nanot et al. (2013) S. Nanot, N. A. Thompson, J.-H. Kim, X. Wang, W. D. Rice, E. H. Hároz, Y. Ganesan, C. L. Pint,  and J. Kono, in Handbook of Nanomaterials, edited by R. Vajtai (Springer, Berlin, 2013) pp. 105–146.
  • Weisman and Kono (2019) R. B. Weisman and J. Kono, eds., Optical Properties of Carbon Nanotubes: A Volume Dedicated to the Memory of Professor Mildred Dresselhaus (World Scientific, Singapore, 2019).
  • Nakanishi and Ando (2001) T. Nakanishi and T. Ando, J. Phys. Soc. Jpn. 70, 1647 (2001).
  • Buldum and Lu (2001) A. Buldum and J. P. Lu, Phys. Rev. B 63, 161403 (2001).
  • Maarouf and Mele (2011) A. A. Maarouf and E. J. Mele, Phys. Rev. B 83, 045402 (2011).
  • Fuhrer et al. (2000a) M. S. Fuhrer, J. Nygård, L. Shih, M. Forero, Y.-G. Yoon, H. J. Choi, J. Ihm, S. G. Louie, A. Zettl,  and P. L. McEuen, Science 288, 494 (2000a).
  • Fuhrer et al. (2000b) M. Fuhrer, A. K. Lim, L. Shih, U. Varadarajan, A. Zettl,  and P. L. McEuen, Physica E: Low-dimensional Systems and Nanostructures 6, 868 (2000b).
  • Perebeinos et al. (2012) V. Perebeinos, J. Tersoff,  and P. Avouris, Phys. Rev. Lett. 109, 236604 (2012).
  • Gao et al. (2021) W. Gao, D. Adinehloo, X. Li, A. Mojibpour, Y. Yomogida, A. Hirano, T. Tanaka, H. Kataura, M. Zheng, V. Perebeinos,  and J. Kono, Carbon 183, 774 (2021).
  • Shim et al. (2015a) W. Shim, Y. Kwon, S.-Y. Jeon,  and W.-R. Yu, Sci. Rep. 5, 1 (2015a).
  • Perebeinos and Tersoff (2009) V. Perebeinos and J. Tersoff, Phys. Rev. B 79, 241409 (2009).
  • Zacharia et al. (2004) R. Zacharia, H. Ulbricht,  and T. Hertel, Phys. Rev. B 69, 155406 (2004).
  • Komori and Makishima (1977) T. Komori and K. Makishima, Textile Research Journal 47, 13 (1977).
  • Shim et al. (2015b) W. Shim, Y. Kwon, S.-Y. Jeon,  and W.-R. Yu, Scientific Reports 5, 16568 (2015b).
  • Tripathy et al. (2020) S. Tripathy, B. Bose, P. P. Chakrabarti,  and T. K. Bhattacharyya, IEEE Transactions on Electron Devices 67, 5676 (2020).
  • Datta (2005) S. Datta, Quantum Transport: Atom to Transistor (Cambridge University Press, 2005).
  • Avriel (2003) M. Avriel, Nonlinear Programming: Analysis and Methods (Courier Corporation, 2003).