On the modification of the Efimov spectrum in a finite cubic box
Abstract
Three particles with large scattering length display a universal spectrum of threebody bound states called “Efimov trimers”. We calculate the modification of the Efimov trimers of three identical bosons in a finite cubic box and compute the dependence of their energies on the box size using effective field theory. Previous calculations for positive scattering length that were perturbative in the finite volume energy shift are extended to arbitrarily large shifts and negative scattering lengths. The renormalization of the effective field theory in the finite volume is explicitly verified. We investigate the effects of partial wave mixing and study the behavior of shallow trimers near the dimer energy. Moreover, we provide numerical evidence for universal scaling of the finite volume corrections.
I Introduction
Fewbody systems with resonant interactions characterized by a large scattering length show interesting universal properties. If is positive, two particles of mass form a shallow dimer with energy , independent of the mechanism responsible for the large scattering length. Examples for such shallow dimer states are the deuteron in nuclear physics, the He dimer in atomic physics, and possibly the new charmonium state in particle physics Braaten:2004rn (); Platter:2009gz (). In the threebody system, the universal properties include the Efimov effect Efimov70 (). If at least two of the three pairs of particles have a large scattering length compared to the range of their interaction, there is a sequence of threebody bound states whose energies are spaced geometrically between and . In the limit , there are infinitely many bound states with an accumulation point at the threebody scattering threshold. These Efimov states or trimers have a geometric spectrum Efimov70 ():
(1) 
which is specified by the binding momentum of the Efimov trimer labeled by . This spectrum is a consequence of a discrete scaling symmetry with discrete scaling factor . In the case of identical bosons, and the discrete scaling factor is . The discrete scale invariance persists if is large but finite, but in this case it connects states corresponding to different values of the scattering length. The scaling symmetry becomes also manifest in the logperiodic dependence of scattering observables on the scattering length Efimov79 (). The consequences of discrete scale invariance and “Efimov physics” can be calculated in an effective field theory for shortrange interactions, where the Efimov effect appears as a consequence of a renormalization group limit cycle Bedaque:1998kg ().
While the Efimov effect was established theoretically already in 1970, first experimental evidence for an Efimov trimer in ultracold Cs atoms was provided only recently by its signature in the threebody recombination rate Kraemer06 (). It could be unravelled by varying the scattering length over several orders of magnitude using a Feshbach resonance. Since this pioneering experiment, there was significant experimental progress in observing Efimov physics in ultracold quantum gases. More recently, evidence for Efimov trimers was also obtained in atomdimer scattering Knoop08 () and in threebody recombination in a balanced mixture of atoms in three different hyperfine states of Li Ottenstein08 (); Huckans08 (), in a mixture of Potassium and Rubidium atoms Barontini09 (), and in an ultracold gas of Li atoms Gross09 (). In another experiment with Potassium atoms Zaccanti09 (), two bound trimers were observed.
The observation of Efimov physics in nuclear and particle physics systems is complicated by the inability to vary the scattering length and one has to focus on the detection of excited states. While twoneutron halo nuclei could be bound due to the Efimov effect, the analysis of known halo nuclei does not show much promise for an unambiguous identification (See Ref. Canham:2008jd () and references therein). Another opportunity to observe Efimov physics is given by lattice QCD simulations of threenucleon systems Wilson:2004de (). A number of studies of the quarkmass dependence of the chiral nucleonnucleon () interaction found that the inverse scattering lengths in the relevant – and channels may both vanish if one extrapolates away from the physical values to slightly larger quark masses Beane:2001bc (); Beane:2002xf (); Epelbaum:2002gb (). This implies that QCD is close to the critical trajectory for an infrared renormalization group limit cycle in the threenucleon sector. It was conjectured that QCD could be tuned to lie precisely on the critical trajectory by tuning the up and down quark masses separately Braaten:2003eu (). As a consequence, the triton would display the Efimov effect. More refined studies of the signature of Efimov physics in this case followed Epelbaum:2006jc (); Hammer:2007kq (). However, a proof of this conjecture can only be given by an observation of this effect in a lattice QCD simulation Wilson:2004de (). The first full lattice QCD calculation of nucleonnucleon scattering was reported in Beane:2006mx () but statistical noise presented a serious challenge. A promising recent highstatistics study of threebaryon systems presented also initial results for a system with the quantum numbers of the triton such that lattice QCD calculations of threenucleon systems are now within sight Beane:2009gs (). For a review of these activities, see Ref. Beane:2008dv (). Since lattice simulations are carried out in a cubic box, it is important to understand the properties of Efimov states in the box. The first step towards this goal is to understand these modifications for a system of three identical bosons.
The corresponding modifications of the Efimov spectrum can be calculated in effective field theory (EFT) since the finite volume modifies the infrared properties of the system. The properties of twobody systems with large scattering length in a cubic box were calculated in Ref. Beane:2003da (). Some properties of threebody systems in a finite volume have also been studied previously. For repulsive and weakly attractive interactions without bound states, Tan has determined the volume dependence of the ground state energy of three bosons up to Tan08 (). In Refs. Beane:2007qr (); Detmold:2008gh (), this result was extended for general systems of bosons. It was used to analyze recent results for three and more boson systems from lattice QCD Beane:2007es (); Detmold:2008fn (); Detmold:2008yn (). For the unitary limit of infinite scattering length, some studies have been carried out as well. The properties of three spin1/2 fermions in a box were investigated in Pricoup07 (). However, this system has no threebody bound states in the infinite volume. The volume dependence of energy levels in a finite volume can also be used to extract scattering phase shifts and resonance properties from lattice calculations Luscher:1990ux (); Luscher:1991cf (). For a recent application of this idea to the resonance, see Refs. UGM1 (); UGM2 ().
In a previous letter, we have investigated the finite volume corrections for a threeboson system with large but finite scattering length Kreuzer:2008bi (). We have studied the modification of the bound state spectrum for positive scattering length using an expansion valid for small finite volume shifts and explicitly verified the renormalization in the finite box within this approximation. In the current paper, we present an extension of this work that avoids the expansion in the energy shift and can be applied to arbitrarily large shifts. Moreover, we extend our previous studies to the case of negative scattering lengths. Our results indicate that the finite volume corrections are subject to a universal scaling relation. We also provide a more detailed discussion of the technical details and our numerical method.
Ii Theoretical Framework
ii.1 Lagrangian and BosonDiboson Amplitude
Three identical bosons interacting via shortrange forces can be described by the Lagrangian Bedaque:1998kg (); Braaten:2004rn ()
(2) 
where the dots indicate higher order terms. The degrees of freedom in the Lagrangian are a boson field and a nondynamical auxiliary field representing the composite of two bosons. Units have been chosen such that . This Lagrangian corresponds to the zero range limit with . Our results will therefore be applicable to physical states whose size is large compared to . Corrections from finite range can be included via higher order terms in the Lagrangian but will not be considered here. The coupling constants and will be determined by matching calculated observables to a twobody datum and a threebody datum, respectively. The central quantity in the threebody sector of this EFT is the bosondiboson scattering amplitude. It determines all threebody observables and is given as the solution of an inhomogeneous integral equation, depicted diagrammatically in Fig. 1. A detailed discussion of the infinite volume case can be found in Braaten:2004rn (). Here we follow the strategy of Refs. Kreuzer:2008bi (); Beane:2003da () and focus on the finite volume case.
The system of three bosons is assumed to be contained in a cubic box with side length and periodic boundary conditions. This leads to quantized momenta . In particular, the possible loop momenta are quantized. As a consequence, the integration over loop momenta is replaced by an infinite sum. The divergent loop sums are regulated by a momentum cutoff similar to the infinite volume case.
The finite volume modifies the infrared physics of the system but does not change the ultraviolet behavior of the amplitudes. Therefore, the renormalization is the same in the infinite and finite volume cases. Of course, this statement only holds if the momentum cutoff is large compared to the momentum scale set by the size of the volume, namely , such that the infrared and the ultraviolet regime are well separated. In the numerical calculations presented in this work, the consistent renormalization of the results will always be verified explicitly. The twobody sector can be completely renormalized by matching the twobody coupling constant to a lowenergy twobody observable, namely the twobody scattering length or the dimer binding energy. Because of the discrete scaling symmetry, the threebody coupling approaches an ultraviolet limit cycle. For convenience, the threebody coupling constant is expressed in the form . The cutoff dependence of the dimensionless function is then given by Braaten:2004rn ()
(3) 
where and is a threebody parameter that can be fixed from a trimer binding energy or another threebody datum.
For the diboson lines in the threebody equation depicted in Fig. 1, the full, interacting diboson propagator has to be used. This quantity corresponds to the exact twobody scattering amplitude. It is obtained by dressing the bare diboson propagator, given by the constant , with bosonic loops as shown in Fig. 2. This leads to an infinite sum that can be evaluated analytically. For a diboson with energy , the propagator is
(4) 
In the limit , this expression reduces to the full diboson propagator in the infinite volume case.
Having obtained the full diboson propagator, the integral equation for the bosondiboson scattering amplitude can be written down explicitly. Using the diagrammatical representation in Fig. 1 and the Feynman rules derived from the effective Lagrangian (2), we obtain:
Here, () are the momenta of the incoming (outgoing) dibosons, while the momenta of the incoming (outgoing) bosons are (). The boson legs have been put on shell but the diboson legs remain offshell. The total energy of the system, , can be treated as a parameter of the equation. The integration over the loop energy can be performed by virtue of the residue theorem. This yields the sum equation
(6) 
with
(7)  
(8) 
If the energy is near a trimer energy , the amplitude exhibits a simple pole and the dependence on the momenta seperates:
(9) 
Matching the residues on both sides of Eq. (6), the boundstate equation
(10) 
is obtained. Values of the energy , for which this homogeneous sum equation has a solution, are identified with the energies of the trimer states.
ii.2 Cubic symmetry
In the infinite volume case, only swave bound states are formed. However, in a finite cubic volume, the spherical symmetry of the infinite volume is broken to a cubic symmetry. In the language of group theory, the infinitely many irreducible representations of the spherical symmetry group are mapped onto the five irreducible representations of the cubic group . The representations of the spherical symmetry group are now reducible and can therefore be decomposed in terms of the five irreducible representations of the cubic group. On the other hand, a quantity transforming according to the irreducible representation of can be written in terms of the basis functions of the spherical symmetry, i.e. the spherical harmonics , via
(11) 
where and is an additional index needed if the representation labeled by appears in the irreducible representation more than once. The linear combinations of spherical harmonics are called “kubic harmonics” BetheVdLage:47 (). The values of the coefficients are known for values of as large as 12 Altmann:65 ().
In order to make contact with the infinite volume formalism, Eq. (10) is rewritten using Poisson’s resummation formula. This identity states that an infinite sum over integer vectors may be replaced by an infinite sum over integer vectors of the Fourier transform of the addend, i.e.
(12) 
where is the Fourier transform of . Applying the identity (12) to Eq. (10) yields
(13) 
The term with gives the corresponding equation in the infinite volume, while the other terms of the infinite sum may be viewed as corrections due to momentum quantization and the breakdown of the spherical symmetry. The explicit recovery of the infinite volume term is useful for bound states, since in this case the analytic structure of the amplitude, i.e. the pole at the binding energy, is identical and only the position of the pole is changed. This approach may be inappropriate for states belonging to the continuous scattering region of the infinite volume.
In order to evaluate the angular integration in Eq. (13), all quantities with angular dependence are expanded in terms of the basis functions of the irreducible representations of , namely in spherical harmonics. The amplitude itself is assumed to transform under the trivial representation of the cubic group, since the representation is solely contained in . The amplitude can therefore be written as
(14) 
The sum runs over those values of associated with the representation of . The first few values are . Since the first occasion where an value appears more than once is , the summation over the multiplicity and hence the index will be suppressed in the following.
The only angular dependence of the quantity is on the cosine of the angle between and . Therefore, can be expanded in Legendre polynomials . These polynomials can in turn be expressed in spherical harmonics via the addition theorem, yielding
(15) 
The exponential function in Eq. (13) can be rewritten using the identity
(16) 
where is the spherical Bessel function of order .
With these expansions, the angular integration in Eq. (13) can be performed. The integral over the three spherical harmonics depending on yields Wigner 3 symbols. Projecting on the th partial wave results in an infinite set of coupled integral equations for the quantities :
(17) 
The sum runs over the partial waves associated with the representation. The quantity can be calculated from Eq. (15) to be
(18) 
Here, is a Legendre function of the second kind. The threebody contact interaction contributes only to the swave, as expected.
Since the bound states in the infinite volume are swave states, Eq. (17) is specialized to the case , yielding
(19) 
The specialization of Eq. (18) to the case reads
(20) 
101  606  96.23  81.33 

