Reader TTS

Download PDF Article Open access Published: 10 June 2020 Unconventional valley-dependent optical selection rules and landau level mixing in bilayer graphene Long Ju1,2 na1 nAff8, Lei Wang 

ORCID: orcid.org/0000-0002-1919-91071,2 na1, Xiao Li  ORCID: orcid.org/0000-0002-1708-43143, Seongphill Moon4, Mike Ozerov  ORCID: orcid.org/0000-0002-5470-11584, Zhengguang Lu4, Takashi

Taniguchi5, Kenji Watanabe  ORCID: orcid.org/0000-0003-3701-81195, Erich Mueller  ORCID: orcid.org/0000-0001-5186-58892, Fan Zhang  ORCID: orcid.org/0000-0003-4623-42006, Dmitry Smirnov 

ORCID: orcid.org/0000-0001-6358-32214, Farhan Rana7 & …Paul L. McEuen  ORCID: orcid.org/0000-0002-6438-57501,2 Show authors Nature Communications volume 11, Article number: 2941 (2020) Cite

this article

8104 Accesses

17 Citations

2 Altmetric

Metrics details

Subjects Materials scienceNanoscience and technologyOptics and photonicsPhysics An Author Correction to this article was published on 24 June 2020

This article has been updated

Selection rules are of vital importance in determining the basic optical properties of atoms, molecules and semiconductors. They provide general insights into the symmetry of the system and

the nature of relevant electronic states. A two-dimensional electron gas in a magnetic field is a model system where optical transitions between Landau levels (LLs) are described by simple

selection rules associated with the LL index N. Here we examine the inter-LL optical transitions of high-quality bilayer graphene by photocurrent spectroscopy measurement. We observed

valley-dependent optical transitions that violate the conventional selection rules Δ|N| = ± 1. Moreover, we can tune the relative oscillator strength by tuning the bilayer graphene bandgap.

Our findings provide insights into the interplay between magnetic field, band structure and many-body interactions in tunable semiconductor systems, and the experimental technique can be

generalized to study symmetry-broken states and low energy magneto-optical properties of other nano and quantum materials.

single-electron bilayer graphene quantum dots Article Open access 02 September 2021 Broken symmetries and excitation spectra of interacting electrons in partially filled Landau levels

Article 20 July 2023 Spin–valley protected Kramers pair in bilayer graphene Article Open access 10 February 2025 Introduction

Bilayer graphene (BLG) has emerged as a two-dimensional semiconductor where the bandgap is tunable by an external electric field. At zero magnetic field, the bandgap hosts strong exciton

resonances1 that obey unusual optical selection rules determined by the electron pseudospin texture. At large magnetic fields, two-dimensional electron gas is expected to form quantized

energy levels that were named after Landau2. Properties of individual Landau levels (LLs) have been studied extensively ever since the birth of graphene using electron transport3,4,5,6,7,

scanning tunneling spectroscopy8, and electronic compressibility measurements9. Inter-LL optical transitions in graphene could provide direct information of the critical energy scales of

SU(4) isospin6 degeneracy lifting through spectroscopy and the wave functions of LLs through selection rules. For semiconductor physics, BLG provides a model system where the interplay

between band dispersion, cyclotron energy, and Coulomb energy of electron-hole pairs can be systematically studied through gating and tuning magnetic field. However, inter-LL transitions in

BLG and the crossover from exciton-dominated to LL-dominated optical response was largely unexplored, let alone the optical selection rules of them.

Previous infrared absorption spectroscopy studies of MLG (both on SiO2/Si10,11 and on hBN12,13), BLG on SiO2/Si substrates14, and thin graphene layers on SiC15,16 suffered from broad peak

width or the lack of gating. And they revealed only optical transitions obeying the conventional Δ|N | = ± 1 selection rule. Further investigations on higher quality graphene samples,

however, have been challenging owing to limited sizes of hBN-encapsulated graphene samples compared with the relevant infrared wavelength. Here we overcome this difficulty by using Fourier

Transformed Infrared photocurrent spectroscopy1 to study high-quality hBN-encapsulated BLG devices. In addition to greatly enhanced signal-to-noise ratio, such spectroscopy method avoids the

background signal owing to absorption in the gates—providing clean data from the layer of interest. Although photocurrent signal is more complicated than optical absorption, it would not

affect the main conclusion of this article as the conversion efficiency should be similar for optical transitions with similar energies. As a result, we observed spectral features with a

linewidth of ~1 meV and many optical transitions that violate the conventional selection rule10,12,14,17,18,19 Δ|N| = ± 1.

Figure 1a shows the band structure of BLG near the K point. Around the electric field-induced bandgap Δ both conduction and valence bands develop three pockets20,21 owing to the trigonal

