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: several different probes have been used to investigate the properties of the two - dimensional electron gas ( 2deg ) in the quantum hall regime , transport and optical experiments , or equilibrium methods including capacitance and magnetization measurements , to name only few . the magnetization of high mobility homogeneous 2deg has been measured by two methods . one method uses sensitive mechanical , torque magnetometers.@xcite more recently , a low - noise superconducting quantum interference device ( squid ) has also been used.@xcite these precision measurements reveal many - body effects such as the exchange enhancement at odd filling factors @xcite and the fractional quantum hall effect and , in addition , an unidentified effect around filling factor @xmath0 that might be connected to skyrmions.@xcite    these experimental techniques are expected to be used soon for measuring the magnetization in lateral superlattices . the magnetization has been calculated for a disordered homogeneous 2deg within the hartree - fock approximation ( hfa)@xcite and within a statistical model for inhomogeneities corresponding to a hartree approximation ( ha ) . @xcite the results show strong exchange effects , already observed in experiments , @xcite and manifestations of the screening properties of the 2deg .    concerning periodic systems , the magnetization and the persistent current have been calculated by kotlyar et . for a finite array of quantum dots using a mott - hubbard model for the electron - electron interactions , both intra - dot and inter - dot.@xcite in an infinite lateral superlattice , defined by a potential periodic in two directions ( electric modulation ) the energy spectrum in the presence of a magnetic field can only be calculated when the ratio of the magnetic flux through a lattice cell to the unit flux quantum , @xmath1 , is a rational number . the unit cell can then be enlarged to have an integer number of flux quanta flowing through it.@xcite in principle , the magnetization of a 2deg in a lateral superlattice can be evaluated in the thermodynamic limit as the negative derivative of the free energy with respect to the magnetic field @xmath2,@xcite since the energy spectrum is always a continuous function of the flux . @xcite the inclusion of the coulomb interaction to the model , within a self - consistent scheme such as the ha @xcite or the hfa,@xcite does severely limit the possibility to effectively change the value of @xmath2 by a small amount in order to take a numerical derivative . nevertheless , such a thermodynamic calculation can be done for a modulation that varies only along one spatial direction , i.  e.  for an array of quantum wires .    in this paper we shall evaluate the magnetization of a 2deg in a periodic potential corresponding to a weak density modulation in both , or in only one spatial direction . in other words our system is either an array of dots or antidots , or of parallel quantum wires , in most cases with a strong overlap . we first consider a finite system with boundaries , and then the unbound system .    for the finite system we are able to calculate the total magnetization . for the infinite system with a 2d potential we shall rather use the definition applicable to a mesoscopic system with a phase coherence length @xmath3 much larger than the spatial period of the square unit cell @xmath4 . it is clear to us that in this manner we are calculating the contribution to the magnetization due to the periodic modulation , neglecting the contribution stemming from the edge of a real system . in an experiment the total magnetization is indeed measured and we would have a hard time arguing that the edge contribution is statistically insignificant in the thermodynamic limit . our way out of this dilemma is to compare several results that we can obtain . first , we compare the results for the finite system of various sizes , by heuristically separating the contribution of the bulk and edge current distributions to the total magnetization . second , we compare quantitatively and qualitatively the bulk contribution in the finite system to the magnetization produced by one unit cell in the infinite system . we point out that experimentally it may be possible to measure only the contribution to the magnetization caused by the periodic modulation by placing the entire squid - loop well within the sample . and third , we compare to the total thermodynamic magnetization of an infinite system which is periodic in only one spatial direction and thus not subject to commensurability difficulties . our calculations show that the bulk contributions to the magnetization are strongly dependent on the presence or absence of the exchange interaction in the models supporting the view of meinel et . that magnetization is an ideal probe of the many - body effects in a 2deg.@xcite magnetization has been calculated for quantum hall systems in the fractional regime with higher order approximations that reproduce more reliably exchange and correlation effects.@xcite here we focus our attention on relatively large structured systems with several landau bands included and have thus to resort to the hfa in order to make the calculation tractable in cpu - time . recent comparison between results from exact numerical diagonalization and the hfa show that in high landau levels the two approaches agree on the formation of charge density waves.@xcite    we shall consider periodic potentials of a short period , 50 nm , which means short with respect to the present technical possibilities , but still realistic . in this case , for gaas parameters and for magnetic fields in the range of few tesla , the screening effects due to the direct columb interaction are weak . however , the exchange effects remain strong , and they amplify the single - particle energy dispersion.@xcite therefore the presence of the periodic potential should become proeminent in the magnetization even for weak amplitudes . and last , but not least , by avoiding strong screening we are benefitted by a shorter computational time . the magnetization is calculated within three models in the paper , self - consistently with respect to the electron - electron coulomb interaction : a finite model using the ha , an infinite model periodic in both spatial directions , using the unrestricted hartree fock approximation ( uhfa ) , and an infinite model periodic in only one spatial direction , using the standard hfa . the model for a finite system consists of a laterally confined 2deg . a hard wall potential ensures that the electrons stay in the square region @xmath5 the wave functions being zero at the boundary . an external modulating potential and a perpendicular magnetic field are applied . the potential has the form @xmath6 where @xmath7 and @xmath8 count the number of dots in @xmath9 and @xmath10 direction respectively , giving in total @xmath11 unit cells . the schrdinger equation is solved by expanding the eigenfunctions in fourier sine - series and the expansion coefficients are found by diagonalizing the hamiltonian matrix . the electron interaction is taken into account in the hartree approximation . the total magnetization can be calculated according to the definition for the orbital and the spin component of the magnetization,@xcite @xmath12 where @xmath13 is the total area of the system . the equilibrium local current is evaluated as the quantum thermal average of the current operator , @xmath14 with velocity operator @xmath15/m^*$ ] , * a * being the vector potential . even though the magnetic field @xmath2 can be varied freely in this model we have used the definition of the orbital magnetization ( [ m_os ] ) rather than evaluating the derivative of the free energy with respect to the magnetic field . the two - dimensional modulation of the system is a square lattice of quantum dots ( @xmath16 ) or antidots ( @xmath17 ) determined by the static external potential @xmath18 or a simple cosine - modulation defined by @xmath19 where @xmath20 is the length of the fundamental inverse lattice vectors , @xmath21 , and @xmath22 . the bravais lattice defined by @xmath23 or @xmath24 has a periodic length @xmath4 and the inverse lattice is spanned by @xmath25 , with @xmath26 . the commensurability condition between the magnetic length @xmath27 and the period @xmath4 requires magnetic - field values of the form @xmath28 , with @xmath29 , and @xmath30 the magnetic flux quantum.@xcite arbitrary rational values can , in principle , be obtained by resizing the unit cell in the bravais lattice .    for this model we evaluate the contribution of the periodic modulation to @xmath31 and @xmath32 , rather than the total magnetization . using in eq.([m_os ] ) the periodicity of the current and spin densities , and the reflection symmetry of the unit cell , we reduce the integrations to a single cell . obviously , in the absence of the modulation @xmath33 and the orbital contribution vanishes . the ground - state properties of the interacting 2deg in a perpendicular homogeneous magnetic field @xmath34 and the periodic potential are calculated within the uhfa for the coulomb interacting electrons at a finite temperature.@xcite the approximation is unrestricted in the sense that the single - electron states do not have to be eigenstates of @xmath35 . the modulation is defined by the potential @xmath36 describing an array of parallel quantum wires . in this case there is no restriction on the magnetic - field values , the magnetic flux through one lattice cell always being infinite . the groundstate is calculated in the hfa , by diagonalizing the hamiltonian in the landau basis , and by expanding the matrix elements in fourier series . therefore we can evaluate directly the total magnetization @xmath37 by the thermodynamic formula appropriate for the canonical ensemble , @xmath38 where @xmath39 is the total energy , and @xmath40 the entropy . we shall assume the temperature is sufficiently low to neglect the second term of eq . ( [ thermodeq ] ) . in view of more realistic results we also assume a small disorder broadening of the landau levels , which we take into account with a gaussian model for the spectral function.@xcite the numerical calculations are performed with gaas parameters , @xmath41 , and @xmath42 . in the case of the infinite periodic modulation , in the uhfa , hfa , or the ha , the bare @xmath20-factor is -0.44 , and is set equal to zero in the model of the finite 2deg in the ha . mostly for numerical reasons we keep a finite temperature , which for the models with 2d potential is 1 k , and for the 1d potential is 0.2 k. in all cases the length of the unit cell is @xmath43 nm . the magnetization for the finite system is shown in fig.[mag_tot_sie ] . the system size is progressively increased , starting from a single cell of 50@xmath4450 nm@xmath45 , to a system of 5@xmath445 cells , keeping the unit - cell size constant . each cell is defined by one period of the modulation potential ( [ v_sq ] ) . when the system consists of more than one cell we , ad hoc , divide the magnetization into an edge part @xmath46 with a contribution only from the first row of cells around the system , and a bulk part @xmath47 with the contribution from the rest of the cells . below we shall see that generally , the magnetization @xmath47 does approach the orbital magnetization expected for a large system as the number of cells @xmath48 is increased . the variable on the @xmath9-axis of the figure , @xmath49 , i.  e.  the number of electrons in a single cell , can approximately be interpreted as a filling factor for the 3@xmath443 and the 5@xmath445 system . this is confirmed by the evolution of the chemical potential @xmath50 through the single electron hartree - states depicted in fig . [ ehns_tot_sie ] . for even - integer values of @xmath49 @xmath50 jumps through gaps of sparsely distributed edge states separating states concentrated into precursors of landau bands . the bulk @xmath47 and the edge @xmath46 contributions to the magnetization as well as the total magnetization are of similar magnitudes . however the oscillations of @xmath47 around zero are more symmetric than those of @xmath46 . this is because the direction of the edge current is more commonly as expected from the classical clockwise motion of the electrons around the sample , thus giving a preferred sign to @xmath46 . modestly increasing the size , from 3@xmath443 to 5@xmath445 unit cells , clearly gives the finite system more of the character of an extended system . the bulk magnetization for the large system is small when @xmath51 is not close to even integers due to the strong screening of the modulation potential away from the edges of the system . this is confirmed by fig . [ mag_tot_0_sie ] showing @xmath31 of the noninteracting 5@xmath445 system . the structures around even integer values of @xmath51 are less steep than for the interacting system . in the presence of the interaction the energy gaps are reduced by screening , which is self - consistently determined by the density of states around @xmath50 . when @xmath50 lies within one band ( i.  e.  @xmath51 is not an even integer ) the coulomb repulsion forces the electron density to spread out more evenly shifting the  effective filling  @xmath49 a small amount . a slight change in the magnetic flux in the finite system only shifts the relation between the number of electrons and the effective filling factor.@xcite    the current distribution for a 5@xmath445 system is shown in fig . [ curr_sie ] which reveals a strong edge current , but also a bulk current structured self - consistently in a complex way by the modulation , the interaction , and the location of @xmath50 with respect to the energy levels . this interplay of complex bulk contributions with the effects of the edge currents opens the question what are the effects of a modulation in an extended electron system on the magnetization . next , we turn our attention to the infinite system , that is modulated in two directions , and calculate the contribution to the magnetization from one unit cell . in fig . [ e_tot_cos_pq2 ] the total energy is shown as a function of the filling factor @xmath51 for the extended periodic 2deg in the uhfa and the ha for the case @xmath52 . two magnetic flux quanta flow through the unit cell and each landau band is split into two subbands which in turn are doubly spin split . filling factor two means thus that both spin states of one landau band are occupied , and in total four subbands are below the fermi level . the modulation with @xmath53 or 0.1 mev is small compared to @xmath54 mev . the minima in the total energy for the uhfa reflect the strong exchange interaction for electrons , added to nearly filled landau bands or subbands thereof . [ m_cos_pq2 ] compares the total magnetization , eq . ( [ m_os ] ) , and its components @xmath31 and @xmath32 calculated within the uhfa , with the total magnetization according to the ha . the main difference between the results of these two approximations is the sharp reduction in the magnetization caused by the exchange interaction around odd integer filling factors . in this region the enhanced spin splitting of the subbands is larger than the subband splitting caused by the modulation . the order of the subbands ( with respect to spin and magnetic subband index ) and their curvature thus leads to @xmath31 being of same sign for @xmath55 and @xmath56 . the behavior is thus different around the even and odd values of @xmath51 . later we see that this is not in the case of @xmath57 . just like for the total energy the different modulation strengths result in minor changes in @xmath31 , because the energy dispersion of the landau bands is determined by the exchange energy rather than by the external or screened potentials . @xcite    the light effective mass and the small @xmath20-factor of electrons in gaas cause the spin contribution @xmath32 to be an order of magnitude smaller than the orbital one . a comparison of @xmath32 in these two approximations can be seen in fig . [ s_mu_pq2 ] . in the case of the uhfa the exchange interaction always leads to the maximum spin polarization . in the ha the spin splitting of the landau bands is only the bare zeeman gap , here about 0.07 mev , i.  e.   much smaller than the intra - band energy dispersion which is of the order of the modulation amplitude , @xmath58 mev . therefore the chemical potential is never able to lie only in one spin subband . new electrons are being added to the system concurrently , to both spin states , resulting in a reduced polarization .    in the case of one flux quantum through a unit cell ( @xmath57 ) each landau band consists of only two subbands , with different spin quantum numbers . this simpler band structure is reflected in the magnetization ( see fig.[m_v01_pq1 ] ) . here the spin - bands hosting @xmath50 for the filling range @xmath59 and @xmath60 have opposite curvature causing a maximum and a minimum in @xmath31 , respectively . for the lower @xmath51 range the fermi level ( @xmath50 ) is in the lower spin subband of the second landau - band pair , see fig . [ ebpq1 ] . here the fermi contours (  fermi surface  ) encircle the energy maxima at the edges of the brillouin - zone ( hole orbits ) , while for the higher @xmath51 range they close around the minimum of the upper spin subband ( electron orbits ) .    in the case of two flux quanta each landau band is fourfold split as can be seen in fig . [ ebpq2 ] . for @xmath61 only the splitting due to the modulation is clearly visible but for @xmath51 not equal to an even integer the spin splitting gets more enhanched due to the strong exchange force . furthermore , we see that the subbands repel each other around the @xmath62 point ( center of the brillouin zone ) resulting in a different occurrence as @xmath51 is swept from 2.8 to to 3.3 . in contrast to the strong @xmath51-dependent interaction effects on the energy spectra we show in fig . [ eb0 ] the  static  energy bands of the noninteracting system that are independent of @xmath51 and the location of @xmath50 .    around @xmath63 the order of the landau subbands is unusual . the states below the fermi level are those of the landau subband @xmath64 , @xmath65 being the orbital and @xmath66 the spin quantum numbers . but the landau subbands above , i.  e.  @xmath67 and @xmath68 , are overlapped due to the strong exchange enhancement of the gaps between all the subbands with the same @xmath65 but opposite @xmath66 . therefore the first states which are populated for @xmath69 belong to the subband @xmath68 , while the subband @xmath67 is populated a bit later ( for slightly higher @xmath51 ) . this situation is well known in the atomic physics of complex atoms , and has also been demonstrated in quantum dots , with a current - spin density - functional approach,@xcite which in principle is more reliable than our uhfa . in our case we observe this spin - state inversion as a small anomaly in @xmath32 around @xmath63 , fig . [ m_v01_pq1 ] , where the maximum of @xmath32 is shifted to @xmath69 . the fact that the magnetization calculated by eq . ( [ m_os ] ) for a unit cell in an infinite doubly modulated system reflects _ only _ the contribution of the modulation is seen in fig . [ m_v5_qad_pq2 ] , where @xmath70 for a dot and an antidot modulation of the same strength are mirror symmetric around zero for low @xmath51 . for higher @xmath51 the mixing of the landau bands due to the coulomb interaction slightly skews the mirror symmetry . for a homogeneous system ( @xmath71 ) the persistent current in the definition of the magnetization ( [ m_os ] ) vanishes and thus @xmath31 . similar effect was seen for the bulk magnetization of the finite system in fig . [ mag_tot_sie ] for noninteger values of @xmath51 . due to the enhancement of subband width caused by the exchange force the transition to zero magnetization with decreasing modulation does not happen in a smooth linear fashion.@xcite      now we come back to the question how the magnetization of the finite and infinite system are related . the magnetization for the finite system can be calculated by either equation ( [ m_os ] ) or ( [ thermodeq ] ) giving the same results for the low temperature assumed here . the orbital magnetization @xmath31 is compared for the finite system of different sizes to the result for a single unit cell in an extended system with the same type of modulation in fig . [ m_samanb_1h ] , but with the interaction accounted for in different ways . if we first look at the results for the infinite periodic system we notice that the largest variation of @xmath72 occurs for the uhfa and the smallest for the ha with the noninteracting case inbetween . this is in accordance with the screening properties mentioned earlier , in the ha the modulation is screened more effectively than in the uhfa . this result is in agreement with the simplistically defined bulk magnetization @xmath47 of the finite system seen in fig . [ mag_tot_sie ] , the main difference being the sharp variation of @xmath47 at even integer values of @xmath51 . their presence indicates that even though the current density is only integrated over the  bulk  of the system the underlying energy spectrum is affected by the chemical potential @xmath50 traversing it s edge states . to be able to get the unit cell of the two different systems to give the same magnetization the finite system has to be even larger , exhausting our means of computation . the magnetization of both systems ( and also of the 1d modulated system , see the next subsection ) compares well with results derived from an older model of statistical inhomogeneities in a 2deg that was used to explain effects caused by oscillating landau level width due to the electrostatic screening.@xcite    the main differences between the magnetization of our finite and infinite periodic systems are : _ i _ ) the asymmetry around the zero line in the case of the finite system , and _ ii _ ) the missing steepness of @xmath31 around even integer filling factors @xmath51 for the infinite system . earlier we saw that the asymmetry is influenced by the contribution from the  edge  of the finite system . the second difference can also be traced to the edge states . in the extended model there are no edge states between the landau subbands . their shape and curvature can thus change sharply with the motion of the chemical potential @xmath50 through them , the self - consistent screening and exchange effects minimize any gaps that might evolve around @xmath50 , which in turn prevents any sharp jumps in @xmath31 . with this in mind it is clear that _ the magnetization for a realistic ( large , but finite ) modulated system is not simply the sum of the magnetization produced by two independent subsystems , the edge and the bulk_. the coulomb interaction makes the separation of the contributions to @xmath70 a nontrivial problem , which can be solved only by an experiment . in addition , we have seen that the self - consistent motion of @xmath50 through the energy bands depends strongly on the approximation used for the electron interaction . we have noticed in our calculations that for the system sizes considered here the magnetization of the finite system is not strongly dependent on whether the modulation is assumed 2d as here , or 1d.@xcite for an infinite system with a 1d modulation we can calculate the thermodynamically defined magnetization ( [ thermodeq ] ) presented in fig . [ 1dmod ] . as mentioned before we see here that @xmath31 for the finite system , especially when it is enlarged , bears a strong similarity with @xmath31 for the infinite 1d modulated 2deg . the calculation of the ground state for a modulated 2deg with arbitrary magnetic fields is impossible for the 2d potential , due to the commensurability restrictions . the problem can be circumvented for a 1d modulated system that we shall now turn our attention to , and formulate predictions of experimental results that can be used to test the importance of the exchange interaction . for such a system we have access to the total magnetization , i.  e.  bulk plus edge , in the thermodynamic limit . in fig . [ 1dmod ] we display results for the infinite system with a 1d potential , obtained with eq . ( [ thermodeq ] ) . first we show the sawtooth profile in the absence of a modulation potential , reflecting the instability of the fermi level in the energy gaps . the exchange interaction determines the spin splitting for odd filling factors , but also amplifies the jumps for even filling factors by almost a factor of two . the reason is the enhancement effect on both the spin gaps and the landau gaps . for the same reason , in the presence of a modulation the jumps are only slightly reduced . similarily , as for the 2d modulation , the exchange interaction also increases the energy dispersion of the landau bands for noninteger filling factors , by lowering the energy of the occupied states , see figs . [ ebpq1 ] and [ ebpq2 ] . hence , the band width depends on the position of the fermi level inside an energy band . this fact prevents the coincidence of the fermi level with a band edge ( top or bottom ) , resulting in sharp cusps for @xmath51 close to integers . the sharpness is an effect of the exchange interaction in the vicinity of a van hove singularity . some amount of disorder may indeed broaden such cusps . in addition , the magnetization jumps may now slightly increase , because of the smearing of the band edges by disorder , which helps the fermi level to enter or to leave a landau band . when the fermi level lies in the middle of a band screening effects are important . in principle screening is important when the modulation period is much bigger than the magnetic length and/or when high landau bands , with extended wave functions , are occupied . however , even here we can see some oscillations in the upper bands , with orbital quantum number @xmath73 . increasing the modulation amplitude , from @xmath74 mev to @xmath75 mev we first see that the magnetization for the bands with @xmath76 is relatively stable . just like in figs . [ e_tot_cos_pq2 ] and [ m_cos_pq2 ] , this is because the exchange amplification of the energy dispersion depends on the filling factor , rather than on the modulation amplitude.@xcite the spin splitting survives now only for @xmath77 , and it is abruptly supressed for @xmath78 . a similar supression occurs for @xmath74 mev , but at a higher @xmath51 , and it can be explained by the inflation of the wave functions in high landau levels which rapidly equilibrates the number of spin - up and spin - down electrons and destroys the enhancement of the spin gap . @xcite such a suppression effect has been recently observed in magnetotransport experiments on short - period modulated systems , in the shubnikov - de haas peaks.@xcite we have calculated the magnetization of periodic systems and discussed with examples various properties of it due to system boundaries , periodicity , and coulomb interaction . we have compared the results of the finite and infinite systems and of the periodic systems in one and two spatial directions . our aim is to provide information for understanding the magnetization measurements in mesoscopic systems , which are expected to become a new direction of experimental investigations .    unlike in other types of experiments , like transport or electromagnetic absorption , the magnetization measurements seem to open a better and more direct access to the intrinsic , quantum electronic structure of the system . in transport experiments this is often intermediated by complicated electron - impurity interactions , and in far - infrared absorption usually the classical collective motion of the electron system is dominant . we have identified in the magnetization several properties of the energy spectrum which are absent or incompletely observed in transport or absorption measuremens , like the exchange enhancement of the energy dispersion or the curvature of the landau bands . according to a recent prediction the exchange effects may also determine hysteresis properties when acting on the energy dispersion , either by varying the modulation amplitude , or by varying the zeeman splitting in tilted magnetic fields , and keeping the filling factor constant.@xcite the magnetization measurements could be the best suited tool for probing such effects . the present calculation further indicates that sought after delicate internal structure of the landau bands , such as the hofstadter butterfly,@xcite in a doubly periodic 2deg might be better observed by magnetization than transport experiments.@xcite this research was supported by the icelandic natural science foundation and the university of iceland research fund . in particular we acknowledge additional support from the student research fund ( s.  e. ) , nato science fellowship ( a.  m. ) , and graduiertenkolleg `` physik nanostrukturierter festkrper '' ( v.  g. ) .
Abstract:

