Density functional theory (DFT) in its original formulation [1,2] and its extended versions [3,4,5,6] has became a universal approach to compute the ground-state (g.s.) and excited [7,8] configurations of many-electron systems held together by an external one-body potential in condensed-matter, atomic, and molecular physics. The DFT strategy has also been used in the area of nuclear structure. The nuclear DFT, a natural extension of the self-consistent mean-field (MF) theory [9,10] is a tool of choice for computations of g.s. properties and low-lying excitations of medium-mass and heavy nuclei.
There are a number of differences between the electronic DFT and nuclear DFT. First of all, because of the absence of external potential, nuclei are self-bound systems, and this creates conceptual problems due to the distinction between intrinsic (symmetry-broken) and laboratory-system densities. This problem can be treated within the extended DFT framework of Levy-Lieb [3,4,5,6,11], see Refs. [12,13,14,15]. The second difference is related to the effective interaction. A complicated form of the effective nucleon-nucleon interaction, as well as complexity of the nuclear matter saturation-mechanism (pertaining to three-body forces), give rise to the presence of higher-order terms in the nuclear energy density functional (EDF). Hence, the EDF cannot be constructed in a model-independent ab initio way based solely on bulk properties of infinite-medium augmented by exchange and gradient corrections, as it is usually done for electronic systems. The current best nuclear EDFs are still phenomenological in nature, either inspired by a MF averaging of the effective interaction such as the zero-range Skyrme  force, or based on systematic expansion involving symmetry-constrained gradient terms . By construction, the nuclear EDF is characterized by a certain number of free parameters which are fitted directly to empirical data.
The third major difference between electronic and nucleonic DFT, addressed in this study, is the presence of two types of fermions and the charge independence of the nuclear interactions giving rise to the isospin symmetry. Since the isospin symmetry is violated essentially by the Coulomb interaction, which is much weaker than the strong interaction between nucleons, many effects associated with the isospin breaking in nuclei can be treated in a perturbative way, making the formalism of isotopic spin a very powerful concept in nuclear structure and reactions [18,19].
Up to electromagnetic effects, the isospin symmetry should be conserved by elementary excitations of nuclei. This is not the case for the elementary excitations modes of the nuclear Hartree-Fock (HF) or Kohn-Sham theory, i.e., proton and/or neutron particle-hole (p-h) excitations. The independent-particle wave function manifestly breaks the isospin symmetry in the g.s. configurations of odd-odd = nuclei and in all other but isoscalar excited states. Two prominent examples are discussed in this work: superdeformed (SD) bands in a doubly magic nucleus Ni  and terminating states in = 40 nuclei .
The paper is organized as follows. We begin in Sec. 2 with general discussion of g.s. configurations of odd-odd = nuclei and p-h excitations in even-even (e-e) nuclei, so as to qualitatively introduce problems faced by MF-based theories around = and to provide motivation for the necessity of isospin restoration. Section 3 contains a detailed presentation of the isospin-projected DFT approach. In particular, we demonstrate that this approach is essentially free from problems related to uncompensated poles plaguing projection techniques ; hence, it requires no further regularization . Applications of the isospin-projected DFT approach to the SD bands in Ni and terminating states in = isotopes of Ca, Sc, Ti, V, and Cr, are presented in Secs. 4.1 and 4.2, respectively. Finally, Sec. 5 contains the conclusions of this work.