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: with recent advances in the field of ultra - cold atomic gases , trapped bose - einstein condensates ( bec s ) have become an important system to study many - body physics such as quantum phase transitions . in particular , the ability to dynamically control the lattice structure and the strength of interaction as well as the absence of impurities in bec s confined in optical lattices , have led to the recent observation of the superfluid to mott - insulator transition @xcite . at the same time , there has been a tremendous interest in studying rotating bec s in harmonic traps ; at sufficient rotation an abrikosov lattice of quantized vortices has been observed @xcite and the realization of strongly correlated quantum states similar to the fractional quantum hall states has been predicted to occur at higher rotation rates @xcite . in these proposals , the rotation can play the role of an effective magnetic field for the neutral atoms and in analogy with electrons , the atoms may enter into a state described by the laughlin wavefunction , which was introduced to describe the fractional quantum hall effect . while this approach yields a stable ground state separated from all excited states by an energy gap , in practice this gap is rather small because of the weak interactions among the particles in the magnetic traps typically used . in optical lattices , the interaction energies are much larger because the atoms are confined in a much smaller volume , and the realization of the fractional quantum hall effect in optical lattices could therefore lead to a much higher energy gap and be much more robust . in a recent paper @xcite , it was shown that it is indeed possible to realize the fractional quantum hall effect in an optical lattice , and that the energy gap achieved in this situation is a fraction of the tunneling energy , which can be considerably larger than the typical energy scales in a magnetic trap . in addition to being an interesting system in its own right , the fractional quantum hall effect , is also interesting from the point of view of topological quantum computation @xcite . in these schemes quantum states with fractional statistics can potentially perform fault tolerant quantum computation . so far , there has been no direct experimental observation of fractional statistics although some signatures has been observed in electron interferometer experiments @xcite . strongly correlated quantum gases can be a good alternative where the systems are more controllable and impurities are not present . therefore , realization of fractional quantum hall states in atomic gases can be a promising resource for topological quantum computation in the future . as noted above the fqh effect can be realized by simply rotating and cooling atoms confined in a harmonic trap . in this situation , it can be shown that the laughlin wavefunction exactly describes the ground state of the many body system @xcite . in an optical lattices , on the other hand , there are a number of questions which need to be addressed . first of all , it is unclear to which extent the lattices modifies the fractional quantum hall physics . for a single particle , the lattice modifies the energy levels from being simple landau levels into the fractal structure known as the hofstadter butterfly @xcite . in the regime where the magnetic flux going through each lattice @xmath0 is small , one expects that this will be of minor importance and in ref . it was argued that the fractional quantum hall physics persists until @xmath1 . in this paper , we extent and quantify predictions carried out in ref .   . whereas ref .   only considered the effect of infinite onsite interaction , we extent the analysis to finite interactions . furthermore , where ref .   mainly argued that the ground state of the atoms was in a quantum hall state by considering the overlap of the ground state found by numerical diagonalization with the laughlin wavefunction , we provide further evidence for this claim by characterizing the topological order of the system by calculating chern numbers . these calculations thus characterize the order in the system even for parameter regimes where the overlap with the laughlin wavefunction is decreased by the lattice structure . in addition to considering these fundamental features of the fqh states on a lattice , which are applicable regardless of the system being used to realize the effect , we also study a number of questions which are of particular interest to experimental efforts towards realizing the effect with atoms in optical lattices . in particular , we show that adding dipole interactions between the atoms can be used to increase the energy gap and stabilize the ground state . furthermore , we study bragg spectroscopy of atoms in the lattice and show that this is a viable method to identify the quantum hall states created in an experiment , and we discuss a new method to generate an effective magnetic field for the neutral atoms in the lattice . the paper is organized as follows : in sec . [ section - no_hardcore ] , we study the system with finite onsite interaction . in sec . [ section - chern ] we introduce chern numbers to characterize the topological order of the system . the effect of the dipole - dipole interaction has been elaborated in sec . [ section - dipole ] . [ nu1/4 ] studies the case of @xmath2 filling factor . [ bragg ] is dedicated to explore bragg spectroscopy of the system . [ magnetic ] outlines a new approach for generating the type of the hamiltonian studied in this paper . the fractional quantum hall effect occurs for electrons confined in a two dimensional plane under the presences of a perpendicular strong magnetic field . if @xmath3 is the number of electrons in the system and @xmath4 is the number of magnetic fluxes measured in units of the quantum magnetic flux @xmath5 , then depending on the filling factor @xmath6 the ground state of the system can form highly entangled states exhibiting a rich behaviors , such as incompressibility , charge density waves , and anyonic excitations with fractional statistics . in particular , when @xmath7 , where @xmath8 is an integer , the ground state of the system is an incompressible quantum liquid which is protected by an energy gap from all other states , and in the landau gauge is well described by the laughlin wavefunction @xcite :    ( z_1,z_2, ... ,z_n)=_j > k^n ( z_j - z_k)^m _ j=1^n e^- y^2_i/2 [ laughlin ] , where the integer @xmath8 should be odd in order to meet the antisymmetrization requirement for fermions . although the fractional quantum hall effect occurs for fermions ( electrons ) , bosonic systems with repulsive interactions can exhibit similar behaviors . in particular , the laughlin states with even @xmath8 correspond to bosons . in this article , we study bosons since the experimental implementation are more advanced for the ultra - cold bosonic systems . we study a system of atoms confined in a 2d lattice which can be described by the bose - hubbard model @xcite with the peierls substitution @xcite , @xmath9 + where @xmath10 is the hopping energy between two neighboring sites , @xmath11 is the onsite interaction energy , and @xmath12 is the phase acquired by a particle going around a plaquette . this hamiltonian is equivalent to the hamiltonian of a u(1 ) gauge field ( transverse magnetic field ) on a square lattice . more precisely , the non - interacting part can be written as -j_<ij > a^_i a_j ( ^j_i ) where @xmath13 is the vector potential for a uniform magnetic field and the path of the integral is chosen to be a straight line between two neighboring sites . in the symmetric gauge , the vector potential is written as @xmath14 . hence , @xmath0 will be the amount of magnetic flux going through one plaquette .    while the hamiltonian in eq . ( [ eq : hamiltonian ] ) occurs naturally for charged particles in a magnetic field , the realization of a similar hamiltonian for neutral particles is not straightforward . as we discuss in sec . [ magnetic ] this may be achieved in a rotating harmonic trap , and this have been very successfully used in a number of experiments in magnetic traps @xcite , but the situation is more complicated for an optical lattice . however , there has been a number of proposals for lattice realization of a magnetic field @xcite , and recently it has been realized experimentally @xcite . @xcite have studied the realization of fractional hall states for a few particles in individual lattice sites . a new approach for rotating the entire optical lattice is discussed in sec . [ magnetic ] . the essence of the above hamiltonian is a non - zero phase that a particle should acquire when it goes around a plaquette . this phase can be obtained for example by alternating the tunneling and adding a quadratic superlattice potential @xcite or by simply rotating the lattice ( sec.[magnetic ] ) . the advantage of confining ultra - cold gases in an optical lattice is to enhance the interaction between atoms which consequently result in a higher energy gap comparing to harmonic trap proposals ( e.g. ref . ) . this enhancement in the energy gap of the excitation spectrum can alleviate some of the challenges for experimental realization of the quantum hall state for ultra - cold atoms .      in order to approximate a large system , we study the system with periodic boundary condition , i.e. on a torus , where the topological properties of the system is best manifested . there are two energy scales for the system : the first is the magnetic tunneling term , @xmath15 , which is related to the cyclotron energy in the continuum limit @xmath16 and the second is the onsite interaction energy @xmath11 . experimentally , the ratio between these two energy scales can be varied by varying the lattice potential height @xcite or by feshbach resonances @xcite . let us first assume that we are in the continuum limit where @xmath17 , i.e. the flux through each plaquette is a small fraction of a flux quantum . a determining factor for describing the system is the filling factor @xmath18 , and in this study we mainly focus on the case of @xmath19 , since this will be the most experimentally accessible regime . we restrict ourself to the simplest boundary conditions for the single particle states @xmath20 , where @xmath21 is the magnetic translation operator which translates the single particle states @xmath22 around the torus . the definition and detailed discussion of the boundary conditions will be elaborated in section [ section - chern ] . the discussed quantities in this section , such as energy spectrum , gap and overlap , do not depend on the boundary condition angles ( this is also verified by our numerical calculation ) .    in the continuum case , for the filling fraction @xmath19 , the laughlin state in landau gauge @xmath23 is given by eq . ( [ laughlin ] ) with @xmath24 . the generalization of the laughlin wavefunction to a torus takes the form @xcite ( z_1,z_2, ... ,z_n)=f_rel(z_1,z_2, ... ,z_n)f_cm(z)e^-_i y^2_i/2 [ laughlin - torus ] , + where @xmath25 is the relative part of the wave function and is invariant under collective shifts of all @xmath26 s by the same amount , and @xmath27 is related to the motion of the center of mass and is only a function of @xmath28 . for a system on a torus of the size @xmath29 , we write the wavefunction with the help of theta functions , which are the proper oscillatory periodic functions and are defined as @xmath30 } ( z|\tau)=\sum_n e^{i \pi \tau ( n+a)^2 + 2\pi i ( n+a)(z+b)}$ ] where the sum is over all integers . for the relative part we have , f_rel=_i < j ( |i ) ^2 .    according to a general symmetry argument by haldane @xcite , the center of mass wave function @xmath27 is two - fold degenerate for the case of @xmath19 , and is given by f_cm(z)=(|2i ) where @xmath31 refers to the two degenerate ground states . this degeneracy in the continuum limit is due to the translational symmetry of the ground state on the torus , and the same argument can be applied to a lattice system when the magnetic length is much larger than the lattice spacing @xmath32 . for higher magnetic filed , the lattice structure becomes more pronounced . however , in our numerical calculation for a moderate magnetic field @xmath33 , we observe a two - fold degeneracy ground state well separated from the excited state by an energy gap . we return to the discussion of the ground state degeneracy in sec . [ section - chern ] .    in the continuum limit @xmath32 , the laughlin wavefunction is the exact ground state of the many body system with a short range interaction @xcite . the reason is that the ground state is composed entirely of states in the lowest landau level which minimizes the magnetic part of the hamiltonian , the first term in eq . ( [ eq : hamiltonian ] ) . the expectation value of the interacting part of the hamiltonian , i.e. the second term in eq . ( [ eq : hamiltonian ] ) , for the laughlin state is zero regardless of the strength of the interaction , since it vanishes when the particles are at the same position .    to study the system with a non - vanishing @xmath0 , we have performed a direct numerical diagonalization of the hamiltonian for a small number of particles . since we are dealing with identical particles , the states in the hilbert space can be labeled by specifying the number of particles at each of the lattice sites . in the hard - core limit , only one particle is permitted on each lattice site , therefore for @xmath3 particles on a lattice with the number of sites equal to @xmath34 , where @xmath35 is the unit lattice side , the hilbert space size is given by the combination @xmath36 . on the other hand , in case of finite onsite interaction , the particles can be on top of each other , so the hilbert space is bigger and is given by the combination @xmath37 . in our simulations the dimension of the hilbert can be raised up to @xmath38 and the hamiltonian is constructed in the configuration space by taking into account the tunneling and interacting terms . the tunneling term is written in the symmetric gauge , and we make sure that the phase acquired around a plaquette is equal to @xmath39 , and that the generalized magnetic boundary condition is satisfied when the particles tunnels over the edge of the lattice [ to be discussed in sec .  [ section - chern ] , c.f . ( [ twist_angles ] ) ] . by diagonalizing the hamilton , we find the two - fold degenerate ground state energy which is separated by an energy gap from the excited states and the corresponding wavefunction in the configuration space . the lauhglin wave function ( [ laughlin - torus ] ) can also be written in the configuration space by simply evaluating the laughlin wave function at discrete points , and therefore we can compared the overlap of these two dimensional subspaces . the energy gap above the ground state and the ground state overlap with the laughlin wavefunction for the case of @xmath19 in a dilute lattice @xmath40 , are depicted in figure . [ tuning_u ] . the laughlin wavefunction remains a good description of the ground state even if the strength of the repulsive interaction tends to zero ( fig . [ tuning_u ] a ) . below , we discuss different limits :    first , we consider @xmath41 : if the interaction energy scale @xmath11 is much larger than the magnetic one @xmath42 , all low energy states lie in a manifold , where the highest occupation number for each site is one , i.e. this corresponds to the hard - core limit . the ground state is the laughlin state and the excited states are various mixtures of landau states . the ground state is two - fold degenerate @xcite and the gap reaches the value in for the hard - core limit at large @xmath43 , as shown in fig . [ fig : gap ] ( a ) . in this limit the gap only depends on the tunneling @xmath10 and flux @xmath0 , and the gap is a fraction of @xmath10 . these results are consistent with the previous work in ref . . secondly , we consider @xmath44 . in this regime , the magnetic energy scale @xmath45 is much larger than the interaction energy scale @xmath11 . for repulsive regime ( @xmath46 ) , the ground state is the laughlin state and the gap increases linearly with @xmath47 , as shown in fig . [ fig : gap ] b.    thirdly , we study @xmath48 where the interaction is absent and the ground state becomes highly degenerate . for a single particle on a lattice , the spectrum is the famous hostadfer s butterfly @xcite , while in the continuum limit @xmath49 , the ground state is the lowest landau level ( lll ) . the single particle degeneracy of the lll is the number of fluxes going through the surface , @xmath50 . so in the case of @xmath3 bosons , the lowest energy is obtained by putting @xmath3 bosons in @xmath50 levels . therefore , the many - body ground state s degeneracy should be : @xmath51 for example , 3 bosons and 6 fluxes gives a 336-fold degeneracy in the non - interacting ground state .    if we increase the amount of phase ( flux ) per plaquette ( @xmath52 ) we are no longer in the continuum limit . the landau level degeneracy will be replace by @xmath53 where @xmath54 is the amount of flux per plaquette and @xmath55 and @xmath56 are coprime @xcite . then , the many - body degeneracy will be : @xmath57 . fourthly , we consider @xmath58 : when u is negative ( i.e. attractive interaction ) in the limit strong interaction regime , the ground state of the system will become a _ pillar _ state . in a pillar state , all bosons in the system condensate into a single site . therefore , the degeneracy of the ground state is @xmath59 and the ground state manifold can be spanned as , _ i ( a^)_i^n | vac . these states will very fragile and susceptible to collapse @xcite .    in a lattice , it is also possible to realize the fractional quantum hall states for attractive interaction in the limit when @xmath60 . assume that the occupation number of each site is either zero or one . since the attraction energy @xmath11 is very high and there is no channel into which the system can dissipate its energy , the probability for a boson to hop to a site where there is already a boson is infinitesimally small . therefore , the high energy attraction will induce an effective hard - core constraint in the case of ultra - cold system . the energy of these state should be exactly equal to their hard - core ground state counterparts , since the interaction expectation value of the interaction energy is zero for the laughlin state . the numerical simulation shows that these two degenerate states indeed have a good overlap with the laughlin wavefunction similar to their repulsive hard - core counterparts and also their energies are equal to the hard - core ground state . these states are very similar to repulsively bound atom pairs in an optical lattice which have recently been experimentally observed @xcite . is varied by changing the size of the lattice ( the size in the two orthogonal directions differ at most by unity ) . the laughlin state ceases to be a good description of the system as the lattice nature of the system becomes more apparent @xmath61 . the overlap is only calculated at the positions shown with dots and the color coding is an extrapolation between the points . ] is varied by changing the size of the lattice ( the size in the two orthogonal directions differ at most by unity ) . the laughlin state ceases to be a good description of the system as the lattice nature of the system becomes more apparent @xmath61 . the overlap is only calculated at the positions shown with dots and the color coding is an extrapolation between the points . ] so far we have mainly considered a dilute lattice @xmath62 , where the difference between a lattice and the continuum is very limited . we shall now begin to investigate what happens for larger values of @xmath0 , where the effect of the lattice plays a significant role . [ high - alpha - overlap ] , shows the ground state overlap with the laughlin wavefunction as a function of the strength of magnetic flux @xmath0 and the strength of the onsite interaction @xmath11 . as @xmath0 increases the laughlin overlap is no longer an exact description of the system since the lattice behavior of the system is more pronounced comparing to the continuum case . this behavior does nt depend significantly on the number of particles for the limited number of particles that we have investigated @xmath63 . we have , however , not made any modification to the laughlin wave function to take into account the underlying lattice , and from the calculations presented so far , it is unclear whether the decreasing wave function overlap represents a change in the nature of the ground state , or whether it is just caused by a modification to the laughlin wave function due the difference between the continuum and the lattice . to investigate this , in the next section , we provide a better characterization of the ground state in terms of chern numbers , which shows that the same topological order is still present in the system for higher values of @xmath0 .    as a summary , we observe that the laughlin wavefunction is a good description for the case of dilute lattice ( @xmath32 ) , regardless of the strength of the onsite interaction . however , the protective gap of the ground state becomes smaller for weaker values of interaction and in the perturbative regime @xmath64 is proportional to @xmath47 for @xmath65 . in the theory of quantum hall effect , it is well understood that the conductance quantization , is due to the existence of certain topological invariants , so called chern numbers . the topological invariance was first introduced by avron _ _ et al.__@xcite in the context of thouless _ et al . _ ( tkndn ) s original theory @xcite about quantization of the conductance . tkndn in their seminal work , showed that the hall conductance calculated from the kubo formula can be expressed into an integral over the magnetic brillouin zone , which shows the quantization explicitly . the original paper of tkndn deals with the single - particle problem and bloch waves which can not be generalized to topological invariance . the generalization to many - body systems has been done by niu _ _ et al.__@xcite and also tao _ _ et al.__@xcite , by manipulating the phases describing the closed boundary conditions on a torus ( i.e. twist angles ) , both for the integer and the fractional hall systems . these twist angles come from natural generalization of the closed boundary condition . to clarify the origin of these phases , we start with a single particle picture . a single particle with charge ( @xmath66 ) on a torus of the size @xmath67 in the presence of a magnetic field @xmath68 perpendicular to the torus surface , is described by the hamiltonian h_s= , where @xmath13 is the corresponding vector potential ( @xmath69 ) . this hamiltonian is invariant under the magnetic translation ,    t_s(*a*)=e^i * * ak**^s / where @xmath70 is a vector in the plane , and @xmath71 is the pseudomomentum , defined by    @xmath72    the generalized boundary condition on a torus is given by the single - particle translation    @xmath73    where @xmath74 and @xmath75 are twist angles of the boundary . the origin of these phases can be understood by noting that the periodic boundary conditions corresponds to the torus in fig . [ torus](a ) . the magnetic flux through the surface arises from the field perpendicular to the surface of the torus . however , in addition to this flux , there may also be fluxes due to a magnetic field residing inside the torus or passing through the torus hole , and it is these extra fluxes which give rise to the phases . the extra free angles are all the same for all particles and all states in the hilbert space , and their time derivative can be related to the voltage drops across the hall device in two dimensions . the eigenstates of the hamiltonian , including the ground state will be a function of these boundary angles @xmath76 . by defining some integral form of this eigenstate , one can introduce quantities , that do not depend on the details of the system , but reveal general topological features of the eigenstates . first we discuss the simplest situation , where the ground state is non - degenerate , and later we shall generalize this to our situation with a degenerate ground state . the chern number is in the context of quantum hall systems related to a measurable physical quantity , the hall conductance . the boundary averaged hall conductance for the ( non - degenerate ) @xmath0th many - body eigenstate of the hamiltonian is @xcite : @xmath77 , where the chern number @xmath78 is given by    [ twist ] ) ( b ) redefining the vector potential around the singularities : @xmath79 is not well - defined everywhere on the torus of the boundary condition . therefore , another vector field @xmath80 with different definition should be introduced around each singularity @xmath81 of @xmath79 . @xmath79 and @xmath80 are related to each other by a gauge transformation @xmath82 and the chern number depends only on the loop integrals of @xmath82 around those singularities regions , c.f . , eq.([chern_gauge]).,title="fig : " ]    [ gauge ] ) ( b ) redefining the vector potential around the singularities : @xmath79 is not well - defined everywhere on the torus of the boundary condition . therefore , another vector field @xmath80 with different definition should be introduced around each singularity @xmath81 of @xmath79 . @xmath79 and @xmath80 are related to each other by a gauge transformation @xmath82 and the chern number depends only on the loop integrals of @xmath82 around those singularities regions , c.f . , eq.([chern_gauge]).,title="fig : " ]    c()= _ 0 ^ 2 d_1 _ 0 ^ 2 d_2 ( _ 1 _ 2^()-_2 _ 1^ ( ) ) [ eq : chern ] , where @xmath83 is defined as a vector field based on the eigenstate @xmath84 on the boundary torus @xmath85 by _ j^()(_1,_2 ) i ^()||^ ( ) . it should be noted that the wave function @xmath76 is defined up to a phase factor on the boundary - phase space . therefore , @xmath76 and @xmath86 are physically equivalent for any smooth function @xmath87 . under this phase change , @xmath83 transforms like a gauge :    @xmath88    hence , the chern number integral is conserved under this gauge transformation and it encapsulates general properties of the system . chern numbers has been used extensively in the quantum hall literature , for characterization of the localized and extended states ( ref .   and refs . therein ) . in this paper , we use the chern number as an indicator of order in the system . moreover , it enables us to characterize the ground state in different regimes , especially where the calculation of the overlap with the laughlin wave function fails to give a conclusive answer .    before explaining the method for calculating the chern number , we clarify some issues related to the degeneracy of the ground state . in some systems , the ground state can be degenerate , this can be intrinsic or it can be as a result of the topology of the system . if the ground state is degenerate , we should generalize the simple integral form of eq . ( [ eq : chern ] ) to take into account the redundancy inside the ground state manifold . for example , in the case of a two - fold degenerate ground state , an extra gauge freedom related to the relative phase between two ground states , and this freedom should be fixed . in other words , as we change the twist angles , we can not keep track of the evolution of both states , since one can not distinguish them from each other . therefore , to uniquely determine the chern number of the ground state(s ) , we should resolve this gauge invariance , which is treated in section [ section - impurity ] and [ gauge - fixing ] . it important to note that the degeneracy in the non - interacting regime is fundamentally different from the degeneracy in the interacting case . in the non - interacting limit , the degeneracy can be lifted by a local perturbation e.g. impurity , while in the hardcore case , the degeneracy remains in the thermodynamic limit @xcite . the latter degeneracy in the ground state is a consequence of the global non - trivial properties of the manifold on which the particles move rather than related to a symmetry breaking which happens in conventional models e.g. ising model . the topological degeneracy is not a consequence of breaking of any symmetry only in the ground state , instead it is the entire hilbert space which is split into disconnected pieces not related by any symmetry transformation . with a general argument , wen @xcite showed that if the degeneracy of a chiral spin system moving on a simple torus is @xmath89 , then it should be @xmath90 on a torus with @xmath91 handles ( riemann surface with genus @xmath91 ) , therefore the topological detergency is an intrinsic feature of the system . in particular , in the context of quantum hall effect , this multicomponent feature of the ground state on a torus has a physical significance , while the single component ground state on a sphere boundary condition gives zero conductance , the torus geometry with multicomponent ground state results in a correct conductance measured in the experiment , since the torus boundary condition is more relevant to the experiment . changing twists angles of the boundary will rotate these components into each other and gives an overall non - zero conductance @xcite . as studied in a recent work by m. oshikawa and t. senthil @xcite , as a universal feature , it has been shown that in presence of a gap , there is a direct connection between the fractionalization and the topological order . more precisely , once a system has some quasiparticles with fractional statistics , a topological degeneracy should occur , which indicates the presence of a topological order . therefore , the amount of the degeneracy is related to the statistics of the fractionalized quasiparticles e.g. in the case of @xmath19 , the two - fold degeneracy is related to 1/2 anyonic statistics of the corresponding quasiparticles . chern number has been also studies for spin @xmath92 system on a honeycomb lattice @xcite for identifying abelian and non - abelian anyons . _ @xcite studies chern number for disordered fermi system using a non - commutative geometry .    to resolve the extra gauge freedom related to the two degenerate ground states , we consider two possibilities : i. lifting the degeneracy by adding some impurities , ii . fixing the relative phase between the two states in the ground state . below , we explore both cases . in the first case , we introduce some fixed impurities to lift the degeneracy in the ground state for all values of the twist angles . this is an artifact of the finite size of the system which we take advantage of . in the presence of perturbation , we show that the system has a topological order in spite of poor overlap with the laughlin state . in the second approach , we use a scheme recently proposed by hatsugai @xcite which is a generalized form for degenerate manifolds .      in this section , we introduce some perturbation in the finite system in form of local potentials ( similar to local impurities in electronic systems ) to split the degeneracy of the ground state and resolve the corresponding gauge invariance , which allows us to compute the chern number . furthermore , the fact that we can still uniquely determine the chern number in the presence of impurities , shows that the system retains its topological order , even when the impurities distort its ground state wavefunction away from the laughlin wavefunction .    in the context of the quantum hall effect , the conventional numerical calculation of various physical quantities such as the energy spectrum , screening charge density profile , wave functions overlaps , and the density - density correlation functions , can not be used for understanding the transport properties of electrons in the presence of impurities ( although useful for studying of isolated impurities @xcite ) . recently , d. n. sheng _ et al . _ @xcite calculated the chern number as an unambiguous way to distinguish insulating and current carrying states in the @xmath93 quantum hall regime which correspond to zero and non - zero chern number , respectively . in this work , a weak random disorder was used to lift the degeneracy of the ground state ( three - fold for @xmath93 ) for a finite number of electrons . the energy splitting between the lowest three states then decreased with increasing number of particles , which indicates the recovery of the degeneracy in the thermodynamic limit . moreover , the mobility gap can be determined by looking at the energy at which the chern number drops towards zero . this energy gap is comparable to the energy gap obtained from the experiment and it is not necessarily equal to the spectrum gap which separates the degenerate ground state from the excited states this shows the significance of chern number calculations for understanding these systems .    in a finite system , the coupling to a single - body interaction , e.g. impurities , can lift the degeneracy and one can uniquely determine the chern number for the individual states by direct integration of eq . ( [ eq : chern ] ) . on one hand , the impurity should be strong enough to split the energy in the ground state ( in this case @xmath94 , where @xmath95 denotes the energy of the @xmath96th energy level ) for all values of the twist angles . on the other hand , these impurities should be weak enough so that the energy splitting in the ground state remains smaller than the thermodynamic gap ( in this case @xmath97 ) .    to calculate the chern number of individual level , as mentioned in the pervious section , we have to fix the phase of the wavefunction . the method that we explore in this section can be considered as a simplified version of the general method developed by hatsugai @xcite which we will explore in the next section . following kohmoto s procedure @xcite , we assume that the ground state @xmath98 may be expanded for all twist angles on a @xmath99-dimensional hilbert discrete space @xmath100 . if @xmath101 in eq . ( [ eq : chern ] ) is a periodic function on the torus of the boundary condition , then by application of stoke s theorem , the chern number will be always zeros . the non - triviality ( non - zero conductance in the case of quantum hall system ) occurs because of the zeros of the wave function , where the phase is not well - defined . therefore , @xmath102 is not well defined everywhere and its flux integral can be non - zero . to uniquely determine the chern number , we assume that the wave function and the vector field are not defined for certain points @xmath81 in @xmath103 regions on the torus of the boundary condition . for simplicity , we first discuss this procedure , in the case of a non - degenerate ground state . for calculating the integral , we should acquire another gauge convention for the wave function inside these @xmath103 regions , e.g. in a discrete system , we may require an arbitrary element of the wave function to be always real , and thereby we can define a new vector field @xmath104 , which is well defined inside these regions . these two vector fields differ from each other by a gauge transformation ( fig . [ torus ] ) : _ j^()(_1,_2)- _ j^()(_1,_2 ) = _ j ( _ 1,_2 ) , [ eq:1d_gauge]and the chern number reduces to the winding number of the gauge transformation @xmath105 over small loops encircling @xmath81 , i.e. @xmath106 ,    c()= _ n _ s_n d [ chern_gauge ] .    the one - dimensional gauge eq .  ( [ eq:1d_gauge ] ) should be resolved by making two conventions . for example , in one convention the first element and in the other the second element of the wavefunction in the hilbert space should be real i.e. transforming the ground state @xmath107 into @xmath108 where @xmath109 is a @xmath99-dimensional vector and @xmath110 is a projection into the ground state and similarly with the other reference vector @xmath111 . since the gauge that relates two vector fields is the same as the one that relates the corresponding wavefunctions ( similar to eq.([eq : gauge_wavefunction ] ) ) , we can uniquely determine the gauge transformation function @xmath82 by evaluating @xmath112 . therefore , the chern number will be equal to the number of windings of @xmath82 around regions where @xmath113 is zero . counting the vorticities has a vigorous computational advantage over the conventional method of direct integration of eq . ( [ eq : chern ] ) . in the direct integration , we need to choose a large number of mesh points for the boundary angles , because of the discrete approximation of derivatives in eq . ( [ eq : chern ] ) , and this makes the calculation computationally heavy . we note that for the system on a lattice , we should exactly diagonalize the hamiltonian which is a sparse matrix as opposed to the continuum case where the hamiltonian is a dense matrix residing on a smaller projected hilbert space ( lowest landau level ) .     and @xmath114 , for the first and the second energy state , respectively . the results are for 3 atoms on a 6x6 lattice ( @xmath115 ) in the hard - core limit . ( b ) @xmath116 for the first level has no vorticity . however , for the second level , as shown in ( c ) , @xmath116 has vorticity equal to one associated with regions where either @xmath117 or @xmath118 vanishes . , width=566 ]    for removing the degeneracy , in our numerical simulations , we add a small amount of impurity which is modeled as delta function potentials located at random sites with a random strength , of the order of the tunneling energy @xmath10 . this is described by a hamiltonian of the form @xmath119 , where @xmath120 numerates the lattice site , @xmath121 is the atom number operator and @xmath122 is the strength of the impurity at site @xmath120 . we choose reference states @xmath123 and @xmath124 to be eigenvectors of the numerically diagonalized hamiltonian at two different twist angles . in fig . [ impurity ] , vorticities of @xmath125 associated with the first and the second energy level is depicted . it is easy to see that the chern number associated to the two ground states is one . number of vortices may vary for the first and second ground states , but their sum is always equal to one . the hard - core limit ( @xmath126 ) is very similar to the case of fractional quantum hall effect , which in the context of hall systems , means a share of @xmath92 ( in @xmath127 unit ) for each ground state @xcite . when the onsite interaction strength is small ( @xmath128 ) , the thermodynamic gap becomes comparable to the ground state energy splitting @xmath129 , the chern number can not be uniquely determined , and the system does nt have topological order . on the other hand , in the limit of strong interaction ( @xmath130 ) , the total chern number associated to the ground states is equal to one , regardless of the impurity configuration . moreover , in the hard - core limit , although the ground state is not described by the laughlin wavefunction , since it is distorted by the impurity ( in our model it can be as low as @xmath131 ) , the chern number is unique and robust . this is an indication that the topological order is not related to the symmetries of the hamiltonian and it is robust against arbitrary local perturbations @xcite . these results indicate the existence of a topological order in the system and robustness of the ground states against local perturbations . the method developed in the previous section has the graphical vortex representation for the chern number which makes it computationally advantageous compared to the direct integration of eq . ( [ eq : chern ] ) . it can not , however , be applied directly to a degenerate ground state , and therefore we had to introduce an impurity potential which lifted the degeneracy . on the other hand , a significant amount of impurity in the system may distort the energy spectrum , so that the underlying physical properties of the lattice and fluxes could be confounded by the artifacts due to the impurities , especially for large @xmath0 . to address this issues , in this section , we explore a generalized method of the previous section based on refs .   and , which works for a degenerate ground state . as arrows for fixed @xmath123 and @xmath124 . ( b ) and ( c ) : surface plots of @xmath132 and @xmath133 ( blue is lower than red ) . @xmath74 and @xmath75 changes from zero to @xmath134 . these plots have been produced for 3 atoms with @xmath135 ( @xmath136 ) in the hard - core limit on a 5x5 lattice . the total vorticity corresponding to each of the reference wave functions ( @xmath123 or @xmath124 ) indicates a chern number equal to one . ]     as arrows for fixed @xmath123 and @xmath124 . ( b ) and ( c ) : surface plots of @xmath132 and @xmath133 ( blue is lower than red ) . @xmath74 and @xmath75 changes from zero to @xmath134 . these plots have been produced for 3 atoms with @xmath135 ( @xmath136 ) in the hard - core limit on a 5x5 lattice . the total vorticity corresponding to each of the reference wave functions ( @xmath123 or @xmath124 ) indicates a chern number equal to one . ]     as arrows for fixed @xmath123 and @xmath124 . ( b ) and ( c ) : surface plots of @xmath132 and @xmath133 ( blue is lower than red ) . @xmath74 and @xmath75 changes from zero to @xmath134 . these plots have been produced for 3 atoms with @xmath135 ( @xmath136 ) in the hard - core limit on a 5x5 lattice . the total vorticity corresponding to each of the reference wave functions ( @xmath123 or @xmath124 ) indicates a chern number equal to one . ]    by generalizing the chern number formalism for a degenerate ground state manifold , instead of having a single vector field @xmath137 , a tensor field @xmath138 should be defined , where @xmath139 for a @xmath140-fold degenerate ground state    _ j^(,)(_1,_2 ) i ^()||^ ( )    similar to the non - degenerate case , when @xmath141 is not defined , a new gauge convention should be acquired for the regions with singularities . this gives rise to a tensor gauge transformation @xmath142 on the border of these regions    _ j^(,)(_1,_2)- _ j^(,)(_1,_2 ) = _ j ^(,)(_1,_2 ) . following hatsugai s proposal@xcite for fixing the ground state manifold gauge , we take two reference multiplets @xmath123 and @xmath124 which are two arbitrary @xmath143 matrices ; @xmath66 is the ground state degeneracy ( equal to 2 in our case ) . in our numerical simulation , we choose the multiplets to be two sets of ground state at two different twist angles far from each other , e.g. @xmath144 and @xmath145 . we define an overlap matrix as @xmath146 where @xmath147 is again the projection into the ground state multiplet , and consider the regions where @xmath148 or @xmath149 vanishes ( similar to zeros of the wave function in the non - degenerate case ) . hence , the chern number for @xmath140 degenerate states , will be equal to the total winding number of @xmath150 for small neighborhoods , @xmath103 , in which @xmath148 vanishes     for fixed @xmath123 and @xmath124 . @xmath74 and @xmath75 changes form zero to @xmath134 . this plot has been produced for 4 atoms with @xmath151 in the hard - core limit on a 5x5 lattice ( @xmath152 ) . although , there are more vortices here compared to fig . [ omega ] , the total vorticity corresponding to each of the trial functions ( @xmath123 or @xmath124 ) indicates a chern number equal to one . ] @xmath153    which is the same as the number of vortices of @xmath154 . it should be noted that the zeros of det @xmath155 and det@xmath156 should not coincide in order to uniquely determine the total vorticity . in fig . [ omega ] , we have plotted @xmath157 , @xmath132 , and @xmath158 , found by numerical diagonalization of the hamiltonian for a mesh ( @xmath159 ) of winding angles @xmath74 and @xmath75 . in this figure , the chern number can be determined be counting the number of vortices and it is readily seen that the winding number is equal to one for the corresponding zeros of @xmath148 ( or @xmath160 ) . .chern number for different configurations in the hard - core limit for fixed filling factor @xmath19 . the laughlin state overlap is shown in the last column . although it deviates from the laughlin state . although the ground state deviates from the laughlin state , the chern number remains equal to one half per state before reaching some critical @xmath161 where the energy gap vanishes . [ cols="^,^,^,^,^,^,^,^",options="header " , ] we have calculated the chern number for fixed @xmath19 and different @xmath0 s by the approach described above . the result is shown in table [ chern ] . for low @xmath32 , we know from sec . [ section - no_hardcore ] that the ground state is the laughlin state and we expect to get a chern number equal to one . for higher @xmath0 , the lattice structure becomes more apparent and the overlap with the laughlin state decreases . however , in our calculation , the ground state remains two - fold degenerate and it turns out that the ground state chern number tends to remain equal to one before reaching some critical @xmath161 . hence , also in this regime we expect to have similar topological order and fractional statistics of the excitations above these states on the lattice . the degeneracy of the ground state is a function of twist angles . the shown results are for 5 atoms on a 5x5 lattice i.e. @xmath162 ( a ) shows first three energy manifolds as a function of the toroidal boundary condition angles , ( b ) shows a cross section of ( a ) at @xmath163 for seven lowest energy levels . the first and the second energy levels get close to each other at @xmath164 . ]     the degeneracy of the ground state is a function of twist angles . the shown results are for 5 atoms on a 5x5 lattice i.e. @xmath162 ( a ) shows first three energy manifolds as a function of the toroidal boundary condition angles , ( b ) shows a cross section of ( a ) at @xmath163 for seven lowest energy levels . the first and the second energy levels get close to each other at @xmath164 . ]    for the arguments above to be applicable , it is essential that we can uniquely identify a two - fold degenerate ground state which is well separated from all higher lying states . for higher flux densities , @xmath165 , the two - fold ground state degeneracy is no longer valid everywhere on the torus of the boundary condition . in this regime , the issue of degeneracy is more subtle , and the finite size effect becomes significant . the translational symmetry argument @xcite , which was used in the section [ section - no_hardcore ] , is not applicable on a lattice and as pointed out by kol _ et al_. @xcite the degeneracy of the ground state may vary periodically with the system size . some of the gaps which appear in the calculation may be due to the finite size and vanish in the thermodynamic limit , whereas others may represent real energy gaps which are still present in the thermodynamic limit . to investigate this , we study the ground state degeneracy as a function of boundary angles ( @xmath166 ) which are not physical observable and therefore the degeneracy in thermodynamic limit should not depend on their value . in particular , fig .  [ level - crossing ] shows the energy levels of five particles at @xmath167 for different values of the twist angles . the first and the second level are split at ( @xmath168 ) , while they touch each other at ( @xmath164 ) . we have observed similar behavior for different number of particles and lattice sizes e.g. 3 and 4 atoms at @xmath169 . in this case , the system seems to not have a two - fold degeneracy . therefore , the ground state enters a different regime which is a subject for further investigation .    for having the topological order , it is not necessary to be in the hard - core limit . even at finite interaction strength @xmath170 , we have observed the same topological order with the help of the chern number calculation . if @xmath11 gets further smaller , the energy gap above the ground state diminishes ( as seen in sec . [ section - no_hardcore ] ) and the topological order disappears . we conclude that the chern number can be unambiguously calculated for the ground state of the system in a regime where laughlin s description is not appropriate for the lattice . the non - zero chern number of a two - fold degenerate ground state , in this case equal to one half per state , is a direct indication of the topological order of the system . in sections [ section - no_hardcore ] and [ section - chern ] above , we have investigated the conditions under which the fractional quantum hall effect may be realized for particles on a lattice . the motivation for this study is the possibility to generate the quantum hall effect with ultra cold atoms in an optical lattice but the results of these sections are applicable regardless of the method used to attain this situation . in this and the following sections , we investigate some questions which are of particular relevance to ultra cold atoms in an optical lattice . first , we introduce a long range , e.g. , dipole - dipole , interaction which turn out to increase the energy gap and thereby stabilizes the quantum hall states . we then turn to the case of @xmath2 and show that in order to realize this state , it is essential to have some kind of long range interaction .      in an experimental realization of the quantum hall effect on a lattice , it is desirable to have as large an energy gap as possible in order to be insensitive to external perturbations . so far , we have studied effect of the short range interaction , and we have shown that the gap increases with increasing interaction strength , but the value of the gap saturates when the short range interaction becomes comparable to the tunneling energy @xmath10 .    in this section , we explore the possibility of increasing the gap by adding a long - range repulsive dipole - dipole interaction to the system . previously , such dipole - dipole interaction has also been studied in ref . + onlinecitecooper2005 as a method to achieve read - rezayi states @xcite of rapidly rotating cold trapped atoms for @xmath171 and as a means to realize fractional quantum hall physics with fermi gases @xcite . the dipole - dipole ( magnetic or electric ) interaction is described by the hamiltonian : h_d - d = u_dipole _ 1i < jn where @xmath172 . @xmath173 are unit vectors representing the permanent dipole moments and the position vectors @xmath174 are in units of the lattice spacing @xmath175 . for simplicity , we assume that all dipoles are polarized in the direction perpendicular to the plane . with time independent dipoles , the strength of the interaction is given by @xmath176 ( or @xmath177 ) where @xmath178 s ( @xmath179 s ) are the permanent magnetic ( electric ) dipole moment . static dipoles will thus give repulsive interaction @xmath180 , but experimentally time varying fields may be introduced which effectively change the sign of the interaction@xcite . for completeness , we shall therefore both investigate positive and negative @xmath181 , but the repulsive interaction corresponding to static dipoles will be the most desirable situation since it stabilizes the quantum hall states .    experimentally the dipole - dipole interaction will naturally be present in the recently realized bose - einstein condensation of chromium@xcite which has a large magnetic moment . however , for a lattice realization , polar molecules which have strong permanent electric dipole moments is a more promising candidate . for typical polar molecules with the electric moment @xmath182 debye , on a lattice with spacing @xmath183@xmath184 m , @xmath185 can be up to few khz , an order of magnitude greater than the typical tunneling @xmath186 which can be few hundreds of hz @xcite .     versus @xmath0 . the results are shown for @xmath187 , @xmath188 , @xmath189 ]     versus @xmath0 . the results are shown for @xmath187 , @xmath188 , @xmath189 ]    to study the effect of the dipole - dipole interaction , we again numerically diagonalize the hamiltonian for a few hardcore bosons ( @xmath130 ) , in the dilute regime @xmath190 , while varying strength of the dipole - dipole interaction . the results of this simulation is shown in fig . [ dipole ] ( a ) . after adding the dipole interaction , the ground state in a dilute lattice remains two - fold degenerate , since the interaction only depends on the relative distance of the particles and keeps the center of mass motion intact . if the dipole interaction becomes too strong @xmath191 , the ground state wave function deviates from the laughlin wavefunction , but the topological order remains the same as for the system without dipole interaction . we verified this by calculating the chern number as explained in sec . [ section - chern ] for different values of the dipole - dipole interaction strength @xmath192 and observed that the total chern number of the two - fold degenerate ground state is equal to one . moreover , as it is shown in fig . [ dipole ] ( b ) adding such an interaction can increase the gap : the lower curve corresponds to the hard - core limit discussed in the previous work @xcite and the upper curve corresponds to the system including the dipole - dipole interaction . this enhancement varies linearly with the flux density @xmath0 in a dilute lattice and does nt depend on the number of particles and consequently , it is expected to behave similarly in the thermodynamic limit . one of the impediment of the experimental realization of quantum hall state is the smallness of the gap which can be improved by adding dipole - dipole interaction . in this section , we showed that this improvement is possible and moreover , by chern number evaluation , we verified that adding dipole interaction does nt change the topological behavior of the ground state manifold . so far we have concentrated on the case of @xmath19 . in this section , we briefly investigate the case of @xmath2 . it is expected that the laughlin wavefunction remains a good description for the ground state of a bosonic system for any even @xmath66 , where @xmath193 . following haldane s argument@xcite , due to the center of mass motion , the degeneracy of the ground state is expected to be @xmath66-fold on a torus . similar to the case of @xmath19 , the laughlin wavefunction should be a suitable description for any @xmath66 provided that the magnetic field is weak so that we are close to the continuum limit , i.e. @xmath194 . also the chern number is expected to be equal to one for the q - fold degenerate ground state , which in the context of quantum hall effect means a share of @xmath195 of the conduction quantum @xmath127 for each state in the q - fold degenerate ground state manifold .    we have done both overlap and the chern number calculation to check those premises . in the case of @xmath2 , significant overlap occurs at low @xmath196 . the average wave function overlap of four lowest energy eigenstates with the laughlin wavefunction is depicted in figure [ overlap_laughlin_1_4 ] , where we have used a generalization of the laughlin wavefunction for periodic boundary conditions similar to eq . ( [ laughlin - torus ] ) @xcite . we observe that the laughlin wavefunction is a reliable description of the system with @xmath2 but only for much more dilute lattices ( @xmath196 ) compared to @xmath19 where significant overlap occurs for @xmath197 . contrary to @xmath19 , where the gap is a fraction of the tunneling energy @xmath10 , the gap for @xmath2 between the 4-fold degenerate ground state and the next energy level is infinitesimal . the reason for the vanishing gap can be understood in the continuum limit from the argument put forward in ref .  : as noted previously the laughlin wavefunction is an exact eigenstate of the hamiltonian with an energy per particle equal to the lowest landau level energy . the energy of the @xmath198 state is thus equal to the @xmath24 state . it thus costs a negligible energy to compress the @xmath2 state to the @xmath19 state , and therefore there is no gap . in an external trap the system will always choose the highest density state which is the @xmath19 state . note , however that this argument only applies to short range interactions . for long range interactions , we expect to see a non vanishing gap .    even though that with short range interactions , the gap is very small in our numerical calculations , it is still sufficiently large that it allows us to unambiguously determine the chern number for the ground state manifold as described in sec . [ gauge - fixing ] . as expected the calculation shows that the chern number is equal to one corresponding to a four - fold degenerate ground state consistent with the generalization of the fermionic case in the fractional quantum hall theory @xcite . in fig . [ overlap_laughlin_1_4 ] , the overlap of the first four lowest energy state with the laughlin wavefunction is depicted . in the absence of the dipole interaction , the ground state overlap is significant only for @xmath199 , however , by adding a moderate dipole interaction ( @xmath200 ) , the overlap becomes more significant for a larger interval of the flux density , i.e. @xmath201 . this is due to the fact that state with lower density become more favorable in the presence of a long - range repulsive interaction . we observed that adding a dipole interaction would lead to an improvement of the gap for @xmath19 and make the laughlin overlap more significant for larger interval of magnetic field strength @xmath0 in the case of @xmath2 . therefore this long - range interaction can be used as an tool for stabilizing the ground state and make the realization of these quantum states , experimentally more feasible .     on a torus . the dashed ( dotted ) line shows the overlap for the system without ( with ) dipole interaction ( @xmath200 ) . the laughlin state is only a good description for a more dilute lattice @xmath202 compared to @xmath19 . the dipole interaction stabilizes the system and the overlap is more significant for higher values of @xmath203 . ] in an experimental realization of the quantum hall states , it is essential to have an experimental probe which can verify that the desired states were produced . in most experiments with cold trapped atoms , the state of the system is probed by releasing the atoms from the trap and imaging the momentum distribution . in ref . , it was shown that this technique provide some indication of the dynamics in the system . this measurement technique , however , only provides limited information , since it only measures the single particle density matrix , and provides no information about the correlations between the particles . in refs .   and a more advanced measurement techniques were proposed , where the particle correlation is obtained by looking at the correlations in the expansion images . in this section , we study bragg scattering as an alternative measurement strategy which reveals the excitation spectrum of the quantum system . in ref . bosonic quantum hall system responses to a perturbative potential is studied . we focus on bragg scattering where two momentum states of the same internal ground state are connected by a stimulated two - photon process @xcite . by setting up two laser beams with frequencies @xmath204 and @xmath205 and wave vectors @xmath206 and @xmath207 in the plane of the original lattice , a running optical superlattice will be produced with frequency @xmath208 and wave vector @xmath209 . ( both frequencies @xmath204 and @xmath205 should be close to an internal electronic dipole transition in the atoms ) . the beams should be weak and sufficiently detuned so that direct photon transitions are negligible , i.e. @xmath210 , where @xmath211 is the frequency of the transition , @xmath212 is the corresponding spontaneous decay rate and @xmath213 are the rabi frequencies related to the laser - atom coupling . in this perturbative regime , the inelastic scattering of photons will be suppressed ; the atom starts from the internal ground state , absorbs one photon from e.g. beam 1 by going to a virtual excited state and then emits another photon into beam 2 by returning to its internal ground state . after this process , the system has acquired an energy equal to @xmath214 and a momentum kick equal to @xmath215 . therefore , the overall effect of the recoil process is a moving ac stark shift as a perturbing moving potential , and the effective hamiltonian represents the exchange of the two - photon recoil momentum and the energy difference to the system and is proportional to the local density i.e. @xmath216 . this process can be used to probe density fluctuations of the system and thus to measure directly the dynamic structure factor @xmath217 and static structure factor @xmath218 . this kind of spectroscopy has been studied for a bec in a trap by blakie _ _ @xcite and zambelli _ et al . _ @xcite and has been realized experimentally in refs .  . ]    . ] also , bragg spectroscopy in an optical lattice is discussed in ref .   in a mean - field approach and also in ref .   as a probe of the mott - insulator excitation spectrum . on the other hand , in the context of quantum hall effect , the static structure factor has been studied for probing the magnetoroton excitations @xcite and charge density waves @xcite . the dynamic and static structure factors are given respectively as    s(,)=_n , 0| n | ^()| 0|^2 ( -e_n+e_0 ) ,    s()=_n,0 | n | ^()| 0|^2 = _ n,0 | n | _ e^i | 0|^2 , where the density fluctuation operator is defined as @xmath219 and the coefficient are defined as fourier transforms of the wannier functions : @xmath220 , where the wannier function @xmath221 is the wave function of an atom localized on a site centered at @xmath222 . below , we focus on deep optical lattices , where @xmath223 .    in the structure factor , there is a sum over the excited states @xmath224 and ground states @xmath225 and the self - term is thus excluded . the ground state on a torus is two - fold degenerate and therefore in our numerics , we add the contribution of both . since we are working on a discrete lattice , there will be a cut - off in the allowed momentum given by the lattice spacing @xmath226 , where @xmath227 is the distance between lattice sites . [ 3atoms - structure](a ) shows the structure factor for the case @xmath19 for a small @xmath0 calculated from our numerical diagonalizations . in the data presented here , we have chosen @xmath228 and but the result should be similar in other directions in the lattice plane . we see that @xmath218 is modulated at a momentum corresponding to the magnetic length . for the parameters that we have investigated , the general features of the structure factor is independent of the size of the system . we obtain the excitation spectrum shown in fig . [ 3atoms - structure ] ( b ) similar to ref .   by the feynman - bijl approximation . in the continuum limit ( @xmath229 ) , we assume that the low - lying excitations are long - wavelength density oscillations and their wave functions can be approximated to have the form @xmath230 . therefore , the variational estimate for the excitation energy is @xmath231 . at zero momentum and at the momentum corresponding to the magnetic length , there are gaps , and we also observe a deviation from the free particle spectrum similar to the magneoroton case as a reminiscent of the wigner crystal . it should be noted that the deviation does not depend on the size or the number of particles in the system . as clearly seen in the fig . [ 3atoms - structure ] the energy spectrum and structure factor deviate from those of free particles , therefore , it could be used as an experimental probe of the system .    the structure factor and excitation spectrum imply some general features that are very different from that of the mott - insulator and superfluid states , and can be used a powerful experimental indication of the quantum hall states . recently , there has been several proposals for producing artificial magnetic field for neutral atoms in optical lattices @xcite , however , the implementation of each of them is still experimentally demanding . recently , there has been an experimental demonstration of a rotating optical lattice @xcite which is equivalent to an effective magnetic field ( see below ) . the technique used in this experiment , however , generates a lattice with a large lattice spacing , because it uses laser beam which are not exactly counter propagating . this longer spacing reduces the energy scale in the lattice and thereby also reduces quantities such as energy gaps . here , we shall now introduce an alternative method for generating a rotating optical lattice , which does not result in an increased lattice spacing . this method consists of rotating the optical lattice by manipulating laser beams .    in a frame of reference rotating with angular velocity @xmath232 around the @xmath233 , the hamiltonian for a particle of mass @xmath8 in an ( planar ) harmonic trap of natural frequency @xmath234 is    @xmath235    at resonance @xmath236 the form is equivalent to the hamiltonian of a particle of charge @xmath66 experiencing an effective magnetic field @xmath237 . therefore , by simply rotating the optical lattice , we can mimic the magnetic field for neutral atoms .    to rotate the optical lattice , we propose to set up four acousto - optic modulators ( aom ) and four focusing composite lenses as shown in the fig . [ lattice_rotation ] . by sweeping the acoustic wave frequency , the beams can be focused and make a rotating optical lattice .        in an aom , for the first order diffracted light we have @xmath238 , where @xmath239 is the wavelength of sound in the medium , @xmath240 is the light wavelength and @xmath241 is half of the angle between a diffracted beam and the non - diffracted beam ( fig . [ lattice_rotation ] ) . by increasing the frequency of the acoustic wave , the diffraction angle increases . however , the beam should be focused by a large aperture lens so that it always passes the region where we want to make the optical lattice . by focusing a similar but counter - propagating beam , we can make a rotating standing wave ( fig . [ lattice_rotation ] ) . by repeating the same configuration in the transverse direction , we can make a rotating optical lattice . in particular , if the aom is far from the composite lenses , @xmath242 , then @xmath243 and @xmath244 where @xmath245 , where the parameters are defined in fig . [ lattice_rotation ] . if we consider a square lattice with dimensions @xmath246 , the number of magnetic flux given by rotation is n_== where @xmath247 is the atomic recoil frequency .    on the other hand , the upper limit for the magnetic field comes from the single bloch band approximation which we made in writing the hamiltonian for the optical lattice . in order for the particles to remain in the first band , the traveling lattice beams should move them adiabatically . from ref . , the adiabaticity condition for a moving lattice with an acceleration @xmath248 equal to @xmath249 at the edge , is @xmath250 , where @xmath251 is the frequency difference between the first and the second band in the lattice . this puts a limit on how large the lattice can become @xmath252 . hence , for @xmath19 with a lattice filling fraction @xmath253 and a typical recoil frequency @xmath254 khz , one can enter the regime of fractional quantum hall effect by rotating the lattice at @xmath255 hz . if a deep optical lattice is used e.g. @xmath256 , the adiabaticity condition is easily satisfied for a lattice of size @xmath257 . the experimentally challenging part will , however , likely be to mitigate the instability of the lattice caused by the thickness of the beam and aberration of the lenses at the turning points i.e. near @xmath258 . an extended study of the realization of the fractional quantum hall states in optical lattices has been presented . we showed that a hamiltonian similar to that of a magnetic field for charged particles can be constructed for neutral ultra cold atoms and molecules confined in an optical lattice . by adding an onsite interaction for the case of @xmath19 , an energy gap develops between the ground state and the first excited state , which increases linearly as @xmath47 and saturates to its value in the hardcore limit @xmath126 . we learned that the laughlin wavefunction is a reliable description of the system for low flux densities @xmath201 . for higher @xmath0 s , the lattice structure becomes more pronounced and a better description of the system can be carried out by investigating the chern number associated to the ground state manifold . the chern number indicates that the system has topological order up to some critical flux density @xmath161 , where the properties of the ground state manifold starts to change . we have also studied @xmath2 , where compared to @xmath19 , the laughlin wavefunction only describes the ground state for lower values of the flux @xmath199 . we showed that a dipole - dipole interaction can enhance the gap and stabilize the system , and therefore make the ground state more experimentally realizable . bragg spectroscopy has been studied as a potential experimental tool to diagnose theses incompressible states . characterization of the ground state by evaluating the chern number , developed in sec . [ gauge - fixing ] , can be generalized to other interesting many - body systems where the conventional overlap calculation fails to work . in particular , this method can be applied to ground states with non - abelian statistics which are appealing candidates for fault - tolerant quantum computation . we thank k. yang , m. greiner , s. girvin and j. i. cirac for fruitful discussions . this work was partially supported by the nsf career award , packard foundation , afsor and the danish natural science research council .
Abstract:

Answer: we analyze a recently proposed method to create fractional quantum hall ( fqh ) states of atoms confined in optical lattices [ a. srensen _ et al . _ , phys . rev . lett . * 94 * 086803 ( 2005 ) ] . extending the previous work , we investigate conditions under which the fqh effect can be achieved for bosons on a lattice with an effective magnetic field and finite onsite interaction . furthermore , we characterize the ground state in such systems by calculating chern numbers which can provide direct signatures of topological order and explore regimes where the characterization in terms of wavefunction overlap fails . we also discuss various issues which are relevant for the practical realization of such fqh states with ultra cold atoms in an optical lattice , including the presence of the long - range dipole interaction which can improve the energy gap and stabilize the ground state . we also investigate a new detection technique based on bragg spectroscopy to probe these system in an experimental realization .