REVIEW 3 major objections 5 minor 2 cited by
Ultrafast Coulomb blockade in an atomic-scale quantum dot
T0 review · 3 major / 5 minor · reviewed 2026-08-11 · deepseek-v4-flash
Pith's one-line read This paper reports direct, real-space, time-domain observation of transient Coulomb blockade at a single selenium vacancy in WSe2 using THz pump–THz probe scanning tunneling microscopy.
desk verdict A real experimental milestone—time-domain Coulomb blockade at a single defect—with a plausible but under-constrained Franck-Condon mechanism; referee it, but insist on sensitivity analysis and data release. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The mechanism that carries the argument is the Franck–Condon blockade, implemented in a time-dependent master equation over three electronic states (neutral ground state, LUMO, LUMO+1) dressed by a single vibrational mode with Franck–Condon factors. The paper defines the parameters — ℏΩ = 8 meV, Huang–Rhys factor S = 2.2 for 2 ML and S = 5 for 1 ML, 3 meV Gaussian broadening, and phonon relaxation time τph = 1 ps — and uses them to compute the rectified charge from the actual THz waveform including reflections. A modified plate-capacitor voltage drop is introduced to model the distance dependence.
What would settle it
Measure the THz pump–probe trace on a selenium vacancy with a substantially weaker electron–phonon coupling (small Huang–Rhys factor): the Franck–Condon blockade predicts that τeff should stay short at ΔV ≈ 0 because back tunneling remains strong, so a persistent rise in τeff would rule out the proposed mechanism.
Extended reading notes
Core claim
Using THz pump–THz probe time-domain sampling, the authors record atomic-scale snapshots of the transient Coulomb blockade of individual selenium vacancies in WSe2 and extract effective charge-state lifetimes. For a vacancy in 1 ML WSe2 the intrinsic lifetime is about 2.2 ps, while in 2 ML WSe2 it ranges up to 86 ps; the measured effective lifetime rises steeply as the dc bias approaches the LUMO from below and peaks near the intrinsic value at ΔV ≈ 0. The paper's key claim is that this behavior is governed by the Franck–Condon blockade: because vibrational relaxation is fast, transitions occur from the vibrational ground state of the charged defect, and at high bias the overlap with accessible low-lying vibrational states of the neutral defect is strongly reduced, so back tunneling to the tip is suppressed. The same master-equation model, with a single 8 meV vibrational mode and Huang–Rhys factor 2.2, reproduces the rectified-charge curves as functions of bias, tip–sample distance, and pump–probe delay.
Load-bearing premise
The interpretation rests on a single-mode vibronic model with fitted parameters — Huang–Rhys factor 2.2, phonon energy 8 meV, 3 meV broadening, and 1 ps phonon relaxation time — and the paper itself states that the one-phonon-mode treatment breaks down at bias offsets near 10 mV.
Editorial extensions
If this is right
- Lightwave-driven STM can now probe localized charge dynamics even when charge-state lifetimes exceed the THz pulse duration, provided back tunneling is suppressed by the Franck–Condon blockade.
- Measuring τeff as a function of dc bias near the defect resonance gives access to the intrinsic charge-state lifetime τ0 from a time-domain experiment.
- The Franck–Condon blockade acts as a charge backflow valve, so the same junction can be operated in weak-injection, strong-injection, and blockade-limited regimes by choosing bias and tip height.
- Tuning bias within one phonon energy of the resonance switches between back-tunneling-dominated and forward-transfer-dominated regimes, enabling unidirectional charge transport.
- The approach can track charge transfer into the substrate with picosecond temporal resolution and atomic spatial resolution.
Reading between the lines
- The sharpness of τeff(ΔV) near resonance is a sensitive measure of electron–phonon coupling, so the technique could be used as a local, time-domain probe of Huang–Rhys factors.
- Because the paper's model uses a single phonon mode and the authors note its breakdown near ΔV ≈ 10 mV, quantitative lifetime extraction may carry systematic errors; a two-mode or temperature-resolved treatment could test whether low-energy interlayer phonons dominate the blockade.
- The same back-tunneling-suppression strategy might be generalizable to other atomic defects in two-dimensional materials and to other selection rules, such as spin or orbital angular momentum, as the authors suggest.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports THz pump–THz probe scanning tunneling microscopy measurements on individual selenium vacancies in monolayer and bilayer WSe2 and claims the first time-domain, atomic-scale observation of transient Coulomb blockade at a single defect. The central observable is the rectified charge per THz transient, measured as a function of dc bias, tip–sample distance, and pump–probe delay; from the delay traces the authors extract an effective charge-state lifetime τeff. They find that τeff grows as the dc gate approaches the defect LUMO and near ΔV ≈ 0 approaches the intrinsic, substrate-limited lifetime τ0. This suppression of back tunneling to the tip is attributed to the Franck–Condon blockade and modeled with a master-equation rate model that includes one vibrational mode (Methods, Eqs. 1–6). The paper also presents orbital-resolved differential images at selected delays, which are interpreted as snapshots of the transient Coulomb blockade.
Significance. If the central claim holds, this is a notable experimental advance: it would be the first real-space, time-domain imaging of Coulomb blockade at an atomic-scale quantum dot, and it would identify Franck–Condon blockade as a practical mechanism for achieving unidirectional charge transfer in lightwave-driven STM. The experimental methodology is strong: multiple defects are studied, several coupling regimes are varied systematically, the THz waveform is calibrated in situ, and the rate-equation model is described in enough detail to be reconstructed. The paper is also commendably transparent about some limitations, such as the acknowledged breakdown of the single-phonon-mode model at ΔV ≈ 10 mV and the presence of trailing THz reflections. However, the attribution of the observed lifetime-versus-bias behavior to Franck–Condon blockade currently rests on several fitted or assumed model parameters, most importantly the unmeasured phonon relaxation time τph, so the quantitative evidence for the mechanism is not yet established independently of those assumptions.
major comments (3)
- [Methods, Eq. (6) and “Implementation of simulation”] The Franck–Condon blockade mechanism requires the charged defect to relax to its vibrational ground state before the back-tunneling window opens, and the model implements this as a decay to the ground state with lifetime τph = 1 ps. This value is an assumption extrapolated from free-carrier thermalization in few-layer MoS2 (ref. 38) to a localized 8 meV defect mode at 5 K. If τph is an order of magnitude longer, transitions from excited vibronic levels of the charged state to low-lying neutral states would remain available, weakening the blockade and reducing τeff near ΔV = 0. No sensitivity analysis is provided. Because this single parameter controls the central mechanism, the manuscript should report how the predicted τeff(ΔV) and QLW(Vdc) curves change when τph is varied over a plausible range (e.g., 0.1–10 ps), or provide an independent experimental bound on τph.
- [Methods, “Implementation of simulation”; Extended Figs. S6b–S6g] The master-equation model is not a parameter-free simulation. The Huang–Rhys factors (S = 2.2 for 2 ML, S = 5 for 1 ML), phonon energy ħΩ = 8 meV, 3 meV Gaussian broadening, tip–LUMO+1 coupling factor Λ, effective capacitor-plane distance zVD = 13 Å, and a 0.4 Å z0 offset are all fitted or chosen to match the same sample's static dI/dV and I(z) data. The claimed “quantitative reproduction” of QLW(Vdc), QLW(z), and τeff(ΔV) therefore has limited evidentiary weight for the Franck–Condon mechanism unless the model is validated out-of-sample or the parameter uncertainties are propagated into the predicted observables. The authors' own statement that the 1-mode treatment breaks down at ΔV ≈ 10 mV further narrows the predictive range over which the comparison in Fig. 3c can be taken as support for the mechanism.
- [Fig. 5 and “Time-domain detection of the ultrafast Coulomb blockade”] The effective lifetimes are extracted by exponential fits over the restricted window 1.5–15 ps, while the claimed saturation near ΔV = 0 is compared with intrinsic lifetimes τ0 = 55–86 ps that exceed the fit window by a factor of 4–50. With trailing THz reflections at up to 30% of the main pulse amplitude (Extended Fig. S2g) and the authors' own statement that fit errors become large for τeff ≳ 20 ps, the data do not robustly establish that τeff quantitatively reaches τ0. The claim requires either a longer-time measurement with the reflection problem controlled, a quantitative uncertainty analysis of the truncated-window fits, or a demonstration that the same fitting protocol applied to the simulation is insensitive to the window choice.
minor comments (5)
- [Section heading] The heading “COULOMB BLOCKADE A T A SINGE SELENIUM V ACANCY” contains several typographical errors and should be corrected.
- [Definitions of ΔV] The sign convention for ΔV is inconsistent across the paper: some passages and Fig. 5 axes use ΔV = Vdc − VLUMO, while other captions and the Methods define the offset as VLUMO − Vdc. A single, clearly stated convention should be used throughout.
- [Methods, “Implementation of simulation”] Extended Fig. S6b shows that a two-mode Franck–Condon model with an additional 20 meV mode (S2 = 0.7) is needed to reproduce the static dI/dV spectrum, but the simulations in Figs. 3c, 4c, and 5 use only the 8 meV mode. The main text should explicitly state that the 20 meV mode is omitted from the dynamical simulations and discuss how its inclusion would affect the predicted back-tunneling asymmetry.
- [Fig. 4b caption] The color scale for the ΔQLW images in Fig. 4b is not defined in the caption; please state the reference used for the differential quantity.
- [Data availability] The data and code availability statements say materials are available “upon reasonable request.” Given the number of fitted parameters in the model, depositing the master-equation code and representative raw pump–probe traces would substantially strengthen reproducibility.
Circularity Check
No significant circularity: the central pump-probe lifetime claim rests on an independent observable; calibrated static parameters and minor self-citations do not force the result.
full rationale
Walking the derivation chain, I find no step in which a claimed prediction reduces by construction to its inputs. The model parameters (S=2.2, hbar-Omega=8 meV, 3 meV Gaussian broadening, tau_0, kappa, tip-LUMO+1 coupling) are calibrated to static STS and I(z) approach curves, and the QLW(z) simulation uses explicitly adjusted voltage-drop parameters (zVD=13 A and a 0.4 A z0 shift). However, these calibrations are disclosed in the Methods as fits rather than disguised predictions, and the central new observable, the THz pump-THz probe decay with its tau_eff(Delta V) trend, is not one of the fitted inputs. The simulations are also cross-defect: S=2.2 is fit on Vac_2.2 and used to model Vac_2.4 in Fig. 3c, and the 1 ML Huang-Rhys factor is rescaled rather than fit to the pump-probe data. tau_eff is compared with an independently measured tau_0 obtained from dc current saturation approach curves, so the near-Delta-V=0 agreement is not a tautology. The self-citations (Refs. 16, 18, 32, 44) supply prior sample characterization, calibration methods, and spin-multiplicity input, but no uniqueness theorem or forbidden-alternative argument rests on them, and they are not used to define the target claim. The assumed phonon relaxation time tau_ph=1 ps is an unmeasured parameter justified by external ultrafast carrier-cooling literature; this is a robustness and sensitivity concern, not circularity. Thus the core claim of observing transient Coulomb blockade and of Franck-Condon suppression of back tunneling has independent, non-circular empirical content.
Assumptions & free parameters
free parameters (8)
- Huang-Rhys factor S for 2 ML WSe2 =
2.2
- Huang-Rhys factor S for 1 ML WSe2 =
5
- Phonon mode energy hbar Omega =
8 meV
- Phonon relaxation time tau_ph =
1 ps
- Gaussian line broadening =
3 meV
- Effective capacitor plane distance zVD =
13 A
- Tip-LUMO+1 coupling factor Lambda =
0.5 to 1.0
- z0 uncertainty for 1 ML simulation =
0.4 A
assumptions (6)
- domain assumption Tunneling dynamics are described by a time-dependent Markov process (master equation), with rates given by first-order perturbation theory
- standard math The Franck-Condon approximation with a single vibrational mode, rigid shift, and no Duschinsky rotation
- domain assumption Fast vibrational relaxation to the ground state of the charged defect (tau_ph ~ 1 ps, few-100 fs to 1 ps)
- domain assumption The THz transient adds linearly to the dc bias at the junction
- domain assumption Voltage division in the tunnel junction follows a plate capacitor model with epsilon_r=6.5 and d=6.5 A for WSe2
- domain assumption The defect is described by three electronic states (neutral ground, LUMO, LUMO+1) with defined multiplicities
Cite this review
Pith. "Pith review of Ultrafast Coulomb blockade in an atomic-scale quantum dot." pith.science (2026). https://pith.science/paper/JWH7QGCN
@misc{pith2026241213718,
author = {Pith},
title = {Pith review of: Ultrafast Coulomb blockade in an atomic-scale quantum dot},
year = {2026},
howpublished = {\url{https://pith.science/paper/JWH7QGCN}},
note = {Machine review of arXiv:2412.13718}
}
abstract
Controlling electron dynamics at optical clock rates is a fundamental challenge in lightwave-driven nanoelectronics. Here, we demonstrate ultrafast charge-state manipulation of individual selenium vacancies in monolayer and bilayer tungsten diselenide (WSe$_2$) using picosecond terahertz (THz) source pulses, focused onto the picocavity of a scanning tunneling microscope (STM). Using THz pump--THz probe time-domain sampling of the defect charge population, we capture atomic-scale snapshots of the transient Coulomb blockade, a signature of charge transport via quantized defect states. We identify back tunneling of localized charges to the tip electrode as a key challenge for lightwave-driven STM when probing electronic states with charge-state lifetimes exceeding the pulse duration. However, we show that back tunneling can be mitigated by the Franck-Condon blockade, which limits accessible vibronic transitions and promotes unidirectional charge transport. Our rate equation model accurately reproduces the time-dependent tunneling process across the different coupling regimes. This work builds on recent progress in imaging coherent lattice and quasiparticle dynamics with lightwave-driven STM and opens new avenues for exploring ultrafast charge dynamics in low-dimensional materials, advancing the development of lightwave-driven nanoscale electronics.
Figures
Forward citations
Cited by 2 Pith papers
-
Ultrafast dynamics of local charge order in a THz-induced metastable quantum state
THz-STM reveals a long-lived metastable state in 1T-TaS2 with local stacking changes and coherent interlayer shear and amplitude mode oscillations at 2.5 and 1.3 THz.
-
Attosecond charge transfer in atomic-resolution scanning tunnelling microscopy
Single-cycle near-infrared pulses drive sub-femtosecond electron transfer across an STM junction, enabling atomic-resolution imaging with attosecond-scale sensitivity.
Reference graph
Works this paper leans on
-
[1]
Weber, J. R. et al. Quantum computing with defects. Proceedings of the National Academy of Sciences 107, 8513–8518 (2010)
work page 2010
-
[2]
Awschalom, D. D., Bassett, L. C., Dzurak, A. S., Hu, E. L. & Petta, J. R. Quantum Spin- tronics: Engineering and Manipulating Atom-Like Spins in Semiconductors. Science 339, 1231364 (2013)
work page 2013
-
[3]
Wolfowicz, G. et al. Quantum guidelines for solid-state spin defects. Nat. Rev. Mater.1–20 (2021)
work page 2021
-
[4]
Maurer, P. C. et al. Room-temperature quantum bit memory exceeding one second. Science 336, 1283–1286 (2012)
work page 2012
-
[5]
et al.Room-temperature entanglement between single defect spins in diamond
Dolde, F. et al.Room-temperature entanglement between single defect spins in diamond. Nat. Phys. 9, 139 (2013)
work page 2013
-
[6]
Heide, C., Keathley, P. D. & Kling, M. F. Petahertz electronics. Nature Reviews Physics6, 648–662 (2024)
work page 2024
-
[7]
Schofield, S. R. et al. Atomically Precise Placement of Single Dopants in Si. Physical Review Letters 91, 136104 (2003)
work page 2003
-
[8]
Fuechsle, M. et al. A single-atom transistor. Nature Nanotechnology7, 242–246 (2012)
work page 2012
Show all 49 references
-
[9]
& Hersam, M
Liu, X. & Hersam, M. C. 2D materials for quantum information science. Nature Reviews Materials 4, 669–684 (2019)
2019
-
[10]
Klein, J. et al. Site-selectively generated photon emitters in monolayer MoS2 via local helium ion irradiation. Nat. Commun. 10, 2755 (2019)
2019
-
[11]
et al.Spin-defect qubits in two-dimensional transition metal dichalcogenides operating at telecom wavelengths
Lee, Y. et al.Spin-defect qubits in two-dimensional transition metal dichalcogenides operating at telecom wavelengths. Nature Communications 13, 7501 (2022)
2022
-
[12]
T., Bray, K., Ford, M
Tran, T. T., Bray, K., Ford, M. J., Toth, M. & Aharonovich, I. Quantum emission from hexagonal boron nitride monolayers. Nat. Nanotechnol.11, 37 (2016)
2016
-
[13]
Mendelson, N. et al. Identifying carbon as the source of visible single-photon emission from hexagonal boron nitride. Nature Materials 20, 321–328 (2021)
2021
-
[14]
Senkpiel, J. et al. Dynamical Coulomb Blockade as a Local Probe for Quantum Transport. Physical Review Letters124, 156803 (2020)
2020
-
[15]
& Gross, L
Kaiser, K., Lieske, L.-A., Repp, J. & Gross, L. Charge-state lifetimes of single molecules on few monolayers of NaCl. Nature Communications 14, 4988 (2023)
2023
-
[16]
Bobzien, L. et al. Layer-Dependent Charge State Lifetime of Single Se Vacancies in WSe2. arXiv arXiv:2407.04508 (2024)
2024 arXiv
-
[17]
Taucer, M. et al. Single-Electron Dynamics of an Atomic Silicon Quantum Dot on the H-Si (100)-(2×1) Surface. Physical Review Letters112, 256801 (2014)
2014
-
[18]
Schuler, B. et al. Large Spin-Orbit Splitting of Deep In-Gap Defect States of Engineered 31 Sulfur Vacancies in Monolayer WS2. Phys. Rev. Lett.123, 076801 (2019)
2019
-
[19]
Cochrane, K. A. et al. Spin-dependent vibronic response of a carbon radical ion in two- dimensional WS2. Nature Communications 12, 7287 (2021)
2021
-
[20]
Stolz, S. et al. Layer-dependent Schottky contact at van der Waals interfaces: V-doped WSe2 on graphene. npj 2D Materials and Applications6, 1–5 (2022)
2022
-
[21]
L., Jelic, V., Hillenbrand, R
Cocker, T. L., Jelic, V., Hillenbrand, R. & Hegmann, F. A. Nanoscale terahertz scanning probe microscopy. Nature Photonics 15, 558–569 (2021)
2021
-
[22]
& Damascelli, A
Boschini, F., Zonno, M. & Damascelli, A. Time-resolved ARPES studies of quantum materials. Reviews of Modern Physics96, 015003 (2024)
2024
-
[23]
Cocker, T. L. et al. An ultrafast terahertz scanning tunnelling microscope. Nature Photonics 7, 620–625 (2013)
2013
-
[24]
Jelic, V. et al. Ultrafast terahertz control of extreme tunnel currents through single atoms on a silicon surface. Nature Physics 13, 591–598 (2017)
2017
-
[25]
Liu, S. et al. Nanoscale coherent phonon spectroscopy. Science Advances8, eabq5682 (2022)
2022
-
[26]
Jelic, V. et al. Terahertz control of surface topology probed with subatomic resolution. arXiv arXiv:2411.07545 (2024)
2024
-
[28]
Roelcke, C. et al. Ultrafast atomic-scale scanning tunnelling spectroscopy of a single vacancy in a monolayer crystal. Nature Photonics 18, 595–602 (2024)
2024
-
[29]
Jelic, V. et al. Atomic-scale terahertz time-domain spectroscopy. Nature Photonics 18, 898– 904 (2024)
2024
-
[30]
L., Peller, D., Yu, P., Repp, J
Cocker, T. L., Peller, D., Yu, P., Repp, J. & Huber, R. Tracking the ultrafast motion of a single molecule by femtosecond orbital imaging. Nature 539, 263–267 (2016)
2016
-
[31]
Imaging surfaces at the space–time limit: New perspectives of time-resolved scanning tunneling microscopy for ultrafast surface science
M¨ uller, M. Imaging surfaces at the space–time limit: New perspectives of time-resolved scanning tunneling microscopy for ultrafast surface science. Progress in Surface Science99, 100727 (2023)
2023
-
[32]
Bobzien, L. et al. Ultrafast state-selective tunneling in two-dimensional semiconductors with a phase- and amplitude-controlled THz-scanning tunneling microscope. APL Materials 12, 051110 (2024)
2024
-
[33]
& Meyer, G
Steurer, W., Gross, L. & Meyer, G. Local thickness determination of thin insulator films via localized states. Applied Physics Letters104, 231606 (2014)
2014
-
[34]
Schuler, B. et al. How Substitutional Point Defects in Two-Dimensional WS 2 Induce Charge Localization, Spin–Orbit Splitting, and Strain. ACS Nano 13, 10520–10534 (2019)
2019
-
[35]
Aghajanian, M. et al. Resonant and bound states of charged defects in two-dimensional semiconductors. Physical Review B101, 081201 (2020)
2020
-
[36]
V., Wu, S
Nazin, G. V., Wu, S. W. & Ho, W. Tunneling rates in electron transport through double- barrier molecular junctions in a scanning tunneling microscope. Proceedings of the National 32 Academy of Sciences102, 8832–8837 (2005)
2005
-
[37]
& Andreev, A
Koch, J., Von Oppen, F. & Andreev, A. V. Theory of the Franck-Condon blockade regime. Physical Review B74, 205438 (2006)
2006
-
[38]
Nie, Z. et al. Ultrafast Carrier Thermalization and Cooling Dynamics in Few-Layer MoS2. ACS Nano 8, 10931–10940 (2014)
2014
-
[39]
Leturcq, R. et al. Franck–Condon blockade in suspended carbon nanotube quantum dots. Nature Physics 5, 327–331 (2009)
2009
-
[40]
Doleˇ zal, J., Sagwal, A., De Campos Ferreira, R. C. &ˇSvec, M. Single-Molecule Time-Resolved Spectroscopy in a Tunable STM Nanocavity. Nano Letters 24, 1629–1634 (2024)
2024
-
[41]
Imai-Imada, M. et al. Orbital-resolved visualization of single-molecule photocurrent channels. Nature 603, 829–834 (2022)
2022
-
[42]
Sheng, S. et al. Terahertz spectroscopy of collective charge density wave dynamics at the atomic scale. Nature Physics 20, 1603–1608 (2024)
2024
-
[43]
Briggs, N. et al. Realizing 2D Materials Via MOCVD. ECS Transactions 75, 725 (2016)
2016
-
[44]
Allerbeck, J. et al. Efficient and Continuous Carrier-Envelope Phase Control for Terahertz Lightwave-Driven Scanning Probe Microscopy. ACS Photonics 10, 3888–3895 (2023)
2023
-
[45]
E., Wei, Y., Everett, N., Jelic, V
Ammerman, S. E., Wei, Y., Everett, N., Jelic, V. & Cocker, T. L. Algorithm for subcycle terahertz scanning tunneling spectroscopy. Physical Review B105, 115427 (2022)
2022
-
[46]
Zhebrak, E. D. A method for calculating the franck–condon factors in the tomographic prob- ability representation of quantum mechanics. J Russ Laser Res37, 123––132 (2016)
2016
-
[47]
Zhao, W. et al. Lattice dynamics in mono- and few-layer sheets of WS2 and WSe2. Nanoscale 5, 9677–9683 (2013)
2013
-
[48]
Jeong, T. Y. et al. Coherent Lattice Vibrations in Mono- and Few-Layer WSe2. ACS Nano 10, 5560–5566 (2016)
2016
-
[49]
et al.Atomically Thin Heterostructures Based on Single-Layer Tungsten Diselenide and Graphene
Lin, Y.-C. et al.Atomically Thin Heterostructures Based on Single-Layer Tungsten Diselenide and Graphene. Nano Letters 14, 6936–6941 (2014)
2014
-
[50]
Krane, N. et al. Mapping the perturbation potential of metallic and dipolar tips in tunneling spectroscopy on MoS2. Physical Review B100, 035410 (2019). 33
2019
Reviewed August 11, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.