Neutron skin thickness and its volume and surface contributions

TL;DR

This study systematically investigates neutron skin thickness in transuranium berkelium isotopes using advanced nuclear theory, revealing how shell effects, deformation, and anisotropy influence neutron distribution and contribute to understanding nuclear symmetry energy.

Contribution

It introduces a detailed decomposition of neutron skin into volume and surface contributions in deformed nuclei, highlighting the effects of shell closures and deformation on neutron skin properties.

Findings

Neutron skin thickness increases with neutron number, with shell effects causing anti-kinks.

Surface contribution dominates near the proton drip line, while volume contribution is generally larger.

Deformation enhances surface diffuseness and neutron skin, with anisotropic effects depending on nuclear shape.

Abstract

Accurate determination of the neutron skin thickness ($\Delta R_{\mathrm{np}}$) in finite nuclei is crucial for constraining the density dependence of the nuclear symmetry energy. In this work, we systematically investigate $\Delta R_{\mathrm{np}}$ in the transuranium berkelium (Bk) isotopic chain using the deformed relativistic Hartree-Bogoliubov theory in continuum (DRHBc). Our results reveal a general increase of $\Delta R_{\mathrm{np}}$ with neutron number $N$, which exhibits anti-kinks at the shell closures $N = 184, 258$ due to the shell effects. By decomposing $\Delta R_{\mathrm{np}}$ into volume and surface contributions through two-parameter Fermi (2pF) fits to angle-averaged DRHBc densities, we find that the volume term accounts for as much as $68\%$ in most nuclei, whereas the surface term dominates only near the proton drip line for $N < 142$. Nuclear deformation is shown to slightly reduce the central radius $R_c$ while significantly enhancing the surface diffuseness $a$, resulting in a notable increase in $\Delta R_{\mathrm{np}}$ that is largely driven by the surface term. Moreover, by extracting 2pF parameters along the symmetry axis ($\theta = 0^\circ$) and perpendicular to it ($\theta = 90^\circ$), we examine the anisotropy of $\Delta R_{\mathrm{np}}$. In prolate deformed nuclei, a pronounced directional dependence emerges: although the nucleus elongates along the symmetry axis, $\Delta R_{\mathrm{np}}$ is substantially larger in the perpendicular direction. This anisotropy is weak for oblate nuclei near shell closures. The anisotropy of $\Delta R_{\rm np}$ is attributed mainly to the volume term, which remains the dominant contribution in most nuclei regardless of direction. These findings provide new insights into the interplay between deformation, shell structure, and the neutron skin in finite nuclei.