⚛️ phenomenology

Bound-state decay of a charged scalar particle in an external Coulomb field

This paper calculates the decay rate of a negatively charged kaon bound to a light nucleus decaying into pions by employing a subtraction method to overcome slow partial-wave convergence and deriving a small- $\alpha Z$ expansion that agrees with all-orders numerical results.

Original authors: Andrzej Czarnecki, Artem O. Davydov

Published 2026-09-30

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Andrzej Czarnecki, Artem O. Davydov

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, particles are rarely solitary travelers. They often find themselves trapped in the electric grip of an atomic nucleus, orbiting like tiny planets around a sun. When these trapped particles decay, or break apart, the process is not the same as it would be for a free particle drifting through empty space. The invisible electric field of the nucleus tugs on the fragments as they fly away, altering the speed and likelihood of the breakup. Physicists have long known that this "Coulomb field" changes the rules of decay, but calculating exactly how much it changes them has been a stubborn mathematical puzzle, especially for atoms with a small number of protons. Understanding these subtle shifts is crucial because they test the limits of our theories about how matter behaves under extreme conditions, bridging the gap between the laws governing single particles and the complex reality of bound systems.

A team of researchers at the University of Alberta has now solved a key piece of this puzzle by creating a simplified model to study how a heavy, negatively charged particle decays while bound to a light nucleus. They focused on a scenario where a particle, which they treated as a type of kaon, sits in the lowest energy orbit around a nucleus and then splits into two new particles: a charged pion and a neutral pion. While real kaons are subject to other forces that complicate their behavior, the researchers stripped the problem down to its electrical core. They treated the nucleus as a fixed, heavy point and the particles as simple, point-like spheres governed by the laws of quantum mechanics. Their goal was to determine the precise rate at which this bound particle decays compared to an identical particle floating freely in space.

The challenge they faced was one of mathematical precision. When the researchers tried to calculate the decay rate using standard methods, the numbers refused to settle. The calculation involves adding up an infinite series of contributions from different angles of motion, but for light nuclei, this series converges so slowly that the computer has to work through thousands of terms just to get a rough answer. It is like trying to hear a whisper in a storm; the signal they are looking for—the tiny difference caused by the nucleus's electric field—is drowned out by the massive background noise of the decay itself. Previous attempts to solve similar problems for muons had to skip over the lightest nuclei entirely because the math became too unstable to trust.

To overcome this, the researchers invented a new way of looking at the problem. Instead of trying to calculate the entire decay process from scratch every time, they realized they could calculate the "free" decay—the part that would happen even without the nucleus—and subtract it out. They treated the interaction with the nucleus as a small correction to a known baseline. By isolating this tiny correction, they could ignore the massive, noisy background and focus only on the subtle distortion caused by the electric field. This subtraction method acted like a filter, allowing them to see the small signal clearly without needing to compute millions of unnecessary terms.

Using this technique, the team calculated the decay rates for nuclei with atomic numbers ranging from one to ten, covering the lightest elements in the periodic table. They found that the bound particle decays slightly slower than a free one, a suppression that grows as the nuclear charge increases. More importantly, they confirmed that their new subtraction method was far more efficient, reducing the computational effort required by about thirty percent while delivering results that matched their theoretical predictions with high precision. They also developed a mathematical expansion to check their work, showing that their numerical results agreed perfectly with the theoretical expectations for light nuclei.

The study revealed that the electric field of the nucleus does not just nudge the decay rate; it introduces a specific, predictable pattern of change. The researchers found that the difference between the bound and free decay rates follows a clear mathematical rule, dominated by a term related to the square of the nuclear charge, with a smaller, more complex correction appearing at the fourth power. This finding aligns with previous theories for muon decay but arises from a different physical mechanism in their simplified model, highlighting how the absence of particle spin changes the underlying dynamics. By successfully navigating the difficult low-charge region where other methods fail, the team has provided a robust tool for future studies. Their work demonstrates that by separating the known from the unknown, scientists can extract precise answers from problems that previously seemed too messy to solve, paving the way for more accurate calculations of how particles behave when trapped in the heart of an atom.

Technical Summary: Bound-State Decay of a Charged Scalar Particle in an External Coulomb Field

Problem Statement
The paper addresses the calculation of the decay rate $\Gamma$ for a particle bound in the Coulomb field of a nucleus, specifically examining the interplay between decay dynamics, binding effects, and final-state interactions. While previous studies have analyzed bound-muon decay, numerical formulations relying on partial-wave expansions of continuum wave functions face significant challenges for light nuclei ( $Z < 6$ ). In this regime, the binding parameter $\alpha Z$ is small, yet the angular-momentum series converges slowly. Direct partial-wave calculations are computationally expensive because each term involves highly oscillatory radial integrals, and reconstructing a decay rate ratio close to unity with sufficient accuracy to resolve small Coulomb corrections is numerically demanding.

