Three-dimensional Pentagon Carbon with a genesis of emergent fermions

Abstract

Carbon, the basic building block of our universe, enjoys a vast number of allotropic structures. Owing to its bonding characteristic, most carbon allotropes possess the motif of hexagonal rings. Here, with first-principles calculations, we discover a new metastable three-dimensional carbon allotrope entirely composed of pentagon rings. The unique structure of this Pentagon Carbon leads to extraordinary electronic properties, making it a cornucopia of emergent topological fermions. Under lattice strain, Pentagon Carbon exhibits topological phase transitions, generating a series of novel quasiparticles, from isospin-1 triplet fermions to triply degenerate fermions and further to Hopf-link Weyl-loop fermions. Its Landau level spectrum also exhibits distinct features, including a huge number of almost degenerate chiral Landau bands, implying pronounced magneto-transport signals. Our work not only discovers a remarkable carbon allotrope with highly rare structural motifs, it also reveals a fascinating hierarchical particle genesis with novel topological fermions beyond the Dirac and Weyl paradigm.

ARTICLE

Received 8 Oct 2016 | Accepted 14 Apr 2017 | Published 5 Jun 2017

Carbon, the basic building block of our universe, enjoys a vast number of allotropic structures. Owing to its bonding characteristic, most carbon allotropes possess the motif of hexagonal rings. Here, with rst-principles calculations, we discover a new metastable three-dimensional carbon allotrope entirely composed of pentagon rings. The unique structure of this Pentagon Carbon leads to extraordinary electronic properties, making it a cornucopia of emergent topological fermions. Under lattice strain, Pentagon Carbon exhibits topological phase transitions, generating a series of novel quasiparticles, from isospin-1 triplet fermions to triply degenerate fermions and further to Hopf-link Weyl-loop fermions. Its Landau level spectrum also exhibits distinct features, including a huge number of almost degenerate chiral Landau bands, implying pronounced magneto-transport signals. Our work not only discovers a remarkable carbon allotrope with highly rare structural motifs, it also reveals a fascinating hierarchical particle genesis with novel topological fermions beyond the Dirac and Weyl paradigm.

DOI: 10.1038/ncomms15641 OPEN

Three-dimensional Pentagon Carbon with a genesis of emergent fermions

Chengyong Zhong1, Yuanping Chen1, Zhi-Ming Yu2, Yuee Xie1, Han Wang3, Shengyuan A. Yang2

& Shengbai Zhang3

1 School of Physics and Optoelectronics, Xiangtan University, Xiangtan, Hunan 411105, China. 2 Research Laboratory for Quantum Materials, Singapore University of Technology and Design, Singapore 487372, Singapore. 3 Department of Physics, Applied Physics, and Astronomy, Rensselaer Polytechnic Institute, Troy, New York 12180, USA. Correspondence and requests for materials should be addressed to Y.C. (email: mailto:[email protected]

Web End [email protected] ) or to S.A.Y. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 8:15641 | DOI: 10.1038/ncomms15641 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15641

Most stable carbon allotropes, including the prominent members such as graphite, diamond, carbon nanotube1 and graphene2, possess a motif of hexagonal rings due

to the dominating sp2 and sp3 bonding character3. In contrast, the motif of pentagonal rings, although can appear in molecules such as fullerene4, is very rare in three-dimensional (3D) allotropes. This is partly due to the unfavourable bonding angles that may easily cause structural instability5, and also because the pentagonal motif does not readily t into the 3D space group symmetry6. It seems unlikely to have a stable 3D carbon allotrope dominated by pentagonal carbon rings.

Meanwhile, the research on topological metals and semimetals (TMs) has been attracting tremendous interest714. In these materials, topologically protected quasiparticles emerge at protected band-crossing points around Fermi Level, giving rise to condensed matter realizations of many peculiar fermions with a variety of distinct electronic properties. For example, in Dirac and Weyl semimetals, the quasiparticles near the protected linear band-crossing points are direct analogues of the Dirac and Weyl fermions in relativistic quantum eld theory, providing unique opportunities to study interesting high-energy physics phenomena such as chiral anomaly and gravitational anomaly in a table-top experiment1518. The emergent topological fermions can be classied by the degeneracy of the band-crossing point, which dictates the number of internal degrees of freedom of the fermion. Dirac and Weyl fermions are associated with crossing points with fourfold and twofold degeneracies19,20. Importantly, the breaking of Lorentz symmetry at lattice scale makes it possible to have TMs that can host new types of fermions with different degeneracies beyond the Dirac and Weyl fermions21-23. These new fermions are now actively searched for, and if realized, they are believed to host a wealth of exotic physical properties unprecedented in high-energy physics.

The search for TMs has been focused on compounds involving heavy elements7,24,25, hoping that the strong spinorbit coupling could help to achieve a nontrivial band structure. On the other hand, light elements such as carbon may offer a distinct alternative. For example, graphene is well-known with two-dimensional (2D) Dirac fermions26. With negligible spinorbit coupling, the real spin in carbon materials can be considered as a dumb degree of freedom27. Then the fundamental time reversal operation (T ) here satises T 2 1, contrasting

