np-Pair Correlations in the Isovector Pairing Model †

: A diagonalization scheme for the shell model mean-field plus isovector pairing Hamiltonian in the O(5) tensor product basis of the quasi-spin SU Λ (2) ⊗ SU I (2) chain is proposed. The advantage of the diagonalization scheme lies in the fact that not only can the isospin-conserved, charge-independent isovector pairing interaction be analyzed, but also the isospin symmetry breaking cases. More importantly, the number operator of the np-pairs can be realized in this neutron and proton quasi-spin basis, with which the np-pair occupation number and its fluctuation at the J = 0 + ground state of the model can be evaluated. As examples of the application, binding energies and low-lying J = 0 + excited states of the even–even and odd–odd N ∼ Z ds -shell nuclei are fit in the model with the charge-independent approximation, from which the neutron–proton pairing contribution to the binding energy in the ds -shell nuclei is estimated. It is observed that the decrease in the double binding-energy difference for the odd–odd nuclei is mainly due to the symmetry energy and Wigner energy contribution to the binding energy that alter the pairing staggering patten. The np-pair amplitudes in the np-pair stripping or picking-up process of these N = Z nuclei are also calculated.


Introduction
It is evident from both theoretical and experimental analysis of available data that, besides neutron-neutron (nn) and proton-proton (pp) pairing, neutron-proton (np) pairing is equally of importance in N∼Z nuclei [1][2][3][4][5][6][7][8][9]. Though isoscalar T = 0 np-pairing in nuclei is of importance in the high-energy regime [10,11], isovector T = 1 np-pairing seems dominating in thelow-energy regime [9,12], where a shell model mean-field plus isovector pairing may provide a simple and clear description of the np-pairing correlations [8,13,14]. Exact solutions of the mean-field plus charge-independent O(5) isovector pairing is available [15,16], in which the same valence proton and neutron mean-fields and the same isovector paring strengths among pp-, nn-and np-pairs are assumed. For nuclei away from the N = Z line, not only the valence neutron and proton mean-fields, but also the isovector pairing strengths among pp-, nn-and np-pairs may differ, for which ananalytic solution has not been reported.
Shell model calculations with effective interactions focusing on the neutron-proton pairing correlations have also been carried out. For example, the pair correlation was investigated by means of the shell model Monte Ca method performed with the modified Kuo-Brown interaction and the pairing plus quadrupole-quadrupole interaction within the fp-shell space [17,18]. The direct diagonalization of the KB3 interaction within the fp-shell showed that the isovector pairing interaction strength seems 2-3 times stronger than the isoscalar one when the total isospin is small, as seems to apply in this case [19,20]. Shell model calculations based on effective interactions with respect to the isovector and isoscalar pairing were also performed within sdfp and f 5 pg 9 subspaces [7,21]. Systematic analysis of N ≈ Z nuclei in various model spaces within the extended pairing plus quadrupole-quadrupole (EPQQ) Hamiltonian has been studied extensively [22][23][24]. Very recently, a distinct quartet strcuture was proposed and applied to describe isovector and isoscalar pairing correlations [25][26][27]. The isovector and isoscalar pairing in N = Z nuclei was also studied through an analysis of the shell-model wavefunctions with effective interactions as reported in [28,29]. Though the agreement of the shell model results with experimental data suggests that the isovector and isoscalar pairing interactions are realistic, the actual interaction strengths are subject to considerable uncertainty due to the fact that the competition of the isovector and isoscalar pairing, deformation and other correlations lead to a very complex picture.
On the other hand, besides calculations based onan M-scheme or a J-scheme with various algorithms [30][31][32][33][34], as shown in the pioneering work of Elliott [35][36][37] and the further work carried out by Moshinsky and many others [38], group theoretical or algebraic descriptions of the full shell model are now feasible [39,40], in which the symmetry adapted bases used are equivalent to the shell model Fock states up to a unitary transformation. It is also well known [13,14] that the isovector pairing Hamiltonian can be built by using generators of the quasispin group O (i) (5) (i = 1, · · · , p), where p is the number of orbits considered in the model space. Thus, the isovector pairing Hamiltonian can be diagonalized within a given irreducible representation of ⊗ p i=1 O (i) (5), of which the basis is simply called the O(5) tensor product basis.
Though the algebraic scheme presented in this article is specifically designed for the nuclear isovector pairing problem, it can also be extended for the O(5) nonlinear σ model that unifies antiferromagnetism and d-wave superconductivity [41]. It is obvious that the scheme in the SU(2) case, where only like-nucleon pairs are considered, can also be applied to Heisenberg spin interaction systems. Therefore, the scheme outlined in this work can also be extended to study spontaneous symmetry breaking processes and the appearance of Nambu-Goldstone modes [42] due to spontaneous symmetry breaking in quantum many-body systems [43,44].
This article is an extension of our previous work [45] and a recent Letter [46], in which neutron-proton pairing in even-even and odd-odd ds-shell nuclei was analyzed by a shell model mean-field plus isovector pairing model reailzed in the O(5) tensor product basis adapted to the quasi-spin SU Λ (2) ⊗ SU I (2) chain. In Section 2, the group theoretical classification of the shell model Fock states, with the isospin related O(5) group as a subgroup in a given j-orbit, is outlined. In Section 3, the relevant canonical and non-canonical basis of the O(5) group are briefly reviewed. In Section 4, based on the results shown in Section 3, the isovector pairing Hamiltonian is diagonalized in the O(5) tensor product basis adapted to the quasi-spin SU Λ (2) ⊗ SU I (2) chain within which some quantities related to the np-pairing corrections can be analyzed. As an example, this diagonalization scheme for the O(5) charge-independent isovector pairing model in a description of even-even and odd-odd N∼Z ds-shell nuclei are presented in Section 5. A short summary is presented in Section 6.

