{"id":"b0c17074-e459-4d6c-bc8b-f32083f85081","arxiv_id":"2608.08267","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"For a harmonic oscillator driven by classical noise, the first-passage-time distribution to an energy threshold under projective measurements is exactly a sum of exponentials with Laguerre-root rates, converging to the classical result only at high thresholds.","lead":"This paper derives the statistics of the time a noisy quantum oscillator first crosses an energy threshold under repeated projective measurements, and compares them to the classical case. It finds the two agree at high thresholds but differ sharply at low ones, and that surviving quantum trajectories develop persistent Wigner negativity.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Appendix C passivity-preservation lemma is algebraically wrong: the correct gap dynamics have a negative diagonal term, so passive states are not generally preserved by (10).","rationale":"The reader's weakest assumption concerns the stochastic-averaging approximation used for the classical FPTD. That is a caveat the authors explicitly acknowledge, and both sides of the comparison assume the same weak-coupling limit; it does not obviously threaten Eq. (11). The passivity error is a concrete internal inconsistency that the reader missed. Eq. (11) itself checks out: small-N_B Laplace transforms match, and the residue expansion yields the stated Laguerre coefficients. The Monte Carlo and master-equation FPTDs agree, though no code is shipped. However, the passivity lemma underpins the assertion that ensemble-averaged Wigner functions are nonnegative for any passive initial state, which supports the ensemble-versus-trajectory interpretation of Wigner negativity. Since the lemma is false, the paper needs a correction or a restriction of the claim (e.g., to ground-state or thermal initial states). This does not overturn the central FPTD and crossover claims, so the reader's CONDITIONAL verdict remains appropriate, but the reasons should include this mathematical error.","tokens_in":19671,"tokens_out":33970,"duration_ms":300069,"concrete_test":"Numerically integrate (10) for N_B = 3 with initial vector (0.5, 0.5, 0). If P_1(t) − P_2(t) becomes negative, which occurs near t = 0.25, the passivity-preservation lemma is false. Alternatively, re-derive the gap equation symbolically from (10) and check for the −2(n+1)D_n term; if present, the paper's RHS is incomplete and the 'any passive initial state' claim must be retracted. For completeness, verify that the ground-state solution remains passive so that Figs. 6 and 8 are unaffected.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In Appendix C the paper claims that (10) preserves passivity, stating ∂_t(P_n − P_{n+1}) = n(P_{n−1} − P_n) + (n+2)(P_{n+1} − P_{n+2}) ≥ 0. This omits a diagonal term. Defining D_n = P_n − P_{n+1}, equation (10) gives ∂_t D_n = nD_{n−1} − 2(n+1)D_n + (n+2)D_{n+1}. The −2(n+1)D_n term lets D_n decrease even when all gaps are nonnegative. Explicit counterexample with N_B ≥ 3 and initial populations P_0 = P_1 = 1/2, P_2 = 0 (a passive state): at t = 0, P_1′ = P_0 − 3P_1 = −1 and P_2′ = 2P_1 = 1, so (P_1 − P_2)′ = −2 and D_1 crosses zero near t = 1/4. Thus the conditioned population dynamics do not preserve the set of passive states, and the conclusion that every passive initial state has a nonnegative ensemble Wigner function is unsupported. The ground-state results may survive because the free solution from |0⟩ is a thermal (geometric) distribution and the absorbing boundary only truncates the tail, but the general lemma is false as stated.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"Ryan et al. study the first-passage-time statistics of a harmonic oscillator driven by classical Gaussian white noise and subject to repeated projective measurements of an energy threshold E_B = N_B + 1/2. The paper defines a stroboscopic protocol, shows numerically that finite-interval results converge as θ → 0, and derives in that limit an exact first-passage-time distribution for initial Fock states: P_0^FPT(t|N0,NB) = Σ_i [L_{N0}(λ_i)/L_{NB-1}^{(1)}(λ_i)] e^{-λ_i t}, where λ_i are the roots of L_{NB}(x). This quantum result is compared with the first-passage distribution of a classically stochastically averaged energy diffusion, and the difference is quantified with a Kullback-Leibler divergence, yielding the claimed quantum-classical crossover. The paper further analyzes conditioned surviving states, including sub-Poissonian statistics, quasi-stationary Wigner functions, and trajectory-resolved Wigner negativity, and verifies that density-matrix and Monte Carlo trajectory descriptions agree.","tokens_in":19911,"tokens_out":35633,"duration_ms":368298,"significance":"If the central result stands, this is a rare exactly solvable quantum first-passage problem with a closed-form distribution, connecting monitored quantum dynamics to the spectral theory of a truncated birth-death chain and to the classical Bessel limit via the Laguerre-to-Bessel correspondence. The paper provides concrete, experimentally testable predictions for trapped-ion platforms and identifies monitored first passage as a route to heralded non-Gaussian resource states. Strengths include the absence of fitted parameters in the central formula, the analytic treatment of the θ → 0 limit, the detailed appendices, and the cross-validation of master-equation and Monte Carlo results in Fig. 12. The main caveats are the approximate classical benchmark used for the crossover claim and a proof in Appendix C that needs correction.","major_comments":[{"comment":"The classical benchmark P_C^FPT is obtained from the stochastically averaged energy diffusion of Eq. (B1), not from the original two-dimensional Langevin equation (1). The text states this approximation, but no quantitative validation is provided for the parameter regimes used in Figs. 4 and 6. Because the KL divergences in Fig. 4b and the quantitative location of the quantum-classical crossover are computed relative to this approximate benchmark, the authors should validate the approximation against direct numerical simulation of Eq. (1), or otherwise bound its error. This does not affect the exact quantum formula (11), but it is load-bearing for the quantitative crossover claim.","section":"Appendix B and Fig. 4"},{"comment":"The derivation of the passivity-preservation lemma is not written correctly. The gap variable D_n = P_n - P_{n+1} satisfies ∂_t D_n = n D_{n-1} - 2(n+1)D_n + (n+2)D_{n+1}; the expression in Appendix C omits the diagonal term -2(n+1)D_n. The argument only needs to evaluate this at the boundary D_n = 0, where the omitted term vanishes, so the conclusion may be salvageable (for example, by noting that the gap generator is Metzler and hence preserves the nonnegative cone), but the printed identity and the inequality following it are not a valid proof as they stand. This matters because the passivity lemma underpins the claim that the ensemble Wigner function remains nonnegative for all passive initial states and the statement that ⟨N(t)⟩_ρ = 0 in Fig. 8.","section":"Appendix C and Sec. III B 2"}],"minor_comments":[{"comment":"The residue expansion should state explicitly that L_{NB-1}^{(1)} is the associated Laguerre polynomial and that the calculation uses L_{NB}'(λ_i) = -L_{NB-1}^{(1)}(λ_i); as written, the sign convention is easy to misread.","section":"Appendix A, Eq. (A5)"},{"comment":"The numerical master-equation calculations should state the Fock-space truncation used, since the infinite-dimensional Hilbert space must be truncated to evaluate e^{θM} and the quasi-stationary states.","section":"Sec. III and Appendix E"},{"comment":"The validity condition for stochastic averaging, described in the text as 'the heating per oscillation is small', should be quantified in terms of ω0 and D; currently no explicit parameter relation is given, which makes it hard to assess the accuracy of the classical benchmark in the plotted regimes.","section":"Appendix B"},{"comment":"The text notes that points with N_B - N_0 = 1 are omitted because of a logarithmic divergence in the KL integrand; a brief comment on how sensitive the displayed KL values are to the short-time behavior would help the reader interpret the crossover quantitatively.","section":"Fig. 4b"}],"recommendation":"major_revision","confidential_remarks":"The paper is well within scope for a quantum-physics journal. The self-citation of the experimental work [31] is used as motivation and testability claim, not as an input to the derivations, so I see no circularity concern. The central exact result, Eq. (11), appears sound and is the main contribution; the requested revisions concern the validation of the classical comparison and the rigor of the passivity proof. If the authors provide the suggested numerical validation and correct the Appendix C derivation, I would support acceptance."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things you should know. First, this is a serious analytic paper with a real new result: the closed-form FPTD for the continuously monitored noisy oscillator (Eq. 11), the Laguerre-to-Bessel crossover, and the trajectory-level Wigner negativity. The derivations in Appendices A and B are careful, and the numeric cross-checks (master equation versus Monte Carlo, finite-θ versus θ→0) are convincing. Second, the passivity-preservation lemma in Appendix C is algebraically wrong, and that error is load-bearing for the general claim that every passive initial state has a nonnegative ensemble Wigner function.\n\nThe issue: with D_n = P_n − P_{n+1}, Eq. (10) gives ∂t D_n = nD_{n−1} − 2(n+1)D_n + (n+2)D_{n+1}. The paper drops the middle term and concludes the right-hand side is nonnegative. The counterexample P_0 = P_1 = 1/2, P_2 = 0 (a passive state for N_B ≥ 3) immediately violates passivity: at t = 0, D_1′ = −2. So the set of passive states is not invariant under the conditioned dynamics, and the proof of nonnegative Wigner functions for arbitrary passive initial states collapses. The ground-state specific result may still hold — thermal states are preserved by the bulk generator and the boundary only truncates — but the authors need to prove that separately or weaken the claim. This is a real, localized error, not a manufactured one.\n\nThe rest of the paper is in good shape. The FPTD formula (11) is derived without fitting, the KL-divergence crossover figure is clean, and the authors are honest that the classical side uses stochastic averaging. No code or data files are shipped, which is a minor reproducibility gap for the numerics, not a correctness issue. The self-citation to the experimental paper [31] is appropriate as motivation; nothing is assumed from it.\n\nWho this is for: anyone working on monitored bosonic modes, first-passage statistics in open quantum systems, or trapped-ion state engineering. It deserves a serious referee. My recommendation: send it to review, but make fixing Appendix C a condition of acceptance. The main result can stand; the overbroad passivity claim cannot.","headline":"Solid new analytic FPTD result for monitored oscillators, but the passivity-preservation lemma in Appendix C is algebraically wrong and needs fixing before acceptance.","tokens_in":20504,"tokens_out":4264,"would_cite":true,"duration_ms":38780,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["03.65.Yz","05.40.-a"],"model":"deepseek-v4-flash","headline":"Repeated energy-threshold measurements give a noisy oscillator an exactly solvable escape-time distribution — Laguerre exponentials — that deviates from classical at low thresholds and merges back as the threshold rises.","keywords":["quantum first-passage time","energy-barrier measurement","monitored harmonic oscillator","Laguerre polynomial","quantum-classical crossover","Wigner negativity","projective measurement","trapped-ion oscillator"],"falsifier":"A direct table-top test: initialize a trapped-ion or cavity oscillator in the ground state, apply the energy-barrier measurement at $N_B=1$ and $N_B=2$ quanta, and record escape times for a sweep of measurement intervals $\\theta$. The quantum prediction is parameter-free: at $N_B=1$ the distribution is geometric with $p=\\theta/(\\theta+1)$ and hence a pure exponential of rate 1 in the continuous limit; at $N_B=2$ with $N_0=0$ it is $(1/\\sqrt{2})(e^{-(2-\\sqrt{2})t}-e^{-(2+\\sqrt{2})t})$. If the observed shape, the two rate constants, or the $\\theta$-scaling disagree with these forms, Eq. (11) is wrong; agreement while the classical Bessel-mode prediction at the same barrier visibly deviates would confirm the low-threshold quantum-classical split.","tokens_in":19419,"feed_emoji":"⏱️","tokens_out":21781,"duration_ms":186492,"temperature":0.7,"pith_summary":"This paper asks whether the timing of a random event — the first time a noisy oscillator crosses an energy threshold — carries quantum information. For a harmonic oscillator driven by classical white noise and probed by repeated projective energy-barrier measurements, it derives an exact first-passage-time distribution in the continuous-measurement limit: a finite sum of exponentials whose rates are the roots of a Laguerre polynomial. The central finding is that at low thresholds these quantum timing statistics differ measurably from the classical prediction, while at high thresholds the difference gradually vanishes. This matters because the timing statistics become effectively classical even while the conditioned quantum states stay nonclassical, making a monitored oscillator a clean testbed for the quantum-classical crossover and for measurement-induced generation of nonclassical states.","feed_headline":"Escape-time statistics go quantum at low energy thresholds","feed_subtitle":"Exact escape-time formula: where measurement timing turns nonclassical, and where it merges back to classical.","key_machinery":"The load-bearing object is the truncated generator $M^{(N_B)}=P_S M P_S$, the master-equation superoperator restricted to the surviving Fock subspace; in the continuous-measurement limit it acts on occupations as the birth-death chain $\\partial_t P_n=(n+1)P_{n+1}-(2n+1)P_n+nP_{n-1}$ with an absorbing boundary $P_{N_B}=0$. Because this generator's eigenvalue equation obeys the same three-term recursion as the Laguerre polynomials, the Laplace transform of the escape probability is the rational function $L_{N_0}(-s)/L_{N_B}(-s)$, and the Heaviside expansion theorem converts it into the finite exponential sum that is the paper's main formula. The long-time surviving state is the kernel of the associated Jacobi matrix, with Fock occupations $P_n^{\\mathrm{qs}}\\propto L_n(\\lambda_1^{N_B})$, a superexponential nonthermal envelope that serves as the reference state for the nonclassicality analysis. On the classical side, the same role is played by the stochastically averaged Fokker-Planck operator $dE=dt+\\sqrt{2E}\\,dW$, whose eigenfunctions are Bessel functions with eigenvalues $(j_0^{(i)})^2/4E_B$. In the trajectory picture, the machinery is the random displacement operator $D(\\alpha)$ with $\\alpha\\sim\\mathcal{N}_\\mathbb{C}(0,\\theta)$ interleaved with the threshold projection, and the projective truncation of Fock support is what produces Wigner negativity in individual surviving trajectories.","core_discovery":"Starting from the Langevin equation $\\ddot{x}+\\omega_0^2 x=\\xi(t)$ and its quantum counterpart, the infinite-temperature amplitude-reservoir master equation, the paper analyzes a stroboscopic protocol: evolve for an interval $\\theta$, then projectively measure whether the oscillator occupies a Fock level below the threshold $N_B$, repeating until the first absorbing outcome. Its central claim is that in the limit of infinitely frequent measurements the first-passage-time density for an initial Fock state $|N_0\\rangle$ is exactly $P_{\\mathrm{FPT}}^0(t|N_0,N_B)=\\sum_i [L_{N_0}(\\lambda_i)/L_{N_B-1}^{(1)}(\\lambda_i)] e^{-\\lambda_i t}$, where $\\lambda_i$ are the roots of the Laguerre polynomial $L_{N_B}$. Equivalently, the escape statistics are exactly the spectral content of a truncated birth-death chain with an absorbing boundary at $n=N_B$. The classical counterpart, obtained from the stochastically averaged energy diffusion $dE=dt+\\sqrt{2E}\\,dW$, is an infinite Bessel-mode expansion with eigenvalues $(j_0^{(i)})^2/4E_B$, and the Mehler–Heine formula connects the two, so the quantum-classical crossover is governed by the Laguerre-to-Bessel correspondence. The paper further establishes that the ensemble density-matrix description and the trajectory-resolved Monte Carlo description give indistinguishable timing statistics, that the quantum and classical mean first-passage times coincide ($N_B-N_0$), and that all quantum-classical deviation lives in the shape of the distribution, concentrated at short and intermediate times when the threshold is within a few quanta of the initial state.","pith_inferences":["Because Eq. (11) has no fitting parameters, the escape-time spectrum at tiny thresholds is a set of universal numbers — for $N_B=2$ and $N_0=0$, rates $2\\mp\\sqrt{2}$ with weights $\\pm 1/\\sqrt{2}$ — so a coarse histogram of escape times from a ground-state oscillator tests the entire derivation quantitatively, not merely its qualitative crossover.","The only approximate leg of the comparison is classical: a numerical first-passage-time distribution from the exact two-dimensional Langevin equation (1) would replace the stochastic-averaging result and could shift the KL-divergence crossover line, while leaving the quantum formula untouched; the exact classical FPTD is therefore the natural next computation.","The mechanism looks like a template rather than a one-off calculation: replacing the Fock-space threshold projector with a different measurement would replace the Laguerre recursion with the corresponding orthogonal-polynomial recursion, so 'escape-time statistics equal the spectrum of a truncated birth-death chain' should hold for other monitored bosonic modes and boundary conditions."],"forward_implications":["At $N_B=1$ the continuous-limit first-passage density is a pure exponential with rate 1, and for finite $\\theta$ it is exactly geometric with $p=\\theta/(\\theta+1)$, because every survival measurement projects the oscillator back into the ground state — a signature no classical diffusion model reproduces at such a low barrier.","The mean first-passage time is identical in the quantum and classical treatments ($N_B-N_0$ in the paper's dimensionless units), so the quantum-classical difference lives entirely in the distribution's shape and higher moments, and it is largest when the threshold is only a few quanta above the initial state.","Ensemble density-matrix evolution and trajectory-resolved Monte Carlo wavefunction simulation yield indistinguishable first-passage-time statistics (verified numerically over a range of $N_0$, $N_B$, and $\\theta$), so the trajectory picture is a faithful microscope on the same process.","Survival conditioning alone creates nonclassical states despite purely classical driving noise: sub-Poissonian number statistics, nonpositivity of the $P$-representation, and Wigner negativity along individual trajectories that persists in the quasi-stationary regime; an observer retaining only survival outcomes sees a nonnegative-Wigner ensemble while an observer with the full noise record sees n","The protocol therefore doubles as a heralded preparation scheme for non-Gaussian continuous-variable resource states, realizable in trapped-ion oscillators and by photonic state truncation, and it demonstrates that operational observables can appear classical before the underlying conditioned state dynamics do."],"supporting_citations":[{"why":"Defines the energy-barrier projective measurement protocol on a trapped-ion oscillator that the model describes and that can test Eq. (11) directly.","marker":"[31]"},{"why":"Supplies the infinite-temperature amplitude-reservoir master equation (2), the quantum analogue of the Langevin dynamics (1).","marker":"[35]"},{"why":"Provides the Laguerre-polynomial recursion that makes the truncated generator's eigenvalue equation exactly solvable.","marker":"[37]"},{"why":"Supplies the Heaviside expansion theorem used to invert the Laplace-domain escape probability into the finite exponential sum.","marker":"[63]"},{"why":"Gives the stochastic-averaging reduction of the two-dimensional Langevin oscillator to the one-dimensional energy diffusion (B1) used for the classical FPTD.","marker":"[66]"},{"why":"Provides the Fokker-Planck and Pontryagin machinery for the classical first-passage eigenfunction expansion and moments.","marker":"[68]"},{"why":"Supplies the inequality bounds on the smallest Laguerre root used for the quantum long-time decay rate $\\lambda_1^{N_B}$ and its $1/E_B$ scaling.","marker":"[62]"},{"why":"Gives the displacement-operator representation of the noise (18)-(19) used in the trajectory-resolved simulation.","marker":"[49]"},{"why":"Supplies the Monte Carlo wavefunction method used to show that individual trajectories reproduce the ensemble first-passage-time statistics.","marker":"[51]"}],"fun_headline_variants":["Exact formula pins where quantum escape times leave classical behind","Measurement timekeeping: quantized escape statistics at low thresholds","Exact escape-time formula: quantum measurements set the clock","Quantized escape times emerge under frequent measurement","Quantum-classical split in first-passage statistics: exact formula"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The classical first-passage-time distribution used for the comparison is not computed from the full two-dimensional Langevin equation but from the stochastically averaged one-dimensional energy diffusion $dE=dt+\\sqrt{2E}\\,dW$, which assumes the oscillator's energy is nearly constant over one oscillation period; if that averaging fails in the low-threshold regime, the quantitative KL-divergence values and the precise location of the crossover would shift, even though the quantum formula itself would remain exact.","fun_headline_variants_meta":{"raw":{"variants":["Exact formula pins where quantum escape times leave classical behind","Measurement timekeeping: quantized escape statistics at low thresholds","Exact escape-time formula: quantum measurements set the clock","Quantized escape times emerge under frequent measurement","Quantum-classical split in first-passage statistics: exact formula"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000823,"raw_usage":{"total_tokens":3663,"prompt_tokens":1069,"completion_tokens":2594,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":685,"completion_tokens_details":{"reasoning_tokens":2515}},"tokens_in":685,"tokens_out":2594,"duration_ms":20077,"temperature":1.0,"reasoning_tokens":2515,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T00:13:22.712005+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct table-top test: initialize a trapped-ion or cavity oscillator in the ground state, apply the energy-barrier measurement at $N_B=1$ and $N_B=2$ quanta, and record escape times for a sweep of measurement intervals $\\theta$. The quantum prediction is parameter-free: at $N_B=1$ the distribution is geometric with $p=\\theta/(\\theta+1)$ and hence a pure exponential of rate 1 in the continuous limit; at $N_B=2$ with $N_0=0$ it is $(1/\\sqrt{2})(e^{-(2-\\sqrt{2})t}-e^{-(2+\\sqrt{2})t})$. If the observed shape, the two rate constants, or the $\\theta$-scaling disagree with these forms, Eq. (11) is wrong; agreement while the classical Bessel-mode prediction at the same barrier visibly deviates would confirm the low-threshold quantum-classical split.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the infinite-temperature amplitude-reservoir master equation (2), the quantum analogue of the Langevin dynamics (1)."},{"cited_title":"Leo´ nski and A","cited_arxiv_id":null,"evidence_quote":"Supplies the Heaviside expansion theorem used to invert the Laplace-domain escape probability into the finite exponential sum."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the stochastic-averaging reduction of the two-dimensional Langevin oscillator to the one-dimensional energy diffusion (B1) used for the classical FPTD."},{"cited_title":"H¨ anggi and P","cited_arxiv_id":null,"evidence_quote":"Provides the Fokker-Planck and Pontryagin machinery for the classical first-passage eigenfunction expansion and moments."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the inequality bounds on the smallest Laguerre root used for the quantum long-time decay rate $\\lambda_1^{N_B}$ and its $1/E_B$ scaling."},{"cited_title":"Van Herstraeten and N","cited_arxiv_id":null,"evidence_quote":"Supplies the Monte Carlo wavefunction method used to show that individual trajectories reproduce the ensemble first-passage-time statistics."}],"review_version":1}