1001  9126  94.05  20.69 
3501  39678  118.67  15.73 
The second line of Eq. (19) indicates admixtures from higher partial waves. Since the leading term in the expansion of the spherical Bessel function is , these contributions are suppressed by at least . They will therefore be small for volumes not too small compared to the size of the bound state. Moreover, contributions from higher partial waves will be suppressed kinematically for shallow states with small binding momentum. This is ensured by the spherical harmonic in the second line of Eq. (19). For higher partial waves, the angular sum yields smaller prefactors relative to the swave () for a given absolute value (see Table 1). Only for small lattices, i.e. when is large, this behavior is counteracted by terms stemming from the spherical Bessel function and higher partial waves may contribute significantly. The first numerical calculations presented in this work were therefore mainly performed neglecting contributions from higher partial waves. For two specific examples, however, we explicitly calculate the contribution from the next partial wave and demonstrate that it is small. Some details on the numerical solution of Eq. (19) are given in the Appendix. In the following section, we present our results.
Iii Results and discussion
By employing the formalism laid out in the previous section, we have calculated energy levels in finite cubic volumes of varying side lengths. The results of these calculations are presented in the following. For convenience, the dependence of the energies on the boson mass is reinstated in this section.
iii.1 Positive scattering length
We first present results for systems with . In this regime, a physical diboson state with a binding energy exists. This energy is therefore identical with the threshold for the breakup of a trimer into a diboson and a single boson. We choose states with different energies in the infinite volume, including shallow as well as deeply bound states:

, ,

, ,

, ,

, ,

, ,
Here, is the trimer energy in the infinite volume. Note that the states Ia, Ib and Ic appear in the same physical system chracterized by . For each of these states, its energy in the finite cubic volume has been calculated for various values of the box side length . In order to check the consistency of our results, the calculation was carried out for several cutoff momenta . For each cutoff, the threebody interaction parameterized by has been adjusted such that the infinite volume binding energies are identical for all considered values of . If our results are properly renormalized, the results for the different cutoffs should agree with each other up to an uncertainty of order stemming from the finiteness of the cutoff.
The results for the states II and III are depicted in Fig. 3 for box sizes between and . The values obtained for different cutoffs indeed agree with each other within the depicted uncertainty bands. Note that these bands do not represent corrections from higher orders of the EFT. For both states, the infinite volume limit is smoothly approached. As the volume becomes smaller, the energy of the states is more and more diminished. This corresponds to an increased binding with decreasing box size.
In the infinite volume, state III is more deeply bound than state II. Naïvely, one therefore expects the former to have a smaller spatial extent than the latter. The size of the state can be estimated via the formula , yielding for state II and for state III. Hence, a given finite volume should affect state II more strongly than the smaller state III. This behavior can indeed be observed. For example, considering a cubic volume with side length , the energy of state II deviates 10% from the infinite volume value, whereas the corresponding difference for state III is less than one percent. On the other hand, the box size for which the energy shift of the state III amounts to 10% is roughly .
For both states, we can now form the dimensionless number , where is the box size at which the energy differs by 10% from the infinite volume value. This yields for state II and for state III. The approximate equality of the two values of may indicate the presence of universal scaling in the finite volume version of the Effective Theory.
In Fig. 4, the three datasets obtained for state II are plotted again, this time in comparison to the results from a calculation using the expansion of the kernel used in Ref. Kreuzer:2008bi () and described in Appendix A. For large values of the box size , the results of both calculations agree with each other. For volumes smaller than , the result of the calculation with the expansion deviates from the result of the full calculation. For this size of the volume, the full result differs by about 20% from the infinite volume binding energy. Accordingly, the expansion employed for the integral kernel can not be applicable any longer. In Fig. 4, results are only depicted for volumes with . For smaller volumes, the results of calculations using the expansion become cutoff dependent, and are hence not are not properly renormalized anymore.
In Fig. 5, we show our results for the two states Ia and Ib in the same physical system characterized by . The volume dependence of the two states is again shown for box sizes between and . The curve corresponding to the more deeply bound state Ib shows a behavior similar to the one observed for the states II and III. The binding energy remains constant until the volume is small enough to affect the state. At this point, the energy of the state is more and more diminished as the volume becomes smaller.
The behavior of the shallow state Ia is different. In the region , the binding is not further increased. The results for smaller volumes show a sharp rise of the threebody energy. The energy of the state becomes positive near . State Ia is close to the threshold for bosondiboson scattering in the infinite volume located at . For comparison, we calculated the energy of the physical diboson according to Beane:2003da (). The resulting curve is the solid line in Fig. 5. Like the energy of the threebody bound states, the energy of the diboson is diminished in finite volumes. For volumes of the size , the energy of the threebody state Ia becomes larger than the diboson energy and starts to grow. This behavior is consistent with the observation that states are always shifted away from the threshold in a finite volume. In the twobody sector, for example, continuum states have been shown to have a power law dependence on the volume, while the volume dependence of bound states is dominated by exponentials Beane:2003da (). The data shown for state Ia can be explained by an exponential for , which characterizes the state as a bound state. For , the data is consistent with a power law, indicating the state indeed behaves like a bosondiboson scattering state if its energy is above the diboson energy.
The other investigated states do not show such a transition since their energy is well below the diboson energy for all considered volumes. It is unclear whether other states would show a similar behavior for smaller box sizes. If this is not the case, the transition from bound to unbound would occur only for states with infinite volume binding energies up to a critical value. For a conclusive analysis of the nature of the described transition, more data is still needed. It would be interesting to see whether such a transition shows up in lattice data for a state that is very close to the diboson threshold in the infinite volume.
,  ,  

