nuclei are quantum many - body systems with both bound and unbound states . a realistic _ ab initio _ description of light nuclei with predictive power must have the capability to describe both classes of states within a unified framework . over the past decade
, significant progress has been made in our understanding of the properties of the bound states of light nuclei starting from realistic nucleon - nucleon ( @xmath11 ) interactions , see e.g. ref .
@xcite and references therein , and more recently also from @xmath11 plus three - nucleon ( @xmath12 ) interactions  @xcite .
the solution of the nuclear many - body problem becomes more complex when scattering or nuclear reactions are considered . for @xmath13 and 4 nucleon systems , the faddeev  @xcite and faddeev - yakubovsky  @xcite as well as the hyperspherical harmonics ( hh ) @xcite or the alt ,
grassberger and sandhas ( ags ) @xcite methods are applicable and successful . however ,
_ ab initio _
calculations for unbound states and scattering processes involving more than four nucleons in total are quite challenging .
the first _ ab initio _ many - body neutron-@xmath0he scattering calculations were performed within the green s function monte carlo method using the argonne @xmath11 potential and the illinois @xmath12 interaction  @xcite .
also , resonances in he isotopes were investigated within the coupled - cluster method using the gamow basis  @xcite .    in a new development ,
we have recently combined the _ ab initio _ no - core shell model ( ncsm )  @xcite and the resonating - group method ( rgm )  @xcite , into a new many - body approach  @xcite ( _ ab initio _
ncsm / rgm ) capable of treating bound and scattering states of light nuclei in a unified formalism , starting from fundamental inter - nucleon interactions .
the ncsm is an _ ab initio _ approach to the microscopic calculation of ground and low - lying excited states of light nuclei with realistic two- and , in general , three - nucleon forces .
the rgm is a microscopic cluster technique based on the use of @xmath14-nucleon hamiltonians , with fully anti - symmetric many - body wave functions built assuming that the nucleons are grouped into clusters .
although most of its applications are based on the use of binary - cluster wave functions , the rgm can be formulated for three ( and , in principle , even more ) clusters in relative motion  @xcite . the use of the harmonic oscillator ( ho ) basis in the ncsm results in an incorrect description of the wave - function asymptotic and a lack of coupling to the continuum . by combining the ncsm with the rgm , we complement the ability of the rgm to deal with scattering and reactions with the use of realistic interactions , and a consistent _ ab initio _ description of the nucleon clusters , achieved via the ncsm .
presently the ncsm / rgm approach has been formulated for processes involving binary - cluster systems only .
however , extensions of the approach to include three - body cluster channels are feasible , also in view of recent developments on the treatment of both three - body bound and continuum states ( see , e.g. , refs .
@xcite ) . as described in detail in refs .
@xcite , the _ ab initio _
ncsm / rgm approach has been already applied to study the @xmath15-@xmath16h , @xmath15-@xmath17he , @xmath15-@xmath18be , and @xmath19-@xmath20he scattering processes , and address the parity inversion of the @xmath21be ground state , using realistic @xmath11 potentials . in that work
, we demonstrated convergence of the approach with increasing basis size in the case of the @xmath22 and @xmath23 scattering .
the @xmath15-@xmath18be calculations were , on the other hand , perfomed only in a limited basis due to the computational complexity of the ncsm calculations of the @xmath24be eigenstates .
it is the purpose of the present paper to expand the applicability of the ncsm / rgm beyond the lightest nuclei by using sufficiently large @xmath25 ho excitations to guarantee convergence of the calculation with the ho basis expansion of both the cluster wave functions and the localized rgm integration kernels .
the use of large @xmath25 values is now feasible due to the recent introduction of the importance truncated ( it ) ncsm scheme  @xcite .
it turns out that many of the basis states used in the ncsm calculations are irrelevant for the description of any particular eigenstate , e.g. , the ground state or a set of low - lying states . therefore ,
if one were able to identify the important basis states beforehand , one could reduce the dimension of the matrix eigenvalue problem without losing predictive power .
this can be done using an importance truncation scheme based on many - body perturbation theory  @xcite .
we make use of the it ncsm wave functions for the cluster eigenstates , in particular the eigenstates of the target nucleus of the binary nucleon - nucleus system , and calculate the one- and two - body densities that are then used to obtain the ncsm / rgm integration kernels .
we benchmark the it approach in basis sizes accessible by the full calculation and apply it within still larger basis sizes until convergence is reached for target nuclei as heavy as @xmath2c or @xmath3o . in this study
, we employ a similarity renormalization group ( srg )  @xcite evolved chiral n@xmath26lo @xmath11 potential  @xcite ( srg - n@xmath26lo ) that is soft enough to allow us reach convergence within about @xmath27 ho excitations in the basis expansion .    in sect .
[ formalism ] , we briefly overview the ncsm / rgm formalism and present for the first time the it ncsm scheme that includes both ground and low - lying excited states in the set of reference states .
next , we present scattering calculation results for the @xmath28-@xmath0he and @xmath29-@xmath0he systems in sect .
in particular , we compare the calculated phase shifts to an r - matrix analysis of experimental data and , further , calculated differential cross sections and analyzing powers in the energy range  6 - 19 mev to the corresponding experimental data .
neutron elastic and inelastic scattering on @xmath1li and proton elastic and inelastic scattering on @xmath1be are investigated in sect .
we present phase shifts , cross sections and scattering lengths .
we predict resonances in @xmath7li and @xmath7be that have not been clearly identified in experiments yet . in sect .
[ n12c ] , we discuss nucleon-@xmath2c results for both bound and unbound states of @xmath30c and @xmath30n , obtained including two @xmath2c bound states , the ground and the first @xmath31 state , in the ncsm / rgm coupled - channel calculations . in sect .
[ n16o ] , we present results for the nucleon-@xmath3o system . in the ncsm / rgm coupled - channel calculations , we take into account the @xmath3o ground state and up to the lowest three @xmath3o negative - parity states .
conclusions are given in sect .
[ conclusions ] .
the _ ab inito _
ncsm / rgm approach was introduced in ref .
@xcite with details of the formalism given in ref .  @xcite . here
we give a brief overview of the main points .    in the present paper ,
we limit ourselves to a two - cluster rgm , which is based on binary - cluster channel states of total angular momentum @xmath32 , parity @xmath33 , and isospin @xmath34 , @xmath35^{(j^\pi t)}\,\frac{\delta(r - r_{a - a , a})}{rr_{a - a , a}}\,.\label{basis}\end{aligned}\ ] ] in the above expression , @xmath36 and @xmath37 are the internal ( antisymmetric ) wave functions of the first and second cluster , containing @xmath38 and @xmath39 nucleons ( @xmath40 ) , respectively .
they are characterized by angular momentum quantum numbers @xmath41 and @xmath42 coupled together to form channel spin @xmath43 . for their parity , isospin and additional quantum numbers we use , respectively , the notations @xmath44 , and @xmath45 , with @xmath46 .
the cluster centers of mass are separated by the relative coordinate @xmath47 where @xmath48 are the @xmath14 single - particle coordinates .
the channel states  ( [ basis ] ) have relative angular momentum @xmath49 .
it is convenient to group all relevant quantum numbers into a cumulative index @xmath50 .
the former basis states can be used to expand the many - body wave function according to @xmath51 as the basis states  ( [ basis ] ) are not anti - symmetric under exchange of nucleons belonging to different clusters , in order to preserve the pauli principle one has to introduce the appropriate inter - cluster anti - symmetrizer , schematically @xmath52 where the sum runs over all possible permutations @xmath53 that can be carried out among nucleons pertaining to different clusters , and @xmath29 is the number of interchanges characterizing them .
the coefficients of the expansion  ( [ trial ] ) are the relative - motion wave functions @xmath54 , which represent the only unknowns of the problem . to determine them one has to solve the non - local integro - differential coupled - channel equations @xmath55
\frac{g^{j^\pi t}_\nu(r)}{r } = 0\,,\label{rgmeq}\ ] ] where the two integration kernels , the hamiltonian kernel , @xmath56 and the norm kernel , @xmath57 contain all the nuclear structure and anti - symmetrization properties of the problem . in particular , the non - locality of the kernels is a direct consequence of the exchanges of nucleons between the clusters .
we have used the notation @xmath58 and @xmath59 to denote the total energy in the center - of - mass frame , and the intrinsic @xmath14-nucleon microscopic hamiltonian , respectively .
the formalism presented above is combined with the _ ab initio _ ncsm in two steps :    first , we note that the hamiltonian can be written as @xmath60 where @xmath61 and @xmath62 are the ( @xmath38)- and @xmath39-nucleon intrinsic hamiltonians , respectively , @xmath63 is the relative kinetic energy and @xmath64 is the sum of all interactions between nucleons belonging to different clusters after subtraction of the average coulomb interaction between them , explicitly singled out in the term @xmath65 ( @xmath66 and @xmath67 being the charge numbers of the clusters in channel @xmath68 ) .
we use identical realistic potentials in both the cluster s hamiltonians and inter - cluster interaction @xmath64 .
accordingly , @xmath36 and @xmath37 are obtained by diagonalizing @xmath61 and @xmath62 , respectively , in the model space spanned by the ncsm @xmath69 ho basis .
note that in the present paper we use soft srg evolved @xmath11 potentials .
therefore , there is no need to derive any further effective interaction tailored to the model space truncation as with these soft interactions our calculations converge in the model spaces we are able to reach .
second , we replace the delta functions in the localized parts of the hamiltonian  ( [ h - kernel ] ) and the norm  ( [ n - kernel ] ) kernels with their representation in the ho model space .
we use identical ho frequency as for the cluster eigenstate wave functions and a consistent model space size ( @xmath70 ) .
we emphasize that this replacement is performed only for the localized parts of the kernels .
the diagonal parts coming from the identity operator in the antisymmetrizers , the kinetic term and the average coulomb potential are treated exactly .    in this paper
, we apply the ncsm / rgm formalism in the single - nucleon projectile basis , i.e. , for binary - cluster channel states  ( [ basis ] ) with @xmath71 ( with channel index @xmath72 ) . as an illustration , let s discuss in more detail the norm kernel that is rather simple in this basis : @xmath73 we can easily recognize a direct term , in which initial and final state are identical ( corresponding to diagram @xmath74 of fig .  [ diagram - norm - pot ] ) , and a many - body correction due to the exchange part of the inter - cluster anti - symmetrizer ( corresponding to diagram @xmath75 of fig .
[ diagram - norm - pot ] ) .
we note that in calculating the matrix elements of the exchange operator @xmath76 we replaced the delta function of eq .
( [ basis ] ) with its representation in the ho model space as discussed above .
this is appropriate as the transposition @xmath76 operator acting on the traget wave function is short - to - medium range . on the contrary ,
the @xmath77-function coming from the identity is treated exactly .
the presence of the inter - cluster anti - symmetrizer affects also the hamiltonian kernel , and , in particular , the matrix elements of the interaction . for a @xmath11 potential one obtains a direct term involving interaction and exchange of two nucleons only ( see diagrams ( @xmath78 ) and ( @xmath79 ) of fig .  [ diagram - norm - pot ] ) , and an exchange term involving three - nucleons .
diagram ( @xmath80 ) of fig .
[ diagram - norm - pot ] describes this latter term , in which the last nucleon is exchanged with one of the nucleons of the first cluster , and interacts with yet another nucleon . for more details on the integration kernels in the single - nucleon projectile basis
we refer the readers to ref .
@xcite .
being translationally - invariant quantities , the norm and hamiltonian kernels can be  naturally " derived working within the ncsm jacobi - coordinate basis . however , by introducing slater - determinant channel states of the type @xmath81^{(j^\pi t ) } r_{n\ell}(r^{(a)}_{\rm c.m.})\ , , \label{sd - basis}\end{aligned}\ ] ] in which the eigenstates of the @xmath82-nucleon fragment are obtained in the sd basis ( while the second cluster is still a ncsm jacobi - coordinate eigenstate ) , it can be easily demonstrated that translationally invariant matrix elements can be extracted from those calculated in the sd basis of eq .
( [ sd - basis ] ) by inverting the following expression : @xmath83 here @xmath84 represents any scalar and parity - conserving translational - invariant operator ( @xmath85 , @xmath86 , etc . ) .
we exploited both jacobi - coordinate and sd channel states to verify our results .
the use of the sd basis is computationally advantageous and allows us to explore reactions involving @xmath29-shell nuclei , as done in the present work . in order to calculate the parts of the integration kernels depicted in fig .
[ diagram - norm - pot ] ( b ) , ( c ) and ( d ) , all information that we need from the sd basis calculation are one - body densities of the target eigenstates .
for the ( e ) part of the integration kernel in fig .
[ diagram - norm - pot ] , we need two - body densities of the target eigenstates obtained in the sd basis .    due to the presence of the norm kernel @xmath87 , eq .
( [ rgmeq ] ) does not represent a system of multichannel schrdinger equations , and @xmath54 do not represent schrdinger wave functions .
the short - range non - orthogonality , induced by the non - identical permutations in the inter - cluster anti - symmetrizers , can be removed by introducing normalized schrdinger wave functions @xmath88 where @xmath89 is the square root of the norm kernel , and applying the inverse - square root of the norm kernel , @xmath90 , to both left and right - hand side of the square brackets in eq .
( [ rgmeq ] ) .
this procedure , explained in more detail in ref .
@xcite , leads to a system of multichannel schrdinger equations @xmath91\frac{\chi^{j^\pi t}_{\nu } ( r)}{r } \nonumber\\[2 mm ] & & + \sum_{\nu^\prime}\int dr^\prime\,r^{\prime\,2 } \,w^{j^\pi t}_{\nu \nu^\prime}(r , r^\prime)\,\frac{\chi^{j^\pi t}_{\nu^\prime}(r^\prime)}{r^\prime } = 0,\label{r - matrix - eq}\end{aligned}\ ] ] where @xmath92 is the energy eigenvalue of the @xmath93-th cluster ( @xmath46 ) , and @xmath94 is the overall non - local potential between the two clusters , which depends on the channel of relative motion , while it does not depend on the energy .
these are the equations that we finally solve to obtain both our scattering and bound - state results .      the primary limitation for the range of applicability of the ncsm in terms of particle number @xmath14 and model spaces size @xmath95 results from the factorial growth of the dimension of the @xmath96 space .
except for light isotopes , it is hardly possible to obtain a converged result using a bare hamiltonian within the @xmath96 spaces that are computationally tractable .    at this point
the importance truncation offers a solution .
the importance truncation in connection with the ncsm was introduced in ref .
@xcite and discussed in detail in ref .
@xcite . in the following
we summarize a few key features of the it - ncsm and generalize the approach to the simultaneous description of excited states .    the motivation for the importance truncation results from the observation that the expansion of any particular eigenstate of the hamiltonian in a full @xmath97-scheme ncsm space typically contains a large number of basis states with extremely small or vanishing amplitudes .
the amplitudes define an adaptive truncation criterion , which takes into account the properties of the hamiltonian and the structure of the eigenstate under consideration . if those amplitudes were known  at least approximately  before actually solving the eigenvalue problem , one could reduce the model space to the most relevant basis states by imposing a threshold on the amplitude .
the amplitude of a particular basis state @xmath98 in the expansion of a specific eigenstate can be estimated using first - order multiconfigurational perturbation theory . in order to set up a perturbation series we need an initial approximation of the target state , the so - called reference state @xmath99 . in practice
this reference state will be a superposition of basis states @xmath100 from a reference space @xmath101 : @xmath102 the reference state and the amplitudes @xmath103 are typically extracted from a previous ncsm calculation . based on @xmath99 as unperturbed state
, we can evaluate the first - order perturbative correction to the target state resulting from basis states @xmath104 .
their first - order amplitude defines the so - called importance measure @xmath105 the energy denominator @xmath106 in a mller - plesset - type partitioning is simply given by the unperturbed harmonic - oscillator excitation energy of the basis state @xmath98 ( see ref .
@xcite for details ) .    imposing an importance threshold @xmath107
, we construct an importance truncated model space including all basis states with importance measure @xmath108 .
since the importance measure is zero for all basis states that differ from all of the states in @xmath101 by more than a two - particle - two - hole excitation , we have to embed the construction of the importance truncated space into an iterative update cycle . after constructing the importance truncated space and solving the eigenvalue problem in that space
, we obtain an improved approximation for the target state that defines a reference state for the next iteration . in order to accelerate the evaluation of the importance measure
, we typically do not use the complete eigenstate as new reference state , but project it onto a reference space spanned by the basis states with @xmath109 , where @xmath110 are the coefficients resulting from the solution of the eigenvalue problem .
the second threshold @xmath111 will be chosen sufficiently small so as not to affect the results for a given @xmath107 threshold .
simple iterative update schemes can be devised for any type of full model spaces , as discussed in refs .
specifically for the @xmath96 space of the ncsm , however , there is an efficient sequential update scheme leading to the it - ncsm(seq ) approach .
it is based on the fact that all states of an @xmath112 space can be generated from the basis states of an @xmath96 space using two - particle - two - hole excitations at most .
thus a single importance update starting from a reference state in an @xmath96 space gives access to all relevant states in an @xmath112 space . making use of this property , in the it - ncsm(seq )
we start with a full ncsm calculation in , e.g. , @xmath113 and use this eigenstate after applying the @xmath111 threshold as reference state for constructing the importance truncated @xmath114 space . after solving the eigenvalue problem for this importance truncated
@xmath114 space we use the resulting eigenstate as reference state to construct the @xmath115 space , and so on .
thus only one importance update is required for each value of @xmath95 , which makes this scheme very efficient computationally . moreover , in the limit of vanishing thresholds , @xmath116 , this scheme recovers the full @xmath96 space at each step of the sequence , i.e. , the it - ncsm(seq ) would recover the full ncsm result .
based on this limiting property , we can obtain a numerical approximation to the full ncsm result by extrapolating the it - ncsm(seq ) observables obtained for a set of different importance thresholds @xmath107 ( and in principle also @xmath111 ) to @xmath117 . through this extrapolation ,
the contribution of discarded basis states , i.e. those with importance measures @xmath118 below the smallest threshold considered , is effectively recovered . because the control parameter @xmath107 is tied to the physical structure of the eigenstate
, we observe a smooth threshold dependence for all observables , which allows for a robust threshold extrapolation . in the case of the energy we can improve the quality of the extrapolation further by considering a perturbative second - order estimate for the energy of the excluded basis states .
while setting up the importance truncated space , all second - order energy contributions @xmath119 for the discarded states with @xmath120 are summed up to provide a correction @xmath121 to the energy eigenvalue . by construction
this correction goes to zero in the limit @xmath117 .
we use this additional information for a constrained simultaneous extrapolation of the energy to vanishing threshold with and without perturbative correction for the excluded states as described in detail in ref .
@xcite .
the whole concept can be generalized to the description of excited states .
for the present application in connection with the ncsm / rgm , we aim at an importance truncated model space that is equally well suited for the description of the lowest @xmath122 eigenstates of the hamiltonian for given parity and total angular momentum projection . instead of using a single reference state ,
we employ different reference states @xmath123 , with @xmath124 , for each of the @xmath122 target states . for each reference state
we define a separate importance measure @xmath125 following eq . .
a basis state @xmath126 is included in the importance truncated space if at least one of the importance measures @xmath127 is above the threshold @xmath107 , i.e. , if it is relevant for the description of at least one of the @xmath122 target states it will be included .
because the different eigenstates are typically dominated by different basis states , the dimension of the importance truncated space grows linearly with @xmath122 .    in the it - ncsm(seq ) scheme
we start with a full ncsm calculation in @xmath113 and use the lowest @xmath122 eigenstates as initial reference states @xmath123 . based on the corresponding importance measures
@xmath125 the importance truncated @xmath114 space is constructed and the lowest @xmath122 eigenvectors in this space serve as new reference states ( after application of the @xmath111 threshold ) for the construction of the @xmath115 space , and so on . from a sequence of it - ncsm(seq ) calculations we obtain a set of @xmath122 eigenvectors for each value of @xmath95 which can be used to evaluate other observables .    by default
we compute the expectation values of @xmath128 and @xmath129 as well as the expectation values of @xmath130 and @xmath131 .
indeed , since we use an importance truncated space in the @xmath97-scheme without explicit angular momentum projection , the eigenstates are not guaranteed to have good angular momentum and isospin .
we therefore monitor the expectation values of @xmath128 and @xmath129 and find values which typically differ by less then @xmath132 from the exact quantum numbers . as in the full ncsm we separate spurious center - of - mass ( cm ) excitations from the physical spectrum by adding a lawson term @xmath133 to the translationally invariant intrinsic hamiltonian @xmath130 ( with the typical choice @xmath134 ) .
the use of this modified hamiltonian provides at the same time a diagnostic for potential cm contaminations of the intrinsic states induced by the importance truncation . as discussed in refs .
@xcite , the independence of the intrinsic energies @xmath135 on @xmath136 and the smallness of @xmath137 demonstate that the it - ncsm(seq ) solutions are free of cm contaminations .    eventually , the wave functions obtained in the it - ncsm(seq ) together with the threshold extrapolated intrinsic energies form the input for the ncsm / rgm calculations discussed in the following .
the purpose of the nucleon-@xmath0he calculations presented in this paper is two - fold .
first , we want to check the predictive power of the srg evolved chiral interaction in the @xmath23 system , where a lot of experimental scattering data exist and where our calculations can be easily converged with respect to the size of the basis expansion .
second , we want to benchmark the importance truncation scheme with the full - space calculations all the way up to very large @xmath69 spaces .    the first _ ab initio _
@xmath23 scattering calculations was reported in ref .
@xcite . the @xmath28-@xmath10 low - lying @xmath138 and @xmath139 @xmath53-wave resonances as well as the @xmath140 @xmath141-wave non - resonant scattering below 5 mev c.m .
energy were obtained using the av18 @xmath11 potential with and without the three - nucleon force , chosen to be either the urbana ix or the illinois-2 model .
the results of these green s function monte carlo ( gfmc ) calculations revealed sensitivity to the inter - nucleon interaction , and in particular to the strength of the spin - orbit force .    soon after , the development of the _ ab initio _ ncsm / rgm approach allowed us to calculate both @xmath28- and ( for the first time ) @xmath29-@xmath10 scattering phase shifts for energies up to the inelastic threshold  @xcite , using several realistic @xmath11 potentials , including the chiral n@xmath26lo  @xcite , the @xmath142  @xcite and the cd - bonn  @xcite @xmath11 potentials .
nucleon-@xmath10 scattering provides one of the best - case scenarios for the application of the ncsm / rgm approach .
this process is characterized by a single open channel up to the @xmath143h threshold , which is fairly high in energy .
in addition , the low - lying resonances of the @xmath0he nucleus are narrow enough to be reasonably reproduced diagonalizing the four - body hamiltonian in the ncsm model space . in the present work we include the first excited state of @xmath0he , the @xmath144 state , as a closed channel in our ncsm / rgm basis space .
we performed extensive nucleon-@xmath0he calculations with the srg - n@xmath26lo @xmath11 potential with a cutoff of 2.02 @xmath145 to check convergence of our ncsm / rgm calculations . in fig .
[ fig : nhe4_phaseconv ] , we present @xmath28-@xmath0he phase shift results for the @xmath141- and @xmath53-waves obtained using an ho basis expansion up to @xmath146 for for the localized parts of the ncsm / rgm integration kernels and for the @xmath0he ground- and the first - excited @xmath144 wave functions ( since these states have positive parity , the @xmath147 expansion is in fact used for the @xmath0he eigenstates ) .
as seen in the figure , the phase - shift convergence is excellent .
in particular , the @xmath146 and the @xmath148 curves lie on top of each other .
the convergence rate demonstrated here is quite similar to that obtained using the @xmath149 @xmath11 potentatial in our erlier study ( compare the present fig .
[ fig : nhe4_phaseconv ] to the left panel of fig . 13 in ref .
@xcite ) .
-@xmath0he phase shifts on the size of the ho basis expansion of the @xmath0he wave functions and the localized parts of the integration kernels .
the @xmath0he ground state and the first @xmath144 excited states were included .
the srg - n@xmath26lo @xmath11 potential with a cutoff of 2.02 @xmath145 and the ho frequency @xmath150 mev were used . ]      as shown in the previous subsection , for the @xmath23 system we are able to reach complete convergence with @xmath0he wave functions obtained within full , non - truncated , ncsm calculations .
we can , therefore , test the performace of the it - ncsm scheme in this system all the way up to very large @xmath70 values and see how well the it - ncsm scheme reproduces the completetly converged results .
it should be noted that for heavier @xmath151 and @xmath9 systems investigated later , full , non - truncated ncsm calculations for the @xmath152 ( @xmath153 ) target nuclei are feasible only up to @xmath154 ( @xmath155 ) .
it is , therefore , desirable and important to benchmark the it - ncsm calculations in a lighter system like @xmath23 in @xmath156 calculations .    in fig .
[ fig : nhe4_phase_full_it ] , we compare @xmath28-@xmath0he phase shifts calculated within the ncsm / rgm with @xmath0he wave functions obtained in a full @xmath157 ncsm calculation and those obtained using @xmath0he wave functions obtained within an @xmath157 it - ncsm calculation .
the agreement of the two sets of phase shifts is excellent .
it should be noted that the dimension of the full @xmath157 @xmath0he ncsm basis is 6344119 .
the dimension of the it - ncsm basis used here to calculate the @xmath0he wave functions was just 992578 , more than a factor of six smaller .
truncation parameters @xmath158 and @xmath159 were used .
the ground state energy from the full ncsm calculation is @xmath160 mev .
the @xmath161 extrapolated ground state energy from the it - ncsm calculation is @xmath162 mev with a difference from the full result less than 10 kev .
the excited @xmath144 energy obtained in the full ncsm calculation was 21.58 mev .
the corresponding extrapolated it - ncsm result was 21.4(1 ) mev .
the slightly lower accuracy of the excited state reproduction in the it - ncsm calculation is manifested in a very small deviation of the @xmath141-wave phase shift at energies above 12 mev ( less than 1 degree at 16 mev ) .
it should be noted that the excited @xmath144 state is not bound .
consequently , it is challenging to reproduce the excited state as well as the ground state in a importance - truncated calculation .
it should be also pointed out that unlike for the energies , no phase shift extrapolation was performed .
the needed one- and two - body densities were calculated from the wave functions obtained in the it - ncsm calculation with the truncation parameters described above .
the excellent agreement of the full and the it - ncsm phase shifts demonstrates that no extrapolation was actually necessary .
obviously , we can check the dependence of observables like phase shifts on the @xmath107 and @xmath111 and perform an extrapolation to vanishing values of these parameters if needed .
-@xmath0he @xmath141- and @xmath53-wave phase shifts .
results obtained with @xmath0he wave functions from full ncsm ( solid lines ) and it - ncsm ( dashed lines ) calculations are compared .
the srg - n@xmath26lo @xmath11 potential with a cutoff of 2.02 @xmath145 , the @xmath146 basis space and the ho frequency @xmath150 mev were used .
see text for details on the it - ncsm calculation . ]
he ( left panels ) and @xmath163he ( right panels ) compared to the r - matrix analysis of experimental data  @xcite .
the ncsm / rgm calculations that included the @xmath0he ground state and the @xmath144 excited state were done using the srg - n@xmath26lo @xmath11 potential with a cutoff of 2.02 @xmath145 . the ho frequency @xmath150 mev and @xmath146 basis space were employed . ]
our calculated @xmath28-@xmath0he and @xmath29-@xmath0he phase shifts are compared to those obtained from an @xmath164-matrix analysis of @xmath165he experimental data  @xcite in fig .
[ fig : n4he_phase ] .
the agreement is quite reasonable for the @xmath141-wave , @xmath166-wave and @xmath167-wave .
the @xmath168 resonance is positioned at higher energy in the calculation and the corresponding phase shifts are underestimated with respect to the @xmath164-matrix results , although the disagreement becomes less and less pronounced starting at about 8 mev .
while the inclusion of negative - parity excited states of the @xmath169particle would likely increase somewhat the @xmath168 phase shifts  @xcite , the observed difference is largely due to a reduction in spin - orbit strength caused by the neglect of the three - nucleon interaction in our calculations .
the importance of the three - nucleon force in reproducing the @xmath164-matrix @xmath168 phase shifts was demonstrated in the gfmc @xmath28-@xmath0he calculations of ref .
overall , the present results obtained with the srg - n@xmath26lo @xmath11 interaction agree better with experiment than our earlier calculations  @xcite with the @xmath149 , n@xmath26lo and cd - bonn @xmath11 potentials .
the only exception is the @xmath141-wave phase shift which is best described using the cd - bonn @xmath11 potential .
the larger spin - orbit strength of the employed srg - n@xmath26lo potential with respect to n@xmath26lo itself is the likely responsible for the improved agreement .
as our calculated phase shifts agree with the experimental ones reasonably well above the center - of - mass energy of 8 mev , we expect a similar behavior for cross section and analysing power in that energy range .
this is indeed the case as shown in fig .
[ fig : n4he_17mevxsay ] , where the calculated differential cross section and analyzing power are compared to experimental data from karlsruhe  @xcite with polarized neutrons of @xmath17017 mev laboratory energy . for the cross section experimental data
see also references in  @xcite .
the cross section is reproduced remarkably well at all angles and the analysing power is in reasonable agreement with the data , particularly at backward angles .
the same quality of agreement can be found for all energies far from the low - lying resonances , as shown in the right panel of fig .
[ fig : n4he_17mevxsay ] for the analysing power at @xmath171 mev and 19 mev
.    a better display of the dependence of our calculated cross section and analysing power upon the incident nucleon energy is provided by fig .
[ fig : p4he_1 ] , where the @xmath163he results for these observables are compared to the data of ref .
@xcite at the proton laboratory energies of @xmath172 , 7.89 , 9.89 , and 11.99 mev . as expected from the behavior of the phase shifts described earlier , for energies relatively close to the resonance region we find a rather poor agreement with experiment , particularly noticeable in the analysing power overall and in the cross section at backward angles .
however , starting at about 10 mev , the agreement improves substantially and data are once again reproduced in a quite satisfactory way at higher energies , as shown in fig .
[ fig : p4he_2 ] , where the ncsm / rgm @xmath163he results are compared to various experimental data sets  @xcite in the energy range @xmath173 mev .
the @xmath1be(@xmath29,@xmath175)@xmath7b capture reaction plays a very important role in nuclear astrophysics as it serves as an input for understanding the solar neutrino flux  @xcite . while the experimental determination of the neutrino flux from @xmath7b has an accuracy of about 9%  @xcite , the theoretical predictions have uncertainties of the order of 20%  @xcite .
the theoretical neutrino flux depends on the @xmath1be(@xmath29,@xmath175)@xmath7b s - factor .
significant experimental and theoretical effort has been devoted to studying this reaction .
the s - factor extrapolation to astrophysically relevant energies depends among other things on the scattering lengths of the proton scattering on @xmath1be .
experimental determination of these lengths was performed recently  @xcite with precision of the order of 30% .
the proton-@xmath1be elastic scattering was also investigated in ref .
@xcite . to benchmark the theoretical calculations used for s - factor extrapolations , an investigation of the mirror capture reaction , @xmath1li(@xmath28,@xmath175)@xmath7li , as well as
the @xmath28+@xmath1li scattering is important .
for example , the @xmath28+@xmath1li scattering lengths are known with a higher accuracy  @xcite .
li ground - state and the @xmath139 and @xmath176 excited state energy dependence on the model - space size @xmath70 , obtained within the importance - truncated ncsm ( solid lines ) , using the srg - n@xmath26lo @xmath11 potential with a cutoff of 2.02 @xmath145 .
the ho frequency @xmath150 mev was employed .
the full - space ncsm results are shown by dashed lines . ]
the first applications of the ncsm approach to the description of the @xmath1be(@xmath29,@xmath175)@xmath7b capture reaction  @xcite required a phenomenological correction of the asympotic behavior of the overlap functions and , further , the scattering @xmath29+@xmath1be wave function was calculated from a phenomenological potential model .
the present investigation within the _ ab initio _
ncsm / rgm approach paves the way for a complete first principles calculation of this capture reaction . here , we limit ourselves to scattering calculations and postpone the capture reaction calculations to a forthcoming paper .    our current limit on the unrestricted ncsm calculations for @xmath1li and @xmath1be is @xmath154 . to improve the convergence of our scattering calculations , we utilize wave functions obtained within the it - ncsm . in that scheme ,
we are able to reach @xmath177 model spaces and calculate both ground as well as low - lying excited states .
this is demonstrated in fig .
[ fig : li7_itncsm ] . with the srg - n@xmath26lo @xmath11 potential with
@xmath178 @xmath145 employed in the present study we reach convergence already around @xmath179 .
also , as seen in the figure , the aggreement between the unrestricted ncsm and the it - ncsm is perfect up to the highest accessible unrestricted space , @xmath154 .
-wave diagonal phase shifts of the @xmath28-@xmath1li elastic scattering ( top panel ) , elastic @xmath1li(@xmath28,@xmath28)@xmath1li cross section ( middle panel ) , and inelastic @xmath1li(@xmath28,@xmath180)@xmath1li(1/2@xmath181 ) cross section ( bottom panel ) .
the ncsm / rgm calculation that included the @xmath1li ground state and the @xmath139 and @xmath176 excited states were done using the srg - n@xmath26lo @xmath11 potential with a cutoff of 2.02 @xmath145 .
wave functions from it - ncsm calculations in the @xmath182 basis and the ho frequency of @xmath150 mev were employed .
experimental data are from ref .  @xcite . ]
-wave diagonal phase shifts of the @xmath28-@xmath1li elastic scattering ( top panel ) , elastic @xmath1li(@xmath28,@xmath28)@xmath1li cross section ( middle panel ) , and inelastic @xmath1li(@xmath28,@xmath180)@xmath1li(1/2@xmath181 ) cross section ( bottom panel ) .
the ncsm / rgm calculation that included the @xmath1li ground state and the @xmath139 and @xmath176 excited states were done using the srg - n@xmath26lo @xmath11 potential with a cutoff of 2.02 @xmath145 .
wave functions from it - ncsm calculations in the @xmath182 basis and the ho frequency of @xmath150 mev were employed .
experimental data are from ref .  @xcite . ]
-wave diagonal phase shifts of the @xmath28-@xmath1li elastic scattering ( top panel ) , elastic @xmath1li(@xmath28,@xmath28)@xmath1li cross section ( middle panel ) , and inelastic @xmath1li(@xmath28,@xmath180)@xmath1li(1/2@xmath181 ) cross section ( bottom panel ) .
the ncsm / rgm calculation that included the @xmath1li ground state and the @xmath139 and @xmath176 excited states were done using the srg - n@xmath26lo @xmath11 potential with a cutoff of 2.02 @xmath145 .
wave functions from it - ncsm calculations in the @xmath182 basis and the ho frequency of @xmath150 mev were employed .
experimental data are from ref .  @xcite . ]
the ncsm / rgm coupled - channel calculations performed for the @xmath4 system include the @xmath1li ( @xmath1be ) ground state , the first excited @xmath139 state as well as the second excited @xmath176 state .
it is essential to include the @xmath176 state in order to reproduce the low - lying @xmath183 resonance in @xmath7li and @xmath7b . using these three states , we are able to reach model spaces up to @xmath182 , which is sufficient concerning the ho basis expansion convergence as can be judged from fig .  [
fig : li7_itncsm ] .
the coupled channel calculation described above gives two bound states for the @xmath28-@xmath1li system , a @xmath31 corresponding to the experimentally observed @xmath7li ground state , bound by 2.03 mev  @xcite , and a @xmath184 corresponding to the @xmath7li first excited state at @xmath185 mev , bound by 1.05 mev  @xcite .
the calculated states are bound by 1.16 mev and 0.17 mev , respectively , i.e. less than in experiment .
this is in part due to the fact that higher excited states of @xmath1li were omitted . in fig .
[ fig : li7_it_317_pwaves ] , we present our results for the diagonal @xmath53-wave phase shifts of the @xmath28+@xmath1li elastic scattering as well as the elastic @xmath1li(@xmath28,@xmath28)@xmath1li and inelastic @xmath1li(@xmath28,@xmath180)@xmath1li(1/2@xmath181 ) cross sections . at low energies , we can identify four resonances two of which can be associated with the experimentally known @xmath7li states : @xmath183 at @xmath186 mev and @xmath184 at @xmath187 mev  @xcite .
the other two resonances , @xmath188 and @xmath31 are not present in the @xmath7li evaluation of ref .
@xcite . they do appear in many theoretical calculations including the gfmc  @xcite , ncsm  @xcite and recoil - corrected continuum shell model ( rccsm )  @xcite .
the @xmath188 resonance also appears in the gcm calculations of ref .
contributions of different resonances to the cross sections can be deduced from fig .
[ fig : li7_it_317_pwaves ] .
the elastic cross section is dominated by the @xmath183 resonance with some contributions from the @xmath31 resonace at higher energy .
the inelastic cross section shows a peak just above the threshold due to the @xmath188 resonance and also a contribution from the @xmath184 resonance .
the appearance of a @xmath188 peak just above threshold of the @xmath1li(@xmath28,@xmath180)@xmath1li(1/2@xmath181 ) reaction was also discussed in ref .
@xcite ( see fig . 10 in that paper ) .
the data of ref .
@xcite seem to rule out a @xmath188 state so close to the threshold .
it is known , however , that the position of the @xmath188 state is sensitive to the strength of the spin - orbit interaction  @xcite .
the three - nucleon interaction , that would increase the strenght of the spin - orbit force , was not included in our present calculations . consequently , our predicted @xmath188 state energy is likely underestimated .
we note that no fit to the experimental treshold was done in the present ncsm / rgm calculations .
still , as seen in the bottom panel of fig .
[ fig : li7_it_317_pwaves ] , the calculated inelastic cross section is very close to the experimental data just above the threshold .
-wave diagonal phase shifts of the @xmath29+@xmath1be elastic scattering ( top and middle panel ) and inelastic @xmath1be(@xmath29,@xmath189)@xmath1be(1/2@xmath181 ) cross section ( bottom panel ) .
the ncsm / rgm calculation that included the @xmath1be ground state and the @xmath139 and @xmath176 excited states were done using the srg - n@xmath26lo @xmath11 potential with a cutoff of 2.02 @xmath145 .
wave functions from it - ncsm calculations in the @xmath182 basis and the ho frequency of @xmath150 mev were employed . in the middle panel , the full - space ncsm ( solid lines ) and the it - ncsm ( dashed lines ) results in the @xmath154 basis are compared . ]
-wave diagonal phase shifts of the @xmath29+@xmath1be elastic scattering ( top and middle panel ) and inelastic @xmath1be(@xmath29,@xmath189)@xmath1be(1/2@xmath181 ) cross section ( bottom panel ) .
the ncsm / rgm calculation that included the @xmath1be ground state and the @xmath139 and @xmath176 excited states were done using the srg - n@xmath26lo @xmath11 potential with a cutoff of 2.02 @xmath145 .
wave functions from it - ncsm calculations in the @xmath182 basis and the ho frequency of @xmath150 mev were employed . in the middle panel , the full - space ncsm ( solid lines ) and the it - ncsm ( dashed lines ) results in the @xmath154 basis are compared . ]
-wave diagonal phase shifts of the @xmath29+@xmath1be elastic scattering ( top and middle panel ) and inelastic @xmath1be(@xmath29,@xmath189)@xmath1be(1/2@xmath181 ) cross section ( bottom panel ) .
the ncsm / rgm calculation that included the @xmath1be ground state and the @xmath139 and @xmath176 excited states were done using the srg - n@xmath26lo @xmath11 potential with a cutoff of 2.02 @xmath145 .
wave functions from it - ncsm calculations in the @xmath182 basis and the ho frequency of @xmath150 mev were employed . in the middle panel , the full - space ncsm ( solid lines ) and the it - ncsm ( dashed lines ) results in the @xmath154 basis are compared . ]    in the mirror system , @xmath29-@xmath1be , we do not find a bound state in the same type of coupled - channel ncsm / rgm calculation as described above for @xmath28-@xmath1li .
as seen in the top and the middle parts of fig .
[ fig : be7_317_pwaves ] , the lowest @xmath31 resonance corresponding to the @xmath7b ground state lies at about 200 kev above the threshold . in experiment , @xmath7b is bound by 137 kev  @xcite .
our calculated lowest @xmath184 resonance appears at about 1 mev .
it corresponds to the experimental @xmath7b @xmath184 state at @xmath190 mev ( 0.63 mev above the @xmath29-@xmath1be threshold ) .
this resonance dominates the inelastic cross section as seen in the bottom part of fig .
[ fig : be7_317_pwaves ] .
the higher lying resonances follow similar patterns as those found in @xmath28-@xmath1li ( fig .
[ fig : li7_it_317_pwaves ] ) .
again , we find @xmath188 and @xmath31 resonances not included in the recent @xmath7b evaluation  @xcite .
we note that experimental efforts are now under way to find these resonances  @xcite .
we further note that our calculated @xmath191 states in @xmath7li and @xmath7b appear at a significantly higher energies than the corresponding @xmath191 states obtained within the microscopic cluster model in ref .
@xcite .    the middle panel of fig .
[ fig : be7_317_pwaves ] demonstrates once again the good accuracy of the importance truncated calculations for a high @xmath25 , @xmath154 , model space .
the it calculation reduced the @xmath1be basis from 43.6 million to 11.9 million in the present case .
be(@xmath29,@xmath29)@xmath1be ( top panel ) and inelastic @xmath1be(@xmath29,@xmath189)@xmath1be(1/2@xmath181 ) ( bottom panel ) differential cross section at @xmath192 calculated within the ncsm / rgm with srg - n@xmath26lo @xmath11 potential with @xmath178 @xmath145 . ]    be(@xmath29,@xmath29)@xmath1be ( top panel ) and inelastic @xmath1be(@xmath29,@xmath189)@xmath1be(1/2@xmath181 ) ( bottom panel ) differential cross section at @xmath192 calculated within the ncsm / rgm with srg - n@xmath26lo @xmath11 potential with @xmath178 @xmath145 . ]
the elastic @xmath29-@xmath1be scattering was measured at @xmath193 and analyzed by the r - matrix approach  @xcite .
cross section calculations within the rccsm at that angle were then published in ref .
@xcite and also in ref .
further , elastic and inelastic cross sections at this angle were analyzed within the time - dependent approach to the continuum shell model ( tdcsm )  @xcite .
our elastic and inelastic differential cross section results at @xmath193 are presented in fig .
[ fig : p_be7_148 ] . in the elastic cross section ,
the first @xmath184 state is visible and beyond the minimum of the cross section , we can see the dominant peak due to the @xmath183 state . at higher energies
, the @xmath31 state contributes as well .
the inelastic cross section at @xmath193 has a similar shape as the reaction cross section shown in fig .
[ fig : be7_317_pwaves ] .
the first @xmath184 state peak dominates at low energy with contributions from the @xmath188 and the second @xmath184 at higher energies .
our findings are in line with the rccsm results .
however , we remind the reader that there is no fitting in our calculations , all results being predictions based on the srg - n@xmath26lo @xmath11 potential . because of this , the positions of our calculated resonances , e.g. , @xmath184 , @xmath183 do not exactly reproduce experiment .
we do not include the experimental data in the figure as they would be shifted compared to the calculated peaks .
there are at least two reasons why our predictions do not match the experimental resonances accurately .
first , our nuclear hamiltonian is incomplete , e.g. no three - nucleon interaction is included .
second , we omitted higher resonances of @xmath1li and @xmath1be due to numerical reasons .
most likely , the omitted resonances would produce some shifts in the calculated peaks .
both these points can and will be addressed in the future . still , our current results contain the bulk of the physics behind the investigated scattering processes .
-wave phase shifts of the @xmath28+@xmath1li ( solid lines ) and the @xmath29+@xmath1be ( dashed lines ) elastic scattering .
the calculations as described in figs .
[ fig : li7_it_317_pwaves ] and  [ fig : be7_317_pwaves ] . ]    .
the @xmath28-@xmath1li and the @xmath29-@xmath1be @xmath141-wave scattering lengths .
theoretical values correspond to calculations as described in figs .
[ fig : li7_it_317_pwaves ] and  [ fig : be7_317_pwaves ] .
experimental values are from refs .
@xcite . [ cols="^,^,^,^,^ " , ]     in fig .
[ fig : be7_li7_317_swaves ] , we present our calculated @xmath28-@xmath1li and the @xmath29-@xmath1be @xmath141-wave phase shifts
. we do not find any evidence for a @xmath194 resonance advocated in ref .
@xcite and discussed in ref .
the corresponding scattering lengths together with the experimental values are given in table  [ tab : s - wave_scatl ] .
with the exception of the @xmath29-@xmath1be @xmath195 , which has a large experimental uncertainty , our calculated scattering lengths do agree with experimental data as to their signs , there are however differences in the absolute values .
again , as dicussed above , the results presented here serve only as a first step towards the _ ab initio _ investigation of the @xmath28-@xmath1li and the @xmath29-@xmath1be reactions .
prospects for a realistic calculation of the @xmath1be(@xmath29,@xmath175)@xmath7b capture are excellent . here
we found the @xmath7b unbound by only 200 kev .
it is quite possible that @xmath7b will become bound ( with the @xmath11 potential employed here : srg - n@xmath26lo with @xmath178 @xmath145 ) by including more excited states of @xmath1be in the coupled - channel ncsm / rgm calculations .
the effect of higher excited states of @xmath1be can be , in fact , most efficiently included by coupling the presently used ncsm / rgm basis with the @xmath7b ncsm eigenstates as outlined in ref .
even if the @xmath7b would not be bound or , most likely , the theshold energy will not agree with the experiment , we have the possibilty to explore a variation of the srg @xmath11 potential evolution parameter @xmath196 and tune this parameter to fit the experimental threshold .
we note that for any @xmath196 the srg - evolved @xmath11 potential will describe all two - nucleon properties as accurately as the original starting @xmath11 potential , here the chiral n@xmath26lo potential of ref .
it should be noted that by adding the three - nucleon interaction , omitted in the present calculations due to computational reasons , the need for a fine - tuning should be significantly reduced , i.e. the results should become @xmath196 independent .
for nucleon scattering calculations on @xmath2c or heavier targets within the ncsm / rgm , the use of the importance truncation becomes essential . for @xmath2c ,
the full - space ncsm calculations are currently limited to @xmath155 ( although successful runs were already performed for @xmath154 on the biggest supercomputers with the latest version of the code mfd  @xcite ) .
this is insufficient for reaching or approaching convergence of the @xmath2c ncsm calculations as seen from fig .
[ fig : c12_itncsm ] and even more so of the ncsm / rgm scattering calculations . the importance - truncated calculations , on the other hand , are feasible up to @xmath177 , where convergence is reached for both the ground state as well excited states .
our @xmath2c calculations are performed with the srg - n@xmath26lo @xmath11 potential with the evolution parameter @xmath197 @xmath145 , a higher value ( i.e. shorter evolution , less soft ) than that used for the lighter nuclei . the use of a small @xmath196 results in large overbinding of heavier nuclei and a significant underestimation of their radii . as seen in fig .
[ fig : c12_itncsm ] , our converged @xmath2c binding energy is about 84.5(8 ) mev , smaller than the experimental value of 92 mev and , further , the agreement of the full - space and importance - truncated results is perfect all the way up to @xmath155 .
state energy dependence on the model - space size @xmath70 for @xmath2c , obtained within the importance - truncated ncsm , using the srg - n@xmath26lo @xmath11 potential with a cutoff of 2.66 @xmath145 .
the ho frequency @xmath198 mev was employed .
the calculation is variational .
no ncsm effective interaction was used .
the full ncsm results were obtained with the code antoine  @xcite . ]
our low - energy @xmath29-@xmath2c phase shift results are shown in fig .
[ fig : pc12 ] .
the comparison of the @xmath157 and @xmath199 results demonstrates good convergence with respect to the ho basis expansion . the @xmath2c ground state and the first @xmath31 state were included in the coupled - channels ncsm / rgm equations .
we note that we also performed a phase shift comparison of the full - space and the importance - truncated calculations up to @xmath200 and found a similarly perfect agreement as presented in fig .
[ fig : nhe4_phase_full_it ] for @xmath28-@xmath0he .
-@xmath2c eigenphase shifts calculated within the ncsm / rgm using the srg - n@xmath26lo @xmath11 potential with a cutoff of 2.66 @xmath145 and the ho frequency @xmath198 mev . full lines ( dotted lines ) corespond to results obtained in the @xmath157 ( @xmath199 ) model space .
the ground state and the first excited @xmath31 state of @xmath2c was included .
the @xmath2c wave functions were obtained within the it ncsm . ]    in the present @xmath29-@xmath2c calculations , we found a single bound state , @xmath139 at -2.98 mev , corresponding to the @xmath30n ground state , bound experimentally by 1.94 mev  @xcite .
the lowest resonance in our calculation is @xmath201 , barely visible at 0.25 mev above threshold . in experiment , this resonance is at 1.56 mev .
our calculated @xmath140 resonance appears at about 1.5 mev above treshold ( in experiment at 0.42 mev above threshold ) and the @xmath202 resonance at about 4.9 mev ( in experiment at 2.61 mev ) .      in the mirror system , @xmath28-@xmath2c , our ncsm /
rgm calculations produce three bound states : @xmath139 at -5.34 mev corresponding to the @xmath30c ground state experimentally bound by 4.95 mev with respect to the @xmath28-@xmath2c threshold , @xmath201 , bound by 2.23 mev ( experimentally bound by 1.26 mev ) , and @xmath140 bound by just 0.03 mev ( experimentally bound by 1.86 mev ) . in experiment , there is also a @xmath202 state bound by 1.09 mev . our present ncsm / rgm calculations including the lowest @xmath188 and and the lowest @xmath31 @xmath2c states do not produce any bound @xmath202 state .
-@xmath2c phase shifts calculated within the ncsm / rgm using the srg - n@xmath26lo @xmath11 potential with a cutoff of 2.66 @xmath145 .
the ho frequency @xmath198 mev and the model - spaces size of @xmath157 were used .
the ground state and the first excited @xmath31 state of @xmath2c was included .
the @xmath2c wave functions were obtained within the it ncsm . ]
our low - energy @xmath28-@xmath2c diagonal phase shifts are shown in fig .
[ fig : nc12 ] .
the @xmath202 resonance is found at 2.8 mev ( experimenally at 1.92 mev with respect to the @xmath28-@xmath2c threshold ) .
the steep drop of the @xmath140 phase shift is due to the presence of the very weakly bound @xmath140 state .
we note that similarly as in the case of @xmath21be , discussed in ref .
@xcite , we observe a significant decrease of the @xmath140 state energy in the @xmath28-@xmath2c ncsm / rgm calculation when compared to the standard ncsm calculation for @xmath30c .
we were able to make these comparisons in model spaces up to @xmath200 where we found this drop to be about 3 mev .
-@xmath2c elastic scattering below and above the calculated @xmath202 resonance .
energies are in the center of mass .
the calculation as described in fig .
[ fig : nc12 ] . ]
analysing powers were measured for proton and neutron scattering on @xmath2c  @xcite and scattering experiments on polarized proton target are under way  @xcite . in fig .
[ fig : nc12_ay ] , we present our calculated analyzing power below and above the energy of the @xmath202 resonance . we note that our calculated @xmath202 resonance appears at 2.8 mev in the center of mass ( experimentally at 1.92 mev ) .
below the resonance , the analyzing power is positive at @xmath203 and negative at @xmath204 . at energies above the resonance ,
the analyzing power reverses its sign .
similar observations were made in calculations performed within the multichannel algebraic scattering ( mcas ) theory  @xcite .
see in particular fig .
5 of ref .
@xcite .
our calculated @xmath30n and @xmath30c bound - state levels and resonances are more spread than the experimental ones .
this is a consequence of an underestimation of the @xmath2c radius found to be 2.05 fm with the srg - n@xmath26lo @xmath11 potential . to remedy this
, one would have to calculate three - nucleon interaction terms induced due to the srg evolution .
this can be done as described in ref .
however , we still need to further develop the ncsm / rgm formalism in order to handle three - nucleon interactions in the scattering calculations .
the calculation of nucleon scattering on @xmath3o is the most challenging among the systems we investigate in this paper .
the @xmath10 clustering plays an important role in the structure of @xmath3o , in particular for the first excited @xmath188 state that is known to be almost impossible to reproduce in ncsm or coupled - cluster calculations .
our present calculations do not include the @xmath10 clustering yet .
as in the case of @xmath2c , we rely on the importance - truncated ncsm calculations for obtaining the @xmath3o wave functions as the full-@xmath70 ncsm calculations are possible only up to @xmath155 . in fig .  [
fig : o16_itncsm ] , we show the ground - state convergence within the it - ncsm and a comparison to the full - space results .
again , up to the largest accessible model space , the agreement between the importance - truncated and the full - space calculations is perfect .       for @xmath3o , obtained within the importance - truncated ncsm , using the srg - n@xmath26lo @xmath11 potential with a cutoff of 2.66 @xmath145 .
the ho frequency @xmath198 mev was employed .
the calculation is variational .
no ncsm effective interaction was used .
the full ncsm results were obtained with the code antoine  @xcite . ]
o phase shifts calculated within the ncsm / rgm using the srg - n@xmath26lo @xmath11 potential with a cutoff of 2.66 @xmath145 and the ho frequency @xmath198 mev in the @xmath177 model space . the ground state and of @xmath3o was included .
the @xmath3o wave functions were obtained within the it ncsm . ]
o phase shifts calculated within the ncsm / rgm using the srg - n@xmath26lo @xmath11 potential with a cutoff of 2.66 @xmath145 . the ho frequency of @xmath198 mev was used .
the @xmath205 channel is shown in the top ( bottom ) panel .
model space sizes up to @xmath177 were considered .
the ground state and of @xmath3o was included .
the @xmath3o wave functions were obtained within the it ncsm . ]
o phase shifts calculated within the ncsm / rgm using the srg - n@xmath26lo @xmath11 potential with a cutoff of 2.66 @xmath145 . the ho frequency of @xmath198 mev was used .
the @xmath205 channel is shown in the top ( bottom ) panel .
model space sizes up to @xmath177 were considered .
the ground state and of @xmath3o was included .
the @xmath3o wave functions were obtained within the it ncsm . ]
it is straightforward to converge nucleon-@xmath3o scattering calculations within the ncsm / rgm using the ho expansion up to @xmath177 .
our calculated @xmath28-@xmath3o phase shifts are shown in fig .
[ fig : no16gs_18 ] and the ho - basis expansion convegence is checked for the @xmath141- and the @xmath166-wave in fig .
[ fig : no16gs_conv ] .
these calculations included the @xmath3o ground state only .
we find two bound states , @xmath140 at -0.88 mev and @xmath202 at -0.41 mev with respect to the @xmath28-@xmath3o threshold . in experiment ,
the @xmath206o ground state is @xmath202 , bound by 4.14 mev , and the @xmath140 state is the first excited state bound by 3.27 mev .
there are also two additional bound states : @xmath139 and @xmath201 .
those are unbound in our calculations .
-@xmath3o phase shifts calculated within the ncsm / rgm using the srg - n@xmath26lo @xmath11 potential with a cutoff of 2.66 @xmath145 and the ho frequency @xmath198 mev in the @xmath207 model space . the ground state and and the lowest @xmath208 , @xmath209 and @xmath194 excited states of @xmath3o
were included .
the @xmath3o wave functions were obtained within the it ncsm . ]    clearly , it is insufficient to consider only the ground state of @xmath3o in the coupled - channel ncsm / rgm scattering calculations .
we , therefore , include in addition the three lowest @xmath3o negative parity states : @xmath208 , @xmath209 , and @xmath194 .
due to computational limitations , in this case we used ho basis expansion up to @xmath70=13 . comparing fig .
[ fig : no16gs312_15 ] to fig .
[ fig : no16gs_18 ] , the @xmath210 negative - parity excited states of @xmath3o generate negative - parity resonances in @xmath206o .
these resonances do appear , however , at much higher energy than in experiment .
the reason for this is the fact that our calculated @xmath3o @xmath210 states have too large excitation energy . in particular , our calculated @xmath208 excited state
has an excitation energy of 15.99 mev while experimentally it lies at just 6.13 mev .
one reason for the discrepancy is the softness of the srg - n@xmath26lo @xmath11 potential we use that results in an overall overbinding of the @xmath3o ground state and in an underestimation of its radius .
another aspect is the challenging problem of the it - ncsm extrapolations of the independent positive and negative - parity state calculations .
the uncertainties of the relative excitation energies are higher than in same - parity calculations . on the positive side our calculation with the negative - parity states , even though with overestimated excitation energies , results in the proper ordering of the @xmath206o bound states .
the ground state is @xmath202 at -1.32 mev and the @xmath140 state gains binding as well , appearing at -1.03 mev .
we also investigated the @xmath29-@xmath3o scattering and @xmath206f states .
when the ncsm / rgm calculations are restricted to the channels involving only the @xmath3o ground state , we find a @xmath140 resonance at 1.0 mev and a @xmath202 resonance at 2.2 mev .
these resonances correspond to the @xmath206f @xmath140 first excited state , bound by 0.105 mev , and the @xmath206f @xmath202 ground state bound by 0.6 mev with respect to the @xmath29+@xmath3o threshold . by coupling channels involving the @xmath210 @xmath3o @xmath208 , @xmath209 and @xmath194 excited states , the calculated @xmath140 and @xmath202 states are still unbound resonances but their energy moves significantly closer to the threshold : the @xmath140 appears at + 0.7 mev and the @xmath202 at + 1.2 mev .
the @xmath206f low - lying states were recently investigated within the coupled - cluster approach with the gamow - hartree - fock basis  @xcite . in those calculations with the n@xmath26lo @xmath11 potential ,
the @xmath140 state is weakly bound while the @xmath202 state remains unbound by about 0.1 mev .
using the srg evolved interaction , the @xmath202 state became bound with the decrease of the cutoff @xmath196 .
we note that our calculated @xmath3o ground state energy , -139.0(8 ) mev ( fig .
[ fig : o16_itncsm ] ) obtained with the srg - n@xmath26lo @xmath11 potential with @xmath197 @xmath145 , compares well with the ccsd coupled - cluster @xmath3o calculations : -137.6 mev with the srg - n@xmath26lo @xmath11 potential with @xmath211 @xmath145  @xcite .
the differences in the positions of the @xmath140 and the @xmath202 are due to deficiencies in our description of the negative parity @xmath210 states , which could be related to the two - body hamiltonian used here as well as the uncertainties of the threshold extrapolations for the excitation energies .
the inclusion of additional @xmath3o excited states would increase the absolute energy of our calculated @xmath206f states . the most efficient way to do
this is by coupling the presently used ncsm / rgm basis with the @xmath206f ncsm eigenstates in as outlined in ref .
by combining the importance truncation scheme for the cluster eigenstate basis with the _ ab initio _ ncsm / rgm approach , we were able to perform many - body calculations for nucleon scattering on nuclei with mass number as high as @xmath212 . with the soft srg - evolved chiral @xmath11 potentials , convergence of the calculations with respect to the ho basis expansion of the target eigenstates and the localized parts of the ncsm / rgm integration kernels
can be reached using @xmath213 .
we first benchmarked the it - ncsm results with the full - space ncsm results for the @xmath23 system .
our neutron-@xmath0he and proton-@xmath0he calculations compare well with an r - matrix analysis of the data in particular at energies above 8 mev , and describe well measured cross sections and analysing powers for those energies .
our calculations of @xmath28-@xmath1li and @xmath29-@xmath1be scattering predict low - lying @xmath188 and @xmath31 resonances in @xmath7li and @xmath7b that have not been experimentally clearly identified yet .
we found that the prospects of a realistic _ ab initio _ calculation of the @xmath1be(@xmath29,@xmath175)@xmath7b capture within our approach are very good . in the present calculations we found the @xmath7b unbound by only 200 kev .
it is quite possible that @xmath7b will become bound ( with the @xmath11 potential employed here : srg - n@xmath26lo with @xmath178 @xmath145 ) by including more excited states of @xmath1be in the coupled - channel ncsm / rgm calculations .
even if the @xmath7b will still not be bound or , most likely , the threshold energy will not agree with the experiment , we have the possibility to explore a variation of the srg @xmath11 potential evolution parameter @xmath196 and tune this parameter to fit the experimental threshold .
the use of the importance - trunacted basis becomes essential in calculations with @xmath2c or @xmath3o targets as the full - space ncsm calculations are limited to @xmath155 .
our @xmath28-@xmath2c and @xmath29-@xmath2c investigations included @xmath2c ground and the first excited @xmath31 states .
we found a single bound state , @xmath140 in @xmath30n as in experiment . in @xmath30c , we found three bound states with the @xmath202 state still unbound contrary to experiment .
our calculated spectrum of @xmath8 states is more spread than in experiment due to the underestimation of the @xmath2c radius , a consequence of the softness of the srg - evolved @xmath11 interaction .
the description of nucleon scattering on @xmath3o within our formalism was the most challenging .
the @xmath10 clustering that plays an important role in the structure of @xmath3o is not yet included in our present calculations .
further , the @xmath210 @xmath3o excited states are more difficult to treat in the it - ncsm approach , as the extrapolations of excitation energies are done from the independent ground state and the negative - parity state calculations . we found a strong impact of the @xmath210 @xmath3o states on the positions of the lowest @xmath9 states .
for example , correct ordering of the @xmath202 and the @xmath140 states in @xmath206o was obtained only when the @xmath210 states were included .
overall , we find that the inclusion of additional excited states of the target nuclei would be beneficial in all studied systems and more significant with the increase of @xmath14 .
coupled - channel ncsm / rgm calculations with many excited states of the traget are computationally challenging .
the most efficient way of including the effects of such states is by coupling the presently used ncsm / rgm basis , consisting of just a few lowest excited states , with the ncsm eigenstates of the composite system as outlined in ref .  @xcite
. work on this coupling is under way .
the use of the srg - evolved @xmath11 interaction facilitates convergence of the ncsm / rgm calculations with respect to the ho basis expansion . on the other hand , due to the softness of these interactions , radii of heavier nuclei
become underestimated . to remedy this
, one would have to calculate three - nucleon interaction terms induced due to the srg evolution .
this can be done as described in ref .
it is essential to further develop the ncsm / rgm formalism in order to handle three - nucleon interactions , both genuine and the srg - evolution induced , in the scattering calculations .
numerical calculations have been performed at the llnl lc facilities and at the nic , jlich . prepared in part by llnl under contract de - ac52 - 07na27344 .
support from the u.  s.  doe / sc / np ( work proposal no .
scw0498 ) , llnl ldrd grant pls-09-erd-020 , and from the u.  s.  department of energy grant de - fc02 - 07er41457 is acknowledged .
this work is supported in part by the deutsche forschungsgemeinschaft through contract sfb 634 and by the helmholtz international center for fair within the framework of the loewe program launched by the state of hesse .
+          r. b. wiringa , s. c. pieper , j. carlson , v. r. pandharipande , phys .
c * 62 * , 014001 ( 2000 ) ; s. c. pieper and r. b. wiringa , ann .
nucl . part .
sci . * 51 * , 53 ( 2001 ) ; s. c. pieper , k. varga and r. b. wiringa , phys .
c * 66 * , 044310 ( 2002 ) .