REVIEW 1 major objections 5 minor 4 references
Tip-induced nitrene generation
T0 review · 1 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read This paper reports that voltage pulses can strip all three azide groups from 2,5,8-triazido-s-heptazine, producing trinitreno-s-heptazine; broken-symmetry DFT and configuration-interaction calculations predict a septet neutral and a…
desk verdict Solid tip-chemistry result with a robust gas-phase high-spin prediction; the on-surface sextet assignment is a theory-based inference, not a measurement. 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 key object is trinitreno-s-heptazine, a heptazine core with three nitrene nitrogen atoms, each carrying an $S=1$ triplet of two unpaired electrons. The reaction pathway is driven by the scanning tip: voltage pulses of 2–2.6 V with tunnelling currents up to about 20 pA detach one $N_2$ molecule per azide group, and the authors find that successive dissociations require slightly stronger pulses. The structural identification uses atomic force microscopy with a CO-functionalized tip, interpreted with the probe-particle model, which translates calculated adsorption geometries into simulated AFM images. The spin-state argument is carried by broken-symmetry DFT and difference-dedicated configuration interaction (DDCI) calculations whose energy differences are mapped onto a Heisenberg–Dirac–van Vleck spin Hamiltonian with an intra-nitrene coupling $J'$ and an inter-nitrene coupling $J$; all fitted couplings are positive. On-surface, PBE DFT geometry optimizations of the anion and dianion in all spin multiplicities, combined with probe-particle AFM simulations, select the anionic sextet as the state that matches experiment and has the lowest computed energy.
What would settle it
Perform a spin-sensitive single-molecule measurement on 4, such as inelastic electron tunnelling spectroscopy in a magnetic field: a sextet ($S=5/2$) anion should show a ladder of spin excitations, while a doublet ($S=1/2$) would show at most a single excitation. Observing the doublet signature, or no spin excitation, would show the molecule is not the anionic sextet.
Extended reading notes
Core claim
The central claim is that tip-induced dissociation of all three azide groups of 2,5,8-triazido-s-heptazine (1) on bilayer NaCl/Au(111) produces trinitreno-s-heptazine (4), and that this molecule is a high-spin organic species: a septet ($S=3$) in the gas phase as a neutral and a sextet ($S=5/2$) as an anion on the surface. The evidence combines atomically resolved AFM images that show the three nitrene nitrogen atoms, an STM dark halo and interface-state scattering that indicate charging of 4, Kelvin-probe force spectra consistent with a negative charge, and quantum-chemical calculations. The computed adiabatic electron affinity of 5.29 eV together with a reduced work function of the NaCl-covered Au(111) surface supports the anionic state. Broken-symmetry DFT and DDCI give positive intra-nitrene couplings $J'$ of about 0.4–1.1 eV and positive inter-nitrene couplings $J$ of about 1–6 meV, implying ferromagnetic alignment of all six unpaired electrons. On the surface, the calculated electron addition puts the extra electron in the heptazine core with opposite spin, reducing the total spin from 3 to 5/2; AFM probe-particle simulations of the relaxed anion match the measured contrast, and the calculations place the sextet below the other spin states.
Load-bearing premise
The conclusion that the observed molecule is the anionic sextet depends on the approximate quantum-chemistry calculation correctly ordering the spin-state energies on the salt film, since the microscope images match both the doublet and the sextet equally well.
Editorial extensions
If this is right
- All three azide groups of one precursor can be dissociated sequentially by voltage pulses, so the same molecule can be converted from triazide to mono-, di-, and trinitrene in place.
- The calculated exchange couplings imply neutral trinitreno-s-heptazine is a septet: three ferromagnetically coupled $S=1$ nitrene centers forming a total spin $S=3$.
- On bilayer NaCl/Au(111), the trinitrene is an anion within at least $\pm0.7$ V of bias, with the extra electron antiferromagnetically coupled to the nitrene spins to give a sextet ($S=5/2$).
- Because the nitrene spin density stays localized on the three nitrogen atoms, the molecule is a candidate building block for on-surface arrays of well-separated spins.
- The strong intra-nitrene coupling compared with the weak inter-nitrene coupling means the local triplet on each nitrene is robust, while ordering between centers is weak and should be sensitive to the connecting heptazine core.
Reading between the lines
- If the electron affinity versus work-function logic is correct, depositing the same molecule on a higher-work-function substrate should leave it neutral; a neutral septet on such a surface would confirm the charge-state switching independently of the spin-state calculation.
- The order-of-magnitude gap between the intra-nitrene coupling (~0.4–1.1 eV) and the inter-nitrene coupling (~1–6 meV) suggests that nitrene-based spin lattices would behave as weakly coupled local triplets at accessible temperatures; long-range order would require stronger linker-mediated exchange than the bare heptazine core provides.
- Because the probe-particle AFM simulations reproduce the heptazine core but not the nitrene atoms, a full charge-density AFM simulation could test whether the calculated nitrene–surface distance is too small; that would directly probe the adsorption geometry that enters the on-surface spin-state assignment.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper reports the generation of mono-, di-, and trinitreno-s-heptazine (2, 3, 4) from 2,5,8-triazido-s-heptazine (1) by voltage pulses applied with an STM/AFM tip on bilayer NaCl/Au(111). AFM and STM are used to identify the precursor and the successive loss of N2 from the azide groups. Gas-phase broken-symmetry DFT (B3LYP and M06-2X) and DDCI calculations are mapped onto Heisenberg spin Hamiltonians, yielding positive intra-nitrene (J') and inter-nitrene (J) exchange couplings; the authors conclude that neutral 4 has a septet ground state with three ferromagnetically coupled S=1 nitrene centers. On the surface, a dark halo, standing-wave pattern, and Kelvin probe force spectroscopy indicate that 4 carries a negative charge. PBE-D3 calculations on a NaCl slab combined with probe-particle AFM simulations lead the authors to assign 4 as an anionic sextet.
Significance. The potential significance is high if the claims hold: the paper would demonstrate on-surface generation of a molecule with three localized spin centers and a high-spin ground state, and it offers a concrete route to nitrene-based molecular building blocks. The work has genuine strengths: the AFM/STM identification of the precursor and its three products is convincing and carefully documented; the gas-phase spin-state prediction is not circular, as the exchange couplings are obtained from ab initio energies and their positive sign is consistent across BS-DFT and DDCI; and the calculated J' values agree with independent experimental aryl nitrene couplings. The DDCI calculation covers the full low-energy spectrum for 4, and the predicted neutral septet is a concrete falsifiable statement that could be tested by matrix ESR. The main weakness is that the experimental assignment of the surface species as an anionic sextet is not pinned by the imaging data; as the authors acknowledge, the probe-particle AFM simulations match the doublet and sextet equally well, so the sextet assignment rests on a PBE-D3 energy ordering that is not benchmarked.
major comments (1)
- [Results and discussion, paragraph 'To compare the AFM measurements of 4 on NaCl with simulations'; Fig. S4] The assignment of 4 on the surface as an anionic sextet is not experimentally pinned. The probe-particle AFM simulations match the experimental contrast for both the anionic doublet and the anionic sextet, so the imaging data cannot discriminate between these spin states. The only evidence favoring the sextet is the PBE-D3 relative energy ordering of spin states on the NaCl slab, but no numerical gaps are reported and no benchmark or error estimate is provided. PBE is known to have uncontrolled errors for spin-state splittings of open-shell organic radicals, so this ordering needs a direct check. A second uncontrolled factor is that the on-surface DFT model in Methods uses a freestanding NaCl(100) slab without the Au(111) substrate, whereas the experiment is performed on bilayer NaCl on Au(111); for an anionic adsorbate the metal will influence screening and spin-state energetics. Because the central claim of a high-spin molecule on the surface rests on this ordering, the authors should either provide a more trustworthy estimate (for example, hybrid or wavefunction-based calculations on a cluster model, or at least report the gaps and validate PBE against a known nitrene system) or explicitly state that the surface spin state is not experimentally determined.
minor comments (5)
- [Methods: Electron affinity calculations] The CASPT2 calculation is described as an adiabatic electron affinity, but the anion is evaluated at the neutral optimized geometry without relaxing the anion. Please either relax the anionic geometry or describe the value as a vertical electron affinity at the neutral geometry.
- [Results and discussion, on-surface DFT paragraph] The statement that 'the DFT calculated energies indicate that the sextet is the anionic ground state' should be accompanied by the actual relative energies of the doublet, quartet, and sextet states; without numerical gaps the reader cannot assess the margin or the robustness of the ordering.
- [Fig. S4 caption] The caption says that energies of the spin multiplicities are indicated, but the values are not visible in the text provided here; please ensure the numerical labels are legible and clearly define the zero of energy.
- [Results and discussion, charge-state argument] The assumed ~1 eV work-function decrease for bilayer NaCl on Au(111) is taken from the literature; since the dark halo, the standing-wave pattern, and the Kelvin probe data directly indicate a negative charge state, please state explicitly that the electron-affinity/work-function comparison is only a consistency check rather than the primary evidence for charging.
- [Notation throughout] The notations '740' and '64–' should be defined at first use (for example, superscript spin multiplicity and subscript charge) to avoid confusion with ordinary chemical formulas.
Circularity Check
No significant circularity: the spin-state predictions are ab initio results benchmarked against independent experimental values, not fits to the measured outcome.
full rationale
The paper's derivation chain is self-contained against its inputs. The central gas-phase prediction of a septet ground state for neutral 4 follows from broken-symmetry DFT (B3LYP and M06-2X) and DDCI calculations mapped onto Heisenberg spin Hamiltonians; the resulting positive intra- and inter-nitrene exchange couplings J' and J are ab initio outputs, not parameters fitted to the observed spin state. The calculated J' values are explicitly compared with independent experimental aryl-nitrene and NH values, providing external validation. The on-surface assignment of 4 as an anionic sextet rests on two distinct lines: experimental charge-state evidence (STM dark halo, standing-wave scattering, KPFS) supported by an external CASPT2 electron affinity and literature work-function values, and DFT relative energies among spin states on the NaCl bilayer. Although the AFM probe-particle simulations match both the doublet and sextet equally well, that fact makes the spin assignment dependent on the accuracy of the PBE spin-state ordering, which is a benchmarking/accuracy concern rather than circularity: the PBE energies are not fitted to the experimental spin state, and the imaging data do not enter the energy calculation. Self-citations appear in methodological contexts (e.g., refs 37, 40, 71), but the load-bearing claims are independently supported by external references or by direct experimental data, so no step reduces by construction to its own inputs.
Assumptions & free parameters
free parameters (3)
- Work function decrease of bilayer NaCl on Au(111) relative to Au(111) =
approx. 1 eV (assumed)
- Exchange coupling parameters J and B from spin-Hamiltonian fits =
J = 1.49 to 6.56 meV, B = 0.28 to 0.87 meV depending on method
- CO-tip probe-particle parameters =
lateral stiffness 0.25 N/m, radial stiffness 30 N/m, O partial charge -0.1 e
assumptions (4)
- domain assumption Each nitrene center carries two unpaired electrons in a triplet (S=1) configuration
- domain assumption The relaxed gas-phase trinitrene has threefold symmetry, so all inter-nitrene couplings J are equal and all intra-nitrene couplings J' are equal
- domain assumption The DDCI and broken-symmetry DFT energy spectra can be faithfully represented by the Heisenberg-Dirac-van Vleck bilinear and biquadratic spin Hamiltonian
- domain assumption PBE DFT relative energies of different spin multiplicities for the anion on the NaCl surface are accurate enough to identify the ground state
Cite this review
Pith. "Pith review of Tip-induced nitrene generation." pith.science (2026). https://pith.science/paper/DBV62B2E
@misc{pith2026250604741,
author = {Pith},
title = {Pith review of: Tip-induced nitrene generation},
year = {2026},
howpublished = {\url{https://pith.science/paper/DBV62B2E}},
note = {Machine review of arXiv:2506.04741}
}
read the original abstract
We generated trinitreno-s-heptazine, a small molecule featuring three nitrene centers, by tip-induced chemistry from the precursor 2,5,8-triazido-s-heptazine on bilayer NaCl on Au(111). The precursor's azide groups were dissociated to form mono-, di- and trinitreno-s-heptazine, yielding molecules with one to three nitrene centers. The precursor and its products are characterized by atomic force microscopy and scanning tunnelling microscopy. Broken-symmetry DFT and configuration interaction calculations of inter- and intra-nitrene exchange couplings suggest a ferromagnetic coupling of the S = 1 nitrene centers, resulting in a high-spin septet ground state for neutral trinitreno-s-heptazine in the gas phase. On bilayer NaCl on Au(111), the combined results of experiments and theory suggest trinitreno-s-heptazine to be an anion with a sextet ground state.
Figures
Figures from the paper (2 more)
Reference graph
Works this paper leans on
-
[97]
(6) Huang, K.; Leung, L.; Lim, T.; Ning, Z.; Polanyi, J
https://doi.org/10.1021/acs.nanolett.5b02974. (6) Huang, K.; Leung, L.; Lim, T.; Ning, Z.; Polanyi, J. C. Single-Electron Induces Double-Reaction by Charge Delocalization. J. Am. Chem. Soc. 2013, 135 (16), 6220–6225. https://doi.org/10.1021/ja400612c. (7) Mohn, F.; Repp, J.; Gross, L.; Meyer, G.; Dyer, M. S.; Persson, M. Reversible Bond Formation in a Gol...
-
[347]
(31) Bastardis, R.; Guihéry, N.; de Graaf, C
https://doi.org/10.1016/S0010-8545(00)80132-0. (31) Bastardis, R.; Guihéry, N.; de Graaf, C. Microscopic Origin of Isotropic Non-Heisenberg Behavior in $S=1$ Magnetic Systems. Phys. Rev. B 2007, 76 (13), 132412. https://doi.org/10.1103/PhysRevB.76.132412. (32) Brown, H. A. A Three-Spin Interaction Model. Journal of Magnetism and Magnetic Materials 1984, 4...
-
[7295]
https://doi.org/10.1021/acsnano.5b02147. (5) Schendel, V.; Borca, B.; Pentegov, I.; Michnowicz, T.; Kraft, U.; Klauk, H.; Wahl, P.; Schlickum, U.; Kern, K. Remotely Controlled Isomer Selective Molecular Switching. Nano Lett. 2016, 16 (1), 93–
-
[7423]
(36) Li, X.-Y.; Lou, F.; Gong, X.-G.; Xiang, H
https://doi.org/10.1021/acs.nanolett.4c01604. (36) Li, X.-Y.; Lou, F.; Gong, X.-G.; Xiang, H. Constructing Realistic Effective Spin Hamiltonians with Machine Learning Approaches. New J. Phys. 2020, 22 (5), 053036. https://doi.org/10.1088/1367- 2630/ab85df. (37) Murrell, J. N.; and Teixeira-Dias, J. J. C. The Dependence of Exchange Energy on Orbital Overla...
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.