Pith. sign in

REVIEW 5 minor 2 cited by

The first-passage time of an underdamped harmonic oscillator splits into three pieces fixed by Q, energy diffusion, and Kramers eigenvalues.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

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

T0 review · grok-4.5

2026-07-12 08:50 UTC pith:OWNDATQ3

load-bearing objection Solid, usable FPT density for underdamped oscillators: three explicit pieces, spectral and energy-diffusion rates, clean numerics; one labeled short-time ansatz.

arxiv 2607.01405 v2 pith:OWNDATQ3 submitted 2026-07-01 cond-mat.stat-mech

First passage time distribution in underdamped harmonic oscillators

Aubin Archambault , Caroline Crauste-Thibierge , Alberto Imparato , Sergio Ciliberto , Ludovic Bellon This is my paper
classification cond-mat.stat-mech
keywords first passage timeunderdamped oscillatorquality factorKramers equationenergy diffusionabsorbing boundaryinstantonnoise pattern
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

This paper gives a closed-form distribution for the first time an underdamped harmonic oscillator, started from equilibrium, reaches a fixed position threshold. Because the oscillator has inertia, the answer depends strongly on the quality factor Q and cannot be taken from the well-studied overdamped case. The authors split the density into an instantaneous piece (initial conditions already above threshold), a short-time piece controlled by pure Hamiltonian rotation of high-energy states, and a long-time piece whose exponential rate is the slowest eigenvalue of the Kramers operator with an absorbing wall (or, for very large Q, the mean energy-crossing time plus half a period). Direct Langevin simulations confirm the formulas across a wide range of Q and barrier heights. They also show that every trajectory that hits the threshold is driven by a characteristic, resonantly amplified noise pattern whose average shape is known exactly from time-reversal of the free relaxation.

Core claim

The first-passage-time density of an underdamped harmonic oscillator is exactly the sum of three contributions: an equilibrium Dirac mass for initial conditions already above threshold, a Hamiltonian flux of high-energy states that empties within one oscillation period, and a long-time survival probability whose rate is the slowest eigenvalue of the absorbing Kramers operator (or the energy-diffusion mean time plus half a period when Q is very large). These expressions match numerical Langevin trajectories for all tested Q and barrier heights.

What carries the argument

The decomposition P(t_fp)=P_I δ(t_fp)+P_II(t_fp)+P_III(t_fp), with P_III built from the slowest eigenvalue λ_1 of the Kramers operator L_xB that includes an absorbing sink at x=x_B (or, equivalently, from the energy-diffusion mean first-passage time τ_Z for large Q).

Load-bearing premise

For intermediate times shorter than one relaxation period the authors replace the true time-dependent escape rate by the long-time eigenvalue evaluated at the moving lowest-contact energy; if that educated guess fails, intermediate plateaus are no longer predicted quantitatively.

What would settle it

A long Langevin trajectory (or experiment) at moderate Q whose first-passage histogram fails to match the three-piece formula, especially the height and duration of the early plateaus controlled by the instantaneous contact energy E†(t).

Watch this falsifier — get emailed when new claim-graph text bears on it.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

0 major / 5 minor

Summary. The manuscript derives the first-passage-time distribution P(t_fp) for the position of an underdamped harmonic oscillator to reach a threshold x_B. The density is decomposed into three contributions (Eqs. 2–3): an instantaneous Dirac piece P_I from initial conditions already above threshold, a short-time Hamiltonian flux P_II of high-energy initial conditions (t_fp ≤ 2π), and a long-time Kramers-escape piece P_III. For moderate Q the long-time rate is identified with the slowest eigenvalue λ_1 of the absorbing Kramers operator L_xB (Sec. III); for large Q it is obtained from the mean first-passage time of the energy-diffusion process plus a half-period delay (Sec. IV, Eqs. 18–26). Short-time dynamics are treated in the Hamiltonian approximation. Direct Langevin simulations, numerical diagonalization of L_xB, and the closed-form energy-diffusion time all agree. The mean trajectories conditioned on first passage are shown to follow an instanton/relaxation path driven by a characteristic noise pattern (Sec. V).

Significance. A quantitative FPT distribution for underdamped oscillators has been missing; most prior analytic work addresses the overdamped limit. The paper supplies closed-form expressions that cover the full range of Q and t_fp, validated by independent numerical methods (eigenvalue spectra of L_xB, long Langevin trajectories, energy-diffusion MFPT) with no free parameters. The companion Letter already uses the mean FPT for an information-engine power estimate that matches experiment; the present article supplies the technical foundation and extends the result to very large Q. The identification of the noise pattern that generates escape events is an additional, potentially reusable insight. These results are of direct interest for stochastic thermodynamics, underdamped information engines, and any setting in which inertial first-passage statistics matter.

