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: recent advances in atom de broglie waveguide technology @xcite and its potential applicability to atom interferometry @xcite and integrated atom optics @xcite create a need for accurate theoretical modelling of such systems in the low temperature , tight waveguide regime where transverse excitations are frozen out and the quantum dynamics becomes essentially one - dimensional ( 1d ) ( tonks - gas limit ) . it has been shown by olshanii @xcite , and also recently by petrov et al . @xcite , that at sufficiently low temperatures and densities , high transverse frequencies @xmath0 , and large positive scattering length , where thermal and longitudinal zero - point energies are small compared with @xmath1 , a bose - einstein condensate ( bec ) in a thin cigar - shaped trap has dynamics which reduce to those of a 1d gas of hard core , or impenetrable , point bosons . this is a model for which the exact many - body energy eigensolutions were found in 1960 using an exact mapping from the hilbert space of energy eigenstates of an _ ideal _ gas of spinless fermions to that of many - body eigenstates of hard core , and therefore _ strongly interacting _ , bosons @xcite . in this limit there are strong short - range pair correlations which are omitted in the gross - pitaevskii ( gp ) approximation , which assumes that all @xmath2 bosons occupy the same orbital ( complete bec , condensed fraction unity ) . in the absence of a trap potential it is known @xcite that the occupation of the lowest orbital is of order @xmath3 where @xmath2 is the total number of atoms , in contrast to @xmath2 for the ideal bose gas as well as the gp approximation . nevertheless , this system exhibits some bec - like behavior such as talbot recurrences following an optical lattice pulse @xcite and dark soliton - like behavior in response to a phase - imprinting pulse @xcite . the case of harmonically trapped , hard core bosons in 1d is more relevant to recent atom waveguide experiments @xcite . the spatial profile of the single - particle density is expressible in closed form , and has recently been shown @xcite to be well approximated by a modified 1d effective field theory , although we have recently shown in a numerically accurate time - dependent calculation @xcite that spatial interference fringes of separated and recombined condensates in the exact many - body solution are much weaker than those predicted by the corresponding time - dependent mean field theory @xcite . although the fermi - bose mapping theorem @xcite implies that all physical properties expressible in terms of spatial configurational probabilities are the same for the actual bosonic system and the fictitious  spinless fermion " system used for the mapping , the momentum distribution of the bosonic system , or more generally its occupation distribution over the relevant orbitals for a given geometry , is very different in the bosonic system . it is known @xcite that for a spatially uniform system of hard core bosons in 1d , the momentum distribution is strongly peaked in the neighborhood of zero momentum , whereas that of the corresponding fermi system is merely a filled fermi sea . in the case of hard core bosons in a 1d harmonic trap , it has been an interesting question whether the system undergoes true bec or merely an attenuated one such as that in the uniform system . ketterle and van druten @xcite have shown that true bec occurs for a finite number of atoms in a 1d harmonic oscillator ( ho ) in the case of an _ ideal _ gas , but the behavior turns out to be different for the exact many - body solutions in the tonks ( 1d ) limit . the most fundamental definition of bec and the condensate orbital is based on the large distance behavior of the one - particle reduced density matrix @xmath4 . if off - diagonal long - range order ( odlro ) is present and hence the largest eigenvalue of @xmath5 is macroscopic ( proportional to @xmath2 ) then the system is said to exhibit true bec and the corresponding eigenfunction , the condensate orbital , plays the role of an order parameter @xcite . although the precise definition of odlro requires a thermodynamic limit not strictly applicable to mesoscopic traps , the gp approximation assumes from the start that odlro and macroscopic occupation of a single orbital are good approximations in a trap , so examination of this assumption is important . we have recently @xcite used the fermi - bose mapping theorem to determine , for this case of 1d impenetrable bosons in a harmonic trap , the exact many - body ground state and its salient features , including the one - particle reduced density matrix and its eigenvalues ( occupation number distribution function ) and eigenfunctions ( natural orbitals ) , as well as the momentum distribution function . we find that the largest eigenvalue of the one - particle density matrix is proportional to lower than the first power of @xmath2 , as is the momentum distribution function at @xmath6 , so that the tonks gas does not show true bec . in this paper we give an overview of our recent work on exact many - body quantum solutions for 1d tonks gases . in particular , secs . ii and iii describe the ground state and dynamic properties of tonks gases , respectively . for the most part we investigate 1d tonks gases assuming the conditions are satisfied for realizing them . in sec . iv we examine the use of bessel optical dipole traps for experimental realization of 1d tonks gases . we first describe the fermi - bose mapping method for obtaining exact solutions of the time - independent many - body schrdinger equation in the tonks limit ( impenetrable point bosons ) and its application to the spatially uniform case ( no trap potential ) . this description is very brief and the original literature @xcite should be consulted for details . some salient features of the recent generalization to the harmonically trapped tonks gas @xcite are then described . our basic model consists of @xmath2 bosonic atoms at zero temperature moving in 1d ( @xmath7 ) , the quantum motion in the other two spatial dimensions having been frozen out by tight transverse confinement via an atomic waveguide . the schrdinger hamiltonian is then assumed to have the structure : @xmath8 where @xmath9 is the one - dimensional position of the @xmath10 particle ( atom ) and @xmath11 is symmetric ( invariant ) under permutations of the particles . the two - particle interaction potential is assumed to contain a hard core of 1d diameter @xmath12 . this is conveniently treated as a constraint on allowed wave functions @xmath13 : @xmath14 rather than as an infinite contribution to @xmath11 , which then consists of all other ( finite ) interactions and external potentials . let @xmath15 be a fermionic solution of @xmath16 which is antisymmetric under all particle pair exchanges @xmath17 , hence all permutations . one can consider @xmath18 to be either the wave function of a fictitious system of `` spinless fermions '' , or else that of a system of real fermions whose spins are all aligned , as in magnetically trapped atomic vapor becs . define a  unit antisymmetric function " @xcite @xmath19 where @xmath20 is the algebraic sign of the coordinate difference @xmath21 , i.e. , it is + 1(-1 ) if @xmath22(@xmath23 ) . for given antisymmetric @xmath18 , define a bosonic wave function @xmath24 by @xmath25 which defines the fermi - bose mapping . @xmath24 satisfies the hard core constraint ( 2 ) if @xmath18 does , is totally symmetric ( bosonic ) under permutations , obeys the same boundary conditions as @xmath18 , and @xmath26 follows from @xmath27 @xcite . in the case of periodic boundary conditions ( no trap potential , spatially uniform system ) one must add the proviso that the boundary conditions are only preserved under the mapping if @xmath2 is odd , but the case of even @xmath2 is easily accomodated by imposing periodic boundary conditions on @xmath18 but _ anti_periodic boundary conditions on @xmath24 . for a trapped system the boundary condition that wave functions vanish at infinity is preserved for all @xmath2 . the mapping theorem leads to explicit expressions for all many - body energy eigenstates and eigenvalues of a 1d scalar condensate ( bosons all of the same spin ) under the assumption that the only two - particle interaction is a zero - range hard core repulsion , represented by the @xmath28 limit of the hard - core constraint . such solutions were obtained in sec . 3 of the original work @xcite for periodic boundary conditions and no external potential . such a situation could be realized experimentally , for example , by trapping the tonks gas around the peak of an tight toroidal optical dipole trap formed using a red - detuned laguerre - gaussian laser beam @xcite . in a system of 1d bosons with a hard - sphere interaction of diameter @xmath12 , it is sufficient at low densities @xcite to consider the case of impenetrable point particles , the zero - range limit @xmath28 . since wave functions of  spinless fermions " are antisymmetric under coordinate exchanges , their wave functions vanish automatically whenever any @xmath29 , the constraint has no effect , and the corresponding fermionic ground state is the ground state of the _ ideal _ gas of fermions , a slater determinant of the lowest @xmath2 single - particle plane - wave orbitals . the exact many body ground state was found to be a pair product of bijl - jastrow form @xmath30| , \ ] ] where @xmath31 is the perimeter of the annular trap . in spite of the very long range of the individual pair correlation factors @xmath32|$ ] , the pair distribution function @xmath33 , which is physically the joint probability density that if one particle is found at @xmath34 a second will be found at @xmath9 , was found to be of short range @xmath35 with @xmath36 the linear number density of the system . clearly , @xmath37 which reflects the hard core nature of the two - particle interaction . by examination of the excited states the system was also found to support propagation of sound with speed @xmath38 @xcite . here we briefly review our recently - obtained exact solution for the many - body ground state of a harmonically trapped tonks gas @xcite . the hamiltonian of n bosons in a 1d harmonic trap is @xmath39   .\ ] ] we again assume that the two - body interaction potential consists only of a hard core of negligible diameter @xmath40 . it follows from the mapping theorem that the exact n - boson ground state @xmath41 is @xmath42 where @xmath43 is the ground state of a fictitious system of @xmath2 spinless fermions with the same hamiltonian and constraint . since the hard - core interaction has no effect in the zero - range limit for fermions , whose wave functions already vanish at contact due to antisymmetry , it follows that the fermionic ground state is a slater determinant of the lowest @xmath2 single - particle eigenfunctions @xmath44 of the harmonic oscillator ( ho ) : @xmath45 the ho orbitals are @xmath46 with @xmath47 the hermite polynomials and @xmath48 being the ground state width of the harmonic trap for a single atom . by factoring the gaussians out of the determinant and carrying out elementary row and column operations , one can cancel all terms in each @xmath49 except the one of highest degree , with the result @xcite @xmath50 substitution into ( 6 ) then yields a simple but exact analytical expression of bijl - jastrow pair product form for the @xmath2-boson ground state : @xmath51 \prod_{1\le j < k\le n}|x_{k}-x_{j}|\ ] ] with @xmath52 and normalization constant @xmath53^{-1/2 }   .\ ] ] it is interesting to note the strong similarity between this exact 1d @xmath2-boson wave function and the famous laughlin variational wave function of the 2d ground state for the quantized fractional hall effect @xcite , as well as the closely - related wave functions for bosons with weak repulsive delta - function interactions in a harmonic trap in 2d found recently by smith and wilkin @xcite . both the single particle density and pair distribution function depend only on the absolute square of the many - body wave function , and since @xmath54 they reduce to standard ideal fermi gas expressions . the single particle density , normalized to @xmath2 , is @xmath55 we shall not exhibit it here since it has recently been calculated by kolomeisky _ _ @xcite ; see also our recent discussion of the time - dependent case @xcite . the pair distribution function , normalized to @xmath56 , is @xmath57 noting that terms with @xmath58 , which vanish by antisymmetry , can be formally added to the summation ( 13 ) , one can write @xmath59 although the hermite polynomials have disappeared from the expression ( 10 ) for the many - body wave function , they reappear upon integrating @xmath60 over @xmath61 coordinates to get the single particle density @xmath62 and over @xmath63 to get the pair distribution function @xmath64 , and the expressions in terms of the ho orbitals @xmath65 are the most convenient for evaluation . figure [ fig : one ] shows a gray scale plot of the dimensionless pair distribution function @xmath66 versus the normalized coordinates @xmath67 for a ) @xmath68 , b ) @xmath69 , and c ) @xmath70 . some qualitative features of the pair distribution function are apparent : in the first place it follows either from the original expression ( 12 ) or from eq . ( 14 ) that @xmath71 vanishes at contact @xmath72 , as it must because of impenetrability of the particles , and we see this to be true in fig . [ fig : one ] . furthermore , the correlation term @xmath73 is a truncated closure sum and approaches the dirac delta function @xmath74 as @xmath75 , as is to be expected since the healing length in a spatially uniform 1d hard core bose gas varies inversely with particle number @xcite . as a result the width of the null around the diagonal @xmath76 decreases with increasing @xmath2 , and vanishes in the limit . away from the diagonal along @xmath77 the pair distribution function rises , exhibits modulations for @xmath78 , due to the oscillatory nature of the ho orbitals , before decreasing back to zero at large distances . for @xmath79 much larger than the healing length , @xmath80 reduces to the uncorrelated density product @xmath81 , so the spatial extent of the pair distribution function is that of the density and varies as @xmath82 @xcite . the reduced single - particle density matrix with normalization @xmath83 is @xmath84 although this multi - dimensional integral can not be evaluated analytically , it can be evaluated numerically by monte carlo integration for not too large values of @xmath2 ( the computing time scales as @xmath85 ) . figure [ fig : two ] shows a gray scale plot of the dimensionless reduced single - particle density matrix @xmath86 versus the normalized coordinates @xmath87 and @xmath88 for a ) @xmath68 , b ) @xmath69 , and c ) @xmath70 . along the diagonal @xmath89 reproduces the single - particle density @xcite . the off - diagonal elements relate to odlro , and it is clear that as @xmath2 increases the off - diagonal elements are decreasing , in contrast with the diagonal . this is an indication that odlro vanishes for a system of hard core bosons in a 1d ho in the thermodynamic limit , although at present only numerical evidence exists , there being no analytical proof generalizing the result of lenard for the untrapped tonks gas @xcite .    in a macroscopic system , the presence or absence of bec is determined by the behavior of @xmath4 as @xmath90 . off - diagonal long - range order is present if the largest eigenvalue of @xmath5 is macroscopic ( proportional to @xmath2 ) , in which case the system exhibits bec and the corresponding eigenfunction , the condensate orbital , plays the role of an order parameter @xcite . although this criterion is not strictly applicable to mesoscopic systems , if the largest eigenvalue of @xmath5 is much larger than one then it is reasonable to expect that the system will exhibit some bec - like coherence effects . thus we examine here the spectrum of eigenvalues @xmath91 and associated eigenfunctions @xmath92 ( `` natural orbitals '' ) of @xmath5 . although natural orbitals are a much - used tool in theoretical chemistry , they have only recently been applied to mesoscopic atomic condensates @xcite . the relevant eigensystem equation is @xmath93 @xmath91 represents the occupation of the orbital @xmath94 , and one has @xmath95 . numerical evaluation of the integral by discretization yields a matrix eigensystem equation giving accurate numerical results for the largest eigenvalues and associated eigenvectors . we remark that for the corresponding problem of @xmath2 free fermions in a 1d ho , in which case @xmath96 is evaluated using the fermion wave function @xmath97 , the natural orbitals are simply the ho orbitals , and @xmath98 , all higher eigenvalues being zero ( filled fermi sea ) . however , the n - boson wave function is the modulus of the fermion wave function , and this leads to significant differences in the spectrum of natural orbitals and eigenvalues for the hard core bose gas . in fig . [ fig : three](a ) we show a log - log plot of the fractional occupation of the lowest orbital @xmath99 versus the total particle number @xmath2 ( solid line ) , along with a best fit power - law @xmath100 ( dashed line ) . this is to be contrasted with the case of a spatially uniform system of hard core bosons for which @xmath101 @xcite . in both cases the fractional occupation decreases with increasing @xmath2 , and thus does not correspond to a true condensate for which @xmath102 . nevertheless , the occupation of the lowest orbital may still be large , @xmath103 , and is larger than the spatially uniform case @xmath104 , so macroscopic quantum coherence effects reminiscent of bec can still result @xcite . figure [ fig : three](b ) shows the distribution of occupations @xmath91 versus orbital number @xmath105 ( the orbitals are ordered according to eigenvalue magnitude , the largest eigenvalue being @xmath106 ) for @xmath68 ( circles ) , @xmath69 ( stars ) , and @xmath70 ( squares ) . this figure shows that as the lowest orbital occupation @xmath107 increases with increasing @xmath2 so does the range of significantly occupied higher - order orbitals with @xmath108 . this means that the dominance of the lowest orbital decreases with increasing @xmath2 , so singling out @xmath109 as a macroscopic wave function for the whole system becomes more problematic with increasing @xmath2 @xcite . next we examine the momentum distribution for the trapped 1d tonks gas . for a spatially uniform system ( no trap ) the natural orbitals are plane waves , so the occupation distribution of the natural orbitals is the same as the momentum distribution . although this is not the case here due to the effect of the harmonic trap potential , the momentum distribution is still physically important . in terms of the boson annihilation and creation operators in position representation ( quantized bose field operators ) the one - particle reduced density matrix is @xmath110 the momentum distribution function @xmath111 , normalized to @xmath112 , is @xmath113 where @xmath114 is the annihilation operator for a boson with momentum @xmath115 . then @xmath116 the spectral representation of the density matrix then leads to @xmath117 where the @xmath118 are fourier transforms of the natural orbitals : @xmath119 . figure [ fig : four ] shows the numerically calculated dimensionless momentum spectrum @xmath120 versus normalized momentum @xmath121 , with @xmath122 , for a ) @xmath68 , b ) @xmath69 , and c ) @xmath70 . the key features are that the momentum spectrum maintains the sharp peaked structure reminiscent of the spatially uniform case @xcite for the 1d ho , and that the peak becomes sharper with increasing atom number @xmath2 . this is to be expected since as the number of atoms increase the many - body repulsion causes the system to become more spatially uniform within the trap interior . by way of contrast , for a 1d fermi gas the corresponding momentum spectrum is a filled fermi - sea , here rounded out by the trapping potential . in particular , the momentum distribution for the 1d trapped fermi gas can be expressed as @xmath123 thus , the momentum spectrum provides a means of distinguishing between the 1d fermi and tonks gases . figure [ fig : five ] shows an example of the tonks ( dashed line ) and fermi gas ( solid line ) momentum spectra for @xmath70 atoms in a harmonic trap , and the difference is clearly seen . furthermore , the distinction only grows larger with increasing number of atoms , as the filled fermi - sea broadens whereas the tonks spectrum becomes narrower . in a recent paper we have devised a scheme to measure the momentum spectrum of trapped gases based on raman outcoupling , and we refer the reader to our paper for details @xcite . the fermi - bose mapping theorem is very easily generalized so as to yield exact solutions of the time - dependent many - body schrdinger equation ( tdmbse ) by noting that since the mapping function @xmath124 of eq . ( 3 ) is independent of time , one can merely replace @xmath125 by @xmath126 , implying that the relationship of eq . ( 4 ) , with time arguments added to the bose and fermi wave functions , is valid for solutions of the tdmbse . in the special case where the only interatomic interaction is that of hard cores of vanishing diameter , but external potentials including a trap potential as well as time - dependent fields may be present , the many - fermion solutions of the tdmbse are slater determinants of solutions of the _ single particle _ tdse in the given external potential , and each many - boson solution is obtained by multiplying the corresponding determinant by the mapping function @xmath124 . some salient features of our recently - obtained dynamical solutions for solitons in a ring geometry and for interference effects in a dynamically split , harmonically trapped tonks gas will be described . the description will again be brief since details are available in the literature @xcite . dark and gray solitons are a generic feature of the nonlinear schrdinger equation with repulsive interactions , and several calculations of their dynamics based on the mean - field gross - pitaevskii ( gp ) equation have appeared @xcite , as well as experiments demonstrating their existence in atomic becs @xcite . since the underlying _ many - body _ schrdinger equation is linear , this raises the question of how observed solitonic behavior arises . here this issue will be examined with the aid of exact many - body solutions for the tonks gas . the model consists of a 1d hard - core bose gas in a toroidal trap , or ring , with cross section so small that motion is essentially circumferencial . the fermi - bose mapping is employed to generate exact solutions for this problem . we identify stationary solutions which reflect some properties of dark solitons from the gp theory when the ring is pierced at a point by an intense blue - detuned laser . we also present dynamical solutions when half of an initially homogeneous ring bec is phase - imprinted via the light - shift potential of an applied laser , leading to gray soliton - like structures whose velocity depends on the imposed phase - shift @xcite . such structures are apparent for times less than the echo time @xmath127 , with @xmath31 the ring circumference and @xmath128 the speed of sound in the bec . on longer time scales the dynamics becomes very complex showing talbot recurrences which are beyond the gp theory . consider @xmath2 bosons in a tight toroidal trap , for example a toroidal optical dipole potential @xcite , and denote their 1d positions measured around the circumference by @xmath9 . this is equivalent to the exactly - solved model @xcite of @xmath2 impenetrable point bosons in 1d with wave functions satisfying periodic boundary conditions with period @xmath31 equal to the torus circumference , and the fundamental periodicity cell may be chosen as @xmath129 . however , the rotationally invariant quantum states of this problem do not reveal any dark soliton - like structures . to proceed we therefore consider the case that a blue - detuned laser field pierces the ring at @xmath130 by virtue of the associated repulsive dipole force : the light sheet then provides a reference position for the null of the dark soliton . assume that the light sheet is so intense and narrow that it may be replaced by a constraint that the many - body wave function ( hence the orbitals @xmath131 ) must vanish whenever any @xmath132 . then the appropriate orbitals @xmath133 are free - particle energy eigenstates vanishing at @xmath130 and periodic with period @xmath31 . the complete orthonormal set of even - parity eigenstates @xmath134 and odd - parity eigenstates @xmath135 are @xmath136   , \nonumber\\ \phi_{n}^{(-)}(x ) & = & \sqrt{2/l}\sin(2n\pi x / l )        , \end{aligned}\ ] ] with @xmath137 running from @xmath138 to @xmath139 . the odd eigenstates are the same as those of free particles with no @xmath130 constraint , since these already vanish at @xmath130 . however , the even ones are strongly affected by the constraint , their cusp at @xmath130 being a result of the impenetrable light sheet at that point . if one bends a 1d box @xmath140 with impenetrable walls into a ring , identifying the walls at @xmath141 , then those particle - in - a - box eigenfunctions which are even about the box center become identical with the @xmath134 , and their cusp results from the nonzero slope of these functions at the walls . the @xmath2-fermion ground state is obtained by inserting the lowest @xmath2 of these orbitals into a slater determinant ( filled fermi sea ) . since @xmath142 , the one - particle density @xmath62 of the corresponding many - boson ground state given by eq . ( 4 ) is the same as that of the @xmath2-fermion ground state . in the thermodynamic limit @xmath75 , @xmath143 , @xmath144 for fixed @xmath7 , @xmath145 , one finds @xcite @xmath146   .\ ] ] @xmath62 vanishes at @xmath130 and approaches the mean density @xmath147 over a healing length @xmath148 with damped spatial oscillations about its limiting value . suppose next that the light - sheet is turned off at @xmath149 by removing the constraint that the wave function vanish at @xmath130 . the solution of the tdmbsb for the many - boson system is then still given by the mapping theorem , but the slater determinant representing the corresponding many - fermion state has to be built from solutions of the time - dependent single - particle schrdinger equation satisfying the initial conditions of eq . ( [ lightsheet ] ) . the odd - parity solutions are stationary in time since they already vanish at the position of the light sheet . the even - parity solutions are nontrivial , but are expressible as sums over the space - periodic , even - parity solutions of the time - independent free - particle schrdinger equation @xcite : @xmath150 where @xmath151 and @xmath152 . the time - dependent density @xmath153 is then the sum of absolute squares of all @xmath2 orbitals in the fermi sea . it is found @xcite that there are two important time scales : one is the poincar recurrence time @xmath154 . noting that @xmath155 is proportional to @xmath156 , one finds that all terms in the sum are time - periodic with period @xmath157 , which is therefore the recurrence time for the density and in fact all properties of our model @xcite . the other important time is the echo time @xmath158 , the time for sound to make one circuit around the torus . recalling that the speed of sound in this system is @xmath38 @xcite , one finds @xmath159 . for @xmath160 after the constraint is removed , the initial density develops sound waves that propagate around the ring , and that we examine below in the context of phase - imprinting . for @xmath161 the evolution is very complex , but complete recurrences occur for times @xmath162 with fractional revivals in between . suppose next that the toroidal impenetrable point bose gas is in its ground state for times @xmath163 ( no light sheet obstacle in this case ) , and a phase - imprinting laser pulse is applied as a delta - function pulse over half the ring at @xmath149 . this is described by the hamiltonian @xmath164\ ] ] where @xmath165 , i.e. , unity for @xmath166 and zero elsewhere . this is the technique used in recent experiments @xcite , here idealized to a delta - function in time and to sharp spatial edges . before the pulse the most convenient free - particle orbitals in ( 5 ) are plane waves @xmath167 where @xmath168 and @xmath169 with @xmath170 . let @xmath171 be the solution of the corresponding single - particle tdse reducing to the above @xmath172 just before the pulse . then the solutions just after the pulse are @xmath173 . the potential gradients at the pulse edges impart momentum kicks to the particles there which induce both compressional waves propagating at the speed , @xmath128 , of sound and density dips ( gray solitons ) moving at speeds @xmath174 . the expansion of @xmath171 in terms of the unperturbed plane waves is evaluated as @xmath175 and the time - dependent density is the sum of the absolute squares of the lowest @xmath2 of these . figure [ fig : six ] shows numerical simulations for @xmath176 , @xmath177 , and @xmath178 ( solid line ) , and @xmath179 ( dashed line ) : due to symmetry we show only half of the ring @xmath180 , the phase shift being imposed at @xmath181 . considering times short compared to the echo time means that the corresponding results are not very sensitive to the periodic boundary conditions , and also therefore apply to a linear geometry . the initial density profile is flat with a value @xmath182 . for both phase shifts two distinct maxima are seen , which travel at close to the speed of sound @xmath128 , and two distinct minima , which are analogous to gray solitons and travel at velocities @xmath183 . in addition , there are also high wavevector oscillations which radiate at velocities greater than @xmath128 , analogous to precursors in electromagnetic wave propagation in a medium . in the case of a phase shift @xmath178 , the density is symmetric about @xmath181 , whereas for a phase - shift other than a multiple of @xmath184 the evolution is not symmetric , see the dashed line where the global minimum moves to the right in reponse to the phase shift . in fig . [ fig : seven ] we plot the calculated velocity of the global density minimum relative to the speed of sound for a variety of phase shifts @xmath185 . the basic trend is that larger phase - shift means lower velocity , in qualitative agreement with recent experiments @xcite , but there is a sharp velocity peak at @xmath186 : this peak results from the cross - over between two local minima in the density . these general features , the generation of gray solitons and density waves , agree with those of the gp theory , but here arise out of the exact many - body calculation .    a detailed comparison between our results and current experiments is not possible as they do not conform to the conditions for a 1d system . however , some estimates are in order to set the appropriate time scales : if we consider @xmath187rb with a ring of circumference @xmath188 @xmath189 m , and a high transverse trapping frequency @xmath190 hz , then we are limited to atom numbers @xmath191 @xcite , so these are small condensates . we then obtain @xmath192 s , and @xmath193 ms for @xmath176 . finally , we remark that since our approach relied on the mapping between the strongly - interacting bose system and a non - interacting  spinless fermi gas " model , this suggests that dark and gray solitons should also manifest themselves in the density for the 1d fermi system . although real fermions have spin , the interactions used here to generate solitons were spin - independent .      mean - field theory ( mft ) has proven remarkably successful at predicting both the static and dynamic behavior of bose - einstein condensates ( becs ) in weakly - interacting atomic vapors @xcite , including the ground state properties @xcite , the spectrum of collective excitations @xcite , four - wave mixing @xcite , matter - wave solitons @xcite , and interference between becs @xcite . the basic notion underlying mft is of a macroscopic wave function @xcite , or order parameter , which defines the spatial mode into which a significant fraction of the atoms condense below the critical temperature . the macroscopic wave function typically obeys a nonlinear schrdinger equation ( nlse ) , the gross - pitaevskii equation , and is most suitable for dilute bose gases . however , the success of mft is not assured in all cases . for example , in one - dimension ( 1d ) a spatially homogeneous ideal gas in its many - body ground state exhibits complete bec into the lowest single - particle state , but no bec at any nonzero temperature . furthermore , previous exact analysis @xcite of a spatially uniform tonks gas by one of us ( mg ) and its extension by lenard @xcite and vaidya and tracey @xcite have shown that in the many - body ground state the occupation of the lowest single - particle state is of order @xmath3 where @xmath2 is the total number of atoms , in contrast to @xmath2 for usual bec , and similar behavior , but with a slightly larger exponent , occurs in the harmonically trapped tonks gas @xcite . we have recently shown that mean - field theory breaks down in the analysis of the dynamics of 1d atom clouds , in that it predicts interference effects that are absent in the exact theory @xcite . _ @xcite have proposed a nonlinear schrdinger equation ( nlse ) with a quartic nonlinearity to extend the usual mean - field theory for 1d atom clouds . for a harmonic trap the ground - state density profiles from their theory show excellent agreement with the exact many - body results ( see fig . 1 of their paper ) . the key question , then , is whether this extended nlse can be used in all circumstances . to address this issue we examine the problem of a 1d atomic cloud in the ground state of a harmonic trap that is split by a blue - detuned laser and recombined , both using an exact many - body treatment based on the fermi - bose mapping and the approximate nlse : the nlse predicts interference whereas the exact analysis does not . consider a tonks gas which is initially in its @xmath2-body ground state in a harmonic trap @xcite . a central gaussian repulsive potential simulating a blue - detuned laser field is turned on quasi - adiabatically at time @xmath149 to split the initial state . after some time @xmath194 both the harmonic trap and repulsive potential are turned off and the two split components allowed to recombine : this is an interference experiment of the cool , cut , interfere variety previously discussed @xcite . the external potential is taken of the specific form @xmath195 and @xmath196 for @xmath197 . figure [ fig : eight ] shows an illustrative example for @xmath70 with a repulsive potential of height @xmath198 , and width @xmath199 , with @xmath200 the ground state harmonic oscillator width , and @xmath201 . this figure shows a gray scale plot of the single - particle density @xmath202 of the @xmath2-boson system as functions of @xmath203 ( horizonatal axis ) and position @xmath204 ( vertical axis ) , with white being the highest density . the potential height was chosen such that @xmath205 , where @xmath189 is the chemical potential of the @xmath2-particle oscillator ground state @xcite , noting that the top of the @xmath2-particle fermi sea is at @xmath206 . as expected , as the repulsive potential turns on it splits the initial ground state into two separated components . upon release at @xmath207 the two components expand and subsequently recombine . what is noteworthy is that although there is some modulation upon recombination there are no strong interference fringes indicative of the interference provisionally expected for bosons : this was a generic finding from our simulations irrespective of the time scale on which the repulsive potential was turned on @xcite . in contrast , the density @xmath208 for each individual orbital @xmath209 can show strong fringes , but with a different period in each case . thus , the minimal inteference is a result of washing out of the individual interferences by summing over n orbitals . thus , the remnant of any interference fringes decreases with increasing @xmath2 and vanishes in the thermodynamic limit . physically , it makes sense that the interference fringes are all but absent since the fermi - bose mapping shows that in this 1d limit the system of bosons acts effectively like a system of free fermions insofar as effects expressible only in terms of @xmath210 are concerned , so interference is not expected @xcite . the lack of interference is therefore a signature of the fermi - bose duality that occurs in 1d systems of impenetrable particles @xcite . we next turn to the mean - field description proposed by kolomeisky _ @xcite for low - dimensional systems . in particular , they introduce an order parameter @xmath211 , normalized to the number of particles @xmath2 , for such systems . using energy functional arguments they deduce the following nlse with quartic nonlinearity for a 1d system of impenetrable bosons : @xmath212\phi   .\ ] ] our goal is to compare the predictions of this nlse for the same cool , cut , interfere simulation in figure [ fig : eight ] , with the initial condition @xmath213 corresponding to the exact many - body solution , all other parameters being equal . figure [ fig : nine ] shows the corresponding gray - scale plot of the density @xmath214 , and two features are apparent : first , during the splitting phase when the repulsive potential is on there is very good overall agreement between the exact many - body theory and the nlse prediction . second , when the split components are released and recombine they produce pronounced interference fringes , in contrast to the exact theory . indeed , this interference in the mft is to be expected on the basis of previous theoretical work @xcite , even though a quartic ( rather than quadratic ) nonlinearity is employed here . thus the mft can not accurately capture the time - dependent dynamics in all situations . so far we have assumed that conditions are satisfied for the tonks limit . however , previous analysis by olshanii @xcite and petrov et al . @xcite show that there are stringent requirements on the temperature , transverse confinement , and linear atom number density to approach the tonks limit . clearly , what is needed first of all it to create a thin cigar shaped atomic trap so that the transverse mode becomes frozen and the atomic motion becomes essentially 1d . recent experimental developments suggest that tonks gases may be realizable in magnetic atom waveguides @xcite , and bongs et al . @xcite have proposed a hybrid trap composed of optical dipole trap formed with a first - order lg beam combined with magnetic longitudinal trapping . here we examine the use of bessel optical dipole traps as an all - optical means to realize 1d tonks gases . ideal bessel beams are solutions of the free - space wave equation which propagate with unchanging beam profile along the propagation axis which we take as @xmath7 in cylindrical coordinates @xmath215 . the electric field of a monochromatic , linearly polarized ideal bessel beam of order @xmath216 and frequency @xmath217 is @xcite @xmath218 where @xmath219 is a scale electric field value , @xmath220 is the @xmath216th order bessel function , @xmath221 is the azimuthal mode number which we take as positive for simplicity in notation , and @xmath222 and @xmath223 are the radial and longitudinal wavevectors such that @xmath224 with @xmath225 . the zeroth - order solution @xmath226 has a central maximum surrounded by concentric rings of roughly equal power while the higher - order solutions @xmath220 have zero on - axis intensity also with concentric rings . the ultility of bessel beams for optical dipole traps lies in the fact that they can produce very elongated , and hence low aspect ratio quasi-1d traps . in practice , of course , one can not produce an ideal bessel optical beams as it carries infinite power . however , using axiconal optics finite power approximations to bessel beams can be produced which can propagate over significant distances . the basic scheme is illustrated in figure [ fig : ten ] where an incident laguerre - gaussian beam of azimuthal mode index @xmath216 is transformed into a bessel beam of the same order . the interested reader is referred to our recent paper @xcite for details , but for an incident gaussian beam on the axiconal optics the resulting zero - order bessel beam intensity profile past the axicon is @xmath227 where @xmath228 is the input spot size , @xmath229 the input power , @xmath230 , @xmath137 and @xmath231 are the index and opening angle of the axicon , and @xmath232 which gives the longitudinal extent of the bessel beam . this intensity profile of the bessel optical beam is readily converted to an effective dipole potential @xmath233 with @xmath234 the laser detuning from the optical transition frequency @xmath235 , @xmath236 the natural linewidth of the optical transition , @xmath237 is the resonant saturation intensity , and @xmath238 . for a red - detuned laser the potential is negative and the atoms are attracted to the regions of high intensity , whereas for a blue - detuned laser the atoms are repelled into the low field regions . here for a red - detuned laser we have approximated the bessel optical dipole potential by a parabolic potential near the axis under conditions of tight confinement , where @xmath239 is the longitudinal position of the peak of the bessel beam , and @xmath240 a red - detuned ( @xmath241 ) @xmath226 optical dipole potential therefore provides confinement in both the radial and longitudinal directions . here @xmath242 is the radial oscillation frequency with corresponding ground state oscillator width @xmath243 , and @xmath244 is the ratio between the longitudinal and radial trap frequencies @xmath245 @xcite , which also determines the aspect ratio between the radial and longitudinal ground state widths @xmath246 ( in the absence of many - body repulsion ) . petrov _ et al . _ @xcite have theoretically studied the diagram of state for a one - dimensional gas of trapped bosons , assuming @xmath247 , and found that a true bec , or at least a quasi - condensate , with concomitant macroscopic occupation of a single state , is only attained for high enough particle numbers @xmath248 with @xmath249 for @xmath250 and temperatures @xmath251 , one obtains a tonks gas of impenetrable bosons for which hard core repulsion between the bosonic atoms prevents them from penetrating through each other in the one - dimensional system , and the system acquires properties reminiscent of a one - dimensional system of fermionic atoms . a highly elongated bessel beam discussed , say with parameters @xmath252 nm , @xmath253 w , @xmath254 cm , @xmath255 @xmath189 m , would be an ideal candidate for the experimental realization of a tonks gas . the low aspect ratio @xmath256 and tight radial confinement @xmath257 nm result in a high upper boundary @xmath258 for the particle number of the tonks gas . for the commonly used @xmath187rb isotope with a scattering length @xmath259 nm one finds @xmath260 . although this is still a fairly low value it should be possible to experimentally realize a small @xmath187rb tonks gas . however , more promising would be the use of the @xmath261rb isotope , where a feshbach resonance can be used to tune the normally negative scattering length to positive values of several hundred nanometers magnitude @xcite . as @xmath258 is proportional to the square of the scattering length even a moderate increase to @xmath262 nm would make it possible to create a larger tonks gas , with say @xmath263 atoms , which should be easily detectable . the bessel beam trap offers some advantages compared to alternative suggested approaches using magnetic waveguides @xcite and a hybrid magnetic - optical trap @xcite . firstly , it involves only a very simple _ all - optical _ system for which the aspect ratio of the trap may be controlled simply by varying the gaussian spot size incident on the axicon . more specifically , being all - optical , it does not involve material surfaces as in magnetic waveguides , which can cause matter - wave decoherence @xcite . furthermore , it allows for the possibility of trapping multiple magnetic sublevels and the investigation of multi - component tonks gases , which would not be possible in the hybrid magnetic - optical trap of bongs _ et al . in this paper we have given an overview of our work on the quantum dynamics of 1d trapped tonks gases . in particular , we hope to have shown that tonks gases display a rich variety of behaviors that are worthy of experimental investigation . the fermi - bose mapping approach to solving for the 1d tonks shows that many features are shared between the tonks gas of impenetrable bosons and a 1d gas of non - interacting fermions : they have the same ground state energies and density profiles , and neither shows true bec . however , the tonks and fermi gases differ markedly in their momentum spectra , and we have devised a scheme to measure this difference @xcite . furthermore , we have shown that the 1d tonks gas can show dark soliton - like solutions which are typically associated with solutions of the mean - field gross - pitaevskii equation . this is so in spite of the fact that mean - field theory greatly overestimates the coherence properties of the tonks gas , as we demonstrated in our simulations of split and recombined tonks gases . our work presented here shows that tonks gases have some degree of coherence but much less than a true bec , and this may have implications for proposed atomic interferometers employing high density , eg . atom laser , sources . in particular , strong transverse confinement is desirable in atomic interferometers to avoid multi - transverse mode effects , and at the same time one would like to keep the atomic density down to avoid many - body shifts , but these are exactly the requirements for the tonks limit ! therefore it is very important to study the bec - tonks transition @xcite and its effect on the performance of atom interferometers . we hope that our work is a first step in this direction . this work was supported at the university of arizona by office of naval research grant n00014 - 99 - 1 - 0806 and also the us army research office . we thank joe triscari for collaborations on the ground state properties of the tonks gas , and the work on bessel optical dipole traps was conducted in collaboration with j. arlt and k. dholakia from st . andrews university , scotland , and j. soneson at the university of arizona . et al . _ , . d. mller _ et al . n.h . dekker _ et al . _ , . thywissen , r.m . westervelt , and m. prentiss , . hinds , m.g . boshier , and i.g . hughes , . , ed . p.r . berman ( academic press , boston , 1997 ) . j. schmiedmayer , . m. olshanii , . d. s. petrov , g. v. shlyapnikov , and j. t. m. walraven , . m. girardeau , . girardeau , . see particularly secs . 2 , 3 , and 6 . a. lenard , . a. g. rojo , g. l. cohen , and p. r. berman , . m.d . girardeau and e.m . k. bongs _ et al . e. b. kolomeisky _ et al . _ . m.d . girardeau and e.m . vaidya and c.a . w. ketterle and n. j. van druten , . o. penrose and l. onsager , . girardeau , e.m . wright , and j.m . triscari , . e.m . wright , j. arlt , and k. dholakia , . aitken , _ determinants and matrices _ ( oliver and boyd , edinburgh and london , 1951 ) , p. 112 . r.b . laughlin , . r.a . smith and n.k . j. l. dubois and h. r. glyde , . m.d . girardeau and e.m . wright , cond - mat/0012365 ( 2000 ) . w.p . reinhardt and c.w . r. dum _ et al . scott , r.j . ballagh , and k. burnett , . jackson , g.m . kavoulakis , and c.j . pethick , . s. burger _ et al . _ , . j. denshlag _ et al . _ . a.e . muryshev _ busch and j.r . for a recent review see f. dalfovo _ et al . _ , . g. baym and c. j. pethick , . y. kagan , g. v. shlapnikov , and j. t. m. walraven , . m. edwards _ et al . s. stringari , . e. v. goldstein , k. plttner , and p. meystre , _ quantum . semiclass . opt . _ * 7 * , 743 ( 1995 ) . _ , _ science _ * 398 * , 218 ( 1999 ) . s. a. morgan , r. j. ballagh , and k. burnett , . w.p . reinhardt and c.w . m. naraschewski _ et al . _ , . j. javanainen and m. wilkens , ; _ ibid . _ * 81 * , 1345 ( 1998 ) . a. j. legget and f. sols , . e. m. wright _ et al . _ , . m.d . girardeau and e.m . k. e. cahill and r. j. glauber . j. durnin , j. j. miceli , and j. h. eberly , _ phys . * 58 * , 1499 ( 1987 ) . , to be published in _ phys . a _ ( 2001 ) . v. m. perez - garcia _ et al . _ , _ phys . a _ * 57 * , 3837 ( 1998 ) . y. s. kivshar and t. j. alexander , cond - mat/9905048 ( 1999 ) . s.  l. cornish _ et  al . _ , * 85 * , 1795 ( 2000 ) . c. henkel and m. wilkens , _ europhys . _ * 47 * , 414 ( 1999 ) . c. henkel , s. poetting , and m. wilkens , _ appl . b _ * 69 * , 379 ( 1999 ) . v. dunjko , v. lorent , and m. olshanii , cond - mat/0103085 ( 2001 ) .
Abstract:

Answer: several experimental groups are currently working towards realizing quasi - one - dimensional ( 1d ) atom de broglie waveguides and loading them with ultracold atoms . the dynamics becomes truly 1d in a regime ( tonks gas ) of low temperatures and densities , and large positive scattering lengths for which the transverse mode becomes frozen , in which case the many - body schrdinger dynamics becomes exactly soluble via a fermi - bose mapping theorem . in this paper we review our recent work on the exact ground state and quantum dynamics of 1d tonks gases and assess the possibilty of approaching the tonks regime using bessel beam optical dipole traps . # 1#2#3#4#1 * # 2 * , # 3 ( # 4 )