Systematically Localizable Operators for Quantum Simulations of Quantum Field Theories

I. INTRODUCTION

• Quantum computers are proposed for accessing real-time dynamics and high-density systems beyond classical polynomial-time calculations, with relevance to nuclear, high-energy, and basic-energy science.The paper frames quantum simulation as a route toward observables that remain difficult for classical approaches.Source: §I · p. 3
• Scalar and pseudoscalar fields are central to nuclear and high-energy physics, motivating efficient quantum algorithms for preparing and evolving lattice field states.The discussion connects pion fields and other scalar degrees of freedom to nuclear interactions and quantum simulation.Source: §I · pp. 3–4
• High-fidelity ground-state preparation is essential because imperfect initialization can generate background particles that obscure subsequent scattering and dynamical measurements.The stated examples include glueballs, pions, and baryon-antibaryon systems in QCD calculations.Source: §I · p. 4

II. LATTICE SCALAR FIELD THEORY AND ITS CORRELATIONS

• Free lattice scalar theory reduces to coupled harmonic oscillators whose spatial correlations and covariance off-diagonal elements decay exponentially, enabling polynomial-in-volume initialization circuits.The ground state is diagonal in momentum space but becomes a correlated Gaussian wavefunction in position space.Source: §I · p. 4

I. INTRODUCTION

• The paper designs localizable circuits for free and interacting lattice scalar theories by organizing controlled operations according to spatial separation and correlation structure.The construction also supports systematic spatial truncations and quantification of the resulting wavefunction-fidelity loss.Source: §I · pp. 4–5

A. Continuous Lattice Scalar Field in One Spatial Dimension

• The lattice formulation discretizes space, uses finite differences for gradients, and diagonalizes the coupled system into momentum-mode harmonic oscillators.Periodic boundary conditions and dimensionless lattice variables provide the setup for the one-dimensional analysis.Source: §II · p. 5, A. Continuous Lattice Scalar Field in… · p. 6

1. Two-Point Correlation Functions and Lattice Artifacts

• The ground-state field correlator is related to the lattice K matrix and has continuum expressions involving modified Bessel functions K0(x) and K1(x).Source: §1 · pp. 7–8
• The lattice two-point function approaches infinite-volume, continuum, and large-separation expressions while exhibiting finite-volume and finite-spacing effects.For ˆm = 0.3 on N = 30 sites, periodic-boundary symmetry appears about the lattice midpoint.Source: §1 · p. 7, Figure 1 · p. 7

2. Entanglement Entropy and Mutual Information

• Mutual information between two lattice sites is computed from reduced-density-matrix partitions of the K matrix and decreases exponentially with site separation.The analysis makes the reduced state approach a tensor product as the sites become infinitely separated.Source: §2 · pp. 8–9
• At large separation, the inter-site parameters γ12 and β12 scale as e^(-ˆmj)/j^p, causing the mutual information to scale as e^(-2ˆmj)/j^p′.Source: §2 · p. 9
• Mutual information falls exponentially faster than the two-point correlation function, with exponents −2ˆm versus −ˆm.Source: §2 · p. 9

• Negativity is the sum of negative eigenvalues of the partially transposed density matrix and therefore measures purely quantum correlations.Source: §3 · p. 9
• Negativity vanishes when the two sites decouple, K12 = 0.Source: §3 · p. 9
• Negativity is localized to the interaction’s extent rather than falling exponentially with spatial separation.Source: §3 · p. 9
• Nearest-neighbor finite differences produce negativities that vanish beyond nearest neighbors, apart from rapidly vanishing finite-size effects.Source: §3 · p. 10
• Changing or smearing the finite-difference operator shifts the nonzero negativities to the corresponding site separations, while mutual information remains exponentially decaying.Source: §3 · p. 10, Figure 2 · p. 10
• Controlled rotations localized according to classical and quantum correlations are the focus of the circuit construction that follows.Source: §3 · p. 11

