Deeply Bound Pionic-Atom States

• Deeply bound pionic-atom states are π⁻-nucleus systems defined by low-lying orbitals with MeV-scale binding energies and narrow widths, serving as precision probes of nuclear matter.
• Their formation through recoil-less reactions like (d,³He) enables high-resolution spectroscopy that clearly separates 1s, 2p, and even 2s states.
• Analysis of these states refines in-medium πN interaction parameters, elucidates neutron density distributions, and supports evidence for partial chiral symmetry restoration.

Deeply bound pionic-atom states are atomic states of a $\pi^-$ bound to a nucleus in which the pion occupies low-lying orbitals such as $1s$, $2p$, or, in some theoretical studies, $2s$, with binding energies of several MeV and widths of order a few hundred keV. They are described by a Klein–Gordon equation in the Coulomb field of the nucleus plus a complex pion–nucleus optical potential, and they are formed not through the ordinary atomic X-ray cascade but through recoil-less nuclear reactions such as $(d,{}^3{\rm He})$. Their quantitative description has made them a precision laboratory for the in-medium $\pi N$ interaction, neutron-density distributions, neutron skins, the renormalization of the isovector amplitude $b_1$, and partial restoration of chiral symmetry in nuclear matter (Gal, 4 Aug 2025).

1. Bound-state definition and optical-potential framework

A pionic atom is described by the Klein–Gordon equation for a $\pi^-$ of reduced mass $\mu$ moving in the Coulomb field $V_c(r)$ of the nucleus plus a complex optical strong-interaction potential $1s$0:

$1s$1

with $1s$2. At leading order in the nuclear density, the $1s$3-wave part of $1s$4 is parametrized as

$1s$5

where $1s$6, $1s$7, $1s$8 is the nucleon mass, and $1s$9 is a phenomenological absorption parameter $2p$0 (Gal, 4 Aug 2025).

The negative isovector scattering length $2p$1 is repulsive. Because of this repulsion, low-lying states such as $2p$2 are pushed partly out of the nucleus, which keeps the absorptive imaginary part of $2p$3 small. The result is that deeply bound states can remain narrow even in very heavy nuclei. The same framework is often written in the Ericson–Ericson form, with explicit $2p$4-wave terms and additional parameters $2p$5, $2p$6, $2p$7, and $2p$8, in analyses that solve the bound-state problem and then embed it in formation-reaction calculations (Ikeno et al., 2022).

Free-space $2p$9 scattering lengths used in this context are quoted as

$2s$0

in one summary of global pionic-atom fits, while a sigma-term analysis uses

$2s$1

within its chosen implementation (Gal, 4 Aug 2025, Ikeno et al., 2022). The medium dependence of the isovector term is commonly expressed through Weise’s formula,

$2s$2

which directly connects the optical potential to finite-density chiral dynamics (Gal, 4 Aug 2025).

2. Formation mechanism and spectroscopic method

Deeply bound $2s$3 and $2s$4 states cannot be reached by atomic X-ray cascade because higher levels such as $2s$5 and $2s$6 are so strongly absorbed that radiative yields vanish. Their observation therefore relies on recoil-less “stripping” reactions at forward angles, which implant a $2s$7 almost at rest in the residual nucleus. The canonical example is

$2s$8

with the kinematics tuned to a momentum transfer $2s$9 (Gal, 4 Aug 2025).

The standard experimental configuration uses deuterons of $(d,{}^3{\rm He})$0 incident on a thin lead or tin target of a few $(d,{}^3{\rm He})$1. The outgoing $(d,{}^3{\rm He})$2 is momentum-analyzed in a high-resolution magnetic spectrometer at $(d,{}^3{\rm He})$3, with typical $(d,{}^3{\rm He})$4, corresponding to $(d,{}^3{\rm He})$5–$(d,{}^3{\rm He})$6 (FWHM). The missing-mass spectrum,

$(d,{}^3{\rm He})$7

shows peaks at the pion binding energies. Variants of this method have been applied at GSI and, more recently, at RIKEN for Sn isotopes at $(d,{}^3{\rm He})$8, with comparable or improved resolution (Gal, 4 Aug 2025).

The reaction is commonly treated in DWIA or effective-number formalisms. In finite-angle studies, different combinations of pion-bound and neutron-hole states dominate at different scattering angles because of the matching condition. Near the recoil-free limit, $(d,{}^3{\rm He})$9 formation is strongest; at finite angles, $\pi N$0 and other configurations become more prominent. Representative calculations for $\pi N$1 at $\pi N$2 show that the dominant structures shift from $\pi N$3 at $\pi N$4 to $\pi N$5-dominated configurations at $\pi N$6 and $\pi N$7 (Ikeno et al., 2011).

