Open access peer-reviewed chapter

Ultrafast Dynamics in Topological Insulators

By Phuoc Huu Le and Chih-Wei Luo

Submitted: November 3rd 2017Reviewed: February 5th 2018Published: March 2nd 2018

DOI: 10.5772/intechopen.74918

Downloaded: 1051


Ultrafast dynamics of carriers and phonons in topological insulator Bi2Se3, CuxBi2Se3 (x = 0, 0.1, 0.125) single crystals were studied by time-resolved pump-probe spectroscopy. The coherent optical phonon (A1g1) is found via the damped oscillation in the transient reflectivity changes (∆R/R) for CuxBi2Se3. The observed red shift of A1g1 phonon frequency suggests the intercalation of Cu atoms between a pair of the quintuple layers of Bi2Se3 crystals. Moreover, the relaxation processes of Dirac fermion near the Dirac point of Bi2Se3 are studied by optical pump and mid-infrared probe spectroscopy through analyzing the negative peak of the ∆R/R. The Dirac fermion-phonon coupling strength was found in the range of 0.08–0.19 and the strength is reduced as it gets closer to the Dirac point. The ultrafast dynamics and fundamental parameters revealed by time-resolved pump-probe spectroscopy are important for designing the optoelectronics in the mid-IR and THz ranges.


  • topological insulators
  • ultrafast dynamics
  • pump-probe spectroscopy
  • Dirac fermion

1. Introduction

Recently, topological insulators (TIs) [1, 2, 3, 4, 5, 6, 7, 8] and two-dimensional (2D) materials such as graphene [9], MoS2, WS2, and MoSe2 [10] are of great interests because of their unique physical properties and applications. These materials have a band structure that is linearly dispersed with respect to momentum, in which the transportation of electrons in these materials is essentially governed by Dirac’s (relativistic) equation with zero rest mass and an effective “speed of light”—c* ≈ 106 m/s [9]. In TIs, a novel electronic state called the topological surface state (TSS) has been predicted and observed [1, 2, 3, 4, 5, 6, 7, 8]. Unlike the trivial insulator, TIs have a spin degenerate and fully gapped bulk state but exhibit a spin polarized and gapless electronic state on the surface [8]. This metallic surface state has a linear energy-momentum dispersion relation in the low-energy region, which is known as a Dirac cone. Unlike the Dirac cone of graphene, the Dirac cone of a TI is protected by the time-reversal symmetry. This robust TSS can survive under time-reversal invariant perturbations, such as surface pollution, crystalline defects, and distortions of the surface [6]. Additionally, because of the fully spin-polarized characteristics of the surface state, TIs have a high potential for the development of spintronic devices and quantum computation [6, 11].

The optoelectronic properties of TIs are important subjects for the development of optoelectronic devices. Therefore, the issues associated with electron–phonon interaction, carrier lifetime, carrier dynamics, energy loss rate, and low-energy electronic responses are very important for optimizing device performance. These ultrafast dynamic properties of the materials can be resolved by pump-probe spectroscopy. This chapter provides a brief introduction to the materials, time-resolved pump-probe spectroscopy, and some ultrafast dynamic properties of Bi-based topological insulators.


2. Bismuth-based topological insulators

Bismuth chalcogenide compounds (Bi2Ch3, Ch = Se, Te) have been extensively investigated in material science and condensed-matter physics because of their intriguing properties regarding thermoelectricity [12, 13, 14] and three-dimensional TIs [15, 16, 17, 18]. Bi2Ch3 is a narrow bandgap semiconductor with a rhombohedral crystal structure belonging to the D3d5R3¯mspace group. The Bi2Ch3 crystal structure is constructed from repeated quintuple layers (QLs) arranged along the c-axis. The unit lattice cell of a Bi2Ch3 crystal is composed of three QLs. Each QL is stacked in a sequence of atomic layers Ch(1)-Bi-Ch(2)-Bi-Ch(1) and is weakly bonded to the next QL via Van der Waals interaction. The crystal structures of Bi2Se3 and Bi2Te3 are shown in Figure 1. For convenience, these crystal structures are also described by a hexagonal lattice, where the a-axis and c-axis lattice constants of Bi2Se3 (Bi2Te3) are 4.138 Å (4.384 Å) and 28.64 Å (30.487 Å), respectively [19].

