Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: the last decade and a half have seen enormous progress in the description of light nuclei as collections of interacting nucleons with the same properties as in vacuum @xcite . after the formulation of potentials that fit nucleon - nucleon scattering data with high accuracy , _ ab initio _ calculations of nuclear structure demonstrated that the energy spectra of nuclei small enough for converged calculations ( mass number @xmath3 ) can be understood as arising from the vacuum nucleon - nucleon interaction @xcite . the addition of a three - nucleon interaction with only two to four parameters produces quantitative agreement with experiment for the nuclear binding energies and spin - orbit splittings , including the correct ground state ( @xmath4 instead of @xmath5 ) for @xmath6b @xcite .    for nuclei with @xmath7 , methods using correlated hyperspherical harmonic bases or faddeev and similar formulations have solved bound - state , scattering , and reaction problems quite successfully @xcite . the _ ab initio _ methods that have been developed for @xmath8 nuclei are suited mainly to treatment of bound states , but there has been significant progress on unbound states in recent years . the green s function monte carlo ( gfmc ) method has been formulated as a particle - in - a - box method and used to compute phase shifts in the @xmath9he system with a full realistic hamiltonian @xcite . the no - core shell model ( ncsm ) has been merged with the resonating group method ( rgm ) to produce phase shifts and reaction cross sections in several systems using effective two - body forces without explicit three - body terms @xcite . the coupled - cluster method has been combined with a gamow shell model basis to compute complex energies @xmath10 of resonances and thus their excitation energies @xmath11 and total widths @xmath12 . it has been applied to resonances of he isotopes and @xmath13 nuclei , but it has not been used to produce cross sections @xcite . all of these approaches to scattering are computationally challenging , and significant human and computer effort is required for each individual system .    even before the development of true scattering / reaction calculations in @xmath8 nuclei , there were accurate energy calculations for many unbound resonance states @xcite , produced by approximating them with bound ( i.e. , square - integrable ) wave functions and real energies . this approach is successful for states sufficiently narrow that important features of their structure can be accommodated within the model space . one way to understand these calculations is as approximations to gamow s decaying complex - energy states @xcite with energies near the real axis ( i.e. , small widths ) . quantum monte carlo ( qmc ) calculations of square - integrable `` pseudobound '' states begin with a variational monte carlo ( vmc ) calculation , in which a complicated but closed - form wave function containing large amounts of correlation is produced by minimizing the energy expectation value with respect to many variational parameters @xcite . the computational effort in the vmc method lies in computing energy expectation values by monte carlo integration over all particle coordinates , and a square - integrable wave function is necessary for both energy and normalization integrals to converge . the variational ansatz incorporates square integrability through particle correlations that decay exponentially at large separation , and reasonable energies result because the long - range tails of the true resonant wave function have small amplitudes relative to the `` interaction region '' where all nucleons interact . the second step in application of the qmc method to nuclei is a gfmc calculation . this step also requires square - integrable wave functions for evaluation of matrix elements . gfmc takes the vmc wave function as a starting point and evolves the schrdinger equation for imaginary values of the time @xmath14 through a series of small steps . this evolution filters high - energy contamination out of the wave function and for large @xmath15 leaves behind the lowest - energy eigenstate contained in the vmc starting point . for bound or narrow states , this procedure `` converges '' at large @xmath15 to a unique energy @xcite . for broad states the gfmc propagation drifts slowly toward the lowest available threshold ; even in very long gfmc propagations , the curve of energy _ versus _ @xmath15 for such states fails to flatten out at a `` converged '' energy . presumably this also occurs for the narrow states , but with a decline too slow to have been noticed in calculations so far . because the calculated energies of narrow resonances in the qmc methods are believed to be accurate , it is natural to ask whether widths could be extracted from these calculations , avoiding the considerable effort of explicit scattering / reaction calculations . attempts to correlate the width with the rate of decline of the energy , @xmath16 , in the late - time gfmc propagation failed , so the rate may depend more on the starting wave function than on the physical width  @xcite .    even in explicit scattering calculations , some widths are effectively too narrow to resolve by qmc methods . direct qmc calculations of widths so far consist of energy calculations for a series of specified boundary conditions @xcite , so the smallest width that can be computed must be larger than the energy resolution of the method . as an extreme example , the ground state of @xmath2b has a total energy of @xmath17 mev and a width of 0.5 kev , roughly 0.01% of the total . since the typical precision claimed for gfmc energy calculations is 1% , this width will not be amenable to calculation in this way .    here i present a method to extract approximate widths from qmc wave functions . the basis of this approach is a relation between the partial width in a specified breakup channel and an integral over the `` interaction region '' where all of the nucleons are close to each other . this relation has been used in the literature , though not ( to my knowledge ) in explicit application to many - body wave functions . it has been used particularly in models of alpha @xcite and proton @xcite decays of heavy nuclei in order to avoid direct integration of the schrdinger equation to impractically large radii . it is also closely related to green s function formulations of scattering theory . there has been considerable interest in applying such formulations to _ ab initio _ calculations in the recent literature : they have been used to extract asymptotic normalization coefficients ( ancs ) from many - body bound states computed in various approximations @xcite and are coming into wider use for scattering problems @xcite . preparation for future reaction calculations of that kind is a primary motivation of the present work . although the integral relation can be used to compute the partial width in any two - body decay channel if wave functions of the parent and daughter nuclei are known , i confine my consideration in the present study to nucleon emission . i compute widths for all narrow ( @xmath18 mev ) one - nucleon decays to bound states in @xmath19 nuclei . i also present results for several specific broad states for which better theoretical information concerning widths would be useful and for some with unbound final states that are well - represented by pseudobound wave functions . the integral relation applied here can be used with either gfmc or vmc wave functions , extending even to use of the same computer routines . application to gfmc requires more computation and additional bookkeeping , so i have chosen in this initial study to use only vmc wave functions . this work represents the extension to unbound states of the methods presented in ref . the integral relation used here may be used with other many - body methods , but it is particularly well - suited to qmc wave functions because of its formulation as a short - range integral and because of the transparent treatment of fermion antisymmetry in the qmc methods . the remainder of this paper is organized as follows : in sec . [ sec : derivation ] i motivate and define the integral method , and i discuss its accuracy and its connection to overlap calculations . in sec . [ sec : vmc ] i describe the application of the integral method to vmc wave functions . in sec . [ sec : results ] i present the results of calculations for specific states , compare computed and experimental widths , and consider how accurately widths may be computed from spectroscopic factors . along the way i present new ( and likely improved ) overlap calculations of vmc wave functions that may prove useful for treatment of transfer reactions . finally in sec . [ sec : conclusion ] i summarize my results and briefly mention future directions for this work . consider a many - body wave function @xmath20 at energy @xmath21 above the threshold for breakup into clusters 1 and 2 that have wave functions @xmath22 and @xmath23 respectively and no internal angular momentum . assume also for the moment that only this breakup channel is open . since @xmath22 and @xmath23 are spinless , the orbital angular momentum @xmath24 of their motion is equal to the total angular momentum @xmath25 of @xmath20 . given appropriate boundary conditions , @xmath20 at large separations @xmath26 of the clusters is a linear combination of an incoming wave proportional to @xmath27 and an outgoing wave proportional to @xmath28 , each normalized to probability flux @xmath29 in its appropriate direction at @xmath30 . ( for all special functions , i follow the conventions of refs . @xcite and @xcite . ) these functions depend on the energy @xmath21 through the wavenumber @xmath31 , where @xmath32 is the reduced mass of the clusters , and through the sommerfeld parameter @xmath33 , where @xmath34 and @xmath35 are the charge numbers of the clusters . @xmath36 and @xmath37 solve the radial coulomb - schrdinger equation @xmath38 with @xmath39 and outward and inward flux , respectively , at @xmath30 . in terms of @xmath37 and @xmath36 , @xmath20 ( assumed to have angular - momentum projection @xmath40 , omitted from subsequent expressions ) is , @xmath41/r_{12}.\ ] ] where @xmath42 is a spherical harmonic . in this single - channel case with both incoming and outgoing waves , the normalization @xmath43 is arbitrary . if @xmath44 and @xmath45 are real , conservation of probability guarantees that the effect of scattering is to multiply the outgoing wave @xmath36 by a complex phase factor relative to the incoming wave @xmath37 , so that the function @xmath46 may be written in terms of a real phase shift @xmath47 as @xmath48 . the function @xmath46 is the single - channel case of the @xmath49-matrix , which gives the amplitude and phase of outgoing waves in terms of specified ingoing amplitudes , and it determines the cross section uniquely . resonances in the scattering of clusters 1 and 2 occur at real energies near poles of @xmath46 and produce peaks in the scattering cross section where @xmath50 around some resonance with energy @xmath11 and width @xmath12 . @xmath49-matrix poles in general occur at complex @xmath44 ( and thus complex @xmath21 ) . it is evident from eq . ( [ eq : inandout ] ) that at such a pole the wave function takes the form @xmath51 with only outgoing flux . this is the same sort of decaying complex - energy state originally formulated in gamow s treatment of alpha decay @xcite . the probability flux out of this state is @xmath52 , so the normalization constant is no longer arbitrary but carries information about the size of the outgoing flux relative to the total wave function . if @xmath20 has been normalized to unit total probability in some finite region , its lifetime is inversely proportional to the outward probability flux . by writing @xmath20 in the form @xmath53^{-1 } i_l -      o_l\right]/r_{12},\ ] ] and defining @xmath20 to have unit norm in some regularization scheme , it can be shown that the residue of the @xmath49-matrix pole is proportional both to the squared normalization constant @xmath54 and to the imaginary part of the pole energy @xcite . the relations among @xmath54 , the lifetime , and the pole location imply that the width @xmath55 . ( derivations that deal rigorously with the normalization of @xmath20 may be found in refs . @xcite . ) the width may thus be computed from the normalization constant @xmath56 . physically realized systems have real energies , so formulation of the width in terms of complex - energy gamow states is often inconvenient . for the qmc methods , not only is the energy real , but the wave functions are stationary waves , sums of incoming and outgoing waves with zero total flux in each channel . again assuming a single open channel of given angular momentum , a standing - wave solution is asymptotically @xmath57/r_{12},\ ] ] for some different @xmath43 than before . @xmath58 is the regular coulomb function that satisfies @xmath59 , @xmath60 satisfies the wronskian relation @xmath61 and the two are related to @xmath36 and @xmath37 by @xmath62 the function @xmath63 is the single - channel case of the @xmath64-matrix . ( [ eq : inandout ] ) , ( [ eq : outgoing ] ) , and ( [ eq : incoming ] ) combine to give @xmath65 so that @xmath66 . ( [ eq : standing - wave ] ) may then be written as @xmath67/r_{12}\ ] ] for some @xmath68 determined by an appropriate normalization of @xmath20 . at a resonance , @xmath69 and @xmath70 . this corresponds to a pole of @xmath63 on the positive real axis , so that @xmath71/r_{12}.\end{aligned}\ ] ] the standing - wave solution at resonance is thus the sum of equal inward and outward fluxes , each of magnitude @xmath72 . this is the rate at which @xmath20 decays through the outgoing wave and is replenished by the ingoing wave . in analogy with the complex - energy case , it can be shown that the residue of the @xmath64-matrix at the resonance pole is proportional to @xmath73 for a wave function normalized to unity in a finite region @xcite , up to corrections arising from the choice of normalization volume . these corrections are small as long as the wave function has a much smaller amplitude in the asymptotic region than in the interaction region @xcite , as is expected for a long - lived resonance state . comparing eqs . ( [ eq : s - matrix - pole ] ) and ( [ eq : gl - norm ] ) , it is unsurprising that the result of a rigorous derivation is @xmath74 . i have sketched the considerations that lead to the connection between resonance width and wave function asymptotic normalization for the case of a single open channel and spinless daughter nuclei . these considerations carry over directly to _ ab initio _ wave functions , but it is necessary to account for complications absent from the toy model . most complications of the _ ab initio _  case amount to additional bookkeeping implied by multiple final - state channels , antisymmetry of the wave function , and daughters with nonzero angular momenta . multiple types of decay products may be emitted , and their non - zero internal angular momenta may couple in multiple ways . the simple right - hand sides of eqs . ( [ eq : inandout ] ) and ( [ eq : standing - wave ] ) are replaced by sums over all breakup and @xmath24 channels of a given total angular momentum @xmath25 and parity @xmath75 . the formalism to describe multi - channel wave functions in this case may be found in many treatments of reaction theory ( e.g. refs . the coefficient @xmath46 in eq . ( [ eq : inandout ] ) is replaced by a matrix @xmath76 that gives the outgoing flux in channel @xmath77 produced by unit incoming flux in channel @xmath78 , and @xmath63 in eq . ( [ eq : standing - wave ] ) is similarly replaced by a matrix @xmath79 that gives the irregular - function amplitude in channel @xmath77 produced by unit regular - function amplitude in channel @xmath78 . if the initial state in a reaction problem has amplitude @xmath80 in channel @xmath77 , then @xmath81_{j_a}\right]_{j } \left\{x_{a}f_{l_a}(\eta_a , k_ar_a ) + y_a g_{l_a}(\eta_a , k_ar_a)\right\}/r_a\ .\ ] ] in this expression , channel @xmath77 is characterized by daughter nuclei with wave functions @xmath82 and @xmath83 , wave number @xmath84 , sommerfeld parameter @xmath85 , and daughter separation @xmath86 . the orbital angular momentum of this channel is @xmath87 , and the square brackets denote coupling of the daughter angular momenta @xmath88 and @xmath87 first to `` @xmath89-coupled '' angular momentum @xmath90 and then with @xmath88 to form total angular momentum @xmath25 . ( the use of this coupling anticipates its later utility in defining channels in qmc wave functions . ) the antisymmetrization operator @xmath91 carries out an antisymmetric sum over all partitions of nucleons into daughters @xmath92 and @xmath93 with mass numbers @xmath94 and @xmath95 and multiplies by the normalization @xmath96 . the index @xmath77 is taken to specify the daughter nuclei as well as the channel quantum numbers @xmath90 , @xmath97 , @xmath88 , @xmath87 , and @xmath98 . away from resonance , the irregular - function amplitudes are given by @xmath99 . at a pole of @xmath79 , all asymptotic channels that are coupled to the resonance have wave functions proportional to @xmath100 . then the corresponding @xmath80 are irrelevant and may be set to zero . if @xmath20 at a pole is normalized to unit probability within some finite volume that includes the whole interaction region , the residue of @xmath101 is proportional to @xmath102 @xcite . the partial width of the resonance in channel @xmath77 is proportional to this residue , just as it was in the single - channel case discussed above @xcite . it can be shown that the imaginary part of the corresponding @xmath49-matrix pole energy is proportional to the sum of these partial widths @xcite , so the pole residues give the partial widths , whereas the @xmath49-matrix pole location gives only the total width . the channel radial functions defined by the summed terms in eq . ( [ eq : kmatrix - multichannel ] ) can be isolated from @xmath20 by projecting it onto the channel functions @xmath103 where the channel function in a given partition @xmath104 of the nucleons into clusters is @xmath105_{j_a}\right]_{j}\ , \\ { \bm\xi}_{ai } & = & { \bm\xi}_{ai}^p \equiv \left\{\mathbf{r}_j - \mathbf{r}_{ai}\right\},~j\in~ai\,.\end{aligned}\ ] ] the @xmath106 are the internal coordinates of @xmath107 , written in terms of nucleon coordinates in the specified partition @xmath104 as differences @xmath108 of the coordinates @xmath109 of nucleons within daughter @xmath110 and the center of mass of daughter @xmath110 . the daughter centers of mass @xmath111 are related to their separation @xmath112 by @xmath113 using bracket notation to indicate inner product in nucleon spin - isospin space and integration over all nucleon coordinates @xmath114 , @xmath115_{j_a}\right]_{j}^\dag      \frac{\delta(r_a - r)}{r_a^2}\psi d^{3a}r,\ ] ] which for an exact solution of the hamiltonian gives the overlap function @xmath116 it is this function that multiplies the daughter wave functions in eq . ( [ eq : kmatrix - multichannel ] ) . ( formally the integral extends over only @xmath117 coordinates , because center - of - mass motion is irrelevant . the wave functions used in nuclear qmc calculations are translationally invariant , allowing the integral to be written as extending over all coordinates ; the effect of center - of - mass motion cancels out when dividing computed quantities by the wave function normalization . ) ( [ eq : gl - norm ] ) suggests that the partial width of a real - energy resonance state may be found by first computing @xmath20 at resonance and then projecting onto the desired channel and dividing by @xmath118 to obtain the partial width @xmath119 this approach is useful when @xmath20 can be computed accurately in the large-@xmath120 region ( as in ref . however , it often happens that @xmath20 is computed accurately in the interaction region , but computation in the asymptotic region described by eq . ( [ eq : kmatrix - multichannel ] ) is difficult or inconvenient . a more robust approach proceeds through a green s function formalism . one begins with the schrdinger equation , @xmath121 with @xmath122 the hamiltonian operator and @xmath21 the energy . to isolate a particular channel @xmath77 , as defined above , it is useful to choose a partition @xmath104 of nucleons into the daughter nuclei in that channel and write @xmath122 as a sum of parts internal to the daughters and parts that depend on their relative motion . working in the center - of - mass frame , @xmath123 where @xmath124 is the kinetic energy of relative motion of the daughters . the operator @xmath125 is the part of the hamiltonian ( kinetic plus potential ) involving only the coordinates @xmath126 within daughter @xmath110 . @xmath127 contains the remaining terms of the potential , consisting of interactions between nucleons in daughter 1 and nucleons in daughter 2 : @xmath128 this is _ not _ an effective interaction but rather the sum of all terms of the many - body potential ( with two - body terms @xmath129 and three - body terms @xmath130 ) linking the two daughters . the point - coulomb interaction between daughters , @xmath131 , is added and subtracted for reasons that will become apparent later ; it is zero if one of the daughters is a neutron . the energy @xmath21 may be similarly broken up into the sum of the daughter internal energies @xmath132 and the channel energy @xmath133 : @xmath134 then eq . ( [ eq : schroedinger ] ) becomes @xmath135 rearranging terms and multiplying by the operator @xmath136^{-1}$ ] , we find @xmath137^{-1 }   \left({\ensuremath{u_\mathrm{rel}}}^{a , p } - v^a_c\right)\psi\\   & & - \left[t_{a , p } + v^a_c - e_a\right]^{-1 }     \left(h_1^{a , p } + h_2^{a , p } - e_{a1 } - e_{a2}\right)\psi\ . \nonumber\end{aligned}\ ] ]    the next step is to project @xmath20 onto channel @xmath77 , as in eq . ( [ eq : projection ] ) . projection of eq . ( [ eq : broken - up ] ) onto a channel function @xmath138 of eq . ( [ eq : channel - single - partition ] ) gives @xmath139^\dag\psi     & = &     -\left[\tilde{\phi}_{a , p}(\mathbf{r}_a^\prime)\right]^\dag       \left[t_{a , p } + v^a_c - e_a\right]^{-1 }       \left({\ensuremath{u_\mathrm{rel}}}^{a , p } - v^a_c\right)\psi\\       & & - \left[\tilde\phi_{a , p}(\mathbf{r}_a^\prime)\right]^\dag        \left[t_{a , p } + v^a_c - e_a\right]^{-1 }        \left(h^{a , p}_1 + h^{a , p}_2 - e_{a1 } - e_{a2}\right)\psi \nonumber\end{aligned}\ ] ] ( abbreviating @xmath140 on the right side by omitting the cluster internal coordinates ) , but since @xmath141 the second line of eq . ( [ eq : channel - projection ] ) is zero , and @xmath142^\dag\psi     = -\left[\tilde\phi_{a , p}(\mathbf{r}_a^\prime)\right]^\dag       \left[t_{a , p } + v^a_c - e_a\right]^{-1 }        \left({\ensuremath{u_\mathrm{rel}}}^{a , p } - v^a_c\right)\psi\ .\ ] ]    now , the operator @xmath143^{-1}$ ] takes functions of nucleon coordinates @xmath144 to functions of nucleon coordinates @xmath145 with different values of the separation @xmath112 but with the @xmath108 untouched . its application to a function @xmath146 and projection onto a second function @xmath147 may be written as an integral over a green s function that contains a product of coulomb wave functions : @xmath148^{-1}\phi_1(\mathbf{r } ) = \frac{2\mu}{\hbar^2k_a}\phi_2^\dag(\mathbf{r}^\prime ) \int d^3 r_a\ \frac{f_{l_a}\!(\eta_a , k_ar_<)g_{l_a}\!(\eta_a , k_ar_>)}{r _ < r _ > } y_{l_a}(\mathbf{\hat{r}}_a^\prime)y_{l_a}^\ast(\mathbf{\hat{r}}_a ) \phi_1(\mathbf{r}),\ ] ] following the usual notation that @xmath149 denotes the smaller of @xmath86 ( a jacobi coordinate specified for partition @xmath104 by @xmath144 ) and @xmath150 ( specified by @xmath104 and @xmath145 ) , while @xmath151 denotes the larger of @xmath86 and @xmath150 . rewriting eq . ( [ eq : direct - overlap ] ) in terms of eq . ( [ eq : better - channel - projection ] ) and applying eq . ( [ eq : greens - function ] ) , the result of integration over @xmath152 and antisymmetrization over partitions @xmath104 is @xmath153^\dag         \left({\ensuremath{u_\mathrm{rel}}}^{a , p } - v^a_c\right)\psi\ , d^{3a}r    \right.\\\nonumber & &    \left . + \frac{f_{l_a}\!(\eta_a , k_ar)}{r }    \left\{b_{a,\infty } + \mathcal{a}_a\int_{r_a > r }    \frac{g_{l_a}\!(\eta_a , k_ar_a)}{r_a } \left[\tilde\phi_{a , p}(\mathbf{r}_a)\right]^\dag         \left({\ensuremath{u_\mathrm{rel}}}^{a , p } - v^a_c\right)\psi\ , d^{3a}r    \right\ }    \right]\,.\end{aligned}\ ] ] the constant of integration @xmath154 can be determined from @xmath20 as in ref . @xcite ; as discussed below , @xmath155 at resonance . finally , i denote the integrals in eq . ( [ eq : full - integral - relation ] ) by @xmath156^\dag    \left({\ensuremath{u_\mathrm{rel}}}^{a , p } - v^a_c\right)\psi\ , d^{3a}r\\ \label{eq : g - norm - integral }    c_a(r ) & = & \frac{2\mu}{\hbar^2k_a } \mathcal{a}_a\int_{r_a < r }    \frac{f_{l_a}\!(\eta_a , k_ar_a)}{r_a } \left[\tilde\phi_{a , p}(\mathbf{r}_a)\right]^\dag    \left({\ensuremath{u_\mathrm{rel}}}^{a , p } - v^a_c\right)\psi\ , d^{3a}r\,,\end{aligned}\ ] ] so that @xmath157f_l(\eta_a , k_ar )   + c_a(r)g_l(\eta_a , k_ar)\right\}/r.\ ] ] many derivations of eq . ( [ eq : integral - schematic ] ) and equivalent ( through analytic continuation ) bound - state expressions may be found in the literature ( e.g. refs . @xcite ) . several observations may be made at this point . the choice to add and subtract the point - coulomb interaction @xmath158 in eq . ( [ eq : h - breakdown ] ) has two important consequences . first , it guarantees that the overlap @xmath159 computed from eq . ( [ eq : integral - schematic ] ) is a linear combination of @xmath160 and @xmath161 , as required for solutions of the schrdinger equation governing @xmath20 . second , the nuclear interaction is short ranged , so that at large separations , the interaction between daughter nuclei is dominated by the monopole term of their coulomb interaction . since this is just equal to @xmath158 , the difference @xmath162 goes rapidly to zero for @xmath120 beyond the range of the nuclear interaction . for the interactions and wave functions discussed below , this typically happens at @xmath163 fm . thus , @xmath164 , and @xmath165 takes on approximately its asymptotic value for all @xmath166 fm . these two properties make eq . ( [ eq : integral - schematic ] ) especially useful for computing asymptotic properties of @xmath20 . indeed , similar integrals tend to appear in scattering theory for exactly this reason . comparison of eqs . ( [ eq : integral - schematic ] ) and ( [ eq : standing - wave - trig ] ) shows that since @xmath167 , @xmath168 as discussed in many places in the literature ( e.g. @xcite ) . in fact , eq . ( [ eq : kmatrix - integral ] ) may be improved upon significantly to yield estimates of @xmath101 that are second - order variational in approximations to @xmath20  @xcite . ( these improvements require considerably more computation , so i do not pursue them here . ) from the discussion in sec . [ sec : plausibility - and - channel - setup ] , at resonance there is no contribution from the regular function so that @xmath155 and @xmath101 has a pole . the residue of @xmath101 at this pole ( for @xmath20 normalized to unity within some radius @xmath169 that contains the interaction region ) is @xmath170 up to corrections of order @xmath171\ , , \ ] ] discussed in refs . it is easy to show from eq . ( [ eq : k - from - s ] ) that the residue of the @xmath49-matrix pole is twice the residue of the corresponding @xmath64-matrix pole , and thus the partial width of the resonance in channel @xmath77 is given by @xmath172 again up to corrections of order @xmath173 . since the integral defining @xmath174 ( eq . ( [ eq : g - norm - integral ] ) ) is short - ranged , eq . ( [ eq : gamma - from - c ] ) for the partial width may be formulated as the square of a straightforward integral over all nucleon coordinates without specifying a boundary . this result for the width is easily generalized to include bound states . these occur at negative energy , @xmath175 , so that @xmath176 and @xmath177 . since they must be square - integrable , their overlap functions are asymptotically @xmath178 for purely imaginary @xmath84 , @xmath36 is proportional to the whittaker function , @xmath179 . working with analytic continuations of the linearly - independent pair of functions @xmath160 and @xmath36 instead of @xmath160 and @xmath161 , an integral relation analogous to eq . ( [ eq : full - integral - relation ] ) may be derived . it yields a result for @xmath180 that is guaranteed to have the correct form ( eq .  ( [ eq : bound - as - outgoing - wf ] ) ) at large @xmath120 even if it is generated from an approximate @xmath20 that does not solve the schrdinger equation exactly . in fact , an early use of integral relations of the form considered here was to produce overlaps with the correct asymptotic behavior from asymptotically - incorrect hartree - fock wave functions @xcite . they have been used more recently to compute the asymptotic normalizations @xmath181 defined by @xmath182 from shell - model and _ ab initio _ wave functions @xcite . it is evident from the discussion above that a variational approximation to the wave function @xmath20 allows two calculations of @xmath180 : one from the definition in eq . ( [ eq : direct - overlap ] ) and one from the integral relation of eq . ( [ eq : full - integral - relation ] ) . if @xmath20 , @xmath183 , and @xmath184 are all eigenstates of their respective hamiltonians , both methods yield the same result . in a typical application of many - body methods , @xmath20 is quite accurate in the interaction region but inaccurate where the daughter nuclei are widely separated . ( [ eq : direct - overlap ] ) yields @xmath180 that is only as accurate as @xmath20 is at @xmath120 and is not guaranteed to have the correct form at large @xmath120 . if @xmath185 is computed from eq . ( [ eq : full - integral - relation ] ) , its value at any radius depends only on the values of @xmath20 within the interaction region ( apart from the integration constant @xmath154 , which is nonzero only for nonresonant open channels ) . at large radius the correct asymptotic shape is guaranteed . it seems likely that since eq . ( [ eq : full - integral - relation ] ) incorporates the hamiltonian directly ( through @xmath186 ) and not just through the approximate @xmath20 , it is also more accurate at smaller @xmath120 than the directly - computed @xmath180 @xcite . in the following calculations i use variational wave functions that were computed by the vmc method . they are approximate solutions for a hamiltonian consisting of the sum of the argonne @xmath187 ( av18 ) two - body @xcite and urbana ix ( uix ) three - body @xcite interactions , which appear both in the computation of the wave functions and in the operator @xmath186 used in the integral relation . in this section i describe enough of the structure of the vmc @xmath104-shell wave functions to discuss their asymptotic properties , i explain the implementation of the integral relation , and i finally comment on the asymptotics of the wave functions . the wave functions are described in detail in ref . @xcite . the vmc wave functions begin with the spin - isospin function @xmath188 , in which the spins and isospins of the first four particles are organized into a slater determinant like those of a filled , @xmath189-particle - like @xmath190 shell . ( the zeroes denote the total angular momentum @xmath25 , its projection @xmath191 , the total isospin @xmath192 , and its projection @xmath193 ; this notation follows conventions in previous qmc work and should not be confused with the channel function @xmath194 ) . the addition of further particles while retaining antisymmetry requires spatial dependence in the wave function , i.e. , placing particles into the @xmath104-shell . this is done using orbitals @xmath195}(r_{\alpha i})$ ] defined in terms of the separation @xmath196 of particle @xmath110 from the center of mass of the `` alpha core . '' the orbitals are derived from woods - saxon potential wells and coupled to form states of definite angular momentum , parity , spatial symmetry , and isospin in an @xmath197-coupled basis , and they fall off exponentially at large @xmath196 . the effects of pairwise interactions between nucleons ( through the av18 potential ) are then accounted for using two - body scalar correlations @xmath198 that mainly account for particles avoidance of the repulsive core of the potential and are functions of the separation @xmath199 of particles @xmath110 and @xmath200 . it has been found advantageous to allow different scalar correlations depending on whether particles @xmath110 and @xmath200 are both in the @xmath201-shell ( @xmath202 ) , both in the @xmath104-shell ( @xmath203 ) , or one in each ( @xmath204 ) . finally , there are analogous three - body scalar correlations ( @xmath205 , @xmath206 ) due to both two- and three - body interactions . all of these elements ( @xmath189 core , @xmath197 coupled orbitals , and scalar correlations ) are antisymmetrized among particle permutations to make up the jastrow wave function : @xmath207 } \big ( \beta_{ls[n ] }        \phi_a(ls[n]jmtt_{z})_{p } \big ) \right\ } \ .\end{aligned}\ ] ] the index @xmath208 denotes a specific permutation of the particles into @xmath201- and @xmath104-shells ( subsequently antisymmetrized by the operator @xmath209 ) and the amplitudes @xmath210}$ ] specify the admixtures of @xmath104-shell @xmath197 states of young diagram @xmath211 $ ] . in turn , the configuration of good @xmath212 , j , m , t,$ ] and @xmath193 is given by @xmath213tt_{z})_p & = &       \phi_{\alpha}(0 0 0 0)_{1234 }       \prod_{5\leq i\leq a }       \phi^{ls[n]}_{p}(r_{\alpha i } )       \nonumber \\   & &   \times \left [ \left[\prod_{5\leq j\leq a}y_{1m_j}(\hat{r}_{\alpha j } )      \right]_{lm_l } \otimes        \left[\prod_{5\leq k \leq a}\chi_k ( \frac{1}{2}m_k)\right]_{sm_s }       \right]_{jm } \nonumber \\   & &   \times \left[\prod_{5\leq l\leq a}\nu_{l}(\frac{1}{2}t_z)\right]_{tt_z } \ , \end{aligned}\ ] ] where the spinors @xmath214 and @xmath215 specify the angular momentum and isospin states of particle @xmath110 . because the hamiltonian includes operator as well as scalar terms ( i.e. , it acts on particle spins and isospins ) , there are operator correlations in addition to central correlations . these are accounted for by writing the variational trial function as @xmath216 where @xmath217 and @xmath218 are two- and three - body operators , and @xmath219 is a symmetrization operator , needed to preserve the antisymmetry of @xmath220 because the @xmath217 and @xmath218 do not commute amongst themselves . the operator correlations are of the form @xmath221 u_q(r_{ij } ) o^q_{ij } \ , \ ] ] where @xmath222 is an operator - independent three - body correlation , and the operators @xmath223 , @xmath224 , @xmath225 , @xmath226 , and @xmath227 ( where @xmath228 are nucleon spin operators , @xmath229 , are isospin operators , and @xmath226 is the tensor operator ) , appear in the largest operator terms in the av18 potential . the @xmath230 , together with the scalar correlation @xmath231 , solve a set of coupled euler - lagrange equations with coefficients that serve as variational parameters , discussed in ref . the central correlation @xmath202 falls off exponentially to reflect the strong binding of the @xmath201-shell particles , as do the @xmath232 functions . the @xmath233 and @xmath234 correlations are constructed to approximate @xmath231 at small @xmath199 but to approach constant values at large @xmath199 . this guarantees that where particles approach each other closely the wave function is governed by the nucleon - nucleon interaction , but that the correlation between widely - separated particles is dominated by binding to a `` mean field '' accounted for in the @xmath195}$ ] orbitals . thus , the asymptotic region of @xmath235 is dominated by the @xmath195}$ ] , which have much longer tails than the @xmath202 . all of the functions appearing so far in this section are specified as functions of variational parameters , either explicitly or in the differential equations solved to compute correlations . the optimum values of those parameters are found by searching the parameter space to minimize the energy expectation value , @xmath236 with both numerator and denominator computed by monte carlo integrals over the particle coordinates . as a final step , the hamiltonian is diagonalized with respect to the @xmath104-shell configurations labeled by @xmath237 $ ] to find the coefficients @xmath210}$ ] . diagonalization both improves the variational energy of the ground state of given @xmath238 and @xmath192 and provides access to excited states . energies of unbound @xmath104-shell states can almost always be lowered by making their wave functions more diffuse ( closer to threshold ) . this is also often true of bound states , where variational energies ( but not gfmc energies ) can lie artificially above breakup thresholds because of shortcomings of the variational ansatz . for both bound and resonance states , this problem is addressed by constraining variational parameter searches to keep charge radii close to experimentally known ground - state charge radii . the values of the variational parameters for all states used in the present calculations were provided by r.  b.  wiringa @xcite . they are the results of calculations in bases of good isospin in which individual nucleons typically can not have definite charges . for convenience in defining the integral relations for neutron or proton decays , i carry out the calculations below in a basis of good nucleon charge ( essentially an @xmath40-scheme in the particle isospins ) so that the emitted nucleon is definitely either a neutron or a proton . for given variational parameters , this is only a change of representation and does not alter observables . the variational minimization has in all cases been carried out for the state of lowest @xmath193 in each isomultiplet . i compute widths of the isobaric analogues of these states by using isospin rotations of the minimized states rather than carrying out independent variational minimizations . calculations of explicit overlap functions @xmath180 ( eq .  ( [ eq : projection ] ) ) and integral relations ( eqs . ( [ eq : f - norm - integral ] ) and ( [ eq : g - norm - integral ] ) ) involve several common elements . while integral relations are much more time - consuming because they contain the operator @xmath186 , most of the computational tasks in organizing the calculation amount to constructing the channel vector @xmath239 given by eq . ( [ eq : channel - definition ] ) and then contracting either @xmath235 or @xmath240 against it . for nucleon emission , i take @xmath82 to be the wave function of the daughter nucleus and @xmath83 to be the spinor of the emitted nucleon . the integrals of eqs . ( [ eq : direct - overlap ] ) and ( [ eq : g - norm - integral ] ) are computed by monte carlo integration over the particle coordinates , using the same sampling algorithm that has long been used to compute the energy expectation value , eq . ( [ eq : variational - energy ] ) . sampling follows the metropolis algorithm , using the weight function @xmath241 where @xmath242 is the variational wave function for the @xmath243-body parent nucleus at particle coordinates @xmath144 . the delta function in eq . ( [ eq : direct - overlap ] ) is accounted for by sampling all particle coordinates and sorting the samples into narrow bins of specified channel radius ; this builds up the entire function @xmath180 from a single monte carlo walk . the normalization integral needed to give @xmath82 unit norm is computed in the same monte carlo walk by which the overlap or integral relation is computed .    only relative coordinates are used in the definitions of @xmath235 and @xmath140 , so the results of the calculations are all explicitly translation invariant , and for operator @xmath244 the quantity computed is @xmath245 the routines used to compute the integral relations were written as modified versions of existing spectroscopic - overlap routines @xcite . the integral - relation routines were used previously to compute bound - state ancs @xcite , and only very minor modification ( replacing regular whittaker functions with regular scattering functions ) was necessary for width calculations .    finally , the operator @xmath186 must be considered . it is just the potential - energy operator of the av18+uix hamiltonian , but with the restriction that only terms involving the @xmath246 ( emitted ) nucleon are considered so that eq . ( [ eq : urel - definition ] ) becomes : @xmath247 its action on @xmath235 is evaluated by calling the potential - energy routines with instructions to omit all terms purely internal to @xmath82 for the given permutation . it is then straightforward to contract @xmath248 with @xmath249^\dag$ ] for a given configuration @xmath144 . it is instructive to examine the integrand of eq . ( [ eq : g - norm - integral ] ) by inserting a delta function @xmath250 into the integral . this is evaluated just like the delta function in eq . ( [ eq : direct - overlap ] ) , by carrying out the monte carlo walk for the full integral and binning the monte carlo samples according to the channel radius instead of summing them all together . a typical result is shown in fig . [ fig : integrand ] , and nearly all cases that i have computed appear similar to this graph apart from an overall scaling . it is evident that the largest contribution comes from @xmath251 fm . it is also evident that the form of the vmc wave function beyond @xmath252 fm is irrelevant for the width , because the factor @xmath162 is zero there . this limits the integral to smaller radii and guarantees convergence of the integral regardless of what lies in the tails of the variational wave function .    ) in two @xmath90 channels is shown as a function of channel radius @xmath120 after integrating over all other coordinates . the case shown is the @xmath253 overlap onto @xmath254 + neutron , but the shapes of the curves are similar in nearly all calculations from vmc wave functions . discrete points show the integrand @xmath255 . the solid curve shows its cumulative integral , @xmath165 , which at large radius is equal to the wave function asymptotic normalization . the integrand goes to zero at small radius because it contains the regular function @xmath256 and at @xmath257 fm because @xmath258 there . the light dotted curve indicates the distribution of samples in the monte carlo integration ( with zero at the bottom of the graph).,width=336 ]    the motivation presented above for the integral relations assumes that @xmath235 , @xmath82 , and @xmath83 are all energy eigenstates . to the extent that the vmc wave functions approximate these states , the integral relations approximate the overlaps and widths _ of the hamiltonian for which they are approximate solutions_. two difficulties then present themselves in applying the integral relations to vmc wave functions : 1 ) comparison of the results to other calculations using the av18+uix hamiltonian is not possible because ( except for @xmath9he ) no such calculations have been done by other methods , and 2 ) comparison to experiment is complicated because experimental resonance energies are not reproduced exactly by the hamiltonian . the mismatch of resonance energies from experiment can be dealt with straightforwardly . the integral relations of eqs . ( [ eq : f - norm - integral ] ) and ( [ eq : g - norm - integral ] ) require an assumed channel energy @xmath133 , from which @xmath84 and @xmath85 are computed . formally @xmath133 should be the channel energy of the av18+uix hamiltonian ( known from gfmc ) , but the mismatch between this and the experimental @xmath133 is often @xmath259 mev . since the potential energy in the @xmath104-shell is much larger ( @xmath260 to @xmath261 mev ) , the experimental @xmath133 is close enough to the av18+uix channel energy that it can plausibly be used in the integrals . ( see sec . [ sec : psky ] . ) this choice accounts for most of the well - known strong dependence of the width on the resonance energy so that comparison with experiment is possible . it was found previously that using the experimental separation energy in integral - relation calculations of bound - state ancs produces results in generally good agreement with experiment @xcite . similarly , the results below indicate that using the experimentally - measured resonance energy as @xmath133 allows prediction of experimental widths . i present the results of using both experimental and gfmc values of @xmath133 in the integral relations . more consistent calculations will require hamiltonians that more precisely reproduce thresholds and resonance energies . such hamiltonians exist in the form of the illinois three - body potentials @xcite , but they have not yet been incorporated into the vmc code used here , and @xmath186 for these interactions requires considerably more computation than for the uix interaction . consider configurations in which particle @xmath243 ( before antisymmetrization ) is far from the first @xmath262 particles . because @xmath203 and @xmath204 approach constants and @xmath202 decays rapidly at large @xmath199 , the shape of @xmath235 in this part of the wave function is dominated by the shapes of the single - particle functions @xmath195}$ ] . this might be expected to give @xmath263 } \gamma_{ls[n ] } \phi_n^{ls[n]}(\omega_{ls[n]}r)\ ] ] for some amplitudes @xmath264}$ ] and jacobian - like factors @xmath265}$ ] that emerge from the correlations of eqs . ( [ eq : jastrow])-([eq : trial - wf ] ) when the overlap integral is computed . the factors @xmath265}$ ] account for the distinction between the channel radius @xmath120 and the distance @xmath196 of a @xmath104-shell particle from the center of mass of the @xmath201-shell core ; these only coincide when @xmath266 and otherwise differ by the mean difference between nucleon - core and nucleon - daughter distances . in general , the @xmath195}$ ] that emerge from the variational procedure do not yield the correct long - range asymptotic shapes for the overlaps @xmath185 . this is most readily seen for open channels , where solutions that should oscillate at large channel radius instead fall off with an assumed exponential dependence . in closed channels , the condition of square integrability gives zero for the analogue of @xmath154 in eqs . ( [ eq : full - integral - relation ] ) and ( [ eq : integral - schematic ] ) , so that eq . ( [ eq : whittaker - asymptotic ] ) holds for the true wave function . because no single @xmath237 $ ] term typically dominates a given @xmath180 , it is in general difficult to construct @xmath195}$ ] to satisfy eq . ( [ eq : whittaker - asymptotic ] ) for all possible channels at once . this difficulty is compounded by the problem that the energy expectation value driving the variational minimization receives very little contribution from the wave function tails , so the variational principle does not constrain the low - probability tails of the wave function very strongly . the application of the integral approach to bound state ancs in ref . @xcite avoided these difficulties and effectively matched the correct asymptotic form onto the better - computed interior of the wave function . pseudobound vmc wave functions approximate resonance states with square - integrable wave functions , in which the @xmath202 and @xmath195}$ ] functions cut off the wave function tails exponentially . this exponential cutoff can be understood at `` medium '' range ( 48 fm ) as forcing a resonance form on the wave function and beyond this range as providing a regularization to normalize the unbound state despite its formally nonzero amplitude at large radius . such a regularization is important for quantities like electromagnetic transition strengths ( e.g.  @xcite ) and for the approximate relation between the asymptotic normalization and the partial width given in eq . ( [ eq : gamma - from - c ] ) @xcite . as long as the resonance wave function is computed reasonably accurately within the region where the integral of eq . ( [ eq : g - norm - integral ] ) is nonzero , and it is normalized to unity over the region where @xmath235 is larger in amplitude than in the asymptotic region , the integral approximates the asymptotic normalization of @xmath161 and thus the partial width . the cutoff of @xmath235 at large nucleon separation is illustrated with overlap functions in fig . [ fig : fixing - asymptotics ] .     of the unbound @xmath267 state of @xmath1b onto @xmath268 proton , computed directly from eq . ( [ eq : direct - overlap ] ) ( discrete points with monte carlo error bars ) and from the integral relation of eq . ( [ eq : full - integral - relation ] ) ( solid curves ) . at large radius , the directly - computed overlap falls off exponentially because the @xmath195}$ ] in the vmc wave function fall off exponentially . however , the integral - relation overlap ( with @xmath155 ) gives the exact proportionality to @xmath161 expected at resonance . the absolute normalization @xmath269 multiplying @xmath161 gives the partial width of @xmath270 in each @xmath271 channel ; the normalized asymptotic is shown as a dashed curve for each channel and manifestly merges with the full overlap at @xmath272 fm.,width=336 ] i now apply the integral relation of eqs . ( [ eq : gamma - from - c ] ) and ( [ eq : g - norm - integral ] ) to compute widths of several unbound energy levels in nuclei of mass numbers @xmath0 , from vmc wave functions . i mostly choose energy levels that empirically have small width ( under 1 mev ) and are dominated by nucleon emission . the integral relation is valid for decay channels in which all three wave functions ( one parent and two daughters ) are known , but i limit this first examination to nucleon - emission channels , where one `` daughter nucleus '' is a neutron or proton . i concentrate on two - body final states , but i also model some three - body decays as sequential processes , e.g. the @xmath2b ground state decaying to a proton and the unbound ground state of @xmath1be . the subsequent decay of @xmath1be to two alpha particles can be neglected because of its small width . ( in principle the widths of unbound daughter states should be integrated over , but i find that in all cases considered this correction is much smaller than either the experimental errors or the widths of omitted decay channels . )    the results of the width calculations are shown in table [ tab : big - width - table ] . each calculation was carried out twice : once assuming the channel energy @xmath133 from the av18+uix hamiltonian ( known from gfmc calculations ) and once using the experimentally - known resonance energy . ( for the second @xmath273 states in @xmath1b and @xmath1li and the @xmath274 states in @xmath2be and @xmath2b , av18+uix energies have never been computed with gfmc ) . where the experimental channel energy is unknown or uncertain , i have used instead the gfmc energy for the av18+il7 hamiltonian , which gives a better overall fit to experimental energies than av18+uix . i have taken the experimental energies and widths in table [ tab : big - width - table ] mainly from data compilations @xcite . nearly all of the widths in the compilations are `` observed '' widths , which coincide with the full width at half maximum ( fwhm ) of cross section peaks and are proportional to the sum of the @xmath64-matrix pole residues . the integral - relation widths should correspond to these quantities . in two cases ( @xmath275 and @xmath276 ) , examination of the source literature indicated that the available experimental widths are `` formal '' widths of the @xmath277-matrix formalism . i have converted the experimental widths listed for those states to `` observed '' widths by the usual relation @xmath278 where @xmath78 ( taken to be 4 fm ) is a matching radius , @xmath279 is the shift function of @xmath277-matrix theory @xcite , prime denotes its derivative , and @xmath280 is defined by the formal width @xmath281 and the penetrability function @xmath282 . ( see the section `` definitions of resonance parameters '' in ref . @xcite . ) the sources of the experimental widths are indicated in table  [ tab : big - width - table ] . the widths presented in the table are sums over all @xmath104- and @xmath283-wave decay channels . the predicted @xmath283-wave contributions are less than 1% of the total in all cases except the decay of @xmath284 to @xmath285 , where it is computed to be 23% of the total . the @xmath104-wave decays are in most cases sums of @xmath286 and @xmath287 contributions . the decompositions into @xmath286 and @xmath287 ( or into a channel - spin coupling scheme ) are available on request to the author . the table only includes cases for which established experimental or gfmc energies are available . i have also computed widths of several states of @xmath2he with varying assumptions about resonance energies as described below , but i have not included these numbers in table [ tab : big - width - table ] . .the results of integral - relation calculations of widths . results are shown from calculations in which the channel energy was assumed equal to its experimental value ( `` from expt energy '' ) and to the av18+uix channel energy known from gfmc ( `` from av18+uix energy '' ) . error estimates are quadrature sums of monte carlo sampling errors and uncertainties from the input energies . where no experimental energy is available , the results in the `` experiment '' columns were computed using the gfmc energy with the av18+il7 hamiltonian , and they are indicated by square brackets . the column `` matches @xmath288 ? '' indicates whether the overlap function appears consistent with a resonance state , as discussed in sec . [ sec : psky ] . see sec . [ sec : psky ] and eq . ( [ eq : zeta ] ) for the definition of @xmath289 . energies are relative to the decay threshold in the center - of - mass frame . experimental energies are taken from data compilations @xcite unless otherwise noted . [ cols="<,<,<,<,<,<,<,^ , < " , ]     the computed widths of those states and channels for which experimental data are available are shown graphically in fig . [ fig : width - width ] . only for those states where partial widths are available from experiment do i show partial widths in the graphs  otherwise total widths are shown . [ fig : width - width ] demonstrates the wide dynamic range of the integral method , extending from 0.0005 to roughly 1 mev . with the exception only of the very broad @xmath290 states , which present problems discussed below , all computed widths are within a factor of three of experiment  within a factor of two if states with uncomputed alpha and direct three - body width are omitted . the error - weighted mean ratio @xmath291 of the integral - method width to the experimental width is @xmath292 , and the @xmath293 statistic for the difference between computed and experimental widths is 5.9 per degree of freedom ( omitting @xmath9he for a total of 19 states ) . restricting consideration to those states with no omitted channels gives @xmath294 and @xmath295 with @xmath296 . the errors used to compute these statistics are mainly in the measured width @xmath297 .    setting aside the resonance energies , the main sources of uncertainty in the theoretical calculation are in the accuracy of the variational wave function and of the potential . because it is unclear how to estimate these errors , the theoretical errors reported in table  [ tab : big - width - table ] are the quadrature sums of monte carlo sampling errors with errors propagated from the input resonance energies . errors propagated from experimental energies are typically small compared with the other errors . however , each gfmc channel energy is the difference of a resonance energy and a threshold energy from separate calculations . the error on this difference is often large compared with the channel energy , and this propagates to a large error on the predicted width .     for those cases included in the averages . the vertical axis shows the ratio @xmath298 of the integral - relation width to the measured width , and the horizontal axis shows the measured width @xmath297 . the horizontal dotted line shows @xmath299 as a guide to the eye . the errors shown are dominated by those on the experimental widths . experimental errors have been added in quadrature to smaller monte carlo statistical errors and errors propagated from the energy used in computing @xmath300 . states for which important alpha or 3-body decays have been neglected are indicated with asterisks next to the state labels , and large ( red ) circles indicate overlaps that seem inconsistent with @xmath288 phase shift . the outlying @xmath301 states of @xmath302be and @xmath302li are complicated by the need to include the small alpha decay channel in multichannel @xmath277-matrix fits when extracting widths from data ; this introduces a considerable spread in reported nucleon widths @xcite.,width=336 ]      the results presented in table [ tab : big - width - table ] and fig . [ fig : width - width ] vary in the degree to which they match the experimental widths . some , like the first @xmath303 state of @xmath1li , are close matches . others , like the two low - lying states of @xmath9he , are very far from agreement with experiment . the @xmath9he cases are particularly interesting , because they are the only ones for which the correct widths for the av18+uix hamiltonian are known from explicit scattering calculations @xcite . the known av18+uix widths are 1.5 mev for the @xmath304 state and 5.0 mev for the @xmath305 state , while applying the integral relation to the pseudobound vmc wave functions gives 0.68 and 0.71 mev , respectively . the @xmath305 state lies above the centrifugal barrier and is so broad that its phase shift does not pass through @xmath288 , so it is unsurprising that the integral method fails for this case . the @xmath304 state does not present these difficulties , so its difference from both experiment and theory is clearly not a shortcoming of the potential but rather of the computational methods . there are several reasons that the application of the integral relation to a pseudobound variational wave function could fail to yield the correct width for the potential . i begin by noting that for the states examined , the function @xmath306 typically does not deviate far from its leading - order dependence on @xmath307 over the range of @xmath86 where @xmath308 differs significantly from zero . the main effect of changing the assumed resonance energy is therefore to change the overall scale of the integrand in eq . ( [ eq : g - norm - integral ] ) without changing its shape much . thus , simply using the experimental channel energy in the integral relation should correct rather accurately for the mismatch between experimental and theoretical channel energies without introducing significant distortions in the integrand . it is possible that the variational minimization with constrained charge radius fails to produce good approximations to some resonance states . it could happen that the variational ansatz is a poor match to a particular resonant wave function or that the vmc wave function , being only approximate , contains `` contamination '' from nonresonant continuum states . contamination more energetic than the desired state is precisely what the gfmc method is intended to remove , and much of the success of the qmc methods lies in the exclusion of low - energy excitations from the variational wave function . any off - resonance contamination produces contributions to the integral relation that do not correspond to the pole residue . the danger of contamination would seem to grow with the resonance width , because a less - peaked density of states implies a wider range of states that are similar to the resonance and difficult to eliminate by variational minimization . the presence of nearby resonances in the same channel may compound this problem by contributing contamination with rather different structure from the resonant wave function being sought . yet another kind of difficulty lies in the normalization of the wave function . the relationship between the width and the asymptotic normalization of standing - wave states found in refs . @xcite depends on the wave function being normalized in a finite volume and its amplitude being small at the boundary of that finite volume as expressed in eq . ( [ eq : humblet - correction ] ) . i have assumed that the normalization volume is effectively defined by the exponential fall - off of the @xmath202 and @xmath195}$ ] correlations in the tails of the variational wave function . if either the true or the computed wave function fails to fall off rapidly enough , the normalization is problematic in ways that can be viewed either as the lack of an effective cutoff or as neglect of a large surface - amplitude correction . the normalization problem is perhaps the most straightforward to examine . the final column of table [ tab : big - width - table ] displays a parameter @xmath309 where @xmath310 is the location of the maximum @xmath311 , and @xmath180 is computed directly from eq . ( [ eq : direct - overlap ] ) . this ratio measures the amplitude of the variational wave function just outside the interaction region relative to that in the interior ( accounting for the @xmath312 dependence of the volume of a spherical shell ) . we may expect to encounter difficulties when @xmath313 , and this occurs for four states : the two @xmath9he states and the @xmath305 and @xmath314 states of @xmath302he . the relevance of @xmath289 may also be viewed in light of eq . ( [ eq : humblet - correction ] ) . if i adopt 8 fm as the effective boundary of the normalization volume and assume ( cf . [ fig : psky ] below ) that typically @xmath315 , then eq .  ( [ eq : humblet - correction ] ) gives @xmath316 for all but a very few of the widths computed here . much larger values of @xmath173 occur for @xmath317 , @xmath318 , and minor decay channels of a few other states . the overlap functions may also be tested for consistency with expectations for a resonance . the integral relation for the width is the @xmath319 limit of @xmath185 as computed from eq . ( [ eq : full - integral - relation ] ) . given a wave function , it is possible to compute the overlap function at all radii using eq . ( [ eq : full - integral - relation ] ) given @xmath154 , and the result is likely more accurate than that from direct calculation of @xmath185 . determining @xmath154 from a variational wave function is a tractable problem @xcite , but i do not pursue it here . since i have tacitly _ assumed _ in computing widths that @xmath155 as required for a resonance state ( cf . ( [ eq : kmatrix - integral ] ) ) , i can compute the full overlap function from the integral relation ( eq . ( [ eq : full - integral - relation ] ) ) with @xmath155 and check that it matches the overlap function computed directly from ( eq . ( [ eq : direct - overlap ] ) ) . if the two overlaps are in rough agreement in the interaction region , then the vmc wave function is consistent with a resonance wave function and may validly be used to compute a width .    as a test of the approach , i apply it to bound states , since their integration constant corresponding to @xmath154 is zero by definition . results are shown for several states in fig . [ fig : bound - overlaps ] , and i have computed them for all of the channels considered in ref . the agreement between the two calculations , especially for @xmath201-shell nuclei where the vmc method is more accurate , is excellent . the deviations of the integral - relation overlap from the direct overlap are likely to be improvements : the integral relation contains more information about the potential than does the vmc wave function by itself , and it guarantees the correct @xmath319 asymptotics . for some nuclei with @xmath320 , 4 , and 7 , gfmc calculations of overlaps exist ( albeit for the av18+il7 hamiltonian , not av18+uix ) @xcite . the results of the gfmc calculations ( dashed curves of fig . [ fig : bound - overlaps ] ) are not severely different from those of applying the integral method to vmc wave functions , and they deviate from the vmc overlaps by similar amounts . this experience with bound states indicates that overlaps computed from integral relations are at least as accurate as directly - computed overlaps of vmc wave functions and are not in conflict with gfmc results . they may therefore be very useful for calculations of spectroscopic factors and transfer and knockout cross sections .    ) ( solid curves ) and from the definition in eq . ( [ eq : direct - overlap ] ) ( squares with monte carlo error bars . ) the @xmath87 and @xmath90 quantum numbers of the `` virtually emitted '' nucleon are indicated by labels near the appropriate curves ( and distinguished by color ) . where they exist , overlaps computed directly from gfmc wave functions  @xcite are shown as ( black ) dashed curves.,width=528 ]    overlap functions are shown for @xmath270 in fig . [ fig : fixing - asymptotics ] and for several other resonance states in fig . [ fig : psky ] . these were computed both directly and by integral relations with @xmath155 . in all cases , the integral relation replaces the artificial exponential fall - off of the vmc wave function with the oscillatory behavior of @xmath321 . based on experience with bound states , the two overlaps should agree at @xmath322 fm . in some cases ( e.g.  @xmath323 ) the match there is quite good , while in others ( e.g.  @xmath324 ) overlaps from the two methods seem to have little to do with each other . i conclude that when results of the two methods are qualitatively very different inside 4 fm , there is an inconsistency with the assumption of @xmath288 phase shift so that the pseudobound vmc wave function may not allow accurate width calculations . the penultimate column of table [ tab : big - width - table ] indicates for each width a qualitative judgment of whether the two methods of computing overlaps agree , and cases where they do not are indicated in fig . [ fig : width - width ] by large ( red ) circles . examination of the table reveals that with the exception of the @xmath325 state in @xmath1b , failures of the `` @xmath326 assumption '' occur only in a few neutron emission channels and in no proton emission channels . considering only states consistent with @xmath327 , i find @xmath328 and @xmath329 , this time for @xmath330 instead of 18 . further restriction to states in which all decay channels are computed gives @xmath331 and @xmath332 for 8 degrees of freedom . these are essentially the same results as for the entire data set . it appears that the best predictor of whether a calculated width will match experiment is simply whether the calculation includes all channels contributing to the measured width .    ) ; points with monte carlo errors ) and from the integral relation of eq . ( [ eq : full - integral - relation ] ) with @xmath155 ( or @xmath333 ; solid curves ) . where more than one angular momentum channel is available , the channels are labeled . the four cases on the left show good agreement between the two methods at small separation ( @xmath322 fm ) , and the corresponding widths agree with experiment . the two states on the right are inconsistent with the @xmath334 assumption since the two types of overlap calculations disagree at @xmath335 fm . computed widths of these two states also disagree with experiment.,width=528 ]      a difficulty presents itself in considering the pair of @xmath4 states of @xmath1be at @xmath337 mev and the pair of @xmath5 states at @xmath338 mev . each of these doublets consists of one @xmath339 and one @xmath340 state that have mixed , with the lower state of each pair predominantly @xmath340 . the vmc wave functions have definite isospin , so incorporating the mixing into calculations of their widths presents ambiguities . simply ignoring the mixing and assigning the @xmath340 widths to the lower states and @xmath339 widths to the higher states does not work . for example , in the @xmath303 doublet , the lower state has measured width @xmath341 kev and the upper has @xmath342 kev ; the unmixed widths from the integral method are 177 kev and 364 kev , respectively . there exist preliminary gfmc calculations of the mixing matrix element of the hamiltonian for each of these doublets @xcite , and this matrix element can be combined with the splitting of the doublet to determine its mixing angle ( e.g. @xcite ) . the splittings of the doublets are small enough that they are not resolved by the existing gfmc calculations , so the experimental splittings must be used to compute mixing angles . since the recommended @xcite 165 kev splitting of the @xmath303 doublet is less than twice the gfmc mixing energy , these two numbers can not be combined to yield a real mixing angle . ( adopting a splitting of 310 kev @xcite produces good agreement with the experimental widths . ) in the @xmath5 doublet , combining the gfmc mixing matrix element with the experimental splitting produces a poor match to the data ( widths of 14.2 and 105 kev _ versus _ @xmath343 and @xmath344 kev from experiment ) .    for each doublet , a mixing angle can be computed by @xmath293 minimization of the difference between theoretical and experimental widths . this produces good agreement with the measured widths , which was not guaranteed . the resulting mixing angle for the @xmath303 doublet has @xmath345 , which is small relative to literature values @xcite . the mixing angle from minimizing @xmath293 for the @xmath5 doublet has @xmath346 , much smaller than both the value of @xmath347 from the gfmc mixing energy and the literature value of @xmath348 @xcite . thus , attempts to compute separate widths for the upper and lower states in each doublet run aground on the problem of finding mixing angles consistent with all available information . for each doublet , the quantity most independent of the mixing angle is the sum of the two widths . the sum is less sensitive to @xmath349 than are the underlying pole residues , but it is not quite independent of @xmath349 because the doublet states are not degenerate . it is the sum for each doublet that is shown in table [ tab : big - width - table ] and fig . [ fig : width - width ] , with the @xmath340 width computed in each case using the lower energy and the @xmath339 width the higher . these sums are in good agreement with experiment . an important use of theoretical estimates of widths is in the identification of observed states . most of the states considered here ( @xmath350 , @xmath18 mev , dominated by nucleon decays ) were found experimentally , and their spins and parities identified , long ago . some exceptions are the first @xmath325 and second @xmath351 states of @xmath1b and the entire low - lying spectra of @xmath2he and @xmath2li . here i attempt to shed light on these systems by calculating widths from vmc wave functions . several theoretical models predict low - lying states of @xmath1b that have not been observed , as discussed in ref . it is possible that the states are simply too broad to be seen easily in experiments , and evidence was recently found for a @xmath325 state at 1.9 mev above the @xmath1b ground state @xcite . an @xmath277-matrix fit to both elastic and inelastic scattering of @xmath302be on protons indicated a @xmath325 state with partial width @xmath353 mev for decay to the @xmath302be ground state and @xmath354 mev to the first excited state . since these widths are within @xmath355 of zero , i have not shown them in table [ tab : big - width - table ] . this state has some support from the calculation of ref . @xcite , where it was found in computed phase shifts at the same energy using the merged no - core shell model and resonating group method . that calculation also indicates reason for caution in applying the integral relation : the phase shift does not approach @xmath288 , and indeed the vmc @xmath325 state fails the @xmath288 consistency test discussed above . i show in table [ tab : big - width - table ] the predicted partial widths of the @xmath325 state at its energies for the av18+uix and av18+il7 hamiltonians . [ fig : b8_0 ] shows the dependence of this prediction on the assumed resonance energy . the claimed experimental widths are consistent with my results , though inconsistency of the overlaps with @xmath326 phase shift makes the significance of this consistency doubtful . i also include in table [ tab : big - width - table ] the av18+il7 energies and vmc - computed partial widths of the unobserved isobaric - analogue state in @xmath1li . the @xmath1li and @xmath302li vmc wave functions are explicitly isospin rotations of the @xmath1b and @xmath302be wave functions , and the @xmath356 decays also fail the @xmath288 condition on the overlaps . the authors of ref . @xcite also find a @xmath351 state at 2.55 mev excitation . because this result has not been confirmed , i omit this state from goodness - of - fit statistics , but i have computed its partial widths to the @xmath302be ground and first - excited states and included the results in table [ tab : big - width - table ] . the overlap functions for these states are compatible with @xmath288 phase shift . however , the computed partial widths do not match those claimed in ref . @xcite : i find @xmath357 mev to the @xmath302be ground state _ versus _ @xmath358 mev measured , and @xmath359 mev to the @xmath302be excited state _ versus _ @xmath360 mev measured . the origin of these differences is not clear . be ground state threshold ) for the first @xmath325 state in @xmath1b . the thicknesses of the bands indicate the errors from monte carlo sampling . the band bounded by solid curves indicates the partial width to the ground state of @xmath302be , and that bounded by dashed curves shows the partial width to the first excited state . the boxes ( offset from each other slightly in the horizontal direction for visibility ) show the experimental results of ref . @xcite without a correction from @xmath277-matrix formal width to observed width . ( this can amount to a @xmath361 reduction at the high end of the allowed width . ) also indicated is the best estimate of the resonance energy from gfmc calculations using the av18+il7 hamiltonian.,width=288 ]      the spectroscopy of @xmath2he remains unclear despite several experimental studies @xcite . the ground state was originally thought to be a @xmath305 resonance state @xcite . subsequently , strong @xmath201-wave @xmath363-@xmath1he interaction was seen near threshold and argued to reflect a @xmath364 virtual state  @xcite . more recent experiments have found a smaller scattering length and thus less support for a virtual state  @xcite . because @xmath201-wave neutrons do not have true resonances passing through @xmath288 phase shift , i do not present a width for this state . several other observations of resonances within a few mev of the @xmath1he+@xmath363 `` threshold '' have been claimed . ( see ref . @xcite for a summary . ) spin and parity assignments for all of these states are uncertain , and matching them to theoretical expectations has proven difficult @xcite . width estimates based on _ ab initio _ calculations could provide useful guidance , so i explore this possibility here . there have been four claims of a state around 1.2 mev above the @xmath1he+@xmath363 `` threshold '' @xcite , and it has additional support from a study of possible analogue states @xcite . this is generally assumed to be the lowest - lying @xmath104-shell state , with @xmath365 , and there is conflicting information concerning its width . @xcite found it to be narrower than their 0.42 mev resolution ; other experiments found 1 mev @xcite , @xmath366 mev @xcite , and 2 mev @xcite . as pointed out particularly by barker @xcite , it is difficult to reconcile widths considerably narrower than 1 mev with the expected strong single - particle character of the @xmath305 resonance . theoretical calculations also place the @xmath305 state a few mev higher ( e.g. at 3 and 4 mev in gfmc calculations with illinois-6 and illinois-2 three - body forces , respectively ) . one possibility is reduction of the @xmath305 energy by an @xmath367-shell component that is missing in the calculations  @xcite . additional resonances have been found at higher energies , one around 2.3 mev @xcite with claimed width @xmath368 mev , one around 4 mev @xcite , and another around 5 mev @xcite . in addition to the @xmath305 state , a @xmath304 state is expected theoretically , though also at higher energy . it is also likely that @xmath367-shell intruder states with @xmath369 and @xmath370 are present in the low - lying spectrum . i computed widths of @xmath104-shell states with @xmath365 and @xmath304 , but the results proved difficult to match with experimental data . the short - range overlaps of the @xmath305 state computed directly and by the integral method are in nice agreement . since it is unclear what resonance energy should be used in the calculation , i have computed widths from the vmc wave function using a range of channel energies in the integral relation . [ fig:9he ] shows these results . they mainly demonstrate the argument of barker @xcite that there is a mismatch between the narrow width of the @xmath305 state claimed by bohlen _ et al . _ @xcite and theoretical expectations of a strongly single - particle state . the resonances claimed by belozerov _ et al . _ @xcite and golovkov _ et al . _  @xcite are consistent with the computed @xmath305 width , but that assignment makes interpretation of the 1.3 and 2.4 mev states of bohlen _ et al . _ @xcite difficult : at least one would have to be an @xmath367-shell state . all channel energies below 6 mev for the @xmath304 state yield computed widths of less than 5 kev . this could match either of the states at 4.3 and 5.3 mev in ref . @xcite , which were found to be narrower than the 100-kev experimental resolution , but the vmc overlaps of this state are inconsistent with @xmath288 phase shift . for this reason , the integral - relation results are probably not reliable predictions of the width . vmc wave functions also exist for states with @xmath371 and @xmath372 , but these wave functions with @xmath367-shell components have not reached the same level of development as the vmc @xmath104-shell states and have not been published . as mentioned above , the @xmath364 state should not be observable as a resonance . the vmc wave functions for the @xmath370 and @xmath372 states indicate that they are mainly made by coupling the @xmath351 state of @xmath1he to an @xmath201-wave neutron and therefore also not true @xmath373 resonances . i did not attempt calculations of partial widths for decay to the @xmath374 state . calculations of the partial widths for decay from positive - parity states to the @xmath1he ground state by emission of a @xmath375-wave neutron yielded partial widths of less than 60 kev and overlaps inconsistent with @xmath288 phase shift . he resonances that have reported widths . the band labeled `` vmc '' shows the predicted width of the @xmath305 state as a function of the resonance energy assumed in the integral relation , and its width reflects the statistical error of the monte carlo integration . where no error was reported for a width ( for example , ref . @xcite has `` @xmath376 mev '' for the width of the state at 2.0 mev ) , i have assigned an error of 0.5 mev and indicated this with broken lines at the upper and lower limits . i have omitted some very broad states with missing or lower - limit errors : a state with @xmath377 mev and @xmath378 mev from ref . @xcite and one with @xmath379 mev and @xmath380 mev from ref . @xcite . see ref . @xcite for additional reports of states without measured widths . the displayed widths are from seth _ et al . _ @xcite , belozerov _ et al . _ @xcite , bohlen _ et al . _ @xcite , and golovkov _ et al . downward arrows indicate widths that include @xmath381 kev instrumental resolution.,width=288 ]      the most recent data compilation for @xmath383 lists five low - lying states of @xmath2li and only assigns firm spins and parities to the lowest two @xcite . a more recent paper @xcite identifies the third state as @xmath384 by comparing spectroscopic factors of vmc wave functions ( older versions of those used here ) with measured @xmath385 cross sections . those authors also assigned @xmath386 and @xmath387 , respectively , to the next two states . this assignment was based partly on the ordering of states in theoretical calculations and partly on the assumption that widths should correlate with computed spectroscopic factors . the results presented in table [ tab : big - width - table ] and fig .  [ fig : width - width ] support these assignments by reproducing the widths of all three unbound states . i conclude by considering other ways to estimate widths from vmc wave functions and comparing them with the integral method . in the absence of integral relations or explicit scattering calculations , widths must be estimated from spectroscopic factors . in applications of the shell model , one often assumes that the width is the product of the spectroscopic factor , @xmath388 ^ 2r^2\,dr\,,\ ] ] and the single - particle width . vmc spectroscopic factors might in fact be more suited to this procedure than those from a shell model , because shell models are typically confined to a single value of the principal quantum number , while the large amount of correlation in vmc wave functions guarantees that eq . ( [ eq : direct - overlap ] ) picks up contributions from all major shells . the crudest estimate of the single - particle width is the `` wigner limit '' @xcite . on the basis of a causality argument , the width of a resonance can be shown to have an approximate upper limit of @xmath389 where @xmath78 is a `` matching radius '' ( typically @xmath390 fm ) defining the edge of the interaction region , and @xmath391 is the penetration factor of ref .  @xcite . since @xmath392 is ( approximately ) an upper limit on the width that a resonance can have , it might approximate a single - particle width . some authors define @xmath392 to include an additional numerical factor multiplying eq . ( [ eq : wigner ] ) , reflecting assumptions about the wave function inside the interaction region . teichmann and wigner  @xcite assumed a constant wave function to arrive at a factor of @xmath393 . other authors make more elaborate assumptions and arrive at a factor of @xmath394 for @xmath395 @xcite . i take eq . ( [ eq : wigner ] ) to define @xmath392 .    in fig . [ fig : wigner ] i use the wigner limit to estimate widths of the states under consideration . in each case , i multiply @xmath392 for a 4 fm radius by the vmc spectroscopic factor from eqs . ( [ eq : direct - overlap ] ) and ( [ eq : spectroscopic - factor ] ) . these estimates plainly do not reproduce measured widths as well as the integral relation . the weighted mean ratio of `` theoretical '' to experimental width for this method , restricted to states consistent with @xmath288 phase shift and purely nucleon - emission decay , is @xmath396 ; the reduced @xmath293 is 1845 for eight degrees of freedom . ( recall that @xmath397 with @xmath398 for the same set of states . ) the mismatch between @xmath399 and @xmath297 can be reduced by choosing a smaller numerical factor to define the wigner limit , but that does not remove the large scatter in @xmath400 . better agreement with experiment is likewise achieved with a smaller radius , but good agreement requires an unphysically small radius in the neighborhood of 2 fm , which again does not remove the large scatter in the ratio @xmath401 . estimates from @xmath392 are typically @xmath402 smaller if they are estimated as `` observed widths '' from eq . ( [ eq : observed - width ] ) , using @xmath403 this helps significantly with neither the overall scale nor the large scatter of the predicted widths .    , and the labels from that figure may be used to identify states here . comparison with fig . [ fig : width - width ] indicates that the integral relation is a significantly better predictor of widths than @xmath399 . note the different vertical scale from fig . [ fig : width - width].,width=336 ]    a better estimate of the single - particle width , and one perhaps more widespread in shell - model studies , is based on woods - saxon potentials . one assumes a potential well of `` standard '' radius and diffuseness and adjusts its depth to produce a resonance at the correct energy . the width @xmath404 of this resonance is then taken as an estimate of the single - particle width . a range of geometric parameters for the potential is usually considered , because the most appropriate values are not known _ a priori_.    i computed estimates of this kind , using diffuseness 0.65 fm and woods - saxon radius @xmath405 , with @xmath243 the mass number of the resonance state ( the defaults in a code provided by b.  a.  brown @xcite ) . i neglected variation of these parameters and estimated the single - particle width @xmath406 from the fwhm of the peak in the computed cross section for the given woods - saxon well . the comparison of @xmath407 with @xmath297 is shown in fig . [ fig : woods - saxon ] . the woods - saxon estimates are systematically low , with a weighted mean @xmath408 for the same eight cases considered above and @xmath409 for eight degrees of freedom . i conclude from this exercise and the similar exercise using wigner limits that vmc widths computed by the integral method are more useful predictors of experimental widths than are the vmc spectroscopic factors .    , divided by experimental widths . symbols are as in fig . [ fig : width - width ] , and the labels from that figure may be used to identify states here . comparison with figs . [ fig : width - width ] and [ fig : wigner ] indicates that this approach predicts experimental widths more accurately than the wigner - limit approach but still not as well as the integral relation.,width=336 ] i have presented plausibility arguments , supported by detailed derivations in the literature , that widths of resonant states can be estimated by evaluating an integral over pseudobound _ ab initio _ wave functions . this approach is approximate , but it avoids a great deal of computation and human labor that would be needed in full - on scattering calculations and would often be complicated by coupled channels . it is nicely suited to quantum monte carlo calculations in that it is insensitive to the difficult - to - compute tails of the many - body wave functions , it involves a short - range integral amenable to monte carlo integration , it uses more information about the hamiltonian than is encoded in the variational wave function , and it can be applied to resonances narrower than the practical energy resolution of the gfmc technique . related integrals yield overlap functions for bound states , and these overlaps are guaranteed to have the correct shapes in their long - range asymptotics even when the variational wave function does not . these may be useful for calculations of spectroscopic factors and of transfer and knockout cross sections . i have implemented integral - method width calculations for one - nucleon emission from wave functions computed by the variational monte carlo method . it yields widths in good agreement with experiment for several states in the @xmath410 mass range . cases of disagreement always involve either open channels for which i have not accounted or a resonant wave function that is not strongly peaked in the interaction region . i have shown that widths predicted in this way are closer matches to experiment than are nave combinations of _ ab initio _ spectroscopic factors with wigner - limit or woods - saxon estimates of single - particle widths . the integral method is thus a useful tool for estimating widths from _ ab initio _ methods that produce pseudobound wave functions .    for the longer term , the calculations presented here represent a learning problem for application of integrals of the type in eqs . ( [ eq : f - norm - integral])-([eq : g - norm - integral ] ) to qmc wave functions . application to gfmc wave functions of the methods used here will be straightforward and mainly involve additional bookkeeping similar to that used in ref . @xcite for direct overlap calculations . integrals of the kind considered here are likely to find their most extensive use in _ ab initio _ calculations of coupled - channel scattering and reactions , and a major goal of the work presented here is to prepare the way for such calculations . i acknowledge useful discussions with c.  a.  bertulani , i.  brida , b.  a.  brown , c.  r.  brune , h.  esbensen , a.  m.  mukhamedzhanov , s.  c.  pieper , and j.  p.  schiffer . i thank r.  b.  wiringa for providing variational wave functions and much guidance in their use . this work was supported by the u.s . department of energy , office of nuclear physics , under contract no . de - ac02 - 06ch11357 . calculations were performed on the fusion computing cluster operated by the laboratory computing resource center at argonne .
Abstract:

Answer: i describe the computation of energy widths of nuclear states using an integral over the interaction region of _ ab initio _ variational monte carlo wave functions , and i present calculated widths for many states . i begin by presenting relations that connect certain short - range integrals to widths . i then present predicted widths for @xmath0 nuclei , and i compare them against measured widths . they match the data more closely and with less ambiguity than estimates based on spectroscopic factors . i consider the consequences of my results for identification of observed states in @xmath1b , @xmath2he , and @xmath2li . i also examine failures of the method and conclude that they generally involve broad states and variational wave functions that are not strongly peaked in the interaction region . after examining bound - state overlap functions computed from a similar integral relation , i conclude that overlap calculations can diagnose cases in which computed widths should not be trusted .