Configuration Mixing in Pt Nuclei within Interacting Boson Model-2 ()
1. Introduction
Some nuclei near closed shells appear to have both the vibrational structure expected for a near-spherical shape, and rotational structure, which is typical of deformed nuclei [1]. This phenomenon of shape coexistence involves two configurations of the nucleus which have different numbers of active nucleons. In an IBM description, the two configurations have different boson numbers; Nv being the same but Nπ different, or vice versa. The most common situation involves a difference in Nπ (Nv) of two bosons between the normal configuration and the so-called intruder configuration, corresponding to a pair excitation across a shell or sub-shell gap [1] [2].
Shape coexistence in atomic nuclei has become a very active field of research during the last decades and clear signals of its existence have been obtained at and near proton or neutron closed shells [1] [2] [3], more in particular in light nuclei with a closed neutron shell at N = 8, 20, 28, and 50 closed shells as well as in heavy nuclei such as the Sn and the Pb nuclei. It seems that, without exception, shape coexistence is associated with the presence of low-lying excited 0+ states.
In the case of the IBM-2, shape coexistence arises including two-particle two-hole (2p-2h) (or even higher np-nh). Excitations across the closed shells, but they consider them as extra bosons, i.e., pairs of nucleons. This extension is called IBM-2 configuration mixing. In two previous papers [4] [5], we used the IBM-CM to extensively study the Pt nuclei. We carried out a detailed analysis of the energy spectra and absolute B (E2) values for states up to an energy of 1.5 MeV. Garc’ıa-Ramos et al., [5], studied shape evolution and shape coexistence in Pt isotopes: Comparing interacting boson model configuration mixing and Gogny mean-field energy surfaces. In this work, the evolution of the total energy surface and the nuclear shape in the isotopic chain 172–194Pt are studied in the framework of the interacting boson model, including configuration mixing. The results are compared with a self-consistent Hartree-Fock-Bogoliubov calculation using the Gogny-D1S interaction and a good agreement between both approaches is found.
In this article, it will be studied to extract the parameters describing the IBM-2 CM Hamiltonian precisely. After showing some theories and laws about IBM-2, it represents the results of the Pt isotopes comparing with some experimental data on energy levels, electromagnetic transition probabilities (B(E2) and B(M1)), monopole transitions, mixing ratios and quadrupole moments.
2. The Interacting Boson Model-2 (IBM-2)
2.1. The Hamiltonian of Interacting Boson Model-2 (IBM-2)
As the IBM-2 model, the protons and neutrons degrees of freedom are taken into account explicitly. Consequently, one can write the Hamiltonian [6] [7] as,
(1)
(2)
Here
is the d-boson energy,
is the strength of the quadrupole interacting between proton and neutron bosons.
The quadrupole moment operator, in the IBM-2 model, has the form:
(3)
where
or
,
is the quadrupole deformation parameter for protons
and neutrons
. Where the terms
and
are the proton-proton and neutron-neutron d-boson interacting only and written by:
(4)
The Majorana interaction is represented the term
, that accounts for the symmetry energy and shits the states with mixed proton-neutron symmetry state with respect to the fully symmetric ones which affects only the relative location of the states with mixed symmetry, it is with respect to the totally symmetric states. For the reason that little experimental information is familiar states with mixed symmetry that has the form:
(5)
2.2. Electromagnetic Transitions and Quadrupole Moments in IBM-2
Generally, the E2 transition operator of one-body in the IBM-2 is
(6)
(7)
which
is in the form of Equation (3). For simplicity, the
has the same value as in the Hamiltonian. Also, one suggests it by the single j-shell microscopy. Generally, the E2 transition results are not sensitive to the choice of
and
, whether
or not. Thus, the reduced electric quadrupole transition rates between
states are given by:
(8)
In IBM-2, the electric quadrupole moment is written as:
(9)
One can calculate IBM-2 eigen functions and energy eigenvalues are usually achieved numerically by the program code NPBOS [8]. Then, the result of eigenvectors can be calculated transition rates and related properties using the program code NPBTRN [8]. The relationship is between the parameters of Equation (2).
2.3. Configuration Mixing in Interacting Boson Model-2 (IBM-2 CM)
Configuration mixing can be treated in the IBM-2 using a technique developed by Duval and Barrett [9]. Separate IBM-2 calculations are achieved for the two configurations and the results are then mixed this is done by the interaction configuration mixing can be remedied in the IBM-2 by a technique method developed by Duval and Barrett [9]. Separate IBM-2 calculations are done for the two configurations and the results are then mixed using the interaction
(10)
where the intruder configuration is assumed to involve the proton shell. There are three parameters in the mixing calculation, the mixing strengths
and
in Equation (23), and the pair excitation energy,
, which gives the relative energies of the two unperturbed configurations.
The total mixing Hamiltonian is then given by
(11)
where H1 (H2) is the IBM-2 Hamiltonian for the first (second) configuration, as given by Equation (11), and an amount fl has been added to the energies of the second configuration.
The mixed wave functions are used to calculate B (E2) values of observed transitions and quadrupole moments. The E2 transition operator is given by [9]:
(12)
where
were defined in Equation (3) and
and
are adjustable parameters. The suffixes 2 and 4 refer to the normal and the intruder configurations respectively.
For simplicity,
and
.
For a mathematical simplicity, the neutron boson and proton boson effective charges are often taken to be equal, and the parameters
and
(in intruder configuration) are taken to be the same as the Hamiltonian parameters
and
respectively. The T(E2) transition operator can then be written as:
(13)
When two configurations are present, this operator becomes [9]:
(14)
where
or
and i denotes the configuration. The effective charges of the two configurations are not the same, in general.
3. Results and Discussion
3.1. Choice of Parameters
The normal configuration for platinum isotopes involves
(sometimes denoted as 2π, two proton boson holes), counting from the Z = 82 closed shell. The neutron configuration for
for example, is
(four neutron boson holes), counting from the N = 126 closed shell. The vibrational spectra can be calculated by diagonalizing the IBM-2 Hamiltonian, (Equation (2)), in the space of two proton and
neutron s and d bosons. In order to describe the rotational states, an alternative configuration must be specified and a separate set of IBM-2 calculations made, based on that configuration. The alternate configuration used for the 186-200Pt isotopes involves a two-particle-four-hole excitation in the shell model proton space [1] [9] [10]. This corresponds to two proton boson particles and two proton boson hole in the IBM-2 space. For simplicity, the proton boson particles and hole are treated equivalently, even though the underlying fermion pair degrees of freedom originate in different major shells.
The IBM-2 calculations have been done in model spaces with
and
to describe the vibrational and rotational states, respectively, the two calculations are combined using Equation (11).
In the phenomenological calculations the parameters appearing in the Hamiltonian ((Equation (2) and Equation (11)) in two configurations may in general depend both on proton (
) and neutron (
) boson number. Guided by the microscopic theory as we have assumed that only
and
depend on
and
i.e.,
,
, while
depends only on
and
on
, i.e.,
,
. Thus a set of isotopes, (constant
) have the same value of
, while a set of isotones, (constant
), have the same value of
. This parameterization allows one to correlate a large number of experimental data. Similarly, when a proton-proton
and neutron-neutron
, interaction is added, the coefficients
are taken as
and
, i.e. the proton boson interaction will only depend on
, and the neutron boson on
. Since there is no information on the location of the states with mixed neutron-proton symmetry we kept the coefficients appearing in the Majorana force
(see Table 1 and Table 2).
The values of the parameters used for the present calculations are given in Table 1. The value of the parameter (boson energy)
for the
configuration,
, is constant for all isotopes and the values of
for the
configuration,
, are nearly constant. The quadrupole-quadrupole interaction strengths
trends for (
) both configurations follow the microscopic predictions [9].
The values of the parameter
used for Pt isotopes are the microscopic predictions from Bijker et al., 1980 [11]. They were reported only for neutron number
. The values of
for Pt isotopes were determined by extrapolating the microscopic trend to larger neutron number. This was done in
Table 1. IBM-2 Hamiltonian parameters for the Pt isotopes (normal configuration
). All energies are in MeV;
and
are dimensionless.
,
,
.
Table 2. IBM-2 Hamiltonian parameters for the Pt isotopes (normal configuration
). (2P-2h), All energies are in MeV;
and
are dimensionless.
,
,
.
a manner consistent with the phenomenologically determined
values which were used for the neighboring platinum and osmium isotones [11]. The same values of
are used for both configurations.
The emphasis in this work is on describing overall trends with constant or smoothly varying parameters of the Hamiltonian (Equations. (2), (11)), rather than obtaining the best possible fit to the experimental data for each nucleus. This is done in an effort to find a set of IBM-2Hamiltonian parameters which is appropriate for the entire isotopic chain.
The normal configuration for platinum involves
(sometimes (
) denoted as (
), two proton boson hole), counting from the Z = 82 closed shell. The neutron configuration for 196Pt for example, is
(four neutron boson holes), counting from the N = 126 closed shell. The vibrational spectra can be calculated by diagonalizing the IBM-2 Hamiltonian, Equation (2), in the space of one proton and
neutron s and d bosons. In order to describe the rotational states, an alternative configuration must be specified and a separate set of IBM-2 calculations made, based on that configuration. The alternate configuration used for the Pt isotopes involves a four-particle-four-hole excitation in the shell model proton space. This corresponds to two proton boson particles and two proton boson hole in the IBM-2 space. The two configurations are depicted schematically energy levels in Table 2. For simplicity, the proton boson particles and hole are treated equivalently, even though the underlying fermion pair degrees of freedom originate in different major shells.
3.2. Energy Spectra
The configuration mixing calculations are done using the computer code NPMIX [10], which calculates the energy eigenvalues and eigenfunctions. The computer codes NPBEMX and BEMIX [10] are subsequently used to calculate matrix elements for transition rates and other properties.
The mixing between the two configurations (
) and (
) is apparent in the experimental data shown in Figures 1-3. Consider, for example, the
,
and
,
states in 186-196Pt isotopes. In the latter, the vibrational (2) state is lower in energy than the rotational (
)
and
states. The states are close in energy and mutually repel, with the vibrational 4+ state being lower in energy and the rotational 4+ state higher in energy than would be expected without configuration mixing.
In the second configuration (
) configuration interaction parameters
,
and
for all isotopes. For most nuclei, the spectral structure comes primarily from the pairing and quadrupole-quadrupole terms in the IBM-2 Hamiltonian, Equation (2), and the
term is relatively unimportant. Since there is only two proton boson to interact through
, the residual interaction among the neutron bosons
interaction is comparatively large.
From Figure 1, in the normal configuration the ground state band the moment of inertia increases, beta band (
) is pushed up, where the state
becomes bandhead (a member of γ-band) of
. In 194-200Pt isotopes, only five
states are observed around 1.25 MeV is produced, some of these states do not have a
symmetry.
In the IBM-2-CM approach, the lightest Pt isotopes, are deformed (172-178Pt). In 186Pt isotope a prolate shape and a γ-soft minimum coexist, but a well-deformed prolate minimum quickly develops in 180Pt isotope, becoming the most pronounced prolate minimum at the mid-shell, i.e., in 182Pt isotope with the prolate shape remaining well pronounced up to 186Pt isotope. Moving towards heavier mass Pt isotopes, γ-flat energy surfaces start to develop. For 188Pt isotope, a much extended energy surface develops in the γ-direction (deviation from symmetry axis), becoming completely γ-unstable when reaching 190–196Pt isotopes, the 198-200Pt isotopes tends to spherical shape because approach to magic number.
We show in Figure 1 the evolution of the energy of the ground state band (
,
,
, and
) increases toward the middle of the major shell with the number of the valence neutrons and remain almost constant for
isotopes. Although these tendencies are well reproduced, the rotational properties or rotational features are somewhat enhanced in the calculated levels for
, which are slightly agreement in energy with the experimental data. From both results, the IBM-2 CM and the experimental data, we observed the fingerprints for structural evolution with a jump between 186Pt and 188Pt isotopes, also the shape of the nucleus changed from prolate (186Pt isotope) to oblate (188Pt isotope) deformation, this case is called phase transition. The 188Pt isotope takes the oblate shape because the yrast states gradually departure to neutron closed shell.
One can also find signatures for a shape or phase transition in the systematics of the quasi-β-band levels shown in Figure 2. In 180-186Pt isotopes configuration mixing, the
baldhead and the state
nearly constant in both IBM-2 data and experimental data. The two levels are pushed up rather significantly from
Figure 1. Ground state band in Pt isotopes. Experimental data are taken from Refs. [12] - [21].
Figure 2. Beta state band in Pt isotopes. Experimental data are taken from Refs. [12] - [21].
186Pt isotope to 188Pt isotope. Consistently with the systematics in the ground-state band. The calculated
and
states are higher than, but still follow, the experimental data. The
state lies close to the
state. However, this calculation predicts that this trend persists even for 188-198Pt isotopes, whereas the energy spacing between the experimental
and
states for these isotopes is larger. Similar deviation occurs for high spin states. This means that our calculations suggest the feature characteristic of the O (6) symmetry (γ-unstable).
The γ-band states shown in Figure 3, we can see the good agreement between IBM-2 calculations and experimental data for
, where the
,
and
states are energetically closed.
The high spin states for example.,
,
,
, states, we observe some deviation occurs for these states. This means that our calculations suggest the feature characteristic of the O (6) symmetry, where the staggering occurs as high spin beta band states such as
,
,
, ..., etc., However, the experimental
Figure 3. Gamma state band in Pt isotopes. Experimental data are taken from Refs. [12] - [21].
levels are lying more regularly, particularly for
, and thus appear to be in between the O(6) symmetry and triaxial rotor symmetry. The deviation of the γ -band structure seems to be nothing more than a consequence of an algebraic nature of the IBM-2, and indeed has also been found in existing phenomenologically calculations of IBM-2.
A characteristic feature of the γ-unstable limit of the IBM-2 is a bunching of γ-band states according to 2+, (3+, 4+), (5+, 6+), ..., that is, 3+ and 4+ are close in energy, etc. This even-even staggering is observed in certain O(6) nuclei but not in all and in some it is, in fact, replaced by the opposite bunching (2+, 3+), (4+, 5+), ... which is typical of a rigid triaxial rotor [22]. From these qualitative observations it is clear that the even–even γ-band staggering is governed by the γ-degree of freedom (i.e., triaxiality) as it changes character in the transition from a γ-soft vibrator to a rigid triaxial rotor. The energy surface becomes rather flat, evolving towards a spherical minimum at 200Pt and beyond. The possibility of triaxial deformation was not considered (tends to spherical shape).
3.3. Electric Transition Probability B(E2)
Calculations of electromagnetic properties give us a good test of the nuclear models prediction. The electromagnetic matrix elements between eigenstates were calculated using the programs NPBTRN and NPBEMX for IBM-2 and IBM-2 CM model respectively.
From Equations (7), we note that an E2 transition mainly depends on identifying proton and neutron bosons effective charges
and
. These isotopes lying in region between SU(3) limit (rotational nuclei) and O(6) limit (γ-soft nuclei), therefore, the relationship between (
,
) and the reduced transition probability B(E2) for rotational limit SU(3) and γ-soft limit O(6) is given in the form [23] [24]:
For O(6) limit
(16)
For SU(3) limit
(17)
where
is the experimental reduced transition probability from the first excited states (
) to the ground state (
), N is the total number of bosons.
The relations (26) and (27), was used to estimate the effective boson charges for proton and neutron bosons (
,
). In these calculations, we use the following criteria to determine the effective charges.
is a constant throughout the whole isotopic chain and the
changes with neutron number. This is true if the neutron (proton) interaction does not depend on the proton (neutron) configurations. The values of
and
are determined by fitting to the five
and
in 194Pt and 196Pt isotopes from the first configuration (
). They are given in Table 3.
For the configuration mixing, the effective charges for bosons is evaluated in the same manner in normal configuration, For simplicity,
and
. The ratio of the two quadrupole effective charges,
in Equation (25) is taken to be the same quadrupole interaction strengths, as the ratio of the corresponding
, for each isotope (see Equation (12) and Table 1, Table 2). This is reasonable, since the effective charge and the strength of the quadrupole interaction are both proportional to the mean square proton radius. Thus, the only new parameter needed to determine the reduced transition rates is
. The value for the Pt isotopes was determined by fitting of
in 196Pt isotope.
It is well known that absolute gamma ray transition probabilities offer the possibility of a very sensitive test of nuclear models and the majority of the information on the nature of the ground state has come from studies of the energy level spacing. The transition probability values of the excited state in the ground state band constitute another source of nuclear information. Yrast levels of even-even nuclei (
) usually decay by E2 transition to the lower lying yrast level with
.
In Table 4 we show the
and
values, which are of the same order of magnitude and display a typical decrease towards the middle of the shell.
As a consequence of possible M1 admixture the
quantity is rather difficult to measure. There is no experimental data to compare the values
Table 3. Effective charge used in E2 transition calculations (
).
Table 4. Electric transition probability
in e2b2 units.
of IBM-2. For 190Pt isotope, we give the agreement with experimental value, from these values seems to decrease for 190-196Pt isotopes.
In Table 4, we show
values. Experimentally the results are radically different for the Ptisotopes. In some Ptisotopes the value seems to decreased towards the middle of the shell, whereas in another 196Pt isotope is increased. Our calculations could not reproduce these contradictory features simultaneously. The results for
values are shown in Table 4. This quantity is rather small since this transition is forbidden in all three symmetries of the IBM-2.
As a consequence of possible M1 admixture the
quantity is rather difficult to measure. For some Ptisotopes, we give the different, conflicting experimental results and we see that no general feature be derived from them, from these values seems to decrease for Ptisotopes.
In 194Pt isotope, there is a good agreement between experimental and calculated
value. This confirms our earlier statement about the nature of the lowest
state in this isotope. Other transitions such as
,
and
are small values in sometimes because these transitions between different bands (cross over transitions) and the selection rules which determine these transition.
In Table 5, the quadrupol moment, qualitatively, for first excited state
, second excited states
and
, for the ground state band, the positive
and
mean a negative intrinsic quadrupole moment (for ground state
). For the beta band, a negative
means a negative intrinsic quadrupole moment
. The negative
implies that the nucleus has an oblate shape; the positive intrinsic quadrupole moment
means that the nucleus has a prolate shape. The overall the IBM-2 and IBM-2 CM results is a good agreement with the experimental data.
4. Concluding Remarks
In the present work, we have analyzed the spectroscopic properties for the Pt isotopic chain 182-200Pt in terms of the interacting boson model-2 (IBM-2) and interacting boson model configuration mixing (IBM-2 CM). Calculations of energy levels for even-even 182-200Pt isotopes were performed with the whole
Table 5. Quadrupole moments
in eb Units for Pt isotopes.
Hamiltonian (Equation (2) and Equation (11)) using NPBOS computer code. 182-200Pt nuclei (Z = 78) have two proton bosons (
normal configuration) and the neutron boson varied from
for 182Pt isotope to
for 200Pt isotope.
The energy levels in Pt isotopes especially (ground state band) are a good agreement with the experimental data. In this study, we see that the excitation energy of the state
in 186-200Pt isotopes (heavier Pt isotopes) is predicted sometimes too high in comparison with experimental data.
In the IBM-2-CM approach, the lightest Pt isotopes are deformed (172-178Pt). In 186Pt isotope a prolate shape and a γ-soft minimum coexist, but a well-deformed prolate minimum quickly develops in 180Pt isotope, becoming the most pronounced prolate minimum at the mid-shell, i.e., in 182Pt isotope with the prolate shape remaining well pronounced up to 186Pt isotope. Moving towards heavier mass Pt isotopes, γ-flat energy surfaces start to develop. The 198-200Pt isotopes tend to spherical shape because of the approach to magic number.
Concerning electric transition probability B(E2) values, we find that in all calculations the overall trend is reproduced reasonably very well in some transitions, but notice some discrepancies present case of the decay of the
states in heavier Pt isotopes, in general, better than the values calculated by Bjjker et al., [11].