Instruction: Write an abstract for the following acadamic paper.

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

Paper: the bak  sneppen ( b  s ) model is an abstract representation of a biological system that evolves according to the darwinian principles of random mutation and natural selection . it was introduced in @xcite in the context of self - organized criticality in systems with spatial interactions .    despite its simplicity , the b  s model captures some of the features that are believed to be characteristic of evolving biological systems . in particular , it predicts evolutionary activity on all time scales with long periods of relative stasis interrupted by bursts of activities , referred to as _ avalanches_. as a consequence of the absence of a characteristic time scale , evolutionary dynamics in the b  s model display long - range dependence in both the temporal and the spatial domain . it is thus suitable as an abstract representation of systems in punctuated equilibrium , a concept that was introduced in @xcite to explain the patterns observed in fossil records . for a more thorough discussion of these ideas , and for applications of the notion of punctuated equilibrium in other scientific disciplines we refer the reader to @xcite . in addition to its usefulness in abstractly representing some key features of palaeontology and macro - evolution , the b  s model has also been employed to analyze the evolution of bacteria in a controlled , competitive environment . in a series of key experiments @xcite , lenski and collaborators cultivated twelve initially identical populations of an e.  coli strain over several years and conserved samples at regular time intervals . they then determined the relative fitness of the conserved samples by putting them into direct competition with a sample taken from the initial populations and measuring their relative growth rates . in @xcite it was shown that the b  s model with random mutations qualitatively reproduces some of the experimental results on relative bacterial fitness obtained in lenski s long - term experimental evolution project @xcite . using an extended multi - trait variant , @xcite extended the explanatory scope of the b  s model to include experimental findings about the interplay of adaptation , randomness and history in bacterial evolution .      the bak  sneppen model characterizes each species in a biological system by a numerical fitness value between zero and one , which represents its degree of adaptedness to its environment and changes as the species evolves . further , each species is assumed to directly interact with exactly two other species , where it is left unspecified if such an interaction represents competition for resources , predator - prey relations or something different entirely . the bak  sneppen model can therefore be visualized as points on a circle , where each point stands for an ecological niche ( or a species occupying that niche ) and neighboring species interact with each other . as time progresses , the fitness landscape evolves in accordance with the following rules , representing in an abstract way the principles of random mutation and natural selection : at each time step , the least adapted species , i.e. the one with the smallest fitness parameter , is removed from the system ( becomes extinct ) and its place is immediately taken by a new species whose fitness is initially modelled as an independent uniformly distributed random variable . in order to take into account the effect of this change on the local environment , the fitness parameters of the two species to either side of the least fit one are also reset to random values . this can be thought of as those two species themselves becoming extinct and superseded by new ones , or as them undergoing mutations in response to their neighbour becoming extinct .      despite its apparently easy definition , the b  s model has withstood most attempts at mathematical analysis in the past . partial results have been obtained , however , in the context of rank - driven processes @xcite and mean - field approximations @xcite . early on it was conjectured based on simulations that the steady - state fitness distribution at a fixed site converges , in the limit of large populations , to a uniform distribution on the interval @xmath0 $ ] , where @xmath2 is approximately equal to 0.667 , but believed to be slightly larger than @xmath3 . there is only a small number of mathematically rigorous result about the bak  sneppen model ; in @xcite it is shown that the steady - state fitness at a fixed site is bounded away from one in expectation , independent of the number of species in the system ; a characterization of the limiting marginal fitness distribution , conditional on a set of critical thresholds , is given in @xcite ( see also @xcite ) . in @xcite , the author proposes to compute the steady - state fitness distribution as the fixed point of the one - step transition equation and uses this method to describe the asymptotic fitness distribution for four species in terms of a compact rational function . in the same paper it is shown that one can not find a similarly simple formula in the b  s model with five species , and that the fitness distribution of a randomly selected species at steady - state in this case is not only not rational , but not even a hypergeometric function .      in this paper , we revisit the bak  sneppen model with five species . in [ thm - main ] , our main theorem , we establish a representation of the steady - state fitness distribution for five species in terms of the solution of an explicit differential equation with hypergeometric coefficients . this steady - state distribution encodes information about the fitness attributes of species in a system that has evolved for a long time . for instance , one can deduce from it how fit , on average , a randomly selected species from the population is expected to be ; this is done in [ cor - marg ] . furthermore , since the steady - state distribution contains information about the joint fitness values of all species in the population , its knowledge allows to draw biologically relevant conclusions about qualitative properties of the system , such as the emergence of one or several dominant species , or the fragmentation of the eco - system into areas of different prevailing fitness . in our simple model , the symmetry of the initial configuration is preserved and no such phenomenon occurs . it is an interesting question whether in systems with a more complicated interaction between species , symmetry can be spontaneously broken . we envisage that similar representations for the asymptotic fitness distribution in larger systems may help pave the way towards a resolution of the question of whether or not , in the limit of infinitely many species , the fitness is indeed asymptotically uniformly distributed on the interval @xmath0 $ ] with @xmath1 . we speculate that the techniques developed in this paper can be generalized to analyze larger systems with more than five species . it seems plausible to expect that the asymptotic joint fitness distribution in such systems can still be characterized as the solution to a certain linear differential equation , even though the coefficients might no longer be hypergeometric functions . we adopt the following formalization of the bak  sneppen model from @xcite . initially , all fitness parameters are independent uniformly distributed and after @xmath4 evolutionary steps the state of the system is represented by the vector @xmath5 ^ 5 $ ] , where the @xmath6th component refers to the fitness of the @xmath6th species . the evolutionary dynamics of the system can be expressed formally by the equation @xmath7 ^ 5,\quad a\in\mathscr{b}([0,1]^5),\ ] ] where the one - step transition kernel @xmath8 encodes the dynamics of the model and is given by @xmath9 here , and in the following , all vector indices are taken modulo five . the sequence @xmath10 is a uniformly ergodic markov chain with absolutely continuous marginal distributions with densities @xmath11 ^ 5\to\mathbb{r}^+$ ] . this means that for any borel set @xmath12 ^ 5)$ ] , @xmath13 and that the random vectors @xmath14 converge in distribution to a steady - state limit @xmath15 . moreover , the @xmath4-step densities @xmath16 satisfy the recursion @xmath17 ^ 3 } { \boldsymbol{1}}_{\left\{\xi_2<\min(\xi_1,\xi_3,{{\mathbf}{x}}_{]\nu[})\right\}}g_k\left({{\mathbf}{x}}_{]\nu[_{{{\mathbf}{\xi}}}}\right){\mathbbm{d}}^3{{\mathbf}{\xi}},\ ] ] where the vectors @xmath18\nu[}\in[0,1]^2 $ ] and @xmath18\nu[_{{{\mathbf}{\xi}}}}\in[0,1]^5 $ ] are obtained from @xmath19 by dropping the @xmath20th , and @xmath21th components , or replacing these components by the components of @xmath22 , respectively . uniform ergodicity of the markov chain @xmath23 implies that the densities @xmath16 converge uniformly to the density @xmath24 of the unique invariant distribution of @xmath23 , which we recognize as the steady - state fitness distribution . we also introduce the notation @xmath25 where @xmath26 denotes the gauss hypergeometric function ( * ? ? ? * section 15.1 ) and @xmath27 is the imaginary unit . the following is the main result of the paper . [ thm - main ] the limiting density @xmath28 is given by @xmath29 ^ 5}({{\mathbf}{x}})\sum_{\nu=1}^5{q\left(\min{\{x_\nu , x_{\nu+1}\}},\max{\{x_\nu , x_{\nu+1}\}}\right)},\ ] ] where @xmath30 . here , @xmath31 and the function @xmath32 is the unique solution of the differential equation @xmath33 with boundary conditions @xmath34 the coefficients @xmath35 , @xmath36 , are hypergeometric functions given by    [ eq - coeff - ode ] @xmath37,\\ c_1(y ) = & -y c_0(y),\\ c_2(y ) = & \frac{6 } { y^3 + 2}\left[y \left(3 y^6 - 38 y^3 - 4\right ) \mathcal{g}'(y)+\left(15 y^6 - 10 y^3 + 4\right ) \mathcal{g}(y)\right],\\ c_3(y ) = & -12 y \left[y \left(4 y^3 - 1\right ) \mathcal{g}'(y)+\left(5 y^3 + 1\right ) \mathcal{g}(y)\right],\\ c_4(y ) = & -3 y^2 \left[y \left(y^3 + 2\right ) \mathcal{g}'(y)+\left(9 y^3 - 2\right ) \mathcal{g}(y)\right],\\ c_5(y ) = & y^3 \left(y^3 + 2\right ) \left[y \mathcal{g}'(y)-\mathcal{g}(y)\right].\end{aligned}\ ] ]    finally , the function @xmath38 is given by @xmath39{\mathbbm{d}}\xi}.\ ] ]    unfortunately , it does not seem possible to evaluate the differential equation @xmath33 in terms of known special functions . as a direct consequence of [ thm - main ] , we obtain the fitness distribution of a single species at steady - state by computing the one - dimensional marginal of [ eq - limdensity ] . [ cor - marg ] in the bak  sneppen model with five species , the fitness distribution of a randomly - selected species at steady - state is absolutely continuous with density @xmath40\mathbf{1}_{[0,1]}(x).\ ] ]    without loss of generality we can compute the distribution of the fitness of a fixed species , say the first one , which we call @xmath41 . averaging over the fitness values of the remaining species and using [ eq - limdensity ] leads to the expression @xmath42 ^ 4}{g(x,\xi_2,\ldots,\xi_5){\mathbbm{d}}^4\xi}=\left(\frac35 + 2\int_0^x{q(y , x){\mathbbm{d}}y}+2\int_x^1{q(x , y){\mathbbm{d}}y}\right)\mathbf{1}_{[0,1]}(x)\ ] ] for the density of @xmath41 . plugging in the explicit formula for @xmath43 stated in [ thm - main ] , differentiating once with respect to @xmath44 and using [ eq - defbcirc0 ] to eliminate @xmath38 results in the equation @xmath45 . the observation that @xmath46 together with @xmath47 completes the proof .    the graphs of the limiting density @xmath48 as well as the marginals of the @xmath4-step densities @xmath16 are depicted in [ fig - margdens ] , illustrating the convergence asserted in [ thm - main ] . moreover , [ fig - comparemargdens ] compares the cumulative distribution functions of the steady - state fitness value at a fixed site in die bak  sneppen model with three ( trivial ) , four ( * ? ? ? * theorem 1 ) and five ( [ thm - main ] ) species , as well as their conjectured limit  assuming @xmath49  in a system with an infinite number of species . clearly , a lot of work remains to be done to fill the gap between five and infinitely many species . 0.45     0.45 it was shown in ( * ? ? ? * proposition 2 ) that the density of the joint distribution of the fitness parameters after @xmath4 steps , starting from a uniform distribution @xmath50)^{\otimes 5}$ ] , is a polynomial given by @xmath51,\end{aligned}\ ] ] where the functions @xmath52 $ ] can be written as @xmath53 the coefficients of @xmath54 for @xmath55 , as computed from [ eq - recgk ] , are tabulated in [ tab - coeff - k1,tab - coeff - k2,tab - coeff - k3,tab - coeff - k4,tab - coeff - k5 ] . in order to identify the function @xmath56 featuring in [ thm - main ] , it is thus sufficient to compute the uniform limit @xmath57 this will occupy most of the rest of this section . first we recall an explicit recursion for the coefficients @xmath58 that was derived in ( * ? ? ? * proposition 5 ) . it corresponds to  and is derived from  the recursion for the densities @xmath16 . it obviates the need to evaluate any integrals and thus allows for the functions @xmath16 to be determined much more quickly and efficiently . [ prop - rec - alpha5 ] the coefficients @xmath58 vanish for @xmath59 and have the following properties :    a.   for @xmath60 , they satisfy @xmath61 , @xmath62 , as well as @xmath63\alpha_{i-2,0,k}+\left[\frac{1}{3}+\frac{1}{(i-2)i}\right]\alpha_{i-3,0,k}\\ & \quad+\frac{i+2}{i}\sum_{p=0}^{3k+1}\frac{\alpha_{i-1,p , k}}{p+1}-\frac{i+4}{2i}\sum_{p=0}^{3k+1}\frac{\alpha_{i-2,p , k}}{p+1}-\frac{i+2}{i}\sum_{p=0}^{i-2}\frac{\alpha_{i-2-p , p , k}}{p+1}\\ & \quad+\frac{i+4}{2i}\sum_{p=0}^{i-3}\frac{\alpha_{i-3-p , p , k}}{p+1}+\sum_{p=0}^{i-2}\frac{\alpha_{i-2-p , p , k}}{i - p } -\frac{1}{2 } \sum_{p=0}^{i-3}\frac{\alpha_{i-3-p , p , k}}{i - p},\quad i{\geqslant}3 . \end{split}\ ] ] b.   for @xmath64 , they satisfy the recursion @xmath65\alpha_{i-2,j , k}+\left[\frac{1}{3}+\frac{1}{i(i-2)}\right]\alpha_{i-3,j , k } ,            & i{\geqslant}3 .                       \end{cases}\ ] ]    the analysis of this three - dimensional recursion is simplified considerably by the fact that for each @xmath66 the sequence @xmath67 becomes eventually constant . more precisely , we have the following result which may be compared to ( * ? ? ? * lemma 2 ) and can be proved along the same lines . [ lemma - coeff - constant ] for each @xmath66 there exists a rational number @xmath68 such that @xmath69 for all @xmath70 . in particular , @xmath71 .    in [ tab - coeff - k1,tab - coeff - k2,tab - coeff - k3,tab - coeff - k4,tab - coeff - k5 ] the coefficients @xmath58 with @xmath70 are printed in bold and the frontier coefficients @xmath72 are marked with boxes . the assertion of [ lemma - coeff - constant ] is thus easily appreciated by visual inspection . in combination with [ lemma - coeff - constant ] allows to recognise the limiting function @xmath43 as the generating function of the limiting coefficients @xmath68 , i.e.  @xmath73 . before proceeding further we define convenient notation for the generating functions of the arrays @xmath58 and @xmath68 along various dimensions and with various indices held fixed .    @xmath74    here , the formal variables @xmath44 , @xmath75 and @xmath76 correspond to indices @xmath6 , @xmath77 and @xmath4 , respectively , and the symbol @xmath78 indicates summation over the index that it replaces . the next result establishes how passing to the limit @xmath79 can be accomplished at the level of generating functions . [ lemma - ggenred ] for each @xmath80 $ ] , it holds that @xmath81 equals @xmath82 . similarly , for non - negative integers @xmath66 , it holds that @xmath83    it suffices to prove the first claim , which follows directly from [ lemma - coeff - constant ] ; it allows us to write @xmath84 after multiplication by @xmath85 the first term vanishes as @xmath76 approaches one , whereas the second one converges to @xmath81 . we now begin analyzing the recursion for @xmath58 in more detail . throughout , we employ the powerful technique of generating functions as described in @xcite . in particular , we make use of the fact that a linear recursion equation for a sequence @xmath86 can be transformed into a differential equation for the generating function @xmath87 by multiplying the original recursion by @xmath88 and summing over @xmath89 . indeed , if the original recursion equation has polynomial coefficients , this property is shared by the resulting differential equation . [ prop - bcircj ] for positive integers @xmath77 , the generating function @xmath90 is given by @xmath91 , where @xmath92.\ ] ]    the last case of [ eq - recn5jgeq1 ] and [ lemma - coeff - constant ] imply that , for @xmath64 , the sequence @xmath93 satisfies the recursion @xmath94\beta_{i-2,j}+\left[\frac{1}{3}+\frac{1}{i(i-2)}\right]\beta_{i-3,j},\quad i{\geqslant}3.\ ] ] after multiplying by @xmath95 and summing over @xmath6 this translates into the differential equation @xmath96 b_{\circ , j}'''(x)=0\ ] ] for the generating functions @xmath97 . the initial conditions are @xmath98 where the last equality follows from the second case of [ eq - recn5jgeq1 ] . the general solution of [ eq - ode - bcircj ] is @xmath99 this can be obtained with the aid of a computer algebra system or checked using the power series representation of hypergeometric functions . the coefficients @xmath100 , @xmath101 , are determined by the initial conditions and are given by @xmath102 , @xmath103 , and @xmath104 , where @xmath105 and @xmath106 are explicit hypergeometric constants . observing that @xmath107 equals @xmath108 , which can be proved via zeilberger s algorithm @xcite , the claim follows . implies the decomposition @xmath109 and it thus only remains to compute the two functions @xmath110 and @xmath38 . unfortunately , the corresponding cases of [ eq - recn5j0,eq - recn5jgeq1 ] can not be analyzed by simply passing to the limit @xmath79 because they involve coefficients @xmath58 with @xmath111 . we therefore need to compute the two - dimensional generating function of the array @xmath112 , a task which is directly modelled after the proof of [ prop - bcircj ] . [ prop - g - circ - circ - j ] for positive integers @xmath77 , the generating function @xmath113 is given by @xmath114 , where @xmath115 . \end{split}\ ] ] in particular , @xmath116 . the proof proceeds by multiplying the @xmath117 case of [ eq - recn5jgeq1 ] by @xmath118 and summing over both @xmath6 and @xmath4 . after some algebraic manipulations one arrives at a differential equation similar to [ eq - ode - bcircj ] which can be solved using the same techniques . our next task will be to transform the first case of the recursion relation into differential form . this corresponds to the case @xmath64 and @xmath119 . the proof proceeds as before by multiplying the recursion equation by @xmath120 , summing over @xmath77 and @xmath4 and interchanging the order of summation where necessary . [ prop - ode1 - 2var ] the generating functions @xmath121 and @xmath122 satisfy the integral equation @xmath123    setting @xmath124 in [ eq - ode1 - 2var ] , applying [ lemma - ggenred ] and differentiating once with respect to @xmath44 yields the following result which will be one of the main ingredients in the proof of our main theorem . [ coro - ode1 - 1var ] the generating functions @xmath110 and @xmath38 satisfy the integro - differential equation @xmath125.\end{aligned}\ ] ]    a second differential equation relating @xmath110 and @xmath38 can be derived from [ eq - recn5j0 ] . the generating functions @xmath110 and @xmath38 satisfy the differential equation @xmath126b_{\circ,0}'(x)\notag\\ & \quad = 3(1-x)\left(1-\int_0^x{g(\xi){\mathbbm{d}}\xi}\right)b_{1,\circ}(x ) + \frac13\left[2+(1-x)^3\right]g(x)b_{1,\circ}'(x).\end{aligned}\ ] ]    using the @xmath119 case of [ eq - recn5jgeq1 ] the two @xmath127 sums can be eliminated from [ eq - recn5j0 ] . this allows to pass to the limit @xmath79 and to replace the coefficients @xmath58 by their steady - state values @xmath68 from [ lemma - coeff - constant ] . multiplying the resulting equation by @xmath95 and summing over @xmath6 yields [ eq - ode2 - 1var ] . the details of the computation are lengthy and hence omitted . the two equations , together with the initial conditions that are tabulated in [ tab - coeff - k1,tab - coeff - k2,tab - coeff - k3,tab - coeff - k4,tab - coeff - k5 ] , are already enough to numerically compute the functions @xmath110 and @xmath38 . it is , however , possible to simplify the problem even further to a single differential equation of order five .    for better readability we introduce the antiderivatives @xmath128 a direct evaluation of the last antiderivative gives [ eq - defg ] . after the substitution @xmath129 , [ eq - ode1 - 1var , eq - ode2 - 1var ] can be rewritten in terms of these functions as @xmath130 differentiating the latter equation three times eliminates the zeroeth and first derivative of @xmath131 ; replacing the remaining second through fifth derivatives with expressions derived from the first equation of the last display produces the fifth - order differential equation @xmath132 . the coefficients @xmath133 are given by @xmath134,\\ c_2(y ) = & 6 \left(2 - 5 y^3\right ) \mathcal{g}(y)-6 y\left(13 y^3 + 2\right ) \mathcal{g}'(y)-9 y^5 \mathcal{g}''(y)+y^3 \left[\left(11 y^3 + 4\right ) \mathcal{g}^{(3)}(y)+\left(y^3 + 2\right ) y \mathcal{g}^{(4)}(y)\right],\\ c_3(y ) = & y \left[-3 \left(19 y^3 + 4\right ) \mathcal{g}(y ) + 3y \left(4 - 9 y^3\right ) \mathcal{g}'(y)+4y^2 \left(4 y^3 - 1\right ) \mathcal{g}''(y)+3 y^3 \left(y^3 + 2\right ) \mathcal{g}^{(3)}(y)\right],\\ c_4(y ) = & 3 y^2 \left[\left(2 - 6 y^3\right ) \mathcal{g}(y ) + 2 y\left(y^3 - 1\right ) \mathcal{g}'(y)+ y^2 \left(y^3 + 2\right ) \mathcal{g}''(y)\right],\\ c_5(y ) = & y^3 \left(y^3 + 2\right ) \left(-\mathcal{g}(y)+y \mathcal{g}'(y)\right).\end{aligned}\ ] ] it follows from [ prop - bcircj ] that the function @xmath135 satisfies the differential equations @xmath136 using these equations to eliminate the higher - order derivatives of @xmath135 from the expressions for the coefficients @xmath137 one obtains the claimed expressions . assertion is a reformulation of [ eq - ode1 - 1var ] . 20    m.  abramowitz and i.  a. stegun . , volume  55 of _ national bureau of standards applied mathematics series_. courier corporation , 1964r . p.  bak and k.  sneppen . punctuated equilibrium and criticality in a simple model of evolution . , 71(24):4083 , 1993 . i.  bose and i.  chaudhuri . bacterial evolution and bak - sneppen model . , 12(5):675685 , 2001 . j.  de  boer , b.  derrida , h.  flyvbjerg , a.  d. jackson , and t.  wettig . simple model of self - organized biological evolution . , 73(6):906 , 1994 . r.  donangelo and h.  fort . model for mutation in bacterial populations . , 89(3):038101 , 2002 . n.  eldredge and s.  j. gould . punctuated equilibria : an alternative to phyletic gradualism . in t. j.  m. schopf , editor , _ models in paleobiology _ , chapter  5 , pages 82115 . freeman copper , san francisco , 1972 . h.  flyvbjerg , k.  sneppen , and p.  bak . mean field theory for a simple model of evolution . , 71(24):4087 , 1993 . a.  gillett , r.  meester , and m.  nuyens . bounds for avalanche critical values of the bak  sneppen model . , 12(4):679694 , 2006 . a.  gillett , r.  meester , and p.  van der  wal . maximal avalanches in the bak  sneppen model . , 43(3):840851 , 2006 . m.  grinfeld , p.  a. knight , and a.  r. wade . bak  sneppen - type models and rank - driven processes . , 84(4):041124 , 2011 . m.  grinfeld , philip  a. knight , and a.  r. wade . rank - driven markov processes . , 146(2):378407 , 2012 . h.  j. jensen . , volume  10 . cambridge university press , 1988 . r.  lenski . the e.  coli long - term experimental evolution project site , 2015 . http://myxo.css.msu.edu/ecoli . richard  e lenski and michael travisano . dynamics of adaptation and diversification : a 10,000-generation experiment with bacterial populations . , 91(15):68086814 , 1994 . r.  meester and d.  znamenski . limit behavior of the bak  sneppen evolution model . , 31(4):19862002 , 2003 . r.  meester and d.  znamenski . critical thresholds and the limit distribution in the bak  sneppen model . , 246(1):6386 , 2004 . m.  petkovek , h.  s. wilf , and d.  zeilberger . . a k peters , ltd . , wellesley , ma , 1996 . e.  schlemm . asymptotic fitness distribution in the bak  sneppen model of biological evolution with four species . , 148(2):191203 , 2012 . h.  s. wilf . . a k peters , ltd . , wellesley , ma , third edition , 2006 . in this appendix we collect the values of the coefficients @xmath58 of the functions @xmath54 for @xmath55 . this serves as an illustration of [ lemma - coeff - constant ] and to determine the initial conditions of the various differential equations encountered in the paper .
Abstract:

Answer: the bak  sneppen model is an abstract representation of a biological system that evolves according to the darwinian principles of random mutation and selection . the species in the system are characterized by a numerical fitness value between zero and one . we show that in the case of five species the steady - state fitness distribution can be obtained as a solution to a linear differential equation of order five with hypergeometric coefficients . similar representations for the asymptotic fitness distribution in larger systems may help pave the way towards a resolution of the question of whether or not , in the limit of infinitely many species , the fitness is asymptotically uniformly distributed on the interval @xmath0 $ ] with @xmath1 .


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: gaining control over the external and the internal degrees of freedom of molecules is of great interest for a range of areas in molecular science . one approach , dating as far back as the 1920s  @xcite , relies on the application of inhomogenous ( static or low frequency ) magnetic and electric fields . this has proven most useful for controlling the full three - dimensional motion of molecules , including control of the velocity distribution , bringing molecules to a stand still , storing them for extended periods of times and trapping them in a confined volume ( for a recent review see  @xcite ) . an alternative and complementary approach relies on the use of non - resonant non - ionizing laser fields typically supplied by pulsed lasers . many studies during the last decade have shown that strong non - resonant laser fields can effectively manipulate the external degrees of freedom of isolated gas phase molecules . the manipulation results from laser - induced forces and torques due to the interaction between the induced dipole moment and the laser field itself . examples of manipulation include deflection  @xcite , focusing  @xcite and slowing  @xcite of molecules through the dependence of the non - resonant polarizability interaction on the intensity distribution in a laser focus . likewise , the dependence of the induced dipole interaction on molecular orientation has proven highly useful for controlling the rotation of a variety of molecules  @xcite . in particular , the spatial orientation of molecules can be sharply confined with respect to axes that are fixed in the laboratory . molecular manipulation by induced dipole forces extends beyond the external degrees of freedom and has also been demonstrated for the internal degrees of freedom such as vibrational motion  @xcite in molecular hydrogen . notably , the electrical field from laser pulses can modify energy potential barriers such that photoinduced bond breakage of a small linear molecule is guided to yield a desired final product  @xcite .    in the case of larger molecules control of the lowest frequency vibrational modes attracts special interest since some of these modes correspond to motion along well - defined reaction coordinates separating two conformational minima ( conformers ) . although many molecules contain large number of conformers it is often just two conformers that dominate important chemical properties , for instance chirality . a particularly important example is found in axially chiral molecules  @xcite such as biaryl systems . in these molecules rotation about a single stereogenic carbon - carbon ( c - c ) bond axis changes the molecule from one enantiomer into the opposite enantiomer ( mirror image ) . recently , we demonstrated that the laser - induced nonresonant polarizability interaction can also be used to influence the internal rotation of an axially chiral molecule around the stereogenic c - c bond axis  @xcite . in particular , we showed that by fixing the c - c bond axis of a substituted biphenyl molecule in space , using laser induced alignment by a long laser pulse , it was possible to initiate torsional motion of the two phenyl rings by a short laser pulse polarized perpendicular to the fixed axis . the purpose of the present work is to extend our recent paper  @xcite . in particular , we discuss our theoretical modelling in detail . a ns laser pulse aligns @xcite the stereogenic axis of 3,5-diflouro-3,5-dibromo - biphenyl ( dfdbrbph ) along its polarization . this enables a much shorter fs pulse , which we will refer to as the kick pulse , to both initiate torsional motion and set the molecule into rotation around the fixed axis . the internal as well as external rotational motion is monitored by fs time - resolved coulomb explosion imaging . the inset of fig . [ fig1 ] shows a model of the dfdbrbph molecule with the stereogenic axis marked by red . according to our quantum chemical calculations , the laser - free torsional potential , illustrated by the red dotted curve in fig . [ fig1 ] , has minima at dihedral angles of @xmath0 corresponding to the r@xmath1 and s@xmath1 enantiomeric forms , and where the dihedral angle is the angle between the two phenyl rings . the twisted equilibrium shape is characteristic of biphenyl compounds . the traditional view is that this non - planarity results from competition between stabilization via conjugation of the @xmath2-orbital systems and steric repulsion between the ortho - positioned atoms @xcite . we might mention that alternative explanations have recently been discussed @xcite . our strategy for controlling @xmath3 relies on a transient modification of the field - free potential curve by a fs kick pulse ( fig . [ fig1 ] ) . the modification , caused by the nonresonant polarizability interaction  @xcite converts the initial stationary quantum states , localized near the minima of the torsional potential , into vibrational wave packets , i.e. , coherent superpositions of several quantum states . the temporal evolution of these wave packets gives rise to time - dependent torsional motion , and the grey and the black curves in fig . [ fig1 ] show the calculated expectation values of @xmath3 for a particular set of laser parameters . the stereogenic c - c bond axis is fixed in the laboratory frame by adiabatic alignment utilizing a moderately intense , linearly polarized ns laser pulse @xcite . the ns pulse is intense enough to keep the c - c axis tightly confined , yet weak enough to modify the torsional potential only slightly . the fs kick pulse is applied with its polarization perpendicular to the aligned c - c bond axis to ensure primarily influence on torsion while avoiding other possible vibrational motion . the dfdbrbph compound is prepared as a racemate and the experiment is carried out on isolated molecules at rotational temperatures of a few kelvin . under these conditions we have equal numbers of molecules initially localized in the @xmath4 or 39@xmath5 conformation ( see fig .  [ fig1 ] ) and no thermally induced transitions between the two conformations occur . the intense probe pulse , sent at time @xmath6 with respect to the kick pulse , removes several electrons from the molecules , thereby triggering coulomb explosion into ionic fragments . in particular , the br@xmath7 and f@xmath7 fragment ions recoil in the planes defined by the br- and f - phenyl rings . by recording the velocities of both ion species with two - dimensional ion imaging @xcite , we determine the instantaneous orientation of each of the two phenyl planes at the time of the probe pulse .     between the br- and f - phenyl planes . the asymmetry in the potential is obtained by orienting the molecules ( here with the br - phenyl plane out of the paper ) , 3d aligning them , and by polarizing the kick pulse at an angle of 13@xmath5 with respect to the second most polarizable axis ( smpa ) ( see fig . [ fig2 ] ) . the red dotted curve at large times illustrates the laser - free time - independent torsional potential . for the s@xmath1 enantiomer , starting out with @xmath8 , the time varying potential induces an oscillatory motion ( grey curve ) corresponding to torsion confined within the initial well . by contrast , due to the induced asymmetry between the two wells , the initial r@xmath1 enantiomer is traversing the central torsional barrier , and ends up as an s@xmath1 enantiomer undergoing internal rotation ( black curve ) . the kick pulse has an intensity of @xmath9 w@xmath10@xmath11 and a duration of 1.0 ps ( fwhm ) . the torsional motion may be monitored by fs time - resolved coulomb explosion imaging . the inset shows a model of the dfdbrbph molecule with the sterogenic axis marked by red ( grey ) . ] to account for the motion of the nuclei of the dfdbrbph molecule we make use of two sets of coordinate systems ( see fig . [ fig2 ] ) : a molecular fixed ( mf ) frame attached to the molecule and a laboratory fixed ( lf ) frame specified by the lasers . the mf coordinates are chosen with the @xmath12 axis along the stereogenic axis pointing from the phenyl ring with the bromines towards the phenyl ring with the flourines , and the @xmath13 axis is chosen along the phenyl ring with the bromines . the lf coordinates are chosen with the @xmath14 axis along the polarization direction of the ns pulse and the @xmath15 axis along the kick pulse polarization direction . in agreement with the experimental observations ( cf . [ sec : results ] ) , we assume that the stereogenic axis of the dfdbrbph molecule is aligned along the @xmath14 axis .     and the laboratory fixed @xmath16 coordinate axes . the dotted line indicates the second most polarizable axis ( smpa ) , which is located 11@xmath5 from the br - phenyl ring and 28@xmath5 from the f - phenyl ring . ] also , we neglect all normal modes , except the lowest one , which corresponds to torsion ( see sec . [ subsec : qmcalc ] for details on quantum chemical calculations , justifying , e.g. , this assumption ) . in this simplified situation and in the absence of the kick pulse the task is reduced to describing the coupled rotations of the two ( rigid ) phenyl rings of the molecule as given in the lf frame by the hamiltonian [ atomic units , @xmath17 , are used throughout unless indicated otherwise ] @xmath18 where @xmath19 , is the angle of the @xmath20 phenyl ring with respect to the kick pulse polarization axis , @xmath21 the moment of inertia for rotation of the @xmath20 phenyl ring around the stereogenic axis ( @xmath22 , @xmath23 ) and @xmath24 is the laser - free torsional potential as obtained from quantum chemical calculations ( fig .  3 ) . by changing the coordinates to the dihedral angle @xmath25 between the two phenyl rings and the weighted azimuthal angle @xmath26 , characterizing the rotation of the molecule , with @xmath27 , we obtain @xmath28 here @xmath29 is the total moment of inertia of the molecule for rotation around the stereogenic axis and @xmath30 is a relative moment of inertia for the two phenyl rings around the axis . a full rotation of either phenyl ring leaves us with the same molecule implying @xmath31-periodic boundary conditions of the eigenfunctions of @xmath32 from eq .  , i.e. , @xmath33 , with @xmath34 and @xmath35 integers . we shall assume that this property also applies to @xmath36 and @xmath3 , so that we simply need to consider eigenfunctions @xmath37 of eq . that separates into rotation of the molecule as described by the @xmath31-periodic function @xmath38 and torsion accounted for by the @xmath31-periodic function @xmath39 . [ rigorously , the bounds for @xmath36 depend on @xmath3  @xcite . at the time scales of interest , however , a molecule generally only makes a small fraction of a full rotation so that the effect of @xmath36-@xmath3 coupling is negligible . ] the separation is physically motivated by considering the energy scales related to rotation and torsion . the energy scale of the prior is given by @xmath40ev . for torsion , on the other hand , the relevant energy is determined by the torsional potential , and a harmonic approximation of the potential near the minimum at @xmath41 yields a frequency corresponding to the energy @xmath42  mev . using that the period of motion is of the order planck s constant divided by the energy , we therefore see that the molecule rotates ( @xmath36 change ) on a nanosecond time scale , whereas the torsion ( @xmath3 change ) is of picosecond duration . the small energy separation of the rotational levels compared to that of the torsional levels implies that many rotational states will be occupied in thermal equillibrium . consequently , we expect the @xmath36 dynamics to behave classically ( see also fig . [ fig7 ] ) , and we only treat the torsional dynamics quantum mechanically . the torsion is described by the stationary schrdinger equation @xmath43 with @xmath44 denoting the energy eigenstates . to solve this equation we expand the hamiltonian onto an orthonormal basis of @xmath31-periodic functions , which we diagonalize . the eigenstates below the torsional barriers exhibit an almost exact fourfold degeneracy , and linear combinations of the approximately degenerate states will also be stable at the time scales of the experiment . in particular , we can produce such stable states that are localized in the wells of the torsional potential . we will denote these states by @xmath45 , with @xmath46 and @xmath47 the smallest @xmath48 of the four degenerate states . in fig . [ fig3 ] we show the first four energy eigenstates and the corresponding localized states along with the torsional potential . mev above the minimum of the torsional potential . from linear combinations of these we obtain the corresponding four localized states shown in ( b ) . the torsional potential is indicated with a dashed line . ] the experiment is generally carried out on a gas of molecules in a mixed state given by the density matrix @xmath49 with the @xmath50 s being weight factors that sum to unity , and where the latter sum includes the states above the torsional barrier . the kick pulse introduces dynamics that we determine from the alignment potential created by the pulse . by definition @xmath51 for the kick pulse , and we use a gaussian envelope @xmath52 ^ 2/2\}$ ] corresponding to an intensity full width at half maximum ( fwhm ) of @xmath53 . we calculate the polarizability tensor in the mf frame as a function of the dihedral angle ( see table  [ tab : table1 ] ) and consequently transform into the lf frame by the application of directional cosine matrices  @xcite . we then arrive at the interaction @xmath54.\end{aligned}\ ] ]    [ cols="<,<,>,>,>,>,>,>,>",options="header " , ] to a round - bottom flask equipped with a magnetic stirring bar under argon we add 3,5-difluorophenyl - boronic acid ( 474 mg , 3.0 mmol ) , 1,3,5-tribromobenzene ( 1134 mg , 3.6 mmol ) , and pd(pph@xmath55)@xmath56 ( 69.6 mg , 2 mol@xmath57 ) . the flask is then evacuated and backfilled with argon twice . to this mixture of solids we add toluene ( 18 ml , degassed ) , etoh ( 6 ml , degassed ) , and 2 m aq . na@xmath58co@xmath55 ( 3 ml , degassed ) . we next heat the mixture to 90@xmath5c under argon for 5.5 h and cool it to room temperature . the mixture is diluted with water and extracted three times with et@xmath58o . we dry the combined organic fractions with mgso@xmath56 , filter and concentrate to give approximately 1.4 g of a crude mixture . based on nmr we find that this mixture contains both the desired mono - suzuki coupling product @xcite , double - coupled products , and a homo - coupling product from the boronic acid . careful column chromatography on sio@xmath58 eluting with n - hexane and c - hexane finally allows us to isolate 182 mg 3,5-dibromo-3,5-difluorobiphenyl ( dfdbrbph ) as a white amorphous solid ( contains @xmath593@xmath57 of the homocoupling product dfdbrbph ) . nmr spectroscopic data ( see also fig . [ fig4 ] ) : @xmath60h nmr @xmath61 ( @xmath62 mhz , cdcl@xmath55 ) 7.69 ( t , j 1.7 hz , 1h ) , 7.61 ( d , j 1.7 hz , 2h ) , 7.05 ( m , 2h ) , 6.85 ( m , 1h ) . @xmath63c nmr @xmath61 ( 100 mhz , cdcl@xmath55 ) 163.5 ( 2c ) , 142.3 , 141.5 , 133.8 , 128.9 ( 2c ) , 123.5 ( 2c ) , 110.1 ( 2c ) , 103.8 . elemental analysis , calculated : c 41.42@xmath57 , h 1.74@xmath57 , br 45.92@xmath57 ; found : c 41.20@xmath57 , 1.66@xmath57 , br 45.90@xmath57 . ( color online ) ( a ) , synthesis of dfdbrbph by pd - catalyzed cross - coupling . ( b)-(d ) , nmr - spectra of synthetic dfdbrbph . integrated signals are indicated by green ( grey ) curves , and the corresponding numbers are situated directly below the respective curves . ] a schematic of the experimental setup is shown in fig . less than hundred milligram of solid dfdbrbph is heated to 130@xmath5c and expanded in 90 bar of he into vacuum using a pulsed even - lavie valve @xcite . the resulting supersonic molecular beam is skimmed and subsequently crossed by three , focused , pulsed laser beams at 90@xmath5 . the first laser beam originates from a q - switched nd : yag laser and consists of 9-ns - long pulses at 1064 nm . it is focused to a spotsize @xmath64 @xmath65 35 @xmath66 m yielding a peak intensity of @xmath67 w/@xmath11 . the purpose of these pulses is to align the stereogenic axis of the dfdbrbph molecules . the second laser beam consists of 0.7-ps - long kick pulses ( @xmath68 = 800 nm ) , obtained by passing part of the output from an amplified ti : sapphire fs laser system through a grating stretcher . the spotsize of @xmath64 = 43 @xmath66 m yields a peak intensity of @xmath69 w/@xmath11 . these kick pulses are used to induce torsion and rotation of the molecules . the third laser beam consists of 25-fs - long pulses , which we achieve by passing another part of the output from the ti : sapphire system through a hollow wave guide compressor setup . the pulses have a spotsize @xmath64 = 25 @xmath66 m and a peak intensity of @xmath70 w/@xmath11 . these probe pulses coulomb explode the molecules and the produced ions are extracted by a weak static electric field , in a velocity imaging geometry , and projected onto a two - dimensional detector , consisting of a micro channel plate ( mcp ) detector backed by a phosphor screen . we detect the f@xmath7 and br@xmath7 ions separately by time gating of the mcp detector around the respective arrival times of the ions . we record the ion images on the phosphor screen by a ccd camera and determine the coordinates of each individual ion hit . figure  [ fig6](a ) displays the ion images of f@xmath7 and br@xmath7 at selected probe times . ( color online ) schematic representation of the experimental setup , showing the vacuum system and the three laser beam paths . the nd : yag ( alignment ) pulses and the probe pulses are polarized vertically , i.e. , perpendicular to the detector plane , whereas the fs kick pulse is polarized in the plane of the detector . ] first , we establish that the ns pulse is able to hold the c - c bond axis of dfdbrbph along its polarization direction . at @xmath71 ps the f@xmath7 image , displayed in fig . [ fig6](a ) , is almost circularly symmetric and the small deviation from circular symmetry is explained by noting that the kick pulse has a finite value at @xmath72 ps . the absence of ions in the innermost region is only compatible with the c - c bond axis being aligned perpendicular to the detector plane and the f - phenyl planes uniformly distributed around the c - c bond axis . the corresponding br@xmath7 image is also essentially circularly symmetric although the hollow - out of the center is not as pronounced as for f@xmath7 . next , we investigate the effect of the kick pulse . already at @xmath73 ps the circular symmetry is broken and both f@xmath7 and br@xmath7 ions start to localize around the polarization direction of the kick pulse . the f@xmath7 ions remain within the radial range of the ring structure , which shows that the kick pulse does not cause any significant distortion of the c - c bond axis alignment but rather initiates an overall rotation of the molecule around this axis . the latter is expected since the torque imparted by the kick pulse will force the smpa ( fig .  2 ) to align along the kick pulse polarization after a delay determined by the kick strength . at @xmath74 ps the localization sharpens for both ions species , but whereas the br@xmath7 ions are simply more localized around the polarization axis the f@xmath7 ion distribution exhibits a four dot structure . this behaviour is compatible with alignment of the smpa along the kick pulse polarization . in the limit of perfect smpa alignment , the br@xmath7 ions would appear as two pairs of strongly localized regions on the detector corresponding to molecules with their br - phenyl plane located either 11@xmath5 clockwise or counterclockwise to the smpa . each pair of br@xmath7 ions would then be angularly separated by 22@xmath5 . in practise and consistent with theory ( see fig .  [ fig8 ] ) the smpa alignment is not strong enough to resolve the two br@xmath7 regions . it is , however , sharp enough to resolve the two pairs of f@xmath7 ion hits due to the much larger offset ( 28@xmath5 ) of the f - phenyl plane from the smpa . at @xmath75 ps the br@xmath7 distribution has localized further showing that the br - phenyl plane has rotated into stronger alignment with the kick pulse polarization . if the torsion was unaffected , i.e. , @xmath3 remained unchanged , the f@xmath7 ion image at @xmath75 ps should exhibit a distinct four - dot structure similar to the image at @xmath76 ps but with a larger angular splitting between the ion regions in each of the two pairs . clearly , the four - dot structure at 2.47 ps is significantly blurred compared to the case at 1.47 ps . thus , we conclude that the kick pulse not only sets the molecule into controlled rotation around the c - c axis , it also initiates torsional motion .    beyond 3.8 ps the angular confinement of both the f@xmath7 and the br@xmath7 ions is gradually lost due to continued rotation around the c - c axis with dihedral dynamics imposed , and the images for both ion species eventually regain their circularly symmetric form ( see images at 12.15 ps in fig .  [ fig6](a ) ) identical to the pre - kick pulse images at -0.87 ps . further insight into the effect of the kick pulse is obtained by analyzing the angular distribution of the f@xmath7 ions as a function of @xmath6 . an example of the f@xmath7 angular distribution at @xmath6 = 1.47 ps is shown in fig . [ fig6](b ) . figure  [ fig6](c ) displays the average angle between the f - phenyl rings and the kick pulse polarization as a function of @xmath6 in the time interval where a clear four - peak structure is visible in the angular distributions . the increase from 22.5@xmath5 at 1.47 ps to 26.5@xmath5 at 3.8 ps shows that the f - phenyl plane gradually moves away from the kick pulse polarization due to the overall rotation of the molecule around the c - c bond axis , and we ascribe the recurrent dips to a periodical variation in @xmath3 . we estimate the period to be @xmath77 ps and the amplitude to @xmath78 for this oscillation .    before presenting the theoretical results on dfdbrbph , we want to confirm the validity of the classical treatment of rotation . to this end we apply the model to a test case , where the simpler 3,5-diflouroiodobenzene ( dfib ) has its most polarizable axis fixed by a ns laser pulse and is then set into rotation about this axis by an orthogonally polarized femtosecond pulse . here , there is no torsion and the rotational dynamics have been observed directly recently  @xcite . in fig . [ fig7 ] , we compare the theoretical and experimental results and find good agreement .        , 1.47 and 2.47 ps . ( c ) expectation value of the dihedral angle for a molecule starting out with the smpa aligned along the kick pulse polarization . the kick pulse is as in fig . [ fig6 ] . ] figure  [ fig8 ] shows the results of a calculation with laser parameters identical to the experimental values and an initial rotational temperature of 0 k. prior to the kick pulse the angular distributions of the br- and f - phenyl rings are isotropic as in the experiment . maximum alignment of the smpa occurs at 1.3 ps and the confinement of the f - phenyl rings at a large angle with respect to the kick pulse polarization ( cf . middle panel , fig . [ fig8](a ) ) explains the distinct four - dot structure observed at @xmath74 ps in the experimental f@xmath7 ion image ( fig . [ fig6](a ) ) . also , at @xmath74 ps the confinement of the br - phenyl rings at a small angle with respect to the kick pulse polarization predicts a much less distinct , if any , four - dot structure in good agreement with the br@xmath7 ion image . at @xmath75 ps the angular localization of the f - phenyl rings has broadened ( right panel , fig .  [ fig8](a ) ) and a blurred four - dot structure is seen , consistent with the experimental result at @xmath75 ps . the distribution of the br - phenyl rings is also broadened ( right panel , fig . [ fig8](b ) ) , but remains localized around the kick pulse polarization fully consistent with the br@xmath7 ion distribution , recorded at 2.47 ps . the theoretical value @xmath79 exhibits oscillations ( fig . [ fig8](c ) ) with a period of @xmath80 ps and amplitude of @xmath81 . the period agrees well with the experimental value ( @xmath82 ps ) , and smaller modulation in @xmath83 is expected in the experiment ( @xmath84 ) partly since here the smpa is not pre - aligned and partly due to deviations from the stated laser peak intensity . the behavior is ascribed to a wave packet of normal modes in the torsional double well potential ( fig . [ fig3](b ) ) for a molecule starting out with the smpa aligned . the qualitative agreement of figs . [ fig6](b ) and [ fig8](c ) corroborates the interpretation of the kick pulse inducing time - dependent torsional motion . we will now discuss some perspectives of the demonstrated laser control of torsion . for one thing , we suggest a time - resolved study of de - racemization  @xcite , where one enantiomer is steered into its mirror form . to do this we break the inversion symmetry of the c - c bond axis by orienting each molecule . this allows us to discriminate between the two enantiomeric forms . theory  @xcite and experiment  @xcite show that orientation can be added to 3d alignment by combining the ns alignment pulse with a static electric field . next , the excitation of torsional motion is optimized by eliminating initial rotation around the c - c bond through alignment of the smpa ( fig . 2 ) prior to the kick pulse using an elliptically rather than a linearly polarized ns pulse @xcite . finally , the interaction strength between the molecule and the kick pulse needs to be increased either through higher intensity , a longer kick pulse or by trains of synchronized kick pulses  @xcite .     or ( b ) an s@xmath1 enantiomer . initially , the molecule is 3d oriented with the br - phenyl end pointing out of the paper and the smpa aligned at an angle of 13@xmath5 with respect to the kick pulse polarization . the kick pulse triggering the torsional motion has a peak intensity of @xmath9 w/@xmath11 and duration ( fwhm ) of @xmath85 ps . the torsional potential is scaled down by 1/4 rather than by increasing the kick strength . in practise reduction of the torsional barrier is possible , for instance , by modifying the aromatic rings or by using halogen substituted biphenylacetylene . ] assuming initial orientation and confinement of the smpa , we have calculated the @xmath3 dynamics for both conformations the dfdbrbph molecule with a reduced torsional barrier ( see caption to fig . [ fig9 ] ) . the results are shown in fig . [ fig9 ] , and clearly , the present method would allow for a time - resolved study of the transition from one enantiomer into the other . a quantitative analysis of the efficiency of the process shows that after the pulse @xmath86 of the molecules starting out as r@xmath1 has changed into s@xmath1 enantiomers , whereas only @xmath87 of the s@xmath1 enantiomers changed into r@xmath1 . the inverse process causing an excess of r@xmath1 enantiomers , is achieved simply by inverting the orientation of the molecules . we believe that the results presented in this paper also fuel the field of molecular junctions  @xcite with new possibilities . the conductivity of molecules like biphenyl will depend on the dihedral angle  @xcite , and as such this type of molecules , if placed between two leads , can function as a molecular junction  @xcite that may be used to control the charge flow and act as , e.g. , a switch  @xcite . importantly , we can realize quantum control of such a system by means of kick pulse schemes and in this way tailor very specific torsional wave packets  @xcite to dictate the detailed time - dependent dihedral motion and thereby the current flow through the molecular junction . seen through the light of controllability , molecular junctions based on laser controlled torsion complements the available schemes , such as the mechanical break junctions  @xcite and the previously suggested resonantly light driven molecular junctions  @xcite . we presented a fs time - resolved study on laser controlled torsion of axially chiral molecules . experimentally , the symmetry axis of dfdbrbph molecules was held fixed in space by a long alignment laser pulse while a fs kick pulse , polarized perpendicular to the fixed axis , was applied . the torsion as well as the overall rotation of the molecule was monitored with fs time resolution using a delayed probe pulse . to substantiate the observations , we developed an original theoretical model , which produces results in good agreement with the experiment . the model divides the dynamics of the molecule into a the rotation with classical behavior and quantum mechanical torsional dynamics . as such the model provides a transparent physical intepretation of essential features of laser controlled torsion of dfdbrbph . based on the theoretical model , we further discussed the extension of the experiment to a setup that will facilitate a fs time - resolved study of de - racemization . finally , we pointed to perspectives of using molecules like dfdbrbph in molecular junctions in order to realize laser controlled charge flow . this work was supported by the carlsberg foundation , the lundbeck foundation , danish natural research foundation , the danish natural science research council , and the danish research agency ( grant no . 217 - 05 - 0081 ) .
Abstract:

Answer: it is demonstrated that strong laser pulses can introduce torsional motion in the axially chiral molecule 3,5-diflouro-3,5-dibromo - biphenyl ( dfdbrbph ) . a nanosecond laser pulse spatially aligns the stereogenic carbon - carbon ( c - c ) bond axis allowing a perpendicularly polarized , intense femtosecond pulse to initiate torsional motion accompanied by a rotation about the fixed axis . we monitor the induced motion by femtosecond time - resolved coulomb explosion imaging . our theoretical analysis corroborates the experimental findings and on the basis of these results we discuss future applications of laser induced torsion , viz . , time - resolved studies of de - racemization and laser controlled molecular junctions based on molecules with torsion .


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: boundary value problems ( bvps ) on infinite intervals arise in several branches of science . the classical numerical treatment of these problems consists in replacing the original problem by one defined on a finite interval , say @xmath0 $ ] where @xmath1 is a truncated boundary . the oldest and simplest approach is to replace the boundary conditions at infinity by the same conditions at the value chosen values as the truncated boundary . this approach was used , for instance , by goldstein @xcite and by howarth @xcite in 1938 to get and tabulate the numerical solution of the blasius problem . however , in order to achieve an accurate solution , a comparison of numerical results obtained for several values of the truncated boundaries is necessary as suggested by fox @xcite , or by collatz @xcite . moreover , in some cases accurate solutions can be found only by using very large values of the truncated boundary . this is , for instance , the case of the fourth branch of the von karman swirling flows , where values of @xmath2 up to @xmath3 were used by lentini and keller @xcite .    to overcome the mentioned difficulties of the classical approach described above , lentini and keller @xcite and de hoog and weiss @xcite suggested to apply asymptotic boundary conditions ( abcs ) at the ; see also the theoretical work of markowich @xcite . those abcs have to be derived by a preliminary asymptotic analysis involving a jacobian matrix of the right - hand side of the governing equations evaluated at infinity . the main idea of this abcs approach is to project the solution into the manifold of bounded solutions . by this approach more accurate numerical solutions can be found than those obtained by the classical approach with the same values of the , because the imposed conditions are obtained from the asymptotic behaviour of the solution . however , we should note that for nonlinear problems highly nonlinear abcs may result . moreover , it has been noticed by j. r. ockendon that unfortunately the analysis is heavy and relies on much previous work ,  see _ math . _ 84c:34201 . on the other hand , starting with the work by beyn @xcite , the abcs approach has been applied successfully to connecting orbits  problems . connecting orbits are of interest in the study of dynamical systems as well as of travelling wave solutions of partial differential equations of parabolic type . however , a truncated boundary allowing for a satisfactory accuracy of the numerical solution has to be determined by trial , and this seems to be the weakest point of the classical approach . hence , a priori definition of the truncated boundary was indicated by lentini and keller @xcite as an important area of research . a free boundary formulation for the numerical solution of bvps on infinite intervals was proposed in @xcite . in this approach the truncated boundary can be identified as an unknown free boundary that has to be determined as part of the solution . as a consequence , the free boundary approach overcomes the need for a priori definition of the truncated boundary . this new approach has been applied to : the blasius problem @xcite , the falkner - skan equation with relevant boundary conditions @xcite , a model describing the flow of an incompressible fluid over a slender parabola of revolution @xcite , and a model describing the deflection of a semi - infinite pile embedded in soft soil @xcite . an application of the free boundary approach to a homoclinic orbit problem can be found in @xcite . moreover , a possible way to extend the free boundary formulation to problems governed by parabolic partial differential equations was proposed in @xcite . it might seem that in order to face numerically a bvp defined on an infinite interval , we have to reformulate it in a way or another . however , recently , we have found that it is also possible to apply directly to the given bvp a finite difference method defined on a quasi - uniform grid . to this end it is necessary to derive special finite difference formulae on the grid involving the given boundary conditions at infinity , but the last grid point value ( infinity ) is not required ; see @xcite for details .    in this paper , for an ocean circulation model , we report a comparison of numerical results obtained by the classical truncated boundary approach with a shooting method , those found by our free boundary approach with a finite difference method , and the ones obtained by a finite difference method with a quasi - uniform grid . a steady - state wind - driven ocean circulation model can be introduced , see ierley and ruehr @xcite , by considering the barotropic vorticity equation @xmath4 in a region defined by @xmath5 $ ] and @xmath6 $ ] with the following boundary conditions @xmath7 and either @xmath8 known as rigid  or no - slip boundary conditions , or @xmath9 known as slippery  or stress - free boundary conditions . equation ( [ pde - model ] ) is written in a non - dimensional form ; @xmath10 is the stream function ; to fix a reference coordinate system is fixed with the @xmath11 axis directed to the east and the @xmath12 axis directed to the north ; @xmath13 is the jacobian of the functions @xmath14 and @xmath15 with respect to @xmath11 and @xmath12 , @xmath16 is the laplacian operator on the @xmath17 plane ; the square @xmath18\times[-1,1]$ ] models a region of the subtropical gyre formation . here the jacobian represents nonlinear advection and the laplacian the viscous drag . we assume that the curl of the wind stress in the region can be approximated by @xmath19 ; @xmath20 and @xmath21 , are non - dimensional parameters characterizing the widths of inertial and viscous boundary layers , respectively . we use impermeability and no - slip conditions ( [ pde - bcs2 ] ) at the coasts and impermeability and slip conditions ( [ pde - bcs3 ] ) at the fluid boundaries . we consider a particular solution to ( [ pde - model ] ) of the form @xmath22 relation ( [ eq : solution ] ) represents the first term in the expansion of a solution of ( [ pde - model ] ) with respect to @xmath12 near boundaries of the region : at @xmath23 and at @xmath24 . substituting ( [ eq : solution ] ) into ( [ pde - model ] ) , using a taylor series expansion near @xmath23 of the wind - stress term and assuming that a steady boundary - layer type solution exists , we obtain the equation for the boundary layer at the western coast , i.e. at @xmath25 , @xmath26 the parameters involved can be reduced to one if we define @xmath27 and introduce the new independent variable @xmath28    the limit of vanishing viscosity ( small values of @xmath21 ) is of particular interest . indeed , the parameter @xmath20 is also small , of the order of @xmath29 . therefore , in terms of the new independent variable @xmath30 , far from the boundaries for asymptotically matching the interior solution @xmath31 , taken of the following form @xmath32 we have to require that @xmath33    our fourth order ordinary differential equation ( [ cmodel4 ] ) can be integrated once , using zero boundary conditions at infinity for the second and third derivative of @xmath34 , to give @xmath35 + u - 1 \ , \qquad \xi \in \left[0,\infty\right ) \ .\ ] ] the boundary conditions follow from ( [ pde - bcs1])-([pde - bcs3 ] ) . in particular , we can have no - slip ( or rigid ) boundary data @xmath36 or stress - free ( or slippery ) boundary conditions @xmath37 therefore , we get the two point bvp defined on an unbounded domain that has been investigated by ierley and ruehr @xcite , mallier @xcite or sheremet et al . @xcite . the parameter @xmath15 in ( [ cmodel ] ) can be used as a measure of the strength of the nonlinearity . in fact , for @xmath38 we get the simple linear model formulated by munk @xcite . ierley and ruehr @xcite discovered an analytical approximation for the relation between the missing initial condition and the parameter @xmath15 . in particular , for rigid conditions they proposed the relation @xmath39 and for slippery conditions they reported the relation @xmath40 the approximations provided by ( [ rigid ] ) and ( [ slippery ] ) will be used for comparison with the corresponding numerical results . in this section we present the numerical methods used in order to solve the ocean model ( [ cmodel ] ) . as a first step we rewrite the ocean equation in ( [ cmodel ] ) as a first order system @xmath41 \label{p_qu } \\[-1.5ex ] & & { \bf g } \left ( { \bf u}(0 ) , { \bf u } ( \infty ) \right ) = { \bf 0 } \ ,   \nonumber\end{aligned}\ ] ] by setting @xmath42 in this way the original bvp ( [ cmodel ] ) specializes into @xmath43 that is , @xmath44 with @xmath45 or @xmath46 in ( [ p_qu ] ) . in the following , in order to set a specific test problem , we consider the ocean model with @xmath47 . it is simple to describe a classical shooting method . we set a value for the truncated boundary @xmath48 . then , we guess a value for the missing initial condition , call it @xmath49 , and integrate the given problem as an initial value problem ( ivp ) . this defines , implicitly , a nonlinear equation @xmath50 , where @xmath49 can be considered as a parameter . in order to get the correct value of @xmath49 , it is possible to apply a root - finder method . we have found that the secant method is particularly suitable to this purpose . the same shooting method can be applied by using the newton s . this requires a more complex treatment involving a system of six differential equations . in both cases , secant or newton s root finder , we set @xmath51 . we consider first the results obtained by applying the secant method for the system of three equations ( [ eq : cmodel : system ] ) . for the problem ( [ cmodel ] ) with no - slip boundary conditions ( [ bc1 ] ) , setting @xmath52 and @xmath53 , we found the missing initial condition @xmath54 by @xmath55 iterations . in the second test , with slip boundary conditions ( [ bc2 ] ) , setting @xmath56 and @xmath57 , we obtained the missing initial condition @xmath58 with @xmath59 iterations .    on condition to provide a value for @xmath60 , it is also possible to implement the newton s iterations . this can be done by differentiating with respect to @xmath49 the governing system . of course , we end up to solve a system of six equations , namely the system in ( [ p_qu ] ) with @xmath61 & \\[-1ex ] & { \bf f } ( \xi , { \bf u } ) = \left(u_2,u_3 , b ( u_2 ^ 2 - u_1 u_3 ) + u_1 -1 , { \ds \frac{\partial f_1}{\partial \beta},\frac{\partial f_2}{\partial \beta } , \frac{\partial f_3}{\partial \beta } } \right)^t \ , \nonumber\end{aligned}\ ] ] with initial conditions , for the no - slip case @xmath62 and for the slip case @xmath63 where @xmath64 , @xmath65 , @xmath66 , and in both cases @xmath67 in both cases @xmath68 .    let us report now the numerical results given by the newton s method with the system of six equations ( [ eq : modnewton ] ) with ( [ eq : modnewton : ics : noslip ] ) or with ( [ eq : modnewton : ics : slip ] ) . for the problem ( [ cmodel ] ) with no - slip boundary conditions ( [ bc1 ] ) , setting @xmath52 , we obtained the missing initial condition @xmath69 by @xmath70 iterations . in the second test , with slip boundary conditions ( [ bc2 ] ) , setting @xmath56 , we found the missing initial condition @xmath71 with @xmath72 iterations . as well known , the secant and newton s methods are convergent provided that we use initial iterates sufficiently close to the root . moreover , the convergence is super - linear and the order of convergence is @xmath73 for the secant method and @xmath74 for newton s one . we notice that , for the root of nonlinear equations , as reported by gautschi @xcite the secant method has an efficiency index higher than that of the newton s method . on the other hand , the newton s method might be preferable since it requires only one initial guess . the computational cost of these two shooting methods is given in table  [ tab : ccost ] . as it is easily seen the shooting - newton method is the less demanding of the two and we can conclude that in relation of the numerical solution of our bvp the newton s method is more efficient than the secant method .    llrcrrrr + bcs & root - finder & @xmath75 & @xmath76 & steps & rejections & _ eval _ & iterations + [ 1.5ex ] no - slip & secant & 1 & 2 & 109111 & 142 & 327771 & 12 + no - slip & newton & 1 & & 1489 & 79 & 4711 & 7 + slip & secant & 0.8 & 1 & 28461 & 208 & 86020 & 13 + slip & newton & 0.8 & & 6263 & 114 & 19139 & 8 +    [ tab : ccost ]    figure  [ ir_shoot_ns ] and [ ir_shoot_s ] display the numerical results obtained by newton s iterations and related , respectively , to model ( [ cmodel ] ) with boundary conditions ( [ bc1 ] ) and ( [ bc2 ] ) ( see also @xcite ) .    [ ] [ ] @xmath77 [ ] [ ] @xmath78 [ ] [ ] @xmath79 [ ] [ ] @xmath80 [ ] [ ] @xmath81 ) with rigid boundary condition ( [ bc1 ] ) by the shooting newton s method . , title="fig : " ]    [ ] [ ] @xmath82 [ ] [ ] @xmath83 [ ] [ ] @xmath79 [ ] [ ] @xmath80 [ ] [ ] @xmath81 ) with slip boundary condition ( [ bc2 ] ) by the shooting newton s method . , title="fig : " ]    in both tests , we have used the following termination criterion @xmath84 with @xmath85 . moreover , the numerical solutions of the ivps were obtained by the * ode23 * solver , from the matlab ode suite written by samphine and reichelt @xcite , with the accuracy and adaptivity parameters defined by default .    as far as the shooting method is concerned , it may be not suitable even for the simple truncated boundary approach . in fact , it would be possible , when a large step size is used , or always for some models , that one obtains floating - point overflows in the calculations . this is exactly the reason for the introduction of the more complex multiple shooting method ( see also @xcite ) .      in order to introduce a free boundary formulation for our problem , we replace the far boundary condition by two boundary conditions at the free boundary @xmath86 @xmath87 where @xmath88 can be considered as a truncated boundary . then we rewrite the resulting free bvp in standard form ( see ascher and russell @xcite ) , defining @xmath89 and using the new independent variable @xmath90 in general , we end up with a bvp belonging to the general class : @xmath91 \ , \nonumber \\[-1.5ex ] \\[-1.5ex ] & & { \bf g } \left ( { \bf u}(0 ) , { \bf u } ( 1 ) \right ) = { \bf 0 } \ , \nonumber\end{aligned}\ ] ] where @xmath92 where , in our case , @xmath93 . in order to simplify notation in ( [ pp])-([ppwhat ] ) and in the following , we omit the dependence of @xmath94 and @xmath95 on @xmath96 .    in order to solve the resulting problem we apply a relaxation method . let us introduce a mesh of points @xmath97 , @xmath98 , for @xmath99 , of uniform spacing @xmath100 and naturally @xmath101 . we denote by the @xmath102dimensional vector @xmath103 the numerical approximation to the solution @xmath104 of ( [ pp ] ) at the points of the mesh , that is for @xmath105 . keller s box scheme for ( [ pp ] ) can be written as follows : @xmath106 \label{boxs } \\[-1.5ex ] & { \bf g } \left ( { \bf v}_0 , { \bf v}_j \right ) = { \bf 0 } \ ,   \nonumber\end{aligned}\ ] ] where @xmath107 . it is evident that ( [ boxs ] ) is a nonlinear system with respect to the unknown @xmath108dimensional vector @xmath109 . following keller , the classical newton s method , along with a suitable termination criterion , is applied to solve ( [ boxs ] ) . let us recall now the main properties of the box scheme proved by keller in the main theorem of @xcite . under the assumption that @xmath110 and @xmath111 are sufficiently smooth , each isolated solution of ( [ pp ] ) is approximated by a difference solution of ( [ boxs ] ) which can be computed by newton s method , provided that a sufficiently fine mesh and an accurate initial guess for the newton s method are used . as far as the accuracy issue is concerned , the truncation error has an asymptotic expansion in powers of @xmath112 .    for the newton s method we used the simple termination criterion @xmath113 , @xmath114 and @xmath115 , is the difference between two successive iterate components and @xmath116 is a fixed tolerance . the key point for the numerical solution of the nonlinear system is that newton s method converges only locally . therefore , some preliminary numerical experiments may be helpful and worth of consideration . however , for the results reported , our initial guess to start the iterations was always as follows : @xmath117    in figures  [ ir_noslip ] and [ ir_slip ] and in tables  [ tt1 ] and  [ tt2 ] we report some of the numerical results , obtained with the free boundary approach , related to different values of @xmath96 obtained by setting @xmath118 and @xmath119 . here and in the following @xmath120 is the standard notation for @xmath121 in simple precision arithmetic . [ ] [ ] @xmath82 [ l][]@xmath83 [ l][]@xmath79 [ ] [ ] @xmath80 [ ] [ ] @xmath81 [ ] [ ] @xmath122 [ ] [ ] @xmath123 ) with no slip conditions ( [ bc1 ] ) by the free boundary approach.,title="fig : " ] [ ] [ ] @xmath124 [ ] [ ] @xmath125 ) with no slip conditions ( [ bc1 ] ) by the free boundary approach.,title="fig : " ] [ ] [ ] @xmath126 [ ] [ ] @xmath127 ) with no slip conditions ( [ bc1 ] ) by the free boundary approach.,title="fig : " ] [ ] [ ] @xmath128 [ ] [ ] @xmath129 ) with no slip conditions ( [ bc1 ] ) by the free boundary approach.,title="fig : " ]    [ ] [ ] @xmath82 [ l][]@xmath83 [ l][]@xmath79 [ ] [ ] @xmath80 [ ] [ ] @xmath81 [ ] [ ] @xmath130 [ ] [ ] @xmath123 ) with slip conditions ( [ bc2 ] ) by the free boundary approach.,title="fig : " ] [ ] [ ] @xmath131 [ ] [ ] @xmath125 ) with slip conditions ( [ bc2 ] ) by the free boundary approach.,title="fig : " ] [ ] [ ] @xmath132 [ ] [ ] @xmath127 ) with slip conditions ( [ bc2 ] ) by the free boundary approach.,title="fig : " ] [ ] [ ] @xmath133 [ ] [ ] @xmath129 ) with slip conditions ( [ bc2 ] ) by the free boundary approach.,title="fig : " ]    cr@.lcc + @xmath134 & @xmath135 & iter & @xmath136 + [ 1.5ex ] 1@xmath1372 & 6 & 485761 & 7 & 0.826184 + 1@xmath1373 & 8 & 792991 & 8 & 0.826141 + 1@xmath1374 & 11 & 098635 & 10 & 0.826141 + 1@xmath1375 & 13 & 402219 & 11 & 0.826142 +    [ tt1 ]    cr@.lcc + [ -1.5ex ] @xmath134 & @xmath135 & iter & @xmath138 + [ 1.5ex ] 1@xmath1372 & 5 & 828307 & 7 & 0.528970 + 1@xmath1373 & 8 & 132813 & 8 & 0.528922 + 1@xmath1374 & 10 & 437875 & 9 & 0.528921 + 1@xmath1375 & 12 & 741323 & 11 & 0.528921 +    [ tt2 ]    we have several reasons to consider a free boundary formulation more effective than the simple truncated boundary approach . first of all , the free boundary conditions are less suitable for the applications of a shooting method forcing us to use a more suitable relaxation ( finite difference ) method . furthermore , we know a priori that the free boundary has to be an increasing function of @xmath134 . moreover , @xmath134 itself can be consider as a continuation parameter . this means that the numerical results obtained for a value @xmath134 can be used as the initial guess for newton s method for the next value of @xmath134 . consequently , except for the first value of @xmath134 , for both kind of boundary conditions the number of iterations of the relaxation method can be reduced to 7 , 6 , 6 , 6 . let us consider the smooth strict monotone quasi - uniform map @xmath139 , the so - called grid generating function , @xmath140 where @xmath141 $ ] , @xmath142 $ ] , and @xmath143 is a control parameter . we notice that for ( [ eq : qu1 ] ) @xmath144 . the problem under consideration can be discretized by introducing a uniform grid @xmath145 of @xmath146 nodes on @xmath147 $ ] with @xmath148 and @xmath149 with @xmath150 , so that @xmath151 defines a quasi - uniform grid on @xmath152 $ ] . the last interval in ( [ eq : qu1 ] ) , namely @xmath153 $ ] , is infinite but the point @xmath154 is finite , because the non integer nodes are defined by @xmath155 with @xmath156 and @xmath157 . the map allows us to describe the infinite domain by a finite number of intervals . the last node of such grid is placed on infinity so that right boundary condition is taken into account correctly .    for the sake of simplicity we consider here the simple scalar case . the finite difference formulae can be applied component - wise to a system of differential equations . we can define the values of @xmath34 on the middle - points of the grid @xmath158 as far as the first derivative is concerned we can apply the following approximation @xmath159 these formulae use the value @xmath160 , but not @xmath161 . both finite difference approximations ( [ eq : y ] ) and ( [ eq : dy ] ) have order of accuracy @xmath162 .    a finite difference scheme on a quasi - uniform mesh for the class of bvps ( [ p_qu ] ) can be defined by using the approximations given by ( [ eq : y ] ) and ( [ eq : dy ] ) . we denote by the @xmath163dimensional vector @xmath164 the numerical approximation to the solution @xmath165 of ( [ p_qu ] ) at the points of the mesh , that is for @xmath105 . we can define a second order finite difference scheme for ( [ p_qu ] ) as @xmath166   & \\[-1ex ] & { \bf g } \left ( { \bf u}_0,{\bf u}_j \right ) = { \bf 0 } \ ,   \nonumber\end{aligned}\ ] ] for @xmath167 , where @xmath168    it is evident that ( [ boxs_qu ] ) is a nonlinear system with respect to the unknown @xmath169dimensional vector @xmath170 . we notice that @xmath171 for all @xmath172 , but when @xmath173 , then @xmath174 and @xmath175 . on the contrary , we choose to set @xmath176 and @xmath177 in order to avoid a suddenly jump for the coefficients of ( [ boxs_qu ] ) . this produces a much smaller error in the numerical solution of the system at @xmath178 .    for the solution of ( [ boxs_qu ] ) we can apply the classical newton s method along with the simple termination criterion @xmath179 , @xmath114 and @xmath180 , is the difference between two successive iterate components and @xmath116 is a fixed tolerance . the results listed in the next section were computed by setting @xmath119 . figures  [ fig : ierley - ruehr_ns ] and [ fig : ierley - ruehr_s ] show the numerical solution of ocean model ( [ cmodel ] ) obtained setting the initial iterate @xmath181 from these figures we notice how the grid is denser close to the origin in comparison with the side of the far boundary at infinity .    [ ] [ ] @xmath30 [ ] [ ] @xmath81 [ ] [ ] @xmath82 [ ] [ ] @xmath182 [ ] [ ] @xmath79 ) with no slip boundary conditions ( [ bc1 ] ) obtained with the map ( [ eq : qu1 ] ) and @xmath183 for @xmath184 . we found a missing value of @xmath185.,title="fig : " ]    [ ] [ ] @xmath30 [ ] [ ] @xmath81 [ ] [ ] @xmath82 [ ] [ ] @xmath186 [ ] [ ] @xmath79 ) with slip boundary conditions ( [ bc2 ] ) obtained with the map ( [ eq : qu1 ] ) and @xmath183 for @xmath184 . we found a missing value of @xmath187.,title="fig : " ] in this paper we describe several methods for the numerical solution a simple wind driven circulation model arising in physical oceanography . our final aim is the comparison of numerical results . this is provided in tables  [ t1][t2 ] , where we used fbf ( free boundary formulation ) and qug ( quasi - uniform grid ) . for the sake of simplicity we limited ourself to compare the computed values of the missing initial condition @xmath49 . llrrr + numerical method & & grid - points & iter & @xmath188 + [ 1.5ex ] shooting - secant & @xmath189 & & 12 & 0.826111 + shooting - newton & @xmath189 & & 7 & 0.826111 + fbf & @xmath190 & @xmath191 & 11 & 0.826142 + fbf & @xmath192 & @xmath193 & 11 & 0.826140 + qug & @xmath194 & @xmath195 & 5 & 0.826180 + qug & @xmath194 & @xmath196 & 5 & 0.826150 +      we applied the simple shooting method to the truncated boundary formulation and a finite difference method to both the free boundary approach and the quasi - uniform grid treatment . for the sake of comparison , all numerical methods used in this study were second order methods . llrrr + [ -4ex ] numerical method & & grid - points & iter & @xmath197 + [ 1.5ex ] shooting - secant & @xmath189 & & 13 & 0.528885 + shooting - newton & @xmath189 & & 8 & 0.528910 + fbf & @xmath198 & @xmath191 & 11 & 0.528921 + fbf & @xmath199 & @xmath193 & 11 & 0.528921 + qum & @xmath194 & @xmath195 & 4 & 0.528927 + qum & @xmath194 & @xmath196 & 4 & 0.528922 +      the reported numerical results allow us to point out that the finite difference method with a quasi - uniform grid is the less demanding approach and that the free boundary approach provides a more reliable formulation than the classical truncated boundary . the shooting method supplemented by the newton s iterations shows that the ocean circulation model can not be considered as a simple test case . in fact , for this method we are forced to use as initial iterate a value close to the correct missing initial condition in order to be able to get a convergent numerical solution .    as a last remark , we want to mention that for the numerical solution of bvps on unbounded domains it is also possible to consider spectral methods that use mapped jacobi , laguerre and hermite functions ( see shen and wang @xcite for a review on this topic ) . w.  j. beyn . global bifurcation and their numerical computation . in d.  rossed , b.  d. dier , and a.  spence , editors , _ bifurcation : numerical techniques and applications _ , pages 169181 . kluwer , dordrecht , 1990 .
Abstract:

Answer: in this paper we report and compare the numerical results for an ocean circulation model obtained by the classical truncated boundary formulation , the free boundary approach and a quasi - uniform grid treatment of the problem . we apply a shooting method to the truncated boundary formulation and finite difference methods to both the free boundary approach and the quasi - uniform grid treatment . using the shooting method , supplemented by the newton s iterations , we show that the ocean circulation model can not be considered as a simple test case . in fact , for this method we are forced to use as initial iterate a value close to the correct missing initial condition in order to be able to get a convergent numerical solution . the reported numerical results allow us to point out how the finite difference method with a quasi - uniform grid is the less demanding approach and that the free boundary approach provides a more reliable formulation than the classical truncated boundary formulation . * key words . * nonlinear boundary value problems , infinite intervals , shooting methods , free boundary formulation , quasi - uniform grid , finite difference methods .    * ams subject classifications . * 65l10 , 65l12 , 34b40 .