# Penta-graphene: A new carbon allotrope

^{a}Center for Applied Physics and Technology, College of Engineering, Peking University, Beijing 100871, China;^{b}Collaborative Innovation Center of Inertial Fusion Sciences and Applications, Ministry of Education, Beijing 100871, China;^{c}Department of Physics, Virginia Commonwealth University, Richmond, VA 23284;^{d}National Laboratory for Infrared Physics, Shanghai Institute of Technical Physics, Chinese Academy of Sciences, Shanghai 200083, China;^{e}Synergetic Innovation Center of Quantum Information & Quantum Physics, University of Science and Technology of China, Hefei, Anhui 230026, China; and^{f}Institute for Materials Research, Tohoku University, Sendai, 980-8577, Japan

See allHide authors and affiliations

Edited by Ho-kwang Mao, Carnegie Institution of Washington, Washington, DC, and approved January 5, 2015 (received for review August 28, 2014)

## Significance

Carbon has many faces––from diamond and graphite to graphene, nanotube, and fullerenes. Whereas hexagons are the primary building blocks of many of these materials, except for C_{20} fullerene, carbon structures made exclusively of pentagons are not known. Because many of the exotic properties of carbon are associated with their unique structures, some fundamental questions arise: Is it possible to have materials made exclusively of carbon pentagons and if so will they be stable and have unusual properties? Based on extensive analyses and simulations we show that penta-graphene, composed of only carbon pentagons and resembling Cairo pentagonal tiling, is dynamically, thermally, and mechanically stable. It exhibits negative Poisson's ratio, a large band gap, and an ultrahigh mechanical strength.

## Abstract

A 2D metastable carbon allotrope, penta-graphene, composed entirely of carbon pentagons and resembling the Cairo pentagonal tiling, is proposed. State-of-the-art theoretical calculations confirm that the new carbon polymorph is not only dynamically and mechanically stable, but also can withstand temperatures as high as 1000 K. Due to its unique atomic configuration, penta-graphene has an unusual negative Poisson’s ratio and ultrahigh ideal strength that can even outperform graphene. Furthermore, unlike graphene that needs to be functionalized for opening a band gap, penta-graphene possesses an intrinsic quasi-direct band gap as large as 3.25 eV, close to that of ZnO and GaN. Equally important, penta-graphene can be exfoliated from T12-carbon. When rolled up, it can form pentagon-based nanotubes which are semiconducting, regardless of their chirality. When stacked in different patterns, stable 3D twin structures of T12-carbon are generated with band gaps even larger than that of T12-carbon. The versatility of penta-graphene and its derivatives are expected to have broad applications in nanoelectronics and nanomechanics.

Carbon is one of the most versatile elements in the periodic table and forms a large number of allotropes ranging from the well-known graphite, diamond, C_{60} fullerene (1), nanotube (2), and graphene (3) to the newly discovered carbon nanocone (4), nanochain (5), graphdiyne (6), as well as 3D metallic structures (7, 8). The successful synthesis of graphene (3) has triggered considerable interest in exploring novel carbon-based nanomaterials. A wealth of 2D carbon allotropes beyond graphene has since been studied (see *SI Appendix*, Table S1 for details). Although some of these polymorphs such as graphdiyne (6) are metastable compared with graphene, they have been successfully synthesized. Moreover, some 2D carbon allotropes are predicted to exhibit remarkable properties that even outperform graphene, such as anisotropic Dirac cones (9), inherent ferromagnetism (10), high catalytic activity (6), and potential superconductivity related to the high density of states at the Fermi level (11). These results demonstrate that many of the novel properties of carbon allotropes are intimately related to the topological arrangement of carbon atoms and highlight the importance of structure–property relationships (12).

