Abstract
A remarkable manifestation of the quantum character of electrons in matter is offered by graphene, a single atomic layer of graphite. Unlike conventional solids where electrons are described with the Schrödinger equation, electronic excitations in graphene are governed by the Dirac hamiltonian1. Some of the intriguing electronic properties of graphene, such as massless Dirac quasiparticles with linear energy–momentum dispersion, have been confirmed by recent observations2,3,4,5. Here, we report an infrared spectromicroscopy study of charge dynamics in graphene integrated in gated devices. Our measurements verify the expected characteristics of graphene and, owing to the previously unattainable accuracy of infrared experiments, also uncover significant departures of the quasiparticle dynamics from predictions made for Dirac fermions in idealized, free-standing graphene. Several observations reported here indicate the relevance of many-body interactions to the electromagnetic response of graphene.
Similar content being viewed by others
Main
We investigated the reflectance R(ω) and transmission T(ω) of graphene samples on a SiO2/Si substrate (inset of Fig. 1a) as a function of gate voltage Vg at 45 K (see the Methods section). We start with data taken at the charge-neutrality point VCN: the gate voltage corresponding to the minimum d.c. conductivity and zero total charge density (inset of Fig. 1c). Figure 1a shows R(ω) of a graphene gated structure (graphene/SiO2/Si) at VCN=3 V normalized by reflectance of the substrate Rsub(ω). Rsub(ω) is dominated by a minimum around 5,500 cm−1 due to interference effects in SiO2. A remarkable observation is that a monolayer of undoped graphene markedly modifies the interference minimum of the substrate leading to a suppression of Rsub(ω) by as much as 15%. This observation is significant because it enables us to evaluate the conductivity of graphene near the interference structure, as will be discussed below.
Both reflectance and transmission spectra of graphene structures can be modified by a gate voltage. Figure 1b,c shows these modifications at various gate voltages normalized by data at VCN: R(V)/R(VCN) and T(V)/T(VCN), where V =Vg–VCN. These data correspond to the Fermi energy EF on the electron side and similar behaviour was observed with EF on the hole side (not shown). At low voltages (<17 V), we found a dip in R(V)/R(VCN) spectra. With increasing bias, this feature evolves into a peak–dip structure and systematically shifts to higher frequency. The T(V)/T(VCN) spectra reveal a peak at all voltages, which systematically hardens with increasing bias. A voltage-induced increase in transmission (T(V)/T(VCN)>1) signals a decrease of the absorption with bias. Most interestingly, we observed that the frequencies of the main features in R(V)/R(VCN) and T(V)/T(VCN) all evolve approximately as .
To explore the quasiparticle dynamics under applied voltages, it is imperative to first discuss the two-dimensional (2D) optical conductivity of charge-neutral graphene, σ1(ω,VCN)+i σ2(ω,VCN), extracted from a multilayer analysis of the devices (see the Methods section). Theoretical analysis6,7,8 predicts a constant ‘universal’ 2D conductivity σ1(ω,VCN)=πe2/2h for ideal undoped graphene. Our R(ω)/Rsub(ω) data are consistent with this prediction. Figure 1a shows a comparison between the experimental R(ω)/Rsub(ω) spectrum and model spectra generated assuming constant σ1(ω,VCN) values. The constant universal conductivity offers a good agreement (within ±15%) with the experimental spectra in the range 4,000–6,500 cm−1. Outside this spectral region, our infrared measurements do not allow us to unambiguously determine the absolute value of σ1(ω,VCN); therefore, the uncertainty of σ1(ω,VCN) increases as shown by the shaded region weighted around the πe2/2h value in Fig. 2a. However, recent infrared studies of graphene revealed a constant conductivity σ1(ω,VCN)=πe2/2h between 2,400 and 24,000 cm−1 (ref. 9 and Mak, K.F. & Heinz, T. 2008 APS March Meeting, Abstract: L29.00006, unpublished). The universal conductivity is only weakly modified in bulk highly ordered pyrolytic graphite10 and extends down to 800 cm−1. Thus, in the following discussion, we will assume σ1(ω,VCN)=πe2/2h throughout the entire range of our data.
Electrostatic doping of graphene introduces two fundamental changes in the optical conductivity σ1(ω,V)+i σ2(ω,V): a strong Drude component formed in the far-infrared with σ1(ω→0)=4–100 πe2/2h accompanied by a shift in the onset of interband transitions at 2EF, as schematically shown in the inset of Fig. 2b. To investigate these effects, we obtained σ1(ω,V)+i σ2(ω,V) (Fig. 2b,c) from voltage-dependent reflectance and transmission spectra (see the Methods section). The key features in the conductivity spectra are independent of uncertainties in σ1(ω,VCN) discussed above. Regardless of the choice of σ1(ω,VCN), under applied biases, we observe a suppression of the conductivity compared with σ1(ω,VCN) and a well-defined threshold structure above which the conductivity recovers the universal value πe2/2h. The energy of the threshold structure systematically increases with voltage, a natural expectation for a transition occurring at 2EF. With a scattering rate 1/τ=30 cm−1 at 71 V independently obtained from transport data, the Drude mode is narrow and confined below the low-ω cutoff of our measurements. We stress that the two voltage-induced transformations of the conductivity, the intraband mode and the onset of interband absorption at 2EF, are interdependent as suggested by our data. Indeed, assuming the intraband component can be described with a simple Drude formula σ1(ω)=σd.c./(1+ω2τ2) using σd.c. and 1/τ obtained from transport measurements, we find that the spectral weight removed from ω<2EF is recovered under the Drude structure, such that the total oscillator strength given by is conserved at any bias with a cutoff frequency Ωc=8,000 cm−1.
Next, we extracted Fermi energy values from the 2EF threshold using two different methods (see the Methods section). We found that the 2EF values (Fig. 3a) are symmetric for biases delivering either holes or electrons to graphene. Moreover, 2EF increases with voltage approximately as (deviations from the square-root law at small biases will be discussed below). Note that EF of Dirac fermions scales with the 2D carrier density N as (refs 2, 3), where vFis the Fermi velocity. In our devices, N=CgV/e, where Cg=115 aF μm−2 is the gate capacitance per unit area. Therefore, the observed dependence of 2EFsubstantiates that graphene samples integrated in gated devices are governed by Dirac quasiparticles.
Interestingly, the 2EF threshold in σ1(ω,V) shows a width of about 1,400 cm−1 that is independent of gate voltage and therefore of carrier density N, irrespective of a sevenfold enhancement of N between 10 and 71 V. This effect is much stronger than the theoretical estimate for thermal smearing of the 2EF feature at 45 K (refs 7, 8), which is about 500 cm−1. A recent theoretical study11 showed that disorder effects and electron–phonon coupling are needed to account for the width of the 2EF threshold in our data. Apart from that, a spatial variation of local EF values observed in graphene on SiO2/Si substrates (ref. 12 and Brar, V. et al., 2008 APS March Meeting, Abstract: U29.00003, unpublished) will inevitably lead to a broadening of the absorption onset at 2EF in σ1(ω), because infrared measurements register the absorption averaged over a large area (a few micrometres in our experiments). The origin of the inhomogeneity of EF in graphene is still an open question12, which needs to be explored using spatially resolved probes such as near-field infrared conductivity studies capable of probing the response of a material with nanometre resolution over a large area13.
Our study has uncovered several new properties of graphene that are beyond the ideal Dirac fermion picture14. First, our study revealed unexpected features of σ1(ω,V) below 2EF. The band structure of ideal graphene implies that the interband transition at 2EF is the lowest electronic excitation in the system apart from the Drude response at ω=0. Therefore, we anticipate finding σ1(ω,V)≈0 up to the 2EF threshold, provided the Drude scattering rate is much smaller than 2EF. This latter condition is fulfilled for all data in Fig. 2, and yet we registered significant conductivity below 2EF (see the Supplementary Information). This result has not been anticipated by theories developed for Dirac fermions6,7,8. Both extrinsic and intrinsic effects may give rise to the residual conductivity in Fig. 2. Among the former, charged impurities and unitary scatterers (edge defects, cracks, vacancies and so on) were shown to induce considerable residual conductivity below 2EF (ref. 11). However, the theoretical residual absorption in ref. 11 is systematically suppressed with voltage, whereas this suppression was not observed in our data. In addition, the magnitude of the theoretical residual absorption is smaller compared with experimental values in Fig. 2. Therefore, it is likely that other mechanisms are also responsible for the residual conductivity in our data. One intriguing interpretation of the residual conductivity is in terms of many-body interactions, which are known to produce a strong frequency-dependent quasiparticle scattering rate 1/τ(ω). It is predicted theoretically that 1/τ(ω) in graphene increases with frequency owing to electron–electron15,16 and electron–phonon interactions17. The energy-dependent scattering rate initiates a marked enhancement of the conductivity compared with the lorentzian form prescribed by the Drude model. Such an enhancement in mid-infrared frequencies has been observed in many systems18,19,20.
A closer inspection of the voltage dependence of the 2EF feature uncovers marked departures from the behaviour anticipated within a single particle picture of graphene. To highlight these deviations, we plot 2EF as a function of the Fermi wave vector kF based on , as shown in Fig. 3a. The EF(kF) plot has a clear physical meaning: it is a direct representation of the band dispersion. We then examine the ratio EF/ℏkF, which is directly related to the Fermi velocity vF. The EF/ℏkF plot as a function of V1/2 and kF in Fig. 3b reveals a departure from linear dispersion with a single value of EF/ℏkF expected within a single particle picture. Moreover, EF/ℏkF increases systematically with decreasing kF values compared with that at high kF values. The observed systematic enhancement of EF/ℏkF at low kF is indicative of many-body interactions14,21,22. Signatures of band renormalization were also observed in a previous magneto-optical study of graphene4. Importantly, even the smallest EF/ℏkF values in Fig. 3b are higher than that of the bulk graphite23 (∼0.9×106 m S−1), which also supports the hypothesis of EF/ℏkF renormalization in graphene. Complementary information on the EF/ℏkF renormalization in graphene can be obtained from photoemission, which is another potent probe of many-body effects in solids. Currently available photoemission data were all collected for epitaxial graphene grown on SiC (refs 24, 25). This complicates a direct comparison with infrared results for exfoliated samples on SiO2/Si substrates reported here. We conclude by noting that the strong deviations of the experimental electromagnetic response from a simple single particle picture of graphene reported in our study challenge current theoretical conceptions of fundamental properties of this interesting form of carbon and also have implications for its potential applications in opto-electronics.
Note added in proof. After the submission of our paper, we became aware of another work on gate tunable infrared properties of graphene26.
Methods
Sample fabrication and infrared measurements
In the graphene devices studied here, monolayer graphene mechanically cleaved from Kish graphite was deposited onto an infrared transparent SiO2(300 nm)/Si substrate2,3, which also serves as the gate electrode. Then, standard fabrication procedures were used to define multiple Cr/Au (3/35 nm) contacts to the sample. The devices studied here exhibit mobility as high as 8,700 cm2 V−1 s−1 measured at carrier densities of ∼2×1012 cm−2. The characteristic half-integer quantum Hall effect is observed in these samples2,3, confirming the single layer nature of our specimen. Infrared experiments were carried out using an infrared microscope operating with a synchrotron source at the Advanced Light Source in the frequency range of 700–8,000 cm−1. The synchrotron beam is focused in a diffraction-limited spot, which is smaller than the sample. We measured the reflectance R(ω) and transmission T(ω) of the graphene devices as a function of gate voltage Vg with simultaneous monitoring of the d.c. resistivity.
Temperature of the graphene sample
Data reported here were obtained in a micro-cryostat with the sample mounted on a coldfinger in vacuum. The temperature of our graphene sample is warmer than that of the coldfinger, owing to thermal radiation from room-temperature KBr optical windows and electrical isolation of the devices from the coldfinger that compromises thermal contact. A sensor mounted in the immediate proximity to the Si substrate of the devices read T=45 K at the lowest temperatures attainable at the coldfinger. Because both the temperature sensor and the device are in a nearly identical environment, we assumed this reading to be accurate for graphene as well.
Extracting the optical constants of graphene
The graphene device contains four layers: (1) graphene with 2D optical conductivity σ(ω)=σ1(ω)+i σ2(ω), (2) SiO2 gate insulator, (3) Si accumulation layer that forms at the interface of SiO2/Si under the applied bias and (4) Si substrate. Properties of layers 2 and 4 are independent of the gate voltage, whereas layers 1 and 3 are systematically modified by Vg. In our analysis of these multilayer structures, we followed the protocol detailed in ref. 27. Specifically, we carried out reflection, transmission and ellipsometric measurements on the Si substrates and SiO2/Si wafers used in our devices and thus obtained the optical constants of layers 2 and 4. We then investigated infrared properties of test devices Ti/SiO2/Si as a function of gate voltage and thus extracted the optical constants of the Si accumulation layer in wafers used for graphene devices. We find that the response of the Si accumulation layer is confined to far-infrared frequencies28 and gives negligible contribution to mid-infrared data in Fig. 1. Finally, we used a multi-oscillator fitting procedure27 to account for the contribution of σ(ω) of graphene to the reflectance and transmission spectra shown in Fig. 1 using standard methods for multilayered structures.
Extracting Fermi energy EF from conductivity spectra
Because of the broadening of the 2EF threshold in σ1(ω,V), the EFvalues can be determined most accurately from the imaginary part of the optical conductivity spectra σ2(ω,V) shown in Fig. 2c. Indeed, these spectra reveal a sharp minimum at ω=2EF in agreement with a previous theoretical prediction29. The minimum in the σ2(ω,V) spectrum is found from the frequency where the derivative of σ2(ω,V) with respect to frequency is zero. The uncertainties of 2EF obtained from this method are related to the accuracy in defining the minimum in the σ2(ω,V) spectrum. Alternatively, 2EF values can be extracted from the centre frequency of the 2EF threshold in σ1(ω,V). The second method has larger uncertainties, as shown in Fig. 3, due to the ambiguity of defining the centre of the 2EF threshold in σ1(ω,V).
References
Semenoff, G. W. Condensed-matter simulation of a three-dimensional anomaly. Phys. Rev. Lett. 53, 2449–2452 (1984).
Novoselov, K. S. et al. Two-dimensional gas of massless Dirac fermions in graphene. Nature 438, 197–200 (2005).
Zhang, Y., Tan, J. W., Stormer, H. L. & Kim, P. Experimental observation of the quantum Hall effect and Berry’s phase in graphene. Nature 438, 201–204 (2005).
Jiang, Z. et al. Infrared spectroscopy of Landau levels of graphene. Phys. Rev. Lett. 98, 197403 (2007).
Deacon, R. S., Chuang, K.-C., Nicholas, R. J., Novoselov, K. S. & Geim, A. K. Cyclotron resonance study of the electron and hole velocity in graphene monolayers. Phys. Rev. B 76, 081406 (2007).
Ando, T., Zheng, Y. & Suzuura, H. Dynamical conductivity and zero-mode anomaly in honeycomb lattices. J. Phys. Soc. Jpn. 71, 1318–1324 (2002).
Peres, N. M. R., Guinea, F. & Castro Neto, A. H. Electronic properties of disordered two-dimensional carbon. Phys. Rev. B 73, 125411 (2006).
Gusynin, V. P. & Sharapov, S. G. Transport of Dirac quasiparticles in graphene: Hall and optical conductivities. Phys. Rev. B 73, 245411 (2006).
Nair, R. R. et al. Fine structure constant defines visual transparency of graphene. Sciencedoi:10.1126/science.1156965 (2008).
Kuzmenko, A. B., van Heumen, E., Carbone, F. & van der Marel, D. Universal optical conductance of graphite. Phys. Rev. Lett. 100, 117401 (2008).
Peres, N. M. R., Stauber, T. & Castro Neto, A. H. The infrared conductivity of graphene. Preprint at <http://arxiv.org/abs/0803.2816> (2008).
Martin, J. et al. Observation of electron–hole puddles in graphene using a scanning single-electron transistor. Nature Phys. 4, 144–148 (2008).
Qazilbash, M. M. et al. Mott transition in VO2 revealed by infrared spectroscopy and nano-imaging. Science 318, 1750–1753 (2007).
Castro Neto, A. H., Guinea, F., Peres, N. M. R., Novoselov, K. S. & Geim, A. K. The electronic properties of graphene. Preprint at <http://arxiv.org/abs/0709.1163> (2007).
González, J., Guinea, F. & Vozmediano, M. A. H. Unconventional quasiparticle lifetime in graphite. Phys. Rev. Lett. 77, 3589–3592 (1996).
Hwang, E. H., Hu, B. Y.-K. & Das Sarma, S. Inelastic carrier lifetime in graphene. Phys. Rev. B 76, 115434 (2007).
Park, C.-H., Giustino, F., Cohen, M. L. & Louie, S. G. Velocity renormalization and carrier lifetime in graphene from the electron–phonon interaction. Phys. Rev. Lett. 99, 086804 (2007).
Basov, D. N., Singley, E. J. & Dordevic, S. V. Sum rules and electrodynamics of high-Tc cuprates in the pseudogap state. Phys. Rev. B 65, 054516 (2002).
Degiorgi, L. The electrodynamic response of heavy-electron compounds. Rev. Mod. Phys. 71, 687–734 (1999).
Basov, D. N. & Timusk, T. Electrodynamics of high-Tc superconductors. Rev. Mod. Phys. 77, 721–779 (2005).
Gonzalez, J., Guinea, F.. & Vozmediano, M. A. H. Marginal-Fermi-liquid behavior from two-dimensional Coulomb interaction. Phys. Rev. B 59, R2474–R2477 (1999).
Das Sarma, S., Hwang, E. H. & Tse, W.-K. Many-body interaction effects in doped and undoped graphene: Fermi liquid versus non-Fermi liquid. Phys. Rev. B 75, 121406 (2007).
Zhou, S. Y. et al. First direct observation of Dirac fermions in graphite. Nature Phys. 2, 595–599 (2006).
Bostwick, A., Ohta, T., Seyller, T., Horn, K. & Rotenberg, E. Quasiparticle dynamics in graphene. Nature Phys. 3, 36–40 (2006).
Zhou, S. Y. et al. Substrate-induced bandgap opening in epitaxial graphene. Nature Mater. 6, 770–775 (2007).
Wang, F. et al. Gate-variable optical transitions in graphene. Science 320, 206–209 (2008).
Li, Z. Q. et al. Light quasiparticles dominate electronic transport in molecular crystal field-effect transistors. Phys. Rev. Lett. 99, 016403 (2007).
Sai, N., Li, Z. Q., Martin, M. C., Basov, D. N. & Di Ventra, M. Electronic excitations and metal–insulator transition in poly(3-hexylthiophene) organic field-effect transistors. Phys. Rev. B 75, 045307 (2007).
Mikhailov, S. A. & Ziegler, K. New electromagnetic mode in graphene. Phys. Rev. Lett. 99, 016803 (2007).
Acknowledgements
Work at UCSD is supported by the DOE (No. DE-FG02-00ER45799). Research at Columbia University is supported by the DOE (No. DE-AIO2-04ER46133 and No. DE-FG02-05ER46215), NSF (No. DMR-03-52738 and No. CHE-0117752), NYSTAR and the Keck Foundation. The Advanced Light Source is supported by the Director, Office of Science, Office of Basic Energy Sciences, of the US Department of Energy under Contract No. DE-AC02-05CH11231.
Author information
Authors and Affiliations
Corresponding author
Supplementary information
Supplementary Information
Supplementary Information and Supplementary Figure 1 (PDF 139 kb)
Rights and permissions
About this article
Cite this article
Li, Z., Henriksen, E., Jiang, Z. et al. Dirac charge dynamics in graphene by infrared spectroscopy. Nature Phys 4, 532–535 (2008). https://doi.org/10.1038/nphys989
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/nphys989
This article is cited by
-
Hybrid Metal Halide Perovskite–Graphene Sensor for Ultra-sensitive X-ray Detection
Arabian Journal for Science and Engineering (2024)
-
Highly directional and carrier density-independent plasmons in quasi-one-dimensional electron gas systems
Communications Physics (2023)
-
Enhancing third harmonic generation using a mid-infrared graphene plasmonic waveguide
Optical and Quantum Electronics (2023)
-
Nanophotonic modulator based on Silicon-ITO heterojunction and slot waveguide with 2D-graphene sheet
Journal of Optics (2023)
-
A numerical investigation study on tunable graphene-squared pixel array-based infrared polarizer
Applied Physics B (2023)