with T 2 1 for the spinful case. Thus in terms of topological

classications, the nontrivial phases in carbon would be fundamentally distinct from those in spinorbit-coupled systems. Indeed, previous studies have identied a few TM carbon allotropes with Weyl points28,29, nodal loops27,2932 and even Weyl surfaces28. Motivated by these observations, one may naturally wonder: is it possible to also realize new types of topological fermions in carbon allotropes?

In this work, based on rst-principles calculations, we propose a new 3D carbon allotrope, named as 3D Pentagon Carbon (3D-C5), which is entirely composed of pentagonal carbon rings. With nicely interlaced structures, the structure shows good stability comparable to other metastable carbon allotropes. We characterize its phononic, mechanical and electronic properties. At the equilibrium state, Pentagon Carbon is a narrow-gap semiconductor. A semiconductormetal quantum phase transition can be driven by an applied strain. Remarkably, at the critical transition point, three bands cross at a single Fermi point in the energy-momentum space, making the emergent low-energy fermions a novel type of isospin-1 triplet fermions. Beyond the transition, the single Fermi point splits into two triply degenerate band-crossing points along the fourfold screw-rotational axis. The material is turned into a novel TM possessing

another type of triply degenerate fermions. When further breaking the screw-rotational symmetry, the triply degenerate fermion points will transform into two inter-connected Hopf-link Weyl loops protected by the remaining mirror symmetry. We show that, for each case, the Landau level (LL) spectrum exhibits distinct features tied to the topological fermions, suggesting pronounced anomalies in magneto-transport experiment. In addition, we discuss a possible reaction pathway for the experimental synthesis of Pentagon Carbon.

ResultsLattice structure and stability. The structure of 3D Pentagon Carbon is shown in Fig. 1a , which can be viewed as formed by interlinking two orthogonal arrays of pentagon-ring nanoribbons and stacking them along the z-direction (Fig. 1b). After linking, the two nanoribbons share one (green-colored) atom. Hence in the Pentagon Carbon lattice, there are two kinds of carbon atoms, that is, the sp2-hybridized (red- and blue-colored) atoms forming the backbones of the pentagon-ring nanoribbons and the sp3-hybridized (green-colored) atoms forming the linkage between ribbons.

The structure has the nonsymmorphic D194h space group symmetry (No. 141, I41/amd). Here an important symmetry element is the screw rotation

C4z, which is a fourfold rotation along z followed by a factional lattice translation of c^z=4, where c is the lattice parameter in the z direction. The detailed structural properties are presented in Table 1 and are compared with several other (meta-) stable 3D carbon allotropes. One observes that the cohesive energy of Pentagon Carbon, although slightly higher than graphite and diamond, is less than T6 (ref. 33) and Bcc-C8 (ref. 34) and is comparable with Carbon Kagome Lattice (CKL)35,36 structure. Due to its porous structure, the mass density of Pentagon Carbon is low (B2.18 g cm 3) comparable with that of graphite. This is also correlated with its relatively low bulk modulus B257 GPa compared with other allotropes (except for graphite), suggesting that external strain can be more easily applied. The phonon spectrum of Pentagon Carbon is also investigated, which shows no soft mode, conrming its dynamical stability (more analysis of lattice structure and stability can be found in Supplementary Fig. 1 and Note 1).

Electronic structure and quantum phase transition. Remarkably, Pentagon Carbon is also a material hosting a series of novel topological fermions beyond the usual Dirac and Weyl fermions. The band structure in Fig. 2a shows that, in equilibrium state, it is a narrow-gap semiconductor with a small energy gap B21.4 meV. The electronic structure shows strong anisotropy between in-plane (xy) and out-of-plane (z) directions in accordance with the crystal structural anisotropy. Focusing on the G point, there are two valence bands (heavy-hole and light-hole)

and a single conduction band at low energy, with the conduction band and the light-hole valence band almost crossing linearly. Note that the degeneracy of the two valence bands at G is dictated by symmetry: the two states belong to a 2D irreducible representation Eg of the D4h group at G. In comparison, the conduction band belongs to the one-dimensional (1D) A1u

representation. Another salient feature is that the two valence bands are degenerate along the GZ path (kz-axis). This is because that there is remaining C4v little group symmetry along this screw axis, such that the 2D representation does not split according to the compatibility relation. From the analysis of wave-function and projected density of states (PDOS), we nd that the low-energy states are mainly from the p orbitals of the sp2 atoms (Fig. 2f, Supplementary Fig. 2).

2 NATURE COMMUNICATIONS | 8:15641 | DOI: 10.1038/ncomms15641 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15641 ARTICLE

a b c

+

C1 b

z

C

y x

d e

Z

P

y

X

N

x

z