warping effect22. Figure 1b presents LL spectra in BLG calculated based on a continuum model with a bandgap of 80 meV (optimized for the spectrum at D = 0.874 V/nm, see Supplementary

Discussion), where each LL is labeled by its orbital index N. Owing to the interlayer potential, LL with N = 0 and 1 have a valley-dependent distribution. Together with electron-hole

asymmetry inherent in the band structure, inter-LL transitions in K and K′ valleys (especially those close to the bandgap) are expected to have strong differences. Electric-dipole-allowed

optical selection rules dictates that initial and final states are different by ± ħ in angular momentum, corresponding to that of circularly polarized photons. In ordinary semiconductors,

this angular momentum conservation is satisfied by interband transitions at zero magnetic field. Therefore, interband LL transitions are allowed only when Δ|N| = 0 as no additional angular

momentum change is induced by LL wave functions. In graphene systems10,12,14, transition metal dichalcogenides17 and other Dirac electron systems18,19, however, an extra pseudospin degree of

freedom introduces an effective angular momentum for electrons. As a result, the conventional optical transition selection rules were shown to be Δ|N| = ± 1 both theoretically and

experimentally10,12,14,23 as illustrated by solid and dashed arrows in Fig. 1b.

a Illustration of the trigonal-warped BLG band structure with a bandgap Δ in the K valley. b Evolution of Landau level energies calculated by a continuum model with Δ = 80 meV. We use red

and green colors to label K and K′ valleys, respectively, throughout the paper. Solid and dashed arrows illustrate inter-LL transitions obeying the usual selection rule of Δ|N| = ± 1. c

Photocurrent spectra with a displacement field of D = 0.874 V/nm at B = 0, 2, and 4 T. The 0 T spectrum features two exciton peaks near the bandgap and a flat spectrum above the bandgap.

Oscillations corresponding to inter-LL transitions emerge at 2 T with an alternated high-low pattern. The contrast between high and low peaks increases at 4 T, where splittings of peaks also

appear. d A 2D color plot of photocurrent spectra as the magnetic field is continuously tuned. Dashed curves trace high peaks in c, which grow stronger towards higher magnetic fields and we

label them by branch A-H (only branch A-D are shown in c). These strong transitions obey the usual optical selection rule of Δ |N| = ± 1. The low peaks in c evolve into weak transitions

between dashed lines, whose optical selection rules will be discussed later.

Our hBN-BLG-hBN sample is similar to that used in ref. 1. Figure 1c presents photocurrent spectra at a displacement field D = 0.874 V/nm with the Fermi level tuned into the bandgap (see ref.

 1 for the measurement technique). At B = 0 T, the spectrum features two sharp peaks, which were previously identified as excitons1. Above excitons, the spectrum starts as a flat line and

gradually develops oscillations with an increasing B field. At 2 T, the oscillations show a pattern of alternating higher and lower peaks. The contrast between these two groups of peaks is

enhanced at B = 4 T. At the same time, both higher and lower peaks split into finer structures (the latter will be shown more clearly in following figures). These spectra present the

continuous evolution of the system behavior from exciton-dominated picture to inter-LL-dominated picture.

To better visualize the evolution of all transition peaks, we plot spectra as a 2D color map in Fig. 1d as B is continuously tuned. High peaks in Fig. 1c evolve into branches in the diagonal

direction as traced by the dashed curves. We label these branches by A-H. The lower peaks in Fig. 1b evolve into branches between the dashed curves. The contrast in intensity of these two

groups of branches increases as B rises. We will call them strong and weak transitions hereafter.

Discrete peaks in photocurrent spectrum correspond to optical transitions between LLs and Fig. 1d reveals a rich collection of spectral information. To better understand the origins of

various transitions, we first focus on branch C as shown in Fig. 2a. We shifted spectra horizontally to align the center of transition peaks to zero energy. This plot clearly shows a pair of

peaks as traced by dashed lines, which split more as B increases.

a Evolution of branch C transitions. For clarity, photocurrent spectra at 2–5 T (0.2 T step) are shifted horizontally to the center and vertically by multiples of 0.2. Dashed lines trace

transitions corresponding to solid arrows in Fig. 1b. Transitions corresponding to dashed arrows in Fig. 1b are too weak to be identified. b Upper panel: the splitting between two peaks in

a, E45. It measures the difference between electron and hole cyclotron energies, ΔEc, which scales linearly with B. Lower panel: the ratio between ΔEc and the average cyclotron energy

\(\overline {E_c} \), indicating an asymmetry of ~18% between electron and hole masses. c Evolution of branch B, which is special among all branches as it involves three transitions. The

extra transition corresponds to transition (0 → 4, K′). Such Δ|N| = ±4 transitions are allowed by the trigonal warping effect. d Evolution of the 2p exciton. Left panel: Fig. 1d near and

below branch A. Right panel: calculated dispersions of LL transitions. When LLs are not well-established at low magnetic fields, the strong 2p exciton dominates the spectrum. At higher

