Consistent description of 12C and 16O using finite range three-body interaction
Consistent description of C and O has been a long standing problem of microscopic cluster models, where the wave function is fully antisymmetrized and the effective interaction is applied not between clusters but between nucleons. When the effective interaction is designed to reproduce the binding energy of O (four ), the binding energy of C (three ) becomes underbound by about 10 MeV. In the present study, by taking into account the coupling with the -coupling shell model components and utilizing Tohsaki interaction, which is phenomenological but has finite-range three-body interaction terms, we show that consistent understanding of these nuclei can be achieved. The original Tohsaki interaction gives small overbound of about 3 MeV for O, and this is improved by slightly modifying three-body Majorana exchange parameter. Also, the coupling with the -coupling shell model wave function strongly contributes to the increase of the binding energy of C. So far the application of Tohsaki interaction has been limited to nuclei, and here, we add Bartlet and Heisenberg exchange terms in the two-body interaction for the purpose of applying it to neutron-rich systems, and it is applied to He.
pacs:21.30.Fe, 21.60.Cs, 21.60.Gx, 27.20.+n
The nuclei C and O are typical light nuclei, and they have been extensively studied based on the cluster approaches (1). Since there is no bound nucleus with mass number 5 or 8, formation of C from three He nuclei ( clusters) is a key process of the nucleosynthesis. Here, the second state at MeV plays a crucial role, which is the second excited state of C and located just above the threshold energy to decay into three He nuclei (2). The existence of three state just at this energy is essential factor in the synthesis of various elements in starts. For also O, cluster structure has been shown to be extremely important. Although the ground state corresponds to the doubly closed shell of the shell model, this configuration can be also interpreted from four and C+ view points, if we take certain limit for the inter cluster distances. Also the first excited state of O at MeV, very close to the threshold to decay into C and He, can be interpreted as C+He cluster state (3), and low-lying cluster states just around the threshold are quite important in the synthesis of O in starts (4). Various cluster models have been proposed and successfully applied to these nuclei.
However consistent description of C and O has been a long standing problem of microscopic cluster models. Here the definition of the microscopic model is that the wave function is fully antisymmetrized and the effective interaction is applied not between clusters but between nucleons. When the effective interaction is designed to reproduce the binding energy of O (four ), the binding energy of C (three ) becomes underbound by about 10 MeV, and on the contrary, when the binding energy of C is reproduced, O becomes overbound by about 20 MeV. We have previously utilized Tohsaki interaction (5), which has finite-range three-body terms, and the obtained result is better than ones only within the two-body terms; however the problem has not been fully solved (6).
One of the clue to solve this problem is the inclusion of the spin-orbit interaction. In most of the cluster models, clusters are defined as simple configuration at some point. These clusters are spin singlet systems and the spin-orbit interaction does not contribute inside clusters and also between clusters. In -coupling shell model, the spin-orbit interaction is quite important and this plays an essential role in explaining the observed magic numbers. According to the -coupling shell model, C corresponds to the subclosure configuration of spin-orbit attractive orbits () and the spin-orbit interaction works attractively, whereas O corresponds to the closure of major shell ( shell), where both spin-orbit attractive () and repulsive () orbits are filled and the contribution of the spin-orbit interaction cancels. Therefore, inclusion of the breaking wave function and taking into account the spin-orbit contribution are expected to decrease the binding energy difference of C and O. To describe the -coupling shell model states with the spin-orbit contribution starting with the cluster model wave function, we proposed the antisymmetrized quasi-cluster model (AQCM) (7); (8); (9); (10); (11); (12); (13). In the AQCM, the transition from the cluster- to shell-model-structure can be described by only two parameters: representing the distance between clusters and , which characterizes the transition of cluster(s) to quasi-cluster(s) and quantifies the role of the spin-orbit interaction.
In nuclear structure calculations, it is quite well known that the central part of the nucleon-nucleon interaction in the calculation should have proper density dependence in order to to satisfy the saturation property of nuclear systems. If we just introduce simple two-body interaction, for instance Volkov interaction (14), which has been widely used in the cluster studies, we have to properly choose Majorana exchange parameter for each nucleus, and consistent description of two different nuclei with the same Hamiltonian becomes a tough work. Thus it is rather difficult to reproduce the threshold energies to decay into different subsystems.
Concerning the density dependence of the interaction, adding zero range three-body interaction term helps better agreements with experiments as in modified Volkov (MV) interaction (15). However, in this case the binding energies become quite sensitive to the choice of the size parameter of Gaussian-type single particle wave function. Especially, the binding energy and radius of He cannot be reproduced consistently. This situation is essentially common in the case of Gogny interaction widely used in mean field studies (16). The nucleus He is a building block of cluster states and it is desired that the size and binding energy are reproduced. The Tohsaki interaction, which has finite range three-body terms, has much advantages compared with the zero range three-body interactions. This interaction is a phenomenological one and designed to reproduce the - scattering phase shift. Also it gives reasonable size and binding energy of the cluster, which is rather difficult in the case of zero-range three-body interaction, and the binding energy is less sensitive to the choice of size parameter of Gaussian-type single particle wave function because of the finite range effect of the three-body interaction. Furthermore, the saturation properties of nuclear matter is reproduced rather satisfactory.
Of course, introducing the term proportional to the fractional power of the density is another possibility to reproduce the saturation properties of nuclear systems as in density functional theories (DFT), instead of introducing three-body interaction terms. However, here we perform parity and angular momentum projections, and we also superpose many Slater determinants based on the generator coordinate method (GCM). In this case, it is desired that the Hamiltonian is expressed in a operator form such as three-body interaction, which enables us to calculate the transition matrix elements between different Slater determinants. From this view point, simplified version of finite range three-body interaction is proposed in Ref. (17).
The purpose of the present work is to combine the use of finite range three-body interaction for the interaction part and AQCM for the wave function part to establish consistent understanding of C and O. The original Tohsaki interaction gives small overbound for O (about 3 MeV) (6), and here, we try to improve by slightly modifying three-body Majorana exchange parameter. For C, subclosure configuration of -coupling shell model, where the spin-orbit interaction plays an important role, is coupled to three model spaced based on AQCM. For O, the closed shell configuration of the shell is the dominant configuration of the ground state, and we apply four model, which coveres the model space of closed shell. Also, the application of Tohsaki interaction has been limited to nuclei, and we add Bartlet and Heisenberg exchange terms in the two-body interaction for the purpose of applying it to neutron-rich systems.
Ii The Model
The Hamiltonian () consists of kinetic energy () and potential energy () terms,
and the kinetic energy term is described as one-body operator,
and the center of mass kinetic energy (), which is constant, is subtracted. The potential energy has central (), spin-orbit (), and the Coulomb parts.
ii.2 Tohsaki Interaction
For the central part of the potential energy (), we adopt Tohsaki interaction. The Tohsaki interaction consists of two-body () and three-body () terms:
where and consist of three terms,
Here, represents the exchange of spatial part of the wave functions of interacting two nucleons, and this is equal to due to the Pauli principle (), where and are spin and isospin exchange operators, respectively. The range parameters are set to be common for the two-body and three-body parts, and the values are listed in Table 1 together with strengths of two-body interaction , three-body interaction , and the Majorana exchange parameters of the two-body interaction. The values of Winger parameters, , are given as . We employ F1 parameter set in Ref. (5).
Until now Tohsaki interaction has been applied only to nuclei, and in this article we extend the application to neutron-rich nuclei. The original Tohsaki interaction does have Wigner and Majorana exchange terms, but spin and isospin exchange terms are missing. Because of this, the interaction gives a weak bound state for a two neutron system, as in the case of Volkov interaction. Therefore, here we add Bartlet () and Heisenberg () exchange terms in Eq. 4 as . The values of and are chosen to be 0.1. By adding these terms, the neutron-neutron interaction (or pron-proton interaction) because weaker than the original interaction, while - scattering phase shift is not influenced by this modification.
ii.3 Spin-orbit interaction
For the spin-orbit part, G3RS (18), which is a realistic interaction originally determined to reproduce the nucleon-nucleon scattering phase shift, is adopted;
where fm fm, and is a projection operator onto a triplet odd state. The operator stands for the relative angular momentum and is the total spin, (). The strength, , has been determined to reproduce the He+ scattering phase shift (19), and MeV has been suggested.
ii.4 Single particle wave function (Brink model)
In conventional cluster models, the single particle wave function has a Gaussian shape (20);
where represents the spin-isospin part of the wave function, and is a real parameter representing the center of a Gaussian wave function for the -th particle. In this Brink-Bloch wave function, four nucleons in one cluster share the common value. Hence, the contribution of the spin-orbit interaction vanishes.
ii.5 Single particle wave function in the AQCM
In the AQCM, clusters are changed into quasi clusters. For nucleons in the quasi cluster, the single particle wave function is described by a Gaussian wave packet, and the center of this packet is a complex parameter;
where and in Eq. (9) represent the spin and isospin parts of the -th single particle wave function, respectively. The spin orientation is governed by the parameters and , which are in general complex, while the isospin part is fixed to be ‘up’ (proton) or ‘down’ (neutron),
The center of Gaussian wave packet is give as
where is a unit vector for the intrinsic-spin orientation, and is a real control parameter describing the dissolution of the cluster. As one can see immediately, the AQCM wave function, which has no imaginary part, is the same as the conventional Brink-Bloch wave function. The AQCM wave function corresponds to the -coupling shell model wave function, such as subshell closure configuration, when and . The mathematical explanation for this is summarized in Ref. (11). For the width parameter, we use the value of fm.
ii.6 AQCM wave function of the total system
The wave function of the total system is antisymmetrized product of these single particle wave functions;
The projections onto parity and angular momentum eigen states can be performed by introducing the projection operators and , and these are performed numerically in the actual calculation.
ii.7 Superposition of different configurations
Based on GCM, the superposition of different AQCM wave functions can be done,
Here, is a set of AQCM wave functions with different values of the and parameters, and the coefficients for the linear combination, are obtained by solving the Hill-Wheeler equation (20). We obtain a set of coefficients for the linear combination for each eigen value of .
We start the application of Tohsaki interaction with He. So far Tohsaki interaction has been applied to 4 nuclei, and here a neutron-rich system is examined. Figure 1 shows the energy convergence of the ground state of He. The model space is ++ and the positions of two neutrons are randomly generated, and we superpose different configurations for the neutrons. The horizontal line at MeV shows the threshold energy of He++ (the experimental value is MeV). Here, the dotted line is the result of original F1 parameter set. We add the Bartlet and Heisenberg terms with the parameters of , and the result is shown as the dashed line. Furthermore, we slightly modify the Majorana exchange parameter of the three-body interaction term and the results of F1’ parameter set is shown as solid line (this is to reproduce the binding energy of O, which will be discussed shortly). For the spin-orbit part, the strength has been determined to be MeV in the analysis of He+ scattering phase shift (19), and here we adopted MeV. The results of dashed and solid lines are similar and difficult to distinguish. Nevertheless, by adding Bartlet and Heisenberg terms, we obtained the binding energy close to the experimental onei (the experimental value of is 0.975 MeV).
For O, according to the -coupling shell model, the ground state corresponds to the closure of major shell ( shell), where both spin-orbit attractive () and repulsive () orbits are filled and the contribution of the spin-orbit interaction cancels. Therefore, breaking configurations are not expected to mix strongly, and here we introduce a four model space, which is known to coincide with the closed shell configuration at the limit of relative distance between clusters equal to zero. The energy of O with a tetrahedron configuration of four ’s is shown in Fig. 2 as a function of the relative distance between clusters. The dotted line is the result of original F1 parameter set, and the result of newly introduced F1’ parameter set is shown as the solid line. Here, F1’ is designed to avoid small overbinding of O when the original F1 parameter set is introduced, and the solid line is less attractive compared with the dotted line by about MeV. The energies of O is obtained by superposing the basis states with - distances of 0.1, 0.5, 1.0, 1.5, 2.0, 2.5, and 3.0 fm based on GCM, and the newly introduced parameter F1’ parameter set gives MeV for the ground state compared to the experimental value of MeV.
The observed root mean square (rms) radius of O is quite large; the charge radius is 2.69 fm (21), and we often underestimate the radius by 0.1-0.2 fm, if we calculate with four models and use only two-body effective interaction such as Volkov interaction (14). The dashed line in Fig. 2 is for the result calculated using Volkov No.2 interaction with . In this case, the energy curve is much shallower and feature of the curve is quite different from the Tohsaki interaction cases. The dashed line shows the energy minimum point around the - distance of 2 fm. Using the Tohsaki interaction with finite-range three-body interaction terms, the solid line in Fig. 2 shows that the lowest energy is obtained with the - distance of 2.5 fm, larger than the result of Volkov interaction (dashed line) by 0.5 fm. The matter rms radius for the state of O with a tetrahedron configuration of four ’s is shown in Fig. 3 as a function of distances between clusters. When the - distance is 2.5 fm, which gives the lowest energy in the Tohsaki interaction cases, the matter radius is 2.49 fm. This matter radius decreases to 2.35 fm if the - distance is 2.0 fm, which gives the lowest energy in the Volkov interaction case. The ground state of O obtained by superposing the basis states with the - distances of 0.1, 0.5, 1.0, 1.5, 2.0, 2.5, and 3.0 fm gives the rms matter radius of 2.49 fm (Tohsaki interaction F1’ parameter set). This value corresponds to the charge radius of 2.64 fm, and the experimental value is almost reproduced.
Here we compare the F1 (original) and F1’ (modified) parameter sets of the Tohsaki interaction. The energy curves of - (Be) calculated with F1 (dotted line) and F1’ (solid line) are compared in Fig. 4. As explained previously, F1’ is designed to avoid small overbinding of O calculated with F1, and the solid line is slightly more repulsive at short - relative distances. However the difference is quite small and less than 1 MeV, and the character of the original F1 that the - scattering phase shift is reproduced is not influenced by this modification.
For C, we prepare -coupling ( breaking) components of the wave functions based on AQCM. We introduce basis states with equilateral triangular configuration of three clusters with the relative distances of = 0.5, 1.0, 1.5, 2.0, 2.5, and 3.0 fm and change clusters to quasi clusters by giving dissolution parameter . The values of are chosen to be 0.2 and 0.4, since the states in between pure three clusters () and -coupling shell model limit () are known to be important for the description of the ground state. In addition to these 12 basis states, we prepare 28 basis sates with various there configurations by randomly generating Gaussian center parameters. This is because, the Hoyle state is a gas-like state without specific shape, and it has been known that not only equilateral triangular configuration, various three cluster configurations couple in this state.
By superposing these 40 basis states based on GCM and diagonalizing the Hamiltonian, energy eigen states are obtained. The F1’ parameter set of the Tohsaki interaction is adopted for the central part. In Fig 5, the ground , first , and second states of C are shown together with the calculated three threshold energy (dotted line). The strength of the spin-orbit interaction, in Eq. 7, is chosen as 0 MeV, 1600 MeV, 1800 MeV, and 2000 MeV in (a), (b), (c), and (d), respectively. The reasonable range of the strength of MeV has been suggested in the He+ scattering phase shift analysis (19). Without the spin-orbit interaction ( 0 MeV), the ground state of C is obtained at MeV in Fig 5 (a) compared with the experimental value of MeV. However, with the spin-orbit effect, the ground state is obtained at MeV in Fig 5 (b) ( 1600 MeV), MeV in (c) ( 1800 MeV), and MeV in (d) ( 2000 MeV). Therefore, the absolute value of the binding energy of C can be almost reproduced with the present interaction and the model, together with the binding energies of He and O. If we measure the energy from the three threshold energy, the ground state of C is MeV in Fig 5 (a), MeV in (b), MeV in (c), and MeV in (d), compared with the experimental value of MeV. Thus the binding energy from the three threshold is also reproduced in the case of MeV (Fig 5 (c)).
The famous Hoyle state, the second state experimentally observed at MeV, appears at MeV in Fig 5 (a), MeV in (b), MeV in (c), and MeV in (d). If we measure from the three threshold energy, these energies correspond to MeV in Fig 5 (a), MeV in (c), MeV in (c), and MeV in (d). Here again MeV gives reasonable agreements with the experiment; however the result implies that we need slightly larger number of basis states to reproduce Hoyle state just above (experimentally 0.38 MeV) the three threshold energy.
The traditional three cluster models have a serious problem that they give very small level spacing between the ground state and the first state, which is the first excited state of C (1). In our case, MeV result in Fig 5 (a) shows that the level spacing is 2.5 MeV and much smaller compared with the experimental value of 4.4389131 MeV. This is improved by the spin-orbit effect, since the excitation from the ground state to the state corresponds to one-particle one-hole excitation to a spin-orbit unfavored orbit from the closure configuration of spin-orbit attractive orbits in the -coupling shell model. The level spacing becomes 3.1 MeV in (b), 3.7 MeV in (c), and 4.6 MeV in (c). Similar trend has been reported in the recent antisymmetrized molecular dynamics (22) and fermionic molecular dynamics (23) calculations.
The appearance of negative parity states in low-lying excitation energy has been considered as the signature of the importance of three cluster structure in C. Experimentally and states have been observed at MeV and MeV, respectively. These and states are reproduced at 9.8 MeV and 12.6 MeV, respectively, when the strength of the spin-orbit interaction is chosen as MeV, which gives reasonable results for the and states.
We have tried to achieve consistent description of C and O, which has been a long standing problem of microscopic cluster model. By taking into account the coupling with the -coupling shell model and utilizing Tohsaki interaction, which is finite-range three-body interaction, we have shown that consistent understanding of these nuclei can be achieved. The original Tohsaki interaction gives small overbound of about 3 MeV for O, and this is improved by slightly modifying three-body Majorana exchange parameter. Also, so far the application of Tohsaki interaction has been limited to nuclei, and here, we added Bartlet and Heisenberg exchange terms in the two-body interaction for the purpose of applying it to neutron-rich systems.
By applying Tohsaki interaction with finite-range three-body interaction terms to O, the lowest energy of the tetrahedron configuration of four ’s is obtained with very large - distance (2.5 fm). After performing GCM, the ground state is obtained with the charge radius of 2.64 fm, compared with the observed value of 2.69 fm. We often underestimate the radius by 0.1-0.2 fm with four models, if we calculate only within the two-body effective interactions, and this is significantly improved.
For C, we prepared various there configurations by randomly generating Gaussian center parameters, and we mix -coupling ( breaking) components based on AQCM. The ground state of C is obtained at MeV with reasonable strength of the spin-orbit interaction compared with the experimental value of MeV. The absolute value of the binding energy of C (and also He and O) can be almost reproduced with the present interaction and the model. If we measure the energy from the three threshold energy, the agreement with the experiment is even more reasonable. The famous Hoyle state (second state) is reproduced just around the three threshold energy. Also, traditional three cluster models give very small level spacing for the ground state and the first state, and this is significantly improved by the spin-orbit effect. The appearance of negative parity states in low-lying excitation energy has been considered as the signature of the three cluster structure of C. Experimentally and states have been observed at MeV and MeV, respectively, and these states are also reproduced within the present framework.
Acknowledgements.The author thank the discussions with Prof. A. Tohsaki. Numerical calculation has been performed at Yukawa Institute for Theoretical Physics, Kyoto University. This work was supported by JSPS KAKENHI Grant Numbers 716143900002.
- Y. Fujiwara et al., Supple. of Prog. Theor. Phys. 68 29 (1980).
- F. Hoyle, D. N. F Dunbar, W. A. Wenzel, W. Whaling, Phys. Rev. 92, 1095c (1953).
- Hisashi Horiuchi and Kiyomi Ikeda, Prog. Theor. Phys. 40, 277 (1968).
- P. Descouvemont, Phys. Rev. C 47, 210 (1993).
- Akihiro Tohsaki, Phys. Rev. C 49, 1814 (1994).
- Naoyuki Itagaki, Akira Ohnishi, and Kiyoshi Katō, Prog. Theor. Phys. 94, 1019 (1995).
- N. Itagaki, H. Masui, M. Ito, and S. Aoyama, Phys. Rev. C 71 064307 (2005).
- H. Masui and N. Itagaki, Phys. Rev. C 75 054309 (2007).
- T. Yoshida, N. Itagaki, and T. Otsuka, Phys. Rev. C 79 034308 (2009).
- N. Itagaki, J. Cseh, and M. Płoszajczak, Phys. Rev. C 83, 014302 (2011).
- T. Suhara, N. Itagaki, J. Cseh, and M. Płoszajczak, Phys. Rev. C 87, 054334 (2013).
- Tadahiro Suhara and Yoshiko Kanada-En’yo, Phys. Rev. C 91, 024315 (2015).
- N. Itagaki, H. Matsuno, and T. Suhara, Prog. Theor. Exp. Phys., in press.
- A.B. Volkov, Nucl. Phys. 74, 33 (1965).
- T. Ando, K. Ikeda, and A. Tohsaki-Suzuki, Prog. Theor. Phys. 64, 1608 (1980)
- J. Decharge and D. Gogny, Phys. Rev. C 21, 1568 (1980).
- Y. Kanada-En‘yo and Y. Akaishi, Phys. Rev. C 69, 034306 (2004).
- R. Tamagaki, Prog. Theor. Phys. 39, 91 (1968).
- Shigeto Okabe and Yasuhisa Abe, Prog. Theor. Phys. 61 1049 (1979).
- D.M. Brink, in Proceedings of the International School of Physics ”Enrico Fermi” Course XXXVI, edited by C. Bloch (Academic, New York, 1966), p. 247.
- I. Angeli and K. P. Marinova, At. Data Nucl. Data Tables 99, 69 (2013).
- Y. Kanada-En’yo, Phys. Rev. Lett. 81, 5291 (1998); Y. Kanada-En’yo, Prog. Theor. Phys. 117, 655 (2007).
- T. Neff and H. Feldmeier, Nucl. Phys. A738, 357 (2004).