Neutron skin thickness and its volume and surface contributions

Using the deformed relativistic Hartree-Bogoliubov theory in continuum, this study systematically analyzes the neutron skin thickness in berkelium isotopes, revealing that while shell closures induce anti-kinks and deformation enhances surface diffuseness, the volume term remains the dominant contributor to the skin thickness, which exhibits significant anisotropy in prolate nuclei.

The Gist

The Big Picture: The "Fuzzy" Edge of an Atomic Nucleus

Imagine an atomic nucleus not as a hard marble, but as a soft, fuzzy ball of dough. Inside this dough, there are two types of ingredients: protons (which are positively charged) and neutrons (which are neutral).

Usually, the protons and neutrons are mixed together pretty evenly in the center. However, in heavy, unstable atoms (like the ones studied in this paper, called Berkelium), the neutrons start to pile up on the outside, creating a "skin" of extra neutrons. This is called the neutron skin.

The thickness of this skin is a huge deal for scientists. It acts like a "thermometer" for the laws of physics that govern how matter behaves inside neutron stars and during supernova explosions. If we can measure how thick this skin is, we can understand the "stiffness" of the nuclear force.

What Did These Scientists Do?

The researchers used a super-computer model called DRHBc (a fancy way of saying they simulated how these fuzzy balls behave when they are squashed or stretched). They looked at a whole chain of Berkelium atoms, adding more and more neutrons to see how the skin changed.

Here are their three main discoveries, explained simply:

1. The "Anti-Kink" Surprise

As you add more neutrons to the atom, the skin generally gets thicker, just like adding more frosting to a cake makes the layer thicker.

• The Twist: However, when the number of neutrons hits specific "magic numbers" (184 and 258), the skin suddenly stops growing as fast. It's like hitting a speed bump.
• Why? At these magic numbers, the neutrons fill up a perfect, stable shell (like a full parking lot). This stability makes the nucleus resist changing shape, causing a temporary pause in the skin's growth.

2. The "Volume" vs. "Surface" Debate

The scientists wanted to know why the skin gets thicker. Is it because the whole ball of dough is getting bigger (Volume), or is it because the fuzzy edge is getting fluffier and more spread out (Surface)?

• The Finding: For most of these atoms, the skin gets thicker because the whole ball is expanding (the Volume contribution). This accounts for about 68% of the skin's thickness.
• The Exception: Only for the very lightest atoms in their study (near the "proton drip line," where the nucleus is barely holding together) does the "fluffiness" of the edge (Surface) become the main reason for the skin.
• The Deformation Effect: Many of these atoms aren't perfect spheres; they are squashed like a rugby ball (prolate) or a pancake (oblate). The study found that when an atom is deformed, it doesn't get much bigger in the center, but its edge becomes much fluffier. This extra fluffiness is what makes the skin thicker in deformed atoms.

3. The "Directional" Skin (The Anisotropy)

This is the most surprising part. Because these atoms are squashed (deformed), the neutron skin isn't the same thickness in every direction.

• The Analogy: Imagine a rugby ball (a prolate nucleus). It is long from top to bottom and short from side to side.
• The Counter-Intuitive Result: You might think the skin would be thickest where the ball is longest (top to bottom). But it's the opposite!
  • The neutron skin is actually thicker on the sides (perpendicular to the long axis) than it is at the tips.
  • Even though the nucleus is stretched out along the long axis, the "fuzz" of neutrons spreads out more on the sides.
• Why? It turns out the "Volume" part of the skin (the main bulk) is responsible for this difference. The way the neutrons and protons are packed inside creates a situation where the skin is naturally thicker on the "equator" of the rugby ball than at the "poles."

Summary in a Nutshell

1. Neutron skins are fuzzy layers of extra neutrons on heavy atoms.
2. As you add neutrons, the skin gets thicker, but it pauses at "magic numbers" where the nucleus is extra stable.
3. The skin gets thicker mostly because the whole nucleus expands, not just because the edge gets fluffier (except for the lightest atoms).
4. Deformation matters: Squashing the nucleus makes the edge fluffier, which thickens the skin.
5. Direction matters: In squashed (rugby-ball shaped) atoms, the neutron skin is surprisingly thicker on the sides than at the tips, driven mostly by how the neutrons are packed inside.

This research helps scientists understand the rules of the universe that apply to everything from the heaviest elements we can make in a lab to the densest stars in the sky.

Technical Summary

Technical Summary: Neutron Skin Thickness and its Volume and Surface Contributions in Transuranium Berkelium Isotopes

