Pith. sign in

REVIEW 2 major objections 3 minor 52 references

Negative differential conductance in triangular molecular assemblies

T0 review · 2 major / 3 minor · reviewed 2026-08-05 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read A triangular cluster of three TBTAP molecules on lead produces negative differential conductance between 0.7 and 0.9 V, driven by inter-molecular Coulomb correlations rather than superconductivity.

desk verdict Solid experimental NDC in a molecular trimer with spatial Coulomb-ring mapping; the PME 'quantitative reproduction' rests on an explicitly questionable weak-coupling assumption that the stress-test correctly identifies as unprotected for the NDC mechanism itself. read the letter →

arxiv 2508.05575 v1 pith:L5TIZH5R submitted 2025-08-07 cond-mat.mes-hall

Chao Li , Vladislav Pokorný , Prokop Hapala , Martin Žonda , Ping Zhou , Silvio Decurtins , Shi-Xia Liu , Fengqi Song
show 2 more authors
Rémy Pawlak Ernst Meyer
This is my paper · ORCID
classification cond-mat.mes-hall PACS 73.23.Hk73.63.-b85.65.+h
keywords negativedifferentialconductanceCoulombblockademoleculartrimerthree-impurityAndersonmodelPaulimasterequationcapacitivecouplingYu–Shiba–Rusinovstatesscanningtunnelingspectroscopy
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper reports a three-molecule device that shows negative differential conductance (NDC): between 0.7 and 0.9 V the current drops as the voltage rises. The device is a C3-symmetric triangular trimer of TBTAP molecules assembled on a Pb(111) surface by lateral STM manipulation, and the authors reproduce its behavior quantitatively with a three-impurity Anderson model solved by Pauli master equations. The mechanism is capacitive: the inter-molecular coupling ($W = 50$ meV) makes the two-electron charge state of the cluster more stable than the three-electron state, and under bias the cluster is driven into a non-equilibrium mixture of singly-occupied states in which the empty site often sits far from the tip, suppressing the tunneling current. The authors argue the effect is pure electron-correlation physics, unrelated to the superconducting substrate, and that the same molecules can be re-programmed by cluster geometry: a linear tetramer behaves as a memory, the triangular trimer as a molecular analogue of a Gunn diode.

What carries the argument

The three-impurity Anderson model (TIAM): three molecular sites at $-90$ meV with on-site repulsion $U = 200$ meV, inter-site capacitive coupling $W = 50$ meV, and negligible direct hopping, weakly coupled to tip and substrate. Large $U$ forbids double occupancy, shrinking the Hilbert space from 64 to 8 states and making the dynamics solvable by Pauli master equations. The tip acts two ways: its electrostatic field (plate capacitor with spherical apex and mirror charge) shifts site energies; its tunneling amplitude decays exponentially with site distance, $t_{ti} \propto \exp(-\beta|r_i - r_t|)$. Molecular orbital shape (the SOMO with a fitted quadrupole moment) reproduces the chiral pattern

What would settle it

Two checks would settle the claim. (1) Re-run the Pauli master equation with the experimentally fitted substrate coupling $\Gamma_s = 10$\u2013$20$ meV and compare the simulated dI/dV maps with the measured discharge rings and NDC patterns; if the rings or the NDC region change shape or vanish, the electrostatic-only assertion fails. (2) Measure the trimer on a normal-metal substrate (e.g., Ag(111)): the claim that NDC is independent of superconductivity predicts the same 0.7–0.9 V NDC there, and its absence would refute that independence.

Watch

Extended reading notes

Core claim

Central claim: the NDC between 0.7 and 0.9 V in a C3-symmetric trimer of TBTAP molecules on Pb(111) is a Coulomb-blockade effect driven by inter-molecular capacitive coupling, purely electron-correlation physics, independent of the superconducting substrate. Many sharp peaks in the dI/dV maps are intra-cluster charge rearrangements, not total-charge changes. The NDC appears between excited-state crossings, where singly occupied states (e.g., $|100\rangle$) fall below doubly occupied states; without inter-site hopping, electrons cannot relax, occupancy spreads across the cluster, and tunneling is suppressed because the empty site is often far from the tip. The non-superconducting model reprod

