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: graphene has drawn a lot of interest not only for their unique electronic properties @xcite but also for their unusual magnetic properties @xcite . magnetic ordering is typically observed in materials with partially filled d or f shell electrons . though there exists a large number of carbon based organic magnets with p electrons , the t@xmath0 is usually below 10 k @xcite . that is why observation of room temperature ferromagnetism in graphene related materials is so surprising . many theoretical works already exist to explain unusual magnetic property of graphene related materials . various reasons have been put forward for the formations of magnetic moments in graphene , such as , existence of zigzag edges @xcite , structural defects ( vacancies , adatoms , etc . , ) @xcite etc . , . zigzag edges induce highly degenerate localized electronic states precisely at the fermi level . these degenerate states are filled with electrons with the same spin in order to minimize the coulomb repulsion energy ( an effective hund s coupling ) . spin magnetic moments can also originate from the dangling bonds localized at the vacancy sites @xcite . j. zhou et . @xcite have predicted that moments can arise in semi - hydrogenated graphene sheets . l. xie et . @xcite observed room temperature ferromagnetism in partially hydrogenated epitaxial graphene . w. li et . @xcite theoretically showed that addition of monovalent and divalent adatom on graphene can also give rise to magnetic moments . this way of generating magnetic moment follows simply from lieb s theorem @xcite , which says that for half - filled bipartitite lattice systems the electrons in sublattice a and b are antiferromagnetically correlated . the ground state magnetic moment for such a system is , m= ( n@xmath1 -n@xmath2)@xmath3@xmath2 , where n@xmath1 and n@xmath2 are the number of a and b sublattice points . much less is known about the interaction and ordering if any between such local magnetic moments . vozmediano et . @xcite studied the interaction between localized moments mediated by rkky - like interaction but their calculated transition temperature was far below the room temperature . possibility of strong long range magnetic coupling via itinerant electrons in defective graphene ( vacancy or hydrogen chemisorption ) was predicted by yazyev and helm @xcite assuming metallicity of the system . l. pisani et . al . @xcite studied more realistic situation of a semiconducting ground state of defective graphene sheet and they suggested that long range spin magnetic ordering is possible in graphene . although there still persists a controversy that the observed magnetism is not its intrinsic property rather magnetism comes from ferromagnetic impurity present in the sample , experimentally ferromagnetism was observed in hopg only after proton irradiation @xcite and hence rules out the role of elemental impurity . y. wang et . @xcite observed room temperature ferromagnetism of partially reduced graphene oxide . but however , they did not stress on the evolution of magnetic hysteresis as a function of temperature , even though they observed an anomalous behaviour of increase in coercivity as a function of temperature . magnetic properties of graphene related materials depend upon their microstructure . structures of graphene oxide ( go ) and reduced graphene oxide ( rgo ) varies depending on the synthesis processes @xcite and consequently the magnetic property also changes @xcite .    here , we have studied the magnetic property of go and partially reduced go from room temperature ( 300k ) down to 10k . we find that both go and rgo show hysteresis at room temperature and as the temperature is lowered the hysteresis loop diminishes . in this paper we will address this anomalous behaviour of diminishing of hysteresis loop upon lowering the temperature for both go and rgo and try to understand the difference in the magnetic behaviour between go and rgo . graphene oxide ( go ) was synthesized using modified hummers method @xcite . to prepare go we used graphite powder , conc . h@xmath4so@xmath5 , nano@xmath6 and kmno@xmath5 . the reaction of h@xmath4so@xmath5 with strong oxidizing agents such as nano@xmath6 and kmno@xmath5 results in increase in the interlayer separation of the carbon layers due to incorporation / intercalation of oxygen atoms . this reaction leads to formation of graphite oxide flakes having single or few layer sheets . the solution was further ultrasonicated and centrifuged . the supernatant was collected and dried in vacuum drying oven at 80@xmath7c to obtain single layer graphene oxide samples . these go sheets were reduced by strong reducing agent hydrazine hydrate to form rgo . the resulting suspension is filtered and washed with water and ethanol and dried in an vacuum oven . the samples were characterised by uv - vis , ftir , raman spectroscopy , and sem . magnetic measurements were carried out using squid mpms xl ( quantum design ) . magnetization data as a function of temperature and magnetic hysteresis data at three different temperatures 10 k , 100 k and 300 k were taken . the magnetic moment values observed are comparable with that in earlier reports @xcite . 1 shows the ftir spectra of both go and rgo . inset of the fig . 1 shows the uv - vis spectroscopy data . we observe that the uv - vis spectra of go shows a sharp absorption peak near 238 nm which is shifted to 252 nm upon reduction . the shift of the peak position of rgo is due to restoration of sp@xmath8 carbon network after removing the functionalized groups as suggested in ref . ftir spectroscopy data ( fig.1 ) also indicates the reduction of oxygen - containing groups from graphene oxide sheet . the characteristic bands of go around 1051 @xmath9 ( c - o ) , 1223 @xmath9 ( c - o - c ) , 1613 @xmath9 ( c = o ) , 1731 @xmath9 ( o = c  oh ) and 3440 @xmath9 ( o - h ) are considerably decreased after the reduction which is consistent with the literature @xcite . in fig.2 we show the raman spectra of go and rgo showing two distinct peaks g and d. raman spectroscopy is a very important tool to study the average size of the sp@xmath8 region in sp@xmath8/sp@xmath10 combined material . graphene oxide shows the g peak near 1605 @xmath9 and d peak near 1370 @xmath9 which agrees with the previous reports @xcite . due to the reconstruction of the sp@xmath11 network at the oxidized region after reduction , the g peak in the rgo shifts to a lower wave number 1595 @xmath9 @xcite . in fig.3 we show the sem images of go and rgo . from the sem image we can clearly see that go shows folded layered structure and for the case of rgo , the sheets appear like unfurled structure . magnetization vs. temperature ( in log - log scale ) data of go and rgo are presented in fig . 4 . in this study a considerable magnetic moment has been observed in both samples . we observed less moment in rgo than go which we shall explain below . in the figure ( 5 @xmath12 6 ) we plot the magnetic hysteresis curve measured at three different temperatures ( 10 k , 100 k and 300 k ) for go and rgo respectively . m - t data on a log - log scale can be fitted with a straight line with a slope of -1 , from 10k upto about 60 - 70 k. almost no coercivity at 10 k supports paramagnetic behaviour at low temperature ( see inset of fig . 5 ) , but with increase in temperature , the hysteresis loop opens up i.e. , coercivity increases with temperature which is quite unusual . another interesting feature we observe is that with reduction of go , the magnetic moment decreases but coercivity increases . go has more paramagnetic contribution and does not saturate at high fields , while reduced graphene oxide ( rgo ) shows approach to saturation at high field along with larger coercivity indicating that upon reduction the magnetism in rgo becomes more cooperative in nature . several theoretical works have tried to explain ferromagnetism in graphene . m. wang et . shows by dft calculation that oxidized graphene can be ferromagnetic @xcite . ftir spectra indicate that graphene oxide has considerable amount of oxygen and oh groups in it . presence of the functionalized groups ( epoxide [ -o- ] , hydroxyl [ -oh ] , carbonyl [ -c = o ] and carboxyl [ -cooh ] ) on the graphene skeleton leads to disruption of sp@xmath8 bonded carbon network at many places . many topological defects are also formed , leading to local crumpling of planar graphene sheets . deviation from the planarity decreases the efficiency of @xmath13-overlap which reduces the connectivity of the @xmath13 network and causes an increase in @xmath14@xmath15 energy gap of the localized double bonds , hence reducing the ring current diamagnetism . therefore magnetic susceptibility / moment can increase upon oxidation of graphene @xcite . on the other hand , disruption of the hexagonal sp@xmath8 bonded carbon network presumably create many topological defects and isolated dangling bonds which carry local moments @xcite and also functionalized oxygenic groups ( -o - h ) can introduce localized magnetic moments @xcite . after on reduction of go , these functionalized groups are removed . this leads to two effects , ( 1 ) the @xmath13 electron with which the functional group was bonded is released as a carrier , and ( 2 ) the local crumpling of the planar carbon network caused due the presence of functional group is cured and locally the planar structure is restored . we can assume that there will be lesser oxygenic groups , defects , ruptures in the rgo compared to go . we find lesser paramagnetic moments in rgo compared to go as can be observed from fig . 4 of m versus t plot . this gives an indication that the magnetization is proportional to the amount of defects such as oxygen and hydroxyl impurities . thus , the magnetic moments in go and rgo is defect in origin . pure graphene is electronically half filled ( one @xmath13 electron per site ) on a bipartite lattice . the @xmath13 electrons on a and b sublattice are antiferromagnetically correlated . go and rgo are electronically insulating and so a description in terms of hubbard type model in real space rather than extended band electron picture is more appropriate . it has been shown that a strong impurity potentials like vacancy or any chemisorbed defect induces short range antiferromagnetic ( ferrimagnetic ) order around itself in a hubbard model on a half filled honeycomb lattice ( just like graphene ) @xcite . the net magnetic moment in a short range antiferromagnetic cluster is given by , m@xmath16@xmath17 -n@xmath18 as mentioned earlier . oxygen atoms bonds with two near neighbor c atoms of opposite sublattices and hence do not give rise to any local moment . oh on the other hand binds with only one c atom ( either a or b sublattice ) on graphene , this can create a local imbalance in n@xmath19 and n@xmath20 leading to a local moment @xcite . in a similar way atomic vacancy ( missing c atom ) can also create magnetic moment . we assume that go sheets have many patches of oh clusters , this patches are structurally distorted and electrically insulating locally because of gap at the local fermi level . nearest neighbour interaction between these @xmath13 electrons within such a patch is antiferromagnetic ( superexchange ) . the magnetic moment of such a small patch is @xmath21 . this uncompensated moment resides at the boundary of the patch . it is like an uncompensated nanosized antiferromagnetic particle with a net moment . most of these patches have very small net moments and behave paramagnetic . go has more paramagnetic moment compared to rgo . this seems natural as reduction process removes many such small moment carrying oh clusters from go . after reduction many oh groups are removed and hence we can consider less clustering of oh groups and isolated oh groups are scattered all over the sheet . since go has more oh groups than rgo , the net magnetic moment of go is higher than rgo . at low temperatures these isolated moments randomly placed in an insulating matrix with very little interaction between them due to highly crumpled structure behaves like paramagnets . it has been observed experimentally @xcite that the crumpling / wrinkling of graphene sheets decreases with increase in temperature . thus , with increase in temperature , crumpling is reduced and interaction between moments increases and hence deviation from paramagnetic behavior is observed . from fig.4 it can be seen that at low temperature due to increase in crumpling both go and rgo are paramagnetic , but at higher temperature(above 60 - 70 k ) we see for go , the moments drops below the paramagnetic slope ( see fig.4 ) and for rgo , the moments increases above the paramagnetic slope . hence go and rgo behaves like antiferromagnetic and ferromagnetic / ferrimagnetic respectively . now we discuss the m - h loop measurements . at 10 k the m - h curve is paramagnetic for both go and rgo showing no coercivity . at 100 k and 300 k both the samples shows clear hysteresis loops . we observe that the coercivity increases with increase in temperature which is very unusual . also the coercivity of rgo is higher than that of go at these temperatures . as we discussed earlier , both go and rgo can be considered to have large number of para moments coming from either defects ( vacancies ) or isolated oh groups . they might also have many clusters of oh groups which can be thought of as nanosized single domain particles . it is these blocked super moments which contributes to irreversibility and hysteresis . here the magnetization curve is known to arise from a competition between the domain anisotropy energy and the zeeman energy ( stoner wolfarth model @xcite ) . the coercive field is proportional to @xmath22 , where @xmath23 is the anisotropy energy density and @xmath24 is the net moment of any patch . if a patch is structurally very much crumpled then the average anisotropy energy of the patch will be small . magnetocrystalline anisotropy energy is very sensitive to the local perpendicular ( to graphene plane ) axis . on a patch having large number of o and oh groups , the variation of this axis from point to point is considerable . the coercive field is dependent on the spatial average ( over the whole patch ) of this anisotropy energy . average anisotropy energy density @xmath23 is smaller in go than in rgo , because go is structurally more disordered than rgo . rgo has much less number of o and oh groups and hence much less local crumpling that leads to higher anisotropy energy density k for the patch . so the coercive field in rgo is larger than go at any temperature as observed . since with increase in temperature the graphene sheets becomes on the average more planar , the average ` k ' of any patches of go and rgo increases , leading to increase in coercivity as temperature increases for both the systems . this we believe to be dominant reason for the increase of coercivity with increase in temperature . so far we have not discussed the interaction between the moments . it may be argued that with increase in temperatures , thermally activated carries could be responsible for interaction between the isolated moments , leading to long range magnetic ordering @xcite . with increasing temperature , more thermally activated electrons will come into play to increase coercivity along with gradual transition from paramagnetic to more nonlinear shape of m - h . since the gap arises due to distortion , defects etc . in the sp@xmath8 sheet , hence lesser gap at the fermi level of rgo than go could be the reason for increase of coercivity upon reduction . to summarize , we synthesized go and rgo chemically and the observed magnetic property of these two system are different . the paramagnetic part of the magnetic moment is larger in go compared to rgo . we argue that the origin of para moments are due to small patches of oh on graphene or other topological defects . reduction process removes many such oh containing regions as well as anneals some local topological defects . this might be the reason for lower paramagnetic moment in rgo compared to go . at higher temperatures there is finite hysteresis loops and the magnetic coercivity anomalously increases with increase in temperatures . we have argued that over and above the small defects and small patches of oh groups which gives paramagnetism , there are also many nanosized oh patches having considerable uncompensated moments . these are blocked along random directions and are responsible for irreversibility , hysteresis , saturation etc . the coercivity of such blocked supermoments are determined by its average local magnetic anisotropy energy . with increase in temperatures the wrinkles on go and rgo decreases on the average , this leads to an increase in the local magnetic anisotropy energy . hence , an increase in coercivity with increase in temperature is observed . larger coercivity in rgo compared to go at any temperatures is due to less number of functional groups than in go and hence less wrinkles on it . our arguments is also supported by the experimental observation of increase in wrinkle formation in graphene with the decrease in temperature . we would like to thank prof . a. roy of iit kharagpur for raman measurement of our samples . we would like to thank amity university , noida for encouraging jk and arc to work in sinp . h. c. neto , f. guinea , n. m. r. peres , k. s. novoselov and a. k. geim ; rev . phys . * 81 * , 109 ( 2009 ) . a.k.geim and k. s. novoselov ; nature mater . * 6 * , 183 ( 2007 ) . a. k. geim and a. h. macdonald ; phys . today * 608 * , 35(2007 )            m. tamura , y. nakazawa , d. shiomi , k. nozawa , y. hosokoshi , m. ishikawa , m.takahashi and m.kinoshita ; chem . * 186 * , 401404(1991 ) . p. m.alemand , k. c.khemani , a.koch , f.wudl , k.holczer , s.donovan , g.gruner and j. d.thompson ; science , * 253 * , 301303(1991 ) . a. rajca , j. wongsriratanakul and s. rajca ; science * 294 * 16(2001 )    h. kumazaki and d. s. hirashima ; j. phys . 77 , no . 4 044705 ( 2008 ) , s. bhowmick and v. b. shenoy ; j. chem . phys.128 , 244717 ( 2008 ) , j. jung and a. h. macdonald ; phys . rev . b 79 , 235433 ( 2009 ) . v. v. nelayev and a. i. mironchik ; materials physics and mechanics 9 26 - 34 ( 2010 ) r. singh and p. kroll ; j. phys . : condens . matter 21 196002 ( 2009 ) . j. j. palacios and j. fernndez - rossier and l. brey ; phys . rev . b 77 , 195428 ( 2008 )    s. choi , b. w. jeong , s. kim and g. kim ; j. phys . : condens . matter * 20 * , 235220 ( 4pp)(2008 ) . v. m. pereira , f. guinea , j. m. b lopes dos santos , m. n. r. peres and a. h. castro neto ; phys . lett . * 96 * 036801 ( 2006 ) . y. h. kim , j. choi , k.j . chang and d. a. tom ; phys . rev.b * 68 * 125420(2003 ) . l. xie , x. wang , j. lu , z. ni , z. luo , h. mao , r. wang , y. wang , h. huang , d. qi , r. liu , t. yu , z. shen , t. wu , h. peng , b. zyilmaz , k. loh , a. t. s. wee , a. and w. chen ; appl . . letts . * 98 * , 193113 ( 2011 )                      c. hontoria - lucas , a. j. lopez - peinado , j. d. d. lopez - gonzalez , m. l. rojas - cervantes and r. m. martin - aranda ; carbon * 33 * , 15851592(1995 ) . g. i. titelman , v. gelman , s. bron , r. l. khalfin , ; y.cohen , bianco - peled , h.carbon * 43 * , 641649(2005 ) . n. b. colthup , l. h. daly and s. e. wiberley ; introduction to infrared and raman spectroscopy , 3rd ed . ; academic press : london , 1990 . s.stankovich , r.d.piner , s.t . nguyen and r. s. ruoff , carbon ( 2006 ) , 44 , 33423347 . m. baraket , s.g . walton , z. wei , e.h . lock , j.t . robinson , p. sheehan , carbon * 48 * , 3382 ( 2010 ) . s. stankovich , d.a . dikin , r.d . piner , k.a . kohlhaas , a. klein - hammes , y. jia , y. wu , s.t . nguyen and r.s . ruoff , carbon * 45 * , 1558 ( 2007 ) . f.tuinstra and j. l. koenig , j. chem . * 53 * , 11261130(1970 ) . * figure 5 * : ( color online ) magnetic hysteresis of go at temperature 10k , 100k and 300k . inset ( top left ) : magnetic hysteresis of go and rgo at 10k . inset ( bottom right ) zoom region of the hysteresis curve . * figure 6 * : ( color online ) magnetic hysteresis of rgo at temperature 10k , 100k and 300k . inset ( top left ) : magnetic hysteresis of go and rgo at 300k . inset ( bottom right ) zoom region of the hysteresis curve .
Abstract:

Answer: in this report we present the temperature evolution of magnetic coercivity of graphene oxide ( go ) and reduced graphene oxide ( rgo ) . we report an anamolous decrease in coercivity of go and rgo with decreasing temperature . we could explain this anamolous behavior invoking the inherent presence of ripple in graphene . we observe antiferromagnetic and ferromagnetic behavior at room temperature for go and rgo respectively , but at low temperatures both shows paramagnetic behavior .


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 existence and uniqueness of solutions to the differential equations of physics is seldom an issue . the very fact that these equations describe physical reality seems to argue that their solutions must exist and be unique . textbook examples of nonuniqueness , for example clairaut s equation ( ref . @xcite , p. 94 ) , seem like exceptional cases , which would not arise in physics in any case .    there is , however , a system of differential equations , arising as a limiting case of a much studied physical problem , in which uniqueness fails for _ every _ solution at _ every _ time . the system is singular laplacian growth @xcite , @xcite . one might call this behavior `` maximally nonunique . '' it is so different from the usual behavior of differential equations that it hardly seems like a differential system at all . the way it occurs is the following : the theory of singular laplacian growth describes the motions of certain singularities of a conformal map , which are all located on the unit circle in the complex plane . they move on the circle , but they can also `` split , '' introducing new singularities , at any time . this behavior , which does not sound very remarkable , essentially implies the maximal nonuniqueness property . in ref . @xcite the term `` fragile '' was suggested for such a system , because its distinguishing feature is that its singularities can break apart . it is natural to ask if other physical systems might reduce to a fragile system in some limit , and if the fragile property manifests itself in the behavior of the system . turbulence is certainly a candidate to be a fragile system : like laplacian growth , turbulence is a phenomenon in which a differential equation has unexpectedly complex solutions . in this paper i argue that turbulence is a system that has a fragile ( maximally nonunique ) limiting case . the limiting case is the burgers equation . ( on rather different grounds a resemblance between the problems of turbulence and laplacian growth has been noted by hastings and levitov @xcite ) . the burgers equation @xcite , @xcite @xmath0 has played a shifting role in the problem of turbulence . it was first introduced as a one - dimensional zero pressure version of the navier - stokes equation ( with @xmath1 the viscosity ) in the hope that it might exhibit in its solutions some of the complexity of turbulence . by the hopf - cole trick , though , it was found to be equivalent to a linear diffusion equation , so that the initial value problem could be solved explicitly : one can write a formula for the solution ! that is surely too simple to be turbulent . attention was thus focussed on the only case which still seemed promising , the @xmath2 limit , the inviscid burgers equation , which is a conservation law for @xmath3 . in this limit the solution may develop discontinuities ( `` shocks '' ) . the distribution of shocks , and their development in time , for random initial data , is a problem called burgers turbulence . the reason that one should think of this as a @xmath2 limit , and not simply as a @xmath4 version of the equation , is that for @xmath4 the solution is ambiguous . to be sure , the solution of the inviscid equation is immediate by the method of characteristics : @xmath3 is constant along lines @xmath5 ( hence straight lines in the x - t plane ) . but if the characteristics cross , which they certainly will do in the case of random initial data @xmath6 , then one has multiple determinations of @xmath7 , as in fig . 1 . the resolution of this difficulty is to take the formula for @xmath7 in case @xmath8 , which is unambiguous , and let @xmath1 approach zero . the result is as shown in fig . 2 : a discontinuity forms at a definite location , separating one determination of @xmath3 from another . the same procedure shows that shocks interact in a definite manner , coalescing as shown in fig . 3 when one overtakes or collides with another . the rules for how these discontinuities form and interact were known long ago in the theory of compressible gas dynamics , where they are called the rankine - hugoniot jump conditions ( see ref . @xcite , pp . 488 - 90 , ref . @xcite , p. 596 ) . they are justified by appeal to the second law of thermodynamics : while @xmath3 is conserved , entropy must increase in the shock . that is what the argument with @xmath8 also does : we resolve the ambiguity in the inviscid burgers equation by putting in a dissipative term with the right sign . the second law of thermodynamics says mechanical energy should be dissipated in the shock , and not created . it is this condition which leads to the interaction rule illustrated in fig . now recall the suggestion of onsager , as emphasized by alexandre chorin , among others , that a statistical theory of the turbulent steady state should be characterized by a negative temperature ( see ref . @xcite , chapter 4 ) . what is meant is that there is a maxwell - boltzmann probability distribution @xmath9 which , instead of giving more weight to the low energy microstates in the ensemble description of the macrostate , as is usual , gives more weight to high energy microstates . ( this statistical `` temperature '' of turbulence has nothing to do with usual thermodynamic temperature , which is only weakly coupled to the mechanical degrees of freedom of interest . ) let us accept this idea for the moment . the second law says that entropy should increase , or perhaps better , in this more general situation , that information should be lost , in irreversible processes . stated in terms of the free energy @xmath10 ( and remembering @xmath11 ) this says that free energy should _ increase _ in the approach to the steady state . this is opposite to what is usual . putting onsager s idea together with the ideas of burgers turbulence thus requires that mechanical energy should be created in the shock rather than dissipated . we see that the limit in the inviscid burgers equation should be @xmath12 , i.e. , we should imagine the viscosity approaching zero through _ negative _ values , from _ below_.    the inviscid burger s equation does not contain temperature or @xmath1 explicitly , of course , but the rules for how the shocks form and interact are now different . the following argument gives the simple idea which is at the heart of this paper . we do not try to give it an elaborate justification , since , like many simple ideas , it may contain some truth even if the arguments are wrong . the argument is that there is a symmetry of the ( general ) burgers equation , @xmath13 we can convert the @xmath14 inviscid burgers equation into the familiar @xmath8 inviscid burgers equation by reversing the signs of @xmath15 and @xmath16 , i.e. , reflecting graphs like fig . 3 through the origin . the result is fig . 4 , in which a single shock has spontaneously split into two , with no reference to initial conditions . this is the fragile property . if shocks can split at any time , as indicated , then the solutions to the differential equation are maximally nonunique . our aim in the previous section was to show that the equations which describe fluid flow , the navier - stokes equations , become maximally nonunique , or `` fragile , '' in some limit . this limit is rather far removed from physical reality , however . one naturally wonders how , if at all , the nonuniqueness property might manifest itself in a real system . the example of laplacian growth encourages one to think that real processes would not completely obscure the underlying `` fragile '' processes @xcite . the system we have imagined is characterized by two temperatures : a `` superhot '' negative temperature , which describes the ensemble of microscopic entities which make up the turbulent state , and the usual temperature , which describes the ensemble at the molecular level . these two ensembles interact only weakly , it has been argued , but , again by the second law of thermodynamics , to the extent that they interact , there must be a flow of energy from the first to the second , from the turbulent ensemble to the molecular ensemble . this picture is reminiscent of the kolmogorov cascade idea ( ref . @xcite , chapter 3 ) , and suggests identifying the turbulent ensemble with the `` energy range , '' ( a hypothetical range in @xmath17 space containing most of the energy ) , and the molecular ensemble with the `` dissipation range , '' ( a range in @xmath17 space , disjoint from the energy range , in which dissipation is important ) , i.e. , assigning a temperature gradient to the kolmogorov picture , with negative temperature at small @xmath17 and positive temperature at large @xmath17 . it is in the energy range , then , at small @xmath17 , that the fragile processes of nonuniqueness would occur . these processes extract mechanical energy from the negative temperature `` bath , '' and this energy must ultimately derive from the forces maintaining the turbulence . thus the onset of fragile processes would appear , on the macroscopic scale , as an increase in resistance to these applied forces . at a more microscopic level it would be the onset of nonuniqueness , allowing splitting and proliferation of microscopic entities as in fig . 4 . at a still more microscopic level , corresponding to the dissipation range , the usual picture would apply , and one would have dissipative processes like fig . 3 . in this view turbulence _ is _ the visible manifestation of nonuniqueness on the intermediate scale of the energy range , and the work done by external forces goes directly into the proliferation of singularities at that scale , and only indirectly into dissipation .    in terms of modelling turbulence , it suggests that the proliferative processes of the energy range , which are still hypothetical , may be like the dissipative processes of the dissipation range , but time reversed ( and on a larger length scale ) . in this way the abstract picture suggested here might persist in a more realistic dynamics . the issue of nonuniqueness may also be relevant to cfd modelling of turbulence . the algorithms of differential equation solvers are not set up for equations which do not have unique solutions . experience with laplacian growth confirms that proximity to a nonunique model may indicate trouble for conventional numerical solutions .    99 r. courant and d. hilbert , _ methods of mathematical physics _ , vol 2 ( interscience publishers , 1962 ) . peterson , * phys 62 * , 284 ( 1989 ) . peterson , cond - mat/9710046 . hastings and l.s . levitov , cond - mat/9607021 . burgers , _ the nonlinear diffusion equation _ , ( riedel , dordrecht , 1974 ) . g. strang , _ introduction to applied mathematics _ , ( wellesley - cambridge press , wellesley ma , 1986 ) . chorin , _ vorticity and turbulence _ , ( springer , 1994 ) .
Abstract:

Answer: the possibility is considered that turbulence is described by differential equations for which uniqueness fails maximally , at least in some limit . the inviscid burgers equation , in the context of onsager s suggestion that turbulence should be described by a negative absolute temperature , is such a limit . in this picture , the onset of turbulence coincides with the proliferation of singularities which characterizes the failure of uniqueness . epsf    2


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 dynamics of fluid particles in an immobile environment has been of fundamental interest to understand anomalous transport processes in confined geometry and porous media  @xcite . examples include different types of lorentz gas models @xcite , polymers in quenched disorder @xcite , hard  sphere mixtures @xcite and ion  conducting silicates @xcite . simple model systems that provide a time  scale separation of transport properties among different species are mixtures of small and large particles at high densities . while it is difficult to experimentally realize such systems on an atomistic scale ( for an exception see ref . @xcite ) , it is possible to carry out experiments on colloidal suspensions that contain disparately  sized particles . about a decade ago , imhof and dhont @xcite performed dynamic light scattering experiments on a binary mixture of colloidal silica particles with a size ratio of 1:9.3 . the effective interactions between colloids are hard  sphere  like , hence phase behavior and transport properties are governed by packing effects . an interesting finding of imhof and dhont is the existence of different phases where the large particles exhibit a structural arrest , yet the small particles are still mobile . the simplest case of such a phase could be one consisting of mobile small fluid - like particles in a crystalline matrix of large particles . it can be realized , at least , in an intermediate regime below the freezing transition .    in this work , we use molecular dynamics simulations to study a binary mixture of charged particles with a size  ratio of 1:5 . similar to the experiments by imhof and dhont @xcite , the system exhibits a phase transition from a fluid to a mixture of crystalline large particles and fluid  like small particles . when temperature is sufficiently lowered , the small particles are localized and crystallize , eventually , in a sublattice with respect to the large  particle lattice . in the temperature regime between the freezing transition and the final crystallization into a sublattice , the movement of the small particles can be characterized by a hopping motion in an external periodic potential ( created by the large particles ) . similar to supercooled liquids , a `` @xmath0 relaxation regime '' can be identified in the incoherent intermediate scattering functions @xmath1 ( with @xmath2 the wavenumber and @xmath3 the time ) of the small particles . these functions show a plateau region between the microscopic regime and the long  time decay to zero . at a given wavenumber , the height of the plateau increases with decreasing temperature . this indicates that the particles are more and more localized . thus , the small particles display a kind of delocalization  to  localization transition around the liquid  to  crystal transition of the large particles . although the effective packing fraction of the small particles , considered in this work , is relatively small , structural correlations among the small particles occur , associated with a complex distribution of single particle potential energies . the system considered in this work shows various nontrivial features that might be generic for the dynamics in confined geometry and porous media , especially when coulomb  like interactions become important . it can be also realized experimentally in charged colloidal suspensions . therefore , with a similar system as the one used by imhof and dhont @xcite the predictions of our model could be rationalized experimentally . on the other hand , our model considers essentially the motion of mobile particles in a periodic potential . it should be much simpler to develop analytic theories for this case than for the transport of particles in a disordered potential energy landscape . charged particles ( colloids ) are modelled by an effective screened coulomb ( or yukawa ) potential , @xmath4 where @xmath5 is the distance between particles @xmath6 and @xmath7 . the parameters @xmath8 ( @xmath9 ) denote the distance between two particles at contact , @xmath10 , with @xmath11 the radius of an @xmath12-particle . the interaction range is characterized by the screening parameter @xmath13 . in the limit @xmath14 , the yukawa potential corresponds to bare coulomb interactions . a shifted lennard  jones ( wca ) potential is used to model the excluded volume of the particles : @xmath15 + \varepsilon , & r\leq2^{\frac16 } \sigma_{\rm { \alpha}{\beta } } , \label{shiftedlj}\\ 0 , & r > 2^{\frac16}\sigma_{\rm { \alpha}{\beta}}. \end{array}\right.\end{aligned}\ ] ] the energy parameter @xmath16 is set to @xmath17 . note that the potential , eq . ( [ shiftedlj ] ) , is purely repulsive . it can be regarded as an approximation to a hard  sphere interaction for two particles at contact . the binary yukawa system of volume fraction @xmath18 consists of @xmath19 small and @xmath20 big particles in three dimensional space . the size ratio between them is @xmath21 . we choose @xmath22 , @xmath23 and @xmath24 , and mass @xmath25 . the yukawa interactions are short  ranged with a fixed value of the inverse screening parameter ( screening length ) @xmath26 . note that a cut  off is introduced at @xmath27 where potential energies @xmath28 are of the order of @xmath29 . for the energy parameters , the values @xmath30 , @xmath31 , and @xmath32 are chosen . this choice implies that small and large particles have the same surface charge density . in the following , all the physical quantities are measured in units of the large particle s mass @xmath33 , diameter @xmath34 and energy @xmath35 . the boltzmann constant is set to @xmath36 . newton s equations of motion for the charged particles are integrated using the velocity form of the verlet algorithm . the time step is chosen to be @xmath37 for @xmath38 , @xmath39 for @xmath40 , and @xmath41 for @xmath42 ( note that @xmath43 is measured in units of the time @xmath44 ) . at each temperature , the system was first equilibrated in the @xmath45 ensemble by coupling it to a stochastic heat bath . the equilibration time was sufficiently longer than the time needed to observe a diffusive motion of the small particles ( typically , a factor of @xmath46 longer for @xmath47 , and a factor of @xmath48 for @xmath49 . ) at the lowest temperatures ( @xmath50 ) , the production runs were over 10 to 20 million time steps . this was not enough to see a diffusive motion of the small particles . note that for the equilibration at a given temperature , we used a final configuration of the next  higher temperature as an initial configuration . microcanonical production runs were performed at the temperatures @xmath51 , 2.5 , 2.0 , 1.0 , 0.8 , 0.5 , 0.4 , 0.2 , 0.15 , 0.1 , 0.05 , 0.04 , 0.03 , 0.025 , 0.02 , 0.015 , 0.012 , 0.01 . at each temperature , eight independent runs were done to improve statistics . first , we consider partial pair correlation functions @xcite @xmath52 ( @xmath53 ) that measure structural correlations between small ( `` s '' ) and large ( `` l '' ) particles . these quantities are proportional to the probability to find a particle of type @xmath12 at a distance @xmath5 from a particle of type @xmath0 .    upon decreasing temperatures the large particles crystallize around @xmath54 into a fcc lattice structure . [ fig1 ] shows the @xmath55 above and slightly below the freezing transition of the large particles at @xmath51 and @xmath56 , respectively . we see that the curve for @xmath57 at @xmath51 shows only tiny differences to that at @xmath56 . this also holds for the functions @xmath58 at the two temperatures . beyond the first peak the amplitude in @xmath59 and @xmath58 is only slightly different at @xmath56 and at @xmath51 this indicates a fluid  like behavior in the small  small and small  large correlations . in contrast , the function @xmath60 changes significantly from @xmath56 to @xmath51 . note that the rather broad peaks at @xmath56 are reminiscent of a pair correlation function of a normal liquid state . however , a closer inspection of other structural quantities such as the static structure factor shows that the large particles form a fcc crystal and the broad peaks in @xmath60 are due to thermal fluctuations at the relatively high temperature @xmath56 . this can be also infered from the snapshot at @xmath56 ( see fig . [ fig1 ] ) . one can clearly identify the crystalline planes , which are strongly distorted by the thermal motion of the large particles .    by decreasing temperature from @xmath56 to @xmath61 , significant structural changes occur . in fig . [ fig2 ] , both @xmath60 and the snapshot clearly show the long  range structural order of the fcc crystal at @xmath61 . the crystalline structure of the large particles is also reflected in the correlations among the small particles at @xmath61 . the function @xmath57 displays peaks at @xmath62 and @xmath63 . the second peak at @xmath64 is at the same position as the first peak in @xmath60 , which corresponds to the nearest neighbor distance in the fcc crystal . this indicates the tendency of the small particles to populate a sublattice with respect to the fcc structure of the large particles . an interesting feature is the rather sharp peak of @xmath57 at @xmath65 . this is due to the coordination of a small particle with another small particle . it means that some of the small particles form pairs in a given local site . note that approximately 75% of the small particles have no nearest neighbor , 20% exhibit a pairing , and 5% of the small particles are two  fold or three  fold coordinated . the occurrence of the small particles with zero neighbour or pairs in local sites is also reflected in @xmath58 . the first peak in this function shows a shoulder around @xmath66 , which is approximately the distance between a large particle and the next  nearest small particle of a small particle pair . the periodic structure of the large particles crystalline matrix is accompanied by fluid  like small particles , at least at high temperatures . to characterize in detail the `` disorder '' in the small  small correlations , we show in fig . [ fig3 ] the small particles distribution of single  particle potential energies ( denoted by @xmath67 ) at four different temperatures . it is remarkable that the main peak in @xmath68 moves to higher energies while temperature decreases . this means that the formation of the fcc structure by the large particles is not energetically favoured by the small particles as it increases their potential energy on average . another significant feature of @xmath67 is the emergence of a high  energy tail . this is magnified in the inset of fig . [ fig3 ] by a double  logarithmic plot of the data . this high  energy tail tends to dissappear towards low temperatures . this is due to the fact that the kinetic energy sets an upper bound for the accessible potential energies . having discussed static properties of the model , we now turn our attention to transport properties of the small particles in the immobile environment of the large particles . fig . [ fig4]a displays the small  particle self - diffusion constant @xmath69 as a function of temperature in a double  logarithmic plot ( note that we determined @xmath69 from the mean  squared displacement of the small particles via the einstein relation @xcite ) . as temperature is decreased from @xmath51 ( i.e.  above the freezing transition of the large particles ) to @xmath61 , the self - diffusion constant decreases by three and a half orders of magnitude . note that the functional form of @xmath69 can neither be described by a power law nor by an arrhenius behavior . a detailed description of the small  particle dynamics can be obtained from the incoherent intermediate scattering functions @xmath1 , the self  part of the time  dependent density  density correlation function for the small particles @xcite . in fig . [ fig4]b this function is shown for different temperatures at @xmath70 . this wavenumber corresponds to the location of the first bragg peak in the static structure factor for the large  large correlations . note that for other wavenumbers a similar behavior of @xmath71 is observed . above the liquid  to  crystal transition of the large particles at @xmath54 , @xmath72 exhibits an exponential one  step decay to zero , as expected for simple liquids . below @xmath73 , a plateau develops at intermediate times , followed by a _ non  exponential _ decay to zero . the emergence of the plateau is associated with a caging of the particles . in this case , the small particles are trapped in local potential basins created by the large particles . the height of the plateaus in @xmath71 ( also called lamb  mssbauer factor @xcite ) measures the localization of the particles with respect to the length scale corresponding to the considered wavenumber @xmath2 @xcite . we can infer from fig . [ fig4]b that around @xmath73 the small particles undergo a transition from a completely delocalized state ( with vanishing plateau ) to a more and more localized state towards low temperatures . in the light of the potential energy distribution ( see fig . [ fig3 ] ) , the delocalization  to  localization transition of the small particles is related to the gradual disappearance of the high  energy tails .    the plateau sets in by an oscillation . with decreasing temperature , this oscillation and the whole microscopic regime , i.e.  the decay of the curves onto the plateau , shift to the right on the time axis this behavior corresponds to a softening of the cage towards lower temperatures and is the subject of a forthcoming publication . further information on the long  time regime and the diffusive motion of the small particles can be extracted from the fourier transform of @xmath71 , the self - part of the van hove correlation function , @xmath74 @xcite . @xmath75 is the probability to find a particle at time @xmath3 at a distance @xmath5 away from the origin at @xmath76 . in fig . [ fig5 ] , this quantity is shown for different times at @xmath56 and @xmath61 in a linear  logarithmic plot . at the high temperature , we observe a regular behavior , which is similar to that of simple liquids . with increasing time , the location of the peaks moves continuously to larger distances . however , the behavior is totally different at low temperatures : at @xmath77 , the correlation function exhibits a single peak . as time goes on , a second peak develops and if one waits until @xmath78 five peaks can be identified . the distance between the peak maxima is around @xmath79 , which corresponds to the lattice constant of the large particle s fcc lattice . we can conclude from this that at low temperatures the small particles do not diffuse in a continuous manner but discontinuously in time by hopping along the sublattice sites of the large particle s fcc lattice . note that this interpretation has been first given for similar features in the glassy behavior of a soft  sphere model by roux _ et al . _ we have presented an extensive md simulation study of a binary yukawa mixture with size ratio 1:5 . in this system , interesting phase behavior is observed where the crystallization of the large particles is accompanied by a dynamic delocalization  to  localization transition of the small particles . we speculate that this kind of phase behavior is generic for broad classes of complex systems , in particular mixtures of disparate  sized ( charged ) colloids ( see , e.g. , the experiments by imhof and dhont @xcite ) . the motion of proteins in lipid cubic phases @xcite might be also very similar to that of the small particles in our yukawa mixture . however , we are not aware of any systematic experiments on fluid  like particles in a crystalline matrix . we hope the present work will stimulate new experimental efforts in this direction .    financial support of the dfg ( sfb 625 , sfb tr6 , and the emmy noether programme , grant no . ho 2231/2 ) and the mwfz mainz are gratefully acknowledged . we thank the nic jlich for a generous grant of computing time on the jump . 0 recent developments can be found in , _ proceedings of the international workshop on dynamics in confinement _ , , _ proceedings of the second international workshop on dynamics in confinement _ , . . . . . . . . . . . . . . . . . .
Abstract:

Answer: molecular dynamics computer simulations are used to investigate the dynamics of a binary mixture of charged ( yukawa ) particles with a size  ratio of 1:5 . we find that the system undergoes a phase transition where the large particles crystallize while the small particles remain in a fluid  like ( delocalized ) phase . upon decreasing temperature below the transition , the small particles become increasingly localized on intermediate time scales . this is reflected in the incoherent intermediate scattering functions by the appearance of a plateau with a growing height . at long times , the small particles show a diffusive hopping motion . we find that these transport properties are related to structural correlations and the single  particle potential energy distribution of the small particles .


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: whereas there are a number of examples of exact analytic solutions of time - independent problems in quantum mechanics , such solutions of time - dependent problems are few . the analytic solutions of both types that do exist tend to provide approximate models to actual physical systems . no doubt such solutions are instructive for gaining insight into the behavior of the physical systems that they describe . nevertheless because the models themselves are often approximations and because one wishes to describe real systems as precisely as possible , one also relies on accurate numerical methods .    the time - dependence of nonrelativistic quantum systems , which is the focus of this paper , has become important in diverse areas of atomic and subatomic physics . examples of these include the study of nuclear processes such as the decay of unstable nuclei and associated phenomena like atomic ionization  @xcite and bremsstrahlung  @xcite , the study of fundamental processes necessary for quantum computing  @xcite , the study of mesoscopic physics or nanophysics devices  @xcite , and the motion of atoms in a trap . a reliable and accurate numerical determination of the time - dependent wave function such as we discuss in this paper will no doubt be necessary and/or helpful in making advances in the understanding of a variety of quantum processes .    in this paper we consider the numerical solution of the time - dependent nonrelativistic schrdinger equation . much has been learned about basic scattering processes from the numerically generated solutions of traveling wavepackets as they pass through a potential region  @xcite , as well as the time - evolution of unstable quantum processes  @xcite . however , the methods used in the past , and still employed currently , are limited in that the solutions often degrade after a certain time interval , so that they reduce to noise . furthermore for processes in which the wave function spreads or travels away from the source one often requires such a large number of space steps that the computation becomes prohibitive . the goal of this study is to improve the existing standard approach by allowing for relatively large step sizes both in time and in space and thus to reduce the number of basic arithmetical calculations while obtaining more accurate solutions . we have been able to make significant improvement to one conventional approach , _ viz . _ , the crank - nicolson ( cn ) implicit integration scheme for the time - dependent schrdinger equation . many years ago the cn approach was shown to be successful in the study of wavepacket scattering in one dimension by goldberg _ et al . _ @xcite . in recent years the cn method continues to be employed for its space- and/or time - development algorithm to study various time - dependent problems . see , for example , refs . the attractive aspect of this method is that the solution is constrained to be unitary at every time step . it is this constraint that makes the solution stable regardless of time- or space - step size . although the evolution of the solution is unitary , the wave function is not correct if the step sizes are too large . the error is of @xmath5 , where @xmath6 and @xmath7 are the spatial and temporal step size , respectively . the method was successfully generalized to two dimensional scattering by galbraith _ et al . _ @xcite and , more recently , to multichannel scattering  @xcite . furthermore alternative methods which are fast computationally were introduced by kosloff and kosloff  @xcite . these involve the fast fourier transform of the kinetic energy operator of the schrdinger equation . variants of this method were discussed in ref .  @xcite . although this approach is fast and is able to handle large time intervals in one pass , it is not unitary and does require a large number of space intervals .    improved cn algorithms have been discussed by a few authors . et al . _ @xcite introduced a seven - point formula for the second - order spatial derivative with error of @xmath8 and an improved time - integration scheme with an error of @xmath9 . they claim to obtain two orders of magnitude improvement in the time advance . moyer  @xcite uses a numerov scheme for the spatial integration method with error of @xmath10 but has a one - stage time evolution giving the cn precision in time , i.e. , with an error of @xmath11 . moyer also introduces transparent boundary conditions for unconfined systems . we have found these very useful for making long - time or large - space problems tractable  @xcite , but we will not discuss such boundary conditions further in this paper . et al . _ @xcite indicate how to generalize the time development to higher order , but do not discuss spatial integration . in this paper we present a generalization of the often - used cn method of obtaining numerical solutions of the time - dependent schrdinger equation . the generalization yields numerical solutions accurate to order @xmath0 in space and @xmath1 in time for any positive integers @xmath2 and @xmath3 , while cn employ @xmath4 . by appropriate choice of @xmath2 and @xmath3 the improvement can be of such a nature that hitherto computationally unfeasible problems become doable , and solutions with low to modest precision can now be obtained extremely accurately .    in the following we consider the generalization of spatial integration in sec . [ spatial ] and the generalization of the time integration in sec . [ time ] . in sec . [ errors ] we discuss errors and a way of dealing with a particular type of boundary condition . we study specific examples to illustrate the improvement of the generalizations over the standard cn procedure in sec . [ examples ] . some general observations and conclusions are made in sec . [ conclusions ] . we describe a general procedure for solving the one - dimensional time - dependent schrdinger equation @xmath12 where @xmath13 and @xmath14 is a given wave function at initial time @xmath15 . in this section we use the standard time - advance procedure of the cn method , but generalize the spatial integration . in sec .  [ time ] we generalize the time - evolution procedure . the time evolution of the system can be expressed in terms of an operator acting on the wave function at time @xmath16 which gives the wave function at a later time @xmath17 according to the equation @xmath18 + the time - evolution operator @xmath19 can be expanded to give a unitary approximation of the operator by setting @xmath20 inserting the approximate form of the operator into eq . ( [ eq:2.3 ] ) , we obtain the equation @xmath21 with an error of @xmath11 . here we focus on the second - order spatial derivative in @xmath22 of eq . ( [ eq:2.2 ] ) and leave improvements with respect to the time derivative to sec . we generalize the usual three - point formula and the seven - point formula of miicu _ et al . _ @xcite , for the second - order derivative to a @xmath23-point formula . such a formula has the form @xmath24 where @xmath25 are real constants . to obtain the coefficients @xmath25 we make expansions    @xmath26    for @xmath27 ; @xmath28 denotes the @xmath29th derivative with respect to @xmath30 . when we add the two equations , the terms with odd - order derivatives cancel , resulting in the equation @xmath31 thus we obtain the system of @xmath2 equations in @xmath2 unknowns , i.e. , @xmath32 for @xmath33 , @xmath34    we solve these equations to obtain @xmath35 . it is evident from the terms on the right side of eqs . ( [ eq:2.8 ] ) that @xmath35 has the form of eq . ( [ eq:2.6 ] ) and the coefficients @xmath25 can be identified . because the equations  ( [ eq:2.8 ] ) are invariant under the change of @xmath36 to @xmath37 , the coefficients satisfy the relation @xmath38 for @xmath39 . for example , the first seven sets of coefficients ( up to the fifteen - point formula ) are given in table  [ table1 ] . coefficients @xmath25 up to @xmath40 . [ cols="^,^,^,^,^,^,^,^,^",options="header " , ]     we list some of the computed results . we observe that the traditional cn method ( @xmath41 and @xmath42 ) has a low precision and that using greater @xmath43 ( smaller @xmath6 ) results in modest gain in precision . by using higher - order time formula one can make significant gain in precision ( seven orders of magnitude ) with no increase in computational time compared with that of ref . ( compare the first row to the last two rows of table  [ table4 ] . ) rows 5 through 9 of table  [ table4 ] illustrate the transition from less precise to more precise solutions . it is consistent with the finding of the authors of ref . @xcite who use an @xmath44 , @xmath45 method and obtain two orders of magnitude improvement of the results of refs . the results are sensitive to surprisingly high orders of @xmath6 and @xmath7 . another test using wavepacket scattering to show the efficacy of the higher - order approach is scattering from a potential . rather than using the square barrier of ref . @xcite , we consider the repulsive pschl - teller type potential  @xcite of the form @xmath46 since this potential does not have discontinuities the improved cn method works well with it . the transmission and reflection coefficients are known analytically . we can also compute them by considering the wavepacket eq . ( [ eq:4.5 ] ) incident on the potential . over a sufficiently long time the wavepacket will have interacted with the potential and transmitted and reflected packets emerge and travel away from the potential region . at that point we can calculate the probabilities of the particle represented by the packet on the left and on the right of the potential ; these probabilities correspond to the transmission and reflection coefficients provided the packet is sufficiently narrow in momentum space . this means that one needs an initial packet which is wide in coordinate space . in our calculation we choose @xmath47 , @xmath48 , @xmath49 , and @xmath50 . this gives a spread in the incident momentum - space wave packet of @xmath51 . the width in momentum space of the reflected and transmitted wavepackets also has this value . the domain of the @xmath30 coordinates is from @xmath52 to @xmath53 and the initial position of the packet is at @xmath54 to ensure that there is no overlap of the initial packet and the potential . we find good agreement between the transmission and reflections probabilities determined by plane - wave scattering approach and the time - dependent calculation as shown in fig . [ figure_1 ] . there are few analytically solvable models of the time evolution of unstable quantum systems  @xcite . realistic systems need to be solved numerically . long - time calculations are required for systems which require both nuclear and atomic time - scales such as ionization and bremsstrahlung due to radioactive decay of the nucleus of an atom  @xcite . to study the short - time anomalous power law behavior of the decay and the long - time inverse - power law behavior the method of this paper is appropriate . this is especially relevant because of the recently observed violation of the exponential - decay law at long times  @xcite .    to illustrate the numerical method discussed in this paper as applied to decaying systems let us consider a variant of the model with a @xmath55-shell potential  @xcite , but with the @xmath55 function replaced by a gaussian . thus in the @xmath56 partial wave of a spherical system the potential is @xmath57},\ ] ] where @xmath58 is the radial coordinate . this potential reduces to the @xmath55-shell potential , @xmath59 when @xmath60 . for small but finite values of @xmath61 this potential leads to scattering results which are good approximations of those of the @xmath55-shell interaction  . initially the quantum system is in the state @xmath62 in our example we take @xmath63 , @xmath64 , @xmath65 , @xmath66 and @xmath67 . using the numerical method of this paper , including the modification of matrix @xmath68 as described in sec . [ bc ] to take care of the boundary conditions at @xmath69 , we determine the wave function at later times , i.e. , @xmath70 . from that we obtain the nonescape probability , as a function of @xmath16 , @xmath71 , which is shown in fig . [ figure_2 ] . it clearly shows the exponential decay - region in time , the inverse power - law behavior for long times , the deviation from exponential decay at short times , and the transition regions  @xcite .    the quadratic short - time behavior is seen in fig . [ figure_2 ] as is the inverse power law behavior at long times  @xcite . remarkably the decaying system can be studied in this manner for a time exceeding thirty half - lives . in fig . [ figure_3 ] we plot the square of the absolute value of the wave function at times @xmath72 . notice that the wave function ( packet ) has three distinct regions : a precursor due to energy components of the initial wave function larger than that associated with the exponential decay , the main packet which corresponds to the exponential decay at the resonance energy , and the follower , which is a small blip that stretches in time ( travels more slowly ) and is due to energy components in the initial wave function which have lower energy that the resonance energy . if the maximum spatial coordinate , which is set at 800 for this calculation , were set at a smaller value , say 400 , then one observes a fuzziness in the precursor of the right - most wave function . this can be attributed to the finite space in which the wavepacket travels so that the fast precursor has been partially reflected from the right boundary and interferes with the wave front of the main wavepacket . the numerical parameters for this calculation are @xmath73 , @xmath74 , @xmath75 and @xmath76 . the generalized cn method that we have presented in this paper gives many orders of magnitude improvement in the precision of the results and several orders of magnitude in the computational time required to obtain the results . clearly , since this method enhances the efficiency of the numerical calculations , it can be a significant tool for studying time - dependent processes . it goes beyond the improvement of ref . @xcite in a systematic way . it also is an advance over the method of ref . @xcite , since the numerov method has an error @xmath10 , and it is difficult to see how it can be generalized systematically to higher order spatial errors . the generalized time evolution algorithm can be applied to moyer s  @xcite method . we have applied the method for @xmath2 and @xmath3 up to and including 20 for both . having achieved significant improvements with these values of @xmath2 and @xmath3 we did not consider larger values although there does not seem to be a practical reason that this can not be done . although the approach seems to saturate at @xmath2 around 10 ( see table  [ table3 ] ) , there is no visible saturation in the time - evolution part of the problem . one would expect higher orders of spatial errors to be significant when the wave function and/or the potential has large spatial fluctuations . even so , at @xmath40 we are using a 15-point formula for the second - order spatial derivative , and it is surprising that a smaller number of points in the formula is not sufficient for optimal results in this case . it should be noted that the higher order method as discussed in this paper are suitable only for well - behaved , sufficiently differentiable solutions ; these occur when the potential function is well behaved . as the authors of ref . @xcite point out for singular functions higher - order methods do not necessarily lead to greater accuracy . it should be noted that in this paper we consider primarily one - dimensional systems , but the method applies equally well to partial - wave equations of two- or three - dimensional systems . the study of the decaying quasi - stable state is an example of the latter . an interesting avenue to investigate further is the impact that this approach may have on two- or three - dimensional systems , where the number of variables involved is equal to the dimension . the peaceman - rachford - type approach  @xcite , also known as the alternating - direction implicit method , of factoring the approximation of the time evolution operator may apply as it did in refs . @xcite or more recently in , for example , refs . @xcite . using the crank - nicolson method the authors of refs . @xcite show that in two dimension the kinetic energy parts of the evolution operator factorizes . whether such factorization can be generalized in the spirit of the method of this paper is under investigation . preliminary studies with the numerov spatial integration scheme  @xcite and the generalized time evolution as described in this work indicate that significant improvements occur if one incorporates appropriate changes in the spatial step size for different regions of space . this is important in the case of discontinuous potentials and potentials that have great variation in some region and little or no variation in other regions . furthermore it is well known that the wavepacket has large fluctuation in a ( short - range ) potential region and little variation in the asymptotic regions . a great savings in computational time can be achieved by using different space - step sizes in the different regions . one needs to investigate whether such variable step size can be incorporated in the generalized spatial integration scheme of this paper . the approach that dealt with the discontinuous first- or second - order derivative in this paper and ref . @xcite is worth exploring . we intend to study this in the future . we are grateful to professor y. nogami for carefully reading the manuscript and the constructive comments , as well as useful and helpful discussions . one of the authors ( wvd ) expresses gratitude for the hospitality of the department of information and communication sciences of kyoto sangyo university where most of this work was completed . he also acknowledges the financial support for this research from the natural sciences and engineering council of canada and the japan society for the promotion of science .
Abstract:

Answer: we present a generalization of the often - used crank - nicolson ( cn ) method of obtaining numerical solutions of the time - dependent schrdinger equation . the generalization yields numerical solutions accurate to order @xmath0 in space and @xmath1 in time for any positive integers @xmath2 and @xmath3 , while cn employ @xmath4 . we note dramatic improvement in the attainable precision ( circa 10 or greater orders of magnitude ) along with several orders of magnitude reduction of computational time . the improved method is shown to lead to feasible studies of coherent - state oscillations with additional short - range interactions , wavepacket scattering , and long - time studies of decaying systems .


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: cellular automata ( ca ) have found many applications in the field of complex systems modeling . there is a number of works devoted to ca models for simulation of real - world phenomena like pedestrian dynamics @xcite , traffic flow @xcite , urban growth @xcite , etc . in most cases , the update rule and neighborhood for ca are determined by human experts that have knowledge of the modeled phenomenon . automatic identification of ca - based models remains an open research issue . several attempts have been made in the literature to develop algorithms for ca identification . however , most of the previous research did not investigate the use of real - world data sets as an input . the available algorithms were designed and tested mainly against synthetic data obtained from ca evolution . the use of real - world data for identification of ca models poses additional challenges due to inherent complexity of modeled phenomena and errors that are made during data acquisition . this paper discusses the possibility of using data mining techniques based on the rough sets theory @xcite to select neighborhood and determine rules for ca models . according to the proposed approach , input data describing observed states of cells are represented in form of a decision table . neighborhood is detected by using algorithms for reducts calculation . cells that do not belong to the neighborhood are removed from the decision table . after that , a rule - learning method is applied to induce a set of decision rules that define the evolution of ca . this approach was tested in identification of both deterministic and probabilistic ca models for synthetic as well as real - world data sets . most of the works related to the ca identification problem use genetic algorithms as a tool to extract update rule and neighborhood from spatio - temporal patterns produced in ca evolution @xcite . in @xcite a genetic algorithm was employed to learn probabilistic rules directly from experimental data for two - dimensional binary - state ca . several methods were proposed that use genetic algorithms to learn ca rules and neighborhood size for image processing tasks @xcite . although application of genetic algorithms was a dominant approach to ca identification , some non - genetic techniques are also available . adamatzky @xcite proposed several approaches to extracting rules for different classes of ca . straatman et al . @xcite developed a form of hill climbing and backtracking to identify ca rules . in @xcite a  deterministic sequential floating forward search method was used to select rules of ca . another approach is based on parameter estimation methods from the field of system identification @xcite . a framework for solving the identification problems for both deterministic and probabilistic ca was presented in @xcite . main drawbacks of the existing ca identification methods are related to the fact that they were either designed for a particular class of ca , or their experimental evaluation was limited to the case of synthetic data . the introductory research results presented in this paper shows that the rough sets approach may be effectively used to develop a universal method for identification of ca models that mimic real - world phenomena . formally , a cellular automaton can be defined as a triple @xmath0 , where @xmath1 is a non - empty set of cell states , @xmath2 is the neighbourhood , and @xmath3 is the update rule . arguments of @xmath3 are the current states of cells in the neighbourhood , while the value of @xmath3 indicates the state of a central cell at the next time step . the problem of ca identification involves finding both the cells neighborhood and the update rule on the basis of a training data set , which includes observed states of the cells . in order to use the rough sets approach for solving this problem , the training set of data has to be represented in the form of a decision table @xmath4 , where @xmath5 is a non - empty set of observations and @xmath6 is a set of cell states : @xmath7 . state of @xmath8-th cell at time step @xmath9 is denoted by @xmath10 . index @xmath11 indicates the central cell , for which neighborhood and update rule have to be found . thus , the cell state @xmath12 is used as the decision attribute . the remaining cell states @xmath13 are condition attributes . the candidate neighbors @xmath14 are the cells for which distance to the central cell @xmath11 is lower than a threshold value . the threshold is determined experimentally . above definition can be easily extended to multidimensional ca by adding the necessary indexes . neighborhood for the @xmath11-th cell can be selected by calculating reducts of the above - defined decision table . a reduct is a subset of the condition attributes , which is sufficient to determine the decision attributes @xcite . taking into account the decision table discussed above , reduct should be defined as a subset of cell states @xmath15 , which preserves discernibility of the observations with regard to the decision @xmath12 , and none of its proper subsets has this ability . observations are discernible if they differ in at least one condition attribute ( cell state ) . each two observations that have different decisions @xmath12 and are discernible by the full set of cell states @xmath16 are also discernible by the reduct @xmath17 .    the neighborhood @xmath2 of a cell is determined as a set of cells , whose states belongs to the reduct @xmath18 , where @xmath8 is a cell index . there may exist multiple reducts for one decision table . selection of the neighborhood is made with regard to the shortest reduct , because size of the neighborhood has to be minimized . if there are several minimal reducts then the one is selected which has the lowest average distance between neighbors and the central cell . when the reduct @xmath17 is found , the condition attributes that do not belong to this reduct are excluded from the decision table . thus , the modified decision table @xmath19 has the following set of attributes : @xmath20 , where @xmath12 remains the decision attribute .    update rule of a cellular automaton is identified as a set of decision rules by taking into account the information from the modified decision table @xmath21 . a particular decision rule @xmath22 has the following form :    @xmath23    where : @xmath24 , @xmath25 , and @xmath1 denotes the set of allowable cell states . two characteristics of decision rules are useful for the proposed method : support and match @xcite . support of rule @xmath22 , denoted by @xmath26 , is equal to the number of observations from @xmath21 for which rule @xmath22 applies correctly , i.e. , premise of the rule is satisfied and the decision given by rule is consistent with the one in decision table . @xmath27 is the number of observations in @xmath21 for which the rule @xmath22 applies , i.e. , premise of the rule is satisfied . based on these two characteristics , a certainty factor is defined for the decision rule @xmath22 :    @xmath28    the certainty factor may be interpreted as a conditional probability that decision of rule @xmath22 is consistent with an observation in the decision table @xmath21 , given that premise of the rule is satisfied . 1 summarizes the main operations that are necessary to identify a cellular automaton by using the rough sets approach . in this study , reducts and decision rules are calculated using algorithms implemented in the rses software @xcite . three algorithms of reducts calculation were examined : exhaustive , genetic , and dynamic reduct algorithm . moreover , the experiments involved application of three algorithms that enable induction of decision rules : exhaustive , genetic , and lem2 . in this section , the proposed approach is applied for identification of three deterministic ca @xcite : elementary cellular automaton with wolfram s rule 184 ( eca-184 ) , deterministic version of nagel - schreckenberg cellular automaton ( nash - d ) , and the conway s game of life cellular automaton ( life ) . eca-184 is a one - dimensional cellular automaton with binary cell states and neighborhood of three cells wide . original definition of the eca-184 update rule is presented in fig . 2 . the upper row in this figure illustrates all possible states of a central cell and its neighborhood . lower row shows states of the central cell after update  in the next time step of the ca evolution . identification of eca-184 was performed using a training data set of 500 observations . each observation in the decision table @xmath29 has covered a group of 21 cells ( candidate neighbors ) and the 11-th cell was considered as the central one . the neighborhood of eca-184 was correctly recognized by each of the reducts calculation algorithms ( exhaustive , genetic , and dynamic ) . the shortest reduct was determined as @xmath30 , thus the modified decision table @xmath21 had four attributes : @xmath31 . tab . 1 presents the decision rules that were calculated from table @xmath21 by using the lem2 algorithm . symbol @xmath32 indicates that a given cell state does not occur in a particular decision rule , e.g. , the rule no . 3 from tab . 1 should be interpreted as : @xmath33 . the set of decision rules in tab . 1 is consistent with the original update rule of eca-184 ( fig . the remaining algorithms of rule induction ( exhaustive and genetic ) have also generated valid sets of decision rules ; however their size was larger ( 6 rules ) . .decision rules generated for eca-184 [ cols="^,^,^,^,^,^",options="header " , ]     the unrealistic model behavior may occur because the decisions for uncertain rules are made randomly , without any coordination . e.g. , the rules 7 and 13 in tab . 5 apply to the same configuration of cells @xmath34 . if for both rules the random decision is 1 then the resulting configuration @xmath35 includes three instead of two occupied cells . this problem was resolved by adding error detection rules to the ca update algorithm . if error is detected by rules ( 4 ) then the update of cells @xmath36 is repeated . @xmath37    above modifications allow the ca to perform realistic simulation of traffic flow . 3 presents time space diagrams for two shock waves . the shock wave a ) was observed in real highway traffic and the shock wave b ) was obtained from evolution of the identified ca . this example illustrates the possibility of using the rough sets approach for identification of ca - based models of real - world phenomena . results of the introductory research show that the data exploration techniques based on rough sets theory enable proper selection of neighborhood and update rule induction for different classes of ca . it was also demonstrated that the rough sets approach is suitable for ca models identification from real - world data sets . the best results were obtained for the proposed ca identification procedure when implemented by using the dynamic reducts algorithm for neighborhood selection and the lem2 algorithm for decision rules induction . nevertheless , in case of stochastic ca identification , some additional actions were necessary ( reduct shortening , adding error detection rules ) . verification of the usability of the proposed approach for a wider set of data as well as comparison with state - of - art methods based on genetic algorithms remain open issues for further research . another interesting topic for future studies is to apply the proposed method in a data exploration system for ca - based image processing @xcite . ladniak m. , piorkowski a. , mlynarczuk m. : the data exploration system for image processing based on server - side operations . in : saeed , k. et al . cisim 2013 . lncs 8104 , pp 168 - 176 . springer , heidelberg ( 2013 )                  sun , x. , rosin , p. l. , martin , r. r. : fast rule identification and neighborhood selection for cellular automata . ieee trans . on systems , man , and cybernetics , part b : cybernetics , 41(3 ) , 749760 ( 2011 )
Abstract:

Answer: in this paper a method is proposed which uses data mining techniques based on rough sets theory to select neighborhood and determine update rule for cellular automata ( ca ) . according to the proposed approach , neighborhood is detected by reducts calculations and a rule - learning algorithm is applied to induce a set of decision rules that define the evolution of ca . experiments were performed with use of synthetic as well as real - world data sets . the results show that the introduced method allows identification of both deterministic and probabilistic ca - based models of real - world phenomena .