Figure 1 | Crystal structure and Brillouin zone of 3D Pentagon Carbon. (a) Structure of 3D Pentagon Carbon, where the purple lines depict the conventional unit cell, and the red and blue Carbon atoms form mutually orthogonal armchair chains which are linked by the green atoms. This structure can also be viewed as composed of orthogonal pentagonal rings sharing the vertex green atoms, as shown in b. (c) Primitive cell of 3D Pentagon Carbon. The red and blue atoms are marked as C1 sites and the green atoms are marked as C2 sites, indicating their different orbital hybridization character. The lattice parameters of the primitive cell are a b c 7.15 , a b 149.86, g 43.11. (d) Top view (from z-axis) of the structure, showing that the red

and blue armchair chains are mutually orthogonal. (e) Brillouin zone of the primitive cell along with the high-symmetry points marked.

Since the semiconducting gap is small, a slight perturbation may close the gap and drive a semiconductor-to-metal phase transition (a tight-binding model analysis is given in Supplementary Figs 3,4 and Note 2). In Fig. 2bd, we show band structures for three applied uniaxial strains e 1, 0.1 and 1% along the z direction. One observes that indeed a semiconductormetal phase transition is induced. A tensile strain drives the system to a semiconductor with wider gap, whereas a compressive strain drives the system towards a metal. Of particular interest is the change of band ordering at the G point.

The A1u singlet state and the Eg doublet switch their positions in energy across the transition. Due to their different symmetry character, A1u and Eg states cannot hybridize in the transition.

Hence at the critical strain, we must have the three bands crossing at a single point, where the conduction and the light-hole bands have linear dispersion in the kxky plane; and all three bands have quadratic dispersion along the screw axis. The 3D band structures in Fig. 2i,j clearly show the anisotropic feature around the triplet point.

Passing the critical point, the band ordering is inverted at G. However, the corresponding states at Z on the Brillouin zone boundary still have the original ordering. This band topology means that the two bands must cross each other between G and Z.

Further, because the two bands still belong to different representations along the screw axis, they must cross in a linear manner with two triply degenerate band-crossing points (one on each side of G), as conrmed in Fig. 2d. Therefore, the metallic phase represents a novel TM phase with a pair of triply degenerate band-crossing points near the Fermi level. In such sense, the quantum phase transition may also be regarded as a topological phase transition from a trivial semiconductor phase to a TM phase.

We can dene a Z2 invariant to characterize this TM phase and the phase transition. Let D and Dz be the values of the energy difference between A band and E band at G point and at Z point,

respectively. Then a Z2 invariant can be dened as w sgn D DZ

2 f 1; 1g. w 1 corresponds to the

trivial case where a band-crossing point does not need to occur, while w 1 is the topologically nontrivial case where a triply

degenerate crossing point must exist between G and Z. The topological phase transition from a trivial semiconductor phase to TM phase (as shown in Fig. 2) is thus characterized by w changing from 1 to 1. The topological character is manifest in the

sense that as long as the Z2 invariant w is well dened (that is, symmetry is preserved) and does not change, the existence of the triply degenerate crossing point and the linear band dispersion along kz near the point are guaranteed and insensitive to system parameters and perturbations.

Emergent fermions and k . p model. Nontrivial band-crossing points are fascinating objects because the fermionic quasiparticles around these points necessarily have a multi-component structure, distinct from the usual Schrodinger fermions. The essential physics of these emergent fermions are characterized by the low-energy effective k . p model expanded around the crossing point.

Around the G point, we can construct the k . p model based on the three basis states of A1u and Eg. Constrained by the D4h symmetry group and the time reversal symmetry for a spinless system (as spinorbit coupling is negligible for carbon), one obtains the k . p model up to k-quadratic order as

H k

C D1k2z D2k2?

2 4

D B1k2z B2k2? iAkx iAky

iAkx 0 B3kxky iAky B3kxky 0

3 5;

1

where the coefcients A, B1, B2, C, D1, D2, and D can be

determined by tting the density functional theory (DFT) result. At transition point, D 0, the three bands cross at a single point,

NATURE COMMUNICATIONS | 8:15641 | DOI: 10.1038/ncomms15641 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15641

Table 1 | Comparison between Pentagon Carbon and other representative carbon allotropes.

System Space Group d (sp2) d (sp3) h (sp2) h (sp3) q B Ecoh CKL F63/MMC 1.53, 1.50 60.0, 115.1, 117.7 2.75 327 8.81

Bcc C8 P4/MMM 1.52, 1.58, 1.59 90.0, 117.6, 139.1 2.97 343 8.49

3D-C5 I41/amd 1.37, 1.44 1.53 109.7, 110.3, 140.1 100.1, 114.3 2.18 257 8.79

T6 P42/MMC 1.34 1.54 115.1, 122.5 106.7, 115.1 2.95 350 8.63

Diamond FD

3M 1.54 109.5 3.55 431 9.09

Graphite F63/MMC 1.42 120 2.24 36 9.22

The table contains information of the space group, the bond length d (), the bond angle y (), the mass density r (g cm 3), the bulk modulus B (GPa) and the cohesive energy E (eV per atom).

a

e f

