{"id":"ecd5446b-4afd-4ac3-a1fa-407166b3f38f","arxiv_id":"2505.12034","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":8.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Exact large deviation functions and cumulants for the WASEP current and tracer statistics are obtained from the ASEP limit, with an explicit Lax pair mapping the MFT to an anisotropic Landau-Lifshitz spin chain.","lead":"This paper derives exact formulas for the rate functions and cumulants of particle current and tracer position in the weakly asymmetric exclusion process, connecting the symmetric and KPZ regimes. It also exhibits an integrable Lax-pair structure for the underlying macroscopic fluctuation theory, opening a new route to exact large deviation results.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Uncontrolled first-cumulant step: Eq. (SE.9) replaces a Fredholm determinant by exp(-Tr(ςA)) with no estimate of the same-order fluctuation terms, so Eq. (5) inherits an unproven self-averaging assumption.","rationale":"The paper's principal claims are that Eq. (5) gives the exact rate function and that the WASEP MFT is integrable. The Lax-pair part is concrete and verifiable: the zero-curvature statement is an algebraic consequence of (33), and it is not where the argument is fragile. The limit checks (SSEP in Section Q, KPZ in Section P) and the matching of κ1-κ3, with κ4 at ν=0, against independent perturbative MFT results are real evidence. The fragile point is exclusively the passage from the finite-q ASEP Fredholm determinant to the ε→0 rate function via (SE.9). The supplement calls this a first-cumulant approximation and asserts self-averaging, but no variance or rigidity estimate is supplied. If the omitted cumulants contribute at O(1/η), the central formula (5) would be off by an additive same-order term; if they contribute at O(1) or O(log(1/η)), the result stands. Because the SSEP and KPZ limits and the low-order cumulant matches could still be consistent with a correction that vanishes in those limits, they do not settle the issue. This is exactly the conditionality that the reader flagged, and it justifies keeping the CONDITIONAL verdict rather than promoting to ACCEPT; the requested second-cumulant computation should be a condition for full acceptance.","tokens_in":84029,"tokens_out":22489,"duration_ms":242084,"concrete_test":"Compute the next term in the cluster expansion of the Fredholm determinant, S_2 = (1/2) Tr((σ_u A)^2), where A = A_2 A_3 A_1 and σ_u(r) = 1/(1+q^{-r}/u), under the same rescaling q = e^{-η}, r = r̃/η, and steepest-descent analysis as Eqs. (SE.12)-(SE.21). If S_2 = O(1/η), then (SE.9) and hence Ψ(u) in Eq. (5) miss a same-order correction; if S_2 = o(1/η), e.g. O(1) or O(log(1/η)), the first-cumulant approximation is justified at leading order. This single analytic check settles whether the central formula stands or needs an additive correction.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central formula (5) is obtained from the \"first cumulant approximation\" of the ASEP Fredholm determinant, Eq. (SE.9) in the supplement. The exact object is log Det(I+K_u) = log E_DPP[∏_ℓ e^{-ς(a_ℓ)}] with ς = log(1+u q^r); replacing it by -Tr(ς A_2 A_3 A_1) discards, at the same formal order, the second and higher cumulants of the linear statistic, equivalently the terms -1/2 Tr((σ A)^2)+... in the cluster expansion, with σ = 1-e^{-ς}. Since η = 2νε is the only small parameter and the computed trace is O(1/η), any fluctuation contribution of order O(1/η) — not merely O(1) — changes Ψ(u) at leading order. The paper asserts that the linear statistics self-averages but gives no estimate of Var(S), no rigidity statement for the limiting point process, and no bound on the boundary terms in the integration by parts (SE.11). The independent checks (SSEP/KPZ limits, cumulants to n=3-4) are strong and make a total failure unlikely, but they do not control the full function Ψ(u). This is the load-bearing assumption of the main result.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript claims to derive exact large-deviation quantities for the weakly asymmetric exclusion process (WASEP) for a two-sided Bernoulli initial condition: a rate function Ψ(u) for the integrated current, the associated cumulants, the tracer-position statistics, and the SSEP/KPZ crossovers. The derivation starts from an ASEP Fredholm determinant formula of Aggarwal, converts it to a WASEP generating function in the limit ε→0, and then extracts cumulants and tails via Legendre-type parametric relations. In a second part, the paper transforms the macroscopic fluctuation theory (MFT) equations into a complex anisotropic Landau-Lifshitz spin system and exhibits an explicit Lax pair, thereby claiming classical integrability of the WASEP MFT and of related quadratic-mobility MFTs. The main quantitative results are benchmarked against independent perturbative and known limiting results. The central formula for Ψ(u) depends, however, on a first-cumulant approximation of a Fredholm determinant in the supplement, for which no error estimate is provided.","tokens_in":84462,"tokens_out":3616,"duration_ms":46150,"significance":"If the main formula (5) is valid, this is a significant advance: exact current and tracer large deviations for the WASEP across the full SSEP-to-KPZ crossover would be a natural benchmark for driven diffusive systems and macroscopic fluctuation theory. The paper ships a large amount of explicit, checkable material, including closed expressions for cumulants, parametric tail formulae, and a concrete Lax pair in Eqs. (36)-(37) whose zero-curvature condition is stated to be verified by direct computation. The consistency checks are strong and appropriately credited: cumulants match perturbative results up to n=3 (n=4 for ν=0), the SSEP limit matches [44,56], the weak-noise KPZ limit matches [50] up to fourth order, and the tracer variance matches [43, Eq. (S135)]. The value of the contribution is therefore high, but the central rate function rests on an unproven self-averaging step, so the paper is currently conditional rather than established.","major_comments":[{"comment":"The replacement Det(I+K_u) ≃ exp(-Tr(ς A_2 A_3 A_1)) is the load-bearing step for the main formula (5). The exact object is an expectation over a determinantal point process of exp(-Σ_ℓ ς(a_ℓ)), and the discarded terms are the second and higher cumulants of this linear statistic. Since η=2νε is the only small parameter and the retained trace is O(1/η), a fluctuation contribution of order O(1/η)—not merely O(1)—would change Ψ(u) at leading order. The text asserts self-averaging but provides no estimate of Var(Σ ς(a_ℓ)), no rigidity statement for the point process, and no bound on the neglected cluster terms. The strong low-order and limit checks do not control the full function Ψ(u). This should be fixed by an explicit error estimate, or by a proof of the self-averaging in the relevant η→0 limit, before the word 'exact' is attached to (5).","section":"Supplement E, Eq. (SE.9)"},{"comment":"The passage from the exact trace formula to the saddle-point integral (SE.21) needs justification of two points: the integration by parts in r states 'there is no boundary term', and the subsequent saddle-point evaluation assumes a unique relevant saddle with contour deformation between v and ζ and vanishing boundary contributions at infinity. The contour C0 is described as containing −1 as an asymptotic point, so boundary contributions are not obviously negligible in this oscillatory regime. A same-order boundary term would alter the O(1/η) result (SE.21) and hence the rate function. The authors should provide the missing estimates or state precisely which contour deformations are used.","section":"Supplement E.2, Eqs. (SE.11)-(SE.21)"},{"comment":"The general stationary cumulant formula in Appendix B is presented as an iterative arbitrary-order result, but its derivation in Appendix L relies on the conjecture (SL.47) and on the assumption (SL.49) that the jump at k=0 is exactly one-half of the jump at infinity 'for any q'. These are verified only for q=1,2 and tested in the KPZ limit. Since the claim 'arbitrary order' depends on these unproved identities, the paper should either prove them or explicitly label the general-order formula as conjectural with the stated support.","section":"Supplement L, Eqs. (SL.47), (SL.49), (SL.55)"},{"comment":"The paper claims the full large-deviation statistics including J<J_c, but the solitonic continuation is derived under a single-valued-branch assumption (SS.13), and the text states that the multi-branch phase-transition case is 'analyzed elsewhere [85]'. Also, the matching of Eq. (30) with the ASEP upper-tail result of [86] is presented as a demonstrated coincidence, while the supporting branch selection for the step initial condition is only sketched. These are limitation statements already contained in the manuscript; they should be elevated: either the missing cases are treated here, or the scope of the 'full statistics' claim should be narrowed accordingly.","section":"Main text, Tails paragraph; Supplement S.4-S.5, T.3"}],"minor_comments":[{"comment":"The variables z and J are used both for random variables and for their saddle-point values; footnote 74 acknowledges this, but the notation is still hard to follow in (10)-(12). A dedicated symbol for saddle-point values would improve readability.","section":"Eqs. (10)-(12) and Supplement H"},{"comment":"The main text says the stationary formulae 'allow iterative calculation of the cumulants to an arbitrary order', while the supplement labels the underlying identity as conjectural. The main text should carry the same caveat.","section":"Appendix B vs. Supplement L"},{"comment":"Panel (c) of Fig. SS.2 is printed in a dot-matrix style that makes the tri-valued branch structure nearly unreadable; a vector plot with labelled branches would be much clearer.","section":"Fig. SS.2"},{"comment":"The symbol ∼ is used for the large-deviation equivalence, but the left-hand side is an integral over ω of an expectation value; it would help to state explicitly which variable(s) are integrated out in the rate function and that the ω-integral is evaluated by saddle point.","section":"Eq. (4)"}],"recommendation":"major_revision","confidential_remarks":"The central result depends on a self-averaging approximation that is a technique developed by the authors in earlier papers, so independent validation is limited to low-order and limiting checks. The editor may wish to ask whether the missing error estimate for the first-cumulant step can be supplied within the Letter format, and whether the in-preparation reference [85] for the solitonic branches should be replaced by a self-contained treatment in the published version. The paper also contains an unusually high density of self-citations ([52-54,58,59,70,85,102-104]); this is not in itself a problem, but it amplifies the need for independent derivations of the load-bearing steps."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Right off: this is the missing WASEP piece. The paper produces exact large deviation functions for the integrated current and tracer in the weakly asymmetric exclusion process, and it gets the two limits right—SSEP and weak-noise KPZ. It also gives an explicit Lax pair mapping the MFT to the anisotropic Landau-Lifshitz chain, which is new and likely useful beyond this problem.\n\nThe strengths are real. The parametric formulas (5), (10)–(12) are concrete and the consistency checks are unusually thorough: cumulants match perturbative results to n=3 (n=4 at ν=0), the SSEP limit reproduces Derrida–Gerschenfeld and Mallick–Moriya–Sasamoto, the KPZ limit matches Le Doussal et al. up to fourth order, and the tracer variance matches Berlioz et al. The Lax pair (36)–(37) is explicit and the zero-curvature computation is given. That is a lot of solid work.\n\nThe soft spot is the step labeled 'first cumulant approximation' in the supplement (SE.9). The Fredholm determinant is replaced by exp(−Tr(ζA)), and the second and higher cumulants of the linear statistic are dropped. The trace is O(1/η), so if the discarded terms are also O(1/η), the central formula (5) changes at leading order. The paper asserts self-averaging but offers no estimate of the variance, no rigidity statement, and no bound on the boundary terms in the integration by parts. The successful limit checks are strong evidence that the final answer is correct, but they exercise only a few cumulants and the two endpoints; they do not control the full function Ψ(u). This is not a takedown—the result may well be exact, and the physics is clearly pointing the right way—but it is a genuine gap in the derivation as written. A referee should ask for an estimate, or for a clear statement that (5) is a conjecture backed by checks.\n\nSecondary, minor: the analytic continuation for J<Jc is deferred for general densities and only worked out for the step initial condition in a restricted parameter range. That is fine for a Letter, but it means the paper overclaims slightly when it says 'exact large deviations' without qualification.\n\nOverall: yes, send it to a serious referee. The paper is important and the checks make it very likely correct; the referee's job should be to close the first-cumulant gap or make its status explicit.","headline":"Exact WASEP large deviations with a clean integrability structure, but the central rate function rests on an unproven self-averaging step that a referee should pin down.","tokens_in":84911,"tokens_out":5866,"would_cite":true,"duration_ms":65319,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Exact large-deviation statistics are derived for the current and tracer in the weakly asymmetric exclusion process, interpolating between SSEP and weak-noise KPZ and revealing the integrability of the underlying macroscopic fluctuation…","keywords":["weakly asymmetric exclusion process","macroscopic fluctuation theory","large deviations","integrated current","tracer particle","Lax pair","Landau-Lifshitz model","Kardar-Parisi-Zhang crossover"],"falsifier":"Compute the next-order correction in $\\eta=2\\nu\\varepsilon$ to the log-determinant in Eq. (SE.9): if corrections appear at order $\\eta^0$ or $\\eta^1$ rather than vanishing, the saddle-point integral (SE.21) is incomplete. Alternatively, simulate the ASEP with rates $1\\pm\\nu\\varepsilon$ up to time $T/\\varepsilon^2$, measure the generating function (4) for several small $\\varepsilon$, and check whether $\\varepsilon\\log(\\cdot)$ converges to $-\\Psi(u)$; a non-vanishing slope in $\\varepsilon$ at fixed $u$ would falsify the claimed rate function.","tokens_in":83822,"feed_emoji":"📈","tokens_out":10698,"duration_ms":101985,"temperature":0.7,"pith_summary":"This paper establishes that the weakly asymmetric exclusion process (WASEP), the driven one-dimensional lattice gas with hopping asymmetry scaled as $1/\\sqrt{T}$, has exact large-deviation statistics rather than only perturbative cumulants. The core object is a rate function $\\Psi(u)$ written as a dilogarithm contour integral over the initial densities; from it, by Legendre-type transforms, the paper obtains the rate functions and all cumulants of the time-integrated current and of a tracer's position. These formulas interpolate between two known worlds: the symmetric exclusion process when the driving $\\nu\\to0$ and the short-time weak-noise Kardar-Parisi-Zhang equation when $\\nu\\to\\infty$, including the crossover of the cubic tail to the $5/2$ and $3/2$ tail exponents. The paper also proves that the macroscopic fluctuation theory of the WASEP is classically integrable by exhibiting an explicit Lax pair through a mapping to the anisotropic Landau-Lifshitz spin chain. If correct, this turns the WASEP into a reference model where full current statistics are available at any driving strength.","feed_headline":"WASEP current and tracer statistics solved exactly from SSEP to KPZ","feed_subtitle":"Current and tracer statistics are now exact for WASEP, with the full SSEP-to-KPZ tail crossover.","key_machinery":"The load-bearing machinery is the first-cumulant reduction of the ASEP Fredholm determinant (SA.1): in the WASEP limit $\\eta=2\\nu\\varepsilon\\to0$, the log-determinant is replaced by the trace of the kernel (SE.9), whose saddle-point evaluation produces the dilogarithm integral (5). Around this sits a calculus that converts the unusual generating function (4) into ordinary cumulants: the auxiliary variable $\\omega$ with weight $F(\\omega)$ is eliminated by a saddle point, and the Legendre system (10)-(12) turns $\\Psi(u)$ into $\\Phi(J)$ and $\\phi(P)$. On the integrability side, the key object is the change of variables (32) taking the MFT fields $(q,p)$ to $(Z,R)$; the matrix (35) realises these fields as a stereographic parametrisation of a complex anisotropic Landau-Lifshitz spin, and the zero-curvature pair (36)-(37) is the Lax representation whose compatibility reproduces the MFT equations.","core_discovery":"On the paper's own terms, the central claim is that the generating function (4), built from the exact ASEP Fredholm determinant, collapses in the WASEP scaling to the simple integral $\\Psi(u)$ in Eq. (5), and that the parametric system (10)-(12) then yields the complete large-deviation functions $\\Phi(J)$ and $\\phi(P)$ for the integrated current, from which the tracer position statistics follow by setting $J=0$. The same formula is shown to carry known results: the $\\nu\\to0$ limit reproduces the SSEP formulas of [44,56] and the $\\nu\\to\\infty$ limit reproduces the weak-noise KPZ rate function for Brownian initial data [50], with an explicit crossover tail function between the $J^3$ and $J^{5/2}$ tails and a $3/2$ tail of the upper current distribution. In addition, the MFT saddle-point equations (31) are shown to be integrable: the change of variables (32) maps them to a complex extension of the anisotropic Landau-Lifshitz model, with Lax pair (36)-(37), so the authors claim integrability for all MFTs with quadratic mobility and their duals.","pith_inferences":["If the first-cumulant self-averaging is exact, the same trace-only reduction should apply to other $q$-deformed Fredholm determinants in the same scaling, giving analogous rate functions for $q$-TASEP-like models.","The explicit crossover tail is a quantitative experimental prediction for driven single-file systems: the rescaled current PDF should collapse onto the universal function $\\Phi_+$ with no intermediate exponent.","The soliton-branch continuation for $J<J_c$ suggests a dynamical phase transition in the WASEP current statistics at large $\\nu^2T$; locating the non-analyticity of $\\Phi(J)$ would test whether the integrable structure, not a saddle-point artefact, controls the transition.","Because the Lax equation is gauge-equivalent to standard integrable nonlinear Schrödinger-type flows, spectral methods may yield Riemann-Hilbert or Fredholm representations of the rate function and connect these MFT large deviations to random-matrix edge statistics."],"forward_implications":["Every cumulant of the integrated current and of the tracer can be written in closed form; the coefficients $\\Psi_n(X,T)$ obey the order-$n$ heat equation (6), so the formulas propagate from $T=0$ data.","The $\\nu\\to0$ limit reproduces the SSEP current and tracer CGFs, including formulas at arbitrary observation point $X$, and the $\\nu\\to\\infty$ limit reproduces the weak-noise KPZ rate function, so the WASEP is a single model connecting both universality classes.","The right tail of the current remains cubic for fixed driving, and as the driving grows the tail crosses over, through the explicit function (28)-(29), to the KPZ $5/2$ tail; the upper tail shows a $3/2$ wall solution.","Integrability extends to the WASIP, weakly asymmetric KMP, and their duals, giving exact current and tracer statistics for all these models at once.","The explicit Lax pair opens the way to inverse-scattering solutions of the WASEP MFT, which would cover more general initial conditions and observables than the two-sided Bernoulli case solved here."],"supporting_citations":[{"why":"Supplies the exact ASEP Fredholm determinant identity for the current, the microscopic starting point of the whole derivation.","marker":"[99]"},{"why":"Gives the alternative kernel representation used in the contour manipulations and in taking the WASEP limit.","marker":"[101]"},{"why":"Develops the first-cumulant method that converts the Fredholm determinant asymptotics into the trace of the kernel, the step that produces Eq. (5).","marker":"[102–104]"},{"why":"Provides the SSEP step-initial-condition current rate function that the $\\nu\\to0$ limit must match.","marker":"[44]"},{"why":"Provides exact SSEP MFT rate functions for arbitrary $X$; the paper checks its $\\Psi_0(u)$ against the result there.","marker":"[56]"},{"why":"Supplies the weak-noise KPZ rate function for Brownian initial conditions used as the $\\nu\\to\\infty$ benchmark.","marker":"[50]"},{"why":"Gives perturbative low-order WASEP cumulants used to check the new cumulant formulas.","marker":"[43]"},{"why":"Proves the integrability of the anisotropic Landau-Lifshitz model from which the explicit Lax pair is imported.","marker":"[92]"},{"why":"Provides the ASEP upper-tail large deviation result that the large-$\\nu$ step-initial-condition tail must coincide with, supporting the $3/2$ exponent claim.","marker":"[86]"}],"fun_headline_variants":["Exact WASEP current/tracer large deviations via integrable MFT","WASEP MFT integrable: Lax pairs from spin chain mapping","One formula for WASEP: full SSEP to KPZ tail crossover","WASEP rates exact: current and tracer from SSEP to KPZ","Integrable WASEP: exact large deviations from MFT"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole result rests on replacing the log of the ASEP Fredholm determinant by the trace of its kernel in the limit $\\varepsilon\\to0$; if that self-averaging step has non-vanishing subleading corrections, the central rate function $\\Psi(u)$ is not exact.","fun_headline_variants_meta":{"raw":{"variants":["Exact WASEP current/tracer large deviations via integrable MFT","WASEP MFT integrable: Lax pairs from spin chain mapping","One formula for WASEP: full SSEP to KPZ tail crossover","WASEP rates exact: current and tracer from SSEP to KPZ","Integrable WASEP: exact large deviations from MFT"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000325,"raw_usage":{"total_tokens":1866,"prompt_tokens":1037,"completion_tokens":829,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":653,"completion_tokens_details":{"reasoning_tokens":734}},"tokens_in":653,"tokens_out":829,"duration_ms":7923,"temperature":1.0,"reasoning_tokens":734,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T20:42:59.156165+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the next-order correction in $\\eta=2\\nu\\varepsilon$ to the log-determinant in Eq. (SE.9): if corrections appear at order $\\eta^0$ or $\\eta^1$ rather than vanishing, the saddle-point integral (SE.21) is incomplete. Alternatively, simulate the ASEP with rates $1\\pm\\nu\\varepsilon$ up to time $T/\\varepsilon^2$, measure the generating function (4) for several small $\\varepsilon$, and check whether $\\varepsilon\\log(\\cdot)$ converges to $-\\Psi(u)$; a non-vanishing slope in $\\varepsilon$ at fixed $u$ would falsify the claimed rate function.","supporting_citations":[],"review_version":1}