Pith. sign in

REVIEW 3 major objections 4 minor 13 references

Disentangling the components of a multiconfigurational excited state in isolated chromophore

T0 review · 3 major / 4 minor · reviewed 2026-08-07 · deepseek-v4-flash

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

Pith's one-line read Submolecular photocurrent maps of a single PTCDA anion are claimed to directly image the two dominant electronic configurations of its multiconfigurational first excited doublet state, with bias choosing which configuration recombines.

desk verdict Bias-switchable photocurrent fingerprints on a single open-shell molecule are a real advance, but the 'prove' language overreaches a non-unique forward model. read the letter →

arxiv 2502.09347 v1 pith:Y7NWPT2E submitted 2025-02-13 physics.atm-clus

Rodrigo Cezar de Campos Ferreira , Amandeep Sagwal , Jiří Doležal , Tomáš Neuman , Martin Švec This is my paper · ORCID
classification physics.atm-clus PACS 68.37.Ef73.40.Gk
keywords light-STMphotocurrentmappingmulticonfigurationalexcitedstatePTCDAanionDysonorbitalsdoubletradicalplasmonicnanocavityrate-equationmodel
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

This paper claims that zero-bias photocurrent maps of a single PTCDA anion molecule, recorded with light-scanning tunneling microscopy in a plasmonic nanocavity, are direct real-space fingerprints of the two dominant electronic configurations of its multiconfigurational first excited doublet state, $D_1^-$. The positive part of the bidirectional current is traced to electron injection into orbital 1 (configuration B), while the negative part is traced mostly to electron extraction from orbital 4 (configuration A); the two orbital signatures appearing in the same map are taken as proof that $D_1^-$ is a superposition of both configurations. A rate-equation model fed by TD-DFT transition densities and Kohn-Sham orbital slices reproduces the contrast and predicts that the bias voltage smoothly selects which recombination pathway dominates, which the experiments confirm by converting the map from an orbital-1 pattern to an orbital-3/4 pattern as the bias crosses zero. If correct, this gives experimental access, at the single-molecule level, to the correlated composition of excited states that previously could only be inferred from ensemble-averaged spectra.

What carries the argument

The central object is a set of Dyson-like orbital maps: for every tunneling transition between many-body states, the tip-mediated rate is taken proportional to the squared wave function of the Kohn-Sham orbital that approximates the transition's Dyson orbital, evaluated on a constant-height plane about 5 Å above the molecule ($z_T = 9.4$ a.u.). Each of the five many-body states ($S_0$, $D_0^-$, $D_1^-$, $D_2^-$, $S_0^{2-}$) is assigned a definite frontier-orbital occupation, so a charge transfer between states is tied to one specific orbital. Those two-dimensional densities are multiplied by a bias-dependent threshold function derived from a Franck-Condon model of the NaCl substrate reorganization ($E_R = 850$ meV), and the steady-state solution of the resulting rate equations gives the position- and bias-dependent net current that is compared, pixel by pixel, with the measured $\mathrm{d}I/\mathrm{d}P$ maps.

What would settle it

Compute true Dyson orbitals for the $D_1^- \to S_0$ and $D_1^- \to S_0^{2-}$ transitions with a correlated many-body method, including the NaCl substrate in the model, and compare their constant-height slices at $z_T = 9.4$ a.u. with the Kohn-Sham orbital 1 and orbital 3/4 slices used here; if the dominant Dyson densities differ appreciably from those slices, the claimed identification of the positive-current lobes with configuration B and the negative-current lobes with configuration A is falsified. A purely experimental second check: re-record the zero-bias photocurrent map at a substantially different tip height, since the orbital fingerprints should persist if they are genuine tunneling channels rather than features of the single chosen plane.

Watch

Extended reading notes

Core claim

