Driving ultrafast spin and energy modulation in quantum well states via photo-induced electric fields

The future of modern optoelectronics and spintronic devices relies on our ability to control the spin and charge degrees of freedom at ultrafast timescales. Rashba spin-split quantum well states, 2D states that develop at the surface of strong spin-orbit coupling materials, are ideal given the tunability of their energy and spin states. So far, however, most studies have only demonstrated such control in a static way. In this study, we demonstrate control of the spin and energy degrees of freedom of surface quantum well states on Bi 2 Se 3 at picosecond timescales. By means of a focused laser pulse, we modulate the band-bending, producing picosecond time-varying electric ﬁ elds at the material ’ s surface, thereby reversibly modulating the quantum well spectrum and Rashba effect. Moreover, we uncover a dynamic quasi-Fermi level, dependent on the Lifshitz transition of the second quantum well band bottom. These results open a pathway for light-driven spintronic devices with ultrafast switching of electronic phases, and offer the interesting prospect to extend this ultrafast photo-gating technique to a broader host of 2D materials.


INTRODUCTION
Quantum well (QW) states in semiconductors are a result of a confining potential that traps electrons (typically in a 2D plane) quantizing their electronic states. These states are easy to produce, have tunable band gaps, and are conceptually simple in terms of quantum systems. For these reasons, they have yielded a plethora of applications such as high-efficiency solar cells, infrared lasers, and qubits [1][2][3][4] . Suitable confinement potentials can be achieved by constructing a multilayer heterostructure, but they can also occur spontaneously on surfaces from an adsorbateinduced localized surface potential [5][6][7] . When these states appear on the surface, the breaking of inversion symmetry removes the spin-degeneracy resulting in spin-polarized QW states [8][9][10][11][12][13] . Combining quantum confinement with strong spin-orbit coupling joins the malleable electronic dispersions of QW states with the robust spin-polarization of Rashba states, making them extremely promising candidates for spin-based devices, such as the Datta Das spin-transistor 14,15 .
Since the Rashba interaction and QW spectrum are both predicated on the underlying electric potential, tuning the surface potential becomes a powerful method to manipulate both their spin and energy degrees of freedom 7,[16][17][18] . Such control has been achieved via in-situ surface chemical doping (irreversible and difficult to implement in devices) 9,12,19,20 and electrostatic gating (ideal for devices but often limited by the magnitude of the QW potential at the surface of bulk samples, 10 7 −10 13 V m −1 ) [21][22][23] . The real downfall of both methods is in the difficulty to dynamically tune the QW states on ultrafast timescales, restricting the ability to study the fundamental mechanisms of electronic relaxation and limiting device application in the time domain.
Here, by combining the surface photovoltage effect (SPV), a well-studied mechanism for producing reversible surface potentials in many semiconductors, with time and angle-resolved photoemission spectroscopy, we demonstrate an ultrafast and reversible control of Rashba splitting and energy level spacing, and subsequent quasi-equilibrium charge storage in correspondence with a photo-induced Lifshitz transition. This is in contrast to proposed methods of spin-current manipulation in topological insulators using the SPV to solely change the local surface potential of the Dirac surface state [24][25][26] . These results present an all-in-one approach for light-driven spin-orbit modulation and charge manipulation by directly controlling the density of states on picosecond timescale. Experiments were performed on p-type Bi 2 Se 3 , a promising candidate for spin-based devices 14,15 , and ideal since it exhibits both a strong SPV effect in the bulk 27,28 and 2D QW Rashba surface states, with the latter readily produced through various surface dopants such as carbon monoxide, water, potassium, rubidium, or even long exposure to vacuum [11][12][13]20,[29][30][31][32] .

