SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

Open Access Nano Express

Dirac cone in α -graphdiyne: a first-principles study

Xiaoning Niu, Xingze Mao, Dezheng Yang, Zhiya Zhang, Mingsu Si* and Desheng Xue

Author Affiliations

Key Laboratory for Magnetism and Magnetic Materials of the Ministry of Education, Lanzhou University, 222 South Tianshui Road, Lanzhou 730000, China

For all author emails, please log on.

Nanoscale Research Letters 2013, 8:469  doi:10.1186/1556-276X-8-469


The electronic version of this article is the complete one and can be found online at: http://www.nanoscalereslett.com/content/8/1/469


Received:12 September 2013
Accepted:26 October 2013
Published:9 November 2013

© 2013 Niu et al.; licensee Springer.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

We investigate the Dirac cone in α-graphdiyne, which is a predicted flat one-atom-thick allotrope of carbon using first-principles calculations. α-graphdiyne is derived from graphene where two acetylenic linkages (-C ≡C-) are inserted into the single bonds (-C-C-). Thus, α-graphdiyne possesses a larger lattice constant which subsequently affects its electronic properties. Band structures show that α-graphdiyne exhibits similar Dirac points and cone to graphene. Further, the tight-binding method is used to exploit the linear dispersion in the vicinity of Dirac points. Thanks to the larger lattice constant, α-graphdiyne yields a lower Fermi velocity, which might make itself an ideal material to serve the anomalous integer quantum Hall effect.

Keywords:
α-graphdiyne; Dirac cone; First-principles calculation

Background

Band theory was first used to study the band structure of graphene over half a century ago [1], and it demonstrated that graphene is a semimetal with unusual linearly dispersing electronic excitations called Dirac electron. Such linear dispersion is similar to photons which cannot be described by the Schrödinger equation. In the vicinity of the Dirac point where two bands touch each other at the Fermi energy level, the Hamiltonian obeys the two-dimensional (2D) Dirac equation [2] as <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/469/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/469/mathml/M1">View MathML</a> with vF being the Fermi velocity, <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/469/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/469/mathml/M2">View MathML</a> the Pauli matrices, and <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/469/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/469/mathml/M3">View MathML</a> the momentum operator. In graphene, the Fermi velocity vF is 300 times smaller than the speed of light. Hence, many unusual phenomena of quantum electrodynamics can be easily detected because of the much lower speed of carriers [3]. Within the framework of tight-binding approximation, the Fermi velocity vF is proved to be dependent on both the lattice constant and the hopping energy. In fact, the hopping energy is also associated with the lattice constant. Thus, the Fermi velocity of Dirac cone materials might be tunable through changing the corresponding lattice constant.

Recently, it was found that Dirac cones not only occur in the 2D carbon allotropes such as graphene, graphyne, and graphdiyne [4], but also can be detected at interfaces of topological insulators [5-11]. It is notable that, in 6,6,12-graphyne [4], the conduction electrons turn out to be superior to that in graphene in one preferred direction over the other, which is due to the rectangular lattice. This is a major step in searching for new Dirac cone materials. Therefore, it is proper to pursue the Dirac cone material with tunable Fermi velocity, which will be the focus of future researches.

In this letter, we predict a novel flat one-atom-thick allotrope of carbon by inserting two acetylenic linkages into the single bonds in graphene. According to the naming method used in [4], we assign it as α-graphdiyne. Up to now, no study has been made on α-graphdiyne both experimentally and theoretically. Thus, theoretical investigation on α-graphdiyne is a must before synthesizing it in experiments. Since α-graphdiyne has a larger lattice constant, it should have potential applications both in quantum tunneling [12] and in anomalous integer quantum Hall effect [13]. In this work, band structures are calculated and a similar Dirac cone to that of graphene is observed. In particular, we introduce a tight-binding model to mimic the hopping energy between the hexagonal vertices, which realizes the linear dispersion of bands near the Dirac points, allowing the Dirac cone to be studied explicitly.

Methods

To simulate the electronic properties, we employ density functional theory with the generalized gradient approximation (GGA) of Perdew-Burke-Ernzerhof (PBE) [14] for the exchange-correlation (XC) potential within the projector augmented wave method, as implemented in VASP [15]. The cutoff energy for plane waves is set to be 500 eV. The vacuum space is at least 15 Å, which is large enough to avoid the interaction between periodical images; 15 ×15×1 and 25 ×25×1 are used for the k-grid of geometry optimization and self-consistent calculation, respectively. During the geometry optimization, all the atoms in the unit cell were allowed to relax and the convergence of force is set to 0.001 eV/Å.

Results and discussion