fields, the 2p exciton evolves into a magneto-exciton, whose dispersion is very close to LL transitions (−2 → 1, K) and (−1 → 2, K′).

By comparing to calculated spectra (see Supplementary Discussion 1B), we identify these two peaks as inter-LL transitions (−5 → 4, K) and (−4 → 5, K′) that correspond to solid arrows in Fig.

 1b. Note that our calculation does not include electron-hole interaction, so excitonic effects are not shown. Transitions indicated by dash arrows in the same energy range as in Fig. 1b are

found to be much weaker by calculations. These four transitions (two represented by solid lines and two represented by dashed lines) form a quartet that obey the conventional optical

selection rule of Δ|N | = ± 1. To the first order, the splitting between transitions (−5 → 4, K) and (−4 → 5, K′) reflects the difference between electron and hole cyclotron energies24. The

upper panel of Fig. 2b shows that the difference ΔEc scales linearly B (the red line). The lower panel of Fig. 2b shows the relative electron-hole asymmetry ΔEc /\(\bar E_c\), where

\(\overline {E_c} \) is the averaged cyclotron energy of electron and hole. The extracted ratio of ~18% is close to previously reported asymmetry in the effective masses of electrons and

holes24. Most strong transition branches behave similarly to branch C (see Supplementary Figure 4).

The only exception is branch B, where three peaks are involved as shown in Fig. 2c. A comparison with our calculation (see Supplementary Figure 3) shows that the extra peak corresponds to (0

 → 4, K′), which is forbidden by the simple selection rules Δ|N| = ± 1. This is a manifestation of the trigonal warping effect in BLG lattice22, which distorts the rotational symmetric

Mexican hat band structure to a three-fold rotational symmetric structure as shown in Fig. 1a. Therefore LLs differed by 3 m (m is an integer) are coupled by this lattice potential,

resulting in an extra ±3 ingredient to the selection rule25,26,27,28,29. This effect results in inter-LL transitions with Δ|N| = ± 4, ±2. Branch B is the only one that has three peaks with

similar energies and oscillator strength involved, owing to the close energies of LL 0 and −3 in the K′ valley as shown in Fig. 1b. Ultimately this is owing to the relatively flat band at

the top of the conduction band in K′ valley. As the LL index increases, energies of LLs are well separated so Δ|N| = ± 4 transitions are far from Δ|N| = ± 1 transitions.

Having understood the origin of branch B, we conclude that branch A corresponds to transitions (−3 → 2, K) and (−2 → 3, K′). LL transitions (−1 → 2, K) and (−2 → 1, K′) are also allowed by

selection rules and should appear at energy below that of branch A, as shown by our calculation in Fig. 2d. However, the 2p exciton peak is in the same energy range and it persists down to

zero magnetic field. As the magnetic field increases, single-particle cyclotron energy gradually surpasses the exciton binging energy. One would expect the nature of this resonance to

gradually change from exciton to inter-LL transition17,30 as the magnetic increases. At high enough magnetic field, the dispersion of this resonance can be almost described by that of the

inter-LL transition, as has been shown in ref. 30.

Next we examine the weak transitions. Figure 3a shows one branch of such transitions between strong transition branch F and G in the range of 2.4–3.0 T with a step size of 0.1 T. The spectra

are shifted to align the center of weak transition peaks to zero energy. Interestingly, three peaks emerge from the single peak at low B and their splitting widens with increased B. Figure 

3b shows a fitting of the spectrum at 2.8 T, where the sum of seven Lorentzian lineshaped peaks with an offset describes the experimental spectrum nicely. The linewidths of all peaks are ~1 

meV, indicating a very low disorder-induced inhomogeneous broadening in our samples. Figure 3c plots the energy splitting of three pairs of peaks extracted from the fitting in Fig. 3b: E78

and E89 are the splittings within branches F and G; E79 represents the splitting between the two side peaks in the branch of weak transitions. We found that E79 is within 13% of E78 + E89.

This observation, together with the fact that weak transition peaks sit at halfway between branches F and G, indicate that the three weak transitions correspond to (−9 → 7, K & K′), (−8 → 8,

K & K′), and (−7 → 9, K & K′) as labeled in Fig. 3b. Other branches of weak transitions behave similarly and more fittings can be found in Supplementary Figure 6. We therefore observed

inter-LL transitions obeying selection rules of Δ|N| = ± 2 and Δ|N| = 0. The former can be understood with the same trigonal warping effect as a combination of Δ|N| = ± 1 and ±3. But the

Δ|N| = 0 transitions goes beyond the single-particle picture in intrinsic BLG.

a Evolution of a typical branch of weak transitions between branch F and G in Fig. 1d. Spectra at 2.4–3.0 T (0.1 T step) are shifted horizontally to the center, and vertically by multiples

of 0.1. Three peaks emerge gradually in the middle. b Fitting of the 2.8 T spectrum using seven Lorentzian peaks. The two side peaks in the weak transition branch correspond to (−9 → 7) and