Pentagons and hexagons are two basic building blocks of carbon nanostructures. From zero-dimensional nanoflakes or nanorings (13) to 1D nanotube, 2D graphene, and 3D graphite and metallic carbon phases (7, 8), hexagon is the only building block. Extended carbon networks composed of only pentagons are rarely seen. Carbon pentagons are usually considered as topological defects or geometrical frustrations (14) as stated in the well-known “isolated pentagon rule” (IPR) (15) for fullerenes, where pentagons must be separated from each other by surrounding hexagons to reduce the steric stress. For instance, C_{60} consists of 12 pentagons separated by 20 hexagons forming the shape of a soccer ball, which is a perfect footnote to IPR. The emergence of carbon pentagons is also found to be accompanied by carbon heptagons in some cases (14), but are separated from each other. Inspired by the synthesis of pure pentagon-based C_{20} cage (16), considerable effort has been made to stabilize fused-pentagon-based and non-IPR carbon materials in various dimensionalities (10, 15). Some non-IPR fullerenes have been experimentally realized (15). A “pentagon-first” mechanism was postulated in the transformation from *sp* carbon chains to *sp*^{2} carbon rings during surface growth of 2D carbon sheets (17). Thus, we conceived the idea of building 2D carbon sheets using fused pentagons as a structural motif. In this work, we show that a 2D carbon allotrope, penta-graphene, consisting entirely of pentagons, can indeed exist. The dynamical, thermal, and mechanical stability of this unique structure is confirmed by a series of state-of-the-art theoretical calculations. In addition, we show that pentagon-based carbon nanotubes, penta-tubes, formed by rolling up the penta-graphene sheet, and 3D twin structures of the recently reported T12-carbon (18) formed by stacking these sheets in different patterns, are both dynamically and thermally stable. We demonstrate that these exotic pentagon-based carbon materials exhibit interesting mechanical and electronic properties.

## Results

### Penta-Graphene Exfoliated from T12-Carbon.

Our search for an all-pentagon-based 2D carbon sheet began by examining the recently proposed T12-carbon phase (18) that can be acquired by heating an interlocking-hexagon–based metastable carbon phase at high temperature (8). We note that there are two kinds of C–C bond lengths, namely, the slightly shorter intralayer bond (*d*_{1}) and the slightly longer interlayer bond (*d*_{2}), as shown in Fig. 1*A*. Here the atoms displayed in ball–stick model and highlighted in yellow form a layered structure which can be chemically exfoliated from the T12-carbon phase. In fact, a monolayer can be fabricated from either a layered structure or a nonlayered structure (19). For van der Waals coupled layered structures like graphite, mechanical exfoliation (3) is sufficient to obtain the monolayer sheet (graphene), whereas for chemically bonded bulk phases like MAX (a family of transition metal carbides or nitrides), chemical exfoliation technique has been developed to extract a single layer (MXene) (20).

The optimized crystal structure of the 2D pentagon-based phase generated by exfoliating a single layer from T12-carbon is shown in Fig. 1*B*. The structure possesses *P*-42_{1}*m* symmetry (space group no. 113) with a tetragonal lattice. The optimized lattice constants are *a* = *b* = 3.64 Å. The top view shows that the new phase is composed entirely of carbon pentagons, forming a beautiful pattern well known as Cairo pentagonal tiling (21). From the side view a buckling (0.6 Å) is observed, leading to a 2D sheet with a total thickness of 1.2 Å. This structure can be considered as a multidecker sandwich, with the 4-coordinated C atoms highlighted in yellow in Fig. 1*B* sandwiched between the 3-coordinated atoms. For convenience of discussion, we hereafter group the *sp*^{3}- and *sp*^{2}-hybridized C atoms as C1 and C2, respectively, and call this new graphene allotrope penta-graphene. The unit cell of penta-graphene contains six atoms as denoted by red dashed lines in Fig. 1*B* in which the C1 to C2 ratio is 1:2. The C1–C2 (1.55 Å) and C2–C2 (1.34 Å) bond lengths show pronounced characters of single and double bonds, respectively, and the bond angle θ_{C2–C1–C2} is 134.2°, indicating the distorted *sp*^{3} character of C1 atoms. Such bond multiplicity (22) of carbon, although absent in the well-known diamond, graphite and graphene, has been found in a number of carbon structures with different dimensionalities (6⇓–8, 22, 23) and is of general chemical interest as it leads to intermediate valency (23). Interestingly, we note that penta-graphene resembles the structure of experimentally identified layered silver azide (AgN_{3}) (24). By replacing the N_{3} moieties and Ag atoms with the triconnected C dimers and tetra-connected C atoms, respectively, the geometry of penta-graphene can be realized.

### Energetic Stability.

