Abstract

We investigated properties of an armchair nanotube attached by specific numbers of diketone groups along the tube axis using density functional theory (DFT) calculations. The results from DFT calculations show that multiple diketone attachments into an armchair nanotube cleave the CC bonds along the tube axis, forming a large hole. Then, the six-membered rings surrounding the hole are planarized, and zigzag edges appear. Due to these geometrical changes, the functionalized armchair nanotubes exhibit properties similar to those in corresponding graphene ribbons with zigzag edges. For example, diketone-attached nanotubes have a spin-polarized ground state with frontier orbitals whose amplitudes are localized at diketone O atoms. As a consequence of the existence of the localized orbitals, unpaired electrons appear only on the diketone O atoms in an armchair nanotube.

1. Introduction

Honeycomb carbon sheets [1, 2] are contained in both single-walled carbon nanotubes [3, 4] and single-layer graphite, also called graphene [5]. The two carbon allotropes, which use the same building blocks, exhibit similar and unique electronic properties resulting from small band gaps [1, 2]. Although the unique properties primarily come from πœ‹ orbitals in the frontier orbital region, there are some differences in their electronic properties. The differences originate from the dimensionality of the nanocarbon materials. In fact, one-dimensional nanotubes exhibit metallic or semi-conducting character, depending on how the honeycomb sheet rolls up [6–9]. The properties of thin strips of graphene (graphene-ribbons) are strongly influenced by their width and edge structures (armchair- and zigzag-edges) [10–23]. In particular, zigzagnanoribbons have radical character, which is localized at their edges.

To demonstrate the ability to connect nanotubes and graphene [24–30], Kosynkin et al. successfully converted a nanotube into a nanoribbon via a chemical treatment [26]. In this intriguing experiment [26], KMnO4 oxidizes a nanotube and unzips the nanotube longitudinally to form a nanoribbon, as shown in Figure 1. In the proposed reaction mechanism, KMnO4 binds initially to a CC bond in the nanotube. After the binding, one diketone group is formed by cleaving the attached CC bond. Because the first diketone formation distorts the neighboring CC bonds, the distorted CC bonds become a suitable site for the second KMnO4 addition. After a sequence of oxidation reactions, the CC bonds of the tube are successfully cleaved along the tube axis. Although their proposal appears plausible and is consistent with the experimental findings, atomic-scale views of the reaction intermediates, which include a nanotube attached by multiple diketone groups, are still required.

For this purpose, computer simulations based on density functional theory (DFT) [31–38] can provide insights into the reaction mechanism. Although several computational studies have focused on unzipping nanotubes with and without H-atom termination [31–35], there are a few reports discussing those created by the addition of oxygen pairs [36, 37]. For example, Rangel et al. used DFT calculations to investigate the sequential addition of oxygen pairs to a model with a finite-length (5,5) tube. The model, which included 100 atoms, has a tube length of ~10 Å. Based on the tube length, the maximum number of oxygen pairs attached to the tube (𝑁op) can be as many as 5. According to the DFT calculations, cleaving the tube CC bonds by diketone formation was found for 𝑁op larger than 3. Although their DFT findings can partially explain the experimental results, the structural features of the oxygen-pair attached nanotubes are artificially affected by the model size. In fact, the separation between the carbon atoms attached by oxygen pairs is influenced by the tube edges. Furthermore, attached oxygen atoms cannot line up along the tube axis at regular intervals. Because these structural features have a strong impact on the electronic properties, detailed insight into the properties of diketone-attached nanotubes is still required.

To improve our understanding of the conversion of nanotubes into nanoribbons by multiple diketone attachments, we employed DFT calculations for a sufficiently long tube model and an infinite-length model using a periodic boundary condition (PBC). In the current study, we focused on how the properties of diketone-attached nanotubes change as a function of the degree of oxidation (i.e., the number of diketones formed on a nanotube surface (𝑁op)). It is expected that the diketone-attached nanotubes have properties in between a pristine nanotube and the resultant graphene-nanoribbon. Thus, we focused on how a nanotube becomes planar via a sequence of diketone formation and how the planarization affects the electronic properties.

