M3Y-P6 Interaction in Nuclear Physics

• M3Y-P6 interaction is a semi-realistic effective nucleon-nucleon force derived from the Michigan-three-range-Yukawa framework, incorporating density-dependent terms and a realistic tensor force.
• It is calibrated to reproduce key properties of finite nuclei and nuclear matter, accurately modeling saturation energy, shell structure, pairing gaps, and magic numbers.
• The interaction drives insights in nuclear structure studies, influencing shell evolution, deformation, isotope shifts, and nucleon scattering, and has been extended via M3Y-P6a.

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 (Nakada, 2012). 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 (Nakada, 2019).

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 (Nakada et al., 2016). 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 (Nakada, 2012).

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 (Nakada, 2012). 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 (Nakada, 2012). 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)}$ (Nakada, 2012). 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 (Nakada et al., 2016).

All density-independent channels use Yukawa radial form factors,

$G$0

and the central, spin-orbit, and tensor terms are expanded in spin-isospin channels through the usual projection operators $G$1 (Nakada, 2012). The density-dependent contact term is

$G$2

with $G$3 and $G$4 in M3Y-P6 (Nakada, 2012).

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 $G$5 to empirically reproduce the single-particle spectra, especially in $G$6Pb (Nakada et al., 2014).

Component	Representative structure	Status in M3Y-P6
Central $G$7	Finite-range Yukawa sum in SE/TE/SO/TO channels	Longest-range central part retained
Spin-orbit $G$8	Finite-range Yukawa LS term	Original M3Y/Paris LS multiplied by $G$9
Tensor $^{100}$0	Finite-range tensor term	Identical to M3Y-Paris
Density-dependent central $^{100}$1 or $^{100}$2	Zero-range contact with $^{100}$3 dependence	Added for saturation and pairing
Density-dependent LS $^{100}$4	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 (Nakada, 2012). In the resulting parameter set, the isotropic spin-saturated symmetric nuclear matter remains stable in the density range as wide as $^{100}$5 (Nakada, 2012). The same work reports the following representative nuclear-matter quantities for M3Y-P6: $^{100}$6, saturation energy $^{100}$7, incompressibility $^{100}$8, effective mass $^{100}$9, and symmetry energy $Z=50$0 (Nakada, 2012).

On the finite-nucleus side, separation energies of proton- or neutron-magic nuclei were shown to be in fair agreement with experimental data (Nakada, 2012). 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 (Nakada et al., 2014). Using this criterion, the predictions with M3Y-P6 were found to correspond well to known data apart from a few exceptions (Nakada et al., 2014).

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 (Nakada et al., 2016). 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 $Z=50$1 and $Z=50$2 interaction could be carried into a practical mean-field parametrization (Nakada et al., 2016).

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

$Z=50$3

or, in a spherical-coupled representation,

$Z=50$4

which makes explicit the orbit-occupancy dependence of the tensor effect (Suzuki et al., 2016, Nakada et al., 2016).

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 $Z=50$5 isotones $Z=50$6Ne, $Z=50$7Mg, $Z=50$8Si and the $Z=50$9 isotones $N=82$0Mg, $N=82$1Si, $N=82$2S (Suzuki et al., 2016). 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 (Suzuki et al., 2016).

Around $N=82$3, the tensor force favors sphericity and acts to maintain the $N=82$4 magic number. The reason given is the $N=82$5-closure of $N=82$6: 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 $N=82$7 to $N=82$8, which the tensor force disfavors because protons partly fill $N=82$9 and thus push up the neutronic $v_{ij}^{(\mathrm{TN})}$0 (Suzuki et al., 2016). In $v_{ij}^{(\mathrm{TN})}$1Si, the spherical minimum is always favored, whereas in $v_{ij}^{(\mathrm{TN})}$2Ne and $v_{ij}^{(\mathrm{TN})}$3Mg deformation is not suppressed entirely but the energy difference between spherical and deformed minima is increased (Suzuki et al., 2016).

