{'fraction_non_alphanumeric': 0.029401579166089487, 'fraction_numerical': 0.014060119130073417, 'mean_word_length': 4.674395755551189, 'pattern_counts': {'":': 0, '<': 2, '<?xml version=': 0, '>': 1, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Light control of surface–bulk coupling by terahertz vibrational coherence in a topological insulator

The demand for disorder-tolerant quantum logic and spin electronics can be met by generating and controlling dissipationless spin currents protected by topology. Dirac fermions with helical spin-locking surface transport offer a way of achieving such a goal. Yet, surface-bulk coupling can lead to strong Dirac electron scattering with bulk carriers and phonons as well as impurities, assisted by such dissipative channel, which results in “topological breakdown”. Here, we demonstrate that coherent lattice vibrations periodically driven by a single-cycle terahertz (THz) pulse can significantly suppress such dissipative channel in topological insulators. This is achieved by reducing the phase space in the bulk available for Dirac fermion scattering into during coherent lattice oscillations in Bi2Se3. This light-induced suppression manifests as a remarkable transition exclusively in surface transport, absent for bulk, above the THz electric fields for driving coherent phonons, which prolongs the surface transport lifetime. These results, together with simulations, identify the critical role of spin–orbit coupling for the “phase space contraction” mechanism that suppresses the surface-bulk coupling. Imposing vibrational quantum coherence into topological states of matter may become a universal light control principle for reinforcing the symmetry-protected helical transport.


INTRODUCTION
Recently, there is increasing evidence that coherently driven states can protect and enhance the stability of quantum systems. A pedagogical example for such dynamic stabilization is the periodically driven Kapitza pendulum 1 where an inverted, yet stable, orientation can be achieved by a sufficiently high frequency vibration of its pivotal point. Applied to topological protected states (TPS), such as those realized in topological insulators (TIs), 2-7 coherent lattice vibrations can have direct and profound effects on surface transport of Dirac fermions, via periodic modulation of electronic states, as illustrated in Fig. 1a. This effect of coherent phonons is in stark contrast to that of incoherent phonons from thermal excitations, which diminishes topological enhancement of surface transport via random lattice scattering, heating, surface-bulk charge transfer and coupling. Ultrafast phononics has been explored as a new avenue to manipulate properties of superconductors, 8 oxides, 9,10 semimetal 11 and photovoltaic semiconductors. 12 Although experimental realizations of coherent phonon generation have been explored in TIs, 7,13,14 the influence of such vibrational coherence on the surface helical transport has not been observed until now.
The prototype TI, such as Bi 2 Se 3 exhibits surface charge transport protected by spin-momentum locking, where the spin-orbit interaction connects the direction of the carrier spin to that of its momentum. 4,[15][16][17] This spin-locking gives rise to coherent transport that manifests as, e.g., lack of current relaxation and prolonged transport lifetime, which is protected against nonmagnetic disorder by symmetry. Here, we explore a nonequilibrium driven TPS, periodically modulated by coherent phonons induced by a noncontact THz push electric field, as a new control knob for further protecting surface carriers from scattering with residual bulk states and/or thermal phonons. Such THz photon driven coherence mechanism complements the spontaneous symmetry protection mechanism in the ground state. However, its observation requires intense THz pump and THz conductivity probe spectroscopy measurements [18][19][20][21][22] of the coherent phonon-dressed TPS, which has never been demonstrated so far.
In this article, we create a periodically driven topologically protected state via coherent lattice vibrations that significantly suppress the surface-bulk coupling in Bi 2 Se 3 . We demonstrate two regimes of post-THz-pump nonequilibrium states, obtained by tuning E THz . In this way we show evidence of nonthermal quantum control of the Dirac fermion transport lifetime τ Dirac , by inducing a transition from incoherent carrier heating to a coherent driving regime. This transition exclusively appears in surface transport and is absent in bulk. It occurs above electric field strengths required for observing coherent phonon oscillations. Together with theoretical simulations, our results reveal a "phase space contraction" mechanism assisted by spin-orbit coupling that is responsible for coherent phonon-enabled suppression of the surfacebulk coupling in dressed topological states out of equilibrium.

RESULTS
We apply a single-cycle THz pulse E THz with amplitude up tõ 500 kV cm −1 , which selectively excites coherent oscillations of E 1 u phonons, shown in Fig. 1b. For E THz = 0, the real part of the static conductivity, σ 1 (ω), at the ground state shows a Drude lineshape superimposed to a well-defined oscillator feature centered at 8 meV as shown in Fig. 1c. This lineshape has been established as arising from joint carrier conductivity responses from both bulk and surface states plus a E 1 u phonon mode, consistent with prior THz electrodynamic studies of Bi 2 Se 3 . 23  E THz ≠ 0, the spectral-temporal characteristics of the complex conductivity response functions are measured by our THz pump-THz probe conductivity spectroscopy that allows us to differentiate surface from bulk contributions as shown in ref. 23 . In this way we measure the transport lifetimes of Dirac currents, τ Dirac , and bulk ones, τ Bulk , which reveals the interaction between topological states and bulk coherent phonons.
We study 50 quintuple layers of n-type Bi 2 Se 3 with bulk insulating bandgap~310 meV and Fermi energy~60 meV above the conduction band edge. This Bi 2 Se 3 film is grown by molecular beam epitaxy on top of a 0.5 mm thick sapphire substrate with caxis perpendicular to surface. The sample is cooled down to 5 K in a He flow cryostat in order to maximize the surface responses by reducing thermal phonons. Figure 1d shows phonon emission after THz pump at 5 K as a function of electric field E THz , up to~500 kV cm −1 . The broadband THz pulses used here are centered at ω pump~4 meV (Fig. 1c) and can excite the E 1 u phonon mode. Interestingly, the phonon emission clearly shows a pronounced multi-cycle oscillation for high driving fields~180 kV cm −1 . The Fourier transform (FT) spectra of these coherent oscillations displays a resonant peak centered at~8 meV, which is consistent with the E 1 u phonon mode (inset, Fig. 1e). The integrated spectral weight (SW) clearly shows a "threshold-like" build-up at E THz~1 80 kV cm −1 (Fig. 1e). Below E th 180 kV cm −1 , the coherent phonon THz emission SW is vanishingly small. This nonlinear emergence of the coherent phonon oscillations is most likely due to the multi-photon excitation of coherent E 1 u phonons~8 meV. The THz pump spectra is broad, from 2 to 6 meV with center~4 meV (Fig. 1c). Therefore, the excitation will be a combination of multiple processes that mainly include two-and three-photon nonlinear processes. This is consistent with the fitted data (Supplementary Note 6), which shows a power law with exponent~2.7 in the E-field dependence. Most intriguingly, this transition in lattice displacement driven outof-equilibrium distinguishes two different THz pumping regimes controlled by the THz electric field: (1) weak pumping can mainly excite incoherent phonons, averaging to zero lattice movement, and mainly accesses a lattice heating regime; (2) intense THz driving allows the establishment of a macroscopic lattice vibration and accesses a periodically driven TPS. The latter regime gives rise to the observed THz radiation, as shown in Fig. 1d. Note that THz pulses will also excite charge carriers by intraband transitions at surface and bulk because of finite doping.
We characterize the ultrafast THz conductivity by measuring the THz fields in the time domain, shown in the inset of Fig. 2a, b. In the inset, we show E(t) after transmission through (i) a reference substrate, E ref (t) (gray) and (ii) the unexcited sample, E Sample (t) (blue). We also compare the above to the pump-induced change ΔE(t) after pump-probe delay Δt pp = 0.5 ps (red, inset). Using FT and Fresnel equations, the full THz dielectric response is obtained and expressed as the real parts of the conductivity, Δσ 1 (ω), Fig. 2a, and the dielectric function, Δε 1 (ω), Fig. 2b. These response functions characterize the change in transport properties of surface Dirac electrons and bulk carriers following THz driving. We observe a clear bipolar frequency-dependence of Δσ 1 (ω), with transition from strong bleaching (Δσ 1 < 0) to enhanced absorption (Δσ 1 > 0) for frequencies between 2 and 10 meV. We also observe a decrease of Δε 1 with increasing probe frequency in the same range. In our previous study of intraband photoexcitation responses in Bi 2 Se 3 , mainly using mid-infrared pulses, 23 we identified a similar bipolar spectral shape of Δσ 1 (ω) dynamics, which arises from the combined responses due to competing surface and bulk states. There we showed that the above behavior can be used to distinguish the intrinsic transport properties of Fig. 1 Vibrational quantum coherence in Bi 2 Se 3 driven by an intense single-cycle THz pulse above a threshold electric field. a Experimental schematics of the THz pump process in Bi 2 Se 3 film. Intense THz pulses drive coherent phonon oscillations of the E 1 u mode that strongly influence surface state transport. b THz electric field profile in the time domain, with peak field as large as 500 kV cm −1 and c broadband spectrum from 1 to 10 meV (pink shade). The static conductivity σ 1 (ω), shown together in c, exhibits Drude transport behavior superimposed to a bulk E 1 u phonon mode resonance at~8 meV (black solid circle). d The post-pump THz emission for various driving THz electric fields demonstrates a threshold-like behavior clearly observed above E th~1 80 kV cm −1 . e Pump electric field dependence of the coherent phonon integrated spectral weight shows a significant increase above threshold E th . Inset: Corresponding Fourier spectrum of timedomain THz emission driven by single-cycle THz pump with E-field of 500 kV cm −1 .
Dirac fermions and bulk electrons. The lineshape can also be explained by quantitative simulation as applied below.
The observed spectral lineshapes as a function of pump-probe time delay can be fitted by using a composite model with two Drude terms and one Lorentz oscillator. These three contributions account for the THz response functions of surface and bulk carriers and phonons, respectivelỹ where ε ∞ is the background electrical permittivity. The plasma frequencies ω 2 p ¼ ne 2 =ε 0 m Ã are proportional to the densities of surface (n S ) and bulk (n B ) free carriers. The subscripts j = S and B in the equation stand for surface and bulk states, respectively. ω, e, and ε 0 denote the frequency, electron charge, and vacuum electrical permittivity, respectively. m* and γ are the effective mass and scattering rates of surface or bulk carriers. The THz response functions measured in the experiment are well reproduced for all time delays by varying the surface and bulk scattering rates γ j , the phonon oscillator strength F, and the phonon damping coefficient Γ. At the same time, we keep the surface and bulk carrier densities fixed, consistent with a physical picture of sub-gap photoexcitation that is decoupled from interband transitions. 23 A unique set of model parameters is obtained by simultaneously fitting conductivity and permittivity in the measured broad frequency range. By only varying the surface τ S and bulk τ B scattering times (i.e., 1/γ), we are able to consistently reproduce the time-dependent experimental results, as shown in Fig. 2a, b, and extract the Dirac (purple line) and bulk (green line) transport properties. This fitting implies that the THz pumping used here gives rise to much smaller changes of surface and bulk quasiparticle densities as compared to the scattering time changes, unlike for high frequency pumping as discussed in ref. 23 . We can thus separate the intrinsic surface and bulk transport properties in both time and frequency domains.
The distinct spectral-temporal characteristics observed here can be used to characterize how THz driving influences carrier transport. In Fig. 2c-f, complete 2D false color plots of Δσ 1 and Δε 1 spectra are presented as a function of THz pump electric field E THz up to 500 kV cm −1 , at short pump-probe delay t pp = 0.5 ps after excitation. The pump-induced THz conductivity spectra Δσ 1 (ω) exhibit a bipolar behavior, demonstrated by the two THz field-cut spectra in Fig. 2c taken at the black dashed lines in Fig. 2d. We observe a pump field-dependent nonlinear negative growth and then decay of Δσ 1 (ω) at 5 K. As shown in Fig. 2, the low frequency negative signals from 2 to 4 meV become dominant for driving fields above~180 kV cm −1 and then progressively decrease with increasing THz pump field. To analyze this behavior, we fit the 2D THz response functions Δσ 1 (ω) and Δε 1 (ω) shown in Fig. 2d, f by using the composite model discussed above. Excellent agreement between this model and the measured transient THz response is achieved for all THz pump fluences, similar to the exemplary fitting in Fig. 2a, b. As we discuss next, these results indicate that the transport properties of Dirac fermions arise from two competing effects: coherent phonon periodic modulation and incoherent phonon heating of the Fermi sea.  180 kV cm −1 , there is a reduction in τ S with increasing field, marked by the black dashed line, which we assign to increased surface electron scattering due to THz heating of the Fermi sea and incoherent thermal phonons. Above~180 kV cm −1 , this reduction in τ S clearly slows down, which indicates suppression of Dirac fermion scattering compared to what is expected for thermal effects. The observation of a clear transition between two excitation regimes indicates the emergence of additional protection absent in the ground state. This stops further suppression of τ S even with increasing THz field pumping, which finally saturates at the highest THz field used. The most salient feature is the coincidence between the transition field observed in the Dirac fermion transport lifetime changes and the nonlinear threshold for emergence of coherent phonon oscillations at E th . This coincidence suggests a strong correlation between the two phenomena, i.e, the emergence of prolonged, new nonequilibrium transport lifetimes when the THz field drives ordered lattice positions via coherent phonon oscillations, instead of disordered lattice positions due to thermal fluctuations. At the same time, the THz pump field dependence of the bulk transport lifetime τ B (inset, Fig. 3a) does not show any transition across E th (Fig. 3b). Instead, τ B displays an opposite field dependence, i.e., an increase in τ B at low pump fields followed by a decrease at high pump fields. This behavior can be understood as a transition from impurity-to phonon-dominated scattering with increasing THz pumping, which is consistent with conventional behaviors seen in the temperature dependence of bulk carrier scattering in a doped semiconductor. 24,25 Figure 3c illustrates our interpretation of these experimental observations as signatures of suppression of the surface-bulk coupling by THz driving, which leads to a nonequilibrium topological state protected and stabilized further by a new dynamic channel from coherent phonon excitations.

DISCUSSION
To explore how coherent lattice displacements driven by E THz influence the topological electronic state, we use first-principles calculations in combination with a microscopic theory of surfacebulk coupling as discussed in details in the Supplementary Note 3. Bi 2 Se 3 has a rhombohedral structure (space group R-3m) with five atoms in the primitive unit cell (Fig. 4a). These atoms are arranged in a quintuple layered structure, with each layer consisting of five planes of Se2-Bi-Se1-Bi-Se2, with lattice parameters a = 4.142 Å and c = 28.635 Å. The details of our first-principles calculation are discussed in the Supplementary Notes 3 and 4. The most important scattering mechanisms of the surface electrons in Bi 2 Se 3 at low temperature are electron-electron (e-e) and electron-impurity (e-imp) scattering. On the one hand, the e-e scattering rate Γ e−e for normal Fermi-liquids is proportional to T 2 . In our THz pump experiment, the surface electronic temperature increases due to intraband scattering after THz excitation, therefore we expect that Γ eÀe / ðE THz Þ 4 leads to reduced carrier transport lifetimes when heating dominates. On the other hand, the e-imp scattering in the surface is weakly temperaturedependent. Next we explore how the e-imp scattering, assisted by surface-bulk coupling, can be reduced nonthermally by intense THz driving of coherent phonons that will compete with Γ e−e above the threshold as demonstrated in Fig. 3a.
The e-imp scattering can proceed via surface-surface (SS) and surface-bulk (SB) channels, whose contribution depends on the available final states, i.e., the surface electron may scatter to either surface or bulk states. According to the microscopic theory of impurity-and phonon-induced SB coupling in TIs, 26 in n-doped Bi 2 Se 3 samples with Fermi level at~60 meV above the conduction band minimum (CBM) similar to ours, the SB scattering channel dominates over the SS one. The SB impurity scattering rate is shown to be proportional to the bulk density of states (DOS) at the Fermi level and is also sensitive to the relative orientation of the surface and bulk orbital pseudospin 〈τ z 〉. 26 Using a model Hamiltonian derived to describe the low-energy electronic structure of Bi 2 Se 3 up to O(k 2 ), 27 we computed the low energy electronic bands and corresponding surface and bulk expectation values of the orbital pseudospin 〈τ z 〉. The results are plotted in Fig.  4b. It can be seen that, at 60 meV above the CBM, the surface and bulk states have parallel orbital pseudospin orientation. Thereby, the SB scattering rate is enhanced and dominates over other channels, consistent with ref. 26 . Since jhτ z i S j 1, jhτ z i B j ' 1, the SB scattering rate of the surface electron at the Fermi level near the conduction band edge takes the following simplified form 26 (Supplementary Note 4) where the summation includes all bulk bands α 0 . g imp α 0 is proportional to the effective electron-impurity coupling and is dependent on the overlap between the surface and the α 0 bulk band wavefunctions. v α 0 ðE F Þ is the DOS of the bulk band α 0 at the Fermi level. The SB scattering rate Γ imp SB is therefore proportional to the bulk DOS at the Fermi level.
Next we analyze the variation of the bulk DOS at the Fermi level, v α 0 ðE F Þ, during the lattice vibration under the E 1 u mode. To determine the lattice dynamics, we calculate the phonon properties with a supercell approach and harmonic approximation using the Phonopy code. 28 In the calculations, a supercell containing 135 atoms, i.e., 3 × 3 × 3 of the primitive unit cell, and a k-point grid of 2 × 2 × 2 are used. The force constants are computed without and with spin-orbit coupling (SOC) using the VASP code. 29,30 Our results show that the frequency of the E 1 u mode is 2.57 THz from the GGA calculation, while the GGA + SOC calculation gives 2.30 THz, closer to the experimental measurement (~2.0 THz). The directions of the lattice vibrations under the E 1 u modes are schematically plotted in Fig. 4c, where the two inequivalent E 1 u modes are illustrated from different views. Note that a hexagonal plane is used to represent all atoms in the quintuple layer without losing generality. The atomic displacement of each layer is determined to be r = (−0.  Fig. 4d, e. Here, the bulk DOS at the Fermi level v α 0 (E F ) is plotted as a function of the lattice vibration amplitude (in units of r) from GGA and GGA + SOC calculations, respectively, i.e., without (Fig. 4d) and with SOC (Fig. 4e). We note that the results obtained for the two inequivalent E 1 u modes are not the same, thus an average curve is also computed and plotted in Fig. 4d, e (black solid circles).
The DOS dependence on phonon amplitude is markedly different with or without SOC, which identifies its critical role for determining the phase space available for the SB scattering rate Γ imp SB . Figure 4d shows that v α 0 (E F ) gradually increases in the regular GGA calculations without SOC, which implies increased Dirac electron scattering. In contrast, the calculations in Fig. 4e with SOC show that v α 0 (E F ) rapidly decreases for initial decrease of lattice vibrations (<0.2 Å from the equilibrium position), and then slowly increases and saturates. The Dirac fermion-impurity scattering rate Γ imp SB , Eq. (2), changes similarly to v α 0 (E F ) as a function of the lattice vibration amplitude, which is proportional to the THz pump field. The simulated behavior clearly explains the slow down in τ S decrease with this new surface protection channel against bulk scattering and saturation at the highest THz field used (Fig. 3a). Therefore, the good agreement between the calculated results when SOC is included and the experimental observations clearly underpins the importance of spin-orbit coupling for the "phase space contraction" mechanism that suppresses the surface-bulk coupling in Bi 2 Se 3 during coherent phonon oscillations. Note that the observed nonequilibrium τ S in Fig. 3a (black circles) is always smaller than the unexcited one mainly due to the unavoidable heating of bulk Fermi sea in highly doped samples like ours. However, the "phase space contraction" mechanism demonstrated in Fig. 4 indicates that further increase in τ S can be achieved, e.g., in samples with the doping below the bulk band edge to remove heating effects of bulk carriers. Our data for the bulk conducting state clearly shows the roadmap and it will be interesting to extend the demonstrated coherent-driving scheme to bulk insulating samples and other systems including semimetals and semiconductors to increase coherent transport. 11,12,31 We expect two major changes with elevated temperature: (1) decrease of coherent phonon amplitude for the same pump field and (2) much stronger surface-bulk coupling, assisted by the population of incoherent phonons. Figure 5 shows these expected Fig. 4 First-principles calculations, carried out by using DFT and microscopic scattering theory, demonstrate the proposed "phase space contraction" mechanism that suppresses the surface-bulk coupling. a Crystal structure of Bi 2 Se 3 . b Low-energy electronic band structure of the surface (dash lines) and bulk (solid lines) states in Bi 2 Se 3 , calculated using a model Hamiltonian. The color map indicates the value of the orbital pseudospin 〈τ z 〉. The Fermi level sits at 60 meV above the conduction band minimum (CBM) and is indicated by the black dotted line. c Schematic representation of the two E 1 u modes calculated in Bi 2 Se 3 at q = 0. The vibration directions of the two inequivalent E 1 u modes are explained with a hexagonal plane that generally represents all five planes in the Bi 2 Se 3 quintuple layer. d, e The variation of the bulk DOS at E F as a function of lattice vibration amplitude (in units of the atomic displacement r) is calculated from d GGA and e GGA + SOC.
behaviors. First, Fig. 5a, b compare THz pump-induced phonon emission at 5 and 300 K for two selected pump fields, above and below E th~1 80 kV cm −1 . One clearly sees coherent phonon oscillation at low temperature above E th , as discussed above. In contrast, the phonon oscillation strength is significantly suppressed at high temperature, due in part to enhanced scattering and thermal population. This result is consistent with our conclusion that the incoherent phonons dominate the THzdriven state at high temperature, unlike for low temperature. This interpretation naturally leads to the conclusion that strong surface-bulk coupling and charge transfer exist at high temperatures as discussed next. Second, THz spectra of Δσ 1 (ω) and Δε 1 (ω) measured at 300 K (red circles) are shown in Fig. 5c, d, which could not be fitted with our composite model (green solid lines) consisting of two Drude contributions with only the scattering rate changes. This is in contrast to the 5 K data, which are fitted very well with the proposed model. This clearly shows that surface and bulk are strongly coupled by the incoherent phonon with charge transfer in this high temperature case, 32 which can only be accounted beyond a linear combination of their individual responses.
There are distinguishing features and robust analysis that show the presence of two transport channels, i.e., both surface and bulk, and reliably determine their transport lifetimes, as shown in Fig. 3a. We have addressed them extensively in our prior THz conductivity study using ultra-broadband, wavelength-selective pumping that extends from mid-IR to near-IR and visible. 23 Applied to the THz pump data, for example, one distinguishing feature is the observation of frequency-dependent cooling times that directly distinguish the faster surface channel from the slower bulk one as shown in the Supplementary Note 5. This cannot be explained by a single state with fixed carrier density. The only physical picture that one can propose to consistently explain this unusual behavior is that there are two relaxation channels, i.e., surface and bulk, with different cooling times. 23 In summery, intense THz-driven coherent phonon modulations provide a previously inaccessible tuning knob for suppressing surface-bulk scattering channels. Our results provide compelling implications of exploring "phononic engineering" of phase space, by modulation of electronic DOS, toward achieving persistent coherent transport in topological device applications.

Experimental schemes
Broadband quasi-single-cycle THz pulses were generated by optical rectification in a MgO-doped LiNbO 3 crystal using the tilted pulse front method. The crystal was pumped by 2 W and 35 fs pulses from a Ti:Sapphire based regenerative amplifier with a 1 kHz repetition rate. A portion of the amplifier beam was used to generate THz probe pulses via optical rectification in a 110-cut 1 mm ZnTe crystal. The probe THz waveforms were measured using electro-optic sampling in another 110-cut 1 mm ZnTe crystal. The pump THz electric field was varied using three wire-grid polarizers placed in the THz beam path. The pump THz and probe THz spot sizes were 1.2 and 0.8 mm, respectively. A 2 mm aperture was placed in front of the sample to ensure uniform illumination of the probed area. The sample was mounted in an optical cryostat along with a bare substrate which served as a reference. More information can be found in Supplementary Note 1.

Sample preparation
The Bi 2 Se 3 thin film sample, 50 quintuple layers (~50 nm thick), is grown by molecular beam epitaxy on a 0.5 mm thick sapphire substrate. The sample is mounted together with a 0.5 mm thick pure sapphire substrate, used as a reference, into a cryostat with temperatures down to T = 5 K. Two copper mounts with the same aperture are placed directly in front of the sample and reference to ensure uniform photoexcitation and accurate comparison of their THz transmission. The sample studied has a Fermi energy E F 60 meV from the bulk conduction band edge, as estimated from the measured plasma frequency in the static THz conductivity spectra. Please note that the exact value of the gap size is not critical for our scheme, nor does it change our conclusion, as long as the interband transition gap in our sample is larger than our pump photon energy~4 meV. Other details can be found in the Supplementary Note 2.

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

Fig. 5
Enhanced surface-bulk coupling and suppression of the coherent phonon oscillations at high temperature. Comparison of THz phonon emission at a 5 K and b 300 K after two selected pump fields below and above E th~1 80 kV cm −1 . This clearly reveals strong coherent phonon emission at low temperature above E th , while the coherent oscillations are suppressed at 300 K at the THz pump fields used. Pump-induced THz spectra of c Δσ 1 (ω) and d Δε 1 (ω) measured at 300 K. The experimental results cannot be fitted by the composite model (green lines) with fixed surface and bulk densities that works well for low temperature data in Fig. 2a, b.