Distinguishing high-harmonic generation from surface and bulk states in topological insulator Bi2Se3

Yang Jiang; Ya Bai; Zeyi Ye; Na Li; Candong Liu; Peng Liu

doi:10.3788/COL202321.043801

Abstract

We demonstrated a scheme to differentiate the high-harmonic generation (HHG) originating from the surface states and bulk states of the topological insulator

{Bi}_{2} {Se}_{3}

. By adopting two-color mid-infrared laser fields on

{Bi}_{2} {Se}_{3}

, we found that the nonlinear response sensitively depends on the relative phase of the driving fields. The even harmonics arise from the surface states with a clear signature, whose modulation period equals the cycle of the second-harmonic generation (SHG) field. We reveal that the weak SHG perturbs the nontrivial dipole phase of the electron-hole pair in surface states, and thus leads to the modulation of HHG. It provides a means to manipulate the ultrafast dynamics in surface states through adopting a weak perturbing laser field.

Keywords

high-harmonic generation topological insulators two-color laser fields

1. Introduction

High-harmonic generation (HHG) is one of the most intriguing topics among the phenomena arising from the interaction of matter with intense laser fields. HHG originating from the ultrafast electron dynamics of atoms or molecules can be used to study the field-driven physics between strong laser pulse and matter. It has irreplaceable applications in many research fields, such as coherent extreme ultraviolet (EUV) to X-ray^[1–3], chemical reaction tracing^[4], molecular orbital imaging^[5–7], and attosecond science^[8–12]. In recent years, the nonperturbative HHG has also been observed in solid materials^[13,14]. It has been proved to be a useful tool to reconstruct the electronic band structures and reveal the ultrafast electron dynamics in condensed matter all-optically^[13–19].

HHG from solids is fundamentally different from that in the gaseous medium, owing to a variety of electronic structures in the lattice of bulk crystals^{[13–17,20–29]}. Both the intraband current and interband polarization of the optically excited carriers can contribute to the emission of high-harmonics, despite the two mechanisms often intertwining^{[20–23,30,31]}. In 2015, Vampa et al. found that the generalized recollision between the field-driven electron-hole pairs is the primary mechanism of HHG in ZnO^[24], which brings new ideas for the study of solid-state HHG. In their study, two-color laser fields composed of the fundamental and second-harmonic generation (SHG) were adopted as a pump. The weak SHG can unbalance the electrons’ trajectories in adjacent half-cycles of driving field, and thus lead to the emission of even harmonics. Latterly, the modulation of high-harmonic emission was found to be sensitive to the symmetry properties of solid materials, which provides a method to analyze the mechanism of HHG^[25–27]. Recently, two-color laser fields driving HHG have been adopted to disentangle the electron-hole dynamics from the surface charge region and the crystal bulk^[27]. This has been used to establish the all-optical spectroscopy of a strongly driven crystal and reveal a laser-induced modification of the band structure^[32].

Three-dimensional (3D) topological insulators (TIs) feature unconventional two-dimensional topological surface states (TSSs). In TSSs, the carriers are massless Dirac fermions, which provide a dissipationless transport channel on the TIs’ surface^[33,34]. Very recently, the observations of HHG from 3D TIs have been reported by several groups, and experimental evidence of HHG originating from strong field dynamics of TSSs has been demonstrated^[35–37]. Despite the exciting progress, strong laser fields can inevitably excite transitions between the bulk states and surface states and, in consequence, affect the nonlinear response. Thus, how to disentangle the bulk and surface contributions of HHG is still an open question.

Sign up for Chinese Optics Letters TOC. Get the latest issue of Chinese Optics Letters delivered right to you！Sign up now

2. Experiments and Methods

In this work, we investigated the HHG from ${Bi}_{2} {Se}_{3}$ crystal, a prototypical 3D TI, in response to two-color laser fields. The modulation of even harmonics sensitively depends on the retarded phase, $ϕ$ , between the fundamental and SHG fields. By adjusting the SHG intensity, two periodic modulations were observed. We found that the harmonic yield originating from the surface states modulates once per cycle of the SHG, while the harmonic yield from the bulk states modulates twice per cycle. The observations are well explained by the generalized recollision model based on the semiconductor Bloch equations^[24,38,39].