High-resolution theoretical studies further indicate that non-yrast states such as $\pi N$8 are expected to be visible in $\pi N$9 spectra. Their simultaneous observation with the ground $b_1$0 state is helpful both for reducing uncertainties associated with neutron wave functions and for reducing experimental uncertainties associated with calibration of the absolute excitation energy (Ikeno et al., 2011).

3. Measured binding energies, widths, and spectral systematics

Representative experimental values for deeply bound states in heavy nuclei are summarized below. These quantities are extracted by fitting missing-mass peaks to Lorentzian line shapes folded with the experimental resolution function (Gal, 4 Aug 2025).

System	$b_1$1 (MeV)	$b_1$2 (MeV)
$b_1$3 $b_1$4	$b_1$5	$b_1$6
$b_1$7 $b_1$8	$b_1$9	$\pi^-$0
$\pi^-$1 $\pi^-$2	$\pi^-$3	$\pi^-$4
$\pi^-$5 $\pi^-$6	$\pi^-$7	$\pi^-$8

Within the Klein–Gordon plus optical-potential framework, theoretical predictions reproduce both $\pi^-$9 and $\mu$0 within $\mu$1 (Gal, 4 Aug 2025). This level of agreement is central to the interpretation of deeply bound states as controlled probes of the $\mu$2-wave $\mu$3 interaction in nuclei.

A separate high-statistics study of the $\mu$4 spectrum observed the atomic $\mu$5 and $\mu$6 states as distinct peak structures in the missing-mass spectrum, and reported that the $\mu$7 state in a Sn nucleus was observed for the first time. The same experiment measured the spectrum at finite reaction angles for the first time and determined the formation cross sections between $\mu$8 and $\mu$9 (Nishi et al., 2017).

Systematics across nuclei reinforce several recurring features. Heavy nuclei support deeply bound $V_c(r)$0 levels with MeV-scale binding energies, yet the widths remain modest because the repulsive $V_c(r)$1 term reduces the pion’s overlap with the absorptive nuclear interior. This is the technical basis for the phrase “deeply bound” in this field: the states are deeply bound in energy, but not strongly broadened by absorption (Gal, 4 Aug 2025).

4. In-medium $V_c(r)$2 dynamics, $V_c(r)$3, and chiral symmetry restoration

The most consequential interpretation of deeply bound pionic-atom spectroscopy is its sensitivity to the density dependence of the isovector $V_c(r)$4-wave amplitude $V_c(r)$5. In the low-density chiral description, the medium dependence is written as

$V_c(r)$6

and the isoscalar channel is supplemented by a density-dependent double-scattering correction,

$V_c(r)$7

in one explicit implementation (Ikeno et al., 2022).

Global fits of “normal” and “deeply bound” pionic-atom X-ray level shifts and widths across the periodic table $V_c(r)$8–$V_c(r)$9 extract

$1s$00

while the isoscalar amplitude is less well constrained but is consistent with

$1s$01

These fits are presented as clear experimental evidence for reduction of the pion decay constant and partial restoration of the QCD chiral condensate $1s$02 in the nuclear medium (Gal, 4 Aug 2025).

The sensitivity of specific observables to $1s$03 has been quantified explicitly. For tin isotopes, numerical solutions show that the $1s$04 binding energy $1s$05 and width $1s$06 vary almost linearly with $1s$07 over $1s$08–$1s$09. For $1s$10, the sensitivities are

$1s$11

per $1s$12 change in $1s$13. For the lighter and more sensitive $1s$14, they are

$1s$15

per $1s$16 change in $1s$17 (Ikeno et al., 2022).

Using experimental precisions of roughly $1s$18, $1s$19, and $1s$20–$1s$21, the corresponding $1s$22 uncertainties are estimated as $1s$23 from $1s$24 alone in $1s$25, $1s$26 from $1s$27 alone, and $1s$28–$1s$29 from the $1s$30–$1s$31 energy gap in $1s$32. This identifies the gap and the $1s$33 width in lighter Sn isotopes as the best observables for a few-MeV determination of $1s$34, provided the neutron density is known accurately (Ikeno et al., 2022).

Large-scale analyses that fit 100 pionic-atom data points, including deeply bound states, reaffirm the density-dependent renormalization of $1s$35. In one such fit, keeping $1s$36 density-independent gives