In a single PTCDA anion on a NaCl/Ag(111) surface inside a plasmonic nanocavity, illumination at 785 nm produces a tunneling current at zero bias whose direction flips between positive and negative over distances of a few bond lengths. The paper explains this bidirectional photocurrent with a five-state model: absorption populates the $D_1^-$ and, through internal conversion, the $D_2^-$ states, and the excited molecule decays back to $D_0^-$ by two competing sequences of single-electron transfers, one passing through the neutral $S_0$ and the other through the doubly negative $S_0^{2-}$. Because $D_1^-$ is a superposition of two configurations (A, with an electron in orbital 4, and B, with a hole in orbital 1), each decay path carries a distinct spatial fingerprint given by the tunneling probability of the orbital involved in its rate-limiting step. The positive current lobes match the tunneling pattern of orbital 1 and the negative lobes match orbitals 3 and 4, so the maps are read as proving that both configurations contribute to the same excited state; applying a few hundred millivolts of bias progressively suppresses one path and enhances the other, turning the visualization into a controlled switch. The authors note that the relative weights of the two configurations are not fixed by this scheme, since a different choice of weights can reproduce the contrast; what is established is the coexistence of both occupations, not their coefficients.

Load-bearing premise

The entire spatial assignment rests on approximating each tunneling transition's Dyson orbital by a single Kohn-Sham orbital evaluated on a flat plane about 5 Å above the molecule, so if the true Dyson orbitals differ substantially from those orbital slices, the matching of map features to configurations A and B would collapse.

Editorial extensions

If this is right

  • A zero-bias photocurrent map on a single open-shell molecule can serve as a real-space readout of which orbitals participate in an excited state's decay, and not merely as a measure of total photocurrent.
  • Positive bias above about 150 mV pushes the decay through the $S_0^{2-}$ path, and the map collapses onto the orbital 1 pattern, isolating configuration B; negative bias below about $-150$ mV favors the $S_0$ path and exposes the orbital 3/4 pattern, isolating configuration A.
  • Under illumination the $\mathrm{d}I/\mathrm{d}V$ maps develop structured, location-dependent features inside the transport gap, showing that the optically populated states open tunneling channels that are absent in the dark.
  • The level scheme extracted from $\mathrm{d}I/\mathrm{d}V$ onsets and TEPL spectra, combined with the reorganization threshold function, reproduces the bias values at which the photocurrent maps switch character.

Reading between the lines

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

  • The paper does not exploit the quantitative side of the fingerprints: if the maps are proportional to Dyson-orbital densities, the relative intensity of positive and negative lobes should encode the configuration weights (TD-DFT here gives 74% for the $1\to 2$ transition and 26% for $2\to 4$), turning the coexistence proof into a measurement of superposition coefficients.
  • Because the rate equations already contain substrate reorganization, the same mapping should be able to track how configuration weights change when the molecule sits on different decoupling layers or couples to different metallic environments.
  • The bias-controlled selection of one decay path over the other is, in miniature, the same competition that governs charge separation at donor-acceptor interfaces in organic photovoltaics, so the technique offers a single-molecule test bed for how configuration mixing steers the direction of photoinduced charge transfer.
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

3 major / 4 minor

Summary. The paper reports light-STM measurements of zero-bias and bias-dependent photon-induced currents on individual PTCDA anions on NaCl/Ag(111), combined with TEPL and dI/dV spectroscopy. The central claim is that submolecular photocurrent maps contain spatial fingerprints of two dominant electronic configurations, denoted A and B, of the multiconfigurational doublet excited state D1^-, and that applying a bias voltage switches the dominant recombination pathway between the two configurations. The interpretation is developed through a rate-equation model in which tip- and substrate-mediated tunneling rates are expressed in terms of Kohn-Sham orbital slices used as approximate Dyson orbitals, with many-body level energies extracted from the same experimental TEPL and dI/dV data. The authors conclude that the correspondence between experimental and simulated maps 'proves' the presence of both configurations in D1^-.

Significance. If the central claim were established, the work would be a substantial advance: it would provide direct real-space visualization of the individual electronic configurations of a multiconfigurational excited state in a single molecule, together with a way to select one recombination pathway over another by bias voltage. The experimental data are novel, the combination of submolecular photocurrent mapping with a rate-equation model is ambitious, and the paper explicitly addresses a difficult open problem in single-molecule spectroscopy. However, the proof-level claim is not supported by the present analysis: the forward model is admitted to be non-unique, the spatial fingerprints rest entirely on an unvalidated Dyson-orbital approximation, and several key parameters are fitted to the same data that the model then reproduces. With a quantitative robustness analysis and appropriately qualified conclusions, the paper could become a valuable contribution; in its current form, the evidentiary gap between 'consistent with' and 'proves' is too wide.

