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 phase structure of quantum chromodynamics(qcd ) in the presence of an external magnetic field has been explored extensively in recent years@xcite . this is of great phenomenological relevance for non - central heavy ion collisions and dense neutron stars as well as cosmological evolutions in its early stage . in particular , the peculiar response of the spontaneously breaking of chiral symmetry to a strong magnetic field is among the central issues and posseses a theoretical challenge . the profound effect of the magnetic field on chiral symmetry breaking had been pointed out in quantum electrodynamics(qed ) very early@xcite . it was found that magnetic fields can enhance the chiral symmetry breaking due to an effective dimension reduction in the lowest landau level(lll ) in the charged fermionic sector . this phenomenon is called _ magnetic catalysis_. in qcd , this problem can be addressed reliably by the lattice monte - carlo simulation without sign problem@xcite . it has shown intriguingly and , however , unexpectedly that the chiral critical temperature decreases with an increasing magnetic field , which is in explicit conflict with magnetic catalysis and thus termed as _ inverse magnetic catalysis _ or _ magnetic inhibition_. many attempts already exist in the literatures trying to explain this contradiction@xcite .    along with the dimension reduction in the fermionic sector , which accounts for the magnetic catalysis , the significance of magnetic fields manifests in another important aspect in a general way in systems with spontanously symmetry breaking . it enhances the fluctuations through the spatial anisotropy of the spectrum of bosonic modes in the system , although the bosons are still in @xmath0 dimension and do not directly suffer from the mermin - wagner - coleman theorem which forbids the spontaneously breaking of continuous symmetry in @xmath1 and @xmath2 dimensions@xcite . this claim follows by the observation that the ginzburg critical window of the chiral phase transition gets widened in the presence of a magnetic field@xcite . the role played by magnetic fields is thus twofold which competes with each other . a simple and straightforward model to investigate this competition is the bose - einstein condensation(bec ) of neutral bosons composed by two equally and oppositely charged fermions@xcite , which is a physical analog of chiral condensate in a magnetic field , yet analytically tractable . the corresponding non - relativistic case had been studied in @xcite , where we found that the critical temperature for the bec was dramatically affected by the magnetic field . in weak coupling domain ( where the coupling fails to support a two fermion bound state in the absence of a magnetic field ) , dimension reduction dominates and we found the usual magnetic catalysis effect . in strong coupling domain ( where a two fermion bound state exists in the absence of a magnetic field ) , however , the fluctuations dominate and the situation gets reversed , we found the inverse magnetic catalysis , i.e. , the critical temperature of bec decreases as an increasing magnetic field .    in the present paper , we shall extend the above analysis to the relativistic fermions . technically , this amounts to replace the non - relativistic fermionic propagator by the dirac one in a magnetic field@xcite . the corresponding one loop diagram underlying the theory will be quadratically divergent . the leading divergence can be removed by the coupling renormalization , similar to the non - relativistic case , where the renormalized coupling is in terms of the scattering length , but the logarithmic subleading divergence remains and a uv cutoff has to be introduced by hand . consequently , the results will be explicitly cutoff dependent , like all relativistic field theories with 4-fermion interactions . nevertheless , the cutoff does not change the results qualitatively . we found that the low lying spectrum of the bound pairs is given by @xmath3 where @xmath4(@xmath5 ) is the momentum component parallel ( perpendicular ) to the magnetic field , @xmath6 is an anisotropicity factor and @xmath7 is the rest energy of a bound pair . this spectrum implies a nonrelativistic like condensation temperature @xmath8^{2/3}\frac{\pi}{\bar\mu}\label{transitiontemperature}.\ ] ] at sufficiently low density . unlike the nonrelativistic case where only the anisotropicity factor @xmath9 is magnetic field dependent , the `` mass '' @xmath10 here depends on the magnetic field as well . in the weak coupling domain , the magnetic field dependence of both @xmath9 and @xmath10 support magnetic catalysis and @xmath11 increases with an increasing magnetic field as in the nonrelativistic case . the situation in strong coupling domain , however , is very different from that in the nonrelativistic case , in which the inverse magnetic catalysis prevails regardless of the strength of the magnetic field . in the present case , the inverse magnetic catalysis only appears when magnetic field is low and the magnetic catalysis effect takes place again when magnetic field is sufficiently high because of opposite roles played by @xmath9 and @xmath10 there . consequently , @xmath11 decreases first and then increases with the magnetic field . the rest of the paper is organized as follows : in section ii we lay out the general formulation and present the mean field approximation . the fluctuations beyond the mean field theory , which is necessary for bec , is calculated under the gaussain approximation in section iii . the numerical solutions to the gap equation and the condensation temperature are presented in section iv . section v is devoted to the conclusions and outlooks . some calculation details and useful formulas are presented in the appendices a , b , c and d. throughout the paper , we will work in euclidean space with the four vector represented by @xmath12 with @xmath13 the matsubara frequency for bosons @xmath14 or for fermions @xmath15 with @xmath16 the temperature . we consider a system consisting of relativistic fermions of mass @xmath17 and chemical potential @xmath18 with opposite electric charge interacting through a short ranged instantaneous attractive interaction . in a quark matter , the attractive interaction stems from the nonperturbative qcd effects and the long range coulomb interaction is perturbative and can be ignored in the leading order . the lagrangian density reads @xmath19 where @xmath20 and @xmath21 with @xmath22 the covariant derivative and @xmath23 the charge magnitude . we also defined the charge conjugate operator @xmath24 and @xmath25 . in the following , we shall consider the situation of a constant magnetic field and work in landau gauge , in which the vector potential @xmath26(the field @xmath27 being along the @xmath28 direction ) . introducing the usual hubbard - stratonovich field @xmath29 coupled to @xmath30 , we obtain the partition function @xmath31.\ ] ] where @xmath32 is the normalization constant and the nambu - gorkov(ng ) spinors are defined as @xmath33 the effective action is given by @xmath34,\ ] ] with    @xmath35    the thermodynamic potential reads @xmath36 and the number density of fermions is given by @xmath37 with @xmath38 the pressure and @xmath39 the spatial volume of the system . integrating out the ng degrees of freedom , the partition function becomes @xmath40\right),\ ] ] with @xmath41-\int d^4 x\frac{|\delta(x)|^2}{g}.\label{effectiveaction}\ ] ] where the trace runs over space , imaginary time , dirac and ng indices . the mean field approximation ignores the fluctuations in @xmath29 field , which amounts to set @xmath42 . in this approximation , we obtain @xmath43 ^ 2\right ) \left((i\omega_n)^2-\left[e_l^-(\delta_0)\right]^2\right)\right],\end{aligned}\ ] ] where @xmath44 and @xmath45 labels the landau levels ( ll ) . notice that , the factor @xmath46 reflects the non - degeneracy in the lowest landau level ( lll ) . we defined @xmath47^{1/2}.\ ] ] which are the energy spectrum of the fermionic exitations in the system . the thermodynamic potential under mean - field approximation is @xmath48 . \label{pressure0}\ ] ]    from the saddle point condition @xmath49}{\delta\delta_0}=0,\ ] ] one obtains the gap equation @xmath50.\label{gapequation}\end{aligned}\ ] ] and @xmath51 at the transition temperatute @xmath52 . the number density under mean field approximation follows from eq.([fermionnumberdensity ] ) and reads @xmath53}{\partial\mu}\right)_t . \label{fermionnumberdensity0}\ ] ]    in bcs limit , the critical temperature coincides with the threshold pairing temperature which , at a given chemical potential , is determined by ( [ gapequation ] ) at @xmath51 first and the chemical potential is in turn solved by the density equation ( [ fermionnumberdensity ] ) . in bec limit , however , the roles of ( [ gapequation ] ) and ( [ fermionnumberdensity ] ) reversed because the critical tempearture becomes separate from the pairing . ( [ gapequation ] ) determines the chemical potential@xcite . the bec critical temperature is determined by eq.([fermionnumberdensity ] ) equals ( [ fermionnumberdensity0 ] ) plus the contribution from the bosonic excitation spectrum , whose low momentum behavior will be given in the next section . there we shall also see that the quantity @xmath54 with @xmath18 the solution of the gap equation measures the binding energy of a relativistic pair and locates the bosonic pole below the two - fermion cut at zero total momentum . to make the formulation analytically tractable , we assume that @xmath55 and @xmath56 the gap equation ( [ gapequation ] ) can be approximated by @xmath57.\label{zerotemgapeq}\end{aligned}\ ] ] where we have set @xmath51 , which amounts to approach the critical point from normal to superconducting phase . obviously , this equation is quadratically divergent and needs to be regularized . to maintain the relativistic invariance , we shall use the pauli - villars regularization scheme . introducing the pv regulators , we obtain @xmath58.\label{pvgapeq}\end{aligned}\ ] ] with @xmath59 and @xmath60 . throughout the paper , we shall assume that the regulator mass @xmath61 but finite . carrying out the integral over the momentum under the constrains imposed by the regularization conditions , we obtain @xmath62.\label{gapeqinappendix}\end{aligned}\ ] ] dropping all terms with negative powers of @xmath63 since @xmath64 by definition , we end up with    @xmath65\right.\\ \nonumber&+\left.2\sum_{l=0}^{\infty}\left[\frac{\mu}{\sqrt{2leb+m^2-\mu^2}}\tan^{-1}\frac{\mu}{\sqrt{2leb+m^2-\mu^2}}-\frac{\mu^2}{2leb+m^2-\mu^2}\right]\right\}\\ & + \frac{1}{4\pi^2}\left[m^2\sum_{s>0}c_s\lambda_s^2\ln\lambda_s^2 - 2\mu^2\sum_{s>0}c_s\ln\lambda_s^2+m^2\left(1+\ln\frac{2eb}{m^2}\right)-2\mu^2\ln\frac{2eb}{m^2}-eb\ln\frac{4\pi eb}{m^2}\right].\label{regularizedgapeq}\end{aligned}\ ] ]    with @xmath66 and @xmath67 the gamma function , and @xmath68 . the detailed steps leading from ( [ gapeqinappendix ] ) to ( [ regularizedgapeq ] ) will be shown in appendix a. the quadratic term in the regulator masses in the last line can be absorbed into the coupling constant as usual and the logarithmic term in the regulator ones remains , leaving the gap equation explicitly depends on an uv cutoff . in the appendix b , we have proved that solutions to this gap equation is a decreasing function of an increasing magnetic field . the gap equation at zero magnetic field can be obtained straightforwardly from ( [ regularizedgapeq ] ) and reads @xmath69 where the critical coupling @xmath70 had absorbed the quadratically divergent term . the strong and weak coupling domain is characterized by a critical coupling , which is the threshold coupling for a bound state to show up in the absence of a magnetic field . it follows then the critical coupling constant is defined as the one for @xmath71 . moreover , we shall denote @xmath72 with @xmath73 an uv cutoff . we thus have @xmath74 this definition establishs an explicit dependence of the dimensionless critical coupling constant @xmath75 on the uv cutoff . in this section , we shall include all nonzero momentum component of the hubbard - stratonovich field to the quadratic order ( gaussian fluctuations ) in the partition function . approching the condensate from the normal phase , we shall reproduce the gap quation ( [ zerotemgapeq ] ) by the thouless critirion and obtain the contribition from the bosonic excitations to the number density ( [ fermionnumberdensity ] ) , which is necessary to determine the critical temperature in bec limit . for the simplicity of notation , we shall not explicitly write out the pauli - villars regulators in the following calculations . but they will be restored whenever the integral is divergent . expanding the bosonic field @xmath29 around the trivial saddle point @xmath51 in the effective action ( [ effectiveaction ] ) , we obtain that @xmath76&\simeq { \cal s}_{\rm eff}[0]-\frac{1}{g}\int d^4 x|\delta(x)|^2\\ & + { \rm tr}\left[g_0 ^ -(x , y)\gamma^5\delta^*(y ) g_0^+(y , x)\gamma^5 \delta(x)\right],\label{guanssianaction}\end{aligned}\ ] ] with @xmath77^{-1}$ ] the diagonal elements in the ng propagator ( [ ngpropagator ] ) , which are nothing but the inverse dirac operators in a magnetic field in coordinate space . the dirac propagator had been firstly obtained by schwinger using the proper time method . here we shall employ the form with explicit landau levels dependence , given in@xcite @xmath78^{-1}\delta(x - y)\\ = & e^{i\phi(x , y)}t\sum_n\int\frac{d^3{\bf k}}{(2\pi)^3}e^{-ik\cdot z}\sum_{l=0}^{\infty}g_l^\pm(\omega_n,{\bf k}),\label{diracpropagator}\end{aligned}\ ] ] with    @xmath79\left[2\pi_+l_l\left(\frac{2{\bf k}_\perp^2}{eb}\right)-2\pi_-l_{l-1}\left(\frac{2{\bf k}_\perp^2}{eb}\right)\right]\right.\\ & + \left.4({\bf k}\cdot\gamma)_\perp l_{l-1}^1\left(\frac{2{\bf k}_\perp^2}{eb}\right ) \right\}.\end{aligned}\ ] ]    where @xmath80 and we also defined the projectors @xmath81 , @xmath82 is the generalized laguerre polynomials(@xmath83 and @xmath84 by definition ) . the gauge dependent phase factor is given by @xmath85 which cancels out in the combination of ( [ guanssianaction ] ) . consequently , the kernel mediating @xmath86 and @xmath29 there is explictly translational invariant , i.e. , depends only on @xmath87 , as expected .    in momentum space , ( [ guanssianaction ] ) becomes @xmath88={\cal s}_{\rm eff}[0]-t\sum_n\int \frac{{d^3\bf p}}{(2\pi)^3}\gamma^{-1}(\omega_n,{\bf p})|\delta(\omega_n,{\bf p})|^2,\label{gaussianeffectiveaction}\ ] ] where the dependence of the coefficient @xmath89 on @xmath90 and @xmath27 is implicit . the instability towards the superconducting phase is then signaled by the condition @xmath91 which is just the thouless criterion for the onset of a long range order . the number density equation is determined by @xmath92 where @xmath93 $ ] with @xmath94 the fermi - dirac distribution function , is the fermionic contribution defined in ( [ fermionnumberdensity ] ) . ( [ thoulesseq ] ) and ( [ gaussiannumbereq ] ) are two basic equations to solve for bcs / bec crossover@xcite . following the convention in@xcite , one can write the number density ( [ gaussiannumbereq ] ) in terms of a phase shift defined by @xmath95 $ ] . we have @xmath96 with @xmath97 $ ] the bose - einstein distrubution function .    to calculate @xmath89 , we make use of the eq.([diracpropagator ] ) and the fourier transformation of the fluctuation field @xmath98 and obtain    @xmath99\left[l_{l^\prime}\left(\frac{2{\bf p}_\perp^2}{eb}\right)l_l\left(\frac{2{\bf q}_\perp^2}{eb}\right)+l_{l^\prime-1}\left(\frac{2{\bf p}_\perp^2}{eb}\right)l_{l-1}\left(\frac{2{\bf q}_\perp^2}{eb}\right)\right]\right.\\ & \left.-8\left({\bf k}\cdot{\bf q}\right)_\perp l^1_{l^\prime-1}\left(\frac{2{\bf p}_\perp^2}{eb}\right)l^1_{l-1}\left(\frac{2{\bf q}_\perp^2}{eb}\right)\right\}-\frac{1}{g}.\end{aligned}\ ] ]    with @xmath100 . working out the sum over matsubara frequecies , one can verify that the condition ( [ thoulesseq ] ) will yield the same gap equation as ( [ zerotemgapeq ] ) . let s proceed by employing the proper - time representation of the propagators introduced by schwinger @xmath101 where the feynmann parametrization scheme had been used . we obtain    @xmath102+(1-x)\left[(\omega'+\omega+2\mu)^2-q_3 ^ 2 - 2leb - m^2 \right ] ) } \\ \nonumber&\times\left\ { \left[\omega^{\prime}(\omega^{\prime}+\omega+2\mu)-k_3q_3-m^2\right]\left[l_{l^\prime}\left(\frac{2{\bf p}_\perp^2}{eb}\right)l_l\left(\frac{2{\bf q}_\perp^2}{eb}\right)+l_{l^\prime-1}\left(\frac{2{\bf p}_\perp^2}{eb}\right)l_{l-1}\left(\frac{2{\bf q}_\perp^2}{eb}\right)\right]\right.\\ & \left.-8\left({\bf k}\cdot{\bf q}\right)_\perp l^1_{l^\prime-1}\left(\frac{2{\bf p}_\perp^2}{eb}\right)l^1_{l-1}\left(\frac{2{\bf q}_\perp^2}{eb}\right)\right\}-\frac{1}{g}.\label{inversegamma}\end{aligned}\ ] ]    where we analytically continuated the matsubara frequency to the continuous one as @xmath103 and shift the frequency from @xmath104 to @xmath105 . the same approximation as in the preceding section , @xmath106 has been made in ( [ inversegamma ] ) that turns the matsubara sum over @xmath107 to the integral over @xmath108    using the generating function of laguerre polynomials @xmath109 one can carry out the sum over landau levels and has    @xmath110}\\ & \times \left\ { \left[k\cdot q - m^2\right]\frac{1+tt'}{(1-t)(1-t')}-8\left({\bf k}\cdot{\bf q}\right)_\perp\frac{tt'}{(1-t)^2(1-t')^2 } \right\}-\frac{1}{g}.\label{covariantgamma}\end{aligned}\ ] ]    with @xmath111 and @xmath112 . we defined @xmath113 , @xmath114 and @xmath115 , by which eq.([covariantgamma ] ) can be written in a covariant form . one may observe that all momentum integrals become gaussian ones , which is the advantage of using the proper - time representations . the integral over @xmath116 can be carried out by shifting the variable to @xmath117 and wick rotating @xmath118 with @xmath119 from @xmath120 to @xmath121 , and we obtain    @xmath122}\\ & \times \left\ { \left[\frac{1}{\lambda}+x(1-x)p^2+m^2\right]\frac{1+tt'}{(1-t)(1-t')}+8\left({\bf k}\cdot{\bf q}\right)_\perp\frac{tt'}{(1-t)^2(1-t')^2 } \right\}-\frac{1}{g}.\end{aligned}\ ] ]    the other integral over momentum @xmath123 can be performed by noticing that @xmath124 and @xmath125 we end up with    @xmath126 -\frac{\tanh(\lambda xeb)\tanh(\lambda(1-x)eb)}{\tanh(\lambda xeb)+\tanh(\lambda(1-x)eb)}\frac{{\bf p}_\perp^2}{eb}}\left\ { \left[\frac{1}{\lambda}+x(1-x)p^2+m^2\right]\coth(\lambda eb)\right.\\ & + \left.\frac{1}{\sinh^2(\lambda eb ) } \left[eb -\frac{\tanh(\lambda xeb)\tanh(\lambda(1-x)eb)}{\tanh(\lambda xeb)+\tanh(\lambda(1-x)eb)}{\bf p}_\perp^2\right ] \right\}-\frac{1}{g}.\label{finalgamma}\end{aligned}\ ] ]    the singularity structure of @xmath127 in the entire complext @xmath128-plane reflects the two - particle spectrum . there will be an isolated real pole representing the two - body bound pair and a branch cut along the real axis representing the continuum of two - particle excitations . for sufficiently strong coupling and low temperature , the contribution to the number density is dominated by the bound pair pole , which is determined by @xmath129 with @xmath130 the solution to the mean field equation ( [ regularizedgapeq ] ) . we henceforth consider the expansion of ( [ finalgamma ] ) around this pole and obtain @xmath131 with @xmath132}\\ & \times\left\{\left[\frac{2}{\lambda}+4x(1-x)\bar{\mu}^2+m^2\right]\coth(\lambda eb)+\frac{eb}{\sinh^2(\lambda eb ) } \right\ } , \label{coeffa}\end{aligned}\ ] ] and @xmath133}\\ \nonumber&\times\frac{\tanh(\lambda xeb)\tanh(\lambda(1-x)eb)}{\tanh(\lambda xeb)+\tanh(\lambda(1-x)eb)}\\ & \times \left\{\left[\frac{1}{\lambda}+4x(1-x)\bar{\mu}^2+m^2\right]\coth(\lambda eb)+\frac{2eb}{\sinh^2(\lambda eb)}\right\}. \label{coeffb}\end{aligned}\ ] ] notice that , the kinetic term in ( [ gammaundergaussian ] ) becomes anisotropic with respect to the directions along and perpendicular to the magnetic field due to the rotational symmetry breaking by the magnetic field . one can readily show that @xmath134 at zero magnetic field as expected . moreover , the coefficients satisfy the inequality @xmath135 regardless of solutions to the gap equation(see appendix d for details ) . the bosonic spectrum at low momentum is determined by the pole of @xmath127 and reads @xmath136 with the anisotropicity factor @xmath137 in contrast to the non - relativistic case , in which the mass of the pair is @xmath138 , the mass is replaced by the @xmath139 , which is magnetic field dependent . once @xmath140 becomes smaller than @xmath17 , stable bosonic pairs appear and they dominate the total number density at sufficienly low temperature . the condensation temperature is determined by setting the chemical potential at the solution of the mean field equation ( [ regularizedgapeq ] ) , i.e. @xmath141 , and the phase factor follows @xmath142 where @xmath143 is the heaviside step function . substituting ( [ phasefactor ] ) into ( [ numberdensity ] ) , we obtain @xmath144^{-1 } , \label{beccondition}\ ] ] where @xmath145 the fermionic contribution had been ignored under the condition @xmath146 of our approximation . solving ( [ beccondition ] ) for @xmath11 , we end up with the formula ( [ transitiontemperature ] ) reported in section i. in this section , we shall explore the bec critical temperature in view of ( [ regularizedgapeq ] ) and ( [ transitiontemperature ] ) . the results display a nontrivial dependence on the strength of the coupling constant and magnetic field . the strong coupling and weak coupling domains correspond to @xmath147 and @xmath148 respectively . for the discussion below , it is convenient to introduce a dimensionless critical temperature @xmath149 where @xmath150 is given by ( [ transitiontemperature ] ) with @xmath151 and @xmath152 . in the strong coupling domain , @xmath153 is the solution of the gap equation at zero magnetic field and @xmath150 is the critical temperature of a bound pair at zero magnetic field . in the weak coupling domain , we set @xmath154 and @xmath150 the critical temperature of a nonrelativistic boson with mass @xmath138 at zero magnetic field . the ( inverse)magnetic catalysis implies increasing(decreasing ) @xmath155 with magnetic field . let us consider first the weak magnetic field limit @xmath156 , where analytic method can be applied . in the strong coupling domain , we expand the rhs of the gap equation to the quadratic order in @xmath157 with the aid of stirling formula for the gamma and di - gamma functions in ( [ regularizedgapeq ] ) and the euler - maclaurin formula for the infinite series there and obtain @xmath158 with the solution @xmath159,\ ] ] where @xmath160 and @xmath161 is the derivative of rhs of ( [ regularizedgapeq ] ) with respect to @xmath18 at @xmath162 and @xmath163 , given by @xmath164 the anisotropicity in this approximation reads @xmath165 with    @xmath166}\left\{\frac{1}{\lambda}(2+x-6x^2 + 3x^3)+x(1-x)^2\left[4x(1-x)\bar{\mu}_0 ^ 2+m^2\right]\right \}\\ & = \frac{3m^2 - 2\bar{\mu}_0 ^ 2}{32\bar{\mu}_0 ^ 2}-\frac{9m^4 - 28m^2\bar{\mu}_0 ^ 2 + 16\bar{\mu}_0 ^ 4}{96\bar{\mu}_0 ^ 3\sqrt{m^2-\bar{\mu}_0 ^ 2}}\tan^{-1}\frac{\bar{\mu}_0}{\sqrt{m^2-\bar{\mu}_0 ^ 2}},\end{aligned}\ ] ]    and @xmath167}\\ \nonumber&\times\left[\frac{3}{\lambda}+4x(1-x)\bar{\mu}_0 ^ 2+m_s^2\right ] \\ = & \frac{1}{2}\ln\lambda^2-\frac{m^2 - 2\bar{\mu}_0 ^ 2}{2\bar{\mu}_0\sqrt{m^2-\bar{\mu}_0 ^ 2}}\tan^{-1}\frac{\bar{\mu}_0}{\sqrt{m^2-\bar{\mu}_0 ^ 2}}+\frac{2}{3}.\end{aligned}\ ] ] it follows that the dimensionless critical temperature @xmath168 with the coefficient @xmath169 the numerical values of the coefficient @xmath170 are tabulated in table i for different couplings and cutoffs in the strong coupling domain and indicate an inverse magnetic catalysis at a weak magnetic field .     the scaled bound state mass versus the dimensionless magnetic field in strong coupling domain . the solid and dashed line corresponds to an uv cutoff @xmath171 and @xmath172 , respectively.,height=264 ]     the scaled bound state mass versus the dimensionless magnetic field in weak coupling domain . the solid and dashed line corresponds to an uv cutoff @xmath171 and @xmath172 , respectively.,height=264 ]    coming to the weak coupling domain , the bound magnetic field is the only catalyst of the bound state and a weak magnetic field triggers a weak bound state , nonrelativistic formulas in @xcite can be carried over in the weak field limit . @xmath173 and @xmath174 with the scattering length @xmath175 extracted from @xmath176 the lll approximation works here and the usual magnetic catalysis emerges ( increasing with @xmath27 ) . the details of the non - relativistic approximation of the gap equation can be found in appendix d.    .[tab : table]numerical values of @xmath170 for different couplings and cutoffs in the strong coupling domain . [ cols="^,^,^",options="header " , ]      the anisotropicity @xmath9 versus the dimensionless magnetic field in strong coupling domain . the solid and dashed line corresponds to solutions to the gap equation with an uv cutoff @xmath171 and @xmath172 , respectively.,height=264 ]     the anisotropicity @xmath9 versus the dimensionless magnetic field in weak coupling domain . the solid and dashed line corresponds to solutions to the gap equation with an uv cutoff @xmath171 and @xmath172 , respectively.,height=264 ]     the ratio @xmath177 versus the dimensionless magnetic field in strong coupling domain . the solid and dashed line corresponds to solutions to the gap equation with an uv cutoff @xmath171 and @xmath172 , respectively.,height=264 ]     the ratio @xmath177 versus the dimensionless magnetic field in weak coupling domain . the solid and dashed line corresponds to solutions to the gap equation with an uv cutoff @xmath171 and @xmath172 , respectively.,height=264 ]     the condensation temperature versus the dimensionless magnetic field in strong coupling domain . the solid and dashed line corresponds to solutions to the gap equation with an uv cutoff @xmath171 and @xmath172 , respectively.,height=264 ]     the condensation temperature versus the dimensionless magnetic field in weak coupling domain . the solid and dashed line corresponds to solutions to the gap equation with an uv cutoff @xmath171 and @xmath172 , respectively.,height=264 ]    beyond the weak field limit , we solve the gap equation ( [ regularizedgapeq ] ) numerically at two different coupling domains that are mildly above and blow the critical coupling with two different sets of uv cutoff and the results are presented in fig.1 and fig.2 . in both cases , the chemical potential at the critical temperature decreases with an increasing magnetic field indicating an increasing binding energy of the bound pairs , consistent with the magnetic catalysis . in strong coupling domain , however , bound pairs may appear at zero magnetic field in contrast to the case in weak coupling domain , in which the pairs can only support by magnetic catalysis . therefore , the curve in fig.1 starts from @xmath178 at @xmath162 while the curve in fig.2 from @xmath179 at @xmath162 . in both cases , the results display an explicit dependence on the uv cutoff , which is intrinsic in the ( 3 + 1 ) dimensional model with contact interaction . but the dependence does not change the results qualitatively . the reason for the cross between the solid and dashed curves at lower magnetic field can be explained as follows . one can show that @xmath180^{-1}.\label{muoflambda}\ ] ] with @xmath181 and @xmath182 the rhs of the gap equation ( [ regularizedgapeq ] ) . it had been proved in the appendix that @xmath183 and ( [ muoflambda ] ) thus changes sign at a critical value @xmath184 . one therefore expects a crossing of the two sets of solutions corresponding to different uv cutoffs at strong coupling domain with @xmath185 as shown in fig.1 . such crossing should be absent in weak coupling since @xmath186 . the magnetic field dependence of the dimensionless critical temperature ( [ dimltc ] ) is determined by that of the anosotropicity factor @xmath9 and the ratio @xmath187 , which are plotted in fig.3 to fig.8 . as we can see , the anisotropicity decreases ( increases ) with magnectic field in strong ( weak ) coupling domains while the ratio @xmath188 increases with magnetic field in both domains . consequently , the competetion between the two opposite field dependence of @xmath189 and @xmath188 in the strong coupling domain yields an inverse magnetic catalysis at lower magnetic field but magnetic catalysis at higher magnetic field . in the weak coupling domain , however , both factors increases with the magetic field and the usual magnetic catalysis emerges . in this paper , we have investigated the bec of bound pairs composed by two equally and oppositely charged relativistic fermions in the presence of a magnetic field . this is inspired by the inverse magnetic catalysis effect found in lattice qcd and is a generalization of our non - relativistic work @xcite . the purpose is to explore the interplay between the pairing dynamics and the fluctuations of bound pairs in a magnetic field , the latter of which underlies the mermin - wigner - coleman theorem on the absence of a long range order in lower dimensions . a nambu - jona - lasinio(njl)-like pairing interaction is considered and the fluctuation is examined within the gaussian approximation , which is valid at a sufficiently low density . unlike the case without a magnetic field where a bound state emerges for the coupling strength above a certain threshold , the bound states occurs also in the weak coupling domain because of the @xmath190 nature of the fermionic spectrum within lll . the binding energy of a bound pair is found to increase with the magnetic field as a consequence of the magnetic catalysis . the dependence of the bec critical temperature on the magnetic field depends on the coupling strength . in weak coupling domain , we found the usual magnetic catalysis . in strong coupling domain , unlike the nonrelativitic case , the inverse magnetic catalysis appears only in small magnetic field region and magnetic catalysis takes place again when magnetic field becomes sufficiently large . the reason for this difference is the dependence of the mass of bound pairs on the magnetic field in relativistic case . the condensation temperature is thus determined by two factors , the anisotropy and the ratio between the mass of bound pairs at zero magnetic field and the one with magnetic field dependence , both depends on the magnetic field . in strong coupling domain , the two factors have opposite field dependence and yield the inverse magnetic catalysis in small magnetic filed region but magnetic catalysis at large field . in weak coupling domain , both factors increase with the magnetic field and magnetic catalysis prevails . this variation indicates that the enhanced fluctuations by the anisotropy in the kinetic terms of the bound pairs is less dominant as compared with the nonrelativistic case . although the system we have studied shares the same physics as the chiral condensate in qcd regarding general prperties of a long range order and its fluctuations and is analytically tractable . the qcd dynamics , however , are expected to be far more complicated and other mechanisms may also contribute the simulated magnetic field dependence of the condensation temperature . for example , the quark loop that is influenced by the magnetic field can impact on the running of qcd coupling . in recently works based on functional renormalization group studies@xcite , inverse magnetic catalysis at a low magnetic field and a delayed magnetic catalysis at a high field had been found . the magnetic field at the turning point ( minimum @xmath52 ) between decreasing and increasing @xmath52 is expected to be of order of several @xmath191 . ads / qcd model suggests the same behavior and similar value of magnetic field at the turning point@xcite . it is shown by one lattice calculation@xcite , however , that even at an extremely strong field @xmath192 , inverse magnetic catalysis still prevails . the reconciliation is still lacking at the present stage . in addition to the theoretical values of our present study , the physics involved may also be applicable in cold dense quark matter , such as the color - flavor - locking superconducting phase and the planar phase in single flavor superconductivity@xcite . they are of interests in dense neutron stars , where magnetic field up to @xmath193 to @xmath194 appears on the surface and even higher magnetic field exists in the interior . this research is partly supported by the ministry of science and technology of china ( mstc ) under the `` 973 '' project no . 2015cb856904(4 ) . b. f. is supported by nsfc under grant no . 11305067 , no . 11535005 and the fundamental research funds for the central universities ,  hust : no . hou and h - c ren are partly supported by nsfc under grant nos . 11375070 , 11135011 and 11221504 . p - p wu is supported by the exploratory fund of henan polytechnic university , no . 13156041 .      in this appendix , the reduction of the gap equation ( [ gapeqinappendix ] ) to the form ( [ regularizedgapeq ] ) shall be presented . the first term inside the brackets of ( [ gapeqinappendix ] ) can be evaluated as @xmath195.\end{aligned}\ ] ] using the infinite product representation of gamma function @xmath196,\ ] ] one has @xmath197 therefore the gap equation ( [ gapeqinappendix ] ) becomes @xmath198\label{intermediate}.\end{aligned}\ ] ] with @xmath199 it can be further simplified for @xmath200 and @xmath201 by eliminating the terms that vanishes in the limit @xmath202 with @xmath203 since @xmath204 . ] . for this purpose , we rewrite ( [ intermediate ] ) as @xmath205\right\}.\end{aligned}\ ] ] we have @xmath206 and the other infinite series converges for each @xmath207 , which vanishes in the limit @xmath208 for @xmath203 and can be ignored . applying the stirling formulas @xmath209 and @xmath210 with @xmath211 the bernoulli numbers , @xmath212 . we thus have , for @xmath203 @xmath213 and @xmath214 we end up with the form of ( [ regularizedgapeq ] ) of the gap equation .      in this appendix , we shall provide analytical proofs of the monotinic dependence of the bound state mass and the anisotropy factors on the magnetic field . because of the complication introduced by pauli - villars regulators , we limit our scope within the approximaton that maintains only the logarithmic terms in the regulator masses .    denoting the rhs of the gap equation ( [ gapeqinappendix ] ) by @xmath215 , we have _ lemma _ : ( 1 ) for a given @xmath27 , @xmath216 is a monotonically increasing function of @xmath18 for @xmath217 ; ( 2 ) for a given @xmath18 within @xmath218 , @xmath216 is a monotonically increasing function of @xmath27 .    proof : ( 1 ) the partial derivative of @xmath216 with respect to @xmath18 reads : @xmath219 with @xmath220 defined in ( [ thetas ] ) and @xmath221 and @xmath222 . let us divide the above derivative into two terms , @xmath223 and @xmath224 , with @xmath224 carrying the logarithmic uv divergence and @xmath223 not . we have @xmath225\\ \simeq&\frac{eb\mu}{\pi^2}\left[\sum_{l=0}^\infty \alpha_l\left(\frac{1}{2leb+m^2-\mu^2}-\frac{1}{2leb+m^2}\right)\right]\ge 0,\end{aligned}\ ] ] and @xmath226.\ ] ] where we have assumed that the terms with negative power of the regulator masses dropped from @xmath223 are too small to offset the direction of the inequality . it follows from the binet formula @xmath227 and its derivative @xmath228 that @xmath229 . consequently @xmath230    \(2 ) taking the derivative of ( [ gapeqinappendix ] ) with respect to @xmath27 , we find that @xmath231 where @xmath232\\ \nonumber\ge&\frac{e\mu}{4\pi^2}\sum_{l=0}^\infty\alpha_l\left[\frac{leb+m^2-\mu^2}{(2leb+m^2-\mu^2)^{\frac{3}{2 } } } \sin2\theta_0\right.\\ \nonumber&-\left.\frac{leb\mu}{(2leb+m^2)(2leb+m^2-\mu^2)}\right]\\ = & \frac{e\mu^2}{2\pi^2}\sum_{l=0}^\infty\alpha_l\frac{m^2-\mu^2}{(2leb+m^2)}\ge 0,\end{aligned}\ ] ] and @xmath233\nonumber\\   = & \frac{1}{2\pi^2}\frac{\partial}{\partial b}\left\{eb\left[\ln\gamma\left(\frac{m^2}{2eb}\right)-\frac{m^2}{2eb}\left(\ln \frac{m^2}{2eb}-1\right)\right.\right.\nonumber\\ & -\left.\left.\frac{1}{2}\ln\frac{4\pi eb}{m^2}\right]\right\}\ge 0.\end{aligned}\ ] ] with the inequality following from ( [ b1 ] ) and ( [ b2 ] ) . the lemma is proved . for the solution of the gap equation , @xmath234 and we arrive at _ theorem _ : the bound state mass is a decreasing function of the external magnetic field . @xmath239 \frac{ebs's\cosh((s'+s)eb)}{s'+s}+\frac{(eb)^2s's}{(s'+s)\sinh((s'+s)eb ) } -\left[\frac{1}{s'+s}+\frac{4s's\bar{\mu}^2}{s'+s}+m^2\right]\\ \nonumber & \times\frac{\tanh ( s'eb)\tanh ( s eb)}{\tanh ( s'eb)+\tanh ( s eb ) } -\frac{2eb\tanh ( s'eb)\tanh ( s eb)}{\sinh((s'+s)eb)(\tanh ( s'eb)+\tanh ( s eb))}\nonumber\\ \nonumber = & \left[\frac{1}{s'+s}+\frac{4s's\bar{\mu}^2}{s'+s}+m^2\right ] \left[\frac{ebs's}{s'+s}\cosh((s'+s)eb)-\frac{\tanh ( s'eb)\tanh ( s eb)}{\tanh ( s'eb)+\tanh ( s eb)}\right ] + \frac{ebs's}{(s'+s)^2}\cosh((s'+s)eb)\\ \nonumber & + \frac{(eb)^2s's}{(s'+s)\sinh((s'+s)eb ) } -\frac{2eb\tanh ( s'eb)\tanh ( s eb)}{\sinh((s'+s)eb)(\tanh ( s'eb)+\tanh ( s eb))}\nonumber\\ \geq & \frac{eb s's}{s'+s}\left[\frac{1}{s'+s}+\frac{4s's\bar{\mu}^2}{s'+s}+m^2\right ] \left[\cosh((s'+s)eb)-1\right ] + \frac{ebs's}{(s'+s)^2}\left[\cosh((s'+s)eb)-\frac{(s'+s)eb}{\sinh((s'+s)eb)}\right]\nonumber \\ \geq & \frac{ebs's}{(s'+s)^2\sinh((s'+s)eb)}\left[\sinh ( 2(s'+s)eb)-2(s'+s)eb\right ] \geq 0.\end{aligned}\ ] ]    a crucial inequality employed above is that @xmath240 for positive @xmath241 and @xmath242 . to prove it , we notice that the function @xmath243 is monotonically increasing for @xmath244 and @xmath245 . setting @xmath246 , it follows from the inequality @xmath247 that @xmath248 repeating the same argument with @xmath249 , we end up with @xmath250 and the inequality is proved .      in this section , we shall take the nonrelativistic ( nr ) approximation of the gap equation ( [ regularizedgapeq ] ) . the conditions for the nr approximation are @xmath251 , @xmath252 with the nr binding energy @xmath253 and cyclotron frequency @xmath254 . let us break the infinite series on rhs of ( [ regularizedgapeq ] ) into two pieces , @xmath255,\ ] ] and @xmath256.\ ] ] with @xmath220 defined in ( [ thetas ] ) and @xmath257 . consequently @xmath258,\end{aligned}\ ] ] and the summation in @xmath259 can be approximated by an integral , i.e. @xmath260 using the formula @xmath261 with @xmath262 a contour on @xmath263-plane going around the positive real axis counterclockwisely , we find @xmath264 it follows from the hermite formula and stirling formula for @xmath265 that @xmath266 and @xmath267 consequently @xmath268 . \label{series}\end{aligned}\ ] ] applying the stirling formula to the gamma function on rhs of ( [ regularizedgapeq ] ) and combining with ( [ series ] ) , we end up with @xmath269.\ ] ] which , upon identifying @xmath270 with @xmath271 , is the nonrelativistic gap equation obtained in @xcite . @xmath274}\left\{\left[\frac{2}{\lambda}+4x(1-x)\bar{\mu}^2+m^2\right]\coth(\lambda eb)+\frac{eb}{\sinh^2(\lambda eb ) } \right\}\\ = & \frac{{\bar{\mu}}}{\pi^2eb}\int_0 ^ 1 dx\int_0^\infty d\xi \xi^{s+\frac{1}{2 } } x(1-x)e^{-\frac{m}{\omega_b}\xi(2x-1)^2 + 8\frac{b}{\omega_b}x(1-x)}\left\{\left[\frac{2eb}{\xi}+4x(1-x)\bar{\mu}^2+m^2\right]\coth(\xi)+\frac{eb}{\sinh^2(\xi ) } \right\},\end{aligned}\ ] ]      @xmath275}\frac{\tanh(\lambda xeb)\tanh(\lambda(1-x)eb)}{\tanh(\lambda xeb)+\tanh(\lambda(1-x)eb)}\\ \nonumber & \times\left\{\left[\frac{1}{\lambda}+4x(1-x)\bar{\mu}^2+m^2\right]\coth(\lambda eb)+\frac{2eb}{\sinh^2(\lambda eb)}\right\}\\ \nonumber = & \frac{{\bar{\mu}}}{\pi^2}\int_0 ^ 1 dx\int_0^\infty d\xi \xi^{s-\frac{1}{2 } } e^{-\frac{m}{\omega_b}\xi(2x-1)^2 + 8\frac{b}{\omega_b}x(1-x ) } \frac{\tanh(\xi x)\tanh(\xi(1-x))}{\tanh(\xi x)+\tanh(\xi(1-x))}\\ & \times\left\{\left[\frac{eb}{\xi}+4x(1-x)\bar{\mu}^2+m^2\right]\coth(\xi)+\frac{2eb}{\sinh^2(\xi)}\right\}.\end{aligned}\ ] ]    and define @xmath272 and @xmath273 as their analytic continuations at @xmath276 , i.e. @xmath277 and @xmath278 . both integrals are convergent for @xmath279 and the main contribution to the integral comes from @xmath280 . in the nonrelativistic limit , @xmath252 and @xmath251 , the integrand is peaked at @xmath281 and only the 2nd and 3rd terms inside the brackets dominate . it follows that @xmath282 upon writing @xmath283 and using the integral representation of the hurwitz zeta function , we find that @xmath284.\ ] ] under the same approximation , we have @xmath285 using the identity @xmath286,\ ] ] and integration by part , we find        d. e. kharzeev , k. landsteiner , a. schmitt and h - u yee , lect . notes phys . * 871 * , 1 ( 2013 ) and _ references therein _ ; v. miransky and i. shovkovy , phys . rep * 576 * ( 2015 ) 1 . k. klimenko , theor . phys . * 89 * , 1161 ( 1991 ) ; z. phys . c * 54 * , 323 ( 1992 ) . v. gusynin , v. miransky and i. shovkovy , phys . lett . * 73 * , 3499 ( 1994 ) ; phys . d * 52 * , 4718 ( 1995 ) ; phys . b * 349 * , 477 ( 1995 ) ; nucl . * b 462 * , 249 ( 1996 ) . g. bali , f. bruchmann , g. endrdi , z. fodor , s.d . katz , s. krieg , a. sch@xmath288fer and k.k . szab@xmath289 , j. high energy phys . 02 , ( 2012 ) 044 ; g. bali , f. bruchmann , g. endrdi , z. fodor , s.d . katz , s. krieg and a. sch@xmath288fer , phys . d * 86 * , 071502 ( 2012 ) . f.preis , a. rebhan and a. schmitt , j. high energy phys . 03 ( 2011 ) 033 . k. fukushima and j. m. pawlowski , phys . d * 86 * , 076013 ( 2012 ) . k. fukushima and y. hidaka , phys . lett . * 110 * , 031601 ( 2013 ) . t. kojo and n. su , phys . lett . * b 720 * , ( 2013 ) . f. bruckmann , g. endrodi , and t. g. kovacs , j. high energy phys . 1304 , 112 ( 2013 ) . j. chao , p. chu and m. huang , phys . rev d * 88 * , 054009 ( 2013 ) ; l. yu , h. liu and m. huang , phys . d * 90 * , 074009 ( 2014 ) ; l. yu , j. van doorsselaere and m. huang , phys . d * 91 * , 074011 ( 2015 ) . r.farias , k.gomes , g. krein and m. pinto , phys . c * 90 * , 025023 ( 2014 ) . b. feng , d.f . hou and h.c . ren , phys . d * 92 * , 065011 ( 2015 ) . e.j.ferrer , v.de la incera and x. wen , phys . d * 91 * , 054006 ( 2015 ) . n. mueller and j. m. pawlowski , phys . d * 91 * , 116010 ( 2015 ) . j. braun , w. a. mian and s. rechenberger , arxiv:1412.6025[hep - ph ] . k. a. mamo , j. high energy phys . 05 , 121 ( 2015 ) . a. flachi , k. fukushima and v. vitagliano , phys . * 114 * , 181601 ( 2015 ) . chen , k. fukushima , x.g . huang and k. mameda , arxiv:1512.08974[hep - ph ] . k. hattori , t. kojo and n. su , arxiv:1512.07361[hep - ph ] . n. mermin and h. wagner , phys . rev . lett . * 17 * , 1133 ( 1966);p . hohenberg , phys . rev . * 158 * , 383 ( 1967 ) ; s.coleman , commun . math . phys . * 31 * , 259 ( 1973 ) . b. feng , d - f hou and h - c ren , _ in preparation_. j. schwinger , phys . rev . * 82 * , 664 ( 1951 ) ; a. chodos , i. everding and d. a. owen , phys . rev . d * 42 * , 2881 ( 1990 ) . c. s de melo , m. randeria and j.r.engelbrecht , phys . lett . * 71 * , 19 ( 1985 ) . p. nozires and s schmitt - rink , j. low temp . phys . * 59 * , 195 ( 1985 ) . y. nishida and h. abuki , phys . d * 72 * , 096004 ( 2005 ) ; h. abuki , nucl . phys . * a 791 * , 117 ( 2007);h . abuki , g. baym , t. hatsuda and y. yamamoto , phys . d * 81 * , 125010 ( 2010 ) . j. deng , a. schmitt and q. wang , phys . d * 76 * , 034013 ( 2007 ) . j. wang , v. de la incera , e.j . efrain and q. wang , phys . d * 84 * , 065014 ( 2011 ) ; e.j . efrain , v. de la incera , j. keith and i. portillo , arxiv:1405.7422[hep - ph ] . g. sun , l. he and p. zhuang , phys . d * 75 * , 096004 ( 2007);l . he and p. zhuang , phys . d * 76 * , 056003 ( 2007 ) . g. endrdi , arxiv:1504.08280[hep - lat ] . ferrer , v. de la incera , phys . * 102 * , 050402 ( 2009 ) ; e.j . ferrer , v. de la incera , nucl . * b 824 * , 217 ( 2010);b . feng , e. j. ferrer and v. de la incera , nucl . * b 853* 213 ( 2011 ) ; phys . * b 706 * , 232 ( 2011 ) ; phys . d * 85 * , 103529 ( 2012 ) . a. schmitt , q. wang and d. h. rischke , phys . lett . * 91 * , 242301 ( 2003 ) ; b. feng , d. f. hou , h. c. ren and p. p. wu , phys . lett . * 105 * , 042001 ( 2010 ) .
Abstract:

