REVIEW 4 major objections 5 minor
Non-equilibrium Effects in Vibrational Modes Pumped by Inelastic Tunneling
T0 review · 4 major / 5 minor · reviewed 2026-08-01 · deepseek-v4-flash
Pith's one-line read Inelastic tunneling can pump hydrogen vibrations on a delafossite surface into long-lived excited states lasting hundreds of picoseconds.
desk verdict A real current-dependent depopulation observation and independent DFT lifetimes are buried under a headline lifetime claim that rests on unmeasured φτ product and a self-contradictory linewidth passage. 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 γ1 in-plane H-Pd bending mode at 42 meV on the hydrogen tiling, and the two-level excitation model used to convert its current-dependent intensity into a lifetime. The model treats the mode as a two-level system with ground and first excited states, an excitation rate Γ=φI/e (quantum yield φ times electron current I), spontaneous relaxation rate 1/τ, and stimulated processes with ratio κ; the steady-state ground-state population is N_g=R/(R+1) with R=κ+(Γτ)^{-1}. Fitting the measured peak height h(I)=h0N_g yields the product φτ, so the absolute lifetime requires an assumption about φ. The paper also uses the observed quantum yield φ≈0.5 from the jump in conductance
What would settle it
Measure the quantum yield independently from the low-current conductance step height, where ground-state depletion is negligible, for each cluster type and check whether the resulting φ values are all near 0.5; alternatively, a two-pulse pump-probe STM measurement of the γ1 peak recovery would measure τ directly and would confirm or refute the hundreds-of-picosecond claim without relying on the two-level fit.
Extended reading notes
Core claim
On the Pd-terminated surface of PdCrO2, hydrogen forms a nonperiodic tiling of clusters, and the H-Pd bonds give sharp inelastic tunneling features. The paper identifies the 42 meV mode γ1 as an in-plane H-Pd vibration and shows that its intensity, normalized to tunneling current, falls monotonically as the current is raised from 50 pA upward. That current dependence is the non-equilibrium signature: the tunneling electrons excite the mode faster than it can relax, depopulating the ground state. Fitting this to a two-level rate-equation model yields values of φτ from 11.3 to 519 ps across five cluster types, together with stimulated de-excitation ratios κ between 0.04 and 0.69. Assuming a qu
Load-bearing premise
The absolute lifetime claim rests on the extra assumption that the quantum yield φ is about 0.5 and the same for every cluster, because the two-level fit fixes only the product φτ; if φ were smaller the inferred lifetimes would be even longer, and if larger, shorter.
Editorial extensions
If this is right
- If the inferred lifetimes are correct, STM-IETS can act as a pump as well as a probe: at currents above about 50 pA the vibrational ground state of a localized mode is measurably depleted, giving access to non-equilibrium vibrational populations on a single adsorbate cluster.
- Long lifetimes mean vibrational energy injected locally into H-Pd bonds is not rapidly dumped into the substrate, so the excitation remains available for subsequent processes such as overtone excitation or stimulated emission, consistent with the observed second harmonic at 84 meV.
- The contrast with H on Pd(111) suggests that surface-bulk electronic decoupling and low carrier density are the controlling factors, so similar long-lived adsorbate modes should appear on other delafossite and layered correlated materials.
- The spatial maps of the 42 meV in-plane and 274 meV out-of-plane modes closely resemble each other, indicating that the two vibrations are coupled, so driving one mode can affect the other; relevant to selective distortion driving of correlated materials.
- The paper links the long-lived vibrations to the high electrocatalytic activity of Pd delafossites in hydrogen evolution: long-lived H-Pd excitations could help overcome kinetic barriers in that reaction.
Reading between the lines
- The fitted quantity is always φτ, never τ alone; if the quantum yield is not constant across clusters, the spread in lifetimes may partly reflect variation in φ rather than variation in relaxation time. A direct per-cluster measurement of φ would separate the two.
- If the long lifetime is really caused by reduced electronic screening from magnetic correlations, then suppressing the antiferromagnetic order—by heating above the ordering temperature or applying a magnetic field—should shorten the vibrational lifetime; this is a testable prediction the paper does not state explicitly.
- The paper's own n-level extension shows the lifetime estimate is stable to adding more levels, but the stimulated-emission ratio κ grows with n, suggesting κ should be read as an effective parameter rather than a microscopic rate.
- The same pump-and-probe scheme could be applied to other strongly correlated surfaces with adsorbate clusters, or used to search for modes whose lifetimes are long enough for coherent control with tailored voltage pulses.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports STM-based inelastic electron tunneling spectroscopy (IETS) on hydrogen clusters on the Pd-terminated PdCrO2 surface. It identifies several vibrational features, most prominently the 42 meV γ1 H–Pd mode, and observes that its normalized IETS peak height decreases monotonically with increasing tunneling current. This decrease is interpreted as depopulation of the vibrational ground state under electron pumping, and is fitted with a two-level model h(I)=h0 N_g with R=κ+(Γτ)^{-1}, Γ=φI/e. The fits yield φτ products between 11.3 ps and 519 ps for different cluster types, together with κ values; assuming φ≈0.5 for all clusters, the authors infer lifetimes up to hundreds of picoseconds. Density-functional-theory calculations of vibrational lifetimes on PdCrO2 and Pd(111) are presented as an independent comparison. The paper's central claim is that these modes have longer lifetimes than on typical metal surfaces, attributed to electronic decoupling and correlation effects in the delafossite substrate.
Significance. If the inferred lifetimes are correct, this would be a notable demonstration of non-equilibrium vibrational pumping in a strongly correlated surface system, with implications for understanding energy dissipation and for hydrogen-evolution catalysis on delafossite oxides. The paper's strengths are its raw experimental observation of a current-dependent decrease in the normalized peak height—a model-independent trend—and the independent DFT lifetime comparison, which provides a plausibility check. However, the absolute lifetime claim is not directly measured: the fit determines the product φτ, and the conversion to τ rests on the unstated accuracy and constancy of φ. The manuscript also contains a direct internal contradiction concerning the linewidth-derived lifetime. These issues affect the central quantitative claim and require substantial revision.
major comments (4)
- [Fig. 3(c)–(d), two-level model (h(I)=h0N_g)] The fit determines the product φτ, not τ. The abstract and conclusion state 'measured lifetimes reaching up to hundreds of picoseconds', but this is only obtained after assuming φ≈0.5 and that φ is identical for all cluster types. The text reports a factor-of-two conductance jump for one T3 cluster, and 'similar step heights' elsewhere, but no per-cluster φ values, no confidence intervals on φτ, and no goodness-of-fit measure. Because the comparison to metal-surface lifetimes depends entirely on this conversion, the manuscript should provide joint confidence regions for (φτ, κ, h0) and a sensitivity analysis for τ versus φ, or explicitly reframe the claim as a φτ product.
- [Paragraphs around Fig. 3(a), linewidth lifetime] The text states first that 'the lifetime of the γ1 mode determined from its apparent energy width (~0.3 ps) is of a similar order...' and immediately afterwards that 'a determination of the lifetime of the vibrational state from its linewidth is not possible, because the width is comparable to the intrinsic resolution limit given by thermal broadening.' These statements are contradictory. If the linewidth is resolution-limited, the 0.3 ps value is not a measured lifetime; if it is meaningful, it would point to a short lifetime inconsistent with the long-lifetime claim. This needs to be resolved, since it bears directly on the validity of the central result.
- [Fig. 2(b), γ2 overtone at 84 meV] The observed γ2 overtone is assigned to progressive pumping but is never used quantitatively. Since the two-level model predicts a specific current dependence of the ground-state population, the current dependence of γ2, or of the γ2/γ1 ratio, should provide an independent cross-check on Γτ. Performing this analysis—or explaining why it cannot be done—would strengthen, or possibly falsify, the proposed pumping model. Without it, the lifetime inference rests on a single fitted product.
- [Fit robustness (Fig. 3(c), T5 cluster; n-level model)] The T5 cluster could not be fitted, and the n-level robustness check is referred to only in the Supplemental Material. The central claim is based on five clusters, one of which fails, and the reported κ values appear to drift with n. Please report the number of data points per cluster, the fit uncertainties, and the n-level results in the main text or provide a reason for excluding T5. As written, the reader cannot assess whether the inferred hierarchy of lifetimes among clusters is statistically meaningful.
minor comments (5)
- [General] Several unit errors: α at '3 mV' and the text 'os comparable' should read 'is comparable'; energies should be meV, not mV. The sentence 'Back to our discussion of the results based on the two-level model, Assuming...' is a fragment and should be rewritten.
- [Fig. 2 callouts] The text says a topographic region comprising ~50 clusters is shown in [Fig. 2(a)], but Fig. 2(a) is a g(V) spectrum; the topographic image is Fig. 2(c). Please correct all figure callouts.
- [Supplemental Fig. S9] Supplemental Fig. S9 is cited both for the central-atom protrusion evidence and for calculated lifetimes; this appears to be two different figures. Please renumber and disambiguate.
- [Quantum yield section] The relationship between the observed factor-of-two increase in g and φ≈0.5 should be stated explicitly. In standard IETS, the step height is not simply twice φ; the definition of φ used here should be given, otherwise the assumption of φ=0.5 is not transparent.
- [Data availability] The data availability statement contains a placeholder '[?]'. This should be completed before publication.
Circularity Check
No circular derivation: the current-dependent peak-height fit is a transparent inversion for the stated φτ product, and the long-lifetime conclusion is independently anchored by DFT; the linewidth contradiction is a correctness issue rather than circularity.
full rationale
The paper's central quantitative claim is not a prediction that is equivalent by construction to its input. The observed monotonic decrease of the normalized γ1 peak height with current (Fig. 3(c)) is a measurement; the two-level fit h(I)=h0·N_g with R=κ+(Γτ)^{-1} and Γ=φI/e is an inversion. The fit output is explicitly reported as 'estimated lifetimes (expressed in ϕτ)' (11.3–519 ps), and the conversion to τ is made under the stated assumption that φ≈0.5 is the same for all cluster types. That is a model-dependent inference, not a definitional identity: the same data plus a different φ would give a different absolute τ, but the paper does not hide the product or the assumption. Independent DFT lifetime calculations for H on Pd/PdCrO2 versus Pd(111) (Fig. 4, 'our results ... show that on the Pd-terminated surface ... lifetimes of the bending modes are up to nine times longer than those of hydrogen on Pd(111)') provide an external benchmark for the qualitative 'longer than typical metals' conclusion. Mode assignments lean on the authors' prior Ref. [19], but the present DFT and external phonon calculations [28] anchor them independently, so this is not a circular self-citation chain. I flag two non-circular defects: (1) the manuscript states both 'The lifetime of the γ1 mode determined from its apparent energy width (~0.3 ps)' and 'A determination of the lifetime ... from its linewidth is not possible'; these conflict and undermine confidence in any single lifetime value, and (2) the observed γ2 overtone (84 meV) is not used quantitatively to cross-check Γτ. Neither reduces the derivation to its inputs by construction.
Assumptions & free parameters
free parameters (3)
- φτ (quantum-yield-times-lifetime) per cluster =
11.3–519 ps for T1, T3, T4, T7; T5 not fitted
- κ (stimulated de-excitation / excitation ratio) =
0.04–0.69 for fitted clusters
- h0 proportionality factor =
not reported individually
assumptions (5)
- domain assumption Two-level excitation model: transitions only between ground and first excited state, and only when the ground state is populated.
- ad hoc to paper The quantum yield φ is the same for all cluster types and independent of current and bias.
- domain assumption Normalized peak height h_norm is proportional to the time-average ground-state population N_g.
- domain assumption The 42 meV mode γ1 is the in-plane H-Pd vibration, 84 meV is its overtone, and the 3/15/59/274 meV features are assigned as stated.
- domain assumption DFT (PBE/GGA) with and without antiferromagnetic ordering of the CrO2 layer captures the relevant vibrational lifetime physics.
Cite this review
Pith. "Pith review of Non-equilibrium Effects in Vibrational Modes Pumped by Inelastic Tunneling." pith.science (2026). https://pith.science/paper/KP7B3BYE
@misc{pith2026260717530,
author = {Pith},
title = {Pith review of: Non-equilibrium Effects in Vibrational Modes Pumped by Inelastic Tunneling},
year = {2026},
howpublished = {\url{https://pith.science/paper/KP7B3BYE}},
note = {Machine review of arXiv:2607.17530}
}
read the original abstract
The properties of strongly correlated electron materials exhibit a surprising sensitivity to small lattice distortions, providing an opportunity for their tuning by selective distortion driving, usually achieved by optical excitations. Using inelastic electron tunneling in scanning tunneling microscopy, we demonstrate that at the surface of a strongly correlated electron material, we can drive vibrational excitations out of equilibrium, by studying the dynamics of localized modes on the Pd-terminated surface of the delafossite oxide PdCrO2. This surface forms a tiling of hydrogen clusters of varying sizes and shapes upon hydrogen adsorption. Our findings reveal that vibrational excitations in the clusters exhibit longer lifetimes than on typical metal surfaces. Detailed analysis of the spectroscopy data reveals signatures of non-equilibrium effects in the excitations which we attribute to the extended lifetimes of these modes. Theoretical calculations support that the long-lived nature of the excitations is related to the unique properties of the substrate.
Figures
Figures from the paper (1 more)
Reviewed August 1, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.