arxiv: 2603.03783 · v1 · submitted 2026-03-04 · ❄️ cond-mat.mes-hall

Recognition: no theorem link

Imaging asymmetric Coulomb blockade phenomena across metallic nanoislands

Junho Bang , Byeongin Lee , Hankyu Lee , Jian-Feng Ge , Doohee Cho

Authors on Pith no claims yet

Pith reviewed 2026-05-15 16:59 UTC · model grok-4.3

classification ❄️ cond-mat.mes-hall

keywords Coulomb blockadenanoislandsscanning tunneling spectroscopywork functionindiumblack phosphorussingle-electron transportorthodox theory

0 comments

The pith

Work function differences between junctions produce shifted and asymmetric Coulomb blockade resonances across indium nanoislands.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper maps how charging resonances in metallic nanoislands vary with position using scanning tunneling spectroscopy on indium islands grown on black phosphorus. The resonances disperse spatially, with their symmetry axis offset in bias voltage and their curvature becoming asymmetric. Orthodox theory calculations reproduce both the offset and the asymmetry once distinct work functions are assigned to the tip-island and island-substrate junctions. This link shows that junction-specific electrostatics, rather than uniform charging, set the detailed shape of the Coulomb blockade spectrum.

Core claim

In indium nanoislands on semiconducting black phosphorus, spatially resolved spectra exhibit charging resonances whose trajectories display a finite shift of the symmetry axis together with pronounced asymmetric curvature. Calculations based on orthodox theory match these features when differences in work function between the two junctions are included, establishing that junction-specific electrostatics control the nanoscale charge transport spectrum.

What carries the argument

Orthodox theory model of Coulomb blockade that incorporates distinct work functions for the tip-nanoisland and nanoisland-substrate junctions.

If this is right

Junction work function must be treated as a tunable parameter when predicting or engineering single-electron spectra.
Spatial STM maps of charging resonances become a direct probe of local junction electrostatics.
Asymmetric curvature implies directional preference in electron tunneling rates across the island.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The same junction asymmetry could be engineered deliberately to create single-electron rectifiers.
In other STM studies of nanoislands, apparent quantum confinement signals may partly reflect unaccounted work function gradients.
Substrate choice offers a route to preset the offset of the Coulomb blockade diamond without gate electrodes.

Load-bearing premise

That the observed spatial dispersion and asymmetry are fully explained by orthodox theory once only work-function differences are introduced.

What would settle it

A set of spectra recorded on islands with deliberately matched work functions in both junctions that still show the same shifts and curvature.

Figures

Figures reproduced from arXiv: 2603.03783 by Byeongin Lee, Doohee Cho, Hankyu Lee, Jian-Feng Ge, Junho Bang.

**Figure 2.** Figure 2: FIG. 2. Spatially resolved Coulomb blockade spectra. (a), Normalized differential conductance [PITH_FULL_IMAGE:figures/full_fig_p013_2.png] view at source ↗

**Figure 3.** Figure 3: FIG. 3. Spatially dispersive Coulomb peaks in an In nanoisland. (a) Spatially resolved normalized [PITH_FULL_IMAGE:figures/full_fig_p014_3.png] view at source ↗

**Figure 4.** Figure 4: FIG. 4. Coulomb blockade asymmetry from work function differences. (a)–(c) Simulated [PITH_FULL_IMAGE:figures/full_fig_p015_4.png] view at source ↗

read the original abstract

Coulomb blockade (CB) arises in nanoscale systems with ultra-small capacitance, where discrete charging effects dictate electron transport, enabling wide-ranging applications based on single-electron transistors. Despite established electrostatic control of charge states in quantum dots and nanoislands, a rigorous quantitative link between junction parameters and the CB spectrum remains elusive. Here, using scanning tunneling spectroscopy, we investigate the spatial variation of CB in indium nanoislands on semiconducting black phosphorus. We observe spatially dispersive charging resonances whose trajectories exhibit a finite shift of the symmetry axis in bias as well as a pronounced asymmetric curvature. By comparing the experimental results with calculations based on orthodox theory, we show that these features originate from work function differences in the junctions, underscoring the importance of junction-specific electrostatics in nanoscale charge transport.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper maps spatially varying asymmetric Coulomb blockade on indium nanoislands and shows that a single work-function offset in orthodox theory reproduces the bias shifts and curvature.