Figure 1.

(a) The crystal structures of (a) Bi2Se3 and (b) Bi2Te3, in which 5-atomic-layer-thick lamellae of –(Se(1)–Bi–Se(2)–Bi–Se(1))– or –(Te(1)–Bi–Te(2)–Bi–Te(1))– is called a quintuple layer (QL).

In 2009, Zhang et al. predicted that the Bi2Ch3 crystal is a strong TI [15]. A calculation of the electronic structure with spin-orbit coupling in the Bi2Se3 crystal has also been performed [15]. By tuning the spin-orbit coupling in the system, band inversion occurred around the Γ point. As these two levels, which are closest to the Fermi energy, have opposite parity, the inversion between them drives the system into a TI phase [15]. Figure 2 shows the calculated energy and momentum dependence of the local density of states (LDOS) for Sb2Se3, Sb2Te3, Bi2Se3, and Bi2Te3. All of these materials have the same rhombohedral crystal structure with the space group D3d5R3¯m. Zhang et al. predicted that Bi2Se3, Bi2Te3, and Sb2Se3 are candidates for a TI, whereas Sb2Se3 is not because the spin-orbit coupling effect of Sb2Se3 is not strong enough to induce band inversion [15]. Following this prediction, Xia et al. [20] and Hsieh et al. [4] investigated the existence of the TSS in Bi2Se3, Bi2Te3, and Sb2Se3 through angle-resolved photoemission spectroscopy (ARPES).

Figure 2.

The calculated energy and momentum dependence of the LDOS for (a) Sb2Se3, (b) Sb2Te3, (c) Bi2Se3, and (d) Bi2Te3 on the (111) surface. The TSSs are clearly seen around the Γ point as a red line in the Sb2Te3, Bi2Se3, and Bi2Se3 graphs. No TSS exists in Sb2Se3 [15].

Figure 3(a) and (b) shows the ARPES results of the surface electronic structure on a Bi2Se3 (111) surface [20]. Around the Γ¯point, the clear V-shaped band is observed to approach the Fermi level. The slopes of this V-shaped band along the Γ¯Μ¯and Γ¯K¯directions are nearly equivalent [20]. The U-shaped bands near the Fermi level and below the V-shaped band are the bulk conduction band (BCB) and bulk valence band (BVB) of Bi2Se3. This result matches the prediction that the surface state exists between the BCB and the BVB. The detail of the surface state is shown in Figure 3(d). The ring-like Fermi surface formed by the Dirac cone-like surface state is centered at the Γpoint. The unique spin-momentum lock behavior can also be observed in this figure.

Figure 3.

The ARPES measurements of the electronic structure on Bi2Se3 near theΓ¯point along the (a)Γ¯Μ¯and (b)Γ¯K¯directions. (c) The bulk 3D Brillouin zone and the surface 2D Brillouin zone of the projected (111) surface. (d) The Fermi surface of the surface state [20].


3. Principle of femtosecond spectroscopy

3.1. Degenerate pump-probe spectroscopy

Highly temporal resolution is one of the unique characteristics in femtosecond optics. By the pump-probe technique, the photoexcited carrier dynamics and phonon dynamics in solid state materials can be clearly resolved. Additionally, the interband and intraband relaxation processes can be also obtained.

The basic understanding of time-resolved pump-probe spectroscopy is introduced as follows. The pump pulses are served as a perturbation which leads to the changes of the electronic population in materials. The probe pulses are used for the detection of the optical property changes of the materials. By controlling the time interval between the pump and probe pulses, the transient changes of the optical properties can be recorded. In pump-probe spectroscopies, the transient reflectivity changes (∆R/R) or transient absorption changes (∆A/A) can be measured.

