HOT Graphene and HOT Graphene Nanotubes
Total Page:16
File Type:pdf, Size:1020Kb
Xu et al. Nanoscale Research Letters (2020) 15:56 https://doi.org/10.1186/s11671-020-3279-1 NANO REVIEW Open Access HOT Graphene and HOT Graphene Nanotubes: New Low Dimensional Semimetals and Semiconductors Lin-Han Xu1* , Shun-Qing Wu1, Zhi-Quan Huang2, Feng Zhang3, Feng-Chuan Chuang2, Zi-Zhong Zhu1* and Kai-Ming Ho3,4 Abstract We report a new graphene allotrope named HOT graphene containing carbon hexagons, octagons, and tetragons. A corresponding series of nanotubes are also constructed by rolling up the HOT graphene sheet. Ab initio calculations are performed on geometric and electronic structures of the HOT graphene and the HOT graphene nanotubes. Dirac cone and high Fermi velocity are achieved in a non-hexagonal structure of HOT graphene, implying that the honeycomb structure is not an indispensable condition for Dirac fermions to exist. HOT graphene nanotubes show distinctive electronic structures depending on their topology. The (0,1) n (n ≥ 3) HOT graphene nanotubes reveal the characteristics of semimetals, while the other set of nanotubes (1,0) n shows continuously adjustable band gaps (0~ 0.51 eV) with tube size. A competition between the curvature effect and the zone-folding approximation determines the band gaps of the (1,0) n nanotubes. Novel conversion between semimetallicity and semiconductivity arises in ultra-small tubes (radius < 4 Å, i.e., n < 3). Keywords: HOT graphene, Nanotubes, Electronic structures, First-principles calculations Introduction in sp2-bonded carbon systems are dependent on the Because of its bonding flexibility, carbon-based systems honeycomb structure. In the lower dimension, the carbon show an unlimited number of different structures with nanotube is a honeycomb structure rolled into a hollow an equally large variety of physical properties. These cylinder with nanometric diameter and μmlength[7–10]. physical properties are, in great part, the result of the di- As there is an infinite number of ways of rolling a sheet mensionality of these structures [1]. Graphene is a single into a cylinder, the large variety of possible helical geom- two-dimensional layer of carbon atoms bound in a hex- etries, defining the tube chirality, provides a family of agonal lattice structure [2] revealing a number of unique nanotubes with different diameters and microscopic struc- properties, such as massless carriers, high Fermi velocity tures [11–13]. The electronic and transport properties are [3], and Dirac cones [4, 5], which are characteristic of certainly among the most significant physical properties of two-dimensional Dirac fermions. The honeycomb lattice carbon nanotubes, and crucially depend on the diameter consisting of two equivalent carbon sublattices plays a and chirality [14–18]. Graphene nanotubes can be either crucial role in forming such intriguing properties [2]. semimetallic [14]orsemiconducting[19–21], with a band Enyashin and Ivanovskii [6] constructed 12 artificial 2D gap varying from zero to a few tenths of an eV, depending carbon networks but found no structures other than the on their diameter and chirality [10, 14, 16]. Furthermore, graphyne allotrope exhibit the graphene-like electronic the band gap of semiconducting tubes can be shown to be behavior. It seems to imply that the Dirac-like fermions simply related to the tube diameter. The semimetallic nanotubes also maintain the unique properties from gra- phene, such as massless carriers, high Fermi velocity [22], * Correspondence: [email protected]; [email protected] 1Department of Physics, Collaborative Innovation Center for Optoelectronic and Dirac cones [23]. Such remarkable results can be ob- Semiconductors and Efficient Devices, Jiujiang Research Institute, Xiamen tained from a variety of considerations, starting from the University, Xiamen, China so-called band-folding approach, based on knowledge of Full list of author information is available at the end of the article © The Author(s). 2020 Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Xu et al. Nanoscale Research Letters (2020) 15:56 Page 2 of 11 the electronic properties of the graphene sheet, to the dir- Method of Calculation ect study of nanotubes using semiempirical tight-binding The present calculations on HOT graphene and HOT approaches [14, 16, 18, 23]. Comparing with more sophis- graphene nanotubes were performed by using a first- ticated ab initio calculations and available experimental principles method based on the density-functional results, finer considerations, such as curvature effects, kF theory (DFT) with the generalized gradient approxi- shifting [24, 25], σ-π hybridization [26] are introduced. mation (GGA) in the form of Perdew-Burke-Ernzerh Graphene and graphene-like materials [6] are considered (PBE) exchange-correlation functional [43], as imple- a revolutionary material for future generation of high- mented in the Vienna Ab initio Simulation Package speed electronic, radio frequency logic devices [27, 28], (VASP) [44, 45]. The wave functions were expanded thermally and electrically conductive reinforced compos- in plane waves up to a cutoff of plane wave kinetic ites [29, 30], catalyst [31], sensors [32–35], transparent energy of 520 eV. The Brillouin zone (BZ) integrals electrodes [27, 36], etc. basing on the unusual properties were performed by using a Monkhorst-Pack [46]sam- all above. Over the past few decades, carbon nanotubes pling scheme with a k-point mesh resolution of 2π × − have also shown great potential in logic circuits, gas stor- 0.03 Å 1. The unit cell basis vectors (representing age, catalysis, and energy storage because of their extraor- unit cell shape and size) and atomic coordinates were dinary electronic, mechanical, and structural properties fully relaxed in each system until the forces on all the [37–39]. Hence, the creation of new carbon allotropes (in- atoms were smaller than 0.01 eV/Å. cluding 2D and 1D) has been the focus of numerous theor- etical and experimental explorations because of their Results and Discussion fundamental scientific and technological importance [40]. Geometric and Electronic Structures of HOT Graphene However, completely clarifying the structures of these excit- The geometric structure of HOT graphene (Fig. 1a) ing carbon phases through current experimental technolo- shows a more complicated bonding situation than gra- gies is usually unrealistic due to their limited quantity, as phene. The variety of carbon polygons in HOT graphene well as the mixture of other phases. Theoretical prediction results in various carbon bonding characters. These is necessary and has yielded great success [31–35, 40–42]. polygons in HOT graphene share common edges with In this study, we designed a new allotrope of graphene each other, and the bonds can be distinguished by the that has two-dimensional Dirac fermions without an ex- two polygons they belong to. Therefore, in our research, clusively hexagonal structure. The new allotrope was con- they are named as 6–8 bonds, 4–8 bonds, 4–6 bonds, 6– structed with interlaced carbon hexagons, octagons, and 6 bonds, and 8–8 bonds. The 4–8 bonds and 6–8 bonds tetragons, and was named HOT graphene. HOT graphene have two different bond lengths: 1.44 Å and 1.47 Å for nanotubes were also constructed by rolling up HOT gra- 4–8 bonds; 1.41 Å and 1.48 Å for 6–8 bonds. The 4–6 phene sheet along with different directions. The electronic bonds, 6–6 bonds, and 8–8 bonds have unique bond property, curvature effect, kF shifting effect, etc. of HOT lengths of 1.44 Å, 1.46 Å, and 1.34 Å, relatively. Figure 1b graphene and nanotubes were calculated using ab initio shows the band structure and density of states (DOS) of calculations based on density function theory (DFT). HOT graphene with the corresponding BZ depicted in Fig. 1 (a) Geometry of HOT graphene; (b) Band structures and DOS of the HOT graphene; (c) the corresponding BZ of HOT graphene Xu et al. Nanoscale Research Letters (2020) 15:56 Page 3 of 11 Fig. 1c. The crossing point of energy bands at the Fermi calculation. Different rolling directions between (1,0) n and (0, level indicates semimetallicity of HOT graphene, which 1) n tubes result in the differences in geometry and bonding is confirmed by the vanishing DOS at the Fermi level. situation. Two tubes, (1,0)6 and (0,1)4, are depicted in Fig. 3d The Dirac point is located at (0, 0.0585, 0) adjacent to Γ. and e to describe the geometric differences between the two The 3D band structure (Fig. 2) presents the band sur- rolling directions. The arrangement of polygons along (1,0) faces near the Fermi level, where one can see the Dirac direction can be divided into two patterns: C4–C6–C8 (orange) cones formed by upper and lower conical surfaces meet- and C8–C6–C4 (blue) which is exactly opposite to each other. ing at two Dirac points exactly at the Fermi surface. The These two opposite patterns alternate along the circumferen- corresponding Fermi velocity (vF) of the Dirac fermions, tial direction of the tube. In direction (0,1), polygons arranging evaluated from the gradient of the linear dispersions of along the tube axis also have two patterns: C4–C8 (blue) and 5 the band structures, is 6.27 × 10 m/s, which is a little C6–C6 (orange). Two C4–C8 patterns alternate with one C6– 5 lower than 8.1 × 10 m/s [22] for graphene nanotube and C6 pattern along the circumferential direction. 5 8.6 × 10 m/s [47, 48] for graphene. The high vF implies To reveal the energy cost in rolling up a sheet into high mobility of carriers in the HOT graphene.