RESULTS
Gating surface QWs with light Figure 1 presents a schematic of the experimental setup and the effect of infrared light pulses on the band-bending and QW state. In a semiconductor, the SPV effect occurs when in-gap surface states lead to a redistribution of charge between the bulk and surface at equilibrium. This non-uniform charge distribution is called the space charge region and produces an internal electric field that 'bends' the bands near the surface, as schematically shown in Fig. 1a. In the case of downward band-bending, illumination of the surface via an infrared pulse generates photoexcited electrons and holes that drift to the surface and bulk, respectively. These carriers partially cancel the internal electric field thereby reducing the band-bending and altering the surface potential (solid yellow region) 33,34 . Equivalently, the downward band-bending potential results in an internal electric field at equilibrium, (E BB ), and this electric field is reduced by an opposing electric field generated by the separated photoexcited charges. The consequence of this is a time-varying electric field in the bulk that decays according to the recombination of separated photoexcited electron holes. Hence, an SPV semiconductor substrate can serve as a time-dependent electric gate to 2D electronic states at the surface, and the infrared pulse can be focused to gate a specific area of interest. The Bi 2 Se 3 samples have metallic surfaces due to a topological surface state and host QW states which can be modified by the pump-induced timedependent charge separation. In this way, the surface of the SPV material itself is the subject of the electric field gating and is modified by the time-varying internal field. An added consequence of the photo-induced charge separation is that the resulting finite dipole at the surface also produces an external electric field between the sample and the detector, E vac , that affects the detected kinetic energy of the photoelectron, resulting in a rigid energy shift of the measured dispersion, allowing us to track the SPV in real-time 27,[35][36][37][38] .
The Bi 2 Se 3 samples in this study are p-type doped with a holelike majority carrier concentration of ≈1 × 10 19 cm −3 (see Supplementary Fig. 3), yielding a fairly conductive bulk that generates a narrow SCR at the surface of~10 nm (see discussion). In panel b, we have modeled the effect of an optical pump on the bandbending potential and QW energy levels. The lower plot shows the distribution of fixed charges that form the equilibrium bandbending potential (black line) and separated charges following photoexcitation (red gradient). Photoexcited electrons are distributed equally in three parts between the surface and two QWs (using the QW wavefunctions to determine the density profile) and an equal number of photoexcited holes piled at the edge of the SCR. The band-bending length-scale is similar to the wavefunction spatial extension near the surface, such that only sufficiently separated charges on either side of the SCR can avoid high recombination rates in order to form a persisting SPV effect 39 . The amplitudes of the excited charges with respect to the fixed charges have been exaggerated by 2× for visibility. The effect of the population of separated photoexcited charges on the QW potential is seen in the upper plot (gray to orange). A secondary well forms at the bulk edge of the SCR (see dotted curve in panel b) that limits the ability to trap holes in the bulk such that the QW potential is not quenched after SPV saturation at the highest fluence (saturation is shown later in Fig. 5b). The squared modulus of the wavefunctions and corresponding energies have been calculated for the equilibrium and maximally pumped condition shown as black and dark orange curves, demonstrating compression of the energy eigenvalues. Figure 1c presents ∂ 2 /∂E 2 ARPES spectra demonstrating the compression of the QWs after pump excitation as also schematically summarized in Fig. 1d. The strong downward band bending at equilibrium leaves the surface n-type doped, filling a large portion of the upper Dirac cone, and the quantized conduction band 10,12,27 , despite the p-doping in the bulk. The first and second QWs, QW1, and QW2 are derived from the bulk conduction band (gray parabolic region in panel d) and exhibit near-parabolic dispersions. The strong spin-orbit coupling inherent to Bi 2 Se 3 as well as inversion symmetry breaking at the cleave-plane generates a Rashba effect that duplicates the QW states into momentumseparated spin-polarized pairs (shown in red and blue in panel d) [8][9][10] . Furthermore, the momentum splitting transforms the band minimum from point-like to a one-dimensional ring, resulting in a van Hove singularity below the intersection of the Rashba bands 40 for both QW1 and QW2 (see DOS in panel d). At equilibrium, with Fig. 1 Tuning quantum well potentials through charge separation. a Diagram of the band-bending in Bi 2 Se 3 inside and outside the illuminated region. Infrared excitation generates excited electron-hole pairs that are separated by the band-bending, reducing the internal electric field. b Photoexcited charge separation decreases the slope of the quantum well potential (gray to orange region). The quantum well energy levels (n = 1, 2) before illumination (solid curves) are compressed in energy after the pump (dotted curves). c ∂ 2 /∂E 2 spectra at equilibrium and after pump pulse demonstrating the spectrum modification. d Illustration of quantum well level shifts in Bi 2 Se 3 before and after illumination. The gray parabolic region represents the conduction band before quantization. The DOS of Rashba-split bands exhibits a singularity near the band bottom. The photo-induced spectrum modification enables toggling of a Lifshitz transition.
no pump beam, only QW1 is populated, exhibiting strong Rashba splitting. After pumping with a 51.8 μJ cm −2 820 nm pulse, QW1 and QW2 shift downwards in energy relative to the bulk conduction band minimum. QW1 shifts downward by~25 meV and QW2 emerges below E F and is subsequently populated resulting in a sudden increase of the DOS. Therefore, the lightdriven modification and subsequent filling of QW2 drives a van Hove singularity below the Fermi level, inducing a Lifshitz transition 41 . Figure 2 demonstrates the evolution of the QW spectrum after ultrafast illumination. Since the SPV decays in time and the effective gating field depends on the SPV strength, a single timeresolved measurement records a continuum of the field effect on the QW. Figure 2a illustrates the field effect on surface-confined QW as a function of delay. The gray region denotes the QW SCR confinement at the surface and the white region represents the near-surface bulk. At equilibrium (I), the band-bending is maximally resulting in an enhanced slope for the QW potential and can be represented by an internal electric field E BB . The maximum reduction in the band-bending is the potential difference between the surface and deep within the bulk at equilibrium, labeled V BB,0 . At τ = 0 (II), when the pump pulse arrives, separated photoexcited carriers generate an opposing electric field to E BB , labeled E eff . The internal field from the separated carriers acts as an effective applied field on the QW potential, driving the time-dependent modification to the spectrum. In the illustrated scenario, the pump pulse is significantly strong as to saturate the available hole trap sites altering the surface potential as a function of delay by ΔV SPV (τ) = V BB,0 − V BB (τ), resulting in reduced but non-zero band-bending. At more positive delays (III-IV), the excited charge recombines thus reducing the effective field E eff , reviving the band-bending potential, and bringing the surface potential back to its equilibrium value.

