# Electron-electron interactions and doping dependence of the two-phonon Raman intensity in graphene

###### Abstract

Raman spectroscopy is a fast, non-destructive means to characterize graphene samples. In particular, the Raman spectra are strongly affected by doping. While the change in position and width of the G peak can be explained by the non-adiabatic Kohn anomaly at , the significant doping dependence of the 2D peak intensity has not been explained yet. Here we show that this is due to a combination of electron-phonon and electron-electron scattering. Under full resonance, the photogenerated electron-hole pairs can scatter not just with phonons, but also with doping-induced electrons or holes, and this changes the intensity. We explain the doping dependence and show how it can be used to determine the corresponding electron-phonon coupling. This is higher than predicted by density-functional theory, as a consequence of renormalization by Coulomb interactions.

## I Introduction

Graphene is the latest carbon allotrope to be discovered, and it is now at the center of a significant research effortNov306(2004) (); GeimRevNM (); Nov438(2005) (); CastroNetoRev (); charlier (); Zhang438(2005) (). Near-ballistic transport at room temperature and high mobilityNov438(2005) (); Zhang438(2005) (); Nov315(2007) (); MorozovNov(2007) (); andrei (); kimmob () make it a potential material for nanoelectronicsHan (); Chen (); Zhang86 (); Lemme (), especially for high frequency applicationsYuming (). Furthermore, its transparency and mechanical properties are ideal for micro and nanomechanical systems, thin-film transistors and transparent and conductive composites and electrodesbunch (); blake1 (); hernandez (); eda ().

Graphene layers can be readily identified in terms of number and orientation by inelastic and elastic light scattering, such as RamanACFRaman () and Rayleigh spectroscopiesCasiraghiNL (); GeimAPL (). Raman spectroscopy also allows monitoring of doping, defects, strain, disorder, chemical modifications and edgesACFRaman (); CasiraghiAPL (); Pisana (); ACFRamanSSC (); DasCM (); Mohi (); ferralis (); YanPrl2007 (); Leandro (); PiscanecPRL (); cancado08 (); Cedge (); Elias (); Dasbila (); Ferrari00 (); Graf2007 (). Indeed, Raman spectroscopy is a fast and non-destructive characterization method for carbonsacftrans (). They show common features in the 800-2000 cm region: the G and D peaks, around 1580 and 1350 cm, respectively. The G peak corresponds to the phonon at the Brillouin zone center ( point). The D peak is due to the breathing modes of six-atom rings and requires a defect for its activationtuinstra (); Ferrari00 (); ThomsenPrl2000 (). It comes from TO phonons around the K point of the Brillouin zonetuinstra (); Ferrari00 (), is active by double resonance (DR)ThomsenPrl2000 (), and is strongly dispersive with excitation energy due to a Kohn Anomaly at KPiscanecPRL (). The activation process for the D peak is inter-valley, and is shown schematically in Fig. 1(d): i) a laser induced excitation of an electron/hole pair; ii) electron-phonon scattering with an exchanged momentum ; iii) defect scattering; iv) electron-hole recombination. DR can also happen as intra-valley process, i. e. connecting two points belonging to the same cone around K (or ), as shown in Fig. 1(b). This gives the so-called D’peak, which is at in defected graphite measured at 514nm.

The 2D peak is the second order of the D peak. This is a single peak in single layer graphene (SLG), whereas it splits in four in bilayer graphene (BLG), reflecting the evolution of the band structureACFRaman (). The 2D’ peak is the second order of the D’ peak. Since both 2D and 2D’ originate from a process where momentum conservation is satisfied by two phonons with opposite wavevectors (q and ), they do not require the presence of defects for their activation, and are thus always present. Indeed, high quality graphene shows the G, 2D and 2D’ peaks, but not D and D’ACFRaman (). Also, under the assumption of electron-hole symmetry, the two-phonon peaks are fully resonantBasko2007 (); BaskoBig (). This means that energy and momentum conservation are satisfied in all elementary steps of the Raman process, as shown schematically in Fig. 1(c,e). Then, all intermediate electronic states are real. As a consequence, two-phonon Raman spectroscopy is sensitive to the dynamics of the photo-excited electron-hole pair, in particular, to the scattering processes it can undergo. This is of crucial importance for the present work.

The effects of doping on the graphene G-peak position [Pos(G)] and Full Width at Half Maximum [FWHM(G)] were reported in Refs. Pisana, ; YanPrl2007, ; DasCM, ; Dasbila, . Pos(G) increases and FWHM(G) decreases for both electron and hole doping. The G peak stiffening is due to the non-adiabatic removal of the Kohn-anomaly at Lazzeri2006 (); Pisana (). The FWHM(G) sharpening is due to Pauli blocking of phonon decay into electron-hole pairs, when the electron-hole gap is higher than the phonon energyLazPRB2006 (); Pisana (), and saturates for a Fermi shift bigger than half phonon energyPisana (); YanPrl2007 (); LazPRB2006 (). A similar behavior is observed for the LO-G peak in metallic nanotubesdasnt (), for the same reasons. In the case of BLG, the different band structure re-normalizes the phonon response to doping differently from SLGYan2008 (); Dasbila (); andobi (). Also in this case the Raman G peak stiffens and sharpens for both electron and hole doping, as a result of the non-adiabatic Kohn anomaly at Dasbila (). However, since BLG has two conduction and valence subbands, with splitting dependent on the interlayer coupling, this changes the slope in the variation of Pos(G) with doping, allowing a direct measurement of the interlayer coupling strengthDasbila (); andobi ().