The experimental setup is shown in Fig. 1. The fundamental pulse was produced at the wavelength of 3820 nm from a home-built system that is composed of two identical three-stage optical parametric amplifiers (OPAs) and a difference frequency generation (DFG) module^[40]. Subsequently, the fundamental pulse passed through an ${AgGaS}_{2}$ (AGS) crystal with the thickness of 400 µm to produce SHG at the wavelength of 1910 nm. Both the fundamental and SHG pulses were synchronized and collinearly focused onto a ${Bi}_{2} {Se}_{3}$ crystal. A zinc sulfide (ZnS) plate was used to compensate for the group delay between the fundamental and SHG pulses. The dual-wavelength wave plate (DWP), a half-wave plate for the fundamental and a full-wave plate for the SHG, was adopted to adjust the polarization of the two-color laser fields. By tuning the angle of the DWP and a polarizer (WP25M-UB, Thorlabs), the SHG intensity was maintained much more weakly than the fundamental intensity ( $< 2 \times 10^{- 2}$ ) as a perturbation factor. Phase delay of the two pulses was tuned by a pair of ${BaF}_{2}$ wedges with an angle of 3°. In the measurements, the phase delay was adjusted by a step of $0.1 π$ . The ${Bi}_{2} {Se}_{3}$ crystal was prepared by using the self-flux method described in Ref. [35]. The emitted high-harmonic spectra were measured by a grating spectrometer equipped with a CCD camera (Andor Shamrock 303i, Newton 920).

Figure 1.Experimental setup. The fundamental and SHG are linearly polarized in s-polarization direction. The mirror plane of Bi₂Se₃ is rotated to be parallel to the polarization of the pump field.

Figure 2 shows a typical HHG spectrum pumped by the two-color laser fields; as a comparison, the HHG spectrum for fundamental pulse only is also shown. In measurements, the SHG intensity is $3.2 \times 10^{- 4}$ times of the fundamental intensity, which is not enough to independently produce HHG. Driven by the fundamental pulse alone, even harmonics should stem from inversion asymmetric states confined in the surface^[35,36]. As the weak SHG is added, notable changes in even harmonics are observed with the increasing yield of the fourth and eighth harmonic (H4, H8), while the sixth harmonic (H6) is opposite. Compared with the even harmonics, the odd harmonics’ intensity remains unchanged. This indicates that even harmonics are sensitive to the weak SHG field, which may originate from the surface states that lack inversion symmetry, while the odd harmonics may come from bulk states with inversion and time-reversal symmetry ^[24,27]. The variation of even harmonics can be interpreted qualitatively by the interference effects of HHG produced from asymmetric two-color laser fields^[41].

$HHG spectrum with and without the SHG, respectively; the intensity of the fundamental is 77 GW/cm2; gray line, the SHG intensity is 0; green line, the SHG intensity is 0.025 GW/cm2; ϕ = π. The small peak at 1.48 eV is the second-order diffraction of H9 from the grating-based spectrometer.$

Figure 2.HHG spectrum with and without the SHG, respectively; the intensity of the fundamental is 77 GW/cm²; gray line, the SHG intensity is 0; green line, the SHG intensity is 0.025 GW/cm²; ϕ = π. The small peak at 1.48 eV is the second-order diffraction of H9 from the grating-based spectrometer.

