{"id":"8d2346c0-4107-4a64-818b-ebb855ab070b","arxiv_id":"2607.28807","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":8.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Exact finite-time density formulas, a finite-dimensional moment algebra, and a rigorous Euler-scale limit with Catalan statistics are derived for the symmetric Dyson exclusion process.","lead":"This paper derives exact formulas for the density evolution of a long-range interacting exclusion process, the symmetric Dyson exclusion process, from deterministic initial states. It proves the large-time Euler-scale limit and the arctic curve, confirming a previous conjecture from the microscopic dynamics.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Euler-scale theorem rests on a sketched two-saddle steepest-descent estimate (Lemma 6.4); the O(N^{-1/2}) remainder is not established.","rationale":"The reader identified two gaps: the Fourier-support claim in Prop. 3.1 and the steepest-descent Lemma 6.4. I agree that Lemma 6.4 is the more load-bearing issue for the central claim, because Theorem 6.2 relies on its uniform error bound. The Fourier-support claim, while sketched, is directly checkable from the trigonometric product (3.1) and appears correct; it is not the principal risk. The steepest-descent lemma, however, is not a routine one-saddle estimate: the phase has two conjugate saddles with equal real parts, and the pole at z=ω coincides with the saddles. The proof's one-sentence invocation of the 'usual construction' does not supply the necessary uniform bounds, particularly for the O(N^{-1/2}) remainder. This does not mean the result is false: the exact formulas (5.20) and (6.15) are credible, and the agreement with the conjectured hydrodynamic arctic curve [23] and with finite-N numerics (Figure 1) provides independent support. But the advertised rigor is incomplete, so the paper should be conditional on a full proof or a more detailed numerical/analytic verification of Lemma 6.4. Hence the reader's CONDITIONAL verdict is appropriate, and I do not recommend a change.","tokens_in":19889,"tokens_out":15703,"duration_ms":161096,"concrete_test":"Evaluate the exact double-contour integral (6.15) numerically for a fixed interior point (ξ,τ) (e.g., ξ=0.5, τ=0.1) for N=20,40,80,160, using high-precision quadrature on the original contours Γω and Γz. Compare K_N with -arg u/π. If |K_N + arg u/π| does not decay like O(N^{-1/2}) or faster, the Euler-scale claim is in question. To isolate Lemma 6.4, independently construct the deformed steepest-descent contours, subtract the crossing contribution ∫Σ dω/(2πiω), and measure the remainder; check that it is uniformly O(N^{-1/2}) over a range of (ξ,τ) in L. An analytic local expansion near (u,u), (u,\\bar u), (\\bar u,u), (\\bar u,\\bar u) after crossing removal would also settle the uniformity claim.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Theorem 6.2, the paper's central Euler-scale result, depends on Lemma 6.4, whose proof is a one-paragraph reference to 'the usual two-contour steepest-descent construction.' The lemma asserts that after deforming the contours in (6.15) and removing the crossing along Σ, the remaining double integral eK_N is uniformly O(N^{-1/2}) on compact subsets of L. This is not a routine single-saddle estimate. The phase Ψ_{ξ,τ}(a)=τ/2(a+a^{-1})-ξ log a+log(a-1) has two conjugate non-degenerate saddles u and \\bar u. Since the coefficients are real, Re Ψ(u)=Re Ψ(\\bar u), so the double integral has four stationary-phase combinations (ω,z)∈{u,\\bar u}^2, not only the diagonal crossing. Near (u,u), the pole at z=ω coincides with the critical point; the local integral with the singular factor 1/(z-ω) contributes at order N^{-1/2}, and it is not shown that this is removed by the crossing subtraction or bounded uniformly on L. The off-diagonal saddle combinations are likewise not analyzed. Thus the asserted uniform error bound, and hence the claimed limit (6.3), is not fully demonstrated. The abstract's promise of a 'rigorous steepest-descent analysis' makes this gap load-bearing. By contrast, Proposition 3.1's Fourier-support claim is checkable from the explicit Laurent form of (3.1) and appears sound; it is not the critical risk.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the symmetric Dyson exclusion process (SDEP), a lattice exclusion process with long-range logarithmic jump rates, using the ground-state transform that maps it to the free-fermion XX chain. The authors derive exact equal-time correlation formulas for deterministic initial configurations, expressing the mixed determinantal kernel through trigonometric Lagrange interpolation. From this they obtain finite- and infinite-volume density evolution equations, a finite-dimensional polynomial algebra for all density moments, Catalan asymptotics for the highest-time coefficients, and an Euler-scale limit shape with an arctic curve for the melting of a compact block. The central claim is that the Euler-scale density profile is ρ = arg u(ξ,τ)/π, where u is the upper-half-plane root of the cubic (6.2), and that the arctic boundary is the discriminant locus (6.4), coinciding with the earlier hydrodynamic conjecture in [23].","tokens_in":20269,"tokens_out":11122,"duration_ms":117148,"significance":"If the claims are fully justified, this is a substantial contribution to the exact theory of long-range interacting exclusion processes. The finite-time formulas (Theorems 4.1 and 5.3) are explicit and parameter-free, and the polynomial moment algebra is a genuinely new structural result. The Euler-scale theorem provides the first microscopic derivation of the conjectured arctic curve of [23] and places the SDEP in the broader context of free-fermion limit-shape problems. The paper is careful in many places, especially the derivation of the moment algebra and the exact double-contour representation. However, the rigorous status of the central Euler-scale theorem currently depends on a sketch rather than a proof, so the work needs revision before it can be regarded as closing the advertised 'rigorous steepest-descent analysis.'","major_comments":[{"comment":"Theorem 6.2, the paper's headline Euler-scale result, rests on Lemma 6.4, whose proof is a one-sentence reference to 'the usual two-contour steepest-descent construction.' This is not a routine estimate. In (6.15) the integrand has a singular factor 1/(z−ω) and the phase Ψ_{ξ,τ} has two conjugate saddles u and ū; because Re Ψ(u)=Re Ψ(ū), the double integral has four stationary combinations (u,u), (u,ū), (ū,u), (ū,ū), not only the diagonal crossing. At (u,u) the pole coincides with the saddle, so the local integral with 1/(z−ω) is a singular saddle-point problem; the text does not show that the crossing subtraction removes the O(N^{−1/2}) contribution, nor does it analyze the off-diagonal combinations. The uniform bound (6.16) is therefore not established, and the limit (6.3) remains conditional on a missing calculation. The authors should give a complete proof of Lemma 6.4.","section":"§6.1.2, Lemma 6.4"},{"comment":"The proof of the finite-volume kernel (3.1) asserts without demonstration that the trigonometric Lagrange function L_y(x) is a Laurent polynomial 'with Fourier modes contained in the Fermi sea.' This is not immediate: writing L_y in terms of q=e^{2π i x/L} produces a prefactor q^{−(N−1)/2} times a polynomial, and the cancellation must be shown explicitly. Since every subsequent formula (Theorems 4.1, 5.3, and 6.2) inherits Prop. 3.1, this step should be proved in the text rather than asserted.","section":"§3.1, Prop. 3.1"},{"comment":"In Prop. 5.8, the combinatorial identity (5.30), which is the only input that makes the average over balanced words vanish, is stated without proof. The sentence on reversal is not enough: one must define h_q for the reversed word and verify the height transform h ↦ 1−h, including words that go below 0. Please supply the calculation; as written, the Catalan claim κ_{r,N}=C_r N^{r+1}+O_r(N^{r−1}) is not fully justified.","section":"§5.5, Prop. 5.8"}],"minor_comments":[{"comment":"The sentence 'Start from the exact representation (6.15).' is duplicated verbatim at the beginning of the subsection.","section":"§6.1.3"},{"comment":"The notation I_n(t) is introduced for the modified Bessel function, whereas Section 4.2 uses I_x(−2wt). After setting w=1/2 the conventions are consistent, but the switch is abrupt; please state explicitly that t in Section 6 is the physical time with w=1/2.","section":"Eq. (6.6)"},{"comment":"The right panel reports 'exact finite-lattice result (5.20) at N=20.' If these are direct numerical evaluations of the exact formula, state so; if any quadrature or truncation is used, it should be described.","section":"Figure 1"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth engaging with. The exact finite-time formulas are genuinely new, and the moment algebra is the cleanest part of the paper. The Euler-scale theorem is probably correct, but the steepest-descent proof is a sketch, and the abstract promises rigor it doesn't yet deliver.\n\nThe core input is Proposition 3.1, the trigonometric Lagrange-interpolation formula for the mixed kernel. That is a nice closed form, and despite a proof that hand-waves the Fourier-support claim, the claim itself is checkable from the explicit Laurent polynomial and appears true. The moment algebra in Section 5 is the strongest section: Theorem 5.3 (moments as traces of (X+sD)^n) is elegant, the degree bound is proved, and the first two universal moments plus the Catalan asymptotics in Proposition 5.8 are clearly derived. I checked the block formulas (5.23)-(5.25); they are consistent with the trace algebra, and the finite representation in Remark 5.6 is a real improvement over an infinite Bessel sum.\n\nThe soft spot is Section 6. Lemma 6.4, which supplies the uniform O(N^{-1/2}) bound for the double contour integral after the crossing contribution is removed, is not proved — the text says it follows from 'the usual two-contour steepest-descent construction' and a short local Gaussian estimate. That is not enough. The phase has two conjugate non-degenerate saddles; the pole at z=ω sits exactly at the (u,u) critical point, and the four saddle combinations (u,u), (u,\\bar u), (\\bar u,u), (\\bar u,\\bar u) all need uniform control. The crossing subtraction along Σ may well take care of the (u,u) contribution, but the paper does not show it, and the off-diagonal combinations are not analyzed at all. Since Theorem 6.2 depends on this lemma, the central Euler-scale result is not fully rigorous as written. I want to emphasize I don't think the theorem is wrong: it matches the companion-paper conjecture and the numerics, and the reduction to crossing gives the right arg/π. But the abstract's claim of a rigorous steepest-descent analysis is premature.\n\nThe citation pattern is fine — the companion paper is used for comparison and context, not fitted. I see no circularity.\n\nWho should read this: anyone working on long-range exclusion processes, determinantal processes, or free-fermion hydrodynamics. It deserves a serious referee. Recommend accept for peer review with a request for a detailed proof of Lemma 6.4 (and a one-line proof or reference for the Fourier-support claim in Prop 3.1).","headline":"Solid exact results for a long-range exclusion process, with the Euler-scale theorem hanging on a steepest-descent lemma that is sketched rather than proved.","tokens_in":20684,"tokens_out":5614,"would_cite":true,"duration_ms":60091,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["60K35","82C22","82B23"],"pacs":["05.40.-a"],"model":"deepseek-v4-flash","headline":"The symmetric Dyson exclusion process is solved exactly from deterministic initial states, yielding the Euler-scale block-melting profile and its arctic curve.","keywords":["symmetric Dyson exclusion process","long-range interacting particle system","free-fermion representation","determinantal process","Lagrange interpolation","arctic curve","Euler scaling","Catalan numbers"],"falsifier":"Evaluate the exact double-contour representation (6.15) numerically at finite N for points (xi,tau) inside the liquid region and compare the result with arg u/pi; if the difference does not decay as N^{-1/2}, the steepest-descent lemma fails. Alternatively, simulate the SDEP block melting at Euler scale and measure the density profile: any consistent deviation from the discriminant curve (6.4) would falsify the derived arctic boundary.","tokens_in":19817,"feed_emoji":"🧊","tokens_out":6109,"duration_ms":64993,"temperature":0.7,"pith_summary":"The paper takes on the symmetric Dyson exclusion process (SDEP), a lattice gas in which particles exclude each other and also repel through a long-range logarithmic Coulomb interaction—equivalently, symmetric random walkers conditioned never to collide. It shows that because the process is the ground-state Doob transform of the free-fermion XX chain, every equal-time observable from a deterministic initial configuration is a determinant built from a mixed kernel. That kernel is evaluated in closed form by trigonometric Lagrange interpolation, which turns the melting of a fully packed block into a finite computation. The paper then derives the exact density profile in the Euler scaling limit as arg u/pi, where u is the upper-half-plane root of an explicit cubic, and identifies the arctic curve—the boundary between frozen and liquid regions—as the locus where two roots coalesce. These results confirm a previously conjectured hydrodynamic description and bring out Catalan numbers as the leading long-time coefficients of the density moments.","feed_headline":"Exact solution fixes the long-range gas's arctic curve","feed_subtitle":"Euler-scale density of a melting block now follows from the microscopic formula, confirming the conjectured hydrodynamic boundary.","key_machinery":"The load-bearing objects are: (i) the ground-state (Doob) transform conjugating the SDEP generator to the XX free-fermion Hamiltonian, which converts SDEP expectations into mixed free-fermion matrix elements; (ii) the trigonometric Lagrange interpolation formula for the initial mixed kernel, which turns deterministic initial data into a closed interpolation problem; (iii) the operator identity e^{sB} X e^{-sB} = X + sD acting on the space of polynomials of degree < N, which makes the whole moment hierarchy finite-dimensional; and (iv) an exact double-contour representation whose phase has critical points exactly the roots of the cubic (6.2). The Euler-scale limit is obtained by a steepest-de","core_discovery":"The central claim is that the block-melting problem for the SDEP admits an exact finite-time formula and a computable macroscopic limit. Starting from the block {1,...,N}, the density at time t is the finite sum (5.20) of products of modified Bessel functions times block Lagrange polynomials. Under Euler scaling x_N = floor(xi N) and t_N = tau N, this density converges uniformly on compact subsets of the liquid region to (1/pi) arg u(xi,tau), where u is the upper-half-plane root of the cubic (6.2). The boundary of the liquid region is the discriminant locus of that cubic, equation (6.4), which separates the frozen phases rho=0 and rho=1 from the liquid region 0