DFT TB

Energy (eV)

1.4

0.7

0.0

1.4

0.7

[afii9797]

A1u

Eg

N X Z P

b

g h

Energy (eV)

0.5

0.0

0.5

A1u

1 %

Eg

A1u

X

X Z

0.5

0.0

0.5

c

Energy (eV)

Eg

i j

k l

0.1 %

N X

X Z

d

Energy (eV)

0.5

0.0

0.5

Eg

A1u

X Z

1 %

N X

N X Z

Figure 2 | Electronic band structure of 3D Pentagon Carbon. (a) Band structure of 3D Pentagon Carbon, as from DFTcalculation (blue solid line) and from tight-binding model (red dashed line). The solid and dashed lines are DFTand tight-binding results, respectively. (bd) The band structures of 3D Pentagon Carbon under uniaxial strain along z-direction with e 1%, 0.1% and 1%, respectively. (e) The zoom-in view of band structure around the G point in a

along NGX. The two degenerate valence bands at G belong to the two-dimensional irreducible representation Eg, and the lowest conduction band belongs to the one-dimensional A1u representation. (f) The charge density distribution for states near the Fermi level in a, indicating these states are mainly from p orbitals of the sp2 hybridized atoms. (gl) The 3D band structures on the NGX and the XGZ planes corresponding to bd, respectively.

as illustrated in Fig. 3a. Then we can further simplify the model

by keeping only the k-linear terms:

H k

Ak k; 2 where k kx; ky; 0, and

lx

2 4

0 i 0

i 0 0

0 0 0

3

5; ly

2 4

0 0 i

0 0 0

i 0 0

3

5; lz

2 4

0 0 0

0 0 i

0 i 0

3 5

p (note that we must have

Do0 and B140 to reproduce the band structure). The emergent triply degenerate fermions around each crossing point are described by the linearized model obtained from equation (1) as

Ht q

2tB1 D1Kcqz iAqx iAqy

iAqx 2tD1Kcqz 0

iAqy 0 2tD1Kcqz

2

4

3

are three of the eight Gell-Mann matrices corresponding to spin-1 (ref. 37). This describes an isospin-1 triplet fermion moving in the xy plane. The fermion is helical, with a well-dened helicity of 1 and 0 corresponding to the eigenvalues of the helicity operator k l/k. The two branches with helicity 1 are massless

while the remaining one with helicity 0 has a at dispersion, and these are constrained by the presence of an emergent chiral symmetry operation C diag1; 1; 1 with H k

; C

f g 0.

The nite mass term for D 6 0 breaks this symmetry, generating

a mass for both 1 branches and gapping the spectrum.

Passing the critical point, the isospin-1 triplet fermion point split into two triply degenerate points located at kz tKc, as

shown in Fig. 3b, where t is an index labelling the two

crossing points and Kc

D=B1

5; 3

where the wave vector q is measured from the crossing point. One notes that in the direction perpendicular to the screw axis, that is, in the qxqy plane, the dispersion is the same as that for the isospin-1 triplet fermion in equation (2). One nds that the triply degenerate fermions basically inherit the structure of the triplet fermions but acquire one more degree of freedom, that is, the motion along the z direction. In the meantime, they lose the well-dened helicity.

4 NATURE COMMUNICATIONS | 8:15641 | DOI: 10.1038/ncomms15641 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15641 ARTICLE

a b c d

0.40.20.0 0.2

0.4

0.4

0.2

0.0

0.2

0.4

0.4

0.40.20.0 0.2

0.4

0.4

0.2

0.0

0.2

0.4

0.4

0.40.20.0 0.2 0.4

0.4

0.2

0.0

0.2

0.4

0.4

Energy (eV)

0.2

0.0

0.2

Z Z

kz

kz

kz

0.4

0.4

0.4

0.2

0.2

0.2

0.0

0.0

0.0

0.2

k x

0.2

k x

0.2

k x

ky

ky

ky

Figure 3 | Emergent fermions in Pentagon Carbon. Schematic of the band-crossing evolution in reciprocal space, from (a) a single isospin-1 triplet fermion point at the Brillouin zone centre, to (b) two triply degenerate fermion points along kz-axis (one on each side of G point), further to (c) two inter-connected (Hopf link) Weyl loops (with one Weyl loop being centered at G point, while the other loop crossing the rst Brillouin zone boundary). In a,b, the red line along the kz-axis marks the (twofold) band degeneracy line protected by the fourfold screw-rotational symmetry. Note that the dispersion transverse to this line (that is, in kxky plane) is of quadratic type (which can be observed from Fig. 2). (d) The DFT band structure of 3D Pentagon Carbon under broken fourfold screw-rotational symmetry (with uniaxial strain of 1% along z-axis followed by 1% tensile strain along x-axis). The red (blue) dot indicates the

crossing point of red (blue) circle in c with the kz-axis.

0.2