Answer: we calculate the magnetization of the two - dimensional electron gas in a short - period lateral superlattice , with the coulomb interaction included in hartree and hartree - fock approximations . we compare the results for a finite , mesoscopic system modulated by a periodic potential , with the results for the infinite periodic system . in addition to the expected strong exchange effects , the size of the system , the type and the strength of the lateral modulation leave their fingerprints on the magnetization .


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: interstellar medium is turbulent with magnetic fields playing a critical role in the key processes , including star formation , propagation and acceleration of cosmic rays as well as regulating heat and mass transfer between different phases . therefore it is essential to have a reliable way of studying magnetic fields . unfortunately , there only a few ways of tracing magnetic fields and their applicability is limited . for instance , tracing magnetic fields with aligned grains ( see @xcite ) requires either the availability of the background stars , if the polarization via extinction is employed , or sensitive far infrared measurements , which frequently imply the use of the space - based instruments . another major technique that employs synchrotron polarization mostly traces magnetic fields in the warm and hot phases of the interstellar medium ( ism ) ( see @xcite for the list of the idealised ism phases ) .    using spectroscopic data for tracing magnetic fields was first suggested and demonstrated in @xcite using the anisotropy of correlation functions of the velocity channel maps . the new approach was motivated by the theoretical advancements establishing that the mhd turbulence is anisotropic ( @xcite , herafter gs95 , @xcite , see @xcite for a review ) and should translate into the anisotropy of the observed correlation fluctuations . later this approach was employed with the use of velocity centroids @xcite and @xcite and the principal component analysis ( pca ) @xcite of the spectroscopic data . in its pca incarnation the idea of tracing magnetic field anisotropies was applied to the observational data in @xcite with the results of anisotropies shown to be consistent with the starlight polarization measurements . further development of the technique resulted in the development of techniques using the anisotropies for obtaining the media magnetization , as it parameterized by the alfven mach number @xcite and the separation of compressible and incompressible turbulence components of observational data as discussed in @xcite . on the basis of the aforementioned improved understanding of the nature of mhd cascade a new technique employing velocity centroid gradients ( vcgs ) was introduced in @xcite ( @xcite , henceforth gl17 ) . instead of appealing to anisotropic velocity correlations as in @xcite the new technique made use of velocity gradients that , according to the mhd turbulence theory tend to be perpendicular to the _ local _ direction of the magnetic field . this idea was further extended in @xcite ( @xcite , henceforth yl17 ) where a robust practical procedure for vcg calculations was suggested . in the latter paper we also applied the new technique to the galfa 21 cm data and demonstrated a good correspondence between the magnetic field tracing with the vcgs and the planck polarization arising from aligned grains . this result motivates the further study of vcgs properties that we undertake in this paper . we note parenthetically that vcgs are not the only way to trace magnetic field using velocity gradients . intensity gradients within the position - position- velocity ( ppv ) channel maps have been recently suggested as an alternative way of exploring the magnetic field topology @xcite . on the basis of the earlier studies of the statistics of channel maps versus the statistics obtained with velocity centroids ( see @xcite and ref . therein ) one can expect that the velocity channel gradients in @xcite ( @xcite , henceforth ly17 ) may have advantages for studying supersonic turbulence as well as tracing the 3d distribution of turbulence . at the same time , the vcgs may have advantages for studying subsonic turbulence . the comparison of the two technique will be done elsewhere . another way of exploring the gradients was introduced in @xcite ( @xcite , henceforth lylc17 ) . this approach employs synchrotron intensity gradients ( sigs ) and provides a way of tracing magnetic field in turbulent environments the synergy of which to the vcg technique we explore in this paper . the theoretical basis for the sigs is the same as for the vcgs however , in the presence of gravity , the response of the sigs and the vcgs is different which provides a way to identify the regions where the matter is subject to the gravitational collapse , as we discuss in the paper . the density statistics is also affected by mhd turbulence . mhd turbulence can imprint some of its properties to density , making its statistics at low sonic mach numbers anisotropic ( see @xcite ) . while both from theory and simulations we expect the density to be a worse tracer of magnetic fields , especially at high sonic mach numbers , the very differences between the velocity and anisotropies can be informative of the properties of the turbulent media . this is the synergy that we also explore in this paper . we note that the first report of the correlation of the magnetic fields and density gradients was done within the numerical study by @xcite reported the alignment of density gradients and magnetic fields with simulation as well as the change of the alignment that occurs after the effects of self - gravity get important starting with a particular density threshold . in our study we make a step further and explore the point - wise tracing of magnetic fields by the intensity gradients ( igs ) as well as the nature of the deviations of the igs directions from those of magnetic field . we show a significant synergy of using the igs , the vcgs and the dust polarimetry data . in addition , we provide a detailed study on how self - gravity lead to the rotation of both vcg and ig vectors and the physics behind it . we also explore the ways of improving the ig tracing of magnetic field by suggesting the procedure filtering for the shocked gas regions .    as a separate but related development it is important to mention the study in @xcite . there a good correlation of directions of the filaments observed in the ppv galfa channel maps and the magnetic field direction as revealed by planck polarimetry was reported . as a result the filaments were identified as a tracer of magnetic field direction in hi high latitude gas . while in the original study the filaments were identified with the actual spacial density enhancements , we believe that the interpretation depends on the thickness of the velocity channel used . indeed , according to @xcite within the sufficiently thin velocity channels the intensity fluctuations are caustics mostly induced by velocity fluctuations . therefore , depending on the velocity channel thickness , the filament tracing is related either to revealing of velocity statistics or density statistics or both . the comparison of the empirical techniques is an interesting problem that we will discuss in another paper . self - gravity is expected to modify the properties of the flow in the vicinity centers of the gravitational collapse . the change of the properties of mhd turbulence within such regions is expected to be the result of vcgs changing their properties , in particular their alignment with the magnetic fields . the density distribution is even more affected by self - gravity . indeed , a recent study by @xcite , reported that the igs rotate with respect to magnetic fields when arriving some density threshold in super - alfvenic high density self - gravitating simulation . while the work from @xcite does not address how self - gravity dynamically change the alignment of igs and vcgs , that encourages a detailed study on how self - gravity can possibly lead to the rotation of gradient vectors and the physics behind it . in this paper , we address the issue of how vcgs and igs behave in the regions of self - dominant gravity for the case of sub- and trans - alfvenic clouds , how does the stage of collapse correlated to the alignment of gradients , and how can the aforementioned phenomenon can be used on tracking self - gravitating regions in observation .    in what follows we briefly summarize the theoretical arguments that justify the use of the vcgs in  [ sec:1 ] , explain our simulation setting in  [ sec:1.5 ] and numerical approach for gradient calculation in  [ sec:2 ] . we provide a comparison of the alignment of vcgs and igs with magnetic fields in different physical conditions in  [ sec:3 ] , and providing methods on improving the alignment . we explore the properties of the vcgs in sub - alfvenic self - gravitating clouds in  [ sec:4 ] and explore their synergistic use with the igs . we illustrate our method in  [ sec:5 ] using observational data in  [ sec:5 ] . we compare our methods to the method of anisotropy of correlation in  [ sec:6 ] . we compare our method to others and discuss our results in  [ sec:7 ] . we summarize our results in  [ sec:8 ] . the theoretical motivation for studying magnetic fields using vcgs are discussed in our earlier papers ( gl17 , yl17 , ly17 ) as well as in the paper on the sigs ( lylc17 ) . in short , in the presence of fast turbulent reconnection @xcite the motions that do not induce magnetic field bending are preferentially excited . these are the motions perpendicular to the local direction of magnetic field . this fact is numerically confirmed in @xcite and is the corner stone of the modern theory of mhd turbulence ( see @xcite ) . due to this alignment of velocity motions , the gradient of velocities are largest perpendicular to the _ local _ direction of magnetic field . it is important that the gradient amplitude for the gs95 turbulence is increasing with the decrease of the eddy scale . therefore the smallest eddies that are most aligned with the local direction of the magnetic field contribute most to the gradients . as a result , by tracing the velocity gradients one can trace the magnetic field direction . the motions at the turbulence injection scale can be super - alfvenic , i.e. faster than the alfven velocity @xmath0 . the gs95 theory is not applicable to such motions . however , as the turbulent velocity decrease with the scale , e.g. for kolmogorov turbulence @xmath1 , for a realistically extended turbulence inertial range at small scales the motions gets subalfvenic and obey the gs95 theory . therefore , by filtering the contributions from larger scales we may still trace the magnetic fields using velocity gradients .    in the presence of gravity , the velocity field is modified . in particular , for subcritical collapse when magnetic field dominate the dynamics , the motions are expected along the stiff magnetic field lines , making velocity gradients aligned with magnetic field . similarly , for supercritical collapse , the collapsing material is expected to drag relatively weak magnetic fields . this should also make magnetic fields and velocity gradients aligned . in other words , for the regions with significant self - gravity effects we expect the change of the direction of magnetic field and velocity gradient alignment from being perpendicular to magnetic field to becoming parallel to magnetic field . ( see section [ subsec:4.1 ] for a detailed discussion )    in gl17 and yl17 the velocity centroids were used as a measure of velocity that is available from observations . the velocity centroids were shown to reflect the statistics of underlying velocity field well , at least for subsonic mhd turbulence @xcite . in the present paper we explore numerically the utility for the vcgs to trace magnetic fields in mhd turbulence for various turbulent driving environments , e.g. different combination of sonic and alfvenic mach number , different compression of fluids in non self - gravitating as well as self - gravitating environment . we shall also establish the guidelines on how to get better alignments numerically even the gradient technique does not provide a good prediction . as there are many factors affecting the alignment of rotated gradients with magnetic fields , we would investigate vcgs and igs through studying mhd simulations in different physical conditions . following yl17 and lylc17 we used a series of compressible , turbulent , isothermal single fluid magnetohydrodynamic ( mhd ) simulations for examination of vcgs and igs from zeus - mp/3d @xcite , a variant of the code zeus - mp @xcite . we set up simulation cubes with several combinations of sonic mach number @xmath2 and alfvenic mach number @xmath3 , which are listed in table [ tab : simualtionparameters ] . for supersonic simulations , we increases the resolution of the simulations to ensure the cube has a sufficiently long inertial range for high @xmath2 cases . the initial magnetic field is along the z - axis and the boundary condition is triply periodic . we injected turbulence using the same method in yl17 . to study the effect of compressibility to alignment of gradients , we also introduce a cube ( set e ) from an incompressible simulation ( see lylc17 for details ) . most of this work is devoted to the studies of diffuse media with the negligible effects of gravity . in section[sec:4 ] we study the case of strongly magnetized media for which self - gravity is important . there we wait until turbulence inside the cube gets saturated and then switch on the self - gravitating module , which employs a periodic fast fourier transform ( fft ) poisson solver . to trace the effect of gravity with respect to time , we take a number of snapshots with even time intervals before and after the of gravity is switch on . we stop our calculation when any pixel in the system violates the truelove criterion @xcite , which informs us when to stop the simulation before having numerical artifacts due to self - gravitating collapse . as self - gravity is a dynamic process , we expect systems having same @xmath2 and @xmath3 but different magnetic criticality , that means a system with gravitational to magnetic energy ratio larger than one , will provide very similar final results . therefore we here only investigate the systems with magnetic criticality @xmath4 to be 2 , which ensures the cube will have physical collapse when self - gravity is activated and before the truelove criterion is violated . the raw data from simulation cube has to be converted to synthetic maps for our gradient studies . assuming the line - of - sight direction is the @xmath5-axis , the intensity @xmath6 , velocity centroid @xmath7 are defined as @xmath8    we also computed polarization by assuming the media is optically thin , faraday rotation is not effective along the line of sight , and the emissitivity is a constant over the line of sight . in other words , the stokes parameters @xmath9 can be expressed in terms of @xmath10 , the y and z direction magnetic fields : @xmath11 the polarized intensity @xmath12 and angle @xmath13 are then defined correspondingly . ( see @xcite for details ) . the three observables @xmath14 will be used to calculate the intensity gradients ( igs ) , velocity centroid gradients ( vcgs ) and synchrotron intensity gradients ( sigs ) , by following the recipe and the _ sub - block averaging_. the _ sub - block averaging _ method , as it introduced in yl17 is based on defining the most probable direction within a block by a local gaussian - fitting peak on the distribution of gradient angle distribution . within the approach , the fitting error provides a quantitative estimate of whether the block size is large enough for providing a correct direction . to compare the performance between vcgs and igs , we employed the _ alignment measure _ that is introduced in analogy with the grain alignment studies ( see @xcite ) : @xmath15 ( see gl17 , yl17 ) with a range of @xmath16 $ ] to measure the relative alignment between _ rotated _ gradients and magnetic fields , where @xmath17 is the relative angle between _ rotated _ gradients to magnetic field . a perfect alignment gives @xmath18 . on the other hand , @xmath19 when vectors tend to be more perpendicular . we shall use am to quantify the performance of the vcgs and igs in various physical conditions . some special values for am can also give the readers how good the alignment it is for certain type of gradient vectors to magnetic fields . for example , @xmath20 indicates the dispersion between _ average _ gradients to magnetic field directions is about @xmath21 , which means the gradient vectors are rather good tracers for the magnetic field direction . on the other hand , @xmath22 indicates a random orientation of gradient vectors relative to magnetic fields . in what follows , we shall study both vcgs and igs and investigate their advantages and synergy for studying magnetic field and other parameters of the ism . shocks are a very important process when the sonic mach number , @xmath2 , which is the ratio of the turbulent injection velocity and the speed of sound , gets large . different phases of the ism are expected to have different sonic mach numbers @xmath2 ( see @xcite ) . the mach number @xmath2 of warm gas is of the order of unity ( see @xcite and ref . therein ) , while the mach number of cold interstellar media ( see @xcite ) as well as molecular clouds can be @xmath23 and even higher . the sonic mach number also characterize the compressibility of turbulent flow . when @xmath2 increases , stronger compressions result in stronger shocks and significantly modify the properties of turbulent gas . we study how shocks change the alignment of gradient vectors and the underlying magnetic field . we start our investigation with the cubes from set `` b '' ( see table [ tab : simulationparameters ] ) to perform an alignment test on both the vcgs and the igs . figure [ fig : vcgig - ms ] shows how the vcgs and the igs behave with various @xmath2 . we note that the alignment of the vcgs in significantly better compared with the igs .     a diagram showing the alignment of vcgs ( red ) , igs ( green ) , and the relative alignment between them ( blue ) as the sonic mach number @xmath2 increases , scaledwidth=44.0% ]    [ tab : simualtionparameters ]    .simulations used in our current work . the magnetic criticality @xmath24 is set to be 2 for all simulation data . [ cols="^,^,^,^,^",options="header " , ]     in a subsequent study on observational studies on the same vela c region ( yuen et al . in prep ) , we not only be able to trace the direction of magnetic field in appropriate scales with very high @xmath25 , but also locating the possible self - gravitating region using the numerical result in this work . that suggests the applicability of our technique for magnetic field tracing even within extreme environments .      in this work , we investigate a number of different regimes for the sonic mach number @xmath2 , but the explored range of the alfven mach numbers @xmath3 is limited to sub- and trans - alfvenic cases only . while our series of studies ( yl17 , lylc17 , this work , @xcite ) shows the power of our gradient technique in different circumstances , one very important regime that is yet to explore is the highly super - alfvenic cases . on the other hand , the expectation from turbulence theory ( gs95 , lv99 ) also make the tracing through gradients more complicated than sub- and trans - alfvenic cases . only sufficient small scales @xmath26 will the turbulent eddies become anisotropic . in fact , the wavenumber removal suggested in section [ subsec:3.2 ] should be able to reveal the alignment of the vcgs to underlying magnetic field provided that @xmath27 is known . the practical procedures for the actual data handling will be discussed elsewhere . our paper explores the use of the velocity centroid gradients ( vcgs ) for tracing magnetic fields and its synergy with other gradient techniques and polarimetry for probing the ism processes . we tested how vcgs , intensity gradients ( igs ) and synchrotron intensity gradients ( sigs ) behave under different physical conditions , including a change of sonic and alfvenic mach numbers and in the presence of self - gravity . we show that in the diffuse media the vcgs trace the magnetic field well and showed that in combination with polarimetry they can also be efficient in identifying the regions dominated by self - gravity . our results in the paper can be briefly summarized :    1 .   in the diffuse media the vcgs trace the magnetic field well and in combination with polarimetry or sigs they can also be efficient in identifying the regions dominated by self - gravity . 2 .   the deviations of the direction of the vcgs and the igs can signify the existence of shocks and can be used as a ways to obtain the turbulence sonic mach number . 3 .   new ways to improve the performance of the vcgs and igs tracing magnetic fields were identified . they include 1 .   removing pixels with higher ig amplitudes . 2 .   removing the contribution of lower wave number . 4 .   for self - gravitating regions the alignment of the vcgs changes : 1 .   in the case of dynamically important magnetic field , the unconstrained falling of matter along magnetic fields lines makes the gradients aligned with magnetic fields . therefore the change of the relative orientation of polarization directions and those determined by the vcgs can provide the observational signature of the self - gravitating collapse within molecular clouds . the time scale for the response of the igs to the gravity is faster than the vcgs in a collapsing cloud , which shed the possibility of determining the stage of collapse by tracing the relative alignment of the igs to the vcgs in the vicinity of high density regions . this opens a way to identify the regions of gravitational collapse by observing the increase of misalignment of the igs and the vcgs for intermediate densities with the return to their alignment for higher densities . 3 .   combining vcgs and sigs explored earlier in @xcite , it is possible to locate self - gravitating regions by examining the relative orientation of the two gradients . the gravity , rather than the density threshold determine the change of the direction of gradient vectors in self - gravitating regions . the vcgs can provide more detailed magnetic field tracing compared to the correlation function anisotropy technique as suggested in lazarian et al . we illustrate our technique by analyzing the observation data : 1 . for the nearby region of galfa - hi used in @xcite to the combination of the vcgs and the igs reveals shocks . 2 .   the full - sky hi4pi hi data together with planck synchrotron foreground intensity data and planck dust polarization data locate the possible self - gravitating regions . the stay of khy at uw - madison is supported by the fulbright - lee hysan research fellowship and department of physics , cuhk . al acknowledges the support of the nsf grant ast 1212096 , nasa grant nnx14aj53 g as well as a distinguished visitor pve / capes appointment at the physics graduate program of the federal university of rio grande do norte , the inct inespao and physics graduate program / ufrn . gaensler , b.  m. , haverkorn , m. , burkhart , b. , et  al . 2011 , http://dx.doi.org/10.1038/nature10446 [ , 478 , 214 ] galtier , s. , nazarenko , s.  v. , newell , a.  c. , & pouquet , a.  2000 , journal of plasma physics , 63 , 447 goldreich , p. , & sridhar , s. 1995 , http://dx.doi.org/10.1086/174600 [ , 438 , 763 ]                          , a. , pogosyan , d. , & esquivel , a. 2002 , in astronomical society of the pacific conference series , vol . 276 , seeing through the dust : the detection of hi and the exploration of the ism in galaxies , ed . a.  r. taylor , t.  l. landecker , & a.  g. willis , 182 lazarian , a. , & pogosyan , d.  2000 , , 537 , 720
Abstract:

Answer: this study proceeds with the development of the technique employing velocity gradients that were identified in ( @xcite , henceforth gl17 ) as a means of probing magnetic field in interstellar media . we demonstrate a number of practical ways on improving the accuracy of tracing magnetic fields in diffuse interstellar media using velocity centroid gradients ( vcgs ) . addressing the magnetic field tracing in super - alfvenic turbulence we introduce the procedure of filtering low spatial frequencies , that enables magnetic field tracing in the situations when the kinetic energy of turbulent plasmas dominate its magnetic energy . we propose the synergic way of of using vcgs together with intensity gradients ( igs ) , synchrotron intensity gradients ( sigs ) as well as dust polarimetry . we show that while the igs trace magnetic field worse than the vcgs , the deviations of the angle between the igs and vcgs trace the shocks in diffuse media . similarly the perpendicular orientation of the vcgs and the sigs or to the dust polarimetry data trace the regions of gravitational collapse . we demonstrate the utility of combining the vcgs , igs and polarimetry using galfa hi and planck polarimetry data . we also provide an example of synergy of the vcgs and the sigs using the hi4pi full - sky hi survey together with the planck synchrotron intensity data .


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 early hi 21 cm - line absorption and emission studies of the galaxy lead to the two - phase model of the interstellar medium in which dense cool hi clouds ( the cold neutral medium ) were in pressure equilibrium with a warm low - density medium ( the warm neutral medium ( wnm ) ) ( clark , radhakrishnan , & wilson 1962 , clark 1965 , radhakrishnan et al . the mean spin temperature of the cold clouds was estimated to be 80 k. only a lower limit to the spin temperature of the wnm of 1000 k was obtained from an upper limit to the optical depth of hi 21 cm - line absorption in the wnm . since then there have been several attempts to measure the optical depth of the wnm ( mebold & hills 1975 , kalberla , mebold , & reich 1980 , liszt , dickey , & greisen 1982 , mebold et al . 1982 , payne , salpeter , & terzian 1983 , kulkarni et al . 1985 ) . using the bonn telescope , mebold & hills ( 1975 ) estimated an optical depth for the wnm based on the lack of hi emission in observations toward cygnus a. they estimated the spin temperature of the wnm to be in the range of 3000 to 8000 k in different velocity ranges toward cygnus a. no absorption in the wnm was detected in any of the other observations mentioned above leading to upper limits on its optical depth and hence lower limits to the spin temperature . the highest of these lower limits to the spin temperature is 10@xmath6 k ( kulkarni et al . 1985 ) . of all the phases of the interstellar medium , the wnm remains the least understood . although hi 21 cm - line emission from the wnm is detected easily in most directions , an accurate value for the temperature of the wnm is lacking . the spin temperature of the wnm is an important parameter with implications for models of the interstellar medium . the most detailed analysis of the equilibrium state of neutral hydrogen in the interstellar medium is the work of wolfire et al . the primary heating mechanism is photo - electric emission from dust grains , while cooling is dominated by fine structure lines from heavy elements in colder regions , and by hydrogen recombination lines in warmer regions . recently , we detected hi 21 cm - line absorption in the wnm using the westerbork synthesis radio telescope ( wsrt ) . the absorption was detected toward cygnus a at lsr velocities of 40 km s@xmath1 and 70 km s@xmath1 ( figs . 3 and 4 , and table 1 in carilli , dwarakanath , & goss 1998 ) . these two velocity ranges were previously identified as being relatively free of cold absorbing clouds . the measured optical depth for the wnm along the line of sight to cygnus a is ( 8.9 @xmath0 1.9)@xmath2 10@xmath7 at 70 km s@xmath1 and ( 8.5 @xmath0 2.0)@xmath2 10@xmath7 at 40 km s@xmath1 , with corresponding spin temperatures of 6000 @xmath0 1700 and 4800 @xmath0 1600 k , respectively . the volume filling factor for the wnm appears to be fairly high and is @xmath8 0.4 ( carilli , dwarakanath , & goss 1998 ) . the measured value of the spin temperature toward cygnus a is consistent with the range of kinetic temperature ( 5500  8700 k ) for the wnm in the model proposed by wolfire et al ( 1995 ) . while the current data toward cygnus a are consistent with this picture , values of the spin temperature of the wnm in more directions are needed to provide a constrained physical model of the interstellar medium . mebold et al . ( 1981 ) have carried out a survey in galactic hi absorption and emission toward a large number of extragalactic radio sources . from this list we have chosen four sources which are strong enough ( @xmath9 10 jy at 20 cm ) to detect hi absorption ( optical depth @xmath8 0.001 ) in an integration time of @xmath8 12 hours . the hi emission profile toward each of these sources shows a few narrow lines plus a broad line component . the existing hi absorption spectrum toward each of these sources , however , has only narrow components corresponding to the narrow emission lines(mebold et al . these absorption components ( optical depth @xmath9 0.05 ) are due to the cold gas ( t@xmath10 @xmath8 80 k ) in the galaxy . the motivation of the present observations is to detect the hi absorption corresponding to the broad line in the hi emission profile arising from the intercloud medium ( warm neutral medium ) . the four sources chosen for the present observations are at b @xmath9 10@xmath11 . these higher latitudes are relatively free of the cold neutral medium , and because of this , hi absorption measurements toward high - latitude sources are less confused by the plethora of narrow absorption lines routinely detected in the galactic plane . the observations are discussed in section 2 . the data analyses and the results obtained are given in section 3 . the implications for the determinations of the temperature and the density of the warm neutral medium based on the current observation are discussed in section 4.1 . a by - product of the current observations is the detection of hi 21 cm - line absorption from high velocity clouds in the outer arm of the galaxy . a discussion on this can be found in section 4.2 . the hi 21 cm - line absorption measurements were carried out toward 3c147 , 3c273 , 3c295 and 3c380 using the westerbork synthesis radio telescope ( wsrt ) in 1998 and 2000 . some details of the observations are given in table 1 . the observations used a bandwidth of 1.25 mhz centered at the frequency of the hi 21 cm - line and 256 spectral channels . the resulting velocity resolution is @xmath8 2.1 km s@xmath1 after hanning smoothing . absolute flux density calibration was obtained by observing 3c286 and 3c48 . since the sources of interest were bright and unresolved at the wsrt resolutions ( @xmath8 15@xmath12 ) , amplitude and phase calibrations of the antennas were carried out through self - calibration . the bandpass calibration was obtained through frequency switching on the source itself . for this purpose , each ` on - line ' scan of 25 minutes duration was sandwiched between two ` off - line ' scans each of 25 minutes duration . these ` off - line ' scans were offset in frequency from the ` on - line ' scan by @xmath0 1.25 mhz respectively . the ` off - line ' scans are unaffected by contamination due to galactic hi and provide an accurate bandpass . the data were analyzed using the astronomical image processing system ( aips ) developed by the national radio astronomy observatory . for each source , a continuum data base was constructed by averaging the line - free channels . the amplitude and phase solutions of antennas obtained through normal and self - calibration of this data set were transferred to the spectral line data . the calibrated and bandpass corrected spectral line data of the source were used in further analysis . the continuum flux density in the source visibilities was removed by making a linear fit to the visibilities in the line - free channels as a function of channel and subtracting it from the source visibilities in all the spectral channels . these continuum - subtracted visibilities were used to make spectral cubes of hi absorption images . the spectral cubes were cleaned . the rms per channel was typically @xmath8 5 mjy / beam . in order to avoid any possible contamination due to galactic hi emission , visibilities in the short spacings ( the inner @xmath8 0.5 km ) were not used in making the hi absorption images . since the aim of the current observations is to detect weak ( optical depth @xmath810@xmath3 ) , and wide ( width @xmath8 50 km s@xmath1 ) lines , the most important requirement is the stability of the bandpass . this stability can be estimated by comparing the two ` off - line ' bandpasses from the same source in a given observing run . such a comparison for 3c147 , 3c273 and 3c295 indicated that the two ` off - line ' bandpasses on a given source are consistent with each other to within 2@xmath210@xmath7 ( 1@xmath5 ) . the calibrated on - line data of these sources is not expected to have any spurious features beyond this level . this stability is greater than , or comparable to the rms noise in the optical depth images toward these three sources ( table 1 ) . in the case of 3c380 , the two ` off - line ' bandpasses differed at the level of 6@xmath210@xmath7 ( 1@xmath5 ) . in addition , in the case of 3c380 , many antennas showed variations across the passband with periods of @xmath8 0.5 mhz ( @xmath8 100 km s@xmath1 ) . due to the limited dynamic range achieved , the 3c380 data were not considered for wide - line analysis . however , these data were suitable for an investigation of narrow hi absorption lines . the hi absorption spectra toward the four sources in table 1 are displayed in the lower panels in figs . 14 respectively . the corresponding hi emission profiles from the leiden dwingeloo sky survey ( hartmann , & burton 1997 ) are in the upper panels of these figures . a gaussian decomposition of these profiles was carried out using the groningen image processing system ( gipsy ) . each optical depth spectrum was hanning smoothed and simultaneously fit with the minimum number of gaussians until the reduced chi - squared value was @xmath13 1 . in each case , the model spectrum was then subtracted from the observed spectrum to obtain the residual spectrum . the best - fit model spectrum shows a featureless residual spectrum . the residual spectra so obtained toward the sources 3c147 , 3c273 , and 3c295 are shown in fig . similar analyses were also carried out on the corresponding hi emission spectra . the parameters of the best - fit gaussians from these analyses are summarized in table 2 . in many cases , the central velocities of hi absorption and hi emission components agree within errors . in such cases , we assume that the absorption and emission arise from the same hi . an estimate of the column density , the spin temperature , and the width expected due to thermal broadening is also given for these components . if the central velocities do not agree within the quoted errors , either the absorption or the emission column is left blank , indicating that they arise from different hi .    in figs . 5 and 6 , the hi absorption spectra toward 3c147 and 3c380 are shown along with the corresponding hi emission spectra respectively . toward the direction of 3c147 note the weak hi absorption ( fig . 5b ) over the velocity range 70 km s@xmath14 v@xmath15 30 km s@xmath1 and a corresponding hi emission ( fig . the weak hi absorption detected toward 3c147 over the velocity range 70 km s@xmath14 v@xmath15 30 km s@xmath1 ( fig . 5b ) has a 10@xmath5 significance compared to the bandpass stability observed in the 3c147 observations . this hi absorption can be combined with the corresponding hi emission ( fig . 5a ) to estimate the spin temperature of the hi gas responsible for them . the hi emission spectrum in t@xmath16 toward 3c147 ( fig . 5a ) is divided by the optical depth values ( fig . 5b ) to obtain spin temperature as a function of v@xmath17 . over the velocity range 70 km s@xmath14 v@xmath15 30 km s@xmath1 where the weak absorption is detected , the spin temperature has an average value of 3000@xmath0500 k ( 1@xmath5 error ) . the high spin temperature and the wide range of velocities over which hi absorption is detected indicates that this hi gas is due to the warm neutral medium of the galaxy . note that the hi absorption toward 3c380 ( fig . 6b ) also shows a similar @xmath18dip@xmath18 over the velocity range 50 km s@xmath14 v@xmath15 20 km s@xmath1 . this feature has less than a 3@xmath5 significance compared to the bandpass stability in the 3c380 observations .    in fig . 7 , the residual hi absorption spectra obtained after subtracting the best - fit models ( table 2 ) from the observed spectra toward 3c147 , 3c273 and 3c295 are shown in top to bottom panels , respectively . in the case of 3c147 ( fig . 7a ) , the velocity range over which the simultaneous fit was carried out was restricted to 110 km s@xmath14 v@xmath15 + 100 km s@xmath1 in order not to exceed the maximum number of gaussians that gipsy can fit simultaneously . the spectrum shown for 3c147 ( fig . 7a ) is restricted to this velocity range . a separate simultaneous fit was carried out with two gaussians over the velocity range 140 km s@xmath14 v@xmath15 90 km s@xmath1 to account for the two high velocity narrow absorption lines detected toward 3c147 ( fig . 5b and table 2 ) . there is a clear separation between these two velocity ranges of the spectrum with no spectral features in the overlapping region 110 km s@xmath14 v@xmath19 -90 km s@xmath1 ( fig . 5b ) enabling independent gaussian analyses in the two velocity ranges . in fig . 7a two spectra corresponding to 3c147 are shown . the solid line indicates the residual absorption spectrum obtained after subtracting from the observed spectrum the model containing only the seven narrow hi absorption lines in the velocity range 110 km s@xmath14 v@xmath15 + 100 km s@xmath1 ( table 2 ) . the broken line indicates the residual spectrum after including the wide hi absorption line ( v@xmath20 29@xmath04 km s@xmath1 and @xmath21v@xmath22 = 53@xmath06 km s@xmath1 ) in addition to the seven narrow lines in the model spectrum . the reduced chi - squared values of the residual spectra are 16 ( solid line ) and 0.6 ( broken line ) respectively . these two residual spectra and the corresponding reduced chi - squared values clearly indicate the existence of a wide hi absorption line toward 3c147 . in addition , a corresponding wide ( @xmath21v@xmath22 = 34@xmath01.1 km s@xmath1 ) hi emission component at a velocity v@xmath20 35.5@xmath00.3 km s@xmath1 is observed in the ldss ( fig . 5a and table 2 ) . the central velocities of these absorption and emission lines agree to within @xmath81@xmath5 . however , their widths are discrepant at the 3@xmath5 level . we identify these two emission and absorption lines as arising from the same hi gas . the estimated spin temperature of this component is 3600 @xmath0 360 k , in agreement with the direct estimate of spin temperature made in the previous paragraph by dividing the hi emission brightness temperature by the corresponding optical depth as a function of v@xmath17 .    the residual spectra shown in figs . 7b , c corresponding to the hi absorption toward 3c273 and 3c295 do not indicate wide hi absorption , although the hi emission spectra toward these sources do show wide lines ( table 2 ) . in these two cases upper limits to the optical depths of wide absorption lines can be obtained . based on these upper limits and the hi column densities of the corresponding emission lines , lower limits to the spin temperatures of these wide emission lines are estimated ( table 2 ) . for example , toward 3c273 , a lower limit ( 3@xmath5 ) of 2600 k is obtained for the emission line at v@xmath17 =  13.2 km s@xmath1 with a width of @xmath21v@xmath22 = 24.7@xmath01.4 km s@xmath1 .    the hi absorption profile toward 3c147 ( fig . 5b ) also shows two narrow , weak hi absorption lines at v@xmath17 = 117.3 and 124.5 km s@xmath1 due to hi in the outer arm of the galaxy . the hi emission profile toward 3c147 ( fig . 5a ) has a line at v@xmath17 = 120.4 km s@xmath1 with a full width at half maximum of 17 km s@xmath1 ( table 2 ) . it appears that this hi emission and the two hi absorption lines may be associated . however , since there is no exact correspondence of the central velocities and the widths of the emission and absorption lines , we can only estimate limits on the spin temperatures . assuming that the widths of these two absorption lines are due to thermal broadening alone , upper limits to their spin temperatures of 1400 k and 270 k are obtained , respectively . the absence of hi absorption corresponding to the hi emission at v@xmath17 = 120.4 km s@xmath1 leads to a lower limit of 1780 k to its spin temperature . similarly , the hi absorption spectrum toward 3c380 ( fig . 6b ) has a discrete narrow line at v@xmath23113.7 km s@xmath1 with an hi emission line in the ldss spectrum at v@xmath17 = 109.7 km s@xmath1 ( fig . 6a ) . however , hi emission line width is more than ten times the hi absorption line width ( table 2 ) . the lower limit to the spin temperature of this wide emission line with no corresponding hi absorption is 800 k. an upper limit to the spin temperature of 210 k can be estimated for the hi absorption at v@xmath23113.7 km s@xmath1 , if the width is due entirely to thermal broadening . the spin temperature of the hi gas is influenced by the intensity of radiation at the frequency of the 21 cm - line , collisions with electrons , protons , and hydrogen atoms , and the intensity of the ly@xmath24 radiation field ( field 1958 , liszt 2001 ) . in the context of the warm neutral medium , the most important factor influencing the spin temperature of the hi gas is the ly@xmath24 radiation field . the hi spin temperature is expected to be close to the kinetic temperature of the gas through the process of resonant scattering of ambient ly@xmath24 photons ( field 1958 ) . the efficacy of this process depends on the ambient density of ly@xmath24 photons . there are no secure measurements of the diffuse galactic ly@xmath24 photon density ; an upper limit of 2.5@xmath210@xmath25 cm@xmath3 in the solar neighborhood has been estimated ( holberg 1986 ) . the areal production rate of ionizing photons from early - type stars in the solar neighborhood is estimated to be 3.7@xmath210@xmath26 photons @xmath27 s@xmath1 ( vacca , garmany , & shull 1996 ) . assuming this to be comparable to the production rate of ly@xmath28 photons leads to a galactic ly@xmath24 flux from early - type stars . liszt ( 2001 ) discusses the relation between the spin temperature and the kinetic temperature as a function of the number density of h - nuclei , the gas pressure of the interstellar medium , and the fraction of galactic ly@xmath24 flux from early - type stars that permeates the diffuse gas . if the fraction of ly@xmath24 flux from early - type stars that permeates the warm gas is @xmath290.0001 , the spin temperature is significantly lower than the kinetic temperature of the hi gas . the spin temperature of 3600 k estimated in the current observations implies a kinetic temperature of 7000 k and a density of 0.4 cm@xmath3 in the two - phase model proposed by wolfire et al ( 1995 ) ( fig . 5 in liszt 2001 ) . on the other hand , if the fraction of ly@xmath24 flux from early - type stars that permeates the warm gas is large ( @xmath90.01 ) , the spin temperature is within 10% of the kinetic temperature . however , hi gas at such a low kinetic temperature as 3600 k does not have a stable phase and is outside the allowed range of values in the model proposed by wolfire et al ( 1995 ) ( fig . 5 in liszt 2001 ) . either we are observing an unstable phase of hi , or very small , if any , fraction of ly@xmath24 radiation from early - type stars permeates the diffuse intercloud medium .    the warm hi gas at v@xmath30 30 km s@xmath1 is associated with the perseus arm of the galaxy at a distance of 2.5 kpc ( l@xmath8 160@xmath31 , b@xmath8 10@xmath31 ) . the mean z of the hi gas is thus @xmath8 460 pc . at this height the electron density is @xmath29 0.02 cm@xmath3 ( taylor & cordes 1993 ) , and combined with the current estimate of n@xmath32 implies that the ionization fraction at this height is @xmath29 0.05 . as discussed earlier in this section , the kinetic temperature of the wnm in the direction of 3c147 can be in the range 3600 - 7000 k. the full width at half maximum of the hi 21 cm emission line from the wnm due to thermal broadening alone is 13 - 18 km s@xmath1 . the observed hi emission line width from the wnm in the direction of 3c147 is 34 km s@xmath1 . if the observed excess width is attributed to turbulence in the wnm , the implied full width at half maximum of the turbulence broadening is 31 - 29 km s@xmath1 . the outer arm is the structure observed in hi in the range 50@xmath33 195@xmath31 , 5@xmath34 35@xmath31 , and 175 @xmath29 v@xmath17 @xmath29 60 km s@xmath1 ( wakker , & van woerden 1991 ) . the outer arm appears to lie between 15 and 20 kpc from the galactic center in this longitude range ( kulkarni , blitz , & heiles 1982 ) . at the latitudes of 3c147 and 3c380 the hi emission images from the ldss show well - delineated structures in the longitude - velocity plane which can be associated with the outer arm . in addition , there are discrete hi emission features in these images at the longitudes of 3c147 and 3c380 . corresponding to these features , the hi emission spectra toward 3c147 and 3c380 from ldss ( figs . 5a and 6a ) clearly show lines at v@xmath17 = 120.4 and 109.7 km s@xmath1 respectively ( table 2 ) . the high velocity hi absorption lines detected in the current observations toward 3c147 ( v@xmath17 = 117.3 and 124.5 km s@xmath1 ) and 3c380 ( v@xmath17 = 113.7 km s@xmath1 ) ( figs . 5b and 6b ) are associated with the respective hi emission features arising in the outer arm . as mentioned in section 3 , in both cases ( i.e. , 3c147 and 3c380 ) , the high velocity clouds have at least two temperature components  one with a spin temperature higher than , or comparable to 1000 k , and the second with a spin temperature lower than 200 k. the higher temperature gas , presumably wide - spread , is observed in hi emission . the narrow line hi absorption may arise from smaller and colder regions ( ` clumps ' ) embedded in this bigger structure . there are probably many such ` clumps ' in the high velocity hi emission features since hi absorption was detected toward two out of the four extragalactic sources observed . the heliocentric distances to these high velocity hi clouds toward 3c147 and 3c380 are 10 kpc and 18 kpc respectively , assuming a mean galactocentric distance of 18 kpc to the outer arm ( kulkarni , blitz , & heiles 1982 , blitz , fich , & stark 1982 ) . at these distances , the heights of these high velocity clouds above the galactic plane are 1.8 kpc and 8 kpc respectively . very little is known about the spin temperature of hi in the outer arm due to paucity of hi absorption measurements . payne , salpeter , & terzian ( 1980 ) reported an hi absorption measurement in the outer arm toward 4c33.48 ( l=66.4@xmath31 , b=8.4@xmath31 ) at v@xmath17 = 128 km s@xmath1 with @xmath35 = 0.011 and t@xmath10 = 69 k. more recently , akeson , & blitz ( 1999 ) detected hi absorption in the outer arm toward the extragalactic sources 1901 + 319 and 2005 + 403 . they report @xmath35 = 0.012 and 0.031 and t@xmath10 = 38 and 147 k for these two lines of sight respectively . the widths of these hi absorption lines and the corresponding hi emission lines are between 5 and 12 km s@xmath1 . although the current hi absorption in the outer arm seems to imply the existence of similar spin temperature components , there are two distinct differences compared to the earlier two measurements : ( 1 ) the optical depths observed in the current observations are about a factor of ten smaller than the previous measurements , and ( 2 ) the current measurements indicate the existence of multiple temperature components in the outer arm hi . the existence of multiple temperature components in other parts of the galaxy has been well - established based on earlier hi observations ( mebold et al 1982 , stark et al 1994 ) . it appears that the less explored outer arm is not too different in this respect . however , the outer arm is different from the inner galaxy in at least two aspects - ( 1 ) the column density of atomic hydrogen can be less than or equal to 10@xmath36 @xmath27 , and ( 2 ) the flux of diffuse ly@xmath28 radiation can be significantly lesser . these two factors play an important role in estimating the allowed range of spin temperatures for an allowed range of kinetic temperatures in the two - phase equilibrium models of wolfire et al ( 1995 ) . it will be interesting to compare these estimates of spin temperatures with those obtained from observations toward many lines of sight in the outer arm . we present observations of the galactic hi 21 cm - line absorption toward four extragalactic sources using the wsrt . towards one of the sources we have detected hi 21 cm - line absorption in the warm neutral medium associated with the perseus arm of the galaxy . the observed optical depth is ( 1.9@xmath00.2)@xmath210@xmath3 leading to a spin temperature estimate of 3600@xmath0360 k. depending on the intensity of the diffuse ly@xmath24 radiation field the kinetic temperature of the wnm is in the range 3600 - 7000 k.    we have also detected hi 21 cm - line absorption from high velocity clouds in the outer arm of the galaxy . we find multiple temperature components ( t@xmath37 200 k and @xmath9 1000k ) in these high velocity clouds similar to that observed in other parts of the galaxy . we thank the referee for a thorough reading of an earlier version of this paper and giving detailed comments to improve the clarity and readability of the paper . the national radio astronomy observatory ( nrao ) is a facility of the national science foundation , operated under cooperative agreement by associated universities , inc .. the westerbork synthesis radio telescope is operated by the netherlands foundation for research in astronomy with financial support from the netherlands organization for scientific research ( nwo ) . we thank rene vermeulen for scheduling and carrying out these observations . cccccccc 3c147 & 161.7 + 10.3 & 21.3 & 16.6 @xmath2 12.1 & 2.4 & 3.3 & 1.5@xmath210@xmath7 & 2 @xmath2 15 3c273 & 290.0 + 64.4 & 26.9 & 511 @xmath2 11.7 & 0.4 & 3.2 & 1.2@xmath210@xmath7 & 1 @xmath2 12 3c295 & 97.5 + 60.8 & 21.3 & 28.9 @xmath2 14.4 & 75.0 & 7.7 & 3.6@xmath210@xmath7 & 2 @xmath2 10 3c380 & 77.2 + 23.5 & 12.3 & 46.1 @xmath2 19.7 & 16.1 & 6.4 & 5.2@xmath210@xmath7 & 1 @xmath2 12    rrrrrrrrrr 3c147 & + 2.8(.5 ) & 5.0(.4 ) & .1(19 ) & + 5.6(2 ) & 11.1(2 ) & 8.5(2)&1.8 & 85 & 2.0 & 0.6(.1 ) & 4.6(.1 ) & .38(21 ) & + 0.7(.04 ) & 6.5(.2 ) & 39.5(4.1 ) & 5.0 & 105 & 2.2 & 8.8(.01 ) & 2.9(.04 ) & .29(5 ) & & & & & @xmath29185 & & 10.7(.06 ) & 7.3(.06 ) & .26(3 ) & 10.8(.04 ) & 9.6(.1 ) & 35.9(.2 ) & 6.7 & 140 & 2.5 & 20.2(.07 ) & 2.4(.8 ) & .009(.7 ) & & & & & @xmath29125 & & 29(4 ) & 53(6 ) & .0019(.2 ) & 35.5(.3 ) & 34(1.1 ) & 6.9(.1 ) & 4.6 & 3600 & 12.8 & 34.7(.8 ) & 2.7(3.9 ) & .0013(1.5 ) & & & & & @xmath29160 & & 43.9(.6 ) & 3.5(1.2 ) & .0016(.5 ) & & & & & @xmath29270 & & & & & 81.8(2.2 ) & 29(6 ) & 0.8(.1 ) & 0.5 & @xmath91780 & & & & & 120.4(1.4 ) & 17(3.6 ) & 0.8(.1 ) & 0.3 & @xmath91780 & & 117.3(.3 ) & 8.0(.9 ) & .0013(.1 ) & & & & & @xmath291400 & & 124.5(.1 ) & 3.5(.2 ) & .003(.1 ) & & & & & @xmath29270 & + 3c273 & + 32.2(.08 ) & 3.0(.2 ) & .0038(.2 ) & & & & & @xmath29200 & & + 24.5(.1 ) & 2.6(.3 ) & .0026(.2 ) & & & & & @xmath29150 & & 5.4(.1 ) & 6.6(.6 ) & .0048(1.1 ) & & & & & @xmath29950 & & 5.8(.03 ) & 3.0(.1 ) & .0147(1.1 ) & 5.7(.03 ) & 6.8(.1 ) & 4.6(.06 ) & 0.6 & 300 & 3.7 & & & & 12.7(1.6 ) & 76.0(8 ) & 0.34(.05 ) & 0.5 & @xmath9900 & & & & & 13.2(.5 ) & 24.7(1.4 ) & 0.95(.06 ) & 0.5 & @xmath92600 & + 3c295 & & & & 0.25(.2 ) & 19.1(.4 ) & 1.2(.02 ) & 0.4 & @xmath91100 & & 22.9(.1 ) & 3.5(.3 ) & .0066(.5 ) & 23.3(.03 ) & 3.9(.08 ) & 2.9(.05 ) & 0.2 & 440 & 4.5 & & & & 41.9(.3 ) & 24.9(.7 ) & 0.9(.02 ) & 0.4 & @xmath9800 & + 3c380 & + 8.4(.1 ) & 4.2(.3 ) & .0098(.7 ) & & & & & @xmath29385 & & + 2.3(.02 ) & 2.8(.05 ) & .0526(.8 ) & & & & & @xmath29170 & & & & & + 0.1(.1 ) & 26.6(.2 ) & 7.6(.1 ) & 3.9 & @xmath94200 & & & & & 3.1(.02 ) & 3.5(.07 ) & 10.1(.2 ) & 0.7 & @xmath29270 & & 6.0(.01 ) & 2.7(.03 ) & .1208(1.2 ) & & & & & @xmath29160 & & 8.3(.1 ) & 8.3(.1 ) & .0356(.1 ) & & & & & @xmath291510 & & & & & 109.7(.6 ) & 38.1(1.4 ) & 1.3(.04 ) & 1.0 & @xmath9800 & & 113.7(.3 ) & 3.1(.6 ) & .0043(.8 ) & & & & & @xmath29210 &
Abstract:

Answer: using the westerbork synthesis radio telescope , we have detected hi 21 cm - line absorption in the warm neutral medium of the galaxy toward the extragalactic source 3c147 . this absorption , at an lsr velocity of 29@xmath04 km s@xmath1 with a full width at half maximum of 53@xmath06 km s@xmath1 , is associated with the perseus arm of the galaxy . the observed optical depth is ( 1.9@xmath00.2)@xmath210@xmath3 . the estimated spin temperature of the gas is 3600@xmath0360 k. the volume density is 0.4 cm@xmath3 assuming pressure equilibrium . toward two other sources , 3c273 and 3c295 , no wide hi 21 cm - line absorption was detected . the highest of the 3@xmath5 lower limits on the spin temperature of the warm neutral medium is 2600 k. in addition , we have also detected hi 21 cm - line absorption from high velocity clouds in the outer arm toward 3c147 and 3c380 at lsr velocities of 117.3 , 124.5 and 113.7 km s@xmath1 respectively . we find two distinct temperature components in the high velocity clouds with spin temperatures of greater than 1000 k and less than 200 k , respectively .