Dynamic QW spectrum and Lifshitz transition
We can monitor the photoexcited charge separation through the energy shift of the ARPES spectrum induced by the in vacuum field produced by ΔV SPV . The change in the surface potential has the effect of increasing the photoelectron kinetic energy at the detector with respect to the equilibrium energy position (ΔKE). Figure 2b presents the effect of the SPV rigid shift on the topological surface state at the Fermi wavevector as a function of delay. The rigid shift of the spectrum extracted from the shift of the Dirac dispersion at fixed momenta well below E F is plotted in white. The Fermi edge as determined by fitting a Lorentzian multiplied by a Fermi-Dirac distribution is plotted in red and denotes the quasi-Fermi level, E Ã F (See Supplementary Fig. 2). The rigid energy shift is a direct measure of the surface potential as a function of delay since ΔKE(τ) =−eΔV SPV (τ), where e is the charge of the electron.
Under the center of charge approach 42,43 , the SPV signal can be described by an average displacement of electron and hole sheet charges induced by the band banding such that: where n sep and d are the charge separation density and separation distance, respectively, and ϵ r ϵ 0 is the dielectric constant of Bi 2 Se 3 .
With the assumption that the separation distance is fixed and approximately equal to the width of the band-bending region, then the measurement of ΔV SPV is proportional to the charge separation concentration n sep , and therefore the charge at the surface 42,44 . To be clear, this is different than the total photoexcited carriers generated by the pump since significant recombination has taken place after 10 ps, with only sufficiently separated charges surviving (see Supplementary Fig. 5). Moreover, if the charge is sufficiently separated and localized (a fit assumption based on the long lifetime of the SPV), then the photo-induced charge separation is analogous to a parallel plate capacitor within the probe beam spot and the effective applied field is also proportional to ΔV SPV and the surface charge. Therefore, the ΔKE shift (panel b for τ > 0) decays proportionally to the separated charge density (panel b for τ > 0), enabling us to map the delay axis to a surface voltage and ultimately an applied field strength (discussed in Fig. 3).
At negatives delays (photoemission occurs before the pump) there is a rigid shift of the recorded ARPES spectrum. This effect is Fig. 2 Time-resolved electric field gating of quantum well spectrum. a Effective applied field, E eff is generated by separated photoexcited carriers, manifesting as a reduction of the band-bending potential. SPV is maximized before quenching the band-bending due to a finite limit on trapped holes. As photoexcited charge recombines, band-bending returns, and the effective field is reduced. b Separated surface charge can be monitored by the rigid energy shift of the ARPES spectrum due to the voltage change at the surface (ΔV SPV ). Black-lined white dots correspond to the SPV-induced rigid shift of the band structure relative to equilibrium E F , red dots correspond to the quasi-Fermi level at k F of the TSS. c DOS at the Fermi level as a function of delay revealing the Lifshitz transition of QW2. d Symmetrized ARPES spectra at equilibrium and three non-equilibrium delays demonstrating the pump-induced quantum well shifts after correcting the SPV shift. e Graphic summarizing quantum well spectrum dynamics corresponding to the spectra in d, presenting a delay-dependent field effect on the quantum wells and long-lived quasi-E F .
due to the photoelectron interacting with an SPV-generated field induced by a pump pulse that arrives after photoemission [35][36][37][38] . The rigid energy shift in the negative delay is therefore shifting the unpumped spectrum and can be ignored for the purposes of this study.
In the positive delay, a maximum ΔKE shift of 72.0 meV is reached at 10 ps after 51.8 μJ cm −2 illumination. Between 10 ps and 125 ps, the SPV-induced shift follows an unusual linear decay profile. After 125 ps the SPV decays exponentially with a time constant of~53 ps. This apparent discrete transition from a linear decay mechanism to an exponential decay is associated with the transit of QW2 above the Fermi level, τ cross . In the absence of surface QWs, Bi 2 Se 3 exhibits a simple exponential decay with no linear region or kink 27,38 . Panel c shows the DOS at the Fermi level, E = (−0.005, 0.005) eV, and at zero momentum, k = (−0.005, 0.005) Å −1 , demonstrating the crossing of the QW2 band bottom and indicating the Lifshitz transition as QW2 relaxes.
The continuous photo-induced modification of the QWs is observed in Fig. 2d which displays four spectra of the Bi 2 Se 3 QW spectrum at equilibrium and at three separate delay values and then schematically simplified in Fig. 2e. The spectra have been shifted to compensate for the extrinsic delay-dependent rigid energy shift associated with the SPV effect (black-lined white dots in panel b) [35][36][37][38] . At equilibrium (far left), only the topological surface state and QW1 are present and occupied below the Fermi level, E F . The full band-bending is present creating a steep QW potential as we described in scenario I of panel a. This results in a large energy spacing for QW1 and QW2, leaving QW2 unpopulated above E F .
For the pumped spectrum, we first look at the data at 10 ps. By this time, nearly all hot carrier recombination has already taken place [45][46][47] , leaving behind long-lived trapped charges (see Supplementary Fig. 4). QW1 shifts downward in energy and QW2 falls below E F as photoexcited charges have reduced the equilibrium band-bending modifying the QW potential; QW2 is subsequently filled below E F and states above E F are populated resulting in a new quasi-equilibrium Fermi level. At τ = 80 ps, the QW spectrum drifts towards higher energies relative to the topological surface state and E F . The competition between the rate of charge being pushed above E F by the upward energy shifts (dotted horizontal lines in Fig. 2e) and the decay rate of charge above E F allows charge to existing above the Fermi level during the course of the spectrum relaxation. This ultimately leads to excess charge remaining in QW2 at τ = 370 ps even when the QW energy spacing has returned to equilibrium. For each delay, we can determine the photo-induced change in the surface potential, ΔV SPV (τ), from ΔKE in panel b, thereby tracking the QW modifications as a function of surface voltage.