Based on first-principles calculation, the lattice structure of α-graphdiyne is predicted for the first time, as shown in Figure 1. It clearly shows that α-graphdiyne has a hexagonal lattice the same as graphene. The optimized lattice constant is 11.42 Å. This is very insightful. On one hand, it has the largest lattice constant compared with currently known carbon allotropes [16] and thus has a much smaller density than graphene and other related carbon allotropes. This makes α-graphdiyne a potential candidate for hydrogen storage [17]. At the same time, the absorbed hydrogen may induce an intrinsic magnetism in the defected system [18,19]. On the other hand, this lattice constant has a very little mismatch with Si(111) surface, that is, one unit cell of α-graphdiyne matches the Si(111)-3 ×3 unit cell well. This suggests that the high possibility is to grow α-graphdiyne epitaxially on Si(111) substrate. After the epitaxial structure is cooled down, one can remove the substrate by chemical etching. In this way, the isolation of monolayer α-graphdiyne might be obtained in experiments.

thumbnailFigure 1. Crystal structure of α-graphdiyne.(a) A unit cell and (b) a 4×4 supercell. (c) A simplified model to mimic the hopping matrix elements along two carbon triple bonds in α-graphdiyne. Carbon atoms 1 and 6 are at vertices of a hexagon in α-graphdiyne. The black balls and blue line represent carbon atoms and the crystalline cell, respectively.

The band structure and density of states (DOS) of α-graphdiyne are shown in Figure 2a,b, respectively. The most important observations from Figure 2a are the linear dispersion near the K point and the zero DOS at the Fermi energy level. However, the corresponding slope of the Dirac cone is obviously smaller than that of graphene and α-graphyne. This has a big effect on the Fermi velocity, as discussed below. The bonding and antibonding orbitals at the Fermi energy level touch each other and develop two slight flat bands as K approaches M, which correspond to the two peaks near the Fermi level in the DOS plot. Similar to the case of graphene and α-graphyne, the Dirac points are located at the K and K, which means that there are even (six) Dirac points in the Brillouin zone, which is in a striking difference from the odd Dirac points observed in topological insulator Bi2Te3[20].

thumbnailFigure 2. Electronic properties of α-graphdiyne.(a) Band structure and (b) DOS. (c) First Brillouin zone with the letters designating high-symmetry points. (d) 2D Dirac cone representing the valence and conduction bands in the vicinity of the K and K points. EF is the Fermi energy.

Due to the breaking symmetry associated with spin-orbit interaction (SOI) in 2D layered materials, a small band gap will be induced at the Dirac points, which can in principle be used to study the quantum spin Hall effect. The energy bands with SOI (not shown for brevity) open a band gap of 22 ×10-3 meV in α-graphdiyne, and the magnitude is close to the value of graphene [21].

To understand the nature of the Dirac cone in α-graphdiyne, we employ the tight-binding method proposed in [22], where an effective hopping parameter is introduced. It is notable that there are six carbon atoms along the effective hopping direction in α-graphdiyne, as shown in Figure 1, while only four in α-graphyne. This makes it more complex to exploit α-graphdiyne than α-graphyne. To simplify the model, two triple bonds with the hopping parameters t1 and t2 for the single and triple carbon bonds are taken. The simplified Hamiltonian equations at the carbon triple bond, i.e., sites 2, 3, 4, and 5, are

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/469/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/469/mathml/M4">View MathML</a>

(1)

where E and V are the electron and on-site energies, respectively. Based on Equation 1, we can easily obtain

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/469/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/469/mathml/M5">View MathML</a>

(2)

When we focus on the electrons near the Dirac cone where E ≈ V, the wave functions are approximately φ2 ≅ -(t1/t2)φ4, φ3 ≅ -(t1/t2)φ1, φ4 ≅ -(t1/t2)φ6, and φ5 ≅ -(t1/t2)φ3. Thus, the hopping term from site 2 to 1 is <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/469/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/469/mathml/M6">View MathML</a>, from site 3 to 4 is <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/469/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/469/mathml/M7">View MathML</a>, from site 4 to 3 is <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/469/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/469/mathml/M8">View MathML</a>, and from site 5 to 6 is <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/469/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/469/mathml/M9">View MathML</a>. With the above four hopping terms, we thus have

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/469/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/469/mathml/M10">View MathML</a>

(3)

which means that the effective direct hopping parameter between sites 1 and 6 is

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/469/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/469/mathml/M11">View MathML</a>

(4)

