Research Article | Open Access

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

**Academic Editor:**Akpan Ikot

#### 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.

#### 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 to [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, 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]where ± are the positions of the two cluster centers and other parameters are shown in Figure 1.

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:where , , and are the potential depth, the parameters , , and are inverse square range, and is cluster centers. We study nucleus that has two identical cluster structures. Hence, we suppose that the parameters and have a symmetric potential form. In Figure 2, the potential is plotted for some values of . The nuclei with the closed shell such as 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 , the nucleons have been observed to cluster structure with a two-center structure. For example, 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).

**(a) fm**

**(b) fm**

**(c) fm**

**(d) fm**

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:where and denote the energy eigenvalue and the mass, respectively. is potential that does not depend on because of the axial symmetry and is separable to and . In the cylindrical coordinates, the Laplace operator takes the following form:

The separation of variables is accomplished by substituting

In the usual way, this leads to the following ordinary differential equations [18]:

The solution of the first equation iswhere 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 and are considered as and . The eigenvectors for three lowest states are shown in Figure 3 and the eigenvectors are plotted for some values of in Figure 4. is the radial quantum number and is the quantum number along -axis having values .

**(a) fm**

**(b) fm**

**(c) fm**

**(d) fm**

The eigenvectors vanish on the boundaries of our system and behave as expected. For the eigenvectors , in addition to boundary conditions, we consider the continuity conditions too:

While we describe nuclei in an axially symmetric model, the orbital angular momentum, , and the intrinsic spin, , are not good quantum numbers because states with different -values with the same parity can mix. In this model, the quantum number is the principal quantum number and the parity is determined by . The magnetic quantum number has the integer values between − and too.

Figures 4(a)–4(d) show that the wave functions associated with -axis have a definite parity with respect to . 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 -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 , there are the principal levels without spin-orbit coupling () and, for large values of , the levels start to form two clusters and their levels overlap with each other.

One of the interesting features of Figure 5 is that states having the same -component of angular momentum, , 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 and but the same value of are degenerate. For example, are degenerate with . 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]where is a dimensionless coupling constant; and are the nucleon and the speed of the light, respectively [15]. The orbital angular momentum operator is and we must transform it into cylindrical coordinates:hence,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 areand the spin-orbit matrix is Now, we use the parabolic approximation [23], and consequently the levels of single particle energies are obtained as a function of by the diagonalization of the potential . The structure of levels is plotted in Figure 6 and labeled them by . In our model, the projection of the total angular momentum on the symmetry axis, , is the good quantum number. Hence, states of different are not coupled by the Hamiltonian but the states with ± are degenerate because the nuclei have reflection symmetry for either of the two possible directions of the symmetry axis.

In Figure 6, the level scheme of the spherical nuclei is obtained for . But, by increasing , two separated clusters are formed. An example of light nuclei by two cluster structures is isotope. The nucleon single particle level scheme is plotted in Figure 6. 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 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 () and () overlap with each other. In our model, the single neutron is in the level at the ground state and its parity is . For the first excited state, the single neutron moves into level . 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 , the level scheme of the spherical nuclei is obtained, and two separated clusters are formed by increasing . 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.

#### Competing Interests

The authors declare that they have no competing interests.

#### References

