Get 20M+ Full-Text Papers For Less Than $1.50/day. Start a 14-Day Trial for You or Your Team.

Learn More →

Implementation of the 8-Nucleon Yakubovsky Formalism for Halo Nucleus 8He

Implementation of the 8-Nucleon Yakubovsky Formalism for Halo Nucleus 8He In order to study the bound-state structure of the Helium halo nuclei, the 8-nucleon Yakubovsky formalism has been implemented for He in a 5-body sub-cluster model, i.e. 𝛼 + 𝑛 + 𝑛 + 𝑛 + 𝑛 . In this case, the 8- nucleon Yakubovsky equations has been obtained in the form of two coupled equations, based on the two independent components. In addition, by removing the contribution interactions of the 8 and 7’s bound nucleons in the formalism, the obtained equations explicitly reduce to the 6-nucleon Yakubovsky equations for He, in the case of effective 3-body model, i.e. 𝛼 + 𝑛 + 𝑛 . In view of the expectation for the dominant structure of He, namely an inert 𝛼 -core and four loosely-bound neutrons, Jacobi configurations of the two components in momentum space have been represented to provide technicalities which were considered useful for a numerical performance, such as bound-state calculations and momentum density distributions for halo-bound neutrons. PACS. No: 21.45.+v Few-body systems; 21.30.Fe Forces in hadronic systems and effective interactions; 21.30.-x Two-nucleon interactions; 21.10.Dr Binding energies and masses. Introduction Studies of rich-neutron light nuclei present stringent tests for their calculations as they probe aspects of the interactions that are less prevalent in nuclei closer to stability. Among them, the physics of strong interactions gives rise to new structures in rich-neutron light nuclei. One of the best prominent cases are the helium halo nuclei, namely 6 8 4 He and He with two and four loosely-bound neutrons, forming an extended halo around the He core. In this regard, 6 8 the rich-neutron helium isotopes, He and He, are the topic of this study. Indeed, both nuclei are Borromean halos [1, 2], they have no bound excited states, and they exhibit the binding energy. Both helium halo nuclei are radioactive and decay quickly that complicates efforts to measure their properties. He has a time-life of only a tenth of a second and in this little time, can be considered as a stable nucleus. Practically, detection of the rare He is a major step for particle physicists. But, it is very important and interesting to understand how helium configures itself with strong interactions after being produced from a particle accelerator. In recent decades, combination of nuclear and atomic physics methods has allowed new courses of precision 6 8 measurements of the ground-state energies and charge radii of He [3, 4] and He [5, 6]. Dealing with the theoretical study of these nuclei is very challenging, because one needs to describe the small separation energy of the halo neutrons and the large radius of the whole system, in addition to the challenges of binding energy calculations. Because they are light-mass nuclei, one can use traditional methods to study them, such as Green’s Function Monte Carlo (GFMC) method [7] and the No-Core Shell Model (NCSM) [8]. In addition, some of the standard methods in nuclear theory has been the progress of Effective Field Theory (EFT) and the Renormalization Group (RG) to nuclear forces. While such helium halo nuclei have been investigated using the NCSM [9], but still, no sense has been obtained for the halo structure properties in such theoretical studies. 8 6 In this paper, to approach a simple physical structure of the observed effects, namely He and He halo-bound nuclei, and to understand the physics that governs the neutron distributions in helium halos, it was felt that by using current/modern computers with high computational speed, beyond previous works [10-14], a further step to study the light halo nuclei, namely He as an 8-nucleon bound system problem, within the reliable Yakubovsky method [15] is now desirable. It is worth mentioning that our obtained results are more difficult to implement bound-state problem E.Ahmady.ph@ut.ac.ir A.A.Rajabi@shahroodut.ac.ir in few-body branch, due to its halo nature and high-dimensional calculations. However, we provide technicalities for numerical performances. This paper is written as follows: In sect. I, the Yakubovsky formalism to the 8-nucleon bound system problem using the standard notation [16] is explicitly derived. In this case, the Lippmann-Schwinger scheme twice step-by- step has been applied on the integral form of Schrödinger equation, and then the identity of the nucleons are implemented. Finally, the 8-nucleon Yakubovsky formalism for effective 𝛼 -core structure in a 5-body sub-cluster model of He leads to a set of two coupled equations related to 2 different components. In sect. II, useful technicalities for numerical applications have been provided. To this aim, corresponding Jacobi configurations of each Yakubovsky component is schematically represented and then the integral form of coupled equations are represented by introducing basis states based on Jacobi momenta in partial-wave (PW) analysis. Finally, a summary in sect. III is given. I. The 8-nucleon Yakubovsky Formalism In branch of few-body systems, each N-body bound system has 𝑁 (𝑁 − 1)/2 different two-body interactions. Therefore, in the 8-nucleon bound system there are 28 different two-nucleon interactions or 28 different cluster decompositions (𝑎 ) having 7 clusters. They are indexed by the only pairing forces in cluster 𝑎 they contain, e.g. 𝑎 = 7 7 7 12 ≡ 12+ 3 + 4 + 5+ 6+ 7 + 8. In order to study the 8-nucleon bound system in the framework of the Yakubovsky method using the sub-cluster notation [16], the method is to first sum up the pair forces in each 7-body fragment (𝑎 ), in a second step among all 6-body fragments (𝑎 ), and then in a third step among all 5-body fragments (𝑎 ). We cut 6 5 off that formalism ending with 5-body sub-clusters in the spirit of the usually used approximating effective 5-body model that is 𝛼 + 𝑛 + 𝑛 + 𝑛 + 𝑛 . In the standard formalism, for an individual 𝛼 -particle, such two-body interactions, namely 𝛼 − 𝑛 and 𝛼 -clusters have not been used, because in this framework 𝛼 -core interacts as a four-nucleon sub- system in the 8-nucleon system. We start with the non-relativistic Schrödinger equation for the 8-nucleon bound system as (1.1) (𝐻 + ∑𝑉 )𝛹 = 𝐸 𝛹 , 0 𝑎 where 𝐻 refers to the free Hamiltonian operator of the 8-nucleon system that will be introduced in the next section, and ∑ 𝑉 ≡ 𝑉 + ⋯+ 𝑉 is the summation of all paring interactions with 28 terms. According to the Lipmann- 𝑎 𝑎 12 78 7 7 Schwinger scheme, Eq. (1.1) is rewritten into an integral equation (1.2) 𝛹 = 𝐺 ∑𝑉 𝛹 , 0 𝑎 where 𝐺 stand for the 8-nucleon free Green’s function operator and in the case of bound states formalism we have −1 𝐺 = [𝐸 − 𝐻 ] . Like previous works for five- and six-nucleon Yakubovsky formalism [12, 13], we start from the 0 0 integral form of the Schrödinger equation, Eq. (1.2), applying Lipmann-Schwinger scheme twice step-by-step. We can define 𝜓 ≡ 𝐺 𝑉 𝛹 by using Eqs. (1.2) and (1.7). Therefore, Faddeev like equation for 8-nucleon Yakubovsky 𝑎 0 𝑎 7 7 amplitudes are given as (1.3) 𝜓 ≡ 𝐺 𝑡 ∑𝛿 𝜓 , 𝑎 0 𝑎 𝑎 𝑏 𝑏 7 7 7 7 7 where 𝑡 is a pairing 𝑡 -matrix operator that follows the Lippmann-Schwinger equation as 𝑡 = 𝑉 + 𝑉 𝐺 𝑡 and 𝑎 𝑎 𝑎 𝑎 0 𝑎 7 7 7 7 7 𝛿 = 1 − 𝛿 . For the next step sub-cluster form the 6-body fragments amplitudes can be written 𝑎 𝑏 𝑎 𝑏 7 7 7 7 𝑎 𝑎 6 𝑎 6 6 ̅ (1.4) 𝜓 = 𝐺 ∑ 𝒯 (𝜓 ) =𝐺 ∑ 𝒯 ∑ 𝛿 𝜓 , 𝑎 𝑎 0 𝑎 0 𝑎 𝑏 𝑏 𝑏 7 6 𝑎 𝑏 7 𝑎 𝑏 6 6 7 6 7 7 7 7 𝑏 ⊂𝑎 𝑏 ⊂𝑎 𝑏 ⊂𝑏 7 6 7 6 7 6 where 𝒯 obeys as 𝑎 𝑏 7 7 (1.5) 𝑎 𝑎 6 6 ̅ ̅ 𝒯 = 𝑡 𝛿 + 𝐺 ∑ 𝑡 𝛿 𝒯 , 𝑎 𝑏 𝑎 𝑎 𝑏 0 𝑎 𝑎 𝑐 𝑐 𝑏 7 7 7 7 7 7 7 7 7 7 𝑐 ⊂𝑎 7 6 next, more decompose the right side of Eq. (1.4) according to 5-body fragments are given as 𝑎 𝑎 5 6 (1.6) 𝜓 = 𝐺 ∑ 𝒯 ∑ 𝛿 𝜓 . 0 𝑎 𝑏 𝑏 𝑏 𝑎 ,𝑎 𝑎 𝑏 7 6 6 6 7 6 7 7 𝑏 ⊂𝑏 𝑏 ⊂𝑎 7 6 7 6 𝑏 ⊂𝑎 6 5 briefly, we work out that method ending with (𝑎 ) in the spirit of the applied though approximate effective 𝛼 + 𝑛 + 𝑛 + 𝑛 + 𝑛 in a 5-body model. Further details could be found in Refs. [12, 13]. In the following, the single nucleons in sub-clusters, namely 7-, 6- and 5-body fragments, i.e. (𝑎 ), (𝑎 ) and (𝑎 ) respectively, will no longer be indexed. 7 6 5 Next, in order to get to 8-nucleon Yakubovsky coupled equations, the identity of the nucleons has been applied. In addition, the number of huge total WF Yakubovsky components has been obtained. It is well-known, the 8-nucleon total WF consists of the summation of all pairing Yakubovsky components with 28 terms which are as 𝛹 = ∑𝜓 ≡ 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 𝑎 12 13 14 15 16 17 18 +𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 23 24 25 26 27 28 +𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 34 35 36 37 38 +𝜓 + 𝜓 + 𝜓 + 𝜓 45 46 47 48 (1.7) +𝜓 + 𝜓 + 𝜓 56 57 58 +𝜓 + 𝜓 67 68 +𝜓 , where the first pair (𝜓 ), having 7-body fragments, consists of 21 components of 6-body fragments (𝑎 ) as 12 6 𝜓 = ∑ 𝜓 𝑎 ≡12 𝑎 ≡12,𝑎 7 7 6 (1.8) 𝑎 ≡12⊂𝑎 7 6 ≡ 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 12,123 12,124 12,125 12,126 12,127 12,128 12,12+34 12,12+35 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 12,12+36 12,12+37 12,12+38 12,12+45 12,12+46 12,12+47 12,12+48 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 , 12,12+56 12,12+57 12,12+58 12,12+67 12,12+68 12,12+78 where (𝑎 ) refers to any 6-body fragments containing the pair restricted to (𝑎 ≡12) and the sum runs over pairs 6 7 𝑎 ⊂ 𝑎 . It means that the sub-clusters (𝑎 ), when broken up, lead to the sub-clusters (𝑎 ). Next first term as 𝜓 7 6 6 7 12,123 is a sub-cluster of 15 components of 5-body fragments (𝑎 ) 𝜓 = ∑ 𝜓 𝑎 ≡12,𝑎 ≡123 7 6 𝑎 ≡12,𝑎 ≡123 7 6 (1.9) 𝑎 ≡123⊂𝑎 6 5 1234 1235 1236 1237 1238 123+45 123+46 123+47 ≡ 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 12,123 12,123 12,123 12,123 12,123 12,123 12,123 12,123 123+48 123+56 123+57 123+58 123+67 123+68 123+78 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 , 12,123 12,123 12,123 12,123 12,123 12,123 12,123 the summation of all 5-body fragments ends up to the 8-nucleon Yakubovsky components for total WF as 𝑎 𝑎 5 5 (1.10) 𝛹 = ∑∑∑𝜓 = 8,820 𝜓 , 𝑎 ,𝑎 𝑎 ,𝑎 7 6 7 6 𝑎 𝑎 𝑎 7 6 5 We find that there are 8,820 components for the total WF of constituent 8 nucleons in a 5-body model. It is worthwhile to mention that to calculate the expectation value of the Hamiltonian operator, the above-mentioned total WF’s have to be used. Clearly, the total WF is anti-symmetrized following the Pauli principle. In the next step, the identity of the nucleons in sub-cluster components are implemented. We start from Eq. (1.4) choosing the case 7-body fragments 𝑎 ≡ 12 with 6-body fragments 𝑎 ≡ 123 7 6 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 12,124 12,125 12,126 12,127 12,128 +𝜓 +𝜓 +𝜓 +𝜓 +𝜓 12,12+34 12,12+35 12,12+36 12,12+37 12,12+38 +𝜓 +𝜓 +𝜓 +𝜓 12,12+45 12,12+46 12,12+47 12,12+48 𝜓 = 𝐺 𝒯 12,123 0 12,12 +𝜓 +𝜓 +𝜓 12,12+56 12,12+57 12,12+58 +𝜓 +𝜓 12,12+67 12,12+68 +𝜓 ( ) 12,12+78 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 23,234 23,235 23,236 23,237 23,238 +𝜓 +𝜓 +𝜓 +𝜓 +𝜓 23,23+14 23,23+15 23,23+16 23,23+17 23,23+18 +𝜓 +𝜓 +𝜓 +𝜓 23,23+45 23,23+46 23,23+47 23,23+48 +𝐺 𝒯 0 12,23 +𝜓 +𝜓 +𝜓 23,23+56 23,23+57 23,23+58 +𝜓 +𝜓 23,23+67 23,23+68 +𝜓 ( ) 23,23+78 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 31,314 31,315 31,316 31,317 31,318 +𝜓 +𝜓 +𝜓 +𝜓 +𝜓 31,31+24 31,31+25 31,31+26 31,31+27 31,31+28 +𝜓 +𝜓 +𝜓 +𝜓 31,31+45 31,31+46 31,31+47 31,31+48 +𝐺 𝒯 , 12,31 +𝜓 +𝜓 +𝜓 31,31+56 31,31+57 31,31+58 +𝜓 +𝜓 31,31+67 31,31+68 +𝜓 ( ) 31,31+78 (1.11) it is easily seen that, by using adequate permutation operators, the huge equation above can be rewritten as 123 123 123 𝜓 = 𝐺 (𝒯 + 𝒯 𝑃 𝑃 + 𝒯 𝑃 𝑃 ) 12,123 0 12,12 12,23 12 23 12,31 13 23 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 12,124 12,125 12,126 12,127 12,128 +𝜓 +𝜓 +𝜓 +𝜓 +𝜓 12,12+34 12,12+35 12,12+36 12,12+37 12,12+38 +𝜓 +𝜓 +𝜓 +𝜓 12,12+45 12,12+46 12,12+47 12,12+48 × , +𝜓 +𝜓 +𝜓 (1.12) 12,12+56 12,12+57 12,12+58 +𝜓 +𝜓 12,12+67 12,12+68 +𝜓 ( ) 12,12+78 by defining the bellow transition operator as 123 123 123 123 𝒯 = 𝒯 + 𝒯 𝑃 𝑃 + 𝒯 𝑃 𝑃 , (1.13) 12,12 12,23 12 23 12,31 13 23 Eq. (1.12) is simplified to 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 12,124 12,125 12,126 12,127 12,128 +𝜓 +𝜓 +𝜓 +𝜓 +𝜓 12,12+34 12,12+35 12,12+36 12,12+37 12,12+38 +𝜓 +𝜓 +𝜓 +𝜓 12,12+45 12,12+46 12,12+47 12,12+48 𝜓 = 𝐺 𝒯 . 12,123 0 +𝜓 +𝜓 +𝜓 12,12+56 12,12+57 12,12+58 +𝜓 +𝜓 (1.14) 12,12+67 12,12+68 +𝜓 ( ) 12,12+78 123 123 123 where 𝒯 follows Lipmann-Schwinger equation as 𝒯 = 𝑡 𝑃 + 𝑡 𝑃 𝐺 𝒯 and 𝑃 = 𝑃 𝑃 + 12 12 0 12 23 𝑃 𝑃 .Subsequently, starting again from Eq. (1.4) to select the case 7-body fragments 𝑎 ≡ 12 with 6-body fragments 13 23 7 𝑎 ≡ 12+ 34 𝜓 + 𝜓 +𝜓 +𝜓 +𝜓 +𝜓 12,123 12,124 12,125 12,126 12,127 12,128 +𝜓 +𝜓 +𝜓 +𝜓 12,12+35 12,12+36 12,12+37 12,12+38 +𝜓 +𝜓 +𝜓 +𝜓 12,12+45 12,12+46 12,12+47 12,12+48 12+34 𝜓 = 𝐺 𝒯 12,12+34 0 12,12 +𝜓 +𝜓 +𝜓 12,12+56 12,12+57 12,12+58 +𝜓 +𝜓 12,12+67 12,12+68 +𝜓 ( ) 12,12+78 𝜓 + 𝜓 +𝜓 +𝜓 +𝜓 +𝜓 34,134 34,234 34,345 34,346 34,347 34,348 +𝜓 +𝜓 +𝜓 +𝜓 34,34+15 34,34+16 34,34+17 34,34+18 +𝜓 +𝜓 +𝜓 +𝜓 34,34+25 34,34+26 34,34+27 34,34+28 12+34 + 𝐺 𝒯 , 0 12,34 +𝜓 +𝜓 +𝜓 34,34+56 34,34+57 34,34+58 +𝜓 +𝜓 34,34+67 34,34+68 (1.15) +𝜓 ( ) 34,34+78 subsequently, by using permutation properties, Eq. (1.15) can be rewritten as 𝜓 + 𝜓 +𝜓 +𝜓 +𝜓 +𝜓 12,123 12,124 12,125 12,126 12,127 12,128 +𝜓 +𝜓 +𝜓 +𝜓 12,12+35 12,12+36 12,12+37 12,12+38 +𝜓 +𝜓 +𝜓 +𝜓 12,12+45 12,12+46 12,12+47 12,12+48 12+34 12+34 𝜓 = 𝐺 (𝒯 + 𝒯 𝑃 𝑃 ) , 12,12+34 0 12,12 12,34 13 24 +𝜓 +𝜓 +𝜓 12,12+56 12,12+57 12,12+58 +𝜓 +𝜓 12,12+67 12,12+68 (1.16) +𝜓 ( ) 12,12+78 then defining 12+34 12+34 12+34 𝒯 = 𝒯 + 𝒯 𝑃 𝑃 , (1.17) 13 24 12,12 12,34 Eq. (1.16) is simplified to 𝜓 + 𝜓 +𝜓 +𝜓 +𝜓 +𝜓 12,123 12,124 12,125 12,126 12,127 12,128 +𝜓 +𝜓 +𝜓 +𝜓 12,12+35 12,12+36 12,12+37 12,12+38 +𝜓 +𝜓 +𝜓 +𝜓 12,12+45 12,12+46 12,12+47 12,12+48 12+34 𝜓 = 𝐺 𝒯 . 12,12+34 0 +𝜓 +𝜓 +𝜓 12,12+56 12,12+57 12,12+58 +𝜓 +𝜓 12,12+67 12,12+68 (1.18) +𝜓 ( ) 12,12+78 12+34 12+34 12+34 ̃ ̃ ̃ where 𝒯 follows the equation 𝒯 = 𝑡 𝑃 + 𝑡 𝑃 𝐺 𝒯 and where 𝑃 = 𝑃 𝑃 . 12 12 0 13 24 The next step is decomposing 𝜓 according to Eq. (1.6). For 7-body fragments 𝑎 ≡ 12 with 6-body 12,123 7 fragments 𝑎 ≡ 123 the all relevant 5-body fragments (𝑎 ) are 1234,1235,1236,1237,1238, 123+ 45,123+ 6 5 46,123+ 47,123+ 48,123+ 56,123+ 57,123+ 58,123+ 67,123+ 68,123+ 78. All 5-body sub-clusters summed up 1234 1235 1236 1237 1238 123+45 123+46 123+47 123+48 𝜓 = 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 +𝜓 + 𝜓 (1.19) 12,123 12,123 12,123 12,123 12,123 12,123 12,123 12,123 12,123 12,123 123+56 123+57 123+58 123+67 123+68 123+78 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 . 12,123 12,123 12,123 12,123 12,123 12,123 Considering Eq. (1.6), first part of the above turns 1234 123 123 (1.20) 𝜓 = 𝐺 𝒯 (𝜓 +𝜓 )+ 𝐺 𝒯 (𝜓 +𝜓 ) 12,123 0 12,12 12,124 12,12+34 0 12,23 23,234 23,23+14 + 𝐺 𝒯 (𝜓 +𝜓 ), 0 12,31 31,134 31,31+24 since (1.21) 𝜓 +𝜓 = 𝑃 𝑃 (𝜓 +𝜓 ), 23,234 23,23+14 12 23 12,124 12,12+34 (1.22) 𝜓 +𝜓 = 𝑃 𝑃 (𝜓 +𝜓 ), 31,134 31,31+24 13 23 12,124 12,12+34 according to Eq. (1.13), Eq. (1.20) leads to 1234 123 (1.23) 𝜓 = 𝐺 𝒯 (𝜓 +𝜓 ). 12,123 0 12,124 12,12+34 Similarly, next we decompose 𝜓 according to Eq. (1.6). For the 7-body fragments 𝑎 ≡ 12 with 6- 12,12+34 7 body fragments 𝑎 ≡ 12+ 34 the possible 5-body fragments (𝑎 ) are 1234,125+ 34,126+ 34,127+ 34,128+ 6 5 34,12+ 345,12+ 346,12+ 347,12+ 348,12+ 34+ 56,12+ 34+ 57,12+ 34+ 58,12+ 34+ 67,12+ 34+ 68,12+ 34+ 78, which are now summarized to 1234 123+45 123+46 123+47 123+48 12+345 (1.24) 𝜓 = 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 12+346 12+348 12+34+56 12+34+57 12+34+58 12+347 + 𝜓 +𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 12+34+67 12+34+68 12+34+78 + 𝜓 + 𝜓 + 𝜓 . 12,12+34 12,12+34 12,12+34 Considering Eq. (1.6), the first part of the above equation is regarded in turn 1234 12+34 12+34 (1.25) 𝜓 = 𝐺 𝒯 (𝜓 +𝜓 )+ 𝐺 𝒯 (𝜓 +𝜓 ), 0 12,123 12,124 0 34,234 34,134 12,12+34 12,12 12,34 since we use 𝜓 +𝜓 = 𝑃 𝑃 (𝜓 +𝜓 ). Eq. (1.25) is simplified according to Eq. (1.17), and 34,234 34,134 13 24 12,123 12,124 rewritten as 1234 12+34 (1.26) 𝜓 = 𝐺 𝒯 (𝜓 +𝜓 ). 0 12,123 12,124 12,12+34 1234 1234 The two amplitudes, 𝜓 and 𝜓 , expressed in Eq. (1.23) and Eq. (1.26) are coupled to each other 12,123 12,12+34 as shown now. The expression Eq. (1.19) can easily be converted to 𝜓 and using in addition Eq. (1.24), with 12,124 usage of Eq. (1.23), one obtains 1234 123 1234 1245 1246 1247 1248 124+35 124+36 124+37 𝜓 = 𝐺 𝒯 ((𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 +𝜓 12,123 0 12,124 12,124 12,124 12,124 12,124 12,124 12,124 12,124 124+38 124+36 124+57 124+58 124+67 124+68 124+78 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 ) 12,124 12,124 12,124 12,124 12,124 12,124 12,124 123+45 123+46 123+48 12+345 1234 123+47 + (𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 12+346 12+347 12+348 12+34+56 12+34+57 12+34+58 + 𝜓 +𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 12+34+67 12+34+68 12+34+78 (1.27) + 𝜓 + 𝜓 + 𝜓 )), 12,12+34 12,12+34 12,12+34 correspondingly Eq. (1.26) yields 1234 12+34 1234 1235 1236 1237 1238 123+45 𝜓 = 𝐺 𝒯 ((𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 12,12+34 0 12,123 12,123 12,123 12,123 12,123 12,123 123+46 123+47 123+48 123+56 123+57 123+58 123+67 + 𝜓 +𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 12,123 12,123 12,123 12,123 12,123 12,123 12,123 123+68 123+78 + 𝜓 + 𝜓 ) 12,123 12,123 1234 1245 1246 1247 1248 124+35 124+36 124+37 + (𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 +𝜓 12,124 12,124 12,124 12,124 12,124 12,124 12,124 12,124 (1.28) 124+38 124+36 124+57 124+58 124+67 124+68 124+78 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 )), 12,124 12,124 12,124 12,124 12,124 12,124 12,124 1234 1234 one separates now the components 𝜓 and 𝜓 in Eq. (1.27) from the rest 12,124 12,12+34 1234 123 1234 1234 𝜓 − 𝐺 𝒯 (𝜓 + 𝜓 ) 12,123 0 12,124 12,12+34 123 1245 1246 1247 1248 124+35 124+36 124+37 = 𝐺 𝒯 (𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 +𝜓 0 12,124 12,124 12,124 12,124 12,124 12,124 12,124 124+38 124+36 124+57 124+58 124+67 124+68 124+78 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 12,124 12,124 12,124 12,124 12,124 12,124 12,124 123+45 123+46 123+47 123+48 12+345 12+346 12+347 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 +𝜓 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 12+348 12+34+56 12+34+57 12+34+58 12+34+67 12+34+68 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 (1.29) 12+34+78 + 𝜓 ), 12,12+34 1234 1234 also, one separates now the components 𝜓 and 𝜓 in Eq. (1.28) from the rest 12,123 12,124 1234 12+34 1234 1234 𝜓 − 𝐺 𝒯 (𝜓 + 𝜓 ) 12,12+34 12,123 12,124 12+34 1235 1236 1237 1238 123+45 123+46 123+47 = 𝐺 𝒯 (𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 +𝜓 0 12,123 12,123 12,123 12,123 12,123 12,123 12,123 123+48 123+56 123+57 123+58 123+67 123+68 123+78 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 12,123 12,123 12,123 12,123 12,123 12,123 12,123 1245 1246 1248 124+35 124+36 124+38 1247 124+37 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 +𝜓 + 𝜓 12,124 12,124 12,124 12,124 12,124 12,124 12,124 12,124 124+56 124+57 124+58 124+67 124+68 124+78 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 ), 12,124 12,124 12,124 12,124 12,124 12,124 (1.30) 1234 1234 with 𝜓 = −𝑃 𝜓 Eqs. (1.29) and (1.30) into a matrix form 12,124 34 12,123 1234 1234 123 123 𝜓 𝒯 (−𝑃 ) 𝒯 𝜓 12;123 12;123 ( ) − 𝐺 ( )( ) 1234 12+34 1234 𝜓 𝒯 (1 − 𝑃 ) 0 𝜓 12;12+34 12;12+34 1245 1246 1248 124+35 124+36 124+38 1247 124+37 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 +𝜓 + 𝜓 12,124 12,124 12,124 12,124 12,124 12,124 12,124 12,124 124+36 124+57 124+58 124+67 124+68 124+78 +𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 12,124 12,124 12,124 12,124 12,124 12,124 123+45 123+46 123+47 123+48 12+345 12+346 12+347 12+348 +𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 +𝜓 + 𝜓 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 (1.31) 12+34+56 12+34+57 12+34+58 12+34+67 12+34+68 12+34+78 +𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 ( 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 ) = 𝐺 , 1235 1236 1237 1238 123+45 123+46 123+47 123+48 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 +𝜓 + 𝜓 12,123 12,123 12,123 12,123 12,123 12,123 12,123 12,123 123+56 123+57 123+58 123+67 123+68 123+78 +𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 12,123 12,123 12,123 12,123 12,123 12,123 12+34 1245 1246 1247 1248 124+35 124+36 124+37 124+38 +𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 +𝜓 + 𝜓 12,124 12,124 12,124 12,124 12,124 12,124 12,124 12,124 124+56 124+57 124+58 124+67 124+68 124+78 +𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 ( ) 12,124 12,124 12,124 12,124 12,124 12,124 ( ) since 1245 1235 124+35 123+45 (1.32) 𝜓 = −𝑃 𝜓 ,𝜓 = −𝑃 𝜓 , 12,124 34 12,123 12,124 34 12,123 the right side of Eq. (1.31) can be factored and achieves the form as 1234 1234 123 123 123 123 𝜓 𝜓 𝒯 (−𝑃 ) 𝒯 𝒯 (−𝑃 ) 𝒯 12;123 12;123 34 34 ( ) − 𝐺 ( )( ) = 𝐺 ( ) 0 0 1234 12+34 1234 12+34 𝜓 𝒯 (1 − 𝑃 ) 0 𝜓 𝒯 (1 − 𝑃 ) 0 12;12+34 34 12;12+34 34 1235 1236 1237 1238 123+45 123+46 123+47 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 12,123 12,123 12,123 12,123 12,123 12,123 12,123 ( ) 123+48 123+56 123+57 123+58 123+67 123+68 123+78 (1.33) +𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 12,123 12,123 12,123 12,123 12,123 12,123 12,123 × , 123+45 123+46 123+47 123+48 12+345 12+346 12+347 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 +𝜓 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 ( ) 12+348 12+34+56 12+34+57 12+34+58 12+34+67 12+34+68 12+34+78 +𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 ( ) 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 the right side can be reduced applying permutation operators and obtains the final form of the 8-nucleon equations as 123 123 𝒯 (−𝑃 ) 𝒯 12;123 (1.34) ( ) = 𝐺 ( ) 1234 12+34 𝜓 𝒯 (1 − 𝑃 ) 0 12;12+34 34 1234 123+45 (1 − 𝑃 −𝑃 −𝑃 −𝑃 )𝜓 + (1 − 𝑃 −𝑃 −𝑃 −𝑃 −𝑃 −𝑃 +𝑃 𝑃 +𝑃 𝑃 +𝑃 𝑃 )𝜓 45 46 47 48 12,123 56 57 58 46 47 48 46 57 46 58 47 57 12,123 × [( )], 1234 125+34 12+345 12+34+56 ( ) 𝜓 + 1 − 𝑃 −𝑃 −𝑃 ( 𝜓 + 𝜓 )+ (1 − 𝑃 −𝑃 −𝑃 −𝑃 +𝑃 𝑃 )𝜓 12;12+34 56 57 58 12;12+34 12;12+34 67 68 57 58 57 68 12,12+34 1234 1234 123+45 we end the formalism with some Yakubovsky independent components as 𝜓 , 𝜓 , 𝜓 , 12;123 12;12+34 12,123 125+34 12+345 12+34+56 (𝜓 + 𝜓 ) and 𝜓 coupled in the equations. According to the 5-body sub-cluster distributions of 12,12+34 12,12+34 12,12+34 above independent components, it is well-known that two first components are related to 4𝑁 + 𝑁 + 𝑁 + 𝑁 + 𝑁 that refer to approximately effective 5-body model (𝛼 + 𝑛 + 𝑛 + 𝑛 + 𝑛 ). Respectively, third component refers to 3𝑁 + 2𝑁 + 𝑁 + 𝑁 + 𝑁 , which is not present in the effective 5-body 𝛼 -core model. Fourth component is a linear combinations that refer again to the third component. This is also beyond the effective 5-body 𝛼 -core model. The last one is related to 2𝑁 + 2𝑁 + 2𝑁 + 2𝑁 sub-cluster which is not contained in the effective 5-body 𝛼 -core model, either. 1234 1234 Regarding the sub-cluster underlying the Yakubovsky components, only 𝜓 and 𝜓 , are related to 12;123 12;12+34 approximately effective configuration of 𝛼 + 𝑛 + 𝑛 + 𝑛 + 𝑛 , where the 𝛼 -core and four loosely bound neutrons approximation are valid. Now, here we convinced enough why we choose two specific components, and those components are related to the effective 𝛼 -core structure. By considering the distinct Jacobi configurations related to 1234 1234 remaining Yakubovsky components, Fig.1, we can choose the first two components, 𝜓 and 𝜓 , and 12;123 12;12+34 obviously other components will not be taken into account. Hint: after removing the contribution interactions of the seventh and eighth nucleons in Eq. (1.34), i.e. removing their permutation operators, the 8-nucleon obtained equations leads to the 6-nucleon ones for halo nucleus He (See Eq. (67) in Ref. [12]). Such reduction, confirms that our implementing the Yakubovsky formalism for the 8- nucleon bound system for He, is appropriate approximation to describe the effective 𝛼 -core structure for Helium halo-bound nuclei, similar to the approximation of the six-nucleon for He in the effective 3-body model [12, 13]. In the above explanations we clearly have discussed why only the two specific components, 𝜓 and 12;123 𝜓 , are relevant to the effective 𝛼 -core structure of He, and the other Yakubovsky components in Eq. (1.34), 12;12+34 will not be taken into account. Therefore, Eq. (1.34) leads to 1234 1234 123 123 𝜓 (1 − 𝑃 −𝑃 −𝑃 −𝑃 )𝜓 𝒯 (−𝑃 ) 𝒯 12;123 45 46 47 48 12,123 ( ) = 𝐺 ( )[( )]. (1.35) 1234 12+34 1234 𝜓 𝒯 (1 − 𝑃 ) 0 𝜓 12;12+34 34 12;12+34 It is worth mentioning that switching off some Yakubovsky independent components changes the Hamiltonian of the comprehensive 8-nucleon system, so the Hamiltonian does not contain all interactions, but we show that in the remaining components the effective interactions of 𝛼 -core as the stable 4-nucleon sub-system is governed in the remaining components (Refer to Fig.1). It is obvious that after removing the contribution interactions of the 5 to 8’s nucleons in the 8-nucleon Yakubovsky equations, Eq. (1.34) straightforwardly reduces to the four- nucleon Yakubovsky equations for He ( See Ref. [10]) as 1234 1234 123 123 𝜓 𝒯 (−𝑃 ) 𝒯 𝜓 12;123 12;123 ( ) = 𝐺 ( )( ), (1.36) 1234 12+34 1234 𝜓 𝒯 (1 − 𝑃 ) 0 𝜓 12;12+34 12;12+34 in order to reduce the elaboration of the 8-nucleon coupled equations, Eq. (1.35), we can present the Yakubovsky 123 12+34 equations for the 8-nucleon problem as a function of two-nucleon 𝑡 -matrices and avoid using 𝒯 and 𝒯 . Consequently, the formalism can be simplified and the numerical solution of integral would be much faster in 123 12+34 comparison with the case where was need to solve first sub-cluster Faddeev-like equation to obtain 𝒯 and 𝒯 by using Pade′ approximation [10]. In addition to Eq. (1.36), the 4-nucleon Yakubovsky coupled equations has been obtained in [10] as 𝜓 𝐺 𝑡 𝑃 (1 − 𝑃 ) 𝐺 𝑡 𝑃 𝜓 1 0 12 34 0 12 1 (1.37) ( ) = ( )( ). ̃ ̃ 𝜓 𝐺 𝑡 𝑃 𝜓 𝐺 𝑡 𝑃 (1 − 𝑃 ) 2 0 12 2 0 12 34 also in Eq. (1.36) the 4-nucleon transition operators obey as 123 (1.38) = 𝑡 𝑃 + 𝑃 𝐺 𝑡 ≡ 𝒯 , 12 0 12 12+34 ̃ ̃ ̃ ̃ ̃ ̃ (1.39) 𝑇 𝑃 = 𝑡 𝑃 + 𝑃 𝐺 𝑡 𝑇 𝑃 ≡ 𝒯 , 12 0 12 According to [10] the 4-nucleon matrix form transition in Eq. (1.36) and (1.37) are similar in action, and according to Eqs. (1.35) and (1.36), the 8-nucleon and the 4-nucleon matrix form transitions are equal in action. At this point, the matrices of Eq. (1.37) can be identified and Eq. (1.35) can be written in terms of two-nucleon t−matrices operator since they are given as 1234 1234 𝜓 (1 − 𝑃 −𝑃 −𝑃 −𝑃 )𝜓 𝐺 𝑡 𝑃 (1 − 𝑃 ) 𝐺 𝑡 𝑃 12;123 45 46 47 48 12,123 0 12 34 0 12 ( ) = ( )[( )], 1234 1234 ̃ ̃ 𝐺 𝑡 𝑃 (1 − 𝑃 ) 𝐺 𝑡 𝑃 𝜓 𝜓 0 12 (1.40) 12;12+34 0 12 34 12;12+34 finally, the linear form of the obtained 8-nucleon Yakubovsky equations in the case of effective 𝛼 -core structure can be written as 1234 1234 1234 (1.41) 𝜓 = 𝐺 𝑡 𝑃 [(1 − 𝑃 )(1 − 𝑃 −𝑃 −𝑃 −𝑃 )𝜓 + 𝜓 ], 12;123 0 12 34 45 46 47 48 12;123 12;12+34 1234 1234 1234 (1.42) 𝜓 = 𝐺 𝑡 𝑃 [(1 − 𝑃 )(1 − 𝑃 −𝑃 −𝑃 −𝑃 )𝜓 + 𝜓 ]. 12;12+34 0 12 34 45 46 47 48 12;123 12;12+34 In the next step we will describe characteristic numerical techniques such as introducing Jacobi configuration of each independent component in momentum space representation to provide corresponding basis states and to understand the halo structure model of the 8-nucleon system with respect to the regarded relevant remaining components; integral representation of the coupled equations; eigenvalue equation form and iteration method for solution of it, and discussions about the halo structure configurations that describe the 8-nucleon halo-bound system in the case of effective 𝛼 -core structure. II. Technicalities for Numerical Applications In this section, first, Jacobi configurations in momentum space are represented to define corresponding momentum basis states for evaluations of the integral kernels. Next, details for a typical algorithm are described to solve coupled integral equations based on PW representation. In order to describe the corresponding momentum basis 8 1234 1234 states, for two independent components of the dominant 𝛼 -core structure of He, namely 𝜓 and 𝜓 the 12;123 12;12+34 standard Jacobi momenta are described in Fig. 1, respectively. 𝑇𝑃 𝑇𝑃 (1) (2) Fig. 1: Schematic representation of the two independent components of 8-nucleon bound system for He (four loosely bound neutrons with an inert 𝛼 -core) in Jacobi configurations. In the bag the 𝛼 -core tightly interacts in a 4-nucleon subsystem. For the first component 𝜓 in terms of the first configuration in Fig. 1 sets corresponding to 𝒖 -set chains 12;123 1/2 −1/2 0 0 0 0 0 0 −1/3 −1/3 2/3 0 0 0 0 0 2 2 0 0 0 0 −1/4 −1/4 −1/4 3/4 𝒑 1/2 −1/2 𝒑 𝒖 0 0 0 0 0 0 4 = , 0 0 5 5 −1/3 −1/3 −1/3 −1/3 2/3 2/3 1/2 −1/2 0 0 0 0 0 0 𝒖 𝒑 7 7 (2.1) −1/4 −1/4 −1/4 −1/4 −1/4 −1/4 3/4 3/4 ( ) (𝒑 ) ( ) 1 1 1 1 1 1 1 1 in the non-relativistic case, we can define the kinetic energy operator by two equivalent forms. The inverse form of the above-mentioned transfer matrix is used to represent the kinetic energy in terms of 𝒖 -set Jacobi momenta. 2 2 2 2 2 2 2 2 𝑝 𝑢 3𝑢 2𝑢 𝑢 3𝑢 𝑢 1𝑢 𝑖 1 2 3 4 5 6 7 𝐻 = ∑ ≡ + + + + + + , (2.2) 2𝑚 𝑚 4 𝑚 3 𝑚 𝑚 8 𝑚 𝑚 3 𝑚 𝑖 =1 where 𝒑 is an individual particle momentum in the center of mass situation (under the condition ∑𝒑 = 0) that is 𝑖 𝑖 described by relative Jacobi momenta 𝒖 ;( 𝑖 = 1,2,3,4,5,6,7,8). Similarly, to 𝜓 in terms of the second 𝑖 12;12+34 configuration in Fig. 1 belongs 0 0 0 0 0 0 1/2 −1/2 1/2 −1/2 0 0 0 0 0 0 𝒑 0 0 0 0 −1/2 −1/2 1/2 1/2 𝒑 1/2 −1/2 𝒑 𝒗 0 0 0 0 0 0 4 = , (2.3) 0 0 5 5 −1/3 −1/3 −1/3 −1/3 2/3 2/3 1/2 −1/2 0 0 0 0 0 0 𝒗 𝒑 7 7 −1/4 −1/4 −1/4 −1/4 −1/4 −1/4 3/4 3/4 ( ) (𝒑 ) ( ) 1 1 1 1 1 1 1 1 according to some appearing Jacobi vectors in the two different configurations in Fig.1, correspondingly, some equal Jacobi momenta have been represented by different notations, i.e. 𝒗 = 𝒖 , because of the standard representation of 7 7 the above transfer matrices, for providing inverse transformations. The kinetic energy form in terms of 𝒗 -set Jacobi momenta, are given as 2 2 2 2 2 2 2 2 𝑝 𝑣 𝑣 1𝑣 𝑣 3𝑣 𝑣 1𝑣 𝑖 1 2 3 4 4 4 4 𝐻 = ∑ ≡ + + + + + + . 2𝑚 𝑚 𝑚 2 𝑚 𝑚 8 𝑚 𝑚 3 𝑚 (2.4) 𝑖 =1 Then, we introduce the two basis-states corresponding two independent components in PW representation. Obviously, these basis are separated into the 3 spaces that are 1-momentum/ 2-angular momentum and spin/ 3-isospin. 1234 1234 The relevant basis-states in PW representation for 𝜓 and 𝜓 respectively are given as 12;123 12;12+34 |𝑢 ; 𝛾 ⟩ ≡ |𝑢 𝑢 𝑢 𝑢 𝑢 𝑢 𝑢 ⟩ 𝑢 1 2 3 4 5 6 7 1 1 | ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ⨂ (𝑙 𝑠 𝑗 (𝑙 )𝑗 𝑗 𝑗 𝐼 (𝑙 )𝑗 𝐼 𝑗 𝐼 𝑙 𝑠 𝑗 𝑙 𝑗 𝐼 𝐼 𝐼 𝐼 𝑙 𝑠 𝑗 𝑙 𝑗 𝐼 𝐼 𝐼 𝐽 ,𝑀 ⟩ 1 12 1 2 2 1 2 3 3 4 3 4 4 5 56 5 4 5 5 4 5 6 6 12 6 7 6 7 6 7 𝑗 2 2 1 1 (2.5) ⨂|(𝑡 )𝑡 (𝑡 )𝑡 (𝑡 𝑡 )𝑡 (𝑡 𝑡 ) 𝑇 ,𝑀 ⟩, 12 3 3 4 4 56 6 6 78 𝑡 2 2 ⟩ ⟩ |𝑣 ; 𝛾 ≡ |𝑣 𝑣 𝑣 𝑣 𝑣 𝑣 𝑣 𝑣 1 2 3 4 5 6 7 ⨂|(𝑙 𝑠 )𝑗 (𝑙 𝑠 )𝑗 (𝑗 𝑗 )𝑆 (𝐿 𝑆 )𝐼 (𝑙 𝑠 )𝑗 (𝑙 𝑗 )𝐼 (𝐼 𝐼 )𝐼 (𝑙 𝑠 )𝑗 (𝑙 𝑗 )𝐼 (𝐼 𝐼 ) 𝐽 ,𝑀 ⟩ 1 12 1 2 34 2 1 2 4 56 5 5 5 5 5 6 4 78 6 7 6 7 6 7 𝑗 (2.6) | ) ( ) ⨂ (𝑡 𝑡 𝑡 (𝑡 𝑡 )𝑡 𝑡 𝑡 𝑇 ,𝑀 ⟩, 12 34 1_4 1_4 56 6 6 78 𝑡 here the first part in the above-introduced basis states are magnitude of momentum vectors in PW representation. In the second part the orbital angular momentum 𝑙 go with the 𝒖 \ 𝒗 , 𝑠 which are two-body spins for nucleons , 𝑗 𝑖 𝑖 𝑖 𝑖 𝑗 𝑖 are total 1- and 2-body angular momenta coupled out of orbital and spin angular momenta, 𝐼 are total 𝑖 -body angular momenta except 𝐼 and 𝐼 . These are total angular momenta of 5, 6 and 7, 8 respectively. Evidently, the third part 5 7 refers to isospin. In this step, in order to implement the numerical techniques, the two coupled equations, Eq. (1.41) and Eq. (1.42), by inserting the completeness relations between the permutation operators, can be projected in momentum space as 1234 ′2 ′ ′′2 ′′ ⟨𝑢 ; 𝛾 |𝜓 ⟩ = ∑∫𝑢 𝑑 𝑢 ∑ ∫𝑢 𝑑 𝑢 12;123 𝛾 𝛾 ′′ 𝑢 ′ ′′ ′′ 1234 ⟨𝑢 ; 𝛾 |𝐺 𝑡 𝑃 |𝑢 ′; 𝛾 ⟩⟨𝑢 ′; 𝛾 |(1 − 𝑃 )(1 − 𝑃 −𝑃 −𝑃 −𝑃 )| 𝑢 ; 𝛾 ′′⟩⟨ 𝑢 ; 𝛾 ′′|𝜓 ⟩ 𝑢 0 12 𝑢 ′ 𝑢 ′ 34 45 46 47 48 12;123 𝑢 𝑢 ′2 ′ ′2 ′ ′ ′ 1234 + ∑∫𝑢 𝑑 𝑢 ∑∫𝑣 𝑑 𝑣 ⟨𝑢 ; 𝛾 |𝐺 𝑡 𝑃 |𝑢 ′; 𝛾 ⟩⟨𝑢 ′; 𝛾 |𝑣 ; 𝛾 ′⟩⟨𝑣 ; 𝛾 ′|𝜓 ⟩, 𝑢 0 12 𝑢 ′ 𝑢 ′ 12;12+34 𝑣 𝑣 𝛾 𝛾 𝑢 ′ 𝑣 ′ (2.7) 1234 ′2 ′ ′2 ′ ⟨𝑣 ; 𝛾 | 𝜓 ⟩ = ∑∫𝑣 𝑑 𝑣 ∑∫𝑢 𝑑 𝑢 𝑣 12;12+34 𝛾 𝛾 𝑣 ′ 𝑢 ′ ′ ′ 1234 ⟨𝑣 ; 𝛾 |𝐺 𝑡 𝑃 |𝑣 ; 𝛾 ′⟩⟨𝑣 ; 𝛾 ′|(1 − 𝑃 )(1 − 𝑃 −𝑃 −𝑃 −𝑃 )|𝑢 ′; 𝛾 ⟩⟨𝑢 ′; 𝛾 |𝜓 ⟩ 𝑣 0 12 34 45 46 47 48 𝑢 ′ 𝑢 ′ 12;123 𝑣 𝑣 (2.8) ′2 ′ ′ ′ 1234 + ∑∫𝑣 𝑑 𝑣 ⟨𝑣 ; 𝛾 |𝐺 𝑡 𝑃 |𝑣 ; 𝛾 ′⟩⟨𝑣 ; 𝛾 ′|𝜓 ⟩, 𝑣 0 12 𝑣 𝑣 12;12+34 𝑣 ′ using techniques like ones presented in Refs. [10-14], it is straightforward to evaluate in PW analysis like the above- mentioned kernels and shifted momenta. It is well-known that the variables for 7 amplitudes on the right side are in general linear combinations of intermediate integration variables which include angles besides momentum magnitudes. After evaluating each term in the above-mentioned coupled integral equations in the standard PW analysis, the obtained equations are the starting point for numerical calculations in the eigenvalue equation form. We can replace the continuous variables in the numerical treatment by a dependence on certain distinct values by using the Gauss-Legendre discretization. The eigenvalue equation can be solved by the iteration method. It can apply Lanczos-like scheme that is efficient for nuclear few-body problems [10]. Also the evaluated coupled integral 𝑖𝑗 equations involve a very large number of interpolations, it can be used the cubic Hermit splines of Ref. [17] for its accuracy and high computational speed. The eigenvalue equation form can be given as 1234 1234 1234 1234 (2.9) 𝜂 (𝐸 ) 𝜓 ( 𝜓 , 𝜓 )= 𝑘 (𝐸 ) 𝜓 ( 𝜓 , 𝜓 ), 12;123 12;12+34 12;123 12;12+34 ( ) where 𝐸 is the energy eigenvalue at which auxiliary Yakubovsky kernel eigenvalue is 𝜂 𝐸 = 1. Because the energy is varied such that one reaches eigenvalue to be one. The Yakubovsky kernel of the equations 𝑘 (𝐸 ) is just depend to 1234 1234 energy E, and 𝜂 (𝐸 ) is eigenvalue with 𝜓 ( 𝜓 , 𝜓 ) as the corresponding eigenvector. In order to solve the 12;123 12;12+34 eigenvalue equation, Eq. (2.9), it can be used the Gaussian quadrature grid points. The coupled equations represent a set of homogenous integral equations, which turn into a high-dimension matrix eigenvalue equation after discretization. Starting from an arbitrary initial 𝜓 ≡ 𝜓 one generates by consecutive applications of 𝑘 (𝐸 ) a sequence of amplitudes 𝜓 , which after orthogonalization form a basis into which 𝜓 is expanded. Then the energy is varied ( ) such that one reaches 𝜂 𝐸 = 1. In order to investigate the loosely interactions of the halo bound neutrons in He, −1 calculations of the momentum probability density as 𝑛 (𝑢 )[𝑓𝑚 ] in terms of 𝑢 [𝑓𝑚 ], correspondingly (𝑣 ) can be 𝑖 𝑖 𝑖 done. Why we just select 4 and 6 Jacobi momenta to calculate the momentum probability density? Because they are the only ones that make up contributions of the halo neutrons (refer to Eqs. (2.1) and (2.3)). By calculating momentum probability densities, we can represent the situation of the halo bound neutrons in momentum space. It is expected that the representation of momentum probability densities approximately are wide and short in comparison with tightly bound neutrons, namely neutrons inside the 𝛼 -core subsystem. Such those momentum distributions confirms that the 8-nucleon in the case of effective 𝛼 -core model is a halo structure system, exactly like halo nucleus He. III. Summary The applications of Faddeev and Yakubovsky equations are traditionally for 3- and 4-nucleon bound and scattering states. Today, after the experiences with 4-, 5- and 6-nucleon bound-state problems within the Yakubovsky approach [13, 14] that the technical expertise has developed in the recent years and the very strong increase of computer power just recently achieved, it is possible to approach the 8-nucleon bound-state problem in that formalism and then numerical solution. For a full solution of the 8-nucleon Yakubovsky equations for He, modern super- computers with parallel algorithms is required. But in order to approach the calculations of the 8-nucleon bound system, for the first time, we have implemented the 8-nucleon Yakubovsky equations for the halo nucleus He in the case of simple effective 𝛼 -core structure. To this end, we stop the sequential sub-clustering with 5 fragments, namely 𝛼 + 𝑛 + 𝑛 + 𝑛 + 𝑛 , though the additional step with 4 fragments could be easily performed, but is not applicable for a physical system. In this method, we restricted the formulation to two-nucleon forces only because loosely bound neutrons with the center of an inert 𝛼 -core tend to have much stronger pairing interactions for large radii neutrons [5]. Therefore, as a simplification in the formalism, and above fact in such a model, three-nucleon forces was ignored. As a result, the implementation of the 8-nucleon bound system in the halo structure model within the Yakubovsky approach leads to two coupled equations in terms of two independent Yakubovsky components. Subsequently, in order to approach the halo structures of He, in addition to the analytical implementations, we have introduced the Jacobi configurations of those components in momentum space representation. Eventually, we provide technicalities which were considered useful for a numerical performance, such as bound-state calculations and proving the halo structure of bound neutrons that are outside of the 𝛼 -core. References [1] M. Brodeur, et. al. Phys. Rev. Lett. 108, 052504 –31 Jan (2012). [2] M.V.Zhukov, et. al. Physics Reports, Volume 231, Issue 4, August (1993), Pages 151-199. [3] S. Bacca, A. Schwenk, G. Hagen, et. al. Eur. Phys. J. A, 42: 553 (2009). [4] L. B. Wang et al., Phys. Rev. Lett. 93, 142501 (2004). [5] P. Mueller et al., Phys. Rev. Lett. 99, 252501 (2007). [6] V. L. Ryjkov et al., Phys. Rev. Lett. 101, 012501 (2008). [7] S. C. Pieper and R. B. Wiringa, Annu. Rev. Nucl. Part. Sci. 51, 53 (2001); S. C. Pieper, arXiv:0711.1500. [8] P. Navratil and W. E. Ormand, Phys. Rev. C 68, 034305 (2003). [9] P. Navratil, V. G. Gueorguiev, J. P. Vary, W. E. Ormand and A. Nogga, Phys. Rev. Lett. 99, 042501 (2007). [10] 10. H. Kamada and W. Glӧckle, Nucl. Phys. A 548, 205 (1992). [11] A. Nogga, H. Kamada and W. Glöckle, Phys. Rev. Lett. 85, 944 (2000). [12] W. Glӧckle and H. Witala, Few-Body Syst. 51, 27-44 (2011). [13] E. Ahmadi Pouya and A. A. Rajabi, Acta, Phys, Pol, B 48: 1279 (2017). [14] E. Ahmadi Pouya and A. A. Rajabi, Eur. Phys. J. Plus, 131: 240 (2016). [15] A. C. Fonseca, Phys. Rev. C 30, 35 (1984). [16] W. Glöckle: The Quantum Mechanical Few-Body Problem. Springer-Verlag, New York (1983). [17] D. Huber, H. Witala, A. Nogga, W. Glӧckle and H. Kamada, Few-Body Syst. 22, 107 (1997) . http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Nuclear Theory arXiv (Cornell University)

Implementation of the 8-Nucleon Yakubovsky Formalism for Halo Nucleus 8He

Pouya, Eskandar Ahmadi; Rajabi, Ali Akbar
Nuclear Theory , Volume 2019 (1904) – Apr 29, 2019
Download PDF
13 pages

Loading next page...
 
/lp/arxiv-cornell-university/implementation-of-the-8-nucleon-yakubovsky-formalism-for-halo-nucleus-wNAPTHjQFS
ISSN
2326-8867
eISSN
ARCH-3340
DOI
10.11648/j.ajmp.20190803.12
Publisher site
See Article on Publisher Site

Abstract

In order to study the bound-state structure of the Helium halo nuclei, the 8-nucleon Yakubovsky formalism has been implemented for He in a 5-body sub-cluster model, i.e. 𝛼 + 𝑛 + 𝑛 + 𝑛 + 𝑛 . In this case, the 8- nucleon Yakubovsky equations has been obtained in the form of two coupled equations, based on the two independent components. In addition, by removing the contribution interactions of the 8 and 7’s bound nucleons in the formalism, the obtained equations explicitly reduce to the 6-nucleon Yakubovsky equations for He, in the case of effective 3-body model, i.e. 𝛼 + 𝑛 + 𝑛 . In view of the expectation for the dominant structure of He, namely an inert 𝛼 -core and four loosely-bound neutrons, Jacobi configurations of the two components in momentum space have been represented to provide technicalities which were considered useful for a numerical performance, such as bound-state calculations and momentum density distributions for halo-bound neutrons. PACS. No: 21.45.+v Few-body systems; 21.30.Fe Forces in hadronic systems and effective interactions; 21.30.-x Two-nucleon interactions; 21.10.Dr Binding energies and masses. Introduction Studies of rich-neutron light nuclei present stringent tests for their calculations as they probe aspects of the interactions that are less prevalent in nuclei closer to stability. Among them, the physics of strong interactions gives rise to new structures in rich-neutron light nuclei. One of the best prominent cases are the helium halo nuclei, namely 6 8 4 He and He with two and four loosely-bound neutrons, forming an extended halo around the He core. In this regard, 6 8 the rich-neutron helium isotopes, He and He, are the topic of this study. Indeed, both nuclei are Borromean halos [1, 2], they have no bound excited states, and they exhibit the binding energy. Both helium halo nuclei are radioactive and decay quickly that complicates efforts to measure their properties. He has a time-life of only a tenth of a second and in this little time, can be considered as a stable nucleus. Practically, detection of the rare He is a major step for particle physicists. But, it is very important and interesting to understand how helium configures itself with strong interactions after being produced from a particle accelerator. In recent decades, combination of nuclear and atomic physics methods has allowed new courses of precision 6 8 measurements of the ground-state energies and charge radii of He [3, 4] and He [5, 6]. Dealing with the theoretical study of these nuclei is very challenging, because one needs to describe the small separation energy of the halo neutrons and the large radius of the whole system, in addition to the challenges of binding energy calculations. Because they are light-mass nuclei, one can use traditional methods to study them, such as Green’s Function Monte Carlo (GFMC) method [7] and the No-Core Shell Model (NCSM) [8]. In addition, some of the standard methods in nuclear theory has been the progress of Effective Field Theory (EFT) and the Renormalization Group (RG) to nuclear forces. While such helium halo nuclei have been investigated using the NCSM [9], but still, no sense has been obtained for the halo structure properties in such theoretical studies. 8 6 In this paper, to approach a simple physical structure of the observed effects, namely He and He halo-bound nuclei, and to understand the physics that governs the neutron distributions in helium halos, it was felt that by using current/modern computers with high computational speed, beyond previous works [10-14], a further step to study the light halo nuclei, namely He as an 8-nucleon bound system problem, within the reliable Yakubovsky method [15] is now desirable. It is worth mentioning that our obtained results are more difficult to implement bound-state problem E.Ahmady.ph@ut.ac.ir A.A.Rajabi@shahroodut.ac.ir in few-body branch, due to its halo nature and high-dimensional calculations. However, we provide technicalities for numerical performances. This paper is written as follows: In sect. I, the Yakubovsky formalism to the 8-nucleon bound system problem using the standard notation [16] is explicitly derived. In this case, the Lippmann-Schwinger scheme twice step-by- step has been applied on the integral form of Schrödinger equation, and then the identity of the nucleons are implemented. Finally, the 8-nucleon Yakubovsky formalism for effective 𝛼 -core structure in a 5-body sub-cluster model of He leads to a set of two coupled equations related to 2 different components. In sect. II, useful technicalities for numerical applications have been provided. To this aim, corresponding Jacobi configurations of each Yakubovsky component is schematically represented and then the integral form of coupled equations are represented by introducing basis states based on Jacobi momenta in partial-wave (PW) analysis. Finally, a summary in sect. III is given. I. The 8-nucleon Yakubovsky Formalism In branch of few-body systems, each N-body bound system has 𝑁 (𝑁 − 1)/2 different two-body interactions. Therefore, in the 8-nucleon bound system there are 28 different two-nucleon interactions or 28 different cluster decompositions (𝑎 ) having 7 clusters. They are indexed by the only pairing forces in cluster 𝑎 they contain, e.g. 𝑎 = 7 7 7 12 ≡ 12+ 3 + 4 + 5+ 6+ 7 + 8. In order to study the 8-nucleon bound system in the framework of the Yakubovsky method using the sub-cluster notation [16], the method is to first sum up the pair forces in each 7-body fragment (𝑎 ), in a second step among all 6-body fragments (𝑎 ), and then in a third step among all 5-body fragments (𝑎 ). We cut 6 5 off that formalism ending with 5-body sub-clusters in the spirit of the usually used approximating effective 5-body model that is 𝛼 + 𝑛 + 𝑛 + 𝑛 + 𝑛 . In the standard formalism, for an individual 𝛼 -particle, such two-body interactions, namely 𝛼 − 𝑛 and 𝛼 -clusters have not been used, because in this framework 𝛼 -core interacts as a four-nucleon sub- system in the 8-nucleon system. We start with the non-relativistic Schrödinger equation for the 8-nucleon bound system as (1.1) (𝐻 + ∑𝑉 )𝛹 = 𝐸 𝛹 , 0 𝑎 where 𝐻 refers to the free Hamiltonian operator of the 8-nucleon system that will be introduced in the next section, and ∑ 𝑉 ≡ 𝑉 + ⋯+ 𝑉 is the summation of all paring interactions with 28 terms. According to the Lipmann- 𝑎 𝑎 12 78 7 7 Schwinger scheme, Eq. (1.1) is rewritten into an integral equation (1.2) 𝛹 = 𝐺 ∑𝑉 𝛹 , 0 𝑎 where 𝐺 stand for the 8-nucleon free Green’s function operator and in the case of bound states formalism we have −1 𝐺 = [𝐸 − 𝐻 ] . Like previous works for five- and six-nucleon Yakubovsky formalism [12, 13], we start from the 0 0 integral form of the Schrödinger equation, Eq. (1.2), applying Lipmann-Schwinger scheme twice step-by-step. We can define 𝜓 ≡ 𝐺 𝑉 𝛹 by using Eqs. (1.2) and (1.7). Therefore, Faddeev like equation for 8-nucleon Yakubovsky 𝑎 0 𝑎 7 7 amplitudes are given as (1.3) 𝜓 ≡ 𝐺 𝑡 ∑𝛿 𝜓 , 𝑎 0 𝑎 𝑎 𝑏 𝑏 7 7 7 7 7 where 𝑡 is a pairing 𝑡 -matrix operator that follows the Lippmann-Schwinger equation as 𝑡 = 𝑉 + 𝑉 𝐺 𝑡 and 𝑎 𝑎 𝑎 𝑎 0 𝑎 7 7 7 7 7 𝛿 = 1 − 𝛿 . For the next step sub-cluster form the 6-body fragments amplitudes can be written 𝑎 𝑏 𝑎 𝑏 7 7 7 7 𝑎 𝑎 6 𝑎 6 6 ̅ (1.4) 𝜓 = 𝐺 ∑ 𝒯 (𝜓 ) =𝐺 ∑ 𝒯 ∑ 𝛿 𝜓 , 𝑎 𝑎 0 𝑎 0 𝑎 𝑏 𝑏 𝑏 7 6 𝑎 𝑏 7 𝑎 𝑏 6 6 7 6 7 7 7 7 𝑏 ⊂𝑎 𝑏 ⊂𝑎 𝑏 ⊂𝑏 7 6 7 6 7 6 where 𝒯 obeys as 𝑎 𝑏 7 7 (1.5) 𝑎 𝑎 6 6 ̅ ̅ 𝒯 = 𝑡 𝛿 + 𝐺 ∑ 𝑡 𝛿 𝒯 , 𝑎 𝑏 𝑎 𝑎 𝑏 0 𝑎 𝑎 𝑐 𝑐 𝑏 7 7 7 7 7 7 7 7 7 7 𝑐 ⊂𝑎 7 6 next, more decompose the right side of Eq. (1.4) according to 5-body fragments are given as 𝑎 𝑎 5 6 (1.6) 𝜓 = 𝐺 ∑ 𝒯 ∑ 𝛿 𝜓 . 0 𝑎 𝑏 𝑏 𝑏 𝑎 ,𝑎 𝑎 𝑏 7 6 6 6 7 6 7 7 𝑏 ⊂𝑏 𝑏 ⊂𝑎 7 6 7 6 𝑏 ⊂𝑎 6 5 briefly, we work out that method ending with (𝑎 ) in the spirit of the applied though approximate effective 𝛼 + 𝑛 + 𝑛 + 𝑛 + 𝑛 in a 5-body model. Further details could be found in Refs. [12, 13]. In the following, the single nucleons in sub-clusters, namely 7-, 6- and 5-body fragments, i.e. (𝑎 ), (𝑎 ) and (𝑎 ) respectively, will no longer be indexed. 7 6 5 Next, in order to get to 8-nucleon Yakubovsky coupled equations, the identity of the nucleons has been applied. In addition, the number of huge total WF Yakubovsky components has been obtained. It is well-known, the 8-nucleon total WF consists of the summation of all pairing Yakubovsky components with 28 terms which are as 𝛹 = ∑𝜓 ≡ 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 𝑎 12 13 14 15 16 17 18 +𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 23 24 25 26 27 28 +𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 34 35 36 37 38 +𝜓 + 𝜓 + 𝜓 + 𝜓 45 46 47 48 (1.7) +𝜓 + 𝜓 + 𝜓 56 57 58 +𝜓 + 𝜓 67 68 +𝜓 , where the first pair (𝜓 ), having 7-body fragments, consists of 21 components of 6-body fragments (𝑎 ) as 12 6 𝜓 = ∑ 𝜓 𝑎 ≡12 𝑎 ≡12,𝑎 7 7 6 (1.8) 𝑎 ≡12⊂𝑎 7 6 ≡ 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 12,123 12,124 12,125 12,126 12,127 12,128 12,12+34 12,12+35 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 12,12+36 12,12+37 12,12+38 12,12+45 12,12+46 12,12+47 12,12+48 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 , 12,12+56 12,12+57 12,12+58 12,12+67 12,12+68 12,12+78 where (𝑎 ) refers to any 6-body fragments containing the pair restricted to (𝑎 ≡12) and the sum runs over pairs 6 7 𝑎 ⊂ 𝑎 . It means that the sub-clusters (𝑎 ), when broken up, lead to the sub-clusters (𝑎 ). Next first term as 𝜓 7 6 6 7 12,123 is a sub-cluster of 15 components of 5-body fragments (𝑎 ) 𝜓 = ∑ 𝜓 𝑎 ≡12,𝑎 ≡123 7 6 𝑎 ≡12,𝑎 ≡123 7 6 (1.9) 𝑎 ≡123⊂𝑎 6 5 1234 1235 1236 1237 1238 123+45 123+46 123+47 ≡ 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 12,123 12,123 12,123 12,123 12,123 12,123 12,123 12,123 123+48 123+56 123+57 123+58 123+67 123+68 123+78 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 , 12,123 12,123 12,123 12,123 12,123 12,123 12,123 the summation of all 5-body fragments ends up to the 8-nucleon Yakubovsky components for total WF as 𝑎 𝑎 5 5 (1.10) 𝛹 = ∑∑∑𝜓 = 8,820 𝜓 , 𝑎 ,𝑎 𝑎 ,𝑎 7 6 7 6 𝑎 𝑎 𝑎 7 6 5 We find that there are 8,820 components for the total WF of constituent 8 nucleons in a 5-body model. It is worthwhile to mention that to calculate the expectation value of the Hamiltonian operator, the above-mentioned total WF’s have to be used. Clearly, the total WF is anti-symmetrized following the Pauli principle. In the next step, the identity of the nucleons in sub-cluster components are implemented. We start from Eq. (1.4) choosing the case 7-body fragments 𝑎 ≡ 12 with 6-body fragments 𝑎 ≡ 123 7 6 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 12,124 12,125 12,126 12,127 12,128 +𝜓 +𝜓 +𝜓 +𝜓 +𝜓 12,12+34 12,12+35 12,12+36 12,12+37 12,12+38 +𝜓 +𝜓 +𝜓 +𝜓 12,12+45 12,12+46 12,12+47 12,12+48 𝜓 = 𝐺 𝒯 12,123 0 12,12 +𝜓 +𝜓 +𝜓 12,12+56 12,12+57 12,12+58 +𝜓 +𝜓 12,12+67 12,12+68 +𝜓 ( ) 12,12+78 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 23,234 23,235 23,236 23,237 23,238 +𝜓 +𝜓 +𝜓 +𝜓 +𝜓 23,23+14 23,23+15 23,23+16 23,23+17 23,23+18 +𝜓 +𝜓 +𝜓 +𝜓 23,23+45 23,23+46 23,23+47 23,23+48 +𝐺 𝒯 0 12,23 +𝜓 +𝜓 +𝜓 23,23+56 23,23+57 23,23+58 +𝜓 +𝜓 23,23+67 23,23+68 +𝜓 ( ) 23,23+78 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 31,314 31,315 31,316 31,317 31,318 +𝜓 +𝜓 +𝜓 +𝜓 +𝜓 31,31+24 31,31+25 31,31+26 31,31+27 31,31+28 +𝜓 +𝜓 +𝜓 +𝜓 31,31+45 31,31+46 31,31+47 31,31+48 +𝐺 𝒯 , 12,31 +𝜓 +𝜓 +𝜓 31,31+56 31,31+57 31,31+58 +𝜓 +𝜓 31,31+67 31,31+68 +𝜓 ( ) 31,31+78 (1.11) it is easily seen that, by using adequate permutation operators, the huge equation above can be rewritten as 123 123 123 𝜓 = 𝐺 (𝒯 + 𝒯 𝑃 𝑃 + 𝒯 𝑃 𝑃 ) 12,123 0 12,12 12,23 12 23 12,31 13 23 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 12,124 12,125 12,126 12,127 12,128 +𝜓 +𝜓 +𝜓 +𝜓 +𝜓 12,12+34 12,12+35 12,12+36 12,12+37 12,12+38 +𝜓 +𝜓 +𝜓 +𝜓 12,12+45 12,12+46 12,12+47 12,12+48 × , +𝜓 +𝜓 +𝜓 (1.12) 12,12+56 12,12+57 12,12+58 +𝜓 +𝜓 12,12+67 12,12+68 +𝜓 ( ) 12,12+78 by defining the bellow transition operator as 123 123 123 123 𝒯 = 𝒯 + 𝒯 𝑃 𝑃 + 𝒯 𝑃 𝑃 , (1.13) 12,12 12,23 12 23 12,31 13 23 Eq. (1.12) is simplified to 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 12,124 12,125 12,126 12,127 12,128 +𝜓 +𝜓 +𝜓 +𝜓 +𝜓 12,12+34 12,12+35 12,12+36 12,12+37 12,12+38 +𝜓 +𝜓 +𝜓 +𝜓 12,12+45 12,12+46 12,12+47 12,12+48 𝜓 = 𝐺 𝒯 . 12,123 0 +𝜓 +𝜓 +𝜓 12,12+56 12,12+57 12,12+58 +𝜓 +𝜓 (1.14) 12,12+67 12,12+68 +𝜓 ( ) 12,12+78 123 123 123 where 𝒯 follows Lipmann-Schwinger equation as 𝒯 = 𝑡 𝑃 + 𝑡 𝑃 𝐺 𝒯 and 𝑃 = 𝑃 𝑃 + 12 12 0 12 23 𝑃 𝑃 .Subsequently, starting again from Eq. (1.4) to select the case 7-body fragments 𝑎 ≡ 12 with 6-body fragments 13 23 7 𝑎 ≡ 12+ 34 𝜓 + 𝜓 +𝜓 +𝜓 +𝜓 +𝜓 12,123 12,124 12,125 12,126 12,127 12,128 +𝜓 +𝜓 +𝜓 +𝜓 12,12+35 12,12+36 12,12+37 12,12+38 +𝜓 +𝜓 +𝜓 +𝜓 12,12+45 12,12+46 12,12+47 12,12+48 12+34 𝜓 = 𝐺 𝒯 12,12+34 0 12,12 +𝜓 +𝜓 +𝜓 12,12+56 12,12+57 12,12+58 +𝜓 +𝜓 12,12+67 12,12+68 +𝜓 ( ) 12,12+78 𝜓 + 𝜓 +𝜓 +𝜓 +𝜓 +𝜓 34,134 34,234 34,345 34,346 34,347 34,348 +𝜓 +𝜓 +𝜓 +𝜓 34,34+15 34,34+16 34,34+17 34,34+18 +𝜓 +𝜓 +𝜓 +𝜓 34,34+25 34,34+26 34,34+27 34,34+28 12+34 + 𝐺 𝒯 , 0 12,34 +𝜓 +𝜓 +𝜓 34,34+56 34,34+57 34,34+58 +𝜓 +𝜓 34,34+67 34,34+68 (1.15) +𝜓 ( ) 34,34+78 subsequently, by using permutation properties, Eq. (1.15) can be rewritten as 𝜓 + 𝜓 +𝜓 +𝜓 +𝜓 +𝜓 12,123 12,124 12,125 12,126 12,127 12,128 +𝜓 +𝜓 +𝜓 +𝜓 12,12+35 12,12+36 12,12+37 12,12+38 +𝜓 +𝜓 +𝜓 +𝜓 12,12+45 12,12+46 12,12+47 12,12+48 12+34 12+34 𝜓 = 𝐺 (𝒯 + 𝒯 𝑃 𝑃 ) , 12,12+34 0 12,12 12,34 13 24 +𝜓 +𝜓 +𝜓 12,12+56 12,12+57 12,12+58 +𝜓 +𝜓 12,12+67 12,12+68 (1.16) +𝜓 ( ) 12,12+78 then defining 12+34 12+34 12+34 𝒯 = 𝒯 + 𝒯 𝑃 𝑃 , (1.17) 13 24 12,12 12,34 Eq. (1.16) is simplified to 𝜓 + 𝜓 +𝜓 +𝜓 +𝜓 +𝜓 12,123 12,124 12,125 12,126 12,127 12,128 +𝜓 +𝜓 +𝜓 +𝜓 12,12+35 12,12+36 12,12+37 12,12+38 +𝜓 +𝜓 +𝜓 +𝜓 12,12+45 12,12+46 12,12+47 12,12+48 12+34 𝜓 = 𝐺 𝒯 . 12,12+34 0 +𝜓 +𝜓 +𝜓 12,12+56 12,12+57 12,12+58 +𝜓 +𝜓 12,12+67 12,12+68 (1.18) +𝜓 ( ) 12,12+78 12+34 12+34 12+34 ̃ ̃ ̃ where 𝒯 follows the equation 𝒯 = 𝑡 𝑃 + 𝑡 𝑃 𝐺 𝒯 and where 𝑃 = 𝑃 𝑃 . 12 12 0 13 24 The next step is decomposing 𝜓 according to Eq. (1.6). For 7-body fragments 𝑎 ≡ 12 with 6-body 12,123 7 fragments 𝑎 ≡ 123 the all relevant 5-body fragments (𝑎 ) are 1234,1235,1236,1237,1238, 123+ 45,123+ 6 5 46,123+ 47,123+ 48,123+ 56,123+ 57,123+ 58,123+ 67,123+ 68,123+ 78. All 5-body sub-clusters summed up 1234 1235 1236 1237 1238 123+45 123+46 123+47 123+48 𝜓 = 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 +𝜓 + 𝜓 (1.19) 12,123 12,123 12,123 12,123 12,123 12,123 12,123 12,123 12,123 12,123 123+56 123+57 123+58 123+67 123+68 123+78 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 . 12,123 12,123 12,123 12,123 12,123 12,123 Considering Eq. (1.6), first part of the above turns 1234 123 123 (1.20) 𝜓 = 𝐺 𝒯 (𝜓 +𝜓 )+ 𝐺 𝒯 (𝜓 +𝜓 ) 12,123 0 12,12 12,124 12,12+34 0 12,23 23,234 23,23+14 + 𝐺 𝒯 (𝜓 +𝜓 ), 0 12,31 31,134 31,31+24 since (1.21) 𝜓 +𝜓 = 𝑃 𝑃 (𝜓 +𝜓 ), 23,234 23,23+14 12 23 12,124 12,12+34 (1.22) 𝜓 +𝜓 = 𝑃 𝑃 (𝜓 +𝜓 ), 31,134 31,31+24 13 23 12,124 12,12+34 according to Eq. (1.13), Eq. (1.20) leads to 1234 123 (1.23) 𝜓 = 𝐺 𝒯 (𝜓 +𝜓 ). 12,123 0 12,124 12,12+34 Similarly, next we decompose 𝜓 according to Eq. (1.6). For the 7-body fragments 𝑎 ≡ 12 with 6- 12,12+34 7 body fragments 𝑎 ≡ 12+ 34 the possible 5-body fragments (𝑎 ) are 1234,125+ 34,126+ 34,127+ 34,128+ 6 5 34,12+ 345,12+ 346,12+ 347,12+ 348,12+ 34+ 56,12+ 34+ 57,12+ 34+ 58,12+ 34+ 67,12+ 34+ 68,12+ 34+ 78, which are now summarized to 1234 123+45 123+46 123+47 123+48 12+345 (1.24) 𝜓 = 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 12+346 12+348 12+34+56 12+34+57 12+34+58 12+347 + 𝜓 +𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 12+34+67 12+34+68 12+34+78 + 𝜓 + 𝜓 + 𝜓 . 12,12+34 12,12+34 12,12+34 Considering Eq. (1.6), the first part of the above equation is regarded in turn 1234 12+34 12+34 (1.25) 𝜓 = 𝐺 𝒯 (𝜓 +𝜓 )+ 𝐺 𝒯 (𝜓 +𝜓 ), 0 12,123 12,124 0 34,234 34,134 12,12+34 12,12 12,34 since we use 𝜓 +𝜓 = 𝑃 𝑃 (𝜓 +𝜓 ). Eq. (1.25) is simplified according to Eq. (1.17), and 34,234 34,134 13 24 12,123 12,124 rewritten as 1234 12+34 (1.26) 𝜓 = 𝐺 𝒯 (𝜓 +𝜓 ). 0 12,123 12,124 12,12+34 1234 1234 The two amplitudes, 𝜓 and 𝜓 , expressed in Eq. (1.23) and Eq. (1.26) are coupled to each other 12,123 12,12+34 as shown now. The expression Eq. (1.19) can easily be converted to 𝜓 and using in addition Eq. (1.24), with 12,124 usage of Eq. (1.23), one obtains 1234 123 1234 1245 1246 1247 1248 124+35 124+36 124+37 𝜓 = 𝐺 𝒯 ((𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 +𝜓 12,123 0 12,124 12,124 12,124 12,124 12,124 12,124 12,124 12,124 124+38 124+36 124+57 124+58 124+67 124+68 124+78 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 ) 12,124 12,124 12,124 12,124 12,124 12,124 12,124 123+45 123+46 123+48 12+345 1234 123+47 + (𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 12+346 12+347 12+348 12+34+56 12+34+57 12+34+58 + 𝜓 +𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 12+34+67 12+34+68 12+34+78 (1.27) + 𝜓 + 𝜓 + 𝜓 )), 12,12+34 12,12+34 12,12+34 correspondingly Eq. (1.26) yields 1234 12+34 1234 1235 1236 1237 1238 123+45 𝜓 = 𝐺 𝒯 ((𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 12,12+34 0 12,123 12,123 12,123 12,123 12,123 12,123 123+46 123+47 123+48 123+56 123+57 123+58 123+67 + 𝜓 +𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 12,123 12,123 12,123 12,123 12,123 12,123 12,123 123+68 123+78 + 𝜓 + 𝜓 ) 12,123 12,123 1234 1245 1246 1247 1248 124+35 124+36 124+37 + (𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 +𝜓 12,124 12,124 12,124 12,124 12,124 12,124 12,124 12,124 (1.28) 124+38 124+36 124+57 124+58 124+67 124+68 124+78 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 )), 12,124 12,124 12,124 12,124 12,124 12,124 12,124 1234 1234 one separates now the components 𝜓 and 𝜓 in Eq. (1.27) from the rest 12,124 12,12+34 1234 123 1234 1234 𝜓 − 𝐺 𝒯 (𝜓 + 𝜓 ) 12,123 0 12,124 12,12+34 123 1245 1246 1247 1248 124+35 124+36 124+37 = 𝐺 𝒯 (𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 +𝜓 0 12,124 12,124 12,124 12,124 12,124 12,124 12,124 124+38 124+36 124+57 124+58 124+67 124+68 124+78 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 12,124 12,124 12,124 12,124 12,124 12,124 12,124 123+45 123+46 123+47 123+48 12+345 12+346 12+347 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 +𝜓 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 12+348 12+34+56 12+34+57 12+34+58 12+34+67 12+34+68 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 (1.29) 12+34+78 + 𝜓 ), 12,12+34 1234 1234 also, one separates now the components 𝜓 and 𝜓 in Eq. (1.28) from the rest 12,123 12,124 1234 12+34 1234 1234 𝜓 − 𝐺 𝒯 (𝜓 + 𝜓 ) 12,12+34 12,123 12,124 12+34 1235 1236 1237 1238 123+45 123+46 123+47 = 𝐺 𝒯 (𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 +𝜓 0 12,123 12,123 12,123 12,123 12,123 12,123 12,123 123+48 123+56 123+57 123+58 123+67 123+68 123+78 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 12,123 12,123 12,123 12,123 12,123 12,123 12,123 1245 1246 1248 124+35 124+36 124+38 1247 124+37 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 +𝜓 + 𝜓 12,124 12,124 12,124 12,124 12,124 12,124 12,124 12,124 124+56 124+57 124+58 124+67 124+68 124+78 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 ), 12,124 12,124 12,124 12,124 12,124 12,124 (1.30) 1234 1234 with 𝜓 = −𝑃 𝜓 Eqs. (1.29) and (1.30) into a matrix form 12,124 34 12,123 1234 1234 123 123 𝜓 𝒯 (−𝑃 ) 𝒯 𝜓 12;123 12;123 ( ) − 𝐺 ( )( ) 1234 12+34 1234 𝜓 𝒯 (1 − 𝑃 ) 0 𝜓 12;12+34 12;12+34 1245 1246 1248 124+35 124+36 124+38 1247 124+37 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 +𝜓 + 𝜓 12,124 12,124 12,124 12,124 12,124 12,124 12,124 12,124 124+36 124+57 124+58 124+67 124+68 124+78 +𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 12,124 12,124 12,124 12,124 12,124 12,124 123+45 123+46 123+47 123+48 12+345 12+346 12+347 12+348 +𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 +𝜓 + 𝜓 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 (1.31) 12+34+56 12+34+57 12+34+58 12+34+67 12+34+68 12+34+78 +𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 ( 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 ) = 𝐺 , 1235 1236 1237 1238 123+45 123+46 123+47 123+48 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 +𝜓 + 𝜓 12,123 12,123 12,123 12,123 12,123 12,123 12,123 12,123 123+56 123+57 123+58 123+67 123+68 123+78 +𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 12,123 12,123 12,123 12,123 12,123 12,123 12+34 1245 1246 1247 1248 124+35 124+36 124+37 124+38 +𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 +𝜓 + 𝜓 12,124 12,124 12,124 12,124 12,124 12,124 12,124 12,124 124+56 124+57 124+58 124+67 124+68 124+78 +𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 ( ) 12,124 12,124 12,124 12,124 12,124 12,124 ( ) since 1245 1235 124+35 123+45 (1.32) 𝜓 = −𝑃 𝜓 ,𝜓 = −𝑃 𝜓 , 12,124 34 12,123 12,124 34 12,123 the right side of Eq. (1.31) can be factored and achieves the form as 1234 1234 123 123 123 123 𝜓 𝜓 𝒯 (−𝑃 ) 𝒯 𝒯 (−𝑃 ) 𝒯 12;123 12;123 34 34 ( ) − 𝐺 ( )( ) = 𝐺 ( ) 0 0 1234 12+34 1234 12+34 𝜓 𝒯 (1 − 𝑃 ) 0 𝜓 𝒯 (1 − 𝑃 ) 0 12;12+34 34 12;12+34 34 1235 1236 1237 1238 123+45 123+46 123+47 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 12,123 12,123 12,123 12,123 12,123 12,123 12,123 ( ) 123+48 123+56 123+57 123+58 123+67 123+68 123+78 (1.33) +𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 12,123 12,123 12,123 12,123 12,123 12,123 12,123 × , 123+45 123+46 123+47 123+48 12+345 12+346 12+347 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 +𝜓 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 ( ) 12+348 12+34+56 12+34+57 12+34+58 12+34+67 12+34+68 12+34+78 +𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 + 𝜓 ( ) 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 12,12+34 the right side can be reduced applying permutation operators and obtains the final form of the 8-nucleon equations as 123 123 𝒯 (−𝑃 ) 𝒯 12;123 (1.34) ( ) = 𝐺 ( ) 1234 12+34 𝜓 𝒯 (1 − 𝑃 ) 0 12;12+34 34 1234 123+45 (1 − 𝑃 −𝑃 −𝑃 −𝑃 )𝜓 + (1 − 𝑃 −𝑃 −𝑃 −𝑃 −𝑃 −𝑃 +𝑃 𝑃 +𝑃 𝑃 +𝑃 𝑃 )𝜓 45 46 47 48 12,123 56 57 58 46 47 48 46 57 46 58 47 57 12,123 × [( )], 1234 125+34 12+345 12+34+56 ( ) 𝜓 + 1 − 𝑃 −𝑃 −𝑃 ( 𝜓 + 𝜓 )+ (1 − 𝑃 −𝑃 −𝑃 −𝑃 +𝑃 𝑃 )𝜓 12;12+34 56 57 58 12;12+34 12;12+34 67 68 57 58 57 68 12,12+34 1234 1234 123+45 we end the formalism with some Yakubovsky independent components as 𝜓 , 𝜓 , 𝜓 , 12;123 12;12+34 12,123 125+34 12+345 12+34+56 (𝜓 + 𝜓 ) and 𝜓 coupled in the equations. According to the 5-body sub-cluster distributions of 12,12+34 12,12+34 12,12+34 above independent components, it is well-known that two first components are related to 4𝑁 + 𝑁 + 𝑁 + 𝑁 + 𝑁 that refer to approximately effective 5-body model (𝛼 + 𝑛 + 𝑛 + 𝑛 + 𝑛 ). Respectively, third component refers to 3𝑁 + 2𝑁 + 𝑁 + 𝑁 + 𝑁 , which is not present in the effective 5-body 𝛼 -core model. Fourth component is a linear combinations that refer again to the third component. This is also beyond the effective 5-body 𝛼 -core model. The last one is related to 2𝑁 + 2𝑁 + 2𝑁 + 2𝑁 sub-cluster which is not contained in the effective 5-body 𝛼 -core model, either. 1234 1234 Regarding the sub-cluster underlying the Yakubovsky components, only 𝜓 and 𝜓 , are related to 12;123 12;12+34 approximately effective configuration of 𝛼 + 𝑛 + 𝑛 + 𝑛 + 𝑛 , where the 𝛼 -core and four loosely bound neutrons approximation are valid. Now, here we convinced enough why we choose two specific components, and those components are related to the effective 𝛼 -core structure. By considering the distinct Jacobi configurations related to 1234 1234 remaining Yakubovsky components, Fig.1, we can choose the first two components, 𝜓 and 𝜓 , and 12;123 12;12+34 obviously other components will not be taken into account. Hint: after removing the contribution interactions of the seventh and eighth nucleons in Eq. (1.34), i.e. removing their permutation operators, the 8-nucleon obtained equations leads to the 6-nucleon ones for halo nucleus He (See Eq. (67) in Ref. [12]). Such reduction, confirms that our implementing the Yakubovsky formalism for the 8- nucleon bound system for He, is appropriate approximation to describe the effective 𝛼 -core structure for Helium halo-bound nuclei, similar to the approximation of the six-nucleon for He in the effective 3-body model [12, 13]. In the above explanations we clearly have discussed why only the two specific components, 𝜓 and 12;123 𝜓 , are relevant to the effective 𝛼 -core structure of He, and the other Yakubovsky components in Eq. (1.34), 12;12+34 will not be taken into account. Therefore, Eq. (1.34) leads to 1234 1234 123 123 𝜓 (1 − 𝑃 −𝑃 −𝑃 −𝑃 )𝜓 𝒯 (−𝑃 ) 𝒯 12;123 45 46 47 48 12,123 ( ) = 𝐺 ( )[( )]. (1.35) 1234 12+34 1234 𝜓 𝒯 (1 − 𝑃 ) 0 𝜓 12;12+34 34 12;12+34 It is worth mentioning that switching off some Yakubovsky independent components changes the Hamiltonian of the comprehensive 8-nucleon system, so the Hamiltonian does not contain all interactions, but we show that in the remaining components the effective interactions of 𝛼 -core as the stable 4-nucleon sub-system is governed in the remaining components (Refer to Fig.1). It is obvious that after removing the contribution interactions of the 5 to 8’s nucleons in the 8-nucleon Yakubovsky equations, Eq. (1.34) straightforwardly reduces to the four- nucleon Yakubovsky equations for He ( See Ref. [10]) as 1234 1234 123 123 𝜓 𝒯 (−𝑃 ) 𝒯 𝜓 12;123 12;123 ( ) = 𝐺 ( )( ), (1.36) 1234 12+34 1234 𝜓 𝒯 (1 − 𝑃 ) 0 𝜓 12;12+34 12;12+34 in order to reduce the elaboration of the 8-nucleon coupled equations, Eq. (1.35), we can present the Yakubovsky 123 12+34 equations for the 8-nucleon problem as a function of two-nucleon 𝑡 -matrices and avoid using 𝒯 and 𝒯 . Consequently, the formalism can be simplified and the numerical solution of integral would be much faster in 123 12+34 comparison with the case where was need to solve first sub-cluster Faddeev-like equation to obtain 𝒯 and 𝒯 by using Pade′ approximation [10]. In addition to Eq. (1.36), the 4-nucleon Yakubovsky coupled equations has been obtained in [10] as 𝜓 𝐺 𝑡 𝑃 (1 − 𝑃 ) 𝐺 𝑡 𝑃 𝜓 1 0 12 34 0 12 1 (1.37) ( ) = ( )( ). ̃ ̃ 𝜓 𝐺 𝑡 𝑃 𝜓 𝐺 𝑡 𝑃 (1 − 𝑃 ) 2 0 12 2 0 12 34 also in Eq. (1.36) the 4-nucleon transition operators obey as 123 (1.38) = 𝑡 𝑃 + 𝑃 𝐺 𝑡 ≡ 𝒯 , 12 0 12 12+34 ̃ ̃ ̃ ̃ ̃ ̃ (1.39) 𝑇 𝑃 = 𝑡 𝑃 + 𝑃 𝐺 𝑡 𝑇 𝑃 ≡ 𝒯 , 12 0 12 According to [10] the 4-nucleon matrix form transition in Eq. (1.36) and (1.37) are similar in action, and according to Eqs. (1.35) and (1.36), the 8-nucleon and the 4-nucleon matrix form transitions are equal in action. At this point, the matrices of Eq. (1.37) can be identified and Eq. (1.35) can be written in terms of two-nucleon t−matrices operator since they are given as 1234 1234 𝜓 (1 − 𝑃 −𝑃 −𝑃 −𝑃 )𝜓 𝐺 𝑡 𝑃 (1 − 𝑃 ) 𝐺 𝑡 𝑃 12;123 45 46 47 48 12,123 0 12 34 0 12 ( ) = ( )[( )], 1234 1234 ̃ ̃ 𝐺 𝑡 𝑃 (1 − 𝑃 ) 𝐺 𝑡 𝑃 𝜓 𝜓 0 12 (1.40) 12;12+34 0 12 34 12;12+34 finally, the linear form of the obtained 8-nucleon Yakubovsky equations in the case of effective 𝛼 -core structure can be written as 1234 1234 1234 (1.41) 𝜓 = 𝐺 𝑡 𝑃 [(1 − 𝑃 )(1 − 𝑃 −𝑃 −𝑃 −𝑃 )𝜓 + 𝜓 ], 12;123 0 12 34 45 46 47 48 12;123 12;12+34 1234 1234 1234 (1.42) 𝜓 = 𝐺 𝑡 𝑃 [(1 − 𝑃 )(1 − 𝑃 −𝑃 −𝑃 −𝑃 )𝜓 + 𝜓 ]. 12;12+34 0 12 34 45 46 47 48 12;123 12;12+34 In the next step we will describe characteristic numerical techniques such as introducing Jacobi configuration of each independent component in momentum space representation to provide corresponding basis states and to understand the halo structure model of the 8-nucleon system with respect to the regarded relevant remaining components; integral representation of the coupled equations; eigenvalue equation form and iteration method for solution of it, and discussions about the halo structure configurations that describe the 8-nucleon halo-bound system in the case of effective 𝛼 -core structure. II. Technicalities for Numerical Applications In this section, first, Jacobi configurations in momentum space are represented to define corresponding momentum basis states for evaluations of the integral kernels. Next, details for a typical algorithm are described to solve coupled integral equations based on PW representation. In order to describe the corresponding momentum basis 8 1234 1234 states, for two independent components of the dominant 𝛼 -core structure of He, namely 𝜓 and 𝜓 the 12;123 12;12+34 standard Jacobi momenta are described in Fig. 1, respectively. 𝑇𝑃 𝑇𝑃 (1) (2) Fig. 1: Schematic representation of the two independent components of 8-nucleon bound system for He (four loosely bound neutrons with an inert 𝛼 -core) in Jacobi configurations. In the bag the 𝛼 -core tightly interacts in a 4-nucleon subsystem. For the first component 𝜓 in terms of the first configuration in Fig. 1 sets corresponding to 𝒖 -set chains 12;123 1/2 −1/2 0 0 0 0 0 0 −1/3 −1/3 2/3 0 0 0 0 0 2 2 0 0 0 0 −1/4 −1/4 −1/4 3/4 𝒑 1/2 −1/2 𝒑 𝒖 0 0 0 0 0 0 4 = , 0 0 5 5 −1/3 −1/3 −1/3 −1/3 2/3 2/3 1/2 −1/2 0 0 0 0 0 0 𝒖 𝒑 7 7 (2.1) −1/4 −1/4 −1/4 −1/4 −1/4 −1/4 3/4 3/4 ( ) (𝒑 ) ( ) 1 1 1 1 1 1 1 1 in the non-relativistic case, we can define the kinetic energy operator by two equivalent forms. The inverse form of the above-mentioned transfer matrix is used to represent the kinetic energy in terms of 𝒖 -set Jacobi momenta. 2 2 2 2 2 2 2 2 𝑝 𝑢 3𝑢 2𝑢 𝑢 3𝑢 𝑢 1𝑢 𝑖 1 2 3 4 5 6 7 𝐻 = ∑ ≡ + + + + + + , (2.2) 2𝑚 𝑚 4 𝑚 3 𝑚 𝑚 8 𝑚 𝑚 3 𝑚 𝑖 =1 where 𝒑 is an individual particle momentum in the center of mass situation (under the condition ∑𝒑 = 0) that is 𝑖 𝑖 described by relative Jacobi momenta 𝒖 ;( 𝑖 = 1,2,3,4,5,6,7,8). Similarly, to 𝜓 in terms of the second 𝑖 12;12+34 configuration in Fig. 1 belongs 0 0 0 0 0 0 1/2 −1/2 1/2 −1/2 0 0 0 0 0 0 𝒑 0 0 0 0 −1/2 −1/2 1/2 1/2 𝒑 1/2 −1/2 𝒑 𝒗 0 0 0 0 0 0 4 = , (2.3) 0 0 5 5 −1/3 −1/3 −1/3 −1/3 2/3 2/3 1/2 −1/2 0 0 0 0 0 0 𝒗 𝒑 7 7 −1/4 −1/4 −1/4 −1/4 −1/4 −1/4 3/4 3/4 ( ) (𝒑 ) ( ) 1 1 1 1 1 1 1 1 according to some appearing Jacobi vectors in the two different configurations in Fig.1, correspondingly, some equal Jacobi momenta have been represented by different notations, i.e. 𝒗 = 𝒖 , because of the standard representation of 7 7 the above transfer matrices, for providing inverse transformations. The kinetic energy form in terms of 𝒗 -set Jacobi momenta, are given as 2 2 2 2 2 2 2 2 𝑝 𝑣 𝑣 1𝑣 𝑣 3𝑣 𝑣 1𝑣 𝑖 1 2 3 4 4 4 4 𝐻 = ∑ ≡ + + + + + + . 2𝑚 𝑚 𝑚 2 𝑚 𝑚 8 𝑚 𝑚 3 𝑚 (2.4) 𝑖 =1 Then, we introduce the two basis-states corresponding two independent components in PW representation. Obviously, these basis are separated into the 3 spaces that are 1-momentum/ 2-angular momentum and spin/ 3-isospin. 1234 1234 The relevant basis-states in PW representation for 𝜓 and 𝜓 respectively are given as 12;123 12;12+34 |𝑢 ; 𝛾 ⟩ ≡ |𝑢 𝑢 𝑢 𝑢 𝑢 𝑢 𝑢 ⟩ 𝑢 1 2 3 4 5 6 7 1 1 | ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ⨂ (𝑙 𝑠 𝑗 (𝑙 )𝑗 𝑗 𝑗 𝐼 (𝑙 )𝑗 𝐼 𝑗 𝐼 𝑙 𝑠 𝑗 𝑙 𝑗 𝐼 𝐼 𝐼 𝐼 𝑙 𝑠 𝑗 𝑙 𝑗 𝐼 𝐼 𝐼 𝐽 ,𝑀 ⟩ 1 12 1 2 2 1 2 3 3 4 3 4 4 5 56 5 4 5 5 4 5 6 6 12 6 7 6 7 6 7 𝑗 2 2 1 1 (2.5) ⨂|(𝑡 )𝑡 (𝑡 )𝑡 (𝑡 𝑡 )𝑡 (𝑡 𝑡 ) 𝑇 ,𝑀 ⟩, 12 3 3 4 4 56 6 6 78 𝑡 2 2 ⟩ ⟩ |𝑣 ; 𝛾 ≡ |𝑣 𝑣 𝑣 𝑣 𝑣 𝑣 𝑣 𝑣 1 2 3 4 5 6 7 ⨂|(𝑙 𝑠 )𝑗 (𝑙 𝑠 )𝑗 (𝑗 𝑗 )𝑆 (𝐿 𝑆 )𝐼 (𝑙 𝑠 )𝑗 (𝑙 𝑗 )𝐼 (𝐼 𝐼 )𝐼 (𝑙 𝑠 )𝑗 (𝑙 𝑗 )𝐼 (𝐼 𝐼 ) 𝐽 ,𝑀 ⟩ 1 12 1 2 34 2 1 2 4 56 5 5 5 5 5 6 4 78 6 7 6 7 6 7 𝑗 (2.6) | ) ( ) ⨂ (𝑡 𝑡 𝑡 (𝑡 𝑡 )𝑡 𝑡 𝑡 𝑇 ,𝑀 ⟩, 12 34 1_4 1_4 56 6 6 78 𝑡 here the first part in the above-introduced basis states are magnitude of momentum vectors in PW representation. In the second part the orbital angular momentum 𝑙 go with the 𝒖 \ 𝒗 , 𝑠 which are two-body spins for nucleons , 𝑗 𝑖 𝑖 𝑖 𝑖 𝑗 𝑖 are total 1- and 2-body angular momenta coupled out of orbital and spin angular momenta, 𝐼 are total 𝑖 -body angular momenta except 𝐼 and 𝐼 . These are total angular momenta of 5, 6 and 7, 8 respectively. Evidently, the third part 5 7 refers to isospin. In this step, in order to implement the numerical techniques, the two coupled equations, Eq. (1.41) and Eq. (1.42), by inserting the completeness relations between the permutation operators, can be projected in momentum space as 1234 ′2 ′ ′′2 ′′ ⟨𝑢 ; 𝛾 |𝜓 ⟩ = ∑∫𝑢 𝑑 𝑢 ∑ ∫𝑢 𝑑 𝑢 12;123 𝛾 𝛾 ′′ 𝑢 ′ ′′ ′′ 1234 ⟨𝑢 ; 𝛾 |𝐺 𝑡 𝑃 |𝑢 ′; 𝛾 ⟩⟨𝑢 ′; 𝛾 |(1 − 𝑃 )(1 − 𝑃 −𝑃 −𝑃 −𝑃 )| 𝑢 ; 𝛾 ′′⟩⟨ 𝑢 ; 𝛾 ′′|𝜓 ⟩ 𝑢 0 12 𝑢 ′ 𝑢 ′ 34 45 46 47 48 12;123 𝑢 𝑢 ′2 ′ ′2 ′ ′ ′ 1234 + ∑∫𝑢 𝑑 𝑢 ∑∫𝑣 𝑑 𝑣 ⟨𝑢 ; 𝛾 |𝐺 𝑡 𝑃 |𝑢 ′; 𝛾 ⟩⟨𝑢 ′; 𝛾 |𝑣 ; 𝛾 ′⟩⟨𝑣 ; 𝛾 ′|𝜓 ⟩, 𝑢 0 12 𝑢 ′ 𝑢 ′ 12;12+34 𝑣 𝑣 𝛾 𝛾 𝑢 ′ 𝑣 ′ (2.7) 1234 ′2 ′ ′2 ′ ⟨𝑣 ; 𝛾 | 𝜓 ⟩ = ∑∫𝑣 𝑑 𝑣 ∑∫𝑢 𝑑 𝑢 𝑣 12;12+34 𝛾 𝛾 𝑣 ′ 𝑢 ′ ′ ′ 1234 ⟨𝑣 ; 𝛾 |𝐺 𝑡 𝑃 |𝑣 ; 𝛾 ′⟩⟨𝑣 ; 𝛾 ′|(1 − 𝑃 )(1 − 𝑃 −𝑃 −𝑃 −𝑃 )|𝑢 ′; 𝛾 ⟩⟨𝑢 ′; 𝛾 |𝜓 ⟩ 𝑣 0 12 34 45 46 47 48 𝑢 ′ 𝑢 ′ 12;123 𝑣 𝑣 (2.8) ′2 ′ ′ ′ 1234 + ∑∫𝑣 𝑑 𝑣 ⟨𝑣 ; 𝛾 |𝐺 𝑡 𝑃 |𝑣 ; 𝛾 ′⟩⟨𝑣 ; 𝛾 ′|𝜓 ⟩, 𝑣 0 12 𝑣 𝑣 12;12+34 𝑣 ′ using techniques like ones presented in Refs. [10-14], it is straightforward to evaluate in PW analysis like the above- mentioned kernels and shifted momenta. It is well-known that the variables for 7 amplitudes on the right side are in general linear combinations of intermediate integration variables which include angles besides momentum magnitudes. After evaluating each term in the above-mentioned coupled integral equations in the standard PW analysis, the obtained equations are the starting point for numerical calculations in the eigenvalue equation form. We can replace the continuous variables in the numerical treatment by a dependence on certain distinct values by using the Gauss-Legendre discretization. The eigenvalue equation can be solved by the iteration method. It can apply Lanczos-like scheme that is efficient for nuclear few-body problems [10]. Also the evaluated coupled integral 𝑖𝑗 equations involve a very large number of interpolations, it can be used the cubic Hermit splines of Ref. [17] for its accuracy and high computational speed. The eigenvalue equation form can be given as 1234 1234 1234 1234 (2.9) 𝜂 (𝐸 ) 𝜓 ( 𝜓 , 𝜓 )= 𝑘 (𝐸 ) 𝜓 ( 𝜓 , 𝜓 ), 12;123 12;12+34 12;123 12;12+34 ( ) where 𝐸 is the energy eigenvalue at which auxiliary Yakubovsky kernel eigenvalue is 𝜂 𝐸 = 1. Because the energy is varied such that one reaches eigenvalue to be one. The Yakubovsky kernel of the equations 𝑘 (𝐸 ) is just depend to 1234 1234 energy E, and 𝜂 (𝐸 ) is eigenvalue with 𝜓 ( 𝜓 , 𝜓 ) as the corresponding eigenvector. In order to solve the 12;123 12;12+34 eigenvalue equation, Eq. (2.9), it can be used the Gaussian quadrature grid points. The coupled equations represent a set of homogenous integral equations, which turn into a high-dimension matrix eigenvalue equation after discretization. Starting from an arbitrary initial 𝜓 ≡ 𝜓 one generates by consecutive applications of 𝑘 (𝐸 ) a sequence of amplitudes 𝜓 , which after orthogonalization form a basis into which 𝜓 is expanded. Then the energy is varied ( ) such that one reaches 𝜂 𝐸 = 1. In order to investigate the loosely interactions of the halo bound neutrons in He, −1 calculations of the momentum probability density as 𝑛 (𝑢 )[𝑓𝑚 ] in terms of 𝑢 [𝑓𝑚 ], correspondingly (𝑣 ) can be 𝑖 𝑖 𝑖 done. Why we just select 4 and 6 Jacobi momenta to calculate the momentum probability density? Because they are the only ones that make up contributions of the halo neutrons (refer to Eqs. (2.1) and (2.3)). By calculating momentum probability densities, we can represent the situation of the halo bound neutrons in momentum space. It is expected that the representation of momentum probability densities approximately are wide and short in comparison with tightly bound neutrons, namely neutrons inside the 𝛼 -core subsystem. Such those momentum distributions confirms that the 8-nucleon in the case of effective 𝛼 -core model is a halo structure system, exactly like halo nucleus He. III. Summary The applications of Faddeev and Yakubovsky equations are traditionally for 3- and 4-nucleon bound and scattering states. Today, after the experiences with 4-, 5- and 6-nucleon bound-state problems within the Yakubovsky approach [13, 14] that the technical expertise has developed in the recent years and the very strong increase of computer power just recently achieved, it is possible to approach the 8-nucleon bound-state problem in that formalism and then numerical solution. For a full solution of the 8-nucleon Yakubovsky equations for He, modern super- computers with parallel algorithms is required. But in order to approach the calculations of the 8-nucleon bound system, for the first time, we have implemented the 8-nucleon Yakubovsky equations for the halo nucleus He in the case of simple effective 𝛼 -core structure. To this end, we stop the sequential sub-clustering with 5 fragments, namely 𝛼 + 𝑛 + 𝑛 + 𝑛 + 𝑛 , though the additional step with 4 fragments could be easily performed, but is not applicable for a physical system. In this method, we restricted the formulation to two-nucleon forces only because loosely bound neutrons with the center of an inert 𝛼 -core tend to have much stronger pairing interactions for large radii neutrons [5]. Therefore, as a simplification in the formalism, and above fact in such a model, three-nucleon forces was ignored. As a result, the implementation of the 8-nucleon bound system in the halo structure model within the Yakubovsky approach leads to two coupled equations in terms of two independent Yakubovsky components. Subsequently, in order to approach the halo structures of He, in addition to the analytical implementations, we have introduced the Jacobi configurations of those components in momentum space representation. Eventually, we provide technicalities which were considered useful for a numerical performance, such as bound-state calculations and proving the halo structure of bound neutrons that are outside of the 𝛼 -core. References [1] M. Brodeur, et. al. Phys. Rev. Lett. 108, 052504 –31 Jan (2012). [2] M.V.Zhukov, et. al. Physics Reports, Volume 231, Issue 4, August (1993), Pages 151-199. [3] S. Bacca, A. Schwenk, G. Hagen, et. al. Eur. Phys. J. A, 42: 553 (2009). [4] L. B. Wang et al., Phys. Rev. Lett. 93, 142501 (2004). [5] P. Mueller et al., Phys. Rev. Lett. 99, 252501 (2007). [6] V. L. Ryjkov et al., Phys. Rev. Lett. 101, 012501 (2008). [7] S. C. Pieper and R. B. Wiringa, Annu. Rev. Nucl. Part. Sci. 51, 53 (2001); S. C. Pieper, arXiv:0711.1500. [8] P. Navratil and W. E. Ormand, Phys. Rev. C 68, 034305 (2003). [9] P. Navratil, V. G. Gueorguiev, J. P. Vary, W. E. Ormand and A. Nogga, Phys. Rev. Lett. 99, 042501 (2007). [10] 10. H. Kamada and W. Glӧckle, Nucl. Phys. A 548, 205 (1992). [11] A. Nogga, H. Kamada and W. Glöckle, Phys. Rev. Lett. 85, 944 (2000). [12] W. Glӧckle and H. Witala, Few-Body Syst. 51, 27-44 (2011). [13] E. Ahmadi Pouya and A. A. Rajabi, Acta, Phys, Pol, B 48: 1279 (2017). [14] E. Ahmadi Pouya and A. A. Rajabi, Eur. Phys. J. Plus, 131: 240 (2016). [15] A. C. Fonseca, Phys. Rev. C 30, 35 (1984). [16] W. Glöckle: The Quantum Mechanical Few-Body Problem. Springer-Verlag, New York (1983). [17] D. Huber, H. Witala, A. Nogga, W. Glӧckle and H. Kamada, Few-Body Syst. 22, 107 (1997) .

Journal

Nuclear Theory – arXiv (Cornell University)

Published: Apr 29, 2019

There are no references for this article.