The obtained effective hopping parameter <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/469/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/469/mathml/M12">View MathML</a> has the same sign as t1, which means that pseudospin in α-graphdiyne has the same direction as in graphene. Thus, many perspectives of graphene can be transferred to α-graphdiyne directly. The magnitude of <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/469/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/469/mathml/M13">View MathML</a> depends on the hopping parameter t2. Remarkably, it equals t1/t2 times the effective hopping parameter in α-graphyne. Thus, the effective hopping parameter should be smaller in α-graphdiyne than in α-graphyne as t1/t2 < 1.

Once we obtain the effective hopping parameter <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/469/mathml/M14','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/469/mathml/M14">View MathML</a>, the standard energy-momentum relation can be obtained directly as [1]

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/469/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/469/mathml/M15">View MathML</a>

(5)

where a is the lattice constant. By fitting the occupied and unoccupied bands in the vicinity of the K point from the first-principles calculations, as illustrated in Figure 2a, the renormalized hopping parameter <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/469/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/469/mathml/M16">View MathML</a> has a value of 0.45 eV. It is much smaller than the value of approximately 3 eV in graphene, which originates from the larger lattice constant in α-graphdiyne. Figure 2c shows the high-symmetry points in the first Brillouin zone. It explicitly shows that the energy bands are degenerate to zero at both K and K points. In Figure 2d, a 2D plot of the Dirac cone of α-graphdiyne is displayed. Due to the same hexagonal lattice as graphene and α-graphyne, the 2D Dirac cone of α-graphdiyne exhibits a similar appearance.

It is known that the Fermi velocity plays a vital role in the photoelectric field and crucially dominates the transport properties. Here, we will focus attention on the study of Fermi velocity of α-graphdiyne. The dispersion close to the K and K points can be expanded as

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/469/mathml/M17','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/469/mathml/M17">View MathML</a>

(6)

where q is the momentum measured relative to the Dirac points, ‘ ±’ the upper and lower Dirac cones, and vF the Fermi velocity, given by <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/469/mathml/M18','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/469/mathml/M18">View MathML</a>. With the lattice constant a = 11.42 Å and the effective hopping parameter <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/469/mathml/M19','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/469/mathml/M19">View MathML</a> = 0.45 eV, the slope of the Dirac cone in α-graphdiyne equals ±4.5 eVÅ compared with ±28 eVÅ in α-graphyne and ±34 eVÅ in graphene [4]. The corresponding Fermi velocity is about 0.11 ×106 m/s, which is much lower than the value in α-graphyne. From this perspective, α-graphdiyne, which has a lower Fermi velocity than other known carbon allotropes, will lead to possible applications in quantum electrodynamics, for example, to observe the anomalous integer quantum Hall effect at room temperature [13].

More information including the helical texture of Dirac cone and Berry’s phase are indeed associated with the detailed wave functions. In this work, we instead calculate the two orbitals at the Dirac point as shown in Figure 3. The charge density explicitly exhibits a 180° rotational symmetry, which is consistent with the theoretical conclusion in the literature [2]. The two orbitals consist of two types of bonds in α-graphdiyne: One is the bonding bonds (Figure 3a) and the other the antibonding bonds (Figure 3b), which are located at the different carbons. As a recent study reported [23], the effective hopping term of the acetylenic linkages is much smaller than the direct hopping between the vertex atoms. This is because the covalent bonds are formed in these acetylenic linkages as illustrated in Figure 3, which subsequently weakens the hopping ability. Thus, the reduced hopping parameter is a natural consequence, which also agrees well with our above tight-binding theory. Future experiments can test this prediction directly.

thumbnailFigure 3. Charge density distributions of two orbitals at the Dirac point. The (a) bonding and (b) antibonding bonds. The isovalues are set to 0.03 Å -3; 3 ×3 supercells are given for the sake of clarity.

Conclusions

In conclusion, we have predicted a novel carbon allotrope called α-graphdiyne, which has a similar Dirac cone to that of graphene. The lower Fermi velocity stems from its largest lattice constant compared with other current carbon allotropes. The effective hopping parameter <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/469/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/469/mathml/M20">View MathML</a> of 0.45 eV is obtained through fitting the energy bands in the vicinity of Dirac points. The obtained Fermi velocity has a lower value of 0.11 ×106 m/s, which might have potential applications in quantum electrodynamics.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

MSS designed the work and revised the paper. XNN calculated the first-principles results. XZM wrote the manuscript. DZY, ZYZ, and DSX have devoted valuable discussion. All authors read and approved the final manuscript.

Acknowledgements

We would like to thank L. Huang (LZU, Lanzhou) for the valuable discussion. This work was supported by the National Basic Research Program of China under no. 2012CB933101, the Fundamental Research Funds for the Central Universities (no. 2022013zrct01), and the National Science Foundation (51202099 and 51372107).