(−7 → 9). These Δ|N| = ± 2 transitions are allowed by the trigonal warping effect. The central peak corresponds to transition (−8 → 8) and its selection rule of Δ|N| = 0 is beyond an ideal

single-particle picture. c Energy splitting between different pairs of peaks indicated by dashed lines in b. The splitting within branch F and G are plotted in pink (E78) and green (E89),

respectively, whereas the sum of them is plotted in red (twice of the difference between electron and hole cyclotron energies). The splitting between (−9 → 7) and (−7 → 9) (E79 plotted in

blue) is very close to E78 + E89—corroborating our assignment of peaks in b.

The weak transitions acquire oscillator strengths through LL mixing by either the trigonal warping Hamiltonian or the exciton binding energy (as we will discuss later). As the magnetic field

increases, such mixing effect is weakened owing to larger LL separation—cyclotron energy. As a result, these weak transitions diminish relative to the strong transitions towards higher

magnetic fields, agreeing with our observations in Fig. 1c, d.

We further explore quantitatively the oscillator strength of weak transitions and the Δ|N| = 0 transitions relative to the integrated oscillator strength of all transitions. Owing to the

complex nature of photocurrent and photoconductivity process and their unknown quantum efficiencies, we limit our discussion to relative oscillator strengths of adjacent transition peaks and

avoid extracting the absolute optical conductivity. Figure 4a plots spectra in the vicinity of the (−5 → 5, K & K′) transition at 1.5–4.8 T. For clarity, we shifted spectra horizontally to

align the center of transition peaks to zero energy. At 1.5 T, the “weak” and “strong” transitions are of similar peak height and oscillator strength. And the only way to differentiate them

is to trace transitions from spectra at higher magnetic fields. As we increase B, the central weak transitions gradually weaken compared to strong transitions on both sides. At 4.8 T, weak

transitions drop to the noise level of the spectrum. We further explore this contrast in oscillator strength at different bandgaps. Figure 4b shows spectra in the vicinity of (−5 → 5, K &

K′) at D = 0.874, 1.02, 1.32 V/nm with the same B = 4 T. Clearly, more oscillator strength is concentrated in the central weak transitions for a bigger bandgap. We can further quantify the

oscillator strength based on fittings as shown in Fig. 3b. Figure 4c plots the oscillator strength of weak transitions (the sum of Δ|N| = 0 and ±2 in red, Δ|N | = 0 in blue) divided by the

total oscillator strength of all transitions. The data from different branches are plotted with different legends at D = 1.32 V/nm. To provide a reference for the decreasing trend of

relative oscillator strength versus cyclotron energy, we plot a solid green curve and a dashed black curve—both are inversely proportional to the squared cyclotron energy.

Relative oscillator strength of weak transitions and Δ|N| = 0 transitions.

a Photocurrent spectra in the vicinity of transition (−5 → 5) at D = 0.874 V/nm. Spectra at 1.5, 2.0, 2.5, 3.0, 3.5, 4.0, 4.5, 4.8 T are shifted horizontally to the center, and vertically by

multiples of 0.6. The contrast between the “strong” and “weak” transitions is negligible at low fields but is appreciable at high fields. At 4.8 T, weak transitions are negligible in the

spectrum. b Spectra in the vicinity of transition (−5 → 5) at D = 0.874, 1.02, and 1.32 V/nm. The relative oscillator strength of weak transitions to strong transitions increases at higher D

for the same magnetic field of 4 T. Even in the scenario where cyclotron energies are similar (1.32 V/nm & 4.5 T, and 0.874 V/nm & 4 T), the oscillator strength is still more concentrated

in weak transitions at higher D. c Relative oscillator strength as a function of the cyclotron energy at D = 1.32 V/nm. Data in blue (red) represent the ratio between the oscillator strength

in Δ|N| = 0 (Δ|N| = 0 & Δ|N| = ± 2) transition and the integrated oscillator strength. The relative oscillator strength in Δ|N| = 0 is proportional to 1/(C.E.)2, where C.E. is the cyclotron

energy. A typical energy scale of LL mixing ~1.8 meV can be extracted from this plot.

The observation of Δ|N| = 0 transitions can be qualitatively explained by a simple model in which adjacent inter-LL transitions are mixed by certain mechanism (represented by a mixing

Hamiltonian Hm). For example, we consider optical transitions (−N→N), (−N→N − 1), and (−N − 1→N) in K valley and label the excited states as |−N, N>, |−N, N − 1>, and |−N − 1, N>,

respectively. With the perturbation by Hm, the state |−N, N > changes into |−N, N>′ = |−N, N + a|−N, N − 1> +b|−N − 1, N>, where the coefficients a and b describe the mixing of excited

states. The oscillator strength of optical transition from the ground state |0> to |−N, N>′ is therefore T = |′<−N, N|H|0>|2 = |a < −N, N − 1|H|0 > + b < −N − 1, N|H|0>|2, where H is the

