arxiv: v2 [nucl-th] 27 Feb 2019

Transcription

1 Beyond Wigner s Isobaric Multiplet Mass Equation: Effect of Charge-Symmetry-Breaking Interaction and Coulomb Polarization J. M. Dong,, J. Z. Gu, Y. H. Zhang, W. Zuo,,, L. J. Wang, Yu.. Litvinov, 5, 6, 7, and Y. Sun Institute of Modern Physics, Chinese cademy of Sciences, Lanzhou 7, China China Institute of tomic Energy, P. O. Box 75), Beijing, China School of Physics, University of Chinese cademy of Sciences, Beijing, China Department of Physics and stronomy, University of North Carolina, Chapel Hill, North Carolina, , US arxiv:6.567v nucl-th] 7 Feb 5 GSI Helmholtzzentrum für Schwerionenforschung, Planckstraße, 6 Darmstadt, Germany 6 School of Physics and stronomy, Shanghai Jiao Tong University, Shanghai, China 7 Collaborative Innovation Center of IFS, Shanghai Jiao Tong University, Shanghai, China Dated: February 8, ) The quadratic form of the isobaric multiplet mass equation IMME), which was originally suggested by Wigner and has been generally regarded as valid, is seriously questioned by recent high-precision nuclear mass measurements. The usual resolution to this problem is to add empirically the cubic and quartic T z -terms to characterize the deviations from the IMME, but finding the origin of these terms remains an unsolved difficulty. Based on a strategy beyond the Wigner s first-order perturbation, we derive explicitly the cubic and quartic T z -terms. These terms are shown to be generated by the effective charge-symmetry breaking and chargeindependent breaking interactions in nuclear medium combined with the Coulomb polarization effect. Calculations for the sd- and lower f p-shells explore a systematical emergence of the cubic T z -term, suggesting a general deviation from the original IMME. Intriguingly, the magnitude of the deviation exhibits an oscillationlike behavior with mass number, modulated by the shell effect. PCS numbers:.8.+y,.75.cs,.65.ef,..dr I. INTRODUCTION Shortly after the discovery of neutron, Heisenberg introduced isospin to describe different charge states of nucleon ]. In this concept, proton p) and neutron n) are treated as an isospin T= / doublet distinguished by different projections T z p)= / and T z n)=+/. s one of the most important predictions in nuclear physics, the isobaric multiplet mass equation IMME) proposed later by Wigner, ] suggests that the mass excesses ME, T, T z ) of the nuclei belonging to an isospin multiplet of mass number and total isospin T follow a simple quadratic equation ME, T, T z )=a+bt z + ct z, ) where T z = N Z)/ is the isospin projection, and the parameters a, b and c are constants for a given multiplet. The elegant IMME, though derived by using the st-order perturbation approximation, has been widely employed to predict the unknown masses of unstable neutron-deficient nuclei. Since its establishment, the IMME is believed to be generally valid ]. With recent advances in radioactive beam facilities, a wealth of exotic masses with increasing precision became available 5]. Unexpectedly large discrepancies between the measured masses and the ones given by the quadratic form of the IMME were observed 6 8]. This calls for an addition of a cubic term dt z or even a quartic term et z to Eq. ) ]. The origin of these higher-order terms, which clearly lies beyond the original IMME of Eq. ), requires explanation. dongjm7@impcas.ac.cn zuowei@impcas.ac.cn sunyang@sjtu.edu.cn

