THz-STM & Ultrafast Pump–Probe STM

• THz-STM and ultrafast pump–probe STM are advanced scanning probe techniques that combine intense, phase-stable THz pulses with STM for time-resolved atomic-scale measurements.
• These methods employ near-field enhancement and pump–probe schemes to map local charge, vibrational, and electronic dynamics with femtosecond to picosecond resolution.
• Applications span probing quantum, electronic, and topological dynamics in low-dimensional materials, enabling direct observation and manipulation at the atomic level.

Terahertz Scanning Tunneling Microscopy (THz-STM) and ultrafast pump–probe STM represent a class of lightwave-driven scanning probe techniques that couple ultrafast THz transients to an STM junction, enabling real-space, time-resolved measurements at the atomic scale. Leveraging the extreme spatial confinement inherent to the STM tip–sample gap and the high field enhancement of antenna effects, these methods provide access to local charge, vibrational, and electronic dynamics on timescales down to a few femtoseconds. Combined with generalized pump–probe schemes where THz and/or optical pulses are used as temporal gates or probes, THz-STM extends the capabilities of conventional STM to interrogate and control a range of quantum, electronic, and structural processes in low-dimensional materials, correlated systems, and molecular junctions (Jelic et al., 2023, Li et al., 2023, Sabanés et al., 2022, Allerbeck et al., 2024, Bae et al., 15 Jul 2025, Rai et al., 9 Jun 2025, Jacobsen et al., 16 Feb 2026).

1. Fundamental Principles and Physical Mechanisms

THz-STM harnesses intense, phase-stable, single-cycle THz pulses to transiently drive the STM junction far from equilibrium. The THz field induces a time-dependent voltage bias $V_{\mathrm{THz}}(t)$ across the nanometer-scale gap, which can be locally enhanced by factors of $10^2$–$10^5$ via the tip’s antenna response. The resulting tunneling current $I(t)$ is governed by the instantaneous barrier potential: $$I(t) \approx I_0 V(t) \exp\left(-\frac{2 z_0}{\hbar} \sqrt{2 m (\Phi - e V(t))} \right),$$ where $\Phi$ is the effective work function, $z_0$ the tip–sample separation, $m$ the electron mass, and $V(t)=F E_{\mathrm{inc}}(t) z_0$ (Jelic et al., 2023).

In pump–probe modalities, a strong “pump” THz transient first excites the system (e.g., charges a localized defect), and at a variable delay $\Delta t$, a second, temporally synchronized “probe” pulse interrogates the evolving state by triggering additional tunneling or current rectification. By mapping the resultant current or transferred charge as a function of delay, one reconstructs ultrafast population dynamics, dephasing rates, or oscillatory modes intrinsic to the sample (Allerbeck et al., 2024, Li et al., 2023).

2. THz Pulse Generation, Delivery, and Near-Field Enhancement

Single-cycle THz pulses for STM are typically generated via optical rectification in nonlinear crystals (e.g., LiNbO₃, ZnTe) using few-μJ, near-infrared femtosecond pulses (Yb-fiber or Ti:Sa lasers). Alternative sources include spintronic THz emitters (STE) driven by ultrafast NIR pulses, which achieve bandwidths up to 30 THz (Müller et al., 2020).

The free-space THz beam is directed and tightly focused onto the tip apex in UHV STM heads using large numerical aperture parabolic mirrors. Crossing the gap, the incident field $10^2$0 is near-field enhanced by the STM tip–sample geometry to $10^2$1, with $10^2$2–$10^2$3 depending on tip radius, gap, and frequency (Jelic et al., 2023, Li et al., 2023). This enhancement enables rectification of μA-level field emission bursts and deterministic ultrafast gating.

Control over the THz waveform’s carrier-envelope phase (CEP), polarization, and amplitude is provided by custom CEP shifters, wire-grid polarizers, and delay lines. Through a Michelson or Mach–Zehnder interferometer, two pulse replicas—pump and probe—can be generated with tunable inter-pulse delay, while pulse tailoring enables transformation of free-space bipolar pulses into unipolar near-field drives, critical for unambiguous time-resolved sampling (Jelic et al., 2023, Li et al., 2023, Allerbeck et al., 2024).

3. Time-Resolved Sampling and Pump–Probe Protocols

In THz-STM and pump–probe STM, the temporal evolution of electronic or structural states is captured by mapping the tunneling response as a function of relative delay between pump and probe pulses. The primary detection modes include:

• THz pump–THz probe: A strong-field pump pulse transiently charges a localized state or alters the local electronic configuration; a weak probe pulse, introduced at a controllable $10^2$4, samples the occupation or continuing dynamics, with the resulting rectified current $10^2$5 reflecting the real-time evolution of the initial excitation (Allerbeck et al., 2024, Jelic et al., 2023).
• Ultrafast gate sampling: In field emission STM, the instantaneous tunneling current spike ($10^2$6100–200 fs) provides a faithful replica of the near-field waveform, enabling direct, in-situ measurement of $10^2$7 at sub-nanometer lateral resolution (Li et al., 2023).
• Photoemission-based detection: A few-femtosecond NIR pulse can trigger multiphoton photoemission; simultaneous or subsequent THz gating then encodes the local THz transient onto the measured photoinduced current (Müller et al., 2020, Sabanés et al., 2022).

Spatial resolution in all these modes is essentially unchanged from standard STM, with full width at half maximum (FWHM) $10^2$80.85 nm for THz-driven field emission, and time resolution limited by the THz pulse envelope or, in photocurrent sampling, by the NIR pulse width (as low as 8–10 fs) (Li et al., 2023, Jelic et al., 2023, Sabanés et al., 2022).

4. Theoretical Models and Data Interpretation

