Nuclear Structure Study of Odd-Odd Yttrium Nuclei within Interacting-Boson Fermi-Fermion Model (IBFFM)


Nuclear structure of odd-mass Yttrium nuclei is studied within the framework of the interacting boson-fermion model (IBFM) with parameters based on the optimization method. In addition to the energy surfaces, both single particle energies and occupation probabilities have been used as a microscopic input for building the IBFM Hamiltonian. Only three strength parameters for the particle-boson-core coupling are fitted to experimental spectra. The IBFM Hamiltonian is then used to compute the energy spectra and electromagnetic transition rates for selected odd-mass Sr nuclei. A reasonable agreement with the available experimental data is obtained for the considered odd-mass nuclei.

Share and Cite:

Mohaisen, A. and Abood, S. (2024) Nuclear Structure Study of Odd-Odd Yttrium Nuclei within Interacting-Boson Fermi-Fermion Model (IBFFM). Journal of Applied Mathematics and Physics, 12, 2020-2031. doi: 10.4236/jamp.2024.126123.

1. Introduction

The neutron deficient nuclei originating from the sector characterized by N and Z approximately equal to 38 to 40 have garnered sustained attention in recent years. Within this domain, the shell model orbitals exhibit discontinuities at significant deformations for nucleon numbers 38 and 40 [1] . Moreover, given that both neutrons and protons populate identical orbitals in this realm, the impacts of such discontinuities are notably pronounced. An “island” of nuclei characterized by a significant prolate deformation in their ground state (with a deformation parameter β2 of approximately 0.4) has been identified at the proton and neutron numbers of 38 and 40, respectively [2] [3] . Furthermore, owing to the relatively low density of shell model orbitals in these nuclei, even minor variations in the number of nucleons lead to remarkable alterations in their shapes. In a similar manner, a correlation can be observed when the spin is enhanced, as various occurrences like alignment of particles and hindrance in their motion are linked with diverse modifications in shape. Furthermore, an additional area of significance within this vicinity adjacent to the N = Z boundary involves the examination of the impact stemming from the pairing of neutrons and protons, with a particular focus on the rivalry between its components characterized by isospin values of T = 0 and T = 1.

The nuclei within this isolated locality generally exhibit relatively small cross-sectional areas in heavy-ion-induced reactions achievable with stable targets and projectiles. For odd-odd nuclei, the challenge is further compounded by their intricate level structures. Hence, when investigating such nuclei, it becomes crucial to employ the IBFFM model.

In the current study, we provide findings from gamma-ray spectroscopy performed in-beam for the odd-odd nucleus 80Y, situated closely to the highly deformed A ≈ 80 nuclei. Previous outcomes related to various rotational bands observed in this nucleus have been documented in Reference [4] . With the introduction of novel measurements, the high-spin level scheme has been fully explored and expanded. Recent supplementary investigations [5] [6] [7] [8] have yielded significant insights into the lower energy segment of the level scheme. The analysis of the detected formations in 80Y is deliberated upon utilizing the interacting boson-fermion-fermion model (IBFFM) and the cranked shell model (CSM).

2. Interacting Boson Fermion-Fermion Model

The model of interacting boson-fermion-fermion is an expansion designed for odd-odd nuclei from the model of interacting boson (IBM) [9] (pertaining to even-even nuclei) and the model of interacting boson-fermion (IBFM) [10] (for odd-mass nuclei). The Hamiltonian of IBFFM is expressible as follows [11] .

H I B F F M ^ = H I B F M ^ ( π ) + H I B F M ^ ( ν ) + H R E S ^ ( π ν ) H I B M ^ (1)

Here, the acronym HIBM is utilized to denote the typical IBM Hamiltonian applicable to even-even core systems [9] . The symbol H π ( ν ) is employed to characterize the interaction involving fermion single-particles, while H B π ( B ν ) is used to elucidate interactions between bosons and fermions. The mathematical expression (1) can be restated as H = H I B F M ( ν ) + H I B F M ( π ) H I B M + H π ( ν ) , wherein H I B F M ( ν , π ) symbolizes the Hamiltonians in the Interacting Boson Fermion Model (IBFM) for nuclei with odd-mass configurations [11] . The final component accounts for a residual interaction between neutrons and protons.