2 Various mechanisms have been proposed to explain the deviations found in individual cases, including the isospin mixing, the high-order Coulomb effect, and the charge-dependent nucleon-nucleon interaction 5]. However, to date there is no consensus as to the origin of the observed large dt z terms. In current shell-model calculations, isospin-nonconversing INC) interactions are determined through fitting to available experimental data 6 ], which are however insufficient to explain the experimental dt z terms. In a recent work ], we have laid out a theoretical framework which considers the contributions of the charge-symmetry breaking CSB) and charge-independence breaking CIB) components in nuclear medium to the effective nucleon-nucleon force. We have found that such effective INC interactions are density-dependent, and thus can no longer be expressed as irreducible tensors as was done by Wigner, ]. This leads us to propose a generalized IMME GIMME) ] to study the Nolen-Schiffer anomaly, which is expressed as ] ME, T, T z ) = a+ b c + nh + a CSB) sym,, T z) ) T z + c c + ) acib) sym,, T z) Tz, ) where nh =.78 MeV is the neutron-hydrogen mass difference. In Eq. ), the two anticipated isospin-symmetry breaking sources, the Coulomb and the nuclear interactions, are clearly separated. The T z -independent b c and c c coefficients are produced solely by the Coulomb interaction, whereas the st nd)-order symmetry energy a CSB) sym, acib) sym, ) originates from the CSB CIB) interaction in nuclear medium ]. In this Rapid Communication we apply the GIMME ] to explore the physics beyond the st-order perturbation in the IMME. We show that the appearance of the high-order T z terms is a more general phenomenon, and that the degree of the deviation to Eq. ), measured by the coefficient of the Tz term, is totally governed by shell effects with a remarkable - dependence. II. T= / ISOBRIC QURTETS Within the st-order perturbation calculation, αtt z is assumed to be eigenstate of the charge-independent Hamiltonian H, withαfor all additional quantum numbers to label this state. The energy produced by the CSB and CIB interactions is given by αtt z H CSB+CIB αtt z =a CSB) sym,, T z)t z + acib) sym,, T z)tz in the absence of the Coulomb force ]. In this case, the a CSB) sym,, T z) and also a CIB) sym,, T z)) is identical for all members of an isobaric multiplet, as seen in the appendix. Thus, the GIMME is reduced to the quadratic form of the IMME. However, if we go beyond the st-order perturbation to calculate αt z H CSB+CIB αt z with inclusion of the Coulomb polarization effect, a CSB) sym,, T z) and also a CIB) sym,, T z)) is no longer a constant for a given isobaric multiplet. One may formally expand them to be a CSB) sym,, T z)t z = a +b T z +c Tz+d Tz and a CIB) sym,, T z)tz= a + b T z + c Tz + d Tz for the T= / quartets. Equation ) can then be rearranged as where the d coefficient is explicitly expressed as d = 8π ME, T, T z )=a+bt z + ct z+ dt z, ) 8π r S CSB) ρ) δρ / δρ / ) dr r ρr) S CIB) ρ) δρ / δρ /) dr. ) In Eq. ),δρ Tz = ρ n r) ρ p r) ] core is the neutron- and proton-density difference in the core of the nucleus with T z. For T= T z T z T z ) nuclei, where there are N Z excess neutrons protons), we call the rest nucleons with an equal number of protons and neutrons) the core in the discussion. The terms S CSB) and S CIB) in Eq. ) are density-dependent symmetry energies

3 characterizing the INC interactions, defined as S CSB) ρ) = ECSB) ρ,β) β=, β 5) S CIB) ρ) = E CIB) ρ,β) β β=. 6) The above results are achieved based on the microscopic Brueckner-Hartree-Fock approach with the V8 together with V) interaction ], in which E CSB) E CIB) ) is the energy associated with the CSB CIB) interaction. The CIB interaction contributes about one order of magnitude smaller than the CSB interaction. From Eq. ), the nonzero d originates primarily from the storder symmetry energy difference between the T z = / core and T z = / core. The nucleonic density distributions inside the nucleus are calculated within a Skyrme energy-density functional approach. Detailed derivation of Eqs., ) are given in the appendix. We note that the dt z term is usually added empirically to account for deviations from the quadratic form of the IMME, but here, it is derived microscopically. The occurrence of the nonzero d coefficient is due to a combined effect of the CSB interaction in nuclear medium together with the treatment of the beyond-st-order perturbation calculation i.e., by including the core polarization induced by the Coulomb force). Both of them are indispensable. We thus conclude that although the high-order Coulomb contribution to the d value is generally believed to be small ], it is necessarily to be included because without it, as we discussed in detail before Eq. ), d term would not appear. d kev) 5-8] ] Exp. Exp.new) Cal. ] ] 8] ] 6] FIG. : Color online) d values for T = / quartets calculated with Eq. ). The presented experimental data are from Refs., ] for =, 8] for,, ] for, 5, 6] for, 7, 8] for, ] for 5, ] for 7, 7] for, 7, ] for 5, and the remaining data are taken from Ref. ]. The numbers in brackets refer to the years of publication. ] In Fig., we show the d values calculated with Eq. ) for T= / isobaric quartets. The SLy interaction ] that satisfies specific constraints defined in our previous work ] is employed to compute the nucleonic density distributions. Remarkably, the results for the d coefficients exhibit a clear -dependence. It is striking that across the sd and lower f p shells, the d values show an oscillation-like behavior with a minimum 6 kev and maximum 7 kev. These occur across the magic numbers 8,, and 8, implying the shell effect behind. Experimental data seem to support the above prediction. In Fig., we also plot the experimentally-extracted d values for the isobaric quartets from the measured mass excesses via d=met z = /) MET z = /) MET z = /)+MET z = /)]/6. These experimental d-values follow the predicted pattern well. s one can see from Fig., the experimental d-values for = 5, 6], 7, 8], and 5 ] multiplets have small uncertainties. Our calculations reproduce these data. With the new mass of S measured with the isochronous mass spectrometry technique in CSRe recently ], the IMME is shown to be revalidated for =, which is supported by our calculation. However, the validity of the IMME represents only special cases of our general conclusion. For the T z = / isobaric quartets, the first test in the f p-shell with =5,, and 5 indicated systematical deviations from the IMME 7]. Calculations based on two INC Hamiltonians, the f 7/ model space ] and the full p f model space with the

