Book Title
Abstract
A general form of the octupole collective Hamiltonian is introduced and analyzed based on fundamental tensors in the sevendimensional tensor space. Possible definitions of intrinsic frames of reference possessing cubic symmetry for the octupole tensor are considered. Cubic intrinsic octupole coordinates or deformations are introduced. Shapes of the octupoloid are investigated. The octupole collective Hamiltonian is expressed in intrinsic coordinates. An intrinsic angular momentum carried by the octupole vibrations is discovered. Small oscillations about an axiallysymmetric pear shape are analyzed. Formulation of a unified quadrupoleoctupole collective model is discussed.
Chapter 1 The octupole collective Hamiltonian.
Does it follow the example of the quadrupole case?
Institute of Theoretical Physics, Faculty of Physics, University of Warsaw, Pasteura 5, 02093 Warsaw, Poland
StanislawG.Rohozinski@fuw.edu.pl
S.G.Rohoziński and L. Próchniak]Leszek Próchniak
Heavy Ion Laboratory, University of Warsaw,
Pasteura 5A, 02093 Warsaw, Poland
prochniak@slcj.uw.edu.pl
1 Introduction
The idea of attributing a definite multipolarity to nuclear collective excitations came from the nuclear liquiddrop model. It was known before the discovery of atomic nuclei that the normal modes of the surface vibrations of a spherical drop of incompressible liquid have definite multipolarities [1]. Many years later, when atomic nuclei were already known, Siegfried Flügge [2] connected the eigenfrequencies of the surface vibrations of a nuclear liquid drop with the excitation energies of lowlying states of eveneven nuclei. These were the beginnings of the collective model.
The most important collective mode in nuclear structure physics is that with multipolarity , because it concerns the lowest excited states in eveneven nuclei. These states are well known experimentally and theoretical methods describing them are well developed. One such method, applicable to eveneven nuclei only, is a description by means of the Schrödinger equation in a collective space. Initiated by Aage Bohr [3] a long time ago, it is very effective and is still in use. The Bohr Hamiltonian became the label of the model. When applying the method one usually starts with a general classical collective Hamiltonian [4, 5]. Its form has been extracted from various microscopic manybody models through methods of the “Adiabatic TimeDependent HartreeFockBogolyubov” (ATDHFB) type (see e.g. Chapt. 12 in Ref.[6]) and then quantized by means of the PodolskyPauli prescription [7, 8, 9]. Such a quasiclassical approach is still used even today. It would be good to use purely quantum methods to describe the collective states. Therefore, the Generator Coordinate Method (GCM) has been recently employed [10, 11]. It leads to integral equations instead of differential ones (see e.g. Chapt. 10 in Ref.[6]). However, in an approximation to the GCM called the Gaussian Overlap Approximation (GOA) the HillWheeler integral equations for collective motion can be reduced back to differential equations with a collective Hamiltonian [12].
Apart from the positiveparity quadrupole states, negativeparity levels have been observed for a long time [13] in the lowenergy spectra of eveneven nuclei. Multipolarity has been attributed to such states. Information on the octupole states has been collected over many years, see Refs. [14] and [15] for reviews. Recently, measurements of a static octupole deformation in radium (Ra) and barium (Ba) isotopes have been reported [16, 17]. However, data on the octupole states are not so rich as in the case of quadrupole excitations. This is perhaps connected with experimental and interpretational difficulties. A common opinion seems to be that the quadrupole collective Hamiltonian stands for a pattern for higher multipolarities. However, it is not so easy. The theory of the octupole degrees of freedom is much more complicated than that of the quadrupole ones, or rather the quadrupole case is exceptionally simple. This will be demonstrated further in subsequent Sections. In Sec. 2 a general form of the octupole collective Hamiltonian is introduced and discussed in comparison with that of the quadrupole one. Possible definitions of intrinsic frames of reference for the octupole tensor in analogy to that for the quadrupole tensor are considered in Sec. 3.1. The octupole deformation parameters are introduced and illustrated. The octupole Hamiltonian is expressed in intrinsic coordinates in Sec. 3.2. Its approximate form for small oscillations around a pear shape is given in Sec. 3.3. In Conclusion, in Sec. 4, a draft of a unified quadrupoleoctupole collective model is recapitulated.
2 General form of collective Hamiltonians
The idea of collective models consists generally in describing some complex (collective) states of a manybody system through the substitution of the coordinates of many particles by a relatively small number of collective variables. An appropriate choice of these collective variables decides the success of the model. Constructing the collective models in question here for the description of lowlying collective states of eveneven nuclei one takes the spherical tensors as collective variables. These tensors are transformed according to irreducible dimensional representations D of the O(3) group of orthogonal transformations in the physical space. The spherical components () of in the laboratory frame, U, are the collective laboratory coordinates. It is assumed that the tensor is electric i.e. has parity and real, which means that its components fulfill the relation . The differential operators ^{a}^{a}aUnits are used here. play the role of the momenta canonically conjugate to the coordinates . The angular momentum operators or the O(3) generators in the collective space fulfill characteristic commutation relations with the coordinates and momenta (see e.g. Appendix A.1. in Ref. [18]). The following vector operator:
(0) 
fulfills such commutation relations and thus plays the role of the collective angular momentum [18, 19].
The nuclear collective system is defined by a collective Hamiltonian . It is assumed that has the following properties:

It is a secondorder differential operator (like the Schrödinger operator),

it is real, i.e. (it describes eveneven nuclei only),

it is a scalar with respect to the orthogonal group O(3) (rotations and inversion in the physical space) or commutes with the angular momentum operators of Eq. (2),

it is Hermitian with a scalar weight ,

it possesses the lowest eigenvalue (positive kinetic energy),

it is an isotropic function of the coordinates, i.e. it does not depend on any other tensor quantities.
No other assumptions are needed at this stage (cf. Ref. [20]). The most general form of with properties (1) – (6) given above reads [18, 12]:
(0) 
The Hamiltonian is determined by three quantities: two real scalar functions — weight and potential , and a symmetric matrix of real socalled inverse inertial functions . Since the matrix is transformed under the O(3) group like a product of two tensors it is called a symmetric bitensor. All these three quantities depend on no other tensors but and therefore they are isotropic tensor fields in the dimensional collective space. They can either be calculated from microscopic manybody models or fitted to experimental data. When the Hamiltonian is obtained by quantization of its classical counterpart the weight is .
In order to investigate possible structures of the inverse inertial bitensor it is convenient to express it as a set of (single) tensors (), namely
(0) 
Assumption (6), that the are isotropic functions of is essential here. Then, an arbitrary tensor can be expressed in the following form:
(0) 
by a number of definite fundamental tensors for given and , and arbitrary scalar coefficients (see Appendix A in Ref. [18]). The original Bohr inverse inertial bitensor [3] has the following form:
(0) 
where is a constant mass parameter.
In the case of a general quadrupole collective Hamiltonian it is well known that the inverse inertial bitensor can have at most six independent components out of fifteen possible ones [5, 18]. It is interesting to know whether there are similar restrictions for the components of the inverse inertial bitensor in the case of an octupole collective Hamiltonian. Unfortunately, the things are much more complicated in that case. The symmetric bitensor has twentyeight components. Can they all be independent? To answer this question one should analyze Eqs. (2) and (2) for .
According to Eq. (2) the even (positiveparity) symmetric octupole bitensor, , can be replaced with four tensors ( ). These tensors should be even isotropic functions of . The tensor fields in the space of the octupole coordinates are built out of some twentysix elementary tensors
(0) 
(square brackets stand for vector coupling to rank ) for and related to each other by about two hundred syzygies or relations in the form of rational integral functions [21]. All the elementary tensors split into two groups with positive and negative spinparity , respectively. Independent fundamental tensors from Eq. (2) are constructed by alignment of the elementary tensors of Eq. (2). The positive spinparity elementary and even (positive parity) fundamental tensors needed to construct the inverse inertial bitensor in question are listed in Sec. A.1. There are twentyeight relevant fundamental tensors altogether and thus, all twentyeight components of the bitensor can be arbitrary tensor fields in the octupole collective space. No additional relations between the components need appear. Twentyeight scalars for and are functions of the four elementary scalars, Eq. (A.1).
3 Intrinsic frame and intrinsic coordinates
3.1 Intrinsic frames of reference
Obviously, the tensor can be represented by different sets of coordinates in different frames of reference. As already stated in Sec. 2, the are the components of in the frame U. Let us take another frame, say U, the orientation of which with respect to the laboratory frame U is given by the Euler angles . We will not use the spherical components of in the frame U. Instead, we shall use their real and imaginary parts, and , respectively, defined in the standard way (cf e.g. Eq. A.7 in Ref. [18]). Then, the transformation rule between the corresponding components of the tensor with respect to U and U, respectively, takes the following form (cf. Ref. [19]):
(0) 
where
(0) 
are the semiCartesian Wigner functions. The BohrMottelson definition of the Wigner functions, , is used [22].
For the frame U to be the intrinsic frame, three properly chosen conditions for the coordinates and should be given, namely
(0) 
for , which determine the three Euler angles, as functions of the laboratory coordinates and, in this way, impart the status of intrinsic coordinates to them. The remaining independent coordinates are usually called deformations.
The concept of an intrinsic (bodyfixed) frame of reference is connected with the descriptions of collective states from the very beginning [3]. The principal axes of the tensor are taken as the intrinsic axes in the quadrupole, , case. According to Eq. (B.1) from Appendix B.2 the wellknown definition of the intrinsic frame is for . It is seen from Appendix 1 that the frame U is then Osymmetric. The two remaining intrinsic coordinates, and , are usually parametrized by the wellknown deformation parameters and .
The problem of the intrinsic frame for the octupole tensor, , has appeared to be less transparent than that for . The tensor has no principal axes and thus no obvious intrinsic frame. Early attempts to define an intrinsic frame failed (cf Ref. [23]). This was, perhaps, the reason why the intrinsic frame of reference was for a long time determined with the quadrupole tensor and octupole coordinates were treated as intrinsic coordinates with respect to that frame (cf. Ref. [19]). And what shall we do when any quadrupole tensor is not at our disposal?
In order to define an intrinsic frame for , having the natural symmetry O, it is convenient to decompose the tensor representation D of the O(3) group onto the O irreps A, F and F, respectively (see Appendix B.1). The corresponding representations will be called the octupole cubic Wigner functions , and , respectively. They are the following combinations of the semiCartesian Wigner functions:
(0) 
The cubic Wigner functions form a unitary set. When using them the transformation rule of Eq. (3.1) takes the following form:
(0) 
where , and are the octupole cubic coordinates in the frame U, proposed in Ref. [24] and defined in Appendix B.2.
It is seen that we have two possible definitions of the Osymmetric frame of reference. We can take Eq. ( ‣ 3.1) in the two following alternative forms: either or for . Here, we shall explore the former definition. Both of them are briefly discussed in Ref. [25]. We see from Eq. ( ‣ 3.1) that in the former case the following relation between the spherical laboratory coordinates and the Euler angles and the octupole deformations holds:
(0) 
The Jacobian of the transformation is equal to
(0)  
The transformation ( ‣ 3.1) is reversible for the deformations contained aside from the hypersurface , where ambiguities appear arising from the symmetries of the shape. The Fcovariant (vector) deformations () are transformed under the O transformations of the intrinsic frame like the Cartesian coordinates of a position vector.
It follows directly from the O symmetry that it is sufficient to consider, for instance, the region of vector deformations forming an infinite triangular pyramid in threedimensional space (see Fig. 1). The fortyeight pyramids obtained by all the O transformations fill up the entire space. Deformation supplements the space of vector deformations to the fourdimensional space. It is invariant under rotations and and changes sign under and inversion together with the corresponding transformations of the vector deformations. Therefore, there are no additional restrictions on values of (see Fig. 4 below). The case when is an exception: it is sufficient then to consider values .
To learn the geometrical meaning of the deformation parameters, and , one can investigate the shapes of the “octupoloid” given by the following equation in the spherical coordinates , , :
(0) 
in accordance with the approach of Ref. [26], bearing in mind that . The shape of the octupoloid is defined by the set of deformations , , and up to the cubic group of transformations i.e. the Bohr rotations and mirror reflections (see Appendix B.1).
Two octupoloids with identical axialsymmetric (pear) shapes but oriented differently are shown in Fig. 2. Obviously, the two sets of vector deformations then lie inside different pyramids of Fig. 1. Changing the sign of the deformation gives inverted octupoloids with the same shape as that in the figure. Examples of octupoloids with asymmetric shapes are shown in Fig. 3. Their orientations and/or handedness can be changed by cubic group transormations of the deformations. For instance, changing the signs of all the deformations gives inverted octupoloids.
However, when deformations belonging to one out of the two irreducible representations of the cubic group are transformed the shape of the octupoloid is changed, as shown in Fig. 4. Changing the sign of without changing the vector deformations gives a change of shape.
3.2 The Hamiltonian in intrinsic coordinates
A reversible transformation between the laboratory and intrinsic coordinates allows us to use interchangeably one or another set of variables. The use of intrinsic coordinates is usually more convenient because it gives the possibility to separate variables, especially the Euler angles. It is evident for potentials depending on the laboratory coordinates through the elementary scalars, the number of which coincides with the number of deformations. Hence, the potential in Eq. (2) is a function of the number of deformation parameters characteristic for a given multipolarity: two for and four for .
To express the kinetic part of the Hamiltonian and angular momenta in intrinsic coordinates we have to convert derivatives with respect to the laboratory coordinates into derivatives with respect to the intrinsic variables. As might be expected (cf. Ref. [27]), independently of the multipolarity of the collective space the components of the angular momentum of Eq. (2) can be expressed as
(0) 
where , , are the Cartesian components of the intrinsic angular momentum depending on the Euler angles and their derivatives only and are given by the following standard formulae: (cf. e.g. Eq. (2.15) in Ref. [18]).
(0)  
For the procedure for converting the derivatives is well known and is presented in detail, e.g., in Ref. [18]. A general quadrupole Hamiltonian expressed in intrinsic variables is divided into the vibrational part depending only on two deformations and the rotational Hamiltonian which contains the angular momenta , , and the deformation dependent moments of inertia. The intrinsic axes are always the principal axes of the tensor of inertia. When we take the kinetic energy with the Bohr inversed inertial bitensor (2) the structure of the Hamiltonian will not change much. Namely, the mixed term in the vibrational kinetic energy will vanish and the remaining five kineticenergy terms — vibrational and rotational — contain one common constant mass parameter instead of different deformationdependent inertial functions.
The procedure for transforming the octupole collective Hamiltonian to intrinsic coordinates is more involved than that for the case of . The first step of the procedure is conversion of the derivatives with respect to the laboratory coordinates into those with respect to the intrinsic variables. To do this one should calculate the derivatives of the intrinsic with respect to the laboratory coordinates. The transformation the reverse of that of Eq. ( ‣ 3.1) can be presented in the following entangled form:
(0) 
for . Derivatives of deformations and with respect to the laboratory variables are equal to
(0) 
Derivatives of the Euler angles with respect to are obtained by solving the following set of linear equations:
(0) 
To calculate all the derivatives from Eqs. (3.2) and(3.2) one should calculate the derivatives of the octupole cubic Wigner functions with respect to the Euler angles. Using Eqs. (3.1) and (3.1) the derivatives of the cubic functions can be expressed by derivatives of the Wigner functions themselves. E.g., handbook [28] contains formulae for the derivatives in question and other relevant properties of the Wigner functions^{b}^{b}bNote that the Wigner functions from Ref. [28] are the complex conjugate of those used here.. Finally, solutions of Eq. (3.2) for derivatives of the Euler angles with respect to are:
(0) 
where
(0) 
and are circular permutations of here and further below. Using Eqs. ( ‣ 3.2), (3.2) and ( ‣ 3.2) we are in a position to express derivatives with respect to the laboratory by derivatives with respect to the intrinsic coordinates in the following two equivalent ways:
(0) 
where the differential operators
(0) 
stand for angular momenta carried by vibrations of the octupole vector deformations and will be called the octupole vibrational angular momenta.
Using both versions of the righthand side of Eq. (3.2) and taking advantage of the unitarity of the octupole cubic Wigner functions we are able to express the Hamiltonian of Eq. (2) for in the intrinsic coordinates. To observe the inherent characteristics of the octupole collective motion an octupole Hamiltonian with the simplest inverse inertial bitensor, namely that of Eq. (2), will be presented here. This Hamiltonian when expressed in the Euler angles and octupole deformations is as follows:
(0) 
where the Cartesian tensor of the moments of inertia is equal to
It is seen that the intrinsic axes of an octupole system are not the principal axes of the moment of inertia as one would expect. The octupole vibrations, contrary to the quadrupole ones, carry their own angular momentum, which interacts by the Coriolis and centrifugal interactions with the total angular momentum (cf. Ref. [30]). This seems to be the most striking feature of the octupole rotations. In conclusion, the kinetic energy part of the Hamiltonian of Eq. (3.2) consists of ten terms: the four separate vibrations and the six rotational terms. In general, the Hamiltonian of Eq. (2) for , when expressed in the intrinsic variables, can contain additionally six mixed vibrational terms of type () and twelve vibrationrotation terms of type .
3.3 Axiallysymmetric deformation
Many eveneven nuclei show a static quadrupole deformation with axial symmetry. In the smalloscillations approximation of the collective Hamiltonian a simple picture of the quadrupole excitations has emerged (consult e.g. Chapter 6 in Ref. [29]). Two separate intrinsic vibrations appear, namely: the vibration of deformation around the equilibrium deformation , and the vibration () strongly coupled to rotations around the symmetry axis. On the other hand, rotations around axes perpendicular to the symmetry axis are weakly coupled to the vibrations and form characteristic rotational bands built on top of every vibrational level.
How is it in the case of a static axiallysymmetric octupole deformation? The equilibrium points in the fourdimensional deformation space are then supposed to be . The equilibrium shape is shown in Fig. 2(a). A rough approximation of the small oscillations around the two (because of the mirror symmetry) equilibrium points for the Hamiltonian of Eq. (3.2) (with the original Bohr kinetic energy) reads as follows:
(0) 
From the form of the Hamiltonian given above the following pattern of the small oscillations around the axialsymmetric octupole shape emerges, namely:

Two harmonic oscillations in coordinates and with stiffnesses and , respectively (x and yvibrations),

The doubleoscillator zvibrations around points with stiffness (cf. Ref. [31]),

The bvibration with stiffness strongly coupled to rotation around the symmetry axis z,

Rotations around the x and yaxes perpendicular to the symmetry axis with constant moment of inertia equal to .
The rotations around the x and y axes form rotational bands on top of vibrational levels. However, the bands are disturbed by the Coriolis interaction, being a kind of the rotationvibration interaction. In turn, centrifugal forces affect the four separate vibrations and can mix them with each other.
4 Conclusion
In the previous Sections a formalism for the octupole collective Hamiltonian has been presented and compared to that for the wellknown quadrupole one. For a few reasons, like a number of degrees of freedom greater by two, negative parity, additional simplifications in the quadrupole case, the theory of the octupole collective Hamiltonian is essentially more complicated, and therefore less developed than that of the quadrupole collective motion. A substantial feature of the octupole motion, which does not seem to be realized, is that the intrinsic vector x, y and zvibrations carry a nonzero angular momentum. This is obviously not the case for the quadrupole  and vibrations.
Obviously, a realistic collective model should take into account both modes, the quadrupole and the octupole together [19]. A separate consideration of the and cases either serves as a tool for developing a formalism and methods of treatment, or is an approximation. When we take, for instance, the kinetic energies of both modes with the Bohr inverse inertial bitensors of Eq. (2), the total quadrupoleoctupole Hamiltonian is the sum of the kinetic energies parametrized by the two mass parameters, and , respectively, and the potential , which can contain a possible quadrupoleoctupole interaction. However, modern collective Hamiltonians are extracted from microscopic theories, which seem to give inverse inertial bitensors dependent on both sets of coordinates for both ’s. Then, for instance, assumption no. (6) from Sec. 2 that the collective Hamiltonians contain isotropic functions of coordinates, is not valid. In consequence, the bitensor can have more than six independent components. Furthermore, it is natural to allow for the appearance of mixed quadrupoleoctupole terms in the total Hamiltonian. These terms would have the following form:
(0)  
Should be invariant under space inversion and Hermitian, is symmetric () and odd. By analogy to Eq. (2) the mixed bitensor can be presented in the following form:
(0) 
for . Tensors should have negative parity. The bitensor of Eq. ( ‣ 4) has 21 components altogether.
When the quadrupoleoctupole collective Hamiltonian is considered assumptions nos. (1) – (6) from Sec. 2 should be extended to the collective space of both tensors, and . The practical role of assumption no. (6) is that no material tensors appear for the nuclear medium. Under these extended assumptions the most general form of a quadrupoleoctupole Hamiltonian reads as follows:
(0)  
It is parametrized by 64 coordinatedependent inertial functions being components of the three inverse inertial bitensors, scalar weight and potential. The weight can possibly be equal to the square root of the determinant of the matrix of components of the inertial bitensors. The potential is a function of the coordinates through nine scalars described as deformations. In order to separate these nine variables from the twelve coordinates, a bodyfixed intrinsic frame of reference and intrinsic coordinates have to be introduced. One can do this in different ways. For instance, the principal axes of tensor oriented by the three Euler angles with respect to the laboratory axes can be treated as the intrinsic axes (cf. Ref. [19]). Then the two remaining intrinsic components of and all the seven intrinsic components of can be considered as deformations. Another way is to exchange the roles of tensors and and consider one of the frames defined in Sec. 3.1 through the octupole tensor as the intrinsic frame. One can also define two intrinsic frames for tensors and separately and treat the rotation of one frame with respect to the other as an intrinsic motion. Then the relative Euler angles have the status of deformations. Finally, in the case of a weak and well separated interaction between the modes one can treat them separately and then diagonalize the interaction within the product basis.
Only recently an attempt to solve a quadrupoleoctupole model, similar to that presented here, however not based on Hamiltonian ( ‣ 4) and with a restricted number of degrees of freedom, has been undertaken [33]. The model has been applied to the positive and negativeparity collective levels of the Gd nucleus. In any case, the problem of the full quadrupoleoctupole collective Hamiltonian is apparently complicated enough and still awaits practical applications to the spectroscopy of nuclear collective excitations.
Acknowledgments
The subject of the present paper was initiated some time ago by the first author (S.G.R.), then a Fellow of the Alexander von Humboldt Foundation, in collaboration with Walter Greiner, the then Director of the
Institut für Theoretische Physik der Johann Wolfgang GoetheUniversität.
The second author (L.P.) gratefully acknowledges support from the Polish National Science Center (NCN), grant no. 2013/10/M/ST2/00427.
Appendix A. Isotropic tensor fields in the octupole collective space
a.1 The positive spinparity elementary tensors
The twelve positive spinparity elementary tensors of Eq. (2) in the octupole collective space are as follows:
(A.1) 
a.2 Independent fundamental even tensors
Sets of fundamental even tensors with even ranks from 0 to 6 are listed in Table A.1. The choice of independent tensors need not be unique. This is because all fundamental tensors of a given rank (all possible alignments of the elementary tensors) are related to each other through a number of syzygies which can eliminate this or that tensor.