it is an important task to visualize , understand and then extract useful information in complex continuous dynamical systems from many science and engineering fields including flight simulation @xcite , wave propagation @xcite , bio - inspired fluid flow motions @xcite , ocean current model @xcite , and etc .
it is , therefore , natural to see many existing tools to study dynamical systems .
for example , a lot of work in dynamical systems focuses on understanding different types of behaviors in a model , such as elliptical zones , hyperbolic trajectories , chaotic attractors , mixing regions , to name just a few .
one interesting tool is to partition the space - time domain into subregions based on certain quantity measured along with the passive tracer advected according to the associated dynamical system . because of such a lagrangian property in the definition , the corresponding partition is named the lagrangian coherent structure ( lcs ) .
one commonly used quantity is the so - called finite time lyapunov exponent ( ftle ) @xcite which measures the rate of separation of a passive tracer with an infinitesimal perturbation in the initial condition , over a finite period of time .
another popular diagnostic of trajectory separation in dynamical systems is the finite - size lyapunov exponent ( fsle ) @xcite which is especially popular in the applications from the oceanography . instead of introducing an infinitesimal perturbation in the initial condition of the fluid particles
, fsle considers finite perturbation .
it does not measure the rate of separation but computes the time required to separate two adjacent tracers up to certain distance .
since all these quantities are long treated as a lagrangian property of a continuous dynamical system , most , if not all , numerical methods are developed based on the traditional lagrangian ray tracing method by solving the ode system using any well - developed numerical integrator .
these approaches , however , require the velocity field defined at arbitrary locations in the whole space depending on the location of each individual particle .
this implies that one has to in general implement some interpolation routines in the numerical code .
unfortunately , it could be numerically challenging to develop an interpolation approach which is computationally cheap , monotone , and high order accurate .
an eulerian approach to the flow map computation for the ftle has been first proposed in @xcite which has incorporated the level set method @xcite .
based on the phase flow method @xcite , we have developed a backward phase flow method for the eulerian ftle computations in @xcite . based on these partial differential equation ( pde ) based algorithms , more recently , we have developed in @xcite an efficient eulerian numerical approach for extracting invariant sets in a continuous dynamical system in the extended phase space ( the @xmath1 space ) . in @xcite , we have proposed a simple algorithm called vials which determines the average growth in the surface area of a family of level surfaces .
this eulerian tool relates closely to the ftle and provides an alternative for understanding complicated dynamical systems .    in this paper
, we are going to further develop these eulerian tools for continuous dynamical systems .
there are several main contributions of this paper .
we will first improve the numerical algorithm for _ forward _ flow map construction .
the papers @xcite have proposed to construct the _ forward _ flow map defined on a fixed cartesian mesh by solving the level set equations or the liouville equations in the _ backward _ direction . in particular , to compute the forward flow map from the initial time @xmath2 to the final time @xmath3 , one needs to solve the liouville equations _
backward _ in time from @xmath3 to @xmath2 with the initial condition given at @xmath3 .
this implementation is numerically inconvenient , especially when incorporating with some computational fluid dynamic ( cfd ) solvers , since the velocity field is loaded from the current time @xmath3 backward in time to the initial time .
this implies that the whole field at all time steps has to be stored in the desk which might not be practical at all . in this paper
, we first propose an eulerian approach to compute the _ forward _ flow map _ on the fly _ so that the pde is solved _ forward _ in time .
developed using this improved algorithm , we propose a simple eulerian approach to compute the isle .
for each of the separation factor @xmath0 in the definition of the isle , typical lagrangian method requires to shoot an individual set of rays and to monitor the variations in these trajectories .
if the scale @xmath0 is changed , these methods will require to restart the whole computations all over again .
in this work , we are going to develop an efficient algorithm which provides the isle field for any arbitrary separation factor @xmath0 .
for an individual value of @xmath0 , the method requires to extract only an isosurface of a volumetric data .
this can be done easily by any well - developed interpolation algorithm or simply the function isosurface in matlab .
another main contribution of this work is to provide a theoretical link between the ftle and the isle .
even though these two quantities are measuring different properties of the flow , it has been reported widely , such as @xcite , that they give visually very similar solutions in many examples .
we are going to show in particular that their ridges can be identified through the ridge of the largest eigenvalue of the associated cauchy - green deformation tensor of the flow . to the best of our knowledge
, this theoretical result is the first one quantitatively revealing the relationship between the ftle ridges and the isle ridges .
indeed , comparisons between ftle and fsle have been made in various studies .
the paper @xcite has pointed out that the ftle might be unable to recognize the boundaries between the chaotic and the large - scale mixing regime .
the article @xcite has argued that the fsle has several limitations in lagrangian coherence detection including aspects from local ill - posedness , spurious ridges and intrinsic jump - discontinuities .
nevertheless , the work @xcite has stated that these two concepts might yield similar results , if properly calibrated , and could be even interchangeable in many flow visualizatons such as oceanography and chemical kinetics . in this paper , however , we would like to emphasize on the computational aspect of using these lyapunov exponents , rather than on comparing the advantages and the disadvantages of these different methods .
the preference for one or the other could be based on the tradition in various fields , or the availability of the computational resources .
this paper is organized as follows . in section [ sec : background ] we will give a summary of several important concepts including the ftle , fsle and isle , and also our original eulerian formulations for computing the flow maps and the ftle . then , our proposed eulerian algorithms will be given in section [ sec : proposal ] .
we will then point out some theoretical properties and relationships between the ftle and the isle in section [ sec : relation ] .
finally , some numerical examples will be given in section [ sec : examples ] .
in this section , we will summarize several useful concepts and methods , which will be useful for the developments that we are proposing .
we first introduce the definition of various closely related lyapunov exponents including the finite time lyapunov exponent ( ftle ) @xcite , the finite size lyapunov exponent ( fsle ) @xcite and also the infinitesimal size lyapunov exponent ( isle ) @xcite .
then we discuss typical lagrangian approaches and summarize the original eulerian approaches as in @xcite .
the discussions here , however , will definitely not be a complete survey .
we refer any interested readers to the above references and thereafter .      in this paper , we consider a continuous dynamical system governed by the following ordinary differential equation ( ode ) @xmath4 with the initial condition @xmath5 .
the velocity field @xmath6 is a time dependent lipschitz function where @xmath7 is a bounded domain in the @xmath8-dimensional space . to simplify the notation in the later sections
, we collect the solutions to this ode for all initial conditions in @xmath9 at all time @xmath10 and introduce the flow map @xmath11 such that @xmath12 represents the arrival location @xmath13 of the particle trajectory satisfying the ode ( [ eqn : ode ] ) with the initial condition @xmath14 at the initial time @xmath15 .
this implies that the mapping will take a point from @xmath16 at @xmath15 to another point @xmath13 at @xmath17 .    the finite time lyapunov exponent ( ftle ) @xcite measures the rate of separation between adjacent particles over a finite time interval with an infinitesimal perturbation in the initial location .
mathematically , consider the initial time to be 0 and the final time to be @xmath18 , we have the change in the initial infinitesimal perturbation given by @xmath19 the leading order term of the magnitude of this perturbation is given by @xmath20^ * \nabla \phi_0^t(\mathbf{x } ) \delta
\mathbf{x}(0 ) \right > } \ , .\ ] ] with the cauchy - green strain tensor @xmath21^ * \nabla \phi_0^t(\mathbf{x})$ ] , we obtain the maximum deformation @xmath22 } } \| \mathbf{e}(0)\| = e^ { \sigma_0^t({\mathbf{x } } ) |t| } \|\mathbf{e}(0 ) \| \ , , \ ] ] where @xmath23 is the eigenvector associated with the largest eigenvalue of the deformation tensor . using this quantity , the finite - time lyapunov exponent ( ftle ) @xmath24
is defined as @xmath25 } = \frac{1}{|t| } \ln \sqrt{\lambda_0^t({\mathbf{x } } ) } \ , .
\label{eqn : ftle}\ ] ] where @xmath26 denotes the largest eigenvalue of the cauchy - green tensor .
the absolute value of @xmath18 in the expression reflects the fact that we can trace the particles either _ forward _ or _ backward _ in time . in the case when @xmath27 , we are measuring the maximum stretch _
backward _ in time and this corresponds to the maximum compression _ forward _ in time . to distinguish different measures , we call @xmath24 the _ forward _ ftle if @xmath28 and the _ backward _ ftle if @xmath27 .
a related concept is a widely used concept in oceanography called the finite size lyapunov exponent ( fsle ) which tries to measure the time it takes to separate two adjacent particles up to certain distance @xcite .
there are two length scales in the original definition of the fsle .
one is the initial distance between two particles which is denoted by @xmath29 .
the other one is the so - called separation factor , denoted by @xmath30 . for any trajectory with initial location @xmath31 with @xmath32 ,
we first determine the shortest time @xmath33 so that @xmath34 the fsle at the point @xmath35 is then defined as @xmath36    numerically , on the other hand , the maximization over the constraint @xmath32 might require special considerations and the corresponding optimization problem might not be easily solved by typical optimization algorithms .
one simple approximation is to consider only @xmath37-neighbors of @xmath35 with each point @xmath38 obtained by perturbing one coordinate of @xmath35 with distance @xmath29 .    to avoid introducing two separate length scales ( @xmath29 and @xmath0 ) and the approximation by @xmath37 neighbors in the final optimization step , one theoretically takes the limit as the length scale @xmath29 tends to 0 .
the fsle can then be reduced to the infinitesimal size lyapunov exponent ( isle ) as defined in @xcite given by @xmath39 where @xmath40 is the shortest time for which @xmath41 .
even with tremendous theoretical development , there are not many discussions on the numerical implementations . because the quantity is defined using lagrangian particle trajectories , most numerical algorithms for computing the ftle or the fsle are developed based on the ray tracing method which tries to solve the flow system ( [ eqn : ode ] ) using a numerical integrator . in the two dimensional case , for example
, the flow map is first determined by solving the ode system with the initial condition imposed on a uniform mesh @xmath42
. then one can apply simple finite difference to determine the deformation tensor and , therefore , its eigenvalues .
the computation of the ftle is a little easier . for a given final time ,
one simply solves the ode system up to that particular time level .
then the ftle at each grid location can be computed directly from these eigenvalues .
the isle , on the other hand , is a little more involved since for a particular separation factor , @xmath0 , one has to continuously monitor the change in these eigenvalues and record the precise moment when it reaches the thresholding value .
computational complexity could be one major concern of these numerical tools .
let @xmath43 and @xmath44 be the number of grid points in each physical direction and in the temporal direction , respectively .
the total number of operations required to compute the flow map @xmath45 is , therefore , @xmath46 .
note however that such complexity is needed for determining the ftle or the isle on one single time level . for all time levels ,
the computational complexity for these usual lagrangian implementation is @xmath47 .
there are several nice ideas to improve the computational time .
one natural idea is to concentrate the computations near the locations where the solutions have large derivatives by implementing an adaptive mesh strategy @xcite .
another interesting approach is to avoid the computation of the overall map @xmath48 but to numerically decompose it into sub - maps @xmath49 @xcite .
then , the extra effort from @xmath50 to @xmath51 , given by @xmath52 are simply some interpolations .
all these proposed improvements are only on the flow map computations , and therefore the ftle . for the fsle or
the isle computations on one single time level and one particular @xmath0 , say @xmath53 , one has to go through all @xmath54 mesh points in the temporal direction . for a different time level or even a different separation factor @xmath0
, we have to go though the same procedure all over again .
we do not aware of many efficient algorithms for the fsle or the isle implementation .
( a ) @xcite .
( a ) lagrangian ray tracing from a given grid location @xmath55 at @xmath2 . note that @xmath56 might be a non - grid point .
( b ) eulerian values of @xmath57 at a given grid location @xmath56 at @xmath3 gives the corresponding take - off location at @xmath2 .
note the take - off location might not be a mesh point.,title="fig:",width=217 ] ( b ) @xcite .
( a ) lagrangian ray tracing from a given grid location @xmath55 at @xmath2 . note that @xmath56 might be a non - grid point .
( b ) eulerian values of @xmath57 at a given grid location @xmath56 at @xmath3 gives the corresponding take - off location at @xmath2 .
note the take - off location might not be a mesh point.,title="fig:",width=217 ]    we briefly summarize the eulerian approach based on the level set method and the louville equation .
we refer interested readers to @xcite and thereafter .
we define a vector - valued function @xmath58 . at @xmath2
, we initialize these functions by @xmath59 these functions provide a labeling for any particle in the phase space at @xmath2 . in particular , any particle initially located at @xmath60 in the extended phase space can be * implicitly * represented by the intersection of @xmath8 codimension-1 surfaces represented by @xmath61 in @xmath62 . following the particle trajectory with @xmath63 as the initial condition in a given velocity field , any particle identity should be preserved in the lagrangian framework and this implies that the material derivative of these level set functions is zero , i.e. @xmath64 this implies the following level set equations , or the liouville equations , @xmath65 with the initial condition ( [ eqn : levelsetequationic ] ) .
the above * implicit * representation embeds all path lines in the extended phase space .
for instance , the trajectory of a particle initially located at @xmath66 can be found by determining the intersection of @xmath8 codimension-1 surfaces represented by @xmath67 in the extended phase space .
furthermore , the forward flow map at a grid location @xmath63 from @xmath2 to @xmath3 is given by @xmath68 where @xmath56 satisfies @xmath69 .
note that , in general , @xmath56 is a non - mesh location .
the typical two dimensional scenario is illustrated in figure [ fig : forwardbackward ] ( a ) .
the solution to ( [ eqn : levelsetequation ] ) contains much more information than what was referred to above .
consider a given mesh location @xmath56 in the phase space at the time @xmath3 , as shown in figure [ fig : forwardbackward ] ( b ) , i.e. @xmath70 in the extended phase space . as discussed in our previous work
, these level set functions @xmath71 defined on a uniform cartesian mesh in fact give the backward flow map from @xmath3 to @xmath2 , i.e. @xmath72 .
moreover , the solution to the level set equations ( [ eqn : levelsetequation ] ) for @xmath73 provides also backward flow maps for all intermediate times , i.e. @xmath74 .    to compute the forward flow map , on the other hand , @xcite has proposed to simply reverse the above process by initializing the level set functions at @xmath3 by @xmath75 and solving the corresponding level set equations ( [ eqn : levelsetequation ] ) backward in time . a typical algorithm of this type is given in algorithm 1 .
note that in the original lagrangian formulation , the attracting lcs or the unstable manifold is obtained by backward time tracing , while the repelling lcs or the stable manifold is computed by forward time integration . in this current eulerian formulation , on the other hand , forward time integration of the liouville equations gives the attracting lcs and backward time marching provides the repelling lcs .
( 4,0 ) ( -0.5,0)(1,0)12     +    algorithm 1 : computing the forward flow map @xmath76 :    1 .
discretize the computational domain to get @xmath77 .
2 .   initialize the level set functions on the last time level @xmath78 by @xmath79 3 .
solve the liouville equations for each individual level set function @xmath80 @xmath81 from @xmath82 down to @xmath2 using any well - developed high order numerical methods like weno5-tvdrk2 @xcite with the boundary conditions @xmath83 where @xmath84 is the outward normal of the boundary .
4 .   assign @xmath85 .
( 4,0 ) ( -0.5,0)(1,0)12
in this section , we introduce two new eulerian algorithms for flow visualizations . the first algorithm is to improve the eulerian numerical method for computing the forward flow map . instead of solving the pde _
backward _ in time as in section [ subsec : oldeulerian ] and @xcite , we develop an eulerian pde algorithm to construct the forward flow map from @xmath2 to the final time @xmath3 so that we do not need to load the time - dependent velocity field from the terminal time level backward to the initial time . in other words , we construct the forward flow map _ on the fly_. this is computationally more natural and convenient .
because of this simple algorithm , we propose to develop an efficient and systematic pde based method for computing the isle field .
typical challenge in obtaining such field is that the choice of the separation factor @xmath0 is in general flow dependent .
one has to vary @xmath0 to extract the necessary information .
usual lagrangian methods unfortunately require shooting different sets of initial rays for each individual separation factor @xmath0 .
this is computationally very inefficient .
the proposed eulerian formulation requires either some simple interpolations or only one single level surface extraction which can be easily implemented using any well - developed contour extraction function such as in .
one disadvantage about the previous eulerian approach for forward flow map computation as discussed in section [ subsec : oldeulerian ] and @xcite is that the level set equation @xmath86 has to be first solved _
backward _ in time from @xmath3 to @xmath2 .
once we have the final solution at the time level @xmath2 , one can then identity the flow map @xmath76 by @xmath87 .
this could be inconvenient especially when we need to access the intermediate forward flow map . in this section
, we propose a new algorithm to construct the forward flow map _ on the fly_.    consider the forward flow map from @xmath2 to @xmath3 .
suppose the time domain @xmath88 $ ] is discretized by @xmath89 discrete points @xmath90 , where @xmath91 and @xmath92 .
then on each subinterval @xmath93 $ ] for @xmath94 , we use the same method as described in section [ subsec : oldeulerian ] to construct the forward flow map @xmath49 .
in particular , we solve the level set equation ( [ eqn : levelset ] ) _ backward _ in time from @xmath95 to @xmath96 with the terminal condition @xmath97 imposed on the time level @xmath95 .
then the forward flow map is given by @xmath98 .
once we have obtained this one step _
forward _ flow map @xmath49 , the _ forward _ flow map from @xmath2 to @xmath3 can be obtained using the composition @xmath99 . and this can be easily done by typical numerical interpolation methods . to prevent extrapolation
, we could enforce @xmath100 \times [ y_{\min},y_{\max}]$ ] .    here
we emphasize that even though we solve equation ( [ eqn : levelset ] ) _ backward _ in time for each subinterval , we access the velocity field data @xmath101 prior to @xmath102 and therefore the forward flow map is indeed obtained _ on the fly_. as a simple example , we consider the first order euler method in the temporal direction .
the forward flow map @xmath49 can be obtained by first solving @xmath103 with the terminal condition @xmath104 , and then assigning @xmath105 .
higher order generalization is straight - forward .
for example , if we use the tvd - rk2 method in the temporal direction @xcite , we can compute the forward flow map @xmath49 by @xmath106 with the terminal condition @xmath104 , and the functions @xmath107 and @xmath108 are two intermediate predicted solutions at the time levels @xmath109 and @xmath90 , respectively .
then the flow map from the time level @xmath96 to @xmath95 is given by @xmath110 .
it is rather natural to decompose the flow map @xmath111 into a composition of maps @xmath112 .
for example , such idea has been used recently in @xcite to improve the computational efficiency of the lagrangian ftle construction between various time levels . define the interpolation operator by @xmath113 , which returns the interpolated mesh values from the non - mesh values .
then the flow map @xmath48 can be decomposed into @xmath114 .
if the maps @xmath115 for @xmath116 are all stored once they are computed , one can form @xmath48 by determining only @xmath117 .
the idea in this work shares some similarities with what we are implementing .
our approach can also re - use all intermediate flow maps , the work @xcite however has concentrated only on improving the computational efficiency of the lagrangian approach .
we are using the idea of flow map decomposition to propose a forward computational strategy for the eulerian formulation .    to obtain a stable evolution in the flow map constructions ,
theoretically we require the interpolation scheme between two local flow maps to be monotone , i.e. the interpolation scheme should preserve the monotonicity of the given data points @xcite .
this is especially important in the context of the phase flow maps for autonomous flows @xcite since the overshooting due to the interpolation could be significantly amplified in the time doubling strategy @xmath118 ^ 2 \ , .\ ] ]      in this section , we propose an efficient eulerian approach to compute the isle function for an arbitrary separation factor @xmath0 based on the techniques developed in the previous section .    to compute the isle at @xmath119 , one has to determine the minimum time @xmath120 at each location @xmath119 for which @xmath121 .
it is , therefore , necessary to keep track of all intermediate values @xmath122 for all @xmath123 s during the flow map construction . at each grid point
@xmath124 and each intermediate time step @xmath90 , we construct a quantity @xmath125 by @xmath126\ ] ] with @xmath127 .
once we determine the flow map @xmath128 , we have also constructed an increasing sequence at each grid point @xmath124 in time @xmath129 , i.e. @xmath130 for each @xmath124 .
now , we interpret this quantity as a spatial - temporal volumetric data in the @xmath1 space .
due to such monotonicity property in the temporal direction , the isosurface @xmath131 forms a graph of @xmath119 , i.e. for each @xmath119 we have a unique @xmath132 so that @xmath133 . and more importantly , the value @xmath134 gives the shortest time so that @xmath135 .
this implies that @xmath136 and so the level surface in the extended phase space , i.e. @xmath137 can be used to define @xmath138 .
finally , the isle function can then be computed using @xmath139 the computational algorithm is summarized in algorithm 2 .
( 4,0 ) ( -0.5,0)(1,0)12     +    algorithm 2 ( our proposed eulerian method for computing the isle ) :    1 .
discretize the computational domain @xmath140 2 .
set @xmath141 .
3 .   for @xmath94 , 1 .   solve the level set equation ( [ eqn : levelset ] ) for @xmath142 with the condition @xmath143 .
2 .   define @xmath144 .
3 .   interpolate the solutions to obtain @xmath145 .
4 .   compute @xmath146 and determine @xmath147 \ , .\ ] ] 4 .
given an arbitrary separation factor @xmath0 .
for each @xmath124 , search for a time step @xmath90 such that @xmath148 . then determine @xmath138 by applying the piecewise linear interpolation using the values @xmath125 and @xmath149 . + the isle function is computed as @xmath150    ( 4,0 ) ( -0.5,0)(1,0)12     + there are several major advantages of the proposed eulerian approach .
it provides an efficient yet simple numerical implementation of @xmath151 for multiple separation factors @xmath0 .
the construction of the quantity @xmath152 can be first done in the background . for each particular given value of @xmath0
, one has to perform only one single isosurface extraction in the extended phase space , i.e. the @xmath1 space . unlike the typical lagrangian implementation where one has to repeatedly keep track of _
all _ particle trajectories for the constraint @xmath153 , the proposed eulerian algorithm provides a systematic way based on a simple thresholding strategy .
note that in our implementation , if it happens that @xmath154 for all @xmath90 , we simply set @xmath155 .
we conclude this section by discussing the computational complexity of our proposed algorithm for a two dimensional flow , i.e. @xmath156 .
let @xmath157 and @xmath54 be the discretization size of one spatial dimension and time dimension respectively .
since the liouville equation is a hyperbolic equation , we have @xmath158 by the cfl condition . at each time
step @xmath90 , a short time flow map @xmath49 is obtained by solving the liouville equations from @xmath159 to @xmath160 , the computational effort is @xmath161 .
then in computing the long time flow map @xmath162 , an interpolation @xmath163 takes @xmath161 operations .
therefore the construction of flow map at each time step requires @xmath164 operations .
computing @xmath165 and @xmath149 also needs @xmath161 operations , the complexity order is kept at @xmath161 . summing up this procedure in all time steps ,
the overall computational complexity is @xmath166 .
this new eulerian method has significantly improved the overall computational complexity in the application to the isle computations .
since the calculations involve the forward flow map from @xmath167 to @xmath168 for all @xmath169 , the original eulerian approach as discussed in section [ subsec : oldeulerian ] requires to solve the liouville equation from each individual time level @xmath168 backward in time to @xmath167 . the overall computational complexity in obtaining all flow maps
is therefore given by @xmath170 which is one order magnitude larger than the newly developed eulerian method .
with algorithms proposed in section [ sec : proposal ] , we can finally develop an efficient numerical approach to compute the _ forward _ ftle field by solving corresponding pdes _
forward _ in time and the isle field for arbitrary separation factor @xmath0 .
the ftle and the isle ( or the fsle ) are indeed measuring different properties of a given flow , even though both quantities depend on the same function @xmath171 which relates the growth of an infinitesimal perturbation .
in particular , ftle measures how much @xmath171 grows over a finite time span @xmath172 $ ] , while the isle measures the time required for the quantity @xmath171 to reach a threshold value prescribed by the separation factor .
they are two different tools to study how rapid @xmath171 grows
. a careful comparison of these two quantities can be found in , for example , @xcite .    in many numerical experiences , on the other hand
, the ftle fields can show striking visual resemblance with the isle fields of certain separation factors .
however , there has not been much studies on the theoretical relationship between these two quantities . in this section
, we are going to demonstrate how these two quantities relate to each other by considering their corresponding _
ridges_. of course , the way to define a _
ridge _ is not unique and one can pick a convenient definition for a particular application @xcite . in this work , we adapt the following simple definition .    for a given scalar function @xmath173 , we define a @xmath174-ridge to be a codimension - one compact and @xmath175 surface , denoted by @xmath176 , such that for every point @xmath177 , @xmath173 is a local maximum along the normal direction to @xmath176 at @xmath119 .    at these locations , @xmath173 is also called a generalized maximum in some literatures .
next , we define a particular neighbourhood of a ridge for later discussion .
suppose @xmath178 is a codimension - one , compact and smooth surface .
we define @xmath179 , the tubular neighbourhood of @xmath178 with radius @xmath180 , as the collection of any point @xmath181 that can be uniquely expressed as @xmath182 , for some point @xmath183 , a real number @xmath184 , and @xmath185 is an unit normal vector of @xmath178 at @xmath119 .    to ensure that any point @xmath181 on @xmath179 can be uniquely expressed as such
, one can equivalently require that the boundary of @xmath179 does not intersect itself .
for a smooth surface , this can always be guaranteed by choosing a small enough positive number @xmath186 .
the definitions of both the ftle and the isle involve @xmath187 .
it is therefore natural to link them through this particular quantity .
the following simple theorem which identifies the @xmath188-ridges with the @xmath189-ridges , will be the very first step in building up our analysis .
[ prop : ftle&lambda ] let @xmath188 and @xmath189 be the ftle function and the largest eigenvalue of the cauchy - green deformation tensor from time @xmath190 to time @xmath18 , respectively . then
, a @xmath188-ridge is also a @xmath189-ridge .
since @xmath191 is a just a scaling factor and the natural logarithm function is strictly increasing , all generalized maxima of @xmath188 as defined in ( [ eqn : ftle ] ) and @xmath189 coincide .    for convenience , in the time window @xmath172
$ ] , we impose @xmath192 if @xmath193 for all
@xmath194 $ ] .
if an isle function @xmath195 on an open set @xmath9 of the spatial domain , there is obviously no visible isle feature on @xmath9 , corresponding to this separation factor @xmath0 .
the following theorem establishes a fundamental result that for certain choices of @xmath0 , the isle function @xmath196 can be positive on a tubular neighbourhood of a @xmath189-ridge .
this explains the non - trivial isle features around the location of a ftle ridge , when one chooses a suitable separation factor @xmath0 .
[ prop : funprop1 ] let @xmath189 be the largest eigenvalue of the cauchy - green deformation tensor from time @xmath190 to time @xmath18 .
suppose @xmath189 is continuous everywhere in the spatial domain , and @xmath176 is a @xmath189-ridge with @xmath197 .
then for any positive number @xmath198 , there is a tubular neighbourhood @xmath199 such that the isle function is positive on @xmath199 .    for each @xmath177 , using the definition of @xmath200 , we have @xmath201
. then by the continuity of @xmath189 at @xmath119 , we can find an open ball @xmath202 with center @xmath119 and radius @xmath203 such that @xmath204 for all @xmath205 . collect all these open balls at every point on @xmath176 into @xmath206 , which is clearly an open cover of @xmath176 . by the compactness of @xmath176
, we can find a finite subcover @xmath207 from @xmath206 .
take @xmath208 . without loss of generality ,
assume @xmath186 is small enough such that the tubular neighbourhood @xmath199 is contained in the union of sets in @xmath209 .
now for any @xmath210 , it is clear that @xmath211 . because the function @xmath212 is a continuous function of @xmath213 , we apply the intermediate value theorem and conclude that it must have reached @xmath0 at some @xmath214 $ ] .
therefore @xmath215 is positive .
the next theorem is to identify a @xmath189-ridge with a @xmath196-ridge of a certain isle function @xmath196 . to justify the condition ,
the points on the location of a @xmath189-ridge , are expected to exhibit highly chaotic behavior .
it is reasonable to assume that we can find a tubular neighbourhood of a @xmath189-ridge , where the time monotonicity of @xmath216 holds for @xmath217 $ ] , and that no other points reach the minimum value of @xmath218 , except those on @xmath176 .
[ prop : funprop2 ] let @xmath189 be the largest eigenvalue of the cauchy - green deformation tensor from time @xmath190 to time @xmath18 .
suppose @xmath176 is a @xmath189-ridge with @xmath219 .
if there exists a tubular neighbourhood @xmath199 of @xmath176 such that    1 .   for any fixed @xmath210 , the quantity @xmath220 is increasing in @xmath217 $ ] , and 2 .
@xmath221 ,    then @xmath176 is also a @xmath222-ridge for the isle function @xmath222 .
let @xmath177 , and @xmath223 be a point in the normal direction of @xmath176 at @xmath119 .
by the second condition , we have @xmath224 . and
since @xmath220 is increasing with @xmath213 , it follows that @xmath225 for @xmath217 $ ] .
therefore @xmath226 by our convention .
we can conclude that no other points on the normal direction of @xmath176 at @xmath119 has positive value of @xmath222 , except @xmath119 itself .
this makes @xmath119 a generalized maximum of @xmath222 .
theorem [ prop : funprop2 ] is a restricted result , in the sense that the whole @xmath189-ridge is preserved for a particular @xmath196 function .
but it is straightforward to extend it to the following corollary , which looks for a certain portion of a @xmath189-ridge that can be conserved in a larger class of isle functions .
[ prop : funcor ] let @xmath189 be the largest eigenvalue of the cauchy - green deformation tensor from time @xmath190 to time @xmath18 .
suppose @xmath176 is a @xmath189-ridge with @xmath219 , and @xmath227 is a connected subset of @xmath176 with @xmath228 .
for any positive number @xmath229 $ ] , if there exists a tubular neighbourhood @xmath230 of @xmath227 such that    1 .
for any fixed @xmath231 , the quantity @xmath220 is increasing in @xmath217 $ ] , and 2 .
@xmath232 ,    then @xmath227 is a @xmath196-ridge for the isle function @xmath196 .
the proof is parallel to that of theorem [ prop : funprop2 ] , with @xmath176 and @xmath200 being substituted by @xmath227 and @xmath0 , respectively .    through theorem [ prop : funprop1 ] to corollary [ prop : funcor ] , we have obtained a close relationship between the @xmath189-ridges and the isle ridges . by the equivalence of the ftle ridges and the @xmath189-ridges
as proved in theorem [ prop : ftle&lambda ] , we can then link the ftle ridges and the isle ridges . in practice , the estimated location of the isle ridges can be inferred by the location of the ftle ridges , while @xmath189 plays an important role for estimating a suitable separation factor defined in the isle ridge . to demonstrate this ,
suppose that there is a ftle ridge @xmath176 satisfying all conditions in theorem [ prop : funprop2 ] .
then @xmath176 is an isle ridge with the separation factor @xmath233 .
although the value of @xmath200 is actually an unknown , we can determine a number @xmath234 that is slightly less than @xmath235 manually . from the definition of the ftle
, we have @xmath236 > e^{lt } \ , .\ ] ]    now , with any separation factor @xmath237 , theorem [ prop : funprop1 ] implies that there is a tubular neighborhood around @xmath176 with positive isle values .
we can then increase the separation factor from @xmath238 approaching to @xmath200 , so that the tubular neighborhood containing @xmath176 becomes narrower , and thus leads to a better approximation of @xmath176 . in practice
, the value of @xmath234 has to be chosen in a trial - by - error fashion and has to be determined case - by - case .
based on the eulerian algorithm we developed in section [ subsec : neweulerian ] and section [ subsec : islecomp ] , we are now able to efficiently construct the isle field of each individual separation factor using one single isocontour extraction .
in this section , we present numerical results to show the feasibility of the proposed eulerian approach .
we will compare the solution from the new approach with the previous eulerian approach discussed in section [ subsec : oldeulerian ] .
it is noteworthy that for a dynamical system in a fixed time window , we can now construct the isle function with an arbitrary choice of the separation factor in only one single numerical flow simulation .
therefore , we can easily verify the theoretical connection between ftle and isle , developed in section [ sec : relation ] .
this first example is a simple flow taken from @xcite to describe a periodically varying double - gyre .
the flow is modeled by the following stream - function @xmath239 \sin(\pi y ) \
, , \ ] ] where @xmath240 in this example , we follow @xcite and use @xmath241 , @xmath242 .
we discretize the domain @xmath243\times[0,1]$ ] using 513 grid points in the @xmath244-direction and 257 grid points in the @xmath245-direction .
this gives @xmath246 .    in figure
[ fig : ftle ] , we compare the ftle field @xmath247 of the double - gyre flow with @xmath248 computed using the lagrangian approach and our proposed eulerian approach , respectively .
these two solutions match extremely well .
we have also checked the @xmath249-error of the forward flow map @xmath250 using different @xmath251 s ranging from @xmath252 to @xmath253 .
since we do not have the exact solution for this flow , the _ exact _ solutions are computed using the lagrangian ray tracing with a very small time step .
figure [ fig : convergence ] shows the errors in the flow map @xmath254 with respect to different @xmath251 s .
we find that the flow map computed using the proposed eulerian approach is approximately second order accurate .
figures [ fig : isle1 ] shows the isle fields @xmath255 with several separation factors @xmath0 . in the implementations ,
the flow map is computed from @xmath2 up to @xmath256 . recall that we impose zero isle value at points whose local separation rate is less than the given separation factor . with the increasing separation factor ,
the regions with non - zero isle values , get smaller .
nevertheless , the regions with non - zero or high isle values , always concentrate near the location of the ftle ridges shown in figure [ fig : ftle ] .
we denote by @xmath176 the most prominent ridge originating from the central bottom of the domain .
note that a rough estimate of the minimum ftle value on @xmath176 is approximately 0.3 .
and so we have the estimate @xmath257 . using theorem [ prop : funprop1 ] ,
we can always find a tubular neighborhood near @xmath176 of non - zero isle using a separation factor @xmath0 smaller than approximately 20 . in figure
[ fig : isle1 ] , we have computed the isle using various separation factors from @xmath258 to @xmath259 . as @xmath0 increases , we can see that there are indeed non - zero isle regions around the major ridge @xmath176 . and
evidently , the isle ridges lie within these regions .    to better compare the ftle field and the isle field , we plot these two solutions together in figure [ fig : isosurface ] . the @xmath244- and the @xmath245- axis represent the computational domain @xmath9 while the @xmath260 axis denotes the temporal direction . in ( a ) , we plot the ftle field @xmath261 at @xmath262 and @xmath263 at @xmath264 . in ( b - d )
, we plot the isosurfaces of the isle field @xmath265 with @xmath266 , respectively .
the @xmath260-values denote the time required for the local separation ratio to achieve @xmath0 for the first time , i.e. the @xmath260-values are the @xmath267 we defined before . near the major ftle ridge
, we can see the values of @xmath120 are significantly lower than the remaining part of the domain .
therefore the corresponding isle values are generalized maxima .
this shows us the isle ridges match with the ftle ridges very well .
we consider a simple analytical field @xcite given by @xmath268 and @xmath269 in this example .
the computational domain is @xmath270 ^ 2 $ ] and we discretize the domain with @xmath271 .
we impose the fixed inflow boundary condition and non - reflective outflow boundary condition as in @xcite .
figure [ fig : ftle2 ] shows the ftle field @xmath272 computed using our original eulerian approach and the proposed eulerian approach .
it is clear that the new approach can successfully handle the inflow and outflow of particles on the boundaries , just like the previous approach does .
next , we are going to apply the results in section [ sec : relation ] to analyze the isle ridges .
figure [ fig : ftle2 ] suggests that there is only one ftle ridge . and
the minimum value of @xmath273 on the ridge is close to @xmath274 . by the definition of ftle
, we have @xmath275 \ , .\ ] ] therefore we can conclude that @xmath276 is a reasonable estimate of the minimum value of @xmath277 on the ridge .
we have tested with several separation factors that are near @xmath278 , and the results can be plotted in figure [ fig : isle2a ] . to better visualize the isle field , we plot in figure [ fig : isle2 ] the time @xmath267 required for different separation factors @xmath0 .
the shorter the time , the bigger the corresponding isle .      [ ! htb ]    we consider the dynamical system governed by a duffing and a van der pol oscillator , as in @xcite .
the system is given by @xmath279    the computational domain is @xmath280\times [ -1.5,1.5]$ ] , with mesh size @xmath281 .
we simulate the flow from @xmath2 to @xmath256 .
the discretization size is rather coarse for such a complicated dynamics , but our new approach can still demonstrate the fine details of the transport barrier as shown in figure [ fig : ftle3 ] .
there are several major ftle ridges that are easily spotted .
the minimum value of @xmath282 on those ftle ridges are very close to @xmath283 .
so @xmath284 .
theorem [ prop : funprop1 ] implies that the separation factor @xmath285 can approximately single out the isle ridges .
we pick several separation factors that are close to @xmath286 , and we find that the isle ridges corresponding to @xmath287 , as shown in figure [ fig : isle3](c ) , have the greatest resemblance to the ftle ridges in figure [ fig : ftle3 ] .      [ !
htb ]    to further demonstrate the effectiveness and robustness of our proposed approach , we test our algorithm on a real life dataset obtained from the ocean surface current analyses real - time ( oscar ) .
the oscar data were obtained from jpl physical oceanography daac developed by esr .
it measures ocean surface flow in the region from @xmath288 to @xmath289 in latitude and from @xmath290 to @xmath291 longitude .
the time resolution is about 5 days and the spatial resolution is @xmath292 in each direction .     in this numerical example
, we consider the area near the line islands in the region from @xmath293 to @xmath294 in latitude and from @xmath295 to @xmath296 longitude . to have a better visualization of both the ftle and isle field , we interpolate the velocity data to obtain a finer resolution of @xmath297 days in the temporal direction and @xmath298 in each spatial direction
. then we look at the ocean surface current within the first 50 days in the year 2014 to obtain the forward ftle field @xmath299 and the isle field @xmath300 .
as demonstrated in figure [ fig : oscar_ftle ] , our new approach works well in extracting fine details of the transport barrier even for real data .
there are a few major ftle ridges spotted in the figure .
according to theorem [ prop : funprop1 ] , we found that the separation factor @xmath0 roughly equals to @xmath301 to @xmath302 can single out the isle ridges .
we have tested several separation factors that are close to @xmath303 , as shown in figure [ fig : oscar_isle ] .
we find that the isle ridges corresponding to @xmath259 , in figure [ fig : oscar_isle](b ) , have the greatest resemblance to the ftle ridges in figure [ fig : oscar_ftle ] . as a comparison
, we have also plotted the isle field @xmath300 for @xmath304 and @xmath305 in figure [ fig : oscar_isle ] ( a ) and ( c ) , respectively .
these results verify the theory developed in section [ sec : relation ] .
to summarize , we emphasize the novelty and importance of both the two eulerian algorithms we developed in this paper .
based on the first algorithm , we are now able to determine the _ forward _ flow map _ on the fly _ so that the pde is solved _ forward _ in time . unlike the original eulerian algorithm
as developed in @xcite , there is no need to store the whole velocity field at all time steps at the beginning of the computations .
this makes the forward flow map computations practical .
more importantly , because of the first algorithm , we are now able to access all intermediate forward flow maps which is necessary in computing the isle field computations .
this eulerian approach is extremely efficient in computing the isle field . to the best of our knowledge ,
this is the first efficient numerical approach for the isle computations .
theoretically , this paper provides a rigorous mathematical analysis to explain the similarity between the ftle ridges and the isle ridges based on @xcite , which has been verified by the proposed eulerian approach .
the work of you was supported by the talents introduction project of nanjing audit university , and the natural science foundation of jiangsu higher education institutions of china ( no.16kjb110012 ) .
the work of leung was supported in part by the hong kong rgc grants 16303114 and 16309316 .                                                                          c.  w. shu .
essentially non - oscillatory and weighted essentially non - oscillatory schemes for hyperbolic conservation laws . in b.  cockburn , c.  johnson , c.w .
shu , and e.  tadmor , editors , _ advanced numerical approximation of nonlinear hyperbolic equations _ , volume 1697 , pages 325432 .
springer , 1998 .
lecture notes in mathematics .
w.  tang , p.w .
chan , and g.  haller .
accurate extraction of lagrangian coherent structures over finite domains with application to flight data analysis over hong kong international airport .
, 20:017502 , 2010 .


