Isospin symmetry breaking in atomic nuclei

The importance of the isospin symmetry and its breaking in elucidating the properties of atomic nuclei is reviewed. The quark mass splitting and the electromagnetic origin of the isospin symmetry breaking (ISB) for nuclear many-body problem is discussed. The experimental data on isobaric analogue states cannot be described only with the Coulomb interaction, and ISB terms in the nucleon-nucleon interaction are needed to discern the observed properties. In the present work, the ISB terms are explicitly considered in nuclear energy density functional and spherical shell model approaches, and a detailed investigation of the analogue states and other properties of nuclei is performed. It is observed that isospin mixing is largest for the $N=Z$ system in the density functional approach.


Introduction
The word "symmetry" is a 16th century Latin derivative from the Greek words for "syn-" (together) and "metron" (measure).Symmetries have played a pivotal role in the development of the physical laws and breaking of these symmetries gives rise to many phenomena occurring in nature.In the development of three monumental theories of physics, the Newtonian mechanics, Einstein's relativity and the Maxwell's electromagnetism, symmetries played a prominent role.The concepts of Galilean and Lorentz invariance principles, equivalence of inertial and gravitational masses and the duality between electric and magnetic fields are the building blocks of these theories.In quantum mechanics, the symmetries are of paramount importance as they lead to conserved quantum numbers which allow the characterization of the quantum mechanical states.The fundamental translational and rotational symmetries in space-time lead to energy, momentum and angular-momentum quantum numbers [1].These symmetries have classical manifestations, but there is another class of symmetries that have no classical analogies and occur only in the quantum mechanical realm.These include exchange symmetry and the symmetries that are broken spontaneously.The former leads to the Pauli exclusion principle for fermions, and the latter is an important class of symmetries associated with phase transitions giving rise to the Nambu-Goldstone boson [2,3].
The main focus of the present work is to review, and also to discuss some new results on isospin symmetry and its breaking in nuclear physics.The isospin symmetry is an approximate symmetry and is broken explicitly by Coulomb and isospin symmetry breaking (ISB) terms in nucleon-nucleon interaction.In this work, we shall investigate the implications of both of these on the nuclear properties.Further, isospin symmetry is also broken spontaneously due to the mean-field approximation which shall be also discussed in this work.The review article is organised in the following manner: Basic ideas regarding the isospin symmetry are presented in Sec. 2. The quantum chromodynamics (QCD) origin of the ISB is discussed in Sec. 3. The classification of charge-dependent interactions between nucleons are presented in Sec. 4. Sec. 5 discusses the empirical data, manifesting the breakdown of the isospin symmetry in hadronic systems.In Sec. 6, the important consequences of the isospin symmetry in nuclear physics, i.e., isobaric mass multiplet equation (IMME) and Nolen-Schiffer anomaly are presented.Sec.7 presents the detailed investigation of the ISB using the density functional approach.The spherical shell model study of the ISB is discussed in Sec. 8. Summary and concluding remarks are included in Sec. 9.

Basic concepts of isospin symmetry
It was in the year 1932 when the neutron was discovered [4] by Chadwick and Heisenberg realised [5] that neutron and proton could be treated as two states of the same particle called "nucleon", considering that the two particles have similar masses and the same intrinsic spin.The two states of the nucleon are distinguished by the values of a new quantum number that was called as "isospin".The name isospin was christened by Wigner in 1937 [6] to express the similarity between this quantity and the spin quantum number.Although the mathematical structure of the two quantities is identical, but in comparison to the spin which is an intrinsic property of a particle, the isospin symmetry is invoked.Nevertheless, the concept of isospin has played a fundamental role in the development of quark model in particle physics [7].Further, in 1959 Weinberg derived the isospin selection rule for particle decay, ∆T ≤ 2, by recognising that Coulomb interaction weakly breaks the isospin symmetry [8].

2/44
In nuclear physics, approximate models have been developed using the concept of isospin symmetry.For instance, the famous isobaric mass multiplet equation (IMME) of Wigner [9] to predict the masses of isobaric chains of isotopes has its roots in the isospin quantum number.In Sec. 6, we shall discuss in detail the importance of this mass formula in predicting the masses of members of a multiplet which are difficult to measure experimentally.The nucleon-nucleon interaction in most of the approaches is also expanded in terms of the isospin quantum number as this allows to evaluate the effects of the ISB terms [10].
In analogy to the spin quantum number of an electron with s = 1 2 , the isospin quantum number of the nucleon is defined as, t = 1 2 with t z = +1/2(−1/2) for a neutron (proton).The total isospin projection for a nuclear system, having a total "A" number of nucleons, is given by with the possible range of isospin quantum number, T ≥ |T z |.The projection quantum number is well defined for a given system, but isospin itself is not conserved in general.For conservation of the isospin quantum number, it is required that the isospin operator must commute with the Hamiltonian of the system [11], i.e., [ where the isospin operator can be written in terms of Pauli matrices with The simplest possible nucleon-nucleon interaction that satisfies the commutation relation of Eq. ( 2) is of scalar form in the isospace, i.e., where the coefficients "a" and "b" depend on space and spin operators.The isospin invariance of a system implies that its Hamiltonian commutes with the rotation operator about any axis defined in the isospace.This is commonly referred to as charge independence (CI) in the literature [11,12].The charge symmetry (CS) operator, another related quantity defined in nuclear physics, indicates a rotation by π about an axis perpendicular to the quantization axis along which the isospin projection quantum number is measured.Conventionally, z-axis is considered as the quantization axis, and y-axis is chosen as the rotational axis.The CS operator is then given by If CS is conserved then the corresponding CS operator, P CS given by Eq. 7, commutes with the Hamiltonian.The above rotation operator will interchange "z" and "−z" axis in the isospace, which means that T z of a state will change into −T z .Therefore, CS symmetry implies the equivalence of neutron-neutron (nn) and proton-proton (pp) forces.It has been demonstrated in Ref. [11] that this equality is a necessary condition, but not sufficient for the validity of CS.The symmetry also has implications for the neutron-proton system and is shown, following the derivation given in Ref. [11], by considering the commutator of T z with the Hamiltonian operator acting on a self-conjugate state |T z = 0⟩ [ P CS , T z ]|T z = 0⟩ = ( P CS T z + P CS T z )|T z = 0⟩ , where we have made use of the following relation along with the fact that P CS is a unitary operator, i.e., P CS P † CS = I, P † CS T z P CS = − T z or T z P CS = − P CS T z .
If P CS commutes with the Hamiltonian for the self-conjugate system, then we can specify the charge (eigenvalue of T z ) and a charge parity quantum number, which is the eigenvalue of P CS , simultaneously.Since the application of P 2 CS brings the system back to its original state, i.e., the eigenvalues of P CS are ±1.For a state with a well defined value of isospin, we obtain The above equation gives the eigenvalues of CS operator in terms of the isospin quantum number, implying isospin conservation for self-conjugate nuclei obeying charge symmetry.
It can be shown that CS operator P CS converts a proton into a neutron or vice-versa only when the bare masses are considered and for physical masses a perturbative correction term needs to be considered [11].The term "bare" is the eigenstate of the Hamiltonian H 0 in the absence of electromagnetic interactions and other charge breaking terms.Denoting the "bare" neutron (proton) of momentum k, and magnetic quantum number m s by |n(p)km s ), such that and where the column vectors denote the isospin part of the wave function.It is noted that momentum and spin remain unaffected due to P CS operator.Therefore, we have where the properties of the Pauli matrices have been used.Therefore, we obtain In a similar manner, we can show that In the relativistic version, P CS in the helicity representation, and is given by [11] As indicated above, for the "bare" neutron-proton system, we have [ H 0 , P CS ] = 0.The corresponding eigen value equation is given by where M 0 is the "bare" nucleon mass.In the presence of the CSB terms, which are considered as the perturbation Hamiltonian H 1 to the charge symmetric Hamiltonian, H 0 , such that the total Hamiltonian resulting the physical neutron and proton states is given by The eigen value equations of this total Hamiltonian are given by and Thus, the "bare" strong Hamiltonian has two-fold degeneracy, resulting from its charge symmetry.This symmetry is broken by the CSB terms resulting into the physical neutron and proton with masses given by M n = 939.6MeVand M p = 938.3MeV.In the above equations, the complete physical states are distinguished by sharp kets | ⟩ instead of the round one | ) of Eqs.(11) to (17) for the "bare" states.Considering H 1 as a perturbative correction, it is shown in appendix (A), how the physical states can be obtained in the lowest order of the perturbation theory.

