Relativistic Feshbach-Villars Equation for Two Spin-$0$ Particles

Abstract: The Feshbach-Villars version of the relativistic quantum mechanics can be extended for two-body systems in such a way that the center-of-mass motion is separated off. The procedure results in an equation of Feshbach-Villars-type in terms of the relative coordinate.

• The paper extends the first-order Feshbach–Villars framework to two interacting spin-0 particles and derives exact center-of-mass separation, producing a relative-motion equation with reduced mass in the kinetic term and total mass in the rest-energy term.
• The authors solve the resulting equation using a Lippmann–Schwinger formulation, Coulomb–Sturmian basis expansion, and matrix continued fractions that treat long-range dynamics accurately while approximating only short-range interactions.
• Applications to pionic hydrogen and pionium produce relativistic binding energies that differ from reduced-mass Klein–Gordon results by approximately 10^-3–10^-2 atomic units for ground states, while highlighting unresolved issues such as retardation and cluster separability.

The Feshbach–Villars (FV) representation of the Klein–Gordon equation has long been recognized as the natural Hamiltonian formulation of relativistic spin-0 quantum mechanics, but it has remained essentially a one-particle formalism. In "Relativistic Feshbach-Villars Equation for Two Spin-0 Particles" (2605.02161), Z. Papp extends the FV framework to two-body systems and demonstrates that the center-of-mass kinetic energy separates exactly, yielding an effective one-body FV equation in the relative coordinate. This is a structurally significant result: center-of-mass separability is precisely the property that most ad hoc relativistic few-body constructions fail to achieve without imposing additional constraints, and its emergence here "naturally" is presented as the key step toward a consistent relativistic few-particle theory.

Motivation: why the Klein-Gordon equation cannot serve as a foundation

The paper opens by restating two well-known deficiencies of the Klein–Gordon (KG0) equation. First, it is second order in time, so specifying $\Psi$ alone does not determine the state; the stationary form contains $E^2$, which blocks any additive composition of single-particle energies into a many-body energy. Second, the KG0 equation admits no consistent extension to interacting few-particle systems for this reason. The FV reformulation of Feshbach and Villars resolves the first-order-in-time problem by splitting the KG wave function into particle ($\phi$) and antiparticle ($\chi$) components, governed by a Schrödinger-like equation

$$i\hbar \frac{\partial}{\partial t}|\psi\rangle = H_{FV0}|\psi\rangle,$$

with

$$H_{FV0} = (\tau_3 + i\tau_2)\frac{p^2}{2m} + \tau_3 mc^2 + (\tau_3+i\tau_2)U + I_2 V,$$

where $U = S + S^2/(2mc^2)$ collects scalar-interaction contributions and $V$ is the time-like component of a Lorentz four-vector potential. The Hamiltonian is Hermitian only in the generalized sense $H_{FV0} = \tau_3 H_{FV0}^\dagger \tau_3$, with real eigenvalues and indefinite normalization $\langle\psi|\tau_3|\psi\rangle = \pm 1$ distinguishing particles from antiparticles.

A central structural feature—and the principal computational obstacle—is that the components are coupled by the kinetic energy operator itself through $E^2$0. Unlike ordinary coupled-channel problems, where inter-channel couplings are short-range potentials vanishing asymptotically, this coupling persists at all distances. Moreover, since $E^2$1, the coupling cannot be eliminated by matrix inversion. Any solution method must therefore confront the kinetic coupling head-on.

The two-body construction

For two free spin-0 particles with masses $E^2$2 and $E^2$3, the additivity of energy in the FV eigenvalue equation allows the unified Hamiltonian to be written directly as a sum of single-particle FV Hamiltonians plus a relative-coordinate interaction:

$E^2$4

with $E^2$5. Because the interaction is assumed to depend only on the relative coordinate $E^2$6—a plausible assumption given that each particle's interaction can originate only from the other—the standard non-relativistic kinetic-energy decomposition applies verbatim inside the $E^2$7 prefactor:

$E^2$8

with reduced mass $E^2$9. Dropping the center-of-mass term leaves the effective relative-motion Hamiltonian

$\phi$0

Two features of this result deserve emphasis. The kinetic term carries the reduced mass $\phi$1, while the rest-energy separation between particle-like and antiparticle-like sectors is set by $\phi$2, not $\phi$3. The author argues explicitly that a $\phi$4 separation would violate the non-relativistic limit, since the rest energy of two particles must be $\phi$5; this is presented as a physical consistency check on the construction rather than a fitted choice.

Solution method: Lippmann-Schwinger equation with matrix continued fractions

The numerical machinery follows the author's earlier work on single-particle FV equations [(2605.02161) references therein]. The Hamiltonian is split into a long-range part $\phi$6 (kinetic energy, rest mass, Coulomb or confining long-range pieces) and a short-range part $\phi$7, and the bound-state problem is cast as a Lippmann–Schwinger equation

$\phi$8

The short-range operator is expanded on a Coulomb–Sturmian basis truncated at $\phi$9 radial functions per partial wave; because the basis is bi-orthonormal, this yields a finite homogeneous algebraic problem whose determinant condition supplies the eigenvalues. Crucially, only the short-range piece is approximated—the long-range Green's operator is treated exactly via the continued-fraction representation of Kónya, Lévai, and Papp for infinite symmetric tridiagonal (Jacobi-matrix) Hamiltonians. In the present setting each element of the tridiagonal matrix is itself a $\chi$0 FV block, so the scalar continued fraction becomes a matrix continued fraction. The author notes that this scheme converges rapidly for bound states and can be analytically continued to scattering energies.