$1s$37

which is $1s$38–$1s$39 more repulsive than the free value $1s$40. Imposing the chiral ansatz for $1s$41 restores full agreement with the free $1s$42 value and reduces $1s$43 by $1s$44–$1s$45 units. The same analysis infers

$1s$46

corresponding to a $1s$47–$1s$48 drop of $1s$49 in nuclear matter (Friedman et al., 2014).

5. Neutron densities, neutron skins, and effective density

Because the isovector term $1s$50 shifts levels in proportion to the local neutron–proton density difference, deeply bound pionic levels in isotope chains serve as probes of neutron density distributions and neutron skins. In Sn, Pb, and Zr isotopes, precise measurements of $1s$51 and $1s$52 level shifts have been presented as sensitive to neutron skins of

$1s$53

(Gal, 4 Aug 2025).

In a 100-point global analysis, neutron densities were parameterized by

$1s$54

with best fits yielding $1s$55 and implying for $1s$56

$1s$57

This result was reported to be in line with parity-violation and proton-scattering results (Friedman et al., 2014).

A complementary perspective is provided by the effective density sampled by the atomic pion. Defining the overlapping density

$1s$58

and the peak position $1s$59 by $1s$60, one obtains $1s$61. Numerical calculations across $1s$62–$1s$63 and for states $1s$64, $1s$65, $1s$66, and others give

$1s$67

This near-constancy means that atomic pions probe nearly the same $1s$68 over a wide range of nuclei and orbitals (Ikeno et al., 2011).

That property has a dual significance. It stabilizes the interpretation of extracted in-medium amplitudes, since different states are not sampling widely different densities. At the same time, it limits direct differential access to the density dependence away from $1s$69. A related study concluded that extracting $1s$70 at densities different from $1s$71 requires $1s$72 accuracy together with simultaneous $1s$73–$1s$74 data in order to break Seki–Masutani correlations (Ikeno et al., 2011).

6. Reaction-theory discrepancies and refined experimental strategies

Although the level energies and widths are well described, the formation reaction still contains unresolved issues. In the angular-distribution measurement of $1s$75, the observed reaction-angle dependence of each state was reproduced in shape by theoretical calculations, but the absolute magnitude of the pionic $1s$76 formation cross section showed a significant discrepancy. At $1s$77,

$1s$78

corresponding to

$1s$79

For the $1s$80 state at the same angle,

$1s$81

giving

$1s$82

Possible origins listed for the $1s$83 suppression include higher-order reaction mechanisms, inadequacies in the DWIA distorted waves or spectroscopic factors, sensitivity to the short-range $1s$84-wave part of the optical potential, and many-body in-medium modifications of the $1s$85 vertex beyond the standard Ericson–Ericson term (Nishi et al., 2017).

Parameter correlations also limit what can be extracted from spectroscopy. The same observables that constrain $1s$86 correlate strongly with the two-nucleon absorption parameter $1s$87. For $1s$88, contours of constant $1s$89 are nearly parallel in the $1s$90 plane, making separation difficult; neutron-density uncertainties can shift the best-fit $1s$91 by several MeV; and the angular ratio $1s$92 in calculated $1s$93 formation spectra is almost flat in $1s$94 and cannot be used alone to determine it (Ikeno et al., 2022).

One refined strategy is to use odd-neutron targets. In $1s$95, the picked-up neutron is the extra $1s$96 neutron, so the final pionic level is a pure $1s$97 configuration and does not suffer additional shifts from residual neutron-hole–pion couplings. Numerical predictions for this case give, for $1s$98,

$1s$99

$2p$00

before experimental resolution. With a $2p$01 FWHM instrumental resolution, the peaks remain well separated, and for a beam intensity of $2p$02 and a $2p$03 Sn target, the expected yield is $2p$04 events/day in the $2p$05 region within $2p$06, sufficient to determine $2p$07 to $2p$08 and $2p$09 to $2p$10 (Ikeno et al., 2013).

Taken together, these developments define the present status of deeply bound pionic-atom states. The basic bound-state phenomenology is robust; the spectroscopy has reached the level at which $2p$11, $2p$12, and potentially $2p$13 states can be separated with sub-MeV resolution; the in-medium renormalization of $2p$14 and the associated partial restoration of chiral symmetry are supported by global fits and isotope studies; and the remaining challenges lie primarily in reaction theory, correlated optical-potential parameters, and precision control of neutron-density inputs (Gal, 4 Aug 2025).