Quantum Chromodynamics and charge symmetry breaking
The microscopic origin of the CSB is the electromagnetic interaction and the mass difference between up-and down-quarks.Therefore, the investigation of CSB in nuclear physics will provide an insight into the quark mass problem and may also give some clues on the flavor dependence of the quark mass spectrum, which is one of the most fundamental unsolved problems in physics [13].In the following, it is demonstrated how mass difference between up-and down-quarks leads to CSB mechanism.We start with the QCD Lagrangian [13] : with the covariant derivative and the gluon field tensor

5/44
In above Eqs.( 21) to (23), "q" denote the quark fields, M the quark mass matrix, g is the strong coupling constant and A µ,a are the gluon fields.λ a are the Gell-Mann matrices and f abc the structure constant of the SU (3) color Lie algebra with a, b, c = 1, 2, ..., 8.The first generation of quark masses at 2 GeV scale in the minimal subtraction (MS) scheme from (2 + 1) flavour lattice QCD with estimated quantum electrodynamics (QED) correction are given by [14] m u = 2.19 (15) MeV , m d = 4.67 (20) MeV , m s = 94(3) MeV .
In the absence of the quark masses, the QCD Lagrangian is given by Defining right-and left-handed quark fields as q R = P R q , q L = P L q , where Lagrangian of Eq. ( 21) can then be rewritten as It is evident from the above equation that left-and right-handed components of the massless quarks do not mix.In the two-flavor case, it is SU (2) R × SU (2) L symmetry, what is called as chiral symmetry, that is preserved.It is shown below that this symmetry is broken by considering the mass term In the above equation, the chiral symmetry is broken even when the two masses are considered to be same.On the other hand, the first term conserves isospin symmetry but the second term breaks it.For m u = m d , the isospin symmetry is preserved by the QCD Lagrangian.Therefore, it is evident that mass difference between up-and down-quarks leads to the ISB, apart from the electromagnetic interactions.

Classification of charge-dependent nucleon-nucleon interaction
Following the Henley and Miller classification scheme [10,11,15], nucleon-nucleon interaction is categorised in terms of the following groups: 1. Class I : Isospin or charge-independent interaction commutes with the isospin operator and can be written as a scalar operator in the isospace where "a" and "b" are isospin independent quantities.6/44 2. Class II : The interaction in this category is charge symmetric, but breaks charge independence and can be written as In the above equation, the charge symmetry operator, e − iπ 2 τy , changes sign for τ 3 (j) and again for τ 3 (i) and therefore first term is invariant under the charge symmetry transformation.However, under an arbitrary rotation in the isospace, this term is not invariant and therefore breaks the charge independence.The second term in the above equation is isoscalar.Coulomb interaction has the form of Eq. (30).

3.
Class III : This class of forces break both charge symmetry and charge independence, but are symmetric under exchange of particles in isospace and is written as where "d" is symmetric under interchange of spin coordinates.This above class of force is different for protons and neutrons, and vanishes for the neutron-proton system.However, it does not lead to isospin mixing as [ V III , T z ] = 0 and [ T 2 , T z ] = 0.The ρ 0 − ω meson mixing gives significant contribution to this class of force and also the neutron-proton mass difference in two-pion-exchange potential.

4.
Class IV : This class of forces break charge symmetry and consequently the charge independence and is generally written as The main contribution to the first term comes from photon and ρ − ω exchange terms, and for the second term, it is from the nucleon mass difference on π and ρ exchanges.It vanishes for neutron-neutron and proton-proton systems, but causes isospin-mixing in the neutron-proton system.

Empirical data on isospin symmetry breaking
In this section, we shall review the empirical evidence on the breaking of the isospin symmetry in nuclear physics.It has been discussed in the introduction section that the idea of the isospin symmetry was conceived when neutron was discovered and Heisenberg [5] realised that since mass of this new particle is similar to that of a proton, the two particles can be treated on equal footing.The mass difference between proton and neutron, a  fundamental constant in physics, is 1.29 MeV.This mass difference is a direct manifestation of the ISB and in fact is very critical for the very existence of the universe [16].In the event that this mass difference was less than that of electron mass of 0.51 MeV, the proton will capture electrons from the atomic shells through the process called as inverse β-decay with the consequence that an atom will cease to exist.In the other scenario that this difference is larger than the observed one, it will be impossible to synthesise the elements beyond hydrogen [16].

Mass differences in mirror and analogue states
The ISB terms of electromagnetic origin and u-d mass splitting give rise to the neutron-proton mass difference.In a major achievement [17], QED and QCD equations have been simultaneously solved on a lattice using powerful computers and the masses of hadrons have been determined using four non-degenerate Wilson fermion flavors.In particular, the neutron-proton mass difference has been accurately reproduced by the calculations and it has shown that it results from the competition between electromagnetic and u-d mass ISB terms.Although the u-d mass difference is a parameter in the model as this quantity is not known empirically, however, using the same u-d mass difference, the masses of a large set of particles having different quark structures have been evaluated and are found to be consistent with the measured values [17].
In Table 1, the masses of light mirror nuclei are listed.The isospin symmetry signifies that the masses of the two mirror nuclei (isospin doublets) should be identical and any difference shall be attributed to ISB terms in the Hamiltonian.The major contributor to the ISB is the Coulomb interaction as one mirror partner has more protons than the other and will have larger Coulomb energy.We shall discuss later that nuclear interaction also breaks the isospin symmetry and contributes to the mass difference.Isospin symmetry, as a 9/44 matter of fact, entails that not just masses but all the properties of the mirror nuclei should be identical [18] apart from the corrections arising from the Coulomb interaction.In Fig. (1), the energy spectrum of the mirror nuclei, 11 B and 11 C are displayed, and the figure depicts a remarkable similarity between the two spectra, a consequence of the isospin symmetry.
The above similarity is not confined to mirror nuclei only, but holds for a broader class of states, so called, the isobaric analogue states.For a given angular-momentum, a state of a nucleus with T greater than T z is called an analogue state as it has identical structure to the state of a neighbouring isobar with |T z | = T .An example of isobaric analogue states of A = 22 triplet are depicted in Fig. (2). 22Mg and 22 Ne having two-protons and two-neutrons as valence particles have T = 1 states only, whereas 22 Na having neutron-proton pair has both T = 1 and 0 states.It is noted from the figure that after shifting the spectrum of 22 Na downwards to match the lowest T = 1 (J π = 0 + ) states of the three nuclei, the other T = 1 states with J π = 2 + , 4 + and 6 + fall in place.The difference in the energies of these states or any other properties is the indication of charge independence breaking (CIB).To further illustrate the the meaning of analogue states, in Fig. (3) the binding of the A = 21 T = 3/2 quadruplet states are depicted.For the two extreme cases, the ground-states are having T = 3/2 and for the two middle cases, the first excited states belong to the T = 3/2 quadruplet.The ground-states for the latter two cases are having T = 1/2.For mirror nuclei, only proton and neutron pairs are exchanged and it is possible to test CSB only.However, for isobaric analogue states of three or more nuclei, CIB is investigated.

