content_image/344411.jpg content_image/344412.jpg
is calculated by microscopic c.m. correction 

<img src='content_image/344411.jpg'>

with A mass number and ˆ P$_{c}$$_{.}$$_{m}$$_{.}$ = ∑ A i ˆ p$_{i}$ the total momentum in the c.m. frame. It has been shown that the microscopic c.m. correction provides more reasonable and reliable results than phenomenological ones [18–20]. 

Therefore, the total energy for the nuclear system becomes 

## III. NUMERICAL DETAILS 

<img src='content_image/344412.jpg'>

In this work, a series of calculations have been performed for both the spherical and deformed nuclei. The Dirac equation for nucleons is solved in a three-dimensional harmonic oscillator basis [21]. For spherical calculations, by increasing the fermionic shells from N$_{f}$ = 20 to N$_{f}$ = 22, the binding energy, charge radius, and neutron skin thickness in $^{208}$Pb change by 0.003%, 0.007%, and 0.1% respectively. For $^{240}$U, the binding energy changes by 0.001% from the axially deformed calculation with N$_{f}$ = 16 to N$_{f}$ = 18. Therefore, a basis of 20 major oscillator shells is used in the spherical calculations and 16 shells in the axially deformed cases. The triaxial calculations are performed with N$_{f}$ = 12, which, for Nd isotopes, provides an accuracy of 0.04% for binding energies in comparison with the calculations with N$_{f}$ = 14. To achieve an accuracy of ∼ 100 keV in the description of both the fission barrier and the energy of the isomer state in $^{240}$Pu, a basis of 20 oscillator shells has been adopted in both the axial and triaxial calculations. 

In order to determine the parameters of Lagrangian density in Eq. (2) and the pairing strength in Eq. (20), a multiparameter fitting to both the binding energies and charge radii for selected spherical nuclei is performed with the Levenberg-Marquardt method [22]. As usual, the masses of neutron and proton are fixed as 939 MeV. The corresponding data [23– 25] for selected spherical nuclei used in the fitting procedure are listed in Table II and III. The empirical neutron pairing gaps for $^{122}$Sn, $^{124}$Sn, and $^{200}$Pb as well as the proton ones for $^{92}$Mo, $^{136}$Xe, and $^{144}$Sm obtained with five-point formula are also employed to constrain the pairing strengths. 