electron–photon coupling Hamiltonian. Here the conventional selection rule of Δ|N| = ± 1 is included as < −N, N|H|0> = 0. With finite a and b, transition (−N→N) could appear in our spectrum

as the system is excited from |0> to |−N, N > ′.

In the first order perturbation theory, the wave function mixing coefficient between two adjacent excited states |−N, N>, and |−N − 1, N> is a, b ~ Em/Ec, where Em = |<−N, N|Hm|−N−1, N>| is

the typical energy scale of the coupling Hamiltonian Hm, whereas Ec is the cyclotron energy. Therefore, the oscillator strength of Δ|N| = 0 transitions relative to the integrated oscillator

strength is |a + b|2 = 4(Em/Ec)2. Qualitatively, (Em/Ec)2 decreases as B increases as long as Em follows a sublinear relation with B, agreeing with the trend we observed in Fig. 4.

Quantitatively, the green curve in Fig. 4c indicates a wave function mixing coefficient of a~b ~ (1.8 meV/Ec) and Em ~ 1.8 meV.

With this energy scale of the mixing Hamiltonian, we discuss possible microscopic mechanisms. Beyond the single-particle picture in ideal BLG, two possible scenarios could lead to our

observation of Δ|N| = 0 transitions. In the first scenario, we include disorders but neglect the electron-hole interaction. Such disorder potential could enable the forbidden optical

transitions31. Experimentally, disordered induced unconventional inter-LL optical transition were indeed observed in GaAs quantum wells32 and Raman spectrum of graphene-like domains on

graphite33. Quantitatively31, the forbidden transition has a relative oscillator strength of ~\(\left( {{\it{\Gamma /{\rm{E}}}}_c} \right)^2\) with respect to the Δ|N| = ± 1 transition,

where Γ is the scattering induced broadening of a single LL and it plays the role of Em. Quantitatively, by setting 4(Em/Ec)2 = \(\left( {{\it{\Gamma /{\rm{E}}}}_c} \right)^{\mathrm{2}}\) to

explain our data, we can obtain Γ = 2Em ~ 3.6 meV. Fitting of the single transition peak in our spectra indicates a Lorentzian linewidth of ~1 meV, which corresponds to 4Γ (two for initial

and final states, whereas another two for the full width at half-maximum). Therefore, the realistic level broadening \({\it{\Gamma }} \sim 0.25\;\mathrm{meV}\) is more than one order of

magnitude smaller than that is needed to explain the observed oscillator strength—indicating possible mechanisms beyond the single-particle picture.

In the second scenario, we include electron-hole interaction. In a general context, interband optical excitations of a 2DEG in magnetic field is described by magneto-excitons34. Owing to

angular momentum conservation, optically active modes of magneto-excitons at k = 0 do not mix34. However, modes at finite k are results of mixing between |−N, N> and |−N−1, N> and the mixing

coefficient a (b) could be significant at klB~1. Here lB is the magnetic length defined as lB = \(\sqrt {\hbar /eB} \) ~26 nm/\(\sqrt {B(\mathrm{Tesla})} \). These modes can be optically

excited when the translational symmetry is broken by disorder or the Moire superlattice. The latter provides a momentum k\( \ge 2\pi /14 \, \mathrm{nm}\), where 14 nm is the largest period

in graphene/hBN Moire superlattice. At 1 Tesla, klB is ~10, easily satisfying the required condition for mixing |−N, N> and |−N − 1, N>. From the energy point-of-view, the Coulomb energy

between electron and hole is expected to play an important role, as the exciton binding energy (~10 meV at D = 0.874 V/nm) is big enough compared with 1.8 meV to explain the data. The

binding energy is expected to increase with enhanced displacement field and giving forbidden transitions more oscillator strength. This picture is supported by our data in Fig. 4b, where the

relative oscillator strength of Δ|N| = 0 transitions at D = 1.32 V/nm is much larger than that at D = 0.874 V/nm. Further theoretical calculation of this specific semiconductor system is

needed to corroborate this second scenario.

Our results present the first measurement of inter-LL transitions in high-quality BLG with a tunable bandgap. The many observations provide a plethora of information about interplay between

the band structure, magnetic field, and many-body interactions. There is no theoretical prediction of the behavior of Δ|N| = 0 transitions in BLG, or its dependence on the bandgap. The

selection rule of inter-LL transitions could also be used to study exotic interaction-driven ground states of materials such as the nematic phase35, where the optical selection rules can be

changed by the broken rotational symmetry. On the other hand, the (magneto-) photocurrent spectroscopy presented here could be extended to study the integer and fractional quantum Hall

states in high-quality graphene devices, shining light on the mechanism of valley and spin degeneracy lifting and excitations such as the magneto-rotons36. Beyond the quantum Hall systems,

