# states in the large boson number limit of the Interacting Boson Approximation model

###### Abstract

Studies of the Interacting Boson Approximation (IBA) model for large boson numbers have been triggered by the discovery of shape/phase transitions between different limiting symmetries of the model. These transitions become sharper in the large boson number limit, revealing previously unnoticed regularities, which also survive to a large extent for finite boson numbers, corresponding to valence nucleon pairs in collective nuclei. It is shown that energies of states grow linearly with their ordinal number in all three limiting symmetries of IBA [U(5), SU(3), and O(6)]. Furthermore, it is proved that the narrow transition region separating the symmetry triangle of the IBA into a spherical and a deformed region is described quite well by the degeneracies , , , while the energy ratio turns out to be a simple, empirical, easy-to-measure effective order parameter, distinguishing between first- and second-order transitions. The energies of states near the point of the first order shape/phase transition between U(5) and SU(3) are shown to grow as n(n+3), in agreement with the rule dictated by the relevant critical point symmetries resulting in the framework of special solutions of the Bohr Hamiltonian. The underlying partial dynamical symmetries and quasi-dynamical symmetries are also discussed.

###### :

21.10.Re, 21.60.Fw, 05.70.Fh6x9

address=Institute of Nuclear Physics, National Centre for Scientific Research “Demokritos”, GR-15310 Aghia Paraskevi, Attiki, Greece

address=Physics Division, Argonne National Laboratory, Argonne, Illinois 60439, USA

address=Wright Nuclear Structure Laboratory, Yale University, New Haven, Connecticut 06520-8124, USA

The Interacting Boson Approximation (IBA) model [1], describing collective phenomena in atomic nuclei in terms of and bosons (of angular momentum 0 and 2, respectively), is known to possess an overall U(6) symmetry, containing three different dynamical symmetries, U(5), SU(3), and O(6), corresponding to near-spherical (vibrational), axially symmetric prolate deformed (rotational), and soft with respect to axial asymmetry (-unstable) nuclei respectively. These limiting symmetries are shown at the vertices of the symmetry triangle [2] of the model, shown in Fig. 1.

An energy functional can be obtained [4] in the classical limit of the model, through the use of the coherent state formalism [5, 6]. Studying this energy functional in the framework of catastrophe theory one can see [7] that a first order phase transition (in the Ehrenfest classification) is predicted to occur between the limiting symmetries U(5) and SU(3), while a second order phase transition is expected between U(5) and O(6). We refer to these transitions as shape/phase transitions. A narrow shape coexistence region is then predicted [4] in the symmetry triangle of the IBA, separating the spherical and deformed phases. The shape coexistence region shrinks into the point of second order phase transition as the U(5)-O(6) line is approached, as shown in Fig. 1.

Shape/phase transitions have been considered recently in the framework of the geometric collective model [8] as well, in which the collective variables and are used. The second order transition between U(5) and O(6) has been described by the E(5) critical point symmetry [9], using in the Bohr Hamiltonian [10] a potential independent of and having the shape of an infinite square well potential in , while the first order transition between U(5) and SU(3) has been described by the X(5) critical point symmetry [11], using in the Bohr Hamiltonian a potential of the form , with having again the shape of an infinite square well potential in , and being a steep harmonic oscillator centered around . E(5) and X(5) are shown in Fig. 1, close to the points of the second and first order phase transitions of the IBA, respectively.

The IBA model predictions have also been studied using different measures of chaotic behavior. It has been found that chaotic behavior prevails over most of the symmetry triangle, with the noticeable exception of two regions, in which highly regular behavior is observed [12]. One of them is located along the U(5)-O(6) line, and its existence is due to the O(5) symmetry known to survive along the whole line [13]. The other regular region is connecting U(5) and SU(3) through a narrow path within the triangle, called the Alhassid–Whelan arc of regularity [12] (also shown in Fig. 1), the symmetry implying its existence being yet unknown.

SU(3) | O(6) | ||||
---|---|---|---|---|---|

Irrep (,) | () | Irrep (,) | () | Irrep () | () |

(2,0) | 0 | () | 0 | ||

(2-4,2) | 1 | (-2) | 1 | ||

(2-8,4) | (4-6)/(2-1) | (2-6,0) | (4-3)/(2-1) | (-4) | 2 |

