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 quasinormal modes ( qnm ) analysis is a central tool in the investigation of gravitational perturbations of stars and black holes . ( for recent comprehensive reviews of the vast literature on the subject , see @xcite and @xcite . ) the qnm analysis of black holes is particularly relevant due to the possibility of observation of the ringdown signals in gravitational wave detectors such as lisa , see @xcite for instance . besides their relevance for future observations , the qnm can also reveal important information about the structure and behavior of solutions of the einstein equations . many generalized solutions have been investigated , including non - asymptotically flat @xcite , charged ( see @xcite for a comprehensive analysis ) , and time - dependent ones @xcite . the ads / cft conjecture , in particular , has motivated many qnm analyses of generalized black hole solutions @xcite . according to the conjecture , the qnm frequencies carry information about the thermal properties of the associated conformal field theory @xcite . the vaidya metric @xcite was originally proposed to describe the spacetime outside a radiating star . it has also been the usual starting point for the study of time dependent black holes qnm @xcite . namely , it corresponds to a time dependent solution of einstein equations for a spherically symmetric body immersed in a unidirectional radial null - fluid flow . it has been also widely used in the analysis of spherically symmetric collapse and the formation of naked singularities for many years . ( for further references , see @xcite and the extensive list of @xcite ) . it is also known that the vaidya metric can be obtained from the tolman metric by taking appropriate limits in the self - similar case@xcite . this result has shed some light on the nature of the so - called shell - focusing singularities @xcite , as discussed in detail in @xcite . the vaidya metric has also proved to be useful in the study of hawking radiation and the process of black - hole evaporation @xcite , in the stochastic gravity program @xcite , and in recent numerical relativity investigations @xcite . the charged version of the metric is also an usual starting point to the study of many aspects of charged black hole physics and naked singularities @xcite . the main purpose of the present paper is to consider the qnm of time dependent backgrounds corresponding to a general vaidya metric , with special emphasis on nonstationary effects , generalizing and further developing the work done in @xcite . this kind of qnm analysis is certainly relevant from the physical point of view since , for instance , it is quite natural to expect that a real black hole be affected by processes which can change its mass such as , for instance , mass accretion or even hawking radiation , which would indeed imply a decreasing mass . any signal coming from a black hole could , in principle , have some nonstationary component .    as a secondary objective , we show that our numerical setup can be used to investigate the qnm frequencies in the case of a vaidya metric evolving towards a reisser - nordstrm naked singularity . one should expect the breakdown of our approach in the presence of a naked singularity , of course . however , we show that it is in principle possible to obtain some information about the singularity from the qnm frequencies , in agreement with the results of ishibashi and hosoya @xcite , who conclude that it is indeed possible to have regular scattering from naked singularities . this paper has three more sections . in the next one we present the main equations describing the most generic vaidya metric in double - null coordinates and write the scalar wave equation in a canonical hyperbolic form , specially suitable for the numerical analysis presented in section iii . in the last section , we present some concluding remarks . the @xmath0-dimensional vaidya metric was first discussed in @xcite . it can be easily cast in @xmath0-dimensional radiation coordinates @xmath5 as done , for instance , in @xcite . the @xmath0-dimensional charged vaidya metric in radiation coordinates , obtained originally in @xcite , reads @xmath6 where @xmath7 , @xmath8 , and @xmath9 stands for the metric of the unit @xmath10-dimensional sphere , assumed here to be spanned by the angular coordinates @xmath11 in the usual way . for the case of an ingoing radial flow , @xmath12 and @xmath1 is a monotonically increasing mass function in the advanced time @xmath13 , while @xmath14 corresponds to an outgoing radial flow , with @xmath1 being in this case a monotonically decreasing mass function in the retarded time @xmath13 . the constant @xmath15 corresponds to the total electric charge . in principle , one can also consider time dependent charges @xmath15 as done , for instance , in @xcite . this situation will of course require the presence of charged null fluids and currents , whose realistic nature we do not address here . we have already reported some preliminary results on qnm for this case in @xcite . we will now extend the approach proposed in @xcite and derive the double - null formulation for the most general vaidya metric : @xmath0-dimensional , in the presence of a cosmological constant , and with varying electric charge . only the main results are presented . the reader can get more details on the employed semi - analytical approach in @xcite and the references cited therein . we recall that the @xmath0-dimensional spherically symmetric line element in double - null coordinates @xmath16 is given by [ uv ] ds^2 = -2f(u , v)dudv + r^2(u , v)d^2_n-2 , where @xmath17 and @xmath18 are non vanishing smooth functions . the energy - momentum tensor of a unidirectional radial null - fluid in the eikonal approximation in the presence of an electromagnetic field @xmath19 is given by [ t ] t_ab = h(u , v)k_a k_b + ( f_acf^c_b -g_ab f_cdf^cd ) , where @xmath20 is a radial null vector and @xmath21 is a smooth function characterizing the null - fluid radial flow . we will consider here , without loss of generality , the case of a flow along the @xmath13-direction .    from the einstein - maxwell equations with metric ( [ uv ] ) and energy - momentum tensor ( [ t ] ) we obtain the following equations for the functions @xmath22 , @xmath23 and @xmath24 : [ b ] f = 2br_,u , [ h ] h = -2 ( ) ( m_,v - ) , @xmath25 where @xmath26 , @xmath1 and @xmath2 are arbitrary integration functions . if we choose @xmath27 constant , we can interpret @xmath1 and @xmath2 as the mass and charge of the solution , respectively . these two functions must be monotonic , and must be chosen in a way that satisfies the null - energy condition @xcite . the details of the derivation of eqs.([b])-([r2 ] ) and the full analysis of the charged vaidya metric in duble null coordinates will be presented in another paper currently under preparation . for our qnm analysis we will consider the evolution of a massless scalar field governed by the klein - gordon equation ( g^ab_,b)_,a = 0 , [ wave1 ] in the background ( [ uv ] ) . taking advantage of the spherical symmetry , we decompose the scalar field in terms of higher - dimensional spherical harmonics = _ , m(u , v)y_m(_1,  ,_n-2 ) , [ psi ] for which ^2_y_m = -(+ n -3)y_m , where @xmath28 stands for the laplacian operator over the @xmath10-dimensional unit sphere . substituting the ansatz ( [ psi ] ) in ( [ wave1 ] ) , we obtain for the metric ( [ uv ] ) ( ( r^n-2_,v)_,u + ( r^n-2_,u)_,v ) + . using now the substitution @xmath29 we get _ , uv + v(u , v)= 0 , [ phi ] where v = f - r_,ur_,v - r_,uv . [ v1 ]    the wave equation ( [ phi ] ) is already written in a canonical hyperbolic form and it will be the starting point for our qnm analysis . nevertheless , the potential ( [ v1 ] ) can still be cast in a more convenient way . from equations ( [ b ] ) and ( [ r2 ] ) , we find r_uv = -f ( - - ) , [ ruv ] and @xmath30 now we can finally write the potential @xmath31 as @xmath32 the conventions we adopted for the mass , electric charge and cosmological constant appearing in the potential ( [ potential ] ) and used in the derivations of this section are the standard ones employed in the vaidya metric literature . in particular , our expressions coincide , in the appropriate limits , with the previous results for the @xmath4 @xcite and for the @xmath33 @xcite cases . however , the commonly employed conventions in the qnm literature are slightly different , see , for instance , @xcite . the wave equation ( [ phi ] ) is already written in a canonical hyperbolic form and it can be numerically integrated by means of a characteristic scheme , _ i.e. _ , one can evolve @xmath34 in @xmath13 knowing @xmath35 . in order to determine @xmath31 , one needs to known @xmath17 and @xmath18 and , consequently , eq . ( [ r2 ] ) must be solved in each evolution step . since only @xmath34 is required to determine @xmath36 , one can implement the algorithm in an efficient way avoiding unnecessary calculations . we use here the same numerical scheme used in @xcite to evolve the system and read the qnm frequencies , which allows us to perform an exhaustive numerical analysis with modest computational resources . all mass functions considered in this work are of the form [ tanh ] 2m(v ) = ( m_f+m_i ) + ( m_f - m_i)_m ( v - v_m ) , where @xmath37 and @xmath38 stand for the initial and final mass , respectively , and @xmath39 controls how fast the change is , with the maximum rate of change at @xmath40 . we use an analogous expression for the time dependent electric charge @xmath2 . our main results do not depend on the exact form of the functions @xmath1 and @xmath2 . choices such as ( [ tanh ] ) , however , are very convenient due to their smooth behavior and the  asymptotically static " limits for @xmath41 .    for the numerical calculation of the @xmath13-evolution , the algorithm requires the evaluation of @xmath34 for all values of @xmath42 corresponding to the exterior region of the black hole , let us say , for @xmath43 , with @xmath44 corresponding to the event horizon . technically , it is easier to control the approach to the regions near the horizon by introducing an appropriated kind of tortoise coordinate @xmath45 such that the external region @xmath43 will correspond to @xmath46 . it is instructive to start by testing our code with the standard reissner - nordstrm ( rn ) case , for which a vast set of results is available in the literature @xcite . the potential ( [ potential ] ) for @xmath47 and @xmath4 becomes simply v(u , v ) = ( + - ) . [ potential_rn ] in order to compare the qnm frequencies obtained from the numerical integration of ( [ phi ] ) for the potential ( [ potential_rn ] ) with well - known results for the rn black hole , we perform the following change of variables u = u - ( - - r_+ ) - ( - - r_- ) , [ u - rn ] where r _ = m , as usual for the rn black hole . this is necessary in order to guarantee that the qnm frequencies be defined with respect to the same scales . moreover , the new tortoise coordinate @xmath45 given by ( [ u - rn ] ) allows a better description of the exterior region of the black hole , notably of the regions very close to the horizon @xmath48 . the numerical grid is , therefore , written in @xmath49 coordinates , and @xmath42 is obtained by inverting the definition ( [ u - rn ] ) .    a comparison between the the numerical results obtained with our code and the values found in the literature is shown in fig . [ fig2b ] .     obtained for a rn black hole with our code ( smooth line ) compared with values found in the literature @xcite . ] we observe a very good agreement with the known values for the rn qnm frequencies . this simple example allows us to calibrate all the algorithm parameters in order to attain a pre - established accuracy . we perform a comprehensive qnm analysis for the wave equation ( [ phi ] ) with the potential ( [ potential ] ) . we are mainly interested in the nonstationary regimes like the one described in @xcite , corresponding to the case @xmath33 , @xmath47 , and @xmath4 in the potential ( [ potential ] ) . as we show in this section , the stationary and nonstationary regimes for qnm are also present for the general situation corresponding to the potential ( [ potential ] ) . moreover , we verify appreciable deviations from the stationary regime whenever @xmath50 or @xmath51 is larger than @xmath52 , where @xmath53 is the imaginary part of the first qnm frequency . this behavior occurs irrespective of the other parameters of the potential ( [ potential ] ) , provided they be consistent with the existence of a black hole . the case of an asymptotically flat ( @xmath4 ) `` time dependent rn black hole '' , _ i.e. _ , a rn black hole with time dependent mass @xmath1 and/or charge @xmath2 functions , is particularly interesting . the determination of the new tortoise variable @xmath45 equivalent to ( [ u - rn ] ) is quite more involved for this case . the problem here is the definition , and the numerical determination , of the event horizon . ( see @xcite for a discussion of the implications of this problem in the semi - analytical approach used here . )    for our purposes here , the event horizon @xmath48 is the last null geodesic ( up to the machine precision ) escaping towards infinity , requiring the full numerical solution of ( [ r2 ] ) prior to the analysis of the wave equation ( [ phi ] ) . the determination of the real cauchy horizon @xmath54 is easier since any null geodesic inside the event horizon tends to @xmath54 along the @xmath13-evolution .    in fig . [ fig_mv ] ,          we present the time variation for the real and imaginary parts of the qnm of the scalar perturbations with @xmath55 for a 4-dimensional rn black hole with time dependent mass and charge functions . the qnm frequencies change in a similar way with the time variation of the mass or the charge , but they are typically more sensitive to the charge variations .    the results from the lower plot of fig . [ fig_mv ] can also be seen in fig . ( [ wrxwi_01 ] ) ,        in the @xmath56 plane , where we see how the variation goes from one stationary state to another through a nonstationary trajectory . in the right plot of fig.[wrxwi_01 ] , we use the minimum value of the variation of @xmath57 in the nonstationary trajectories of the left plot to quantify the deviation from the stationary trajectory . these values are plotted against @xmath58 , which can be related to a measure of how fast the variations are . we can see that the behavior can be described extremely well by a linear fit at first , deviating from this fit for @xmath59 ( faster variations ) . the case of asymptotically extremal @xmath60 black holes is specially notable in our approach due to numerical technicalities . our results for this case are presented in fig . [ wrxwi_03 ] .        in the right plot we present a convergence test , essential in this case , with the results obtained for different resolutions , while in the left plot we show how the results depend on the manner in which @xmath61 . we remind here that an extremal rn black hole is , of course , yet a regular black hole and not a naked singularity , and its qnms are well defined . it is only due to numerical limitations that most works in the literature do not reach this limit .    the left plot of fig . [ wrxwi_03 ] shows once again that the deviation from the stationary behavior increases as we consider faster variations of the background , which are quantified by @xmath58 . as this deviation becomes smaller , however , it is possible to note a nontrivial trajectory of the complex frequencies in the plane ( more easily discernible in the right plot ) , whose origin is still unclear . we can use our formalism and numerical setup to probe the formation of a naked singularity in the spacetime . the idea that it would be possible to have regular scattering from naked singularities is rather old @xcite . in fact , as it was first noticed by gibbons @xcite , minimally coupled scalar fields have the remarkable property of being regular at the origin of a reissner - nordstrm solution , where the spacetime manifold is irremediably singular . the non - singular behavior of scalar @xcite and other @xcite fields around a reissner - nordstrm naked singularity has been investigated since then . despite being mathematically well posed , the field dynamics near a naked singularity are typically ambiguous from the physical point of view since it is not clear which boundary condition one needs to impose for the field at the singularity . the main motivation of this part of our investigation is the work of ishibashi and hosoya @xcite , which demonstrated explicitly that is indeed possible to have unambiguous and regular scattering from naked singularities . after all , the presence of a naked singularity might be less harmful than originally conceived . this has been explored recently in a quantum field theory scenario in the series of papers @xcite . we apply our numerical setup to simulate a charged black hole losing mass while keeping its charge constant , leading to a reissner - nordstrm naked singularity .        for a fixed charge value @xmath15 , we set a mass function ( [ tanh ] ) between two values @xmath37 and @xmath38 , and take the final mass to be @xmath62 . at the point @xmath63 when @xmath64 the evaluation of @xmath65 is problematic since for the case of a naked singularity the full range of @xmath42 must include the singularity at @xmath66 . despite the fact that the scalar field can indeed be finite at @xmath66 @xcite , our code can not deal properly with the terms like @xmath67 in the vicinity of the singularity . this point is now under investigation , and a better numerical code is being developed for this case . our objective here is to investigate how the qnm behave until the very last moment before the formation of a naked singularity . the results obtained for the @xmath68 part of the qnms are presented in fig . ( [ mc_of_v ] ) . all simulations end at the point when an extremal rn black hole is formed , but the frequencies are all different ! this is the result of the different nonstationary trajectories followed by the qnm in each case . we remark here that in the cases with a more rapid variation @xmath69 it seems there is not enough time before the simulation stops to see the maximum deviation of the frequency and the subsequent relaxation to a possible asymptotic value . our results confirm , for the general vaidya case , the same nonstationary behavior corresponding to the inertia of the qnm frequencies identified in @xcite . in particular , for situations with @xmath70 , the qnm frequencies will not follow in the @xmath56 plane a trajectory corresponding to the instantaneous frequency associated to a rn black hole of a given @xmath71 ratio ( see figs . [ wrxwi_01 ] and [ wrxwi_03 ] ) , and the inertial behavior of the qnm can be identified . moreover , we see that the faster the change in @xmath1 or @xmath2 is , the bigger the deviation from the stationary regime will be . we could determine that the deviation from the stationary behavior is proportional to @xmath58 for masses and charges varying according to ( [ tanh ] ) , at least for values of @xmath72 , see fig . [ wrxwi_01 ] .    similarly to the uncharged case @xcite , the imaginary parts of the qnm frequencies are typically more sensitive to mass and charge changes than the real parts . however , we have identified a new behavior whose physical origin is still unclear for us : the qnm frequencies are typically more sensitive to charge than mass variations ( compare the upper and lower graphics in fig . [ fig_mv ] ) . this leads to a curious fact and a further conclusion . consider two situations such that , in the first , electric charge is constant and mass is increasing , let us say according to ( [ tanh ] ) . in the second , mass is constant and electric charge is varying in a way that both cases have the same quotient @xmath71 as function of @xmath13 . for these two situations , the instantaneous qnm frequencies will evolve differently in the @xmath73 . according to our results , the case of varying charge will typically show a bigger deviation from the stationary curve @xmath71 . this new behavior demonstrates explicitly that the ratio @xmath71 is not enough to characterize a non - stationary black hole . ( this does not represent any challenge to the no - hair theorems , though , since they typically deal with stationary solutions . )    as for the @xmath33 case considered in @xcite , an interesting extension of this work would be analysis of the highly damped qnm . since for such overtones the ratio @xmath74 is always larger than for the fundamental ( @xmath75 ) qnm considered here , including , for sufficient large @xmath0 , cases for which @xmath76 , it would be interesting to check if the stationary behavior could be improved for @xmath77 . in particular , it would be very interesting to test the high sensitivity of the instantaneous qnm overtones to electric charge variations . we remark that the numerical analysis presented here can not be directly extended to the @xmath77 case since one can not identify the overtones numerically with sufficient accuracy . we believe , however , this could be attained by means of the wkb approximation . e.  berti , v.  cardoso and a.  o.  starinets , class . * 26 * , 163001 ( 2009 ) [ arxiv:0905.2975 [ gr - qc ] ] . r.  a.  konoplya and a.  zhidenko , rev . phys . * 83 * , 793 ( 2011 ) . e. berti , v. cardoso and c.m . will , phys . rev . * d73 * , 064030 ( 2006 ) .      c. molina , d. giugno , e. abdalla and a. saa , phys . rev . * d69 * , 104013 ( 2004 ) . du , b. wang and r.k . su , phys . rev . * d70 * , 064024 ( 2004 ) ; e. abdalla , b. wang , a. lima - santos and w.g . qiu , phys . lett . * b538 * , 435 ( 2002 ) ; e. abdalla , k.h.c . castello - branco and a. lima - santos , phys . * d66 * , 104018 ( 2002 ) . e. abdalla , c.b.m.h . chirenti and a. saa , phys . rev . * d74 * , 084029 ( 2006 ) . e. abdalla , c.b.m.h . chirenti and a. saa , jhep * 0710 * , 086 ( 2007 ) . x. he , b. wang , s.f . wu and c.y . lin , phys . b673 * , 156 ( 2009 ) . m.  k.  parikh and f.  wilczek , phys . b * 449 * , 24 ( 1999 ) [ arxiv : gr - qc/9807031 ] . a.  ishibashi and a.  hosoya , phys . d * 60 * , 104028 ( 1999 ) [ arxiv : gr - qc/9907009 ] . iyer and c.v . vishveshwara , pramana * 32 * , 749 ( 1989 ) .      c.  chirenti and a.  saa , arxiv:1012.5110 [ gr - qc ] ; to apper in j. phys . conf . ser . f. girotto and a. saa , phys . rev . * d70 * , 084014 ( 2004 ) . a.  saa , phys . d * 75 * , 124019 ( 2007 ) [ arxiv : gr - qc/0701153 ] . konoplya , phys . rev . * d68 * , 124017 ( 2003 ) . v.  d.  sandberg , phys . d * 12 * , 2226 ( 1975 ) . g.  w.  gibbons , commun . * 45 * , 191 ( 1975 ) . m.  martellini , c.  reina and a.  treves , phys . d * 17 * , 2573 ( 1978 ) . f.  belgiorno , phys . d * 58 * , 084017 ( 1998 ) ; d.  batic , d.  chin and m.  nowakowski , eur . j.   c * 71 * , 1624 ( 2011 ) [ arxiv:1005.1561 [ gr - qc ] ] . j.  p.  m.  pitelli and p.  s.  letelier , phys . d * 77 * , 124030 ( 2008 ) [ arxiv:0805.3926 [ gr - qc ] ] . j.  p.  m.  pitelli and p.  s.  letelier , phys . d * 80 * , 104035 ( 2009 ) [ arxiv:0911.2626 [ gr - qc ] ] . p.  s.  letelier and j.  p.  m.  pitelli , phys . d * 82 * , 043514 ( 2010 ) [ arxiv:1008.0991 [ gr - qc ] ] . p.  s.  letelier and j.  p.  m.  pitelli , phys . d * 82 * , 104046 ( 2010 ) .
Abstract:

Answer: we consider non - stationary spherically symmetric @xmath0-dimensional charged black holes with varying mass @xmath1 and/or electric charge @xmath2 , described by generic charged vaidya metrics with cosmological constant @xmath3 in double null coordinates , and perform a comprehensive numerical analysis of the fundamental quasinormal modes ( qnm ) for minimally coupled scalar fields . we show that the `` instantaneous '' quasinormal frequencies exhibit the same sort of non - stationary behavior reported previously for the four - dimensional uncharged case with @xmath4 . such property seems to be very robust , independent of the spacetime dimension and of the metric parameters , provided they be consistent with the existence of an event horizon . the study of time dependent reissner - nordstrm black holes allows us to go a step further and quantify the deviation of the stationary regime for qnm with respect to charge variations as well . we also look for signatures in the quasinormal frequencies from the creation of a reissner - nordstrm naked spacetime singularity . even though one should expect the breakdown of our approach in the presence of naked singularities , we show that it is possible , in principle , to obtain some information about the naked singularity from the qnm frequencies , in agreement with the previous results of ishibashi and hosoya showing that it would be indeed possible to have regular scattering from naked singularities .


Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: in the risk minimisation of portfolio optimisation it is very important to consider how the returns of different companies correlate with each other . for this purpose the study of the equal time cross correlations between stocks has attracted much interest  @xcite : the clustering properties and the comparison between the time and ensemble averages have provided much useful information in this respect .    in statistical physics _ time dependent _ correlations are also of major interest . due to their role in the fluctuation dissipation theorem they constitute as the main tool for determining transport coefficients . the famous onsager reciprocity relations have their roots in the symmetry properties of the time - dependent cross correlations  @xcite . obviously , in an economic system there is no reason to assume that the time reversal symmetry or detailed balance is maintained . nevertheless , it is of interest to investigate the time dependent cross correlations between stock returns because they contain information about the way how the prices influence each other , which are the dominant stocks and to what extent this dominance is reflected in the price changes under the conditions of an efficient market . the time dependent correlations between the indices of different stock exchanges were already studied empirically by @xcite and with a microscopic model by @xcite . they showed that there exists a time shift in the cross - correlations which arises from the fact that the different stock markets are open in different time cycles during the day as the earth rotates . in this article we study time - dependent cross correlation functions of the returns of different stocks taken from the new york stock exchange ( nyse ) . as we will show , in many cases the maximum of the correlation as a function of time is not at zero but shifted , meaning that there exists some `` pulling '' effect between the companies , _ i.e. , _ one of them influences the price behaviour of the others . however , this effect can not be strong and the shift should be small , otherwise the effect could be utilized for arbitrage purposes , which is excluded from an efficient market . in fact , the investigated cases do not contradict with these criteria . the time dependent correlation between the stocks was studied before by @xcite . the results of @xcite seems to contradict with ours because they studied weakly returns , and they found significant cross - correlations on the weakly scale . however , the results agree considering the `` pulling '' effect , namely that the cross - correlations are asymmetric . the paper is organized as follows : in the next section we give a short description of the data set . in section 3 we present the method of analysis and show how it works on an artificial set of data . section 4 is devoted to the presentation of the results . the paper terminates with a discussion . one of the stylized facts of markets is that the auto - correlation of stock returns decays exponentially with a very short characteristic correlation time , which is in the range of a few minutes @xcite . this is understood as a signature of market efficiency@xcite . since cross correlations could also be used for arbitrage , one should not expect effects much beyond the above scale and therefore high frequency data are needed . we have analyzed the trade and quote ( taq ) database for @xmath0 days over the time period from 01.12.1997 to 09.03.1998 , which includes tick - by - tick data for about 10000 companies . since this is quite a short time period we selected only those companies which were traded more than 15000 times such that the number of companies reduced to 195 . having these 195 time series we have to face the following problem : since the tradings do not happen simultaneously , the values of the returns have also to be defined for the time intervals between the tradings . according to the rules of the stock exchanges we have considered the price as constant between two changes the whole trading time @xmath1 during one day is divided into @xmath2 small intervals or windows of size @xmath3 . if the trading happens in the interval @xmath4 the return takes the value @xmath5 \over \ln\left [ p(t-\delta t ) \right]}\   \nonumber , \end{aligned}\ ] ] where for simplicity the day index @xmath6 is not indicated ; otherwise it is zero .    in order to avoid the problem of major return values stemming from the differences between opening and previous day s closing prices we simply took the days as independent , _ i.e. _ , the averaging is separated into two steps : over the intraday trading time @xmath1 and over the trading days . the data prepared in this way were then analysed from the point of view of time - dependent cross correlations . as mentioned in the introduction we want to investigate the correlation of returns as the function of the time shift between pairs of stocks return time series . the definition of the time - dependent correlation function @xmath7 is @xmath8 where @xmath9 is the variance of the return . the notation @xmath10 means averaging over the whole trading time @xmath1 and important details of this process will be given in the following . since the smallest interval between two tradings is one second , then @xmath11 seems to be a natural choice . however , for such a short window it quite often happens that at a given time step there is no transaction for one of the stocks ( or for both ) such that the return results in a zero contribution to the total correlations . since the number of non - zero contributions is small , the correlation coefficients as a function of the time shift , @xmath12 , will strongly fluctuate . to avoid this problem one has to enlarge the time difference , @xmath13 , and average the correlations over the starting points of the returns . in this way the average in eq . ( [ correl_def ] ) means the following @xmath14 where the first sum runs over the starting points of the returns and the second one runs over the @xmath13 wide windows of the returns . in order to illustrate the effect that by taking larger time difference it is easier to identify the peaks in the correlation function  in other words to locate the time which gives the maximal correlation  we simulated two series of artificial data sets the first is a one dimensional persistent random walk ( prw )  @xcite , which deviates from a normal rw by the fact that the probability , @xmath15 , that it jumps in the same direction as in the previous step is higher than @xmath16 , _ i.e. _ , the random walker remembers its history . the probability of an increment @xmath17 at time @xmath4 is @xmath18 the other time process is simply generated from the first one by shifting it by @xmath19 and adding to it gaussian random noise with zero mean and width @xmath20 : @xmath21 the advantage of this model is that the correlation function can be calculated analytically and the position of the maximum correlations can be adjusted at @xmath19 : @xmath22 after generating the two data sets we randomly drop points from both sets and keep only the fraction @xmath23 of the points in order to have the same problem as with the original data sets that the jumps do not occur at the same time in the different time series . it is apparent from fig.[artificial ] . that increasing the time difference , @xmath13 helps identifying the time of maximum correlation . the fact that we dropped random points from the original data changes slightly the position of the maximum correlation as compared to eq.([anal_artific ] ) . fig.[artificial ] . shows that the decay of the correlation function is not exponential as in eq.([anal_artific ] ) but it decays approximately linearly down to the noise level . this is due to the averaging procedure we use with the increased time difference , @xmath13 . the correlation corresponding to larger time difference @xmath24 can be written as the weighted sum of the one - step correlation functions , @xmath25 , which belongs to @xmath26 : @xmath27 where @xmath28 is the return belonging to one second time difference . changing @xmath12 in eq.([correl_atiras ] ) means changing the weights of the one - step correlation functions . since the correlation function of the original data sets , see eq.([anal_artific ] ) , decays exponentially , the maximum , @xmath29 , will give the main contribution to the sum in eq.([correl_atiras ] ) and because its weight is linear in @xmath12 , then @xmath24 will decay approximately linearly . ( it should be noted that the normalization factor in eq.([correl_def ] ) does not change this consideration since it is independent of @xmath12 . ) there is only one question left namely how can we choose a smaller value for @xmath12 than for @xmath13 ? the time dependent cross correlation of the returns contains a product of the _ return _ of company a with that of company b shifted by @xmath12 . as the return is defined with the window @xmath30 the values of @xmath31 could only be multiples of @xmath13 . the solution is simply that one shifts the starting point of the return of company b by @xmath12 , as evident in eq . ( [ correl_average ] ) , i.e. we make the time shift in the price function and in this way allow any time shift larger than the minimum trading time . the above arguments of averaging give support to choose a value for @xmath13 that is larger than the minimum trading time . however , it should not be too large since the averaging leads then to the smearing out of the maximum . as the width of the one - step correlation should be a few minutes , much larger time difference would mean that in the sum of eq.([correl_atiras ] ) we mainly have terms , which are only due to noise . this suggests that the optimal choice for @xmath13 is of the order of magnitude of one minute . as mentioned in sec . [ data ] we have studied the correlation of 195 companies , which were traded during the available 54 days more than 15000 times . in accordance with the arguments presented in the previous section we have used @xmath32 but checked that the results are quite robust within the range @xmath33 . as already mentioned , we averaged over the starting points of the returns . for the maximum of the time shift we choose @xmath34 . this is definitely beyond any reasonable characteristic time for correlations in return values because of market efficiency . in fact , using such a large value for the time shift allows us to measure the noise level , which the possible effect should be compared with .    for the resulting 195 * 194/2 correlation functions we measured the maximum value , @xmath35 , the position @xmath36 at which time shift this maximum was found , and the ratio , @xmath37 , of the maximum and strength of the noise defined as the variance of the correlation values for time shift values between @xmath38 and @xmath39 . we looked at those pairs of companies for which these three values exceeded a prescribed threshold values , which we defined for @xmath40 as : @xmath41 , @xmath42 , @xmath43 . one example of the measured correlation function can be seen in fig . [ esv - xon ] . in this case the company xon ( exxon )  which is a large oil company  `` pulls '' the esv ( ensco international ) which provides drilling service to oil and gas companies . this effect is quite weak but the large value of @xmath37 shows that it is significant . the maximal value of the correlations turn out to be quite small , in average less than 0.1 , ( _ e.g. _ see fig . [ esv - xon ] ) , although the generally quoted equal time cross correlations have much larger values . the root of this effect lies in the choice of the time difference , @xmath13 . increasing @xmath13 increases the values of the equal time correlations @xcite .    in some cases the position of the maximum correlation was found at values much larger than few minutes which would be inconsistent with the efficient market behavior . a closer inspection revealed that in such cases the peak in the correlation function is caused by two major return values in the considered time series . the contribution of their product to the correlation  at appropriate value of the time shift  dominates the maximum of the correlation function . these are not the effects we are looking for , therefore we did not take them into account . in order to check whether the peak in the correlation is due to some single large return value or due to persistent influence of one of the stocks on the other we also studied how the correlation changes if the analyzed time window changes . we measured the shifted time correlation also for the first and for the second half of the given 3 month period and studied whether the correlation function remains qualitatively the same . we also measured the correlation for shorter and for larger time difference , _ i.e. _ @xmath44 and @xmath45 , respectively , because it may happen that by changing the time difference also the position of the maximal correlation value changes due to the averaging procedure described in eq.([correl_atiras ] ) . this can happen if the time dependent correlation function for @xmath13 has an asymmetric peak ; see fig . [ fre - tmx ] . let us suppose that the left hand side is higher than the right hand one . for @xmath46 the maximum will be shifted towards left as it can be shown through simple examples using eq.([correl_atiras ] ) . in the case of fig . [ esv - xon ] the correlation function is also asymmetric but not at its peak ( not near the maximum ) , which means that the maximum will not be shifted by increasing the time difference , @xmath13 . the results show that the characteristic time shift is around 100 sec . which is consistent with the effective market hypothesis . a time shift larger than the characteristic time of the decay of the return auto - correlations would contradict with the efficient market picture and could be used to arbitrage . in general the more frequently traded companies are influencing ( `` pulling '' ) the less frequently traded ones . this is not surprising since obviously the more frequently traded companies are more important . it is therefore more likely that they influence a smaller company than the other way around . although this is the generic situation , there are a few exceptions when a less often traded company `` pulls '' the other one .    in this study we found that in general one `` small '' company is influenced by many `` large '' companies and one `` large '' company pulls many `` small '' ones . as can be seen in fig . [ arrow ] this behaviour can be represented as a graph of directed links , where there are nodes from which many links go out ( meaning that this node is influenced by many others ) and there are other nodes where many links go in ( these are the big companies influencing the less important ones ) . in this paper we have analyzed the time dependent cross correlation functions of the returns of stocks at the nyse . we have studied whether there exists any pulling effect between stocks , i.e. , whether at a given time the return value of one stock influences that of another stock at a different time .    in general we can see two types of mechanisms to generate significant correlation between two stocks :    * some external effect ( e.g. economic , political news etc . ) that influences both stock prices simultaneously . in this case the change for both prices appears at the same time , and the maximum of the correlation is at zero time shift . * one of the companies has an influence to the other ( e.g. one of the company s operation depends on the other . ) . in this case the price change of the influenced stock appears later in time because it needs some time to react on the price change of the first stock , in other words one of the stocks pulls the other . this pulling effect has been the main focus of our study in this article . since the correlation between stocks was expected to be small and the available set of data was somewhat limited we had to do a careful analysis . for this reason and test purposes we generated an artificial data set with which we showed that by increasing the time window of the returns and by averaging over their starting points the detection of the correlation effect gets easier .    with the real data we saw that it is possible to find pairs of stocks where the pulling effect exists , though it turned out to be small . in addition , the characteristic time shift - given by the position of maximum correlation - was found to be of the order of a few minutes . these findings are compatible with the efficient market picture .    as for the pulling effect we found that generically the more traded , and thus more important companies pull the relatively smaller companies . this result is consistent with that of @xcite . in this light it is not surprising that in the study of the time - dependent cross correlation functions of pairs of companies in the dow jones industrial average index no pulling effect was found . this underlines the fact that the dow jones companies are indeed among the most important stocks of the new york stock exchange .    finally we would like to propose that although the observed pulling effect was small , our careful analysis could show that it is significant for a considerable set of pairs of companies . we think that this property of the stock market should be added to the so called stylized facts . of course , further analysis on more extensive data is needed to clarify further details of the time dependent cross correlations . * acknowledgment : * we are grateful to a. l. barabsi for discussions , to soon - hyung yook and hawoong jeong for their help in processing the data . support by otka ( t029985 ) and the academy of finland , project no . 1169043 ( finnish centre of excellence programme 2000 - 2005 ) are gratefully acknowledged . 99 r. n. mantegna , eur.phys . j. b. , * 11 * , 193 ( 1999 ) . l. kullmann , j. kert ' esz , r. n. mantegna , physica a , * 287 * , 412 , ( 2000 ) . l. e. reichl , _ a modern course in statistical physics _ , j. wiley and sons , new york 2nd ed . lin , r. f. engle , t. ito , rev . financial studies , * 7 * , 507 ( 1997 ) . n. vandewalle , ph . boveroux , f. brisbois , eur . j. b , * 15 * , 547 , ( 2000 ) . c. schulze , int . c , 13 , 207 ( 2002 ) . a. lo , a.c mackinlay , rev . financial studies , * 3 * , 175 ( 1990 ) . j. campbell , a. lo , a.c . mackinlay , _ the econometrics of financial markets _ , princeton university press ( 1997 ) . r. mantegna and h.e . stanley , _ econophysics _ , cambridge up , ( 2000 ) . y. liu , p. gopikrishnan , p. cizeau , m. meyer , c. peng , h. eugene stanley , phys . e , * 60 * , 1390 , ( 1999 ) . e. f. fama , j. finance , * 25 * , 383 ( 1970 ) . r. furth , ann . ( leipzig ) * 53 * , 177 ( 1917 ) ; g.h . weiss , _ aspects and applications of the random walk _ , north - holland , amsterdam , ( 1994 ) . g. bonanno , f. lillo , r. n. mantegna , quantitative finance , * 1 * , 96 , ( 2001 ) . xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx= + arc atlantic richfield co. petroleum refining + aud automatic data processing data communications and information services . + bls bellsouth corp . telephone communication + bmy bristol - myers squibb co. pharmaceutical preparations + cci citicorp banking . + cd cendant corp . travel , real estate , vehicle , and financial services + cdg cliffs drilling co. international drilling company + chv chevron corp . energy and chemical company , petroleum refining + cpq compaq computer corp . electronic computers + dis walt disney co. entertainment company + ek eastman kodak co. photography + esv ensco international inc . drilling oil and gas wells + evi energy ventures inc . oil and gas field machinery + fcn first chicago nbd corp . banking + flc falcon drilling co. inc . marine - based drilling + ge general electric co. electronics , machinery + glm global marine inc . drilling oil and gas wells + glx glaxo wellcome plc . pharmaceutical preparations + gtw gateway 2000 inc . electronic computers + hal halliburton co. oil field services + hd home depot inc . home improvement retailer + hwp hewlett - packard co. computers + ibm international business machines corp . computers + ip international paper co. paper + jnj johnson & johnson health care products + jpm morgan j.p . banking + ko coca - cola co. soft drinks + leh lehman brothers holdings financial services + mob mobil corp . petroleum refining + mot motorola inc . semiconductor technology + mrk merck & co inc . pharmaceutical preparations + mu micron technology inc . semiconductor technology + ne noble drilling corp . drilling oil and gas wells + nn newbridge networks corp . telephone and telegraph apparatus + noka nokia corp . mobile phones + pfe pfizer inc . pharmaceutical preparations + pg procter & gamble co. soap and other detergents + rd royal dutch petroleum comp . petroleum refining + sbh smithkline plc pharmaceutical preparations + slb schlumberger limited ltd oil and gas field services + sub summit bank corp . banking + tbr telecomunicacoes brasileiras s.a . telecommunications + ter teradyne inc . electrical instruments + tmx telefonos de mexico telephone communication + trv travelers group inc . fire , marine and casualty insurance + ual ual corp . air transportation + vrc varco international inc . oil and gas field services + wdc western digital corp . computer storage devices + wla warner lambert co. pharmaceutical preparations + wmt wal - mart stores inc . retail - variety stores + xon exxon corp . petroleum refining +
Abstract:

