Best Student Paper Award Efficient Auger scattering in Landau-quantized graphene F. Wendler∗a, H. Funka, M. Mittendorffb, S. Winnerlc, M. Helmc,d, A. Knorra and E. Malice aInstitute of Theoretical Physics, Technical University Berlin, Berlin, Germany; bInstitute for Research in Electronics and Applied Physics, University of Maryland, College Park, Maryland, USA; cHelmholtz-Zentrum Dresden-Rossendorf, Dresden, Germany; dTechnische Universit¨atDresden, Dresden, Germany; eDepartment of Applied Physics, Chalmers University of Technology, G¨oteborg, Sweden ABSTRACT We present an analytical expression for the differential transmission of a delta-shaped light field in Landau- quantized graphene. This enables a direct comparison of experimental spectra to theoretical calculations reflect- ing the carrier dynamics including all relevant scattering channels. In particular, the relation is used to provide evidence for strong Auger scattering in Landau-quantized graphene. Keywords: Landau-quantized graphene, Auger scattering, carrier dynamics, differential transmission spec- troscopy 1. INTRODUCTION Since its discovery in 2004,1 graphene has attracted an enormous amount of interest2 being considered as a novel material for optoelectronic devices.3{8 To exploit the fascinating electronic and optical properties of this atomically thin material, it is of crucial importance to understand the non-equilibrium carrier dynamics. There has been a lot of research in this field,9{27 revealing that carrier-carrier and carrier-phonon scattering dominate the ultrafast relaxation dynamics.28, 29 Remarkably, Auger scattering { which is generally suppressed in ordinary semiconductors { turns out to be of particular importance in graphene, where it gives rise to carrier multiplication occurring at certain experimentally controllable conditions.23, 30{33 In an external magnetic field, the electronic properties change dramatically, as the energy is quantized into non-equidistant Landau levels.34 Unlike in the case without a magnetic field, there are only a few studies of the carrier dynamics in Landau-quantized graphene.35{41 Here, optical selection rules42 allow the excitation of specific inter-Landau level transitions opening the possibility to address the carrier dynamics of selected levels in a pump-probe experiment. Plochocka et al. have investigated the carrier relaxation in highly energetic Landau levels (with an index n ∼ 100) which was found to be slowed down by the magnetic field as a consequence of the energy quantization.35 In a more recent joint experiment theory study, we have investigated the carrier dynamics within the energetically lowest Landau levels.36 Based on a pump probe experiment and microscopic modeling of the carrier dynamics, we were able to provide evidence for strong Auger-type scattering in Landau-quantized graphene: Including this scattering channel in the theoretical model, an unexpected behavior in the measured differential transmission could be explained. In this Article, we derive the analytic expression that was used in Ref. [36] to compare the experimental differential transmission spectra (DTS) with the calculations of the Landau level occupations. To this end, we introduce the Bloch equations for Landau-quantized graphene in Section 2, calculate the absorbance for a delta- shaped probe pulse in Section 3, and ultimately derive an analytic expresion for the differential transmission in Section 4. ∗fl[email protected], Telephone: +49 (0)30 314 23029 Ultrafast Phenomena and Nanophotonics XIX, edited by Markus Betz, Abdulhakem Y. Elezzabi, Kong-Thon Tsen, Proc. of SPIE Vol. 9361, 936105 · © 2015 SPIE CCC code: 0277-786X/15/$18 · doi: 10.1117/12.2075458 Proc. of SPIE Vol. 9361 936105-1 2. BLOCH EQUATIONS FOR LANDAU-QUANTIZED GRAPHENE The many-particle Hamilton operator for charge carriers in Landau-quantized graphene subject to an optical light field is given by X y X y H = iai ai + i~ Ωif af ai; (1) i i;f p λ where the first term is the free energy part that is determined by the Landau level dispersion n = λvF 2~ne0B 43 which depends on the band index λ = ±1, the Fermi velocity vF ≈ 1 nm=fs, the Landau level index n = 0; 1; 2;:::, the elementary charge e0, and the magnetic field B. The second term describes the carrier-light e0 coupling expressed through the Rabi frequency Ωif (t) = Mif ·A(t) with the free electron mass m0, the vector m0 potential A(t), and the optical matrix element38, 44 αni αnf m0vF − + ∗ Mif = iδξi,ξf δmi;mf p λi^ δnf ;ni−1 + λf ^ δnf ;ni+1 = −Mfi: (2) 2 2~ Here, ξ = ±1 is the valley index that distinguishes the two equivalent Dirac cones of graphene.34 Furthermore, 2 m = 0; 1;:::;L e0B=(2π~) − 1 is a quantum number expressing thep Landau level degeneracy with the areap of 2 34 ± graphene L , αn=0 is a constant that can take the values αn=0 = 2 and αn6=0 = 1, and ^ = (e^x ∓ ie^y)= 2 are Jones vectors describing left (+) and right (−) circularly polarized light.45 y y The temporal evolution of the microscopic polarization pif (t) = haf aii(t) (where aj and aj are fermionic d creation and annihilation operators of the state j) is obtained using the Heisenberg equation of motion i~ dt O(t) = [O;H] which provides an equation of motion for the operator O.28, 46 Using the Hamilton operator from Eq. 1, this yields ∗ p_if (t) = (i4!if − γ) pif (t) − Ωif [ρf (t) − ρi(t)] ; (3) with the energy difference ~4!if = f − i, a phenomenologically introduced dephasing γ, and the occupation probability ρj = pjj in the state j. The dephasing is equivalent to a Landau level broadening and can be caused by disorder47 or, alternatively by interaction with acoustic phonons.48, 49 3. ABSORBANCE OF A DELTA-SHAPED PULSE For a two-dimensional nanostructure in the linear optics regime, the absorbance reads26, 28, 29 !Imχ±(!) α±(!) = ; (4) c j1 + i!χ±(!)=2cj2 where c is the light velocity and ± ± j (!) χ (!) = 2 ± (5) 0! A (!) is the linear optical susceptibility, with the current density j, and the permittivity 0. The superscript ± denotes the σ±-circularly polarized components of the respective quantity in the Jones vector basis (^±). The components 1 h i A±(!) = p A+ei(!∓∆ωif )τ + A−ei(!±∆ωif )τ ; (6) 2 π 0 0 of a δ-shaped probe pulse that is in resonance with the transition i ! f is obtained from + cos 4!if t − cos 4!if t A(t) = δ (t + τ) A0 + A0 ; (7) sin 4!if t − sin 4!if t Proc. of SPIE Vol. 9361 936105-2 ± where τ describes the temporal delay of the probe pulse, and A0 are the amplitudes of the left (+) and right (−) circularly polarized contributions. The current density is given by28 ± ~e0 X h ± ± i j (!) = −i 2 Mif pif (!) + Mfipfi(!) : (8) 2m0L if ± rot i4!if t where a term proportional to A was neglected. Using the substitution pif = pif e , we transform Eq. 3 to the rotating frame rot rot ∗ −i4!if t p_if (t) = −γpif (t) − Ωif [ρf (t) − ρi(t)] e ; (9) h rot i rot and obtain its Fourier transform (FT) exploiting the relation FT p_if (t) = i!pif (!) rot e0vF − + FT [δ (t − τ)(ρf − ρi)] (!) pif (!) = δξi,ξf δmi;mf αni αnf λiδnf ;ni−1A0 + λf δnf ;ni+1A0 : (10) 4~ ! − iγ Here, the rotating wave approximation (RWA) was applied by omitting terms e±i24!if t in comparison to 1. ∗ rot −i4!if t Since the symmetry relation pif = pfi effects an inverse rotation of pfi compared to pif , cf. pfi = pfi e , rot rot rot the RWA for pfi (!) is not obtained by exchanging i $ f in Eq. 10, but the relation pfi (!) = −pif (!) must be used instead. Exploiting the convolution theorem, we can write FT [δ (t + τ)(ρf − ρi)] = FT [δ (t − τ)] ∗ FT [ρf − ρi] 1 = ei!τ [ρ (τ) − ρ (τ)] : (11) 2π f i 1 0 0 0 where ∗ denotes the convolution of two functions f and g:(f ∗ g)(!) = −∞ d! f(! )g(! − ! ), Using Eqs. 2, 8, rot rot 10, 11, as well as the relations pif (!) = pif (! − 4!if ) and pfi(!) = −p´if (! + 4!if ), the current density reads 3 2 " i(!−∆ω )τ i(!+∆ω )τ # e v B X 2 e if e if j±(!) = 0 F α α [ρ (τ) − ρ (τ)] A±δ − A∓δ 64π2 2 ni nf f i 0 nf ;ni±1 imp 0 nf ;ni∓1 imp ~ if ! − ∆ωif − iγif ! + ∆ωif − iγif (12) Here, we performed the summations over the valley and m degrees of freedom. In due consideration of the 2 Kronecker deltas δξi,ξf δmi;mf in Eq. 10 this yields a factor of L e0B=(π~). Omitting terms with a non vanishing rotation ei!τ~ 6= 1, we obtain the linear optical susceptibility (Eq. 5) using the Eqs. 6 and 12 p 3 2 ± πe0vFB X 2 ρf (τ) − ρi(τ) χ (!) = 2 2 2 αni αnf imp δnf ;ni±1: (13) 32π 0! ~ if ! − ∆ωif − iγif Neglecting the term i!χ(!)=2c in the denominator of Eq. 4 which is usually small, we find the final expression for the absorbance of a delta-shaped pulse p 3 2 ± πe0vFB X 2 γ α (!) = 2 2 αni αnf [ρf (τ) − ρi(τ)] 2 δnf ;ni±1: (14) 32π 0c! 2 ~ if (! − ∆ωif ) + γ 4. DIFFERENTIAL TRANSMISSION The idea of the differential transmission spectroscopy is to optically excite the system under investigation using a strong pump pulse. Then, the differential transmission of a much weaker probe pulse incident after a certain time delay τ is measured. Based on the assumption that the occupations of the system are not affected by the weak probe pulse, the differential transmission allows to draw conclusions about the carrier dynamics in the system.