⚛️ nuclear theory

Direct Boundary Matching: A Bound-State Technique for Nuclear Scattering with Lagrange-Legendre Functions

This paper introduces a Direct Boundary Matching Method (DBMM) that solves nuclear scattering problems, including coupled-channel cases, by incorporating outgoing wave boundary conditions directly into a Lagrange-Legendre $L^2$ matrix formulation, thereby eliminating the need for Bloch operators or complex coordinate scaling while achieving accuracy comparable to traditional Numerov integration.

Original authors: Jin Lei

Published 2026-08-27

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Jin Lei

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

In the subatomic world, atoms are not static islands but dynamic systems where particles constantly collide, scatter, and interact. To understand how these collisions happen, physicists must solve a difficult mathematical puzzle: describing how a particle moves when it is far away from the target it is hitting. When a particle is trapped inside an atom, its behavior is relatively simple to calculate because it stays in one place, fading away quickly as it moves outward. However, when a particle is scattered, it travels forever, oscillating like a wave that never stops. This infinite, wavy nature makes it incredibly hard to use standard computer tools, which are designed for things that stay put. For decades, scientists have had to use complex, multi-step workarounds to force these endless waves into a format their computers could handle, often requiring them to split the problem into separate regions or rotate their mathematical coordinates into strange, imaginary spaces just to make the numbers work.

A researcher at Tongji University has now proposed a more direct way to solve this problem, bypassing the usual detours. By using a specific set of mathematical building blocks called Lagrange-Legendre functions, the new method treats the scattering problem as a single, unified equation. Instead of splitting the calculation or using complex tricks to manage the infinite waves, the researcher simply writes the rule for how the wave should behave at the very edge of the calculation zone directly into the final line of the computer's equation. This approach, called Direct Boundary Matching, allows the computer to solve for the scattering result in one go, keeping the entire calculation in real, physical space rather than twisting it into complex mathematical territory.

The method was tested by simulating the scattering of protons off carbon-12 nuclei at an energy of 30 million electron volts. The results matched almost perfectly with the established, highly trusted numerical methods that have been used for years to solve these problems. The new technique reproduced the probability of the particles bouncing off and the specific shifts in their wave patterns with extreme precision, showing differences so small they are barely measurable. The study confirms that this direct approach is not just a theoretical idea but a practical tool that works as well as the traditional, more complicated methods.

One of the most significant advantages of this new technique is how naturally it handles situations where multiple paths are possible. In many nuclear collisions, a particle might bounce off, or it might excite the target nucleus into a different state, creating a web of interconnected possibilities. Traditional methods often struggle to keep these different paths organized without adding layers of extra complexity. The new method, however, builds these connections directly into the structure of the equation, treating the different paths as parts of a single, larger system. This makes it particularly useful for creating fast, simplified models, or "emulators," that can predict the outcomes of complex collisions without needing to run the full, heavy calculation every time.

While the new method does not necessarily run faster than the old ones for a single, simple calculation, its strength lies in its clarity and flexibility. It removes the need for special operators and separate matching steps that have been standard for decades, replacing them with a straightforward set of rules that anyone can see and understand. The researcher has already made the code for this method available to the public, allowing other scientists to use it immediately. By providing a clear, direct path to solving these scattering problems, this work offers a simpler way to explore the fundamental interactions that shape the atomic nucleus, potentially speeding up future discoveries in nuclear physics without sacrificing accuracy.

Technical Summary: Direct Boundary Matching with Lagrange-Legendre Functions

Problem Statement
The theoretical description of quantum scattering presents a fundamental challenge distinct from bound-state calculations. While bound-state wave functions decay exponentially and can be accurately represented within finite domains using square-integrable ( $L^2$ ) bases, scattering states exhibit oscillatory asymptotic behavior extending to infinity. This oscillatory nature precludes straightforward $L^2$ representation, as truncated bases of decaying functions cannot capture infinite-range oscillations. Historically, this has necessitated specialized methods for continuum problems, creating a dichotomy between nuclear structure and nuclear reaction theory. Traditional approaches, such as numerical integration followed by asymptotic matching, become increasingly cumbersome for coupled-channel problems. While "bound-state techniques for the continuum" (e.g., R-matrix, Complex Scaling, Lorentz Integral Transform) have emerged to reformulate scattering within localized bases, they often rely on complex formal machinery, such as Bloch operators or coordinate rotations into the complex plane.