major comments (3)
  1. [Main text, final paragraph of the photocurrent section; SI, 'A different parametrization ...'] The statement that the signatures of orbital 1 in the positive current and orbital 4 in the negative current 'prove the presence of both configurations A and B in the state D1^-' is load-bearing and is not justified by the analysis. The SI explicitly concedes that 'a different parametrization of the model could yield similar results even if different contributions of the two configurations were considered,' and no quantitative fit metric or uniqueness study is provided for the comparison in Fig. 2e. Moreover, the negative current is also generated by orbital 3 of D2^-, so the negative-lobe pattern does not by itself single out orbital 4. The data are consistent with the two-configuration model, but they do not uniquely determine the configuration weights; the word 'prove' should be replaced by a weaker claim, or supplemented by a sensitivity analysis over the fitted parameters and configuration weights.
  2. [SI, 'Finally, to obtain the spatial dependency ...', Eq. (4), z_T = 9.4 a.u.] The spatial fingerprints used to identify configurations A and B rest entirely on approximating each Dyson orbital by a single Kohn-Sham orbital evaluated on a constant-height plane and normalized to its own maximum. The manuscript provides no validation of this approximation and no test of robustness to the choice of plane height, the DFT functional, the fractional molecular charge (-0.5e) used in the Octopus calculation, or the expected difference between true Dyson orbitals and single Kohn-Sham orbitals. Since the identification of the negative-current feature with orbital 4 and configuration A depends on this approximation, a robustness check is required before the assignment can be regarded as established.
  3. [Table S1 and SI Eqs. (3)-(14); Fig. 2d,e] The model re-uses the same experimental data from which its input energies and rates are derived, so the agreement in Fig. 2e and the bias-dependent maps is not an independent validation. The many-body energies in Table S1 are extracted from the TEPL and dI/dV spectra that the model then simulates, and eta0 is explicitly selected to match the ratio of photocurrent maxima to sequential-tunneling-current maxima. This does not invalidate the model as an interpretation tool, but it means that the correspondence between experiment and simulation cannot, by itself, carry the burden of proving the configuration assignment. The authors should clearly separate fitted parameters from predicted quantities and, where possible, test whether a single parameter set held fixed across all maps and bias voltages is sufficient.
minor comments (4)
  1. [Fig. S5 caption] The caption uses the symbol omega2 for both the bias modulation frequency and the laser modulation frequency; one of the two should be labeled omega1 to avoid ambiguity.
  2. [Table S1] The many-body energies are quoted without uncertainties; since they are derived from experimental spectra, an estimate of the error bars would help assess how well-determined the level alignment is.
  3. [SI Eq. (4)] Normalizing the DOS by its maximum removes the absolute orbital amplitude from the tunneling rate; this choice should be stated explicitly and justified, because it affects the relative weights of the different tunneling channels.
  4. [References] Some references contain incomplete bibliographic data, for example reference [15] lacks a final page/article number and several Nature-family entries list only volume and starting page; these should be completed for publication.

Circularity Check

0 steps flagged · score 1.0 of 10

No circular reduction in the central inference; configuration weights come from TD-DFT, not from the photocurrent maps, though the 'prove' wording overstates a non-unique model.

full rationale

The central claim—that the D1- state contains configurations A and B—is not obtained by fitting the photocurrent maps to those configurations. The composition is taken from TD-DFT transition analysis (SI: 'the composition of the D0-D1 transition that dominantly consists of the transitions 2->4 (26%) and 1->2 (74%)'), and the spatial fingerprints (orbitals 1 and 4) are computed Kohn-Sham/Dyson-orbital slices used in the rate-equation model. The experimental maps are then compared with the simulated maps, so the match is a consistency check with independent TD-DFT input, not a reduction of the conclusion to the input by construction. Several parameters (eta0, alpha, gamma0_T, gamma0_S, kappa0, kappa_pl, gamma_IC, E_R) are fitted to the same experimental spectra or to the current-amplitude ratio, but these fits calibrate the model; they do not by themselves determine the spatial A/B assignment, which rests on the independently calculated orbital shapes. The SI candidly states the model's non-uniqueness ('a different parametrization ... could yield similar results even if different contributions of the two configurations were considered'), and the Dyson-orbital approximation is an assumption; these are limitations on proof strength and correctness, not circularity. Self-citations (e.g., the threshold-function line-shape model from 'our earlier work', SI Ref. [3]) are used as modeling tools but are not load-bearing uniqueness constraints. No equation or fitted parameter is renamed as a prediction of the central claim, so no specific circular reduction can be exhibited.