Another significant result is that in SLG the ratio of the heights of the 2D and G peaks, I(2D)/I(G), and their areas, A(2D)/A(G), is maximum for zero dopingACFRaman (); heinz (); berciaud (), and decreases for increasing doping. On the other hand, this shows little dependence on doping for BLGDasCM (); Dasbila (). Fig. 2 plots the combined data for SLG and BLG from Refs. DasCM, ; Dasbila, ; ACFRaman, ; heinz, ; berciaud, ; tan1, . Note that Refs. DasCM, ; Dasbila, reported height ratios, while here, as discussed later, we analyze the area ratio A(2D)/A(G), which encompasses both trends of I(2D)/I(G) and FWHM(2D)/FWHM(G).

Due to residual disorder, the energy of the Dirac point can fluctuate across the sample on a scale smaller than the laser spot, which leads to spatial inhomogeneity of the doping levelCasiraghiAPL (); yacoby (). We attribute the difference in the behavior of the two SLG curves in Fig.2 to a different degree of residual charge inhomogeneity in the polymeric electrolyte experiments of Refs. DasCM, ; Dasbila, . On the other hand, the use of this electrolyte enabled probing a very large doping range, because the nanometer-thick Debye layer gives a much higher gate capacitance compared to the usual 300nm SiO back gate DasCM (); Dasbila (); Pisana (). Note as well that A(2D)/A(G) for the most intrinsic samples measured to date is 12–17 ACFRaman (); heinz (); berciaud (); tan1 (), much higher than the zero gating values in Refs. DasCM, ; Dasbila, , as shown in Fig. 2. This points again to sources of disorder in the gated samples of Refs. DasCM, ; Dasbila, , while the absence of a significant D peak excludes large amounts of structural defects.

Here, we show that the 2D intensity doping dependence results from its sensitivity to the scattering of the photoexcited electron and hole. Assuming the dominant sources of scattering to be phonon emission and electron-electron collisions, we note that, while the former is not sensitive to doping, the latter is. Then, the 2D doping dependence can be used to estimate the corresponding electron-phonon coupling (EPC).

## Ii Doping Dependence of Two Phonon Raman Intensity

### ii.1 Theoretical Dependence

Raman scatteringRaman () is an electron-mediated process where electromagnetic radiation exchanges vibrational quanta (phonons) with a crystal. A complete description requires the detailed knowledge of (i) electronic structure, (ii) phonon dispersions, (iii) mutual interactions between electrons and phonons (i.e. electron-electron, electron-phonon and phonon-phonon scattering).

The Raman spectrum of graphene consists of a set of distinct peaks. Each characterized by its position width, height, and area. The frequency-integrated area under each peak represents the probability of the whole process. It is more robust with respect to various perturbations of the phonon states than width and height. Indeed, for an ideal case of dispersionless undamped phonons with frequency the shape of the -phonon peak is a Dirac distribution , with zero width, infinite height, but well-defined area. If the phonons decay (e. g, into other phonons, due to anharmonicity, or into electron-hole pairs, due to electron-phonon coupling), the lineshape broadens into a Lorentzian, but the area is preserved, as the total number of phonon states cannot be changed by such perturbations. If phonons have a weak dispersion, states with different momenta contribute at slightly different frequencies. This may result in an overall shift and a non-trivial peak shape, but frequency integration across the peak means counting all phonon states, as in the dispersionless case. Thus, the peak area is preserved, as long as the Raman matrix element itself is not changed significantly by the perturbation. The latter holds when the perturbation (phonon broadening or dispersion) is smaller than the typical energy scale determining the matrix element. Converting this into a time scale using the uncertainty principle we have that, if the Raman process is faster than the phonon decay, the total number of photons emitted within a given peak (i. e., integrated over frequency across the peak), is not affected by phonon decay, although their spectral distribution can be. Although the graphene phonons giving rise to the D and D’ peak are dispersive due to the Kohn Anomalies at K and PiscanecPRL (), their relative change with respect to the average phonon energy is at most a few %, thus we are in the weakly dispersive case discussed above. The phonon decay in graphene is in the picosecond timescale, while the Raman process is faster, in the femtosecond timescalePisana (); Bonini2007 (); LazzeriHP (). Then, we will analyze the area ratio, A(2D)/A(G), which encompasses both variations in height ratio, I(2D)/I(G), and width: FWHM(2D)/FWHM(G).

We first consider the G peak. For the one-phonon process, allowed by momentum conservation, which gives rise to the G peak, the picture is entirely different from the two-phonon case. As shown in Fig.1a, the process responsible for the G peak is determined by virtual electron-hole pairs with energy , where is the laser excitation energy (for a typical visible Raman measurement ). If the Fermi energy, , stays below , as in Refs. DasCM, ; Dasbila, , these electronic states are not strongly affected. Only the final phonon state is influenced by doping, which manifests itself in a change of Pos(G) and FWHM(G)DasCM (); Dasbila (); Pisana (); YanPrl2007 (). However, the area of the peak is determined by the total spectral weight of the phonon state, which is preserved. Thus, we do not expect any significant dependence of A(G) on doping, as long as the doping is not too strong, so that . We can then take the measured doping dependence of A(2D)/A(G) as representative of the A(2D) trend. Note that A(G) can change as a function of other external parameters, such as the Raman excitation energyVidano (); Pocsik1998 (); Cancado2007 (); ACFRaman (); Ferrari00 (). However, for fixed excitation, such as in the experiments discussed here, the above argument holds.