- W. von Oertzen, M. Freer, and Y. Kanada-En'yo, “Nuclear clusters and nuclear molecules,”
*Physics Reports*, vol. 432, no. 2, pp. 43–113, 2006. View at: Publisher Site | Google Scholar - Y. Kanada-En'yo, T. Suhara, and F. Kobayashi, “Cluster formation, breaking, and excitation in light nuclei,” in
*Proceedings of the International Nuclear Physics Conference (INPC '13)*, Florence, Italy, June 2013. View at: Publisher Site | Google Scholar - M. Freer and 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. View at: Publisher Site | Google Scholar - M. R. Shojaei and 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. View at: Publisher Site | Google Scholar - Y. Kanada-En'yo and H. Horiuchi, “Structure of light unstable nuclei studied with antisymmetrized molecular dynamics,”
*Progress of Theoretical Physics Supplement*, no. 142, pp. 205–263, 2001. View at: Google Scholar - M. Freer, “Clusters in nuclei,”
*Scholarpedia*, vol. 5, no. 6, article 9652, 2010. View at: Publisher Site | Google Scholar - D. M. Brink, “History of cluster structure in nuclei,”
*Journal of Physics: Conference Series*, vol. 111, no. 1, Article ID 12001, 2008. View at: Publisher Site | Google Scholar - H. Feldmeier, “Fermionic molecular dynamics,”
*Nuclear Physics A*, vol. 515, no. 1, pp. 147–172, 1990. View at: Publisher Site | Google Scholar - A. Tohsaki, H. Horiuchi, P. Schuck, and G. Röpke, “Alpha cluster condensation in
^{12}C and^{16}O,”*Physical Review Letters*, vol. 87, no. 19, Article ID 192501, 4 pages, 2001. View at: Publisher Site | Google Scholar - Y. Funaki, A. Tohsaki, H. Horiuchi, P. Schuck, and 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, Article ID 051306, 2003. View at: Google Scholar - T. Yamada and P. Schuck, “Dilute multi-
*α*cluster states in nuclei,”*Physical Review C—Nuclear Physics*, vol. 69, no. 2, p. 24309, 2004. View at: Google Scholar - 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. View at: Publisher Site | Google Scholar - C. L. Zhang, B. Schuetrumpf, and W. Nazarewicz, “Nucleon localization in light and heavy nuclei,” https://arxiv.org/abs/1607.00422. View at: Google Scholar
- 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. View at: Publisher Site | Google Scholar - M. Mirea, “Two center shell model with Woods-Saxon potentials,”
*Romanian Reports in Physics*, vol. 59, no. 2, pp. 523–531, 2007. View at: Google Scholar - J.-P. Ebran, E. Khan, T. Nikšić, and D. Vretenar, “How atomic nuclei cluster,”
*Nature*, vol. 487, no. 7407, pp. 341–344, 2012. View at: Publisher Site | Google Scholar - D. Aranghel and A. Sandulescu, “Shell effects in the fragmentation potential for superheavy elements,”
*Romanian Journal of Physics*, vol. 60, no. 1-2, pp. 147–155, 2015. View at: Google Scholar - J. Maruhn and W. Greiner, “The asymmetrie two center shell model,”
*Zeitschrift für Physik*, vol. 251, no. 5, pp. 431–457, 1972. View at: Publisher Site | Google Scholar - T. Dongjiao, Y. Ye, and M. A. Dewanto, “Generalized Matrix Numerov Solutions to the Schrödinger Equation,” 2014. View at: Google Scholar
- G. Avdelas, A. Konguetsof, and 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. View at: Publisher Site | Google Scholar | Zentralblatt MATH - R. del Río and C. Villegas-Blas, Eds.,
*Spectral Theory of Schrödinger Operators*, vol. 340 of*Contemporary Mathematics*, American Mathematical Society, Providence, RI, USA, 2004. View at: Publisher Site | MathSciNet - W. Greiner and J. A. Maruhn,
*Nuclear Models*, Springer, Berlin, Germany, 1996. View at: Publisher Site - S. Gomez and R. Romero, “Few-electron semiconductor quantum dots with Gaussian confinement,”
*Open Physics*, vol. 7, no. 1, pp. 12–21, 2009. View at: Publisher Site | Google Scholar - B. Pritychenko, A. A. Sonzogni, D. F. Winchell et al., “Nuclear reaction and structure data services of the national nuclear data center,”
*Annals of Nuclear Energy*, vol. 33, no. 4, pp. 390–399, 2006. View at: Publisher Site | Google Scholar - E. N. D. File, B. Part, and V. Version,
*National Nuclear Data Center*, Brookhaven National Laboratory, Upton, NY, USA, 1990.

#### Copyright

Copyright © 2017 Nafiseh Roshanbakht and Mohammad Reza Shojaei. 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^{3}.