---
title: M3Y-P6 Interaction in Nuclear Physics
url: https://www.emergentmind.com/topics/m3y-p6-interaction
type: topic
---

# M3Y-P6 Interaction in Nuclear Physics

The M3Y-P6 interaction is a semi-realistic effective nucleon-nucleon interaction in the M3Y-type family, developed for self-consistent mean-field calculations so as to preserve a close linkage to the bare nucleonic interaction while remaining usable for finite nuclei and nuclear matter. It derives from the Michigan-three-range-Yukawa framework based on the \(G\)-matrix, retains the tensor channels and the longest-range central channels, and supplements them with phenomenological modifications, especially density-dependent terms, to reproduce saturation, shell structure, pairing, and neutron-matter properties [1206.0794]. In subsequent work it became a reference interaction for Hartree-Fock, Hartree-Fock-Bogolyubov, quasiparticle random-phase approximation, angular-momentum projection, and related calculations of shell evolution, deformation, isotope shifts, halos, rotational energies, and positive-energy single-particle potentials [1912.05756].

## 1. Origins and design objectives

The M3Y-type interactions were developed to bridge the gap between microscopic and phenomenological treatments: they start from the microscopic M3Y \(G\)-matrix, especially the M3Y-Paris version, but allow phenomenological modification to better reproduce nuclear saturation and specific splittings [1606.03169]. Within this program, new parameter sets M3Y-P6 and M3Y-P7 were introduced by modifying the M3Y interaction while maintaining the tensor channels and the longest-range central channels [1206.0794].

The specific motivation for M3Y-P6 was broader than a fit to finite nuclei alone. Its parameters were adjusted so as to reproduce microscopic results of neutron-matter energies, the measured binding energies of doubly magic nuclei including \(^{100}\)Sn, and the even-odd mass differences of the \(Z=50\) and \(N=82\) nuclei in self-consistent mean-field calculations [1206.0794]. This construction preserved the realistic tensor force and the central one-pion-exchange component, while altering shorter-range central strengths and introducing density dependence for saturation and pairing.

This design makes M3Y-P6 “semi-realistic” in a specific sense. The tensor force \(v_{ij}^{(\mathrm{TN})}\) is kept identical to M3Y-Paris, and the longest-range central part from the original M3Y-Paris is also retained, while the density-dependent term and selected other channels are phenomenologically readjusted [1206.0794]. A plausible implication is that M3Y-P6 was intended not merely as a global fit, but as an interaction in which changes in shell structure or deformation can be traced to identifiable pieces of the effective force.

## 2. Formal structure

The effective Hamiltonian used with M3Y-P6 is
\[
H = H_N + V_C - H_{\mathrm{c.m.}},
\qquad
H_N = K + V_N,
\qquad
K = \sum_i \frac{\mathbf{p}_i^2}{2M},
\qquad
V_N = \sum_{i<j} v_{ij},
\]
with the two-body interaction decomposed as
\[
v_{ij} = v_{ij}^{(\mathrm{C})} + v_{ij}^{(\mathrm{LS})} + v_{ij}^{(\mathrm{TN})} + v_{ij}^{(\mathrm{DD})},
\]
where later presentations also write the density-dependent central term as \(v_{ij}^{(\mathrm{C}\rho)}\) [1206.0794]. In the broader M3Y framework, a density-dependent LS term \(v_{ij}^{(\mathrm{LS}\rho)}\) is introduced in the related variant M3Y-P6a rather than in the original M3Y-P6 [1606.03169].

All density-independent channels use Yukawa radial form factors,
\[
f_n^{(\mathrm{X})}(r)=\frac{e^{-\mu_n r}}{\mu_n r},
\]
and the central, spin-orbit, and tensor terms are expanded in spin-isospin channels through the usual projection operators \(P_{\mathrm{SE}}, P_{\mathrm{TE}}, P_{\mathrm{SO}}, P_{\mathrm{TO}}\) [1206.0794]. The density-dependent contact term is
\[
v_{ij}^{(\mathrm{DD})}
=
\Big(
t_\rho^{(\mathrm{SE})} P_{\mathrm{SE}} [\rho(\mathbf{r}_i)]^{\alpha^{(\mathrm{SE})}}
+
t_\rho^{(\mathrm{TE})} P_{\mathrm{TE}} [\rho(\mathbf{r}_i)]^{\alpha^{(\mathrm{TE})}}
\Big)\delta(\mathbf{r}_{ij}),
\]
with \(\alpha^{(\mathrm{SE})}=1\) and \(\alpha^{(\mathrm{TE})}=1/3\) in M3Y-P6 [1206.0794].