Answer: the bose - einstein condensation of bound pairs made of equally and oppositely charged fermions in a magnetic field is investigated using a relativistic model . the gaussian fluctuations have been taken into account in order to study the spectrum of bound pairs in the strong coupling region . we found , in weak coupling reagion , the condensation temperature increases with an increasing magnetic field displaying the magnetic catalysis effect . in strong coupling region , the inverse magnetic catalysis appears when the magnetic field is low and is replaced by the usual magnetic catalysis effect when magnetic field is sufficiently high , in contrast to the nonrelativistic case where the inverse magnetic catalysis prevails in strong coupling region regardless of the strength of the magnetic field . the resulting response to the magnetic field is the consequence of the competition between the dimensional reduction by landau orbitals in pairing dynamics and the anisotropy of the kinetic spectrum of the bound pairs . we thus conclude that dimensional reduction dominates in weak domain and strong coupling one except the small magnetic field region , where the enhanced fluctuations dominates .


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: in recent years a considerable improvement in the study of large - scale structure has been obtained with the development of three - dimensional codes able to follow simultaneously the evolution of dark and baryonic matter . it is important to include the gas in n - body simulations for at least two reasons : first we must estimate the large amount of gas that is observed in clusters of galaxies through x - ray emission ; second , we must remove any ambiguity in the identification of galaxies . in the numerical simulations , the hydrodynamical equations have been solved either using mesh based methods @xcite or particle methods like the smoothed particle hydrodynamics @xcite . the latter , being lagrangian methods , have the advantage of covering a larger range of clustering scales , but they are not very reliable in low density regions for the calculation of thermodynamical quantities . eulerian methods , on the contrary , give a better estimate of thermodynamical quantities , but are constrained by the limited grid resolution in their capability of following the collapse of high density regions . a large dynamical range can be attained only by increasing the number of grid - points , but this is obviously limited by the available computer memory . at present the largest mesh used in eulerian calculations has @xmath4 grid - points . the analysis of the evolution at different cosmological scales is obtained by changing the physical normalization of the box size , and sometimes the results obtained at a given scale are used for setting up initial and boundary conditions at another scale ( e.g. cen & ostriker 1992 ) . since three - dimensional calculations pose severe limits on the number of grid - points used in eulerian methods , there has been in recent years an attempt to make up for the lack of resolution with the use of more accurate algorithms . one is looking for numerical methods able to well describe both low density regions and strong gradients in the fluid flow . very interesting is , for example , the possibility of introducing mesh refinement schemes in grid - based codes which , associated with shock capturing methods , would considerably improve the resolution in high density regions @xcite . the limited resolution reachable in three - dimensional codes makes sometimes difficult to separate , in the discussion of results , those effects which are only numerical from those which contains real physical information . analysis like the one of kang et al . ( 1994 ) are very useful for discriminating between different methods in three dimensions , but , because of the limitation in computer memory affecting all of the methods , a simultaneous one - dimensional analysis of the same phenomena appears quite convenient . working only in a single dimension makes possible to span a large dynamical range and to follow structure formation at different scales . approximate analytic solutions and their behaviour in highly non - linear regimes can also be better tested in one - dimensional calculations ( e.g. williams et al . recently , quilis , ibnez & saez ( 1994 ) have studied the applicability of modern high - resolution shock - capturing methods to the study of cosmological structures in presence of pressure forces . however , they do not follow at the same time the dynamics of the dark matter and so they are not able to compare the clustering properties of the two different components . in our work , instead , we study the evolution of cosmological one - dimensional perturbations when both baryonic and dark matter are considered . a detailed study of the evolution of gas and collisionless matter in a single pancake was presented by bond et al . ( 1984 ) and shapiro & struck - marcell ( 1985 ) . in both of these works one - dimensional calculations of the coupled evolution of baryonic and dark component were performed including the effects of radiative and compton cooling terms and thermal conductivity of the gas . the main difference between these two analyses is in the choice of the numerical algorithm used for solving the hydrodynamical equations which was lagrangian in one case @xcite and eulerian in the other one @xcite . a similar analysis has been recently done by thoul & weinberg ( 1994 ) , for a spherical configuration . they also have used a lagrangian approach which has the advantage of giving a higher resolution where needed without having , in one - dimension , the same problems of grid distortion present in more dimensions . all of the previous one - dimensional works following the coupled evolution of dark and baryonic matter consider isolated perturbations in the expanding universe and start their simulations when hydrodynamical effects begin to become important in the evolution . the information concerning the formation epoch and mass of these isolated structures is derived from n - body calculations assuming some model of large - scale structure formation .    in this work we do not limit our attention to isolated perturbations but we want to explore the differences in the statistical properties of the one - dimensional clustering of baryonic and dark matter ( hereafter bm and dm , respectively ) components in an expanding background . we analyse the effect of pressure forces and adiabatic heating in the dynamics of the gas and we do not include , at this stage , cooling terms in our computations . our aim is to understand first the evolution of one - dimensional structures only in the presence of `` compressible '' effects . in addition , the inclusion of cooling processes requires the specification of characteristic time and length scales in the problem and in one - dimensional cosmological calculations it is not obvious how one can choose these quantities in connection with theory and observations . therefore , we have preferred to make simulations which are scale - free ; the absence of features in the primordial power - spectrum together with the higher spatial resolution implied by the use of one - dimensional simulations would permit to better discriminate between numerical artifacts and physical effects in the results . at the end , physical information on the computed quantities can easily be extracted by fixing the normalization quantities , like the box size and the final time . the plan of the paper is as follows . in section 2 we introduce the dynamical equations both for the hydrodynamical variables and for collisionless matter . the numerical methods used for the solution of these equations are also presented together with some numerical tests . the results obtained in numerical simulations of cosmological one - dimensional structures are shown in section 3 . discussion and conclusions are drawn in section 4 . in this section we present the dynamical equations that we use for evolving the bm and dm components in the case of one - dimensional perturbations . we assume a flat universe ( @xmath5 ) in which the bm component accounts for @xmath6 of the total mass ( @xmath7 ) and the rest is in the form of dm ( @xmath8 ) . in the present analysis we neglect the effects of radiative and compton coolings and any possible external heating . the dm component is approximated as a pressureless fluid and the bm component as a perfect fluid . we introduce the following set of dimensionless variables : @xmath9 where @xmath10 is the comoving coordinate , @xmath11 is the time , @xmath12 is the cosmological expansion factor , @xmath13 is the peculiar gravitational potential , @xmath14 and @xmath15 are the matter and internal energy densities , respectively . we have three basic units of normalization : the final time @xmath16 , the comoving cell size @xmath17 and the mean baryon density @xmath18 at time @xmath16 . for a cosmological application we can take @xmath16 and @xmath18 equal to the present age and density of the universe , respectively . the evolution of the collisional component is described by the following set of hydrodynamical equations : @xmath19 = 0\ , \ ] ]    @xmath20 = -{\dot a\over a}\tilde \varrho \tilde v-{1\over   a}\tilde \varrho{\partial\tilde \phi\over\partial\tilde x}\ , \ ] ]    @xmath21 = -2{\dot a\over a}\tilde \epsilon-{1\over a } \tilde p{\partial \tilde v \over\partial \tilde x}\ , \ ] ]    where @xmath22 is the comoving pressure . we use the equation of state of an ideal gas with adiabatic index @xmath23 , so that @xmath24    the physical temperature @xmath25 is related to the previous normalized quantities by @xmath26 where @xmath27 is boltzman s constant and @xmath28 the mass of the proton .      in the dynamical equations for the collisionless particles we use as time integration variable the scale factor @xmath29 : @xmath30 @xmath31    the peculiar gravitational potential due to both the bm and dm component is computed by solving the poisson equation : @xmath32 where @xmath33 and @xmath34 is the newton s constant . the hydrodynamical equations have been integrated using the flux corrected transport ( fct ) method @xcite , a hybrid shock capturing method . in this method two difference schemes are blended together : a second order lax - wendroff scheme is used in regions of smooth flow , while a first order lax scheme is used near discontinuities . then second order accuracy is ensured everywhere except near flow jumps where the dissipation introduced by the low order scheme guarantees monotonicity in the behaviour of flow variables . the blending between the two schemes is controlled by a monotonicity constraint that leads to the sharpest possible discontinuity profiles . however , if the grid is not sufficiently fine , very strong gradients can be represented as a sequence of discontinuous jumps . this effect usually does not hinder the convergence of the method and can be reduced or completely avoided refining the grid @xcite .    for the integration of the collisionless matter we have used a particle - mesh code @xcite . we use the same mesh as in the hydrodynamical part of the code , while the number of particles is a multiple of the number of mesh points in order to have a good resolution also in low density regions . the interpolation used for computing the mass density and the forces acting on each particle is obtained by a tsc scheme ( e.g. hockney & eastwood 1981 ) , which ensures a good accuracy in the estimates of previous quantities without leading to an excessive slowing down of the computation . for the time integration we use a second order leap - frog method so that the accuracy in the n - body is comparable to that attained in the hydrodynamical part of the code . the peculiar gravitational potential is computed using the standard fast fourier transform technique . the accuracy of the code has been tested against known analytical solutions . for the hydrodynamical part we performed the shock tube test comparing the analytical and numerical solutions . we found a very good agreement in smooth regions , while the shock is numerically diffused over about five grid - points . the evolution of an initial sinusoidal perturbation has been instead used for testing the performance of the n - body code : in this case the numerical solution has been compared with the results obtained using the zeldovich approximation @xcite which provides the exact one - dimensional solution up to shell crossing . in figure 1 we compare , just before shell crossing , the resulting dm distributions for the zeldovich and n - body solutions for two different choices of the number of grid - points @xmath35 ( @xmath36 and @xmath37 ) . the number of particles is instead kept fixed and set equal to @xmath38 . we notice that already with the smallest number of grid - points the n - body solution is in very good agreement with the zeldovich result . the difference in the height of the peak in the two cases is only due to the different numerical resolution implied by the adopted grids . we used a sinusoidal perturbation also for testing the coupled evolution of dark and baryonic components and to see how far the code is able to follow the collapse of a baryonic fluctuation . in figure 2 we show the density profile of both components after the shell crossing time , obtained using a grid of @xmath39 points . we can see that fluctuations in the baryonic matter as larger as @xmath40 can be easily resolved with this resolution . some typical problems of the fct method start to appear after the formation of strong shocks which are not well resolved using the previous grid , but they do not destroy the solution and a local spatial refinement would permit to continue the evolution . since situations as this one are never found in the following cosmological simulations , we believe that a grid of @xmath39 is sufficient for our purposes . in this section we present the results of the evolution of one - dimensional perturbations in an expanding background . we analyse the distribution of bm and dm , paying particular attention to the clustering properties of the two components and to the effects due to their coupling . we consider a grid of comoving length @xmath41 , whose physical size grows according to the expansion parameter @xmath29 . we subdivide our grid into @xmath39 nodes and adopt periodic boundary conditions . the number of particles used in the n - body part of the code is @xmath38 . the initial conditions for a given model are determined as follows . the dm component is perturbed according to a gaussian random density field characterized by a power - law spectrum of the general form @xmath42 where @xmath43 is a normalization constant and @xmath1 the spectral index . the short wavelength cut - off at the scale @xmath44 grid - points ensures that the results are not affected by the sampling of modes whose size is close to that of the resolution of the simulation . the normalization constant @xmath43 is fixed requiring that the one - dimensional variance @xmath45 at the final time @xmath16 ( corresponding to the expansion factor @xmath46 ) for linearly evolved perturbations is equal to unity . in computing the variance we assume a filtering radius @xmath47 so that @xmath48 the function @xmath49 is a one - dimensional top - hat filter : @xmath50    initially the dm particles are evolved according to the zeldovich ( 1970 ) algorithm which is known to give the exact solution up to the shell crossing . we start our numerical computation when one of the following requirements is satisfied : i ) the largest density fluctuation is equal to unity ; ii ) the first shell crossing has occurred : this typically corresponds to epochs with @xmath51 for all models . the numerical evolution of bm starts when @xmath52 , i.e. just after the recombination epoch . before this time we assume that the distribution of bm is constant , as it would be in comparison with dm distribution , because of the coupling between radiation and baryonic matter before hydrogen recombination . the internal energy is fixed assuming that the cooling of the gas is due only to adiabatic expansion between the epoch of recombination and the initial time of the simulations . in the subsequent evolution the perturbation in the baryonic matter is induced by gravitational coupling with the dm component . as already pointed out , our simulations are completely scale - free : the final amplitude of the fluctuations depend only on the scale @xmath53 relative to which the initial normalization is fixed . the final thermodynamical quantities are completely determined once the amplitude of initial perturbations and the initial value of @xmath54 are given . for example , if we think of @xmath53 as corresponding of @xmath55 ( @xmath56 is the hubble constant in units of 100 km s@xmath57 mpc@xmath57 ) , where observationally the variance is found to be unity , the whole box would correspond to @xmath58 mpc ( we use @xmath59 ) and then we would be looking the evolution at very large scales . in order to study the evolution on smaller scales but fixing the normalization of the initial spectrum at the same physical scale , we should use a larger value for @xmath53 in units of @xmath41 and , consequently , we should start with larger density fluctuations . these simulations would produce larger final temperatures on these smaller scales as we expect . we consider three different models with primordial spectral index @xmath60 : in three dimensions , these would correspond to @xmath61 ( the so - called harrison - zeldovich spectrum ) respectively , covering in this way the range of values usually adopted for cosmological models . for each model we run three simulations with different realizations of the initial conditions in order to obtain more accurate estimates of the relevant quantities . here we consider various statistical tests , such as the distribution of density perturbations , peaks and thermodynamical quantities , their correlations and the density power - spectrum . the results referring to the dm component have been previously smoothed by a gaussian filter with a radius of ten grid - points , roughly corresponding to twice the numerical spreading of the shock front in the hydrodynamical calculations .      in figures 3a and 3b we show the behaviour of @xmath62 and @xmath63 in a realization for each model at two different epochs : @xmath64 and the final one @xmath46 , respectively . the results refer to about a tenth of the whole grid . different runs of the same model show a very similar qualitative behaviour . the density fluctuations are computed with respect to the mean value of the corresponding component . as expected , increasing the spectral index @xmath1 , we have less power on large scales and , consequently , peaks in density contrast are more frequent , but they have a smaller density contrast , while underdense regions are less extended . we notice that only for @xmath2 peaks in bm are higher than those in dm and their contribution to the local gravitational field starts to be comparable to that of dm .    for the @xmath65 case there is a correspondence between the distribution of dm and bm although there are always substructures in the collisionless component which are not present in the other one ( see also the following discussion on relative bias , power - spectrum and number of peaks ) . this behaviour is enhanced as we decrease the spectral index : bm appears more and more clumped while dm stays more spread with small substructures . this is particularly evident in the structures appearing in the model with @xmath2 , where collisionless matter is spread over a region which is several times larger than that occupied by the peak in bm . it is interesting to notice also that there is not an exact correspondence between the positions of the largest bm and dm peaks .    comparing the figures at two different times we see the tendency to the merging of substructures which produces higher peaks and larger voids . this behaviour is again more evident in simulations with more power on large scales . for example , in the panels referring to the @xmath2 simulation , it is possible to note that the two central peaks in the bm component at @xmath64 quickly merge forming a unique high - density region at @xmath46 . this phenomenon is less pronounced for dm , which at @xmath46 continues to be divided into small substructures . a quantitative estimate of the spatial correspondence between the bm and dm density distributions can be obtained studying the cross - correlation coefficient ( see e.g. coles , melott & shandarin 1993 ) @xmath66 where @xmath62 and @xmath63 are respectively the density contrasts in the bm and dm components , smoothed with a gaussian window of radius @xmath67 and @xmath68 are the corresponding r.m.s .. the mean is computed over all grid - points . definition ( 18 ) implies @xmath69 . the limit @xmath70 corresponds to @xmath71 , where @xmath72 is a constant : in this case there is a perfect agreement in the positions of the structures in the two components . in figure 4 , the cross - correlation coefficient @xmath73 is plotted as a function of the gaussian filtering radius @xmath67 for the three models when @xmath64 and @xmath46 . the error bars , shown for clarity only at the final time , represent the scatter r.m.s . between the three different realizations of each model . at earlier times these errors are found to be always smaller . as expected , decreasing the filtering radius and/or allowing for a longer evolution , the differences between the two components increase and , consequently , the coefficient @xmath73 decreases . this is particularly true for the model with @xmath2 , where , for small @xmath67 ( up to @xmath41/400 ) , the cross - correlation is even less than 0.4 , denoting the small correspondence between bm and dm distributions which was already clear in figure 3 . on the contrary , a value of @xmath73 close to unity is obtained in the simulations with @xmath65 , even when a small filtering radius is adopted . the study of high - density regions in bm is particularly interesting because they are expected to be related to the positions where galaxy formation occurs . we define as peaks all the grid - points which are local maxima of the density fluctuation field . in figure 5 we show , at two different times , the number of peaks @xmath74 having an height greater than a given threshold @xmath75 , both for bm and dm components . the comparison between the models with different spectral indices shows once again a different behaviour . only for the model with @xmath2 , the number of high peaks in bm is larger than that of dm at both considered times . the latter model is also the one which shows the largest time evolution in the number of peaks . in the other cases , @xmath74 does not changes significantly and the peak number in the bm component is always less than that in dm one .      in figure 6 we show the behaviour of the baryonic quantities @xmath76 , @xmath22 and @xmath25 for the case @xmath2 at @xmath46 ( the results refer to the same realization shown in the upper panel of figure 3b but in this case the whole grid is displayed ) . this is the model in which the highest values in thermodynamical quantities are obtained , thus some effects , which are also present in the other models , are better visible . the baryonic density and the pressure are given in normalized units , while the temperature is shown in kelvin degrees . in absence of heating from external sources and cooling processes , we notice a strong correspondence among peaks in density , pressure and temperature . however , in correspondence of the largest density fluctuations we observe a double peak on the top of the temperature profile which shows , as expected , that heating mainly occurs in the regions in which the velocity gradient is large . the temperature distribution in the different models is better illustrated by the cumulative volume filling factor and mass fraction @xmath77 as a function of the temperature @xmath25 . figure 7 shows the behaviour of these quantities at @xmath46 as obtained by summing over the contribution of all three realizations for each model . we notice that the fraction of matter at high temperature ( e.g. @xmath78 k ) increases as we decrease the spectral index : it is approximately 15% and 2% for @xmath2 and @xmath79 , respectively , while no grid - point has such a high temperature in the @xmath65 model ( in this case the highest temperature obtained is about 25 k ) . at the same time , most of the volume is at low temperatures : this is particularly evident in the case @xmath2 , where , due to the dominance of empty regions , more than 90% of the volume has @xmath80 k. on the contrary , the two cumulative functions are more similar in the model with @xmath65 , where structures are smaller and more uniformly distributed . these results suggest that , in absence of other dissipative or heating effects ( viscosity , thermoconduction , starbursts , etc . ) , it is difficult to obtain from one - dimensional cosmological large - scale perturbations temperatures as high as those seen in x - ray observations . figure 8 shows contour plots for the number of zones characterized by a given temperature and baryonic density . the additional information which we can extract from these graphs is the correlation between these two quantities in all the three models . this correlation is particularly strong in the case @xmath2 and refers to the large empty regions characterizing this model . on the contrary , more scattered distributions are presented in the models with @xmath79 and @xmath65 .      in order to study the time evolution of clustering at different scales , we have computed the power - spectrum @xmath81 . if the density fluctuations remain linear at large scales , @xmath81 is expected to grow as @xmath82 . in figure 9 we show the behaviour of the power - spectrum both for dm and bm for the three models at different times . the units of the wavenumbers @xmath83 are such that @xmath84 corresponds to the fundamental wavelength of the computational grid . considering the models with @xmath65 and @xmath79 , we see that at small wavenumbers ( i.e. at large scales ) the growth of the perturbations in the dm component is in good agreement with linear theory . the situation is different for the model with @xmath2 where the growth of perturbations is non - linear even at these scales . at larger wavenumbers ( i.e. at small scales ) deviations from linear theory are found , as expected . in particular , we notice that there is a faster evolution for the models with @xmath2 and @xmath79 : the non - linear growth produces the coupling of high- and low-@xmath83 modes . in the model with @xmath65 , where initially the large - scale power is small , this effect does not appear . then , our simulations , even if starting with a power - law spectrum , seems to exclude the possibility of having self - similar evolution when large scale power is present . only the model with @xmath65 presents a self - similar evolution in the dm at least for @xmath85 . the behaviour of the bm is very different . at the beginning , due to the assumption of uniform distribution on all scales , the power - spectrum is zero . at early times all models show a very slow growth produced by the coupling with the dm fluctuations ( the curves for the model @xmath2 are not visible because @xmath81 is too low ) . only at late epochs ( @xmath86 ) , the fluctuations in the bm grow considerably and quickly reach those in the dm . at the final time @xmath46 the power - spectra in the two components are different at small scales , where the growth of gas perturbations is delayed by the pressure : this is particularly true in the models with @xmath2 and @xmath79 , for which higher values of the pressure are reached .    an interesting way to measure the relative growth of the perturbations is the bias factor @xmath87 , defined as @xmath88 where @xmath89 is the variance of the density fluctuation field in the @xmath90-component , smoothed using a gaussian window of radius @xmath67 . figure 10 shows the behaviour of @xmath87 for the three models at @xmath64 and @xmath46 . the scatter r.m.s . between the different realizations of each model is also plotted only at the final time . the models with @xmath79 and @xmath65 behave in a similar way : the dm appears more clustered than the bm at small scales , while at larger scales they tend to have the same clustering properties . as expected , the pressure in the bm prevents shell - crossing and a collapse comparable to that occurred in the dm . at larger scales where dissipation is negligible , the two distributions are similar . in the case @xmath2 instead , the bm at the final time @xmath46 has a larger variance than the dm component at all scales ; this is in agreement with the density distribution shown in figure 3b . the motivation for the different behaviour of this model is the large peculiar velocities of the dm component . in fact , the absence of dissipative effects delays the formation of more compact structures . in the bm , instead , higher temperatures and pressures are reached , but they are not sufficient to overcome the gravitational force and avoid the collapse . in addition the conversion of the kinetic into internal energy allows the formation of more compact baryonic structures .    from figure 10 it is possible to extract also information about the time evolution of the bias factor . again we can note a different behaviour between the models . in the case @xmath2 , the evolution of @xmath87 from @xmath64 to @xmath46 is rapid at all scales : @xmath87 changes by a factor @xmath91 . on the contrary , the increase of @xmath87 in the same time interval is smaller for @xmath79 and @xmath65 ( approximately 15% and 5% , respectively ) but the final value is always below unity . in this paper we have studied the evolution of one - dimensional perturbations in a medium composed by bm and dm . initial gaussian perturbations with a power - law spectrum are assumed in the dominant dm component and then the perturbations in the bm component , initially uniformly distributed , are induced through gravitational coupling . we observe that bm does not fall immediately into the potential wells created by the distribution of the collisionless component , but when this happens , the amplitude of density perturbations in the gas becomes quickly comparable to that of the dm . at large scales the final distribution of the two components , as shown by the cross - correlation coefficient and by the bias factor , is very similar . the main differences are present at small scales where the dm component presents a lumpiness which is not found in the bm . a quite important feature shown by all simulations is the non - exact correspondence between largest dm and bm peaks .    in our simulations peaks in the dm are always higher than those in bm except for @xmath2 , i.e. for the spectrum with more power on large scales . in this case the variance in baryon density distribution is larger than that of dm at all scales . moreover , in the highest baryon peaks the self - gravity of baryons is comparable to the gravitational field due to the collisionless component and this can have an important effect on the evolution of the latter . in general , many peaks in the dm component are associated to a single peak in the bm . these are usually embedded in the corresponding baryonic peak which stays more spread in space because of the effect of pressure working against the collapse . for the highest bm peaks present in the case @xmath2 the situation is reversed : the pressure , although higher than in the previous cases , is not able to oppose sensibly the infall of bm and we end up with a distribution in the dm which is more spread than that in the bm . this behaviour is confirmed by the value of the bias which , for @xmath2 is larger than one at all scales . our analysis is completely scale - free and the final distribution of the gas temperature depends only on the assumed power spectrum . the highest values of the temperature are obtained for the case @xmath2 and they correspond to the high density regions . as the spectral index increases , large - scale motion is reduced and this decreases the height of the baryonic peaks and therefore also the final heating of the gas is smaller . dissipative processes present in the collisional component allow the transformation of kinetic into internal energy and this makes possible the formation of narrow baryonic peaks even in the absence of cooling processes not included in the present work . in our simulations pressure forces have never been able to equal gravitational forces and so the motion of infall is expected to continue even in presence of high temperatures . the inclusion of cooling terms would amplify this effect and possibly cause a fragmentation of baryonic peaks . work is presently in progress for testing the effect of cooling and heating processes in the evolution of one - dimensional cosmological perturbations . this , however , will require the specification of the physical scales involved in the problem and therefore the scale - free character of the present study will then be lost . we thank francesco lucchin and sabino matarrese for discussions . this work was partially supported by italian murst . anninos p. , norman m.l . , clarke d.a . , 1994 , apj , 436 , 11    bond j.r . , centrella j. , szalay a.s . , wilson j.r . , 1984 , apj , 210 , 515    boris j.p . , book d.l . , 1973 , j. comput . 11 , 38    bryan g.l . , norman m.l . , stone j.m . , , ostriker j.p . , 1994 , comput . comm . , submitted    cen r. , 1992 , apjs , 78 , 341    cen r.y . , ostriker j.p . , 1992 , apj , 393 , 22    chiang w. , ryu d. , vishniac e.t . , 1989 , apj , 339 , 603    coles p. , melott a.l . , shandarin s.f . , 1993 , mnras , 260 , 765    evrard a.e . , 1990 , apj , 363 , 349    hernquist l. , katz n.s . , 1989 , apjs , 64 , 715    hockney r.w . , eastwood j.w . , 1981 , computer simulations using particles . mcgraw - hill , new york , ny    kang h. , ostriker j.p . , cen r. , ryu d. , hernquist l. , evrard a.e . , bryan g.l . , norman m.l . , 1994 , apj , 430 , 83    katz n. , gunn j.e . , 1991 , apj , 377 , 365    quilis v. , ibnez j.m . , saez d. , 1994 , a&a , 286 , 1    ryu d. , ostriker j.p . , kang h. , cen r. , 1993 , apj , 414 , 1    shapiro p.r . , struck - marcell c. , 1985 , apj , 57 , 205    thoul a.a . , weinberg d.h . , 1994 , preprint , astro - ph/9410009    williams b.g . , heavens a.f . , peacock j.a . , shandarin s.f . , 1991 , mnras , 250 , 458    woodward p. , colella p. , 1984 , j. comput . phys . , 54 , 115    zalesak s.t . , 1979 , j. comput . 31 , 335    zeldovich ya . b. , 1970 , a&a , 5 , 84 * figure 1 . * sinusoidal test for the single dark matter component : distribution of density fluctuations just before shell crossing . comparison between the results obtained from n - body ( solid line ) and zeldovich approximation ( open squares ) , when two different numbers of grid - points @xmath35 are adopted : @xmath36 ( left ) and @xmath37 ( right ) . * figure 2 . * sinusoidal test for the coupled evolution of baryonic and dark matter components : distribution of density fluctuations after shell crossing . the dotted and solid lines refer to the baryonic and dark matter components , respectively .    * figure 3a . * the density fluctuation distribution at @xmath64 for the three different models : @xmath2 ( top panel ) , @xmath79 ( central panel ) and @xmath65 ( bottom panel ) . the dotted and solid lines refer to the baryonic and dark matter components , respectively . only about a tenth of the whole grid is displayed . * figure 4 . * the cross - correlation coefficient @xmath73 as a function of the gaussian filtering radius @xmath67 ( in units of the whole grid @xmath41 ) at the time @xmath64 ( dotted line ) and @xmath46 ( solid line ) for the different models : @xmath2 ( left ) , @xmath79 ( centre ) and @xmath65 ( right ) . error bars ( displayed for clarity only at @xmath46 ) represent the r.m.s . scatter between the different simulations . * figure 5 . * the number of peaks @xmath74 in dark ( solid line ) and baryonic ( dotted line ) components as a function of their height @xmath75 when @xmath64 ( left column ) and @xmath46 ( right column ) for the different models : @xmath2 ( top ) , @xmath79 ( centre ) and @xmath65 ( down ) . * figure 6 . * the distribution of baryonic density fluctuation @xmath92 , pressure @xmath22 and temperature @xmath25 for the case @xmath2 at @xmath46 . the whole grid is displayed . pressure is in arbitrary units while temperature is in kelvin degrees . * figure 7 . * the cumulative volume filling factor ( solid line ) and mass fraction ( dotted line ) as a function of the temperature @xmath25 at the final time @xmath46 for the different models : @xmath2 ( left ) , @xmath79 ( centre ) and @xmath65 ( right ) . * figure 8 . * contour plots for the number of grid - points characterized by a given temperature @xmath25 and baryonic density @xmath93 at the time @xmath46 for the different models : @xmath2 ( left ) , @xmath79 ( centre ) and @xmath65 ( right ) . the contour levels are defined as follows : @xmath94 , where @xmath95 is a positive integer ; the outermost contour corresponds to @xmath96 and contours inside it have gradually increasing @xmath95 .    * figure 9 . * the time evolution of the power - spectrum @xmath81 for baryonic ( left column ) and dark matter ( right column ) component for the different models : @xmath2 ( top ) , @xmath79 ( centre ) and @xmath65 ( down ) . the curves refer to different epochs : @xmath97 ( dotted lines from down to top ) and @xmath46 ( solid line ) . * figure 10 . * the relative bias @xmath98 as a function of the filtering scale @xmath67 at @xmath64 ( dotted line ) and @xmath46 ( solid line ) for the different models : @xmath2 ( left ) , @xmath79 ( centre ) and @xmath65 ( right ) . error bars ( displayed for clarity only at @xmath46 ) represent the r.m.s . scatter between the different simulations .
Abstract:

Answer: we present the results of a series of one - dimensional n - body and hydrodynamical simulations which have been used for testing the different clustering properties of baryonic and dark matter in an expanding background . initial gaussian random density perturbations with a power - law spectrum @xmath0 are assumed . we analyse the distribution of density fluctuations and thermodynamical quantities for different spectral indices @xmath1 and discuss the statistical properties of clustering in the corresponding simulations . at large scales the final distribution of the two components is very similar while at small scales the dark matter presents a lumpiness which is not found in the baryonic matter . the amplitude of density fluctuations in each component depends on the spectral index @xmath1 and only for @xmath2 the amplitude of baryonic density fluctuations is larger than that in the dark component . this result is also confirmed by the behaviour of the bias factor , defined as the ratio between the r.m.s of baryonic and dark matter fluctuations at different scales : while for @xmath3 it is always less than unity except at very large scales where it tends to one , for @xmath2 it is above 1.4 at all scales . all simulations show also that there is not an exact correspondence between the positions of largest peaks in dark and baryonic components , as confirmed by a cross - correlation analysis . the final temperatures depend on the initial spectral index : the highest values are obtained for @xmath2 and are in proximity of high density regions . galaxies : clustering  galaxies : formation  large - scale structure of universe