In this particular model, the 80Y nucleus is delineated as a composite system comprising two fermions: one proton fermion and one neutron-hole fermion, interacting with an 80Sr core. The scope of this model is restricted to the depiction of states involving two quasiparticles. The initial phase involves the determination of the phenomenological parameters in the first five terms of Equation (1) through the characterization of the 80Sr core nucleus and the two neighboring odd-mass nuclei, namely 81Y and 79Sr. The Hamiltonian (Equation (1)) will subsequently undergo diagonalization within a limited space formed by linking proton and neutron states to the core space of bosons. The parameters from Ref. [12] were employed for the Sr core within the IBM framework. In the case of the positive-parity states in 81Y, the IBFM parameters established for the odd-mass Y isotopes in Ref. [13] were utilized for the analysis. In the present calculations, the eccentric proton was permitted to populate the orbitals 1g9/2 and 2d5/2 within the framework of the spherical shell model. Extensive information can be found in the references provided (Refs. [12] [13] ). Subsequently, both states of positive and negative parity have been computed for the nucleus 79Sr, following a methodology akin to the one expounded in reference [13] . The odd fermion was permitted to inhabit the spherical orbitals 1f5/2, 2p3/2, and 2p1/2 for the negative-parity states. The boson-fermion interaction terms parameters were A0 = −0.11, Γ0 = 0.29, and Γ0 = 8.20 MeV2, while the values for the negative-parity states were −0.11, 0.055, and 8.20, respectively. These parameters determine the strengths of the boson-fermion interaction’s exchange, quadrupole-quadrupole, and monopole-monopole terms, respectively [9] . The codes PHINT and FBEM [14] , and ODDA and PBEM [15] , for even-even and odd-mass nuclei, respectively, have been used to calculate the energy levels and electromagnetic decay rates. Scholten IBFFM algorithm has been used to conduct the computations for the odd-odd nucleus [16] . The program’s description may be found in Ref. [17] , and for information on how the Hamiltonian (Equation (1)) is diagonalized, we refer to this work. The neutron-proton residual interaction in this code is made up of a spin-spin term and a quadrupole-quadrupole interaction:

H π ν = V q Q π Q ν + V s σ π σ ν (2)

The IBFFM code includes the residual interactions spin-spin, multipole-multipole, spin-spin-delta, and tensor interaction, which are represented as [11] :

H R E S ^ π ν = H δ ^ + H σ σ σ ^ + H σ σ ^ + H T ^ + H M M ^ (3)


H δ ^ = 4 π V δ ( r π r ν ) δ ( π r R 0 ) δ ( r ν R 0 ) (4)

is a surface delta-function interaction between the proton and neutron,

H σ σ σ ^ = 4 π V σ σ σ ( r π r ν ) ( σ π σ ν ) δ ( π r R 0 ) δ ( r ν R 0 ) (5)

is a surface contact spin-spin interaction,

H σ σ ^ = 3 V σ σ σ π σ ν (6)

is a spin-spin interaction,

H T ^ = V T [ 3 ( σ π r π ) ( σ ν r ν ) r π ν 2 σ π σ ν ] (7)

is a tensor interaction, and

H M M ^ = 4 π δ ( r π r ν ) r π r ν k m V k Y k m * ( θ π , ϕ π ) Y k m ( θ ν , ϕ ν ) (8)

the residual interaction Hamiltonian was calculated as follows:

H R E S ^ = 4 π V δ δ ( r π r ν ) δ ( r π R 0 ) 3 V σ σ ( σ π σ ν ) + V T [ 3 ( σ π r π ) ( σ ν r ν ) r π ν 2 σ π σ ν ] (9)

where R 0 = 1.2 A 1 / 3 , the parameters V δ , V σ σ σ , V σ σ , V T and V R are constants, and determined the strengths of the interactions.

3. Results and Discussion