minor comments (5)
  1. Sec. VI: the short-time rate ansatz Γ(t)=λ_1(E†(t)) for π<t<2π is clearly labeled an educated guess and works well for Q=7, but a one-sentence quantitative statement of residual discrepancy (or a brief comparison against the full spectrum of L_xB) would help readers judge its accuracy for intermediate Q.
  2. Fig. 7 and the accompanying text note the plateau structure of the transient for large Q; a short remark that these plateaux arise from the imaginary parts of the higher eigenvalues of L_xB (already mentioned in Sec. VI) would make the connection fully explicit.
  3. Eq. (5) and the surrounding paragraph: the rough estimate Γ_K=e^{-B}/Q is useful, but a brief cross-reference to the more accurate expressions Γ_B=λ_1 and Γ_Z later in the paper would avoid any impression that Γ_K is the final result.
  4. Appendix A: the reconstruction of the FPT histogram from residence times is elegant; a single sentence clarifying that the same long trajectory can be re-analyzed for any x_B would further emphasize the numerical efficiency.
  5. Minor typographical points: “tfp” vs “t_fp” consistency, and a few missing spaces after commas in the abstract and Sec. II.

Circularity Check

0 steps flagged

No significant circularity: rates from independent operator eigenvalues or energy-diffusion MFPT; simulations are external validation only.

full rationale

The claimed FPT density decomposes into three explicit pieces (Eqs. 2–3). P_I and P_II follow from the equilibrium measure and Hamiltonian phase-space rotation (Fig. 1), independent of any rate. For long times the rate Γ is obtained either as the slowest eigenvalue λ_1 of the absorbing Kramers operator L_xB (numerically diagonalized by FEM, Sec. III) or, at large Q, as the closed-form energy-diffusion mean first-passage time τ_Z (Eqs. 18–20) plus a half-period delay (Eq. 26). Both constructions are parameter-free and do not use the Langevin histograms that later serve as validation (Figs. 2–3, 7). The only uncontrolled step is the short-time ansatz Γ(t)=λ_1(E†(t)) for π<t<2π (Sec. VI), which is explicitly labeled an “educated guess” and does not feed back into the long-time exponential or the large-Q formula. Self-citations to the companion Letter [17] merely point to the intermediate-Q summary; the present manuscript supplies the detailed derivations. No fitted parameter is renamed a prediction, no uniqueness theorem is imported, and no definitional identity is disguised as a result. The derivation chain is therefore self-contained against external numerical benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 4 axioms · 0 invented entities

The central claim rests only on the standard underdamped Langevin equation, the associated Kramers operator, and two well-known asymptotic reductions (Hamiltonian rotation at short times, energy diffusion at large Q). No numerical constants are fitted to produce the distribution; all free parameters of the model (Q, B) are inputs that are scanned, not adjusted.

axioms (4)
  • domain assumption Position and velocity obey the underdamped Langevin equations (1a,b) with additive white noise of intensity fixed by the fluctuation-dissipation relation.
    Stated at the opening of Sec. II; all subsequent operators and simulations start from this SDE.
  • domain assumption For Q ≫ 1 the energy is a slow variable whose marginal density satisfies the energy-diffusion Fokker–Planck equation with drift A = −E/Q and diffusion D = E/Q.
    Invoked in Sec. IV; standard reduction of Kramers’ equation for weak friction.
  • standard math The long-time survival probability inside the absorbing domain is dominated by the eigenvalue of largest real part of the non-self-adjoint Kramers operator with absorbing boundary.
    Used in Sec. III; follows from the spectral expansion of a linear evolution operator with negative spectrum.
  • domain assumption At short times (t ≲ 2π) the phase-space density rotates rigidly under the Hamiltonian flow before diffusion becomes appreciable.
    Basis of the P_I and P_II contributions (Sec. II and Fig. 1).

pith-pipeline@v1.1.0-grok45 · 17451 in / 2448 out tokens · 34113 ms · 2026-07-12T08:50:50.391292+00:00 · methodology

0 comments
read the original abstract

We derive the distribution of the first passage time $t_{fp}$ for the position $x$ of an underdamped harmonic oscillator to overcome a threshold $x_B$. As the $t_{fp}$ distribution depends on the oscillator quality factor $Q$ different approaches are used. At very large quality factor ($Q\gg 100$) and intermediate and long $t_{fp}$ the proof is based on an energy diffusion model, whereas at medium quality factor ($Q\sim 10$) the proof is based on the study of the eigenvalues of the Kramers linear differential operator with absorbing boundary conditions. For all $Q$ and short $t_{fp}$ we use a Hamiltonian approximation. The theoretical predictions are in excellent agreement with direct numerical simulations of underdamped oscillator dynamics. Finally we show that the mean of the trajectories ending at $t_{fp}$ presents a particular shape driven by a specific noise pattern.