Reduction of QW spacing and prolonged charge storage
We further investigate the delay dependence of the QW level spacing and charge occupation in order to better understand the pump-induced QW state modification. In Fig. 2, we demonstrated that the surface potential is a function of delay, and so we are able to track the QW modification with respect to surface voltage. We can take this one step further and estimate the effective field, E eff , from the energy positions of the QW1 and QW2 band bottoms. Figure 3a displays momentum integrated ∂ 2 /∂E 2 spectra for QW1 and QW2 band bottoms as a function of delay. The momentum integration windows for QW1 and QW2 band bottoms are centered at k = 0.000 Å −1 and k = 0.023 Å −1 , respectively, with a width of Δk = 0.002 Å −1 . The QW1 and QW2 band bottom positions, E 1 and E 2 , (black-lined white dots) track the evolution of the energy levels as the QW spectrum is modified . The effective field on the quantum well system is extrapolated from the average best fit line and the triangular-well approximation. c Integrated spectral intensity of the QW2 band bottom using integration window in a over a large fluence range. E F crossing time, (τ cross ), +10 ps is used as a cutoff (minimum set at 10 ps) for exponential fits (red). d QW2 decay rate, τ decay , (black open diamonds) and τ cross (red dots) as a function of maximum SPV at each fluence (and extrapolated effective field on QWs) indicating that the lifetime of QW2 states above E F is fixed but quasiequilibrium lifetime is set by peak field strength. d y axis error bars are from fit covariance in c (τ decay ) and linear fit in Supplementary Fig. 9 (τ cross ), fluence error is 10% conservatively.
S.T. Ciocys et al. by the field effect and subsequent relaxation. The positions of E 1 and E 2 are displayed in Fig. 2d for comparison with the full dispersion. After pump illumination, QW1 and QW2 shift toward higher binding energies, dropping QW1 by~25 meV and QW2 bỹ 35 meV and bringing QW2~25 meV below E F . Both QWs shift back towards lower binding energies linearly as a function of delay at different rates. For QW2, 126 ps corresponds to τ cross (see vertical dashed line in panel a), the delay time in which the band bottom crosses E F . After this point, the spectral intensity of the band bottom decreases at an accelerated rate.
An important property of QW systems is the intersubband-gap, defined as the difference in energy between the QW2 and QW1 band bottom energies ΔE = E 2 − E 1 , which sets crucial material properties including the optical response. The delay dependence of the QW band bottom positions, and therefore the surface potential dependence of the intersubband-gap, demonstrates ultrafast-field tunability. Panel b shows the potential dependence of the intersubband-gap ΔE for the two highest fluences. We focus on the band-gap values for τ < τ cross where QW2 is fully populated and the position of the band bottom is not strongly influenced by the Fermi level.
The effective applied field on the QW spectrum can be extrapolated from the intersubband gap. The triangular-well potential is a widely used approximation to describe surface confinement, especially for low-energy wavefunctions, and is suited for capturing the effective electric field experienced by the QWs (see Supplementary Fig. 7d). A more precise model of the confinement potential may be necessary for the condition of higher-energy wavefunctions (n > 2) 12,48-50 .
The triangular-well potential is described by, V(z) = eFz, in which F represents the steepness of the well with units of the electric field. As such, the QW wavefunctions take on the form of Airy functions with energy eigenvalues 48,50 : The energy eigenvalues, ϵ n , are referenced to the bulk conduction band minimum before quantization, which represents the bottom of the QW potential. By subtracting ϵ 1 from ϵ 2 , the equation can be rearranged to solve for F using only the intersubband-gap energy, and the known Airy coefficients, c 2 ≈ 2.338 and c 1 ≈ 4.088.
Utilizing the above equation and the average best fit lines from panel b, the x axis can be rescaled in terms of the effect applied field on the QW states. We use the average effective mass from dispersion fits of m * = 0.45 m e . The total electric field in the QW region at equilibrium (ΔV SPV = 0) is 42.0 mV nm −1 and 55.7 mV nm −1 at the maximal SPV (ΔV SPV = 78.5), reasonable electric field values for surface QW systems 23 . This is equivalent to applying a field, E eff between 0 mV nm −1 (ΔV SPV = 0 mV) and 13.7 mV nm −1 (ΔV SPV = 78.5 mV) opposite to the equilibrium band-bending field E BB .
The field effect on the intersubband-gap shown in Fig. 3b exhibits a nearly linear dependence since the QW electric field, F, is always far from 0, existing in the range of (40-60 mV nm −1 ). The applied field axis ticks are determined by the y axis values of the best fit line and are independent of the surface voltage axis. The even spacing of both x axis ticks demonstrates that the intersubband gap is near-linear in both surface voltage and effective electric field, as expected from the center of charge separation model in Equation (2).
In addition to enabling an ultrafast-field effect on the QW energy positions and intersubband-gap, the ultrafast gating grants access to the charge dynamics of the QW spectrum. In particular, we uncover quasi-equilibrium charge occupation and determine recombination timescales, not accessible by any static gating technique. Furthermore, we observe that the crossing of the QW2 band bottom above E F (vertical dashed line in Fig. 3a, τ cross = 126 ps at 59.1 μJ cm −2 ) corresponds to the transition of the SPV relaxation from a linear to an exponential decay mechanism as observed in Fig. 2b. Figure 3c displays the integrated spectral intensity of the QW2 band bottom as a function of delay (shaded region in panel a) for fluences between 1.5 and 59.1 μJ cm −2 , showing stable charge storage before fluence-dependent τ cross values that then depopulate exponentially. At the lowest fluences, there is no stable region and only exponential behavior. The decay region in red is fitted with exponential fits, with the start of the fitting range (red dots) set by τ cross + 10 ps and a minimum cutoff of 10 ps.
As a function of fluence, the resulting timescales from the exponential fit, τ decay , are plotted in Fig. 3d (diamonds) alongside τ cross (red dots). The decay rate of charge in QW2, τ decay , is independent of fluence with an average decay time of 53 ps. On the other hand, τ cross increases logarithmically for fluences above 4 μJ cm −2 , marking a critical fluence for the Lifshitz transition of QW2. For fluences in which the τ decay = 0, QW2 never dips below E F after the pump, and therefore does not undergo a Lifshitz transition. All together, this demonstrates that the charge stored in QW2 is stable while it is below E F and charge above E F decays with a timescale of~53 ps. Moreover, the timescale for stable stored charge in QW2, τ cross , is a consequence of the dynamic Lifshitz transition and can be controlled by varying the pump fluence, equivalent to increasing the SPV strength.