We took into consideration negative-parity configurations in 79Sr and positive-parity configurations in 81Y for the negative-parity ones we were able to calculate the positive-parity states in 80Y by permitting positive-parity configurations in both linked odd-mass nuclei.

The other two potential combinations, which consist of configurations with negative parity in 81Y, were not examined, as it is anticipated that they would contribute to elevated excitation levels. Calculation of B(E2) transition probabilities employed an operator that was deliberately selected to be identical to the one utilized in characterizing the E2 characteristics of the even-even core and the two corresponding odd-mass isotopes. Conversely, the M1 characteristics were assessed using an operator that bore resemblance to the one employed for the odd-mass isotopes [17] .

Having established the IB(F)M parameters for the core and odd-mass nuclei, the computations pertaining to the odd-odd nucleus were uncomplicated, necessitating only the selection of the neutron-proton interaction parameters. The impact of the quadrupole-quadrupole interaction was deemed minimal on the outcomes, ultimately leading to its exclusion. It is comprehensible, given that the majority of quadrupole interactions have already been addressed by the interaction between bosons and fermions. The interaction between spins has implications on the arrangement of the lowest energy states; however, a trade-off was necessary to describe both positive and negative states effectively, requiring very low values for Vs. Consequently, the final outcomes were minimally affected by this parameter. Here, we present the results obtained using Vs = 0.03 MeV.

Figure 1 and Figure 2 depict the outcomes of the computations for the positive-parity states, reaching spin 14, and the negative-parity states, extending up to spin 15, respectively. The derived states have been organized into groupings that are tentatively linked to those observed in experimental settings. The correlation between the computed states and the ones witnessed has been established by examining their decay mechanisms. The computed branching ratios are juxtaposed with the empirical ones in Table 1, revealing that the computations satisfactorily depict the data. Evident from Table 1 is the assignment of the 5+ state

Figure 1. Comparison of the experiment [18] for the positive-parity states with the predictions from IBFFM.

Figure 2. Comparison of the experiment [18] for the negative-parity states with the predictions from IBFFM.

Table 1. A comparison is made between the IBFFM and experimental [18] branching ratios in 80Y. The initial section of the table illustrates the in-band transitions, while the subsequent part delineates the interrelations among the lowest states of the positive-parity bands. In the case of each particular level, the most robust branch is standardized to 100.

emerging from band 5 as the second computed ( 5 2 + ) state. The calculated 5 1 + state could possibly correspond to the level situated immediately below the 6+ state of band 3. The three states positioned beneath the 6+ state are influenced by the feeble transitions of 0.383, 0.232, and 0.186 MeV, thus complicating their proper characterization, with uncertainties even surrounding their sequential arrangement. However, supposing that the state at 0.675 MeV (which experiences decay from the 1.059 MeV, 6+ state through the 383 keV transition) represents the 5 1 + state, it becomes apparent that the 6 1 + state primarily decays into this particular state, whereas the 7 1 + and 6 2 + states (belonging to band 5) exhibit minimal connections to this state, a phenomenon corroborated by experimental observations.

The IBFFM wave-functions for the positive-parity states are dominated by the g9/2 proton and neutron orbitals, while in the negative-parity states the dominating orbitals are πg9/2 and νf5/2. In some cases, this confirms the Nilsson assignments made on the basis of the two-qp-plus-rotor model. On the other hand, the present calculations cannot provide a strong contribution of the d5/2 orbital to the wave-functions. The d5/2 orbital in a spherical shape is positioned at a considerable distance in comparison to the g9/2 orbital. Consequently, it becomes involved in the wave-functions of odd-A isotopes with contributions that are less than 10%. As a result, we anticipate that the depiction of the well deformed scenario will not be very accurate, particularly when the [431] 1/2 orbital (originating from d5/2) significantly competes.

Figure 3 illustrates the prediction of the B(M1)/B(E2) ratio for the designated bands. The estimation indicates an overestimation by a factor of 3 - 4 at the lower spins within band 5, and an underestimation for band 6. Notably, the contribution of the [431] 1/2 orbital is deemed significant for both bands. In addition,