Scattering lengths
The 1 S 0 scattering length provides a direct evidence of CSB and CIB effects in the nuclear force [19].In Table 2, neutron-neutron, proton-proton and neutron-proton singlet scattering length parameters are listed after correcting for the electromagnetic effects.To measure CIB effects, one defines the quantity, ∆a CD , as and for CSB the quantity, ∆ CSB , is defined as The above differences in the scattering lengths translate into 2.5% and 0.6% changes in the potential [19].These differences are significant and need to be considered to describe the nuclear processes accurately.The present level of understanding is that there are primarily two mechanisms behind the measured ISB effects in nuclear force.The first is the mass splitting observed for hadrons in the isospin channel, and the second is the meson mixing and the irreducible meson-photon exchanges.This leads to a smaller kinetic energy for neutrons as compared to protons and increases the neutron-neutron scattering length by 0.25 fm in comparison to the proton-proton one.Besides this, the nucleon mass asymmetry also contributes to the meson exchange diagrams of 2π and N ∆ types that contribute to the nuclear force [11].
The up-and down-quark mass difference and the electromagnetic interaction also leads to mixing of neutral mesons with same spin and parity, but different isospin [11].The most important mixing is between ρ 0 and ω neutral mesons and it has been shown [11] that this mixing can completely explain the observed CSB in the singlet state.However, it has been remarked that only on-shell mixing matrix elements have been extracted, and the calculations are not quite accurate [11].
In an attempt to unravel the source of CSB in nuclear force, three different models have been studied [20].These models include, ρ − ω mixing, nucleon mass splitting and the phenomenological Argonne potential (V 18 ).It has been shown that all the three models give similar values for the 1 S 0 scattering length.Although considerable departures were noted in these models for ISB in 3 P J partial waves, but the contribution of these to some observable properties, like the binding energy of the mirror pair 3 H-3 He were found to be of the order of only 6%.Further studies are being planned to probe the cause of the ISB effects in nuclear force [20].

Isobaric mass multiplet equation (IMME) and Nolen-Schiffer anomaly
In this section, the IMME and the Nolen-Schiffer anomaly are discussed.

Isobaric mass multiplet equation
In this section, we shall discuss the IMME, which is one of the most important consequences of the isospin symmetry in nuclear physics [12].It was demonstrated by Wigner [9] that the mass difference of the isospin analogue states can be expressed in terms of a quadratic polynomial in T z .This was shown using the first order perturbation correction arising from the Coulomb interaction, which explicitly breaks the isospin symmetry in nuclei.The correction was evaluated using the basis functions having well defined isospin quantum number.To begin with, it may be noted that the Coulomb interaction has the following tensorial decomposition in the isospace : 11/44 The first, second and third terms in Eq. ( 35) are isoscalar, isovector and isotensor components of the Coulomb interaction, respectively.In an analogous manner, the total Hamiltonian of the nuclear many-body system is considered to have the following isospin expansion in terms of T = 0 and 1 channels [21] : where the operator separates the explicit isospin dependence from the radial and spin components of the nucleon-nucleon interaction and are given by [21] In Eq. ( 36), the T = 1 component of nucleon-nucleon interaction, involving the radial and the spin components can be written as v where the summation over "µ" includes Coulomb, Yukawa and the isoscaler terms of the nucleon-nucleon interaction, and the quantities S and V (r) represent the spin and radial dependence of the interaction terms.
In the following, we shall express the interaction in terms of neutron-proton, neutron-neutron and protonproton terms, which allows to delineate the terms in Eq. ( 38) corresponding to ISB and CSB mechanisms.As an illustration, we consider the evaluation of the following matrix element : In the above equation "C" denotes the Clebsch-Gordan (CG) coefficients [22].The other matrix elements can be similarly calculated.Thus, the proton-proton, neutron-neutron, and the parts of proton-neutron matrix elements with T = 1 can be written as where . Similarly, 12/44 and Adding and subtracting ( 40), ( 41) and ( 42), we obtain It is now evident from Eqs. ( 44) and ( 45) that v (1) and v (2) are responsible for CSB and CIB in the nucleon-nucleon interaction, respectively.The Hamiltonian given by Eq. ( 36) can be expressed in terms of charge-independent and charge-violating parts as [21] and the binding energy is given by If the structure of H ′ CV , involves only two-body interactions, we can write it as a tensor upto rank two as where k = 0, 1 and 2 correspond to H

H
(2) In the following, we shall now obtain IMME of Wigner by considering the contribution of H ′ CV in the first order perturbation theory.This contribution gives the energy splitting of the isobaric mass multiplets, i.e.,

