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: many astrophysical observations provide evidences for the existence of non - baryonic dark matter from the galactic scale to the cosmological scale  @xcite . its existence in our universe has been revealed only by its gravitational effects , while its mass and interactions to the standard model ( sm ) particles are still unclear . among many scenarios that had been proposed to account for non - baryonic dark matter , the weakly interacting massive particle ( wimp ) scenario , where the particle with weak scale mass and interactions could be thermally frozen out and its abundance would be the observed dark matter relic density  @xcite , is one of the most attractive scenarios because it may be produced and detected at the current or future colliders . actually many works in the literature have dealt with this issue  @xcite , but they have mainly focused on the large hadron colliders ( lhc ) .    the dark matter search at colliders is based on the production of sm particle(s ) with missing energies . the typical signals are mono - jet  @xcite , mono - photon  @xcite , mono-@xmath3  @xcite or mono-@xmath4  @xcite with missing energies . in principle , the particles which are produced together with the dark matter particle(s ) could be any sm fundamental particles or their composite states . for example , the dark matter production associated with a heavy - quark - and - heavy - anti - quark ( @xmath5 ) pair may also play an important role in searching for dark matter . it is known that the dark matter pair production associated with the @xmath5 pair could be more effective in the case that the dark matter interaction with the sm particles depends on the mass of the relevant sm quark  @xcite .    in this paper , we propose to use the process of dark matter production associated with a heavy quarkonium to search for dark matter at colliders . a heavy quarkonium is a meson which is a bound state of a @xmath5 pair . naively speaking , the heavy - quark - pair - associated production of dark matter would be more effective for the dark matter detection than the light - quark - pair - associated production if the effective coupling of dark matter to the quark pair is proportional to the quark mass . in a similar way , the dark matter production associated with a heavy quarkonium would play a more essential role in the probe of such operators than dark matter production associated with a light meson . the dark matter production associated with the sm particle(s ) can be investigated in any collider experiments . in this work , we focus on the dark matter production associated with a heavy quarkonium at @xmath0 factories , whose center - of - mass ( cm ) energy is @xmath6 gev . however , we note that this search will easily be applied to the international linear collider ( ilc ) and large hadron collider ( lhc )  @xcite .    there is no evidence for dark matter from the collider experiments so far . however , there are some reports for dark matter candidate signals from direct detection of dark matter , which is carried out in underground experiments , and also from indirect detection of dark matter from astrophysical observations , in particular , in the dark matter mass region of @xmath7 gev . on the other hand , the dark matter mass region below a few gev has not been well investigated , especially , in the direct detection experiments of dark matter because the very light dark matter can not hit a nucleon with significant recoil . one of the strongest bounds on the dark matter signals in the indirect detection of dark matter comes from the @xmath8-ray observation of dwarf spheroidal galaxies by the fermi - lat satellite  @xcite . however , the constraints become weakened below the dark matter mass region less than about 4 gev because of the uncertainty of hadronization of final state particles in the dark matter annihilation and the limit of the photon energy threshold of the detector in the fermi gamma - ray space telescope  @xcite .    since such light dark matter mass region could be investigated at colliders  @xcite , the collider detection of dark matter would be complementary to the direct detection of dark matter and the indirect detection of dark matter .    at @xmath0 factories , dark matter can not only be produced directly , but also produced from the decays of a heavy quarkonium . for example , the invisible decay of @xmath9 can give constraints on the properties of light dark matter through the decay of @xmath10 followed by @xmath11 nothing at @xmath0 factories  @xcite . also , the charm factories can contribute to search for dark matter by investigating the invisible decay of a charmonium , @xmath1  @xcite . it is found that the @xmath9 and @xmath1 decays at @xmath0 and charm factories are more suitable to probe light dark matter interactions than the mono - jet search at high energy hadron colliders when the mass of mediator is not large  @xcite . the dark matter production associated with a heavy quarkonium at @xmath0 factories would also be complementary to the invisible decay of a heavy quarkonium as well as the mono - jet searches at high energy colliders . this paper is organized as follows . in sec . [ model ] , we describe the relevant effective operators for the dark matter interactions with the sm quarks . as a candidate for dark matter , we take into account two cases : dirac fermionic dark matter and real scalar dark matter . we adopt the effective field theory ( eft ) approach for the interaction of dark matter with the sm quarks . in sec . [ calc ] , we calculate the cross sections for the dark matter production at @xmath0 factories as well as for the sm backgrounds . we also interpret the constraints on the properties of the dark matter interaction to the upper limit on the dark matter - nucleon scattering cross section . finally , we summarize our results in sec . [ summary ] . for a model - independent search for dark matter , we assume that @xmath12 , which stands for either a dirac fermion or real scalar , is the only component of dark matter and singlet under the sm gauge group . the extension to the majorana fermion or complex scalar dark matter would be straightforward and not so much different from the dirac fermion or real scalar studied here . we adopt the effective field theory approach , where the mediator connecting the dark sector to the sm particles is heavy enough to be integrated out  @xcite . after symmetry breaking , the dark matter interactions to the sm fields can be described by higher - dimensional operators . since we are interested in the dark matter production associated with a heavy quarkonium , the relevant effective operators must contain the heavy quark fields . the dimension-6 operators in the dirac fermion case which we will consider in this paper are  @xcite @xmath13 where @xmath14 is the heavy - quark mass and @xmath15 represents the energy scale of new physics . effectively @xmath16 , where @xmath17 is the mass of the mediator and @xmath18 are the coupling of the mediator to the dark sector and sm sector , respectively . for the ( pseudo)scalar interaction operators @xmath19 , an additional helicity suppression factor is taken into account .    in the real scalar case , the relevant dimension-6 operators are given by  @xcite @xmath20 there are more dimension-6 operators which would be necessary for complete study of the dark matter production . for example , one can consider the operators which include the gluonic field strengths ( @xmath21 ) . they may contribute to the production of dark matter particles associated with a heavy quarkonium at the one - loop level so that the dark matter production process would be suppressed . for the probe of the operators which include the electromagnetic field strengths ( @xmath22 ) or a lepton pair , the dark matter production associated with a photon would be more efficient at the lepton colliders like @xmath0 factories or ilc . it is well known that the eft approach is valid when the particle mediating the interaction is heavier than the typical energy scale of the process , which is the maximum momentum transfer in the process . at high - energy hadron collider where the energy scale of the fundamental process is not fixed , the eft approach may break down unless the mediator is not too heavy compared to the cm energy of the collider . then , the eft approach may over- or under - estimate the result in a uv - complete theory  @xcite . it is worthwhile to note that @xmath0 factories are relatively free from this issue because the cm energy is fixed as well as not so high .    some operators in eqs . ( [ effopfermion ] ) and ( [ effopscalar ] ) do not respect the sm gauge group symmetry  @xcite . in particular , the scalar operators @xmath23 require a @xmath24- or @xmath25 channel scalar exchange for uv completion . in ref . @xcite , it is shown that in the @xmath24-channel scalar - exchange model , the eft approach prediction is completely different from that of the uv complete model . therefore , one should be more cautious about the interpretation of the bounds on such operators obtained in the eft approach . the production of a heavy quarkonium is described by non - relativistic qcd ( nrqcd ) , which is an effective field theory of qcd  @xcite . as usual , the factorization between short - distance and long - distance physics is assumed . while the factorization between high - energy and low - energy physics for the inclusive heavy - quarkonium production has not been proved , the factorization for its exclusive production at @xmath0 factories was given in  @xcite . the schematic expression for the production cross section is given by @xmath26 , where @xmath27 is the long - distance matrix element ( ldme ) , which is the probability of the evolution of a heavy - quark pair into a heavy quarkonium @xmath28 , and @xmath29 is the corresponding short - distance coefficient , which is responsible for the production of a heavy - quark and anti - quark pair . the coefficient @xmath29 is obtained by integrating the squared amplitude over the phase space after averaging the spin states of initial particles and summing the spin states of final particles . the amplitude for the production of @xmath28 can be obtained by projecting the production amplitude of a @xmath5 pair onto @xmath28 via the spin and color projection operators of @xmath28 ( generically denoted by @xmath30 here ) : @xmath31,\ ] ] with @xmath32 being a matrix that acts on spinors of the @xmath5 pair in the amplitude of the @xmath5 pair production  @xcite .    in nrqcd , a heavy - quark and anti - quark pair can be produced in a color - singlet state or color - octet state . the color - octet state can evolve into a color - singlet state by emitting or absorbing soft gluons . typically , the color - octet ldme is suppressed by @xmath33 or higher in the heavy - quark velocity @xmath34 compared to the color - singlet ldme . therefore , unless there exists a kinematical enhancement factor in the short - distance coefficient for the color - octet amplitude , the color - singlet amplitude dominates . in the @xmath0 factories , the initial colliding beams are the electron and positron , which are neutral under @xmath35 gauge symmetry therefore , the production of the heavy quarkonium in the color - singlet state can occur only when the other particles are color - neutral or in the color - singlet state . however , the production of a heavy quarkonium in a color - octet state must accompany with at least a color - octet particle , for instance , an additional gluon or a pair of quark and anti - quark in the color - octet state . in the dark matter production or off - shell @xmath4 boson production at @xmath0 factories , the color - singlet heavy quarkonium production is dominant . the color - octet heavy quarkonium production has a wavefunction suppression factor of @xmath36 or higher . furthermore , the color - octet process has an additional suppression factor due to the additional color - octet gluon or a pair of quark and anti - quark . therefore , for the dark matter production associated with a heavy quarkonium at @xmath0 factories , it is sufficient to take into account only the color - singlet contribution . the color - singlet ldme is determined by the electromagnetic decay of the heavy quarkonium . explicitly , for charm quark mass @xmath37 gev and to leading - order in @xmath34 , we take the color - singlet ldmes to be  @xcite @xmath38    first we calculate the sm backgrounds . in the sm , the feynman diagrams contributing to the heavy quarkonium production with missing energy have at least one weak - boson exchange , which is connected to a pair of neutrinos . there are also diagrams with two weak - boson exchanges , but they are suppressed by @xmath39 so that we ignore those diagrams . then , we left with seven diagrams in the sm , whose representative diagrams are shown in fig . [ diagramsm ] . the irreducible background for the dark matter production associated with a heavy quarkonium is the @xmath40 process followed by @xmath41 , as depicted in the first two diagrams in fig . [ diagramsm ] . there is also a @xmath3 fusion diagram , which is required by gauge invariance , as shown in the last diagram in fig . [ diagramsm ] . other diagrams not shown are obtained by reversing the flow of fermions .    in the monojet production at hadron colliders , except the irreducible sm background @xmath42 followed by @xmath43 , there is another type of the sm background , for example , the process @xmath44 , where both @xmath45 and @xmath46 are missed . we note that this kind of background can be ignored at @xmath0 factories because the relevant process is @xmath47 , which is highly suppressed . we note that , at the @xmath48 machines , there may be continuum background coming from @xmath49 . we calculate this continuum background by choosing the invariant mass @xmath50 of the @xmath51 pair in the region @xmath52 , where @xmath53 and @xmath54 are the mass and decay width of @xmath28 respectively for @xmath55 or @xmath2 . we find that the cross section for the continuum background is less than @xmath56 ab , which can be neglected at @xmath0 factories . there may also be a continuum background coming from process like @xmath57 ( @xmath58 ) since @xmath1 may be detected experimentally via the decay @xmath59 . nevertheless the cross section from this continuum contribution is also negligible . we take the electromagnetic coupling to be @xmath60 at the scale @xmath61 and the charm quark mass to be @xmath37 gev . then , the cross sections for a heavy quarkonium production associated with a pair of neutrinos in the sm are calculated as ( 1 ab = 1 attobarn = @xmath62 @xmath63 ) @xmath64 the cross section for @xmath65 production is much smaller than that for the @xmath66 production . this is mainly because of the dominated contributions from photon - fragmentation diagrams , where a virtual photon evolves into a @xmath1  @xcite . all the diagrams in fig . [ diagramsm ] contribute to the @xmath67 production while only the first diagram contributes to the @xmath68 production . however , the first diagram does not correspond to the photon - fragmentation contribution and is suppressed by @xmath69 as compared to other photon - fragmentation diagrams . the expected event number for the @xmath66 production is about @xmath70 @xmath71 with an integrated luminosity of @xmath70 @xmath72 ab@xmath73 without any cut . since the expected event number for the @xmath65 production is entirely negligible , even one event of the @xmath2 plus missing energy measured at @xmath0 factories would imply the existence of new physics .      in this section , we calculate the cross section for the dirac fermion dark matter production associated with a heavy quarkonium . there are two feynman diagrams , one of which is shown in fig . [ diagram ] . the other diagram is obtained by reversing the flow of the charm quark . the bulb in the figure represents vertices of various effective operators . the relevant operators to the heavy quarkonium production depend on the quantum number of the heavy quarkonium , where @xmath74 for @xmath1 and @xmath75 for @xmath2 . for example , the operators @xmath76 contribute only to the @xmath77 process , while the operators @xmath78 only to the @xmath79 . we note that only the @xmath78 operators can contribute to the invisible decays of the heavy quarkonium states @xmath9 and @xmath1  @xcite . thus , the dark matter production associated with a heavy quarkonium could provide constraints on a larger set of operators , which can not be probed in the invisible decays of @xmath9 and @xmath1 . figure  [ figurexsec ] shows the cross sections ( a ) for @xmath80 and ( b ) for @xmath81 as a function of the dark matter mass @xmath82 for @xmath83 gev ( dashed line ) , @xmath84 gev ( solid line ) , and @xmath85 gev ( dotted line ) , respectively . the effective operators for the red and blue lines are @xmath86 ( @xmath87 ) and @xmath88 ( @xmath89 ) for the @xmath1 ( @xmath2 ) production , respectively . the horizontal line in fig .  [ figurexsec](a ) is the sm background in the framework of nrqcd    in fig . [ figurexsec](a ) , the cross sections for the @xmath90 production can reach about @xmath91 fb for @xmath92 gev and about @xmath93 ab for @xmath94 gev respectively in the case of the operator @xmath88 . however , for the @xmath86 operator , the cross section reaches about @xmath95 ab for @xmath83 gev , but the production is negligible for @xmath94 gev due to a suppression factor @xmath96 . the curves in the figure drop rapidly as the dark matter mass grows , indicating the kinematical limit of the dark matter mass of @xmath97 . for other operators @xmath98 , similar features can be obtained .    in fig . [ figurexsec](b ) , the cross section for the @xmath99 production can reach about @xmath100 fb for @xmath83 gev and about @xmath101 ab for @xmath94 gev in the case of the @xmath87 operator . for the @xmath89 operator , the cross section reaches about @xmath102 fb for @xmath83 gev and about @xmath103 ab for @xmath94 gev , respectively . effects from the kinematical limit of the dark matter mass of @xmath104 are also evident . for operators @xmath105 , similar features are obtained .    in fig . [ figurelam ] , we plot bounds on the scale @xmath15 ( a ) for @xmath77 and ( b ) for @xmath106 in the cases of the integrated luminosity ( @xmath107 ) of @xmath70 @xmath108 ( dashed line ) and @xmath109 @xmath108 ( solid line ) , respectively . in the @xmath90 production , we set the bound on @xmath15 to be the value where the signal - to - background ratio , which is defined by @xmath110 , is equal to 5 . in the @xmath99 production , we set the bound to be the value where the number of events is one because the sm background is negligible . in the @xmath111 production , the bound can reach about @xmath112 ( @xmath113 ) gev with @xmath114 ( @xmath109 ) @xmath108 for the @xmath88 operator . for the @xmath86 operator , the bound is much less due to the chirality suppression factor in the operator . in the @xmath99 production , the bound can reach about @xmath112 ( @xmath115 ) gev for @xmath87 operator and about @xmath116 ( @xmath117 ) gev for @xmath89 operator with @xmath114 ( @xmath109 ) @xmath108 , respectively . a couple of comments are in order here . first , as the dark matter mass approaches the kinematical limit , the bound on @xmath15 which can be achieved at @xmath0 factories becomes quite tiny as clearly shown in fig . [ figurelam ] . we note that for the region below some value of @xmath15 , our limits become untrustworthy . recall that @xmath118 . a very tiny @xmath15 can be achieved only when @xmath119 becomes large since the mediator mass @xmath17 should be larger than @xmath120 for validity of eft approach . on the other hand , @xmath121 can not be too large in order to perform perturbative calculation . for @xmath122 , the eft description breaks down in the region @xmath123 . second , we note that the number of events is not large . so the best thing one can hope for is to detect inclusive signals rather than exclusive ones . that is , the heavy quarkonium would be identified by the invariant mass distribution of all the final state particles in the decay products ( except for missing energy ) instead of using the exclusive mode like @xmath124 .      in this section , we consider the real scalar dark matter production associated with a heavy quarkonium at @xmath0 factories . the relevant feynman diagrams are the same as in the dirac dark matter case with the replacement of the dark matter fermion lines by the dark matter scalar lines . the relevant effective operators are shown in eq . ( [ effopscalar ] ) .    in fig . [ figurexsec](a ) , the cross section for the @xmath125 production in the @xmath126 operator case is depicted for @xmath83 gev ( green dashed line ) and @xmath94 gev ( green solid line ) . the cross section can reach about @xmath127 ab for @xmath83 gev , but for @xmath94 gev the cross section is below @xmath128 @xmath108 . similar feature is observed for the @xmath129 operator . we note that the cross section for the @xmath130 production vanishes for real dark matter . the bound on the scale @xmath15 for @xmath131 is shown in fig . [ figurelam](a ) for the integrated luminosity of 1 @xmath108 ( green dashed line ) and 50 @xmath108 ( green solid line ) , respectively . again the bound is set to be the signal - to - background ratio @xmath132 . the bound on @xmath15 can reach about @xmath133 @xmath134 gev for @xmath114 @xmath72 @xmath108 respectively for the operator @xmath135 . the bound on @xmath15 deduced above may be converted into the bound for the dark matter - nucleon scattering cross section for the direct detection of dark matter . we note that , however , this conversion holds with the assumption that the effective operators in eqs . ( [ effopfermion ] ) and ( [ effopscalar ] ) are either flavor blind or proportional to the quark masses . in case that dark matter has charm - philic property , the bound from the dark matter production associated with a heavy quarkonium can not be converted into the cross section for the dark matter - nucleon scattering .    for the effective operators , the dark matter - nucleon scattering cross sections are shown in refs . @xcite : @xmath136 where @xmath137 is the reduced mass of the dark matter and nucleon and @xmath138 is the mass of the nucleon . @xmath139 stands for the spin - independent ( si ) and spin - dependent ( sd ) cross section for the dark matter - nucleon scattering . the scalar form factor of the nucleon @xmath140 where @xmath141 . here , we use @xmath142 , @xmath143 , and @xmath144 , for the scalar form factors , @xmath145 , @xmath146 , for the vector form factors , and @xmath147 , @xmath148 , and @xmath149 for the axial - vector form factors , respectively  @xcite . we note that if the dark matter particles interact only with the charm quarks , the dark matter - nucleon scattering cross sections are negligible or significantly reduced .    in fig . [ figuredirect ] , we show ( a ) the spin - independent cross section and ( b ) the spin - dependent cross section versus the dark matter mass with the bounds of @xmath15 extracted for each dark matter mass from fig .  [ figurelam ] for the corresponding operators and luminosities and applied to the appropriate formulas given in eqs . ( [ dmnucleona ] ) to ( [ dmnucleond ] ) . the black and gray lines at the lower right corners are 90% c.l . exclusion limits of cmdslite  @xcite and supercdms  @xcite , respectively . also at the lower right corners of these two plots in fig . [ figuredirect ] , the orange dotted ( solid ) lines are 90% c.l . exclusion limits of cms for @xmath150 ( @xmath151 ) from 1 to 10 gev  @xcite , while the cyan dotted ( solid ) lines are for @xmath150 ( @xmath151 ) from atlas and the purple dotted ( solid ) lines are for @xmath152 ( @xmath153 ) from atlas , respectively  @xcite . since the bounds from cms and atlas are obtained for the effective operators ( [ effopfermion ] ) , they strongly depend on the uv completion of the operators if other particles in the uv completion model are not heavy enough to be integrated out at the lhc energy scale . the red and blue lines correspond to the @xmath154 and @xmath155 operators in the @xmath90 ( @xmath156 ) production , respectively . the green line corresponds to the @xmath126 operator in the @xmath157 production . the dashed and solid lines represent the cross section corresponding to @xmath114 @xmath108 and @xmath109 @xmath108 at @xmath0 factories , respectively . the bound on the dark matter - nucleon scattering cross section could reach about @xmath158 @xmath63 in the spin - independent case and @xmath62 @xmath63 in the spin - dependent case for the very light @xmath82 . for @xmath159 and @xmath160 @xmath108 at @xmath0 factories , the bound could be about @xmath161 @xmath63 in the spin - independent case and about @xmath162 @xmath63 in the spin - dependent case . the bounds obtained from the mono - jet search at the lhc might be stronger than those from @xmath0 factories . however , the bounds from the lhc are restricted to the large mass region of the mediator , but those from @xmath0 factories would be applied to broader range of the mediator mass ( @xmath163 ) .    as shown in fig . [ figuredirect ] , the lhc bounds are more efficient than the bounds obtained from the dark matter production associated with a heavy quarkonium at @xmath0 factories by a few orders of magnitude for @xmath164 gev . however , this search is meaningful to provide complementary search in the lepton colliders for the universal couplings of dark matter . if the dark matter particle dominantly couples to the charm quark , the lhc bound would become weaker for @xmath165 since the current bounds arise from the light quark interactions in the @xmath51 or @xmath166 ( @xmath167 ) collisions . furthermore , the bounds from the direct detection could be negligible since they do not give any constraints on such operators at leading order . for the scalar operators @xmath168 , the lhc bounds are still powerful , but the direct detection bounds become weaker by @xmath169 . we note that one should be careful for the comparison of the results of dark matter search at colliders with those from direct detection . as shown in ref . @xcite , the eft approach may be inadequate for the hadron collider dark matter search because the typical energy scale in the process is not fixed . in addition , some of the effective operators in eqs . ( [ effopfermion ] ) and ( [ effopscalar ] ) do not respect the sm gauge symmetry , in particular , @xmath170 symmetry . if one imposes the sm gauge symmetry for the effective operators , the uv complete model , which is responsible for generating these operators , may not approach the eft even in the limit of infinite mass of the mediator . this is true for the scalar operators in eq . ( [ effopfermion ] )  @xcite . the contact interaction for the scalar@xmath171scalar operator can not be realized at the energy scale of the lhc because the sm higgs boson which is needed to impose the sm @xmath170 gauge symmetry is still a propagating mode . in this case , the comparison in fig . [ figuredirect ] should be re - interpreted with the analysis in the full theory . however , at @xmath0 factories , one can obtain the contact interaction by integrating out both the heavy degree of freedom of the mediator and sm higgs boson . the effective scale @xmath15 will depend on the masses of the mediator as well as the sm higgs boson . the effective operators thus obtained would respect just the gauge symmetries of qed and qcd , which is necessary . in this paper , we investigate light dark matter production associated with a heavy quarkonium at @xmath0 factories . for the interaction of dark matter to the sm sector , the effective field theory approach is adopted , but it would be straightforward to extend this study to a realistic model with uv completion . so far the dark matter search at colliders has been focused on hadron colliders . in this paper , we showed that the @xmath0 factories may play a role in searching for dark matter . especially , for the light mediator mass , the @xmath0 factories would be more effective . we took into account the @xmath90 and @xmath99 production , but this would easily be extended to the charmonium in the higher resonances like @xmath172 , @xmath173 , and @xmath174 , as well as the bottomonium production . combining all the results , one may obtain stronger bounds on the effective energy scale @xmath15 , which together with the dark matter mass determines the size of the dark matter - nucleon scattering cross section for direct detection . the eft approach holds only when the mediator mass is larger than the typical energy scale of the relevant processes . therefore the dark matter search via the eft approach would be reliable only for @xmath175 . for @xmath176 , one must consider a uv complete model or another eft model including the particles which can not be integrated out at the energy scale . another important point for the eft approach is to impose the full sm gauge symmetry  @xcite . for example , the @xmath170 symmetry may require the sm higgs boson exchange between the dark sector and sm sector for the scalar@xmath171scalar interaction , which can not be integrated out at the energy of the lhc . searches for dark matter at the @xmath0 factories are relatively free from such complications .    the search for light dark matter via its production associated with a heavy quarkonium at @xmath0 factories would provide both the alternative and complementary ways to the hadron colliders and the heavy quarkonium invisible decays . g.  jungman , m.  kamionkowski and k.  griest , phys . * 267 * ( 1996 ) 195 [ hep - ph/9506380 ] . g.  bertone , d.  hooper and j.  silk , phys . * 405 * ( 2005 ) 279 [ hep - ph/0404175 ] . j.  einasto , [ arxiv:0901.0632 [ astro-ph.co ] ] . a.  de simone , g.  f.  giudice and a.  strumia , jhep * 1406 * ( 2014 ) 081 [ arxiv:1402.6287 [ hep - ph ] ] and references therein . y.  bai , p.  j.  fox and r.  harnik , jhep * 1012 * ( 2010 ) 048 [ arxiv:1005.3797 [ hep - ph ] ] . j.  goodman , m.  ibe , a.  rajaraman , w.  shepherd , t.  m.  p.  tait and h.  b.  yu , phys . d * 82 * , 116010 ( 2010 ) [ arxiv:1008.1783 [ hep - ph ] ] . p.  j.  fox , r.  harnik , j.  kopp and y.  tsai , phys . d * 85 * ( 2012 ) 056011 [ arxiv:1109.4398 [ hep - ph ] ] . y.  bai and t.  m.  p.  tait , phys . b * 723 * ( 2013 ) 384 [ arxiv:1208.4361 [ hep - ph ] ] . l.  m.  carpenter , a.  nelson , c.  shimmin , t.  m.  p.  tait and d.  whiteson , phys . d * 87 * ( 2013 ) 7 , 074005 [ arxiv:1212.3352 ] . t.  lin , e.  w.  kolb and l.  t.  wang , phys .  rev . d * 88 * ( 2013 ) 6 , 063510 [ arxiv:1303.6638 [ hep - ph ] ] . wang , c.  yu , and t.c .  yuan , work in progress . m.  ackermann _ et al . _ [ fermi - lat collaboration ] , phys .  rev . d * 89 * ( 2014 ) 042001 [ arxiv:1310.0828 [ astro-ph.he ] ] . n.  fernandez , j.  kumar , i.  seong and p.  stengel , phys . d * 90 * , no . 1 , 015029 ( 2014 ) [ arxiv:1404.6599 [ hep - ph ] ] ; n.  fernandez , i.  seong and p.  stengel , arxiv:1511.03728 [ hep - ph ] . e.  izaguirre , g.  krnjaic , p.  schuster and n.  toro , phys . d * 88 * ( 2013 ) 114015 [ arxiv:1307.6554 [ hep - ph ] ] . k.  cheung , p.  y.  tseng , y.  l.  s.  tsai and t.  c.  yuan , jcap * 1205 * ( 2012 ) 001 [ arxiv:1201.3402 [ hep - ph ] ] . o.  buchmueller , m.  j.  dolan and c.  mccabe , jhep * 1401 * ( 2014 ) 025 [ arxiv:1308.6799 [ hep - ph ] ] . n.  f.  bell , y.  cai , j.  b.  dent , r.  k.  leane and t.  j.  weiler , arxiv:1503.07874 [ hep - ph ] . s.  baek , p.  ko , m.  park , w.  i.  park and c.  yu , arxiv:1506.06556 [ hep - ph ] . g.  t.  bodwin , e.  braaten and g.  p.  lepage , phys . d * 51 * , 1125 ( 1995 ) [ erratum - ibid . d * 55 * , 5853 ( 1997 ) ] [ hep - ph/9407339 ] . g.  t.  bodwin , x.  garcia i tormo and j.  lee , phys . * 101 * ( 2008 ) 102002 [ arxiv:0805.3876 [ hep - ph ] ] ; phys . d * 81 * ( 2010 ) 114014 [ arxiv:1003.0061 [ hep - ph ] ] . e.  braaten and j.  lee , phys . d * 67 * ( 2003 ) 054007 [ phys .  rev . d * 72 * ( 2005 ) 099901 ] [ hep - ph/0211085 ] . g.  t.  bodwin , j.  lee and c.  yu , phys . d * 77 * ( 2008 ) 094018 [ arxiv:0710.0995 [ hep - ph ] ] . g.  t.  bodwin , h.  s.  chung , d.  kang , j.  lee and c.  yu , phys . d * 77 * , 094017 ( 2008 ) [ arxiv:0710.0994 [ hep - ph ] ] ; h.  s.  chung , j.  lee and c.  yu , phys . b * 697 * , 48 ( 2011 ) [ arxiv:1011.1554 [ hep - ph ] ] . y.  fan , j.  lee and c.  yu , phys .  rev . d * 87 * , no . 9 , 094032 ( 2013 ) [ arxiv:1211.4111 [ hep - ph ] ] . g.  t.  bodwin , e.  braaten , j.  lee and c.  yu , phys .  rev . d * 74 * ( 2006 ) 074014 [ hep - ph/0608200 ] . g.  belanger , f.  boudjema , a.  pukhov and a.  semenov , comput . commun . * 180 * ( 2009 ) 747 [ arxiv:0803.2360 [ hep - ph ] ] . u.  haisch , f.  kahlhoefer and j.  unwin , jhep * 1307 * ( 2013 ) 125 [ arxiv:1208.4605 [ hep - ph ] ] . h.  y.  cheng and c.  w.  chiang , jhep * 1207 * ( 2012 ) 009 [ arxiv:1202.1292 [ hep - ph ] ] . a.  crivellin , m.  hoferichter and m.  procura , phys . d * 89 * ( 2014 ) 054021 [ arxiv:1312.4951 [ hep - ph ] ] ; m.  hoferichter _ et al . _ , phys . * 115 * ( 2015 ) 092301 [ arxiv:1506.04142 [ hep - ph ] ] . r.  agnese _ et al . _ [ supercdms collaboration ] , [ arxiv:1509.02448 [ astro-ph.co ] ] . r.  agnese _ et al . _ [ supercdms collaboration ] , phys . * 112 * ( 2014 ) 24 , 241302 [ arxiv:1402.7137 [ hep - ex ] ] . v.  khachatryan _ et al . _ [ cms collaboration ] , eur . j.  c * 75 * ( 2015 ) 5 , 235 [ arxiv:1408.3583 [ hep - ex ] ] . g.  aad _ et al . _ [ atlas collaboration ] , eur . j.  c * 75 * ( 2015 ) 7 , 299 [ arxiv:1502.01518 [ hep - ex ] ] .
Abstract:

Answer: we investigate light dark matter production associated with a heavy quarkonium at @xmath0 factories in a model - independent way by adopting the effective field theory approach for the interaction of dark matter with standard model particles . we consider the effective operators for the dark matter - heavy quark interaction , which are relevant to the production of dark matter associated with a heavy quarkonium . we calculate the cross sections for dark matter production associated with a @xmath1 or @xmath2 to compare with the standard model backgrounds . we set bounds on the energy scale of new physics for various effective operators and also obtain the corresponding limits for the dark matter - nucleon scattering cross sections for light dark matter with mass of the order of a few gev .


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: it has been recognized shortly after the discovery of pulsars that neutron stars have velocities much in excess of those of any other normal stellar populations in our galaxy ( gunn & ostriker 1970 ; minkowski 1970 ) . however , it is only in the last few years that a significant body of evidence has come into place to support the view that type ii supernovae are nonspherical and neutron stars receive large kick velocities at birth . lyne and lorimer ( 1994 ) analyzed pulsar velocities in light of new proper - motion measurements ( harrison et al . 1993 ) and up - to - date pulsar distance scale ( taylor & cordes 1993 ) , and concluded that pulsars were born with a mean speed of @xmath9  km  s@xmath7 , much larger than previously thought . more recent studies of pulsar velocities have adopted more sophisticated satistical methods and included better treatment of selection effects and uncertainties ( lorimer et al . 1997 ; hansen & phinney 1997 ; cordes & chernoff 1997 ) . they all yielded a result for the pulsar birth velocity in qualitative agreement ( although not in quantitative details ) with that of lyne & lorimer ( 1994 ) . in particular , cordes and chernoff ( 1997 ) found that the three - dimensional space velocities of their sample of 47 pulsars have a bimodal distribution , with characteristic speeds of @xmath10  km  s@xmath7 and @xmath11  km  s@xmath7 ( corresponding to @xmath12 and @xmath13 of the population ) . they also estimated that pulsars with velocities greater than @xmath14  km  s@xmath7 may be underrepresented owing to selection effects in pulsar surveys [ the uncertainty in the high - velocity end of the pulsar velocity distribution function has also be emphasized by hansen & phinney ( 1997 ) ] . concrete evidence for the existence of pulsars with velocities of @xmath15  km  s@xmath7 has come from the observation of the guitar nebula pulsar ( b2224 + 65 ) , which produces a bow shock when plowing through the interstellar medium ( cordes , romani & lundgren 1993 ) . in addition , studies of pulsar - supernova remnant associations have uncovered a number of pulsars having velocities much greater than @xmath14  km  s@xmath7 ( frail et al .  1994 ) , although in many cases the associations are not completely secure . compelling evidence for supernova asymmetry and pulsar kicks also comes from the detection of geodetic precession in the binary pulsar psr 1913 + 16 ( cordes et al .  1990 ; arzoumanian et al .  1996 ) , and the orbital plane precession in the psr j0045 - 7319/b star binary ( lai et al . 1995 ; kaspi et al . 1996 ) and its fast orbital decay ( which indicates retrograde rotation of the b star with respect to the orbit ; see lai 1996 ) . these results demonstrate that binary break - up ( as originally suggested by gott , gunn & ostriker 1970 ; see iben & tutukov 1996 ) can not be _ solely _ responsible for the observed pulsar velocities , and that _ natal kicks are required_. in addition , evolutionary studies of neutron star binary population imply the existence of pulsar kicks ( e.g. , deway & cordes 1987 ; fryer & kalogera 1997 ; fryer et al .  1998 ; see also brandt & podsiadlowski 1995 ) . finally , there are a large number of direct observations of nearby supernovae ( e.g. , cropper et al .  1988 ; trammell et al .  1993 ; mccray 1993 ; utrobin et al .  1995 ) and supernova remnants ( e.g. , morse , winkler & kirshner 1995 ; aschenbach et al .  1995 ) in radio , optical , and x - ray bands which support the notion that supernova explosions are not spherically symmetric . the origin of the pulsar velocities is still unknown . two classes of mechanisms for the _ natal kicks _ have been suggested . the first class relies on convective instabilities in the collapsed stellar core and within the rebounding shock ( e.g. , burrows & fryxell 1992 ; burrows et al .  1995 ; janka & mller 1994 ,  1996 ; herant et al . the asymmetries in the matter and temperature distributions associated with the instabilities naturally lead to asymmetric matter ejection and/or asymmetric neutrino emission . numerical simulations indicate that the local , post - collapse instabilities are not adequate to account for kick velocities higher than a few hundred km  s@xmath7 ( janka & mller 1994 ) . a variant of this class of models therefore relies on the global asymmetric perturbations seeded in the presupernova cores ( goldreich et al . 1996 ; see also bazan & arnett 1994 ) . clearly , the magnitude of kick velocity depends on the degree of initial asymmetry in the imploding core ( burrows & hayes 1996 ) . due to various uncertainties in the presupernova stellar models , it is not clear at this point whether sufficiently large initial perturbations can be produced in the precollapse core ( lai & goldreich 1998 ) .    in this paper we focus on the second class of models in which the large pulsar kick velocities arise from asymmetric neutrino emission induced by strong magnetic fields . since @xmath16 of the neutron star binding energy ( a few times @xmath17  erg ) is released in neutrinos , tapping the neutrino energy would appear to be an efficient means to kick the newly - formed neutron star . magnetic fields are naturally invoked to break the spherical symmetry in neutrino emission . but the actual mechanism is unclear . we first briefly comment on previous works on this and related subjects . a number of authors have noted that parity violation in weak interactions may lead to asymmetric neutrino emission from proto - neutron stars ( chugai 1984 ; dorofeev et al .  1985 ; vilenkin 1995 ; horowitz & piekarewicz 1997 ) . however , their studies are largely unsatisfactory for a number of reasons : they either failed to identify the most relevant neutrino emission / interaction processes or the relevant physical conditions in proto - neutron stars , or stopped at estimating the magnetic field effects on neutrino opacities . chugai ( 1984 ) and vilenkin ( 1995 ) ( see also bezchastnov & haensel 1996 ) considered neutrino - electron scattering and concluded that the effect is extremely small ( e.g. , to obtain @xmath18  km  s@xmath7 would require a magnetic field of at least @xmath8  g ) . however , neutrino - electron scattering is less important than neutrino - nucleon scattering in determining the characteristics of neutrino transport in proto - neutron stars . similarly , dorofeev et al . ( 1985 ) considered neutrino emission by urca processes in strong magnetic fields . but as shown by lai & qian ( 1998 ) , in the bulk interior of the neutron star , the asymmetry associated with neutrino emission is cancelled by the asymmetry associated with neutrino absorption . concerning the parity violation effect in proto - neutron stars , horowitz & li ( 1997 ) recently suggested that the asymmetry in neutrino emission may be enhanced due to multiple scatterings of neutrinos by nucleons which are slightly polarized by the magnetic field . they estimated that a field strength of a few times @xmath19  g is adequate to account for kick velocities of a few hundred km  s@xmath7 . however , the paper by horowitz & li ( 1997 ) only discusses the idealized situations of scattering media , and ignores various essential physics that is needed in a proper treatment of neutrino transport in proto - neutron stars . in particular , it does not consider the effect of neutrino absorption which one might suspect to wash out the cumulative effect from multiple scatterings . our own study of the parity violation effect ( lai & qian 1998 ) was flawed in the treatment of scattering terms in the neutrino transport equation . as shown by arras & lai ( 1998 ) , detailed balance requires that there be no cumulative effect from multiple scatterings in the bulk interior of the proto - neutron star where local thermodynamic equilibrium applies to a good approximation . enhancement of neutrino emission asymmetry from multiple scatterings obtains only after neutrinos thermally decouple from proto - neutron star matter , and therefore is insignificant . bisnovatyi - kogan ( 1993 ) attributed pulsar kicks to asymmetric magnetic field distribution in proto - neutron stars . using the fact that neutron decay rate can be modified by the magnetic field , he inferred that neutrino emissions from opposite sides of the neutron star surface are different . however , neutron decay is not directly relevant for neutrino emission from a newly - formed neutron star . roulet ( 1997 ) ( whose paper was posted while this paper was being written ) considered the relevant neutrino absorption processes . but his study was restricted to calculating the neutrino cross sections as functions of neutrino energy ( corresponding to  2.1 of this paper ) , and therefore was insufficient for addressing the issue of asymmetric neutrino emission .    finally , it is worthwhile to mention a speculative idea on pulsar kicks which relies on nonstandard neutrino physics . kusenko & segr ( 1996 ) suggested that asymmetric @xmath20 emission could result from the mikheyev - smirnov - wolfenstein flavor transformation between @xmath20 and @xmath0 inside a magnetized proto - neutron star because a magnetic field changes the resonance condition for the flavor transformation . with a @xmath20 mass of @xmath21  ev , they claimed that magnetic fields of @xmath22  g can give the pulsar a kick velocity of a few hundred km  s@xmath7 . however , their treatment of neutrino transport was oversimplified : e.g. , they ignored that neutrinos of different energies have different resonance surfaces . furthermore , a simple geometric effect and realistic proto - neutron star conditions easily reduce their estimated kick velocity by an order of magnitude . most likely , a magnetic field strength of @xmath23  g is needed to produce the observed average pulsar kick velocity via their mechanism ( see qian 1997 ) . similarly strong magnetic fields are required in variants of the kusenko & segr mechanism , e.g. , those considered by akhmedov et al . ( 1997 ) ( which relies on both the neutrino mass and the neutrino magnetic moment to facilitate the flavor transformation ) . as should be clear from the brief review in  1.1 , in spite of quite a number of suggestions , there is at present no consensus on the magnitudes of the magnetic field induced asymmetry in neutrino emission from proto - neutron stars and the resulting kick velocities . in this paper we study the effect of asymmetric magnetic field topology on pulsar kicks . since the energy levels of @xmath2 and @xmath3 in a magnetic field are quantized , the @xmath0 and @xmath5 absorption opacities near the neutrinosphere and the neutrino - matter decoupling region depend on the local magnetic field strength . if the _ magnitude _ of magnetic field in the north pole is different from that in the south pole ( the field does not need to be ordered ) , then asymmetric neutrino flux may be generated . here we ignore the effects of magnetic fields on the proto - neutron star structure through the equation of state  these are secondary effects as far as neutrino transport is concerned . the organization of this paper is as follows . in  2 we calculate the @xmath0 and @xmath1 absorption cross sections as functions of neutrino energy in strong magnetic fields . section 3 summarizes the basic features of neutrino transport near the neutron star surface . in  4 we derive and evaluate the `` rosseland mean '' neutrino opacities which are directly related to the local neutrino flux . we then consider in  5 the change in neutrino flux due to the modified opacities , and estimate the kick velocity resulting from an asymmetric magnetic field topology . in  6 we briefly consider how the emergent neutrino spectra may be modified by the strong magnetic field . in  7 we present our conclusion . unless noted otherwise , we shall use units in which @xmath24 , @xmath25 , and the boltzmann constant @xmath26 are unity . in this section we consider how magnetic fields modify the ( total ) neutrino absorption ( @xmath27 and @xmath28 ) cross sections for specific neutrino energies . we present our calculation of the @xmath0 absorption cross section in detail . the @xmath1 absorption cross section can be obtained in a similar manner .    to begin with , it is useful to review various energy scales involved in the problem . near the neutrinosphere and the neutrino - matter decoupling region , the matter density is typically @xmath29@xmath19 g@xmath30@xmath31 , and the temperature is about several mev . the electrons are extremely relativistic and highly degenerate , with fermi energy @xmath32mev , where @xmath33 is the electron fraction ( the number of electrons per nucleon ) and @xmath34 . the neutrons and protons have fermi energies @xmath35^{2/3}~$]mev and @xmath36mev , respectively , and thus they are essentially nondegenerate .      in strong magnetic fields the transverse motion of the electron is quantized into landau levels . the electron energy dispersion relation is e_e=^1/2 ,  n=0 ,  1 ,  2 ,  where @xmath37 is the longitudinal momentum ( along the magnetic field ) , and @xmath38 is the landau level index . the dimensionless field strength @xmath39 is defined as b ,  with  b_c = m_e^2c^3e = 4.41410 ^ 13  g. the quantization effect on the proton is extremely small due to the large proton mass , and is neglected throughout this paper . since the cyclotron radius ( characterizing the size of the landau wavefunction ) @xmath40^{1/2}= [ \hbar/(m_e c)]b^{-1/2}=386\,b^{-1/2}$ ] fm is much greater than the range of the weak interaction [ @xmath41 fm , where @xmath42 gev is the mass of w boson ] , one can use the @xmath43 theory with contact interaction . locally ( at the point of interaction ) the landau wavefunction resembles a plane wave . thus we expect that the spin - averaged matrix element is the same as in the case without magnetic fields . this expectation has been borne out by detailed calculations ( e.g. , fassio - canuto 1969 ; matese & oconnell 1969 ) . therefore we only need to consider the effect of magnetic field on the electron phase space . since the neutrino energy ( a few to tens of mev ) of interest is much less than the nucleon mass , we first consider the case where the nucleon mass is taken to be infinity , an approximation adopted in the standard zero - field calculations ( e.g. , tubbs & schramm 1975 ) . the cross section for @xmath44 can be written as a sum over the ( spatial and spin ) states of the electron ( per unit volume ) ^_b(e_)=a_e(1-f_e)(e_+q - e_e ) , where we have defined ag_f^2 ^2_c ( c_v^2 + 3c_a^2 ) , and the other symbols have their usual meanings : @xmath45 is the universal fermi constant , @xmath46 refers to the cabbibo angle , @xmath47 and @xmath48 are the coupling constants , @xmath49 mev is the neutron - proton mass difference , @xmath50 is the @xmath0 energy , and @xmath51 is the pauli blocking factor . the electron occupation number @xmath52 is f_e=1 + 1 , where @xmath53 is the electron energy , @xmath54 is the electron chemical potential , and @xmath55 is the matter temperature . for zero field , the sum becomes @xmath56 , and we have ^_b=0(e_)=a^2(e_+q ) ^1/2(1-f_e ) , with @xmath57 . for @xmath58 , knowing that the degeneracy of the landau level ( per unit area ) is @xmath59 , with @xmath60 for @xmath61 and @xmath62 for @xmath63 , we have ^_b(e_)=a _ -^_n g_nm_e^2b2(1-f_e ) ( e_+q - e_e)=am_e^2b2 ^ 2 _ n=0^n_max g_n e_ep_z , n(1-f_e ) , [ absb]where @xmath57 , and @xmath64^{1/2}$ ] . the upper limit @xmath65 for the sum is the maximum value of @xmath38 for which @xmath66 is meaningfully defined , i.e. , n_max = int , [ nmax]where int@xmath67 $ ] stands for the integral part of @xmath68 . henceforth , @xmath65 is similarly defined , and we shall not explicitly write out the lower and upper limits for the sum . figure 1 depicts @xmath69 for @xmath70 and @xmath71 as functions of @xmath50 for typical conditions ( @xmath72 , @xmath73 mev , and @xmath74 mev ) in a proto - neutron star . note that ideally , one should compare @xmath69 and @xmath71 at the same matter density . for @xmath75 ( so that @xmath76 ) , the matter density is related to @xmath54 via @xmath77 in practice , we have found that for a given @xmath54 , the dependence of @xmath78 on the field strength is rather weak for the range of parameters of interest in this paper . thus in all of our calculations here and below , we compare the cross sections at a given @xmath54 . a prominent feature of @xmath69 shown in fig . 1 is the presence of `` resonances , '' which occur when @xmath79 ( cf . eq .  [ [ absb ] ] ) , i.e. , at the neutrino energies for which the final state electron lies in a `` stationary '' landau level ( with zero longitudinal momentum ) . these resonances reflect the divergence of the density of states of electrons at these stationary landau levels . the cross section for @xmath80 can be similarly calculated . one only needs to replace @xmath81 with @xmath82 and @xmath54 with @xmath83 in the expressions for @xmath0 absorption to obtain the desired expressions for @xmath1 absorption . the numerical results for @xmath1 absorption cross section are shown in fig .  2 for the same physical parameters as used in fig .  1 . the divergent cross sections at the `` resonances '' shown in figs .  1 and 2 are clearly unphysical . they are smoothed out by the thermal motion of the neutron . to the leading order , the width of the resonance is @xmath84 , corresponding to the doppler shift of the neutrino energy in the neutron s rest frame . the recoil energy , of order @xmath85 , is much less effective in smoothing out the resonances , and will be neglected . with this approximation , the smoothed cross section can be written as ^_b(e_)v_n^2 dv_nf_nd_n _ b^(e_ ) , [ aver]where @xmath86 and @xmath87 is the velocity distribution function for the neutron . if we replace the factor @xmath51 [ which depends on @xmath88 in @xmath89 with the average value @xmath90 , which can be approximately taken as @xmath51 evaluated at @xmath57 , then the angular integral in eq . ( [ aver ] ) can be carried out analytically . for @xmath91 , we have ^_b=0(e_)1-f_ev_n dv_nf_n13 ^ 2e _ ( e_e^2-m_e^2)^3/2|^e_e = e_(1+v_n)+q _ e_e = e_(1-v_n)+q , [ abs0aver]which is reduced to ^_b=0(e_)(e_+q)^21-f_e[abs0aver0]if we neglect the electron mass @xmath92 . for @xmath58 , we find has an integrable singularity at @xmath93 . ] ^_b(e_)1-f_ev_n dv_nf_nm_e^2b2 ^ 2e _ _ n g_n^1/2|^e_e = e_(1+v_n)+q _ e_e = e_(1-v_n)+q . [ absaver]we could have approximately included the neutron recoil effect by replacing @xmath81 with @xmath94 in eqs . ( [ abs0aver])([absaver ] ) . but we choose to leave out this effect since as mentioned previously , it is ineffective in smoothing out the resonances in @xmath69 . we have numerically evaluated eqs . ( [ abs0aver ] ) and ( [ absaver ] ) , and the results are depicted in fig .  1 together with the cross sections before averaging . similar results for the @xmath1 absorption cross section are shown in fig . as expected , we find that for @xmath91 , the averaged cross section @xmath95 is almost identical to @xmath71 ( see eq . [ [ abs0aver0 ] ] ) . however , for @xmath58 , the thermal averaging is essential in obtaining the smoothed physical absorption cross sections . the oscillatory behaviors of @xmath96 in figs .  1 and 2 can be demonstrated analytically . using the poisson summation formula ( see e.g. , ziman 1979 ) : _ n=0^f(n+12)= _ s=-^(-1)^s_0^f(x)e^2i s xdx , we can rewrite the sum in eq . ( [ absaver ] ) as @xmath97^{1/2}\nonumber\\ & & = { 2\over 3}(e_e^2-m_e^2)^{3/2 } + ( 2m_e^2b)^{3/2}{\rm re}\sum_{s=1}^\infty { 1-i\over 4\pi i s^{3/2}}\exp\left[2\pi i s \left({e_e^2-m_e^2\over 2m_e^2b } \right)\right ] . \label{summ}\end{aligned}\ ] ] the first term on the right - hand side of eq . ( [ summ ] ) corresponds to the zero - field result , and the second term ( involving the sum of many oscillatory terms ) corresponds to the correction due to the quantization effect . the sum is clearly convergent , and thus for small @xmath98 , the averaged cross section @xmath96 is reduced to the zero - field value @xmath95 .    for sufficiently large @xmath98 at a given @xmath50 , or , equivalently , for sufficiently small @xmath50 at a given @xmath98 , only the ground state landau level is filled by the electron , i.e. , @xmath99 ( see eq . [ [ nmax ] ] ) . in this limit , @xmath100 , and only the first term in eq . ( [ absb ] ) needs to be retained . we then have ^_b(e_)am_e^2b2 ^ 2 ( e_+q)^1/2(1-f_e ) [ sigbapp]for @xmath0 absorption . a similar expression for @xmath1 absorption is obtained by replacing @xmath81 with @xmath82 and @xmath54 with @xmath83 in eq . ( [ sigbapp ] ) . the physical conditions which figs .  1 and 2 are based on ( @xmath74  mev , @xmath72 , and @xmath73  mev , corresponding to a matter density of @xmath101 g  @xmath31 ) are typical of the neutrinosphere of the proto - neutron star . in fig . 3 we plot the @xmath0 and @xmath1 absorption cross sections for a number of slightly different physical parameters , specified by @xmath102 and @xmath55 ( all with @xmath72 ) . the cross section for elastic neutral - current scattering @xmath103 , @xmath104 ( e.g. , tubbs & schramm 1975 ) , is also shown for comparison . absorption dominates the transport of @xmath0 for high energies , with the ratio of absorption to scattering cross section per nucleon @xmath105 approaching 5 for @xmath106 . for smaller energies ( @xmath107 ) , absorption is severely suppressed by the pauli blocking of the final electron states , and scattering becomes more important . in the case of @xmath1 , pauli blocking is negligible , and scattering almost always dominates the transport , with the ratio of scattering to absorption cross section per nucleon @xmath108\simeq 2\,(0.1/y_e)$ ] for @xmath109  mev .    as we can see from figs .  1 - 3 , for a given set of physical parameters of the medium ( @xmath110 or @xmath78 ) , the magnetic field effect on the absorption opacity is more prominent for lower - energy neutrinos than for higher - energy ones . this is because for a given field strength , as the neutrino energy increases , more and more landau levels become available for the electron in the final state ( see eq . [ [ nmax ] ] ) . consequently , the electron phase space and hence the cross section should approach the zero - field values according to the correspondence principle . we found that typically when the number of landau levels filled by the electron ( @xmath65 ) is greater than 23 , the magnetic field effect on the absorption opacity becomes negligible . the `` critical '' neutrino energy is approximately given by e_^(crit)m_e(2bn_max+1)^1/2 6(n_max3)^1/2 ( b10 ^ 15g)^1/2   mev . for @xmath111 , one can approximate @xmath96 by the zero - field value @xmath95 with good accuracy ( within a few percent ) . alternatively , we can define a threshold magnetic field , @xmath112 , above which the @xmath0 and @xmath1 absorption opacities are modified ( by more than a few percent ) by the magnetic field : b_thb_c 310 ^ 15(e_10  mev)^2 ( 3n_max )  g , where @xmath113 mev is the average neutrino energy . in the deeper layers of the neutron star , the typical neutrino energy is higher . therefore we expect the magnetic field effect on the absorption opacities to be smaller there . in addition , for @xmath0 , absorption becomes unimportant as compared to scattering in the high density region ( although not necessarily in the deepest interior , where @xmath0 can be degenerate ) . thus , in order to assess the importance of the magnetic field effect on the neutrino emission , we only need to focus on the surface region , the region between the neutrinosphere and the neutrino - matter decoupling layer . we note that because the absorption opacities in the magnetic field for different neutrino energies oscillate around the the zero - field values , without a detailed calculation of the mean opacity ( appropriately averaged over the neutrino energy spectrum ) one can not even answer the qualitative question such as whether the magnetic field increases or decreases the local neutrino energy flux .    in   36 , we study how the modified absorption cross sections affect the emergent @xmath0 and @xmath1 energy fluxes and spectra from the proto - neutron star . as discussed in  2.4 , the quantization effect of magnetic fields on neutrino absorption is important in the region only near the neutron star surface ( between the neutrinosphere and the neutrino - matter decoupling layer ) . in this section we consider neutrino transport in this relatively low density region .    because @xmath0 and @xmath1 with typical energies of 10 mev have comparable absorption and scattering opacities ( see fig . 3 ) , the emergent @xmath0 and @xmath1 energy fluxes and spectra are determined essentially in the same region where these two neutrino species just begin to thermally decouple from the proto - neutron star matter . the @xmath0 and @xmath1 occupation numbers in this region are then approximately given by @xmath114 + 1 } , \label{occnue}\ ] ] and @xmath115 + 1 } , \label{occnueb}\ ] ] respectively , where @xmath116 , and @xmath117 is the @xmath0 chemical potential . from these occupation numbers , we obtain @xmath118 and @xmath119 where for example , @xmath120 and @xmath121 are the number and energy densities , respectively , for the @xmath0 . we can estimate the magnitude of @xmath122 as follows . the net deleptonization rate of the proto - neutron star is @xmath123 where @xmath124 is the initial electron fraction of the proto - neutron star , @xmath125 is the proto - neutron star mass , @xmath126 is the nucleon mass , @xmath127 is the neutrino diffusion timescale , and @xmath128 is the proto - neutron star radius . the sum of the @xmath0 and @xmath1 luminosities is @xmath129 from eqs . ( [ nuno])([nulum ] ) , we have @xmath130 where the numerical value is obtained for @xmath131 and @xmath132 km . therefore , to a good approximation , the @xmath0 and @xmath1 in the decoupling region have fermi - dirac energy distributions with zero chemical potentials . let @xmath133 be the specific intensity of neutrinos of a given species , where @xmath134 is the spatial position , and @xmath135 specifies the direction of propagation . the general transport equation takes the form i_e=_e^i_e^(fd ) + ^_e 14d i_e - i_e , [ inu]where @xmath136 , and the subscript `` @xmath137 '' indicates that the relevant physical quantities depend on the neutrino energy @xmath137 . the opacities @xmath138 and @xmath139 are related to the corresponding cross sections as , e.g. , @xmath140 and @xmath141 for @xmath0 . note that we have included the effect of stimulated absorption for neutrinos ( e.g. , imshennik & nadezhin 1973 ) by introducing the corrected absorption opacity @xmath142 $ ] ( @xmath55 is the matter temperature ) . hereafter , we shall suppress the superscript `` @xmath143 '' in @xmath138 , and the notation @xmath144 should be understood as having included the correction due to stimulated absorption . the quantity @xmath145 in eq . ( [ inu ] ) is the equilibrium neutrino intensity , given by the fermi - dirac distribution i_e^(fd)=e^3c^2h^31(e / t)+1 . note that the time derivative term @xmath146 on the left - hand - side of eq . ( [ inu ] ) has been neglected . this corresponds to an instantaneous redistribution of emission and absorption sources ( i.e. , an instantaneous redistribution of matter temperature )  the timescale for the redistribution is of order the mean free path divided by @xmath25 , much smaller than the neutrino cooling time of the star . the zeroth order moment of the transport equation can be written as _ e=_e^ c(u_e^(fd)-u_e ) , [ dfnu]where @xmath147 is the ( spectral ) neutrino energy flux , @xmath148 is the energy density , and @xmath149 is the corresponding equilibrium energy density . with the eddington closure relation , the first moment of eq . ( [ inu ] ) becomes _ e =- c3_e^(t)u_e , [ fnu]where @xmath150 . equations ( [ dfnu ] ) and ( [ fnu ] ) completely specify the neutrino radiation field . here we are concerned with the total neutrino energy flux emitted from a local surface region of the neutron star . the magnetic field strength in this local region is assumed to be constant . since the scale heights for physical quantities such as the matter density and temperature are much smaller than the local radii , we adopt the plane - parallel geometry . integrating eq . ( [ fnu ] ) over the neutrino energy , we find that the total neutrino energy flux ( of a specific species ) @xmath151 is given by f =- c3_rdudr , [ ftot]where @xmath152 is the total neutrino energy density , and @xmath153 is the `` rosseland mean opacity '' defined as = de(1/_e^(t))(d u_e / dr)d e(du_e / dr ) . [ ross]since neutrino emission and absorption rates are extremely fast even in the regions near the proto - neutron star surface , a radiative equilibrium , @xmath154 , @xmath20 , and @xmath155 is much weaker , the overall radiative equilibrium in the region near the proto - neutron star surface mainly involves matter and @xmath0 and @xmath1 . since the coupling between matter and @xmath0 is comparable to that between matter and @xmath1 , we further assume that an approximate radiative equilibrium holds separately for @xmath0 and @xmath1 . ] approximately holds , i.e. , _ e^ ( u_e^(fd)-u_e)de0 . > from eq .  ( [ dfnu ] ) , we have @xmath156 , or @xmath157 is approximately constant . if we introduce the rosseland mean optical depth measured from the stellar surface via @xmath158 , with @xmath157 constant , eq . ( [ ftot ] ) can be integrated to yield cu=3f(_r+23 ) , [ profile]where we have used the boundary condition @xmath159 at the surface ( @xmath160 ) . equation ( [ profile ] ) resembles the standard eddington profile . note that outside the decoupling region the neutrino energy density u may not be equal to @xmath161 specified by the matter temperature @xmath55 . similarly , @xmath162 may differ substantially from @xmath163 ( see  6 ) . however , we can still approximate @xmath162 by the fermi - dirac function at the neutrino temperature @xmath164 , which is equal to the matter temperature at the decoupling layer . equation ( [ ross ] ) can then be rewritten as = 307 ^ 4t_^5de1_e^(t ) e^4 ( e / t_)^2 . [ rosst]with the cross sections given in  2 , we can calculate the rosseland mean opacity @xmath165 in magnetic fields using eq . ( [ rosst ] ) . note that in eq . ( [ rosst ] ) @xmath166 depends on the local matter temperature @xmath55 . only in the regions not far from the decoupling layer , do we have @xmath167 . figure 4 shows the ratio of the rosseland mean opacity @xmath165 of @xmath0 to the zero - field value @xmath168 as a function of the magnetic field strength for typical conditions near the proto - neutron star surface . we choose these conditions to be @xmath169  mev , and @xmath170  mev , all with @xmath72 . figure 5 shows similar results for @xmath1 with similar physical parameters ( @xmath171  mev ; note that the results are insensitive to @xmath54 ) . we see that deviation of @xmath165 from the zero - field value by more than @xmath172 generally requires field strength @xmath173 . the behavior of the ratio @xmath174 as a function of @xmath39 is certainly not obvious a priori . since the rosseland mean integral is dominated by neutrinos with the smallest opacities , and since @xmath175 scales linearly with @xmath98 in the regime where only the ground landau level is filled ( see eq . [ [ sigbapp ] ] ) , we expect that for sufficiently large @xmath98 , the rosseland mean opacity increases linearly with increasing @xmath98 . this is indeed what is indicated in figs .  4 and 5 . for `` intermediate '' field strength , the @xmath175 curve for @xmath0 ( see fig .  1 ) has a significant portion which is below the zero - field value at low neutrino energies . for such intermediate fields , we would expect @xmath165 to be less than the zero - field value . this explains the non - monotonic behavior of @xmath174 as a function of @xmath98 shown in fig . 4 . for @xmath1 , the non - monotonic behavior is much less pronounced ( see the curve for @xmath176 mev in fig . 5 ) since the energy range in which @xmath177 is smaller . imagine a proto - neutron star with an asymmetric magnetic field topology , e.g. , the north pole and the south pole have different _ magnitudes _ of magnetic field ( say , by a factor of a few ) . what is the difference in the neutrino fluxes from these two regions due to the different absorption opacities of @xmath0 and @xmath1 , and what is the resulting kick velocity of the neutron star ? we recall that the fractional asymmetry @xmath178 in the total neutrino emission ( of all species ) required to generate a kick velocity @xmath179 is @xmath180 \,[e_{\rm tot}/(3\times 10^{53}~{\rm erg})]$ ] , where @xmath181 is the total neutrino energy radiated from the neutron star ( the neutron star mass has been set to @xmath182 ) . to quantitatively understand how the magnetically - modified opacities change the local neutrino flux , we will have to evaluate these opacities at different depths ( corresponding to different physical conditions ) below the neutron star surface , and to perform complete cooling calculations which take into account these opacities . however , figs . 4 and 5 clearly indicate that , even near the stellar surface , where one expects the modification to the opacities to be most prominent , and even with the most favorable conditions for the magnetic field to be effective in changing the opacities ( i.e. , low density and low temperature ) , a @xmath172 change in the rosseland mean opacity ( for @xmath0 or @xmath1 ) requires a magnetic field with @xmath183 , or @xmath184 g. in the following , we shall content ourselves with merely estimating the change in the local neutrino flux from the modified opacities . we first consider the case in which the matter temperature throughout the neutron star is unaffected by the magnetic field . one might imagine this as a possibility if convection dominates energy transfer between different parts of the star . under the neutrino - matter decoupling layer the specific neutrino energy density @xmath162 is the same as the fermi - dirac function @xmath163 specified by the matter temperature . the emergent neutrino flux is then inversely proportional to the rosseland mean optical depth @xmath185 from the decoupling layer to the stellar surface ( see eq .  [ [ profile ] ] ) . since @xmath185 for @xmath0 or @xmath1 at the decoupling layer is of order a few , the flux @xmath157 is inversely proportional to @xmath153 evaluated near the decoupling region . if we use the curve for @xmath74  mev and @xmath186  mev in fig .  4 ( which has the most pronounced magnetic correction ) , we find that the fractional change in the @xmath0 flux is @xmath187 at @xmath70 ( if we used the curve for @xmath188 or 20  mev , which more closely represents the physical condition at the decoupling region , the modification @xmath189 would be smaller by a factor of a few or more ) . similarly , using the curve for @xmath74  mev in fig .  5 to maximize the magnetic correction , we find @xmath190 at @xmath70 . since the fluxes of other types of neutrinos are unchanged and @xmath191 depend on the @xmath192-@xmath193 scattering cross sections , which can be affected by the magnetic field . however , the layers are located much deeper than those of @xmath0 and @xmath5 ( since @xmath192-@xmath193 scattering opacity is small ) , the modification to the @xmath192-@xmath193 scattering opacity due to electron landau levels is therefore smaller than the modification to the absorption opacities of @xmath4 and @xmath5 considered in this paper . ] , and since different neutrino species more or less carry away the same amount of energy , the change in the total neutrino flux is only @xmath194 at @xmath70 . if we now imagine that the north pole of the neutron star has @xmath195 while the south pole has a field strength a factor of a few smaller , then , taking into account the geometric reduction ( @xmath196 ) , the flux asymmetry is about @xmath197 , just enough to give a kick velocity of @xmath21  km  s@xmath7 to the neutron star . with @xmath153 evaluated at more realistic conditions , we expect the resulting @xmath178 to be even smaller . the assumption that the matter temperature is unaffected by the magnetic field ( and therefore spherically symmetric ) is unlikely to be valid . in reality , when a region of the star ( say the north pole ) emits more flux , it also cools faster . with reduced temperature , this region will tend to emit less flux . in the extreme case , one would find @xmath198 , with @xmath199 the neutrino energy density in the stellar core , and @xmath185 the rosseland mean optical depth from the core to the surface . since in the bulk interior of the star , modification to @xmath153 due to the magnetic field is negligible , we would conclude that the change in the neutrino flux is much smaller than what is estimated in the last paragraph .    to summarize , the modifications to the absorption opacities of @xmath0 and @xmath5 due to quantized landau levels give rise to rather small change in neutrino fluxes from the neutron star . to generate appreciable kick velocities ( @xmath6  km  s@xmath7 ) from this effect requires the difference of field strengths at the two opposite poles to be at least @xmath8  g ( and possibly much larger ) . we now consider how the oscillatory behaviors ( see figs . 13 ) of the absorption opacities of @xmath0 and @xmath1 may manifest themselves in the radiated neutrino energy spectra . for an absorption dominated medium , neutrinos with energy @xmath137 decouple from the matter and start free - streaming at a column depth ( in g  @xmath200 ) @xmath201 . thus the neutrino energy spectrum is given by @xmath202 with @xmath203 evaluated at the depth @xmath204 . in the presence of significant scattering , decoupling occurs at depth @xmath205 with the effective opacity @xmath206 , whereas the neutrinosphere from which neutrinos start free - streaming is located at @xmath207 . near the decoupling region we have @xmath208 , and thus from eq . ( [ fnu ] ) the neutrino spectral energy flux is given by f_e~_e^eff_e^(t ) cu_e^(fd)(y~1/_e^eff ) . the above qualitative argument can be made more rigorous by directly solving eqs . ( [ dfnu ] ) and ( [ fnu ] ) . consider a plane - parallel atmosphere and define the ( energy - dependent ) optical depth via @xmath209 . equations ( [ dfnu ] ) and ( [ fnu ] ) can be combined to give u_e=_e^2 [ u_e - u_e^(fd ) ] ,   with  _ e^2=3_e^_e^(t ) . with a constant @xmath210 , the above equation can be solved exactly . for the boundary condition @xmath211 at @xmath212 , the neutrino energy flux at the surface is given by f_e=_e/31 + 2_e/3_0^d_e_ecu_e^(fd)(_e)(-_e_e ) cu_e^(fd)(_e_e^-1 ) . [ fnu2 ]    figure 6 shows the spectral energy fluxes of @xmath0 calculated using eq . ( [ fnu2 ] ) for different field strengths ( @xmath70 and @xmath213 ) . the matter temperature profile is obtained from eq . ( [ profile ] ) assuming @xmath214 . for a total energy flux of @xmath215 ( where @xmath216 is the effective temperature ) , we have @xmath217^{1/4}$ ] . we choose @xmath218 mev in these examples . for simplicity , we have neglected the variations of the opacities as a function of @xmath219 . thus the temperature at decoupling ( @xmath220 ; see eq . [ [ fnu2 ] ] ) is given by @xmath221^{1/4}$ ] , where @xmath222 and we take the relevant opacities to have the corresponding values for @xmath72 , @xmath73 mev , and @xmath74 mev . this simplification should not introduce qualitative changes in the spectrum . as seen from fig . 6 , the spectrum is pinched at the high energies because of the @xmath223-dependence of the opacities . the magnetic field introduces modulations in the spectrum . for lower magnetic field strengths , the modulations occur at the lower energy band . as the magnetic field gets stronger , the modulations shift to the higher energy part of the spectrum . the distortions and modulations of the @xmath0 and @xmath1 energy spectra by strong magnetic fields may offer an interesting possibility to limit the magnetic field strength in proto - neutron stars through the detection of neutrinos from a galactic supernova . we note , however , that the @xmath0 spectra shown in fig .  6 are based on very approximate calculations . they serve only as illustrative examples of the magnetic field effect on the neutrino spectra . more detailed calculations are needed to obtain realistic spectra which include the magnetic field effect . in this paper , we have studied the issue of whether asymmetric magnetic field topology in a proto - neutron star can induce asymmetric neutrino emission from the star through the modifications of the neutrino absorption opacities by the magnetic field . these modifications arise from the quantized landau levels of electrons and positrons produced in strong magnetic fields . by calculating the appropriate mean neutrino opacities in the magnetic field , we demonstrate that this mechanism is rather inefficient in generating a kick to the neutron star : to obtain appreciable kick velocities ( @xmath6  km  s@xmath7 ) , the difference in the field strengths at the two opposite poles of the star must be at least @xmath8  g.    this work was started while d.l . held the richard c. tolman postdoctoral fellowship at caltech . additional support was provided by nasa grant nag 5 - 2756 and nsf grant ast-9417371 , as well as a start - up fund for new faculty at cornell university . is supported by the david w. morrisroe postdoctoral fellowship at caltech .
Abstract:

Answer: in proto - neutron stars with strong magnetic fields , the cross section for @xmath0 ( @xmath1 ) absorption on neutrons ( protons ) depends on the local magnetic field strength due to the quantization of energy levels for the @xmath2 ( @xmath3 ) produced in the final state . if the neutron star possesses an asymmetric magnetic field topology in the sense that the magnitude of magnetic field in the north pole is different from that in the south pole , then asymmetric neutrino emission may be generated . we calculate the absorption cross sections of @xmath4 and @xmath5 in strong magnetic fields as a function of the neutrino energy . these cross sections exhibit oscillatory behaviors which occur because new landau levels for the @xmath2 ( @xmath3 ) become accessible as the neutrino energy increases . by evaluating the appropriately averaged neutrino opacities , we demonstrate that the change in the local neutrino flux due to the modified opacities is rather small . to generate appreciable kick velocity ( @xmath6  km  s@xmath7 ) to the newly - formed neutron star , the difference in the field strengths at the two opposite poles of the star must be at least @xmath8  g. we also consider the magnetic field effect on the spectral neutrino energy fluxes . the oscillatory features in the absorption opacities give rise to modulations in the emergent spectra of @xmath0 and @xmath1 .