Using Bogoliubov-de Gennes formalism for inhomogeneous systems, we have studied superconducting properties of a bundle of packed carbon nanotubes, making a triangular lattice in the bundle's transverse cross-section. The bundle consists of a mixture of metallic and doped semiconducting nanotubes, which have different critical transition temperatures. We investigate how a spatially averaged superconducting order parameter and the critical transition temperature depend on the fraction of the doped semiconducting carbon nanotubes in the bundle. Our simulations suggest that the superconductivity in the bundle will be suppressed when the fraction of the doped semiconducting carbon nanotubes will be less than 0.5, which is the percolation threshold for a two-dimensional triangular lattice.

Single wall carbon nanotubes (SWCNTs) represent a unique class of quasione-dimensional nanoscale systems exhibiting various interesting phenomena. Among other exciting features, it was demonstrated that individual single wall carbon nanotubes may have intrinsic superconducting properties [1]. However, because of their extremely small diameter (just few nanometers), and thus strongly one-dimensional character, the superconducting order parameter may have significant “phase slips” due to thermal and quantum fluctuations, leading to a finite conductivity in the system below the critical temperature [2, 3]. Carbon nanotubes can also form bundles and ropes [4], with tens and hundreds of individual SWCNTs in the bundle, coupled to each other by dispersive Van der Waals forces. Such kind of system may exhibit reduced “phase slips” effects and as a result, much stronger conductivity drop below the critical temperature. The overall length of a SWCNT in the bundle also plays a significant role. For example, reducing the bundle’s length to 300 nm may destroy the superconductivity in the system due to increasingly high quantum fluctuations [2]. Generally speaking, for nanoscale systems with the quantum level spacing approaching the superconducting gap energy , the superconductivity vanishes [5].

It is expected that doping of SWCNTs in a bundle by, for example, boron, may significantly improve their superconducting properties [6]. At a proper level of doping, the Fermi level may be at a one-dimensional singularity of the energy spectrum that gives a higher density of states (DOSs), that will lead to a higher critical temperature . In particular, we assume here that such kind of mechanism of doping enhanced may be much better pronounced in the case of semiconducting SWCNTs, which may have higher DOS due to lower energy of the Van Hove singularities. This is in contrast to metallic SWCNTs, where singularities in the DOS are at higher energies and can be filled at much higher doping levels only (according to the Kataura plot [7]). This means that a bundle consisting of doped semiconducting nanotubes could be a much better superconductor compared to a bundle made of metallic SWCNTs.

However, synthesis of SWCNTs by currently known methods usually results in a mixture of semiconducting and metallic nanotubes. Since the nanotubes after the synthesis initially are not doped (or unintentionally slightly p-type doped, e.g., by oxygen of atmosphere), those are only metallic tubes, which may have superconducting transition, while semiconducting tubes will be “diluting” superconductivity in the bundle by the inverse proximity effect [8]. Upon doping (i.e., by electrochemical methods), the semiconducting tubes can become superconducting with a higher superconducting gap and thus a higher than in metallic nanotubes.

One may try to estimate a spatially averaged order parameter and the corresponding effective critical temperature for a bundle consisting of a mixture of these two types of SWCNTs. From an experimentalist’s point of view, it is even more important to solve a bit more complex problem: for a given fraction of doped semiconducting SWCNTs in the bundle and the experimentally determined critical temperature , to predict the critical temperature for a bundle, consisting only of doped semiconducting SWCNTs. It will be also interesting to know whether it is possible to obtain the critical temperature higher than in other carbon based nanostructures, like in alkali metal doped fullerenes.

