Next Article in Journal
Recent Advances of Upconversion Nanomaterials in the Biological Field
Previous Article in Journal
Preparation and Characterization of Sludge-Based Magnetic Biochar by Pyrolysis for Methylene Blue Removal
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Essential Electronic Properties of Silicon Nanotubes

1
Department of Physics/QTC/Hi-GEM, National Cheng Kung University, Tainan 701, Taiwan
2
Center of General Studies, National Kaohsiung University of Science and Technology, Kaohsiung 811, Taiwan
*
Author to whom correspondence should be addressed.
Nanomaterials 2021, 11(10), 2475; https://doi.org/10.3390/nano11102475
Submission received: 26 August 2021 / Revised: 16 September 2021 / Accepted: 20 September 2021 / Published: 22 September 2021

Abstract

:
In this work, the various electronic properties of silicon nanotubes (SiNTs) were investigated by the density functional theory. The cooperative and competitive relationships between the chiral angle, periodic boundary conditions, and multi-orbital hybridizations create unusual narrow gaps and quasi-flat bands in the ultra-small armchair and zigzag tubes, respectively. The features varied dramatically with tube radii. Armchair SiNTs (aSiNTs) have an indirect-to-direct band gap transition as their radius is increased to a particular value, while zigzag SiNTs (zSiNTs) present a metal-semiconductor transition. The projected density of states was used to elucidate the critical transitions, and the evolution of p and s orbital mixing states during the process are discussed in detail. The information presented here provides a better understanding of the essential properties of SiNTs.

1. Introduction

Carbon nanotubes (CNTs, single wall) contain a honeycomb lattice with various nanometer diameters and are regarded as a quasi one-dimensional (1D) tubular material [1,2,3,4,5,6,7]. The curvature effect leads to the misorientation of 3 p z orbitals and induces significant s p 3 hybridizations on the cylindrical surface, unlike the pure s p 2 hybridization of planar graphene. Therefore, CNTs act as metallics or semiconductors depending on their diameter and chiral vector. The large variety of electronic properties in different CNT structures ensure their potential application in many areas [8,9,10,11,12,13,14,15,16,17,18,19,20,21].
Silicon is a group-IV element with a similar electronic configuration to carbon atoms. However, the size of Si compared to C causes it to fill the s p 3 hybridization state preferentially. Early on, Fagan and Baierle et al. proposed hypothetical silicon nanotubes and studied their electronic and structural properties [22]. They considered two Si nanotube geometry types: hexagonal (h-NTs) and gear-like nanotubes (g-NTs). The two types differ in their Si atomic hybridization, where one is a rolled planar silicene (a monolayer sheet of Si atoms) and the other a buckled silicene that forms a nanotube. At the same time, Si g-NTs are more energetically stable than Si h-NTs [23].
SiNTs were first successfully synthesized in 2001 [24]. Different SiNT growth processes have been reported since then [25,26,27,28,29,30,31]. Sha et al. used a nano channel Al 2 O 3 [26] to fabricate SiNTs by chemical vapor deposition (CVD) [26,32,33]. After that, a variety of techniques have been used to create SiNTs, such as molecular beam epitaxy (MBE) [27], hydrothermal synthesis [29], and gas phase synthesis [30,34]. Ishai and Patolsky demonstrated that the diameter and wall thickness could be controlled between 1.5 to 500 nm [35]. The following are several possible applications for SiNTs. Yoo et al. used a surface sol-gel reaction to obtain SiNTs for energy storage devices [36]. Some other groups enhanced SiNT fabrication techniques to improve their properties for biological, chemical, and electrical applications in sensors, nanoscale electronics, and optoelectronics [37,38,39].
The electronic properties of SiNTs were theoretically studied by using the tight-binding model [23,40,41,42,43,44,45] and the first-principles calculations [22,46,47,48]. Most papers focused on the conducting properties, i.e., the band gap variations, but the detailed geometric parameters, a whole feature of energy bands, spatial charge density, the sophisticated projected density of states (PDOSs), and their dependence on chirality and diameter are still lacking. In this work, we used the first-principles density functional theory (DFT) under the Vienna ab initio Simulation Package (VASP) to thoroughly explore the essential properties of SiNTs. A whole feature of energy bands is discussed, including two groups of σ parabolic bands, frequent anti-crossings between the σ and π bands, a peculiar quasi-flat band in ultra-small zSiNTs near the Fermi level, the initial/final energies of σ and π bands and their mixing strongly dependent on the tube radius and chirality, etc. The above features of energy bands were not presented or discussed well before. Further, we display the spatial charge density distribution to clearly show the transition from multi-orbital to single-orbital hybridizations when increasing the radius. A significant difference in the transition process between armchair and zigzag SiNTs exists. The sophisticated PDOSs reveal the equal and nonequal contributions of p x and p y orbitals to the highest valence band (HVB) edges with s p / s p 2 / s p 3 bonding in the indirect-/direct-gap semiconductors for various radii. The above information is helpful to understand the important features in the energy bands, such as the unusual small energy gaps, nonmonotonous energy dispersions, and the existence of quasi-flat and anti-crossing bands.

2. Materials and Methods

The first-principles calculations based on the density functional theory (DFT) were carried out using the Vienna ab initio Simulation Package (VASP) [49,50]. The exchange-correlation energy derived from many-particle Coulomb interactions was evaluated using the Perdew–Burke–Ernzerhof functional (PBE) [51] under generalized gradient approximation. The one-dimensional periodic boundary condition was set along the y-direction, and the vacuum distance in the x- and z-space was set to 15 Å to make sure no interaction with neighboring SiNTs. The convergence of the Helmann–Feynman force was less than 0.01 eV/Å during ionic relaxations. The maximal cutoff energy of the wave function expanded by the plane wave was 500 eV. The Brillouin zone (BZ) and geometric optimization were implemented with 1 × 1 × 15 k-point mesh and 1 × 1 × 500 k-point mesh for the PDOSs calculations via the Gamma scheme.