2. Method of Calculation

The DFT calculations used in the current study are based on the Perdew-Wang (PW91) functional [39]. When some numbers of diketones are formed on a CC bond perpendicular to the axis of a (5,5) armchair nanotube, the CC bonds at the binding sites are proposed to be cleaved [26]. After cleaving the attached CC bonds along the axis, zigzag edges appeared at the binding sites. Previous theoretical studies indicate that a spin-polarized solution for the zigzag nanoribbons with attached ketone groups is energetically stable relative to its spin-unpolarized solution [18–20]. We performed both spin-polarized and spin-unpolarized calculations to obtain local minima of the diketone-attached nanotubes.

In the PBC calculations implemented in the Vienna ab initio Simulation Package (VASP v.4.6) [40], we employed a hexagonal supercell in Figure 2. The supercell contains 200 C atoms in the (5,5) nanotube and some number of oxygen pairs. In the current study the number of attached oxygen pairs (𝑁op) ranges from 1 to 7. The kinetic energy cutoff of the plane-wave basis set is 349.5 eV with ultrasoft Vanderbilt-type pseudopotentials [41]. We allowed full geometry relaxation in the axial direction, but separations between adjacent tubes are fixed at 15.5 Å (Figure 2) (In the PBC calculations (Figure 2), we relaxed cell parameter in c vector as well as position of included atoms. However two cell parameters in a and b vectors were fixed, and eventually intertube separations are fixed). The intertube separations were large enough to avoid significant intertube interactions [42–44]. The cell size is also large enough to avoid interactions between the guests located on the neighboring unit cells. We employed a 1Γ—1Γ—5π‘˜-point mesh for geometry optimization of the nanotube systems. In PBC calculations of metallic systems with no band gap, we sometimes faced a difficulty in converging their total energy. To avoid the difficulty, the convergence was enhanced in the VASP program using fractional occupancy generated by Gaussian broadening (𝜎GB=0.03 eV) of the one-electron energy levels. As a result, the energy levels close to the Fermi level have fractional occupations. Based on the PW91 calculations, the CC bond length in graphene was determined to be 1.421 Å, which is consistent with the experimental value [1, 2] suggesting that the method of our choice is appropriate for the present study.

In the finite-length model calculation performed with the Gaussian 03 program [45], we fully optimized a (5,5) nanotube terminated by H atoms (C230H20) plus 𝑁op oxygen pairs. The length of the tube model was approximately 27 Å, which was large enough to model the corresponding infinite-length tube [46–49]. We used the 6–31G* basis set [50] for the attached oxygen pairs as well as the neighboring 58 C atoms, and the 3–21 G basis set for all of the other atoms [51]. Similar to previous theoretical studies [18–20], single and triplet spin states of the diketone-attached nanotubes were obtained for the finite-length calculations.

3. Results and Discussion

3.1. Formation of One Diketone Group on a (5,5) Nanotube

As proposed by Kosynkin et al. [26], the diketone formation would be an indispensable reaction to cleave the tube CC bonds at the binding sites. In this section, we describe how DFT calculations were used to investigate relationships between the diketone formation and the CC bond cleavage. For this purpose, we constructed two types of initial geometries for the addition of two oxygen atoms into the (5,5) nanotube, which can be distinguished by whether the attached C–O bonds are longer or shorter. As a result of PW91 optimization, we obtained two local minima, as shown in Figure 3. Figure 3(a) displays the optimized structure for the formation of one diketone group on the (5,5) nanotube. As shown in Figure 3(a), the optimized CO bond length is 1.24 Å, which indicated the presence of two CO double bonds (The optimized geometry was obtained by using the initial geometry of the nanotube attached by one oxygen-atom pair with shorter C–O bonds). Then, the CC bond at the binding site is broken (CC separation, 2.73 Å). Another type of optimized structure was also obtained, and it is shown in Figure 3(b) (The optimized geometry was obtained by using the initial geometry of the nanotube attached by one oxygen-atom pair with longer C–O bonds). The optimized structure is 31.3 kcal/mol more unstable compared to the structure in Figure 3(a). A similar value was obtained in [37]. In the oxygen-pair addition shown in Figure 3(b), two single CO bonds (1.47 Å) are formed. Then, the OO bond (1.52 Å) stretches, but it does not completely break. The optimized structure retains the tube CC bond at the binding site.