Methodology
To investigate these convergence issues, the authors employ a simplified scalar model: a negatively charged scalar particle ( $K^-$ ) bound to a static, infinitely heavy, point-like nucleus decays into charged ( $\pi^-$ ) and neutral ( $\pi^0$ ) scalar particles via a local three-scalar interaction. All particles are treated as point-like, and the final pions are massless. The system is described using the Klein-Gordon equation in an external Coulomb potential.

The core methodological innovation is a subtraction method designed to overcome the slow convergence of the partial-wave expansion:

Separation of Contributions: The total decay rate is decomposed into a "free pion" contribution (neglecting final-state Coulomb interactions) and a "Coulomb-distortion" correction ( $\Delta C$ ).
Analytic Reduction: The free pion contribution is evaluated separately using a one-dimensional integral representation derived from the Fourier transform of the bound-state wave function. This avoids the need for a partial-wave sum for the dominant term.
Partial-Wave Application: The partial-wave expansion is applied only to the difference between the interacting pion wave function and the free pion wave function (the Coulomb distortion). Since this correction is significantly smaller than the total rate, the partial-wave sum converges much faster and requires a lower cutoff ( $l_{max}$ ).
Analytic Verification: The numerical results are cross-verified using a diagrammatic small- $\alpha Z$ expansion through order $(\alpha Z)^4$ . This involves calculating the imaginary part of the kaon self-energy in the external field, considering zero, one, and two Coulomb photon insertions on the pion line, including seagull and ladder diagrams.

Key Results

Numerical Efficiency: For light nuclei ( $1 \le Z \le 10$ ), the subtraction method reduces the required partial-wave cutoff by approximately 25–30% compared to direct calculations. For $Z=1$ , the cutoff is reduced from $l_{max}=1136$ (direct) to $802$ (subtracted). The results for $1 \le Z \le 5$ from both methods agree to four significant digits.
Decay Rate Ratio: The ratio of bound-to-free decay rates is determined numerically. The paper presents the rate suppression $1 - \Gamma/\Gamma_0$ for $Z$ up to 10.
Analytic Expansion: The diagrammatic expansion yields the following series for the decay rate ratio:
$\frac{\Gamma}{\Gamma_0} = 1 - \frac{1}{2}(\alpha Z)^2 + \frac{17}{24}(\alpha Z)^4 + O((\alpha Z)^5)$
The leading term $-(\alpha Z)^2/2$ arises entirely from the free-pion contribution. The Coulomb distortion of the outgoing pion does not contribute at order $(\alpha Z)^2$ but provides the leading correction at order $(\alpha Z)^4$ .
Coulomb Distortion Coefficient: The coefficient of the $(\alpha Z)^4$ term arising from the Coulomb distortion is found to be $4/3$ . This is confirmed numerically, as the quantity $\Delta C / (\alpha Z)^4$ approaches $4/3$ as $Z$ decreases.
Cancellation Mechanism: At order $(\alpha Z)^4$ , the two-photon contributions (seagull graph and two linear insertions) cancel exactly ( $+2(\alpha Z)^4 - 2(\alpha Z)^4 = 0$ ), leaving the one-photon contribution as the sole source of the distortion correction at this order.

Significance and Claims
The paper claims that the subtraction method effectively resolves the numerical instability associated with low- $Z$ bound-state decay calculations, a problem that previously led researchers to omit light nuclei ( $Z < 6$ ) from bound-muon decay analyses. By isolating the small Coulomb distortion correction, the method avoids the numerical difficulty of reconstructing a ratio near unity from large, nearly cancelling terms.

The authors emphasize that while the leading coefficient of the suppression ( $-1/2$ ) coincides with the result for bound-muon decay (the Überall result), the physical origin differs: in this scalar model, the Coulomb distortion of the outgoing particle only enters at order $(\alpha Z)^4$ , whereas in the spin-1/2 muon case, it contributes at order $(\alpha Z)^2$ . The work provides a validated framework and tools (the subtraction strategy and diagrammatic expansion) for future, more realistic calculations of bound-state decays, including the full spin-1/2 problem of bound-muon decay.

Technical Summary: Bound-State Decay of a Charged Scalar Particle in an External Coulomb Field

More like this