Around $v_{ij}^{(\mathrm{TN})}$4, the pattern reverses. The tensor force facilitates deformation and acts to erode the $v_{ij}^{(\mathrm{TN})}$5 magic number because the full occupation of $v_{ij}^{(\mathrm{TN})}$6 produces a $v_{ij}^{(\mathrm{TN})}$7-closure with maximal repulsive tensor energy at sphericity; deformation moves the configuration toward spin saturation and reduces that repulsion (Suzuki et al., 2016). In this picture, $v_{ij}^{(\mathrm{TN})}$8Si and $v_{ij}^{(\mathrm{TN})}$9Mg lose the energetic preference for sphericity, while $$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},$$0S exhibits competing near-spherical and well-deformed minima compatible with shape coexistence (Suzuki et al., 2016).

These deformation studies are consistent with the broader magic-number program of M3Y-P6. The interaction reproduces the experimentally observed inversion of $$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},$$1 and $$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},$$2 neutron orbits from $$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},$$3Ca to $$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},$$4Ca and the observed trend in their spacing, whereas Skyrme and Gogny interactions, especially without tensor force, fail to do so (Nakada et al., 2016). More generally, the tensor force and the spin-isospin channel originating from the one-pion exchange potential were identified as key drivers of the $$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},$$5- and $$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},$$6-dependence of shell gaps and hence of the appearance and disappearance of magic numbers (Nakada et al., 2014).

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

$$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},$$7

or equivalently

$$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},$$8

with

$$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},$$9

The parameter $$v_{ij} = v_{ij}^{(\mathrm{C})} + v_{ij}^{(\mathrm{LS})} + v_{ij}^{(\mathrm{TN})} + v_{ij}^{(\mathrm{DD})},$$0 is chosen so as not to alter the empirical $$v_{ij} = v_{ij}^{(\mathrm{C})} + v_{ij}^{(\mathrm{LS})} + v_{ij}^{(\mathrm{TN})} + v_{ij}^{(\mathrm{DD})},$$1 splitting from M3Y-P6, and $$v_{ij} = v_{ij}^{(\mathrm{C})} + v_{ij}^{(\mathrm{LS})} + v_{ij}^{(\mathrm{TN})} + v_{ij}^{(\mathrm{DD})},$$2 is used to avoid instability at high density (Nakada et al., 2014).

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 $$v_{ij} = v_{ij}^{(\mathrm{C})} + v_{ij}^{(\mathrm{LS})} + v_{ij}^{(\mathrm{TN})} + v_{ij}^{(\mathrm{DD})},$$3 orbits while making the $$v_{ij} = v_{ij}^{(\mathrm{C})} + v_{ij}^{(\mathrm{LS})} + v_{ij}^{(\mathrm{TN})} + v_{ij}^{(\mathrm{DD})},$$4 functions distribute more broadly (Nakada et al., 2014). In lead isotopes, this broadening of the neutron $$v_{ij} = v_{ij}^{(\mathrm{C})} + v_{ij}^{(\mathrm{LS})} + v_{ij}^{(\mathrm{TN})} + v_{ij}^{(\mathrm{DD})},$$5 orbit strengthens the kink in the isotope shifts at $$v_{ij} = v_{ij}^{(\mathrm{C})} + v_{ij}^{(\mathrm{LS})} + v_{ij}^{(\mathrm{TN})} + v_{ij}^{(\mathrm{DD})},$$6; the Hartree-Fock calculations show that switching from M3Y-P6 to M3Y-P6a increases the mean-square radius of $$v_{ij} = v_{ij}^{(\mathrm{C})} + v_{ij}^{(\mathrm{LS})} + v_{ij}^{(\mathrm{TN})} + v_{ij}^{(\mathrm{DD})},$$7 by $$v_{ij} = v_{ij}^{(\mathrm{C})} + v_{ij}^{(\mathrm{LS})} + v_{ij}^{(\mathrm{TN})} + v_{ij}^{(\mathrm{DD})},$$8 (Nakada et al., 2014).