(2-12,6) | (6-15)/(2-1) | (2-10,2) | (6-10)/(2-1) | (-6) | 3-(3/) |

(2-16,8) | (8-28)/(2-1) | (2-14,4) | (8-21)/(2-1) | (-8) | 4-(8/) |

States with zero angular momentum are of particular interest, since centrifugal effects are absent, facilitating the detection of underlying symmetries. In the U(5) limit of the IBA, the energies of states increase linearly with the number of bosons, corresponding to their phonon number. In the SU(3) limit, the energies of bandheads are determined by the eigenvalues of the second order Casimir invariant of SU(3), the results shown in Table 1. In the O(6) limit, they are determined by the second order Casimir invariant of O(6), the results also shown in Table 1. From Table 1 it is clear that in the limit of large boson numbers , the energies of bandheads in all three dynamical symmetries of the IBA grow linearly, .

IBA numerical calculations have been performed using the recently developed code IBAR [14, 15], which can handle large boson numbers. The standard two-parameter IBA Hamiltonian [16] has been used, depending on the parameters and . Examining the degeneracy , a hallmark of the X(5) critical point symmetry, we find [17] that its locus in the IBA triangle is a straight line, falling withing the coexistence region for large boson numbers, as seen in Fig. 2. Similar results are obtained for the degeneracies and , also shown in Fig. 2.

Plotting the ratio [17], related to the first of the degeneracies mentioned above, we see in Fig. 3 that it exhibits the behavior expected [18] for an order parameter of a first (second) order phase transition for (). It is remarkable that experimental data around the isotones, which are known to be the best empirical examples of X(5) [19, 20, 21, 22, 23], shown in Fig. 4, exhibits a behavior very similar to the one expected for a first order phase transition, while data around Ba, the best example of E(5) [24, 23], shows the behavior expected for a second order transition.

E(5) | Z(5) | X(5) | Norm | IBA-Norm | Z(4) | Norm | X(3) | Norm | |
---|---|---|---|---|---|---|---|---|---|

0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | |

3.03 | 3.91 | 5.65 | 1.0 | 1.0 | 2.95 | 1.0 | 2.87 | 1.0 | |

7.58 | 9.78 | 14.12 | 2.50 | 2.48 | 7.60 | 2.57 | 7.65 | 2.67 | |

13.64 | 17.61 | 25.41 | 4.50 | 4.62 | 13.93 | 4.71 | 14.34 | 5.00 | |

21.22 | 27.39 | 39.53 | 7.00 | 7.13 | 21.95 | 7.43 | 22.95 | 8.00 | |

30.31 | 39.12 | 56.47 | 10.00 | 9.85 | 31.65 | 10.72 | 33.47 | 11.67 |

Model | D | Model | D | ||||
---|---|---|---|---|---|---|---|

E(5) | + | 5 | |||||

X(5) | 5 | Z(5) | 5 | ||||

X(3) | 3 | Z(4) | 4 | 1 |

Now we turn attention to energies of states within the critical point symmetries E(5) and X(5), mentioned above, as well as within Z(5) [25], a similar solution of the Bohr Hamiltonian, also using an infinite square well potential in , with . As seen in Table 2, although the relevant energies look different in each model when normalized to the energy of the state, they become identical if normalized to the energy of [3]. Moreover, they increase as . Similarly, energies of states within Z(4) [26] [similar to Z(5), but with fixed to ] increase as . Within X(3) [27] [similar to X(5), but with fixed to 0], they increase as , as also shown in Table 2. These results can be easily interpreted [3] by taking into account the order of the Bessel functions appearing as eigenfunctions in these models, given in Table 3, as well as the fact that the spectra of the Bessel functions are found to increase approximately as , this result being exact only for [3]. Consideration of the second order Casimir operator of E() [28, 26], the Euclidean algebra in dimensions, shows [3] that the present situation is a partial dynamical symmetry [29] of Type I [30], in which some of the states (the states in the present case) preserve all the relevant symmetry. The recent conjecture [31] of a partial SU(3) dynamical symmetry underlying the Alhassid–Whelan arc of regularity is also receiving attention.

It is a nontrivial result [3] that the IBA near the critical point of the first order transition also yields energies of states increasing as , i.e. in the same way as critical point symmetries based on infinite square well potentials in 5 dimensions (degrees of freedom) predict, as shown in Table 2.

