Abstract
Originating from relativistic quantum field theory, Dirac fermions have been invoked recently to explain various peculiar phenomena in condensed-matter physics, including the novel quantum Hall effect in graphene1,2, the magnetic-field-driven metalâinsulator-like transition in graphite3,4, superfluidity in 3He (ref. 5) and the exotic pseudogap phase of high-temperature superconductors6,7. Despite their proposed key role in those systems, direct experimental evidence of Dirac fermions has been limited. Here, we report the first direct observation of relativistic Dirac fermions with linear dispersion near the Brillouin zone (BZ) corner H, which coexist with quasiparticles that have a parabolic dispersion near another BZ corner K. In addition, we also report a large electron pocket that we attribute to defect-induced localized states. Thus, graphite presents a system in which massless Dirac fermions, quasiparticles with finite effective mass and defect states all contribute to the low-energy electronic dynamics.
Similar content being viewed by others
Main
For most condensed-matter systems, the physics is formulated in terms of the non-relativistic Schrödinger equation, and the low-energy excitations are quasiparticles with finite effective mass. For some special systems, for example, graphene/graphite, where the dispersion is expected to be linear near the Fermi energy, EF, and only touch EF at one point (Dirac point), the physics is described by the relativistic Dirac equation with the speed of light replaced by the Fermi velocity vF. The low-energy excitations in this case are Dirac fermions, which have zero effective mass and a vanishing density of states at the Dirac point. The Dirac fermions are proposed to be responsible for various anomalous phenomena observed in these systems1,2,3,8. So far, transport experiments in graphene have shown results in agreement with the presence of Dirac fermions1,2. Phase analysis of quantum oscillations in graphite has also suggested coexistence of both Dirac fermions and quasiparticles with finite effective mass9. Here, we report the first direct observation of massless Dirac fermions coexisting with quasiparticles with finite effective mass in graphite, by using angle-resolved photoemission spectroscopy (ARPES). ARPES provides the advantage of directly probing the electronic structure with both energy and momentum information, not accessible by any other measurement.
Figure 1a shows an ARPES intensity map measured near the Brillouin zone (BZ) corner H. The out-of-plane momentum, kz, is 0.5c*, where c* is the reciprocal lattice constant (see the Methods section for the extraction of kz values). Following the maximum intensity in this map, a linear Î-shaped dispersion can be clearly observed. The dispersion can be better extracted by following the peak positions in the momentum distribution curves (MDCs), momentum scans at constant energies, shown in Fig. 1b. Here, two peaks in the MDCs disperse linearly and merge near EF. The Fermi velocity extracted from the dispersion is 0.91±0.15Ã106âmâsâ1, similar to a value 1.1Ã106âmâsâ1 reported by a magnetoresistance study of graphene2. We note that at low energy near EF, this linear dispersion is also observed along other cuts through the H point, with similar Fermi velocity. This linear and isotropic dispersion is in agreement with the behaviour of Dirac fermions.
a, ARPES intensity map taken near the H point (photon energy h ν=140âeV, kzâ0.50c*), along a cut through H and perpendicular to kx (see the red line in the BZ shown in c). The inset shows a schematic diagram of the Dirac cone dispersion near EF in the three-dimensional Eâkxâky space. b, MDCs from EF to â2.0âeV. The MDCs are normalized to have the same amplitude and are displaced by the same amount so that the dispersion can be directly viewed by following the peak positions at each energy. The dotted lines are guides to the eyes for the linear-dispersing peaks in the MDCs. c, Three-dimensional BZ for graphite with high-symmetry directions relevant for this paper highlighted with red, green and blue lines.
Another way of probing the linear and isotropic dispersion is to study the intensity maps at constant energy. At EF (Fig. 2a), the intensity map shows a small object near H. The details of this small object will be discussed later. With increasing binding energy, this object expands and shows a circular shape (Fig. 2b,c). We note that only the circular shape in the first BZ is clearly observed. This is attributed to the dipole matrix element10, which suppresses or enhances the intensity in different BZs. However, taking the three-fold symmetry of the sample, this circular shape in the first BZ is expected to extend to other BZs and thus the electronic structure is isotropic near H. As the energy changes to â0.9âeV, the constant energy map deviates slightly from the circular shape (see the arrow in Fig. 2d). This deviation increases with binding energy and a trigonal distortion is clearly observed at â1.2âeV (Fig. 2e). This trigonal distortion is determined by the relevant tight-binding parameters for graphite and further studies to analyse this trigonal distortion are in progress. Overall, Fig. 2 shows that from EF to â0.6âeV, the electronic structure is isotropic in the kxâky plane. Similarly, the Fermi velocity measured within the first BZ is 1.0Ã106âmâsâ1 with a â¤10% variation along different directions, consistent with the circular constant energy maps shown here. Combining the results of Figs 1 and 2, we conclude that from EF to â0.6âeV, the dispersion shows a cone-like behaviour near each BZ corner H, similar to that expected for graphene (Fig. 2f).
aâe, ARPES intensity maps near H (h ν=140âeV, kzâ0.50c*) taken at energies from EF to â1.2âeV. The circles are guides for the circular intensity pattern near the H point. The arrows in d and e point to deviation from the circle. f, Schematic diagram of the dispersion for graphene near six BZ corners in the three-dimensional Eâkxâky space.
To resolve the details of the low-energy dispersion and the small object at EF (Fig. 2a), we show in Fig. 3a an intensity map measured near H with lower photon energy to give better energy and momentum resolution. In the intensity map, two bands dispersing linearly towards EF can be distinguished, as is also clear in the MDCs (Fig. 3b) where two peaks can be observed for all binding energies. The extracted dispersion (open circles in Fig. 3a) from MDCs (Fig. 3b) shows two bands dispersing linearly towards EF, with a minimum separation of 0.020±0.004âà â1 at EF. This linear dispersion near the H point, as well as the isotropic electronic structure shown in Fig. 2 from EF to â0.6âeV, is a basic characteristic of Dirac quasiparticles, which points to the presence of Dirac quasiparticles in the low-energy excitations near the H point in graphite. Furthermore, from the extracted dispersions, the Dirac point is extrapolated to be 50±20âmeV above EF, and thus the small object observed at EF is a hole pocket, in agreement with previous studies of the three-dimensional band structure of graphite9,11. Assuming an ellipsoidal shape for the hole pocket12, the hole concentration is estimated to be 3.1±1.3Ã1018âcmâ3, from the 0.020âà â1 separation of the peaks at EF. This hole concentration is in agreement with reported values12,13,14. We note that, given the current resolution of the ARPES technique, we are not able to resolve the two hole pockets at the H point reported by the other experimental probe15. In fact, according to ref. 15, the difference in energy between these two hole pockets at the H point is â¤1âmeV, which is beyond the current resolution of the ARPES technique. The presence of holes with Dirac fermion dispersion is further supported by the angle-integrated intensity (Fig. 3c), which is proportional to the two-dimensional density of states, barring the matrix element. In this energy range, a linear behaviour, similar to that expected for Dirac fermions, is observed. In addition, the energy intersect is at â50âmeV above EF, in agreement with the Dirac point energy extrapolated from the dispersions.
a, ARPES intensity map near the H point (h ν=65âeV, kzâ0.45c*) along the AHLâ² direction (green line in the BZ shown in Fig. 1c). The open blue circles show the MDC dispersion. The dotted lines are guides for the linear dispersion. b, MDCs at energies from EF to â0.2âeV for the data shown in a. Note that, similar to Fig. 2, the intensity of the Ï band is strongly enhanced in the first BZ (HâA direction), due to the dipole matrix element10. c, Intensity obtained by integrating over both kx and ky for data taken near H (h ν=140âeV, kzâ0.50c*). The intensity has been symmetrized with respect to the Dirac point energy (EDâ50âmeV) to compare directly with the expected density of states (DOS) for Dirac fermions (inset). An overall linear behaviour is observed with some weak additional intensity around 100 meV from the Dirac point energy. The origin of this additional weak intensity is unclear and needs further investigation.
Figures 1â3 show that near the H point, the low-energy excitations in graphite are Dirac fermions characterized by linear and isotropic cone-like dispersion, in agreement with transport measurements in graphite where Dirac fermions are suggested to coexist with quasiparticles that have finite effective mass9. To gain direct insight into the different types of quasiparticles, ARPES can provide a unique advantage by directly measuring the effective mass as well as accessing its momentum dependence. Figure 4a shows the intensity map near another high-symmetry point in the BZ corner, the K point. The dispersion (open circles) shows a parabolic behaviour, in sharp contrast to the linear behaviour observed near the H point (Figs 1, 3). This parabolic dispersion points to the presence of quasiparticles with finite effective mass. To determine the effective mass, we first extract the low-energy dispersion, then fit the MDC and energy distribution curve (EDC) dispersions (Fig. 3b,c) with a parabolic function. In both cases, the effective mass is determined to be 0.069±0.015me, where me is the free-electron mass. This effective mass measured by ARPES is different from the values reported by transport measurements13,14,16, 0.052 and 0.038me for electrons and holes respectively. This difference may be due to the fact that transport measurements are not momentum selective, and therefore the mass measured is the average mass over all kz values. On the other hand, ARPES is momentum selective and the value for the effective mass is only for this specific kz value. Taking this into account, the agreement between these measurements is reasonable. In summary, the data presented so far show that the low-energy excitations in graphite change from massless Dirac fermions with linear dispersion near H (Figs 1â3) to quasiparticles with parabolic dispersion and finite effective mass near K (Fig. 4a).
a, ARPES intensity map near K (h ν=50âeV, kzâ0.08c*) along the ÎKMâ² direction (blue line in the BZ shown in Fig. 1c). The open circles are the dispersions extracted from MDCs. b, MDCs from EF to â50âmeV for the data in a. The open circles mark the peaks clearly resolved in the data. The inset shows the MDC dispersion from â10 to â50âmeV, with the parabolic fit used to extract the effective mass. c, EDCs from k0 to k12, as indicated in e. The open circles are the peak positions for the large electron pocket. d, Closer view of the data shown in a with the EDC dispersion (dotted line). e, Intensity map near the K point measured in some parts of the sample, which shows an additional large electron pocket. The open circles are dispersions extracted from EDCs shown in c. f, MDC at EF from the data shown in e. The black arrows point to the peaks from the large electron pocket, which are separated by â0.1âà â1, whereas the grey arrow points to the peak from the Ï band.
We now discuss another interesting feature observed in graphite, that is, a large electron pocket near EF. Figure 4e shows the intensity map measured in the same experimental conditions as Fig. 4d except at a different spot. In Fig. 4e, a strong and large electron pocket within 50âmeV below EF is the dominant feature. Weak signatures of the parabolic Ï band (as in Fig. 4d) can still be observed, as the MDC at EF (Fig. 4f) demonstrates. Here, in addition to the two main peaks (black arrows) corresponding to the electron pocket, a central weak peak (grey arrow) corresponding to the top of the parabolic band can also be distinguished. From the separation (â0.1âà â1) between the two main peaks at EF, the electron concentration is determined to be 8.0±0.7Ã1019âcmâ3, which is an order of magnitude higher than the values reported13,14. Moreover, from the dispersion (Fig. 4e), the effective mass is extracted to be 0.42±0.07me, which is also much larger than any mass reported by transport measurements13,14,16.
This large electron pocket is observed in most of the samples measured, and thus it represents an important feature associated with graphite. We note that a similar large electron pocket has been reported recently and proposed to be associated with either edge states or dangling bonds12. Here we provide detailed characterization of this large electron pocket, which is important in revealing its origin. We propose defect-induced localized states as a possible explanation for this large electron pocket, on the basis of the following reasons. First, the electron concentration and effective mass measured for this electron pocket are much larger than reported values. Second, although the parabolic Ï band in Fig. 4a,d is observed in all the samples measured and in different spots (averaged over â100âμm) within the same sample, this large electron pocket strongly depends on the spot position within the same sample. Third, scanning tunnelling microscopy shows that zigzag edges can induce a peak in the local density of states at an energy (ââ0.03âeV) similar to the electron pocket observed here17. In fact, it has been shown that a low concentration of defects (for example, edge states, vacancies and so on) can induce self-doping to the sample18,19. If this interpretation is correct, then further studies on this large electron pocket may provide insights into the magnetic properties of nano-graphite ribbons, because it has been proposed that some defect-induced localized states are magnetic20,21.
Methods
High-resolution ARPES data have been taken at Beamlines 12.0.1 and 7.0.1 of the Advanced Light Source in Lawrence Berkeley National Laboratory with photon energies from 20 to 155âeV with energy resolution from 15 to 65âmeV. Data were taken from single-crystal graphite (Kish graphite for Fig. 3a,b and natural graphite for other figures) at a temperature of 25âK. Throughout this paper, the kz values are estimated using the standard free-electron approximation of the ARPES final state22,23,24. The inner potential needed for extracting the kz value was determined from the periodicity of the detailed dispersion at the Brillouin zone centre ÎA using a wide range of photon energies from 34 to 155âeV. The consistency of the kz values is confirmed by the degeneracy of the Ï bands near H and a maximum splitting near K (Fig. 5), both in agreement with band structure11,25,26. Note that although this large energy range is good for discussing the splitting of the Ï bands, it cannot capture the deviation from linear behaviour at the K point expected at low energy (as discussed in Fig. 4).
References
Zhang, Y. B., Tan, Y.-W., Stormer, H. L. & Kim, P. Experimental observation of the quantum Hall effect and Berryâs phase in graphene. Nature 438, 201â204 (2005).
Novoselov, K. S. et al. Two-dimensional gas of Dirac fermions in graphene. Nature 438, 197â200 (2005).
Kopelevich, Y., Torres, J. H. S. & da Silva, R. R. Reentrant metallic behavior of graphite in the quantum limit. Phys. Rev. Lett. 90, 156402 (2003).
Du, X., Tsai, S.-W., Maslov, D. L. & Hebard, A. F. Metal-insulator-like behavior in semimetallic bismuth and graphite. Phys. Rev. Lett. 94, 166601 (2005).
Volovik, G. E. Field theory in superfluid 3He: what are the lessons for particle physics, gravity and high-temperature superconductivity? Proc. Natl Acad. Sci. 96, 6042â6047 (1999).
Rantner, W. & Wen, X.-G. Electron spectral function and algebraic spin liquid for the normal state of underdoped high Tc superconductors. Phys. Rev. Lett. 86, 003871 (2001).
Franz, M. & Tesanovic, Z. Algebraic Fermi liquid from phase fluctuations: âtopologicalâ Fermions, votex âBerryonsâ, and QED3 theory of cuprate superconductors. Phys. Rev. Lett. 87, 257003 (2001).
González, J., Guinea, F. & Vozmediano, M. A. H. Unconventional quasiparticle lifetime in graphite. Phys. Rev. Lett. 77, 003589 (1996).
Lukâyanchuk, I. A. & Kopelevich, Y. Phase analysis of quantum oscillations in graphite. Phys. Rev. Lett. 93, 166402 (2004).
Shirley, E. L., Terminello, L. J., Santoni, A. & Himpsel, F. J. Brillouin-zone-selection effects in graphite photoelectron angular distributions. Phys. Rev. B 51, 013614 (1995).
Dresselhaus, M. S. & Dresselhaus, G. Intercalation compounds of graphite. Adv. Phys. 51, 1â186 (2002).
Sugawara, K., Sato, T., Souma, S., Takahashi, T. & Suematsu, H. Fermi surface and edge-localized states in graphite studied by high-resolution angle-resolved photoemission spectroscopy. Phys. Rev. B 73, 045124 (2006).
Soule, D. E. Magnetic field dependence of the Hall effect and magnetoresistance in graphite single crystals. Phys. Rev. 112, 698â707 (1958).
Zhang, Y. B., Small, J. P., Amori, M. E. S. & Kim, P. Electric field modulation of galvanomagnetic properties of mesoscopic graphite. Phys. Rev. Lett. 94, 176803 (2005).
Toy, W. W., Dresselhaus, M. S. & Dresselhaus, G. Minority carriers in graphite and the H-point magnetoreflection spectra. Phys. Rev. B 15, 4077â4090 (1977).
Galt, J. K., Yager, W. A. & Dail, H. W. Jr. Cyclotron resonance effects in graphite. Phys. Rev. 103, 1586â1587 (1956).
Kobayashi, Y., Fukui, K.-I., Enoki, T., Kusakabe, K. & Kaburagi, Y. Observation of zigzag and armchair edges of graphite using scanning tunneling microscopy and spectroscopy. Phys. Rev. B 71, 193406 (2005).
Peres, N. M. R., Guinea, F. & Castro Neto, A. H. Electronic properties of disordered two-dimensional carbon. Phys. Rev. B 73, 125411 (2006).
Pereira, V. M., Guinea, F., Lopes Dos Santos, J. M. B., Peres, N. M. R. & Castro Neto, A. H. Disorder induced localized states in graphene. Phys. Rev. Lett. 96, 036801 (2006).
Nakada, K., Fujita, M., Dresselhaus, G. & Dresselhaus, M. S. Edge state in graphene ribbons: nanometer size effect and edge shape dependence. Phys. Rev. B 54, 017954 (1996).
Wakabayashi, K., Fujita, M., Ajiki, H. & Sigrist, M. Electronic and magnetic properties of nanographite ribbons. Phys. Rev. B 59, 008271 (1999).
Law, A. R., Johnson, M. T. & Hughes, H. P. Synchrontron-radiation-excited angle-resolved photoemission from single-crystal graphite. Phys. Rev. B 34, 004289 (1986).
Himpsel, F. J. Angle-resolved measurements of the photoemission of electrons in the study of solids. Adv. Phys. 32, 1â51 (1983).
Hüfner, S. Photoelectron Spectroscopy (Springer, Berlin, 1995).
Charlier, J.-C., Gonze, X. & Michenaud, J.-P. First principles study of the electronic properties of graphite. Phys. Rev. B 43, 004579 (1991).
Nilsson, J., Castro Neto, A. H., Guinea, F. & Peres, N. M. R. Electronic properties of graphene multilayers. Preprint at <http://www.arxiv.org/abs/cond-mat/0604106> (2006).
Acknowledgements
We thank A. Castro Neto, V. Oganesyan, A. Bill, K. McElroy, C. M. Jozwiak and D. Garcia for useful discussions and E. Domning and B. Smith for beam line 12.0.1 control software. This work was supported by the National Science Foundation through Grant No. DMR03-49361, the Director, Office of Science, Office of Basic Energy Sciences, Division of Materials Sciences and Engineering of the US Department of Energy under Contract No. DEAC03-76SF00098 and by the Laboratory Directed Research and Development Program of Lawrence Berkeley National Laboratory under the Department of Energy Contract No. DE-AC02-05CH11231.
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Rights and permissions
About this article
Cite this article
Zhou, S., Gweon, GH., Graf, J. et al. First direct observation of Dirac fermions in graphite. Nature Phys 2, 595â599 (2006). https://doi.org/10.1038/nphys393
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/nphys393
This article is cited by
-
Enhanced optical conductivity and many-body effects in strongly-driven photo-excited semi-metallic graphite
Nature Communications (2023)
-
Nanophotonic modulator based on Silicon-ITO heterojunction and slot waveguide with 2D-graphene sheet
Journal of Optics (2023)
-
A ferromagnetic hybrid Weyl semimetal in two dimensions: the monolayer AgCrS2
Journal of Materials Science (2023)
-
Observation of interaction-induced phenomena of relativistic quantum mechanics
Communications Physics (2021)
-
Computational Study of Electronic and Thermoelectric Properties of ZnO/Graphene Heterostructures
International Journal of Thermophysics (2021)