On the X-Ray Edge Problem in Graphene

Show more

1. Introduction

The so-called X-ray edge problem is related to the absorption of high-energy electromagnetic radiation. In the process of x-ray absorption, a core electron absorbs an x-ray photon and is injected from the material. In an insulator, the energy of the ejected electron is unequivocally determined by that of the absorbed photon.

In metal, however, the presence of the conduction electrons makes this process much more interesting. As a result of the core electron ejection, there is a localized change of screened Coulomb potential seen by the conduction electrons. This sudden change of potential leads to transient perturbation of electron gas that is to creation of electron-hole pairs. This process changes the spectral function of a core hole with respect to the case of insulator. It is this spectral function which determines the energy dependence of the cross-section for ejecting photoelectrons from the core state [1] - [6] .

Starting from the discovery of graphene, large attention attracts the systems which occupy a place intermediate between an insulator and a metal. The gap between the conduction and the valence band in graphene is zero, which allows creation of electron-hole pairs with the energies down to zero, like in a metal. On the other hand, the density of electron states at the Fermi level in the pristine graphene is zero, like in an insulator. So the present paper studies core-hole spectral function in such systems, noticing however, that the results can be applied to any Dirac semimetal.

The X-ray edge problem in graphene [7] was addressed previously [8] [9] . (See also closely connected works, which studied the theory of dynamical spin response for fullerene molecules [10] , bilayer graphen [11] , topological insulator [12] , Kitaev honeycomb model [13] and intercalated graphite [14] .) Thus in Ref. [8] , both the spectral density of the core hole and the X-ray response function were calculated using the method by Nozieres and de Dominicis [3] . The calculations were done for the local interaction (see below).

In our paper, we calculate the spectral density of the core hole using the method developed by Mahan [1] . For the local interaction, we reproduce the result obtained previously [7] . (We only consider the case of pristine graphene.) For the unscreened Coulomb interaction, we obtain the result which is equivalent to the classical result for the metal with the constant density of states of itinerant electrons but with the local interaction [1] . Thus lower density of the electron- hole excitations in graphene in comparison to metal is compensated by singularity of the unscreened Coulomb interaction. As far as the calculation, the X-ray response function is concerned; we only present the general approach to calculation of a one-electron transient Green’s function L, leaving the exposition of the results until the next publication.

2. Spectral Density of the Core Hole

2.1. Spectral Density

The model Hamiltonian is [1]

$H={\zeta}_{h}{d}^{\u2020}d+{\displaystyle \underset{k}{\sum}{\zeta}_{c}}\left(k\right){c}_{k}^{\u2020}{c}_{k}+\frac{1}{V}{\displaystyle \underset{k{k}^{\prime}}{\sum}V}\left(k{k}^{\prime}\right){c}_{k}^{\u2020}{c}_{{k}^{\prime}}{d}^{\u2020}d,$ (1)

where ${d}^{\u2020},d$ are the hole creation and annihilation operators, ${c}^{\u2020},c$ are the conduction (valence) creation and annihilation operators.

The hole Green’s function is determined as

${G}_{h}\left(t\right)=-i\langle Td\left(t\right){d}^{\u2020}\left(0\right)\rangle =-i\Theta \left(t\right)\langle {\text{e}}^{iHt}d{\text{e}}^{-iHt}{d}^{\u2020}\rangle \equiv -i\Theta \left(t\right){\text{e}}^{-i{\zeta}_{h}t}\rho \left(t\right).$ (2)

The spectral density is

${A}_{h}\left(\Omega \right)={\displaystyle {\int}_{-\infty}^{\infty}\frac{\text{d}t}{2\text{\pi}}{\text{e}}^{i\left(\omega -{\zeta}_{h}\right)t}\rho \left(t\right)}\text{\hspace{0.05em}}\text{\hspace{0.05em}}.$ (3)

In the absence of interaction between the core hole and the itinerant electrons

${A}_{h}\left(\Omega \right)=\delta \left(\Omega -{\zeta}_{h}\right).$ (4)