Load-bearing premise

The load-bearing premise is that the discharge-ring geometry and the NDC pattern depend only on the local electrostatic environment, so the simulations can use a substrate coupling of $\Gamma_s \sim 0.01$ meV even though the YSR fit gives $\Gamma_s = 10$\u2013$20$ meV. The paper states this without derivation in the transport-calculations section, calling the weak-coupling assumption 'questionable'; if the electrostatic-only assertion is wrong, the quantitative reproduction o

Editorial extensions

If this is right

  • The same molecular building block implements different devices by geometry alone: a linear tetramer acts as a bistable memory, while the triangular trimer acts as a molecular analogue of a Gunn diode, with no change in chemistry.
  • Because the NDC is independent of the superconducting substrate, the trimer should function on normal-metal surfaces, cleanly separating the Coulomb-blockade physics from YSR and superconductivity physics.
  • The eight-state Hilbert space of the trimer makes it a fully simulable unit, offering a tractable building block for designing molecular logic and charge-state computing architectures.
  • Charge-state readout of such clusters cannot rely on counting dI/dV peaks as total-charge toggles: many peaks signal charge redistribution among sites at constant total charge.
  • The hexamer demonstrates electrostatic gating of inner sites by outer molecules, extending the design space from three-site clusters toward larger, programmable arrays.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the mechanism is right, the NDC voltage window (0.7–0.9 V) should be tunable by modifying the inter-site capacitance $W$ — for example, by changing molecular spacing, linker chemistry, or substrate screening — since $W$ sets the separation between doubly and triply occupied configurations.
  • The non-equilibrium-occupancy mechanism predicts similar NDC in other frustrated geometries of strongly coupled charge sites with negligible inter-site hopping, such as triangular arrays of different radical molecules; the essential ingredients are capacitive coupling and threefold frustration, not the specific molecule.
  • The paper's assertion that discharge-ring geometry is independent of the substrate coupling $\Gamma_s$ is directly testable: at higher temperatures the ring edges should broaden with $k_B T$ rather than with $\Gamma_s$, and cotunneling signatures would appear if the weakly coupled assumption breaks down.
  • The YSR-validated parameters imply a doubly occupied, spin-carrying trimer ground state; the same platform might therefore combine NDC switching at high bias with spin-sensitive transport at low bias, extending toward the hybrid superconducting devices the authors mention as an outlook.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

2 major / 3 minor

Summary. The paper reports the assembly of a triangular trimer of TBTAP radical molecules on Pb(111) and its characterization by low-temperature STM/STS. The central experimental observation is a negative differential conductance (NDC) region between roughly 0.7 and 0.9 V, accompanied by voltage-dependent discharge rings in dI/dV maps. The authors model the trimer with a three-impurity Anderson model, solve the equilibrium low-energy physics with NRG, and simulate the finite-bias transport with a Pauli master equation (PME). They attribute the NDC to non-equilibrium occupation of singly occupied states caused by strong inter-site capacitive coupling and negligible direct hopping. A hexamer structure is also presented as a topological variant. The claim is that the NDC is purely electronic, independent of superconductivity, and quantitatively reproduced by the model.

Significance. If the theoretical account is correct, the work would be a valuable demonstration of a molecular-scale, gate-tunable NDC element whose functionality is controlled by cluster geometry rather than by the molecule-lead interface. The experimental data are of high quality: the manipulation is documented, the discharge-ring evolution is systematic, the chiral patterns are reproducible, and the finite Hilbert space of the trimer makes the modeling transparent. The authors also provide code links and use open-source packages, which is a strength. However, the manuscript's central quantitative claim rests on a PME calculation performed in a coupling regime that the authors themselves concede is questionable and that is inconsistent with the coupling used in their own YSR/NRG analysis. This weakens, but does not necessarily invalidate, the proposed NDC mechanism.