A comparison between the two optimized structures elucidated the roles of the oxygen-pair addition in the CC bond cleavage. To understand the different chemical bonds between the two C2O2 networks, the valence concept was utilized. Conceptually, the two C2O2 networks are connected by an electron transfer, which is indicated by the arrows in Figure 4. The green arrow in Figure 4 indicates a single-electron transfer from the OO bond to the CO bond, while a blue arrow indicates a single-electron transfer from the attached CC bond to the CO bond. As a result of the electron transfer, the CC and OO bonds are cleaved, and two double CO bonds are formed (Figure 4). This process can be rationalized with the Mulliken charge populations, and the calculated values are also presented in Figure 4. Thus, the DFT results support the experimental prediction posited in [26].

3.2. A Series of Diketone-Attached Infinite-Length (5,5) Nanotube

In this section, we discuss a series of diketone-attached infinite-length (5,5) nanotube where the number of attached diketone groups (𝑁op) ranges from 2 to 7. Their local minima were obtained from spin-polarized and spin-unpolarized DFT calculations. Table 1 shows the relative energy between the spin-polarized and spin-unpolarized solutions of each diketone-attached nanotube (Δ𝐸state=𝐸polarizedβˆ’πΈunpolarized). The negative Δ𝐸state values in Table 1 suggest that a spin-polarized solution is energetically preferable relative to a spin-unpolarized solution. In particular, 3-and 4-diketone-attached nanotubes have substantial Δ𝐸state values (βˆ’4.2 andβˆ’6.2 kcal/mol, resp.).

Figure 5 provides an overview of the structures for the optimized 𝑁op-diketone-attached nanotubes in the spin-polarized ground state. For comparison, the optimized structures for the corresponding H-attached nanotubes are also displayed in Figure 6. In Figure 5, the attached oxygen (green) atoms are aligned along the tube axis irrespective of the number of attached oxygen atoms. It is important to note that the O···O separations along the tube axis are nearly constant (~2.5 Å) for larger 𝑁op values. Such geometrical features cannot be found by Rangel et al. [36]. As in the case of the 1-diketone-attached nanotube, the CC bonds break as the multiple diketone groups are formed on the nanotube. With an increase in 𝑁op, larger CC separations are found. Accordingly, a larger hole is generated in the 𝑁op-diketone-attached nanotube with a larger 𝑁op value, as shown in Figure 5. For example, the 7-diketone-attached nanotube has a hole of 19.4Γ—6.6 Å. Interestingly, the hole width (6.6 Å) is comparable to the diameter of the pristine (5,5) nanotube (6.9 Å). In terms of the hole formation, similarities were found in the H-attached nanotubes in Figure 6. However, the H-attached nanotube has a slightly smaller hole than the corresponding diketone-attached nanotube case. The different hole sizes can be reproduced in the finite-length calculations. According to Mulliken population analyses in the finite-length calculations, the attached oxygen atoms have substantial negative charges (βˆ’0.38e~βˆ’0.29e) (nanotubes attached by smaller number of oxygen pairs have more significant charges with negative signs). Between the paired oxygen atoms attached to the nanotubes, repulsive forces exert, forming the larger holes in the diketone-attached nanotubes than those in the corresponding H-atom-attached nanotubes where such negative charges do not appear on the attached H-atoms.

