Interrelation between $ar{p}$-Ca Atom Spectra and Nuclear Density Profiles

This work studies $\bar{p}$-Ca atom spectra in light of the strong shifts and level widths, using the optical model with several types of parametric coefficients. The spectroscopic quantities are obtained as the eigenvalues of the Dirac equation, where the nuclear densities computed via nuclear Density Functional Theory and the effect of the anomalous magnetic moment are incorporated. The results indicate that the isovector term’s contribution to the optical potential is crucial for explaining the systematical differences in the strong shifts between $^{40}$Ca and $^{48}$Ca. Furthermore, it is found that both the strong shifts and the level widths exhibit significant dependence on the nuclear density profiles. These findings provide critical insights into the nuclear structures, particularly in the context of Calcium isotopes, by offering a more comprehensive understanding of the underlying nuclear-hadron properties.

💡 Research Summary

This paper investigates the strong interaction–induced energy shifts and level widths observed in antiprotonic calcium atoms (¯p‑Ca), focusing on the isotopes ^40Ca and ^48Ca. The authors employ a relativistic framework based on the Dirac equation, which naturally incorporates the large anomalous magnetic moment of the antiproton and the spin‑orbit coupling. The electromagnetic interaction is treated with a Coulomb potential derived from the nuclear charge density, supplemented by the leading‑order vacuum polarization correction.

The strong interaction is introduced via an optical potential V_opt = U – iW, expressed in the linear‑density approximation as a sum of isoscalar (b₀) and isovector (b₁) terms multiplied by the total (ρ₀ = ρ_n + ρ_p) and difference (ρ₁ = ρ_n – ρ_p) nucleon densities, respectively, plus gradient (p‑wave) terms (c₀, c₁). Three parameter sets are defined: Type I (b₀ only), Type II (b₀ + b₁), and Type III (b₀ + imaginary c₀). All sets share the same b₀ value taken from recent global fits (b₀ = 1.3 + 1.9 i fm). The isovector coefficient is reduced relative to scattering‑length estimates (b₁ = –8.0 + 1.7 i fm) to account for density folding, while the imaginary p‑wave strength is set to c₀ = 0 + 1.2 i fm³.

Nuclear density profiles are supplied in two ways. First, a three‑parameter Fermi (3pF) model reproduces the known proton charge radius and assumes simple neutron‑skin thicknesses (e.g., –0.16 fm for ^40Ca, –0.03 fm for ^48Ca). Second, self‑consistent Skyrme‑DFT calculations (SLy4, SLy5, SkM*) provide more realistic proton and neutron distributions. The DFT results predict a larger neutron‑skin for ^48Ca (≈0.14–0.15 fm) than the 3pF parametrizations, reflecting recent experimental indications of a “kink” in charge radii between the two isotopes.

The Dirac radial equations are discretized on a radial grid, assembled into a matrix, and diagonalized to obtain bound‑state energies. The strong shift ΔE is defined as the difference between the binding energy with and without the optical potential; the width Γ is proportional to the imaginary part of the optical potential evaluated at the bound‑state wave function. The calculations assume spherical symmetry for both isotopes, which is justified by their doubly‑magic nature.

Results show that the isoscalar‑only model (Type I) predicts almost identical shifts for ^40Ca and ^48Ca (≈3 eV), far below the experimentally observed difference of several eV. Inclusion of the isovector term (Type II) generates a systematic shift difference of about 5 eV, bringing the theoretical prediction much closer to the measured values (≈10 eV). Adding an imaginary p‑wave term (Type III) markedly increases the calculated widths, reproducing the experimentally observed broadening of tens of eV.

A key finding is the sensitivity of both ΔE and Γ to the underlying nuclear density. When DFT densities replace the simple 3pF profiles, the calculated shifts and widths change by 10–20 %, with larger neutron skins leading to more positive shifts and broader widths. This demonstrates that antiprotonic atom spectroscopy can serve as a high‑precision probe of subtle nuclear‑structure features such as neutron‑skin thickness and surface diffuseness.

The authors emphasize that their parameter sets are not intended to replace global optical‑potential fits but rather to explore nucleus‑specific effects in calcium isotopes. They argue that forthcoming high‑resolution X‑ray measurements, possibly employing transition‑edge sensor microcalorimeters, will be capable of resolving the spin‑orbit doublet of the (n,l) = (5,4) level. Such data would allow a stringent test of the isovector and p‑wave contributions and could refine the extraction of nuclear density parameters from antiprotonic atom spectra.

In conclusion, the study establishes that (i) the isovector b₁ term and the p‑wave c₀ term are essential for reproducing the systematic differences in strong shifts and widths between ^40Ca and ^48Ca, and (ii) realistic nuclear density profiles, especially those reflecting neutron‑skin variations, have a pronounced impact on the spectroscopic observables. These insights reinforce the role of antiprotonic atoms as a valuable tool for nuclear‑structure investigations and suggest a pathway toward more precise determinations of neutron‑skin thicknesses in medium‑mass nuclei.