Perturbation Scheme for the Effective Nuclear Force

Based on a systematic Hartree-Fock+BCS calculation, new perturbation scheme is proposed to modify the existing Skyrme-type effective nuclear force. Much attention is paid to have a better description of fission property of heavy nuclei. In terms of incorporating the nuclear medium effects on the fission barrier height more correctly, two typical Skyrme parameter sets (SkM* and SLy4) are modified. In conclusion, according to the comparison to experiments, the calculated fission barrier heights are improved by almost 90$\%$.


Constrained Hartree-Fock+BCS theory
For carrying out the density functional calculation, the constrained Hartree-Fock+BCS theory (CHF+BCS) is utilized to impose a constraint on the quadrupole deformation. The master equation is obtained by the variational principle: where H means the Hamiltonian operator of many nucleon systems, the quadrupole parameter β plays a role of the Lagrangian multiplier for the quadrupole constraint βQ, and the trial function ψ is taken as the Slater determinant. The BCS-type pairing interaction is included in H together with the nuclear and the Coulomb interactions. Each deformed state and the corresponding energy surface are obtained by choosing the value of β.
The calculation is performed using the SkyAX code [1] in which the quadrupole deformation is given on the three-dimensional Cartesian coordinate. In the SkyAX code, the octupole moment is optimized by adding a small octupole moment to the initial wave functions under the quadrupole constraint. Although the axial symmetry is assumed for the SkyAX calculations, it does not require anything more for the quadrupole constraint calculations. Indeed, the quadrupole-deformed nuclei can be fully described within the axial symmetric framework.

Doubly-constrained Skyrme perturbation scheme
Let ρ be the density. For Skyrme-type effective nuclear interaction [2], the Hamiltonian density reads We focus on the competition between t 0 ρ 2 (two-body force) and t 3 ρ 2+α (medium effect) to have a better description of fission and a better prediction power for heavy nuclear physics in general. The perturbation scheme optimizes the nuclear medium effects (medium effect); more precisely, nuclear density-dependent term (t 3 -term) of the Skyrme interaction. The refit protocol of Skyrme perturbation is shown as follows: • adding a perturbation δα to fit the theoretically calculated fission barrier height to the experiment; • changing the parameter t 3 to keep the quality of the original interaction in terms of reproducing the static properties; A whole t 3 -term is optimized as t ′ 3 ρ 2+α+δα to reproduce the experimental fission property, where δa and t ′ 3 are the perturbation and the modified-value for t 3 , respectively. One of the two reference quantities is the fission barrier height (to be compared to [3]), and the other is the binding energy (to be compared to [4]). The Skyrme perturbation scheme is judged to be working correctly only if is satisfied. This condition of permitting 5% difference is rather reasonable, because almost a few % differences in binding energy already exist in the original interactions. The proposed Skyrme perturbation theory makes use of the different mathematical behavior due to the power index (α) and the multiplication parameter (t 3 ).

Results
Two Skyrme-type effective nuclear interactions mSkM* and mSLy4 interactions are proposed by modifying SkM* [5] and SLy4 [6], respectively The differences between the existing theoretical values and the experimental values (5.00 MeV [3]) are 2.36 MeV for SLy4 and 2.53 MeV for SkM*. The experimental fission barrier height for 236 U (inner barrier) is utilized to modify the effective nuclear force parameters, and the binding energy of 236 U is also fitted to keep the original quality of static calculations.
The obtained interaction is shown in Table I. For the validity check, the binding energies and the radii of 16 O, 40 Ca, 90 Zr, 120 Sn, and 208 Pb are also utilized (Table II and III), as they are referred to in most of Skyrme parameter fittings. For the quality check, nuclear matter properties are compared in Table IV. The fission barrier height is systematically calculated for Uranium to Curium isotopes (Table V). Calculated isotopes are selected based on whether the comparable experimental data [3] exist or not. Possibly due to the inclusion of the pairing interaction, the quality of barrier height calculations for odd-odd, odd-even and even-odd nuclei are as good as those for even-even nuclei.

Discussion
The Skyrme-type effective nuclear forces typically overestimate the fission barrier heights, which prevents us to have a precise prediction/control of nuclear fission processes. The proposed interaction shows a significant improvement of the calculated fission barrier height not only for Uranium isotopes but also for the other isotopes. As an average, for both two interactions, the description of fission barrier height is improved by almost 90%.