Methodology: Direct Boundary Matching Method (DBMM)
This paper presents the Direct Boundary Matching Method (DBMM), a technique that reformulates scattering problems into a localized, $L^2$ representation using Lagrange-Legendre basis functions. The method operates entirely in real coordinate space and is characterized by two primary features:

Basis Selection: The method employs Lagrange-Legendre functions defined on a finite interval $[0, R]$ . Unlike Lagrange-Laguerre functions, which possess an inherent exponential decay factor suitable only for bound states or complex scaling, Lagrange-Legendre functions are polynomial in nature and impose no asymptotic behavior. This neutrality allows the basis to represent both decaying and oscillatory functions, with the asymptotic character determined entirely by the boundary conditions.
Direct Boundary Incorporation: The core innovation of DBMM is the direct imposition of the outgoing wave boundary condition into the matrix equation, bypassing the need for Bloch operators or two-step matching procedures typical of R-matrix methods.
- The radial Schrödinger equation is discretized at interior mesh points using the Lagrange-Legendre basis and Baye's exact analytical derivative matrices.
- The outgoing wave condition, $\psi'_{sc}(R) = \gamma_s \psi_{sc}(R)$ , where $\gamma_s$ is the complex logarithmic derivative of the outgoing Coulomb-Hankel function, is directly encoded in the last row of the matrix equation.
- This results in a complex, non-symmetric $N \times N$ linear system where the first $N-1$ rows represent the Schrödinger equation at interior points, and the final row enforces the boundary condition.

Coupled-Channel Extension
The formalism is extended to coupled-channel problems without requiring additional formal machinery like generalized Bloch operators. The wave function is decomposed into an incident wave (present only in the entrance channel) and a scattered wave. This decomposition naturally leads to an effective source potential that distinguishes between the entrance channel (driven by the short-range part of the diagonal potential) and other channels (driven by the full coupling potential, including Coulomb multipoles). The resulting system is a block matrix equation of dimension $n_c N \times n_c N$ , where $n_c$ is the number of channels. The diagonal blocks follow the single-channel structure, while off-diagonal blocks contain the coupling potentials.

Results
Benchmark calculations were performed for proton scattering on $^{12}\text{C}$ at a laboratory energy of $E_{lab} = 30$ MeV, utilizing the Koning-Delaroche global optical potential. The results were compared against the established Numerov integration method:

Accuracy: The S-matrix elements computed by DBMM showed excellent agreement with the Numerov method across partial waves $\ell = 0$ to $10$. The differences in the modulus $|S_\ell|$ were below $2.5 \times 10^{-5}$ , and phase differences were less than $0.01^\circ$ .
Wave Functions: The radial wave functions (real and imaginary parts) computed via DBMM matched the Numerov solutions precisely throughout the computational domain, from the nuclear interior to the asymptotic region.
Convergence: A convergence study indicated that achieving high phase-shift accuracy ( $\Delta \delta < 0.01^\circ$ ) requires a specific number of mesh points $N$ depending on the boundary radius $R$ . The author notes that DBMM may require approximately 1.5–2 times more mesh points than variational R-matrix formulations to reach the same accuracy, leading to a potential runtime increase of a factor of 3–8 due to the $O(N^3)$ scaling of the matrix inversion.

Significance and Claims
The paper claims that the primary contribution of DBMM is methodological simplification and conceptual transparency rather than superior numerical efficiency.

Conceptual Transparency: Each row of the matrix equation corresponds to a clear physical statement (Schrödinger equation at interior points, boundary condition at the surface), making the method accessible for teaching and implementation.
Implementation Simplicity: The method requires only standard linear algebra routines, eliminating the need for the Bloch operator formalism or complex coordinate rotations.
Practical Utility:
- Emulator Construction: Unlike R-matrix methods which yield boundary values requiring reconstruction, DBMM provides the scattered wave function coefficients directly as part of the solution. This is critical for reduced-basis emulator applications, enabling the construction of unified reduced bases across coupled channels to achieve massive speedups ( $10^3$ – $10^4$ ) in continuum-discretized coupled-channels (CDCC) calculations.
- Few-Body Scattering: DBMM offers an alternative real-coordinate boundary-matching scheme that can interface with complex-scaled interior solutions, potentially extending few-body scattering calculations to lower energies where complex scaling faces numerical difficulties.

The work concludes that while DBMM does not outperform the R-matrix method in raw computational speed for individual calculations, its direct formulation and natural extension to coupled channels make it a valuable tool for nuclear reaction theory, particularly in the context of emulator development and pedagogical clarity. The method is implemented in the open-source Julia package SLAM.jl.

More like this