the @xmath0he three - atomic system is of considerable interest in various fields of physical chemistry and molecular physics .
studies of the helium dimer and trimer represent an important step towards understanding the properties of helium liquid drops , superfluidity in @xmath0he films , and so on ( see , for instance , refs .
besides , the helium trimer is probably a unique system where a direct manifestation of the efimov effect  @xcite can be observed since the binding energy @xmath2 of the @xmath0he dimer is extremely small ( @xmath3mk  @xcite ) even in the molecular scale . for this reason ,
the helium trimer is certainly of interest for nuclear physicists , too .
moreover a theoretical study of the @xmath0he trimer is based just on the same methods of the theory of few - body systems that are used in solving three  body nuclear problems .    from the standpoint of the general theory of few
body systems , the @xmath0he trimer belongs to three  body systems that are most difficult for a specific investigation , first , owing to its efimov nature , and second , because it is necessary to take into account the practically hard core in the interatomic he  he - interaction  @xcite . at the same time
the problem of three helium atoms can be considered as an example of an ideal three  body quantum problem since the @xmath0he atoms are identical neutral bosons with zero spin and the analysis of this problem is complicated neither by separation of spin  isospin variables nor by taking into account the coulomb interaction .
there is a great number of experimental and theoretical studies of @xmath0he clusters .
however , most of the theoretical investigations consist merely in computing the ground - state energies of clusters of that sort , mainly on the basis of variational methods @xcite . besides , the methods based on hyperspherical expansions of the schrdinger and faddeev equations @xcite in the coordinate representation were used . also , the faddeev integral equations in the momentum representation were employed in refs .
@xcite while the results of ref .
@xcite are based on a direct solving the two - dimensional faddeev differential equations in configuration space . from the experimental studies we would like to mention those of refs .
@xcite where clusters consisting of a small number of noble gase atoms were investigated .
though much effort was undertaken for studying molecular clusters various problems related to the @xmath0he trimer remained beyond the scope of thorough consideration .
in particular , the elastic scattering phases of a helium atom on a helium dimer and breakup amplitudes ( at ultralow energies ) have been calculated only recently  @xcite .
these computations were preceded only by the computation of characteristics of the he
he@xmath4 scattering at zero energy  @xcite and estimation of the recombination rate @xmath5  @xcite .    as a matter of fact
, we have already pointed out basic reasons for computations of excited states and scattering being especially difficult in the @xmath0he@xmath6 system .
first , this is a low energy of the dimer @xmath2 which necessitates to consider very large domains in the configuration space with a characteristic size of hundreds of  .
second , a very strong repulsive component in the he  he interaction produces large errors in the standard approximation of the three
atomic hamiltonian at short distances between atoms .
the capacities of modern computers do not yet allow one to reach dimensions of grids that would remove both the above - mentioned reasons and would provide stable results with the use of the conventional methods .    the present paper is a sequel of studies of the @xmath0he@xmath6 system undertaken in the papers  @xcite within an approach that is capable , as we think , to resolve both the above - mentioned numerical problems .
in these papers the repulsive component of the he  he interaction at short distances between atons is approximated by a hard core .
this allows one to investigate the @xmath0he@xmath6 system within a mathematically rigorous method of solving a three - body problem in the boundary - condition model developed in  @xcite .
an important advantage of such an approach that essentially diminishes computational errors is the necessity to approximate , inside the core domains , only the laplacian operator instead of the sum of this operator and a huge repulsive components of the he  he - potentials ( see @xcite ) . in  @xcite ,
such an approach has been successfully applied for calculating not only scattering but also binding energies of the ground and excited states of the helium trimer .
investigation made in  @xcite has shown that the method proposed in  @xcite is well suited for performing three
body molecular computations in the case where repulsive components of interatomic interactions are of a hard core nature .
there is a series of works @xcite showing that the excited state of the @xmath0he trimer is initiated indeed by the efimov effect  @xcite . in these works the various versions of the aziz @xmath0he@xmath0he potential were employed ( hfdhe2  @xcite , hfd - b  @xcite , and lm2m2  @xcite ) . however , the basic result of refs .
@xcite on the excited state of the helium trimer is the same : this state disappears when the interatomic potential is multiplied by the  amplification factor " @xmath7 of order 1.2 .
more precisely , if this potential is multiplied by the increasing factor @xmath8 then the following effect is observed .
first , the difference @xmath9 between the dimer energy @xmath10 and the energy of the trimer excited state @xmath11 increases
. then the behavior of this difference radically changes and with further increase of @xmath7 it monotonously decreases . at @xmath12 the level @xmath13 disappears .
it is just such a nonstandard behavior of the energy @xmath11 as the coupling between helium atoms becomes more and more strengthening , points to the efimov nature of the trimer excited state . and vice versa , when @xmath7 slightly decreases ( no more than 2% ) , the second excited state @xmath14 appears in the trimer  @xcite .
this paper is aimed at elucidating the fate of the trimer excited state upon its disappearance in the physical sheet when @xmath8 and at studying the mechanism of arising of new excited states when @xmath15 .
as the interatomic he  he potential , we use the potential hfd - b  @xcite .
we have established that for such he  he - interactions the trimer excited level @xmath13 merges with the threshold @xmath2 at @xmath16 and with further decreasing @xmath7 it transforms into a virtual level of the first order ( a simple real pole of the analytic continuation of the scattering matrix ) lying in the unphysical energy sheet adjoining the physical sheet along the spectral interval between @xmath2 and the three  body threshold .
we trace the position of this level for @xmath7 increasing up to 1.5 .
besides , we have found that the excited ( efimov ) levels for @xmath15 also originate from virtual levels of the first order that are formed in pairs .
before a pair of virtual levels appears , there occurs a fusion of a pair of conjugate resonances of the first order ( simple complex poles of the analytic continuation of the scattering matrix in the unphysical sheet ) resulting in the virtual level of the second order .    as it will be clear from the further exposition ( see sect .  [ results ] ) ,
the above - mentioned resonances are not , generally speaking , genuine resonances of the @xmath0he@xmath6 trimer since they are situated outside of the energy domain for which we can rigorously prove the applicability of the method we are using for computing the resonances .
we will call the resonances found outside the range of guaranteed applicability of the method the _ ( quasi)resonances_.    the paper is organized as follows .    in sect .
[ difffs ] , we describe the method of search for resonances in a three  body system on the basis of the faddeev differential equations .
the idea of the method consists in calculating the analytic continuation of the component @xmath17 ( see formula  ( [ s0 ] ) ) of the scattering matrix corresponding to the ( @xmath18 ) process , in the physical sheet with the use of these equations .
a particular attention in this section is paid to the description of the parabolic domain on the physical sheet where one can analytically continue the function @xmath19 by numerical solving the coordinate space faddeev partial equations . for the potentials we use , the three  body resonances ( including virtual levels ) lying in the unphysical sheet of energy @xmath20 plane adjoining the physical sheet along the interval @xmath21 are the roots of the function @xmath19 in the physical sheet .
we have earlier employed this method for computing resonances as roots of @xmath17 in the three  nucleon problem  @xcite .    in sect .
[ results ] , we first briefly describe the numerical method we use to solve the @xmath1 scattering problem for the @xmath0he@xmath6 system with going out into the domain of complex energies .
then we describe the results of our calculations .
some notation used throughout the paper is as follows : by @xmath22 we denote the complex plane ; @xmath23 stands for the main branch of the function @xmath24 , @xmath25 for any @xmath26 ; the symbol @xmath27 is used for the quadrant @xmath28 , @xmath29 ; by @xmath30 we understand the hilbert space of complex  valued functions which are integrable on @xmath27 with the absolute value squared ; the symbol @xmath31 stands for the complex number conjugated to @xmath20 .
in this paper , we will consider the @xmath0he@xmath6 system in the state with the total angular momentum @xmath32 .
first we consider the case where the interatomic interactions are described by conventional smooth potentials that include no hard - core component . in this case , the angular partial analysis reduces the initial faddeev equation for three identical bosons to a system of coupled two - dimensional integro - differential equations  @xcite @xmath33f_l(x , y)=-v(x)\psi_l(x , y)\,.\ ] ] here , @xmath34 stand for the standard jacobi variables , @xmath28 and @xmath29 , and @xmath35 for the partial component of the kinetic energy operator .
functions from the domain of @xmath36 are assumed to obey the boundary conditions @xmath37 which are quite standard when the expansions over bispherical basis are used .
the potential @xmath38 is assumed to be central . in our paper
, the energy @xmath20 can get both real and complex values . at @xmath32
the partial angular momentum @xmath39 corresponds both to the dimer and an additional atom .
the momentum @xmath39 can assume only even values , @xmath40    the partial wave functions @xmath41 are expressed through the faddeev partial components @xmath42 by the relations @xmath43 where @xmath44 and @xmath45 .
the explicit form of the functions @xmath46 can be found , e.g. , in refs .
@xcite ( see also  @xcite ) .
here we only deal with a finite number of equations  ( [ fadpart ] ) , assuming that @xmath47 where @xmath48 is a certain fixed even number , @xmath49 .
the condition @xmath50 is equivalent to the supposition that the potential @xmath38 only acts in the two - body states with @xmath51 .
the spectrum of the schrdinger operator for a system of three identical bosons with such a potential is denoted by @xmath52 .
its is well known ( see , e.g. , ref .
@xcite ) that if the potential @xmath38 is smooth and decreasing as @xmath53 together with its derivatives not slower than @xmath54 , @xmath55 , then the asymptotic conditions as @xmath56 and/or @xmath57 for the partial faddeev components of the @xmath58 scattering wave functions '' .
the asymptotics of these functions in the total three - body configuration space @xmath59 contains , apart from the incident wave , only the so - called outgoing spherical waves ( see , e.g. , @xcite ) . ] for @xmath60 , @xmath61 , read@xmath62\right\ } \\        & & +    \displaystyle\frac{\exp({\rm i}\sqrt{z}\rho)}{\sqrt{\rho } }                  \left[a_l(z,\theta)+o\left(1\right)\right]\ , .
\end{array}\ ] ] we assume that the @xmath0he@xmath4 dimer has an only bound state with an energy @xmath2 , @xmath63 , and wave function @xmath64 .
this function is assumed to be normalized so that for all @xmath65 values of @xmath64 are real .
the notations @xmath66 , @xmath67 , and @xmath68 , @xmath69 , are used for the hyperradius and hyperangle .
the coefficient @xmath70 , @xmath60 , for @xmath71 is the elastic scattering amplitude .
the functions @xmath72 provide us , at @xmath61 , the corresponding partial faddeev breakup amplitudes . note that for @xmath60 the correction terms @xmath73 in coefficients of outgoing waves @xmath74 , @xmath71 , and @xmath75 , @xmath61 , in  ( [ asbcparts ] ) are of the form , respectively , @xmath76 and @xmath77 .
this property ensures uniqueness of the solution of the boundary value problem  ( [ fadpart][asbcparts ] ) for real scattering energies @xmath71  @xcite .
the component of the @xmath78-wave partial scattering matrix for a system of three helium atoms is given for real @xmath60 , @xmath71 , by the expression @xmath79 while the scattering phases read @xmath80 where @xmath81 stands for the momentum conjugated to the jacobi variable  @xmath82 .
our goal is to study the analytic continuation of the scattering matrix @xmath19 into the complex plane ( the physical sheet ) . as it follows from the results of refs .
@xcite , roots of the function @xmath17 in the physical sheet of energy @xmath20 plane correspond to the location of the three - body resonances situated in the unphysical sheet connected with the physical sheet by crossing the spectral interval @xmath21 .
this statement is a particular case of more general statements regarding the three - body resonances obtained in  @xcite for the case of two - body potentials decreasing in the coordinate space not slower than exponentially .
we assume that @xmath38 is just a potential which falls off exponentially and , thus , for all @xmath28 @xmath84 with some positive @xmath85 and @xmath86 .
for the sake of simplicity we even assume sometimes that @xmath38 is finite , i.e. , for @xmath87 , @xmath88 .
looking ahead , we note that , in fact , in our numerical computations of the @xmath0he@xmath6 system at complex energies we make a  cutoff " of the interatomic he  he - potential at a sufficiently large radius @xmath89 .
it is well known that different representations of the same holomorphic function ( for instance , either by a series or by an integral ) allow one to describe this function only in some parts of its riemann surface .
the description  @xcite of the holomorphy domains for different truncations of the total three - body scattering matrix in the physical sheet was based on the use of the faddeev integral equations in the momentum representation . in this paper , we make use of the faddeev equations in the configuration space .
therefore it is necessary to perform an investigation , independent of  @xcite , of domains in the physical sheet where we can analytically continue the faddeev components @xmath90 and the amplitudes @xmath91 and @xmath92 just with the use of the configuration space techniques .
let us list briefly the main results of this investigation obtained by us for the @xmath0he@xmath6 system under the assumption  ( [ vexp ] ) .
to formulate these results we distinguish the following three domains in the complex plane @xmath93 .
the domain @xmath95 where the faddeev components @xmath90 ( and , hence , the wave functions @xmath96 ) can be analytically continued in @xmath20 so that the differences @xmath97 at @xmath98 turn out to be elements of @xmath99 .
the domain @xmath95 is described by the inequality @xmath100 for fixed @xmath34 the functions @xmath101 are continuous in @xmath20 up to the rims of the cut along the continuous spectrum @xmath102 .
the domain @xmath103 where both the elastic scattering amplitude @xmath70 and the faddeev breakup amplitudes @xmath92 can be analytically continued in @xmath20 , @xmath104 , provided that the functions @xmath90 obey the asymptotic formulas  ( [ asbcparts ] ) .
this domain is described by the inequalities @xmath105    3@xmath94 . and
finally , we distinguish the domain @xmath106 , most interesting for us , where the analytic continuation in @xmath20 , @xmath104 , can be only done for the elastic scattering amplitude @xmath70 ( and consequently , for the scattering matrix @xmath19 ) ; the analytic continuabilty of the amplitudes @xmath92 in the whole domain @xmath106 is not required . the set @xmath106 is a geometric locus of points obeying the inequality @xmath107    for the domains @xmath95 , @xmath103 , and @xmath106 , the following chain of inclusions @xmath108 is valid .
note that the type  ( [ pipsi ] ) or ( [ pis ] ) condition , @xmath109 is equivalent to the inequality @xmath110 therefore , for @xmath111 the set @xmath95 is the domain bounded by the parabola @xmath112 for @xmath113 this set coincides with the domain @xmath114 analogously , if @xmath115 then the domain @xmath106 is described by the inequality  ( [ parabmu ] ) ; whereas for @xmath116 by the inequality @xmath117 as to the curves bounding the domains  ( [ piamain ] ) and ( [ piamu ] ) , we only notice that their order with respect to the variables @xmath118 and @xmath119 is higher than the second order .
it is easy to check that each of these curves is connected , symmetric with respect to the @xmath118 axis and crosses the latter only once .
for the first curve this intersection occurs at , the slope angle of the tangent at the point of intersection being independent of @xmath2 , . as @xmath120 , the boundaries  ( [ piamain ] ) and  ( [ piamu ] ) are asymptotically approximated by the type  ( [ parabola ] ) parabolas with coefficients @xmath121 and @xmath122 which can be computed explicitly .    to prove the assertion1@xmath94 concerning the domain
@xmath95 we note that the functions @xmath101 given by the formulas ( [ flphil ] ) satisfy the equations @xmath123\phi_l(x , y;z )   + v(x)\sum_{l'}\int_{-1}^{+1 }   d\eta\,h_{l l'}(x , y,\eta)\,\phi_{l'}(x',y',z)=               \chi_l(x , y;z)\,\ ] ] where @xmath124 obviously , for @xmath125 the functions @xmath126 fall off exponentially as @xmath127 .
moreover , for all the directions @xmath128 the uniform estimate @xmath129 is valid with @xmath130 consequently , if the condition  ( [ pipsi ] ) holds then the inhomogeneous terms @xmath126 considered as functions of the variables @xmath131 and @xmath82 at fixed @xmath20 , are elements of @xmath132 . at the same time
, the vectors @xmath133 turn out to be holomorphic functions of @xmath125 with respect to the @xmath30 norm .    in the problem under consideration , the spectrum of the faddeev matrix operator defined by the l.h.s . of eqs .
( [ faddxi ] ) and by the boundary conditions ( [ bcstandard ] ) in the hilbert space constituted of the vectors , @xmath134 , coincides with the spectrum @xmath52 of the corresponding three - boson schrdinger operator with two - body potentials @xmath38 only acting in the states with @xmath51 .
this means that for any energy @xmath125 lying outside of the spectrum @xmath52 , the inhomogeneous system  ( [ faddxi ] ) is uniquely solvable in the class of the functions @xmath135 , @xmath51 . since outside of the set @xmath52 the resolvent of the faddeev operator is a holomorphic operator - valued function of the variable @xmath20 , each of the components @xmath136 of the solution of eqs .
( [ faddxi ] ) also is a holomorphic function of @xmath98 .
the bound - state energies of the three - boson system under consideration turn out to be poles of the first order for @xmath136 .
thus , the faddeev partial components @xmath90 admit the analytic continuation in @xmath20 in the form  ( [ flphil ] ) into the domain @xmath137 .    the proof of the assertions of 2@xmath94 and 3@xmath94 regarding the domains @xmath106 and @xmath103 is rather cumbersome .
this is why we here only outline its main steps .
note that the proof is based on the integral equations method and it is quite standard ( see , for instance , ref .
@xcite , chapter  v ) .
first , the equations  ( [ faddxi ] ) are rewritten in the form of the faddeev partial integral equations . to do this
, it suffices to reverse the operators @xmath138 in ( [ faddxi ] ) .
since the variables @xmath131 and @xmath82 in @xmath36 are separated , the kernels @xmath139 , @xmath140 , @xmath141 , of the respective resolvents @xmath142 are explicitly expressed in terms of the two - body problem . analytic properties in the variable @xmath20 and coordinate asymptotics of the kernels @xmath139 are well known ( see ref .
@xcite , chapters  iv and  v ) .
iterations first  improve " and then stabilize the asymptotic properties of the iterated kernels and inhomogeneous terms of the faddeev equations .
( in the case under consideration , this stabilization requires only three iterations . )
further , it turns out that , for @xmath143 , the iterated kernels are represented by sums of exponentially decreasing terms admitting , in certain domains of the configuration space , an explicit asymptotic factorization with respect to @xmath144 and @xmath145 .
since we are working in the domain where @xmath146 , the corresponding asymptotic factors of these terms , along with the asymptotics of the iterations of the inhomogeneous term , determine the coordinate asymptotics of the functions @xmath101 .
therefore , finally we are able to determine the geometric locus of the points @xmath20 in the complex plane for which there exists a ( non - empty ) set in the configuration space such that the leading term of the coordinate asymptotics of the function @xmath101 in this set represents a term of the form @xmath147 , and thus , for these @xmath20 the scattering matrix @xmath19 is well defined .
this geometrical locus is just the domain @xmath106 . in this domain , as @xmath148 and/or @xmath127 , the functions @xmath101 admit the asymptotic representation @xmath149+{\rm      f}_0(y;z)\right\ } \\ \nonumber & + &   \displaystyle\frac{\exp({\rm i}\sqrt{z}\rho)}{\sqrt{\rho } }       \left[a_l(z,\theta)+o(1)\right ] + { \rm f}_{1,l}(x , y;z ) \ ] ] with @xmath150 where @xmath151 in a parabolic neighborhood of the @xmath82-axis , the functions @xmath152 are also subjected to the asymptotic estimates @xmath153 where @xmath154 is an arbitrary fixed number smaller than unity , @xmath155 .    as to the domain @xmath103 , the leading asymptotic term of each of the functions @xmath152 for @xmath156 is a spherical wave @xmath75 with the amplitude @xmath157 being a differentiable function of the angle @xmath68 .
therefore , for @xmath158 the term @xmath152 in the r.h.s . of the formula  ( [ phipis ] ) can be added to the asymptotic term with a spherical wave preceding @xmath159 . in the domain @xmath106 , and hence , in a narrower domain @xmath103 the condition @xmath160 holds .
consequently , for @xmath156 , the faddeev components @xmath90 do obey the standard asymptotic conditions like ( [ asbcparts ] ) .
therefore , for any @xmath155 the dominant term of the asymptotics of the function @xmath161 , @xmath162 , in the domain @xmath163 reads as @xmath147 as @xmath148 .
this means that , for @xmath162 , it is always possible by solving the equations  ( [ fadpart ] ) to separate explicitly the elastic scattering amplitude @xmath70 and , thus , to construct the analytic continuation of the scattering matrix @xmath19 .    outside of the domain @xmath106 the numerical construction of @xmath19 by solving the faddeev differential equations is , in general , impossible since for @xmath163 and @xmath155 both functions @xmath164 and
@xmath165 , @xmath166 , include terms decreasing slower than @xmath74 as @xmath148 .      in the case of potentials with hard core , the partial faddeev differential equations for a system of three identical bosons at @xmath32 acquire the form @xmath167f_l(x , y)=\left\ {              \begin{array}{cl } -v(x)\psi_l(x , y ) , & x > c \\                      0 ,                   & x < c\ , , \end{array}\right.\ ] ] where @xmath168 , @xmath169 , is the core size .
the partial wave functions @xmath41 are expressed via faddeev partial components @xmath42 by the formulas  ( [ ftconn ] ) . the components @xmath42 satisfy the standard boundary conditions ( [ bcstandard ] ) .
the two - body central potential @xmath38 acts only beyond the core domain , i.e. only where @xmath170 .
we assume as before that @xmath38 falls off not slower than exponentially as and , hence , it satisfies the condition ( [ vexp ] ) for some @xmath171 and @xmath172 .
a main difference between the model with hard core and those with smooth potentials is that the functions @xmath42 in this model satisfy the auxiliary boundary conditions @xmath173 requiring that the wave functions @xmath174 vanish on the boundary @xmath175 of the core domain .
it can be shown that in fact the conditions ( [ bccorepart ] ) force the wave functions ( [ ftconn ] ) to vanish also inside the core domain at all energies except for a certain countable set of real values of @xmath20 ( see ref .
@xcite and references therein ) .
asymptotic conditions for the partial faddeev components @xmath90 of the @xmath58 scattering wave functions as @xmath56 and/or @xmath57 are again of the form ( [ asbcparts ] ) .
the only difference is that the dimer wave function @xmath64 is considered as zero in the core domain , i.e. for @xmath176 .    in the hard - core model ,
all the assertions of sect .
[ holomorphydomains ] regarding the holomorphy domains of the functions @xmath101 and the scattering matrix @xmath19 in @xmath20 still hold true .
we have already noticed that the roots of @xmath19 in the physical sheet of energy @xmath20 plane correspond to the location of the three - body resonances in the unphysical sheet adjoining the physical sheet along the spectral interval @xmath21 . in the case under consideration , this statement is an immediate consequence of the unitarity of the scattering matrix @xmath17 for @xmath60 , @xmath177 , @xmath178 indeed , as we have established , the functions @xmath101 are holomorphic functions of @xmath98 .
since the boundary value problem  ( [ fadpart][asbcparts ] ) is uniquely solvable , one easily verifies that the boundary values and for these functions on the rims of the cut along are related to each other as @xmath179 since , on the one hand , their asymptotics  ( [ phipis ] ) as @xmath148 and/or @xmath127 has the same structure and , on the other hand , @xmath180 since @xmath181 consequently , @xmath182 and @xmath183 therefore , it follows from eq .
( [ s0unitar ] ) that for @xmath177 @xmath184^{-1 } \qquad\mbox{and}\qquad     { \rm s}_0(e-{\rm i}0)=[{\rm s}_0(e+{\rm i}0)]^{-1}\,.\ ] ] this means that the function @xmath185 is continued through the cut @xmath186 $ ] into the domain @xmath187 as @xmath188 . in a similar manner
, @xmath189 is continued into the domain @xmath190 , again as @xmath188 .
all this signifies that the scattering matrix @xmath19 admits analytic continuation at least into the domain @xmath106 of the unphysical energy sheet connected with the physical sheet by crossing the interval @xmath186 $ ] , the value of the continued function @xmath17 at @xmath158 in the unphysical sheet coinciding with the value of @xmath188 at the same @xmath20 but in the physical sheet .
recall that those points @xmath20 on unphysical sheets are called resonances where the analytically continued scattering matrix possesses poles .
the resonances @xmath20 with zero imaginary part @xmath191 and @xmath192 are called the virtual levels .
thus , we have here presented a simple proof of the fact that the resonances including the virtual levels corresponding to poles of the analytic continuation of the scattering matrix @xmath17 in the unphysical sheet connected with the physical one by crossing the spectral interval @xmath186 $ ] are the roots of this matrix in the physical sheet . at the same time , the poles of the function @xmath19 in the physical sheet correspond to bound states of the three - boson system under consideration .    concluding the subsection
, we note that it follows from eq .
( [ aadj ] ) that and , hence , @xmath193 for any @xmath158 .
this means that the roots of the function @xmath19 are situated symmetrically with respect to the real axis .
in the present work we make use of the faddeev equations ( [ fadpartcor ] ) considered together with the boundary conditions ( [ bcstandard ] ) , ( [ asbcparts ] ) and ( [ bccorepart ] ) to calculate the values of the @xmath0he@xmath6 scattering matrix @xmath19 in the physical sheet
. we search for the resonances including the virtual levels as roots of @xmath19 and for the bound - state energies as positions of poles of @xmath19 .
all the results presented below are obtained for the case @xmath194 .    in all our calculations ,
as the interatomic he  he - interaction we employed the widely used semiempirical potential hfd - b constructed by r.a.aziz and co - workers  @xcite .
this potential is of the form @xmath195f(\zeta )   \right \}\,\ ] ] where @xmath196 .
the function @xmath197 reads @xmath198 ^ 2 } ,      & \mbox{if   $ \zeta\leq d$ }    \cr 1 ,   & \mbox{if   $ \zeta >   d $ } \ , . } \ ] ] for completeness the parameters of the potential hfd - b are given in table  i.    the value of the parameter @xmath168 ( the core  diameter " of particles ) is chosen to be so small that its further decrease does not appreciably influence the dimer binding energy @xmath2 and the energy of the trimer ground state @xmath199 . unlike papers  @xcite , where @xmath168 was taken to be equal 0.7 , now we take @xmath200 .
we have found that such a value of @xmath168 provides at least six reliable figures of @xmath2 and three figures of @xmath199 .    since the statements of sect .
[ difffs ] are valid , generally speaking , only for the potentials decreasing not slower than exponentially , we cut off the potential hfd - b setting @xmath201 for @xmath87 . we have established that this cutoff for provides the same values of @xmath2 ( @xmath202mk ) , @xmath199 ( @xmath203k ) and phases @xmath204 which were obtained in our earlier calculations @xcite performed with the potential hfd - b .
comparison of these results with results of other researchers can be found in refs .
in all the calculations of the present work we take . note that if the formulas from sect .
[ difffs ] including the parameter @xmath86 are used for finite potentials , one should set @xmath205 .    before making numerical approximation of the system of equations ( [ bcstandard ] ) , ( [ fadpartcor ] ) , ( [ bccorepart ] ) at @xmath194 we rewrite it in terms of a new unknown function @xmath161
that is expressed via the faddeev component @xmath206 by the relation  ( [ flphil ] ) .
note that for @xmath98 the function @xmath161 is square integrable in @xmath34 ( see sect .
[ holomorphydomains ] ) .
therefore , this function is uniquely determined by the asymptotic condition @xmath207 that can be easily approximated and programmed .
one could , for instance , require @xmath208 at a sufficiently large @xmath209 and look for a numerical solution of the system ( [ bcstandard ] ) , ( [ fadpartcor ] ) , ( [ bccorepart ] ) satisfying this condition .
further , for @xmath158 , one could , going sufficiently far from @xmath209 into the domain of smaller ( but nevertheless , providing the asymptotics  ( [ phipis ] ) ) values of @xmath66 , separate the elastic scattering amplitude @xmath70 , putting , e.g. , @xmath210 , where the value of @xmath131 corresponds to the maximum of the function @xmath64 . such an approach is , however , not effective in view of a relatively slow decrease of the exponentials @xmath211 and @xmath212 as well as of the function @xmath213 in the energy domain of interest for us in @xmath106 . for a proper approximation of the condition  ( [ phias0 ] ) ,
very large values of @xmath209 are to be taken .
this is just a reason why one should take into account the asymptotics of the function @xmath161 as @xmath53 and/or @xmath148 .
though the asymptotic formula  ( [ asbcparts ] ) only holds for @xmath156 , we employ it also for @xmath214 . indeed , when @xmath214 , the leading term of the asymptotics of @xmath161 as @xmath148 and @xmath163 , @xmath155 , is given by the same expression ( see sect .
[ holomorphydomains ] ) as in eq .
( [ asbcparts ] ) . outside of the parabola @xmath163 , it suffices to require the condition ( [ phias0 ] ) to be satisfied .
the presence , in eq .
( [ asbcparts ] ) , of the spherical wave does not contradict this requirement .
therefore , the use of asymptotic condition ( [ asbcparts ] ) is justified even if @xmath214 .
when solving the boundary - value problem ( [ bcstandard ] ) , ( [ asbcparts ] ) , ( [ fadpartcor ] ) , ( [ bccorepart ] ) written in terms of the function @xmath215 , we carry out its finite - difference approximation in polar coordinates @xmath66 and @xmath68 .
the grid is chosen in such a way that the points of intersection of arcs @xmath216 , @xmath217 , and rays @xmath218 , @xmath219 , with the line @xmath175 turn out automatically to be its knots .
the @xmath220 points are chosen according to the formulas @xmath221 where @xmath222 stands for the number of arcs inside the core domain and @xmath223 the nonlinear monotonously increasing function @xmath224 , @xmath225 , satisfying the conditions @xmath226 and @xmath227 is chosen in the form @xmath228\\
\alpha_1\tau+\tau^\nu & , & \tau\in(\tau_0,1 ]          \end{array }          \right .. \ ] ] the values of @xmath229 , @xmath230 and @xmath231 , @xmath232 are determined via @xmath233 and @xmath154 from the continuity condition for @xmath224 and its derivative at the point @xmath233 . as a rule
, we took values of @xmath233 within 0.1 and 0.2 .
the value of the power @xmath154 depends on the cutoff radius @xmath234@xmath235 its range being within 2 and 4 in our calculations .    the knots @xmath236 at @xmath237
are taken according to @xmath238 .
the rest knots @xmath236 , @xmath239 are chosen equidistantly .
such a choice of the grid is prescribed by the need to have a higher density of points in the domain where the functions @xmath101 are most rapidly changing , i.e. for small values of @xmath66 and/or @xmath131 and lower in the asymptotic domain . in this work
, we used the grids of dimension .
the number of the last arc knots in @xmath68 lying in the core domain was usually equal to @xmath240 .
the finite - difference approximation of the integro - differential equations ( [ fadpartcor ] ) and boundary conditions ( [ bcstandard ] ) , ( [ bccorepart ] ) for @xmath194 reduces the problem to a system of @xmath241 linear algebraic equations .
the finite - difference equations corresponding to the arc @xmath242 include initially the values of the unknown function @xmath243 from the arc @xmath244 . to eliminate them , we express these values through the values of @xmath243 on the arcs @xmath242 and @xmath245 by using the asymptotic formula ( [ asbcparts ] ) , just in the manner described in the concluding part of appendix a of ref .  @xcite . in @xcite ,
this approach was only used for computing the energies of bound states .
now we extend it also on the scattering problem .
( note that the formulas ( a10 ) and ( a11 ) in @xcite related to the described approach contain misprints .
the values @xmath246 in these formulas should be replaced with inverse values @xmath247 . )
the matrix of the resultant system of equations has a block - three - diagonal form ( see ref .
@xcite , appendix a ) .
every block has the dimension @xmath248 and consists of the coefficients standing at unknown values of the function @xmath243 in the grid knots belonging to a certain arc @xmath216 .
the main diagonal of the matrix consists of @xmath249 such blocks .
in contrast to @xcite , in the present paper we solve the block - three - diagonal algebraic system on the basis of the matrix sweep method .
this allows us to dispense with writing the system matrix on the hard drive and to carry out all the operations related to its inversion immediately in ram .
besides , the matrix sweep method reduces almost by one order the computer time required for computations on the grids of the same dimensions as in  @xcite .
we searched for the resonances ( roots of the function @xmath17 on the physical sheet ) and bound - state energies ( roots of the function @xmath188 for real @xmath250 ) of the helium trimer by using the complex version of the secant method . within this method ,
the approximation @xmath251 to a root of a holomorphic function @xmath252 is constructed from the two previous approximations @xmath253 and @xmath254 according to the formula @xmath255    as the relationship ( [ s0sym ] ) implies the symmetry of properties of the scattering matrix @xmath19 with respect to the real axis , we performed all the calculations for @xmath19 only for @xmath256 ( except the tests of the code ) .
we start with a study of graph surfaces of the real and imaginary parts of the scattering matrix @xmath19 in the domain of its holomorphy @xmath257 .
the root lines of the functions @xmath258 and @xmath259 obtained in the case of the grid parameters and  are depicted in fig .  1 .
both resonances ( roots of @xmath19 ) and bound - state energies ( poles of @xmath19 ) of the @xmath0he trimer are associated with the intersection points of the curves @xmath260 and @xmath261 . when the roots or poles are simple , these curves intersect each other at the right angle .
note that for real @xmath262 the function @xmath19 is real and , thus , @xmath261 . in fig .  1 , along with the root lines we also plot the boundaries of the domains @xmath106 , @xmath103 and @xmath95 .
one can observe that a  good " domain @xmath106 includes none of the points of intersection of the root lines @xmath260 and @xmath263 .
nevertheless , as we will see below , the going beyond this domain is of an interest , even though the asymptotic formula ( [ asbcparts ] ) is not valid for @xmath264 and the function @xmath19 calculated there can not be interpreted as the scattering matrix .
the caption for fig .
1 points out positions of the four  resonances " , the roots of @xmath19 , found immediately beyond the boundary of the domain @xmath106 .
as one could expect , the values of the function @xmath19 at @xmath264 and positions of its roots in @xmath265 turn out to be unstable and strongly depend on the value of the cutoff radius @xmath266 , whereas the dependence on the number of knots is weak . in particular , for @xmath267 , a ( quasi)resonance , closest to the real axis , is situated at the point if at the point if and at the point if the same ( quasi)resonance in fig .  1 ( calculated for ) is situated at the point @xmath268  mk . if @xmath269 is fixed , the increase of @xmath266 up to 800 shifts this point to the point @xmath270mk .
all the aforesaid regarding the instability of the function @xmath19 values and positions of its roots beyond the domain @xmath106 bears no relation to its pole at the point , corresponding to a trimer excited - state energy , even though this energy does not belong to @xmath106 .
the point is that the position of the pole of @xmath17 is only determined by the position of the root of the determinant of the linear algebraic system we solve , whereas the inhomogeneous term of the system plays no role .
therefore , the search for the poles of the grid function @xmath19 is equivalent to the search for the binding energies of the trimer .
the grids we have used turn out to be quite sufficient for this purpose . the convergence of our results for @xmath13 with respect to the parameters @xmath271 and their accuracy can be judged from the values of the difference @xmath272 obtained with different grids and shown in table  ii .
we would like to stress that we do not consider the roots of function @xmath19 drawn in fig .  1 as genuine resonances for the @xmath0he@xmath6 system since they are situated beyond the domain @xmath106 where the applicability of our method is proved .
we should rather consider them as artifacts of the method .
however it is remarkable that the  true " ( i.e. , getting inside @xmath106 ) virtual levels and then the energies of the excited ( efimov ) states appear just due to these ( quasi)resonances when the potential @xmath38 is weakened . this is the object of our further consideration .
following  @xcite , instead of the initial potential @xmath273 , we will consider the potentials @xmath274 to establish the mechanism of formation of new excited states in the @xmath0he trimer , we have first calculated the scattering matrix @xmath19 for @xmath15 . in table
iii for some values of @xmath7 from the interval between 0.995 and 0.975 , we present the positions of roots and poles of @xmath19 , we have obtained at real @xmath275 . we have found that for a value of @xmath7 slightly smaller than @xmath276 , the ( quasi)resonance closest to the real axis ( see fig .
1 ) gets on it and transforms into a virtual level ( the root of @xmath19 ) of the second order corresponding to the energy value where the graph of @xmath19 , @xmath277 , @xmath250 , is tangent to the axis @xmath20 .
this virtual level is preceded by the ( quasi)resonances mk for @xmath278 and mk for @xmath279 .
the originating virtual level is of the second order since simultaneously with the root of the function @xmath17 , also the conjugate root of this function gets on the real axis . with a subsequent decrease of @xmath7 the virtual level of the second order splits into a pair of the virtual levels @xmath280 and @xmath281 , @xmath282 of the first order which move in opposite directions . a characteristic behavior of the scattering matrix @xmath19 when resonances transform into virtual levels is shown in fig .
the virtual level @xmath281 moves towards the threshold @xmath2 and `` collides '' with it at @xmath283 . for @xmath284
the function @xmath19 has no longer the root corresponding to @xmath281 . instead of the root
, it acquires a new pole corresponding to the second excited state of the trimer with the energy @xmath14 . note that though the virtual levels @xmath280 and @xmath281 appear beyond the domain @xmath106 , already at @xmath285 the point @xmath281 turns out to be inside this domain .
therefore , it should be considered as a `` true '' virtual level of the trimer .
we expect that the subsequent efimov levels originate from the virtual levels just according to the same scheme as the level @xmath14 does .
the other purpose of the present investigation is to determine the mechanism of disappearance of the excited state of the helium trimer when the two - body interactions become stronger owing to the increasing coupling constant @xmath8 .
it turned out that this disappearance proceeds just according to the scheme of the formation of new excited states ; only the order of occurring events is inverse .
the results of our computations of the energy @xmath13 when @xmath7 changes from 1.05 to 1.17 are given in table  iv .
in the interval between @xmath286 and @xmath287 there occurs a  jump " of the level @xmath13 on the unphysical sheet and it transforms from the pole of the function @xmath17 into its root , @xmath288 , corresponding to the trimer virtual level .
the results of calculation of this virtual level where @xmath7 changes from 1.18 to 1.5 are presented in table  v. for all the values of @xmath7 presented in tables  iv and  v , the dimer possesses an only bound state .
we have found that the first excited state of the dimer appears only at @xmath289 .
note that in the case of finite potentials the geometric characteristics of the domain @xmath106 where the function @xmath19 can be calculated reliably , are only determined by the value of @xmath290 ( see formula  ( [ pipsi ] ) for @xmath205 ) .
when @xmath290 increases , the domain @xmath106 is enlarged .
it is easy to check that the energies of the excited state level @xmath11 and of the virtual level @xmath291 given in tables  iv and  v belong to the corresponding domains @xmath292 . for @xmath293 ,
this results in a weak dependence of the calculated values of @xmath11 and @xmath291 on the parameters @xmath294 , @xmath249 and ( this is especially important ) on the parameter @xmath209 .    in essential , we chose the values of the cutoff hyperradius @xmath209 given in tables  iii
v from the scaling considerations . as a matter of fact
, we took the value of @xmath209 following the formula @xmath295 where the  constant " @xmath296 corresponds to an appropriate choice of @xmath209 at @xmath297 .
it has been established in  @xcite that such a choice is ensured if @xmath298400600 . in determining the values of @xmath299 ,
indicated in table  iii , we followed the formula  ( [ rhomax ] ) literally .
as the  constant " @xmath300 , we took its value corresponding to the base value of @xmath301600 .
the values of @xmath302 presented in tables  iv and v correspond to the choice of @xmath303 in the interval within 400 and 800 .
all the results presented in tables  iii  v have been obtained with the grids parameters @xmath269 .
the authors are grateful to prof .
v.b.belyaev and prof .
h.toki for help and assistance in calculations at the supercomputer of the research center for nuclear physics of the osaka university , japan . also , the autors thank mrs .
t.dumbrais for her help in translation of the text into english .
one of the authors ( a.k.m . ) is much indebted to prof .
w.sandhas for his hospitality at the universitt bonn .
the support of this work by the deutsche forschungsgemeinschaft and russian foundation for basic research is gratefully acknowledged .
99 rama krishna m.v .
, whaley k.b .
rev . lett .
p.  1126 .
lehman k.k . ,
scoles g. // science .
p.  2065 .
grebenev s. , toennies j.p .
, vilesov a.f .
// science .
p.  2083 .
efimov v. //
luo f. , mcbane g.c .
, kim g. , giese c.f . , gentry w.r . //
j. chem . phys .
p.  3564 .
luo f. , giese c.f .
, gentry w. r. // j. chem . phys . 1996 . v.  104 .
p.  1151 .
schllkopf w. , toennies j.p .
// science .
1994 . v.  266 .
p.  1345 .
aziz r.a . ,
nain v.p.s .
, carley j.s . , taylor w.l .
, mcconville g.t .
// j. chem . phys .
1979 . v.  79 .
p.  4330 .
aziz r.a . ,
mccourt f.r.w .
, wong c.c.k .
1987 . v.  61 .
p.  1487 .
, slaman m.j .
// j. chem .
1991 . v.  94 .
p.  8047 .
tang k.t . ,
toennies j.p . , yiu c.l . //
1995 . v.  74 .
p.  1546 .
mcmillan w.l .
pandharipande v.r . ,
zabolitzky j.g . ,
pieper s.c . ,
wiringa r.b . , helmbrecht u.
1983 . v.  50 .
p.  1676 .
usmani n. , fantoni s. , pandharipande v.r .
p.  6123 .
pieper s.c . ,
wiringa r.b .
, pandharipande v.r .
// phys . rev . b. 1985 . v.  32 .
p.  r3341 .
rick s.w . ,
lynch d.l .
, doll j.d .
// j. chem . phys .
1991 . v.  95 .
p.  3506 .
levinger j.s .
// yadernaya fizika .
, lin c.d . , greene c.h . //
a. 1996 . v.  54 .
nielsen e. , fedorov d.v .
, jensen a.s .
// lanl e - print physics/9806020 .
nakaichi - maeda s. , lim t.k .
a. 1983 . v.  28 .
cornelius th .
, glckle w.
// j. chem . phys .
1986 . v.  85 .
p.  3906 .
carbonell j. , gignoux c. , merkuriev s.p .
// few  body systems
schllkopf w. , toennies j.p . //
j. chem . phys .
1996 . v.  104 .
p.  1155 .
kolganova e.a . ,
motovilov a.k .
, sofianos s.a .
a. 1997 . v.  56 .
p.  1686 .
motovilov a.k . ,
sofianos s.a .
, kolganova e.a .
phys . lett .
1997 . v.  275 .
kolganova e.a . ,
motovilov a.k .
, sofianos s.a . //
j.  phys .
p.  1279 .
fedichev p.o .
, reynolds m.w .
, shlyapnikov g.v .
1996 . v.  77 .
p.  2921 .
merkuriev s.p .
, motovilov a.k .
1983 . v.  7 .
merkuriev s.p . ,
motovilov a.k .
, yakovlev s.l .
kolganova  e.a .
, motovilov  a.k .
// phys . atom .
faddeev l.d .
, merkuriev s.p .
quantum scattering theory for several particle systems .
kluwer academic publishers , doderecht , 1993 .
merkuriev s.p . ,
gignoux c. , laverne a. // ann . phys .
1976 . v.  99 .
motovilov  a.k .
1996 . v.  107 .
motovilov a.k .