---
title: Weakly Relativistic Lattice Hamiltonians
url: https://www.emergentmind.com/topics/weakly-relativistic-lattice-hamiltonians
type: topic
---

# Weakly Relativistic Lattice Hamiltonians

Weakly relativistic lattice Hamiltonians are first-quantized finite-grid Hamiltonians obtained by starting from the positive-energy relativistic kinetic operator and truncating its low-momentum expansion at order $\hat p^4$. On a one-dimensional finite domain, this construction requires boundary-consistent lattice realizations of both the second and fourth momentum moments, so the Hamiltonian depends explicitly on the boundary condition. In the formulation presented in "First-Quantized Relativistic Quantum Simulation with Periodic and Dirichlet Boundary Conditions" [2606.21794], the relevant cases are periodic boundary conditions (PBC) and Dirichlet boundary conditions (DBC), with the corresponding weakly relativistic kinetic operator written as
\[
\hat T_\tau=\frac{\hat P_\tau^2}{2m}-\frac{\hat P_\tau^4}{8m^3c^2},\qquad \hat P_\tau^4\equiv (\hat P_\tau^2)^2,\qquad \tau\in\{\mathrm P,\mathrm D\}.
\]
The central structural point is that weakly relativistic simulation on a lattice requires reconstruction of both $\langle \hat P_\tau^2\rangle$ and $\langle \hat P_\tau^4\rangle$, and that these operators are boundary-dependent rather than universal.

## 1. Relativistic origin and lattice definition

The continuum starting point is the positive-energy single-particle kinetic operator
\[
\hat T_{\rm rel}=mc^2\left(\sqrt{1+\frac{\hat p^2}{m^2c^2}}-1\right).
\]
In the weakly relativistic regime, the expansion is truncated at fourth order,
\[
\hat T_{\rm rel}=\frac{\hat p^2}{2m}-\frac{\hat p^4}{8m^3c^2}+O\!\left(\frac{\hat p^6}{m^5c^4}\right),
\]
so the lattice problem is reduced to constructing boundary-consistent versions of $\hat p^2$ and $\hat p^4$ on a finite grid [2606.21794].

On the lattice, the total Hamiltonian is
\[
\hat H_\tau^{\rm tot}=\hat T_\tau+\hat V_\tau,
\]
with diagonal potential operator
\[
\hat V_\tau=\sum_{j=0}^{N-1}V\!\left(x_j^{(\tau)}\right)|j\rangle\langle j|.
\]
The paper emphasizes that the boundary condition is part of the Hamiltonian definition, not an afterthought: PBC and DBC correspond to different lattice momentum operators. This is decisive for the relativistic correction, because the leading correction depends on $\langle \hat P_\tau^4\rangle$ in addition to $\langle \hat P_\tau^2\rangle$.

## 2. Finite-domain discretization and first-quantized encoding

The wavefunction is encoded in a first-quantized $L$-qubit register,
\[
\psi=\sum_{j=0}^{N-1}c_j\,|j\rangle,\qquad N=2^L.
\]
The finite domain is discretized differently for the two boundary conditions [2606.21794].

For PBC, the grid is
\[
x_j^{(\mathrm P)}=j\Delta_{\mathrm P},\qquad \Delta_{\mathrm P}=\frac{R}{N},
\]
so the lattice is cyclic and the last point is connected back to the first.

For DBC, the grid is
\[
x_j^{(\mathrm D)}=(j+1)\Delta_{\mathrm D},\qquad \Delta_{\mathrm D}=\frac{R}{N+1},
\]
with virtual endpoints
\[
\psi(0)=0,\qquad \psi(R)=0.
\]
In this case the computational basis stores only the interior points of an open chain; the physical boundary points are not represented as basis states.

This difference in grid definition is not merely notational. Under PBC the natural discrete kinetic structure is generated by cyclic translation. Under DBC the desired operator is the open-chain finite-difference Laplacian, and the finite register must therefore exclude the unphysical wrap-around link that would otherwise couple the two ends of the chain.