our method is particularly suitable for studying strongly correlated electron physics in engineered van der Waals heterostructures of 2D materials37,38,39,40,41, e.g., the Moire gap and

correlated insulating gaps in twisted graphene/TMD systems as well as ABC-stacked trilayer graphene/hBN superlattices.

The basic sample configuration and measurement scheme are similar to that in ref. 1. The data at magnetic fields higher than 5 Tesla was collected at the SCM3 in the National High Magnetic

Field Laboratory (NHMFL). Sample temperature in the Cornell setup is ~10 K, whereas that in the NHMFL setup is ~5 K owing to the helium exchange gas environment. A Bruker IFS 66 v

spectrometer was used at NHMFL. The illumination is focused by a ZnSe lens located outside (inside) the cryostat at Cornell (NHMFL). The thickness of hBN flakes is ~20 nm for both top and

bottom gates. All the measurement is done at a total charge density of zero to minimize the dark current noise. The displacement field is calculated by using a dielectric constant of boron

nitride of 3.9. The data in this manuscript are mostly based on two separate high-quality graphene samples.

The data that support the findings of this study are available from the corresponding author upon request.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

References Ju, L. et al. Tunable excitons in bilayer graphene. Science 358, 907 (2017).

Article  ADS  CAS  Google Scholar 

Landau, L. D. & Lifshitz, E. Course of Theoretical Physics Vol 3 Quantum Mechanics. (Pergamon Press, 1958).

Zhang, Y., Tan, Y.-W., Stormer, H. L. & Kim, P. Experimental observation of the quantum Hall effect and Berry’s phase in graphene. Nature 438, 201 (2005).

Article  ADS  CAS  Google Scholar 

Novoselov, K. S. et al. Two-dimensional gas of massless Dirac fermions in graphene. Nature 438, 197 (2005).

Article  ADS  CAS  Google Scholar 

Novoselov, K. S. et al. Unconventional quantum Hall effect and Berry’s phase of 2π in bilayer graphene. Nat. Phys. 2, 177 (2006).

Article  Google Scholar 

Young, A. F. et al. Spin and valley quantum Hall ferromagnetism in graphene. Nat. Phys. 8, 550 (2012).

Article  CAS  Google Scholar 

Hunt, B. et al. Massive Dirac Fermions and Hofstadter Butterfly in a van der Waals Heterostructure. Science 340, 1427 (2013).

Article  ADS  CAS  Google Scholar 

Miller, D. L. et al. Observing the quantization of zero mass carriers in graphene. Science 324, 924–927 (2009).

Article  ADS  CAS  Google Scholar 

Feldman, B. E., Krauss, B., Smet, J. H. & Yacoby, A. Unconventional sequence of fractional quantum Hall states in suspended graphene. Science 337, 1196–1199 (2012).

Article  ADS  CAS  Google Scholar 

Jiang, Z. et al. Infrared spectroscopy of Landau levels of graphene. Phys. Rev. Lett. 98, 197403 (2007).

Article  ADS  CAS  Google Scholar 

Deacon, R., Chuang, K.-C., Nicholas, R., Novoselov, K. & Geim, A. Cyclotron resonance study of the electron and hole velocity in graphene monolayers. Phys. Rev. B 76, 081406 (2007).

Article  ADS  Google Scholar 

Chen, Z.-G. et al. Observation of an intrinsic bandgap and Landau level renormalization in graphene/boron-nitride heterostructures. Nat. Commun. 5, 4461 (2014).

Article  ADS  CAS  Google Scholar 

Russell, B. J., Zhou, B., Taniguchi, T., Watanabe, K. & Henriksen, E. A. Many-particle effects in the cyclotron resonance of encapsulated monolayer graphene. Phys. Rev. Lett. 120, 047401

(2018).

Article  ADS  CAS  Google Scholar 

Henriksen, E. A. et al. Cyclotron resonance in bilayer graphene. Phys. Rev. Lett. 100, 087403 (2008).

Article  ADS  CAS  Google Scholar 

Sadowski, M. L., Martinez, G., Potemski, M., Berger, C. & de Heer, W. A. Landau level spectroscopy of ultrathin graphite layers. Phys. Rev. Lett. 97, 266405 (2006).

Article  ADS  CAS  Google Scholar 

Orlita, M. et al. Magneto-optics of bilayer inclusions in multilayered epitaxial graphene on the carbon face of SiC. Phys. Rev. B 83, 125302 (2011).

Article  ADS  Google Scholar 

Wang, Z., Shan, J. & Mak, K. F. Valley-and spin-polarized Landau levels in monolayer WSe 2. Nat. Nanotechnol. 12, 144 (2017).

Article  ADS  CAS  Google Scholar 

Schafgans, A. et al. Landau level spectroscopy of surface states in the topological insulator Bi 0.91 Sb 0.09 via magneto-optics. Phys. Rev. B 85, 195440 (2012).