4 GPFX interaction 5, 6] plus the Ormand-Brown INC Hamiltonian 6], could not reproduce the experimental data 7]. Our calculated d values for =, 5, and shown in Fig. agree qualitatively with experiment. For =5, Ref. 7] initially reported a very large d=) kev. With a later remeasurement of the IS of 5 Co ], it was reduced and agrees now well with our prediction see Fig. ). From our calculations, significant d-values are expected in the p f -shell. Precision mass measurements of the relevant nuclei are required. Especially interesting would be the confirmation of the maximum d-value at =, for which the identification and mass determination of the T= / IS in Ti are needed. The = multiplet has been taken in Ref. 8] as an example to show violations of the IMME in the sd-shell nuclei. The universal sd USD and USDB isospin-conserving Hamiltonians supplemented with an INC part yield too small d values d=. kev USD) and. kev USDB) for the =, J π = 5/ + quartet 8]). The valence-space calculations based on the low-momentum two-nucleon and three-nucleon forces derived from the chiral effective field theory 7] give d = 8 kev for the = quartet 8], which disagrees with the experimental value 6.7) kev. Our result d=. kev is close to the experimental data. f 7/ d / s / a) 8] ] 6] - fm - ) 6 5 = b) = = - =7 - r fm) 5 6 FIG. : Color online) a) Schematic illustration highlighting the filling pattern of excess neutrons for the T z = T nucleus. b) The core density differenceδρ / δρ / for T= / isobaric quartets. =7 and = correspond to the left and right panels of a), respectively. We stress that the occurrence of nonzero d s, which marks deviations from the original IMME, is fundamental. The variation of the d coefficient is driven by the shell effect. We find that once the excess neutrons in the T z = / of a multiplet fill a level below above) a large shell gap, as schematically illustrated in the left right) panel of Fig. a), the smallest largest) d value appears see Fig. ). For instance, the T z = T member of the =7 =) multiplet, 7 Cl Ca), has three excess neutrons filling below above) the N = shell gap. When changes from 7 to, the excess neutrons protons) for T z = / T z = /) member gradually occupy the upper f 7/ orbit. The neutrons protons) in the core tend to be more loosely tightly) bound if they and the excess neutrons protons) fill the same different) orbits), leading to the -dependent differences in neutron proton) density of the core. s seen in Fig. b), when changes from 7 to,δρ / δρ / = ρ core n r) ρ p r) ] + ρ core T z =/ p r) ρ n r) ] increases significantly, particularly in the region of r= fm. This shell T z = / effect is brought into the d coefficient via the integral in Eq. ). In the T= / multiplets, the excess neutrons protons) in the T z = / T z = /) nuclei occupy a level above the p / d 5/ or d / f 7/ shell gap for the = or =5 5 quartets, respectively, leading to a relatively large violation of the IMME. This mechanism holds also true for the strong breakdown of the = quartet since the excess neutrons protons) in the T z = / T z = /) nuclei occupy a level above the s / p / shell gap. However, the = quartet is excluded from Fig. for discussion because generally, the Skyrme functional does not quantitatively apply to such light nuclei. We conclude that the magnitude of the d coefficient, which measures the degree of deviation from the original IMME, depends on the shell filling. s shown in Fig., we predict a local maximum of d for the = quartet. The existing experimental measurements are divergent for the excitation energy of the IS in Ne, ], leading to completely different conclusions regarding the IMME. n experimental confirmation of this energy is necessary.

