1 Introduction
Carbon is an extremely versatile element with multiple chemical bonding possibilities, resulting in the existence of a large number of carbon allotropes [
1]. Graphite and diamond are two classic carbon-based materials. Graphite is a good conductor, and hence, is typically used as a standard electrode [
2-
5]. Diamond is a well-known superhard material, and is widely used as hard tips in experiments [
6-
8]. In 1985, C
60 and other buckyballs were discovered and several novel properties of the buckyballs were explored. After the gold rush for the buckyball, carbon nanotubes have attracted intense research interest since their discovery in 1991. In 2004, single-layer graphene was experimentally obtained for the first time, and has gained considerable attention in the carbon family [
9]. Till date, the discovery and study of carbon allotropes have led to several Nobel Prizes, including the discovery of the buckyball C
60 in 1996 and experiments with graphene in 2010 [
10,
9]. There are other proposed structures such as the recently erethynylated dehydroannulenes, expanded radialenes, radiaannulenes [
11], sparse fullerene structures [
12], conjugated zinc porphyrin nanoball [
13], belt-shaped carbon structures [
14] and others [
15-
17].
Carbon allotropes have wide applications in several industrial fields. Graphite, with a reasonably high Li ion capacity and being strong enough to withstand mechanical strains during the insertion process, has been used as an efficient anode material in Li ion batteries produced by SONY [
18-
20] Diamond, with ultrastrong C-C sp3 bonds, is extremely robust and can be used as a hard tip for various indentation experiments [
21-
24]. Carbon fiber, with a low density and certain advanced mechanical properties, is an important high-quality composite used in the aircraft field [
25-
28].
Due to the industrial and academic significance of carbon allotropes, there has been increasing interest in recent years in discovering more two-dimensional (2D) carbon allotropes. These newly discovered carbon allotropes have certain different characteristic properties that have been discussed separately in the respective original studies. However, a comparative survey of these new carbon allotropes discovered in the last few years is lacking, and is the focus of the present review. A large number of review articles have been published on the well-studied carbon allotropes, such as buckyballs [
10], graphene [
29-
40], diamond [
41-
43], carbon nanotubes [
44-
48]. Hence, we are not going to discuss these well-known carbon allotropes. Instead, we will primarily focus on the 2D and three-dimensional (3D) carbon allotropes explored in the past decade.
In this review, we have presented a comparative survey of the physical and mechanical properties of a variety of carbon allotropes, including seven 2D and eight 3D materials. The results of the comparison, summarized in Tables 2 and 3, will be valuable for further studies on other properties of these carbon allotropes. By analyzing the relationship between the structural topology of these carbon allotropes, we propose a systematic route to synthesize these carbon allotropes by manipulating the carbon atoms or C-C bonds in the two mother structures, graphene and diamond. This theoretical study can facilitate the discovery of new carbon allotropes in a more systematic way.
2 2D semiconductor carbon
Semiconductors play a fundamental role in everyday life and are of key importance for the development of the present electric technique. Several semiconduction 2D carbon allotropes with finite electronic band gaps have been proposed, which will be discussed in this section.
2.1 Graphenylene
Graphenylene, also called biphenylene carbon, was first reported by Balaban
et al. [
62] and was successfully synthesized in the experiment [
63,
64]. First-principles simulations proposed to synthesize graphenylene from porous graphene via dehydrogenation [
65].
Graphenylene has a lattice structure of P
6/mmm symmetry as shown in Fig. 1A [
60]. It is a 2D sp
2-hybridized carbon network with cyclohexatriene unit cells containing twelve carbon atoms. There are three inequivalent C-C bond lengths as listed in Table 1.
Graphenylene is energetically less stable than graphene by about 0.62 eV per atom, but more stable than graphdiyne by about 0.14 eV per atom [
66].
The Young’s modulus for graphenylene is 648.80 GPa (217.3 N/m) (assuming an effective thickness of 0.335 nm) and the Poisson’s ratio is 0.259 [
60]. Anisotropic effects were found on mechanical properties for graphenylene nanoribbons [
67]. The Young’s modulus is 68.31 N/m and 79.42 N/m along the armchair and zigzag directions, respectively. The ultimate strength is 429.43 nN (9.59 N/m) and 1140.59 nN (25.47 N/m) for the armchair and zigzag graphenylene nanoribbons, respectively.
Theoretical results for the electronic band gap of the graphenylene are quite diverse. Small band gaps of 0.025 eV and 0.034 eV have been predicted by two recent works [
56,
68], while a band gap of in the range 0.8-1.08 eV were also predicted [
60,
69,
70]. The electronic band gap is tunable by hydrogenation and halogenation [
59]. Graphenylene in tubal forms have also been studied, and the electronic band gap is about 0.7 eV which can be tuned by axially strain [
71].
Graphenylene has been proposed to be efficient gas separators, by taking advantage of its porous configuration. By stretching graphenylene with different strains, it is able to seperate several gases, such as CO
2, N
2, CO, CH
4, or
3He [
56,
58,
61]. The porous nature of the graphenylene enables it to be a promising candidate as anode material for Li ion battery of high mobility and high capacity [
66,
72].
2.2 Penta-graphene
Penta-graphene was proposed theoretically by Zhang
et al. in 2015 [
50], which was predicted to be a metastable structure that can undergo phase transition via strain or temperature [
73,
74]. Penta-graphene has not been synthesized in the experiment, and may not easy to observe this material according to the works by Ewels
et al. [
75] The cohesive energy per atom is
eV [
76]. Compared with the cohesive energies of graphene (
eV/atom), penta-graphene is less stable by about 0.9 eV per atom.
Figure 1B shows the lattice structure of penta-graphene in square unit cell with lattice constant of 3.64 Å. There are six atoms in the unit cell, two C1 carbon atoms and four C2 carbon atoms. Penta-graphene is a closely resembling of the single bond C1-C2 (1.55 Å) and double bond C2-C2 (1.34 Å) [
77,
78]. It is a quasi 2D material with buckling height of 1.2 Å. Penta-graphene consists 5-carbon rings, which can be obtained through graphene by the reduction of some C-C bonds.
Mechanical properties of peta-graphene have attracted intense research interest in recent years. As a direct result of the buckling configuration, penta-graphene exhibits a negative Poisson’s ratio of
[
50]. The Poisson’s ratio can be increased to 0.243 or 0.236 through hydrogenation or fluorination [
79]. The in-plane Young’s modulus for penta-graphene is about 263.8 GPa·nm (263.8 N/m) in Ref. [
50], or 376 N/m in Ref. [
80]. Cranford
et al. predicted the yield strength to be 12.2 N/m at strain 5.1%, while the ultimate strength to be 23.8 N/m at strain 55.8%, and the elastic toughness to be 5.4 J/m
2 [
80]. Sun
et al. obtained similar value of 23.51 N/m for the ultimate strength at strain 0.18 [
81]. Hydrogenation, fluorination and the hydroxyl functional groups can efficiently tune a variety of mechanical properties of penta-graphene, including Young’s modulus, yield strength, stability, and the buckling height [
74,
79,
82,
83].
Penta-graphene has a rather large electronic band gap typically around 3.2 eV [
50,
84-
87]. Even larger values were also predicted for the electronic band gap of penta-graphene, eg., 4.1–4.3 eV [
77] or 4.48 eV [
76] which are close to the value of insulators. The electronic band gap can be manipulated via various approaches, including stacking style [
84], doping or hydrogenation [
79,
88,
89], vacancy or Stone–Wales defects [
89], and layer number [
76]. Penta-graphene in the nanoribbon form was predicted to exhibit some interesting electronic and magnetic properties [
90,
91].
Theoretical works have predicted high thermal conductivity for penta-graphene, which is dominately contributed by the phonon transport as the electron thermal conductivity is neglectable due to large band gap. The thermal conductivity at room temperature was predicted to be 167 W·m
-1·K
-1 [
92], 197.85 W·m
-1·K
-1 [
93], 645 W·m
-1·K
-1 [
94], and 350 W·m
-1·K
-1 [
95]. The thermal conductivity is insensitive to the layer number of penta-graphene [
96], but it can be tuned through hydrogenation [
95] or oxidization [
97].
There are some potential applications for penta-graphene owing to its quasi-2D and slightly buckling configuration and the large electronic band gap. Theoretical studies have suggested penta-graphene to be applied as a low-cost catalyst for CO oxidation [
98], hydrogen storage [
88], and anode material for Li ion battery [
99].
2.3 Twin graphene
Twin graphene was proposed theoretically in 2017 through first-principles calculations and MD simulations [
51]. There are 18 atoms in the unit cell with two inequivalent carbon atoms C1 and C2 as shown in Fig. 1C, leading to three bond lengths
Å,
Å, and
Å. The lattice structure of the twin graphene can be obtained by substituting the hexagonal carbon ring in the graphyne by a pair of hexagonal carbon rings (see Figs. 1C and D).
From an energetic point of view, a larger cohesive energy means better structural stability. The cohesive energy of twin graphene is smaller than that of graphene by about 0.8–0.9 eV based on different approaches [
51].
Twin graphene has an electronic band gap about 1.0 eV depending on the first-principles method, which is very close to the band gap of silicon (1.17 eV), one of the most important industrial semiconducting materials. It is thus expected that the twin graphene is one of the most promising candidates for the electric devices.
Benefit from its “twin” structure, twin graphene inherits some advanced mechanical properties from graphene. For example, the in-plane stiffness is about 186 N/m (20.9 eV/atom), and the Poisson’s ratio is 0.32. The bending stiffness 1.39 eV is very close to that of graphene [
100]. In particular, twin graphene is able to sustain a large strain of 17% and 16% along the zigzag and armchair directions, respectively. These excellent mechanical properties can support the application of twin graphene in the flexible electric devices.
2.4 Graphyne
The existence of graphyne was conjectured before 1960 [
62], and only Graphdiyne has been synthesized in the experiment [
101]. A large piece of graphdiyne film was grown on copper [
102] and silver substrates [
103]. As stimulated by the success in the experiment, lots of theoretical works have been done to explore mechanical, electric, and possible applications for graphyne. There have been some review articles on a variety of properties of graphyne [
104-
108].
The crystal structure of the graphyne family is shown in Fig. 1D, which displays a flat porous configuration for this material [
109,
110]. Note that the length of the carbon chain (acetylene link) connecting two hexagonal carbon rings can be varied, resulting in unlimited number of possible structures in the graphyne family. The four graphyne members with shortest acetylene links are displayed in Fig. 1D, which are called graphyne, graphdiyne, graphyne-3, and graphyne-4. The carbon chain has four carbons for the specific graphyne shown in Fig. 1D. There are four inequivalent carbon-carbon bonds in graphyne, with bond length 1.432, 1.396, 1.233, and 1.339 Å [
111].
The maximum binding energies for graphyne and graphdiyne are obtained at 7.95 and 7.78 eV/atom, respectively. While the values are 8.66 eV/atom for graphite and 8.22 eV/atom for fullerene C
60, as conpared [
112]. The calculated energy [
113],
E, for graphdiyne is 0.803 eV/atom with respect to graphene, whereas E values for diamond, graphite, (6, 6) nanotube, fullerene C
60, and carbyne are
,
, 0.114, 0.364, and 1.037 eV/atom, respectively.
The Young’s modulus for the armchair and zigzag directions in graphyne was predicted to be 170.4 N/m and 224.0 N/m by Buehler
et al [
114]. The ultimate stresses along the armchair and zigzag directions for graphyne are 48.2 GPa and 107.5 GPa in Ref. [
114], or 14.437 N/m and 20.471 N/m in Ref. [
115]. The ultimate strains for the armchair and zigzag directions for graphyne are 0.0819 and 0.1324 in Ref. [
114], or 0.112 and 0.177 in Ref. [
115]. Peng
et al. demonstrated the elastic deformation of the graphyne in a wide strain range [
116]. Defects are found to be important for various mechanical properties of the graphyne [
117]. Graphyne sheet undergos three continuous phase transitions as the timperature increases from 2800 K to 5000 K, where some defects such as 4, 5, 7, 8 and 9-membered rings and large holes start to appear at 1700 K [
118].
Due to its porous configuration, the Young’s modulus depends on the length of the acetylene link in graphyne members. For instance, Yue
et al. found that the in-plane stiffness decreases from 166 N/m to 88 N/m with the number of carbon atoms in the acetylene line increasing from 1 to 4 [
119]. Similar size dependence was also found in other works [
120,
121]. Pei
et al. found that the in-plane stiffness is 7.60 eV/Å
2 and the Poisson’s ratio is 0.453 for graphdiyne, respectively [
122]. The ultimate stresses along the armchair and zigzag directions for graphdiyne are 9.538 N/m and 20.835 N/m in Ref. [
115]. The ultimate strains for the armchair and zigzag directions for graphyne are 0.109 and 0.208 in Ref. [
115].
The graphyne members have intrinsic electronic band gap from 0.46 eV to 1.22 eV depending on the first-principles approaches and the length of the acetylene line [
104,
111,
112,
123-
125]. The electronic properties can be tuned by hydrogenation [
126]. Long
et al. predicted that the electron mobility and hole mobility for graphdiyne are
cm
2·V
-1·s
-1 and
cm
2·V
-1·s
-1, respectively [
127,
128]. Similar results are obtained by Xi
et al. [
129]. The electronic properties for graphyne members can be efficiently tuned via adsorption of transion metal [
128] or strain engineering [
130].
Taking advantage of their porous structures, recent studies have proposed some possible applications for the graphyne members in fields such as supercapacitors [
131], anode materials for lithium batteries [
132] or lithium ion battery [
133-
137], the Li ion storage [
138,
139], hydrogen purification [
125], catalysts [
140,
141], hydrogen production in a photoelectrochemical water splitting cell [
142], membrane for hydrogen gas separation [
143],
3 2D metallic carbon
3.1 Phagraphene
In 2015, Wang
et al. [
53] proposed phagraphene to be a new stable carbon allotrope, using the systematic evolutionary structure searching. This planar carbon allotrope is energetically comparable to graphene and more favorable than other carbon allotropes proposed in previous works due to its sp
2-hybridization and dense atomic packing. While it has not been snthesized in lab. The total energy of phagraphene is calculated to be
eV/atom. Compared with graphene (
eV/atom), phagraphene is less stable by just 0.2 eV/atom [
53].
The structure of phagraphene is shown in Fig. 1E with Pmg space group, which is constructed by 5-6-7 carbon rings. From the structural point of view, this structure can be derived from graphene by rotating a C-C bond in the 6-carbon ring. There are five inequivalent C-C bonds in phagraphene with bond length of 1.52 Å (C1–C6), 1.44 Å (C2–C3), 1.41 Å (C3–C4), 1.43 Å (C4–C5), and 1.40 Å (C5–C6) [
53].
The elastic modulus of phagrphene is calculated to be
GPa (
N/m) along the armchair direction and
GPa (
N/m) along the zigzag direct (assuming an effective thickness of 0.335 nm). And the tensile strength of phagraphene is predicted to be around
GPa, independent of loading direction [
144]. The Young’s modulus is 292.9 N/m and the Poisson’s ratio is 0.255 for phagraphene [
81]. The ultimate tensile strength is 25.39 N/m at the strain of 0.18 and 25.57 N/m at the strain of 0.16 in the armchair and zigzag directions, respectively [
81]. The mechanical properties for phagraphene can be tuned by hydrogenation or fluorination [
145].
Phagraphene has a special electronic band structure with distorted Dirac cones, and these Dirac cones are robust against external strain [
146]. Phagraphene is a conductor with zero band gap. The electronic band gap can be opened through hydrogenation and fluorination, which can induce band gaps of 4.29 eV and 3.23 eV for the fully hydrogenation and fluorination, respectively [
145]. The electronic band structure for the phagraphene nanoribbon can be manipulated by B-, N-, and BN-doping [
147], or free edges [
148,
149].
The thermal conductivity is found to be anisotropic with
W·m
-1·K
-1 and
W·m
-1·K
-1 along the armchair and zigzag directions, respectively [
144]. First-principles calculations predict phagraphene to be a promising candidate for anode materials in Li ion battery, with higher mobility and Li ion storage capacity [
150].
3.2 Biphenylene
Biphenylene was first synthesized by Lothrop in 1941, which is a planar configuration containing 4-, 6-, and 8-carbon rings [
151]. As shown in Fig. 1F, we can obtain the structure from graphene by rotating two C-C bonds in the 6-carbon ring to get the 4- and 8-carbon rings. there are two types of carbon atoms in biphenylene, with C1 from the 4-carbon ring and C2 from other carbon rings. These two types of atoms form four different C-C bond lengths in biphenylene, i.e., 1.458 Å, 1.454 Å, 1.407 Å, and 1.447 Å, and four different bond angles.
The total energy of biphenylene is calculated to be
eV/atom and the corresponding volues for graphene and phagraphene are
and
eV/atom, respectively. So biphenylene is only 0.47 eV/atom less stable than graphene. Biphenylene is even stable upto a high temperature of 5000 K [
73].
The Young’s modulus is
GPa (
N/m) and
GPa (
N/m) along the armchair and zigzag, respectively (assuming an effective thickness of 0.335 nm) [
73]. Yedla
et al. found that the Young’s modulus is as large as 3.4 TPa (1139 N/m) along the armchair direction, and is 652 GPa (218.4 N/m) in the zigzag direction [
152]. The ultimate stress is 138 GPa at strain 0.064 for the armchair direction, and 23 GPa at strain 0.15 for the zigzag direction [
152].
First-principles calculations have predicted biphenylene to be metallic [
153], and biphenylene nanotubes are also metallic, irrespective to their chiralities [
154]. The biphenylene nanoribbons of zigzag edges are also metallic, while biphenylene nanoribbons with armchair terminations have finite electronic band gap that decreases with increasing ribbon width [
153]. The electronic and excitonic optical absorption properties have also been investigated for the biphenylene [
155]. The electronic band gap for biphenylene can be opened by hydrogenation or fluorination, and a metal to insulator transition was realized by full hydrogenation [
154].
Biphenylene has a large Li ion storage capacity and high Li ion mobility [
150]. Furthermore, the volume expansion during the Li insertion is about 11%. These properties enable biphenylene to be a suitable candidate for the anode of Li ion battery. Biphenylene also has a reasonablly large capacity for the storage of hydrogen [
156].
4 3D metallic carbon
Owing to the outstanding properties of metallic carbon as well as their great potential applications, design and synthesis of three-dimensional metallic carbons are currently one of the hot issues.
4.1 H18 carbon
H
18 carbon was proposed by Zhao
et al. by first-principles calculations in 2015 [
55]. This 3D lattice essentially consists of covalently bonded graphene layers, so the configuration of H
18 carbon can be regarded as a derivative from the graphene as shown in Fig. 1G. This new carbon allotrope has not been synthesized in the experiment, but it may exist in recent detonation experiment as revealed by the XRD simulation [
164]. It is a metallic carbon allotrope with 3D sp
2-sp
3 hybridized bonding network structure in the P6/mmm
symmetry as shown in Fig. 1G. We can obtain the structure through graphene by bond rotation. The lattice parameters are
Å and
Å. There are three inequivalent carbon atoms, forming four bond lengths: 1.591, 1.631, 1.475, and 1.317 Å. The density of H
18 carbon is 3.135 g/cm
3.
The cohesive energy is 8.75 eV/atom, so H
18 carbon is energetically stable. There is no imaginary frequency in the phonon dispersion, which implies the thermal stability for H
18 carbon. The stability of H
18 carbon is further verified by computing the independent elastic constants, which satisfy the Born stability condition [
165]. The bulk modulus and shear modulus of H
18 carbon are 360 GPa and 361 GPa, respectively.
4.2 Hex-C18
Hex-C
18 is a 3D metallic carbon allotrope that was proposed by Liu
et al. with the global structure search [
157]. Although this new carbon allotrope has not been synthesized in the experiment, XRD calculations have provided some clues for the existance of Hex-C
18 in detonation soot [
164,
166]. Figure 2A shows the hexagonal lattice structure for Hex-C
18 with space group P6
3/mcm
. We can also obtain the structure through graphene by some bonds rotations. There are 18 atoms in the unit cell, which are categrized into the sp
3 hybridized carbon C
1 and sp
2 hybridized carbon C
2. The lattice constants are
Å and
Å. Hex-C
18 is a porous structure with low density of 2.41 g/cm
3.
The bulk modulus and shear modulus for Hex-C
18 are 282.9 GPa and 254.6 GPa, respectively. The Vickers hardness of Hex-C
18 was predicted to be 42.2 GPa with Chen’s model [
167] or 80.8 GPa with Gao’s model [
168], so Hex-C
18 is a superhard material. Hex-C
18 was suggested to be a possible anode material, with low Li diffusion energy barriers and high specific capacity.
4.3 Tri-C9
In 2016, Wen
et al. proposed a new 3D metallic carbon allotrope termed as Tri-C
9, by distorting the sp
3 hybridization bond in the diamond structure [
158]. This carbon allotrope has not been verified experimentally, but simulations suggested that it is possibe to produce Tri-C
9 by compressing graphite at pressure above 70 GPa. The lattice structure for Tri-C
9 is shown in Fig. 2B, which has the R32 (No.155) space group. The lattice constants are
Å and
Å, respectively. The bond lengths are 1.526 Å and 1.528 Å, which are close to each other. The six bond angles in tetrahedron of Tri-C
9 are 60°, 105.778°, 105.778°, 122.944°, 122.944°, and 123.807°, which are substantially different from that of diamond (109.5°), indicating the distortion of sp
3 bonds in Tri-C
9.
The stability of Tri-C9 was supported by examing the phonon dispersion and the elastic constants. There is no imaginary frequency in the full phonon dispersion at ambient pressure, so Tri-C9 is thermaldynamically stable. The independent elastic constants satisfy the mechanical stability criteria, which further suggests the mechanical stability of Tri-C9.
First-principles calculations have also predicted some mechanical properties for Tri-C
9. The bulk modulus and shear modulus for Tri-C
9 are 365 GPa and 272 GPa, respectively. The Vickers hardness according to the model by Chen [
167] is 34.8 GPa.
4.4 O-type and T-type carbon
In 2017, Liu
et al. proposed a series of carbon allotropes by self-assembling ultrathin diamond nanostripes, termed O-type and T-type carbon [
159]. Figure 2C shows the configuration for two sampling structures. The allotrope is orthorhombic (Cmmm) or tetragonal (P4
2/mmc) corresponding to the odd or even number of diamond nanostripes that are assembled. The allotrope of orthorhombic symmetry is called O-type carbon, while the allotrope of tetragonal symmetry is called T-type carbon. T10 carbon was found in the experiment [
169]. The thermal and mechanical stability of the O-type and T-type carbons were confirmed by examining phonon dispersions and elastic constants. The cohesive energies of these two allotropes are larger than some other carbon allotropes, like the T-carbon [
160]. A distinct feature for the T-type and O-type carbons is their superhard property that is inherit from diamond.
5 3D superhard carbon
There has been a long-standing search for 3D metallic, superhard carbon, stable under ambient conditions.
5.1 T-carbon
T-carbon was predicted theoretically in 2011 by Sheng
et al. [
160] and has been synthesized in the experiment in 2017 [
170]. It is a semiconductor with electronic band gap 2.25 eV or 2.968 eV [
160]. As shown in Fig. 2D, T-carbon can be regarded as a derivative of the diamond structure by substituting each carbon atom by a carbon tetrahedron. The resultant lattice constant is 7.52 Å and the equilibrium density is 1.5 g/cm
3. The cohesive energy per atom is predicted to be 6.573 eV.
The bulk and shear modulus of T-carbon are 169 GPa and 70 GPa, respectively [
160]. Gao
et al. [
168] found that T-carbon is a superhard material with hardness of 61.1 GPa, while Chen
et al. [
171] predicted the hardness to be less than 10 GPa.
Due to its low density, T-carbon may be used for hydrogen storage, with the gravimetric and volumetric hydrogen capacities around 7.7 wt% and 0.12 kg H
2/L, respectively [
160]. In the meantime, the thermal conductivity of T-carbon is about 33.06 W·m
-1·K
-1, which is much lower than the value of 2388.69 W·m
-1·K
-1 for diamond, because the interaction within T-carbon is weaker than diamond [
172].
5.2 C14-diamond
In 2017, Wu
et al. performed first-principles calculations to predict C
14-diamond [
161]. It is a derivative structure from the diamond as shown in Fig. 2E, and has many similar properties as diamond. For instance, the cohesive energy for C
14-diamond is 8.983 eV/atom, which is very close to the cohesive energy of diamond, so C
14-diamond is thermodynamical stable. The stability of this structure was also confirmed by the generalized elastic stabiliity criteria with nin independent elastic constants. The density of C
14-diamond is 3.37 g/cm
3, also very close to the density of diamond. C
14-diamond is a superhard material with hardness of 55.8 GPa.
5.3 Tetragonal C64
In 2016, Wei
et al. proposed a new semiconductor with band gap of 1.32 eV – tetragonal C
64. [
162]. This carbon allotrope has a tetragonal symmetry (space group I4
1/amd No. 141) as shown in Fig. 2F. The lattice constants are
Å and
Å for the tetragonal unit cell. The density of tetragonal C
64 is 2.562 g/cm
3. The stability of the tetragonal crystal is verified using the criteria for the mechanical stability [
173] with six independent elastic constants. Furthermore, the there is no imaginary frequency for the phonon dispersion over the whole Brillouin zone, indicating the dynamical stability of tetragonal C
64.
According to the model by Lyakhov
et al. [
174], the hardness for tetragonal C
64 is 60.2 GPa, so it is a superhard material. But the hardness is 34.0 GPa from the model by Chen
et al. [
167]. The Young’s modulus is 510 GPa and the Poisson’s ratio is 0.178. The bulk modulus B and the shear modulus G are 264 GPa and 217 GPa within the Voigt–Reuss–Hill approximation, so the B/G ratio is 1.22. This B/G ratio is lower than 1.75, indicating a brittle character [
175]. The weakest tensile strength is 48.1 GPa in the [
111] direction, while the strongest tensile strength is 78.8 GPa in the [
110] direction. The lowest shear strength is 48.2 GPa in the (100)[010] direction and the highest shear strength is 60.1 GPa in the (111)[-1-12] direction.
5.4 C20 T-carbon
Using first-principles calculations, Wang
et al. proposed the structure for C
20 T-carbon and discussed its stability [
163]. It is an insulator with electronic band gap of 5.44 eV. The structure of C
20 T-carbon is shown in Fig. 2G, which has a space group P2
13 (No. 198). The lattice constant is 4.495 Å. There are twenty carbon atoms in the primitive unit cell with three inequivalent carbon atoms, resulting in four different C-C bond length: 1.541 Å, 1.548 Å, 1.550 Å, and 1.541 Å.
C20 T-carbon was predicted to be stable, as there is no imaginary frequency in the whole phonon dispersion. The stability of C20 T-carbon is also verified by computing these three independent elastic constants in cubic lattice, which satisfy the Born stabiligy condition. The bulk and shear modulus is 395 GPa and 427 GPa according to the Voigt–Reuss–Hill approximation. The B/G ratio 0.93 is lower than 1.75, so C20 T-carbon is brittle. The Poisson’s ratio of C20 T-carbon is 0.11. One advanced property for C20 T-carbon is its high hardness of 72.76 GPa, so this is a superhard material. The weakest tensile strength of C20 T-carbon is about 71.1 GPa in the direction and The lowest shear strength of C20 T-carbon is 54.5 GPa.
5.5 P-carbon
Recently, P-carbon was proposed by Pan
et al. using the evolutionary partical swarm structural search approach [
183]. It is a semiconductor with direct band gap of 3.52 eV [
183]. We can design new carbon allotropes from not only graphene and diamond but also carbon nanotube. P-carbon can be obtained by compressing large-diamenter single-walled carbon nanotubes, and the stability of the structure has been verified theoretically [
173]. Figure 3 shows the configuration of P-carbon, which has a tetragonal symmetry (space group I4/mmm, No. 139) with lattice parameters
Å and
Å.
The average bond length is 1.54 Å. As a result of the microporous structure, P-carbon has a low density of 3.1 g/cm
3 [
183]. P-carbon is stable at high pressure and can be even more stable than graphite at pressure above 115 GPa [
183].
The bulk moduli for P-carbon is 334.2 GPa, while the shear moduli is 360.1 GPa. Specifically, P-carbon is a superhard material with hardness predicted to be 86.7 GPa. The tensile strengths along [001], [
100], and [
110] directions are 109.33 GPa at strain 0.32, 37.39 GPa at strain 0.24 and 42.46 GPa at strain 0.24, respectively [
183].
6 Future prospects and summary
6.1 Searching for more carbon allotropes
6.1.1 Theoretical design of more carbon allotropes
By comparing the structural topology of existing carbon allotropes in the above, we have learn how to design new carbon allotropes by manipulating the elemental structures (carbon atoms or C-C bonds) in the two most fundamental carbon materials – graphene and diamond. The derivative structure typically inherits the dimensionality of the mother structure; i.e., allotropes derived from graphene (diamond) are always in 2D (3D) configuration. New carbon allotropes can be predicted in a systematic way along this line with the assistance of the first-principles calculations. For example, Enyashin and Ivanovskii have systematically designed twelve 2D carbon allotropes based on graphene by extending the primitive unit cell [
17].
In particular, it is possible to construct more complex structures by recursively using this designing technique. As an example, the twin-graphene can be theoretically regarded as the covalent sticking of two graphene single layers. Repeating this twin process will yield a new structure with a finite thickness about four-fold of the thickness of graphene. As a result, the theoretical designing technique can lead to endless carbon based allotropes.
6.1.2 More experimental studies needed
Sometimes, new allotropes are first observed in the experiment and then analyzed theoretically, such as the discovery of C
60 buckyball. Several decades ago, it is still a challenge for theoreticalists to imagine the exact atomic structure of a newly discovered allotrope. For example, Kroto and heath made a great effort to figure out the novel structure of the buckyball C
60 [
10]. However, with the fast development of the observation techniques, nowadays it is able to examine the structure of new materials on the atomic level, so it becomes an easier task for theoreticalists to determine the actual configuration of the new allotrope after it is synthesized in the labotory.
On the other hand, theoreticalists are now capable to design new carbon allotropes in a systematic manner, mainly thanks to the increase of the computational speed and the fast development of the first-principles calculations approach. Many carbon allotropes of novel imaginary configurations have been predicted theoretically. A key challenge for these predicted carbon allotropes is the experimental verification, because there are usually lots of different allotropes in the product of the experiment. One successful example is the T-carbon, which was first predicted theoretically in 2011 and verified experimentally in 2017. More experiments are needed to verify the theoretical predictions of these new carbon allotropes.
6.2 Properties and applications of carbon allotropes
The properties and applications of each carbon allotrope were discussed separately in original works, which have been summarized in Tables 2 and 3. We have compared most fundamental properties for carbon allotropes in this review, which can facilitate a comparative investigation for possible applications of these carbon allotropes. The first meaningful task is to complete these blank fields in Table 2 by examining these properties that have not been studied.
Another important direction is to compare properties or applications for these carbon allotropes. For instance, most 2D carbon allotropes are proposed to be good electrode materials for the Li ion battery. A natural question arises: who is the best? A systematic study is needed to address such questions.
Furthermore, 2D materials are usually stacked into finite-thickness functional devices in practical applications. Some 2D carbon allotropes have quite different electronic and optical properties, and these varying physical properties can be integrated by stacking in the vertical van der Waals heterostructure form. This stacking technique will lead to the creation of full carbon multi-functional devices.
6.3 Summary
We have summarized some major properties discussed in the present review in Tables 2 and 3. There are general features that can be extracted from the comparison of some properties listed in these tables. Carbon allotropes can be obtained from the two most fundamental carbon materials – graphene and diamond, by manipulating the elemental structures (carbon atoms or C-C bonds); eg., bond rotation, bond extension, or bond reduction for 2D allotropes. More specifically, phagraphene, biphenylene and graphenylene have the lowest relative energy to graphene, and they are obtained through bond rotation. Graphynes and twin graphene can be obtained through the bond extension, which have the medium energy. The penta-graphene is obtained through the bond reduction, which has the highest energy. According to the cohesive energy of these 2D allotropes, the bond rotation is the most feasible manipulation energetically while bond reduction is the most unfeasible.
There is a strong correlation between the Young’s modulus and the density of the 2D allotropes in Table 2; i.e., the Young’s modulus is linearly proportional to the density as displayed in Fig. 4. Note that twin graphene deviates from the correlation, because the twin graphene essentially has two layers in the out-of-plane direction. In other words, twin graphene is not an ideal 2D allotrope. Similar correlation is also found for the 3D carbon allotropes listed in Table 3; i.e., the Young’s modulus and Shear modulus are both linearly proportional to the density of the material as displayed in Figs. 5 and 6. The phenomenon of increasing modulus (
E) with increasing density (
ρ) has also been found in other materials, like the bones [
184], where the power factor (
α) varies in different models
. According to Figs. 4, 5 and 6, the power factor is
for the 2D and 3D carbon allotropes.
Higher Education Press and Springer-Verlag GmbH Germany, part of Springer Nature 2019