In Ref.BaskoBig, the following expressions for the 2D and 2D’ areas were obtained:

(1a) | |||

(1b) |

where is the electron charge, is the speed of light, is the fine structure constant, and is the electron velocity (its experimental value is Jiang2007 (); Rotenberg2007SSC (); Zhou2008 ()). is the scattering rate of the photoexcited electron and hole. Note that we define as the imaginary part of the energy, so it determines the decay of the amplitude, while the decay of the probability is determined by . This includes all sources of inelastic scattering. Assuming the two main mechanisms for electron scattering to be the emission of phonons and electron-electron collisions, we write:

(2) |

Here we include the phonons near and K, responsible for D and D’. The corresponding emission rates, and , enter the numerators in Eqs. (1a), (1b).

Two points regarding Eqs. (1a), (1b) should be emphasized. First, the scattering rates depend on the electron energy, , which is defined by half the laser energy, [see Eq. (III.1) in the next section]. Second, if impurity scattering is significant compared to other scattering mechanisms, the corresponding elastic scattering rate cannot be simply included in and Eqs. (1a), (1b). The whole Raman intensity calculation should be done differently. Eqs. (1a), (1b) thus neglect impurity scattering. For short-range impurities this assumption is justified by the absence of a large D peak in the spectra of Refs. DasCM, ; Dasbila, . Long-range disorder is efficiently screened (even though the vanishing density of states at the Dirac point requires the screening to be nonlinearFogler (); Novikov (); DasSarma (); Polini ()); it is precisely this screening that gives rise to the inhomogeneous concentration of electrons/holes and spatial fluctuations of the Dirac point energy.

In principle, there are no reasons for a strong dependence of on carrier density. However, does exhibit such a dependence. Indeed, in undoped graphene at low temperatures, the photoexcited electron finds itself in a state with some momentum, , measured from the Dirac point, in the empty conduction band. To scatter into a state with a different momentum , it has to give away some energy and momentum to another electron in the full valence band. This second electron would have to be promoted to the conduction band (as there are no available empty states in the valence band) into a state with momentum , leaving a hole in the valence band with . Momentum and energy conservation require:

(3a) | |||

(3b) |

where is the quasiparticle dispersion, assumed the same for electrons and holes. For Dirac particles, , the only possibility to satisfy both conservation laws is to have all four momenta parallel. If the spectrum is convex, , the two equations can be satisfied by a set of momenta with non-zero measure, i. e. the phase space is finite. If it is concave, , they are incompatible. In SLG the spectrum is Dirac to a first approximation, resulting in an uncertaintyGuinea1996 (). This can be resolved by taking into account corrections from electron-electron interactions, which make the spectrum concave,AbrikosovBeneslavskii (); Guinea94 () and the interband process forbidden.

As new carriers are added to the system, intraband electron-electron collisions become allowed. The momentum and energy conservation become:

(4) | |||

(5) |

which can be satisfied for any quasiparticle dispersion. These collisions give a contribution to which increases with carrier concentration. As a consequence, the total in Eq.1a increases, leading to an overall decrease of A(2D), consistent with the experimental trend in Fig. 2.

The above arguments essentially use the non-convexity of the electronic spectrum in the conduction band, and thus apply to SLG only. In BLG, the spectrum is parabolic near the Dirac point, so that , and the phase-space restrictions are absent. Thus, electron-electron collisions are allowed even at zero doping, and the collision rate has a much weaker dependence on , which, in first approximation, can be neglected. Thus, A(2D) is expected to have a weak dependence on , as seen in Fig. 2, where the experimental A(2D)/A(G) for BLG shows a negligible variation with dopingDasbila ().

To quantify the doping effects on the SLG A(2D), we first calculate the electron-electron scattering rate, , in the random-phase approximation, analogously to Refs. DasSarma2008, ; Polini2008, . is given by the imaginary part of the on-shell electronic self-energy, for , with and p counted from the Dirac pointGuinea1996 (). Here we consider the limiting case, when the energy of the photoexcited electron () far exceeds . The carrier concentration is . In this case, the collisions are dominated by small momentum transfers, , so does not depend on and is proportional to , the proportionality coefficient depending only on the dimensionless Coulomb coupling constant ( being the dielectric constant):

(6) |

where the function is given by:

and , , are:

(8a) | |||

(8b) | |||

(8c) | |||

(8d) |

Fig. 3 plots , calculated numerically.

Thus, we expect A(2D) to change with as:

(9) |

with a constant. Note that a variation of the dielectric constant will affect A(2D). Given the negligible dependence of A(G) on doping, Eq. (9) can be rewritten as

(10) |

where is another constant.

### ii.2 Fit to Experiments

Fig. 4 plots as a function of . This dependence, according to Eq. (10), should correspond to two symmetric straight lines joining at . As noted in Sec. I, close to the data from the two polymer electrolyte gating experiments do not converge to the same value. However, for both a linear rise of is seen at higher energies. Also, while the data represented by open circles in Fig.4 are almost symmetric, a significant asymmetry is seen for electron doping in the set represented by the half-filled circles, while the two sets are in good agreement for hole doping.

