Instruction: Write an abstract for the following acadamic paper.

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

Paper: an optical lattice is a crystal structure of off - resonant laser light , wherein atoms are trapped due to the optical dipole force . this creates a periodic potential in which the motional states of a single atom exhibits many of the usual textbook solid state characteristics , such as energy band structure , bloch oscillations and wannier - stark ladders  @xcite . in an interacting many - body system the periodic potential constitutes an ideal setting for realizing a wide range of condensed matter phenomena such as the superfluid to mott insulator transition  @xcite . likewise controlled atomic collisions in an optical lattice potential is a promising approach to implement quantum gate operations  @xcite .    in the context of ultracold atomic collisions a feshbach resonance arises , when a bound state of a closed scattering channel is coupled to the scattering continuum in the entrance channel  @xcite . in an applied magnetic field the separated atoms are described by eigenstates of the hyperfine and zeeman - hamiltonian , and the coupling between the different channels is a result of the atomic interaction , which is not diagonal in the hyperfine basis . the position of the resonance is tunable by varying the strength of the applied magnetic field , since the open and closed channels have different magnetic moments , and therefore experience different zeeman shifts . when the closed channel bound state is tuned to be degenerate with the colliding atoms in the entrance channel the scattering cross section is resonantly enhanced . a feshbach resonance may thus be used to regulate the atomic interactions , which in the ultracold regime is characterized entirely by the @xmath0-wave scattering length . in addition , since the feshbach resonance provides an adiabatic connection between bound dimer states and the continuum , a pair of unbound atoms may be converted into a molecule by slowly ramping the magnetic field across the resonance position .    to study two - particle physics on a lattice experimentally , a state with either zero or two atoms at each lattice site may be prepared by adiabatic magnetic field sweeps across a feshbach resonance starting from a pure atomic cloud  @xcite . this prepares a state with at most two atoms in each lattice well , since sites with more than two atoms are eliminated by inelastic collisions . the unbound atoms may be subsequently removed by an optical purification pulse . the sample thus prepared contains only isolated pairs , whose binding can then be controlled by further adjustments of the magnetic field . a special situation arises in an optical lattice due to the band structure . the periodic potential leads to a structured continuum with scattering states grouped in energy bands separated by band gaps . the depth of the lattice controls the width of the bands and their mutual separation . in this paper we concentrate on the lowest band , where the band structure then implies that the continuum has both a lower and an upper edge . this has fundamental implications for both the bound states and the scattering states of the system . a counter - intuitive effect of the lattice band structure is the existence of bound pairs composed of atoms , which repel each other  @xcite . such repulsively bound pairs exist , when the repulsive interaction is strong enough to lift the state of two atoms on the same lattice site out of the continuum and into the band gap , thereby preventing the pair from dissociating and making it a bound state of the system  @xcite . in one dimension repulsively bound pairs exist for an arbitrarily small repulsion . near a feshbach resonance further intriguing two - body effects become possible . in particular , the feshbach resonance gives rise to motionally bound molecules , which are stable only at finite center of mass momentum with respect to the lattice , while decaying if at rest . it is thus possible to tune the location of the scattering resonance at a fixed magnetic field by controlling the molecular motion  @xcite . here we derive analytical results fully characterizing both the bound and scattering states of an atom pair in an optical lattice near a feshbach resonance . we represent the lattice physics in a discrete model with nearest neighbor hopping . to describe the feshbach resonance we adopt a two - channel model with one open channel and one closed channel , which is adequate to describe the physics of both the scattering continuum and the bound states of the system for an isolated resonance  @xcite . contrary to  @xcite we include the background interaction in the entrance channel , which plays an important role both in scattering and in the bound state dynamics . the direct scattering process in the open channel interferes with the resonant scattering proceeding via the closed channel resonance state leading to a fano profile in the transmission probability . in the bound state spectrum the presence of an additional bound state leads to avoided crossings that are important when considering adiabatic transitions induced by ramping the magnetic field . the same fundamental physics of a discrete state embedded in a structured continuum occurs in the fano - anderson model , which describes autoionizing states in atoms  @xcite and localized electron states in solids  @xcite . from a many - body perspective the combination of an optical lattice with resonant control of the atomic interactions makes a variety of spin models attainable  @xcite .    in a previous paper @xcite we described the basic results of the discrete lattice model of the feshbach resonance physics , and suggested how the salient features of the system may be probed in experiments . here we present the full analytic theory of feshbach scattering in a one - dimensional lattice ( sec . [ sec : scattering ] ) . we introduce a generalized , one - dimensional scattering length , which is defined at both the upper and the lower threshold of the lowest continuum band in the lattice . this generalized scattering length shows a characteristic feshbach divergence when a bound state of the system crosses either band edge , and we identify the magnetic field positions for these divergences as a function of the lattice depth for a specific feshbach resonance . in sec . [ sec : bound - states ] we consider the bound state solutions of the combined lattice - feshbach problem and derive analytic expressions for the binding energies and the wavefunctions of the dimer states . we also show that the molecules exhibit cooperative tunneling with a tunneling rate consistent with that measured in experiments  @xcite . a spectral analysis of the unbound and dimerized atom pair states is presented in sec . [ sec : spectral - analysis ] , where we demonstrate that threshold effects lead to a non - trivial peak structure of the molecular spectral function inside the band . our discrete lattice model includes only the lowest continuum band and neglects couplings to higher bands in the lattice . in sec . [ sec : validity ] we discuss when the lowest band approximation is valid and how this restriction may be relaxed . we find that a moderate interband coupling only modifies our results on a quantitative level . throughout we give quantitative results for a specific feshbach resonance for which our model is valid  @xcite . we consider atoms with mass @xmath1 in a cubic optical lattice potential @xmath2 created by off - resonant laser light with wavelength @xmath3 . the lattice period is , @xmath4 , and the transverse , @xmath5 , and longitudinal depths , @xmath6 , are measured in units of the photon recoil energy , . for a cubic lattice the motion separates along the three spatial directions , and the hamiltonian , @xmath7 , is diagonalized by products of three bloch functions , @xmath8 , one for each spatial direction @xmath9 , where @xmath10 is the quasi - momentum in the range @xmath11 $ ] , and @xmath12 is the band index . the energy eigenstates form the usual band structure , splitting the continuum into bands separated by energy gaps that decrease with increasing @xmath12 . the band gaps grow and the widths of the bands decrease , as the lattice is made deeper by increasing the laser intensity . we will assume that the lattice depth along the transverse directions , @xmath5 , is much larger than the depth of the longitudinal lattice , @xmath6 . this allows us to only include the ground state along the two perpendicular directions , leaving an effectively one - dimensional lattice . as the scattering energy is increased consecutive transverse channels open , however we only consider scattering at energies much smaller than the transverse level splitting , thus justifying the neglect of the transverse degrees of freedom . a confinement induced resonance occurs at the lower edge of the continuum when the ( free space ) scattering length @xmath13 and the width of the transverse confining potential @xmath14 are related by @xmath15  @xcite . however , this resonance may be avoided by not choosing the transverse confinement too tight , and for the quantitative results we present here , it plays no role . the bloch waves are delocalized , but can be combined into orthonormal localized basis functions , wannier functions , which in the longitudinal direction take the form @xmath16 in addition to the band index they are labeled by a site index , @xmath17 , specifying at which lattice sites , @xmath18 , the wannier functions are localized . a similar wannier basis can be constructed in the two transverse directions . when it is not clear from the context , which wannier basis we refer to , we will add a superscript @xmath19 or @xmath20 to the wannier functions indicating to which lattice direction they belong . the dynamics of a single atom in the @xmath12th longitudinal band , _ i.e. _ the kinetic energy and lattice potential , is characterized by the tunneling amplitude , , between lattice sites separated by a distance @xmath21 . here @xmath22 only depends on the coordinate along the weak lattice direction , and @xmath1 is an integer . in the transverse directions tunneling can be disregarded , due to the much larger lattice depth .    in this work we consider only the motion in the lowest band of the lattice , and we disregard tunneling between non - adjacent sites , assuming @xmath23 for all . in the basis of the wannier states the single particle hamiltonian for motion in the longitudinal direction then becomes @xmath24 with @xmath25 . we have set the zero of energy to be the expectation value of the energy for an atom in a lowest band wannier state , _ i.e. _ @xmath26 . the hilbert space is now represented by a discrete basis of the localized wannier functions , and accordingly we have translated the continuum description of the system into a discrete lattice model . by our definition @xmath27 corresponds to the center of the lowest longitudinal bloch band , and the single - particle energies of the discrete lattice model are @xmath28 . the approximations made are the omission of higher bands , and the neglect of tunneling beyond nearest neighbor sites . both of these restrictions may be relaxed in a systematic way  @xcite . we use the longitudinal wannier functions as a basis for the two - particle states on the lattice . hence the two - particle wavefunction @xmath29 will denote the amplitude of finding the atoms in the wannier orbitals @xmath30 and @xmath31 localized around the discrete positions @xmath32 and @xmath33 . considering only the lowest bloch band we have @xmath34 .    to analyze the motion of the two particles in the lattice we introduce center of mass , @xmath35 , and relative , @xmath36 , coordinates . for the center of mass motion a plane - wave ansatz @xmath37 is possible , with the first brillouin zone for the center of mass momentum , @xmath38 , running from @xmath39 to @xmath40 . this separation is pertinent , since states with distinct center of mass momenta are not coupled by the interactions or the lattice , and we can thus obtain separate results for the relative motion for each @xmath38 . taking the interaction potential @xmath41 to be purely on - site with strength @xmath42 ( defined in appendix  [ sec : calc - param ] ) , the relative motion of two atoms in the open channel is then described by the hamiltonian , @xmath43 , where the relative motion lattice hamiltonian , @xmath44 , contains the discrete laplacian  @xcite @xmath45 , \label{lattice_laplacian}\ ] ] and the center of mass energy , @xmath46 . the eigenstates of @xmath47 are discrete plane waves , @xmath48 . the corresponding energies , @xmath49 , are the sum of the single - particle energies of atoms with momenta @xmath50 , @xmath51 $ ] . the ( retarded ) green s function for non - interacting atoms on the lattice is @xmath52^{-1}$ ] , where @xmath53 is a positive infinitesimal . in momentum space it is diagonal @xmath54 while its coordinate space form may be found from @xmath55 due to the delta function interactions we only need to consider @xmath56 . for energies inside the band ( @xmath57 ) the solution is propagating @xmath58 with @xmath59 , while the the solution outside the band ( @xmath60 ) falls off exponentially : @xmath61^{z / a}. \label{g0_z_out}\ ] ] here @xmath62 . for energies above the band the sign of @xmath63 alternates between lattice sites .    the green s function related to the open channel hamiltonian , including the interaction , @xmath41 , is found by explicitly solving the dyson equation @xmath64 in coordinate space the result has the simple form @xmath65 defining the diagonal elements of the non - interacting momentum space green s function to be @xmath66^{-1}$ ] the interacting open channel green s function has the momentum space form @xmath67 which follows from the matrix elements of ( [ dyson ] ) in the relative momentum eigenbasis , using that for an on - site interaction the matrix element , @xmath68 , is constant for all @xmath69 .      with the expressions for the green s functions the formal solution of the schrdinger equation for the relative motion becomes explicit . below we derive the transmission probability for open channel collisions and analyze the near threshold scattering in terms of a generalized scattering length . this section introduces the concepts and quantities central to the ensuing discussion of two - channel scattering . we first present results relevant to scattering of distinguishable particles , and then discuss how to treat the case of identical bosons or fermions . choosing the situation with the incident wave entering from the left with quasi - momentum @xmath59 , the scattering state of two particles colliding under the influence of the entrance channel potential , @xmath41 , is given by the solution of the lippmann - schwinger equation @xmath70 for our on - site interaction of strength , @xmath42 , the scattering wavefunction becomes @xmath71 where we have used that @xmath72 . we note that for the contact potential @xmath73 equals its asymptotic form @xmath74 everywhere , hence we may readily identify the scattering amplitude for the background interaction : @xmath75 which is unit - less in a one - dimensional system . for the one - dimensional problem and distinguishable particles it is natural to recast this in terms of transmission and reflection amplitudes . these are defined from the left and right asymptotic form of the wavefunction @xmath76 by comparing with ( [ psi_e_bg_asymp ] ) it follows that the transmission amplitude is @xmath77 , while the reflection amplitude coincides with the scattering amplitude , @xmath78 . the resulting transmission probability is @xmath79 as one might expect , the transmission vanishes in the limit where @xmath80 , and approaches unity as the strength of the on - site interaction is diminished . as usual , the reflection coefficient is @xmath81 , which is unity at the edges of the band , where propagation ceases . the minimum value of the reflection probability , @xmath82 , is attained in the middle of the band , @xmath27 . we emphasize that the transmission is independent of the sign of the interaction . this is a special property of the delta - function potential .    as usual the scattering amplitude may also be related to the phase shift , @xmath83 , of the scattered wave . in one dimension the relation is  @xcite @xmath84 in general there will be two partial waves , corresponding to even and odd solutions of the relative motion problem . for the delta - function interaction the odd partial wave has vanishing phase shift due to symmetry , and the reflection and transmission probabilities are then given by @xmath85 and @xmath86 , respectively . the preceding derivation of the one - dimensional scattering in the lattice assumed that the two colliding atoms are distinguishable . if the collision involves identical bosons or fermions the scattering state @xmath87 must be an even or an odd function of the separation @xmath88 , respectively , to make the wavefunction either symmetric or anti - symmetric under interchange of the two atoms . for fermionic atoms in the same internal state this implies that the wavefunction vanishes for @xmath89 and the contact interaction does not produce any scattering .    for identical bosons the properly symmetrized scattering state may be written as @xmath90 when identical particles collide in one dimension the transmission and reflection coefficients have no meaning , since the incoming and outgoing fluxes are identical , however the scattered wave still experiences a phase shift given by ( [ eq : f_and_delta ] ) . even though the transmission profile calculated above does not play a role in the scattering of identical bosons , it may nonetheless be mapped out spectroscopically . consider transitions from a deeper bound state @xmath91 in the molecule to an energy eigenstate inside the band induced by applied rf- or photodissociation fields , which are described by some transition operator @xmath92 . if the initial state is well localized with respect to the lattice spacing , such dissociative transitions towards the state ( [ psi_bg_sym ] ) occur with a probability @xmath93 , obtained above .      when the relative quasi - momentum approaches the center or the edges of the brillouin zone , the two - particle energies @xmath94 tend to @xmath95 and the scattering amplitude attains a limiting form @xmath96 with @xmath97 at the bottom of the band and @xmath98 at the top of the band . these limits allow a natural definition of a generalized one - dimensional scattering length in the lattice @xmath99 from the limit @xmath100^{-1}$ ] of the scattering amplitude as @xmath101  @xcite . an equivalent definition is through the derivative of the phase shift with respect to the relative momentum @xmath102 in analogy with the usual free space scattering in three dimensions . our lattice scattering length is a generalization of the usual concept , in that it is defined for collision energies at both the lower _ and _ the upper edge of the continuum . the sign of @xmath103 at the upper edge of the continuum is a convention , which will be justified for the two - channel case below . our generalized scattering length depends on the center of mass motion of the pair , which is a crucial feature in the lattice , especially when we discuss the two - channel case below .    for @xmath104 a bound state is situated below the continuum , making @xmath103 positive as expected . in this case the pole of ( [ f_bg_lim ] ) lies along the positive imaginary axis , @xmath105 , and the analytic continuation of the plane wave @xmath106 is the dying exponential of the bound state . conversely , for a repulsively bound pair state lying above the continuum , @xmath107 , and @xmath103 is negative . the pole of the scattering amplitude is then at @xmath108 , corresponding to a wave function decaying exponentially with the separation between the two atoms and a phase factor @xmath109 , which alternates between @xmath110 and @xmath111 from one lattice site to the next , in accordance with the behavior of the green s function for the relative motion at @xmath112 ( [ g0_z_out ] ) .    as @xmath113 the scattering length diverges , and the bound state approaches the edge of the continuum . our intuition from scattering theory in three dimensions about the relation between weakly bound states and the scattering length thus holds both below and above the continuum band however , since @xmath42 is proportional to the free space background scattering length , @xmath114 , the generalized one - dimensional scattering length and its three - dimensional equivalent are inversely related  @xcite . we now consider scattering in the presence of a feshbach resonance mixing the entrance channel with an energetically closed channel through a coupling , @xmath115 . this gives rise to a set of two coupled equations for the relative motion of an atom pair    @xmath116    @xmath117    due to the large splitting between the bound states in the closed channel the resonance physics is faithfully represented by considering only a single bound state , the resonance state @xmath118 , of @xmath119 with an energy @xmath120 with respect to the center of the open channel band . this amounts to a single - pole approximation for the closed channel green s function @xmath121 the resonance state energy @xmath120 is composed , in part , of the free space resonance energy @xmath122 , which is linearly tunable with an applied magnetic field , due to the difference @xmath123 between the separated atoms and the closed channel magnetic moments , and , in part , of the displacement of the closed channel relative to the open channel bloch band as illustrated in fig . [ fig : energydiagram ] . at all magnetic field strengths we take the zero of energy to be the center of the open channel continuum band . in free space the bare resonance energy can be related to the magnetic field position of the zero energy resonance , @xmath124 , as outlined in appendix  [ resonance_shift ] . is the energy of @xmath125 relative to the open channel threshold . applying an optical lattice changes the position of the open and closed channel by @xmath126 and @xmath127 , respectively , where @xmath128 is the average atomic energy and @xmath129 is the molecular spectrum of the first bloch band both along direction @xmath91 . note that the position of the open channel thresholds in the lattice depends on the quasi - momentum , @xmath38 , and therefore we choose @xmath27 to be the center of the open channel band . hence , the thresholds are located at @xmath95 . the difference in energy between the closed channel bound state and the center of the open channel band in the presence of the lattice is @xmath130 . ]    as a function of continuous position variables the longitudinal lattice potential for each atom is @xmath131 with @xmath132 , and therefore the combined lattice potential for a pair of atoms in the resonance state is @xmath133 where we have introduced the center of mass , @xmath134 , and relative , @xmath135 , coordinates . since the size of the closed channel bound state is much smaller than the lattice spacing , we may consider the limit of @xmath136 for this state . consequently , the resonant state has a center of mass motion corresponding to a lattice potential twice as deep as that governing the motion of the individual atoms , while its mass is the sum of the two atomic masses . this leads to a tunneling rate for the resonant state , @xmath137 , which is much reduced in comparison with the rate at which the free atoms tunnel through the lattice , @xmath138 . this bare molecular tunneling leads to a modulation of the resonance energy with the center of mass momentum of the pair , which we may parametrize as @xmath139 the band averaged resonance energy , @xmath140 , incorporates the relative energy shifts of the atomic and molecular bands in the lattice .      with the explicit expression for the closed channel green s function ( [ g_cl_approx ] ) the closed channel part of the scattering problem ( [ coupled_eq2 ] ) may be solved formally : @xmath141 where we have introduced a coupling matrix element , @xmath142 , incorporating the structure of the resonance state . for a typical resonance in alkali collisions the coupling arises due to the difference between the triplet and singlet molecular potentials , and it is therefore inherently of short range compared with the lattice spacing . it is therefore reasonable that the coupling only acts between atoms localized at the same lattice site . the open channel component of the scattering state in the presence of the resonance may then be expressed as @xmath143 \label{psi_op_res}\ ] ] for all @xmath144 . again the relative quasi - momentum of the scattering atoms is @xmath59 , and @xmath145 is the molecular self - energy , which for the onsite coupling may be written as : @xmath146 .    from the form of the scattered wave ( [ psi_op_res ] ) the total transmission coefficient including both background and resonant contributions is found to be of the fano form  @xcite @xmath147 the fano parameters @xmath148 and @xmath149 depend on the resonance shift and width functions , which are related to the real and imaginary parts of the molecular self - energy , @xmath150 where the real part @xmath151 describes the self consistently determined shift of the resonance position from @xmath152 to @xmath153 , and @xmath154 is the decay width of the resonance state due to the coupling to the continuum . outside the continuum @xmath155 , and the molecules are stable against dissociation . the real part of the molecular self - energy outside the band @xmath156 determines the energies of the bound states for a given resonance energy as described in section  [ sec : bound - states ] . the width and shift of the resonance are illustrated in fig . [ fig : shift_and_width ] along with the fano parameters for @xmath157 . throughout the paper we illustrate our analytical results for a set of parameters where @xmath158 and @xmath159 . as outlined in appendix  [ sec : calc - param ] this corresponds to @xmath160rb atoms in a lattice with longitudinal and transverse depths of @xmath161 and 30@xmath161 , respectively , in an applied magnetic field tuned close to the feshbach resonance near 414 g  @xcite .     and width , @xmath162 of the resonance for @xmath157 where the band edges are at @xmath163 ( a ) . the real part of the molecular self - energy is continuous across the band edges and has a pole at the uncoupled bound state energy , @xmath164 ( see section  [ sec : bound - states ] ) . the fano parameters , @xmath165 and @xmath10 , parametrize the shape of the transmission resonance ( b ) . here @xmath166 . ] we note that @xmath167 is independent of the coupling , @xmath168 . it is related to the background transmission probability through the relation @xmath169 , which is a consequence of inversion symmetry in the problem  @xcite . the reflection maximum of the fano profile occurs at @xmath170 and is nearly independent of the center of mass motion , due to the weak tunneling of the closed channel molecules . this point of total reflection only occurs , when the bare resonance state lies inside the continuum . just as in the non - resonant case the reflection probability is unity at the band edges . hence the reflection profile has two minima , when @xmath171 . when the bare resonance state lies outside the band , @xmath172 has a single minimum . the two - channel reflection profiles are plotted in figs . [ fig : r_vs_k](a)-(d ) and [ fig : r_vs_eres](a ) as colormaps and in figs . [ fig : r_vs_k](e)-(h ) and [ fig : r_vs_eres](c ) as solid lines . for the parameters used in the quantitative calculations throughout this paper @xmath173 , and hence the resonance energy appears nearly independent of @xmath38 in figs . [ fig : r_vs_k](a)-(d ) .     for a fixed @xmath157 . the horizontal dashed line is the repulsively bound pair state existing above the band for @xmath174 , while the diagonal dashed curve is the band averaged resonance energy , @xmath175 . ( b ) closed channel population of the bound states . ( c ) cuts showing the reflection probability at @xmath176 0 , -3 , 10 , and @xmath177 ( solid lines from top to bottom at @xmath27 ) . the dashed curve indicate the reflection profile in the absence of the coupling . ] the condition for total transmission is @xmath178 , corresponding to @xmath179 . as @xmath152 is ramped upwards ( downwards ) starting below ( above ) the band , this point of destructive interference between the direct and resonant scattering amplitudes first appears at the magnetic field , where a bound state of the coupled channels system enters the continuum to become a scattering resonance . as we shall see below , this occurs at two critical values where @xmath180 - exactly the points delineating the region , where the condition @xmath178 can be satisfied . the reflection profile will thus only have both a fano maximum and a fano minimum if the conditions @xmath171 and @xmath181 are simultaneously met .    neglecting the background interaction corresponds to the limit @xmath182 , or equivalently @xmath183 . in this case the resonance shift vanishes and the reflection probability assumes a simple breit - wigner - like profile  @xcite @xmath184 ^ 2+\frac{\hbar^2\gamma_0 ^ 2(e , k)}{4}},\ ] ] where @xmath185 . the opposite limit of @xmath186 on the other hand corresponds to @xmath187 and @xmath188 , giving rise to a breit - wigner - like form of the transmission coefficient @xmath189 ^ 2+\frac{\hbar^2\gamma_\infty^2(e , k)}{4}},\ ] ] corresponding to a dip in the reflection probability at the resonance , which now occurs at @xmath190 . the resonance shift , @xmath191 , is independent of energy and momentum , while the width is given by @xmath192 . both the shift and the width approach zero as the strength of the onsite interaction is increased . the observations regarding indistinguishable particles in section  [ sec : identical ] can be directly carried over to the situation with more than one scattering channel . again , the transmission coefficient @xmath193 has no straightforward meaning in the scattering of identical bosons . but following fano  @xcite it is clear that the transition probability from an initial state @xmath91 to the stationary scattering state @xmath194 under the action of some transition operator @xmath92 is @xmath195 . therefore , the resonance profile is accessible spectroscopically by scanning the energy of the final state across the continuum . starting from a deeper bound state of the system the dissociated atom signal will have the same energy variation as @xmath193 . the scattering amplitude in the coupled channels case consists of a background and a resonant contribution , @xmath196 , where the resonant part @xmath197 ^ 2}{e - e_{\rm{res}}(b , k)-\sigma_m(e , k ) } \label{f_res}\ ] ] is obtained by comparing ( [ psi_op_res ] ) with @xmath198 for @xmath199 . the phase shift is related to the scattering amplitude in the same way as in the single - channel case .    as in the case of single - channel scattering we define the one - dimensional scattering length , @xmath200 , from the limit of the scattering amplitude , @xmath201^{-1}$ ] , at the extrema of the relative energy , @xmath202 . after some algebra we arrive at the simple expression @xmath203 \label{a_1d}\ ] ] for the scattering length , with the upper and lower sign corresponding to the top and the bottom edge of the continuum , respectively . this lattice scattering length is of the standard feshbach form @xcite , and reduces to the non - resonant expression ( [ a_bg ] ) , when the coupling @xmath168 is set to zero . however , the generalized , magnetic field dependent scattering length at the top and the bottom of the continuum differ , in that they do not diverge at the same value of the resonance energy . the motivation for the chosen sign convention of the scattering length at the upper band edge is guided by the rationale that the two scattering lengths then become equal in the limit of a deep lattice , where @xmath204 goes to zero .    ) . solid lines show @xmath205 at @xmath206 0 , 0.4 , and 0.7 @xmath40 ( from the left to the right ) , while the dashed lines correspond to @xmath207 at the same @xmath38-values ( increasing from the right to the left ) . ( b ) the corresponding binding energies below the band ( solid lines ) and above the band ( dashed lines ) and the limiting form ( [ e_bind_universal ] ) , valid near the feshbach resonances in @xmath208 ( thin dot - dashed lines ) . ]    it is convenient to define a pair of detunings @xmath209 , where @xmath210 . we thus get two resonance positions , one for scattering near the upper band edge , and one for scattering near the lower band edge , as shown in fig . [ fig : a_1d ] . we caution that the relative position of @xmath211 and @xmath212 depends on the sign of the magnetic moment difference , @xmath123 . if @xmath213 ( @xmath214 ) , then @xmath211 lies above ( below ) @xmath212 . the lattice resonance positions @xmath211 and @xmath212 are shown in fig . [ fig : b0 ] as function of the lattice depth , @xmath215 . they are separated by @xmath216 and approach each other as @xmath215 is increased , since @xmath217 decreases when the lattice becomes deeper . furthermore , in the limit where both the transverse and the longitudinal lattice depths are large , _ i.e. _ , @xmath218 , the wannier orbitals approach harmonic oscillator states , and we have that @xmath219 , while @xmath220 and @xmath221 . this gives the following asymptotic behavior of the magnetic field resonance position @xmath222   + \frac{e_{\rm        r}}{\delta\mu}\sqrt{v_0^{\parallel } }   + \frac{2e_{\rm{r}}}{\delta\mu}\sqrt{v_0^\perp}\ ] ] where @xmath223 is the feshbach resonance position in the absence of the lattice , and the dimensionless function @xmath224 is defined in appendix [ resonance_shift ] . the asymptotic behavior is indicated by the black dash - dotted curves in fig . [ fig : b0 ] .    . @xmath211 and @xmath212 are plotted for @xmath225 ( dashed curves ) and @xmath226 ( solid curves ) . in the example @xmath227 is positive so both @xmath211-curves lie above their respective @xmath212-curves . the left axis shows @xmath228 relative to the feshbach resonance position in the absence of the lattice , and the right axis shows the actual magnetic field position for the specific feshbach resonance considered  @xcite . the black dash - dotted curves are given by eq . ( [ eq : asymptoticb ] ) . ]    as shown in sec . [ sec : binding_energy ] the divergence of the lattice scattering length is associated with a bound state entering or leaving the continuum . hence the lattice gives rise to a @xmath38-dependent shift of the magnetic field threshold for the formation of two - body bound states away from the resonance position in free space . this fact , illustrated in fig . [ fig : a_1d ] , was also observed by orso _ et al . _ @xcite who solved the two - body physics in a one - dimensional optical lattice within a single - channel model . the zero - crossing of the scattering length occurs for @xmath229 , which corresponds to magnetic field values , @xmath230 , shifted from the resonance positions . the bound state solutions of the coupled - channels problem have energies , @xmath231 , coinciding with the poles of the scattering amplitude @xmath232 . the pole locations satisfy the equation @xmath233=\frac{{\mathcal{w}}^2{\rm{sgn}}(e_b)}{e_b - e_{\rm{res}}(b , k ) } , \label{e_bind}\ ] ] and can be determined as a subset of the roots of a quartic polynomial . the bound state energies are plotted as the solid lines in figs . [ fig : r_vs_k](a)-(d ) and [ fig : r_vs_eres](a ) . if @xmath234 the entrance channel bound state energy is @xmath235 with a wavefunction proportional to @xmath236 . this single - channel bound state is indicated by the sinusoidal dashed line in fig  [ fig : r_vs_k](a)-(d ) and the horizontal dashed line in fig . [ fig : r_vs_eres](a ) .    an important consequence of the structured continuum in the lattice potential is that the stability of a molecular bound state may depend on its center of mass momentum . for a given magnetic field , the scattering resonance only exists for a range of center of mass motion bloch states , as shown in fig . [ fig : r_vs_k](c)-(d )  @xcite . outside this range there is instead a true bound state of the system . such motionally bound states are stabilized by their kinetic energy , which displaces the pair state into the band gap for the relative motion , where the pair can not disintegrate . hence , for a fixed magnetic field , where a molecule at rest is unstable , accelerating the lattice provides a new avenue for tuning the feshbach resonance . the bound state pairs can tunnel together through the lattice with an effective tunneling rate , which depends on their binding energy . we now determine the tunneling matrix element of the pairs directly from their dispersion relation .    in the absence of atomic tunneling processes , @xmath237 , and the energies of the bound states are given by @xmath238 ^ 2}{4}+{\mathcal{w}}^2}. \label{ebound_no_tunneling}\ ] ] with the inter - channel coupling turned off these zero tunneling bound state energies reduce to the uncoupled pair state energies for @xmath239 ; @xmath175 and @xmath42 . with a small non - vanishing single - particle tunneling amplitude we may expand ( [ e_bind ] ) to lowest order in @xmath240 and find @xmath241 comparing the momentum dependent part of the bound state energy with the dispersion relation for non - interacting pairs , @xmath242 , we see that in the limit of a deep lattice the bound atom pairs tunnel through the lattice with an effective tunneling rate @xmath243 as one would expect from arguments based on second order perturbation theory  @xcite . we emphasize that the molecule dispersion depends on the applied magnetic field . cooperative tunneling of atom pairs consistent with an effective pair hopping matrix element given by @xmath244 has been demonstrated experimentally by measurement of the frequency of population transfer in a double well  @xcite and by demonstrating that the transport properties of an attractive fermi gas in an optical lattice depends strongly on the formation of local pairs  @xcite .      instead of the bound state energy it is often instructive to consider the _ binding energy _ , defined as the distance of the bound state energy from the ( upper or lower ) band edge : @xmath245 the binding energy is positive for a bound state lying above the band and negative for a bound pair below the continuum . from fig . [ fig : r_vs_k](a)-(d ) we observe that the binding energy increases with increasing @xmath246 in agreement with the single - channel calculation in  @xcite , which only discussed bound states with negative binding energy as no transverse confinement was included .    as @xmath231 approaches either of the band edges , the binding energy vanishes in a particularly simple way . to see this we observe that @xmath247 for @xmath248 . replacing @xmath231 with @xmath249 in the denominator on the right hand side in ( [ e_bind ] ) it then follows that @xmath250 ^ 2},\ ] ] as the bound state energy approaches @xmath249 . the denominator is exactly the resonance shape of the generalized scattering length ( [ a_1d ] ) , and the binding energy therefore approaches the universal form @xmath251 when the scattering length is tuned to be large with respect to the lattice spacing @xmath252 in the vicinity of the feshbach resonance . we have introduced the effective _ reduced _ mass @xmath253 of the atom pair with center of mass quasi - momentum , @xmath38 , in the lattice . its limits at the top and the bottom of the band are @xmath254 , with the upper ( lower ) sign applying for @xmath255 ( @xmath256 ) . these limits are plotted in fig . [ fig : m_eff ] as a function of @xmath38 . from the expression for the effective reduced mass it is clear that the sign of the universal binding energy expression ( [ e_bind_universal ] ) adheres to the definition above . the criteria of applicability for the quadratic approximation is @xmath257 , or @xmath258 . the universal expression ( [ e_bind_universal ] ) is compared with the exact binding energies in fig . [ fig : a_1d](b ) both above and below the continuum . due to the proximity of the zero crossing of @xmath259 , where ( [ e_bind_universal ] ) diverges , the agreement below the band is worse than above the band .    from eqs . ( [ a_1d ] ) and ( [ e_bind_universal ] ) it is clear that the binding energy vanishes when @xmath260 . this condition for the disappearance of the feshbach bound state coincides with the critical resonance energy at which the fano profile develops a point of total transmission , as anticipated above . experimentally the binding energy of the dimer states can be measured by modulating the depth of the lattice . the bound pairs dissociate , when the frequency of modulation satisfies a resonance condition , such that the dimers can exchange energy with the undulating lattice . the resonance is marked by the appearance of unbound atoms  @xcite . the formal solution of the bound state problem consistent with the single - pole approximation for the closed channel green s function is @xmath261 where the closed channel population , @xmath262 , is determined by the normalization condition : @xmath263 . in coordinate space the open channel component of the bound state wavefunction therefore becomes @xmath264 while the closed channel component of the bound state is described by @xmath265 . here we have used that the coupling is of zero range , and that the bare resonance state is localized on a single lattice site . the closed channel weight of a bound atom pair with energy @xmath231 is @xmath266 ^ 2}\frac{|e_b|}{(e_b^2-e_k^2)^{3/2}}\right]^{-1 } , \label{z}\ ] ] which is plotted in figs  [ fig : r_vs_k](i)-(l ) and  [ fig : r_vs_eres](b ) . as the bound state approaches the upper or the lower band edge , it follows from the limiting form of the binding energy that @xmath267 , and consequently @xmath262 vanishes as @xmath268 in the limit where @xmath269 .    in momentum space the open channel part of the bound state wavefunction is given by @xmath270 which is peaked near the center of the brillouin zone ( @xmath271 ) when the bound state energy is negative , and peaked near the edges of the brillouin zone ( @xmath272 ) for @xmath273 , making it possible to discern the bound states with @xmath274 from those with @xmath273 in time of flight imaging  @xcite . this reflects that @xmath275 is predominantly comprised of the @xmath276-states with energies , @xmath94 , in proximity of @xmath231 . for @xmath273 ( @xmath274 ) this corresponds to the states near the top ( bottom ) of the band . the height of the peaks is reduced as @xmath38 approaches @xmath277 , where all relative motion quasi - momentum states are degenerate in our model . the closed channel component of the bound state is independent of the relative momentum : @xmath278 .    , where a weakly bound molecular state ( solid line , narrow bars ) at rest is situated just below the band . the repulsively bound pair state lying above the band ( dashed line , thick bars ) , is peaked at the edge of the brillouin zone , and its sign alternates from one site to the next . the first column shows the momentum space wavefunction ( [ phi_mom ] ) in units of @xmath279 , the second column @xmath280 . ]    the size of the bound state @xmath281 ^ 2}\frac{e_k^2}{(e^2_b - e^2_k)^2},\ ] ] reflects the magnitude of the binding energy of atom pairs . in particular , it diverges when the bound state enters the continuum from above or from below : @xmath282 in this limit the open channel part of the bound state wavefunction assumes a universal form @xmath283^{z / a}.\ ] ]    the bound state wavefunctions are shown in fig . [ fig : bound_state_wf ] in both coordinate and momentum space for @xmath284 . at this resonance energy a bound molecular state exists just below the band for @xmath157 ( lower bound state branch in fig .  [ fig : r_vs_k](b ) ) . with a small binding energy , this state extends over several lattice sites as shown in fig  [ fig : bound_state_wf ] . for the repulsively bound pair state lying above the band ( upper bound state branch in fig . [ fig : r_vs_k](b ) ) the sign of the coordinate space wavefunction alternates between lattice sites . we have previously shown that a molecule has a finite probability of surviving a sweep through the continuum band when the applied magnetic field is varied across the feshbach resonance  @xcite . the survival probability increases the faster the sweep and the weaker the coupling to the continuum . given the existence of motionally bound states , an alternative experiment to investigate the passage of a molecule through the continuum may then be conducted by starting with the applied magnetic field held at a value such that @xmath285 , and molecules at rest prepared in the bound state with energy below @xmath175 ( represented by the the lower branch of bound states in figs .  [ fig : r_vs_k ] and  [ fig : r_vs_eres ] ) . accelerating the lattice imparts momentum @xmath286 to the molecules . if the magnetic field is then ramped adiabatically to a new value , where @xmath287 the moving molecules remain stable , while molecules at rest are now subject to dissociative decay . with a further acceleration of the lattice , either in the opposite direction , or in the same direction , but letting @xmath38 cross one or several brillouin zones , the molecules may be swept through the continuum to a final state with center of mass momentum @xmath288 , where the molecules are stable provided @xmath289 . the probability of a molecule surviving the center of mass momentum ramp through the continuum depends on the rate of lattice acceleration and the strength of the coupling to the continuum states . in this section we will perform a spectral analysis based on the full green s functions of both the open channel atoms and closed channel molecules . in general , given a retarded green s operator , @xmath290 , the corresponding spectral function is the imaginary part of the diagonal green s function , @xmath291 , in a basis @xmath292 of orthonormal states .      for the unperturbed system ( @xmath42=@xmath168=0 ) the spectral function is @xmath293 where the density of states is @xmath294 for @xmath57 . the unperturbed system has no bound states , and hence the spectral function is zero at all energies outside the continuum band . inside the band the scattering states are equally distributed over the whole lattice , making the spectral function translationally invariant and directly proportional to the density of states , which diverges at the band edges , @xmath295 .      both for the single - channel and the two - channel case we will need the entrance channel green s function , @xmath296 , _ i.e. _ , the propagator from one relative coordinate in the lattice , @xmath297 , to another , @xmath88 . when we consider just a single scattering channel this can be obtained from the dyson equation ( [ dyson ] ) : @xmath298 where we have used that the unperturbed green s function is translationally invariant . this expression for @xmath299 reduces to eq . ( [ gu_z ] ) when @xmath88 or @xmath297 is equal to zero .    for the spectral function we have to treat @xmath57 and @xmath60 separately . outside the continuum @xmath63 is real and we obtain the spectral function @xmath300 where @xmath301 is the decay constant for the bound state . this gives the probability distribution for finding a bound state atom pair @xmath88 lattice sites apart , when the bound state energy is @xmath302 . inside the band we find @xmath303 where @xmath73 is the scattering wave defined in eq . ( [ psi_e_bg ] ) . changing the sign of @xmath88 is equivalent to changing the direction of the incoming wave in the scattering wave from @xmath276 to @xmath304 . these two waves are degenerate and therefore it is natural to expect the average of the two for a given energy in the spectral function . we have used that the density of states is unchanged for elastic scattering , _ i.e. _ @xmath305 is also the density of states for atoms scattering under the influence of @xmath41 . it is instructive to calculate the total spectral weight residing inside the continuum band : @xmath306 with the addition of the bound state contribution ( [ eq : a_u_bound ] ) the sum rule , @xmath307 , is then seen to be satisfied .      when we have two coupled channels there is a spectral density for both the open channel atoms and the closed channel molecules . we start by addressing the latter . the dressed green s function for the closed channel is found by summing all the terms in the dyson series based on the bare propagator ( [ g_cl_approx ] ) . the result is @xmath308 where @xmath309 is the molecular self - energy ( [ sigma_cl ] ) . if @xmath310 is outside the continuum , @xmath311 is real , and the spectral function is only non - zero at the poles of @xmath312 , which coincide with those of the scattering amplitudes , _ i.e. _ the bound states of the coupled channels problem ( [ e_bind ] ) . the spectral weight at the poles is given by the corresponding residues , which are just the probability of finding the bound state pair in the closed channel , @xmath262 : @xmath313 for energies inside the continuum , we obtain the usual structure , @xmath314 ^ 2 + \frac{\hbar^2\gamma^2(e , k)}{4 } } , \label{eq : a_m_continuum}\ ] ] where the energy dependent shift , @xmath315 , and decay rate , @xmath316 , are defined in eqs . ( [ eq : delta ] ) and ( [ eq : gamma ] ) , respectively . in the middle of the band the shift and width are slowly varying functions of energy , and the molecular spectral function has an approximately lorentzian shape with a peak position given approximately by the resonance energy @xmath152 . however , near the band edges significant deviations from this simple form arise due to the strong interference with the continuum threshold . this leads to a double peak structure of @xmath317 , when the bound state is about to enter the continuum from below , as is evident from fig . [ fig : a_m ] , which shows @xmath317 as a function of the energy and @xmath152 for a molecule at rest .    . inside the band @xmath317 is a continuous function given by ( [ eq : a_m_continuum ] ) . outside the band the delta function peaks at the bound state poles of the propagator , @xmath312 , are indicated by vertical lines with a height of @xmath318 . offset below is a contour plot of @xmath317 inside the band , with bright colors indicating high spectral density . the solid lines in this plane are the bound state energies , while @xmath175 is indicated by the dashed line . ] the molecule lifetime also follows from fermi s golden rule , @xmath319 . though the density of states diverges at the edges of the continuum ( van hove singularities ) the decay rate goes to zero at both thresholds , due to vanishing overlap of the background scattering state with the bare resonance state when @xmath320 .      for the open channel the full green s function in the resonant case reads @xmath321 where @xmath322 $ ] is the self - energy for the relative motion of an atom pair in the open channel . this gives the spectral function @xmath323 ^ 2\delta(e - e_b)\ ] ] for energies outside the band . the open channel green s functions has exactly the same poles as the molecular propagator ( [ eq : g_m ] ) , but the residues at the atomic poles depend on the separation between the atoms in the pair . however , using the identity @xmath324 one may show that @xmath325 such that the combined spectral weight of the open channel and closed channel at a pole is unity . we must sum over all possible separations for the open channel weight , since the component of the bound state wavefunction in the open channel is distributed over the lattice according to ( [ eq : psi_b_op ] ) . thus far we have only investigated the dynamics in the lowest bloch band of the lattice potential . in this section we analyze the neglected coupling to higher bands . first we estimate the size of the second order energy shift from the excluded higher bands and derive a criterion for the validity of our theoretical model . next , in sec . [ sec : incl - high - bloch ] we discuss how a modest coupling to higher atomic bands is expected to modify the central conclusions of this work .      in this section we set up a criterion for using the theory , by comparing the second order energy shift from the omitted atomic and molecular exited bloch bands to the relevant energy scales . the dominant contribution to the energy shift comes from the bare resonance state coupling to atom pairs in higher bands , while the shift from open channel atoms coupling to higher molecular bloch bands are several orders of magnitude smaller for two reasons . firstly , the coupling matrix elements between an atom pair in the lowest band and an excited band bare molecule are much smaller than the coupling to a molecule in the lowest band ( see appendix  [ sec : mol - couplings ] ) . secondly , the energy splitting between the ground and excited molecular bands is larger than the atomic band gap , since the molecules experience a deeper lattice potential . hence we approximate the unperturbed wavefunction by @xmath326 , the resonance state in the lowest molecular band . for simplicity we will furthermore assume the width of the bloch bands to be small compared to the separation from the lowest band , so that the energy can be approximated by the average energy of the band .    under these assumptions the second order energy shift is @xmath327 where the sum is over pairs of atoms in the @xmath12th and @xmath1th bloch band , @xmath328 is the inter - band coupling matrix element defined in appendix [ sec : calc - param ] , and @xmath329 is the average energy of the @xmath12th band . if the atoms are identical bosons , we can take the sum to be over @xmath330 to avoid double counting .    we will require the shift to be small compared to the relevant energy scales . if @xmath331 is comparable to the gap between the ground and first excited band @xmath332 , the higher bands will become populated , and the lowest band approximation of our model will break down  @xcite . on the other hand a stronger condition can be derived by requiring that the energy shift be small in comparison with the width of the ground band , @xmath333 , which sets the relevant energy scale for the dynamics in our model . by fitting the respective energies from numerical band structure calculations at different lattice depths we arrive at the condition @xmath334 where @xmath335 when comparing @xmath331 with the band gap , and @xmath336 when requiring the energy shift to be much smaller than the width of the lowest band  @xcite . the scaling factors are the bohr radius , @xmath337 , the bohr magneton , @xmath338 , gauss ( g ) , the atomic mass unit ( @xmath339 ) , and nanometers . if both criteria are met , coupling to higher bands can be safely neglected . otherwise , if the energy shift incurred from the coupling to atoms in higher bloch bands is significant with respect to the band width yet negligible with regards to the band gap , we might expect quantitative changes to our results , but as we discuss below , the qualitative features of our model will remain . in the case where the lowest band approximation fails coupling to higher bands can be systematically accounted for in the tight binding limit  @xcite , leading to an effective single - band hamiltonian with a non - local coupling between dressed molecules and unbound atom pairs  @xcite . both of the above stated conditions can be fulfilled for the feshbach resonance studied by syassen _ _ in @xmath160rb near 414 g  @xcite . in the entrance channel both atoms are in the @xmath340 hyperfine state . the resonance is narrow , with @xmath341  mg , and has a very small difference in magnetic moments between the two channels , @xmath342  khz / g , which is preferable if eq . ( [ eq : validity_condition ] ) is to hold . a final condition for our results to be experimentally relevant is that the variation of the resonance energy due to magnetic field noise @xmath343 is much smaller than the bandwidth @xmath333 . with a lattice wavelength of 830.44 nm  @xcite the recoil energy is @xmath344  khz ( equivalent to @xmath345 ) and for the particular resonance considered here @xmath34632 mg . for a longitudinal lattice depth of @xmath161 we therefore require @xmath347 according to table  [ tab : parameters ] . in ref . @xcite the magnetic field noise was found to be less than 4 mg .      as discussed in the previous subsection the major correction to our model comes from the coupling to higher atomic bands . we now include these couplings into our model . for each atomic band we apply the analysis of sec . [ sec : entr - chann - greens ] to obtain band specific green s functions @xmath348 and @xmath349 for the relative motion of one atom in bloch band @xmath12 and another in band @xmath1 . as for the first band @xmath350 is the green s function for noninteracting particles and @xmath351 is the green s function for particles with an on - site interaction , @xmath41 . we determine @xmath348 and @xmath349 by replacing @xmath94 by @xmath352 $ ] in eq . , and replacing @xmath42 by the band specific interaction strength @xmath353 in eq . ( for the definition of @xmath353 see appendix [ sec : calc - param ] ) . we now introduce the coupling @xmath115 between the first molecular band and all the atomic bands . this gives the molecular self - energy @xmath354 which reduces to ( [ sigma_cl ] ) , when we exclude contributions from higher bands . because we only consider energies well outside the excited bands , @xmath355 is real and slowly varying with energy for @xmath356 . therefore we conclude that the lifetime @xmath316 of the molecules is unchanged and that the change in the resonance shift @xmath357 must vary slowly . a change of the resonance shift by an amount @xmath358 will affect some of the results in secs . [ sec : two - chann - scatt ] , [ sec : bound - states ] and [ sec : resonant - case ] . the changes amount to simply replacing @xmath359 by @xmath360 , a renormalized resonance energy .    for @xmath27 and @xmath157 we obtain for the parameters in @xcite an approximate relationship for the change in resonance shift @xmath361 due to the higher atomic bands  @xcite . for the lattice depth used throughout this paper @xmath358 is only 2% of the bandwidth @xmath333 , and it is reasonable to neglect higher bands .    a direct coupling between ground and excited atomic bands may be included by renormalizing the local interactions to include the effect of virtual transitions to excited levels through band changing collisions . this leads to an energy dependent correction to the coupling parameters  @xcite . but since the correction is inversely proportional to the gap separating the ground and excited bands , it is also expected to be small . we have assumed onsite interactions and only nearest - neighbor tunneling in all atomic bands which is a questionable assumption for the higher bands , but eq . ( [ eq : deltachange ] ) still provides an estimate of the magnitude of the correction to the real part of the molecular self - energy . we have presented a complete description of the feshbach physics of an atom pair in an optical lattice represented by a discrete two - channel hamiltonian . we have demonstrated how the structured continuum due to the periodic potential modifies the resonance line shape and position as well as the properties of the bound states outside the continuum . the lattice potential facilitates repulsively bound pair states in the gap between the two - particle bloch bands , and near a feshbach resonance a further class of exotic bound states may be created , which are only stable when their center of mass momentum exceeds a critical value , thus enabling us to tune the feshbach physics by controlling the molecular motion . ignoring the small modulation of the resonance energy due to tunneling of the bare molecules , the feshbach molecules are stable bound states of the system if their center of mass quasi - momentum in the first brillouin zone satisfies @xmath362 , where @xmath363 if @xmath364 the pair state lies inside the continuum and decays with a rate @xmath162 , which we have determined analytically . for magnetic fields such that @xmath365 the feshbach molecules are stable for all @xmath38 . note that while the decay rate near threshold is proportional to the density of states , which diverges at the band edges , @xmath162 goes continuously to zero as @xmath310 approaches the edge of the continuum , since the overlap between the bare resonance state and the open channel continuum states is suppressed due to the background interactions in the open channel . an additional effect of the background interactions is to shift and modify the resonance profile . with a non - zero onsite interaction @xmath42 the position of the resonance is shifted away from the bare resonance energy @xmath120 , and the transmission of atoms in the open channel is described by a fano profile , with a shape determined by the energy dependent shift and width of the resonance . the energy dependence is most pronounced near the band edges , where the coupling to the continuum is strongest . this gives rise to threshold effects , when the resonance energy is tuned close to the band edge , as evidenced e.g. in the double peak structure of the molecule spectral function . we have outlined how the energy dependence of the scattering can be mapped out spectroscopically by inducing dissociative transitions from molecular bound states to two - particle scattering states inside the band .    from the expansion of the scattering amplitude as the collision energy approaches the band edges we have defined a generalized scattering length , which exhibits a feshbach resonance at both the lower and the upper threshold of the continuum . the diverging scattering length coincides with the appearance / disappearance of the feshbach bound state . in the limit of a deep lattice potential the lowest bloch band narrows , and these two resonances coalesce . since the bare molecules and the entrance channel atoms experience different lattice induced energy offsets , the position of both the upper and the lower threshold feshbach resonance is shifted relative to the position of the feshbach resonance in free space . for the specific resonance we have considered here , this shift exceeds the width of the resonance . for the interaction between the atoms in the open channel we assume that @xmath41 is a contact potential of strength @xmath366 , where @xmath367 is the ( free space ) background scattering length . therefore the on - site interaction between a pair of atoms in bands @xmath12 and @xmath1 becomes @xmath368 ^ 4        \nonumber \\      & & \times\int dx_3 \ , \big [ w_{n0}^\parallel(x_3)w_{m0}^\parallel(x_3 )      \big]^2,\end{aligned}\ ] ] under the assumption that only the lowest band in the transverse lattice is occupied . since all lattice wells are identical , the integrals can be performed over wannier functions centered on any site . the coupling to the closed channel is of short range , since it is proportional to the difference between the atomic triplet and singlet potentials , which have the same long range form  @xcite . in addition , the closed channel bound state has a size ( typically tens of  ) , which is small compared with the lattice spacing ( hundreds of nm ) . we thus assume a purely on - site coupling , which in free space would couple molecular and atomic plane wave states with a strength @xmath369 . here @xmath123 is the difference between the closed and open channel magnetic moments and @xmath370 is the magnetic field width of the feshbach resonance in the absence of the lattice . in the discrete lattice model this results in the matrix element @xmath371 for the coupling of a pair of atoms in bands @xmath12 and @xmath1 to a bare molecule in band @xmath372 . the matrix elements are zero whenever @xmath373 is even due to parity of the wannier functions . a molecule experiences twice as deep a lattice as the atoms , see eq . ( [ v_lat_mol ] ) , and has twice the atomic mass . consequently , separate wannier functions , @xmath374 , have to be calculated for the center of mass motion of the resonance state .    with the exception of sec . [ sec : validity ] we will consider exclusively the lowest bloch band for the atoms and molecules . therefore , we will only need three parameters to characterize the dynamics in the lattice : @xmath375 , @xmath376 and @xmath377 .    as discussed in section  [ sec : validity ] our discrete lattice model is valid for a narrow feshbach resonance . as a particular example we use the resonance in @xmath160rb near 414 g , which has been studied by syassen _ et al . _ @xcite . for this resonance @xmath341 mg and @xmath342 khz / g , while the background scattering length is 100.8 bohr radii . for all the plots in this paper we consider a lattice with @xmath378 and @xmath379 . this results in interaction parameters @xmath380 and @xmath381 . table  [ tab : parameters ] gives the relevant model parameters for this resonance and a range of lattice depths . .discrete lattice parameters for the @xmath160rb feshbach resonance near 414 g. the recoil energy is @xmath344 khz . [ cols="^,^,^,^,^,^",options="header " , ] in free space the coupling between the collision channels induces a shift from the magnetic field @xmath382 , where the resonance state is degenerate with the entrance channel threshold [ @xmath383 , to the position @xmath384 of the scattering length divergence . this offset may be approximately evaluated by @xmath385 here @xmath386 $ ] is a dimensionless function of the ratio @xmath387 between the background ( free space ) scattering length and the so - called mean scattering length @xmath388 , which is a characteristic length scale of the long range @xmath389 part of the molecular potential  @xcite . in the lattice the energy of the closed channel resonance state relative to the zero of energy ( center of the lowest band ) is then given by @xmath390+e_m(k)-e_a,\ ] ] where the lattice induced energy shifts @xmath391 and @xmath392 of the open and closed channel , respectively , are defined in fig . [ fig : energydiagram ] . the molecular wannier states , @xmath393 , where @xmath394 and @xmath395 refers to the molecular bloch band index and lattice site index along direction @xmath91 , respectively , constitute a complete set on which we can expand the atomic function @xmath396\times w_{n0}^\parallel(x_3)w_{m0}^\parallel(x_3)$ ] . if we denote the expansion coefficients @xmath397 , we have @xmath398 and @xmath399 , where @xmath400 is the @xmath401-norm . now a simple norm argument gives the following identity for the neglected molecular couplings to higher bands @xmath402 where @xmath403 for the 414 g resonance in @xmath160rb . for a lattice with @xmath404 and @xmath226 this results in the estimate @xmath405 for an atom pair in the lowest band ( @xmath406 ) . we therefore neglect the coupling of ground band atoms to molecules in higher bands in comparison with the matrix elements , which couple molecules in the lowest bloch band with excited atom pairs . m. ben dahan , e. peik , j. reichel , y. castin , and c. salomon , phys . lett . * 76 * 4508 ( 1996 ) . s. r. wilkinson , c. f. bharucha , k. w. madison , q. niu , and m. g. raizen , phys . rev . lett . * 76 * , 4512 ( 1996 ) . m. greiner , m. o. mandel , t. esslinger , t. hnsch , and i. bloch , nature ( london ) * 415 * , 39 ( 2002 ) . i. bloch , j. dalibard , and w. zwerger , rev . * 80 * , 885 ( 2008 ) . o. mandel , m. greiner , a. widera , t. rom , t. w. hnsch , and i. bloch , nature ( london ) * 425 * , 937 ( 2003 ) . t. khler , k. gral , and p. s. julienne , rev . phys . * 78 * , 1311 ( 2006 ) . g. thalhammer , k. winkler , f. lang , s. schmid , r. grimm , and j. hecker denschlag , phys . . lett . * 96 * , 050402 ( 2006 ) . t. volz , n. syassen , d. m. bauer , e. hansis , s. drr , and g. rempe , nature phys . * 2 * , 692 ( 2006 ) . k. winkler , g. thalhammer , f. lang , r. grimm , j. hecker denschlag , a. j. daley , a. kantian , h. p. bchler , and p. zoller , nature ( london ) * 441 * , 853 ( 2006 ) . c. ospelkaus , s. ospelkaus , l. humbert , p. ernst , k. sengstock , and k. bongs , phys . lett . * 97 * , 120402 ( 2006 ) . p. a. f. da veiga , l. ioriatti , and m. ocarroll , phys . e * 66 * , 016130 ( 2002 ) . s. m. mahajan and a. thyagaraja , j. phys . a * 39 * , l667 ( 2006 ) . j. hecker denschlag and a. j. daley , proceedings of the international school of physics `` enrico fermi '' , course clxiv , ultra - cold fermi gases , e - print arxiv : cond - mat/ 0610393 . m. grupp , r. walser , w. p. schleich , a. muramatsu , and m. weitz , j. phys . b. * 40 * , 2703 ( 2007 ) . n. nygaard , b. i. schneider , and p. s. julienne , phys . a * 73 * , 042705 ( 2006 ) . u. fano , phys . rev . * 124 * , 1866 ( 1961 ) . p. w. anderson , phys . rev . * 124 * , 41 ( 1961 ) . duan , phys . lett . * 95 * , 243202 ( 2005 ) . l. d. carr and m. j. holland , phys . rev . a * 72 * , 031604(r ) ( 2005 ) . f. zhou , phys . b * 72 * , 220501(r ) ( 2005 ) . n. nygaard , r. piil , and k. mlmer , phys . a * 77 * , 021601(r ) ( 2008 ) . s. flling , s. trotzky , p. cheinet , m. feld , r. saers , a. widera , t. mller , and i. bloch , nature ( london ) * 448 * , 1029 ( 2007 ) . n. strohmaier , y. takasu , k. gnter , r. jrdens , m. khl , h. moritz , and t. esslinger , phys . lett * 99 * , 220601 ( 2007 ) . n. syassen , d. m. bauer , m. lettner , d. dietze , t. volz , s. drr , and g. rempe , phys . lett . * 99 * , 033201 ( 2007 ) . j. u. nckel and a. d. stone , phys . b * 50 * , 17415 ( 1994 ) . m. olshanii , phys . lett . * 81 * , 938 ( 1998 ) . t. bergeman , m. g. moore , and m. olshanii , phys . * 91 * , 163201 ( 2003 ) . v. peano , m. thorwart , c. mora , and r. egger , new . j. phys . * 7 * , 192 ( 2005 ) . j. i. kim , j. schmiedmayer , and p. schmelcher , phys . a * 72 * , 042711 ( 2005 ) . v. a. yurovsky , phys . a * 71 * , 012709 ( 2005 ) . v. a. yurovsky , phys . rev . a * 73 * , 052709 ( 2006 ) . r. piil and k. mlmer , phys . a * 76 * , 023607 ( 2007 ) . s. a. gurvitz and y. b. levinson , phys . b * 47 * , 10578 ( 1993 ) . j. eberly , am . j. phys . * 33 * , 771 ( 1965 ) . a. j. moerdijk , b. j. verhaar , and a. axelsson , phys . a * 51 * , 4852 ( 1995 ) . g. orso , l. p. pitaevskii , s. stringari , and m. wouters , phys . * 95 * , 060402 ( 2005 ) . r. micnas , j. ranninger , and s. robaszkiewicz , rev . * 62 * , 113 ( 1990 ) . a. b. kuklov and b. v. svistunov , phys . * 90 * , 100401 ( 2003 ) . duan , e. demler , and m. d. lukin , phys . lett * 91 * , 090402 ( 2003 ) . d. petrosyan , b. schmidt , j. r. anglin , and m. fleischhauer , phys . a * 76 * , 033606 ( 2007 ) . t. miyakawa and p. meystre , j. low . phys * 148 * , 429 ( 2007 ) . m. khl , h. moritz , t. stferle , k. gnter , and t. esslinger , phys . lett . * 94 * , 080403 ( 2005 ) . the fit was performed for @xmath407 and @xmath408 , and the sum in ( [ eq:2nd_order_shift ] ) was taken over the lowest 35 bands . d. b. m. dickerscheid , u. al khawaja , d. van oosten , and h. t. t. c. stoof , phys . a * 71 * , 043604 ( 2005 ) . k. b. gubbels , d. b. m. dickerscheid , and h. t. c. stoof , new . j. phys * 8 * , 151 ( 2006 ) . r. b. diener and t .- l . ho , phys . rev . lett . * 96 * , 010402 ( 2006 ) . s. drr , t. volz , and g. rempe , phys . a * 70 * , 031601(r ) ( 2004 ) . this fit was performed for @xmath409 and @xmath410 , and the contribution from the first 35 bands were included in ( [ eq : sigma_m_higher ] ) .
Abstract:

Answer: we analyze the scattering and bound state physics of a pair of atoms in a one - dimensional optical lattice interacting via a narrow feshbach resonance . the lattice provides a structured continuum allowing for the existence of bound dimer states both below and above the continuum bands , with pairs above the continuum stabilized by either repulsive interactions or their center of mass motion . inside the band the feshbach coupling to a closed channel bound state leads to a fano resonance profile for the transmission , which may be mapped out by rf- or photodissociative spectroscopy . we generalize the scattering length concept to the one - dimensional lattice , where a scattering length may be defined at both the lower and the upper continuum thresholds . as a function of the applied magnetic field the scattering length at either band edge exhibits the usual feshbach divergence when a bound state enters or exits the continuum . near the scattering length divergences the binding energy and wavefunction of the weakly bound dimer state acquires a universal form reminiscent of those of free - space feshbach molecules . we give numerical examples of our analytic results for a specific feshbach resonance , which has been studied experimentally .