Because of the large-hole formation, the 7-diketone-attached nanotube cannot retain a tubular structure. However, some of the six-membered rings seem to be planarized and are shown as circles in Figure 5. To determine how the hole formation results in planarizing of the six-membered rings, we investigate in detail the structural features of the 2-, 3-, 5-, and 7-diketone-attached nanotubes, as shown in Figure 7. The six-membered rings attached by a ketone group are labeled to identify the binding sites. For example, the binding sites in the 7-diketone-attached nanotube have blue carbon atoms in the six-membered rings with labels ranging from β€œ3” to β€œ9.” The dependence of the CC bond cleavage on the location of the attached six-membered rings is displayed in Figures 7 and 8(a). In Figures 7 and 8(a), the two carbon atoms are significantly separated in the nanotubes with larger numbers of attached diketone groups. In particular, the longest CC separation was found in the middle of the hole in the 7-diketone-attached nanotube and is labeled β€œ6” in Figure 7. On the other hand, smaller CC separations were found in the nanotubes attached by a smaller number of diketone groups (𝑁op=2 or 3).

Planarization of the six-membered rings surrounding the hole is quantified in Figure 8(b) using the dihedral angles (D) of the four C atoms (a, b, c, and d (aβ€², bβ€², cβ€², and dβ€²) in Figure 7) neighboring the attached site. It is important to note that the pristine (5,5) tube and the planer graphene have absolute D values of 19.8 (Dihedral angles of a, b, c, and d atoms have positive D values, whereas those of aβ€², bβ€², cβ€², and dβ€²β€² have negative values.) and 0, respectively, as shown in Figure 9. Thus, the dihedral angles are useful parameters for determining if a six-membered ring is planar in the diketone-attached tube. In fact, their D values would approach 19.8 degrees, when the diketone-attached nanotube exhibits a tubular structure. In contrast, the planarized six-membered rings appear in the diketone-attached nanotube with D values close to 0.

As shown in Figure 8(b), we observe an increase in the D values near the binding sites of the 2-diketone-attached nanotube. In addition to the D values larger than 19.8 degrees, two six-membered rings have a smaller D value. Much smaller D values are observed for the 3-diketone-attached nanotube where two diketone-attached rings, labeled 4 and 6, have D values of approximately 10 degrees. Decreasing D values indicate that the curvature of the ring is substantially reduced. When larger numbers of diketone groups are attached on a tube, nearly planar six-membered rings are present. For example, some six-membered rings have D values ranging from 0 to 10 degrees in the 7-diketone-attached nanotube. More importantly, Figures 8(a) and 8(b) exhibit a clear correlation between the planarity of a six-membered ring and the dissociation of a CC bond via diketone formation. When carbon atoms are more than 5 Å apart, their D values are less than 10 degrees. Because the planarized six-membered rings have zigzag edges, the edge structures perturb the electronic properties of the armchair nanotube.

3.3. Electronic Properties of a Finite-Length Nanotube Attached by Diketone Groups

In this section, we discuss the electronic properties of the 7-diketone-attached nanotube using a finite-length model. In the finite-length calculations, we obtained results similar to the infinite-length cases in terms of their geometrical features. Furthermore, the 7-diketone-attached nanotube in the triplet state is 5 kcal/mol more stable than that the singlet state. Because of the stability of the triplet state, the 7-diketone-attached nanotube has radical character. According to the DFT calculations, the unpaired electrons are primarily localized on the diketone O atoms, as shown in Figure 10, where isosurface 𝛼- and 𝛽-spin densities are given by pink and blue, respectively. To understand the unique radical character, we investigated representative 𝛼-spin orbitals in the frontier orbital region of the 7-diketone-attached nanotube in Figure 11. The HOMO-2 and HOMO-4 are only 0.19 and 0.34 eV below the HOMO. These orbitals have amplitudes localized on diketone-O atoms, which is not the case for pristine and H-attached [32–35] nanotubes. Similar orbital features are observed in the magnetic orbitals of ketone-attached nanoribbons [20]. Therefore, the multiple diketone attachments yield an armchair nanotube orbital properties similar to those in the corresponding zigzag nanoribbon.