read the letter

The main takeaway is that STM spectroscopy on these indium nanoislands reveals charging resonances that disperse across the island and show a clear offset in bias plus asymmetric curvature. The authors run orthodox Coulomb blockade calculations that include different work functions for the tip-island and island-substrate junctions and find that this single adjustment lines up with the observed trajectories. That spatial mapping plus the direct comparison to theory is the concrete new piece; similar asymmetries have been noted before but not imaged this way on nanoislands with explicit junction parameters. The data presentation looks clean and the modeling is straightforward, which is useful for anyone building single-electron devices where junction electrostatics matter. The soft spot is that the work-function difference is fitted to the data rather than measured independently, so the match shows consistency but does not exclude other spatially varying effects such as local capacitance gradients or substrate band bending. Without error bars on the fit or a test that rules out alternatives, the attribution is reasonable but not uniquely diagnostic. This is the kind of incremental but practical result that people in mesoscopic transport and 2D-material devices would want to see. It is solid enough to deserve a serious referee who can check the full quantitative comparisons and any additional controls in the manuscript.

Referee Report

2 major / 1 minor

Summary. The manuscript reports scanning tunneling spectroscopy measurements of Coulomb blockade (CB) resonances in indium nanoislands on black phosphorus. It observes spatially dispersive charging lines whose trajectories in bias voltage exhibit a finite offset of the symmetry axis together with pronounced asymmetric curvature. The authors compare these data to orthodox-theory calculations and conclude that the observed spatial dispersion and asymmetry originate from work-function differences between the tip-island and island-substrate junctions.

Significance. If the attribution is robust, the work supplies a concrete illustration that junction-specific electrostatic parameters (here, work-function offsets) can dominate the detailed shape of CB spectra even in nominally metallic islands. The spatial mapping capability adds a useful experimental handle for separating local electrostatic effects from intrinsic island properties, which is relevant for the design and modeling of single-electron devices.

major comments (2)

[Abstract and theoretical-modeling section] Abstract and theoretical-modeling section: the claim that orthodox-theory calculations 'reproduce the features once work-function differences are included' is presented without quantitative fit metrics (e.g., rms deviation of resonance positions, error bars on extracted offsets, or goodness-of-fit statistics). Because the work-function difference is the sole adjustable parameter invoked to match both the bias-axis shift and the curvature asymmetry, the absence of these diagnostics leaves open whether the agreement is predictive or achieved by construction.
[Results section (comparison of experiment and theory)] Results section (comparison of experiment and theory): the manuscript does not report an independent experimental determination of the work-function offset (e.g., from separate Kelvin-probe or photoemission data on the same junctions). Without such a constraint, alternative spatially varying electrostatic contributions (local capacitance gradients, substrate band bending, or tip-induced potentials) cannot be quantitatively excluded as sources of equivalent trajectories under orthodox theory.

minor comments (1)

[Figures] Figure captions and axis labels should explicitly state the range of work-function differences explored in the calculations and whether any other parameters were held fixed.

Simulated Author's Rebuttal

2 responses · 1 unresolved

We thank the referee for the constructive comments, which have prompted us to strengthen the quantitative presentation and discussion of alternative explanations in the manuscript. We address each major comment below.

read point-by-point responses

Referee: [Abstract and theoretical-modeling section] Abstract and theoretical-modeling section: the claim that orthodox-theory calculations 'reproduce the features once work-function differences are included' is presented without quantitative fit metrics (e.g., rms deviation of resonance positions, error bars on extracted offsets, or goodness-of-fit statistics). Because the work-function difference is the sole adjustable parameter invoked to match both the bias-axis shift and the curvature asymmetry, the absence of these diagnostics leaves open whether the agreement is predictive or achieved by construction.

