{"id":"2d728722-132f-4dd7-83bc-ccf3af7e0676","arxiv_id":"2501.13069","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A quantum algorithm based on Haag-Ruelle theory and LCU proposes to prepare hadron wave packets from the vacuum in 3D lattice QCD, with a success probability that shrinks polynomially with lattice spacing, energy, and momentum width.","lead":"This paper designs a quantum-computing recipe for creating moving particles, like protons and pions, out of empty space in a lattice model of QCD. It is a step toward simulating hadron collisions in real time on quantum computers, which classical methods cannot do.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central claim hinges on unverified interpolator overlaps (Eq. 2): after the staggered mapping (Eqs. 72-76) with spin-isospin mixing, the operators are asserted, not shown, to couple to pions/nucleons; Section V.B concedes only O(a) control.","rationale":"The reader's weakest assumption is the same one I identify: unverified interpolator overlaps satisfying Eq. (2). I agree. After a careful reading, I did not find an internal inconsistency or a stronger objection: the Haag-Ruelle framework, the LCU implementation, and the qubit encodings are presented coherently; the assumptions about vacuum preparation and spectral knowledge are explicitly stated; and the extension to spin-1/2 particles is plausible. The genuinely load-bearing unproven step is the lattice interpolator's non-vanishing matrix element. The paper's own Section V.B caveat about spin-isospin mixing and O(a) errors makes this concern concrete rather than speculative. This does not require changing the reader's verdict: CONDITIONAL with moderate confidence is appropriate. I would not move to REJECT because the overlap is likely nonzero for standard lattice interpolators and the failure mode is an unproven assumption, not a demonstrated contradiction. The proposed concrete test would settle the concern by checking the overlaps in the free staggered theory first and then on a small interacting lattice with exact diagonalization.","tokens_in":25224,"tokens_out":10276,"duration_ms":121484,"concrete_test":"Diagonalize the free staggered Hamiltonian (68) on a finite lattice with gauge links set to identity and compute, on the one-particle shells, the overlaps <pi^+(k)|O_pi^+(0)|0> and <p(k)|O_p(0)|0> for the operators obtained from (90)-(92) via the mapping (70)-(76). Check that neither overlap vanishes for any k, and compare the two inequivalent spinor identifications of Section V.A: if one identification gives a zero overlap or the two give O(1) different overlaps, the interpolator is not a reliable creation operator. If the free overlaps are nonzero, repeat on a small (e.g., 2^3) Kogut-Susskind lattice with truncated SU(3) links using Lanczos exact diagonalization to verify the interacting overlap remains nonzero and obeys a power-law in a.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The method's efficiency is governed by the success probability (9), which is nonzero only if the lattice interpolating operators O_pi, O_p satisfy the condition <alpha|O|Omega> != 0 for the physical one-particle states. In Section V.B the authors import continuum interpolators (90)-(91), map them through the staggered identification (72)-(76), and insert gauge links to restore gauge invariance. This is where the central claim becomes conditional. The staggered theory mixes spin and isospin and has no symmetry that cleanly separates pion from nucleon quantum numbers on the lattice; the authors themselves state that 'identification of particles ... is possible only near the continuum limit, and up to O(a) errors.' No computation or argument is given that the resulting operator, after the gauge-link insertions and the O(a) redefinitions, has non-vanishing overlap with the target one-particle states. It could, in principle, couple dominantly to opposite-parity partners, additional tastes, or multi-hadron states. Moreover, the quoted LCU norms C_pi ~ 432/a^3 and C_p ~ 11500/a^4.5 are stated only approximately, without the full expansion, so the claimed polynomial scaling of rho is not independently verifiable from the manuscript. The efficiency statement in the abstract therefore inherits an unproven ingredient: without the overlap condition, the success probability could be zero or exponentially small in 1/a, not merely polynomially suppressed.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a strategy for preparing hadronic wave packets from the interacting vacuum in lattice gauge theories, based on Haag-Ruelle creation operators implemented via linear combinations of unitaries (LCU). It extends the authors' previous work [56] from scalar field theories to particles of general spin, with detailed treatments of U(1), SU(2), and SU(3) gauge theories in one spatial dimension, and a proposal for two-flavor QCD in three dimensions using Kogut-Susskind staggered fermions. The paper gives staggered-to-Dirac field identifications, gauge-invariant interpolating operators for pions and nucleons, qubit encodings (including an extension of generalized superfast encoding to odd fermionic operators), and estimates of LCU norms and success probabilities. The central claim is that wave-packet preparation succeeds with a probability that vanishes only polynomially in the lattice spacing, wave-packet energy, and momentum narrowness.","tokens_in":25447,"tokens_out":6669,"duration_ms":76392,"significance":"If the central claim holds, the paper would supply a concrete initial-state preparation protocol for real-time QCD scattering simulation on digital quantum computers, a genuinely useful step. The paper is commendably explicit about several technical building blocks: the spin-1/2 adaptation of the Haag-Ruelle construction, the LCU decomposition of interpolating operators in one-dimensional gauge theories, and the extension of generalized superfast encoding to odd operators. It also states its limitations honestly, including the staggered-fermion spin-isospin mixing. However, the efficiency claim is conditional on two ingredients that are not established in the manuscript: the nonzero overlap of the lattice interpolating operators with the intended one-particle states, and the exact values of the LCU norms entering the success probability. These are load-bearing for the abstract's polynomial-scaling statement, and they need concrete support before the claimed significance is realized.","major_comments":[{"comment":"Section V.B, Eqs. (90)-(91) and condition (2): the nonzero overlap of the staggered, gauge-link-inserted interpolators with the physical pion and nucleon one-particle states is asserted rather than established. Since the success probability in Eq. (9) is proportional to the squared norm of a†_ψ|Ω⟩, a vanishing or exponentially suppressed overlap would make the preparation fail or become exponentially hard, and the authors' own admission that staggered fermions mix spin and isospin makes this a real possibility. The manuscript should supply a concrete argument or calculation, for example a tree-level or lattice-perturbative evaluation of ⟨α|O_π,p|Ω⟩, or existing spectroscopy results for the exact staggered operators, before the polynomial-efficiency claim can be accepted.","section":"V.B"},{"comment":"The quoted LCU norms, C_π± ∼ 432/a^3 and C_p ∼ 11500/a^4.5, are stated only approximately and without the full coefficient sums required by Eq. (8). These numbers enter α in Eq. (9) and hence the success probability, so the claimed polynomial scaling of ρ in a and Ē is not independently verifiable from the manuscript. In particular, the conclusion that the staggered-fermion overheads are ∼10^4 for pions and ∼10^8 for nucleons does not follow from the displayed C values. A full expansion of the complete gauge-invariant interpolating operators, or a reproducible counting argument, is needed.","section":"V.B"},{"comment":"The complete gauge-invariant staggered interpolating operators are never written out: Eq. (92) gives only one baryonic link-insertion term, and the text states that the full expressions are 'extremely cumbersome'. Because these operators are the central object of the algorithm and determine both the overlap condition (2) and the LCU norms, the manuscript should provide them explicitly or give a precise constructive algorithm, rather than leaving their structure implicit.","section":"V.A/V.B"}],"minor_comments":[{"comment":"The symbol δ_p is introduced as the linear size of the momentum-space support of ψ̃, but the precise normalization and the dependence of the exponent d on the spatial dimension are not stated; please define δ_p in terms of ψ̃ so that the scaling can be checked.","section":"Section II, Eq. (10)"},{"comment":"The even-site summation with step 2 is introduced for one-dimensional theories, but it is not stated whether the same even-site restriction applies to the three-dimensional interpolators constructed in Section V; please clarify.","section":"Section III, Eq. (18)"},{"comment":"The ζ_αβ(x) strings for SU(3) are said to be 'analogous' to the SU(2) case, but they are not written out; since they are needed to reproduce the C values and to implement the operators, they should be defined explicitly in the text or in an appendix.","section":"Section IV, Eqs. (65)-(67)"},{"comment":"The phase factors A(y) and D(x,z) are introduced without explanation; a brief derivation of why this particular choice yields the continuum Hamiltonian in Eq. (77) would improve readability and allow readers to assess the O(a) ambiguities.","section":"Section V.A, Eqs. (70)-(76)"}],"recommendation":"major_revision","confidential_remarks":"The paper is largely an extension of the authors' own PRX Quantum paper [56], and the novelty claim of a 'first description' for QCD should be checked against Refs. [52] and [53], which also address hadron wave-packet preparation in gauge theories. The main risk is the unverified overlap condition after the staggered mapping; if the authors can add a concrete overlap calculation or reframe the claim as conditional on this assumption, a revised version may be publishable."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's the quick take: this is a serious theory proposal, not a numerical demonstration. If you care about initial-state preparation for real-time QCD simulation on quantum computers, it's worth reading and worth citing; just don't treat the efficiency claim as a proven theorem.\n\nWhat's new: The Haag-Ruelle-plus-LCU idea comes from their earlier scalar-field paper [56], but this paper does the real work of adapting it to spin-1/2 and higher-spin particles, gives explicit qubit-level interpolating operator decompositions for U(1), SU(2), SU(3) in 1D, extends the generalized superfast encoding to odd fermionic operators, and writes down pion and nucleon creation operators in 3D two-flavor QCD in Kogut-Susskind form. That last step—first explicit wave-packet protocol for hadrons in 3D QCD—is genuinely new. The paper is also unusually honest: it lists the three ingredients it needs (vacuum preparation, spectral knowledge, interpolating operators), and it explicitly concedes that staggered fermions mix spin and isospin, so particle identification holds only in the continuum limit up to O(a) errors, and that the overheads are large (~10^4 for pions, ~10^8 for nucleons).\n\nThe soft spots are real but not hidden. The main one is the overlap condition in Eq. (2): the whole success probability is proportional to |<alpha|O|Omega>|^2, and for the staggered-lattice operators with gauge links inserted, the paper asserts rather than shows that this overlap is nonzero for the physical pion/nucleon states. The stress-test note's worry that the operators could couple to opposite-parity partners, extra tastes, or multi-hadron states is legitimate. There is also no numerical check, and the quoted LCU norms C_pi ~ 432/a^3 and C_p ~ 11500/a^4.5 are stated only approximately, so the polynomial scaling of the success probability isn't fully verifiable from the manuscript. The formulas for N and rho come from the previous paper [56]; that's not circular—it's a published source—but it does mean the independence of the efficiency claim rests partly on trust.\n\nWhere does that leave us? The central construction is coherent, the caveats are acknowledged, and the missing piece—non-vanishing overlap—is a standard assumption in lattice hadron spectroscopy, usually checked by computing two-point correlators. This paper doesn't do that check, so it's a conditional result. But it's exactly the kind of work that a serious referee should engage with, and the odd-operator GSE extension alone is a useful technical contribution. I'd send it to review, and I'd cite it if I were writing about wave-packet preparation or superfast encodings. It's not a paper you should desk-reject; it's a paper where you ask the authors to either prove or numerically support the overlap assumption, and to give the full C expansions.","headline":"A serious theory proposal for hadron wave-packet preparation via Haag-Ruelle and LCU, genuinely new in scope but with the load-bearing overlap assumption unproven; worth a serious referee.","tokens_in":26065,"tokens_out":3399,"would_cite":true,"duration_ms":30334,"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 claims that wave packets of pions and nucleons can be created from the interacting vacuum of three-dimensional two-flavor QCD using Haag-Ruelle creation operators implemented with LCU, with a success probability that vanishes…","keywords":["wave packet preparation","quantum simulation of QCD","Haag-Ruelle scattering theory","linear combination of unitaries","lattice gauge theory","staggered fermions","initial state preparation","hadron scattering on quantum computers"],"falsifier":"On a small three-dimensional lattice with two-flavor staggered fermions, compute the overlaps $\\langle \\alpha|\\hat O_\\pi(x)|\\Omega\\rangle$ and $\\langle \\alpha|\\hat O_p(x)|\\Omega\\rangle$ for the pion and proton interpolators of Section V.B: if the overlap is zero, or decays exponentially in $1/a$ instead of polynomially, then the creation operators cannot prepare the claimed wave packets. A cheaper proxy is to check the two-point correlation function for a one-particle pole with the expected quantum numbers on a modest-size lattice.","tokens_in":24966,"feed_emoji":"⚛️","tokens_out":9180,"duration_ms":84544,"temperature":0.7,"pith_summary":"Real-time scattering simulations in quantum field theory need initial states made of two wave packets of stable particles moving on the vacuum, and for confined theories like QCD these states are hard to prepare. This paper extends the Haag-Ruelle scattering-theory approach to lattice gauge theories and, for the first time, to three-dimensional two-flavor QCD, giving explicit creation operators for pions and nucleons. The operators act on the interacting vacuum and are implemented on a digital quantum computer with the linear-combination-of-unitaries technique; the preparation succeeds upon measuring an ancilla with probability that vanishes polynomially, not exponentially, in the lattice spacing, the wave-packet energy, and the momentum narrowness. If the required vacuum preparation and spectral knowledge are available, this supplies the missing initial-state ingredient for simulating hadron collisions step by step.","feed_headline":"New recipe creates hadron wave packets from the QCD vacuum","feed_subtitle":"Only polynomially vanishing success probability keeps real-time hadron collisions within reach.","key_machinery":"The load-bearing object is the Haag-Ruelle creation operator $\\hat a^\\dagger_\\psi$, a spacetime-smeared interpolating operator $\\hat O(x)$ satisfying $\\langle \\alpha|\\hat O|\\Omega\\rangle \\ne 0$ for the one-particle state to be created. The interpolator selects the particle species by symmetry quantum numbers, while the Schwartz function $\\psi$, a smooth rapidly decaying function whose Fourier transform has support only on the corresponding mass hyperboloid, cuts off multiparticle contamination. On the lattice the creation operator becomes a finite sum over even lattice translations and time samples, implemented by LCU with an ancillary register; the same LCU decomposition handles the non-unitary interpolators, whose norm constants $C$ enter the success probability. The other essential piece is the continuum-to-staggered-fermion identification of Eqs. (72)--(76), which maps pion and nucleon interpolators into gauge-invariant products of staggered fields and link operators, and the generalized superfast encoding, extended to odd operators, which maps those products to qubit strings.","core_discovery":"The central claim is that wave packets of composite hadrons can be produced directly from the interacting vacuum of lattice QCD, without adiabatically deforming free-theory states. The construction uses Haag-Ruelle creation operators $\\hat a^\\dagger_\\psi = \\sum_x a^d \\int dt\\, \\psi(t,x) e^{iHt} \\hat O(x) e^{-iHt}$, where $\\hat O$ is an interpolating operator carrying the particle's quantum numbers and $\\psi$ is smeared to have momentum support only on the one-particle mass hyperboloid. For two-flavor QCD in three dimensions the paper gives explicit interpolators for pions and nucleons, made gauge invariant by inserting link operators, translated from staggered fermions into Dirac spinors through the identification of Eqs. (72)--(76), and encoded as qubit operators using the generalized superfast encoding extended to odd fermionic operators. The resulting success probability is estimated as $\\rho \\sim Z/(C^2 \\delta_p^d \\bar E)$ with spin-one-half modifications, vanishing polynomially; the main quantitative cost appears in the constants $C$, with pion interpolators giving $C\\sim 432/a^3$ and proton interpolators $C\\sim 11500/a^{4.5}$, and the paper notes that particle identification in staggered fermions is reliable only near the continuum limit up to $O(a)$ errors.","pith_inferences":["Editorial inference: the decisive numerical test of the method is to compute $\\langle \\alpha|\\hat O|\\Omega\\rangle$ for the staggered pion and nucleon interpolators on small lattices; if these overlaps vanish or decay exponentially toward the continuum, the polynomial-success-probability claim fails even though the formalism is intact.","Editorial inference: because the cost is dominated by repeated vacuum preparation, the practical feasibility of this initial-state strategy will be set by the quality of vacuum preparation and spectral inputs, so progress on those two problems translates directly into progress on scattering simulation.","Editorial inference: the large constants for nucleons ($C\\sim 11500/a^{4.5}$) suggest that a first demonstration is more likely with mesons or with the one-dimensional SU(3) toy model, where the same building blocks appear with much smaller $C$ values.","Editorial inference: if staggered-fermion spin-isospin mixing proves too severe near current lattice sizes, porting the interpolators to a Wilson-fermion or loop-string-hadron formulation, which the paper says is straightforward, would cut the overhead by orders of magnitude."],"forward_implications":["Pion-pion and proton-proton collisions could be simulated step by step, not just as S-matrix amplitudes, once two such wave packets are placed on the vacuum.","The success probability vanishes only polynomially in lattice spacing, wave-packet energy, and momentum narrowness, so the preparation cost stays manageable in principle rather than being exponentially suppressed.","Preparing $n$ wave packets requires on average $O(1/\\rho^n)$ repetitions of vacuum preparation and creation, with $n=2$ the relevant case; amplitude amplification can quadratically improve this.","The same construction extends to other particles and other lattice formulations; with Wilson fermions the authors estimate the interpolator constants would be $C_\\pi = 12/a^3$ and $C_p = 24/a^{4.5}$, reducing the overhead.","The only theoretical limitation stated is the presence of massless particles, since an isolated mass shell is needed."],"supporting_citations":[{"why":"Supplies the Haag-Ruelle wave-packet creation method for scalar theories that this paper extends to gauge theories and QCD.","marker":"[56]"},{"why":"Provides the original Haag-Ruelle scattering theory that justifies creation operators of the form $\\hat a^\\dagger_\\psi = \\int dx\\, \\psi \\hat O$.","marker":"[38, 39]"},{"why":"Formalizes Haag-Ruelle scattering for gapped lattice systems, giving the lattice-theoretic foundation for applying the construction on a spatial lattice.","marker":"[64]"},{"why":"Provides the LCU technique used to implement the non-unitary creation operator on a quantum computer.","marker":"[67]"},{"why":"Supplies the Kogut-Susskind Hamiltonians and the SU(2)/SU(3) gauge-link decompositions used to express interpolators as linear combinations of unitaries.","marker":"[70]"},{"why":"Give the staggered-fermion identification of continuum Dirac fields used to map pion and nucleon interpolators onto the lattice.","marker":"[73, 74]"},{"why":"Introduces the generalized superfast encoding of fermions onto qubits, which the paper extends to odd operators for baryons.","marker":"[75]"},{"why":"Uses generalized superfast encoding for lattice gauge theory simulation, providing the resource context for the encoded operators.","marker":"[76]"},{"why":"Provide the traditional-lattice pion and nucleon interpolating operators that are imported and rewritten in staggered-fermion form.","marker":"[77-80]"}],"fun_headline_variants":["Quantum computers create QCD wave packets in 3D for first time","Haag-Ruelle method prepares hadron wave packets from vacuum","Direct vacuum-to-wave-packet route for lattice QCD","3D QCD wave packet creation enabled on quantum hardware","Wave packets from QCD vacuum: polynomial cost recipe"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that the lattice operators written down for pions and nucleons really do overlap with the corresponding one-particle states (equation (2)); if that overlap is zero or strongly suppressed after the staggered-fermion mapping, the created wave packets are not the intended hadrons.","fun_headline_variants_meta":{"raw":{"variants":["Quantum computers create QCD wave packets in 3D for first time","Haag-Ruelle method prepares hadron wave packets from vacuum","Direct vacuum-to-wave-packet route for lattice QCD","3D QCD wave packet creation enabled on quantum hardware","Wave packets from QCD vacuum: polynomial cost recipe"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000168,"raw_usage":{"total_tokens":1290,"prompt_tokens":1003,"completion_tokens":287,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":619,"completion_tokens_details":{"reasoning_tokens":203}},"tokens_in":619,"tokens_out":287,"duration_ms":3415,"temperature":1.0,"reasoning_tokens":203,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T16:28:33.799129+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"On a small three-dimensional lattice with two-flavor staggered fermions, compute the overlaps $\\langle \\alpha|\\hat O_\\pi(x)|\\Omega\\rangle$ and $\\langle \\alpha|\\hat O_p(x)|\\Omega\\rangle$ for the pion and proton interpolators of Section V.B: if the overlap is zero, or decays exponentially in $1/a$ instead of polynomially, then the creation operators cannot prepare the claimed wave packets. A cheaper proxy is to check the two-point correlation function for a one-particle pole with the expected quantum numbers on a modest-size lattice.","supporting_citations":[],"review_version":1}