Contents lists available atScienceDirect
Physics
Letters
B
www.elsevier.com/locate/physletb
Chiral
geometry
and
rotational
structure
for
130
Cs
in
the
projected
shell
model
F.Q. Chen (
)
a,
J. Meng (
)
a,
b,
c,
∗
,
S.Q. Zhang (
)
aaStateKeyLaboratoryofNuclearPhysicsandTechnology,SchoolofPhysics,PekingUniversity,Beijing100871,People’sRepublicofChina bYukawaInstituteforTheoreticalPhysics,KyotoUniversity,Kyoto606-8502,Japan
cDepartmentofPhysics,UniversityofStellenbosch,Stellenbosch,SouthAfrica
a
r
t
i
c
l
e
i
n
f
o
a
b
s
t
r
a
c
t
Articlehistory:
Received18June2018
Receivedinrevisedform30July2018 Accepted22August2018
Availableonline24August2018 Editor: W.Haxton
Keywords:
Nuclearchirality Projectedshellmodel Chiralgeometry
K plot Azimuthal plot
The projectedshellmodel withconfiguration mixingfornuclearchiralityisdeveloped andapplied to theobservedrotationalbandsinthechiralnucleus130Cs.Forthechiralbands,theenergyspectraand
electromagnetictransitionprobabilitiesarewellreproduced.Thechiralgeometryillustratedinthe
K plot
andtheazimuthal plot is
confirmedtoberobustagainsttheconfigurationmixing.Theotherrotational bandsarealsodescribedinthesameframework.©2018TheAuthor(s).PublishedbyElsevierB.V.ThisisanopenaccessarticleundertheCCBYlicense (http://creativecommons.org/licenses/by/4.0/).FundedbySCOAP3.
Spontaneouschiralsymmetrybreakinginatomicnucleihas at-tractedintensivetheoreticalandexperimentalstudiessinceitsfirst predictionbyFrauendorfandMeng[1] in1997.Ittakesplacewhen the angular momentum of a triaxial nucleus has non-vanishing components on all of the three intrinsic principal axes. These three components, contributed by the valenceproton(s), the va-lence neutron(s), and the rotating core, respectively, can form a left- oraright-handedconfiguration.
The picture of chiral geometry is described in the intrinsic frame, and it is manifested in the laboratory system by the ob-servation of the chiral doublet bands, which are a pair of near-degenerate
I
=
1 bands with the same parity. Intensive efforts havebeendevotedtothesearchofchiraldoubletbandsinvarious massregions.Sofararound50pairsofchiraldoubletshavebeen discoveredinA∼
80 [2,3],100 [4–9],130 [10–16] and190 [17,18] mass regions, see Refs. [19–21] for reviews. The chiral doublet bandsobservedexperimentallyuptonowhasbeencompiledvery recentlyinRef. [22].Theoreticalstudiesofthechiraldoubletbands havebeencarriedoutbytheparticlerotormodel(PRM)[1,23–26], thetilted axiscrankingmodel(TAC) [1,27–31], theTAC plusran-*
Correspondingauthorat:SchoolofPhysics,PekingUniversity,Beijing100871, People’sRepublicofChina.E-mailaddresses:mengj@pku.edu.cn(J. Meng),sqzhang@pku.edu.cn
(S.Q. Zhang).
domphaseapproximation[32],thecollectiveHamiltonian[33,34], the interactingboson–fermion–fermion model(IBFFM) [12,35,36], andthegeneralizedcoherentstatemodel[37].
Theangularmomentum projection(AMP)approach,which re-storesthebrokenrotationalsymmetryinthemeanfield wavefunc-tions, is a promising tool for the microscopic description to the nuclearsystem. TheAMPapproach couldbe basedonmeanfield descriptions fromNilsson+BCS [38] tovarious densityfunctionals [39–43],andhasbeenappliedtovariousproblemsassummarized inarecentreview[44].AsoneoftheimplementationsoftheAMP approach,theprojectedshellmodel[38] wasused inattemptsto understand the chiral doublet bands [45,46] in Cs isotopic chain by Bhat et al. The observed energy spectra and electromagnetic transitionsforthedoubletbandsarewellreproducedinthese cal-culations.However,itwasfounddifficulttogiveanillustrationfor theunderlying chiralgeometry.The difficultyliesinthe factthat theangularmomentumgeometryisdefinedintheintrinsicframe, whiletheangularmomentumprojectedwavefunctionsarewritten inthelaboratoryframe.Thisdifficultyisovercomerecentlybythe introductionoftheK plot and theazimuthal plot inRef. [47].The
K plot isdefinedastheprobabilitydistributionsofthecomponents oftheangularmomentumonthethreeaxesoftheintrinsicframe, whiletheazimuthal plot isdefinedastheprobabilitydistributions ofthe polarandtheazimuthal anglesof theangularmomentum inthe intrinsicframe. The K plot andthe azimuthal plot provide
probabilitydistributionsoftheorientationrelatedquantitiesinthe https://doi.org/10.1016/j.physletb.2018.08.039
0370-2693/©2018TheAuthor(s).PublishedbyElsevierB.V.ThisisanopenaccessarticleundertheCCBYlicense(http://creativecommons.org/licenses/by/4.0/).Fundedby SCOAP3.
intrinsic frame, as demonstrated by the chiral doublet bands in 128Cs [47]. With the help of the K plot and the azimuthal plot, theintrinsicorientationofthe angularmomentumcanbe shown explicitly. Consequently, the existence of aplanar orientation for chiralbandscouldbeconfirmedwithoutambiguity.Intheangular momentumprojectedframework,thechiralgeometrycanalsobe illustratedbytheroot-mean-squarevaluesoftheangular momen-tumcomponentsintheintrinsicframe[48].
Itisnoted, however,that theprojectedcalculationinRef. [47] islimitedtotheconfigurationofthelowest
π
h11/2 andthefourthν
h11/2orbitalsonly,andtheconfigurationmixingisneglected.For a betterdescription tothe nuclearsystem, the theoretical frame-workinRef. [47] needstobe generalizedtoincludevarious con-figurations, i.e., to the projected shell model. In order to have a betterunderstandingofthechiraldoubletbands,itisnecessaryto confirm whetherthe chiralgeometryis robust against configura-tionmixing. Another advantageis that theprojected shellmodel canprovidea simultaneousdescriptionforall rotationalbandsin onenucleusonthesamefoot.Inthisworktheprojectedshellmodelwithconfiguration mix-ingfornuclearchiralityisdevelopedandappliedto130Cs[10].The chiralgeometryforthechiraldoubletbandsisillustratedinterms ofthe K plot andtheazimuthal plot,inasimilarwayasRef. [47]. Theeffectoftheconfigurationmixingonthechiralgeometryand thedescriptionoftheotherrotationalbandsarediscussed.
The framework of the projected shell model is based on the standardpairingplusquadrupoleHamiltonian[49],
ˆ
H= ˆ
H0−
χ
2 μˆ
Qμ+Qˆ
μ−
GMPˆ
+Pˆ
−
GQ μˆ
Pμ+Pˆ
μ,
(1)which includes a sphericalsingle-particle partand the separable twobodyinteractions,i.e.,thequadrupole–quadrupoleinteraction, the monopole pairing, and the quadrupole pairing. The intrinsic ground state of an odd–odd nucleus can be denoted as
|
ν0π0,where
ν
0 andπ
0 are single-particle orbitals blocked. The state|
ν0π0 is determined by the variation principle with thecon-straintsonthequadrupolemoments andtheaverageparticle num-bers:
δ
ν0π0
| ˆ
H− λ
q0Qˆ
0− λ
q2Qˆ
2− λ
NNˆ
− λ
ZZˆ
|
ν0π0=
0,
(2)andcanbewrittenasatwo-quasiparticlestate ontopofa quasi-particlevacuumoftheeven-evencore
|
0:|
ν0π0= β
ν+0β
+
π0
|
0.
(3)By the variation in Eq. (2), the state
|
ν0π0 and the vacuum|
0areobtained,aswellascorrespondingquasiparticleoperators{β
+ν
,
β
+π}
,withν
andπ
thesingle-particle orbitals. Varioustwo-quasiparticle states can be constructed as
{|
κ}
= {β
ν+β
π+|
0}
, in whichκ
specifies different two-quasiparticle configurations. The effect of the configuration mixing, which was neglected in Ref. [47],canbethentakenintoaccount.Thetwo-quasiparticle states
|
N,Z,κ withgoodparticlenum-berN and Z canbeprojectedfrom
|
κ:|
N,Z,κ≡ ˆ
PNPˆ
Z|
κ.
(4)The symmetry restored basis is constructed by the angular mo-mentumprojection:
{ ˆ
PIM K|
N,Z,κ}.
(5)The diagonalization ofthe Hamiltonian in the projected basis (5) leads to theHill-Wheeler equation (Notethat theseprojected basisvectorsarenotorthogonal[38]):
Kκ(
N,Z,κ
| ˆ
HPˆ
K KI|
N,Z,κ−
EIσN,Z,κ
| ˆ
PK KI|
N,Z,κ)
fKIσκ=
0,
(6)inwhich
σ
specifiesdifferenteigenstatesofthesamespin I.By solvingEq. (6) theeigenenergiesEIσ andthewavefunctions|
σI M
=
Kκ
fKIσκP
ˆ
M KI|
N,Z,κ (7)are obtained. With the wave functions (7) the electromagnetic transitionsandotherphysicalquantitiescanbecalculated.
In order to examine the configuration mixing, it is necessary to knowthe weightsof differenttwo-quasiparticle configurations
|
κinthewavefunction(7).Theseweightscouldbeobtainedbyresorting the concept ofcollective wave functions inthe genera-tor coordinate method(GCM), andregarding
{
K,
κ
}
inEq. (7) as generatorcoordinates.Thecorresponding generatingfunctionsare theprojectedbasisinEq. (5),andthecorrespondingnormmatrix elementswriteN
I(K,
κ
;
K,
κ
)
≡
N,Z,κ| ˆ
PIK K|
N,Z,κ.
(8)The collective wave functions gIσ
(
K,
κ
)
can be obtained by thesquarerootofthenormmatrix(8) andthefunction fKIσκ :
gIσ
(
K,κ
)
=
K,κN
1/2I
(K
,
κ
;
K,
κ
)
fKIσκ,
(9)and they are proved to satisfy the normalizationandorthogonal relation
KκgIσ1∗
(
K,κ
)g
Iσ2(K,
κ
)
= δ
σ1σ2
.
(10)IntheGCM,thecollectivewavefunctionsareunderstoodas prob-ability amplitudes of the generator coordinates. Therefore, the weightoftheconfiguration
|
κinthestate|
σI Mcouldbewrit-tenas:
Wκ
=
K
|
gIσ(K
,
κ
)
|
2.
(11)Thecollectivewave functions(9) arealsousedinthecalculations ofthe K plot andtheazimuthal plot,whichhavebeenintroduced inRef. [47].
In the following, the five rotational bands observed in 130Cs [50–52],includingthechiraldoubletbands,areinvestigatedbythe presentapproach.TheparametersintheHamiltonian(1) aretaken fromRef. [53].The quadrupoledeformationparameters(
β,
γ
) are constrainedtobe(0.
20,
30◦). Thischoice of(β,
γ
)
isthesameas thatadoptedinthecalculationfor128CsinRef. [47].Itagrees rea-sonably with thedeformation parameters (0.
19,
39◦) usedin the projectedshellmodelcalculationinRef. [53].Different from the calculation in Ref. [47], where the lowest
π
h11/2 andthe fourthν
h11/2 orbitalsare blocked,herethe two-quasiparticle configurationsβ
ν+β
π+|
0 are constructed using or-bitalsν
andπ
fromthe N=
4 andN=
5 majorshells.The quasi-particle energycutoff Eν+
Eπ≤
3.
5MeV isadoptedforthe two-quasiparticleconfigurations,inwhichEν and Eπ arethe quasipar-ticle energies ofthe orbitalsν
andπ
,respectively. There are 60 two-quasiparticle configurations withpositiveparity and68with negative paritytakenintoaccountintheconfigurationspace.Half of themcan be projected ontoeven K valuesandthe other half canbeprojectedontooddK values[54].Thereforethedimension oftheprojected basis(5) forspin I is30(
2I+
1)
forthe positive paritysubspaceand34(
2I+
1)
forthenegativeparitysubspace.Fig. 1. (Coloronline)(Upperpanels)Experimentalandcalculatedrotationalbands in130Cs.Panel(a)showsbandswithpositiveparitywhilepanel(b)showsthose
withnegativeparity.ExperimentaldataaretakenfromRefs. [50–52].Thecalculated energylevelsareshiftedbytakingthelevel9+ inbandAasareference.(Lower panels)ExperimentalandcalculatedvaluesofS(I) = [E(I) −E(I−1)]/2I forbands A,B(c)andbandsC,D(d).
Thecalculatedenergyspectraoftherotationalbands(denoted asbandsA–E)in130CsareshowninFigs.1(a)and(b),in compar-isonwiththeavailabledata[50–52].
As shown in Fig. 1 (a), the calculated yrast and yrare bands withpositive parityagreewell withtheexperimental chiral dou-blet bands, band A and band B, including the near degeneracy between the partner bands. Note that the configuration mixing, which is absent in the calculation in Ref. [47], is takeninto ac-countin thepresentcalculation. Thepresentresultsdemonstrate thattheneardegeneracybetweenthechiraldoubletbandsis per-sistentevenwiththeconfigurationmixing.
Thecalculatedyrastandyrarebandswithnegativeparityagree reasonably with the experimental bands D and C, as shown in Fig.1(b). Thelevels with I
≥
18h in¯
band Dare overestimated, whichmightbeduetoacrossingwithfour-quasiparticle configu-rationsandisbeyondthe presenttwo-quasiparticleconfiguration space.Thecalculatedlowest band withthedominatingconfiguration
ν
g7/2π
h11/2 suggestedin Ref. [51] lies around 1 MeV above the experimentalband E.Thedominatingconfiguration(s) willbe ex-plainedinthe following.Onemaynoticethat themomentof in-ertia for band E is larger than those for bands C and D, which suggestsa largerdeformation forband E.Inthissense,band Eis beyondthe presentscope. Thereforeweexclude bandEfromthe followingdiscussions.Thecompositionofeachstate canbecalculatedfromEq. (11), andthen thedominating configuration(s)can be recognized. The dominating configurations for bands A–D at the band head are given in Table 1, together with those suggested in the previous studies.Thedominatingconfigurationsfoundbythepresent calcu-lationcoincidewiththepreviousassignments,exceptbandC.The dominatingconfigurationfoundforbandC is
ν
h11[5th]/2π
d[2nd]5/2 while thepreviously assignedoneisν
h11/2π
g7/2.Thisdifferenceisdue to the strong mixing between the g7/2 and the d5/2 orbitals as showninFig.2.InFig.2,thecompositionsofconfigurationsforbandsA–D are shownasfunctionsof spin.ForbandsA andB, theconfiguration
ν
h[115th/2]π
h11[1st/2] isdominantuntilI=
16h,¯
afterwhichstrong config-urationmixingoccurs.ForbandC,theconfigurationsν
h[115th/2]π
d[52nd/2]Table 1
ThedominatingconfigurationsforbandsA–Datthebandheadinthepresent cal-culation,incomparisonwiththosesuggestedinthepreviousstudies.Theordinal numbersinthesecondcolumndenotespecificsingle-particleorbitals.For exam-ple,thenotationπh11[1st/]2representsthefirst(lowest)single-particleorbitalinthe
πh11/2subshell.
Bands Present configuration Previous configuration References A, B νh[115th/2]πh[111st/2] νh11/2πh11/2 [10,52,51]
C νh[115th/2]πd5[2nd/2] νh11/2πg7/2 [50,51]
D νh[115th/2]πd
[2nd]
5/2 νh11/2πd5/2 [50,51]
Fig. 2. (Coloronline)CompositionofconfigurationforbandsA–D.Foreachband, contributionsofthedominatingconfigurationsasfunctionsofspinarepresented bycoloredlineswiththeircorrespondingconfigurationlabeledbythesamecolor. Contributionsofotherconfigurations(iflargerthan1%)arepresentedbydark yel-lowlines.
and
ν
h[115th/2]π
g7[2nd/2] stronglycompete witheachother.The configu-rationν
h[115th/2]π
d5[2nd/2] winsatthebandhead,buttheconfigurationν
h[115th/2]π
g7[2nd/2] takesoverat I=
8h.¯
Forband D,theconfigurationν
h[115th/2]π
d5[2nd/2] isdominantuntil I=
13h,¯
afterwhichstrong con-figurationmixingoccurs.The amplitudes of signature splitting of bands A–D, reflected bythequantity S
(
I)
= [
E(
I)
−
E(
I−
1)
]/
2I,areshowninFigs.1(c) and(d).ForbandsAandB,theangularmomentumofthevalence neutronorientatesalongthelongaxis(l-axis),andthatofthe va-lenceprotonorientatesalongtheshortaxis(s-axis).Bothofthem areperpendicular approximatelytotheangularmomentumofthe collective rotation.Thevaluesof S(
I)
stay almostindependent of spin,asisseeninFig.1(c).Infact,thespinindependenceofS(
I)
has beenrecognized as one ofthe criteriafor the chiraldoublet bands[4,55]. ForbandsCandD, the strongconfigurationmixing leadstotheincreasingsignaturesplittingwithspin[56] asshown inFig.1(d).Moreover,it isnotedinFig.1that theamplitudesofS
(
I)
increase withspin, whichisduetotheincreasing configura-tionmixing.In thefollowing we focus ourdiscussion on theproperties of thechiraldoubletbandsAandB.
Thecalculated intrabandB
(
E2)
andB(
M1)
,andtheinterbandB
(
M1)
forbandsAandBareshowninFig.3,incomparisonwith thedata[57].Figure3(a)showsthattheintrabandB(
E2)
are simi-larforbandsAandB,inbothexperimentalandtheoreticalresults, as expectedfor a pair ofchiral doublets [24,55]. The B(
E2)
val-uesaresomewhatoverestimatedbythecalculation.Thesimilarity oftheintrabandB(
M1)
betweenthechiraldoublets,togetherwithFig. 3. (Coloronline)(a)TheintrabandB(E2),(b)theintrabandB(M1),and(c)the interbandB(M1)forthechiraldoubletbandsin130Cscalculatedbytheprojected
shellmodel,comparedwiththeexperimentaldatatakenfromRef. [57].
the staggering behavior of both intra- and interband B
(
M1)
, are found in Figs. 3(b) and(c), in both experimental andtheoretical results.Thesefeaturesarealsoexpectedassignaturesofthechiral modesaccordingtoRefs. [24,55].Thecharacteristicfeaturesofthe electromagnetictransitions expectedforthe chiral doublet bands showninFig.3,whichhavebeenreproducedwithasingle config-uration[47],arepersistentevenwhentheconfigurationmixingis takenintoaccount.TheangularmomentumgeometryofthechiraldoubletsAand B can be illustrated by the K plot and the azimuthal plot, as in Ref. [47].TheK plot,i.e.,theK distributionspIσ
(
|
K|)
forthe angu-larmomentumonthethreeprincipalaxes,areshowninFig.4at thespins I=
10,
11,
14,
15,
18,
19h.¯
Theseresults couldbe com-paredwiththecorresponding resultsobtainedbythePRM,which havebeendiscussedinRefs. [26,58].AsseenintheK plot inFig.4,theevolutionofthechiralmodes canbeexhibited.
Forspins I
=
10 and11h,¯
the K plots shownin Fig.4are in accordancewiththe expectationforchiral vibrationwithrespect to the s–l plane,which takesplace nearthe band headwith in-sufficientcollectiverotation.TheprobabilityatKi=
0 issignificantforband A, which indicates a wave function symmetric with
re-specttoKi
=
0 andcorrespondstoa0-phononstate.Ontheotherhand,the vanishingprobability at Ki
=
0 forband B indicatesanantisymmetricwavefunctioncorrespondingtoa1-phononstate. Forspins I
=
14 and15¯
h,themostprobablevalue forKiap-pears atKi
∼
13h for¯
bandA, whichsuggeststhatthe collectiverotation aroundthei-axisdevelopsandtheangularmomentafor band Adeviate fromthes–l plane withtheincrease ofspin. The peaks of the Kl- and Ks-distributions also correspond to K
val-uessubstantiallyawayfromzero,sothemostprobableorientation oftheangularmomentumisaplanar,indicatingtheoccurrenceof staticchirality.TheK plots forbandBaresimilartothoseforband A, and the occurrence of static chirality is thus supported. The tunneling betweenthe left-handedand the right-handed config-urations isresponsiblefortheenergyseparationbetweenthetwo bands,whichexplainswhythedoublets areclosestto eachother inenergyatI
=
14 and15h.¯
For spins I
=
18 and 19 h,¯
the Ks- and Kl-distributionsbe-come broad,andthe peaksofthe Ki-distributions become sharp
at Ki
∼
I. Both features suggest that the angular momentamoveclosetothei-axis.Thereforethestaticchiralitydisappearsandthe aplanarrotationisreplacedbytheprincipalaxisrotation.
InFig.5,theprofilesfortheorientationoftheangular momen-tumonthe
(θ,
φ)
plane,theazimuthal plot,areshownforthesame spinsasinFig.4.Thedefinitionsoftheangles(θ,
φ)
canbefound inRef. [47].ForspinsI
=
10 and11h,¯
theprofilesfortheorientationofthe angularmomentumforbandAhaveasinglepeakat(θ
∼
60◦,
φ
=
90◦)
, which suggests that the angular momentum stays within the s–l plane,inaccordancewiththeexpectationfora 0-phonon state. On the other hand, the profiles for band B show a node at(θ
∼
60◦,
φ
=
90◦)
, withtwo peaks at(θ
∼
65◦,
φ
∼
45◦)
and(θ
∼
65◦,
φ
∼
135◦)
, respectively. The existence of the node and thetwopeakssupportstheinterpretationofa1-phononvibration. Thereforetheinterpretationofchiralvibrationisdemonstrated.ForspinsI
=
14 and 15h,¯
theazimuthal plots forbandsAand Baresimilar.Twopeakscorrespondingtoaplanarorientationsare found,i.e.(θ
∼
75◦,
φ
∼
35◦)
and(θ
∼
75◦,
φ
∼
145◦)
forband A, while(θ
∼
70◦,
φ
∼
40◦)
and(θ
∼
70◦,
φ
∼
140◦)
forbandB.These featurescouldbeunderstoodasarealizationofstaticchirality.The non-vanishing distribution forθ
=
90◦ andφ
=
90◦ reflects the tunnelingbetweentheleft- andright-handedconfigurations.For spins I
=
18 and 19 h,¯
the peaks for the azimuthal plotfor band A move toward (
θ
∼
80◦,
φ
∼
20◦) and (θ
∼
80◦,
φ
∼
160◦), namely closeto the i-axis.This isin accordance withtheFig. 4. (Color online) The K plot,i.e., the K distributions for the angular momentum on the short (s),intermediate (i) and long (l)axes, calculatedat spins I=
Fig. 5. (Coloronline)Theazimuthal plot, i.e.,profilefor the orientationofthe angularmomentum onthe (θ, φ)plane,calculatedat spins I=10, 11, 14, 15, 18, 19¯h,
respectively,forthechiraldoubletbandsin130Cs.
dominance of rotation around the i-axis reflected in the K plot
discussed above, suggesting the disappearance of chiral geome-try and the onset of principal axis rotation. The peaks for the
azimuthal plot for band B locate at
(θ
∼
75◦,
φ
∼
35◦)
and(θ
∼
75◦,
φ
∼
145◦)
,whicharesimilartothose atI=
14,
15h but¯
ap-proaching the i-axis. Thus the chiral geometry is weakened. In general,theevolutionoftheangularmomentumgeometryintheazimuthal plot isconsistentwiththoseinthe K plot.
Byblockingthe lowest
π
h11/2 and thefourthν
h11/2 orbitals, the chiral geometry in K plot and azimuthal plot has been dis-cussedinRef. [47].Figs.4and5
confirmthatthechiralgeometry in K plot and azimuthal plot canbe robust against configuration mixing.Insummary,theprojectedshellmodelwithconfiguration mix-ing for nuclear chirality is developed and applied to the ob-servedrotationalbandsinthechiralnucleus130Cs.Boththechiral andachiral bandsare described on the samefoot. Forthe chiral bands, the energy spectra, S
(
I)
, B(
M1)
and B(
E2)
are well re-produced.Thechiralgeometryisdemonstratedinthe K plot andtheazimuthal plot. The robustness ofthe chiral geometryagainst theconfiguration mixingisconfirmed. As another advantage,the otherrotationalbandsaredescribedsimultaneouslywiththesame Hamiltonian.
Acknowledgements
Thiswork waspartlysupported bytheNationalKeyR&D Pro-gramofChina No.2018YFA0404400, theNationalNaturalScience FoundationofChina(GrantsNo.11335002andNo.11621131001), and the China Postdoctoral Science Foundation under Grant No. 2017M610688.
References
[1]S.Frauendorf,J.Meng,Nucl.Phys.A617(1997)131.
[2]S.Y.Wang,B.Qi,L.Liu,S.Q.Zhang,H.Hua,X.Q.Li,Y.Y.Chen,L.H.Zhu,J.Meng, S.M.Wyngaardt,etal.,Phys.Lett.B703(2011)40.
[3]C.Liu,S.Y.Wang,R.A.Bark,S.Q.Zhang,J.Meng,B.Qi,P.Jones,S.M.Wyngaardt, J.Zhao,C.Xu,etal.,Phys.Rev.Lett.116(2016)112501.
[4]C.Vaman,D.B.Fossan,T.Koike,K.Starosta,I.Y.Lee,A.O.Macchiavelli,Phys.Rev. Lett.92(2004)032501.
[5]P.Joshi,M.P.Carpenter,D.B.Fossan,T.Koike,E.S.Paul,G.Rainovski,K.Starosta, C.Vaman,R.Wadsworth,Phys.Rev.Lett.98(2007)102501.
[6]D.Tonev,M.S.Yavahchova,N.Goutev,G.deAngelis,P.Petkov,R.K.Bhowmik, R.P.Singh,S. Muralithar,N.Madhavan,R.Kumar,et al.,Phys.Rev.Lett.112 (2014)052501.
[7]E.O.Lieder,R.M.Lieder,R.A.Bark,Q.B.Chen,S.Q.Zhang,J.Meng,E.A.Lawrie, J.J.Lawrie,S.P.Bvumbi,N.Y.Kheswa,etal.,Phys.Rev.Lett.112(2014)202502.
[8]N.Rather,P.Datta,S.Chattopadhyay,S.Rajbanshi,A.Goswami,G.H.Bhat,J.A. Sheikh,S.Roy,R.Palit,S.Pal,etal.,Phys.Rev.Lett.112(2014)202503.
[9]I.Kuti,Q.B.Chen,J.Timár,D.Sohler,S.Q.Zhang,Z.H.Zhang,P.W.Zhao,J.Meng, K.Starosta,T.Koike,etal.,Phys.Rev.Lett.113(2014)032501.
[10]K.Starosta,T.Koike,C.J.Chiara,D.B.Fossan,D.R.LaFosse,A.A.Hecht,C.W. Beau-sang,M.A.Caprio,J.R.Cooper,R.Krücken,etal.,Phys.Rev.Lett.86(2001)971.
[11]S.Zhu,U.Garg,B.K.Nayak,S.S.Ghugre,N.S.Pattabiraman,D.B.Fossan,T.Koike, K.Starosta,C.Vaman,R.V.F.Janssens,etal.,Phys.Rev.Lett.91(2003)132501.
[12]D.Tonev,G.deAngelis,P.Petkov,A.Dewald,S.Brant,S.Frauendorf,D.L. Bala-banski,P.Pejovic,D.Bazzacco,P.Bednarczyk,etal.,Phys.Rev.Lett.96(2006) 052501.
[13]C.M. Petrache,G.B.Hagemann,I.Hamamoto,K. Starosta,Phys. Rev.Lett.96 (2006)112502.
[14]E.Grodner,J.Srebrny,A.A.Pasternak,I.Zalewska,T.Morek,C.Droste,J. Mierze-jewski, M. Kowalczyk,J. Kownacki,M. Kisieli ´nski,et al., Phys.Rev.Lett. 97 (2006)172501.
[15]S.Mukhopadhyay,D.Almehed,U.Garg,S.Frauendorf,T.Li,P.V.M.Rao,X.Wang, S.S.Ghugre,M.P.Carpenter,S.Gros,etal.,Phys.Rev.Lett.99(2007)172501.
[16]A.D.Ayangeakaa,U.Garg,M.D.Anthony,S.Frauendorf,J.T.Matta,B.K.Nayak, D.Patel,Q.B.Chen,S.Q.Zhang,P.W.Zhao,etal.,Phys.Rev.Lett.110(2013) 172504.
[17]D.L.Balabanski,M.Danchev,D.J.Hartley,L.L.Riedinger,O.Zeidan,J.Y.Zhang, C.J.Barton,C.W.Beausang,M.A.Caprio,R.F.Casten,etal.,Phys.Rev.C70(2004) 044305.
[18]E.A. Lawrie,P.A.Vymers,J.J.Lawrie,C.Vieu,R.A.Bark,R.Lindsay, G.K. Ma-bala, S.M.Maliage,P.L.Masiteng,S.M.Mullins, etal.,Phys.Rev.C78(2008), 021305(R).
[19]J.Meng,S.Q.Zhang,J.Phys.G,Nucl.Part.Phys.37(2010)064025.
[20]J.Meng,Q.B.Chen,S.Q.Zhang,Int.J.Mod.Phys.E23(2014)1430016.
[21]J.Meng,P.W.Zhao,Phys.Scr.91(2016)053008.
[22]B.W.Xiong,Y.Y.Wang,arXiv:1804.04437,2018.
[23]J.Peng,J.Meng,S.Q.Zhang,Phys.Rev.C68(2003)044324.
[24]T.Koike,K.Starosta,I.Hamamoto,Phys.Rev.Lett.93(2004)172502.
[25]S.Q.Zhang,B.Qi,S.Y.Wang,J.Meng,Phys.Rev.C75(2007)044307.
[26]B.Qi,S.Q.Zhang,J.Meng,S.Y.Wang,S.Frauendorf,Phys.Lett.B675(2009) 175.
[27]V.I.Dimitrov,S.Frauendorf,F.Dönau,Phys.Rev.Lett.84(2000)5732.
[28]H. Madokoro, J. Meng, M. Matsuzaki, S. Yamaji, Phys. Rev. C 62 (2000), 061301(R).
[29]P. Olbratowski, J.Dobaczewski, J.Dudek, W.Płóciennik, Phys. Rev.Lett. 93 (2004)052501.
[30]P.Olbratowski,J.Dobaczewski,J.Dudek,Phys.Rev.C73(2006)054308.
[31]P.W.Zhao,Phys.Lett.B773(2017)1.
[32]D.Almehed,F.Dönau,S.Frauendorf,Phys.Rev.C83(2011)054308.
[33]Q.B.Chen,S.Q.Zhang,P.W.Zhao,R.V.Jolos,J.Meng,Phys.Rev.C87(2013) 024314.
[34]Q.B.Chen,S.Q.Zhang,P.W.Zhao,R.V.Jolos,J.Meng,Phys.Rev.C94(2016) 044301.
[35]D.Tonev,G.deAngelis,S.Brant,S.Frauendorf,P.Petkov,A.Dewald,F.Dönau, D.L.Balabanski,Q.Zhong,P.Pejovic,etal.,Phys.Rev.C76(2007)044313.
[36]S.Brant,D.Tonev,G.deAngelis,A.Ventura,Phys.Rev.C78(2008)034301.
[37]A.A.Raduta,A.H.Raduta,C.M.Petrache,J.Phys.G,Nucl.Part.Phys.43(2016) 095107.
[38]K.Hara,Y.Sun,Int.J.Mod.Phys.E4(1995)637.
[39]B.Sabbey,M.Bender,G.F.Bertsch,P.-H.Heenen,Phys.Rev.C75(2007)044305.
[40]T.R.Rodríguez,J.L.Egido,Phys.Rev.C81(2010)064323.
[41]T.R.Rodríguez,A.Arzhanov,G.Martínez-Pinedo,Phys.Rev.C91(2015)044315.
[42]T.Nikši ´c,D.Vretenar,G.A.Lalazissis,P.Ring,Phys.Rev.Lett.99(2007)092502.
[43]J.M.Yao,J.Meng,P.Ring,D.Vretenar,Phys.Rev.C81(2010)044311.
[44]Y.Sun,Phys.Scr.91(2016)043005.
[45]G.H.Bhat,J.A.Sheikh,R.Palit,Phys.Lett.B707(2012)250.
[46]G.H.Bhat,R.N.Ali,J.A.Sheikh,R.Palit,Nucl.Phys.A922(2014)150.
[47]F.Q.Chen,Q.B. Chen,Y.A.Luo,J. Meng,S.Q. Zhang, Phys.Rev.C96 (2017) 051303.
[48]M.Shimada,Y.Fujioka,S.Tagami,Y.R.Shimizu,Phys.Rev.C97(2018)024319.
[49]P.Ring,P.Schuck, TheNuclearManyBody Problem,Springer-Verlag,Berlin, 1980.
[50]P. Sala, N. Blasi, G.L. Bianco, A. Mazzoleni, R. Reinhardt, K. Schiffer, K. Schmittgen,G.Siems,P.V.Brentano,Nucl.Phys.A531(1991)383.
[51]R. Kumar, D. Mehta, N. Singh, H. Kaur, A. Görgen, S. Chmel, R. Singh, S. Murlithar,Eur.Phys.J.A11(2001)5.
[52]A.J.Simons,P.Joshi,D.G.Jenkins,P.M.Raddon,R.Wadsworth,D.B.Fossan,T. Koike,C.Vaman,K. Starosta,E.S.Paul,etal.,J.Phys. G,Nucl.Part.Phys.31 (2005)541.
[53]Z.C.Gao,Y.Chen,Y.Sun,Phys.Lett.B634(2006)195.
[54]J.A.Sheikh,G.H.Bhat,Y.Sun,G.B.Vakil,R.Palit,Phys.Rev.C77(2008)034313.
[55]S.Y.Wang,S.Q.Zhang,B.Qi,J.Meng,Chin.Phys.Lett.24(2007)664.
[56]S.Frauendorf,J.Meng,Z.Phys.A356(1996)263.
[57]L.L.Wang,X.G.Wu,L.H.Zhu,G.S.Li,X.Hao,Y.Zheng,C.Y.He,L.Wang,X.Q.Li, Y.Liu,etal.,Chin.Phys.C33(2009)173.
[58]B. Qi, S.Q. Zhang, S.Y. Wang, J.M. Yao, J. Meng, Phys. Rev. C 79 (2009), 041302(R).