Ultrafast reversible reduction of Rashba splitting in QW1
Another appealing capability of time-resolved photo-gating is the ultrafast control and tunability of the Rashba splitting since the Rashba effect depends on the inversion symmetry breaking electric field at the surface and can therefore be tuned by an applied field [51][52][53] . The Rashba effect is an energy splitting that is, in the simplest form, linear in momentum with ΔE R = α R k where α R k is known as the Rashba parameter. In the case of parabolic bands, this leads to a rigid splitting in momentum k 0 = α R m/ h 229,54 . By directly fitting the bands, we can obtain an accurate determination of α R . In Fig. 4a, we present the band fitting of the Rashba-split QW1 band at equilibrium and τ = 30 ps after the pump. The momentum positions of the Rashba bands as a function of energy (red and blue dots) were determined by peak fitting the momentum distribution curves along delay-dependent binding energy ranges and avoiding the Rashba crossing at Γ.
In panel b, the delay dependence of the Rashba splitting is clearly seen in the QW1 band bottom momentum positions captured by integrating the second derivative in energy spectra from panel a over a narrow energy range 15 meV below the dispersion, avoiding the increased intensity at the Rashba intersection at Γ. The black dashed curve serves as a guide to the eye, indicating a 15% reduction in k 0 at τ = 0 that returns to equilibrium. The peak positions for the negative-momentum band bottom are plotted in pink and closely follow the 15% reduction. Figure 4c presents the resulting change in Rashba parameter (yellow dots) as a function of delay at high pump fluence demonstrates a linear relationship of the Rashba effect with respect to the surface photovoltage and therefore the applied effective field, and subsequently dependent on the dynamic Lifshitz transition. Figure 4d plots the percent change of the Rashba parameter as function of the SPV at each delay as well as the extrapolated effective field on the QW. The curves are spaced artificially by 5% intervals for clarity. The field effect behavior of the Rashba parameter is linear for all fluences, decreasing as the applied field is driven to more negative values with~15% reduction in Rashba splitting at a maximal applied field of 13.5 mV nm −1 . From the average of the fluence measurements, the response of the Rashba parameter versus field strength is 1.26 ± 0.21% nm mV −1 . This ultimately demonstrates precise control of the Rashba splitting of the QW states with ultrafast photo-gating, with the Rashba parameter tuned by the pump pulse intensity and dependent on the Lifshitz mediated relaxation.
Lifshitz-dependent surface charge trapping Finally, in Fig. 5 we report a non-monotonic dependence of the quasi-Fermi level, E Ã F , with respect to fluence that demonstrates the significance of the Lifshitz transition on the infrared response of the QW system.
The delay dependence of E Ã F for all measured fluences is plotted in Fig. 5a. E Ã F is measured in reference to the equilibrium Fermi level (i.e. the difference between the red fit and white fit from Fig. 2b). For fluences above 4 μJ cm −2 , E Ã F is offset from equilibrium starting in a plateau region for delays less than τ cross preceding a small peak after τ cross that leads into exponential decay. Curiously, the offset at early delays is non-monotonic with fluence, peaking at~15 meV for 5.7 μJ cm −2 and leveling off towards~5 meV at high fluence. The encircled points indicate the values that are averaged to generate the E Ã F offset values in panel e. The hump that occurs near τ cross , may be indicative of a large charge transfer from QW2 to the TSS following the reverse Lifshitz transition of the QW2 singularity. This may additionally suggest that the sustained E Ã F plateau region is also from the charge in the QW states leaking into the TSS during relaxation. For the purposes of this study, we focus on E Ã F at early delays. The relatively small offset of E Ã F is in contrast to previous time-resolved studies on p-type Bi 2 Se 3 -related samples without QWs 24 , emphasizing the influence of QW1 and QW2 on the relaxation dynamics.
The behavior of E Ã F at early delays with respect to pump fluence can be explained within the framework of the stored equilibrium charge in QW1 and QW2 described in Fig. 3. To start, in Fig. 5b we display the fluence dependence of the maximal SPV, |ΔV SPV (τ = 10 ps)|, which is proportional to the initial net charge at the surface after hot carrier relaxation, Δn surf = n sep . Next, the fraction of charge stored below E F in QW1 and QW2 is calculated. To do this, we use linear extrapolations of the QW1 and QW2 band bottom positions as a function of SPV from the high fluence data in Fig. 3a and the cumulative DOS of Rashba-split parabolic bands shown as dashed blue curves (see Supplementary Fig. 8f). The band bottom positions, and therefore the cumulative DOS, are functions of the SPV strength and are calculated from the phenomenological fit of the SPV maximal values (dashed red curve). The sum of the QW1 and QW2 cumulative DOS corresponds to the total available quasi-equilibrium charge states for photoexcited charge (solid blue curve). Figure 5c depicts two photoexcitation regimes that summarize the role of the Lifshitz transition on E Ã F . The brackets on the left indicate the pump-induced QW energy shift from equilibrium positions, the blue and red regions correspond to excited charge below and above E F , respectively. At moderate fluences of 3−7 μJ cm −2 QW2 undergoes a Lifshitz transition, however, only QW1 has significant available states below E F leaving excess separated charge to populate above E F . At high fluence, QW2 has shifted well below E F enabling a large degree of quasi-equilibrium charge storage reducing the burden on states above E F .
To further quantify the effect, it is more useful to reparameterize the curves of panel b in terms of the surface voltage to avoid reliance on the SPV phenomenological fit. In this way, the total net charge concentration is a linear function and the QW1 and QW2 band bottoms are linearly dependent on the x axis. Because the QW1 and QW2 band bottom positions are empirically determined, the only consequential free parameter is the relative amplitude of the total available quasi-equilibrium charge states (blue regions) with respect to n surf (dashed black line). The difference between the aforementioned populations represents the excess photoexcited charge at the surface (light-red region), which describes the charge remaining after QW1 and QW2 are filled below E F and that must occupy above-E F states. From panel d, it is clear that the emergence of the QW2 band bottom below E F , which acts to enhance the DOS contribution of QW2 below E F , necessitates a peaked behavior in the excess charge population near the Lifshitz transition, regardless of the fluence dependence of the SPV effect.
The proportionality factor between the charge stored in QW1 and QW2 with respect to n surf is chosen such that the excess charge best fits the E Ã F offset plotted in panel e. The excess charge from the model in panel d closely captures the non-monotonic behavior in the measured E Ã F values.

