Skyrme Interactions in Nuclear Modeling

• Skyrme interactions are phenomenological, zero-range forces with density and momentum dependence used to model nucleon-nucleon interactions in nuclei and astrophysical environments.
• They underpin both static mean-field and advanced beyond-mean-field methods, enabling detailed calculations of nuclear structure, reaction dynamics, and correlation effects.
• Recent extensions incorporate higher-order derivatives and refined density dependencies to address ultraviolet divergences and improve simulations of neutron-star matter.

The Skyrme interaction is a phenomenological, zero-range effective force widely used to model nucleon-nucleon interactions in atomic nuclei and infinite nuclear matter. It is formulated as a local pseudopotential with explicit density and momentum dependence, supporting both static mean-field approaches and dynamic beyond-mean-field frameworks. Its parameter sets—often called Skyrme "forces" or "functionals"—nest within energy density functionals and are crucial for calculations spanning nuclear structure, reactions, and astrophysical environments.

1. General Formulation and Parametrization

The canonical Skyrme interaction encompasses central, momentum-dependent, spin-orbit, tensor, and density-dependent components. In coordinate space, the two-body and three-body terms take the following standard forms:

$$\begin{aligned} V_2(1,2) &= t_0\,(1+x_0P_\sigma)\,\delta(\mathbf r_1-\mathbf r_2) \ &+ \tfrac{1}{2}t_1\,(1+x_1P_\sigma)\,[\,\mathbf k'^2\,\delta(\mathbf r)\;+\;\delta(\mathbf r)\,\mathbf k^2\,] \ &+ t_2\,(1+x_2P_\sigma)\,\mathbf k'\!\cdot\!\delta(\mathbf r)\,\mathbf k \ V_{LS}(1,2) &= i\,W_0\,(\bm\sigma_1+\bm\sigma_2)\!\cdot\!\bigl[\mathbf k'\times\delta(\mathbf r)\,\mathbf k\bigr] \ V_3(1,2,3) &= t_3\,\delta(\mathbf r_1-\mathbf r_2)\,\delta(\mathbf r_2-\mathbf r_3) \end{aligned}$$

where $\mathbf k, \mathbf k'$ are relative-momentum operators, $P_\sigma$ is the spin exchange operator, and $t_0, t_1, t_2, t_3, W_0$ and $x_0, x_1, x_2$ are coupling constants and exchange parameters. For density-dependent effects, the three-body term is commonly mapped to an effective two-body contact interaction $t_3\,\rho^\alpha\,\delta(\mathbf r)$, with exponent $\alpha$ controlling compressibility and surface energy.

Skyrme parameter sets vary in their balance and explicit formulation of these terms, with classical examples (SII, SIII, SIV, SV, SVI) displaying systematic differences in $t_3$ (density dependence) and $t_1, t_2$ (momentum dependence) (Tohyama, 2021).

2. Treatment in Mean-Field and Beyond-Mean-Field Theories

Mean-Field Theory

In Hartree–Fock (HF) and energy density functional approaches, Skyrme interactions form the backbone for calculating nuclear ground-state properties and excitation spectra. The HF energy per nucleon in infinite matter is analytically tractable, with explicit dependence on density, momentum, and asymmetry:

$$\frac{E^{(1)}}{A}(\rho,\delta) = \frac{3\hbar^2}{10m}\Bigl(\tfrac{3\pi^2}{2}\rho\Bigr)^{2/3}G_{5/3}(\delta) + \ldots$$

with $$G_\beta(\delta) = \tfrac12\left[(1+\delta)^\beta+(1-\delta)^\beta\right]$$ (Moghrabi et al., 2012).

Beyond Mean Field and Correlations

Moving beyond the mean field, zero-range Skyrme interactions cause ultraviolet divergences in second-order many-body perturbation theory, notably scaling as $\Lambda^5$ with cutoff $\Lambda$ due to velocity-dependent terms (Moghrabi et al., 2012, Kaiser, 2015). This necessitates cutoff–regularization and simultaneous refitting of all Skyrme parameters to preserve empirical properties.

Time-dependent density-matrix theory (TDDM) utilizes the Skyrme interaction for both mean-field and residual (correlation-generating) channels, demanding a consistent mapping of the three-body term (and hence density dependence) to the residual two-body force for correlated ground states (Tohyama, 2021).

3. Density and Momentum Dependence: Effects on Correlations

Systematic investigation reveals that the magnitude of ground-state correlations induced by the Skyrme force—particularly in closed-shell nuclei such as ${}^{16}$O and ${}^{40}$Ca—depends critically on the relative strengths of density-dependent ($t_3$) and momentum-dependent ($t_1, t_2$) terms.

Param. Set	$t_3$	$t_1$	$t_2$	$\Delta E_{\rm corr}$ (${}^{16}$O)	$\Delta E_{\rm corr}$ (${}^{40}$Ca)
SII	small	large	large	–6.2 MeV	–7.5 MeV
SIV	small	large	large	–3.4 MeV	–3.6 MeV
SV	small	large	large	–7.1 MeV	–5.6 MeV
SIII	int.	int.	int.	–14.1 MeV	–22.4 MeV
SVI	large	small	small	–28.0 MeV	–56.5 MeV

Sets with strong momentum dependence and weak density dependence yield minimal correlation energies, while those with dominant density dependence (notably SVI) induce robust correlations, aligning with shell-model and perturbative benchmarks (Tohyama, 2021). This reflects cancellation of attractive $t_0,t_3$ and repulsive $t_1,t_2$ in high-momentum transfer residual interactions, crucial for beyond-HF consistency. A modest downward readjustment of SVI restores HF binding energies without suppressing correlations.

4. Applications in Nuclear Structure, Reactions, and Astrophysics

Nuclear Structure

Skyrme functionals underpin static and dynamical mean-field calculations:

• Accurate reproduction of binding energies, charge radii, and spin-orbit splittings across semi-magic nuclei.
• Treatment of collective excitation modes (GMR, GDR, Gamow–Teller resonance), where improved spin-isospin and tensor structures (as in the SAMi family) yield sum rule exhaustion and realistic response centroids (Roca-Maza et al., 2012, Zheng et al., 2016).
• Implementation of density-dependent spin–orbit terms (Skyrme–ddso) enhances the charge-radius "kink" at N=126 in Pb isotopes and resolves surface density discrepancies (Kanada-En'yo, 2022).

Nuclear Reactions

Time-dependent Hartree–Fock (TDHF) simulations establish the form and parametrization of the Skyrme interaction as critical to barrier heights, fusion cross sections, and dissipation in heavy-ion collisions. Spin–orbit and tensor components decisively affect fusion thresholds and damping patterns (Stevenson et al., 2018). Semiclassical ETF2-based models confirm SVI, SII, and SIII as optimally accurate for fusion-barrier systematics (Ghodsi et al., 2015).

Astrophysics

The extension to neutron-star matter and core-collapse environments relies on Skyrme EOS predictions. Extended functionals (with advanced momentum and density dependence, e.g. $t_4,t_5$ terms) eliminate finite-density instabilities, fit both finite-nucleus and neutron-star constraints, and match observed mass-radius curves and Urca thresholds (Zhang et al., 2015, Duan et al., 2024, Wang et al., 2023, Wang et al., 2024). New families allow systematic variation of symmetry energy slope ($L$), effective mass splitting, and high-momentum behavior, supporting unified modeling across laboratory and astrophysical regimes.

5. Advanced Extensions: High-Momentum and Flexible Density Terms

Recent work generalizes the Skyrme pseudopotential to next-to-next-to-next-to-leading order (N$^3$LO), N$^4$LO, and N$^5$LO forms, incorporating derivative operators up to tenth order and density-dependent terms as Fermi-momentum expansions. This yields:

• Saturation of nucleon optical potentials up to 2 GeV.
• Enhanced ability to match both empirical nuclear optical data and neutron-star observables.
• Systematic control over the symmetry energy ($E_{\rm sym}$), slope ($L$), and its higher derivatives.
• Efficient implementation in transport models for heavy-ion collisions, accommodating the required momentum dependence for multi-GeV dynamics (Wang et al., 2018, Wang et al., 2024).

6. Open Issues: Regularization, Tensor Breakdown, Future Directions

Ultraviolet Divergences and Regularization

Zero-range Skyrme interactions inherently produce divergences in second and higher-order MBPT calculations. Cutoff regularization and refitting per Lambda have been successfully implemented, but a unique, cutoff-independent scheme (e.g. dimensional regularization) remains desirable for universality (Moghrabi et al., 2012, Kaiser, 2015).

Tensor Components and Spin–Isospin Channels

It is demonstrated that the standard zero-range tensor force contributes zero energy at the mean-field level and is thus invalid in conventional Skyrme functionals; central spin–spin terms can partially compensate by tuning shell evolutions, especially in open-shell chains (Dong et al., 2020). Optimized parameterizations improve Gamow–Teller and spin-dipole resonances, critical for weak-interaction and neutrino processes (Roca-Maza et al., 2012, Zheng et al., 2016).

Future Prospects

Outstanding endeavors include:

• Improving the matching between Skyrme effective masses and non-relativistic BHF results at high density.
• Global re-optimization of traditional Skyrme (and Skyrme–ddso) parameters in tandem with new density-dependent gradient and spin–orbit terms.
• Further developing extended pseudopotentials to N$^n$LO for transport codes and EOS calculations in ultradense and highly asymmetric matter.

7. Summary Table: Representative Skyrme Parametrizations and Key Attributes

Name	Density-dep. ($t_3$, $\alpha$)	Momentum-dep. ($t_1$, $t_2$)	Notable Feature	Ref.
SVI	Large $t_3$, $\alpha$	Small $t_1$, $t_2$	Strong correlations	(Tohyama, 2021)
SII	Small $t_3$, $\alpha$	Large $t_1$, $t_2$	Weak correlations	(Tohyama, 2021)
SIII	Intermediate	Intermediate	Moderate correlations	(Tohyama, 2021)
SAMi	Refined spin-isospin, $J^2$	Realistic GTR/SDR/IAR	Spin–isospin accuracy	(Roca-Maza et al., 2012)
Sky3	Extended $t_4$, $t_5$ ($\rho^{1/3}$)	Non-monotonic $m^*(\rho)$	EOS & neutron stars	(Duan et al., 2024)
eMSL07/08/09	$t_4$, $t_5$ for 3N momentum	Eliminates instabilities	Heavy neutron stars	(Zhang et al., 2015)
N3LO...N5LO	Up to $p^{10}$ momentum-dep.	Fermi momentum expansion	2 GeV transport validity	(Wang et al., 2024)

The Skyrme interaction, in both standard and extended forms, remains the backbone of nuclear EDF theory. Its ongoing development—driven by regularization, improved spin–isospin properties, flexible density/momentum dependence, and integration with transport models—continues to address longstanding challenges in nuclear structure, reactions, and compact-object physics.