## Perturbation Scheme for the E ff ective Nuclear Force

Kun Ratha K ean 1 2 , , Takashi N ishikawa 2 , and Yoritaka I wata 3

1 Japan Agency of Atomic Energy, Ibaraki, Japan

2 Institute of Innovative Research, Tokyo Institute of Technology, Tokyo, Japan

3 Kansai University, Osaka, Japan

E-mail: iwata phys@08.alumni.u-tokyo.ac.jp

(Received July 14, 2019)

Based on a systematic Hartree-Fock + BCScalculation, new perturbation scheme is proposed to modify the existing Skyrme-type e ff ective nuclear force. Much attention is paid to have a better description of fission property of heavy nuclei. In terms of incorporating the nuclear medium e ff ects 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%.

KEYWORDS: Hartree-Fock+BCS, nuclear Skyrme interaction, new perturbation scheme

## 1. 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:

$$\delta \left < \psi | \mathcal { H } - \beta Q | \psi \right > = 0,$$

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.

## 2. Doubly-constrained Skyrme perturbation scheme

Let ρ be the density. For Skyrme-type e ff ective nuclear interaction [2], the Hamiltonian density reads

$$H = t _ { 0 } \rho ^ { 2 } + \frac { t _ { 3 } } { 6 } \rho ^ { 2 + \alpha } + \cdots.$$

We focus on the competition between t 0 ρ 2 (two-body force) and t 3 ρ 2 + α (medium e ff ect) 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 e ff ects (medium e ff ect); more precisely, nuclear

Table I. Skyrme parameter sets (10 parameters profiling the Skyrme-type e ff ective nuclear force; for the detail, e.g., see [6]).

|                           | SLy4         | SkM*         | mSLy4        | mSkM*        |
|---------------------------|--------------|--------------|--------------|--------------|
| t 0 (MeV · fm 3 )         | - 2488 . 913 | - 2645 . 000 | - 2488 . 913 | - 2645 . 000 |
| t 1 (MeV · fm 5 )         | 486.818      | 410.000      | 486.818      | 410.000      |
| t 2 (MeV · fm 5 )         | - 546 . 395  | - 135 . 000  | - 546 . 395  | - 135 . 000  |
| t 3 (MeV · fm 3(1 + α ) ) | 13777.000    | 15595.000    | 14645.710    | 16004.450    |
| x 0                       | 0.834        | 0.090        | 0.834        | 0.090        |
| x 1                       | - 0 . 344    | 0.000        | - 0 . 344    | 0.000        |
| x 2                       | - 1 . 000    | 0.000        | - 1 . 000    | 0.000        |
| x 3                       | 1 . 354      | 0.000        | 1 . 354      | 0.000        |
| W 0 (MeV · fm 5 )         | 123 . 000    | 130.000      | 123 . 000    | 130.000      |
| α                         | 0.166667     | 0.166667     | 0.195967     | 0.179000     |

Table II. Binding energy per nucleon [MeV]. The calculated values are compared to the experimental value.

|        | Experiment [4]   | SLy4   |   SkM* | mSLy4   | mSkM*   |
|--------|------------------|--------|--------|---------|---------|
| 208 Pb | 7 . 87           | 7 . 87 |   7.85 | 7 . 86  | 7 . 80  |
| 120 Sn | 8 . 50           | 8 . 50 |   8.48 | 8 . 52  | 8 . 45  |
| 90 Zr  | 8 . 71           | 8 . 73 |   8.48 | 8 . 80  | 8 . 45  |
| 40 Ca  | 8 . 55           | 8 . 66 |   8.63 | 8 . 81  | 8 . 67  |
| 16 O   | 7 . 98           | 8.10   |   8.14 | 8.37    | 8.25    |

Table III. Nuclear radius (proton radius - neutron radius - total radius) [fm]. The calculated values are compared to the experimental value (proton radius).