Time-dependent tunneling under THz drive is fundamentally a nonequilibrium process. Formulations span:

• Rate equation/master equation models: For charge-state dynamics in defects, vibronic ladders are populated with transitions modeled by time-dependent rates incorporating Franck–Condon factors, tunnel coupling, and phonon relaxation (Allerbeck et al., 2024).
• Anderson–Holstein frameworks: Electron–electron and electron–phonon correlations in single-molecule or quantum-dot systems are treated using variational non-Gaussian states to capture dynamics such as Kondo resonance collapse, phonon sidebands, and long-lived vibrational coherences after THz excitation (Shi et al., 2019).
• Nonequilibrium Green’s function (NEGF) theory: For addressing time-dependent local density of states and rectified currents in the presence of Floquet driving or inhomogeneous fields, NEGF expressions compute $10^2$9, $10^5$0, and their relation to transient or steady-state LDOS (Jacobsen et al., 16 Feb 2026).

The applicable phenomenology includes adiabatic tunneling (when junction response time $10^5$1 field variation), ultrafast nonthermal tunneling, and delayed thermionic emission by hot electrons with relaxation times extracted from two-temperature models (Sabanés et al., 2022).

5. Applications to Nanoscale Charge, Lattice, and Topological Dynamics

THz-STM and ultrafast pump–probe STM have been applied to a variety of quantum material systems:

• Atomic-scale quantum dots and charge localization: Ultrafast Coulomb blockade effects in WSe₂ defects were directly visualized, with characteristic charge-state lifetimes ($10^5$2) and back-tunneling dynamics resolved on ps timescales. The Franck–Condon blockade was shown to selectively suppress back-tunneling by restricting vibronic transitions and favoring unidirectional charge flow (Allerbeck et al., 2024).
• Coherent phonon excitation in 2D semiconductors: Selective, defect-modulated driving of forbidden breathing and shear phonon modes in MoTe₂ was achieved, with mode lifetimes $10^5$330–50 ps and amplitudes tunable by local band bending at defects. Local dipoles induced by band bending modify the field–phonon coupling and allow controllable excitation on the atomic scale (Rai et al., 9 Jun 2025).
• Ultrafast observation of collective modes: Pump–probe trSTM captured the dynamic competition between phasons and amplitudons in charge density wave insulators such as (TaSe₄)₂I, with sub-ps oscillations at 0.22 THz and a parametrically-amplified 0.11 THz “daughter” phason mode observable through damped sinusoidal fits and Fourier analysis (Bae et al., 15 Jul 2025).
• Real-space Floquet engineering: Ultrafast THz-STM has been proposed and modeled as a local probe of Floquet gaps and edge states in graphene and nanoribbons under strong mid-IR circular driving, exploiting the STM's Å-scale spatial resolution to map dynamic topological states inaccessible to bulk transport or conventional photoemission (Jacobsen et al., 16 Feb 2026).

6. Technical Limitations, Energy and Temporal Resolution

Key technical constraints of THz-STM and ultrafast pump–probe STM include:

• Tip-antenna low-pass filtering: The frequency response of the tip–sample junction limits the attainable THz bandwidth; current state-of-the-art resolves up to 15 THz, with half-cycle temporal features as short as 115 fs (Müller et al., 2020).
• Unipolar pulse requirements: Accurate time-domain mapping depends on generating near-unipolar driving fields and isolating single current bursts, necessitating careful waveform engineering and self-consistent calibration (Jelic et al., 2023).
• Sample conductivity: THz rectification requires a conducting or degenerate junction; nonmetallic samples cannot be measured directly by THz-CC methods (Jelic et al., 2023).
• Ultimate time resolution: Set by the shorter of the NIR pump (in photoemission-based STM) or the duration of the THz field emission spike; values $10^5$410–20 fs are standard, but may be driven to the attosecond regime with further advances (Sabanés et al., 2022).
• Energy resolution: For differential THz-STS, the weak-probe amplitude $10^5$5 determines the energy window, with practical values yielding $10^5$6 a few meV (Jelic et al., 2023).

7. Outlook and Future Directions

Advances in THz-STM instrumentation and methodology are driving the field towards broader applicability and finer resolution. Expected developments include:

• Higher field enhancement and bandwidth: Plasmonic and antenna-engineered tips will increase $10^5$7 and thus attainable bias and THz bandwidth, allowing access to even faster charge and spin phenomena (Li et al., 2023, Müller et al., 2020).
• Simultaneous multimodal readout: Integration with local optical, vibrational, or luminescence detection (“THz–STML”) and implementation of multispectral pump–probe protocols will expand the accessible observables (Li et al., 2023, Jelic et al., 2023).
• Emergent quantum and topological dynamics: Extension to ultrafast control and probing of Floquet phenomena, nonlinear excitations, and light-induced ordering in correlated, topological, and low-dimensional systems (Jacobsen et al., 16 Feb 2026, Bae et al., 15 Jul 2025).
• Single-molecule and few-electron dynamics: By leveraging attosecond gating and carrier-envelope phase stabilization, future THz-STM platforms may resolve single-electron tunneling, vibrational wavepacket evolution, and nonadiabatic molecular switching at atomic resolution (Allerbeck et al., 2024, Sabanés et al., 2022, Shi et al., 2019).

THz-STM and ultrafast pump–probe STM thus establish a versatile and robust platform for time-resolved, atomic-scale quantum measurement, bridging the gap between femtosecond optics and STM’s intrinsic spatial precision (Jelic et al., 2023, Sabanés et al., 2022, Müller et al., 2020, Allerbeck et al., 2024, Li et al., 2023, Bae et al., 15 Jul 2025, Rai et al., 9 Jun 2025, Jacobsen et al., 16 Feb 2026).