∆BE(αT
CV + H CV + H (2) The explicit T z dependence of the energy splitting of the multiplet can be factored out using the Wigner-Eckart theorem [22] as CV ||αT ⟩ CV ||αT ⟩ , where the algebraic expressions of the Clebsch-Gordan coefficients have been used [23].The double-bars in the above equation denote the isospin reduced matrix elements that depend on the total isospin and any other quantum number necessary to define the states uniquely.The expression on the right side of Eq. ( 53) acquires the following quadratic form by rearranging the terms: The coefficient "a" depends on the isoscalar part of the interaction with a small contribution from isotensor component, coefficient "b" is related to the isovector and "c" on isotensor parts of the nucleon-nucleon interaction.These parameters are fitted to the empirical data and include the effects of the Coulomb term and also the ISB term of the nucleon-nucleon interaction.In the derivation of Eq. ( 53), it is assumed that ISB term of the nucleon-nucleon interaction has the same structure as that of the Coulomb interaction and it is an interesting problem to derive the isospin dependence of the nucleon-nucleon interaction starting from QCD inspired microscopic theories [24].It is possible to have a rough estimate of the three parameters in Eq. ( 55) by assuming a uniform charged spherical distribution for the nucleus, and the Coulomb energy is then simply given by [12] with Recent advancements in nuclear physics have led to a reevaluation of the IMME, particularly, in light of newly available data regarding nuclei with T ≥ 3/2 [12].These experimental endeavours, while promising, present formidable challenges due to the necessity for exceptionally precise mass measurements of analog nuclei that are notably deficient of neutrons.Despite decades of effort, achieving the requisite level of precision has proven to be a gradual process.
Fig. 4 displays results of some recent studies where researchers have meticulously fitted known data to the IMME, employing up to cubic term in T z to explore the potential necessity for higher-order correction terms in order to test the predictive power of the model [12].The analysis reveals interesting insight that with the exception of the isobar A = 9, the data agrees remarkably well with the quadratic IMME expression originally proposed by Wigner [9].Moreover, coefficients derived for a uniform spherical charge distribution, as depicted in Fig. (4), exhibit a notable agreement with the fitted values.
The IMME has been derived by considering tensorial representation of the Coulomb interaction.However, since the parameters of the mass formula are fitted to the data, it also includes the contributions from the ISB terms of the nuclear potential.It is, therefore, implicitly assumed that ISB terms of the nuclear potential have also the same tensorial form as that of the Coulomb interaction.In a seminal investigation [25], the author tried to deduce the ISB contribution to the effective two-body interaction using a microscopic approach by employing the bare Argonne AV18 potential.The medium effects were incorporated by solving the Bethe-Goldstone equation within the Brueckner approach.In the following, we discuss some important results from this work [25].In nuclear matter approximation, energy density with the inclusion of CSB and CIB terms is given by It is evident from above that CSB term of the nucleon-nucleon interaction gives the first order symmetry energy coefficient, whereas the CIB contributes to the even orders only.In the above equation, the density of nuclear matter is ρ = ρ n + ρ p and the parameter β = (ρ n − ρ p )/ρ denotes the isospin asymmetry.The derivation of Eq. ( 54) is based on Wigner's assumption of |αT T z ⟩ as the eigenstate of the charge-independent Hamiltonian H 0 , where α encompasses all additional quantum numbers to uniquely specify this state.All charge-violating two-body interactions, including the Coulomb interaction H C among protons and H CSB+CIB of CSB and CIB interactions, are treated as first-order perturbations.The total binding energy is expressed as where the first-order perturbation terms H C and H CSB+CIB are expressed as tensor operators of up to rank two.The matrix elements of tensor operators can be written in terms of reduced matrix elements using the Wigner-Eckart theorem, with the coefficients involving solely depending on T and T z and also on radial and spin components.However, in a nuclear medium, H CSB+CIB becomes a density-dependent effective interaction, and can no longer be represented as an irreducible tensor [25].The perturbation energy term corresponding to CSB and CIB components, i.e., ⟨αT T z | H CSB+CIB |αT T z ⟩, lacks analytic forms unlike the Coulomb interaction.In the presence of effective CSB and CIB interactions, the perturbed energy can be represented as [25] ⟨αT where the coefficients "a" contain the zeroth order symmetry contribution and I = (N − Z)/A = 2T z /A .The generalized IMME (GIMME), incorporating the perturbation correction from CIB and CSB terms, is then expressed as 15/44 where ∆ nH = 0.782 MeV is the mass difference between a neutron and a hydrogen nucleus.The above derivation disentangles the effects of effective charge-violating nuclear interactions from those occurring due to the Coulomb potential.The terms like b c and c c in Eq. ( 60) are solely induced by the Coulomb interaction and are independent of the total isospin projection T z .However, novel components such as a (CSB) sym,1 (A, T z ) and a (CIB) sym,2 (A, T z ), originating from CSB and CIB interactions in the nuclear medium involve T z dependence.This dependence of the coefficients signifies that the microscopically derived IMME cannot be expressed as a simple polynomial in T z .
The coefficients associated with the T z and T 2 z terms of the IMME, deduced from the measured masses of T = 3/2 isobaric quartets [26], are depicted in Fig. 5, compared with those obtained from the model of a non-uniformly charged sphere [27].Furthermore, Fig. 5 illustrates the contributions stemming from CSB and CIB effects in Eq. ( 60) for nuclei with T z = T as illustrative examples.The analysis reveals that the contribution attributed to the CSB effect to the coefficient of the T z term rises approximately from −80 keV to −220 keV as the mass number A increases from 17 to 53 [25].This observation aligns well with estimations for T = 1 multiplets documented in Table 5.4 of Ref. [28].Broadly speaking, the CSB effect leads to a diminution of the T z term coefficient by 2.0% − 3.1%, while the CIB effect augments the coefficient of the T 2 z term by 1.6% − 4.4%.It is noteworthy that although the energy splitting within the isobaric multiplet primarily arise from the Coulomb interaction, the corrections to the IMME clearly have CSB and CIB origins.z extracted from the experimental data [26], along with the calculated coefficients of the Coulomb contribution (plus ∆ np ) from a simple non-uniformly charged sphere (dashed curves) are shown for comparison.(The figure is adopted from Ref. [25]).

Nolen-Schiffer anomaly
The energy splitting in an isobaric mass multiplet for a given angular-momentum and parity should be equal to the difference in the Coulomb energies and be described by the Wigner's IMME.However, it has been known for more than fifty year that even after including a string of correction terms, which includes the electromagnetic spin-orbit interactions, an overall shift from a few to ten percent in the Coulomb displacement energies (CDE) still remains unresolved.This difference is known as Nolen-Schiffer anamoly [29,30].Using IMME, CDE for adjacent members of a multiplet is given by [Ref.[25]] where T z is taken for the isobar with larger proton number.With the GIMME, the above expression is modified as where ∆ N SA is arising from CSB and CIB and is given by [Ref.[25]] with T z> = T z + 1.It is evident from above that ∆ N SA has contribution from CSB and CIB apart from the usual Coulomb perturbation [25].∆ N SA accounts for about 2% to 3% of the CDE for isobaric quartets.

7
Nuclear energy density functional approach and isospin symmetry breaking In this section, we shall employ the nuclear energy density functional (NEDF) approach to investigate the isospin mixing for the nuclear many-body system.This section is subdivided into five subsections.
In subsection (7.1), a general outline of the NEDF approach is presented from historical perspective.In subsection (7.2), generalized NEDF model is presented that preserves isospin-symmetry at the DFT level and includes densities in all isospin channels.An approximate isospin projection method is discussed in subsection (7.3).In subsection (7.4), ISB terms are introduced in NEDF method and finally in subsection (7.5) numerical results of a few isobaric chains of isotopes are presented and discussed.

Timeline of the DFT in nuclear physics
DFT has become an indispensable and a valuable framework to investigate the properties of many-body quantum systems [31].It is now widely used in physics, chemistry and material science to simulate many-body correlations.The DFT approach is based on the two theorems of Hohenberg-Kohn [32].The first of these theorems states that every system has a unique local density that describes the ground-state energy of a system, and the second theorem provides a procedure on how to obtain this density.However, the major impediment in the application of the DFT approach is how to construct the energy density functional, in the first place, in terms of the local density that will lead to the exact ground-state energy of the Schrodinger wave function approach.Hohenberg-Kohn approach only establishes that there is a local density for each system that gives the same ground-state energy as that of the Schrodinger approach.The ground-state energy of a system in an external potential, v ext , is given by [33] E(ρ) = ⟨Ψ| H|Ψ⟩ , In the first line of the above equation, the energy is expressed in terms of the standard Schrodinger wave function.The second line is based on the Hohenberg-Kohn theorem, which expresses the ground-state energy completely in terms of the spatial density.v ext in the above equation is the external potential and in the case of atomic and molecular systems, it is the Coulomb interaction.The quantity F (ρ) is the universal energy density functional for a given species of particles, and the problem with the Hohenberg-Kohn theory is that it doesn't provide any procedure on how to construct it.A commonly adopted scheme proposed by Kohn and Sham [34] is to express the density in terms of the single-particle orbits, ϕ i (⃗ r) as The kinetic energy in terms of single-particle orbits is given by The functional F (ρ) corresponds to the Hamiltonian of the system with the kinetic energy contribution, evaluated in Eq. ( 66).For atomic systems, the direct component of Coulomb energy is approximated by the classical electrostatic energy and only the unknown component of the energy functional is the exchange term.
In nuclear physics, DFT has its roots in the 1970s and 80s pursuit to develop an effective interaction that will elucidate the properties of nuclei all across the periodic table [35].In this class of nucleon-nucleon forces, the interaction potential containing several terms is completely parameterized.The most successful of these effective interactions is the Skyrme force introduced in late fifties [36][37][38].It was shown that this short-range expansion of the two-body potential leads to a simple structure of the Hamiltonian density in terms of the nuclear and kinetic densities in the Hartree-Fock (HF) approach.For the three-body force, Skyrme also introduced a zero-range force and it was shown that it corresponds to a linear density-dependent two-body force in HF approximation for time even systems.The most important advantage of the Skyrme potential is that HF equations can be derived directly without explicitly calculating two-and three-body matrix elements of the potential.The other advantage of the Skyrme interaction is that analytic expressions of the parameters are obtained in the nuclear matter approximation.In the pioneering work [39], HF calculations were performed for doubly-closed shell systems and it was demonstrated that binding energies and radii are reproduced quite accurately with two different parameter sets obtained in this work.The original parameter sets of the Skyrme interaction, referred to as SI and SII, were obtained by fitting the binding energy and density of nuclear matter, and also the binding energy and radii of 16 O and 208 Pb closed shell nuclei.
In a significant development, new parameter sets of the Skyrme interaction were obtained by fitting the coupling constants to a larger body of data on closed and semi-closed nuclei [40].In this work, binding energies and charge radii of 16 O, 40 Ca, 48 Ca, 56 Ni, 90 Zr, 140 Ce and 208 Pb were used in the fitting procedure and several new sets of Skyrme parameters were obtained.These parameter sets called as SIII, SIV, SV and SVI correspond to varying density dependencies of the interaction potential.In this exhaustive work, around 200 nuclei were studied in the spherical HF+BCS approximation and it was shown that binding energies and radii could be reasonably reproduced with all four new parameter sets.However, it was shown that calculated single-particle spectra deviate in the four sets as stronger density dependence leads to a higher level density [40].
The above sets of parameters from SI to SVI, referred to as the first generation, are rooted in the nuclear Hamiltonian with linear density dependence.It was realised very early that although these sets describe the bulk properties quite well, but completely overestimate the energy of the isoscalar giant resonance [35].This mode related to the incompressibility modulus, K nm [41] and for SKIII, its value is more than 300 MeV and the empirical estimate from the isoscalar mode is about 220 MeV.It has been shown that in order to reproduce the empirical value, the density dependence of the Skyrme interaction should be between 1/6 and 1/3 and a new Skyrme interaction set SkM was developed to reproduce the isoscalar monopole resonance [42,43].The introduction of the fractional density dependence is a complete departure from the linear dependence in the first generation of Skyrme force as this dependence cannot be derived from the Hamiltonian of the system.However, it can be justified from the Hohnberg-Kohn DFT formulation as the energy density functional of a quantum many-body system can be expressed in terms of any power (integer or fractional) of the local density.The Skyrme energy density in nuclear physics ,therefore, falls within the realm of DFT formulation of quantum mechanics as soon as the fractional power of the density is introduced.
In most of the DFT formulations in nuclear physics, beginning from 1990s, the energy density functional is constructed disregarding the Hamiltonian origin of the new terms introduced in the functional.This new path in nuclear DFT was laid in the seminal work [44] where the spin-orbit interaction of the original Skyrme Hamiltonian was altered with the objective that anomaly found in the description charge radii of the Pb-isotopes could be resolved.It was noted that Skyrme HF+BCS calculations with the standard parameter set fails to reproduce the kink found in the charge radii of Pb-isotopes at N=126 [45].However, the HF+BCS analysis performed with the relativistic density function, referred to as RMF approach [46], was capable of reproducing the kink observed in the data.In the work [44], effective non-relativistic energy density was derived from the RMF approach and it was identified that spin-orbit term contained in this functional was different from that of spin-orbit interaction employed in the original Skyrme density functional.An auxiliary term in the spin-orbit part of the density functional was introduced and by refitting the parameters, the kink observed in the charge radii of Pb-isotopes was replicated in the Skyrme approach [44].
There is now a plethora of Skyrme energy density functional parameter sets and are obtained by a fitting coupling constants to ground-state properties of magic and semi-magic nuclei, energies of giant resonances, fission barrier heights and the properties of symmetric and asymmetric nuclear matter.In a comprehensive investigation [47], 240 Skyrme functionals were scrutinized against 11 benchmark empirical properties of nuclear matter ranging from low density dilute Fermi gas to a high density neutron matter.It was shown that out of 240 sets, only 16 were able to satisfy all the proposed constraints [47].
In recent years, modern optimization techniques have been developed to fit the Skyrme density functional coupling constants to a larger set of data, including those of deformed nuclei [48].In this work, masses, charge radii and odd-even mass differences of 72 nuclei ( 28 spherical and 44 deformed ) were selected in the fitting procedure.The derivative free optimization algorithm POUNDerS [49] was used and it was found to be much superior to the standard fitting algorithms.It was shown that using the newly obtained parameter set, called as UNEDF0, returns root mean square (rms) deviations of 0.45 MeV and 1.2 Mev for two-neutron separation energies and masses, respectively, of nuclei above A=80.However, the UNEDF0 functional was not designed to provide a proper description of the single-particle energies as none of the properties related to shell structure were used in the fitting procedure.In a subsequent work [50], a new Skyrme functional, referred to as UNEDF1, was developed by also including the excitation energies of fission isomers, observed in the actinides in the fitting procedure apart from the observables employed in UNEDF0 stage.
In the third stage of the development of Skyrme functional following the UNEDF pathway, a new parameterization set so called UNEDF2 [51] was obtained by adding the tensor terms and including the data on single-particle energies of doubly magic nuclei.Further, the data set used in UNEDF1 was extended in order to improve the pairing properties in heavy mass region.However, it was realized [51] that no significant improvement in the overall description of the observed properties was achieved by employing more data in the fitting protocol.The conclusion of the work [51] is reproduced here "However, after adding new data aiming to better constrain the nuclear functional, its quality has improved only marginally.These results suggest that the standard Skyrme energy density has reached its limits, and significant changes to the form of the functional are needed".
It is quite evident that Skyrme density functional needs to be enriched to achieve the quality of predictions of spectroscopic level.This statement is also applicable to Gogny [52,53] and relativistic [46] density functionals as quality of agreement with data is similar to that achieved with Skyrme functional.There are several pathways to extend the Skyrme density functional approach and quite obvious generalization is to include the higher order derivative terms beyond the traditional quadratic form [54].There are also attempts to augment the nuclear density functionals with long-range pion exchange contributions, obtained from the chiral effective field theory [55].
Further, it is known that existing NEDF have deficiencies that can be overcome in order to improve the quality of predictions.For instance, the NEDF breaks the isospin symmetry as can be easily checked that energy functionals are expressed in terms of isoscalar and only z-component of the isovector density, x-and y-components that exchange neutrons and protons are missing.The density functional is therefore not scalar in the isospin space and breaks the isospin symmetry.This symmetry is quite important for nuclei close to the N ∼ Z line and it is known that a phenomenological Wigner term needs to be added in the microscopic mass models to correctly describe the extra binding energy measured for these systems [56,57].Further, in several studies, neutron-proton pairing term have been included in mean-field studies [58][59][60][61] and the importance of the T=0 and T=1 pairing modes have been debated.In these studies, the two pairing modes are mutually exclusive, however, in some exactly solvable models, it turns out that both modes can coexist [62,63].
Nevertheless, in all the above investigations, neutron-proton coupling has been considered only in the pairing channel and is disregarded in the particle-hole channel which clearly is an inconsistent approach that breaks the isospin symmetry.In an attempt to develop a generalized DFT approach that obeys isospin symmetry in the particle-hole channel, the nuclear energy density functional has been constructed [64] that contains terms up to quadratic order in terms of local isoscalar and the three components of the isovector densities.The energy density functional has been constructed by considering all possible terms that leave the functional invariant under spatial rotations, isospin rotations and space inversion.In the following, we shall provide some details about this generalized approach and shall closely follow the Refs.[64][65][66].

Isospin invariant nuclear DFT
In the generalized isospin invariant approach, neutrons and protons are treated as two states of the same particle and the ground-state wavefunction in the Kohn-Sham (Hartree-Fock) approximation for a system of "A" nucleons is given by where the Kohn-Sham orbitals are related to the original fermion states through a unitary transformation The particle density matrix in the isospin representation is given by and the trace of the above density gives the total particle number "A".The non-local scalar and vector densities in the coordinate space can be written in terms of the above density matrix as where σ and τ are the Pauli matrices in the spin and isospin spaces, respectively.In the above equation, m = 0 corresponds to the isoscalar component, and m = 1, 2 and 3 are the three isovector components.The energy of the nuclear system in terms of Skyrme density functional and the Coulomb terms is given by The Skyrme energy density functional is given by where the first term is the kinetic energy density.The second and third terms are isoscalar and isovector parts of the Skyrme interaction, respectively and when expressed as a functional of local densities and considering only terms that are invariant under spatial and isospin rotations, these terms can be written as where details on various terms and densities can be found in the Refs.[64][65][66].The Kohn-Sham or Hartree-Fock potential is obtained through variation of the energy density functional of Eq. ( 72), i.e., where Kohn-Sham orbitals are obtained through the solution of the following mean-field equations It is known that the mean-field approximation breaks the isospin symmetry spontaneously and leads to spurious isospin mixing in the nuclear states.This is especially true for systems with neutron excess and necessitates that isospin projection of the HF states need to be performed.As a matter of fact isospin projection was performed in early eighties [67,68] to investigate the discrepancies reported in the HF approach to describe the Coulomb displacement energies and differences in proton and neutron radii.In order to estimate the isospin mixing due to Coulomb interaction, the Hamiltonian was diagonalized in the isospin projected states.
A procedure similar to the above diagonalization of the Hamiltonian has been adopted in the Skyrme DFT approach [69,70].It has been shown that in comparison to the angular-momentum and particle-number projection formalism in the DFT which are plagued with the singularity problem [71], the isospin projection is free from this divergence.There are essentially two sources of divergences in nuclear DFT, one is due to reason that nucleon vertices are not anti-symmetrized and secondly because of the application of the generalized Wick's theorem [72][73][74].These problems don't exist for nuclear density functionals derived from the Hamiltonian based approaches, for instance, in the case of first generation of the Skyrme functionals.As already discussed earlier, in all the modern nuclear DFT approaches, the functional terms are added with no linkages to the nuclear Hamiltonian.It has been shown that nuclear DFT with integer power of density dependence, a regularization approach is possible to rectify the divergence problem.However, for nuclear DFT with fractional density dependence, regularization method is not possible and divergence problem persists [72,75].
The divergence occurs for the transition density as these are multiplied by the inverse of rotated norm kernel, which has a pole structure for some rotational angles.In the case of isospin projection, the pole structure has been investigated analytically and it has been shown that for the density dependence of the standard Skyrme functional which has order less than three, the isospin projected energy density is free from spurious divergences [70].
Using the isospin projection approach, a chain of Ca-isotopes was studied with SLy4 interaction and the isospin results were compared before rediagonalization (BR) and after rediagonalization (AR) [69].It was noted that in the absence of the Coulomb potential, the two results were quite different with BR results showing no isospin mixing for the N = Z, but then increased with neutron excess.This isospin mixing is spurious and occurs due to the spontaneous ISB by the mean-field approximation.The results of AR depict no isospin mixing in the absence of Coulomb potential as expected in an accurate treatment.In the presence of the Coulomb potential, the two results again differ with AR showing quite large mixing for the N = Z system and then dropping exponentially with neutron excess.The BR results show smaller isospin mixing for the N = Z system, but then remains non-zero even for the extreme neutron excess system of 60 Ca for which the AR shows zero mixing [69].

Approximate isospin projection method
In the present work, we have followed an approximate isospin projection approach as outlined in Refs.[65,66].In this scheme, the following auxiliary Kohn-Sham equation is solved : In the above equation, three Lagrangian parameters (λ x , λ y , λ z ) have been introduced that fix the average values of the three components of the isospin, and is similar to the cranking model approach used in high-spin nuclear physics to generate the high angular momentum rotational states.It is known that the cranking approach is a first order approximation to angular-momentum projection in the Kamlah expansion [76] of the projection operator.In Eq. ( 89), the cranking term in introduced in the isospace and, therefore, it will approximately project out the isospin quantum number.In the spherical polar coordinate system, (λ x , λ y , λ z ) = (λsinθ cosϕ, λsinθ sinϕ, λcosθ) and it has been shown that results are independent of the azimuthal angle (ϕ) [65] and, therefore, all the calculations can be performed with two parameter of λ and θ.This prescription ensures that ⟨ t y ⟩ = 0 and corresponds to the invariance of the time-reversal symmetry in the calculations.The interpretation of the results becomes quite simple in the absence of the Coulomb interaction as for this case, neutron-proton mixed DFT preserves the isospin and the states will have well defined isospin quantum number.This symmetry invariance implies that a given "T " state should be degenerate for all the isobaric analogue states which have different T z = (N − Z)/2.On the two extremes of the isobaric multiplet, the lowest states will have T = T z = |(N − Z)|/2.Each single-particle Kohn-Sham energy state 22/44 will have neutron-proton and Kramers degeneracy.For the system having T z = (N − Z)/2, the lowest "Z/2" single-particle energies will be four-fold degenerate, and the remaining higher (N − Z)/2 single-particle energies will be two-fold degenerate.This system is solved using the standard DFT approach with fixed neutron and proton numbers and the λ parameters is chosen equal to λ = λ (Tz=T ) np = (λ n − λ p ) and θ is varied to obtain the solutions for analogue states of other isobars [64].For θ = 180 • , protons and neutrons are interchanged and for values of θ between 0 • and 180 • , the states are mixed.θ = 90 • has equal neutron and proton mixing in the absence of ISB terms.Although neutron-proton mixing of the states varies with θ, but energies remain constant.
The inclusion of the Coulomb interaction favours the states with larger < T z > as these states will have lower contribution to the electrostatic energy.The single-particle energies now become "θ" dependent and the level crossings occur.This leads to ping-pong divergence problem in the HF iterative solution and it is difficult to find a converged solution.To avoid this divergence, a shift in the λ value is introduced [65] in the following manner :

Isospin symmetry breaking in nuclear DFT
The ISB terms have been introduced in Skyrme density functional approach [77] to shed light on the discrepancies obtained in the evaluation of mirror displacement energies (MDEs) and triplet displacement energies (TDEs).It was observed that using three different standard Skyrme functionals of SV T , SkM * and SLy4, the MDEs for lighter mass nuclei were underestimated by about 10%, and the observed staggering in the TDEs was completely missing in the calculated values.It was shown that by introducing ISB terms of class II and III, the discrepancies in the evaluation of MDEs and TDEs could be clarified [77].These terms were introduced in the three dimensional Skyrme DFT solver HFODD with neutron-proton mixing included at the HF level [78].In the present work, we have included these terms in the HFBTHO solver [66] and in the following we shall discuss this new development.The ISB terms introduced in the neutron-proton mixed HFBTHO program are : The corresponding contributions to NEDF are Using the above definition, Eq. ( 94) become Similarly, where ρ and s are scalar and spin-vector densities, respectively.The expressions of particle and spin densities in the cylindrical coordinate system are provided in appendix (B).The cylindrical system is used in the HFBTHO solver as it employs the axial symmetry and the DFT problem reduces to two dimensions.The two coupling constants of the ISB terms have been separately fitted to the experimental values of MDEs and TDEs [79] since their contributions have been shown to be of the order of few MeVfm 3 in comparison to the isoscalar coupling constant of the Skyrme functional, t 0 , which is about 1000 MeVfm 3 .It is noted from Eq. (96) that contribution from the class II force depends on the mixed neutron-proton densities and the generalized HFBTHO code is required to be used [66].We have performed numerical calculations for a few isobaric chains of isotopes with the inclusion of the ISB terms in order to investigate the modifications on the nuclear properties.In particular, it would quite interesting to study the additional isospin mixing admitted by the ISB terms apart from the mixing with the Coulomb potential.The investigation has been undertaken for the isotopic chains of A = 48 (T = 4) and A = 78 (T = 11).These isobaric chains have already been studied in Ref. [66] and in the present work we shall explore the modifications on these already published results with the inclusion of the ISB terms.The Skyrme functional SkM * has been employed and the coupling constants t II 0 = 7 MeVfm 3 and t III 0 = −5.6MeVfm 3 have been adopted from Ref. [77].As already stated, these parameters have been fitted to reproduce the MDEs and TDEs in lighter mass region and may not be optimum one's to study other properties or other regions.Nevertheless, the purpose of the present work is to have an overall picture of the modifications on the properties and these parameters will suffice.

Results and Discussion
The Lagrangian parameters used in the isocranking model are listed in Table (3) for the studied isobaric chains.In Fig. (6), the total energy as a function of the isocranking angle, θ are displayed for the isobaric chains of A = 48 (left panel) and 78 (right panel).The calculations have been carried out for four sets : (i) without Coulomb and without ISB , (ii) without Coulomb but with ISB , (iii) with Coulomb but without ISB and (iv) with Coulomb and also with ISB terms.The total energy without Coulomb and without ISB is independent of the θ, which means that it has the same value for all the isobars in a given chain.This is expected as nuclear DFT of Eq. ( 96) is isoscalar and changing neutrons into protons will not alter the nuclear properties.There would be some spurious isospin mixing due to the mean-field approximation, but it is evident from the calculations that this is not very significant.Total energy calculated with ISB terms but without Coulomb is lower in energy before θ = 90 • and above this angle, it is higher.This means that contribution before θ = 90 • is attractive and above it has a repulsive behaviour.This can easily understood by noting that the ISB term of class II has "t z " dependence and since its coupling constant is negative, the contribution of this term will be attractive for N > Z, zero for N = Z and negative for N < Z.This is the behaviour noted for the total energy in Fig. (6) with the inclusion of the ISB terms.The contribution of ISB terms is about 2 MeV for the A = 48 isobars of 48 Ni and 48 Ca, and is about 4 and 8 MeV for the A = 78 isobars of 78 Ni and 78 Sn, respectively.The expectation value of T 2 is displayed in Fig. (7) for the two isobaric chains.In the absence of the Coulomb potential and ISB, this value should be exactly equal to 20 (132), but this value is 20.03 (132.04) for A = 48 (78) isobars.This small deviation is caused by the spurious isospin mixing of the mean-field approximation.In the presence of the ISB terms but without Coulomb, ⟨ T 2 ⟩ is slightly higher than without ISB before θ = 90 • and after this angle, it becomes lower.In the presence of the Coulomb potential, the results are different with ⟨ T 2 ⟩ showing the maximum deviation around θ = 90 • , i.e., for N ∼ Z system.The inclusion of the ISB terms, ⟨ T 2 ⟩ depicts further deviations.The expectation values of ⟨T x ⟩ and ⟨T z ⟩ are depicted in Figs.(8) and ( 9) for the isobars of A=48 and 78, respectively.It is evident from the figures that there are minor deviations with and without ISB terms.In particular for the case without Coulomb interaction, the inclusion of ISB term for θ = 90 • has no impact on ⟨J x ⟩ since for this angle, there is equal mixing of neutrons and protons.8 Spherical shell model approach to isospin symmetry breaking The well-known phenomenological interactions like Cohen-Kurath interaction [80] for p-shell and USD interactions [81][82][83] for sd-shell are usually constructed as isospin-conserving interactions.While these 27/44 interactions are widely successful in reproducing low-energy nuclear observables, the low-energy spectra obtained for different nuclei within the same isospin multiplet would be degenerate for such isospin-conversing interactions.The isospin symmetry is broken at the nuclear level due to the following three reasons : 1. Coulomb interactions among protons, 2. mass differences between proton and neutron, and 3. the charge-dependent nature of nuclear forces.
The isospin non-conserving interactions are required to study the difference in low-energy spectra between mirror pair nuclei whose proton (Z) and neutron (N ) numbers having the same mass number (A) are exchanged.The mirror energy difference (MED) between the isospin analogue state (IAS) is a measure of the ISB in atomic nuclei, and it is defined as : Here, E I are the excitation energies of analogous states with angular momentum I in a mirror pair with T z = ±T .
To discuss isospin-related properties, we have performed calculations using the no-core shell model [84] and the shell model for a few selected pairs of nuclei across the sd-shell region.Since the no-core shell model is computationally expensive for heavier nuclei, we consider only one t z = 1 mirror pair ( 20 Na -20 F) in the lower sd-shell region.Recently, t z = 1/2 mirror pair ( 21 Na -21 Ne) has been studied in Ref. [85] within the no-core shell model formalism.In order to study the effect of ISB with the shell model, we have selected four pairs of nuclei in the mid sd-shell region corresponding to t z = ± 1/2 ( 25 Al/ 25 Mg), t z = ± 1 ( 26 Si/ 26 Mg), t z = ± 3/2 ( 23 Ne/ 23 Al), and t z = ± 2 ( 24 Ne/ 24 Si).We have calculated the low-energy spectra for the above-mentioned pairs of nuclei using the newly developed USDC [86] interaction and compared them to the corresponding experimental data.We have performed ab initio no-core shell model calculations for t z = 1 mirror pair ( 20 Na -20 F) to study the effect of ISB.The low-energy spectra for these two nuclei are shown in Fig. 10.We have used two N N interactions, namely, the inside nonlocal outside Yukawa (INOY) [87] and the chiral next-to-next-next-toleading order (N3LO) [88] interactions.The INOY interaction is a phenomenological N N interaction whose short-range part (r < 1 fm) is modified by a controllable non-local part, and the long-range part (r > 3 fm) is similar to the Yukawa tail of the Argonne V18 (AV18) [89] interaction.The N3LO interaction used in this work is a chiral interaction obtained from QCD using chiral perturbation theory (χPT).The first step in the no-core shell model is to decide an optimum frequency for each interaction and it is obtained by plotting ground-state energies with harmonic oscillator (hω) frequency for different model spaces controlled by the parameter, N max .The no-core shell model results reported in this work correspond to N max = 4, and optimum frequencies 20 and 14 MeV, respectively, for INOY and N3LO interactions.From Fig. 10, it can be seen that the ordering of the ground and first excited states for both nuclei are well reproduced by INOY and N3LO interactions.However, these two states are obtained to be degenerate up to N max = 4 model space calculation.Better results for the overall spectra can be expected for higher N max calculation.Experimentally, the ( 1 + 2 ) and 1 + 3 states have large MEDs of -527 and -487 keV, respectively.The calculated MEDs using INOY (N3LO) are -3.115(-0.852) and -4.234 (-2.956)MeV, respectively, for 1 + 2 and 1 + 3 states.In Fig. 11, we have plotted the occupancies of single nucleon harmonic oscillator orbitals up to N = 2 major shell for the 1 + 2 and 1 + 3 states.From the occupancy plots, it can be seen that the large MEDs in the NCSM results are due to the difference in occupancies of π(1s 1/2 ) of 20 Na and ν(1s 1/2 ) of 20 F. Additionally, the difference in two body matrix elements (TBME) involving 1s 1/2 orbitals for proton and neutron are also responsible for the large MEDs for these two states.

Shell model results
To study the effect of ISB within the shell model formalism, we have used the isospin non-conserving USDC interaction.This interaction is fitted to 854 states of 117 nuclei across the sd-shell.The USDC interaction is constrained by a renormalized G-matrix Hamiltonian, and an analytic Coulomb part is added along with a 2.2 % increase in the proton-neutron (T = 1) TBMEs.The full USDC Hamiltonian is provided in proton-neutron formalism comprising of three parts : 1) the isospin-conserving part, 2) the isospin-breaking part, and 3) the Coulomb part.We applied this newly developed interaction to calculate MEDs for Fig. 12 shows the comparison between the calculated and experimental low-energy spectra of t z = ± 1/2 mirror pair ( 25 Al -25 Mg).The shell model results for USDC interaction are in good agreement with the experimental low-energy spectra, except the ordering of 3/2 + 2 and 7/2 + 2 are reversed for 25 Mg compared to the experimental data.Out of the nine low-lying states considered in this work, three states, namely 1/2 + 1 , 5/2 + 2 , and 3/2 + 2 show MEDs greater than -100 keV.The experimental (calculated) MEDs of these states are -133 (-115), -175 (-130), and -128 (-95) keV, respectively.The three blue circles in Fig. 16 represent these three states.Lastly, Fig. 15 shows the comparison between the calculated and experimental low-energy spectra of T z = ± 2 mirror pair ( 24 Si -24 Ne).All three excited states for this pair of nuclei show large MEDs : 2 + 1 (-110 keV), 2 + 2 (-426 keV), and 0 + 2 (-1.298MeV).The shell results for MEDS of these states are too small compared to the experimental data.The large MEDs in the 0 + 2 state can be explained by considering the Thomas-Ehrman shift [90].The two cyan pentagons in Fig. 16 represent the other two states except for 0 + 2 , which show larger than 100 keV MED.In Fig. 17, we have plotted the proton (P) and neutron (N) occupancies of 1/2 + 1 , 5/2 + 2 , and 3/2 + 2 states of 25 Al/ 25 Mg pair along with the 0 + 2 state of 24 Si/ 24 Ne mirror pair.From the figure, it can be seen that the proton occupancies of the proton-rich partner can be compared to the neutron occupancies of the neutron-rich partner and vice-versa, and the MEDs are solely due to the difference in pp, nn, and np part of the USDC interaction.

