Charge radii of Sn isotopes in the relativistic mean field approximation

This paper investigates the kink in the charge radii of Sn isotopes around $N=82$ using the relativistic mean-field framework, revealing that while the small components of Dirac spinors for neutron states near the Fermi level significantly influence the proton central potential and contribute to the kink, they alone are insufficient to fully explain its observed magnitude.

The Gist

Imagine an atomic nucleus as a bustling city. The protons are the citizens who carry a positive charge, and the neutrons are the neutral neighbors who help hold the city together. The size of this city is measured by its charge radius—essentially, how far out the protons spread.

Scientists have been watching these cities grow as they add more neutrons. Usually, the city grows smoothly, like a balloon being inflated. But in the city of Tin (Sn), something strange happens when the population of neutrons hits exactly 82. The city doesn't just grow; it suddenly "kinks" or jolts outward. It's like a balloon that suddenly gets a weird, sharp bulge at a specific point.

This paper is a detective story about why that bulge happens, specifically looking at the city through the lens of Relativistic Mean Field (RMF) theory. Think of RMF as a high-tech, super-accurate map that accounts for Einstein's theory of relativity (where things get weird at high speeds).

Here is the breakdown of their investigation in simple terms:

1. The Mystery of the "Kink"

For decades, scientists knew about this "kink" in Tin isotopes. It's a famous puzzle.

• The Problem: Old, standard maps (non-relativistic models) couldn't draw this kink correctly. They predicted a smooth curve, missing the sharp bulge entirely.
• The Hope: The new, high-tech RMF maps were better at drawing the kink in Lead (Pb) cities, so scientists hoped they could solve the Tin mystery too.

2. The Secret Ingredient: The "Small Shadow"

In the world of quantum physics, particles like neutrons are described by mathematical objects called Dirac spinors. You can think of a neutron as having two parts:

• The Big Part (Large Component): This is the main body of the neutron, the part we usually see and count.
• The Small Part (Small Component): This is a tiny, subtle "shadow" or echo of the neutron that only appears when you look at it through the lens of relativity.

The Big Discovery:
The authors found that this "Small Shadow" is the secret sauce.

• When neutrons fill up specific "rooms" (orbitals) in the nucleus, their Small Shadows interact with the protons.
• It turns out that neutrons in certain rooms (specifically those with a specific spin direction, called $j = l - 1/2$) cast a much stronger shadow than their neighbors in the opposite room ($j = l + 1/2$).
• This strong shadow pushes the protons outward, creating the "kink" in the city's size.

The Analogy:
Imagine two people standing in a hallway.

• Person A (the standard neutron) is wearing a normal coat. They cast a normal shadow.
• Person B (the special neutron) is wearing a coat that, due to relativistic effects, casts a giant, distorted shadow that pushes the walls of the hallway outward.
• The paper shows that in Tin, the "Person B" type of neutrons are the ones causing the city to bulge.

3. The Twist: It's Not the Whole Story

The researchers were excited. They thought, "Aha! The Small Shadow explains the kink!"
But then they ran the numbers, and a problem appeared.

• The Result: While the Small Shadows did create a kink, it wasn't big enough. The theoretical kink was too small compared to what experiments actually see.
• The Real Culprit: The model failed to get the size of the Tin city before the kink (for lighter isotopes) right. Because the starting size was wrong, the final jump looked too small.

It's like trying to measure a jump. If you start from the wrong spot on the ground, even if you jump perfectly, you won't land where you expected.

4. Why This Matters

Even though the model didn't get the exact size right, the paper proves a fundamental point:

• Relativity is real and necessary. The "Small Shadow" (the small component of the Dirac spinor) is a genuine relativistic effect. Without it, you can't even get a kink at all.
• Non-relativistic models are blind. Standard physics models ignore this "Small Shadow," which is why they fail to see the kink. They are like trying to see a ghost with your eyes closed.

The Conclusion

The paper concludes that while the Relativistic Mean Field model is a powerful tool that correctly identifies why the kink happens (the interaction of the "Small Shadows" of neutrons with protons), it still needs some tuning to get the exact numbers right.

In a nutshell:
The "kink" in Tin atoms is caused by a subtle, relativistic effect where certain neutrons cast a "shadow" that pushes the atomic nucleus outward. While this explains the existence of the kink, the current maps still need a little more work to get the size of the kink perfectly right. It's a victory for understanding the mechanism, even if the final measurement isn't perfect yet.

Technical Summary

1. Problem Statement

The paper investigates the "kink effect" (KE) observed in the nuclear charge radii ($R_c$) of Tin (Sn) isotopes around the neutron magic number $N=82$.

• The Phenomenon: Experimental data shows a pronounced change in the slope of the charge radius evolution when crossing $N=82$. Specifically, the radius increases more rapidly for $N > 82$ than for $N < 82$.
• The Challenge: While the Relativistic Mean-Field Approximation (RMFA) successfully reproduces the kink in Lead (Pb) isotopes ($N=126$), it struggles to accurately reproduce the magnitude of the kink in Sn isotopes. Standard non-relativistic models (like Skyrme-Hartree-Fock) also fail to capture this effect.
• Specific Issues:
  1. Theoretical models fail to reproduce the "arch-like" behavior of charge radii for $A < 132$ (lighter than the magic nucleus $^{132}\text{Sn}$).
  2. Even with improvements, the magnitude of the kink for $A > 132$ is significantly underestimated compared to experimental data.
• Goal: To determine if the mechanisms driving the kink in Pb isotopes (specifically the role of Dirac spinor small components) apply to Sn isotopes and to clarify the origin of the discrepancy in the RMFA framework.

