{"id":"031d336d-fbe2-46bc-81c3-c9897630d02f","arxiv_id":"2607.13229","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":3.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A finite harmonic chain with noise and end thermostats converges, under n^{3/2} time scaling, to a fractional heat equation with a Neumann fractional Laplacian and nonlocal boundary terms.","lead":"The paper derives, for a one-dimensional chain of atoms with random momentum exchanges and heat baths at both ends, the correct mathematical boundary conditions for superdiffusive heat flow. It shows the average energy profile converges to a fractional heat equation with a special nonlocal boundary correction.","discovery_kind":"review","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Boundary covariance closure (Props 4.4–4.5) is cited to companion [14], not proved here; the nonlocal boundary terms and c_bd depend entirely on that deferred argument.","rationale":"The paper's central claim is the rigorous derivation of macroscopic boundary conditions for superdiffusive heat transport, specifically the convergence to the fractional heat equation with a Neumann fractional Laplacian and nonlocal boundary terms. The proof outline is coherent: the bulk term asymptotics (Prop 4.1) are sketched in detail, and the stochastic current contribution is shown to vanish. The critical step is the boundary covariance closure, which converts the open equation (4.40) into the closed fractional equation (2.7). The reader's weakest_assumption identified exactly this, and I agree. Notably, the sub-assumption T*T=pi^2 I, highlighted by the reader, is actually correct: T equals pi times the isometry RHE (restriction of the Hilbert transform of the even extension), so T*T=pi^2 I follows from the isometry property. Thus the real concern is not that identity but the unproved Propositions 4.4 and 4.5, which are deferred to the companion paper [14]. This is not an internal inconsistency but a verification gap: the advertised theorem is not self-contained, and the crucial closure is untested here. The paper is honest about this, stating it is a review of [14], and no ad hominem is warranted. Given the reader's CONDITIONAL verdict, my read does not change it, so verdict_should_be is UNCHANGED. The concrete test I propose—a direct numerical simulation of the microscopic model plus an independent re-derivation of the deferred propositions—would settle whether the concern lands.","tokens_in":18812,"tokens_out":14345,"duration_ms":126955,"concrete_test":"Independently verify the two deferred propositions: (1) Re-derive [14, Props 9.1–9.2] from the stationary covariance equations of the half-line limit, checking that (4.42) and (4.45) hold with the stated constants. (2) As a numerical falsification, simulate the stochastic harmonic chain (1.2) with n=1024, gamma=gamma_tilde=1, T_L=1, T_R=2, and initial local equilibrium at T_ini=1; at time t=2*n^{3/2}, compare the measured energy profile to the solution of (2.7) with coefficients (2.17). A deviation significantly larger than the expected n^{-1/2} fluctuations would indicate the closure is wrong.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central theorem (2.6) reduces to closing the evolution equation (4.40) for the energy profile. Closure is achieved through Propositions 4.4 and 4.5, which express the boundary-bulk covariance fields b(p,v) and b(pr,v) in terms of the bulk energy and bath temperatures via the operator T. Neither proposition is proved in this manuscript; both are cited to [14, Props 9.1–9.2, Thm C.1]. The combination of these propositions (together with T*T=pi^2 I) is what produces the nonlocal boundary terms in (2.7) and fixes the coefficient c_bd = sqrt(2 gamma_tilde) c_bulk/(1+gamma_tilde)^2. If the L^2 bounds (4.37)–(4.39) only give weak convergence and the asserted pointwise relations for b(p,v) and b(pr,v) fail, the closed equation would not be (2.7) and the advertised boundary conditions would be incorrect. Since the manuscript is explicitly a review of [14], the correctness of the advertised result is untestable from this paper alone. The operator identity T*T=pi^2 I itself is consistent with an isometry (T=pi times the Hilbert transform composed with even extension), so the weak point is not that identity but the unshown asymptotic relations.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper studies a finite one-dimensional unpinned harmonic chain with stochastic nearest-neighbor momentum exchange and Langevin heat baths at the endpoints. Under the superdiffusive scaling t ∼ n^{3/2}, the authors claim that the averaged microscopic energy profile converges to the solution of a fractional heat equation on [0,1] whose generator is a Neumann fractional Laplacian plus nonlocal boundary terms, with explicit coefficients c_bulk = (2^3 γ)^{-1/2} and c_bd = √(2 γ̃) c_bulk/(1+γ̃)^2. The argument is presented as an outline: the bulk term is computed in detail in §4.4, while the boundary closure, which converts the open evolution equation (4.40) into the closed equation (2.7), is deferred to the companion paper [14]. The paper also states entropy, covariance, and current estimates; several of these are proved only in special cases or cited to [14].","tokens_in":19120,"tokens_out":9169,"duration_ms":76824,"significance":"If the main theorem holds, this is a significant contribution: it provides a rigorous derivation of macroscopic boundary conditions for superdiffusive energy transport in open chains from a microscopic dynamics, and it identifies concrete nonlocal boundary terms for fractional Laplacians that appear to be new. The paper has genuine strengths: the bulk computation in §4.4 is spelled out, the coefficients are derived from microscopic parameters with no fitted constants, and the probabilistic interpretation in Remark 2.5 is illuminating. No circularity or hidden free parameters are apparent. However, the central boundary covariance closure is not proved in this manuscript, so the advertised result is conditional on the companion paper.","major_comments":[{"comment":"These results are the load-bearing step. Equation (4.40) contains the boundary-bulk covariance fields b_n^{(p,v)} and b_n^{(pr,v)}. Propositions 4.4 and 4.5 assert asymptotic relations expressing these fields in terms of the bath temperatures and the bulk energy, and Theorem 4.6 combines them (using T*T = π² I) to yield the formula that, inserted into (4.40), produces the nonlocal boundary terms and the coefficient c_bd in (2.17). All three are cited to [14, Props 9.1–9.2, Thm C.1], with no proof or proof scheme in this manuscript. Since the main theorem (2.6) reduces exactly to this closure, the present text does not by itself establish the advertised result. The authors acknowledge this: Section 2.2 states 'The detailed argument is given in [14]'.","section":"Section 4.5, Propositions 4.4–4.5 and Theorem 4.6"},{"comment":"The abstract claims 'we prove that the averaged microscopic energy profile converges...', but the proof of Theorem 2.6 is an outline and essential ingredients are deferred to [14]: Theorem 2.2 (well-posedness of the limiting equation), Theorem 2.8 (current bound), the general case of Theorem 3.1, and Theorem 3.3. In addition, Proposition 4.3, used to derive the boundary term in (4.40), is stated without proof or a reference. A reader cannot verify the core claim from this paper alone. The paper should either be reframed explicitly as an announcement/review of [14], with 'prove' replaced by 'announce', or the omitted arguments should be included.","section":"Section 2.2 / Theorem 2.6"},{"comment":"Proposition 4.3 gives the asymptotic formula for the boundary term π_p^{(pr)} that enters the open equation (4.40). No proof is given, and no citation to [14] appears at this point. Since the formula fixes the prefactor 2^{1/2} γ̃ / π and the kernel 1/((πℓ)² + γ² ϱ⁴), it is load-bearing for the boundary coefficient. Please add a proof or a precise reference to the corresponding statement in [14].","section":"Section 4.5, Proposition 4.3"}],"minor_comments":[{"comment":"The Introduction says 'In this paper we review our recent result [14]', while the Abstract claims 'we prove'. This inconsistency should be resolved; the paper's role (research announcement versus self-contained proof) must be stated consistently.","section":"Introduction / Abstract"},{"comment":"The displayed definition reads 'where φ̂'_s(j) = √2 ∫_0^1 φ(u) sin(πju) du', but the subsequent use in (4.36), 'φ̂'_s(ℓ) = -πℓ φ̂_c(ℓ)', indicates the sine transform of φ', not of φ. The definition likely should involve φ'. Please fix the notation.","section":"Equation (4.28)"},{"comment":"The article number '0650002' appears to be a typo; the standard J. Stat. Mech. formatting would be '065002'.","section":"Reference [25]"},{"comment":"The parity argument used to conclude that θ_{pr}^{(e)} = 0 is only sketched. A sentence explaining the symmetry or a pointer to [14] would help the reader.","section":"Section 4.4"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is essentially an announcement of the authors' companion paper [14], and several core results (Props 4.4–4.5, Thm 4.6, Thm 2.2, Thm 2.8) are deferred to that preprint. The editor should verify the status of [14] before deciding. If [14] is not yet accepted, acceptance of the present paper is risky. The abstract's 'prove' claim is too strong for the content as written."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is an outline of the authors' own result from [14], not a new result standing alone. The introduction says so plainly. Its value is that it puts the macroscopic equation (2.7) with the nonlocal bath terms and the explicit coefficients c_bulk and c_bd in one place, and it gives a readable path from the microscopic model to that equation. The bulk part of the argument (Prop 4.1) is actually spelled out in Section 4.4, including the discrete Fourier analysis and the constant (2^3 gamma)^(-1/2). That is a concrete, checkable computation, and it looks right. The entropy and covariance bounds are standard and cleanly presented. The soft spot is exactly where the stress test lands: the closure of the boundary terms. Propositions 4.4 and 4.5, plus Theorem 4.6, are cited to [14] and not proved here. Everything rests on those: they convert the open equation (4.40) into the closed equation (2.7) and they deliver the coefficient c_bd. The operator identity T*T = pi^2 I is stated with a reference; it is consistent with an isometric structure, so I don't doubt that, but the asymptotic relations for b_n are the load-bearing piece. Without seeing those proofs, this manuscript alone does not fully support the advertised conclusion. That is not a fatal flaw because the companion exists and the dependence is stated honestly; it does mean a referee needs to read [14] before signing off. There is nothing fitted here. The coefficients come from gamma and gamma_tilde, and no free parameters appear. The citation to [14] is appropriate, not a dodge. The one thing I'd want in revision is at least a sketch of how Props 4.4 and 4.5 are obtained, or a statement of the main estimate behind them, so the survey stands a little more on its own. Who is this for? Anyone working on fractional hydrodynamics or open harmonic chains who wants a concise statement of the boundary conditions without digging through the full proofs. It deserves a serious referee, but the referee should get [14] as part of the package. If [14] holds up, this is a useful summary; if not, the present paper inherits the problem.","headline":"A clear, self-described outline of the authors' own companion theorem; the decisive boundary-closure step is deferred to [14], so this paper stands or falls with that companion.","tokens_in":750,"tokens_out":877,"would_cite":true,"duration_ms":42886,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35R11","82C70","60K35"],"pacs":[],"model":"deepseek-v4-flash","headline":"A finite unpinned harmonic chain with stochastic momentum exchange and Langevin baths is proved to converge, under superdiffusive scaling, to a fractional heat equation on [0,1] with a Neumann fractional Laplacian and nonlocal boundary term","keywords":["fractional heat equation","superdiffusion","harmonic chain","Langevin baths","nonlocal boundary conditions","Neumann fractional Laplacian","energy transport","hydrodynamic limit"],"falsifier":"Compute numerically the left and right sides of (4.42) in a finite chain for increasing n and several values of γ̃; if the difference fails to vanish as n grows—or if a direct calculation shows T*T≠π²I on L²[0,∞)—the nonlocal boundary terms would not emerge and the limit equation would differ. Alternatively, measure the steady-state temperature profile of a long open chain and compare it with the stationary solution of (2.7); a mismatch in the boundary layer would indicate the closure is wrong.","tokens_in":18694,"feed_emoji":"🔥","tokens_out":4807,"duration_ms":60245,"temperature":0.7,"pith_summary":"The paper establishes the macroscopic boundary conditions for superdiffusive energy transport in an open one-dimensional chain. For a finite unpinned harmonic chain with stochastic nearest-neighbor momentum exchanges and Langevin heat baths at both ends, the authors prove that under the time scaling t∼n^{3/2} the averaged energy profile converges to a temperature field solving a fractional heat equation on [0,1]. The generator is a Neumann fractional Laplacian plus nonlocal boundary terms induced by the baths, with explicit coefficients cbulk=(2^3γ)^{-1/2} and cbd=√(2γ̃)cbulk/(1+γ̃)^2. This provides a rigorous derivation of boundary conditions for fractional heat equations from a microscopic model and predicts boundary layers typical of unpinned chains.","feed_headline":"Two coefficients fix the boundary conditions of superdiffusive heat flow","feed_subtitle":"A rigorous limit of an open harmonic chain yields a fractional heat equation whose bath-induced boundary terms are now explicit.","key_machinery":"The key object is the covariance closure for the boundary-bulk fields b(p,v)_n and b(pr,v)_n, stated in Propositions 4.4–4.5 and Theorem 4.6. These propositions assert that boundary-bulk correlations can be expressed through the bath temperatures and the Fourier cosine coefficients of the bulk energy field, via the integral operator T on L²[0,∞) for which T*T=π²I. This closure converts the non-closed open evolution (4.40) into the closed fractional equation (2.7). The bulk part is carried by the Fourier/cosine decomposition of the fractional Laplacian and by the energy equipartition Proposition 4.2, which together yield the coefficient cbulk.","core_discovery":"The central discovery is that, in the hydrodynamic limit, the open chain's energy profile does not follow a simple Dirichlet or Neumann fractional heat equation. Instead, the end baths generate additional nonlocal terms that enter the limiting equation as boundary absorption and creation rates, enforcing the bath temperatures at the endpoints while producing boundary layers. The proof works by converting the open evolution, which is not closed, into a closed fractional equation for the energy field: boundary-to-bulk covariance fields are expressed in terms of the bath temperatures and the Fourier coefficients of the bulk energy, using an integral operator T satisfying T*T=π²I. The result is","pith_inferences":["Because the proof of the boundary covariance closure is cited to the companion paper rather than carried out here, a direct numerical test of the identity (4.42) at moderate system sizes would show whether the predicted nonlocal boundary terms are robust or only asymptotic.","The same nonlocal boundary structure may transfer to other low-dimensional nonlinear chains, where simulations already display superdiffusion; if so, the ratio cbd/cbulk provides a testable signature of how the bath coupling mechanism changes the boundary layer.","The creation-and-annihilation reading of the limiting equation suggests constructing a Markov jump process on [0,1] with those rates; its invariant measure should reproduce the steady-state temperature profile predicted by the fractional equation, offering an independent check."],"forward_implications":["If the theorem is correct, open superdiffusive heat transport is fully characterized macroscopically by a well-posed fractional heat equation with the bath temperatures imposed at the endpoints.","The bath-induced boundary conditions are nonlocal, not classical Dirichlet or Neumann conditions, and they produce boundary layers characteristic of unpinned chains, with a natural interpretation in terms of creation and annihilation rates for a Lévy-type process.","The ratio cbd/cbulk is maximal at γ̃=1 and vanishes as γ̃→0 or γ̃→∞, predicting how the relative strength of the baths modifies the boundary layer.","The energy currents are bounded uniformly in time by C/√n, so the macroscopic limit has controlled fluctuations and a well-defined temperature in the bulk via equipartition."],"fun_headline_variants":["Superdiffusive heat flow gets explicit boundary terms","Fractional heat equation gains physical boundary conditions","Open chain limit reveals nonlocal bath effects on heat","Bath-induced boundary layers in superdiffusive transport","Rigorous limits fix fractional Laplacian boundary terms"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The whole boundary picture rests on the covariance closure in Propositions 4.4–4.5 and Theorem 4.6: that boundary-bulk correlations can be exactly expressed through bath temperatures and the bulk energy field via the operator T with T*T=π²I, a proof that is only cited to the companion paper, not given here.","fun_headline_variants_meta":{"raw":{"variants":["Superdiffusive heat flow gets explicit boundary terms","Fractional heat equation gains physical boundary conditions","Open chain limit reveals nonlocal bath effects on heat","Bath-induced boundary layers in superdiffusive transport","Rigorous limits fix fractional Laplacian boundary terms"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000164,"raw_usage":{"total_tokens":1064,"prompt_tokens":706,"completion_tokens":358,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":450,"completion_tokens_details":{"reasoning_tokens":284}},"tokens_in":450,"tokens_out":358,"duration_ms":4362,"temperature":1.0,"reasoning_tokens":284,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-02T05:47:48.615366+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute numerically the left and right sides of (4.42) in a finite chain for increasing n and several values of γ̃; if the difference fails to vanish as n grows—or if a direct calculation shows T*T≠π²I on L²[0,∞)—the nonlocal boundary terms would not emerge and the limit equation would differ. Alternatively, measure the steady-state temperature profile of a long open chain and compare it with the stationary solution of (2.7); a mismatch in the boundary layer would indicate the closure is wrong.","supporting_citations":[],"review_version":1}