• Each spatial site is mapped to an nQ-qubit register with Hilbert-space dimension ns = 2^nQ.Source: B. Digitized Lattice Scalar Field · p. 11
• The digitization uses a symmetric field-eigenstate basis spanning discretized values between −φmax and φmax.Source: B. Digitized Lattice Scalar Field · p. 11
• Starting from an equal-weight lattice wavefunction, a normalized non-unitary construction reproduces the continuum wavefunction on the digitized Hilbert space.Source: B. Digitized Lattice Scalar Field · p. 11
• The resulting wavefunction is real and invariant under φ → −φ, reducing the independent real amplitudes by reflection symmetry.Source: B. Digitized Lattice Scalar Field · p. 11

III. QUANTUM CIRCUITS FOR GROUND STATE WAVEFUNCTIONS

• The restructuring begins by at most doubling quantum operations while making long-distance circuit elements removable.Source: §III · p. 11
• For lattice scalar ground states, truncating the localized circuit yields exponentially improvable wavefunction fidelity because two-point correlations decay exponentially with separation.Source: §III · pp. 11–12

A. Controlled Rotations: the θ-Angles

• The generic circuit for an arbitrary real, positive wavefunction uses 2^nQ − 1 rotations and 2^nQ − 2 CNOT gates.Source: A. Controlled Rotations: the θ-Angles · p. 12
• At level ℓ, the θℓ vector contains 2^ℓ angles indexed by the binary values of the control qubits.Source: A. Controlled Rotations: the θ-Angles · p. 12
• The rotation angles are obtained from ratios of wavefunction amplitudes across binary-defined bipartitions.Source: A. Controlled Rotations: the θ-Angles · pp. 12–13
• For computational-basis-symmetric wavefunctions, symmetry constraints reduce the independent θ-angle degrees of freedom.Source: A. Controlled Rotations: the θ-Angles · pp. 12–13

B. Systematically Localizable Operators: the α-Angles

• The α-angles decompose each controlled-rotation level from local interactions outward, isolating sensitivity to long-distance controls.Source: B. Systematically Localizable… · pp. 13–14
• For free scalar ground states, the α-representation exposes exponentially suppressed rotation angles as control extent increases, enabling high-fidelity locality truncation.Source: B. Systematically Localizable… · p. 13
• The α indices identify the originating level, number of controls, and binary-interpreted control value, while recursive relations translate between α- and θ-angles.Source: B. Systematically Localizable… · pp. 13–14
• Partial decompositions can retain selected control counts, including site-wise controls, by setting selected α-angles to zero and solving subsequent angles recursively.Source: B. Systematically Localizable… · p. 14
• Symmetries reduce the number of independent angles, including to 2^(nQ−1) − 1 for symmetric wavefunctions.Source: B. Systematically Localizable… · pp. 14–15
• GHZ entanglement is captured by localized nearest-neighbor operations, whereas W-state preparation retains one increasingly nonlocal θ-angle per level because its two-qubit negativity remains nonzero across separations.Source: B. Systematically Localizable… · p. 15

Circuit localization

• α-angles scale with the K matrix and inverse two-point correlation function, rather than mutual information or negativity, with distance set by the maximally separated controlled qubits.The perturbative hierarchy follows spatial separation through powers of K_ij and asymptotically decays as e^(-m r)/r^(3/2).Source: §IV · p. 16

Numerical scaling

• Numerical free-field results show α-angle magnitudes bounded by an inverse-two-point-function form, while θ-angles lack the same spatial correlation pattern.For N = 10, nQ = 2, φmax = 3.5, and m̂ = 0.3, the extracted Mα is approximately 2m̂ at this lattice volume.Source: §IV · p. 16, Figure 3 · p. 17

Truncation and fidelity