Figures

Figures reproduced from arXiv: 2607.01405 by Alberto Imparato, Aubin Archambault, Caroline Crauste-Thibierge, Ludovic Bellon, Sergio Ciliberto.

Figure 1
Figure 1. Figure 1: FIG. 1. Phase space ( [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. (a) Eigenvalues [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. (a) Phase space density P(x, v, t) at large times ( [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Distribution of phase [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. (a) Mean trajectories of the position [PITH_FULL_IMAGE:figures/full_fig_p006_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Distribution of first passage time [PITH_FULL_IMAGE:figures/full_fig_p007_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: We’re only left with normalizing the histogram to estimate the pdf P(tfp), with a very high precision: we typically run simulation long enough to have N ∼ 106 . The benefit of this approach is also that we can compute P(tfp) for any value of xB on the same dataset, without running the simulation again. Appendix B: Distributions P(x) and P(v) at E = B In section Sec. IV we have shown that the pdf of the ang… view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. The pdfs of [PITH_FULL_IMAGE:figures/full_fig_p008_9.png] view at source ↗

Review history (2 revisions) →

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. First passage time for an underdamped harmonic oscillator and application to the power of an information engine

    cond-mat.stat-mech 2026-07 conditional novelty 7.0

    Closed-form first-passage-time distribution for an underdamped harmonic oscillator (short-time Hamiltonian + long-time Kramers) agrees with micro-cantilever data and yields the power of an information engine.

  2. First passage time for an underdamped harmonic oscillator and application to the power of an information engine

    cond-mat.stat-mech 2026-07 unverdicted novelty 5.0

    Derives first passage time distribution for underdamped harmonic oscillator via eigenvalues and approximations, applies to information engine power with micro-cantilever experiment.

Reference graph

Works this paper leans on

29 extracted references · 3 linked inside Pith · cited by 1 Pith paper

  1. [2]

    The mean response of the linear system thus corre- sponds to the relaxation of an harmonic oscillator start- ing inx(t fp) =x B with velocity⟨v(t fp)⟩. In Fig. 6(b), we plot the distribution of velocities att fp, which is inde- pendent of the quality factorQand perfectly described byP tfp(v) =ve −v2/2. Such expression is expected by noticing thatJ x =v P(...

  2. [3]

    t1 "t2 "t3

    BesidesP I andP II, dealing with intrinsically short time contributions, an educated guess is used for contri- butionP III described by Eq. 3c. Indeed, the escape rate Γ(t) converges towardsλ 1(B) in the large time limit, and we know it is zero fort < πsince there is no contact be- tween area III and the absorbing boundaryx=x B in the 7 0: 4: 8: 12: 16: 2...

  3. [4]

    Redner,A Guide to First-Passage Processes(Cam- bridge University Press, 2001)

    S. Redner,A Guide to First-Passage Processes(Cam- bridge University Press, 2001)

    2001

  4. [5]

    Chandrasekhar, Stochastic problems in physics and astronomy, Rev

    S. Chandrasekhar, Stochastic problems in physics and astronomy, Rev. Mod. Phys.15, 1 (1943)

    1943

  5. [6]

    A. J. Bray and G. Majumdar, S. N.and Schehr, Persis- tence and first-passage properties in nonequilibrium sys- tems, Adv. Phys.62, 225 (2013)

    2013

  6. [8]

    Majumdar, Brownian functionals in physics and com- puter science, Curr

    S. Majumdar, Brownian functionals in physics and com- puter science, Curr. Sci.89, 2076 (2005)

    2076

  7. [9]

    Godec and R

    A. Godec and R. Metzler, First passage time statistics for two-channel diffusion, J. Phys. A: Math. Theor.50, 084001 (2017)

    2017

  8. [10]

    Shin and A

    J. Shin and A. B. Kolomeisky, Target search on dna by interacting molecules: First-passage approach, J. Chem. Phys.151, 125101 (2019)

    2019

  9. [11]

    S. N. Majumdar and R. M. Ziff, Universal record statis- tics of random walks and l´ evy flights, Phys. Rev. Lett. 101, 050601 (2008)

    2008

  10. [12]

    M. R. Evans and S. N. Majumdar, Diffusion with optimal resetting, J. Phys. A: Math. Theor.44, 435001 (2011)

    2011

  11. [13]

    M. R. Evans, S. N. Majumdar, and G. Schehr, Stochastic resetting and applications, J. Phys. A: Math. Theor.53, 193001 (2020)

    2020

  12. [14]

    Admon, S

    T. Admon, S. Rahav, and Y. Roichman, Experimental realization of an information machine with tunable tem- poral correlations, Phys. Rev. Lett.121, 180601 (2018)

    2018

  13. [15]

    Tal-Friedman, A

    O. Tal-Friedman, A. Pal, A. Sekhon, S. Reuveni, and Y. Roichman, Experimental realization of diffusion with stochastic resetting, J. Phys. Chem. Lett.11, 7350 (2020)

    2020

  14. [16]

    Vatash and Y

    R. Vatash and Y. Roichman, Many-body colloidal dy- namics under stochastic resetting: Competing effects of particle interactions on the steady state distribution (2025), arXiv:2504.10015 [cond-mat.soft]

    Pith/arXiv arXiv 2025

  15. [17]

    Besga, A

    B. Besga, A. Bovon, A. Petrosyan, S. N. Majumdar, and S. Ciliberto, Optimal mean first-passage time for a brow- nian searcher subjected to resetting: Experimental and theoretical results, Phys. Rev. Res.2, 032029 (2020)

    2020

  16. [18]

    Besga, F

    B. Besga, F. Faisant, A. Petrosyan, S. Ciliberto, and S. N. Majumdar, Dynamical phase transition in the first- passage probability of a brownian motion, Phys. Rev. E 104, L012102 (2021)

    2021

  17. [19]

    Faisant, B

    F. Faisant, B. Besga, A. Petrosyan, S. Ciliberto, and S. N. Majumdar, Optimal mean first-passage time of a brow- nian searcher with resetting in one and two dimensions: experiments, theory and numerical tests, J. Stat. Mech. 2021, 113203 (2021)

    2021

  18. [20]

    Archambault, C

    A. Archambault, C. Crauste-Thibierge, A. Imparato, S. Ciliberto, and L. Bellon, First passage time for an underdamped harmonic oscillator and application to the power of an information engine (2026), companion letter, with a focus on short time behavior and an application to the power of an information engine, arXiv:2607.01404 [cond-mat.stat-mech]

    Pith/arXiv arXiv 2026

  19. [21]

    Archambault, C

    A. Archambault, C. Crauste-Thibierge, A. Imparato, S. Ciliberto, and L. Bellon, Ancillary movies obtained by direct numerical simulations of the Langevin Eq. 1, show- ing the phase space evolution forB= 1 andB= 2, with two examples of quality factors:Q= 7 andQ= 100, arxiv.org/src/2607.01405/anc (2026)

    Pith/arXiv arXiv 2026

  20. [22]

    For energy that is quadratic in position and speed, the relaxation time is half ofτ r:τ E r =τ r/2 =Q

  21. [23]

    C. W. Gardiner,Handbook of Stochastic Methods for Physics, Chemistry, and the Natural Sciences, 3rd ed., Springer Series in Synergetics (Springer, 2004)

    2004

  22. [24]

    Risken,The Fokker-Planck Equation: Methods of So- lution and Applications, Springer Series in Synergetics (Springer Berlin Heidelberg, 1996)

    H. Risken,The Fokker-Planck Equation: Methods of So- lution and Applications, Springer Series in Synergetics (Springer Berlin Heidelberg, 1996)

    1996

  23. [25]

    Kramers, Brownian motion in a field of force and the diffusion model of chemical reactions, Physica7, 284 (1940)

    H. Kramers, Brownian motion in a field of force and the diffusion model of chemical reactions, Physica7, 284 (1940)

    1940

  24. [26]

    Zwanzig,Nonequilibrium Statistical Mechanics(Ox- ford University Press, New-York, 2001)

    R. Zwanzig,Nonequilibrium Statistical Mechanics(Ox- ford University Press, New-York, 2001)

    2001

  25. [27]

    H¨ anggi, P

    P. H¨ anggi, P. Talkner, and M. Borkovec, Reaction-rate theory: fifty years after kramers, Rev. Mod. Phys.62, 251 (1990)

    1990

  26. [28]

    Saadat, I

    E. Saadat, I. Latella, and S. Ruffo, Lifetime of locally sta- ble states near a phase transition in the thirring model, J. Stat. Mech.2023, 083207 (2023). 10

    2023

  27. [29]

    Bouchet, C

    F. Bouchet, C. Nardini, and T. Tangarifec, Non- equilibrium statistical mechanics of the stochastic Navier-Stokes equations and geostrophic turbulence, in 5th Warsaw School of Statistical Physics, edited by B. Ci- chocki, M. Napi´ orkowski, and J. Piasecki (Warsaw Uni- versity Press, 2014) p. 3

    2014

  28. [30]

    M. L. Rosinberg, G. Tarjus, and T. Munakata, Heat fluc- tuations for underdamped langevin dynamics, EPL (Eu- rophysics Letters)113, 10007 (2016)

    2016

  29. [31]

    Barros, S

    N. Barros, S. Whitelam, S. Ciliberto, and L. Bellon, Learning efficient erasure protocols for an underdamped memory, Phys. Rev. E111, 044114 (2025)

    2025