{"id":"7e8c3e58-32a5-4c94-a8d1-78b1ef124ace","arxiv_id":"2411.09589","paper_version":2,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"In a damped quantum harmonic oscillator, any initial population distribution whose first r moments match those of the thermal equilibrium relaxes at rate 2γ(r+1), producing a pronounced Mpemba effect.","lead":"This paper shows that a quantum oscillator coupled to a warm bath can relax to equilibrium faster when its initial energy distribution matches certain statistical moments of the equilibrium, a quantum version of the Mpemba effect. The result gives a simple, exact criterion for when such accelerated thermalization happens, which could help design faster qubit resets and state preparations.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"General criterion in Sec. 4 is misstated: Eq. (35) requires population moments to vanish rather than match thermal moments; vacuum with nth=2 and h=3 is a counterexample.","rationale":"The reader's weakest_assumption (Markovianity of the Lindblad master equation) is a recognized model limitation rather than an internal mathematical error, and the paper is explicit about that modeling framework. The central spectral/moment derivation for populations checks out, and the population Mpemba construction is sound. The load-bearing flaw is in the generalized coherence theorem, which the abstract advertises as a general criterion for super-accelerated thermalization: Eq. (35) demands zero population moments for s=0 instead of matching the thermal moments. This is an internal inconsistency, demonstrable by the vacuum-state counterexample, and it is repeated in Appendix B.2. Because the main Sec. 3.3 result and the first Sec. 4 theorem (conditions 1 and 2) remain valid, the paper should not be rejected; it should be accepted conditionally after correcting the statement of the general theorem and verifying the fix against the counterexample. The numerical simulations and the core physics are unaffected by this correction.","tokens_in":20786,"tokens_out":20118,"duration_ms":188616,"concrete_test":"Verify the counterexample: set nth=2, h=3, ρ_i=|0⟩⟨0|. Check that Eq. (35) holds for all allowed (s,l). Compute the exact population solution of Eq. (4); since C1=1−⟨n⟩/nth=1, the spectral decomposition contains e^{−2γt}, hence the relaxation rate is 2γ<3γ, falsifying the theorem as stated. Re-run the same check with the corrected condition ∑ n^l(ρ_{n,n}(0)−P_n^{(S)})=0: vacuum is excluded, and the theorem becomes consistent with Sec. 3.3.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The main Sec. 3.3 result—initial populations with first r moments equal to thermal moments relax at spectral rate 2γ(r+1)—appears correct. However, the general theorem in Sec. 4 (Eq. 35) and its Appendix B.2 version are internally inconsistent. For s=0, l≥1, Eq. (35) demands ∑ n^l ρ_{n,n}(0)=0, i.e. vanishing population moments. The correct moment-matching condition is ∑ n^l ρ_{n,n}(0)=∑ n^l P_n^{(S)} (equivalently, vanishing moments of the deviation ρ−P^{(S)}), as proven in Sec. 3.3 and Appendix A.2. Counterexample: nth=2, h=3, initial vacuum |0⟩⟨0|. All conditions in Eq. (35) are satisfied (mean=0, coherences zero), but the population block has C1=1≠0, so the slowest decaying mode has rate 2γ, not ≥3γ as the theorem claims. Appendix B.2 repeats the same error when it says Q_l^{(0)}(0)=0 reproduces Appendix A.2. The fix is to state the s=0 condition as matching thermal moments, or to restrict the general theorem to s≥1 and treat populations separately.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper analyzes thermalization of a damped quantum harmonic oscillator coupled to a thermal bath, described by the standard Markovian Lindblad master equation. Using generating-function methods for the spectral problem and a moment hierarchy, it shows that an initial diagonal state whose first r population moments equal those of the thermal equilibrium distribution does not excite the r slowest population eigenmodes and therefore relaxes at the accelerated rate 2γ(r+1). The same moment matching is used to construct explicit Mpemba-effect examples, comparing thermal states at nearby temperatures with specially prepared Fock-state mixtures. The paper then extends the discussion to initial states with coherences and states a general theorem, Eq. (35), intended to give a complete criterion for super-accelerated thermalization of the full density matrix. Numerical diagonalization of the Pauli master equation supports the population-subspace results.","tokens_in":21023,"tokens_out":8745,"duration_ms":85388,"significance":"The population-subspace result is a clean, parameter-free analytical criterion: matching the first r moments to the thermal distribution removes the r slowest population eigenmodes, giving a relaxation rate growing linearly in r. This is derived self-consistently and is confirmed numerically, and the constructed Mpemba examples are explicit and convincing. The paper therefore makes a useful contribution to the quantum Mpemba literature. However, the advertised general criterion in Sec. 4 is not correct as stated for the population block, and since this criterion is highlighted in the abstract and conclusions as one of the central results, the error is load-bearing and must be repaired.","major_comments":[{"comment":"The general theorem is stated incorrectly for the s=0 block. Eq. (35) requires, for s=0 and l=1,2,... with 2l<h, that sum_n n^l rho_{n,n}(0)=0. This imposes vanishing initial population moments, not equality to the thermal moments. The correct population condition, proven in Sec. 3.3 and Appendix A.2, is sum_n n^l rho_{n,n}(0)=sum_n n^l P_n^{(S)}. Consequently the theorem as written is false. A concrete counterexample is n_th=2, h=3, and the initial vacuum state |0><0|: all conditions in Eq. (35) are satisfied because all coherences vanish and the only population condition is the vanishing first moment, but Eq. (26) gives C1=1, so the slowest decaying population mode has rate 2 gamma, not the claimed rate >=3 gamma. The theorem should either be restricted to s>=1, with the population block treated separately by the moment-matching condition, or Eq. (35) should be replaced by the condition that the moments of rho_{n,n}(0)-P_n^{(S)} vanish.","section":"Sec. 4, Eq. (35); Appendix B.2, Eq. (82)"},{"comment":"The statement that the coherence theorem 'also holds for s=0 and l=1,2,3,..., i.e. for populations, reproducing the result given in Appendix A.2' is directly contradicted by Appendix A.2 itself. In Appendix A.2 the accelerated-relaxation condition is Eq. (63), Q_l(0)=Q_l^{(S)} for l=1,...,r, not Q_l(0)=0. The vanishing-moment condition used for coherences cannot reproduce the population result, because the thermal distribution has nonzero population moments. This sentence must be corrected together with Eq. (35).","section":"Appendix B.2, final paragraph"},{"comment":"Because Eq. (35) is presented as the paper's general criterion for super-accelerated thermalization, the error in the s=0 case affects the central claim of the paper as summarized in the abstract, not merely a peripheral lemma. The population result remains valid, but the generality claimed for the coherence theorem should be restated carefully, for instance by formulating separate conditions for the population block and for the coherence blocks s>=1.","section":"Abstract and Sec. 4"}],"minor_comments":[{"comment":"The notation P_n^{(S)}(0) on the right-hand side is misleading because the stationary distribution does not depend on time; it should read P_n^{(S)}.","section":"Sec. 3.3, Eq. (28)"},{"comment":"The footnote says 'the condition given by Eq. (25) is never satisfied' when discussing distributions without finite moments; the intended reference appears to be Eq. (27), the condition C1=0, rather than Eq. (25).","section":"Sec. 3.3, footnote 2"},{"comment":"The numerical truncation at Nmax=1800 is adequate for the reported nth values, but the paper does not state how the truncation error was checked for the slowest decay exponents; a brief convergence statement would be helpful.","section":"Sec. 3.4"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The main population result is right and worth knowing: for the damped harmonic oscillator, if the first r moments of the initial population distribution match those of the thermal equilibrium, the relaxation rate jumps to 2γ(r+1). The generating-function derivation is clean, the moment equations check out, and the numerics back the claims. This is a genuine subfield-level advance—it turns scattered observations about quantum Mpemba effects into an exact criterion in a paradigmatic model.\n\nBut the general theorem in Sec. 4, Eq. (35), overreaches. As written, for s=0 it demands that the population moments ∑ n^l ρ_{n,n}(0) vanish, not that they match the thermal moments. The correct condition is matching, as proven in Sec. 3.3 and App. A.2. The counterexample is real: take nth=2, h=3, and the vacuum state. All conditions in Eq. (35) are satisfied (mean zero, no coherences), yet the populations decay at rate 2γ, not ≥3γ, because C1=1. The bullet-point theorem just above Eq. (35) is fine, since it states the matching condition explicitly and restricts coherences to |m-n|>2r+1. The error is in the \"general criterion\" and the corresponding claim in App. B.2 that it reproduces App. A.2—it does not.\n\nThis is a localized but load-bearing flaw for the coherence part. It does not sink the population result, and the fix is straightforward: state the s=0 case as moment matching, or restrict Eq. (35) to s≥1 and treat populations separately. The paper should not be rejected for it, but the theorem should be corrected. The heavy-tailed caveat is honestly acknowledged in Sec. 3.3, and the Markovian assumption is standard for the model.\n\nFor readers working on quantum Mpemba effects or open-system thermalization, this paper is useful and deserves a serious referee. I would cite the population result, and I'd bring it to reading group. With the Sec. 4 correction, it's publishable.","headline":"Solid population moment-matching result, but the general coherence theorem in Eq. (35) is misstated for s=0 and needs a fix before it can be used as stated.","tokens_in":21559,"tokens_out":2366,"would_cite":true,"duration_ms":21459,"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.30.-d","42.50.-p"],"model":"deepseek-v4-flash","headline":"Matching the first r moments of the thermal state makes a damped quantum oscillator relax r+1 times faster.","keywords":["Mpemba effect","quantum thermalization","damped quantum harmonic oscillator","method of moments","Lindblad master equation","super-accelerated relaxation","population coherences","bosonic mode"],"falsifier":"Numerically integrate the Pauli master equation starting from a two-Fock mixture $p|0\\rangle\\langle 0|+(1-p)|n_1\\rangle\\langle n_1|$ with $p=n_{\\rm th}/n_1$ so the mean is exactly $n_{\\rm th}$; if the Kullback-Leibler distance to equilibrium decays with exponent $4\\gamma t$ rather than $8\\gamma t$, the vanishing of the slow-mode amplitude predicted by the theory is disproved.","tokens_in":20565,"feed_emoji":"⚛️","tokens_out":11945,"duration_ms":97283,"temperature":0.7,"pith_summary":"The paper studies the simplest open quantum system, a harmonic oscillator damped by a thermal bath, and asks when it thermalizes faster than expected. Its central claim is that if the initial population distribution has its first r moments exactly equal to those of the equilibrium thermal distribution, the oscillator relaxes to equilibrium at rate $2\\gamma(r+1)$, one step faster for every matched moment. Because a distribution can sit far from equilibrium and still match a few moments, this produces a genuine quantum Mpemba effect inside a paradigmatic model: farther-from-equilibrium states can overtake nearer ones. The statement is proven exactly by the method of moments and the Lindblad spectrum, and it extends to a broad class of initial states with coherences under explicit conditions.","feed_headline":"Matching r equilibrium moments speeds relaxation by factor r+1","feed_subtitle":"Quantum Mpemba effect appears in a damped harmonic oscillator whenever its first r population moments match the bath's.","key_machinery":"The carrier of the argument is the moment hierarchy of the Pauli master equation for the populations: the $l$-th moment $Q_l(t)$ evolves only through lower moments, so if the first $r$ moments start at their equilibrium values they remain locked there while the $(r+1)$-th moment relaxes at rate $2\\gamma(r+1)$. Spectrally, this is backed by the exactly solvable Lindbladian spectrum ($\\lambda_\\alpha=-2\\gamma\\alpha$ in the population sector and $\\lambda^{(s)}_\\alpha=-2\\gamma(\\alpha+|s|/2)+i\\omega_0 s$ for coherences), whose left eigenoperators are polynomials in the number operator — polynomial eigenfunctions are exactly why moment matching cancels the slow modes.","core_discovery":"For the damped quantum harmonic oscillator in a Markovian thermal bath at temperature $T$, the population spectrum of the Lindbladian is $\\lambda_\\alpha=-2\\gamma\\alpha$ and the coherence spectrum is $\\lambda^{(s)}_\\alpha=-2\\gamma(\\alpha+|s|/2)+i\\omega_0 s$. The paper proves that any initial population distribution $P_n(0)$ whose first $r$ moments match those of the thermal distribution $P^{(S)}_n$ has zero amplitude in the $r$ slowest population modes, so it relaxes with leading rate $|\\lambda_{r+1}|=2\\gamma(r+1)$. With coherences present, the same super-accelerated relaxation survives provided $\\rho_{n,m}(0)=0$ for $|m-n|\\le 2r+1$; a general moment theorem on the transformed coherences guarantees relaxation no slower than $\\gamma h$ for any positive integer $h$. This is the quantum Mpemba effect in a foundational model.","pith_inferences":["Editorial inference: the same triangular moment hierarchy should appear in any Markovian thermalization map whose slow eigenmodes are polynomials of a conserved observable; the acceleration mechanism is therefore more general than the harmonic-oscillator model itself.","Editorial inference: an experiment could verify the mechanism by monitoring photon-number moments rather than the full distribution — matched low-order moments should be frozen in time while the first unmatched moment decays at the predicted rate.","Editorial inference: because the argument requires finite moments, distributions with power-law tails cannot be accelerated; a truncated experimental approximation of such a distribution should show the acceleration degrade as the cutoff is lowered.","Editorial inference: the crossing time in the paper's $n_{\\rm th}=2$ example, near $\\gamma t \\simeq 0.28$, offers a sharp quantitative benchmark for cavity or trapped-ion implementations of the effect."],"forward_implications":["A thermal initial state at any temperature $T'$ close to $T$ always excites the slowest population mode, so it relaxes at rate $2\\gamma$; any state with the first $r$ moments matched relaxes at $2\\gamma(r+1)$, making the Mpemba crossing generic.","Coherence does not destroy the speedup: if the low-lying coherences up to $|m-n|\\le 2r+1$ are absent initially, the full density matrix still relaxes at a rate no slower than $2\\gamma(r+1)$.","At zero bath temperature, population-only acceleration disappears because all thermal moments vanish, so the zero-temperature bosonic Mpemba effect must be carried by coherences rather than by population moments.","Moment-matched initial states provide a practical fast-reset mechanism for quantum devices: a state engineered to match a few thermal moments of the surrounding bath will reset faster than a nearby thermal state."],"supporting_citations":[{"why":"Supplies the Lindblad master equation for the damped harmonic oscillator coupled to a thermal bath, the model from which everything follows.","marker":"[39]"},{"why":"Provides the Davies-map criterion for accelerated thermalization that motivates rotating to diagonal states and analyzing the population subspace.","marker":"[29]"},{"why":"Establishes the spectral-amplitude framework for exponentially accelerated approach to stationarity in Markovian open quantum systems.","marker":"[16]"},{"why":"Defines the Kullback-Leibler distance used as the equilibrium-distance measure in the Mpemba comparisons.","marker":"[15]"},{"why":"Underpins the spectral theory of the semi-infinite Jacobi matrix used to solve the population relaxation problem.","marker":"[57]"},{"why":"Connects the birth-death spectral problem to orthogonal polynomials, supporting the analytic form of the eigenfunctions.","marker":"[58]"},{"why":"Provides the zero-temperature bosonic Mpemba counterpart that sets the contrast for the finite-temperature population result.","marker":"[35]"}],"fun_headline_variants":["r matched moments make quantum relaxation superfast","Quantum Mpemba: match r moments, relax with rate r+1","Super-accelerated thermalization via moment matching","Damped oscillator's Mpemba effect: faster when further from equilibrium","Moment matching speeds up quantum thermalization"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The result assumes the oscillator is weakly coupled to a memoryless thermal bath so the Lindblad equation is valid, and that the initial population distribution has finite moments.","fun_headline_variants_meta":{"raw":{"variants":["r matched moments make quantum relaxation superfast","Quantum Mpemba: match r moments, relax with rate r+1","Super-accelerated thermalization via moment matching","Damped oscillator's Mpemba effect: faster when further from equilibrium","Moment matching speeds up quantum thermalization"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000216,"raw_usage":{"total_tokens":1477,"prompt_tokens":1038,"completion_tokens":439,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":654,"completion_tokens_details":{"reasoning_tokens":359}},"tokens_in":654,"tokens_out":439,"duration_ms":4059,"temperature":1.0,"reasoning_tokens":359,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T20:28:26.323476+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Numerically integrate the Pauli master equation starting from a two-Fock mixture $p|0\\rangle\\langle 0|+(1-p)|n_1\\rangle\\langle n_1|$ with $p=n_{\\rm th}/n_1$ so the mean is exactly $n_{\\rm th}$; if the Kullback-Leibler distance to equilibrium decays with exponent $4\\gamma t$ rather than $8\\gamma t$, the vanishing of the slow-mode amplitude predicted by the theory is disproved.","supporting_citations":[{"cited_title":"Breuer and F","cited_arxiv_id":null,"evidence_quote":"Supplies the Lindblad master equation for the damped harmonic oscillator coupled to a thermal bath, the model from which everything follows."},{"cited_title":"Thermodynamics of the Quan- tum Mpemba Effect","cited_arxiv_id":null,"evidence_quote":"Provides the Davies-map criterion for accelerated thermalization that motivates rotating to diagonal states and analyzing the population subspace."},{"cited_title":"Exponentially Accelerated Approach to Stationarity in Markovian Open Quantum Systems through the Mpemba Effect","cited_arxiv_id":null,"evidence_quote":"Establishes the spectral-amplitude framework for exponentially accelerated approach to stationarity in Markovian open quantum systems."},{"cited_title":"Equidistant quenches in few-level quantum systems","cited_arxiv_id":null,"evidence_quote":"Defines the Kullback-Leibler distance used as the equilibrium-distance measure in the Mpemba comparisons."},{"cited_title":"Spectral Theory for the Differential Equations of Simple Birth and Death Processes","cited_arxiv_id":null,"evidence_quote":"Underpins the spectral theory of the semi-infinite Jacobi matrix used to solve the population relaxation problem."},{"cited_title":"Linear birth and death models and asso- ciated Laguerre and Meixner polynomials","cited_arxiv_id":null,"evidence_quote":"Connects the birth-death spectral problem to orthogonal polynomials, supporting the analytic form of the eigenfunctions."},{"cited_title":"Bosonic Mpemba effect with non-classical states of light","cited_arxiv_id":null,"evidence_quote":"Provides the zero-temperature bosonic Mpemba counterpart that sets the contrast for the finite-temperature population result."}],"review_version":1}