Authors: We agree that quantitative fit metrics are essential to demonstrate the robustness of the modeling. In the revised manuscript we have added the root-mean-square deviation between measured and calculated resonance positions for multiple spatial cuts, together with uncertainty estimates on the single fitted work-function offset obtained from a global fit across the dataset. These metrics confirm that one offset value reproduces the observed shifts and curvature asymmetries with low RMS error, indicating the agreement is not achieved by construction. revision: yes
Referee: [Results section (comparison of experiment and theory)] Results section (comparison of experiment and theory): the manuscript does not report an independent experimental determination of the work-function offset (e.g., from separate Kelvin-probe or photoemission data on the same junctions). Without such a constraint, alternative spatially varying electrostatic contributions (local capacitance gradients, substrate band bending, or tip-induced potentials) cannot be quantitatively excluded as sources of equivalent trajectories under orthodox theory.

Authors: We acknowledge that an independent experimental constraint on the work-function offset would be desirable. Such measurements are not feasible within the present STM setup on the identical nanoscale junctions. In the revision we have expanded the discussion to show that the principal alternatives (capacitance gradients or substrate band bending) would produce qualitatively different spatial trajectories or bias asymmetries than those observed; the data are instead consistent with a uniform offset across islands. We therefore maintain that the work-function interpretation remains the most parsimonious. revision: partial

standing simulated objections not resolved

Independent experimental determination of the work-function offset (e.g., via Kelvin-probe or photoemission on the same junctions) is not available in the current dataset.

Circularity Check

1 steps flagged

Work-function difference fitted to reproduce CB asymmetry; agreement with orthodox theory is by construction

specific steps

fitted input called prediction [Abstract; theory-comparison paragraphs]
"By comparing the experimental results with calculations based on orthodox theory, we show that these features originate from work function differences in the junctions"

The work-function offset is introduced as the sole adjustable parameter in the orthodox CB model and is tuned to reproduce the measured trajectories, finite bias shift, and asymmetric curvature. Because the model contains only this one free parameter and no external calibration of the offset is given, the match is achieved by construction; the data are reproduced once the input is chosen to fit them.

full rationale

The paper's central claim—that asymmetric CB features originate from junction work-function differences—is supported only by inserting an adjustable offset into orthodox-theory calculations and tuning it to match the observed spatial dispersion, bias-axis shift, and curvature. No independent measurement of the offset is reported, so the reproduction demonstrates consistency with some value of the parameter rather than an independent prediction. This matches the 'fitted_input_called_prediction' pattern exactly.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The interpretation relies on orthodox single-electron theory plus the assumption that work-function differences are the dominant source of asymmetry; no new entities are introduced.

free parameters (1)

work function difference between tip-island and island-substrate junctions
Adjusted to reproduce the observed bias-axis shift and curvature asymmetry.

axioms (1)

domain assumption Orthodox theory of Coulomb blockade applies without significant quantum or substrate effects varying across the island
Invoked when comparing experimental trajectories to calculated spectra.

pith-pipeline@v0.9.0 · 5437 in / 1276 out tokens · 22459 ms · 2026-05-15T16:59:04.894406+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

38 extracted references · 38 canonical work pages

[1]

M. A. Nielsen and I. L. Chuang,Quantum computation and quantum information(Cambridge univer- sity press, 2010)

work page 2010
[2]

M. H. Devoret and R. J. Schoelkopf, Amplifying quantum signals with the single-electron transistor, Nature406, 1039 (2000)

work page 2000
[3]

J. P. Pekola, O.-P. Saira, V . F. Maisi, A. Kemppinen, M. Möttönen, Y . A. Pashkin, and D. V . Averin, Single-electron current sources: Toward a refined definition of the ampere, Reviews of Modern Physics 85, 1421 (2013)

work page 2013
[4]

R. S. Shaikhaidarov, K. H. Kim, J. W. Dunstan, I. V . Antonov, S. Linzen, M. Ziegler, D. S. Golubev, V . N. Antonov, E. V . Il’ichev, and O. V . Astafiev, Quantized current steps due to the ac coherent quantum phase-slip effect, Nature608, 45 (2022)

work page 2022
[5]

M. H. Devoret, D. Esteve, and C. Urbina, Single-electron transfer in metallic nanostructures, Nature 360, 547 (1992). 8