## 3. Periodic boundary conditions and cyclic-translation reconstruction

For PBC, the basic operator is the unitary cyclic translation
\[
\hat A|j\rangle=|j+1\ \mathrm{mod}\ N\rangle,\qquad \hat A^\dagger|j\rangle=|j-1\ \mathrm{mod}\ N\rangle,
\]
with
\[
\hat A^\dagger\hat A=\hat A\hat A^\dagger=\hat{\mathds 1},\qquad \hat A^N=\hat{\mathds 1}.
\]
The standard second-order finite-difference momentum-squared operator is
\[
\hat P_{\rm P}^2=-\frac{\hbar^2}{\Delta_{\rm P}^2}\left(\hat A+\hat A^\dagger-2\hat{\mathds 1}\right),
\]
and its expectation value reduces to a single real translation moment,
\[
\langle \hat P_{\rm P}^2\rangle=\frac{2\hbar^2}{\Delta_{\rm P}^2}\left(1-\mathrm{Re}\langle \hat A\rangle\right).
\]

The fourth moment is obtained by squaring the same operator,
\[
\hat P_{\rm P}^4=\frac{\hbar^4}{\Delta_{\rm P}^4}\left(\hat A^2+\hat A^{\dagger 2}-4\hat A-4\hat A^\dagger+6\hat{\mathds 1}\right),
\]
which implies
\[
\langle \hat P_{\rm P}^4\rangle=\frac{\hbar^4}{\Delta_{\rm P}^4}\left(2\,\mathrm{Re}\langle \hat A^2\rangle-8\,\mathrm{Re}\langle \hat A\rangle+6\right).
\]
Defining translation moments
\[
m_l=\mathrm{Re}\langle \hat A^l\rangle,\qquad l=1,2,
\]
the PBC kinetic sector is completely determined by $m_1$ and $m_2$. The corresponding weakly relativistic kinetic expectation value is
\[
\langle \hat T_{\rm P}\rangle=\frac{1}{2m}\langle \hat P_{\rm P}^2\rangle-\frac{1}{8m^3c^2}\langle \hat P_{\rm P}^4\rangle.
\]

Operationally, this means that under PBC the relativistic correction can be reconstructed entirely from cyclic-translation measurements. That feature gives the PBC case a particularly compact estimator structure: no boundary-specific observables are required.

## 4. Dirichlet boundary conditions and boundary-local correction terms

For DBC, the desired kinetic operator is the open-chain finite-difference Laplacian rather than the cyclic one. The construction in [2606.21794] expresses the open-chain operator using the same cyclic translation $\hat A$, supplemented by boundary-local terms that remove the unphysical wrap-around link.

The boundary coupling operator is defined as
\[
\hat E_0=|N\!-\!1\rangle\langle 0|+|0\rangle\langle N\!-\!1|.
\]
This isolates the coupling between the first and last interior points. The DBC momentum-squared operator is then
\[
\hat P_{\rm D}^2=-\frac{\hbar^2}{\Delta_{\rm D}^2}\left(\hat A+\hat A^\dagger-2\hat{\mathds 1}-\hat E_0\right).
\]
If one introduces the cyclic finite-difference operator on the DBC grid spacing,
\[
\hat P_{\rm cyc}^2(\Delta_{\rm D})=-\frac{\hbar^2}{\Delta_{\rm D}^2}\left(\hat A+\hat A^\dagger-2\hat{\mathds 1}\right),
\]
then
\[
\hat P_{\rm D}^2=\hat P_{\rm cyc}^2(\Delta_{\rm D})+\frac{\hbar^2}{\Delta_{\rm D}^2}\hat E_0.
\]
Its expectation value is therefore
\[
\langle \hat P_{\rm D}^2\rangle=\langle \hat P_{\rm cyc}^2(\Delta_{\rm D})\rangle+\frac{\hbar^2}{\Delta_{\rm D}^2}\langle \hat E_0\rangle,
\]
with
\[
\langle \hat E_0\rangle=c_{N-1}^*c_0+c_0^*c_{N-1}.
\]
The correction is explicitly a boundary coherence.