5 5 d kev) a) 5] 8 6 Exp. Exp.new) Cal. ] b) 5] ] e kev) - -6 FIG. : Color online) d and e values for T = quintets calculated with Eqs. 8, ). The presented experimental data are taken from Refs. 8, 8] for =, and the remaining data are from Ref. ]. The numbers in brackets refer to the years of publication. III. T= ISOBRIC QUINTETS With the similar derivation for T= / quartets, the GIMME of Eq. ) is rewritten for T= quintets as with the d and e coefficients given as whereρ exc n ME, T, T z )=a+bt z + ct z + dt z + et z, 7) d = π π e = π 6 r S CSB) ρ) δρ δρ ) dr r ρr) S CIB) ρ) δρ ) δρ dr, 8) dr r ρr) S CIB) ρ) δρ δρ ) δρ δρ )+8 ρ exc n Tz =+ρ exc p Tz = )], ) T z = ρ exc p T z = ) is the density of the N Z excess neutrons protons) in the T z = T T z = T) nucleus. The e coefficient just originates from the CIB interaction, and the detailed derivation is presented in the appendix. In Fig., we compare the calculated d and e values for T= quintets with the available experimental data. The calculations suggest that the d values are overall small, and the magnitude does not exceed kev. The calculated d values for = 6 in Fig. a) show a similar pattern as the one for T= / quartets with a maximum at =. The occurrence of this variation shares the same physical origin as that in the T= / quartets, namely the excess neutrons in the T z = T= member occupy a level above the p / d 5/ shell gap, resulting in a relatively larger d value as compared with those of its neighbors. In, a significant violation of the IMME for the = quintet was reported 8]. However, with a later measurement of the excitation energy of the lowest T= state in Na, the IMME was revalidated 8]. Our calculated value of d=.7 kev is in agreement with the later experimental measurement see Fig. a)). The calculation presented in Ref. 8] deviates substantially from both measurements. The data point for = shows a deviation to our calculated value. The studies of this quintet were carried out by several collaborations ], and some of the measured data are controversial ]. We note that both, the experimental and our predicted values, are small. The calculated and experimental e values are compatible with zero, and they are in agreement to each other.

6 6 IV. SUMMRY Based on our recently-proposed GIMME ], we have established an isobaric multiplet mass equation that includes highorder terms within a new strategy beyond the st-order perturbation approximation. The explicit expression of d and also e) coefficient which quantifies the deviation of the quadratic form of the original IMME has been derived. The emergence of the nonzero cubic T z -term, and hence the violation of the quadratic IMME, have basic roots. We have found that the chargesymmetry-breaking interaction in nuclear medium, characterized in our theory by the st-order symmetry energy, combined with the core polarization effect primarily induced by Coulomb force, are responsible for the breakdown of the quadratic IMME. The effective charge-symmetry-breaking and charge-independence-breaking interactions were extracted by employing an ab initio method, i.e., the Brueckner theory with bare V8 and V interactions without any phenomenological adjustments. Remarkably, we found that the calculated d values for quartets and quintets follow an oscillation-like behavior throughout the sd- and lower f p-shell regions as a consequence of the shell effect, and the experimental d values extracted from the measured masses agree the oscillation pattern qualitatively. If all excess neutrons in the T z = T nucleus of a multiplet fill a level above a large shell gap, the breakdown of the quadratic IMME tends to be strong. Therefore, it is straightforward to predict the nuclei where the strong deviations from the original IMME are expected, which is essential for guiding future experimental efforts. cknowledgements This work was supported by the National Natural Science Foundation of China under Grants No , No. 5, No. 5, No , and No. 575, by the 7 Program of China under Grants No. CB8 and No. CB85, by the National Key Program for S&T Research and Development No. 6YF5, 6YF5), by the Youth Innovation Promotion ssociation of Chinese cademy of Sciences, by the Helmholtz-CS Joint Research Group HCJRG-8, and by the European Research Council ERC) under the European Union s Horizon research and innovation program Grant agreement No. 688 STRUm ). Y.H.Z. acknowledges support by the ExtreMe Matter Institute EMMI at the GSI Helmholtzzentrum fĺźr Schwerionenforschung, Darmstadt, Germany. PPENDIX The st-order symmetry energy a CSB) sym,, T z) and the nd-order one a CIB) sym,, T z) for finite nuclei in Eq. ) originate from the CSB and CIB interactions in nuclear medium, respectively, which can be expressed as ] a CSB) sym,, T z) = I a CIB) sym,, T z) = I πr ρr)s CSB) ρ)βr)dr, ) πr ρr)s CIB) ρ)β r)dr, ) where I= N Z)/=T z / is the isospin asymmetry of a given nucleus. βr)=ρ n r) ρ p r))/ρr) is the local isospin asymmetry in whichρ p r) andρ n r) are the proton and neutron density distributions inside the nucleus. Here the symmetry energy coefficient is called simply as the symmetry energy. The spherical-nuclei approximation is employed to achieve a concise result. To achieve the a CSB) sym,, T z) and a CIB) sym,, T z), we should gain the nucleonic density distributionsρ n andρ p of an isobaric analog state IS) whose T is larger than T z. We assume αtt z is the eigenstate of the charge-independent Hamiltonian H, withαfor all additional quantum numbers to label this state. In the first-order perturbation, the energy produced by the CSB and CIB interactions is given by αtt z H CSB+CIB αtt z =a CSB) sym,, T z)t z + acib) sym,, T z)tz in the absence of Coulomb force ]. In this case, the wave function of the IS with T z = T ) with N neutrons and Z+ protons N> Z) is obtained with IS = T, T z = T = T T ] rigidly, where T = t j) is the isospin lowering operator, j excess neutron orbits in. is the ground state of the parent nucleus with N neutrons and Z protons belonging to a multiplet with T= T z. Thus, ρ n +ρ p ) IS = ρ n +ρ p ) parent and ) IS = T )ρ n ρ p ) parent are obtained. The same situation also applies to the IS with T z = T ). ccordingly,