One of the ways to find the Green’s function in the presence of interaction is the linked cluster expansion

$\rho \left(t\right)=\langle T\mathrm{exp}\left[-i{\displaystyle {\int}_{0}^{t}\text{d}{t}_{1}V\left({t}_{1}\right)}\right]\rangle =\mathrm{exp}\left[{\displaystyle \underset{l=1}{\overset{\infty}{\sum}}{F}_{l}}\left(t\right)\right],$ (5)

where

${F}_{l}\left(t\right)=\frac{{\left(-1\right)}^{l}}{l}{\displaystyle {\int}_{0}^{t}\text{d}{t}_{1}}\cdots {\displaystyle {\int}_{0}^{t}\text{d}{t}_{l}}{\langle TV\left({t}_{1}\right)\cdots V\left({t}_{l}\right)\rangle}_{\text{connected}}$ (6)

The first term in the series $\underset{l=1}{\overset{\infty}{\sum}}{F}_{l}}\left(t\right)$ just renormalizes the core level energy. Asuming the electron-hole interaction is weak, we can restrict ourselves by the second term [1]

${F}_{2}\left(t\right)=\frac{1}{{V}^{2}}{\displaystyle \underset{{k}_{1}{k}_{2}}{\sum}{\left|V\left({k}_{1}{k}_{2}\right)\right|}^{2}}\left[\frac{it{n}_{1}}{{\zeta}_{1}-{\zeta}_{2}}-\frac{{n}_{1}\left(1-{n}_{2}\right)}{{\left({\zeta}_{1}-{\zeta}_{2}\right)}^{2}}\left(1-{\text{e}}^{it\left({\zeta}_{1}-{\zeta}_{2}\right)}\right)\right].$ (7)

Using assumption $V\left({k}_{1}{k}_{2}\right)=V\left({k}_{1}-{k}_{2}\right)$ and again ignoring the shift of the hole energy due to the electron-hole interaction, described by the first term in brackets in Equation (7), one obtains

${F}_{2}\left(t\right)=-{\displaystyle {\int}_{0}^{\infty}\frac{\text{d}u}{{u}^{2}}{R}_{e}\left(u\right)}\left(1-{\text{e}}^{-iut}\right),$ (8)

where

${R}_{e}\left(u\right)=-\frac{2}{\text{\pi}V}{\displaystyle \underset{q}{\sum}{\left|V\left(q\right)\right|}^{2}\text{Im}\left[{P}^{\left(1\right)}\left(q,u\right)\right]}\text{\hspace{0.05em}}\text{\hspace{0.05em}},$ (9)

and

$-\text{Im}\left[{P}^{\left(1\right)}\left(q,u\right)\right]=\frac{2\text{\pi}}{V}{\displaystyle \underset{k}{\sum}{n}_{F}}\left({\zeta}_{k}\right)\left[1-{n}_{F}\left({\zeta}_{k+q}\right)\right]\delta \left(u+{\zeta}_{k}-{\zeta}_{k+q}\right).$ (10)

is the imaginary part of the polarization operator (in a single-bubble approximation). Here and further on we take into account spin and valley degeneracy. For the monolayer graphene [3] [4] [5]

${P}^{\left(1\right)}\left(q,u\right)=\frac{{q}^{2}}{4\sqrt{{v}_{F}^{2}{q}^{2}-{u}^{2}}}.$ (11)

2.2. Local Interaction

Let the interaction of graphene electrons with the core hole is local $V\left(q\right)=V=$ $\text{const}$ . We obtain

${R}_{e}\left(u\right)=\frac{{V}^{2}}{4{\text{\pi}}^{2}}{\displaystyle {\int}_{0}^{u/{v}_{F}}\frac{{q}^{3}\text{d}q}{\sqrt{{u}^{2}-{v}_{F}^{2}{q}^{2}}}}=\frac{{V}^{2}{u}^{3}}{6{\text{\pi}}^{2}{v}_{F}^{4}}$ (12)

Hence

