Abstract
Under strong laser fields, electrons in solids radiate high-harmonic fields by travelling through quantum pathways in Bloch bands in the sub-laser-cycle timescales. Understanding these pathways in the momentum space through the high-harmonic radiation can enable an all-optical ultrafast probe to observe coherent lightwave-driven processes and measure electronic structures as recently demonstrated for semiconductors. However, such demonstration has been largely limited for semimetals because the absence of the bandgap hinders an experimental characterization of the exact pathways. In this study, by combining electrostatic control of chemical potentials with HHG measurement, we resolve quantum pathways of massless Dirac fermions in graphene under strong laser fields. Electrical modulation of HHG reveals quantum interference between the multi-photon interband excitation channels. As the light-matter interaction deviates beyond the perturbative regime, elliptically polarized laser fields efficiently drive massless Dirac fermions via an intricate coupling between the interband and intraband transitions, which is corroborated by our theoretical calculations. Our findings pave the way for strong-laser-field tomography of Dirac electrons in various quantum semimetals and their ultrafast electronics with a gate control.
Similar content being viewed by others
Introduction
Strong laser fields give rise to extreme nonlinear optical phenomena in atoms â high harmonic generation (HHG)1,2. In the gas phase, lightwave-driven (or laser-field-driven) electrons radiate intense high-order harmonics through the three-step re-collision processes of tunnel-ionization, free acceleration, and recombination3,4, which has offered the fundamental basis for attosecond science and technology for the last decades5. More recently, since the first discovery of solid-state HHG in a ZnO crystal6, high harmonic radiation in solids has emerged as a fascinating new route to probe electronic states in sub-laser-cycle timescale7,8,9,10,11,12,13,14,15,16,17,18 and to realize compact, coherent light sources with high photon energy (up to ~40âeV)8. In particular, the lightwave-driven quantum pathways of electrons in the momentum space can be visualized by analyzing the spectrum and polarization of high harmonics, which allows to reconstruct band structures9,10, measure electronic properties such as Berry curvature11,12,13 and Chern number14,15, and observe coherent quasi-particle dynamics in sub-cycle timescale16,17,18. However, unlike the gas phase, delocalized electrons in solids can travel complicated multiple quantum pathways via interband and intraband transitions7,19,20,21, which sensitively depends on chemical potential and bandgap size under a given laser excitation condition22. In semiconductors with a large bandgap, the complication can be effectively addressed by high-order sideband generation10,16,17,18 where one laser pulse prepares coherent interband excitation of electron-hole pairs (or excitons) and another laser pulse at terahertz frequencies strongly drives their intraband motion. Unfortunately, it is challenging to avoid the complication for semimetals and narrow bandgap semiconductors with high-order sideband generation because the charge carriers can be readily excited via the interband transition across a small bandgap under strong laser fields even at terahertz frequencies13,23,24,25.
Graphene is a two-dimensional semimetal (or a zero-bandgap semiconductor) with linear energy dispersion hosting massless Dirac fermions26. Electrons in graphene exhibit strong light-matter interaction due to its unusual band structure. As a linear optical response, the interband and intraband transitions lead to a high absorption coefficient over broad spectral range from terahertz to ultraviolet frequencies27,28. The optical resonances from the interband transition enhance nonlinear optical responses29,30,31,32, including third harmonic generation and four-wave mixing. As laser intensity increases over ~1âGWâcmâ2, collective oscillation of massless Dirac fermions via the intraband transition enables remarkably efficient HHG at terahertz frequencies33,34. According to a recent experimental study under extreme laser intensity (~ 1 TW cmâ2)23, elliptically polarized excitation can further enhance HHG via efficient pair generation mechanism from light-induced transient band structure modification, which is not possible in the standard re-collision model for the gas phases. Theoretical studies suggest other intriguing pathways available in graphene for HHG processes22,35,36,37,38,39,40,41,42. On the other hand, the interband and intraband transitions in graphene can be drastically controlled by electrostatically injecting carriers28,43,44 because the two-dimensional linear band with the high Fermi velocity (~ 106âmâsâ1) allows large modification of the chemical potential by small amount of carrier injection. Therefore, graphene provides an ideal platform to resolve the interplay between the interband and intraband transitions in quantum pathways of carriers under strong laser fields.
In this work, we investigate widely tunable HHG in graphene with electrostatic control of the chemical potential to resolve quantum pathways of massless Dirac fermions under strong laser fields. Under linearly-polarized excitation with laser intensity of ~1âGWâcmâ2, harmonic intensity exhibits non-monotonic dependence on the chemical potential, which arises from the destructively interfering excitation channels via the multi-photon interband transitions. Under elliptically polarized excitation, as the light-matter interaction deviates beyond the perturbative regime, charge carriers excited by the multi-photon transitions efficiently radiate the high harmonics with the drastically rotated polarization axis, which is completely different from the previous study23. Interestingly, the helical response is substantially suppressed by blocking the interband transitions. Our full quantum mechanical theory identifies a characteristic mechanism for massless Dirac fermions where the intricate coupling between interband and intraband transitions efficiently drives charge carrier motion under elliptically polarized excitation. Our study paves the way for visualizing sub-cycle dynamics of massless charge carriers in quantum materials including graphene, topological insulators and Weyl semimetals, and realizing their ultrafast electronics and attosecond photonics with an electrostatic control.
Results and discussion
HHG measurement with a gate control of chemical potentials
Figure 1aâc schematically illustrate the electronic processes in graphene under strong laser fields. In the real space (Fig. 1a), laser field (E(t)) centered at the frequency of Ï excites and transports charge carriers to generate an ultrafast anharmonic current (J(t)) which radiates high harmonics (I(nÏ))7,45, where n is a positive integer to represent harmonic orders. In the momentum space (Fig. 1b, c), J(t) can be generated from charge carriers following the multiple quantum pathways which are composed of the interband and intraband transitions. The control of the chemical potential (μ) leads to a profound impact on the individual quantum pathway. For the charge-neutral case (μâ=â0) in Fig. 1b, where μ is at the Dirac point or the charge-neutral point, charge carriers in the pathways can be excited and recombined via linear and nonlinear optical transitions between the valence band and the conduction band (black solid arrow). The photo-excited carriers can also be driven via the intraband transition simultaneously (red and blue solid arrows). For the hole-doped case (μâ<â0) in Fig. 1c, the interband transition (black dashed arrow) is forbidden for the empty states between the Dirac point and μ, which does not create photo-excited carriers28,43,44. Subsequently, the laser field cannot drive charge carriers via the connected intraband transitions (red and blue dashed arrows). For the electron-doped case (μâ>â0), the same pathways are not available either since the interband transition is blocked for the filled states between the Dirac point and μ. On the other hand, the laser field can drive electrostatically injected carriers to generate J(t) solely via intraband transition without interband transition34,46. Therefore, control of μ allows to understand distinct contributions of different quantum pathways to the HHG process in graphene.
a Schematic of HHG measurement. An intense mid-infrared femtosecond laser pulse \({{{{{\bf{E}}}}}}(t)\) generates an anharmonic ultrafast current \({{{{{\bf{J}}}}}}(t)\) in the graphene device with ion-gel gating, radiating high harmonics \({{{{{{\bf{I}}}}}}}^{{{{{{\boldsymbol{(}}}}}}{{{{{\boldsymbol{n}}}}}}{{{{{\boldsymbol{\omega }}}}}}{{{{{\boldsymbol{)}}}}}}}\). b, c Schematic of chemical potential (\(\mu\)) dependent HHG process in momentum space. For the charge-neutral case (b), \({{{{{\bf{J}}}}}}(t)\) is generated simultaneously by the interband transition (black solid arrow), intraband transition along the x-direction (red solid arrow), and intraband transition along the y-direction (blue solid arrow). For the highly doped case, the interband transition (black dashed arrow) and connected intraband transitions (red and blue dashed lines) are blocked due to Pauli blocking. d HHG spectrum. The fifth and seventh harmonics are observed at 1.38 and 1.93âeV, respectively, under excitation photon energy of 0.28âeV. e Intensity of the fifth harmonics as a function of the laser peak intensity \({I}_{{{{{{\rm{exc}}}}}}}\). The harmonic intensity (black circles) initially scales with the fifth power of \({I}_{{{{{{\rm{exc}}}}}}}\) (blue line), but it gradually saturates to the second power (red line) of \({I}_{{{{{{\rm{exc}}}}}}}\) as increasing \({I}_{{{{{{\rm{exc}}}}}}}\). Black arrow indicates the intensity \({I}_{{{{{{\rm{exc}}}}}}}\) of 3.1âGW/cm2. f Laser polarization profile with ellipticity \({\varepsilon }_{{{{{{\rm{exc}}}}}}}\) of 0.3. g Normalized polarization profile of the fifth harmonics under elliptically polarized excitation (f) with different \({I}_{{{{{{\rm{exc}}}}}}}\). h Polarization profile of the fifth harmonics with gate voltage control \({V}_{{{{{{\rm{G}}}}}}}\) under \({\varepsilon }_{{{{{{\rm{exc}}}}}}}\)â=â0.3. For clarity, polarization plot for \({V}_{{{{{{\rm{G}}}}}}}\)â=ââ1.5âV is multiplied by a factor of 5. As \({V}_{{{{{{\rm{G}}}}}}}\) decreases from 0.9âV to â1.5, the harmonic intensity is reduced by ~50 times while the polarization axis is rotated by 52.7°.
The experimental configurations are schematically described in Fig. 1a. Graphene grown by chemical vapor deposition method (CVD) is transferred on a sapphire substrate. The device structure of a field-effect transistor is fabricated on graphene to control and characterize μ. The source electrode (S) and the drain electrode (D) measure static current along the graphene channel under the bias voltage (Vb). Ion-gel covers the gate electrode (G) and graphene to electrostatically inject electrons or holes into graphene by applying the gate voltages (VG). Intense mid-infrared femtosecond laser pulses (E(t)) are irradiated on graphene with elliptical polarization. The polarization major axis (the x-direction) and the minor axis (the y-direction) are assigned parallel and perpendicular, respectively, to the direction along the source and the drain electrodes of the devices. The ratio of two orthogonal electric fields (Ey/Ex), defined as ellipticity (εexc), is controlled by a liquid crystal variable retarder. Combination of the polarization analysis optics, the spectrometer, and EMCCD measures and analyzes detailed characteristics of transmitted harmonics (I(nÏ)) including their spectra, intensity, and polarization. More detailed configuration is provided in Supplementary Note.
High harmonic spectrum (Fig. 1d) shows the fifth and seventh harmonics at 1.38âeV and 1.93âeV, respectively, under laser excitation at 0.28âeV (see Supplementary Fig. 1 for subtraction of background luminescence). Laser peak intensity (Iexc) dependence indicates that the non-perturbative response emerges in generation of the harmonic signals under strong laser fields. Under εexcâ=â0 (i.e., linearly-polarized along the x-direction), intensity of the fifth harmonics (black circles in Fig. 1e) scales initially to the fifth power of Iexc (blue line), which is expected from the perturbative response47, but saturating roughly to the second power of Iexc (red line) as Iexc increases from 1.8 to 55âGWâcmâ2. In addition, elliptically polarized excitation with εexcâ=â0.3 reveals unusual polarization profile of harmonics. Figure 1g displays a polar plot of the fifth harmonic intensity which is measured after a linear polarizer as a function of the polarizer angle (Ï). The orientation of the sample and the major axis of the laser polarization are fixed along the x-direction (Ïâ=â0°) as described in Fig. 1a, f, respectively. The polarization axis of the harmonic signals in Fig. 1g is rotated counterclockwise by 24.9° for 2.4âGWâcmâ2 in comparison to the laser polarization axis (Fig. 1f). Higher Iexc rotates the polarization axis of the harmonic signals even further, which cannot be expected from the perturbative response47(see Supplementary Note for Iexc - dependence of the fifth harmonic expected from perturbative nonlinear optics). The polar plots of harmonics intensity and laser intensity are normalized to clearly visualize the rotation angles. Despite low Iexc of a few GWâcmâ2, we observe similar behavior of the fifth harmonics to the experimental results in the previous study23 with high Iexc of ~1âTWâcmâ2. Surprisingly, however, we observe that application of \({V}_{{{{{{\rm{G}}}}}}}\) drastically modifies the fifth harmonics (Fig. 1h). As \({V}_{{{{{{\rm{G}}}}}}}\) is varied from 0.9âV to â1.5âV for the same \({I}_{{{{{{\rm{exc}}}}}}}\) of 3.1âGWâcmâ2 (marked with black arrow in Fig. 1e), harmonic intensity is reduced by ~50 times while polarization axis is rotated from 52.7° to nearly 0°. The widely tunable harmonic signals with \({V}_{{{{{{\rm{G}}}}}}}\) suggest that charge carriers indeed undergo different pathways of harmonic generation depending on \(\mu\).
In order to calibrate our devices from \({V}_{{{{{{\rm{G}}}}}}}\) to \(\mu\), we measure both electrical transport and infrared transmittance spectra of graphene. Blue line in Fig. 2a indicates that the maximum resistance is achieved at \({V}_{{{{{{\rm{G}}}}}}}\)â=â0.95âV, identifying as the charge-neutral point (\(\mu\)â=â0âeV). Resistance sharply drops for \({V}_{{{{{{\rm{G}}}}}}}\)â<â0.95âV (or \({V}_{{{{{{\rm{G}}}}}}}\)â>â0.95âV) due to carrier injection of holes (or electrons). \(\mu\) away from the Dirac point can be identified by infrared transmission spectroscopy which probes modulation of interband optical transition28,43,44. In the linear optical process, interband transition is forbidden by Pauli blocking for photon energy below \(2\left|\mu \right|\). Real part of optical conductivity shows a step-like absorption edge at photon energy of \(2\left|\mu \right|\), and Fig. 2b demonstrates such step-like profiles with broadening from scattering processes. The absorption edge gradually moves to higher energy as \(\mu\) is lowered from the Dirac point. Fitting to the model based on Kubo formula provides \(\mu\) which are summarized with black empty squares in Fig. 2a. \(\mu\) measured from the electrical and optical measurements shows good agreement for ion-gel gating of graphene with a capacitance 4.8âμFâcmâ2 (see Supplementary Note for more details on the experimental determination of \(\mu\)). In this study, we focus on the hole-doped region (\(\mu\) < 0) which provides wider electro-chemical window for our graphene devices without chemical degradation by ion-gel. Nevertheless, the physical mechanism should be identical for the electron-doped region (\(\mu\) > 0) as discussed above.
a Electrical transport measurement of the graphene channel. Resistance (blue solid line) shows the maximum at the gate voltage \({V}_{{{{{{\rm{G}}}}}}}\) of 0.95âV, identifying the charge-neutral case (\(\mu\)â=â0âeV). b Infrared transmission measurement of the graphene device. As \({V}_{{{{{{\rm{G}}}}}}}\) decreases from 0.1âV to â1.3âV, a step-like absorption edge moves to higher energy. Fitting to the model based on Kubo formula (gray solid lines) determines \(\mu\) which is summarized in a with black squares. câe Harmonic intensity as a function of \(2\left|\mu \right|\) (or \(2\left|\mu \right|/{E}_{{{{{{\rm{ph}}}}}}}\) where \({E}_{{{{{{\rm{ph}}}}}}}\) is laser photon energy). \({I}_{{{{{{\rm{x}}}}}}}^{(5\omega )}\) (red circles) and \({I}_{{{{{{\rm{y}}}}}}}^{(5\omega )}\) (blue circles) represent harmonic intensity along the x-direction and y-direction, respectively. \({I}_{{{{{{\rm{x}}}}}}}^{(5\omega )}\)â+â\({I}_{{{{{{\rm{y}}}}}}}^{(5\omega )}\) (black squares) represents total intensity of the harmonics. Under linearly-polarized excitation (c), \({I}_{{{{{{\rm{x}}}}}}}^{(5\omega )}\) shows a resonance-like profile with the maximum intensity as \(2\left|\mu \right|\)â=â0.87âeV while \({I}_{{{{{{\rm{y}}}}}}}^{(5\omega )}\) is absent beyond the noise level. Under elliptically polarized excitation with ellipticity \({\varepsilon }_{{{{{{\rm{exc}}}}}}}\)â=â0.2 (d) and \({\varepsilon }_{{{{{{\rm{exc}}}}}}}\)â=â0.4 (e), the fifth harmonics shows \(I\)x(5Ï) and \(I\)y(5Ï) with a resonance-like profile and a step-like profile, respectively.
Interband transition channels of HHG
As the simplest case, we examine harmonic generation in graphene under linearly-polarized laser excitation (\({\varepsilon }_{{{{{{\rm{exc}}}}}}}\) = 0). \({I}_{{{{{{\rm{exc}}}}}}}\) is kept at 3.1âGWâcmâ2 to avoid any possible laser-induced damage during a long-term exposure for accurate optical measurement (see Supplementary Fig. 2 for the laser-induced damage on graphene). In Fig. 2c, the fifth harmonics are observed for the polarization along the x-direction (red circles labeled with \({I}_{{{{{{\rm{x}}}}}}}^{(5\omega )}\)) with the non-monotonic dependence on the chemical potential. For the y-direction, no signal is observed beyond the noise level (blue circles labeled with \({I}_{{{{{{\rm{y}}}}}}}^{(5\omega )}\)). We find that the fifth harmonics are linearly-polarized along the parallel direction of the laser polarization regardless of the sample orientation (see Supplementary Fig. 3 for the harmonic signals from different sample orientations). As the chemical potential is lowered, \({I}_{{{{{{\rm{x}}}}}}}^{(5\omega )}\) shows a resonance-like profile with the maximum intensity at \(2\left|\mu \right|\)â=â0.87âeV. Dramatically small signal for \(2\left|\mu \right|\)â>â1.5âeV suggests that the intraband transition of electrostatically injected carriers alone provide minor contribution to HHG in highly doped graphene under \({I}_{{{{{{\rm{exc}}}}}}}\) of 3.1âGWâcmâ2. Instead, the interband carrier excitation plays an important role in the pathway to generate the fifth harmonics.
Enhancement of \({I}_{{{{{{\rm{x}}}}}}}^{(5\omega )}\) at the specific \(\mu\) indicates that the interband excitation of charge carriers is dominated by the multi-photon transition. The multi-photon transition rate is determined by the multiple resonant channels with the energy separation of \({{mE}}_{{{{{{\rm{ph}}}}}}}\), where \(m\) and \({E}_{{{{{{\rm{ph}}}}}}}\) denote the number of photons involved in the transition and the individual photon energy of laser, respectively47. According to the recent studies on the perturbative nonlinear optical process in graphene29,40,48, these resonant channels can destructively interfere each other while their magnitudes of transition rate are in the similar order of magnitude. As the resonant channels are sequentially switched off by Pauli blocking, the harmonic signal can increase as the destructive interference is partly eliminated. The harmonic signal eventually disappears as all the relevant resonant channels are switched off. Therefore, enhancement of \({I}_{{{{{{\rm{x}}}}}}}^{(5\omega )}\) in Fig. 2c can be explained by partial elimination of destructive interference among available interband transition channels via multi-photon transition40.
Our interpretation can be further supported by measurement with the different laser excitation energies. For the laser excitation energy at 0.31âeV and 0.35âeV, \({I}_{{{{{{\rm{x}}}}}}}^{(5\omega )}\) shows similar resonance-like profiles, but the entire profiles are precisely shifted with the laser excitation energies (see Supplementary Fig. 4 for the excitation energy dependence). In principle, harmonics enhanced by multi-photon transition exhibits series of multiple sharp resonance-like profiles at \(2\left|\mu \right|/{E}_{{{{{{\rm{ph}}}}}}}\)â=ân, which is expected by theoretical calculation in the perturbative regime40. However, we observe a broad resonance-like profile (Fig.2c) with the maximum intensity at \(2\left|\mu \right|/{E}_{{{{{{\rm{ph}}}}}}}\) slightly above 3. This can be due to ultrafast electronic scattering process which drastically broaden and merge series of sharp resonance-like profiles together. (see Supplementary Fig. 5 for the variation of \(2\left|\mu \right|\) for the maximum \({I}_{{{{{{\rm{x}}}}}}}^{(5\omega )}\) over devices and Supplementary Note for the theoretical calculation of \({I}_{{{{{{\rm{x}}}}}}}^{(5\omega )}\) as a function of \(2\left|\mu \right|\) and discussion on the detailed characteristics).
HHG under elliptically polarized excitation
Harmonic generation in graphene is further examined under elliptically polarized laser excitation. Elliptically polarized excitation allows to probe the interplay between microscopic processes in the generation of harmonics. According to the standard re-collision model for atomic gases3, rotation of electric field steers tunnel-ionized electrons away from parent ions during the acceleration process, which decreases probability of the recombination process. Harmonic generation efficiency consequently drops with laser ellipticity49. For solids with crystalline symmetry, elliptically polarized excitation can enhance HHG efficiency as recently demonstrated for MgO, Si, Bi2Se3, and graphene19,20,23,50. Theoretical studies suggest several possible models where the enhancement is due to the increased recombination probability at the neighboring atomic sites20 or due to the tightly coupled dynamics of the intraband and interband transitions22,51. Alternatively, strong laser fields can modify a band structure of solids from a semiconducting structure to a semimetallic structure, which facilitates ultrafast current generation under elliptically polarized excitation23,52.
Figure 2d and e show the fifth harmonic signals in graphene with laser ellipticity \({\varepsilon }_{{{{{{\rm{exc}}}}}}}\)= 0.2 and 0.4, respectively. While \({\varepsilon }_{{{{{{\rm{exc}}}}}}}\) is modified, the pulse energy is set identically to the case of linearly-polarized excitation. The harmonic signals are observed for both polarization directions along the x-direction (red circles) and the y-direction (blue circles). For the charge-neutral case (\(2\left|\mu \right|\)â=â0âeV), total intensity of the harmonic signal (black squares) is noticeably increased under the laser excitation of \({\varepsilon }_{{{{{{\rm{exc}}}}}}}\,\)= 0.2 in comparison to \({\varepsilon }_{{{{{{\rm{exc}}}}}}}\) = 0. This enhancement can be attributed to efficient generation of harmonic signals for \({I}_{{{{{{\rm{y}}}}}}}^{(5\omega )}\). Despite small laser fields along the y-direction, elliptically polarized excitation generates surprisingly strong \({I}_{{{{{{\rm{y}}}}}}}^{(5\omega )}\). For example, \({I}_{{{{{{\rm{y}}}}}}}^{(5\omega )}\) reaches nearly 60 percent of \({I}_{{{{{{\rm{x}}}}}}}^{(5\omega )}\) despite 25 times smaller laser intensity along the minor axis than the major axis. Remarkably, \({I}_{{{{{{\rm{y}}}}}}}^{(5\omega )}\) is even greater than \({I}_{{{{{{\rm{x}}}}}}}^{(5\omega )}\) for laser excitation of \({\varepsilon }_{{{{{{\rm{exc}}}}}}}\) = 0.4. We do not observe any noticeably different behavior for other sample orientations (see supplementary Fig. 6 for the sample orientation dependence of harmonic signals for elliptically polarized excitation). Efficient generation of \({I}_{{{{{{\rm{y}}}}}}}^{(5\omega )}\) under elliptically polarized excitation suggests an intriguing pathway in graphene where laser fields \({E}_{{{{{{\rm{x}}}}}}}\) and \({E}_{{{{{{\rm{y}}}}}}}\) are strongly coupled for the generation of harmonics.
The chemical potential dependence shows that this pathway can be switched off with the \(\mu\) control. \({I}_{{{{{{\rm{y}}}}}}}^{(5\omega )}\) disappears a step-like profile where intensity decreases roughly across \(2\left|\mu \right|\)â=â0.87âeV for both \({\varepsilon }_{{{{{{\rm{exc}}}}}}}\) = 0.2 and 0.4. On the other hand, \({I}_{{{{{{\rm{x}}}}}}}^{(5\omega )}\) shows a resonance-like profile with maximum at \(2\left|\mu \right|\,\)=â0.87âeV as in the case of linearly-polarized excitation (Fig. 2c). Under laser excitation energies at 0.31âeV and 0.35âeV, we observe similar profiles for \({I}_{{{{{{\rm{x}}}}}}}^{(5\omega )}\) and \({I}_{{{{{{\rm{y}}}}}}}^{(5\omega )}\), but the entire profiles are precisely shifted with the laser excitation energies (see Supplementary Fig. 7 for the excitation energy dependent HHG given by elliptically polarized light). This means that the pathway for the efficient generation of \({I}_{{{{{{\rm{y}}}}}}}^{(5\omega )}\) also requires the multi-photon interband excitation of charge carriers. However, the two contrasting dependences on chemical potential imply that the interplays of microscopic electronic processes are markedly different in the pathways generating \({I}_{{{{{{\rm{x}}}}}}}^{(5\omega )}\) and \({I}_{{{{{{\rm{y}}}}}}}^{(5\omega )}\).
The pathways of charge carrier in real space can be probed by measuring polarization profile of harmonic signals. Ultrafast current (\({{{{{{\bf{J}}}}}}}^{(5{{{{{\rm{\omega }}}}}})}\)) induced by laser field generates harmonic radiation field (\({{{{{{\bf{E}}}}}}}^{(5{{{{{\rm{\omega }}}}}})}\)) satisfying the following equation:
\({J}_{\hat{\varphi }}^{(5\omega )}\) and \({I}_{\hat{\varphi }}^{(5\omega )}\) are ultrafast current and harmonic radiation intensity along direction \(\hat{{{{{{\boldsymbol{\varphi }}}}}}}\) at frequency of \(5\omega\), respectively. In our experiment, \({I}_{\hat{\varphi }}^{(5\omega )}\) can be measured by recording harmonic intensity after a linear polarizer as a function of polarizer angle (\(\varphi\)). Black empty circles in polar plots of Fig. 3a show polarization of fifth harmonics for the charge-neutral case (column labeled with \(2\left|\mu \right|\)â=â0âeV) and the highly doped case (column labeled with \(2\left|\mu \right|\)â=â1.4âeV) as the two representative cases with control of \({\varepsilon }_{{{{{{\rm{exc}}}}}}}\) (see Supplementary Fig. 8 for more measurement at different chemical potentials). The first column (labeled with Excitation) is laser intensity profile which is shown as a guidance to describe the direction of the major axis (the x-direction, \(\varphi\)â=â0°), minor axis (the y-direction, \(\varphi\)â=â90°), and ellipticity of laser excitation. Laser excitation helicity is set as counterclockwise (left-handed) rotation with respect to laser propagation direction. Under linearly-polarized excitation (\({\varepsilon }_{{{{{{\rm{exc}}}}}}}\)â=â0), the harmonic signals for both cases (\(2\left|\mu \right|\)â=â0âeV and \(2\left|\mu \right|\)â=â1.4âeV) exhibit linear polarization along the x-direction, meaning that laser excitation induces \({{{{{{\bf{J}}}}}}}^{(5{{{{{\rm{\omega }}}}}})}\) in graphene along the parallel direction of the laser polarization.
a Polarization profiles of the fifth harmonics. Black circles show the fifth harmonic intensity measured after a linear polarizer as a function of polarizer angle \(\varphi\). Polarization profile of laser excitation with ellipticity \({\varepsilon }_{{{{{{\rm{exc}}}}}}}\) = 0, 0.12, 0.3, and 0.5 (left column labeled with Excitation) is plotted as a guidance. For the charge-neutral case (middle column labeled with \(2\left|\mu \right|\) = 0âeV), the harmonic polarization sensitively rotates with \({\varepsilon }_{{{{{{\rm{exc}}}}}}}\). For the highly doped case (right column labeled with \(2\left|\mu \right|\) = 1.4âeV), the harmonic polarization nearly does not rotate up to \({\varepsilon }_{{{{{{\rm{exc}}}}}}}\)â=â0.5. Pink solid lines show the theoretically calculated polarization, which exhibits good agreement with experimental results. b, c Direction of harmonic polarization rotation under elliptically polarized excitation. Under laser excitation with left-handed helicity in (b) (or right-handed helicity in (c)), the polarization axis rotates in counterclockwise direction (or clockwise direction) with \({\varepsilon }_{{{{{{\rm{exc}}}}}}}\). d, e Harmonic intensity as a function of \({\varepsilon }_{{{{{{\rm{exc}}}}}}}\). Both \({I}_{{{{{{\rm{x}}}}}}}^{(5\omega )}\) (red circles) and \({I}_{{{{{{\rm{y}}}}}}}^{(5\omega )}\) (blue circles) are observed for the charge-neutral case whereas \({I}_{{{{{{\rm{y}}}}}}}^{(5\omega )}\) is nearly absent for the highly doped case. Theoretical calculation results (solid lines) show a good agreement with experimental results.
As laser ellipticity \({\varepsilon }_{{{{{{\rm{exc}}}}}}}\) increases, harmonic signals show completely different behaviors. For the charge-neutral case, harmonic polarization axis sensitively rotates with ellipticity. Under \({\varepsilon }_{{{{{{\rm{exc}}}}}}}\) = 0.5, polarization axis is rotated by 58.7° with respect to the laser polarization axis. The rotation direction of polarization axis is determined by laser excitation helicity. Under laser excitation with the left-handed helicity, polarization axis rotates in the counterclockwise direction, which is also shown in Fig. 3b with additional data (\({\varepsilon }_{{{{{{\rm{exc}}}}}}}\)= 0.2 and 0.4). Under right-handed helicity excitation, the polarization axis rotates by exactly the same angles but in clockwise direction (Fig. 3c). This indicates that elliptical laser excitation can rotate the orientation of \({{{{{{\bf{J}}}}}}}^{(5{{{{{\rm{\omega }}}}}})}\) following the laser field helicity. However, polarization axis nearly does not rotate up to \({\varepsilon }_{{{{{{\rm{exc}}}}}}}\) = 0.5 for the highly doped case (Fig. 3a). The two contrasting behaviors upon \({\varepsilon }_{{{{{{\rm{exc}}}}}}}\) are summarized in Fig. 3d, e where red circles, blue circles, and black squares represent \({I}_{{{{{{\rm{x}}}}}}}^{(5\omega )}\), \({I}_{{{{{{\rm{y}}}}}}}^{(5\omega )}\), and \({I}_{{{{{{\rm{x}}}}}}}^{(5\omega )}\)â+â\({I}_{{{{{{\rm{y}}}}}}}^{(5\omega )}\), respectively. For \(2\left|\mu \right|\)â=â0âeV (Fig. 3d), elliptically polarized excitation efficiently generates \({I}_{{{{{{\rm{y}}}}}}}^{(5\omega )}\) in addition to \({I}_{{{{{{\rm{x}}}}}}}^{(5\omega )}\). For example, substantially smaller laser field along the y-direction efficiently generates \({I}_{{{{{{\rm{y}}}}}}}^{(5\omega )}\) at \({\varepsilon }_{{{{{{\rm{exc}}}}}}}\) = 0.12, which leads to harmonic signals with remarkably large polarization rotation and ellipticity. On the other hand, \({I}_{{{{{{\rm{y}}}}}}}^{(5\omega )}\) is nearly absent for \(2\left|\mu \right|\)â=â1.4âeV (Fig. 3e) for all ellipticities, which orients the harmonic signals nearly along the x-direction. Considering the fact that multi-photon interband transition are suppressed upto the five photon transition (1.38âeV) for \(2\left|\mu \right|\)â=â1.4âeV due to Pauli blocking, charge carrier excitation via the interband transition also plays an important role for charge carriers to generate \({J}_{{{{{{\rm{y}}}}}}}^{(5\omega )}\) efficiently.
In the previous study on graphene23, the efficient generation of \({I}_{{{{{{\rm{y}}}}}}}^{(5\omega )}\) is attributed to transient band structure modification under intense laser excitation (~1âTWâcmâ2). According to the suggested model for a semiconductor52,53, strong laser field closes a bandgap by lowering a conduction band and raising a valence band, which increases overlap between the two bands. Pairs of electrons and holes in this semimetallized structure can be readily generated without the interband transition. This pair generation mechanism can be particularly facilitated in graphene because the conduction band and valence band overlaps at Dirac point without a bandgap. Their following analysis on charge carrier trajectory shows that elliptically polarized excitation with finite ellipticity can exhibit larger recombination probability of charge carriers than linearly-polarized excitation. However, in our experiment with orders of magnitude lower laser intensity (~1âGWâcmâ2), systematic chemical potential dependence of harmonic signals (Fig. 2câe) suggests that charge carriers are dominantly excited via the multi-photon transitions instead of the pair generation in the transiently semimetallized band structure. Therefore, we conclude that graphene hosts an additional efficient pathway to generate high harmonics under elliptically polarized excitation in the multi-photon excitation regime.
Theoretical analysis of HHG in graphene
In order to understand the underlying microscopic mechanism, we solve Liouville-von Neumann equation (see Supplementary Note on Liouville-von Neumann equation) to obtain photo-excited charge carrier population for the charge neutral case where graphene is excited by linearly-polarized laser field along x-direction (\({E}_{{{{{{\rm{x}}}}}}}\)). Population relaxation rate is assumed negligibly small to visualize photo-excited carrier population under strong laser fields. Figure 4a shows the conduction band population (\({\rho }_{{{{{{\rm{CB}}}}}}}\)) immediately after laser excitation with intensity of 3.1âGWâcmâ2. In the momentum space, population is mostly concentrated around the Dirac points at the K and Kâ² points where broad low-energy interband transition is available from terahertz to ultraviolet frequencies along Dirac cones. The detailed profile of carrier population around the K point (region marked with a white square) is shown in Fig. 4b. Dashed lines are constant energy contours on a Dirac cone, which describes the conduction band states vertically separated by energy \({{mE}}_{{{{{{\rm{ph}}}}}}}\) from the valence band. Carriers are mostly distributed from the Dirac point to the states on the contour of \({5E}_{{{{{{\rm{ph}}}}}}}\) along \({k}_{{{{{{\rm{y}}}}}}}\)-axis while carriers are absent along the \({k}_{{{{{{\rm{x}}}}}}}\)-axis. For single-photon transition, the pronounced nodal line can be explained by optical selection rule of massless Dirac fermion. Direction of transition dipole moment circulates around the Dirac point, which prohibits interband transition for states with \({k}_{{{{{{\rm{y}}}}}}}\) = 0 under linearly-polarized laser field along the x-direction (\({E}_{{{{{{\rm{x}}}}}}}\)) (Fig. 4a). Such characteristic excitation profile from the optical selection rule is also observed in ARPES measurements54 and ultrafast pump-probe studies55. Our numerical calculation in Fig. 4a and b indicates that the multi-photon transition as well as the single-photon transition is also blocked for the states on the \({k}_{{{{{{\rm{x}}}}}}}\)-axis.
a Photo-excited charge carrier population in the conduction band \({\rho }_{{{{{{\rm{CB}}}}}}}\) under linearly-polarized laser field along the x-direction (\({E}_{{{{{{\rm{x}}}}}}}\)). b Detailed profile of carrier population around K point (region marked with a white square in a). Dashed lines are constant energy contours on a Dirac cone, which describes conduction band states vertically separated by energy \({{mE}}_{{{{{{\rm{ph}}}}}}}\). c Transition dipole moment profile. Color scale and arrow show the magnitude and direction of transition dipole moment. d, e Microscopic channels for carrier dynamics at the states on the \({k}_{{{{{{\rm{y}}}}}}}\)-axis. Under linearly-polarized excitation along the x-direction (d), the interband transition (black arrows labeled with \({M}_{{{{{{\rm{x}}}}}}}^{{{{{{\rm{inter}}}}}}}\)) creates and recombine photo-excited electrons and holes (red filled and empty circles, respectively). Simultaneously, the intraband transition drives carriers along the \({k}_{{{{{{\rm{x}}}}}}}\)-direction (orange arrows labeled with \({M}_{{{{{{\rm{x}}}}}}}^{{{{{{\rm{intra}}}}}}}\)). Under elliptically polarized excitation (e), the intraband transition also drives carriers along the \({k}_{{{{{{\rm{y}}}}}}}\)-direction (blue arrows labeled with \({M}_{{{{{{\rm{y}}}}}}}^{{{{{{\rm{intra}}}}}}}\)). f Polarization profile of the fifth harmonics under laser ellipticity \({\varepsilon }_{{{{{{\rm{exc}}}}}}}\)â=â0.2 from selected pathways. Black, orange, and blue solid lines show the polarization profile from quantum pathways via \({M}_{{{{{{\rm{x}}}}}}}^{{{{{{\rm{inter}}}}}}}\)â+â\({M}_{{{{{{\rm{y}}}}}}}^{{{{{{\rm{inter}}}}}}}\), \({M}_{{{{{{\rm{x}}}}}}}^{{{{{{\rm{inter}}}}}}}\)â+â\({M}_{{{{{{\rm{x}}}}}}}^{{{{{{\rm{intra}}}}}}}\), and \({M}_{{{{{{\rm{x}}}}}}}^{{{{{{\rm{inter}}}}}}}\)â+â\({M}_{{{{{{\rm{y}}}}}}}^{{{{{{\rm{intra}}}}}}}\), respectively. Black empty circles show experimental results for the same laser excitation condition.
The population relaxation processes quickly relax the anisotropic population (Fig. 4a, b) to a thermalized isotropic population within ~30âfs54,55,56, which typically generates the same photocurrent regardless of direction of an external electric field57. However, the anisotropic excitation profile plays an important role in high harmonic generation because charge carriers under strong laser fields radiate harmonic fields by interband and intraband transitions in sub-laser-cycle timescale (<15âfs in our experiment)45. Figure 4d, e schematically describes the microscopic channels for carrier dynamics at the states on the \({k}_{{{{{{\rm{y}}}}}}}\)-axis (i.e. \({k}_{{{{{{\rm{x}}}}}}}\) \(=\) 0, \({k}_{{{{{{\rm{y}}}}}}}\) \(\ne\) 0) as the representative cases. Under linearly-polarized excitation along x-direction (Fig. 4d), the interband transition (black arrows labeled with \({M}_{{{{{{\rm{x}}}}}}}^{{{{{{\rm{inter}}}}}}}\)) creates and recombines photo-excited electrons and holes (red filled and empty circles, respectively). Simultaneously, laser field \({E}_{{{{{{\rm{x}}}}}}}\) drives photo-excited carriers via the intraband transition along the \({k}_{{{{{{\rm{x}}}}}}}\)-direction (orange arrows labeled with \({M}_{{{{{{\rm{x}}}}}}}^{{{{{{\rm{intra}}}}}}}\)). Under elliptically polarized excitation (Fig. 4e), laser field \({E}_{{{{{{\rm{y}}}}}}}\) also drives carriers via the intraband transition along the \({k}_{{{{{{\rm{y}}}}}}}\)-direction with a phase delay of Ï/2 (blue arrows labeled with \({M}_{{{{{{\rm{y}}}}}}}^{{{{{{\rm{intra}}}}}}}\)). The interband transition from laser field \({E}_{{{{{{\rm{y}}}}}}}\) (\({M}_{{{{{{\rm{y}}}}}}}^{{{{{{\rm{inter}}}}}}}\)) is not allowed for the states as shown in Fig. 4b. The intricate coupling of \({M}_{{{{{{\rm{x}}}}}}}^{{{{{{\rm{inter}}}}}}}\) and \({M}_{{{{{{\rm{y}}}}}}}^{{{{{{\rm{intra}}}}}}}\) provides an additional pathway to generate ultrafast current that is not available under linearly-polarized excitation. In particular, the population profile from \({M}_{{{{{{\rm{x}}}}}}}^{{{{{{\rm{inter}}}}}}}\) (Fig. 4b) can exhibit highly anisotropic response in terms of ultrafast current generation along the two orthogonal directions via \({M}_{{{{{{\rm{x}}}}}}}^{{{{{{\rm{intra}}}}}}}\) and \({M}_{{{{{{\rm{y}}}}}}}^{{{{{{\rm{intra}}}}}}}\).
The harmonic signals from graphene can be theoretically calculated by solving density matrices with quantum master equation in the Houston basis (see Supplementary Note for more details on theoretical calculation). Under the same laser excitation condition in our experiment, the theoretical calculation (pink solid lines in Fig. 3a, solid lines in Fig. 3d, e) reproduces our experimental results with an excellent agreement. Detailed chemical potential dependence of \({I}_{{{{{{\rm{x}}}}}}}^{(5\omega )}\) and \({I}_{{{{{{\rm{y}}}}}}}^{(5\omega )}\) also shows the characteristic features from charge carrier excitation via the multi-photon transition, which is consistent with our experimental observation in Fig. 2câe (see Supplementary Note for calculation of chemical potential dependent HHG).
In order to resolve the individual contributions from different pathways under elliptically polarized excitation, we employ the uniform velocity gauge where Bloch functions and gauges are fully decomposed. In the velocity gauge, light-electron interaction can be written in Hamiltonian, \(p\cdot {{{{{\bf{A}}}}}}\left(\tau \right)/c\). \({{{{{\bf{A}}}}}}(\tau )\) and \(c\) are the vector potential of laser field and the speed of light, respectively, and \(p\) is the optical matrix element of Dirac electrons in graphene. This representation allows to distinguish charge carrier dynamics naturally via the interband and intraband transitions induced by the time-dependent vector potential of a single external laser pulse. \(p\cdot {{{{{\bf{A}}}}}}\left(\tau \right)/c\) can be decomposed into the arbitrarily model gauge based on the interband and intraband transitions under laser fields along the x- and y-directions (i.e., \({p}_{{{{{{\rm{x}}}}}}}^{{{{{{\rm{inter}}}}}}}{A}_{{{{{{\rm{inter}}}}}}}\left(\tau \right)/c,\) \({p}_{y}^{{{{{{\rm{inter}}}}}}}{A}_{{{{{{\rm{inter}}}}}}}\left(\tau \right)/c\), \({p}_{{{{{{\rm{x}}}}}}}^{{{{{{\rm{intra}}}}}}}{A}_{{{{{{\rm{intra}}}}}}}\left(\tau \right)/c,\) and \({p}_{{{{{{\rm{y}}}}}}}^{{{{{{\rm{intra}}}}}}}{A}_{{{{{{\rm{intra}}}}}}}\left(\tau \right)/c\)) so that the interaction part of the Hamiltonian can be rewritten by the microscopic channels of \({M}_{{{{{{\rm{x}}}}}}}^{{{{{{\rm{inter}}}}}}},\) \({M}_{{{{{{\rm{y}}}}}}}^{{{{{{\rm{inter}}}}}}}\), \({M}_{{{{{{\rm{x}}}}}}}^{{{{{{\rm{intra}}}}}}}\), and \({M}_{{{{{{\rm{y}}}}}}}^{{{{{{\rm{intra}}}}}}}\). One can examine harmonic signals generated from individual decomposed pathways selectively by constructing a desired Hamiltonian from linear combination of microscopic channels.
Solid lines in Fig. 4f show polar plots of fifth harmonic signals as a function of linear polarizer angle (\(\varphi\)), which are calculated for representative pathways under laser excitation with \({\varepsilon }_{{{{{{\rm{exc}}}}}}}\) = 0.2. Black circles in Fig. 4f show the fifth harmonic signals measured in experiment as a reference. The pathway of \({M}_{{{{{{\rm{x}}}}}}}^{{{{{{\rm{inter}}}}}}}\)â+â\({M}_{{{{{{\rm{x}}}}}}}^{{{{{{\rm{intra}}}}}}}\) (orange solid line) generates a linearly-polarized harmonic signal along the x-direction despite the presence of both laser fields \({E}_{{{{{{\rm{x}}}}}}}\) and \({E}_{{{{{{\rm{y}}}}}}}\), indicating that \({I}_{{{{{{\rm{y}}}}}}}^{(5\omega )}\) cannot be generated without \({M}_{{{{{{\rm{y}}}}}}}^{{{{{{\rm{inter}}}}}}}\) or \({M}_{{{{{{\rm{y}}}}}}}^{{{{{{\rm{intra}}}}}}}\). The pathway of \({M}_{{{{{{\rm{x}}}}}}}^{{{{{{\rm{inter}}}}}}}\)â+â\({M}_{{{{{{\rm{y}}}}}}}^{{{{{{\rm{inter}}}}}}}\) (black solid line) generates near linearly-polarized harmonic along a direction slightly rotated from x-direction. Significantly weaker \({E}_{{{{{{\rm{y}}}}}}}\), which is five times smaller than \({E}_{{{{{{\rm{x}}}}}}}\), generates negligibly small but non-zero \({I}_{{{{{{\rm{y}}}}}}}^{(5\omega )}\). However, the pathway of \({M}_{{{{{{\rm{x}}}}}}}^{{{{{{\rm{inter}}}}}}}\)â+â\({M}_{{{{{{\rm{y}}}}}}}^{{{{{{\rm{intra}}}}}}}\) (blue solid line) generates surprisingly strong \({I}_{{{{{{\rm{y}}}}}}}^{(5\omega )}\) with polarization major axis rotated by 35°. This clearly shows that charge carriers in graphene exhibit the intricate coupling between \({M}_{{{{{{\rm{x}}}}}}}^{{{{{{\rm{inter}}}}}}}\) and \({M}_{{{{{{\rm{y}}}}}}}^{{{{{{\rm{intra}}}}}}}\) under elliptically polarized excitation. Although calculation of accurate harmonic signals requires full Hamiltonian with all four microscopic channels (\({M}_{{{{{{\rm{x}}}}}}}^{{{{{{\rm{inter}}}}}}}\), \({M}_{{{{{{\rm{y}}}}}}}^{{{{{{\rm{inter}}}}}}}\), \({M}_{{{{{{\rm{x}}}}}}}^{{{{{{\rm{intra}}}}}}}\), and \({M}_{{{{{{\rm{y}}}}}}}^{{{{{{\rm{intra}}}}}}}\)), it is remarkable to note that the pathway of \({M}_{{{{{{\rm{x}}}}}}}^{{{{{{\rm{inter}}}}}}}\)â+â\({M}_{{{{{{\rm{y}}}}}}}^{{{{{{\rm{intra}}}}}}}\) alone reproduces our experiment result (black circles) with a good agreement. Our theoretical calculation result agrees with recent theoretical results22 where the strong coupling is attributed to large intraband transition strength along the \({k}_{{{{{{\rm{y}}}}}}}\)-direction (blue arrow in Fig. 4e (\({M}_{{{{{{\rm{y}}}}}}}^{{{{{{\rm{intra}}}}}}}\))) than along the \({k}_{{{{{{\rm{x}}}}}}}\)-direction (orange arrow in Fig. 4e (\({M}_{{{{{{\rm{y}}}}}}}^{{{{{{\rm{intra}}}}}}}\))) for the states excited by interband transition (black arrow in Fig. 4e (\({M}_{{{{{{\rm{x}}}}}}}^{{{{{{\rm{inter}}}}}}}\))).
In conclusion, the chemical potential control in our study revealed the quantum pathways of massless Dirac fermions in graphene for HHG under strong laser fields. Charge carriers are excited via the multi-photon interband transitions, which forms destructively interfering quantum pathways for HHG. As the light-matter interaction deviates beyond the perturbative regime, the highly anisotropic carrier excitation profile of massless Dirac fermions in momentum space leads to the intricate coupling between the interband and intraband transitions, which allows to control HHG in terms of intensity and polarization state under elliptically polarized excitation. Considering the fact that the circulating profile of the transition dipole moment (Fig. 4c) universally presents in solids hosting gapless charge carriers with topologically-protected band crossings (see Supplementary note for transition dipole moment of massless Dirac fermions), the underlying principles of HHG in graphene can be also utilized for broad range of quantum semimetals including topological insulators and Weyl semimetals to resolve quantum pathways of Dirac electrons under strong laser fields. In addition, the intensity and polarization modulation scheme of high harmonics provides a microscopic mechanism to enable a novel light source based on HHG in quantum semimetals with a fast and efficient control via an electrostatic gating.
Methods
Device fabrication
Polycrystalline graphene film is grown by chemical vapor deposition on Cu foil (Nilaco corporation, #CU-113213, 30âµm thick, 99.9% purity). The Cu foil is annealed at 1030â°C for 4-hour under flow of H2 at 70 sccm with a total pressure of 5âTorr; then graphene film is grown at 1040â°C for 2-h 30âmin by additionally introducing methane in Ar carrier gas at 0.15 sccm. During the cooling process after growth, the methane injection is maintained until the film had cooled to 600â°C, then only H2 is introduced until the temperature reached <150â°C to allow unloading of the sample. Poly (methyl methacrylate) (PMMA, 996âK, 8% in Anisole) is spun onto as grown graphene. To etch the Cu substrate, the PMMA/graphene/Cu foil is floated on a Cu etchant of FeCl3 aqueous solution (Sigma-Aldrich, #667528) for 30âmin. After complete etching of Cu, the PMMA/graphene film is transferred onto a surface of ultrahigh purity deionized water (DIW) by scooping with a SiOx substrate, then and releasing the ultrathin film. The process is repeated twice in a clean DIW to rinse the surface, then the film is transferred onto a double-side polished sapphire substrate. The sample is dried by heating at 60â°C in air for 10âmin, followed by annealing at 165â°C for 10âmin and soaking in acetone at 60â°C for 10âmin to remove the PMMA. Finally, annealing is performed at 350â°C in air for 15âmin to remove the PMMA residue.
After transfer process, drain, source, and side gate contacts are deposited by e-beam evaporator (Cr 5ânm/Au 50ânm) to construct ion-gel-based field-effect transistor and distance between drain and source is about a few millimeters. For ion-gel gate dielectric, 0.1âg of poly(vinylidene fluoride-co-hexafluoropropylene) (PVDF-HFP) and 0.195âg of 1-ethyl-3-methylimidazolium bis (trifluoromethylsulfonyl)imide ([EMIM][TFSI]) are dissolved in the mixture of 2âml of 2-butanon and 0.25âml of propylene carbonate. [EMIM][TFSI] and PVDF-HFP were purchased from Sigma Aldrich and ARKEMA, respectively. The solution is stirred overnight at 45â°C, then drop-casted onto the device to cover gate electrodes and graphene channel. The sample is dried in argon filled glove box. Source, drain, and gate contacts are electrically connected by pair of voltage source meter (Keithley 2400, 2450) to apply gate voltage and read channel resistance.
HHG measurements
Based on femtosecond laser system (Light Conversion PHAROS), mid-infrared pulses are prepared using optical parametric amplifier (ORPHEUS) and difference frequency generator (LYRA). The output serves wavelength-tunable multi-cycle pulses with repetition rate of 100âkHz. The spectral linewidth of the pulse is 15.4âmeV in full-width half-maximum and the pulse duration is estimated to be 120âfs assuming a Fourier-transform-limited pulse. To control its ellipticity, liquid crystal retarder (Thorlabs LCC1111-MIR) is employed, whose optical axis is oriented at an angle of 45° with respect to the laser polarization. Then, mid-infrared pulses are focused at roughly center of the graphene device by ZnSe focusing objectives with spot size of 150 \(\mu\)m. Emitted HHG has been collected by 50X objective lens on transmission geometry, and its polarization is analyzed by half-wave plate mounted on motorized stage and fixed Glan-Taylor polarizer. The HHG spectra are recorded by an electron-multiplying charge-coupled device detector (ProEM, Princeton instruments) and grating spectrometer (SP-2300, Princeton instruments) at Materials Imaging & Analysis Center of POSTECH.
Data availability
All of the data that support the findings of this study are available in the main text or Supplementary Information. Source data are available from the corresponding authors on request.
Code availability
Custom codes used in this work can be provided by the corresponding author on request.
References
Ferray, M. et al. Multiple-harmonic conversion of 1064 nm radiation in rare gases. J. Phys. B 21, L31 (1988).
Corkum, P. B. & Krausz, F. Attosecond science. Nat. Phys. 3, 381â387 (2007).
Keldysh, L. Ionization in the field of a strong electromagnetic wave. Sov. Phys. JETP 20, 1307â1314 (1965).
Corkum, P. B. Plasma perspective on strong field multiphoton ionization. Phys. Rev. Lett. 71, 1994â1997 (1993).
Krausz, F. & Stockman, M. I. Attosecond metrology: from electron capture to future signal processing. Nat. Photon. 8, 205â213 (2014).
Ghimire, S. et al. Observation of high-order harmonic generation in a bulk crystal. Nat. Phys. 7, 138â141 (2011).
Ghimire, S. & Reis, D. A. High-harmonic generation from solids. Nat. Phys. 15, 10â16 (2019).
Luu, T. T. et al. Extreme ultraviolet high-harmonic spectroscopy of solids. Nature 521, 498â502 (2015).
Vampa, G. et al. All-optical reconstruction of crystal band structure. Phys. Rev. Lett. 115, 193603 (2015).
Borsch, M. et al. Super-resolution lightwave tomography of electronic bands in quantum materials. Science 370, 1204â1207 (2020).
Luu, T. T. & Wörner, H. J. Measurement of the Berry curvature of solids using high-harmonic spectroscopy. Nat. Commun. 9, 916 (2018).
Liu, H. et al. High-harmonic generation from an atomically thin semiconductor. Nat. Phys. 13, 262â265 (2017).
Schmid, C. P. et al. Tunable non-integer high-harmonic generation in a topological insulator. Nature 593, 385â390 (2021).
Silva, R. E. F., Jiménez-Galá, Ã., Amorim, B., Smirnova, O. & Ivanov, M. Topological strong-field physics on sub-laser-cycle timescale. Nat. Photon. 13, 849â854 (2019).
Chacón, A. et al. Circular dichroism in higher-order harmonic generation: Heralding topological phases and transitions in Chern insulators. Phys. Rev. B 102, 134115 (2020).
Zaks, B., Lium, R. B. & Sherwin, M. S. Experimental observation of electronâhole recollisions. Nature 483, 580â583 (2012).
Costello, J. B. et al. Reconstruction of Bloch wavefunctions of holes in a semiconductor. Nature 599, 57â61 (2021).
Langer, F. et al. Lightwave-driven quasiparticle collisions on a subcycle timescale. Nature 533, 225â229 (2016).
Tancogne-Dejean, N., Mücke, O. D., Kärtner, F. X. & Rubio, A. Ellipticity dependence of high-harmonic generation in solids originating from coupled intraband and interband dynamics. Nat. Commun. 8, 745 (2017).
You, Y. S., Reis, David, A. & Ghimire, S. Anisotropic high-harmonic generation in bulk crystals. Nat. Phys. 13, 345â349 (2017).
Vampa, G. et al. Theoretical analysis of high-harmonic generation in solids. Phys. Rev. Lett. 113, 073901 (2014).
Sato, S. A., Hirori, H., Sanari, Y., Kanemitsu, Y. & Rubio, A. High-order harmonic generation in graphene: Nonlinear coupling of intraband and interband transitions. Phys. Rev. B 103, L041408 (2021).
Yoshikawa, N., Tamaya, T. & Tanaka, K. High-harmonic generation in graphene enhanced by elliptically polarized light excitation. Science 356, 736â738 (2017).
Bai, Y. et al. High-harmonic generation from topological surface states. Nat. Phys. 17, 311â315 (2021).
Uchida, K., Pareek, V., Nagai, K., Dani, K. M. & Tanaka, K. Visualization of two-dimensional transition dipole moment texture in momentum space using high-harmonic generation spectroscopy. Phys. Rev. B 103, L161406 (2021).
Novoselov, K. S. et al. Two-dimensional gas of massless Dirac fermions in graphene. Nature 438, 197â200 (2005).
Mak, K. F. et al. Measurement of the optical conductivity of graphene. Phys. Rev. Lett. 101, 196405 (2008).
Horng, J. et al. Drude conductivity of Dirac fermions in graphene. Phys. Rev. B 83, 165113 (2011).
Jiang, T. et al. Gate-tunable third-order nonlinear optical response of massless Dirac fermions in graphene. Nat. Photon. 12, 430â436 (2018).
Soavi, G. et al. Broadband, electrically tunable third-harmonic generation in graphene. Nat. Nanotechnol. 13, 683â588 (2018).
Hong, S.-Y. et al. Optical third-harmonic generation in graphene. Phys. Rev. X 3, 021014 (2013).
Hendry, E., Hale, P. J., Moger, J., Savchenko, A. K. & Mikhailov, S. A. Coherent nonlinear optical response of graphene. Phys. Rev. Lett. 105, 097401 (2010).
Hafez, H. A. et al. Extremely efficient terahertz high-harmonic generation in graphene by hot Dirac fermions. Nature 561, 507â511 (2018).
Kovalev, S. et al. Electrical tunability of terahertz nonlinearity in graphene. Sci. Adv. 7, eabf9809 (2021).
Feng, Y. et al. Semiclassical analysis of ellipticity dependence of harmonic yield in graphene. Phys. Rev. A 104, 043525 (2021).
Chizhova, L. A., Libisch, F. & Burgdörfer, J. High-harmonic generation in graphene: Interband response and the harmonic cutoff. Phys. Rev. B 95, 085436 (2017).
Mrudul, M. S. & Dixit, G. High-harmonic generation from monolayer and bilayer graphene. Phys. Rev. B 103, 094308 (2021).
Liu, C., Zheng, Y., Zeng, Z. & Li, R. Driving-laser ellipticity dependence of high-order harmonic generation in graphene. Phys. Rev. A 97, 063412 (2018).
Al-Naib, I., Sipe, J. E. & Dignam, M. M. Nonperturbative model of harmonic generation in undoped graphene in the terahertz regime. N. J. Phys. 17, 113018 (2015).
Hipolito, F., Dimitrovski, D. & Pedersen, T. G. Iterative approach to arbitrary nonlinear optical response functions of graphene. Phys. Rev. B 99, 195407 (2019).
Zhang, Y. et al. Orientation dependence of high-order harmonic generation in graphene. Phys. Rev. A 104, 033110 (2021).
Wang, X. âQ. & Bian, X. âB. Model-potential method for high-order harmonic generation in monolayer graphene. Phys. Rev. A 103, 053106 (2021).
Wang, F. et al. Gate-variable optical transitions in graphene. Science 320, 206â209 (2008).
Li, Z. Q. et al. Dirac charge dynamics in graphene by infrared spectroscopy. Nat. Phys. 4, 532â535 (2008).
Chang, Z. Fundamentals of attosecond optics (CRC, 2011).
Cheng, B. et al. Efficient terahertz harmonic generation with coherent acceleration of electrons in the Dirac semimetal Cd3As2. Phys. Rev. Lett. 124, 117402 (2020).
Shen, Y. R. The Principles of Nonlinear Optics (John Wiley and Sons, 1984).
Cheng, J. L., Vermeulen, N. & Sipe, J. E. Third order optical nonlinearity of graphene. N. J. Phys. 16, 053014 (2014).
Burnett, N. H., Kan, C. & Corkum, P. B. Ellipticity and polarization effects in harmonic generation in ionizing neon. Phys. Rev. A 51, R3418âR3421 (1995).
Baykusheva, D. et al. All-optical probe of three-dimensional topological insulators based on high-harmonic generation by circularly polarized laser fields. Nano Lett. 21, 8978â8970 (2021).
Li, L. et al. Reciprocal-space-trajectory perspective on high-harmonic generation in solids. Phys. Rev. Lett. 122, 193901 (2019).
Tamaya, T., Ishikawa, A., Ogawa, T. & Tanaka, K. Higher-order harmonic generation caused by elliptically polarized electric fields in solid-state materials. Phys. Rev. B 94, 241107 (2016).
Tamaya, T., Ishikawa, A., Ogawa, T. & Tanaka, K. Diabatic mechanisms of higher-order harmonic generation in solid-state materials under high-intensity electric fields. Phys. Rev. Lett. 116, 016601 (2016).
Aeschlimann, S. et al. Ultrafast momentum imaging of pseudospin-flip excitations in graphene. Phys. Rev. B 96, 020301 (2017). (R).
Mittendorff, M. et al. Anisotropy of excitation and relaxation of photogenerated charge carriers in graphene. Nano Lett. 14, 1504â1507 (2014).
Heide, C. et al. Electronic coherence and coherent dephasing in the optical control of electrons in graphene. Nano Lett. 21, 9403â9409 (2021).
Koppens, F. H. L. et al. Photodetectors based on graphene, other two-dimensional materials and hybrid systems. Nat. Nanotechnol. 9, 780â793 (2014).
Acknowledgements
J.K. and M.K. acknowledge help from D.P. and Prof. U.J. for graphene device preparation. S.C., M.K., M.-H.J., and J.K. acknowledge the support from Institute of Basic Science (IBS-R034-D1). J.K. and M.K. acknowledge the support from the National Research Foundation of Korea grants (NRF-2020R1A4A1018935, and 2020R1A2C2103166). This study was supported by Brain Korea 21 FOUR project for Education and research center for future materials (F21YY7105002). This research was also supported by the MSIT (Ministry of Science and ICT), Korea, under the ITRC (Information Technology Research Center) support program (IITP-2022-RS-2022-00164799) supervised by the IITP (Institute for Information & Communications Technology Planning & Evaluation). Y.K. and J.L. acknowledge the support from the DGIST R&D program (22-CoE-NT-01). J.K. and J.L. also acknowledge the support from the national research foundation (NRF) of Korea (2022R1A2B5B01001582). G.Y.C. acknowledges the support of the National Research Foundation of Korea (NRF) funded by the Korean Government (Grant No. 2020R1C1C1006048) and Grant No. IBS-R014-D1. G.Y.C. acknowledges the Air Force Office of Scientific Research under Award No. FA2386-20-1-4029 and No. FA2386-22-1-4061. G.Y.C. also thanks the support by Samsung Science and Technology Foundation under Project Number SSTF-BA2002-05. C.-J.K. acknowledges the support from the National Research Foundation of Korea grants (2020R1C1C1014590 and 2022M3H4A1A01012718).
Author information
Authors and Affiliations
Contributions
J.K. conceived the project. S.C. and M.K. built optical set-up and obtained HHG spectra. S.C., M.K., S.C., and S.K. fabricated graphene devices with a gate control. Y.K., G.C., and H.K. obtained theoretical calculation results of HHG in graphene. C.-J.K., M.J.P., S.C., M.K., S.C., and S.K. analyzed gating behavior of graphene with an ion-gel under strong laser field. J.K., M.-H.J., J.D.L., B.J.K., G.C., S.C., M.K., Y.K., and H.K. analyzed HHG data. All authors discussed and wrote the manuscript together.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Communications thanks the anonymous reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.
Additional information
Publisherâs note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights 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/.
About this article
Cite this article
Cha, S., Kim, M., Kim, Y. et al. Gate-tunable quantum pathways of high harmonic generation in graphene. Nat Commun 13, 6630 (2022). https://doi.org/10.1038/s41467-022-34337-y
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-022-34337-y
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.