Here, we explain more experimental details about the detection of ∆R/R. As shown in Figure 4, the pump-induced reflectivity changes are plotted as R(t). The fmodulation is the modulation frequency of the chopper. The flaser is the repetition rate of femtosecond laser pluses. The period of the R(t) is correlating to the period of the Ipump(t) pulse train. Since the reflectivity of material, R(t), is modulated by the Ipump(t) pulse train, the intensity of reflective probe beam IR,probe(t) is also modulated by the Ipump(t) pulse train. Thus, the IR,probe(t) can be described by the superposition of the DC intensity signal IR,0 and AC intensity signal ΔIRwith the specific frequency f(= fmodulation). Typically, the ΔIRis much smaller than the IR,0 in the order of 10−3–10−6. By using the phase-lock technique, the amplitude of AC intensity signal ΔIRcan be extracted out by providing the reference frequency fmodulation for the lock-in amplifier. Because the ΔIRis small compared to the IR,0, the IR,probet=IR,0+ΔIR/2IR,0. Thus, the IR,0can be obtained by using a multimeter for the measurements of IR,probet. The measured ∆R/Rcan directly relate to the ΔIR/IR,0via the following relationship.


Figure 4.

The scheme for the principle of the pump-probe technique. The details are described in the text.

where the Ii,probeis the intensity of the incident probe beam. Thus, the time evolution of ΔRt/Rcan be measured by swapping the time interval Δtbetween the pump pulse and the probe pulse.

3.2. Optical pump and mid-infrared probe spectroscopy

The plasma edge of the doped n-type semiconductor usually lies in the mid-infrared (MIR) regime. By measuring the reflectivity around the plasma edge, many characteristics of carriers such as scattering rate and carrier concentration can be obtained [21]. The development of a pulsed mid-infrared light source provides the opportunities for understanding the dynamics of carriers. The mid-infrared pump-probe spectroscopy has been already applied on various materials (i.e., oxides, semiconductors, superconductors, graphene, and topological insulators) [22, 23, 24, 25, 26, 27, 28]. In the reflection-type mid-infrared pump-probe spectroscopy, the effect of multiple reflections should be considered in the analysis, and the dynamical characteristics of carriers can be further obtained through modeling the measured data with the Drude-Lorentz model.

Figure 5 shows a schematic diagram of our optical pump and mid-infrared probe (OPMP) spectroscopy. The light source of the pump-probe system is a regenerative amplifier with 800 nm central wavelength, 5 kHz repetition rate, and 30 fs pulse duration. The beam is split into a pump beam (40% of the incident light) and a probe beam (60% of the incident light). The probe beam passes through a 0.7-mm-thick GaSe crystal to generate mid-infrared (MIR) pulses, in which the MIR wavelength can be tuned from 9.0 μm (138 meV) to 14.1 μm (88 meV) through differential frequency generation (DFG). The optical pump beam with the fluence of 68 μJ/cm2 and a spot size of 485 μm (in diameter) is focused on the sample using a 150 mm lens. An Au-coated off-axis parabolic mirror with f = 200 mm is employed to focus the probe beam on the sample surface with a spot diameter of 392 μm. It is ensured that the spot size of the pump beam is larger than that of the probe beam. The probe beam is further collimated and refocused onto a MIR detector (e.g., liquid nitrogen-cooled HgCdTe) using an Au-coated off-axis parabolic mirror (f = 50 mm).

Figure 5.

Schematic diagram of the optical pump and mid-infrared probe (OPMP) system.


4. Ultrafast dynamics in topological insulators

4.1. Time-resolved spectroscopy in a topological insulators