Article  ADS  Google Scholar 

Orlita, M. et al. Observation of three-dimensional massless Kane fermions in a zinc-blende crystal. Nat. Phys. 10, 233 (2014).

Article  CAS  Google Scholar 

Varlet, A. et al. Anomalous sequence of quantum hall liquids revealing a tunable lifshitz transition in bilayer graphene. Phys. Rev. Lett. 113, 116602 (2014).

Article  ADS  Google Scholar 

McCann, E. & Fal’ko, V. I. Landau-level degeneracy and quantum hall effect in a graphite bilayer. Phys. Rev. Lett. 96, 086805 (2006).

Article  ADS  Google Scholar 

Dresselhaus, M. S. & Dresselhaus, G. Intercalation compounds of graphite. Adv. Phys. 51, 1–186 (2002).

Article  ADS  CAS  Google Scholar 

Gusynin, V., Sharapov, S. & Carbotte, J. J. Anomalous absorption line in the magneto-optical response of graphene. Phys. Rev. Lett. 98, 157402 (2007).

Article  ADS  CAS  Google Scholar 

Zou, K., Hong, X. & Zhu, J. Effective mass of electrons and holes in bilayer graphene: electron-hole asymmetry and electron-electron interaction. Phys. Rev. B 84, 085408 (2011).

Article  ADS  Google Scholar 

Cao, T., Wu, M. & Louie, S. G. Unifying optical selection rules for excitons in two dimensions: band topology and winding numbers. Phys. Rev. Lett. 120, 087402 (2018).

Article  ADS  CAS  Google Scholar 

Zhang, X., Shan, W.-Y. & Xiao, D. Optical selection rule of excitons in gapped chiral fermion systems. Phys. Rev. Lett. 120, 077401 (2018).

Article  ADS  CAS  Google Scholar 

Abergel, D. & Fal’ko, V. I. Optical and magneto-optical far-infrared properties of bilayer graphene. Phys. Rev. B 75, 155430 (2007).

Article  ADS  Google Scholar 

Dresselhaus, M. S. & Mavroides, J. G. Energy band parameter determination in graphite. Carbon 3, 465–473 (1966).

Article  CAS  Google Scholar 

Doezema, R., Datars, W., Schaber, H. & Van Schyndel, A. Far-infrared magnetospectroscopy of the Landau-level structure in graphite. Phys. Rev. B 19, 4224 (1979).

Article  ADS  CAS  Google Scholar 

Stier, A. V. et al. Magnetooptics of exciton Rydberg states in a monolayer semiconductor. Phys. Rev. Lett. 120, 057405 (2018).

Article  ADS  CAS  Google Scholar 

Briskot, U., Dmitriev, I. & Mirlin, A. J. Quantum magneto-oscillations in the ac conductivity of disordered graphene. Phys. Rev. B 87, 195432 (2013).

Article  ADS  Google Scholar 

Lyo, S., Jones, E. & Klem, J. Breaking of the usual selection rule for magnetoluminescence in doped semiconductor quantum wells. Phys. Rev. Lett. 61, 2265 (1988).

Article  ADS  CAS  Google Scholar 

Kühne, M. et al. Polarization-resolved magneto-Raman scattering of graphenelike domains on natural graphite. Phys. Rev. B 85, 195406 (2012).

Article  ADS  Google Scholar 

Sári, J. & Tőke, C. Theory of inter-Landau-level magnetoexcitons in bilayer graphene. Phys. Rev. B 87, 085432 (2013).

Article  ADS  Google Scholar 

Fradkin, E., Kivelson, S. A., Lawler, M. J., Eisenstein, J. P. & Mackenzie, A. P. Nematic Fermi fluids in condensed matter physics. Annu. Rev. Condens. Matter Phys. 1, 153–178 (2010).

Article  ADS  CAS  Google Scholar 

Girvin, S., MacDonald, A. & Platzman, P. Magneto-roton theory of collective excitations in the fractional quantum Hall effect. Phys. Rev. B 33, 2481 (1986).

Article  ADS  CAS  Google Scholar 

Cao, Y. et al. Correlated insulator behaviour at half-filling in magic-angle graphene superlattices. Nature 556, 80 (2018).

Article  ADS  CAS  Google Scholar 

Chen, G. et al. Evidence of a gate-tunable Mott insulator in a trilayer graphene moiré superlattice. Nature 15, 237 (2019).

CAS  Google Scholar 

Liu, X. et al. Spin-polarized correlated insulator and superconductor in twisted double bilayer graphene. Preprints at https://arxiv.org/abs/1903.08130

Regan, E. C. et al. Mott and generalized Wigner crystal states in WSe2/WS2 moiré superlattices. Nature 579, 359–363 (2020).

Article  ADS  CAS  Google Scholar 

Tang, Y. et al. Simulation of Hubbard model phsics in WSe2/WS2 moiré superlattices. Nature 579, 353–358 (2020).