The two kinds of fermionic quasiparticles, the isospin-1 triplet fermions and the triply degenerate fermions, are new types of emergent fermions. Unlike the Dirac and Weyl fermions studied before, these new fermions do not have direct analogues in the relativistic quantum eld theory. The isospin-1 triplet fermion point is not topologically protected: it marks the critical point in the quantum phase transition. In contrast, the two triply degenerate fermion points are protected by the nontrivial band topology and the fourfold screw axis. These features are schematically illustrated in Fig. 3a,b.

When further breaking the screw-rotational symmetry, the double degeneracy of the E band will be lifted; hence the triply degenerate fermion point will disappear. However, note that the original A and E bands have different eigenvalues under the mirror operations Mxz and Myz. Hence crossings between each split E band and A band in the mirror-invariant kxkz and kykz planes will be protected if the respective mirror symmetry remains. For example, by an additional uniaxial strain along x direction or along diagonal directions in the xy plane, we nd that, interestingly, the resulting band crossings form two orthogonal concatenated Weyl loops, which has the topology of a Hopf link. As shown in Fig. 3c, for an additional tensile strain along the x direction, one Weyl loop is centred at the G point lying in the kxkz plane, while the other loop is in the kykz plane and crosses the rst Brillouin zone boundary. Around each point on the loop, the energy dispersion along two transverse directions is of the Weyl-like linear type as shown in Fig. 3d. It needs to be emphasized that the band crossing forms the Weyl-loop fermion, not Weyl fermion, because the inversion symmetry is still preserved. It is particularly interesting that the two Weyl loops form an inter-connected Hopf-link pattern, and these two loops are protected by the Mxz and Myz mirror symmetries, respectively.

It should also be noted that there is no symmetry to pin these Weyl loops at a constant energy, and, indeed, the energy varies along each loop.

Landau level spectrum. Topological character of the band structure can manifest when subjecting to a magnetic eld3840. For a 3D system, the electron motion in the plane perpendicular to the eld is quantized into LLs. The dispersion along the eld connects the levels to form dispersive Landau bands. For example, at each Weyl point, there is a gapless chiral Landau band, such that an additional parallel electric eld will drive charge creation or depletion depending on the chirality. Since Weyl points appear in pairs of opposite chirality, the process

leads to charge pumping between the Weyl points. This non-conservation of chiral charge at a single Weyl point is the well-known chiral anomaly14,18,40. Experimentally, it manifests as a negative magnetoresistance, which has recently been demonstrated in Weyl semimetals15,16.

For Pentagon Carbon, considering a magnetic eld along the screw axis such that the screw rotational symmetry is preserved. Focusing on the low-energy physics around the G point where the quantum phase transition occurs, we nd that the LL spectrum exhibits drastic difference across the transition. As illustrated in Fig. 4a, one observes that in the semiconductor phase, the LL spectrum shows a spectral gap, typical for a semiconducting state. Remarkably, in the TM phase, there appear two chiral Landau bands crossing the spectral gap and crossing the location of the original triply degenerate fermion points (Fig. 4b). This is like that of a Weyl point; however, there is a marked difference: the chiral band here is in fact with a huge number of Landau bands nearly degenerate onto a single curve21. This behaviour is due to the small values of the in-plane quadratic terms such as B2 and B3. If one articially increases these parameters, there will be a visible splitting between the almost degenerate Landau bands. Similar to chiral anomaly, a parallel E-eld would pump charges from one triply degenerate point to the other. Due to the large number of the bands, it is expected to generate pronounced signals in magneto-transport experiment.

0.2

a b c

Energy

0

3

0

kz

3 0

3 3 0

3 3 kz kz

Figure 4 | Landau level spectra. Landau level spectra for (a) the semiconductor phase, (b) the TM phase (corresponding to Fig. 2d) and (c) the TM phase with further screw-rotational symmetry breaking (corresponding to Fig. 3c). The magnetic eld is along the z-direction. The red circles in b indicate the location of the two triply degenerate fermion points (Fig. 3b).

NATURE COMMUNICATIONS | 8:15641 | DOI: 10.1038/ncomms15641 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15641

When further breaking the screw-rotational symmetry, the triply degenerate fermion points disappear, giving rise to inter-connected Hopf-link Weyl loops. In the LL spectrum, we see that the main effect of the symmetry breaking is to split the chiral Landau bands, as indicated in Fig. 4c. These interesting LL features reect the evolution of the band-crossing points and will be useful for probing the quantum phase transition and the emergent topological fermions.

DiscussionOur work discovers a new metastable 3D carbon allotropethe Pentagon Carbon, which exhibits a beautiful lattice structure entirely composed of pentagonal ringsa highly rare lattice motif for 3D crystals. Originating from this unique structure, there appear interesting quantum phase transitions and the associated novel emergent fermions. Remarkably, we reveal a fascinating transformation from isospin-1 triplet fermion point to triply degenerate fermion point and nally to inter-connected Hopf-link Weyl loops.