A(2D)/A(G) for intrinsic samples measured at 514.5 nm excitation, the same used in Refs.DasCM, ; Dasbila, , is in the range 12-17ACFRaman (); heinz (); berciaud (), represented by the star in Fig. 4 at 14.5. This is in good agreement with the ratio measured for carbon whiskerstan1 (). These show a 2D peak very similar to graphene, being composed of mis-oriented graphene layerstan1 (); latil (). However, their Raman spectra are much less susceptible to charged impurities or surface doping, being bulk materialstan1 (). This corresponds to , which we use to eliminate the effect of doping inhomogeneity, by constraining at zero doping. We also need to consider the dielectric constant of the polymer electrolyte DasCM (), , giving . Thus, we fit the data with a one-parameter expression:

(11) |

We fit separately each branch of the two data-sets, as shown by solid and dotted lines in Fig. 4. We get : 18, 21, 29, 65 meV, with an average meV.

## Iii Raman Intensities and Electron-Phonon Coupling

### iii.1 Theoretical Background and Electron-Phonon Coupling Definitions

Even though graphite and other -hybridized materials have been investigated for more than 50 yearsWallace1947 (); tuinstra (), all the fundamental physical properties needed for the interpretation of the Raman spectra have undergone an intense debate, which seems to be just beginning to converge. Interestingly, several features of both phonon dispersions and band structure of graphene are determined by the EPC. For example, in the Kohn anomalies around or KPiscanecPRL () the correction to the phonon frequencies due to EPC results in a linear slope of the optical phonon branches as the wave vector approaches or K. The EPC and phonon dispersions calculations of Ref. PiscanecPRL, have been confirmed at the point by inelastic X-ray scatteringMaultzsch2004 (), and by the measurement of FWHM(G) in graphite, graphene and nanotubesLazPRB2006 (); piscaPRB (); ACFRaman (); Pisana (), once an-harmonic effects are taken into accountBonini2007 (); ACFRaman (); Pisana (). For the K point, the precise slope of the anomaly is debated Lazzeri2008 (); Graf2007 (); Gruneis2009 (). Another EPC effect is the kink in the electron dispersion, 200 meV below , seen by angle-resolved photoemission spectroscopy (ARPES)Rotenberg2007 (); Zhou2008 (). This is attributed to a correction to the electron energy due to EPCZhou2006 (); Rotenberg2007 (); Zhou2008 (), although alternative explanations also existOlevano (). Thus, a correct EPC determination is a fundamental step for an accurate description of the physical properties of graphene, and nanotubes, being rolled up graphene sheets.

To link the 2D intensity to the EPC we first consider the rate of phonon emission by the photoexcited electron/hole, . This is obtained from the imaginary part of the electron self-energy, . For , as in the case of the Raman measurements at 2.41 eV excitation of Refs.DasCM, ; Dasbila, , we haveBaskoBig ():

Then, from Eq. (2):

(13) |

The dimensionless coupling constants correspond to phonons close to and K, respectively, and determine their rate of emission. We define them as:

(14) |

Here Gruneis2009 () and ,ACFRaman () is the mass of the carbon atom, is the unit cell area. and have the dimensionality of a force and are the proportionality coefficients between the change in effective hamiltonian and the lattice displacement along the corresponding phonon mode. Strictly speaking, the relevant phonon states are not exactly at and K, as shown in Fig. 1. However, the corresponding deviation, , is small compared to the K-K’ distance, and is neglected. All observables depend on the dimensionless EPCs, and .

Eq. (14) follows the notation of Ref. BaskoBig, . Since different EPC definitions are used in the literature, it is quite useful to give here matching rules for all of them, which will be necessary when comparing the EPC values obtained here with previous (and future) reports. The EPCs can be conveniently matched by either relating them to the nearest-neighbor tight-binding model, where the constants are expressed in terms of a single parameter: , the derivative of the nearest-neighbor electronic matrix element with respect to the interatomic distance, or by comparing expressions for various observables. For example, doping leads to a G peak shift due to EPC. This is expressed in terms of asPisana (); YanPrl2007 (); Lazzeri2006 ():

(15) |

The corrections to the phonon dispersions as function of wavevector q, measured from or K, arePiscanecPRL (); piscaPRB (); BaskoBig ():

(16a) | |||

(16b) | |||

(16c) |

Note that the mode splits into longitudinal () and transverse () at finite . Note also that due to analytical properties of the logarithm and square root, Eq. (15) at and Eqs. (16a)–(16c) at acquire imaginary parts, which correspond to the phonon decaying into a continuum of electron-hole pairsLazPRB2006 (). In this case gives the FWHM of the corresponding Lorentzian profile. At Eqs. (16a) and (16c) give the profile of the Kohn anomalies.

In Refs. PiscanecPRL, ; Pisana, ; Calandra2007, ; Lazzeri2008, the EPCs are defined as the matrix elements of the Kohn-Sham potential, differentiated with respect to the phonon displacements. What enters the observables are their squares, averaged over the Fermi surface in the limit . The matching rule is then:

(17a) | |||

In Ref. YanPrl2007, the dimensionless coupling constant is defined as the proportionality coefficient in Eq. (15). Thus,