3. Results

3.1. Geometric Structures

SiNTs are akin to rolling up a silicene sheet from the origin to the lattice vector (chiral vector) C h = n a 1 + m a 2 ≡ (n,m), where a 1 and a 2 are primitive hexagonal lattice vectors. We focus on axial and diagonal chiral vectors, namely armchair (n,n) and zigzag (n,0) types, for which cross-section and side views are presented in Figure 1. The buckling, enhanced by s p 3 hybridization, results in a gear-like structure on the cylindrical surface. Two inequivalent silicon atoms have smaller and larger distances to the tube axis, R 1 and R 2 , respectively. The periodic length of the 1D unit cell is 3 u 2 and 3 u 2 for armchair and zigzag, respectively, where u 2 is the Si-Si bond length “along” the axis.
The optimized geometric parameters, including the optimal lattice constant, the ground state energy E 0 , Si-Si bond lengths u 1 (“perpendicular” to the axis) and u 2 (“along” the axis), radii R 1 and R 2 , and buckling distance Δ z, are summarized in Table 1 for the armchair (n,n) and zigzag ( n ,0) tubes with n ∈ (2,25) and n ∈ (4,21), respectively. The reported radius of the SiNTs was adjustable in the range of 7.5–25 Å using ultra-high vacuum chemical vapor deposition (UHV-CVD) [35]. The following are the results of R 1 ( R 2 ) in the range of 1.731–26.428 (1.826–12.659) and 2.688–26.617 (2.881–13.129) Å for aSiNTs and zSiNTs, respectively. The buckling distance Δ z, determined by R 2 R 1 , is inversely proportional to the radius, roughly. In addition, the difference between the two bond lengths u 1 and u 2 (in units of the standard Si-Si bond length, u 0 2.28 Å) reflects the geometric distortion. When a silicene sheet is rolled up to form a tube, it could induce stress ( u 1 , 2 < 1 ) and strain ( u 1 , 2 > 1 ) between Si atoms. Table 1 shows more stress in the armchair tubes and strain in zigzag tubes. The ground state energy E 0 decreases as the radius increases according to the reduction in the curvature energy [52] and buckling.

3.2. Energy Bands

The aSiNT band structures with chiral vectors of (2,2), (5,5), (10,10), (15,15), (20,20), and (25,25) are shown in Figure 2a–f. There are two high-symmetry points in the 1D BZ of the nanotube, namely, the Γ point at k = 0 and X point at k = ± π / a , where a is the lattice constant of the nanotubes. Another crucial point is located between Γ and X, namely K at k = 2 π /3a. According to the BZ folding scheme, the Dirac cone at (2 π /3a, 2 π / 3 a) in the 2D BZ of hexagonal monolayers is folded into the K point (k = 2 π /3a) of the 1D BZ [53,54,55]. Therefore, the Γ and K points significantly influence the magnitude of band gaps. Unlike the metallic armchair CNTs, semiconducting behaviors are found in all aSiNTs. The low-energy bands of ultra-small aSiNTs are strongly dependent on the radius (Figure 2a–c). Si(2,2) presents an indirect energy gap ( E g ) of about 0.35 eV. The HVB at point Γ is associated with the misorientation of the 2 p z orbitals, which induces hybridization of the p and s states. On the other hand, the lowest conduction band ( L C B ) is between points Γ and K. The conduction and valence bands are asymmetric to E F = 0 due to low symmetry in the geometric structure (the deformed honeycomb structure). The large curvature suppresses the π and π * bands. Accordingly, it is impossible to use single-orbital bonding (2 p z orbital) to investigate the essential low-energy properties of the ultra-small tubes. When the radius increases, a valence and a conduction π band approach the Fermi level, leading to a direct gap around the K point (with a slight red shift with respect to the exact K point). The low-energy bands (15,15)–(25,25) have similar features, including a tiny direct gap close to the exact K point and the symmetric π and π * bands against the E F = 0, as shown in Figure 2d–f. Two groups of σ parabolic bands are initiated at the Γ point ranging from E v = −1.25 to −2.7 eV, associated with the σ -bonds between the neighboring p x and p y orbitals, respectively. Their non-degenerate energies arise from the different boundary conditions along the circumference and tube axis. The deep initial energy of the σ bands (≈−1.25 eV) reveals strong p-p bondings on a relatively flat surface. These properties contrast with the ultra-small tubes, e.g., Si(2,2), where the weak p-s bonding states have almost zero energy. The σ bands have frequent anti-crossings with the π bands in the E v = 2 2.5 eV region, where a lot of extra band-edge states and some quasi partial flat bands exist. This area may have exceeding van Hove singularities (discussed later). The anti-crossing behaviors denote that the buckling and curvature cause complex orbital hybridizations. The s bands, associated with the 3 s orbitals, come into existence at a deeper energy ( E v < −3 eV) due to the largest ionization energy (the lowest one-site energy in the tight-binding model [56]).
The zSiNTs energy dispersions of Si(4,0)–Si(21,0) are shown in Figure 3a–i. The zSiNTs present a wider variety of low-energy properties than the aSiNTs when changing the radius. The chiral indexes (n,0) for n = 4–9 exhibit metallic features, while n = 10–21 are semiconductors [45,57,58,59]. In ultra-small zSiNTs, e.g., Si(4,0), there are two quasi-flat bands located between E v = 0 0.25 eV. The dispersionless feature is attributed to the lengthened Si-Si bonds in the strongly distorted honeycomb structure (with additional bondings inside). The HVB (LCB) in Si(6,0) is upward (downward) parabolic dispersive along the Γ -X direction crossing E F = 0 between points Γ and K. The 1D band crossing denotes the metallicity. When the radius is increased, the band crossing is replaced by a tiny band gap at the Γ point, as shown in Si(10,0)–Si(21,0) (the band gap size ( E g ) is discussed later and summarized in Table 1 and illustrated in Figure 4). Only very small zSiNTs belong to the 1D metals. This unusual metallic behavior is not dominated by the periodic boundary conditions and π bonding; it mainly comes from the significant s p 3 -orbital hybridization on a highly curved and buckled surface. The σ and π band mixing may also play a significant role in curved silicene nanoribbons. These materials are in contrast to planar systems, such as graphene, where the σ and π orbitals appear orthogonal and cannot exist in a mixed state [60,61]. It is worth noting that a peculiar quasi-flat band exists in all Si(n = even,0) tubes that moves from E v 0.25 eV (Si(4,0)) to E v 1 eV (for Si( n > 10 ,0)). The appearance or absence of the quasi-flat band in Si(n = even/odd,0) tubes are related to the D n h /D n d symmetry [56].

