TY - JOUR AU - Liu, Zhirong AB - Abstract Inspired by the great development of graphene, more and more research has been conducted to seek new two-dimensional (2D) materials with Dirac cones. Although 2D Dirac materials possess many novel properties and physics, they are rare compared with the numerous 2D materials. To provide explanation for the rarity of 2D Dirac materials as well as clues in searching for new Dirac systems, here we review the recent theoretical aspects of various 2D Dirac materials, including graphene, silicene, germanene, graphynes, several boron and carbon sheets, transition-metal oxides (VO2)n/(TiO2)m and (CrO2)n/(TiO2)m, organic and organometallic crystals, so-MoS2, and artificial lattices (electron gases and ultracold atoms). Their structural and electronic properties are summarized. We also investigate how Dirac points emerge, move, and merge in these systems. The von Neumann–Wigner theorem is used to explain the scarcity of Dirac cones in 2D systems, which leads to rigorous requirements on the symmetry, parameters, Fermi level, and band overlap of materials to achieve Dirac cones. Connections between existence of Dirac cones and the structural features are also discussed. graphene, silicene, graphyne, two-dimensional, Dirac cone, Dirac materials INTRODUCTION Two-dimensional (2D) crystal was thought to be unstable in nature for a long time until 2004 when graphene, a one-atom-thick honeycomb structure composed of carbon atoms, was successfully prepared [1]. The discovery of graphene has attracted great interest because of the promising prospects of graphene in both basic and applied research [2,3]. In particular, the Dirac-cone structure gives graphene massless fermions, leading to half-integer [4,5]/fractional [6,7]/fractal [8–10] quantum Hall effects (QHE), ultrahigh carrier mobility [11], and many other novel phenomena and properties [12,13]. The great development of graphene also opens a door for discovering more 2D materials [14,15]. Up to now, over hundreds of 2D materials have been found, including group-IV compounds, binary systems of group III-V elements, metal chalcogenides, complex oxides, and so on [14,15]. But among them, only graphene [4,5,16], silicene and germanene (graphene-like silicon and germanium, respectively) [17], several graphynes (sp-sp2 carbon allotropes) [18,19], and some other systems [20–31] have been predicted to be Dirac materials (summarized in Fig. 1). Moreover, only Dirac cones in graphene have been truly confirmed experimentally [4,5]. Why are the 2D Dirac materials so rare? What are the special structure and properties for them? And how can we find new massless Dirac systems? Figure 1. Open in new tabDownload slide 2D materials with Dirac cones, which contain various elementary substances (Es) and compounds (Cp) with the elements of B, C, O, Si, Ge, Mo, V, Pb, etc. In detail, there are Pmmn boron, graphene, α-graphyne, 6,6,12-graphyne, 14,14,18-graphyne, S-graphene, square carbon SG-10b, silicene, germanene, CO on Cu(111), (VO2)3/(TiO2)5, so-MoS2, and Pb2(C6H4)3 in the clockwise direction. Figure 1. Open in new tabDownload slide 2D materials with Dirac cones, which contain various elementary substances (Es) and compounds (Cp) with the elements of B, C, O, Si, Ge, Mo, V, Pb, etc. In detail, there are Pmmn boron, graphene, α-graphyne, 6,6,12-graphyne, 14,14,18-graphyne, S-graphene, square carbon SG-10b, silicene, germanene, CO on Cu(111), (VO2)3/(TiO2)5, so-MoS2, and Pb2(C6H4)3 in the clockwise direction. This paper reviews the theoretical studies on the 2D materials with Dirac cones. We start with a survey on the known 2D Dirac materials up to date, where emphases are put on their structural and electronic properties. Then, we discuss the underlying mechanism for the rarity of 2D Dirac systems, i.e. how Dirac cones emerge and merge in a 2D crystal, and what conditions are required for a 2D Dirac system. THE KNOWN 2D DIRAC MATERIALS Graphene Graphene is the most typical 2D material with Dirac cones [12]. It can be regarded as a single layer of graphite, which has two C atoms per unit cell arranged in a hexagonal lattice (Fig. 2a). The C−C bond length in graphene is 1.42 Å, and the lattice constant is 2.46 Å. Graphene is rather stable because the C atoms are sp2-hybridized and bind together with both σ and π bonds. Figure 2. Open in new tabDownload slide (a) Atomic structure and (b) energy bands of graphene. The energy bands were calculated by first-principles [32] (black line) and TB model with Equation (2) (green line) and Equation (3) (red line). Figure 2. Open in new tabDownload slide (a) Atomic structure and (b) energy bands of graphene. The energy bands were calculated by first-principles [32] (black line) and TB model with Equation (2) (green line) and Equation (3) (red line). The electronic structure of graphene can be described by a tight-binding (TB) approach [12,16]. Under nearest-neighbor approximation, the Hamiltonian is simplified into a 2 × 2 matrix for each wave vector k: \begin{equation} H\!\left( {\bf k} \right) = \left[ {\begin{array}{cc} {\varepsilon _0 } &\quad { - \sum\limits_{i = 1,2,3} {t_i \exp \left( {i{\bf k} \cdot {\bf d}_i } \right)} } \\ {{\rm c}{\rm .c}{\rm .}} &\quad {\varepsilon _0 } \end{array}} \right], \end{equation} (1) where di (i = 1, 2, 3) are vectors that connect a C atom to its three nearest neighbors, ti are the corresponding hopping energies, and ε0 is the on-site energy. The energy bands are thus solved to be \begin{equation} E_ \pm \!\left( {\bf k} \right) = \pm \left| {\sum\limits_{i = 1,2,3} {t_i \exp \left( {i{\bf k} \cdot {\bf d}_i } \right)} } \right| \end{equation} (2) by setting the Feimi level to be ε0 = 0. For the equilibrium structure, ti ≡ t0 (≈ 2.7 eV) and di ≡ r0 (≈ 1.42 Å), and the valence and conduction bands contact at K and K′ points of the hexagonal Brillouin zone. Expanding the energy bands around K (or K′) gives \begin{equation} E_ \pm \left( {\bf q} \right) = \pm \hbar \upsilon _{\rm F} \left| {\bf q} \right|, \end{equation} (3) where k = K + q, and υF = 3t0r0/2|${\hbar}$| is the fermi velocity (∼106 m/s). Equation (3) shows that graphene has a cone-like band structure with linear dispersion near K (or K′) points, similar to a relativistic particle. The TB results are consistent with first-principles calculations at low-energy regions (Fig. 2b) [32]. The density of states (DOS) per unit cell (with a degeneracy of 4 included) near Fermi level is expressed as \begin{equation} \rho (E) = \frac{{4\left| E \right|}}{{\sqrt 3 \pi t_0 ^2 }}.\end{equation} (4) Thus, graphene is a gapless semiconductor with zero DOS at Fermi level. The Hamiltonian near K (or K′) points can be also transformed [12,33] into \begin{equation} H = \upsilon _{\rm F} {\bf p} \cdot {\boldsymbol \sigma }, \end{equation} (5) where p = −iℏ∇ is the momentum operator and σ is the Pauli matrice. Equation (5) is identical to the massless Dirac equation (or Dirac–Weyl equation with spin S = 1/2) by replacing c (speed of light) with υF. Therefore, the K (or K′) points are also called Dirac points, and the linear band structure is named as Dirac cone. Graphene possess massless Dirac fermions with pseudospins of ±1/2. Many novel physical phenomena and electronic properties have been found in graphene due to the Dirac-cone structure. For example, a particular Landau level forms in graphene when a uniform perpendicular magnetic field B is applied [5,34], \begin{equation} E(N) = {\mathop{\rm sgn}} (N)\sqrt {2e\hbar \upsilon _{\rm F} ^2 \left| N \right|B} \end{equation} (6) while E(N)∝N in normal semiconductors. The Landau levels in 2D systems can be verified by Shubnikov-de Haas oscillations or QHE (if the magnetic field is strong enough). Different from the conventional 2D systems, graphene shows a half-integer QHE [4,5] with \begin{equation} {\boldsymbol \sigma }_{{\rm x}y} {\rm = }(4e^2 {/}h)(N + 1/2),\end{equation} (7) which can be explained by the pseudospins and degeneracy of Dirac fermions. When interactions between electrons and magnetic flux quanta become very strong, new quasiparticles with fractional electronic charge can form and lead to correlated electron phenomena. Fractional QHE has been observed in suspended graphene [6,7] and graphene on hexagonal BN (h-BN) [35] by reducing external perturbations, and the conductance G can be expressed as \begin{equation} G = ve^2 /h,\end{equation} (8) where v = p/q is the filling factor with p and q integers. Phase transitions between different fractional quantum Hall states have also been observed in suspended graphene, suggesting changes in spin and/or valley polarization [36]. In 1976, Hofstadter predicted a recursive electronic structure for 2D-confined electrons in magnetic field, which is called Hofstadter's butterfly [37]. This physical picture has been verified by the observation of fractal QHE in monolayer [9,10] and bilayer [8] graphene on h-BN when the magnetic length is comparable to the size of the superlattice. Besides for the various QHEs, ultrahigh carrier mobility has also been found in graphene due to the massless Dirac-cone structure [13]. For graphene on Si/SiO2 substrate, charge impurities are main scattering source which limit the mobility to be the order of 104 cm2 V−1 s−1 at low temperature [38]. The mobility of graphene on h-BN can be ∼6 × 104 cm2 V−1 s−1 which is three times larger than that on Si/SiO2, because the former substrate is flatter and has less charge impurities [39]. If all extrinsic scatterings are excluded, graphene has intrinsic mobility of 2 × 105 cm2 V−1 s−1 at room temperature [40,41]. Castro et al. [42] have studied the scattering mechanism of suspended graphene at different temperatures. They found that flexural phonons limit the intrinsic mobility for temperature higher than 10 K. Various electron–phonon couplings have been investigated in a recent work which emphasizes the contribution of high-energy, optical, and zone-boundary phonons to electrical resistivity at room temperature [43]. Li et al. [44] has also revealed that both longitudinal acoustic (LA) and transverse acoustic (TA) phonons are important in determining the intrinsic mobility of graphene. The electronic properties of graphene are affected by many factors in reality, like ripples and substrates. Suspended graphene shows intrinsic ripples with height of ∼1 nm and size of ∼10 nm due to the thermal fluctuation [45,46]. Mid-gap states, non-zero DOS at fermi energy, and charge inhomogenity (electron-hole puddles) can form due to the ripple-induced modulation of hopping terms and localized states [47]. Graphene on a substrate also exhibits corrugations because of the substrate roughness and the lattice mismatch. For example, the typical corrugation is 0.2 nm in height for graphene on SiO2 but only 0.02 nm for graphene on mica, which is in agreement with the substrate morphology [48]. For graphene on substrates, the charge inhomogenity mainly arises from substrate-induced charge impurites [49] and can be improved by using h-BN instead of SiO2 [50]. The interaction between graphene and substrates also contains strain effect, charge transfer, orbital hybridization, etc. A large uniaxial strain up to 24% might open the bandgap of graphene, but the Dirac cones keep robust for small and moderate uniform deformations [51]. The uniaxial and shear strains cause a moving of Dirac points [51,52] and anisotropic Fermi surface for graphene [51,53]. Pereira and Castro Neto have proposed strain engineering of graphene by patterning the substrates [54]. The local-strain-induced gauge fields can be tailored to generate tunnel suppression, transport gap, and electron confinement in graphene [54]. Especially for a local triangular strain in graphene, the gauge field acts as a pesudomagnetic field which satisfies a simple model derived by Guinea et al. [55], \begin{equation} B_{\rm s} = 8\beta c/a,\end{equation} (9) where a is the lattice constant, β = −∂ln t/∂a ≈ 2 (t is hopping energy), c = Δm/D is the largest strain Δm for a disc of diameter D. When Δm = 10% and D = 100 nm, Bs is estimated to be ∼40 T and the largest Landau gap ∼0.25 eV, resulting an observable QHE [55]. Landau levels have been detected in graphene nanobubbles with a strain-induced pesudomagnetic field greater than 300 T [56]. The graphene-metal interfaces can be divided into two classes according to the binding energies Eb [57]. Graphene on Al, Ag, Cu, Au, and Pt surfaces which has weak adsorption (physisorption, Eb < 0.5 eV) preserves the Dirac-cone structure and shows n- or p-doping because of charge transfer. However, graphene interacts strongly with Co, Ni, Pd, and Ti (chemsorption, Eb > 0.8 eV) by p-d hybridization, leading to an opened band gap of graphene submerged under the conducting states from metals [57]. Silicene and germanene The great success on graphene has inspired people to search for other 2D Dirac materials. Since carbon belongs to group IVA, silicon and germanium in the same group attract more and more attention. The critical problems are whether graphene-like silicon (silicene) and germanium (germanene) exist and what their electronic properties are. Using first-principles calculations, Cahangirov et al. [17] have predicted that silicene and germanene have stable low-buckled honeycomb structures (Fig. 3a and b). The bond lengths in silicene and germanene are 2.25 and 2.38 Å, respectively, which are much longer than the C−C length in graphene. The enlarged bond lengths weaken the π−π interaction and cause distinct coupling of σ and π bonds to form buckled structures. The buckling distance is 0.44 Å for silicene and 0.64 Å for germanene. Both low-buckled silicene and germanene are semimetals with Dirac cones at K points (Fig. 3c and d), and the Fermi velocity was estimated to be ∼106 m/s, very close to that of graphene. Figure 3. Open in new tabDownload slide Atomic and electronic structures of (a, c) silicene and (b, d) germanene [17] (Copyright 2009 The American Physical Society). Figure 3. Open in new tabDownload slide Atomic and electronic structures of (a, c) silicene and (b, d) germanene [17] (Copyright 2009 The American Physical Society). However, the buckled structures greatly affect the electronic properties of silicene and germanene and bring about new physics beyond graphene. Liu et al. [58] have reported that the spin–orbit coupling (SOC) increases with the buckling degree and opens a band gap of 1.55 meV in silicene and 23.9 meV in germanene at Dirac points (Fig. 4a). The energy dispersion near Dirac points is written as \begin{equation} E({\bf k}) = \pm \sqrt {(v_{\rm F} \hbar {\bf k})^2 + \Delta _{{\rm so}} ^{\rm 2} } , \end{equation} (10) where Δso is the effective SOC. The spin–orbit band gap and non-trivial topological properties might result in detectable quantum spin Hall effect. Considering the buckled structures, Ni et al. [59] have proposed to open the band gap of silicene or germanene by applying a vertical electric field Ez to break the atomic symmetry (Fig. 4b). Drummond et al. [60] have considered the effect of SOC and vertical electric field simultaneously and have derived \begin{equation} E_{{\boldsymbol \sigma }\xi , \pm } \left( {\bf k} \right) = \pm \sqrt {\left( {\hbar \upsilon _{\rm F} {\bf k}} \right)^2 + \left( {\Delta _{{\rm so}} + {\boldsymbol \sigma }\xi \Delta _{\rm z} } \right)^2 } , \end{equation} (11) where ξ = ±1 distinguishes K and K′, σ = ±1 distinguishes spin-up and -down, and Δz is a half of the electric-field-induced band gap. Thus, silicene becomes a semimetal at a critical field Ec ≈ 20 mV Å−1 when Δso = ±Δz. And it undergoes a transition from a topological insulator (TI) when |Ez| < Ec to a simple band insulator (BI) when |Ez| > Ec. Similar results (Fig. 4c) have been obtained by Ezawa with further consideration of Rashba SOC, and a coexistence of TI, BI, and metal regions in a silicene sheet by application of an inhomogeneous electric field [61]. Ezawa also depicted two phase diagrams of silicene with exchange field M and electric field Ez and with light A and electric field Ez, which include various insulating, semimetal, and metallic states (Fig. 4d and e) [62,63]. Recently, a new quantum state, valley-polarized quantum anomalous Hall state (Fig. 4f) has been predicted in silicene through tuning the Rashba SOC [64]. The intrinsic carrier mobilities of silicene, germanene, and graphene have the same order of magnitude under deformation potential theory (only the electron–acoustic phonon interaction is considered) [65,66], but the electron–optical phonon coupling in silicene or germanene is ∼25 times smaller than that in graphene [67]. Based on the state transition of silicene under electric fields, a silicene-based spin filter has been proposed and demonstrated to exhibit 98% spin polarization [68], and a thermal-induced pure valley and spin currents have also been found in silicene junctions [69]. Figure 4. Open in new tabDownload slide (a) The evolution of SOC-induced gap with buckling for silicene [58] (Copyright 2011 American Physical Society). (b) The band gap of silicene as a function of vertical electric field [59] (Copyright 2011 American Chemical Society). (c) The band structures of a silicene nanoribbon when Ez = 0 and Ez = 2Ec [61]. (d) Topological phase diagram of silicene in Ez − M plane [62] (Copyright 2012 American Physical Society). (e) Topological phase diagram of silicene in Ez − A plane [63] (Copyright 2013 American Physical Society). (f) The electronic structure of a zigzag-terminated silicene with the valley-polarized QAHE [64] (Copyright 2014 American Physical Society). Figure 4. Open in new tabDownload slide (a) The evolution of SOC-induced gap with buckling for silicene [58] (Copyright 2011 American Physical Society). (b) The band gap of silicene as a function of vertical electric field [59] (Copyright 2011 American Chemical Society). (c) The band structures of a silicene nanoribbon when Ez = 0 and Ez = 2Ec [61]. (d) Topological phase diagram of silicene in Ez − M plane [62] (Copyright 2012 American Physical Society). (e) Topological phase diagram of silicene in Ez − A plane [63] (Copyright 2013 American Physical Society). (f) The electronic structure of a zigzag-terminated silicene with the valley-polarized QAHE [64] (Copyright 2014 American Physical Society). The predicted novel properties make silicene and germanene promising in electronics, but can silicene or germanene be obtained experimentally? Lalmi et al. [70] have reported an epitaxial growth of a silicene sheet on Ag(111). The scanning tunneling microscope (STM) and low energy electron diffraction (LEED) images show a (2|${\sqrt {3}}$| × 2|${\sqrt {3}}$|⁠) R30° superstructure (with respect to Ag), and the Si−Si distance is 0.18–0.19 nm. Later studies revealed that silicene on Ag(111) has various superstructures [71–74], i.e. (⁠|${\sqrt{3}}$| × |${\sqrt {3}}$|⁠)R30°-Si (with respect to Si), (4 × 4), (⁠|${\sqrt {13}}$| × |${\sqrt {13}}$|⁠) R13.9°, and (⁠|${\sqrt {7}}$| × |${\sqrt {7}}$|⁠) R19.1° phases, etc. The superstructures are sensitive to the silicon coverage and substrate temperature [71,72], and the Si−Si distance varies between 2.28 and 2.5 Å [74]. Although several groups declared that they observed Dirac-cone structure of silicene on Ag(111) [75,76], the conclusion is controversial [77,78], especially after Lin et al. [78] showed that the Dirac cones disappear for (4 × 4) silicene on Ag(111). Lin et al. [78] suggested that the linear dispersion of silicene on Ag(111) measured by Angular-resolved photoelectron spectroscopy [76] originate from the Ag bulk sp-band. More and more works support the absence of Dirac cone in silicene on Ag(111) due to the strong interaction between Si and Ag atoms and symmetry breaking [79–83]. Silicene has also been successfully grown on ZrB2(0001) [84], Ir(111) [85], and MoS2 [86]. Silicene on ZrB2(0001) is a (⁠|${\sqrt {3}}$| × |${\sqrt {3}}$|⁠)-reconstructed Si honeycomb sheet containing three kinds of Si atoms, which has a direct band gap [84]. For silicene on Ir(111), a buckled (⁠|${\sqrt {3}}$| × |${\sqrt {3}}$|⁠)-Si sheet with an undulation is formed [85]. Quhe et al. [87] predicted that the absence of Dirac cones for silicene on various metal substrates is common because of the strong band hybridizations. But the Dirac cones might be restored by intercalating alkali metal atoms between silicene and substrates [87]. Recently, a highly buckled silicene nanosheet has been reported to form on MoS2 [86]. First-principles calculations have shown that silicenes on MoS2, MoSe2, and GaTe are metallic while those on MoTe2, GaS, and GaSe are gapless semiconductors [88]. Theorists have also proposed several substrates with weak silicene–substrate interaction to reserve the Dirac fermions, like h-BN and hydrogen-passivated Si(111) [80,81,89]. Up to now, no freestanding silicene has been obtained, and the Dirac cone of silicene has not been confirmed. A new substrate is needed for both growing silicene and preserving the Dirac cones. Since germanene has similar electronic properties to silicene, the growth of germanene is also attractive. In 2014, Li et al. [90] first reported the fabrication of germanene sheet on Pt(111) which has a distorted (⁠|${\sqrt {19}}$| × |${\sqrt {19}}$|⁠) superstructure with respect to the substrate. Calculations have shown that germanene on Ag(111)/h-BN has an adsorption energy of −464∼−428/−130 meV per Ge atom and is metallic/semiconducting [91]. Germanene on MoS2 has also been discussed theoretically which is a p-doped semiconductor with a band gap of 24 meV [92]. More honeycomb sheet of group IVA elements and III–V binary compounds have been investigated theoretically to give semiconducting properties except SiGe [93,94]. Stanene which is a honeycomb sheet of tin atoms has been predicted to be a sizable-gap QSH insulator [95] but might be too unstable to exist in reality. Graphynes Graphynes are a series of 2D carbon allotropes composed of sp and sp2 hybridized atoms, which were proposed by Baughman et al. in 1987 [96]. Graphynes have various structures based on the ratio and arrangement of sp and sp2 C atoms and possess high thermal stability [96]. The C–C bonds in most graphynes fall into three types: the triple bond (C≡C) with a length of 1.21–1.24 Å, the double bond (C=C) with a length of 1.31–1.35 Å, and the conjugated double bond (C=C) with a length of 1.39–1.45 Å [97,98]. Early first-principles calculations showed that graphynes can be semiconductors or semimetals or metals [97–99]. In 2012, Malko et al. [18] reported that α-, β-, and 6,6,12-graphyne (Fig. 5a–c) could have Dirac-cone structures (Fig. 5d–f). Both α-graphyne and β-graphyne exhibit hexagonal symmetry, but the Dirac points in the former are at K points while they locate along Γ-M line in the latter. In addition, 6,6,12-graphyne has a rectangular lattice and possesses two kinds of distorted Dirac cones with self-doping. These results not only indicate that graphynes might have more versatile properties than graphene but also break the stereotype that a honeycomb or hexagonal structure is a prerequisite for Dirac cones. Figure 5. Open in new tabDownload slide Atomic and band structures of (a, d) α-graphyne, (b, e) β-graphyne, and (c, f) 6,6,12-graphyne [18] (Copyright 2012 American Physical Society). Figure 5. Open in new tabDownload slide Atomic and band structures of (a, d) α-graphyne, (b, e) β-graphyne, and (c, f) 6,6,12-graphyne [18] (Copyright 2012 American Physical Society). Why are some graphynes semimetals with Dirac cones while others are semiconductors? Kim et al. [100] have first noticed that α-, β-, and γ graphyne are topologically equivalent to graphene where the role of triple bonds can be simplified as effective hopping terms. Liu et al. [101] have also used a TB model with effective hopping terms to fit the band structures of a few graphynes very well. Huang et al. [19] then derived a criterion for the existence of Dirac cones in graphynes, \begin{equation} \frac{{t_r }}{{t_b }} = 1\;{\rm or}\; - 2 \le \frac{{t_r }}{{t_b }} \le - 1,\end{equation} (12) where tr and tb are the two kinds of effective hopping energies in a graphene-like TB model to describe the band structures of graphynes. Equation (12) indicates that whether Dirac cones exist or not in graphynes is determined by the combination of hopping energies. With tr and tb values derived from first-principles calculations, this criterion successfully explains the properties of various graphynes, including α-, β-, γ- and 6,6,12-graphyne. It also suggests that more graphynes with Dirac cones can emerge by tailoring the hopping terms. Actually, they have predicted that two graphynes (14,14,14-graphyne and 14,14,18-graphyne) without hexagonal symmetry have Dirac cones [19]. Three other rectangular graphynes have recently been designed, and two of them (6,6,18-graphyne and h-12,12,20-graphyne) have anisotropic Dirac cones along Γ-X′ lines [102]. Hexagonal δ-graphyne which is more stable than α- or β-graphyne shows Dirac-cone structure, too [103]. More importantly, δ-graphyne becomes a TI under SOC with an induced gap of 0.59 meV and a Z2 topological invariant of v = 1 [103]. A TB approach to investigate SOC effect in more graphynes has been developed, and rich phase transitions were predicted [104,105]. Besides designing new graphynes, another way to construct Dirac systems is tuning the structure of known graphynes. Cui et al. [106] have observed a transition from semiconductor to semimetal with Dirac cones in graphdiyne under uniaxial tensile strain. The reverse way is also possible, e.g. a transition from gapless to finite gap system has also been achieved by an application of tensile strain to 6,6,12-graphyne [107]. The room-temperature intrinsic carrier mobilities of α-, β-, and 6,6,12-graphyne calculated under deformation potential theory vary from 0.8 × 104 to 5.4 × 105 cm2 V−1 s−1, which can be even higher than that of graphene [108]. Further consideration of various electron–phonon interactions has indicated that the mobilities of α- and γ graphynes are limited by LA phonon scatterings when temperatures are lower than 600 K and are ∼104 cm2 V−1 s−1 at room temperature [109]. Recently, an analytic formula for the intrinsic mobility of 2D Dirac systems has been derived by Li et al. [44] considering the influences of LA and TA phonon modes and anisotropy. It reveals that the high mobility of graphynes originates from a suppression of TA phonon scatterings. Directional transport properties of 6,6,12-graphyne have been predicted, which can be further manipulated by strain [110]. Although novel electronic properties have been proposed in graphynes, the experimental progress is still in its infancy [111]. Among the various graphynes, only graphdiyne-like films have been reported to be synthesized by homocoupling of monomers [112]. The related reaction mechanisms have been investigated both experimentally and theoretically [113–115], and the main difficulty to grow graphynes is how to reduce the side reactions. Rectangular carbon and boron allotropes Carbon is the element with the most allotropes. More 2D carbon lattices with Dirac cones are expected after the discovery of graphene and graphynes. Liu et al. [116] proposed a buckled carbon sheet with tetrarings ‘T-graphene’ and thought it as a 2D Dirac material. Although T-graphene was demonstrated to be metallic soon later [117], it drives people to notice the 2D carbon allotropes with multi-member rings. In 2014, Xu et al. [23] designed a series of 2D carbon allotropes by reconstructing graphene. Moreover, three stable rectangular systems named as S-graphene (Fig. 6a), D-graphene (Fig. 6b), and E-graphene (Fig. 6c) were predicted to have Dirac cones [23]. S-graphene consists of eight C atoms in a unit cell and contains four- and six-member rings. It has two self-doping Dirac cones along M-X and Y-Γ lines separately. D-graphene is composed of sp and sp2 C atoms while E-graphene composed of sp3 and sp2 C atoms. Both of them have distorted Dirac cones, but the Dirac points are on different positions. More carbon sheets with different rings and hybridized forms have also been suggested to be Dirac materials [23]. Figure 6. Open in new tabDownload slide Atomic and electronic structures of (a–c) S-graphene, D-graphene, E-graphene [23] (Copyright 2014 Royal Society of Chemistry), and (d) Pmmn boron [24] (Copyright 2014 American Physical Society). Figure 6. Open in new tabDownload slide Atomic and electronic structures of (a–c) S-graphene, D-graphene, E-graphene [23] (Copyright 2014 Royal Society of Chemistry), and (d) Pmmn boron [24] (Copyright 2014 American Physical Society). Boron—the neighbor of carbon in the periodic table—is a fascinating element because of its chemical and structural complexity. The energy landscape of boron clusters is glass-like, being in sharp contrast to carbon and boron nitride systems [118]. Boron might be the second element that can possess freestanding flat monolayer structures, evidenced by the experimental synthesis of single-walled and multiwalled boron nanotubes and the stability calculation of boron sheets [119]. Boron sheets are usually composed of triangular and hexagonal motifs arising from the competition between two- and three-center bonding [120]. Recently, a novel 2D boron sheet (Pmmn boron, Fig. 6d) has been predicted to have distorted Dirac cones [24]. It has eight atoms in a rectangular cell with the lattice constants of 4.52 and 3.26 Å. The B atoms form a buckled triangular layer that can be classified into two sublattices (buckled chains and hexagons, illustrated in Fig. 6d by different colors). The predicted Dirac points locate at the Γ-X line, and the hybrid of in-plane (px orbitals from the buckled boron chains) and out-of-plane states (pz orbitals from the buckled irregular boron hexagons) are responsible for the emergence of Dirac cone [24]. Transition-metal oxides: (VO2)n/(TiO2)m and (CrO2)n/(TiO2)m Dirac cones also exist in multilayer nanostructures of transition-metal oxides. The system (VO2)n/(TiO2)m, formed by n VO2 monolayers and m TiO2 monolayers grown along the rutile (001) direction, was predicted to possess Dirac cones for n = 3 and 4, where the spectrum in the vicinity of the Fermi level is dominated by V 3d orbitals [30]. The in-plane cell of this pseudo-2D system is square, but it does not hold the full square symmetry and thus escape the restriction of the full square symmetry on the existence of Dirac cones [121]. Remarkably, although the dispersion is linear along the diagonal (Γ-M) direction, it is quadratic perpendicular to the diagonal [30,122]: \begin{equation} E\!\left( {\bf k} \right) = \pm \sqrt {\left( {\upsilon _{\rm F} k_1 } \right)^2 + \left( {{{k_2^2 } \Big/ {2m}}} \right)^2 } , \end{equation} (13) where υF and m are the Fermi velocity and effective mass, respectively. It is different from the usual Dirac cones as that in graphene, and is thus referred to as semi-Dirac cones [30]. The strongly anisotropic dispersion is expected to give rise to peculiar transport and thermodynamic properties. Similar semi-Dirac point has also been found in a photonic crystal consisting of a square array of elliptical dielectric cylinders [123] and a microwave graphene analog composed of coupled dielectric cylindrical resonators [124]. Recently, a close inspection has revealed that the semi-Dirac cones in (VO2)n/(TiO2)m are formed by merging three usual Dirac cones at each semi-Dirac point, which results in a non-trivial Chern number in this system (H. Q. Huang, et al. Chern insulator in TiO2/VO2 nanostructures, unpublished). In addition, usual Dirac cones have been predicted in (CrO2)n/(TiO2)m for n = 4 [3,125]. Organic and organometallic crystals Another class of materials which affords Dirac cones are organic conductors. A realized example is the quasi-2D organic conductor α-(BEDT-TTF)2I3 under high pressure [28,126]. α-(BEDT-TTF)2I3 is composed of conductive layers of BEDT-TTF molecules and insulating layers of |${\rm I}_3^ -$| anions, and the in-plane cell is rectangular [126]. According to the band calculation, α-(BEDT-TTF)2I3 possesses two Dirac points at the Fermi level, which, in contrast to the case of graphene, locate at non-symmetric positions in the k-space [28]. In addition, the Dirac cones are heavily tilted and the Fermi velocity along k is not equal to −k, i.e. the inversion symmetry of the Dirac cones is lost. Experimentally, the zero-mode Landau level expected to appear at the Dirac points has been successfully detected in the magnetoresistance measurement [127]. The molecular interactions in organic crystals such as α-(BEDT-TTF)2I3 are much weaker than the interactions in atomic crystals such as graphene and silicene. As a result, the energy scale in organic crystals is small and their applications at the room temperature are limited. To remedy this shortcoming, metal atoms can be introduced to glue organic compounds. Organometallic crystals of Pb2(C6H4)3, Ni2(C6H4)3, and Co2(C6H4)3 with a hexagonal lattice have recently been proposed to possess Dirac cones [22,26]. Due to the incorporation of the magnetic Ni and Co atoms, Ni2(C6H4)3 and Co2(C6H4)3 are half-metallic, which may extend the applications of 2D Dirac materials in spintronics [22]. Systems with a pseudospin of S= 1 Dirac cones in most Dirac materials can be described by Dirac–Weyl equation with pseudospin S = 1/2 as Equation (5). However, Dirac cones with unusual spin are also possible. Wang et al. [21] have constructed various carbon allotropes with square symmetry by enumerating the carbon atoms in a unit cell up to 12. One of them, SG-10b (Fig. 7a), was predicted to possess Dirac cones with a pseudospin of S = 1 in a TB calculation [21]. Two cone-like bands and one flat band contact at the M point at the Fermi level (Fig. 7a), which can be described by a Dirac–Weyl Hamitonian: \begin{equation} H = \upsilon _F {\bf q} \cdot S, \end{equation} (14) where S is a pseudospin of the quantum number 1 with three eigenvalues −1, 0 and +1. The Massless Dirac–Weyl fermions with a pseudospin S = 1 have interesting properties that are different from S = 1/2. The Berry's phase enclosing the Dirac point of S = 1 vanishes, so the Dirac point of S = 1 can exist solely, unlike those of S = 1/2 that appears in pairs. The enlarged pseudospin also leads to an enhanced Klein tunneling where the barrier is transparent for all incident angles [128]. In addition, the presence of the flat band under S = 1 serves as a good starting point to study correlated electron systems [129]. Figure 7. Open in new tabDownload slide Atomic and band structures for (a) SG-10b [21] (Copyright 2013 AIP Publishing LLC), (b) a single-walled hexagonal graphene antidote lattices (GALs) [21] (Copyright 2013 AIP Publishing LLC), and (c) so-MoS2 [25] (Copyright 2014 American Physical Society). Figure 7. Open in new tabDownload slide Atomic and band structures for (a) SG-10b [21] (Copyright 2013 AIP Publishing LLC), (b) a single-walled hexagonal graphene antidote lattices (GALs) [21] (Copyright 2013 AIP Publishing LLC), and (c) so-MoS2 [25] (Copyright 2014 American Physical Society). SG-10b is not the only atomic crystal predicted to possess Dirac cones with S = 1. The odd-walled hexagonal graphene antidote lattices also possess massless Dirac fermions with S = 1 at the Fermi level determined by first-principles calculations (Fig. 7b) [20,21]. Another example is a MoS2 allotrope [25]. Although normal monolayer MoS2 with a hexagonal lattice (h-MoS2) is a direct-gap semiconductor [130], the square MoS2 sheet (so-MoS2) has been shown to possess massless Dirac fermions from d electrons (Fig. 7c) [25]. The single Dirac point locates at the Γ point, so it can be recognized to be S = 1. Recently, a modified TB model based on T3-lattice has provided a continuous interpolation between the fermions with S = 1/2 and the fermions with S = 1 [131]. More 2D Dirac systems with S = 1 can be expected. Artificial lattices: electron gases and ultracold atoms The extraordinary properties of graphene due to the existence of Dirac cones have stimulated the search of the ‘quantum simulators’ of Dirac cones in artificial systems with distinct mechanisms, e.g. 2D electron gases, ultracold atoms, and photonic crystals [132]. Ultracold atoms trapped in optical lattices serve to mimic condensed matter phenomena very well. Realizing Dirac cones with cold atoms loaded in a honeycomb optical lattice was first proposed in 2007 [29]. Cold atom systems with orbital degrees of freedom can even afford new features which are impossible to realize in graphene, e.g. px,y-orbital counterpart of graphene with flat bands and Wigner crystallization [133,134]. Cold atoms with other structures of symmetries can also afford Dirac cones. For example, atoms trapped in a T3 lattice or a line-centered-square lattice behave as the massless Dirac fermions with pseudospin S = 1 [135,136]. In 2012, the experimental creation of Dirac cones at K and K′ was reported in an ultracold gas of 40K atoms trapped in a 2D honeycomb optical potential lattice of interfering laser beams [137]. Moreover, with changing the lattice parameters, the positions of the two Dirac points inside the Brillouin zone move accordingly, and finally merge and annihilate each other when the parameters exceed a critical limit—a situation that is extremely challenging to observe in solids. 2D electron gases are gases of electrons free to move in two dimensions but tightly confined in the third. They are usually found in nanostructures made from semiconductors. A theoretical analysis on the electronic structure of 2D electron gases under an external periodic potential of hexagonal symmetry was conducted in 2009, which revealed that Dirac cones are generated at K and K′ as those in graphene [138,139]. The required parameters were estimated to be achievable under laboratory conditions for semiconductors. In 2012, by combining both experimental and theoretical studies, it was revealed that the main obstacle preventing the realization of Dirac-like physics in 2D electron gases for semiconductors is related to the interplay between the electron density and the lattice constant [140]. In the same year, a different scheme was developed to successfully achieve Dirac cones at a Cu surface [27]. Cu(111) surface is characterized by its nearly free electron gases with very long coherence lengths. Adsorbed CO molecules on clean Cu(111) were assembled into a lattice by positioning the molecules individually using the tip of a STM, which exert periodic potential over the free electron gases. With CO molecules in a hexagonal lattice, the existence of linearly dispersing, massless quasiparticles, and a DOS characteristic of graphene have been clearly shown in the experiments [27]. More importantly, the CO lattice sites can be arbitrarily tuned globally or locally to mimic various effects, providing a versatile means of synthesizing exotic electronic phases in condensed matter. Photonic crystals are periodic optical nanostructures that affect the motion of photons in a way similar to how ionic lattices affect electrons in solids. Photonic crystals offer a distinctive route to obtain dispersion relations with characteristic Dirac cones [132,141,142]. However, photons are boson, being markedly different from the other systems discussed in this review which are fermion, and photons rarely exist in the equilibrium state. Therefore, advances of photonic crystals with Dirac cones will not be addressed here. Readers are directed to reference [132] which provides an excellent review in this field. New 2D Dirac systems keep springing up [143], but they are rather rare compared to the numerous 2D materials. WHY ARE DIRAC MATERIALS SO RARE? Moving and merging of Dirac points Dirac cones are rather robust under perturbation. For example, when a uniaxial or shear strain is applied, the band structure of graphene keeps gapless and the Dirac point moves to a new k location near the original one [51,52]. A TB analysis on the honeycomb lattice with the nearest neighboring interactions indicats that Dirac cones exist in the system as long as three hopping integrals (t1, t2, t3) can form a triangle [144]: \begin{equation} \left| {t_2 - t_1 } \right| \le \left| {t_3 } \right| \le \left| {t_2 + t_1 } \right|. \end{equation} (15) For a general 2D system that contains two atoms of the same species in each unit cell, the criterion for the existence of Dirac cones becomes \begin{equation} \left| {t_1^2 + t_2^2 - t_3^2 - t_4^2 } \right| \le \left| {2t_1 t_2 - 2t_3 t_4 } \right|, \end{equation} (16) and the locations of Dirac cones are given explicitly in a function of four hopping integrals (t1, t2, t3, t4) [121]. The effect of strain on Dirac cones can be well explained in terms of the variation of hopping integrals [52]. When there is no external magnetic field, the system is invariant under the time-reversal operation, where a Dirac point at k0 is always accompanied by its pair at −k0, e.g. K and K′ for graphene. The Berry phase around a Dirac point is either +π or −π, and the Berry phase around a Dirac point at k0 is always opposite to that around its pair at −k0. When two Dirac points with opposite Berry phases move in the k-space under any perturbation and arrive at the same point, they merge and their Berry phases annihilate each other [121,145,146]. At the merging point, the energy dispersion is semi-Dirac which remains linear in one direction but becomes quadratic in the other. Upon further perturbation, a gap is induced and the Dirac points disappear. If there is only a single pair of Dirac points in the system, they can be merged only at high symmetric k-points. For example, the Dirac point observed near X in 6,6,12-graphyne is actually formed by merging two Dirac points [18,19]. An in-plane AC electric field can also drive the moving and merging of Dirac cones [147]. On the other hand, the opposite processes are also possible, i.e. the energy gap may disappear and a pair of Dirac points emerges at high symmetric k-points. For example, DFT calculation showed that a uniaxial strain of 9% induced a pair of Dirac points at Γ in graphdiyne [106]. To achieve Dirac cone, merging or emerging is not an easy task in experiments. The critical strain required to annihilate Dirac cones in graphene is as high as 24% due to the well separated Dirac points originally at K and K′ [51]. Such a large strain is experimentally impractical. A possible way to reduce the modulation difficulty is to replace some single bonds in graphene with acetylenic linkages (–C≡C–), i.e. to adopt graphynes or graphdiynes. For example, one pair of Dirac points in 6,6,12-graphyne locates very close to the high symmetric point X, which is possible to annihilate at small strain. Unfortunately, the controlled preparation of graphynes and graphdiynes is difficult by itself. At present, Dirac cone merging is achieved only in artificial honeycomb lattices where parameters are much more adjustable. By patterning CO molecules on clean Cu(111), the hexagonal potential lattice of electron gases was effectively modulated to demonstrate a transition from massless to massive Dirac fermions in the system [27]. In an ultracold gas of 40K atoms trapped in a 2D honeycomb optical potential lattice, the merging and annihilation of two Dirac points were clearly recorded when the lattice anisotropy exceeded a critical limit [137]. Interestingly, more than two Dirac points can also merge together. When three Dirac points merge, their Berry phases cannot completely annihilate each other and thus possess a non-trivial topology. Quasi-2D (VO2)3/(TiO2)m possesses semi-Dirac cones as described by Equation (13) in its band structure, but the Berry phase around the semi-Dirac points is ±π and the Chern numbers are non-zero, which were actually caused by the merging of three Dirac points (H. Q. Huang, et al. Chern insulator in TiO2/VO2 nanostructures, unpublished). When four Dirac points merge together, the energy dispersion becomes parabolic in any direction near the merging points. This is what happens in bilayer graphene. Its nature of four-Dirac-point-merging can be clearly demonstrated by exerting the next-nearest-neighbor interlayer hopping: each merging point will split into four conventional Dirac points with linear dispersions [12]. The Berry phase around the merging points in bilayer graphene is +2π (−2π), suggesting that they are merged by three +π (−π) and one −π (+π). An illustration for the various situations is given in Fig. 8. Although the merging of two and three Dirac cones both lead to a semi-Dirac cones near which the energy dispersion is linear in one principal axis and quadratic in the other, the band structure in the merging of three Dirac cones is more twisted, appearing in a banana shape (H. Q. Huang, et al. Chern insulator in TiO2/VO2 nanostructures, unpublished). Figure 8. Open in new tabDownload slide The merging of (a) one, (b) two, (c) three, and (d) four Dirac points (H.Q. Huang, et al. Chern insulator in TiO2/VO2 nanostructures, unpublished). 3D band structure and 2D energy contor are shown at top and bottom panels, respectively. A simplified 2 × 2 Hamiltonian of |$\left[{ { 0 \quad H_{12}} \atop {{\rm c.c.} \quad 0}} \right]$| is adopted with the expression of H12 given at the bottom for each case. Figure 8. Open in new tabDownload slide The merging of (a) one, (b) two, (c) three, and (d) four Dirac points (H.Q. Huang, et al. Chern insulator in TiO2/VO2 nanostructures, unpublished). 3D band structure and 2D energy contor are shown at top and bottom panels, respectively. A simplified 2 × 2 Hamiltonian of |$\left[{ { 0 \quad H_{12}} \atop {{\rm c.c.} \quad 0}} \right]$| is adopted with the expression of H12 given at the bottom for each case. Existent conditions of Dirac cones Dirac points are not just energy band crossings but also singularities in the spectrum of Hamiltonians [148]. Dirac cones are usually topologically protected and have inherent connection to the Z2 quantization defined in terms of Berry phases [149,150]. However, why Dirac cones are not ubiquitous but rather rare in 2D materials? Briefly speaking, the rarity of Dirac cones in 2D systems can be explained with the von Neumann–Wigner theorem [151,152]. Considering a general two-band system whose Hamiltonian is given as \begin{equation} H\!\left( {\bf k} \right) = \left[ {\begin{array}{*{20}c} {H_{11} \left( {\bf k} \right)} &\quad {H_{12} \left( {\bf k} \right)} \\ {{\rm c}{\rm .c}{\rm .}} &\quad {H_{22} \left( {\bf k} \right)} \\ \end{array}} \right], \end{equation} (17) and we obtain \begin{eqnarray} &&\!\!\!\!\!\! E_ \pm \left( {\bf k} \right) = \frac{{H_{11} \left( {\bf k} \right) + H_{22} \left( {\bf k} \right)}}{2} \nonumber\\ &\!\!\!\! \pm &\sqrt {\left( {\frac{{H_{11} \left( {\bf k} \right) - H_{22} \left( {\bf k} \right)}}{2}} \right)^2 + \left| {H_{12} \left( {\bf k} \right)} \right|^2 } .\nonumber\\ \end{eqnarray} (18) The conditions to determine the Dirac point are thus \begin{equation} \left\{ \begin{array}{l} H_{11} \left({\bf k} \right) - H_{22} \left({\bf k} \right) = 0 \\ {\mathop{\rm Re}\nolimits} \left[ {H_{12} \left({\bf k} \right)} \right] = 0 \\ {\mathop{\rm Im}\nolimits} \left[ {H_{12} \left({\bf k} \right)} \right] = 0 \end{array} \right., \end{equation} (19) where |${\mathop{\rm Re}\nolimits} \left[ {H_{12} \left({\bf k} \right)} \right]$| and |${\mathop{\rm Im}\nolimits} \left[ {H_{12} \left({\bf k} \right)} \right]$| denote the real and imaginary part of the complex number H12. The three conditions in Equation (19) must be simultaneously fulfilled to have an energy degeneracy, which is known as the von Neumann–Wigner theorem [151]. In 2D systems, we have two variables (kx, ky) to be solved from the above three equations. Because the number of variables is less than the number of equations, the problem is overdetermined and we usually get no solution. This is the main reason why Dirac cones are so rare. To make the existence of Dirac points possible, a constraint on H(k) is required, which would reduce the number of conditions in Equation (19) by one [152]. The required constraint is usually attributed to the symmetries of the system. A representative example is the space–time inversion symmetry [152,153]. The invariance under the time reversal, T, gives \begin{equation} T:H\left( {\bf k} \right) = H^* \left( { - {\bf k}} \right), \end{equation} (20) and the spatial inversion gives \begin{equation} I:H\!\left( {\bf k} \right) = {\boldsymbol \sigma }_x H\left( { - {\bf k}} \right){\boldsymbol \sigma }_x , \end{equation} (21) where σx is a Pauli matrix. Under a combination of T and I, \begin{equation} TI:H\!\left( {\bf k} \right) = {\boldsymbol \sigma }_x H^* \left( {\bf k} \right){\boldsymbol \sigma }_x , \end{equation} (22) which leads to H11(k) = H22(k) and cancels the first equation in Equation (19), making the existence of Dirac cones possible. Most Dirac materials, if not all, are actually protected from this kind of space–time inversion symmetry. Not all symmetries are effective in providing constraint on H(k). Generally, a symmetry operation would relate the Hamiltonian at a point k to that at another point k′, e.g. Equations (20) and (21). Although they are useful properties, what we need here is a constraint on the form of H(k). Therefore, to provide a constraint, the k-points should be kept unchanged after the symmetry operation (invariant k-point). The combination of T and I meets such a requirement for a general k-point as shown in Equation (22), but the separate T or I does not. Similarly, a sole rotation cannot provide an effective constraint either and is not essential for the existence of Dirac cones. For example, when a shear strain is applied upon graphene, although the rotation symmetries are destroyed, Dirac cones exist in the system before the strain exceeds a large threshold [51,52]. Although the achieved constraint from the system symmetries makes the existence of Dirac cones possible, it does not guarantee it. The major reason is that kx and ky are real numbers and appear in H(k) in the form of a sine or cosine function. The solution of Equation (19) may be still absent for certain parameter values even with the constraint [121]. Equations (15) and (16) are two examples of the required parameters to guarantee Dirac cones. Another excellent example comes from a comparison between β- and γ-graphynes. β- and γ-graphynes share completely identical symmetries, but β-graphyne is semimetal with Dirac cone while γ-graphyne is semiconductor. A detailed analysis reveals that the hopping integrals in β-and γ-graphynes are responsible for their existence/absence of Dirac cones [19]. Hopping integrals are affected by atomic geometry such as bond length, which is further related to the crystal lattice. A simplified analysis on a general 2D system, which contains two atoms of the same species in each unit cell, demonstrates that a hexagonal cell is the most favorable for the existence of Dirac cones, and the favorableness gradually diminishes when the cell evolves into a square one (Fig. 9) [121]. This is consistent with the fact that many Dirac materials are observed in hexagonal lattice. Figure 9. Open in new tabDownload slide Hexagonal cell is favorable for the existence of Dirac cones. A system containing two atoms in each unit cell is considered with the first atom located at corners. Location regions of the second atom, which enable the existence of Dirac cones, are marked by the filled areas [121]. Figure 9. Open in new tabDownload slide Hexagonal cell is favorable for the existence of Dirac cones. A system containing two atoms in each unit cell is considered with the first atom located at corners. Location regions of the second atom, which enable the existence of Dirac cones, are marked by the filled areas [121]. To observe the novel properties particular to the Dirac cones, the Fermi level should lie at the Dirac points and no other bands go through the Fermi level. There should not have any other bands than Dirac points overlap at the Fermi level, otherwise the low-energy excitation would be dominated by the carriers in the hole and electron pockets caused by the overlap. For example, although quasi-2D Dirac cones were predicted in LaAgBi2 [154], they may be difficult to detect experimentally. The band overlap can be avoided when the diagonal elements in Equation (17), H11(k) and H22(k), are k-independent. This condition is automatically satisfied in a bipartite system under a nearest-neighbor approach. Actually, all Dirac materials listed in Fig. 1 are bipartite except Pmmn boron. It is noted that being bipartite is not sufficient condition to possess Dirac cones. For example, the single-layer h-BN sheet has a bipartite lattice as graphene, but its two k-independent H11 and H22 are not equal to each other due to the different onsite energies of B and N atoms, so it is semiconductor with a large energy gap [155]. Symmetries are usually needed to produce Dirac cones, but too high symmetries may be disadvantageous. If there is too much constraint on H(k), the number of independent equations may be smaller than the number of independent variables, and thus the obtained solutions are not discrete points, but a continuous line in the k-space. In such a case, the carriers near the Fermi level cannot be described by massless Dirac equation, and the corresponding novel properties are absent. For 2D carbon allotropes with full square symmetry (wallpaper group p4mm), many systems have band contact at the Fermi level appearing in lines [21]. Therefore, they are metal, and the absence of Dirac cones in them is not caused by the absence of contacts between valence and conductance bands, but caused by the existence of too much contacts. The only carbon allotrope with full square symmetry recognized to possess Dirac cone, SG-10b, has three bands contacting at the Dirac point, but not two, and cannot be described by the Equation (17) [21]. A study on a general 2D atomic crystal containing two atoms in each unit cell and a 2D electron gas under a periodic muffin-tin potential also suggested that conventional Dirac cones with S = 1/2 (such as those observed in graphene) are difficult, if not impossible, to be achieved under full square symmetry [121]. As a support for this opinion, it is noted that although the Dirac material (VO2)n/(TiO2)m and (CrO2)n/(TiO2)m have a pseudo-2D square lattice, they do not possess full square symmetry. Based on the analyses given above, we can have a brief summary on why Dirac cones are so rare. The central origin is the von Neumann–Wigner theorem: the number of variables (kx and ky) is usually less than the number of equations to determine the Dirac points. In more details, to achieve Dirac materials, at least three conditions are required: (1) Symmetry. Specific symmetries are required to reduce the number of equations to be solved. The k-points should be unchanged after the symmetry operation (invariant k-point). Too low or too high symmetries are both disadvantageous. (2) Parameters. Even when the number of equations is equal to the number of variables, the solution does not necessarily exist since the variables (kx and ky) are real numbers and appear in the equations in the form of a sine or cosine function. Therefore, proper parameters are required. This is usually described as a phase diagram in the parameter space. (3) Fermi level and band overlap. The Fermi level should lie at the Dirac points and there should not be any other band than Dirac points overlap at the Fermi level. The rigorous conditions may explain the scarcity of Dirac materials. It is noted that the energy dispersion near the Dirac point solved from the above conditions is linear in most cases, because a quadratic dispersion needs more constraints to satisfy dE/dk = 0 at the crossing point. The understanding on the conditions for Dirac materials also provides a guide to search/design new Dirac materials in the future. Structural features of the known 2D Dirac materials After discussing the conditions for Dirac cones to exist, we now briefly analyze the common structural features of the known atomic Dirac materials to see how they are beneficial for the existence of Dirac cones. First, most Dirac materials have spatial inversion symmetry. As demonstrated above, in combination with time reversal, spatial inversion provides an effective constraint on H(k) to reduce the number of equation. The exceptions are transition-metal oxides of (VO2)n/(TiO2)m and (CrO2)n/(TiO2)m. Their Dirac points are protected by mirror symmetry. As a result, the Dirac points locate along the Γ-M lines which are invariant under mirror operation. Secondly, all atomic Dirac materials have even number of atoms in a unit cell. Many of them are bipartite and composed of only one element (mainly belongs to group IVA). The bipartite feature forces the diagonal Hamiltonian elements to be constants, and thus is favorable not only in reducing the equation number but also in avoiding band overlap at the Fermi level. Single element, on the other hand, is helpful for the on-site energies at different sites to be equal. The different properties of graphene and h-BN can be explained in terms of relation between the element number and the on-site energy. We can transform the materials into binary systems by introducing the concept of superatoms. In detail, we replace half atoms in a unit cell of each system in Fig. 1 except (VO2)n/(TiO2)m by one superatom and the other half by another superatom (Fig. 10c), and redraw the structures. Then, graphene, α-graphyne, silicene, germanene, CO on Cu(111) and Pb2(C6H4)3 can be transformed to a graphene-like structure (Fig. 10a), and the 6,6,12-graphyne, 14,14,18-graphyne, S-graphenen, SG-10b, Pmmn boron, and so-MoS2 can be transformed to a rectangular or square structure constituted by two superatoms in spatial inversion (Fig. 10b). Figure 10. Open in new tabDownload slide (a) Honeycomb and (b) rectangular structure for 2D Dirac materials by introducing (c) the concept of superatoms. Topological structures for (d) graphene, silicene, germanene, graphynes, and (e) Pmmn boron. Figure 10. Open in new tabDownload slide (a) Honeycomb and (b) rectangular structure for 2D Dirac materials by introducing (c) the concept of superatoms. Topological structures for (d) graphene, silicene, germanene, graphynes, and (e) Pmmn boron. Lastly, hexagonal honeycomb structure is common in atomic Dirac materials. This feature may relate to the requirement of parameters. As illustrated in Fig. 9, hexagonal cell is more favorable in providing appropriate parameters for the existence of Dirac cones. We take an effective bond method to analyze the topological structure of the systems. Assume that C–C, Si–Si, Ge–Ge, and C–C≡C–C bonds can be expressed as one effective bond, and then the elementary substances in Fig. 1 have the same honeycomb topology (Fig. 10d). For Pmmn boron, since it has two sublattices which contribute the valence and conduction band separately [24], we connected the atoms in the same sublattice and found that two honeycomb structures were formed (Fig. 10e). CONCLUSIONS Dirac cones are not only the linear energy dispersion around discrete points but also singularities in the spectrum of Hamiltonians and are topologically protected. Dirac-cone structures give materials unusual properties and promising prospects in both fundamental and applied research. However, the probability of Dirac cones in 2D systems is rather low. Herein, we have systematically surveyed the known 2D Dirac materials and discussed how Dirac cones emerge and merge in the systems. Rigorous conditions on symmetry and parameters are required to achieve Dirac cones in 2D systems, which provide an explanation for the rarity of Dirac materials. Looking forward, we believe that more and more 2D Dirac materials will be discovered, and a thorough understanding on the existing conditions of Dirac cones is greatly helpful in seeking/designing new systems. We are grateful to Wenhui Duan, Huaqing Huang, Shuqing Zhang, and Zhenzhu Li for discussions. FUNDING This work was supported by the National Natural Science Foundation of China (21373015). REFERENCES 1. Novoselov KS , Geim AK , Morozov SV , et al. Electric field effect in atomically thin carbon films , Science , 2004 , vol. 306 (pg. 666 - 9 ) Google Scholar Crossref Search ADS PubMed WorldCat 2. Novoselov KS , Fal’ko VI , Colombo L , et al. A roadmap for graphene , Nature , 2012 , vol. 490 (pg. 192 - 200 ) Google Scholar Crossref Search ADS PubMed WorldCat 3. Abergel DSL , Apalkov V , Berashevich J , et al. Properties of graphene: a theoretical perspective , Adv Phys , 2010 , vol. 59 (pg. 261 - 482 ) Google Scholar Crossref Search ADS WorldCat 4. Novoselov KS , Geim AK , Morozov SV , et al. Two-dimensional gas of massless Dirac fermions in graphene , Nature , 2005 , vol. 438 (pg. 197 - 200 ) Google Scholar Crossref Search ADS PubMed WorldCat 5. Zhang YB , Tan YW , Stormer HL , et al. Experimental observation of the quantum Hall effect and Berry's phase in graphene , Nature , 2005 , vol. 438 (pg. 201 - 4 ) Google Scholar Crossref Search ADS PubMed WorldCat 6. Bolotin KI , Ghahari F , Shulman MD , et al. Observation of the fractional quantum Hall effect in graphene , Nature , 2009 , vol. 462 (pg. 196 - 9 ) Google Scholar Crossref Search ADS PubMed WorldCat 7. Du X , Skachko I , Duerr F , et al. Fractional quantum Hall effect and insulating phase of Dirac electrons in graphene , Nature , 2009 , vol. 462 (pg. 192 - 5 ) Google Scholar Crossref Search ADS PubMed WorldCat 8. Dean CR , Wang L , Maher P , et al. Hofstadter's butterfly and the fractal quantum Hall effect in moire superlattices , Nature , 2013 , vol. 497 (pg. 598 - 602 ) Google Scholar Crossref Search ADS PubMed WorldCat 9. Ponomarenko LA , Gorbachev RV , Yu GL , et al. Cloning of Dirac fermions in graphene superlattices , Nature , 2013 , vol. 497 (pg. 594 - 7 ) Google Scholar Crossref Search ADS PubMed WorldCat 10. Hunt B , Sanchez-Yamagishi JD , Young AF , et al. Massive Dirac fermions and Hofstadter butterfly in a van der Waals heterostructure , Science , 2013 , vol. 340 (pg. 1427 - 30 ) Google Scholar Crossref Search ADS PubMed WorldCat 11. Bolotin KI , Sikes KJ , Jiang Z , et al. Ultrahigh electron mobility in suspended graphene , Solid State Commun , 2008 , vol. 146 (pg. 351 - 5 ) Google Scholar Crossref Search ADS WorldCat 12. Castro Neto AH , Guinea F , Peres NMR , et al. The electronic properties of graphene , Rev Mod Phys , 2009 , vol. 81 (pg. 109 - 62 ) Google Scholar Crossref Search ADS WorldCat 13. Weiss NO , Zhou H , Liao L , et al. Graphene: an emerging electronic material , Adv Mater , 2012 , vol. 24 (pg. 5782 - 825 ) Google Scholar Crossref Search ADS PubMed WorldCat 14. Butler SZ , Hollen SM , Cao LY , et al. Progress, challenges, and opportunities in two-dimensional materials beyond graphene , ACS Nano , 2013 , vol. 7 (pg. 2898 - 926 ) Google Scholar Crossref Search ADS PubMed WorldCat 15. Xu MS , Liang T , Shi MM , et al. Graphene-like two-dimensional materials , Chem Rev , 2013 , vol. 113 (pg. 3766 - 98 ) Google Scholar Crossref Search ADS PubMed WorldCat 16. Wallace PR . The band theory of graphite , Phys Rev , 1947 , vol. 71 (pg. 622 - 34 ) Google Scholar Crossref Search ADS WorldCat 17. Cahangirov S , Topsakal M , Akturk E , et al. Two- and one-dimensional honeycomb structures of silicon and germanium , Phys Rev Lett , 2009 , vol. 102 pg. 236804 Google Scholar Crossref Search ADS PubMed WorldCat 18. Malko D , Neiss C , Vines F , et al. Competition for graphene: graphynes with direction-dependent Dirac cones , Phys Rev Lett , 2012 , vol. 108 pg. 086804 Google Scholar Crossref Search ADS PubMed WorldCat 19. Huang HQ , Duan WH , Liu ZR . The existence/absence of Dirac cones in graphynes , New J Phys , 2013 , vol. 15 pg. 023004 Google Scholar Crossref Search ADS WorldCat 20. Ouyang FP , Peng SL , Liu ZF , et al. Bandgap opening in graphene antidot lattices: the missing half , ACS Nano , 2011 , vol. 5 (pg. 4023 - 30 ) Google Scholar Crossref Search ADS PubMed WorldCat 21. Wang JY , Huang HQ , Duan WH , et al. Identifying Dirac cones in carbon allotropes with square symmetry , J Chem Phys , 2013 , vol. 139 pg. 184701 Google Scholar Crossref Search ADS PubMed WorldCat 22. Ma YD , Dai Y , Li XR , et al. Prediction of two-dimensional materials with half-metallic Dirac cones: Ni2C18H12 and Co2C18H12 , Carbon , 2014 , vol. 73 (pg. 382 - 8 ) Google Scholar Crossref Search ADS WorldCat 23. Xu LC , Wang RZ , Miao MS , et al. Two dimensional Dirac carbon allotropes from graphene , Nanoscale , 2014 , vol. 6 (pg. 1113 - 8 ) Google Scholar Crossref Search ADS PubMed WorldCat 24. Zhou XF , Dong X , Oganov AR , et al. Semimetallic two-dimensional boron allotrope with massless Dirac fermions , Phys Rev Lett , 2014 , vol. 112 pg. 085502 Google Scholar Crossref Search ADS WorldCat 25. Li WF , Guo M , Zhang G , et al. Gapless MoS2 allotrope possessing both massless Dirac and heavy fermions , Phys Rev B , 2014 , vol. 89 pg. 205402 Google Scholar Crossref Search ADS WorldCat 26. Wang ZF , Liu Z , Liu F . Organic topological insulators in organometallic lattices , Nat Commun , 2013 , vol. 4 pg. 1471 Google Scholar Crossref Search ADS PubMed WorldCat 27. Gomes KK , Mar W , Ko W , et al. Designer Dirac fermions and topological phases in molecular graphene , Nature , 2012 , vol. 483 (pg. 306 - 10 ) Google Scholar Crossref Search ADS PubMed WorldCat 28. Katayama S , Kobayashi A , Suzumura Y . Pressure-induced zero-gap semiconducting state in organic conductor alpha-(BEDT-TTF)(2)I-3 salt , J Phys Soc Jpn , 2006 , vol. 75 pg. 054705 Google Scholar Crossref Search ADS WorldCat 29. Zhu SL , Wang BG , Duan LM . Simulation and detection of Dirac fermions with cold atoms in an optical lattice , Phys Rev Lett , 2007 , vol. 98 pg. 260402 Google Scholar Crossref Search ADS PubMed WorldCat 30. Pardo V , Pickett WE . Half-metallic semi-Dirac-point generated by quantum confinement in TiO2/VO2 nanostructures , Phys Rev Lett , 2009 , vol. 102 pg. 166803 Google Scholar Crossref Search ADS PubMed WorldCat 31. Cai TY , Li X , Wang F , et al. Emergent topological and half semimetallic Dirac fermions at oxide interfaces arXiv:1310.2471 OpenURL Placeholder Text WorldCat 32. Wang JY , Zhao RQ , Liu ZF , et al. Widely tunable carrier mobility of boron nitride-embedded graphene , Small , 2013 , vol. 9 (pg. 1373 - 8 ) Google Scholar Crossref Search ADS PubMed WorldCat 33. Divincenzo DP , Mele EJ . Self-consistent effective-mass theory for intralayer screening in graphite-intercalation compounds , Phys Rev B , 1984 , vol. 29 (pg. 1685 - 94 ) Google Scholar Crossref Search ADS WorldCat 34. Zheng YS , Ando T . Hall conductivity of a two-dimensional graphite system , Phys Rev B , 2002 , vol. 65 pg. 245420 Google Scholar Crossref Search ADS WorldCat 35. Dean CR , Young AF , Cadden-Zimansky P , et al. Multicomponent fractional quantum Hall effect in graphene , Nat Phys , 2011 , vol. 7 (pg. 693 - 6 ) Google Scholar Crossref Search ADS WorldCat 36. Feldman BE , Levin AJ , Krauss B , et al. Fractional quantum hall phase transitions and four-flux states in graphene , Phys Rev Lett , 2013 , vol. 111 pg. 076802 Google Scholar Crossref Search ADS PubMed WorldCat 37. Hofstadter DR . Energy-levels and wave-functions of bloch electrons in rational and irrational magnetic-fields , Phys Rev B , 1976 , vol. 14 (pg. 2239 - 49 ) Google Scholar Crossref Search ADS WorldCat 38. Adam S , Hwang EH , Galitski VM , et al. A self-consistent theory for graphene transport , Proc Natl Acad Sci USA , 2007 , vol. 104 (pg. 18392 - 7 ) Google Scholar Crossref Search ADS PubMed WorldCat 39. Dean CR , Young AF , Meric I , et al. Boron nitride substrates for high-quality graphene electronics , Nat Nanotechnol , 2010 , vol. 5 (pg. 722 - 6 ) Google Scholar Crossref Search ADS PubMed WorldCat 40. Chen J-H , Jang C , Xiao SD , et al. Intrinsic and extrinsic performance limits of graphene devices on SiO2 , Nat Nanotechnol , 2008 , vol. 3 (pg. 206 - 9 ) Google Scholar Crossref Search ADS PubMed WorldCat 41. Morozov SV , Novoselov KS , Katsnelson MI , et al. Giant intrinsic carrier mobilities in graphene and its bilayer , Phys Rev Lett , 2008 , vol. 100 pg. 016602 Google Scholar Crossref Search ADS PubMed WorldCat 42. Castro EV , Ochoa H , Katsnelson MI , et al. Limits on charge carrier mobility in suspended graphene due to flexural phonons , Phys Rev Lett , 2010 , vol. 105 pg. 266601 Google Scholar Crossref Search ADS PubMed WorldCat 43. Park CH , Bonini N , Sohier T , et al. Electron-phonon interactions and the intrinsic electrical resistivity of graphene , Nano Lett , 2014 , vol. 14 (pg. 1113 - 9 ) Google Scholar Crossref Search ADS PubMed WorldCat 44. Li ZZ , Wang JY , Liu ZR . Intrinsic carrier mobility of Dirac cones: the limitations of deformation potential theory , J Chem Phys , 2014 , vol. 141 pg. 144107 Google Scholar Crossref Search ADS PubMed WorldCat 45. Meyer JC , Geim AK , Katsnelson MI , et al. The structure of suspended graphene sheets , Nature , 2007 , vol. 446 (pg. 60 - 3 ) Google Scholar Crossref Search ADS PubMed WorldCat 46. Xu P , Neek-Amal M , Barber SD , et al. Unusual ultra-low-frequency fluctuations in freestanding graphene , Nat Commun , 2014 , vol. 5 pg. 3720 Google Scholar Crossref Search ADS PubMed WorldCat 47. Guinea F , Katsnelson MI , Vozmediano MAH . Midgap states and charge inhomogeneities in corrugated graphene , Phys Rev B , 2008 , vol. 77 pg. 075422 Google Scholar Crossref Search ADS WorldCat 48. Lui CH , Liu L , Mak KF , et al. Ultraflat graphene , Nature , 2009 , vol. 462 (pg. 339 - 41 ) Google Scholar Crossref Search ADS PubMed WorldCat 49. Zhang YB , Brar VW , Girit C , et al. Origin of spatial charge inhomogeneity in graphene , Nat Phys , 2009 , vol. 5 (pg. 722 - 6 ) Google Scholar Crossref Search ADS WorldCat 50. Burson KM , Cullen WG , Adam S , et al. Direct imaging of charged impurity density in common graphene substrates , Nano Lett , 2013 , vol. 13 (pg. 3576 - 80 ) Google Scholar Crossref Search ADS PubMed WorldCat 51. Pereira VM , Neto AHC , Peres NMR . Tight-binding approach to uniaxial strain in graphene , Phys Rev B , 2009 , vol. 80 pg. 045401 Google Scholar Crossref Search ADS WorldCat 52. Li Y , Jiang XW , Liu ZF , et al. Strain effects in graphene and graphene nanoribbons: the underlying mechanism , Nano Res , 2010 , vol. 3 (pg. 545 - 56 ) Google Scholar Crossref Search ADS WorldCat 53. Choi SM , Jhi SH , Son YW . Effects of strain on electronic properties of graphene , Phys Rev B , 2010 , vol. 81 pg. 081407(R) Google Scholar Crossref Search ADS WorldCat 54. Pereira VM , Neto AHC . Strain engineering of graphene's electronic structure , Phys Rev Lett , 2009 , vol. 103 pg. 046801 Google Scholar Crossref Search ADS PubMed WorldCat 55. Guinea F , Katsnelson MI , Geim AK . Energy gaps and a zero-field quantum Hall effect in graphene by strain engineering , Nat Phys , 2010 , vol. 6 (pg. 30 - 3 ) Google Scholar Crossref Search ADS WorldCat 56. Levy N , Burke SA , Meaker KL , et al. Strain-induced pseudo-magnetic fields greater than 300 Tesla in graphene nanobubbles , Science , 2010 , vol. 329 (pg. 544 - 7 ) Google Scholar Crossref Search ADS PubMed WorldCat 57. Khomyakov PA , Giovannetti G , Rusu PC , et al. First-principles study of the interaction and charge transfer between graphene and metals , Phys Rev B , 2009 , vol. 79 pg. 195425 Google Scholar Crossref Search ADS WorldCat 58. Liu CC , Feng WX , Yao YG . Quantum spin Hall effect in silicene and two-dimensional germanium , Phys Rev Lett , 2011 , vol. 107 pg. 076802 Google Scholar Crossref Search ADS PubMed WorldCat 59. Ni Z , Liu Q , Tang K , et al. Tunable bandgap in silicene and germanene , Nano Lett , 2012 , vol. 12 (pg. 113 - 8 ) Google Scholar Crossref Search ADS PubMed WorldCat 60. Drummond ND , Zolyomi V , Fal’ko VI . Electrically tunable band gap in silicene , Phys Rev B , 2012 , vol. 85 pg. 075423 Google Scholar Crossref Search ADS WorldCat 61. Ezawa M . A topological insulator and helical zero mode in silicene under an inhomogeneous electric field , New J Phys , 2012 , vol. 14 pg. 033003 Google Scholar Crossref Search ADS WorldCat 62. Ezawa M . Valley-polarized metals and quantum anomalous Hall effect in silicene , Phys Rev Lett , 2012 , vol. 109 pg. 055502 Google Scholar Crossref Search ADS PubMed WorldCat 63. Ezawa M . Photoinduced topological phase transition and a single Dirac-cone state in silicene , Phys Rev Lett , 2013 , vol. 110 pg. 026603 Google Scholar Crossref Search ADS PubMed WorldCat 64. Pan H , Li ZS , Liu CC , et al. Valley-polarized quantum anomalous Hall effect in silicene , Phys Rev Lett , 2014 , vol. 112 pg. 106802 Google Scholar Crossref Search ADS PubMed WorldCat 65. Shao ZG , Ye XS , Yang L , et al. First-principles calculation of intrinsic carrier mobility of silicene , J Appl Phys , 2013 , vol. 114 pg. 093712 Google Scholar Crossref Search ADS WorldCat 66. Ye X-S , Shao Z-G , Zhao H , et al. Intrinsic carrier mobility of germanene is larger than graphene's: first-principle calculations , RSC Adv , 2014 , vol. 4 (pg. 21216 - 20 ) Google Scholar Crossref Search ADS WorldCat 67. Roome NJ , Carey JD . Beyond graphene: stable elemental monolayers of silicene and germanene , ACS Appl Mater Interfaces , 2014 , vol. 6 (pg. 7743 - 50 ) Google Scholar Crossref Search ADS PubMed WorldCat 68. Tsai WF , Huang CY , Chang TR , et al. Gated silicene as a tunable source of nearly 100% spin-polarized electrons , Nat Commun , 2013 , vol. 4 pg. 1500 Google Scholar Crossref Search ADS PubMed WorldCat 69. Niu ZP , Dong SH . Valley and spin thermoelectric transport in ferromagnetic silicene junctions , Appl Phys Lett , 2014 , vol. 104 pg. 202401 Google Scholar Crossref Search ADS WorldCat 70. Lalmi B , Oughaddou H , Enriquez H , et al. Epitaxial growth of a silicene sheet , Appl Phys Lett , 2010 , vol. 97 pg. 223109 Google Scholar Crossref Search ADS WorldCat 71. Feng BJ , Ding ZJ , Meng S , et al. Evidence of silicene in honeycomb structures of silicon on Ag(111) , Nano Lett , 2012 , vol. 12 (pg. 3507 - 11 ) Google Scholar Crossref Search ADS PubMed WorldCat 72. Jamgotchian H , Colignon Y , Hamzaoui N , et al. Growth of silicene layers on Ag(111): unexpected effect of the substrate temperature , J Phys-Condens Mat , 2012 , vol. 24 pg. 172001 Google Scholar Crossref Search ADS WorldCat 73. Liu Z-L , Wang M-X , Xu J-P , et al. Various atomic structures of monolayer silicene fabricated on Ag(111) , New J Phys , 2014 , vol. 16 pg. 075006 Google Scholar Crossref Search ADS WorldCat 74. Enriquez H , Vizzini S , Kara A , et al. Silicene structures on silver surfaces , J Phys-Condens Mat , 2012 , vol. 24 pg. 314211 Google Scholar Crossref Search ADS WorldCat 75. Chen L , Liu CC , Feng BJ , et al. Evidence for Dirac fermions in a honeycomb lattice based on silicon , Phys Rev Lett , 2012 , vol. 109 pg. 056804 Google Scholar Crossref Search ADS PubMed WorldCat 76. De Vogt P , Padova P , Quaresima C , et al. Silicene: compelling experimental evidence for graphenelike two-dimensional silicon , Phys Rev Lett , 2012 , vol. 108 pg. 155501 Google Scholar Crossref Search ADS PubMed WorldCat 77. Arafune R , Lin CL , Nagao R , et al. Comment on ‘Evidence for Dirac fermions in a honeycomb lattice based on silicon’ , Phys Rev Lett , 2013 , vol. 110 pg. 229701 Google Scholar Crossref Search ADS PubMed WorldCat 78. Lin CL , Arafune R , Kawahara K , et al. Substrate-induced symmetry breaking in silicene , Phys Rev Lett , 2013 , vol. 110 pg. 076801 Google Scholar Crossref Search ADS PubMed WorldCat 79. Cahangirov S , Audiffred M , Tang P , et al. Electronic structure of silicene on Ag(111): strong hybridization effects , Phys Rev B , 2013 , vol. 88 pg. 035432 Google Scholar Crossref Search ADS WorldCat 80. Guo ZX , Furuya S , Iwata J , et al. Absence of Dirac electrons in silicene on Ag(111) surfaces , J Phys Soc Jpn , 2013 , vol. 82 pg. 063714 Google Scholar Crossref Search ADS WorldCat 81. Guo ZX , Furuya S , Iwata J , et al. Absence and presence of Dirac electrons in silicene on substrates , Phys Rev B , 2013 , vol. 87 pg. 235435 Google Scholar Crossref Search ADS WorldCat 82. Chen MX , Weinert M . Revealing the substrate origin of the linear dispersion of silicene/Ag(111) , Nano Lett , 2014 , vol. 14 (pg. 5189 - 93 ) Google Scholar Crossref Search ADS PubMed WorldCat 83. Yuan Y , Quhe R , Zheng J , et al. Strong band hybridization between silicene and Ag(111) substrate , Phys E , 2014 , vol. 58 (pg. 38 - 42 ) Google Scholar Crossref Search ADS WorldCat 84. Fleurence A , Friedlein R , Ozaki T , et al. Experimental evidence for epitaxial silicene on diboride thin films , Phys Rev Lett , 2012 , vol. 108 pg. 245501 Google Scholar Crossref Search ADS PubMed WorldCat 85. Meng L , Wang YL , Zhang LZ , et al. Buckled silicene formation on Ir(111) , Nano Lett , 2013 , vol. 13 (pg. 685 - 90 ) Google Scholar Crossref Search ADS PubMed WorldCat 86. Chiappe D , Scalise E , Cinquanta E , et al. Two-dimensional Si nanosheets with local hexagonal structure on a MoS2 surface , Adv Mater , 2014 , vol. 26 (pg. 2096 - 101 ) Google Scholar Crossref Search ADS PubMed WorldCat 87. Quhe R , Yuan YK , Zheng JX , et al. Does the Dirac cone exist in silicene on metal substrates? , Sci Rep , 2014 , vol. 4 pg. 5476 Google Scholar Crossref Search ADS PubMed WorldCat 88. Scalise E , Houssa M , Cinquanta E , et al. Engineering the electronic properties of silicene by tuning the composition of MoX2 and GaX (X = S, Se, Te) chalchogenide templates , 2D Mater , 2014 , vol. 1 pg. 011010 Google Scholar Crossref Search ADS WorldCat 89. Liu H , Gao J , Zhao J . Silicene on substrates: a way to preserve or tune its electronic properties , J Phys Chem C , 2013 , vol. 117 (pg. 10353 - 9 ) Google Scholar Crossref Search ADS WorldCat 90. Li LF , Lu SZ , Pan JB , et al. Buckled germanene formation on Pt(111) , Adv Mater , 2014 , vol. 26 (pg. 4820 - 4 ) Google Scholar Crossref Search ADS PubMed WorldCat 91. Li L , Zhao M . First-principles identifications of superstructures of germanene on Ag(111) surface and h-BN substrate , Phys Chem Chem Phys , 2013 , vol. 15 (pg. 16853 - 63 ) Google Scholar Crossref Search ADS PubMed WorldCat 92. Li XD , Wu SQ , Zhou S , et al. Structural and electronic properties of germanene/MoS2 monolayer and silicene/MoS2 monolayer superlattices , Nanoscale Res Lett , 2014 , vol. 9 pg. 110 Google Scholar Crossref Search ADS PubMed WorldCat 93. Liu CC , Jiang H , Yao YG . Low-energy effective Hamiltonian involving spin-orbit coupling in silicene and two-dimensional germanium and tin , Phys Rev B , 2011 , vol. 84 pg. 195430 Google Scholar Crossref Search ADS WorldCat 94. Garcia JC , de Lima DB , Assali LVC , et al. Group IV graphene- and graphane-like nanosheets , J Phys Chem C , 2011 , vol. 115 (pg. 13242 - 6 ) Google Scholar Crossref Search ADS WorldCat 95. Xu Y , Yan BH , Zhang HJ , et al. Large-gap quantum spin Hall insulators in tin films , Phys Rev Lett , 2013 , vol. 111 pg. 136804 Google Scholar Crossref Search ADS PubMed WorldCat 96. Baughman RH , Eckhardt H , Kertesz M . Structure-property predictions for new planar forms of carbon - layered phases containing sp2 and sp atoms , J Chem Phys , 1987 , vol. 87 (pg. 6687 - 99 ) Google Scholar Crossref Search ADS WorldCat 97. Enyashin AN , Ivanovskii AL . Graphene allotropes , Phys Status Solidi B , 2011 , vol. 248 (pg. 1879 - 83 ) Google Scholar Crossref Search ADS WorldCat 98. Narita N , Nagai S , Suzuki S , et al. Optimized geometries and electronic structures of graphyne and its family , Phys Rev B , 1998 , vol. 58 (pg. 11009 - 14 ) Google Scholar Crossref Search ADS WorldCat 99. Yin WJ , Xie YE , Liu LM , et al. R-graphyne: a new two-dimensional carbon allotrope with versatile Dirac-like point in nanoribbons , J Mater Chem A , 2013 , vol. 1 (pg. 5341 - 6 ) Google Scholar Crossref Search ADS WorldCat 100. Kim BG , Choi HJ . Graphyne: hexagonal network of carbon with versatile Dirac cones , Phys Rev B , 2012 , vol. 86 pg. 115435 Google Scholar Crossref Search ADS WorldCat 101. Liu Z , Yu GD , Yao HB , et al. A simple tight-binding model for typical graphyne structures , New J Phys , 2012 , vol. 14 pg. 113007 Google Scholar Crossref Search ADS WorldCat 102. Yang DZ , Si MS , Zhang GP , et al. Crystal momentum-dependent anisotropy of the Dirac cone in the rectangular carbon allotropes , Europhys Lett , 2014 , vol. 107 pg. 20003 Google Scholar Crossref Search ADS WorldCat 103. Zhao MW , Dong WZ , Wang AZ . Two-dimensional carbon topological insulators superior to graphene , Sci Rep , 2013 , vol. 3 pg. 03532 Google Scholar Crossref Search ADS WorldCat 104. van Miert G , Juricic V , Smith CM . Tight-binding theory of spin-orbit coupling in graphynes , Phys Rev B , 2014 , vol. 90 pg. 195414 Google Scholar Crossref Search ADS WorldCat 105. van Miert G , Smith CM , Juricic V . High-Chern-number bands and tunable Dirac cones in beta-graphyne , Phys Rev B , 2014 , vol. 90 pg. 081406 Google Scholar Crossref Search ADS WorldCat 106. Cui HJ , Sheng XL , Yan QB , et al. Strain-induced Dirac cone-like electronic structures and semiconductor-semimetal transition in graphdiyne , Phys Chem Chem Phys , 2013 , vol. 15 (pg. 8179 - 85 ) Google Scholar Crossref Search ADS PubMed WorldCat 107. Wang GX , Si MS , Kumar A , et al. Strain engineering of Dirac cones in graphyne , Appl Phys Lett , 2014 , vol. 104 pg. 213107 Google Scholar Crossref Search ADS WorldCat 108. Chen JM , Xi JY , Wang D , et al. Carrier mobility in graphyne should be even larger than that in graphene: a theoretical prediction , J Phys Chem Lett , 2013 , vol. 4 (pg. 1443 - 8 ) Google Scholar Crossref Search ADS PubMed WorldCat 109. Xi JY , Wang D , Yi YP , et al. Electron-phonon couplings and carrier mobility in graphynes sheet calculated using the Wannier-interpolation approach , J Chem Phys , 2014 , vol. 141 pg. 034704 Google Scholar Crossref Search ADS PubMed WorldCat 110. Padilha JE , Fazzio A , da Silva AJR . Directional control of the electronic and transport properties of graphynes , J Phys Chem C , 2014 , vol. 118 (pg. 18793 - 8 ) Google Scholar Crossref Search ADS WorldCat 111. Ivanovskii AL . Graphynes and graphdyines , Prog Solid State Chem , 2013 , vol. 41 (pg. 1 - 19 ) Google Scholar Crossref Search ADS WorldCat 112. Li GX , Li YL , Liu HB , et al. Architecture of graphdiyne nanoscale films , Chem Commun , 2010 , vol. 46 (pg. 3256 - 8 ) Google Scholar Crossref Search ADS WorldCat 113. Gao HY , Franke JH , Wagner H , et al. Effect of metal surfaces in on-surface glaser coupling , J Phys Chem C , 2013 , vol. 117 (pg. 18595 - 602 ) Google Scholar Crossref Search ADS WorldCat 114. Zhang YQ , Kepcija N , Kleinschrodt M , et al. Homo-coupling of terminal alkynes on a noble metal surface , Nat Commun , 2012 , vol. 3 pg. 1286 Google Scholar Crossref Search ADS PubMed WorldCat 115. Wang JY , Zhang SH , Zhou JY , et al. Identifying sp-sp2 carbon materials by Raman and infrared spectroscopies , Phys Chem Chem Phys , 2014 , vol. 16 (pg. 11303 - 9 ) Google Scholar Crossref Search ADS PubMed WorldCat 116. Liu Y , Wang G , Huang QS , et al. Structural and electronic properties of t graphene: a two-dimensional carbon allotrope with tetrarings , Phys Rev Lett , 2012 , vol. 108 pg. 225505 Google Scholar Crossref Search ADS PubMed WorldCat 117. Huang HQ , Li YC , Liu ZR , et al. Comment on ‘Structural and electronic properties of t graphene: a two-dimensional carbon allotrope with tetrarings’ , Phys Rev Lett , 2013 , vol. 110 pg. 029603 Google Scholar Crossref Search ADS PubMed WorldCat 118. De S , Willand A , Amsler M , et al. Energy landscape of fullerene materials: a comparison of boron to boron nitride and carbon , Phys Rev Lett , 2011 , vol. 106 pg. 225502 Google Scholar Crossref Search ADS PubMed WorldCat 119. Wu X , Dai J , Zhao Y , et al. Two-dimensional boron monolayer sheets , ACS Nano , 2012 , vol. 6 (pg. 7443 - 53 ) Google Scholar Crossref Search ADS PubMed WorldCat 120. Tang H , Ismail-Beigi S . Novel precursors for boron nanotubes: the competition of two-center and three-center bonding in boron sheets , Phys Rev Lett , 2007 , vol. 99 pg. 115501 Google Scholar Crossref Search ADS PubMed WorldCat 121. Liu ZR , Wang JY , Li JL . Dirac cones in two-dimensional systems: from hexagonal to square lattices , Phys Chem Chem Phys , 2013 , vol. 15 (pg. 18855 - 62 ) Google Scholar Crossref Search ADS PubMed WorldCat 122. Banerjee S , Singh RRP , Pardo V , et al. Tight-binding modeling and low-energy behavior of the semi-Dirac point , Phys Rev Lett , 2009 , vol. 103 pg. 016402 Google Scholar Crossref Search ADS PubMed WorldCat 123. Wu Y . A semi-Dirac point and an electromagnetic topological transition in a dielectric photonic crystal , Opt Express , 2014 , vol. 22 (pg. 1906 - 17 ) Google Scholar Crossref Search ADS PubMed WorldCat 124. Bellec M , Kuhl U , Montambaux G , et al. Topological transition of Dirac points in a microwave experiment , Phys Rev Lett , 2013 , vol. 110 pg. 033902 Google Scholar Crossref Search ADS PubMed WorldCat 125. Li X , Zhang F , Niu Q , et al. Superlattice valley engineering for designer topological insulators , Sci. Rep. , 2014 , vol. 4 pg. 6397 Google Scholar Crossref Search ADS PubMed WorldCat 126. Tajima N , Kajita K . Experimental study of organic zero-gap conductor alpha-(BEDT-TTF)(2)I-3 , Sci Technol Adv Mat , 2009 , vol. 10 pg. 024308 Google Scholar Crossref Search ADS WorldCat 127. Tajima N , Sugawara S , Kato R , et al. Effect of the zero-mode Landau level on interlayer magnetoresistance in multilayer massless Dirac fermion systems , Phys Rev Lett , 2009 , vol. 102 pg. 176403 Google Scholar Crossref Search ADS PubMed WorldCat 128. Urban DF , Bercioux D , Wimmer M , et al. Barrier transmission of Dirac-like pseudospin-one particles , Phys Rev B , 2011 , vol. 84 pg. 115136 Google Scholar Crossref Search ADS WorldCat 129. Green D , Santos L , Chamon C . Isolated flat bands and spin-1 conical bands in two-dimensional lattices , Phys Rev B , 2010 , vol. 82 pg. 075104 Google Scholar Crossref Search ADS WorldCat 130. Mak KF , Lee C , Hone J , et al. Atomically thin MoS2: a new direct-gap semiconductor , Phys Rev Lett , 2010 , vol. 105 pg. 136805 Google Scholar Crossref Search ADS PubMed WorldCat 131. Raoux A , Morigi M , Fuchs JN , et al. From Dia- to paramagnetic orbital susceptibility of massless fermions , Phys Rev Lett , 2014 , vol. 112 pg. 026402 Google Scholar Crossref Search ADS PubMed WorldCat 132. Polini M , Guinea F , Lewenstein M , et al. Artificial honeycomb lattices for electrons, atoms and photons , Nat Nanotechnol , 2013 , vol. 8 (pg. 625 - 33 ) Google Scholar Crossref Search ADS PubMed WorldCat 133. Wu CJ , Bergman D , Balents L , et al. Flat bands and Wigner crystallization in the honeycomb optical lattice , Phys Rev Lett , 2007 , vol. 99 pg. 070401 Google Scholar Crossref Search ADS PubMed WorldCat 134. Wu CJ , Das Sarma S . px,y-orbital counterpart of graphene: cold atoms in the honeycomb optical lattice , Phys Rev B , 2008 , vol. 77 pg. 235107 Google Scholar Crossref Search ADS WorldCat 135. Bercioux D , Urban DF , Grabert H , et al. Massless Dirac-Weyl fermions in a T-3 optical lattice , Phys Rev A , 2009 , vol. 80 pg. 063603 Google Scholar Crossref Search ADS WorldCat 136. Shen R , Shao LB , Wang BG , et al. Single Dirac cone with a flat band touching on line-centered-square optical lattices , Phys Rev B , 2010 , vol. 81 pg. 041410 Google Scholar Crossref Search ADS WorldCat 137. Tarruell L , Greif D , Uehlinger T , et al. Creating, moving and merging Dirac points with a Fermi gas in a tunable honeycomb lattice , Nature , 2012 , vol. 483 pg. 302 Google Scholar Crossref Search ADS PubMed WorldCat 138. Gibertini M , Singha A , Pellegrini V , et al. Engineering artificial graphene in a two-dimensional electron gas , Phys Rev B , 2009 , vol. 79 pg. 241406 Google Scholar Crossref Search ADS WorldCat 139. Park CH , Louie SG . Making massless dirac fermions from a patterned two-dimensional electron gas , Nano Lett , 2009 , vol. 9 (pg. 1793 - 7 ) Google Scholar Crossref Search ADS PubMed WorldCat 140. Nadvornik L , Orlita M , Goncharuk NA , et al. From laterally modulated two-dimensional electron gas towards artificial graphene , New J Phys , 2012 , vol. 14 pg. 053002 Google Scholar Crossref Search ADS WorldCat 141. Parimi PV , Lu WT , Vodo P , et al. Negative refraction and left-handed electromagnetism in microwave photonic crystals , Phys Rev Lett , 2004 , vol. 92 pg. 127401 Google Scholar Crossref Search ADS PubMed WorldCat 142. Peleg O , Bartal G , Freedman B , et al. Conical diffraction and gap solitons in honeycomb photonic lattices , Phys Rev Lett , 2007 , vol. 98 pg. 103901 Google Scholar Crossref Search ADS PubMed WorldCat 143. Zhang LZ , Wang ZF , Du SX , et al. Prediction of a Dirac state in monolayer TiB2 , Phys Rev B , 2014 , vol. 90 pg. 161402 Google Scholar Crossref Search ADS WorldCat 144. Hasegawa Y , Konno R , Nakano H , et al. Zero modes of tight-binding electrons on the honeycomb lattice , Phys Rev B , 2006 , vol. 74 pg. 033413 Google Scholar Crossref Search ADS WorldCat 145. Montambaux G , Piechon F , Fuchs JN , et al. Merging of Dirac points in a two-dimensional crystal , Phys Rev B , 2009 , vol. 80 pg. 153412 Google Scholar Crossref Search ADS WorldCat 146. Montambaux G , Piechon F , Fuchs JN , et al. A universal Hamiltonian for motion and merging of Dirac points in a two-dimensional crystal , Eur Phys J B , 2009 , vol. 72 (pg. 509 - 20 ) Google Scholar Crossref Search ADS WorldCat 147. Delplace P , Gomez-Leon A , Platero G . Merging of Dirac points and Floquet topological transitions in ac-driven graphene , Phys Rev B , 2013 , vol. 88 pg. 245422 Google Scholar Crossref Search ADS WorldCat 148. Kaufmann RM , Khlebnikov S , Wehefritz-Kaufmann B . Singularities, swallowtails and Dirac points. An analysis for families of Hamiltonians and applications to wire networks, especially the Gyroid , Ann Phys-NY , 2012 , vol. 327 (pg. 2865 - 84 ) Google Scholar Crossref Search ADS WorldCat 149. Kariyado T , Hatsugai Y . Symmetry-protected quantization and bulk-edge correspondence of massless Dirac fermions: application to the fermionic Shastry-Sutherland model , Phys Rev B , 2013 , vol. 88 pg. 245126 Google Scholar Crossref Search ADS WorldCat 150. Hatsugai Y . Topological aspect of graphene physics , J Phys: Conf Ser , 2011 , vol. 334 pg. 012004 Google Scholar Crossref Search ADS WorldCat 151. Landau LD , Lifshitz EM . , Quantum Mechanics: Non-Relativistic Theory (Third Edition) , 1981 Oxford Butterworth-Heinemann pg. 79 Google Scholar Google Preview OpenURL Placeholder Text WorldCat COPAC 152. Asano K , Hotta C . Designing Dirac points in two-dimensional lattices , Phys Rev B , 2011 , vol. 83 pg. 245125 Google Scholar Crossref Search ADS WorldCat 153. Manes JL , Guinea F , Vozmediano MAH . Existence and topological stability of Fermi points in multilayered graphene , Phys Rev B , 2007 , vol. 75 pg. 155424 Google Scholar Crossref Search ADS WorldCat 154. Wang KF , Graf D , Petrovic C . Quasi-two-dimensional Dirac fermions and quantum magnetoresistance in LaAgBi2 , Phys Rev B , 2013 , vol. 87 pg. 235101 Google Scholar Crossref Search ADS WorldCat 155. Zhao RQ , Wang JY , Yang MM , et al. BN-embedded graphene with a ubiquitous gap opening , J Phys Chem C , 2012 , vol. 116 (pg. 21098 - 103 ) Google Scholar Crossref Search ADS WorldCat © The Author(s) 2015. Published by Oxford University Press on behalf of China Science Publishing & Media Ltd. All rights reserved. For Permissions, please email: journals.permissions@oup.com This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0/), which permits non-commercial reuse, distribution, and reproduction in any medium, provided the original work is properly cited. for commercial re-use, please contact journals.permissions@oup.com © The Author(s) 2015. Published by Oxford University Press on behalf of China Science Publishing & Media Ltd. All rights reserved. For Permissions, please email: journals.permissions@oup.com TI - The rare two-dimensional materials with Dirac cones JO - National Science Review DO - 10.1093/nsr/nwu080 DA - 2015-03-01 UR - https://www.deepdyve.com/lp/oxford-university-press/the-rare-two-dimensional-materials-with-dirac-cones-49vDK0NCEN SP - 22 EP - 39 VL - 2 IS - 1 DP - DeepDyve ER -