Nontrivial band-crossing point is also predicted in a few recent works. Particularly, Bradlyn et al.21 classied all possible degeneracy points at high-symmetry k-points for the 230 space groups in spinorbit-coupled systems, and triple-degenerate points have also been found in a few materials with strong spinorbit coupling22,41,42. The triply degenerate points here are distinct from those previous examples in the following aspects: rst, due to the negligible SOC, Pentagon Carbon is essentially a spinless system, and, as we mentioned at the beginning, spinful and spinless systems belong to different symmetry/topological classes (For example, the time reversal operation has distinct properties between the two cases. And rotations for spinful systems will also operate on the spin degree of freedom, requiring a double group representation, distinct from spinless systems.); second, while the degeneracy points in ref. 21 are at high-symmetry points, here they are on symmetry lines (GZ). Hence, without altering the symmetry, the latter is free to move along the symmetry line by applying a strain with observable physical consequences, whereas the former is pinned at the high-symmetry k-points. A work by Heikkila and Volovik discussed possible triply degenerate points in Bernal graphite27, where each point connects four Dirac band-crossing lines with each line protected by a Z2 invariant N1 1. In contrast, the

triply degenerate points in Pentagon Carbon are on a single-band degeneracy line along kz-axis (as in Fig. 3b); it is not a Dirac line because the dispersion is quadratic along the transverse directions, so the corresponding Z2 invariant, as dened in ref. 27, is trivial (N1 0). The protection of the band

degeneracy line and triply degenerate point here has a qualitatively different physical origin, which is originated from the fourfold screw-rotational symmetry and band ordering topology, not present in Bernal graphite. We note that the different congurations of degeneracy lines in Bernal graphite and Pentagon Carbon precisely reect their distinct symmetries: the arrangement of four Dirac lines with Z2 indices ( 1, 1, 1, 1) observed in Bernal graphite is allowed by its

threefold rotational symmetry, but it is not compatible with the fourfold screw-rotation in Pentagon Carbon.

Here, we focus on the effect of uniaxial strain in driving the phase transition. Other types of lattice deformation can also produce a similar effect, as long as the fourfold screw-rotational symmetry is preserved. For example, isotropic lattice strain or hydrostatic pressure can drive the phase transition with qualitatively similar features: the critical value is 1% for

isotropic strain and 4.75 GPa for hydrostatic pressure. We also mention that like other carbon allotropes, Pentagon Carbon has

excellent mechanical property. From rst-principles calculation, it shows a linear elastic regime up to 10% strain (Supplementary Fig. 5). Hence the strain-induced phase transitions discussed here should be readily achievable.

To guide experiment, we simulate X-ray diffraction peaks of Pentagon Carbon and contrast them with those of representative carbon allotropes (Supplementary Fig. 6). We also suggest a possible reaction route for the chemical synthesis. Molecules with pentagonal rings can be synthesized, for example, in cyclopentadiene with different substitutes. Starting from 1,4-dibromocyclopentadiene, the cyclopentadiene building units can be linked into a 1D ribbon through Yamamoto polymerization. The 1D ribbons could be functionalized with bromides and boronic acids which are connected into the 3D framework through Suzuki coupling (the details are described in Supplementary Fig. 7). Each step in this process is either already achieved or have a high chance to be achieved because of the existence of similar reactions. Given its energetic and dynamic stabilities and in view of the rapid progress in experimental techniques which have realized a number of carbon structures in recent decade34,43, we expect that Pentagon Carbon could also be realized in the near future.

Methods

Computational methods. We performed rst-principles calculations within the DFT formalism as implemented in VASP code44,45. The potential of the core electrons and the exchange-correlation interaction between the valence electrons are described, respectively, by the projector augmented wave46 approach and the generalized gradient approximation in the PerdewBurkeErnzerhof47 implementation. The kinetic energy cutoff of 500 eV is employed. The atomic positions were optimized using the conjugate gradient method. The energy and force convergence criteria are set to be 10 5 eV and 10 2 eV 1, respectively.

The band structure and projected density of states (PDOS) were calculated using the primitive cell. The Brillouin zone was sampled with a 7 7 7 Monkhorst

Pack special k-point grid for geometrical optimization. Phonon spectra are calculated using force-constants method, and the dynamical matrices are computed using the nite differences method in a 2 2 2 supercell to eliminate errors in

the low frequency modes (calculation with 3 3 3 supercell is also carried out for

verication). The force constants and phonon spectra were obtained using the Phonopy package48. Calculation with the HeydScuseriaErnzerhof (HSE06) hybrid functional49 is also performed to verify the band structure. The values of the band gap and the critical transition strain become larger compared with Perdew BurkeErnzerhof results (the HSE06 gap is about 471 meV and the critical strain is about 2%), but the qualitative features of the bands and phase transitions remain

the same (Supplementary Fig. 8 and Note 3). The effect of van der Waals interaction is also tested and is found to be weak, for example, by including the van der Waals interaction, the cohesive energy is decreased from 8.79 to 8.89 eV

and the lattice constants are decreased by about 0.2% (Supplementary Fig. 9).

Data availability. The data that support the ndings of this study are available from the corresponding authors upon reasonable request.

