# Two-Center Gaussian Potential Well for Studying Light Nucleus in Cluster Structure

Authors: Nafiseh Roshanbakht; Mohammad Reza Shojaei
Journal: Advances in High Energy Physics
 (2017)

---

## Abstract

The clustering phenomena are very important to determine structure of light nuclei and deformation of spherical shape is inevitable. Hence, we calculated the energy levels of two-center Gaussian potential well including spin-orbit coupling by solving the Schrödinger equation in the cylindrical coordinates. This model can predict the spin and parity of the light nuclei that have two identical cluster structures.


---

## Body

## 1. Introduction

In the light nucleus, deformation plays an important role in determining nuclear structure. These deviations from spherical structures are found in the axial deformations and the clustering because numerous experimental studies have revealed clustering phenomena in them [1, 2]. Freer and Merchant in 1997 studied the role of clustering and cluster models in nuclear reactions and examined the evidence for α-cluster chain configurations in the even-even light nuclei from Be8 to Si28 [3, 4]. In 2001, Kanada-En’yo and Horiuchi showed that the valence neutrons in nucleus with cluster structure can form a covalent bond between the clusters and named it “nuclear molecules” [5]. For example, Be9 has a neutron that exchanges between the two alpha clusters. They studied nuclear molecules by the one-center deformed harmonic oscillator potential [5, 6].For studying light nucleus in cluster model, many physicists used various methods such as antisymmetrized molecular dynamics (AMD) [7] and fermionic molecular dynamics (FMD) [8] and condensates and the THSR Wave-Function [9–11]. These models have many important advantages for studying nuclei in cluster models but do not make assumptions about the cluster [12].The two-center shell model is more appropriate for modeling nuclei with clustered nature. This model assumes a clustered structure, and each cluster is represented by its own potential. This potential is the two-center harmonic oscillator (TCHO) [13]. It has a complete analytical solution but does not have a sufficiently sharp edge to satisfy finite separation energies.In this article, we consider axially symmetric nuclear parameterization, so the nuclear surface equation is given in the cylindrical coordinates [14, 15]. In this case, the deformation removes the degeneracy of single nucleon levels associated with spherical symmetry [16]. In the next section, we introduce new potential for studying the nucleus with the cluster structure.


## 2. 2D Potential in Cylindrical Coordinates