|        | Experiment [7]   | SLy4               | SkM*               | mSLy4              | mSkM*              |
|--------|------------------|--------------------|--------------------|--------------------|--------------------|
| 208 Pb | 5 . 50           | 5.46 - 5.62 - 5.55 | 5.53 - 5.62 - 5.56 | 5.70 - 5.84 - 5.78 | 5.58 - 5.74 - 5.68 |
| 120 Sn | 4 . 65           | 4.59 - 4.73 - 4.68 | 4.58 - 4.73 - 4.67 | 4.79 - 4.92 - 4.86 | 4.68 - 4.83 - 4.77 |
| 90 Zr  | 4 . 27           | 4.23 - 4.28 - 4.26 | 4.23 - 4.28 - 4.26 | 4.40 - 4.47 - 4.43 | 4.32 - 4.37 - 4.35 |
| 40 Ca  | 3 . 48           | 3.45 - 3.38 - 3.41 | 3.46 - 3.40 - 3.43 | 3.57 - 3.50 - 3.53 | 3.53 - 3.47 - 3.50 |
| 16 O   | 2 . 70           | 2.79 - 2.72 - 2.76 | 2.77 - 2.71 - 2.74 | 2.83 - 2.76 - 2.80 | 2.84 - 2.78 - 2.81 |

density-dependent term ( t 3-term) of the Skyrme interaction. The refit protocol of Skyrme perturbation is shown as follows:

first adding a perturbation δα to fit the theoretically calculated fission barrier height to the experiment; second 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

Table IV. Nuclear matter properties calculated by [8] (for the detail, see [9]). From left side, energy per nucleon E A / , imcompressivity K , symmetry energy S , slope of the symmetry energy L are shown.

| parameter   |   E / A [MeV] |   K [MeV] |   S [MeV] |   L [MeV] |
|-------------|---------------|-----------|-----------|-----------|
| SkM*        |        -15.77 |       216 |        30 |        46 |
| SLy4        |        -15.97 |       231 |        32 |        46 |
| mSkM*       |        -15.39 |       234 |        30 |        44 |
| mSLy4       |        -15.21 |       269 |        31 |        32 |

correctly only if

$$\left | 1 - \frac { < \psi | t _ { 3 } ^ { \prime } \rho ^ { 2 + \alpha + \delta \alpha } | \psi > } { < \psi | t _ { 3 } \rho ^ { 2 + \alpha } | \psi > } \right | < 0. 0 5$$

is satisfied. This condition of permitting 5% di ff erence is rather reasonable, because almost a few % di ff erences in binding energy already exist in the original interactions. The proposed Skyrme perturbation theory makes use of the di ff erent mathematical behavior due to the power index ( α ) and the multiplication parameter ( t 3).

## 3. Results

Two Skyrme-type e ff ective nuclear interactions mSkM* and mSLy4 interactions are proposed by modifying SkM* [5] and SLy4 [6], respectively The di ff erences 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 e ff ective 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, where odd nuclei are well calculated without blocking.

## 4. Discussion

The Skyrme-type e ff ective 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%. The proposed interactions somewhat underestimate the saturation density, which will be addressed in the future research.

## References

- [1] P.-G. Reinhard, D. J. Dean, W. Nazarewicz et al. , Phys. Rev. C 60 014316 (1999).
- [2] D. Vautherin, and D. M. Brink, Phys. Rev. C 5 626 (1972).
- [3] RIPL-2: https: // www-nds.iaea.org RIPL-2 / /
- [4] NuDat 2: https: // www.nndc.bnl.gov nudat2 / /

Table V. Fission barrier height of heavy nuclei [MeV]: inner barrier (left) and outer barrier (right).