References

1. Iijima, S. Helical microtubules of graphitic carbon. Nature 354, 5658 (1991).2. Geim, A. K. & Novoselov, K. S. The rise of graphene. Nat. Mater. 6, 183191 (2007).

3. Hirsch, A. The era of carbon allotropes. Nat. Mater. 9, 868871 (2010).4. Kroto, H. W., Heath, J. R., OBrien, S. C., Curl, R. F. & Smalley, R. E. C60: Buckminsterfullerene. Nature 318, 162163 (1985).

5. Ewels, C. P. et al. Predicting experimentally stable allotropes: instability of penta-graphene. Proc. Natl Acad. Sci. USA 112, 1560915612 (2015).

6. Dresselhaus, M. S., Dresselhaus, G. & Jorio, A. Group Theory: Application to the Physics of Condensed Matter (Springer Science & Business Media, 2007).

7. Wan, X., Turner, A. M., Vishwanath, A. & Savrasov, S. Y. Topological semimetal and Fermi-arc surface states in the electronic structure of pyrochlore iridates. Phys. Rev. B 83, 205101 (2011).

8. Soluyanov, A. A. et al. Type-II Weyl semimetals. Nature 527, 495498 (2015).9. Xu, S. Y. et al. Discovery of a Weyl fermion semimetal and topological Fermi arcs. Science 349, 613617 (2015).

10. Turner, A. M., Vishwanath, A. & Head, C. O. Beyond band insulators: topology of semimetals and interacting phases. Topol. Insul. 6, 293324 (2013).

6 NATURE COMMUNICATIONS | 8:15641 | DOI: 10.1038/ncomms15641 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15641 ARTICLE

11. Hasan, M., Xu, S. & Neupane, M. Topological insulators, topological crystalline insulators, topological Kondo insulators, and topological semimetals. Preprint at https://arxiv.org/abs/1406.1040

Web End =https://arxiv.org/abs/1406.1040 (2014).

12. Yang, S. A., Pan, H. & Zhang, F. Dirac and Weyl superconductors in three dimensions. Phys. Rev. Lett. 113, 046401 (2014).

13. Yang, S. A., Pan, H. & Zhang, F. Chirality-dependent Hall effect in Weyl semimetals. Phys. Rev. Lett. 115, 156603 (2015).

14. Aji, V. Adler-Bell-Jackiw anomaly in Weyl semimetals: application to pyrochlore iridates. Phys. Rev. B 85, 241101 (2012).

15. Huang, X. et al. Observation of the chiral-anomaly-induced negative magnetoresistance in 3D Weyl semimetal TaAs. Phys. Rev. X 5, 031023 (2015).

16. Parameswaran, S. A., Grover, T., Abanin, D. A., Pesin, D. A. & Vishwanath, A. Probing the chiral anomaly with nonlocal transport in three-dimensional topological semimetals. Phys. Rev. X 4, 031035 (2014).

17. Nielsen, H. B. & Ninomiya, M. The Adler-Bell-Jackiw anomaly and Weyl fermions in a crystal. Phys. Lett. B 130, 389396 (1983).

18. Son, D. T. & Spivak, B. Z. Chiral anomaly and classical negative magnetoresistance of Weyl metals. Phys. Rev. B 88, 104412 (2013).

19. Young, S. M. et al. Dirac semimetal in three dimensions. Phys. Rev. Lett. 108, 140405 (2012).

20. Yang, S. A. Dirac and Weyl materials: fundamental aspects and some spintronics applications. SPIN 06, 1640003 (2016).

21. Bradlyn, B. et al. Beyond Dirac and Weyl fermions: unconventional quasiparticles in conventional crystals. Science 353, aaf5037 (2016).

22. Chang, G. et al. New fermions on the line in topological symmorphic metals. Preprint at https://arxiv.org/abs/1605.06831

Web End =https://arxiv.org/abs/1605.06831 (2016).

23. Orlita, M. et al. Observation of three-dimensional massless Kane fermions in a zinc-blende crystal. Nat. Phys. 10, 233238 (2014).

24. Dai, X. Weyl semimetals: a group family picture. Nat. Mater. 15, 56 (2015).25. Liu, Z. K. et al. Evolution of the Fermi surface of Weyl semimetals in the transition metal pnictide family. Nat. Mater. 15, 2731 (2016).

26. Novoselov, K. S. et al. Two-dimensional gas of massless Dirac fermions in graphene. Nature 438, 197200 (2005).

27. Heikkila, T. T. & Volovik, G. E. Nexus and Dirac lines in topological materials. New J. Phys. 17, 093019 (2015).

28. Zhong, C. et al. Towards three-dimensional Weyl-surface semimetals in graphene networks. Nanoscale 8, 72327239 (2016).

29. Chen, Y. et al. Nanostructured carbon allotropes with Weyl-like loops and points. Nano Lett. 15, 69746978 (2015).

30. Weng, H. et al. Topological node-line semimetal in three-dimensional graphene networks. Phys. Rev. B 92, 045108 (2015).

