The Triaxial Rotation Vibration Model in the Xe-Ba Region

11 Pages·0.12 MB·English

by  U. Meyer

The Triaxial Rotation Vibration Model in the Xe-Ba
Region
U. Meyer1, Amand Faessler, S.B. Khadkikar
Institute for Theoretical Physics, University of Tu¨bingen
Auf der Morgenstelle 14, D–72076 Tu¨bingen, Germany
8 Abstract
9
9
Collective quadrupoledegrees of freedomgive riseto vibrations
1
and rotations in nuclei. The axial Rotation Vibration Model
n
a (RVM) is here extended to describe also triaxial equilibrium
J shapes with β and γ vibrations allowing for the interaction

9
between vibrations and rotations. This Triaxial Rotation Vi-
1
bration Model (TRVM) is applied to Ba and Xe isotopes with

1 A 120 to 130. This area has recently been pointed out to be
v ≈
a region for the O(6) limit of the Interacting Boson Approx-
5
3 imation (IBA). The present work shows that the TRVM can
0 equally well describe these nuclei concerning their excitation
1
energies and E2 branching ratios.
0
8
9
/
h
t
-
l PACS numbers: 21.60 Ev, 23.20.-g, 23.20.Js, 27.60.+j
c
u
n
:
v Keywords: Collective nuclear models, triaxiallity, Rotation Vibration Model, Xe-Ba nu-
i
X clei, energy levels, E2 branching ratios.
r
a
1
Supported by the DFG Graduiertenkolleg MU 705/3
1 Introduction
The collective quadrupole degrees in nuclei yield rotational and vibrational excitations.
They are described by the Bohr-Mottelson Hamiltonian [1]. The wavefunctions for ro-
tations and β and γ vibrations in strongly deformed axially symmetric nuclei have been
greatly improved in the Rotation Vibration Model (RVM) by Faessler and Greiner [2, 3,
4, 5] taking the interaction between the rotations and vibrations into account with an
axially symmetric equilibrium deformation.
Recently, the Ba and Xe region with mass numbers A 120 to 130 has been studied
≈
experimentally [6, 9-15] and interpreted by several models: They compare different ap-
proaches of the Interacting Boson Approximation (IBA) in the U(5) and the O(6) limits
(with different E2 operators) [7] and the axial Rotation Vibration Model (RVM) and the
Asymmetric Rotor Model (ARM) [8]. They exclude the ARM [9, 11] because it cannot
describe a K=0 band built on γ vibrations [11] and because of the wrong staggering of
the 2+, 3+, 4+, 5+, ... excitation energies in the K=2 band. The argument with the K=0
band built on γ vibrations is trivial: The ARM does not contain γ vibrations and thus
cannot describe it. Because the staggering in the (K=2) quasi γ band is not described
correctly in the ARM, but given correctly in the RVM, the staggering seems to be due to
rotation vibration interaction. Although the RVM [2, 3, 4, 5] does quite well in describing
the data, the Cologne group concludes [11] that the O(6) IBA with the additional param-
eter χ in the E2 transition operator does in average agree better with the data. Here, we
want to show that this is connected with the restriction of the RVM [2, 3, 4, 5] to axial
symmetry while IBA allows also for triaxiallity (O(6) = γ instable limit). We present
here an extension of the axial RVM to triaxiallity (Triaxial Rotation Vibration Model =
TRVM). This model has the same number of parameters as the IBA used [9, 11]. We
obtain an equally good agreement with the data as in IBA.
InSection 2we present the modelwhich isanextension ofthe RVM[2,3,4,5]allowing
for a triaxial equilibrium deformation (TRVM). In Section 3 we compare the excitation
energies and the branching ratios in the best measured Xe and Ba isotopes. Section 4
summarizes the main results.
2 The Triaxial Rotation Vibration Model
In the Rotation Vibration Model [2, 3, 4, 5] one characterizes as in the Bohr-Mottelson
Model [1] the surface of the nucleus by quadrupole deformations:
R(θ,φ) = R (1+ α Y (θ,φ)) (1)
0 2µ 2µ
µ
X
The deformation parameters α are considered as dynamical variables and depend clas-
2µ
sically on time. Their behaviour is governed by the Hamiltonian:
1
H = B α˙† α˙ +V(α ) (2)
2 2µ 2µ 2µ
µ
X
The information about the equilibrium shape of the nucleus is contained in the potential
energy V(α ). The five quadrupole degrees of freedom α can be replaced by the three
2µ 2µ
Euler angles (φ,θ,ψ) for rotations and two deformation parameters
α = 2 (φ,θ,ψ)(β +a′(t))+( 2 + 2 )(a +a′(t)) . (3)
2µ Dµ0 0 0 Dµ2 Dµ−2 2 2
Here β and a give the equilibrium shape of the nucleus while a′(t) and a′(t) describe
0 2 0 2
vibrations around this shape. They are connected with the Bohr-Mottelson parameters
through:
β +a′(t) = βcos(γ)
0 0
1
a +a′(t) = βsin(γ) (4)
2 2
√2
In deriving the Hamiltonian for the TRVM we follow very closely the derivation for the
axial RVM in Ref.[4]. With the transformation (3) the Hamiltonian (2) can be written
(see in Ref.[4] eqs. (3) and (4)):
H = T +V = T +T +T +V (a′,a′) (5)
rot vib rotvib a0a2 0 2
The different terms are obtained by the straight forward transformation (3) and the
expansion in powers of a′/b and a′/a up to second order. This assumes that the β and
0 0 2 2
γ deformation is large compared to the vibrational amplitudes. This approach therefore
does not allow to describe spherical or nearly spherical nuclei. The several terms in the
TRVM Hamiltonian are given by:
ˆI2 Iˆ2 Iˆ2
3 3
T = − +
rot 2I 16Ba2
0 2
h¯2 ∂2 1 ∂2
T = ( + )
vib −2B ∂a′2 2∂a′2
0 2
2
ˆI2 Iˆ
T = − 3 f (β ,a ,a′,a′) (6)
rotvib 0 0 2 0 2
2I
0
Iˆ2 +Iˆ2
+ + −f (β ,a ,a′,a′)
1 0 2 0 2
2I
0
Iˆ2 a′
+ 3 f (a ,a′)+2ǫ 0
16Ba2 2 2 2 β
2 0
The functions f , f , and f are obtained by the expansion mentioned above:
0 1 2
a′ a′2 2
f = 2 0 +3 0 + (a2 +2a a′ +a′2)
0 − β β2 β2 2 2 2 2
0 0 0
1 1 1
f = √6 (a +a′) √6 a′(a +a′)
1 3 β 2 2 − β2 0 2 2
0 0
a′ a′2
2 2
f = 2 +3
2 − a a2
2 2
Forthe potentialenergy we assume a harmonicoscillator potential aroundtheequilibrium
shape:
1
V(a′,a′) = C a′2 +C a′2 (7)
0 2 0 0 2 2
2
For the diagonalization of the Hamiltonian we chose a basis of eigenstates of the Hamilto-
nian H = T +T +V. The eigenstates of the unperturbed Hamiltonian H are easily
0 rot vib 0
obtained and labelled by the total angular momentum I, the projection on the intrinsic
3-axis K and the number of phonons for the β-(n ) and γ vibrations (n ):
0 2
1
IK,n n = 2I +1 1 2 ( I +( )I I ) 1 βˆ† n0 0 1 βˆ† n2 0 (8)
| 2 0i (cid:18) 16π2 1+δK0(cid:19) DMK − DM−K sn0! (cid:16) 0(cid:17) | isn2! (cid:16) 2(cid:17) | i
Due to the rotation vibration part of the Hamiltonian T several eigenstates of H are
rotvib 0
mixed. The unperturbed energies as eigenvalues of H may easily be obtained. They are
0
given by:
h¯2 K2 1 1
EIK = (I(I +1) K2) + +(n + )E +(n + )E
n2n0 − 2I 16Ba2 0 2 β 2 2 γ
0 2
0,2,4,6,............ for K=0
with I =
K,K +1,K +2,... for K=0
(
6
C C
0 2
E = h¯ ; E = h¯ (9)
β γ
sB sB
2
h¯
ǫ =
I
0
2
I = 3Bβ
0 0
The operators βˆ† are creation operators for harmonic oscillations. We diagonalize the
Hamiltonian (5), (6), (7) in the complete basis (8). For numerical reasons the Hilbert
space has to be truncated. We chose the 31 lowest basis states with K 6 and up to two
≤
phonons (n +n 2). In notation (8) the most important states IK,n n are:
2 0 2 0
≤ | i
I0,00 ...... ground state band
| i
I2,00 ...... K=2 quasiγ band
| i
I0,10 ...... one-phononγ band (10)
| i
I0,01 ...... one-phononβ band
| i
I4,00 ...... K=4 quasiγ band
| i
The TRVM has as the IBA four independent parameters. The vibration energies E and
β
E , the inverse moment of inertia ǫ = 1/I and the triaxial equilibrium deformation a .
γ 0 2
The Hilbert space is too limited to describe the full variation of the moment of inertia.
As in the competing IBA [13] we make for the energies the Lipas ansatz
E = E /(1+α E ) (11)
0 L 0
E is the excitation energy obtained after diagonalization. The Lipas parameter α is
0 L
quite small ( 10−4[keV−1]). It describes a variable moment of inertia. We would like to
≈
note that the number of parameters including the Lipas parameter is the same in TRVM
and in IBA (including the effective boson charge in IBA). In the O(6) limit the branching
ratios of IBA depend only on the parameter χ. But similarly most of the branching ratios
in TRVM depend essentially only on the effective triaxiallity a + a′. In addition, one
2 2
could omit practically the β band and therefore the parameter E in TRVM, since the
β
results for the branching ratios do almost not depend on the β band and the agreement
would be equally good by omitting this band. But on the other side, this means also that
our prediction for the β bandhead is not reliable. Both models therefore effectively have
the same number of parameters.
3 Comparison with the Data
3.1 Energy spectra
In Figure 1 we display our results for the energy spectrum of 130Ba. The free parameters
which give an overall best fit of the experimental data from Ref.[13] are listed in Table 1.
For comparison, in Figure 1 furthermore the calculated IBA values and the experimental
results are shown. For the low spin states in the ground state band both theoretical
results are in nice agreement with experiment. As a result of the applied Lipas fitting
procedure lower spin states as the 4+, 6+ and also the 8+ state lie slightly higher while
the 10+-state is little below the experimental values. This is due to the stronger effect of
the Lipas procedure on higher spin states.
For the excited K=2 band (quasi γ band) we obtain in TRVM the experimentally
observed staggering of even-odd angular momentum states. A pure triaxial rotor model
without vibrations but with a γ deformation of the nuclear surface does not exhibit this
feature. The staggering thus seems to be the result of the rotation vibration coupling.
We would like to note that the staggering increases with excitation energies in TRVM.
In contrast, in IBA the staggering seems to be constant with growing excitation energy.
Both models agree qualitatively with the data in the γ band but not in detail.
The K=0 band is interpreted in the TRVM as the γ one-phonon band. The bandhead
(0+) at 1179 keV is fitted through the vibration energy E . For the TRVM the 0+–2+
γ
splitting is nicely reproduced and in better agreement with experiment than IBA. The β
vibrational bandhead is at around 1150 keV. Up to now, the levels of this band have not
been identified in the experimental analysis which might be difficult because the two 0+
energies are very close. Since our values for the experimentally observed energy levels are
almost not affected by the β band, we do not display it in Figure 1. Table 1
In Figure 2 we show our results for the 126Xe spectrum. There a different set of Figure1
parameters must be used in order to fit the experimental data best. This set is displayed
in Table 2. The agreement with the experimental level spacings is in the two theoretical
models about the same. This holds for the ground state band as well as for the quasi
γ band. In TRVM the staggering in low lying states is rather low. The free parameters
a /β and E have to be increased relative to these in 130Ba to get the best agreement.
2 0 β
Theexcited0+ bandisingoodagreement withexperiment andthelevel spacingsobtained
in TRVM are in better agreement with the experiment than the IBA calculations.
Table 2
Figure 2
3.2 E2 Branching Ratios
TheelectricquadrupoleoperatorcanagainbeobtainedfromtheoneinRef.[4]byreplacing
a′ a +a′. In collective intrinsic coordinates, it is given by:
2 2 2
→
1 1
3Z 2 5 2 4 5 2
m(E2,µ) = R2 2 β 1+ β + 2 a′ 1+ β +
4π 0Dµ0 0 7 π 0 Dµ0 0 7 π 0
(cid:18) (cid:19) (cid:18) (cid:19)
 1    
2 5 2
2 (a 2 2(a +a′)2)+ (12)
Dµ07 π 0′ − 2 2
(cid:18) (cid:19)
1 1
4 5 2 4 5 2
( 2 + 2 ) 1 β (a +a′) a′(a +a′)
Dµ2 Dµ−2  − 7 π 0 2 2 − 7 π 0 2 2 
(cid:18) (cid:19) (cid:18) (cid:19)
  
The reduced E2 transition probabilities are defined by [2, 3, 4, 5]:
2I +1
f 2
B(E2;I I ) = I m(E2) I (13)
i f i f
→ 2I +1|h || || i|
i
We have calculated the measured B(E2) ratios to be able to compare them with the
experimental results for both 130Ba and 126Xe. In Table 3 the results for 130Ba are given.
We see a fair agreement between the results of TRVM and the experiment as well as those
+
of IBA except for the transition from the initial 4 state. Here our result is about two
3
orders of magnitude larger while the IBA result is two orders of magnitude lower than the
+ + + +
experimental ratio B(E2,4 2 )/ B(E2,4 4 ). A similar discrepancy exists also
3 2 3 2
+ + → →
for the transition (4 2 ) as can be seen from the last entry in Table 3. Any B(M1)
3 1
+→ +
contribution to the 4 4 transition cannot cure this disagreement. Table 3
3 2
→
Table 4 shows the ratios of certain B(E2) transitions in 126Xe . Here the agreement
is much better between the models and the experiment; and both Interacting Boson
Approximation and the Triaxial Rotation Vibration Model agree nicely with the data.
This might be because only pure quadrupole transitions are used in Xe to define the
ratios, while the uncertainty in the E2/M1 ratio in the normalising states is ignored in
the case of Ba and transitions are assumed to be pure E2 transitions.
Table 4
4 Conclusions
The Triaxial Rotation Vibration Model (TRVM) is successfully applied to Ba and Xe
nuclei to obtain low lying energy levels. These results compare well with those obtained
in the O(6)-limit of IBA and with the experimental values both for the spectra and for
relative E2 transition rates. The bandhead of the β vibrational band are not seen in
experiment so far. Unfortunately, we are not able to specify very exactly the energy
region for the bandhead or other levels of the β band since our calculated values are not
very sensitive to the parameter E .
β
Acknowledgment: We would like to thank Dr. Ingo Wiedenh¨over and Prof. Peter von
Brentano for providing us the IBA results and for fruitful discussions.
References
[1] A. Bohr, B. Mottelson, Dan. Mat. Fys. Medd. 27 (1953) No.16
[2] A. Faessler, W. Greiner, Z. Phys. 168 (1962) 425
[3] A. Faessler, W. Greiner, Z. Phys. 170 (1962) 105 and 177 (1964) 190
[4] A. Faessler, W. Greiner, R.K. Sheline, Nucl. Phys. 70 (1965) 33
[5] A. Faessler, W. Greiner, R.K. Sheline, Phys. Rev. 135 (1964) B591 and Nucl. Phys.
62 (1965) 241
[6] R.F. Casten, P. von Brentano, Phys. Lett. B152 (1985) 22
[7] F. Iachello, A. Arima, ”The Interacting Boson Model”, Cambridge Univ. Press, 1987
[8] A.S. Dawydow, V.S. Rostovsky, Nucl. Phys. 12 (1959) 58, A.S. Dawydow, Nucl.
Phys. 24 (1961) 682 and A.S. Dawydow, A.A. Chaban Nucl. Phys. 20 (1960) 499
[9] W. Lieberz, A. Dewald, W. Frank, A. Gelberg, W. Krips, D. Lieberz, R. Wirowski
and P. von Brentano, Phys. Lett. B240 (1990) 38
[10] A. Dewald, P. Sala, R. Wrzal, G. B¨ohm, D. Lieberz, G. Siems, W. Wirowski, K.O.
Zell, A. Gelberg, P. von Brentano, P. Nolan, A.J. Kirwan, P.J. Bishop, R. Julin, A.
Lampinen, J. Hattula, Nucl. Phys. A545 (1992) 822
[11] F.Seiffert, W.Lieberz, A.Dewald, S.Freund, A.Gelberg, A.Granderath,D.Lieberz,
R. Wirowski and P. von Brentano, Nucl. Phys. A554 (1993) 287
[12] G. Siems, U. Neuneyer, I. Wiedenh¨over, S. Albers, M. Eschenauer, R. Wirowski, A.
Gelberg, P. von Brentano, Phys. Lett. B320 (1994) 1
[13] K. Kirch, G. Siems, M. Eschenauer, A. Gelberg, R. Ku¨hn, A. Mertens, U. Neuneyer,
O. Vogel, I. Wiedenh¨over, P. von Brentano, T. Otsuka, Nucl. Phys. A587 (1995) 211
[14] R. Ku¨hn, K. Kirch, I. Wiedenh¨over, O. Vogel, M. Wilhelm. U. Neuneyer, M. Luig,
A. Gelberg, P. von Brentano, Nucl. Phys. A597 (1996) 85
[15] A. Dewald, D. Weil, R. Kru¨cken, R. Ku¨hn, R. Peusquens, H. Tiesler, O. Vogel, K.O.
Zell, P. von Brentano, D. Bazzacco, C. Rossi-Alvarez, P. Pavan, D. DeAcun˜a, G. De
Angelis, M. De Poli, ”A Dynamic Shape Effect in 126Ba at Low-Spin”, preprint 1996
130
Ba
3500.0
+ +
10 TRVM exp.[13] 10 IBA
+
10 +
3000.0 + + 8
8 8
+
2500.0 8+ 7+ 8+ 7+ 8+ 7
+
V] + + 6
e 6 6
k
gy [ 2000.0 + + +
n Ener 6+ 5+ 2+ 6+ 54+ 2+ 6+ 54+ 2
atio 1500.0 4+ +
Excit + 3+ + + 3+ 0 + 3+ 0+
1000.0 4 0 4 4
+
+
2 +
2 2
+ + +
500.0
2 2 2
0.0
Fig.1: Ground state, quasi γ and K=0 band for 130Ba in keV for the
IBA calculation and the Triaxial Rotation Vibration Model (TRVM) compared
with experimental values from Ref.[13].
126
Xe
3500.0
+
10 TRVM + exp.[13] + IBA
+ 10 + 8
10
3000.0 8 +
8
+
7
+ + + +
2500.0 8 7 8+ 7 8
+ +
V] + 6 6
e 4 +
k + 4
gy [ 2000.0 6 5+ 5+
Ener 6+ 5+ 2+ 6+ + 2+ 6+ + 2+
n 4 4
citatio 1500.0 4+ + + 0+ 3+ +
Ex + 3+ 0 + 3 + 0
1000.0 4 4 4
2+ + +
2 2
+ + +
500.0
2 2 2
0.0
Fig.2: Ground state, quasi γ and K=0 band for 126Xe in keV for the
IBA calculation and the Triaxial Rotation Vibration Model (TRVM) compared
with experimental values from Ref.[13].
E [keV] E [keV] ǫ[keV] a /β α [keV 1]
β γ 2 0 L −
1168 1179 80.5 0.329 10 4
−
Table 1: Parameters of the TRVM in the present calcula-
tions for 130Ba. E and E are the vibration energies, ǫ is
β γ
the inverse moment of inertia and a2/β0 the relative ratio of
the static γ and β deformations. α is the Lipas parameter.
L
E [keV] E [keV] ǫ[keV] a /β α [keV 1]
β γ 2 0 L −
1769 1314 81.3 0.333 10 4
−
Table 2: Parameters of the TRVM in the present calcula-
tions for 126Xe. E and E are the vibration energies, ǫ is
β γ
the inverse moment of inertia and a2/β0 the relative ratio of
the static γ and β deformations. α is the Lipas parameter.
L
I I IBA–CQF exp.[13] TRVM
i f
→
+ +
2 0 6.2 6.2 0.7 4.5
2 1
→ + ±
2 100.0 100.0 100.0
1
+ +
3 2 100.0 100.0 100.0
1 2
→ +
4 29 22.0 3.0 16
1
+ ±
2 6.1 4.5 0.6 2.6
1
+ + ±
4 2 100.0 100.0 100.0
2 2
→ +
4 58 54.0 10.0 57
1
+ ±
2 0.12 2.3 0.4 2.75
1
+ + ±
0 2 100.0 100.0 100.0
2 2
→ +
2 2.5 3.3 0.2 3.0
1
+ + ±
2 0 100.0 100.0 100.0
3 2
→ +
2 0.11 21.0 4.0 0.04
2
+ ±
4 1.3 2.7 0.5 0.11
1
+ ±
2 0.0 3.3 0.6 0.06
1
+ ±
0 0.24 0.017 0.003 0.04
1
+ ±
3 112 18
1
+ +
4 2 0.02 2.9(5) 670
3 2
→ +
3 108 97(17) 214
1
4+ 100.0 100.0∗ 100.0
2
4+ 0.01 3.4(6)∗ 7.2
1
+
2 0.0 0.30(6) 23.3
1
Table 3: B(E2) branching ratios of the TRVM
for 130Ba in comparison with IBA and with ex-
periment from Ref.[13].