|        | Experiment (RIPL-2 [3])   | SLy4            | SkM*            | mSLy4           | mSkM*           |
|--------|---------------------------|-----------------|-----------------|-----------------|-----------------|
| 231 U  | 4 . 40 , 5 . 50           | 6 . 90 , 6 . 51 | 5 . 34 , 5 . 50 | 3 . 63 , 4 . 72 | 3 . 36 , 3 . 37 |
| 232 U  | 4 . 90 , 5 . 40           | 7 . 10 , 8 . 66 | 5 . 63 , 5 . 48 | 3 . 84 , 4 . 84 | 3 . 65 , 3 . 45 |
| 233 U  | 4 . 35 , 5 . 55           | 7 . 49 , 8 . 09 | 6 . 01 , 5 . 55 | 4 . 11 , 4 . 99 | 4 . 01 , 3 . 56 |
| 234 U  | 4 . 80 , 5 . 50           | 8 . 02 , 8 . 31 | 6 . 45 , 6 . 60 | 4 . 38 , 5 . 28 | 4 . 35 , 3 . 76 |
| 235 U  | 5 . 25 , 6 . 00           | 8 . 77 , 8 . 75 | 6 . 99 , 5 . 86 | 4 . 68 , 5 . 56 | 4 . 74 , 3 . 97 |
| 236 U  | 5 . 00 , 5 . 67           | 7 . 36 , 8 . 29 | 7 . 53 , 6 . 18 | 5 . 00 , 5 . 89 | 5 . 09 , 4 . 27 |
| 237 U  | 6 . 40 , 6 . 15           | 9 . 86 , 8 . 52 | 7 . 93 , 6 . 46 | 5 . 30 , 6 . 15 | 5 . 37 , 4 . 50 |
| 238 U  | 6 . 30 , 5 . 50           | 10 . 2 , 8 . 92 | 6 . 37 , 6 . 97 | 5 . 60 , 6 . 51 | 5 . 66 , 4 . 78 |
| 239 U  | 6 . 45 , 6 . 00           | 10 . 7 , 9 . 28 | 8 . 70 , 8 . 21 | 5 . 88 , 6 . 90 | 5 . 88 , 5 . 05 |
| 236 Np | 5 . 90 , 5 . 40           | 8 . 87 , 8 . 27 | 7 . 11 , 5 . 49 | 4 . 63 , 4 . 89 | 4 . 79 , 3 . 40 |
| 237 Np | 6 . 00 , 5 . 40           | 9 . 52 , 8 . 39 | 7 . 66 , 7 . 76 | 4 . 98 , 5 . 21 | 5 . 15 , 3 . 67 |
| 238 Np | 6 . 50 , 5 . 75           | 10 . 0 , 8 . 35 | 8 . 16 , 6 . 18 | 5 . 29 , 5 . 51 | 5 . 48 , 3 . 93 |
| 237 Pu | 5 . 10 , 5 . 15           | 8 . 92 , 7 . 69 | 7 . 17 , 4 . 50 | 4 . 55 , 4 . 17 | 4 . 77 , 2 . 74 |
| 238 Pu | 5 . 60 , 5 . 10           | 9 . 55 , 8 . 14 | 7 . 75 , 5 . 34 | 4 . 89 , 4 . 47 | 5 . 16 , 2 . 98 |
| 239 Pu | 6 . 20 , 5 . 70           | 10 . 1 , 8 . 67 | 8 . 25 , 5 . 62 | 5 . 25 , 4 . 85 | 5 . 55 , 3 . 31 |
| 240 Pu | 6 . 05 , 5 . 15           | 10 . 6 , 8 . 79 | 8 . 76 , 6 . 08 | 5 . 59 , 5 . 25 | 5 . 87 , 3 . 62 |
| 241 Pu | 6 . 15 , 5 . 50           | 11 . 2 , 8 . 81 | 9 . 19 , 6 . 34 | 5 . 90 , 5 . 64 | 6 . 20 , 3 . 96 |
| 242 Pu | 5 . 85 , 5 . 05           | 11 . 6 , 8 . 92 | 9 . 48 , 6 . 69 | 6 . 20 , 6 . 06 | 6 . 39 , 4 . 21 |
| 243 Pu | 6 . 05 , 5 . 45           | 11 . 8 , 9 . 07 | 9 . 56 , 6 . 97 | 6 . 46 , 6 . 47 | 6 . 52 , 4 . 41 |
| 244 Pu | 5 . 70 , 4 . 85           | 11 . 9 , 9 . 34 | 9 . 49 , 7 . 19 | 6 . 65 , 6 . 60 | 6 . 59 , 4 . 50 |
| 245 Pu | 5 . 85 , 5 . 25           | 11 . 7 , 10 . 2 | 9 . 37 , 7 . 52 | 6 . 73 , 6 . 43 | 6 . 53 , 4 . 16 |
| 239 Am | 6 . 00 , 5 . 40           | 9 . 51 , 7 . 50 | 7 . 74 , 4 . 75 | 4 . 80 , 3 . 74 | 5 . 15 , 2 . 29 |
| 240 Am | 6 . 10 , 6 . 00           | 10 . 1 , 8 . 00 | 8 . 29 , 5 . 05 | 5 . 19 , 4 . 15 | 5 . 59 , 2 . 66 |
| 241 Am | 6 . 00 , 5 . 35           | 10 . 7 , 8 . 66 | 8 . 87 , 5 . 43 | 5 . 56 , 4 . 53 | 5 . 96 , 2 . 99 |
| 242 Am | 6 . 32 , 5 . 78           | 11 . 3 , 8 . 54 | 9 . 35 , 5 . 89 | 5 . 90 , 4 . 97 | 6 . 27 , 3 . 29 |
| 243 Am | 6 . 40 , 5 . 05           | 11 . 9 , 8 . 57 | 9 . 74 , 6 . 26 | 6 . 22 , 5 . 35 | 6 . 52 , 3 . 59 |
| 244 Am | 6 . 25 , 5 . 90           | 12 . 1 , 8 . 74 | 9 . 81 , 8 . 52 | 6 . 49 , 5 . 75 | 6 . 68 , 3 . 78 |
| 241 Cm | 7 . 15 , 5 . 50           | 10 . 1 , 7 . 39 | 8 . 35 , 4 . 45 | 5 . 12 , 3 . 39 | 5 . 56 , 1 . 92 |
| 242 Cm | 6 . 65 , 5 . 00           | 10 . 7 , 8 . 04 | 9 . 00 , 4 . 94 | 5 . 49 , 3 . 79 | 5 . 96 , 2 . 28 |
| 243 Cm | 6 . 33 , 5 . 40           | 11 . 5 , 8 . 25 | 9 . 53 , 5 . 32 | 5 . 84 , 4 . 21 | 6 . 30 , 2 . 61 |
| 244 Cm | 6 . 18 , 5 . 10           | 12 . 0 , 8 . 21 | 9 . 91 , 9 . 74 | 6 . 18 , 4 . 59 | 6 . 54 , 2 . 87 |
| 245 Cm | 6 . 35 , 5 . 45           | 12 . 3 , 8 . 51 | 10 . 0 , 6 . 02 | 6 . 46 , 5 . 02 | 6 . 76 , 3 . 09 |
| 246 Cm | 6 . 00 , 4 . 80           | 12 . 4 , 9 . 14 | 10 . 0 , 6 . 25 | 6 . 68 , 4 . 91 | 6 . 83 , 3 . 12 |
| 247 Cm | 6 . 12 , 5 . 10           | 12 . 2 , 9 . 86 | 9 . 92 , 6 . 58 | 6 . 77 , 4 . 63 | 6 . 79 , 2 . 85 |
| 248 Cm | 5 . 80 , 4 . 80           | 11 . 8 , 10 . 7 | 9 . 68 , 6 . 02 | 6 . 77 , 4 . 23 | 6 . 72 , 2 . 52 |
| 249 Cm | 5 . 63 , 4 . 95           | 11 . 5 , 10 . 4 | 9 . 40 , 5 . 50 | 6 . 65 , 3 . 79 | 6 . 56 , 2 . 14 |

- [5] J. Bartel, P. Quentin, M. Brack, C. Guet, and H.B. Hakansson, Nucl. Phys. A 386 , 79 (1982).
- [6] E. Chabanat, P. Bonche, P. Haensel, J. Meyer, and F. Schae ff er, Nucl. Phys. A 635 , 231 (1998).
- [7] IAEA database: https: // www-nds.iaea.org radii / /
- [8] J. Stone, private communication.
- [9] M. Dutra, O. Lourenc ¸o, J. S. S´ Martins, and A. Delfino, Phys.Rev. a C 85 035201 (2012).
- [10] P. Kl¨pfel, P.-G. Reinhard, T. J. B¨rvenich, and J. A. Maruhn, Phys. Rev. u u C 79 034310 (2009).
Kun Ratha Kean Takashi Nishikawa Yoritaka Iwata
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$\%$.