The dynamic properties of photoexcited TIs have attracted a great deal of attention. For example, the relaxation behavior of a carrier near the Fermi surface has been observed by the time-resolved angle-resolved photoemission spectroscopy (Tr-ARPES) [29, 30, 31, 32]. Figure 6(c) shows that the 1.55 eV photons excite the electrons from the bulk valence band to a higher-lying state in the bulk materials. Then, the photoexcited carriers fall into the bulk conduction band (BCB) and the surface state within 1 ps [31]. In Figure 6(a), we can see the rise time of curve 10 is ~1 ps. This means that after photoexcitation, the carriers in the higher lying band are rapidly relaxed into the BCB, then cooled to the bottom of the BCB via intraband scattering. These interband transitions and intraband scattering are shown in Figure 6(d) and (e) [31].

Figure 6.

(a) The normalized population within the integration windows indicated in (b). (b) The integration windows over the BCB and the surface state. (c)–(f) Schematics of carrier dynamics over the transition energy range [31]. (c)t = 0 ps. (d)t~ 0.5 ps. (e)t~ 2.5 ps. (f)t > 5 ps.

Furthermore, the relaxation time of curve 10 in Figure 6(a) is longer than 10 ps. This slow relaxation indicates the metastable behavior of the population of carriers in the BCB [30, 31]. Meanwhile, as curves 6–9 shown in Figure 6(b), the population of surface states also exhibits an unusually long-lived existence [31]. Here, the relaxation bottleneck is attributed to the scattering processes between the BCB and the surface state [31]. As Figure 6(f) shows, the photoexcited carriers first relax via surface-bulk scattering and then cooling via surface-state intraband scattering. This scattering channel is mainly in response to the acoustic phonon-mediated surface-bulk coupling and the acoustic phonon scattering of the surface-state Dirac fermions [32]. The Tr-ARPES can directly deliver information about the population changes of the electronic state near the Fermi level. However, reports on the transition processes occurring in the early stages after photoexcitation are rare. To fully understand the photoexcited carrier dynamics, studies for the interband transition and the intraband cooling are needed, which can be revealed using optical pump/optical probe spectroscopy (OPOP) and optical-pump/mid-infrared probe (OPMP) spectroscopy.

4.2. Interband relaxations in topological insulators

The interband relaxation of photoexcited carriers in topological insulator (TI) single crystals is examined by the optical pump and optical probe spectroscopy [33]. In this section, we present the phonon and carrier dynamics in doped TI CuxBi2Se3 (x = 0, 0.1, 0.125) single crystals. Figure 7(a) shows the typical ΔR/Rsignals as a function of delay time for CuxBi2Se3 crystals at room temperature. Generally, different energy-transfer processes can be unambiguously extracted from the time evolution of ΔR/Rcurve. After pumping, the thermalization between electrons and optical phonons which occurred in a sub-picosecond timescale is characterized by the fast component in ΔR/R. A subsequent slow component in a timescale of several picoseconds is assigned to the thermalization between electrons and acoustic phonons [34]. After these electron-lattice relaxation processes, the heat diffusion out of the illuminated area on the sample is further revealed by the quasi-constant component in ΔR/R[35]. Furthermore, all of the ΔR/Rcurves show two damped oscillation components with different periods.

Figure 7.

Temporal variations of ΔR/Rsignals for CuxBi2Se3 crystals (x = 0, 0.1, 0.125) at room temperature by using the 1.55 eV degenerate pump-probe spectroscopy, shown (a) in the full timescale and (b) in short timescale.