Assumptions & free parameters 11 free parameters · 3 assumptions · 0 invented entities

No new physical entities are introduced. The main burdens on the reader are the many fitted parameters and the modeling assumptions that convert molecular orbital shapes into simulated current maps. These parameters and assumptions carry the weight of the comparison between theory and experiment, and the paper does not provide independent quantitative validation of them.

free parameters (11)
  • gamma0^T (tip tunneling rate scale) = 1 µeV
    Sets the overall magnitude of tip-molecule tunneling; chosen so that tip tunneling dominates over substrate tunneling.
  • gamma0^S (substrate tunneling rate scale) = 0.1 µeV
    Chosen relative to gamma0^T without independent measurement.
  • alpha (fraction of voltage drop across molecule-tip gap) = 0.8
    Justified by 'a good match of the simulation with experimental data, particularly the dI/dV curves' (SI after Eq. 5).
  • kappa0 (background plasmon-mediated decay rate) = 0.1 meV
    Selected consistently with previous studies; no independent measurement.
  • kappa_pl (tip-position-dependent decay rate) = 0.5 meV
    Selected consistently with previous studies.
  • z1, z2 (plasmon potential localization heights) = 30 a.u., 40 a.u.
    Described as 'suitable chosen parameters approximating the localization of the plasmonic field' (SI Eq. 10).
  • eta0 (laser pumping rate) = 3e-7
    Selected so that the ratio of photocurrent maxima to sequential tunneling current maxima corresponds to the experiment (SI after Eq. 11).
  • gamma_IC (internal conversion rate) = 2 meV
    Chosen to reproduce the experimentally observed population of D2^-.
  • E(S0) and E(S0^2-) (many-body state energies) = 0.200 eV and 0.380 eV
    Extracted from experimental dI/dV spectra using the threshold function (SI Table S1 and text).
  • E(D1^-) and E(D2^-) = 1.330 eV and 1.456 eV
    Taken from TEPL emission peak energies, i.e., from the same experiment being modeled.
  • E_R (reorganization energy) = 850 meV
    Carried from previous work (Sellies et al.), not fitted here.
assumptions (3)
  • domain assumption The five-state rate-equation model with no spin degeneracy and with triplet states, photon-assisted tunneling, and hot electrons omitted
    Stated in the main text: 'Within the framework of this simplified picture, in which the neutral or doubly negative triplet states, photon-assisted tunneling, and role of hot electrons are omitted'. If any of these omitted channels contributes significantly, the assignment of photocurrent contrast to specific configurations A and B could be wrong.
  • ad hoc to paper Tunneling rates factor as gamma0 x DOS(x,y) x Theta(Delta E), with DOS approximated by the squared Kohn-Sham orbital evaluated on a plane above the molecule
    SI Eqs. (3)-(4). This factorization and the plane evaluation are specific modeling choices that directly determine the spatial fingerprints attributed to each configuration.
  • domain assumption The substrate reorganization is described by a threshold function built from an optical-phonon spectral density with E_R = 850 meV and phonon bounds 18-31 meV
    SI Eqs. (12)-(14). This model is carried from earlier work and is a critical component for reproducing the gradual bias dependence of the photocurrent, but it is not independently validated here.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Disentangling the components of a multiconfigurational excited state in isolated chromophore." pith.science (2026). https://pith.science/paper/Y7NWPT2E

@misc{pith2026250209347,
  author       = {Pith},
  title        = {Pith review of: Disentangling the components of a multiconfigurational excited state in isolated chromophore},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/Y7NWPT2E}},
  note         = {Machine review of arXiv:2502.09347}
}
read the original abstract