• Fidelity converges exponentially as control distance is increased, with the N = 10 fit giving Z = 0.97(10) and η = 0.98(2).The result supports high-fidelity preparation even with limited spatial connectivity.Source: §IV · p. 18, Figure 6 · p. 19
• Truncating α-angles, controls, or both follows spatial and correlation hierarchies and permits systematic fidelity improvement by restoring omitted operations.This makes α-angle circuits preferable to θ-angle circuits for initializing scalar-field simulations.Source: §IV · pp. 18–19

Limitations

• The free-field evidence is limited by small classically simulated systems, finite-volume effects, lattice artifacts, and a discrepancy between fitted and lightest spectral masses.Larger classical calculations are required to establish the finite-volume interpretation firmly.Source: §IV · p. 19

V. INTERACTING λφ4 SCALAR FIELD THEORY

• The interacting calculation constructs and diagonalizes the digitized Hamiltonian classically, then extracts θ- and α-angles from the resulting ground-state wavefunction.Systems with up to nine spatial sites were treated using Mathematica and a Lanczos algorithm.Source: §V · p. 20
• For the interacting theory, Mα = 1.86(2) agrees with M2pt = 1.82(1) and MG = 1.86(1), supporting α-angle suppression by correlation functions.The scalar Hamiltonian mass is Mφ = 2.22, while finite-volume and digitization effects remain present.Source: §V · p. 21
• In interacting λφ4 theory, α-angles decrease exponentially with control distance, enabling exponentially precise ground-state preparation at a fixed control distance.The demonstration uses N = 8, nQ = 2, m = 1.6, φmax = 1.3, and λ = 32.Source: Figure 7 · p. 20, §V · p. 21

VI. SUMMARY AND CONCLUSIONS

• The demonstrated free and interacting systems have substantial finite-volume effects and significant two-qubit digitization errors, which are reduced with three or more qubits per site.The results are insufficient to inform continuum simulations with fully quantified uncertainties and require much larger classical calculations.Source: §VI · p. 22
• The approach is anticipated to apply beyond one-dimensional scalar theory to systems with a mass gap or localized spatial correlations, including QCD and the standard model.The paper specifically identifies efficient preparation of pion fields as a nuclear-physics application.Source: §VI · pp. 22–23

Appendix A: Continuous Lattice Scalar Field Theory with Four Spatial Sites

• The four-site periodic lattice uses allowed momenta to construct a position-space circulant coupling matrix K.The first row contains onsite, nearest-neighbor, and next-nearest-neighbor couplings.Source: Appendix A · p. 23

Appendix B: Mutual Information

• The reduced density matrix is obtained by partitioning target-site fields from traced degrees of freedom, exploiting Gaussian integrations.Its spectrum is organized through eigenvalues of the reduced matrix and oscillator eigenvalue towers.Source: Appendix B · pp. 23–24
• The reduced-density-matrix analysis does not assume particular lattice boundary conditions or correlation structure.The diagonal terms can be decoupled by a change of variables.Source: Appendix B · p. 23
• At large separation, inter-site reduced-density-matrix parameters scale as γ12, β12 ∼ e^−mj/jp, while the sites approach a tensor-product state.The off-diagonal couplings vanish asymptotically, leaving independently described site density matrices.Source: Appendix B · p. 24
• The mutual information scales as ∼e−2mj/jp′ at large separation, falling exponentially faster than the two-point correlations.Perturbative eigenvalue shifts contribute quadratically to the mutual information.Source: Appendix B · pp. 24–25

Appendix C: Negativity

• Negativity complements correlation functions and mutual information by probing entanglement through the partial transpose of a subsystem.For two sites, the first site's Hilbert space is partially transposed without tracing out degrees of freedom.Source: Appendix C · p. 25, §1 · p. 25

1. A Lattice with Two Sites