The slow oscillation components, as shown in Figure 7(a), are attributed to the coherent acoustic phonons (CAPs) generated by ultrafast laser pulses. This damped slow oscillation in ΔR/Ris generated by the interference between two probe beams, respectively, reflected from the sample surface and the strain pulse that propagate longitudinally with the sound velocity. The relationship between the period τCAPof the slow oscillation and the longitudinal sound velocity vsis τCAP=λ/vsn2sin2θ, where λis the probe wavelength, nis the refractive index at λ, and θis the incident angle of the probe beam [33]. Consequently, the sound velocity can be estimated by measuring the CAP oscillations when the refractive index of the material is known. The frequency of the CAP for the Bi2Se3 crystals is ~0.033 THz (~ 30 ps in period). Additionally, it completely decays within ~60 ps. The disappearance of the CAP (slow oscillation) around 60 ps, according to the strain pulse model, is determined by the penetration depth of an 800 nm probe beam in Bi2Se3 crystals. Taking the refractive index of Bi2Se3 crystals reported in [36], the sound velocity is estimated to be 1996 m/s at room temperature [33]. Figure 7(a) also reveals that the periods of the slow oscillations in CuxBi2Se3 (x = 0, 0.1, 0.125) crystals vary slightly from 29.9 to 30.2 ps.

The fast oscillation components of CuxBi2Se3 crystals are presented in Figure 7(b), which can be extracted by removing the relaxation background from the ΔR/Rsignals. The results are presented in Figure 8(a). The frequency of the component is 2.148 THz, which can be further assigned as the A1g1 coherent optical phonon (COP) mode of Bi2Se3, based on comparison with the steady-state Raman spectroscopy [37]. Interestingly, the frequencies of the fast oscillations considerably vary with Cu content (x) of the CuxBi2Se3 samples and are associated with the changes in the chain length of the QL and in the lattice constant of c-axis.

Figure 8.

(a) High-frequency temporal variations of ΔR/Rsignals for CuxBi2Se3 (x = 0, 0.1, 0.125) crystals extracted fromFigure 7(b). (b) Phonon frequency (red squares), dephasing time of A1g1 phonon mode (blue circles) and changes in lattice constant ofc-axis (green triangles) in CuxBi2Se3 crystals as a function of Cu doping concentrations.

Figure 8(a) shows the fast oscillation component for CuxBi2Se3 (x = 0, 0.1, 0.125) crystals. In order to quantitatively analyze these oscillations, a damped oscillation function, Aosccos2πfosct+ϕet/τdephasing, was used to fit the original data in Figure 8(a) to get the dephasing time (τdephasing) and the phonon frequency (fosc) for the CuxBi2Se3 crystals. As shown in Figure 8(b), both dephasing time and phonon frequency shrink as Cu concentrations increase (x). This indicates that an additional Cu atom deforms the Se-Bi-Se-Bi-Se chain in CuxBi2Se3 crystals. Furthermore, the lattice constant of c-axis increases slightly with increasing Cu concentrations [Figure 8(b)], implying that the QL chain in CuxBi2Se3 is stretched by introducing Cu atoms. Thus, the scenario of stretching the QL chain length is that the Cu atoms (form a mediated layer) are intercalated between QLs to strengthen the interaction between QLs. Moreover, the QLs are further deformed by these intercalated Cu atoms.

4.3. Intraband relaxations in topological insulators