7 the a CSB) sym,, T z) and also a CIB) sym,, T z)) is identical for all members of an isobaric multiplet. Therefore, they can be merged into the b c and c c, respectively. Namely, the GIMME is reduced to the quadratic form of the IMME, and hence the IMME is not breakdown. Here we treat this energy produced by CSB and CIB interactions beyond the first-order perturbation approximation, i.e., calculate the αt z H CSB+CIB αt z instead of αtt z H CSB+CIB αtt z, with the inclusion of Coulomb interaction. ccordingly, the a CSB) sym,, T z) and also a CIB) sym,, T z)) is no longer a constant for a given isobaric multiplet. Yet, we can expand them as a CSB) sym,, T z)t z = a + b T z + c Tz+ d Tz and a CIB) sym,, T z)tz= a + b T z + c Tz+ d Tz for T= / isobaric quartets, and thus Eq. ) is written as ME, T, T z ) = a+b c + nh ) T z + ) a + b T z + c Tz+ d Tz 7 = +c c Tz+ a+a + a + c c + c + c ) a + b T z + c Tz+ d Tz, ) + b c + nh + b + b ) T z ) Tz+ d + d ) T z. ) Therefore, if we introduce a new a coefficient to replace the a+a + a, and define the b, c, d coefficients by then the GIMME for quartets is rewritten as b = b c + nh + b + b, c = c c + c + c, d = d + d ) ME, T, T z )=a+bt z + ct z + dt z. ) The central task is to calculate the d coefficient. It should be stressed that, because of the core polarization induced by the Coulomb force, the state of the T / T is not the exact description of the physical analog state T z = T ). The excess neutron densityρ exc. n,parent, instead ofρ n ρ p should be used in the transition density, as discussed in Ref. 5]. Therefore, one obtains ) Tz=T = ) ) T Tz=T + T δρ T z =T, 5) whereδρr)=ρ core n r) ρ p r) is the neutron- and proton-density difference in the core with an equal number of protons and neutrons) in the T z = T nucleus. The same situation also applies to the IS with T z = T ). In order to improve the accuracy, we use the ) Tz=T = ) ) T Tz=T + ) δρtz =T+δρ Tz =T. 6) T In Fig., we give two examples for T = / doublet to show the validity of such a treatment. Eq. 6) works rather well, and is much better than Eq. 5). In a word, for the T= / isobaric quartets, the density difference is summarized as ), T T z z= =/ ) = + ) T z =/ δρ/ +δρ /, Tz = ) + ) 7) T z = / δρ / +δρ /, Tz = ), T T z z= = / withδρ Tz = ρ core n r) ρ p r) ] = ρ n r) ρ p r) ] core. T z T z

