Centrifugal-corrected harmonic oscillator model for spherical proton emitters

This paper proposes a centrifugal-corrected harmonic oscillator model, calibrated with experimental data and integrated with relativistic mean field theory, to accurately predict and analyze proton radioactivity half-lives in spherical nuclei while providing a robust framework for future nuclear structure research.

The Gist

Imagine an atomic nucleus as a crowded, chaotic dance floor. Usually, the dancers (protons and neutrons) hold hands tightly, forming a stable circle. But sometimes, the floor gets too crowded with one type of dancer (protons), and the circle becomes unstable. To regain balance, the nucleus might kick one dancer out. This ejection of a proton is called proton radioactivity.

This paper is like a team of physicists trying to build a better crystal ball to predict exactly how long it takes for that proton to escape.

Here is the breakdown of their work using simple analogies:

1. The Problem: The "Slippery Slope"

In physics, we know that for a proton to escape, it has to tunnel through a "wall" of energy (a barrier). Think of this like a ball trying to roll out of a deep bowl.

• The Old Model: Previous models tried to predict how fast the ball rolls out using a standard "harmonic oscillator" (like a perfect spring). It was okay, but it missed a crucial detail.
• The Missing Piece: When the proton is kicked out, it often spins (has orbital angular momentum). Imagine the proton isn't just rolling straight out; it's spinning like a top. This spin creates a "centrifugal force" (like the force pushing you to the side when a car turns sharply). This force makes the energy wall higher and harder to climb.
• The Issue: The old models ignored this spinning effect. For protons, this spin matters a lot more than it does for other types of radioactive decay. Ignoring it was like trying to predict a race time without knowing if the runner was wearing heavy boots.

2. The Solution: The "Centrifugal-Corrected" Model

The authors, led by Xiao-Yan Zhu, built a new, improved model. They added a specific "correction factor" to their equations to account for that spinning force.

• The Analogy: Imagine you are calculating how long it takes a skier to slide down a hill.
  • Old Model: You only looked at the steepness of the hill.
  • New Model: You realized the skier is also doing a triple-twist spin in the air. You added a new variable to your math to account for the extra energy needed to spin while sliding.
• The Result: They found a specific number (a "centrifugal parameter") that, when plugged into their formula, made their predictions match real-world experiments almost perfectly. They could now predict the "escape time" (half-life) within a factor of 2.4, which is a huge improvement in the world of nuclear physics.

3. The "Spectroscopic Factor": The "Ticket" to the Dance Floor

To calculate the time, you also need to know how likely the proton is to even try to leave. This is called the spectroscopic factor.

• The Analogy: Think of the nucleus as a theater. The proton is an actor wanting to leave the stage. But the actor can only leave if their seat (orbit) is empty in the "daughter" nucleus (the group left behind).
• The Method: The authors used a sophisticated computer simulation (called Relativistic Mean Field theory) to check the "seating chart" of the nucleus. They calculated the probability that the proton's seat was actually empty and ready for it to exit. This made their prediction much more accurate than just guessing.

4. The "Secret Code": Connecting Structure to Escape

The team discovered a beautiful, simple relationship between the structure of the nucleus and the difficulty of escaping.

• The Analogy: They found a linear "code." If you know how "stiff" the energy wall is (the fragmentation potential), you can predict how "wide" the door is for the proton to escape.
• Why it matters: This confirms that the internal architecture of the atom (how the protons and neutrons are arranged) is directly linked to how they behave when they try to escape. It's like realizing that the shape of a keyhole determines exactly how fast a key can turn.

5. The Future: Predicting the Unknown

Finally, they used their new, super-accurate crystal ball to look into the future.

• The Application: There are many unstable atoms that scientists know exist (or suspect exist) but haven't measured yet because they are so rare or hard to catch.
• The Prediction: The authors used their model to predict the "escape times" for these mysterious atoms. They checked their predictions against other known laws and found they matched up well. This gives other scientists a reliable map to guide their future experiments in finding and studying these exotic elements.

Summary

In short, this paper is about fixing a broken map.

1. The Map: A mathematical model to predict how long unstable atoms last before spitting out a proton.
2. The Flaw: The old map ignored the "spin" of the proton, making it inaccurate.
3. The Fix: They added a "spin correction" and a better way to calculate the proton's "ticket" to leave.
4. The Reward: They now have a tool that predicts these events with high precision, helping us understand the very limits of how matter can exist.

Technical Summary

Here is a detailed technical summary of the paper "Centrifugal-corrected harmonic oscillator model for spherical proton emitters."

1. Problem Statement

Proton radioactivity is a quantum tunneling phenomenon where unstable, proton-rich nuclei emit protons to reach a more stable state. While various theoretical models exist (e.g., Universal Decay Law, Generalized Liquid Drop Model), accurately predicting proton emission half-lives ($T_{1/2}$) remains challenging due to two primary factors:

1. Spectroscopic Factors ($S_p$): The probability of the emitted proton pre-formation within the nucleus is difficult to calculate precisely using phenomenological methods.
2. Centrifugal Potential Sensitivity: Unlike $\alpha$-decay or cluster radioactivity, proton emission half-lives are extremely sensitive to the orbital angular momentum ($l$) of the emitted proton. Existing models often fail to adequately correct for the centrifugal barrier, leading to significant deviations (orders of magnitude) for high-$l$ transitions.

