{"id":"4bb949ef-2c06-44dd-a2e2-fe305d9c22d4","arxiv_id":"2505.10256","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"For a harmonic oscillator chain with exchange noise, the diffusive limit is two independent heat equations for weak springs, a coupled parabolic system at the critical strength, and volume fluctuations satisfy an Edwards-Wilkinson-style identity.","lead":"A noisy chain of coupled springs is shown to behave, at large scales, like either plain heat diffusion or a coupled nonlinear diffusion system, depending on how strongly the springs push. The paper also proves precise statistical properties of the density fluctuations around this large-scale behavior.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Theorem 2.11 is proved for bounded measurable initial profiles, but its proof of item (b) invokes the C∞-profile energy estimate (2.14); the needed energy hydrodynamics is not available under the stated assumptions.","rationale":"I read the paper in good faith and find the core hydrodynamic results plausible: Theorems 2.5, 2.7, and 2.8 have a coherent entropy/block-structure architecture, the fourth-moment bound in Theorem 2.13 is a genuine technical contribution, and the correlation estimate in Theorem 2.14 supplies the main quantitative input. My stress-test does not identify a flaw in the random-walk local-time argument itself: the apparent boundary-Laplacian notation in Lemma 6.2 is consistent with a second difference along the boundary line, and the final absorption for κ≥1 is sound. The load-bearing concern is narrower but directly hits the advertised fluctuation characterization: Theorem 2.11 part (b) is the only mechanism that identifies the limit martingale's quadratic variation, yet its proof is not valid under the stated hypotheses because it imports Theorem 2.8's C∞-smooth, product-measure energy-profile estimate. For κ=1, the energy hydrodynamics for general bounded measurable e0 is explicitly not obtained, so the expected quadratic variation formula (2.16) is unsupported for the full statement of Theorem 2.11. This does not overturn the paper: it suggests the main fluctuation theorem should be restricted to smooth product-type initial data, or the authors should supply a proof of (2.16) under Assumption 2.1 alone. The reader's weakest assumption (initial two-point correlation decay) is a stated and essential hypothesis, and I agree it limits scope, but the unstated smoothness gap in the proof of part (b) is more immediately load-bearing for the central claim. I therefore recommend keeping the CONDITIONAL verdict, with the condition being that the hypotheses of Theorem 2.11 be reconciled with the energy-profile estimate actually used.","tokens_in":37523,"tokens_out":34306,"duration_ms":331436,"concrete_test":"Attempt to re-prove the convergence in (2.16) for a discontinuous initial energy profile, e.g. e0(u) = 1 + 1_{[0,1/2)}(u), with v0 smooth and µN the corresponding product Gaussian measure satisfying Assumption 2.1. In the proof of Section 4.4, replace the use of (2.14) by only the convergence statements that are actually established under Assumption 2.1 (Theorems 2.5 and 2.7, the latter leaving energy in probability open for κ=1). If the weighted energy sum (1/N)Σ_x E[e^N_s(x)+e^N_s(x+1)-v^N_s(x)^2-v^N_s(x+1)^2] (T^-_{f_N s} ∇G)^2 cannot be shown to converge to ∫ 2χ(s,u)(∇G)^2 du without the C∞ product-measure hypotheses of Theorem 2.8, then Theorem 2.11 must either be restricted to those hypotheses or supplied with a new energy-profile estimate for bounded measurable e0.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central fluctuation result, Theorem 2.11, is stated for bounded measurable v0 and e0 with e0-v0^2 ≥ 0, under Assumption 2.1. The proof of part (b) in Section 4.4 computes the limit of the expected quadratic variation of the volume fluctuation martingale and needs to replace the microscopic energy profile e^N_s by the macroscopic e(s, ·) in the combination e^N_s(x)+e^N_s(x+1)-v^N_s(x)^2-v^N_s(x+1)^2. The text does this 'using (2.14)', i.e. using the quantitative bound sup_{0≤t≤T} max_x |e^N_t(x)-e(t,x/N)| ≤ C log N / N. But (2.14) is Theorem 2.8, whose hypotheses include v0,e0 ∈ C^∞_b and an initial product measure of the form (2.9); these are not part of Theorem 2.11. Moreover, for κ=1 the energy hydrodynamic limit in probability is explicitly left open in Theorem 2.7, and the only energy-profile statement, Theorem 2.8, is conditional on smoothness and product structure. Thus, for the generality claimed in Theorem 2.11, the convergence needed for (2.16) is not proven. If Theorem 2.11 is intended only for smooth product initial measures, that restriction is missing; if it is intended for general initial measures, the proof has a genuine gap in the derivation of part (b), which is the only quantitative handle on the limit martingale.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the Bernardin-Stoltz model of a harmonic chain perturbed by an exchange noise, with the Hamiltonian strength scaled by α_N = α N^{-κ} and in the diffusive time scale N^2. The main results are: (i) for κ>1, a hydrodynamic limit in probability for both volume and energy empirical measures, converging to two autonomous heat equations (Theorem 2.5); (ii) for κ=1, a hydrodynamic limit in probability for the volume only, with a drift term (Theorem 2.7), and a quantitative L∞-rate estimate for the energy profile under stronger smooth/product initial data (Theorem 2.8); (iii) a non-equilibrium fluctuation result for the volume field, tightness and partial characterization of the limit martingale (Theorem 2.11). The proofs are based on entropy-method tightness, Dynkin martingales, a uniform fourth-moment bound (Theorem 2.13), and a decay estimate for the two-point volume correlation (Theorem 2.14), the latter being the main technical ingredient.","tokens_in":37818,"tokens_out":9456,"duration_ms":88515,"significance":"If the stated results are correct, the paper would extend rigorous hydrodynamic limit and non-equilibrium fluctuation results for the Bernardin-Stoltz model in a regime where the two conserved quantities do not evolve autonomously, and it would provide quantitative energy-profile estimates. The proof strategy is clearly organized and uses standard tools (entropy method, Duhamel representations, random-walk estimates), with explicitly stated assumptions and no fitted parameters. The paper also honestly identifies several open problems, such as the lack of a full characterization of the limit martingale and of energy correlation decay. These strengths make the work potentially useful to the stochastic-particle-systems community. However, two load-bearing technical issues in the correlation estimate and in the fluctuation theorem need to be resolved before the central claims can be accepted.","major_comments":[{"comment":"The estimate (6.4) cannot hold as stated for κ>1. For a non-constant smooth initial profile v0, the term -(∇_N v^N_0(x))^2 in the definition of g^N_0(x), Eq. (6.1), is generically of order 1, while α_N^2 N^2 = α^2 N^{2-2κ} tends to 0. Even if the α_N N^2 h^N_0 term is small, this leaves the left-hand side of (6.4) of order 1 and the right-hand side tending to 0, a contradiction. The proof of Lemma 6.2 does not repair this: from the bound on h^N in the proof one obtains ‖g^N‖ bounded by C + Cα_N N + Cα_N^2 N^{8/3}‖φ^N‖, not by (6.4). Since (6.4) is the input to the bootstrap (6.5) proving Theorem 2.14, the proof of the key correlation decay estimate is invalid as written. A corrected bound of the form ‖g^N‖ ≤ C(1 + N^{2/3}‖φ^N‖) might still imply Theorem 2.14, so the result may be salvageable, but the present argument does not establish it.","section":"Section 6, Lemma 6.2, Eq. (6.4)"},{"comment":"The proof of item (b) explicitly invokes (2.14), i.e. the quantitative energy-profile estimate of Theorem 2.8, to replace e^N_s by e(s,·) in the quadratic variation. However Theorem 2.8 requires v0,e0 ∈ C∞_b, a product initial measure of the form (2.9), and initial data (2.13); none of these assumptions appear in Theorem 2.11, whose stated hypotheses are only Assumption 2.1, bounded measurable profiles, and convergence of the initial fluctuation field. Moreover, Theorem 2.7 explicitly leaves the energy hydrodynamic limit open for κ=1 in the general setting. Consequently, the convergence needed for (2.16) is not proved under the hypotheses of Theorem 2.11. The statement must either restrict Theorem 2.11 to the smooth product setting of Theorem 2.8 or supply a new energy-profile estimate valid under Assumption 2.1.","section":"Section 4.4, proof of Theorem 2.11(b), Eq. (2.16)"},{"comment":"Related to the previous comment, the bootstrap in (6.5) also relies on the specific form of (6.4). If (6.4) is replaced by the correct order-1 bound for ‖g^N‖, the factor (1/N)‖g^N‖ would give a contribution of order 1/N, which is still compatible with the desired conclusion; but the manuscript does not provide that argument. As written, the proof of Theorem 2.14 is incomplete, and Theorem 2.14 is used in the proofs of Theorems 2.5, 2.7, 2.8, and 2.11.","section":"Section 6, proof of Theorem 2.14"}],"minor_comments":[{"comment":"There is a typo: 'weal solution' should read 'weak solution'.","section":"Theorem 2.5"},{"comment":"The proof refers to 'Theorem 3.5' when the intended statement is Lemma 3.5.","section":"Corollary 3.6"},{"comment":"The phrase 'general initial measures' overstates the hypotheses: Assumption 2.1 requires the initial two-point volume correlation to be O(1/N) pointwise and imposes uniform bounds on discrete gradients of the initial profiles; Example 2.2 exhibits only a subclass of such measures. The wording should be adjusted to match the stated assumptions.","section":"Abstract and Introduction"},{"comment":"The text refers to 'Theorem 6.1' and 'Theorem 6.2' where the objects are labelled as Lemma 6.1 and Lemma 6.2; this should be corrected for consistency.","section":"Section 6"},{"comment":"In the Aldous-criterion estimate, the displayed probability incorrectly mixes the volume martingale and the energy martingale: it should read |M^{e,N}_{τ+θ}(f)-M^{e,N}_τ(f)|, not |M^{v,N}_{τ+θ}(f)-M^{e,N}_τ(f)|.","section":"Lemma 3.11 proof"}],"recommendation":"major_revision","confidential_remarks":"The paper contains substantial original work and the overall research program is relevant for the journal. The two major issues identified are load-bearing: the false/inapplicable estimate in Lemma 6.2 and the unjustified use of Theorem 2.8 in the proof of Theorem 2.11(b). Both appear addressable by either adding hypotheses or developing additional estimates, so I recommend major revision rather than rejection. I saw no issues with novelty disclosure or citation practices."},"author_rebuttal":null,"desk_editor":null,"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["60H15","60J27","82C22","60K50"],"pacs":[],"model":"deepseek-v4-flash","headline":"For a stochastic harmonic oscillator chain, this paper proves two autonomous heat equations when the Hamiltonian part decays faster than 1/N, a coupled system at the critical rate, and a martingale characterization of volume fluctuations.","keywords":["Bernardin–Stoltz model","harmonic oscillators","hydrodynamic limit","non-equilibrium fluctuations","Ornstein–Uhlenbeck process","Edwards–Wilkinson universality","two-point correlation decay","fourth-moment estimates"],"falsifier":"Start the $\\kappa=1$ dynamics from a smooth profile with nonzero $\\nabla v_0$ and measure the energy profile: the coupled equation predicts a visible $2\\alpha\\nabla v^2$ source term, so its absence would disprove Theorem 2.8. Independently, begin from a product Gaussian with $O(1)$ long-range covariance: Theorem 2.14 predicts $\\sup_{0\\le t\\le T}\\|\\varphi^N_t\\|_{\\ell^\\infty}\\le C/N$, so a numerical violation of that decay would falsify the correlation estimate and the tightness arguments built on it.","tokens_in":37304,"feed_emoji":"⚛️","tokens_out":9430,"duration_ms":87456,"temperature":0.7,"pith_summary":"The paper studies the Bernardin–Stoltz model, a harmonic chain of oscillators with random nearest-neighbor exchanges, under diffusive scaling and with Hamiltonian strength $\\alpha_N=\\alpha N^{-\\kappa}$. It proves that when $\\kappa>1$ the two conserved macroscopic fields—volume and energy—each converge to the solution of a plain heat equation, $\\partial_t v=\\Delta v$ and $\\partial_t e=\\Delta e$, so the conservation laws decouple. At the critical value $\\kappa=1$ the limiting system is coupled: $\\partial_t v=\\Delta v+2\\alpha\\nabla v$ and $\\partial_t e=\\Delta e+2\\alpha\\nabla v^2$, with volume convergence in probability and energy convergence in expectation, plus an $O(\\log N/N)$ bound for smooth profiles. It also characterizes non-equilibrium volume fluctuations for general initial states satisfying a short-range correlation condition: limit points have the form $V_t(f)=V_0(T_t f)+M_t(f)$, where $M_t$ is a mean-zero martingale whose mean quadratic variation is $\\int_0^t\\int 2\\chi(s,u)(\\nabla f)^2\\,du\\,ds$, with $\\chi=e-v^2$. A sympathetic reader would care because this gives a rigorous, parameter-dependent derivation of macroscopic heat and drift equations from a microscopic unbounded-state-space dynamics, together with the correlation estimates needed to control fluctuations out of equilibrium.","feed_headline":"Weak Hamiltonian dynamics yields two decoupled heat equations","feed_subtitle":"At the critical Hamiltonian strength the equations couple, yet volume fluctuations retain a martingale-plus-heat limit.","key_machinery":"The argument is carried by the two-point volume correlation function $\\varphi^N_t(x,y)=E[\\eta_t(x)\\eta_t(y)]-v^N_t(x)v^N_t(y)$, defined off the diagonal of $\\mathbb{T}_N^2$. Its evolution equation is governed by a two-dimensional random walk that is reflected on the diagonal lines $y=x\\pm1$; via Duhamel's principle the correlation is written as an expectation of that walk, and the key estimate bounds the walk's local time at the diagonal by $O(1/N)$, yielding $\\sup_t\\|\\varphi^N_t\\|_{\\ell^\\infty}\\le C/N$ for $\\kappa\\ge1$. A second load-bearing object is the discrete $H^{-1}$ norm $\\|\\eta^2\\|_{-1,N}$, whose generator computation produces the uniform fourth-moment bound $\\int_0^T N^{-1}\\sum_x E[\\eta_s(x)^4]\\,ds\\le C(1+\\alpha_N N)$. Together these estimates replace the one- and two-block replacement steps that are unavailable for unbounded real-valued variables, and they power both the hydrodynamic-limit proofs and the fluctuation field.","core_discovery":"The central discovery is a decoupling–coupling transition controlled by $\\kappa$, the exponent in $\\alpha_N=\\alpha N^{-\\kappa}$. For $\\kappa>1$, Theorem 2.5 shows the volume and energy empirical measures converge in probability to $v$ and $e$ solving the two autonomous heat equations. For $\\kappa=1$, Theorem 2.7 shows the volume field converges in probability to $\\partial_t v=\\Delta v+2\\alpha\\nabla v$, while Theorem 2.8 shows the energy profile converges at rate $\\sup_t\\max_x |e^N_t(x)-e(t,x/N)|\\le C\\log N/N$ to $\\partial_t e=\\Delta e+2\\alpha\\nabla v^2$. For fluctuations, Theorem 2.11 proves tightness of the volume fluctuation field and identifies every limit point as $V_t(f)=V_0(T_t f)+M_t(f)$ with a mean-zero martingale $M_t$ of mean quadratic variation $\\int_0^t\\int 2\\chi(s,u)(\\nabla f)^2\\,du\\,ds$; from invariant Gaussian product measures the limit is the Ornstein–Uhlenbeck equation $dV_t=\\Delta V_t\\,dt+\\sqrt{2/\\beta}\\,\\nabla\\dot W_t$.","pith_inferences":["If the correlation decay extends to four-point functions, the volume fluctuation limit should be fully characterized rather than merely tight, giving the noise term in the volume equation; the paper explicitly leaves this as an open problem.","The same machinery suggests that for $\\kappa<1$, where the proof's absorption step fails, the volume–energy coupling may produce superdiffusive or anomalous transport; that is an implicit prediction, not a theorem of this paper.","One could test the role of the initial-state assumption by running the dynamics from a correlated Gaussian initial measure with $O(1)$ long-range covariances: the predicted $C/N$ decay of $\\varphi^N_t$ should fail, separating the contribution of the initial condition from the dynamics.","The $O(\\log N/N)$ energy-profile rate is likely not optimal; a sharper estimate would require better short-time control of the diagonal local time and could be checked numerically against the actual sup-norm discrepancy."],"forward_implications":["For $\\kappa>1$, volume and energy evolve macroscopically as two independent heat equations; no coupling term survives in the limit.","At $\\kappa=1$, energy transport is slaved to volume: a nonzero volume gradient produces a source term $2\\alpha\\nabla v^2$ in the energy equation.","Out-of-equilibrium volume fluctuations are universal in form for the allowed initial measures: every limit point is a heat-evolved initial field plus a mean-zero martingale, with the noise intensity fixed by the local static compressibility $\\chi=e-v^2$.","Two-point volume correlations decay as $C/N$ uniformly up to a fixed time for $\\kappa\\ge1$, so second-order structure is tractable even though the state space is unbounded.","For smooth data at $\\kappa=1$, the discrete energy profile is within $O(\\log N/N)$ of the solution of the coupled PDE uniformly in space and time."],"supporting_citations":[{"why":"Introduces the Bernardin–Stoltz oscillator model with exchange noise and its two conserved quantities; this is the dynamics under study.","marker":"[8]"},{"why":"Derives the hydrodynamic limit for harmonic oscillators with continuous noise, giving the earlier heat-equation benchmark this paper extends to exchange-type noise.","marker":"[2]"},{"why":"Supplies the entropy-method toolbox—Sobolev-space tightness, Aldous criterion, and $H^{-1}$ estimates—used throughout the hydrodynamic and fluctuation proofs.","marker":"[21]"},{"why":"Establishes equilibrium volume fluctuations of the harmonic BS model as an Ornstein–Uhlenbeck process, the stationary case recovered in Theorem 2.11(d).","marker":"[4]"},{"why":"Studies the BS model under generic potentials and derives coupled or anomalous behavior, providing the background for interpreting the $\\kappa=1$ coupling.","marker":"[13]"},{"why":"Derives coupled SPDEs from the BS model, showing how volume and energy can couple beyond the purely harmonic case and motivating the coupled limiting system.","marker":"[1]"}],"fun_headline_variants":["Stochastic oscillator chain: decoupling transition in heat equations","Noise and Hamiltonian strength decide heat equation coupling","Volume fluctuations: martingale limit with heat kernel","Oscillator chain scaling: from decoupled to coupled parabolic systems","Heat decoupling for strong driving in stochastic oscillator chains"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The proof needs the initial two-point volume correlations to decay as $C/N$, together with an initial $H^{-1}$ bound on the squared configuration; an initial state with order-one long-range correlations is not covered, and the correlation estimate and fluctuation characterization would likely break down.","fun_headline_variants_meta":{"raw":{"variants":["Stochastic oscillator chain: decoupling transition in heat equations","Noise and Hamiltonian strength decide heat equation coupling","Volume fluctuations: martingale limit with heat kernel","Oscillator chain scaling: from decoupled to coupled parabolic systems","Heat decoupling for strong driving in stochastic oscillator chains"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000706,"raw_usage":{"total_tokens":3160,"prompt_tokens":904,"completion_tokens":2256,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":520,"completion_tokens_details":{"reasoning_tokens":2188}},"tokens_in":520,"tokens_out":2256,"duration_ms":16983,"temperature":1.0,"reasoning_tokens":2188,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T21:14:35.541097+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Start the $\\kappa=1$ dynamics from a smooth profile with nonzero $\\nabla v_0$ and measure the energy profile: the coupled equation predicts a visible $2\\alpha\\nabla v^2$ source term, so its absence would disprove Theorem 2.8. Independently, begin from a product Gaussian with $O(1)$ long-range covariance: Theorem 2.14 predicts $\\sup_{0\\le t\\le T}\\|\\varphi^N_t\\|_{\\ell^\\infty}\\le C/N$, so a numerical violation of that decay would falsify the correlation estimate and the tightness arguments built on it.","supporting_citations":[{"cited_title":"Bernardin","cited_arxiv_id":null,"evidence_quote":"Derives the hydrodynamic limit for harmonic oscillators with continuous noise, giving the earlier heat-equation benchmark this paper extends to exchange-type noise."},{"cited_title":"Kipnis and C","cited_arxiv_id":null,"evidence_quote":"Supplies the entropy-method toolbox—Sobolev-space tightness, Aldous criterion, and $H^{-1}$ estimates—used throughout the hydrodynamic and fluctuation proofs."},{"cited_title":"Bernardin, P","cited_arxiv_id":null,"evidence_quote":"Establishes equilibrium volume fluctuations of the harmonic BS model as an Ornstein–Uhlenbeck process, the stationary case recovered in Theorem 2.11(d)."}],"review_version":1}