Elsevier

Ultramicroscopy

Volume 184, Part A, January 2018, Pages 134-142
Ultramicroscopy

Analytical modeling of electron energy loss spectroscopy of graphene: Ab initio study versus extended hydrodynamic model

https://doi.org/10.1016/j.ultramic.2017.08.014Get rights and content

Highlights

  • We present an analytical modeling of the STEM-EELS data for free-standing graphene.

  • We use an ab initio method and an extended hydrodynamic model with Dirac correction.

  • The extended hydrodynamic model is calibrated to optical limit of the ab initio data.

  • The results of modeling compare favorably with three independent EELS experiments.

  • Some EELS data likely reveal Dirac contribution after subtraction of zero-loss peak.

Abstract

We present an analytical modeling of the electron energy loss (EEL) spectroscopy data for free-standing graphene obtained by scanning transmission electron microscope. The probability density for energy loss of fast electrons traversing graphene under normal incidence is evaluated using an optical approximation based on the conductivity of graphene given in the local, i.e., frequency-dependent form derived by both a two-dimensional, two-fluid extended hydrodynamic (eHD) model and an ab initio method. We compare the results for the real and imaginary parts of the optical conductivity in graphene obtained by these two methods. The calculated probability density is directly compared with the EEL spectra from three independent experiments and we find very good agreement, especially in the case of the eHD model. Furthermore, we point out that the subtraction of the zero-loss peak from the experimental EEL spectra has a strong influence on the analytical model for the EEL spectroscopy data.

Introduction

With its unique electrical and optical properties (relatively low loss, high confinement, mechanical flexibility, and good tunability) graphene is an ideal material for plasmonic applications covering a wide frequency range from terahertz up to infrared, even stretching into the visible regime [1], [2], [3], [4]. Electron energy loss spectroscopy (EELS) is a commonly used experimental technique for investigating electronic and plasmonic properties of materials, including graphene sheets [1]. High-energy single-particle inter-band excitations in graphene, which are often misnomered as π and π + σ plasmons [5], have been studied recently by EELS experiments using high-energy electron beams (∼100 keV) in scanning transmission electron microscope (STEM) on samples consisting of free-standing, single-layer graphene (SLG) [5], [6], [7], [8], [9], [10], [11], [12], and multi-layer graphene (MLG) [8], [9], [10], [11], [12].

Theoretical modeling of the EELS data of SLG and MLG is an active field of research [13], [14], [15], [16], [17], [18], [19]. In our previous publication [20], we treated the MLG as layered electron gas with in-plane polarizability modeled by a two-dimensional (2D), two-fluid hydrodynamic (HD) model [21] for the inter-band transitions of π and σ electrons of SLG, yielding good agreement with the experimental EEL spectra [11] for N < 10 graphene layers in STEM. We have also used the same version of the HD model for graphene's π and σ electrons in conjunction with an empirical Drude-Lorentz model for metal substrate to reproduce the momentum-resolved experimental EELS data for low-energy electron reflection (∼10 eV) from monolayer graphene supported by Pt(111), Ru(0001), and Ni(111) substrates [19], as well as for high-quality graphene grown on peeled-off epitaxial Cu(111) foils [22].

It should be mentioned that, while the agreement of the HD model with the experimental EEL spectra for MLG from Ref. [11] covered the regions around the principal π and π + σ peaks, there was no experimental data for energy losses below  ≈ 3 eV, which is a consequence of the subtraction of the zero-loss peak (ZLP). At the same time, the HD model does not incorporate the Dirac physics of low-energy excitations in graphene [20]. However, in the meantime, several EELS experiments were performed with high-energy electron beams in STEM, showing intriguing increase in spectral intensity as the energy loss decreases below  ≈ 3 eV, even after the ZLP subtraction [5], [7], [12]. In that respect, we pose a question whether the low-energy, inter-band excitations of π electrons in intrinsic graphene play any detectable role in the low-loss range of EELS, and if so, whether the new generation of monochromators can open up possibility to explore the Dirac physics of graphene in STEM. In order to address this question, we attempt to reproduce the STEM-EELS data from those experiments [5], [7], [12] by formulating an extended HD (eHD) model, which includes a Dirac correction. This correction treats the low-energy contribution of graphene's π electron inter-band transitions in a manner consistent with the Dirac-cone approximation for graphene's π electron bands near the K point in the Brillouin zone (BZ).

On the other hand, Despoja et al. have used an ab initio method to calculate the energy-loss rate of a point blinking charge in the vicinity of a graphene monolayer [23] and have obtained values for the energies of π and π + σ peaks in the EELS spectra, which are in good agreement with the reported experimental values [11]. In addition, they have calculated the so-called loss function, Im[1/ɛ(q,ω)], where ɛ(q, ω) is the dielectric function of SLG obtained by ab initio methods [24], and have obtained a very good agreement with the experimental STEM-EELS data for SLG [11]. Moreover, those authors were able to implement their ab initio method in the optical, or the long wavelength limit (q → 0), and hence compute a universal, frequency dependent 2D conductivity of SLG, σ(ω), in a broad range of energies of interest for EELS in STEM [25]. At the same time, it was observed in Ref. [5] that using a response function of graphene in the optical limit represents an excellent approximation for analytical modeling of the EELS data with the electron beam under normal incidence and for small collection angles. Taking advantage of that situation, we propose here an analytical expression for the optical conductivity σ(ω) of SLG within the eHD model, containing several free parameters which are fixed via direct comparison with the optical conductivity obtained by the ab initio method.