The same mechanism was then extended to other proton-magic chains. Using spherical HFB calculations with M3Y-P6a, almost equal charge radii between $$v_{ij} = v_{ij}^{(\mathrm{C})} + v_{ij}^{(\mathrm{LS})} + v_{ij}^{(\mathrm{TN})} + v_{ij}^{(\mathrm{DD})},$$9Ca and $v_{ij}^{(\mathrm{C}\rho)}$0Ca are reproduced, a kink is predicted at $v_{ij}^{(\mathrm{C}\rho)}$1 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 (Nakada, 2015). Later summaries emphasized that M3Y-P6a reproduces the kink at $v_{ij}^{(\mathrm{C}\rho)}$2 in Pb isotope shifts even without degeneracy between the $v_{ij}^{(\mathrm{C}\rho)}$3 and $v_{ij}^{(\mathrm{C}\rho)}$4 levels, improves the $v_{ij}^{(\mathrm{C}\rho)}$5-dependence of Sn charge radii, and predicts $v_{ij}^{(\mathrm{C}\rho)}$6 (Nakada et al., 2016).

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 $v_{ij}^{(\mathrm{C}\rho)}$7-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 $v_{ij}^{(\mathrm{C}\rho)}$8-dependence of matter radii for $v_{ij}^{(\mathrm{C}\rho)}$9Mg, identified a halo in $v_{ij}^{(\mathrm{LS}\rho)}$0Mg, and introduced the mechanism called “unpaired-particle haloing,” in which pair correlation enhances halos in odd-$v_{ij}^{(\mathrm{LS}\rho)}$1 nuclei (Nakada et al., 2018). The halo in $v_{ij}^{(\mathrm{LS}\rho)}$2Mg was predicted to have peanut shape in its intrinsic state, reflecting $v_{ij}^{(\mathrm{LS}\rho)}$3-wave contribution (Nakada et al., 2018).

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 (Abe et al., 2022). When pairing is included, the pair correlations significantly change these contributions even for well-deformed heavy nuclei (Abe et al., 2023).

For quadrupole collectivity in tin isotopes, spherical HFB plus QRPA calculations with M3Y-P6 reproduced $v_{ij}^{(\mathrm{LS}\rho)}$4 and $v_{ij}^{(\mathrm{LS}\rho)}$5 well in $v_{ij}^{(\mathrm{LS}\rho)}$6, while constrained-HFB calculations indicated that neutron-deficient $v_{ij}^{(\mathrm{LS}\rho)}$7Sn are soft against quadrupole deformation and have almost flat potential-energy curves in the range $v_{ij}^{(\mathrm{LS}\rho)}$8 (Omura et al., 2023). The analysis attributed this softness to the near degeneracy of $v_{ij}^{(\mathrm{LS}\rho)}$9 and $G$00 together with pairing (Omura et al., 2023).

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 $G$01 incident energy (Nakada et al., 2024). 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 (Nakada et al., 2024).

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 $G$02-closed magicity, in better agreement with data than Gogny-D1S and comparable to shell-model results with $G$03EFT interaction (Nakada et al., 18 Aug 2025). 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 $G$04 in $G$05-shell nuclei (Yoshinaga et al., 8 Dec 2025).

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$06-matrix, that the effective mass is a bit low at about $G$07, and that the framework contains no explicit full $G$08-force term, because $G$09 effects are folded into effective density-dependent two-body terms (Nakada, 2019). 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 (Davesne et al., 2016). The same work nevertheless found that the central part of M3Y globally reproduces Brueckner-Hartree-Fock results in all four $G$10 channels up to and beyond saturation density, and that special combinations of partial waves can be used to constrain the tensor parameters (Davesne et al., 2016).

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 (Nakada et al., 2016). 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.