We suggest a new Gogny-type finite-range effective interaction including a third Gaussian in the central term. Based on simple arguments valid for an arbitrary radial form factor, the three ranges are obtained in connection with physical grounds, relating them to one-boson exchange interactions. Moreover, some parameters of the longest range are fixed through the G-matrix elements of the One Pion Exchange Potential. On top of giving a fairly good description of atomic nuclei properties comparable with other existing parametrisations, the resulting interaction leads to a remarkable improvement of some infinite matter properties that are relevant for astrophysical calculations.

We present a new Gogny-type finite-range effective interaction including a third Gaussian in the central term. Based on simple arguments valid for an arbitrary radial form factor, the three ranges are fixed from physical grounds, relating them to one-boson exchange interactions. Moreover, some parameters of the longest range are fixed through the G-matrix elements of the One Pion Exchange Potential. On top of giving a fairly good description of atomic nuclei properties, the resulting interaction leads to a remarkable improvement of some infinite matter properties that are relevant for astrophysical calculations.

We describe the new version (v3.06h) of the code HFODD that solves the universal nonrelativistic nuclear DFT Hartree-Fock or Hartree-Fock-Bogolyubov problem by using the Cartesian deformed harmonic-oscillator basis. In the new version, we implemented the following new features: (i) zero-range three- and four-body central terms, (ii) zero-range three-body gradient terms, (iii) zero-range tensor terms, (iv) zero-range isospin-breaking terms, (v) finite-range higher-order regularized terms, (vi) finite-range separable terms, (vii) zero-range two-body pairing terms, (viii) multi-quasiparticle blocking, (ix) Pfaffian overlaps, (x) particle-number and parity symmetry restoration, (xi) axialization, (xii) Wigner functions, (xiii) choice of the harmonic-oscillator basis, (xiv) fixed omega partitions, (xv) consistency formula between energy and fields, and we corrected several errors in the previous versions.

We discuss advanced statistical methods to improve parameter estimation of nuclear models. In particular, using the Liquid Drop Model for nuclear binding energies, we show that the area around the global $\chi^2$ minimum can be efficiently identified using Gaussian Process Emulation. We also demonstrate how Markov-chain Monte-Carlo sampling is a valuable tool for visualising and analysing the associated multidimensional likelihood surface.

We perform a detailed analysis of the parameters of the Skyrme SN2LO interaction. By mean of a covariance matrix analysis, we have been able to provide error bars and build the complete covariance matrix. The latter will be used in future applications to properly study the propagation of errors on nuclear observables calculated with this pseudo-potential.

We discuss advanced statistical methods to improve parameter estimation of nuclear models. In particular, using the Liquid Drop Model for nuclear binding energies, we show that the area around the global χ minimum can be efficiently identified using Gaussian Process Emulation. We also demonstrate how Markov-chain Monte-Carlo sampling is a valuable tool for visualising and analysing the associated multidimensional likelihood surface.

We present recent developments obtained in the so-called N2LO extension of the usual Skyrme pseudo-potential.In particular, we discuss the isovector splitting mass in infinite nuclear matter and the pairing gaps of selected semi-magic even-even nuclei.

We present the development of the extended Skyrme N2LO pseudopotential in the case of spherical even-even nuclei calculations. The energy density functional is first presented. Then we derive the mean-field equations and discuss the numerical method used to solve the resulting fourth-order differential equation together with the behavior of the solutions at the origin. Finally, a fitting procedure for such an N2LO interaction is discussed and we provide a first parametrization. Typical ground-state observables are calculated and compared against experimental data.

By considering infinite matter constraints only, we suggest in this paper that the Gogny interaction should benefit from a third Gaussian in its central part. A statistical analysis is given to select the possible ranges which are compatible with these constraints and which minimize a $\chi^2$ function.

We discuss some infinite matter properties of two finite-range interactions widely used for nuclear structure calculations, namely Gogny and M3Y interactions. We show that some useful informations can be deduced for the central, tensor and spin-orbit terms from the spin-isospin channels and the partial wave decomposition of the symmetric nuclear matter equation of state. We show in particular that the central part of the Gogny interaction should benefit from the introduction of a third Gaussian and the tensor parameters of both interactions can be deduced from special combinations of partial waves. We also discuss the fact that the spin-orbit of the M3Y interaction is not compatible with local gauge invariance. Finally, we show that the zero-range limit of both families of interactions coincides with the specific form of the zero-range N3LO Skyrme interaction and we emphasize from this analogy the benefits of N3LO.

We perform a detailed analysis of the properties of the finite-range tensor term associated with the Gogny and M3Y effective interactions. In particular, by using a partial-wave decomposition of the equation of state of symmetric nuclear matter, we show how we can extract their tensor parameters directly from microscopic results based on bare nucleon-nucleon interactions. Furthermore, we show that the zero-range limit of both finite-range interactions has the form of the next-to-next-to-next-leading-order (N3LO) Skyrme pseudopotential, which thus constitutes a reliable approximation in the density range relevant for finite nuclei. Finally, we use Brueckner-Hartree-Fock results to fix the tensor parameters for the three effective interactions.

We discuss the contributions to the equation of state for the NlLO Skyrme pseudopotential (l = 2,3). We show that by adding fourth- and sixth-order gradient terms, it is possible to fairly reproduce the spin/isospin decomposition of an equation of state obtained from ab initio methods. Moreover, by inspecting the partial-wave decomposition of the equation of state, we show for the first time a possible way to add explicit constraints on the sign of the tensor terms of the Skyrme interaction.

The possibility of adding a D-wave term to the standard Skyrme effective interaction has been widely considered in the past. Such a term has been shown to appear in the next-to-next-to-leading order of the Skyrme pseudo-potential. The aim of the present article is to provide the necessary tools to incorporate this term in a fitting procedure: first, a mean-field equation written in spherical symmetry in order to describe spherical nuclei and second, the response function to detect unphysical instabilities. With these tools it will be possible to build a new fitting procedure to determine the coupling constants of the new functional.