Bateman Method: Algebraic & Numerical Insights

• The Bateman method is a set of mathematical constructions that combines algebraic techniques for quadratic Hamiltonians with interpolation-based separable expansions for multidimensional integral equations.
• It provides a framework for constructing ladder operators that yield explicit eigenvalues and spectra for non-Hermitian quantum systems, as seen in the Bateman Hamiltonian formulation.
• The approach also encompasses the Bateman–Hillion ansatz, which generates exact solutions for scalar wave equations and underpins high-precision numerical methods in quantum scattering.

The term "Bateman method" refers to a set of mathematical constructions originating from work by H. Bateman on dissipative systems and advanced through numerous extensions. It encompasses both algebraic techniques for quadratic Hamiltonians—most notably the Bateman Hamiltonian for damping—and interpolation-based separable expansions for multidimensional integral equations, as well as solution-generating techniques for wave equations via the Bateman–Hillion ansatz. These diverse applications reflect the Bateman method’s foundational role in non-Hermitian quantum mechanics, numerical quantum scattering, and exact solutions to the wave equation.

1. Algebraic Structure and the Bateman Hamiltonian

The canonical Bateman Hamiltonian is given, in dimensionless units ($\hbar=1$, $\omega=1$), by:

$$H = (p_x^2 - p_y^2) + (x^2 - y^2) - b(xp_y + y p_x),$$

with $[x, p_x] = [y, p_y] = i$, all other commutators vanishing, and $b = \gamma / (m \omega)$ the dimensionless damping parameter. In this formulation, $x$ and $y$ are rescaled coordinates, and $p_x$, $p_y$ are conjugate momenta.

The adjoint (regular) matrix representation is constructed by evaluating commutators $[H, O_i]$ for the operator basis $\omega=1$0. The resulting $\omega=1$1 adjoint matrix yields, via its characteristic equation, four eigenvalues—two pairs of positive and negative "natural frequencies":

$\omega=1$2

This matrix-algebraic method provides the explicit frequencies and corresponding ladder operators of the system (Fernández, 2020).

2. Ladder Operators, Spectra, and Commutation Relations

For each eigenvalue $\omega=1$3, explicit ladder operators $\omega=1$4 are constructed as nontrivial solutions of $\omega=1$5:

• $\omega=1$6, $\omega=1$7
• $\omega=1$8, $\omega=1$9
• $$H = (p_x^2 - p_y^2) + (x^2 - y^2) - b(xp_y + y p_x),$$0, $$H = (p_x^2 - p_y^2) + (x^2 - y^2) - b(xp_y + y p_x),$$1
• $$H = (p_x^2 - p_y^2) + (x^2 - y^2) - b(xp_y + y p_x),$$2, $$H = (p_x^2 - p_y^2) + (x^2 - y^2) - b(xp_y + y p_x),$$3

The commutation relations are:

• $$H = (p_x^2 - p_y^2) + (x^2 - y^2) - b(xp_y + y p_x),$$4
• $$H = (p_x^2 - p_y^2) + (x^2 - y^2) - b(xp_y + y p_x),$$5; all other commutators vanish.

Two towers of formal eigenfunctions can be constructed, each as a family generated from a vacuum annihilated by a pair of lowering operators. For example, the "shifted-oscillator" vacuum $$H = (p_x^2 - p_y^2) + (x^2 - y^2) - b(xp_y + y p_x),$$6 leads to eigenstates

$$H = (p_x^2 - p_y^2) + (x^2 - y^2) - b(xp_y + y p_x),$$7

with spectra

$$H = (p_x^2 - p_y^2) + (x^2 - y^2) - b(xp_y + y p_x),$$8

However, none of these formal eigenfunctions are $$H = (p_x^2 - p_y^2) + (x^2 - y^2) - b(xp_y + y p_x),$$9-normalizable since $[x, p_x] = [y, p_y] = i$0 is Hermitian but possesses complex eigenvalues (Fernández, 2020).

3. Bateman Interpolation Scheme for the Lippmann–Schwinger Equation

A distinct application of the Bateman method is in the construction of separable expansions for potentials in the two-variable Lippmann–Schwinger equation, bypassing partial-wave decomposition, especially advantageous in high-energy or multinucleon problems (Kuruoglu, 2013). The essential steps are:

• Construct a grid: Choose $[x, p_x] = [y, p_y] = i$1 nodes in the $[x, p_x] = [y, p_y] = i$2 (momentum modulus) variable and $[x, p_x] = [y, p_y] = i$3 nodes in $[x, p_x] = [y, p_y] = i$4, forming a Cartesian product grid $[x, p_x] = [y, p_y] = i$5.
• Build the interpolant: Form the Bateman rank-$[x, p_x] = [y, p_y] = i$6 interpolant for a potential $[x, p_x] = [y, p_y] = i$7 as

$[x, p_x] = [y, p_y] = i$8

where $[x, p_x] = [y, p_y] = i$9 is the inverse of the sampled potential matrix.

• The separable form enables transforming the Lippmann–Schwinger equation into a finite-dimensional linear system for the $b = \gamma / (m \omega)$0-matrix elements. This structure supports decoupling the interpolation (matrix) grid from the quadrature (integration) grid, thus flexibly balancing computational effort and accuracy.

This Bateman interpolation approach achieves multi-digit precision for the nucleon-nucleon $b = \gamma / (m \omega)$1-matrix with only a few hundred interpolation points, as direct convergence studies in the literature confirm (Kuruoglu, 2013).

4. Exact Beam Solutions by the Bateman–Hillion Ansatz

A further facet of the Bateman method involves constructing exact solutions of the scalar wave equation via the Bateman–Hillion ansatz (Ducharme, 2014). The general construction is:

$b = \gamma / (m \omega)$2

where $b = \gamma / (m \omega)$3 solves the standard paraxial equation in variable $b = \gamma / (m \omega)$4. To ensure the physical beam exhibits intensity fall-off as $b = \gamma / (m \omega)$5, a space-time constraint $b = \gamma / (m \omega)$6 is enforced, typically in integrals via a Dirac delta function. This endows the solution with spherical wavefronts at large distances.

In the paraxial regime, this formalism recovers the equivalence between the optical paraxial wave and non-relativistic Schrödinger equations under the identification $b = \gamma / (m \omega)$7, highlighting a deep correspondence in the functional forms of Gouy phase evolution in optics and matter waves (Ducharme, 2014).

5. Numerical Implementation, Convergence, and Limitations in the Bateman Interpolation Approach

For the Bateman interpolation scheme, grid construction and numerical integration details are essential for accuracy and stability. Key features include:

• Partitioning $b = \gamma / (m \omega)$8 (momentum axis) and $b = \gamma / (m \omega)$9 ($x$0-axis) into finite elements; using Gauss–Legendre quadrature points for numerical integration.
• Typical convergence rates show 3–4 digit agreement with Nystrom benchmarks for $x$1 grid points. Stability to $x$2 is achieved for $x$3 (Kuruoglu, 2013).
• When the interpolation and quadrature grids coincide, the Bateman scheme reduces to the Nystrom method. The Bateman method's flexibility lies in decoupling grid choice for interpolation and accuracy control.
• The approach generalizes to three- and four-body problem kernels, but matrix sizes grow exponentially with variable number, making large-scale application challenging unless iterative techniques or domain decomposition are used.

6. Physical Interpretation, Generalizations, and Context

The Bateman method unifies a number of algebraic and numerical approaches to non-Hermitian systems, separable expansions in quantum scattering, and exact solution families for wave equations. The Bateman–Hillion solution framework, in particular, yields exact Hermite–Gaussian beams with non-paraxial validity, structured via space–time constraints to mirror spherical wavefront properties.

Historically, Bateman’s original work on dissipative systems led to several quantization attempts and algebraic treatments of non-Hermitian Hamiltonians. Subsequent works by Hillion and others (notably Deguchi–Fujiwara and Fernández) have clarified the spectral, algebraic, and functional analytic structure of the Bateman Hamiltonian and its ladder-operator framework (Fernández, 2020), while applications in numerical quantum scattering demonstrate the method's practical impact (Kuruoglu, 2013). The Bateman method thus persists as a key technique in mathematical physics, with applications ranging from foundational studies of non-Hermitian quantum systems to efficient computational algorithms for high-dimensional integral equations.

---

References:

• Algebraic treatment of the Bateman Hamiltonian (Fernández, 2020)
• Bateman method for two-body scattering without partial-wave decomposition (Kuruoglu, 2013)
• Constrained Bateman-Hillion Solutions for Hermite-Gaussian Beams (Ducharme, 2014)