3.3. Band Gap Variation with Tube Radius

The size of the band gap ( E g ) versus the tube mean radius ( R m e a n = R 1 + R 2 2 ) deserves further discussion (Figure 4). The indirect E g increases and decreases as the radius increases for Si(2,2)–Si(4,4). The band gap size variance occurs according to the movement of the HVB-edge state from the Γ to K points. E g increases again when the direct E g is formed in Si(5,5). After that, E g decreases monotonically as the radius increases further. The narrow gaps induced in aSiNTs are attributed to the combined effects of buckling and curvature, which considerably contrasts with metallic aCNTs [62]. The buckling (summarized in Table 1) and curvature decrease as the radius increases, causing a reduction in E g . The zigzag tube zero gaps in Si(4,0)–Si(9,0) correspond to the existence of quasi-flat bands or 1D band crossings. Non-zero E g occurs in Si(10,0) and reaches the first local maximum in Si(14,0). After that, E g presents an oscillatory behavior that includes the global maximum in Si(16,0) and the local minima (nonzero) in Si(15,0), Si(18,0), and Si(21,0), i.e., when n = 3q (q is an integer) [48]. Some publications discussed hexagonal zSiNTs (in analogy to carbon nanotubes) [23,45]. In this condition, n varied from 5 to 12 are metals. Band gap is opened at n = 13 and then presents an oscillatory behavior with the zero gaps occurring at the chiral vector (n = 3q,0). That is, the gear-like zSiNTs with large radii (n = 10–21) are always semiconductors, while the hexagonal zSiNTs have similar electronic properties with carbon nanotubes, which can be metal or semiconductors. The semiconducting features in the gear-like zSiNTs may be attributed to the two inequivalent silicon atoms that have the different distances to the tube axis (Table 1).

3.4. Spatial Charge Density Distribution

The low-lying energy bands are closely related to the type of orbital bonding. Spatial charge distributions, which are sensitive to the surface profile, are used for identifying and examining the multi- and single-orbital hybridizations and the charge transfer [63]. Armchair and zizag SiNTs present similar and distinct charge distributions, as do the orbital hybridizations, demonstrated by the (x,z)-plane projection shown in Figure 5 and Figure 6, respectively. The aSiNTs possess clear σ and π bonds (Figure 5). The (3 s, 3 p x , 3 p y ) orbitals of Si interact with adjacent Si through σ bonds (red rectangles), while the perpendicular 2 p z orbitals show significant π bonding (green rectangles). Both bond types survive in the ultra-small tubes but become non-orthogonal to each other. The σ bond is strengthened with the truncated bond length as the radius increases and the buckling decreases, while the opposite is true for the observable sp 3 -orbital bonds. With a sufficiently large radius, e.g., Si(15,15)–Si(25,25), the σ -bond strength is similar to planar monolayer silicene [64], i.e., the tubular structure does not affect the charge distribution much, as revealed in the bond lengths u 1 and u 2 that are close to 1 (Table 1). The significant s p 3 -orbital bonds in the small tubes result from the misorientations in the (3 p x , 3 p y , 3 p z ) orbitals. The misorientations lead to drastic changes in the low-energy band structures, e.g., the unusual small energy gaps and the nonmonotonous energy dispersions (Figure 2 and Figure 3). To fully comprehend the essential properties of the gear-like tubular structures, non-negligible sp 3 multi-orbital hybridizations might be simultaneously included in the tight-binding model. These four-orbital hybridizations are more pronounced in zigzag than armchair systems (Figure 6) due to the longer Si-Si distance along the azimuthal direction. The charge density is more localized around the Si atoms, especially with a small radius, e.g., Si(4,0). The yellow region between the neighboring Si atoms is dominant in Si(4,0), which corresponds to an sp 3 -orbital bond. The red region (highest density) in the semi-ellipse shape gradually transforms into a triangle as the radius is increased. This change is caused by the decreased buckling that causes significant charge transfer from the sp 3 -orbital bond to the sp 2 -orbital bond.

3.5. The Projected Density of States