31. Wang, J. T. et al. Body-centered orthorhombic C_{16}: a novel topological node-line semimetal. Phys. Rev. Lett. 116, 195501 (2016).

32. Cheng, Y. et al. Body-centered tetragonal C16: a novel topological node-line semimetallic carbon composed of tetrarings. Small 13, 1602894 (2017).33. Zhang, S., Wang, Q., Chen, X. & Jena, P. Stable three-dimensional metallic carbon with interlocking hexagons. Proc. Natl Acad. Sci. USA 110, 1880918813 (2013).

34. Johnston, R. L. & Hoffmann, R. Superdense carbon, C8: supercubane or analog of .gamma.-silicon? J. Am. Chem. Soc. 111, 810819 (1989).

35. Chen, Y. et al. Carbon kagome lattice and orbital-frustration-induced metal-insulator transition for optoelectronics. Phys. Rev. Lett. 113, 085501 (2014).

36. Zhong, C., Xie, Y., Chen, Y. & Zhang, S. Coexistence of at bands and Dirac bands in a carbon-Kagome-lattice family. Carbon 99, 6570 (2016).

37. Georgi, H. & Jagannathan, K. Lie algebras in particle physics. Am. J. Phys. 50, 10531053 (1982).

38. Jeon, S. et al. Landau quantization and quasiparticle interference in the three-dimensional Dirac semimetal Cd3As2. Nat. Mater. 13, 851856 (2014).

39. Hosur, P. & Qi, X. Recent developments in transport phenomena in Weyl semimetals. C. R. Phys. 14, 857870 (2013).

40. Yu, Z. M., Yao, Y. & Yang, S. A. Predicted unusual magnetoresponse in type-II Weyl semimetals. Phys. Rev. Lett. 117, 077202 (2016).

41. Zhu, Z., Winkler, G. W., Wu, Q., Li, J. & Soluyanov, A. A. Triple point topological metals. Phys. Rev. X 6, 031003 (2016).

42. Weng, H., Fang, C., Fang, Z. & Dai, X. Topological semimetals with triply degenerate nodal points in y-phase tantalum nitride. Phys. Rev. B 93, 241202 (2016).

43. Liu, P., Cui, H. & Yang, G. W. Synthesis of body-centered cubic carbon nanocrystals. Cryst. Growth Des. 8, 581586 (2008).

44. Kresse, G. & Furthmller, J. Efciency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. Mater. Sci. 6, 1550 (1996).

45. Kresse, G. & Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. Rev. B 47, 558561 (1993).

46. Blchl, P. E. Projector augmented-wave method. Phys. Rev. B 50, 1795317979 (1994).

47. Perdew, J. P. et al. Atoms, molecules, solids, and surfaces: applications of the generalized gradient approximation for exchange and correlation. Phys. Rev. B 46, 66716687 (1992).

48. Togo, A., Oba, F. & Tanaka, I. First-principles calculations of the ferroelastic transition between rutile-type and CaCl2-type SiO2 at high pressures. Phys. Rev.

B 78, 134106 (2008).49. Heyd, J., Scuseria, G. E. & Ernzerhof, M. Hybrid functionals based on a screened Coulomb potential. J. Chem. Phys. 118, 8207 (2003).

Acknowledgements

We thank the helpful discussions with Man-Fung Cheung, D.L. Deng, Shuai Yuan, Weiyu Xie and Damien J. West. This work was supported by the National Natural Science Foundation of China (Nos. 51376005 and 11474243). S.Z. acknowledges the support by US DOE under Grant No. DE-SC0002623. Z.-M.Y. and S.A.Y. are supported by the Singapore MOE Academic Research Fund Tier 1 (SUTD-T1-2015004) and Tier 2 (MOE2015-T2-2-144).

Author contributions

Y.C. proposed the Pentagon Carbon. Y.C., S.A.Y. and S.Z. conceived the original ideas. C.Z., Y.C., Z.-M.Y. and S.A.Y. conducted the calculations. Y.C. and S.A.Y. wrote the manuscript. All authors discussed the results and commented on the manuscript at all stages.

Additional information

Supplementary Information accompanies this paper at http://www.nature.com/naturecommunications

Web End =http://www.nature.com/ http://www.nature.com/naturecommunications

Web End =naturecommunications

Competing interests: The authors declare no competing nancial interests.

Reprints and permission information is available online at http://npg.nature.com/reprintsandpermissions/

Web End =http://npg.nature.com/ http://npg.nature.com/reprintsandpermissions/

Web End =reprintsandpermissions/

How to cite this article: Zhong, C. et al. Three-dimensional Pentagon Carbon with a genesis of emergent fermions. Nat. Commun. 8, 15641 doi: 10.1038/ncomms15641 (2017).

Publishers note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional afliations.

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articles Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

Web End =http://creativecommons.org/licenses/by/4.0/

r The Author(s) 2017

NATURE COMMUNICATIONS | 8:15641 | DOI: 10.1038/ncomms15641 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 7