In the spin-orbit channel, the standard LS force in M3Y-type interactions is finite-range, and in M3Y-P6 the LS channel from the original M3Y/Paris force is multiplied by an overall factor \(2.2\) to empirically reproduce the single-particle spectra, especially in \(^{208}\)Pb [1412.1558].

| Component | Representative structure | Status in M3Y-P6 |
|---|---|---|
| Central \(v^{(\mathrm{C})}\) | Finite-range Yukawa sum in SE/TE/SO/TO channels | Longest-range central part retained |
| Spin-orbit \(v^{(\mathrm{LS})}\) | Finite-range Yukawa LS term | Original M3Y/Paris LS multiplied by \(2.2\) |
| Tensor \(v^{(\mathrm{TN})}\) | Finite-range tensor term | Identical to M3Y-Paris |
| Density-dependent central \(v^{(\mathrm{DD})}\) or \(v^{(\mathrm{C}\rho)}\) | Zero-range contact with \(\rho^\alpha\) dependence | Added for saturation and pairing |
| Density-dependent LS \(v^{(\mathrm{LS}\rho)}\) | Zero-range LS term with density dependence | Introduced in M3Y-P6a |

## 3. Calibration in nuclear matter and finite nuclei

M3Y-P6 was constructed to improve neutron-matter properties while preserving desirable finite-nucleus performance [1206.0794]. In the resulting parameter set, the isotropic spin-saturated symmetric nuclear matter remains stable in the density range as wide as \(\rho \lesssim 6\rho_0\) [1206.0794]. The same work reports the following representative nuclear-matter quantities for M3Y-P6: \(k_{F0} \approx 1.34~\mathrm{fm}^{-1}\), saturation energy \(\approx -16.24~\mathrm{MeV}\), incompressibility \(\approx 239.7~\mathrm{MeV}\), effective mass \(M_0^\ast/M \approx 0.596\), and symmetry energy \(a_{t0} \approx 32.14~\mathrm{MeV}\) [1206.0794].

On the finite-nucleus side, separation energies of proton- or neutron-magic nuclei were shown to be in fair agreement with experimental data [1206.0794]. In wider self-consistent surveys, magic numbers were identified by vanishing pair correlations in spherical Hartree-Fock-Bogolyubov calculations, with submagic numbers assigned when the energy gain due to pairing is sufficiently small [1401.8033]. Using this criterion, the predictions with M3Y-P6 were found to correspond well to known data apart from a few exceptions [1401.8033].

A later overview sharpened that assessment, stating that M3Y-P6’s predictions of magic and submagic numbers match experiments for almost all known nuclei and outperform the Gogny D1M and M3Y-P7 parameterizations [1606.03169]. In the same line of work, the interaction was presented as furnishing a new theoretical instrument for advancing nuclear mean-field approaches, precisely because realistic interaction derived from the base \(2N\) and \(3N\) interaction could be carried into a practical mean-field parametrization [1606.03169].

## 4. Shell evolution, magic numbers, and deformation

