{"id":"5acb01f1-e593-49bd-8a8e-54516ce908ed","arxiv_id":"2509.22579","paper_version":1,"verdict":"REJECT","confidence":"HIGH","novelty_score":4.0,"correctness_risk":"high","formal_verification":"none","parameter_count":1,"one_line_summary":"A first-quantized quantum simulation framework approximates relativistic kinetic energy via a perturbative expansion of finite-difference momentum operators under periodic and Dirichlet boundary conditions.","lead":"This paper develops a circuit-based recipe for simulating relativistic kinetic energy corrections on a discretized quantum grid, handling both periodic and Dirichlet boundary conditions. It expresses the relativistic Hamiltonian as combinations of quantum translation operator expectations, intended for variational ground-state energy estimation on near-term quantum computers.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"DBC correction is algebraically unsound: the claimed identities for E1 and E2 in Eq. (20) drop diagonal boundary terms, so Eqs. (20)-(26) do not give the stated finite-difference ⟨p^4⟩_D.","rationale":"The reader's weakest assumption identifies precisely the same load-bearing algebraic error in the Dirichlet-boundary section: the definitions/identities for E1 and E2 are false, and the dropped boundary-diagonal terms alter ⟨p^4⟩_D and hence the DBC relativistic kinetic energy in Eq. (26). This is not a matter of convention or approximation; it is a direct miscalculation of the operator whose expectation value is supposed to be estimated. The PBC part of the paper is more secure, but the abstract and conclusions advertise DBC as a key contribution, so the central claim is damaged. I recommend keeping the reader's REJECT verdict: the algebraic error is concrete, reproducible, and load-bearing, and no numerical results are provided that could accidentally compensate for it. My own reading confirms the reader's finding rather than introducing a new objection.","tokens_in":7722,"tokens_out":3839,"duration_ms":29249,"concrete_test":"Take N=4 (L=2), define A as the cyclic shift and E0=|3⟩⟨0|+|0⟩⟨3|. Compute M=(A+A†−2I−E0)^2 explicitly and compare its expectation on a uniform state |ψ⟩=(1/2)(|0⟩+|1⟩+|2⟩+|3⟩) with the value obtained from Eq. (20) using the paper's E1=|1⟩⟨3|+|3⟩⟨1| and E2=|0⟩⟨2|+|2⟩⟨0|. If the two values differ, the discrepancy is exactly the dropped diagonal boundary terms, confirming that Eq. (26) does not compute the stated DBC kinetic energy. A purely symbolic check also suffices: verify that A E0 + E0 A† − (|1⟩⟨N−1|+|N−1⟩⟨1|) = 2|0⟩⟨0|.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing step is the Dirichlet-boundary expansion of Eq. (20). With A = Σ|j+1⟩⟨j| (cyclic shift) and E0=|N−1⟩⟨0|+|0⟩⟨N−1|, direct algebra gives A E0 + E0 A† = 2|0⟩⟨0| + |1⟩⟨N−1| + |N−1⟩⟨1|, not |1⟩⟨N−1| + |N−1⟩⟨1| as claimed in Eq. (21). Similarly, E0 A + A† E0 = |N−1⟩⟨0| + |0⟩⟨N−2| + |N−2⟩⟨0| + |N−1⟩⟨N−1|, not |0⟩⟨N−2| + |N−2⟩⟨0| as claimed in Eq. (22). The dropped diagonal terms contribute to the expectation value of (A(1) − 2I − E0)^2, so Eq. (20)'s expression for ⟨p^4⟩_D misses, for example, a 2(mcL_m)^4|c0|^2 contribution from the |0⟩⟨0| term in E1. For a uniform four-point state this changes the fourth-momentum term by a factor of three. Since Eq. (26) is explicitly built from these identities for E1 and E2, the DBC relativistic kinetic energy formula is not the expectation value of the finite-difference operator the paper defines. The PBC construction appears internally consistent, but the new DBC contribution—the paper's distinctive extension—rests on this false algebra. No numerical validation or independent derivation is supplied to mitigate the algebraic error.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a first-quantized quantum algorithm for estimating relativistic kinetic-energy corrections via the perturbative expansion of √(m²c⁴+p²c²). The wavefunction is discretized on an L-qubit grid, and the squared momentum operator is replaced by a finite-difference operator built from cyclic shift operators. For periodic boundary conditions, the even powers of momentum are expressed through expectation values of A^(l) = A_x^l + (A_x†)^l, leading to a decomposition of the relativistic kinetic energy with precomputable coefficients β_l. For Dirichlet boundary conditions, an additional operator E0 is introduced to remove the periodic wrap-around, and the paper derives correction terms E1 and E2, together with a circuit protocol for estimating ⟨p⁴⟩_D. The final variational recipe combines kinetic and potential energy estimates to obtain the total ground-state energy.","tokens_in":8172,"tokens_out":22045,"duration_ms":144298,"significance":"If correct, the manuscript would offer a simple, near-term-friendly route to relativistic kinetic-energy corrections in first-quantized quantum simulation. The PBC construction is transparent: the finite-difference operator is standard, the coefficients are derived rather than fitted, and the Hadamard-test circuit for ⟨A^(l)⟩ is plausible. The DBC extension, however, is the paper's distinctive contribution, and its central algebraic identities are incorrect. The error is load-bearing: Eqs. (20)–(26) do not compute the expectation value of the finite-difference operator defined in Eq. (15), and the Table I measurement protocol does not estimate the quantities it claims. The paper contains no numerical example or independent verification that would catch this error.","major_comments":[{"comment":"The operator identities for E1 and E2 are false. With A = Σ|j+1⟩⟨j| (cyclic) and E0 = |N−1⟩⟨0|+|0⟩⟨N−1|, direct algebra gives E1 = A E0 + E0 A† = 2|0⟩⟨0| + |1⟩⟨N−1| + |N−1⟩⟨1|, and E2 = E0 A + A† E0 = |0⟩⟨N−2| + |N−2⟩⟨0| + 2|N−1⟩⟨N−1|. The terms 2|0⟩⟨0| and 2|N−1⟩⟨N−1| are missing from Eqs. (21)–(22). Consequently, the expectation values in Eqs. (23)–(24) should read ⟨E1⟩ = 2P0 + c1*c_{N−1}+c_{N−1}*c1 and ⟨E2⟩ = 2P_{N−1} + c0*c_{N−2}+c_{N−2}*c0. This error is not cosmetic: for a uniform four-point state, Eq. (20) with the printed forms gives ⟨p⁴⟩_D = 1.5 (mcL_m)^4, whereas the true expectation of (A(1)−2I−E0)² is 0.5 (mcL_m)^4. Since Eq. (26) is explicitly built from these identities, the DBC relativistic kinetic-energy formula is not the expectation value of the finite-difference operator the paper defines.","section":"Sec. II.B, Eqs. (21)-(22)"},{"comment":"The reference-state protocol in Table I is based on the false E1/E2 identities. The measurement Ps − Pf − Pg estimates only the off-diagonal overlaps such as c1*c_{N−1}+c_{N−1}*c1; it cannot capture the new diagonal terms 2P0 and 2P_{N−1} that the true E1 and E2 contain. Thus the circuit protocol does not produce the DBC kinetic-energy correction in Eq. (26) as stated. A corrected protocol must either add boundary probability measurements or redefine the Fd set so that the diagonal boundary contributions are explicitly included. This is fixable, but it changes the core of Section III.","section":"Sec. III, Table I"},{"comment":"The boundary condition implemented by D = A(1)−2I−E0 is not standard Dirichlet. This operator removes only the cyclic coupling between site 0 and site N−1, leaving boundary sites as dynamical variables with self-energy −2 and a single neighbor. It does not enforce ψ(0)=ψ(1)=0. If the authors intend a genuine Dirichlet problem, the finite-difference operator must be different (e.g., boundary rows/columns removed or constrained). If they intend an open-boundary or nonperiodic condition, the terminology should be changed throughout the abstract and text. This distinction is load-bearing for the physical claim of the paper.","section":"Sec. II.B, Eq. (15) and terminology"}],"minor_comments":[{"comment":"The dimensionless parameter L_m = λ_m/δx is dimensionally inconsistent unless the physical domain length R is set to 1. Since the physical coordinate is x' = R x, the actual grid spacing is Rδx, and the correct ratio is λ_m/(Rδx). Please clarify the normalization of R or define L_m accordingly.","section":"Sec. II.A, Eq. (7)"},{"comment":"The formal operator identity in Eq. (20) is correct if E1 and E2 are kept as exact operators. The problem is the immediate substitution of the incorrect closed forms. Please present the exact expansion and the corrected expectation values together, so that the source of the error is visible.","section":"Sec. II.B, Eq. (20)"},{"comment":"The manuscript would benefit from a small numerical example, e.g., a uniform four-point state, to verify the PBC and DBC formulas. Such a check would have exposed the factor-of-three discrepancy in ⟨p⁴⟩_D and would greatly increase confidence in the corrected version.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The PBC section is sound but incremental; the DBC part is the claimed novelty and is currently algebraically incorrect. I am not recommending immediate rejection because the error is localized and can in principle be repaired by using the exact E1/E2 identities and updating the expectation-value/circuit protocol. However, acceptance should be contingent on a complete re-derivation of Section II.B and III with corrected operators, plus a numerical check on a small state. If the DBC framework cannot be repaired, the paper would not meet the bar for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper is a first-quantized recipe for relativistic kinetic energy on a lattice, with circuits for periodic and Dirichlet boundary conditions. The periodic-boundary part is a clean, correct extension of the earlier first-quantized framework. The Dirichlet part, which is the paper's distinctive contribution, has a load-bearing algebraic error: the claimed identities for E1 and E2 in Eqs. (21)-(22) drop diagonal boundary terms. I checked directly: with A the cyclic shift and E0 as defined, A E0 + E0 A† equals 2|0⟩⟨0| + |1⟩⟨N−1| + |N−1⟩⟨1|, not the two-term expression in Eq. (21). Similarly E0 A + A† E0 has a 2|N−1⟩⟨N−1| term missing from Eq. (22). These dropped terms contribute to ⟨p^4⟩_D, so Eq. (20) and the final kinetic formula Eq. (26) do not give the expectation value of the finite-difference operator they define. For a uniform four-point state the fourth-momentum term is overestimated by a factor of three. The error cannot be waved away by assuming the wavefunction vanishes at the boundaries: the paper explicitly allows |ψ(0)⟩ ≠ |ψ(1)⟩ and the ansatz states are unrestricted. This is not a nitpick: the DBC section is the paper's reason to exist, and it rests on this algebra.\n\nWhat the paper does well: the PBC decomposition into translation-operator expectation values is correctly derived, and the circuit for estimating those expectation values is standard and plausible. The coefficients β_l are precomputable, which is nice. There is no numerical or hardware demonstration, so even the correct part is only a recipe, not a validated method. The error is fixable—add the missing boundary projectors and the formulas would change—but as written the DBC result is wrong.\n\nThis paper deserves a serious referee because the mistake is specific and the PBC part may be useful to some readers, but my recommendation is reject in current form. If you send it to review, the referee should be asked to verify Eqs. (21)-(22) explicitly. I would not cite it in its present state.","headline":"The DBC section has a concrete algebraic error that breaks the paper's main new result; the PBC part is correct but a straightforward extension.","tokens_in":8608,"tokens_out":16138,"would_cite":false,"duration_ms":104450,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["81P68","81-08"],"pacs":[],"model":"deepseek-v4-flash","headline":"A first-quantised, grid-based quantum circuit recipe estimates relativistic ground-state energies as perturbative sums of translation-operator expectation values, under periodic and Dirichlet boundary conditions.","keywords":["relativistic quantum simulation","first quantisation","finite-difference method","periodic boundary condition","Dirichlet boundary condition","translation operator","quantum adder","variational quantum simulation"],"falsifier":"Evaluate the operators E1 = A E0 + E0 A^† and E2 = E0 A + A^† E0 on the boundary basis states |0⟩ and |N−1⟩ for a small grid (for example, L=3 qubits): the exact expressions contain additional diagonal terms 2|0⟩⟨0| and 2|N−1⟩⟨N−1|. Comparing ⟨p^4⟩_D computed from Eq. (20) with the exact finite-difference expectation value of (A^(1) − 2I − E0)^2 on a state with nonzero c0 or c_{N−1} will show a discrepancy proportional to |c0|^2 + |c_{N−1}|^2, settling whether the DBC circuit estimates the intended discretised Hamiltonian.","tokens_in":7651,"feed_emoji":"⚛️","tokens_out":5152,"duration_ms":120800,"temperature":0.7,"pith_summary":"This paper tries to establish that relativistic ground-state energies of a single particle in one dimension can be estimated on a quantum computer without second quantisation. The key move is to replace the square-root relativistic kinetic Hamiltonian by a perturbative expansion in powers of p^2, and to express each power as a finite-difference combination of translation-operator expectation values on an L-qubit position grid. The authors give explicit formulas and circuits for periodic and Dirichlet boundary conditions, and propose a variational minimisation of kinetic plus potential energy. If correct, this offers a concrete route to including relativistic corrections in near-term quantum simulations of simple quantum systems.","feed_headline":"Translation gates estimate relativistic ground-state energies","feed_subtitle":"A finite grid and momentum expansion let near-term devices test relativistic effects under two boundary conditions.","key_machinery":"The load-bearing identity is the finite-difference representation of squared momentum on an L-qubit grid: p̂^2 = −(mc L_m)^2 (A^(1) − 2I), where A^(1) = A_x + A_x^† is the sum of the quantum adder and subtractor translation gates, and L_m = λ_m/δx is the reduced Compton wavelength divided by the grid spacing. The relativistic kinetic energy is then expanded as a power series in p̂^2 with coefficients α_l from the Klein-Gordon expansion. The circuits measure the real parts of controlled-translation expectation values with a single control qubit, and the DBC boundary terms are extracted through reference-state overlap measurements. Variational ansatz states with tunable parameters complete the","core_discovery":"The central claim is that the relativistic kinetic energy can be computed as a linear combination of expectation values of translation operators A^(l) = A_x^l + (A_x^†)^l acting on a discretised L-qubit wavefunction. Under periodic boundary conditions, ⟨p^{2l}⟩_P = (−(mc L_m)^2)^l ⟨(A^(1) − 2I)^l⟩, so the first relativistic correction only requires ⟨A^(1)⟩ and ⟨A^(2)⟩. Under Dirichlet boundary conditions, the same quantities are expressed through the same translation operators plus additional edge operators E0, E1, E2, and E0^2, with a table of reference states that allow each expectation value to be measured. The total ground-state energy is obtained by variational optimisation of the sum o","pith_inferences":["An implication the authors leave implicit is that the same translation-based finite-difference construction could be extended to higher-order Laplacian stencils: the ⟨A^(2)⟩ circuit already supplies the data needed for a second-order-precision second derivative, potentially improving accuracy without additional qubits.","A reader checking the DBC derivation will find that the identities for E1 and E2, as written in Eqs. (21) and (22), omit diagonal boundary terms that appear in the exact operator products; correcting this would change the boundary contribution to ⟨p^4⟩_D and therefore to the relativistic energy estimate in Eq. (26).","The paper assumes real probability amplitudes for the simplified circuits; extending to complex amplitudes would require measuring both real and imaginary parts of the controlled-translation expectation values, doubling the number of circuit runs on near-term hardware.","A natural testable extension is to apply the PBC part of the recipe to a known exactly solvable potential (for example, a harmonic trap) and compare the variational relativistic ground-state estimate against perturbation theory in a regime where L_m is small."],"forward_implications":["For periodic boundary conditions, estimating the relativistic kinetic energy to l-th order requires only the translation expectation values ⟨A^(1)⟩, ⟨A^(2)⟩, ..., ⟨A^(l)⟩, all of which are accessible from the same class of controlled-translation circuits.","For Dirichlet boundary conditions, the nonrelativistic kinetic energy can be obtained by reusing the periodic-boundary result plus a single extra expectation value ⟨E0⟩, so existing PBC data can be recycled.","The first-order relativistic Dirichlet correction (through p^4) is assembled from the PBC relativistic kinetic energy plus four edge expectation values: ⟨E0⟩, ⟨E1⟩, ⟨E2⟩, and ⟨E0^2⟩.","If the method works, near-term variational quantum simulators could compute relativistic corrections for one-dimensional model systems without needing expensive second-quantised mappings.","Adding one qubit halves the grid spacing and improves the finite-difference accuracy, provided the perturbative parameter L_m stays small enough for the expansion to converge."],"fun_headline_variants":["Translation gates measure relativistic kinetic energy on quantum grids","Relativistic energies from translation operators under two boundary conditions","First-quantized recipe for relativistic simulation on near-term quantum devices","Translation gates estimate relativistic energies with a finite-grid expansion","Relativistic ground-state energies via translation operator expectations"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The Dirichlet boundary-condition derivation depends on the operator identities E1 = |1⟩⟨N−1| + |N−1⟩⟨1| and E2 = |0⟩⟨N−2| + |N−2⟩⟨0| following from their definitions, but the true operators contain extra diagonal boundary terms, so the paper's p^4 Dirichlet formula does not actually match the finite-difference operator it claims to implement.","fun_headline_variants_meta":{"raw":{"variants":["Translation gates measure relativistic kinetic energy on quantum grids","Relativistic energies from translation operators under two boundary conditions","First-quantized recipe for relativistic simulation on near-term quantum devices","Translation gates estimate relativistic energies with a finite-grid expansion","Relativistic ground-state energies via translation operator expectations"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000768,"raw_usage":{"total_tokens":3206,"prompt_tokens":675,"completion_tokens":2531,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":419,"completion_tokens_details":{"reasoning_tokens":2466}},"tokens_in":419,"tokens_out":2531,"duration_ms":97090,"temperature":1.0,"reasoning_tokens":2466,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T14:49:44.986400+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Evaluate the operators E1 = A E0 + E0 A^† and E2 = E0 A + A^† E0 on the boundary basis states |0⟩ and |N−1⟩ for a small grid (for example, L=3 qubits): the exact expressions contain additional diagonal terms 2|0⟩⟨0| and 2|N−1⟩⟨N−1|. Comparing ⟨p^4⟩_D computed from Eq. (20) with the exact finite-difference expectation value of (A^(1) − 2I − E0)^2 on a state with nonzero c0 or c_{N−1} will show a discrepancy proportional to |c0|^2 + |c_{N−1}|^2, settling whether the DBC circuit estimates the intended discretised Hamiltonian.","supporting_citations":[],"review_version":1}