Total energy calculations are performed to investigate the thermodynamic stability of penta-graphene. Although this phase is metastable compared with graphene and previously reported 2D carbon allotropes (6, 11, 14) due to its violation of the IPR, it is more stable than some nanoporous carbon phases such as 3D T-carbon (25), 2D α-graphyne (6), and (3, 12)-carbon sheet (26) (Fig. 2*A*). We also note that penta-graphene is energetically preferable over some experimentally identified carbon nanostructures such as the smallest fullerene, C_{20} and the smallest carbon nanotube, implying that the 2D penta-graphene sheet could be synthesized. Although C_{20} cage and penta-graphene share the structural motif of fused pentagons, unlike the highly curved C_{20} cage where all of the C atoms exhibit distorted *sp*^{2} hybridization leading to a large strain energy, in penta-graphene the onset of *sp*^{3} hybridization lowers the curvature of fused carbon pentagons, thus partially releasing the strain.

### Dynamic Stability.

Next we study the lattice dynamics of penta-graphene by calculating its phonon dispersion. The results are presented in Fig. 2*B*. The absence of imaginary modes in the entire Brillouin zone confirms that penta-graphene is dynamically stable. Similar to the phonons of graphene (27, 28), there are three distinct acoustic modes in the phonon spectra of penta-graphene. The in-plane longitudinal and in-plane transverse modes have linear dispersion near the Γ point, whereas the out-of-plane (ZA) mode has quadratic dispersion when ** q** approaches

**0**. The quadratic ZA mode in the long-wavelength region is closely associated with the bending rigidity and lattice heat capacity of the nanosheet, which is discussed in detail in

*SI Appendix*, text S1. A remarkable phonon gap is observed in the phonon spectra. Detailed analysis of the atom-resolved phonon density of states (PhDOS) reveals that the double bonds between the

*sp*

^{2}hybridized C2 atoms are predominant in the dispersionless high-frequency modes (Fig. 2

*B*), which is quite similar to the phonon modes in earlier reported

*sp*

^{2}–

*sp*

^{3}hybrid carbon structures (7, 8).

### Thermal Stability.

The thermal stability of penta-graphene is examined by performing ab initio molecular dynamics (AIMD) simulations using canonical ensemble. To reduce the constraint of periodic boundary condition and explore possible structure reconstruction, the 2D sheet is simulated by (4 × 4), (6 × 6), and (4√2 × 4√2) *R*45° supercells, respectively. After heating at room temperature (300 K) for 6 ps with a time step of 1 fs, no structure reconstruction is found to occur in all of the cases. Furthermore, we find that the penta-graphene sheet can withstand temperatures as high as 1,000 K, implying that this 2D carbon phase is separated by high-energy barriers from other local minima on the potential energy surface (PES) of elemental carbon. The snapshots of atomic configurations of penta-graphene at the end of AIMD simulations are shown in *SI Appendix*, Fig. S1. The effect of point defects or rim atoms on the stability of the penta-graphene sheet is also studied by introducing mono- and di-vacancies, Stone–Wales-like defect, adatoms, and edge atoms. The results are presented in *SI Appendix*, text S2, Figs. S2–S4, where one can see that the stability and structure of penta-graphene is robust, despite the defects.

### Mechanical Stability.

As we fix the supercell during all of the MD simulations, it is necessary to assess the effect of lattice distortion on structural stability. To guarantee the positive-definiteness of strain energy following lattice distortion, the linear elastic constants of a stable crystal has to obey the Born–Huang criteria (29). We calculate the change of energy due to the in-plane strain to determine the mechanical stability of penta-graphene. For a 2D sheet, using the standard Voigt notation (26), i.e., 1-*xx*, 2-*yy*, and 6-*xy*, the elastic strain energy per unit area can be expressed as_{11}, C_{22}, C_{12}, and C_{66} are components of the elastic modulus tensor, corresponding to second partial derivative of strain energy with respect to strain. The elastic constants can be derived by fitting the energy curves associated with uniaxial and equi-biaxial strains. The curves are plotted in *SI Appendix*, Fig. S5. For a mechanically stable 2D sheet (29), the elastic constants need to satisfy *C*_{11}*C*_{22} − *C*_{12}^{2} > 0 and *C*_{66} > 0. Due to the tetragonal symmetry of penta-graphene, we have *C*_{11} = *C*_{22}. Thus, in this case we only need to satisfy *C*_{11}>|*C*_{12}| and C_{66} > 0. Under uniaxial strain, *ε*_{yy} = 0, *C*_{11} = 265 GPa·nm. Under equi-biaxial strain, *ε*_{xx} = *ε*_{yy}, we have *C*_{11} + *C*_{12} = 247 GPa·nm, hence, *C*_{12} = −18 GPa·nm. Thus, the calculated elastic constants satisfy *C*_{11} > |*C*_{12}|, and the calculated *C*_{66} is positive, indicating that the 2D penta-graphene sheet is mechanically stable.