8 8 n - p ) Tz=-/ - fm - ) =7, T=/ Eq. 5) Eq. 6) Direct Cal. 5 6 r fm) =, T=/ 5 6 r fm) FIG. : The density difference between the neutron and proton, i.e.,, for =7 and = isobaric doublets as examples. The Eqs. 5, 6) are employed to calculate the of T z = / nuclei in the framework of the Skyrme-Hartree-Fock method, and compared with the direct calculation with the Skyrme-Hartree-Fock method. We assume thatδρ /,δρ /,δρ /,δρ / are equidistant since the proton number is equidistantly increasing from T z = / to T z = / nuclei, and hence ) ) δρ / δρ / = δρ/ δρ /. Since we expand the first-order and the second-order symmetry energy produced by the CSB and CIB interactions as a CSB) sym,, T z)t z = a + b T z + c Tz+ d Tz and a CIB) sym,, T z)tz= a + b T z + c Tz + d Tz, the d and d are given as d = = = 6 a CSB) sym,, )+acsb) sym,, ) acsb) sym,, ] ) acsb) sym,, ) πr S CSB) ρ) )T z =/ dr πr S CSB) ρ) )T dr z = / πr S CSB) ) ρ) + ) ] δρ/ +δρ / dr T z =/ + πr S CSB) ) ρ) + ) ] δρ / +δρ / dr T z = / = { drπr S CSB) 6 ρ) )Tz=/ ) T z = / ) + ) ] ) δρ/ +δρ / + T z =/ + ) ] } δρ / +δρ / T z = / drπr S CSB) ρ) ) )] δρ / δρ / + δρ/ δρ / = 6 = 8π drπr S CSB) δρ/ ) ρ) δρ / + δρ/ δρ / ) ] drr S CSB) ρ) δρ / δρ / )]. 8)

9 a CIB) sym,, ) acib) sym,, ) acib) sym,, ] )+acib) sym,, ) d = = { 8 ) ) ) + ) πr ρr) S CIB) ρ) ) T z =/ dr πr ρr) S CIB) ρ) ) T z = / dr πr ρr) S ) CIB) ρ) + ) ] δρ/ +δρ / dr T z =/ πr ρr) S CIB) ) ρ) + ) ] } δρ / +δρ / dr T z = / = drπr 8 ρr) S ) CIB) ρ) { ) T z =/ T z = / ) + ) ] ) δρ/ +δρ / + T z =/ + ) ] } δρ / +δρ / T z = / = drπr 8 ρr) S CIB) 8 ) ρ) { 8 ) T z =/ T z = / δρ/ ) +δρ / + ) ) ] δρ/ +δρ / T z =/ + δρ / ) +δρ / + ) ) ]} δρ / +δρ / T z = / = drπr 8 ρr) S CIB) 8 ) ρ) { 8 ) T z =/ δρ T z = / 7 / δρ /) + 8 ρn ) 5 ρ p T z = / δρ /+ ) δρ / ) 5 T z =/ δρ /+ )]} δρ / = drπr 8 8 ρr) S CIB) ρ) { δρ/ +δρ / ] ρn ρ p )Tz=/ ) ] δρ T z = / 7 / δρ /) + 8 ρn ) 5 ρ p T z = / δρ /+ ) δρ / ) 5 T z =/ δρ /+ )]} δρ / = drπr { 8 ρr) S CIB) ρ) δρ 7 / δρ /) 6 ) ) } δρ/ δρ / δρ/ +δρ / 7 = drπr ρr) S CIB) ρ) δρ / δρ /). ) Finally, the d coefficient for the T = / isobaric quartets takes the form of d = d + d = 8π r S CSB) ρ) δρ / δρ / ) dr 8π r S CIB) ρ) ) δρ ρr) / δρ / dr. )

10 For T= isobaric quintets, similarly, the is summarized as ), T T z = z= ) + T z = δρ +δρ ), T z = = δρ, T z = ) + T z = δρ +δρ ), T z = ) T z= T z =. We expand the symmetry energy a CSB) sym,, T z)t z and a CIB) sym,, T z)tz as acsb) sym,, T z)t z = a + b T z + c Tz + d Tz + e Tz and a CIB) sym,, T z)tz = a + b T z + c Tz + d Tz + e Tz, and thus Eq. ) is written as ME, T, T z ) = a+b c + nh ) T z + ) a + b T z + c Tz+ d Tz+ e Tz = +c c Tz+ a+a + a + c c + c + c ) a + b T z + c Tz+ d Tz+ e Tz, ) + b c + nh + b + b ) T z ) Tz+ d + d ) T z+ e + e ) T z. ) The GIMME for quintets is reduced as ME, T, T z )=a+bt z + ct z+ dt z+ et z. ) with d = d + d, ) e = e + e. ) The d, d, e, e coefficients are given by d = = = a CSB) sym, = = = π, )+acsb) sym,, ) acsb) sym,, ) acsb) sym,, )] πr S CSB) ρ) )T dr πr S CSB) z = ρ) )T dr z = ) πr S CSB) ρ) + ] T z = δρ +δρ ) dr ) πr S CSB) ρ) + ] T z = δρ +δρ ) dr { drπr S CSB) ρ) )Tz= ) ) + ] T z = δρ +δρ ) + T z = drπr S CSB) ρ) ] δρ +δρ δρ δρ drπr S CSB) ρ) δρ δρ + ] δρ δρ ) ) + ]} T z = δρ +δρ ) r S CSB) ρ) δρ δρ ] dr. 5)