work page 1992
[6]

Shnirman and G

A. Shnirman and G. Schoen, Quantum measurements performed with a single-electron transistor, Physical Review B57, 15400 (1998)

work page 1998
[7]

Schoelkopf, P

R. Schoelkopf, P. Wahlgren, A. Kozhevnikov, P. Delsing, and D. Prober, The radio-frequency single- electron transistor (RF-SET): A fast and ultrasensitive electrometer, Science280, 1238 (1998)

work page 1998
[8]

M. Yoo, T. Fulton, H. Hess, R. Willett, L. Dunkleberger, R. Chichester, L. Pfeiffer, and K. West, Scanning single-electron transistor microscopy: Imaging individual charges, Science276, 579 (1997)

work page 1997
[9]

Vlaic, S

S. Vlaic, S. Pons, T. Zhang, A. Assouline, A. Zimmers, C. David, G. Rodary, J.-C. Girard, D. Roditchev, and H. Aubin, Superconducting parity effect across the Anderson limit, Nature Com- munications8, 14549 (2017)

work page 2017
[10]

Trivini, J

S. Trivini, J. Ortuzar, K. Vaxevani, B. Viña-Bausá, F. S. Bergeret, and J. I. Pascual, Local control of parity and charge in nanoscale superconducting lead islands, Physical Review Letters135, 226001 (2025)

work page 2025
[11]

Flensberg, A

K. Flensberg, A. A. Odintsov, F. Liefrink, and P. Teunissen, Towards single-electron metrology, Inter- national Journal of Modern Physics B13, 2651 (1999)

work page 1999
[12]

J. P. Pekola, J. J. Vartiainen, M. Möttönen, O.-P. Saira, M. Meschke, and D. V . Averin, Hybrid single- electron transistor as a source of quantized electric current, Nature Physics4, 120 (2008)

work page 2008
[13]

F. A. Zwanenburg, A. S. Dzurak, A. Morello, M. Y . Simmons, L. C. Hollenberg, G. Klimeck, S. Rogge, S. N. Coppersmith, and M. A. Eriksson, Silicon quantum electronics, Reviews of Modern Physics85, 961 (2013)

work page 2013
[14]

Elzerman, R

J. Elzerman, R. Hanson, L. Willems van Beveren, B. Witkamp, L. Vandersypen, and L. P. Kouwen- hoven, Single-shot read-out of an individual electron spin in a quantum dot, Nature430, 431 (2004)

work page 2004
[15]

T. A. Fulton and G. J. Dolan, Observation of single-electron charging effects in small tunnel junctions, Physical Review Letters59, 109 (1987)

work page 1987
[16]

C. W. Beenakker, Theory of Coulomb-blockade oscillations in the conductance of a quantum dot, Physical Review B44, 1646 (1991)

work page 1991
[17]

Lafarge, P

P. Lafarge, P. Joyez, D. Esteve, C. Urbina, and M. Devoret, Two-electron quantization of the charge on a superconductor, Nature365, 422 (1993)

work page 1993
[18]

Takahashi and S

S. Takahashi and S. Maekawa, Effect of Coulomb blockade on magnetoresistance in ferromagnetic tunnel junctions, Physical Review Letters80, 1758 (1998)

work page 1998
[19]

Van Bentum, R

P. Van Bentum, R. Smokers, and H. Van Kempen, Incremental charging of single small particles, Physical Review Letters60, 2543 (1988). 9

work page 1988
[20]

Deshpande, W

A. Deshpande, W. Bao, Z. Zhao, C. Lau, and B. LeRoy, Imaging charge density fluctuations in graphene using Coulomb blockade spectroscopy, Physical Review B83, 155409 (2011)

work page 2011
[21]

M. T. Woodside and P. L. McEuen, Scanned probe imaging of single-electron charge states in nanotube quantum dots, Science296, 1098 (2002)

work page 2002
[22]

I.-P. Hong, C. Brun, M. Pivetta, F. Patthey, and W.-D. Schneider, Coulomb blockade phenomena observed in supported metallic nanoislands, Frontiers in Physics1, 13 (2013)

work page 2013
[23]