One of the defining uses of M3Y-P6 is the study of shell evolution through its explicit tensor force and its realistic spin-isospin content. In spherical calculations, the tensor-force contribution to a single-particle energy can be written as
\[
\varepsilon^{(\mathrm{TN})}(\nu)
=
2\sum_{\nu'>0} n_{\nu'}\langle\nu\nu'|v^{(\mathrm{TN})}|\nu\nu'\rangle,
\]
or, in a spherical-coupled representation,
\[
\varepsilon^{(\mathrm{TN})}(j)
=
\frac{1}{2j+1}
\sum_{j'J}
n_{j'}(2J+1)\langle jj'J|v^{(\mathrm{TN})}|jj'J\rangle,
\]
which makes explicit the orbit-occupancy dependence of the tensor effect [1604.03202; 1606.03169].

This mechanism was examined in detail in the first application of an M3Y-type interaction to deformed nuclei, namely axially symmetric constrained Hartree-Fock calculations for the \(N=20\) isotones \(^{30}\)Ne, \(^{32}\)Mg, \(^{34}\)Si and the \(N=28\) isotones \(^{40}\)Mg, \(^{42}\)Si, \(^{44}\)S [1604.03202]. The main conclusion was that the tensor force mainly causes a configuration-dependent energy shift, with little impact on the locations of the quadrupole minima themselves [1604.03202].

Around \(N=20\), the tensor force favors sphericity and acts to maintain the \(N=20\) magic number. The reason given is the \(\ell s\)-closure of \(N=20\): at sphericity the spin-saturated configuration leads to small net tensor-force effects, while intruder deformed configurations require breaking the closure by moving neutrons from \(0d_{3/2}\) to \(0f_{7/2}\), which the tensor force disfavors because protons partly fill \(0d_{5/2}\) and thus push up the neutronic \(0f_{7/2}\) [1604.03202]. In \(^{34}\)Si, the spherical minimum is always favored, whereas in \(^{30}\)Ne and \(^{32}\)Mg deformation is not suppressed entirely but the energy difference between spherical and deformed minima is increased [1604.03202].

Around \(N=28\), the pattern reverses. The tensor force facilitates deformation and acts to erode the \(N=28\) magic number because the full occupation of \(n0f_{7/2}\) produces a \(jj\)-closure with maximal repulsive tensor energy at sphericity; deformation moves the configuration toward spin saturation and reduces that repulsion [1604.03202]. In this picture, \(^{42}\)Si and \(^{40}\)Mg lose the energetic preference for sphericity, while \(^{44}\)S exhibits competing near-spherical and well-deformed minima compatible with shape coexistence [1604.03202].

These deformation studies are consistent with the broader magic-number program of M3Y-P6. The interaction reproduces the experimentally observed inversion of \(p0d_{3/2}\) and \(p1s_{1/2}\) neutron orbits from \(^{40}\)Ca to \(^{48}\)Ca and the observed trend in their spacing, whereas Skyrme and Gogny interactions, especially without tensor force, fail to do so [1606.03169]. More generally, the tensor force and the spin-isospin channel originating from the one-pion exchange potential were identified as key drivers of the \(Z\)- and \(N\)-dependence of shell gaps and hence of the appearance and disappearance of magic numbers [1401.8033].

## 5. M3Y-P6a and the density-dependent spin-orbit extension

A major development related to M3Y-P6 was the construction of M3Y-P6a, which replaces part of the enhanced two-body LS force by a density-dependent LS term motivated by three-nucleon physics. The added interaction is
\[
v_{ij}^{(\mathrm{LS}\rho)}
=
2i\, D[\rho(\mathbf{R}_{ij})]\,
\mathbf{p}_{ij}\times\delta(\mathbf{r}_{ij})\,\mathbf{p}_{ij}\cdot(\mathbf{s}_i+\mathbf{s}_j),
\]
or equivalently
\[
v_{ij}^{(\mathrm{LS}\rho)}
=
D[\rho(\mathbf{R}_{ij})]\,
\left(-\nabla_{ij}^2\delta(\mathbf{r}_{ij})\right)\,
\mathbf{L}_{ij}\cdot(\mathbf{s}_i+\mathbf{s}_j),
\]
with
\[
D[\rho(\mathbf{r})]
=
-w_1\frac{\rho(\mathbf{r})}{1+d_1\rho(\mathbf{r})}.
\]
The parameter \(w_1\) is chosen so as not to alter the empirical \(\ell s\) splitting from M3Y-P6, and \(d_1=1.0~\mathrm{fm}^3\) is used to avoid instability at high density [1412.1558].

The physical consequence of this modification is orbit dependent. With the total splitting kept constant, the density-dependent LS term tends to shrink the wave functions of the \(j=\ell+1/2\) orbits while making the \(j=\ell-1/2\) functions distribute more broadly [1412.1558]. In lead isotopes, this broadening of the neutron \(0i_{11/2}\) orbit strengthens the kink in the isotope shifts at \(N=126\); the Hartree-Fock calculations show that switching from M3Y-P6 to M3Y-P6a increases the mean-square radius of \(n0i_{11/2}\) by \(0.49~\mathrm{fm}^2\) [1412.1558].

The same mechanism was then extended to other proton-magic chains. Using spherical HFB calculations with M3Y-P6a, almost equal charge radii between \(^{40}\)Ca and \(^{48}\)Ca are reproduced, a kink is predicted at \(N=82\) in Sn isotope shifts, and the overall isotope-shift data for long Pb and Sn chains are reproduced more successfully than without the density-dependent LS term [1504.07445]. Later summaries emphasized that M3Y-P6a reproduces the kink at \(N=126\) in Pb isotope shifts even without degeneracy between the \(n1g_{9/2}\) and \(n0i_{11/2}\) levels, improves the \(N\)-dependence of Sn charge radii, and predicts \(\langle r^2\rangle_p(^{40}\mathrm{Ca}) \approx \langle r^2\rangle_p(^{48}\mathrm{Ca})\) [1606.03169].

This body of work does not replace M3Y-P6; rather, it supplements it. A plausible interpretation is that M3Y-P6 established the semi-realistic central, tensor, and standard LS framework, while M3Y-P6a incorporated a specific \(3N\)-inspired modification in the LS channel that proved decisive for certain radial observables.

## 6. Later applications and extensions

After its initial calibration and shell-structure applications, M3Y-P6 was used in a wide range of self-consistent studies. In neutron-rich magnesium isotopes, axial HFB calculations with M3Y-P6 reproduced the measured \(N\)-dependence of matter radii for \(^{34-40}\)Mg, identified a halo in \(^{37}\)Mg, and introduced the mechanism called “unpaired-particle haloing,” in which pair correlation enhances halos in odd-\(N\) nuclei [1804.08852]. The halo in \(^{37}\)Mg was predicted to have peanut shape in its intrinsic state, reflecting \(p\)-wave contribution [1804.08852].

In rotational spectroscopy, angular-momentum projection on axial Hartree-Fock solutions with M3Y-P6 showed that, except for light or weakly deformed nuclei, the ratios of the individual Hamiltonian terms to the total Peierls-Yoccoz rotational energy are insensitive to nuclides and deformation; kinetic contributions are large and close to rigid-rotor values, central-force contributions are sizable, and noncentral-force contributions are not negligible [2203.00954]. When pairing is included, the pair correlations significantly change these contributions even for well-deformed heavy nuclei [2312.14347].

For quadrupole collectivity in tin isotopes, spherical HFB plus QRPA calculations with M3Y-P6 reproduced \(E_x(2_1^+)\) and \(B(E2;0_1^+\to2_1^+)\) well in \(N\geq 64\), while constrained-HFB calculations indicated that neutron-deficient \(^{106-110}\)Sn are soft against quadrupole deformation and have almost flat potential-energy curves in the range \(|q_0|\lesssim 200~\mathrm{fm}^2\) [2303.16483]. The analysis attributed this softness to the near degeneracy of \(n0g_{7/2}\) and \(n1d_{5/2}\) together with pairing [2303.16483].

M3Y-P6 has also been carried beyond bound-state structure. A 2024 study constructed self-consistent positive-energy single-particle potentials from M3Y-P6 and showed that, when the same interaction is used for the target mean field and the single-folding potential, nucleon-nucleus elastic-scattering differential cross sections are reproduced almost comparably to empirical optical potentials up to \(80~\mathrm{MeV}\) incident energy [2403.19961]. The same paper stressed that a single energy-independent effective interaction can then generate a single-particle potential compatible with available experimental data over this range [2403.19961].

More recent applications have targeted observables with enhanced sensitivity to spin-isospin and tensor channels. Self-consistent mean-field calculations allowing time-reversal breaking found that M3Y-P6 reproduces magnetic dipole moments particularly well in nuclei adjacent to \(jj\)-closed magicity, in better agreement with data than Gogny-D1S and comparable to shell-model results with \(\chi\)EFT interaction [2508.12550]. In a different many-body framework, shell-model calculations with density-dependent interactions adapted from M3Y-P6, Gogny-D1S, and Gogny-GT2 concluded that only the M3Y-P6 functional properly describes the magicity of \(N=28\) in \(pf\)-shell nuclei [2512.07376].

## 7. Theoretical issues and limitations

Despite its successes, M3Y-P6 is not presented as a fully microscopic interaction. Review discussions emphasize that its density-dependent central and LS terms are fitted rather than derived directly from the \(G\)-matrix, that the effective mass is a bit low at about \(0.6\), and that the framework contains no explicit full \(3N\)-force term, because \(3N\) effects are folded into effective density-dependent two-body terms [1912.05756]. These are limitations internal to the semi-realistic program rather than incidental technicalities.

A more specific formal issue was identified in analyses of infinite matter and the zero-range limit of finite-range interactions. There it was argued that the spin-orbit term of the M3Y interaction is not compatible with local gauge invariance, and the recommendation was to prefer a zero-range spin-orbit term to maintain local gauge invariance and compatibility with the continuity equation [1607.00835]. The same work nevertheless found that the central part of M3Y globally reproduces Brueckner-Hartree-Fock results in all four \((S,T)\) channels up to and beyond saturation density, and that special combinations of partial waves can be used to constrain the tensor parameters [1607.00835].

These tensions define the present status of M3Y-P6. On one side, it has been repeatedly validated as a practical semi-realistic interaction that describes magic numbers, shell evolution, deformation, isotope shifts, halos, rotational energies, scattering observables, and selected electromagnetic moments across a broad mass range [1606.03169]. On the other side, its low effective mass, fitted density dependence, and the gauge-invariance issue of the finite-range LS term indicate that M3Y-P6 is best viewed as a physically informed effective interaction rather than a final reduction of the bare nuclear force.

Source: https://www.emergentmind.com/topics/m3y-p6-interaction