{"id":"a1c2ad32-f1eb-4db1-b9a6-165f8e25381b","arxiv_id":"2412.01436","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"high","formal_verification":"none","parameter_count":2,"one_line_summary":"A weakly damped Brownian particle is said to show Kibble-Zurek-like scaling, but the scaling exponents are asserted without derivation and the method's validity is explicitly violated.","lead":"This paper claims that a driven underdamped Brownian particle, in the limit of very weak friction and stiffness, shows a freezing effect similar to the Kibble-Zurek mechanism of quantum phase transitions. It reports scaling exponents for the irreversible work, but the key derivations are missing and the justification relies on the author's own prior work.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Claimed KZ exponents rest on a linear-response calculation the paper itself says is invalid, backed only by self-citation; Eq. 12 under the stated f(τ)>τ condition gives τ^0, not τ^{-2}, and Eq. 3 is not normalized at t=0.","rationale":"The paper's advertised novelty is the new Kibble-Zurek scaling, Eqs. (13)-(14). Those exponents are the load-bearing output; without them the paper is only an observation that in a singular double limit the relaxation time diverges. The derivation path is: linear-response formula Eq. (12), a freeze-out time f(τ), and an explicit admission that f(τ) exceeds the plotted τ range, so the linear-response validity condition is not respected, followed by a self-citation to Ref. [4] as the only justification. The reader's weakest assumption identifies exactly this step. My analysis strengthens it: under the stated inequality f(τ)>τ, Eq. (12) with a bounded normalized relaxation function gives an O(1) impulse work, not the plotted τ^{-2} scaling, so the internal consistency of the sudden-process claim is also in doubt. The additional normalization failure of Eq. (3) at t=0 reinforces that the response function entering the calculation is not reliable. None of this is an attack on the author; it is a description of the argument's exposed load-bearing step. An exact Langevin simulation would settle whether the exponents survive outside linear response, but as it stands the central scaling claim is unsupported. The reader's REJECT verdict is therefore appropriate; my read does not change it.","tokens_in":3097,"tokens_out":17055,"duration_ms":156956,"concrete_test":"Simulate the exact Langevin equation (1) with the same parameters as Figs. 1-2 (γ=ω0=0.01, δω0=0.001, m=1, β=1/γ) for τ spanning 0.002 to 0.01 and 2×10^5 to 10^6, averaging the mean work over enough trajectories to extract d ln W_irr/d ln τ. If the sudden slope is not -2 or the slow slope is not -1, the claimed exponents fail; if they match, the self-citation rescue is empirically confirmed.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central new results are the exponents ηKZ=-2 and -1 in Eqs. (13)-(14), extracted from the impulse-work expression Eq. (12). The paper's own caveat in Sec. C says: 'the solution f(τ) is always greater than the range of τ used in these graphics, which indicates that the range of validity of linear response is not respected.' The only rescue offered is 'previous studies illustrate that the Kibble-Zurek exponent is calculated correctly [4]'—a self-citation with no derivation shown for this system. That assumption is load-bearing: if the exponents are wrong, the central scaling claim collapses. The problem is not merely a vague regime issue. If f(τ)>τ for the plotted ranges, then in Eq. (12) the lower integration limit τ-f(τ) clips to 0; since ρ(t-t')=Ψ0/Ψ0(0) is bounded (about 3 in the γ/ω0=1 limit), W_I ≈ (δω0^2/τ^2)∫_0^τ∫_0^t ρ dt' dt, which is O(1), not τ^{-2}. Thus the sudden-process slope in Fig. 1 is not reproducible from Eq. (12) as written. Additionally, Eq. (3) fails a trivial normalization test: at t=0 the right-hand side equals ω^2/ω0^2=3 for γ/ω0=1, not 1, so the 'frozen value 3' of Eq. (7) is partly an artifact of an incorrect constant in the response function.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript claims that a driven underdamped Brownian motion exhibits the Kibble-Zurek mechanism in the simultaneous limit γ,ω₀→0⁺ with γ/ω₀=1. Using linear response theory, it derives a divergent relaxation time (Eq. 6), a frozen relaxation function (Eq. 7), a high-temperature initial condition (Eq. 9), new scaling exponents η_KZ=-2 and -1 for the impulse part of the irreversible work (Eqs. 13-14), and a pausing effect in the optimal protocol (Eq. 15). The results are presented in a very brief form, with several key steps asserted rather than derived.","tokens_in":3456,"tokens_out":7571,"duration_ms":57900,"significance":"If established, the work would provide a simple classical analog of the Kibble-Zurek mechanism and could be of interest to the stochastic thermodynamics community. However, the paper does not provide a self-contained derivation of its scaling claims; the only support for the central exponents is the author's own prior work, and the manuscript explicitly acknowledges that the linear response calculation is used outside its range of validity. The normalization error in the response function further undermines the quantitative claims. The significance is therefore contingent on a substantial revision that the present manuscript does not provide.","major_comments":[{"comment":"Equation (3) is not a normalized response function: at t=0 its right-hand side equals ω²/ω₀² = 3 when γ/ω₀=1, rather than 1 as required by the definition Ψ₀(0)/Ψ₀(0)=1. This error propagates into the 'frozen' value 3 reported in Eq. (7) and into the impulse-work integrand in Eq. (12), since all subsequent results use Ψ₀(t)/Ψ₀(0). Please re-derive Eq. (3) from the Langevin equation (1) and correct the normalization before any scaling claims can be assessed.","section":"I.A, Eq. (3)"},{"comment":"The paper states that for the plotted ranges f(τ)>τ, so the lower integration limit τ-f(τ) in Eq. (12) is negative and the integral effectively runs from 0 to τ. In that regime the normalized response (Eq. (3) with γ,ω₀ small and γ/ω₀=1) is approximately constant (≈3) over the integration domain, so the double integral is O(τ²) and W_irr^I is O(δω₀²), i.e., independent of τ. This contradicts the claimed η_KZ=-2 scaling in Eq. (13) and the slope in Fig. 1. Please clarify how the numerical results were obtained and whether Eq. (12) is the correct expression in the regime of the figures.","section":"I.C, Eq. (12)"},{"comment":"The scaling exponents are asserted without derivation. No closed form for f(τ) is given, no asymptotic analysis of Eq. (12) is shown, and no direct numerical integration of the original Langevin dynamics is provided. The only justification offered is a citation to the author's previous work [4]. Given that the paper itself admits that the linear response validity condition is violated in the plotted regime, the exponents cannot be considered established unless they are derived within this manuscript or verified by a calculation that does not rely on the invalid approximation.","section":"I.C, Eqs. (13)-(14)"},{"comment":"The 'high-temperature condition' T∝1/γ is stated without derivation or a clear argument from the fluctuation-dissipation theorem. As presented, it is an additional assumption that is singular in the γ→0 limit and places the system in a regime where the validity of linear response around the initial equilibrium state must be re-examined. Please provide a derivation or explicitly label this as a condition that defines the KZ limit rather than a consequence of linear response.","section":"I.B, Eq. (9)"}],"minor_comments":[{"comment":"The paper never defines the Kibble-Zurek exponent for the quantum Ising chain or states the comparison quantitatively. A few sentences defining the standard KZ scaling and how it is generalized here would place the results in context.","section":"Abstract/Introduction"},{"comment":"The axis labels are unclear (e.g., 'τ-2 Wirr I' in Fig. 1). Please use standard notation, e.g., W_irr^I on the y-axis, and state explicitly whether the figures show W_irr^I versus τ on a log-log scale or some rescaled quantity.","section":"Fig. 1 and Fig. 2"},{"comment":"The optimal protocol result g*(t)≈1/2 is taken directly from reference [5] and is not derived in this manuscript. Since the connection to the Kibble-Zurek mechanism is the point of this section, a more explicit explanation of how the pausing effect follows from the preceding analysis would be helpful.","section":"I.D, Eq. (15)"},{"comment":"Reference [1] is missing its article number/page (likely J. Chem. Phys. 129, 024114 (2008)). Please verify all references for completeness.","section":"References"}],"recommendation":"reject","confidential_remarks":"The manuscript is extremely brief for the scope of its claims: the scaling exponents are not derived, the central calculation is acknowledged to be outside the validity regime, and the only support for the key exponents is a self-citation. The normalization error in Eq. (3) suggests a fundamental problem with the response-function derivation. In my view this cannot be fixed by a minor revision; the paper would need a complete rewrite with rigorous derivations or numerical validation. I therefore recommend rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper has a genuinely interesting observation buried under a pile of unsubstantiated claims. The limit γ,ω0→0 with γ/ω0=1 does make the relaxation time diverge, and that is a legitimate curiosity. But the two new exponents, ηKZ=-2 and -1, are asserted, not derived. The only justification given is a self-citation (ref. [4]) after the paper itself concedes that the linear-response regime is not reached. That alone would worry me; what worries me more is that the equations as written do not support the plotted slopes. If f(τ)>τ, the lower limit in Eq. (12) is negative; clipping it to zero gives an O(1) result, not τ^{-2}. And Eq. (3) fails a basic check: for γ/ω0=1, the right-hand side at t=0 equals 3, not 1, so the frozen value in Eq. (7) is at least partly a normalization artifact. These are load-bearing problems, not minor typos. The high-temperature condition T∝1/γ is stated without derivation and looks ad hoc. What the paper does well is identify a corner of the problem worth thinking about: a classical Brownian system whose relaxation time diverges in a particular limit, with a possible analogy to Kibble-Zurek. The idea is not crazy, but the execution is a sketch, not a paper. I would not send this to a referee in its current form. The author needs to supply a derivation of the exponents, fix the response function normalization, and show that the plotted scaling actually follows from the equations. If that can be done, this might become a solid contribution. Right now, it is not. Recommendation: reject, but do not close the door; a rewritten version with the missing steps could deserve a serious look.","headline":"Interesting limit, but the central exponents are asserted without derivation, the response function is misnormalized, and the stated equations contradict the claimed scaling.","tokens_in":3974,"tokens_out":3168,"would_cite":false,"duration_ms":26797,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper argues that a single driven underdamped Brownian particle, in the double limit of vanishing friction and vanishing natural frequency taken at the same rate, exhibits the Kibble-Zurek mechanism: its relaxation time diverges, its…","keywords":["Kibble-Zurek mechanism","underdamped Brownian motion","linear response theory","relaxation time divergence","irreversible work scaling","sudden process","optimal protocol","pausing effect"],"falsifier":"Directly integrate the driven Langevin equation (1) for $\\gamma=\\omega_0=0.01$, $\\delta\\omega_0=0.001$, over switching times $\\tau$ in the same range as Figures 1 and 2, computing the exact irreversible work from the full work distribution without the linear-response approximation, and compare the exponents of $W_{\\rm irr}^I(\\tau)$; if they are not $-2$ (sudden) and $-1$ (slowly varying), the central scaling claim fails.","tokens_in":2879,"feed_emoji":"🧊","tokens_out":5489,"duration_ms":45470,"temperature":0.7,"pith_summary":"This paper argues that a single driven underdamped Brownian particle, when the friction and the natural frequency are made small at the same rate, behaves like a Kibble-Zurek system: its relaxation time diverges, its relaxation function freezes to a constant, and the irreversible work obeys a power law in the switching time. The paper derives the exponents $\\eta_{KZ}=-2$ for sudden drives and $\\eta_{KZ}=-1$ for slowly varying drives using linear response theory. This matters because it moves a paradigmatic nonequilibrium scaling from quantum chains at zero temperature to a minimal classical stochastic system, potentially making the mechanism easier to observe and control.","feed_headline":"Underdamped Brownian motion freezes like a Kibble-Zurek system","feed_subtitle":"Relaxation time diverges and irreversible work obeys power-law scalings, bringing a quantum-chain effect to a single classical oscillator.","key_machinery":"The central object is the relaxation function $\\Psi_0(t)/\\Psi_0(0)$ of the underdamped Brownian motion, computed by linear response theory from the Langevin equation, together with the associated relaxation time $\\tau_R=(\\gamma^2+\\omega_0^2)/(2\\gamma\\omega_0^2)$. The mechanism is the limit $\\gamma,\\omega_0\\to 0^+$ at equal rates: this makes the relaxation time diverge while the relaxation function freezes at a constant, putting the dynamics in a critical-like frozen state. For the scaling, the key expression is Eq. (12), the linear-response formula for the impulse-part irreversible work, evaluated against the freezing point $\\hat{t}=f(\\tau)$ defined by Eq. (11).","core_discovery":"In the simultaneous limits $\\gamma\\to 0^+$ and $\\omega_0\\to 0^+$ with $\\gamma/\\omega_0=1$, the underdamped Langevin equation (1) has a relaxation time $\\tau_R=(\\gamma^2+\\omega_0^2)/(2\\gamma\\omega_0^2)$ that diverges (Eq. 6), while the normalized relaxation function $\\Psi_0(t)/\\Psi_0(0)$ tends to the constant 3 (Eq. 7) and the initial susceptibility $\\Psi_0(0)$ diverges (Eq. 8). These are the signatures of Kibble-Zurek behavior: the system freezes because the initial equilibrium temperature must scale as $T\\propto 1/\\gamma$, implying large thermal fluctuations that prevent equilibration. For a linear driving $\\omega(t)=\\omega_0-\\delta\\omega_0\\,t/\\tau$, the irreversible work in the impulse part (Eq. 12) scales as $\\tau^{-2}$ for sudden processes and $\\tau^{-1}$ for slowly varying processes, matching the known quantum-chain exponent in the second case and giving a new exponent in the first. The optimal protocol then reduces to a pause at half the driving with two jumps, which the author attributes to the generic Kibble-Zurek limit rather than to symmetry about a critical point.","pith_inferences":["A natural extension not pursued in the paper is that the same equal-rate vanishing limit may transfer the $\\eta_{KZ}=-2$ sudden-process exponent to other underdamped linear models, such as trapped ions or levitated nanoparticles, where both $\\gamma$ and $\\omega_0$ can be tuned toward zero simultaneously.","If the scaling survives beyond linear response, a direct integration of the exact Langevin dynamics would show a collapse of $W_{\\rm irr}^I(\\tau)$ onto the predicted power laws; checking that collapse would separate the linear-response extrapolation from genuine Kibble-Zurek physics.","The pausing optimal protocol suggests a practical control strategy: to avoid large irreversible work, one should hold the trap frequency fixed in the middle of a fast drive rather than sweeping it through the whole interval.","The paper fixes the ratio $\\gamma/\\omega_0=1$; exploring nearby ratios could reveal a crossover between exponents or a family of scaling curves, a question the paper leaves open."],"forward_implications":["If the claimed limit is realized, the driven underdamped Brownian particle becomes a minimal classical system in which the relaxation time diverges, so the Kibble-Zurek mechanism is not restricted to quantum chains at zero temperature.","For sudden switching, the irreversible work in the impulse part decays as $\\tau^{-2}$, a new scaling distinct from the usual $\\tau^{-1}$ of the transverse-field Ising chain; the author attributes the difference to the different nature of the relaxation time.","For slowly varying processes, the scaling exponent is $\\eta_{KZ}=-1$, the same as in the quantum Ising chain.","The optimal protocol that minimizes average work and its fluctuations becomes approximately constant at $1/2$, with two jumps, meaning the driving pauses in the middle; the author notes this is a coincidence with the critical-point pause, since the driving here is not symmetric about a critical point.","The initial equilibrium temperature must be high ($T\\propto 1/\\gamma$), so the frozen state is associated with strong thermal fluctuations rather than low-temperature criticality."],"supporting_citations":[{"why":"Supplies the underdamped Langevin equation and white-noise statistics used to write Eq. (1).","marker":"[1]"},{"why":"Provides the linear-response formalism and fluctuation-dissipation theorem that produce the relaxation function and the $T\\propto 1/\\gamma$ condition.","marker":"[2]"},{"why":"Defines the Kibble-Zurek scaling baseline (transverse-field Ising chain) that the slowly-varying exponent $\\eta_{KZ}=-1$ is compared with.","marker":"[3]"},{"why":"The previous study cited to justify that the linear-response exponents remain correct even outside the formal range of validity.","marker":"[4]"},{"why":"Gives the optimal-protocol solution $g^*(t)\\approx 1/2$ used for the pausing effect.","marker":"[5]"},{"why":"Supplies the comparison that in the traditional Kibble-Zurek driving the pause happens at the symmetric critical point, which the paper calls a coincidence.","marker":"[6]"}],"fun_headline_variants":["Quantum-chain criticality from a single driven oscillator","Divergent relaxation time in driven Brownian motion","New Kibble-Zurek scaling for sudden processes","Optimal pause emerges in Brownian Kibble-Zurek limit"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"Everything rests on the assumption that the linear-response impulse-work formula, applied outside its formal range of validity (the paper states that $f(\\tau)$ is always greater than the $\\tau$ range used), still gives the correct Kibble-Zurek exponents because an earlier calculation with the same approximation worked.","fun_headline_variants_meta":{"raw":{"variants":["Quantum-chain criticality from a single driven oscillator","Divergent relaxation time in driven Brownian motion","New Kibble-Zurek scaling for sudden processes","Optimal pause emerges in Brownian Kibble-Zurek limit"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000322,"raw_usage":{"total_tokens":1805,"prompt_tokens":932,"completion_tokens":873,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":548,"completion_tokens_details":{"reasoning_tokens":809}},"tokens_in":548,"tokens_out":873,"duration_ms":8228,"temperature":1.0,"reasoning_tokens":809,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T04:23:08.103543+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Directly integrate the driven Langevin equation (1) for $\\gamma=\\omega_0=0.01$, $\\delta\\omega_0=0.001$, over switching times $\\tau$ in the same range as Figures 1 and 2, computing the exact irreversible work from the full work distribution without the linear-response approximation, and compare the exponents of $W_{\\rm irr}^I(\\tau)$; if they are not $-2$ (sudden) and $-1$ (slowly varying), the central scaling claim fails.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the underdamped Langevin equation and white-noise statistics used to write Eq. (1)."},{"cited_title":"Gomez-Marin, T","cited_arxiv_id":null,"evidence_quote":"Provides the linear-response formalism and fluctuation-dissipation theorem that produce the relaxation function and the $T\\propto 1/\\gamma$ condition."},{"cited_title":"Deffner, Physical Review E 96, 052125 (2017)","cited_arxiv_id":null,"evidence_quote":"The previous study cited to justify that the linear-response exponents remain correct even outside the formal range of validity."},{"cited_title":"Naz´ e, M","cited_arxiv_id":null,"evidence_quote":"Gives the optimal-protocol solution $g^*(t)\\approx 1/2$ used for the pausing effect."},{"cited_title":"Naz´ e, Journal of Statistical Mechanics: Theory and E x- periment 2024, 073205 (2024)","cited_arxiv_id":null,"evidence_quote":"Supplies the comparison that in the traditional Kibble-Zurek driving the pause happens at the symmetric critical point, which the paper calls a coincidence."}],"review_version":1}