The authors aim to develop a more robust theoretical framework that systematically evaluates proton radioactivity half-lives for spherical nuclei by integrating advanced nuclear structure calculations with a refined treatment of centrifugal effects.

2. Methodology

The authors propose an Improved Harmonic Oscillator Model (HOPM) that incorporates three key technical components:

A. Theoretical Framework for Half-Life

The half-life is calculated using the standard relation $T_{1/2} = \hbar \ln 2 / \Gamma$, where the decay width $\Gamma$ depends on the spectroscopic factor ($S_p$), a normalization factor ($F$), and the barrier penetration probability ($P$).

• Potential Model: The nuclear potential is modeled using a modified harmonic oscillator potential ($V_N = -V_0 + V_1 r^2$) truncated at the classical turning point, combined with Coulomb and centrifugal potentials.
• Centrifugal Correction: A new term, $d \cdot l(l+1)$, is added to the logarithmic half-life equation to account for the centrifugal barrier effects, which are more pronounced in proton emission than in $\alpha$-decay.

B. Calculation of Spectroscopic Factors ($S_p$)

Instead of using phenomenological fits, the authors calculate $S_p$ using a microscopic approach:

• Relativistic Mean Field (RMF) Theory: Combined with the BCS (Bardeen-Cooper-Schrieffer) method to handle pairing correlations.
• Interaction: The DD-ME2 effective interaction is used to determine the single-particle structure and the probability ($u_j^2$) that the proton orbit is empty in the daughter nucleus.

C. Analytical Derivations

• Fragmentation Potential ($V_{frag}$): The authors verify a linear relationship between the logarithm of the reduced width ($\log_{10} \gamma^2$) and the fragmentation potential ($V_{frag} = V_C - Q_p$).
• Centrifugal Parameter ($d$): An analytical expression for the centrifugal parameter is derived based on the harmonic oscillator approximation: $d_{Ae} \approx \frac{\hbar}{\mu r_1^2 \omega \ln 10}$. This is compared against an empirically fitted value.

3. Key Contributions

1. Centrifugal-Corrected Model: The introduction of the adjustable parameter $d$ in the term $d \cdot l(l+1)$ specifically addresses the high sensitivity of proton emission to orbital angular momentum, a limitation in previous HOPM applications.
2. Microscopic Spectroscopic Factors: The integration of RMF+BCS (DD-ME2) provides a more physically grounded calculation of $S_p$ compared to previous semi-microscopic or purely phenomenological estimates.
3. Structure-Dynamics Connection: The paper establishes and verifies an analytical link between nuclear structure (via fragmentation potential) and tunneling dynamics (via reduced width), confirming that nuclear structure effects directly influence the decay probability.
4. Predictive Framework: The model is extended to predict half-lives for proton emission candidates in the NUBASE2020 database that are energetically allowed but not yet experimentally quantified.

4. Results

• Parameter Optimization: By fitting experimental data for 32 spherical nuclei, the authors obtained:
  • Nuclear potential depth: $V_0 = 62.4$ MeV.
  • Centrifugal parameter: $d = 0.143$.
• Accuracy Improvement:
  • The improved model (Eq. 17) reduces the standard deviation ($\sigma$) between calculated and experimental half-lives to 0.380.
  • This outperforms existing models: UDL ($\sigma=0.403$), N-GNL ($\sigma=0.525$), and UFM ($\sigma=0.705$).
  • The model controls errors within a factor of 2.4 for experimental data.
• Centrifugal Effect Validation: For high angular momentum cases (e.g., $l=5$), the uncorrected model (Eq. 16) underestimated half-lives by nearly 6 orders of magnitude. The inclusion of the $d \cdot l(l+1)$ term corrected this, bringing predictions in line with experimental data (e.g., for $^{155}$Ta, the deviation improved from $\Delta = -0.223$ to $\Delta = -0.16$).
• Analytical vs. Empirical $d$: The analytically derived $d_{Ae} \approx 0.167$ yielded larger deviations ($\sigma=0.511$) compared to the fitted $d=0.143$. The authors attribute this to simplifications in the harmonic oscillator assumption (neglecting deformation and higher-order corrections), validating the necessity of empirical fitting for high precision.
• Predictions: The model successfully predicts half-lives for candidates like $^{111}$Cs, $^{169}$Ir$^m$, and $^{172}$Au$^m$, showing consistency with the modified Geiger-Nuttall law.

5. Significance

This work significantly advances the theoretical understanding of proton radioactivity by:

• Enhancing Precision: Providing a robust tool that achieves better agreement with experimental data than current standard models, particularly for high-$l$ transitions.
• Bridging Theory and Experiment: Demonstrating that microscopic nuclear structure calculations (RMF+BCS) can be effectively coupled with tunneling models to predict decay rates.
• Guiding Future Research: Offering reliable predictions for unobserved proton emitters near the proton drip line, which is crucial for guiding future experiments in nuclear physics and understanding the limits of nuclear stability.
• Framework Extension: Laying the groundwork for future studies on deformed nuclei and multi-proton emission systems.