Figure 3. IBFFM calculations and experimental data [18] for B(M1)/B(E2) for the bands in 80Y. The dashed-dotted line represents the IBFFM forecasts for each band, and the cranked shell model (CSM) is used to compute the other curves for the given combinations.

the B(M1) strengths are also projected to be underestimated for the odd-spin members of the positive-parity band 3, a trend reflected in the predicted branching ratios presented in Table 1.

The limitations of these computations seem to lie in their incapacity to precisely delineate the relative positioning of band-heads and the spacings of the lowest levels within the bands. Enhanced positioning accuracy of negative and positive-parity band-heads would have necessitated adjustments in the parameter Vs. Conversely, while there could have been some enhancement in the spacing of the lowest states, its impact on the computed decay characteristics was minimal. These deficiencies may be ascribed to the insufficiency of the neutron-proton interaction model employed. The significance of this interaction in elucidating intricacies of the observed level structure was underscored in both the two-qp-plus-rotor model [6] and in other IBFFM computations. Particularly in the IBFFM calculations, superior outcomes have been achieved by incorporating a surface-δ interaction [11] . Due to the limitation of the utilized code to only incorporate interaction, the examination of such effects was precluded. Further scrutiny was conducted on our IBFFM computations pertaining to the g 9 / 2 g 9 / 2 band, which was also identified in the more massive odd-odd Y isotopes, 82Y [19] [20] [21] and 84Y [22] [23] . Illustrated in Figure 4 is the manner in which this band is delineated in 82Y and 84Y, employing calculations akin to those expounded upon earlier for 80Y. The fundamental characteristics of these bands are effectively expounded upon until the point of backbending (interaction with other bands, as elaborated below). While not perfectly in phase with the observed data, the computations also exhibit an inversion of the signature splitting at lower angular momenta, and at higher angular momenta (exceeding 9ħ for the lighter isotopes, and 11ħ for 84Y), the signature splitting is accurately depicted.

4. Conclusion

High-spin states have been identified in the neutron-deficient nucleus 80Y through IBFFM. The unambiguous association of prompt gamma rays with this nucleus has been achieved through their calculations within this model. The existing level arrangement [4] has been expanded by the inclusion of fresh bands and transitions. Numerous band configurations exhibiting rotational properties have been detected up to angular momenta of at least 17, suggesting a nucleus characterized by significant quadrupole deformation. The scrutinized band structures are deliberated based on the interacting boson-fermion-fermion model,

Figure 4. IBFFM predictions for the positive-parity band 5 and its equivalents in 82Y [19] [20] [21] [22] and 84Y [22] [23] are compared in an experiment.

thereby facilitating identifications of their principal configuration. The determination of absolute electromagnetic decay transition probabilities (level lifetimes), which can now be achieved using the most extensive Ge detector arrays, would be crucial in enhancing the characterization of this nucleus structure and the progression of its shapes as a function of rotational frequency.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.