Summary and conclusions
The concept of isospin symmetry, originally introduced by Heisenberg [5] and Wigner [6], has played a pivotal role in the development of particle and nuclear physics models.It contributed to the development of eight-fold way for the classification of elementary particles.The utility of the isospin symmetry in nuclear physics was demonstrated very early by Wigner [9].In this work, the nucleon-nucleon interaction was considered to be isoscalar and the only isospin violating term was the Coulomb interaction, which was treated using the perturbation theory.It was shown that mass excess of the isobaric analogue states can be expressed as a quadratic polynomial in T Z with the coefficients depending on the isospin quantum number.This is what is referred to as the IMME and it known to describe the masses of the analogue states quite accurately.As a matter of fact, this model has been found to be quite useful to predict the nuclear masses of the analogue states for which it is difficult to perform the measurements [26].The masses of the isobaric analogue states are needed in modelling the stellar evolution.
The isospin symmetry has provided a useful framework for the classification of the nucleon-nucleon interaction.The nuclear potential is expressed in terms of isoscalar, isovector and isotensor components, what is referred to as Henley-Miller classification [10,11,15].The dominant part of the nuclear interaction is the isoscalar component and the isospin symmetry is weakly broken by isovector and isotensor components.This is evident from the empirical data available from nucleon-nucleon scattering.ISB has two manifestations : CSB and CIB.CSB breaking is a subgroup of the CIB as it involves a special rotational invariance in the isospin space.In the Henley-Miller classification scheme, class I force is scalar in isospin space, Class II force maintains charge symmetry, but breaks charge independence and class III and IV break both charge symmetry and charge independence.
The ab initio no-core shell model and shell model are two powerful approaches to study the small difference in isobaric analogue states due to ISB.While the no-core shell model results overpredict the corresponding MEDs, the shell model results with USDC interaction underpredict the experimental MEDs.Better results for the no-core shell model can be expected for higher N max calculations, and some improved results from shell model calculations can be obtained by using USDCm and USDIm interactions [86] where the Coulomb terms are further constrained.
In an attempt to describe MDEs and TDEs using the nuclear DFT approach, the ISB terms of class II and III have been introduced in the nuclear density functional [77].The conventional Skyrme DFT is unable to illustrate these properties, and, in particular, the staggering observed in TDEs is completely absent in the DFT calculations.It has been demonstrated [77] that with the inclusion of the ISB terms in the Skyrme DFT, staggering observed in the data is reproduced.The two parameters have been fitted to the masses of the lighter mass nuclei independently without tweaking the original coupling constants of the Skyrme functional.We have investigated the isobaric chains of isotopes using these newly introduced ISB terms in the energy functional.The study has been performed using the isospin cranking approach, which is an approximation to the isospin projection.It has been shown that these terms are attractive in nature for N > Z and become repulsive for Z < N .For the N = Z system, there is no contribution of the ISB term to the energy.
For an accurate treatment of ISB terms in DFT, isospin projection needs to be performed.The ISB term of class II involves neutron-proton mixed densities and the proton and neutron particle numbers are not well defined in the HF approach.Therefore, three-dimensional isospin projection is required to be carried out.In the absence of the ISB terms, one-dimensional isospin projection has been performed in Skyrme DFT approach and it has been shown that this particular projection is free from the divergence problem that particle-number and angular-momentum projection is beset with.For this special case, the singularity problem has been analyzed analytically and it would be interesting to investigate whether the three-dimensional isospin projection is also free from the divergence problem.
The coupling constants of the ISB terms have been separately fitted to the masses in the lighter mass region.Although the contribution of ISB terms is quite small, nevertheless, in an accurate treatment the parameters of the Skyrme functional need to be refitted by also involving the coupling constants of the ISB terms in the fitting protocol.Further, the ISB terms introduced in the energy functional have been adopted from the shell model perespective.It is very essential to derive the ISB terms of the nuclear potential from the ab initio methods, for instance, following the chiral effective field theory approach.