1.25  200  4.30392  4.24545  3.53099 
200  4.57097  4.57753  
1  300  4.60298  4.60872  5.07581 
400  4.59579  4.60056  
0.9  200  4.31223  4.27927  6.068 
Since the size of the finite volumes where the state crosses the diboson energy is comparable to the size of the state itself, the breaking of the spherical symmetry may be a relevant effect here. To assess the influence of the higher partial waves, we extract from Eq. (17) two coupled equations for the swave amplitude () and the amplitude . These coupled equations are then solved in a coupled channel approach. The energies of state Ia obtained by this method are summarized and compared to the swave only results in Table 2. For , the state is still below the dimer state. The inclusion of the higher partial wave leads to a small downward shift in the energy. For , we have done calculations using 3 different cutoffs. The results for different cutoffs agree to two significant digits indicating the results are properly renormalized, but the binding is slightly reduced by the contribution. For , the effect of the higher partial wave is again a small downward shift. All results show only a deviation of about 1% from the swave only result. In summary, we find that the correction from the admixture is extremely small. Moreover, the correction is of the same order of magnitude as the finite cutoff corrections in our calculation and a more quantitative statement requires improved numerical methods.
,  ,  

1  200  11.86  11.15 
400  11.79  11.08  
0.7  200  19.06  20.70 
400  18.97  20.64 
In order to establish an estimate of typical corrections from higher partial waves, we have also performed calculations including the contributions for the more deeplybound state II. The resulting numbers are given in Tab. 3. The investigated box sizes are about three times larger than the state itself. The contribution of the higher partial wave is now several percent. This is still a small correction but considerably larger than the finite cutoff uncertainty. This suggests that the extremely small corrections for state Ia are related to its unusual behavior. The dominance of the swave my be associated with the closeness of the state to the threshold. A more detailed investigation of higher partial waves including a description of the numerical methods will be the subject of a future publication.
Since the shallower state Ia is more affected by a given finite volume than the deeper state Ib, the ratio of the energies of the states is changing. In the infinite volume, this ratio is 23.08. For , just before the shallow state crosses the dimer energy, this ratio has decreased to 7.4. Note that this ratio differs from the discrete scaling factor even in the infinite volume limit. This behavior is expected for shallow states close to the bound state threshold Braaten:2004rn (). The ratio will be approached when deeper states are considered. For example, the infinite volume ratio of state Ib and the much more deeply bound state Ic is 344 and already closer to the discrete scaling factor 515.
If we assume that the combination is indeed a universal number, we are able to predict for the states Ia and Ib. The results are for state Ib and for state Ia. The energies calculated for these volumes are , corresponding to an 8% shift, for state Ib and , also corresponding to an 8% shift. These findings support the assumption that the finite volume corrections obey universal scaling relations.
The states Ia and Ib appear in the same physical system characterized by . In this system, an even more deeply bound state, denoted as state Ic, is present. In the infinite volume, the energy of this state is . This corresponds to a binding momentum of . Since this is already comparable to the momentum cutoffs of a few hundred inverse scattering lengths employed before, we also performed calculations for this state using a much larger cutoff of . The threebody force for this cutoff has been fixed such that the energy of the most shallow state Ia is reproduced. The resulting energy of state Ic is then . This differs from the energy obtained using the smaller cutoff by 0.4%. This difference can be attributed to effects stemming from the finiteness of the cutoff. We have calculated the energy of this state in finite volumes of the sizes , and . The results of this calculations are summarized in Table 4. The values obtained using the large cutoff show no effect of the finite volume at all. From the infinite volume binding energy, the size of the state can be estimated via to be . So, we do not expect any visible effect since the finite volume is fifty times as large as the state itself. However, for the smaller cutoff , there are very small deviations from the infinite volume energy. But these deviations are smaller than the uncertainty stemming from the finiteness of the cutoff, which is estimated to be of order .
9440.91  –  9401.32  –  
1  9440.91  0%  9401.36  
0.75  9440.91  0%  9400.53  
0.5  9440.91  0%  9399.38  0.02% 
iii.2 Negative scattering length
Now we turn to systems with . In this regime, the twobody interaction is attractive but no diboson bound state exists in the infinite volume limit. The only possible breakup process for a threeboson bound state is therefore the breakup into three single bosons. The threshold for this process is . As for the case with positive scattering length, we choose states with different energies in the infinite volume:

, ,

, ,

, ,
Here, is the trimer energy in the infinite volume. For each of these states, its energy in the finite cubic volume has been calculated for various values of the box side length . In order to check the consistency of our results, the calculation was carried out for several cutoff momenta . For each cutoff, the threebody interaction parameterized by has been adjusted such that the infinite volume binding energies are identical for all considered values of .
The results for the states NI and NII are shown in Fig. 6 for box sizes between and . The results for the state NIII with box sizes between and are depicted in Fig. 7. The values obtained with different cutoffs all agree with each other within the uncertainty bands indicating proper renormalization. Our findings are similar to those in the positive scattering length regime. All three states smoothly approach the infinite volume limit. As the box size becomes smaller the energy of the state is more and more diminished. The overall behavior is identical to the one in the positive scattering length case described in the previous section.
State  

NI  29%  1.2  0.7  
NII  118%  1.8  1.05  
NIII  3250%  7.2  4.2 
The more deeply bound a state is in the infinite volume, the smaller is its spatial extent. Estimating the size via the formula yields for state NI, for state NII, and for state NIII. A given finite volume should therefore affect state NII more than state NI, and state NIII should be the most affected. When considering a cubic volume with side length , we find the energies given in Table 5. The relative deviation of state NI is indeed four times smaller than the shift for state NII and a hundred times smaller than the shift for the shallow state NIII. On the other hand, the box length , for which the energy of each states deviates 10% from its infinite volume value, is smaller the more deeply bound a state is. These values are also given in Table 5.
As in the case of positive scattering length, we have investigated the dimensionless combination . From state NI, we obtain . This leads to the predictions for state NII and for state NIII in good agreement with the explicitly calculated values. Additionally, we form another dimensionless combination, , where is the box side length where the energy of the state is twice as large as the infinite volume value. For the three investigated states, this box size is also given in Table 5. From the value for the state NI, we get . From this value, we predict for state NII and for state NIII. Again, the value for state NII is as predicted, whereas the value for state NIII is about 10% off the prediction from universal scaling.
For states far away from the threshold, we expect that the regimes of negative and positive scattering length are governed by the same scaling factor. Therefore, the dimensionless combination should, for such deeply bound states, have a common value for both signs of . For positive scattering lengths, state Ib is an example of a rather deeply bound state. For this state, we find and from that . For negative scattering lengths, we chose a state with by setting . The energy of this state is shifted by 10% in a volume with side length , yielding . The values of for the two signs of the scattering length are indeed close to each other for these two states. This behavior provides numerical evidence for the universality of finite volume effects for both positive and negative scattering lengths.
Iv Conclusions
In this paper, we have extended our earlier studies of the Efimov spectrum in a cubic box with periodic boundary conditions Kreuzer:2008bi (). The knowledge of the finite volume modifications of the spectrum is important in order to understand results from future 3body lattice calculations. Using the framework of EFT, we have derived a general set of coupled sum equations for the Efimov spectrum in a finite volume.
Specializing to , we have calculated the spectrum for both positive and negative scattering lengths and verified the renormalization in the finite box explicitly. We have removed the expansion for small finite volume shifts used in Ref. Kreuzer:2008bi () and presented an extension that can be applied to arbitrarily large shifts. Typically, the binding of all threebody states increases as the box size is reduced. Moreover, we provided a more detailed discussion of the technical details and our numerical method. We have investigated the breakdown of the linear approximation in detail and find that the expansion is applicable as long as the finite volume shift in the energy is not larger than 15–20% of the infinite volume energy.
We have studied the spectrum for positive and negative scattering lengths in detail and provided numerical evidence for universal scaling of the finite volume effects. The scaling properties can be quantified by the dimensionless product of the threebody binding energy in the infinite volume limit and the square of the box length corresponding to a finite volume shift of ten percent of the infinite volume energy. For sufficiently deep states, we obtained numerical evidence that approaches a universal number for both signs of the scattering length. These findings suggest that the properties of deeper states in the finite volume can be obtained from a simple rescaling and do not require explicit calculations. A more detailed analysis of this issue along the lines of Braaten:2002sr () is in progress.
For positive scattering lengths, we have investigated a spectrum of three states in the same physical system and studied the behavior of the shallowest state in the vicinity of the dimer energy which specifies the scattering threshold in the infinite volume. We found that this state drastically changes its behavior as a function of the box length once its energy becomes equal to the dimer energy. At this point the energy of the shallowest threebody state starts to grow and eventually becomes positive. The observed behavior is consistent with exponential suppression of the finite volume corrections below the dimer energy and power law suppression above.
The effect of an admixture of the partial wave has been investigated for two different threebody states for box sizes of the order of the size of the state. For a generic state well separated from threshold, we found corrections of the order of a few percent for volumes that are about three times larger than the state itself. For the state closest to threshold, the effect of this admixture turns out to be surprisingly small for volumes about twice as large as the state and was found to be less than . As the finite cutoff corrections for this state are about the same size, a more quantitative study requires an improved treatment of these corrections. Such studies are in progress.
Finally, our method should be extended to the threenucleon system. This requires also the inclusion of higher order corrections in the EFT and finite temperature effects as lattice calculations will inevitably be performed at a small, but nonzero, temperature. Work in these directions is in progress. With high statistics lattice QCD simulations of threebaryon systems within reach Beane:2009gs (), the calculation of the structure and reactions of light nuclei appears now feasible in the intermediate future. Our results demonstrate that the finite volume corrections for such simulations are calculable and under control. This also opens the possibility to test the conjecture of an infrared limit cycle in QCD for quark masses slightly larger than the physical values Braaten:2003eu ().
Acknowledgements.
We would like to thank D. Lee and B. Metsch for helpful discussions. This research was supported by the DFG through SFB/TR 16 “Subnuclear structure of matter” and the BMBF under contracts No. 06BN411 and 06BN9006.Appendix A Numerical treatment
In this appendix, we give some details on the numerical solution of Eqs. (17, 19). The starting point for a first numerical treatment of the formalism is the homogeneous integral equation
(21) 
The first step is to transform this equation into a finitedimensional problem. Due to the oscillatory nature of the integrand, it is not sensible to use a finite number of sampling points. Therefore, a set of basis functions will be used. The choice of the basis functions is guided by the knowledge of the boundstate amplitude in the infinite volume case. Asymptotically, the amplitude in the infinite volume behaves like Bedaque:1998kg (); Braaten:2004rn ()
with the universal number and a momentum scale . The boundstate amplitude obtained by the infinite volume formalism has precisely such a form, with only a few zeros in the interval (See, e.g., Ref. Bedaque:1998kg ()). Therefore, the basis functions are chosen as Legendre functions with logarithmic arguments as
(22) 
These basis functions are orthogonal with respect to a suitably chosen scalar product:
(23) 
The amplitude in Eq. (21) is replaced by its expansion in the basis functions and the th component is projected out using the scalar product. The resulting equation
(24) 
with
(25) 
can be interpreted as a matrix equation when truncating the set of basis functions.
A nontrivial solution of Eq. (24) can be found if the condition
(26) 
is fulfilled. Values of , for which this is the case, are then considered as trimer energies in the finite volume. This is only a valid interpretation if the truncation of the basis and the neglect of the higher partial waves induce only small corrections for the result.
The value of the parameter for which the condition (26) is satisfied is found via a root finding algorithm. This requires the calculation of in each iteration. To save numerical effort, it is possible to expand the kernel around the binding energy in infinite volume if the shift in the binding energy is small. The expansion is done up to first order:
(27) 
This expansion has been used to obtain the results in Kreuzer:2008bi (). It is applicable as long as the finite volume shift in the energy is not larger than 15–20% of the infinite volume energy.
In the following, details on the numerical methods used to calculate the kernel matrix elements as defined in Eq. (25) are presented. The integration over is performed using a logarithmically distributed GaussLegendre quadrature with 64 points.
The integrand of the integration is strongly oscillating due to the term. Therefore, this integral is evaluated using the Fast Fourier Transformation technique (FFT). The computation of Fouriertype integrals via FFT is explained in detail in NumRep:07 (). The sampling for the FFT is the most timeconsuming part of the calculation. The number of points to sample, , is determined by the highest “frequency” that is to be accessed. The frequencies in the present case are the values of . The highest possible frequency is connected to and the integration range by
(28) 
For realistic values of these parameters, namely box sizes of a few , an of about 3000 and cutoff values of several hundred , a typical value for is or about one million sampled points.
A summation over threedimensional integer vectors of a quantity depending only on the absolute values involves substantial double counting when naïvely done. To circumvent this problem, an ordered list with all absolute values that are possible for vectors in and their multiplicity has been created. The result of the summation can be viewed as converged when going to absolute values of about (see the example in Fig. 8). For radii from about 2000 on, the intermediate sums oscillate around the converged result. To reduce runtime, the result of the summation has been calculated as a mean of 15 intermediate sums for radii from 2750 to 3500.
The calculation of the matrix elements has been performed on the cluster of the HISKP at the University of Bonn. As already stated above, the number of sampling points for the FFT is the parameter with the largest influence on the runtime. When using the expanded integral kernel, the matrix and the derivative matrix have to be calculated only once. For a typical sample size of , this takes about 100 minutes. When going to smaller volumes, the shifts in the binding energy become larger and the Taylor expansion of the integral kernel breaks down. In this case, it is inevitable to recalculate the kernel matrix in each iteration of the root finding algorithm. This amounts to 10 to 15 evaluations of the kernel matrix and a runtime of 8 to 10 hours per data point.
References
 (1) E. Braaten and H.W. Hammer, Phys. Rept. 428, 259 (2006) [arXiv:condmat/0410417].
 (2) L. Platter, Few Body Syst. 46, 139 (2009) [arXiv:0904.2227 [nuclth]].
 (3) V. Efimov, Phy. Lett. 33B, 563 (1970).
 (4) V. Efimov, Sov. J. Nucl. Phys. 29, 546 (1979).
 (5) P.F. Bedaque, H.W. Hammer, and U. van Kolck, Phys. Rev. Lett. 82, 463 (1999) [arXiv:nuclth/9809025]; Nucl. Phys. A 646, 444 (1999) [arXiv:nuclth/9811046].
 (6) T. Kraemer, M. Mark, P. Waldburger, J.G. Danzl, C. Chin, B. Engeser, A.D. Lange, K. Pilch, A. Jaakkola, H.C. Nägerl, and R. Grimm, Nature 440, 315 (2006).
 (7) S. Knoop, F. Ferlaino, M. Mark, M. Berninger, H. Schoebel, H.C. Naegerl, R. Grimm Nature Physics 5, 227 (2009) [arXiv:0807.3306 [condmat]].
 (8) T. B. Ottenstein, T. Lompe, M. Kohnen, A. N. Wenz, S. Jochim, Phys. Rev. Lett. 101, 203202 (2008) [arXiv:0806.0587 [condmat]].
 (9) J. H. Huckans, J. R. Williams, E. L. Hazlett, R. W. Stites, K. M. O’Hara, Phys. Rev. Lett. 102, 165302 (2009) [arXiv:0810.3288 [physics.atomph]].
 (10) G. Barontini, C. Weber, F. Rabatti, J. Catani, G. Thalhammer, M. Inguscio, F. Minardi, Phys. Rev. Lett. 103, 043201 (2009) [arXiv:0901.4584v1 [condmat.other]].
 (11) N. Gross, Z. Shotan, S. Kokkelmans, L. Khaykovich, Phys. Rev. Lett. 103, 163202 (2009) [arXiv:0906.4731v1 [condmat.other]].
 (12) M. Zaccanti, B. Deissler, C. D’Errico, M. Fattori, M. JonaLasinio, S. Müller, G. Roati, M. Inguscio, G. Modugno, Nature Physics 5, 586 (2009) [arXiv:0904.4453v1 [condmat.quantgas]].
 (13) D. L. Canham and H.W. Hammer, Eur. Phys. J. A 37, 367 (2008) [arXiv:0807.3258 [nuclth]] and references therein.
 (14) K.G. Wilson, Nucl. Phys. Proc. Suppl. 140, 3 (2005) [arXiv:heplat/0412043].
 (15) S.R. Beane, P.F. Bedaque, M.J. Savage, and U. van Kolck, Nucl. Phys. A 700, 377 (2002) [arXiv:nuclth/0104030].
 (16) S.R. Beane and M.J. Savage, Nucl. Phys. A 717, 91 (2003) [arXiv:nuclth/0208021]; Nucl. Phys. A 713, 148 (2003) [arXiv:nuclth/0206113].
 (17) E. Epelbaum, U.G. Meißner, and W. Glöckle, Nucl. Phys. A 714, 535 (2003) [arXiv:nuclth/0207089].
 (18) E. Braaten and H.W. Hammer, Phys. Rev. Lett. 91, 102002 (2003) [arXiv:nuclth/0303038].
 (19) E. Epelbaum, H.W. Hammer, U.G. Meißner and A. Nogga, Eur. Phys. J. C 48, 169 (2006) [arXiv:hepph/0602225].
 (20) H.W. Hammer, D.R. Phillips and L. Platter, Eur. Phys. J. A 32, 335 (2007) [arXiv:0704.3726 [nuclth]].
 (21) S. R. Beane, P. F. Bedaque, K. Orginos and M. J. Savage, Phys. Rev. Lett. 97, 012001 (2006) [arXiv:heplat/0602010].
 (22) S. R. Beane et al., Phys. Rev. D 80 (2009) 074501 arXiv:0905.0466 [heplat].
 (23) S. R. Beane, K. Orginos and M. J. Savage, Int. J. Mod. Phys. E 17, 1157 (2008) [arXiv:0805.4629 [heplat]].
 (24) S. R. Beane, P. F. Bedaque, A. Parreno and M. J. Savage, Phys. Lett. B 585, 106 (2004) [arXiv:heplat/0312004].
 (25) S. Tan, Phys. Rev. A 78, 013636 (2008) [arXiv:0709.2530 [condmat]].
 (26) S. R. Beane, W. Detmold and M. J. Savage, Phys. Rev. D 76, 074507 (2007) [arXiv:0707.1670 [heplat]].
 (27) W. Detmold and M. J. Savage, Phys. Rev. D 77, 057502 (2008) [arXiv:0801.0763 [heplat]].
 (28) S. R. Beane, W. Detmold, T. C. Luu, K. Orginos, M. J. Savage and A. Torok, Phys. Rev. Lett. 100 (2008) 082004 [arXiv:0710.1827 [heplat]].
 (29) W. Detmold, M. J. Savage, A. Torok, S. R. Beane, T. C. Luu, K. Orginos and A. Parreno, Phys. Rev. D 78, 014507 (2008) [arXiv:0803.2728 [heplat]].
 (30) W. Detmold, K. Orginos, M. J. Savage and A. WalkerLoud, Phys. Rev. D 78, 054514 (2008) [arXiv:0807.1856 [heplat]].
 (31) L. Pricoupenko and Y. Castin, J. Phys. A 40, 12863 (2007) [arXiv:0705.1502 [condmat]].
 (32) M. Lüscher, Nucl. Phys. B 354, 531 (1991).
 (33) M. Lüscher, Nucl. Phys. B 364, 237 (1991).
 (34) V. Bernard, U.G. Meißner and A. Rusetsky, Nucl. Phys. B 788, 1 (2008) [arXiv:heplat/0702012].
 (35) V. Bernard, M. Lage, U.G. Meißner and A. Rusetsky, JHEP 08, 024 (2008) [arXiv:0806.4495 [heplat]].
 (36) S. Kreuzer and H.W. Hammer, Phys. Lett. B 673, 260 (2009) [arXiv:0811.0159 [nuclth]].
 (37) F. C. von der Lage and H. A. Bethe, Phys. Rev. 71, 612 (1947).
 (38) S. L. Altmann and A. P. Cracknell, Rev. Mod. Phys. 37, 19 (1965).
 (39) E. Braaten, H. W. Hammer and M. Kusunoki, Phys. Rev. A 67 (2003) 022505 [arXiv:condmat/0201281].
 (40) W. H. Press, S. A. Teukolsky, W. T. Vetterling, B. P. Flannery. Numerical Recipes: The Art of Scientific Computing, Third Edition, Cambridge University Press, 2007.