major comments (2)
  1. [Supplementary Text, Eq. (S7)-(S8); main text 'Negative differential conductance emerging from non-equilibrium occupancy] The PME calculation sets Γs≈0.01 meV ≪ kBTexp, explicitly outside the Γs=10–20 meV inferred from YSR fits [27] and used in the NRG calculation of Figs. S10G/S11. The authors state that the ring size is independent of Γ, but that statement concerns the equilibrium charge-state locus, not the NDC mechanism. The NDC in Fig. 3B is a dynamical steady-state effect: it requires the population to be spread over the singly occupied states |100⟩, |010⟩, |001⟩ via Eq. (S7), and those populations depend on Γs and Γt. No calculation is shown for Γs≈10–20 meV; at that coupling the substrate relaxation rate exceeds kBT/ℏ and the non-equilibrium occupancy should be strongly suppressed. Please provide a calculation in the broadened-coupling regime (or a convincing argument, e.g., an adiabatic elimination or finite-Γ master equation) showing that the NDC survives, or reframe the claim accordingly.
  2. [Main text 'Triangular trimer assemblies and their discharging behavior'; Fig. 2 and Fig. S7] The phrase 'quantitatively reproduces the experimental observations' is stronger than the evidence supports. The simulation maps in Fig. 2A’-H’ are in arbitrary units and the model parameters (εid=-90 meV, W=50 meV, t=0, the tip electrostatic parameters rtip, ztip, zC, zS, and the quadrupole moment Qxx) are fitted to the same discharge-ring data that the model then reproduces. The YSR NRG check (Fig. S10G/S11) is not an independent validation because it uses the same U, W, t and ε≈-104 meV inferred from the charging-ring analysis. I recommend presenting the comparison as a self-consistent parametrization or semi-quantitative reproduction and adding a sensitivity analysis for at least W and the tip parameters.
minor comments (3)
  1. [Figure 4 caption and panels] The caption lists negative sample voltages (Vs = −0.37 V, −0.39 V, ...) but the panel labels show 'VS = 0.37 V' etc. The missing minus signs should be corrected for clarity.
  2. [Abstract and main text] The phrase 'quantitatively reproduces' should be reconciled with the arbitrary-unit simulations in Fig. 2 and Fig. S7; otherwise the claim is misleading.
  3. [Throughout] Several typographical/spacing errors remain, e.g., 'an molecular system', 'e fficiently', 'di fferential'. A careful proofread is needed.

Circularity Check

1 steps flagged · score 6.0 of 10

Ring/NDC 'reproduction' is a fit to the same Fig. S7 maps; the PME uses Γs≈0.01 meV against the authors' own YSR value of 10–20 meV, and the stated Γ-independence justification covers ring geometry, not the NDC mechanism.

  1. fitted input called prediction [Supplementary Text, 'Pauli master equations' (p. S6); see also main text p. 4 and Figs. 2, S7]
    "We obtained the parameters of the model by fitting the experimental data plotted in Fig. S7. The fitted parameters of the electrostatic model read rtip = 0.3 nm, ztip = 0.6 nm, zC = 2 nm and zS = −0.09 nm."

    The same Fig. S7 data set—experimental dI/dV maps from Vs = 0.47 to 0.87 V, including the NDC regions—is later presented as being reproduced by the PME calculation (Figs. 2A'–H' and Fig. S7 columns 2 and 4). Since rtip, ztip, zC, zS, the quadrupole, W = 50 meV, and t = 0 (also fitted: 'The best fit of the experimental data was obtained for the zero value of inter-site hopping t') were all tuned to these maps, the agreement is a goodness-of-fit, not a prediction. The NDC pattern and the t = 0 mechanism invoked to explain it are consequences of the fitted model, so the central claim that the calculation 'quantitatively reproduces' the NDC reduces to the fit.

full rationale

The paper contains one genuine circular step: the PME parameters (tip electrostatics rtip, ztip, zC, zS, quadrupole, W = 50 meV, t = 0) are fitted to the experimental dI/dV maps of Fig. S7, and the same maps are then exhibited (Fig. 2, Fig. S7/S8) as a 'quantitative reproduction' and used to 'pinpoint physical parameters.' The central NDC pattern is part of the fitted dataset, so the agreement is training-set performance, not an independent prediction. This is a fitted-input-called-prediction circularity. The YSR NRG calculation is a genuine out-of-sample check on independent low-energy spectra, but because it uses parameters inferred from the ring fit (ε≈−104 meV, U=200, W=50, Γ=20 meV), it is a consistency check rather than an independent derivation of the NDC mechanism. The authors explicitly flag the PME coupling assumption as questionable: Γs∼0.01 meV is two to three orders below Γs,exp=10–20 meV inferred from their own YSR fits [27]. Their defense that ring size/shape depend only on electrostatics addresses the locus of charge transitions, not the non-equilibrium NDC mechanism, which is governed by the master-equation populations (Eq. S7) that depend on Γs and Γt. This is a robustness/correctness concern, not itself circularity, but it compounds the fitted nature of the central simulation. Nevertheless, the model has independent content—explicit many-body Hamiltonian, master-equation dynamics, hexamer transfer—so the paper is not wholly derivative. Score 6 reflects one central prediction reducing to a fit; it is not 8–10 because no self-citation chain forces the result and the YSR channel provides partial external constraint.

Assumptions & free parameters 9 free parameters · 6 assumptions · 0 invented entities

The paper introduces no new physical particles or forces. All additional structure is in the form of fitted model parameters and an ad hoc electrostatic tip model, which are captured as free parameters and domain assumptions.

free parameters (9)
  • epsilon_id (impurity level) = -90 meV
    Assumed from the broadening peak at -90 mV in Fig. S2, so it is set from experimental data, not derived.
  • U (on-site Coulomb) = 200 meV (infinite in PME)
    Taken from previous TBTAP studies and then approximated as infinite in the master equation; not independently derived here.
  • W (inter-site capacitive coupling) = 50 meV
    Fitted to reproduce the experimental discharge-ring data in Fig. S7; also said to be consistent with previous dimer fits.
  • t (direct hopping) = 0
    Obtained as the best fit to the experimental data, simplifying the many-body basis.
  • Gamma_s (substrate coupling in PME) = 0.01 meV
    Chosen to stay in the sequential tunneling regime, but inconsistent with the experimental value 10-20 meV from YSR fits; the authors acknowledge the assumption is questionable.
  • beta (tunneling decay constant) = not specified
    Tunable parameter in the exponential model t_ti proportional to exp(-beta |ri-rt|); no fixed value given.
  • tip electrostatic parameters (rtip, ztip, zC, zS) = rtip=0.3 nm, ztip=0.6 nm, zC=2 nm, zS=-0.09 nm
    Fitted to match the evolution of charging rings with voltage; zS does not correspond to the physical surface position.
  • Qxx (molecular quadrupole moment) = not specified
    Fitted to reproduce the elliptical shape of the discharge rings.
  • epsilon (local energy in YSR NRG) = -104 meV
    Used in the NRG calculation of YSR spectra (Fig. S10), close to but not identical to the -90 meV used in PME.
assumptions (6)
  • domain assumption Each molecule is described by a single correlated quantum level with spin.
    The three-impurity Anderson model in Eq. S2 reduces each TBTAP molecule to one orbital; vibrational and multi-orbital effects are neglected.
  • domain assumption Double occupancy is prohibited (infinite U) in the PME calculation.
    Reduces the Hilbert space from 64 to 8 states; justified by comparing with finite U=200 meV results in Fig. S6, not by a microscopic derivation.
  • domain assumption Direct inter-site hopping t is negligible.
    Set to zero as the best fit; the physical justification is weak and the model relies on the substrate for charge rearrangement.
  • ad hoc to paper System is in the sequential tunneling regime (Gamma_s, Gamma_t << kBT).
    Underlies the Pauli master equation; the authors themselves note this is questionable given the experimental Gamma_s = 10-20 meV.
  • domain assumption The tip electrostatic field is modeled as a parallel-plate capacitor with a single spherical adatom and an electrostatic mirror at zS.
    The model is introduced in the supplement to fit the ring evolution; the mirror position zS is unphysical.
  • domain assumption Tunneling amplitudes decay exponentially with tip-molecule distance and the tip wavefunction is s-like.
    Standard approximation in STM theory; used to compute the spatially dependent tunneling rates.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Negative differential conductance in triangular molecular assemblies." pith.science (2026). https://pith.science/paper/L5TIZH5R

@misc{pith2026250805575,
  author       = {Pith},
  title        = {Pith review of: Negative differential conductance in triangular molecular assemblies},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/L5TIZH5R}},
  note         = {Machine review of arXiv:2508.05575}
}
read the original abstract