By changing the phase delay of the two-color fields, periodic modulations of HHG were observed experimentally, as shown in Figs. 3(a)–3(c). For weak SHG intensity of $0.025 {GW / cm}^{2}$ , the yield of H4 and H6 oscillates once per cycle of the SHG, that is, the modulation period equals $T_{2 ω}$ . However, the modulation of H5 is absent in the same measurement. As the SHG intensity is increased to $1.2 {GW / cm}^{2}$ ( $1.6 \times 10^{- 2}$ times that of the fundamental), the modulation of H5 appears. The intensity of H5 modulates twice per cycle of the SHG, that is, the modulation period is $T_{2 ω} / 2$ . The modulation depth of H6 is enhanced; however, the modulation period does not change. At the same time, the intensity modulation of H4 shows a distinctive feature: two extrema emerge in one SHG cycle. By Fourier transform, two kinds of modulation periods can be extracted: one is $T_{2 ω}$ ; the other is $T_{2 ω} / 2$ . Therefore, we speculate that for very weak SHG, even harmonics with $T_{2 ω}$ modulation originate from asymmetric surface states. With the increase of the SHG intensity, the contribution from the bulk states of ${Bi}_{2} {Se}_{3}$ begins to play a role. Therefore, the multiperiodic modulation of even harmonics can be attributed to the interplay between the surface states and bulk states.

Figure 3.Harmonic yield versus the phase delay of two-color laser fields. (a) Modulation of H4 for the SHG intensity of 0.025 GW/cm² (blue square and line) and 1.2 GW/cm² (orange square and line); (b), (c) Modulation of H5 (b) and H6 (c) with the same experimental condition as (a).

In order to understand why there are two kinds of modulations for the even harmonics, we discussed the HHG mechanisms by using semiconductor Bloch equations (SBEs)^[28]. It is known that the modulations of HHG are closely related to the symmetry of solids^[25-27]. According to Ref. [28], the complex transition dipole moment (TDM) $D_{c v} (k)$ that describes the interband coherence contributes to even harmonics from crystals that lack inversion symmetry^[28]. For the 3D TI, the symmetry is different for the electronic states inside the bulk and confined at the surface^[34,42]. For the bulk states with inversion symmetry, TDM follows the restrictions of $D_{| |, c v} (- k) = D_{| |, c v}^{*} (k)$ and $D_{| |, c v} (- k) = \pm D_{| |, c v} (k)$ , which should be purely real or imaginary. For the surface states, however, the TDM that lacks inversion symmetry should take a complex value, $D_{| |, c v} (k) = | D_{| |, c v} (k) | e^{i θ (k)}$ ; here $θ (k)$ is the nontrivial transition dipole phase (TDP)^[28].

3. Results and Discussion