Moreover, taking further advantage of working in the limit of optical response of SLG, we derive an analytical expression for the probability density for losing energy ω, P(ω), of fast electrons traversing graphene under normal incidence, which takes frequency dependent conductivity σ(ω) as input. The resulting formula may be readily applied to model the EELS of any isotropic 2D material, which can be described by a scalar conductivity given in local form. Hence, we use both the eHD and ab initio results for σ(ω) of SLG to obtain probability densities P(ω) which are then directly compared with the experimental EELS data from three independent experiments.

Finally, using the eHD model with and without the Dirac correction, we explore the possible role of Dirac physics in the experimental STEM-EELS setup and its effects on the ZLP subtraction from those spectra.

Section snippets

Theoretical methods

In a typical (S)TEM-EELS experiment operating at the voltage on the order of several tens of kV (for example, 40 kV in Ref. [7], 60 kV in Ref. [5], and 100 kV in Ref. [11]) the momentum transfer of the incident electron is close to zero, so we shall use a straight-line trajectory while neglecting relativistic effects [20], [26]. We use a Cartesian coordinate system with r={R,z} and assume that SLG occupies the plane z=0, where R={x,y} is the in-plane position and z the distance from it.

Results and discussion

As regards the ab initio study, the first part of the calculation consists of determining the KS ground state of SLG and the corresponding wave functions ϕnK(R,z) and energies En(K). We use the experimental value of a=4.651 a.u. [29] for the graphene unit cell parameter in parallel direction. For the unit cell in the perpendicular (or z) direction (separation between periodically repeated graphene layers) we take L=5a=23.255 a.u. For calculating KS wave functions and energies we use a

Conclusions

We have presented an analytical modeling of the EELS data for free-standing graphene using an ab initio method and the 2D, two-fluid eHD model. An analytical expression for the probability density for energy loss of fast electrons traversing graphene under normal incidence is derived in the long wavelength limit in terms of an arbitrary form of the conductivity of graphene given in the local, i.e., frequency-dependent form. We found a very good agreement between the results obtained by ab initio

Acknowledgements

This work is supported by the Ministry of Education, Science and Technological Development of the Republic of Serbia (Project No. 45005). T.Dj. also acknowledges support from the Ministry of Education, Science and Technological Development of the Republic of Serbia (Project No. 171023). V.D. acknowledges support from the QuantiXLie Center of Excellence. T.Dj., I.R., V.D., and D.B. wish to acknowledge the support of the COST Action MP1306 “Modern Tools for Spectroscopy on Advanced Materials: a

References (36)

  • P. Wachsmuth et al.

    High-energy collective electronic excitations in free-standing single-layer graphene

    Phys. Rev. B

    (2013)
  • P. Wachsmuth et al.

    Plasmon bands in multilayer graphene

    Phys. Rev. B

    (2014)
  • J. Lu et al.

    Plasmon dispersion on epitaxial graphene studied using high-resolution electron energy-loss spectroscopy

    Phys. Rev. B

    (2009)
  • U. Bangert et al.

    STEM plasmon spectroscopy of free standing graphene

    Phys. Status Solidi A

    (2008)
  • T. Eberlein et al.

    Plasmon spectroscopy of free-standing graphene films

    Phys. Rev. B

    (2008)
  • J.C. Idrobo et al.

    A short story of imaging and spectroscopy of two-dimensional materials by scanning transmission electron microscopy

    Ultramicroscopy

    (2017)
  • M.B. Roman et al.

    Low energy intraband plasmons and electron energy loss spectra of single and multilayered graphene

    Plasmonics

    (2017)
  • D.J. Mowbray

    Theoretical electron energy loss spectroscopy of isolated graphene

    Phys. Status Solidi B

    (2014)
  • Cited by (13)

    • Conductivity models for electron energy loss spectroscopy of graphene in a scanning transmission electron microscope with high energy resolution

      2020, Ultramicroscopy
      Citation Excerpt :

      Inspired by these developments, we also explore here the effects of reducing the maximum collection angle of a monochromated electron beam and propose a computational procedure that can potentially be used to deduce the real and imaginary parts of the optical conductivity of doped graphene from experimental EELS data for a broad frequency range. We note that using the optical conductivity σ(ω) in Eqs. (2)–(4) makes it possible to complete the q integrals in Eq. (1) in an analytical form for arbitrary qc, so that the resulting expressions for Pohm(ω) and Prad(ω) simply take σ(ω) as an input function [6,34] This may be taken available from empirical analytical models, tabulated data of numerical computations, or even digitized experimental data. Finally, a few comments are in order regarding possible nonlocal effects, i.e., the q dependence of a model conductivity σ(q, ω) of graphene.

    • Backscattered electron spectra from graphite

      2018, Physics Letters, Section A: General, Atomic and Solid State Physics
      Citation Excerpt :

      The optical properties of graphite have been measured by energy-loss spectroscopy (EELS) techniques [16–18]. Also, in some of our previous study of graphene simple hydrodynamic model were used to explain HREELS experiments [19–23]. In this work the backscattered electron spectra from graphite sample are investigated both experimentally and theoretically at impact energies between 500 and 5000 eV.

    View all citing articles on Scopus
    View full text