REVIEW 3 major objections 5 minor 42 references
Vibrational excitations in magnetic triangular nanographenes
T0 review · 3 major / 5 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read The inelastic tunneling steps at ±22.5 mV and ±35 mV in a single phenalenyl radical on Au(111) are off-resonant vibrational excitations of three out-of-plane modes, not spin excitations.
desk verdict Clear evidence that the IETS steps in phenalenyl are vibrational rather than spin excitations, with a plausible mode assignment that needs one honest paragraph about a 2.5 meV discrepancy. 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 central object is the vibrationally perturbed SOMO wavefunction, $\delta\Psi_{m,k} = \delta Q_k (\partial \Psi_m/\partial Q_k)$, whose overlap with an s-wave tip orbital enters the tunneling matrix element. For each mode $k$, the equilibrium structure is displaced by $\pm\delta Q_k/2$, the resulting change in the SOMO is computed with DFT, and the tunneling matrix element gives the inelastic intensity. The mechanism is off-resonant vibration-assisted tunneling, with strong intensity only when the mode's out-of-plane displacement has the same local symmetry as the SOMO phase: the degenerate $k=1,2$ modes (20.0 meV) create a threefold triangular ring with a node in the center, while the $k=4$ mode (35.4 meV) is centered on the molecule.
What would settle it
Measure phenalenyl on a weakly coupled decoupling layer, such as NaCl on Cu(111), and repeat the below/above-step map subtraction. If the step energies shift by more than a few meV or the maps lose the predicted triangular-ring and center patterns, the gas-phase mode assignment and the vibrational mechanism would fail; conversely, if the maps were to follow the Kondo/spin-density positions, a spin origin would be indicated.
Extended reading notes
Core claim
The paper establishes that the symmetric differential-conductance steps observed at $\pm 22.5$ mV and $\pm 35$ mV on a single phenalenyl radical on Au(111) are off-resonant vibrational excitations, not spin excitations. A single $S=1/2$ radical has no spin-flip excited state in zero magnetic field, and the maps of step intensity do not follow the singly occupied molecular orbital (SOMO) or its Kondo signature; instead, DFT-based simulations of the vibrationally perturbed wavefunction reproduce the measured maps. The steps are assigned to the degenerate out-of-plane modes $k=1,2$ at 20.0 meV and to mode $k=4$ at 35.4 meV, and the agreement of the below/above-step map subtraction with simulation is the decisive evidence.
Load-bearing premise
The calculations treat the molecule as isolated in the gas phase, so the argument assumes that sitting on Au(111) does not significantly change the vibrational modes or their coupling to the electron, even though the surface strongly screens the unpaired spin.
Editorial extensions
If this is right
- In IETS of a single $S=1/2$ nanographene radical, symmetric conductance steps cannot be spin excitations, and the spatial maps shown here identify them as vibrational.
- The inelastic intensity is controlled by the vibrationally perturbed wavefunction, not the elastic orbital density, so IETS maps can be concentrated where the elastic signal is suppressed.
- Only vibrational modes whose out-of-plane motion shares the local symmetry of the SOMO phase couple strongly; for phenalenyl these are exactly the 20.0 meV and 35.4 meV modes.
- Choosing a tip position at the molecule's edge suppresses vibrational background, which is a practical guideline for measuring spin excitations in triangulene-based systems.
- The below/above-step map-subtraction protocol yields directly the spatial distribution of the inelastic tunneling probability and can be applied to other molecular adsorbates.
Reading between the lines
- Editorial inference: the same map-subtraction protocol should separate vibrational from magnetic IETS steps in triangulene dimers and chains, where both types of excitation can appear at similar energies.
- Editorial inference: if substrate hybridization does not strongly renormalize the modes, the 20.0 meV and 35.4 meV energies serve as a fingerprint for identifying phenalenyl on other metal surfaces.
- Editorial inference: the symmetry-matching rule suggests that for larger triangulenes the number of strongly coupled vibrational modes will remain small, which should simplify the interpretation of spin IETS on those systems.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports low-temperature STM/IETS measurements on the phenalenyl radical (S = 1/2) adsorbed on Au(111). The authors observe symmetric conductance steps at bias voltages of ±22.5 mV and ±35 mV. Using the absence of spin excitations for a spin-1/2 center, the spatial distribution of the step intensities, and DFT-based Bardeen simulations of vibration-assisted tunneling, they assign the features to three out-of-plane vibrational modes: a degenerate pair at 20.0 meV (k = 1, 2) and a mode at 35.4 meV (k = 4). They conclude that the steps are off-resonant vibrational excitations and that the approach can discriminate vibrational from magnetic IETS features in nanographenes.
Significance. If the conclusions hold, this work provides a practical method for distinguishing vibrational from spin excitations in IETS of magnetic nanographenes, which is an important issue for the field of pi-magnetism. The paper has clear strengths: the spin-1/2 argument rigorously rules out conventional spin-flip excitations; the DFT-based Bardeen simulations use no free parameters fitted to the IETS maps, with fixed tip parameters and standard B3LYP/6-311g(d,p); and the simulated mode energies and spatial maps are compared directly with experiment. The central vibrational-origin conclusion is well supported. However, the quantitative assignment of the 22.5 mV feature to the 20.0 meV modes is not adequately justified, and the gas-phase DFT assumption is not discussed in light of the strong molecule-substrate hybridization evidenced by the Kondo resonance.
major comments (3)
- [Main text, discussion of Fig. 1d and Fig. 3] The manuscript assigns the experimental step at ±22.5 mV to the calculated degenerate modes k = 1, 2 at 20.0 meV, while the ±35 mV step matches the k = 4 mode at 35.4 meV almost exactly. The 2.5 meV (12%) discrepancy for the lower-energy pair is never commented on. Since the IETS step position is a direct measure of the excitation energy of the adsorbed molecule, this mismatch is quantitatively material. A global rigid shift due to the substrate or DFT error cannot explain why one pair of modes matches to within 0.4 meV while the other is off by 2.5 meV. The authors should discuss the expected accuracy of the gas-phase vibrational energies, consider anharmonicity or surface-induced renormalization, or provide additional evidence that the 22.5 mV feature indeed corresponds to the k = 1, 2 pair rather than to some other mode. Without this discussion, the quantitative mode assignment for one of the two features is not yet established.
- [Methods (first paragraph) and Fig. 4 comparison] The DFT calculations are performed for an isolated gas-phase molecule, whereas the experiment is performed on a molecule strongly hybridized with Au(111), as evidenced by the pronounced Kondo resonance. The manuscript does not state what effect the surface is expected to have on the vibrational energies, the SOMO shape, or the electron-phonon coupling. This is a load-bearing assumption because both the mode energies and the simulated spatial maps are used for the assignment. The authors should either justify the transferability of gas-phase results to the adsorbed system or show that the conclusions are robust to, for example, different DFT functionals, basis sets, or a simplified substrate model. Such a test would directly address the correctness risk of the assignment.
- [Fig. 4c-e and Conclusion] The statement that the "conformity between simulated and measured maps accordingly proofs" that the steps are caused by off-resonant vibrational excitations is stronger than the evidence presented. The comparison is visual and qualitative, and for the 22.5 mV step the authors acknowledge that the subtraction procedure overestimates the Kondo background and creates six darker areas. While the central dark region is a genuine signal, a quantitative comparison metric (e.g., a normalized cross-correlation or a line-cut overlay) would strengthen the claim. This is not a fatal flaw, but it should be addressed so that the spatial evidence is as convincing as the spin-1/2 argument.
minor comments (5)
- [Text near Eq. (1)] There is a typo: "integrataion plane" should be "integration plane."
- [Paragraph discussing Fig. 2 and Kondo resonance] "Deconstructive interference" should be "destructive interference."
- [Section on k = 1, 2 modes (Fig. 4a discussion)] "Dislocations of the carbon atoms" should likely be "displacements of the carbon atoms," since the text refers to normal-mode motion, not lattice defects.
- [Conclusion] "Proofs" should be "proves" (or "demonstrates"). Also, "exctitation" in the sentence before Eq. (1) should be "excitation."
- [Fig. 4c caption and main text] The subtraction procedure for the maps (constant-height maps at different biases) is described only briefly. Stating whether the maps were normalized or drift-corrected before subtraction, and how the bias offsets were chosen, would improve reproducibility.
Circularity Check
No significant circularity: the mode assignment is tested against experimental spatial maps, and the DFT/Bardeen inputs are not fitted to the measured IETS data.
full rationale
The paper's central claim is that the IETS steps at ±22.5 mV and ±35 mV arise from three out-of-plane vibrational modes calculated with gas-phase DFT (B3LYP/6-311G(d,p)) and Bardeen tunneling matrix elements. The calculation inputs—work function 5 eV, tip height 7 Å, integration plane 1.6 Å—are fixed and not adjusted to reproduce the experimental maps. The measured spatial distributions in Fig. 4c are the external benchmark against which the simulated maps are compared; no fitted parameter is renamed as a prediction. The energy comparison (20.0 meV vs 22.5 mV and 35.4 meV vs 35 mV) also provides an independent quantitative check, even if the 2.5 meV offset for the lower mode is a possible accuracy concern rather than circularity. The only relevant self-citation (Ref. 41, Reecht et al., which shares an author) is used as a methodological reference for the electron-phonon coupling simulation; the method is stated in the paper and validated by the experimental comparison, so it is not load-bearing as a self-citation. Ref. 12, also with overlapping authors, establishes the spin-1/2 ground state but is independently supported by the Kondo signal and by Lieb/Ovchinnikov-type spin rules. No step in the derivation defines the predicted quantity in terms of the experimental observable or imports a uniqueness theorem from the authors' prior work.
Assumptions & free parameters
free parameters (3)
- Tip workfunction phi =
5 eV
- Tip height z_t =
7 A
- Integration plane position z_S =
1.6 A
assumptions (4)
- domain assumption Gas-phase DFT (B3LYP/6-311g(d,p)) describes the vibrational modes and SOMO wavefunction of phenalenyl on Au(111) accurately enough for mode assignment.
- domain assumption Bardeen's tunneling formalism with an s-wave tip at 7 A height and a planar integration surface captures the relative inelastic tunneling probabilities.
- standard math The off-resonant vibrational excitation probability is proportional to the tunneling matrix element of the vibrationally perturbed molecular orbital delta Psi_m,k (Lorente-Persson theory).
- domain assumption A single S=1/2 spin has no excited spin state in zero magnetic field, so the IETS steps cannot be spin excitations.
Cite this review
Pith. "Pith review of Vibrational excitations in magnetic triangular nanographenes." pith.science (2026). https://pith.science/paper/UG5HTMWU
@misc{pith2026241119670,
author = {Pith},
title = {Pith review of: Vibrational excitations in magnetic triangular nanographenes},
year = {2026},
howpublished = {\url{https://pith.science/paper/UG5HTMWU}},
note = {Machine review of arXiv:2411.19670}
}
read the original abstract
Inelastic electron tunneling spectroscopy (IETS) is a powerful measurement technique often used in scanning tunneling spectroscopy to probe excited states of various nanostructures, e.g., the magnetic properties of complex spin systems. The observed excited states can be of magnetic and vibrational origin and it is therefore necessary to differentiate between these two excitation mechanisms. Here, we investigate the spin S = 1/2 phenalenyl radical on Au(111). IETS measurements feature inelastic excitations, whereas the spatial distribution of their intensity excludes any spin excitations. Comparison to theoretical simulations proves the vibrational origin of those excitations and allows us to assign the observed features to distinct vibrational modes.
Figures
Reference graph
Works this paper leans on
-
[1]
Yazyev, O. V. Hyperfine Interactions in Graphene and Related Carbon Nanostructures. Nano Letters 2008, 8, 1011--1015
work page 2008
-
[2]
Min, H.; Hill, J. E.; Sinitsyn, N. A.; Sahu, B. R.; Kleinman, L.; MacDonald, A. H. Intrinsic and Rashba spin-orbit interactions in graphene sheets. Phys. Rev. B 2006, 74, 165310
work page 2006
-
[3]
Yazyev, O. V. Emergence of magnetism in graphene materials and nanostructures. Reports on Progress in Physics 2010, 73, 056501
2010
-
[4]
K.; Tretyakov, E.; Baumgarten, M.; Ardavan, A.; Sadeghi, H.; Lambert, C
Slota, M.; Keerthi, A.; Myers, W. K.; Tretyakov, E.; Baumgarten, M.; Ardavan, A.; Sadeghi, H.; Lambert, C. J.; Narita, A.; Müllen, K.; Bogani, L. Magnetic edge states and coherent manipulation of graphene nanoribbons. Nature 2018, 557, 691--695
work page 2018
-
[5]
de Oteyza, D. G.; Frederiksen, T. Carbon-based nanostructures as a versatile platform for tunable -magnetism. Journal of Physics: Condensed Matter 2022, 34, 443001
work page 2022
-
[6]
Pavliček, N.; Mistry, A.; Majzik, Z.; Moll, N.; Meyer, G.; Fox, D. J.; Gross, L. Synthesis and characterization of triangulene. Nature Nanotechnology 2017, 12, 308--311
work page 2017
-
[7]
A.; Müllen, K.; Fasel, R.; Feng, X.; Ruffieux, P
Mishra, S.; Beyer, D.; Eimre, K.; Liu, J.; Berger, R.; Gröning, O.; Pignedoli, C. A.; Müllen, K.; Fasel, R.; Feng, X.; Ruffieux, P. Synthesis and Characterization of -Extended Triangulene. Journal of the American Chemical Society 2019, 141, 10621--10625
work page 2019
-
[8]
Su, J. et al. Atomically precise bottom-up synthesis of -extended [5]triangulene. Science Advances 2019, 5, eaav7717
work page 2019
Show all 42 references
-
[9]
Triangulenes: From Precursor Design to On-Surface Synthesis and Characterization
Su, J.; Telychko, M.; Song, S.; Lu, J. Triangulenes: From Precursor Design to On-Surface Synthesis and Characterization. Angewandte Chemie International Edition 2020, 59, 7658--7668
2020
-
[10]
A.; Fasel, R.; Feng, X.; Ruffieux, P
Mishra, S.; Xu, K.; Eimre, K.; Komber, H.; Ma, J.; Pignedoli, C. A.; Fasel, R.; Feng, X.; Ruffieux, P. Synthesis and characterization of [7]triangulene. Nanoscale 2021, 13, 1624--1628
2021
-
[11]
On- Surface Synthesis and Characterization of [7] Triangulene Quantum Ring
Su, J.; Fan, W.; Mutombo, P.; Peng, X.; Song, S.; Ondráček, M.; Golub, P.; Brabec, J.; Veis, L.; Telychko, M.; Jelínek, P.; Wu, J.; Lu, J. On- Surface Synthesis and Characterization of [7] Triangulene Quantum Ring . Nano Letters 2021, 21, 861--867
2021
-
[12]
Observation of the Magnetic Ground State of the Two Smallest Triangular Nanographenes
Turco, E.; Bernhardt, A.; Krane, N.; Valenta, L.; Fasel, R.; Juríček, M.; Ruffieux, P. Observation of the Magnetic Ground State of the Two Smallest Triangular Nanographenes. JACS Au 2023, 3, 1358--1364
2023
-
[13]
Ovchinnikov, A. A. Multiplicity of the ground state of large alternant organic molecules with conjugated bonds. Theoretica chimica acta 1978, 47, 297--304
1978
-
[14]
Lieb, E. H. Two theorems on the Hubbard model. Phys. Rev. Lett. 1989, 62, 1201--1204
1989
-
[15]
A.; Fasel, R.; Feng, X.; Ruffieux, P
Mishra, S.; Beyer, D.; Eimre, K.; Ortiz, R.; Fernández-Rossier, J.; Berger, R.; Gröning, O.; Pignedoli, C. A.; Fasel, R.; Feng, X.; Ruffieux, P. Collective All - Carbon Magnetism in Triangulene Dimers **. Angewandte Chemie 2020, 132, 12139--12145
2020
-
[16]
L.; Franke, K
Krane, N.; Lotze, C.; Reecht, G.; Zhang, L.; Briseno, A. L.; Franke, K. J. High-Resolution Vibronic Spectra of Molecules on Molybdenum Disulfide Allow for Rotamer Identification. ACS Nano 2018, 12, 11698--11703
2018
-
[17]
Atomically Precise Synthesis and Characterization of Heptauthrene with Triplet Ground State
Su, X.; Li, C.; Du, Q.; Tao, K.; Wang, S.; Yu, P. Atomically Precise Synthesis and Characterization of Heptauthrene with Triplet Ground State . Nano Letters 2020, 20, 6859--6864
2020
-
[18]
Magnetic Excitations in Ferromagnetically Coupled Spin-1 Nanographenes
Turco, E.; Wu, F.; Catarina, G.; Krane, N.; Ma, J.; Fasel, R.; Feng, X.; Ruffieux, P. Magnetic Excitations in Ferromagnetically Coupled Spin-1 Nanographenes. 2024; https://arxiv.org/abs/2407.00728
2024 arXiv
-
[19]
On-surface synthesis of triangulene trimers via dehydration reaction
Cheng, S.; Xue, Z.; Li, C.; Liu, Y.; Xiang, L.; Ke, Y.; Yan, K.; Wang, S.; Yu, P. On-surface synthesis of triangulene trimers via dehydration reaction. Nature Communications 2022, 13, 1705
2022
-
[20]
Orbital-symmetry effects on magnetic exchange in open-shell nanographenes
Du, Q.; Su, X.; Liu, Y.; Jiang, Y.; Li, C.; Yan, K.; Ortiz, R.; Frederiksen, T.; Wang, S.; Yu, P. Orbital-symmetry effects on magnetic exchange in open-shell nanographenes. Nature Communications 2023, 14, 4802
2023
-
[21]
R.; Sanz, S.; Rey, D.; Corso, M.; Frederiksen, T.; Pascual, J
Hieulle, J.; Castro, S.; Friedrich, N.; Vegliante, A.; Lara, F. R.; Sanz, S.; Rey, D.; Corso, M.; Frederiksen, T.; Pascual, J. I.; Peña, D. On- Surface Synthesis and Collective Spin Excitations of a Triangulene - Based Nanostar . Angewandte Chemie International Edition 2021, 6...
2021
-
[22]
A.; Feng, X.; Ruffieux, P.; Fernández-Rossier, J.; Fasel, R
Mishra, S.; Catarina, G.; Wu, F.; Ortiz, R.; Jacob, D.; Eimre, K.; Ma, J.; Pignedoli, C. A.; Feng, X.; Ruffieux, P.; Fernández-Rossier, J.; Fasel, R. Observation of fractional edge excitations in nanographene spin chains. Nature 2021, 598, 287--292
2021
-
[23]
Zhao, C.; Catarina, G.; Zhang, J.-J.; Henriques, J. C. G.; Yang, L.; Ma, J.; Feng, X.; Gröning, O.; Ruffieux, P.; Fernández-Rossier, J.; Fasel, R. Tunable topological phases in nanographene-based spin-1/2 alternating-exchange Heisenberg chains. 2024; https://arxiv.org/abs/2402.13590
2024 arXiv
-
[24]
G.; Fischer, F
Delgado, A.; Dusold, C.; Jiang, J.; Cronin, A.; Louie, S. G.; Fischer, F. R. Evidence for excitonic insulator ground state in triangulene Kagome lattice. 2023; http://arxiv.org/abs/2301.06171
2023 arXiv
-
[25]
Catarina, G.; Henriques, J. C. G.; Molina-Sánchez, A.; Costa, A. T.; Fernández-Rossier, J. Broken-symmetry magnetic phases in two-dimensional triangulene crystals. Physical Review Research 2023, 5, 043226
2023
-
[26]
J.; Gupta, J
Heinrich, A. J.; Gupta, J. A.; Lutz, C. P.; Eigler, D. M. Single-Atom Spin-Flip Spectroscopy. Science 2004, 306, 466--469
2004
-
[27]
C.; Lambe, J
Jaklevic, R. C.; Lambe, J. Molecular Vibration Spectra by Electron Tunneling. Phys. Rev. Lett. 1966, 17, 1139--1140
1966
-
[28]
Lambe, J.; Jaklevic, R. C. Molecular Vibration Spectra by Inelastic Electron Tunneling. Phys. Rev. 1968, 165, 821--832
1968
-
[29]
C.; Rezaei, M
Stipe, B. C.; Rezaei, M. A.; Ho, W. Single-Molecule Vibrational Spectroscopy and Microscopy. Science 1998, 12, 1732--1735
1998
-
[30]
J.; Lutz, C
Heinrich, A. J.; Lutz, C. P.; Gupta, J. A.; Eigler, D. M. Molecule Cascades . Science 2002, 298, 1381--1387
2002
-
[31]
Single-molecule chemistry
Ho, W. Single-molecule chemistry. The Journal of Chemical Physics 2002, 117, 11033--11061
2002
-
[32]
Frisch, M. J. et al. Gaussian16 R evision C .01. 2016; Gaussian Inc. Wallingford CT
2016
-
[33]
C.; Krane, N.; Drost, R.; Sobrino, N.; Bernhardt, A.; Jur\' c c ek, M.; Fasel, R.; Ruffieux, P.; Liljeroth, P.; Jacob, D
Turco, E.; Aapro, M.; Ganguli, S. C.; Krane, N.; Drost, R.; Sobrino, N.; Bernhardt, A.; Jur\' c c ek, M.; Fasel, R.; Ruffieux, P.; Liljeroth, P.; Jacob, D. Demonstrating Kondo behavior by temperature-dependent scanning tunneling spectroscopy. Phys. Rev. Res. 2024, 6, L022061
2024
-
[34]
J.; Champagne, A
Parks, J. J.; Champagne, A. R.; Hutchison, G. R.; Flores-Torres, S.; Abru\ na, H. D.; Ralph, D. C. Tuning the Kondo Effect with a Mechanically Controllable Break Junction. Phys. Rev. Lett. 2007, 99, 026601
2007
-
[35]
J.; Pascual, J
Fern\'andez-Torrente, I.; Franke, K. J.; Pascual, J. I. Vibrational Kondo Effect in Pure Organic Charge-Transfer Assemblies. Phys. Rev. Lett. 2008, 101, 217203
2008
-
[36]
S.; Anders, F
Eickhoff, F.; Kolodzeiski, E.; Esat, T.; Fournier, N.; Wagner, C.; Deilmann, T.; Temirov, R.; Rohlfing, M.; Tautz, F. S.; Anders, F. B. Inelastic electron tunneling spectroscopy for probing strongly correlated many-body systems by scanning tunneling microscopy. Phys. Rev. B 20...
2020
-
[37]
Many-body effects in magnetic inelastic electron tunneling spectroscopy
Koryt\'ar, R.; Lorente, N.; Gauyacq, J.-P. Many-body effects in magnetic inelastic electron tunneling spectroscopy. Phys. Rev. B 2012, 85, 125434
2012
-
[38]
Spin excitations and correlations in scanning tunneling spectroscopy
Ternes, M. Spin excitations and correlations in scanning tunneling spectroscopy. New Journal of Physics 2015, 17, 063016
2015
-
[39]
Theory of Single Molecule Vibrational Spectroscopy and Microscopy
Lorente, N.; Persson, M. Theory of Single Molecule Vibrational Spectroscopy and Microscopy. Phys. Rev. Lett. 2000, 85, 2997--3000
2000
-
[40]
Symmetry Dependence of Vibration-Assisted Tunneling
Pavlicek, N.; Swart, I.; Niedenf\"uhr, J.; Meyer, G.; Repp, J. Symmetry Dependence of Vibration-Assisted Tunneling. Phys. Rev. Lett. 2013, 110, 136101
2013
-
[41]
L.; Franke, K
Reecht, G.; Krane, N.; Lotze, C.; Zhang, L.; Briseno, A. L.; Franke, K. J. Vibrational Excitation Mechanism in Tunneling Spectroscopy beyond the Franck-Condon Model. Phys. Rev. Lett. 2020, 124, 116804
2020
-
[42]
before" the step) were subtracted from the higher energy maps (taken
Bardeen, J. Tunnelling from a many-particle point of view. Phys. Rev. Lett. 1961, 6, 57 mcitethebibliography main.tex0000664000000000000000000005321014722334372011234 0ustar rootroot [ journal=jacsat, manuscript=article, ] achemso graphicx siunitx nicefrac dI/dV Vibrational ex...
1961
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.