The femtosecond snapshots of the relaxation processes and Dirac fermion-phonon coupling strength of 3D TI Bi2Se3 were revealed by OPMP spectroscopy [26]. In this study, several selected Bi2Se3 single crystals with a wide range of carrier concentrations (n) from 51.5 × 1018 to 0.25 × 1018 cm−3 were studied. Table 1 summarizes the doping levels of samples (#1: n = 51.5 × 1018 cm−3, #2: n = 13.9 × 1018 cm−3, #3: n = 5.58 × 1018 cm−3, and #4: n = 0.25 × 1018 cm−3). The OPMP spectra and the corresponding ARPES images of the samples are shown in Figure 9(a) and (b) [26]. The OPMP spectra clearly show a positive ΔR/Rpeak for high n ≥ 13.9 × 1018 cm−3 (Bi2Se3 #1 and #2). In contrast, this positive peak gradually diminishes as ndecreases, while an additional negative peak appears for the cases of n = 5.58 × 1018 cm−3 and n = 0.25 × 1018 cm−3.

CodeEF-EDirac point
Carrier concentrationnsurface /(nsurface+ nbulkd)
(1018 cm−3)
(1013 cm−2)
#1422−51.5 ± 0.84−1.450.11
#2325−13.9 ± 0.26−0.830.20
#3284−5.58 ± 0.25−0.720.35
#4260−0.25 ± 0.01−0.470.89

Table 1.

Fermi energy and carrier concentrations of bulk and surface states of various Bi2Se3 single crystals. All samples are n-type. “d = 23.5 nm” is the penetration depth of 800 nm pumping light.

Figure 9.

Carrier concentration (n) dependence of the transient changes in reflectivity ΔR/Rin Bi2Se3 single crystals. (a) ΔR/Rof samples #1 (n = 51.5×1018 cm−3), #2 (n = 13.9×1018 cm−3), #3 (n = 5.58×1018 cm−3), and #4 (n = 0.25×1018 cm−3) with a fluence of 34 μJ/cm2 for pumping and a photon energy of 141 meV for probing. (b) ARPES images on samples of (a) [26].

Based on the ARPES image and the energy band structure of TI Bi2Se3, a model is proposed [in Figure 10(a)] for the optical pumping (1.55 eV) and mid-infrared probing processes to elucidate the origins of both positive and negative signals. The band gap of Bi2Se3 is ~300 meV, as shown in the ARPES images of Figure 9(b), which is much larger than the probe photon energy (87~153 meV) of the mid-infrared (mid-IR). Thus, it does not allow the occurring of the interband transitions between the valence band (VB) and the conduction band (CB) of the bulk. Meanwhile, the free-carrier absorption in the CB [the probe (1) in Figure 10(a)] and Dirac cone surface state [the probe (2) in Figure 10(a)] will dominate the probe processes, which can be assigned to the positive and negative peaks in ΔR/R, respectively. To reveal the physical meanings of the positive peak in ΔR/R, the photon energy dependence of ΔR/Rfor #1 sample is studied and shown in Figure 10(b). Clearly, ΔR/Rgradually changes from positive to negative as decreasing the photon energy. At around 136 meV (1100 cm−1), it appears that intermediate signals mixed with both positive and negative peaks, corresponding to deep in the Fourier transform infrared (FTIR) reflectance spectrum [the inset of Figure 10(b)]. The excited carriers after pumping suffer the so-called intervalley scattering, leading to the red shift of the reflectance spectra. Thus, the reflectivity increases as a function of time when probing photon energy is higher than the position of the 136 meV deep. In contrast, the reflectivity decreases as a function of time when probing photon energy is smaller than the 136 meV deep. Similar results were also observed in a typical semiconductor n-type GaAs [23].

Figure 10.

Schematic energy band structure and photon energy-dependent ΔR/Rin a bulk state. (a) Schematic band structure of TIs based on the ARPES images in Figure 9(b) and the pump-probe processes. CB: Conduction band, VB: Valence band, SS: Surface state.Ri,Dirac &Rf,Dirac: The circumferences of initial/final states in Dirac cone for probing. (b) With a fluence of 38 μJ/cm2 for pumping, the ΔR/Rof Bi2Se3 #1 at various photon energies (wavenumber) from 87 to 153 meV (700–1234 cm−1). Inset: The Fourier transform infrared (FTIR) reflectance spectrum of Bi2Se3 #1. The gray area indicates the range of the mid-IR photon energy used in this study [26].

As found in Figure 9 and Table 1, the amplitude of positive peak in ΔR/Rgradually decreases as bulk carrier concentrations reduce. Meanwhile, the negative peak of ΔR/Rincreases while reducing the bulk and surface carrier concentrations. Intriguingly, the negative peak increases substantially with an increasing ratio of the surface to total carrier concentration [nsurface /(nsurface + nbulkd) in Table 1], implying a close relation between the negative peak of ΔR/Rand Dirac fermions. In addition, Figure 11(a) shows the ΔR/Rsignal as a function of the pumping fluences. The positive peak exhibits a stronger dependence on the pumping fluences than the negative peak does. For a pumping fluence of 3.3 μJ/cm2, the maximum photo-induced carrier density Δnis around 2.54 × 1018 cm−3. Indeed, if one absorbed photon generates one photo-induced carrier, the maximum photo-induced carrier density can be estimated by Δn = (1−RF/(E×δ), where R = 0.55 is the reflectance, F = 3.3 μJ/cm2 is pumping fluence, E = 2.48 × 10−19 J (= 1.55 eV) is the pumping photon energy, δ = 23.5 nm is the penetration depth. Consequently, the negative peak still subsists at the low pumping fluence of 3.3 μJ/cm, while the positive peak almost vanishes [see Figure 11(b)]. Namely, the process (1) associated with the positive peak can be suppressed and the process (2) associated with the negative peak can be preserved by reducing the pumping fluences. To quantitatively certify the relation between the negative peak and Dirac fermions, the amplitude of the negative peak dependence of probing photon energy is studied using low nsamples #3 and #4 to avoid disturbance of the positive peak [Figure 11(b)]. According to Fermi’s golden rule, the amplitude of the negative peak should be proportional to the transition probability (Tif) between the initial and final density of states in the Dirac cone. Indeed, in Figure 11(c), the Ri,Dirac × Rf,Dirac presents linear relation with probing photon energy, reflecting the proportional relation to the transition rate between the initial and final density of states for the mid-IR probe process (2) in the Dirac cone [Figure 10(a)]. This confirms that the negative peak of ΔR/Ris predominantly attributed to the mid-IR probe process (2) in the Dirac cone. Consequently, the ultrafast dynamics of the Dirac fermions can be clearly disclosed by the negative peak of ΔR/R.

Figure 11.

Pumping fluence and photon energy dependence of ΔR/Rand its amplitude and rising (decay) time in the surface state. (a) With probing photon energy of 141 meV, the ΔR/Rof Bi2Se3 #4 at various pumping fluences from 3.3 to 105 μJ/cm2. (b) With pumping fluence of 3.3 μJ/cm2, the ΔR/Rof Bi2Se3 #4 at various photon energies from 90 to 152 meV. (c) The photon energy dependence of negative peak amplitude of ΔR/Rin (b). The photon energy-dependent normalized absorption probability [dashed line, i.e.Ri,Dirac ×Rf,Dirac inFigure 10 (a)] of the mid-IR probe beam in the Dirac cone. (d) The photon energy dependence of the rising time (τr) and decay time (τd) of ΔR/Rin (b) [26].

As shown in Figure 11(d), both rising time (τr) and decay time (τd) of the negative peak of ΔR/Rstrongly depend on the probing photon energy. The τr becomes longer when the probed regime is closer to the Dirac point (or smaller probing energy). The ultrafast relaxation picture for Dirac fermions in TIs can be established. The major process right after the 1.55 eV pumping is that the carriers in the bulk valence band (BVB) are excited to the bulk conduction band (BCB). The carrier recombination between the BCB and BVB can be ignored in this study because of the large timescale (typically >> 1 ns) for such a process. Consequently, the unoccupied states in BVB would mainly be refilled by carriers in the upper Dirac cone. Carriers in the Dirac cone can be easily transferred into the unoccupied states in BVB due to the overlapping between the Dirac cone and BVB [see Figure 9(b)], leading to the increase in the number of the unoccupied states near the Dirac point and thus enhancing the absorption channel for process (2) in the Dirac cone [Figure 10(a)]. Therefore, the reflectivity of the mid-IR probing light decreases within 1.47~3.60 ps, that is, the rising time of the negative peak in Figure 11(b) and (d). Once the carriers in the Dirac cone relax into BVB, the excited carriers in the BCB are subsequently injected into the unoccupied states in the Dirac cone to diminish the absorption channel for the mid-IR process (2) [Figure 10(a)] and consequently lead to the increased mid-IR reflectivity within 14.8~87.2 ps. The timescale (τd) of this process is several tens of picoseconds, which is much longer than the τr of several picoseconds, because the carriers in BCB cannot directly transfer into the top of the Dirac cone without overlaps occurring between them and other auxiliaries, for example, phonons.

The relaxation of Dirac fermions has been demonstrated via phonon medium [38, 39]. The coupling strength (λ) between Dirac fermions and phonons varies at different positions of the Dirac cone, which can be revealed from the photon energy-dependent rising time. Based on the second moment of the Eliashberg function [40], the coupling strength (λ) is inversely proportional to the relaxation time (τe) of excited electrons:


where ωis the phonon energy which couples with the electrons. For the estimate of ω2, some vibrational modes are more efficiently coupled to Dirac fermions than the others. For Bi2Se3, the symmetric A1g1mode of ~8.9 meV is coherently excited by photoexcitation and efficiently coupled. Taking τe = τr in Figure 11(d) and Te = 370 K (obtained from [32] at the aforementioned low pumping fluence) to estimate the coefficient of (πkBTe/3) in Eq. (2), photon energy dependence of the Dirac fermion-phonon coupling strength is λ= 0.08–0.19. The Dirac fermion-phonon coupling strength measured by the present OPMP becomes significantly smaller near the Dirac point (the point of K//=0). As getting closer to the Dirac point, Dirac fermions will possess a weaker coupling with the phonons to reduce the scatterings with phonons. In addition, the effective mass of Dirac fermions in the surface state gradually decreases when approaching the Dirac point, which is consistent with the results in graphene [9].


5. Conclusion

We report the ultrafast dynamics of carriers and phonons in topological insulator Bi2Se3, CuxBi2Se3 (x = 0, 0.1, 0.125) single crystals. By time-resolved pump-probe spectroscopy, one damped fast oscillation was clearly observed in the transient reflectivity changes (∆R/R) for CuxBi2Se3, which is assigned to the coherent optical phonon (A1g1). The frequency of A1g1phonon decreases considerably with increasing Cu contents, suggesting the intercalation of Cu atoms between quintuple layers of Bi2Se3. The schematic illustration of the direct transitions and subsequent relaxation processes induced by optical excitation in Bi2Se3 single crystals is also reported here. The femtosecond snapshots of the relaxation processes were revealed by optical pump and mid-infrared probe spectroscopy. Especially, the Dirac fermion dynamics in the Dirac cone surface state near the Dirac point of Bi2Se3 was unambiguously revealed through the negative peak of ∆R/R. The Dirac fermion-phonon coupling strength was found in range of 0.08–0.19 and the strength is reduced as getting closer to the Dirac point. These results are extremely crucial to the design of Dirac fermion devices and optoelectronics, especially in the mid-IR and THz ranges.



Financial support from Vietnam National Foundation for Science and Technology Development (NAFOSTED) under grant number 103.99–2015.17, the Ministry of Science and Technology of the Republic of China, Taiwan (Grant No. 103-2628-M-009-002-MY3, 103-2119-M-009-004-MY3, 106-2119-M-009-013-FS, 106-2628-M-009-003-MY3) and the Grant MOE ATU Program at NCTU are gratefully acknowledged.

© 2018 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Phuoc Huu Le and Chih-Wei Luo (March 2nd 2018). Ultrafast Dynamics in Topological Insulators, Two-dimensional Materials for Photodetector, Pramoda Kumar Nayak, IntechOpen, DOI: 10.5772/intechopen.74918. Available from:

chapter statistics

1051total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Graphene, Transition Metal Dichalcogenides, and Perovskite Photodetectors

By Zhi Yang, Jinjuan Dou and Minqiang Wang

Related Book

First chapter

Synthesis Strategies about 2D Materials

By Jianghao Wang, Guangshe Li and Liping Li

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us