The projected density of states (PDOS) provides valuable information about the cooperative and competitive relationships between the chiral angle, the periodic boundary conditions, and the multi-orbital hybridizations, as shown in Figure 7 and Figure 8. Most of the peaks are asymmetric, except for a few symmetric peaks that correspond to quasi-flat bands. A pair of asymmetric peaks in aSiNTs (Figure 7) near the Fermi level is associated with the tiny band gap and low-lying parabolic energy bands (Figure 2). However, their compositions strongly depend on the radius. The HVB edge compositions in the indirect-gap semiconductor Si(2,2) come from the p z (brown curve) and s (red curve) orbitals equally. However, for the direct-gap semiconductor Si(5, 5), the HVB edge compositions come from the p z - and p x -orbitals (green curve). The p x -orbital component gradually decreases as the radius increases. Subsequently, the PDOSs in Si(15,15)–Si(25,25) around the E F = 0 are dominated by the p z -orbital only, similar to monolayer silicene. The more significant contribution of p x compared to p y (blue curve) to the HVB edges in small tubes is ascribed to the strong curvature-induced hybridization along the circumference direction. For large tubes, e.g., Si(10,10)–Si(25,25), the p z -orbital begins at the HVB edge (close to the zero energy), while the p x and p y orbitals start at E v ≈−1.2 eV. The former and latter have a strong hybridization at E v −2.5 eV, where the most substantial peaks exist. Many anti-crossing bands occur at this energy (Figure 2c–f). Though the compositions of the HVB edges in zSiNTs with small radii are similar to aSiNTs, the PDOS has a finite value at E F = 0, e.g., in Si(4,0)–Si(6,0), reflecting metallic behavior due to the quasi-flat bands (Figure 3a) and the 1D band crossings (Figure 3b). The similarities in PDOS features between aSiNTs and zSiNTs include a strong dominance of the π -orbital around E F = 0 in large raddi conditions and the initial energy of E v −1.2 eV for the p x and p y orbitals. On the other hand, large aSiNTs, e.g., Si(15,15), Si(20,20), and Si(25,25) show two strong symmetric peaks at E c −1 and −3.2 eV. The former corresponds to the highly degenerate states at the X point, while the latter corresponds to the quasi-flat band induced by the s p 3 -orbital hybridization. Only the latter one exists in zSiNTs. Generally, independent 3 p z orbitals are not observable in the small armchair or zigzag SiNTs, but the opposite is true in large ones. In the large SiNTs, the π -orbital is in the range of E v = 0 1 eV, similar to monolayer silicene.

4. Conclusions

In summary, the geometric structure and electronic properties of SiNTs with different chirality and diameters were computed and analyzed in detail using first-principles calculations. The results demonstrate that strong multi-orbital hybridizations induce semiconducting and metallic properties in the small armchair and zigzag tubes, respectively. There is an indirect-to-direct band gap transition in aSiNTs when the tube diameter increases, while a metal-semiconductor transition occurs in zSiNTs. The transitions are attributed to the decreasing curvature and buckling, which transfers the complex orbital hybridization to single-orbital bonding. The orbital hybridization creates various electronic properties, including unusual small energy gaps, nonmonotonous energy dispersions, the existence of quasi-flat bands, and frequent anti-crossings between σ and π bands. The calculated PDOS reveals the equal and nonequal contributions of p x and p y to HVB edges for sp/sp 2 /sp 3 bonding and zero and nonzero PDOS at E F = 0 depending on the chirality and diameters. In addition, a strong symmetric peak located at the middle energy level ( E v = −3.2 eV) demonstrates the significant sp 3 bond. These findings shed light on the potential applications of SiNTs in nanoelectronic devices.

Author Contributions

