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: recent studies of high frequency ( hf ) data for foreign exchange ( fx ) rates @xcite regarded the many daily realizations of asset returns as constituting a statistical ensemble of histories . the single hf time series from which such ensembles can be extracted are known to present clear periodic patterns , with period of one day ( see for instance * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * and therein cited references ) . for example , in the eur / usd case , the volatility over successive 10 minutes intervals monotonically increases or decreases during specific time windows within the day . this behaviour is observed by averaging the volatility on both a daily and a weakly basis . the idea in @xcite is to endogenize the periodic , nearly deterministic behavior of the eur / usd exchange rate hf volatility into a time inhomogeneous stochastic process for the evolution of the asset . the analysis of the ensemble of daily histories produced by such a process reveals also a peculiar form of scaling obeyed by the probability density functions ( pdf ) of the returns , once they are aggregated over intervals of variable span within a fixed intraday window ( see also * ? ? ? we recall that the scaling symmetry is a relation linking the marginal pdf of a stochastic process with the duration of sampling or observation interval . in the simplest case of time intervals starting at the beginning of the daily window , the scaling property observed for eur / usd exchange rate implies that the pdf @xmath0 for a return @xmath1 over an interval of width @xmath2 , once multiplied by a factor @xmath3 , yields a specific function of @xmath4 : @xmath5 where @xmath6 is the scaling function , and @xmath7 is a suitable scaling exponent . this form leads to @xmath8 , for any arbitrary rescaling @xmath9 of the interval duration @xmath2 , and means that the pdf of the returns has a known structure , independently from the duration of sampling intervals . a process with such a property is thus called _ self - similar_. as a consequence , if we observe an ensemble of histories generated by a self - similar stochastic process and data are available for instance at a 30 minutes frequency , the scaling relation makes it possible to construct the densities of the stochastic process at both finer and time - aggregated durations . moreover , an hf self - similar process allows ( unconditional ) density forecasts for any future horizon ( say from few minutes to several hours ) . such a possibility is relevant for structural analysis of hf data , for financial forecasting , and , as we argue in the following , for the identification of efficient trading rules . the hf scaling symmetry ( [ eq_scaling ] ) addressed in @xcite has further interesting aspects : @xmath6 is non - gaussian and @xmath7 is definitely lower than @xmath10 . on the basis of the stability of the gaussian density under time aggregation , one would expect @xmath11 and a gaussian @xmath6 for independent returns . on the other hand , a scaling with @xmath12 , but a non - gaussian @xmath6 , is generally found in the empirical analysis of the pdf of aggregated returns for intervals sliding along a single , long historical time series . the deviation from gaussianity found for the pdf of returns sampled in this way is an anomalous scaling feature generally ascribed to long range dependence ( as revealed , e.g. , by the presence of volatility clustering and heteroskedasticity with persistent behaviours ) . however , for efficient markets one still has @xmath13 @xcite . in the hf eur / usd case mentioned above , besides non - gaussianity the strong deviation of @xmath7 from @xmath10 emphasizes the anomalous character of the scaling and implies the time - inhomogeneity ( non - stationarity ) of the returns process . in particular , in @xcite the time inhomogeneous scaling has been attributed to a non - markovian , strong dependence of the intraday returns , as , e.g. , manifested by the analysis of the ensemble - averaged volatility autocorrelation function . this strong dependence amounts to an effective volatility clustering phenomenon , partially hidden by the fact that the average volatility reduces in time when @xmath14 . by imposing consistency with the anomalous scaling of the aggregated return pdf , a mathematical model has been then introduced for generating the histories in the ensemble . within this model , the pdf of each return retains memory of the previous ones . once properly calibrated , the model replicates very well the statistical properties of the empirical hf ensemble in the case of eur / usd data . in particular , it is able to reproduce the patterns observed on autocovariance and autocorrelation functions of returns , absolute returns , and squared returns obtained from the hf time series .    in the present paper we apply the non - markovian model introduced in @xcite to the description of the hf data of a different financial asset , namely the s & p 500 index . a first goal is to show that the model is able to replicate the features of an asset different from fx exchange rates . more importantly , we will show how it might be used to develop an efficient intra - day trading rule , exposing a small arbitrage opportunity . this empirical application builds on the scaling property and is made under the assumption of correct model specifications . the process proposed in @xcite , conditionally to a minimum information input , i.e. the opening price and eventually the first few 10-minutes returns , can be used to provide long - range density forecasts ( for several steps ahead within a given daily time range ) . we compute density forecasts also using a simpler and more standard garch process @xcite ; namely , the asymmetric gjr model ( see * ? ? ? the density forecasts obtained from our model and those of the garch benchmark are used to define price bounds whose violation gives a trading signal . the bounds might be interpreted as predicted supports and resistances , or they could be read as a simulated price range with a given confidence interval . our approach thus belongs to the large plethora of technical analysis - based trading methods whose performances have been studies by different authors , @xcite , among others . however , the indicators used to derive the trading signals are model - based and not derived from a pure _ chartist _ approach . in - sample and out - of - sample tests demonstrate the existence of trading opportunities , leading , anyhow , to relatively small average margins of profit . as a matter of fact , those profits vary over time , reaching interesting values in specific periods and always beating the average profits associated with the garch model . the paper proceeds as follow . in subsection 1.1 we briefly describe the data used within this study , while section 2 is devoted to the presentation of the model and to parameter calibration . section 3 deals with the construction of density forecasts , and section 4 describes the trading strategy and provides the empirical results . section 5 concludes .      in this work we focus on the analysis of the s&p 500 index . we consider a dataset ranging from september 30th 1985 to october 20th 2010 . after excluding those days for which the records are not complete ( _ e.g. _ , for holidays or stock exchange anticipated closures ) , the whole dataset includes @xmath15 trading days . for each single day @xmath16 ( @xmath17 ) , we extract the index values , @xmath18 , every 10 minutes between 9:40 a.m. ( when we set @xmath19 ) and 13:00 a.m. ( @xmath20 ) , new york time . our time unit , @xmath21 minutes , is fixed by the characteristics of the dataset at our disposal . from one side it ensures that linear correlation effects among consecutive returns are negligible ; on the other side , however , it poses a bias on the effectiveness of the strategy proposed in the following , since trading actions can be executed only with a certain time - delay induced by the @xmath21-minute database discreteness . the empirical returns of the @xmath16-th day on the time - scale @xmath2 are then defined as @xmath22 , @xmath23 , @xmath24 , and are regarded as specific realizations of stochastic variables @xmath25 . in this way , our ensemble amounts to @xmath15 realizations of the stochastic process @xmath26 , with logarithmic unit increments @xmath27 . below , @xmath28 refers to a generic realization of the random variable @xmath29 . wherever appropriate , to simplify our formulas we will switch to the simpler notation @xmath30 and @xmath31 . similarly , @xmath32 will indicate the @xmath33-th return on a @xmath21-minute - scale in the @xmath16-th history realization of our ensemble . our first task is thus to identify a proper analytical model for the joint pdf of the returns @xmath34 , @xmath35 , for any @xmath23 , which correctly reproduces the statistics of the ensemble @xmath36 . in @xcite it has been suggested that hf financial time - series for the fx exchange rates offer the opportunity to deal with many realizations of the same stochastic process every day . the set of daily histories , restricted to suitable time - windows , can thus be regarded as an ensemble for testing the statistical properties of the underlying process . in fact these histories are not completely independent , because heteroskedasticity and volatility correlations , if estimated by time averages along the whole time series @xmath37 from which the ensemble @xmath38 is extracted , definitely exceed the intraday range . however , for a large enough total number of days @xmath39 in the time series one can expect inter - day correlations effects to vanish , allowing for a reliable statistics of the postulated ensemble .    the general ideas at the basis of the mathematical model adopted here to describe such an ensemble can be traced back in @xcite . financial market returns in the hf ensemble relative to the eur / usd exchange rate are linearly uncorrelated if taken over intervals of 10 minutes or more . however , they are also dependent , as shown , e.g. , by the nonzero volatility autocorrelation function . based on these evidences and on the existence of time inhomogeneous anomalous scaling properties ( see below ) , it has been proposed that the joint pdf s of these returns have the form of convex combinations of the joint pdf s of processes with independent gaussian increments of different , time dependent widths . this implies that the scaling function @xmath6 of the aggregated return pdf is also a convex combination of gaussians of variable width @xcite , as often assumed in phenomenological studies of anomalous scaling ( see * ? ? ? * and cited references ) . the main feature of the time - inhomogeneous model described in @xcite is the construction of the joint pdf @xmath40 , on the basis of the scaling symmetry valid for the pdf of the aggregated return @xmath41 the pdf of @xmath42 , @xmath43 , satisfies a scaling symmetry with scaling exponent @xmath7 and scaling function @xmath6 , if , for every @xmath33 , @xmath44 [ this equation , in more precise notations , is equivalent to eq . ( [ eq_scaling ] ) ] . under the premises of this scaling symmetry , @xmath40 is then reconstructed as @xmath45 where @xmath46^{1/2},\ ] ] @xmath47 is a pdf for a mixture of gaussian distributions with different widths @xmath48 ( see , e.g. , * ? ? ? * ; * ? ? ? * ) , and the scaling function is given by @xmath49 the coefficients @xmath50 determine the time - inhomogeneity of the variables @xmath51 s , which is also manifest in a peculiar scaling form for its pdf @xmath52 , namely , @xmath53 only for @xmath11 the variables @xmath51 s become identically distributed . see @xcite for additional details on the derivation of the joint pdf . the next step is the identification of a proper parametrization for the scaling function @xmath6 . since @xmath48 is a measure of volatility , a convenient way of representing this fat - tailed scaling function in finance ( see , @xmath54 , * ? ? ? * ) is through an _ inverse - gamma _ density @xmath55 : @xmath56 where @xmath57 and @xmath58 are a scale and a form parameter , respectively . as a simplifying assumption , we disregard the weak skewness of the return distribution , whose inclusion in the modeling framework is currently under development . the resulting scaling function @xmath59 has then a power - law decay for large @xmath60 with exponent @xmath61 , whereas @xmath62 simply sets the scale of its width . within this parametrization , @xmath6 is thus identified in terms of the two parameters @xmath62 and @xmath58 .    in this paper we extend the application of this model to the s&p 500 index , working with a larger set of realizations ( @xmath15 ) with respect to the eur / usd dataset used in @xcite . the larger database enables us to reduce the statistical uncertainties of the empirical quantities . interestingly , we find that the same general features of the dynamics of the eur / usd exchange rate also characterize the s&p 500 index and that they are both very well reproduced by the above model . in the appendix we give an account of these results , which include the observation of a systematic decrease of the empirical second moment of the elementary intraday returns ( over intervals of @xmath21 minutes ) as a function of time ( time - inhomogeneity ) , the vanishing of their linear correlation at and above the elementary time - scale , and the presence of strong non - linear correlations . in the next section we report instead the empirical evidence used to calibrate the three parameters @xmath63 , and @xmath58 . the scaling exponent @xmath7 and the scaling function @xmath6 play here a central role as they define respectively the @xmath50 s and @xmath64 appearing in the joint pdf of a given daily realization of the process . for the purposes of this study , we adopt a two step calibration procedure . at first , we determine @xmath7 , and then we use @xmath7 in order to obtain a data - collapse which allows identifying @xmath55 .    a quantitative way to calibrate @xmath7 is offered by the analysis of the moments of @xmath42 . a simple integration shows that if @xmath43 satisfies eq . ( [ eq_simple_scaling ] ) , independently of @xmath6 we have @xmath65\sim t^{qd}$ ] , @xmath66 . in fig . [ fig : nl_moments ] we plot the results of the linear regression , @xmath67 , of the logarithm of the time evolution of the empirical moments @xmath68 the resulting regression slope @xmath69 , identifies the value of @xmath7 and gives a proof of the anomalous scaling symmetry [ eq . ( [ eq_simple_scaling ] ) ] satisfied by the empirical s&p 500 data . quite remarkably , this value is very close to that estimated for fx exchange rates in refs . it is important to point out that we limit our scaling analysis to @xmath70 because of an evident multiscaling behaviour for @xmath71 ( see also * ? ? ? taking advantage of the knowledge of the scaling exponent @xmath72 , a @xmath73 minimization algorithm shows that the values @xmath74 and @xmath75 adequately fit ( see fig .  [ fig : scaling_fcn ] ) the data - collapse plot obtained via eq . ( [ eq_simple_scaling ] ) . these three parameter values , calibrated upon the entire available history , are used for in - sample evaluations of the model and subsequent trading strategies . however , we are interested in performing out - of - sample analyses also . in such cases , starting from our 25-years database , we decided to use the first 15 years of data ( from 1985 to 1999 ) to identify specific values for the parameters @xmath76 to be used to realize the trading strategy for the year 2000 . we then repeatedly shifted year by year , always using the 15 previous years to calibrate the model , until 2009 ( the last complete year ) . so , @xmath54 , we used the data from 1994 to 2008 to fit the scaling function for the strategy to be used for year 2009 . [ fig : parameters ] shows the optimal parameter values used for the implementation of the out - of - sample strategies from 2000 to 2009 . the scaling parameter @xmath7 slightly increases over time and stabilize from 2005 to 2009 . in contrast , the parameters characterizing the _ inverse - gamma _ density show evidence of a clear drop once the crisis period of 2008 is introduced in the calibration sample ( see , @xmath54 , the comparison between the _ inverse - gamma _ s of year 2005 and 2009 presented in fig .  [ fig : inverse_gammas ] ) . the model proposed here and in @xcite has implicit density forecasting capabilities . indeed , assuming correct model specifications , given the parameter values and the first daily return , it determines the pdf of the daily evolution of returns . this determination is not limited to each single return , but _ extends to any aggregation of returns _ within the validity of the scaling symmetry . such a possibility allows the construction of a trading strategy based on density forecasts . furthermore , at each point in time within the chosen scaling window , the forecasts for the remaining part might be updated using the additional intra - daily information . for this reason , we consider in the following two different trading approaches : the first makes use of the opening value of the traded asset at the beginning of the chosen time - window , @xmath77 , and then bases only on this value the density forecasts for the remaining part of it ( we call it _ unconditional trading _ ) ; in the second approach the density forecasts are updated using the information of the most recent returns ( the _ conditional trading _ , since the forecast is conditioned to the fresh intraday information ) . with both trading approaches it is possible to extract trading signals from density forecasts : given a certain quantile @xmath78 , if at a certain time within the intra - daily range under study the observed market price is above ( below ) the @xmath79 ( @xmath80 ) quantile of the predicted price density we have a buy ( sell ) signal . in this case , for each day @xmath16 , we calculate the quantiles for the expected price at time @xmath33 , conditioned to the opening value @xmath77 only . note that our approach , in computing the density forecasts , does not rely on specific previous - days asset information . it is assumed that all the information contained in the past observations has been compressed in the model parameters , namely @xmath76 . according to eq . ( [ eq_joint_probability ] ) , the pdf for the return @xmath51 , @xmath52 , is given by : @xmath81 the lower - bound values of the expected returns , @xmath82 , @xmath83 for the lower quantiles @xmath80 are then easily obtained by numerically solving the following equation with respect to @xmath82 @xmath84 due to the parity of the scaling function @xmath59 , the corresponding upper - bound values @xmath85 are simply obtained via sign flip : @xmath86 . the lower and upper expected values of @xmath26 can then be easily calculated , because @xmath26 is a monotonic function of the @xmath51 s , namely , @xmath87 and so the quantiles of @xmath26 are directly related to the quantiles of @xmath51 , once the daily opening value is given . in order to simplify our notations , we will drop the explicit dependence of the quantiles of @xmath26 on the trading day , and indicate them as @xmath88 and @xmath89 .    summarizing , for every choice of @xmath80 with @xmath78 , for every time @xmath33 from 1 to 20 within each trading day , two price barriers are obtained , @xmath88 and @xmath89 . with probability @xmath90 the price at time @xmath33 is placed between these values . for instance , in fig . [ fig : day_1 ] the results of the in - sample analysis for the first day of the dataset ( september 30th , 1985 ) are shown . as explained below in details , the comparison between these barrier values and the actual real market price lead us to the definition of a buy or sell action . the empirical analysis considers both in - sample and out - of - sample cases . the difference between the two is that for the former a unique set of values @xmath76 is used , calibrated with the whole 25-year dataset , while in the latter case the parameters @xmath76 are calibrated each year , on the basis of the previous 15-year history . for both approaches we will use the quantiles @xmath91 ( @xmath92 , respectively ) .      to exploit the non - markovian character of our model , besides the opening value of the day , @xmath77 , we can use the value of , say , the first return of the day @xmath93 to condition the subsequent expected evolution of the index . using the model prediction for the joint pdf of the returns [ eq . ( [ eq_joint_probability ] ) ] , we can build a conditioned scheme where , for example , the first return is used to condition the successive probability densities : @xmath94 more generally , we can write the following expression , conditioning to the first @xmath95 returns the pdf of the successive @xmath96 returns , with @xmath97 and @xmath98 : @xmath99 the construction of the price barriers @xmath88 and @xmath89 proceeds then as described in the previous section , simply replacing the joint pdf @xmath35 with the conditional pdf @xmath100 .    in fig . [ fig : day_2101 ] some interesting features of conditional trading signals can be detected . in particular , it is clear the influence of the nearest past returns to determine the expected amplitude of the next ones . conditional trading will thus be based on bounds which are not fixed within the day but may vary according to the new information included in the most recent market prices . when the index value @xmath101 breaks through the lower @xmath88 or the upper @xmath89 quantile , this violation can be used to define a trading strategy . note that the ensemble property we observe lasts from 9:30 a.m. to 1 p.m. as a consequence , we define a trading strategy which operates within this time lapse . to avoid the impact of news arrivals between 1 p.m. and the opening of the next day , the trading strategy opens and closes positions within the day . furthermore , since we are using a @xmath21-minute dataset , density forecasts ( and therefore the quantiles at level @xmath80 ) will be available from 9:40 a.m ..    within a certain day @xmath16 , given its specific quantiles at level @xmath80 , @xmath88 and @xmath89 , the trading signals and the trading activity are defined as follows :    * if there are no open positions * * buy if @xmath102 * * sell if @xmath103 * if there are open positions * * close a long position if @xmath104 * * close a short position if @xmath105 * * close long or short positions if they are still in place at 1 p.m.    by construction , multiple trades are possible within volatile days . differently , in trending days , single operations will take place , while during stable days , no positions will be taken . our purpose is to monitor the performances of the trading strategy defined above . therefore , we simulate the time evolution of a trader using our strategy and having an initial cash amount equal to 1 million . given the previous remarks about the validity of the ensemble properties and their link with the trading strategy , at the beginning of the day and at 1 p.m. the simulated portfolio is entirely composed by cash . furthermore , in order to avoid losses larger than the portfolio value when implementing short positions , we limit the investment value to 90% of the overall portfolio value . for symmetry , we apply the same rule also for long positions . as a result , when trades are created , 10% of the portfolio remains in cash . once a signal is observed , the trade is executed at the price @xmath101 .    as mentioned above , three different quantile levels will be considered : @xmath106 , @xmath107 , @xmath108 , and we simulate both unconditional and conditional quantiles for all these three levels . specifically , we calculate price barriers using the opening price only ( no conditioning ) , or conditioning also on the first 3 , 6 or 9 returns . higher numbers of conditioning points are not used , given that , within the day , the model proposed works only from the opening up to 1 p.m at a @xmath21-minute frequency , thus with a total of 21 observations per day . cccc & 25% & 10% & 5% +   + 0 & 2.619 & 1.602 & 0.590 + 3 & 3.788 & 2.842 & 1.862 + 6 & 3.449 & 3.027 & 2.313 + 9 & 3.243 & 3.070 & 2.184 +   + 0 & 1.599 & 1.180 & -0.483 + 3 & 3.409 & 2.748 & 1.840 + 6 & 3.094 & 2.942 & 2.564 + 9 & 2.529 & 2.652 & 1.886 +   + 0 & 3.577 & 1.977 & 1.470 + 3 & 4.173 & 2.932 & 1.882 + 6 & 3.824 & 3.110 & 2.087 + 9 & 3.996 & 3.469 & 2.455 +    the first column reports the conditioning elements ( 0 stands for no conditioning ) . the first row shows the quantile level . [ taba ]      we first analyse the trading strategy in - sample , in order to evaluate its abilities in terms of yearly profits and average return by trade ( in basis points ) . at this stage , we also verify the impact of the different number of conditioning returns . table ( [ taba ] ) reports the average return of the trades generated by the strategy during the range from october 1985 to october 2010 . profits are indicated as average basis points per trade , and distinguishing also between long and short trades . this allows to verify if the strategy better identifies signals in a specific direction . three elements clearly appear . the first one is that the profit decreases with decreasing quantile values , as if stronger signals provide smaller performances irrespective of the trade sign . we explain such an unexpected result by the fact that , due to the 10-minute discreteness of our database , trades are not valued precisely at the price at which the barrier violation occurs , but with a certain delay . this time - delay ( which could be diminished using an higher - frequency dataset ) reduces the potential profit coming from relevant ( and unexpected ) price movements . those , in fact , are the price changes detected by the model - based quantiles at the 5% level . the second relevant comment we make refers to the relation between average profit and conditioning information . we observe that conditioning the price barriers to the first returns of the day increases the profit . this is somewhat expected as during the first part of the day the model adapts its behaviour to the most recent data . and it provides therefore a better fit to the ensemble property within the remaining part of the day . again , such a result holds independently of the trade sign . regarding the number of conditioning values , we note that the relation between conditioning points and profits is not monotonic , and has a maximum between one and six conditioning points . as a consequence of this result and of the previous one about the quantile / profit relation , the out - of - sample analyses discussed in the next section have been performed with three conditioning points , which is the value with maximum profit for 25% quantiles . table ( [ tabb ] ) reports the same quantities as table ( [ taba ] ) , but by year ( we drop 1985 and 2010 where only part of the year was available ) . in general , the use of a 25% quantile provides the higher average profit . furthermore , we observe that the largest part of negative values are identified in the range 2004 - 2007 , where the market was clearly upward trending and in a low volatility phase . this is an expected outcome , since the model detects violations which are associated with high price fluctuations . the largest profits are located during 1986 and 1987 , in a high volatility period . overall , the average profit per trade is relatively small , and occasionally even negative . on the other hand , during some specific market phases the trading strategy provides large average profits per trade , see for instance the 90 s . .average profit per trade in basis points ( yearly values ) . [ cols="^,^,^,^,^,^,^,^,^,^ " , ]     [ tabf ]    overall , the results suggest the existence of some trading opportunities with respect to the time series and sample size used . furthermore , they show how the combination of advanced statistical methodologies could be used for the development of trading rules or trading schemes which have a large resemblance with those commonly used in technical analysis . in this study we realized the first practical application of a non - markovian model for asset evolution @xcite , deduced by the scaling symmetry revealed by empirical analyses . the calibration of the parameters ( @xmath7 , @xmath62 and @xmath58 in the present case ) has been facilitated by the opportunity of considering the intraday evolution of a given asset within specific time - windows , as an ensemble of different histories repeating each day @xcite . while previous studies demonstrated this property for fx exchange rates , in the present context we showed that it is possible to extend this interpretation and to verify that the model proposed in @xcite correctly reproduces the empirical evidence also for the s&p 500 hf data .    in particular , performing both in - sample and out - of - sample analyses we demonstrated how the scaling properties of the stochastic process can be used to derive long - term ( intraday ) density forecasts . in turns , these density forecasts can be used to define trading signals and implement an intra - day trading strategy which exposes a small arbitrage opportunity . by comparing the trading outcomes to those obtained from a standard garch model , namely , the gjr @xcite , we showed better performances for the trading strategy based on the proposed model . the average trade profit is limited over the entire time span , even though local levels might be higher . further studies aiming at improving the trading strategy and the empirical application of the model are thus required and under development . summarizing our findings , we can say that the proposed model has some potential for the development of trading strategies aimed at hedging the volatility risk , since their performances are positive during high market volatility , and characterized by a lower risk compared to the market index . signals extracted from the model could also be considered as confirmatory signals for other strategies working with high frequency data , or could be used to detect relevant market movements .    in this empirical example we do not consider several elements that could have an impact on the trading strategy profits . we motivate this by the need of evaluating the model in comparison to a simple benchmark . across the elements we did not include , we have the trading costs . once those are introduced , the profits reported in the previous tables would be sensibly reduced . however , the trading strategy we implement is based on a fixed frequency database , using a 10-minute interval . this has a relevant impact on the trading outcomes . in fact , if a quantile violation is observed at time @xmath33 , we execute the trade with the price observed at this same point in time . however , the violation could have taken place in any instant in the ten minutes before @xmath33 . a trader using our approach would produce quantiles to be used for each period of 10 minutes , but would immediately detect the violation , and operate in the market soon after it ( assuming she / he fully trusts the signal ) . on the contrary , working with a fixed time span of 10 minutes , we lose part of the potentially relevant content of the signal , since the price at time @xmath33 might be significantly different from the real price observed at the trade execution just after the violation occurred . another element not included in our trading example is the remuneration of the bank account . in addition , overnight liquidity operations could be introduced given that the portfolio is entirely into cash from 1 p.m. of day @xmath33 up to 9:39 a.m. of day @xmath109 . finally , we note that even the trading strategy could be improved , for instance introducing stop - loss and take - profit bounds on the implemented orders . we thank m. zamparo for useful discussions . this work is supported by `` fondazione cassa di risparmio di padova e rovigo '' within the 2008 - 2009 `` progetti di eccellenza '' program . 9    admati , a. and pfleiderer , p. , a theory of intraday patterns : volume and price variability . _ review of financia studies , _ 1988 , * 1 * , ( 1 ) , 340 .    allez , r. and bouchaud , j - p . , individual and collective stock dynamics : intra - day seasonalities . _ new journal of physics , _ 2011 , * 13 * , 25010 .    andersen , t.g . and bollerslev , t. , intraday periodicity and volatility persistence in financial markets . _ journal of empirical finance , _ 1997 , * 4 * , 115158 . andersen , t.g . and bollerslev , t. , heterogeneous information arrivals and return volatility dynamics : uncovering the long run in high volatility returns . _ journal of finance , _ 1997 , * lii 3 * , 9751005 . baldovin , f. , bovina , d. , camana , f. and stella , a.l . , modeling the non - markovian , non - stationary scaling dynamics of financial markets . in _ econophysics of order - driven markets _ ( 1st edn ) , edited by f. abergel , b. k. chakrabarti , a. chakraborti and m. mitra , part iii , pp . 239252 , 2011 ( new economic windows , springer ) . baldovin , f. and stella , a.l . , scaling and efficiency determine the irreversible evolution of a market . _ pnas , _ 2007 , * 104 * , ( 50 ) , 1974119744 . bassler , k.e . , mccauley , j.l . and gunaratne , g.h . , nonstationary increments , scaling distributions , and variable diffusion processes in financial markets . _ pnas , _ 2007 , * 104 * , ( 44 ) , 1728717290 . bollerslev , t. , generalized autoregressive conditional heteroskedasticity . _ journal of econometrics , _ 1986 , * 31 * , 307327 . bouchaud , j - p . and potters , m. , theory of financial risk and derivative pricing ( 2nd edn ) , cambridge university press , 2009 . clark , p.k . , a subordinated stochastic process model with finite variance for speculative prices . _ econometrica , _ 1973 , * 41 * , ( 1 ) , 135155 . dacorogna , m.m . , gencay , r. , muller , a.u . , olsen , r.b . , and pictet , o.v . , an introduction to high frequency finance , academic press , san diego ( 2001 ) . engle , r.f . , autoregressive conditional heteroskedasticity with estimates of the variance of u.k . inflation . _ econometrica , _ 1982 , * 50 * , 9871008 . glosten , l. , jagannathan , r. and runkle , d. , relationship between the expected value and the volatility of the nominal excess returns on stocks . _ journal of finance , _ 1993 , * 48 * , 17791801 . guillaume , d.m . , dacorogna , m.m . , dave , r.r . , muller , a.u . , olsen , r.b . , and pictet , o.v . , from the bird s eye to the microscope : a survey of new stylized facts of the intra - daily foreign exchange market , _ finance and stochastics _ 1997 , * 1 * , 95 - 129 . mantegna , r.n . and stanley , h.e . , scaling behaviour in the dynamics of an economic index . _ nature , _ 1995 , * 376 * , 4649 . mueller , u. a. , dacorogna , m. , olsen , r.b . , pictet , o.v . , schwarz , m. and morgenegg , c. , statistical study of foreign exchange rates , empirical evidence of a price change scaling law , and intraday analysis . _ journal of banking & finance , _ 1990,*14 * , ( 6 ) , 11891208 . neely , c.j . , and weller , p.a . , intraday technical trading in the foreign exchange market , journal of international money and finance , volume 22 , issue 2 , april 2003 , pages 223 - 237 . lo , a.w . , h. mamaysky , and j. wang , foundations of technical analysis : computational algorithms , statistical inference , and empirical implementation , journal of finance , 55(4 ) , 17051770 .    park , c.h , and scott , h. , what do we know about the profitability of technical analysis ? journal of economic surveys , volume 21 , issue 4 , pages 786826 , 2007 . seemann , l. , mccauley , j.l . , gunaratne g.h . intraday volatility and scaling in high frequency foreign exchange markets , international review of financial analysis volume 20 , issue 3 , june 2011 , pages 121 - 126 , 2011 . stella , a.l . , baldovin , f. , anomalous scaling due to correlations : limit theorems and self - similar processes . j. stat . p02018 ( 2010 ) . wang , f. , yamasaki , k. , havlin , s. , stanley , h.e . , indication of multiscaling in the volatility return intervals of stock markets . e 77 , 016109 ( 2008 ) . wirjanto , t.s . , and xu , d. , the applications of mixture of normal distributions in finance : a selected survey . working paper , university of waterloo , 2009 . first of all we show that the linear correlations @xmath110\over \sqrt{m_2(t,1)m_2(1,1)}}\ ] ] already decay to an approximate zero value after 10 min ( see fig .  [ fig : lin_corr ] ) . indeed , we could have chosen a 5 min scale for the last few years , but in the years before 2000 this quick decay to zero of the correlation is not verified . to avoid dealing with linearly correlated sets , we assumed a conservative attitude and decided to work at the 10-minute scale . we then observe that the empirical second moment @xmath111 ^ 2/m$ ] systematically decreases as a function of @xmath33 in the interval considered ( see fig . [ fig:2nd_moment ] ) . this is a clear indication of return non - stationarity of the underlying process at this time scale . the model reproduces such a decay , and the empirical shape is pretty well fitted . the value @xmath112 , which best fits the shape of the curve , is in good accordance with the value @xmath72 derived in the main text via analysis of non - linear moments . a widely studied indicator of non - markovianity and of presence of non - linear correlations is the volatility autocorrelation , defined as @xmath113 - \mathbb e[|r_1|]\;\mathbb e[|r_t| ] }    \over        { \mathbb e[|r_1|^2]- \mathbb e [ |r_1|]^2}},\ ] ] whose empirical counterpart is @xmath114-{{1\over m } \sum_{l=1}^{m } |r_1^l| \sum_{l'=1}^{m } |r_t^{l'}|}}\over { \sum_{l=1}^{m } |r_1^l|^2 -{{1\over m } \sum_{l=1}^{m } |r_1^l| \sum_{l'=1}^{m } |r_1^{l'}|}}}.\ ] ] fig . [ fig : volatility_autoc ] pictures again the substantial agreement between the two ( since we average over the whole ensemble , the parameters derived from the in - sample analysis have been used ) .
Abstract:

Answer: we demonstrate that a stochastic model consistent with the scaling properties of financial assets is able to replicate the empirical statistical properties of the s&p 500 high frequency data within a window of three hours in each trading day . this result extends previous findings obtained for eur / usd exchange rates . we apply the forecast capabilities of the model to implement an explicit trading strategy . trading signals are model - based and not derived from chartist criteria . in - sample and out - of - sample tests indicate that the model performs better than a benchmark asymmetric garch process , and expose the existence of small arbitrage opportunities . we discuss how to improve performances and why the trading strategy is potentially interesting to hedge volatility risk for s&p index - based products . anomalous scaling ; memory ; intraday returns ; intraday strategy .


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: black hole evaporation@xcite was one of the most surprising predictions of the field of quantum field theory in curved backgrounds . since the phenomenon has not been observed experimentally , it is of crucial importance that the assumptions underlying the prediction be examined especially carefully so as to try to understand the process as deeply as possible . one of the most unsettling features of the derivation is the dependence of the derivation on the behaviour of the fields at arbitrarily high frequencies . in fact , if we trace , in the usual derivation , the origin of the thermal emission from a solar mass sized black hole a second after the formation of the black hole , that thermal radiation has its origin in frequencies in the incoming ( vacuum ) radiation of the order of @xmath0 ( where any units make only trivial changes in the exponent ) . since @xmath1 for such frequencies correspond to masses vastly larger that the mass of the universe , it seems certain that the quantum gravitational effects of such frequencies would completely alter the behaviour of the field at such scales . would such quantum gravitational effects also destroy the thermal emission from the black hole ? since we have no quantum theory of gravity , this question is difficult to answer . however , the thermal emission is not only characteristic of black holes , it is also characteristic of dumb holes , the sonic analog of black holes@xcite@xcite . a dumb hole forms when the velocity of the fluid exceeds the velocity of sound at some closed surface . that surface forms a sonic horizon , an exact sonic analog of a black hole horizon . as i showed in 1981 , the propagation of sound waves in such a hypersonic fluid flow are exactly the same as the propagation of scalar waves in a black hole spacetime . since the derivation is not commonly known , it worthwhile repeating it here . the isentropic equation for fluid flow in a medium which is has an equation of state @xmath2 is given by ( + * v * ) = - p ( ) + ( ) = 0    i will assume that the flow is irrotational , so that we can write = the first equation can then be written as ( + 12 * vv * + h ( ) ) = 0 where @xmath3 , from which we get bernoulli s equation + 12*vv * + h ( ) = 0 ( the time dependent constant of integration can be absorbed into @xmath4 ) . we now introduce the small perturbations of @xmath4 and @xmath5 about the background flow , defining @xmath6 and @xmath7 . we get the equations + * v*+ c^2 = 0 + * v*+^2 + ln()=0    where @xmath8 . dividing the first equation by @xmath9 and adding taking the derivative with respect to @xmath10 we finally obtain ( ddt + * v * + ) c^2 ( ddt + * v * ) - 1 = 0 now , this is exactly the equation of motion for a scalar wave in the metric    @xmath11 or in four dimensions @xmath12 if the flow is stationary , we can define the killing vector ^=_0^and find the killing horizon at the point where @xmath13 . if @xmath14 is an integrable vector , we can also define @xmath15 and get the metric ds^2= ( ( c^2-v^2)d^2 - ( _ ij + v^iv^jc^2-v^2)dx^idx^j    if we furthermore assume that the flow is radial , ( in which case @xmath14 is always integrable ) the metric becomes @xmath16 where @xmath17 the usual angular metric .    by exactly the same arguments as for black holes , one expects that this hypersonic flow will produce thermal radiation with a temperature given by the  surface gravity " at the black hole horizon , which in this case , results in a temperature t= 12 dvdr_v = c    however this prediction is again made on the assumption that these sonic equations are valid for arbitrarily short wavelengths , a clearly ridiculous assumption . at wavelengths shorter than of order the intermolecular spacing , the equations for sound waves deviate significantly from the low frequency continuum fluid approximation . does this change in the fluid equations at high frequency alter the prediction of thermal emission ? if it does , one begins to loose faith in the black hole predictions . if it does not , perhaps it will give us a clue as to the real origin of the radiation .    while an exact solution of the fluid equations which fully took into account the atomic nature of the fluid would be most welcome , i am unable to supply it . however one of the most significant of the deviations from continuum fluid flow is the change in the dispersion relation for sound waves . typically at wavelengths approaching the inter atomic spacing , the group velocity of the sound waves drops to values much less than the low frequency value . can this change in the dispersion relation change the prediction of a thermal spectrum of sound waves emitted from such a dumb hole ?    thus the model i want to test is a model for sound waves within a fluid flow with a sonic horizon , or for scalar waves in a background metric , in which the dispersion relation at high frequencies differs from the usual low frequency form . in particular the model will contain a high frequency cut off , i.e. , a frequency above which there are no normal modes for the propagation of sound waves . as jacobson has argued , a naive cutoff leads to not only a suppression of the thermal emission , but also a suppression of all waves and vacuum fluctuations from the hole . in our model however the cutoff is to arises naturally in the theory .    because i am interested in the effects of changes at high frequencies to the thermal emission , not in the detailed behaviour of a real fluid , i will simplify the model so as to make the solution as easy as possible while still retaining the features responsible for the thermal emission in the naive theory . thus , to make the calculations easy , i will look at a two dimensional model ( i.e. , one space and one time ) of the propagation of waves in a flowing medium and one in which the high frequency dispersion relation is non - standard . i will also assume a constant velocity of sound ( c=1 ) , and i will assume that the background density @xmath5 is a constant . ( note that although a two dimensional background flow can not have a constant density and a varying velocity , this assumption does not change the location of the horizon , nor does it significantly alter the propagation equation of sound waves except at the very lowest frequencies . ) the model equation of motion of the scalar waves i will use in examining the effects of alteration of the high frequency behaviour of the theory is    ( _ t + _ x v ) ( _ t + v_x ) - f^2(_x ) = 0 where the function @xmath18 is such that @xmath19 is an odd analytic function of @xmath20 , @xmath21 for small @xmath22 and @xmath23 for large @xmath22 . i.e. , the dispersion function for the fluid ( taking v=0 ) is @xmath24 note that if @xmath21 for all @xmath22 , then we have the simple equation = 0 for which we can make the change of variables @xmath25 to get = 0 which has as solution ( , x ) = f(-z ) + g(+z ) where @xmath26 this is just the equation for a scalar field in a two dimensional black hole spacetime , where near the horizon , we can write ds^2=2dvdx_v=1 ( x - x_0 ) d^2 - dx^22 dvdx_v=1 ( x - x_0 )    just as for the black hole spacetime this leads to the conclusion that the temperature of the emitted radiation is given by t=12 dvdx_v=1    alternatively , one can use the analytic continuation arguments to imaginary @xmath27 time and demand regularity of the green s function at the horizon in the regular euclidean metric ds^2 = ( 1-v^2 ) d^2 + 11-v^2 dx^2 with the identification + 2= dvdx[identification ] to arrive at the same conclusion .    however , as soon as we alter the dispersion function @xmath19 , this analyticity argument fails . under the analytic continuation of @xmath28 , the equation of motion for the field becomes complex , rather than real . one can draw no conclusion from the behaviour of the green s function for this complex metric , nor does one have any reason to demand that one make the identification of the euclidean @xmath29 coordinate of eqn . [ identification ] , since the equation of motion no longer has a metric that this identification makes regular . in order to determine the consequence to the thermal emission process of the altered dispersion relation , we must go back to basic principles and calculate the bogoliubov coefficients for the conversion of ingoing positive frequency modes to outgoing negative frequency modes . i have carried out this calculation numerically for a couple of dispersion functions @xmath19 . as is usual , the particle production rate is equal to the overlap integral between a negative frequency final mode and a positive frequency initial mode . one can either start with a positive frequency initial mode , propagate this through the region containing the black or dumb hole , and calculate the negative frequency component of the final mode , or start with a negative frequency final mode and calculate the positive frequency initial component of this by propagation backwards in time . it turns out that the latter is the easier to do numerically . thus i start with some outgoing purely negative frequency wave packet . i then propagate it backwards in time toward the horizon until it is converted to an initial pulse , and calculate the positive frequency component of that initial packet . it is here necessary to state , in the case of the non - trivial dispersion relation , what i mean by positive and negative frequency . this is important because the dispersion relation implies that the notion of positive frequency does not refer to positive or negative sign of the frequencies _ per se_. the dispersion relation for constant velocity of fluid flow looks like _ = vk |f(k)| [ disp ] even if @xmath30 is less than unity , there are values of positive @xmath20 for which @xmath31 both are positive , and values of negative @xmath20 for which both are positive . in particular if @xmath32 , where @xmath33 is the solution to @xmath34 , then both possible values , @xmath31 , are positive . ( @xmath33 represents the wave vector at which the phase velocity of the wave is zero . ) thus , in the sense of the sign of @xmath35 , there are no positive frequencies at all for @xmath36 and no negative frequencies at all for @xmath37 . however , the term ` positive frequencies ' or ` negative frequencies ' refers not to the sign of @xmath35 ( although there is a correlation for free fields ) but to fixed frequency modes which have positive or negative norm . the field equations for the simplified two dimensional model for sound waves ( with constant density ) can be derived from the action i= |_t + v(t , x ) _ x |^2 - |(f(_x ) |^2 dx which leads to a momentum ( t , x)= _ t+ v(t , x)_xand an inner product = -i2(^*(t , x ) (t , x ) - ^*(t , x)(t , x ) ) dx if we assume a solution with a definite frequency and wave number @xmath38 and with constant @xmath30 , we find that = f(k ) |_,k|^2 dx where @xmath35 is given by equation [ disp ] . the definition of a positive frequency mode is that the mode have a definite frequency and that the inner product be positive , and thus that @xmath39 be positive . this selects the positive sign for @xmath40 and the negative for @xmath41 and independent of the sign of @xmath35 . one can rephrase the above discussion by stating that one is definig positive frequency by the sign of the frequency in the co - moving frame . defining @xmath42 , the mode @xmath43 becomes @xmath44 . in the case of the usual black hole situation , or the case in which the dispersion relation is linear for all frequencies for the dumb hole , we must propagate the packet backward in time until it enters the region of black ( or dumb ) hole formation . during that back propagation , the packet is constantly blue shifting , so that the frequency of the wave by the time it reaches the formation time is absurdly high . this is clearly not calculable numerically , but the process has been well understood analytically since hawking s original calculation . in the case in which the dispersion relation is such that @xmath35 becomes a constant at high frequencies , however , the calculation can be done numerically , and the positive and negative frequency overlap integral can be computed numerically . the dispersion function @xmath19 is taken to be given by f(k ) = k_0((kk_0)^n ) ^1n where @xmath45 is the spatial frequency at which the dispersion function changes from linear to constant . this gives the length scale over which the dispersion function leads to a non - locality in the propagation equations of about @xmath46 . the parameter @xmath47 determines the sharpness of the transition , with large @xmath47 implying a very fast transition between the two regimes . i have studied what happens with @xmath48 and @xmath49 in most detail . a calculation with @xmath50 was also carried out and the results were essentially the same as those for @xmath49 . i will first describe the @xmath48 case in some detail , and then briefly present the results for @xmath49 . 1the soft dispersion function ( @xmath48 ) and corresponding group velocity ( dotted line ) at zero fluid velocity . the transition wave number @xmath45 is 512 . a non - linear dispersion function means that the group velocity of the wave when the fluid velocity,@xmath30 is @xmath51 , namely v_g = dd k , is not a constant , but depends on frequency . in figure [ fig1 ] is plotted @xmath52 and @xmath53 , the group velocity , as a functions of @xmath20 for the value of @xmath48 . one can now give a hand waving description of the effect of this non - trivial dispersion function on the behaviour of a wave packet . i will describe ( and in the next section analyze ) the problem by propagating an outgoing packet backward in time . consider a long wavelength packet which is an outgoing packet at late times . it travels outward away from the sonic horizon and against the fluid flow . as we go backward in time , it approaches nearer and nearer the horizon , until when it is sufficiently close it begins to blue shift exponentially . eventually the frequency becomes sufficiently high that the group velocity of the wave is less that the velocity of the fluid at the horizon . the fluid therefore drags this outgoing packet in toward the hole ( or , since we are examining the system backward in time , the packet at earlier times is further and further from the horizon ) . the packet seems to bounce off the horizon because of the chage in the relation of the group and the fluid velocity of the packet . the naive picture is however complicated by the fact that the blue shift of the wave near the horizon is spatially very non - uniform . the front edge of the wave , the edge nearest the horizon , is blue shifted to a much larger extent ( at any given time ) than the trailing edge of the wave . in fact the leading edge is blue shifted to the final frequency long before the tail of the pulse has come near the sonic horizon . thus , it seems dangerous to simply analyze the process either in position space , or in momentum space . 2 the velocity profile of the fluid used . the system is assumed to have periodic boundary conditions . note that @xmath54 corresponds to the horizon . @xmath55 is a cauchy horizon . however , the numerical calculation seems to support this naive picture . in the numerical calculations , i have chosen my space to have periodic boundary conditions with @xmath56 identified with @xmath57 . this was done primarily because of the difficulty of specifying any other boundary conditions with the non - trivial wave equation with its unusual dispersion relations . the velocity field is chosen to be given by the function v(x)= 1- .2 ( 24(2(x-.25))+23 ) which is plotted in figure [ fig2 ] . this function , which is constant with value of .8 ( in units where the low wave number velocity of sound in a still fluid is unity ) over most of the range , was chosen precisely to have a broad range of constant value with a smooth transition to a region where @xmath58 . although the region where @xmath58 is small , this turns out to be irrelevant because the wave never penetrates that region . the horizon occurs where @xmath59 , which is where @xmath60 , or @xmath61 . it is the lower value ( @xmath54 ) which represents the horizon , while the other represents an irrelevant cauchy horizon . another value of interest is the derivative of velocity with respect to position , @xmath62 at the horizon , since this gives the temperature . this is given by _ v=1 = .2(48(1-(23/24)^2)=.4(47 ) 8.6 giving a temperature of t=12dvdx= 1.37 . ( recall that our units are such that @xmath63 , so that the thermal wave number is @xmath64 and the thermal wave length is thus approximately 8 , which is longer than the whole spatial domain . the final ( late in time ) wave packets ( recall that we will be propagating the packets backward in time ) , chosen to be purely negative frequency packets , are thus complex , and are traveling away from the horizon . they are confined as much as is possible to the region where the velocity is constant . they now travel ( backwards in time ) toward the horizon of the dumb hole , where they begin to blue shift . eventually their frequency becomes sufficiently high that the group velocity drops to less than unity and they are dragged backward out of the region of the hole to the region where the fluid velocity is constant . in this region we can again analyze the wave into its positive and negative frequency parts , and determine the bogoliubov transformation quantity for this particular wave packet . this analysis must be carried out in the constant velocity regions . whether or not a wave is `` positive frequency '' depends on its time dependent behaviour . however , both because of the finite size of the system , and because of the finite time for which one can reliably carry out the integration , the temporal behaviour is difficult to use in determining whether or not the wave is positive frequency . however , in regions where the fluid velocity is constant , one can directly compare the relation between the field and its first time derivative to determine its positive and negative frequency components . we have ( _ t - v _ x)= + ( _ t - v_x)= f^2(_x)fourier transforming under the assumption that @xmath30 is a constant , we have ( _ -vk ) _ ( k)=_(k ) + ( _ -vk)_(k)= f^2(k ) _ ( k ) we can solve for @xmath35 and choose that _ = vk |f(k)| with @xmath65 the positive frequency branch for all values of @xmath20 ( independent of the sign of @xmath65 ) . now , given @xmath66 and @xmath67 we can split it into its positive and negative frequency parts by @xmath68 and @xmath69 , and the @xmath70 subscript denotes the positive and negative frequency component . solving for the positive frequency part , we get _ + ( k)= ( ( -i(k)/|f(k)| + ( k ) ) + _ + ( k)= ( ( ( k ) -i|f(k)|(k ) )    the procedure now is to set up the initial ( late time ) wave as a purely negative frequency wave , propagate it backward in time until the resultant wave is again entirely ( or as much as is possible ) in the constant velocity regime . the resultant wave is then analyzed into its positive and negative frequency parts , and the inner product of each of these parts separately is calculated . note that the inner product of the full wave itself is conserved , and this can be used as a test of the integration routine . we finally get the actual bogoliubov coefficients for the wave from the future to the past ( and thus also from the past to the future ) and can compare this with the thermal hypothesis . one can also calculate the expected bogoliubov coefficient under the assumption that the outgoing mode is thermally distributed , as would be expected from the naive application of the hawking procedure to the wave . the wave packets used are not plane wave states , and thus one must integrate the thermal factor over the spectrum of the outgoing wave to find the expected amplitude of the thermal bogoliubov coefficient . it this thermal factor equals the bogoliubov coefficient calculated from the evolution equations , then the thermal radiation for dumb holes hypothesis is validated even with this non standard high frequency dispersion function . 3 the real ( solid line ) and imaginary ( dotted line ) parts of a final pulse emerging from the dumb hole . the dispersion function has @xmath48 and @xmath71 .    4 an intermediate stage in the evolution backward in time of the real part of the pulse of figure [ fig3 ] . note that the front of the wave has crashed into the horizon , and been blue shifted to a high frequency pulse whose group velocity is less than the fluid velocity . the dotted line is the same calculation but at half the spatial resolution as the solid line .    5 the pulse as it must have been initially to produce the outgoing pulse of figure [ fig3 ] . note that it is an outgoing pulse but is being dragged in toward the horizon by the inflowing fluid . again the dotted line is at half the resolution as the solid line to give an estimate of the numerical errors . figures 3 gives the real ( solid line ) and imaginary ( dotted line ) parts of the wave packet at the final time . in figures [ fig4 ] and [ fig5 ] i show the real part of the wave at two separate earlier times . in figure [ fig4 ] the wave is just interacting with the horizon . part of the wave has been blue shifted ( going backwards in time ) and part is still at the final low frequency . figure 6 gives the incoming wave in the constant velocity region as it must have been in order to produce the outgoing wave of figure [ fig3 ] . we note in figure [ fig4 ] that the front edge of the wave has already crashed into the  horizon " ( @xmath72 ) , been blue - shifted and begun to be dragged back out of the dumb hole by the inflowing ( which in the reverse time direction is outflowing ) fluid because its group velocity has decreased at the higher frequencies . at the same time the trailing edge of the pulse is still traveling toward the hole .    in figures [ fig4 ] and [ fig5 ] i also show the outcome of the numerical evolution for two different grid spacings ( the dotted is at half the resolution of the solid line ) to give a feeling for the accuracy of the integration . the key difference is in the phase velocity of the wave being dragged out of the hole . this is expected , because of the finite differencing of the term @xmath73 term in the equations . the centered differencing scheme converts this term to @xmath74 instead of the exact @xmath75 . this means that at large @xmath20 , the effective velocity of the fluid is decreased by @xmath76 . i.e. , at high frequencies , the effective fluid velocity is less leading to a slightly smaller dragging velocity for the high frequency wave , as is seen in the simulation . 6 the spectrum of the final pulse of figure [ fig3 ] . note that it is purely positive wave number and negative frequency . 7 the spectrum of the initial pulse of figure [ fig5 ] . the solid line is positive frequency components while the dotted line is negative frequencies . since positive frequencies correspond to positive wave number and vice versa , this still represents a purely outgoing wave .    in figure [ fig6 ] i give the spectrum ( in k space ) of the final outgoing wave packet . this wave has been designed to be purely negative frequency , left traveling and concentrated in the region where the velocity is constant . it represents a wave packet traveling away from the hole against the inflowing fluid flow . in figure [ fig7 ] we have the spectrum corresponding to figure [ fig6 ] , divided into positive and negative frequency parts , again in @xmath20 space . the dotted line is the negative frequency components while the solid line represents the positive frequency components . we note that to very good accuracy , the interaction with the horizon has left the wave as a wave purely traveling to the left ( even though the fluid velocity is now sufficiently high to drag the left traveling wave to the right ) . the positive @xmath20 modes are negative frequency while the negative @xmath20 modes are positive frequency . 8 a narrower final pulse but otherwise the same as in figure [ fig3 ] . figures [ fig8 ] and [ fig9 ] give the final and initial wave packets for another higher frequency pulse while figures [ fig10 ] and [ fig11 ] give the spectra of these pulses again split into positive and negative frequency parts . we see in figure [ fig11 ] a much smaller conversion of positive to negative frequency , as would have been expected from the higher frequency of the initial pulse . all of these cases are for the @xmath48 dispersion function . let us now compare the above with the value for the total bogoliubov transformation expected on the thermal hypothesis . given the final mode represented by figure [ fig3 ] or [ fig8 ] , the expected number of particles resident in this mode in a thermal state of temperature @xmath77 is given by _ t=1 ( ^*(k)(k)-^*(k)(k ) ) dk where @xmath78 is the normalization factor = ( ^*(k)(k ) -^*(k)(k ) ) dk on the other hand , the expected number of particles in this mode , @xmath79 , if we assume that the input state was in the vacuum state is just given by the positive frequency portion of the initial wave form ( the bogoliubov coefficient ) . this is just the norm of the positive frequency part of the normalized wave .    comparing the two values , we find for the low frequency pulse in the soft ( @xmath48 ) dispersion function case that _ t= .02302 ;   = .02336 while for the high frequency pulse _ t=.0005364 ;   = .0005545 in each case the bogoliubov estimate of the particle production in that particular pulse is essentially the same as the thermal estimate . the change in high frequency behaviour of the theory makes little or no difference to the expected thermal output from the dumb hole . ( i suspect the error to be related to the fact that the final wave pulse is not entirely contained within the constant velocity region . ) moreover , in addition to the integrated bogoliubov transformation , one can test whether the full spectrum itself of the outgoing radiation is thermal . since the velocity in the simulation is time independent , the frequency of any wave in the continuum , infinite spatial section case is conserved . consider an final left moving pulse , contained entirely in the constant velocity region , with components @xmath80 . for any value of k , there are two possible values for the frequency @xmath35 as given in figure [ fig12 ] . since the pulses are left moving pulses , the relation between the frequency @xmath35 and @xmath20 is given by the left - moving branch labeled as such in the figure . since the final pulse is chosen to be a purely negative frequency pulse , it is , as mentioned previously , composed of purely positive values of k. now evolve the wave backward in time until one has evolved it past its interaction with the horizon and back out into the region where the velocity is a constant . because the problem is time independent , the frequency @xmath35 is conserved . thus the amplitude of the wave over some finite frequency regime must be the same . furthermore , we have already noted that the left moving wave remains a left moving wave , and thus remains on the same left moving branch of the dispersion relation . this amplitude is given by ^*(t , k ) ( t , k ) dk however , on the left moving branch , we have a definite relation between @xmath35 and @xmath20 , and between @xmath67 and @xmath81 . in particular for the left moving branch , @xmath82 . thus , we have that the magnitude of the wave over a frequency interval @xmath83 is f(k ) |(t , k)|^2 dkddwhere the integral is over all values of @xmath20 for which @xmath84 has the desired value . from the plot we see that along the left moving branch there are three values of @xmath20 for any value of @xmath35 near zero . if we choose our initial pulse to have @xmath20 positive , two of these lie on the negative frequency branch , and one lies on the positive frequency branch . now the initial pulse is chosen to be non- zero only for small values of @xmath20 . after the interaction with the horizon , the pulse has been blue shifted , and has non - zero values for @xmath20 lying near @xmath85 , the wave number where the phase velocity of the waves in the constant velocity regions ( @xmath86 ) is zero . let us define the three values of @xmath20 at which the value of @xmath35 is the same as @xmath87 , for the value of k nearest @xmath88 , @xmath89 for the nearest value of @xmath20 to @xmath90 and @xmath91 for the value nearest @xmath92 . under the thermal hypothesis , the initial pulses , obtained by transporting the final pulse backward in time , will be composed of positive and negative frequency components , with the two magnitudes in a frequency interval related by    @xmath93 is the group velocity at wave vector @xmath20 . furthermore , we should have    ( k_i)^2/v_g(k_i)| only if these relations hold will the bogoliubov transformation satisfy the thermal hypothesis for all final wave forms .    in figure [ fig13 ] i have plotted the three functions @xmath94 , the positive frequency part of the initial pulse , @xmath95 , the predicted positive frequency part of the inital pulse from the final pulse and @xmath96 , the prediction of the positive frequency part of the final pulse from the negative frequency part of the final pulse . these are plotted against the frequency @xmath35 . these are for the final pulse of figure [ fig3 ] and associated initial pulse of figure [ fig6 ] . note that if the first two agree , the last one s agreeing is a test of the conservation of the norm by the evolution . because of the finite size of the region and the periodic boundary conditions , only a finite number of @xmath20 values are calculated , but it is clear that these three curves are identical to numerical accuracy . note that the temperature is of order 1.3 , and this graph thus represents about 10 e - foldings of the thermal function . thus we can conclude that not only are the two particular final wave forms chosen for the soft dispersion relation thermal , but that any final wave form would be thermal . thus despite the change in the high frequency dispersion function , the output of the dumb hole will continue to be thermal . the spectrum , and not just the integrated intensity , is a thermal spectrum .    the dispersion function with @xmath48 present a relatively gentle transition from the low frequency linear regime to the high frequency constant regime . if we alter the dispersion function to give a much more sudden transition , namely the dispersion function with @xmath49 ( see fig . [ fig14 ] ) , does this alter the thermal spectrum of the emitted radiation ? the answer is no . the sharp dispersion function has a much longer initial pulse for the same situation as used above for the @xmath48 case . in order to allow the initial pulse to be contained entirely within the constant velocity region , i have scaled the problem by a factor of 2 . thus the transition @xmath45 in figure [ fig14 ] is larger by 2 to 1024 , the velocity profile of figure [ fig15 ] has been rescaled by a factor of 2 , increasing the temperature to approximately 2.7 . the final pulse of figure [ fig16 ] is just figure [ fig3 ] rescaled by a faxtor of 2 . in figure [ fig17 ] i present the initial pulse which must have produced the final pulse of figure [ fig16 ] in the case in which @xmath49 . despite the far more sudden transition between the regimes in the dispersion function , the relation between the final and the initial pulses is again accurately given by a thermal relation as plotted in figure [ fig18 ] . again the final positive and negative frequency components of the initial pulse can be accurately predicted from the final pulse using the thermal hypothesis . why are the results thermal , and what lessons can we learn from this simple problem for the more interesting one of black hole evaporation ? the crucial lessons seems to be that , despite the naive derivations , black hole evaporation depends neither on the ultra high frequency behaviour of the fields near the black hole , nor on the behaviour of the field theory in imaginary time . the process of the thermal particle creation appears to be primarily a low , not high , energy phenomenon . the black hole and dumb hole have natural time scales- the the mass in the former case , and the inverse gradient of velocity in the latter . the behaviour of the fields near the horizon of the holes is governed by these relatively ( compared with plank or atomic ) time scales . in the vicinity of the hole , the fields experience changes on these scales . however , as far as the high frequency parts of the field are concerned , these time scales are very long and the changes are slow . the high frequency phenomena simply adjust adiabatically to these slow changes . if the state is the vacuum state at high frequencies , it remains the vacuum . one would expect this to remain true no matter what the form of the theory was at those high frequencies and energies . it is only in the regime where the natural time scale of the hole is of the same order as the time scale of the field theory that these slow changes near the hole can make their influence felt . i.e. , it is only for those wave whose frequencies are of the order of 1/m for the black hole that anything happens .    on that time scale , the structure near the black hole is that the null geodesics diverge exponentially . the time between neigbouring outgoing geodesics is exponential in the affine parameter along the geodesics . it is this exponential divergence ( shared by both the black hole and the dumb hole ) that leads to the thermal spectrum emitted by the hole . it is this exponential divergence at the low frequency time scale , and not at the highest frequencies of the theory that matters for the thermal emission . one can encapsulate this freedom of the low frequency thermal emission from traces of the high freqeuncy behaviour of the theory by the phrase that dumb holes ( and thus by hypothesis also black holes ) have no quantum hair . just as the red shift at the horizon wipes out all traces of the rich complexity of classical material which could have formed the black hole , it also seems to wipe out all traces of the quantum structures at high frequencies wich are red shifted to form the thermal radiation from the black hole . i would like to thank ted jacobson who revived my interest in dumb holes , and pointed out in conversations and in ref.@xcite that one of the crucial effects of the atomic nature of matter is to alter the high frequency dispersion relation of the sound waves . conversations with him have been very helpful to me in understanding dumb holes . i would also like to thank matt choptuik who insisted through the years on thinking about the stability of the solution scheme before coding rather than when the numerical instabilities become obvious . without the lessons he taught me i would not have solved this problem . in addition a conversation with him on the numerical techniques needed for this problem gave me the courage to attack it . this research was carried out under a fellowship from the canadian institute for advanced research and with an nserc research grant . the equations of motion for the system are solved via a set of implicit , symmetric in space and time , finite difference equations . the spatial grid @xmath97 is taken as a uniform grid , with @xmath98 . the term @xmath99(x_i)$ ] is evaluated by taking the fft of @xmath100 , @xmath101 , multiplying by @xmath102 , and then taking the inverse fft . the temporal grid is staggered , with @xmath103 defined on the time slices @xmath104 with uniform differences @xmath105 . the @xmath106 on the other hand are evaluated on the staggered sites @xmath107 . the grid points run from 0 to @xmath108 where @xmath109 . the periodic boundary conditions are enforced by taking @xmath110 and @xmath111 , and similarly for @xmath106 . the differencing step equations are taken to be _ m+1,i & - & _ t2_x ( v(x_i+1_m+1,i+1 -v(x_i-1_m+1,i-1 ) + & = & _ m , i+_t2_x ( v(x_i+1_m , i+1 - v(x_i-1)_m , i-1 ) + _ t _ m , i + _ m+1,i & - & _ t2_xv(x_i ) ( _ m+1,i+1 - _ m+1,i-1 ) + & = & _ m , i+_t2_xv(x_i ) ( _ m , i+1 - _ m , i-1 ) + _ t [ f(_x^2 ) ] _ m+1,i    this differencing technique was chosen because , if we take @xmath30 as a constant , these equations can be shown to be stable , and non - damping . both of these characteristics are crucial for the success of the integration . the stability is important because of the severe frequency shifting that takes place on the horizon , a frequency shifting which could easily incite any instabilities present in the code . the non - damping nature is also crucial because of the necessity of extracting small effects ( the bogoliubov coefficients ) from the result of the evolution . in particular , since we are testing precisely the high frequency effects , it is important not to hide these behind artificial viscosity at high frequencies .    to see the stability , assume that the velocity is constant , and fourier analyze the difference equations in both space and time . we get the characteristic equation ( 2i ( _ t/2 ) - iv_t_x ( k_x))^2 = _ t^2 f(-k^2 ) or ( _ t/2 ) = _ t ( + v_x ( 2 k _ x ) ) this will have real solutions for @xmath35 for all @xmath20 as long as @xmath112 is chosen to be sufficiently small . in particular the rhs of the equation must be less than unity for all values of @xmath20 . since @xmath113 is bounded by @xmath114 , and @xmath115 , and since @xmath116 , the right hand side is bounded by @xmath117 so we require _ t < _ xk_0_x+v_max in which case the @xmath52 are all real the equations are stable , and non - damping . this would not have been true had i employed an explicit scheme for the time difference equations .    because of the use of an implicit scheme , the time stepping requires the inversion of matrix equation . to solve this , i used a double sweep technique . to illustrate the technique , let us solve the equation v^0_iu_i - v^p_iu_i+1 + v^n_iu_i-1 = s_i [ defu ] both of the equations for the time stepping are of this form , with different coefficient vectors @xmath118 . these coefficients are the same at all times as they depend only on @xmath30 and @xmath119 . the solution will be found by assuming the recursion relation u_i = a_i u_i-1+b_i first solve the recursion relation by substituting for @xmath120 and @xmath121 using this recursion relation into eqn [ defu ] . setting the coefficients of @xmath122 to zero , we get recursion relations for the @xmath123 and @xmath124 : a_i= -v^n_ia_i+1 v^p_i+v^0_i + b_i= -v^p_i b_i+1a_i+1 v^p_i+v^0_i with @xmath125 . now find two solutions to the homogeneous equations u^0_i+1= a_i+1 u^0_i + b_i + u^1_i+1=a_i+1 u^1_i with @xmath126 and @xmath127 . neither of these are periodic .    to find a periodic solution to the non - homogeneous equations , first recursively solve b_i= s_i - v^p_ib_i+1 a_i+1 v^p_i+v^0_i with @xmath128 . then recursively solve u_i= a_i u_i-1 + b_i with @xmath129 this will not satisfy the boundary conditions . however , one can now subtract a suitable linear combination of @xmath130 and @xmath131 to enforce the periodic boundary conditions .
Abstract:

Answer: the naive calculation of black hole evaporation makes the thermal emission depend on the arbitrary high frequency behaviour of the theory where the theory is certainly wrong . using the sonic analog to black holes dumb holes i show numerically that a change in the dispersion relation at high frequencies does not seem to alter the evaporation process , lending weight to the reality of the black hole evaporation process . i also suggest a reason for the insensitivity of the process to high frequency regime .    # 1#2    figure # 1    _ # 2 _