Relativistic Quantum Simulation under Periodic and Dirichlet Boundary Conditions: A First-Quantised Framework for Near-Term Devices

Jaewoo Joo; Jeongho Bang; Kyunghyun Baek; Timothy P. Spiller

REVIEW 3 major objections 3 minor 1 cited by

Relativistic Quantum Simulation under Periodic and Dirichlet Boundary Conditions: A First-Quantised Framework for Near-Term Devices

T0 review · 3 major / 3 minor · reviewed 2026-08-04 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read 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.

desk verdict 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. read the letter →

arxiv 2509.22579 v1 pith:7ZGGIQ2A submitted 2025-09-26 quant-ph

Jaewoo Joo , Timothy P. Spiller , Kyunghyun Baek , Jeongho Bang This is my paper · ORCID

classification quant-ph MSC 81P6881-08

keywords relativisticquantumsimulationfirstquantisationfinite-differencemethodperiodicboundaryconditionDirichlettranslationoperatoraddervariational

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

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.

What carries the argument

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

What would settle it

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.

Watch

Extended reading notes

Core claim

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

Load-bearing premise

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.

Editorial extensions

If this is right

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.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

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.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the 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.

What 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.

This 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.

Referee Report

3 major / 3 minor

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.

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 (3)

[Sec. II.B, Eqs. (21)-(22)] 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.
[Sec. III, Table I] 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.
[Sec. II.B, Eq. (15) and terminology] 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.

minor comments (3)

[Sec. II.A, Eq. (7)] 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.
[Sec. II.B, Eq. (20)] 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.
[General] 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.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity; the derivation chain is self-contained and no fitted parameter is relabeled as a prediction.

full rationale

The paper derives the PBC kinetic-energy expectation values directly from the finite-difference definition in Eq. (7) and the textbook perturbative expansion in Eq. (2). The DBC extension in Eqs. (15)-(26) is an algebraic manipulation of the same operators; no parameter is fitted to data and no external result is invoked as a uniqueness proof. The self-citations (Refs. [10], [11]) are used only to point to prior first-quantised circuits, while the relevant formulas are restated in the paper, so they are not load-bearing. The reader's flagged issue—the operator identities for E1 and E2 in Eqs. (21)-(22) appear to drop diagonal boundary terms—is a correctness/derivation error, not a circularity, because the paper does not define those operators in terms of the claimed final expectation values. The PBC portion is also self-contained. Therefore, by the hard rules, no circular step can be exhibited.

Assumptions & free parameters 1 free parameters · 3 assumptions · 0 invented entities

The central claim rests on standard quantum mechanics, a finite-difference approximation, and a boundary operator. No new physical entities are postulated. The DBC identities are where the derivation breaks.

free parameters (1)

L_m = lambda_m / delta_x = not fitted; user-selected scale
Combines Compton wavelength with grid spacing. The paper requires L_m sufficiently small for the perturbation expansion to converge; it is a chosen simulation parameter, not fitted to data.

assumptions (3)

domain assumption Relativistic kinetic energy expansion sqrt(1+p^2/(m^2 c^2))-1 converges and can be truncated
Section I, Eq. (2). Valid only when <p^2> << m^2 c^2, which the authors acknowledge.
domain assumption Finite-difference approximation p^2 = -(mcL_m)^2 (A^(1)-2I) is accurate enough at first-order precision
Section I, Eq. (7). Higher-order Laplacians are deferred to future work.
domain assumption Dirichlet boundary correction is exactly encoded by E0 = |N-1><0| + |0><N-1| with zero ghost-point boundary
Section II.B, Eqs. (15)-(16). The operator is asserted rather than derived from a physical Dirichlet constraint.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Relativistic Quantum Simulation under Periodic and Dirichlet Boundary Conditions: A First-Quantised Framework for Near-Term Devices." pith.science (2026). https://pith.science/paper/7ZGGIQ2A

@misc{pith2026250922579,
  author       = {Pith},
  title        = {Pith review of: Relativistic Quantum Simulation under Periodic and Dirichlet Boundary Conditions: A First-Quantised Framework for Near-Term Devices},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/7ZGGIQ2A}},
  note         = {Machine review of arXiv:2509.22579}
}

read the original abstract