We report the creation and characterization of a molecular-scale negative differential conductance (NDC) device by assembling a triangular trimer of 4,5,9,10-tetrabromo-1,3,6,8-tetraazapyrene (TBTAP) molecules on a superconducting Pb(111) substrate. Using low-temperature scanning tunneling spectroscopy, we observe robust NDC behavior manifesting as a decrease in current with increasing voltage between 0.7-0.9 V arising from the interplay of Coulomb blockade and strong inter-molecular capacitive coupling within the molecular cluster. Gate-controlled charging and discharging processes are directly visualized via two-dimensional differential conductance mapping, which reveals the emergence of Coulomb rings and spatial regions of NDC. Theoretical modeling using a three-impurity Anderson model and master equation approach quantitatively reproduces the experimental observations and demonstrates that the NDC emerges purely from electron correlations, independent of the underlying superconductivity. By tuning the geometry to a hexamer structure, we further show that cluster topology provides versatile control over electronic properties at the molecular scale. These results establish a functional platform for implementing multifunctional molecular devices and highlight a strategy toward programmable and scalable nanoelectronics.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

52 extracted references · 52 canonical work pages

  1. [27]

    Li, et al., Individual Assembly of Radical Molecules on Superconductors: Demonstrating Quantum Spin Behavior and Bistable Charge Rearrangement

    C. Li, et al., Individual Assembly of Radical Molecules on Superconductors: Demonstrating Quantum Spin Behavior and Bistable Charge Rearrangement. ACS Nano 19 (3), 3403–3413 (2025)

    work page 2025

  2. [1]

    Xiang, X

    D. Xiang, X. Wang, C. Jia, T. Lee, X. Guo, Molecular-Scale Electronics: From Concept to Function. Chem. Rev. 116 (7), 4318–4440 (2016)

    work page 2016

  3. [2]

    Evers, R

    F. Evers, R. Koryt ´ar, S. Tewari, J. M. van Ruitenbeek, Advances and challenges in single- molecule electron transport. Rev. Mod. Phys. 92, 035001 (2020)

    work page 2020

  4. [3]

    Famili, et al., Self-Assembled Molecular-Electronic Films Controlled by Room Temper- ature Quantum Interference

    M. Famili, et al., Self-Assembled Molecular-Electronic Films Controlled by Room Temper- ature Quantum Interference. Chem 5 (2), 474–484 (2019)

    work page 2019

  5. [4]

    Bedrossian, D

    P. Bedrossian, D. M. Chen, K. Mortensen, J. A. Golovchenko, Demonstration of the tunnel- diode effect on an atomic scale. Nature 342 (6247), 258–260 (1989)

    work page 1989

  6. [5]

    I.-W. Lyo, P. Avouris, Negative Differential Resistance on the Atomic Scale: Implications for Atomic Scale Devices. Science 245 (4924), 1369–1371 (1989)

    work page 1989

  7. [6]

    B. Xu, Y . Dubi, Negative di fferential conductance in molecular junctions: an overview of experiment and theory. J. Phys. Condens. Matter 27 (26), 263202 (2015)

    work page 2015

  8. [7]

    Single-photon detection enabled by negative differential conductivity in moir\'e superlattices

    K. Nowakowski, et al., Single-photon detection enabled by negative differential conductivity in moir´e superlattices (2025), arXiv:2505.13637

    work page Pith review arXiv 2025

