-
Unstable nuclei far from the stability line in nuclear charts exhibit numerous novel phenomena, for instance, the emergence of new magic shells and the disappearance of conventional ones [1–7], the island of inversion [8–10], and dilute matter distributions – halo structures [11–13], which are attracting the interest of many researchers in the fields of nuclear physics and related disciplines. While enriching the knowledge of nuclear physics, the emergence of novelties in unstable nuclei has also subverted our traditional understanding of atomic nuclei.
Beryllium (Be) isotopes are perfect candidates for exploring the novelty appearing in unstable nuclei. The novel parity inversion of 11Be is one of the significant issues. As deduced from the conventional shell model, the ground state (GS) of 11Be must have odd-parity because the odd neutron should populate the lower
p1/2 orbit. However, an early experiment with the reaction 11B(n,p ) 11Be indicated thatJπ=1/2− is not the GS [14], which was further proved to be even-parity1/2+ by later experiments [15–17]. At a rather early stage, it was suggested that the even- and odd-parity states1/2± of 11Be should be very close [18]. After a long time of silence, there were a number of attempts to interpret the GS parity of 11Be under the particle-core coupling picture [19–24], indicating that quadruple core excitation and pairing blocking can be essential [20, 21]. Besides, the deformation [25–27], the continuum [28] and effects beyond the mean field [22, 29] might also be significant for correctly giving the even-parity GS of 11Be.Despite the ingenious GS, special interests were also devoted to the halo occurrence in 11Be. Experimentally, 11Be has been confirmed to be of the one-neutron halo structure [30–32]. Theoretically, the energy density functional (EDF) theory presents dilute matter distribution for 11Be, while it fails to reproduce the even-parity GS [26, 27, 33]. Based on the variational shell model, it was indicated that the weak couplings between the core and odd neutron can be responsible for the one-neutron halo in 11Be [19]. Recently, the halo effects of 11Be on the one-neutron knockout reactions were studied using the halo effective field theory [34, 35]. As a consensus, nuclear halo structures are contributed to by neutrons or protons populating the weakly bound low-l states, and the neighboring high-l orbits do not contribute to the halo but stabilize the dilute nucleons via pairing correlations [36].
However, up to now both novelties in 11Be, i.e., the parity inversion and neutron halo, still cannot be interpreted in a unified framework, and the halo occurrence may essentially rely on the even-parity GS. The EDF theory, such as that adopted by the widely used Skyrme Hartree-Fock models, cannot uniformly describe the even-parity GS and halo occurrence in 11Be, as aforementioned. The relativistic mean field (RMF) theory [36–41], that stands on the microscopic meson exchange picture of the nuclear force, succeeds in describing the halo structures of spherical and deformed nuclei [36, 42], but fails to reproduce the parity inversion of 11Be. Implemented with Fock diagrams, the relativistic Hartree-Fock (RHF) theory [43–47] and an extension – the relativistic Hartree-Fock-Bogoliubov (RHFB) theory [48, 49] have achieved many successes in describing various nuclear phenomena [50], including those of unstable nuclei [51–54] and superheavy ones [55].
In particular, the degrees of freedom associated with the π-pseudo-vector (π-PV) and ρ-tensor (ρ-T) couplings, important ingredients of the effective nuclear force, play significant roles in providing self-consistent descriptions of nuclear shell evolution [56–58], the pseudo-spin symmetry restoration [47, 59] and the emergences of new magicity [60, 61]. Recently, considering nuclear deformation that may be significant for 11Be, the axially deformed RHF [62] and RHFB [63] models, respectively called the D-RHF and D-RHFB models, have been established by utilizing an expansion of the spherical Dirac Woods-Saxon (DWS) base [64]. It was revealed that the π-PV and ρ-T couplings, which contribute mainly via the Fock diagrams, evidently enhance the deformation effects of nuclei [62, 63].
In particular, the D-RHFB model has the advantage of a unified treatment of spin-orbit coupling, tensor force, deformation, pairing correlations and the continuum [63], the significant mechanisms which are essential for unstable nuclei. This encourages us to explore the underlying mechanism related to the even-parity GS and halo occurrence in 11Be. The paper is organized as follows. In Sec. II, the decomposition of the canonical single-particle energy and its physical significance are introduced. Afterwards, the parity inversion and the halo phenomenon are discussed in Sec. III. Finally, we give a summary in Sec. IV.
-
To understand the microscopic picture of 11Be, the language of the basis expansion is rather helpful. In the D-RHFB model, the spherical DWS base is used to expand the quasi-particle wave functions
ψV andψU [63],ψVνπm=∑aCVνπm,aψam,ψUνπm=∑aCUνπm,aψam,
(1) where m and π denote the angular momentum projection and parity, ν stands for the index of the orbits in the
πm block. Imposing axial symmetry and reflection symmetry, m and π remain as good quantum numbers. For a spherical DWS base,ψam represents the wave function and the index a =(nκ) carries the principle number n and κ-quantity, namelyκ=±(j+1/2) withj=l∓1/2 , where l is the orbital angular momentum. In the following, the indexi=(νπm) is used to denote deformed orbits, while the index a specifies the states in the spherical DWS base. It is worth noting that the DWS base has the advantage of providing appropriate asymptotic behaviors for the wave functions due to the nature of the Woods-Saxon type potential [65]. For unstable nuclei, e.g., 11Be with a halo structure, this advantage is essential to provide a reliable description.For the convenience of analysis, the canonical transformation from the quasi-particle space is performed in general to get canonical s.p. states by diagonalizing the density matrix. In practice, the canonical s.p. space is taken to be similar to the Hartree-Fock space, for both of which the density matrices are diagonal. Utilizing the orthogonality of the DWS base, there is a straightforward way to build the density matrix, which reads as
ρπmaa′=∑νCVνπm,aCVνπm,a′.
(2) Notice that for the contributions from the odd-nucleon in the blocked orbit, the expansion coefficient
CV shall be replaced by theCU one in Eq. (1). The diagonalization of the above density matrixρπmaa′ gives the eigenvaluesv2i , namely the occupation probability of orbiti=mπν , and the eigenvectorsˆDπmν , the set of the expansion coefficientsDi,nκ of the canonical wave functions in the spherical DWS base.Thus, in terms of the spherical DWS base, the canonical s.p. energy
εi and wave functionψi for orbit i can be expressed as:εi=∑aa′Di,ahπmaa′Di,a′,ψνπm=∑nκDi,nκψnκm,
(3) where
haa′ is the s.p. Hamiltonian in the spherical DWS base, the index a is for spherical DWS basis states. Following the deformation of a nucleus, orbits which originally had the same sphericalnκ quantities cannot degenerate any more, leading to a mixture of sphericalnκ -components in deformed s.p. orbits. As seen from Eq. (3), the couplings between various spherical DWS components a and a' give the canonical s.p. energy. Thus, after summing the principle numbers n andn′ , they can be decomposed into the diagonal termsEκd,i and off-diagonal onesEκκ′c,i as:εi=∑κκ′Eκκ′i=∑κp2iκEκd,i+∑κ≠κ′piκpiκ′Eκκ′c,i,
(4) where
p2iκ describes the sum proportion of the κ-component of the orbitmπν , andEκκ′i corresponds to the coupling effects between the spherical κ andκ′ components,pik=∑nD2i,nκ,Eκκ′i=∑nn′Di,nκhπmaa′Di,n′κ′.
(5) Using the above conventions, the diagonal term
Eκd,i and the off-diagonal oneEκκ′c,i can be defined asEκd,i≡Eκκi/p2iκ,Eκκ′c,i≡Eκκ′i/(piκpiκ′).
(6) Further, considering the sum over κ components, the diagonal contribution
Ed,i and the off-diagonal oneEc,i to the canonical s.p. energyεi read asEd,i=∑κEκκi,Ec,i=∑κ≠κ′Eκκ′i.
(7) As aforementioned, more spreading of the proportion
p2iκ over different κ-components is in general due to enhanced deformation effects. Thus, using the language of the spherical DWS basis expansion, theEκd,i -term can be taken as the centroid energy of the κ-component of the deformed orbit i, or taken as the single-particle-like energy of the spherical κ-fragment of orbit i. In contrast, theEκκ′c,i -term describes the inherent correlations between the κ- andκ′ -components, which in fact reflects the shape evolution of orbit i. Following deformation, it is not hard to deduce that orbit i becomes more bound with the negativeEc,i value, and vice versa. -
Aiming at a unified mechanism behind the even-parity GS and neutron halo of 11Be, we focus on the D-RHFB calculations with PKA1 [47], one of the most complete RHF Lagrangians which contains the π-PV and ρ-T couplings. For comparison, the RHF Lagrangians PKO3 [56] containing π-PV coupling and PKO2 [56], and the RMF one DD-LZ1 [66] are considered as well. In the pairing channel, the finite-range Gogny force D1S [67] is utilized as the pairing force, due to the advantage of a finite range, i.e., natural convergence with the configuration space. The details of the calculations can be found in Ref. [63].
To find the global minimum for an odd-A nucleus using the D-RHFB model, the blocking effects shall be treated carefully [68]. Similar to a popular treatment in the RMF calculations, the equal filling approximation is adopted [69, 70]. Within the Bogoliubov scheme, we attempted to block the first quasi-particle level for various
πm -blocks to determine the ground state. According to the mapping relation between the canonical single-particle and the Bogoliubov quasi-particle states, the blocked quasi-particle orbits are mapped in general with the s.p. ones closest to the Fermi level. For 11Be, we tried to block the statesmπ= 1/2− ,1/2+ ,3/2+ and5/2+ , and only show the results ofmπ=1/2− and1/2+ , from which the ground state can be determined by selected models. -
We first performed the shape constrained calculations for 11Be using PKA1, PKO2, PKO3 and DD-LZ1. Figures 1(a) and 1(b) show the binding energy
EB (MeV) of 11Be as a function of the quadruple deformation β, which corresponds to the blocking of even- and odd-parity statesmπ=1/2± , respectively. It is shown that only PKA1 results in the same even-parity GS as the experiments [14–16]. Moreover, a rather large prolate deformation (β≃1.2 ) is obtained by PKA1 for the GS, while the others give an odd-parity GS with spherical or near spherical shapes. Compared with the experimental value, which reads as –65.478 MeV [71], the binding energy given by PKA1 is overestimated by roughly 3 MeV. In fact, the calculations of nearby nuclei show that the binding energies for light nuclei are overestimated by PKA1in general.Figure 1. (color online) Binding energies
EB (MeV) as functions of the quadruple deformation β for 11Be with even [plot (a)] and odd [plot (b)] parity. The results are calculated by PKA1, PKO2, PKO3 and DD-LZ1. Referring toβ=0 , plots (c) and (d) show the PKA1 contributions ofEB respectively with even and odd parity, including the total, the π-PV and ρ-T ones, and the sum of the other channels.As an example, Table 1 shows the calculated binding energies
EB of 10,11,12Be, as well as the one-neutron separation energySn for 11Be and the two-neutron oneS2n for 12Be, in comparison with the experimental data [72]. It can be found that the systematics of the binding energies, here from 10Be to 12Be, are different between PKA1 and the other selected models. Specifically, allEB values given by PKA1 are larger than the experimental ones by1.4∼3.5 MeV. While the binding energy of 10Be is underestimated by all other selected Lagrangians, the one of 12Be is overestimated by different extents. However, for theSn andS2n values, all selected models provide a similarly accurate result as the experimental values. Thus, deviations in the binding energy of a few MeV do not influence the conclusions made in this work substantially.EB Sn S2n 10Be 11Be 12Be 11Be 12Be Exp. −64.98 −65.48 −68.65 0.50 3.67 PKA1 −66.40 −68.50 −72.21 2.10 5.81 PKO3 −63.63 −65.82 −70.13 2.19 6.50 PKO2 −63.50 −65.73 −69.00 2.23 5.50 DD-LZ1 −63.61 −66.09 −69.67 2.47 6.06 Table 1. Binding energies
EB (MeV) of 10,11,12Be, along with the one-neutron separation energy (MeV)Sn for 11Be and two-neutron oneS2n for 12Be calculated by PKA1, PKO2, PKO3, and DD-ME2, as compared with the experimental (Exp.) data [72].In order to clarify the mechanism behind the GS of 11Be, Figs. 1 (c) and 1 (d) show the contributions of
EB referring toβ=0 given by PKA1 for even- and odd-parity cases, respectively, namely the total, the π-PV one that is enclosed in both PKA1 and PKO3, the ρ-T one which is uniquely considered by PKA1, and the sum from other channels, which corresponds to all degrees of freedom adopted in PKO2 and DD-LZ1. It is interesting to see that both the π-PV and ρ-T contributions are evidently enhanced following the enlarged deformation. Compared with the odd-parity one [plot (d)], more remarkable enhancements are presented by the π-PV and ρ-T couplings for the even-parity case [plot (c)], which are seen to be crucial for reproducing the even-parity GS of 11Be.In fact, referring to other selected models that fail to reproduce the GS of 11Be, it is also helpful to understand the role played by the π-PV and ρ-T couplings. In Figs. 1 (a) and 1(b), both PKO2 and DD-LZ1, which do not contain either π-PV or ρ-T couplings, present the odd-parity GS with a fairly stable spherical shape and even-parity local minima with a rather soft prolate shape. For PKO3, which contains the π-PV coupling, an odd-parity GS with a rather soft near-spherical shape and an even-parity local minimum with a solid prolate shape are obtained, and the energy difference between the even- and odd-parity minima are significantly reduced, as compared with PKO2 and DD-LZ1. Thus, combined with the PKA1 results, it can be concluded that both π-PV and ρ-T couplings, which evidently enhance the deformation effects, play a key role in determining the even-parity GS of 11Be.
To understand why the even-parity GS of 11Be can be correctly reproduced at large deformation in terms of the spherical DWS base, the expansion proportions of the even-parity orbit
1/2+2 occupied by the odd neutron are shown as functions of the deformation β in Fig. 2 (a). Figure 2 (b) shows the canonical s.p. energiesEi (MeV) with occupation probabilities (in thick bars) for spherical 11Be and the GS, and theEκd,i values deduced from Eq. (6) with the expansion proportions (in thick bars) of the dominant κ fragments of the orbits1/2±2 for the GS. The results are extracted from the PKA1 calculations by blocking the statemπ=1/2+ . As shown in Fig. 2 (a), the1d5/2 component in the orbit1/2+2 increases gradually with a simultaneously reduced2s1/2 component when 11Be becomes prolately deformed. Meanwhile, the inversion of the even- and odd-parity orbits1/2±2 from the spherical 11Be to the GS is found in Fig. 2(b) and is consistent with an even-parity GS.Figure 2. (color online) Expansion proportions of neutron orbit
1/2+2 with respect to the deformation β [plot(a)], and the canonical s.p. energiesEi (MeV) and occupation probabilities (in thick bars) of spherical (Sph.) 11Be and the GS [plot (b)]. The results are extracted from the PKA1 calculations by blocking the statemπ=1/2+ , and theEκd,i values (MeV) and expansion proportions (in thick bars) of the neutron orbits1/2±2 of the GS are also shown in plot (b).Moreover, as deduced from the high-lying κ fragments of the
1/2+2 orbit in Fig. 2 (b), theEc,i value given by Eq. (7) shall be negative and large, indicating strong attractive inherent correlations between the κ fragments of the1/2+2 orbit. In contrast, the κ fragments of the odd-parity orbit1/2−2 are deeply bound, and the deduced large, positiveEc,i value reveals strong repulsive inherent correlations between the κ fragments of the1/2−2 orbit. Thus, from the opposite effects carried by the correlation termsEc,i , it is not hard to understand the inversion of the even- and odd-parity orbits1/2±2 from spherical 11Be to the deformed GS, in which the enhancements due to the ρ-T and π-PV couplings are very important as deduced from the nature of the channels [63].In fact, the PKA1 results are consistent with the indication by Bohr and Mottelson that for a p-shell nucleus 11Be the even-parity GS can be a consequence of a strong preference for the orbit
[220,1/2] with prolate deformation [73]. Supplementary, it should be emphasized that the expansion proportions in the theoretical calculation are different from the spectroscopic factor, which is a basic quantity characterizing the single-particle nature of nuclear excitations in both experiments and theories. In this work, the calculations of 11Be are performed by restrictions on the mean-field level. Thus, one cannot take the expansion proportions in Fig. 2 (a) as some equivalent concept of the spectroscopic factor, or even compare it directly with the measured data of the transfer reaction (see Ref. [32] and references therein).To enhance the understanding of the largely deformed GS and the role played by the π-PV and ρ-T couplings, Fig. 3 shows the interaction matrix elements (IMEs) between selected neutron states and core orbits of 11Be, including the total (upper plots) and the sum of the ρ-T and π-PV couplings (lower plots). The results are given by the spherical RHFB calculations [48] with PKA1 by blocking the even-parity state
s1/2 . In fact, the selected even-parity statesν1s1/2 ,ν2s1/2 , andν1d5/2 are the dominant orbital components of the deformed orbit1/2+2 , as seen from Fig. 2 (a).Figure 3. (color online) Interaction matrix elements (MeV) between a selected neutron (ν) and the core orbits of 11Be. The left and right plots correspond to the couplings with proton (π) and neutron (ν) cores, respectively, including the total (upper plots) and the sum of the ρ-T and π-PV couplings (lower plots). The results are obtained from spherical RHFB calculations with PKA1 by blocking the neutron s state.
As the general feature of the IMEs in Fig. 3, the proton
(π) core (left plots) shows stronger couplings with selected neutron orbits than the neutron(ν) core (right plots). This can be explained by the effects of the Fock terms [74]. Specifically, the Fock terms arising from the isoscalar σ- and ω-couplings do not contribute to neutron-proton (np ) interactions but present net repulsive neutron-neutron (nn ) and proton-proton (pp ) ones, thus partly canceling their net attractive Hartree terms in thenn andpp channels. In contrast, as shown in the lower plots of Fig. 3, the isovector ρ-T and π-PV couplings present rather strong attractive contributions mainly through the Fock terms, particularly for thenp interactions (left-lower plot). Compared with other RHF and RMF models, thenp interactions are evidently enhanced due to the strong ρ-T coupling in PKA1 [51, 59].Combined with Fig. 2, the IMEs in Fig. 3 can help us understand the evident deformation effects given by PKA1 in Fig. 1. As seen from Fig. 3, the neutron state
ν2s1/2 , namely the1/2+2 orbit atβ=0 , shows rather weak couplings with both neutron and proton cores due to the nodal differences, and thus a loosely bound2s1/2 orbit is found in Fig. 2(b), which indicates an unfavorable even-parity configuration for spherical 11Be. Compared withν2s1/2 , theν1d5/2 state shows much stronger couplings with the cores. Combined with the gradually enlarged1d5/2 proportion and the simultaneously reduced2s1/2 one in Fig. 2 (a), the couplings between the1/2+2 orbit and the cores are then notably enhanced following the enlarged deformation, which eventually leads to a fairly deep bound1/2+2 orbit in the GS of 11Be. In addition, as shown in Fig. 2 (a), the1s1/2 component is also enfolded slightly into the orbit1/2+2 in the GS, which also plays some role in deepening the orbit1/2+2 , due to its rather strong couplings with the core orbits in Fig. 3.In contrast, for odd-parity 11Be, which is not shown in detail, the mixing of the
1p1/2 - and1p3/2 -components also happens in the populated odd-parity orbits1/2−1 and1/2−2 following an increase in the deformation. However, the mixing of high-lying odd-parity components such as2p and1f are hindered by thesd -shell between the p- andpf -shells. This indicates that the sum of the1p3/2 and1p1/2 proportions in these orbits will not change much, which is confirmed by the D-RHFB calculations. As a result, the enhancement of the binding energy from spherical 11Be to the prolate minimum with odd parity is not so remarkable; see Fig. 1 (b). In contrast, the nearby1d5/2 component can be easily enfolded into the even-parity orbit1/2+2 , following the deformation from spherical 11Be to the GS. Thus, due to strong couplings of the orbitals1d5/2 and1s1/2 with the core orbits, more evident enhancement following the deformation is found in even-parity 11Be than that in the odd-parity case, which is crucial for correctly determining the even-parity GS of 11Be. -
On the basis of correctly reproducing the parity of the GS, we further explore the halo phenomenon in 11Be. Figure 4 shows the neutron (N) and proton (P) density profiles of the GS of 11Be given by PKA1, including the logarithm ones along the symmetric z-axis (left plot) and the contour maps (right plot), in which the dashed ellipsoid defined by the deformation β and radius r is shown as a reference for the size and shape of 11Be. As shown in Fig. 4, the neutrons in 11Be are distributed rather extensively along the z-axis. Different from the decoupling picture [42, 75], the neutron dispersion in 11Be shows a consistent shape with the whole nucleus. Moreover, it can be seen from the left plot of Fig. 4 that the neutron dispersion is almost fully a result of the even-parity orbit
1/2+2 (in pattern of horizontal lines), as referred to in the proton density profiles. It is worth noting that such results are consistent with the experimental indication of a one-neutron halo [11, 30–32].Figure 4. (color online) Neutron (N) and proton (P) densities (fm-3) of 11Be given by PKA1, including logarithm ones (left plot) along the z-axis and contour maps (right plot) in which the dashed ellipsoid is defined by the deformation β and radius r.
However, the orbit
1/2+2 in Fig. 2 (b) is somewhat deeply bound for the GS of 11Be, which seems to conflict with the extensive neutron distributions in Fig. 4. Combined with the results in Fig. 2, a picture becomes apparent. As shown in Fig. 2 (b), the neutron orbit2s1/2 is loosely bound for spherical 11Be because of the nodal difference, and for the same reason, thes1/2 fragment (mainly the2s1/2 one) remains loosely bound in the largely deformed GS. Due to the vanishing centrifugal barrier, this leads to the halo occurrence in the GS of 11Be. In contrast, as hindered by the centrifugal barrier, the neighboringd5/2 fragment does not contribute to the halo, but plays an essential role in stabilizing the halo via the attractive inherent correlations with thes1/2 fragment, as indicated by the negativeEc,i value of the halo orbit1/2+2 aforementioned.In fact, the enhancements due to the ρ-T and π-PV couplings are not only crucial for the parity inversion but also significant for the formation and stabilization of the halo in 11Be. Following the deformation from spherical 11Be to the GS, the spherical fragment
2s1/2 remains weakly bound to form the halo, and the intruding fragment1d5/2 in the orbit1/2+2 enhances the couplings with the core orbits of 11Be, resulting in an even-parity GS, in which the ρ-T and π-PV couplings play an important role. This also stabilizes the halo via an attractive inherent correlation with the2s1/2 fragment. Moreover, due to the nodal difference between the2s1/2 and1d5/2 components, such an inherent correlation is not strong enough to hinder halo formation. It is worth noting that this stability mechanism is different from that of halos in spherical unstable nuclei, which are stabilized in general via enhanced pairing correlations by neighboring high-l orbits [36]. -
In summary, employing the D-RHFB model, we explored the mechanism behind the parity inversion and neutron halo of 11Be, which is uniformly interpreted by the RHF Lagrangian PKA1. Using the language of the spherical DWS basis expansion, it is illustrated that evidently enhanced deformation effects by the ρ-T and π-PV couplings are crucial for describing the even-parity ground state and the neutron halo of 11Be correctly. Moreover, a new microscopic picture of the halo formation and stability is revealed in this work, which paves a way to understanding the halo in deformed unstable nuclei. For light nuclei, the effects of the rotation and vibration corrections might be important. In the future, further implementations of the D-RHFB model are expected by considering the angular momentum projection, the number projection, and other beyond mean field effects.
-
The authors want to thank Prof. S.-G. Zhou, Prof. F. Q. Chen, and Dr. X. X. Sun for fruitful discussions and the Supercomputing Center of Lanzhou University and the Southern Nuclear Science Computing Center for the computing resources provided.