### Mechanical Properties.

Having confirmed the stability of penta-graphene, we next systematically study its mechanical properties. The in-plane Young’s modulus, which can be derived from the elastic constants by *C*_{12} is negative for this nanosheet, leading to a negative Poisson’s ratio (NPR), viz., ν_{12} = ν_{21} = *C*_{12}/*C*_{11} = −0.068. To confirm this unusual result, we calculated the lateral response in the *y* direction when the lattice endures a tensile strain in the *x* direction. We examine cases with *ε*_{xx} = 5%, 6%, and 7%. As expected, we find that the equilibrium lattice constant in the *y* direction is expanded in all of the cases (Fig. 3*A*). This confirms the NPR of penta-graphene. It is well known that Poisson’s ratio is defined as the negative ratio of the transverse strain to the corresponding axial strain. Normally, this ratio is positive as most solids expand in the transverse direction when subjected to a uniaxial compression. Although the continuum mechanics theory does not rule out the possibility of emergence of NPR in a stable linear elastic material, it is fairly rare to find such NPR material in nature. However, it has been found that some artificial materials have NPR and exhibit excellent mechanical properties (31, 32). Such materials, usually referred to as auxetic materials or mechanical metamaterials, are of broad interest in both scientific and technological communities (33). Thus, penta-graphene with such unusual mechanical property may have multiple applications such as a tension activatable substrate, a nanoauxetic material, or a deformable variable-stiffness material.

Besides in-plane stiffness, ideal strength is also a very important mechanical property of a nanomaterial. We study the ideal strength of penta-graphene by calculating the variation of stress with equi-biaxial tensile strain using different cells. The results are plotted in Fig. 3*B*, which shows that the strain at the maximum stress is 21%. Such an ultrahigh ideal strength is exciting. However, we should note that phonon instability might occur before mechanical failure. Such failure mechanism has been well studied in graphene where the phonon softening induced by Kohn anomaly occurs before the stress reaches its maximum in the primitive cell (27, 28). To check whether similar phonon-dominant failure mechanism exists in penta-graphene, we compute the phonons under increasing equi-biaxial strain. The results at the critical point of phonon softening are plotted in Fig. 3*C*. We find that phonon softening does not arise until the magnitude of equi-biaxial strain reaches 17.2%, which is smaller than the magnitude of 21% obtained from the stress–strain curve. For comparison, we also calculate phonons of the equi-biaxially stretched graphene. The observed softening of the K_{1} mode at the Dirac point under equi-biaxial tensile strain of 14.8% (Fig. 3*D*) is in excellent agreement with previous work (27, 28). This indicates that the critical strain of penta-graphene is significantly larger than that of pristine graphene. It is also comparable to that of carrier-doping–strengthened graphene (28). At the critical strained state, the single bond lengths between C1 and C2 atoms reach ∼1.77 Å, which is comparable with the experimentally (34) and theoretically (35) reported longest C–C bond length. Detailed analyses on the eigenvectors corresponding to the imaginary modes reveal that the structure fracture stems from the breakdown of some of the *σ* bonds between C1 and C2 atoms.

### Electronic Properties.

To probe the electronic properties of penta-graphene, we calculate its band structure and corresponding total and partial density of states (DOS). As shown in Fig. 4*A*, penta-graphene is an indirect band-gap semiconductor with a band gap of 3.25 eV [computed using the Heyd–Scuseria–Ernzerhof (HSE06) functional] (36, 37), because the valance band maximum (VBM) lies on the Γ–X path whereas the conduction band minimum is located on the M–Γ path. However, due to the existence of the sub-VBM on the M-Γ path, which is very close to the true VBM in energy, penta-graphene can also be considered as a quasi-direct–band-gap semiconductor. Analysis of its partial DOS reveals that the electronic states near the Fermi level primarily originate from the *sp*^{2} hybridized C2 atoms, which is further confirmed by calculating the band-decomposed charge density distributions, as shown in Fig. 4 *B*–*E*. A simplified tight-binding model is used to understand the underlying physics behind the band-gap opening feature in the band structure of penta-graphene (see *SI Appendix*, text S3 for details). We argue that it is the presence of the *sp*^{3}-hybridized C1 atoms that spatially separates the *p*_{z} orbitals of *sp*^{2}-hybridized C2 atoms, hindering full electron delocalization and thus giving rise to a finite band gap. The dispersionless, partially degenerate valance bands lead to a high total DOS near the Fermi level, lending to the possibility that Bardeen–Cooper–Schrieffer superconductivity can be achieved in this nanosheet through hole doping (38).

### Penta-Tubes: Rolled-Up Penta-Graphene.

It is well known that the electronic properties of carbon nanotubes are closely related to graphene according to the zone folding approximation (39). Due to the gapless semimetallic feature of graphene, the electronic properties of carbon nanotubes are highly chirality-dependent: a carbon nanotube is metallic only when its chiral vector (*n*, *m*) satisfies *n* − *m* = 3*l*, where *l* is an integer. The difficulty in fabricating and separating carbon nanotubes with certain conductance (metallic or semiconducting) greatly hinders its application in nanoelectronics. A previous study proposed a family of metallic carbon nanotubes based on metallic Heackelite sheet (11). It is therefore natural to expect that the penta-graphene–based nanotubes could be semiconducting regardless of chirality. To test this hypothesis we have constructed a series of pentagon-based carbon nanotubes by rolling up the penta-graphene sheet along the (*n*, *m*) chiral vectors, where *n = m* range from 2 to 8 (Fig. 5*A*). The tubes with other chiralities (*n ≠ m*) failed to converge to stable tubular structures. We name these pentagon-based carbon nanotubes penta-tubes. The optimized geometry of a (3, 3) penta-tube is illustrated in Fig. 5*B*. The dynamic and thermal stability of this nanotube is confirmed by carrying out phonon calculations and AIMD simulations, respectively. The results are presented in Fig. 5*C* and *SI Appendix*, Fig. S7, respectively. We find that not only all of the (*n*, *n*) penta-tubes are dynamically robust (*SI Appendix*, Fig. S8) but also they are thermally stable up to 1,000 K. Analysis of their band structures and DOS reveals that all of the stable penta-tubes are semiconducting. The calculated results are summarized in *SI Appendix*, Table S2. Except for the highly curved (2, 2) penta-tube, the band gaps of the (*n*, *n*) penta-tubes are not sensitive to their diameters. Thus, chirality-independent semiconducting carbon nanotubes can be produced for application in nanoelectronics. The semiconducting behavior of penta-tubes can be attributed to the electronic structure of penta-graphene, which resembles the case of other semiconducting monolayers such as h-BN monolayer and the corresponding nanotubes (40) which inherit the semiconducting feature.

### 3D Carbon Structures: Stacked Penta-Graphene Layers.

To further explore the structural versatility of penta-graphene, we have altered the stacking patterns of the penta-graphene layers, leading to a 3D structure as shown in Fig. 6 *A* and *B*. Following the nomenclature used to analyze the structural character of fused-pentagon-based Pentaheptite (14), we define the layer stacking in T12-carbon as AB type. The stacking of the designed structure shown in Fig. 6 *A* and *B* is then termed as AA type. It can be viewed as a twin structure of the T12-carbon phase. The calculated phonon spectra of AA-T12 are presented in Fig. 6*C*, confirming its dynamical stability. Indeed, more complicated structures are expected to be built, akin to the pentaheptite modification of graphene (14). An example of such structures containing four penta-graphene layers (24 atoms) per unit cell, termed ABAA-T24, is presented in *SI Appendix*, Fig. S9. The detailed structural information of these new 3D phases as well as T12-carbon, for comparison, is given in *SI Appendix*, Table S3. The calculated cohesive energies (averaged on each carbon atom) are −8.87, −8.92, and −8.98 eV for AA-T12, ABAA-T24, and T12-carbon, respectively, indicating that the AA-T12 and ABAA-T24 phases are nearly as stable as T12-carbon. Band structure calculations suggest that both AA-T12 and ABAA-T24 are semiconducting with the energy band gaps of 5.68 and 5.33 eV, respectively, which are even larger than that of T12-carbon (4.56 eV), indicating that these phases can be highly electrically resistant and optically transparent like those products of cold compressed graphite (41). The bulk moduli of these polymorphs are calculated by fitting the third-order Birch Murnaghan equation of states (42). Although the bulk moduli of AA-T12 (359 GPa) and ABAA-T24 (380 GPa) are slightly smaller than that of T12-carbon (403 GPa), they are comparable with that of cubic BN, suggesting their potential applications in machining. We note that a very recent theoretical work identified the AA-T12 structure by considering its Si analog (43). It is also pointed out that, like T12-carbon, AA-T12 is a universal structure shared by C, Si, Ge, and Sn, suggesting this family of tetragonal structures may be ubiquitous in elemental allotropes of group IVA elements.

## Discussion

We have demonstrated via AIMD simulations that the metastable penta-graphene structure can withstand very high temperature. We note that the experimentally synthesized dodecahedral C_{20}, the smallest carbon fullerene consisting of only carbon pentagons, is metastable, but possesses outstanding thermal stability (up to 3,000 K) (44). These results imply that thermodynamic criteria may not be the deterministic factor in the synthesis of new carbon-based materials. In fact, due to the bonding versatility of carbon, the PES might be fairly complicated with numerous local minima (corresponding to metastable phases) separated from each other by considerable energy barriers. Graphite and cubic diamond are energetically superior to almost all other carbon polymorphs, namely, they correspond to the two lowest-lying valleys on the PES. The high-energy barrier between graphite and diamond makes it possible for both graphite and diamond to coexist in nature. However, recent experimental (41) and theoretical (45) advances have identified many intermediate phases between graphite and diamond during cold compression of graphite. Some of these phases are pressure-recovered, i.e., they can exist when the external pressure is removed. This stability is also ascribed to the considerable energy barrier. In fact, even in some surfaces of carbon structures, different structural reconstruction patterns are separated by appreciable kinetic energy barriers (35). These findings highlight the vital role that kinetics (22, 23) plays in carbon structure evolution, and fuel the exploration of new metastable carbon phases as functional materials.

We now reflect on the relationship between the special atomic configuration and exotic mechanical properties of penta-graphene. When the structure is under uniaxial tension, the expansive lateral response has two kinds of impact on the total energy: on one hand, it elongates the bond and weakens the binding energy; on the other hand, it significantly reduces the difference between the two lattice constants *a* and *b*, thus helping the structure to get close to its original tetragonal symmetry, reducing the bond distortion around the *sp*^{3}-bonded C1 atoms and lowering the strain energy. The structure thus evolves to its equilibrium as a result of compromise between these two competing factors. Detailed illustrations of atom evolution under uniaxial lattice stretch are presented in *SI Appendix*, text S4. Such regime is reminiscent of an earlier work (46) which argued that the combination of “chemical criteria” and “crystallographic criteria” in carbon materials can lead to exceptional mechanical properties. The ultrahigh critical strain of penta-graphene is also intimately related to its atomic structure. The buckled structure slows down the bond elongation and hence the structure is highly stretchable. Besides, graphene has topologically protected point-like Fermi surface, and the coupling between the electron states near the Fermi level and certain phonon mode leads to fast phonon softening under biaxial tension (27), whereas penta-graphene does not suffer from Kohn anomaly, because it is semiconducting.

One practical issue in the synthesis of penta-graphene is how to selectively break the interlayer covalent bonds in T12-carbon (18). To address this challenge, we point to a similar strategy where hydrogen intercalation was successfully used to decouple a graphene layer from the H–SiC (0001) surface (47, 48). Details of our exfoliation scheme are given in *SI Appendix*, text S5.

In summary, we showed that a 2D carbon sheet, penta-graphene, composed entirely of pentagons can be obtained by chemically exfoliating a single layer from the T12-carbon phase. Although penta-graphene is energetically metastable compared with graphene, it is dynamically stable and can withstand temperatures up to 1,000 K. Due to its special atomic configuration, penta-graphene has unusual properties, such as (*i*) it exhibits NPR, similar to that recently reported in a single-layer black phosphorus sheet (32); (*ii*) it exhibits ultrahigh ideal strength that can even outperform graphene; (*iii*) it is semiconducting, thus, there is no need to functionalize penta-graphene for opening the band gap as is the case with graphene. In addition, penta-graphene can be rolled up to form a 1D pentagon-based nanotube that is semiconducting regardless of its chirality. Therefore, there is no need to develop special techniques to separate semiconducting nanotubes from the metallic ones as is the case with conventional carbon nanotubes. Penta-graphene can also be stacked to form 3D stable structures displaying different properties from those of the mother-phase T12-carbon. Thus, penta-graphene sheet not only possesses exotic properties by itself but also can be used to build new structures. We hope that these findings will motivate experimental efforts. Once synthesized, these new carbon allotropes may not only enrich carbon science but also may lead to an untold number of applications.

## Methods

First-principles calculations and AIMD simulations within the framework of density functional theory are performed using Vienna Ab initio Simulation Package (VASP) (49). The 2D system is separated from its periodic images by a vacuum distance of 20 Å in the perpendicular direction. Projector augmented wave (PAW) (50) method is used to treat interactions between ion cores and valance electrons. Plane waves with a kinetic energy cutoff of 500 eV are used to expand the valance electron (2*s*^{2}2*p*^{2}) wavefunctions. The exchange-correlation potential is incorporated by using the generalized gradient approximation (51) due to Perdew–Burke–Ernzerhof in most of our calculations whereas a hybrid HSE06 (36, 37) functional is used for high-accuracy electronic structure calculations. The first Brillouin zone is represented by K points sampled using the Monkhorst–Pack scheme (52) with a grid density of 2π × 0.02 Å^{-1}. For geometry relaxation, the convergence thresholds for total energy and atomic force components are set at 10^{−4} eV and 10^{−3} eV Å^{-1}, respectively. In AIMD simulations the convergence criterion of total energy is set as 1 meV. Temperature control is achieved by Nosé thermostat (53). Structure relaxations are performed without any symmetry constraint. Phonon properties are calculated using finite displacement method implemented in Phonopy (54). A (4 × 4) supercell is constructed to calculate the atomic forces by using VASP, with a very high accuracy (stringent energy convergence criteria 10^{−8} eV per unit cell). Phonon calculations using a larger supercell (6 × 6) yield consistent results. For phonon calculations in graphene, an (8 × 8) supercell is used, which has been found in previous work (27) to be sufficient to take into account the long-range interatomic interactions.

## Acknowledgments

The authors thank the crew of the Center for Computational Materials Science, the Institute for Materials Research, Tohoku University (Japan), for their continuous support of the HITACHSR11000 supercomputing facility. This work is partially supported by grants from the National Natural Science Foundation of China (NSFC-51471004, NSFC-11174014, NSFC-10990104, and NSFC-11334008), the National Grand Fundamental Research 973 Program of China (Grant 2012CB921404), and the Doctoral Program of Higher Education of China (20130001110033). S.Z. acknowledges funding from the Graduate School of Peking University that enabled him to visit P.J.'s group at Virginia Commonwealth University, where the present work is partially conducted. P.J. acknowledges support of the US Department of Energy, Office of Basic Energy Sciences, Division of Materials Sciences and Engineering under Award DE-FG02-96ER45579.

## Footnotes

- ↵
^{1}To whom correspondence should be addressed. Email: qianwang2{at}pku.edu.cn.

Author contributions: Q.W. designed research; S.Z. and J.Z. performed research; S.Z., J.Z., Q.W., X.C., Y.K., and P.J. analyzed data; and S.Z., J.Z., Q.W., X.C., Y.K., and P.J. wrote the paper.

The authors declare no conflict of interest.

This article is a PNAS Direct Submission.

This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1416591112/-/DCSupplemental.

## References

- ↵
- ↵
- ↵.
- Novoselov KS, et al.

- ↵
- ↵
- ↵
- ↵
- ↵.
- Zhang S,
- Wang Q,
- Chen X,
- Jena P

- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵.
- Lee C,
- Wei X,
- Kysar JW,
- Hone J

- ↵.
- Burns S

- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵.
- Mao WL, et al.

- ↵.
- Murnaghan FD

- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵

## Citation Manager Formats

## Article Classifications

- Physical Sciences
- Physics