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 phase - space formulation of quantum mechanics has received a renewed attention in the last decades @xcite because it allows to employ notions and tools of classical physics in the quantum realm . the central object in this approach is _ wigner s function _ @xmath0 which is merely a phase - space @xmath1 representation of the density matrix @xmath2 . in a closed system , the dynamical equation for the wigner function is equivalent to _ schrdinger s ( or von neumann ) equation _ while , in the @xmath3 limit , it recovers _ liouville s equation _ for the phase - space distribution in classical mechanics . another remarkable property of the wignerian representation is the _ average property_. the expectation value of a quantum operator @xmath4 is obtained from the corresponding classical variable @xmath5 ( related with @xmath4 by weyl s rule ( * ? ? ? * ch .  1 ) ) by a prescription analogous to the classical one , namely @xmath6 thus , in spite of not being necessarily positive , the wigner function can be considered as a sort of quantum - mechanical `` distribution '' . besides , the marginal distributions of @xmath7 give the true quantum probabilities of @xmath8 , as @xmath9 , or of @xmath10 , as @xmath11 . therefore , the wigner formalism provides a natural quantum - classical connection . this has been specially exploited in the search for quantum analogues of various classical effects ( e.g. , chaos ) or in the problem of the fuzzy borders between the quantum and classical worlds @xcite .    in contact with the surrounding medium , the ( open ) system experiences dissipation , fluctuations , and decoherence @xcite . typically one is interested in a system consisting of a few relevant degrees of freedom coupled to its environment , or bath , which has a very large number of degrees of freedom ( e.g. , phonons , photons , nuclear spins , etc . ) . the system is not necessarily microscopic , but it can be a mesoscopic system described by some collective variables which under appropriate conditions can behave quantum mechanically @xcite . because of the genericity of these situations and the fundamental issues raised ( e.g. , approach to thermal equilibrium , dissipation in tunnelling and coherence , measurement problem ) , the study of _ quantum dissipative systems _ is of interest in several areas of physics and chemistry @xcite . the dynamics of an open system can in many cases be formulated in terms of a _ quantum master equation _ for the ( reduced ) density matrix . in the wigner representation this can be compactly written as @xmath12 , with @xmath13 a certain evolution operator . taking the classical limit , the master equation for a particle of mass @xmath14 in a potential @xmath15 goes over the klein  kramers equation @xcite @xmath16 \wf % % % % % % % % % % % % % % % % % % % \;,\ ] ] which is simply a _ planck equation _ in phase space . the first two terms ( poisson bracket ) generate the classical reversible evolution ( liouville equation ) . the last term ( `` collision operator '' ) accounts for irreversible effects due to the coupling to the environment ( dissipation and fluctuations ) . the damping parameter @xmath17 measures the strengh of the coupling to the bath , which is at temperature @xmath18 . it is worth recalling that phase - space dynamics also includes problems reducible to `` mechanical '' analogues : josephson junctions , certain electrical circuits , chemical reactions , etc . a suitable non - perturbative technique to solve classical fokker  planck equations of systems with few degrees of freedom is the _ continued - fraction method _ @xcite . this is a special case of the expansion into complete sets ( grad s ) method to solve kinetic equations in statistical mechanics @xcite . in brief , the distribution is expanded into a series of orthogonal polynomials and , as a result , the kinetic equation is transformed into an _ infinite set of coupled equations _ for the expansion coefficients @xmath19 . approximate solutions are obtained by truncating the hierarchy of equations at various levels . to get manageable expressions , however , the truncation needs to be performed at a low level . in the continued - fraction method , instead of truncating directly , one seeks for a basis in which the range of index coupling is short ( ideally , the equation for @xmath19 involves @xmath20 , @xmath19 , and @xmath21 ) . then , the ( differential ) recurrence relations among the @xmath19 can be solved by iterating a simple algorithm , the structure of which is like that of a continued fraction ( [ app : rr - cf ] ) . for the klein  kramers equation  ( [ kk ] ) the tridiagonal chain of coupled equations for the @xmath19 is called the _ brinkman hierarchy _ it has been solved by continued - fraction methods to study classical particles subjected to dissipation and fluctuations @xcite . the method has been also extended to rotational brownian motion problems involving classical spins and dipoles ( vd . @xcite and references therein ) .    to deal with quantum dissipative systems , however , is a more delicate task @xcite . to begin with , phenomenological quantisation of dissipative systems poses fundamental problems ( e.g. , with the uncertainty and superposition principles ) . a rigorous route is to model the bath in a simple way , quantise it together with the system , and eventually trace over the environmental variables . however , the resulting theoretical frameworks are technically involved . _ quantum master equations _ , except for simple cases , are difficult to solve . an alternative is provided by _ quantum state diffusion _ methods ( vd . @xcite ) , where stochastic evolution equations for state vectors in hilbert space are introduced as computational tools . nevertheless , its implementation seems to be restricted to systems with discrete spectrum ( e.g. , oscillators , 2-state systems , etc . ) . similarly , _ quantum langevin equations _ @xcite are of limited use beyond nearly harmonic systems . for an arbitrary system , there exist exact _ path - integral _ expressions for the evolution of the density matrix ( involving the so - called _ influence functional _ that incorporates environmental effects ) @xcite . however , those expressions are difficult to evaluate , even numerically , because the propagating function is highly oscillatory , rendering numerical methods unstable at long times . finally , _ quantum monte carlo _ simulations can in principle be used . nevertheless , in spite of ongoing progress ( vd . @xcite ) , they are computationally complex and suffer from the ( dynamical ) sign problem .    this situation strongly motivates the development of _ alternative methods _ for quantum dissipative systems . inspired on their suitability for classical systems , continued - fraction techniques have been developed for a number of problems . shibata and co - workers @xcite applied them to solve a c - number quantum fokker  planck equation for a spin in a dissipative environment . vogel and risken @xcite employed continued - fraction methods to solve master equations in quantum nonlinear optics ( vd.also @xcite ) . recently @xcite this method has been extended to study genuine phase - space problems within the wigner formalism . a generalisation of the brinkman hierarchy for quantum master equations of the caldeira  leggett type was presented . it was shown that the continued - fraction method for the classical problem can , in principle , be adapted to solve this hierarchy , yielding a promising technique to study several classes of nonlinear quantum systems subjected to environmental effects .    in this article we first present the details necessary for the derivation of the quantum hierarchies and discuss their solution by continued fractions . the approach is then implemented for the problem of quantum brownian motion in periodic potentials ( a demanding problem with ( partly ) continuous spectrum ) . both the cosine and simple ratchet potentials ( lacking inversion symmetry ) are considered . for the former a number of previous results are recovered ( so testing the method ) and extended , while the physical interpretations are revisited within the wigner formalism . for particles in ratchet potentials we study the effects of finite damping ( inertia ) on the rectified velocities . the phenomenology is interpreted in terms of the interplay of thermal hopping , overbarrier wave reflection , and tunnelling . results for non - equilibrium dynamics under oscillating forcing are also discussed , which show the combined effect of quantum phenomena , thermal activation , and dissipation in a nonlinear system . a number of technical issues are consigned to the appendices . we start from the following generic form for the quantum master equation in the wigner representation @xcite : @xmath22 \wf % % % % % % % % % % % % % % % % % % \;. \;\ ] ] the first three terms , identifying @xmath23 , give the classical klein  kramers equation  ( [ kk ] ) . the mixed diffusion term @xmath24 is heuristically related with the colour of the quantum noise @xcite . the series of derivatives @xmath25 ( wigner  moyal term ) gives the quantum contribution to the unitary evolution of the closed system . conditions under which this series can be truncated are sometimes discussed , to avoid dealing with an infinite - order partial differential equation . however , we shall fully keep this term because it can be specially important in nonlinear systems @xcite . conditions of validity for the master - equation description are discussed in @xcite . equations of the type of ( [ wkk : dxp ] ) can be derived modelling the particle surroundings as a _ bath of oscillators _ representing the normal modes of the environment . however , a number of assumptions are typically involved like semiclassical or high - temperature bath ( rendering the kernels of the influence functional local in time ) , weak system - bath coupling ( born  markov approximations ) , etc . then , one would heuristically expect equation ( [ wkk : dxp ] ) being valid for small enough @xmath26 . notwithstanding this , the _ structure _ of the equation seems to be quite generic . thus , a quantum master equation recently derived for strong coupling @xcite , and valid for all @xmath26 , involves terms similar to those in ( [ wkk : dxp ] ) [ with @xmath8-dependent coefficients and renormalised @xmath15 ] . besides , the phase - space representation of the celebrated lindblad master equation is obtained by simply adding to ( [ wkk : dxp ] ) terms of the form @xmath28 and @xmath29 @xcite . most of these terms can be readily incorporated in the treatment below ( indeed , the term @xmath24 , absent in the original caldeira  leggett equation , is included to illustrate this ) . the same , in principle , would apply to possible extensions of the quantum master equation  ( [ wkk : dxp ] ) . first it is convenient to introduce appropriate scaled units . these are based on a characteristic length @xmath30 ( e.g. , the distance between the minima in a double - well oscillator , the period in a periodic potential ) and a characteristic energy @xmath31 ( e.g. , a barrier height ) . then one scales frequencies by @xmath32 , time by @xmath33 , forces by @xmath34 , the action by @xmath35 , etc . on the other hand , @xmath36 ( @xmath37 in a high-@xmath18 bath ) is handled as an effective temperature @xmath38 and we introduce some convenient thermal rescalings via @xmath39 . thus , in the equations @xmath10 is scaled by @xmath40 , the potential and @xmath41 appear divided by @xmath42 , @xmath17 enters in the combination @xmath43 , time is multiplied by @xmath44 , while the thermal de broglie wave length @xmath45 enters divided by @xmath30 . omitting any marks for the scaled quantities ( but keeping in mind specially the thermal rescaling of @xmath46 , @xmath15 and @xmath10 ) , the master equation is simply written as @xmath47 \wf % % % % % % % % % % % % % % % % % % % \ ; , \quad\ ] ] with the coefficient in the quantum sum given by @xmath48 planck s constant @xmath49 is introduced in terms of the characteristic action @xmath35 via the quantum parameter @xmath50 ( denoted @xmath51 in @xcite ) @xmath52 the second relation gives @xmath53 in terms of the kondo parameter @xmath54 , related by the third one with our @xmath50 . the classical limit is approached as @xmath55 , i.e. , letting @xmath56 . let us now introduce a splitting of the evolution operator @xmath13 that will facilitate the calculation of the matrix elements required when expanding @xmath57 into complete sets . on inspecting equations  ( [ wkk : scaled ] ) and  ( [ kappa : s ] ) , we see that extending the quantum sum down to @xmath58 , we obtain the part @xmath59 of the classical term . thus , we can decompose @xmath13 into the irreversible , kinetic , and potential parts : @xmath60 % % % % % % % % % % % % % % % % % % \lkin & = & -\big(\p-\bdxp\pp\big)\ , \px \\[0.5ex ] % % % % % % % % % % % % % % % % % % \lv & = & \sumiqo \qcoef^{(\iq ) } \ , \bv^{(2\iq+1 ) } \ , \pp^{(2\iq+1 ) } \end{array } \right . % % % % % % % % % % % % % % % % % % \;.\ ] ] this natural splitting has the advantage of dealing with @xmath15 ( the most problem dependent part ) in a unified way . besides , we have grouped the term @xmath61 with the kinetic term @xmath62 because they are structurally similar . in the expansion into complete sets approach to solve kinetic equations , the distribution is expanded in an orthonormal basis @xmath63 and the coupled equations for the coefficients @xmath64 derived @xcite . for the klein  kramers equation  ( [ kk ] ) one uses hermite functions of @xmath10 as basis ( * ? ? ? * sec  4.4 ) @xmath65^{-1/2 } % % % % % % % % % % % % % % % % % % \;.\ ] ] one starts with the expansion in the momentum because , while the @xmath10 dependence of the hamiltonian is fully specified ( kinetic energy @xmath66 ) , the potential part @xmath15 depends on the problem . thus , explicit manipulations can be done on the parts involving @xmath10 , which are valid for any system . on the other hand , the hermite basis has a number of advantages @xcite ; not the least is the handling of the derivatives @xmath67 in the dynamical equation by means of the associated creation and annihilation operators @xmath68 and @xmath69 ( this will prove very convenient in the quantum case ) .    for the klein  kramers equation the resulting equations for the expansion coefficients @xmath64 are called the _ brinkman hierarchy _ @xcite . this plays an important rle in classical systems , both for analytical treatments  derivations of @xmath70 expansions , etc. and to develope efficient non - perturbative numerical methods . in fact , it has the structure of a @xmath71-term recurrence relation which , after expansion in the position basis , can be solved by continued fractions . in what follows we shall proceed in the quantum case following as closely as possible the steps of the classical limit . we start expanding the wigner function as @xmath72 where we have extracted the boltzmann - type factor @xmath73 . this involves the auxiliary parameter @xmath74 and potential @xmath75 , which is frequently chosen proportional to @xmath76 , i.e. , @xmath77 . the results should not depend on the @xmath78 or @xmath79 used , but these may ( i ) put parts of the evolution operator @xmath13 in self - adjoint form and/or ( ii ) improve the stability and convergence of the eventual numerical implementation .    from the orthonormality of the @xmath80 , the expansion `` coefficients '' can be written as @xmath81 . then , differentiating with respect to @xmath46 and using the evolution equation @xmath12 , we get the dynamical equations for @xmath64 @xmath82 to get the matrix elements @xmath83 , which are still operators on the @xmath8-dependence of @xmath84 , one needs the @xmath85 transform of @xmath13 . this can be split as @xmath86 , with @xmath87 , @xmath88 , and @xmath89 given by ( [ wkk : lir : lkin : lv ] ) . the transformation of these operators is done in [ app : lb ] , taking advantage of results for _ normal ordering _ of @xmath90 and @xmath91 ( which we use to express the @xmath67 and @xmath10 s in @xmath13 ) . then , with the @xmath92 , @xmath93 , and @xmath94 obtained we compute their matrix elements between the @xmath95 , getting the @xmath83 . the details of this part of the calculation are given in [ app : qnm ] . introducing the results obtained for @xmath83 into @xmath96 one gets the following _ quantum ( brinkman ) hierarchy _ @xmath97 } \ , \big [ \g_{\ip}^{\iq,- } \ , \bv^{(2\iq+1 ) } \big ] \ , \ec_{\ip-(2\iq+1 ) } + \sqrt{(\ip-1)\ip}\ ; \gammamm \ , \ec_{\ip-2 } \nonumber\\ % % % % % % % % % % % % % % % % % % & & { } + \sqrt{\ip}\ ; \dm \ec_{\ip-1 } + \gammao_{\ip } \ , \ec_{\ip } + \sqrt{\ip+1}\ ; \dp \ec_{\ip+1 } \\ % % % % % % % % % % % % % % % % % % & & { } + \sqrt{(\ip+1)(\ip+2)}\ ; \gammapp \ , \ec_{\ip+2 } + \sum_{\iq=0}^{\infty } \big [ \g_{\ip+(2\iq+1)}^{\iq,+ } \ , \bv^{(2\iq+1 ) } \big ] \ , \ec_{\ip+(2\iq+1 ) } \nonumber % % % % % % % % % % % % % % % % % % \;.\end{aligned}\ ] ] the auxiliary ( damping ) parameters introduced read @xmath98 % [ ( \etap-\etapro)+2\ip(\etab-\etapro ) ] % % % % % % % % % % % % % % % % % % % \;,\ ] ] which involve the following @xmath78-related parameters ( note the sign exchange ) @xmath99 in equation ( [ qbh ] ) , the operators on the @xmath8 dependences are @xmath100 here @xmath101 , @xmath102 is the confluent hypergeometric ( kummer ) function ( * ? ? ? * ch . 13.6 ) , and @xmath103 is given by equation ( [ kappa : s ] ) . in is important to note that the @xmath104 act only on @xmath15 not on the @xmath105 ( [ app : lb ] ) . the quantum hierarchy ( [ qbh ] ) is equivalent to the caldeira  leggett master equation ( [ wkk : dxp ] ) . previously , several hierarchies had been derived @xcite , but the treatment was for closed ( hamiltonian ) systems @xcite or involved unsuitable bases like @xmath106 @xcite . in contrast , our hierarchy of _ hermitian moments _ constitutes a direct quantum generalisation of the brinkman hierarchy . to verify this , note that in the classical case only the @xmath58 terms survive in the sums of ( [ qbh ] ) ( [ app : qnm ] ) . then , setting the usual @xmath107 , equation  ( [ qbh ] ) reduces to @xmath108 % \;,\ ] ] which is indeed the celebrated hierarchy associated to the klein  kramers equation .    for potentials obeying @xmath109 , @xmath110 , the quantum hierarchy ( [ qbh ] ) is a recurrence relation with coupling to a finite number of terms . examples include the harmonic oscillator @xmath111 and the double - well ( duffing ) oscillator @xmath112 . finite - coupling recurrences are always reducible to @xmath71-term vector recurrences ( [ app : rr - cf ] ) and solvable by continued fractions , as in the classical case . however , for non - polynomial @xmath15 this approach is prevented by the infinite coupling range in the index @xmath113 , consequence of the infinite series of @xmath10-derivatives @xmath114 in the wigner  moyal term . in what follows we shall exploit the expansion in the position basis to show how this problem can be circumvented . recall that the coefficients @xmath64 are still functions of @xmath8 . inserting their expansion in an orthonormal basis @xmath115 into the quantum hierarchy ( [ dcdt : qtot : lbar ] ) , we get the generic form @xmath116 with matrix elements @xmath117 . to put ( [ qbh ] ) into the above form we just need to ( i ) substitute the operators ( and @xmath8-dependent coefficients ) by their matrix elements @xmath118 , ( ii ) attach position superscripts to the @xmath64 , and , finally , ( iii ) sum over them . on so doing we get the _ expanded quantum hierarchy _ @xmath119_{\ix\ixp}\ , \ec_{\ip-(2\iq+1)}^{\ixp } + \sqrt{(\ip-1)\ip } \sum_{\ixp } \gammamm \ , \delta_{\ix\ixp } \ , \ec_{\ip-2}^{\ixp } \nonumber\\ % % % % % % % % % % % % % % % % % % \fl & { } + \sqrt{\ip } \sum_{\ixp } \dm_{\ix\ixp}\ , \ec_{\ip-1}^{\ixp } + \sum_{\ixp } \gammao_{\ip}\ , \delta_{\ix\ixp } \ , \ec_{\ip}^{\ixp } + \sqrt{\ip+1 } \sum_{\ixp } \dp_{\ix\ixp}\ , \ec_{\ip+1}^{\ixp } \\ % % % % % % % % % % % % % % % % % % \fl & { } + \sqrt{(\ip+1)(\ip+2 ) } \sum_{\ixp } \gammapp \ , \delta_{\ix\ixp } \ , \ec_{\ip+2}^{\ixp } + \sumiqo \sum_{\ixp } \big [ \g_{\ip+(2\iq+1)}^{\iq,+ } \bv^{(2\iq+1 ) } \big]_{\ix\ixp}\ , \ec_{\ip+(2\iq+1)}^{\ixp } \nonumber % % % % % % % % % % % % % % % % % % \;.\end{aligned}\ ] ] considering that @xmath120 and @xmath77 , we see that to get all matrix elements involved we just need those of @xmath121 and of @xmath122 , namely @xmath123 because the @xmath104 only increase the order of the potential derivatives . the expanded hierarchy has the form of a system of ordinary differential equations , @xmath124 . if @xmath125 and @xmath126 are some large truncation indices for the bases @xmath127 and @xmath115 , we have a problem with @xmath128 equations . one may be tempted to try direct runge  kutta integration or numerical diagonalization of the associated @xmath129 matrix . however , the dimensions involved are typically very large , so that it is worth to pursue a continued - fraction treatment .    as mentioned before the @xmath10-recurrence has a finite coupling range for polynomial potentials . for other potentials , with an appropriate choice of the basis @xmath115 , the matrix elements ( [ matr - elem : der : v ] ) may vanish when the second index @xmath130 lays at a certain `` distance '' of the first @xmath131 . then , the expanded hierarchy would be a recurrence relation in the position indices with a finite coupling range , also amenable for a continued - fraction treatment . to proceed in these cases , we just need to transform the two - index extended hierarchies into ordinary recurrence problems ( i.e. , with one index ) . in order to apply the continued - fraction method when finite coupling range in one of the indices is attained , we need to convert first the two - index recurrences into ordinary recursions . in the resulting one - index recurrences the coefficients will be matrices acting on appropriate vectors formed with the @xmath132 . the generic expanded form ( [ dcdt : generic ] ) can be converted into an _ one - index _ recurrence relation in the momentum index introducing the following vectors and matrices @xmath133 where @xmath134 . that is , for a fixed @xmath113 one forms the column vector with all the @xmath132 with @xmath135 ( this index may take positive and negative values ) . similarly , for given @xmath136 we build the matrix @xmath137 with all indices @xmath138 .    before giving explicit expressions for the matrix coefficients we introduce the following notation : @xmath139 , @xmath140 , @xmath141 , and @xmath142 . then , the matrix quantum hierarchy ( [ dcdt : vectors : matrices : p ] ) can be written as @xmath143 % % % % % % % % % % % % % % % % % % & & { } + \mq_{\ip}^{- } \mc_{\ip-1 } + \mq_{\ip } \mc_{\ip } + \mq_{\ip}^{+ } \mc_{\ip+1 } \\[+0.5ex ] % % % % % % % % % % % % % % % % % % & & { } + \mq_{\ip}^{++ } \mc_{\ip+2 } + { \textstyle\sum}_{\iq\ge1 } \mq_{\ip}^{+(2\iq+1 ) } \mc_{\ip+(2\iq+1 ) } \nonumber % % % % % % % % % % % % % % % % % % \;.\end{aligned}\ ] ] we have accounted for @xmath144 , @xmath145 [ vd .  equation  ( [ qbh : pre - matrix : i ] ) ] and incorporated the @xmath58 terms into the central @xmath146 ones . the matrix coefficients explicitly read @xmath147_{\ix\ixp}\ , \\ % % % % % % % % % % % % % % % % % % \big ( \mq_{\ip}^{-- } \big)_{\ix\ixp } & = & -\sqrt{(\ip-1)\ip}\ ; \gammamm\ , \delta_{\ix\ixp } \\ % % % % % % % % % % % % % % % % % % \big ( \mq_{\ip}^{- } \big)_{\ix\ixp } & = & - \sqrt{\ip}\ , \dm_{\ix\ixp } - \big [ \g_{\ip}^{0,- } \bv ' \big]_{\ix\ixp}\ , \\ % % % % % % % % % % % % % % % % % % \big ( \mq_{\ip } \big)_{\ix\ixp } & = & - \gammao_{\ip}\,\delta_{\ix\ixp } \\ % % % % % % % % % % % % % % % % % % \big ( \mq_{\ip}^{+ } \big)_{\ix\ixp } & = & - \sqrt{\ip+1}\ , \dp_{\ix\ixp } - \big [ \g_{\ip+1}^{0,+ } \bv ' \big]_{\ix\ixp}\ , \\ % % % % % % % % % % % % % % % % % % \big ( \mq_{\ip}^{++ } \big)_{\ix\ixp } & = & -\sqrt{(\ip+1)(\ip+2)}\ ; \gammapp\ , \delta_{\ix\ixp } \\ % % % % % % % % % % % % % % % % % % \big ( \mq_{\ip}^{+(2\iq+1 ) } \big)_{\ix\ixp } & = & - \big [ \g_{\ip+(2\iq+1)}^{\iq,+ } \bv^{(2\iq+1 ) } \big]_{\ix\ixp}\ , % % % % % % % % % % % % % % % % % % \end{array } \right . % \;.\ ] ] for finite coupling range in @xmath113 ( e.g. , in polynomial potentials ) , the hierarchy ( [ qbh : matrix : p ] ) is a recurrence relation to which continued - fraction methods can now be applied . this amounts to convert the underlying large @xmath129 problem into @xmath125 problems with associated @xmath148 matrices [ or @xmath149 . as the recursion `` coefficients '' are matrices , one would need _ matrix continued - fraction _ methods ( [ app : rr - cf ] ) .      when the @xmath10-recurrence has an infinite coupling range ( in periodic potentials , the morse potential , etc . ) , we can convert the @xmath150-index hierarchy ( [ dcdt : generic ] ) into matrix form , but introducing recurrences in the position index @xmath131 @xmath151 that is , @xmath152 , with @xmath153 . thus , for a fixed @xmath131 one constructs the column vector with the @xmath132 for all @xmath154 , and for given @xmath155 , the matrix @xmath137 with all @xmath156 . comparing the matrix in ( [ dcdt : vectors : matrices : x ] ) with equation  ( [ q : matr - elem : p : i ] ) for @xmath157 , we can construct @xmath158 diagonal by diagonal @xmath159_{\ix\ixp}\ , \\ % % % % % % % % % % % % % % % % % % \big ( \mq_{\ix\ixp } \big)_{\ip,\ip-2 } & = & -\sqrt{(\ip-1)\ip}\ ; \gammamm\ , \delta_{\ix\ixp } \\ % % % % % % % % % % % % % % % % % % \big ( \mq_{\ix\ixp } \big)_{\ip,\ip-1 } & = & - \sqrt{\ip}\ , \dm_{\ix\ixp } - \big [ \g_{\ip}^{0,- } \bv ' \big]_{\ix\ixp}\ , \\ % % % % % % % % % % % % % % % % % % \big ( \mq_{\ix\ixp } \big)_{\ip,\ip } & = & - \gammao_{\ip}\,\delta_{\ix\ixp } \\ % % % % % % % % % % % % % % % % % % \big ( \mq_{\ix\ixp } \big)_{\ip,\ip+1 } & = & - \sqrt{\ip+1}\ , \dp_{\ix\ixp } - \big [ \g_{\ip+1}^{0,+ } \bv ' \big]_{\ix\ixp}\ , \\ % % % % % % % % % % % % % % % % % % \big ( \mq_{\ix\ixp } \big)_{\ip,\ip+2 } & = & -\sqrt{(\ip+1)(\ip+2)}\ ; \gammapp\ , \delta_{\ix\ixp } \\ % % % % % % % % % % % % % % % % % % \big ( \mq_{\ix\ixp } \big)_{\ip,\ip+(2\iq+1 ) } & = & - \big [ \g_{\ip+(2\iq+1)}^{\iq,+ } \bv^{(2\iq+1 ) } \big]_{\ix\ixp}\ , % % % % % % % % % % % % % % % % % % \end{array } \right . % \;.\ ] ] all other diagonals being zero @xmath160 , @xmath145 .    to recognise better the structure of these matrices , we decompose them into the `` free '' and potential contributions @xmath161 . for @xmath162 , including the kinetic and irreversible terms , we have the pentadiagonal structure ( tridiagonal for @xmath107 ) @xmath163 the part due to @xmath15 , on the other hand , has the following alternate dense structure , reflecting the odd powers of @xmath67 in the wigner  moyal term : @xmath164_{\ix\ixp } \!\!\!\!&\!\!\!\ ! 0 \!\!\!\!&\!\!\!\ ! [ \g_{3}^{1,+ } \bv^{(3)}]_{\ix\ixp } \!\!\!\!&\!\!\!\ ! 0 \!\!\!\!&\!\!\!\ ! [ \g_{5}^{2,+ } \bv^{(5)}]_{\ix\ixp } \!\!\!\!&\!\ ! \ddots \\[0.ex ] % % % % % % % % % % % % % % % % % % line-1 [ \g_{1}^{0,- } \bv']_{\ix\ixp } \!\!\!\!&\!\!\!\ ! 0 \!\!\!\!&\!\!\!\ ! [ \g_{2}^{0,+ } \bv']_{\ix\ixp } \!\!\!\!&\!\!\!\ ! 0 \!\!\!\!&\!\!\!\ ! [ \g_{4}^{1,+ } \bv^{(3)}]_{\ix\ixp } \!\!\!\!&\!\!\!\ ! 0 \!\!\!\!&\!\ ! \ddots \\[0.ex ] % % % % % % % % % % % % % % % % % % line-2 0 \!\!\!\!&\!\!\!\ ! [ \g_{2}^{0,- } \bv']_{\ix\ixp } \!\!\!\!&\!\!\!\ ! 0 \!\!\!\!&\!\!\!\ ! [ \g_{3}^{0,+ } \bv']_{\ix\ixp } \!\!\!\!&\!\!\!\ ! 0 \!\!\!\!&\!\!\!\ ! [ \g_{5}^{1,+ } \bv^{(3)}]_{\ix\ixp } \!\!\!\!&\!\ ! \ddots \\[0.ex ] % % % % % % % % % % % % % % % % % % line-3 [ \g_{3}^{1,- } \bv^{(3)}]_{\ix\ixp } \!\!\!\!&\!\!\!\ ! 0 \!\!\!\!&\!\!\!\ ! [ \g_{3}^{0,- } \bv']_{\ix\ixp } \!\!\!\!&\!\!\!\ ! 0 \!\!\!\!&\!\!\!\ ! [ \g_{4}^{0,+ } \bv']_{\ix\ixp } \!\!\!\!&\!\!\!\ ! 0 \!\!\!\!&\!\ ! \ddots \\[0.ex ] % % % % % % % % % % % % % % % % % % line-4 0 \!\!\!\!&\!\!\!\ ! [ \g_{4}^{1,- } \bv^{(3)}]_{\ix\ixp } \!\!\!\!&\!\!\!\ ! 0 \!\!\!\!&\!\!\!\ ! [ \g_{4}^{0,- } \bv']_{\ix\ixp } \!\!\!\!&\!\!\!\ ! 0 \!\!\!\!&\!\!\!\ ! [ \g_{5}^{0,+ } \bv']_{\ix\ixp } \!\!\!\!&\!\ ! \ddots \\[0.ex ] % % % % % % % % % % % % % % % % % % line-5 [ \g_{5}^{2,- } \bv^{(5)}]_{\ix\ixp } \!\!\!\!&\!\!\!\ ! 0 \!\!\!\!&\!\!\!\ ! [ \g_{5}^{1,- } \bv^{(3)}]_{\ix\ixp } \!\!\!\!&\!\!\!\ ! 0 \!\!\!\!&\!\!\!\ ! [ \g_{5}^{0,- } \bv']_{\ix\ixp } \!\!\!\!&\!\!\!\ ! 0 \!\!\!\!&\!\ ! \ddots \\[0.ex ] % % % % % % % % % % % % % % % % % % line-6 \ddots \!\!\!\!&\!\!\!\ ! \ddots \!\!\!\!&\!\!\!\ ! \ddots \!\!\!\!&\!\!\!\ ! \ddots \!\!\!\!&\!\!\!\ ! \ddots \!\!\!\!&\!\!\!\ ! \ddots \!\!\!\!&\!\ ! \ddots \end{array } \right ) % % % % % % % % % % % % % % % % % % \ ] ] with an appropriate choice of the basis @xmath115 the matrix elements @xmath118 may couple only a finite number of basis functions . then , the recurrence  ( [ dcdt : vectors : matrices : x ] ) would be solvable by continued fractions when the @xmath113-recurrence is infinite and this approach is seemingly precluded . even for finite coupling in @xmath113 , the @xmath131-recurrence may be preferable in certain situations ( e.g. , at low @xmath18 ) , as indicated by risken for a cosine potential in the classical limit ( * ? ? ? 11.5.6 ) ) . in these cases , one has been able to reduce the underlying large @xmath129 problem to @xmath126 problems with associated @xmath165 matrices . once the @xmath166 or @xmath167 are solved , we can reconstruct the wigner function from its expansion coefficients @xmath132 in the double basis @xmath168 @xmath169 obtaining the _ full solution of the quantum master equation_. if preferred , one can switch to the familiar density matrix @xmath170 by the inverse of transformation ( [ wigner : def ] ) @xmath171 we can get any observable inserting the above @xmath57 in equation  ( [ average ] ) for the averages . nevertheless , common observables can many times be extracted _ directly _ from the @xmath132 . let us illustrate this with transport observables , which are characterised by the averages @xmath172 . for an arbitrary function of @xmath10 only , we get from ( [ w : expansion : xp ] ) @xmath173 where we have used @xmath174 , with @xmath175 [ equations  ( [ hermite ] ) and ( [ w : expansion ] ) ] . next we introduce the auxiliary integrals ( for their calculation vd.[app : auxint ] ) @xmath176 then , the moments ( case @xmath177 ) can be written in terms of the expansion coefficients as @xmath178 . taking into account that @xmath179 has the parity of @xmath113 and carrying this to @xmath180 , we get for the first two @xmath181 the first moment is the _ current _ and the second characterises _ fluctuations _ around it . the zeroth moment corresponds to the _ normalisation _ of @xmath7 , and imposses on the @xmath132 the condition @xmath182 . finally , the quantum probabilities of @xmath8 and @xmath10 , given by the marginal distributions @xmath183 and @xmath184 , can be expressed in a similar fashion @xmath185 henceforth we shall apply the method described to incorporate fluctuations and dissipation in the problem of _ quantum transport in periodic potentials_. the simplest model consists of a particle evolving in a cosine potential subjected to external forces ( tilted `` washboard '' potential ) , which is also a paradigm of _ quantum brownian motion _ others include potentials lacking inversion symmetry ( ratchets ) , which have been used to model rectification of current and _ directional motion _ in several systems @xcite . these problems are demanding because periodic potentials have ( partly ) continuous spectra @xcite , rendering inappropriate methods devised for systems with discrete levels . in this section we compute the matrix coefficients of our recurrences for arbitrary periodic @xmath15 . with them we can implement the continued - fraction method to solve the quantum master equation , which we first check by regaining the classical dissipative and quantum hamiltonian limits .      in periodic potentials the @xmath10-recurrence can not be used because of its infinite coupling range ( the derivatives of @xmath15 neither vanish , nor decrease , as they are essentially the potential itself ) . however , one can guess that the @xmath8-recurrence may have a coupling range related with the number of harmonics in @xmath15 . to exploit this we introduce plane waves as basis functions and the fourier expansion of the potential derivative : @xmath186 to preserve the periodicity of @xmath7 , one extracts the external force @xmath187 from @xmath15 including it in a generalised @xmath188 ; cf .  equation  ( [ dmdp ] ) . then , the auxiliary potential @xmath189 is set proportional to the periodic part @xmath77 .    to compute the matrix elements in @xmath158 , we only need to do the integrals ( [ matr - elem : der : v ] ) with plane waves for @xmath190 . using @xmath191 , the integral @xmath192 , and the expansion @xmath193 , we find @xmath194_{\ix\ixp}=\v_{\ix-\ixp}^{(\iq ) } % % % % % % % % % % % % % % % % % % \;.\ ] ] from these elements and @xmath195 we get @xmath196_{\ix\ixp } & = & 0 & \delta_{\ix\ixp } & - & \etapm^{2\iq+1 } \gs_{\ip}^{\iq}(-\etapro \qab^{2}\ldb^{2 } ) & \v_{\ix-\ixp } ' \\ % % % % % % % % % % % % % % % % % % \dpm_{\ix\ixp } & = & \big(\iu\,\dplmi\ix-\etapm\bf\big ) & \delta_{\ix\ixp } & - & \dplmi\epspot & \bv_{\ix-\ixp } ' % % % % % % % % % % % % % % % % % % \end{array } % % % % % % % % % % % % % % % % % % \ ] ] with @xmath197 , @xmath101 and [ cf . equation ( [ gamma : qcoeff ] ) ] @xmath198 that is , @xmath199 is the coefficient in ( [ gamma : qcoeff ] ) with @xmath200 . inserting equation  ( [ dmdpvper : matr - elem ] ) into the general matrices @xmath158 of section  [ sec : mqbh ] we explicitly get the matrix coefficients of our recurrences ( [ app : matrices ] ) . for the cosine potential @xmath201 the coupling range is @xmath202 , because @xmath203 for @xmath204 and we have a @xmath71-term recurrence @xmath205 . for a @xmath150-harmonic ratchet potential @xmath206 % % % % % % % % % % % % % % % % % % % \;,\ ] ] the range is @xmath150 ( @xmath207-term recurrence ) because @xmath203 for @xmath208 , and so on . as we restricted @xmath74 , we have @xmath209 , so that @xmath210 . thus , the argument of @xmath211 is positive and @xmath212 can act as a regularisation factor [ an advantage of allowing @xmath213 in ( [ w : expansion ] ) ] . this factor reduces the weight of the off - diagonal terms inside @xmath158 , enhancing the numerical stability when going into the deep quantum regime . then we typically use @xmath214@xmath215 , while for classical calculations @xmath107 performs better . concerning the auxiliary potential @xmath189 , for the problems considered below we use @xmath216 . convergence is achieved with @xmath217 hermite functions and @xmath218 plane waves [ then @xmath219 . finally , a word on the scaled quantities ( section  [ sec : wkk : scaled ] ) . to get period @xmath220 , we have set @xmath221 with @xmath222 the original period . the characteristic energy @xmath31 is given by the potential amplitude @xmath223 ; then @xmath224 is scaled by @xmath225 and the other characteristic quantities are : action , @xmath226 , frequency , @xmath227 , and force , @xmath228 .    ) for a particle in a tilted sinusoidal potential ( two periods are displayed along @xmath8 ) . here @xmath229 , and @xmath230 , while @xmath231 , @xmath232 , @xmath233 , and @xmath234 ( left to right , top to bottom ) . , title="fig : " ] ) for a particle in a tilted sinusoidal potential ( two periods are displayed along @xmath8 ) . here @xmath229 , and @xmath230 , while @xmath231 , @xmath232 , @xmath233 , and @xmath234 ( left to right , top to bottom ) . , title="fig : " ] ) for a particle in a tilted sinusoidal potential ( two periods are displayed along @xmath8 ) . here @xmath229 , and @xmath230 , while @xmath231 , @xmath232 , @xmath233 , and @xmath234 ( left to right , top to bottom ) . , title="fig : " ] ) for a particle in a tilted sinusoidal potential ( two periods are displayed along @xmath8 ) . here @xmath229 , and @xmath230 , while @xmath231 , @xmath232 , @xmath233 , and @xmath234 ( left to right , top to bottom ) . , title="fig : " ]      before going into the quantum regime , let us check that for small enough @xmath235 we recover the classical results . recall that in the deterministic limit a particle in a periodic potential has two critical forces @xmath236 ( retrapping force ) and @xmath237 ( force at which the barrier disappears ) . for @xmath238 the attractors of the dynamical system are solutions `` locked '' around the potential minima , whereas for @xmath239 the `` running '' solution is an attractor globally stable . in the range @xmath240 the system exhibits bistability . finally , @xmath236 decreases with @xmath17 [ for the cosine potential @xmath241 , whereas @xmath236 tends to @xmath237 as the damping increases , narrowing the bistability range . thermal fluctuations are incorporated by the klein  kramers equation , which for classical particles in periodic potentials was solved by continued fractions in @xcite ( using the @xmath10-recurrence ) . here we solve the caldeira  leggett quantum master equation ( [ wkk : dxp ] ) using a large @xmath50 . as shown before , when @xmath242 we have a @xmath71-term recurrence which , for the static response , can be written as @xmath243 with @xmath244 , @xmath245 , and @xmath246 . solving it with the continued - fraction method we reobtain risken s impressive graphs for the classical distribution @xmath57 ( * ? ? ? some of them are displayed in figure  [ fig : risken ] for weak damping and various external forces ( here @xmath247 and @xmath248 ) . at low @xmath249 the distribution , always periodic along @xmath8 , has maxima at the potential minima , while the profile in momenta is a maxwellian envelope @xmath250 . as @xmath249 is increased we see the evolution from these _ locked _ solutions ( @xmath251 ) to bistability between the locked and _ running _ solutions , and eventually purely running states @xmath252 $ ] are obtained at large forces . having checked the retrieval of the classical results , we now enter the quantum regime by decreasing @xmath253 . firstly , it will be useful to get a more quantitative feeling of the relation of our parameter @xmath50 and `` how quantal is the system '' . to this end we study the bands of the corresponding hamiltonian problem ( without force ) . inserting the _ bloch _ function @xmath254 in the schrdinger equation , we arrive at the following eigenvalue problem ( in scaled units ; section  [ sec : wkk : scaled ] ) @xmath255 u(\x ) = e\ , u(\x ) % % % % % % % % % % % % % % % % % % \ ; , \quad % % % % % % % % % % % % % % % % % % -\half\leq \q\leq\half % % % % % % % % % % % % % % % % % % \;.\ ] ] expanding @xmath256 in plane waves ( [ vper ] ) and truncating @xmath257 at some large @xmath126 one gets a @xmath258 matrix which is diagonalised numerically . various bands @xmath259 pertaining to problems discussed below are displayed in figure  [ fig : bands ] . we see that the rule of thumb is that @xmath260 is of the order of the number of bands laying below the barrier . finally , the quantum parameter used by kandemir @xcite is @xmath261 , while chen and lebowitz employed @xmath262 @xcite . besides , for a josephson junction problem we have @xmath263 , with @xmath264 and @xmath265 the josephson and coulomb energies .     ( left ) and @xmath150 ( right ) . the `` islands '' correspond to zones of negative @xmath7 . , title="fig : " ]    ( left ) and @xmath150 ( right ) . the `` islands '' correspond to zones of negative @xmath7 . , title="fig : " ]      let us continue for a while in the hamiltonian case to familiarise ourselves with the structure of the schrdinger eigenstates in a phase - space representation . kandemir @xcite obtained analytical approximations for the eigenfunctions of a particle in a cosine potential in terms of mathieu functions . they were represented through the wigner function associated to the density matrix of pure states @xmath266 . ( in @xcite various dirac deltas occurring were regularised with gaussians @xmath267 and in the plots @xmath268 was used ; we attain the same setting @xmath230 since our @xmath10 is thermally rescaled . )    using a very small damping ( @xmath269 ) in the caldeira  leggett master equation we reobtain kandemir s hamiltonian results for the ground state ( figure  [ fig : kandemir ] ) . we see that at low @xmath50 the state is extended or delocalised , becoming more localised at the potential wells when @xmath270 is increased . according to @xcite , the regions of negative @xmath7 are associated to brillouin zone boundaries . we have represented a larger range of @xmath10 to see how these @xmath271 islands are repeated in the @xmath10 direction ( @xmath7 can be as small as @xmath272 bordering the limits of our numerical accuracy ) .    to see how these structures arise , note that for small @xmath50 the potential can be treated perturbatively [ equation  ( [ schrodinger : bloch ] ) ] . then , the wave functions have the form @xmath273 , with @xmath274 . the first term is a pure plane wave , and the rest first order corrections due to the sinusoidal @xmath15 . performing the wigner transform ( [ wigner : def ] ) of @xmath275 , we get for the ground state @xmath276 ( cf . @xcite ) @xmath277 \nonumber\\ % % % % % % % % % % % % % % % % % % & - 2|\lambda|^{2 } \cos(2\x)\ , \delta(\p ) + 2\iu\lambda \sin(\x ) \big [ \delta(\p-\half ) + \delta(\p+\half ) \big ] % % % % % % % % % % % % % % % % % % \;.\end{aligned}\ ] ] the first line gives the momentum localisation at @xmath278 ( unperturbed ) and at @xmath279 . the first term in the second line arises from the interference of @xmath280 and @xmath281 , and is responsible for the weak modulation of the height of the @xmath278 bell along @xmath8 . this modulation is more intense for larger @xmath282 ( i.e. , larger @xmath50 ) . the last term accounts for the interference between the plane wave centred at @xmath278 with those at @xmath279 , and can produce the negative islands . therefore , the above simple functional form captures most features of the wigner funtions of figure  [ fig : kandemir ] . after checking the connection with the classical and hamiltonian quantum limits , we finally include dissipation and temperature in the quantum case . here examples of particles in cosine and ratchet potentials will be discussed . transport properties of weakly damped particles in a cosine potential were studied by chen and lebowitz @xcite . they started from the system - plus - bath density matrix , traced out the bath variables , and performed a perturbative calculation in the potential height ( followed by a resummation ) to get @xmath283 . for low forces a free - particle like behaviour @xmath284 was obtained . increasing @xmath249 , the wave - vector associated to @xmath10 approaches the first zone boundary ( figure  [ fig : bands ] ) . there , while landau  zener tunnelling can bring the particle to the next band , bragg scattering reduces the velocity ( @xmath285 ) . eventually , at larger forces , @xmath10 corresponding to states inside the next band become favoured , and the free - particle behaviour is progressively recovered .     vs.  @xmath249 for various dampings at @xmath286 and @xmath287 . right : marginal distribution of momenta @xmath288 for @xmath289 at @xmath290 , @xmath291 , @xmath292 , @xmath293 , and @xmath207 . the arrows mark the zone boundaries @xmath294 . bottom panels : two views of @xmath57 for @xmath289 and @xmath290 . , title="fig : " ]   vs.  @xmath249 for various dampings at @xmath286 and @xmath287 . right : marginal distribution of momenta @xmath288 for @xmath289 at @xmath290 , @xmath291 , @xmath292 , @xmath293 , and @xmath207 . the arrows mark the zone boundaries @xmath294 . bottom panels : two views of @xmath57 for @xmath289 and @xmath290 . , title="fig : " ]   vs.  @xmath249 for various dampings at @xmath286 and @xmath287 . right : marginal distribution of momenta @xmath288 for @xmath289 at @xmath290 , @xmath291 , @xmath292 , @xmath293 , and @xmath207 . the arrows mark the zone boundaries @xmath294 . bottom panels : two views of @xmath57 for @xmath289 and @xmath290 . , title="fig : " ]    vs.  @xmath249 for various dampings at @xmath286 and @xmath287 . right : marginal distribution of momenta @xmath288 for @xmath289 at @xmath290 , @xmath291 , @xmath292 , @xmath293 , and @xmath207 . the arrows mark the zone boundaries @xmath294 . bottom panels : two views of @xmath57 for @xmath289 and @xmath290 . , title="fig : " ]    since high @xmath295 approximations were involved in their calculations , we set @xmath23 and @xmath296 in the master equation  ( [ wkk : dxp ] ) . solving it with the continued - fraction method we reobtain the effect just described ( figure  [ fig : cheleb92 ] ; here we work in the regime @xmath297 ) . note that this quantum slowing effect is reduced as @xmath17 increases , since the coupling to the bath `` broadens '' the effective energy levels . this broadening makes less relevant the presence of the band gap ( `` bridges '' it ) and brings the curves progressively closer to the free - damped - particle behaviour @xmath284 .    with our method we get , in addition , the wigner function ( i.e. , the full density matrix ) . the distribution of velocities ( [ wp ] ) , @xmath184 , when the curves start to rise again ( @xmath298 ) , shows two peaks separated by a minimum corresponding to the wave vector of the first zone boundary . loosely speaking the system shows coexistence of _ two _ quantum `` running '' solutions . classical running solutions have a wiggling structure along @xmath8 , as the particle slows down near the potential maxima ( figure  [ fig : risken ] ) . the full @xmath57 shows that the quantum running solutions have a nearly straight structure ( substrate insensitive ) . we also recognise the familiar crescent islands of negative @xmath7 at large @xmath10 ( there @xmath299 ) . around @xmath300 a complex structure is developed along @xmath8 ( `` locked part '' ) , with the maxima deformed and the minima becoming slightly negative ( @xmath301 ) . these features may reflect some interference between the locked and the nearby running solution . however , simply adding a plane wave @xmath302 to the @xmath303 used to analyse the wigner functions of the hamiltonial problem [ equation  ( [ wigner : turkeys ] ) ] we have not been able to reproduce such structure . anyway , it is erased when integrating over @xmath8 to get @xmath304 . we now turn our attention to periodic potentials without spatial inversion symmetry . this feature combined with out - of - equilibrium conditions can give rise to directional motion with net unbiased driving ( _ ratchet effect _ ) . this directional motion has been invoked to explain the behaviour of molecular motors and , with the emergence of nanoscience , opens the way to build mesoscopic devices and engines based on it . the theoretical work has concentrated on the classical behaviour and in the high - damping regime @xcite . the few quantum studies also addressed strong system - bath coupling in two situations : semiclassical @xcite and a extreme quantum case disregarding thermal effects @xcite . an exception to the large @xmath17 studies is @xcite , but there the substrate potential ( and the force ) were treated perturbatively . here , we shall solve the quantum master equation ( [ wkk : dxp ] ) with the ( non - perturbative ) continued - fraction method for particles in ratchet potentials , taking into account finite damping ( @xmath305 non - negligible inertia ) . for the ratchet potential ( [ vrat ] ) , the fourier coefficients of @xmath306 are @xmath307 and @xmath308 ( recall the thermal scaling of @xmath309 ; section [ sec : wkk : scaled ] ) . we use @xmath310 , which smooths a small shoulder that this potential exhibits on the `` easy '' side ( figure  [ fig : bands ] ) . as discussed before , the range of index coupling is @xmath150 , which gives a @xmath207-term recurrence . for the stationary response it can be written as [ cf.equation  ( [ rr:3:cosine ] ) ] @xmath311 which can be folded onto a canonical @xmath71-term recurrence by introducing appropriate ( block ) vectors and matrices ( [ app : rr - cf ] ) . the potential considered is the minimal extension of the cosine potential lacking inversion symmetry , so that the average velocities can be different for positive and negative forces ( vd . figure [ fig : cl - ratchet ] ) . we excite the system with a square - wave force switching alternatively between @xmath312 , and compute the `` rectified '' current of particles @xmath313 here @xmath314 are the correponding stationary velocities since we consider adiabatic conditions ( @xmath315 ) . in what follows we briefly investigate the classical limit ( understudied for finite damping ) and then proceed to study quantum corrections . figure  [ fig : cl - ratchet ] displays @xmath316 as a function of @xmath18 , showing the appearance of non - zero rectified velocities ( ratchet effect ) . the rectification is optimum at some intermediate temperature ( at too low @xmath18 there is hardly a response , and at too high @xmath18 the potential is irrelevant , so its asymmetry plays no rle in either case ) . at @xmath317 the results coindice with those obtainable from the analytical @xmath283 for overdamped classical particles ( * ? ? ? 11.3 ) ( vd . also @xcite ) @xmath318 with @xmath319/t$ ] . inertia , however , broadens the @xmath316 curves and shifts the maxima to higher @xmath18 ( for the lowest @xmath229 the maximum moves to a slightly lower @xmath18 ; vd . this broadening is accompanied by a slower decrease of @xmath316 with @xmath18 . this is consistent with the asymptotic behaviour of @xmath316 when @xmath320 @xcite . for overdamped particles it goes as @xmath321 , while for weak damping it decreases only with @xmath322 ( @xmath323 ) . we have checked that our results recover these dependences ( curves not shown due to the smallness of the asymptotic @xmath324 ) .    ) with @xmath229 and @xmath230 . two forces are used @xmath325 ( left ) and @xmath326 ( right ) , to show the lack of symmetry in the response ( running part more developed to the easy side ) and the locked regions reflecting the potential profile . bottom left : rectified classical velocity vs.  temperature for various values of the damping and force @xmath327 . right : @xmath328 vs.  @xmath18 for the positive force . , title="fig : " ] ) with @xmath229 and @xmath230 . two forces are used @xmath325 ( left ) and @xmath326 ( right ) , to show the lack of symmetry in the response ( running part more developed to the easy side ) and the locked regions reflecting the potential profile . bottom left : rectified classical velocity vs.  temperature for various values of the damping and force @xmath327 . right : @xmath328 vs.  @xmath18 for the positive force . , title="fig : " ]    in our case @xmath329 ( the `` retrapping '' force ; section  [ sec : risken ] ) , so that at @xmath330 the attractors are locked solutions . then , the absolute @xmath328 vs.  @xmath18 curves start from @xmath331 at @xmath330 , depart from zero as @xmath18 is increased and evolve towards @xmath332 for @xmath320 ( free damped particle ) . in the overdamped case this evolution is nearly a step , whereas the slope of the initial raising decreases with decreasing @xmath17 ( if @xmath238 ) . to understand this , let us assimilate @xmath283 to the escape rate @xmath333 ( modified by some mean free path ) , so that @xmath334 . kramers theory @xcite shows that @xmath335 increases monotonously with @xmath17 . therefore , the @xmath328 vs.  @xmath18 curves for finite @xmath17 must lay ordered below the overdamped result . on the other hand , the maximum rectification occurs in the temperature region where @xmath328 transits from @xmath336 to @xmath249 . this intermediate region is narrower ( and the slopes steeper ) the higher the damping is ( figure  [ fig : cl - ratchet ] ) , while the curves start to raise from zero at a lower @xmath18 . therefore , increasing the damping the @xmath337 curves will be narrower and the maxima shift to lower @xmath18 , as observed . eventually , this also allows to explain the anomalous curve showing the return ( @xmath229 ) . here @xmath338 , so that at @xmath330 the system starts close to the zone of bistability . then the particle can perform incursions into the running solution , increasing the initial slope of @xmath339 but now when the damping is lowered , shifting the maximum to a lower @xmath18 .      after the mandatory exploration of the classical limit , we proceed to make the system quantal , by decreasing @xmath50 . figure  [ fig : q - ratchet ] shows the rectified current with @xmath340 ( with @xmath341 bands below the barrier ; cf . figure [ fig : bands ] ) for various @xmath17 , together with the reference classical curves . we have set @xmath342 is all cases , to have roughly the same amount of locked and running components in the solutions . at high temperatures ( @xmath343 ) the classical and quantum curves coincide , as expected . however , decreasing @xmath18 down to @xmath344 , the ratchet effect gets reduced or enhanced , with respect to the classical result , depending on the dissipation . finally , at low temperatures ( @xmath345 ) the rectification is always reduced . with @xmath346 ( using @xmath342 ) for @xmath340 and the reference classical curve ( @xmath347 ) . the insets of @xmath348 and @xmath215 show enlargements of the peaks ( together with results for @xmath349 and @xmath350 ) . bottom right panel : kramers classical rate vs.  @xmath17 at various forces ( from melnikov  meshkov formula @xcite ) . the crosses mark the points with @xmath351 . , title="fig : " ]   with @xmath346 ( using @xmath342 ) for @xmath340 and the reference classical curve ( @xmath347 ) . the insets of @xmath348 and @xmath215 show enlargements of the peaks ( together with results for @xmath349 and @xmath350 ) . bottom right panel : kramers classical rate vs.  @xmath17 at various forces ( from melnikov  meshkov formula @xcite ) . the crosses mark the points with @xmath351 . , title="fig : " ]   with @xmath346 ( using @xmath342 ) for @xmath340 and the reference classical curve ( @xmath347 ) . the insets of @xmath348 and @xmath215 show enlargements of the peaks ( together with results for @xmath349 and @xmath350 ) . bottom right panel : kramers classical rate vs.  @xmath17 at various forces ( from melnikov  meshkov formula @xcite ) . the crosses mark the points with @xmath351 . , title="fig : " ]   with @xmath346 ( using @xmath342 ) for @xmath340 and the reference classical curve ( @xmath347 ) . the insets of @xmath348 and @xmath215 show enlargements of the peaks ( together with results for @xmath349 and @xmath350 ) . bottom right panel : kramers classical rate vs.  @xmath17 at various forces ( from melnikov  meshkov formula @xcite ) . the crosses mark the points with @xmath351 . , title="fig : " ]    to understand qualitatively the underlying physics we turn to some semiclassical results for various related problems . first , let us consider the _ reflection / transmission coefficient _ for an asymmetric saw - tooth potential @xcite . the results depend on the energy of the incident quantum particle . the reflection coefficient for energies higher than the barrier , @xmath352 , and the underbarrier transmission coefficient , @xmath353 , read @xmath354 % % % % % % % % % % % % % % % % % % % \;,\ ] ] where @xmath355 , @xmath356 are the slopes at both sides of the cusp , @xmath357 the particle energy , and @xmath358}$ ] . suppose now that we have an asymmetric potential , with an easy direction to the right , and we deform it with forces @xmath312 , computing the corresponding reflection coefficients @xmath359 and @xmath360 . for energies higher than the barrier , the above @xmath352 gives @xmath361 ( [ app : semicl ] ) . that is , _ overbarrier wave reflection _ is less intense when the slope of the potential is smaller . thus overbarrier transmission ( as thermal hopping ) is favoured along the easy direction for energies above the barrier @xcite . in contrast , using @xmath353 for energies below the barrier , we find @xmath362 , so that @xmath363 . that is , the reflection is more intense in the easy direction or , equivalently , tunnelling events are favoured in the hard direction @xcite .    for a distribution of particle energies one would expect a competition between these two types of quantum effects . then , the enhancement of the rectified velocity under @xmath312 forcing for the @xmath348 curve follows from the increase of the thermal escape with @xmath17 in the intermediate - to - low damping range @xmath364 ( figure  [ fig : q - ratchet ] ) . at that damping more escape events occur and the particles launched over the barrier are subjected to the discussed phenomenon of wave reflection , which favours net motion along the easy side and hence increases @xmath316 . lowering the damping , on the other hand , less escape events are produced and wave - reflection becomes progressively dominated by tunnel events ( favouring the hard direction ) . at @xmath344 one would expect the tunnel to be thermally assisted through the higher bands , which are wider . a simple estimation of channel probabilities comes from multiplying a band thermal population ( using kramers rate ) by the probability of tunnel ( given by the semiclassical transmission coefficient @xmath353 ) . here the exponential reduction of @xmath333 with @xmath357 compites with the exponential increase of @xmath353 for the upper bands ( [ app : semicl ] ) . at @xmath344 and for @xmath340 we find that tunnel is indeed favoured through the upper bands .    to check that the behaviours found at the maxima are systematic , we increase @xmath50 to @xmath349 ( @xmath365 bands below the barrier ) and reduce it to @xmath366 ( @xmath291 bands ; here the temperatures around the peak are still inside the range of validity of the master equation , small @xmath26 ) . the results ( insets ) show clearly the tendency to amplify either the enhancement or the reduction when the system becomes more quantal . finally , at low temperatures ( @xmath367 ) thermal activation becomes ineficient for all dampings . then a reduction of the rectification is expected , due to dominance of below - barrier energies and hence of tunnel events ( through the lowest bands ) . note that the decrease of @xmath316 relative to the classical results is comparable in the curves shown . this suggests that we have similar crossover temperatures @xmath368 in this damping range . ( @xmath368 , below which quantum effects dominate , is useful as a measure of quantum corrections at a given @xmath18 ) . this quantity is known exactly in the problem of the escape from a `` quadratic - plus - cubic '' potential @xmath369 $ ] ( * ? ? ? 14 ) @xmath370 in our units @xmath371 and in the damping range studied the factor in brackets is close to one ( @xmath372 ) . this gives crossovers @xmath373 and hence comparable reductions of the ratchet effect , as observed . this reduction of the rectification is the precursor of the current reversals ( @xmath374 ) found in these systems when tunnelling completely dominates . although we find small @xmath374 in some cases , the results can not be fully trusted . for example , at @xmath348 , @xmath340 and @xmath375 we have @xmath376 , which is bordering the limits of validity of the quantum master equation employed here . we conclude with the study of the ( non - adiabatic ) dynamical response to oscillating forces @xmath377 . then , the coefficients of the expansion ( [ w : expansion : xp ] ) of the wigner function are periodic functions of @xmath46 , so they can be fourier expanded as follows : @xmath378 to lowest order in the probing field , the static part of the corresponding vectors [ equation  ( [ dcdt : vectors : matrices : x ] ) ] obeys a recurrence relation of the type ( vd . [ app : pert : ac ] ) @xmath379 the equations for the harmonics have the form ( @xmath380 is the identity matrix ) @xmath381 % % % % % % % % % % % % % % % % % % % \;,\ ] ] with the `` forcing '' @xmath382 involving the previous order results . specifically , @xmath383 , with @xmath384 the part of @xmath385 corresponding to @xmath386 , namely @xmath387 % % % % % % % % % % % % % % % % % % line-1 \sqrt{1}\,\etam\delta f \!&\ ! 0 \!&\ ! \sqrt{2}\,\etap\delta f \!&\ ! 0 \!&\ ! \ddots \\[0.ex ] % % % % % % % % % % % % % % % % % % line-2 0 \!&\ ! \sqrt{2}\,\etam\delta f \!&\ ! 0 \!&\ ! \sqrt{3}\,\etap\delta f \!&\ ! \ddots \\[0.ex ] % % % % % % % % % % % % % % % % % % line-3 0 \!&\ ! 0 \!&\ ! \sqrt{3}\,\etam\delta f \!&\ ! 0 \!&\ ! \ddots \\[0.ex ] % % % % % % % % % % % % % % % % % % line-6 ( 4 ) \ddots \!&\ ! \ddots \!&\ ! \ddots \!&\ ! \ddots \!&\ ! \ddots \end{array } \right ) % % % % % % % % % % % % % % % % % % \ ] ] [ vd . [ app : matrices ] for the general @xmath158 ; the @xmath384 actually entering in @xmath382 is divided by @xmath386 , cf . equations ( [ eom : generic ] ) and ( [ coefficients : a ] ) ] . equations ( [ 3rr : statics ] ) and ( [ 3rr : harmonics ] ) can be solved sequentially with the continued - fraction method . figure  [ fig : suscep - w ] shows the linear susceptibility ( @xmath388 ) computed in this way for a particle in a _ cosine _ potential as a function of @xmath389 . ( here we return to safe ground because @xmath390@xmath391 ; for related dielectric and kerr relaxation curves of a quantum rotator vd . @xcite and references therein . ) for the largest damping the results correspond to the classical limit . there the line - shape @xmath392 broadens and extends to frequencies lower than the oscillation frequency near the bottom of the potential wells ( @xmath393 in our units ) . for low damping , this classical _ spreading of oscillation frequencies _ has an important contribution from the amplitude dependence of the period of oscillation in anharmonic potentials @xcite . this non - dissipative contribution to @xmath392 is also known in classical spin problems , where the precession frequency depends on @xmath394 due to the magnetic anisotropy @xcite .     and @xmath395 with @xmath289 ( @xmath305 classical ) , @xmath396 and @xmath397 ( less to more peaked curves ) . left panel : real part . right panel : imaginary part ( dashed line , @xmath398 but halving @xmath399 ) . vertical lines : loci of the transition frequencies of the associated quadratic - plus - quartic anharmonic oscillator [ equation  ( [ delta - pert - fey ] ) ] . , title="fig : " ]   and @xmath395 with @xmath289 ( @xmath305 classical ) , @xmath396 and @xmath397 ( less to more peaked curves ) . left panel : real part . right panel : imaginary part ( dashed line , @xmath398 but halving @xmath399 ) . vertical lines : loci of the transition frequencies of the associated quadratic - plus - quartic anharmonic oscillator [ equation  ( [ delta - pert - fey ] ) ] . , title="fig : " ]    the quantum regime is approached by decreasing the kondo parameter @xmath400 [ equation  ( [ kondos ] ) ] . reducing @xmath17 with a fixed @xmath50 we find that the @xmath401 curves develop a multipeaked structure . these peaks could be vaguely seen as the mentioned nonlinear oscillations becoming quantised . to give content to this statement , we investigate perturbatively the effect of the anharmonic terms of a cosine potential ( those beyond @xmath402 ) on the energy levels of the harmonic oscillator part . since the main contribution should come from the quartic term , we consider the potential @xmath403 . first - order perturbation theory gives the eigenvalues ( * ? ? ? 6.3 ) @xmath404 expanding the cosine potential @xmath405 , we identify @xmath406 and @xmath407 ( @xmath408 , soft anharmonicity ) . then we can compute the level spacing @xmath409 % % % % % % % % % % % % % % % % % % % \;,\ ] ] which has acquired a dependence on the level index @xmath410 absent in the harmonic case . the associated transition frequencies @xmath411 are plotted in figure  [ fig : suscep - w ] [ scaled as @xmath412 . they agree well with the location of the main peaks of the quantum @xmath392 curves , supporting their relation with the nonlinearity of the potential . furthermore , the contribution from transitions between the higher levels should be reduced when they are less populated . decreasing @xmath18 , we indeed see that the corresponding peaks ( located at lower @xmath389 ) are substantially reduced . in the suggestive language of nonlinear oscillations we would say that lowering the temperature the large amplitude oscillations ( having longer period ) become less probable . the applicability of the quantum nonlinear oscillator picture is grounded on the flattness of the lowest bands at the @xmath50 considered , which are well approximated by constant levels ( vd.figure  [ fig : bands ] ) . we have computed the semiclassical width of the bands in a cosine potential in [ app : semicl ] . an upper bound for the width [ equation  ( [ wb : approx : top ] ) ] reads in dimensionless units @xmath413 $ ] , with @xmath414 the band mean energy relative to the potential amplitude . for moderate @xmath50 this bound gives in fact very flat bands , specially the deepest ones ( @xmath415 ) . this also implies that the width of the @xmath392 peaks has a negligible contribution from the intrinsic level width , being dominated by the damping and thermal broadening . traditionally , the continued - fraction method has been employed to solve fokker  planck equations for few - variable classical systems ( translational and rotational problems ) . when compared with numerical simulations of the langevin equation , the method has several _ shortcomings _ : ( i ) it is quite specific of the problem to be solved ( one needs to recalculate the recurrence coefficients for each potential ) , ( ii ) its convergence and stability depend on the parameters of the problem ( and can be poor in some ranges ) , and ( iii ) it does not return `` trajectories '' , which in the simulations provide helpful insight . nevertheless , when the method can be used its _ advantages _ are most valuable : ( i ) it is free from statistical errors , ( ii ) is essentially nonperturbative , ( iii ) specially apt to get stationary solutions ( static and dynamic ) , ( iv ) high efficiency ( allowing to explore parameter ranges out of the reach of the simulation ) , and ( v ) it gives the distribution @xmath7 , which certainly compensates for the lack of trajectories . for these reasons , and the lack of quantum langevin simulations , it was worth to develope continued - fraction methods for quantum dissipative systems . this had been done for problems of spins in a thermal bath and in quantum nonlinear optics . in this article we have discussed in detail the adaptation of the continued - fraction method to tackle quantum master equations in phase - space problems , taking advantage of the parallels with the classical case provided by the wigner formalism . we have seen that the quantum extension of the continued - fraction method is more problem dependent ( except for polynomial potentials ) , due to the necessity of using the @xmath8-recurrence of the quantum hierarchies and finding suitable bases . nevertheless , it inherits from the classical methods most of the abovementioned advantages , in particular , the obtaining of the wigner function , from which any observable can be computed . furthermore , the eigenvalue spectrum of the system , although helpful in understanding the physics , is not necessarily required . this is advantageous when dealing with non - bounded hamiltonians , continuous spectra , etc . finally , by changing the appropriate quantum parameter , the connection with the classical results is attained in a natural way . note that `` continued - fraction '' is a quite generic term ( like `` series expansion '' ) , which can be found in many different contexts , and in particular in works dealing with quantum dissipative systems ( see , e.g. , @xcite and references therein ) . however , in those cases one usually obtains some continued - fraction expression for a certain quantity ( e.g. , the linear susceptibility , memory functions , etc . ) , whereas we have obtained the _ complete solution of the quantum master equation _ by continued - fraction methods . the approach allows to study the interplay of quantum phenomena , nonlinearity , thermal fluctuations , and dissipation in several classes of systems . this is important , because methods optimised for one of the sides ( e.g. , nonlinearity or quantal behaviour ) tend to perform poorly for the others ( e.g. , fluctuations or dissipation ) . in addition , the visualization of @xmath57 plus the knowledge of the classical phase - space structure ( orbits , separatrices , attractors , etc . ) can provide valuable insight in difficult problems . finally , the main limitations of the approach are those that the starting master equations may have . we have considered the celebrated caldeira  leggett quantum master equation , but the method can be applied to a class of equations of this type , and to possible generalisations .    we have implemented the method in the canonical problem of quantum brownian motion in periodic potentials . we have considered both the cosine and ratchet potentials . for the former we have recovered and extended a number of classic results and interpreted them on the basis of the wigner formalism . results for the dynamics under oscillating forcing have also been obtained , illustrating the modification of the nonlinear effect of spreading of oscillation frequencies by quantum phenomena . for particles in ratchet potentials we have studied , under adiabatic conditions , the effects of finite damping on the rectified velocity . taking into account the competition of thermal hopping , overbarrier wave reflection and tunnel events , together with semiclassical analytical results , has allowed us to understand in great detail the physics of these systems . this work was supported by dges , project bfm2002 - 00113 , and dga , grant no .  b059/2003 . we acknowledge f.  falo , j.  j.  mazo , and a. p. piazzolla for useful discussions and valuable inspiration . here we briefly discuss the solution of recurrence relations by continued - fraction methods . we first consider @xmath71-term recurrences and then differential recurrences , vector - matrix cases , folding of larger - coupling recursions into @xmath71-term ones , and the problem of the intial conditions . ( for a brief description of the relation of continued fractions , series expansions , recursions , and orthogonal polynomials see @xcite . ) suppose we have the simplest case of a @xmath71-term recurrence relation of the form @xmath416 with @xmath417 and @xmath418 given quantities . to obtain the @xmath19 one inserts in ( [ rr ] ) the anstze @xcite ( upper and lower signs for @xmath419 and @xmath420 ) @xmath421 getting the following ladder coefficients @xmath422 and shifts @xmath423 @xmath424 if the recurrence is finite or the @xmath19 decrease quickly with @xmath425 ( e.g. , when they are coefficients of some judiciosly chosen expansion ) , we can truncate at some large @xmath426 , by setting @xmath427 and @xmath428 . next we _ iterate downwards _ in ( [ s : a ] ) getting @xmath422 and @xmath423 , down to @xmath429 . then , to get the @xmath19 from ( [ c : risken : com ] ) , we need @xmath430 , which is obtained from @xmath431 combination of equation  ( [ rr ] ) at @xmath429 and @xmath432 . thus , starting from @xmath430 , we _ iterate @xmath433 upwards _ , getting the solution of the recursion . for @xmath434 ( homogeneous recurrences ) , we have @xmath435 and the solution reads @xmath436 . note that @xmath422 is given in terms of @xmath437 in the denominator , which can in turn be written as a fraction with @xmath438 in the denominator , and so on . this has the structure of a _ continued fraction _ , naming the method @xmath439    the iterative solution can be cast into explicit form as a series of products of continued fractions @xcite . the above form , however , has the virtue of being easy to implement in a computer . convergence is checked increasing @xmath440 and repeating the procedure . solving equation ( [ rr ] ) in this way requires a computational effort of order @xmath441 , instead of the order @xmath442 of an arbitrary linear algebra problem . this reduction arises from the tridiagonal structure of the matrix associated to ( [ rr ] ) . ( for analytical inversion of tridiagonal matrices see @xcite and references therein . ) we can solve analogously a differential recurrence of the type @xmath443 laplace transformation converts this equation into [ @xmath444 @xmath445 % % % % % % % % % % % % % % % % % % % \;,\ ] ] where @xmath446 has been used . then , introducing @xmath447 and @xmath448 , the above equation has the structure of the ordinary recurrence ( [ rr ] ) . the quantities involved in the above recurrences need not to be scalars , but they can be @xmath449-vectors ( @xmath19 and @xmath418 ) and the coefficients @xmath417 are @xmath450-matrices . then we speak of _ matrix continued fractions_. the only change in the above solution is that the fraction bars then mean matrix inversion ( `` from the left '' , @xmath451 ) . this is important , because a recurrence relation involving more than three coefficients can be `` folded '' into a @xmath71-term recurrence by introducing vector and matrix quantities . let us show how to do so for the following @xmath207-term recursion @xmath452 defining the @xmath150vectors @xmath453 and the @xmath454matrices @xmath455 it can be easily seen that equation  ( [ rr:5 ] ) for @xmath456 and @xmath457 is equivalent to @xmath458 then , insertion of @xmath459 gives the matrix version of ( [ s : a ] ) . note that the @xmath19 and @xmath418 can already be vectors and the @xmath417 matrices . then , the above @xmath460 , @xmath461 , and @xmath462 are `` block '' vectors and matrices . finally , to iterate upwards @xmath459 , the `` initial condition '' @xmath463 is obtained from the matrix version of ( [ condini ] ) ( which we write @xmath464 ) . as in the absence of forcing ( @xmath465 ) the general solution involves an overall multiplicative constant , one adds to this system an extra equation to fix it ( e.g. , a normalisation contidion on the coefficients ; vd .  section  [ sec : obs ] ) . alternatively , we can fix arbitrarily one element of @xmath463 and rescale the solution at the end . in either case , one has an extra equation which added to @xmath466 yields the extended @xmath467 system @xmath468\hline\\[-2.ex ] % % % % % % % % % % % % % % % % % % \mbox{\,l.h.s.\ eq.\ , } \end{array } \right ) } \hspace{-0.5em } \begin{array}{c } \uparrow \\ { \scriptstyle ( j+1 ) } \\ \downarrow \end{array } % % % % % % % % % % % % % % % % % % \quad % % % % % % % % % % % % % % % % % % \mathbf{b}^{\ast } = \left ( \begin{array}{c } b_{1 } \\[-1.5ex ] \vdots \\[-1.5ex ] b_{j } % % % % % % % % % % % % % % % % % % \\[0.ex]\hline\\[-2.5ex ] % % % % % % % % % % % % % % % % % % \mathrm{\!\!r.h.s.\!\!\ ! } \end{array } \right ) % % % % % % % % % % % % % % % % % % \ ] ] this system can be solved by any method appropriate for cases with more equations than unknowns ( e.g. , using @xmath469- or svd - decomposition @xcite ) , getting @xmath463 . here we derive @xmath470 , the @xmath85 transform of the operator @xmath13 in the master equation with @xmath471 $ ] [ equation  ( [ w : expansion ] ) ] . we split @xmath470 as @xmath86 , with @xmath87 , @xmath88 , and @xmath89 given by ( [ wkk : lir : lkin : lv ] ) . for the `` product '' of operators we use that the transformation of the product is equal to the product of the transformations : @xmath472 which follows by sandwiching `` identities '' @xmath473 between the factors . besides , the multiplication by functions of @xmath1 commutes with @xmath85 , so that only the differential terms need to be transformed . from these two properties we see that _ only @xmath474 and @xmath475 need to be calculated _ ( and then raised to some power or multiplied ) . _ calculation of @xmath474_. in order to transform @xmath121 , we apply @xmath476 , to an arbitrary function @xmath477 . the @xmath10 dependent part of @xmath73 cancels out and we obtain @xmath478 which has the form of a `` displaced '' differential operator . _ calculation of @xmath475_. applying @xmath479 to a function @xmath480 , we similarly get @xmath475 : @xmath481 this @xmath475 can be expressed in terms of the following creation and annihilation operators @xmath482 thus , using the shifted @xmath78 parameters @xmath483 [ equation  ( [ etam : etap ] ) ] we arrive at @xmath484 for @xmath107 ( the choice in the classical case ) , @xmath485 and @xmath486 , so that @xmath487 .      using the product property @xmath489 , we have @xmath490 . then , replacing @xmath67 by its `` bar '' form ( [ ppb ] ) and writing @xmath491 [ from  ( [ bm - bp ] ) ] , we find @xmath492 here we have also used @xmath493 and the commutation rule @xmath494=1 $ ] to get a _ normally ordered _ form . if we introduce some compact notation for the coefficients @xmath495 we can finally write @xmath496 ( the choice @xmath107 gives the self - adjoint form @xmath497 @xcite . )      using equations  ( [ pxb ] ) and  ( [ ppb ] ) for @xmath474 and @xmath475 , as well as @xmath491 and the definition ( [ dmdp ] ) of @xmath499 we get for the kinetic part @xmath500 : @xmath501 the coefficients of @xmath90 and @xmath91 are @xmath502 and @xmath503 , also defined in ( [ dmdp ] ) , and we have @xmath504      to obtain @xmath506 [ vd .  equation  ( [ wkk : lir : lkin : lv ] ) ] , we use the product property ( [ bar : product ] ) to raise @xmath475 to @xmath507 , getting @xmath508 to get later the matrix elements of @xmath94 between hermite functions it is convenient to cast ( [ lvb ] ) in a normally ordered form . to this end we use pathak s theorems @xcite in the following generalised form @xmath509 ! } \ ; \frac{(\etapro/2)^{q}}{q ! } \;.\ ] ] here @xmath510 with @xmath511 $ ] ( @xmath512 is the result of moving the @xmath513 to the left as if they were scalars ) . the generalisation resides in letting @xmath514 and @xmath513 have a constant but different from one commutator . in our case @xmath515=\etapro % % % % % % % % % % % % % % % % % % % \;,\ ] ] follows from the commutator @xmath494=1 $ ] with @xmath516 . note that @xmath514 and @xmath513 are not the adjoint of each other , but the notation is symbolic .    with help from equation  ( [ pathak:2:gral : odd : alt : rep ] ) and @xmath517 [ equation  ( [ kappa : s ] ) ] we can write ( [ lvb ] ) as follows @xmath518 where we have introduced the notations @xmath519 ! } \frac{(\etapro/2)^{q}}{q ! } \bv^{(2\iq+1 ) } % % % % % % % % % % % % % % % % % % \ ; , \qquad % % % % % % % % % % % % % % % % % % \hat{x}^{\ell } = \nol\big(\aplus+\aminus\big)^{\ell}\nor % % % % % % % % % % % % % % % % % % \;.\ ] ] for a given @xmath520 different powers of @xmath521 appear . we take a common factor @xmath522 and add coefficients , e.g. , @xmath523 with @xmath521 , @xmath524 , with @xmath525 , etc . @xmath526 for @xmath527 in @xmath528 we write @xmath529}=\px^{2q}\bv^{(2\ell+1)}$ ] and introduce the operator @xmath530 @xmath531 now the @xmath532 and @xmath533 dependences factorise , so we can sum the series of @xmath532 in @xmath534 , getting an exponential . then , recalling the definition of @xmath103 we finally obtain the normally ordered form for @xmath94 we were looking for @xmath535 the operator @xmath536 does not act on what could appear after @xmath15 , since @xmath530 was introduced as some abreviated notation for the high - order derivatives of the potential . here we calculate the matrix elements @xmath537 , with the transformed operator @xmath86 obtained in [ app : lb ] . we split in parts the calculation by introducing the corresponding decomposition : @xmath538 , @xmath539 , and @xmath540 , where we have used the custom `` bra''``ket '' notation . the calculation is simplified because we have all operators @xmath470 expressed as normally ordered forms of @xmath90 and @xmath91 . then we can take advantage of the ladder actions @xmath541 , @xmath542 , and the number property @xmath543 .      from  ( [ lirb ] ) for @xmath92 and the orthonormality of the @xmath545 , we find for @xmath546 @xmath547 then , we introduce the notation @xmath548 for the coefficient of the diagonal term [ cf.equation  ( [ gammas ] ) ] , multiply by @xmath84 and sum over @xmath549 , getting @xmath550      using equation  ( [ lkinb ] ) for @xmath93 and the `` ladder '' action of @xmath90 and @xmath91 , we have @xmath552 where we have used the kronecker delta to interchange @xmath113 and @xmath553 . then , multiplying by @xmath84 and summing over @xmath549 , one obtains @xmath554      this will take longer . first , the operator @xmath556 appearing in equation ( [ lvb : norord : el ] ) for @xmath94 can be put in explicit form by using the binomial formula @xmath557 the @xmath514 and @xmath513 are proportional to @xmath90 and @xmath91 [ equation  ( [ am - ap ] ) ] , so the contribution of each summand to @xmath558 involves @xmath559 . from orthonormality only one term of the sum ( [ sum : pathak ] ) will contribute , namely @xmath560 , while @xmath561 will be odd , @xmath562 with @xmath563 ( reflecting the odd powers of @xmath67 in @xmath89 ) . then , in terms of @xmath520 , the index @xmath564 is restricted to @xmath565    [ [ relating - the - results - for - ippipmp2iq1 . ] ]   relating the results for @xmath562 . + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    the two cases can be easily related . first , the combinatorial coefficients are equal @xmath566 , because @xmath567 . second , the @xmath78 factors , coming from the proportionality of the @xmath568 , to @xmath569 , are simply obtained from each other by exchanging pluses and minuses ( @xmath516 ) : @xmath570 third and last , the matrix elements are also relatable . for @xmath101 , we have @xmath571 , while for @xmath572 @xmath573 which is formally identical to the former by replacing @xmath574 and @xmath575 . thus , from now on we consider only the case @xmath101 ( factor @xmath576 ) . the result for @xmath572 , will be readily obtained from this by replacing @xmath577 from all these considerations , we can write the matrix elements of @xmath556 as @xmath578 where we omit a kronecker delta ensuring @xmath101 . the rest of the calculation consists of multiplying this by @xmath579 and sum over all @xmath533 to get @xmath580 [ equation  ( [ lvb : norord : el ] ) ] .    [ [ restrictions - on - the - sum - on - ell . ] ]   restrictions on the sum on @xmath533 . + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    the contribution of the infinite sum ( [ lvb : norord : el ] ) to @xmath580 is fortunately cut - off . first , the lower index in @xmath581 should be positive and , introducing the shifted index @xmath582 , we obtain ( we rename @xmath583 at the end ) @xmath584 here we have already used that the sum is also cut - off from above at @xmath585 , since @xmath586 , for @xmath587 . on the other hand , for @xmath588 [ equation  ( [ lvb : norord : el ] ) ] we have @xmath589 ! } ( -\ldb^{2}\px^{2})^{\ell } \frac{(-1)^{\iq}\ldb^{2\iq}}{(2\iq+1 ) ! } \bv^{(2\iq+1 ) } % % % % % % % % % % % % % % % % % % \;.\ ] ] combining this with @xmath590 , gathering @xmath591 with @xmath592 to form the operator @xmath530 [ equation  ( [ aqlq ] ) ] , and recalling the definition  ( [ kappa : s ] ) of @xmath103 , we obtain @xmath593 ! } \ ; \big\langle\ip\big| ( \bp)^{(2\iq+1)+\ell}\bd^{\ell } \big|\ipp\big\rangle \ ; \frac{(-\dlpxx)^{\ell}}{\ell ! } \ ; \bv^{(2\iq+1 ) } % % % % % % % % % % % % % % % % % % \;.\ ] ] note that the quotient of factorials equals @xmath594 , with the pochhammer symbol defined by ( * ? ? ? 13.5 ) @xmath595    [ [ computation - of - the - matrix - element . ] ]   computation of the matrix element . + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    to conclude , we only need the matrix element @xmath596 in ( [ lvb : matele:1 ] ) . we first pass @xmath507 times the action of @xmath91 to @xmath597 by taking its adjoint @xmath90 and then we use repeteadly its downward ladder action , getting @xmath598 with @xmath101 . then , in terms of the pochhammer symbol ( [ pochhammer ] ) , we can write @xmath599 now , taking into account that the factor @xmath600 cancels that of @xmath601 and using the modified coefficient @xmath602 [ equation  ( [ gamma : qcoeff ] ) ] , we arrive at @xmath603 \bv^{(2\iq+1 ) } % \;\delta_{\ip-(2\iq+1),\ipp } % % % % % % % % % % % % % % % % % % \;.\ ] ] comparing with the series definition of the confluent hypergeometric ( kummer ) function @xmath102 ( * ? ? ? 13.6 ) we get @xmath604 \bv^{(2\iq+1 ) } % % % % % % % % % % % % % % % % % % \quad % % % % % % % % % % % % % % % % % % \kum(a , c\,;z ) = \sum_{\ell=0}^{\infty } \frac{(a)_{\ell}}{(c)_{\ell } } \frac{z^{\ell}}{\ell ! } % % % % % % % % % % % % % % % % % % \;.\ ] ] the operator in front of @xmath76 is simply @xmath605 [ equation  ( [ gamma : qcoeff ] ) ] , so that we can finally write @xmath606 note that @xmath607 is a polynomial because the first negative index cuts the upper pochhammer symbol . indeed the `` integer - plus - one '' second argument @xmath608 tells us that it is an associated laguerre polynomial ( * ? ? ? 13 ) @xmath609 in our case @xmath610 and @xmath611 , so that @xmath612 .    [ [ final - expressions - for - sum_ippq_ipipprm - vec_ipp . ] ]   final expressions for @xmath555 . + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    once we have computed @xmath613 , the matrix elements of @xmath94 , we multiply by @xmath84 and sum over @xmath549 ( i.e. , over @xmath520 ) , to get @xmath614 } \big [ \g_{\ip}^{\iq,- } \ , \bv^{(2\iq+1 ) } \big ] \ , \ec_{\ip-(2\iq+1 ) } % % % % % % % % % % % % % % % % % % \;.\ ] ] note that this @xmath520-sum is restricted because @xmath615 . this yields the upper limit @xmath616 $ ] , with @xmath617 $ ] the integer part of @xmath618 . the contribution from @xmath619 is easily obtained recalling the transformation rules ( [ rule ] ) ( here the sum is not restricted ) @xmath620 \ , \ec_{\ip+(2\iq+1 ) } % % % % % % % % % % % % % % % % % % % \;,\ ] ] with @xmath621 given by ( [ gamma : qcoeff ] ) . we have enclosed into square brackets the action of the operators @xmath622 to recall that they only act on the @xmath8-dependence of the potential .    for @xmath107 , the choice in the classical case , @xmath485 , @xmath486 , and @xmath623 . then , @xmath624 and @xmath625 , so that only the @xmath626 terms survive . the chain of equations then adquires an unbalanced structure ( like in the old hierarchies @xcite ) , which may result in a poor stability when the quantum terms are important . in constrast , @xmath214 gives @xmath627 , and hence the weight of the terms at both sides of @xmath113 ( @xmath305 above and below the diagonals in the matrices @xmath158 ) is similar , improving the stability when going into the deep quantum regime . for special @xmath78 and @xmath189 , the integrals @xmath628 and @xmath180 [ equation  ( [ ial : knl ] ) ] appearing in the expressions ( [ j : ial : knl ] ) for the observables , are easily done . for @xmath107 , we simply have @xmath629 , @xmath630 , and @xmath631 . if @xmath632 ( i.e. , @xmath216 ) , we get @xmath633 . then the first moments reduce to @xmath634 ( normalisation ) and to @xmath635 and @xmath636 . in the general case , one can derive recurrences relating the @xmath180 with lower order ( in @xmath533 ) ones . thus , only @xmath637 is needed , which can be found analytically . concerning the @xmath628 , we mostly use @xmath216 , and hence @xmath633 . anyway , for periodic potentials one can derive recurrence relations for them which can also be solved by continued fractions .      to get the recurrence for the @xmath180 we start from their definition ( [ ial : knl ] ) with index @xmath638 , express the last @xmath10 as @xmath491 [ equation  ( [ bm - bp ] ) ] , and then use @xmath639 , getting @xmath640 this recurrence can be used to get the @xmath641 from the @xmath180 . we can get analytical expressions for the first few @xmath180 with help from the tabulated integral ( * ? ? ? * equation  ( 7.374 - 4 ) ) @xmath642 for @xmath643 , we use this result with @xmath644 , obtaining @xmath645 although @xmath646 and @xmath647 can now be obtained from ( [ knl : rr ] ) , they can also be done with help from equation  ( [ hermite : func : int:1 ] ) , getting : @xmath648 , @xmath649 , and @xmath650      for a general periodic potential we can derive a recurrence relation for the integrals @xmath651\,u_{\ix}(\x)$ ] [ equation  ( [ ial : knl ] ) ] . we integrate by parts the expression for @xmath652 and use @xmath653 with @xmath654 , getting @xmath655 for a finite number of harmonics ( @xmath656 , if @xmath657 ) , the recurrence ( [ rr : ial : vper ] ) has a finite coupling range . for instance , for the ratchet potential ( [ vrat ] ) , we have @xmath307 and @xmath308 and the above recurrence reduces to @xmath658 % % % % % % % % % % % % % % % % % % \ ; , \qquad % % % % % % % % % % % % % % % % % % \epspott=\epspot\,\vo/2\kt % % % % % % % % % % % % % % % % % % \;.\ ] ] with @xmath659 as input ( obtained numerically , e.g. , by simpson rule ( * ? ? ? * app . 2 ) ) , this recursion can be solved by continued fractions ( [ app : rr - cf ] ) . here we write expliticly the matrices @xmath158 for arbitrary periodic potentials . we shall give the results for a slightly generalised basis , obtained by shifting the origin of momenta by an amount @xmath660 [ cf.equation  ( [ w : expansion : xp ] ) ] @xmath661 this momentum shift is convenient to `` catch '' solutions centred far from zero @xmath10 , where the ordinary basis would need many @xmath80 to reconstruct the distribution . ( a useful choice is , when generating a curve , to set @xmath662 of the precedent point . ) the momentum shift is handled as a change of variable @xmath663 . then , @xmath664 , but the term @xmath62 makes the force to enter in the combination ( scaled units ) @xmath665 as mentioned in section  [ sec : vper ] , we extract the force @xmath187 from @xmath15 and set the auxiliary potential @xmath189 proportional to the periodic part @xmath77 . the @xmath158 are obtained by inserting the formulae for @xmath666 , and @xmath667_{\ix\ixp}$ ] [ equation  ( [ dmdpvper : matr - elem ] ) ] , into the general matrices of section  [ sec : mqbh ] . instead of giving a single expression , we write separately the matrices for @xmath668 and @xmath669 . the central matrix @xmath670 is determined by the terms involving @xmath671 , and reads @xmath672 % % % % % % % % % % % % % % % % % % line-1 \sqrt{1}(\iu\,\dmi\ix\!-\!\etam\bfsh ) \!\!\!\!&\!\!\!\ ! \gammao_{1 } \!\!+\!\ ! \pshift\,\iu\ix \!\!\!\!&\!\!\!\ ! \sqrt{2}(\iu\,\dpl\ix\!-\!\etap\bfsh ) \!\!\!\!&\!\!\!\ ! \sqrt{2\cdot3}\ , \gammapp \!\!\!\!&\!\!\!\ ! 0 \!\!\!\!&\!\ ! \ddots \\[0.ex ] % % % % % % % % % % % % % % % % % % line-2 \sqrt{1\cdot2}\ , \gammamm \!\!\!\!&\!\!\!\ ! \sqrt{2}(\iu\,\dmi\ix\!-\!\etam\bfsh ) \!\!\!\!&\!\!\!\ ! \gammao_{2 } \!\!+\!\ ! \pshift\,\iu\ix \!\!\!\!&\!\!\!\ ! \sqrt{3}(\iu\,\dpl\ix\!-\!\etap\bfsh ) \!\!\!\!&\!\!\!\ ! \sqrt{3\cdot4}\ , \gammapp \!\!\!\!&\!\ ! \ddots \\[0.ex ] % % % % % % % % % % % % % % % % % % line-3 0 \!\!\!\!&\!\!\!\ ! \sqrt{2\cdot3}\ , \gammamm \!\!\!\!&\!\!\!\ ! \sqrt{3}(\iu\,\dmi\ix\!-\!\etam\bfsh ) \!\!\!\!&\!\!\!\ ! \gammao_{3 } \!\!+\!\ ! \pshift\,\iu\ix \!\!\!\!&\!\!\!\ ! \sqrt{4}(\iu\,\dpl\ix\!-\!\etap\bfsh ) \!\!\!\!&\!\ ! \ddots \\[0.ex ] % % % % % % % % % % % % % % % % % % line-4 0 \!\!\!\!&\!\!\!\ ! 0 \!\!\!\!&\!\!\!\ ! \sqrt{3\cdot4}\ , \gammamm \!\!\!\!&\!\!\!\ ! \sqrt{4}(\iu\,\dmi\ix\!-\!\etam\bfsh ) \!\!\!\!&\!\!\!\ ! \gammao_{4 } \!\!+\!\ ! \pshift\,\iu\ix \!\!\!\!&\!\ ! \ddots \\[0.ex ] % % % % % % % % % % % % % % % % % % line-6 ( 5 ) \ddots \!\!\!\!&\!\!\!\ ! \ddots \!\!\!\!&\!\!\!\ ! \ddots \!\!\!\!&\!\!\!\ ! \ddots \!\!\!\!&\!\!\!\ ! ! \ddots \end{array } \right ) % % % % % % % % % % % % % % % % % % \ ] ] with @xmath673 and @xmath674 given by ( [ gammas ] ) , @xmath675 , and @xmath676 . here the periodic potential does not appear ; only @xmath187 enters in the matrices @xmath670 . besides , they do depend on the index @xmath131 ( cf . below ) . the matrices @xmath677 are determined by the terms involving @xmath678 in equation  ( [ dmdpvper : matr - elem ] ) and read @xmath679 % % % % % % % % % % % % % % % % % % line-1 \etam \gs_{1}^{0 } \!-\ ! \sqrt{1}\,\epspot \!\!\!&\!\!\ ! -\pshift\,\epspot \!\!\!&\!\!\ ! \etap \gs_{2}^{0 } \!-\ ! \sqrt{2}\,\epspot \!\!\!&\!\!\ ! 0 \!\!\!&\!\!\ ! \etap^{3 } \gs_{4}^{1 } \!\!\!&\!\!\ ! 0 \!\!\!&\!\!\ ! \ddots \\[0.ex ] % % % % % % % % % % % % % % % % % % line-2 0 \!\!\!&\!\!\ ! \etam \gs_{2}^{0 } \!-\ ! ! -\pshift\,\epspot \!\!\!&\!\!\ ! \etap \gs_{3}^{0 } \!-\ ! \sqrt{3}\,\epspot \!\!\!&\!\!\ ! 0 \!\!\!&\!\!\ ! \etap^{3 } \gs_{5}^{1 } \!\!\!&\!\!\ ! \ddots \\[0.ex ] % % % % % % % % % % % % % % % % % % line-3 \etam^{3 } \gs_{3}^{1 } \!\!\!&\!\!\ ! 0 \!\!\!&\!\!\ ! \etam \gs_{3}^{0 } \!-\ ! \sqrt{3}\,\epspot \!\!\!&\!\!\ ! -\pshift\,\epspot \!\!\!&\!\!\ ! \etap \gs_{4}^{0 } \!-\ ! \sqrt{4}\,\epspot \!\!\!&\!\!\ ! 0 \!\!\!&\!\!\ ! \ddots \\[0.ex ] % % % % % % % % % % % % % % % % % % line-4 0 \!\!\!&\!\!\ ! \etam^{3 } \gs_{4}^{1 } \!\!\!&\!\!\ ! 0 \!\!\!&\!\!\ ! \etam \gs_{4}^{0 } \!-\ ! \sqrt{4}\,\epspot \!\!\!&\!\!\ ! -\pshift\,\epspot \!\!\!&\!\!\ ! \etap \gs_{5}^{0 } \!-\ ! \sqrt{5}\,\epspot \!\!\!&\!\!\ ! \ddots \\[0.ex ] % % % % % % % % % % % % % % % % % % line-5 \etam^{5 } \gs_{5}^{2 } \!\!\!&\!\!\ ! 0 \!\!\!&\!\!\ ! \etam^{3 } \gs_{5}^{1 } \!\!\!&\!\!\ ! 0 \!\!\!&\!\!\ ! \etam \gs_{5}^{0 } \!-\ ! \sqrt{5}\,\epspot \!\!\!&\!\!\ ! -\pshift\,\epspot \!\!\!&\!\!\ ! \ddots \\[0.ex ] % % % % % % % % % % % % % % % % % % line-6 \ddots \!\!\!&\!\!\ ! \ddots \!\!\!&\!\!\ ! \ddots \!\!\!&\!\!\ ! \ddots \!\!\!&\!\!\ ! \ddots \!\!\!&\!\!\ ! ! \ddots \end{array } \right ) % % % % % % % % % % % % % % % % % % \ ] ] where @xmath211 is evaluated at @xmath680 , with @xmath516 . these matrices enjoy properties opposite to those of @xmath670 . the substrate potential does appear ( multiplying all elements in front of the matrices ) while the force is absent . besides , the matrices @xmath677 are independent of the index @xmath131 , which permits to write them simply as @xmath681 . the only technical complication with respect to the common choice @xmath107 is that we need to compute the kummer functions @xmath102 included in @xmath211 , which in our case are simple polynomials [ equation  ( [ assoc : laguerre : poly : kummer ] ) ] . note that to compute the observables we must eventually undo the momentum shift @xmath682 . in this appendix we derive some semiclassical formulae discussed in the main text . in section  [ ratqbm ] we invoked the current through a saw - tooth barrier deformed by @xmath683 ( figure [ fig : saw - width ] ) . ( b ) and by @xmath684 ( c ) . right panel : width of the energy bands @xmath685 in a cosine potential ( exact and various approximations ) . , title="fig : " ]   ( b ) and by @xmath684 ( c ) . right panel : width of the energy bands @xmath685 in a cosine potential ( exact and various approximations ) . , title="fig : " ]    we introduced the reflexion coefficients @xmath359 and @xmath360 for @xmath312 . one has to analyse separately the cases of incident particles with energies above or below the barrier . in the first case , the reflexion coefficient by a potential with a cusp ( with slopes @xmath355 and @xmath356 ) is given by the first equation  ( [ rs ] ) @xcite . owing to the right slope in figure  [ fig : saw - width ] is infinity , we make it finite and finally take the vertical limit , getting @xmath686 next we consider the case of tunnelling through the barrier . the semiclassical transmission coefficient is also written in ( [ rs ] ) . defining @xmath687 , we calculate @xmath688 and @xmath689 for positive and negative forces . after simple integrals , we find @xmath690 with @xmath357 the particle energy ( cf . reference  @xcite ) . for moderate forces , @xmath691 , we get for the transmission coefficients @xmath692 $ ] @xmath693 thus , for @xmath357 below the barrier @xmath363 , while @xmath361 for @xmath357 above the barrier .      in order to get the semiclassical transmision coefficient through a potential @xmath405 , we need to do the integral @xmath694 } \big| % % % % % % % % % % % % % % % % % % % \;,\ ] ] with turning points obeying @xmath695 . due to @xmath696 below the barrier , we change the sign inside the square root extracting an @xmath697 which is absorbed by the modulus . simple tranformations bring the integral into a tabulated form ( * ? ? ? * equation  ( 2.576 ) ) , and we obtain @xmath698 % % % % % % % % % % % % % % % % % % \ ; , \qquad % % % % % % % % % % % % % % % % % % \varepsilon = e/\vo \ ; , \quad s=\sqrt{\half(1-\varepsilon ) } % % % % % % % % % % % % % % % % % % % \;,\ ] ] where @xmath699 while @xmath700 and @xmath701 the first and second complete elliptic integrals ( * ? ? ? insertion in @xmath692 $ ] gives the transmision coefficient through a cosine potential . this @xmath353 also gives the _ semiclassical energy bands _ @xcite @xmath702 with @xmath703 the oscillation frequency in the wells . the width of the bands is therefore @xmath704 \big\ } % % % % % % % % % % % % % % % % % % \ ; , \qquad % % % % % % % % % % % % % % % % % % \varepsilon = e/\vo % % % % % % % % % % % % % % % % % % \;.\ ] ] the parameter @xmath705 goes from @xmath202 when @xmath706 ( potential bottom ) to @xmath336 at the barrier top @xmath707 . for @xmath708 we can expand the elliptic integrals , getting @xmath709 \big\ } \nonumber\\ % % % % % % % % % % % % % % % % % % & + ( 1+\varepsilon)^{2 } \big\ { 1 + \case{1}{8 } \ln \big [ \case{1}{32 } ( 1+\varepsilon ) \big ] \big\ } \big ) \big ] % % % % % % % % % % % % % % % % % % \;.\end{aligned}\ ] ] at the potential bottom @xmath706 this gives @xmath710 $ ] . the band width increases exponentially with the energy . near the top , @xmath711 , one has @xmath712 , and we can use the series expansion of the elliptic integrals , obtaining @xmath713 % % % % % % % % % % % % % % % % % % \;.\ ] ] figure  [ fig : saw - width ] shows that equations ( [ wb : approx : bottom ] ) and ( [ wb : approx : top ] ) approximate well the exact results in their respective ranges ( the later also serves as un upper bound ) . if we disregard the terms with @xmath714 in ( [ wb : approx : bottom ] ) , the result happens to fit remarkably @xmath685 in all the range . in this last appendix we shall solve the following generic differential equation @xmath715 perturbatively in the forcing @xmath716 . ( for cosine forcing @xmath717 . ) comparing with section  [ sec : dynamics ] , we see that @xmath718 corresponds to the vectors @xmath719 , @xmath720 to the static part of the matrices @xmath158 , and the perturbation @xmath721 to @xmath722 , the part of @xmath670 involving @xmath386 .    a periodic solution can be fourier expanded as @xmath723 with @xmath724 the unperturbed response , while @xmath725 for complex @xmath718 . separating terms oscillating with @xmath726 and @xmath727 , the left - hand side of ( [ eom : generic ] ) reads @xmath728 \e^{+\iu\omega_{\iw } t } \\ % % % % % % % % % % % % % % % % % % & + & \sum_{\iw=1}^{\infty } \dfm^{\iw}\ , \big [ -\iu(\omega_{\iw}\tau)\tilde{\ec}^{(\iw ) } + \rel\,\tilde{\ec}^{(\iw ) } \big ] \e^{-\iu\omega_{\iw } t } % % % % % % % % % % % % % % % % % % \;.\end{aligned}\ ] ] to obtain the right - hand side , we use @xmath729 , redefine the indices ( keeping the same names ) to get all the oscillating factors at @xmath730 , and introduce the definition @xmath731 , arriving at @xmath732 \e^{+\iu\omega_{\iw } t } \\ % % % % % % % % % % % % % % % % % % & + & \sum_{\iw=1}^{\infty } \dfm^{\iw } \big [ \vpert\,\tilde{\ec}^{(\iw-1 ) } + \left(\dfp\dfm\right ) \vpert\,\tilde{\ec}^{(\iw+1 ) } \big ] \e^{-\iu\omega_{\iw } t } % % % % % % % % % % % % % % % % % % \;.\end{aligned}\ ] ] equating terms with the same oscillating factors at both hand sides ( _ uniqueness _ of the fourier expansion ) , we obtain @xmath733 \\[1.ex ] % % % % % % % % % % % % % % % % % % \iu(\omega_{\iw}\tau)\ec^{(\iw ) } + \rel\,\ec^{(\iw ) } & = & \vpert\ , \big[\ ; \ec^{(\iw-1 ) } & + & \left(\dfp\dfm\right ) \ec^{(\iw+1 ) } \big ] % % % % % % % % % % % % % % % % % % \end{array } % \;,\ ] ] the @xmath734 are included as @xmath735 with @xmath736 . the contribution @xmath737 comes from products of  rotating " terms @xmath738 . for instance , the oscillating factor @xmath739 in @xmath740 when multiplied by the part @xmath741 of the field raises the harmonic from @xmath742 to @xmath743 . the term @xmath744 comes from products of  counter - rotating " terms @xmath745 . for example , the product of the oscillating factor @xmath746 and the @xmath747 part of the field lowers the order from @xmath748 to @xmath743 . this contribution ( the  contamination " from higher harmonics ) has an order higher in @xmath716 than the contribution of @xmath740 , which is therefore the leading one . indeed , expanding the @xmath749 in powers of @xmath716 we get the leading terms at each harmonic @xmath750 the @xmath751 equation gives the static response [ corresponding to ( [ 3rr : statics ] ) ] , @xmath388 the linear dynamical response , etc . [ corresponding to equation  ( [ 3rr : harmonics ] ) ] .
Abstract:

Answer: the continued - fraction method to solve classical fokker planck equations has been adapted to tackle quantum master equations of the caldeira  leggett type . this can be done taking advantage of the phase - space ( wigner ) representation of the quantum density matrix . the approach differs from those in which some continued - fraction expression is found for a certain quantity , in that the full solution of the master equation is obtained by continued - fraction methods . this allows to study in detail the effects of the environment ( fluctuations and dissipation ) on several classes of nonlinear quantum systems . we apply the method to the canonical problem of quantum brownian motion in periodic potentials both for cosine and ratchet potentials ( lacking inversion symmetry ) .