DISCUSSION
Using the reported dielectric constant of Bi 2 Se 3 of a nominal ϵ r ≈ 30−100 55,56 , we can extract an approximate SCR width, describing the spatial extent of the band-bending utilized in Fig.  1b, d ¼ ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi 2 ΔV BB;0 ϵ r ϵ 0 = e n A p % 10 nm 34 , where the acceptor ion concentration, n A , is approximately equal to the bulk majority carrier density in the SCR and V BB,0~3 00 mV based on the observation of the QW3 energy level at short delay and the zeropoint energy of the QW spectrum (see Supplementary Fig. 6). A second independent method is to use charge conservation at the SPV saturation, since band-bending cancellation occurs when the total surface charge is equal to the volumetric bulk charge times the depth. The maximal surface charge can be calculated by Luttinger's theorem 57 on the QW1, QW2, and TSS bands giving 5 × 10 12 cm −2 (see Supplementary Fig. 8), and since only~1/4 (from ΔV SPV /V BB,0 ) of the charge is screened in the SCR we get, d = 5 × 10 12 cm −2 /(1 × 10 19 cm −3 /4)~10 nm. These estimates ignore possible variations in both the dopant concentration and dielectric constant between surface and bulk, as well as hexagonal deformation and possible non-isotropic charge storage at the Fermi surface 20,58 . However, the approximate SCR width is consistent with the band-bending length-scale observed previously in Bi 2 Se 3 59 . Initial reports of the Rashba effect on metallic surfaces suggested that applied fields on the order of 1000 mV nm −1 are insufficient to produce the observed Rashba effect, and that the nuclear core potentials are responsible for large splittings 52,60 . In our experiment we observe a 15% reduction on an~1 eV Å Rashba effect with fields near 10 mV nm −1 only. There is a large body of evidence demonstrating small changes to surface or interface fields generate large changes to the Rashba parameter 51,61,62 . A reasonable explanation that has been presented is that the applied field plays an indirect role on the Rashba effect by shifting electrons to and from the nuclear cores at the surface plane 53,61,[63][64][65] . With respect to Bi 2 Se 3 , Eremeev et al. presents a debated hypothesis on the origin of the QWs from surface modulated van der Waal's gaps 66 . Although the origin of the QWs has been further argued to be from surface band-bending 20,30 , their study indirectly demonstrates the role of wavefunction Outlined points indicate averaged points used in e. b Maximal SPV as a function of fluence indicating separated charge concentration at the surface. Labeled blue curves demonstrate QW1 and QW2 states below E F determined from the phenomenological fit (dashed red) and the SPV dependence of the band bottom positions in Fig. 3a. c Graphical comparison of quasi-Fermi level shift between medium and high fluence. Medium fluence (3−7 μJ cm −2 ), QW2 photo-carrier generation is significantly larger than available states in shifted QWs. At high fluence, photo-carrier generation has saturated and QW2 provides increased storage below E F . d Representation of b with the x axis parameterized by SPV surface potential. Excess carriers after filling QW trap states (light-red) contribute to raising E Ã F (black-lined white dots). eE Ã F at short delays from a and total surface charge minus trapped surface charge (light-red curve)(black dashed curve minus solid blue curve in d) as a function of surface voltage. Error bars for Δn and ΔV SPV are conservative 10% and 2.5 meV, respectively. Y axis error in e determined by difference between E Ã F measured at QW1 and TSS k F .

