{"id":"7fff257d-ac2b-4afd-b06f-a30d4d69e9fd","arxiv_id":"2505.06353","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A duality between a coalescing population model and a noisy long-range FKPP equation, supported by large-scale simulation, shows identical butterfly light cone scaling for alpha>0.5 in 1D.","lead":"This paper connects two mathematical descriptions of how quantum operators spread in long-range interacting systems, showing via a duality that the discrete population picture and a noisy continuous reaction-diffusion equation produce the same butterfly light cone. The result gives a unified tool for predicting chaos propagation in quantum simulators where interactions decay as a power law.","discovery_kind":"unification","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Conclusion hinges on numerically demonstrated but unproven soft/hard universality, with the claimed α>0.5 range untested below α=0.85; the Conclusion's Lévy index (µ=2α−2) contradicts Eq. (4) and Appendix D (µ=2α−1).","rationale":"The reader's verdict correctly identifies the weakest assumption: the equivalence between the hard-constraint population model and the soft-coalescence model is numerically demonstrated but not proven. My reading of the manuscript confirms this is the place where the central claim is least secure. The exact duality between the soft model and the FKPP equation is a rigorous construction for a given Markov process, but the hard model is the physically motivated operator-growth dynamics, and the leap from hard to soft is the unproven step. The numerical evidence is substantial in the tested regimes, but the paper claims α>0.5 while showing no data below α=0.85; that is a concrete, falsifiable gap rather than a stylistic concern. The Lévy index inconsistency is a separate, real defect: Eq. (4) and Appendix D.5.b say µ=2α−1, the Conclusion says µ=2α−2, and the paper even claims the latter contrasts with previous work giving µ=2α−1, which would contradict its own derivation. This needs correction, but it does not change the verdict because the scaling comparison is performed on the population models, not on the continuum PDE. A CONDITIONAL verdict with high confidence remains appropriate: the paper should either prove or sharply extend the numerical evidence for soft/hard universality, and fix the index typo, before the unification claim is stated as a firm conclusion.","tokens_in":26175,"tokens_out":8438,"duration_ms":84474,"concrete_test":"Run the paper's sparse-window Gillespie algorithm for both the hard (N=2) and soft (equilibrium height≈2) constraint models at α = 0.60, 0.70, and 0.80, extending t as far as the no-finite-size data structure allows (target t ~ 10^3–10^4 in the stretched-exponential regime). If the log ℓ(t) curves for the two models fail to collapse under the predicted exp(const × (ln t)^2) scaling, the soft/hard universality assumption is violated and the central claim must be narrowed. In parallel, re-derive the fractional Laplacian order from the dual SDE in Appendix D.5.b to determine whether Eq. (4)'s µ=2α−1 or the Conclusion's µ=2α−2 is the correct continuum statement.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that the two perspectives share identical butterfly light cone scaling for α>0.5 in 1D rests on a chain: the hard-constraint model (the actual operator-growth picture) is equivalent to the soft-coalescence model, which is exactly dual to the stochastic long-range FKPP equation. The first link is the load-bearing one. The paper explicitly states it cannot implement the hard constraint in the dual model ('We are not able to implement this constraint exactly in the dual model') and that no dual process is known for hard-constraint arrivals (Appendix D, Section f). The soft-hard equivalence is therefore an unproven universality assumption, supported only by numerical evidence. That evidence in Fig. 4 covers α = 0.85–1.05 and α = 1.3–2.0, while the conclusion announces α>0.5; the interval 0.5<α<0.85, where the predicted stretched-exponential exponent diverges most rapidly, is not tested. Fig. 5 further shows that other front statistics (Np, Ns, M) do not exactly overlap between the two models, so the ℓ(t) overlap in Fig. 4 is the only direct evidence for the claimed identity. Separately, the paper is internally inconsistent about the Lévy index: Eq. (4) and Appendix D.5.b give µ = 2α−1, but the Conclusion states µ = 2α−2. This inconsistency does not by itself invalidate the numerical scaling comparison, but it must be resolved before the FKPP equation is unambiguously identified.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes to unify two complementary descriptions of operator spreading in chaotic long-range interacting systems: the discrete stochastic population-dynamics picture with a hard local height cap, and the continuum noisy long-range FKPP equation. The authors introduce an intermediate soft-constraint model in which the hard cap is replaced by a coalescence process that fixes the equilibrium population. They prove, in Appendix D, an exact duality between this soft-constraint model and a coupled system of diffusions that, in the continuum limit, becomes the stochastic long-range FKPP equation. The light-cone equivalence follows from the duality. Since the hard-constraint model is the one directly tied to the operator-growth picture, the authors then compare hard- and soft-constraint models numerically in one dimension, finding overlapping butterfly light-cone curves for the tested parameters, and they derive a third-order perturbative scaling around α=1. The central claim is that, for α>0.5 in 1D, the discrete operator-growth picture and the noisy long-range FKPP description share identical butterfly light-cone scaling.","tokens_in":26516,"tokens_out":3329,"duration_ms":34106,"significance":"The exact duality in Appendix D is a genuine contribution: it is mathematically clean, holds for arbitrary population size N, and gives a direct bridge between a discrete Markov process and a stochastic PDE with multiplicative noise. The numerical algorithm is also a strength: it simulates up to 10^9 sites with event selection in constant time and no finite-size space cutoff, which is necessary for the stretched-exponential regime near α=1. The third-order perturbative solution of the iterative light-cone equation around α=1 is a useful benchmark for a regime that is extremely hard to access asymptotically. If the soft-hard equivalence were proven or convincingly established over the full claimed range, the unification would be an important conceptual step for the operator-spreading literature. However, the paper currently overreaches in its headline range α>0.5 and contains an internal inconsistency in the Lévy index, so the significance is conditional on those points being fixed.","major_comments":[{"comment":"There is an internal inconsistency in the Lévy index of the claimed FKPP equation. Equation (4) states μ=2α−1, and Appendix D.5.b derives the same index by writing ∂t f = −(1−f) D (−Δ)^{α−1/2} f + λ f(1−f) + sqrt(β f(1−f)) η. The Conclusion, however, states that the resulting long-range FKPP equation has Lévy index μ=2α−2 and contrasts this with Ref. [34], which supposedly gives μ=2α−1. The last sentence would mean the paper's own derivation gives the same index as the work it claims to contradict. This must be resolved before the FKPP equation is unambiguously identified, because the title and abstract present that equation as one of the two unified perspectives.","section":"Conclusion, Eq. (4), and Appendix D.5.b"},{"comment":"The load-bearing link in the unification is the equivalence between the hard-constraint model (the operator-growth picture) and the soft-constraint coalescence model. The authors state in the Introduction that they cannot implement the hard constraint exactly in the dual model, and Appendix D.4.f explicitly concludes that no dual process exists for arrival-dependent hard-constraint births. The soft-hard equivalence is therefore a numerical conjecture, not a theorem. The numerical evidence covers α=0.85–1.05 and α=1.3–2.0, while the Conclusion announces the result for all α>0.5 in 1D. The interval 0.5<α<0.85, where the predicted stretched-exponential exponent diverges most quickly, is not tested. Moreover, Fig. 5 shows that other front statistics (Np, Ns, M) do not exactly overlap between the two models, so the ℓ(t) overlap in Fig. 4 is the only direct scaling-level evidence. I recommend either restricting the claim to the tested α range or adding simulations and, if possible, an analytic argument for the missing interval.","section":"Introduction and Appendix D.4.f; Figs. 4 and 5"},{"comment":"The validation of the α=1 prediction is partly circular. The constant c1 in Eq. (7) is computed from K, and K is read from the plateau in Fig. 4(c), which is obtained from the same simulation data used to evaluate the prediction in Fig. 4(b); in addition, the constant C is adjusted by a vertical shift. The authors disclose this, and the test is still meaningful for the functional form, but the text's phrase 'no fitting of unknown constants' (near Fig. 4) is too strong. The agreement should be described as a consistency check of the predicted scaling form with a constant determined from the same data, not as an independent parameter-free confirmation.","section":"Fig. 4(b), Eq. (7), and Fig. 4(c)"}],"minor_comments":[{"comment":"The word 'varitable' appears in the sentence introducing the change of variable τ=t/2+ut; it should be 'variable'.","section":"Appendix E.1"},{"comment":"The caption says the dashed line is Eq. (7) with an adjusted value for C, while the text says the only fitting is a global vertical shift; please reconcile these descriptions so the reader knows exactly which constants are free.","section":"Fig. 4 caption and text near Eq. (7)"},{"comment":"The observation that Np, Ns, and M do not exactly overlap between soft and hard constraints is important and should be moved or at least highlighted in the main text, because it delimits the claimed equivalence to ℓ(t) scaling rather than to full statistical equivalence of the two stochastic processes.","section":"Appendix A and Fig. 5"},{"comment":"The passage from the lattice SDE (D38)–(D39) to the continuum fractional Laplacian equation (D40) is heuristic; a sentence clarifying the sense in which this limit is taken (e.g., Fourier/normalization of the kernel near the origin) would help the reader assess the status of the resulting FKPP equation.","section":"Appendix D.5.b"},{"comment":"The scaling functions in Fig. 2 for the hard constraint model are stated without derivation; since the paper relies on them for the phase boundaries, a brief derivation or a precise citation to the exact results of Refs. [8,31,33] for each regime would improve self-containedness.","section":"Fig. 2"}],"recommendation":"major_revision","confidential_remarks":"The exact duality for the soft-constraint model is the strongest part of the paper and is likely correct. The main risk is overclaiming: the α>0.5 statement goes beyond the numerically tested range, and the Lévy-index inconsistency in the Conclusion must be fixed. I believe these are fixable within the manuscript's scope, hence major revision rather than reject."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Punchline: the exact duality is real, and the numerical work is serious, but the conclusion overreaches. The paper proves an exact duality between the soft-constraint population model and the stochastic long-range FKPP equation, and that is the part worth building on. The numerical comparison between soft and hard constraint models is extensive and well done—the 10^9-site sparse algorithm is a genuine technical achievement.\n\nWhat is new: extending the Doering-Mueller-Smereka duality to birth-at-distance interactions, the third-order perturbative expansion near alpha=1, and the large-scale evidence that the soft and hard models have the same light-cone scaling where tested. The derivation in Appendix D is clean, and the honest admission that the hard constraint has no dual counterpart is to the authors' credit.\n\nSoft spots, in proportion: the hard-soft equivalence is numerically demonstrated, not proven, and the paper says so plainly. That makes the unification claim conditional, not established. The conclusion's \"for alpha>0.5 in 1D\" goes beyond the data—the simulations cover roughly alpha=0.85 to 2.0, and the 0.5 to 0.85 window is where the stretched exponential exponent changes fastest. That range simply has not been tested. There is also a genuine internal inconsistency: Eq. (4) and Appendix D give the Levy index as 2 alpha - 1, while the Conclusion states 2 alpha - 2. That is not a trivial typo—it changes the identification of the PDE—and it must be resolved before the result is cited as a unification. The fitting of K from the same data used to validate the prediction is a milder concern; it is disclosed and the plateau in Fig. 4(c) makes it reasonably convincing.\n\nBottom line: the duality result deserves to be published and cited. The soft-hard universality is a strong numerical claim that is very plausible but not proven; the paper should say \"for alpha in the tested range\" and fix the Levy index. This would benefit from a serious referee—send it out.","headline":"The exact soft-constraint duality is real and the numerics are serious, but the hard-soft universality is numerically supported only, and the Levy-index inconsistency in the Conclusion needs fixing before the unification claim is taken as established.","tokens_in":27038,"tokens_out":2428,"would_cite":true,"duration_ms":22396,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper establishes that discrete operator-growth population dynamics and the continuum noisy long-range FKPP equation describe the same butterfly light cone in power-law quantum systems, by proving an exact duality through a…","keywords":["operator spreading","butterfly light cone","long-range FKPP equation","population dynamics","duality","coalescence","quantum chaos","power-law interactions"],"falsifier":"Simulate hard- and soft-constraint models in one dimension at several $\\alpha$ values in $(0.5,1.5)$ with the reported no-finite-size algorithm, run until the iterative integral equation plateaus, and compare the inferred exponents $d\\log\\ell/d\\log t$; any difference outside numerical error would falsify the claimed equivalence.","tokens_in":1802,"feed_emoji":"🦋","tokens_out":5244,"duration_ms":97012,"temperature":0.7,"pith_summary":"This paper tries to close the gap between two standard ways of describing operator spreading in chaotic quantum systems: a discrete population process that counts how many nontrivial operators live on each site, and a continuum stochastic reaction-diffusion equation of FKPP type. The authors construct an exact duality between a population process and the stochastic FKPP equation by replacing the hard cap on local population with a soft constraint, namely pairwise coalescence that keeps the equilibrium population at $N$. Under that soft model the duality is a mathematical identity, a time-independent combined correlator that forces the butterfly light cones of the two descriptions to coincide. The hard-cap model is not exactly dualizable, so the final step is numerical: in one dimension the soft and hard models show matching light-cone scalings for $\\alpha>0.5$. A sympathetic reader would care because if the chain holds, the discrete growth picture and the noisy PDE are equivalent descriptions of the same physics, and one can choose whichever is easier to compute.","feed_headline":"Two views of quantum chaos share the same light cone","feed_subtitle":"A soft-constraint population model is exactly dual to the noisy long-range FKPP equation, unifying the two views.","key_machinery":"The load-bearing object is the duality identity between the soft-constraint population model and the stochastic long-range FKPP equation. On one side is an integer height $h(x,t)$ with birth-at-distance rate $G(\\|y-x\\|)h(x)$ and on-site coalescence rate $\\tfrac{\\beta}{2}h(h-1)$; on the other side is a continuum field $f(x,t)\\in[0,1]$ satisfying $\\partial_t f = -(1-f)\\tilde D(-\\Delta)^{\\alpha-1/2}f + \\tilde\\lambda f(1-f) + \\sqrt{\\tilde\\beta f(1-f)}\\,\\eta$. The identity $\\langle \\prod_x (1-f(x,\\tau))^{h(x,t-\\tau)}\\rangle = \\text{const}$ transfers the light cone of one process to the other. It is the exact, parameter-free link between the discrete and continuum perspectives; the hard constraint has no such dual, so the soft model is the hinge of the whole argument.","core_discovery":"The central claim is that the butterfly light cone of operator growth in a system with power-law interactions is the same object whether computed from a discrete birth-at-distance population process or from a stochastic long-range FKPP equation. The bridge is an intermediate population model in which the hard local maximum $N$ is replaced by a coalescence process $A+A\\to A$ with equilibrium population $N$. For this soft model the paper proves a duality: a combined correlator $\\langle \\prod_x (1-f(x,\\tau))^{h(x,t-\\tau)}\\rangle$ is independent of $\\tau$, which equates the event that spreading has reached distance $L$ with the event that the dual field $f$ has traveled distance $L$. The paper then shows numerically, with a no-finite-size-effect algorithm that can track up to $10^9$ sites, that the soft and hard constraint models have overlapping light-cone curves for $d=1$ and $\\alpha>0.5$, including the critical crossover near $\\alpha=1$, and concludes that the two perspectives share identical butterfly light-cone scaling.","pith_inferences":["Beyond the paper: the same duality should extend to higher dimensions and to the dissipative, non-unitary settings the paper mentions as future directions.","Beyond the paper: the duality gives a microscopic reading of the FKPP noise term as coalescence noise, so finite-$N$ corrections to the light cone could be tested in power-law quantum simulators and should scale as $1/\\sqrt{N}$.","Beyond the paper: the fact that soft and hard constraint light cones overlap even at early times hints at a stronger coupling between the two processes, possibly a domination or coupling that would convert the numerical agreement into a proof."],"forward_implications":["The butterfly light-cone scalings of discrete population dynamics and the noisy long-range FKPP equation are identical for one-dimensional power-law interactions with $\\alpha>0.5$.","The dual FKPP equation has a Levy index $\\mu=2\\alpha-2$, rather than $2\\alpha-1$ as obtained by cutoff mean-field approximations, showing that noise and discreteness alter the effective transport.","The iterative integral equation $\\ell(t)^{2\\alpha}\\sim \\int_0^t \\ell(\\tau)\\ell(t-\\tau)\\,d\\tau$ organizes the stretched-exponential, power-law, and linear regimes, including the critical point $\\alpha=1.5$ and the crossover near $\\alpha=1$.","The higher-order analytic scaling near $\\alpha=1$ can serve as a benchmark for times where asymptotic scaling is not numerically or experimentally reachable.","The sparse Gillespie algorithm with no finite-size effect can probe light cones over $10^9$ sites and supplies a practical tool for both constraint models."],"supporting_citations":[{"why":"Supplies the Markov-generator-transpose duality construction that the paper adapts to birth-at-distance processes.","marker":"[30]"},{"why":"Gives the cluster-seeding iterative scaling argument and the leading-order butterfly light-cone scalings that the paper refines.","marker":"[31]"},{"why":"Provides the hard-constraint operator-growth population model whose scaling functions are the numerical benchmark.","marker":"[8]"},{"why":"Gives rigorous scaling results for the hard-constraint model that anchor the claimed stretched-exponential and power-law regimes.","marker":"[33]"},{"why":"Establishes the local short-range stochastic FKPP and population-dynamics consistency that this paper extends to long range.","marker":"[26]"},{"why":"Earlier mean-field cutoff treatment of the long-range FKPP equation whose Levy index $2\\alpha-1$ is contrasted with the exact dual index $2\\alpha-2$.","marker":"[34]"},{"why":"Supplemental derivations of the soft-constraint duality and of the scaling function near $\\alpha=1$ that underpin the analytic results.","marker":"[32]"}],"fun_headline_variants":["Duality links discrete and continuum chaos views","Soft-constraint model unifies operator growth","Noisy FKPP matches population dynamics exactly","Two chaos views share one butterfly light cone","Equilibrium population reveals chaos equivalence"],"cache_read_input_tokens":29056,"weakest_assumption_plain":"The chain of equivalence assumes that replacing the hard local population cap by a soft coalescence constraint with the same equilibrium population leaves the butterfly light-cone scaling unchanged; the paper proves the soft model is dual to FKPP but cannot construct a dual for the hard model.","fun_headline_variants_meta":{"raw":{"variants":["Duality links discrete and continuum chaos views","Soft-constraint model unifies operator growth","Noisy FKPP matches population dynamics exactly","Two chaos views share one butterfly light cone","Equilibrium population reveals chaos equivalence"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000233,"raw_usage":{"total_tokens":1477,"prompt_tokens":910,"completion_tokens":567,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":526,"completion_tokens_details":{"reasoning_tokens":502}},"tokens_in":526,"tokens_out":567,"duration_ms":5939,"temperature":1.0,"reasoning_tokens":502,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T22:44:42.869415+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Simulate hard- and soft-constraint models in one dimension at several $\\alpha$ values in $(0.5,1.5)$ with the reported no-finite-size algorithm, run until the iterative integral equation plateaus, and compare the inferred exponents $d\\log\\ell/d\\log t$; any difference outside numerical error would falsify the claimed equivalence.","supporting_citations":[{"cited_title":"Chen and T","cited_arxiv_id":null,"evidence_quote":"Supplies the Markov-generator-transpose duality construction that the paper adapts to birth-at-distance processes."},{"cited_title":"In this step, we do not yet select on which site the event starts","cited_arxiv_id":null,"evidence_quote":"Provides the hard-constraint operator-growth population model whose scaling functions are the numerical benchmark."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplemental derivations of the soft-constraint duality and of the scaling function near $\\alpha=1$ that underpin the analytic results."}],"review_version":1}