All calculations presented in this work were obtained by using isospin projection implemented in the DFT solver HFODD [57,58]. We used the harmonic-oscillator basis of =12 shells. We considered SIII [41], SLy4 [55], and SkP [56] Skyrme functionals, as well as for their Landau and tensor versions denoted by subscripts and , respectively, see Ref. [59] for details. Pairing was neglected.