In conclusion, states in the large boson number limit of the IBA exhibit many interesting properties, as well as degeneracies to non-zero angular momentum states, which invite further investigations into determining of the symmetries underlying these regularities.

###### Acknowledgements.

Work supported in part U.S. DOE Grant No. DE-FG02-91ER-40609 and under Contract DE-AC02-06CH11357.## References

- [1] F. Iachello, and A. Arima, The Interacting Boson Model, Cambridge University Press, Cambridge, 1987.
- [2] R. F. Casten, Nuclear Structure from a Simple Perspective, Oxford University Press, Oxford, 1990.
- [3] D. Bonatsos, E. A. McCutchan, and R. F. Casten, Phys. Rev. Lett., accepted.
- [4] F. Iachello, N. V. Zamfir, and R. F. Casten, Phys. Rev. Lett. 81, 1191–1194 (1998).
- [5] J. N. Ginocchio, and M. W. Kirson, Phys. Rev. Lett. 44, 1744–1747 (1980)
- [6] A. E. L. Dieperink, O. Scholten, and F. Iachello, Phys. Rev. Lett. 44, 1747–1750 (1980).
- [7] D. H. Feng, R. Gilmore, and S. R. Deans, Phys. Rev. C 23, 1254–1258 (1981).
- [8] A. Bohr, and B. R. Mottelson, Nuclear Structure, Vol. II, Benjamin, New York, 1975.
- [9] F. Iachello, Phys. Rev. Lett. 85, 3580–3583 (2000).
- [10] A. Bohr, Mat. Fys. Medd. K. Dan. Vidensk. Selsk. 26, no. 14 (1952).
- [11] F. Iachello, Phys. Rev. Lett. 87, 052502 (2001).
- [12] Y. Alhassid, and N. Whelan, Phys. Rev. Lett. 67, 816–819 (1991).
- [13] A. Leviatan, A. Novoselsky, and I. Talmi, Phys. Lett. B 172, 144–148 (1986).
- [14] R. J. Casperson, IBAR code (unpublished).
- [15] E. Williams, R. J. Casperson, and V. Werner, Phys. Rev. C 77, 061302(R) (2008).
- [16] V. Werner, N. Pietralla, P. von Brentano, R. F. Casten, and R. V. Jolos, Phys. Rev. C 61, 021301(R) (2000).
- [17] D. Bonatsos, E. A. McCutchan, R. F. Casten, and R. J. Casperson, Phys. Rev. Lett. 100, 142501 (2008).
- [18] F. Iachello, and N. V. Zamfir, Phys. Rev.Lett. 92, 212501 (2004).
- [19] R. F. Casten, and N. V. Zamfir, Phys. Rev. Lett. 87, 052503 (2001).
- [20] R. Krücken, et al., Phys. Rev. Lett. 88, 232501 (2002).
- [21] D. Tonev, et al., Phys. Rev. C 69, 034334 (2004).
- [22] A. Dewald, et al., Eur. Phys. J. A 20, 173–178 (2004).
- [23] R. F. Casten, and E. A. McCutchan, J. Phys. G: Nucl. Part. Phys. 34, R285–R320 (2007).
- [24] R. F. Casten, and N. V. Zamfir, Phys. Rev. Lett. 85, 3584–3586 (2000).
- [25] D. Bonatsos, D. Lenis, D. Petrellis, and P. A. Terziev, Phys. Lett. B 588, 172–179 (2004).
- [26] D. Bonatsos, D. Lenis, D. Petrellis, P. A. Terziev, and I. Yigitoglu, Phys. Lett. B 621, 102–108 (2005).
- [27] D. Bonatsos, D. Lenis, D. Petrellis, P. A. Terziev, and I. Yigitoglu, Phys. Lett. B 632, 238–242 (2006).
- [28] A. O. Barut, and R. Raczka, Theory of Group Representations and Applications, World Scientific, Singapore, 1986.
- [29] Y. Alhassid, and A. Leviatan, J. Phys. A: Math. Gen. 25, L1265–L1271 (1992).
- [30] A. Leviatan, Phys. Rev. Lett. 98, 242502 (2007).
- [31] M. Macek, P. Stránský, P. Cejnar, S. Heinze, J. Jolie, and J. Dobeš, Phys. Rev. C 75, 064318 (2007).