• The partial-transpose eigenvectors are symmetric and antisymmetric oscillator combinations, with negative contributions appearing in opposite-sign eigenvalue pairs.Positive eigenvalues do not contribute to the negativity.Source: §1 · p. 25
• For two coupled sites, the negativity is determined by the summed negative eigenvalues and depends on the inter-site coupling K12.The sign of K12 changes which symmetry sector receives positive or negative eigenvalues, but not the negativity itself.Source: §1 · pp. 25–26
• When K12 = 0, the negativity vanishes because the two sites decouple into independent fields.This provides the explicit zero-entanglement limit of the two-site calculation.Source: §1 · p. 26

2. Lattices with more than Two Sites

• For sites 0 and j in larger lattices, negativity follows the mutual-information calculation after transforming γ12 and β12 into partial-transpose matrices.The transformation exchanges γ12 and −β12 before diagonalizing the reduced problem.Source: §2 · p. 26
• When negativity is nonzero, one of the two decoupled oscillators contributes no negativity while the other contains alternating-sign eigenvalues.The total negativity can therefore be reduced to the negativity of the second oscillator.Source: §2 · p. 26

Appendix D: Digitized Lattice Scalar Field Theory with Two Spatial Sites

• The two-site digitized theory uses N = 2 spatial sites and nQ = 2 qubits per site, giving each scalar field a four-dimensional Hilbert space.The basis of the full field theory is the tensor product of the single-site digitized bases.Source: Appendix D · p. 26
• The digitized operator generates an entangled wavefunction with explicit φ →−φ and site-exchange symmetries.The normalization and amplitudes are determined by the coupling-matrix elements K11 and K12.Source: Appendix D · pp. 26–27
• The two-site example extends straightforwardly to more sites but does not capture the richer correlation and entanglement structure, especially its falloff with separation.This limits what can be inferred about spatial localization from the two-site calculation.Source: Appendix D · p. 27

Appendix E: θ- and α-Distributions at Fixed Control Distance h

• The α-angle hierarchy follows the coupling matrix K and decays exponentially with physical separation and scalar mass.The maximum α-angle magnitude decreases exponentially as the control distance increases.Source: Appendix E · p. 27
• At h = 16, 17, approximately half of the θ-angles at each level ℓ = h lie near θ ∼π 2.These fixed-distance distributions are representative of the behavior at other control distances.Source: Appendix E · p. 28, Figure 9 · p. 28
• The number of O(1) θ-rotation operators grows as 2h−1, so highly nonlocal operators dominate the count of entangling operations.Unlike θ-angles, α-angles decrease continuously within fixed h as well as across h.Source: Appendix E · p. 28

Recursive angle construction

• The resulting α-angles are generated recursively from the θ-angles using Eq. (33), with α55,k determined from θ5,k, α51,k, and α53,k modulo the relevant control-space sizes.The displayed relation is α55,k = θ5,k − α51,k mod 2 − α53,k mod 8.Source: Appendix F · pp. 28–29

Site-wise decomposition

• The example uses N = 3 lattice sites with nQ = 2 qubits per site and selects one largest-h operator at each physical distance r.This choice incorporates entire neighboring sites and motivates the site-wise α-angle protocol.Source: Appendix F · p. 29

Controlled-operator decompositions

• For the sixth qubit, a five-qubit-controlled θ5,k operator is decomposed into controlled operators with h = 1, 3, and 5, excluding h = 0, 2, and 4.The decomposition is described as physically motivated by the underlying quantum field theory.Source: Appendix F · p. 29
• For the fifth qubit, the chosen decomposition contains one uncontrolled single-qubit rotation and controlled operators with h = 2 and h = 4.These components arise from decomposing the four-control h = 4 operator defined by θ4,k angles.Source: Appendix F · p. 29

Spatial grouping and numerical example

• The angles are grouped by spatial separation as r = 0, 1, and 2, and Table I reports θ- and α-angle values for the free-field three-site example.The numerical setup uses open boundary conditions, m = 0.3, and φmax = 3.5.Source: Appendix F · p. 30, Table I · p. 31