2. Methodology

The authors employ the Relativistic Mean-Field Approximation (RMFA) using the NL3* parameter set.

• Theoretical Framework:
  • The nucleons are described by Dirac spinors $\psi_a(r)$ satisfying the time-independent Dirac equation with scalar ($S$) and vector ($V$) mean fields.
  • The "no-sea" approximation is used (excluding nucleon-antinucleon pair creation).
  • Pairing correlations are treated using both the standard BCS approximation and the Particle-Number Projected BCS (PBCS) to handle particle-number fluctuations near closed shells more accurately.
• Analysis Techniques:
  • Schrödinger-like Equation: The Dirac equation is transformed into an equivalent Schrödinger-like equation to isolate the central potential ($V_{cent}$) and spin-orbit (SO) potential. This allows for the separation of the large ($G$) and small ($F$) components of the Dirac spinors.
  • Modified Spin-Orbit Interaction: A density-dependent SO interaction ($\tilde{V}_{SO}$) is introduced to test if reducing the SO strength in the nuclear surface region improves the results.
  • Configuration Studies: The authors artificially re-order neutron single-particle levels to test different filling sequences (e.g., $1h_{9/2}$ vs. $2f_{7/2}$) and isolate the contribution of specific orbitals.
  • Component Isolation: Calculations are performed where the small components ($F$) of the valence neutron spinors are set to zero ($F=0$) to quantify their specific contribution to the charge radius and the central potential.

3. Key Contributions and Findings

A. The Role of Small Components and Spin-Orbit Partners

The study confirms that the small components ($F$) of the Dirac spinors for valence neutrons are crucial for the formation of the kink, but their effect depends heavily on the orbital type:

• Mechanism: The small components contribute significantly to the proton central potential ($V_{cent}$). In RMFA, the scalar and vector potentials are large and opposite in sign ($|S| \sim V \gg |S-V|$). This amplifies the contribution of the small components to the potential, which is negligible in non-relativistic models.
• $j = l - 1/2$ vs. $j = l + 1/2$:
  • Neutrons in orbitals with $j = l - 1/2$ (e.g., $1h_{9/2}$) have small components that interfere constructively in the nuclear interior, creating a significant repulsive effect on the proton potential, pushing protons outward, and increasing the charge radius.
  • Neutrons in orbitals with $j = l + 1/2$ (e.g., $2f_{7/2}$) have small components that interfere destructively in the interior, resulting in a much weaker effect on the radius.
• Conclusion: The kink arises because, as neutrons fill the $N>82$ shell, the occupancy of the $1h_{9/2}$ ($j=l-1/2$) orbital increases. The large difference in the radial structure of the small components between the $1h_{9/2}$ and its spin-orbit partner $1h_{11/2}$ drives the sharp increase in radius.

B. Limitations of the Standard RMFA

Despite identifying the relativistic mechanism, the standard RMFA with NL3* fails to fully reproduce the experimental data:

• Underestimation of the Kink: The calculated kink magnitude ($\Delta^2 R_c$) is significantly smaller than experimental values.
• Failure for $A < 132$: The model cannot reproduce the arch-like behavior of charge radii for isotopes lighter than $^{132}\text{Sn}$. This failure in the $A < 132$ region is the primary reason the kink magnitude is underestimated, as the kink is defined relative to the slope change at $N=82$.
• Pairing Sensitivity: The PBCS approach yields better agreement with experiment than standard BCS for $A > 132$ because it correctly handles the occupation probabilities of the $1h_{9/2}$ and $2f_{7/2}$ orbitals, which standard BCS misrepresents (overestimating $2f_{7/2}$ and underestimating $1h_{9/2}$).

C. Impact of Modified Spin-Orbit Interaction

Introducing a density-dependent SO interaction ($\tilde{V}_{SO}$) that weakens the SO force in the nuclear surface region:

• Increases the occupancy of the $1h_{9/2}$ orbital (which is crucial for the kink).
• Improves the agreement with experimental charge radii for $A > 132$.
• However, it does not solve the problem for $A < 132$, nor does it fully close the gap in the kink magnitude.

4. Significance and Conclusions

• Relativistic Origin of the Kink: The paper establishes that the kink in Sn isotopes is a genuine relativistic effect. It stems from the specific structure of the Dirac spinors (small components) and their interaction with the proton central potential. This mechanism is intrinsic to the RMFA and is difficult to replicate in non-relativistic models without ad-hoc adjustments.
• Distinction from Pb Isotopes: While the mechanism (small components of $j=l-1/2$ orbitals) is similar to that in Pb isotopes, the Sn case is more complex because the kink involves the filling of the $1h_{9/2}$ orbital, which is close in energy to the $2f_{7/2}$ orbital. The specific interplay of these orbitals and the failure to describe the $A < 132$ region highlights limitations in the current NL3* parameterization.
• Future Directions: The authors note that the Relativistic Hartree-Fock (RHF) approach with density-dependent couplings (specifically including pion-nucleon tensor interactions) has shown better success in reproducing the Sn kink. This suggests that while the RMFA captures the qualitative relativistic mechanism (small components), a more complete treatment of the tensor force and density dependence is required for quantitative accuracy.

Summary: The paper successfully isolates the relativistic mechanism (small Dirac components) responsible for the charge radius kink in Sn isotopes but concludes that the standard RMFA is insufficient to reproduce the full experimental magnitude, primarily due to its inability to describe the charge radii of isotopes lighter than the magic nucleus $^{132}\text{Sn}$.