{'fraction_non_alphanumeric': 0.042434715821812595, 'fraction_numerical': 0.022625448028673834, 'mean_word_length': 4.322602623062511, 'pattern_counts': {'":': 0, '<': 12, '<?xml version=': 0, '>': 12, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 1, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 5, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Microscopic three-cluster study of light exotic nuclei

I develop a microscopic three-cluster model for exotic light nuclei. I use the hyperspherical formalism, associated with the Generator Coordinate Method. This model is well adapted to halo nuclei, since the long-range part of the radial wave functions is accurately reproduced. The core wave functions are described in the shell model, including excited states. This technique provides large bases, expressed in terms of projected Slater determinants. Matrix elements involve seven-dimension integrals, and therefore require long calculation times. I apply the model to 11Li, 14Be, 15B, and 17N described by two neutrons surrounding a 9Li, 12Be, 13B and 15N core, respectively. The 17Ne (as 15O+p+p) and 15Ne (as 13O+p+p) mirror nuclei are briefly discussed. I present the spectra and some spectroscopic properties, such as r.m.s. radii or E2 transition probabilities. I also analyze the importance of core excitations.


I. INTRODUCTION
Exotic nuclei represent one of the major interests in modern nuclear physics [1]. These nuclei are located close to the driplines (neutron or proton) and, owing to the low binding energy, present specific properties, such as an anomalously large radius or modifications of the shell structure. Recent development of experimental facilities have provided a large number of new data, which need to be understood by theory.
The main property of exotic nuclei being their low binding energy, the relative wave functions extend to large distances. Theoretical models therefore need to reproduce this property. A widely used approach is the hyperspherical method [2], where the Jacobi coordinates are replaced by a set of angles, and by a single length, referred to as the hyperradius. The three-body equation is then replaced by a set of coupled differential equations depending on the hyperradius. The hyperspherical formalism can be extended to systems involving more than three particles [3]. Many works in atomic and in nuclear physics have been performed within this method.
In nuclear physics, most applications are carried out in non-microscopic models. In other words, the nucleus is seen as a three-body system, with a structureless core, and two external nucleons. Typical applications are the 6 He and 11 Li nuclei, modeled by α+n+n and 9 Li+n+n three-body structures. This approach therefore relies on nucleon+nucleon and nucleon+core potentials. It simulates the Pauli principle by an appropriate choice of these interactions. Core excitations are in general neglected.
This non-microscopic model can be extended to microscopic theories, where the system is described by a A-body Hamiltonian. The core nucleus is defined in the shell model, and a full antisymmetrization is taken into account. The use of the hyperspherical formalism guarantees the correct long-range behaviour of the wave func- * pdesc@ulb.ac.be tion. The model only relies on a nucleon-nucleon interaction, and on the description of the core wave functions. In contrast with non-microscopic models, it allows to include core excitations without any further parameter.
Microscopic cluster models are used in nuclear physics for many years (see reviews in Refs. [4,5]) but most applications deal with the two-cluster variant, much easier than multicluster approaches. The application of the hyperspherical formalism to microscopic cluster theories is recent. The first works were focused on 6 He [6,7] which involves the α particle as a core. The α particle can be accurately described by a (0s) 4 shell model configuration which makes the calculations relatively easy. The model was extended to three-body continuum states [8], and microscopic 6 He wave functions have been used in CDCC calculations of elastic scattering on heavy targets [9].
The main limitation of microscopic three-cluster models in hyperspherical coordinates is that the matrix elements involve the numerical calculation of (many) sevendimension integrals [7]. Consequently, applications are essentially limited to light systems, most of them involving an α core. In addition to the 6 He system mentioned above, a recent application focuses on 8 Li, described by an α+ 3 H+n three-cluster configuration [10].
The aim of the present work is to go beyond this limitation. The development of the computing technology, and in particular of the parallelization possibilities permits to consider heavier systems, involving p-shell nuclei as the core. I analyze four nuclei ( 11 Li, 14 Be, 15 B, 17 N) which are modelled by a p-shell core and two surrounding neutrons. Compared to previous applications, the present calculations face two difficulties: (1) the presence of p-shell orbitals, and (2) the need of several Slater determinants for the core.
The paper is organized as follows. In Sec. II, I present a general overview of the microscopic three-cluster model. In Sec. III, I apply the method to various exotic nuclei. Concluding remarks are presented in Sec. IV.

A. General overview
My main goal is to solve a A-body problem, where the Hamiltonian is given by In this definition, t i is the kinetic energy of nucleon i, and v ij a nucleon-nucleon interaction which contains nuclear and Coulomb components. Recent ab initio models (see, for example, Ref. [11]) are developed to find exact solutions of the Schrödinger equation associated with (1). These models, however, are in general not well adapted to halo nuclei, where the long-range part of the wave functions plays a crucial role. My approach, in contrast, is based on the cluster approximation [4,12]. The total wave function, solution of the Schrödinger equation is written schematically as where φ i are the internal wave functions of the clusters (in practice, they are defined in the shell model), and G is a radial function, depending on the relative coordinates between the clusters. The antisymmetrization operator A ensures that the wave function is completely antisymmetric.
The cluster approximation permits a strong simplification of the calculations, in comparison with ab initio models. It is also well adapted to halo nuclei or, more generally, to states presenting a strong cluster structure. Owing to this approximation, however, effective nucleonnucleon interactions v ij , such as the Volkov [13] or the Minnesota [14] force, must be used. These forces simulate missing effects, such as the tensor component, by an effective central interaction. Ideally, three-body forces should be introduced, and developments have been achieved in that direction for realistic forces (see Ref. [15] for a review). For effective interactions, however, even if some work has been done [16,17], the use of three-body forces is essentially limited to well bound nuclei, such as 12 C or 16 O. In most cluster calculations, three-body forces are therefore neglected.
Several variants of microscopic cluster models exist. Of course, the simplest version is a two-cluster model, which is used since many years [18]. Three-cluster models have been developed in various directions: frozen triangular configurations [19], 2+1 configurations essentially aimed for nucleus-nucleus scattering such as 7 Be+p [20], and, more recently, genuine three-body models using the hyperspherical formalism [7]. The present work is based on the hyperspherical approach, which is described in more detail in the next subsections.

B. Core wave functions
Let me consider Eq. (3) where the internal wave functions φ i are associated with three clusters with nucleon numbers (A 1 , A 2 , A 3 ) and charge numbers (Z 1 , Z 2 , Z 3 ). I assume that clusters 2 and 3 are s-shell nuclei and, more specifically in the present work, that they are neutrons. I also assume that the oscillator parameter b is common to the three clusters. Until now, the use of the hyperspherical formalism in microscopic models was limited to an α core. This is well adapted, for example, to 6 He[7], 6 Li [7], or 8 B [10]. The main reason for this limitation is that, as it will be discussed later, matrix elements involve seven-dimensional integrals.
This limitation, however, restricts the applications to a few light nuclei. In the present work, I go beyond this limitation, and extends the model to p-wave orbitals. The consequences in terms of computer times are twofold: (i) the matrix elements involve p-shell orbitals which has a strong impact on the matrix elements of the nucleonnucleon interaction (quadruple sums over the individual orbitals); (ii) p-shell nuclei such as 9 Li or 12 Be involve several Slater determinants (up to 90 for 9 Li), whereas the α particle is described by a single Slater determinant.
I remind here the main properties. Let me consider a Slater determinantΦ i involving A 1 single-particle orbitals where ϕ n (s) are harmonic oscillator functions and where m s , m t = ±1/2 are associated with the spin and the isospin, respectively. The first step is to define the list of N S Slater determinants consistent with the Pauli principle. Assuming that the s-shell is filled, the number of Slater determinants for p-shell nuclei is , where C j i is the number of combinations of j elements among i elements. For 9 Li, 12 Be, 13 B and 15 N, I have N S = 90, 15, 20, 6, respectively. The N S Slater determinantsΦ i must be projected on the various angular momenta I, L, S, T from a diagonalization of the I I I 2 , L L L 2 , S S S 2 and T T T 2 operators (I I I = L L L + S S S). This procedure provides basis functions with good quantum numbers as Finally, basis functions (5) are used to diagonalize the Hamiltonian (1) for cluster 1, and I get the core wave functions as This technique corresponds to a standard shell-model approach, and can be extended to higher shells (the summation may involve an additional quantum number, associated with the degeneracy). The specificity of the cluster model is that (6) is only the very first step of the calculations. These internal wave functions must then be introduced in the three-body wave functions (3).

C. Three-cluster wave functions
Let me come back to the three-cluster wave functions (3). In the Generator Coordinate Method (GCM, see Refs. [4,18]), the clusters are located at three points, as represented in Fig. 1. The generator coordinates R R R 1 is associated with the external clusters (neutrons in the present applications) and R R R 2 with the relative motion between the core and the c.m. of clusters 2 and 3. They are variational parameters, and are not associated with the nucleon coordinates (in other words, the antisymmetrization operator does not act on R R R 1 and R R R 2 ). A GCM basis state is defined as where Φ ν1 c1 are the core wave functions (6), and Φ ν2 2 and where Φ ν3 3 correspond to the external clusters. In this equation, A 12 = A 1 +A 2 and A 23 = A 2 +A 3 . If the internal wave functions are defined in the harmonic oscillator model with a common oscillator parameter b, these basis functions are Slater determinants, and the c.m. motion is factorized exactly. For the sake of clarity, I do not write the spins I i of the clusters. A first angular-momentum coupling is performed on the spins I 2 and I 3 of the external clusters with Then I introduce the hyperspherical formalism, well known in non-microscopic models [21], and extended more recently to microscopic theories [7,22]. I define scaled Jacobi coordinates as which provide the hyperradius R and hyperangle α from Both coordinates are complemented by the angles associated with X X X and Y Y Y as From the basis functions (8), I project on the total angular momentum J and parity π. This is achieved by a double projection on the orbital momenta ℓ x and ℓ y ; a projected basis function is written as which represents a four-dimension integral. This doubleprojection technique has been used in three-body models aimed at studying nucleus-nucleus scattering such as 7 Be+p [23] for example. In the hyperspherical formalism, one introduces the hypermoment K, which is a generalization of the angular momentum in three-body systems. An hyperspherical basis function reads where index γ stands for γ = (c 1 I 23 Iℓ x ℓ y L). In this equation, function ϕ ℓxℓy K (α) is defined by P a,b n (x) being a Jacobi polynomial. The normalization coefficient N ℓxℓy K is given by where n = (K −ℓ x −ℓ y )/2 is a positive integer. The main advantage of the hyperspherical formalism is to reduce the number of degrees of freedom to one. The configuration space is therefore spanned by a single generator coordinate R. This approach is also well adapted to threebody continuum states [8]. Of course, matrix elements between basis states (13) are more demanding in terms of computer times, since they involve seven-dimensional integrals (see next subsection) which must be performed numerically.
The total wave function of the system is given by a superposition of basis functions (13) as where N is the number of generator coordinates (typically N ≈ 10). In this expansion, the summation over the hypermoment K is limited to a maximum value K max . As usual in hyperspherical models, K max must be large enough to ensure the convergence of the physical quantities (energies, radii, etc.). Coefficients f Jπ γKi represent the generator function, and are obtained from a diagonalization of the Hamiltonian kernel.

D. Matrix elements
The main part of the calculation concerns matrix elements between projected basis functions (13). The Hamiltonian kernel reads wherẽ withΦ Jπ γ (R 1 , R 2 ) given by Eq. (12). The projection over the hypermoment therefore represents a double integration over α and α ′ . This quadrature is performed numerically. From Eq. (12), the matrix elements (18) involve eight angular integrals. In fact, owing to the rotation invariance, three angles can be fixed, which leads to fivedimensional integrals. The calculation provides where R R R 2 is along the z axis and R R R ′ 2 is in the xz plane. This expression is also valid for other rotation-invariant operators. A generalization to operators, such as the electromagnetic operators, is straightforward.
The calculation of the matrix elements (18) is therefore performed in several steps:

Matrix elements between Slater determinants
AΦ i Φ ν2 2 Φ ν3 3 are first computed. One-body operators involve double sums over the individual orbitals, whereas two-body operators involve quadruple sums. When dealing with three s clusters, the quadruple sums contain 3 4 terms, but this amounts to 6 4 when one cluster belongs to the p shell. Compared to previous works on 6 He or 12 C this property makes the calculations 16 times longer. (7) are then constructed with the transformation coefficients (6). The number of three-cluster Slater determinants is given by N S (2I 2 + 1)(2I 3 + 1). For 3. The next step is to compute the projected matrix elements (19) which involve five numerical quadratures (typically ∼ 16 − 20 points are used for each angle).

Matrix elements between wave functions
4. Finally the projection over the hypermomentum is performed with (17). For a given set of generator coordinates (R, R ′ ), all integrals are evaluated simultaneously.
This process must be repeated for all sets of (R, R ′ ) values. For systems such as 11 Li which involve many Slater determinants, the full calculation is extremely time consuming. This can be achieved with an optimization of the codes, and using modern computing facilities. In practice, a parallelization is performed over the hyperangles α and α ′ in (17). Let me also mention that the total number of basis functions (i.e. the size of the matrices) can be as large as 20000. This raises precision issues since the basis is (highly) non orthogonal.

III. APPLICATION TO LIGHT EXOTIC NUCLEI
A. Conditions of the calculations The calculations are performed with the Volkov V2 nucleon-nucleon interaction [13]. For the oscillator parameter, I use b = 1.60 fm, a standard value for p-shell nuclei. The optimal b value stems from a compromise: to reproduce the core radius, and to minimize the binding energy of the core. Small variations of b can be compensated by a slight readjustment of the nucleon-nucleon interaction. The first step is to determine the core wave functions (6) Values of the generator coordinate R are selected from 1.5 to 15 fm with a step of 1.5 fm. For the maximum hypermomentum, I use K max = 16. These conditions are sufficient to ensure the convergence of the energies and r.m.s. radii.
The central V2 interaction is complemented by a zerorange spin-orbit force with amplitude S 0 [7]. This amplitude is fixed to 40 MeV.fm 5 , except for 11 Li, where I use S 0 = 50 MeV.fm 5 . The Volkov potential contains the admixture parameter M , whose standard value is M = 0.6. This value, however, can be slightly modified without changing the fundamental properties of the interaction. The aim is to reproduce exactly the groundstate energies. For weakly bound nuclei, the properties are sensitive to the long-range part of the wave function, and therefore to the binding energies. With M = 0.7750, 0.6119, 0.5945 and 0.6177, I reproduce the experimental two-neutron separation energies S 2n of 11 Li, 14 Be, 15  In Fig. 2, I present the convergence of the ground-state energies with the maximum hypermomentum K max . In all cases, K max = 16 guarantees a good convergence. In this general overview, I also evaluate the importance of core excitations. In Fig. 2, the dotted lines represent the energies obtained by neglecting core excitations. This is done in Eq. I discuss now the properties of each nucleus by increasing mass, except for 11 Li which I present as last application since, by far, it is the most complicated. B. 14 Be as a 12 Be+n+n system Before considering a full diagonalization of the basis, I first display the energy curves, where a single value of the generator coordinate R is considered. The energy curves E Jπ i (R) are therefore obtained from the eigenvalue problem where N Jπ γK,γ ′ K ′ (R, R ′ ) is the overlap kernel. They provide a useful overview of the system. From the attractive or repulsive character, one can predict the existence of bound states (or of narrow resonances). In all cases, the three-body threshold is subtracted.
The 12 Be+n+n energy curves are displayed in Fig. 3. In positive parity, there is a minimum for J = 0 + and J = 2 + , which correspond to the 14 Be ground state and to the 2 + resonance, respectively. In negative parity, the 0 − and 2 − curves are repulsive. There is a shallow minimum for J = 1 − which might be associated with a broad resonance. A deeper analysis of such resonances, however, would require a specific formalism for continuum states [8], and is beyond the scope of the present work.
An interesting characteristic of the system is provided by the energy convergence as a function of the maximum R value, which I denote as R max . In Fig. 4, I show 14 Be energies obtained by increasing the N value in Eq. (16)  or, in other words, by increasing R max . Two different behaviours can be clearly observed. The 0 + and 2 + energies are almost stable above R max ≈ 10 fm, and present a plateau. The 2 + energy is 0.25 MeV, in nice agreement with experiment (0.28 ± 0.01 MeV [25]). This result emphasizes the importance of a microscopic theory; a non-microscopic three-body model, based on 12 Be+n and on n + n phenomenological interactions, does not reproduce the 0 + and 2 + energies simultaneously [26]. The 12 Be(g.s.)+n+n component is 87% in the ground state, and 67% in the 2 + resonance. This means that core excitations play a role, and may explain why non-microscopic models, which ignore core excitations, cannot predict the 2 + energy accurately. The other curves in Fig. 4 do not present a plateau, which means that no further narrow resonance can be expected. In particular, the 1 − curve is typical of a continuum state where, according to the variational principle, the minimum energy is zero. Another useful method to distinguish between narrow and broad resonances is to analyze the convergence of the r.m.s. radius r 2 as a function of R max . This is displayed in Fig. 5. The 0 + and 2 + low-lying states reach convergence near R max ≈ 12 fm. In contrast, the second eigenvalues present a diverging behaviour at large R max . This is expected since, strictly speaking, matrix elements involving continuum states diverge. Although I cannot specifically address continuum states, this technique allows a clear distinction between narrow and broad resonances. It is consistent with the energy convergence shown in Fig. 4, and confirms that r.m.s. radii in the continuum should be considered very carefully. The ground-state radii are presented in Table I, and compared with experiment. The neutron radius is in nice agreement with experiment. The experimental proton radius, however, is significantly larger than the predicted value. This is difficult to understand from a 12 Be+n+n model, where the 14 Be proton radius should be close to the core proton radius. Experimental values are partly model dependent, and a reanalysis of the 14 Be charge radius would be welcome.  [27]. b Ref. [28]. c Ref. [29]. d Ref. [30]. e Ref. [31]. f Ref. [32].
Notice that the radii are slightly sensitive to the oscillator parameter b. From the total matter radius √ < r 2 > m , one can define the expectation value of the hyperradius < R 2 > from where < r 2 1 > is the core radius. In the shell model, this quantity is proportional to b and takes the values < r 2 1 >= 49b 2 /24 for 12 Be, < r 2 1 >= 27b 2 /13 for 13 B, < r 2 1 >= 32b 2 /15 for 15 N and 15 O, and < r 2 1 >= 17b 2 /9 for 9 Li. In contrast < R 2 >, which is associated with the external neutrons, weakly depends on the oscillator parameter. For compact states, < R 2 > is small and the matter radius approximately varies linearly with b. For halo states, however, < R 2 > is the dominant term, and the matter radius is almost insensitive to the oscillator parameter.
C. 15 B as a 13 B+n+n system The 15 B ground state is known to be bound by 3.77 MeV. Although the spin assignment is no definite, there are strong indications for a spin J = 3/2 − . Two excited states at E x = 1.336 MeV and E x = 2.743 MeV have been reported in Ref. [33]. On the theoretical side, shellmodel [34] and Antisymmetrized Molecular Dynamics (AMD, see Ref. [35]) calculations have been performed. Non-microscopic calculations are unavailable until now, essentially due to the lack of reliable 13 B+n potentials.
In Fig. 6, I present the energy curves of 15 B. As expected, the lowest energy is obtained for J = 3/2 − . I may also expect other bound states, corresponding to minima in the energy curves. The GCM spectrum, including all generator coordinates, is presented in Fig. 7. The theoretical spectrum shown inf Fig. 7 is remarkably supported by experiment [33], although there is no spin assignment. The amount of 13 B(g.s.)+n+n component is shown for the GCM calculation. As suggested by Fig. 6, the role of core excitation is small in 15 B. As for 14 Be, the positiveparity energy curves are not completely repulsive, but do not support narrow resonances. A more reliable study of positive-parity states in 15 B would require introducing sd-shell components in the 13 B wave functions. This would considerably increase the computer times, and is not feasible at the moment.
The proton and matter radii, given in Table I are in reasonable agreement with experiment. The theoretical radius of the 13 B core is 2.31 fm, which means an increase of 0.33 fm for 15 B.
Let me briefly discuss the 15 Ne mirror nucleus, which has been addressed experimentally in Ref. [36] and theoretically in Ref. [37]. The ground state and first excited state have been observed in two-neutron knockout reactions from a 17 Ne beam at 2.5 and 4.4 MeV above the 13 O+p+p threshold. Both states are expected to be broad.
In the present model, the 15 B and 15 Ne mirror nuclei are studied in the same conditions. The energy curves are shown in Fig. 6(b). The 3/2 − and 5/2 − curves present minima which are associated with the ground and first excited states. These minima are close to the Coulomb barrier, and clearly suggest continuum states. According to the higher centrifugal barrier, however, the 5/2 − excited state could be narrower than the 3/2 − ground state. The 15 Ne spectrum is shown in Fig. 7. Although the GCM energies are obtained in the bound-state approximation, the results are close to the energies observed experimentally [36]. Microscopic calculations involving a 15 N or a 15 O core are relatively simple since only the ground state 1/2 − and the first excited state 3/2 − are present in the 0 ω shell model. A previous three-cluster microscopic calculation was performed in Ref. [38], where the main goal was to address the possible existence of a proton halo in 17 Ne. Figure 8 shows the energy curves, which present a minimum near R = 0 for negative parity. This suggests a compact structure for negative-parity states, in contrast with the other systems considered here.
The 17 N spectrum is displayed in Fig. 9. The general agreement with experiment is reasonable but, as in  Ref. [38] the ordering of the two first excited states is incorrect. For all states the 15 N(g.s.)+n+n component is dominant. No positive parity state is found, as suggested by the repulsive energy curves of Fig. 8.
The 17 Ne spectrum is studied with the same nucleonnucleon interaction; the Majorana parameter is the same as in 17 N. The GCM reproduces the binding energy very well (−1.05 MeV in the model, to be compared to the experimental value −0.93 MeV), which shows that the Coulomb shift is accurately described.
Electromagnetic transitions have been suggested to be a valuable tool to investigate the structure of 17 Ne [39,40]. The E2 transition probabilities have been studied by relativistic Coulomb excitation with a 17 Ne radioactive beam [41][42][43]. The B(E2) values, computed without any effective charge, are presented in Table II. The B(E2, 1/2 − → 5/2 − ) GCM value is in excellent agreement with the latest data of Marganiec et al. [43]. The experimental value of Chromik et al. [42] is larger, but is probably influenced by nuclear effects which have been dismissed in the analysis [43]. The transition to the 3/2 − state is also well reproduced, which shows that the GCM wave functions are reliable. The matter radius of the 17 N ground state, given in Table I, is an excellent agreement with experiment [29]. For 17 Ne, however, even if the binding energy is lower, the GCM does not support the large difference with 17 N. I rather confirm the conclusion of Ref. [38], that there is no evidence for a proton halo in 17 Ne.
FIG. 9. 17 N and 17 Ne energy spectra, compared with experiment [44]. Only negative-parity states are shown. Numbers at the left of GCM states correspond to the amount of 15 N(g.s.)+n+n configuration (in %).

E. 11 Li as a 9 Li+n+n system
The 11 Li nucleus has been studied in many experimental and theoretical works. The very low binding energy (S 2n = 0.367 MeV [24]) is responsible for a remarkable halo structure, as suggested by the large r.m.s. radius. The present calculation is made difficult because of the large number of Slater determinants (90) involved in the 9 Li core. A previous microscopic three-cluster study was performed in Ref. [45], but with a frozen triangular geometry. In Ref. [46], the authors describe the core in an α + t + n + n multicluster configuration. The 9 Li+n+n energy curves are displayed in Fig. 10. A minimum is obtained for J = 3/2 − and, to a lesser extend, for J = 1/2 − . The 3/2 − minimum corresponds to the 11 Li ground state. The role of core excitations is illustrated in Fig. 2. When core excitations are neglected, the state is unbound. The importance of core excitations was already pointed out in Ref. [46]. In contrast with the previous examples, where the neutron number of the core N = 8 corresponds to a closed shell, the various 9 Li+n+n configurations are not orthogonal to each other. Consequently, a ground-state component in 11 Li cannot be estimated. The proton, neutron and matter radii are presented in Table I. The calculation confirms the large enhancement of the matter radius, compared to the 9 Li core (the shell-model value is √ < r 2 > m ( 9 Li) = b 17/9 = 2.20 fm). The GCM value is, however, slightly smaller than experiment. For the proton radius, the experimental 11 Li value (2.467 fm) is significantly larger than the 9 Li value (2.217 fm), which suggests that the neutron halo of 11 Li affects the core [31,47]. In the present model, the 9 Li proton radius is < r 2 > p ( 9 Li) = b 4/3 = 1.84 fm. The difference between 11 Li and 9 Li is therefore 0.13 fm, which is smaller than experiment (0.25 fm) (see a detailed discussion in Refs. [31,47]).

IV. CONCLUSION
The main goal of the present work is to extend the hyperspherical formalism to microscopic three-cluster models. The hyperradial functions are expanded over a Gaussian basis using the Generator Coordinate method. Then the wave functions are expressed in terms of projected Slater determinants. This approach is well adapted to multicluster systems. With a single generator coordinate, it provides an accurate description of the wave functions, even at long distances. The calculations of the matrix elements, however, require very long computer times, owing to the seven-dimension integrals necessary for the angular momentum projection. The extension to p-shell cores raises additional difficulties due to (i) the quadruple sums involved in the matrix elements of two-body interactions; (ii) the presence of several Slater determinants; (iii) the introduction of core excited states. The calculations are made possible thanks to an efficient parallelization of the computer code.
The model has been applied to some exotic light nuclei: 14 Be, 15 B, 17 Ne and 11 Li. In all cases, the only parameter (the admixture parameter M involved in the Volkov nucleon-nucleon interaction) is adjusted on the ground-state energy. For 14 Be, the 2 + excitation energy is well reproduced, in contrast with non-microscopic models. The GCM spectrum of 15 B is in nice agreement with the experimental energies. An exploratory study of the 15 Ne mirror system, which is unbound, is consistent with the experimental energies, and suggests that the first excited state is narrower than the ground state. The 17 N spectrum presents many states which, in general, are fairly well reproduced by the GCM. The ground state is confirmed to be 1/2 − , but the order of the first two excited states in incorrect, as in Ref. [38]. In the 17 Ne mirror nucleus, the ground-state energy, as well as E2 transition probabilities are in good agreement with experiment. For 11 Li, calculations are extremely long, owing to the 90 Slater determinants involved in the core. The neutron and matter radii are explained fairly well, but the proton radius is smaller than the experimental value.
The present model could be extended to deal with three-body continuum states [8]. On the other hand, the wave functions could be used as an input to determine scattering cross sections, as it was done for 6 He for example [9].