Fock States in a Given j-Orbit with the Group Theoretical Classification
In the shell model, let {a † j m j , t m t , a j m j , t m t } be a set of the (valence) nucleon creation and annihilation operators in the j-orbit, where m j is the quantum number of the angular momentum projection, t = 1/2 is the quantum number of the isospin, and m t = 1/2 or −1/2 is the quantum number of the isospin projection, respectively. It is well known that the total number of many-particle product (Fock) states, Γ, provided by {Ξ † φ n |0 = ∏ n m j m t a † j m j ,t m t |0 }, where |0 is the (valence) nucleon vacuum state, and φ n , up to the permutations among n creation operators, stands for the n unequal sub-indices involved, is given by which is due to the fact that the maximal number of creation operators involved in the nonzero many-particle product states is 4j + 2 restricted by Pauli exclusion. It is obvious [47] that the set of operators {Q φ n ,φ n = Ξ † φ n Ξ φ n , 1 ≤ n, n ≤ 4j + 2} generate the unitary group U(2 4j+2 ). The set of the many-particle product (Fock) where the irreducible representation (irrep) ( 1 2 , · · · , 1 2 , There are several important subgroup chains useful to provide various complete basis vectors of the irreps of O(8j + 4), among which the following chain is used to label the complete basis vectors with n j valence nucleons in the j-orbit, where T is the quantum number of the total isospin, J is the quantum number of the total angular momentum and N (j) =n j /2 − Ω j with Ω j = j + 1/2, in whichn j is the number operator of valence nucleons in the j-orbit. The generators of the O(5) group in this case are (J = 0, T = 1) pair creation operators A † µ (j), pair annihilation operators A µ (j), the number operator of valence nucleons in the j-orbitn j and isospin operators T µ (j), with corresponding to m t = 1/2 or −1/2, of which the commutation relations among the above O(5) generators were explicitly shown in [15]. The generators of Sp(2j + 1) are given by with k = 1, 3, · · · , 2j, and ρ = k, k − 1, · · · , −k for a given k, whereã jm j ,tm t = (−) j+m j a j −m j ,tm t , and A j B j (k) ρ stands for the angular momentum coupling with j ⊗ j ↓ k. It should be stated that the generators of Sp(2j + 1) shown in (7) are commutative with the O(5) generators shown in (4)- (6).
Thus, for a fixed number of valence particles n j , the labels of the O(8j + 4) irrep are redundant; the complete basis vectors of (3) may be denoted as where in (8), v 1 and v 2 , being positive integers or positive half-integers simultaneously, are used to label a possible irrep of O(5) with v 1 ≥ v 2 ≥ 0, which are related to the O(5) seniority number of nucleons v and the reduced isospin t with v 1 = Ω j − v/2 and v 2 = t. v and t indicate that there are v nucleons coupled to the isospin t, which are free from the angular momentum J = 0 and T = 1 pairs. v and t also label the corresponding irrep of Sp(2j + 1) simultaneously represented by a two-column Young diagram 2 v , 1 2t with v + 2t boxes in the first column and v boxes in the second column. β and α in (8) are multiplicity labels for given T and J needed in the reduction O(5) ↓ O T (3) ⊗ O N (2) and Sp(2j + 1) ↓ SU J (2), respectively. For the O(5) seniority zero case corresponding to v = t = 0 and J = 0 discussed in this work, the quantum numbers of Sp(2j + 1) and SU J (2) are thus neglected. In this case, for a given number of valence nucleons n j , the basis vectors (8) can be constructed by using n j (J = 0, T = 1) pair creation operators A + µ (j) coupled to isospin T as shown in [48].

Diagonalizaing the Mean-Field Plus Isovector Pairing Hamiltonian
In this section, the spherical mean-field plus isovector pairing model is diagonalized in the O(5) tensor product basis. In general, the model Hamiltonian with p orbits may be written as [49] where are collective pairing operators, π,i and ν,i are valence proton and neutron single-particle energies in the i-th orbit, G π > 0, G ν > 0 and G πν > 0 are pp-, nn-and np-pairing interaction strength, respectively.
The Hamiltonian (19) is digonalized in the subspace of tensor product ⊗ p i=1 O (i) (5) basis when p j-orbits of the shell model are considered, in which each copy of the O(5) irreducible representation (irrep) is adapted to the O(5)⊃SU Λ (2) ⊗ SU I (2)⊃U Λ (1) ⊗ U I (1) chain. Though the procedure for seniority nonzero cases is the same, in this article, only seniority-zero configuration with total angular momentum J = 0 constructed from the tensor product of p copies of the O(5) irrep (Ω i , 0) is considered, in which only equal proton and neutron quasispin I i = Λ i in the i-th orbit is allowed according to (17). Moreover, though the proton and neutron single-particle energies with π,i = ν,i can be considered in the same way, only π,i = ν,i = i (i = 1, · · · , p) case, which is a good approximation for N∼Z nuclei, is considered in the following for simplicity, with which the Hamiltonian (19) diagonalized in the seniority-zero configuration is suitable to describe low-lying J = 0 + states of both even-even and odd-odd N∼Z nuclei. Eigenstates of (19) within the seniorityzero J = 0 + subspace are denoted as where the eigenstate |ξ; n, M T with total number of valence nucleons n = ∑ p i=1 n i and total isospin projection according to the relations shown in (11), C ξ; Λ 1 ,··· , Λ p n 1 m T (1),··· , n p m T (p) is the corresponding expansion coefficient, and ξ labels the ξ-th eigenstate with the same n and M T . The number of nppairs in the i-th orbit for a given Λ i can be expressed in terms of the neutron (proton) quasi-spin as for i = 1, 2, · · · , p. Matrix elements of each terms involved in (19) labeling scheme needed in the diagonalization can be evaluated according to the results shown in the previous section, of which the analytical expressions are provided in Appendix A.
It is obvious that this diagonalization scheme is equivalent to the M T -scheme realized The results produced from this scheme have been checked against the exact solution of the model for up to three pairs shown in [15,50], which shows that the results produced from this scheme are exactly the same as those obtained from the formalism provided in [15,50].
It can be verified that the ground-state isospin projection of the model for the seniorityzero case is always M T = 0 when the number of valence nucleon pairs k = n/2 is even, while it is always M T = 1 when k = n/2 is odd. The total isospin T is a good quantum number of the system only when x = 1. In the single-j case, The possible total isospin T of the system when x = 1 is given by One can check that the dimension of the subspace spanned by the above pair states is exactly equal to the dimension of the O(5) irrep (Ω, 0) with Hence, when G π = G ν = G πν , the lowest isospin is T = 1 when k is odd, and T = 0 when k is even. Let G π = G ν = G and G πν = x.  (19) for n = 12 and n = 10 valence particles, respectively, in the ds-shell with j 1 = 1/2, j 2 = 3/2 and j 3 = 5/2 orbits, for which the experimentally deduced singleparticle energies above the 16

Model Applications to Even-Even and Odd-Odd ds-Shell Nuclei
As an application of this diagonalization scheme, some low-lying J = 0 + level energies of even-even and odd-odd A = 18-28 nuclei up to the half-filling in the ds-shell outside the 16 O core are fitted by the Hamiltonian (19) in the charge-independent approximation with G π = G ν = G πν = G, of which some results have been reported in our recent paper [46]. In order to fit binding energies of these nuclei, in addition to the mean-field plus isovector pairing, the Coulomb energy and the symmetry energy with the isospin-dependent part of the Wigner energy contribution to the binding are considered with the expression of the model Hamiltonian, the same as that used in [46,50]: whereĤ 0 is given by (19) with π,i = ν,i = i , which is the mean-field plus isovector pairing Hamiltonian (19) in the charge-independent form, BE( 16 O) = 127.619 MeV is the binding energy of the 16 O core taken as the experimental value, (n) is the average binding energy per valence nucleon in the ds-shell with j 1 = 1/2, j 2 = 3/2 and j 3 = 5/2 orbits, of which the number of valence nucleons dependent form is determined from a best fit to binding energies of all ds-shell nuclei considered.
is the Coulomb energy [51] and is the parameter of the symmetry energy and the isospin-dependent part of the Wigner energy contribution, of which the first term is taken to be the empirical global symmetry energy paramter provided in [51], while δα(A) is adjusted according to the experimental binding energy of the nucleus with a given mass number A needed to account for local deviation from the first term when the Hamiltonian (26) is used; the experimentally deduced single-particle energies above the 16 O core with 1 = 1s 1/2 = −3.27 MeV, 2 = 0p 3/2 = 0.94 MeV, 3 = 0d 5/2 = −4.14 MeV [50] are used for the mean-field, and in order to achieve a better fit for low-lying J = 0 + level energies, the overall isovector pairing strength is taken as G = 1 MeV for all the nuclei considered. The best fit yields of which the first constant is very close to the value of the average binding energy per valence nucleon with avg = −2.301 MeV used in [50]; the contribution from the second term to the binding related to the two-body interaction becomes smaller because relatively larger pairing strength is used in the present calculation, while the third term is related to the three-body interaction as further correction. The parameter δα(A) obtained from the fitting procedure is provided in Table 1. Table 1. The value of δα(A) (in MeV) of (28) obtained by fitting to the binding energies and some low-lying J = 0 + level energies of even-even and odd-odd A = 18-28 nuclei in this model. Since the binding energies and a few low-lying J π = 0 + level energies of even-even and odd-odd A = 18-28 nuclei were fit together, deviations remain between the fitted values and experimental binding energies shown in Table 2 with a root mean square deviation σ BE = 0.32 MeV, except 22 F and 22 Al, for which J π = 0 + level energies are not available experimentally. Table 3 shows the lowest experimentally known J = 0 + level energies (in MeV) of these even-even and odd-odd ds-shell nuclei fit by (26) with the same model parameters as used in fitting the binding energies, in which the corresponding shell model results (SM) obtained by using the KSHELL code [34] with the USD (W) interaction [52] are also provided for comparison. The root mean square deviation of the fit values to these excited J = 0 + level energies is σ level = 1.20 MeV, while the average deviation of the excited Exp appears to be φ = 17.4%, where the sum runs over all the excited level energies of these nuclei. In addition, when the ground state of the nucleus is not a J π = 0 + state, which cannot be determined from present calculation for J = 0 + states only, the eigen-energy of (26) is given by where E ex (ξ, T, J = 0) > 0 is the excitation energy of the ξ-th excited state with isospin T and J = 0. The theoretical value of BE(Z, N) is adjusted to reproduce a reasonable value of the excitation energy E ex (ξ, T, J = 0). Due to the Coulomb energy contribution and the freedom in adjusting the binding energy with a reasonable value of the excitation energy in this case, there is about a few hundreds of keV energy difference in these excitation energies of mirror nuclei with J = 0 ground state, as shown in Table 3. Panel (a) of Figure 2 shows the double binding-energy difference defined as [54] calculated by using the binding energies of both the related even-even and odd-odd nuclei Included in the fitting procedure, which shows that the experimental data are well-fit by the model Hamiltonian (26). Moreover, it is clearly shown in panel (a) of Figure 2 that δV pn becomes comparatively smaller for odd-odd N = Z nuclei with A = 22 and 26 in this case. Since the one-and two-body interaction dominating average binding energy term and the Coulomb energy term of (26) only contribute a Z and N independent constat to δV pn , the symmetry energy and the isospin-dependent part of the Wigner energy term seem to be the main source that alters the usual pairing gap staggering pattern, which is consistent to the claim made in [6] that the double-binding energy difference (31) actually reveals the evidence for the Wigner energy contribution to the binding, where the usual staggering pattern of the pairing gaps disappears [50]. Alternatively, instead of BE(Z, N), we calculate the double pairing energy difference defined as where in which |ξ = 1, n, TM T is the lowest eigenstate of the model, with eitherĤ P = G ∑ ρ A + ρ A ρ orĤ P = G A + 0 A 0 , of which the first one is the total pairing energy contribution, while the second one is the np-pairing energy contribution to the binding. By substituting BE(Z, N) used in (31) with which removes the symmetry energy and the isospin-dependent part of the Wigner energy contribution to the binding energy, the resultant δV pn obtained from (31) should be close to the double-pairing energy difference (32). And indeed, as shown in panel (b) of Figure 2, the value of δV pn is very close to one of the δE values calculated with the total pairing energy contribution and that with the np-pairing energy contribution. Most noticeably, the staggering pattern appears, and the actual np-pairing energy in the odd-odd N = Z nuclei turns to be comparatively strong. Table 4 shows actual nn-, pp-and np-pairing contribution at the ground state or the lowest eigenstate of (26) for these nuclei defined by and the percentage of the np-pairing energy contribution to the binding η np = E pp ). It can be observed that E (1) nn in the N = Z + 2 nuclei is the same as E (1) pp in the Z = N + 2 mirror nuclei, while E (1) pp in the N = Z nuclei due to the chargeindependent isovector pairing is adopted. However, E pp in odd-odd N = Z nuclei, which shows that the np-pairing energy contribution to the binding is the largest in odd-odd N = Z nuclei.
Since the total number of valence nucleons n and the total isospin projection M T are good quantum numbers of the system, the number of valence protons and that of valence neutrons are certainly fixed in each nucleus with n π = (n + 2M T )/2 and n ν = (n − 2M T )/2, respectively. As shown in Section 3, the number of np-pairs in the i-th orbit can be defined as q i = Ω i − 2Λ i within the seniority-zero configuration, where the neutron (proton) quasi-spin Λ i is a good quantum number in the O(5) tensor product basis adapted to the SU Λ (2) ⊗ SU I (2) ⊃ U Λ (1) ⊗ U I (1) chain. Therefore, the average number of the np-pairs in the lowest J = 0 + state of the model can be defined as Thus, the average number of nn-pairs and that of pp-pairs are given by k nn = (n ν − k np )/2, k pp = (n π − k np )/2.
k np , k nn , and k pp values for each nucleus at the lowest J = 0 + state are shown in Table 4. Since the number of np-pairs is not a conserved quantity, its fluctuation in the lowest J = 0 + state of these nuclei defined as is also provided. It can be observed from Table 4 that the k np value is a definite integer for nuclei with less than or equal to one valence neutron or proton, for which the k np value is also easily countable, while (38) must be used for evaluating k np for nuclei with more valence neutrons and protons. It is obvious that the k np value is indeed relatively large in the odd-odd N = Z nuclei, which is consistent with the larger np-pairing energy contribution to the binding shown in Table 4, while the average number of the np-pairs k np in the even-even nuclei is considerably small with very large fluctuation. For example, k np = 0.711 with ∆k np = 1.342 in 24 Mg, while k np = 0.585 with ∆k np = 1.066 in 28 Si. The ∆k np value in these even-even N = Z nuclei is almost two times of the corresponding average value. Though the np-pair occupation number defined as ζ np = k np /k in the even-even N = Z nuclei is small, the np-pairing energy contribution is still comparable to the nnor pp-pairing energy contribution. Using the data shown in pp /k pp in 24 Mg and 28 Si, respectively. As shown in Figure 3, the distribution of either the percentage of the np-pairing energy η np or the np-pair occupation number ζ np is symmetric with respect to the N = Z line owing to the charge-independent approximation with the clear even-odd staggering pattern along the lines parallel to the N = Z line on the N-Z plane, where the nuclei with fixed M T are connected by each line. In short, the np-pairing is more favored in odd-odd N = Z nuclei as described by the charge-independent isovector pairing model. Table 3. A few lowest J = 0 + level energies (in MeV) of the 22 even-even and odd-odd ds-shell nuclei fit by (26) (Th), where T ξ denotes the ξ-th excited levels with isospin T, the label g denotes the ground state, the experimental data (Exp) are taken from [53], '-' denotes the corresponding level is not observed experimentally, and the shell model results (SM) are obtained by using the KSHELL code [34] with the USD (W) interaction [52] and the parameters of (26) are the same as those used in fitting the binding energies.