Numerical results: pionic hydrogen and pionium

The method is applied to the $\chi$1–$\chi$2 (pionic hydrogen) and $\chi$3–$\chi$4 (pionium) Coulomb systems using physical masses ($\chi$5, $\chi$6). Representative binding energies in atomic units are:

System	State	KG0 (reduced mass)	FV0 (this work)
$\chi$7–$\chi$8	$\chi$9	$$i\hbar \frac{\partial}{\partial t}|\psi\rangle = H_{FV0}|\psi\rangle,$$0	$$i\hbar \frac{\partial}{\partial t}|\psi\rangle = H_{FV0}|\psi\rangle,$$1
$$i\hbar \frac{\partial}{\partial t}|\psi\rangle = H_{FV0}|\psi\rangle,$$2–$$i\hbar \frac{\partial}{\partial t}|\psi\rangle = H_{FV0}|\psi\rangle,$$3	$$i\hbar \frac{\partial}{\partial t}|\psi\rangle = H_{FV0}|\psi\rangle,$$4	$$i\hbar \frac{\partial}{\partial t}|\psi\rangle = H_{FV0}|\psi\rangle,$$5	$$i\hbar \frac{\partial}{\partial t}|\psi\rangle = H_{FV0}|\psi\rangle,$$6
$$i\hbar \frac{\partial}{\partial t}|\psi\rangle = H_{FV0}|\psi\rangle,$$7–$$i\hbar \frac{\partial}{\partial t}|\psi\rangle = H_{FV0}|\psi\rangle,$$8	$$i\hbar \frac{\partial}{\partial t}|\psi\rangle = H_{FV0}|\psi\rangle,$$9	$$H_{FV0} = (\tau_3 + i\tau_2)\frac{p^2}{2m} + \tau_3 mc^2 + (\tau_3+i\tau_2)U + I_2 V,$$0	$$H_{FV0} = (\tau_3 + i\tau_2)\frac{p^2}{2m} + \tau_3 mc^2 + (\tau_3+i\tau_2)U + I_2 V,$$1
$$H_{FV0} = (\tau_3 + i\tau_2)\frac{p^2}{2m} + \tau_3 mc^2 + (\tau_3+i\tau_2)U + I_2 V,$$2–$$H_{FV0} = (\tau_3 + i\tau_2)\frac{p^2}{2m} + \tau_3 mc^2 + (\tau_3+i\tau_2)U + I_2 V,$$3	$$H_{FV0} = (\tau_3 + i\tau_2)\frac{p^2}{2m} + \tau_3 mc^2 + (\tau_3+i\tau_2)U + I_2 V,$$4	$$H_{FV0} = (\tau_3 + i\tau_2)\frac{p^2}{2m} + \tau_3 mc^2 + (\tau_3+i\tau_2)U + I_2 V,$$5	$$H_{FV0} = (\tau_3 + i\tau_2)\frac{p^2}{2m} + \tau_3 mc^2 + (\tau_3+i\tau_2)U + I_2 V,$$6

The KG0 values shown are computed with the reduced mass $$H_{FV0} = (\tau_3 + i\tau_2)\frac{p^2}{2m} + \tau_3 mc^2 + (\tau_3+i\tau_2)U + I_2 V,$$7 in the effective FV form, which the author explicitly labels as lacking physical basis; they are included only for comparison. The FV0 results differ from them at the level of $$H_{FV0} = (\tau_3 + i\tau_2)\frac{p^2}{2m} + \tau_3 mc^2 + (\tau_3+i\tau_2)U + I_2 V,$$8–$$H_{FV0} = (\tau_3 + i\tau_2)\frac{p^2}{2m} + \tau_3 mc^2 + (\tau_3+i\tau_2)U + I_2 V,$$9 atomic units in the ground states—a small but resolvable relativistic correction attributable to the correct $U = S + S^2/(2mc^2)$0 rest-energy scale. The robustness of the matrix continued fraction is what makes such fine distinctions numerically meaningful.

Limitations and open questions

Several caveats are stated plainly in the paper. The interaction model assumes instantaneous, relative-coordinate-dependent potentials; retardation effects—by analogy with the Liénard–Wiechert problem—are acknowledged as a major challenge for relativistic few-body theory, though the author suggests they may be less relevant for stationary eigenvalue problems. The numerical demonstration covers only pure Coulomb binding without strong-interaction corrections (e.g., pion-nucleon interactions in pionic hydrogen), and the comparison against KG0 uses a reduced-mass variant the author himself deems unphysical, so the tables serve more as internal consistency checks than as benchmarks against experiment. The claimed generalization to spinning particles and to genuine few-body (three-plus particle) systems remains a conjecture, not a derivation. Cluster separability beyond two bodies—ensuring isolated subsystems behave independently—is identified as the most crucial outstanding problem, of which the present center-of-mass separation is a necessary but not sufficient ingredient.

Conclusion

This work establishes that the Feshbach–Villars formalism, being a genuine first-order energy eigenvalue equation, permits exact separation of center-of-mass motion for two interacting spin-0 particles, producing an effective FV equation in the relative coordinate with reduced mass in the kinetic term and total mass in the rest-energy term. Combined with the matrix-continued-fraction solution technique, this yields precise spectra for pionic hydrogen and pionium and provides a concrete, constraint-free route toward relativistic few-body wave mechanics—an extension whose viability for spin and for three or more particles the paper identifies but does not yet demonstrate.