arXiv:1012.3484v2 [cond-mat.str-el] 26 Feb 2012 Contents Stat Current Graphene: in Interactions Electron-Electron I.Qaiatce nGraphene in Quasiparticles III. V h olm rbe n hre Impurities Charged and Problem Coulomb The IV. I hreplrzto n ierscreening linear and polarization Charge II. .Introduction I. .Pyia Observables Physical D. screening and modes Collective D. .Lweeg eairna h ia point Dirac the near behavior Low-energy A. spectrum Tight-binding A. .Fnt est em-iudregime Fermi-liquid density Finite C. function Polarization C. .Sotnosms generation mass Spontaneous B. Hamiltonian fermion Dirac B. .Ifiiesako layers of stack Infinite E. .Oeve fmi results main of Overview E. .fsmrule f-sum F. ngahn,w hwta ltoao neetn susre st is issues research interesting graphene of of plethora field the a that that and show experimentally, as br under we main fully also the the graphene, We reviewing of in to some physics. addition review In mesoscopic als and bilayers. layers as is graphene single invariance well in Lorentz as fields spe effects emerging magnetic whole edge an the and of span a that size effect effects is The We to graphene metals. leads infrared. in and ordinary role in adatoms fundamental effects and/or a many-body impurities of ana models of pro the classical The effect analyze problem. also Kondo impurity We Coulomb and the regime. and coupling problem reg strong many-body coupling the e weak in the the discuss states in we quasi-particles systems, invariant these in Lorentz in interactions interactions range electron-electron long of of problem the review We 2018) 22, October (Dated: Nation Physics, of Department 117542 and Singapore Centre Commonwealt Research 590 University, Graphene Boston Physics, of Department Neto In´e Castro Juana H. Sor Madrid, A. de Materiales de Ciencia de Instituto Guinea Nation F. Physics, of Department 117542 and Singapore Centre Research Graphene Pereira M. Urbana-Ch Vitor at Illinois of University Physics, of Department Uchoa Universit Bruno 82 Vermont, of University Physics, of Department Kotov N. Valeri .Conductivity response spin 2. and Charge 1. generation mass Excitonic mass 2. explicit Finite 1. lifetime Quasiparticle 3. analysis Strong-coupling/RPA 2. analysis Weak-coupling 1. 22 22 21 19 19 18 16 16 16 15 13 10 10 10 9 9 7 6 5 4 4 2 to set ftemn-oyproblem many-body the of aspects stood m,adsrnl orltdelectronic correlated strongly and ime, itneo neegn ia iudof liquid Dirac emerging an of xistence rpee trigfo h screening the from Starting graphene. l xiigadvibrant. and exciting ill mag,11 etGenSre,Ubn,Ilni 61801 Illinois Urbana, Street, Green West 1110 ampaign, icse ntecneto finite of context the in discussed o lc,Brigo,Vrot05405 Vermont Burlington, Place, y et ftemn-oypolmin problem many-body the of pects hwta oet naineplays invariance Lorentz that show anoe,bt hoeial and theoretically both open, main el rz3 -84 ard Spain Madrid, E-28049 3, Cruz la de s I.Itrcineet nmssoi systems mesoscopic in effects Interaction VII. oyadcnetosbtenthe between connections and logy tu,fo h lrvoe othe to ultraviolet the from ctrum, I neatosa onaisadltiedefects lattice and boundaries at Interactions VI. eydsusteeet fstrong of effects the discuss iefly vne otn ascuet 21 and 02215 Massachusetts Boston, Avenue, h .Srn orltosi graphene in correlations Strong V. lmo h antcinstability magnetic the of blem s icse naaoywith analogy in discussed lso lUiest fSnaoe cec rv 3, Drive Science 2 Singapore, of University al 3, Drive Science 2 Singapore, of University al .Knoeffect Kondo D. Experiments in Physics Supercritical D. .Mgeimi unu dots quantum in Magnetism A. states Surface A. lattice honeycomb the in gaps Mass A. Problem Coulomb the of Solution Exact A. .Mda ttsadRno ag Fields Gauge Random and States Midgap C. moments magnetic Local C. Interactions Particle Many to Single From C. .Sae tvcnisadcracks and vacancies at States B. instabilities magnetic and Charge B. Screening and Charge Induced B. .RK interaction RKKY E. .Superconductivity F. .Ectn n pnaeu asGeneration Mass Spontaneous and Excitons problem Body 2. Two Interacting 1. Mass Finite ( 3. Coupling Strong ( 2. Coupling Weak Sections 1. Cross Transport and Scattering DOS, 3. Instabilities Supercritical Spectrum and 2. Equations Wave 1. sadPerspectives and us g < g g > g c ) c ) 44 44 43 43 42 42 38 38 36 34 32 32 32 31 31 31 31 30 29 26 26 25 23 23 23 2
B. Charging effects. Coulomb blockade 45 ergy given by
2 VIII. Interactions in strong magnetic fields 46 p K0 = , (1.1) 2m∗ IX. Interactions in bilayers 48 A. Charge polarization 50 where p is the electron momentum and m∗ is a free pa- B. Quasiparticles 51 rameter of the theory that, for lack of a better name, C. Many-body instabilities 52 is called effective mass. Fermi-Dirac statistics implies that electrons carry spin 1/2 and that, in the ground X. Conclusions 53 state, all states with energy below the so-called Fermi energy, EF , are occupied, and all the states above it are Acknowledgments 55 empty. With these two basic assumptions and simple considerations about electron scattering by defects, the References 55 Drude-Sommerfeld model was capable of describing ex- perimental data of several generations of scientists. The understanding of why these two assumptions are I. INTRODUCTION valid for a strongly interacting problem, such as electrons in a metal, had to wait for the development of two major One of the most important problems in theoretical concepts: (i) the band structure theory that explains that physics is the understanding of the properties of quantum the interaction of the electrons with a periodic lattice of systems with an infinitely large number of interacting de- ions produces states that, as the plane waves described by grees of freedom, the so-called many-body problem. In- (1.1), are extended over the entire lattice (Bloch, 1928); teractions are present in almost all areas of physics: soft and (ii) the theory of screening, that is, that metals and hard condensed matter, field theory, atomic physics, are dynamically polarizable materials and that electrons quantum chemistry, nuclear physics, astrophysics, and so act collectively to screen electric fields in their interior on. Interactions between particles are responsible for a (Lindhard, 1954). Hence, long range Coulomb interac- plethora of effects and many-body states, from the band tions become effectively short ranged and weak enough structure of crystals to superconductivity in metals, from to give substance to Drude’s assumptions. In this case the quark-gluon plasma in heavy ion collisions to asymp- the effective mass m∗ reflects the change in the inertia of totic freedom in quantum chromodynamics (QCD). It is the electron as it moves around in an effective medium. the competition between the kinetic energy of the parti- Nevertheless, there are situations when these assump- cles, that is, their inertia, and interactions among them tions fail even in crystalline systems, and that is when that leads to the richness and complexity of these differ- interesting things happen, namely, the free electron pic- ent phases. For these reasons, many-body interactions ture breaks down. are very specific, and the hardest to describe theoreti- In fact, there are many instances where the Fermi liq- cally. uid ground state becomes unstable. Electrons not only One of the greatest theoretical achievements of the interact with static ions but also with their vibrations, last century, the Landau theory of the Fermi liquid the phonons. Electron-phonon interactions, in the pres- (Baym and Pethick, 1991), asserts something very sim- ence of strong screening, can lead to an effective at- ple but, at the same time, very deep: that the excita- tractive interaction between electrons producing a catas- tions of a large (indeed, infinite) collection of strongly trophic Fermi surface instability towards a superconduct- interacting particles can be described as an equally large ing ground state (Tinkham, 1996). Fermi surface insta- collection of weakly-interacting quasi-particles that carry bilities also happen in special situations in the presence the same quantum numbers as the original particles. This of Fermi surface nesting which can lead to charge and statement is far from trivial. Consider, for instance, the spin density wave ground states (Gruner, 1994). Crys- behavior of electrons in a metal. The electrons inter- tals with inner shell electrons, such as transition metals, act among themselves and with the ions in the crystal can also have many-body instabilities due to the strong via strong long-range Coulomb interactions. It is not local interactions between the electrons, leading to insu- at all clear what is the outcome of this complex inter- lating states with magnetic properties as in the case of acting problem. Without having any deep theoretical re- Mott insulators (Mott, 1949). Another important case of sources to treat this problem, except an extraordinary in- Fermi liquid breakdown is when the electron density is tuition, visionaries like Paul Drude (Drude, 1900a,b) and very low and the screening disappears. Arnold Sommerfeld (Hoddeson et al., 1987) settled the Notice that in quantum mechanics the momentum of foundations for the understanding of this complex prob- the particle relates to its wavelength, λ, by p = ~/λ lem by postulating, shamelessly, that (1) electrons prop- and hence the kinetic energy (1.1) behaves as K = ~2 2 agate freely in a non-relativistic (Galilean invariant) way /(2m∗λ ). If the average distance between electrons is (Drude’s contribution), and (2) electrons obey Fermi- ℓ we see that the average kinetic energy per electron has ~2 2/d d Dirac statistics (Sommerfeld’s contribution). Galilean to be of the order EK nd (2m∗) where nd =1/ℓ is invariance dictates that the electrons have a kinetic en- the average electron density≈ in d spatial dimensions. On 3 the other hand, the Coulomb interaction is given by: graphene’s electron density as n. As the electronic prop- erties of graphene are sensitive to environmental condi- e2 tions, they will be modified by the presence of other lay- V (r)= , (1.2) ǫ0r ers. In fact, as we are going to show, bilayer graphene has properties which are rather different than its mono- where e is the electron charge, and ǫ0 the dielectric con- layer counterpart. Furthermore, due to the same peculiar stant of the medium. Notice that the Coulomb energy per dispersion relation, the electronic density of states, ρ(E), 2 1/d 2 electron is of the order EC e n /ǫ0. Thus, the ratio vanishes at the Dirac point, ρ(E) E /v , and hence ≈ d ∝ | | F of Coulomb to kinetic energy is given by rs = EC /EK graphene is a hybrid between an insulator and a metal: 1/d 2 2 d ∝ (n0/nd) , where n0 = [m∗e /(~ ǫ0)] depends only on neutral graphene is not a metal because it has vanishing material properties. Therefore, at high electron densities, density of states at the Fermi energy, and it is not an n n , the kinetic energy dominates over the Coulomb insulator because it does not have a gap in the spectrum. d ≫ 0 energy, which can be disregarded, and the Fermi liq- This means that pristine (or lightly doped) graphene can- uid description is safe. At low densities, n n the not screen the long range Coulomb interaction in the d ≪ 0 Coulomb energy is dominant and new electronic phases, usual (metallic) way, although it is possible to produce such as ferromagnetism and Wigner crystallization, can electronic excitations at vanishingly small energy. This become stable (Ceperley, 1978). Therefore, the rela- state of affairs makes of graphene a unique system from tive strength of the kinetic to Coulomb interactions in the point of view of electron-electron interactions. The Galilean invariant systems is completely controlled by long-range interactions lead to non-trivial renormaliza- the electron density. Notice that in all the cases dis- tion of the Dirac quasiparticle characteristics near the cussed above the Galilean invariance was kept intact and charge neutrality point, and the resulting electronic state the driving force for the many-body instabilities was the can be called Dirac liquid, to be distinguished from the enhancement of the Coulomb relative to the kinetic en- Fermi liquid behavior at finite chemical potential (away ergy. from the Dirac point, where conventional screening takes With the advent of graphene (Novoselov et al., 2004a), place.) a two dimensional crystal of pure carbon, this picture has The unusual relation between kinetic and Coulomb en- changed and a new example of Fermi liquid breakdown ergies not only affects the electron-electron interactions, has emerged in a big way. In graphene, due to its peculiar but also the interactions of the electrons with charged lattice structure, the electrons at the Fermi energy are impurities, the so-called Coulomb impurity problem. In described in terms of an effective Lorentz invariant theory a metal described by a Galilean invariant theory of the where the kinetic energy is given by the Dirac dispersion form (1.1), screening also makes the interaction with the (Castro Neto et al., 2009a) impurity short ranged, and hence the scattering problem effectively reduces to the one of a short range impurity. K = v p , (1.3) G ± F | | In graphene, because of the lack of screening the situ- ation is rather different, and one has to face the prob- where v is the Fermi-Dirac velocity, and the signs F lem of the effect of the long range part of the poten- refer to two linearly dispersing bands. If we take± (1.3) at tial. Scattering by long range interactions has a long face value and reconsider the argument given above on history in physics and it leads to the issue of logarith- the relevance of the Coulomb interactions we reach very mic phase shifts (Baym, 1969). In graphene, because different conclusions. For one, the form of the Coulomb of its emergent Lorentz invariance, this issue is magni- interaction remains the same as in (1.2), since v is a ma- F fied. Since Coulomb interactions between electrons and terial’s property and hence much smaller than the speed electron scattering by Coulomb impurities are closely re- of light, c. This means that the photons which mediate lated issues, one expects that many of the anomalies of the Coulomb interaction are still much faster than the one problem are also reflected in the other. electrons and, thus, the electron-electron interaction can be considered as instantaneous. Therefore, the Coulomb Another interesting consequence of the scaling of the interaction (1.2) actually breaks the Lorentz invariance kinetic energy with momentum is related to the issue of of (1.3). Secondly, because of the linear scaling of the ki- electron confinement. If electrons are confined to a re- netic energy with momentum, we see that the average ki- gion of size L the energy of the states is quantized, no 1/2 netic energy per electron has to scale like EG ~vF n matter whether the electrons obey Galilean or Lorentz and consequently the ratio of Coulomb to kinetic≈ energy invariance. However, the quantization of energy is rather is given by different in these two cases. In a Galilean invariant sys- tem, like the one described by (1.1) the energy levels 2 2 EC e are spaced as ∆E0 1/L while in graphene Lorentz α = = , (1.4) invariance, (1.3), implies∝ ∆E 1/L. Hence, the size EG ǫ0~vF G dependence of the energy levels in∝ sufficiently small sam- and is independent of the electronic density n, depend- ples of graphene is rather different than one would find ing only on material properties and environmental con- in normal metals. Moreover, since the Coulomb energy ditions, such as ǫ0. Here, and from now on, we refer to scales like 1/L we expect Coulomb effects to be stronger 4 in nanoscopic and mesoscopic graphene samples. a) b)
4 −
Furthermore, the fact that graphene is a two dimen- 0 sional (2D) system has strong consequences for elec- −4 − tronic motion in the presence of perpendicular magnetic Energy (eV) fields. Since a perpendicular magnetic field B leads to Γ M K Γ a quantization of the energy in terms of Landau lev- els, and the electrons cannot propagate along the di- rection of the field, its effect is singular, in the sense c) d) E that the problem has a massive degeneracy. So, strong M magnetic fields can completely quench the kinetic energy K K of the electrons that become dispersionless. The elec- tronic orbits are localized in a region of the size of the Γ magnetic length: ℓB = ~c/(eB). For a Galilean in- variant system, such as the one described by (1.1), for K ’ p p ~/ℓB the kinetic energy per electron is of order K ≈ ~ω B where ω = ~/(m ℓ2 ) is the cyclotron C C ∗ B FIG. 1 (Color online) a) Honeycomb lattice with the two frequency.≈ ∝ On the other hand, for graphene, using (1.3), sublattices in graphene. The red arrows are nearest neighbor one has E ~ω √B where ω = √2v /ℓ , which ∗ G ≈ G ∝ G F B vectors. b) Tight-binding spectrum for the π π bands. The is a consequence of the Lorentz invariance. Notice that horizontal line intersecting the K point corresponds− to the in both cases the Coulomb energy per electron scales like Fermi level at half-filling. c) Brillouin zone centered around 2 EC e /(ǫ0ℓB) √B. Hence, in a Galilean invariant the Γ point. d) Dirac cone resulting from the linearization of system∝ the Coulomb∝ energy is smaller than the kinetic the tight-binding spectrum around the K points (blue circles). energy at high fields while for Lorentz invariant systems they are always comparable. Thus, one expects Coulomb interactions to be hugely enhanced in the presence of II. CHARGE POLARIZATION AND LINEAR SCREENING these magnetic fields. In the 2D electron gas (2DEG) this unusual state of affairs is what leads to the fractional A. Tight-binding spectrum quantum Hall effect (FQHE) (Laughlin, 1983). In isolated form, carbon has six electrons in the orbital configuration 1s22s22p2. When arranged in the honey- comb crystal shown in Fig.1(a), two electrons remain in Given all these unusual circumstances, many questions the core 1s orbital, while the other orbitals hybridize, come to mind: How does screening of the long range 2 2 forming three sp bonds and one pz orbital. The sp Coulomb interaction work in graphene? Can graphene orbitals form the σ band, which contains three localized be described in terms a Lorentz invariant theory of quasi- electrons. The bonding configuration among the pz or- particles? Is the Coulomb impurity problem in graphene bitals of different lattice sites generates a valence band, or the same as in a normal metal? In what circumstances π-band, containing one electron, whereas the antibonding is graphene unstable towards many-body ground states? configuration generates the conduction band (π∗), which Are there quantum phase transitions (Sachdev, 1999) in is empty. the phase diagram of graphene? Do magnetic moments From a kinetic energy point of view, the electronic sin- form in graphene in the same way as they do in normal gle particle dispersion in graphene is essentially defined metals? What is the ground state of graphene in high by the hopping of the electrons between nearest neighbor magnetic fields? carbon sites in the honeycomb lattice. Unlike square or triangular lattices, the honeycomb lattice is spanned by two different sets of Bravais lattice generators, forming The objective of this review is not to cover the basic a two component basis with one set for each triangular aspects of graphene physics, since this was already cov- sublattice. Defining a label for electrons sitting in each ered in a recent review (Castro Neto et al., 2009a), but of the two sublattices, say A and B, the free hopping to try to address some of these questions while keeping Hamiltonian of graphene is others open. The field of many-body physics will always be an open field because a seemingly simple question al- 0 = t a† (Ri)bσ(Rj ) + h.c. µ nˆσ(Ri), H − σ − ways leads to another question even more profound and σ, ij σ,i Xh i X harder to answer in a definitive way. In many ways, what (2.1) we have done here is to only scratch the surface of this where aσ(Ri), bσ(Ri) are fermionic operators for sublat- rich and important field, and leave open a large number tices A and B respectively,n ˆσ(Ri) is the number op- of interesting and unexplored problems. erator, σ = , labels the spin and ij means summa- ↑ ↓ h i 5 tion over nearest neighbors. The two energy scales in around the valleys centered at K gives rise to an the Hamiltonian are t 2.8 eV, which is the hopping effective low energy description± of the electrons that energy between nearest≈ carbons, and µ, the chemical mimics the spectrum of massless Dirac particles. In potential away from half-filling [see Fig.1(b)]. In a ho- this effective theory, the elementary excitations around mogeneous system, deviations from half-filling (µ = 0) the Fermi surface are described by a Dirac Hamiltonian are routinely induced either by charge transfer from a (Semenoff, 1984), substrate (Giovannetti et al., 2008), by application of a back gate voltage (Novoselov et al., 2005, 2004a,b), or else by chemical doping (Calandra and Mauri, = Ψ† [vk γ µτ σ ]Ψk , (2.6) H0 kσ · − 0 ⊗ 0 σ 2007; Gr¨uneis et al., 2009; McChesney et al., 2010; σk Uchoa et al., 2008b). X In momentum space the free Hamiltonian of graphene where is
Ψkσ = (aK+k,σ,bK+k,σ,b K+k,σ,a K+k,σ) (2.7) µ tφp − − 0 = Ψp† ,σ − − Ψp,σ , (2.2) H tφp∗ µ p,σ − − X is a four component spinor for sublattice and valley de- grees of freedom. In this representation, γi = τ3 σi where Ψp,σ = (ap,σ,bp,σ) is a two component spinor and , where τ and σ are the usual Pauli matrices, which⊗ 3 operate in the valley and sublattice spaces respectively ip ai (i =1, 2, 3 correspond to x, y and z directions, and τ =1 φp = e · (2.3) 0 i=1 and σ0 = 1 are identity matrices). The form of the spec- X trum mimics the relativistic cone for massless fermions is a tight-binding function summed over the nearest (Wallace, 1947), neighbor vectors E (k)= v k µ (2.8) a √3 a √3 ± ± | |− a1 = ax,ˆ a2 = xˆ + a y,ˆ a3 = xˆ a yˆ , − 2 2 − 2 − 2 where the Fermi velocity v = (3/2)ta 6eVA˚ is nearly (2.4) 300 times smaller than the speed of light,≈ i.e. v 1 where a 1.42A˚ is the carbon-carbon spacing. The di- 106m/s. From now on we set ~ = k = 1 everywhere,≈ × agonalization≈ of Hamiltonian (2.2) yields the spectrum B except where it is needed. For simplicity of notation, we of the two π-bands of graphene in tight-binding approx- call the Fermi velocity v (i.e. v v) throughout this imation (Wallace, 1947), F review. ≡
E (p)= t φp µ . (2.5) The Hamiltonian (2.6) is invariant under a pseudo- ± ± | |− time reversal symmetry operation, = i(τ0 σ2) , 1 S ⊗ C H − = H, ( is the complex conjugation operator), The +( ) sign in the spectrum corresponds to the con- whichS S is equivalentC to a time reversal operation for each duction− (valence) band. valley separately. It is also invariant under a true time The hexagonal Brillouin zone (BZ) of graphene shown reversal symmetry (TRS) operation, which involves an in Fig.1(c) has three high symmetry points: the Γ point, additional exchange between the valleys, = (τ σ ) . T 1 ⊗ 1 C located at the center of the BZ, the M point, which In the absence of back scattering connecting the two indicates the position of the Van Hove singularities of valleys, the Hamiltonian can be decomposed in two in- the π-π∗ bands, where the density of states (DOS) is dependent valley species of Dirac fermions with opposite logarithmically divergent, and the K points, where the chiralities: π-bands touch, and the DOS vanishes linearly. An ex- tensive description of the band structure of graphene = Ψ† [vk σ µ]Ψ k , (2.9) H0,+ +,kσ · − +, σ and its electronic properties is reviewed in detail by σ,k Castro Neto et al., 2009a. X 0, = Ψ† ,kσ [ vk σ∗ µ]Ψ ,kσ, (2.10) H − − − · − − k Xσ, B. Dirac fermion Hamiltonian where Ψ ,kσ = (a K+k,σ,b K+k,σ) are two component The topology of the Fermi surface in undoped spinors.± In this review,± unless± otherwise specified, we graphene is defined by the six K points where the con- will arbitrarily choose one of the two cones and assume duction and valence bands touch, E (K) = φK = 0. an additional valley degeneracy in the Hamiltonian. So These special points form two sets± of nonequivalent±| | valley indexes will be generically omitted unless explicitly points, K and K′, with K = K′ and K =4π/(3√3a), mentioned. A more detailed description of the symmetry which cannot be connected− by the generators| | of the properties of the graphene Hamiltonian can be found in reciprocal lattice. The linearization of the spectrum (Gusynin et al., 2007). 6
FIG. 2 Diagram for¡ the polarization bubble corresponding to eq. (2.12).
C. Polarization function
The Green’s function of graphene is a 2 2 matrix represented in the sublattice basis by ×
G G Gˆ(k, τ)= aa ab , Gba Gbb where Gaa = T [ak(τ)ak† (0)] and so on, with τ as the imaginary time.−h In the low energyi sector of the spectrum, close to the Dirac points, the non-interacting Green’s 1 function is Gˆ(0)(k,iω) = [iω + µ vk σ]− , or equiv- alently, in a chiral representation, − ·
1 1+ sσˆ Gˆ(0)(k,iω)= k , (2.11) 2 iω + µ sv k s= X± − | | whereσ ˆk = σ k/ k is twice the quantum mechanical helicity operator· for| a| Dirac fermion with momentum k, and s = labels the two branches with positive and negative energy± in one cone. It is clear that the positive and negative branches within the same cone have also opposite helicities. The polarization function in one loop is calculated di- FIG. 3 (Color online) Polarization bubble Π(1)(q,ω) for rectly from the bubble diagram shown in Fig. 2, graphene, within the Dirac approximation. Panels (a) and (b) show, respectively, a density plot of the real and imagi- (1) Π (q,iω)= N s,s′ (p, q) nary parts of the polarization bubble, Π(1)(q,ω), defined in p ′ F × X Xs,s eq. (2.12), and normalized to the DOS at the Fermi level, f[Es′ (p + q)] f[Es(p)] ρ(µ). Panels (c) and (d) present constant frequency cuts at − , (2.12) ω/µ = 0.5, 1.0, 1.5, 2.0, 2.5, 3.0. In panel (e) we show the Es′ (p + q) Es(p) iω − − static limit, Π(1)(q, 0), whose closed form expression is writ- 1 ten in eq. (2.16). Notice the transition from a constant value where f(E)= eE/T +1 − is the Dirac-Fermi distribu- (q < 2kF ) to the linear in q dependence at large momenta. tion, with T as temperature, N = 4 is the degeneracy for The derivative of the polarization is shown in the same panel, two spins and two valleys, and and can be seen to vary continuously. In (f) we plot the real (black/solid) and imaginary (red/dashed) parts of the uni- 1 form limit (2.17). ′ (p, q)= tr(1 + ss′σˆpσˆp q) (2.13) Fs,s 4 + are the matrix elements due to the overlap of wavefunc- the screening of charge is completely suppressed, and the tions for intraband (s = s′) and interband (s = s′) tran- − polarization function describes the susceptibility of the sitions. ’tr’ means trace over the sublattice indexes. In vacuum to particle-hole pair production, exactly as in a more explicit form, s,s′ (p, q) = [1 + ss′ cos θp,p+q] /2, the diagonal time component of the polarization tensor where θ is the angle betweenF p and p + q. The full mo- in massless Quantum Electrodynamics (QED), QED2+1 mentum, frequency, and chemical potential dependence (Appelquist et al., 1988; Gonz´alez et al., 1994; Pisarski, of (2.12) is shown in panels (a-d) of Fig. 3. 1984), In metals, screening is a many-body property directly related to the polarizability of the electrons around the 1 q2 Π(1)(q,ω)= . (2.14) Fermi surface. In graphene, because the density of −4 v2q2 ω2 states (DOS) vanishes linearly around the Dirac points, − ρ(E) E µ /v2, exactly at the neutrality point (µ = 0) Here we have performed a Wickp rotation to real frequen- ∝ | − | 7 cies, iω ω +0+. Since the Fermi surface in this case for ω < vq, which is due only to intraband transitions. is just a→ point, there is no phase space for intraband ex- The polarization function in graphene is a regular func- citations at zero temperature due to the Pauli principle. tion everywhere except at ω = vq, where it has an on- The process of creation of particle-hole pairs involves in- shell singularity delimiting the| | border of the particle-hole coherent excitations of electrons from the lower to the continuum. upper band. The continuum of particle-hole excitations The polarization was derived originally by Shung, is well defined for all virtual transitions with ω >vq. 1986a and later rederived by a number of authors (Ando, For finite µ there is a crossover in the behavior of the 2006; Barlas et al., 2007; Hwang and Das Sarma, 2007; polarization function. The DOS around the Fermi level is Wunsch et al., 2007). These results rely on the cone ap- finite and the intraband excitations dominate the infrared proximation, which ignores contributions coming from behavior of the polarization. For vq µ and ω µ , ≪ | | ≪ | | the non linear part of the spectrum. In addition, the band the leading term in the polarization function is (Shung, width is assumed to be infinite. Although the charge po- 1986a) larization for Dirac fermions in 2D is well behaved and does not require cut-off regularization in the ultraviolet, (1) 2 µ ω the physical cut-off of the band, D, generates small cor- Π (q,ω) | 2| 1 . (2.15) ≈− πv − ω2 v2q2 ! rections that vanish only in the D limit. In this − sense, the ‘exact’ expression for the→ static ∞ polarization As in a Fermi liquid, there isp a particle-hole continuum function (ω = 0) for arbitrary momentum is
2 (1) 2kF q 2kF 2kF 1 2kF π Π (q, 0) = + θ(q 2k ) 1 + sin− , (2.16) − πv − F 2πv q − q q − 2 s
where kF = µ /v is the Fermi momentum, and θ(x) interband transitions is ω > 2 µ vq for q < 2kF , as is a step function.| | The static polarization is plotted in shown schematically in Fig. 5. | |− Fig. 3(e). At q 2k the static polarization exhibits a crossover ≈ F from a two dimensional electron gas (2DEG) to Dirac D. Collective modes and screening fermion behavior. For details of the polarization func- tion in the 2DEG please refer to Fig. 4. As in the 2DEG, The Coulomb interaction among the electrons in the polarization of graphene is constant for q < 2kF . graphene gives rise to collective modes and metallic For q > 2kF , it eventually becomes linear in q for large screening when the Fermi level is shifted away from the momenta. At the crossover, the static polarization and Dirac points. In a 2D system, the bare Coulomb interac- its first derivative are continuous at q = 2kF . The dis- tion is given by continuity only appears in the second derivative. This is distinct from the 2DEG case, where the first derivative is 2πe2 V (q)= , (2.18) discontinuous. The difference will affect the spacial de- ǫ q pendence of the Friedel oscillations in the two systems. 0 In the opposite limit, for arbitrary ω and q 0, the where e is the charge of the electron and ǫ is the effec- → 0 polarization function becomes tive dielectric constant of the medium. For graphene in contact with air and a substrate with dielectric constant q2 2 µ 1 2 µ ω κ, ǫ = (1+ κ)/2. In most of the experiments, graphene Π(1)(q 0,ω)= | | + ln | |− , 0 → 2πω ω 2 2 µ + ω lies on top of some substrate like SiO2 or SiC, where di- | | (2.17) electric effects are moderate (for instance, the dielectric which is shown in Fig. 3(f). The presence of a pocket constant of SiO2 is κ 4). The background dielectric ≈ of electrons (holes) around the Dirac points opens a gap constant can be significantly enhanced in the presence of in the particle hole continuum for interband excitations substrates in contact with strong dielectric liquids such (ω > vq). From Eq. (2.17), it is clear that the imagi- as ethanol (κ 25) or water (κ 80) (Jang et al., 2008; ≈ ≈ nary part of the polarization function at small momen- Ponomarenko et al., 2009). tum is zero unless ω > 2 µ [Fig. 3(b)]. This is so be- As usual, the collective modes follow from the zeros of cause the phase space for| vertical| interband excitations the dielectric function is Pauli blocked for ω < 2 µ , generating a gap for optical absorption in the infrared.| | At finite q, the threshold for ǫ(q,ω)= ǫ [1 V (q)Π(1)(q,ω)] , (2.19) 0 − 8
ω ω
interband 2µ
intraband
q q 2kF
FIG. 5 (Color online) Colored regions represent the particle- hole continuum of graphene due to interband (gray area) and intraband (green) transitions. On the left: half-filled case; right: finite µ case, away from half filling. Dashed line: acous- tic plasmon for the single layer (ωp √µq). ∝
tric constant (Shung, 1986a),
qe2 2 µ ω (q) ǫ (q) ǫ ln | |− p . (2.21) 0 ≈ 0 − 2ω (q) 2 µ + ω (q) p | | p As in the 2DEG, the screened Coulomb interaction for q< 2kF is V (q) 1 2πe2 = (2.22) ǫ(q, 0) ǫ0 q + qT F
2 where qT F =4πe kF /(vǫ0) is the Thomas-Fermi momen- tum (kF = µ /v), which sets the size of the screening cloud. In the| | presence of an external charged impurity Ze, the induced charge, δZ, has a non-oscillatory com- FIG. 4 (Color online) Polarization bubble Π(1)(q,ω) for the ponent coming from the q 0 limit of the polarization conventional 2DEG. Panels (a) and (b) show, respectively, a 3 1 → that decays as (kF r )− (as in a 2DEG), and an oscilla- density plot of the real and imaginary parts of the polarization tory part which corresponds to the Friedel oscillations at bubble, Π(1)(q,ω), normalized to the DOS at the Fermi level. q = 2kF . The Friedel oscillations in graphene decay as Panels (c) and (d) present constant frequency cuts at ω/µ = 3 cos(2kF r)/(kF r ), differently from the 2DEG case, where 0.5, 1.0, 1.5, 2.0, 2.5, 3.0. In panel (e) we plot the static limit, 2 Π(1)(q, 0), and in (f) the uniform limit, Π(1)(0,ω). the decay is of the form cos(2kF r)/r . The difference is caused by the fact that the static polarization function in the 2DEG has a cusp at q = 2kF , whereas in graphene, the first derivative is continuous [cfr. Figs. 3(e) and 4(e)]. (1) 2 calculated here in the Random Phase Approximation For undoped graphene, V (q)Π = (π/2)[e /(vǫ0)] − (RPA). Since graphene is a 2D system, the collective [see Eq. (2.14)], and the static dielectric function is a plasmon mode is gapless. The leading term in the po- constant. The effective Coulomb interaction in this case larization for small frequency and momenta (compared is to k ) is shown in Eq. (2.15). From it one can easily F V (q) 1 2πe2 extract the infrared dependence of the plasmon, = , (2.23) ǫ(q, 0) ǫRPA q
2 2 where ǫRPA = ǫ0 + (π/2)(e /v) is the effective back- ωp(q)= (2µe /ǫ0)q , (2.20) ground dielectric constant, renormalized by the inter- p band transitions. Additional many body effects resulting which follows the same dispersion as the plasmon en- from self-energy insertions in the bubbles logarithmically countered in the 2DEG. The √q dependence of the plas- renormalize this correction to zero in the q 0 limit, mon was recently confirmed by a high resolution en- as will be clear in Sec. III of this review. On→ the dy- ergy loss spectroscopy (EELS) measurement in graphene namical side, inserting Eq. (2.14) into Eq. (2.19), one (Liu et al., 2008). Additional corrections due to the in- can easily see that no collective modes are allowed in un- terband excitations (which are absent in the 2DEG) can doped graphene, at zero temperature, within the RPA be absorbed into the definition of the background dielec- framework. At half-filling, RPA is justified in the limit 9 of large number of fermionic species, N, which favors di- with k = (q, kz), q is an in-plane momentum, and agrams with maximal number of bubbles at each order of (Fetter, 1974) perturbation theory. In graphene, the physical number of species is N = 4, and additional corrections beyond sinh(qd) S(q, kz)= (2.28) RPA coming from the exciton channel near the on-shell cosh(qd) cos(kzd) singularity of the bubble, ω vq, were shown to gener- − ate a new acoustic plasmon| mode|∼ (Gangadharaiah et al., is the structure factor for a stack with an infinite num- 2008). In the static limit (ω 0), vertex corrections in ber of layers. In the limit when the distance between the the bubble are perturbatively→ small and RPA can be jus- layers d is small, Eq. (2.27) recovers the isotropic case 2 2 2 tified in the calculation of the dielectric function even at V (k)=4π(e /ǫ0)/(q +kz ), whereas in the opposite limit 2 et al. (d ) one gets the 2D case, V (k)=2π d (e /ǫ0)/q. half-filling (Kotov , 2008b). The structure of per- → ∞ turbation theory in graphene will be discussed in detail In any case, the polarization function must be integrated in Sec. III. over a cylindrical Fermi surface of height 2π/d, and so (1) In addition to the low energy acoustic mode due to Π (q,ω) acquires an additional factor of 1/d compared intraband transitions, graphene has also two high en- to the single layer case. The extension of this problem to ergy optical plasmons generated by interband excitations include the interlayer hopping dispersion in the polariza- around the Van-Hove singularities of the π π∗ bands, tion was considered by Guinea, 2007. − Away from the neutrality point (µ = 0), instead of and also by optical transitions between σ π∗ and π σ∗ 6 bands (Eberlein et al., 2008; Kramberger− et al., 2008− ). a single acoustic mode as in the monolayer, the ze- The measured optical gaps of the π and π σ band roes in the dielectric function of the multilayer gener- plasmons in graphene are 4.5 eV and 15 eV, respectively.− ate a plasmon band, where the modes are labeled by k [ π/d,π/d]. For q 1/d, the plasmon dispersion Similar modes were also observed in graphite, where they z ∈ − ≪ appear blue shifted to 7eV and 24 eV respectively, ac- is (Shung, 1986a) cording to optical data (Taft and Philipp, 1965), X-ray 2 measurements (Shulke et al., 1988), and ab-initio calcu- 2 2µe ωp(q, kz)= qS(q, kz). (2.29) lations (Marinopoulus et al., 2004). ǫ0
In the kz = 0 mode, the charge fluctuations be- E. Infinite stack of layers tween different layers are in-phase, and the result- 2 2 ing plasmon mode is optical, ωp(q, 0) (4µe /ǫ0d) + 3 2 ≈ In the case of an infinite stack of graphene layers, the 4 (vq) . For ωp(q) > 2µ, this mode is damped by Hamiltonian term for the Coulomb interaction among all the particle-hole continuum due to interband transitions the electrons can be written in real space as (see Fig. 6), in agreement with energy loss spectroscopy data (Laitenberger and Palmer, 1996). The out-of-phase 2 e 3 3 1 modes (for k = 0) are acoustic. At the edge of the plas- C = d rd r′ nˆ(r) nˆ(r′) , (2.24) z 6 H ǫ0 r r′ mon band, the mode k = π/d disperses linearly with Z | − | z ± the in-plane momentum, ω (q, π/d) = µe2d/ǫ q, in wheren ˆ(r) is the 3D particle density operator. In the p ± 0 contrast with the 2DEG dispersion (ωp √q) present absence of interlayer hopping, as in the case for exam- ∝p ple of several graphite intercalated compounds, the elec- in the single layer. Except for the lack of an inter- trons remain confined in each layer, but the unscreened band particle-hole continuum and the associated damp- Coulomb lines fill the entire space in between the layers, ing, similar plasmon band features are also expected in coupling all the electrons in the system. In that case we the 2D layered electron gas, for fermions with quadratic may constrain the local density operatorn ˆ to be in the dispersion (Hawrylak, 1987). form (Visscher and Falikov, 1970)
∞ F. f-sum rule nˆ(r) d nˆ(r)δ(z ld) (2.25) → − l= X−∞ The f-sum rule is a generic statement about conser- where l is an integer labeling the layers, and d is the vation of the number of particles and results from the distance between layers. In momentum space, making analytical properties of the retarded charge susceptibil- a discrete sum over the layers, the Coulomb interaction ity. It can be generically defined as (Nozi`eres, 1964) between all the electrons is ∞ R 2 dωωImχ (k,ω)= π [[ , nˆ( k)] , nˆ(k)] , (2.30) e 3 h H − i C = d k nˆ( k)V (k)ˆn(k) , (2.26) Z−∞ H ǫ0 − Z where is the Hamiltonian,n ˆ is the particle density op- where erator,HχR is a retarded charge susceptibility, χ(k, τ) = e2 T [ˆn(k, τ)ˆn( k, 0)] , and ... is an expectation value cal- V (k)=2π d S(q, kz) (2.27) h − i h i ǫ0q culated in some basis. 10
ω for a Dirac Hamiltonian, where D is the ultraviolet cut- off. A similar dependence with the cut-off also occurs ω in the true 3D relativistic problem, where the sum rule p reflects the number of particles contained in the vacuum of the theory, which is formally divergent (Ceni, 2001; Goldman and Drake, 1982). In graphene, as in any two band semi-metal or semiconductor, the validity of the f- Intraband sum rule is physically recovered when the periodicity of q the electronic spectrum is restored back into the Hamil- 2kF tonian. FIG. 6 (Color online) Plasmon band (hatched region) for an infinite stack of graphene layers. Red line: optical mode III. QUASIPARTICLES IN GRAPHENE kz = 0. Dashed line: acoustic mode kz = π/d, ωp √µq, with linear dispersion, at the edge of the band. All the∝ other modes in between are acoustic. Adapted from Shung, 1986a. The quasiparticle properties of graphene are modi- fied by the presence of long-range Coulomb interactions. Their effects are especially pronounced when the Fermi As in any solid, the exact electronic Hamiltonian of energy is close to the Dirac point (µ 0), and can re- graphene can be decomposed into a Hamiltonian of free sult in strong renormalization of the Dirac≈ band struc- electrons, plus a periodic potential due to the lattice, and ture (the Fermi velocity v), and the quasiparticle residue interactions. If the interactions depend only on densities, (Z). Consequently, many physical characteristics, such the commutators in Eq. (2.30) can be calculated exactly, as the compressibility, spin susceptibility and the spe- and the only term that survives is the Kinetic energy due cific heat can be strongly affected by interactions. Even to the free electrons, when the Fermi surface is large and the system is a Fermi k2 liquid, there are strong modifications of the physics near [[ , nˆ( k)] , nˆ(k)] = Ne , (2.31) the Dirac point due to the presence of additional peaks in h H − i m the quasiparticle decay rate, related to plasmon-mediated where m is the bare electron mass and Ne is the number decay channels. Even reconstruction of the Dirac cone of fermions in the band. Choosing, for example, a basis structure near the charge neutrality point appears possi- of non-interacting fermions, the sum rule in graphene is ble, as indicated by recent Angle-Resolved Photoemission 2 Spectroscopy (ARPES) measurements. All these effects ∞ Nek dωωImΠ(1)(k,ω)= π , (2.32) are sensitive to the value of the Coulomb interaction con- m Z−∞ stant in graphene, α. as in metals, where Π(1)(k,ω) is the bare polarization bubble, calculated using the full non-interacting spec- trum (dictated by the lattice symmetry). The validity A. Low-energy behavior near the Dirac point of the f sum-rule does not require Galilean invariance of the quasiparticles, but of the free electrons, which are not 1. Weak-coupling analysis relativistic and hence obey the Schrodinger equation. For low energy effective Hamiltonians, such as the The interaction parameter which characterizes the Dirac Hamiltonian in graphene (which do not include strength of the Coulomb interaction in graphene is the periodicity of the spectrum in the Brillouin zone), (Eq. (1.4)) the f-sum rule above is still formally satisfied when ap- e2 plied for the electrons (holes) in the conduction (valence) α = . (3.1) band only, as can be explicitly checked by direct substi- ǫ0v tution of the polarization due to intraband transitions, At k = 0 screening is absent, and the interaction po- Eq. (2.15), into Eq. (2.32). The number of electrons F 2 2 tential in momentum space: (holes) in this band, Ne = kF A/π, where A = 3√3a /2 is the unit cell area, is set by the size of the Fermi surface, 2πe2 and the verification of the sum rule follows as in a Fermi V (p)= . (3.2) ǫ p liquid. 0 The Dirac Hamiltonian, however, violates the f-sum The value of α = 2.2/ǫ0 depends on the dielectric envi- rule (2.32) when interband transitions are taken into ac- ronment since, as previously discussed, ǫ0 = (1+ κ)/2 count. In that case, the left hand side of Eq. (2.32) be- for graphene in contact with air and a substrate with comes independent of the chemical potential, consistent dielectric constant κ. In vacuum, α =2.2. with the fact that (Sabio et al., 2008) In the case of small coupling, α 1, we can employ ≪ D standard perturbation theory, involving the perturbative [[ , nˆ( k)] , nˆ(k)] = k2 (2.33) h H − i 4 computation of the self-energy Σ(k,ω), which enters in a 11
a) b) A) B) ¡
c) ¡ d) FIG. 8 (a) Self-energy and (b) Vertex corrections to the po-
larization bubble ¡
¡ While the linearity of the spectrum in graphene was FIG. 7 Self-energy diagrams: (a) First order Hartree-Fock, realized a long time ago (Wallace, 1947), in the con- (b) Second order loop diagram (first diagram in the RPA se- text of studying graphite formed by layers of graphene, ries), (c) Second order exchange (vertex correction) diagram, the self-energy correction Eq. (3.8) due to interac- (d) Rainbow diagram. tions was first investigated perturbatively much later by Gonz´alez et al., 1994. The non-trivial velocity standard way the Dirac fermion Green’s function (GF), renormalization is due to the unscreened, long-range for a given valley: Coulomb interactions. Similar logarithmic divergencies were also found in gapless 3D semiconductors, where 1 the Dirac spectrum originated from special symmetries k G( ,ω)= + . ωσ0 vσ k Σ(k,ω)+ iσ00 sign(ω) (Abrikosov and Beneslavskii, 1971). − · − (3.3) The above calculation forms the basis of the Renor- malization Group (RG) analysis. In the RG spirit one It is convenient to decompose the self-energy into two integrates out the high momentum degrees of freedom, pieces with different pseudo-spin structure i.e. regions of momenta Λ > p > Λ1, and the results vary with the quantity ln(Λ/Λ |) | l. Here we denote by l 1 ≡ Σ(k,ω)=Σ0(k,ω)+Σv(k,ω), Σ0 σ0, Σv σ k, the RG parameter, so that the infrared limit corresponds ∝ ∝ · (3.4) to l (i.e. one integrates down to the infrared scale → ∞ where σ0 = 1 is the unit matrix, which from now on will k 0, l = ln(Λ/k)). From Eq. (3.9) we obtain not be written explicitly. Then we have → dv α e2 Z = v = . (3.10) G(k,ω)= , (3.5) dl 4 4ǫ0 ω Z(vσ k +Σv) − · This equation has to be supplemented with an additional where Z is the quasiparticle residue equation reflecting the absence of charge (e2) renormal-
1 ization: Z− =1 ∂Σ0/∂ω , (3.6) − de2 =0 . (3.11) and Σv is responsible solely for the velocity renormaliza- dl tion. The first order diagram shown in Fig. 7(a) is the There are several ways to understand this. It was ar- Hartree-Fock exchange contribution, and can be readily gued early on that the vertex function does not acquire evaluated (we denote by G(0) the non-interacting GF): any divergent contributions, which is related to the ex- pected regular behavior of the polarization operator to d2pdε all orders in graphene (Gonz´alez et al., 1994). More re- Σ(1)(k,ω)= i G(0)(k + p,ω + ε)V (p) , (3.7) (2π)3 cently, explicit calculations up to two loop order were Z performed (de Juan et al., 2010; Kotov et al., 2008b); it which at low external momenta exhibits a logarithmic was confirmed that the vertex function is finite in the low- singularity energy limit. In addition, direct examination of the po- larization function at two loop level (Kotov et al., 2008b) α Σ(1)(k,ω)=Σ(1)(k)= vσ k ln(Λ/k), Λ/k 1. found that the self-energy correction, Fig. 8(a), acquires v 4 · ≫ (3.8) a logarithmic divergence which can be absorbed into the renormalized velocity v(k) (Eq. (3.9)), while the vertex At this order we have Σ0 = 0, i.e. Z = 1 due to the frequency independence of the interaction potential, and correction of Fig. 8(b) is finite: the quasiparticle velocity increases: 2πe2 Π(2b)(q, 0) = finite = 0.53α2 . (3.12) α ǫ q − v(k)= v 1+ ln(Λ/k) , Λ/k 1. (3.9) 0 4 ≫ Incidentally, this contribution leads to enhancement of The ultraviolet cutoff Λ 1/a represents the momentum the dielectric static screening (i.e. the dielectric constant ∼ π 2 scale up to which the spectrum is Dirac-like. beyond linear (RPA) order becomes ǫ = 1+ 2 α+0.53 α .) 12
Alternatively, one can argue that in two-dimensional possible at αc 0.8. This fixed point is infrared unstable ≈dv field theories with Coulomb interactions the charge since near αc, dl = C(α αc)v, C> 0, i.e. for α > αc, 2 − − e does not flow because it appears as a coefficient v flows towards zero (α flows to ) while for α < αc, in a nonanalytic term in the action (Herbut, 2006; v flows towards (α flows to zero.)∞ Of course it is not Ye and Sachdev, 1998). The conclusion then is that only clear that this estimate∞ is reliable since the fixed point the quasiparticle velocity and residue (see below) are value αc is not small, and we used perturbation theory renormalized. In particular, at first order we can combine (α 1) to derive this result. On the other hand, a flow ≪ Eqs. (3.10),(3.11) into a single one reflecting the renor- towards strong coupling for α > αc is consistent with malization (running) of the coupling α: the formation of an excitonic insulator (mass generation), for which strong evidence has accumulated by now, as dα α2 = . (3.13) we discuss in Section III.B. Recent numerical simulations dl − 4 give the value αc 1 (see Section III.B). Finally, we also∼ find that Z is renormalized at second Therefore we have an infrared stable fixed point at α = 0, order, since the self-energy is frequency dependent. From and the flow towards it is logarithmic: Eq. (3.6) we can expand to second order of bare pertur- 2 4 Nα bation theory Z 1 24 ln(Λ/k), which would lead us α(k) , k 0 . (3.14) ≈ − 2 ∼ ln(Λ/k) → dZ Nα to an RG equation for Z: dl = 24 Z, to be solved together with Eq. (3.13), or Eq. (3.17− ), depending on the Thus the Coulomb interactions are marginally irrelevant. desired level of approximation. Alternatively, Eq. (3.6) This is equivalent to a logarithmically divergent velocity: is already written in a “nonperturbative” way. Ignoring v(k) (e2/4) ln(Λ/k), k 0. ∼ → for the moment the running of α, we have at low energies 1 24 a. Two-loop results. It is instructive to examine Z = 2 , k/Λ 0. 1+ Nα ln(Λ/k) → Nα2 ln(Λ/k) → corrections beyond first order (Mishchenko, 2007; 24 (3.18) Vafek and Case, 2008), since additional effects appear, (2b) such as renormalization of Z. For example the first This result, along with the previous one for Σ0 , brings diagram in the RPA series shown in Fig. 7(b) is us to the infrared behavior (we use ω and k interchange- ably in the infrared limit): d2pdε Σ(2b)(k,ω)=i G(0)(k+p,ω+ε)(V (p))2Π(1)(p,ε) . 1 (2π)3 Z , Σ α2ω ln(ω) , ω 0. (3.19) Z ∼ α2 ln(ω) 0 ∼ | | → (3.15) | | An explicit evaluation at low energies and momenta gives This is characteristic of a marginal Fermi liquid a single logarithmic divergence (Das Sarma et al., 2007; Gonz´alez et al., 1994). How- 2 ever, this regime is never achieved if the running of α is (2b) Nα Σ (k,ω)= (ω + vσ k) ln(Λ/k), k/Λ 0, taken into account, as is intuitively clear from the above − 24 · → (3.16) equations. As we will see later from the solution of the (2b) 2 (2b) RG equations for Z and α, in fact Z tends to level off in i.e. Σ0 = (Nα /24)(ω) ln(Λ/k), and Σv = (Nα2/24)vσ k−ln(Λ/k). Because the polarization bub- the infrared, and the system has well-defined quasiparti- −ble is proportional· to the number of fermion flavors cles. N = 4 (valley+spin), we have explicitly written the It is interesting to note that trigonal distortions, N dependence. By comparing with Eq. (3.5), we find which change the band structure away from the Dirac that the velocity is changed by an amount ( N/24 equation, are modified by the electron-electron interac- N/24)α2v ln (Λ/k). − − tion, and their irrelevance at low energies is enhanced In addition, other diagrams at second order have to be (Foster and Aleiner, 2008). As a result, the linear dis- added, such as the vertex correction of Fig. 7(c). Most persion becomes an even more robust feature of graphene importantly, this diagram is also proportional to ln Λ. (Rold´an et al., 2008). Collecting all contributions one finds the RG equation for the velocity flow (Vafek and Case, 2008) b. Influence of disorder. Before we proceed, let us briefly dv α N address the effect of disorder. Two major sources of = v δ α2v , (3.17) dl 4 − 12 − disorder are scalar potential random fluctuations (e.g. formation of electron-hole puddles), and vector gauge with δ 0.03. One observes that the contribution of the field randomness, related to formation of ripples. Start- “RPA”≈ diagram is numerically dominant at second order ing with the latter, i.e. a gauge field coupled to the (it is larger than the rest by a factor of 10 for N = 4.) Dirac fermion pseudospin σ A, and characterized by · In addition, the second order tendency is a decrease of variance ∆, Aµ(r1)Aν (r2) = ∆δµν δ(r1 r2), one can the velocity. Consequently a finite coupling fixed point is readily deriveh the correspondingi RG equations− in the 13
α a)
b)
= ¡+
¡ ¡
FIG. 10 RPA self-energy,¡ which includes an infinite resum- mation of polarization bubbles.
The RPA potential is given by Disorder ∆ 2πe2 RPA p V ( ,ε)= 2 (1) . (3.22) FIG. 9 An attractive line of fixed pints for interactions and ǫ0p 2πe Π (p,ε) gauge field disorder. − Quite remarkably, at low momenta one can evaluate the singular contribution to the self-energy analytically weak disorder and interactions limit (Herbut et al., 2008; Stauber et al., 2005) 8 Σ(RPA)(k,ω)= [ F (λ)ω + F (λ)vσ k] ln(Λ/k) , Nπ2 − 0 1 · d∆ dα α2 ∆ =0, = + α . (3.20) (3.23) dl dl − 4 π where we have defined Gauge field disorder itself is not renormalized, while the π λ = Nα. (3.24) interplay of disorder and interactions leads to a line of 8 4 attractive fixed points located at: α∗ = π ∆, as shown in Fig. 9. Physically the variance is related to the char- This parameter is measuring the importance of polariza- acteristic height h, and length L of the corrugations of tion loop contributions relative to the bare Coulomb term the surface, ∆ h4/(L2a2). Thus weak disorder generi- (i.e. the ratio of the second term to the first in the de- cally shifts the∼ fixed point away from α = 0, while strong nominator of Eq. (3.22)). The RPA is generally expected disorder can have an even more profound effect (Section to be valid when the loops dominate over other diagrams, i.e. N 1. Provided this condition is satisfied, we can VI.C). ≫ In addition, for weak interactions, the inclusion of also analyze the strong-coupling regime λ 1, and the crossover toward the weak-coupling one (λ≫ 1), i.e. we scalar (density fluctuations) disorder turns out to be ≪ a relevant perturbation which grows under renormal- can hope to cover a wide range of α values. ization, and thus away from the perturbative regime The calculated functions F0 and F1 in Eq. (3.23) are (Aleiner and Efetov, 2006). Moreover, gauge field disor- √1 λ2 π der, when combined with strong-enough interactions, can − arccos λ 1+ , λ< 1, F (λ)= − λ − 2λ cause the interactions to grow (Vafek and Case, 2008). It 1 √λ2 1 π has been argued that the strong-coupling regime for dis- − ln λ + λ2 1 1+ ,λ> 1, order and interactions generically occurs when all types λ − − 2λ p (3.25) of disorder consistent with graphene’s symmetries are in- cluded (Foster and Aleiner, 2008). 2 λ2 π − arccos λ 2+ , λ< 1, A detailed analysis of this complex situation is beyond − √ 2 − λ F (λ)= λ 1 λ the scope of this work, and from now on we continue our 0 λ2 −2 π discussion of clean graphene. − ln λ + λ2 1 2+ ,λ> 1. λ√λ2 1 − − λ − p (3.26) This leads to the system of RG equations for v and Z, to 2. Strong-coupling/RPA analysis leading order in 1/N The full RPA treatment was performed by many au- dv 8 et al. = (F1(λ) F0(λ)) v , (3.27) thors (Das Sarma , 2007; Foster and Aleiner, 2008; dl Nπ2 − Gonz´alez et al., 1999; Kotov et al., 2009; Polini et al., 2007; Son, 2007). Here we mostly follow Son, 2007. The dZ 8 RPA self-energy is shown diagrammatically in Fig. 10, = F0(λ)Z. (3.28) and corresponds to the equation dl −Nπ2 d2pdε At strong-coupling, λ 1, one finds Σ(RPA)(k,ω)= i G(0)(k + p,ω + ε)V RPA(p,ε) . ≫ (2π)3 dv 8 Z = v , (3.29) (3.21) dl Nπ2 14
dZ 8 5 = ln (2λ)Z. (3.30) dl −Nπ2 4 The first equation, after integration, leads to the low- energy result (k 0) → 3 Λ η 8 λ v(k)/v = , η = , (3.31) 2 k Nπ2 1 which implies that the quasiparticle dispersion is of the form 0 0 2 4 6 8 8 ω(k) kz, z =1 . (3.32) l ∼ − Nπ2 The existence of the anomalous velocity dimension, η, 1 and consequently z = 1, is characteristic of the strong- coupling regime Nα6 (Son, 2007). However this 0.8 strongly-coupled fixed→ point ∞ is infrared unstable, since, due to the velocity increase, the RG for α flows towards 0.6 weak coupling. (One also expects that for certain N < Z 0.4 Nc and α 1 an excitonic gap can appear, which will be discussed≫ in Section III.B.) In this regime Z can be approximated perturbatively (in 1/N) as 0.2
8 1 0 Z 1 ln(Nαπ/4) ln(Λ/k), Nα 1, (3.33) 0 5 10 15 20 25 30 ≈ − π2 N ≫ l which can be obtained from Eq. (3.30) by ignoring the FIG. 11 RG flow of the coupling λ and the quasiparticle scale dependence of λ. residue Z as a function of the RG scale l; the infrared limit is at l . From Gonz´alez et al., 1999. In the weak-coupling limit λ 1, it is easy to verify →∞ that we recover the previous result≪ (3.10) for the velocity v (leading to a flow for α towards zero), and the previ- ously encountered perturbative equation for Z At higher energies however (away from the fixed point but still much lower than the bandwidth vΛ), the system dZ 8 λ2 N exhibits marginal Fermi liquid behavior. = Z, i.e. Z 1 α2 ln(Λ/k). (3.34) dl −Nπ2 3 ≈ − 24 At finite (but still small) density away from the Dirac point, i.e. k = 0, the logarithmic behavior in The last formula is written to first order in Nα. the infrared is cut-off6 by the Fermi momentum, i.e. Eqs. (3.33),(3.34) allow us to have a qualitative under- ln(Λ/kF ), kF /Λ 0, and the RG stops away from standing of the behavior of Z as a function of the RG the fixed point. For→ comparison with experiments, the scale l. If the initial value of α is large, at the initial flow toward this stable fixed point should be stopped at RG steps Z decreases logarithmically fairly fast (due to a scale set by the (small) density, temperature, or fre- the weak ln(α) dependence in Eq. (3.33), even though α quency, whichever is higher. itself decreases). Eventually, when α has decreased sub- One can also perform a numerical evaluation of the 1 stantially (α (ln(Λ/k))− ), Z is governed by Eq. (3.34), main RPA equation Eq. (3.21) (Polini et al., 2007). For meaning that∼Z will stop decreasing, and will level off for small density, and with logarithmic accuracy (ln(Λ/kF )), l = ln(Λ/k) . → ∞ this is equivalent to evaluating, by using the notation of A numerical evaluation of the system of equations Eq. (3.23), and taking into account Eqs. (3.4),(3.5),(3.6) (3.27),(3.28) confirms the anticipated behavior and is shown in Fig. 11, (Gonz´alez et al., 1999). (The equa- 1 Z = (1 ∂Σ(RPA)/∂ω) 1 = , tion for the coupling λ = π Ne2/(ǫ v) is obtained by − 8 8 0 − 1+ 2 F0(λ) ln(Λ/kF ) 2 π 2 Nπ observing that (dλ/dl) = ( 1/v ) 8 N(e /ǫ0)(dv/dl), due (3.35) to the lack of charge renormalization.)− We conclude that the flow of λ is towards weak coupling, no matter how 8 v∗/v = Z 1+ F (λ) ln(Λ/k ) . (3.36) large its initial value is. Z does not renormalize to zero at Nπ2 1 F low energy due to the RG decrease of λ. Thus, near the weak-coupling infrared fixed point, the marginal Fermi Here v∗ is the renormalized velocity. At any finite den- liquid (Eq. (3.19)) is ultimately not reached, and the sity the numerical evaluation of Σ(RPA) also picks up system behaves as a Fermi liquid (although the quasi- finite (subleading) contributions, while it can be shown particle decay rate is non-Fermi liquid like, see below.) (Polini et al., 2007) that the leading perturbative results 15
1 (a) (b) Λ = 101 Λ = 102 Λ = 103 4 0.8 Λ = 10 Λ = 105 2DES Z
0.6
(a) 0.4 0 1 2 3 4 5 f FIG. 13 (Color online) (a) Density dependence of the velocity 1.8 for suspended graphene, from (Elias et al., 2011). The solid (b) line is the result of RG treatment within RPA (Eq. (3.27)). 1.6 (b) Reshaping of the Dirac cone due to the interaction-driven renormalization (increase) of the Fermi velocity at low mo- 1.4 menta. The outer cone represents the linear Dirac spectrum without many-body effects. /v ⋆
v 1.2 logarithmic renormalization of the velocity predicted by 1 theory fits the data fairly well, and thus offers a direct proof that the Dirac cones can be reshaped by long-range 0.8 electron-electron interactions near the Dirac point, as 0 1 2 3 4 5 schematically shown in Fig. 13(b). Finally, ARPES mea- f surements of quasi-freestanding graphene grown on the FIG. 12 (Color online) Exact evaluation of the RPA equations carbon face of SiC have also detected logarithmic velocity for (a) the quasiparticle residue, and (b) the Fermi velocity. renormalization (Siegel et al., 2011). On the horizontal axis f is defined as f Nα. Λ is in units 2 1 ≡ of kF . Values of Λ from 10 to 10 correspond to density 11 −2 ∼ 13 −2 3 n from n 10 cm to n 10 cm (while Λ 10 is 3. Quasiparticle lifetime ultra low density∼ n 109cm∼−2). The values of Λ (in∼ units ∼ of kF ) can be converted into density n via: Λ/kF 220/√n˜, 10 −2 ≈ The inverse quasiparticle lifetime (decay rate) due n˜ = n/(10 cm ). The curves labeled 2DES refer to the to electron-electron interactions, 1/τ , is an important √ ee case of 2DEG with parabolic bands, where f = 2rs, and quantity which is relevant to many properties of graphene rs 1/√n. From (Polini et al., 2007). ∼ (and Fermi systems in general). In particular the depen- dence of 1/τee on energy (or temperature) determines the importance of the electron-electron interaction con- such as Eqs. (3.33),(3.34) are readily reproduced. The tribution, relative to other processes, to transport, and RPA results are shown in Fig. 12, and exhibit the natu- interpretation of spectroscopic features, such as ARPES. ral density dependence tendency, i.e. the strongest renor- The decay rate is determined by the imaginary part of malization occurs at the lowest densities. Similar RPA the self-energy, ImΣ(k,ω). The first diagram which has results have been obtained by Das Sarma et al., 2007. energy dependence, and thus a non-zero imaginary part, A significant velocity enhancement was observed in the is the one bubble diagram of Fig. 7(b), whose real part infrared conductivity (Li et al., 2008), which reported is given by Eq. (3.16), i.e. behaves as in Eq. (3.19) at around 15% increase of the Fermi velocity, having value low energies. We can therefore deduce, for energies and 6 as high as v∗ 1.25 10 m/s at the lowest densities momenta close to the mass shell (Gonz´alez et al., 1996), (compared to v≈ 1.1 ×106m/s at higher density). The ≈ × system is at a finite Fermi energy µ 0.2eV. However ImΣ(2b)(k,ω) α2θ(ω vk) ω, ω vk , (3.37) the velocity renormalization is not logarithmic,≈ and it is ∼ − ≈ not clear what is the origin of this effect. i.e. the decay rate is linear in energy. In addition, there A recent study of suspended graphene which measures is an on-shell (“light cone”, ω = vk) discontinuity, where the cyclotron mass (Elias et al., 2011) has detected sig- the rate experiences a jump. This on-shell behavior is nificant logarithmic renormalization of the Fermi veloc- due to the fact that, for ω
The above behavior is valid at the Dirac point and T = by confining the electrons into finite-size configurations, µ = 0, while for small T,µ, it is valid for energies of order such as quantum dots (Ponomarenko et al., 2008). In max(T,µ). Notice also that the linear energy behavior of these cases the gap generation mechanism is not intrinsic Eq. (3.37) is very different from the conventional Fermi to graphene, and the value of the gap depends strongly on liquid result ImΣ ω2 (Das Sarma et al., 2007), which the external conditions. However even in such situations would occur for a finite∼ Fermi surface (µ = 0) and is due interactions can play an important role by increasing the to intra-band particle-hole excitations. 6 gap. The on-shell discontinuity present at the one-loop Consider a gap arising from an external potential that level Eq. (3.37) disappears when the full RPA self- alternates between the two sublattices energy is evaluated (Fig. 10). In this case one obtains (Khveshchenko, 2006) mass = ∆0 nσ(Ri) ∆0 nσ(Ri) . (3.40) H − σ,i A σ,i B X∈ X∈ ImΣ(RPA)(k,ω) ln (πα)θ(ω vk)(ω vk), ω vk . ∼ − − ≈(3.38) Consequently an additional pseudospinor structure re- Away from the mass shell, the energy dependence is nat- lated to σ3 is generated, and the new Green’s function urally linear: has the form 1 ImΣ(RPA)(k,ω) ln (πα) ω, ω vk . (3.39) G(k,ω)= . (3.41) ∼ ≫ ω vσ k ∆0σ3 Σ(k,ω) − · − − (RPA) k The full dependence ImΣ ( ,ω) has to be evalu- Here ∆ is the explicit “mass” of the graphene elec- ated numerically (Das Sarma et al., 2007), and the re- 0 (RPA) trons (while Σ(k,ω) contains the information about sults confirm the smooth rise of ImΣ from the point interactions, assumed to perturbatively renormalize ω = vk. all the other terms.) The new spectrum is then In the limit of zero doping µ 0, when the system ap- E(k) = v2k2 + ∆2, with a gap of 2∆ . Com- proaches the fixed point α = 0,→ we argued previously that 0 0 puting the± Hartree-Fock interaction correction to ∆ the residue Z does not approach zero (i.e. the marginal 0 leads to a renormalizedp mass ∆˜ (Kane and Mele, 2005; Fermi liquid behavior ultimately does not manifest itself.) 0 Kotov et al., 2008a) On the other hand the marginal Fermi-liquid behavior is expected to be much more robust as far as the inverse α ∆˜ /∆ 1+ ln (D/∆ ). (3.42) lifetime, ImΣ ω, is concerned, because the running of 0 0 ≈ 2 0 the coupling α∼(ω) only introduces logarithmic variation on top of a much stronger linear energy dependence. The above enhancement can be substantial. For ex- The linear decay rate discussed above is consis- ample for a bare gap due to spin-orbit coupling ∆0 3 ∼ tent with ARPES experiments (Bostwick et al., 2007; 10− meV (Min et al., 2006; Yao et al., 2007), and taking Zhou et al., 2008), and STM measurements of graphene into account the bandwidth D = vΛ 7eV, the logarith- ≈ on graphite (Li et al., 2009a) (see also the discussion in mic factor is around 15. In fact one should integrate the ˜ Grushin et al., 2009). RG equation for the renormalized mass ∆0 as a function of ln(Λ) simultaneously with the equation for the running coupling α(ln(Λ)), Eq. (3.13), down to the lowest infrared B. Spontaneous mass generation scale ∆0 (bare gap). This leads to the stronger de- ∼ ˜ α 2 pendence ∆0/∆0 = (1+ 4 ln (D/∆0)) (Kane and Mele, It is an intriguing possibility that graphene can un- 2005), and the perturbative expansion of this result is dergo a metal-insulator transition for strong enough Eq. (3.42). It is interesting to note that the logarithmic Coulomb interaction α, due to an excitonic pairing mech- mass renormalization formula in graphene Eq. (3.42) is anism. We restrict ourselves to the charge neutrality similar to the well-known expression for the electromag- point µ = 0 since the excitonic pairing tendency de- netic mass of the electron (accounting for radiative cor- creases quickly beyond that. rections) in 3D relativistic QED (Weisskopf, 1939).
1. Finite explicit mass 2. Excitonic mass generation
Before we outline the main results, let us mention that We now turn to the possibility of spontaneous gap gen- an explicit gap can also open in graphene under certain eration due to long-range Coulomb interactions (we set conditions that depend on graphene’s environment. For the explicit gap ∆0 = 0 in Eq. (3.41)). In relativistic example there are suggestions of a detectable gap in situ- QED in two space (plus one time) dimensions, QED2+1, ations when graphene is on a substrate with specific sym- the study of this phenomenon, called chiral symmetry metry, creating sublattice asymmetry in the graphene breaking, started quite a while ago (Appelquist et al., plane, and thus making the graphene electrons massive 1986; Pisarski, 1984), and is still going strong today. (gapped) (Zhou et al., 2007). Gaps can also be produced Graphene is actually different from QED2+1 because only 17 α = 8 (Khveshchenko, 2009). At strong coupling α the gap is non-zero only below a critical number of→ fermion ∞ Excitonic flavors (since the effective interaction scales as 1/N in this limit); for example N 7.2 (Khveshchenko, 2009), α Insulator c ≈ Nc 7 (Liu et al., 2009). Near≈ the critical coupling the low-momentum gap α scales as c C ∆(0) D exp , (3.44) ∝ −√α α Gapless semimetal eff − eff,c where C is a constant, the critical αeff,c = 1/2, and the form of the effective coupling αeff depends on the level of approximation used — for example an improve- ment over the static RPA potential leads to: αeff = 0 4 N N α/(1 + Nπα/8√2) (which gives N 7.2, α 1, and c c ≈ ≫ αc = 1.13, N = 4 (Khveshchenko, 2009)). The form of FIG. 14 Schematic phase diagram in the α N plane. Eq. (3.44) suggests that the transition is of infinite or- − der (Berezinskii-Kosterlitz-Thouless type). Even though Eq. (3.44) is only valid near the critical coupling, numer- the fermions are confined to a 2D plane, while the field ical results find that the gap in units of the bandwidth, lines extend through the whole 3D space. In addition, the ∆(0)/D, is exponentially small in a wide range of cou- Coulomb interaction in graphene can be considered in- plings (Khveshchenko and Leal, 2004). Since D 7 eV, stantaneous since the speed of light c is much larger than this implies ∆(0) meV, i.e. a rather small gap≈ value. ∼ the Fermi velocity (v c/300). Hence, Lorenz invari- Finally, recent work that takes into account the renor- ance is not respected, which≈ reflects the non-relativistic, malization of the coupling constant and the quasiparticle purely band origin of the Dirac quasiparticles. The anal- residue suggests that αc could be much larger than pre- ysis in relativistic QED reveals that dynamical mass can viously found (Gonz´alez, 2010; Sabio et al., 2010a). be generated below a critical number of fermion flavors The above results are based on various approxima- Nc, with the mass scale set by the coupling itself, which tion schemes and it is therefore important to com- has dimension of energy in pure QED2+1. A transition is pare them with direct numerical simulations of the lat- also found in non-relativistic graphene, where the gener- tice field theory model. Recent Monte Carlo calcula- ated mass scale is related to the ultraviolet energy cutoff tions (Drut and L¨ahde, 2009a,b,c; Hands and Strouthos, (bandwidth D = vΛ) since the coupling α is dimension- 2008) provide strong evidence that spontaneous mass less in this case. generation does occur, and give comparable values The gap equation can be obtained as a self-consistent for the critical couplings: N 9.6, α 1 c ≈ ≫ solution for the self-energy within RPA (i.e. vertex (Hands and Strouthos, 2008), αc = 1.1, N = 4 corrections are neglected), and is referred to as the (Drut and L¨ahde, 2009b). Unfortunately the Monte Schwinger-Dyson equation. It has the form Carlo simulations do not allow for an exact determina- tion of the gap size, and for that we can only rely on the d2kdω V RPA(p k,ε ω)∆(k,ω) previously described Schwinger-Dyson equation (leading ∆(p,ε)= i − − . (2π)3 ω2 v2k2 ∆2(k)+ i0+ to small gaps). For graphene deposited on SiO the value Z − − (3.43) of αSiO2 0.79 and is therefore not enough to generate The structure of the solution has been analyzed ex- a gap; only≈ experiments on ultrahigh mobility suspended tensively (Gamayun et al., 2010; Gorbar et al., 2002; samples can potentially reveal the insulating state. Khveshchenko, 2009; Khveshchenko and Leal, 2004; The overall phase diagram of graphene in the α N Liu et al., 2009) at different levels of approximation. plane is expected to look as shown in Fig. 14,− with The equation is simplified significantly if the static RPA αc 1 and Nc 7 9. At finite temperature one RPA ≈ ≈ − potential is used V (p, 0) (Khveshchenko and Leal, expects the existence of a critical temperature Tc 2004), while the dynamical equation has also been stud- ∆(0), while finite doping µ very quickly destroys the∼ ied on-shell (∆(p,ε = vp)) (Khveshchenko, 2009), as well gap (Liu et al., 2009). Application of magnetic field as numerically (Liu et al., 2009). perpendicular to the graphene layer leads to enhance- The mass gap ∆(p) has strong momentum dependence, ment of the excitonic instability due to the formation of due to the long-range nature of the Coulomb interaction. Landau levels (Gorbar et al., 2002; Gusynin et al., 2006; ∆(p) decreases at large momenta and reaches maximum Khveshchenko, 2001a). In addition, it has been suggested value at small momenta where it levels off. For fixed that an in-plane magnetic field favors a gapped excitonic physical value of N = 4, a transition to a gapped state is state (Aleiner et al., 2007), due to the instability of a found above a critical coupling αc. Some of the calculated system of electrons and holes polarized in opposite direc- values are: αc =0.92 (Gamayun et al., 2010), αc =1.13 tions. 18
power of energy (or temperature), as in a Fermi liquid, while the quasiparticle residue is finite at the Fermi sur- face.
The existence of a plasmon-related peak in the quasi- particle decay rate, which originates from intraband transitions in which an electron can decay into a plas- mon, was pointed out in the context of intercalated graphite, where the physics is dominated by graphene layers (Lin and Shung, 1996; Shung, 1986b). For n-doped graphene (µ > 0), which is relevant to ARPES ex- periments, a double-feature is found in the decay rate ImΣ: a peak at positive energies, signaling an on- set of plasmon emission, and a sharp spectral feature at negative energies, below the Dirac point, and sepa- rated from it by an amount proportional to the plasmon frequency (Hwang and Das Sarma, 2008b; Polini et al., 2008a). This is the so-called “plasmaron” — a resonance which consists of a quasiparticle strongly coupled to plas- mons (Lundqvist, 1967). Plasmaron features have been previously detected for example in optical measurements FIG. 15 (Color online) ARPES data from (Bostwick et al., of Bismuth (Tediosi et al., 2007). 2010), showing strong features at the Dirac point, which is below the Fermi energy (at 0). The splitting shown in (H) is The above calculations were done within RPA theory. attributed to the presence of “plasmarons” — quasiparticles Line widths have also been analyzed via ab-initio many- strongly bound to plasmons — and depends on the value of body methods (Park et al., 2009; Trevisanutto et al., α (α 0.5 fits the data.) ≈ 2008). Experiments generally show a well-pronounced linear quasiparticle spectrum (Bostwick et al., 2007; Sprinkle et al., 2009; Zhou et al., 2007, 2008), with ad- The physical structure of the gapped state depends on ditional features near the Dirac point which seem to de- the nature of pairing between the valleys — for example pend on the way graphene is prepared, and its purity. one can have charge density wave states (Khveshchenko, For example, gap-like features have been observed near 2001b) with modulation of the electronic density around the Dirac point (Zhou et al., 2007), and attributed to the two sublattices (which corresponds to intravalley par- external, substrate-related factors. Bending of the Dirac ing), or Kekule dimerization (Hou et al., 2007) which spectrum (kink-like feature) was attributed to plasmons corresponds to tripling of the unit cell (intervalley pair- (Bostwick et al., 2007). Most recently manifestations of ing). One generally expects that interactions beyond the sharp plasmaron spectral intensities have been ob- the long-range Coulomb potential, such as short-range served in quasi-freestanding graphene (Bostwick et al., repulsion, would favor particular states, including time- 2010), where a reconstruction of the Dirac point crossing reversal symmetry broken (spin) states. Further discus- seems to take place, as shown in Fig. 15. A diamond- sion appears in Section V.A. like shape appears due to crossing of charge and plas- maron bands. Comparison of the RPA calculation for the energy splitting with experiment leads to the value of C. Finite density Fermi-liquid regime α 0.5 (Fig. 15.) Bostwick et al., 2010 also suggest that the≈ plasmaron features were obscured in earlier measure- As the density increases above half-filling, i.e. ments on non free-standing graphene (Bostwick et al., graphene is at a finite, not necessarily small, chemical po- 2007), due to the several times stronger screening (and tential µ, with a finite Fermi surface, a crossover towards consequently smaller α.) Perhaps most importantly, all a Fermi liquid regime takes place. In this case the lower the current activity in ARPES on different graphene sam- (hole) band becomes irrelevant and the physics near the ples reveals that the electron-electron interactions can af- Fermi surface is dominated by intra-band transitions in fect strongly the physics around the Dirac point, even for the conduction (upper) band (assuming µ> 0). However relatively large density (Fermi energy). the physics near the Dirac point can still be very strongly affected due to the presence of plasmon and “plasmaron” Tunneling spectroscopy measurements, combined with features in the quasiparticle spectral function. ab-initio calculations, have also found evidence for The quasiparticle width near kF is quite simi- density-dependent interactions effects in the tunneling lar to the case of an ordinary 2D electron gas current (Brar et al., 2010) which arise from the sharp (Das Sarma et al., 2007; Hwang and Das Sarma, 2008b; spectral features in the quasiparticle decay rate below Polini et al., 2008a), and is proportional to the second the Dirac point, as discussed above. 19
this gives (1/κ0)= v π/(4n) — behavior which can be clearly seen in experiment Fig. 16. The interaction effectsp in the ground state energy ac- quire divergent contributions in the limit of small den- sity kF /Λ 0, similarly to the previously discussed self- energy (velocity)≈ renormalization. Ignoring any finite (non-diverging) terms, one finds (Barlas et al., 2007) α E /n = (vk ) ln(Λ/k ), (k /Λ) 0, (3.45) ex 6 F F F →
Nα2 E /n = G(α)(vk ) ln(Λ/k ), (3.46) RPA − 6 F F where the function G(α) is defined as G(α) = ∞ 2 2 √ 2 1 (1/2) 0 dx(1+x )− ( x +1+Nπα/8)− , and, in par- FIG. 16 (Color online) Inverse compressibility, measured by ticular, at zero coupling G(0) = 1/3. The above results Martin et al., 2008. The red line is the compressibility of non- exactlyR follow the velocity renormalization, i.e. are equiv- interacting Dirac fermions. alent to the substitution v v(k ) in the free compress- → F ibility (1/κ0) = v π/(4n), where v(kF ) is the running velocity calculated within RPA at the infrared scale kF . D. Physical Observables The result is particularlyp simple at the Hartree-Fock (ex- change) level (when the velocity follows Eq. (3.9)): The interaction-driven singular logarithmic structure near the Dirac point (for µ 0) encountered in the 1 π α ≈ = v 1+ ln(Λ/k )+ O(α2) , (3.47) fermion self-energy, and in particular the renormaliza- κ 4n 4 F tion of the Fermi velocity, can manifest itself in numer- r ous physical observables, such as the charge compressibil- and was obtained by a number of authors (Barlas et al., ity and the spin susceptibility, which exhibit non Fermi- 2007; Hwang et al., 2007; Sheehy and Schmalian, 2007). liquid behavior. Interactions can also affect the conduc- The above results are valid at zero temperature. We tivity near the Dirac point, leading to deviations from also point out that exactly at zero density kF = 0, but 2 1 2 the celebrated quantized value σ0 = e /4~ expected for T = 0, the compressibility behaves as: κ− (v /T )(1+ 6 2 ∼ free Dirac fermions (Castro Neto et al., 2009a). (α/4) ln(T0/T )) , where T0 is the temperature related to the ultraviolet cutoff; since Λv 7 eV, then T0 8 104K. This is easily understood≈ since in the infrared≈ 1. Charge and spin response limit× near the “critical point” n = T = 0 it’s the larger scale, either vkF , or the temperature T , which enters the a. Compressibility. First we discuss the compressibility physical observables (Sheehy and Schmalian, 2007). κ, which was recently measured (Martin et al., 2008), Of course Eqs. (3.45),(3.46) are valid only asymptoti- Fig. 16, and it was concluded that no interaction effects cally (kF 0), and at any finite density the compressibil- were clearly visible in those samples. Theory predicts ity should→ be calculated numerically. This was achieved significant (α dependent) deviations from the free elec- by expressing the ground state energy via the charge re- tron behavior (Barlas et al., 2007; Hwang et al., 2007; sponse function (Barlas et al., 2007). Polini et al., 2008b; Sheehy and Schmalian, 2007). Fig. 17, upper panel, illustrates the variation of 1/κ The computation of the compressibility requires knowl- with density for fixed interaction. Most notably, 1/κ edge of the ground state energy, which contains the first is larger than the free value 1/κ0. Also, the full order Hartree-Fock exchange contribution Eex, and the RPA implementation weakens the first order Hartree- correlation energy ECorr, describing all the higher order Fock (exchange) result, due to the different signs in effects. Keeping in mind applications of the theory for Eqs. (3.45),(3.46). For example, at α = 0.8 the RPA fairly strong coupling (α 1), the contribution of ECorr term is approximately 1/2 of the exchange, and thus has can be substantial. The correlation∼ energy can be readily to be taken into account (although the RPA effects be- 1 1 calculated within the RPA approximation, i.e. we take come weaker for α 0). Asymptotically, (κ− /κ0− ) E = E . The total ground state energy E, per ln(Λ/k ), as k /Λ → 0. The lower panel gives the vari-∼ Corr RPA F F → unit area, is the sum E = Ekin +Eex+ERPA. The kinetic ation κ/κ0 as a function of the interaction for different 2 energy Ekin = (2/3)vkF n, and n = (kF ) /π is the parti- densities; naturally the deviation from the free limit in- cle density. The inverse compressibility is then calculated creases with increasing interaction and decreasing den- as 1/κ = ∂2E/∂n2, which is equivalent to the usual def- sity. inition involving the variation of the chemical potential The increase of the inverse compressibility, κ0/κ, as with density, (1/κ) = ∂µ/∂n. For free Dirac particles a function of the interaction α (at fixed density), and 20
1.5 5 of the full RPA analysis led us to conclude that α< 0.1. 4 α = 0.8 It has also been argued that exchange and correlation ef- ] 2 RPA 3 fects vanish and do not manifest themselves at all in the 1 (leading logs) 2 compressibility (Abergel et al., 2009). These discrepan- cies indicate that the issue is still unsettled, while it’s eV cm 1
-12 0 also possible (indeed, quite probable) that interaction ef- 0.5 0.001 0.01 0.1 1 [10 fects are obscured by charge inhomogeneities (electron-
κ) Hartree-Fock hole puddles) in these samples. Nevertheless theory pre-
(1/ dicts strong systemic (albeit logarithmic) deviations from 0 Kinetic energy Fermi-liquid theory, and it would be important to test these predictions in cleaner, more uniform, high-mobility, 0.1 0.2 0.3 0.4 0.5 0.6 12 -2 low-density samples. density n [10 cm ]
1 b. Spin susceptibility. The paramagnetic spin suscepti- bility, χs, shows behavior very similar to the charge com- 0.9 1 10 pressibility, i.e. (χs/χs,0) decreases as the interaction in-
0.8 2 creases (Barlas et al., 2007). This is again related to the 10 fact that χ 1 is calculated via the ground state energy, 0 s− 0.7 and is proportional to the Fermi velocity v. It was also 3 κ/κ 10 pointed out that the same effect, i.e. the logarithmic 0.6 growth of the exchange energy, Eq. (3.45), can lead to
4 suppression of ferromagnetism in graphene at low densi- 0.5 10 ties (Peres et al., 2005). The full calculation of χs within 0.4 RPA was carried out by Barlas et al., 2007. 0 1 2 3 4 5 Nα On the other hand the orbital diamagnetic susceptibil- 2 ity, χdia, is proportional to v , because the quasiparti- FIG. 17 Upper panel: Inverse compressibility calculated at cle current that couples to the vector potential contains different levels of approximation as a function of density. The v (the magnetic field is perpendicular to the graphene inset enlarges the low-density region. Lower panel (adapted plane). Therefore interaction corrections lead to an in- from Barlas et al., 2007): Compressibility calculated within RPA, relative to the free level for different couplings and den- crease of χdia (Sheehy and Schmalian, 2007) and, conse- sities. Here N = 4 is the Dirac fermion degeneracy. The quently, orbital effects are expected to dominate in the numbers refer to the values of Λ/kF , which can be converted susceptibility. At the Dirac point, kF = 0, one finds at 10 −2 into density n via: Λ/kF 220/√n˜,n ˜ = n/(10 cm ). This finite temperature 2 ≈ 11 −2 implies (Λ/kF ) 10 for n 10 cm , and (Λ/kF ) 10 for − n 1013cm 2. ∼ ∼ ∼ α 2 ∼ χdia/χdia,0 = 1+ ln(T0/T ) , (3.48) 4 where the non-interacting χ = e2v2/(6πc2T ) with decreasing density (for fixed interaction), represents dia,0 − non-Fermi liquid behavior, and reflects the lack of screen- (Ghosal et al., 2007). Here c is the speed of light. At T = 0, n = 0, we have χ e2v/(c2√n), and interaction ing. By contrast, in a 3D (and 2D) Fermi liquid with 6 dia,0 ∼− a screened potential κ0/κ decreases; for example within corrections readily follow from the v dependence. This re- sult is, strictly speaking, valid for T B µ = v√πn, Hartree-Fock, κ0/κ 1 rs/6 < 1, and eventually goes ≪ ≪ through zero, signaling≈ an− instability (Mahan, 2000) (al- whereas for B = 0 the orbital susceptibility is zero for µ =0as T 0, and is finite only when the Fermi energy though the critical value of rs depends strongly on the 6 → level of approximation.) Such an instability does not oc- is at the Dirac point. It has been suggested that an in- cur in graphene, which is related to the impossibility of teraction driven positive (paramagnetic) contribution to Wigner crystallization (Dahal et al., 2006). It should be the orbital susceptibility can therefore become dominant in doped graphene, χ [e2v2/(µc2)]α ln α , α 1 noted that for larger densities (larger than the density orb ∼ | | ≪ range shown in Fig. 17) the logarithmic corrections be- (Principi et al., 2010). come unimportant and the system recovers the Fermi liq- uid behavior, i.e. eventually κ/κ0 becomes larger than 1. c. Specific heat. The specific heat is logarithmically sup- 2 Fits of the experimental data for κ with adjusted pressed due to the suppression of the DOS v− . Con- 6 2 ∼ (slightly larger) velocity v = 1.1 10 m/s show that sequently CV CV,0/(ln(T0/T )) ,T/T0 1, where × 6 2 2 ∼ ≪ α 0 (Fig. 16), while the use of v = 10 m/s by CV,0 T /v is the free Dirac fermion specific heat. Sheehy≈ and Schmalian, 2007 at the Hartree-Fock level The full∼ RPA calculation, valid also for large coupling, produced α 0.4. On the other hand, the application was carried out by Vafek, 2007. ≈ 21 d. Graphene as a quantum critical system. A uni- ductivity should have the form fied view of the above behavior is presented in Sheehy and Schmalian (2007), where it was stressed that Cα˜ σ(ω)/σ0 =1+ α , (3.49) the logarithmic corrections are manifestations of scal- 1+ 4 ln (Λv/ω) ing behavior around the quantum critical point at n = 0,T = 0. As discussed previously, at finite chemical where the constant C˜ 0.01, as argued by Mishchenko, potential, T = 0,n = 0, graphene behaves as a Fermi 2008; Sheehy and Schmalian≈ , 2009. The smallness of liquid, whereas at T6 = 0, a quantum critical region C˜ reflects the near cancellation of self-energy and ver- fans out of the point n6 = 0,T = 0. In the critical re- tex corrections, and thus the effect of interactions is gion it is natural to call graphene a Dirac liquid, where small. This value is also consistent with optical measure- the proximity to the Dirac point is important for phys- ments on suspended samples (Nair et al., 2008), as well ical phenomena at finite T . This puts graphene’s be- as graphene on a substrate (Li et al., 2008), which find havior into the general framework of quantum critical σ(ω) to be very close to σ0, and frequency independent phenomena (Sachdev, 1999). In practical terms, it im- in a wide range of energies. plies that the logarithmically divergent velocity contri- In the strict DC limit ω = 0, the presence of disorder, butions are cut-off by the largest scale: temperature T , in combination with interactions, can alter the conduc- kF √n, or magnetic field. Computing physical quanti- tivity. For example, for weak gauge field disorder (∆) ties∼ in perturbation theory (Hartree-Fock or RPA) nat- where an attractive line of fixed points exists (Fig. 9) 4 urally involves these infrared scales. The separation be- with α∗ = π ∆, calculations show that the conductivity tween the Dirac liquid and the Fermi liquid regimes in (on the fixed line) increases relatively to the free limit the n T plane is defined by the crossover tempera- (Herbut et al., 2008): σ = [π/2 +(4 π)∆]e2/h. For − − ture T ∗(n) = vkF (1 + (α/4) ln(Λ/kF )), kF = √πn, and stronger scalar and vector disorder/interactions where thus the temperature dependencies quoted previously, the couplings run away to infinity the problem is non- are valid for T0 >T >T ∗(n). The ultraviolet tem- perturbative, and a complex variety of behavior is ex- 4 perature scale T0 8 10 K, while for typical graphene pected (Foster and Aleiner, 2008). 12≈ ×2 2 densities n . 10 cm− , T ∗(n) 10 K. For clean graphene at µ = vkF = 0 it was pointed out ∼ (Fritz et al., 2008; Kashuba, 2008; M¨uller et al., 2008) that at high temperature (compared to the frequency), the conductivity is expected to have the form:
2. Conductivity 2 0.76 e 2 σ = 2 , T α ω, (3.50) The behavior of the electrical conductivity in graphene α h ≫ has been extensively reviewed (Das Sarma et al., 2011; where α(T )=4/ ln(Λv/T ) is the running Coulomb cou- Peres, 2010). It is believed that charged impurities pling. This form reflects electron-electron inelastic colli- 2 and resonant scatterers are the main sources of scat- sions with scattering rate 1/τee α T . The linear tem- tering away from the Dirac point, and to extent the perature dependence is characteristic∼ for Dirac particles. long- or short-range part of the Coulomb potential con- The above formula is valid as long as 1/τee is the dom- tributes to scattering is a matter of ongoing debate inant scattering mechanism (collision-dominated trans- (Chen et al., 2008; Monteverde et al., 2010; Ni et al., port), and implies that clean graphene at the neutral- 2010; Ponomarenko et al., 2009; Reed et al., 2010). ity point should exhibit a universal, interaction-limited Here we will only mention effects related to long-range conductivity, reflecting essentially the quantum critical electron-electron interactions near the Dirac point. Inter- behavior of graphene in this regime (T µ). With ≫ action corrections to the minimum metallic conductivity increased doping (µ/T ), a crossover takes place to a 2 π 2 of free Dirac fermions, σ0 = e /(4~) = e /h (Fradkin, Fermi liquid regime with screened interactions, where 2 1 2 2 1986; Lee, 1993), are more involved, because this expres- τ − α T /µ, (M¨uller et al., 2008) and the conductiv- ee ∼ sion does not contain the quasiparticle velocity, while the ity is dominated by charged impurity scattering. electric charge is not renormalized. The debate was fu- It has also been pointed out that for µ = 0 graphene eled in part by electrical measurements of the minimum behaves as an almost “perfect” fluid, in a sense that conductivity (at the Dirac point) which turned out to its shear viscosity, η, relative to the entropy density 2 ~ be somewhat larger than σ0 (Geim and Novoselov, 2007; s is anomalously small: η/s = (0.13/α (T ))( /kB) Tan et al., 2007). Theoretically, at T = 0 (or T ω (M¨uller et al., 2009). This ratio measures how strongly where ω is the external frequency), it is expected≪ that the excitations in a fluid interact. At room temperature any interaction effect should have sub-leading character, η/s of graphene is smaller than η/s of any known cor- and the frequency can enter only through the running related quantum fluid, and is close to the lower bound of 1 ~ proposed to exist for a large class of strongly of the coupling α(ω). Even though some debate still ex- 4π kB ists (Herbut et al., 2008; Juricic et al., 2010; Mishchenko, interacting quantum field theories (Kovtun et al., 2005). 2008; Sheehy and Schmalian, 2009) as to the implemen- Therefore, due to its quantum critical nature near the tation of the cut-off regularization procedure, the con- Dirac point, graphene is suggested to behave as a strongly 22 correlated quantum liquid and should exhibit signatures found to be amplified by excitonic effects, improv- of electronic turbulence (M¨uller et al., 2009). ing screening of interactions between quasiparti- cles. This analysis leads to values of α ranging from α 1/7 in the static limit to α 2 at high E. Overview of main results frequencies.≈ Very recent measurements≈ of the cy- clotron mass in suspended graphene (Elias et al., Before we proceed with further topics related to inter- 2011) have found logarithmic velocity renormaliza- actions in graphene, let us broadly summarize the main tion and extract, within the RPA scheme, an ef- findings and questions raised so far: fective value of graphene’s dielectric constant ǫG 3.5. One can also expect that near the Dirac point,≈ 1. For clean graphene at the neutrality point µ = 0, where interactions lead to singular effects, addi- interactions are not screened and are marginally tional factors can be important such as disorder, irrelevant; the fixed point α = 0 is approached ∗ inhomogeneities, rippling, etc., and thus obscure logarithmically (or, equivalently, the quasiparticle the clean behavior. velocity increases logarithmically). From a theory standpoint, the approach towards this fixed point 5. In the Fermi-liquid regime, where interactions are is well understood both from weak and strong- screened, the physics near the Dirac point can still coupling (RPA) perspectives. Since in graphene be strongly affected — this is due to resonant fea- one can have α 1 under rather conventional ex- ∼ tures in the quasiparticle self-energy, reflecting in- perimental conditions, our understanding of RPA teractions of quasiparticles with plasmons. calculations is important. RPA is justified only in the limit of large number of fermion species (N 1), while for N = 4 it should work for IV. THE COULOMB PROBLEM AND CHARGED weak≫ to moderate coupling; however there are indi- IMPURITIES cations, coming mostly from two-loop calculations, that vertex corrections are numerically small, and The consideration of non-interacting Dirac electrons thus RPA should work well. Disorder generally in 2D under a Coulomb field is of paramount relevance drives the system away from the clean fixed point, for graphene, and for several reasons. First of all, the towards finite or even strong coupling, depending Coulomb problem for relativistic fermions has many fea- on disorder type. tures that are unfamiliar in condensed matter systems, and which resemble long standing predictions made in the 2. The resulting behavior near the Dirac point is that context of QED in strong fields. As such, and given that of a non-Fermi-liquid with a quasiparticle decay having α 1 makes graphene intrinsically strongly cou- rate which is linear in energy, and decreasing quasi- pled, it can∼ provide the first experimental ground for test- particle residue. All physical characteristics related ing many elusive predictions from strong-coupling QED. to the quasiparticle velocity (which increases loga- On the other hand, the single particle Coulomb prob- rithmically) are affected, and predicted to exhibit lem constitutes the first step in addressing nontrivial fea- systemic, interaction dependent, deviations from tures of the full, many-body interacting problem. Char- their non-interacting values as the Dirac point is acteristics like non-linear screening, or the supercritical approached, either as a function of density or tem- instabilities, provide valuable insight in grasping some perature. proposed many-body effects, like exciton condensation, 3. Can graphene be driven into an excitonic insulating or spontaneous mass generation in graphene. state? At the Dirac point the long-range Coulomb Historically, however, the motivation for studying the interactions can lead to bound electron-hole pairs, Coulomb problem comes from the seminal experimental creating a gap. There has been intense debate observations (Novoselov et al., 2004a) that the field effect whether this can happen under realistic conditions in graphene prepared on SiO2 is characterized by carrier — since the critical interaction strength appears to mobilities that do not depend on the Fermi energy or carrier density (the DC conductivity, σ = me n , with be αc 1, it seems possible to occur in suspended ∼ m const.), and that carriers are chiral Dirac| fermions| samples (α = 2.2). So far no experimental indica- ≃ tions have been observed. in 2D (Novoselov et al., 2005; Zhang et al., 2005). Early semiclassical investigations (Adam et al., 2007; Ando, 4. What is the value of the interaction α? Clearly, 2006; Nomura and MacDonald, 2007, 2006) showed that since α = 2.2/ǫ0 is dielectric constant depen- such linear-in-density conductivity could be explained by dent, working with different substrates could scattering of unscreened Coulomb impurities, which are 10 2 lead to changes in interaction-dependent effects typically seen in silica in concentrations of 10 cm− (Jang et al., 2008). There are also suggestions that (Ando et al., 1982). As a result, transport∼ in the pres- graphene has an “intrinsic” value of α (Reed et al., ence of charged impurities rapidly became one of the 2010), arising from dynamical dielectric screen- most studied topics in the quest for the ultimate mobil- ing. The polarizability of the Dirac fermions was ity in graphene. Since, as we saw before, Coulomb’s law 23 is exactly preserved in undoped graphene, and approx- (Khalilov and Ho, 1998; Novikov, 2007a) imately preserved for small and moderate doping, the scattering processes are essentially governed by the bare m ε g/r Fj (r) + [∂r + j/r]Gj (r) = 0 (4.2a) − − Coulomb problem, unlike conventional metals, where ∂r j/rFj (r)+ m + ε + g/r Gj (r)=0. (4.2b) screening is perfect. A thorough understanding of this − problem is therefore important not only for its theoreti- This coupled pair of first order equations can be straight- cal relevance and its import on electron-electron interac- forwardly reduced to two decoupled second order equa- tions, but also for its experimental implications, and our tions. Free solutions (g =0)of (4.1) exist when ε > m , understanding of transport in graphene. and are simple spherical waves whose k-normalized| | | ver-| Finally, it is highly significant that this is an exactly sion reads solvable problem. This means that most quantities can be obtained exactly, allowing us to unveil many interact- k ε + m Jj 1/2(kr)Φj 1/2 Ψj = | | − − (4.3) ing and non-interacting effects that are not within reach s2 ε isε ε m Jj+1/2(kr)Φj+1/2 of the perturbative approaches already discussed. We | | p | − | p proceed to show several such features. On account of (sx sgn(x)). For nonzero g, one readily sees from (4.2) the long range nature of the Coulomb field, inter-valley that≡ the solutions at r 0 behave as processes are not relevant, and hence we will solve the ∼ γ 2 2 problem within each (independent) valley in the Dirac F (r), G(r) r± , γ = j g . (4.4) description of fermions in graphene. ∼ − The general exact solution is given inp terms of confluent hypergeometric, or Whittaker’s functions, both in the massive (Gamayun et al., 2009; Gupta and Sen, 2008; A. Exact Solution of the Coulomb Problem Gupta et al., 2010; Khalilov and Ho, 1998; Novikov, 2007a; Pereira et al., 2008a), and massless cases 1. Wave Equations and Spectrum (Gupta and Sen, 2009; Pereira et al., 2007; Shytov et al., 2007b). In the massless case, one can map (4.2) into A Coulomb center of charge Z e generates the poten- the familiar Coulomb radial Schrodinger equation in 3D 2 | | tial U(r) = Ze /(ǫ0r) for the electrons. Without any (Pereira et al., 2007): loss of generality let us consider Z > 0. The electronic 2 2 2 dynamics is governed by the wave equation ∂r f + ε +2gε/r γ(γ 1)/r f (r)=0, (4.5) ± − ∓ ± g where the f are linear combinations of F and G, ε2 v iσ ∇ + σ Mv Ψ(r)= EΨ(r). (4.1) ± − · − r 3 takes the place of the Schrodinger energy, and γ plays the role of angular momentum. Since the solution is for- 2 2 mally the same, the appearance of ε instead of ε means Here we use g = Zα = Ze /(ǫ0v), with ǫ0 reflect- ing the effective dielectric constant of the embedding that the massless case admits no bound solutions, as we medium, and the mass M accounts for the more gen- expect on account of the absence of a spectral (mass) eral possibility of a symmetry breaking gap. Through- gap. The massive case, however, has a well defined infi- nite spectrum of bound solutions when ε < m , given out this chapter we shall use the scaled energy and | | | | mass ε = E/v, m = Mv, and k = √ε2 m2. by (Khalilov and Ho, 1998) Even though m = 0 for ideal graphene without− in- n + j2 g2 teractions, nonzero m can be induced in many ways. εn,j = sgm − , (4.6) 2 One of them is through interaction with suitable sub- g2 + np+ j2 g2 strates, of which some experimental hints have been − q reported (Gr¨uneis and Vyalikh, 2008; Li et al., 2009a; p lowest level is given by ε ε = s m 1 (2g)2. Martinazzo et al., 2010; Zhou et al., 2007). In terms of G ≡ 0,1/2 g − the original tight-binding Hamiltonian, the mass M aris- p ing from a sublattice symmetry is related to the parame- 2 2. Supercritical Instabilities ter ∆0 introduced in eq. (3.40) via Mv = ∆0. The axial symmetry of the potential allows us to use the eigen- Consideration of eq. (4.4) immediately reveals a com- states of the total pseudo angular momentum, Jz = Lz + plication if g>gc = 1/2, because γ becomes imaginary σz/2, which is conserved (DiVincenzo and Mele, 1984). for the lowest angular momentum channels (j = 1/2). 1/2 ± We write Ψj† = r− [Fj (r)Φj 1/2(φ), iGj (r)Φj+1/2(φ)], The solution (4.4) is neither regular nor divergent, but − where j = 1/2, 3/2,... are the eigenvalues of Jz, rather oscillates endlessly towards r = 0. This is patho- ± ± ipφ and the cylindrical harmonics read Φp(φ) = e /√2π. logical because the space of solutions is of dimension 2, A detailed derivation of the 2D Dirac equation for gen- and we can no longer discard an irregular contribution eral radial potentials is given by Novikov, 2007a. In our since both linearly independent solutions are square in- case, Eq. (4.1) reduces to the following radial equations tegrable. In other words, there is no boundary condition 24
the pure Coulomb case, εG(g) decreases towards zero in a singular way at g = gc. In a regularized poten- tial, εG depends also on the cutoff radius R, and is al- lowed to monotonically penetrate the negative energy re- gion, until eventually touching the lower continuum at ε = m. If g is further increased, εG dives into the hole (positron)− continuum and becomes a resonance. Other levels will sequentially follow at higher g. The diving point for εG(g) defines a renormalized critical coupling, g˜c > gc that is characterized by a log singularity at 2 2 mR 0:g ˜c gc + π / log (mR) (Gamayun et al., 2009; Khalilov∼ and≃ Ho, 1998; Pereira et al., 2008a; Zhu et al., FIG. 18 (Color online) Schematic drawing of the level diving 2009), strongly depending on the regularization. process in the supercritical regime, and of the resulting quasi- This diving of bound levels entails a complete re- spectrum of levels for massive and massless fermions. structuring of the vacuum. If the level was empty, an electron-hole pair will be immediately created: the elec- tron remains tightly bound and shielding the center, at the origin to univocally select the solution. Secondly, while the hole is ejected to infinity (Greiner et al., 1985; in the massive case the level ε becomes imaginary, sig- G Zeldovich and Popov, 1972). The supercritical regime naling a loss of self-adjointness of the Dirac Hamiltonian is thus characterized by spontaneous pair creation, or for g > 1/2. a spontaneous Schwinger mechanism (Schwinger, 1951). Physically, both effects are a symptom that the po- One expected consequence is a strong signature of these tential has such a strong divergence that particles are resonances in the hole sector of the scattering and trans- inexorably attracted and “fall” into the origin, leading port cross sections. to a collapse of the system (for example, the endless os- An essential detail is that these resonances are not cillations can be read as an infinite phase shift). This usual bound levels diluted inside a continuum, where “fall to the center” is a general characteristic of diverg- their lifetime essentially disappears. One consequence ing potentials in any dimension of space. For power law of the chiral nature of Dirac fermions, combined with potentials, one particular power signals the threshold of the long range tail of the Coulomb potential, is that the criticality. The Coulomb potential is the marginal case supercritical levels in the relativistic Coulomb remain for the Dirac equation (both in 2D and 3D), just like the sharply defined, with diverging lifetime. For example, potential 1/r2 is the marginal case of the 3D Schrodinger for S states (j = 1/2), one shows that these resonances equation (Landau and Lifshitz, 1981). This, of course, follow (Gamayun et al., 2009) begs the question of regularization. Regularizing the po- tential introduces an additional boundary condition at 3π π/√2ξ 3π(β βc) some short distance R, which allows a formal solution, εn m 1+ξ+i e− , ξ = − , (4.7) ≈− 8 8ββc and cures the total collapse of the system (Case, 1960; 2 Perelomov and Popov, 1970). In graphene the lattice is when g & g˜c, and where β = iγ, βc = g˜c 1/4. In the natural regulator and there are no ultraviolet issues. real space the localization of the supercritical− levels is p But the physics in the supercritical regime depends ex- controlled by the reduced Compton wavelength: λC = plicitly on the short range details. 1/(mv). The modulus squared of their wavefunction de- This supercritical collapse has a long history in the cays as Ψ†Ψ exp( 8gr/λC ) and, consequently, even context of QED, where the Dirac equation stands as the inside the continuum,∝ − such levels retain a highly localized basis for understanding the stability of matter. In QED nature, which is whyp they are so relevant, in particular the collapse would occur for ZαQED > 1, which lead to in their potential for screening (Pereira et al., 2008a). extensive investigations regarding the stability of heavy nuclei having Z>Zc = 137 (Case, 1960; Greiner et al., 1985; Popov, 1971a,b; Zeldovich and Popov, 1972). Af- b. Massless Electrons. The spectrum in this case is con- ter regularization Zc 170, which makes the problem tinuous everywhere, and thus there is no sequential div- highly academic, and→ QED’s predictions untestable. In ing and restructuring of the hole continuum as described graphene, on the contrary, Zc 1, which opens the real above. But the pathology associated with Eq. (4.4) still possibility of testing the supercritical∼ instability in a con- exists. Physically, the massless situation is rather more densed matter setting. catastrophic since the solution in a regularized potential reveals an infinite number of quasi-localized resonances in the hole sector (Gamayun et al., 2009; Pereira et al., a. Massive Electrons. To understand the physics in the 2007; Shytov et al., 2007b). This is a highly non-trivial supercritical regime we can follow the level εG as effect for several reasons: (i) in the massless case there the coupling increases (Fig. 18) (Greiner et al., 1985; is no natural length scale in the problem to characterize Pereira et al., 2008a; Zeldovich and Popov, 1972). For such localized states; (ii) the system abruptly develops 25 an infinite quasi-bound spectrum at g >gc, when its spectral fingerprint is rather featureless for g a + ib πn/√g2 g2 ε e− − c , (a,b) (g), (4.8) n ≈− R ∼ O which has an essential singularity at gc, an energy scale/lower bound set explicitly by the regularization dis- tance, R, and diverging lifetimes close to the critical point. Since the width of these states vanishes linearly, they are practically bound states (hence the designation quasi-bound states). In real space, the localization scale is determined by the regularization distance R itself. Since meso and nanoscopic devices are of high interest, it is pertinent pointing out that massless Dirac fermions FIG. 19 (Color online) (a) LDOS, N(ε,r) at r = a for several in a finite-sized system mimic in all aspects the physics couplings g 2 lation at ε = 0 is confirmed in Fig. 19(d) where we show Pj 2gsε 2 2 the supercritical contributionn ¯j (ε, r) as a function of nj(ε, r)= Fγ 1 + Fγ + Fγ Fγ 1 (4.9) 2π2γ2r − j − log( ε ). At positive energies the LDOS exhibits peri- | | odically| | decaying oscillations in εr [inset of Fig. 19(c)], for g < g , and F represents the Coulomb func- c l with extrema separated by nπ, within logarithmic ac- tion F ( gs , ε r) (Abramowitz and Stegun, 1964). The l ε curacy (Shytov et al., 2007a≈). When directly measured function−N(ε,| r)| is plotted in Fig. 19(a) for different cou- in STM such oscillations can be used to extract the elec- plings and distances. Apart from the evident particle- tronic dispersion, as done by Ouyang et al., 2002. hole asymmetry, the LDOS remains rather featureless, We point out that, since the solution of the supercrit- even at the shortest distances. If g >gc the corre- sponding analytical expression obtained in the regular- ical problem involves a nontrivial ad-hoc regularization, ized potential is more complex, but still has a closed form these results have been checked numerically against ex- (Pereira et al., 2007). In this case, supercritical channels act solution of the full tight-binding problem in the hon- ( j < 1/2) need to be isolated from undercritical ones eycomb lattice, being found that the analytical Dirac re- (|j| > 1/2), yielding two contributions to the LDOS: sults reproduce the full lattice problem down to distances | | as small as the lattice scale (Pereira et al., 2007). N(ε, r)= n¯j (ε, r)+ nj (ε, r). (4.10) The striking differences between the two regimes and j < g j > g the violent modification of the ground state at strong cou- | X| | | | X| | | pling are likewise evident in the behavior of the scattering The total LDOS for this case is shown in Fig. 19(c) for phase-shits, δj (ε). They admit closed formed expressions g = 1.0, and at different distances to the impurity. It is at both g Shytov et al., 2007a) and g > gc (Castro Neto et al., recognized that conventional procedures of the theory of 2009b; Shytov et al., 2007b). For example, the under- metals, like self-consistent screening, linear response or critical S-matrix reads (Novikov, 2007a) Friedel sum rules, are not straightforward in this sys- tem. For example, within the Dirac (effective mass) ap- iπ(j γ) 2iδj (ε) je − Γ(1 + γ igsε) proximation, the ultraviolet cutoff scale enters explic- Sj (ε)= e = − , (4.11) γ igsε Γ(1 + γ + igsε) itly in Friedel’s sum rule, and Levinson’s theorem is − modified (Lin, 2006) (Levinson’s theorem is one of the which is energy independent, but considerably asymmet- fundamental results in quantum scattering theory, as- ric with respect to the sign of g. The corresponding δj serting that in the Schr¨odinger’s equation with a non- are shown in Fig. 19(b)(inset) as a function of coupling singular spherically symmetric potential the zero en- strength. Note how δ1/2 (the most important partial ergy scattering phase-shift exactly counts the number of wave) behaves rather differently from the others: only bound states: δl(0) = Nlπ). One consequence is that δ1/2 shows the expected sign for the attractive/repulsive a na¨ıve application of Friedel’s sum rule can yield di- situations. On the other hand, in the supercritical regime vergent displaced charges (DiVincenzo and Mele, 1984). there is a strong ε-dependence of δj . In the top row of Even though these divergences are artificial in the target Fig. 19(d) we present (δj mod π) as a function of log(ε). lattice problem, they point, already at a single particle In the attractive sector (ε< 0 if g > 0) the abrupt steps level, to the anomalous screening properties of graphene. centered around π/2 mark the position of the infinite quasi-bound spectrum (which, as per (4.8), accumulates exponentially at ε = 0), whereas in the attractive sector 1. Weak Coupling (g