Due to its flexible shape, the sextic oscillator seems to be a suitable candidate to serve as the
-potential
in the Bohr Hamiltonian. It is the lowest-order polynomial potential that admits two minima. With a potential minimum at
, it can describe spherical nuclear shapes, whereas with a minimum at
, it can account for a nucleus with a deformed shape. Furthermore, the transition between the two shape phases can be implemented smoothly within an exact analytical framework simply by varying the potential parameters
a and
b in (
11). These features have been employed in a variety of models describing a wide range of nuclear shapes and transitions between them. Several review papers have discussed solvable models within the framework of the Bohr Hamiltonian, paying special attention to the shape-phase transitions and critical-point solutions related to them [
19,
30,
39]. In the present review, we focus on models in which the sextic oscillator was considered as the
potential and report on recent developments in the field, placing them in the context of earlier results.
3.1.
The first application of the sextic oscillator was proposed in 2004 in Ref. [
20] to describe the transition between spherical and
-unstable nuclear shapes. In this case, the
shape variable does not play a role. Based on the corresponding critical-point symmetry, this approach can be called
-sextic. The formalism was first developed for solutions with
and 1, encompassing ten states altogether: three
, three
, one
, two
, and one
levels. These model states can typically be assigned to the first four levels of the ground-state band
, the first three levels of the quasi-
band
, and the bandheads of two further
bands.
The formalism of the sextic oscillator was developed in Ref. [
20], where the energy eigenvalues, wavefunctions, and
transitions were derived. The evolution of the energy spectrum was studied by varying the
a and
b model parameters, and the structure of the phase space was established. It was found that the phase transition occurs on the critical parabola
. This follows from selecting
in Equation (
12), which corresponds to
, i.e., to the combinations
,
and
,
. Closed-form benchmark numbers were calculated for the energy and
ratios and were compared to the benchmark numbers of the
critical-point symmetry and other models. As an example, the spectroscopic data for the
Ba nucleus were calculated and compared to the experimental values.
In the next application of the model [
21], the formalism was further developed: matrix elements with the
wave functions were calculated analytically in terms of confluent hypergeometric functions. The model was then systematically applied to the chain of even-even Ru (
98–108), Pd (
102–110), and Cd (
106–116) isotopes. The potential parameters were determined from the experimental energy eigenvalues, and the trajectories of the chains were plotted on the
phase space. Several nuclei were identified as potential critical-point candidates, i.e., those close to the critical parabola. Examples include
Ru,
Pd,
Cd, and
Cd, which were proposed previously as candidates using other models, and
Cd as a new one.
In the applications in [
20,
21], the key step was assigning appropriate experimental states to the excited
model states. These states have characteristic decay patterns in electric quadrupole transitions: the one denoted by
is expected to decay strongly to the second excited
state assigned to
, whereas the node-excited
state
is expected to decay to the first excited
state
(see
Figure 3 for the notation).
Further applications of this version of the model have been presented. In Ref. [
40], seven Te isotopes (from
to 139) and the chain of six Xe isotopes (
122–132) were discussed. All the nuclei were found to lie in the domain of a deformed minimum (above the critical parabola), with
Te closest to the position of a critical-point nucleus. In Ref. [
41], the chains of Zr (
88–100) and Sn (
110–122) nuclei were investigated. Most of these nuclei were also found to belong to the domain of a deformed minimum; however,
Zr and
Sn were proposed as nuclei close to the critical point.
The sextic oscillator model was developed further in Ref. [
22], where the exact analytic formalism was extended to
, i.e., to a higher degree of the QES solution. In this case, the roots of a cubic algebraic equation had to be calculated (see the discussion in
Section 2.2). The considerably larger model space (30 states compared to 10) allowed for a more elaborate comparison with the experimental data but also raised new questions. The enlarged basis included the ground-state band (up to
); states that could be assigned to a
band (up to
); and three excited
states, which act as bandheads. Two of them were analogous to those appearing in the simple version of the model (
and (2,0)0), whereas the (3,0)0 state was a higher nodally excited configuration. Identifying the experimental counterpart of the excited
states remained the key step in the application of the extended model.
Furthermore, a new question arose. In the original version of the model [
20,
21], the states belonging to the same
multiplet were considered degenerate in energy. The pattern of states with larger
L values appearing in the extended model revealed that the spectrum did not follow this degeneracy; rather, typically, the levels with higher
L appeared higher in the spectrum. For this reason, in Ref. [
22], the Hamiltonian was extended by a phenomenological rotational term
. The justification for this rotational term was presented in Ref. [
22].
In the first application of the extended model [
22], the Ru and Pd chains were revisited. With the larger
M, the
numbers (in Equation (
12)) also changed (
was used instead of
), which slightly modified the location of the critical parabola. The conclusions for the Ru chain were essentially the same as in the limited model: the critical point was again close to
Ru. The situation was somewhat different for the Pd chain, perhaps due to the enlarged model space and the rotational term. Most of the Pd nuclei were located in the domain of the spherical minimum (below the critical parabola); however,
Pd, which had a rich enough collection of
values for comparison with the experiment, showed marked signs of a critical-point nucleus. The results for the energies and
B(
E2) in
Pd are shown in
Figure 4. Besides this nucleus,
Pd was also found to be located close to the critical parabola.
In the next application [
23] of the extended model, Pt and Os nuclei were studied close to the closure of the
shell. A thorough analysis of the decay pattern of the excited
states was carried out in order to assign the model states to their experimental counterparts. It was found that the members of the node-excited band built on the
state systematically became lower with increasing
A, and eventually, the (2,0)0 bandhead state fell below the (1,3)0 bandhead state for
Pt. This is shown in
Figure 5.
This change was accompanied by a drastic change in the potential shape, indicating the realization of the phase transition from the
-unstable to the spherical domain. Similar trends were also observed for the Os chain. However, the last isotope with a reasonably comprehensive spectroscopic data set,
Os, remained clearly in the domain of a
-unstable deformed shape. Indirect signs indicated that the shape-phase transition could occur near
Os, but the spectroscopic information was insufficient to prove this unambiguously. This study also confirmed that
Pt and
Pt were rather good examples of the
-unstable configuration (see
Figure 6, taken from Ref. [
23]), whereas this was less pronounced for the Os isotopes, which, especially the lighter ones, showed signs of a triaxial configuration.
There have been further studies estimating the performance of the sextic oscillator as a
potential by abandoning exact solutions. In Ref. [
32], higher QES solutions were considered using numerical calculations. A comparison with a large array of experimental spectroscopic data for the nuclei discussed in Ref. [
21] indicated that the order
seems to be the optimal choice in most calculations. This corresponds to
, i.e., the extended scheme with
M up to 2 presented in Refs. [
22,
23], providing strong support for this fully analytical model.
Ref. [
42] focused on an alternative approach to addressing the problem of degeneracies in the sextic oscillator model. This problem was addressed in Refs. [
22,
23] by directly introducing a phenomenologic rotational term in the Hamiltonian. In Ref. [
42], a rotational term
was directly incorporated into the
potential to resolve the degeneracy of the multiplets. Numerical methods were used to calculate the spectroscopic quantities, and the model was used to reproduce experimental data for a number of Xe isotopes. The results showed good agreement with the data, especially in the case of the staggering pattern appearing in the quasi-
bands.
3.2. w Non-Zero
As discussed in
Section 2, exact solutions for a
potential can be obtained if the potential can be separated into
and
components. This requires certain assumptions and approximations. To obtain a solvable problem in the
variable, further approximations are necessary. Generally, it is assumed that
has a minimum at either
or
, and the potential is expanded around these minima. The two choices lead to different solutions and different physical systems, i.e., different nuclear shape configurations.
When the expansion of
is carried out around
, the Bohr Hamiltonian describes nuclear shapes that evolve from an axially symmetric shape to a triaxial one. By combining this system with a
potential that allows a spherical or deformed minimum, phase transitions from the spherical to the axially deformed shapes can be described with an
critical-point symmetry [
7]. A key development in this area was presented in Ref. [
27], where the sextic oscillator was used as
), while the
potential was defined as
. With an appropriate variable change, the solutions of the differential equation obtained for the
component were written in terms of the spheroidal functions. This model was coined the sextic and spheroidal approach (SSA) and was later also referred to as
-sextic. Calculations were carried out for nuclei typically considered candidates for the
symmetry: Os (
176–180, 188–190),
Nd,
Dy, and Hf (
166–68). The predictions were compared to those of alternative models with different
used in the Bohr Hamiltonian, i.e., the infinite square well (ISW) and the Davidson (D) potential, and also to those of the coherent state model (CSM), which is based on a boson operator formalism. It was found that the best agreement with the experimental energy eigenvalues was found for the D, SSA, and CSM models, whereas for the
values, the SSA and the CSM performed slightly better. Among the nuclei considered,
Os was found to come closest to the critical point of the phase transition.
By applying the expansion of the
around
, transitions between the prolate and oblate shape phases can be described, and the corresponding critical-point solution is
[
9]. The sextic oscillator was applied in this setting as a
, while the
potential was defined as
[
26]. Solutions to the differential equation of the
component were obtained after transforming the latter into the Mathieu equation. This model was referred to as the sextic and Mathieu approach (SMA) or
-sextic. Predictions of the model were calculated for five isotopes—
Os,
Os,
Os,
Th, and
Th—which partially overlap with the examples studied in the
-sextic model. Since transitions among well-deformed nuclei (from prolate to oblate) through a triaxial critical point were studied, regions where this transition was expected, such as Os and Th isotopes, were considered. Predictions for the energy spectrum and the electric quadrupole transitions in the ground,
, and
bands were compared to the experimental values, as well as to the results obtained from the coherent state model (CSM). Signatures for triaxiality were analyzed in the energy spectra. Reasonable agreement was found between the experimental and the SMA and CSM values. A similar analysis (with similar results) was carried out for the nuclei
Hf and
W [
43]. Besides the application of the SMA to certain nuclei, this work also focused on exploring the relation between the SMA and the CSM. It was found that the SMA formalism could be obtained by quantizing the classical CSM equations, and the SMA represented the strong coupling limit of the CSM.
3.3. The -Rigid Cases
A special kind of solvability arises when the variable is “frozen” to some value, typically or , i.e., the minimum values of in the and models. This ensures the separation of variables while reducing the number of active variables.
As described in
Section 2.1, when
, the
model reduces to
—the model that can account for axially symmetric
-rigid shapes. Combining this with the sextic oscillator as
, the
-sextic model arises [
29]. Predictions for some characteristic energy ratios were calculated, and compared to the corresponding values obtained for polynomial versions of the
model:
,
-
,
-
, and
-
. A simplified version of this model, without a quartic term, was considered in Ref. [
44] and compared to the same models. The energy spectrum of the ground and the first two excited
bands was numerically calculated for a large number of nuclei: Ru (
98–108), Mo (
, 102), Xe (
116–130), Ce (
, 134), Nd (
146–150), Sm (
, 152), Gd (
, 154), Dy (
, 156),
Os, Pt (
180–196),
Hg, and
Ra. The energy eigenvalues and
ratios were reproduced reasonably well. A first-order shape-phase transition between a
-rigid prolate harmonic vibrator and
-rigid prolate anharmonic vibrator was identified, and candidates for the corresponding critical point were proposed. The best candidates were found to be
Ru,
Xe,
Xe, and
Nd, while
Xe,
Os, and
Pt were also close to the critical point. We note that
Ru was also considered a critical-point nucleus in the spherical to
-unstable shape-phase transition, i.e., in the
-sextic model.
The
-sextic model has been further developed along several lines. In Ref. [
45], the minimal length (ML) formalism was incorporated into the model, whereas in Ref. [
46], the deformation-dependent mass (DDM) formalism was considered. As a further option, in Ref. [
47], the order of the QES solutions was varied, and the results were compared. Calculations were carried out in all three studies for the same set of nuclei as discussed in Ref. [
29]. The quality of the fits was generally improved due to the extra parameters, but the conclusions concerning the critical-point nuclei essentially remained the same.
Another
-rigid model can be obtained from
-sextic by fixing
to
[
28], leading to
-sextic. In this case, exact solutions could be obtained for the ground and the
bands, whereas for the
band, approximations had to be applied. Numerical calculations were carried out for the Xe (
128-132) and Pt (
192-196) nuclei. The results for the energy spectrum and the
transitions were compared to the experimental values, as well as to the predictions of the
model. Qualitative agreement of the data sets was found, and a shape-phase transition was identified near
Xe, where a deformed minimum gave way to a spherical one. It is notable that these regions were also identified as candidates for phase transition in the
and
models.
3.4. Unrestricted Sextic Oscillator
The QES formalism restricts the parameters of the sextic oscillator, meaning that the allowed potential shapes are also restricted. It is known, for example, in [
20,
21,
22], that in the case of a double-well structure (left of the critical parabola in
Figure 2), the deformed minimum is always deeper than the spherical one. Double-well structures different from this one can be studied for the general sextic oscillator; however, in this case, the exact analytical calculations have to be replaced with numerical approaches.
In Ref. [
31], the general double-well sextic potential was used as a
potential to describe
-unstable nuclei. The model states were constructed using a numerical diagonalization procedure in a basis classified by the
group chain. The model parameters were determined by fitting the low-lying energy spectra of Mo isotopes (
, 98, and 100). Reasonable agreement was found between the theoretical and experimental energy eigenvalues. The difference between the energy minima showed a decreasing trend with an increasing mass number, similar to the height of the potential barrier between the minima. The
values calculated with the resulting wave functions also showed reasonable agreement with the experimental values.
Mo was found to exhibit properties expected near the critical point.
A similar study was carried out in Ref. [
48] within the
formalism to describe the first-order phase transition between spherical and axially deformed nuclear shapes. The
potential was a one-parameter double-well sextic oscillator with two degenerate minima (one spherical and one deformed). The calculated spectrum, as a function of the free parameter, exhibited predictions in accordance with the
and
-
models in certain limits of the parameter. Depending on the relative energy of the barrier compared to the low-energy states, the study identified cases of coexisting shapes in the ground state and the
-excited state. The paper discussed critical situations of shape coexistence from a phenomenological perspective, employing the density of probability distribution for the deformation and
transitions. Experimental evidence supporting these critical phenomena was mentioned for a few selected nuclei (
Nd,
Hf, and
Pu). Similar critical potentials with shape coexistence are expected for higher excited states. In addition, in Ref. [
49], along the same lines, the sextic potential was used to show the connections between shape coexistence and shape transitions.
Lastly, it is worth highlighting the connection between the sextic potential and Catastrophe Theory (CT) [
50]. This theory provides a framework for understanding abrupt transitions in various natural and artificial systems. It offers a unique perspective on the dynamics of complex systems and the potential for abrupt, transformative changes within them. CT is relevant in multiple fields, including science, economics, psychology, engineering, ecology, etc. Specifically, the sextic potential is associated with the A
butterfly catastrophe, whose germ is represented by the function
. In this context, the potential used is typically expressed as
, where
and
are free parameters. Importantly, for
-independent potentials, this setup not only accounts for transitions from spherical to deformed shapes through a second-order phase transition but also accommodates the critical point at which the minima corresponding to spherical and deformed shapes become degenerate, giving rise to a first-order phase transition. In a recent publication [
51], this relationship was explored within the context of an algebraic approach, leading to the derivation of the general phase diagram for the butterfly catastrophe. This comprehensive phase diagram includes analytical expressions for key features such as the tricritical point, the critical lines for both first- and second-order transitions, and additionally, the delineation of spinodal and anti-spinodal lines that demarcate the regions of coexistence between spherical and deformed shapes. This connection to Catastrophe Theory provides valuable insights into the complex dynamics of nuclear shape transitions.
3.5. Further Generalizations
The sextic oscillator is not the only quasi-exactly solvable potential. In fact, polynomial potentials of degree
belong to this class. Following the sextic oscillator (
), the decadic oscillator (
) also admits similar mathematical formalism. The coefficients of the potential terms also obey certain constraints. In Ref. [
52], the decadic oscillator was employed to describe nuclei that were supposed to lie near the
critical point, characterizing the phase transition from the spherical to the
-unstable shape phases. The spectroscopic results were compared to the experimental data and the results obtained from various models (
,
-sextic, etc.). It was found that the quality of the fit to the experimental energy eigenvalues (in terms of
) was better for the sextic oscillator and the decadic oscillator for three and two nuclei, respectively (
Pd,
Ru, and
Cd vs.
Cd and
Ba), whereas it was of similar quality for one case (
Pd). In the case of six Xe isotopes, the performance of the decadic oscillator was better for four nuclei. In all cases, both models performed better than the original
model based on the infinite square as well as
. However, for the
ratios, the
model almost always performed best, and the sextic oscillator was the worst. This may be due to the difference in the asymptotic tail of the wave functions. It is important to note that higher-order polynomial potentials lead to more intricate recurrence relations and, consequently, more complex truncation conditions. Exact analytical solutions for such potentials cannot be attained as straightforwardly, as demonstrated in Equation (
12). This significantly reduces the number of algebraic analytical solutions. With regard to these exact solutions, due to the intricacies involved in truncating the series expansion, interested readers should be mindful of the complications outlined in Ref. [
53].
Another generalization of the sextic oscillator was discussed in Ref. [
54]. The authors discussed odd-
A nuclei in the
-soft region in a model where the
-sextic Hamiltonian is combined with a single nucleon on a
orbit. In the resulting
model, the single-particle features are represented by a five-dimensional spin–orbit interaction and an angular momentum degeneracy-breaking term proportional to
. The energy spectrum and the
rates were precisely calculated for odd Ir isotopes with
to 195. Good agreement was found between the experimental data and the model predictions.