S.T. Ciocys et al.
overlap with the surface-most quintuple layer as a key ingredient to the Rashba effect. We have shown through our simplified bandbending calculation of the QW wavefunctions in Supplementary  Fig. 7e that the shift in QW position is field-linear and nearly 2 Å over 10 mV nm −1 , a significant fraction of the quintuple layer thickness. In totality, our results are consistent with a picture that the origin of the Rashba splitting is largely controlled by wavefunction overlap with asymmetric nuclear potentials within the first quintuple layer, controlled linearly by the field. One consequence of the ultrafast QW modification is a large rapid and reversible change in the DOS at the Fermi level. This is made possible by the emergence of QW2 below E F after pump excitation; changing the topology of the Fermi surface and triggering a reversible Lifshitz transition. The kink in the relaxation of the SPV as well as the Rashba modification is correlated with τ cross , indicating that the reversal of the Lifshitz transition plays a key role in the relaxation. Furthermore, the new states near the Fermi level imply an increase in conductivity at the surface affecting spin-dependent transport properties, and occupation of the QW2 band bottom would influence intersubband optical absorption 67 . As such, this ultrafast-field effect device provides many opportunities for applications that are sensitive to the DOS at the chemical potential and requires further attention exploring the DOS-dependent scattering dynamics.
The SPV photo-gating method can be extended to other semiconductor systems with 2D electron gases. In fact, a similar gating method to the one presented here has been utilized successfully in black phosphorous to drive gap renormalization at the surface 68,69 . An example of a future application is a heterojunction system in which an exfoliated or MBE-grown 2D material is placed on top of an SPV substrate separated by a thin insulating layer. In this case, a pump pulse focused on the grounded top layer would result in transient gating of the illuminated surface region. Moreover, SPV time constants can vary significantly from picoseconds to milliseconds depending on the system allowing for a large dynamic range of study 36,69,70 .
In conclusion, we have exposed complex QW dynamics using a time-resolved SPV-based gate. The phenomenology of the photoinduced QW modification is general for surface QW systems because it relies on the surface band-bending of not only topological insulators but many semiconducting systems 34 . For this system, the unique decay profile of both the SPV and Rashba reduction is contingent on the dynamic DOS near the Fermi level and the existence of a van Hove singularity in the Rashba-split QWs. Accordingly, we have demonstrated local control of the near-E F DOS; a useful tool for driving electronic phases such as superconductivity 71 , magnetism 72 , and, more broadly, correlated phases 73 . Furthermore, the ultrafast gating mechanism also revealed dynamically stored charge in QW1 and QW2 that ultimately controls the resulting quasi-Fermi level, impossible to observe using static gating techniques. Lastly, we have demonstrated ultrafast reversible modulation of the Rashba effect on 10-100 picosecond timescales opening up the possibility for GHz spin-current manipulation. The temporary photo-induced Rashba tuning creates a platform for spin-orbit modulation enabling lightdriven spintronic devices such as ultrafast spin-transistors and photo-controlled Rashba circuitry 14,74,75 .