[1] Nazarewicz, W., Dudek, J., Bengtsson, R., Bengtsson, T. and Ragnarsson, I. (1985) Microscopic Study of the High-Spin Behaviour in Selected A ≃ 80 Nuclei. Nuclear Physics A, 435, 397-447.
[2] Lister, C.J., Ennis, P.J., Chishti, A.A., Varley, B.J., Gelletly, W., Price, H.G. and James, A.N. (1990) Shape Changes in N = Z Nuclei from Germanium to Zirconium. Physical Review C, 42, R1191(R).
[3] Lister, C.J., Campbell, M., Chishti, A.A., Gelletly, W., Goettig, L., Moscrop, R., Varley, B.J., James, A.N., Morrison, T., Price, H.G., Simpson, J., Connel, K. and Skeppstedt, O. (1987) Gamma Radiation from the N = Z Nucleus 8040Zr40. Physical Review Letters, 59, 1270-1273.
[4] Bucurescu, D., Ur, C.A., Bazzacco, D., Rossi Alvarez, C., Spolaore, P., Petrache, C.M., et al. (1995) Rotational Structures in the Odd-Odd Nucleus 80Y. Zeitschrift für Physik A, 352, 361-362.
[5] Regan, P.H., Chandler, C., Pearson, C.J., Blank, B., Grzywacz, R., Lewitowicz, M., et al. (1997) Spectroscopy of Odd-Odd Nuclei within the Interacting Boson-Fermion-Fermion Model Based on the Gogny Energy-Density Functional. Acta Physica Polonica B, 28, 431.
[6] Döring, J., Schatz, H., Aprahamian, A., de Haan, R.C., Göres, J., Wiescher, M., et al. (1998) New Isomer in 80Y. Physical Review C, 57, 1159.
[7] Chandler, C., Regan, P.H., Blank, B., Pearson, C.J., Bruce, A.M., Catford, W.N., et al. (2022) Levels in 130Cs and IBFFM. Physical Review C, 61, Article ID: 044309.
[8] Ressler, J.J., Walters, W.B., Grzywacz, R., Batchelder, J.C., Bingham, C.R., Gross, C.J., et al. (2023) IBFFM Yrast Band in Odd-Odd Nuclei Associated with O(6) and SU(3) Limits. Physical Review C, 63, Article ID: 067303.
[9] Iachello, F. and Arima, A. (1987) The Interacting Boson Model. Cambridge University Press, Cambridge.
[10] Iachello, F. and Scholten, O. (1979) Interacting Boson-Fermion Model of Collective States in Odd-A Nuclei. Physical Review Letters, 43, 679-682.
[11] Paar, V. (1984) Extended Parabolic Rules for Odd-Odd Nuclei. In: Dombrádi, Z. and Fenyes, T., Eds., In-Beam Nuclear Spectroscopy, Vol. 2, Akadémiai Kiadó, Budapest, 675.
[12] Bucurescu, D., Cata, G., Cutoiu, D., Constantinescu, G., Ivascu, M. and Zamfir, N.V. (1983) Description of the Neutron Deficient Sr and Zr Isotopes in the Interacting Boson Model. Nuclear Physics A, 401, 22-40.
[13] Bucurescu, D., Cata, G., Ivascu, M., Zamfir, N.V., Liang, C.F. and Paris, P. (2020) Collectivity of the Neutron-Deficient Odd Yttrium Isotopes. Journal of Physics G: Nuclear Physics, 14, L175.
[14] Scholten, O. (1979) Computer Codes PHINT and FBEM. KVI Report No. 63.
[15] Scholten, O. (1982) Computer Codes ODDA and PBEM. KVI Report No. 252.
[16] Scholten, O. (1988) IBFFM Code.
[17] Chou, W.-T., McHarris, W.C. and Scholten, O. (1988) Deformed Odd-Odd 180,182,184Re Isotopes in the Interacting-Boson Model. Physical Review C, 37, 2834-2851.
[18] Singh, B. (2005) Nuclear Data Sheets for A = 80. Nuclear Data Sheets, 105, 223-418.
[19] Womble, P.C., Döring, J., Glasmacher, T., Holcomb, J.C., Johns, G.D., Johnson, T.D., et al. (1993) Nuclear Structure of Odd-Odd 82Y. Physical Review C, 47, 2546-2559.
[20] Paul, S.D., Jain, H.C., Chattopadhyay, S., Jhingan, M.L. and Sheikh, J.A. (1995) High-Spin Structure of 82Y. Physical Review C, 51, 2959-2972.
[21] Mukai, J., Odahara, A., Tomura, H., Suematsu, S., Mitarai, S., Kuroyanagi, T., et al. (1994) High-Spin States in the Odd-Odd Nucleus 82Y. Nuclear Physics A, 568, 202-220.
[22] Chattopadhyay, S., Jain, H.C., Sheikh, J.A., Agarval, Y.K. and Jhingan, M.L. (1993) Rotational Structure and Signature Inversion in Odd-Odd 84Y. Physical Review C, 47, R1(R).
[23] Chattopadhyay, S., Jain, H.C., Paul, S.D., Sheikh, J.A. and Jhingan, M.L. (1995) High Spin Structure of 84Y. Physical Review C, 49, 116-130.

Copyright © 2024 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.