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

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.

Nafiseh Roshanbakht 1 and Mohammad Reza Shojaei 1

Academic Editor:Akpan Ikot

Department of Physics, Shahrood University of Technology, P.O. Box 36155-316, Shahrood, Iran

Received 24 October 2016; Revised 29 December 2016; Accepted 16 January 2017; 12 February 2017

This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The publication of this article was funded by SCOAP³ .

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] [figure omitted; refer to PDF] 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.

[figure omitted; refer to PDF]

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: [figure omitted; refer to PDF] 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 of _z0 .

(a) _{z 0} = 0 fm

[figure omitted; refer to PDF]

(b) _{z 0} = 3 fm

[figure omitted; refer to PDF]

(c) _{z 0} = 4 fm

[figure omitted; refer to PDF]

(d) _{z 0} = 5 fm

[figure omitted; refer to PDF]

If 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: [figure omitted; refer to PDF] where E and M denote the energy eigenvalue and the mass, respectively. _V0 (ρ,z) is potential that does not depend on [varphi] because of the axial symmetry and is separable to _V0 (ρ) and _V0 (z). In the cylindrical coordinates, the Laplace operator takes the following form: [figure omitted; refer to PDF]

The separation of variables is accomplished by substituting [figure omitted; refer to PDF]

In the usual way, this leads to the following ordinary differential equations [18]: [figure omitted; refer to PDF]

The solution of the first equation is [figure omitted; refer to PDF] 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 [Planck's over 2pi] are considered as M=1 and [Planck's over 2pi]=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 eigenvectors _Rnρ (ρ) of three lowest states.

[figure omitted; refer to PDF]

Figure 4: The eigenvectors Z(z) of three lowest states for some values of _z0 .

(a) _{z 0} = 0 fm

[figure omitted; refer to PDF]

(b) _{z 0} = 3 fm

[figure omitted; refer to PDF]

(c) _{z 0} = 4 fm

[figure omitted; refer to PDF]

(d) _{z 0} = 5 fm

[figure omitted; refer to PDF]

The eigenvectors _Rnρ (ρ) 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: [figure omitted; refer to PDF]

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 ^(-1)N . The magnetic quantum number has the integer values between -_nρ and _nρ too.

Figures 4(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 in z-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 in z-axis with the distance between the two cluster centers.

[figure omitted; refer to PDF]

One of the interesting features of Figure 5 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 of _nρ 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] [figure omitted; refer to PDF] 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[arrow right]=∇[arrow right]V×p[arrow right] and we must transform it into cylindrical coordinates: [figure omitted; refer to PDF] hence, [figure omitted; refer to PDF] 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 [figure omitted; refer to PDF] and the spin-orbit matrix is [figure omitted; refer to PDF] 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.

[figure omitted; refer to PDF]

In Figure 6, 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 ^(-1)N =-1. For the first excited state, the single neutron moves into level _Jz =1/²⁺ . 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. For _z0 =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.

References

[1] W. von Oertzen, M. Freer, Y. Kanada-En'yo, "Nuclear clusters and nuclear molecules,", Physics Reports , vol. 432, no. 2, pp. 43-113, 2006.

[2] Y. Kanada-En'yo, T. Suhara, F. Kobayashi, "Cluster formation, breaking, and excitation in light nuclei," in Proceedings of the International Nuclear Physics Conference (INPC '13), Florence, Italy, June 2013.

[3] M. Freer, A. C. Merchant, "Developments in the study of nuclear clustering in light even-even nuclei,", Journal of Physics G: Nuclear and Particle Physics , vol. 23, no. 3, pp. 261-323, 1997.

[4] M. R. Shojaei, N. R. Bakht, "Calculation of energy spectrum of 12C isotope with modified Yukawa potential by cluster models,", Pramana , vol. 87, no. 4, article 54, 2016.

[5] Y. Kanada-En'yo, H. Horiuchi, "Structure of light unstable nuclei studied with antisymmetrized molecular dynamics,", Progress of Theoretical Physics Supplement , no. 142, pp. 205-263, 2001.