Sample synthesis and preparation
The p-type Bi 2 Se 3 single crystals were grown via directional slow solidification with ≈1 at.% Mg substitution with bismuth in order to bring E F close to the bulk valence band maximum. The crystals were coated in graphite spray and then cleaved in a vacuum along the (111) plane at a base pressure of 5 × 10 −11 Torr. The sample experienced rapid dosing triggered by initial exposure to the infrared pump resulting in the formation of a clear QW spectrum similar to the observations by Lin et al. 8 , with no further dosing occurring after subsequent pumping. This dosing is perhaps due to the presence of trapped gasses in the graphite spray or the crystal itself, and produced results similar to the effects of carbon monoxide and water surface dosing 12,13 . We have observed this dosing effect multiple times but only after using graphite spray on specifically p-type Bi 2 Se 3 samples.

Measurement details
For the pump-probe measurement, we used a pump of 1.48 eV p-polarized photons generated by a cavity-dumped Ti:sapphire oscillator operating with repetition rates of 1357.50 kHz (51.8 μJ cm −2 ) and 678.75 kHz (all other fluences). The probe beam consists of 5.94 eV p-polarized photons generated by frequency quadrupling of the 1.48 eV source in two BBO stages 76 . The total time resolution from the two beams is~300 fs. The pump and probe beam spot full-width-half-maxima was 120 μm and 30 μm, respectively. The delay between the pump and the probe is achieved with a mechanical translation stage. Due to the long delay range, careful consideration was taken to ensure that the maximum pointing drift due to the delay stage was minimized to <15 μm. A variable neutral density filter applied to the pump beam enabled fluence-dependent measurements. All reported fluence values indicate incident fluence on the Bi 2 Se 3 surface. All spectra were taken at the system's base temperature with liquid nitrogen of~80 K. Data were analyzed using PyARPES analysis framework 77 .

DATA AVAILABILITY
All data and code used for analysis are made available upon request to the corresponding author.