The fourth moment is more intricate because $\hat P_{\rm D}^4=(\hat P_{\rm D}^2)^2$ generates additional boundary terms. The derived expression is
\[
\hat P_{\rm D}^4=\hat P_{\rm cyc}^4(\Delta_{\rm D})+\frac{\hbar^4}{\Delta_{\rm D}^4}\left(4\hat E_0-\hat E_1-\hat E_2+\hat E_0^2\right),
\]
where
\[
\hat E_1=\hat A\hat E_0+\hat E_0\hat A^\dagger,\qquad \hat E_2=\hat E_0\hat A+\hat A^\dagger\hat E_0.
\]
The explicit endpoint forms are
\[
\hat E_0^2=|0\rangle\langle 0|+|N-1\rangle\langle N-1|,
\]
\[
\hat E_1=2|0\rangle\langle 0|+|1\rangle\langle N-1|+|N-1\rangle\langle 1|,
\]
\[
\hat E_2=2|N-1\rangle\langle N-1|+|0\rangle\langle N-2|+|N-2\rangle\langle 0|.
\]
Hence
\[
\langle \hat P_{\rm D}^4\rangle=\langle \hat P_{\rm cyc}^4(\Delta_{\rm D})\rangle+\frac{\hbar^4}{\Delta_{\rm D}^4}\left(4\langle \hat E_0\rangle-\langle \hat E_1\rangle-\langle \hat E_2\rangle+\langle \hat E_0^2\rangle\right).
\]

In measurable form, the paper writes
\[
\langle \hat P_{\rm D}^2\rangle=\frac{\hbar^2}{\Delta_{\rm D}^2}\left(2-2m_1+b_0\right),
\]
\[
\langle \hat P_{\rm D}^4\rangle=\frac{\hbar^4}{\Delta_{\rm D}^4}\left(2m_2-8m_1+6+4b_0-b_1-b_2+b_{00}\right),
\]
with
\[
b_0=\langle \hat E_0\rangle,\qquad b_{00}=\langle \hat E_0^2\rangle,\qquad b_1=\langle \hat E_1\rangle,\qquad b_2=\langle \hat E_2\rangle.
\]
The DBC kinetic energy is therefore
\[
\langle \hat T_{\rm D}\rangle=\frac{1}{2m}\langle \hat P_{\rm D}^2\rangle-\frac{1}{8m^3c^2}\langle \hat P_{\rm D}^4\rangle.
\]

The physical interpretation of the correction terms is localized: $\hat E_0$ is endpoint-to-endpoint coherence, $\hat E_0^2$ gives endpoint probabilities, and $\hat E_1,\hat E_2$ involve endpoints and near-endpoints. Their role is to remove the unphysical periodic wrap-around link and restore the open-chain DBC Laplacian exactly.

## 5. Measurement primitives and energy-estimation workflow

The energy-estimation protocol separates kinetic and potential contributions [2606.21794]. For PBC, kinetic estimation requires only the translation moments
\[
m_1=\mathrm{Re}\langle \hat A\rangle,\qquad m_2=\mathrm{Re}\langle \hat A^2\rangle.
\]
These are measured using a Hadamard-test-style controlled-unitary circuit: ancilla in $|0\rangle$, apply Hadamard, apply controlled-$\hat A^l$, apply Hadamard, and measure ancilla in $Z$. The circuit yields
\[
\langle Z_C\rangle=\mathrm{Re}\langle \psi|\hat U|\psi\rangle,
\]
and for $\hat U=\hat A^l$ gives $m_l$.