lattice simulations are used in many areas of quantum physics , ranging from nuclear and particle physics to atomic and condensed matter physics  @xcite . in such calculations
the system is solved numerically using a discrete spacetime lattice over a finite volume . in practice
this finite volume is usually taken to be a cubic box with periodic boundaries .
when simulating composite objects such as bound states , the boundaries of the periodic box will modify quantum wave functions .
this leads to finite volume shifts in the binding energies , and detailed knowledge of these finite volume effects is necessary to improve high precision lattice calculations .    in ref .
@xcite , lscher derived a formula for the finite volume mass shift of s - wave bound states of two particles with reduced mass @xmath1 interacting via a potential with finite range @xmath2 .
when a bound state with energy @xmath3 is put in a periodic cubic box of length @xmath4 , its energy in the rest frame is shifted by an amount @xmath5 where @xmath6 is the binding momentum and @xmath7 is the asymptotic wave function normalization defined by @xmath8 for @xmath9 . for potentials with exponential fall - off , @xmath10 for large @xmath11 ,
the formula is modified by exponentially small corrections when the binding momentum @xmath12 is smaller than @xmath13 .
bound states moving in a finite periodic volume also have a topological phase correction to the energy  @xcite .
this topological phase correction contains information about the number and mass of the constituents of the bound states .
it must be included when determining scattering phase shifts for composite objects in a finite volume .    in ref .
@xcite , we briefly discussed the generalization of lscher s formula ( [ eq : delta_m - luscher ] ) for the finite volume mass shift to higher partial waves . in this paper , we present the full derivation of these results . moreover , we give the general mass shift for states with angular momenta up to @xmath14 . in general ,
the mass shift for a given state depends on its transformation properties with respect to the symmetry group of the cubic box .
in addition to reducing finite volume effects from precision lattice calculations , our finite volume results can also be used as a diagnostic tool to probe the angular momentum and radial structure of the bound state wave function .
a summary of our results is given by eq .
together with table  [ tab : results ] at the end of sec .
[ sec : results ] .
we also derive a general mass shift formula for two - dimensional systems in a finite area with periodic boundaries .    in our discussion here
we only consider stable bound states .
however , resonances above continuum thresholds can also be described from the volume dependence of avoided level crossings  @xcite .
our results are universal and can be applied to a wide range of systems .
we will discuss a few examples below .    in particle physics
, there is some interest in hadronic molecules with angular momentum  @xcite . in the case of s - waves , the deuteron and some exotic weakly bound states such as the h - dibaryon
were recently studied in lattice qcd  @xcite .
similar investigations for exotic bound states with angular momentum appear feasible in the future . in atomic physics ,
several experiments have investigated strongly - interacting p - wave feshbach resonances in @xmath15li and @xmath16k  @xcite , which can be tuned to produce bound p - wave dimers .
if such systems are simulated in a finite volume , our results can be used to describe the finite volume dependence of the dimers .
halo nuclei constitute a special class of weakly - bound nuclei with molecular character .
they consist of a compact core and one or more nucleons with low separation energy  @xcite . among halo nuclei
there are some systems with nonzero orbital angular momentum .
a well - known example of a p - wave halo state is the @xmath17 excited state in @xmath18be .
the electromagnetic properties of the low - lying states in this nucleus can be well described in a two - body halo picture of a @xmath19be core and a neutron  @xcite .
a similar strategy was applied to radiative neutron capture on @xmath20li  @xcite . a related class of systems is given by nuclei with an @xmath21-cluster structure such as @xmath22be and excited states of @xmath23c  @xcite .
the hoyle state in @xmath23c was recently calculated for the first time in an _ ab initio _ approach using chiral effective field theory and nuclear lattice calculations  @xcite .
understanding the finite volume corrections is crucial in such a calculation . due to the @xmath21-cluster structure
, the volume dependence of the hoyle state is governed by the volume dependence of a three - body system .
finally , we note that the asymptotic normalization coefficient of the bound state wave function appears in the mass shift formula .
hence , our results can be used to extract this quantity from lattice calculations at finite volume .
the asymptotic normalization is directly connected to zero - energy capture reactions  @xcite , which play an important role in nuclear astrophysics .
the paper is organized as follows .
we start with some prerequisites and a general discussion of the finite volume mass shift in section  [ sec : boundstates ] .
lscher s result for s - waves is reviewed in section [ sec : s - wave ] , while our extension to higher partial waves is given in section [ sec : p - wave - and - higher ] .
in particular , we discuss the mass shift for the irreducible representations of the cubic group , relate the sign of the shift to the leading parity , and derive a trace formula for the multiplet - averaged mass shift for states with arbitrary angular momentum @xmath24 . in section [ sec : numtest ] , we verify our results numerically for two model systems .
the case of two spatial dimensions is treated in section  [ sec:2d ] .
finally , a brief summary and outlook is given in section [ sec : conclusion ] .
in the appendix , we discuss the relation of the asymptotic normalization constant to scattering parameters and derive explicit shallow binding relations between the effective range , the binding momentum and the asymptotic normalization coefficient .
as a starting point , we first review several results from  @xcite . we closely follow lscher s derivation , but consider a more general system with arbitrary angular momentum and non - local interactions .
we start with some definitions and basic identities about bound states in rotationally symmetric potentials @xmath25 . this will also help to define our notation used throughout the discussion . for a bound state with energy @xmath26 and
angular - momentum quantum numbers @xmath27 , we separate the wave function @xmath28 as a product of the radial wave function @xmath29 and spherical harmonics @xmath30 , @xmath31 this leads to the radial schrdinger equation @xmath32u_\ell(r ) = 0 \ , , \label{eq : sg - rad}\ ] ] with @xmath33 and @xmath34 .
the mass is denoted by @xmath1 since in later applications it will be the reduced mass for a two - particle system .
the normalized radial wave function satisfies @xmath35    for vanishing potential , we have two linearly independent solutions of the free wave equation .
these are the riccati  bessel functions @xmath36 and riccati  neumann functions @xmath37 , which can be expressed in terms of ordinary bessel functions .
it is often more convenient to work with the riccati  hankel functions @xmath38 these have the asymptotic form @xmath39 as @xmath40 . for future reference , we give the explicit expressions for @xmath41 for @xmath42 :    @xmath43      in order to discuss the asymptotic form of bound state wave functions we define solutions @xmath44 of that fulfill the condition @xmath45 the bound - state solution normalized according to can then be written as @xmath46 where @xmath47 and @xmath48 is the asymptotic normalization coefficient .
the factor @xmath49 adjusts the phase such that @xmath29 is a real function . if the potential has a finite range @xmath2 , _
i.e. _ , @xmath50 for @xmath9 , we have the exact identity @xmath51    the asymptotic normalization @xmath7 is itself an interesting quantity . in low - energy astrophysical reactions
it determines the capture rate of a single proton or neutron at zero relative energy  @xcite . in the shallow binding limit
it is also directly related to low - energy scattering parameters . for @xmath52 ,    [ eq : gamma - r ] @xmath53 whereas for angular momenta @xmath54 one has @xmath55    we note that @xmath7 will in general depend upon the binding momentum @xmath12 .
we will derive the relations   and   in appendix  [ sec : shallowbinding ] .
the s - wave formula is well known and given in  @xcite , although without the @xmath56 correction estimate . for higher angular momentum
the situation is qualitatively different due to the divergence of @xmath57 at @xmath58 .
we now consider a system of two spinless particles with reduced mass @xmath1 with zero total momentum .
we will work with relative coordinates @xmath59 .
the interaction is given by a potential @xmath60 , which we in general allow to be non - local in configuration space .
it is assumed to be hermitian , rotationally symmetric , and to have a finite range @xmath2 , _
i.e. _ , @xmath61    we consider the case when the system has a bound state @xmath62 with angular quantum numbers @xmath27 .
the schrdinger equation , @xmath63 can be written as @xmath64 in configuration space .
we note that for a local potential , @xmath65 eq .
reduces to the familiar form @xmath66\,\psi_b({\mathbf{r } } ) = -{e_b}\,\psi_b({\mathbf{r } } ) \ , .
\label{eq : sg - local}\ ] ] regardless of the locality of the interaction , the wave function @xmath67 has the asymptotic form @xmath68 where @xmath6 is the binding momentum .
we now consider what happens when the two - body system is put into a cubic periodic box with length @xmath69 . for this problem
it is convenient to define a periodic extension of the potential @xmath70 we take @xmath71 to be an exact periodic solution of the finite volume schrdinger equation , @xmath72 with the finite - volume hamiltonian @xmath73 .
it is clear that @xmath74 approaches the infinite volume eigenvalue @xmath75 and @xmath76 as @xmath77 .
we now derive a formula for the finite volume mass ( energy ) shift , @xmath78 let us define a state @xmath79 by adding together periodic copies of the infinite volume wave function in , @xmath80 this clearly satisfies the periodicity condition . acting upon this state with the finite volume hamiltonian , we get @xmath81 the final result can be written as @xmath82 where we have defined @xmath83 as @xmath84 with the substitution @xmath85 for each term in the sum , this can be rewritten as @xmath86 due to the finite range of the potential we only get contributions from the domain @xmath87 .
we note that @xmath88 when @xmath89 and @xmath90 .
therefore we can use the asymptotic form of the wave function and find that @xmath91 .
this means that @xmath79 is an approximate solution of the finite volume schrdinger equation   for large @xmath4 .
motivated by this , we write the exact finite volume solution @xmath71 explicitly as @xmath92 we take @xmath71 to be unit - normalized per volume @xmath93 .
the same is true of @xmath79 up to corrections of order @xmath94 .
we will choose @xmath21 such that @xmath95    consider now the matrix element @xmath96 .
acting with @xmath97 on @xmath79 , we get @xmath98 according to and . on the other hand ,
acting with @xmath97 on @xmath99 yields @xmath100 combining these two results we find @xmath101    we first consider the numerator in this expression .
obviously , @xmath102 we note that the factor of @xmath21 here will cancel the @xmath21 in the denominator of  .
we can now simplify further starting with @xmath103 for each @xmath104 we can make the substitutions @xmath105 and @xmath85 .
these leave the integrals invariant , and we get @xmath106 setting @xmath107 and @xmath108 yields @xmath109 where @xmath110 counts the number of repeated periodic copies .
the fact that @xmath110 diverges simply reflects the fact that we are working with periodic wave functions with normalization measured per volume @xmath93 , and @xmath110 will cancel in the final result . for the integral to be non - zero
, both @xmath111 and @xmath112 have to be close to @xmath113 due to the finite range of the potential . from the assumption
@xmath69 it then follows that all terms with @xmath114 are suppressed by at least a factor of @xmath115 , and we have @xmath116 the possible nonvanishing values of @xmath117 are @xmath118 .
we therefore arrive at @xmath119 for the denominator in , an analogous procedure yields @xmath120 \label{eq : psi0-psi0-final}\ ] ] with the same constant @xmath110 as above . combining   and  , the constant @xmath110 cancels and we get @xmath121 where we have renamed @xmath122 back to @xmath104 .
is a general result valid for any angular momentum .
the dependence of the mass shift upon quantum numbers @xmath27 will emerge from the wave function @xmath123 and the resulting overlap integrals in  . in the following ,
we explore this dependence in detail and denote the mass shift as @xmath124 .
for @xmath52 the asymptotic wave function   is given as    @xmath125    with @xmath126    due to the finite range @xmath90 of the potential we only have contributions with @xmath127 in  .
hence we can insert the asymptotic form for @xmath128 and get @xmath129 we can furthermore use the schrdinger equation   to eliminate the potential . doing this and then renaming @xmath130 , we get @xmath131   \psi_b^*\big(|{\mathbf{r}}|\big)\right\ }   \frac{{\mathrm{e}}^{-\kappa|{\mathbf{r}}+{\mathbf{n}}l|}}{|{\mathbf{r}}+{\mathbf{n}}l| }   + { \mathcal{o}}\big({\mathrm{e}}^{-{\sqrt2\kappa l}}\big ) \\   & = \frac\gamma{\sqrt{4\pi}}\sum\limits_{|{\mathbf{n}}|=1 }   \int{\mathrm{d}^3r}\,\psi_b^*\big(|{\mathbf{r}}-{\mathbf{n}}l|\big )   \frac{1}{2\mu}\left[{\delta}_r-\kappa^2\right]\frac{{\mathrm{e}}^{-\kappa r}}{r }   + { \mathcal{o}}\big({\mathrm{e}}^{-{\sqrt2\kappa l}}\big ) \ , . \end{split}\ ] ] in the second line we have shifted the integration variable and used partial integration to let the laplacian act on @xmath132 . finally , we use the fact that @xmath133 is a green s function for the operator @xmath134 , @xmath135\frac{{\mathrm{e}}^{-\kappa r}}{4\pi r }   = -{\delta^{(3)}}({\mathbf{r } } ) \ , .
\label{eq : helmholtz - gf}\ ] ] this allows us to perform the integral and arrive at @xmath136 in the last step we have inserted the asymptotic form of the wave function for @xmath137 , and the sum yields a factor of six .
this is just lscher s result   as given in the introduction .
we now discuss the generalization of the mass shift formula to arbitrary angular momentum . the general form for the asymptotic wave function
is @xmath138 inserting this into   and performing steps analogous to those presented above for the s - wave case , we find @xmath139\!\psi_b^*({\mathbf{r}}-{\mathbf{n}}l)\right\}\ ,   { y}_\ell^m(\theta,\phi)\frac{{\mathrm{i}}^\ell\gamma\hat{h}_\ell^+({\mathrm{i}}\kappa r)}r   + { \mathcal{o}}\big({\mathrm{e}}^{-{\sqrt2\kappa l}}\big)\ , .
\label{eq : delta_m - int - asympt}\ ] ]    the crucial ingredient is the relation @xmath140 \ , , \label{eq : sh - hankel}\ ] ] where @xmath141 are the solid harmonics defined via @xmath142 .
the derivation of follows from lemma b.1 in ref .
@xcite , which proves that @xmath143 for any smooth function @xmath144 .
we obtain   by using the relation , @xmath145 and noting that @xmath146 and @xmath147 , where @xmath148 is a spherical hankel function of the first kind .
we can illustrate using an example .
for the case @xmath149 and @xmath150 we have @xmath151 and @xmath152 .
a straightforward calculation shows that indeed @xmath153 \,,\ ] ] with @xmath154 .    using   to rewrite
, we get @xmath155\!\psi_b^*({\mathbf{r}}-{\mathbf{n}}l )   \vphantom{\left[\frac{{\mathrm{e}}^{-\kappa r}}{r}\right]}\right\ }   \left\{r_\ell^m\left(-\frac1\kappa{\boldsymbol{\nabla}}_r\right )   \left[\frac{{\mathrm{e}}^{-\kappa r}}{r}\right]\right\ }   + { \mathcal{o}}\big({\mathrm{e}}^{-{\sqrt2\kappa l}}\big ) \,.\ ] ] we now integrate by parts and pass the laplacian through the differential operator @xmath156 . since the operators both consist of partial derivatives , this is not a problem when the wave function is smooth .
we assume that this is the case , with the possible exception of a measure zero region that can be omitted from the integral .
the partial integrations give a factor @xmath157 .
we can now proceed in exactly the same way as for s - waves .
we perform one more integration by parts so that the laplacian acts on @xmath132 .
this yields a delta function times a factor of @xmath158 , and the final result is then @xmath159 for @xmath160 we can insert the asymptotic form   since it is evaluated in the asymptotic region .      for @xmath149 , we find the same result for all three p - wave states , @xmath161 when compared to the s - wave case , the sign of the p - wave mass shift is reversed while the magnitude is the same .
qualitatively , this means that s - wave bound states are more deeply bound when put in a finite volume while p - wave bound states are less bound .
this will be discussed in more detail later .
we next discuss the results for @xmath162 . from
, we find @xmath163 where @xmath164    we note that the size and even the sign of the mass shift depends on the quantum number @xmath165 .
to understand this effect , we need to take into account that our cubic finite volume breaks the rotational symmetry group down to a cubic subgroup .
the cubic symmetry group @xmath166 is a finite subgroup of @xmath167 with 24 elements .
there are five irreducible representations of @xmath166 .
they are conventionally called @xmath168 , @xmath169 , @xmath170 , @xmath171 and @xmath172 , and their dimensionalities are @xmath173 , @xmath173 , @xmath174 , @xmath175 and @xmath175 , respectively .
irreducible representations @xmath176 of the rotation group @xmath167 are reducible with respect to @xmath166 for @xmath177 . for further details about the decomposition
see , for example , ref .
@xcite .
in our discussion we assume that the infinite volume system has no partial wave mixing , such that orbital angular momentum @xmath24 is a good quantum number .
we also assume that there are no accidental degeneracies in the bound state spectrum , so we can use @xmath24 as a label for the family of cubic representations split apart at finite volume .
parity invariance remains unbroken by the cubic volume , and we have @xmath178 just as in the infinite volume case
. for clarity , however , we will indicate parity explicitly with @xmath179 superscripts in the following .
with our assumptions , an s - wave state in infinite volume will map onto an @xmath180 state at finite volume .
also a p - wave triplet will map onto the three elements of the @xmath181 representation at finite volume . for d - waves , however , the five d - wave states are split into a @xmath182 triplet and an @xmath183 doublet , @xmath184    in the following we use the notation @xmath185 , @xmath186 , for the basis vectors of the irreducible cubic representations .
we can rewrite the finite volume mass shift in eq .
as @xmath187 where @xmath188 is the translation operator by displacement @xmath189 . we can also calculate the matrix elements of @xmath190 in the @xmath27 basis . in this case
there will be off - diagonal matrix elements connecting @xmath27 and @xmath191 when @xmath165 and @xmath192 are equivalent modulo 4 .    according to ref .
@xcite , the unitary transformation between the two basis sets for the five d - wave states is    @xmath193    and    @xmath194    so , for example , we have @xmath195 where we have defined @xmath196 as a straightforward generalization of  .
as expected from cubic symmetry , the mass shift is the same for all three @xmath182 states , and the same for both @xmath183 states . to summarize our results ,
we write the mass shift for a state belonging to irreducible representation @xmath197 with angular momentum @xmath24 as @xmath198 we list the coefficients @xmath199 for @xmath200 in table  [ tab : results ] .
.coefficient @xmath201 in the expression for the finite volume mass shifts for @xmath200 .
@xmath197 indicates the corresponding representation of the cubic group . [ cols="^,^,^",options="header " , ]      the sign of the finite volume mass shift can be understood in terms of the parity of the wave function . in infinite volume
the tail of each bound state wave function must vanish at infinity . in the finite volume , however , the bound state wave functions with even parity along a given axis can remain nonzero everywhere . only the derivative needs to vanish , and the kinetic energy
is lowered by broadening the wave function profile . on the other hand ,
a wave function with odd parity along a given axis must change sign across the boundary . in this case
the wave function profile is compressed and the kinetic energy is increased .
we have illustrated both cases for a one - dimensional square well potential in fig .
[ fig : wf ] .        in three dimensions ,
the situation is slightly more complicated , which can be seen from the fact that for @xmath162 the sign of the mass shift depend on the representation of the cubic group even though the parity is just @xmath202 for all states . in order to understand this
, we consider the basis polynomials for the cubic representations .
these basis polynomials are obtained by decomposing the cubic basis vectors in terms of solid harmonics which are homogeneous polynomials in @xmath203 , @xmath204 and @xmath205 . for @xmath206
the basis polynomials are also given explicitly in  @xcite .
for a given polynomial @xmath207 , we define its _ leading parity _ as @xmath208 where @xmath209 is the maximum degree of @xmath210 with respect to any one of the three variables .
it is this leading parity that determines the asymptotic behavior of the mass shift as @xmath211 .
more precisely , we have @xmath212 for the @xmath199 in eq .
[ eq : delta_m - gamma ] .    it can easily be checked that this relation holds for all results presented in table  [ tab : results ] . for @xmath162 , for example
, we have the basis polynomials    @xmath213    and hence @xmath214 for the @xmath182 representation and @xmath215 for the @xmath183 representation .
the expressions for the finite volume mass shift become simpler when we sum over all @xmath165 for a given @xmath24 .
we can rewrite as @xmath216 inserting the asymptotic form of the wave function , @xmath217^ *   \bigg|_{{\mathbf{r}}={\mathbf{n}}l } \,,\ ] ] and using a second time yields @xmath218   \bigg|_{{\mathbf{r}}={\mathbf{n}}l } \\   + { \mathcal{o}}\big({\mathrm{e}}^{-{\sqrt2\kappa l}}\big ) \,.\ ] ] now , from the well - known relation @xmath219 and @xmath220 we get an analogous expression for the solid harmonics , which then carries over to @xmath221 for any sufficiently smooth function @xmath144 .
finally , we have @xmath222 which follows from eq .  .
putting everything together , we arrive at @xmath223   \bigg|_{{\mathbf{r}}={\mathbf{n}}l } + { \mathcal{o}}\big({\mathrm{e}}^{-{\sqrt2\kappa l}}\big ) \\   & = ( -1)^{\ell+1}(2\ell+1)\cdot3|\gamma|^2\frac{{\mathrm{e}}^{-\kappa l}}{\mu l }   + { \mathcal{o}}\big({\mathrm{e}}^{-{\sqrt2\kappa l}}\big)\ , , \end{split } \label{eq : avshift-1}\ ] ] where the sum just yields a factor of six . dividing by @xmath224 , we obtain the average mass shift for states with angular momentum @xmath24 , @xmath225 apart from the alternating sign
, this average shift is independent of @xmath24 .
eq .   can be verified explicitly for the the results presented in sec .
[ sec : results ] ( _ cf .
_  table  [ tab : results ] ) . for @xmath162 , for example , one has to average over the three - dimensional representation @xmath182 and the two - dimensional representation @xmath183 .
in order to verify our predictions numerically , we put the schrdinger equation   on a discrete spatial lattice such that the hamiltonian becomes an ordinary matrix .
we then calculate the corresponding energy eigenvalues and eigenvectors .
we use a hat symbol to denote dimensionless lattice units .
for example , we have @xmath226 where @xmath227 denotes the lattice spacing .
the free lattice hamiltonian is given by @xmath228 \label{eq : h - latt}\ ] ] where @xmath229 and @xmath230 are creation and annihilation operators for a lattice site @xmath231 and @xmath232 is a unit vector in the @xmath233-direction .
the corresponding lattice dispersion relation is @xmath234 with the lattice function @xmath235\ ] ] and the lattice momenta @xmath236 the binding momentum for a bound state with energy @xmath237 is determined by @xmath238 the lattice green s function for the hamiltonian   is @xmath239    we impose periodic boundary conditions by defining the distance @xmath240 to the origin as @xmath241      we calculate the mass shift using three different methods :    1 .   as a direct difference in energies , eq .  , where we use a very large volume ( @xmath242 ) to approximate the infinite volume result
2 .   from the overlap formula .
3 .   using discretized versions of and  , which we obtain by replacing @xmath132 with the lattice green s function .
more precisely , we write the asymptotic bound state wave function   as @xmath243 and replace the continuum green s function @xmath244 with the lattice version @xmath245 effectively , this amounts to the replacement @xmath246 in the mass shift formula .
the lattice green s function is also used to calculate the asymptotic normalization @xmath7 from the lattice data .
this procedure has the advantage of avoiding large lattice discretization errors .      in the following we report physical quantities in units where the reduced mass @xmath1 is set to 1 .
we first use a gaussian potential , @xmath247 with @xmath248 and @xmath249 .
this potential does not have a finite range in a strict mathematical sense , but the range corrections can be entirely neglected in comparison with other errors in our numerical calculation . the smoothness of the gaussian potential helps to minimize lattice discretization artifacts . in fig .
[ fig : shifts - sp - gauss ] we show the s- and p - wave mass shifts obtained with the three methods described in sec .
[ sec : methods ] above .
the results from the three different methods described above agree well for both s- and p - waves . in order to compare the dependence on the box size @xmath4 with the predicted behavior , we have plotted @xmath250 against @xmath4 ( we use the absolute value of @xmath190 since the s - wave mass shift is negative ) . for both s- and p - waves ,
the expected linear dependence is clearly visible .     as functions of the box size @xmath4 ( in lattice units ) for a gaussian potential .
we show the results obtained from the direct difference eq .
( crosses ) , evaluation of the overlap integral eq .
( squares ) , and discretized versions of eqs .  , ( circles ) .
the dashed lines show linear fits to the overlap integral results.,scaledwidth=50.0% ]    when we perform a linear fit to the overlap integral data ( dashed lines in fig .
[ fig : shifts - sp - gauss ] ) we obtain @xmath251 , @xmath252 for the s - wave results and @xmath253 , @xmath254 for the p - wave results .
the values for the asymptotic normalization are in good agreement with the results @xmath255 ( s - wave ) and @xmath256 ( p - wave ) that are obtained directly from the @xmath257 data . inserting the corresponding energy eigenvalues into the lattice dispersion relation
, we find @xmath258 ( s - wave ) and @xmath259 ( p - wave ) , again in quite good agreement with the fit results .
the remaining small discrepancies can be attributed to the mixing with higher partial waves induced by the lattice discretization and the fact that we have not performed a continuum extrapolation to vanishing lattice spacing .      for a simple step potential , @xmath260 which we use with @xmath261 and @xmath262 ,
the numerical calculation becomes more difficult since the discontinuous shape introduces considerable lattice artifacts .
yet we discuss it here because it has a strict finite range and we find that for a small lattice spacing of @xmath263 the results are satisfactory . in fig .
[ fig : shifts - sp - step ] , we show a plot analogous to the one presented for the gaussian potential .
again , the results from the different methods agree well and the expected linear behavior is clearly visible . furthermore , the results from the three methods agree well with each other already for smaller @xmath4 ( compared to the results for the gaussian potential ) , as expected from the fact that the step potential does not have a tail .     as functions of the box size @xmath4 ( in lattice units ) for a simple step potential .
the symbols are as in fig .
[ fig : shifts - sp - gauss].,scaledwidth=50.0% ]    from fitting to the overlap integral data ( dashed lines in fig .
[ fig : shifts - sp - step ] ) we obtain @xmath264 , @xmath265 for the s - wave results and @xmath266 , @xmath267 for the p - wave results . from the @xmath257 data we find @xmath268 , @xmath269 ( s - wave ) and @xmath270 , @xmath271 ( p - wave ) .
given that we do not have error estimates for the @xmath257 results , the overall agreement is quite good .    finally , we also check our result for the d - wave mass splittings using the step potential with @xmath263 . in fig .
[ fig : shift - d - t2-e - step ] , we show the mass shift for the d - wave states in both the @xmath182 and the @xmath183 representation . due to the polynomial coefficients @xmath199 ( see eq .   and table  [ tab : results ] )
one does not expect a linear dependence on @xmath4 for @xmath250 .
hence , we simply plot @xmath190 as a function of @xmath4 directly and do not perform a fit .
nevertheless , we see that ( except for very small @xmath4 , where obviously the condition @xmath272 is not satisfied ) the agreement between the three methods to calculate @xmath190 is very good and hence conclude that our mass shift formula indeed gives the right result also for @xmath162 .     for @xmath182 rep .
( left panel ) and @xmath183 rep .
( right panel ) as a function of the box size @xmath4 ( in lattice units ) for a simple step potential .
the symbols are as in fig .
[ fig : shifts - sp - gauss].,title="fig:",scaledwidth=49.0% ]   for @xmath182 rep .
( left panel ) and @xmath183 rep .
( right panel ) as a function of the box size @xmath4 ( in lattice units ) for a simple step potential .
the symbols are as in fig .
[ fig : shifts - sp - gauss].,title="fig:",scaledwidth=49.0% ]
in this section we derive a formula for the finite - volume ( or rather finite - area ) mass shift of bound states in two - dimensional systems .
the results can be used , for example , in lattice simulations of cold atomic systems , which can be prepared experimentally to be effectively two - dimensional  @xcite .    in two dimensions ,
the schrdinger equation is @xmath273 with @xmath274\psi_b({\mathbf{r}})\ ] ] in polar coordinates .
states are described by a single angular momentum quantum number @xmath275 , and for the wave function we have the separation @xmath276 with @xmath277 the two linearly independent solutions of the free radial equation @xmath278 are just the bessel and neumann functions @xmath279 and @xmath280 . for a bound state , we have @xmath281 , and the wave function has the asymptotic form @xmath282 where @xmath283 is the modified bessel function of the second kind .
it is related to the hankel function of the first kind , @xmath284 via @xmath285 as in the three - dimensional case , @xmath7 is the asymptotic normalization . inserting   into   yields a form which is more similar to the three - dimensional expression , eq .  .
to make the analogy to the calculations in sec .
[ sec : p - wave - and - higher ] as explicit as possible , we will use the hankel function in the following intermediate steps and only express the final results in terms of the modified bessel function .
nearly all of the three - dimensional calculation carries over and we just replace all exponential terms by hankel functions .
the overlap integral for the mass shift is now @xmath286 from the asymptotic form of the hankel function , @xmath287 it is clear that in principle we still have an exponential behavior . in deriving eq .
we have used this to write @xmath288 in the following , we will simply write the correction terms as @xmath289 , as in the three - dimensional case .    the two - dimensional analogue of the relation   is @xmath290 where @xmath291 . this follows from @xmath292 and @xmath293 the proof for eq .   can be carried out in the same manner as the three - dimensional proof in lemma b.1 in  @xcite , where one uses the expansion of @xmath294 ( 2d vectors ) in terms of bessel functions . as the final ingredient
, we have @xmath295   \frac{{\mathrm{i}}}4 h_0^{(1)}({\mathrm{i}}\kappa r ) = -\delta^{(2)}({\mathbf{r } } ) \ , .
\label{eq : helmholtz - gf-2d}\ ] ] using all this in steps completely analogous to those in three dimensions , we get @xmath296 for @xmath150 ( two - dimensional s - waves ) , this directly yields @xmath297    in fact , eq .   can be simplified further .
inserting the asymptotic form for the wave function for @xmath298 and using   a second time gives @xmath299   \bigg|_{{\mathbf{r}}={\mathbf{n}}l } \\
+ { \mathcal{o}}\big({\mathrm{e}}^{-{\sqrt2\kappa l}}\big ) \,.\end{gathered}\ ] ] from   and @xmath300 it is clear that @xmath301which then yields @xmath302 for any sufficiently smooth @xmath144 .
this is essentially the same relation that we used to derive the trace formula in the three - dimensional case , only that here we do not have to sum over different @xmath165 .
together with the two - dimensional analogue of  , @xmath303 we then get @xmath304\bigg|_{{\mathbf{r}}={\mathbf{n}}l }   + { \mathcal{o}}\big({\mathrm{e}}^{-{\sqrt2\kappa l}}\big ) \\   & = ( -1)^{m+1}\cdot\frac{2|\gamma|^2}{\mu } k_0(\kappa l )   + { \mathcal{o}}\big({\mathrm{e}}^{-{\sqrt2\kappa l}}\big ) \ , .
\end{split } \label{eq : delta_m - general-2d-2}\ ] ] as we shall see in the following , this is the final result for @xmath150 and any odd @xmath165 , whereas for even @xmath305 things become slightly more complicated .    in general , we have to take into account that the finite volume breaks the original planar rotational symmetry of the system down to the symmetry group of a square .
we find that states with the same absolute value of @xmath165 may mix to form good eigenstates in the finite volume .
more precisely , we have the symmetric and antisymmetric combinations @xmath306 for @xmath305 .
when we calculate the mass shift for these states ( in the same way as described in sec .
[ sec : results ] ) , we get mixing terms of the form @xmath307    since the condition for the mixing of states is @xmath308 they do not play a role for odd @xmath165 ( in fact , they vanish in this case ) . for even @xmath165
, however , we have to take them into account and find @xmath309 as our final result . as an illustration
, we give the explicit results for @xmath310 :    @xmath311   + { \mathcal{o}}\big({\mathrm{e}}^{-{\sqrt2\kappa l}}\big ) \,,\ ] ]    @xmath312   + { \mathcal{o}}\big({\mathrm{e}}^{-{\sqrt2\kappa l}}\big ) \,.\ ] ]
in this work we have derived explicit formulae for the mass shift of p- and higher - wave bound states in a finite volume and discussed their decomposition into states transforming according to the representations of the cubic group .
we have compared our numerical results for @xmath313 with numerical calculations of the finite - volume dependence for lattice gaussian and step potentials and found good agreement with the predictions . for @xmath314 , the mass shift of a given
state @xmath27 depends on the angular momentum projection @xmath165 due to the breaking of rotational symmetry .
when we average over all @xmath165 in a multiplet , however , the absolute value of the mass shift is even independent of @xmath24 . the mass shift for states in representations of the cubic group , however , is the same for all states .
the sign of the mass shift can be understood from the leading parity of the representations .
finally , we have derived corresponding expressions for the finite volume mass shift in two - dimensional systems .
lattice calculations also provide a method to extract asymptotic normalization coefficients , which are of interest in low - energy astrophysical capture reactions .
we have pointed out how the asymptotic normalization and binding momentum of a shallow bound state can be used to extract the effective range from a simulation .
our work provides a general framework for future lattice studies of molecular states with angular momentum in systems with short - range interactions .
applications to nuclear halo systems and molecular states in atomic and hadronic physics appear promising .
an important next step would be to include coulomb effects into the framework in order to investigate the volume dependence of bound states of charged particles , which are much easier to treat experimentally .
in particular , this extension of the formalism is important if one wants to describe proton - halo nuclei .
another interesting direction would be to analyze the volume dependence of resonances along the lines of refs .
@xcite .
this research was supported in part by the dfg through sfb / tr 16 `` subnuclear structure of matter '' , the bmbf under contract no . 06bn9006 , and by the us department of energy under contract no .
de - fg02 - 03er41260 .
s.k . was supported by the `` studienstiftung des deutschen volkes '' and by the bonn - cologne graduate school of physics and astronomy .
in this section , we derive the relation between the asymptotic normalization of the bound state wave function and the effective range of the corresponding two - particle scattering process that we have quoted in sec .  [
sec : prereq ] .
it was already pointed out by lscher in  @xcite that the asymptotic normalization @xmath7 of the bound state wave function is related to scattering parameters .
more precisely , the analytically - continued elastic scattering amplitude in forward direction has a pole at the bound state energy , and the residue of this pole is proportional to @xmath315 .    in the limit of shallow bound states , _
i.e. _ , @xmath316 , it is possible to make a more direct connection to the effective range in the corresponding scattering channel . the crucial ingredients for this are given in a general investigation of wigner causality bounds  @xcite .
we express the elastic scattering amplitude @xmath317 in terms of the scattering phase shift @xmath318 , @xmath319 } \,,\ ] ] for which we have the well - known effective range expansion @xmath320 where @xmath321 and @xmath322 are the scattering and effective range parameters , respectively .
a scattering state is described by a wave function @xmath323 that is a solution of the radial schrdinger equation   for positive real center - of - mass momentum @xmath324 . as it is done in  @xcite
, we choose the normalization such that outside the range of the potential ( @xmath9 ) we have @xmath325 \,.\ ] ] for a bound state , we have @xmath47 and @xmath326 .
hence , the wave function is @xmath327 \\
= { \mathrm{i}}^\ell \kappa^\ell \hat{h}_\ell^+({\mathrm{i}}\kappa r)\ ] ] in this case .
we immediately see that @xmath328 where @xmath29 is the radial bound state wave function from the previous sections ( normalized to @xmath173 ) .
note that the phase convention is chosen such that the wave functions are real . from  @xcite
we have the relation @xmath329 ^ 2 \label{eq : er-1}\ ] ] with @xmath330 for @xmath52 , we get @xmath331 note that @xmath332 as @xmath316 , and since we only consider finite - range potentials , we can in fact write @xmath333 .    using the normalization of the wave function we can rewrite   as @xmath334 ^ 2\right\ } \ , .
\label{eq : er-2}\ ] ] the remaining integral can be expressed as @xmath335 ^ 2   = \kappa^{2\ell}\int\limits_r^\infty\mathrm{d}r\ ,   \left[{\mathrm{i}}^\ell\hat{h}_\ell^+({\mathrm{i}}\kappa r ) \right]^2 \ , .
\label{eq : int - exterior - intermediate}\ ] ] for @xmath52 , we get @xmath335 ^ 2   = \frac{{\mathrm{e}}^{-2\kappa r}}{2\kappa }   = \frac1{2\kappa } - r + { \mathcal{o}}(\kappa ) { { \ \ } \text{as}{\ \ } } \kappa\to0 \,.\ ] ] together with   this yields @xmath336 up to the given order this is equivalent to the relation @xmath337 from  @xcite . for @xmath54 ,
the integral is @xmath335 ^ 2   = \frac{\gamma\left(\ell-\frac12\right)\gamma\left(\ell+\frac12\right)}\pi   \left(\frac r2\right)^{-2\ell+1 } + { \mathcal{o}}(\kappa )   { { \ \ } \text{as}{\ \ } } \kappa\to0 \ , .
\label{eq : int - exterior}\ ] ] we see that the leading term exactly cancels the one in   such that we arrive at @xmath338    in all expressions above , the asymptotic normalization of course depends on the binding momentum , _ i.e. _ , @xmath339 , only in general this dependence can not be calculated analytically .
using the results that we have derived in this paper , it is possible to determine both @xmath12 and @xmath7 from the volume dependence of a given bound state .
if this state is sufficiently shallow , the formulae derived in this chapter can then be used to get an estimate for the effective range of the interaction .
d.  lee , prog .  part .  nucl .
* 63 * ( 2009 ) 117 [ arxiv:0804.3501 [ nucl - th ] ] .
a.  bazavov _ et al .
_ , rev .
* 82 * ( 2010 ) 1349 [ arxiv:0903.3598 [ hep - lat ] ] . s.  r.  beane , w.  detmold , k.  orginos and m.  j.  savage , prog .
* 66 * ( 2011 ) 1 [ arxiv:1004.2935 [ hep - lat ] ] . m.  lscher , commun .
* 104 * ( 1986 ) 177 .
s.  bour , s.  knig , d.  lee , h .- w .
hammer , and u .- g .
meiner , [ arxiv:1107.1272 [ nucl - th ] ] .
c.  a.  regal , c.  ticknor , j.  l.  bohn and d.  s.  jin , phys .
* 90 * ( 2003 ) 053201 [ arxiv : cond - mat/0209071 ] . c.  h.  schunck _ et al . _ , phys .
a * 71 * ( 2005 ) 045601 [ arxiv : cond - mat/0407373 [ cond-mat.soft ] ] .
s.  typel and g.  baur , phys .
* 93 * ( 2004 ) 142502 [ arxiv : nucl - th/0406068 ] .
w .  hammer and d.
r.  phillips , nucl .
phys .  a  * 865 * ( 2011 ) 17 [ arxiv:1103.1087 [ nucl - th ] ] .
g.  rupak and r.  higa , phys .
* 106 * ( 2011 ) 222501 [ arxiv:1101.0207 [ nucl - th ] ] .
a.  tohsaki , h.  horiuchi , p.  schuck and g.  rpke , phys .  rev .  lett .
* 87 * ( 2001 ) 192501 [ arxiv : nucl - th/0110014 ] .
j.  r.  taylor , _ scattering theory : the quantum theory of nonrelativistic collisions _ , dover publications , inc .
m.  l.  goldberger and k.  m.  watson , _ collision theory _ , john wiley & sons ,