⚛️ quantum physics

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

This paper proposes a first-quantised framework for simulating relativistic quantum systems on near-term devices by discretizing wavefunctions on a grid, utilizing finite-difference methods for momentum operators, and employing perturbative expansions to enable variational estimation of both non-relativistic and relativistic ground-state energies under periodic and Dirichlet boundary conditions.

Original authors: Jaewoo Joo, Timothy P. Spiller, Kyunghyun Baek, Jeongho Bang

Published 2026-06-30

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Jaewoo Joo, Timothy P. Spiller, Kyunghyun Baek, Jeongho Bang

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine you are trying to predict how a tiny, fast-moving particle behaves. In the everyday world, we use standard rules of physics. But when particles move incredibly fast—close to the speed of light—those standard rules break down, and we need "relativistic" rules to get the answer right.

This paper presents a new "recipe" for teaching a quantum computer how to solve these fast-particle puzzles, specifically for two different types of "playgrounds" (boundary conditions).

Here is the breakdown using simple analogies:

1. The Problem: The "Square Root" Nightmare

In physics, the energy of a fast-moving particle involves a complicated math formula that looks like a giant square root. It's very hard for a computer to calculate a square root directly, especially when you are trying to simulate a quantum system.

The Paper's Solution: Instead of trying to calculate the square root directly, the authors use a "trick" called perturbation theory. Think of it like approximating a curved hill by stacking flat blocks on top of each other.

The first block is the standard, slow-moving energy (what we already know).
The next blocks are small "correction" layers that account for the speed.
By adding just a few of these correction layers, the computer can get a very good estimate of the total energy without needing to solve the impossible square root equation.

2. The Setup: The Digital Grid

To simulate this on a quantum computer, the authors turn the continuous space where the particle moves into a digital grid.

Imagine a ruler divided into tiny, tiny segments.
The quantum computer uses "qubits" (quantum bits) to represent these segments. If you have more qubits, you get more segments, making the ruler more precise.
The particle's wave (its position and movement) is represented as a series of numbers sitting on these grid points.

3. The Two Playgrounds: Periodic vs. Dirichlet

The paper tests this method on two different types of "walls" that the particle might hit:

Periodic Boundary Conditions (PBC): The Pac-Man World
Imagine a video game like Pac-Man. If the character walks off the right edge of the screen, they instantly reappear on the left edge. The world is a loop.
- In the paper: The math is simpler here because the wave just flows smoothly from one end to the other. The computer uses a "quantum adder" (a tool that shifts the particle's position by one grid step) to calculate the energy.
Dirichlet Boundary Conditions (DBC): The Trapped Bird
Imagine a bird trapped in a cage with solid walls. If the bird hits the wall, it bounces back or stops; it doesn't wrap around. The wave must be zero at the walls.
- In the paper: This is trickier. The "wrap-around" tool doesn't work perfectly here. The authors had to invent a special "correction tool" (an extra mathematical operator) to fix the edges where the particle hits the wall. They showed that you can still use the results from the "Pac-Man" world and just add a few extra calculations to get the "Trapped Bird" answer.

4. How the Computer Does It: The "Shift and Measure" Game

The authors designed a specific set of instructions (a quantum circuit) for the computer to follow:

Prepare the State: Set up the particle's wave on the grid.
The "Shift" Move: Use a special gate to "shift" the particle one step to the right (or left) on the grid.
The Interference: By shifting the wave and then measuring how it overlaps with its original position, the computer can figure out how "curvy" the wave is.
The Result: In physics, how curvy the wave is tells you the momentum (how fast it's moving). Once you have the momentum, you can calculate the energy.

The paper shows that by combining these "shifts" with the "correction layers" (the perturbation theory mentioned earlier), the computer can estimate the energy of the particle in both the "Pac-Man" world and the "Trapped Bird" world.

5. Why This Matters (According to the Paper)

Near-Term Use: This method is designed for "near-term" quantum devices. These are the early, somewhat noisy quantum computers we have today, not the perfect, massive ones of the future.
Efficiency: It allows scientists to use variational optimization. Think of this as a "trial and error" game where the computer tweaks its settings to find the lowest possible energy state (the ground state) for the particle.
Versatility: It works for both the "looping" world and the "walled" world, giving researchers a flexible toolkit.

In Summary:
The authors have built a bridge between complex relativistic physics and the limited hardware of today's quantum computers. They turned a hard square-root problem into a series of simple "shift and measure" steps, allowing the computer to estimate the energy of fast-moving particles in different environments with reasonable accuracy.

Technical Summary: Relativistic Quantum Simulation under Periodic and Dirichlet Boundary Conditions

Problem Statement
Simulating quantum systems is computationally demanding due to the exponential growth of complexity with system size, a challenge exacerbated when relativistic effects must be incorporated to accurately describe physical phenomena. While perturbation theory is a standard approach for studying relativistic effects beyond the non-relativistic framework, implementing these simulations on near-term quantum devices remains difficult. Existing methods often rely on second quantization, which is standard for quantum chemistry but may lack intuitive handling for systems with non-linearities. There is a need for a practical, first-quantized framework capable of estimating relativistic ground-state energies for systems under different boundary conditions on near-term hardware.

Methodology
The authors propose a first-quantized framework for relativistic quantum simulation that discretizes the wavefunction across a finite grid represented by system qubits. The core methodology involves the following components:

Discretization and Translation Operators:
The continuous 1D spatial wavefunction is mapped to a discrete grid of $2^L$ points using $L$ qubits, where the position $x$ is scaled to the domain $0 \le x < 1$ . The squared momentum operator $\hat{p}^2$ is approximated using the finite-difference method based on quantum translation operators. The translation operator $\hat{A}_x$ (quantum adder) shifts the position state $|j\rangle$ to $|j+\delta x\rangle$ , while its conjugate $\hat{A}_x^\dagger$ (quantum subtractor) shifts in the opposite direction. The squared momentum is expressed as $\hat{p}^2 = -(mc L_m)^2 (\hat{A}^{(1)} - 2\hat{1})$ , where $L_m = \lambda_m / \delta x$ is a dimensionless parameter involving the reduced Compton wavelength.
Perturbative Expansion of Kinetic Energy:
The relativistic kinetic energy Hamiltonian $\hat{H}_{re}^K$ is derived from the square-root form of the relativistic energy-momentum relation. It is expanded into a perturbative series of higher-order momentum terms ( $\hat{p}^2, \hat{p}^4, \hat{p}^6, \dots$ ) valid in the regime where $\langle \hat{p}^2 \rangle \ll m^2 c^2$ . The total relativistic kinetic energy is approximated by truncating this series, incorporating coefficients derived from the Klein-Gordon equation.
Boundary Condition Handling:
The framework addresses two distinct boundary conditions:
- Periodic Boundary Conditions (PBC): The translation operators preserve unitarity by wrapping the grid ( $\hat{A}_x |2^L-1\rangle = |0\rangle$ ). The expectation values of even powers of momentum are expressed as combinations of $\langle \hat{A}^{(l)} \rangle$ .
- Dirichlet Boundary Conditions (DBC): Here, the wavefunction does not wrap ( $\psi(0) \neq \psi(1)$ ). The authors introduce a correction operator $\hat{E}_0 = |2^L-1\rangle\langle 0| + |0\rangle\langle 2^L-1|$ to account for the boundary. The expectation values for DBC are formulated as the PBC terms plus additional terms involving $\hat{E}_0$ and related operators ( $\hat{E}_1, \hat{E}_2$ ), which depend on the probability amplitudes at the grid boundaries.
Quantum Simulation Algorithm:
The authors design quantum circuits to estimate the required expectation values:
- Kinetic Energy: Controlled-translation gates are applied iteratively to a control qubit. Measuring the control qubit in the Z-basis yields the real part of $\langle (\hat{A}_x^\dagger)^l \rangle$ , which corresponds to the necessary momentum terms.
- Potential Energy: A density matrix $\hat{\rho}_V$ representing the potential is used in a circuit to estimate $\langle \hat{V} \rangle$ via the probabilities of the system state.
- Variational Optimization: The total energy is estimated by minimizing the sum of kinetic and potential energy expectation values over a set of variational parameters $\{\bar{\kappa}\}$ in an ansatz state $|\psi(\{\bar{\kappa}\}, x)\rangle$ .

Key Contributions

First-Quantized Recipe: The paper provides a specific recipe for relativistic simulation using first quantization, offering an alternative to the widely used second-quantization methods, particularly beneficial for systems with non-linearities.
Unified Boundary Framework: It presents a unified approach to handle both Periodic and Dirichlet boundary conditions. Crucially, it demonstrates that the data obtained for PBC can be reused to estimate DBC energies by adding specific boundary correction terms, reducing computational overhead.
Near-Term Device Compatibility: The method utilizes standard quantum circuit primitives (controlled-adders, Hadamard gates, and measurements) and variational optimization, making it suitable for implementation on near-term quantum devices.
Explicit Circuit Designs: The authors provide detailed circuit designs for computing the expectation values of the translation operators and the boundary correction terms required for the perturbative expansion.

Results and Claims
The paper claims to offer a "practical route" for simulating relativistic effects on near-term quantum devices. The theoretical results show that:

Relativistic kinetic energy terms up to a specific order of perturbation can be decomposed into linear combinations of expectation values of translation operators ( $\langle \hat{A}^{(l)} \rangle$ ) and boundary operators ( $\langle \hat{E}_d \rangle$ ).
The coefficients for these expansions can be pre-calculated based on physical parameters ( $m, c, \hbar, \delta x$ ).
The method allows for the variational optimization of ansatz states to estimate both non-relativistic and relativistic ground-state energies.

Significance
The authors position this work as a foundational step toward supporting future developments in quantum physics and chemistry. By enabling the simulation of relativistic effects on near-term hardware, the framework addresses a gap in current quantum simulation capabilities. The paper emphasizes that while increasing the number of qubits improves grid resolution, care must be taken to ensure the dimensionless parameter $L_m$ remains small enough for the perturbative expansion to remain valid. The work does not claim to solve the full Dirac equation or handle ultra-relativistic limits but focuses specifically on the regime where relativistic corrections to the non-relativistic framework are relevant. Future work is suggested to explore higher-order Laplacian operators and more accurate ansatz construction methods.