In order to illustrate the influence of tensor densities on single-particle and global nuclear properties, in Fig. 2 are shown the radial components, = , of the neutron vector SO densities , calculated for the nickel isotopes between =28 and 50. Calculations have been performed for the Skyrme SLy4 interaction,[19] by using the Hartree-Fock-Bogolyubov (HFB) method with the spherical symmetry assumed.[20]
One can see that the SO densities are mostly positive and peaked near the surface. At =28, the SO density is large and its major part comes from the occupied f orbital (SUS system). By adding neutrons in the shell above =28, this part is gradually cancelled by an increasing in magnitude, negative contribution from the SO partner f. At the same time contributions from the p and p orbitals appear. When both pairs of the SO partners are occupied around =38, and when the g orbital is still empty (SS system), the SO density is rather small. Beyond =40, it increases again until the g orbital becomes fully occupied at =50. Note the shift of the SO densities to larger distances, which occurs at the point of the switch-over between the dominating 1 and 1 contributions.
A similar pattern of varying SO densities is valid for all shells. For SS systems, one obtains small SO densities, while for SUS systems, the SO densities are large. Therefore, the SO densities are small at magic shells =2, 8, and 20, and large at magic shells =28, 50, 82, and 126.
Since the partners are always occupied first, the SO densities are mostly positive. Note that the derivatives of particle densities are mostly negative, and also peaked at the surface; therefore, for positive coupling constants, the SO and tensor forces split the SO partners in opposite directions, cf. Eq. (4).