Studying the excited states of doublets is challenging for their typically multiconfigurational character. We employ light-scanning-tunneling microscopy (light-STM) to investigate photon-induced currents on a single open-shell PTCDA anion molecule placed into a plasmonic nanocavity between a tip and a substrate, irradiated by laser. Submolecular mapping reveals a zero-bias bidirectional photocurrent strongly varying with the lateral position of the tip apex above the molecule. We elucidate the mechanism in terms of a theoretical model in which a multiconfigurational doublet state is excited and decays back to the anion ground state through sequential electron transfers with the tip and the substrate. The correspondence of the experimental and theoretical contrast proves the correlated character of the excited state which can be described as a superposition of two dominating electronic configurations. By applying bipolar voltage on the junction with the molecule, we switch the dominant recombination pathway from one of the configurations to the other, effectively disentangling the multiconfigurational state individual components through visualization of their Dyson orbitals, as corroborated by theoretical modelling.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

13 extracted references · 11 canonical work pages

  1. [1]

    Kaiser, S

    K. Kaiser, S. Jiang, M. Romeo, F. Scheurer, G. Schull, and A. Rosławska, Gating Single-Molecule Fluorescence with Electrons, Physical Review Letters 133, (2024)

    work page 2024

  2. [2]

    Jiang, T

    S. Jiang, T. Neuman, R. Bretel, A. Boeglin, F. Scheurer, E. Le Moal, and G. Schull, Many-Body Description of STM-Induced Fluorescence of Charged Molecules, Phys. Rev. Lett. 130, (2023)

    work page 2023

  3. [3]

    Vasilev, S

    K. Vasilev, S. Canola, F. Scheurer, A. Boeglin, F. Lotthammer, F. Chérioux, T. Neuman, and G. Schull, Exploring the Role of Excited States’ Degeneracy on Vibronic Coupling with Atomic-Scale Optics, ACS Nano 18, 28052 (2024)

    work page 2024

  4. [4]

    Sellies, J

    L. Sellies, J. Eckrich, L. Gross, A. Donarini, and J. Repp, Controlled Single-Electron Transfer Enables Time-Resolved Excited-State Spectroscopy of Individual Molecules, Nature Nanotechnology (2024)

    work page 2024

  5. [5]

    A. D. Becke, Density-functional thermochemistry. III. The role of exact exchange, The Journal of Chemical Physics 98, 5648 (1993)

    1993

  6. [6]

    M. J. Frisch, G. W. Trucks, H. B. Schlegel, G. E. Scuseria, M. A. Robb, J. R. Cheeseman, G. Scalmani, V. Barone, G. A. Petersson, H. Nakatsuji, X. Li, M. Caricato, A. V. Marenich, J. Bloino, B. G. Janesko, R. Gomperts, B. Mennucci, H. P. Hratchian, J. V. Ortiz, A. F. Izmaylov, J. L. Sonnenberg, D. Williams-Young, F. Ding, F. Lipparini, F. Egidi, J. Goings...

    work page 2016

  7. [7]

    Grimme, S

    S. Grimme, S. Ehrlich, and L. Goerigk, Effect of the damping function in dispersion corrected density functional theory, J Comput Chem 32, 1456 (2011)

    work page 2011

  8. [8]

    Lu, A comprehensive electron wavefunction analysis toolbox for chemists, Multiwfn, The Journal of Chemical Physics 161, (2024)

    T. Lu, A comprehensive electron wavefunction analysis toolbox for chemists, Multiwfn, The Journal of Chemical Physics 161, (2024)

    work page 2024

Show all 13 references
  1. [9]

    N. Tancogne-Dejean et al., Octopus, a computational framework for exploring light-driven phenomena and quantum dynamics in extended and finite systems, The Journal of Chemical Physics 152, (2020)

    2020

  2. [10]

    J. P. Perdew and A. Zunger, Self-interaction correction to density-functional approximations for many-electron systems, Phys. Rev. B 23, 5048 (1981)

    1981

  3. [11]

    P. A. M. Dirac, Note on Exchange Phenomena in the Thomas Atom, Math. Proc. Camb. Phil. Soc. 26, 376 (1930)

    1930

  4. [12]

    Slater, A Simplification of the Hartree-Fock Method, Phys

    C. Slater, A Simplification of the Hartree-Fock Method, Phys. Rev. 81, 385 (1951)

    1951

  5. [13]

    Imai-Imada, H

    M. Imai-Imada, H. Imada, K. Miwa, J. Jung, T. K. Shimizu, M. Kawai, and Y. Kim, Energy-level alignment of a single molecule on ultrathin insulating film, Phys. Rev. B 98, (2018)

    2018

Pith tools

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