Answer: we study the time dependent cross correlations of stock returns , i.e. we measure the correlation as the function of the time shift between pairs of stock return time series using tick - by - tick data . we find a weak but significant effect showing that in many cases the maximum correlation appears at nonzero time shift indicating directions of influence between the companies . due to the weakness of this effect and the shortness of the characteristic time ( of the order of a few minutes ) our findings are compatible with market efficiency . the interaction of companies defines a directed network of influence .


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: be / x - ray binaries are x - ray sources composed of a be star and a neutron star . most of these systems are transient x - ray pulsars displaying strong outbursts in which their x - ray luminosity increases by a factor @xmath1 ( see negueruela 1998 ) . in addition , those systems in which the neutron star does not rotate fast enough for centrifugal inhibition of accretion to be effective ( see stella et al . 1986 ) display persistent x - ray luminosity at a level @xmath2 erg s@xmath3 . the high - energy radiation is believed to arise due to accretion of material associated with the be star by the compact object . it has long been known that accretion from the fast polar wind that is detected in the uv resonance lines of the be primaries can not provide the observed luminosities , even for detections at the weakest level ( see waters et al . 1988 and references therein ; see also the calculations made for x persei by telting et al . therefore it is believed that the material accreted comes from the dense equatorial disc that surrounds the be star . waters et al . ( 1989 ) modelled the radial outflow as a relatively slow ( @xmath4 km s@xmath3 ) dense wind . however most modern models for be stars consider much slower outflows , due to strong evidence for rotationally dominated motion ( quasi - keplerian discs ) . this is due not only to the line shapes ( see hanuschik et al . 1996 ) , which set an upper limit on the bulk motion at @xmath5 ( hanuschik 2000 ) , but also to the success of the global one - armed oscillation model ( which can only work in quasi - keplerian discs ) at explaining v / r variability in be stars ( hummel & hanuschik 1997 ; okazaki 1997ab ) . the viscous decretion disc model ( okazaki 1997b ; porter 1999 ; negueruela & okazaki 2000 ) considers material in quasi - keplerian orbit with an initially very subsonic outflow velocity that is gradually accelerated by gas pressure and becomes supersonic at distances @xmath6 , i.e. , much further than the orbits of neutron stars in be / x - ray transients . the transient a0535 + 26 is one of the best studied be / x - ray binaries ( clark et al . 1998 and references therein ) . it contains a slowly rotating ( @xmath7 ) neutron star in a relatively wide ( @xmath8 ) and eccentric ( @xmath9 ) orbit around the b0iiie star v725 tau ( see finger et al . 1996 ; steele et al . 1998 ) . after its last giant outburst in february 1994 ( clark et al . 1998 ; negueruela et al . 1998 ) , the strength of the emission lines in the optical spectrum of v725 tau has declined steadily . the last normal ( periodic ) outburst took place in september 1994 and the source has since not been detected by the batse experiment on board _ v725 tau , the optical counterpart to a0535 + 26 , was observed on november 7th 1998 , using the 4.2-m william herschel telescope , located at the observatorio del roque de los muchachos , la palma , spain . the telescope was equipped with the utrecht echelle spectrograph using the 31.6 lines / mm echelle centred at h@xmath0 and the site1 ccd camera . this configuration gives a resolution @xmath10 over the range @xmath11  . the data have been reduced using the _ starlink _ packages ccdpack ( draper 1998 ) , echomop ( mills et al . 1997 ) and dipso ( howarth et al . 1997 ) . a detailed analysis of the whole spectrum is left for a forthcoming paper . in figure [ fig : opt ] , we show the shape of h@xmath0 , h@xmath12 and @xmath136678 . when in emission , these three lines sample most of the radial extent of the circumstellar envelope . however , it is apparent that the lines seen in fig . [ fig : opt ] correspond to photospheric absorption from the underlying star . the emission contribution from circumstellar material , if any , is certainly very small . the asymmetry in the shape of h@xmath0 and @xmath136678suggests that some fast - moving material is present close to the stellar surface ( see hanuschik et al . 1993 ; rivinius et al . 1998 for the discussion of low - level activity in disc - less be stars ) , but the circumstellar disc is basically absent . h@xmath12 , which , when in emission , is typically produced at distances of a few @xmath14 , looks completely photospheric . in be stars h@xmath0 probes a region extending to @xmath15 , as measured from line - peak separation ( hummel & vrancken 1995 ) and direct imaging ( quirrenbach et al . 1997 ) . again , circumstellar material seems to be almost absent from this region . there is weak emission emission in - filling at the line centre  transient emission components have been seen in this star during the disc - less state ( haigh et al . 2000 ) , a behaviour typical of disc - less be stars ( rivinius et al . 1998 and references therein ) . ( 250,250 ) ( 0,0 )      observations of the source were taken using the proportional counter array ( pca ) on board _ rossixte _ on 1998 august 21 and 1998 november 12 for a total on - source time of 4170 s and 2250 s , respectively .    in both observations there is an excess of @xmath16 counts / s/(5 pcu ) in the 2.515 kev range of the standard2 data compared to the faint source background model . fits to power - law models with interstellar absorption result in flux estimates of @xmath17 and @xmath18 ( 210 kev ) for the two observations respectively . these can only be considered as upper limits on the flux , due to the uncertain contribution of diffuse galactic disc emission to the count rates ( valinia & marshall 1998 ) .      the issue of whether the source of high energy radiation was active or not during the low activity optical phase can be solved by searching for the previously reported x - ray pulsations at @xmath19103.5-s spin period ( e.g. , finger et al . 1996 ) . in order to improve the signal - to - noise , we accumulated events from the top anode layer of the detectors . we also used the latest version of the faint background model .    for the power spectral analysis we selected a stretch of continuous data and divided it into intervals of 309 s. a power spectrum was obtained for each interval and the results averaged together . given that the pulse frequency ( @xmath20 hz ) lies on a region dominated by red noise , we have to correct for such noise if the statistical significance of the pulsations are to be established . first , we fitted the poisson level by restricting ourselves to the frequency range 0.20.4 hz , that is far away from the region where the red noise component may contribute appreciably . the strongest peak in the power spectrum corresponds to @xmath21 . we also searched for periodicities in the light curves by folding the data over a period range and determining the @xmath22 of the folded light - curve ( epoch - folding technique ) . in this case we used 20 phase bins ( 19 degrees of freedom ) and a range of 100 periods , around the expected period . this method has the advantage that the result is not affected by the presence of gaps in the data , hence a longer baseline can be considered than with fourier analysis . times in the background subtracted light - curve were converted into times at the solar - system barycentre . the result for the 1998 november observation is shown in fig  [ fig : epoch ] . we found that the peak at @xmath23 is significant at @xmath24 , confirming that the source was active during the observations . ( 250,280 ) ( 0,0 )    it is worth mentioning that the detection levels shown in fig . [ fig : epoch ] were obtained without _ a priori _ knowledge of the frequency of pulsations . in other words , we searched for pulsations in the frequency / period range shown in the figure . if we take into account the fact that we are interested in the pulse period _ at _ @xmath25 , the peak becomes still more significant . the analysis of the 1998 august observation provides a much less significant detection . a peak at the expected frequency ( @xmath26 hz ) is seen in the power spectrum . however , epoch - folding analysis gave a significance of @xmath27 only when we considered one single period , that is , the number of trials is one . a search for pulsations in a period range did not yield any maximum above the 3-@xmath28 detection level although a peak at the 103-s period is present ( see fig . [ fig : epoch2 ] )    ( 250,280 ) ( 0,0 )      the pulse shape ( see figure  [ fig : pulse ] ) is nearly sinusoidal , as expected from the absence of second or higher harmonics in the power density spectra . the amplitude of the modulation is @xmath29 count s@xmath3 in the 320 kev energy range , which implies a pulse fraction of @xmath30% . given the unknown contribution from galactic disc diffuse emission , this represents only a lower limit to the pulsed fraction in the signal from the source . we have divided the november 1998 observation into two sections , corresponding to the peak of the pulse ( phase bins 0.6 to 1.0 ) and the interpulse minimum ( phase bin 0.10.5 ) . an absorbed power - law fit in the energy range 2.710 kev to the two spectra gave @xmath31 , @xmath32 , @xmath33 ( 18 dof ) for pulse maximum and @xmath34 , @xmath35 , @xmath36 ( 18 dof ) for pulse minimum . the two values are consistent with each other within the error margins . the lack of spectral changes with phase requires any significant component of the detected flux due to diffuse emission to have a spectrum similar to that of the pulsar . ( 250,280 ) ( 0,0 )      formally the x - ray spectra are equally well represented by an absorbed power - law , blackbody and bremsstrahlung models . table [ tab : models ] shows the spectral fit results . all these models gave fits of comparable quality , which means that we are unable to distinguish meaningfully between the different spectral models of table [ tab : models ] , even though the blackbody fit is unlikely to have any physical meaning , because of very small emitting area and the fact that it does not require any absorption ( introducing @xmath37 does not improve the fit )  see rutledge et al . ( 1999 ) for a discussion of the physical inadequacy of this model for neutron stars . the value of the hydrogen column density ( @xmath37 ) , which is consistent for the power - law and bremsstrahlung fits , is too high to be purely interstellar . the interstellar reddening to the source must be smaller than the measured @xmath38 ( steele et al . 1998 ) , which is the sum of interstellar and circumstellar contribution from the disc surrounding the be star . according to the relation by bohlin et al . ( 1978 ) , @xmath39 implies @xmath40 and therefore there must be a substantial contribution of local material to the absorption .    from the spectral fits , we estimate the 320 kev x - ray flux to be @xmath41 and @xmath42 for the august 1998 and november 1998 observations respectively , assuming a distance of 2 kpc ( steele et al . although the values of the spectral parameters are consistent with each other within the error margins , they all show the same trend , namely , a harder spectral state during the 1998 august observations ( lower photon index and higher blackbody and bremsstrahlung temperatures ) . lc +   + @xmath43 & 2.6@xmath440.2 + @xmath37 ( 10@xmath45 atoms @xmath46 ) & 5.6@xmath442.2 + @xmath47(dof ) & 1.46(43 ) +   + kt@xmath48 ( kev ) & 1.45@xmath440.05 + r ( km ) & 0.09@xmath440.01 + @xmath47(dof ) & 1.73(44 ) +   + @xmath49 ( kev ) & 6.4@xmath441.3 + @xmath37 ( 10@xmath45 atoms @xmath46 ) & @xmath50 + @xmath47(dof ) & 1.31(44 ) +   +   + @xmath43 & 3.2@xmath440.3 + @xmath37 ( 10@xmath45 atoms @xmath46 ) & 10.3@xmath442.6 + @xmath47(dof ) & 1.09(43 ) +   + @xmath51 ( kev ) & 1.40@xmath440.05 + @xmath52 ( km ) & 0.11@xmath440.01 + @xmath47(dof ) & 1.01(44 ) +   + @xmath49 ( kev ) & 4.4@xmath440.7 + @xmath37 ( 10@xmath45 atoms @xmath46 ) & 5.4@xmath442.1 + @xmath47(dof ) & 0.96(44 ) +    [ tab : models ]    ( 250,280 ) ( 0,0 ) our november 1998 x - ray observations represent a clear detection of a0535 + 26 at a time when the optical counterpart showed no evidence for the presence of circumstellar material . moreover , haigh et al . ( 2000 ) present spectroscopy showing that the disc was already absent as early as late august 1998 , when our first observation was taken . the observed luminosities in the 210 kev range ( @xmath53 ) are definitely smaller than the quiescence luminosity observed in other occasions when the equatorial disc surrounding the be star was present . for example , motch et al . ( 1991 ) observed the source on several occasions at a level @xmath54 in the 120 kev range ( correcting their value to the adopted distance of 2 kpc ) using _ exosat_.    the _ exosat _ observation , as well as several other quiescence detections of be / x - ray binaries with @xmath55 , have always been interpreted in terms of accretion on to the surface of the neutron star from the equatorial outflow from the be star . however , on this occasion we have observed the source when the disc of the be star had been absent for several months and at an x - ray luminosity two orders of magnitude lower . therefore we can not consider as an _ a priori _ assumption that the emission mechanism at operation is the same .    in order to explain the difference in the luminosity by two orders of magnitude , we can invoke accretion from a far less dense outflow or assume that some other emission mechanism is at work . we first consider the possibility that the observed luminosity is due to accretion on to the surface of the neutron star . assuming an efficiency @xmath56 in the conversion of gravitational energy into x - ray luminosity , @xmath57 translates into an accretion rate @xmath58 .    following stella et al . ( 1986 ) , we define the corotation radius as that at which the pulsar corotation velocity equals keplerian velocity . the corotation radius is given by @xmath59 where @xmath60 is the gravitational constant , @xmath61 is the mass of the neutron star ( assumed to be @xmath62 ) and @xmath63 is the spin period . for a0535 + 26 , @xmath64 . the magnetospheric radius at which the magnetic field begins to dominate the dynamics of the inflow depends on the accretion rate and can be expressed as @xmath65 where @xmath66 is the neutron star magnetic moment and @xmath67 is the accretion rate . @xmath68 is a constant that in the case of a0535 + 26 has been determined to be @xmath69 when an accretion disc is present ( finger et al . 1996 ) and from theoretical calculations is expected to have a similar value for wind accretion . following finger et al . ( 1996 ) , we will assume a magnetic dipolar field @xmath70 , resulting in @xmath71 .    for the accretion rate @xmath72 derived above , the magnetospheric radius would be @xmath73 . therefore @xmath74 and the neutron star must be in the centrifugal inhibition regime . in order to estimate the solidity of this result , we point out that if the 110 kev cyclotron line detected in the spectrum of a0535 + 26 ( kendziorra et al . 1994 ; grove et al . 1995 ) is the second harmonic instead of the first , as has been suggested , the magnetic field ( and magnetic moment ) would be smaller by a factor 2 . however , this would require a higher value for @xmath68 in order to fit the observations of a qpo in this system ( finger et al . 1996 ) , leaving the value of @xmath75 unaffected . an efficiency in the conversion of gravitational energy into radiation as low as @xmath76 will translate into a reduction in @xmath75 by only a factor @xmath77 . therefore we conclude that the neutron star is certain to be in the inhibited regime .    according to corbet ( 1996 ) , when the source is in the inhibition regime , a luminosity comparable to that observed could be produced by release of gravitational energy at the magnetospheric radius ( for short , accretion onto the magnetosphere . however , even assuming that the magnetosphere is at the corotation radius and an efficiency @xmath56 ( i.e. , best case ) , in order to produce @xmath57 , the accretion rate needed is @xmath78 . with such an accretion rate , the magnetosphere would be driven in well within the corotation radius , and produce a luminosity of @xmath79 by accretion on to the surface of the neutron star . therefore we conclude that the observed luminosity is not due to accretion on to the magnetosphere . therefore we are left with the following possibilities for the origin of the x - ray emission :    * accretion on to the neutron star through some sort of leakage through the magnetospheric barrier . this could adopt two forms . either directly from the be star outflow and only through a fairly limited region near the spin axis or mediated by an accretion torus , supported by the centrifugal barrier , with a small amount of material managing to penetrate the magnetosphere . * thermal emission from the heated core of the neutron star . brown et al . ( 1998 ) and rutledge et al . ( 1999 ) have studied thermal emission from x - ray transients in quiescence . they predict a thermal luminosity in quiescence @xmath80 where @xmath67 here represents the long term average mass accretion rate . from the number of observed type ii outbursts in a0535 + 26 , we assume one giant outburst every 510 years , which translates into a long - term average of @xmath81 . this would imply quiescence thermal emission in the range @xmath82 , which is consistent with our observations . it is very difficult to estimate the fraction of the signal that actually comes from the source , though the pulsed component is evidently a lower limit to it . the diffuse emission from the galactic disc is not well described at high galactic longitudes ( for a0535 + 26 , @xmath83 ) , but if the assumption by valinia & marshall ( 1998 ) that its latitude distribution should be similar to that in the galactic ridge can be held , then it should not be very strong at the position of a0535 + 26 ( @xmath84 ) . in any case , the total ( source + diffuse ) flux detected is lower than the average diffuse emission from the galactic ridge , which is @xmath85 in the 210 kev band ( valinia & marshall 1998 ) . given that the fitted spectra are much softer than the model fits to galactic ridge diffuse emission by valinia & marshall ( 1998 ) , which have photon indexes @xmath86 , and the similitude between the pulse - peak and pulse - minimum spectra , it seems likely that most of the detected signal comes actually from the source . the high value of @xmath37 obtained in all the non - thermal fits argues for the presence of large amounts of material in the vicinity of the neutron star . this could be caused by the pile - up of incoming material outside the magnetosphere . the observed spectrum is much softer than the spectra of be / x - ray binaries at low luminosity during quiescence states ( motch et al . 1991 quote a photon index of 1.7 for their observations of a0535 + 26 ) and could favour the thermal emission interpretation . brown et al . ( 1998 ) proposed that thermal reactions deep within the crust of an accreting neutron star would maintain the neutron star core at temperatures of @xmath87 , similar to that in a young radio pulsar . during the quiescent state of transient accretors , the conduction of thermal energy from the core should result in a detectable thermal spectrum from the neutron star atmosphere . in a high magnetic field pulsar , this thermal emission should be pulsed due to both the anisotropic surface temperature distribution caused by the dependence of the thermal conductivity on the magnetic field ( shibanov & yakovlev 1996 ) , and the anisotropic radiation pattern from the neutron star atmosphere resulting from the magnetic field ( zavlin et al . 1995 ) . our blackbody fit to the a0535 + 26 spectrum resulted in an emission radius much smaller than the neutron star . rutledge at al . ( 1999 ) show that fits to hydrogen atmosphere spectral models result in larger emission areas and lower effective temperature than blackbody fits . the spectra may therefore be consistent with thermal emission from the pulsar . however , if this interpretation is correct , most of the emitted luminosity should be below the energy band that we observe . in that case , the bolometric luminosity may well exceed that predicted from our estimates of the long - term average accretion rates . our detection of a0535 + 26 can also be used to set limits on the outflow from the be companion . the condition for centrifugal inhibition is @xmath88 . therefore the minimum accretion rate at which there is no inhibition corresponds to that at which @xmath89 . using the values above , we obtain an accretion rate onto the magnetosphere @xmath90 . if the observed luminosity is due to accretion on to the surface of the neutron star , the rate of mass flow into the vicinity of the neutron star is then constrained to be in the range @xmath91 ( of course , if it is due to thermal emission , only the upper limit holds ) . the lower limit represents a small fraction of the mass lost from the be star ( which should be @xmath92 ) , but the upper limit is close to the mass - loss values derived in the decretion disc model ( okazaki 1997b ; porter 1999 ) . this value is also close to the accretion rate needed to sustain the quiescence luminosity observed by motch et al . ( 1991 ) , which is @xmath93 . such an accretion rate would represent a substantial fraction of the stellar mass - loss , though still only a small fraction of the disc mass ( estimated to be @xmath94 ) . as has been calculated above , most of the long - term accretion rate is due to the giant outbursts , where a substantial fraction of the be disc mass must be accreted .    from all the above , it is clear that the amount of material reaching the vicinity of the neutron star during the disc - less phase of the companion is smaller than during previous quiescence states . this can not be due to an orbital effect because motch et al . ( 1991 ) observed the source at different orbital times and also because our two observations took place close to periastron ( orbital phases @xmath95 for the august observation and @xmath96 for the november observation , according to the ephemeris of finger et al . 1996 ) . existing evidence seems to argue against the existence of a persistent accretion disc surrounding the neutron star in a0535 + 26 . the statistical analysis of clark et al . ( 1999 ) showed that there is no significant contribution from an accretion disc to the optical / infrared luminosity of a0535 + 26 . this does not rule out the presence of an accretion disc ( which could , for example , be too small to radiate a significant flux in comparison to the be circumstellar disc ) . it is believed that a0535 + 26 does indeed form an accretion disc around the neutron during type ii outbursts , since very fast spin - up and quasi - periodic oscillations have been observed ( finger et al . 1996 ) . however , the lack of spin - up during type i outbursts led motch et al . ( 1991 ) to conclude that no persistent accretion disc was present . in contrast , the be / x - ray binary 2s1845@xmath97024 shows large spin - up during every type i outburst ( finger et al . 1999 ) . if no accretion disc is present , the reduced amount of material reaching the neutron star must be directly due to a change in the parameters of the outflow from the be star . within the framework of modern models for be star discs , considering very subsonic outflow velocities , such a change can only be due to a lower outflow density . unfortunately , since the details of the magnetic inhibition process are poorly understood , we can only constrain the mass rate reaching the neutron star to be below that corresponding to the transition at which inhibition occurs , which is very close to the rate deduced from previous quiescence observations in which the source was not in the inhibition regime . a0535 + 26 was active at a time when the optical counterpart v725 tau showed no evidence for a circumstellar disc . the luminosity was two orders of magnitude lower than in previous quiescence detections and accretion was centrifugally inhibited . given that the observed luminosity can not be due to accretion onto the magnetosphere , we are observing either some material leaking through the magnetosphere or thermal emission from the heated core of the neutron star . in any case , this detection represents a state of an accreting x - ray pulsar that had not been observed before . further observations of be / x - ray binaries in a similar state ( when their companions have lost their discs and very little material can reach the vicinity of the neutron star ) are needed . observations with _ chandra _ or _ xmm _ , which combine much higher sensitivities with more adequate energy ranges , could determine whether the observed spectrum is compatible with thermal emission models . we thank jean swank and the _ rossixte _ team for granting a target of opportunity observation and carrying it out on very short notice . simon clark is thanked for his help in preparing the proposal . the wht is operated on the island of la palma by the royal greenwich observatory in the spanish observatorio del roque de los muchachos of the instituto de astrofsica de canarias . the observations were taken as part of the ing service observing programme . this research has made use of data obtained through the high energy astrophysics science archive research center online service , provided by the nasa / goddard space flight center .
Abstract:

Answer: using the _ rossixte _ experiment , we detect weak x - ray emission from the recurrent be / x - ray transient a0535 + 26 at a time when the optical counterpart v725 tau displayed h@xmath0 in absorption , indicating the absence of a circumstellar disc . the x - ray radiation is strongly modulated at the 103.5-s pulse period of the neutron star , confirming that it originates from a0535 + 26 . the source is weaker than in previous quiescence detections by two orders of magnitude and should be in the centrifugal inhibition regime . we show that the x - ray luminosity can not be due to accretion on to the magnetosphere of the neutron star . therefore this detection represents a new state of the accreting pulsar . we speculate that the x - ray emission can be due to some matter leaking through the magnetospheric barrier or thermal radiation from the neutron star surface due to crustal heating . the observed luminosity is probably compatible with recent predictions of thermal radiation from x - ray transients in quiescence . the detection of the x - ray source in the inhibition regime implies a reduced density in the outflow from the be companion during its disc - less phase .