Show all 52 references
  1. [8]

    B. W. Knight, G. A. Peterson, Theory of the Gunn E ffect. Phys. Rev. 155, 393–404 (1967)

    1967

  2. [9]

    Esaki, P

    L. Esaki, P. J. Stiles, New Type of Negative Resistance in Barrier Tunneling. Phys. Rev. Lett. 16, 1108–1111 (1966)

    1966

  3. [10]

    M. C. Rogge, F. Cavaliere, M. Sassetti, R. J. Haug, B. Kramer, Negative conductance patterns of quantum dots: experiment and theory. New J. Phys. 8 (12), 298 (2006)

    2006

  4. [11]

    Kostyrko, B

    T. Kostyrko, B. R. Bulka, Symmetry-controlled negative di fferential resistance e ffect in a triangular molecule. Phys. Rev. B 79, 075310 (2009)

    2009

  5. [12]

    Wrze ´sniewski, I

    K. Wrze ´sniewski, I. Weymann, Dark states in spin-polarized transport through triple quantum dot molecules. Phys. Rev. B 97, 075425 (2018)

    2018

  6. [13]

    K. S. Kim, et al. , Visualizing Atomic-Scale Negative Di fferential Resistance in Bilayer Graphene. Phys. Rev. Lett. 110, 036804 (2013)

    2013

  7. [14]

    He, et al

    G. He, et al. , Negative Differential Conductance & Hot-Carrier Avalanching in Monolayer WS2 FETs. Sci. Rep. 7 (1), 11256 (2017)

    2017

  8. [15]

    Zhao, et al., Electric-Field-Driven Negative Differential Conductance in 2D van der Waals Ferromagnet Fe3GeTe2

    M. Zhao, et al., Electric-Field-Driven Negative Differential Conductance in 2D van der Waals Ferromagnet Fe3GeTe2. Nano Lett. 21 (21), 9233–9239 (2021)

    2021

  9. [16]

    Grobis, A

    M. Grobis, A. Wachowiak, R. Yamachika, M. F. Crommie, Tuning negative differential resis- tance in a molecular film. Appl. Phys. Lett. 86 (20), 204102 (2005)

    2005

  10. [17]

    Chen, et al., Mechanism for Negative Differential Resistance in Molecular Electronic De- vices: Local Orbital Symmetry Matching

    L. Chen, et al., Mechanism for Negative Differential Resistance in Molecular Electronic De- vices: Local Orbital Symmetry Matching. Phys. Rev. Lett. 99, 146803 (2007). 12

    2007

  11. [18]

    Fern ´andez-Torrente, D

    I. Fern ´andez-Torrente, D. Kreikemeyer-Lorenzo, A. Str ´o˙zecka, K. J. Franke, J. I. Pascual, Gating the Charge State of Single Molecules by Local Electric Fields. Phys. Rev. Lett. 108, 036801 (2012)

    2012

  12. [19]

    H. B. Heersche, et al. , Electron Transport through Single Mn 12 Molecular Magnets. Phys. Rev. Lett. 96, 206801 (2006)

    2006

  13. [20]

    M. L. Perrin, et al., Large negative differential conductance in single-molecule break junc- tions. Nat. Nanotechnol. 9 (10), 830–834 (2014)

    2014

  14. [21]

    C. S. Lent, B. Isaksen, M. Lieberman, Molecular Quantum-Dot Cellular Automata. J. Am. Chem. Soc. 125 (4), 1056–1063 (2003)

    2003

  15. [22]

    Bandyopadhyay, R

    A. Bandyopadhyay, R. Pati, S. Sahu, F. Peper, D. Fujita, Massively parallel computing on an organic molecular layer. Nat. Phys. 6 (5), 369–375 (2010)

    2010

  16. [23]

    Berger, et al

    J. Berger, et al. , Quantum dissipation driven by electron transfer within a single molecule investigated with atomic force microscopy. Nat. Commun. 11 (1) (2020)

    2020

  17. [24]

    C. Li, Z. Wang, Y . Lu, X. Liu, L. Wang, Conformation-based signal transfer and processing at the single-molecule level. Nat. Nanotech. 12 (11), 1071–1076 (2017)

    2017

  18. [25]

    Li, et al., Strong signature of electron-vibration coupling in molecules on Ag(111) trig- gered by tip-gated discharging

    C. Li, et al., Strong signature of electron-vibration coupling in molecules on Ag(111) trig- gered by tip-gated discharging. Nat. Commun. 14 (1), 5956 (2023)

    2023

  19. [26]

    Pawlak, et al., Towards Gate-Tunable Topological Superconductivity in a Supramolecular Electron Spin Lattice (2025), in review

    R. Pawlak, et al., Towards Gate-Tunable Topological Superconductivity in a Supramolecular Electron Spin Lattice (2025), in review

    2025

  20. [28]

    M. R. Wasielewski, et al., Exploiting chemistry and molecular systems for quantum informa- tion science. Nat. Rev. Chem. 4, 490–504 (2020)

    2020

  21. [29]

    Bairamkulov, G

    R. Bairamkulov, G. De Micheli, Superconductive Electronics: A 25-Year Review [Feature]. IEEE Circuits Syst. Mag. 24 (2), 16–33 (2024)

    2024

  22. [30]

    S. Wu, G. Nazin, X. Chen, X. Qiu, W. Ho, Control of relative tunneling rates in single molecule bipolar electron transport. Phys. Rev. Lett. 93 (23), 236802 (2004)

    2004

  23. [31]

    Kumar, C

    D. Kumar, C. Krull, Y . Yin, N. V . Medhekar, A. Schiffrin, Electric field control of molecular charge state in a single-component 2D organic nanoarray. ACS Nano 13 (10), 11882–11890 (2019)

    2019

  24. [32]

    N. A. Pradhan, N. Liu, C. Silien, W. Ho, Atomic Scale Conductance Induced by Single Im- purity Charging. Phys. Rev. Lett. 94, 076801 (2005)

    2005

  25. [33]

    G. V . Nazin, X. H. Qiu, W. Ho, Charging and Interaction of Individual Impurities in a Mono- layer Organic Crystal. Phys. Rev. Lett. 95, 166103 (2005). 13

    2005

  26. [34]

    Wong, et al., Characterization and manipulation of individual defects in insulating hexag- onal boron nitride using scanning tunnelling microscopy.Nat

    D. Wong, et al., Characterization and manipulation of individual defects in insulating hexag- onal boron nitride using scanning tunnelling microscopy.Nat. Nanotechnol. 10 (11), 949–953 (2015)

    2015

  27. [35]

    Koci ´c, et al

    N. Koci ´c, et al. , Periodic Charging of Individual Molecules Coupled to the Motion of an Atomic Force Microscopy Tip. Nano Lett. 15 (7), 4406–4411 (2015)

    2015

  28. [36]

    Li, et al., Imaging local discharge cascades for correlated electrons in WS 2/WSe2 moir´e superlattices

    H. Li, et al., Imaging local discharge cascades for correlated electrons in WS 2/WSe2 moir´e superlattices. Nat. Phys. 17 (10), 1114–1119 (2021)

    2021

  29. [37]

    Manikandan, P

    M. Manikandan, P. Nicolini, P. Hapala, Computational Design of Photosensitive Polymer Templates To Drive Molecular Nanofabrication.ACS Nano 18 (14), 9969–9979 (2024)

    2024

  30. [38]

    S. G. Balasubramani, et al., TURBOMOLE: Modular Program Suite for Ab Initio Quantum- Chemical and Condensed-Matter Simulations. J. Chem. Phys. 152 (18), 184107 (2020)

    2020

  31. [39]

    Kaufmann, K

    J. Kaufmann, K. Held, ana cont: Python Package for Analytic Continuation. Comput. Phys. Commun. 282, 108519 (2023)

    2023

  32. [40]

    Bergeron, A.-M

    D. Bergeron, A.-M. S. Tremblay, Algorithms for Optimized Maximum Entropy and Diag- nostic Tools for Analytic Continuation. Phys. Rev. E 94, 023303 (2016)

    2016

  33. [41]

    ˇZitko, NRG Ljubljana (2021), doi:10.5281 /zenodo.4841076, https://doi.org/10

    R. ˇZitko, NRG Ljubljana (2021), doi:10.5281 /zenodo.4841076, https://doi.org/10. 5281/zenodo.4841076

    2021

  34. [42]

    Kir ˇsanskas, J

    G. Kir ˇsanskas, J. N. Pedersen, O. Karlstr ¨om, M. Leijnse, A. Wacker, QmeQ 1.0: An open- source Python package for calculations of transport through quantum dot devices. Comput. Phys. Commun. 221, 317–342 (2017)

    2017

  35. [43]

    Hapala, et al., ppafm: Probe Particle AFM Simulation Code — PhotonMap Branch (2025), https://github.com/Probe-Particle/ppafm/tree/PhotonMap, accessed: 2025-07- 02

    P. Hapala, et al., ppafm: Probe Particle AFM Simulation Code — PhotonMap Branch (2025), https://github.com/Probe-Particle/ppafm/tree/PhotonMap, accessed: 2025-07- 02

    2025

  36. [44]

    Argento, R

    C. Argento, R. H. French, Parametric tip model and force–distance relation for Hamaker constant determination from atomic force microscopy. J. Appl. Phys. 80 (11), 6081–6090 (1996)

    1996

  37. [45]

    Baier, C

    R. Baier, C. Leendertz, M. C. Lux-Steiner, S. Sadewasser, Toward quantitative Kelvin probe force microscopy of nanoscale potential distributions. Phys. Rev. B 85 (16) (2012)

    2012

  38. [46]

    D. Z. Gao, et al. , Using Metallic Noncontact Atomic Force Microscope Tips for Imaging Insulators and Polar Molecules: Tip Characterization and Imaging Mechanisms. ACS Nano 8 (5), 5339–5351 (2014)

    2014

  39. [47]

    Gross, et al., Investigating atomic contrast in atomic force microscopy and Kelvin probe force microscopy on ionic systems using functionalized tips

    L. Gross, et al., Investigating atomic contrast in atomic force microscopy and Kelvin probe force microscopy on ionic systems using functionalized tips. Phys. Rev. B 90 (15) (2014)

    2014

  40. [48]

    Krej ˇc´ı, P

    O. Krej ˇc´ı, P. Hapala, M. Ondr´aˇcek, P. Jel´ınek, Principles and simulations of high-resolution STM imaging with a flexible tip apex. Phys. Rev. B 95 (4) (2017)

    2017

  41. [49]

    C. J. Chen, Tunneling matrix elements in three-dimensional space: The derivative rule and the sum rule. Phys. Rev. B 42 (14), 8841–8857 (1990). 14

    1990

  42. [50]

    Dole ˇzal, et al., Real Space Visualization of Entangled Excitonic States in Charged Molec- ular Assemblies

    J. Dole ˇzal, et al., Real Space Visualization of Entangled Excitonic States in Charged Molec- ular Assemblies. ACS Nano 16 (1), 1082–1088 (2021)

    2021

  43. [51]

    Oinonen, et al., Advancing scanning probe microscopy simulations: A decade of develop- ment in probe-particle models

    N. Oinonen, et al., Advancing scanning probe microscopy simulations: A decade of develop- ment in probe-particle models. Comput. Phys. Commun. 305, 109341 (2024)

    2024

  44. [52]

    Goldozian, F

    B. Goldozian, F. A. Damtie, G. Kir ˇsanskas, A. Wacker, Transport in serial spinful multiple- dot systems: The role of electron-electron interactions and coherences. Sci. Rep. 6 (1), 22761 (2016). Acknowledgments V .P., M.ˇZ and P.H. thank Tom´aˇs Novotn´y for many fruitful di...

    2016

Pith tools

Reviewed August 5, 2026 · model on record in the stance chip above.