Moments of Inertia, Magnetic Dipole Moments, and Electric Quadrupole Moments of the Lithium Isotopes ()
1. Introduction
In the shell model, there is a core made up of paired nucleons. This core may be spherically symmetric in which case it gives rise to the spherical symmetric of the independent particle model or axially symmetric, as in the Nilsson model [1] which also be referred to as the deformed independent particle model. The analysis of nuclear spectra within the framework of the shell model is feasible only for relatively few nuclei: those which are close to magic numbers. The shell model has a limited range of useful applicability, covering nuclei in the vicinity of closed shells only. The complexity of nuclear spectra increases very rapidly as we go farther and farther away from closed shell [2] . The basic ideas concerning non spherical nuclei have been most completely described by A. Bohr and B. Mottelson [2] . A non-spherical nucleus is characterized by the moment of inertia about the axis perpendicular to the symmetry axis of the nucleus, its magnetic dipole moment, and its electric quadrupole moment. The elongation of the nucleus is related to the interaction between the surface and the nucleons outside closed shells.
A description of deformed nuclei has been given by a model proposed and developed by A. Bohr and B. Mottelson [2] . The success of the independent-particle approximation for spherical nuclei near closed shells naturally suggests adopting a similar procedure for deformed nuclei. Thus, as a first guess for the deformed nucleus internal wave function, one must take an independent-particle wave function, generated from a deformed potential. One of the most successful models for generating realistic intrinsic single particle states of deformed potentials is that first proposed by Nilsson [1] . This model was limited to nuclei with axially symmetric quadruple deformations. Positive values of the deformation parameter correspond to prolate deformation and negative values to oblate deformation. The success of the description of many nuclei by means of deformed potential can be taken as an indication that by distorting a spherical potential in this manner we automatically obtain the right combination of spherical eigenfunctions that makes the corresponding Slater determinant a better approximation to the real nuclear wave function [2] . From this point of view, the deformed potential is a definite prescription for a convenient mixing of various configurations of the spherical potential. Considerable evidence has accumulated for rotational structure for the p and s-d shell nuclei. The absolute values of the rotational energies or equivalently the moments of inertia require a knowledge of the fine details of the intrinsic nuclear structure. Consequently, the investigation of the nuclear moments of inertia is a sensitive check for the validity of the nuclear structure theories [3] - [11] .
The investigation of the different characteristics of the deformed nuclei, whether having axes of symmetry or not, is a very interesting subject in the theory of nuclear structure. Many authors have used several models to deal with this problem [12] - [19] . A common feature of systems that have rotational spectra is the existence of a deformation, by which is implied a feature of anisotropy that makes it possible to specify an orientation of the system. The moments of inertia require knowledge of the fine details of the intrinsic nuclear structure.
The large quadrupole moments observed in some nuclei, which do not belong to closed shells, implied a collective deformation and thereby a rotational degree of freedom. The most central parameter of collective rotation is the quadrupole moment and the moment of inertia of deformed nuclei [4] [13] [14] . The study of the velocity fields for the rotational motion of the axially symmetric deformed nuclei led to the formulation of the Schrödinger fluid [4] [5] [6] [7] [8] [14] [15] . Since the Schrödinger-fluid theory is an independent particle model, the cranking model approximation for the velocity fields and the moments of inertia play the dominant role in this theory.
The single particle Schrödinger fluid is one of the very interesting models which is created directly from the time-dependent Schrödinger wave equation by a suitably chosen type of complex wave functions [14] [15] . The single particle Schrödinger fluid is a concept which is used to describe the collective motions of the nucleons in an axially symmetric deformed nucleus. This concept can be applied to study the rotational motion of a deformed nucleus. This model makes it possible to formulate the well-known equation of continuity, Euler's equation, and Navier-Stokes equations of fluid mechanics [20] as results from the separation of the real and the imaginary parts of the time-dependent Schrödinger wave equation. The single-particle potential that represents the residual interparticle interaction inside the nucleus is taken in the form of a three-dimensional anisotropic oscillator.
Until now, the best description of the nuclear moments of inertia can be obtained within the framework of the so-called nuclear superfluidity model [16] [17] [18] [19] . The method adopted for the calculation of the eigenvalues and the eigenfunctions of the non-axial Hamiltonian makes it therefore possible to use the Belyaev formula for the determination of the nuclear moments of inertia.
It is well known that the nucleons inside the nucleus occupy approximately 1/50 of the volume of the nucleus. It is not surprising to find that nucleon properties are maintained inside the nucleus. In particular, this situation is responsible for the fact that the magnetic dipole moments of nucleons inside nuclei are the same as for free nucleons. In accordance with the above, we describe the motion of each nucleon individually in a common field. Because the surface of deformed nuclei is distorted at some moment, the potential felt by the nucleons is not spherically symmetric [18] [19] .
The next step was to include the effect of deformation on the single particle motion. Hartree-Fock calculations revealed that ground states of alpha-like nuclei in the p-shell are axially deformed [2] . The Nilsson model was limited to nuclei with axially symmetric quadrupole deformations, where the deformation is measured by the deformation parameter
.
In treating the internal motion in the nucleus, it is assumed that the individual nucleons move independently in a certain fixed non-spherical field of the nucleus. The Hamiltonian of the internal motion can then be represented, as in the ordinary model, in the form of a sum of one-particle Hamiltonians. According to Nilsson’s model [1] the nucleons inside the nucleus are moving independently in an averaging field in the form of anisotropic oscillator, with
, added to it a spin-orbit term and a term proportional to the square of the orbital angular momentum of the nucleon. The nucleon energy eigenvalues and eigenfunctions are then obtained by solving the time-independent Schrödinger wave equation in spherical polar coordinates and applying the method of diagonalizing the matrices [9] .
Because the surface is distorted at some moment, the potential felt by particles is not spherically symmetric, the particles will move in orbits appropriate to an aspherical shell-model potential. In the case of deformed nuclei, the theoretical question to be settled first is whether the nucleus has an axis of symmetry. Historically, several applications of the theory were made [5] [9] [11] on the assumption that the deformed nucleus does have such an axis of symmetry. Most of the work in the heavier nuclei is made on this assumption. In a field with axial symmetry, only the component of the angular momentum along the axis of symmetry is conserved. Historically, much after the application of the axially symmetric-rotor model, a systematic attempt was made by several authors, especially Davydov and his collaborators [12] [19] to check the consequences of the general rotor Hamiltonian that has no axis of symmetry (usually called an asymmetric rotor). The total motion of the nucleons is thus composed of two parts: an internal motion with respect to the body-fixed reference frame, described by an internal wave function, and the motion of the body-fixed reference frame itself. The nuclear superfluidity model [17] [18] [19] and the cranked Nilsson model [13] provide us with powerful methods for calculating the nuclear moment of inertia, the electric quadrupole moment, and other characteristics of deformed nuclei which have no axes of symmetry.
Accordingly, it is important to know whether the nucleus is deformed and has an axis of symmetry or not in order that we can choose the suitable model which produces correct values for its deformation characteristics, such as its moment of inertia, magnetic dipole moment, and electric quadrupole moment, which are in good agreement with the corresponding well-known experimental values.
In this paper, we carry out the derivations of the four models: the Nilsson model, the single-particle Schrödinger fluid model, the Cranked Nilsson model, and the nuclear superfluidity model, which are different in their structures and therefore in their applications, and accordingly clarify how the moments of inertia, the magnetic dipole moment, and the electric quadrupole moment of a deformed nucleus can be obtained in framework of these models. Our choice of the six lithium isotopes; 6Li, 7Li, 8Li, 9Li, 10Li, and 11Li, is due to the fact that their complete deformation structures are not yet completely known. Accordingly, we applied the above mentioned four models to the six lithium isotopes in order to calculate their moments of inertia, their magnetic dipole moments and their electric quadrupole moments. The variations of these moments with respect to the deformation parameter
, which describes the deviation from the spherical case, the non-axiality parameter
, which shows that the nucleus does not have an axis of symmetry, and the non-deformed oscillator parameter
, are also given in this paper.
2. The Nilsson Model
2.1. Formulation of the Model
In the Nilsson model [1] , all nucleons are assumed to move in the field of the following potential:
(2.1)
where:
(2.2)
(2.3)
In Equation (2.1), M is the nucleon mass,
is its orbital angular momentum vector, and
is its spin vector. Hence, the single particle Hamiltonian in the Nilsson model is given by:
(2.4)
where:
, (2.5)
In Equations (2.2) and (2.3)
is given by
(2.6)
where
is the deformation parameter. The parameter
is related to well-known deformation parameter
by the relation:
(2.7)
Transforming to the dimensionless variables
(2.8)
we get:
(2.9)
We rewrite Equation (2.9) as follows:
(2.10)
where:
, (2.11)
and
. (2.12)
Finally, we have
(2.13)
Schrödinger wave equation for the unperturbed Hamiltonian
is given by:
, (2.14)
with solutions [4] [10]
, (2.15)
(2.16)
where:
,
,
(2.17)
From the state
we construct the function
, as usual, in the form
. For Fermions, like nucleons,
. The nuclear state has definite values of the parity π, even or odd value of
, and Ω, where:
(2.18)
For the ground state,
, which is positive so that Ω is not a negative value. Hence, the least possible value of this quantum number is
. Finally, applying the variational method to the perturbed Hamiltonian
, with respect to the nuclear state
we obtain the final energy eigenvalues and eigenfunctions.
2.2. The Magnetic Dipole Moment
Using the total spin vector
and the total rotational vector
we can write the magnetic dipole moment as [1]
(2.19)
The complete method of calculating the magnetic dipole moment is given by Nilsson [1] .
2.3. The Electric Quadrupole Moment
Assuming a charge distribution in accordance with the Thomas-Fermi statistical model applied to the oscillator potential, one obtains, for the case of the axially symmetric nuclei, the intrinsic quadrupole moment, to the second order in the deformation parameter δ [1]
(2.20)
where Z is the number of protons and R is to be taken equal to the radius of charge of the nucleus. The relation between the measured quadrupole moment, denoted by
, and
is given by:
(2.21)
where I is the spin-quantum number of the specified nuclear state and K is its component along the body-fixed z-axis. Calculating the charge radius of the nucleus, the measured quadrupole moment for a nucleus with an axis of symmetry is then obtained as function of the deformation parameter δ.
3. The Single-Particle Schrödinger Fluid
3.1. Formulation of the Model
Let us consider a nucleus consisting of A nucleons. We assume that each nucleon in this nucleus (proton or neutron) has mass M and is moving in a single-particle potential
, which is deformed with time t, through its parametric dependence on a classical shape variable
. Here,
is assumed to be an externally prescribed function of t. Thus, the Hamiltonian for the present problem is given by [4] [7] [14] [15]
(3.1)
The single-particle time-dependent wave function
which satisfies the time-dependent Schrödinger wave equation that describes the motion of a nucleon, is defined as
(3.2)
To obtain a fluid dynamical description of the wave function
, we use the polar form of the wave function. We first isolate the explicit time dependence in the form:
(3.3)
where
is the energy density which depends on the time through the parameter
. Then, we write the complex wave function
in the following polar form:
(3.4)
where
and
are assumed to be real functions of
and
. Finally, we assume that the function
is positive definite. In the case of rotation, the parameter
becomes the angle of rotation,
, where Ω is the angular velocity. Substituting Equations (3.1), (3.3) and (3.4) into (3.2) and separating the real and the imaginary parts we get, when multiplying by 2Φ, a pair of coupled equations for Φ and S as follows:
. (3.5)
and
, (3.6)
We may call Equation (3.5) a modified Schrödinger equation because it differs from the usual time-independent Schrödinger equation
by an added term which we refer to as the “dynamical modification potential”:
(3.7)
Hence, we obtain two equations the first is:
(3.8)
where
is the irrotational velocity and
is the density. This equation is the well-known equation of continuity in fluid mechanics [20] . It can be rewritten in the form:
(3.9)
where
and
.
The second equation is:
(3.10)
which is a modified Schrödinger equation with:
(3.11)
Equation (3.4) can be written simply as
, so that:
. (3.12)
and
Therefore,
(3.13)
The current of the single particle state is defined by
[14] , so that:
, (3.14)
where
, and we finally get the current density:
(3.15)
Euler’s equation for the non-viscous fluid flow is given by [20]
(3.16)
where p is the pressure on the fluid at a point
at an instant of time t. For an ideal fluid,
is related to the enthalpy per unit mass, w, of the fluid by the following manner:
(3.17)
Therefore, Euler’s equation can be rewritten as
(3.18)
After integration and using
we get
(3.19)
so that
(3.20)
where S is the velocity potential for v, and we can write:
, (3.21)
So that the modified Schrödinger equation takes the form
(3.22)
where w is now the “enthalpy” of the single-particle Schrödinger fluid [14] [15] . Hence, we have a set of fluid dynamical equations completely analogous to those which describe a classical fluid.
For the modified potential we take the anisotropic oscillator, as proposed by Nilsson [1] with angular frequencies given by:
(3.23)
(3.24)
The parameter
depends on the deformation parameter
in the following way [1] [9]
, (3.25)
where
is the value of
for
. The deformation parameter
is related to the well-known deformation parameter
by (2.7).
The parameter
can vary in the range
.
Also, we assume that the adiabatic approximation is valid for this fluid, that is the angle of rotation
is constant of time. Hence, the collective kinetic energy T of the nucleus is given by [14] [15]
(3.26)
where
is the total density distribution of the nucleus and
is the total velocity field:
(3.27)
3.2. Moments of Inertia from Fluid Dynamical Viewpoint
We now examine the cranking moment of inertia in terms of the velocity fields. Bohr and Mottelson [2] showed that for harmonic oscillator case at the equilibrium deformation, where:
(3.28)
the cranking moment of inertia is identically equal to the rigid moment of inertia:
(3.29)
In terms of expression (3.19) involving the velocity fields, this result asserts the equality of the collective kinetic energy of the Schrödinger fluid and that of rigidly rotating classical fluid:
(3.30)
at the equilibrium deformation. We emphasize that the above equations hold for any number of nucleons occupying any set of single particle harmonic oscillator states at the deformation defined by equilibrium condition (3.28). It holds for a one particle state. For this case, Equation (3.30) becomes:
(3.31)
at the equilibrium deformation of the single particle state
(3.32)
Equation (3.31) is a remarkable identity. The scalar product of
and
which occurs on the left side is replaced on the right side, by the absolute square of
. It forces one to inquire whether the irrotational field
is equal to
. The answer, of course, is no. For,
posses compressible line vortices. It could be impossible to equal the velocity field for rigid rotation
, which has no singularity and is everywhere incompressible and rotational. Despite this qualitative difference between
and the other velocity in Equation (3.31), this shows that, as regards their effects under the integral upon the overall kinetic energy (or the internal parameter), these two velocity fields are equivalent at the equilibrium deformation. We note that the cranking moment of inertia
and the rigid moment of inertia
are equal only when the harmonic oscillator is at the equilibrium deformation. At other deformations, they can, and do, deviate substantially from one another [14] .
The following-expressions for the cranking moment of inertia,
, and the rigid-body moment of inertia,
, hold [14] [15] :
(3.33)
(3.34)
where E is the total nuclear ground-state energy defined by
(3.35)
and q is the ratio of the summed single particle quanta in the y-and z-directions:
(3.36)
The quantity q is known as the anisotropy of the configuration.
σ is the deformation of the potential [14] [15] . Concerning the magnetic dipole moment of an axially deformed nucleus, we can apply the same method as given by [21] .
4. The Cranked Nilsson Model
It is well known that nearly all fully microscopic theories of nuclear rotation are based on or related in some way to the cranking model, which was introduced by Inglis [18] [19] in a semi classical way, but it can be derived fully quantum mechanically, at least in the limit of large deformations, and not too strong K-admixtures (
). The cranking model has the following advantages [13] :
1) In principle, it provides a fully microscopic description of the rotating nucleus. There is no introduction of redundant variables, therefore, we can calculate the rotational inertia parameters microscopically within this model and get a deeper insight into the dynamics of rotational motion.
2) It describes the collective angular momentum as a sum of single-particle angular momentum. Therefore, collective rotation as well as single-particle rotation, and all transitions in between such as decoupling processes, are handled on the same footing.
3) It is correct for very large angular momenta, where classical arguments apply.
A simple and widely used way to describe the change of the single-particle structure with rotation is given by the cranked Nilsson model (CNM) [11] [13] . It is the method of calculating the shell correction energy that made it possible to do large-scale calculations where the nuclear potential-energy surface was explored in detail as a function of different deformation degrees of freedom. Important achievements in this field include the prediction of super deformed high-spin states and terminating bands.
4.1. The Single-Particle Hamiltonian
The single particle Hamiltonian in this model assumes the form [13]
, (4.1)
where:
. (4.2)
Here, the oscillator parameters
and
assume the form [13]
,
, (4.3)
,
where β (or ε) and γ are the quadrupole deformation degrees of freedom. The second term in the right-hand side of Equation (4.1) is given by:
, (4.4)
where the stretched square radius
is written in the form
, (4.5)
The hexadecapole potential is defined to obtain a smooth variation in the γ-plane so that the axial symmetry is not broken for
and
. It is of the form [13] :
, (4.6)
where the
parameters are chosen as
,
,
.
and
. (4.7)
In Equation (4.7) t refers to the stretched coordinates
, etc., and
in Equation (4.4) refers to the hexadecapole deformations degree of freedom.
4.2. Derivations
4.2.1. The Hamiltonian H(0)
The angular frequencies, Equations (4.3), can be simplified to
, (4.8)
, (4.9)
. (4.10)
Taking the squares of these equations and adding we get
(4.11)
Accordingly, the Hamiltonian
takes the form
(4.12)
The two deformation parameters
and
are equal and they are related to the well-known deformation parameter
by the following relation:
. (4.13)
Hence, the Hamiltonian
takes to the first order in
the form:
. (4.14)
4.2.2. The Hamiltonian H(1)
Direct substitution for the different quantities in the operator
gives:
(4.15)
4.2.3. The Term
Substituting for
and
in this term, we obtain
(4.16)
4.3. The Single Particle Energy Eigenvalues and Eigenfunctions
The method of finding the energy eigenvalues and eigenfunctions of the Hamiltonian H can be summarized as follows:
(i) Solving the Schrödinger’s equation
, (4.17)
exactly.
(ii) Modifying the functions
to become eigenfunctions for the solutions of the corresponding equation for
.
(iii) Using the functions obtained in step (ii) to construct the complete function
, the eigenfunction of the Hamiltonian H, in the form of linear combinations of the above functions, as basis functions, with given total angular momentum j and parity
.
(iv) Constructing the Hamiltonian matrix H by calculating its matrix elements with respect to the basis functions defined in step (iii).
(v) Diagonalizing the Hamiltonian matrix H to find the energy eigenvalues
and eigenfunctions
as functions of the non-deformed oscillator parameter
and the parameters of the potentials.
4.3.1. The Solutions of Equation (4.17)
The solutions of the equation
, are given, with the usual notations, by [8] [9] [10] [11]
, (4.18)
, (4.19)
where
are the normalized spherical harmonics with
and
is the nucleon orbital angular momentum quantum number.
The radial wave functions
are given by
, (4.20)
where
,
and the number of quanta of excitation N is related to the orbital angular momentum quantum number
by
or 1.
The last function in the right-hand side of Equation (4.20) is the associated Laguerre polynomial. Since the nucleon has spin
and intrinsic spin wave functions
, where
, the single particle wave functions of the Hamiltonian
are, then, given by
. (4.21)
4.3.2. The Eigenfunctions of the Hamiltonian
Wave functions with given values of the number of quanta of excitations N, the orbital angular momentum quantum number
, the total angular momentum quantum number J, and the parity π can be constructed from the functions (4.21), in the usual manner [8] [9] [10] [11] , as follows:
. (4.22)
The functions
are used as basis functions for the construction of the single particle nuclear wave functions with given total angular momentum J and parity π, in the usual manner, as follows:
. (4.23)
Accordingly, we obtain 15 wave functions, states, namely:
and
.
The matrix elements of the Hamiltonian
with respect to the functions (4.23) are given by:
(4.24)
4.4. The Hamiltonian Matrix Elements
4.4.1. The Matrix Elements of the Operator r2
The matrix elements of the operator r2 with respect to the basis functions
are given, with the usual notations [8] [9] [10] [11] , by
(4.25)
where
and
.
4.4.2. The Matrix Elements of the Operator
The matrix elements of the spherical harmonics
with respect to the basis functions
are given by [9]
. (4.26)
From the condition that
to be satisfied by the 3j-symbols in Equation (4.26),
, we get:
, (4.27)
, (4.28)
. (4.29)
4.5. Total Nuclear Quantities
We define the total quantities [13]
, (4.30)
, (4.31)
with the summation over the occupied orbitals in a specific configuration of the nucleus. The shell energy is now calculated from [13]
, (4.32)
where
is the smoothed single-particle sum evaluated according to the Strutinsky prescription. The detailed formulas for
are discussed for
and
in [13] .
The pairing energy is an important correction that should decrease with increasing spin and become essentially unimportant at very high spins. To obtain an (
) average pairing gap ∆, which varies as
, the pairing strength G is chosen as [13]
, (4.33)
with g1/g0 ≈ 1/3. Furthermore, the number of orbitals included in the pairing calculation should vary as
and
for protons (p) and neutrons (n), respectively.
The total nuclear energy is now calculated by replacing the smoothed single-particle sum by the rotating-liquid-drop energy and adding the pairing correction,
, (4.34)
or
(4.35)
where
,
is liquid drop energy,
and B is the smooth moment of inertia factor,
. The shell and pairing energies are evaluated separately for protons and neutrons at
, while the renormalization of the moment of inertia introduces a coupling when evaluating
for
. In the computer program,
is included only for
. The protons and neutrons are also coupled through the requirement that the shape of the respective potentials and the rotational frequencies are identical.
In the liquid drop model [3] [13] , the nuclear mass is given by:
(4.36)
(4.37)
In this formula,
and
,are the surface and Coulomb energies of a nucleus with a sharp surface in units of their corresponding values for spherical shape. The second term in the Coulomb energy is a (shape-independent) diffuseness correction with d being the diffuseness. The Coulomb energy constant is often defined as
. When calculating the nuclear mass, one should note that the average pairing energy should be subtracted from
.
The calculation of the Coulomb correction in particular, is somewhat involved: the original six-dimensional integral can be simplified only to four dimensions [13] . Furthermore, the use of stretched coordinates leads to complicated expressions. The radius expressed as a function of the angles in the stretched-coordinate system is obtained by requiring a constant value for the potential in (4.1) (neglecting the
term):
, (4.38)
where the spherical harmonics are functions of the angles
and
and the harmonic-oscillator part of the potential is expressed in
and
. From the definition of the stretched coordinates, it is straightforward to express the angles
and
in the corresponding angles in the spherical system,
and
.
Because of the incompressibility of nuclear matter, the nuclear volume is kept constant when the nucleus is deformed. This is achieved by varying the frequency
from its value for a spherical shape,
. The integration of the nuclear volume is most easily performed in the stretched-coordinate system and then multiplied with the corresponding Jacobian, a constant proportional to
. From the single-particle wave functions, the electric (or mass) quadrupole moment may be calculated as
, (4.39)
where
for protons and 0 (or 1) for neutrons. However, as we go further away from closed shells some new, very simple and systematic features start to show up for some nuclei. This is true for nuclei with mass number A in the range
, for A > 225, for nuclei in the s-d shell
and for p shell nuclei in
. Odd-A nuclei in these regions are characterized by exceptionally large positive quadruple moments, even-even nuclei in the same region all have a rather low-lying first excited state with
and electric quadruple radiation are strongly enhanced.
5. The Nuclear Superfluidity Model
5.1. Single Nucleon in a Deformed Nucleus
For a nucleus with quadruple deformation, one can write the nuclear radius as [2] [3]
(5.1)
where
is the radius of the sphere having the same volume and
are the spherical harmonics. Since the body-centered frame was selected as the principal axes, we have for the set of
in this frame
and
.
The Hamiltonian of single nucleon in this average field is given by [18] [19]
(5.2)
where:
(5.3)
(5.4)
The single-particle Hamiltonian then becomes:
(5.5)
The coefficients
can be expressed, for
, as
,
(5.6)
The frequencies
and
are connected with the deformation parameters
and
by [9]
, (5.7)
, (5.8)
. (5.9)
The single particle Hamiltonian then becomes:
(5.10)
This Hamiltonian does not produce the experimental single-particle energy levels so that we add to it two terms, suggested by Nilsson [1] , proportional to the spin-orbit term and the square of the orbital-angular momentum vector of the nucleon as follows:
(5.11)
The constants C and D are given by [1]
,
, (5.12)
The Hamiltonian (5.10) then becomes:
(5.13)
with
(5.14)
. (5.15)
(5.16)
The solutions of the Schrödinger equation corresponding to the Hamiltonian H0 are given, with the usual notations, by:
, (5.17)
. (5.18)
In Equation (5.17),
are the single-particle spin eigenfunctions with
,
are the normalized spherical harmonic functions with
, and the radial wave functions
are given by (4.20).
The Schrödinger wave equation
, (5.19)
can then be solved by applying the perturbation method [9] to obtain the perturbed energy eigenvalues
and the corresponding eigenfunctions
, to the second order of the approximation for
as follows:
, (5.20)
(5.21)
5.2. Calculations of the Moment of Inertia
5.2.1. Even-Even Nuclei
In the calculations, the chemical Potentials
and
are determined from the relation [9]
(5.22)
The expression for the moment of an even-even nucleus is given by [9] :
(5.23)
where
is given by the relation:
(5.24)
The only non-vanishing matrix elements
are:
(5.25)
5.2.1. Odd-Mass Nuclei
It is known that moments of inertia for the odd-mass nuclei are almost larger than that of neighboring doubly even nuclei. The odd-mass nucleus may be represented as a superfluid doubly even core plus the residual unpaired nucleons. In the present work, the moment of inertia
of the odd-mass nucleus is obtained as:
(5.26)
is the superfluid core contribution and
is the unpaired nucleon contribution. Of course, the contribution of even-even nuclei is not identical with that of the equivalent doubly even nucleus. This is clear since the unpaired particle will occupy one of the available levels [9] . Such partially filled levels correspond to unfilled ones with respect to the doubly even equivalent nucleus:
(5.27)
The unpaired nucleon contribution
was calculated with the help of the cranking model [9]
(5.28)
5.3. Calculations of the Electric Quadrupole Moment
For the non-axial case, the intrinsic quadrupole moment of a nucleus consisting of Z protons is given by [11]
(5.29)
where the single-particle operator
is given by
(5.30)
Carrying out the integration in Equation (5.30) with respect to the wave functions
which is evaluated in terms of the functions
, one then obtains
(5.31)
Filling the single-particle wave functions
for the given nucleus in its ground-state it is then possible to calculate the quadrupole moment by calculating the necessary matrix elements of Equation (5.31) and evaluating the expansion coefficients of the functions
in terms of the functions
as obtained from the variational and the perturbation methods.
6. Results and Discussions
According to previous works [9] [10] [11] , the Nilsson model parameters
and
are allowed to take on the values
and 0.08,
, for
and 2; and
for
,
takes values in the interval
with a step 0.01. Accordingly, the magnetic dipole moments and the electric quadrupole moments for the six p-shell nuclei 6Li, 7Li, 8Li, 9Li, 10Li and 11Li are calculated.
By assigning suitable values for the quantum numbers
and
, in accordance with the Nilsson model, we constructed the ground states of the six p-shell nuclei 6Li, 7Li, 8Li, 9Li, 10Li and 11Li by filling their ground states with successive single-particle states as given by Equation (3.32). For more details concerning this filing, see Appendix-1 in ref. [8] . Accordingly, the ground state in each nucleus is filled with the corresponding wave functions. As a first result,
and
are calculated for each nucleus. The corresponding reciprocal
moments of inertia
and
are calculated.
In Figures 1-6, we present the calculated values of the reciprocal moments of inertia according to the cranking model and the rigid-body model as functions of the deformation parameter
for the nuclei 6Li, 7Li, 8Li, 9Li, 10Li and 11Li, respectively.
Figure 1. Reciprocal moments of inertia of the nucleus 6Li.
Figure 2. Reciprocal moments of inertia of the nucleus 7Li.
Figure 3. Reciprocal moments of inertia of the nucleus 8Li.
Figure 4. Reciprocal moments of inertia of the nucleus 9Li.
Figure 5. Reciprocal moments of inertia of the nucleus 10Li.
Figure 6. Reciprocal moments of inertia of the nucleus 11Li.
In Table 1, we present the calculated values of the reciprocal moments of inertia of the nuclei 6Li, 7Li, 8Li, 9Li, 10Li and 11Li, which are in good agreement with the corresponding experimental values [22] . The experimental moments of inertia are also given in this table. The values of the nondeformed oscillator parameter
, for the six nuclei, are given. The values of the deformation parameter
are also given in this table.
It is seen from Table 1 that the calculated values of the cranking-model reciprocal moments of inertia are in better agreement with the corresponding experimental values rather than the other values. The disagreement between the values of the rigid-body reciprocal moments of inertia and the corresponding experimental values are due to the fact that the pairing correlation is not taken in concern in this model [9] [11] .
For the comparison between the calculated values of the magnetic dipole moments for the cases of axial symmetry and nonaxial (asymmetrical case), we only consider one nucleus as a sample case. For this purpose, we present in Table 2, the calculated values of the magnetic dipole moment of the nucleus ⁶Li, for both symmetric and asymmetric cases. Experimental values are also given.
For the electric quadrupole moments, we also did the same as for the magnetic
Table 1. The calculated values of the reciprocal moments of inertia of the nuclei 6Li, 7Li, 8Li, 9Li, 10Li and 11Li which are in good agreement with the corresponding experimental values are given. The values of
and
are also given in this table.
Table 2. The magnetic dipole moment of ⁶Li.
dipole moments. Accordingly, the calculated values of the electric quadrupole moment of 6Li (in e m barns) are given in Table 3 for the axially symmetric case (
) and for the non-axial case corresponding to
.
In Table 4, we present the calculated values of the reciprocal moment of inertia of the nucleus ⁶Li by using the concept of the single-particle Schrödinger fluid, for axially symmetric case, for both cranking-and rigid-body models, and the nuclear superfluidity model for the non-axial case. Also, we present in Table 4 the corresponding experimental value. The values of the deformation parameter
, the non-axiality parameter
and the oscillator parameter
are also given in Table 4.
In Table 5, we present the calculated values of the magnetic dipole moments of the six nuclei by applying the Nilsson model.
In Table 6, we present the calculated values of the electric quadrupole moments of the six nuclei by applying the cranked Nilsson model.
It is seen from Figures 1-6, and Table 1 that each nucleus has two values of the deformation parameter
which produce good agreement between the calculated and the corresponding experimental reciprocal moments of inertia for
Table 3. The electric quadrupole moment of ⁶Li.
Table 4. Reciprocal moment of inertia of ⁶Li.
Table 5. Calculated values of the magnetic dipole moments of the six nuclei by applying Nilsson model. The calculated values are also given [23] [24] .
Table 6. Calculated values of the electric quadrupole moments of the six nuclei by applying the cranked Nilsson model. The calculated values are also given [23] [24] .
the cranking model. Concerning the rigid-body model, the calculated values are not in good agreement with the corresponding experimental values, a result which always occurs with this model for most of the deformed nuclei [9] .
Finally, we see that the single particle Schrödinger fluid has been applied successfully to the six p-shell nuclei 6Li, 7Li, 8Li, 9Li, 10Li and 11Li with a suitable choice and filling of the single-particle anisotropic harmonic oscillator states especially for the cranking-model moments of inertia of deformed p-shell nuclei.
Concerning the calculation of the magnetic dipole moments of the six nuclei, it was found that the Nilsson model is an accurate framework for calculating these moments. The Cranking Nilsson model is the right model for the calculation of the electric quadrupole moments of the six nuclei. Furthermore, the nuclear superfluidity model also produced good results for the calculation of the nuclear moments of inertia for the nuclei which do not possess axes of symmetry.
According to the obtained results for the deformation characteristics, the assumption that the nucleus 6Li is deformed and has an axis of symmetry seems to be the more reliable assumption, a result which has been also obtained for some nuclei in the s-d shell [6] [11] .
7. Conclusion
To show that whether a nucleus is deformed and has an axis of symmetry, it is important to calculate its deformation properties, such as its moment of inertia, its magnetic dipole moment, and its electric quadrupole moment. These properties depend on the correct choice of the nuclear model which upon its application provides us with results in good agreement with the corresponding experimental values. For this reason, we choose in the present paper the well-known four models of deformed nuclei; two of them assume that the nucleus does have an axis of symmetry and the others assume that the nucleus does not possess an axis of symmetry. The six lithium isotopes 6Li, 7Li, 8Li, 9Li, 10Li and 11Li are not known to have axes of symmetry or not. Accordingly, it was our goal to choose these isotopes and apply the four well-known models of deformed nuclei: the single-particle Schrödinger fluid model, the Nilsson model, the cranked Nilsson model, and the nuclear superfluidity model to decide whether these nuclei have axes of symmetry or not. According to our calculations, the nucleus 6Li is the only nucleus which has an axis of symmetry.