(18) |

Note that the expression linking EPC to FWHM(G) in Ref. YanPrl2007, underestimates FWHM(G) by a factor 2.

The dimensionless EPC reported in the ARPES analysis of Refs. Rotenberg2007, ; Zhou2008, ; Rotenberg2008, ; GruneisARPES, and in the scanning tunneling spectroscopy (STS) experiment of Ref. EvaAndrei, was measured from the ratio of the electronic velocities below and above the kink in the electron dispersion. This ratio is determined by the derivative of the real part of the electronic self-energy due to the EPC. The latter can be calculated if one takes the Dirac spectrum for electrons and a constant dispersion for phonons. For Calandra2007 ():

(19) | |||||

Here is the ultraviolet cutoff, of the order of the electronic bandwidth. We then get the matching rule:

(20) | |||||

However, we note that is subject to Coulomb renormalizationsBaskoAleiner (). This implies that depends on the electronic energy scale, such as the electron energy , the Fermi energy , or the temperature , whichever is larger: . This dependence is shown in Fig. 6 of Ref. BaskoAleiner, in the semi-logarithmic scale. In a Raman measurement this scale is given by the energy of the photo-excited electron: , as long as . Thus, in Eq. (13) . On the other hand, to estimate the EPC effects on the phonon dispersions in the intrinsic graphene, the relevant electron energy is of the order of the phonon energy. Thus, in Eq. (16c) . From Fig. 6 of Ref. BaskoAleiner, we estimate that for and 1.2 for (taking eV to represent Raman measurements in the visible range).

The situation with Eq. (20) is more complicated, since the cutoff appears explicitly. The logarithmic term is determined by all energy scales from down to . Thus, the proper expression is

(21) | |||||

### iii.2 Experimental Electron-Phonon Coupling

From Eq.(12), our overall average , derived from a fit to all the data in Fig. 4, gives:

(22) |

On the other hand, the hole doping side of Fig. 4 shows two data sets very consistent with each other. We can thus get another estimate taken from the average for just the hole doping side. This would give:

(23) |

Based on measurementsPisana (); YanPrl2007 () and DFT calculationsPiscanecPRL (), the value of can be reliably taken . Indeed, DFT gives PiscanecPRL () and , corresponding, from Eqs. (14), (17a) to . Even though and are subject to Coulomb renormalization, , which contains their ratio, is not.BaskoAleiner () The experimental extracted from FWHM(G) in graphene and graphiteACFRaman (); LazPRB2006 () according to Eq. (16a) and from the dependence of Pos(G) on Fermi energy according to Eq. (15), give YanPrl2007 () and Pisana ().

On the other hand, the value of is still debatedCalandra2007 (); BaskoAleiner (); Lazzeri2008 (). The calculated DFT , together with the DFT (both taken from Ref. PiscanecPRL, ) gives . However, Ref. BaskoAleiner, suggested this should be enhanced by Coulomb renormalization by up to a factor 3, depending on the background dielectric constant. In order to compare with our fits, we need consider that the corrections to the phonon dispersion are determined by electronic states with energies lower than those contributing to the Raman signal. As discussed in Sec. IIIA, for . Our fit in Eq. (22) corresponds to , while Eq. (23) gives , resulting in and , respectively. These are bigger than DFT by a factor of about 3.5 and 1.7, respectively.

A recent GW calculation gave Lazzeri2008 (). Combining this with the GW Gruneis2008 (), we get , a factor greater than DFT, in good agreement with our fitted average on the hole side.

Ref. Gruneis2009, reported inelastic x-ray scattering measurements of the phonon dispersions near K more detailed than those originally done in Ref. Maultzsch2004, , now giving a phonon slope at K of . Using Eq. (16c) at and taking the experimental value Zhou2008 () (the bare electron velocity, i. e. below the phonon kink), we obtain , a factor higher than DFT, again in good agreement with our fitted average on the hole side.

Another EPC estimate can be derived from the 2D and 2D’ area ratio. Combining Eqs. (1a),(1b),12,13 we get:

(24) |

For intrinsic SLG and graphite whiskers, the experimental A(2D)/A(2D’) is ACFRaman (); heinz (); berciaud (); tan1 (), which gives and . Since in this case , this results in , a factor 4.5 higher than DFT, in agreement with our upper estimate from Eq. (22).

We finally consider the EPC derived from ARPES and STS. For an estimate, we approximate the dependence as linear in . We take , as given by DFT (assumed to be valid at high energies), and leave as the only free parameter determining this linear dependence:

(25) |

Taking Rotenberg2007 (); Zhou2008 (); EvaAndrei (); GruneisARPES (), , and substituting Eq. (25) in Eq. (21), we get:

(26) |

Note that the dependence on the precise value of is weak: setting changes the first coefficient to , and the second (more important as it multiplies the larger coupling constant) varies only by 2%. The measurements in Refs. GruneisARPES, ; Rotenberg2007, ; EvaAndrei, ; Rotenberg2008, ; Zhou2008, gave , respectively. The smallest of these values, , from Eq. (26) corresponds to , while the highest to . Even the smallest is almost twice our upper bound fit of Eq. (22) and would imply an EPC renormalization of almost one order of magnitude. Resolution effects could play a role in this overestimationCalandra2007 ().