[6] M. Freer, "Clusters in nuclei,", Scholarpedia , vol. 5, no. 6, article 9652, 2010.

[7] D. M. Brink, "History of cluster structure in nuclei,", Journal of Physics: Conference Series , vol. 111, no. 1, 2008.

[8] H. Feldmeier, "Fermionic molecular dynamics,", Nuclear Physics A , vol. 515, no. 1, pp. 147-172, 1990.

[9] A. Tohsaki, H. Horiuchi, P. Schuck, G. Röpke, "Alpha cluster condensation in ¹² C and ¹⁶ O,", Physical Review Letters , vol. 87, no. 19, 2001.

[10] Y. Funaki, A. Tohsaki, H. Horiuchi, P. Schuck, G. Röpke, "Analysis of previous microscopic calculations for the second 0+ state in 12C in terms of 3-α particle bose-condensed state,", Physical Review C , vol. 67, no. 5, 2003.

[11] T. Yamada, P. Schuck, "Dilute multi-α cluster states in nuclei,", Physical Review C--Nuclear Physics , vol. 69, no. 2, pp. 24309, 2004.

[12] M. Freer, "The clustered nucleus--cluster structures in stable and unstable nuclei,", Reports on Progress in Physics , vol. 70, no. 12, pp. 2149-2210, 2007.

[13] C. L. Zhang, B. Schuetrumpf, W. Nazarewicz, "Nucleon localization in light and heavy nuclei," https://arxiv.org/abs/1607.00422

[14] M. Mirea, "Realistic orbital momentum operators for the superasymmetric two-center shell model,", Nuclear Physics A , vol. 780, no. 1-2, pp. 13-33, 2006.

[15] M. Mirea, "Two center shell model with Woods-Saxon potentials,", Romanian Reports in Physics , vol. 59, no. 2, pp. 523-531, 2007.

[16] J.-P. Ebran, E. Khan, T. Niksic, D. Vretenar, "How atomic nuclei cluster,", Nature , vol. 487, no. 7407, pp. 341-344, 2012.

[17] D. Aranghel, A. Sandulescu, "Shell effects in the fragmentation potential for superheavy elements,", Romanian Journal of Physics , vol. 60, no. 1-2, pp. 147-155, 2015.

[18] J. Maruhn, W. Greiner, "The asymmetrie two center shell model,", Zeitschrift für Physik , vol. 251, no. 5, pp. 431-457, 1972.

[19] T. Dongjiao, Y. Ye, M. A. Dewanto, "Generalized Matrix Numerov Solutions to the Schrödinger Equation," 2014

[20] G. Avdelas, A. Konguetsof, T. E. Simos, "Generalization of Numerov's method for the numerical solution of the Schrodinger equation in two dimensions,", Computers and Chemistry , vol. 24, no. 5, pp. 577-584, 2000.

[21] R. del Río, C. Villegas-Blas, Spectral Theory of Schrödinger Operators , vol. 340, of Contemporary Mathematics, American Mathematical Society, Providence, RI, USA, 2004.

[22] W. Greiner, J. A. Maruhn, Nuclear Models , Springer, Berlin, Germany, 1996.

[23] S. Gomez, R. Romero, "Few-electron semiconductor quantum dots with Gaussian confinement,", Open Physics , vol. 7, no. 1, pp. 12-21, 2009.

[24] B. Pritychenko, A. A. Sonzogni, D. F. Winchell, V. V. Zerkin, R. Arcilla, T. W. Burrows, C. L. Dunford, M. W. Herman, V. McLane, P. Oblozinský, Y. Sanborn, J. K. Tuli, "Nuclear reaction and structure data services of the national nuclear data center,", Annals of Nuclear Energy , vol. 33, no. 4, pp. 390-399, 2006.

[25] E. N. D. File, B. Part, V. Version, National Nuclear Data Center , Brookhaven National Laboratory, Upton, NY, USA, 1990.