11 e = a CSB), ) acsb) sym, sym,, ) acsb) sym,, )+acsb) sym,, )+ acsb) sym,, T z)t z Tz = = πr S CSB) ρ) )T dr+ πr S CSB) z = ρ) )T dr z = ) πr S CSB) ρ) + ] T z = δρ +δρ ) dr ) πr S CSB) ρ) + ] T z = δρ +δρ ) dr+ T z Tz = πr S CSB) T ρ)δρ dr z = { drπr S CSB) ρ) )Tz= + ) T z = ) ] + δρ T z +δρ ) ) } + δρ = T z +δρ ) ]+6δρ = { ) ] = π drr S CSB) ρ) + δρ ) } T z +δρ ) + + δρ = T z +δρ ) ] 6δρ = { = π drr S CSB) ρ) )Tz= + ) ]+ δρ +δρ )+δρ +δρ ) ] } 6δρ = π =. T z = drr S CSB) ρ) δρ +δρ )+δρ +δρ )+δρ +δρ ) 6δρ ] 6) 6d = a CIB) sym,, ) acib) sym,, ) acib) sym,, )+ acib) sym,, ) = { πr ρr) S CIB) ρ) ) dr πr T z = ρr) S CIB) ρ) ) dr T z = πr ρr) S CIB) ) ρ) + ] T z = δρ +δρ ) dr + πr ρr) S CIB) ) ρ) + ] T z = δρ +δρ ) dr} = drπr { ρr) S ) CIB) ρ) ) ρn ρ T z = p T z = ) + ] T z = δρ ) +δρ ) + + ] } T z = δρ +δρ ) = drπr { 6 ρr) S ) CIB) ρ) ) ] T z = T z = δρ +δρ ) δρ +δρ ) ] )Tz= δρ +δρ )+ ) } δρ T z +δρ ) = = { drπr 6 ρr) S CIB) ρ) } δρ δρ ) δρ +δρ ) δρ δρ ) = { drπr 6 ρr) S CIB) ρ) ) } δρ δρ = drπr ρr) S CIB) ρ) δρ ) δρ. 7)

12 6e = a CIB) sym,, )+ acib) sym,, ) acib) sym,, ) acib) sym,, )+ acib) sym,, T z)tz T z = = { πr ρr) S CIB) ρ) ) dr+ πr T z = ρr) S CIB) ρ) ) dr T z = πr ρr) S ) CIB) ρ) + ] T z = δρ +δρ ) dr πr ρr) S CIB) ) ρ) + ] T z = δρ +δρ ) dr+ T Tz z T z = πr } ρr) S CIB) ρ)δρ dr = drπr { 6 ρr) S ) CIB) ρ) + ) T z = T z = ρn ) ρ p + ] ρn T z = δρ ) +δρ ) ρ p + ] } T z = δρ +δρ ) + 6δρ = drπr { 6 ρr) S CIB) ρ) δρ +δρ ) + ) δρ T z +δρ ) = + δρ +δρ ) + } )T δρ z +δρ ) 6δρ = = drπr { 6 ρr) S CIB) ρ) δρ δρ ) δρ + ) δρ δρ ) δρ T z = + δρ + δρ ) δρ + ) δρ + δρ ) } δρ δρ T z = = drπr { 5 6 ρr) S CIB) ρ) 8 δρ δρ + 7 ) δρ δρ + ) δρ +δρ )+ ] δρ δρ ) + ) δρ +δρ ) ] } δρ δρ ) δρ = 6 T z = T z = { drπr ρr) S CIB) ρ) 8 δρ δρ ) + δρ δρ ) )Tz= ) Finally, we obtain the d, e coefficients for the T= quintets taking the form of whereρ exc n d = d + d = π e = e + e = π 6 r S CSB) ρ) δρ δρ ) dr π dr r ρr) S CIB) ρ) δρ δρ ) δρ δρ )+8 ρ exc n T z = ]}. 8) r ρr) S CIB) ρ) δρ δρ ) dr, ) T z =+ρ exc p T z = )], ) Tz = ρ exc p Tz = ) is the density of the N Z excess neutrons protons) in the T z = T T z = T) nucleus. ] W. E. Heisenberg, Z. Phys. 77, ). ] E. P. Wigner, in Proc. of the R.. Welch Foundation Conf. on Chemical Research, Houston, edited by W. O. Millikan R.. Welch Foundation, Houston, 57), Vol.. ] S. Weinberg and S. B. Treiman, Phys. Rev. 6, 65 5). ] W. Benenson and E. Kashy, Rev. Mod. Phys. 5, 57 7). 5] W. J. Huang, G. udi, M. Wang, F. G. Kondev, S. Naimi, and X. Xu, Chin. Phys. C, 7). 6] E. Kashy, W. Benenson, D. Mueller, and R. G. H. Robertson, Phys. Rev. C, 5 75). 7] Y. H. Zhang et al., Phys. Rev. Lett., 5 ).

