In our SHF+BCS model, we have used the energy density functional with the Skyrme interaction SLy4 and the pairing force strengths defined as in Ref.. Additionally, for the deformed SHE with , the pairing strengths have been scaled to reproduce the experimental neutron (0.696 MeV) and proton (0.803 MeV) pairing gaps in Fm. For the superheavy isotones with =184, whose experimental nuclear masses are unknown, the pairing strengths have been scaled to reproduce the pairing gaps of the finite-range droplet model (FRDM).
For all nuclei considered, a self-consistent total binding energy ( ) is computed with a quadratic constraint on the mass quadrupole moment . Our study covers the prolate deformations b (barns) with a step of 10b for prolate deformed SHE and the oblate/prolate deformations b in the case of spherical SHE. To fix the position of the nucleus center of mass, an additional constraint on the mass dipole moment, , is assumed.