The patterns of the localized orbitals are reminiscent of an iso-distant one-dimensional chain, whose molecular orbital can be given by the following equation:πœ“π‘—=ξ‚™2𝑛+1π‘›ξ“π‘ž=1ξ‚΅sinπœ‹π‘—π‘žξ‚Άπœ’π‘›+1π‘ž,(1) where 𝑛 is the number of atoms contained in the chain, πœ’π‘ž is an atomic orbital on the π‘žth atom, and Ψ𝑗 is the 𝑗th molecular orbital. As mentioned above, the diketone-attached nanotube has two chains with seven oxygen atoms occurring at nearly regular intervals. Thus, their orbital features can be roughly understood from (1). In these systems, πœ’π‘ž corresponds to oxygen p (p(O)) orbitals parallel to or perpendicular to a plane containing a newly formed zigzag edge, which is depicted schematically in Figure 12. For example, the HOMO-2 consists of perpendicular p(O) orbitals, whereas the HOMO-4 consists of parallel p(O) orbitals. It is interesting to note that the localized orbitals have nodes resulting from the sine terms in (1). For example, the HOMO-2 has three nodes along an oxygen-chain, and the HOMO-4 has one node.

Based on the p(O) orbital orientation in the localized orbitals, the two types of orbitals have different distributions of πœ‹ orbitals in the tube moiety. In fact, the HOMO-2 has carbon p orbitals slightly distributed into the tube. On the other hand, the HOMO-4 does not have p orbitals on the tube C atoms, except for carbon atoms that are next-nearest neighbor to a diketone oxygen atom. The different behaviors can be understood from Figure 12 [52]. When parallel p(O) orbitals are perpendicular to the tube πœ‹ orbitals, they cannot mix with the tube πœ‹ orbitals. As a result, localized orbitals containing parallel p(O) orbitals do not exhibit amplitudes on the tube moiety. In contrast, perpendicular p(O) orbitals can interact with the tube πœ‹ orbitals. Because of the interactions, some amplitudes appear in the tube moiety in the frontier orbital containing perpendicular p(O) orbitals. Consequently, the planarization of the six-membered rings with multiple diketone attachments plays a significant role in determining the orbital features in the frontier orbital region.

Detailed analyses of the frontier orbitals revealed that HOMO-4 is occupied by an 𝛼 electron, but it is not occupied by a 𝛽 electron. Therefore, the 𝛼 HOMO-4 is one of the key orbitals that determines the spin density distribution in the 7-diketone-attached nanotube. However, there are many frontier orbitals with amplitudes that exist only on the diketone groups. In fact, nine 𝛼-orbitals lie less than 0.5 eV below the HOMO. These orbitals can be distinguished by the number of nodes as well as the orientation of the p(O) orbitals involved (e.g., perpendicular or parallel orientation). Because many magnetic orbitals exist in the frontier orbital region, we cannot identify all of the orbitals that determine the spin density distributions in the diketone-attached nanotube. However, localized orbitals resulting from p oxygen atoms, such as the HOMO-4, can contribute to the localized radical character.

4. Conclusions

Density functional theory (DFT) PW91 calculations were employed to investigate the properties of a (5,5) nanotube attached to some number of diketone groups. According to the DFT calculations, the formation of two double CO bonds on the nanotube cleaves one CC bond at the attachment site. For the high degree of diketone formation on the armchair nanotube, the CC bonds at the binding sites lengthen significantly, and, eventually, a large hole is formed. Then, the six-membered rings surrounding the large hole are substantially planarized. Because the planarized six-membered rings have zigzag edges attached to diketone groups, the functionalized armchair nanotubes exhibit electronic properties similar to those in nanographene ribbons with zigzag edges. In fact, orbitals with amplitudes primarily localized at diketone groups appear in the frontier orbital region. Based on the orbital features, diketone-attached nanotubes have a spin-polarized ground state, with the unpaired electrons localized on the oxygen atoms

Acknowledgments

The project is partially supported by a Grant-in-Aid for Young Scientists (B) from the Japan Society for the Promotion of Science (JSPS) (T. Yumura at the Kyoto Institute of Technology) (no. 22710088) and a Grant-in-Aid for Scientific Research (C) from JSPS (T. Yumura at the Kyoto Institute of Technology) (no. 23560934).