If the nucleus structure consists of two identical clusters, the nuclear surface in the cylindrical coordinates is defined by [15, 17](1)ρ=b1-z+c1a2z≤z1ρ2-R2-z-c22z1≺z≺z2b1-z-c1a2z2≤z,where ±c1 are the positions of the two cluster centers and other parameters are shown in Figure 1.Figure 1
Parameterizations of nuclear shape in two cluster structures.Hence, we use two-center Gaussian potential well for calculating the energy levels of nucleus which are nearly spherical or deformed because of having cluster structure. Our potential in cylindrical coordinates can be split into two parts:(2)V0ρ,z=Aexp⁡-αz-z02+Bexp⁡-βρ2,z≥0Cexp⁡-λz+z02+Bexp⁡-βρ2,z<0,where A, C, and B are the potential depth, the parameters α, β, and λ are inverse square range, and z0 is cluster centers. We study nucleus that has two identical cluster structures. Hence, we suppose that the parameters A=C and α=λ have a symmetric potential form. In Figure 2, the potential is plotted for some values of z0. The nuclei with the closed shell such as He4 nucleus have a spherical shape and are exceptionally stable. Hence, the spherical Gaussian potential well (Figure 2(a)) can describe energy levels of the nucleus that has a spherical shape. But in light nucleus with almost half full shells N=2, the nucleons have been observed to cluster structure with a two-center structure. For example, Be8 isotope at the ground state has two alpha cluster structures and decays to alpha particles by 92 keV energy. In this case, the potential in Figure 2(d) can be used to describe it. Likewise, for light neutron-rich nucleus that have covalent binding between two clusters due to valence neutrons, the potentials are shown in Figures 2(b) and 2(c).Figure 2
The Gaussian potential for some values ofz0.
(a)
z 0 = 0 fm (b)
z 0 = 3 fm (c)
z 0 = 4 fm (d)
z 0 = 5 fmIf the potential of the quantum system to be examined is axially symmetric, then the Schrödinger equation in the cylindrical coordinates can be employed. For a three-dimensional problem, the Laplacian in cylindrical coordinates is used to express the Schrödinger equation in the following form:(3)-ħ22M∇2ψρ,ϕ,z+V0ρ,zψρ,ϕ,z=Eψρ,ϕ,z,where E and M denote the energy eigenvalue and the mass, respectively. V0(ρ,z) is potential that does not depend on ϕ because of the axial symmetry and is separable to V0(ρ) and V0(z). In the cylindrical coordinates, the Laplace operator takes the following form:(4)∇2=∂2ψ∂ρ2+1ρ∂ψ∂ρ+1ρ2∂2ψ∂ϕ2+∂2ψ∂z2.The separation of variables is accomplished by substituting(5)ψρ,ϕ,z=RρZzΦϕ.In the usual way, this leads to the following ordinary differential equations [18]:(6)d2Φϕdϕ2+m2Φϕ=0,(7)d2Zzdz2-2Mħ2V0z-EZz+γ2Zz=0,(8)d2Rρdρ2+1ρdRρdρ-γ2+m2ρ2Rρ-2Mħ2V0ρRρ=0.The solution of the first equation is(9)Φϕ=12πexp⁡imϕ,where m is the magnetic quantum number and integer. The next step is to obtain the solutions of (7) and (8) which cannot be solved by analytical methods. But the Hamiltonian is diagonal, and there are some numerical methods to solve it. We use a well-known Numerov’s numerical method to solve ordinary differential equations of second order [19, 20]. In this method, we first organize the square matrix Hamiltonian of the system consisting of the kinetic energy matrix and the potential energy matrix; then we obtain the eigenvectors of Hamiltonian and the eigenvalues for the stationary states of the time-independent Schrödinger equation. The values of M and ħ are considered as M=1 and ħ=1. The eigenvectors Rnρ(ρ) for three lowest states are shown in Figure 3 and the eigenvectors Znν(z) are plotted for some values of z0 in Figure 4. nρ is the radial quantum number and nz is the quantum number along z-axis having values 0,1,2,….Figure 3
The eigenvectorsRnρ(ρ) of three lowest states.Figure 4
The eigenvectorsZ(z) of three lowest states for some values of z0.
(a)
z 0 = 0 fm (b)
z 0 = 3 fm (c)
z 0 = 4 fm (d)
z 0 = 5 fmThe eigenvectorsRnρ(ρ) vanish on the boundaries of our system and behave as expected. For the eigenvectors Znν(z), in addition to boundary conditions, we consider the continuity conditions too:(10)
Z
n
ν
z
⟶
0
=
Z
n
ν
z
⟶
0
z
<
0
z
>
0
,
(11)
∂
Z
n
ν
z
⟶
0
∂
z
=
∂
Z
n
ν
z
⟶
0
∂
z
z
<
0
z
>
0
.While we describe nuclei in an axially symmetric model, the orbital angular momentum,l, and the intrinsic spin, s, are not good quantum numbers because states with different l-values with the same parity can mix. In this model, the quantum number N=nρ+nz is the principal quantum number and the parity is determined by -1N. The magnetic quantum number has the integer values between −nρ and nρ too.Figures4(a)–4(d) show that the wave functions associated with z-axis have a definite parity with respect to z=0. For finite separations, the wave functions do not feel the presence of the other cluster across the barrier in low energy. The energy levels with the same quantum number which located in the first or in the second potential well have the same value.The eigenvalues from the solution of the Hamiltonian equation inz-axis are obtained for the different values of the distance between the two cluster centers, and the variation of them is given in Figure 5. In this figure, the shell structure of the spherical nucleus is reproduced for the two separate clusters, in higher energy. At z0=0, there are the principal levels without spin-orbit coupling (N,m) and, for large values of z0, the levels start to form two clusters and their levels overlap with each other.Figure 5
Variations of the eigenvalues inz-axis with the distance between the two cluster centers.One of the interesting features of Figure5 is that states having the same z-component of angular momentum, m, repulse each other and prevent crossing. This is consistent with the theory of Neumann and Wigner. They showed that such level crossings are highly improbable [21].The levels with different values ofnρ and nz but the same value of N are degenerate. For example, (nz,nρ,m)=(1,0,0) are degenerate with (nz,nρ,m)=(0,1,0). We calculated the spin-orbit coupling to eliminate this degeneracy.


