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 crab pulsar ( psr b0531 + 21 ) is a young radio pulsar with a spin frequency @xmath2 and frequency derivative @xmath3 ( lyne et al . its characteristic magnetic field is about @xmath4 at the magnetic poles ] . the crab pulsar has been monitored continuously for decades years by various telescopes since it was discovered in 1968 ( lyne et al . 1993 ; lyne et al . 2015 ; wang et al . its braking index is @xmath5 ( lyne et al . different braking index values @xmath6 , @xmath7 and @xmath8 are also reported between glitches ( wang et al . 2012 ; lyne et al . observational details are given in table [ parameters ] . long - term observations have found that almost all pulsars ( except accreting x - ray pulsars in binary systems ) are spinning down . the spin - down behavior can be described by a power law ( lyne et al . 1993 ) : @xmath9 here , @xmath10 is usually taken as a constant and @xmath11 is the braking index . the braking index and second braking index are defined respectively ( livingstone et al . 2005 ) : @xmath12 @xmath13 where @xmath14 is the pulsar spin frequency , @xmath15 , @xmath16 and @xmath17 are the first , second and third frequency derivative , respectively . the spinning down of pulsars is usually assumed to be braked by the magnetic dipole radiation ( shapiro & teukolsky 1983 ) . in the magnetic dipole radiation model , a neutron star rotates uniformly in vacuo at a frequency @xmath14 and possesses a magnetic moment @xmath18 . the corresponding slowdown rate is : @xmath19 where @xmath20 is the magnetic dipole moment ( @xmath21 is the polar magnetic field and @xmath22 is the radius of neutron star ) , @xmath23 is the moment of inertia , @xmath24 is the speed of light , and @xmath25 is the angle between the rotational axis and the magnetic axis ( i.e. , the inclination angle ) . the spin down behavior of pulsar in this model can be described as @xmath26 . the braking index is exactly three if @xmath27 , and @xmath25 are constant . to date , only eight pulsars , including the crab pulsar , have measured the meaningful braking indices for they are young and own relatively larger @xmath15 ( espinoza et al . 2011 ; lyne et al . their braking indices are all smaller than three . it means that there are other physical processes needed to slow down the pulsar . many mechanisms have been proposed in order to explain the braking index observations , e.g. , the pulsar wind model ( xu & qiao 2001 ; wu et al . 2003 ; contopoulos & spitkovsky 2006 ; yue et al . 2007 ) , a changing magnetic field strength ( lin et al . 2004 ; chen & li 2006 ; espinoza et al . 2011 ) , a changing inclination angle ( lyne et al . 2013 ) , and additional torques due to accretion ( liu et al . 2014 ) . however , these models are more or less not consistent with observations or can not simulate the long - term evolution of pulsars . and these models can not explain the different braking indices detected between glitches ( wang et al . 2012 ; lyne et al . furthermore , the effect of pulsar death ( death line , ruderman & sutherland 1975 ; or death valley , chen & ruderman 1993 ; zhang et al . 2000 ) should be considered in modeling the long term rotational evolution of the pulsar ( contopoulos & spitkovsky 2006 ) . the pulsar wind model considers both the pulsar spin - down and the particle acceleration in the magnetosphere ( xu & qiao 2001 ; wu et al . 2003 ) . in this paper , an updated pulsar wind model is built based on previous researches ( xu & qiao 2001 ; yue et al . 2007 ; contopoulos & spitkovsky 2006 ; li et al . it includes : ( 1 ) both magnetic dipole radiation and particle outflow , and their dependence on the inclination angle ; ( 2 ) the particle outflow depends on the specific acceleration model ; ( 3 ) the primary particle density may be much larger than the goldreich - julian charge density ; ( 4 ) the effect of pulsar death is considered in modeling the rotational evolution of the pulsar . this model can calculate both the short term and long term evolution of pulsars . it is applied to the crab pulsar which has the most detailed timing observations . possible applications to other sources are also illustrated . the pulsar wind model and model calculations of the crab pulsar are presented in section 2 . discussions and conclusions are given in section 3 and section 4 , respectively . [ parameters ]    .observations of the crab pulsar . [ cols="<,<",options="header " , ]     notes :    ( a ) : @xmath28 is the coefficient of primary particle density @xmath29 .    ( b):@xmath25 and @xmath30 are respectively the inclination angle and the magnetic field . ( c ) : @xmath31 is the initial period . ( d):@xmath32 and @xmath33 are the transition period and corresponding age . at this transition point the braking index is @xmath34 . because the minimum braking index of model sclf(ics ) is @xmath35 , there is no such a transition point , which can be seen in the evolution in @xmath1 diagram ( see figure [ figevo&models ] ) . ( e ) : @xmath36 and @xmath37 are the period and age when pulsar stops radio emission . the death period is calculated by equation ( [ pdeath ] ) . a braking index @xmath8 has been monitored after the removal of the effects of glitches during an epoch when there are many glitches occured ( lyne et al . 2015 ) . it is different from the long term underlying braking index @xmath38 of the crab pulsar . previously , wang et al . ( 2012 ) have measured two different values of braking index @xmath39 and @xmath40 between glitches . they attribute it to the effect of varying particle wind strength ( wang et al . it is generally accepted that glitch is caused by the inner effect but may lead to some effect in the outer magnetosphere , for example : the changing of primary particle density during this epoch . in the pulsar wind model , it can be denoted by the changing coefficient of particle density @xmath41 . then equation ( [ nwind ] ) should be modified : @xmath42 equation ( [ nkappa ] ) can be rewritten as : @xmath43 where @xmath44 is the characteristic age , @xmath45 can be viewed as the typical time scale of the changing particle density . when the magnetosphere is in equilibrium ( i.e. , when the glitch activities of the pulsar is not very active ) , the time scale of particle density variation ( @xmath46 ) may be very large ( i.e. , larger than the characteristic age @xmath47 ) . the third term in equation ( [ ntauk ] ) does not contribute . but , if they are comparable with each other , the braking index changes substantially . in other words , when the outflow particle density increases and the effect of particle wind becomes stronger , means @xmath48 or @xmath49 , the braking index will be smaller than @xmath38 . but when the outflow particle density decreases ( @xmath50 or @xmath51 ) , the braking index will be larger than @xmath38 . as the out - flow particle density tending to a certain value ( larger or smaller than previous value ) , the braking index tends to its underlying value ( @xmath38 ) . generally , the increased ( or decreased ) component of the outflow particle density ( may be @xmath52 ) is so much small that it can be ignored . figure [ fignk ] shows the braking index as function of @xmath46 for the crab pulsar in the vg(cr ) model . in order to better understand the observations of @xmath53 ( lyne et al . 2015 ) , we consider the @xmath48 only . as we can see in this figure , the braking index is insensitive to @xmath46 when it is larger than @xmath54 . but , when it is comparable with the characteristic age ( about @xmath0 yr ) , the braking index decreases sharply . when braking index is @xmath8 , the changing rate of the particle density is @xmath55 . in this interval ( from mjd 51000 to mjd 53000 , about @xmath56 years ) , the particle density has changed by @xmath57 . the change of particle density will result in changes of radiation energy loss rate and period derivative . their changes respectively are : @xmath58 and @xmath59 . the pulsar wind model can be tested with observations ( jodrell bank monthly ephemeris of the crab pulsar ) . the effect of glitches is not taken into consideration in the pulsar wind model . the changes of period and period derivative caused by glitches should be added ( lyne et al . 2015 ) . figure [ figob&model ] is the model calculation compared with observations ( without modelling the transient variation caused by glitches ) . when the amplitude of glitches is relatively small , the change of particle density is not obvious , then the model calculations ( the blue points ) fit well with observations ( the black points ) . however , when the amplitude of glitches is large , the effect of glitches must be added ( the red points ) . moreover , the effect of particle density change should be taken into consideration . a factor @xmath60 ( which is calculated above ) is added to the red points ( from mjd 51804.75 afterwards ) , shown as the green points . the green points can match the general trend of the observations .     in the vg(cr ) model . the dashed line is the underlying braking index @xmath38 . the dotted line is the smaller braking index @xmath8 measured during glitch activities ( lyne et al . 2015 , figure 7 there).,scaledwidth=45.0% ]      the pulsar period and its fist derivative evolution with time in the vg(cr ) model are shown in the first and second panel of figure [ figparameters ] . the curves evolve slowly in its early age and then rise sharply , indicating that the crab pulsar is mainly braked by magnetic dipole radiation firstly and then mainly by particle wind . the transition period is @xmath61 and age of the crab pulsar will be @xmath62 calculated in the vg(cr ) model ( see table [ tab_calculations ] ) . the braking index evolution with time in the vg(cr ) model is shown in the third panel of figure [ figparameters ] . braking index decreases from @xmath63 to @xmath64 different minimum braking indices for different acceleration models . the bottom panel of figure [ figparameters ] shows the second braking index evolution with time . the curve gradually decreases from @xmath65 to 30/49 . figure [ figevo&short ] shows the evolution of the crab pulsar in @xmath1 diagram from birth to its @xmath66 years old . these points are respectively taken when the crab pulsar is : @xmath67 , @xmath68 , @xmath62 , @xmath69 and @xmath66 years old . period derivative evolves with period by a function of @xmath70 ( espinoza et al . 2011 ) . in the log - log plot , the evolution curve evolves with a slop ( @xmath71 ) . as the braking index falling from @xmath63 to @xmath67 , the pulsar evolve from magnetic dipole radiation dominated case ( the left points ) to the particle wind dominated case ( the right points ) and the bottom is the transition point ( @xmath72 ) . clearly , the curve evolves more sharply after the transition point . the asymptotic behaviors are discussed in the appendix b.             +            -@xmath73 diagram . five points are chosen when the pulsar is : @xmath67 , @xmath68 , @xmath62 , @xmath69 and @xmath66 years old.,scaledwidth=45.0% ]      as shown is the figure [ figevo&short ] , the line evolves up right with a slope about @xmath67 after the transition point just like psr j1734 - 3333 whose braking index is @xmath74 ( espinoza et al . if the pulsar continually evolve in this trend , it would not dead and might end itself as a magnetar . it will be hard to explain the observations for : ( 1 ) the number of magnetars are smaller than radio pulsars and most of the pulsars occupy the central region in the @xmath75-@xmath73 diagram ( see figure [ figevo&long ] ) ; ( 2 ) only a small portion of magnetars are radio emitters ( olausen & kaspi 2014 ) ; ( 3 ) the properties of the crab pulsar , psr j1734 - 3333 and other high magnetic field pulsars are significantly different from magnetars ( ng & kaspi 2011 ) . as pulsars slow - down , the density of outflow particles may decrease and the effect of particle wind will recede ( for the sclf cases , the number of particles ( or @xmath28 ) may close to constant , but the particle density still decreases because a reducing goldreich - julian charge density ) . therefore , the effect of pulsar death must be included in modeling the long - term rotational evolution of pulsars . the essential condition of radio emission for a pulsar is its pair production ( sturock 1971 ; ruderman & sutherland 1975 ) . the so called pulsar `` death '' means the stopping of pulsar emission ( ruderman & sutherland 1975 ; chen & ruderman 1993 ; zhang et al . 2000 ; contopoulos & spitkovsky 2006 ) . the death defined by ruderman & sutherland ( 1975 ) is that : a pulsar dead when the maximum potential drop @xmath76 available from the pulsar is smaller than the acceleration potential needed to accelerate particles . because this death line is model dependent , chen & ruderman ( 1993 ) proposed death `` valley '' by modifying the boundary conditions . and this death `` valley '' was updated by considering different particle acceleration and photon radiation processes ( zhang et al . however these works are just to separate pulsars with radio emission from those without . contopoulos & spitkovsky ( 2006 ) included the effect of pulsar death when modeling the rotational evolution of pulsars . however , the braking indices there are always larger than three . considering particle acceleration and pulsar death simultaneously may give a comprehensive interpretation for both short and long - term rotational evolution of pulsars . the visual explanation for pulsar death in the pulsar wind model is that the particles in the magnetosphere are exhausted and there are no primary particles generated . for the physical understanding , we can refer to the equation ( [ edottotal ] ) , expressions @xmath77 and @xmath78 ( i.e. , @xmath79 , xu & qiao 2001 ) , with the pulsar spinning down , the maximum potential of the pulsar ( @xmath80 ) drops , when the maximum potential can not meet the need to accelerate particles ( @xmath81 ) , the pulsar death . besides , the acceleration potential is weakly dependent on the @xmath82 , that is why we can make calculations under the constant acceleration potential ( the cap model ) ( yue et al . 2007 ) . figure [ figevo&long ] shows the long - term evolutions of the crab pulsar in the vg(cr ) model . the dashed curve is the line shown in figure [ figevo&short ] and the solid is the one with pulsar death ( according to equation ( [ eqn_edotdeath ] ) ) . these two lines evolve similar at the early age , and divide obviously after the transition point . the second line falls down afterwards and moves down - right with a slop about @xmath83 finally . the initial period and transition period can be calculated according the second line , which are same with the values given by the first line within precision . indicating that the effect of pulsar death can be ignored in the early time . as mentioned before , the particle density decreases as pulsar spin - down . the gap between these two lines is caused by the reducing of particle density . when the spin period of crab approaches to the death period @xmath84 , the radio emission tends to stop and the slop of the second line is very large . crab evolution in various acceleration models are given in figure [ figevo&models ] . the bottom line is the evolution curve in sclf(ics ) model . there is no transition point , indicating that the particle wind effect of this model is very weak . the top line is the evolution of the crab pulsar in the og model . this model is for high - energy radiation , the death of this model may differ from radio emission models . the results here for the og model are just crude approximations to show the effect of pulsar death . evolution of the crab pulsar in other models are similar with each other . diagram has listed all known radio pulsars , magnetars , and millisecond pulsars ( updated from figure 1 in tong & wang 2014 ) . the dot - dashed line represents the fiducial death line ( ruderman & sutherland 1975 ) . pulsars with meaningful braking indices ( espinoza et al . 2011 ; lyne et al . 2015 ) are defined by large square . the arrows represent their evolution directions and each arrow indicate its motion in the next 10000 yr ( for the vela pulsar , 20000 yr is used ; 5000 yr for psr j1846 - 0258 and psr j1119 - 6127).,scaledwidth=45.0% ]     for explanations.,scaledwidth=45.0% ] in the magnetic dipole radiation model , a pulsar is assumed as a magnetic dipole rotating in vacuum . it is just a low - order approximation of the true magnetosphere , and there must be particles acceleration and generation process so that the radio emission can be received . the actual magnetosphere surrounding the pulsar must be filled with plasma and corotate with pulsar because of magnetic freeze . according the special relativity , the corotation magnetosphere can not infinitely extend but only exist within the light cylinder . the magnetosphere is divided into two parts by the light cylinder : the closed magnetic field lines and the open magnetic field lines ( including polar gap and outer gap ) . in the open field line region , particles are accelerated , generate radiation and flow out ( thus taken away some of the rotational energy of the central neutron star ) . in the pulsar wind model , the pulsar is braked down by magnetic dipole radiation and particle wind . but in the electric current loss model , the electric current flow generates a spin - down torque ( harding et al . 1999 ; beskin et al . based on the same magnetospheric physics , the pulsar wind model and the current loss case are just different point of views , will produce the same results ( e.g. section 3.2 in tong et al . 2013 ) . and they can be used to check each other . it is meaningful to take advantage of the existing information and explain more observations . the particle wind worked on pulsar spin - down was first proven by observations of intermittent pulsar psr b1931 + 24 whose radio emission switches between `` on '' and `` of '' state ( kramer et al . corresponding calculations for intermittent pulsars in pulsar wind model were made by li et al . ( 2014 ) .    as one of the best observed pulsar , observations of the crab pulsar are relatively comprehensive ( see table [ parameters ] @xmath85 timing results from jodrell bank website ) . it will be the best target to test radiation models . the pulsar parameters can be calculated by equations ( [ edot ] ) , ( [ nwind ] ) and ( [ mwind ] ) when model dependent acceleration potential ( or @xmath86 ) is given ( see table [ expressions ] ) . the goldreich - julian charge density is taken as the primary particle density in these models of xu & qiao ( 2001 ) and wu et al ( 2003 ) . however , yue et al . ( 2007 ) constrained the primary particle density by observations of several pulsars and predicted that it should be @xmath87 times of goldreich - julian charge density . in our calculations , the primary particle density is at least @xmath88 times of goldreich - julian charge density . in the vg(cr ) model , a particle density @xmath89 is chosen . when applying the pulsar wind model to intermittent pulsars , a goldreich - julian charge density is assumed ( li et al . as the pulsar spin - down , its out - flow particle density decreases . in the @xmath1 diagram , these intermittent pulsars are on the right of those young pulsars ( see figure [ figevo&long ] ) . this may explain the different particle density in the crab pulsar and intermittent pulsars . other studies also found the need for a high particle density in the pulsar magnetosphere , e.g. , when modeling the eclipse of the double pulsar system , giant pulses of the crab pulsar , and magnetar x - ray spectra ( lyutikov 2008 and references therein ) , optical / ultraviolet excess of x - ray dim isolated neutron stars ( tong et al . 2010 , 2011 ) . magnetohydrodynamical simulations of the crab nebula also require a much higher electron flow from the magnetosphere ( buhler & blandford 2014 and references therein ) . besides , a super goldreich - julian current has been proposed in a series theoretical works about the pairs creations in the acceleration gap ( timokhin & arons , 2013 ; zanotti et al . 2012 ; hirotani 2006 ) . combined with the results in this paper , it is possible that there are high density plasmas ( e.g. , @xmath90-@xmath91 times the goldreich - julian density ) in both open and closed field line regions in the pulsar magnetosphere .      here the updated pulsar wind model is employed to compare with the observations of the crab pulsar . it can also be applied to other sources . for eight pulsars with meaningful braking index measured ( espinoza et al . 2011 ; lyne et al . 2015 ) , their magnetospheric parameters can also be calculated , including magnetic field , inclination angle , particle density etc . but the higher second frequency derivatives of some pulsars ( hobbs et al . 2010 ) may be attributed to the fluctuation of magnetosphere ( contopoulos 2007 ) . during this process , more information will be helpful to constrain the model parameters , e.g. , inclination angle , second braking index . their evolution on the @xmath1 diagram can be calculated , which will be similar to figure [ figevo&long ] . the pulsar psr j1734@xmath923333 has braking index @xmath93 ( espinoza et al . its rotational energy loss rate will be dominated by the particle wind . also from figure [ figevo&long ] , the braking indices of young pulsars will naturally be divided into two groups . braking indices of pulsars in the first group are close to three , their evolution directions are down to the right in the @xmath1 diagram ( @xmath94 ) . braking indices in the second group are close to one and their evolution directions are up to the right ( @xmath95 ) . generally speaking , sources in the second group will be older than the ones in the first group . at present , there are some hints for the evolution of pulsar braking index ( espinoza 2013 ) . future more observations will deepen our understanding on this aspect . glitch is a phenomenon of the pulsar suddenly spinning up and then slowing down by a larger rate . in the pulsar wind model , the lower braking index @xmath8 during glitches ( of the crab pulsar ) is due to the larger out - flow particle density . the increasing particle density is a possible performance associated with glitches . glitches induced magnetospheric activities are very common in magnetars ( kaspi et al . 2003 ; dib & kaspi 2014 ) . for the high magnetic field pulsar psr j1846@xmath920258 , it also shows a smaller braking index after glitch ( livingstone et al . the cause of a smaller braking index may be similar to the crab pulsar case . due to a larger particle outflow , the pulsar should experience some net spindown after glitch . this has already been observed ( livingstone et al . 2010 ) . some primary calculations for psr j1846@xmath920258 has been done previously ( tong 2014 , based on a wind braking model designed for magnetars ) . the smaller braking index of psr j1846@xmath920258 after glitch is consistent with the results here .    from figure [ figevo&long ] and figure [ figevo&models ] , the pulsar will go up - right in the @xmath1 diagram in the wind braking dominated case . therefore , for pulsars in the up right corner of the @xmath1 diagram ( e.g. , high magnetic field pulsars and magnetars ) , their true dipole magnetic field may be much lower than the characteristic magnetic field . for magnetar sgr 1806@xmath9220 , its characteristic magnetic field is about @xmath96 ( tong et al . this dipole magnetic field is too high to understand comfortably ( vigano et al . 2013 ) . using the physical braking mechanism in this paper , the dipole magnetic field of sgr 1806@xmath9220 will be much lower . in the late stage of a pulsar , it will bent and go downward in the @xmath1 diagram . the crab pulsar and other pulsars will not evolve into the cluster of magnetars . furthermore , for magnetars , they also will go downward in the @xmath1 diagram when they are aged enough . this may corresponds to the case of low magnetic field magnetars ( tong & xu 2012 ) . an updated pulsar wind model which considered the effect of particle density and pulsar death is developed . it is employed to calculate parameters and simulate the evolution of the crab pulsar . the braking index @xmath97 is the combined effect of magnetic dipole radiation and particle wind . the lower braking indices measured between glitches are caused by the increasing particle density . by adding the glitch parameters and a change of period derivative caused by the change of particle density , the theoretical evolution curve is well fitted with observations ( figure [ figob&model ] ) . this may be viewed as glitch induced magnetospheric activities in normal pulsars . giving the model dependent acceleration potential drop , the magnetic field , inclination angle and particle density of the crab pulsar can be calculated . the evolution of braking index and the crab pulsar in @xmath1 diagram can be obtained . in the @xmath1 diagram , the crab pulsar will evolve towards the death valley , not to the cluster of magnetars ( figure [ figevo&long ] ) . different acceleration models are also considered . the possible application of the present model to other sources ( pulsars with braking index measured , and the magnetar population ) is also mentioned . the authors would like to thank the referee barsukov d.p . for helpful comments and r.x.xu , j.b.wang for discussions . h.tong is supported xinjiang bairen project , west light foundation of cas ( lhxz201201 ) , qing cu hui of cas , and 973 program ( 2015cb857100 ) . the treatment of contopoulos & spitkovsky ( 2006 ) is employed in the following . assuming that the open field line regions rotate at a angular velocity @xmath98 but not the angular velocity of the neutron star . @xmath99 the particle density in this acceleration gap is @xmath28 times the goldreich - julian charge density : @xmath100 the energy taken away by the out - flowing particles is : @xmath101 the perpendicular component and parallel component of the magnetic moment are respectively related to the magnetic dipole radiation and particle acceleration . the rotation energy loss rate can be written as : @xmath102 with the expression of @xmath86 @xmath103 for a constant acceleration potential @xmath104 , correspondingly , @xmath105 for the physical acceleration models , @xmath86 is model dependent . its expression in the vg(cr ) case is calculated in the following . the potential drop of open field lines is ( ruderman & sutherland 1975 ) : @xmath106 where @xmath98 is the angular frequency of the open field lines ( contopoulos & spitkovsky 2006 ) , @xmath107 is the thickness of the gap , @xmath108 and @xmath109 is the curvature radius of the open field lines ( @xmath110 ) ( ruderman & sutherland 1975 ) . the expression of @xmath86 can be rewritten as : @xmath111 the exponent of @xmath112 is the minimum braking index in the vg(cr ) model , @xmath113 ( li et al . 2014 ) . in the wind braking model , the minimum braking index is always around one ( li et al . 2014 ) . therefore , it is taken as one in equation ( [ eqn_edotdeath ] ) ( for all the acceleration models ) . during the numerical calculations , the neutron star radius is taken as @xmath114 ( i.e. , @xmath115 ) . in the pulsar wind model , the pulsar is braked down by the combination of magnetic dipole radiation and particle wind . the effect of particle wind is more and more important with pulsar spinning down . as shown in figure [ figparameters ] , the period and period derivative evolve up sharply with age . if the pulsar braking is dominated by magnetic dipole radiation , equation ( [ nudotpowerlaw ] ) can be written as : @xmath116 the evolution of spin frequency is @xmath117 ( assuming initial spin frequency is very large ) . correspondingly , the period and period derivative evolution with time are respectively : @xmath118 , and @xmath119 . if the pulsar is braked down mainly by particle wind , the braking index is about @xmath67 . the spin - down power law is : @xmath120 the evolution of spin frequency is @xmath121 . correspondingly , the period and period derivative evolution with time are respectively : @xmath122 , and @xmath123 . clearly , the period and period derivative evolve faster by the exponential form .
Abstract:

Answer: the pulsar wind model is updated by considering the effect of particle density and pulsar death . it can describe both the short term and long term rotational evolution of pulsars consistently . it is applied to model the rotational evolution of the crab pulsar . the pulsar is spun down by a combination of magnetic dipole radiation and particle wind . the parameters of the crab pulsar , including magnetic field , inclination angle , and particle density are calculated . the primary particle density in acceleration region is about @xmath0 times the goldreich - julian charge density . the lower braking index between glitches is due to a larger outflowing particle density . this may be glitch induced magnetospheric activities in normal pulsars . evolution of braking index and the crab pulsar in @xmath1 diagram are calculated . the crab pulsar will evolve from magnetic dipole radiation dominated case towards particle wind dominated case . considering the effect of pulsar `` death '' , the crab pulsar ( and other normal pulsars ) will not evolve to the cluster of magnetars but downwards to the death valley . different acceleration models are also considered . applications to other sources are also discussed , including pulsars with braking index measured , and the magnetar population . [ firstpage ]    pulsars : general  pulsars : individal ( psr b0531 + 21 )  stars : magnetar  stars : neutron


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: large number of pions with energies below or of the order of a few hundred mev are commonly produced in high energy hadronic or nuclear collisions experiments . very soon , it has been proposed to interpret this phenomenon as resulting from classical radiation , in analogy with the classical nature of electromagnetic waves corresponding to a large number of photons . this idea has received only a marginal attention until the early @xmath0 s where it has been rediscovered and further developed , on solid theoretical grounds  @xcite . in ref . @xcite , blaizot and krzywicki studied the dynamics of a classical pion field configuration described by the non - linear sigma model , in the context of a simple idealized model of a high energy collision . they found a solution where the classical pion field oscillates coherently in a given direction in isospin space . in parallel , bjorken proposed that , in such collisions , the chiral quark condensate could be momentarily misaligned from its vacuum orientation in isospin space , resulting in a non - trivial long wavelength configuration of the pion field  @xcite . this intuitive picture corresponds to blaizot and krzywicki s solution . the `` disoriented chiral condensate '' ( dcc ) subsequently decays toward ordinary vacuum by coherent ( classical ) radiation of soft pions . the emitted pions can be though of as being in a coherent state  @xcite , schematically : r @xmath1 where  @xmath2 is the creation operator of a free pion with isospin  @xmath3 and momentum  @xmath4 and where @xmath5 is the classical source emitting pions ( summation over repeated indices is implied ) . the latter has the following form : @xmath6 that is , its orientation in isospin space is _ momentum independent_. this reflects the fact that the disorientation of the pion field in isospin space is assumed to be constant in space ( in the volume occupied by the dcc ) . the classical pion field configuration corresponding to this state is therefore a superposition of waves , which oscillate in the very same direction in isospin space . one can speak of a ( linearly ) polarized field configuration  @xcite . this property is at the origin of the most striking signature of the dcc , namely the anomalously large event - by - event fluctuations of the neutral fraction of the _ total _ number of radiated pions @xmath7 where @xmath8 is the number of pions of type @xmath9 . indeed , exploiting isospin symmetry of strong interactions , it is an easy exercise to show that the probability distribution of the above ratio is given by the inverse square root law ( see e.g.  @xcite ) : @xmath10 this is a direct manifestation of the coherence of the dcc state and is to be contrasted with the sharply peaked binomial distribution one obtains in the case of incoherent pion production . in particular , the probability that less than @xmath11 of the total number of radiated pions be neutral is expected to be as large as @xmath12 . this striking property is at the basis of all present strategies for experimental searches ( see e.g.  @xcite and references therein ) . although it had been understood that the dcc field configuration is a solution of the low energy dynamics of strong interactions , a microscopic understanding of its possible formation was still lacking . an important progress has been made in this respect in 1993 , when k.  rajagopal and f.  wilczek realized that a classical long wavelength pion field configuration could indeed be formed during the out of equilibrium chiral phase transition in heavy ion collisions : the rapid expansion of the system results in a sudden suppression of initial fluctuations ( quenching ) which in turn triggers a dramatic amplification of soft pion modes  @xcite .      to model the dynamics of the far from equilibrium chiral phase transition , rajagopal and wilczek considered the classical @xmath13 linear sigma model , with action @xmath14 the parameters @xmath15 , @xmath16 and @xmath17 are related to physical quantities via : @xmath18 m_^2 = m_^2 - 2 f_^2 = ( f_^2 - v^2 ) , h = f_m_^2 . @xmath18and are chosen so that @xmath19  mev , @xmath20  mev and @xmath21  mev . the corresponding equations of motion are evolved numerically on a lattice for a given initial field configuration . the dramatic effect of expansion is modeled by assuming an instantaneous quench from above to below the critical temperature at initial time . in practice , the values of the field @xmath22 and its first time derivative @xmath23 are chosen as independent random variables at each site of a lattice . they are sampled from gaussian distributions centered at @xmath24 , corresponding to a high temperature symmetric phase , and with variances ( in lattice units ) @xmath25 and @xmath26 , which are _ smaller _ than what they would be in the high temperature phase , thereby describing quenched fluctuations . = 7.cm    the main result of ref . @xcite is reproduced in fig . [ fig : rajwil ] , which shows the squared amplitude of different field modes as a function of time : due to the far from equilibrium initial condition , one observes a dramatic amplification of low momentum modes of the pion field at intermediate times . this phenomenon is analogous to that of domain formation after quenching a ferromagnet , or to the growth of quantum fluctuation during the preheating of the early universe after a period of inflation . it can be qualitatively understood in a mean - field approximation  @xcite : at early times , the local curvature of the effective potential ( effective mass squared ) is negative and soft modes , for which the associated effective frequency is imaginary , experience amplification . this is the so - called spinodal instability , well known in condensed matter or nuclear physics . at intermediate times , the periodic oscillations of field averaged value around the minimum of its effective potential trigger a further amplification via the mechanism of parametric resonance  @xcite .    in terms of the underlying quantum dynamics , the amplification of fourier modes amplitudes is to be interpreted as the formation of a classical configuration of the pion field . because of expansion , which eventually causes interactions to freeze out , the system may be left in such a configuration and subsequently decay through coherent pion emission . this proposal has attracted much attention and the quench scenario has been extensively studied and further developed . in particular , quantum effects in the mean field approximation have been included ( see e.g.   @xcite ) and the drastic quench approximation has been abandoned in favor of a description taking explicitely expansion into account  @xcite . in this context , the importance of initial conditions has been pointed out  @xcite . it has been realized that , in most cases , initial fluctuations are not suppressed efficiently enough and expansion hardly drive the system into a region of instability . this observation led to the conclusion that the formation of a strong pion field is a rare phenomenon . exploiting the assumption of local equilibrium at initial time , it has been possible to estimate the probability that a potentially observable signal is produced in this scenario  @xcite ( see also ref . @xcite ) . for central pb - pb collisions at cern sps energies , this analysis yields the upper bound @xmath27 where `` observable '' means that the number of produced dcc pions is required to be greater than a typical multiplicity fluctuation  @xcite . note that this estimate is obtained within the quench scenario and , therefore , corresponds to a conditional probability , which assumes that the conditions leading to a quench are met . thus the dcc is expected to be a fairly rare phenomenon , a fact which has important implications . for instance , it is very unlikely that more than one `` observable '' dcc domain be formed in a single collision . moreover , it is clear that the phenomenon should be looked for on an event - by - event basis ( see also  @xcite ) . the quench scenario has been widely accepted as a microscopic description of dcc formation in heavy ion collisions and , as such , has been extensively used for dcc phenomenology  @xcite . paradoxically , the expected large event - by - event fluctuations of the neutral ratio ( cf . eq . ( [ neutralfluct ] ) ) have never been observed in actual simulations of the out - of - equilibrium phase transition  @xcite . of course , deviations from the perfect law ( [ neutralfluct ] ) were expected on the basis of a microscopic model underlying the original picture of a perfectly polarized dcc , cf . ( [ dccstate ] ) and ( [ polarized ] ) . however , the complete absence of large fluctuations in the quench scenario suggests that this idealized picture is far from being realized in the context of the far - from - equilibrium chiral phase transition  @xcite . whether the original picture is realized or not in a realistic model is a question of great phenomenological importance and requires a careful investigation . if the rapid suppression of initial fluctuations provides a robust mechanism leading to the formation of a strong coherent pion field , it is not clear whether it can explain the hypothetical polarization ( [ polarized ] ) originally proposed . the key point is to realize that the usual assumption of a ( locally ) thermalized initial state implies that the field modes are completely uncorrelated at initial time : they are independent ( gaussian ) random numbers . therefore , in order to generate a dcc configuration , the microscopic mechanism at work needs not only to be efficient in amplifying the modes amplitudes , but it should also build correlations between amplified modes .    to study this question , we have performed a detailed statistical analysis of generic pion field configurations produced after a quench , with particular emphasis on their isospin structure  @xcite . for this purpose , it is sufficient to consider the original model of rajagopal and wilczek  @xcite , which exhibits all the relevant physical features . in particular , the classical field approximation allows one to take into account the full nonlinearities , which are crucial for the study of correlations . we sample initial field configurations in the gaussian statistical ensemble described in sec . [ sec : quench ] , and we follow the exact time evolution of each configuration by solving numerically the equations of motion corresponding to the classical action ( [ action ] ) . in this way , we can study the produced configurations on an event - by - event basis . to analyze the isospin orientations of distinct field modes , a sensitive observable is the neutral fraction of pions in each mode * k * , defined as @xmath28 where @xmath29 are the averaged occupation numbers corresponding to the classical field configuration at the time @xmath30 of measurement : @xmath31 here @xmath32 and @xmath33 denote the fourier component of the field and its time derivative respectively and @xmath34 is the free pion dispersion relation .    in the quench approximation , where the initial fluctuations are suppressed by hand , the system is initially unstable and all initial configurations undergo a dramatic amplification , as in fig . [ fig : rajwil ] . this procedure therefore selects the interesting events . to distinguish between different mechanisms responsible for amplification , we analyze the properties of the statistical ensemble at various final times @xmath30 . the plots presented here are obtained for @xmath35  fm , corresponding to the end of the spinodal instability , but the results to be discussed below are found to be insensitive to this choice . = 6.5 cm    let us concentrate on the amplified long wavelength modes , which are the one we are interested in . the event - by - event distribution of the neutral fraction ( [ ratio ] ) for the zero mode ( the most amplified ) is shown in fig . [ fig : mode ] . at a given time , all the amplified modes exhibit a similar distribution . in particular , we observe that nothing significant happens concerning their isospin structure : the neutral fraction distribution is essentially the same as the corresponding one in the initial ensemble . the latter can be computed exactly  @xcite and is represented by the dashed line in fig . [ fig : mode ] . for a given mode @xmath4 , it reads @xmath36 \ , ,   \label{indist}\ ] ] where @xmath37    \left ( \frac{\omega + 1}{\omega - ( 1 - 2f ) } \right)^{3/2 } \ , , \ ] ] and , where @xmath38 and @xmath39 are the variances of the fourier components of the field and its time derivative respectively . ] @xmath40 although not exactly  @xmath41 , the distribution is very broad , exhibiting large fluctuations around the mean value @xmath42 , which is the relevant point for phenomenology . as emphasized above , these fluctuations are a direct manifestation of the classical nature of the field modes . the small deviations from the perfect law ( [ polarized ] ) come from the fact that the latter are not strictly linearly polarized  @xcite .    however , we find that the that distinct modes have completely independent polarizations in isospin space , as can be seen on fig . [ fig : correl ] , which shows the correlation between the neutral fractions in different ( amplified ) modes . these modes oscillate coherently , but their directions of oscillation in isospin space are completely random . in other words , different modes act as independent dccs . = 6.5 cm    this has the important phenomenological consequence that the large event - by - event fluctuations of the neutral fraction are rapidly washed out when the contributions of several modes are added in a momentum bin  as commonly done in experimental analysis  even when one limits one s attention to soft modes only . this is illustrated on fig . [ fig : bin ] , where we show the neutral fraction distribution in a bin containing @xmath43 modes . the expected signal is considerably reduced , already for such a small bin . this explains the absence of large fluctuations reported in previous works  @xcite , where the authors typically considered the contribution of a large number of modes . in conclusion , contrarily to what is usually assumed , the non - linear dynamics does not build the required correlation between modes : the state produced in the simplest form of the quench scenario , where no correlations are present initially , is not the originally proposed dcc . our detailed analysis reveals that a more realistic picture is that of a superposition of waves having independent orientations in isospin space : an `` unpolarized '' dcc configuration , as depicted on fig . [ fig : dccs ] . = 6.5 cm    of course , one can not exclude the possibility that the required correlations between modes be indeed ( partially ) formed in actual nuclear or hadronic collisions by means of some other mechanism . the resulting configuration would lie between the two extreme pictures represented on fig . [ fig : dccs ] . our point here is that there exits , at present , no microscopic justification for the polarized dcc state which has been looked for so far in experiments  @xcite . instead , within our present theoretical understanding of dcc formation in heavy ion collisions , one expects an unpolarized state to be produced . it is , therefore , important to take this information into account in further theoretical as well as experimental investigations . since it has been proposed in the early 1990 s , the idea that a dcc might be formed in high energy collisions experiments has triggered an intense activity , both theoretical and experimental . attempts to detect this phenomenon include the t-864 minimax experiment , a dedicated search in proton - antiproton collisions at the fermilab tevatron  @xcite , or event - by - event analysis of pb - pb collisions by the wa98 and na49 collaborations at cern sps  @xcite . no clear experimental evidence for dcc formation has been reported so far . this is a bit disappointing in view of the fact that the phenomenon seems so natural , at least in the context of heavy ion collisions . indeed , the growth of long wavelength fluctuations after a quench is a generic phenomenon , observed in many actual situations , for example in condensed matter experiments , and which is likely to happen in heavy ion collisions as well . in particular , it has been shown to occur in various different models  @xcite . the absence of experimental evidence may indicate that the appropriate conditions for dcc formation have not been met yet ( see the interesting discussion in ref . @xcite , where it is argued that semi - peripheral heavy ion collisions at rhic energies might provide more favorable conditions for efficient quenching ) .    however , we would like to point out that the present experimental situation can actually be simply understood within our present theoretical understanding . indeed , the upper bound ( [ dccproba ] ) for the probability of dcc formation at sps energies is at the edge of the current experimental limit probed by the wa98 and na49 collaborations  @xcite . moreover , our results concerning the ( iso-)polarization of the dcc indicates that the phenomenon we are seeking may actually be more difficult to observe than originally thought . in particular , it would be of first interest to analyze existing data in the light of the present results and to design new strategies for future dcc searches at rhic and lhc  @xcite , taking into account its unpolarized nature .    besides providing a useful tool for the experimental study of the qcd chiral phase transition  @xcite , the detection of a dcc is of fundamental interest on its own . indeed , the possible formation of a coherent classical pion field is to be expected on very general grounds as a direct consequence of the bosonic nature of pions . this phenomenon is analogous to that of bose - einstein condensation , currently observed in condensed matter experiments . the fact that such a macroscopic wave of qcd matter may be produced is a very exciting possibility . i think it is worth the effort . i acknowledge a fruitful collaboration with a.  krzywicki about the present subject as well as many interesting conversations about physics and other topics . i thank j .- p . blaizot for stimulating discussions . see the review by j.  p.  blaizot and a.  krzywicki , acta phys . polon .  * b27 * ( 1996 ) 1687 and references therein . j.  p.  blaizot , a.  krzywicki , phys .  rev . * d46 * ( 1992 ) 246 ; phys .  rev .  * d50 * ( 1994 ) 442 . j.  d.  bjorken , slac - pub-5673 ; acta phys .  polon . * b23 * ( 1992 ) 637 . g.  amelino - camelia , j.  d.  bjorken , s.  e.  larsson , phys . * d56 * ( 1997 ) 6942 . d.  horn , r.  silver , annals of phys . * 66 * ( 1971 ) 509 ; j.  c.  botke , d.  j.  scalapino , r.  l.  sugar , phys . * d9 * ( 1974 ) 813 . j.  serreau , phys . rev .  * d63 * ( 2001 ) 054003 ; see also j.  serreau , phd thesis , lpt - orsay-01 - 27 , [ hep - ph/0104023 ] . b.  mohanty , d.  p.  mahapatra , t.  k.  nayak , phys . * c66 * ( 2002 ) 044901 ; b.  mohanty , t.  k.  nayak , d.  p.  mahapatra , y.  p.  viyogi , z.  huang , m.  suzuki , x.  n.  wang , phys . * d50 * ( 1994 ) 2277 ; d.  boyanovsky , h.  j.  de vega , r.  holman , s.  prem kumar , phys . * d56 * ( 1997 ) 5233 ; y.  kluger , v.  koch , j.  randrup , x.  n.  wang , phys . * c57 * ( 1998 ) 280 ; j.  i.  kapusta , s.  m.  wong , phys . * 86 * ( 2001 ) 4251 ; r.  c.  hwa , c.  b.  yang , phys . * b534 * ( 2002 ) 69 . k.  rajagopal , f.  wilczek , nucl . * b404 * ( 1993 ) 577 . s.  mrwczyski , b.  mller , phys . * b363 * ( 1995 ) 1 ; d.  i.  kaiser , phys . * d59 * ( 1999 ) 117901 . a.  dumitru , o.  scavenius , phys . * d62 * ( 2000 ) 076004 . d.  boyanovsky , h.  j.  de vega , r.  holman , phys . * d51 * ( 1995 ) 734 . f.  cooper , y.  kluger , e.  mottola , j.  p.  paz , phys . rev .  * d51 * ( 1995 ) 2377 ; m.  a.  lampert , j.  f.  dawson , f.  cooper , phys . * d54 * ( 1996 ) 2213 . a.  bialas , w.  czyz , m.  gmyrek , phys .  rev . * d51 * ( 1995 ) 3739 ; j.  randrup , phys . * 77 * ( 1996 ) 1226 . a.  krzywicki , j.  serreau , phys . * b448 * ( 1999 ) 257 . t.  s.  bir , c.  greiner , phys . * 79 * ( 1997 ) 3138 . see e.g. t.  c.  petersen , j.  randrup , phys .  rev .  * c61 * ( 2000 ) 024906 . s.  gavin , a.  gocksch , r.  d.  pisarski , phys . * 72 * ( 1994 ) 2143 ; j.  randrup , nucl . * a616 * ( 1997 ) 531 ; k.  rajagopal , talk given at hirschegg 97 , hep - ph/9703258 . g.  holzwarth , j.  klomfass , phys .  rev . * d66 * ( 2002 ) 045032 . m.  asakawa , h.  minakata , b.  mller , phys . * d58 * ( 1998 ) 094011 . minimax collaboration , t.  c.  brooks _ et al . _ , phys . rev .  * d61 * ( 2000 ) 032003 . wa98 collaboration , m.  m.  aggarwal _ et al . .  lett . * b420 * ( 1998 ) 169 ; see also nucl - ex/0206017 ; na49 collaboration , h.  appelshauser _ et al . _ , phys . * b459 * ( 1999 ) 679 . see e.g.   the event - by - event program of the star detector at rhic , ` http://www.star.bnl.gov/star ` ; searches of centauro - like events are planned at lhc e.g.   with the castor detector , see e.  gladysz - dziadus _ et al . _ , hep - ex/0209008 . c.  m.  lattes , y.  fujimoto , s.  hasegawa , phys . * 65 * ( 1980 ) 151 . f.  bedaque , a.  k.  das , mod . * a8 * ( 1993 ) 3151 ; a.  barducci , l.  caiani , r.  casalbuoni , m.  modugno , g.  pettini , r.  gatto , phys . * b369 * ( 1996 ) 23 ; j.  schaffner - bielich , j.  randrup , phys . * c59 * ( 1999 ) 3329 . k.  rajagopal , nucl . * a680 * ( 2000 ) 211 .
Abstract:

Answer: one of the main aims of present and upcoming high energy heavy ion collision experiments is to study new phases of matter at extreme temperature and density . it is expected that a nontrivial classical pion field configuration can occasionally form during the out - of - equilibrium chiral phase transition . we have recently shown that , contrarily to what has been assumed so far , this configuration is not identical to the so - called disoriented chiral condensate ( dcc ) , proposed in the early 1990 s . a detailed analysis reveals that a more realistic picture is that of an `` unpolarized '' dcc , where the fourier modes of the field have completely independent orientations in isospin space instead of being aligned with each other as in the original dcc . this has important implications concerning the possible detection of the phenomenon . in particular , the main expected signature of the original dcc , which is used in most experimental searches , is absent in the unpolarized case . we point out that the fact that no evidence of dcc formation has been reported so far in nuclear collisions actually agrees with our present theoretical understanding . new experimental strategies should be designed to look for the unpolarized dcc in existing data from sps as well as in future searches at rhic and lhc .


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 single stable isotope of beryllium , @xmath1be , is a pure product of cosmic - ray spallation of heavy ( mostly cno ) nuclei @xcite . analyses of be abundances in metal poor stars @xcite have found a relationship between [ fe / h ] = log [ n(a)/n(b)]@xmath2 - log [ n(a)/n(b)]@xmath3 and log(be / h ) with slope close to one , and between [ o / h ] and log(be / h ) with slope between 1 and 1.5 , depending on the oxygen indicator used . independently of the behaviour of the [ o / fe ] ratio at lower metallicities , these results suggest a primary production of be in the early galaxy @xcite .    as a primary element , and assuming cosmic - rays to be globally transported across the early galaxy , be may show a smaller scatter than the products of stellar nucleosynthesis @xcite at a given time , suggesting its potential use as a cosmochronometer . so far , the linear relations appear to be very tight , showing a surprisingly low scatter comparable to the measurement errors . this picture , however , might change with the increase of the samples analysed , as hinted by the results of @xcite . nevertheless , turn - off stars of the globular clusters ngc 6397 and ngc 6752 were found @xcite to have the same be abundance of field stars of the same metallicity . this strongly support the production of be to be a global process . ages derived from these abundances , in a comparison with a model of the evolution of be with time @xcite , show an excellent agreement with ages derived from theoretical isochrones , supporting the use of be as a cosmochronometer . moreover , @xcite showed that be abundances could be used to study the differences in the time scales of star formation in the halo and the thick disc of the galaxy .    in this letter , we report the discovery of an extremely be enriched halo star , hd 106038 , with an abundance 1.2 dex higher than stars of similar metallicity . this unique star deviates considerably from the observed relations of be with fe and o. it was identified during the analysis of a large sample containing near to one hundred halo and thick disc stars ( smiljanic et al . 2008 , in preparation ) . neither the standard scenario for be production , involving spallation of cosmic - rays on nuclei of the interstellar medium @xcite , nor the superbubbles ( sbs ) scenario @xcite seem to be able to produce such be enriched objects . the sbs model predict a scatter in the be abundance @xcite that may explain the stars found by @xcite and @xcite which have similar atmospheric parameters but be abundances differing by @xmath4 0.5 dex . the very high be abundance in hd 106038 , however , would require an extremely poor mixing of the sne ejecta with the ism which seems to be difficult to justify @xcite .        to the best of our knowledge , there is only one other case of extremely be enhanced star in the literature . the star j37 of the open cluster ngc 6633 was found by @xcite to have log ( be / h ) = @xmath09.0 @xmath5 0.5 . the chemical peculiarities of star j37 might be best explained by the accretion of rocky material similar to chondritic meteorites @xcite . as we shall see below , the accretion of such a material is unlikely for our population ii star . the science raw data and calibration files of hd 106038 were retrieved from the eso science archive facility . the spectra were originally obtained in 2000 april 12 with the uves spectrograph @xcite of the eso vlt at cerro paranal , chile . the data of hip 7459 ( cd@xmath061 282 ) , a halo star used as comparison in this work , were obtained in 2005 september 22 with the same instrument . the spectra have r @xmath4 45000 and a final s / n @xmath4 70 in the be region .    for both stars , we adopt the atmospheric parameters derived by @xcite . the parameters were calculated with the standard spectroscopic approach using fe i and fe ii lines ( table [ tab : par ] ) . we refer the reader to the original work for more details .    . the stars hd 106038 and hip 7459 ( filled squares ) are compared to the linear relation defined by the stars from @xcite ( open circles ) . the starred symbols are the stars from @xcite and @xcite . two of them deviate from the linear relation by @xmath4 0.50 dex.,width=245 ]    abundances ( table [ tab : par ] ) were derived through the synthesis of the spectrum around the be ii resonance lines at 3131.065  and 3130.420  . the codes for calculating synthetic spectra are described in @xcite . grids of model atmospheres without overshooting calculated by the atlas9 program @xcite were adopted . the list of atomic lines is that compiled by @xcite and the molecular line list is described in @xcite . a solar be abundance was derived using the uves solar spectrum . we estimate the total uncertainty from atmospheric parameters , continuum placement , and synthetic fit affecting the be abundance to be @xmath6 = @xmath5 0.13 dex . zero point errors might also be present due , for example , to the adopted model atmosphere . however , we are conducting a differential analysis and these should cancel out in a comparison between similar stars .    . symbols are as in fig . [ fig : boe ] . the stars of @xcite are not shown since o abundances were not derived by these authors.,width=245 ]    a comparison between the spectra of the two stars is shown in figure [ fig : be ] , in which the extreme enhancement of the be lines of hd 106038 when compared to the normal hip 7459 is clear , confirming its very high be overabundance . we show in figures [ fig : boe ] and [ fig : boe2 ] these two stars in the [ fe / h ] vs. log ( be / h ) and [ o / h ] vs. log ( be / h ) diagrams , respectivelly , together with the stars from @xcite and @xcite . of particular interest are the two stars from @xcite that deviate from the linear trend . the anomalous position of hd 106038 clearly stands out . .the adopted atmospheric parameters and beryllium abundances derived using synthetic spectra for hd 106038 , hip 7459 , a comparison star , and for the sun . [ cols="^,^,^,^,^,^ " , ]      information on other chemical abundances might help understanding the origin of the be overabundance . an overabundance of cno elements , for example , would offer a good explanation for the enhancement , since these elements are dominant in the production of be by spallation processes . @xcite determined a lithium abundance of a(@xmath7li ) = 2.48 and claimed a detection of @xmath8li , @xmath8li/@xmath7li = 0.031 , compatible with the other 8 detections out of a sample of 26 stars . the high @xmath7li abundance , however , results in a high @xmath8li abundance , a(@xmath8li ) = 0.97 , while the mean for the other detections is 0.80 . its @xmath7li is particularly remarkable since it is larger than the spite plateau @xcite . the plateau as found by @xcite is at 2.22 , which means hd 106038 has an abundance excess of @xmath9(@xmath7li ) = 2.13 . given that lithium is also expected to be produced by cosmic - ray spallation it seems safe to conclude that both @xmath7li and be overabundances have the same origin . abundances of other elements available in the literature are shown in figure [ fig : pat ] . abundances of na , ca , mg , si , ti , cr , ni , y , and ba are from @xcite . the carbon abundance is taken from @xcite and oxygen from @xcite . abundances of mn and sc are from @xcite and of s and zn from @xcite . in the same figure we also show the mean abundances of the samples analysed in these same works for comparison . in addition to be and li , the star also shows clear enhanced abundances of si , ni , and of the neutron capture elements y and ba . an enrichment in s - process elements may explain the enhanced zr line in figure [ fig : be ] . we also note the possibly larger amounts of c , s , mg , and zn , though these remain marginally compatible with the mean abundances of the samples .    , of hd 106038 determined in the literature . the abundances are represented by the full squares . in this case the error bar denotes the actual uncertainty quoted by the original work . the starred symbols indicate the mean abundance of that element for the remaining sample as found in the same work . in this case the error bar indicates the range of abundances for that given element in the original work . hd 106038 is clearly overabundant in si , ni , y , and ba , and show slightly larger abundances of c , s , mg , and zn when compared to the mean of the original samples.,width=245 ] since the standard scenario for cosmic - ray spallation does not explain the enhancement of be in hd 106038 , a peculiar and/or rare event may be related to its formation . a combination of two or more rare events to produce the observed features is unlikely , we therefore concentrate on single events .    to reproduce the very particular chemical pattern of hd 106038 , a nucleosynthetic site must be able to overproduce si and ni without overproducing other @xmath10 and iron group elements . elements in normal halo stars with the same metallicity as hd106038 come mostly from sne ii . it is therefore unlikely that the same sne ii may produce the observed enhanced [ si / fe ] and [ ni / fe ] ratios . moreover , it has about 16 times more be than what models involving sn ii predict for its metallicity @xcite . sne ia produce large amounts of fe , ni , and other iron group elements but are not expected to produce large amounts of o and si @xcite . therefore , in case of a contribution from sne ia , ratios such as [ o / fe ] and [ si / fe ] would fall below normal halo stars , contrary to what is observed . moreover , sne ia are expected to produce about one order of magnitude less spallation products than sn ii @xcite . all these features , however , may possibly be found in the ejecta of a hypernova ( hne ) , which can be enriched in both intermediate mass elements ( as s and si ) as in iron group elements @xcite . moreover , hne can produce large amounts of be and li by spallation @xcite . therefore , we suggest that the material which formed this star was probably contaminated by the nucleosynthetic products of a hne . hypernovae are core - collapse sne ( usually of type ic ) with exceptionally large kinetic energy production , resulting in spectra dominated by very broad absorption line blends @xcite . the energy released in the explosion can be one order of magnitude larger than that of normal core - collapse sne @xcite . some hypernovae , typically the most massive and energetic events , are linked to gamma - ray bursts @xcite . @xcite and @xcite calculated the yields of spallation products resulting from hne explosions . while @xcite calculate the energy distribution of the ejecta with a hydrodynamic code and solve the cosmic - ray transfer equation , @xcite use an empirical formula for the energy distribution and do not solve the transfer equation but adopt an approximation to have the mass fraction of the ejecta that produces the spallation products . @xcite claim the simplifications adopted by @xcite to overestimate the yields by a factor @xmath4 3 . the yield of be per hne can be one or two orders of magnitude larger than the one per sne ii @xcite . however , as a rare event , they are not major contributors of be in the galaxy . @xcite predict @xmath7li/@xmath1be @xmath4 8.6 and be / o @xmath4 @xmath11 ( both ratios by number ) . the calculations by @xcite predict @xmath7li/@xmath1be @xmath4 4.2 , also by number . both predictions are close to what is observed in hd 106038 ; the ratio between the observed excess of @xmath7li and @xmath1be is @xmath7li/@xmath1be = 5.6 ( while the hne is expected to have produced all the observed be abundance , it is responsible only for the excess of li with respect to the primordial plateau ) . we can not estimate the contribution of the possible hne on the observed oxygen abundance , thus only a lower limit can be placed , be / o > @xmath12 , given by the assumption that all the observed oxygen has been produced by the hne .    both models , however , predict much more @xmath8li than observed , @xmath7li/@xmath8li @xmath4 1.9 by @xcite and @xmath7li/@xmath8li @xmath4 1.2 by @xcite , by number . the observed ratio between the excess of @xmath7li and the @xmath8li abundance is @xmath7li/@xmath8li @xmath13 15 . we consider this ratio an upper limit since , given its fragility , some @xmath8li has probably been destroyed in previous evolutionary phases . the production of be without a corresponding production of @xmath8li would be extremely difficult to understand . nucleosynthetic calculations by @xcite find the ejecta of hne to have smaller amounts of c and o and larger amounts of si , s , and ar , when compared to normal sne . @xcite and @xcite also note larger [ ( zn , co)/fe ] and smaller [ ( mn , cr)/fe ] ratios . an overabundance of zn , which is not observed , can be avoided with a deeper mas cut , which would also result in a larger [ ni / fe ] @xcite . these are in qualitatively agreement with the observations , supporting the hne hypothesis . the weak s - process in massive stars seems to efficiently produce only elements with a mass number 136 . ] up to 90 @xcite . the y overabundance may require an enhanced flux of neutrons , which would also contribute to the production of ni . the ba overabundance , however , is more difficult to understand . it is not clear whether this same mechanism would result in the overproduction of ba . moreover , a significant amount of ba is expected to be produced only by the main s - process in agbs or by the r - process , usually associated to massive stars . since pollution by agbs is not possible ( see below ) , the ba overabundance is likely a product of a massive star . for example , although not expected , @xcite found ba to be overabundant by a factor of 5 in the spectrum of sn 1987a . although ba might pose a problem for our scenario , we recall that theoretical predictions for the r - process elements in hne are not available . therefore , whether this scenario is able to explain the ba abundance is still an open question .    the hne scenario , at least qualitatively , is able to explain most features observed in hd 106038 within a single peculiar event . more work , however , is still needed to show whether the scenario still holds quantitatively . a detailed comparison with nucleosynthetic predictions of theoretical models is necessary to validate or not the hne hypothesis .      in this subsection we present some alternative scenarios for the origin of the be enhancement in hd 106038 , which were discarded for the reasons presented below .    a pollution by agb stars , or any kind of evolved star , can at once be discarded . although these may explain the s - process elements , and maybe li , they do not explain the be overabundance . on the contrary , the material ejected by an agb or by a massive star , after successive mixing events , would be depleted in be . the sbs scenario @xcite would be able to reproduce the observed li / be and be / o ratios only with an extreme model where particles are accelerated from pure sne ejecta . the sbs evolution models , however , predict that most material inside of a sb comes from the ism . moreover , the remaining chemical peculiarities are not typical of sne ii ejecta and thus can not be explained within the same scenario . another possibility is the engulfing of a sub - stellar object , a planet or planetesimals debris as in star j37 of the open cluster ngc 6633 @xcite . this , however , can be excluded with a robust quantitative argument . the accreted material would be confined to the surface convective layer of the star . in a metal poor star this layer is much shallower than in a solar metallicity star . with the equation given by @xcite we estimate the surface convective layer of hd 106038 to have @xmath4 @xmath14 m@xmath15 . the mass of be in this layer is @xmath4 @xmath16 m@xmath15 while in a star with normal abundance of be it would be @xmath4 @xmath17 m@xmath15 . assuming the accreted material to have a composition similar to chondrites meteorites @xcite a mass of fe of @xmath18 m@xmath15 would also be accreted with the required mass of be . however , the total mass of fe in the convective layer of the star is @xmath4 @xmath19 m@xmath15 . in this scenario all , or almost all , fe in the convective layer of the star would come from the accreted material . hd 106038 would then originally be a population iii metal free star ; an extremely unlikely possibility . in addition , we note that the most metal poor star found to host planets has [ fe / h ] = @xmath00.68 @xcite .    a longer exposure to eps would be the natural explanation if hd 106038 , for some reason , was younger than halo stars of the same metallicity . its be abundance would be a result of the accumulated action of eps in a cloud where star formation was , somehow , delayed . its be abundance higher than the solar photospheric one could be a sign of a solar or younger age , in agreement with the suggestion of @xcite that be abundances could be used as a cosmic clock . the position of the star in the hr diagram , although not favourable for a good age determination , favours an older age and argues against this hypothesis . if the star originated in or near the galactic centre , the enhanced star formation and supernovae events could provide an enhanced eps flux . this flux might also originate from a non - stellar source such as the central black hole . a bulge origin for the star , however , seems unlikely . in particular the abundances of ni , y , and ba are not compatible with the ones of bulge stars @xcite . since it requires another source for the ni and s - process overabundances , we also discard this hypothesis . this work was developed during the visit of r.s . to eso made possible by a capes fellowship ( 1521/06 - 3 ) and support from the eso dgdf . ashwell , j. f. , jeffries , r. d. , smalley , b. , deliyannis , c. p. , steinhauer , a. , king , j. r. , 2005 , mnras , 363 , l81 asplund , m. , lambert , d. l. , nissen , p. e. , primas , f. , smith , v. v. , 2006 , apj , 644 , 229 boesgaard , a. m. , deliyannis , c. p. , king , j. r. , ryan , s. g. , vogt , s. s. , beers , t. c. , 1999 , aj , 117 , 1549 boesgaard , a. m. , novicki , m. c. , 2006 , apj , 641 , 1122 castelli , f. , kurucz , r. l. , 2003 , in proceedings of the iau symp . 210 , a20 cochran , w. d. , endl , m. , wittenmyer , r. a. , bean , j. l. , 2007 , apj , 665 , 1407 coelho , p. , barbuy , b. , melndez , j. , schiavon , r. p. , castilho , b. v. , 2005 , a&a , 443 , 735 dekker , h. , dodorico , s. , kaufer , a. , delabre , b. , kotzlowski , h. , 2000 , spie , 4008 , 534 fabbian , d. , asplund , m. , carlsson , m. , kiselman , d. , 2006 , a&a , 458 , 899 . fields , b. d. , daigne , f. , cass , m. , vangioni - flam , e. , 2002 , apj , 581 , 389 iwamoto , k. , brachwitz , f. , nomoto , k. , kishimoto , n. , umeda , h. , hix , w. r. , thielemann , f. k. , 1999 , apjs , 125 , 439 iwamoto , k. , mazzali , p. a. , nomoto , k. , umeda , h. , et al . , 1998 , nature , 395 , 672 iwamoto , k. , nomoto , k. , mazzali , p. a. , et al . , 2003 , lecture notes in physics , 598 , 243 king , j. r. , 2001 , aj , 122 , 3115 lodders , k. , 2003 , apj , 591 , 1220 mazzali , p. a. , iwamoto , k. , nomoto , k. , 2000 , apj , 545 , 407 mazzali , p. a. , lucy , l. b. , butler , k. , 1992 , a&a , 258 , 399 mcwilliam , a. , rich , r. m. , 1994 , apjs , 91 , 748 melndez , j. , shchukina , n. g. , vasiljeva , i. e. , ramrez , i , 2006 , apj , 642 , 1082 molaro , p. , bonifacio , p. , castelli , f. , pasquini , l. , 1997 , a&a , 319 , 593 murray , n. , chaboyer , b. , arras , p. , hansen , b. , noyes , r. w. , 2001 , apj , 555 , 801 nakamura , k. , shigeyama , t. , 2004 , apj , 610 , 888 nakamura , t. , umeda , h. , iwamoto , k. , nomoto , k. , hashimoto , m. , hix , w. r. , thielemann , f. k. , 2001 , apj , 555 , 880 nissen , p. e. , akerman , c. , asplund , m. , fabbian , d. , kerber , f. , 2007 , a&a , 469 , 31 nissen , p. e. , chen , y. q. , schuster , w. j. , zhao , g. , 2000 , a&a , 353 , 722 nissen , p. e. & schuster , w. j. , 1997 , a&a , 326 , 751 nomoto , k. , tominaga , n. , umeda , h. , kobayashi , c. , maeda , k. , 2006 , nuclear physics a , 777 , 424 parizot , e. , 2000 , a&a , 362 , 786 parizot , e. , drury , l. , 2000 , a&a , 356 , l66 pasquini , l. , bonifacio , p. , randich , s. , galli , d. , gratton , r. g. , 2004 , a&a , 426 , 651 pasquini , l. , bonifacio , p. , randich , s. , galli , d. , gratton , r. g. , wolff , b. , 2007 , a&a , 464 , 601 pasquini , l. , galli , d. , gratton , r. g. , bonifacio , p. , randich , s. , valle , g. , 2005 , a&a , 436 , l57 podsiadlowski , p. , nomoto , k. , maeda , k. , nakamura , t. , mazzali , p. , schmidt , b. , 2002 , apj , 567 , 491 primas , f. , duncan , d.  k. , pinsonneault , m.  h. , deliyannis , c.  p. , thorburn , j.  a. , 1997 , apj , 480 , 784 rayet , m. , hashimoto , m. , 2000 , a&a , 354 , 740 richard , o. , michaud , g. , richer , j. , 2002 , apj , 580 , 1100 reeves , h. , fowler , w. a. , hoyle , f. , 1970 , nature , 226 , 727 spite , f. , spite , m. , 1982 , a&a , 115 , 357 suzuki , t. k. , yoshii , y. , 2001 , apj , 549 , 303 valle , g. , ferrini , f. , galli , d. , shore , s. n. , 2002 , apj , 566 , 252
Abstract:

Answer: as part of a large survey of halo and thick disc stars , we found one halo star , hd 106038 , exceptionally overabundant in beryllium . in spite of its low metallicity , [ fe / h ] = @xmath01.26 , the star has log(be / h ) = @xmath010.60 , which is similar to the solar meteoritic abundance , log(be / h ) = @xmath010.58 . this abundance is more than ten times higher the abundance of stars with similar metallicity and can not be explained by models of chemical evolution of the galaxy that include the standard theory of cosmic - ray spallation . no other halo star exhibiting such a beryllium overabundance is known . in addition , overabundances of li , si , ni , y , and ba are also observed . we suggest that all these chemical peculiarities , but the ba abundance , can be simultaneously explained if the star was formed in the vicinity of a hypernova . [ firstpage ]    stars : abundances  stars : chemically peculiar  stars : individual : hd 106038