B Densities in Cylindrical Basis
In the case of the axially symmetric even-even nuclei with z-axis as symmetry axis, the third component of the total angular momentum j z is a good quantum number for the representation of a single-particle state.Here, we introduce the axial coordinate representation (rϕz) for the HFB wave functions.If Ω k is the eigenvalue of j z , associated with the single-particle state k, then we represent the quasiparticle wave functions for Ω k > 0 as where In the explicit derivation of densities and currents in the cylindrical coordinates, the following expressions are quite useful, where σ r ss ′ , σ ϕ ss ′ , σ z ss ′ are Pauli spin matrices in cylindrical polar coordinate system and σ 0 ss ′ = δ ss ′ .
3. For the Laplacian, we have The particle-hole (ph) isoscalar and isovector non-local densities are expressed in terms of the quasiparticle wave functions and isospin Pauli matrices as Using the cylindrical expansion of the quasiparticle wave function Eq. (B1), we obtain where τ 0 tt ′ = δ tt ′ and τ 1 tt ′ , τ 2 tt ′ , τ 3 tt ′ are the isospin Pauli matrices.

Vector Spin Density
The ph spin isoscalar and isovector densities are expressed as

Figure 1 .Figure 2 .
Figure 1.(Color online) Approximate charge symmetry of nuclear levels in mirror nuclei.

