-
The atomic nucleus is a microscopic quantum many-body system, and its shape provides an intuitive understanding of nuclear spatial density distributions. Stable quadrupole-deformed nuclear shapes have been known to exist in the intrinsic frame for a long time [1]. With reflection symmetry breaking, some nuclei appear to have octupole deformation in the intrinsic frame, such as a pear shape [2]. The study of the pear-shaped nuclei has been linked to the frontiers of particle physics, nuclear physics, and atomic physics, since it provides a distinctive probe to test the charge-parity (CP) violation beyond the Standard Model [3].
The odd-mass octupole-deformed nuclei have attracted special interest as they have enhanced nuclear Schiff moments owing to the presence of the large octupole collectivity and nearly degenerate parity doublets [4]. The best candidates of such nuclei have been found in the light-actinide region, such as
223Th [5] and225Ra [6]. For a lighter mass region withA∼ 150, continuous efforts have also been devoted owing to its octupole double-magic character with proton and neutron numbersZ∼56 andN∼88 [2]. Indeed, the Coulomb excitation experiments of the octupole double-magic nucleus144Ba and its neighbor146Ba yield enhancedE3 matrix elements and provide direct evidence of octupole deformation [7, 8]. A further experimental study on level structures shows the coexistence of reflection asymmetric and symmetric shapes in the neutron-rich nucleus144Ba [9]. For odd-mass nuclei in this mass region, although no direct experimental evidence of octupole deformation was found, the characteristics of octupole collectivity such as parity doublet bands have been systematically observed in nuclei such as143,145Ba [10, 11],145,147La [12, 13],147Pr [14],151Pm [15, 16], and153,155Eu [17−20].The
N=90 isotone153Eu is a stable isotope. The excited level structure of153Eu has been extensively studied experimentally by Coulomb excitations [19, 21, 22],154Sm (d,3n ) [23], and150Nd (7Li ,4n ) [17, 18] reactions. The positive-parity band built on theK=5/2+ ground state and the negative-parity band built on theK=5/2− state at 97.5 keV constitute two nearly degenerate bands, which combined with the strongE1 transitions linking them, show the prominent characteristics of parity doublet bands [17−19]. The bands have been extended up to spins as high as39/2+ and43/2− [18]. However, the difference between the orbital and rotational gyromagnetic factors, i.e., thegK−gR values, extracted from theB(M1)/B(E2) ratios [17, 18] and the measured lifetimes [19] were found to be different for the opposite parity states of doublet bands, indicating that the bands are not based on an intrinsic parity-mixed state. As the preponderance ofE1 transitions and large electric-dipole moments seem unable to be interpreted in the frame of reflection-symmetric single-particle models, explaining these experimental observations in theory remains a challenge.A reflection-asymmetric triaxial particle rotor model (RAT-PRM) was recently developed, which included both triaxial and octupole degrees of freedom [24, 25]. The model has been applied to investigate the multiple chiral doublet bands with octupole correlations [24, 26, 27]. By analyzing the RAT-PRM Hamiltonian, a new symmetry for an ideal chirality-parity violation system was found and the corresponding selection rules of the electromagnetic transitions were derived [28]. Moreover, the RAT-PRM has been applied to investigate the energy spectra and electromagnetic properties of the possible parity doublet bands in some odd mass nuclei, such as
73Br ,79Se inA∼ 80 [29, 30],143Ba ,151Pm inA∼ 150 [10, 31], and217Ra ,223Th inA∼ 220 [32, 33] mass regions.This study will investigate the octupole correlations in the stable nucleus
153Eu using RAT-PRM. The model is briefly introduced in Sec. II and the numerical details are presented in Sec. III. The calculated results for the positive- and negative-parity bands, such as the energy spectra and electromagnetic transitions are discussed in Sec. IV. A summary is given in Sec. V. -
The total RAT-PRM Hamiltonian [24] is the sum of the intrinsic Hamiltonian
ˆHp(n)intr. for valence protons (neutrons) and the collective HamiltonianˆHcore of a reflection-asymmetric triaxial rotor. The intrinsic HamiltonianˆHp(n)intr. for valence nucleons is expressed as [34−36],ˆHp(n)intr.=∑ν>0(εp(n)ν−λp(n))(a†νaν+a†ˉνaˉν)−Δp(n)2∑ν>0(a†νa†ˉν+aˉνaν),
(1) with
λp(n) representing the Fermi energy andΔp(n) representing the pairing gap parameter. The single-particle energyεp(n)ν is obtained by diagonalizing a reflection-asymmetric triaxial single-particle HamiltonianˆHp(n)s.p. [24], which has the form of the Nilsson Hamiltonian [37]. The core HamiltonianˆHcore is expressed as,ˆHcore=3∑k=1ˆR2k2Jk+12E(0−)(1−ˆPc),
(2) with
ˆRk representing the angular momentum operators of the core andJk=J0sin2(γ−2kπ/3) representing the moments of inertia (MoIs) of the irrotational flow type. In the second term, the core parity splitting parameterE(0−) is treated as a free parameter, describing the excitation energy of the virtual0− state [38], and the core parity operatorˆPc is the product of the single-particle parity operatorˆπ and the total parity operatorˆP .The total RAT-PRM Hamiltonian is diagonalized numerically in the symmetrized strong coupling basis with good parity and angular momentum, which gives rise to the eigenvalues and eigen wavefunctions [24]. For the electric multipole transition calculations, the corresponding operators contain two terms [1, 24],
ˆM(Eλ,μ)=ˆq(c)λμ+ˆq(p)λμ=3Ze4πRλ0βλμ+eeffn∑i=1rλiY∗λμ,
(3) which consider contributions from core and valence particles. Here,
R0=1.2A1/3 fm is the nuclear radius, andeeff is the effective charge. For the electric quadrupole (E2 ) transition calculations, the valence particle term in Eq. (3) is neglected, since it is much smaller than the core term [39]. For the magnetic dipole transition calculations, the corresponding operator is given asˆM(M1,μ)=√34πeℏ2Mc[(gp(n)−gR)ˆjp(n)1μ],
(4) where
gp andgR are the effective gyromagnetic ratios for the valence proton (neutron) and collective core, respectively, andˆj1μ denotes the spherical tensor in the laboratory frame. -
The deformation parameters in the present RAT-PRM calculation of
153Eu are based on the microscopic covariant density functional theory (CDFT) in a 3D lattice [40] with the successful point coupling density functional PC-PK1 [41]. Figure 1 shows the potential energy surface (PES) in the(β20,β30) plane for153Eu calculated by the constrained CDFT in 3D lattice calculations. The global minimum of the PES is found to be located atβ20=0.30 andβ30=0.00 . All the axial asymmetric deformations, i.e.,β2μ,β3μ,β4μ,(μ≠0) , are found to be zero, which means that the ground state of153Eu corresponds to an axial and reflection symmetric shape. The PES exhibits a relatively soft character around the minimum. With the above adopted deformation parameters, a unified description for the positive- and negative-parity bands and the electromagnetic transitions between them in153Eu can be obtained. In addition, the RAT-PRM calculations withβ30=0.03,0.05 , and0.07 are also performed to investigate the effect of octupole deformation.For the intrinsic part, the reflection asymmetric triaxial Nilsson Hamiltonian with the parameters
κ,μ in Ref. [42] is solved in the harmonic oscillator basis [43]. The proton Fermi energy is chosen asλp=43.13 MeV and the pairing gap is calculated using the empirical formulaΔp=12/√A MeV. The single-particle space available to the odd proton is truncated to 13 levels, with six above and six below the Fermi level. Increasing the size of the single-particle space does not influence the band structure in the calculations.For the core part, it turns out that a spin-dependent MoI, i.e.,
J(I)=(a+bI)ℏ2/ MeV is necessary to reproduce the experimental energy spectra. This spin-dependent MoI can be attributed to the increasing behavior versus spin for the effective MoI extracted from the experimental data of positive- and negative-parity bands in153Eu [19]. The parameter a is 30 for the positive-parity band and 18 for the negative-parity band, while the parameter b is 0.85 for both bands. The core parity splitting parameterE(0−)=2.0 MeV is used. An ad hoc Coriolis attenuation factor ξ is generally needed in PRM description [44], which is 0.7 for positive-parity band and0.8 for the negative-parity band in the present calculation.The intrinsic electric-quardupole moment
Q0=6.75 eb, determined empirically in Ref. [19], is adopted tor the electric-multipole transition calculations, while the intrinsic electric-dipole momentQ10=34πR0Zβ10 is calculated withR0=1.2A1/3 fm. For the magnetic-dipole transition calculations, the gyromagnetic ratios for the collective rotor and valence proton are given bygR=Z/A andgp=gl+(gs−gl)/(2l+1) , respectively [1, 39]. -
The calculated physical observables for the positive- and negative-parity bands in
153Eu , including the excitation energies, the energy staggering parametersS(I)=[E(I)−E(I−1)]/(2I) , and the intrabandE2 andM1 transition probabilities, are shown in Fig. 2 in comparison with the data available [19].Figure 2. (color online) The energies
E(I) (a), the energy staggering parametersS(I)=[E(I)−E(I−1)]/(2I) (b), and the intrabandE2 (c) andM1 (d) transition probabilities for the positive-parity bandπ=+ and negative-parity bandπ=− in 153Eu calculated by RAT-PRM (lines), in comparison with the experimental data [19] (symbols).As shown in Fig. 2(a), the experimental positive- and negative-parity bands are nearly degenerate in energy. The negative-parity band is situated at a slightly higher level in the bandhead vicinity. As the spin increases, it undergoes a downward shift with respect to the positive-parity band, reaching a point of intersection at
I=11/2ℏ . The calculated excitation energies well reproduce the data for both the positive-parity bandπ=+ and negative-parity bandπ=− . For the observed spin range of5/2ℏ≤I≤37/2ℏ , the calculated average energy difference between the positive- and negative-parity bands is 10 keV, which is close to the experimental value of 14 keV. Figure 2(b) shows the behavior of signature splitting for the positive- and negative-parity bands, as represented by the staggering parametersS(I) . In contrast to the nearly vanished splitting observed for the positive-parity band, theS(I) values exhibit pronounced signature splitting for the negative-parity band. The RAT-PRM calculations well reproduce the experimental values and behaviors in the spin regionI≥11/2ℏ , and the differentS(I) behaviors for the two bands are attributed to their different configurations. The configuration for the positive-parity band is found to be dominated byπg7/2[Ω=5/2] , andπh11/2[Ω=5/2] for the negative-parity band. For the positive-parity band, the componentπg7/2[Ω=5/2] with a relatively-high Ω is always predominant as the spin increases, leading to the nearly vanished signature splitting. For the negative-parity band, the componentπh11/2[Ω=5/2] decreases rapidly while the componentπh11/2[Ω=3/2] with a relatively-low Ω increases as the spin increases, which leads to the pronounced signature splitting.As shown in Fig. 2(c), the experimental
B(E2) values for the positive-parity band increase as the spin increases, whereas for the negative-parity band, the first two data follow the increasing tendency but the last datum atI=21/2ℏ has an obvious drop. The calculatedB(E2) values well reproduce the increasing trend and are close to the experimental values in magnitude. More experimental data are necessary to pin down theB(E2) tendency for the negative-parity band. TheM1 transition strengths are sensitive to the single-particle components of the intrinsic wave functions. As shown in Fig. 2(d), the experimentalB(M1) values exhibit a pronounced staggering behavior for the negative-parity band, while this behavior is not shown in the positive-parity band. Although the calculatedB(M1) values differ by a factor of two from the experimental results, the trends of the experimentalB(M1) values are well reproduced for both positive- and negative-parity bands, indicating the proper intrinsic wave functions in the present calculations.The observation of enhanced interband E1 transitions connecting the positive- and negative-parity bands is an important signal of octupole correlations in atomic nuclei. Figures 3(a) and (b) show the calculated interband
B(E1) values from the positive-parity to negative-parity bands (π+→π− ) and those from the negative-parityto positive-parity bands (π−→π+ ), respectively, in comparison with the available experimental data [19]. The calculatedB(E1) values withβ30=0.00 generally underestimate the experimental data. By introducing the octupole deformation in the RAT-PRM calculations, theB(E1) values are found to depend sensitively on the value of the octupole deformation parameterβ30 . The calculatedB(E1) values withβ30=0.07 are two orders of magnitude higher than those withβ30=0.00 . The calculated values show good agreement with the experimental data forβ30=0.05 . This result indicates that the single-particle contributions in Eq. (3) alone are insufficient to account for the substantialE1 transitions observed in153Eu . By contrast, no notable influences on the excitation energies andB(E2) andB(M1) values were observed whenβ30 was changed from 0.00 to 0.07. Furthermore, considering the soft behavior of PES inβ20 , RAT-PRM calculations withβ20=0.35 andβ30 changing from 0.00 to 0.07 have been performed to investigate the sensibility ofB(E1) onβ30 at differentβ20 . Based on the reproduction of the energy spectra, the sensibility ofB(E1) toβ30 atβ20=0.35 showed similar behavior to that atβ20=0.30 .Figure 3. (color online) The calculated
B(E1) values of the interbandE1 transitions in 153Eu in comparison with the available experimental data [19] for (a) the positive-parity to negative-parity bands (π+→π− ) and (b) the negative-parity to positive-parity bands (π−→π+ ). The RAT-PRM calculations withβ30=0.00,0.03,0.05, and0.07 are denoted with different lines.To determine the roles of the two terms in Eq. (3), Figs. 4(a) and (b) show the calculated
B(E1) values for the interbandE1 transitions from the positive-parity to negative-parity bands in153Eu forβ30=0.00 andβ30=0.05 , respectively. As shown in Fig. 4(a), theE1 transitions are completely due to the contribution of the intrinsic valence particle part in Eq. (3) since the collective dipole momentq(c)10=0 withβ30=0.00 . Forβ30=0.05 , as shown in Fig. 4(b), theE1 transitions are primarily due to the contribution of the collective part in Eq. (3). The contribution of the collective part increases and shows a staggering behavior at high spins as the spin increases, as does the calculatedB(E1) value. This staggering behavior may result from variations in the main components of the intrinsic wave functions.Figure 4. (color online) The calculated
B(E1) values for the interbandE1 transitions from the positive-parity to negative-parity bands in 153Eu for (a)β30=0.00 and (b)β30=0.05 . CalculatedB(E1) values using only the first term of Eq. (3) (labeled as coll.) and only the second term of Eq. (3) (labeled as s.p.) are denoted by the upside down and up triangles, respectively. The totalB(E1) values (labeled as Total) are denoted by open squares. Forβ30=0.00 , the collective part vanishes.The main components of the RAT-PRM wave functions in terms of the strong coupled basis
|IMK⟩χνp (denoted as|Kν⟩p for short) are shown in Fig. 5. Here,|IMK⟩ is the Wigner function, with I, M, and K denoting the quantum numbers of the total angular momentum and its projections along the third axis in the laboratory and intrinsic frames, andχνp representing the intrinsic wave function of the νth proton single-particle level|ν⟩p . The main component of the positive-parity band is|5/2,1⟩p , and that for negative-parity band is|5/2,4⟩p , i.e., the positive- and negative-parity bands are built on one-proton configurations|1⟩p and|4⟩p , respectively. As shown in Figs. 5(a) and (b), for the positive-parity band, the largest component|5/2,1⟩p plays an overwhelming role, with the amplitude always larger than 0.80, whereas the largest component of the negative-parity band|5/2,4⟩p decreases rapidly with spin, i.e., its amplitude decreases from 0.94 (5/2ℏ ) to 0.43 (35/2ℏ ) and 0.56 (37/2ℏ ). In comparison, the main components of the intrinsic wave functions for the positive- and negative-parity bands are nearly unchanged withβ30=0.05 , as shown in Figs. 5(c) and (d).Figure 5. (color online) The main components of the intrinsic wave functions expanded in the strong-coupled basis
|IMK⟩χνp (denoted as|K,ν⟩p for short) for (a) & (c) the positive-parity bandπ=+ and (b) & (d) negative-parity bandπ=− within spin region52ℏ≤I≤372ℏ withβ30=0.00 (left plane) andβ30=0.05 (right plane).|IMK⟩ is the Wigner function with I, M, and K being the quantum numbers of the total angular momentum and its projections along the third axis in the laboratory frame and intrinsic frame, respectively. Theχνp represent the intrinsic wave functions of the proton single-particle level|ν⟩p obtained by diagonalizing the reflection-asymmetric Nilsson Hamiltonian.To understand why the octupole deformation
β30 inherently results in a different performance ofB(E1) , as shown in Fig. 4, we further examine the main spherical harmonic oscillator components of the single-particle levels|1⟩p and|4⟩p of the proton forβ30=0.00 andβ30≠0.00 . Whenβ30=0.00 , the parity is a good quantum number and the spherical components with different parities cannot mix. The positive-parity level|1⟩p has the dominant componentg7/2 , mixed with components such asi11/2 andd5/2 . The negative-parity level|4⟩p has the dominate componenth11/2 , mixed with components such asj15/2 andf7/2 . Whenβ30≠0.00 , the parity is no longer a good quantum number and the spherical components with different parities can mix with each other. The level|1⟩p dominated byg7/2 is mixed with the opposite parity componenth11/2 , while the level|4⟩p dominated byh11/2 is mixed withg7/2 . Although only the intrinsic valence particle part in Eq. (3) contributes to theE1 transitions forβ30=0.00 , the intrinsic and collective parts contribute to the transitions forβ30≠0.00 . For the intrinsic part, the intrinsic components that can contribute nonzeroE1 single-particle matrix elements change a little asβ30 increases, resulting in the nearly unchanged single-particle contribution ofE1 transitions, as shown in Fig. 4. For the collective part, theE1 transitions can be enhanced by the matrix element between the largest componentg7/2 in|1⟩p , the same componentg7/2 mixed in|4⟩p , that between the largest componenth11/2 in|4⟩p , and the same componenth11/2 mixed in|1⟩p , etc. Since the probabilities of these matrix elements become significant asβ30 increases, the contributions of the collective part will dominate theB(E1) for a largeβ30 value, as shown in Fig. 4(b) forβ30=0.05 .In Ref. [31], the octupole correlations of the observed low-lying
K=5/2± positive- and negative-parity structures in nucleus151Pm , theN=90 isotonic neighbor of153Eu , have been investigated. The observed doublet bands in151Pm and153Eu show quite similar characteristics in both the energy spectra of positive- and negative-parity bands, the enhanced interbandE1 transitions, and in their different g factors of the opposite-parity states in the two bands. The similarity of the two nuclei suggests the same origin for these experimental characteristics. Based on the RAT-PRM calculations, the observed nearly-degenerate positive- and negative-parity bands in153Eu and151Pm can be interpreted as two separate bands based on a substantial reflection-asymmetric core and two individual proton configurations. This interpretation differs with the parity doublet bands built on a single parity-mixed configuration, e.g., that observed in223Th [32]. Further systematic studies of interest include the observed parity doublet bands and comparisons of the parity doublet bands built on a single parity-mixed configuration and those on two different configurations in theA∼150 mass region. Notably, although the deformation parameters from the microscopic CDFT are adopted, some free parameters still exist in the present RAT-PRM, such as the core parity splitting parameterE(0−) and the moment of inertia for reproducing the experimental data. Further combination with the microscopic theory, such as the cranking CDFT and the beyond mean-field approach in the future, is relevant to constrain these parameters. -
In summary, a reflection-asymmetric triaxial particle rotor model has been used to investigate the octupole correlations of the observed low-lying
K=5/2± positive- and negative-parity bands in the stable nucleus153Eu . The energy spectra, energy staggering parameters, and intrabandE2 andM1 transition probabilities are well reproduced by the RAT-PRM calculations. The calculated interbandB(E1) values withβ30=0.00 generally underestimate the experimental data. The calculatedB(E1) values are found to depend sensitively on the octupole deformation parameterβ30 , which can be enhanced up to approximately two orders of magnitude by changingβ30 from 0.00 to 0.07. With increasing octupole deformationβ30 , the calculatedB(E1) values increase, showing good agreement with the experimental data for theβ30=0.05 case, which indicates that only the single-particle contributions are unable to account for the largeE1 transitions observed in153Eu . By analyzing the main components of the RAT-PRM wave functions, the nearly degenerate positive- and negative-parity bands can be described by the bands based on two individual proton configurations, i.e.,πg7/2[Ω=5/2] andπh11/2[Ω=5/2] , respectively. This interpretation is the same as that for the isotonic neighbor151Pm but differs from the parity doublet bands built on a single parity-mixed configuration.