Spatial variations of the superconducting order parameter are significant for nanoscale systems, including nanotubes [9, 10]. In this work, we use a microscopic theory based on inhomogeneous Bogoliubov-de Gennes equations to establish how the superconducting properties of a bundle depend on the fraction of doped semiconducting nanotubes. We assume that the nanotubes in the bundle are approximately of the same radii and tightly packed making a triangular lattice in the bundle’s transverse cross-section, with the primitive vectors , . Here are the unit basis vectors, and is the average intertube distance. The lattice can be enumerated by indexes , which correspond to the position of a nanotube , but in this work we prefer to enumerate nanotubes in a bundle using a single index through the mapping , where is the number of nanotubes in the raw. The doped semiconducting nanotubes are assumed to occupy the fraction of the sites in the bundle, and the metallic nanotubes occupy the fraction of the sites. A full three-dimensional description of a bundle taking into account the band structure and chirality of individual nanotubes would make simulations of the Bogoliubov-de Gennes equations too complex. On the other hand, assumption about the translation invariance along the longitude direction of the nanotubes would make the simulations marginally simpler, because one still needs to deal with quantum problem in three dimensions, but at the same time, this assumption would throw away the band structure difference between semiconducting and metallic nanotubes. In this study, we adopted a simplified two-dimensional picture which allows us to repeat our simulations several dozens of times for different arrangements of the nanotubes in the bundle. In the two-dimensional picture, the conduction electrons can stay in a nanotube or can hop to the neighboring nanotubes (sites). In principle, there may be three different hoping constants, with the hopping matrix elements equal to either ,  ,  and , corresponding to the hopping between metallic-metallic , metallic-semiconducting , or semiconducting-semiconducting nanotubes. Moreover, these parameters may significantly fluctuate from one site to another, due to variations of the intertube distance or the mismatch between SWCNTs of different radii and chirality, but in this work, we neglect all these fluctuations. In the superconducting regime, Cooper pairs can be formed within nanotubes and can also hop from one tube to another.

For the description of the system, we utilized the Hamiltonian form: where the quantum-mechanical operator creates an electron on site (using the single indexing), the operator eliminates an electron from the site , and represents the electron density on site with the spin polarization . The electron spin, , can point up or down. is the onsite interaction potential. This term in a case of attractive interaction may lead to pairing in the nanotube . is a strength of the coupling between electrons localized at neighboring tubes and .

Using the Bogoliubov transformation, which diagonalizes the Hamiltonian equation (1), we arrive to inhomogeneous Bogoliubov-de Gennes equations for the quasiparticle amplitudes on the lattice sites [11]: where the kinetic operator and the superconducting order parameter can be represented as where are the nearest neighbor vectors for a triangular lattice, is the mean field (Hartree) potential, is the chemical potential. is the conventional, -type order parameter. One should solve (2) together with the self-consistency conditions as follows: where the pairing strength may depend on the type of CN at the sites and . The -type order parameter (within a given nanotube ) is simply Note that the summation in (4), (5) is done over the positive eigenvalues only. The amplitudes obey the constraints for any (normalization) and for any , being the site index of the triangular lattice.

Here, we assumed the same constant hopping parameter between any type of nanotubes. In this work, we assumed that the pairing may only happen between electrons in the same nanotube, therefore neglecting a much weaker pairing mechanism between neighboring nanotubes. In principle, a weak attraction mechanism may stimulate the formation of a Cooper pair with one electron in one nanotube and the second electron in one of its nearest neighbors. This may result in coexisting order parameters in the system. The coexistence of order parameters with different symmetries, was studied, for example, for uranium-based superconducting materials [12]. The possibility of such pairing on a triangular lattice may result in unconventional superconducting properties. For example, a 2D triangular lattice was recently considered as a test-bed for a possibility of -wave spin-triplet superconductivity [13].

We studied how the spatially averaged superconducting order parameter depends on the fraction of doped semiconducting SWCNTs in a bundle at different temperatures. For this purpose, we generated realizations for a given number of randomly placed semiconducting nanotubes in the bundle. Other nanotubes in the bundle are assumed to be metallic. To model different superconducting pairing strengths for different types of nanotubes, we set for the doped semiconducting nanotubes and for the metallic ones, so the doped semiconducting nanotubes have a higher critical temperature. For simplicity, we set (half-filled band). In our simulations we considered 16 × 16 nanotubes in the bundle, forming a triangular lattice.

The results of simulations are shown in Figure 1. At , the spatially averaged order parameter scales approximately as a square root of the fraction of doped semiconducting SWCNTs in the bundle . Here and are the order parameters for pure metallic and doped semiconducting bundles. Note the convexity of the dependence. However, at finite temperatures, which are between the critical temperature of a pure metallic and pure doped semiconducting bundles, the averaged order parameter vanishes much faster with the decreasing of . For example, at , which is close to the critical temperature of a pure semiconducting SWCNT bundle , the order parameter decreases very quickly with the decreasing of and almost vanishes at (see Figure 1).

In Figure 2, we plot as a function of temperature for several values of . One can clearly see how the order parameter vanishes above the critical temperature. Note that with the lowering of the fraction of semiconducting nanotubes , the temperature dependence of the order parameter shows less pronounced phase transition because of the “dirty” nature of inhomogeneous spatial distribution of the pairing properties that is similar to dirty superconducting transition in case of large concentration of impurities and in alloys. We also would like to emphasize that the systems under consideration are finite; therefore, the order parameter vanishes not abruptly as in the BCS theory for an infinite system, even for the homogeneous case . Because of the presence of the temperature tails in Figure 2, we define the critical temperature as a temperature at which the order parameter drops to just 5% of its value at .