O Exp
4.46 5.07 5.04 0 + (T ξ = 2 1 ) 6.52 6.80 6.52 0 + (T ξ = 0 2 ) 6.73 5.      Finally, np-pair stripping reactions, such as (α, d) or ( 3 He, p), or the corresponding picking-up process, should be sensitive tests for np-pairing correlations in nuclei, for which the np-pair amplitude determined by the matrix elements of the np-pair operator in target and product states are of importance. Here, we only calculate the np-pair amplitude contributed from the q-th orbit defined as with M T = 0 for N = Z nuclei, where |ξ i , n; T i M T and |ξ f , n + 2; T f M T are eigenstate of the isovector pairing model for target nucleus with n valence nucleons and that of product nuclei with n + 2 valence nucleons of total angular momentum J = 0, respectively, and the expression is consistent with that shown in (4.15) of [55]. The matrix elements of A † 0 (j q ) in the O(5) tensor product basis used in the expansion of (20) are given by (Ω 1 , 0) · · · (Ω q , 0) · · · (Ω p , 0) Λ 1 · · · Λ q · · · Λ p n 1 , m T (1) · · · n q , m T (q) · · · n p , m T (p) (Ω 1 , 0) · · · (Ω q , 0) · · · (Ω p , 0) Λ 1 · · · Λ q · · · Λ p n 1 , m T (1) · · · n q , m T (q) · · · n p , m T (p) which is used in calculating the np-pair amplitude contributed from the q-th orbit shown in (39). It is obvious that the selection rule for (39) is ∆T = +1, −1. The np-pair amplitudes in the ds-shell with |∆T| = 1 are shown in Table 5. As estimated in [55], in which the single-particle energies of the mean-field were taken to be degenerate, and the isospin was not considered explicitly; the amplitude for the reactions among the J = 0 + ground states or the lowest J = 0 + excited states only depends on j withB(j) = (j 2 + 2j + 3/4)/(2π), from which one achievesB(1/2) = 0.318,B(3/2) = 0.955 andB(5/2) = 1.910. TheB(5/2) value is indeed very close to the values of (39) for the reactions among the lowest J = 0 + states shown in Table 5 because the j = 5/2 orbit is the lowest in energy. Nevertheless, the amplitude for the j = 3/2 orbit shown in Table 5 are one or two orders of magnitude smaller thanB(3/2), mainly due to the fact that there is a 0.87 MeV gap between the j = 3/2 and the j = 5/2 orbit. Since the j = 1/2 orbit is the highest in energy, the corresponding amplitudes shown in Table 5 are also systematically smaller thanB(1/2).  20