References

  1. Wallace PR: The band theory of graphite.

    Phys Rev 1947, 71:622-634. Publisher Full Text OpenURL

  2. Neto AHC, Guinea F, Peres NMR, Novoselov KS, Geim AK: The electronic properties of graphene.

    Rev Mod Phys 2009, 81:109-162. Publisher Full Text OpenURL

  3. Neto AHC, Guinea F, Peres NMR: Drawing conclusions from graphene.

    Phys World 2006, 19:33-37. OpenURL

  4. Malko D, Neiss C, Vines F, Görling A: Competition for graphene graphynes with direction-dependent dirac cones.

    Phys Rev Lett 2012, 108:086804. PubMed Abstract | Publisher Full Text OpenURL

  5. Fu L, Kane CL, Mele EJ: Topological insulators in three dimensions.

    Phys Rev Lett 2007, 98:106803. PubMed Abstract | Publisher Full Text OpenURL

  6. Takahashi R, Murakami S: Gapless interface states between topological insulators with opposite Dirac velocities.

    Phys Rev Lett 2011, 107:166805. PubMed Abstract | Publisher Full Text OpenURL

  7. Kane CL, Mele EJ: Quantum spin hall effect in graphene.

    Phys Rev Lett 2005, 95:226801. PubMed Abstract | Publisher Full Text OpenURL

  8. Kane CL, Mele EJ: Z2 topological order and the quantum spin hall effect.

    Phys Rev Lett 2005, 95:146802. PubMed Abstract | Publisher Full Text OpenURL

  9. Bernevig BA, Zhang SC: Quantum spin hall effect.

    Phys Rev Lett 2006, 96:106802. PubMed Abstract | Publisher Full Text OpenURL

  10. Moore JE, Balents L: Topological invariants of time-reversal-invariant band structures.

    Phys Rev B 2007, 75:121306(R). OpenURL

  11. Zhang H, Liu CX, Zhang SC: Spin-orbital texture in topological insulators.

    Phys Rev Lett 2013, 111:066801. PubMed Abstract | Publisher Full Text OpenURL

  12. He WY, Chu ZD, He L: Chiral tunneling in a twisted graphene bilayer.

    Phys Rev Lett 2013, 111:066803. PubMed Abstract | Publisher Full Text OpenURL

  13. Novoselov KS, Geim AK, Morozov SV, Jiang D, Katsnelson MI, Grigorieva IV, Dubonos SV, Firsov AA: Two-dimensional gas of massless Dirac fermions in graphene.

    Nature 2005, 438:197-200. PubMed Abstract | Publisher Full Text OpenURL

  14. Perdew JP, Burke K, Ernzerhof M: Generalized gradient approximation made simple.

    Phys Rev Lett 1996, 77:3865-3868. PubMed Abstract | Publisher Full Text OpenURL

  15. Kresse G, Furthmüller J: Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set.

    Phys Rev B 1116, 54:9-11186. OpenURL

  16. Ivanovskii AL: Graphynes and graphdyines.

    Prog Solid State Chem 2013, 41:1-19. Publisher Full Text OpenURL

  17. Koo J, Hwang HJ, Kwon Y, Lee H: Graphynes and graphdyines.

    J Phys Chem C 2012, 116:20220. Publisher Full Text OpenURL

  18. Si MS, Li JY, Shi HG, Niu XN, Xue DS: Divacancies in graphitic boron nitride sheets.

    Europhys Lett 2009, 86:46002. Publisher Full Text OpenURL

  19. Li J, Gao D, Niu X, si M, Xue D: g-B3N3C: a novel two-dimensional graphite-like material.

    Nanoscale Res Lett 2012, 7:624. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  20. Chen YL, Analytis JG, Chu JH, Liu ZK, Mo SK, Qi XL, Zhang HJ, Lu DH, Dai X, Fang Z, Zhang SC, Fisher IR, Hussain Z, Shen ZX: Experimental realization of a three-dimensional topological insulator, Bi2Te3.

    Science 2009, 325:178-181. PubMed Abstract | Publisher Full Text OpenURL

  21. Gmitra M, Konschuh S, Ertler C, Ambrosch-Draxl C, Fabian J: Band-structure topologies of graphene: spin-orbit coupling effects from first principles.

    Phys Rev B 2009, 80:235431. OpenURL

  22. Kim BG, Choi HJ: Graphyne: hexagonal network of carbon with versatile Dirac cones.

    Phys Rev B 2012, 86:115435. OpenURL

  23. Huang H, Duan W, Liu Z: The existence/absence of Dirac cones in graphynes.

    New J Phys 2013, 15:023004. Publisher Full Text OpenURL