We present a new recipe for relativistic quantum simulation using the first quantisation approach, under periodic (PBC) and Dirichlet (DBC) boundary conditions. The wavefunction is discretised across a finite grid represented by system qubits, and the squared momentum operator is expressed using the finite-difference method based on quantum translation operations. The relativistic kinetic energy is approximated through a perturbative expansion of the total kinetic Hamiltonian, incorporating higher-order momentum terms. The approach would allow variational optimisation of appropriate ansatz states to estimate both non-relativistic and relativistic ground-state energies on a quantum computer. This work offers a practical route to simulating relativistic effects on near-term quantum devices, supporting future developments in quantum physics and chemistry.

Figures

Figures reproduced from arXiv: 2509.22579 by the authors.

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

Quantum Gravity Simulation: Quantum simulation with a minimum length based on the generalised uncertainty principle
quant-ph 2026-07 reject novelty 4.0 of 10

Grid discretization in qubit simulations generates a generalized uncertainty principle; matching high-energy momentum to special relativity implies a mass-dependent minimum position uncertainty.

Reference graph

Works this paper leans on

19 extracted references · 2 linked inside Pith · cited by 1 Pith paper

[1]

Feynman,International Journal of Theoretical Physics,21, 467 (1982)

R.P. Feynman,International Journal of Theoretical Physics,21, 467 (1982)

1982
[2]

J. Li, B. A. Jones, and S. Kais,Science Advances9, eadg4576 (2023)

2023
[3]

Edery and P

A. Edery and P. Laporte,J. Phys. Commun.2, 025024 (2018)

2018
[4]

J. J. Sakurai and J. Napolitano, Modern Quantum Me- chanics,Cambridge: Cambridge University Press(2020)

2020
[5]

Fauseweh,Nat

B. Fauseweh,Nat. Comms,15, 2123 (2024)

2024
[6]

A. J. McCaskeyet al.,npj Quantum Inf.5, 99 (2019)

2019
[7]

H. H. S. Chanet al.,Science Advances9, eabo7484 (2023)

2023
[8]

Suet al.,PRX Quantum2, 040332 (2021)

Y. Suet al.,PRX Quantum2, 040332 (2021)

2021

Show all 19 references

[9]

Lubaschet al.,Phys

M. Lubaschet al.,Phys. Rev. A101, 010301(R) (2020)

2020
[10]

Joo, and H

J. Joo, and H. Moon, arXiv:2109.09216

arXiv
[11]

Joo, and T

J. Joo, and T. Spiller,New J. Phys.25, 083041 (2023)

2023
[12]

Vedral, A

V. Vedral, A. Barenco, and A. Ekert,Phys. Rev. A54, 147 (1996)

1996
[13]

J. Joo, M. Elliott, D. Oi, E. Ginossar, and T. P. Spiller, New J. Phys.18, 023028 (2016)

2016
[14]

P. C. S. Costa, S. Jordan, and A. Ostrander,Phys. Rev. A99, 012323 (2019)

2019
[15]

P. Overet. al.,Computers and Fluids288, 106508 (2025)

2025
[16]

Park and N

C.-Y. Park and N. Killoran,Quantum8, 1239 (2024)

2024
[17]

Zylberman and F

J. Zylberman and F. Debbasch, arXiv:2307.08384

arXiv
[18]

A. M. Childs, J.-P. Liu, and A. Ostrander,Quantum5, 574 (2021)

2021
[19]

Arseniev, D

B. Arseniev, D. Guskov, R. Sengupta, and I. Zacharov, Phys. Rev. A1114, 042625 (2025)

2025

Pith tools

Reviewed August 4, 2026 · model on record in the stance chip above.

[1] [1]

Feynman,International Journal of Theoretical Physics,21, 467 (1982)

R.P. Feynman,International Journal of Theoretical Physics,21, 467 (1982)

1982

[2] [2]

J. Li, B. A. Jones, and S. Kais,Science Advances9, eadg4576 (2023)

2023

[3] [3]

Edery and P

A. Edery and P. Laporte,J. Phys. Commun.2, 025024 (2018)

2018

[4] [4]

J. J. Sakurai and J. Napolitano, Modern Quantum Me- chanics,Cambridge: Cambridge University Press(2020)

2020

[5] [5]

Fauseweh,Nat

B. Fauseweh,Nat. Comms,15, 2123 (2024)

2024

[6] [6]

A. J. McCaskeyet al.,npj Quantum Inf.5, 99 (2019)

2019

[7] [7]

H. H. S. Chanet al.,Science Advances9, eabo7484 (2023)

2023

[8] [8]

Suet al.,PRX Quantum2, 040332 (2021)

Y. Suet al.,PRX Quantum2, 040332 (2021)

2021