Problem Statement
Accurate determination of the neutron skin thickness ($\Delta R_{np}$), defined as the difference between the root-mean-square radii of neutron and proton density distributions, is critical for constraining the density dependence of the nuclear symmetry energy. While $\Delta R_{np}$ is well-studied in spherical nuclei like $^{208}\text{Pb}$, its behavior in deformed, transuranium nuclei remains less explored. Specifically, the relative contributions of volume and surface effects to the skin thickness, and the potential anisotropy of the skin in deformed systems, require systematic investigation. Furthermore, experimental extraction of neutron radii is challenging and often model-dependent, necessitating robust theoretical frameworks to provide reliable predictions.

Methodology
This study employs the Deformed Relativistic Hartree-Bogoliubov theory in continuum (DRHBc) to systematically investigate the neutron skin thickness in the transuranium berkelium (Bk) isotopic chain. The DRHBc framework provides a unified description of nuclear deformation, pairing correlations, and continuum effects.

• Theoretical Framework: The calculations utilize the PC-PK1 density functional for the particle-hole channel and a zero-range pairing force for the particle-particle channel. The Dirac Woods-Saxon (DWS) basis is used to solve the DRHB equations, properly treating continuum effects.
• Decomposition Strategy: To separate $\Delta R_{np}$ into volume and surface contributions, the angle-averaged DRHBc density distributions are fitted using a two-parameter Fermi (2pF) distribution via the Levenberg-Marquardt algorithm. The neutron skin is decomposed using the surface expansion formalism, where $\Delta R_{np} = \Delta R_{np}^{\text{vol}} + \Delta R_{np}^{\text{surf}}$. The volume term is derived from the difference in equivalent sharp radii ($R_n - R_p$), while the surface term arises from differences in surface diffuseness.
• Anisotropy Analysis: To examine directional dependence, 2pF parameters are extracted separately along the symmetry axis ($\theta = 0^\circ$) and perpendicular to it ($\theta = 90^\circ$).
• Comparative Analysis: Results from the deformed DRHBc calculations are compared with those from spherical Relativistic Continuum Hartree-Bogoliubov (RCHB) calculations to quantify the specific influence of nuclear deformation.

Key Contributions and Results

1. Evolution of Neutron Skin: The study confirms that $\Delta R_{np}$ generally increases with neutron number $N$ across the Bk isotopic chain. However, "anti-kinks" (local decreases or flattening of the growth trend) are observed at the neutron shell closures $N = 184$ and $N = 258$, attributed to shell effects.
2. Volume vs. Surface Dominance: By decomposing the skin thickness, the authors find that the volume term ($\Delta R_{np}^{\text{vol}}$) is the dominant contributor in most nuclei, accounting for up to 68% of the total thickness in neutron-rich isotopes. The surface term ($\Delta R_{np}^{\text{surf}}$) only becomes the dominant contribution near the proton drip line for nuclei with $N < 142$.
3. Impact of Deformation: Nuclear deformation significantly enhances the neutron skin thickness compared to spherical calculations. This enhancement is primarily driven by the surface term. Deformation slightly reduces the central radius ($R_c$) and central density ($\rho_0$) but significantly increases the surface diffuseness parameter ($a$), leading to a larger overall skin.
4. Anisotropy of the Neutron Skin: A pronounced anisotropy is identified in prolate deformed nuclei. Despite the nucleus being elongated along the symmetry axis ($\theta = 0^\circ$), the neutron skin thickness is substantially larger in the perpendicular direction ($\theta = 90^\circ$). This anisotropy is weak in oblate nuclei near shell closures.
   • Origin of Anisotropy: The directional dependence is attributed mainly to the volume term. Although the equivalent sharp radii ($R_n$ and $R_p$) are larger along the symmetry axis, the difference ($R_n - R_p$) is smaller there than in the perpendicular direction. The surface term also contributes to this anisotropy, with diffuseness being larger in the perpendicular direction for prolate nuclei.
5. Shell Effects on Anisotropy: Near magic numbers ($N=184, 258$), where deformation is minimal, the anisotropy of $\Delta R_{np}$ is significantly reduced, though a small residual anisotropy persists due to slight residual deformations.

Significance
The paper claims to provide new insights into the interplay between nuclear deformation, shell structure, and the neutron skin in finite nuclei. By systematically decomposing the neutron skin into volume and surface components for deformed systems, the study highlights that deformation-induced enhancements in skin thickness are largely driven by surface effects. Furthermore, the discovery of significant anisotropy in the neutron skin thickness—specifically that the skin is thicker perpendicular to the symmetry axis in prolate nuclei—challenges simple geometric intuitions and underscores the necessity of considering directional dependence in deformed nuclei. These findings offer refined constraints for the density dependence of the nuclear symmetry energy, which has broad implications for understanding exotic nuclei, nuclear masses, and astrophysical phenomena such as neutron star structure.