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: ultracold atomic gases in optical lattices may be used as analog quantum simulators of condensed matter models @xcite . simulating the fermi - hubbard model is particularly interesting . this model is nearly impossible to solve using traditional techniques , yet is amenable to simulation with cold atoms . most importantly , this model is relevant for understanding a wide class of strongly correlated electron systems : for example , many believe that it captures the essential physics responsible for high - temperature superconductivity in cuprates @xcite . the major technical challenges to this program are developing improved cooling and measurement techniques . here we explore how light scattering can be used to detect  magnetic " order in a gas of fermionic atoms in an optical lattices , focussing on the antiferromagnetic order found at low temperatures in the insulating phase of the hubbard model . producing ( and detecting ) such an antiferromagnetic state is a key step on the road to exploring superfluidity in the hubbard model . the hubbard hamiltonian for a two - spin - component fermi gas is @xmath1 where @xmath2 label lattice sites , angle brackets denote sums over nearest neighbors , @xmath3 is the spin label , and the occupation of site @xmath4 is @xmath5 . the coefficient @xmath6 is the hopping energy between adjacent sites , and @xmath7 is the on - site interaction energy . jaksch et al . @xcite show how a gas of atoms in an optical lattice reduces to this model at low temparatures , and provided expressions for @xmath6 and @xmath7 in terms of atomic properties . we consider the case @xmath8 .    at half - filling ( one atom per lattice site ) and strong repulsion ( @xmath9 ) , the hamiltonian reduces to an antiferromagnetic heisenberg model : @xmath10 with @xmath11 and @xmath12 is the spin operator for site @xmath4 . at lower lattice depths , @xmath13 is renormalized to a slightly lower value by the direct interaction between atoms on neighboring sites @xcite . an infinite system described by @xmath14 undergoes a second order phase transition to an antiferromagnetically ordered nel state at a temperature @xmath15 . achieving these temperatures is far from trivial . experiments on cold fermions @xcite have reached sufficiently low temperatures to observe characteristics of the mott insulator , but so far @xmath16 . much lower temperatures will be reached in the near future , when the next generation of cooling protocols are implemented @xcite . many of these approaches are interesting not only for their utility , but also for the insight they provide about fundamental issues of thermodynamics and quantum statistics . the proposals for detecting antiferromagnetic ordering in an atomic gas mostly rely on the fact that appropriately tuned light couples differently to the two spin components @xcite . for example , in an experiment where one can directly image individual lattice sites @xcite , one can use this selectivity to directly resolve the antiferromagnetic ordering . an alternative proposal involves measuring the spatial noise correlations in images of the density profile after turning off the trapping and lattice potentials @xcite . the experimental noise contains information about the density - density correlation functions of the system , and in principle encodes the antiferromagnetism . this method has been used to detect imposed lattice order in a bosonic mott insulator @xcite and a fermionic band insulator @xcite . its disadvantage is sensitivity to technical noise and that averaging over many experimental shots under similar conditions is needed to obtain sufficient statistics . here we propose detecting magnetic order through _ in - situ _ bragg diffraction of _ light _ ( as opposed to bragg diffraction of _ atoms _ , another common technique @xcite ) . this method is analogous to neutron scattering @xcite and magnetic x - ray scattering @xcite , as used in solid state physics . optical bragg diffraction has been used previously in cold atomic gases to confirm the crystalline ordering of nondegenerate atoms in optical lattices @xcite . compared to the other proposed techniques , the primary advantages of bragg diffraction are its simplicity , speed , and relatively large signal . in this work , we report a detailed theoretical analysis of bragg scattering from an array of atoms trapped in an optical lattice and demonstrate its usefulness under typical experimental conditions . we consider atoms confined in a three dimensional simple cubic optical lattice potential , @xmath17 , $ ] where the lattice constant @xmath18 is half of the lattice laser wavelength and @xmath19 is the depth of the optical potential , assumed to be large enough that the tight - binding approximation holds . although the transverse gaussian intensity profile of the laser beams causes @xmath19 to be position dependent , in the relevant part of the cloud it can be taken to be constant . two collisionally closed atomic hyperfine sublevels ( @xmath20 , @xmath21 ) may be treated as a pseudospin-@xmath22 system , where there is no mechanism for spin relaxation . we assume that the two states are separated in energy by a splitting @xmath23 , and have optical resonance transitions suitable for imaging with line width @xmath24  concrete examples using @xmath25li will be given later . an external magnetic field is applied to tune the interactions between atoms via a feshbach resonance . this field also defines the quantization axis for the zeeman sublevels of the hyperfine states and , hence , of the pseudo - spin states . we initially consider bragg scattering in the limit of low probe intensity and low optical density . in this limit the born approximation may be used . the more general case will be considered in sec . [ sec : multscatt ] and the appendix . our system consists of atoms in an optical lattice described by operators for position @xmath26 and occupation number @xmath27 . in the born approximation , the total amplitude for elastic scattering reduces to @xcite @xmath28 where the difference between incoming and outgoing wavevectors of the photons is @xmath29 , and @xmath30 is the scattering amplitude for individual spin-@xmath3 atoms . here , @xmath31 is the polarization vector of the incoming photon with wavevector @xmath32 and polarization @xmath33 ( generally elliptical ) , and @xmath34 is the polarization which couples to the resonant optical transition . the phase shifts @xmath35 are related to the linewidth and detuning by @xmath36 . note that eq . is only valid when the optical density of the sample is small compared to one . optically dense samples will be considered in sec . [ sec : multscatt ] . to achieve spin dependent scattering the probe light frequency must be near an atomic resonance . hence the interaction between light and the atoms is neither purely diffractive nor purely absorbtive . this is clearly illustrated by considering the lattice free case and assuming the atoms are uniformly distributed with densities @xmath37 and @xmath38 . the electric susceptibility @xmath39 of the bulk sample is then @xcite , @xmath40 where @xmath41 is the angular frequency of the light , @xmath42 is the speed of light in vacuum . the real and imaginary components of the index of refraction @xmath43 and @xmath44 are then given by the usual relationship : @xmath45 . the dependence of these quantities on the probe laser frequency for @xmath25li is show in fig . [ fig : ref ] in the case where @xmath46 . ( color online ) real ( @xmath47 , solid red line ) and imaginary ( @xmath44 , dashed blue line ) parts of the index of refraction for a uniform gas of @xmath25li atoms as a function of laser detuning . the atoms are taken to be in the two lowest hyperfine states with densities @xmath48=3.3\times 10^{12}\,\mathrm{cm^{-3}}$ ] in an external magnetic field of @xmath49 , yielding a zeeman splitting of @xmath50 . the laser excites the @xmath51 line with linewidth @xmath52 . thin vertical lines denote the two resonances and their midpoint . ]    as seen in the figure , when the probe detuning is halfway between the levels , @xmath43 is unity and @xmath44 has a local minimum . for this detuning the scattering phase shifts of the two states obey @xmath53 . this will be the operating point of our proposed measurement . there are three advantages to working at this point : ( 1 ) we can neglect bulk dispersive effects because @xmath54 , ( 2 ) absorption is small , and ( 3 ) we are maximally sensitive to the spin order . experiments measure the differential scattering cross - section @xmath55\\ & = & \frac{9}{4 k^2 } \sum_{\lambda_f}\left|(\mathbf{e}^{*}_{\mathbf{k}_f,\lambda_f } \cdot \mathbf{e}_m ) ( \mathbf{e}^{*}_m \cdot \mathbf{e}_{\mathbf{k}_i,\lambda_i}\!)\right|^2 \\ & & \times \sum_{\sigma,\sigma',j , j ' } \big [ \left\langle \hat{n}_{j\sigma } \hat{n}_{j'\sigma ' } e^{i \mathbf{k}(\hat{\mathbf{r}}_j-\hat{\mathbf{r}}_{j ' } ) } \right\rangle \\ & & \times e^{i ( \delta_{\sigma}-\delta_{\sigma ' } ) } \sin(\delta_{\sigma } ) \sin(\delta_{\sigma'})\big].\end{aligned}\ ] ] the density matrix @xmath56 describes the initial state of the system ( where @xmath57 is the thermal density matrix of the atoms ) . tracing over outgoing polarizations gives a factor @xmath58 , where @xmath59 @xmath60 is the angle between the magnetization axis and the incoming ( outgoing ) wavevector , assuming the incoming elliptical polarization is chosen to maximize the signal . onto the 2d space @xmath61 . defining the unit vectors @xmath62 and @xmath63 , the optimal incoming polarization is @xmath64 . ] the thermal average @xmath65 factorizes with @xmath66 , where we assume half - filling @xmath67 , @xmath68 , and @xmath69 ) . the factor @xmath70 , where @xmath71 is the debye - waller factor due to the zero - point motion of the atoms around the lattice sites @xmath72 in terms of the harmonic oscillator length @xmath73 of the individual wells . the differential cross - section becomes @xmath74 where we introduce the crystal structure factor and spin structure factor @xmath75 . the coefficients are @xmath76 and @xmath77 with @xmath78 . as described previously , for a laser frequency where @xmath79 ( this is achieved when @xmath80 ) , the coefficients simplify : @xmath81    the factor @xmath82 in eq . represents the well - known result that ( spin - selective ) bragg scattering measures the spin structure factor . in solids one uses neutron scattering for the same purpose and the coupling between the magnetic moment of the neutrons and electron spin causes different scattering amplitudes for @xmath83 electrons @xcite . in cold atoms the situation is even more favorable as this result shows : by choosing detuning with opposite phase shifts for the two spin states one can basically turn off the signal from the underlying lattice and probe only the spin order ( i.e.  @xmath84 ) . the predicted locations of bragg diffraction peaks ( i.e.  maxima in @xmath85 ) are given by the laue condition : the difference in wavevector between the outgoing and incoming scattered light ( @xmath29 ) is equal to a reciprocal lattice vector @xmath86 of the crystal @xcite . rewriting this condition in terms of incoming and outgoing unit propagation vectors @xmath87 , resp . and the wavelength of the scattered light @xmath88 : @xmath89 where the miller indices of the scattering plane are defined for a cubic lattice by @xmath90 . additionally , the magnitude of the reciprocal lattice vectors available for bragg scattering is further restricted by the the vector triangle inequality : @xmath91 where the second line assumes elastic scattering . the consequence of the above expression is that for the configurations typical of optical lattice experiments where @xmath92 only a handful of choices for @xmath93 will produce bragg scattering . the structure factors @xmath82 and @xmath94 in eq .   are sums that may be performed exactly for a simple cubic lattice with length per side of @xmath95 sites ( @xmath96 total sites ) . by recognizing that the sum in the crystal structure factor is that of a geometric sequence in this configuration , it evaluates to @xmath97 if @xmath93 corresponds to a reciprocal lattice vector of the crystal then the sum is maximized with @xmath98 , demonstrating the coherent enhancement of the scattering when the laue condition is satisfied . as @xmath93 moves away from a reciprocal lattice vector , @xmath94 falls off approximately quadratically . the spin structure factor @xmath82 may be evaluated in a similar way . magnetically ordered phases may be described by a spin ordering vector @xmath99 such that the average staggered spin is @xmath100 , having @xmath101-component @xmath102 . in the mean field approximation the spin dependent factor in @xmath82 is @xmath103 . in particular , a ferromagnetic ( fully polarized ) phase has @xmath104 and a nel afm has @xmath105 . with this notation @xmath106}{\sin^2[(k_j+q_j ) a /2]}.\ ] ] when @xmath107 corresponds to a reciprocal lattice vector , @xmath82 is maximized with value @xmath108 . the locations of the maxima in @xmath94 and @xmath82 are well - known results from solid state physics . as explained in standard textbooks @xcite , bragg scattering for a simple cubic lattice occurs only at planes where the miller indices are integers in the absence of magnetic ordering . if the crystal contains antiferromagnetic ordering , half - integer miller indices are also possible , indicating a doubling of the crystal unit cell . planes with mixtures of integer and half - integer miller indices have zero scattering amplitude .    in the case of a paramagnetic phase ( e.g.  when @xmath16 ) , the spins @xmath109 appearing in the sum @xmath82 do not show long range order . we can then estimate the diffuse background by assuming that all the spins point in random directions : @xmath110 , independent of @xmath93 , i.e.  this term describes incoherent scattering . for a real experiment at finite temperature , a portion of the sample is likely to be in the paramagnetic phase . the incoherent scattering from these disordered spins will yield a diffuse background . to demonstrate this technique for a typical experimental configuration , we simulate a crystal of @xmath111 atoms evaporatively cooled to fermi degeneracy in a far off - resonance optical trap , then loaded into a 3d simple - cubic lattice created from three orthogonal far red detuned laser plane - wave standing waves at @xmath112 . this experiment is currently underway at rice university . the @xmath113 atoms are prepared in an incoherent balanced mixture of the @xmath114 , @xmath115 zeeman sublevels of the @xmath116 atomic ground state ( denoted @xmath21 and @xmath20 , respectively ) . these two sublevels play a role analogous to the spin-@xmath117 electrons in crystalline solid . an externally applied uniform magnetic field controls interatomic interactions via a feshbach resonance @xcite . we choose the direction of the magnetic field to be along one of the lattice axes for experimental convenience . this field also determines the quantization axis of the atomic states . the probe light is tuned near the @xmath118 resonance of @xmath119 with wavelength @xmath120 . because @xmath88 is similar in magnitude to the lattice constant @xmath18 , only a few lattice planes will produce diffraction . applying eq . : @xmath121 the only valid nonmagnetic scattering planes for our system are @xmath122 , @xmath123 and permutations of these indices . ( the @xmath124 case represents unscattered light ) . the allowed magnetic scattering planes are @xmath125 . the basic geometry for two of the scattering planes is depicted in fig . [ fig : cube ] , with the values of the scattering angles shown . ( color online ) geometry for scattering from @xmath126 plane ( left ) and @xmath127 plane ( right ) for a cubic lattice with lattice spacing @xmath18 . the scattering planes are shaded . the thin solid line is normal to the scattering plane . the color of the lattice sites indicates spin ; this figure shows afm ordering . the arrows represent the directions of the incoming and outgoing scattered light . the labeled angles are between the outgoing light and the normal , which is the same value as the angle between the incoming light and the normal . the values of the angles are given by the bragg condition : @xmath128 , assuming lattice spacing of @xmath129 and probe wavelength of @xmath130.,title="fig:"](color online ) geometry for scattering from @xmath126 plane ( left ) and @xmath127 plane ( right ) for a cubic lattice with lattice spacing @xmath18 . the scattering planes are shaded . the thin solid line is normal to the scattering plane . the color of the lattice sites indicates spin ; this figure shows afm ordering . the arrows represent the directions of the incoming and outgoing scattered light . the labeled angles are between the outgoing light and the normal , which is the same value as the angle between the incoming light and the normal . the values of the angles are given by the bragg condition : @xmath128 , assuming lattice spacing of @xmath129 and probe wavelength of @xmath130.,title="fig : " ]    in an experiment one could probe a nonmagnetic peak such as @xmath126 to confirm the cubic ordering of the sample and a magnetic peak such as @xmath127 to detect the afm ordering . the ratio of these two peak intensities is a measure of the staggered magnetization , proportional to @xmath131 and independent of @xmath132 . the zeeman splitting between states @xmath20 and @xmath21 is @xmath133 when the external magnetic field is tuned near the feshbach resonance at @xmath49 . the bulk optical properties of the sample are those shown earlier in fig . [ fig : ref ] for one atom per site and equal spin populations .    to estimate the optical signal strength produced by the scattering we consider a detector ( ccd array or photodiode ) placed @xmath134 away from the atoms along the direction of the outgoing wavevector where we expect a peak . we then numerically evaluate eq . , assuming a perfectly ordered sample of @xmath135 atoms and a lattice depth of @xmath136 ( about 7 lattice photon recoil energies ) . the value of the spin structure factor varies depending on the phase we wish to calculate . for a fully polarized phase the spin ordering vector is chosen to be @xmath104 . for the afm phase , we use @xmath137 . finally , for the paramagnetic state we assume the spin structure factor is independent of direction with a value of @xmath138 , as mentioned previously in sec . [ sec : cubic ] . the maximum light intensity for the bragg peak is given by @xmath139 , where @xmath140 denotes the differential cross - section in direction @xmath141 . we take the probe laser to be a plane wave with intensity @xmath142 and beam waist much larger than the sample size . the total scattered power is estimated by numerically integrating over the area of the detector from the center of the peak radially to the first minimum . the results of the calculations for the various phases and the two scattering planes of interested are given in table [ tab : pred ] . for the maximal scattering configurations we find approximately @xmath143 of power in the scattered beam spread over about @xmath144 of angular divergence . using commerically available photodiodes ( with an appropriate collecting lens ) and integration times of @xmath145 should be sufficient to detect this power . for higher sensitivity , inexpensive ccd cameras may be used . and a noise floor of about 80 photons at this wavelength . ]    .[tab : pred]predicted experimental values of maximum scattering cross - section , peak intensity , and scattered power for various magnetic states ( paramagnet [ pm ] , antiferromagnet [ afm ] and polarized [ pol . ] ) and scattering planes using typical parameters , calculated as described in the text [ cols="^,^,^,^",options="header " , ]     by varying the direction of the incoming beam , we find that samples of this size are sensitive to the direction of the incoming beam on the @xmath146 scale . given the size of our simulated crystal , this is consistent with the expected diffraction limit . the previous analysis assumed that the sample was completely ordered . at finite temperature , espectially for @xmath147 near @xmath148 , the afm ordering will be incomplete . to the estimate effect of finite temperature on the bragg scattering intensity , we evaluate eqs .  , , and for a system with @xmath149 , assuming the staggered magnetization as calculated by a quantum monte - carlo calculation for a smaller system with @xmath150 . while providing an estimate of the role of thermal fluctuations , our procedure of scaling up the 216 site result to 68921 sites will overestimate the finite - size smearing of the phase boundary . we used the alps ( algorithms and libraries for physical simulations ) package to calculate this magnetization @xcite . as shown in fig .  [ fig : pint2 ] this extrapolation predicts that the @xmath151 bragg peak begins to appear as the temperature goes below @xmath152 , which is consistent with the bulk behavior of the system  @xcite . by measuring the ratio of the two bragg peak intensities we gain a measurement of @xmath102 independent of shot - to - shot fluctuations in @xmath132 . ( color online ) staggered spin per atom of the system ( red squares ) calculated using a quantum monte carlo simulation of @xmath153 atoms , and maximum scattering cross - sections of the magnetic and nonmagnetic bragg peaks ( blue circles and green triangles , resp . ) as a function of temperature . the cross - sections are calculated using eq .   with the above value of staggered magnetization for a system of @xmath154 atoms . ] for optically dense samples multiple scattering can not be neglected when calculating the differential cross - section . we calculate the effect of multiple scattering by numerically solving the lippmann - schwinger equation , which is derived in the appendix . the key result is the scattering amplitude ( also eq .  ) @xmath155 where @xmath156 are numerical coefficients found by solving eq .  . using these results , we qualitatively study the consequences of finite optical density ( od ) , in particular the interplay with the position uncertainty of the atoms . ( the debye - waller factor in eq . is an approximation of this effect @xcite . ) in this section we consider a spherical sample and assume that all lattice sites within a radius @xmath157 are occupied . position uncertainty caused by zero - point motion of the atoms is modeled as random displacements of the ( point - like ) atoms by a gaussian distribution with standard deviation @xmath158 . for the ground state of a sinusoidal well in the tight - binding regime @xmath159 , where @xmath19 is the lattice depth in units of the lattice recoil energy . for a typical experimental value of @xmath160 , @xmath161 . the lattice spacing is @xmath162 and the wavelength of the probe beam is @xmath163 . ( a  c ) : differential cross - section of the @xmath164 bragg peak for a fully polarized spherical sample with @xmath165 atoms ( @xmath166 ) , resonant probe light , and varying mean displacements ( disorder ) @xmath167 ( corresponding to ( a ) , ( b ) , ( c ) ) calculated from eqs .   and ( blue triangles ) , as a function of outgoing direction @xmath168 relative to the peak center ( directions shown in fig . [ fig : cube ] ) . for comparison we also show the born approximation result ( red circles ) . the lines are guides to the eye . ( d ) : intensity of the @xmath126 bragg peak as a function of detuning @xmath169 in the fully polarized case @xmath170 ( solid red , dashed green , dotted blue lines ) . the resonance is at @xmath171 . ( e  f ) : intensity of the @xmath127 bragg peak as a function of @xmath169 for a nel afm state with zeeman splitting between spin states of @xmath172 for varying @xmath170 ( solid red , dashed green , dotted blue lines ) and particle numbers ( e ) : @xmath165 ( @xmath173 at @xmath171 ) and ( f ) : @xmath174 ( @xmath175 at @xmath171 ) . the data in ( d  f ) have been symmetrized about @xmath176 to remove the effect of residual polarization caused by the small sample size . the vertical dotted lines in ( e , f ) indicate the atomic optical resonance . ] _ fully polarized : _ we first consider a model for a fully polarized state ( all atoms in the same pseudospin state ) with @xmath177 ( @xmath165 atoms ) . the angular distribution of the differential cross - section is shown in fig . [ fig : multscatt](a  c ) for various @xmath178 using both multiple - scattering and born approximation methods . the width of the peak is consistent with the diffraction limit for a sample of this size . as shown in fig . [ fig : multscatt](d ) , we find that for high optical density , the intensity of the @xmath179 bragg peak develops a local minimum at resonance ( zero detuning ) . on resonance the optical density is highest , and even though the optical scattering cross - section for individual atoms is maximal , the bragg signal has a local minimum because of absorption . this effect is similar to that seen in the experiment by birkl et al . @xcite , except in that case the disorder was caused by temperature . _ afm : _ the antiferromagnetic cloud is modeled by assigning lattice sites different pseudospin states according to nel order . we define @xmath169 as the laser frequency detuning relative to the center of the two atomic resonances . we assume a small zeeman splitting @xmath180 in the simulation so that the small sample of 925 atoms is optically dense when the probe laser is tuned directly between the resonances ( @xmath181 ) . this splitting is much smaller than the experimental parameters for @xmath25li near the feshbach resonance , where @xmath182 . the optical density of the simulation is comparable to what one expects with @xmath183 @xmath25li atoms , , where @xmath184 is the thickness of the sample . ] about @xmath185 times the sample size expected in the experiment . the results are shown for different radii @xmath157 ( i.e.  different particle numbers ) and disorder @xmath178 in fig . [ fig : multscatt](e , f ) . we find that the maximum intensity of the @xmath127 bragg peak scattered from an optically dense sample ( e.g.  fig . [ fig : multscatt](e ) ) is not found on one of the resonances ( as one finds in the born approximation or samples with low optical density ) , but rather around the midpoint in between the @xmath83 states ( @xmath186 ) . disorder increases incoherent scattering , reducing the amount of light coherently scattered into the bragg peaks  the disorder effectively enhances the optical density of the sample in a nonlinear way . under typical experimental conditions ( finite temperature in a harmonic trap ) the outer shell of the atomic cloud will be disordered and we expect the afm order to develop only at the center . the disordered atoms could then shield the afm from the probe light if the optical density in the shell becomes too high . the above calculation demonstrates that in this case it is favorable to detune away from the @xmath83 states to increase the penetration depth of the light and maximize the signal at the @xmath187 bragg peak . the physical picture is that for sufficiently large detuning ( but between the two resonances ) the disordered portion of the sample only contributes weak incoherent scattering , but the afm portion still exhibits detectable coherent bragg scattering from the magnetic planes . the ordered state of a classical antiferromagnet is characterized by having a nonzero staggered magnetization @xmath188 . in the absence of a symmetry breaking field , @xmath189 will point in a random direction on the bloch sphere ( see fig . [ fig : bloch ] ) , and sequential experiments will find this order pointing in a new direction . quantum mechanically the physics is similar , with the added feature that the staggered magnetization may be in a coherent superposition pointing in all directions . from a practical standpoint , this distinction is minimal : a measurement will yield a definite ( and random ) direction for @xmath189 . this ( classical or quantum ) uncertainty means that the the scattered intensity @xmath190 of a magnetic bragg peak will fluctuate from experiment to experiment . the scattered light only interacts with the @xmath20 and @xmath21 eigenstates , and hence is only sensitive to the @xmath101-component of the order @xmath102 . assuming an isotropically distributed @xmath188 , the probability distribution of the intensity will be @xmath191 , with a mean value of @xmath192 , where @xmath193 is the intensity when the spins are parallel to @xmath101 . this purely geometric effect necessitates repeating the experiment multiple times . ( color online ) cartoon of the bloch sphere representations of the various afm configurations . left : isotropic afm state , where the spins ( large red arrows ) may lie anywhere on the bloch sphere . the @xmath101-component of the staggered spin vector is shown by the small blue arrow . center : canted afm state , where the spins are restricted to be near the @xmath194-plane but lift off of it by some small canting angle ( exaggerated here for clarity ) . in this case the @xmath101-projection of the staggered spin vector is zero . right : canted afm state after applying a @xmath195-pulse . now the staggered spin has a nonzero @xmath101-component ( small blue arrow ) , which may be detected by bragg scattering . ]    recent numerical calculations point out that a sufficient polarization of the sample ( @xmath196 ) favors a staggered magnetization lying in the @xmath194-plane of the bloch sphere @xcite . this response to a spin imbalance is analogous to the response of an electronic antiferromagnet to a magnetic field and can be explained in terms of a n anisotropic spin susceptibility . as illustrated in fig . [ fig : bloch ] , it is relatively easy to add a finite polarization to @xmath194-antiferromagnetic order : one simply cants each spin slightly towards the @xmath101-axis of the bloch sphere . on the other hand , it is much more energetically costly to polarize a @xmath101-antiferromagnet : one must completely flip individual spins , adding @xmath197 energy to the system for each spin that is flipped . thus , the lowest energy configuration of a polarized afm has the nel order in the @xmath194-plane . one expects that even at finite temperature the thermodynamic ensemble will be dominated by configurations with @xmath188 lying in the @xmath194-plane . a simple model is to take the spins of the two sublattices to be @xmath198 $ ] , where @xmath199 is uniformly distributed over @xmath200 radians of a circle and @xmath201 is a small _ fixed _ canting angle relative to the @xmath194-plane , which is set by the polarization .    in the canted configuration , the @xmath101-component of the staggered magnetization is zero ( center frame of fig . [ fig : bloch ] ) , making the afm invisible to bragg scattering or optical imaging . however , the spins may be reoriented to partially project along the @xmath101-axis by an rf @xmath202-pulse at the zeeman splitting frequency , giving a non - zero bragg signal ( right panel of fig . [ fig : bloch ] ) . the resulting shot - to - shot intensity probability distribution of the bragg scattered light is then @xmath203 with a mean of @xmath204 . the intensity distributions for bragg scattering of the isotropic afm and canted afm ( after applying the @xmath202-pulse ) are shown in fig . [ fig : ipdf ] . note that the @xmath202-pulse does not change the probability distribution of the isotropic case because of symmetry . measuring these probability distributions provides a method to experimentally distinguish between the isotropic and canted afm states , particularly because of the different behavior of the two distributions when @xmath205 . ( color online ) shot - to - shot probability distributions for the intensity of the bragg scattering signal for two cases . the solid red line shows the case when the staggered magnetization may point in any direction . the dashed blue line shows the case for a canted afm , after applying a @xmath195-pulse to rotate the spins from the @xmath194-plane to the @xmath206-plane . the vertical lines indicate the mean values of the respective distributions . the vertical scale is normalized such that the areas under both curves are one . ] these calculations confirm that optical bragg diffraction is a viable method to observe antiferromagnetic ordering of atoms in an optical lattice . the experimental measurement consists of using a ccd or photodiode to measure the intensities of the peaks . the presence of the @xmath126 bragg peak indicates cubic ordering , while the @xmath127 bragg peak indicates antiferromagnetic ordering . in the actual experiment , only a portion of the atoms near the center of the crystal will be in the antiferromagnetic state . the ratio of the @xmath127 intensity to the @xmath126 intensity provides a measure of the degree of ordering . additionally , the shot - to - shot intensity distribution provides a method to determine whether the antiferromagnetic state is isotropic or canted . here we derive the full transition matrix ( t - matrix ) , including multiple scattering events . we consider scattering of a photon off of @xmath132 2-level atoms located at positions @xmath207 described by the matter - light hamiltonian @xmath208 where @xmath209 is the energy splitting between the two levels , @xmath210 is the @xmath101-component spin operator for the @xmath211-th atom in terms of the corresponding ground state @xmath212 and excited state @xmath213 , @xmath214 is the frequency of a photon with wavevector @xmath215 and polarization @xmath88 , @xmath216 and @xmath217 are creation and annihilation operators of the same photon , @xmath218 is the dipole operator ( @xmath219 is the electron charge ) and @xmath220 is the quantized electric field at the location of the atom at @xmath221 . in the rotating wave approximation ( rwa ) the matter - light interaction becomes @xcite @xmath222,\end{aligned}\ ] ] where @xmath223 , @xmath224 , @xmath225 , @xmath226 , @xmath227 and @xmath228 ( @xmath229 denotes the volume of the sample ) . in the low intensity limit ( much less than the saturation intensity of the atomic transition ) we may consider scattering of a single photon . then within the rwa the lippmann - schwinger equation for the outgoing photon @xmath230 ( where @xmath231 is an incoming plane wave photon with wavevector @xmath232 and polarization @xmath33 , @xmath233 is a positive infinitesimal ) closes with the ansatz @xmath234 and one obtains a linear equation for effective amplitudes @xmath156 ( related to @xmath235 by @xmath236 ) @xmath237 where @xmath238 is the ( renormalized ) detuning at site @xmath211 and @xmath239 is the linewidth of the transition . the dimensionless function @xmath240 has @xcite @xmath241,\\ \gamma_{jl}&=&\frac{3}{2 } \left[p \frac{\sin(k r_{jl})}{k r_{jl}}+q \left ( \frac{\cos(k r_{jl})}{(k r_{jl})^2}-\frac{\sin(k r_{jl})}{(k r_{jl})^3 } \right ) \right],\end{aligned}\ ] ] where @xmath242 and @xmath243 , is only valid when @xmath244 , which is the experimentally relevant limit . ] and @xmath245 and @xmath246 depend on whether one considers a @xmath247 or @xmath248 transition . one has @xmath249 for @xmath247 and @xmath250 for @xmath248 @xcite . the transition matrix is @xmath251 eq . includes the effects of multiple scattering to all orders , provided one can solve for @xmath156 . in the limit of low optical density of the sample , multiple scattering can be neglected and one obtains the born approximation result @xmath252 or @xmath253 we finally note that the scattering amplitude @xmath254 is related to the t - matrix through @xmath255 in the born approximation this reduces to eq .  . we thank david huse , han pu , and pedro duarte for suggestions . this work was supported under aro award w911nf-07 - 1 - 0464 with funds from the darpa ole program . additional funding was provided by nsf , onr , and the welch and keck foundations .
Abstract:

Answer: antiferromagnetism of ultracold fermions in an optical lattice can be detected by bragg diffraction of light , in analogy to the diffraction of neutrons from solid state materials . a finite sublattice magnetization will lead to a bragg peak from the @xmath0 crystal plane with an intensity depending on details of the atomic states , the frequency and polarization of the probe beam , the direction and magnitude of the sublattice magnetization , and the finite optical density of the sample . accounting for these effects we make quantitative predictions about the scattering intensity and find that with experimentally feasible parameters the signal can be readily measured with a ccd camera or a photodiode and used to detect antiferromagnetic order .


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: the physical properties of the interstellar medium ( ism ) are coupled to the star life cycles and affect the dynamics of the galaxy ; therefore , the understanding of the ism composition and phases has broad relevance . the ism has been studied by a variety of observational techniques , e.g. 21  cm radio emission @xcite , ir @xcite , and uv absorption spectroscopy @xcite . x - ray high - resolution spectroscopy has lately become an attractive method as it provides data that are unique and complementary to those gathered with other techniques . by using a bright continuum - dominated x - ray source as a backlight and grating spectrometers , e.g. the _ chandra _ high energy transmission grating spectrometer ( hetgs ) or the _ xmm - newton _ reflection grating spectrograph ( rgs ) , it is possible through the modeling of the observed absorption features to unravel key ism physical and chemical conditions . due to the penetrating power of x - rays through the gas and dust , high - resolution x - ray spectroscopy allows the direct measurement of spectral features such as absorption lines and edges that lead to the identification of multiple atomic ionization species , molecules , solid composites , and to the determination of binding energies @xcite . reliable templates are thus needed to diagnose the different absorption signatures in such a diverse and complex plasma . useful photoabsorption models to render the ism features have been progressively refined and incorporated in the x - ray spectral fitting packages . wabs , a model included in the xspec x - ray fitting package , was based on a compilation @xcite of the photoabsorption cross sections by @xcite , and included the solar abundances of @xcite ; however , it neglected contributions from ions , molecules , grains , and abundance enhancements . due to this latter limitation , @xcite energy fitted these cross sections for 17 elements in a model referred to as phabs , allowing the treatment of abundances as adjustable parameters . within the xspec spectral fitting package @xcite , tbabs considered improved cross sections and revised abundances as well as contributions from grains and the h@xmath0 molecule @xcite . in order to comply with _ chandra _ and _ xmm - newton _ spectral requirements , an updated version of this model ( tbnew ) has been available , which implements finer cross - section meshes to delineate the resonance structures of the oxygen and neon k edges and the iron l edge . it must be noted that all these models only consider photoabsorption by neutral atomic species . on the other hand , @xcite have analyzed ism absorption in low - mass x - ray binary ( lmxb ) spectra with the spex package using three absorption models that assume ionization and energy balance : hot for collisionally ionized plasmas ; slab for plasmas with arbitrary compositions ; and amol that takes into account the effects of dust and molecules . lrlcllc 4u  163653 & 105 & 1999 oct 20 & 29 & timed & @xmath1 & @xmath2 + & 1939 & 2001 mar 28 & 27 & timed & & + & 6635 & 2006 mar 22 & 23 & continuous & & + & 6636 & 2007 jul 02 & 26 & continuous & & + 4u  173544 & 704 & 2000 jun 09 & 24 & timed & @xmath3 & @xmath4 + & 6637 & 2006 aug 17 & 25 & continuous & & + & 6638 & 2007 mar 15 & 23 & continuous & & + 4u  182030 & 1021 & 2001 jul 21 & 9.6 & timed & @xmath5 & @xmath6 + & 1022 & 2001 sep 12 & 11 & timed & & + & 6633 & 2006 aug 12 & 25 & continuous & & + & 6634 & 2006 oct 20 & 26 & continuous & & + & 7032 & 2006 nov 05 & 47 & continuous & & + cygnus  x1 & 3407 & 2001 oct 28 & 17 & continuous & @xmath7 & @xmath8 + & 3724 & 2002 jul 30 & 8.8 & continuous & & + cygnus  x2 & 8170 & 2007 aug 25 & 65 & continuous & @xmath9 & @xmath10 + & 8599 & 2007 aug 23 & 59 & continuous & & + & 10881 & 2009 may 12 & 66 & continuous & & + & 1102 & 1999 sep 23 & 28 & timed & & + gx  9 + 9 & 703 & 2000 aug 22 & 20 & timed & @xmath11 & @xmath12 + & 11072 & 2010 jul 13 & 95 & timed & & + sco  x1 & 3505 & 2003 jul 21 & 16 & continuous & @xmath13 & @xmath14 + xte j1817330 & 6615 & 2006 feb 13 & 18 & continuous & @xmath15 & @xmath16 + & 6616 & 2006 feb 24 & 29 & continuous & & + & 6617 & 2006 mar 15 & 47 & continuous & & + & 6618 & 2006 may 22 & 51 & continuous & & +    @xcite have carried out a study of ism oxygen k absorption using _ chandra _ spectra of the lmxb xte  j1817330 . spectral fits were performed with xspec using a physical model referred to as warmabs , which is based on precalculated tabulations of ionic fractions and level populations from the xstar photoionization code to generate synthetic emission and transmission spectra . in this work , a model benchmark of the theoretical atomic data for oxygen ions by @xcite required wavelength shifts for both the k lines and photoionization cross sections in order to fit the observed spectra accurately . moreover , an attempt to develop a definitive absorption model for neutral oxygen with a new computed cross section revealed standing discrepancies between the laboratory and observed line positions @xcite . the warmabs model has been subsequently used with the adjusted cross sections of @xcite and @xcite to analyze ism oxygen absorption toward eight lmxbs @xcite . oxygen abundances and ionization parameters along different lines of sight were therein determined , finding only small abundance differences due to the theoretical atomic data , and estimating the @xmath17(o  ii)/@xmath17(o ) column - density ratio to be around 10% . in the case of magnesium , new photoabsorption cross sections by @xcite in warmabs were used to fit the mg  k edge in _ xmm - newton _ spectra of the lmxb gs  1826238 predicting that most of this element was in ionized form . the warmabs model includes the most accurate atomic data for atomic neutrals and ions , but it is constrained by the assumption of ionization equilibrium . this may not be practical owing to deviations in the ionizing energy distribution , other ionization mechanisms , or the superposition of multiple components along the line of sight . hence , we present a model of ism x - ray absorption , to be referred to as ismabs , that takes into account both neutrals and charged ions allowing their column densities to vary as free parameters . we benchmark ismabs by means of fits to high - resolution spectra from x - ray binary sources obtained with the _ chandra _ hetgs and comparisons with the tbnew and warmabs models . this procedure has involved a detailed study of the oxygen and neon k edges as well as the iron l edge . the outline of this paper is as follows . in section  [ sec_dat ] the data reduction process is briefly summarized , and in section  [ sec_mod ] we describe the selected atomic database and the x - ray absorption model developed to render the ism absorption features . in section  [ sec_fit_mo ] we show fit comparisons of the tbnew and ismabs models . results obtained from the broadband fit are reviewed in section  [ sec_broad ] , with particular attention to the neon k- and iron l - edge regions , followed by a detailed discussion of the fits in section  [ sec_results ] . finally , in appendix  [ appendix ] we give an account of the importance of ionization balance in abundance determinations . the most suitable sources for x - ray ism studies are bright lmxbs as they do not usually show intrinsic spectral features , and preferably with photon fluxes with high signal - to - noise ratios to enable the measurement of the ism contribution . to hone the present photoabsorption model , we have selected eight candidates from the catalog of lmxbs @xcite and the transmission grating catalog and archive ( tgcat ) such that the total number of counts in each hetgs spectrum is greater than @xmath18 . we have also included spectra from the high - mass x - ray binary cygnus  x1 due to its high photon flux , giving a total of 25 observations . observational details are listed in table  [ tab1 ] , where it is worth adding that the high - energy resolution of the _ chandra _ gratings is certainly indispensable in the present modeling developments . the spectra have been reduced with the standard ciao threads , the zero - order position being estimated in each case with the findzo algorithm . spectral fitting was carried out with the xspec package ( version 12.8.1 ) . llll h & @xcite@xmath19 & mg  i , mg  ii , mg  iii & @xcite@xmath20 + he  i & opacity project@xmath21 & si  i & @xcite + he  ii & @xcite & si  ii , si  iii&@xcite@xmath22 + c  i , c  ii , c  iii & @xcite@xmath20 & s  i & @xcite + n  i , n  ii , n  iii & @xcite@xmath22 & s  ii , s  iii & @xcite + o  i & @xcite@xmath23 & ar  i & @xcite + o  ii , o  iii & @xcite@xmath22 & ar  ii , ar  iii & @xcite@xmath22 + ne  i & @xcite@xmath24 & ca  i & @xcite + ne  ii , ne  iii & @xcite@xmath25 & ca  ii , ca  iii & @xcite + & & fe & @xcite@xmath26 +    ism photoabsorption causes significant modifications to the observed x - ray spectra , and must therefore be taken into account in any high - resolution spectroscopic analysis . for this purpose , @xcite developed the tbabs and tbnew absorption models that have been widely used in the last decade . however , they only include neutral atomic species compromising their ability to render complex ism spectral structures such as the o , ne , and mg k edges @xcite . given this limitation , we have devised a new ism x - ray absorption model  to be referred to as ismabs  that compiles the most accurate atomic data available for neutrally , singly , and doubly ionized species of h , he , c , n , o , ne , mg , si , s , ar and ca ; in the case of iron , an experimental cross section has been chosen to model the l - edge absorption region ( see section  [ db ] ) .    photoabsorption in ismabs takes the form @xmath27 where @xmath28 is the observed spectrum , @xmath29 the spectrum emitted by the source , and @xmath30 is the absorption coefficient . the optical depth @xmath31 is usually defined as @xmath32 @xmath33 being the total photoionization cross section comprising all the atoms and @xmath34 is the hydrogen column density . if a neutral gas is considered , @xmath35 may be written @xmath36 where @xmath37 is the photoionization cross section of a chemical element with atomic number @xmath38 , and @xmath39 gives its abundance relative to hydrogen . it is worth noting that , due to the abundance dependency on column density , there is a degeneracy between @xmath39 and @xmath34 , making it possible by varying one or the other to obtain multiple solutions in the calculation of @xmath31 . consequently , we prefer to adopt the following definition for @xmath31 : @xmath40 here @xmath41 and @xmath42 are respectively the cross section and column density for the @xmath43-th ion . we therefore manage the ion column densities , including @xmath34 , as model parameters in order to determine @xmath31 , thus avoiding any parametric ambivalence . ismabs always assumes that @xmath17 ( ) @xmath44()@xmath45 , and as a default , @xmath17 ( ) @xmath46 . but @xmath17 ( ) is a free parameter which the user can modify if desired . it must be emphasized that there are no changes to the cross sections during fitting except for interpolation ; i.e. they are read just once by the code when the model is invoked . the energy range covered by the model is @xmath47  kev . the redshift @xmath48 is also a model parameter , and we will include turbulent broadening in the next revision to ismabs ( with the turbulence velocity @xmath49 as the model parameter in units of kms@xmath50 ) ; however , observations point to a turbulence ( @xmath51  kms@xmath50 ) that can not be detected by _ chandra _ @xcite . furthermore , for this velocity range , the natural broadening of inner - shell resonances dominate over turbulence broadening . ismabs has been specially developed for x - ray fitting packages such as xspec @xcite , isis @xcite and sherpa @xcite . in addition to high signal - to - noise x - ray spectra from a back - light source , reliable atomic data are required to model the ism absorption signatures . the atomic databases used by many of the currently available x - ray absorption models are not complete and of sufficient accuracy to model grating - resolution spectra . this situation was illustrated in the work of @xcite , who used four high - resolution _ chandra _ spectra of the lmxb xte  j1817330 to benchmark the oxygen theoretical atomic database @xcite in the xstar photoionization code . the photoabsorption cross sections of the lowly ionized species required several adjustments : wavelength shifts of the k@xmath52 resonances in o  i ( @xmath53 m ) and o  ii ( @xmath54 m ) and of the whole cross sections for both of these ions . moreover , a new and more accurate theoretical cross section for o  i by @xcite showed significant differences in the oscillator - strength distribution , and brought to light a discrepancy in the wavelength scale , namely @xmath55  m , between the observed k@xmath52 line position and the laboratory measurement @xcite . the quantitative differences in photoabsorption modeling brought about by these two o  i cross sections , i.e. @xcite and @xcite , were examined by @xcite in a study of ism absorption along several lines of sight , constraining them to less than 5% in the average ionization parameter and oxygen abundance . a further point of interest is the presence of condensed matter in the ism , which can be detected through the analysis of x - ray absorption fine structure in the form of spectral modulations near the photoabsorption edges , and by the absence of strong inner - shell resonance transitions below the ionization threshold @xcite ; however , their proper characterization would require beforehand the mastering of the atomic absorption edges which is the central topic of the present work .      in table [ tab2 ] we list the reference sources for the cross sections incorporated in ismabs . we include the theoretical atomic photoabsorption data for  by @xcite ; and by @xcite for the singly and doubly ionized species of si , s , ar , and ca , which were calculated with the @xmath56-matrix method using an optical potential to account for auger damping . for the neutral species of the latter four elements , we adopt , in the absence of better data , the theoretical cross sections of @xcite even though they do not take into account the inner - shell resonances that give rise to the k - edge structures . these cross sections are also used for the neutral , singly , and doubly ionized species of c. for oxygen ions , we select the theoretical cross section for by @xcite and those by @xcite for and , including the corrections specified by @xcite as previously mentioned . in the case of neon , we take the @xmath56-matrix data ( including auger damping ) by @xcite for and by @xcite for and . for the magnesium ions ,  , we take the photoabsorption cross sections computed by @xcite with the @xmath56-matrix method that account for orbital relaxation effects to reduce a previous over - estimate of the k - edge position @xcite .    in the case of iron , @xcite fitted the l - edge region using the experimental cross section of metallic iron by @xcite . on the other hand , @xcite adjusted the l  edge using the atomic data for the neutral by @xcite , but included a molecular solid component since the edge can not be fitted with a pure gas model . using _ hetg spectra from five x - ray binaries , @xcite adjusted the iron l  edge in the tbnew absorption model with the experimental cross section @xcite . in the present work , we adopt the metallic iron measurements @xcite but with a wavelength correction of @xmath57  m . this shift was originally suggested by @xcite , who argued that the energy resolution of the laboratory data in combination with the hetg absolute instrumental wavelength accuracy would allow it . we also considered the widely used atomic fe l - shell cross section of @xcite , but the resulting spectral fit was inferior to that obtained with the shifted cross section of @xcite . finally , it is worth pointing out that both warmabs and ismabs use the same atomic database except for fe since warmabs uses the fe cross - section of @xcite . lclllll ne  iii & k@xmath58 & 13.453 & @xmath59 & @xmath60 ( ne  ix k@xmath52 ) & @xmath61 & @xmath62 + ne  iii & ? & & @xmath63 & & & + ne  iii & k@xmath64 & 13.709 & @xmath65 & @xmath66 & & @xmath67 + ne  ii & ? & & @xmath68 & @xmath69 ( ne  vii k@xmath52 ) & & + ne  ii & ? & & @xmath70 & & & + ne  ii & k@xmath58 & ? & @xmath71 & @xmath72 & & + ne  ii & k@xmath64 & 13.992 & @xmath73 & @xmath74 & & @xmath75 + & & & @xmath76 & @xmath77 ( ne  v k@xmath52 ) & & + ne  i & k@xmath64 & 14.298 & @xmath78 & @xmath79 & @xmath80 & @xmath81 + ne  iii & k@xmath52 & 14.516 & @xmath82 & @xmath83 & @xmath84 & @xmath85 + ne  ii & k@xmath52 & 14.605 & @xmath86 & @xmath87 & @xmath88 & @xmath89 +      following the benchmark of the oxygen cross sections by @xcite , we proceed to compare the theoretical and observed absorption lines in the ne - edge region of spectra with high signal - to - noise fluxes , namely cygnus  x2 and xte  j1817330 . figure  [ fig1 ] shows the best fits in the 1315    interval using a simple functional model , i.e. a powerlaw and several gaussian profiles . for each source , all observations are fitted simultaneously , but the spectra are combined for illustrative purposes ; the black data points correspond to the observations while the red solid lines represent the model best fit for each case . the bottom panels show fit residuals in units of @xmath90 . several absorption lines are clearly observed including those due to the k@xmath52 transitions of and . identifications of the spectral features and line wavelength positions are listed in table  [ tab3 ] where the estimates by @xcite , @xcite , and @xcite are also tabulated for comparison . we have identified absorption lines due to transitions in , , and , and the line at @xmath76   presumably belongs to the k@xmath52 transition of . a comparison between the theoretical resonance positions with those determined from the spectral fits reveal the inaccuracies of the theoretical atomic data . the left panel of figure  [ fig2 ] shows with solid lines the theoretical cross sections for  @xcite that are included in the ismabs model . the vertical dashed lines are placed at the positions of the absorption features determined from the functional fits . in table [ tab3 ] we give the line assignments based on a comparison between our functional fits and the theoretical atomic data . identifications made in previous work using astrophysical data are also listed : cygnus  x2 was analyzed by @xcite and @xcite while xte  j1817330 was included by @xcite in their analysis of 11 lmxbs . it may be noted that the k@xmath52 position as well as those with @xmath91 are displaced with respect to the observations , and the cross sections of and also show similar displacements . in this respect , @xcite have previously indicated a similar wavelength offset using high - resolution _ chandra _ spectra . in what follows we have adjusted the k@xmath52 resonance positions , and shifted the cross sections for the continuum and higher resonances for these three species in order to obtain the best possible agreement with the observed lines . the shifts are , respectively , @xmath92 m  and @xmath93 m  for the ne  ii and ne  iii k@xmath52 resonances ; and @xmath94 m  , @xmath95 m , and @xmath96 m  for the ne  i , ne  ii and ne  iii cross sections . the right panel of figure  [ fig2 ] shows the cross sections after the application of the wavelength corrections to be finally included in the ismabs model . in this section we compare two different x - ray absorption models , namely tbnew and ismabs , the former including photoionization cross sections for elements with atomic number @xmath97 @xcite and the contribution from molecules and interstellar grains . it therefore favors the modeling of the observed extinction and depletions in the diffuse ism . model parameters are the hydrogen column density ( @xmath34 in units of @xmath98 ) , elemental abundances ( @xmath99 relative to solar ) , and the redshift ( @xmath48 ) . however , tbnew only takes into account neutral atomic species , and therefore , the warm or ionized ism phases can not be adequately modeled . the ismabs model has been described in detail in section  [ sec_mod ] . lcccccc @xmath34 & @xmath100 & @xmath101 + @xmath102 & @xmath103 & @xmath104 + @xmath105 & @xmath106 & @xmath107 + @xmath108 & @xmath109 & @xmath110 + reduced @xmath111 & @xmath112 & @xmath113 +    figure [ fig3 ] shows the ne k- , fe l- , and o k - edge regions in the broadband fit of xte  j1817330 with both tbnew and ismabs . in each panel , the black data points correspond to the observations while the solid lines represent the best fits with tbnew ( blue line ) and ismabs ( red line ) . the bottom panels depict fit residuals in units of @xmath90 . it is clear that the ismabs residuals are smaller particulary around the neon and oxygen edges . in the fe l - edge region the residuals are similar since both models use the same experimental atomic data , although small differences are noted perhaps due to the continuum . it is found that the inclusion of charged ions in addition to neutrals leads to a significant improvement in the fits . the difference in the statistics is of @xmath114 if the fit is performed with ismabs , with 4 new free parameters . moreover , as shown in table  [ tab4 ] , if only the column densities for neutrals are considered ( i.e. considering the same number of free parameters ) , both models give similar results , however the ismabs model performs a better fit , with @xmath115 , due to the higher accuracy in the atomic data included . following @xcite , we have performed a broadband fit ( 1124  ) of all the spectra listed in table  [ tab1 ] using ismabs and a power - law model . we rebin the data to at least 20 counts per channel in order to use @xmath111 statistics . for each source , all observations are fitted simultaneously varying the power - law normalization and photon index , thus accounting for possible continuum variations between observations that are always detected . the h , o , ne , and fe column densities , as well as those of o and ne ions , are taken as free parameters in ismabs . figures  [ fig4 ] and  [ fig5 ] show the best ismabs fits in the 1124    interval . in the case of sco  x1 the fit was performed in the 1524    interval as there are no counts below @xmath116   . for each source the observations have been combined for illustrative purposes . it must be noted that , even though this figure looks similar to figure  1 of @xcite , it is now obtained with our new photoabsorption model . table  [ tab5 ] lists the fit results where it may be appreciated that ismabs yields the ion column densities directly , and a statistically acceptable fit is obtained for each source . a comparison of the present hydrogen column densities with previous results is given in table  [ tab6 ] finding satisfactory agreement , i.e. within the error bars . as mentioned in @xcite , the column - density variations around cygnus  x1 are not resolved in the 21  cm survey ; therefore , the values derived by means of spectral fits should be more reliable . a zoom of the broadband fit in the interval 1315    is depicted in figures  [ fig6 ] and  [ fig7 ] ( sco  x1 is excluded for lack of counts below @xmath116   ) . it may be seen that the absorption features in this wavelength region are more complex than a simple edge : k@xmath52 absorption lines of and are detected at @xmath117   and @xmath118  , respectively , as well as k@xmath64 from , , and ( @xmath119  , @xmath120  , and @xmath121  , respectively ) . the line at @xmath122  , attributed to k@xmath58 , is observed in all sources . therefore , the occurrence of such lines justifies the inclusion of accurate cross sections for and in ismabs in order to procure reliable spectral fits ; for instance , in figures  [ fig6 ] and  [ fig7 ] low residuals are found around the k@xmath64 and k@xmath58 lines at @xmath120   and @xmath123  , respectively . furthermore , a similar zoom is shown in figures  [ fig8 ] and  [ fig9 ] of the fe  l edge ( 1718.5   ) which benchmarks the cross section adopted in ismabs ( see section  [ db ] ) . it is worth recalling that the selected cross section for was obtained from measurements of solid iron , and it seems to fit the observations adequately . this finding supports previous assertions that most of the fe in the ism is in solid form @xcite . we have developed a new spectral model , referred to as ismabs , to fit ism x - ray absorption that takes into account photoabsorption cross sections for neutrally , singly , and doubly ionized species of cosmically abundant elements , namely h , he , c , n , o , ne , mg , si , s , ar , ca , and fe . in other words , ismabs enables the determination of the ion column densities directly . although the predominant component of the ism is the neutral gas , the inclusion of singly and doubly ionized species provides improved fits when compared with a model ( e.g. tbnew ) that only accounts for neutral absorption . this shows that the understanding of the relevant atomic systems and processes , as well as the accuracy and completeness of the relevant atomic data , are both crucial to avoid incorrect interpretations of the absorption structures detected in x - ray high resolution spectra . this may occur in studies of ism x - ray absorption that rely on fits with gaussian profiles . llllll 4u  163653 & @xmath124 & 3.30 & @xmath125 & @xmath126 + 4u  173544 & @xmath127 & 3.03 & @xmath128 & @xmath129 + 4u  182030 & @xmath130 & 1.52 & @xmath131 & @xmath132 + gx  9 + 9 & @xmath133 & 2.10 & @xmath134 & + cygnus  x1 & @xmath135 & 8.10 & @xmath136 & @xmath137 @xmath23 + & & & & @xmath138 + cygnus  x2 & @xmath139 & 2.20 & @xmath140 & @xmath141 + sco  x1 & @xmath142 & 1.47 & & @xmath143 + & & & & @xmath144 + xte  j1817330 & @xmath145 & 1.58 & & +    figure  [ fig10 ] shows for all sources a comparison of the oxygen , neon , and iron total column densities obtained from the ismabs broadband fit with those derived by @xcite . for an atomic element @xmath146 , say , @xmath147 . the respective column densities by @xcite were obtained through a fit of eight galactic sources which includes those of the present report except for sco  x1 and xte  j1817330 . the average values for the oxygen column densities are @xmath148  @xmath98 ( ismabs ) and @xmath149  @xmath98 @xcite . in the case of neon we obtained an average value @xmath150  @xmath98 ( ismabs ) compared to @xmath151  @xmath98 @xcite . lastly , for iron we estimated an average value of @xmath152  @xmath98 ( ismabs ) compared to @xmath153  @xmath98 by @xcite .    from such data we have estimated the @xmath154 and @xmath155 column - density ratios derived from the ismabs broadband fit . neon is chosen as a reference since , as a noble gas , it is neither depleted to molecular nor to solid forms . the present average ratio is @xmath156 , while that derived by @xcite is @xmath157 . ratio differences may by due to systematic effects such as line saturation as quoted by @xcite for the case of xte  j1817330 . the ratio derived by @xcite is higher than the present value ; however , as discussed by these authors , the inclusion of black - hole systems ( e.g. cygnus  x1 ) may modify results due to the presence of intrinsic emission or absorption features . we also obtain the average ratio @xmath158 which is lower than solar , namely @xmath159 @xcite and @xmath160 @xcite . this discrepancy may be due to either ism fe depletion or ne overabundance . in fact , as previously discussed , we have found a ne abundance greater than solar in the lines of sight toward 4u  163653 , 4u  173544 , gx  9 + 9 , and xte  j1817330 . the @xmath155 ratio by @xcite ( @xmath161 ) , on the other hand , is in a fair agreement with ismabs . a further comparison of the iron and oxygen column densities finds the ismabs ratio @xmath162 in close agreement with the solar ratio of @xcite ( @xmath163 ) but somewhat lower than that by @xcite ( @xmath164 ) . future analysis of the fe k - edge absorption region would indeed help to constrain the ism fe abundance fractions in order to settle this issue . the high - resolution x - ray spectra toward eight lmxbs have been analyzed . this study involved a broadband fit with ismabs to constrain the hydrogen column density , finding good agreement with previous work . then , a simultaneous fit of the neon- and iron - edge regions was performed for each source , where the resulting column densities indicate the presence of a lowly ionized gas along the different lines of sight . we intend to carry out in the near future a more extensive analysis of ism photoabsorption using all the grating spectra from x - ray binaries available from the _ xmm - newton _ and _ chandra _ archives . our target is to create a library of column densities and abundances for all the detectable atomic species in these spectra . we believe that once the atomic k edges are fully analyzed and modeled , the search for molecular and dust spectral features can be pursued in earnest , a topic of current interest in ism astrophysics . finally , the present ismabs model can be downloaded for use in x - ray data analysis packages such as xspec @xcite , isis @xcite and sherpa @xcite . this work was supported by grant 15400673 of the chandra theory program . lllllllllc 4u  163653 & @xmath165 & @xmath166 & @xmath167&@xmath168 & @xmath169 & @xmath170&@xmath171&@xmath172 & @xmath173 + 4u  173544 & @xmath174 & @xmath175 & @xmath176&@xmath177 & @xmath178 & @xmath179&@xmath180&@xmath181 & @xmath182 + 4u  182030 & @xmath183 & @xmath184 & @xmath185&@xmath186 & @xmath187 & @xmath188&@xmath189&@xmath190 & @xmath191 + gx  9 + 9 & @xmath192 & @xmath193 & @xmath194&@xmath195 & @xmath196 & @xmath197&@xmath198&@xmath199 & @xmath200 + cygnus  x1 & @xmath201 & @xmath202 & @xmath203&@xmath204 & @xmath205 & @xmath206&@xmath207&@xmath208 & @xmath209 + cygnus  x2 & @xmath210 & @xmath211 & @xmath212&@xmath213 & @xmath214 & @xmath215&@xmath216&@xmath217 & @xmath218 + sco  x1&@xmath219 & @xmath220 & @xmath221&@xmath222 & ( fixed ) & ( fixed ) & ( fixed)&@xmath223 & @xmath224 + xte  j1817330&@xmath225 & @xmath226 & @xmath227&@xmath228 & @xmath229 & @xmath230&@xmath231&@xmath232 & @xmath233 + here we discuss the differences in fitting spectra with the ismabs and warmabs models , the latter as implemented in previous work @xcite . an important complication affecting our modeling concerns the smooth background photoelectric absorption which can not be identified uniquely by spectral features . this background is always present , and arises from hydrogen , helium , and valence - shell absorption by the same ions responsible for the edges and resonance absorption lines . since all these mechanisms produce a smooth opacity , it is impossible to distinguish them uniquely . if the gas is assumed to be in ionization balance , then the relative columns of neutral h , , and are determined ; if not , it is important to be aware of the degeneracy between the smooth absorption produced by h , , and . this degeneracy is illustrated graphically in figure  [ fig11 ] , which shows the @xmath111 contours of the fit to the xte  j1817330 spectrum for various combinations of the h and column densities . the solid curves show contours of constant @xmath111 spaced at @xmath234 . this clearly illustrates that the combinations of these variables lying along the diagonal are equivalent while those that fit the data have an almost exact inverse relationship . this obvious behavior must be understood and accounted for when deriving abundances relative to h or he in x - ray spectral fitting . in the results reported in the main body of this paper , which are based on ismabs , we have assumed that all the smooth opacity in the x - ray band comes from h + and that the column densities are negligible . the warmabs model includes the results of a photoionization equilibrium calculation , and therefore , the h and column densities are tied together . the value taken from the warmabs fit of the xte  j1817330 spectrum is shown with a dot in figure  [ fig11 ] . this corresponds to an ionization parameter of @xmath235 for which the column density of neutral h + is 10.8 times that of ; i.e. , the equilibrium assumption removes the degeneracy between the two variables . a quantity of interest associated to the fitting of interstellar x - ray absorption spectra is the abundance of elements such as oxygen and neon relative to hydrogen . if the smooth background absorption in the x - ray band were solely due to neutral h , then the h and o or ne column densities derived from the fits could be used to derive the relevant abundance ratios by simply taking column - density ratios . the degeneracy between and h complicates this procedure since the total h column density can not be uniquely determined . this is shown graphically in figure  [ fig11 ] with the vertical dashed curves which are curves of constant o abundance relative to h for the fit to the spectra of xte  j1817330 . the curves are labeled according to the conventional notation : @xmath236=0 $ ] corresponds to the solar oxygen abundance @xcite , @xmath236=1 $ ] to @xmath237 solar , etc . this shows that , although acceptable fits can be obtained for any value of @xmath34 in the range up to @xmath238  @xmath98 , the inferred abundance of o relative to h can take any value from 0.75 solar to 10 times solar .    finally , we mention the unavoidable limitation that different numerical model implementations , even if they are based on the same atomic data , inevitably differ due to numerical details . for example , calculations in ismabs use cross sections tabulated on a grid that is finer than most instrumental grids , but which inevitably requires interpolation or smoothing to produce a model for xspec . warmabs , on the other hand , adopts such tabulations for some cross sections while for others analytic expressions are employed . examples of the latter are the all - important cross sections for neutral o , ne , and mg of @xcite . also , is important to note that the fe cross - sections are different between both models , ismabs uses @xcite cross - section while warmabs uses @xcite cross section . small numerical differences in the cross section or opacity are exponentially amplified in the transmission coefficient which is most closely related to the observed spectrum . thus , in spite of efforts to minimize numerical differences between codes , they are unavoidable at some level and it is useful to attempt to quantify them . to summarize , we identify the three principal issues that affect the quantitative results of fitting to interstellar absorption and model comparisons : ( i ) the importance of ions and whether they are treated in equilibrium ratios or not ; ( ii ) the ambiguity associated to the smooth background ; and ( iii ) numerical implementation issues . we now proceed to illustrate them with quantitative examples .    in figure  [ fig12 ] we show a comparison of the ion column densities derived from fitting to xte  j1817330 with warmabs and ismabs . in the case of warmabs , we adopt two ionization components with @xmath239 and @xmath235 , both with a column density of 10@xmath240  @xmath98 . the elemental abundances for the two components are allowed to vary freely in order to obtain the best fit . the four observations of this source are fitted simultaneously , and all the absorption parameters ( i.e. the warmabs ionization parameters , column densities , elemental abundances , turbulent velocities , and redshifts ) are forced to be the same for all the different observations . the continuum is taken to be the sum of a black body and a power law , and the parameters of these components are allowed to differ between the various observations to reflect source variability . the best fit has a @xmath241 . warmabs can output a list of ion column densities ( even though these are not free parameters ) , and in figure  [ fig12 ] we show a comparison of the warmabs ion column densities with those from an ismabs fit . in the latter the h and column densities are fixed at the values taken from the warmabs fit . we take the column densities from all the ions from the warmabs best fit and put them into ismabs . this includes the h , , and column densities . using the same continuum parameters we obtain @xmath242 . this @xmath111 difference is likely to be numerical since all other parameters of the two fits are identical . we then allow the ion column densities in ismabs and the continuum parameters to vary freely in order to obtain the best fit . in doing so we do not allow the h , , and column densities to vary , and we preserve the restriction that the absorption model is the same for all observations . this procedure produces a fit with @xmath243 . figure  [ fig12 ] shows the ratio of the column densities derived from these two different methods . this represents both an astrophysical experiment on whether the ratios show photoionization equilibrium , and in some sense , a test of model accuracy ( note that the si and s column densities are negligible in both fits , so they are omitted ) . the two fits agree to within the error bars for most ions , a notable exception being for which the warmabs model predicts a significantly lower column density than ismabs . this indicates the existence of a component in the absorbing gas toward this source which is more highly ionized than the equilibrium gas that can actually fit most neutrals and singly ionized species . oxygen shows such an effect more clearly than other elements due to its relatively high abundance . a further step in exploring the effects of ionization equilibrium in absorption fitting is to sum the ion column densities derived from the fits with and without equilibrium ( i.e. using warmabs or ismabs ) to get elemental abundances . results are shown in figure  [ fig13 ] where the ratio of elemental abundances are obtained by summing up the ismabs ion column densities to the elemental abundances obtained from the warmabs column densities . the abundances shown here for ismabs make use of the h and column densities derived by warmabs . also shown are the ratio of the ismabs abundances to the values in the compilation by @xcite . this shows that , in spite of the disagreement for , both fitting methods yield abundances that are consistent with each other and with the compilation . however , the ne abundance in this source is an exception ; ismabs manifests a tendency to produce larger abundances since it is capable of accounting for ion absorption that is not tied up to the ionization equilibrium assumption .
Abstract:

Answer: we present an x - ray absorption model for the interstellar medium , to be referred to as ismabs , that takes into account both neutral and ionized species of cosmically abundant elements , and includes the most accurate atomic data available . using high - resolution spectra from eight x - ray binaries obtained with the _ chandra _ high energy transmission grating spectrometer , we proceed to benchmark the atomic data in the model particularly in the neon k - edge region . compared with previous photoabsorption models , which solely rely on neutral species , the inclusion of ions leads to improved spectral fits . fit parameters comprise the column densities of abundant contributors that allow direct estimates of the ionization states . ismabs is provided in the appropriate format to be implemented in widely used x - ray spectral fitting packages such as xspec , isis and sherpa .