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: strongly correlated many - body systems are challenging due to their highly non - trivial interplay between different energy scales . in experiment , many of their properties are probed by measuring their linear response to weak external driving forces . an additional branch of physical phenomena is currently being explored by driving quantum many - body systems far out of equilibrium and studying such new states of matter . in particular , recent progress in atomic physics has made it possible to tune systems of ultracold atoms at will between different interacting regimes . for example , a seminal experiment by greiner et al . @xcite shows collapse and revival phenomena in atomic bose gases after an interaction quench , where excellent isolation from the environment allows the observation of such non - equilibrium behavior at remarkably long time scales .    from the theoretical side , non - equilibrium quantum systems are very challenging since many well - established methods from equilibrium physics can not be directly applied . significant progress has been achieved in low - dimensional systems and for quantum impurity systems with numerical methods like nrg and dmrg : their newly developed extensions called time - dependent density matrix renormalization group ( td - dmrg ) @xcite and time - dependent numerical renormalization group ( td - nrg ) , were , e.g. , successfully applied to interaction quenches in the bose - hubbard model @xcite and to impurity spins coupled to bosonic or fermionic baths . more recent applications of these methods also include interaction quenches in fermionic lattice systems in one dimension @xcite and the calculation of steady state currents through nano - devices @xcite . the extension of conventional analytical techniques seems to be more difficult . for example , the renormalization group approach , which is designed to construct effective low - energy hamiltonians , is not ideally suited to understand non - equilibrium situations , where energy scales well above the low - energy sector can be excited . the situation is not much better in weakly interacting systems , where diagrammatic expansions are usually the first tool to analyze ground states and their elementary excitations , and often yield reliable results in thermal equilibrium . nevertheless , even in weakly interacting systems perturbation theory is often not capable of describing the long - time behaviour of non - equilibrium observables , since truncation errors can become uncontrolled at large time scales .    in this paper , we present a new analytical approach to the real - time evolution problem which merges the advantages of perturbation theory and renormalization group theory , and at the same time , leaves behind their shortcomings mentioned above . in two short previous publications , this approach has already been introduced and was applied to two different problems . @xcite here we give a somehow more pedagogical introduction , provide more details and discuss its general applicability to the field of dissipative quantum systems . in the last twenty years , quantum physics in a dissipative environment has played an important role in solid state physics , quantum optics , quantum computing , chemical and even biological systems ( for a review , see refs . in addition , this field is especially suitable for checking the accuracy of our approximation scheme since many results from other approaches are known . however , our approach is applicable to a much wider class of problems , including also lattice models , see for example ref . @xcite . our approach is based on an analogy to canonical perturbation theory in classical mechanics . we give a simple illustration of canonical perturbation theory and show how canonical transformations can improve perturbative expansions in real - time evolution problems . we show that an analogous implementation for quantum many - body systems is possible , based on wegner s flow equation approach @xcite . independently , the same approach has been developed in the field of high energy physics @xcite . using this approach , we reproduce the exact solution for a quantum dissipative oscillator and show that efficient and precise numerical solutions of the analytical equations can be obtained . we also illustrate the failure of naive perturbation theory in this simple quantum mechanical system .    for the spin - boson model , we obtain results that are in excellent agreement with known results . in the regime of weak coupling to the bath , we reproduce exact results for the decoherence of a spin without tunnel splitting . for finite tunnel splitting , we calculate the non - equilibrium correlation function of the spin projection and also obtain the quantitatively correct coherent decay of the spin polarization in the weakly damped ohmic regime . in section [ feq ] we give an introduction into the basic ideas of our method and its technical details . first we motivate our approach using an analogy from classical mechanics . then we discuss the method in detail . in section [ dho ] , we apply this approach to a simple exactly solvable model which is useful to understand the technical details of our method . in section [ tss ] , we apply our method to the non - trivial spin - boson model to analyze the reliability of our approximation scheme . a brief summary and outlook concludes the paper in section [ disc ] . there exist well - established methods to handle time - dependent perturbations in classical mechanics , and early attempts to handle these problems date back already to newton , who considered the small distortion of the moon orbit caused by the gravitational force of the sun . much progress in this field has been motivated by the ever increasing accuracy of observational data for planetary motion and satellites , and the need to make accurate predictions based on this . the common approaches to these problems are collectively summarized as  canonical perturbation theory " . the basic idea behind canonical perturbation theory can be simply illustrated by the hamilton function of a weakly perturbed classical harmonic oscillator ,    @xmath0    where a quartic anharmonic term with a small coupling @xmath1 perturbs the trajectory . we consider the initial conditions @xmath2 and @xmath3 in the following .    by exploiting the smallness of the quartic perturbation , one might be tempted to employ naive perturbation theory which uses a series in powers of @xmath4 as an ansatz for the perturbed solution @xmath5 ,    @xmath6    the trajectory @xmath7 is just the solution of the unperturbed problem ,    @xmath8    from hamilton s equations ,    @xmath9    we obtain the equation of motion    @xmath10    which has the solution    @xmath11    this result already reveals the caveat of this approach , since the so called secular term @xmath12 yields an error growing unbounded in time . in fact , it is a well known general result from classical mechanics that such secular terms invalidate naive perturbation theory for large times . naive perturbation theory can be much improved in the framework of canonical perturbation theory . this approach first transforms the hamilton function to a suitable normal form , and after solving the equations of motion in this canonical basis , the normal coordinates are reexpressed through the old coordinates . in this manner secular terms can be avoided . to implement this idea , we first look for a canonical transformation of variables    @xmath13    that brings the hamilton function to the following normal form , denoted by @xmath14 ,    @xmath15    it is easy to see that the poisson bracket @xmath16 vanishes , hence the equations of motion for @xmath17 and @xmath18 can be solved trivially . these variables obey the initial conditions @xmath19 and @xmath20 . in our example , the corresponding transformation of variables is    @xmath21    which is performed perturbatively in the parameter @xmath4 . the normal form ( [ normal ] ) has been chosen such that the equation of motion for the new variables @xmath22 and @xmath23 can now be solved exactly , without producing any secular term . using this strategy , the final result is    @xmath24    where @xmath25 and @xmath26 . we show a comparison of naive perturbation theory and canonical perturbation theory in fig . [ canonicalpic ] , which demonstrates the usefulness of canonical perturbation theory . we compare the different approaches to solve the equations of motion for the anharmonic oscillator from eq . ( [ hamilton ] ) . the difference between the numerically exact solution and canonical perturbation theory according to eq . ( [ canonicalresult ] ) can hardly be noticed . naive perturbation theory yields large errors already after a few oscillations , with an error that grows linear in time @xmath27 . our parameters are @xmath28 and coupling strength @xmath29 . , width=302 ]    canonical perturbation theory yields renormalized parameters of the unperturbed problem , but does not lead to secular terms . this improvement becomes directly visible by expanding the contribution @xmath30 in powers of  @xmath4 ,    @xmath31    this expansion generates the secular term from eq . ( [ naiveresult ] ) occuring in naive perturbation theory , making it obvious that canonical perturbation theory contains a summation over secular terms which in total yield a much improved result in comparison to naive perturbation theory . notice that this is true in spite of the fact that _ both approaches have been expanded to the same power in the small coupling constant_. +      in quantum many - body systems , the canonical way to evaluate the real time evolution of observables starting from some non - thermal initial state is the keldysh technique @xcite , which defines a contour ordered s - matrix in order to develop perturbative expansions for non - equilibrium greens functions . just as in our example of naive perturbation theory , secular terms can occur in any finite order of perturbation theory , and there is no universal solution for how to sum up these secular terms . these difficulties can make it very difficult or even impossible to study the transient evolution of observables into a steady state . one might wonder why an analogue of canonical perturbation theory for quantum many - body systems had not been developed earlier . the basic reason is the notorious difficulty to transform quantum - many body systems into normal form . in addition , the continuum of energy scales often causes often non - perturbative effects in coupling constants . since the advent of renormalization group techniques , the much easier problem of constructing effective low - energy theories has attracted most of the attention . the more general problem of constructing normal forms of interacting hamiltonians has been successfully treated only during the last years . how to make use of the flow - equation method to implement the analogue of canonical perturbation theory in quantum mechanics . u denotes the full unitary transformation that relates the @xmath32 to the @xmath33 basis . , width=120 ]    a general description how to transform interacting quantum many body systems into a non - interacting normal form has been given in 1994 by f. wegner.@xcite let us briefly review the basic ideas of the flow equation approach ( for more details see ref.@xcite ) . a many - body hamiltonian  @xmath34 is diagonalized through a sequence of infinitesimal unitary transformations with an anti - hermitean generator  @xmath35 , @xmath36 \ , \label{eqdhdb}\ ] ] with @xmath37 the initial hamiltonian . canonical " generator @xcite is the commutator of the diagonal part  @xmath38 with the interaction part  @xmath39 of the hamiltonian , @xmath40 $ ] . it can be formally shown that the choice of the canonical generator is by construction suitable to eliminate interaction matrix elements with energy transfer @xmath41 . under rather general conditions , an increasingly energy diagonal hamiltonian is obtained . for @xmath42 the hamiltonian will be energy - diagonal and we denote parameters and operators in this basis by  @xmath43 , e.g.  @xmath44 .    usually , this scheme can not be implemented exactly . the generation of higher and higher order interaction terms in ( [ eqdhdb ] ) makes it necessary to truncate the scheme in some order of a suitable systematic expansion parameter ( usually the running coupling constant ) . still , the infinitesimal nature of the approach makes it possible to deal with a continuum of energy scales and to describe non - perturbative effects . this had led to numerous applications of the flow equation method where one utilizes the fact that the hilbert space is not truncated as opposed to conventional scaling methods .    in ref . @xcite , these features have been exploited in order to develop an analogue of canonical perturbation theory in classical mechanics for quantum many - body systems . the general setup is described by the diagram in fig . [ forwardback ] , where @xmath45 is some initial non - thermal state whose time evolution one is interested in . however , instead of following its full time evolution it is usually more convenient to study the real time evolution of a given observable a that one is interested in . this is done by transforming the observable into the diagonal basis in fig . [ forwardback ] ( _ forward transformation _ ) :    @xmath46 \label{operatorflow}\ ] ]    with the initial condition @xmath47 . the key observation is that one can now solve the real time evolution with respect to the energy - diagonal @xmath48 exactly , thereby avoiding any errors that grow proportional to time ( i.e. , secular terms ) : this yields @xmath49 . now since the initial quantum state is given in the @xmath50 basis , one undoes the basis change by integrating ( [ operatorflow ] ) from @xmath33 to @xmath32 ( _ backward transformation _ ) with the initial condition @xmath51 . one therefore effectively generates a new non - perturbative scheme for solving the heisenberg equations of motion for an operator , @xmath52 , in exact analogy to canonical perturbation theory . the dissipative harmonic oscillator is a widely used toy model in the field of quantum optics and is also used in many other contexts.@xcite it describes a quantum oscillator of frequency @xmath53 coupled linearly to a heat bath consisting of bosonic normal modes @xmath54 . @xmath55    the operators @xmath54 fulfill canonical commutation relations    @xmath56= \delta_{kk^\prime}\ ] ]    and their influence on physical properties of the quantum oscillator can be fully described by the spectral function    @xmath57    in experiment , the high frequency part of @xmath58 only affects the short time- response of the system , and thus , this part of the spectrum is usually cut off from @xmath58 . in consequence , the function @xmath59 is commonly approximated by a power - law behavior @xmath60 . three different regimes of the exponent are distinguished , for @xmath61 the bath is called  subohmic " , for @xmath62 the bath is called  ohmic " and for @xmath63  superohmic " . we imagine that the system is prepared in a well - defined quantum state at some time @xmath64 and the subsequent real time evolution of physical quantities is then defined by the hamiltonian ( [ dhomodel ] ) . the flow equation method does not restrict the class of possible initial states , but in this work we will consider the system - bath complex at time @xmath65 to be prepared in a product state    @xmath66    with the quantum oscillator in a coherent state    @xmath67    and the heat bath in the bosonic vacuum state @xmath68 . in such a state , the displacement @xmath69 will be finite and the effects of decoherence and dissipation will manifest themselves in the real time evolution of the observable @xmath70 . the flow equation method , outlined in section [ feq ] , can solve this problem exactly without approximations . first , the flow equations for the hamiltonian are derived . then the real time evolution of the operators @xmath71 and @xmath72 is implemented by their time - dependent flow equations . this will allow us to study the forward - backward transformation scheme in the context of an exactly solvable model . as the first step of our program , the coupled form of the hamiltonian ( [ dhomodel ] ) will be diagonalized by infinitesimal unitary transformations . due to the quadratic nature of the hamiltonian ( [ dhomodel ] ) , the flow equations for the hamiltonian can be derived in closed form , as shown in ref . . commuting the interaction part with the non - interacting boson part of ( [ dhomodel ] ) yields the canonical generator @xmath73 $ ] of unitary transformations :    @xmath74    this generator leads to additional terms in the flowing hamiltonian that violate the form invariance of equation ( [ dhomodel ] ) . the flowing hamiltonian preserves its initial form if the generator is extended by additional terms . in the following , we omit the explicit dependence of coefficients on the parameter @xmath75 .    @xmath76    with the coefficients    @xmath77    which are chosen such that they leave the flowing hamiltonian form invariant . for this purpose the function @xmath78 is still arbitrary , except for the obvious requirement @xmath79 . for our numerical evaluations of the flow equations later on , we will chose @xmath80 , which leads to good convergence properties in the limit @xmath81    the flowing parameters of the transformed hamiltonian are governed by the coupled differential equations    @xmath82    the renormalization of the bath frequencies @xmath83 will have a vanishing effect on time - dependent observables in the thermodynamic limit @xmath84 . therefore , the flow equations of these energies can be ignored for our purposes . in the limit @xmath85 , the hamiltonian is diagonalized and the tunneling matrix element @xmath53 is renormalized to some value @xmath86 . @xmath87    in order to proceed with the real - time evolution of observables , the second step of the program outlined in section [ feq ] requires the analogous transformation of the system operators . the flow of the bosonic operator @xmath88 is determined by the generator ( [ gendho ] ) and yields the following structure    @xmath89    with the flow equations    @xmath90    the initial conditions are @xmath91 . during the flow towards @xmath85 , the operator @xmath71 changes its structure into a complicated superposition of bath operators . it is now easy to formulate an exact transformation that yields the operator @xmath92 , which is time - evolved with respect to the hamiltonian ( [ dhomodel ] ) . first , the operator @xmath93 is trivially time evolved with respect to the hamiltonian ( [ hinfty ] ) as @xmath94 this operation endows the coefficients in ( [ bansatz ] ) with trivial phase factors    @xmath95    leading to the operator    @xmath96    the second step is to obtain the operator b(t ) from the operator @xmath97 by reverting the unitary transformation @xmath98 can be expressed as a @xmath75-ordered exponential , @xmath99 . however , this expression is only formally useful since it can not be evaluated without additional approximations . ] used to diagonalize the hamiltonian ( [ dhomodel ] ) , formally represented by the relation @xmath100 . for this purpose , we again make an ansatz for the flow of the operator @xmath101    @xmath102    where all coefficients have both real and imaginary part , since the initial conditions at @xmath103 are given by the complex valued expressions from ( [ initialb ] ) . since the ansatz of eq . ( [ bbt ] ) is formally identical to that of ( [ bansatz ] ) , the unitary flow of @xmath104 can again be calculated by using the flow equations ( [ bfeq ] ) . the operator @xmath105 is now obtained by integrating the flow equations ( [ bfeq ] ) from @xmath103 to @xmath50 , using the parameters from ( [ bbt ] ) and the initial conditions for @xmath103 , as given by eq . ( [ initialb ] ) . since all transformations are unitary , the operator @xmath106 is readily obtained by the hermitean conjugate of eq . ( [ bbt ] ) . all transformations used up to now do not depend on the initial state of the quantum system . in calculations of time - dependent physical quantities , the operator @xmath105 can be evaluated with respect to equilibrium heat baths at finite temperature as well as arbitrary non - equilibrium ensembles .      _ 1 . naive perturbation theory _ + in analogy to classical perturbation theory as discussed in the introduction , a perturbative result for quantum evolution can be obtained by directly expanding the heisenberg equation of motion for an operator . the exact time evolution of the operator @xmath107 is    @xmath108    where @xmath34 is given by eq . ( [ dhomodel ] ) . a perturbative expansion of eq . ( [ heisenx ] ) is typically performed in the interaction picture , which we define by    @xmath109    and @xmath110 is the non - interacting part of the hamiltonian . ( [ diracx ] ) leads to the equation of motion    @xmath111\ ] ]    which can then be expanded as    @xmath112 \nonumber\\                & + & i^2   \int_0^t d\tau_1 \int_0^{\tau_1 } d \tau_2[h _   { int}^i(\tau_1),[h_{int}^i(\tau_2),\hat{x } ] ] \nonumber\\                & + & o(h_{int}^{i})^3 \ . \label{expanded}\end{aligned}\ ] ]    here @xmath113 is the interaction part of the hamiltonian in the interaction picture . neglecting all contributions of @xmath114 or higher in the perturbative expansion will yield a result that is correct at least to @xmath115 . using the inital state ( [ coherentstate ] ) , the first order contribution from eq . ( [ expanded ] ) vanishes , since @xmath116 . likewise , all odd orders of perturbation theory vanish . evaluation of the two commutators needed for the second order contribution to @xmath117 yields    @xmath118    in analogy to our example of naive perturbation theory in classical mechanics from section [ feq ] , a secular term occurs in the expansion from eq . ( [ xtnaive ] ) which leads to large errors on long time scales ( see fig . [ expansion ] ) . exact solution _ + we next derive a closed analytical expression for the quantity @xmath119 as a benchmark for a full numerical solution of the flow equations . several other exact results for the dissipative quantum oscillator have been obtained by haake and reibold.@xcite according to eq . ( [ bbt ] ) , we have @xmath120 . the coefficient @xmath121 is also given by the commutator    @xmath122_- , \label{commute}\ ] ]    and using the invariance of ( [ commute ] ) under the unitary flow ( [ bansatz ] ) , we obtain    @xmath123\nonumber\\ & = & 2 \sum_k ( \tilde{\alpha}_k^2 + \tilde{\bar{\alpha}}_k^2 ) \cos   ( \omega_k t).\end{aligned}\ ] ]    it is possible to evaluate this sum in closed form by defining    @xmath124    what leads to the exact result for the dynamics    @xmath125    for an ohmic bath with @xmath126 , the function @xmath127 can be evaluated as    @xmath128 ^ 2 \biggr)^{-1 } .\end{aligned}\ ] ]      in order to illustrate the potential applications of our diagonalization scheme , we calculate the expectation value @xmath129 by a numerical integration of the flow equations eq . ( [ bfeq ] ) . for our numerical calculation , we specify the spectral function of the bath as an ohmic one with a sharp cutoff frequency @xmath130 . @xmath131    this spectral function is discretized with @xmath132 states with a constant energy spacing of @xmath133 . the systems of coupled differential equations ( 23 ) and ( 26 ) have to be solved separately for each point in time , since the initial conditions of eq . ( [ initialb ] ) depend on time . for a finite number @xmath134 of bath modes , the initial conditions for the bath modes have a recurrence time of @xmath135 , and it is expected that the error blows up rapidly at this time scale . indeed , we could confirm this behavior numerically . nevertheless , for times @xmath136 , a number of about @xmath132 bath modes is sufficient to agree with the exact result within an error of less than @xmath137 . + summing up , numerical simulations with only @xmath138 bath states provide excellent agreement with analytical solutions , with finite size effects occuring only at time scales of @xmath139 . in [ feqs ] , we briefly discuss how to efficiently implement our transformation scheme for finite size systems . comparison of the time - dependent displacement @xmath129 of the dissipative harmonic oscillator , obtained from the analytical result of eq . ( [ analytical ] ) ( continuous lines ) and numerical integration of the flow equations ( crosses ) . different damping strengths @xmath140 have been chosen , and the tunneling matrix element @xmath53 is almost renormalized to zero at @xmath141 , leading to a much slower oscillation period . the curve has been normalized to 1 by chosing @xmath142 . the comparison demonstrates that a numerical solution of real time observables with about 1000 bath states reproduces the analytical result already excellent with an error below @xmath143 . , width=302 ] the dissipative two - level system    @xmath144    is a fundamental model for the description of decoherence and dissipation in quantum systems.@xcite the two- state system is represented by pseudospin operators @xmath145 and the effect of dissipation is caused by a linear coupling to a bosonic bath . all bath properties can again be modeled by the spectral function @xmath59 defined in eq . ( [ spectrum ] ) .    typically , the effect of decoherence manifests itself if the two- level system is prepared in an eigenstate and the coupling to the environment is switched on subsequently . the observable @xmath146 describes then the tunneling dynamics of the initially pure quantum state . this problem turns out to be non - trivial and nearly all known results rely on approximations . the approximation scheme of the flow equation approach can be applied to this problem in a controlled way as shown below . this example demonstrates the ability of our method to treat real - time evolution in non - integrable models without any problems due to secular terms .      in order to approximately diagonalize the hamiltonian ( [ tls ] ) , we employ the following generator for the unitary flow:@xcite    @xmath147    with @xmath75-dependent coefficients    @xmath148    this generator has an important conceptional property that makes it different from the generator ( [ gendho ] ) used in section [ dho ] . ) is again not the canonical generator in order to leave the flowing hamiltonian form invariant up to higher orders in the coupling @xmath149 . it also makes use of the approximation @xmath150 and neglects small fluctuating parts of this expectation value . ] it does not leave the hamiltonian form invariant during the flow , since it generates additional interactions . in this case , these newly generated interactions are of @xmath115 , which we neglect since we consider the couplings @xmath149 as a small expansion parameter . for more details about the flow equation approach to the spin - boson model we refer the reader to ref . @xcite . due to the expansion in the couplings @xmath149 , the flow equation calculation in equilibrium is reliable ( meaning errors less than 10% ) for values @xmath151 for an ohmic bath .    in the above expressions , normal ordering is denoted by @xmath152 , which ensures that the truncated higher order interaction terms have vanishing expectation values with respect to the quantum state used for normal - ordering . in equilibrium , normal - ordering is therefore performed with respect to the equilibrium ground state , @xmath153 , where @xmath154 is the bose - einstein distribution . this procedure is not ideal for real - time evolution of physical observables out of a non - thermal initial state @xmath155 . in order to minimize our truncation error , we define a more general normal - ordering procedure    @xmath156    here , the correct non - thermal initial state is used for normal - ordering . the flow equations for the hamiltonian then read :    @xmath157    in the limit @xmath85 , the interaction part of the hamiltonian decays completely and the transformed hamiltonian is given by    @xmath158    with a renormalized tunneling matrix element @xmath86 . it has been shown in ref . @xcite that @xmath159 obeys the correct universal scaling behavior for ohmic baths@xcite , @xmath160 .      the truncation scheme for the flow of the hamiltonian ( [ tls ] ) can be employed in the same way for the transformation of the spin operators @xmath161 . in this section , only the transformations of the operator @xmath162 are presented . details of the transformations of the operators @xmath163 and @xmath164 are given in appendix [ feqs ] . an ansatz for the flow of the operator @xmath162 is formally given by the commutator @xmath165 $ ] , which contains all contributions to the flow of @xmath162 which are of first order in the couplings @xmath149 . for convenience , we parametrize this ansatz for the flowing operator @xmath166 as    @xmath167    where @xmath168 and @xmath169 are related to @xmath170 and @xmath171 by complex conjugation . all newly generated terms in the differential equation @xmath172 $ ] are of @xmath173 in the coupling constants . these truncated terms are normal - ordered according to the convention ( [ normalorder ] ) , which finally determines the differential flow of the coupling constants as    @xmath174    with the initial conditions @xmath175 , @xmath176 . in the thermodynamic limit , the observable @xmath162 decays completely , @xmath177 .    for the application to real - time evolution problems , it is again straightforward to obtain the time - evolved operator @xmath178 by evaluating @xmath179 with the diagonal hamiltonian @xmath14 from ( [ freetls ] ) . it is easy to see that the transformed observable ( [ sigmaxb ] ) remains invariant under time evolution , and only its coefficients change to time- dependent functions    @xmath180    where @xmath181 and @xmath182 remain unchanged since these contributions commute with @xmath14 . the operator @xmath178 with coefficients ( [ coefft ] ) can be regarded as an effective operator for the calculation of observables in real - time with an error of @xmath173 . although this effective operator relies on a perturbative expansion of the flow equations , it is able to correctly describe observables on all time scales , e.g. the error remains controlled also for long times . this property can be understood from the analogy to canonical perturbation theory of classical mechanics , where the hamiltonian function is first transformed to normal form . ( [ feqsigmax ] ) together with the initial conditions ( [ coefft ] ) , the effective operator @xmath178 can be integrated back into the inital basis of the problem , thereby inducing a non - perturbative solution of the heisenberg equation of motion for the operator @xmath162 , as illustrated in fig . [ forwardback ] . formally , this solution for the operator @xmath183 is given as    @xmath184    after the coefficients @xmath185 have been obtained , e.g. by a numerical solution of the flow equations , any desired correlation function of the operator @xmath186 can be calculated . although it is in principle possible to obtain analytical approximations to flow equations like ( [ feqsigmax ] ) and ( [ feqtls ] ) , we solve the flow equations for the spin operators numerically in this paper . however , these numerical solutions rely on the truncation scheme we introduced above . in order to check the accuracy of this truncation scheme , we first analyze it by comparing it against an exact result . comparison against exact results _ + an exact analytical solution of the spin - boson model can be easily obtained at the so called pure dephasing point , where @xmath187 . since @xmath164 is a conserved quantitity in this case , the environment can be traced out analytically . an initial state that contains locally a pure state @xmath188 will become entangled with the environment during the course of time . a measure for this process is given by the decay of the off - diagonal matrix element of the reduced density matrix , @xmath189 , and it was shown @xcite that this can be written as @xmath190 with    @xmath191    note that the quantity @xmath192 is identical to the observable @xmath193 , which can be directly obtained from eq . ( [ effsigx ] ) . at zero temperature this observable shows a sluggish decay to the ground state expectation value @xmath194 . we now compare this exact result with the flow equation solution at zero temperature . for the given initial state , the quantity @xmath192 is in the flow equation scheme given by the real function @xmath195 from eq . ( [ effsigx ] ) . we evaluate this quantity for ohmic baths with a cutoff frequency @xmath130 ( superohmic baths are also possible , subohmic baths have not yet been treated successfully within our approach ) . the calculations show that in the regime of low damping strengths @xmath196 the agreement with the exact result is excellent , see fig . [ decoherence ] . deviations from the exact result are on all time scales bounded by corrections of @xmath197 . the numerical results also show the small transient oscillations of the analytical result , which oscillate with approximately the cutoff frequency . this can be interpreted as a band edge effect which would vanish if a smooth cutoff function like @xmath198 were used for the bath spectrum . notice that a numerical solution of the flow equations at finite temperatures ( not shown here ) is straightforward and also in good agreement with the exact result ( [ exact ] ) . +     decay of the off - diagonal matrix element @xmath192 of the reduced density matrix of the spin - boson model without tunnel splitting ( @xmath187 ) at @xmath199 . the dashed and dotted lines represent the exact analytical result given by ( [ exact ] ) , and the full lines the numerical solution of the flow equations . the result demonstrates that the truncation scheme for the flow equations indeed is valid on all time scales , with an error that is systematically controlled in @xmath197 . the bath spectral function has been discretized with @xmath200 states and cut off at @xmath201.,width=302 ]    _ 2 . decay of a polarized spin _ +     real time evolution of the spin expectation values @xmath202 starting from a polarized spin in z - direction with a relaxed ohmic bath . we used a damping strength @xmath203 and the parameters @xmath204 , @xmath205 , @xmath206 and @xmath207.,width=302 ]    the formulation of the spin - boson model was originally motivated by the so - called quantum tunneling problem . @xcite in this problem , the two - state system is initially prepared in an eigenstate @xmath208 of the pseudospin - operator , corresponding to a localized state in a double - well potential , which can be reduced to an effective two - state system . the bath , which we again take as ohmic with a frequency cutoff @xmath130 , is initially in thermal equilibrium and denoted by @xmath209 . thus , the initial state @xmath155 of the total system is given by the product state    @xmath210    after switching on the coupling to the dissipative environment , the time - dependent tunneling dynamics between the two possible states of the two - state system is described by the observable    @xmath211    within the flow equation approach this is given by the coefficient @xmath212 from the ansatz ( [ flowyz ] ) in appendix [ feqs ] . we have numerically solved the flow equations ( [ feqyz ] ) corresponding to the observable ( [ tunneling ] ) and depicted the result in fig . [ magnetization ] . as expected , our solutions show long - time stability without secular terms analogous to canonical perturbation theory . for intermediate times , our curves agree well with the well - established niba approximation , @xcite and for long - time scales @xmath213 the expectation value @xmath214 vanishes as expected . + _ 3 . non - equilibrium correlation functions _ +     the two - time non - equilibrium correlation function @xmath215 for two different damping strengths @xmath216 and @xmath203 . the full and the dashed line show the result for @xmath217 and the dotted lines correspond to @xmath218 . parameters : @xmath205 , @xmath219 , @xmath206 , @xmath220.,width=302 ]    without any additional effort we can also calculate two - time correlation functions based on our non - perturbative solution of the heisenberg equations of motion for the spin operators . as an example we discuss the non - equilibrium correlation function    @xmath221    where @xmath222 is the waiting time for the first measurement after switching on the dynamics . notice that in thermal equilibrium this correlation function is time translation invariant and does not depend on the waiting time @xmath223 . it is the relevant quantity for , e.g. , neutron scattering experiments @xcite . experimental results for this correlator were for example obtained for hydrogen trapped by oxygen in niobium @xcite , where protons can tunnel between two trap sites of the @xmath224 sample .    for studying the non - equilibrium dynamics of the spin - spin correlation function , we again use the quantum state @xmath225 from ( [ initialstate ] ) with the bath at zero temperature as the initial state . within the flow equation approach , @xmath226 is readily evaluated as    @xmath227.\nonumber\\\end{aligned}\ ] ]    all coefficients are explained in appendix [ feqs ] . it turns out from numerical calculations ( see fig . [ correlation ] ) that the dependence on @xmath223 is only weak in the limit of small damping strengths and the correlations are very close to the equilibrium behavior . our study of real - time dynamics in dissipative quantum systems was motivated by developing an analogous method to canonical perturbation theory as known from classical mechanics . our results show : ( i ) the unitary transformation scheme of the flow - equation approach reproduces the well - known advantages of canonical perturbation theory from classical mechanics , especially the absence of secular terms in the time evolution . for example in the spin - boson model , physical observables can be calculated reliably on all time scales since perturbation theory is performed in a unitarily transformed basis . ( ii ) the identification of a suitable expansion parameter is crucial to obtain reliable results . notice that the implementation of the flow equation scheme is not restricted to bosonic baths or to impurity models , see e.g. refs . ( iii ) our method is particularly suitable for studying different initial states , since all transformations and approximations are performed on the operator level . we acknowledge the financial support through sfb / tr  12 of the deutsche forschungsgemeinschaft . a.h . also acknowledges support through sfb  608 and s.k . acknowledges support through the center for nanoscience cens munich , and the german excellence initiative via the nanosystems initiative munich nim . _ 1 . derivation of flow equations _ + we briefly provide the flow equation transformations for the spin operators @xmath163 and @xmath164 , that are constructed in complete analogy to the example from section [ tss ] . using the generator ( [ gentss ] ) , the anstze for the flowing spin components @xmath163 and @xmath164 read    @xmath228    the flow equations for these operators are readily derived as    @xmath229    and    @xmath230    next , we solve the heisenberg equations of motion for the operators @xmath231 and @xmath232 , which have the formal solution @xmath233 , with the hamiltonian @xmath14 given by eq . ( [ freetls ] ) . the solutions read :    @xmath234    these operators differ formally from those of eq . ( [ ansatze ] ) , therefore we chose a different ansatz for the transformation of these operators ,    @xmath235    which is for both @xmath163 and @xmath164 identical . this ansatz yields the time - dependent flow equations    @xmath236    with the initial conditions    @xmath237    for @xmath238 and    @xmath239    for @xmath240 . the effective solutions @xmath241 of the respective heisenberg equations of motion are now obtained by determining eq . ( [ flowyz ] ) in the case @xmath50 . numerical implementation of flow equations _ + in order to numerically integrate a set of @xmath242 coupled differential equations , an adaptive step size fourth order runge kutta algorithm is a fast and accurate choice . discretizing a bosonic or fermionic bath with equal energy spacing and using about @xmath132 bath states yields very accurate solutions up to time scales of @xmath243 . in order to determine initial conditions like eq . ( [ initials ] ) , it is sufficient to integrate the flow equations up to @xmath244 , where only a fraction of @xmath245 of the couplings @xmath149 has not decayed exponentially yet . in order to integrate flow equations with inital conditions like eq . ( [ initials ] ) , a standard runge kutta algorithm works not properly for large values of the parameter @xmath75 , since the exponential smallness of the couplings exceeds floating point precision . in our implementation , we stored therefore the flow of the couplings from the diagonalization of the hamiltonian and supplemented it to the integration of the flow equation with time - dependent inital conditions . although this procedure can not use an adaptive stepsize in order to control the error during integration , it turned out to be very precise in all cases we used it .
Abstract:

Answer: we discuss a new analytical approach to real - time evolution in quantum many - body systems . our approach extends the framework of continuous unitary transformations such that it amounts to a novel solution method for the heisenberg equations of motion for an operator . it is our purpose to illustrate the accuracy of this approach by studying dissipative quantum systems on all time scales . in particular , we obtain results for non - equilibrium correlation functions for general initial conditions . we illustrate our ideas for the exactly solvable dissipative oscillator , and , as a non - trivial model , for the dissipative two - state system .


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: results obtained in the solar , atmospheric and reactor experiments are interpreted as the evidence of the neutrino oscillations and the existence of a small ( @xmath2  ev ) neutrino mass . but absolute values of neutrino masses and their hierarchy can not be determined solely from oscillation experiments . experimental searches for the neutrinoless double beta decay enable to check up the law of lepton charge conservation , to define the neutrino mass nature ( dirac or majorana ) and to investigate the absolute scale of neutrino masses  @xcite . evidence of the observation of @xmath3ge @xmath4-decay  @xcite gave rise to a great discussion . thus appears the necessity to check this result and its interpretation with other isotopes and with higher sensitivity . a lot of new projects of the experiments with few tons isotopes are under consideration and development at present , for example genius  @xcite with @xmath3ge and exo  @xcite with @xmath0xe . the goal of the future experiments is to reach of the sensitivity @xmath5 years and to a neutrino mass of the order @xmath6  ev . however , for unequivocal definition on a neutrino mass ( or limit on mass ) the knowledge of a nuclear matrix element is necessary . there are many publications , devoted to theoretical calculations of the nuclear matrix elements both for @xmath7 and @xmath4 decay modes . results differ more than 2 times . therefore , experimental observations of @xmath7-decay for different isotopes with maximal accuracy are extremely important . possible test of the theoretical calculations has been proposed in  @xcite . where it was stated that final test of the theoretical calculations can be accomplished if @xmath4-decay of three ( or more ) isotopes will be observed . @xmath3ge @xmath7-decay was measured in the few experiments  @xcite , though some small disagreement exists . @xmath0xe @xmath7-decay is not observed up to now , so existing limits for half life time greatly differ , from @xmath8  y to @xmath9  y  @xcite . large track detector devis ( detector visualizing ) is the time projection chamber ( tpc ) in the magnetic field and it is operating at itep  @xcite . the first goal of the experiment is the investigation of @xmath0xe @xmath7-decay . schematic view of the setup one can see on the fig . [ pic_setup ] . tpc is placed at the magnet with diameter @xmath10 m. magnetic field can be changed in the range @xmath11 kilogauss with homogeneity not worse than 2% . source of @xmath12-decay is @xmath0xe , located in the central part of tpc , separated by thin mylar films ( 50  mkm ) out of two adjacent parts , filled by methane . detector is operating at atmospheric pressure . tpc detector has two drift sections , upper and lower ( high voltage drift electrodes are located on the upper and lower plates ) . uniformly distributed potential from 35  kv to 0  kv is fixed on the all - lateral surfaces ( including mylar films ) using special electrodes , creating the homogeneous electric drift field . ionization from passed charge particles in the gas is drifted in direction of an electric field to the proportional chambers , arranged in the center of tpc . three multiwire chambers ( for each gas volume ) have two rows of sensitive wires registering the events in the upper and lower drift sections . signals from two sensitive wires ( the upper and the lower ) are summarized on the input of preamplifier . detector has @xmath13 sensitive wires . additional readout from stripes ( cathode wire electrodes formed @xmath14 related to the sensitive wires ) is made for methane volumes  @xcite . the coordinate measurements are carried out with tdc ( time digital converter ) with number of the sensitive wire being coordinate along the magnetic field while the drift time being the coordinate along the electric field . the trajectory of the electron in magnetic field has a spiral form . therefore in the basic projection ( number of sensitive wire  drift time ) track looks like a sinusoid , and in the additional projection ( number of cathode strip  drift time ) like a cycloid . the momentum and motion direction are defined by fitting of the trajectory . block scheme of the trigger and recording electronics one can see on fig .  [ pic_el ] . signals from preamplifiers of sensitive wires and cathode strips [ a ] with 13  m cables go to discriminators [ d ] with threshold 30  mv ( 15  mv ) with standard nim ( @xmath15  v ) output . the discriminator can work in three modes : duration above threshold , formation of the duration , start and end of the signal from the wire or strip . we mainly used duration above the threshold mode . thus duration on an output of the discriminator corresponds to the duration of a signal exceeding the threshold . dead time of the device was equal to 60  ns . discriminators have two outputs . each wire and strip was connected to tdc using twisted pairs and the time of each hit from wires or strips was fixed . other outputs from discriminators of sensitive wires were used to form trigger . parameters of tdc ( full time of tdc circle  80  mks , one bin  20  ns ) were chosen to minimize the errors in spatial resolution due to electronics , compared to errors by multiple scattering and diffusion . since it is unknown when an event happens the setup worked in mode `` run '' . the clock generator was started by computer , and stopped by trigger . trigger started information recording from tdc to computer . the signal `` fast dump '' [ fd ] cleared the memory of all tdc if trigger did not occurred during 320  mks , 4 cycles of tdc . fd was blocked if it was coincided with trigger , so we had not losses due to fd . on the first stage we used 320 tdc blocks  @xcite . the basic part of trigger electronic consists of two trigger processors m@xmath16 and m@xmath17 connected to left and right volumes of tpc . the trigger signal from processor arises , when more than 20 sensitive wires from 32 wires on the inputs of m@xmath18 , coincide in the 6 mks gate . this time corresponds to the size of the sinusoid track and drift velocity in the ch@xmath19 volume . the signals from m@xmath16 and m@xmath17 coincide at the scheme [ cc1 ] . at the first stage we used as a trigger coincidences between 4 fixed wires ( 2 wires from both left and right ch@xmath19 volumes ) . use of the processors increased the trigger efficiency from large drift distance because the register efficiency of separate sensitive wires is equal to 8592% depended on drift pass .    a signal from another processor [ pac ] arising for events with many tracks in methane volumes included as anti - coincidence to the basic trigger [ cc3 ] . usage of this processor allowed skipping directly in the trigger a part of events with many tracks , thus we reduced the time of data read - out from tdc and increased an alive operating time of setup up to 85% . signals of counters of the active protection cck from charge cosmic radiation were also included in anti - coincidences to basic trigger [ cc2 ] . the description of active protection is resulted in part 5 .    at an initial stage we selected events with 2 electrons , directed in the different sides . the rate of triggers was 45/s . as tpc detector consists of three separate not pumped out volumes operating at atmospheric pressure filling with working gases and maintenance of necessary purity of gases appeared difficult enough . gas supply system consists of 3 separate cycles ( see fig . [ pic_gaz ] ) :    1 .   methane cycle  for two large side volumes , in sum about 10 m@xmath20 , 2 . xenon cycle  for central target volume of 2 m@xmath20 , 3 . subsystem for the proportional counters of active protection from cosmic radiation . two methane volumes were filled by gradual replacing of air with methane , at the first stage by pumping - out with submission of methane , then by flowing up of necessary quantity of volumes . filling of central target part with mixture xe+ch@xmath19 was essentially more difficult problem , because xe ( even with the natural composition of isotopes ) is a very expensive gas so flowing up becomes impossible . the technique with the chemical reaction described in  @xcite was used . first , the central volume was filled by co@xmath21 . then co@xmath21 gas was absorbed in special reactor with alkali ( naoh ) with a compensation of the pressure decrease by ( xe+ch@xmath19 ) mixture . it was necessary to use the admixture of methane to increase the drift velocity as such in pure xe at our drift field v@xmath22  kv / cm was too small . the gas composition was controlled in all volumes by chromatograph axt-002 - 01 ( russian ) with accuracy @xmath23  2% .    after filling the volumes the gas purification is provided in the closed cycles . to have high efficiency for large ( l @xmath24 m ) drift distances deep purification of the working gases from electronegative impurities ( o@xmath21 and h@xmath21o vapor ) is necessary up to a level better than 10 ppm . reactors with ni / sio@xmath21 adsorbent used at the early stage posses fine purification properties , but should be used at high temperature @xmath25c . at such temperature the small part of methane dissociates and hydrogen appears in the detector volumes what could noticeably change the drift velocity . but the main lack of this adsorbent is that it contains the big admixture of @xmath26u and @xmath27th . so the isotopes @xmath28rn and @xmath29rn appear in the detector due to diffusion from a reactor . the decay of their daughter isotopes @xmath30bi and @xmath31tl could imitate the events of double beta decay . we investigated various getters from the point of view of their cleaning properties , capacity on o@xmath21 ( it was necessary to purify large volumes of the gas @xmath32  m@xmath20 ) and the diffusion of rn isotopes . a getter of the active metals alloy was developed specially for us at the institute of chemistry and technology of organic compound . two reactors for methane and xenon volumes operating at room temperature were prepared . using of new reactors allowed to reduce rn concentration more than in 50 times . calibrating measurements are described below in section 10 . gas electronegative impurities were controlled using amplitude measurements with the radioactive source @xmath33fe located in windows in lateral tpc wall on drift distances of 24 and 129 cm  @xcite , and by the trigger rate of setup during measurements . tpc is surrounded by the gas proportional counters of active protection from cosmic radiation on all sides . counters were made of aluminum pipe square cross - section @xmath34  @xmath35 with 50 mkm central wire . the length of counters varies from 1.5 up to 3.5 m. four in fold joins them in groups with one amplifier per fold . total active surface of 6 panels of counters is @xmath36  m@xmath37 . the counters are pumped and filled by ch@xmath19/ar mixture and works 3 months without purification , then gas changes . the counters efficiency checked with the aid of the radioactive source @xmath38cd was not worse than 92% . all 6 panels are included in the coincedence scheme cck ( see fig . [ pic_el ] ) . it gives out a veto signal veto on full time of drift in methane in case of signals from any two panels after cosmic particle passage through tpc . in addition , the separate segments of protection register @xmath39-mesons crossed xe - volume also giving a veto signal in the trigger and registering in tdc . such including gave the additional suppression of events concerned with cosmic particles . accidental coincidences are equal 4% and are controlled regularly . to know the drift velocity in methane volume is necessary to calculate the electron kinetic energy using a form of the trajectory . the knowledge of drift velocity in xe volume is necessary , because the time difference in ch@xmath19 and xe volumes defines the coordinate on drift distance and also for restoration of event vertex . measurements of the drift velocities and checks of its dependence on drift length are carried out using the cosmic trigger . we look for cosmic particles ( for example , coincidence of the counters pointed on fig . [ pic_setup ] ) crossing the volume of tpc . the delay time of each sensitive wire depends on drift distance at known cosmic trigger . the centers of the delay time distributions for all sensitive wires approximate with a straight line well . the inclination of a straight line gives drift velocity as the distances between wires are fixed . to test this measurement we also use the time distributions for the central wires in ch@xmath19 , because we know the drift distance at fixed trigger ( with not so good accuracy ) . the drift velocity in ch@xmath19 coincides for both cases and equals w@xmath40  cm / sec . it weakly depends on small admixture of nitrogen in volume ( accumulated from the atmosphere ) and was measured regularly during the runs . changes were not more than 10% . the drift velocity in xe is determined from time distribution for central wire in xe volume in case of the cosmic trigger and ranges within @xmath41  cm / sec with xe concentration . the difference in drift times in ch@xmath19 and xe volumes determines the absolute event coordinate along an electric drift field @xmath42 = ( t_{xe } - t_{ch_4})\cdot(w_{ch_4}\cdot{}w_{xe})/(w_{xe}-w_{ch_4})\ ] ] and is fixed on tdc . at the correct work of the detector the times ( coordinates ) should be distributed uniformly . distribution of this difference between drift times in ch@xmath19 and xe volumes is resulted on fig . [ pic_dmexe ] . such distributions are checked few times in a day and allow watching the tpc condition . one can see that small efficiency fall takes place along the drift . a few measurement series were made using different xe concentration with the natural isotope composition for the optimal quantity of xe determination . tpc showed good stability with 50% concentration of xe and small ( 4% ) admixture of ic@xmath19h@xmath43 in order to cancel the `` excitation '' of xe wires . on the other hand monte - carlo simulations of the tpc efficiency for @xmath12-decay showed the weak dependence on xe concentration . so our measurements were made with the following gas composition : xe  ( 50% ) , ch@xmath19  ( 46% ) , ic@xmath19h@xmath43  ( 4% ) . such concentration corresponds to 5 kg of xe in central part tpc . drift velocity in the central volume decreased up to @xmath44  cm / mks . so to determine the coordinate along the drift pass ( measuring of the drift time difference in ch@xmath19 and xe volumes ) special block tdc with total time of 160 mks was used .    during 1000 hours of the exposition broken into two intervals of 515 and 485 hours because of failure with a high - voltage cable ( 10  kv supply of transformers for the sources of the stabilized current for a magnet ) about @xmath45 triggers were recorded . additional expositions with a radioactive source @xmath46bi for determination of the energy resolution and checking of the efficiency on - line programs were also made . in other runs with enlarged rn isotopes concentration in the central volume of the detector ( reactors ni / sio@xmath21 used ) we checked up the quality of the new reactors and correctness of work both off - line and modeling programs ( see section 10 ) simultaneously . reduction of the rn concentration by a factor of 50 due to use of reactors on the basis of an alloy of metals resulted in the number of double events reducing only @xmath47 times . at this level it became clear that an additional source of events with two electrons was the cosmic particles crossing the dead zones in xe volume . active protection from charged cosmic was reinforced with few new fold cosmic counters . it decreased quantity of double events by @xmath48 times . the large statistics of events from lateral covers without a direction selection has been recorded for the determination of the efficiency of separate sensitive wires in ch@xmath19 and xe volumes depended on drift distances . the efficiency appeared to be high enough , except for nearest to mylar films and to lateral covers wires . these results were used in monte - carlo calculations . -event , title="fig : " ] + _ number of sensitive wire and drift time coordinates _ -event , title="fig : " ] + _ number of cathode strip and drift time coordinates _ events with electrons escaped from xe in the two different sides in the right and left volumes filled with ch@xmath19 are selected . such selection twice reduces efficiency of the detector to @xmath12-decay , but considerably simplifies procedure of events distinguishing as it is not required to disentangle electron tracks escaped in one side . projection of a track in number of a sensitive wire and drift time coordinates is sinusoid , while in number cathodic strip and drift time coordinates is cycloid as cathodic strips are located under @xmath14 angle to sensitive wires . the loop or narrow part of cycloid is directed in specific direction related to drift depended on a motion direction of the electron . so one can select the electrons escaped from xe volume or from lateral covers . events candidates in double - beta decay had two electrons from xe volumes .    the on - line program uses two simple algorithms for events distinguishing . first one for selection of single tracks is restriction on number of hits n in each of methane volumes @xmath49  n  @xmath50 on 56 sensitive wires . second one is fast algorithm for determination of an electrons direction . in a time window related the trigger it is determined as @xmath51 , where @xmath52 is average time , @xmath53 , @xmath54  maximal and minimal drift times for each of methane volumes . then for each methane volume in this time window the numbers of cathode strip hits above and below this average drift time are calculated . if the number of strip hits above @xmath52 exceeds number of strip hits below @xmath52 it is considered that the narrow part of cycloid ( loop ) is directed downwards to drift that corresponds to one electron direction and in the other case to the opposite direction . such simple algorithms have reduced the detector rate from 45 event / s after the electronic trigger up to 2.2 event / s . the factor of background suppression at this stage is equal to 20 . the example of @xmath12-event is presented on fig . [ pic_2b_ev ] .    after on - line program s procession the information is stored on a disk and the further off - line analysis of events is proceeded . further this rate has been additionally reduced by @xmath55% by removal of the events with an additional track in methane volume ( 3-rd track ) , by more careful fitting and by removal of events with large energy cosmic particles ( straight tracks ) . residuary events passes full fit in both number of a sensitive wire and drift time coordinates and in number of cathode strip and drift time additional coordinates . it should be noted that fitting programs used only tracks in methane volumes . in the basic coordinates a projection of electrons tracks in each methane volume is fitted by function @xmath56 with four parameters , where x is coordinate along which number of a sensitive wire varies and z is coordinate along drift . in this case @xmath57 is the angle between a momentum of the electron and a magnetic field , @xmath58 is larmor radius , a is the mean line of a sinusoid and @xmath59 is a phase of a sinusoid . it is supposed that the momentum and corresponding larmor radius of electron in methane volume changes weakly . losses in 60  cm of ch@xmath19 gas make tens kev , and momentum of electron from the double beta decay is in a range from 100  kev up to 2  mev . to take into account multiple scattering in methane two breaks on a trajectory are supposed during the fitting with larmor radius being the same but angles @xmath57 and @xmath59 changed . this is especially essential for low energies when scattering grows . with one break the number of fitting parameters grows up to 8 , with two breaks  up to 12 for 56 points of the data in each of methane volume . software package of minimization minuit and software package hbook from library of programs for processing of physical measurements from cern  @xcite were used in fit .    since cathode strips are directed relative to sensitive wires under @xmath14 the coordinate along which varies number of a cathodic strip is @xmath60 , where x is earlier coordinate along which number of a sensitive wire varied , and y is coordinate along this wire . taking into account that projections of a electron trajectory are sinusoids in coordinates z , x and y , x are displaced on an angle @xmath61 depending on a direction of electron movement , we obtain within a constant @xmath62 the equation of two cycloids with loops directed upwards or downwards depending on a sign prior @xmath63 . a constant responsible for shift of all cycloid along the @xmath64 axis is obtained at fitting events in coordinates number of a strip , coordinate along the drift . -event ]    from two hypotheses of a direction of electron the one with the best @xmath65 is chosen . on fig . [ pic_fit ] experimental points and fit for the candidate of @xmath12-decay event when two electrons escaped from xe to the left and right methane volumes are shown . below is presented only one of the cycloids corresponding to the correct for the candidate in @xmath12 event electron direction and on the right side the loop of a cycloid is directed upwards while at the left is directed downwards . projections of a trajectories of the each electron in number of a sensitive wire and drift time coordinates ( and number of a cathode strip and drift time coordinates ) are fitted separately . spatial precision along drift at sensitive wires and cathodic strips differs . the signal amplitude from a cathodic strip is 6 times less than that from a sensitive wire and on one cathode strip there can be added signals from several wires , this leaded to shift of the time .    as a result of work of the fitting program the quantity of events had decreased from @xmath66  events / s up to @xmath67  events / hour , i.e. the factor of suppression of a background at this stage ( it is mainly electrons from lateral covers ) is  @xmath23  200 . then events are reviewed visually analyzing information from both ch@xmath19 and xe volumes .    for determination efficiency of the fitting program the additional runs with the increased concentration of @xmath29rn in xe volume were used where events with two electrons with a vertex in xe arose from decay of @xmath30bi  daughter nucleus from @xmath29rn chain . all of the recorded data were viewed without work of fitting programs and after their work . efficiency of the fitting and direction definitions programs is  @xmath23  70% for events with two electrons ( see section 10 ) . the basic tool of the software for monte carlo simulation development is package geant ( versions 3.21 ) from cern library  @xcite .    on the basis of this package the program model of the detector was built . it includes the exact geometry of the detector , the full trigger and the emulation of the on - line events selection . simulation results in a file formated completely identically to the data received in real experiment . then this file is analyzed in the same way as the experimental data .    crossing sensitive wires and strips of the detector forms in perpendicular to drift plane a matrix of virtual cells , such that in each of them there is a part of one wire and one strip . according to the borders of these cells the track of the charged particle breaks into slices . for each of them the drift distance is recalculated to drift time taking in account diffusion , displacement in a magnetic field and efficiency of wires . received times are written in the corresponding channels of wires and strips . the real hardware restrictions are taken into account when final times are obtained and transferred to the further processing . direct viewing of the events generated in the simulation has not revealed appreciable differences from real ones .        ' '' '' xe , % & 20 & 40 & 50 & 60 +    ' '' ''    @xmath68 & 1.50 & 0.98 & 0.79 & 0.64 +    ' '' ''    @xmath69 & 0.30 & 0.39 & 0.395 & 0.384 +    ' '' '' @xmath70 & 12.3 & 10.8 & 9.9 & 9.3 +    ' '' ''    @xmath71 & 2.46 & 4.32 & 4.95 & 5.58 +    efficiency of the detector has been determined using simulation depending on xe concentration . product of efficiency and number of xe atoms is the value , which is used to obtain the half - life . it shows weak dependence from xe concentration . the results of modeling are listed in the table 1 . efficiency for @xmath7-decay at chosen 50% xe concentration @xmath72 , while for @xmath4-decay @xmath73 at work with both drift sections . small efficiency to @xmath7-decay is connected with too soft electrons spectrum and small probability for electrons to escape from xe volume and to give the trigger . a separate measurements were made to check existing programs modeling correctness and to determine efficiency of on - line programs . we used a radioactive source @xmath46bi , which conversion electrons have energy @xmath74 and @xmath75  mev and lay at energy region of the electrons of @xmath0xe double beta decay . the intensity of a radioactive source @xmath46bi was calibrated with the aid of nai(tl ) detector . the only change in selection of events during the work with a radioactive source consisted in changing selection by a direction to search events with electrons moving from a source through the detector . to reduce background it was used the fact that position of a radioactive source is known . events in the window on the strips corresponding to position of a source with a correct value of the time difference of signals in ch@xmath19 and xe volume ( ` a methane - xenon difference ' ) were chosen . measurements with the detector were carried out for two positions of the source : near drift ( drift distance of 24  cm ) and far drift ( distance of 129  cm ) . efficiency of registration differed  @xmath23  1% from the modelling data for near drift and  @xmath23  10% for far drift . such discrepancy can be explained by ` cancellation ' of a source by a background ( double tracks ) and for far drift the additional contribution is brought with distortions of a drift field at the large drift distances .    to calibrate the detector and to check all procedures of event selection and monte - carlo simulation of detector work the additional run  @xmath23  17 hours of live time was made . in this run xe cleaning was carried out by ni / sio@xmath21 reactor resulting in @xmath29rn diffusion from the reactor to tpc that produced two - electron events with a vertex in xe volume due to @xmath30bi decay  a daughter nucleus of @xmath29rn chain . quantity of rn atoms in xe was determined by measurement of the delayed @xmath76 coincidences from decay of a cascade @xmath30bi @xmath77 @xmath30po with characteristic time @xmath78  mks ( a half - life of @xmath30po ) . since sensitive wires worked in a proportional mode , @xmath79-particle with energy  @xmath23  8  mev was selected by the amplitude . we used special tdc , operating in the mode `` start '' with full time @xmath80 ms . tdc was started by electron escaped from xe volume in one side  @xmath81 and stopped by a signal with the large amplitude in xe  @xmath79 . distribution on the arrival time of a signal @xmath79 answered to a half - life @xmath30po . measurement of such coincidences after cleaning xe volume by ni / sio@xmath21 reactor during @xmath48 hours gave number of coincidences @xmath82/h and when used the reactor on the basis of active alloys @xmath83/h . measurements were carried out during @xmath84 and @xmath85  hours . reduction of the @xmath86 number with time corresponded to a half - life time of @xmath29rn @xmath87  days . the similar method was used to determine the rn concentration in ch@xmath19 volume . efficiency to @xmath88 coincidences was calculated by monte carlo method using full geometry and imitation of work of programs and was @xmath89% . taking into account the efficiency concentration of rn in xe volume was equal to @xmath90 working with getter ni / sio@xmath21 and @xmath91 using new getter from the active metal alloy . efficiency for @xmath12-events from @xmath30bi decay has been calculated too , assuming that after decay of daughters rn ions @xmath30bi @xcite was formed and fixed on charged films separating gas volumes . full enough cascade of decays @xmath30bi  @xcite has been used to simulate the events with two electrons that can mimicry the @xmath12-decay . these events appeared as @xmath81-decay of @xmath30bi with the subsequent emission of the second electron due to internal conversion or a birth of compton or mller electron inside xe - volume . both electrons escaped to ch@xmath19-volumes and gave the trigger .    in 17 hours exposition 109 @xmath12-events were found . having concentration of @xmath29rn determined from the delayed @xmath88 coincidences in view of registration efficiency of events from @xmath30bi one could expect 119 @xmath12-events . observed and expected number of events coincided sufficiently well , that gave some confidence of the correctness of all selection procedures . after off - line processing resulted files were viewed visually . selection was carried out by the following criteria . 1 .   events with a wrong direction of the electron ( from a lateral cover ) were thrown out . a great bulk made such events . the fitting program usually made mistakes due to presence of the additional hits ( sparks ) in wires and in strips , and due to the distortion of the tracks as a result of scattering and distortion of the electric drift field near the mylar films . strengthening of selection criteria led to the further loss of efficiency . events with additional tracks in any of methane volumes both in top and bottom were thrown out . events of ( 1 + 2 ) types were  @xmath23  36/h . events with two or more tracks in xe - volume and events when the track in xe was absent (  @xmath23  1/h ) were thrown out also . events with a seen track of a cosmic particle in any of volumes were thrown out (  @xmath23  3/h ) . double viewing of the whole material was made . the efficiency of viewing was close to 100% . after visual viewing there were 954 events . in addition 45 events have been removed from the two - electron events sample with small quantity of the worked wires . events with actuated wires in xe volume n@xmath92 out of total 17 sensitive wires were thrown out . in ch@xmath19 volumes actuation of n@xmath93 wires out of 30 nearest to xe was required . efficiency has not decreased practically , but events with severely different number of actuated wires in l and r volumes have been removed , that most likely was connected to accidental coincidence in drift of the left and right tracks .    909 events recorded during 1000 hours remained for the analysis . -events , title="fig : " ] + distribution along the drift of remaining 909 @xmath12-events showed substantial growth of the number of the events from small drift distances . this distribution differed from the same distributions after on - line ( see fig . [ pic_dmexe ] ) and off - line processions . similar concentration of events was observed for the events concerned with @xmath30bi decay ( see fig . [ pic_dmexe2 ] ) . we assume that these events can be concerned with an interaction of the cosmic particles in multiwire detectors frames and also with possible concentration growth of ions @xmath30bi near the detector as a result of diffusion rn from a frame s material . events from small drift distances with the corresponding corrections of efficiency have been thrown out . number of @xmath12-candidate events was reduced more than twice while efficiency decreased only by 10% . energy spectrum of the single electron for @xmath12-candidate events was more rigid than spectrum of monte - carlo events passed full processing . therefore , energy restrictions were chosen according to monte carlo distributions .    on fig . [ pic_e_mc ] distributions of energy for each electron and total energy of two electrons for @xmath4 and @xmath4 monte - carlo events are presented . energy cuts for each electron @xmath94  mev for @xmath7 decay ( see fig . [ pic_e_mc ] ) were introduced with small efficiency loss . assuming the same statistics to be obtained during an exposition with @xmath0xe and all 254 remaining events to be background we estimate sensitivity of the detector in 1000 measurement hours as @xmath95  years .    to estimate sensitivity for @xmath4 decay energy cuts @xmath96  mev for every electron energy and for the sum energy @xmath97  mev were used . such wide interval for the sum energy was used because the electron energy was measured in the methane volumes and losses in xe were not taken into account . further measurement of ionization losses in xe volume and corresponding narrowing of an interval for the two - electron sum energy is planned . with such cuts there were 152 events , that led to limit estimation @xmath98  years for 1000 measurement hours . in the near future the exposition of tpc detector with @xmath0xe will be started . sensitivity of setup can be increased using both the top and bottom drift sections . the results discussed were obtained with half of detector ( only the top drift section ) . measurement ionization losses will give additional criteria in suppression of background events . including the second half of detector , measurement of ionization losses in xe and increase of measurements time will allow to improve results . authors are grateful to m.v.danilov for constant interest , support and discussion of results and possible ways of their improvement . work is supported by the grant of the russian basic research fond 02 - 02 - 16481 . 99 s.r.elliot , p.vogel , ann . nucl . part . sci . * 52 * , 115  ( 2002 ) . s.m.bilenky , hep - ph/0403245 . h.v.klapdor-kleingrothaus , a.dietz , h.l.harley , i.v.krivosheina , mod . a16 * , 2409  ( 2001 ) . h.v.klapdor-kleingrothaus , m.hirsh , z.phys . * a359 * , 361  ( 1997 ) . m.danilov , r.de voe , a.g.dolgolenko _ et al . b480 * , 12  ( 2000 ) . s.m.bilenky , s.t.petkov , hep - ph/0405237 . c.e.aalseth ( igex coll . ) , nucl . ( proc . suppl . ) * b87 * , 236  ( 1999 ) . m.bakalyarov ( heidelberg - moscow coll . ) , phys . rev . * d55 * , 54  ( 1997 ) . v.artemiev , e.brakchman , m.ivanovsky _ et al . _ , phys . lett . * b280 * , 159  ( 1992 ) . c.vuilleumier , j.busto , j.farine _ et al . _ , phys . rev . * d48 * , 1009  ( 1993 ) . r.luescher , j.farine , f.boehm _ et al . _ , phys . lett . * b434 * , 407  ( 1998 ) . r.bernabei ( dama coll . ) , nucl . ( proc . suppl . ) * b110 * , 88  ( 2001 ) . yu.gavriljuk , a.gangapshov , v.kuzminov _ et al . _ , yadernaya fizika , intern . conf . on non - accelerator new physics , ( dubna , 2003 ) . v.a.artemiev , e.v.brachman , m.a.ivanovsky _ et al . _ , yadernaya fizika v.63 , 1312  ( 2000 ) . v.a.artemiev , e.v.brachman , m.a.ivanovsky _ et al . _ , nim a v.477 , 414  ( 2002 ) . v.a.artemiev , e.v.brachman , m.a.ivanovsky _ et al . _ , nim a v.303 , 309  ( 1991 ) . cern program library ( http://cernlib.web.cern.ch/cernlib ) m.k.moe and d.d.lowenthal , phys . rev . * c22 * , 2186  ( 1980 ) . lbnl isotopes project ( http://ie.lbl.gov/databases/databases.html )
Abstract:

Answer: detector devis is tpc in the magnetic field . it is dedicated to the investigation of the double - beta decay of xe . setup sensitivity was estimated in the series of measurements with xe with natural isotopes composition . detector allows measuring two - neutrino double - beta decay of @xmath0xe with half - life less than @xmath1 years . institute for theoretical and experimental physics