Conceptualization, H.-Y.L., M.-F.L. and J.-Y.W.; methodology, H.-Y.L., M.-F.L. and J.-Y.W.; software, H.-Y.L., M.-F.L. and J.-Y.W.; validation, H.-Y.L. and J.-Y.W.; formal analysis, H.-Y.L. and J.-Y.W.; investigation, H.-Y.L. and J.-Y.W.; resources, H.-Y.L. and J.-Y.W.; data curation, H.-Y.L. and J.-Y.W.; writing—original draft preparation, H.-Y.L. and J.-Y.W.; writing—review and editing, H.-Y.L. and J.-Y.W.; visualization, H.-Y.L. and J.-Y.W.; supervision, M.-F.L. and J.-Y.W.; funding acquisition, M.-F.L. and J.-Y.W. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Science Council of Taiwan (Grant No. 107-2112- M-194-010-MY3) (Grant No. 109-2112-M-992-001-).

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Iijima, S. Helical microtubules of graphitic carbon. Nature 1991, 354, 56–58. [Google Scholar] [CrossRef]
  2. Journet, C.; Maser, W.K.; Bernier, P.; Loiseau, A.; de la Chapelle, M.L.; Lefrant, S.; Deniard, P.; Lee, R.; Fischer, J.E. Large-scale production of single-walled carbon nanotubes by the electric-arc technique. Nature 1997, 388, 756–758. [Google Scholar] [CrossRef]
  3. Thess, A.; Lee, R.; Nikolaev, P.; Dai, H.; Petit, P.; Robert, J.; Xu, C.; Lee, Y.H.; Kim, S.G.; Rinzler, A.G.; et al. Crystalline Ropes of Metallic Carbon Nanotubes. Science 1996, 273, 483–487. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Yudasaka, M.; Komatsu, T.; Ichihashi, T.; Iijima, S. Single-wall carbon nanotube formation by laser ablation using double-targets of carbon and metal. Chem. Phys. Lett. 1997, 278, 102–106. [Google Scholar] [CrossRef]
  5. Ebbesen, T.W.; Ajayan, P.M. Large-scale synthesis of carbon nanotubes. Nature 1992, 358, 220–222. [Google Scholar] [CrossRef]
  6. Ebbesen, T.W.; Hiura, H.; Fujita, J.; Ochiai, Y.; Matsui, S.; Tanigaki, K. Patterns in the bulk growth of carbon nanotubes. Chem. Phys. Lett. 1993, 209, 83–90. [Google Scholar] [CrossRef]
  7. Seraphin, S.; Zhou, D.; Jiao, J.; Withers, J.C.; Loutfy, R. Effect of processing conditions on the morphology and yield of carbon nanotubes. Carbon 1993, 31, 685–689. [Google Scholar] [CrossRef]
  8. Iijima, S.; Ichihashi, T. Single-shell carbon nanotubes of 1-nm diameter. Nature 1993, 363, 603–605. [Google Scholar] [CrossRef]
  9. Lan, A.; Iqbal, Z.; Aitouchen, A.; Libera, M.; Grebel, H. Growth of single-wall carbon nanotubes within an ordered array of nanosize silica spheres. Appl. Phys. Lett. 2002, 81, 433–435. [Google Scholar] [CrossRef]
  10. Altan, H.; Huang, F.; Federici, J.F.; Lan, A.; Grebel, H. Optical and electronic characteristics of single walled carbon nanotubes and silicon nanoclusters by tetrahertz spectroscopy. J. Appl. Phys. 2004, 96, 6685–6689. [Google Scholar] [CrossRef] [Green Version]
  11. Schnorr, J.M.; Swager, T.M. A Emerging Applications of Carbon Nanotubes. Chem. Mater. 2011, 23, 646–657. [Google Scholar] [CrossRef] [Green Version]
  12. Zhang, D.W.; Zhao, Y.B.; Goodenough, J.B.; Lu, Y.H.; Chen, C.H.; Wang, L.; Liu, J.W. A Exfoliation from carbon nanotubes versus tube size on lithium insertion. Electrochem. Commun. 2011, 13, 125–128. [Google Scholar] [CrossRef]
  13. Zhang, H.; Cao, G.P.; Yang, Y.S. A Carbon nanotube arrays and their composites for electrochemical capacitors and lithium-ion batteries. Energy Environ. Sci. 2009, 2, 932–943. [Google Scholar] [CrossRef]
  14. Wang, Y.; Wu, M.; Jiao, Z.; Lee, J.Y. Sn@CNT and Sn@C@CNT nanostructures for superior reversible lithium ion storage. Chem. Mater. 2009, 21, 3210–3215. [Google Scholar] [CrossRef]
  15. Han, Y.; Wang, Y.J.; Wang, Y.P.; Jiao, L.F.; Yuan, H.T.; Liu, S.X. A Hydrothermal synthesis and electrochemical properties of cobalt-carbon nanotubes nanocomposite. Electrochim. Acta. 2011, 56, 3258–3263. [Google Scholar] [CrossRef]
  16. Welna, D.T.; Qu, L.T.; Taylor, B.E.; Dai, L.M.; Durstock, M.F. Vertically aligned carbon nanotube electrodes for lithium-ion batteries. J. Power Sources 2011, 196, 1455–1460. [Google Scholar] [CrossRef]
  17. Wang, Y.; Lee, J.Y. One-step, confined growth of bimetallic tin-antimony nanorods in carbon nanotubes grown in situ for reversible Li+ ion storage. Angew. Chem. Int. Ed. 2006, 45, 7039–7042. [Google Scholar] [CrossRef]
  18. Benvidi, A.; Kakoolaki, P.; Zare, H.R.; Vafazadeh, R. Electrocatalytic oxidation of hydrazine at a Co(II) complex multi-wall carbon nanotube modified carbon paste electrode. Electrochim. Acta. 2011, 56, 2045–2050. [Google Scholar] [CrossRef]
  19. Su, Q.; Liang, Y.; Feng, X.; Mullen, K. Towards free-standing graphene/carbon nanotube composite films via acetylene-assisted thermolysis of organocobalt functionalized graphene sheets. Chem. Commun. 2010, 46, 8279–8281. [Google Scholar] [CrossRef]
  20. Liu, S.; Zhu, J.; Mastai, Y.; Felner, I.; Gedanken, A. Preparation and characteristics of carbon nanotubes filled with cobalt. Chem. Mater. 2000, 12, 2205–2211. [Google Scholar] [CrossRef]
  21. Duran-Valdeiglesias, E.; Zhang, W.; Alonso-Ramos, C.; Serna, S.; Le Roux, X.; Maris-Morini, D.; Caselli, N.; Biccari, F.; Gurioli, M.; Filoramo, A.; et al. Tailoring carbon nanotubes optical properties through chirality-wise silicon ring resonators. Sci. Rep. 2018, 8, 11252. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Fagan, S.B.; Baierle, R.J.; Mota, R.; da Silva, A.J.R.; Fazzio, A. Ab initio calculations for a hypothetical material: Silicon nanotubes. Phys. Rev. B 2000, 61, 9994–9996. [Google Scholar] [CrossRef]
  23. Chegel, R.; Behzad, S. Bandstructure modulation for Si-h and Si-g nanotubes in a transverse electric field: Tight binding approach. Superlattices Microstruct. 2013, 63, 79–90. [Google Scholar] [CrossRef]
  24. Kiricsi, I.; Fudala, A.; Konya, Z.; Hernadi, K.; Lentz, P.; Nagy, J.B. The advantages of ozone treatment in the preparation of tubular silica structures. Appl. Catal. A Gen. 2000, 203, L1. [Google Scholar] [CrossRef]
  25. Schmidt, O.G.; Eberl, K. Nanotechnology—Thin solid films roll up into nanotubes. Nature 2001, 410, 168. [Google Scholar] [CrossRef]
  26. Sha, J.; Niu, J.; Ma, X.; Xu, J.; Zhang, X.; Yang, Q.; Yang, D. Silicon Nanotubes. Adv. Mater. 2002, 14, 1219–1221. [Google Scholar] [CrossRef]
  27. Jeong, S.Y.; Kim, J.Y.; Yang, H.D.; Yoon, B.N.; Choi, S.H.; Kang, H.K.; Yang, C.W.; Lee, Y.H. Synthesis of Silicon Nanotubes on Porous Alumina Using Molecular Beam Epitaxy. Adv. Mater. 2003, 15, 1172–1176. [Google Scholar] [CrossRef]
  28. Chen, Y.W.; Tang, Y.H.; Pei, L.Z.; Guo, C. Self-Assembled Silicon Nanotubes Grown from Silicon Monoxide. Adv. Mater. 2005, 17, 564–567. [Google Scholar] [CrossRef]
  29. Tang, Y.H.; Pei, L.Z.; Chen, Y.M.; Guo, C. Self-assembled silicon nanotubes under supercritically hydrothermal conditions. Phys. Rev. Lett. 2005, 95, 116102. [Google Scholar] [CrossRef]
  30. De Crescenzi, M.; Castrucci, P.; Scarcelli, M.; Diociauti, M.; Chaudhari, P.S.; Balasubramanian, C.; Bhave, T.M.; Bhoraskar, S.V. Experimental imaging of silicon nanotubes. Appl. Phys. Lett. 2005, 86, 231901. [Google Scholar] [CrossRef]
  31. Xie, M.; Wang, J.S.; Fan, Z.Y.; Lu, J.G.; Yap, Y.K. Growth of p-type Si nanotubes by catalytic plasma treatments. Nanotechnology 2008, 19, 365609. [Google Scholar] [CrossRef] [Green Version]
  32. Taghinejad, M.; Taghinejad, H.; Abdolahad, M.; Mohajerzadeh, S. A Nickel-Gold Bilayer Catalyst Engineering Technique for Self-Assembled Growth of Highly Ordered Silicon Nanotubes (SiNT). Nano Lett. 2013, 13, 889–897. [Google Scholar] [CrossRef] [PubMed]
  33. Morata, A.; Pacios, M.; Gadea, G.; Flox, C.; Cadavid, D.; Cabot, A.; Tarancon, A. Large-area and adaptable electrospun silicon-based thermoelectric nanomaterials with high energy conversion efficiencies. Nat. Commun. 2018, 9, 4759. [Google Scholar] [CrossRef] [PubMed]
  34. Castrucci, P.; Scarselli, M.; De Crescenzi, M.; Diociaiuti, M.; Chaudhari, P.S.; Balasubramanian, C.; Bhave, T.M.; Bhoraskar, S.V. Silicon nanotubes: Synthesis and characterization. Thin Solid Films 2006, 508, 226–230. [Google Scholar] [CrossRef]
  35. Ishai, M.B.; Patolsk, F. Shape- and Dimension-Controlled Single-Crystalline Silicon and SiGe Nanotubes: Toward Nanofluidic FET Devices. J. Am. Chem. Soc. 2009, 131, 3679–3689. [Google Scholar] [CrossRef] [PubMed]
  36. Yoo, H.; Park, E.; Bae, J.; Lee, J.; Chung, D.J.; Jo, Y.N.; Park, M.S.; Kim, J.H.; Dou, S.X.; Kim, Y.J.; et al. Si Nanocrystal-Embedded SiOx nanofoils: Two-Dimensional Nanotechnology-Enabled High Performance Li Storage Materials. Sci. Rep. 2018, 8, 6904. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Kasavajjula, U.; Wang, C.S.; Appleby, A.J. Nano- and bulk-silicon-based insertion anodes for lithium-ion secondary cells. J. Power Sources 2007, 163, 1003–1039. [Google Scholar] [CrossRef]
  38. Chan, C.K.; Peng, H.; Liu, G.; McIlwrath, K.; Zhang, X.F.; Huggins, R.A.; Cui, Y. High-performance lithium battery anodes using silicon nanowires. Nat. Nanotechnol. 2008, 3, 31–35. [Google Scholar] [CrossRef]
  39. Zhu, J.; Yu, Z.; Burkhard, G.F.; Hsu, C.M.; Connor, S.T.; Xu, Y.; Wang, Q.; McGehee, M.; Fan, S.; Cui, Y. Optical Absorption Enhancement in Amorphous Silicon Nanowire and Nanocone Arrays. Nano lett. 2009, 9, 279–282. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Ahmadi, N.; Shokri, A.A.; Elahi, S.M. Optical Transition of Zigzag Silicon Nanotubes Under Intrinsic Curvature Effect. SILICON 2016, 8, 217–224. [Google Scholar] [CrossRef]
  41. Behzad, S.; Chegel, R. Magnetic Field-Induced Splitting of Optical Spectra in Silicon Nanotubes: Tight Binding Calculations. SILICON 2016, 8, 43–55. [Google Scholar] [CrossRef]
  42. Chegel, R.; Behzad, S. Electronic Properties of SiNTs Under External Electric and Magnetic Fields Using the Tight-Binding Method. J. Electron. Mater. 2014, 43, 329–340. [Google Scholar] [CrossRef]
  43. Zhang, R.Q.; Lee, H.L.; Li, W.K.; Teo, B.K. Investigation of Possible Structures of Silicon Nanotubes via Density-Functional Tight-Binding Molecular Dynamics Simulations and ab Initio Calculations. J. Phys. Chem. B 2005, 109, 8605–8612. [Google Scholar] [CrossRef] [PubMed]
  44. Seifert, G.; Koehler, T.; Urbassek, H.M.; Hernandez, E.; Frauenheim, T. Tubular structures of silicon. Phys. Rev. B 2001, 63, 193409. [Google Scholar] [CrossRef] [Green Version]
  45. Yang, X.; Ni, J. Electronic properties of single-walled silicon nanotubes compared to carbon nanotubes. Phys. Rev. B 2005, 72, 195426. [Google Scholar] [CrossRef]
  46. Barnard, A.S.; Russo, S.P. Structure and Energetics of Single-Walled Armchair and Zigzag Silicon Nanotubes. J. Phys. Chem. B 2003, 107, 7577–7581. [Google Scholar] [CrossRef]
  47. Matheus, P.L. Double-walled silicon nanotubes: An ab initio investigation. Nanotechnology 2018, 29, 075703. [Google Scholar]
  48. Wang, C.; Fu, X.; Guo, Y.; Guo, Z.; Xia, C.; Jia, Y. Band gap scaling laws in group IV nanotubes. Nanotechnology 2017, 28, 115202. [Google Scholar] [CrossRef] [Green Version]
  49. Kresse, G.; Furthmueller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 1996, 54, 11169–11186. [Google Scholar] [CrossRef]
  50. Kresse, G.; Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. Rev. B 1999, 59, 1758–1775. [Google Scholar] [CrossRef]
  51. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 1996, 77, 3865–3868. [Google Scholar] [CrossRef] [Green Version]
  52. Guelseren, O.; Yildirim, T.; Ciraci, S. Systematic ab initio study of curvature effects in carbon nanotubes. Phys. Rev. B 2002, 65, 153405. [Google Scholar] [CrossRef] [Green Version]
  53. Saito, R.; Fujita, M.; Dresselhaus, G.; Dresselhaus, U.M. Electronic structure of chiral graphene tubules. Appl. Phys. Lett. 1992, 60, 2204–2206. [Google Scholar] [CrossRef]
  54. Hamada, N.; Sawada, S.; Oshiyama, A. New One-Dimensional Conductors-Graphitic Microtubules. Phys. Rev. Lett. 1992, 68, 1579–1581. [Google Scholar] [CrossRef] [PubMed]
  55. Mintmire, J.W.; White, C.T. Universal density of states for carbon nanotubes. Phys. Rev. Lett. 1998, 81, 2506–2509. [Google Scholar] [CrossRef] [Green Version]
  56. Saito, R.; Dresselhaus, G.; Dresselhaus, M.S. Physical Properties of Carbon Nanotubes; Imperial College Press: London, UK, 1998; p. 272. [Google Scholar]
  57. Ouyang, M.; Huang, J.L.; Cheung, C.L.; Liebe, C.M. Energy gaps in “metallic” single-walled carbon nanotubes. Science 2001, 292, 702–705. [Google Scholar] [CrossRef] [Green Version]
  58. Kane, C.L.; Mele, E.J. Size, shape and low energy electronic structure of carbon nanotubes. Phys. Rev. Lett. 1997, 78, 1932–1935. [Google Scholar] [CrossRef] [Green Version]
  59. Blase, X.; Benedict, L.X.; Shirley, E.L.; Louie, S.G. Hybridization effects and metallicity in small radius carbon nanotubes. Phys. Rev. Lett. 1994, 72, 1878–1881. [Google Scholar] [CrossRef] [PubMed]
  60. Lin, S.Y.; Lin, Y.T.; Thuy Tran, N.T.; Su, W.P.; Lin, M.F. Feature-rich electronic properties of aluminum-adsorbed graphenes. Carbon 2017, 120, 209–218. [Google Scholar] [CrossRef]
  61. Lin, S.Y.; Chang, S.L.; Thuy Tran, N.T.; Yang, P.H.; Lin, M.F. H-Si bonding-induced unusual electronic properties of silicene: A method to identify hydrogen concentration. Phys. Chem. Chem. Phys. 2015, 17, 26443–26450. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Umeno, Y.; Kitamura, T.; Kushima, A. Metallic-semiconducting transition of single-walled carbon nanotubes under high axial strain. Comput. Mater. Sci. 2004, 31, 33–41. [Google Scholar] [CrossRef]
  63. Jacobsen, H. Chemical Bonding in View of Electron Charge Density and Kinetic Energy Density Descriptors. J. Comput. Chem. 2009, 30, 1093–1102. [Google Scholar] [CrossRef] [PubMed]
  64. Liu, H.; Lin, S.Y.; Wu, J. Stacking-configuration-enriched essential properties of bilayer graphenes and silicenes. J. Chem. Phys. 2020, 153, 154707. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Cross-section and side views for (a) armchair (9,9) and (b) zigzag (9,0) SiNTs.