Jarillo-Herrero, S

P. Jarillo-Herrero, S. Sapmaz, C. Dekker, L. P. Kouwenhoven, and H. S. Van Der Zant, Electron-hole symmetry in a semiconducting carbon nanotube quantum dot, Nature429, 389 (2004)

work page 2004
[24]

D. L. Klein, R. Roth, A. K. Lim, A. P. Alivisatos, and P. L. McEuen, A single-electron transistor made from a cadmium selenide nanocrystal, Nature389, 699 (1997)

work page 1997
[25]

Negishi, T

R. Negishi, T. Hasegawa, K. Terabe, M. Aono, H. Tanaka, T. Ogawa, and H. Ozawa, IV characteris- tics of single electron tunneling from symmetric and asymmetric double-barrier tunneling junctions, Applied Physics Letters90, 223112 (2007)

work page 2007
[26]

S. Wu, G. Nazin, X. Chen, X. Qiu, and W. Ho, Control of relative tunneling rates in single molecule bipolar electron transport, Physical Review Letters93, 236802 (2004)

work page 2004
[27]

Reiner, A

J. Reiner, A. K. Nayak, A. Tulchinsky, A. Steinbok, T. Koren, N. Morali, R. Batabyal, J.-H. Kang, N. Avraham, Y . Oreg,et al., Spectroscopic visualization of a robust electronic response of semicon- ducting nanowires to deposition of superconducting islands, Physical Review X10, 011002 (2020)

work page 2020
[28]

Hanna and M

A. Hanna and M. Tinkham, Variation of the Coulomb staircase in a two-junction system by fractional electron charge, Physical Review B44, 5919 (1991)

work page 1991
[29]

Wilkins, E

R. Wilkins, E. Ben-Jacob, and R. Jaklevic, Scanning-tunneling-microscope observations of Coulomb blockade and oxide polarization in small metal droplets, Physical Review Letters63, 801 (1989)

work page 1989
[30]

Schönenberger, H

C. Schönenberger, H. Van Houten, and H. Donkersloot, Single-electron tunnelling observed at room temperature by scanning-tunnelling microscopy, Europhysics Letters20, 249 (1992)

work page 1992
[31]

Averin and K

D. Averin and K. Likharev, inMesoscopic Phenomena in Solids(Elsevier, 1991) pp. 173–271

work page 1991
[32]

See Supplemental Material for further analysis and more detailed description of methods

work page
[33]

Schmeidel, H

J. Schmeidel, H. Pfnür, and C. Tegenkamp, Coulomb blockade effects in Ag/Si (111): The role of the wetting layer, Physical Review B80, 115304 (2009)

work page 2009
[34]

K. Gong, L. Zhang, W. Ji, and H. Guo, Electrical contacts to monolayer black phosphorus: A first- principles investigation, Physical Review B90, 125441 (2014)

work page 2014
[35]

L. P. Kouwenhoven, G. Schön, and L. L. Sohn, inMesoscopic Electron Transport(Springer, 1997) pp. 10 1–44

work page 1997
[36]

N. A. Pradhan, N. Liu, C. Silien, and W. Ho, Atomic scale conductance induced by single impurity charging, Physical Review Letters94, 076801 (2005)

work page 2005
[37]

Cockins, Y

L. Cockins, Y . Miyahara, S. D. Bennett, A. A. Clerk, S. Studenikin, P. Poole, A. Sachrajda, and P. Grutter, Energy levels of few-electron quantum dots imaged and characterized by atomic force microscopy, Proceedings of the National Academy of Sciences107, 9496 (2010)

work page 2010
[38]

J. Qian, B. I. Halperin, and E. J. Heller, Imaging and manipulating electrons in a one-dimensional quantum dot with Coulomb blockade microscopy, Physical Review B81, 125323 (2010). 11 CIS CTI RIS RTI Tip In BP Black Phosphorus In Islands STM Tip VB (V) (dI/dV)/(I/V) (arb. units) (a) (b) (c) −1−1 −0.5−0.5 00 0.50.5 110 0.01 0.02 0.03 0.04 −0.2−0.2 −0.1−0.1...

work page 2010