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: stochastic modeling plays an important role in many areas of science , such as physics , ecology or chemistry @xcite . stochasticity may appear due to the lack of complete knowledge about all the relevant variables , the precise dynamics of the system or the interactions with the environment . in some cases , one can obtain a compact description of a complicated system considering only a few relevant variables but at the expense of losing deterministic predictability . often , probabilities for some fundamental processes can be assigned on the basis of symmetries and other considerations , or on empirical analyis , and the dynamics of the process can be derived bottom - up . stochasticity appears together with delay terms in many situations of interest , such as gene regulation @xcite , physiological processes @xcite or postural control @xcite . the combined effects of stochasticity and delay are , however , not completely understood . from the mathematical point of view , stochastic processes including delay are difficult to analyze due to the non - markovian character . most of the previous approaches have focused on stochastic differential equations , that consider continuous variables @xcite , or random walks in discrete time @xcite , where delay can be taken into account increasing the number of variables . models with discrete variables but continuous time are the natural description of many systems such as chemical reactions , population dynamics or epidemic spreading . in some cases , discreteness can be a mayor source of fluctuations , not well captured by continuous models @xcite . the approach with discrete variables and continuous time was used in @xcite . most often , the delay time is taken to be a constant with zero fluctuations . this is not very realistic in the applications , since it is unusual to have a deterministic delay when the rest of the dynamics is stochastic . we will take this consideration into account by allowing the delay times to be random variables with arbitrary probability density functions .    in this work we study some simple , yet general , stochastic birth and death processes including delay . we will develop tree different approaches to the analytical study of this kind of non - markovian processes , in the general case of stochastically distributed delay : a direct approach in subsection ( [ sec : simple ] ) , an effective markovian reduction in subsections ( [ sec : elaborated ] ) and ( [ sec : fullfeedback ] ) , and a master equation approach , together with a time - reversal invariance assumption , in section ( [ dc ] ) . the first direct approach method is interesting for its simplicity , but its application is limited to systems with first order reactions and without feedback . the second one , effective markovian reduction , is rather flexible and general and its development is one of the main advances of this paper . the last master equation approach complements the previous , giving information about the full probability distribution . the main limitation of all the approaches is the need to assume that completion times for delayed reactions are independent random variables ( independent of each other and of other variables of the system ) , although the initiation rates may depend on the state of the system , allowing , for example , for feedback and crowding effects , so we do not consider this limitation to be very relevant for practical applications . although our methodology is rather general , we present it here using specific examples that have been grouped in two categories : delay in the degradation ( section [ dd ] ) and delay in the creation ( section [ dc ] ) . we end the paper with a brief discussion and comments in section [ comments ] . some more technical details are left for the two appendices . we will start by studying simple stochastic birth and death processes that include delay in the degradation step . a process of this type was proposed in @xcite as a model for protein level dynamics with a complex degradation pathway . we consider first the simplest possible process including delayed degradation : @xmath0 that is , a particle @xmath1 is created at a rate @xmath2 and disappears (  dies " or  degrades " ) a time @xmath3 after created . we allow the delay time @xmath3 to be randomly distributed i.e. the lifetimes @xmath3 of the created particles are random variables , that for simplicity we consider independent and identically distributed , with probability density @xmath4 . although not considered in this paper , the case of non - identically distributed delay times , in particular a probability density that depends on the time from birth , can also be treated . however , as commented above , the case of non - independent delay times does not seem to be tractable with the methods we present below . we note first that distributed delay is completely equivalent to degradation at a rate that depends on the  age " @xmath5 ( time form creation ) of the particle , i.e. , processes @xmath6 are equivalent if the rate @xmath7 and the probability density of the delay @xmath4 are related by : @xmath8 with @xmath9 being @xmath10 the cumulative distribution of the delay - time . this is so because @xmath11 is the probability of dying at the time interval @xmath12 , if the particle is still present at @xmath5 , and so it is nothing but the probability @xmath13 that the delay time @xmath3 belongs to that same interval conditioned to the particle still being alive at time @xmath5 , an event with probability @xmath14 . in the notation of @xcite , @xmath7 is nothing but the _ conditional failure rate_. we take @xmath15 as the time origin , so the number of alive particles at time @xmath16 is @xmath17 for @xmath18 . let @xmath19 the probability of @xmath20 particles being alive at time @xmath16 . in the remaining of this subsection we assume that there is no feedback , in the sense that the creation rate @xmath2 is independent on the number of particles @xmath20 , but , for the sake of generality , we do allow it to be a function of time @xmath21 . the non - feedback assumption allows us to obtain a full analytical solution . as shown in the appendix , independently of the form of the delay distribution , @xmath19 follows a poisson distribution @xmath22 with average @xmath23 . if the creation rate , @xmath21 , is independent of time , a steady state is reached , in which the average number of particles is @xmath24 , again independently of the form of the delay distribution . we will now compute the time correlation function . we shall see that its analytical expression does depend on the form of the delay distribution . we start from the relation : @xmath25 with @xmath26 particles created after ( before ) @xmath16 . @xmath27 can be computed exactly as before ( now taking @xmath16 as the time origin ) , so we have : @xmath28 the evolution of the number of particles already present at @xmath16 depends on the age @xmath5 of these particles . their survival probability until time @xmath29 can be written as : @xmath30 where we used @xmath31 , so we find : @xmath32 from this , one easily obtains the correlation function : @xmath33(t , t)=\langle n(t)\langle n(t+t)|n(t ) \rangle\rangle-\langle n(t)\rangle\langle n(t+t)\rangle$ ] @xmath34(t , t)=\int_0^tdt'c(t')\hat{f}(t+t - t')\label{corrsimple},\ ] ] if @xmath35 , independent of time , a steady - state can be reached with correlation function @xmath36(t)=\lim_{t\to\infty}k[n](t , t)$ ] . for a constant rate @xmath37 , which would be equivalent to an exponential delay distribution @xmath38 , it has the usual exponential decay @xmath36(t)=(c/\gamma ) e^{-\gamma t}$ ] . for a fixed delay time @xmath39 , corresponding to @xmath40 , the correlation function is a straight line @xmath36(t)=c(\tau_0-t)$ ] for @xmath41 and @xmath36(t)=0 $ ] for @xmath42 . for other distributions of delay time , the correlation function adopts different forms , but it is always monotonically decreasing . in figure ( [ correlationsimple ] ) we plot the correlation function for two different types of distribution of delay , for different values of the variance of the delay . we see that the distribution with fatter tail displays a slower asymptotic decay , and that the decay is slower as the variance of the delay increases . numerical simulations , performed with a conveniently modified version of the gillespie algorithm @xcite , are in perfect agreement with this exact result , providing a check of its correctness . we remark that the functional form of the decay of the correlation function depends on the delay distributed and can differ from the exponential decay found in systems without delay .    ) , as a function of time , plotted in logarithmic scale , for two different types of delay distribution , gamma and lognormal , for two values of the variance of the delay : @xmath43 ( left panel ) and @xmath44 ( right panel ) ; in both cases the average delay is @xmath45 and the creation rate is @xmath46 . we also plot a exponential decay with exponent one ( dot - dashed line ) , for comparison . note that delay distributions with larger variance and fatter tayls display slower asymptotic decay . ( online version in colour . ) [ correlationsimple ] ]      we now consider a process including both instantaneous and delayed degradation steps : @xmath47 this is , particles are created at a rate @xmath2 and each particle can be eliminated by two processes : i ) instantaneous degradation at a rate @xmath37 ; ii ) delayed degradation , initiated at a rate @xmath48 but completed only a time @xmath3 after initiation . again , we will allow the delay - degradations times to be random variables that , for simplicity , will be independent and identically distributed with probability density function @xmath4 .    for the process to be completely defined , one has to specify if a particle that initiates delayed - degradation at time @xmath16 and thus will disappear at @xmath49 ( this kind of particles will be called  infected " ) , can also disappear before the completion of this reaction , through instantaneous degradation . in the most general case , this can happen at a rate @xmath50 , not necessarily equal to @xmath37 . note that , in the case of first - order degradation ( @xmath50 not dependent on the number of particles @xmath20 ) , this instantaneous degradation is completely equivalent to a system with @xmath51 , after modifying the distribution of the delayed - degradation times in the following way : @xmath52 that is , when instantaneous degradation is added to infected particles , the probability that the lifetime is equal to @xmath3 has two contributions : ( i ) a particle initially has a lifetime @xmath3 ( probability density @xmath4 ) and survives up to this time ( an event with probability @xmath53 ) ; ( ii ) a particle has a lifetime larger than @xmath3 ( probability @xmath54 ) , but survives up to @xmath3 ( probability @xmath53 ) and then undergoes instantaneous degradation ( at rate @xmath50 ) . the consideration of these two contributions leads straightforwardly to eq.([ftau2 ] ) . we see that omitting first order instantaneous degradation of infected particles comprises no loss of generality , given that the treatment is valid for general distributions of delay . if @xmath48 and @xmath37 are independent of @xmath20 , the process is equivalent to the one - variable system discussed in the previous subsection [ sec : simple ] ) with a conveniently modified distribution of delay : @xmath55 this comes from the fact that a particle may disapear at time @xmath3 because it did not disapeared or was infected before and is degraded instantaneously ( probability density @xmath56 ) or because it got infected at some previous time ( @xmath57 ) with an appropriate lifetime ( @xmath58 , probability density @xmath59 ) . this includes as particular cases the ones studied in @xcite . the results of subsection ( [ sec : simple ] ) allows us to obtain the full solution also in the general case of distributed delay . if @xmath48 or @xmath37 depend on @xmath20 the processes are not anymore equivalent , two variables are necessary and a new approach is needed for the analysis . in the following we develop this method . we will also consider the case in which the creation rate @xmath2 depends on the number of particles . the full process corresponds to the following two - variable system : @xmath60 where we have split the proteins into two types : @xmath61 are infected particles that will die precisely at a time @xmath3 ( itself a stochastic variable ) after being infected and @xmath62 are non - infected (  active " ) particles ( so @xmath63 ) . we allow the rates to depend on @xmath64 , the number of @xmath65 , active , particles , but not on @xmath66 , the number of @xmath67 , infected , particles which are considered to be  inert " ; this condition will be realxed in the next subsetion . following @xcite , we have introduced the auxiliary particles @xmath68 whose number is given by the stochastic variable @xmath69 . the introduction of @xmath68 will allow us to obtain the properties of @xmath66 by using the relation : @xmath70 where the discrete process @xmath69 is a sequence of step ( heaviside ) functions and its derivative must be understood as a series of dirac - delta functions . here we have introduced the family of  survival " stochastic processes @xmath71 defined in the following way : first , for each @xmath57 we obtain a value of @xmath72 independently drawn from the distribution @xmath4 . next , we set @xmath73 , if @xmath74 , and @xmath75 , otherwise . this can be considered as the indicator function of a virtual is defined for all @xmath57 , regardless if a particle is actually infected a time @xmath57 . however it only contributes to ( [ nonmarkov ] ) if a particle is actually infected at time @xmath57 , since only then @xmath76 particle that is infected at @xmath57 and survives up to a time @xmath77 . it follows from the definition that : @xmath78    expressions ( [ nonmarkov]-[corrsurvival ] ) are the main advances of this section and provide us with the necessary tools to derive the main properties of the stochastic process ( [ cx0 ] ) . in the case considered in @xcite there is a fixed delay ( @xmath40 ) and no instantaneous degradation of infected particles ( @xmath51 ) , so one has simply @xmath79 . the inclusion of the survival process @xmath71 allows us to consider the general case of distributed delay and rates depending on the state of the system . note that the process followed by @xmath80 is markovian as the delay only appears in variable @xmath67 , so the properties of @xmath81 can be obtained using markovian methods , and the properties of the variable @xmath66 can be derived afterwards using ( [ nonmarkov]-[corrsurvival ] ) . in particular , the first moments follow : @xmath82 using standard markovian methods @xcite , one can prove that the process @xmath80 is described by the master equation : @xmath83 with @xmath84 the step operator , @xmath85 . in this section , we allow the creation rate @xmath2 to depend on the number of @xmath65-particles , constituting a feedback term on the number of `` active '' particles . from the master equation one easily derives the equations for the moments , the first of them read : @xmath86 in the case that @xmath87 is a linear function of @xmath64 and @xmath37 and @xmath48 do not depend on @xmath64 ( and none of them depend on @xmath66 or @xmath81 ) , the system of equations is closed and can be solved . for non - linear systems , we will make use of van kampen s expansion @xcite . this is a standard systematic expansion of the master equation , that consists on assuming a deterministic , and a stochastic part for the variables that scale differently with a large parameter , @xmath88 ( typically the volume or system size ) i.e. @xmath89 , @xmath90 . one can then write the master equation for the new variables @xmath91 , @xmath92 and expand in powers of @xmath93 . the method is generically valid , provided that the rates depend on the variables only trough @xmath94 ( plus higher orders in @xmath95 ; a common factor depending on @xmath88 multiplying all rates is also acceptable ) , which is fulfilled by most systems of interest , and that the macroscopic equations have a steady state as a single attractor . the equations for the macroscopic components are : @xmath96\phi_a,\label{fia}\\ \frac{d\phi_z}{dt}&=&d(\phi_a)\phi_a.\end{aligned}\ ] ] the stochastic contributions , to first order in @xmath93 , read : @xmath97\langle\xi_a\rangle\label{xia},\\ \frac{d\langle\xi_z\rangle}{dt}&=&\widetilde{d}\langle\xi_a\rangle\label{xiz},\\ \frac{d\langle\xi_a^2\rangle}{dt}&=&-2\left[\widetilde{\gamma}+\widetilde{d}-c'(\phi_a)\right]\langle\xi_a^2\rangle+\left(\widetilde{\gamma}+\widetilde{d}\right)\phi_a+c(\phi_a),\\ \frac{d\langle\xi_z^2\rangle}{dt}&=&2\widetilde{d}\langle\xi_a\xi_z\rangle+\widetilde{d}\phi_a,\\ \frac{d\langle\xi_a\xi_z\rangle}{dt}&=&-\left[\widetilde{\gamma}+\widetilde{d}-c'(\phi_a)\right]\langle\xi_a\xi_z\rangle+\widetilde{d}(\langle\xi_a^2\rangle-\phi_a),\label{xiaxiz}\end{aligned}\ ] ] with @xmath98 , @xmath99 . usually , for the ansatz about the scaling of the variables to work ( and so the expansion ) , the equations for the macroscopic components must have a single stable fixed point . in this case , however , the equation for @xmath100 does not have a fixed point , and @xmath101 and @xmath102 grow without bound . this grow , nevertheless , is consistent with @xmath103(t)}=o(\omega^0)$ ] , and the expansion can still be applied . ( [ xia]-[xiaxiz ] ) is a system of closed linear equations and so can always be solved . to compute the time correlations of @xmath66 from eq.([correlationnz ] ) we need the time correlations of @xmath81 . we note that : @xmath104 and that @xmath105 ( for @xmath106 ) can be obtained integrating ( [ na]-[nz ] ) or ( [ xia]-[xiz ] ) . in the general , non - linear , case , using first order van kampen s expansion , one obtains , over the steady state : @xmath107,\ ] ] with @xmath108 and @xmath109 the solution of @xmath110 . the derivative that appears in ( [ correlationnz ] ) is : @xmath111,\ ] ] with @xmath112 . putting all the pieces together , one finally obtains : @xmath113(t)&=&\langle n_i(t_0)n_i(t_0+t)\rangle_\textrm{st}-\langle n_i\rangle_\textrm{st}^2\\ & = & \omega\widetilde{d}\phi_{a , st}\int_0^\infty dt'\hat{f}(t+t')+\omega\widetilde{d}u\langle\xi_a\xi_z\rangle_\textrm{st}\int_0^\infty ds\int_0^\infty dr\hat{f}(s)\hat{f}(r)e^{-u|t+s - r|}.\label{proteincorr}\nonumber\end{aligned}\ ] ] proceeding in a similar way , one can derive : @xmath114(t)&=&\omega\langle\xi_a\xi_z\rangle_\textrm{st}u\int_{0}^{\infty}dt'e^{-u(t+t')}\hat{f}(t')\\ k_\textrm{st}[n_i , n_a](t)&=&\omega\langle\xi_a\xi_z\rangle_\textrm{st}u\int_{0}^{\infty}dt'e^{-ut'}\hat{f}(t+t')+\omega\widetilde{d}\langle\xi_a^2\rangle_\textrm{st}\int_{0}^{t}dt'e^{-ut'}\hat{f}(t - t'),\\ k_\textrm{st}[n_a](t)&=&\omega\langle\xi_a^2\rangle_\textrm{st}e^{-ut},\end{aligned}\ ] ] with @xmath115(t)\equiv\langle n_u(t_0+t)n_v(t_0)\rangle_\textrm{st}-\langle n_u\rangle_\textrm{st}\langle n_v\rangle_\textrm{st}$ ] . this finally allows to express the correlation function for the total number of particles , @xmath116 , as : @xmath117(t)=k_\textrm{st}[n_i](t)+k_\textrm{st}[n_a , n_i](t)+k_\textrm{st}[n_i , n_a](t)+k_\textrm{st}[n_a](t).\label{proteincorrtot}\ ] ] in this case , the average of @xmath20 again depends only on the average delay , @xmath118 , but the second moment depends on the delay distribution in a more complicated way , through factors involving the integral of @xmath119 .    in figure ( [ figure : correlationprot ] ) this result is compared with numerical simulations , showing a very good agreement . note that the treatment of the delayed reactions is exact , the only approximation coming from the use of van kampen s expansion , which is needed when non - linearities are present , but whose error scales as @xmath93 . like in the previous case , the process in which the distribution of delay has fatter tail shows slower decay for the correlation function . the insets show the time correlation for the number of `` infected '' particles , @xmath67 , which gives the largest contribution to the difference between different distributions . symbols come from numerical simulations and lines from the theoretical analysis eqs.([proteincorr]-[proteincorrtot ] ) . the creation rate is @xmath120 , parameters values are : @xmath121 and @xmath122 . ( online version in colour . ) [ figure : correlationprot ] ]      we now consider the case in which the creation rate depends on all present particles @xmath123 with @xmath20 the total ( inert+active ) number of @xmath1-particles . as noted before , this single - variable model can account for instantaneous plus delayed degradation , in the case that the degradation and  contagion " rates , @xmath37 and @xmath48 before , do not depend on the state of the system . for simplicity , we restrict our attention to this case . this process can be treated with the approach of the previous subsection introducing the additional variable @xmath68 , @xmath124 with @xmath69 the corresponding random variable giving the number of @xmath68 particles . we see that : @xmath125 with @xmath71 the same as in the previous section . the probability distribution for @xmath126 follows a master equation of the form : @xmath127 details of the derivation of the master equation in systems with delay are given in the appendix . here , @xmath128 , with @xmath129 the probability density of the delay distribution , although , since we are only interested in the properties of variable @xmath20 , we will not be using this expression . the key step in this case is to note that eq.([mez ] ) allows us to derive the statistical properties ( moments and correlations ) of @xmath69 as a function of those of @xmath130 . then , using ( [ nfeedback ] ) we will be able to self - consistently derive the properties of @xmath20 . more specifically , the approach proceeds as follows : summing eq.([mez ] ) over @xmath20 , we can obtain an equation for the evolution of @xmath131 , but that still depends on @xmath20 ( in this step the contribution of the second term in eq.([mez ] ) vanishes ) : @xmath132 the two times probability distribution @xmath133 follows a similar equation . conditioning carefully , summing over the variable @xmath20 and considering separately the case @xmath134 ( which turns out to be singular ) , we find : @xmath135\langle c(n_a)|n_{z_1}\rangle p(n_{z_1},t_1)\nonumber % \left.\frac{d^2p(n_{z_1},t_1;n_{z_2},t_2)}{dt_1dt_2}\right|_{t_2=t_1}&=&\delta(0)\left[(1-e_{z_2})\delta_{n_{z_1},n_{z_2}}e_{z_1}^{-1}+(1-e_{z_2}^{-1})\delta_{n_{z_1},n_{z_2}}\right]\langle c(n_a)|n_{z_1}\rangle p(n_{z_1},t_1)\nonumber\end{aligned}\ ] ] from ( [ nfeedback ] ) we easily obtain : @xmath136 while ( [ mez2 ] , [ mez2 t ] ) imply : @xmath137 and we finally obtain the following set of integral equations for the moments : @xmath138 in the case of linear feedback , @xmath139 , this system of equations is closed . for non - linear systems , one can use van kampen s expansion as explained above . in the steady state , one finds : @xmath140(t)&=&\int_{0}^{\infty}dx\left[\int_{0}^{t+x}dyk_\textrm{st}[n](t+x - y)\hat{f}(x)\hat{f}(y)+\int_{t+x}^{\infty}dyk_\textrm{st}[n](t+x - y)\hat{f}(x)\hat{f}(y)\right]\nonumber\\ & + & \omega c(\phi_\textrm{st})\int_{0}^{\infty}dx\hat{f}(t+x)\label{corfullfeedback}\end{aligned}\ ] ]    eq . ( [ averagefullfeedback ] ) shows that the steady state number of particles depends only on the average delay . ( [ corfullfeedback ] ) shows that the correlations depend on the delay distribution in a non - trivial way . the analysis of this equation is left for future work . we now turn our attention to the case in which the creation reaction , that is initiated stochastically , takes a finite time to be completed . for simplicity we assume that the degradation reaction is instantaneous . schematically , we have : @xmath141 in this case , if the creation rate does not depend on the number of particles , @xmath20 , then the delay in the creation is completely irrelevant , since the probability that a new particle appears at time @xmath16 is equal to the probability that its creation started at a time @xmath142 , but this equal to the probability that a particle starts its creation at time @xmath16 ( with a shift in the time if @xmath2 is time - dependent ) , so the process is completely equal to one with instantaneous creation .    following @xcite we will adopt here an approach different from that of the previous sections , that , besides the moments , will allow us to obtain an expression for the full probability distribution . for completeness we will explain here the method in some detail . for additional considerations , the reader is referred to @xcite . in the appenix it is shown that the master equation of the process ( [ eq : processcreation ] ) is : @xmath143+(e^{-1}-1)\left[\sum_{n'=0}^{\infty}\int_0^\infty d\tau c(n')p(n',t-\tau;n , t)f(\tau)\right]\label{master_delaydistrib}.\end{aligned}\ ] ] the master equation ( [ master_delaydistrib ] ) can be written as : @xmath144+(e^{-1}-1)[\tilde c(n , t)p(n , t)],\ ] ] where the effective creation rate , @xmath145 , is given by : @xmath146 the conditional probability @xmath147 follows a master equation identical to ( [ master_delaydistrib ] ) with all the probabilities conditioned to @xmath148 at time @xmath149 . from it , and using that @xmath150 , we obtain the following evolution equation for the conditional average : @xmath151 for @xmath152 , with initial condition @xmath153 . the knowledge of the steady value @xmath154 , allows the calculation of the steady - state probabilities @xmath155 , obtained by imposing @xmath156 in eq.([master_effdelay ] ) , as @xcite : @xmath157 @xmath158 is fixed by the normalization condition . all is left to do now is to compute the effective creation rate @xmath159 . the effective creation rate will be computed using expression ( [ averagedelaygen ] ) . in the general case of nonlinear creation rate , we will use van kampen s expansion to linearize @xmath160 around the macroscopic component of @xmath20 . we have : @xmath161 , so @xmath162 using ( [ averagedelaygen ] ) we obtain : @xmath163 equation ( [ vkfluc2 ] ) is in general a non - linear integro - differential equation , that can be difficult to solve . here , however , we will focus on the cases in which ( [ vkfluc2 ] ) has a stable steady state as a single attractor , which is the solution of @xmath164 . this is the regimen in which the validity of van kampen s expansion is guaranteed . we reach now a delicate point . eq.([vkflucgen ] ) is a ( linear ) integro - differential equation . to solve it , we would need an initial condition in the whole interval @xmath165 but we only know a one - time condition @xmath166 . we will circumvent this difficulty by assuming that , over the steady state , the system is statistically invariant under time - inversion , which implies @xmath167 . this condition , together with the value of @xmath168 at time @xmath16 , allows to find the solution of ( [ vkflucgen ] ) . the time - reversal invariance assumption in the steady state is fulfilled by any markovian system that follows detailed balance . our system follows detailed balance ( as any one - step process @xcite ) , but , due to the presence of delay , it is not markovian . so the time - reversal invariance is an assumption , whose validity needs to be checked . it was shown in @xcite that in this system the assumption is approximately valid .    in the case of constant delay , @xmath169 , the time - reversal symmetric solution of ( [ vkflucgen ] ) is @xcite : @xmath170 and using ( [ cexpanded ] ) we finally obtain : @xmath171+c'(\phi_\textrm{st})h(\tau)n,\ ] ] where @xmath172 is the steady state solution of ( [ vkfluc2 ] ) . from eq.([probgeneral ] ) one can obtain the steady - state probabilities @xmath155 . the mean value and variance are given by : @xmath173 from ( [ vksigma ] ) one can see that , interestingly , in the case of negative feedback ( @xmath174 ) , as the delay is increased the fluctuations change from sub - poissonian ( @xmath175 ) to super - poissonian ( @xmath176 ) . this is illustrated in figure ( [ figure : phasespace ] ) , where we also show the line of the hopf bifurcation for the deterministic system . it is usually obtained that a negative feedback reduces the magnitude of the fluctuations @xcite , when delay present we see that this negative feedback can change totally its effect , giving rise to an increase of the fluctuations .    the time correlation function can also be obtained from ( [ xiresult ] ) , as : @xmath117(t)=\sigma^2_\textrm{st}h(t)\ ] ] in the case of negative feedback it becomes non - monotonic , developing peaks of alternating sign at approximately multiples of the delay , signaling the presence of stochastic oscillations . for positive feedback , the time correlation is always positive , but not necessarily monotonic . we will finish by noting that the `` effective markovian reduction '' method used in the previous section can also be used for the case of delay in the creation with feedback . to be completely general , we allow two delays , one in the creation ( with probability density @xmath177 ) , and one in the degradation ( with probability density @xmath178 ) . the process is schematized as follows : @xmath179 with @xmath180 random variables distributed according to @xmath181 . with the addition of two new variables , the process can be rewritten as : @xmath182 which allows us to note that : @xmath183 in this case , the survival function @xmath184 is defined as : @xmath185 , if @xmath186 , and @xmath187 , otherwise , being @xmath188 and @xmath189 random times obtained from the corresponding pdf s @xmath190 and @xmath191 . @xmath184 is equal to one if a virtual particle that initiated its creation at time @xmath57 finished it at some intermediate time @xmath192 and since then had a lifetime greater that @xmath193 , so that it is still alive at @xmath16 , being zero otherwise . it follows that:@xmath194 in the case that the creation rate @xmath160 does not depend on the number of @xmath1-particles , the number of @xmath68-particles follows a markovian process ( poisson process ) , and the properties of @xmath20 can be derived from ( [ nt ] ) . if the creation rate depends on the number of @xmath1-particles i.e. if feedback is present , the properties of @xmath81 can be derived formally as a function of @xmath20 and then the properties of @xmath20 can be derived self - consistently trough ( [ nt ] ) , as done in subsection ( [ sec : fullfeedback ] ) . plane , for creation with constant delay and a negative feedback given by a creation rate of the form @xmath195 ( note that @xmath196 is the strength of the negative feedback ) . the `` poissonian line '' , @xmath197=\langle n\rangle$ ] , obtained through the approximation ( [ vksigma ] ) , marks the transition from sub - poissonian to super - poissonian fluctuations , while the hopf line marks the hopft transition into oscillatory behavior in the deterministic system . parameters values are : @xmath198 @xmath122 . ( online version in colour . ) [ figure : phasespace ] ] in this paper we have analyzed general stochastic birth and death models that include delay . we have presented three different methods that together constitute a general toolbox to study stochastic models including delay .    in sub - section ( [ sec : simple ] ) we have shown that when the creation rate is independent of the state of the system ( no feedback ) and the initiation of the delayed degradation and the instantaneous degradation are first order reactions ( rate not depending on the state of the system ) , the process can be solved fully in an exact fashion for general distributions of delay , showing always poissonian character and a monotonically decreasing time correlation function given by ( [ corrsimple ] ) .    in sub - sections ( [ sec : elaborated ] ) , ( [ sec : fullfeedback ] ) we have considered a more general process with delay in the degradation step , allowing the initiation of the delay degradation and the instantaneous degradation to be higher order reactions , as well as the presence of feedback in the creation rate . the method allows to reduce the system to a markovian one , where usual techniques can be used . explicit expressions for the time correlation for general delay distributions were obtained . in this case the correlation might be non - monotonic , if feedback is present , but typically decreases monotonically . section ( [ dc ] ) shows that when the delay appears in the creation reaction and feedback is present , the delay typically has more dramatic consequences . in the case of fixed delay , it is shown that for negative feedback , the fluctuations are amplified as the delay increases , going beyond the level found when no feedback is present , and the time correlation function becomes oscillatory , alternating positive and negative values at approximately multiples of the delay . in the positive feedback case , again for fixed delay , the fluctuations are reduced with increased delay and the time correlation function remains always positive . we start by considering the case @xmath199 . for the sake of simplicity , we focus on the case with creation rate , @xmath2 , independent of time , but the generalization to time - dependent @xmath2 is straightforward . since the time origin is taken at @xmath15 , the probability of observing zero particles at time @xmath200 is equal to the following limit : @xmath201,\label{p0tini}\ ] ] with @xmath202 playing the role of a small time - increment and @xmath203 . this expression follows from the fact that , in order to find the system with zero particles at time @xmath16 , in every previous infinitesimal time interval ( @xmath204 ) one of the following two ( incompatible ) events must take place : either a particle is not created ( probability @xmath205 ) or a particle is created with a lifetime smaller that @xmath206 ( probability @xmath207 ) . we now have : @xmath208\delta t =- c\int_{0}^tdt'\hat{f}(t - t'),\ ] ] with @xmath209 , so we find @xmath210 following a similar line of reasoning , @xmath19 can be computed as : @xmath211\prod_{{0\le j\le m-1}\atop{j\ne i_1,i_2,\dots , i_n}}\left[1-c\delta t\hat{f}(t - t_{i_j})\right]\ ] ] this expression results from the consideration of choosing the times @xmath212 at which the @xmath20 particles are created and survive up to @xmath16 . the @xmath213-th particle is created with probability @xmath214 and survives up to @xmath16 with probability @xmath215 . the other factor comes from the fact that at the other time intervals either a particle is not created or it is created but dies before @xmath16 .    using @xmath216=e^{-c\int_0^tdt'f(t - t')}\ ] ] and replacing the sums by integrals in the limit @xmath217 @xmath218^n\ ] ] we finally obtain : @xmath219^n}{n!},\ ] ] that is , a poisson distribution with average @xmath220 . in the steady state ( found as the limit @xmath221 ) , the average becomes @xmath222 . remarkably , this poissonian character is completely independent of the form of the delay distribution . as commented above , this result can be easily generalized to the case in which the creation rate depends on time , @xmath223 , obtaining again a poisson distribution with average @xmath224 . here we derive the master equation of the process ( [ eq : processcreation ] ) . we consider first the case of fixed delay @xmath3 . we start with the following identity : @xmath225 it is immediate to see that @xmath226 . in the case of fixed delay , the second sum can be evaluated introducing a tree - times probability as : @xmath227 now , @xmath228 . expanding in a similar way the term @xmath229 , and taking the limit @xmath230 , we can obtain the master equation of the process : @xmath143+(e^{-1}-1)\left[\sum_{n'=0}^{\infty}c(n')p(n',t-\tau;n , t)\right]\label{master_delay}.\end{aligned}\ ] ] in the case of distributed delay , we start considering a discrete distribution of delays i.e. @xmath231 with corresponding probabilities @xmath232 . the continuum limit can then be obtained making @xmath233 . the creation term in ( [ identitycreation ] ) can be written as : @xmath234 now , @xmath235 , that is , the probability that a particle started its creation at time @xmath236 with a creation time equal to @xmath237 . replacing in the previous equation and performing the appropriate sums we obtain : @xmath238 that in the continuum limit reduces to @xmath239 . considering in a similar way the other terms in ( [ identitycreation ] ) and taking the limit @xmath240 one can obtain the master equation for distributed delay ( [ master_delaydistrib ] ) . we acknowledge financial support by the micinn ( spain ) and feder ( eu ) through project fis2007 - 60327 . l.f.l . is supported by the jae predoc program of csic . m. barrio , k. burrage , a. leier and t. tian . oscillatory regulation of hes1 : discrete stochastic delay modelling and simulation . _ plos comput biol _ , * 2*(9 ) : e117 1017 - 1030 , 2006 . doi:10.1371/journal.pcbi.0020117        j. g. milton , j. l. cabrera , t. ohira , s. tajima , y. tonosaki , c. w. eurich , s. a. campbell the time - delayed inverted pendulum : implications for human balance control _ chaos _ , * 19 * 026110 ( 1 - 12 ) , 2009 , doi : 10.1063/1.3141429    j. boulet , r. balasubramaniam , a. daffertshofer and a. longtin stochastic two delay - differential model of delayed visual feedback effects on postural dynamics _ a _ , * 368 * , 1911 423 - 438 , 2010 . doi : 10.1098/rsta.2009.0214        t. d. frank multivariate markov processes for stochastic systems with delays : application to the stochastic gompertz model with delay _ phys . rev . e _ , * 66 * 011914 ( 1 - 8 ) , 2002 . doi:10.1103/physreve.66.011914    t. d. frank , p. j. beek and r. friedrich fokker - planck perspective on stochastic delay systems : exact solutions and data analysis of biological systems _ phys . e _ , * 68 * 021912 ( 1 - 10 ) , 2003 . doi : 10.1103/physreve.68.021912            t. galla intrinsic fluctuations in stochastic delay systems : theoretical description and application to a simple model of gene regulation _ phys . * 80 * 021909 ( 1 - 9 ) , 2009 . doi : 10.1103/physreve.80.021909
Abstract:

Answer: we study general stochastic birth and death processes including delay . we develop several approaches for the analytical treatment of these non - markovian systems , valid , not only for constant delays , but also for stochastic delays with arbitrary probability distributions . the interplay between stochasticity and delay and , in particular , the effects of delay in the fluctuations and time correlations are discussed .


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: after the experimental observation of bose - einstein condensate ( bec ) @xcite , there have been many experimental studies to investigate different quantum phenomena in a laboratory previously not accessible for study in a controlled environment , such as , quantum phase transition @xcite , quantum collapse of matter wave under attraction @xcite , four - wave mixing of matter waves @xcite , formation of vortex lattice @xcite , interference of matter waves @xcite , josephson tunneling @xcite , anderson localization @xcite etc . the generation and the dynamics of self - bound large quantum waves have drawn much attention lately . there have been some studies of self - bound matter waves or solitons in one ( 1d ) @xcite or two ( 2d ) @xcite space dimensions . a ( self - bound ) bright soliton travels at a constant velocity in 1d , due to a cancellation of nonlinear attraction and defocusing forces @xcite . the collision between two analytic 1d bright solitons is always elastic @xcite with the solitons emerging after collision without a change of velocity and shape . the 1d soliton has been observed in nonlinear optics @xcite and in bose - einstein condensates @xcite . however , a two- or three - dimensional ( 3d ) soliton can not be realized for two - body contact attraction alone due to collapse @xcite . the 1d bec solitons studied experimentally @xcite are really quasi solitons behaving like real solitons at large velocities . at low velocities a deformation of their shapes is expected upon collision . there have been a few proposals for creating a self - bound 3d matter - wave state which we term a quantum ball . some of these proposals involve an engineering of the atomic scattering length using the feshbach resonance technique to generate dynamically stabilized solitons in 2d and 3d @xcite . others consider extra interactions usually neglected in dilute bec of alkali atoms to create a stationary localized state . in the presence of an additional nonlocal dipolar interaction a 2d bec soliton can be generated @xcite . it has been suggested by maucher _ _ @xcite that in the case of rydberg atoms , off - resonant dressing to rydberg nd states can provide a nonlocal long - range attraction which can form a quantum ball . the collapse instability can be stopped in this case by a repulsive contact interaction . in this paper we demonstrate that a _ very small _ repulsive three - body interaction in the presence of an attractive two - body contact interaction can generate a stable quantum ball . the collapse is stopped in this case by the repulsive three - body interaction . although , some theoretical suggestions for generating a quantum ball , in the presence of a repulsive core in the two - body atomic interaction , seem viable @xcite , so far there has not been success in their experimental realization . nevertheless , there are questions about the dynamics of a quantum ball which are very intriguing . the study of the dynamics of a quantum elementary particle like an electron to determine simultaneously its position and velocity is doomed to failure due to the heisenberg uncertainty relation . on the other hand , for a quantum ball the uncertainty relation is not of concern due to its large mass and it can be traced like a classical object by its position and velocity at each instant . as the quantum ball is self bound it can move like a classical ball obeying newton s first law of motion . however , very little is known about the interaction dynamics of two quantum balls and that of a quantum ball with other objects and we address these questions in this paper . we consider the mean - field gross - pitaevskii ( gp ) equation with the inclusion of a three - body interaction for the study of statics and dynamics of a quantum ball . the quantum balls are bound by an attractive two - body contact interaction in the presence of a repulsive three - body contact interaction . we use the realistic parameters of @xmath0li atoms , in this study , with a negative scattering length corresponding to two - body attraction @xcite . the effect of atom loss due to three - body recombination is included in the study of dynamics . the three body loss rate is not accurately known for this system @xcite for the parameter domain used in this study . we have chosen three - body loss rates that ensure our system does not decay significantly during our dynamical simulations . it is expected that the dynamics of quantum balls will be independent of the details of the mechanism responsible for self binding and we do not believe that the results obtained here are so peculiar as to have no general validity . in fact , a preliminary study revealed similar dynamics for quantum balls made of dipolar atoms @xcite . a stationary quantum ball can be formed for the two - body attraction above a critical value in the 3d gp equation for any finite three - body repulsion . the statical properties of the quantum ball are studied using a variational analysis and a numerical solution of the 3d gp equation . the variational and numerical results are found to be in good agreement with each other . the stability of the quantum ball is established numerically under a small perturbation introduced by changing the three - body interaction by a small amount , while the quantum ball is found to execute sustained breathing oscillation . a quantum ball can move freely without deformation along any direction with a constant velocity . we study the frontal and angular collisions between two quantum balls . only the collision between two integrable 1d solitons is truly elastic @xcite . as the dimensionality of the soliton is increased such collision is expected to become inelastic with loss of energy in 2d and 3d . in the present numerical simulation of frontal collision between two quantum balls , at sufficiently large velocities the collision is found to be quasi elastic when the two quantum balls emerge after collision with practically no deformation and without any change of velocities . unlike classical balls , obeying the newton s laws , there is no change in the directions of motion of the quantum balls after collision . however , upon impact with a rigid impenetrable plane the quantum ball bounces like a classical elastic ball obeying the usual laws of reflection . at small velocities the collision between two quantum balls is inelastic and the quantum balls form a single bound entity in an excited state executing breathing oscillation , which we call a quantum - ball breather . we present the 3d gp equation used in this study in sec . [ ii ] and a variational analysis of the same for an analytic understanding of the formation of the quantum ball . in sec . [ iii ] we present the numerical results for stationary profiles of a quantum ball . we present numerical tests of stability of a quantum ball under a small perturbation . the quasi - elastic nature of collision of two quantum balls at large velocities and the formation of a quantum - ball breather at small velocities are demonstrated by real - time simulation . we end with a summary of our findings in sec . we consider a quantum ball in the presence of a three - body interaction and the mean - field model appropriate for this study . the mean - field gp equation for @xmath1 atoms of mass @xmath2 is @xcite @xmath3 }   \phi({\bf r},t ) , \label{eq1}\end{aligned}\ ] ] where @xmath4 is the scattering length , and @xmath5 is the three - body interaction term . the negative scattering length @xmath4 represents two - body attraction and the positive @xmath5 to three - body repulsion .    for an analytic understanding of the formation of a quantum ball convenient variational approximation of eq . ( [ eq1 ] ) can be obtained with the following gaussian ansatz for the time - independent stationary wave function @xcite @xmath6,\end{aligned}\ ] ] where @xmath7 , @xmath8 is the width . the energy density corresponding to eq . ( [ eq1 ] ) is given by @xmath9 consequently , the total energy per atom @xmath10 becomes @xmath11 the width @xmath8 of a stationary quantum ball with negative energy corresponds to a global minimum of energy @xmath12 : @xmath13 without the quintic term ( @xmath14 ) the quantum ball of width @xmath15 is tantamount to an unstable towne s soliton @xcite . for stability a non - zero quintic term ( @xmath16 ) is necessary .     versus @xmath8 in @xmath17 m for different @xmath1 from eq . ( [ eq6 ] ) for @xmath18 and @xmath19 m@xmath20/s . variational ( b ) @xmath21 and ( c ) @xmath22 phase plots for @xmath18 and @xmath19 m@xmath20/s , respectively , illustrating the regions of formation of a stable and metastable quantum ball obtained from eq . ( [ eq6 ] ) . , title="fig : " ]   versus @xmath8 in @xmath17 m for different @xmath1 from eq . ( [ eq6 ] ) for @xmath18 and @xmath19 m@xmath20/s . variational ( b ) @xmath21 and ( c ) @xmath22 phase plots for @xmath18 and @xmath19 m@xmath20/s , respectively , illustrating the regions of formation of a stable and metastable quantum ball obtained from eq . ( [ eq6 ] ) . , title="fig : " ]   versus @xmath8 in @xmath17 m for different @xmath1 from eq . ( [ eq6 ] ) for @xmath18 and @xmath19 m@xmath20/s . variational ( b ) @xmath21 and ( c ) @xmath22 phase plots for @xmath18 and @xmath19 m@xmath20/s , respectively , illustrating the regions of formation of a stable and metastable quantum ball obtained from eq . ( [ eq6 ] ) . , title="fig : " ] in the numerical calculation , we use the parameters of @xmath0li atoms , e.g. , @xmath23 @xcite and @xmath24 amu , where @xmath25 is the bohr radius . unlike the 1d case , the 3d gp equation ( [ eq1 ] ) does not have analytic solution and different numerical methods , such as split - step crank - nicolson @xcite and fourier spectral @xcite methods , are used for its solution . we solve the 3d gp equation ( [ eq1 ] ) numerically by the split - step crank - nicolson method using both real- and imaginary - time propagation in cartesian coordinates using a space step of @xmath26 @xmath17 m and a time step of @xmath27 ms in all calculations @xcite . all imaginary - time simulations were performed in a box of size @xmath28 unless otherwise stated . imaginary - time simulation is employed to get the lowest - energy bound state of a quantum ball , while the real - time simulation is to be used to study the dynamics using the initial profile obtained in the imaginary - time propagation @xcite . there are different c and fortran programs for solving the gp equation @xcite and one should use the appropriate one . in the imaginary - time propagation the initial state was taken as in eq . ( [ eq1 ] ) and the width @xmath8 set equal to the variational solution obtained by solving eq . ( [ eq6 ] ) . the convergence will be quick if the guess for the width @xmath8 is close to the final width .    , ( b ) energy @xmath29 and ( c ) peak atom density @xmath30 versus the number of @xmath0li atoms @xmath1 in a quantum ball for three different three - body interactions @xmath31 m@xmath20/s , @xmath32 m@xmath20/s , @xmath33 m@xmath20/s . , title="fig : " ] , ( b ) energy @xmath29 and ( c ) peak atom density @xmath30 versus the number of @xmath0li atoms @xmath1 in a quantum ball for three different three - body interactions @xmath31 m@xmath20/s , @xmath32 m@xmath20/s , @xmath33 m@xmath20/s . , title="fig : " ] , ( b ) energy @xmath29 and ( c ) peak atom density @xmath30 versus the number of @xmath0li atoms @xmath1 in a quantum ball for three different three - body interactions @xmath31 m@xmath20/s , @xmath32 m@xmath20/s , @xmath33 m@xmath20/s . , title="fig : " ]    the variational width of a stationary quantum ball can be obtained from a solution of eq . ( [ eq6 ] ) . if a minimum of energy ( [ eq5 ] ) exists , it can be either a global minimum with negative energy corresponding to a stable state or a local minimum of positive energy corresponding to a metastable state . the energy @xmath34 at @xmath35 even for a very small non - zero @xmath5 , and @xmath36 as @xmath37 . hence collapse is not allowed in the presence of a very small three - body repulsion . for certain values of the parameters there is a negative energy region between these two limiting values where the global minimum corresponding to a stable stationary state is located . for other sets of parameters , the energy changes monotonically between the above two limiting values without ever becoming negative or may have a local minimum with positive energy corresponding to a metastable state . this is illustrated in fig . [ fig1](a ) for @xmath0li atoms by a plot of @xmath12 versus @xmath8 for different @xmath1 from eq . ( [ eq5 ] ) for @xmath18 and @xmath19 m@xmath20/s . for @xmath38 there is no minimum of @xmath12 and there can not be a quantum ball . for @xmath39 there is a minimum at positive energy corresponding to a metastable state . finally , for @xmath40 and 10000 there are minima at negative energies corresponding to stable states . the parameter domains for the formation of stable and metastable states are shown in @xmath21 and @xmath41 phase plots for @xmath18 and @xmath19 m@xmath20/s , respectively , in figs . [ fig1](b ) and ( c ) . an interesting scaling relation @xmath42 is noted in fig . [ fig1](b ) . although there is a lower limit on the number of atoms @xmath1 for the formation of a stable quantum ball , viz . figs . [ fig1](b ) and ( c ) , there is no upper limit on @xmath1 . in the following we will only be concerned with the global minimum with negative energy corresponding to a stable stationary state . next we compare in fig . [ fig2](a ) the numerical and variational root - mean - square ( rms ) radius @xmath43 of a @xmath0li quantum ball versus number of atoms @xmath1 for three different values of the three - body term : @xmath31 m@xmath20/s , @xmath32 m@xmath20/s , @xmath33 m@xmath20/s . the variational result for the rms radius is given by : @xmath44 , where @xmath8 is the equilibrium variational width . for small @xmath1 , the quantum balls are well localized with small size and the agreement between numerical and variational results is better . in fig . [ fig2](b ) we show the numerical and variational energies @xmath29 of a quantum ball versus @xmath1 for different @xmath5 . the energy of a bound quantum ball is negative in all cases and its absolute value is plotted . in fig . [ fig2](c ) the numerical and variational peak atom density @xmath45 of the quantum ball versus @xmath1 is illustrated for different @xmath5 .     of a @xmath0li quantum ball for different @xmath1 and @xmath5 . ]    to study the density distribution of a @xmath0li quantum ball we calculate the reduced 1d density defined by @xmath46 in fig . [ fig3 ] we plot this reduced 1d density as obtained from variational and numerical calculations for different @xmath1 and @xmath5 . for a fixed three - body term @xmath5 , the quantum ball is more compact with the decrease of number of atoms @xmath1 . for a fixed number of atoms @xmath1 , the quantum ball is more compact for a small three - body term @xmath5 . the agreement between the two densities is better for the compact quantum balls of smaller size as in fig . [ fig2 ] .     and @xmath47 m@xmath20/s moving with velocity @xmath48 cm / s calculated with different space steps @xmath49 . ] we performed numerical tests of stability of the @xmath0li quantum balls under a small perturbation ( details not reported here ) . we considered a quantum - ball wave function as calculated by imaginary - time propagation and performed real - time propagation with the imaginary - time wave function under a small perturbation introduced at @xmath50 upon changing the scattering length by less than @xmath51 . after this sudden perturbation the quantum ball starts a breathing oscillation . the continued oscillation of the quantum ball over a long interval of time establishes its stability .    to study the dynamics of the quantum balls , we need to set these in motion . this can be achieved by multiplying the imaginary - time wave function by a phase @xmath52 with @xmath53 , where @xmath54 and @xmath55 are momentum and velocity , and real - time simulation is then performed using these wave functions for the study of moving quantum balls with velocity @xmath55 . however , to achieve the desired velocity numerically an accurate wave function calculated over a large space domain and small space and time steps are needed ; otherwise the numerically generated velocity is always less than the expected velocity @xmath55 , except when @xmath55 is very small . this is illustrated in fig . [ fig5 ] where we plot the displacement versus time for the moving quantum ball , along the @xmath56 axis , of @xmath57 and @xmath47 m@xmath20/s after multiplying its wave function by @xmath58 , so that @xmath59 @xmath17m@xmath60 m@xmath61 leading to @xmath48 cm / s . the numerically generated dynamics was obtained with space steps @xmath62 @xmath17 m , 0.05 @xmath17 m , and 0.025 @xmath17 m and illustrated in fig . [ fig5 ] . the ideal theoretical result for @xmath63 is also shown . the result for @xmath64 @xmath17 m is satisfactory and this value of space step has been used in all calculations .    in the following we study the collision dynamics of quantum balls , where we use a three - body term @xmath5 with dissipation corresponding to a loss of atoms from the quantum ball due to molecule formation . to test the solitonic nature of the present quantum balls , we study the frontal head - on collision of two quantum balls at large velocity including a three - body term with absorption : @xmath65 m@xmath20/s . the imaginary - time profiles of the quantum balls shown in fig . [ fig3 ] with @xmath57 and @xmath47 m@xmath20/s each are used as the initial wave functions in the real - time simulation of collision , with two identical quantum balls placed at @xmath66 @xmath17 m at @xmath50 . to set the quantum balls in motion along the @xmath56 axis in opposite directions the respective imaginary - time wave functions are multiplied by @xmath67 and real - time simulation is then performed using these wave functions for the study of dynamics . the corresponding velocity is @xmath48 cm / s . to illustrate the dynamics , we plot the isodensity contour of the colliding quantum balls in fig . [ fig5 ] at different times obtained by real - time simulation over a box of size @xmath68 . the initial profiles of the balls are shown in fig . [ fig5](a ) . the balls come close to each other in fig . [ fig5](b ) , coalesce to form a single entity in figs . [ fig5](c ) and ( d ) , form two separate balls in fig . [ fig5](e ) , and are well separated in fig . [ fig5](f ) with identical profiles as in fig . [ fig5](a ) . during collision , in figs . [ fig5](c ) and ( d ) the identity of the two separate balls give rise to a larger object which eventually breaks up into two quantum balls . considering the three - dimensional nature of collision , the distortion in the profile of the quantum balls after collision is found to be negligible , recalling that their identities were lost during collision . it is useful to contrast the collision shown in fig . [ fig5 ] with the corresponding elastic frontal collision of two classical balls . in the quantum collision the identity of the balls is lost during interaction and a quantum ball can not be followed during collision like a classical ball . apart from that , the position and velocity of the quantum balls before and after collision are identical with those of the classical balls in this process . besides the frontal collision considered in fig . [ fig5 ] we also consider two types of nonfrontal collisions including a three - body term with absorption . first , we consider the collision between two quantum balls moving along the @xmath56 axis but on laterally displaced tracks . at @xmath50 two balls of @xmath57 and @xmath65 m@xmath20/s each are placed at @xmath66 @xmath17 m , @xmath69 @xmath17 m , respectively , and set into motion along @xmath56 axis in opposite directions by multiplying the respective imaginary - time wave functions by @xmath70 and performing real - time simulation in a box of size @xmath71 with these wave functions for the study of dynamics . the collision is illustrated in fig . [ fig6](a)-(f ) through successive snapshots of isodensity contours of the system before , during and after collision . in this case in figs . [ fig6](a ) and ( b ) the balls approach along the @xmath56 axis , in figs . [ fig6](c ) and ( d ) they interact by losing identity , and in figs . [ fig6](e ) and ( f ) they eventually come out of the interaction region undeformed while moving along the @xmath56 axis maintaining their original trajectories and conserving their velocities . this collision has no classical analogue . in the elastic collision of two classical balls in this case , the balls will be deflected from their original trajectories conserving energy and momentum . another type of collision of interest is the angular collision of two quantum balls which we now study . for this purpose , at @xmath50 two balls of @xmath57 and @xmath65 m@xmath20/s each are placed at @xmath66 @xmath17 m , @xmath72 @xmath17 m , respectively , and set into motion towards the origin @xmath73 with equal velocities by multiplying the respective imaginary - time wave functions by @xmath74 to set the quantum balls in motion with a velocity of 18.9 cm / s and performing real - time simulation in a box of size @xmath75 with these wave functions for the study of dynamics . again the isodensity profiles of the quantum balls before , during , and after collision are shown in figs . [ fig7](a)-(b ) , ( c)-(d ) , and ( e)-(f ) , respectively . in this case the balls again come out after collision undeformed maintaining their original trajectories and conserving their velocities . if we contrast this collision with the corresponding elastic collision of two classical balls , the position and velocity of the quantum balls before and after collision are identical with those of the classical balls in this process . however , again the quantum balls lose their identity during the collision . the elastic interaction of a quantum ball with external objects is also of interest . for this purpose we consider its interaction with a rigid elastic plane upon perpendicular and angular impacts including a three - body term with absorption . to study the vertical impact with a rigid elastic plane we place a quantum ball of @xmath57 and @xmath65 m@xmath20/s at @xmath76 and q@xmath77 at @xmath50 and set it in motion along the @xmath56 axis with a velocity of 18.18 cm / s by multiplying the imaginary - time wave function by @xmath78 and performing real - time simulation with this wave function in a box of size @xmath79 . the reflection of the ball from the wall is achieved by just imposing reflecting boundary condition at the surface in the crank - nicolson algorithm @xcite . the interaction dynamics in this case is illustrated by successive snapshots of isodensity contour before , during and after interaction with the rigid plane in figs . [ fig8](a)-(f ) . the quantum ball moves without any deformation in figs . [ fig8](a)-(b ) , gets deformed in proximation of the rigid plane in figs . [ fig8](c)-(d ) , bounces off without any deformation and without any change of speed and with the direction of motion reversed in figs . [ fig8](e)-(f ) . the dynamics of the quantum ball is the same as that of an elastic classical ball except near the rigid plane when the quantum ball gets deformed . the interaction of a quantum ball upon angular impact with a rigid plane is studied next . a quantum ball of @xmath57 and @xmath65 m@xmath20/s is placed at @xmath80 @xmath17 m and set into motion towards @xmath81 @xmath17 m , @xmath77 by multiplying the imaginary - time wave function by @xmath82 to set the quantum ball in motion with an initial speed of 18.9 cm / s and performing real - time simulation with this wave function in a box of dimension @xmath83 . the reflection of the ball from the wall is achieved by imposing reflecting boundary condition at the surface in the crank - nicolson algorithm @xcite . the uniform motion of the quantum ball without deformation before and after the collision are shown in figs . [ fig9](a)-(b ) and [ fig9](e)-(f ) , respectively , while its deformation in the proximity of the rigid plane is shown in figs . [ fig9](c)-(d ) . the quantum ball bounces like an elastic classical ball with the same speed obeying the classical law of reflection . the results presented so far demonstrate beyond doubt that for large velocities the quantum ball interacts elastically with another quantum ball or an external rigid plane with the conservation of kinetic energy . however , the collision is inelastic at small velocities and large deformation of the quantum ball is possible with the nonconservation of kinetic energy . we performed simulation at small velocities ( about one tenth of the velocities considered so far ) for the interaction with a rigid plane . we find that the quantum ball bounces off the surface at such low incident velocities but with some deformation in shape and nonconservation of kinetic energy resulting in a reduction in the final velocity .    to study the inelastic collision dynamics we consider two quantum balls each with @xmath57 and a three - body term with dissipation , @xmath84 m@xmath20/s , place them at @xmath85 @xmath17 m , @xmath86 and set them in motion in opposite directions along the @xmath56 axis with a small velocity : @xmath87 cm / s . in this case we have used a smaller dissipative term to avoid a large loss of atoms over long - time collision dynamics for a small velocity . the dynamics is started by multiplying the respective imaginary - time wave functions by @xmath88 ) and performing real - time simulation for the study of dynamics . the imaginary and real - time simulations were performed in boxes of size @xmath89 and @xmath90 , respectively . the dynamics is illustrated by a plot of the time evolution of 1d density @xmath91 in fig . [ fig10 ] ( a ) and the corresponding two - dimensional contour plot is shown in fig . [ fig10 ] ( b ) . the two quantum balls come close to each other at @xmath92 and coalesce to form a quantum - ball breather and never separate again . the combined bound system remain at rest at @xmath92 continuing small breathing oscillation because of a small amount of liberated kinetic energy which creates the quantum - ball breather . hence at sufficiently small incident velocities the collision of two quantum balls lead to the formation of a quantum - ball breather and at large velocities one has the quasi - elastic collision of two quantum balls . we demonstrated the creation of a stable , stationary bec quantum ball ( a self - bound bec ) under attractive two - body and repulsive three - body contact interactions employing a variational and full numerical solution of the 3d gp equation . the statical properties of the quantum ball are studied by the variational approximation and a numerical imaginary - time solution of the 3d gp equation . the dynamical properties are studied by a real - time solution of the gp equation including an absorptive three - body term . the quantum ball can move with a constant velocity without deformation . at large velocities , the collision between the two quantum balls is quasi elastic with no visible deformation of the final quantum balls . we studied head - on and angular collisions of two quantum balls . in all cases , unlike classical balls , the balls come out of the collision region without deformation maintaining their velocities and directions of motion unchanged . in elastic collision of two classical balls one can have a change in the direction of their motion subject to energy and momentum conservation . at small velocities , the collision between two quantum balls is inelastic with the formation of a quantum - ball breather after collision .     as collapse is not allowed , even in the presence of a very small three - body repulsion , the present suggestion of realizing a trapless bec quantum ball seems to be attractive from an experimental point of view . the size of a trapped dipolar bec is determined by the harmonic oscillator lengths of the trap , whereas the size of the present quantum ball is determined by the two - body and three - body interactions . one should start with a tapped dipolar bec for @xmath93 where no droplet can be formed , viz . [ fig1](c ) . now using the feshbach resonance technique , one should make the scattering length more attractive to enter the stable domain from unstable domain in fig . [ fig1](c ) . if the harmonic trap is weak then initial size of the trapped bec should be large , and by increasing the two - body attraction the size of the quantum ball could be made much smaller : the sudden change in size will identify the trapless quantum ball . we thank the fundao de amparo  pesquisa do estado de so paulo ( brazil ) ( project : 2012/00451 - 0 and the conselho nacional de desenvolvimento cientfico e tecnolgico ( brazil ) ( project : 303280/2014 - 0 ) for support . m. h. anderson _ et al . _ , science * 269 * , 198 ( 1995 ) ; k. b. davis , m. o. mewes , m. r. andrews , n. j. vandruten , d. s. durfee , d. m. kurn , and w. ketterle , , 3969 ( 1995 ) ; c. c. bradley , c. a. sackett , j. j. tollett , and r. g. hulet , phys . lett . 75 , 1687 ( 1995 ) .                b. d. josephson , phys . lett . * 1 * , 251 ( 1962 ) ; f. s. cataliotti , s. burger , c. fort , p. maddaloni , f. mi- nardi , a. trombettoni , a. smerzi , m. inguscio , science * 293 * , 843 ( 2001 ) ; s. k. adhikari , phys . rev . a * 72 * , 013619 ( 2005 ) ; eur . . j. d * 25 * , 161 ( 2003 ) . p. w. anderson , phys . rev . * 109 * , 1492 ( 1958 ) ; j. billy , v. josse , z. zuo , a. bernard , b. hambrecht , p. lugan1 , d. clment , l. sanchez - palencia , p. bouyer , and a. aspect , nature ( london ) * 453 * , 891 ( 2008 ) ; g. roati , c. derrico1 , l. fallani1 , m. fattori , c. fort , m. zaccanti , g. modugno , m. modugno , and m. inguscio , nature ( london ) * 453 * , 895 ( 2008 ) ; s. k. adhikari and l. salasnich , phys . rev . a * 80 * , 023606 ( 2009 ) . v. m. perez - garcia , h. michinel , and h. herrero , phys . a * 57 * , 3837 ( 1998 ) ; f. k. abdullaev , a. gammal , a. m. kamchatnov , and l. tomio , int . j. mod b * 19 * , 3415 ( 2005 ) ; k. e. strecker , g. b. partridge , a. g. truscott , and r. g. hulet , nature ( london ) * 417 * , 150 ( 2002 ) ; l. khaykovich , f. schreck , g. ferrari , t. bourdel , j. cubizolles , l. d. carr , y. castin , and c. salomon , science * 256 * , 1290 ( 2002 ) . s. k. adhikari , , 063613 ( 2004 ) , 036608 ( 2004 ) ; h. saito and m. ueda , phys . lett . * 90 * , 040403 ( 2003 ) ; f. k. abdullaev , j. g. caputo , r. a. kraenkel , and b. a. malomed , phys . rev . a * 67 * , 013605 ( 2003 ) . b. satari , v. slavni , a. beli , a. bala , p. muruganandam , and s. k. adhikari , comput . * 200 * , 411 ( 2016 ) ; v. lonar , a. bala , a. bogojevi , s. skrbi , p. muruganandam , and s. k. adhikari , comput . . commun . * 200 * , 406 ( 2016 ) ; r. kishor kumar , l. e. young - s . , d. vudragovi , a. bala , p. muruganandam , and s. k. adhikari , comput . . commun . * 195 * , 117 ( 2015 ) ; v. lonar , l. e. young - s . , s. krbi , p. muruganandam , s. k. adhikari , and a. bala , comput . phys . commun . * 209 * , 190 ( 2016 ) ; l. e. young - s . , d. vudragovic , p. muruganandam , s. k. adhikari , a. balaz , comput . * 204 * , 209 ( 2016 ) .
Abstract:

Answer: we study the statics and dynamics of a stable , mobile , three - dimensional matter - wave spherical quantum ball created in the presence of an attractive two - body and a _ very small _ repulsive three - body interaction . the quantum ball can propagate with a constant velocity in any direction in free space and its stability under a small perturbation is established numerically and variationally . in frontal head - on and angular collisions at large velocities two quantum balls behave like quantum solitons . such collision is found to be quasi elastic and the quantum balls emerge after collision without any change of direction of motion and velocity and with practically no deformation in shape . when reflected by a hard impenetrable plane , the quantum ball bounces off like a wave obeying the law of reflection without any change of shape or speed . however , in a collision at small velocities two quantum balls coalesce to form a larger ball which we call a quantum - ball breather . we point out the similarity and difference between the collision of two quantum and classical balls . the present study is based on an analytic variational approximation and a full numerical solution of the mean - field gross - pitaevskii equation using the parameters of @xmath0li atoms .