Thus, our fits to the doping dependent Raman area ratios point to a significant renormalisation, by a factor 1.7-3.5, of the TO mode close to K, responsible for the Raman D and 2D peaks. Our lower bound estimate is consistent with recent GW calculations and phonon measurements, but our upper bound is much lower than the smallest estimate derived by ARPES.

## Iv Conclusions

We have shown that the 2D intensity dependence on doping can be explained considering the influence of electron-electron interactions on the total scattering rate of the photogenerated electrons (holes). We have given a simple formula linking 2D peak area to the Fermi level shift. Fitting this to the available experimental data we got an estimate for the EPC value of the TO phonons close to K, responsible for the Raman D and 2D peaks. This is larger than that from DFT calculations, due to renormalisation by Coulomb interactions. However, our fitted EPC is still significantly smaller than those reported in ARPES or STS experiments.

## V Acknowledgments

We acknowledge A. Das, S. Berciaud, A. Bonetti, P.H. Tan for useful discussions. A.C.F. acknowledges funding from the Royal Society and the European Research Council grant NANOPOTS.

## References

- (1) K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. Zhang, S. V.Dubonos, I. V. Grigorieva, A. A. Firsov; Science, 306, 666 (2004).
- (2) A. K. Geim, K. S. Novoselov; Nature Mater. 6, 183 (2007).
- (3) A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov, A. K. Geim; Rev. Mod. Phys. 81, 109 (2009).
- (4) J. C. Charlier, P.C. Eklund, J. Zhu, A.C. Ferrari, Topics Appl. Phys.111, 673 (2008).
- (5) K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, M. I. Katsnelson, I. V. Grigorieva, S. V. Dubonos, A. A. Firsov; Nature (London), 438, 197 (2005).
- (6) Y. Zhang, Y.W. Tan, H. L. Stormer, P. Kim; Nature (London), 438, 201 (2005).
- (7) K. S. Novoselov, Z. Jiang, Y. Zhang, S. V. Morozov, H. L. Stormer, U. Zeitler, J. C. Maan, G. S. Boebinger, P. Kim, A. K. Geim; Science, 315, 1379 (2007).
- (8) S. V. Morozov, K. S. Novoselov, M. I. Katsnelson, F. Schedin, D. C. Elias, J. A. Jaszczak, A. K. Geim; Phys. Rev. Lett., 100, 016602 (2008).
- (9) X. Du, I. Skachko, A. Barker, E. Y. Andrei, Nature Nano. 3, 491 (2008)
- (10) K. I. Bolotin, K. J. Sikes, J. Hone, H. L. Stormer, P. Kim Phys. Rev. Lett. 101, 096802 (2008) ; K. I. Bolotin, K. J. Sikes, Z. Jiang, G. Fundenberg, J. Hone, P. Kim, H. L. Stormer, Solid State Comm. 146, 351 (2008).
- (11) M. Y. Han, B. Oezylmaz, Y. Zhang, P. Kim; Phys. Rev. Lett., 98, 206805 (2007).
- (12) Z. Chen, Y.M. Lin, M. Rooks, P. Avouris; Physica E, 40, 228 (2007).
- (13) Y. Zhang, J. P. Small, W. V. Pontius, P. Kim; Appl. Phys. Lett., 86, 073104 (2005).
- (14) M. C. Lemme, T. J. Echtermeyer, M. Baus, H. Kurz; IEEE El. Dev. Lett.,28, 4 (2007).
- (15) Y.M. Lin, K. A. Jenkins, A. Valdes-Garcia, J. P. Small, D. B. Farmer, P. Avouris, Nano Lett. 9, 422 (2009)
- (16) J. S. Bunch, A. M. van der Zande, S. S. Verbridge, I. W. Frank, D. M. Tanenbaum, J. M. Parpia, H. G. Craighead, P. L. McEuen, Science 315, 490 (2007).
- (17) P. Blake, P. D. Brimicombe, R. R. Nair, T. J. Booth, D. Jiang, F. Schedin, L. A. Ponomarenko, S. V. Morozov, H. F. Gleeson, E. W. Hill, A. K. Geim, K. S. Novoselov, Nano Lett. 8, 1704 (2008)
- (18) Y. Hernandez, V. Nicolosi, M. Lotya, F. Blighe, Z. Sun, S. De, I. T. McGovern, B. Holland, M. Byrne, Y. Gunko, J. Boland, P. Niraj, G. Duesberg, S. Krishnamurti, R. Goodhue, J. Hutchison, V. Scardaci, A. C. Ferrari, J. N. Coleman, Nature Nano. 3, 563 (2008)
- (19) G. Eda, G. Fanchini, M. Chhowalla, Nature Nano. 3, 270 (2008).
- (20) A. C. Ferrari, J. C. Meyer, V. Scardaci, C. Casiraghi, M. Lazzeri, F. Mauri, S. Piscanec, Da Jiang, K. S. Novoselov, S. Roth, A. K. Geim; Phys. Rev. Lett., 97, 187401 (2006).
- (21) C. Casiraghi, A. Hartschuh, E. Lidorikis, H. Qian, H. Harutyunyan, T. Gokus, K. S. Novoselov, A. C. Ferrari; Nano. Lett., 7, 2711 (2007).
- (22) P. Blake, E. W. Hill, A. H. Castro Neto, K. S. Novoselov, D. Jiang, R. Yang, T. J. Booth, A. K. Geim; Appl. Phys. Lett., 91, 063124 (2007).
- (23) C. Casiraghi, S. Pisana, K. S. Novoselov, A. K. Geim, A. C.Ferrari; Appl. Phys. Lett., 91, 233108 (2007).
- (24) L. M. Malard, J. Nilsson, D. C. Elias, J. C. Brant, F. Plentz, E. S. Alves, A. H. Castro Neto, M. A. Pimenta, Phys. Rev. B 76, 201401 (1999)
- (25) S. Pisana, M. Lazzeri, C. Casiraghi, K. S. Novoselov, A. K. Geim, A. C. Ferrari, F. Mauri, Nat. Mat. 6, 198 (2007).
- (26) S. Piscanec, M. Lazzeri, F. Mauri, A. C. Ferrari, J. Robertson; Phys. Rev. Lett., 93, 185503 (2004).
- (27) L. G. Cançado, R. Beams, L. Novotny Cond.Mat. 0802.2709
- (28) C. Casiraghi, A. Hartschuh, H. Qian, S. Piscanec, C. Georgi, A. Fasoli, K. S. Novoselov, D. M. Basko, A. C. Ferrari, Nano Lett. 9, 1433 (2009)
- (29) D. C. Elias, R. R. Nair, T. M. G. Mohiuddin, S. V. Morozov, P. Blake, M. P. Halsall, A. C. Ferrari, D. W. Boukhvalov, M. I. Katsnelson, A. K. Geim, K. S. Novoselov, Science 323, 610 (2009).
- (30) A. C. Ferrari; Solid State Comm., 143, 47 (2007).
- (31) A. Das, S. Pisana, S. Piscanec, B. Chakraborty, S. K. Saha, U. V. Waghmare, R. Yang, H. R. Krishnamurhthy, A. K. Geim, A. C. Ferrari, A. K. Sood Nature Nano. 3, 210 (2008)
- (32) A. Das, B. Chakraborty, S. Piscanec, S. Pisana, A. K. Sood, A. C. Ferrari, Phys Rev B 79, 155417 (2009)
- (33) N. Ferralis, R. Maboudian, C. Carraro, Phys. Rev. Lett. 101, 156801 (2008)
- (34) T. M. G. Mohiuddin, A. Lombardo, R. R. Nair, A. Bonetti, G. Savini, R. Jalil, N. Bonini, D. M. Basko, C. Galiotis, N. Marzari, K. S. Novoselov, A. K. Geim and A. C. Ferrari, Phys. Rev. B 79, 205433 (2009)
- (35) J. Yan, Y. Zhang,P. Kim, A. Pinczuk, Phys.Rev.Lett. 98, 166802 (2007).
- (36) D. Graf, F. Molitor, K. Ensslin, C. Stampfer, A. Jungen, C. Hierold, and L. Wirtz, Nano Lett. 7, 238 (2007)
- (37) A.C. Ferrari, J. Robertson Phys. Rev. B 61, 14095 (2000); ibid. 64, 075414 (2001).
- (38) A. C. Ferrari, J. Robertson (eds), Raman spectroscopy in carbons: from nanotubes to diamond, Theme Issue, Phil. Trans. Roy. Soc. A 362, 2267-2565 (2004).
- (39) C. Thomsen, S. Reich, Phys. Rev.Lett. 85, 5214 (2000).
- (40) F.Tuinstra,J.L. Koenig, J. Chem. Phys.53, 1126 (1970).
- (41) D. M. Basko, Phys. Rev. B 76, 081405(R) (2007).
- (42) D. M. Basko, Phys. Rev. B 78, 125418 (2008).
- (43) M. Lazzeri, F. Mauri, Phys. Rev. Lett. 97, 266407 (2006).
- (44) M. Lazzeri, S. Piscanec, F. Mauri, A. C. Ferrari, J. Robertson, Phys. Rev. B 73, 155426 (2006).
- (45) A. Das, A. K. Sood, A. Govindaraj, A. M. Saitta, M. Lazzeri, F. Mauri, C. N. R Rao, Phys Rev Lett. 99, 136803 (2007).
- (46) J.Yan, E. A. Henriksen, P. Kim, A. Pinczuk, Phys. Rev. Lett. 101, 136804 (2008).
- (47) T. Ando, J. Phys. Soc. Jpn. 76, 104711 (2007).
- (48) S. Berciaud, S. Ryu, L. E. Brus, T. F. Heinz, Nano Lett., 9, 346 (2009)
- (49) S. Berciaud, S. Ryu, L. E. Brus, T. F. Heinz, APS March L26.00006, (2009)
- (50) P. H. Tan, C. Y. Hu, J. Dong, W. C. Shen, B. F. Zhang, Phys. Rev B, 64, 214301 (2001)
- (51) J. Martin, N. Akerman, G. Ulbricht, T. Lohmann, J. H. Smet, K. von Klitzing, A. Yacoby, Nat. Phys. 4, 148 (2008).
- (52) L.I. Mandelshtam, G. S. Landsberg, Z. Phys. 50, 169 (1928); C. V. Raman, K. S. Krishnan, Nature 121, 501 (1928); ibid. 121, 619 (1928).
- (53) N. Bonini, M. Lazzeri, N. Marzari, F. Mauri, Phys. Rev. Lett. 99 176802 (2007).
- (54) M.Lazzeri, S. Piscanec, F. Mauri, A. C. Ferrari, J. Robertson, Phys Rev. Lett. 95, 236802 (2005)
- (55) R. P. Vidano, D. B. Fischbach, L. J. Willis, T. M. Loehr, Solid State Commun. 39, 341 (1981).
- (56) I. Pócsik, M. Hundhausen, M. Koos, L. Ley, J. Non-Cryst. Solids 227–230, 1083 (1998).
- (57) L. G. Cançado, A. Jorio, M. A. Pimenta, Phys. Rev. B 76, 064304 (2007).
- (58) Z. Jiang, E. A. Henriksen, L. C. Tung, Y.-J. Wang, M. E. Schwartz, M. Y. Han, P. Kim, H. L. Stormer, Phys. Rev. Lett. 98, 197403 (2007).
- (59) A. Bostwick, Solid State Commun. 143, 63 (2007).
- (60) S. Y. Zhou, D. A. Siegel, A. V. Fedorov, A. Lanzara, Phys Rev. B 78, 193404 (2008).
- (61) L. M. Zhang, M. M. Fogler, Phys. Rev. Lett. 100, 116804 (2008).
- (62) M. M. Fogler, D. S. Novikov, L. I. Glazman, and B. I. Shklovskii, Phys. Rev. B 77, 075420 (2008).
- (63) E. Rossi, S. Das Sarma, Phys. Rev. Lett. 101, 166803 (2008).
- (64) M. Polini, A. Tomadin, R. Asgari, and A. H. MacDonald, Phys. Rev. B 78, 115426 (2008).
- (65) J. González, F. Guinea, M. A. H. Vozmediano, Phys. Rev. Lett. 77, 3589 (1996). To be precise, the imaginary part of the electronic self-energy has a discontinuity on the mass shell. This discontinuity is resolved by taking into account the real part as well, which makes the electronic spectrum concave.
- (66) A. A. Abrikosov, S. D. Beneslavskii, Zh. Eksp. Teor. Fiz 59, 1280 (1970) [Sov. Phys. JETP 32, 699 (1971)].
- (67) J. González, F. Guinea, M. A. H. Vozmediano, Mod. Phys. Lett. B 7, 1593 (1994); Nucl. Phys. B 424, 595 (1994); J. Low Temp. Phys. 99, 287 (1994).
- (68) E. H. Hwang, B. Y.-K. Hu, S. Das Sarma, Phys. Rev. B 76, 115434 (2007).
- (69) M. Polini, R. Asgari, G. Borghi, Y. Barlas, T. Pereg-Barnea, A. H. MacDonald, Phys. Rev. B 77, 081411 (2008).
- (70) S. Latil, V. Meunier, L. Henrard, Phys. Rev. B 76, 201402 (2007)
- (71) P. R. Wallace, Phys. Rev. 71 622 (1947).
- (72) J. Maultzsch, S. Reich, C. Thomsen, H. Requardt, P. Ordejón, Phys. Rev. Lett. 92, 075501 (2004).
- (73) S. Piscanec, M. Lazzeri, J. Robertson, A. C. Ferrari, F. Mauri, Phys. Rev. B 75, 035427 (2007)
- (74) M. Lazzeri, C. Attaccalite, L. Wirtz, and F. Mauri, Phys. Rev. B 78, 081406 (2008).
- (75) A. Grüneis, J. Serrano, A. Bosak, M. Lazzeri, S.L. Molodtsov, L. Wirtz, C. Attaccalite, M. Krisch, A. Rubio, F. Mauri, T. Pichler, arXiv:0904.3205 (2009)
- (76) A. Bostwick, T. Ohta, T. Seyller, K. Horn, E. Rotenberg, Nature Phys. 3, 36 (2007).
- (77) S.Y. Zhou, G. H. Gweon, A. Lanzara, Ann. Phys. (N.Y.)., 321, 1730 (2006)
- (78) P. E. Trevisanutto, C. Giorgetti, L. Reining, M. Ladisa, V. Olevano, Phys. Rev. Lett. 101, 226405 (2008).
- (79) M. Calandra, F. Mauri, Phys. Rev. B 76, 205411 (2007).
- (80) C. S. Leem, B. J. Kim, Chul Kim, S. R. Park, T. Ohta, A. Bostwick, E. Rotenberg, H. -D. Kim, M. K. Kim, H. J. Choi, and C. Kim, Phys. Rev. Lett. 100, 016802 (2008).
- (81) A. Grüneis, C. Attaccalite, A. Rubio, D. Vyalikh, S.L. Molodtsov, J. Fink, R. Follath, W. Eberhardt, B. Büchner, and T. Pichler, Phys. Rev. B 79, 205106 (2009).
- (82) G. Li, A. Luican, E. Y. Andrei, Phys. Rev. Lett. 102, 176804 (2009).
- (83) D.M. Basko, I.L. Aleiner, Phys. Rev. B 77, 041409(R) (2008).
- (84) A. Grüneis, C. Attaccalite, T. Pichler, V. Zabolotnyy, H. Shiozawa, S. L. Molodtsov, D. Inosov, A. Koitzsch, M. Knupfer, J. Schiessling, R. Follath, R. Weber, P. Rudolf, L. Wirtz, A. Rubio, Phys. Rev. Lett. 100, 037601 (2008).