Article  ADS  CAS  Google Scholar 

Download references

L.J. and L.W. acknowledge support from a Kavli Postdoctoral Fellowship. The work at Cornell was supported in part by the Cornell Center for Materials Research with funding from the NSF MRSEC

program (DMR-1719875), by the Nanoelectronics Research Corporation (NERC), a wholly owned subsidiary of the Semiconductor Research Corporation (SRC), through the Institute for

Nanoelectronics Discovery and Exploration (INDEX), by the Air Force Office of Scientific Research (MURI: FA9550-16-1-0031), and by the Office of Naval Research under grant number

N00014−12-1-0072. Sample fabrication was performed at the Cornell Nanoscale Science & Technology Facility, a member of the National Nanotechnology Infrastructure Network, which is supported

by the National Science Foundation (ECCS-1542081). X.L. acknowledges the support from City University of Hong Kong (Project no. 9610428), and also thanks Erwin Schrodinger International

Institute for Mathematics and Physics of the University of Vienna for its hospitality during the finalization of this work. F.Z. is supported by ARO under Grant No. W911NF-18-1-0416 and NSF

under Grant No. DMR-1921581 (DMREF)..

Author informationAuthor notes Long Ju

Present address: Department of Physics, Massachusetts Institute of Technology, Cambridge, MA, 02139, USA

These authors contributed equally: Long Ju, Lei Wang.

Authors and Affiliations Kavli Institute at Cornell for Nanoscale Science, Ithaca, NY, 14853, USA

Long Ju, Lei Wang & Paul L. McEuen

Laboratory of Atomic and Solid State Physics, Cornell University, Ithaca, NY, 14853, USA

Long Ju, Lei Wang, Erich Mueller & Paul L. McEuen

Department of Physics, City University of Hong Kong, Kowloon, Hong Kong SAR

Xiao Li

National High Magnetic Field Laboratory, Tallahassee, FL, 32312, USA

Seongphill Moon, Mike Ozerov, Zhengguang Lu & Dmitry Smirnov

National Institute for Materials Science, 1-1 Namiki, Tsukuba, 305-0044, Japan

Takashi Taniguchi & Kenji Watanabe

Department of Physics, University of Texas at Dallas, Richardson, TX, 75080, USA

Fan Zhang

School of Electrical and Computer Engineering, Cornell University, Ithaca, NY, 14853, USA

Farhan Rana

AuthorsLong JuView author publications You can also search for this author inPubMed Google Scholar

Lei WangView author publications You can also search for this author inPubMed Google Scholar

Xiao LiView author publications You can also search for this author inPubMed Google Scholar

Seongphill MoonView author publications You can also search for this author inPubMed Google Scholar

Mike OzerovView author publications You can also search for this author inPubMed Google Scholar

Zhengguang LuView author publications You can also search for this author inPubMed Google Scholar

Takashi TaniguchiView author publications You can also search for this author inPubMed Google Scholar

Kenji WatanabeView author publications You can also search for this author inPubMed Google Scholar

Erich MuellerView author publications You can also search for this author inPubMed Google Scholar

Fan ZhangView author publications You can also search for this author inPubMed Google Scholar

Dmitry SmirnovView author publications You can also search for this author inPubMed Google Scholar

Farhan RanaView author publications You can also search for this author inPubMed Google Scholar

Paul L. McEuenView author publications You can also search for this author inPubMed Google Scholar

L.J. conceived the experiment and performed the photocurrent measurement with the help from S.M., M.O., Z.L. and D.S. L.W. fabricated the device. X.L., E.M., and F.Z. contributed to the

theoretical calculation. T.K. and K.W. grew the BN crystal. L.J. wrote the manuscript with input from all. P.M. supervised the whole project.

Corresponding author Correspondence to Paul L. McEuen.

The authors declare no competing interests.

Peer review information Nature Communications thanks Sergei Sharapov and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

InformationRights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or

format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or

other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in

the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the

copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this articleCite this article Ju, L., Wang, L., Li, X. et al. Unconventional valley-dependent optical selection rules and landau level mixing in bilayer graphene. Nat Commun 11, 2941

(2020). https://doi.org/10.1038/s41467-020-16844-y

Download citation

Received: 30 December 2019

Accepted: 01 April 2020

Published: 10 June 2020

DOI: https://doi.org/10.1038/s41467-020-16844-y

Share this article Anyone you share the following link with will be able to read this content:

Get shareable link Sorry, a shareable link is not currently available for this article.

Copy to clipboard Provided by the Springer Nature SharedIt content-sharing initiative

This article is cited by Physics of 2D Materials for Developing Smart Devices Neeraj GoelRahul Kumar Nano-Micro Letters (2025)

Engineering symmetry breaking in 2D layered materials Luojun DuTawfique HasanZhipei Sun Nature Reviews Physics (2021)