Processing math: 100%

Number of states with given spin I and isospin T for n fermions in a j orbit

  • In this study, we investigate formulas of the number of states with a given total spin I and isospin T for n nucleons in a single-j shell denoted by DIT(j,n). Talmi's recursion formulas for the number of states with a given z-axis projection of total spin are generalized by further considering the isospin couplings and are applied to derive explicit formulas of DIT(j,n).
  • Determining the number of states with a given spin I for identical nucleons in a single-j shell [denoted by DI(j,n)] is a common practice in nuclear structure theory and atomic physics. In the nuclear shell model, DI(j,n) is usually tabulated for relatively small j and n. In cases with larger j or n, DI(j,n) might be obtained by subtracting the number of states with total angular momentum projection M=I+1 from that with M=I [1], where M equals m1+m2++mn, and mi [i=1,2,,n] is the projection of j on the z-axis for the i-th nucleon. Other recipes to extract DI(j,n) include Racha's formulas in terms of the seniority scheme [2] and generating function method studied extensively by Katriel et al. [3] and Sunko and collaborators [46].

    The first explicit and compact formula of DI(j,n) was obtained specifically for I=0 and n=4 by Ginocchio and Haxton while studying the fractional quantum Hall effect [7]. In Ref. [8], DI(j,n) was constructed empirically for n=3 and 4, and for a few I values with n=5. The formulas for n=3 was soon proved by Talmi, who introduced his recursive formulas [9]. These formulas were later derived based on the reduction rule from SU(4) to SO(3) group [10]; the formula for n=4 was derived by the reduction rule from SU(5) to the SO(3) group [11], with a demonstration that DI(j,n) can be actually derived based on the reduction rule from SU(n+1) to SO(3) group. The Ginocchio-Haxton formula of D0(j,4) was also revisited by Zamick and Escuderos [12]. An explicit recursion formula from n1 particles to n particles was obtained in Ref. [13] by introducing "pseudo" particles which allow fermions take integer spins. In Refs. [14, 15], the number of states, denoted by DI(j,n), was applied to derive sum rules of six-j and nine-j symbols, some of which were also revisited in Ref. [ 16]. In the last decade, Pain and collaborators extensively studied the odd-even staggering of DI(j,n) [17, 18], and compact formulas of DI(j,n) for n=3,4, and 5 [19]. The enumeration of the number of states with given spin was also extended to boson systems in Ref. [20]. The study of DI(j,n) motivated a number of studies related to a single-j shell; here, we mention Refs. [2134] without providing further details for completeness.

    Given that the nuclear shell model Hamiltonian respects the isospin symmetry, it is natural to generalize the enumeration of DI(j,n) to the number of states with given spin I and isospin T, denoted by DIT(j,n), for nucleons in a single-j shell. In Ref. [35], Zamick and Escuderos found a few simple relations between DIT(j,4) with T=0 and DIT(j,4) with T=2. In Ref. [36], compact and explicit formulas of DIT(j,n) for n=3 and 4 were derived in terms of sum rules of six-j and nine-j symbols involving the expression of DI(j,n) given in Ref. [11].

    Similar to the enumeration of DI(j,n), DIT(j,n) can be obtained in terms of the number of states with given spin projection M and isospin projection MT. This number is denoted by NMMT(j,n). For n fermions in a single-j shell,

    DIT(j,n)=[NM=I,MT=T(j,n)NM=I+1,MT=T(j,n)][NM=I,MT=T+1(j,n)NM=I+1,MT=T+1(j,n)]=NM=I,MT=T(j,n)+NM=I+1,MT=T+1(j,n)NM=I+1,MT=T(j,n)NM=I,MT=T+1(j,n).

    (1)

    In this study, we investigate DIT(j,n) in terms of NMMT(j,n), which is obtained by generalizing Talmi's recursion formula of NM(j,n), in which isospin is not considered. Compact formulas of DIT(j,n) are derived for n=3 as an exemplification. The rest of this paper is organized as follows. In Sec. II, we generalize Talmi's formulas for NM(j,n) by further considering the isospin symmetry, and present recursion formulas of NMMT. In Sec. III, we apply our generalized recursion formulas of NMMT and extract NMMT(n,j) for n=3. In Sec. IV, DIT(j,3) is derived by using Eq. (1). In Sec. V, we summarize our study. In Appendix A, we provide an alternative mathematical proof of our generalized Talmi's recursion formula by using the generating function method. In Appendix B, we tabulate DIT(j,n) with n=3 and 4, j=5/2,7/2, and 9/2.

    We use the convention that the z-axis projection of the isospin for a neutron and a proton equals 1/2 and 1/2, respectively. For n nucleons in a single-j shell, we denote the projections of total spin I and total isospin T, respectively, as M and MT. With this convention, the neutron number equals n/2+MT, and the proton number equals n/2MT. Let us denote the z-axis spin projections of neutrons and protons by using (mν1,mν2,...,mν2j+1) and (mπ1,mπ2,...,mπ2j+1), respectively, with the convention that j=mν1>mν2>...>mν2j+1=j for neutrons and j=mπ1>mπ2>...>mπ2j+1=j for protons. We define nνi=1 (or nπi=1) for i=1,...,2j+1 if the orbit of mνi (or mπi) is filled; otherwise, it equals zero.

    According to the above conventions, we have

    M=nν1mν1+nν2mν2+...+nν2j+1mν2j+1+nπ1mπ1+nπ2mπ2+...+nπ2j+1mπ2j+1,

    and the maximum of I equals

    Imax=Mmax=(j+(j1)+...+(jn/2MT+1))+(j+(j1)+...+(jn/2+MT+1))=n(2j+1)/2n2/4M2T.

    (2)

    We denote the number of states with given M and MT for n nucleons in a single-j shell by using NMMT(j,n). Clearly, we can divide the values of (nν1,nν2j+1,nπ1,nπ2j+1) into sixteen cases, as tabulated in Table 1. Here, we exemplify them by using ( nν1=1,nν2j+1=1,nπ1=1,nπ2j+1=1). In this case, there are (n4) nucleons distributed in orbits for which the absolute values of mi values equal or are below (j1). Consequently, the number of states with given M and MT, i.e., NMMT(j,n), equals the number of states of the same M and MT, but with (n4) nucleons in a (j1) shell, namely, NMMT(j1,n4). According to this classification, we are able to obtain NMMT(j,n) by using the values of NMMT(j,n), where either j is smaller than j, or n is smaller than n, or both (j,n) are smaller than (j,n). In other words, from Table 1 we have

    Table 1

    Table 1.  Number of states corresponding to the available values of nν1,nν2j+1,nπ1,nπ2j+1.
    nν1nν2j+1nπ1nπ2j+1Number of states
    1111NM,MT(j1,n4)
    1110NMj,MT1/2(j1,n3)
    1101NM+j,MT1/2(j1,n3)
    1100NM,MT1(j1,n2)
    1011NMj,MT+1/2(j1,n3)
    1010NM2j,MT(j1,n2)
    1001NM,MT(j1,n2)
    1000NMj,MT1/2(j1,n1)
    0111NM+j,MT+1/2(j1,n3)
    0111NM,MT(j1,n2)
    0101NM+2j,MT(j1,n2)
    0100NM+j,MT1/2(j1,n1)
    0011NM,MT+1(j1,n2)
    0010NMj,MT+1/2(j1,n1)
    0001NM+j,MT+1/2(j1,n1)
    0000NM,MT(j1,n)
    DownLoad: CSV
    Show Table

    NMMT(j,n)={nν1,nν2j+1,nπ1,nπ2j+1}NMj(nν1nν2j+1+nπ1nπ2j+1),MT12(nν1+nν2j+1nπ1nπ2j+1)(j1,nnν1nν2j+1nπ1nπ2j+1).

    (3)

    Here, the summation over {nν1,nν2j+1,nπ1,nπ2j+1} refers to the sixteen cases listed in Table 1. Note that the above formula holds also for n=1,2,3,4, with the convention that NMMT(j,n)=0 if n<0 and NMMT(j,0)=δM,0δMT,0. Note also that Eq. (3) can be proved alternatively based on the generating-function technique described in Ref. [3], as presented in Appendix A.

    The value of NMMT(j,n) with n=1 or 2 is trivially determined as follows. For MT=±1/2 and |M|j, NM,MT=±12(j,1)=1; otherwise, it equals zero. For n=2 with MT=±1, the highest value of M equals Mmax=2j1 and for MT=0, Mmax=2j. For MT=±1 and M2j1,

    NM,MT=±1(j,2)=2j+1|M|2,

    (4)

    where ""and"" means taking the largest integer without exceeding the value inside; for MT=0 and M2j,

    NM,MT=0(j,2)=2j+1|M|.

    (5)

    Otherwise, NMMT(j,2)=0. With these results, we can construct explicit formulas of NMMT(j,3) by using Eq. (3).

    We first discuss NMMT(j,3) with the requirement that MT=12 and j<MMmax. We recursively apply Eq. (3) to this term and obtain the value of NMMT(j,n):

    NMMT(j,n)=NMMT(j1,n)+{nν1,nν2j+1,nπ1,nπ2j+1}NMj(nν1nν2j+1+nπ1nπ2j+1),MT12(nν1+nν2j+1nπ1nπ2j+1)(j1,nnν1nν2j+1nπ1nπ2j+1)=NMMT(jl,n)+l1i=0{nν1,nν2j+1,nπ1,nπ2j+1}NM(ji)(nν1nν2j+1+nπ1nπ2j+1),MT12(nν1+nν2j+1nπ1nπ2j+1)(ji1,nnν1nν2j+1nπ1nπ2j+1),

    (6)

    where denotes summations that exclude the case {nν1,nν2j+1,nπ1,nπ2j+1}={0,0,0,0} and decouple the term NM,MT(j1,n) [the last term in the summation of the sixteen terms of Table 1]; l is fixed with the requirement that NMMT(jl,n)=0 while NMMT(jl+1,n)0, which means that

    M>n(2j+1)2n24M2Tln,Mn(2j+1)2n24M2T(l1)n,

    or equivalently,

    l=2j+12n4M2T+Mn+1.

    (7)

    For n=3 and MT=1/2, the value of l becomes

    l=jM+13+1.

    (8)

    According to Eq. (6), with the conditions NMMT(jl,n)=0 and NMMT(jl+1,n)0, we have

    NM,MT=12(j,3)=l1i=0{nν1,nν2j+1,nπ1,nπ2j+1}NM(ji)(nν1nν2j+1+nπ1nπ2j+1),1212(nν1+nν2j+1nπ1nπ2j+1)(ji1,3nν1nν2j+1nπ1nπ2j+1).

    (9)

    On the right hand side of the above formula, there are many terms that break the requirement of Eq. (2) for n2. Let us exemplify this by using two cases: {nν1,nν2j+1,nπ1,nπ2j+1}={0,0,0,1} and {1,1,1,0}. For the first case, {nν1,nν2j+1,nπ1,nπ2j+1}={0,0,0,1}, the corresponding contribution to NM,MT=12(j,3) on the right hand of Eq. (10) is NM+(ji),MT=1(ji1,2). Given that M>j here, the value of M+(ji) is larger than the maximum Mmax of two nucleons in a single ji1 shell, i.e., Mmax=2j2i2. Similarly, for {nν1,nν2j+1,nπ1,nπ2j+1}={1,1,1,0}, the corresponding contribution on the right hand side is NM(ji),MT=0(ji1,0). Given that M(ji)>0, NM(ji),MT=0(ji1,0)=δM(ji),0=0. As a result, we obtain

    NM,MT=12(j,3)=l1i=0(NM(ji),MT=0(ji1,2)+NM(ji),MT=1(ji1,2)+NM2(ji),MT=12(ji1,1))=l1i=0(3j3i1M+3j3i1M2+1)=(3j1M)l3l(l1)2+l1i=03j3i1M2+l.

    (10)

    The summation of 3j3i1M2 is tedious but straightforward. One has

    l1={2kifwhen3jM=6k+1,6k+2,6k+3,2k+1ifwhen3jM=6k+4,6k+5,6k+6.

    When 3jM=6k+1, the summation becomes

    l1i=03j3i1M2=2ki=0(3k2i+i2)=2ki=0(3k2i)+2k1i=0i+k=3k2+k.

    This summation for the case 3jM=6k+2,3,4,5,6 can be obtained similarly, as follows:

    l1i=03j3i1M2={3k2+kif3jM=6k+1,3k2+2kif3jM=6k+2,3k2+3k+1if3jM=6k+3,3k2+4k+1if3jM=6k+4,3k2+5k+2if3jM=6k+5,3k2+6k+3if3jM=6k+6.

    (11)

    Substituting these results into Eq. (10), we obtain

    NM,MT=12(j,3)={9k2+6k+1if3jM=6k+1,9k2+9k+2if3jM=6k+2,9k2+12k+4if3jM=6k+3,9k2+15k+6if3jM=6k+4,9k2+18k+9if3jM=6k+5,9k2+21k+12if3jM=6k+6=14(3jM)2+12(3jM)+η(3jM),

    (12)

    where η(3jM)=(0,14) if mod(3jM,2)=(0,1), respectively.

    Let us now address NM,MT=12(j,3) with Tz=1/2 and 0<Mj. In this case, Eq. (3) is applied recursively for jM+1 times. As a result, we have

    NMMT(j,n)=NMMT(j1,n)+{nν1,nν2j+1,nπ1,nπ2j+1}NMj(nν1nν2j+1+nπ1nπ2j+1),MT12(nν1+nν2j+1nπ1nπ2j+1)(j1,nnν1nν2j+1nπ1nπ2j+1)=NMMT(M1,n)+jMi=0{nν1,nν2j+1,nπ1,nπ2j+1}NM(ji)(nν1nν2j+1+nπ1nπ2j+1),MT12(nν1+nν2j+1nπ1nπ2j+1)(ji1,nnν1nν2j+1nπ1nπ2j+1),

    (13)

    where, again, denotes summation among the sixteen sets of {nν1,nν2j+1,nπ1,nπ2j+1} excluding {nν1,nν2j+1,nπ1,nπ2j+1}={0,0,0,0}. For MT=12 and n=3, the above formula is reduced to

    NM,MT=12(j,3)=1+NM,MT=12(M1,3)+jMi=0(NM(ji),MT=0(ji1,2)+NM+(ji),MT=0(ji1,2)+NM(ji),MT=1(ji1,2)+NM+(ji),MT=1(ji1,2)+NM,MT=12(ji1,1)+NM2(ji),MT=12(ji1,1)+NM,MT=12(ji1,1)+NM,MT=12(ji1,1)+NM+2(ji),MT=12(ji1,1)+NM,MT=32(ji1,1))=1+NM,MT=12(M1,3)+jMi=0(ji+M1)+jM2i=0(jiM1)

    +jMi=0ji+M12+jM3i=0jiM12+3jM1i=01,

    (14)

    where the first term " 1" is given by NM=0,MT=0(ji1,0)=1 with {nν1,nν2j+1,nπ1,nπ2j+1,i}={1,1,1,0,jM}; the second term NM,MT=12(M1,3) is given in Eq. (12) with j=M1, i.e.,

    NM,MT=12(M1,3)=M22M+34.

    (15)

    The summations in Eq. (14) equal zero when the upper limit of the given summation index is larger than its lower limit. On the right hand side, the summations of ji+M12 and jiM12 are classified into two cases, jM=2k and 2k+1. For jM=2k, we obtain

    jMi=0ji+M12=2ki=02ki+2M12=2ki=0(ki+2M12+i2)=2ki=0(ki+2M12)+2k1i=0i+k=(jM+1)(M12)+(jM2)2.

    The summation with jM=2k+1 is obtained similarly. The results of these two cases are unified as follows:

    jMi=0ji+M12=(jM+14)(M12)+(jM2)2+{0ifmod(jM,2)=0,14ifmod(jM,2)=1.

    (16)

    Following a similar procedure, we obtain

    jM3i=0jiM12=(jM2)2jM2+{0ifmod(jM,2)=0,+14ifmod(jM,2)=1.

    (17)

    Substituting Eqs. (15)–(17) into Eq. (14), we obtain NM,MT=12(j,3) with 0<Mj in a compact form:

    NM,MT=12(j,3)=32j2+j+14M22.

    (18)

    The results for NMMT(j,n) are reduced to NM(j,n) for identical particles, and Eq. (3) is reduced to Talmi's recursion formula. For this case, NM,MT=32(j,3) was obtained in Ref. [19], which we cite for completeness. For j<MMmax, according to Eqs. (2.11a)–(2.11b) in Ref. [19],

    NM,MT=32(j,3)=112(3jM)2+α(3jM),

    (19)

    where

    α(3jM)=0,112,13,14,13,112,

    for mod(3jM,6)=0,1,2,3,4,5, respectively. For 0<Mj, according to Eqs. (2.19a)-(2.19b) in Ref. [19],

    NM,MT=32(j,3)=12j216M2+β(M12),

    (20)

    where

    β(3jM)=112,14,112,

    for mod(M12,3)=0,1,2, respectively.

    Provided the results of NMMT(j,3), i.e., Eqs. (12), (18), (19), and (20), we can readily substitute them into Eq. (1) to obtain the number of states with given I and T, i.e., DIT(j,3).

    We first consider the case T=1/2. When I>j and T=1/2, we obtain

    DI,T=12(j,3)=NI,12(j,3)+NI+1,32(j,3)NI+1,12(j,3)NI,32(j,3).

    (21)

    According to Eqs. (12) and (19), the above formula is reduced to

    DI,T=12(j,3)=13I+j+13+η(3jI)η(3jI1)+α(3jI1)α(3jI)=13I+j+{0ifmod(3jI,3)=0,23ifmod(3jI,3)=1,13ifmod(3jI,3)=2=3jI+23.

    (22)

    For I=j and T=1/2,

    DI,T=12(j,3)=NI,12(j,3)+NI+1,32(j,3)NI+1,12(j,3)NI,32(j,3).

    (23)

    By using Eqs. (12), (18), (19), and (20), we have

    DI,T=12(j,3)=23j+712+α(2j1)η(2j1)β(j12)=23j+{23ifmod(2j,3)=1,0ifmod(2j,3)=3,13ifmod(2j,3)=5=2j+23.

    (24)

    Similarly, for 0<I<j and T=1/2,

    DI,T=12(j,3)=NI,12(j,3)+NI+1,32(j,3)NI+1,12(j,3)NI,32(j,3).

    (25)

    By using Eqs. (18) and (20), the above formula is reduced to

    DI,T=12(j,3)=23J+13+β(I+12)β(I12)={23I+23ifmod(2I,3)=1,23Iifmod(2I,3)=3,23I+13ifmod(2I,3)=5=2I+23.

    (26)

    Eqs. (24)–(26) can be rewritten as follows. For Ij and T=1/2,

    DI,T=12(j,3)=3jI+23,

    (27)

    and for 0<Ij and T=1/2,

    DI,T=12(j,3)=2I+23.

    (28)

    These results [Eqs. (27)–(28) in this paper] are consistent with Eqs. (25)–(26) in Ref. [ 36], which were obtained in terms of sum rules of six-j symbols, except that the formulas derived in this paper are expressed in a more transparent and understandable form.

    The T=3/2 case is much simpler than the T=1/2 case. One uses Talmi's recursion formulas for NM=I,32(j,n) and obtains DI,T=32(j,3) straightforwardly. The resulting formulas are consistent with those given in Ref. [ 8] and Ref. [19]. We show them next for completeness. According to Eq. (1) in Ref. [ 8], for 0<Ij and T=3/2,

    DI,T=32(j,3)=2I+36;

    (29)

    and according to Eq. (2) in Ref. [ 8], for Ij and T=3/2,

    DI,T=32(j,3)=3j3I6+δI.

    (30)

    where

    δI={0if mod(3j3I,6)=1,1otherwise.

    In principle, one can follow the above procedure as explained in Eqs. (1)–(3) for larger n values. However, this procedure becomes formidably tedious with lengthy formulas and tables for larger n values. Instead of those formulas, we calculate DIT(j,n) for n=3 and 4 and j=5/2, 7/2, and 9/2 and tabulate them in Table B1 and Table B2 in Appendix B.

    To summarize, in this paper we generalize Talmi's recursion formula of the number of states with given spin I by further consideration of the isospin symmetry. This generalization is also proved alternatively based on the generating function method. We also exemplify the generalized Talmi's recursion formulas obtained in this paper by applying our generalized Talmi's recursion formulas in constructing the number of states with given spin and isospin of three nucleons in a single-j shell.

    In this Appendix, we present an alternative proof of the recursion formula presented in Eq. (3) based on the generating-function method.

    A generating function is a polynomial function for which the expansion coefficients are related to the number of states with a given requirement. We define the generating function corresponding to NMMT(j,n), similar to that in Ref. [17],

    fj(x,y,z)=n=0M=MT=znxMyMTNMMT(j,n),

    which yields

    fj(x,y,z)=n=0M=MT=znxMyMT{nν1,...,nν2j+1,nπ1,...,nπ2j+1}(δn,nν1+...+nπ2j+1δMT,12(nν1+...nπ2j+1)δM,nν1mν1+...+nπ2j+1mπ2j+1)=2j+1i=1(1+zxmνiy12)(1+zxmπiy12)=fj1(x,y,z)gi(x,y,z)

    with

    fj(x,y,z)=2ji=2(1+zxmνiy12)(1+zxmπiy12)gj(x,y,z)=(1+zxjy12)(1+zxjy12)×(1+zxjy12)(1+zxjy12).

    The generating function is expanded in powers of z as follows:

    fj(x,y,z)=n=0znfj,n(x,y)

    with

    fj,n(x,y)=1n!nznfj(x,y,z)|z=0=M=MT=NMMT(j,n)xMyMT.

    Using the Leibniz formula, one obtains

    fj,n(x,y)=1n!nk=0(nk)nkznkfj1(x,y,z)×kzkgj(x,y,z)|z=0.

    By using Eqs. (A1), (A5), the above relation yields

    M=MT=NMMT(j,n)xMyMT=1n!nk=0(nk)M=MT=NMMT(j1,nk)xMyMT×kzkgi(x,y,z)|z=0.

    Given that gj(x,y,z) is a fourth-order polynomial, the only derivatives that are nonzero correspond to

    k=0:0z0gj(x,y,z)|z=0=1,k=1:1z1gj(x,y,z)|z=0=xjy12+xjy12+xjy12+xjy12,k=2:2z2gj(x,y,z)|z=0=y+y1+x2j+x2j+1+1,k=3:3z3gj(x,y,z)|z=0=xjy12+xjy12+xjy12+xjy12,

    k=4:4z4gj(x,y,z)|z=0=1.

    Eq. (A7) holds for any independent variables x and y; therefore, the coefficients of homogeneous terms on both sides are equal, which immediately leads to Eq. (3).

    In this Appendix, we tabulate DIT(j,n) with j=3/2– 9/2 and n=3–4, for convenience.

    Table B1

    Table B1.  DIT(j,n) for j=5/2-9/2, with n=3. ID in this Table represents that there are D states with spin I for the jn configuration.
    jTID
    5212(12)1,(32)1,(52)2,(72)2,(92)1,(112)1,(132)1
    5232(32)1,(52)2,(92)1
    7212(12)1,(32)1,(52)2,(72)3,(92)2,
    (112)2,(132)2,(152)1,(172)1,(192)1
    7232(32)1,(52)2,(72)3,(92)1,(112)1,(152)1
    9212(12)1,(32)1,(52)2,(72)3,(92)3,(112)3,
    (132)3,(152)2,(172)2,(192)2,(212)1,(232)1,(252)1
    9232(32)1,(52)1,(72)1,(92)2,(112)1,(132)1,(152)1,(172)1,(212)1
    DownLoad: CSV
    Show Table

    Table B2

    Table B2.  The same as Table B1 except for n=4.
    jTID
    52002,23,31,43,51,62,81
    52112,22,33,42,52,61,71
    52201,21,41
    72003,24,32,45,52,65,72,83,91,102,121
    72113,23,35,44,55,64,74,82,92,101,111
    72201,22,42,51,61,81
    92003,26,32,47,54,67,74,
    87,93,105,112,123,131,142,161
    92114,24,37,46,58,67,78,
    86,96,104,114,122,132,141,151
    92202,22,31,43,51,63,71,82,91,101,121
    DownLoad: CSV
    Show Table
    [1] A. de-Shalit and I. Talmi, Nuclear Shell Model Theory (Academic, New York, 1963)
    [2] G. Racah, Phys. Rev. 63, 367 (1943) doi: 10.1103/PhysRev.63.367
    [3] J. Katriel, R. Paunce, and J. J. C. Mulder, Int. J. Quantum Chem. 23, 1855 (1983); J. Katriel and A. Novoselsky, J. Phys. A 22, 1245 (1989)
    [4] D. K. Sunko and D. Svrtan, Phys. Rev. C 31, 1929 (1985)
    [5] D. K. Sunko, Phys. Rev. C 33, 1811 (1986) doi: 10.1103/PhysRevC.33.1811
    [6] D. K. Sunko, Phys. Rev. C 35, 1936 (1987) doi: 10.1103/PhysRevC.35.1936
    [7] J. N. Ginocchio and W. C. Haxton, in Symmetries in Science VI, edited by B. Gruber and M. Ramek (Plenum Press, New York, 1993), p. 263
    [8] Y. M. Zhao and A. Arima, Phys. Rev. C 68, 044310 (2003) doi: 10.1103/PhysRevC.68.044310
    [9] I. Talmi, Phys. Rev. C 72, 037302 (2005) doi: 10.1103/PhysRevC.72.037302
    [10] H. Jiang, F. Pan, Y. M. Zhao et al., Phys. Rev. C 87, 034313 (2013) doi: 10.1103/PhysRevC.87.034313
    [11] Y. M. Zhao and A. Arima, Phys. Rev. C 71, 047304 (2005) doi: 10.1103/PhysRevC.71.047304
    [12] L. Zamick and A. Escuderos, Phys. Rev. C 71, 054308 (2005) doi: 10.1103/PhysRevC.71.054308
    [13] M. Bao, Y. M. Zhao, and A. Arima, Phys. Rev. C 93, 014307 (2016)
    [14] Y. M. Zhao and A. Arima, Phys. Rev. C 70, 034306 (2004)
    [15] Y. M. Zhao and A. Arima, Phys. Rev. C 72, 054307 (2005) doi: 10.1103/PhysRevC.72.054307
    [16] J. -C. Pain, Phys. Rev. C 84, 047303 (2011) doi: 10.1103/PhysRevC.84.047303
    [17] J. -C. Pain, Phys. Rev. C 97, 064311 (2018) doi: 10.1103/PhysRevC.97.064311
    [18] J. -C. Pain, Phys. Rev. C 99, 054321 (2019) doi: 10.1103/PhysRevC.99.054321
    [19] M. Poirier and J. -C. Pain, Phys. Rev. C 104, 064324 (2021) doi: 10.1103/PhysRevC.104.064324
    [20] L. H. Zhang, Y. M. Zhao, L. Y. Jia et al., Phys. Rev. C 77, 014301 (2008)
    [21] L. Zamick and A. Escuderos, Phys. Rev. C 71, 014315 (2005) doi: 10.1103/PhysRevC.71.014315
    [22] L. Zamick and A. Escuderos, Ann. Phys. (NY) 321, 987 (2006) doi: 10.1016/j.aop.2005.08.006
    [23] A. Escuderos and L. Zamick, Phys. Rev. C 73, 044302 (2006)
    [24] L. Zamick, Phys. Rev. C 75, 064305 (2007)
    [25] L. Zamick and P. Van Isacker, Phys. Rev. C 78, 044327 (2008) doi: 10.1103/PhysRevC.78.044327
    [26] L. Zamick and A. Escuderos, Phys. Rev. C 86, 047306 (2012) doi: 10.1103/PhysRevC.86.047306
    [27] L. Zamick and A. Escuderos, Phys. Rev. C 87, 044302 (2013) doi: 10.1103/PhysRevC.87.044302
    [28] P. Van Isacker and S. Heinze, Phys. Rev. Lett. 100, 052501 (2008) doi: 10.1103/PhysRevLett.100.052501
    [29] A. Leviatan, Prog. Part. Nucl. Phys. 66, 93 (2011) doi: 10.1016/j.ppnp.2010.08.001
    [30] C. Qi, X. B. Wang, Z. X. Xu et al., Phys. Rev. C 82, 014304 (2010)
    [31] C. Qi, Phys. Re. C 83, 014307 (2011) doi: 10.1103/PhysRevC.83.014307
    [32] C. Qi, Z. X. Xu, and R. J. Liotta, Nucl. Phys. A 884-885, 21 (2012) doi: 10.1016/j.nuclphysa.2012.04.007
    [33] X. B. Wang and F. R. Xu, Phys. Rev. C 85, 034304 (2012) doi: 10.1103/PhysRevC.85.034304
    [34] Y. Y. Cheng, J. J. Shen, G. J. Fu et al., Phys. Rev. C 100, 014318 (2019) doi: 10.1103/PhysRevC.100.014318
    [35] L. Zamick and A. Escuderos, Phys. Rev. C 72, 044317 (2005) doi: 10.1103/PhysRevC.72.044317
    [36] Y. M. Zhao and A. Arima, Phys. Rev. C 72, 064333 (2005) doi: 10.1103/PhysRevC.72.064333
  • [1] A. de-Shalit and I. Talmi, Nuclear Shell Model Theory (Academic, New York, 1963)
    [2] G. Racah, Phys. Rev. 63, 367 (1943) doi: 10.1103/PhysRev.63.367
    [3] J. Katriel, R. Paunce, and J. J. C. Mulder, Int. J. Quantum Chem. 23, 1855 (1983); J. Katriel and A. Novoselsky, J. Phys. A 22, 1245 (1989)
    [4] D. K. Sunko and D. Svrtan, Phys. Rev. C 31, 1929 (1985)
    [5] D. K. Sunko, Phys. Rev. C 33, 1811 (1986) doi: 10.1103/PhysRevC.33.1811
    [6] D. K. Sunko, Phys. Rev. C 35, 1936 (1987) doi: 10.1103/PhysRevC.35.1936
    [7] J. N. Ginocchio and W. C. Haxton, in Symmetries in Science VI, edited by B. Gruber and M. Ramek (Plenum Press, New York, 1993), p. 263
    [8] Y. M. Zhao and A. Arima, Phys. Rev. C 68, 044310 (2003) doi: 10.1103/PhysRevC.68.044310
    [9] I. Talmi, Phys. Rev. C 72, 037302 (2005) doi: 10.1103/PhysRevC.72.037302
    [10] H. Jiang, F. Pan, Y. M. Zhao et al., Phys. Rev. C 87, 034313 (2013) doi: 10.1103/PhysRevC.87.034313
    [11] Y. M. Zhao and A. Arima, Phys. Rev. C 71, 047304 (2005) doi: 10.1103/PhysRevC.71.047304
    [12] L. Zamick and A. Escuderos, Phys. Rev. C 71, 054308 (2005) doi: 10.1103/PhysRevC.71.054308
    [13] M. Bao, Y. M. Zhao, and A. Arima, Phys. Rev. C 93, 014307 (2016)
    [14] Y. M. Zhao and A. Arima, Phys. Rev. C 70, 034306 (2004)
    [15] Y. M. Zhao and A. Arima, Phys. Rev. C 72, 054307 (2005) doi: 10.1103/PhysRevC.72.054307
    [16] J. -C. Pain, Phys. Rev. C 84, 047303 (2011) doi: 10.1103/PhysRevC.84.047303
    [17] J. -C. Pain, Phys. Rev. C 97, 064311 (2018) doi: 10.1103/PhysRevC.97.064311
    [18] J. -C. Pain, Phys. Rev. C 99, 054321 (2019) doi: 10.1103/PhysRevC.99.054321
    [19] M. Poirier and J. -C. Pain, Phys. Rev. C 104, 064324 (2021) doi: 10.1103/PhysRevC.104.064324
    [20] L. H. Zhang, Y. M. Zhao, L. Y. Jia et al., Phys. Rev. C 77, 014301 (2008)
    [21] L. Zamick and A. Escuderos, Phys. Rev. C 71, 014315 (2005) doi: 10.1103/PhysRevC.71.014315
    [22] L. Zamick and A. Escuderos, Ann. Phys. (NY) 321, 987 (2006) doi: 10.1016/j.aop.2005.08.006
    [23] A. Escuderos and L. Zamick, Phys. Rev. C 73, 044302 (2006)
    [24] L. Zamick, Phys. Rev. C 75, 064305 (2007)
    [25] L. Zamick and P. Van Isacker, Phys. Rev. C 78, 044327 (2008) doi: 10.1103/PhysRevC.78.044327
    [26] L. Zamick and A. Escuderos, Phys. Rev. C 86, 047306 (2012) doi: 10.1103/PhysRevC.86.047306
    [27] L. Zamick and A. Escuderos, Phys. Rev. C 87, 044302 (2013) doi: 10.1103/PhysRevC.87.044302
    [28] P. Van Isacker and S. Heinze, Phys. Rev. Lett. 100, 052501 (2008) doi: 10.1103/PhysRevLett.100.052501
    [29] A. Leviatan, Prog. Part. Nucl. Phys. 66, 93 (2011) doi: 10.1016/j.ppnp.2010.08.001
    [30] C. Qi, X. B. Wang, Z. X. Xu et al., Phys. Rev. C 82, 014304 (2010)
    [31] C. Qi, Phys. Re. C 83, 014307 (2011) doi: 10.1103/PhysRevC.83.014307
    [32] C. Qi, Z. X. Xu, and R. J. Liotta, Nucl. Phys. A 884-885, 21 (2012) doi: 10.1016/j.nuclphysa.2012.04.007
    [33] X. B. Wang and F. R. Xu, Phys. Rev. C 85, 034304 (2012) doi: 10.1103/PhysRevC.85.034304
    [34] Y. Y. Cheng, J. J. Shen, G. J. Fu et al., Phys. Rev. C 100, 014318 (2019) doi: 10.1103/PhysRevC.100.014318
    [35] L. Zamick and A. Escuderos, Phys. Rev. C 72, 044317 (2005) doi: 10.1103/PhysRevC.72.044317
    [36] Y. M. Zhao and A. Arima, Phys. Rev. C 72, 064333 (2005) doi: 10.1103/PhysRevC.72.064333
  • 加载中

Cited by

1. Poirier, M., Pain, J.-C. Exact expressions of the distributions of total magnetic quantum number and angular momentum in single- j orbits: A general technique for any number of fermions[J]. Physical Review C, 2024, 109(2): 024306. doi: 10.1103/PhysRevC.109.024306

Tables(3)

Get Citation
X. Yin and Y. M. Zhao. Number of states with given spin I and isospin T for n fermions in a j orbit[J]. Chinese Physics C. doi: 10.1088/1674-1137/ac7eb3
X. Yin and Y. M. Zhao. Number of states with given spin I and isospin T for n fermions in a j orbit[J]. Chinese Physics C.  doi: 10.1088/1674-1137/ac7eb3 shu
Milestone
Received: 2022-06-11
Article Metric

Article Views(1415)
PDF Downloads(45)
Cited by(1)
Policy on re-use
To reuse of subscription content published by CPC, the users need to request permission from CPC, unless the content was published under an Open Access license which automatically permits that type of reuse.
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Email This Article

Title:
Email:

Number of states with given spin I and isospin T for n fermions in a j orbit

    Corresponding author: Y. M. Zhao, ymzhao@sjtu.edu.cn
  • 1. Shanghai Key Laboratory of Particle Physics and Cosmology, School of Physics and Astronomy, Shanghai Jiao Tong University, Shanghai 200240, China
  • 2. Collaborative Innovation Center of IFSA (CICIFSA), Shanghai Jiao Tong University, Shanghai 200240, China

Abstract: In this study, we investigate formulas of the number of states with a given total spin I and isospin T for n nucleons in a single-j shell denoted by DIT(j,n). Talmi's recursion formulas for the number of states with a given z-axis projection of total spin are generalized by further considering the isospin couplings and are applied to derive explicit formulas of DIT(j,n).

    HTML

    I.   INTRODUCTION
    • Determining the number of states with a given spin I for identical nucleons in a single-j shell [denoted by DI(j,n)] is a common practice in nuclear structure theory and atomic physics. In the nuclear shell model, DI(j,n) is usually tabulated for relatively small j and n. In cases with larger j or n, DI(j,n) might be obtained by subtracting the number of states with total angular momentum projection M=I+1 from that with M=I [1], where M equals m1+m2++mn, and mi [i=1,2,,n] is the projection of j on the z-axis for the i-th nucleon. Other recipes to extract DI(j,n) include Racha's formulas in terms of the seniority scheme [2] and generating function method studied extensively by Katriel et al. [3] and Sunko and collaborators [46].

      The first explicit and compact formula of DI(j,n) was obtained specifically for I=0 and n=4 by Ginocchio and Haxton while studying the fractional quantum Hall effect [7]. In Ref. [8], DI(j,n) was constructed empirically for n=3 and 4, and for a few I values with n=5. The formulas for n=3 was soon proved by Talmi, who introduced his recursive formulas [9]. These formulas were later derived based on the reduction rule from SU(4) to SO(3) group [10]; the formula for n=4 was derived by the reduction rule from SU(5) to the SO(3) group [11], with a demonstration that DI(j,n) can be actually derived based on the reduction rule from SU(n+1) to SO(3) group. The Ginocchio-Haxton formula of D0(j,4) was also revisited by Zamick and Escuderos [12]. An explicit recursion formula from n1 particles to n particles was obtained in Ref. [13] by introducing "pseudo" particles which allow fermions take integer spins. In Refs. [14, 15], the number of states, denoted by DI(j,n), was applied to derive sum rules of six-j and nine-j symbols, some of which were also revisited in Ref. [ 16]. In the last decade, Pain and collaborators extensively studied the odd-even staggering of DI(j,n) [17, 18], and compact formulas of DI(j,n) for n=3,4, and 5 [19]. The enumeration of the number of states with given spin was also extended to boson systems in Ref. [20]. The study of DI(j,n) motivated a number of studies related to a single-j shell; here, we mention Refs. [2134] without providing further details for completeness.

      Given that the nuclear shell model Hamiltonian respects the isospin symmetry, it is natural to generalize the enumeration of DI(j,n) to the number of states with given spin I and isospin T, denoted by DIT(j,n), for nucleons in a single-j shell. In Ref. [35], Zamick and Escuderos found a few simple relations between DIT(j,4) with T=0 and DIT(j,4) with T=2. In Ref. [36], compact and explicit formulas of DIT(j,n) for n=3 and 4 were derived in terms of sum rules of six-j and nine-j symbols involving the expression of DI(j,n) given in Ref. [11].

      Similar to the enumeration of DI(j,n), DIT(j,n) can be obtained in terms of the number of states with given spin projection M and isospin projection MT. This number is denoted by NMMT(j,n). For n fermions in a single-j shell,

      DIT(j,n)=[NM=I,MT=T(j,n)NM=I+1,MT=T(j,n)][NM=I,MT=T+1(j,n)NM=I+1,MT=T+1(j,n)]=NM=I,MT=T(j,n)+NM=I+1,MT=T+1(j,n)NM=I+1,MT=T(j,n)NM=I,MT=T+1(j,n).

      (1)

      In this study, we investigate DIT(j,n) in terms of NMMT(j,n), which is obtained by generalizing Talmi's recursion formula of NM(j,n), in which isospin is not considered. Compact formulas of DIT(j,n) are derived for n=3 as an exemplification. The rest of this paper is organized as follows. In Sec. II, we generalize Talmi's formulas for NM(j,n) by further considering the isospin symmetry, and present recursion formulas of NMMT. In Sec. III, we apply our generalized recursion formulas of NMMT and extract NMMT(n,j) for n=3. In Sec. IV, DIT(j,3) is derived by using Eq. (1). In Sec. V, we summarize our study. In Appendix A, we provide an alternative mathematical proof of our generalized Talmi's recursion formula by using the generating function method. In Appendix B, we tabulate DIT(j,n) with n=3 and 4, j=5/2,7/2, and 9/2.

    II.   NUMBER OF STATES WITH GIVEN z-COMPONENT PROJECTIONS OF SPIN AND ISOSPIN
    • We use the convention that the z-axis projection of the isospin for a neutron and a proton equals 1/2 and 1/2, respectively. For n nucleons in a single-j shell, we denote the projections of total spin I and total isospin T, respectively, as M and MT. With this convention, the neutron number equals n/2+MT, and the proton number equals n/2MT. Let us denote the z-axis spin projections of neutrons and protons by using (mν1,mν2,...,mν2j+1) and (mπ1,mπ2,...,mπ2j+1), respectively, with the convention that j=mν1>mν2>...>mν2j+1=j for neutrons and j=mπ1>mπ2>...>mπ2j+1=j for protons. We define nνi=1 (or nπi=1) for i=1,...,2j+1 if the orbit of mνi (or mπi) is filled; otherwise, it equals zero.

      According to the above conventions, we have

      M=nν1mν1+nν2mν2+...+nν2j+1mν2j+1+nπ1mπ1+nπ2mπ2+...+nπ2j+1mπ2j+1,

      and the maximum of I equals

      Imax=Mmax=(j+(j1)+...+(jn/2MT+1))+(j+(j1)+...+(jn/2+MT+1))=n(2j+1)/2n2/4M2T.

      (2)

      We denote the number of states with given M and MT for n nucleons in a single-j shell by using NMMT(j,n). Clearly, we can divide the values of (nν1,nν2j+1,nπ1,nπ2j+1) into sixteen cases, as tabulated in Table 1. Here, we exemplify them by using ( nν1=1,nν2j+1=1,nπ1=1,nπ2j+1=1). In this case, there are (n4) nucleons distributed in orbits for which the absolute values of mi values equal or are below (j1). Consequently, the number of states with given M and MT, i.e., NMMT(j,n), equals the number of states of the same M and MT, but with (n4) nucleons in a (j1) shell, namely, NMMT(j1,n4). According to this classification, we are able to obtain NMMT(j,n) by using the values of NMMT(j,n), where either j is smaller than j, or n is smaller than n, or both (j,n) are smaller than (j,n). In other words, from Table 1 we have

      nν1nν2j+1nπ1nπ2j+1Number of states
      1111NM,MT(j1,n4)
      1110NMj,MT1/2(j1,n3)
      1101NM+j,MT1/2(j1,n3)
      1100NM,MT1(j1,n2)
      1011NMj,MT+1/2(j1,n3)
      1010NM2j,MT(j1,n2)
      1001NM,MT(j1,n2)
      1000NMj,MT1/2(j1,n1)
      0111NM+j,MT+1/2(j1,n3)
      0111NM,MT(j1,n2)
      0101NM+2j,MT(j1,n2)
      0100NM+j,MT1/2(j1,n1)
      0011NM,MT+1(j1,n2)
      0010NMj,MT+1/2(j1,n1)
      0001NM+j,MT+1/2(j1,n1)
      0000NM,MT(j1,n)

      Table 1.  Number of states corresponding to the available values of nν1,nν2j+1,nπ1,nπ2j+1.

      NMMT(j,n)={nν1,nν2j+1,nπ1,nπ2j+1}NMj(nν1nν2j+1+nπ1nπ2j+1),MT12(nν1+nν2j+1nπ1nπ2j+1)(j1,nnν1nν2j+1nπ1nπ2j+1).

      (3)

      Here, the summation over {nν1,nν2j+1,nπ1,nπ2j+1} refers to the sixteen cases listed in Table 1. Note that the above formula holds also for n=1,2,3,4, with the convention that NMMT(j,n)=0 if n<0 and NMMT(j,0)=δM,0δMT,0. Note also that Eq. (3) can be proved alternatively based on the generating-function technique described in Ref. [3], as presented in Appendix A.

    III.   RECURSION FORMULA OF NMMT(j,3)
    • The value of NMMT(j,n) with n=1 or 2 is trivially determined as follows. For MT=±1/2 and |M|j, NM,MT=±12(j,1)=1; otherwise, it equals zero. For n=2 with MT=±1, the highest value of M equals Mmax=2j1 and for MT=0, Mmax=2j. For MT=±1 and M2j1,

      NM,MT=±1(j,2)=2j+1|M|2,

      (4)

      where ""and"" means taking the largest integer without exceeding the value inside; for MT=0 and M2j,

      NM,MT=0(j,2)=2j+1|M|.

      (5)

      Otherwise, NMMT(j,2)=0. With these results, we can construct explicit formulas of NMMT(j,3) by using Eq. (3).

    • A.   MT=1/2 and j<MMmax

    • We first discuss NMMT(j,3) with the requirement that MT=12 and j<MMmax. We recursively apply Eq. (3) to this term and obtain the value of NMMT(j,n):

      NMMT(j,n)=NMMT(j1,n)+{nν1,nν2j+1,nπ1,nπ2j+1}NMj(nν1nν2j+1+nπ1nπ2j+1),MT12(nν1+nν2j+1nπ1nπ2j+1)(j1,nnν1nν2j+1nπ1nπ2j+1)=NMMT(jl,n)+l1i=0{nν1,nν2j+1,nπ1,nπ2j+1}NM(ji)(nν1nν2j+1+nπ1nπ2j+1),MT12(nν1+nν2j+1nπ1nπ2j+1)(ji1,nnν1nν2j+1nπ1nπ2j+1),

      (6)

      where denotes summations that exclude the case {nν1,nν2j+1,nπ1,nπ2j+1}={0,0,0,0} and decouple the term NM,MT(j1,n) [the last term in the summation of the sixteen terms of Table 1]; l is fixed with the requirement that NMMT(jl,n)=0 while NMMT(jl+1,n)0, which means that

      M>n(2j+1)2n24M2Tln,Mn(2j+1)2n24M2T(l1)n,

      or equivalently,

      l=2j+12n4M2T+Mn+1.

      (7)

      For n=3 and MT=1/2, the value of l becomes

      l=jM+13+1.

      (8)

      According to Eq. (6), with the conditions NMMT(jl,n)=0 and NMMT(jl+1,n)0, we have

      NM,MT=12(j,3)=l1i=0{nν1,nν2j+1,nπ1,nπ2j+1}NM(ji)(nν1nν2j+1+nπ1nπ2j+1),1212(nν1+nν2j+1nπ1nπ2j+1)(ji1,3nν1nν2j+1nπ1nπ2j+1).

      (9)

      On the right hand side of the above formula, there are many terms that break the requirement of Eq. (2) for n2. Let us exemplify this by using two cases: {nν1,nν2j+1,nπ1,nπ2j+1}={0,0,0,1} and {1,1,1,0}. For the first case, {nν1,nν2j+1,nπ1,nπ2j+1}={0,0,0,1}, the corresponding contribution to NM,MT=12(j,3) on the right hand of Eq. (10) is NM+(ji),MT=1(ji1,2). Given that M>j here, the value of M+(ji) is larger than the maximum Mmax of two nucleons in a single ji1 shell, i.e., Mmax=2j2i2. Similarly, for {nν1,nν2j+1,nπ1,nπ2j+1}={1,1,1,0}, the corresponding contribution on the right hand side is NM(ji),MT=0(ji1,0). Given that M(ji)>0, NM(ji),MT=0(ji1,0)=δM(ji),0=0. As a result, we obtain

      NM,MT=12(j,3)=l1i=0(NM(ji),MT=0(ji1,2)+NM(ji),MT=1(ji1,2)+NM2(ji),MT=12(ji1,1))=l1i=0(3j3i1M+3j3i1M2+1)=(3j1M)l3l(l1)2+l1i=03j3i1M2+l.

      (10)

      The summation of 3j3i1M2 is tedious but straightforward. One has

      l1={2kifwhen3jM=6k+1,6k+2,6k+3,2k+1ifwhen3jM=6k+4,6k+5,6k+6.

      When 3jM=6k+1, the summation becomes

      l1i=03j3i1M2=2ki=0(3k2i+i2)=2ki=0(3k2i)+2k1i=0i+k=3k2+k.

      This summation for the case 3jM=6k+2,3,4,5,6 can be obtained similarly, as follows:

      l1i=03j3i1M2={3k2+kif3jM=6k+1,3k2+2kif3jM=6k+2,3k2+3k+1if3jM=6k+3,3k2+4k+1if3jM=6k+4,3k2+5k+2if3jM=6k+5,3k2+6k+3if3jM=6k+6.

      (11)

      Substituting these results into Eq. (10), we obtain

      NM,MT=12(j,3)={9k2+6k+1if3jM=6k+1,9k2+9k+2if3jM=6k+2,9k2+12k+4if3jM=6k+3,9k2+15k+6if3jM=6k+4,9k2+18k+9if3jM=6k+5,9k2+21k+12if3jM=6k+6=14(3jM)2+12(3jM)+η(3jM),

      (12)

      where η(3jM)=(0,14) if mod(3jM,2)=(0,1), respectively.

    • B.   Tz=1/2 and 0<Mj

    • Let us now address NM,MT=12(j,3) with Tz=1/2 and 0<Mj. In this case, Eq. (3) is applied recursively for jM+1 times. As a result, we have

      NMMT(j,n)=NMMT(j1,n)+{nν1,nν2j+1,nπ1,nπ2j+1}NMj(nν1nν2j+1+nπ1nπ2j+1),MT12(nν1+nν2j+1nπ1nπ2j+1)(j1,nnν1nν2j+1nπ1nπ2j+1)=NMMT(M1,n)+jMi=0{nν1,nν2j+1,nπ1,nπ2j+1}NM(ji)(nν1nν2j+1+nπ1nπ2j+1),MT12(nν1+nν2j+1nπ1nπ2j+1)(ji1,nnν1nν2j+1nπ1nπ2j+1),

      (13)

      where, again, denotes summation among the sixteen sets of {nν1,nν2j+1,nπ1,nπ2j+1} excluding {nν1,nν2j+1,nπ1,nπ2j+1}={0,0,0,0}. For MT=12 and n=3, the above formula is reduced to

      NM,MT=12(j,3)=1+NM,MT=12(M1,3)+jMi=0(NM(ji),MT=0(ji1,2)+NM+(ji),MT=0(ji1,2)+NM(ji),MT=1(ji1,2)+NM+(ji),MT=1(ji1,2)+NM,MT=12(ji1,1)+NM2(ji),MT=12(ji1,1)+NM,MT=12(ji1,1)+NM,MT=12(ji1,1)+NM+2(ji),MT=12(ji1,1)+NM,MT=32(ji1,1))=1+NM,MT=12(M1,3)+jMi=0(ji+M1)+jM2i=0(jiM1)

      +jMi=0ji+M12+jM3i=0jiM12+3jM1i=01,

      (14)

      where the first term " 1" is given by NM=0,MT=0(ji1,0)=1 with {nν1,nν2j+1,nπ1,nπ2j+1,i}={1,1,1,0,jM}; the second term NM,MT=12(M1,3) is given in Eq. (12) with j=M1, i.e.,

      NM,MT=12(M1,3)=M22M+34.

      (15)

      The summations in Eq. (14) equal zero when the upper limit of the given summation index is larger than its lower limit. On the right hand side, the summations of ji+M12 and jiM12 are classified into two cases, jM=2k and 2k+1. For jM=2k, we obtain

      jMi=0ji+M12=2ki=02ki+2M12=2ki=0(ki+2M12+i2)=2ki=0(ki+2M12)+2k1i=0i+k=(jM+1)(M12)+(jM2)2.

      The summation with jM=2k+1 is obtained similarly. The results of these two cases are unified as follows:

      jMi=0ji+M12=(jM+14)(M12)+(jM2)2+{0ifmod(jM,2)=0,14ifmod(jM,2)=1.

      (16)

      Following a similar procedure, we obtain

      jM3i=0jiM12=(jM2)2jM2+{0ifmod(jM,2)=0,+14ifmod(jM,2)=1.

      (17)

      Substituting Eqs. (15)–(17) into Eq. (14), we obtain NM,MT=12(j,3) with 0<Mj in a compact form:

      NM,MT=12(j,3)=32j2+j+14M22.

      (18)
    • C.   MT=3/2 and 0<M

    • The results for NMMT(j,n) are reduced to NM(j,n) for identical particles, and Eq. (3) is reduced to Talmi's recursion formula. For this case, NM,MT=32(j,3) was obtained in Ref. [19], which we cite for completeness. For j<MMmax, according to Eqs. (2.11a)–(2.11b) in Ref. [19],

      NM,MT=32(j,3)=112(3jM)2+α(3jM),

      (19)

      where

      α(3jM)=0,112,13,14,13,112,

      for mod(3jM,6)=0,1,2,3,4,5, respectively. For 0<Mj, according to Eqs. (2.19a)-(2.19b) in Ref. [19],

      NM,MT=32(j,3)=12j216M2+β(M12),

      (20)

      where

      β(3jM)=112,14,112,

      for mod(M12,3)=0,1,2, respectively.

    IV.   NUMBER OF STATES WITH GIVEN I AND T
    • Provided the results of NMMT(j,3), i.e., Eqs. (12), (18), (19), and (20), we can readily substitute them into Eq. (1) to obtain the number of states with given I and T, i.e., DIT(j,3).

      We first consider the case T=1/2. When I>j and T=1/2, we obtain

      DI,T=12(j,3)=NI,12(j,3)+NI+1,32(j,3)NI+1,12(j,3)NI,32(j,3).

      (21)

      According to Eqs. (12) and (19), the above formula is reduced to

      DI,T=12(j,3)=13I+j+13+η(3jI)η(3jI1)+α(3jI1)α(3jI)=13I+j+{0ifmod(3jI,3)=0,23ifmod(3jI,3)=1,13ifmod(3jI,3)=2=3jI+23.

      (22)

      For I=j and T=1/2,

      DI,T=12(j,3)=NI,12(j,3)+NI+1,32(j,3)NI+1,12(j,3)NI,32(j,3).

      (23)

      By using Eqs. (12), (18), (19), and (20), we have

      DI,T=12(j,3)=23j+712+α(2j1)η(2j1)β(j12)=23j+{23ifmod(2j,3)=1,0ifmod(2j,3)=3,13ifmod(2j,3)=5=2j+23.

      (24)

      Similarly, for 0<I<j and T=1/2,

      DI,T=12(j,3)=NI,12(j,3)+NI+1,32(j,3)NI+1,12(j,3)NI,32(j,3).

      (25)

      By using Eqs. (18) and (20), the above formula is reduced to

      DI,T=12(j,3)=23J+13+β(I+12)β(I12)={23I+23ifmod(2I,3)=1,23Iifmod(2I,3)=3,23I+13ifmod(2I,3)=5=2I+23.

      (26)

      Eqs. (24)–(26) can be rewritten as follows. For Ij and T=1/2,

      DI,T=12(j,3)=3jI+23,

      (27)

      and for 0<Ij and T=1/2,

      DI,T=12(j,3)=2I+23.

      (28)

      These results [Eqs. (27)–(28) in this paper] are consistent with Eqs. (25)–(26) in Ref. [ 36], which were obtained in terms of sum rules of six-j symbols, except that the formulas derived in this paper are expressed in a more transparent and understandable form.

      The T=3/2 case is much simpler than the T=1/2 case. One uses Talmi's recursion formulas for NM=I,32(j,n) and obtains DI,T=32(j,3) straightforwardly. The resulting formulas are consistent with those given in Ref. [ 8] and Ref. [19]. We show them next for completeness. According to Eq. (1) in Ref. [ 8], for 0<Ij and T=3/2,

      DI,T=32(j,3)=2I+36;

      (29)

      and according to Eq. (2) in Ref. [ 8], for Ij and T=3/2,

      DI,T=32(j,3)=3j3I6+δI.

      (30)

      where

      δI={0if mod(3j3I,6)=1,1otherwise.

      In principle, one can follow the above procedure as explained in Eqs. (1)–(3) for larger n values. However, this procedure becomes formidably tedious with lengthy formulas and tables for larger n values. Instead of those formulas, we calculate DIT(j,n) for n=3 and 4 and j=5/2, 7/2, and 9/2 and tabulate them in Table B1 and Table B2 in Appendix B.

    V.   SUMMARY
    • To summarize, in this paper we generalize Talmi's recursion formula of the number of states with given spin I by further consideration of the isospin symmetry. This generalization is also proved alternatively based on the generating function method. We also exemplify the generalized Talmi's recursion formulas obtained in this paper by applying our generalized Talmi's recursion formulas in constructing the number of states with given spin and isospin of three nucleons in a single-j shell.

    APPENDIX A: GENERATING-FUNCTION TECHNIQUE
    • In this Appendix, we present an alternative proof of the recursion formula presented in Eq. (3) based on the generating-function method.

      A generating function is a polynomial function for which the expansion coefficients are related to the number of states with a given requirement. We define the generating function corresponding to NMMT(j,n), similar to that in Ref. [17],

      fj(x,y,z)=n=0M=MT=znxMyMTNMMT(j,n),

      which yields

      fj(x,y,z)=n=0M=MT=znxMyMT{nν1,...,nν2j+1,nπ1,...,nπ2j+1}(δn,nν1+...+nπ2j+1δMT,12(nν1+...nπ2j+1)δM,nν1mν1+...+nπ2j+1mπ2j+1)=2j+1i=1(1+zxmνiy12)(1+zxmπiy12)=fj1(x,y,z)gi(x,y,z)

      with

      fj(x,y,z)=2ji=2(1+zxmνiy12)(1+zxmπiy12)gj(x,y,z)=(1+zxjy12)(1+zxjy12)×(1+zxjy12)(1+zxjy12).

      The generating function is expanded in powers of z as follows:

      fj(x,y,z)=n=0znfj,n(x,y)

      with

      fj,n(x,y)=1n!nznfj(x,y,z)|z=0=M=MT=NMMT(j,n)xMyMT.

      Using the Leibniz formula, one obtains

      fj,n(x,y)=1n!nk=0(nk)nkznkfj1(x,y,z)×kzkgj(x,y,z)|z=0.

      By using Eqs. (A1), (A5), the above relation yields

      M=MT=NMMT(j,n)xMyMT=1n!nk=0(nk)M=MT=NMMT(j1,nk)xMyMT×kzkgi(x,y,z)|z=0.

      Given that gj(x,y,z) is a fourth-order polynomial, the only derivatives that are nonzero correspond to

      k=0:0z0gj(x,y,z)|z=0=1,k=1:1z1gj(x,y,z)|z=0=xjy12+xjy12+xjy12+xjy12,k=2:2z2gj(x,y,z)|z=0=y+y1+x2j+x2j+1+1,k=3:3z3gj(x,y,z)|z=0=xjy12+xjy12+xjy12+xjy12,

      k=4:4z4gj(x,y,z)|z=0=1.

      Eq. (A7) holds for any independent variables x and y; therefore, the coefficients of homogeneous terms on both sides are equal, which immediately leads to Eq. (3).

    APPENDIX B: DIT(J,N) WITH J=3/29/2 AND N=34
    • In this Appendix, we tabulate DIT(j,n) with j=3/2– 9/2 and n=3–4, for convenience.

      jTID
      5212(12)1,(32)1,(52)2,(72)2,(92)1,(112)1,(132)1
      5232(32)1,(52)2,(92)1
      7212(12)1,(32)1,(52)2,(72)3,(92)2,
      (112)2,(132)2,(152)1,(172)1,(192)1
      7232(32)1,(52)2,(72)3,(92)1,(112)1,(152)1
      9212(12)1,(32)1,(52)2,(72)3,(92)3,(112)3,
      (132)3,(152)2,(172)2,(192)2,(212)1,(232)1,(252)1
      9232(32)1,(52)1,(72)1,(92)2,(112)1,(132)1,(152)1,(172)1,(212)1

      Table B1.  DIT(j,n) for j=5/2-9/2, with n=3. ID in this Table represents that there are D states with spin I for the jn configuration.

      jTID
      52002,23,31,43,51,62,81
      52112,22,33,42,52,61,71
      52201,21,41
      72003,24,32,45,52,65,72,83,91,102,121
      72113,23,35,44,55,64,74,82,92,101,111
      72201,22,42,51,61,81
      92003,26,32,47,54,67,74,
      87,93,105,112,123,131,142,161
      92114,24,37,46,58,67,78,
      86,96,104,114,122,132,141,151
      92202,22,31,43,51,63,71,82,91,101,121

      Table B2.  The same as Table B1 except for n=4.

Reference (36)

目录

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return