${F}_{2}=-\frac{{V}^{2}}{6{\text{\pi}}^{2}{v}_{F}^{4}}{\displaystyle {\int}_{0}^{{\xi}_{0}}\left(1-{\text{e}}^{-iut}\right)u\text{d}u}\text{\hspace{0.05em}}\text{\hspace{0.05em}}.$ (13)

The upper limit of u ( ${\xi}_{0}$ ) can be taken as the band width. Actually, Equations ((9) and (10)) can be joined together to give a simple formula [15]

$\begin{array}{c}{R}_{e}\left(u\right)=4{V}^{2}{\int}_{\zeta (k>0}\frac{{\text{d}}^{2}k}{{\left(2\text{\pi}\right)}^{2}}{\int}_{\zeta ({k}^{\prime}<0}\frac{{\text{d}}^{2}{k}^{\prime}}{{\left(2\text{\pi}\right)}^{2}}\delta \left(u-\zeta \left(k\right)+\zeta \left({k}^{\prime}\right)\right)\\ =\frac{{V}^{2}}{4}{\displaystyle {\int}_{0}^{\infty}\rho}\left(\u03f5\right)\text{d}\u03f5{\displaystyle {\int}_{-\infty}^{0}\rho}\left({\u03f5}^{\prime}\right)\text{d}{\u03f5}^{\prime}\delta \left(u-\u03f5+{\u03f5}^{\prime}\right)\end{array}$ (14)

(where the second equation is correct for a pristine graphene). For monolayer graphene with the dispersion law $\u03f5\left(k\right)=\pm {v}_{F}k$ the density is

$\rho \left(\u03f5\right)=\frac{2\left|\u03f5\right|}{\text{\pi}{v}_{F}^{2}}.$ (15)

Thus we recover the final result of Equation (12).

The influence of the electron-hole interaction on the core-hole spectral density in graphene for local interaction is much weaker than in the case of Coulomb interaction. In particular, ${F}_{2}$ defined by Equation (13) converges to the same constant for $t\to \pm \infty $ . Hence δ-function singularity of the spectral density, which we had in the absence of the electron-hole scattering survives in the pristine graphene. However, the interaction gives a background in addition to this singularity, which we’ll study in details. For the case of monolayer graphene the background up to a numerical factor can be presented as

${A}_{b.g}\left(\Omega \right)={\displaystyle {\int}_{-\infty}^{\infty}\frac{\text{d}t}{2\text{\pi}}{\text{e}}^{\left(\Omega -{\zeta}_{h}\right)t}}\left({\text{e}}^{\frac{{V}^{2}}{6{v}_{F}^{4}}{\displaystyle {\int}_{0}^{{\xi}_{0}}{\text{e}}^{-iut}}u\text{d}u}-1\right).$ (16)

Calculating integral with respect to $\text{d}u$ we obtain

${A}_{b.g}\left(\stackrel{\u02dc}{\Omega}\right)={\displaystyle {\int}_{-\infty}^{\infty}\frac{\text{d}t}{2\text{\pi}}{\text{e}}^{i\stackrel{\u02dc}{\Omega}t}}\left[{\text{e}}^{g\left[i{\text{e}}^{-it}/t+\left({\text{e}}^{-it}-1\right)/{t}^{2}\right]}-1\right],$ (17)

where

$g=\frac{{V}^{2}{\xi}_{0}^{2}}{6{v}_{F}^{4}}.$ (18)

Recalling that we assumed from the beginning that g is small, we can expand the exponent in Equation (17). In lowest order approximation we obtain

${A}_{b.g}\left(\stackrel{\u02dc}{\Omega}\right)=g{\displaystyle {\int}_{-\infty}^{\infty}\frac{\text{d}t}{2\text{\pi}}{\text{e}}^{i\stackrel{\u02dc}{\Omega}t}}\left[\frac{i{\text{e}}^{-it}}{t}+\frac{{\text{e}}^{-it}-1}{{t}^{2}}\right]=g\stackrel{\u02dc}{\Omega}\theta \left(\stackrel{\u02dc}{\Omega}\right)\theta \left(1-\stackrel{\u02dc}{\Omega}\right).$ (19)

The result of Equation (19) coincides with that obtained previously [7] .

2.3. Unscreened Coulomb Interaction

Now consider the case of unscreened Coulomb interaction

$V\left(q\right)=\frac{2\text{\pi}{e}^{2}}{q}.$ (20)

Substituting imaginary part of Equation (11) and Equation (20) into Equation (9) we obtain

${R}_{e}\left(u\right)={e}^{4}{\displaystyle {\int}_{0}^{u/{v}_{F}}\frac{q\text{d}q}{\sqrt{{u}^{2}-{v}_{F}^{2}{q}^{2}}}}=\frac{{e}^{4}u}{{v}_{F}^{2}}.$ (21)

Hence

${F}_{2}=-g{\displaystyle {\int}_{0}^{{\xi}_{0}}\left(1-{\text{e}}^{-iut}\right)}\frac{\text{d}u}{u},$ (22)

where $g={e}^{4}/{v}_{F}^{2}$ . The upper limit of u is taken as the band width ${\xi}_{0}$ . It is interesting that for a monolayer graphene with Coulomb interaction that we obtain the result coinciding with that for a metal with the local interaction $V\left(q\right)=V=\text{const}$ [3] , the only difference is that for a metal $g=2{\left({N}_{F}V\right)}^{2}$ . Like in the case of metal, ${F}_{2}$ defined by Equation (22) logarithmically diverges for $t\to \pm \infty $ . This leads to disappearance of the δ-function singularity in the spectral density.

${A}_{h}\left(\stackrel{\u02dc}{\Omega}\right)={\displaystyle {\int}_{-\infty}^{\infty}\frac{\text{d}t}{2\text{\pi}}{\text{e}}^{i\stackrel{\u02dc}{\Omega}t-g{\displaystyle {\int}_{0}^{1}\left(1-{\text{e}}^{-iut}\right)}\frac{\text{d}u}{u}}},$ (23)

which is equivalent to the classical result for the metal with the constant density of states of itinerant electrons but with the local interaction [1]

${A}_{h}\left(\stackrel{\u02dc}{\Omega}\right)=2\theta \left(\stackrel{\u02dc}{\Omega}\right)\mathrm{sin}\left(\text{\pi}g\right)\Gamma \left(1-g\right)\frac{{\text{e}}^{-\stackrel{\u02dc}{\Omega}}}{{\stackrel{\u02dc}{\Omega}}^{1-g}},$ (24)

where $\stackrel{\u02dc}{\Omega}=\left(\Omega -{\zeta}_{h}\right)/{\xi}_{0}$ . Thus lower density of the electron-hole excitations in graphene in comparison to metal is compensated by singularity of the unscreened Coulomb interaction.

3. Absorption of X-Rays by Core Level to Conduction Band Transition

The net transition rate in the photoemission experiment in the time representation is expressed as the product of two factors: a one-electron transient Green’s function, and the deep-level Green’s function ${G}_{h}$ [3] . In this Section we formulate analog of the Nozieres-de Dominicis equation for the former for the case when the density of states of itinerant electrons goes linearly to zero at the Fermi level.

A one-body transient Green’s function ${L}_{k{k}^{\prime}}$ is given by

${L}_{k{k}^{\prime}}\left(t-{t}^{\prime}\right)=-{\phi}_{k{k}^{\prime}}\left({t}^{\prime},t;t,{t}^{\prime}\right),$ (25)

and one electron Green’s function (calculated in the presence of transient potential) $\phi $ is defined by equation

${\phi}_{k{k}^{\prime}}\left(\tau ,{\tau}^{\prime};t,{t}^{\prime}\right)={G}_{k{k}^{\prime}}\left(\tau -{\tau}^{\prime}\right)-i{\displaystyle {\int}_{t}^{{t}^{\prime}}\text{d}}{\tau}^{\u2033}{\displaystyle \underset{q{q}^{\prime}}{\sum}{G}_{k{q}^{\prime}}}\left(\tau -{\tau}^{\u2033}\right){V}_{q{q}^{\prime}}{\phi}_{{q}^{\prime}{k}^{\prime}}\left({\tau}^{\u2033},{\tau}^{\prime};t,{t}^{\prime}\right),$ (26)

where

${G}_{k{k}^{\prime}}\left(t\right)=\{\begin{array}{l}{\delta}_{k{k}^{\prime}}{\text{e}}^{-i{\u03f5}_{k}t}\theta \left(t\right)\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{if}\text{\hspace{0.17em}}{\u03f5}_{k}>0\\ -{\delta}_{k{k}^{\prime}}{\text{e}}^{-i{\u03f5}_{k}t}\theta \left(-t\right)\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{if}\text{\hspace{0.17em}}{\u03f5}_{k}<0\end{array}$ (27)

Consider the case of local interaction and assume that electron-hole can be taken as [3]

${V}_{q{q}^{\prime}}=-Vu\left(q\right)u\left({q}^{\prime}\right),$ (28)

where u is a cut-off function. In this case the moments are decoupled, and we can write down closed equation for the quantity

$\phi ={\displaystyle \underset{k{k}^{\prime}}{\sum}{\phi}_{k{k}^{\prime}}}.$ (29)

Retaining only the “active” variable one obtains [3]

$\phi \left(\tau \right)=G\left(\tau -{\tau}^{\prime}\right)+iV{\displaystyle {\int}_{t}^{{t}^{\prime}}\text{d}}{\tau}^{\u2033}G\left(\tau -{\tau}^{\u2033}\right)\phi \left({\tau}^{\u2033}\right).$ (30)

The Green’s function is

$G\left(t\right)=\{\begin{array}{l}{\displaystyle {\int}_{0}^{\infty}\rho}\left(\u03f5\right){u}^{2}\left(\u03f5\right){\text{e}}^{-i\u03f5t}\text{d}\u03f5\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\left(t>0\right)\\ -{\displaystyle {\int}_{-\infty}^{0}\rho}\left(\u03f5\right){u}^{2}\left(\u03f5\right){\text{e}}^{-i\u03f5t}\text{d}\u03f5\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\left(t<0\right)\end{array}$ (31)

Taking the cut-off function as

$u\left(\u03f5\right)={\text{e}}^{-\left|\u03f5\right|/\xi},$ (32)

one obtains

$G\left(t\right)=\frac{2}{\text{\pi}{v}_{F}^{2}}\frac{1}{{\left(i\left|t\right|+1/\xi \right)}^{2}}.$ (33)

Equation (30) for $t=0$ and ${t}^{\prime}\to \infty $

$\phi \left(\tau \right)=G\left(\tau -{\tau}^{\prime}\right)+iV{\displaystyle {\int}_{0}^{\infty}\text{d}}{\tau}^{\u2033}G\left(\tau -{\tau}^{\u2033}\right)\phi \left({\tau}^{\u2033}\right),$ (34)

After being Fourier transformed becomes a Wiener-Hopf equation

$\phi \left(\omega \right)=G\left(\omega \right){\text{e}}^{i\omega {\tau}^{\prime}}+iVG\left(\omega \right){\Phi}_{+}\left(\omega \right).$ (35)

The functions ${\Phi}_{+}\left(\omega \right)$ and ${\Phi}_{-}\left(\omega \right)$ , analytic in the upper and lower half planes respectively, are defined as follows

${\Phi}_{\pm}\left(\omega \right)=\pm \frac{1}{2\text{\pi}i}{\displaystyle {\int}_{-\infty}^{+\infty}\frac{\phi \left({\omega}^{\prime}\right)\text{d}{\omega}^{\prime}}{{\omega}^{\prime}-\omega \mp i\delta}}.$ (36)

Using the decomposition

$\phi \left(\omega \right)={\Phi}_{+}\left(\omega \right)+{\Phi}_{-}\left(\omega \right),$ (37)

we rewrite Equation (35) as

$\left[1-iVG\left(\omega \right)\right]{\Phi}_{+}\left(\omega \right)+{\Phi}_{-}\left(\omega \right)=G\left(\omega \right){\text{e}}^{i\omega {\tau}^{\prime}}.$ (38)

Solution of Equation (38) with the Green’ function (33) will be a subject of further research.

4. Conclusions

We calculated the spectral density of the core hole for pristine graphene for two models of electron-hole interaction: unscreened Coulomb and local interaction.

The calculations are performed using the method developed by Mahan. We also formulate analog of Nozieres-de Dominices equation for calculating the net transition rate in the photoemission experiment in graphene.

Acknowledgements

E.K. cordially thanks Max-Planck-Institut fur Physik Komplexer Systeme for the hospitality extended to him during his visit there.

References

[1] Mahan, G. (1967) Exitons in Metals: Infinite Hole Mass. Physical Review, 163, 612-617.

[2] Anderson, P.W. (1967) Infrared Catastrophe in Fermi Gases with Local Scattering Potentials. Physical Review Letters, 18, 1049-1051.

http://doi.org/10.1103/Phys.Rev.Lett.18.1049

[3] Nozieres, P. and de Dominicis, C.T. (1969) Singularities in the X-Ray Absorption and Emission of Metals. III. One-Body Theory Exact Solution. Physical Review, 178, 1097-1107.

[4] Schotte, K.D. and de Schotte, U. (1969) Tomonaga’s Model and the Threshold Singularity of X-Ray Spectra of Metals. Physical Review, 182, 479-482.

[5] Doniach, S. and Sondheimer, E.H. (1998) Green’s Functions for Solid State Physicists. Imperial College Press, London.

https://doi.org/10.1142/p067

[6] Mahan, G.D. (2000) Many-Particle Physics. 3rd Edition, Kluver Academics/Plenum Publishers, New York.

https://doi.org/10.1007/978-1-4757-5714-9

[7] Castro Neto, A.H., Guinea, F., Peres, N.M.R., Novoselov, K.S. and Geim, A.K. (2009) The Electronic Properties of Graphene. Reviews of Modern Physics, 81, 109-162.

http://doi.org/10.1103/RevModPhys.81.109

[8] Yang, S.-R.E. and Lee, H.C. (2007) X-Ray Edge Problem of Graphene. Phys. Rev. B, 76, 245411-1-245411-7.

[9] Gonzalez, J., Guinea, F. and Vozmediano, M.A.H. (1993) Electrostatic Screening in Fullerene Molecules. Modern Physics Letters B, 7, 1593-1599.

http://doi.org/10.1142/S0217984993001612

[10] Roder, G., Tkachov, G. and Hentschel, M. (2011) Photoabsorption Spectra and the X-Ray Edge Problem in Graphene. Europhysics Letters, 94, 67002-p1-67002-p6.

https://doi.org/10.1209/0295-5075/94/67002

[11] Lee, H.C. (2010) X-Ray Edge Singularity of Bilayer Graphene. Phys. Rev. B, 82, 153410-1-153410-10.

[12] Lee, H.C. (2010) Theoretical Study of X-Ray Absorption of Three-Dimensional Topological Insulator Bi2Se3. Phys. Rev. B, 83, 193107-1-193107-4.

[13] Knolle, J., Kovrizhin, D.L., Chalker, J.T. and Moessner, R. (2015) Dynamics of Fractionalization in Quantum Spin Liquids. Phys. Rev. B, 92, 115127-p1-115127-p20.

[14] Shung, K.V.K. (1986) Dielectric Function and Plasmon Structure of Stage-1 Intercalated Graphite. Phys. Rev. B, 34, 979-993.

[15] Gogolin, A.O., Nersesyan, A.A. and Tsvelik, A.M. (1998) Bosonization and Strongly Correlated Systems. Cambridge University Press, Cambridge.