Figure 1. Cross-section and side views for (a) armchair (9,9) and (b) zigzag (9,0) SiNTs.
Nanomaterials 11 02475 g001
Figure 2. Band structures of aSiNTs with different chiral vectors: (a) (2,2), (b) (5,5), (c) (10,10), (d) (15,15), (e) (20,20), and (f) (25,25).
Figure 2. Band structures of aSiNTs with different chiral vectors: (a) (2,2), (b) (5,5), (c) (10,10), (d) (15,15), (e) (20,20), and (f) (25,25).
Nanomaterials 11 02475 g002
Figure 3. Band structures of zSiNTs with different chiral vectors: (a) (4,0), (b) (6,0), (c) (10,0), (d) (12,0), (e) (14,0), (f) (16,0), (g) (18,0), (h) (20,0), and (i) (21,0).
Figure 3. Band structures of zSiNTs with different chiral vectors: (a) (4,0), (b) (6,0), (c) (10,0), (d) (12,0), (e) (14,0), (f) (16,0), (g) (18,0), (h) (20,0), and (i) (21,0).
Nanomaterials 11 02475 g003
Figure 4. The dependence of band gap E g on the radius for aSiNTs (black curve) and zSiNTs (red curve).
Figure 4. The dependence of band gap E g on the radius for aSiNTs (black curve) and zSiNTs (red curve).
Nanomaterials 11 02475 g004
Figure 5. The spatial charge densities on the (x,z)-plane for aSiNTs with different chiral vectors: (a) (2,2), (b) (5,5), (c) (10,10), (d) (15,15), (e) (20,20), and (f) (25,25).
Figure 5. The spatial charge densities on the (x,z)-plane for aSiNTs with different chiral vectors: (a) (2,2), (b) (5,5), (c) (10,10), (d) (15,15), (e) (20,20), and (f) (25,25).
Nanomaterials 11 02475 g005
Figure 6. The spatial charge densities on the (x,z)-plane for zSiNTs with different chiral vectors: (a) (4,0), (b) (6,0), (c) (10,0), (d) (12,0), (e) (14,0), (f) (16,0), (g) (18,0), (h) (20,0), and (i) (21,0).
Figure 6. The spatial charge densities on the (x,z)-plane for zSiNTs with different chiral vectors: (a) (4,0), (b) (6,0), (c) (10,0), (d) (12,0), (e) (14,0), (f) (16,0), (g) (18,0), (h) (20,0), and (i) (21,0).
Nanomaterials 11 02475 g006
Figure 7. The orbital-projected densities of states for aSiNTs with different chiral vectors: (a) (2,2), (b) (5,5), (c) (10,10), (d) (15,15), (e) (20,20), and (f) (25,25).
Figure 7. The orbital-projected densities of states for aSiNTs with different chiral vectors: (a) (2,2), (b) (5,5), (c) (10,10), (d) (15,15), (e) (20,20), and (f) (25,25).
Nanomaterials 11 02475 g007
Figure 8. The orbital-projected densities of states for zSiNTs with different chiral vectors: (a) (4,0), (b) (6,0), (c) (10,0), (d) (12,0), (e) (14,0), (f) (16,0), (g) (18,0), (h) (20,0), and (i) (21,0).
Figure 8. The orbital-projected densities of states for zSiNTs with different chiral vectors: (a) (4,0), (b) (6,0), (c) (10,0), (d) (12,0), (e) (14,0), (f) (16,0), (g) (18,0), (h) (20,0), and (i) (21,0).
Nanomaterials 11 02475 g008
Table 1. Geometric properties of aSiNTs and zSiNTs: lattice constant, ground state energy E 0 , Si-Si bond lengths u 1 (“perpendicular” to the axis) and u 2 (“along” the axis), radii R 1 and R 2 , buckling distance Δ z, and band gaps E g .
Table 1. Geometric properties of aSiNTs and zSiNTs: lattice constant, ground state energy E 0 , Si-Si bond lengths u 1 (“perpendicular” to the axis) and u 2 (“along” the axis), radii R 1 and R 2 , buckling distance Δ z, and band gaps E g .
aSiNTsLattice con. (Å) E 0 (eV/atom) u 1 / u 0 u 2 / u 0 R 1 (Å) R 2 (Å) Δ z (Å) E g (eV)
Si(2,2)3.54−4.7461.0031.0431.7312.6880.9570.350
Si(3,3)3.722−4.7861.0011.0152.5003.1810.6810.785
Si(4,4)3.751−4.8121.0001.0093.9674.5690.6020.360
Si(5,5)3.764−4.8280.9991.0054.7965.3250.5290.485
Si(6,6)3.766−4.8370.9991.0036.1186.6480.5300.382
Si(7,7)3.772−4.8430.9981.0027.0017.5010.5000.316
Si(9,9)3.775−4.8500.9971.0049.1499.6490.5000.237
Si(10,10)3.779−4.8510.9971.00410.36210.8510.4890.208
Si(12,12)3.780−4.8540.9991.00012.48812.9700.4810.168
Si(15,15)3.777−4.8560.9970.99915.59515.9090.3140.131
Si(20,20)3.782−4.8580.9960.99920.94321.1510.2080.099
Si(25,25)3.782−4.8590.9950.99826.42826.6170.1890.078
zSiNTsLattice con. (Å) E 0 (eV/atom) u 1 / u 0 u 2 / u 0 R 1 (Å) R 2 (Å) Δ z (Å) E g (eV)
Si(4,0)3.806−4.7031.1321.0331.8262.8811.0550
Si(5,0)3.814−4.6941.0221.0252.2823.1630.8810
Si(6,0)3.757−4.7081.0081.0113.4074.0750.6680
Si(7,0)3.760−4.7201.0061.0083.9294.5050.5760
Si(8,0)3.770−4.7311.0041.0074.6615.2220.5610
Si(9,0)3.773−4.7401.0041.0065.1995.7330.5340
Si(10,0)3.773−4.7471.0041.0055.9066.4340.5280.186
Si(11,0)3.773−4.7511.0031.0046.4636.9730.5100.190
Si(12,0)3.776−4.7551.0031.0047.1457.6530.5080.202
Si(13,0)3.779−4.7581.0031.0037.7118.2060.4950.320
Si(14,0)3.776−4.7601.0021.0028.3918.8830.4920.325
Si(15,0)3.780−4.7621.0011.0028.9459.4280.4830.248
Si(16,0)3.777−4.7641.0011.0019.62110.1040.4830.237
Si(17,0)3.779−4.7651.0011.00110.19110.6690.4780.299
Si(18,0)3.782−4.7661.0011.00210.84511.3230.4780.202
Si(19,0)3.781−4.7671.0011.00111.42611.8990.4730.282
Si(20,0)3.780−4.7681.0011.00112.08812.5640.4760.254
Si(21,0)3.781−4.7691.0011.00112.65913.1290.4700.170
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, H.-Y.; Lin, M.-F.; Wu, J.-Y. Essential Electronic Properties of Silicon Nanotubes. Nanomaterials 2021, 11, 2475. https://doi.org/10.3390/nano11102475

AMA Style

Liu H-Y, Lin M-F, Wu J-Y. Essential Electronic Properties of Silicon Nanotubes. Nanomaterials. 2021; 11(10):2475. https://doi.org/10.3390/nano11102475

Chicago/Turabian Style

Liu, Hsin-Yi, Ming-Fa Lin, and Jhao-Ying Wu. 2021. "Essential Electronic Properties of Silicon Nanotubes" Nanomaterials 11, no. 10: 2475. https://doi.org/10.3390/nano11102475

APA Style

Liu, H. -Y., Lin, M. -F., & Wu, J. -Y. (2021). Essential Electronic Properties of Silicon Nanotubes. Nanomaterials, 11(10), 2475. https://doi.org/10.3390/nano11102475

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop