-
The normally deformed (ND) and super deformed (SD) states of
36 Ar have been the research topics of several theoretical models [1, 2] since the discovery of the SD rotational band [3]. In recent years, properties of the hypernuclear system,37Λ Ar, have often been discussed [4–7], although this hypernuclear system has not been experimentally observed. The most popular topics associated with36 Ar and37Λ Ar are the SD states; the ground states are ND. For theΛ hyperon occupying the s-shell orbital, the relativistic mean field (RMF) model [4] and beyond-RMF calculations [7] indicate that the SD state of37Λ Ar gives a largerΛ separation energy,SΛ , than the ground state. However, antisymmetrized molecular dynamics (AMD) [5] and the Skyrme-Hartree-Fock (SHF) model [6] show that theSΛ of the SD state is smaller than that of the ground state. Compared to AMD [5] and the beyond-RMF calculations [7], the SHF model [6] merely provides mean-field results, which breaks rotational symmetry and thus hinders further uses of the rotational bands and the transition probabilities observed in laboratories. Since 2015, the beyond-mean-field SHF model forΛ hypernuclei has gradually improved [8–12]. The angular momentum projection (AMP) technique is used to restore rotational symmetry, and the generator coordinate method (GCM) is employed to deal with shape mixing. Therefore, the beyond-mean-field SHF model is used to research the rotational bands corresponding to the configuration36 Ar⊗ s Λ , which is the first of the three motivations of this paper.Beside the s-shell orbital, the
Λ hyperon may be excited onto the p-shell orbitals, and such aΛ hyperon coupled with an ND or SD nuclear core can provide various configurations. These configurations correspond to low-lying energy levels with negative parity, and the beyond-mean-field SHF model can be used to predict them. For example, in a previous study [10], the rotational energy levels corresponding to the configurations8 Be⊗ p∥Λ and8 Be⊗ p⊥Λ [13, 14] were precisely reproduced. Therefore, the second motivation of this paper is to obtain predictions for configurations such as36 Ar⊗ pΛ , which will help identify the negative-parity energy levels of37Λ Ar in future experiments.For several typical hypernuclear systems such as
9Λ Be,13 Λ C,21 Λ Ne, and25Λ Mg,Λ hyperons on the s-shell orbital or excited onto the p-shell orbitals have been studied using various theoretical models [10–12,15–18] in recent years. However, the sd-shellΛ hyperon in these hypernuclei does not provide a bound state because they are too light to bind such aΛ hyperon. Some (K− ,π− ) and (π+ , K+ ) reactions [19–21] reveal that the configuration39 Ca⊗sdΛ gives bound energy levels (that is, energy levels below the separation threshold of39 Ca+Λ ), whereas31 S⊗sdΛ gives unbound energy levels.37Λ Ar is between32Λ S and40Λ Ca in the hypernuclear chart, which is probably to provide bound states for configurations such as36 Ar⊗ sdΛ . Hence, the third motivation of this paper is to theoretically investigate the existence of bound states for36 Ar⊗ sdΛ and obtain the corresponding configurations with Nilsson quantum numbers.In this paper, a beyond-mean-field SHF model that considers AMP techniques and the GCM is employed to study the low-lying energy levels corresponding to various configurations of
37Λ Ar. These configurations include an ND or SD nuclear core coupled to aΛ hyperon onto the orbitals of s, p, or sd shells. This paper is organized as follows. In Sec. II, the formalism of the beyond-mean-field calculation is introduced. Sec. III presents the results and discussions. In Sec. IV, we draw conclusions about the study. -
In a hypernuclear system, the physical quantities we are most concerned with are its eigenstates, eigenenergies, and transition rates. To obtain these physical quantities, the beyond-mean-field SHF model begins with a mean-field state of the
Λ hypernucleus in the intrinsic frame of reference [22]|Φ(NΛ)(β)⟩=|ΦN(β)⟩⊗|ΦΛ⟩,
(1) where
|ΦN(β)⟩ represents a nuclear core with the deformation parameter β, and|ΦΛ⟩ denotes the single-particle wave function of aΛ hyperon.A general eigenstate of the
Λ hypernuclear system is given by linear superposition:|ΨJMα⟩=∑βFJα(β)ˆPJMK|Φ(NΛ)(β)⟩,
(2) where β plays the role of the generator coordinate, and the AMP operator,
ˆPJMK , restores the rotational symmetry of the system. The mean-field state in the intrinsic frame of reference,|Φ(NΛ)(β)⟩ , is derived by a constrained SHF model [6], and a density-dependent delta interaction (DDDI) [23] is considered for the nucleons asG(r)=V0[1−ρ(r)ρ0].
(3) The weight function
FJα(β) in Eq. (2) and the eigenenergies are determined by the Hill-Wheeler-Griffin (HWG) equation [24],∑β′[H′JKK(β,β′)−EJαNJKK(β,β′)]FJα(β′)=0,
(4) in which the Hamiltonian and norm elements are given by
H′JKK(β,β′)=⟨Φ(NΛ)(β′)|^H′ˆPJKK|Φ(NΛ)(β)⟩,
(5) NJKK(β,β′)=⟨Φ(NΛ)(β′)|ˆPJKK|Φ(NΛ)(β)⟩.
(6) The corrected Hamiltonians
^H′ is^H′=ˆH−λp(ˆNp−Z)−λn(ˆNn−N),
(7) where the final two terms on the right-hand side account for the fact that the projected wave function does not provide the correct number of particles on average [25–27], and
ˆH is given by the energy density functional (EDF) in Ref. [28].For a certain configuration determined by the intrinsic state
|Φ(NΛ)(β)⟩ , the projected energy is derived asE(β,J,K)=H′JKK(β,β)NJKK(β,β).
(8) Given two general eigenstates,
|α;J⟩ and|α′;J′⟩ , the reducedE2 transition rate between them is derived asB(E2,Jα→J′α′)=12J+1|⟨α′;J′||ˆQ2||α;J⟩|2,
(9) where
⟨α′;J′||ˆQ2||α;J⟩=√2J′+1∑Mμββ′FJ′α′∗(β′)FJα(β)×CJ′K′JM2μ⟨Φ(NΛ)(β′)|ˆQ2μˆPJMK|Φ(NΛ)(β)⟩.
(10) In the above equation,
ˆQ2μ=r2Y2μ(φ,θ) is the electric quadrupole operator [29], andCJ′K′JM2μ denotes the Clebsh-Gordon coefficients.Parameters: In this paper, the SLy4 force is used for the
NN interaction, and the strength of the pairing force isV0=−410 MeV fm3 for both protons and neutrons [30], with a smooth pairing energy cutoff of 5 MeV around the Fermi level [31, 32]. For theNΛ interaction, we use the SLL4 force [28], which offers the best fit for theΛ separation energy in the spherical SHF model [28].Notations: Single-particle orbitals of the
Λ hyperon are denoted by Nilsson quantum numbers[Nn3ml]Ωπ ; the same notation is also adopted in Ref. [4]. Then, a certain configuration of theΛ hypernuclear system is given asA−1 Z⊗Λ[Nn3ml]Ωπ , in whichA−1 Z is the nuclear core, andΛ[Nn3ml]Ωπ represents aΛ hyperon occupying the orbital[Nn3ml]Ωπ . To emphasize the deformation,36 ArN and36 ArS are used to represent the ND and SD nuclear cores of37Λ Ar, respectively. For example, the configuration36 ArS⊗Λ [110]12− denotes an SD nuclear state of36 Ar coupled with aΛ hyperon occupying the orbitalΛ [110]12− .Model space: Because intrinsic wave functions are kept axially symmetric, the quadrupole deformation parameter β of the nuclear core is the exclusive generator coordinate. If the
Λ orbital is specified, the range of β and the number of basis functions determine the model space jointly. In this paper, β is between−1.4 and2.8 , and 90 basis functions are evenly spaced over this range for both36 Ar and37Λ Ar. -
The potential energy surfaces (PESs) of
36 Ar and37Λ Ar derived using the mean-field calculation are shown in Fig. 1, and we can see that there are three energy minima for36 Ar. Two of these minima, located atβ=−0.16 andβ=0.12 , both give ND configurations, whereas the third, located atβ=0.63 , exhibits an SD prolate shape. For37Λ Ar, the addition of oneΛ hyperon onto the orbitals of thes, p, or sd shells causes various impurity effects. In our current mean-field calculation, the hyperon Λ [000]12+ makes|β| of the three energy minima slightly smaller, which is already known from early research [33]. Such a reduction in deformation is not clear in Fig. 1 because the nuclear core36 Ar is sufficiently heavy against theNΛ attraction. For the p-shellΛ hyperon, Fig. 1 clearly shows thatΛ [110]12− in37Λ Ar makes the prolate minima significantly deeper and nearly eliminates the oblate minimum, whereasΛ [101]32− andΛ [101]12− make the prolate energy minima shallower and enhance the oblate minimum. These opposite phenomena are due to the density distributions of the p-shellΛ orbitals. Fig. 2 shows the density distributions of several p-shell and sd-shellΛ s. It must be emphasized that the density distributions ofΛ [101]32− andΛ [101]12− are nearly identical and share the same contour map in Fig. 2, andΛ [202]52+ andΛ [202]32+ do the same. In the top-right panel of Fig. 2, we can see thatΛ [110]12− is prolately distributed; thus, its coupling to a prolate nuclear core leads to a lower binding energy. However, the top-left panel of Fig. 2 shows thatΛ [101]32− andΛ [101]12− are oblately distributed; hence, they prefer oblate nuclear cores. For the sd-shellΛ hyperon, Fig. 1 shows the PESs of configurations36 Ar⊗ Λ [220]12+ ,36 Ar⊗ Λ [202]52+ , and36 Ar⊗ Λ [202]32+ because only these three configurations give bound states, the properties of which will be discussed in detail in the final two paragraphs of this section.Figure 1. (color online) PESs obtained from the mean-field calculations for certain configurations indicated in the legend. The horizontal axis represents the quadrupole deformation parameter β.
Figure 2. (color online) Density distributions of the
Λ hyperon on the orbitalsΛ [101]32− orΛ [101]12− (top-left pannel),Λ [110]12− (top-right panel),Λ [202]52+ orΛ [202]32+ (bottom-left panel), andΛ [220]12+ (bottom-right pannel). The deformation parameter β of each corresponding nuclear core is constrained to the average,ˉβ , of the bandhead given by the GCM calculation.Figure 3 gives the projected PESs,
E(β,J,K) , on each angular momentum J for36 Ar. In this figure, all the ND and SD energy minima of theJπ=0+ PES are more obvious than those of the mean-field PES, which is due to the energy gained from the restoration of rotational symmetry. The energy levels derived by the GCM are also given in the same figure, and we can clearly see that there are two rotational bands for the ND states and the SD states of36 Ar. In this current calculation, the bandhead of the ND band corresponds to the ground state (g.s.), whereas that of the SD band is the 4th0+ state. The bandheads of the ND and SD bands are located at –308.02 MeV and –300.96 MeV, respectively, which indicate the binding energy of these two states. Figure 3 also shows that the g.s. and SD bands cross atJπ=6+ , and these calculated results agree with the observed data in Ref. [3]. Besides the ND and SD rotational bands, there are other observed low-lying energy levels such as the 4.4 MeV 2+ level [3]; however, these levels may involve two-quasi-particle (2-qp) excitations [2]. 2-qp excitations are beyond the basis space of this current model; therefore, only the properties of the g.s. and SD bands are discussed hereafter.Figure 3. (color online) Projected PESs,
E(β,J) , and the GCM energy levels of36 Ar. The angular momentum and parity for each projected PES are given in the legend, and the mean-field PES labeled by 'MF' is also shown for comparison. The solid bullets and horizontal bars indicate the GCM energy levels, which are plotted at their average deformationˉβ .The projected PESs of
36 Ar⊗ Λ [000]12+ are given in Fig. 4, and the GCM energy levels of the g.s. and SD bands are also shown in the same figure. The addition ofΛ [000]12+ turns theJπ=0+ states of36 Ar into those withJπ=12+ and turns the other states into spin doublets. The band-head energies of the g.s and SD bands are –327.22 MeV and –319.79 MeV, respectively. Therefore, we can deduce that theΛ separation energy,SΛ , for the ground state is 19.20 MeV, whereasSΛ for theJπ=12+ SD state is 18.83 MeV, which indicates thatSΛ of the SD state is smaller than that of the ground state. This is in agreement with the mean-field SHF calculation using the SKI4 parameters [6] but conflicts with the RMF calculations [4, 7]. This conflict stems from the fact that the RMF calculations give an SD state with a localized density, which leads to a larger overlap betweenΛ and the SD nuclear core. This does not occur in the SHF calculation.Figure 4. (color online) Same as Fig. 3, but for
37Λs Ar. The notationΛs denotes aΛ hyperon on the s-shell orbital, that is,Λ [000]12+ , and only the GCM energy levels of the g.s. and SD bands are preserved for simplification.Figures 5 and 6 give the collective wave functions
gJα(β) of the ND-band and SD-band states for36 Ar and37Λ Ar, respectively. These two figures reproduce the observed features of the ND and SD bands qualitatively. For each state of the ND band, it is shown that the maximum ofgJα(β) corresponds to β within the ND region; for the SD-band states, all the collective wave functions reach their maxima atβ≈0.65 , which indicates an SD shape. It is also shown that, compared to those of36 Ar, the s-shellΛ slightly pulls thegJα(β) of the hypernuclear states in37Λ Ar toward the spherical shape, which leads to a reduction in deformation.Figure 5. (color online) Collective wave functions
gJα(β) for ND-band states in36 Ar and37Λ Ar. The angular momenta, J, are shown in the legend of each panel.Figure 6. (color online) Same as Fig. 5, but for the SD-band states in
36 Ar and37Λ Ar.To discuss the reduction in deformation more comprehensively, the g.s. and SD bands of
36 Ar and37Λ Ar in Figs. 3 and 4 are extracted and shown in Fig. 7, where the observed data [34, 35] on36 Ar are also given for comparison. UntilJπ=4+ , the current model reproduces the observed g.s. band of36 Ar very well. Moreover, the intervals between the energy levels of the calculated SD band are in good agreement with the observed data. However, the current model gives a significantly higherJπ=6+ energy level for the g.s. band than the observed value, and the predicted SD band is nearly 3 MeV greater than the experimental data. This disagreement with the observed data likely stems from the fact that quasi-particle excitations are absent in the GCM space of our model. In Fig. 7, it is also shown that the addition ofΛ [000]12+ shifts the excited energy levels of37Λ Ar slightly upward compared to those of36 Ar because the moment of inertia of the nuclear core is reduced owing to the reduction in defromation. The reduction in deformation is also indicated by the intra-bandB(E2) values of37Λ Ar in comparison with those of36 Ar. As shown in Fig. 7,B(E2,52+→12+) in the g.s. band of37Λ Ar is reduced by9% compared toB(E2,2+→0+) in the g.s. band of36 Ar, whereas for the SD band of37Λ Ar,B(E2,52+→12+) is reduced by approximately5% . The contrast between the reductions in theB(E2) values of the g.s. band and those of the SD band indicates that deformation of the SD band is more stable than that of the g.s. band. For the g.s. band, the stability of deformation is also indicated by the PESs given in Fig. 1, which show that there are two ND energy minima and the barrier between them is less than 1.5 MeV. Furthermore, the collective wave functions in Fig. 5 show that the two ND energy minima compete with each other. Finally, it is deduced that shape deformation in the g.s. band is relatively soft and unstable. Note that the beyond-mean-field RMF calculation provides a clearer reduction in theB(E2) values in37Λ Ar [7]. This is because the relativisticNΛ interaction is stronger than the nonrelativistic interaction, as verified in Refs. [6, 9].Besides
Λ [000]12+ , the ND and SD nuclear cores, denoted as36 ArN and36 ArS in this paper, may couple with aΛ hyperon on the p- or sd-shell orbitals, and each specific combination corresponds to a certain configuration of37Λ Ar. Energy levels corresponding to these configurations, which are lower than the separation threshold of36 Ar+Λ , are given in Fig. 8. It is shown that the configurations36 ArN⊗Λ [110]12− ,36 ArN⊗Λ [101]32− , and36 ArN⊗Λ [101]12− give three rotational bands with bandheads at 11.13 MeV, 9.02 MeV, and 8.80 MeV, respectively. Because the ground state of36 Ar is oblate, the configurations36 ArN⊗Λ [101]32− and36 ArN⊗Λ [101]12− give lower energies than the configuration36 ArN⊗Λ [110]12− . However, for the SD state, the configuration36 ArS⊗Λ [110]12− gives a rotational band with the bandhead at 13.70 MeV, which is approximately 5 MeV lower than the energy levels given by the configurations36 ArS⊗Λ [101]32− and36 ArS⊗Λ [101]12− . This is because the SD state of36 Ar is prolately deformed, and its coupling withΛ [110]12− produces lower energies. Figure 8 also shows that, although the configurations36 ArS⊗Λ [101]32− and36 ArS⊗Λ [101]12− both give bound states, their energy levels are near the separation threshold of36 Ar+Λ . In our earlier research [10], the beyond-mean-field SHF model successfully reproduced the negative-parity energy levels of9Λ Be, which were denoted as a genuine hypernuclear state and9 Be-analogue. Furthermore, the same model predicted the negative-parity levels of13 Λ C, which was in good agreement with the observed values [11]. Therefore, the beyond-mean-field SHF model is powerful, and the predicted energy levels in Fig. 8 for the p-shellΛ are reliable and can help identify the negative-parity levels of37Λ Ar in future experiments.Figure 8. (color online) Comprehensive energy spectra for
37Λ Ar. Configurations corresponding to each rotational band are indicated by the solid box and dashed arrow. Energy levels higher than the separation threshold of36 Ar +Λ are neglected for simplification.As deduced in Sec. I,
37Λ Ar may provide bound states for configurations such as36 Ar⊗ sdΛ . Thus, in future experiments, theoretical predictions for the energy levels of36 Ar⊗ sdΛ are crucial for identifying the bound state near the separation threshold of36 Ar+Λ . The current calculation finds bound states for configurations36 Ar⊗ Λ [220]12+ ,36 Ar⊗ Λ [202]52+ , and36 Ar⊗ Λ [202]32+ , the energy levels of which are shown in the top right corner of Fig. 8. For36 ArN⊗Λ [220]12+ and36 ArN⊗Λ [202]52+ , only the bandheads withJπ=12+ and52+ are lower than the separation threshold, and their energies are 18.45 MeV and 17.91 MeV, respectively. For36 ArN⊗Λ [202]32+ , the bandhead is withJπ=32+ , and its energy is 17.94 MeV. For the same configuration, there is another bound state withJπ=52+ ; however, the energy level is only 0.24 MeV lower than the separation threshold.The mean-field PESs of the three configurations discussed above are shown in Fig. 1. It is found that
Λ [220]12+ enhances the energy minima on the prolate side of the PES and weakens the oblate one, whereasΛ [202]52+ andΛ [202]32+ exhibit opposite effects. Similar to cases of the p-shellΛ , impurity effects of those in the sd-shell mainly stem from their density distributions, which are shown in Fig. 2. ForΛ [220]12+ , the main component of its wave function isY2,0(θ,φ) and it is elongated along the symmetry axis (shown in the bottom-right panel of Fig. 2); hence, it prefers a prolate nuclear core. Conversely, forΛ [202]52+ andΛ [202]32+ , the main components of their wave functions areY2,±2(θ,φ) ; therefore, they are both compressed in the direction of the symmetry axis (shown in the bottom-left panel of Fig. 2) and prefer oblate nuclear cores. As shown in Fig. 8, this type of effect also causes the GCM energy levels of the configurations36 ArN⊗Λ [202]52+ and36 ArN⊗Λ [202]32+ to be approximately 0.5 MeV lower than that of36 ArN⊗Λ [220]12+ . This is because, for all three configurations, the nuclear core36 ArN is oblately distributed, and its coupling withΛ [202]52+ orΛ [202]32+ leads to relatively lower energies. -
In summary, a beyond-mean-field SHF model that includes AMP techniques and the GCM is introduced in this paper. Based on this model, a comprehensive energy spectrum of
37Λ Ar is given and researched in detail. Attention is mainly paid to the bound states that are formed by an ND or SD nuclear core coupled with aΛ hyperon occupying one of the s-, p-, or sd-shell orbitals. For the s-shellΛ hyperon, ND and SD bands corresponding to the configurations36 ArN⊗Λ [000]12+ and36 ArS⊗Λ [000]12+ are given. The intra-bandB(E2) values of these two bands decreased compared to those of36 Ar, which is due to the reduction in deformation. For the p-shellΛ hyperon, three ND bands and one SD band are shown, which correspond to the configurations36 ArN⊗Λ [110]12− ,36 ArN⊗Λ [101]32− ,36 ArN⊗Λ [101]12− , and36 ArS⊗Λ [110]12− , respectively. Moreover, the configurations36 ArS⊗Λ [101]32− and36 ArS⊗Λ [101]12− give two bound states, the energy levels of which are below the separation threshold of36 Ar+Λ within 0.5 MeV. For the sd-shellΛ hyperon, three configurations36 ArN⊗Λ [220]12+ ,36 ArN⊗Λ [202]52+ , and36 ArN⊗Λ [202]32+ give bound states with energy levels lying at 18.45 MeV, 17.91 MeV, and 17.94 MeV above the ground state, respectively. Finally, it is expected that the predictions in this paper will aid future experiments in identifying the observed states of37Λ Ar. -
One of the authors thanks Xian-Ye Wu for informative discussions.
Beyond-mean-field study of 37Λ Ar based on the Skyrme-Hartree-Fock model
- Received Date: 2022-01-10
- Available Online: 2022-07-15
Abstract: We present the hypernuclear states of