For DBC, the same $m_1,m_2$ are needed, together with the boundary terms $b_0,b_1,b_2,b_{00}$. The boundary coherences are obtained from overlap probabilities via
\[
\langle |f\rangle\langle g|+|g\rangle\langle f|\rangle = 2P_{fg}^+-P_f-P_g,
\]
where
\[
P_f=|\langle f|\psi\rangle|^2,\qquad P_g=|\langle g|\psi\rangle|^2,\qquad P_{fg}^+=|\langle s_{fg}^+|\psi\rangle|^2,
\]
and
\[
|s_{fg}^+\rangle=\frac{|f\rangle+|g\rangle}{\sqrt 2}.
\]
Consequently, DBC requires only a small number of endpoint and near-endpoint overlap measurements in addition to the translation measurements.

For the potential term, since $\hat V_\tau$ is diagonal in the computational basis,
\[
\langle \hat V_\tau\rangle=\sum_{j=0}^{N-1}V(x_j^{(\tau)})|c_j|^2,
\]
it is estimated by sampling the position register in the computational basis and classically post-processing the sampled positions. The total-energy estimate is
\[
\langle \hat H_\tau^{\rm tot}\rangle=\langle \hat T_\tau\rangle+\langle \hat V_\tau\rangle.
\]

This workflow makes the estimator structure boundary-aware. Under PBC the kinetic sector is reconstructed from translation measurements alone, whereas under DBC the same cyclic estimator is supplemented by boundary-local terms.

## 6. Benchmarks, validation, and error structure

The framework is validated in four benchmark classes [2606.21794]. For the free particle with PBC, the analysis uses Fourier modes $q_n$ and the exact lattice dispersion
\[
p_{\rm lat}^2(n)=\frac{4\hbar^2}{\Delta_{\rm P}^2}\sin^2\!\left(\frac{\pi n}{N}\right).
\]
The benchmark demonstrates convergence of the lattice relativistic energy to the continuum result as the grid is refined and separates discretization error, defined as the difference between lattice and continuum, from weak-relativistic truncation error, defined as the difference between the exact square-root relativistic energy and the $\hat P^4$ expansion.

For the infinite square well with DBC, the states are sine modes
\[
\psi_s(j)=\sqrt{\frac{2}{N+1}\sin\!\left[\frac{\pi s(j+1)}{N+1}\right].
\]
The open-chain spectrum is confirmed as
\[
p_{\rm D,lat}^2(s)=\frac{4\hbar^2}{\Delta_{\rm D}^2}\sin^2\!\left(\frac{\pi s}{2(N+1)}\right).
\]
Most importantly, the boundary-corrected cyclic estimator exactly reproduces the open-chain matrix moments $\langle \hat P_{\rm D}^2\rangle$ and $\langle \hat P_{\rm D}^4\rangle$. This is the direct validation of the boundary-local correction scheme.

The smooth-potential benchmarks use the periodic cosine potential
\[
V_{\rm P}(x)=V_0\left[1-\cos\left(\frac{2\pi x}{R}\right)\right]
\]
for PBC and the harmonic confining potential
\[
V_{\rm D}(x)=\frac12 m\omega^2\left(x-\frac{R}{2}\right)^2
\]
for DBC. The procedure diagonalizes the nonrelativistic Hamiltonian and then evaluates the relativistic correction on the ground state. The comparison between direct matrix evaluation and estimator reconstruction shows near-perfect agreement, validating the full workflow, including kinetic terms and diagonal potential sampling.

A finite-shot measurement test simulates the statistical behavior of the measurement primitives. The kinetic-energy estimator RMSE scales as
\[
O(M^{-1/2}),
\]
as expected for sampling error. PBC uses only translation moments, while DBC has a slightly larger constant prefactor because it includes the extra boundary overlap measurements. The reported agreement between estimator reconstruction and direct matrix evaluation, together with the explicit separation of finite-grid discretization, weak-relativistic truncation, and measurement errors, delineates the scope of the method: it is a boundary-consistent first-quantized framework for weakly relativistic quantum simulation on a one-dimensional finite grid.

Source: https://www.emergentmind.com/topics/weakly-relativistic-lattice-hamiltonians