To find out how TDM affects the observed HHG modulation in our experiments, we focus on the interband polarization between conduction and valence bands in SBEs. The interband current generated by a single quantum orbital satisfies^[28] $J_{inter} (t) \propto D_{| |, c v}^{*} (k) D_{| |, c v} (k^{'}) F (t^{'}) e^{- i Ω t - i S (k, t, t^{'}) - \frac{t - t^{'}}{T_{2}}},$ (1)where $k^{'} = k + A (t) - A (t^{'})$ is the crystal momentum of electrons excited from valence band to conduction band, $k$ is the crystal momentum of electrons recombining back to the valence band, $A (t)$ is the vector potential of the laser field $F (t)$ , $S (k, t, t^{'})$ is the classical action, $Ω = N ω$ is the frequency of the $N$ th harmonic with $ω$ the fundamental frequency, and $T_{2}$ is the dephasing time. In generalized recollision, the electron-hole pair is created at a time of $t^{'}$ and accelerated in the laser field and then recollides with each other at a time of $t$ , thus, bringing the emission of high harmonics.

Considering that the SHG is very weak and only perturbs the dynamical phase of the electron wave packet, the classical action is given by^[38] $S_{2} (t, ϕ) = S_{1} (t) - \int_{t^{'}}^{t} d τ v (τ, t^{'}) \cdot A_{2} (τ, ϕ),$ (2)where $A_{2} (τ, ϕ) = A_{2} \sin (2 ω τ + ϕ)$ is the vector potential of the SHG, $S_{1}$ the action related to the fundamental pulse only, $v (k) = \nabla_{k} ε_{g} (k)$ the group velocity, and $ε_{g} (k)$ the momentum dependent bandgap. The second term in Eq. (2) is an additional dipole phase induced by the SHG, $σ (t, ϕ) = \int_{t^{'}}^{t} d τ v (τ, t^{'}) \cdot A_{2} (τ, ϕ)$ , where $t^{'}$ and $t$ are the excitation and the recombination time of electron-hole pairs, respectively.

We calculated the interband current by summing up the transient current from adjacent half-cycles, which is similar to the derivation in Ref. [28], $J_{| |, inter} (t) \propto [e^{i σ (t, ϕ)} - e^{- i N π} \cdot e^{- i σ (t, ϕ)}] \cdot D_{| |, c v}^{*} (k) D_{| |, c v} (k^{'}) F (t^{'}) e^{- i N ω t - i S (k, t, t^{'}) - \frac{t - t^{'}}{T_{2}}},$ (3)where the additional dipole phase is the origin of even harmonics, and can be written in the form of $σ (t, ϕ) = σ_{a} \cos ϕ + σ_{b} \sin ϕ = Φ (t) \cos [ϕ - α (t)]$ ^[24,38], where $Φ (t) = \sqrt{σ_{a}^{2} + σ_{b}^{2}}$ and $α (t) = a \tan (σ_{b} / σ_{a})$ . Clearly, the additional dipole phase oscillates with the relative phase delay, $ϕ$ . The amplitude of the $N$ th harmonic is proportional to $Δ_{N} = [e^{i σ (t, ϕ)} - e^{- i N π} \cdot e^{- i σ (t, ϕ)}] D_{| |, c v}^{*} (k_{s}) D_{| |, c v} (k_{s}^{'})$ . Then, the high-harmonic yield can be calculated by $S_{HHG} (Ω) \propto {| Δ_{N} |}^{2}$ . For HHG from bulk states of ${Bi}_{2} {Se}_{3}$ with real-valued TDM, the spectral intensity of odd harmonics is $S_{odd} (Ω) \propto \cos^{2} σ$ , and that of even harmonics is $S_{even} (Ω) \propto \sin^{2} σ$ . The additional dipole is far less than 1, due to the very weak SHG field. Thus, the spectral intensity of odd harmonics can be simplified as $S_{odd} \propto 1 - Φ^{2} (t) \cos^{2} [ϕ - α (t)] .$ (4)The spectral intensity of even harmonics is $S_{even} \propto Φ^{2} (t) \cos^{2} [ϕ - α (t)] .$ (5)Therefore, the modulation periods of both even and odd harmonics from bulk states are all $T_{2 ω} / 2$ .

For the surface states, the dipole phase plays an indispensable role. Because the SHG is very weak, its influence on TDM can be ignored^[24,28,38]. Taking the expression of complex TDM, we can get $Δ_{N} = {e^{i [Δ θ + σ (t, ϕ)]} - e^{- i N π} \cdot e^{- i [Δ θ + σ (t, ϕ)]}} \cdot | D_{| |, c v} (k) | \cdot | D_{| |, c v} (k^{'}) |$ , where $Δ θ = θ (k^{'}) - θ (k)$ . The spectral intensity of the odd harmonics is $S_{odd} (Ω) \propto \cos^{2} (Δ θ + σ)$ , and that of the even harmonics is $S_{even} (Ω) \propto \sin^{2} (Δ θ + σ)$ . Considering $σ$ is very small, we can get $S_{odd} \propto \cos^{2} Δ θ - \sin (2 Δ θ) \cdot σ (t, ϕ) - \cos (2 Δ θ) \cdot σ^{2} (t, ϕ),$ (6)and $S_{even} \propto \sin^{2} Δ θ + \sin (2 Δ θ) \cdot σ (t, ϕ) + \cos (2 Δ θ) \cdot σ^{2} (t, ϕ) .$ (7)

In both Eqs. (6) and (7), the first terms are contributed from the nontrivial dipole phase, and the second and third terms describe the modulation due to the SHG field that oscillates as a function of phase delay $ϕ$ . The modulation periods of HHG from the surface states are all $T_{2 ω}$ . Similar modulation has been observed in recent experiments; however, no in-depth explanation has been given^[26,27]. Our theoretical analysis further explains the relationship between the modulation of even harmonics and the symmetry restrictions. We summarize these symmetry restrictions of even harmonics in Table 1.

	Without SHG	With SHG
Inversion symmetry (bulk)	σ = 0, Δθ = 0; no even harmonics	σ ≠ 0, Δθ = 0; even harmonics; modulation, T_2ω/2
Inversion symmetry breaking (surface)	σ = 0, Δθ ≠ 0; even harmonics	σ ≠ 0, Δθ ≠ 0; even harmonics; modulation, T_2ω

Table 1. Summary of Modulation Results Dictated by Symmetry

View all Tables

The analytical expression of the intraband current is relatively complex; however, it has a form similar to that of the interband polarization^[28]. Through the above analytical analysis, we interpret the origin of the two kinds of HHG modulation.

Our analysis above focuses on the contributions of surface states and bulk states to HHG. It confirms our speculation that for very weak SHG, even harmonics with $T_{2 ω}$ modulation originate from asymmetric surface states, and with the increase of the SHG intensity, the contribution from the bulk states of ${Bi}_{2} {Se}_{3}$ begins to play a role. It is known that the contribution from surface states of TIs contains TSSs and two-dimensional electron gas (2DEG) in the process of HHG^[35]. Here we carried out the simulation by solving the SBEs with the Hamiltonian of TSSs and bulk states derived from the tight-binding model (TBM)^[43]. We calculated the HHG from TSSs and bulk states, respectively. The grid points of two-dimensional momentum grid are set to $400 \times 400$ in the length gauge. The wavelength of the fundamental is set to 3820 nm, the intensity is set to $40 GW / {cm}^{2}$ , and the wavelength of the SHG is set to 1910 nm. The dephasing times are set at $T_{1} = 200 fs$ and $T_{2} = 5 fs$ .

Figure 4(a) shows the high-harmonic spectrum from TSSs with and without the SHG in pump pulses, respectively. The intensity of H4 increases, and the intensity of H6 decreases in the presence of the SHG compared with the absence of the SHG, which is consistent with the experimental observations displayed in Fig. 2. Therefore, we speculate that the even harmonics in the experiment mainly come from TSSs. For comparison, in Fig. 4(b) we plot the high-harmonic spectrum from bulk states driven by the fundamental and two-color laser fields. Unsurprisingly, even harmonics vanish when the SHG is absent, whereas even harmonics can be produced by adding very weak SHG. These indicate that the weak SHG breaks the symmetry of the driven field and induces asymmetry in electron-hole trajectory; thus, the even harmonics are produced^[23].

Figure 4.Calculated high-harmonic spectra. (a) HHG from the surface states with and without the SHG; blue line, harmonic spectrum without the SHG; orange line, harmonic spectrum for the SHG intensity equal to 40 MW/cm² and ϕ = 0.6π. (b) HHG from the bulk states with and without the SHG; blue line, harmonic spectrum without the SHG; orange line, harmonic spectrum for the SHG intensity equal to 40 MW/cm² and ϕ = 0.6π.

We also calculated the high-harmonic yield that modulates relative to the phase delay $ϕ$ , as shown in Figs. 5(a)–5(d). For HHG from the TSSs, the modulation periods of H4 and H5 are all equal to $T_{2 ω}$ . The modulation depth of harmonic yield scales linearly with the increase of the SHG intensity [Figs. 5(a) and 5(b)]. Figures 5(c) and 5(d) show the modulation of HHG from the bulk states; however, the modulation periods of H4 are $T_{2 ω} / 2$ and the modulation depth is proportional to the SHG intensity. Besides, the modulation of H5 begins to emerge as the SHG intensity reaches $40 {MW / cm}^{2}$ [Fig. 5(d)]. As the SHG intensity further increases, the total yield of H5 decreases, which is qualitatively consistent with the experimental measurements in Fig. 3(b). Therefore, the HHG modulation in the numerical calculation verifies the analytical analysis above.

Figure 5.High-harmonic yield modulated with the phase delay of two-color laser fields. (a), (b) Modulation of H4 (a) and H5 (b) from the surface states, for the SHG intensity of 4 MW/cm² (red line) and 40 MW/cm² (blue line), respectively; (c), (d) modulation of H4 (c) and H5 (d) from the bulk states, for the SHG intensity of 4 MW/cm² (red line), 40 MW/cm² (blue line), and 400 MW/cm² (purple line), respectively.

4. Conclusion

To conclude, the experimentally observed even harmonics with the modulation period of $T_{2 ω}$ can be attributed to the coupling between the nontrivial dipole phase of the surface states and the additional dynamical phase accumulated in the acceleration of the electron-hole pairs. The $T_{2 ω} / 2$ modulation emerges in even harmonics originating from the bulk states as the SHG intensity further increases. Our findings are helpful in distinguishing the strong field-driven nonperturbative dynamics inside topologically protected surface states and bulk states. High-harmonic emission provides a novel ultrafast spectroscopy for probing the strong-field physics in a topological matter. In addition, it also gives new inspiration to manipulate the transportation of helical Dirac fermions by using a weak perturbing laser field. At the same time, we notice that the surface states of TIs contain rich physical properties such as TSSs and 2DEG, so how to distinguish their contribution to HHG is still a problem to be explored^[35,42].

References

[1] T. Popmintchev, M. C. Chen, D. Popmintchev, P. Arpin, S. Brown, S. Alisauskas, G. Andriukaitis, T. Balciunas, O. D. Mucke, A. Pugzlys, A. Baltuska, B. Shim, S. E. Schrauth, A. Gaeta, C. Hernandez-Garcia, L. Plaja, A. Becker, A. Jaron-Becker, M. M. Murnane, H. C. Kapteyn. Bright coherent ultrahigh harmonics in the keV X-ray regime from mid-infrared femtosecond lasers. Science, 336, 1287(2012).

[2] T. Popmintchev, M. C. Chen, P. Arpin, M. M. Murnane, H. C. Kapteyn. The attosecond nonlinear optics of bright coherent X-ray generation. Nat. Photonics, 4, 822(2010).

[3] Z. H. Chang, A. Rundquist, H. W. Wang, M. M. Murnane, H. C. Kapteyn. Generation of coherent soft X rays at 2.7 nm using high harmonics. Phys. Rev. Lett., 79, 2967(1997).

[4] H. J. Worner, J. B. Bertrand, D. V. Kartashov, P. B. Corkum, D. M. Villeneuve. Following a chemical reaction using high-harmonic interferometry. Nature, 466, 604(2010).

[5] J. Itatani, J. Levesque, D. Zeidler, H. Niikura, H. Pepin, J. C. Kieffer, P. B. Corkum, D. M. Villeneuve. Tomographic imaging of molecular orbitals. Nature, 432, 867(2004).

[6] C. Vozzi, M. Negro, F. Calegari, G. Sansone, M. Nisoli, S. De Silvestri, S. Stagira. Generalized molecular orbital tomography. Nat. Phys., 7, 822(2011).

[7] H. Niikura, N. Dudovich, D. M. Villeneuve, P. B. Corkum. Mapping molecular orbital symmetry on high-order harmonic generation spectrum using two-color laser fields. Phys. Rev. Lett., 105, 053003(2010).

[8] P. M. Paul, E. S. Toma, P. Breger, G. Mullot, F. Auge, P. Balcou, H. G. Muller, P. Agostini. Observation of a train of attosecond pulses from high harmonic generation. Science, 292, 1689(2001).

[9] G. Sansone, E. Benedetti, F. Calegari, C. Vozzi, L. Avaldi, R. Flammini, L. Poletto, P. Villoresi, C. Altucci, R. Velotta, S. Stagira, S. De Silvestri, M. Nisoli. Isolated single-cycle attosecond pulses. Science, 314, 443(2006).

[10] P. B. Corkum, F. Krausz. Attosecond science. Nat. Phys., 3, 381(2007).

[11] F. Krausz, M. Ivanov. Attosecond physics. Rev. Mod. Phys., 81, 163(2009).

[12] L. Gallmann, C. Cirelli, U. Keller. Attosecond science: recent highlights and future trends. Annu. Rev. Phys. Chem., 63, 447(2012).

[13] S. Ghimire, A. D. DiChiara, E. Sistrunk, P. Agostini, L. F. DiMauro, D. A. Reis. Observation of high-order harmonic generation in a bulk crystal. Nat. Phys., 7, 138(2011).

[14] O. Schubert, M. Hohenleutner, F. Langer, B. Urbanek, C. Lange, U. Huttner, D. Golde, T. Meier, M. Kira, S. W. Koch, R. Huber. Sub-cycle control of terahertz high-harmonic generation by dynamical Bloch oscillations. Nat. Photonics, 8, 119(2014).

[15] B. Zaks, R. B. Liu, M. S. Sherwin. Experimental observation of electron-hole recollisions. Nature, 483, 580(2012).

[16] G. Vampa, T. J. Hammond, N. Thire, B. E. Schmidt, F. Legare, C. R. McDonald, T. Brabec, D. D. Klug, P. B. Corkum. All-optical reconstruction of crystal band structure. Phys. Rev. Lett., 115, 193603(2015).

[17] S. Ghimire, D. A. Reis. High-harmonic generation from solids. Nat. Phys., 15, 10(2019).

[18] J. Y. Cao, N. Li, Y. Bai, P. Liu, R. X. Li. Even-order high-harmonic generation from solids in velocity gauge. Chin. Opt. Lett., 19, 043201(2021).

[19] C. X. Zhang, X. H. Li, E. Chen, H. R. Liu, P. P. Shum, X. H. Chen. Hydrazone organics with third-order nonlinear optical effect for femtosecond pulse generation and control in the L-band. Opt. Laser Technol., 151, 108016(2022).

[20] A. F. Kemper, B. Moritz, J. K. Freericks, T. P. Devereaux. Theoretical description of high-order harmonic generation in solids. New J. Phys., 15, 023003(2013).

[21] T. Higuchi, M. I. Stockman, P. Hommelhoff. Strong-field perspective on high-harmonic radiation from bulk solids. Phys. Rev. Lett., 113, 213901(2014).

[22] P. G. Hawkins, M. Y. Ivanov, V. S. Yakovlev. Effect of multiple conduction bands on high-harmonic emission from dielectrics. Phys. Rev. A, 91, 013405(2015).

[23] G. Vampa, C. R. McDonald, G. Orlando, D. D. Klug, P. B. Corkum, T. Brabec. Theoretical analysis of high-harmonic generation in solids. Phys. Rev. Lett., 113, 073901(2014).

[24] G. Vampa, T. J. Hammond, N. Thire, B. E. Schmidt, F. Legare, C. R. McDonald, T. Brabec, P. B. Corkum. Linking high harmonics from gases and solids. Nature, 522, 462(2015).

[25] T. T. Luu, H. J. Worner. High-order harmonic generation in solids: a unifying approach. Phys. Rev. B, 94, 115164(2016).

[26] T. T. Luu, H. J. Worner. Observing broken inversion symmetry in solids using two-color high-order harmonic spectroscopy. Phys. Rev. A, 98, 041802(R)(2018).

[27] G. Vampa, H. Z. Liu, T. F. Heinz, D. A. Reis. Disentangling interface and bulk contributions to high-harmonic emission from solids. Optica, 6, 553(2019).

[28] S. C. Jiang, H. Wei, J. G. Chen, C. Yu, R. F. Lu, C. D. Lin. Effect of transition dipole phase on high-order-harmonic generation in solid materials. Phys. Rev. A, 96, 053850(2017).

[29] S. C. Jiang, J. G. Chen, H. Wei, C. Yu, R. F. Lu, C. D. Lin. Role of the transition dipole amplitude and phase on the generation of odd and even high-order harmonics in crystals. Phys. Rev. Lett., 120, 253201(2018).

[30] D. Golde, T. Meier, S. W. Koch. High harmonics generated in semiconductor nanostructures by the coupled dynamics of optical inter- and intraband excitations. Phys. Rev. B, 77, 075330(2008).

[31] J. Y. Cao, F. S. Li, Y. Bai, P. Liu, R. X. Li. Inter-half-cycle spectral interference in high-order harmonic generation from monolayer MoS2. Opt. Express, 29, 4830(2021).

[32] A. J. Uzan-Narovlansky, A. Jimenez-Galan, G. Orenstein, R. E. F. Silva, T. Arusi-Parpar, S. Shames, B. D. Bruner, B. H. Yan, O. Smirnova, M. Ivanov, N. Dudovich. Observation of light-driven band structure via multiband high-harmonic spectroscopy. Nat. Photonics, 16, 428(2022).

[33] D. Hsieh, Y. Xia, D. Qian, L. Wray, J. H. Dil, F. Meier, J. Osterwalder, L. Patthey, J. G. Checkelsky, N. P. Ong, A. V. Fedorov, H. Lin, A. Bansil, D. Grauer, Y. S. Hor, R. J. Cava, M. Z. Hasan. A tunable topological insulator in the spin helical Dirac transport regime. Nature, 460, 1101(2009).

[34] H. J. Zhang, C. X. Liu, X. L. Qi, X. Dai, Z. Fang, S. C. Zhang. Topological insulators in Bi2Se3, Bi2Te3 and Sb2Te3 with a single Dirac cone on the surface. Nat. Phys., 5, 438(2009).

[35] Y. Bai, F. C. Fei, S. Wang, N. Li, X. L. Li, F. Q. Song, R. X. Li, Z. Z. Xu, P. Liu. High-harmonic generation from topological surface states. Nat. Phys., 17, 311(2021).

[36] C. P. Schmid, L. Weigl, P. Grossing, V. Junk, C. Gorini, S. Schlauderer, S. Ito, M. Meierhofer, N. Hofmann, D. Afanasiev, J. Crewse, K. A. Kokh, O. E. Tereshchenko, J. Gudde, F. Evers, J. Wilhelm, K. Richter, U. Hofer, R. Huber. Tunable non-integer high-harmonic generation in a topological insulator. Nature, 593, 385(2021).

[37] D. Baykusheva, A. Chacon, J. Lu, T. P. Bailey, J. A. Sobota, H. Soifer, P. S. Kirchmann, C. Rotundu, C. Uher, T. F. Heinz, D. A. Reis, S. Ghimire. All-optical probe of three-dimensional topological insulators based on high-harmonic generation by circularly polarized laser fields. Nano Lett., 21, 8970(2021).

[38] N. Dudovich, O. Smirnova, J. Levesque, Y. Mairesse, M. Y. Ivanov, D. M. Villeneuve, P. B. Corkum. Measuring and controlling the birth of attosecond XUV pulses. Nat. Phys., 2, 781(2006).

[39] P. B. Corkum. Plasma perspective on strong-field multiphoton ionization. Phys. Rev. Lett., 71, 1994(1993).

[40] Y. Bai, C. J. Cheng, X. L. Li, P. Liu, R. X. Li, Z. Z. Xu. Intense broadband mid-infrared pulses of 280 MV/cm for supercontinuum generation in gaseous medium. Opt. Lett., 43, 667(2018).

[41] X. He, J. M. Dahlstrom, R. Rakowski, C. M. Heyl, A. Persson, J. Mauritsson, A. L’Huillier. Interference effects in two-color high-order harmonic generation. Phys. Rev. A, 82, 033410(2010).

[42] D. Hsieh, J. W. McIver, D. H. Torchinsky, D. R. Gardner, Y. S. Lee, N. Gedik. Nonlinear optical probe of tunable surface electrons on a topological insulator. Phys. Rev. Lett., 106, 057401(2011).

[43] D. Baykusheva, A. Chacon, D. Kim, D. E. Kim, D. A. Reis, S. Ghimire. Strong-field physics in three-dimensional topological insulators. Phys. Rev. A, 103, 023101(2021).