13 8]. T. Gallant et al., Phys. Rev. Lett., 85 ). ] Y. H. Lam et al., t. Data Nucl. Data Tables, 68 ). ] M.. Bentley, S. M. Lenzi, Prog. Part. Nucl. Phys. 5, 7 7). ] M. MacCormick, G. udi, Nucl. Phys. 5, 6 ). ] M. Brodeur, T. Brunner, S. Ettenauer,. Lapierre, R. Ringle, B.. Brown, D. Lunney, and J. Dilling, Phys. Rev. Lett. 8, 5 ). ]. Signoracci and B.. Brown, Phys. Rev. C 8, R) ). ] E. M. Henley and C. E. Lacy, Phys. Rev. 8, 8 6). 5] G. Bertsch and S. Kahana, Phys. Lett. B, 7). 6] W. E. Ormand and B.. Brown, Nucl. Phys., 8). 7]. P. Zuker, S. M. Lenzi, G. Martĺłnez-Pinedo, and. Poves, Phys. Rev. Lett. 8, 5 ). 8] K. Kaneko, Y. Sun, T. Mizusaki, and S. Tazaki, Phys. Rev. Lett., 755 ). ] Y. H. Lam, N.. Smirnova, and E. Caurier, Phys. Rev. C 87, 5 ). ] K. Kaneko, Y. Sun, T. Mizusaki, and S. Tazaki, Phys. Rev. C 8, R) ). ] J. M. Dong, Y. H. Zhang, W. Zuo, J. Z. Gu, L. J. Wang, and Y. Sun, Phys. Rev. C 7, R) 8). ] N. uerbach and. Lev, Nucl. Phys. 8, 65 7). ] J. Britz,. Pape and M. S. ntony, t. Data Nucl. Data Tables 6, 5 8). ] C. Y. Fu et al., Phys. Rev. C 8, 5 8). 5]. Kankainen et al., Phys. Rev. C, R) 6). 6] M. B. Bennett et al., Phys. Rev. C, 6R) 6). 7] F. Herfurth et al., Phys. Rev. Lett. 87, 5 ). 8] M. C. Pyle et al., Phys. Rev. Lett. 88, 5 ). ] C. Yazidjian et al., Phys. Rev. C 76, 8 7). ] R. Ringle et al., Phys. Rev. C 75, 555 7). ] J. Su et al., Phys. Lett. B 756, 6). ] E. Chabanat, P. Bonche, P. Haensel, J. Meyer, R. Schaeffer, Nucl. Phys. 65, 8). ] J. Dong, W. Zuo, and J. Gu, Phys. Rev. C, 5 5). ] B.. Brown and R. Sherr, Nucl. Phys., 6 7). 5] M. Honma, T. Otsuka, B.. Brown, and T. Mizusaki, Phys. Rev. C 65, 6R) ). 6] M. Honma, T. Otsuka, B.. Brown, and T. Mizusaki, Phys. Rev. C 6, 5 ). 7] E. Epelbaum, H.-W. Hammer, and U.-G. Meißner, Rev. Mod. Phys. 8, 77 ). 8] B. E. Glassman et al., Phys. Rev. C, 5R) 5). ] E. G. delberger et al., Phys. Rev. Lett. 8, ); 8, ). ] K. Blaum et al., Phys. Rev. Lett., 68 ). ] S. Triambak et al., Phys. Rev. C 7, 5 6). ].. Kwiatkowski et al., Phys. Rev. C 8, 5R) ). ]. Kankainen et al., Phys. Rev. C 8, 55R) ). ] N. uerbach, Phys. Rep. 8, 7 8). 5] N. uerbach, N. V. Giai, Phys. Rev. C, 78 8).