Using the data plotted in Figure 2, we calculated how the critical temperature depends on the concentration of doped semiconducting SWCNTs. We used  eV to be consistent with the data in [1, 6], so will correspond to K. In Figure 3, we plot the averaged critical temperature as a function of the fraction of semiconducting CNs. First, one can clearly see a nonlinear dependence of the . Second, in Figure 3 one can observe that for , the critical temperature changes with the change of with a steeper slope than for . For , the weaker superconducting metallic CNs are arranged in bigger islands, which surround and isolate doped semiconducting nanotubes from each other. The inverse proximity effect of low metallic tubes will significantly reduce the averaged order parameter in the system. Note, corresponds to the site percolation threshold on the plane triangular lattice [14]. If one would go beyond the mean field description and would take into account the phase fluctuations of the order parameter, the dependence would be even steeper, because of the enhanced phase fluctuations in relatively well-isolated doped semiconducting nanotubes, surrounded by metallic nanotubes. For , this will result in a rapid decrease of the spatially averaged order parameter.

According to our model for an optimally doped bundle consisting of 100% semiconducting SWCNTs, the should increase to the unsuppressed of K (see Figure 3). The effect of suppression similar to the one discussed here has been observed in alkali doped fullerene molecular alloys of [15], and adding non-superconducting component, that is, molecules, which do not show any superconducting pairing due to symmetry reasons and probably due to weaker electron-phonon coupling, strongly suppressed from 19 K in 100% , that is, in to K in 20% substituted alloy. The experiments with selectively separated metallic and semiconducting SWCNTs which now become available by new methods of effective separation will allow to check the validity simple model of presented here and to clarify the role of quantum fluctuations, which has not been accounted here.

One has to note that even for relatively small bundles 16 × 16, the physical properties have relatively small variations for different realizations of the spatial distributions of semiconducting nanotubes. To support this observation, we plot a whole set of the spatially averaged order parameters for 50 random configurations of the nanotubes for different temperatures. We choose the equal number of the metallic and semiconducting nanotubes in the bundle . In Figure 4, one can see that there are relatively limited variations of the order parameter corresponding to different configurations of the nanotubes. The overall dispersion is limited to approximately at zero temperature and becomes relatively significant at temperatures approaching the critical temperature. For experimentalists, this means that one needs to take extra care about the measurements of small samples near the critical temperature.

In Figure 5, we plot a particular realization of the spatial distribution of the superconducting order parameter (in units of ) at the percolation regime () and zero temperature. In Figure 5, red dots mark the triangular lattice of a 16 × 16 bundle. One can see the formation of big clusters of nanotubes with the same superconducting properties. For the measurements of the transverse conductivity of such bundles at finite temperatures, one should expect to see a dramatic drop of the conductivity for , that is, above the percolation threshold for the triangular lattice. At this concentration, the semiconducting nanotubes will likely create a connected network within a finite bundle of nanotubes. We propose that this prediction will be tested experimentally.

In conclusion, we introduced a mean field microscopic model to describe superconductivity in a bundle of a mixture of carbon nanotubes of different superconducting properties. We have studied the dependence of the spatially averaged superconducting gap on the fraction of doped semiconducting SWCNT (with a higher pairing strength) in the bundle at different temperatures. Note that for inhomogeneous nanoscale systems, the dependence for different concentrations may be nonlinear, as a manifestation of the breakdown of the bulk BCS theory. Indeed, our calculations of show that for the dependence has less steep slope. The reason is that the bundle is a highly inhomogeneous system with a pronounced inverse proximity effect. At , below the site percolation threshold for a 2D triangular lattice, the bundle can be seen as a collection of finite islands of “good” superconductors (doped semiconducting nanotubes), diluted by normal or weakly superconducting material (metallic nanotubes). Such islands demonstrate significantly the suppressed superconductivity, even in the mean field description, due to the enhanced inverse proximity effect.

Note that our mean field BdG model is unable to predict and properly describe quantum phase fluctuations of the order parameter in isolated doped semiconducting nanotubes, where the superconductivity will be suppressed even stronger. Future research using, for example, Ginzburg-Landau inhomogeneous equations [2, 1618] is necessary to describe such kind of effects. Because the dynamics of Cooper pairs in doped carbon nanotubes can be more close to the diffusive regime, the Usadel equations can be applied to calculate the finite conductivity at [19].