Figure 3 .
Figure 3. (Color online) The experimental binding energies of the A = 21, T z = ± 1 2 , ± 3 2 nuclei -both for ground states and for the T = 32 isobaric analogue states.The Coulomb displacement energy (CDE) between adjacent members of the T = 3 2 quadruplet is indicated.(Adopted from Ref.[12]).

Figure 4 .
Figure 4. (Color online) Coefficients of the IMME of the form ∆BE(T, T z ) = a + bT z + cT 2 z + dT 3 z obtained from fits to experimental data for T = 3 2 quadruplets .(a) The experimentally determined d coefficients, expected to be zero.(b) The experimentally determined b, c and d coefficients along with predictions of the simple charged sphere model.(Adopted from Ref. [12]).

Figure 5 .
Figure 5.For the T = 3/2 quartets, coefficients of (a) T z and (b) T 2z extracted from the experimental data[26], along with the calculated coefficients of the Coulomb contribution (plus ∆ np ) from a simple non-uniformly charged sphere (dashed curves) are shown for comparison.(The figure is adopted from Ref.[25]).

Figure 6 .
Figure 6.(Color online) Total HF energy of the A = 48 and 78 for the isospin states of T = 4 and 11, respectively.

Figure 9 .
Figure 9. (Color online) Expectation values ⟨ T x ⟩ and ⟨ T z ⟩ of the A = 78 for the isospin state T = 11.