## 3. The Effect of Spin-Orbit Coupling on Levels

After the presentation of the solutions of the Schrödinger equation with the potential in (2), we want to explain the effect of spin-orbit coupling on the single particle levels for spherical nuclei and for nuclei with two cluster structures separated at an infinite distance.The spin-orbit coupling is defined in the form of [15, 22](12)VL·S=-λ12m0c2L→·S→,where λ=35 is a dimensionless coupling constant; m0 and c are the nucleon and the speed of the light, respectively [15]. The orbital angular momentum operator is L→=∇→V×p→ and we must transform it into cylindrical coordinates:(13)L^±=∓ħe±iφ∂V∂ρ∂∂z-∂V∂z∂∂ρ±i∂V∂z1ρ∂∂φ;L^z=iħ∂V∂ρ1ρ∂∂φ,hence,(14)L→·S→=12L+s-+L-s++Lzsz.The shift of the energy levels due to the spin-orbit coupling was calculated using the orthogonality of the wave functions. The matrix elements of the spin operators are(15)ms′s^+ms=ħδms′,ms+1,ms′s^-ms=ħδms′,ms-1,ms′s^zms=ħmsδms′,ms,and the spin-orbit matrix is (16)nz′,nρ′,m′,ms′12L+s-+L-s++Lzsznz,nρ,m,ms=ħ2nz′,nρ′,m′L+nz,nρ,mδms′,ms-1+ħ2nz′,nρ′,m′L-nz,nρ,mδms′,ms+1+ħmsnz′,nρ′,m′Lznz,nρ,mδms′,ms.Now, we use the parabolic approximation [23], and consequently the levels of single particle energies are obtained as a function of z0 by the diagonalization of the potential V0+VL·S. The structure of levels is plotted in Figure 6 and labeled them by jz. In our model, the projection of the total angular momentum on the symmetry axis, jz=m+ms, is the good quantum number. Hence, states of different jz are not coupled by the Hamiltonian but the states with ±jz are degenerate because the nuclei have reflection symmetry for either of the two possible directions of the symmetry axis.Figure 6
Variations of the eigenvalues (nz,nρ,m) in z-axis including spin-orbit coupling with the distance between the two cluster centers.In Figure6, the level scheme of the spherical nuclei is obtained for z0=0. But, by increasing z0, two separated clusters are formed. An example of light nuclei by two cluster structures is Be9 isotope. The nucleon single particle level scheme is plotted in Figure 6. Be9 isotope is stable and has 4 protons and 5 neutrons. Alpha cluster is in the first level and in the second level there will be other alpha clusters and the remaining single nucleon located in the next level. Experimental evidence for Be9 shows that it has a ground state with spin and parity 3/2− and a 1/2+ excited state [24, 25]. For predicting these results, we must use the potential well in Figure 4(b). When two separate centers in potential well are formed, the levels (1,0,0) and (0,1,0) overlap with each other. In our model, the single neutron is in the level Jz=3/2 at the ground state and its parity is -1N=-1. For the first excited state, the single neutron moves into level Jz=1/2+. We will predict other excited states by considering the rotational band and Coulomb potential between the clusters.


## 4. Conclusion

In light nuclei, close to decay threshold, the cluster phenomenon is the favored mode. Hence, we use two-center Gaussian potential well for calculating the energy levels of nucleus which are nearly spherical or deformed because of having two identical cluster structures. We solved the Schrödinger equation in the cylindrical coordinates and plotted the level scheme without considering spin-orbit coupling. Then we modified results with it. Forz0=0, the level scheme of the spherical nuclei is obtained, and two separated clusters are formed by increasing z0. In principle, a two-center potential well can be used to study light nuclei with two cluster structures and make assumptions about the cluster in them.



---

*Source: 9309636-2017-02-12.xml*