Summary
In this work, a diagonalization scheme for the shell model mean-field plus isovector pairing Hamiltonian in the O(5) tensor product basis is adapted to accommocate quasispin, which means the scheme is equivalent to a M T -scheme realized in the O(5) ⊃ SU Λ (2) ⊗ SU I (2) ⊃ U Λ (1) ⊗ U I (1) basis. Additionally, the scheme conserves a chargeindependent isovector pairing interaction while accommodating isospin symmetry breaking and provides for a number operator that counts the effective number of np-pairs that can be realized in a neutron-proton quasi-spin basis.
To illustrate the vaue of the theory, the scheme was used to determine binding energies and low-lying J = 0 + excited states of even-even and odd-odd N∼Z ds-shell nuclei by fitting to the charge-independent isovector pairing interaction, within which the np-, nnand pp-pairing energy contributions to the binding in the ds-shell nuclei were estimated. The results show that the np-pairing contribution to the binding energies of the odd-odd N = Z nuclei is systematically larger than that in the even-even nuclei.
It is also shown that the decrease in the double binding-energy difference for the odd-odd nuclei is mainly due to the symmetry energy and Wigner energy contributions that alter the pairing staggering pattern. In particular, the average number of the np-pairs in the J = 0 + ground state or the lowest J = 0 + excited state of the even-even and odd-odd ds-shell nuclei were evaluated. The resuts serve to show that the average number of the np-pairs k np in the even-even N = Z nuclei is considerably smaller with large fluctuation in comparison to that in the odd-odd N = Z nuclei, which leads to the conclusion that the isovector np-pairing is more favored in odd-odd N = Z nuclei. np-pair amplitudes of given single-particle orbits of the model useful in evaluating he neutron-proton transfer reaction rates among these N = Z nuclei are also calculated.
It should be stated that the computation of the isovector np-pair number is demonstrated for the even-even and odd-odd N∼Z ds-shell nuclei described by the isovector pairing model restricted within the O(5) seniority-zero subspace only, where the isoscalar np-pairs are not involved. In order to reveal the actual np-pair contents in these N∼Z nuclei, other O(5) seniority-nonzero configurations must be considered, for which an alternative O(8) model [56][57][58] should be more convenient. Nevertheless, as has been shown in our recent work on the O(8) model [59], not only the binding energies and the low-lying J π = 0 + level energies shown in Tables 2 and 3, but also the isovector pairing energy contributions to the binding energies provided in Table 4 are the same as those calculated from the O(8) model, where the isoscalar np-pairs are also involved. Therefore, the conclusion of the present work on the isovector pairing energy contribution to the binding energies of these ds-shell nuclei is still valid even in the presence of isoscalar np-pairs.
Further applications of this scheme to the model with more j-orbits or in other major shells are straightforward, for which detailed analysis will be a part of our future work. Extensions of the algebraic scheme outlined may also be applicable to study the spontaneous symmetry breaking process in the O(5) nonlinear σ model [41] and Heisenberg spin interaction systems [43].

Conflicts of Interest:
The authors have no conflict of interest.