Figure 10 .
Figure 10.Comparison of no-core shell model results with the experimental data for low-lying spectra of 20 Na and 20 F isotopes using INOY and N3LO interactions.

Figure 16 .
Figure16.The experimental MEDs are compared to the calculated ones using USDC interaction for a few states whose observed MEDs are more than 100 keV.

Table 1 .
The multiplet mass splitting of the lightest nuclear isodoublets.

Table 2 .
Charge independence of nuclear Hamiltonian imposes the equalities of the NN scattering lengths;
23Ne/23Al), and T z = ± 2 ( 24 Ne/ 24 Si) mirror pairs in the mid sd-shell region and the results are shown in Figs. 12 -17.Figure 12.Comparison of shell model results with the experimental data corresponding to low-lying states for t z = ± 1/2 mirror pair ( 25 Al -25 Mg).
Figure 13.Comparison of shell model results with the experimental data corresponding to low-lying states for t z = 1 mirror pair ( 26 Si -26 Mg).A comparison between the calculated and experimental spectra consisting of seven low-lying states of T z = ± 1 mirror pair ( 26 Si -26 Mg) is shown in Fig.13.The shell model results are in good agreement with experimental spectra except the ordering of 2 + 3 and 4 + 1 states are reversed for the case of26Si.Out of these seven low-lying states, five states show more than 100 keV MEDs. ) keV, respectively.A large MED of 477 keV is observed experimentally between the (4 + ) state of26Si and 4 + state of26Mg.The calculated MED for this state is 46 keV compared to the experimental MED of 477 keV.Further investigation of this state would be interesting using different charge-dependent interactions.The five red triangles in Fig.16represent those five states with MEDs more than -100 keV.Comparison of shell model results with the experimental data corresponding to low-lying states for t z = 3/2 mirror pair ( 23 Al-23 Ne).The low-energy spectra of T z = ± 3/2 mirror pair ( 23 Al -23 Ne) is shown in Fig.14.The calculated results for six low-lying states are in good agreement with the experimental data.Out of these states, 1/2 + 1 , 3/2 + 1 and 3/2 + 2 show large MEDs of -467, -347, and -432 keV, respectively.The calculated results are -146, -137, and -37 keV, respectively, corresponding to those three states.A comparison of calculated and experimental MEDs for these three states is shown as three green squares in Fig.16.Comparison of shell model results with the experimental data corresponding to low-lying states for t z = ± 2 mirror pair ( 24 Si-24 Ne).