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 path integral ( pi ) has been used in quantum physics since the revolutionary work of feynman  @xcite , although the basic observation goes back to dirac  @xcite who appreciated the rle of the lagrangian in short - time evolution of the wave function , and even suggested the time - slicing procedure for finite , i.e. , non - infinitesimal , time lags . since then the pi approach yielded invaluable insights into the structure of quantum theory  @xcite and provided a viable alternative to the traditional operator - formalism - based canonical quantization . during the second half of the 20th century , the pi became a standard tool in quantum field theory  @xcite and statistical physics  @xcite , often providing the easiest route to derivation of perturbative expansions and serving as an excellent framework for ( both numerical and analytical ) non - perturbative analysis  @xcite . feynman pi has its counterpart in pure mathematics , namely , in the theory of continuous - time stochastic processes  @xcite . there the concept of integration over a space of continuous functions ( so - called fluctuating paths or sample paths ) had been introduced by wiener  @xcite already in 1920 s in order to represent and quantify the brownian motion . interestingly enough , this so - called wiener integral ( or integral with respect to wiener measure ) was formulated 2 years before the discovery of the schrdinger equation and 25 years before feynman s pi formulation . the _ local time _ for a brownian particle ( in some literature also called _ sojourn time _ ) has been of interest to physicists and mathematicians , since the seminal work of paul  lvy in 1930 s  @xcite . in its essence , the local time characterizes the time that a sample trajectory @xmath1 of a given stochastic process spends in the vicinity of an arbitrary point @xmath2 . this in turn defines a sample trajectory @xmath3 of a new stochastic process . a rich theory has been developed for local - time processes that stem from diffusion processes ( see , e.g. , ref . @xcite and citations therein ) . for later convenience , we should particularly highlight the ray  knight theorem which states that the local time of the wiener process can be expressed in terms of the squared bessel process  @xcite . in contrast to mathematics , the concept of the local time is not uniquely settled in physics literature . various authors define essentially the same quantity under different names ( local time , occupation time , traversal time , etc . ) , and with different applications in mind . for example , in ref . @xcite the _ traversal time _ is used to study quantum scattering and tunneling processes , in  @xcite the small - temperature behavior of the equilibrium density matrix is analyzed with a help of the _ occupation time _ , while in ref . @xcite the large - time behavior of path integrals that contain functionals of the _ local time _ is discussed . the aim of this paper is to derive a local - time pi representation of the bloch density matrix , i.e. , the matrix elements @xmath4 of the gibbs operator . this can serve not only as a viable alternative to the commonly used feynman  kac representation but also as a powerful tool for extracting both large and small - temperature behavior . apart from the general theoretical outline , our primary focus here will be on the low - temperature behavior which is technically more challenging than the large - temperature regime . in fact , the large - temperature expansion was already treated in some detail in our previous paper  @xcite . last but not least , we also wish to promote the concept of the local time which is not yet sufficiently well known among the path - integral practitioners . the structure of the paper is as follows . to set the stage we recall in the next section some fundamentals from the feynman pi which will be needed in later sections . in section [ sec : heuristic ] we provide motivation for the introduction of a local time , and construct a heuristic version of the local - time representation of pi s . the key technical part of the article is contained in section  [ sec : lt ] , where we derive by means of the replica trick the local - time representation of the bloch density matrix . relation to the sturm  liouville theory is also highlighted in this context . a local - time analog of the feynman  salam formula  @xcite is presented in section  [ sec : functionals ] and its usage is illustrated with a computation of the one - point distribution of the local time . since a natural arena for local - time pi s is in thermally extremal regimes , we confine our attention in section  [ sec : asymp ] on large- and small-@xmath5 asymptotic behavior of the bloch density matrix . there we also derive an explicit leading - order behavior in large-@xmath5 ( i.e. , low - temperature ) expansion . the analysis is substantially streamlined by using the rayleigh  ritz variational principle . finally , section  [ sec : conclusion ] summarizes our results and discusses possible extensions , applications , and future developments of the present work . for the reader s convenience the paper is supplemented with appendix which clarifies some finer technical details . consider a non - relativistic one - dimensional quantum - mechanical system described by a time - independent hamiltonian @xmath6 where @xmath7 . throughout this paper we will study the matrix elements @xmath8 of the gibbs operator @xmath9 , where @xmath10 is the inverse temperature and @xmath11 is the boltzmann constant . the matrix @xmath12 , known also as the bloch density matrix , is a fundamental object in quantum statistical physics , as the expectation value of an operator @xmath13 at the temperature @xmath14 can be written in the form @xmath15 where @xmath16 is the partition function of the system . in case of need , ensuing quantum mechanical transition amplitudes can be obtained from ( [ rhodef ] ) via a wick rotation which formally amounts to the substitution @xmath17 , converting thus the gibbs operator @xmath9 to the quantum evolution operator @xmath18 . the matrix elements in eq.([rhodef ] ) can be represented via the path integral as  @xcite @xmath19 \right\}\ , .\ ] ] this represents a  sum " over all continuous trajectories @xmath20 , @xmath21 $ ] , connecting the initial point @xmath22 with the final point @xmath23 . it should be noted that the integral @xmath24 $ ] is the classical euclidean action integral along the path @xmath20 with @xmath25 . in the following we will denote the euclidean action as @xmath26 . the integrand in @xmath26 , i.e. , @xmath27 , can be identified with the classical hamiltonian function , in which the momentum @xmath28 is substituted for @xmath29 . one can also regard ( [ pix ] ) as an expectation value of the functional @xmath30 $ ] over the ( driftless ) brownian motion with the diffusion coefficient @xmath31 , and duration @xmath32 , that starts at point @xmath33 , and terminates at @xmath34 . the latter stochastic process is also known as brownian bridge . the purpose of this section is twofold . firstly , we would like to motivate a need for reformulation of pi s in the language of local - time stochastic process . in particular , we point out when such a reformulation can be more pertinent than the conventional `` sum over histories '' prescription . secondly , we wish to outline a heuristic construction of the local - time representation of pi s . more rigorous and explicit ( but less intuitive ) formulation of pi s over ensemble of local times will be presented in the subsequent section .    to provide a physically sound motivation for the local - time representation of pi s we follow exposition of paulin _ .  @xcite . to this end we first consider the diagonal elements of the bloch density matrix , i.e. , @xmath35 ( often referred to as the boltzmann density ) . upon shifting @xmath36 , and setting @xmath37 , @xmath38 ( @xmath39 is the thermal de  broglie wavelength ) , the pi ( [ pix ] ) can be reformulated in terms of dimensionless quantities @xmath40 and @xmath41 as @xmath42 \right\ } .\ ] ] note in particular , that in contrast to @xmath43 the measure @xmath44 does not explicitly depend on @xmath5 , and thus @xmath5-dependent parts in the pi are under better control . such a rescaled representation is particularly useful when discussing large- and/or small-@xmath5 behavior of the path integral in question . path fluctuations in the potential are controlled by @xmath45 , and the factor @xmath5 quantifies the significance of the potential @xmath46 with respect to the kinetic term .    for small @xmath5 ( i.e. , high temperature ) , typical paths @xmath47 stay in the vicinity of the point @xmath33 , as depicted in fig . [ fig : paths ] ,     are plotted as functions of the dimensionless time @xmath40 . the solid green path , representing a typical trajectory with a high value of @xmath5 , exhibits large fluctuations , whereas the dashed red path , corresponding to small @xmath5 , stays in the vicinity of the initial and final point @xmath33 . on the right , two local - time profiles @xmath48 are shown . the broad one ( solid green ) arises from the violently fluctuating path @xmath49 , whereas the narrow one ( dashed red ) corresponds to the path with small fluctuations . on the left , we depict a generic potential @xmath50 . ]    and therefore a systematic wigner  kirkwood expansion can be readily developed by taylor - expanding the potential part of the action  @xcite . when @xmath5 is large ( i.e. , low temperature ) , the trajectories @xmath49 fluctuate heavily around the value @xmath33 , and the potential term @xmath46 dominates over the kinetic one . from the statistical physics point of view , the most important contribution to the low temperature behavior of the path integral ( [ pixresc ] ) should come from those paths that spend a sizable amount of time near the _ global _ minimum of the potential @xmath50 . for this reason , it is important to be able to keep track of the time which a given path spends in an infinitesimal neighborhood of an arbitrary point @xmath0 . let us define , for each wiener trajectory @xmath20 present in the feynman path integral ( [ pix ] ) the ensuing _ local time _ as @xmath51 since the local time @xmath52 is a functional of the random trajectory @xmath53 for @xmath54 , it represents a random variable . from the definition ( [ ltdef ] ) we can immediately see that @xmath55 for all @xmath56 , @xmath57 and that @xmath3 has a compact support . in addition , it can be proved  @xcite that local - time trajectories @xmath3 are , with probability one , continuous curves which ( similarly as trajectories in the underlying wiener process ) are nowhere differentiable . in fig . [ fig : paths ] we depict two examples of representative local - time trajectories . an extensive mathematical discussion of properties of the local time can be found , e.g. in refs . @xcite .    with the definition ( [ ltdef ] ) the potential part of the euclidean action can be recast into form @xmath58 . a local - time representation of the bloch density matrix @xmath12 is then given by @xmath59 \ \!\delta \ ! \left ( \int_\mathbb{r } dx \ \ ! l^x   - \beta\hbar \right ) \exp \left[-\frac{1}{\hbar } \int_\mathbb{r }   dx \ \ ! l^x v(x ) \right ] , \ ] ] where the pi `` sum '' is taken over all local - time trajectories @xmath60 , with @xmath0 being the independent variable ( not to be mistaken with the wiener trajectory @xmath20 ) . the @xmath61-function enforces the normalization constraint mentioned above . basically , transition to the local - time description represents a change ( or a functional substitution ) of stochastic variables @xmath62 . the weight factor @xmath63 appearing in ( [ heurform ] ) can be formally written in the form @xmath64   \ = \   \exp \left[-\frac{1}{\hbar}\int_{0}^{\beta\hbar}\!\ ! d\tau \ \ ! \frac{m}{2}\dot{x}^2 \right ] \ \ ! \det\left(\frac{\delta l^x(\beta\hbar)}{\delta x(\tau ) } \right)^{-1}\ , . \label{iii.7.a}\end{aligned}\ ] ] it is a function of @xmath33 , @xmath34 ( which are implicitly present in @xmath65 ) and @xmath5 , and a functional of the local time @xmath60 . of course , these cavalier manipulations do not have more than a heuristic nature , and it is , indeed , a non - trivial task to determine @xmath63 directly from ( [ iii.7.a ] ) . for this reason , we will in the following section tackle this problem indirectly . in this section , we present a derivation of the local - time representation of the bloch density matrix ( [ rhodef ] ) . initially , we limit ourselves to considering the case of the diagonal part , @xmath66 , and arrive at the key result ( [ diaglapl ] ) , which expresses the matrix elements in the ensuing laplace picture with respect to @xmath5 . this intermediate outcome is shown to agree with the sturm  liouville theory . in the next step , we generalize the latter result to the off - diagonal elements ( @xmath67 ) . the inverse laplace transform will then yield the sought local - time representation of pi [ cf . ( [ offdiag ] ) ] . it follows from the definition ( [ rhodef ] ) that the function @xmath12 satisfies the heat equation @xmath68 \rho(x_a , x_b,\beta ) \ = \ 0\ , , \ ] ] with the initial condition @xmath69 . this is merely a wick - rotated ( @xmath70 ) analogue of the schrdinger equation . the feynman  kac formula  @xcite then ensures that the pi ( [ pix ] ) can be calculated by solving corresponding parabolic differential equation ( [ diffusioneq ] ) .    in the laplace picture eq . ( [ diffusioneq ] ) takes the form @xmath71 \widetilde{\rho}(x_a , x_b , e ) \ = \ \delta(x_a - x_b)\ ,   , \ ] ] with @xmath72 . ( [ lapdifeq ] ) implies that @xmath73 is nothing but the green function of the operator @xmath74 . with the benefit of hindsight , we represent the green function @xmath73 as a path integral over fluctuating fields  so - called functional integral  @xcite . this is rather standard strategy in quantum field theory  @xcite , and in our case it yields @xmath75 } } { \int_{\psi(x_{-})=0}^{\psi(x_{+})=0 } \mathcal{d}\psi(x ) \ \ ! e^ { -\frac{1}{2 } \mathcal{a}^e[\psi ]   } } \ ,   , \ ] ] where @xmath76 & \ \equiv \ \int_{x_{-}}^{x_{+ } } \!\!\!dx \ \ ! \psi(x ) \left[-\frac{\hbar^2}{2 m } \frac{d^2}{d x^2 } + v(x ) + e \right ] \psi(x ) \nonumber\\[2 mm ] & \ = \ \int_{x_{-}}^{x_{+ } } \!\!\!dx   \left [ \frac{\hbar^2}{2 m } \psi'(x)^2 + ( v(x ) + e ) \psi(x)^2 \right],\end{aligned}\ ] ] is the euclidean action functional of the field - theoretic path integral . the super - index `` @xmath77 '' in @xmath26 indicates the shift in the potential @xmath50 by the amount @xmath77 . here , we have confined our quantum - mechanical system within a finite box @xmath78 $ ] , with @xmath79 and @xmath80 . a real scalar field @xmath81 satisfies dirichlet boundary conditions @xmath82 so as to ensure the validity of the operations being performed . since we will ultimately want to invert the laplace transform to regain from @xmath83 the original bloch density matrix @xmath12 , we can not treat the denominator in ( [ fieldrepre ] ) as an irrelevant normalization constant ( which is the usual practice in quantum field theory ) but we have to take care of its @xmath77-dependence . to this end we take advantage of the formula @xmath84 which is a simple version of the so - called _ replica trick_. [ the usual replica - trick formula ( cf . e.g. , ref . @xcite ) can be obtained from ( [ rep.1a ] ) by integrating both sides with respect to @xmath85 and subsequently dividing by @xmath86 . ] with the help of ( [ rep.1a ] ) we can rewrite ( [ fieldrepre ] ) as a multidimensional functional integral @xmath87 } \ ,   , \ ] ] where the multiplet @xmath88 is a @xmath89-component  replica " field in @xmath90 dimensions , @xmath91 denotes the scalar product @xmath92 , and we have rescaled the fields by a factor of @xmath93 in passing . the factor @xmath94 results from a pi generalization of the well know mean - value identity @xmath95 valid for any two vectors in @xmath89-dimensional statistically isotropic environments .    as a side remark , note that we may now invert the laplace transform , using the trivial identity @xmath96 to obtain the representation @xmath97 } \ , .\ ] ] upon scaling @xmath98 , this agrees with the formula ( 2.9 ) in ref .  @xcite . though the result ( [ replrepre ] ) generalizes to arbitrary number of dimensions of the @xmath0-space , i.e. , @xmath99 , our further development will be illustrated ( for simplicity s sake ) only on the one - dimensional case .      in order to derive the weight factor @xmath63 ( cf . ( [ iii.7.a ] ) ) we have arrived at the representation ( [ replrepre ] ) with @xmath89 replica fields . this form is still not very transparent , and a further simplification step is needed to get rid of an explicit dependence of the measure on @xmath89 . to this end , we note that @xmath100 $ ] is in fact the action of a @xmath89-dimensional harmonic oscillator with the `` time '' variable @xmath0 , `` position '' variable @xmath101 , mass @xmath102 , and time - dependent `` frequency '' @xmath103 . considering for a moment only diagonal matrix elements , @xmath66 , spherical symmetry in the replica field space allows to reduce the path integral ( [ afterrepl ] ) to its radial part . due to the boundary conditions , @xmath104 , only the zero - angular - momentum ( @xmath40-wave ) contribution is non - vanishing ( generally a weighted sum over radial pi s with different angular momenta would be required ) . this will be rigorously justified at the end of this section . the corresponding _ radial _ pi representation for ( [ afterrepl ] ) reads  @xcite @xmath105 } \ , .\ ] ] here , the radial part @xmath106 of the @xmath89-dimensional replica field @xmath107 is always non - negative , i.e. , @xmath108 ; the area of a unit sphere in @xmath89 dimensions , @xmath109 , may be replaced by its small-@xmath89 asymptotic form @xmath110 ; and @xmath111 have been introduced to regularize the origin of the @xmath107-space . the new action functional @xmath112 \ \equiv \ \mathcal{a}^e[\eta ] \ + \ \int_{x_-}^{x_+ } \!\ ! dx \ \ ! \frac{(d\!-\!1)(d\!-\!3)}{8 \eta^2(x ) } \ , , \ ] ] is the euclidean action functional of the _ radial harmonic oscillator _ it contains an additional centrifugal potential term ( edwards  gulyaev or langer term  @xcite ) , which emerges from bessel function @xmath113 present in the finite sliced form of the radial pi ( [ diaglapl ] ) . at this point we should stress that in contrast to the quantum - mechanical radial pi , one can use safely the asymptotic expansion for the bessel function @xmath113 ( see , e.g. ref . @xcite ) : @xmath114 > 0)\ , , \nonumber \\[2 mm ] & & y_j \ = \ ( m/\varepsilon \hbar)r_jr_{{j-1}}\ , , \label{asymptot - b}\end{aligned}\ ] ] in the euclidean pi sliced form . here the infinitesimal `` time '' slice @xmath115 is related to the number of slices @xmath116 via the relation @xmath117 . in quantum mechanic this is a problematic step because ( [ asymptot - b ] ) requires @xmath118 > 0 $ ] while there @xmath118 = \mbox{re}[(m / i\varepsilon \hbar)r_jr_{{j-1 } } ] = 0 $ ] . fortunately , the pi for radial harmonic oscillator is exactly solvable even in the case of @xmath0-dependent oscillator frequency . the solution reads  @xcite @xmath119 \right\ } \nonumber \\[2 mm ] & = & \ \frac{\hbar^2}{m } \frac{\sqrt{\eta_1 \eta_2}}{g(x_1 ) } i_{d/2 - 1}\left(\frac{\hbar^2}{m } \frac{\eta_1 \eta_2}{g(x_1)}\right ) \exp\left[- \frac{\hbar^2}{2 m } \left(\frac{f'(x_2)}{f(x_2 ) } \eta_2 ^ 2 - \frac{g'(x_1)}{g(x_1 ) } \eta_1 ^ 2 \right ) \right ] . \label{radialhosol}\end{aligned}\ ] ] the functions @xmath120 and @xmath121 are two independent solutions of the differential equation @xmath122y(x ) \ = \ \left [ -\frac{\hbar^2}{2 m } \frac{d^2}{dx^2 } + v(x ) + e \right ] y(x ) \ = \ 0\ , , \ ] ] with the initial conditions @xmath123 and @xmath124 , and @xmath125 and @xmath126 . in addition , the wronskian @xmath127 is independent of @xmath0 , as can be proved by differentiation and by using the fact that @xmath128 and @xmath129 both satisfy eq . ( [ diffeqfg ] ) . by equating the values of @xmath130 at points @xmath131 and @xmath132 , and taking into account the initial conditions for @xmath128 and @xmath129 , we find a useful identity @xmath133 . now , the pi in ( [ diaglapl ] ) can be sliced at point @xmath33 , and expressed as @xmath134 the limits in eq . ( [ diaglapl ] ) are readily carried out with the help of the asymptotic formulas @xmath135 , and @xmath136 , valid for @xmath137 . subsequent integration over @xmath138 brings ( [ diaglapl ] ) to form @xmath139 where @xmath140 solves eq . ( [ diffeqfg ] ) with initial conditions @xmath141 and @xmath142 , and @xmath143 solves the same equation with @xmath144 and @xmath145 . the denominator in ( [ diaggreenexpl ] ) is the wronskian @xmath146 . the full derivation is given in appendix  [ sec : app ] . although rather explicit , eq . ( [ diaggreenexpl ] ) is not well suited for the laplace transform inversion , since functions @xmath120 and @xmath121 contain @xmath77 in a non - trivial way , which , in addition , significantly hinges on the actual form of @xmath50 . for formal manipulations it is still better the employ the pi representation ( [ diaglapl ] ) . for instance , using eq . ( [ inverselap ] ) we can easily invert the laplace transform to return from @xmath77 back to the @xmath5-variable , namely @xmath147 } \ , . \label{boltzmann.ii}\end{aligned}\ ] ] note that we have utilized the asymptotic form @xmath110 which holds for @xmath148 . we shall see shortly that ( [ boltzmann.ii ] ) can be straightforwardly related to the local - time pi representation of the boltzmann density matrix . let us finally comment on the higher - angular - momentum terms which , as claimed , should not contribute to the expression ( [ diaglapl ] ) . for arbitrary angular momentum @xmath149 , we employ formula ( [ radialhosol ] ) with a slight modification @xmath150 . now , for example , in the limit @xmath151 , this goes like @xmath152 , which , multiplied by the prefactor @xmath153 , implies the behavior @xmath154 . that is , only the ( @xmath155)-term can give a non - vanishing contribution .      consider again eq . ( [ lapdifeq ] ) and a finite interval @xmath156 $ ] . the corresponding green function of the operator @xmath157 can be easily constructed ( at least formally ) with the help of the sturm  liouville theory  @xcite . an immediate consequence of the latter is that for @xmath158 the green function has the form @xmath159 where the functions @xmath120 and @xmath121 satisfy eq . ( [ diffeqfg ] ) with the initial conditions @xmath160 and @xmath161 , and @xmath162 and @xmath163 , respectively . in addition , the above green function should be symmetric due to the hermitian nature of @xmath164 .    the sturm  liouville theory ensures that the solution to the second - order differential equation ( [ diffeqfg ] ) is unique , when specifying the values of @xmath165 and @xmath166 at some point @xmath167 . therefore , the functions @xmath128 and @xmath129 must coincide with @xmath168 and @xmath169 of eq . ( [ diaggreenexpl ] ) , and the diagonal part of ( [ sturmgreen ] ) , i.e. @xmath170 , reduces to expression ( [ diaggreenexpl ] ) . this is an important consistency check of our representation ( [ diaglapl ] ) . let us now generalize the pi representation ( [ diaglapl ] ) to the full bloch density matrix , i.e. , we wish to include also the off - diagonal matrix elements , @xmath67 . if we go back to the replica representation ( [ afterrepl ] ) , we realize that the requirement @xmath67 spoils rotational symmetry in the replica field space , and thus precludes straightforward reduction to a radial path integral . instead of refining the reduction procedure , we simply make a guess , which , as we prove in appendix  [ sec : app ] , coincides with the well - established sturm  liouville formula ( [ sturmgreen ] ) . our guess is based on mathematical results presented in  @xcite . in particular , we claim that the extension of the representation ( [ diaglapl ] ) to off - diagonal matrix elements should read @xmath171 } \ , , \ ] ] with the action functional @xmath172 \ \equiv \ \mathcal{a}^e[\eta ] \ + \ \int_{x_-}^{x_+ } \!\ ! dx \ \ ! \eta^2(x)}\ ,   , \ ] ] where @xmath173 $ ] is defined in eq . ( [ actionho ] ) , and @xmath174 is a piecewise constant function @xmath175 \\    ( d-1)(d-3 ) & { \rm otherwise } . \end{cases}\ ] ] at this point we can invert the laplace transform with the help of eq . ( [ inverselap ] ) . as a result , we obtain the sought local - time pi representation of the bloch density matrix ( [ rhodef ] ) , namely @xmath176 } \ , .\ ] ] here , integrations over @xmath177 run from @xmath178 to @xmath179 , i.e. , the paths @xmath177 are non - negative . comparing this result with the anticipated heuristic form ( [ heurform ] ) , we can identify @xmath180 . representation ( [ offdiag ] ) allows us to identify the weight factor ( [ iii.7.a ] ) with @xmath181 \ = \ \frac{2}{d^2 } ( \eta_- \eta_+)^{\frac{1-d}{2 } } \eta(x_a ) \eta(x_b ) \exp \left\ { -\int_{x_-}^{x_+ } \!\!\ ! dx \left[\frac{\hbar^2}{2 m}{\eta'}^2 + \frac{m}{\hbar^2}\frac{\delta(x)}{8 \eta^2 } \right ] \right\ } .\ ] ] contrary to expectation , the right - hand - side of this expression does not depend on @xmath5 . sub - index @xmath89 in @xmath182 indicates that the weight factor must be regularized when we pull it out of the pi ( [ offdiag ] ) . by analogy with quantum mechanics one can represent ( [ offdiag ] ) in the discretized time - sliced form . in such a case the weight @xmath182 would be a product of terms involving the bessel functions @xmath113 , if @xmath183 , or @xmath184 , if @xmath185 ( see ref . @xcite ) . last but not least , expressions ( [ offdiag])-([weights ] ) indicate that the square root of @xmath3 is ( at least from physicist s point of view ) more convenient variable to describe local - time trajectories than @xmath3 alone . formula ( [ offdiag ] ) provides a way of rewriting the pi ( [ pix ] ) in terms of the local time . in this section , we consider more general scenario , in which the initial path integral is of the form @xmath186 \ \ ! \exp\left\{-\frac{1}{\hbar } \int_{0}^{\beta\hbar } \!\!\!d\tau \left[\frac{m}{2}\dot{x}^2 + v(x ) \right ] \right\ } , \ ] ] where @xmath128 is an arbitrary functional of the local time @xmath187 , which itself is ( as seen in section  [ sec : heuristic ] ) a functional of the paths @xmath20 . relation ( [ pifunc ] ) represents a local - time analog of the feynman  salam formula  @xcite . to bring it into more manageable form , we may observe that for any @xmath2 , the action of @xmath3 in the pi ( [ pifunc ] ) can be taken over by the functional derivative @xmath188 , acting on the exponential . this becomes transparent after rewriting the potential part as @xmath189 . the entire functional @xmath190 $ ] can be therefore pulled out of the path integral , which then allows to write @xmath191 \rho(x_a , x_b,\beta ) \ , .\ ] ] when we employ the local - time representation of pi for @xmath12 ( cf . ( [ offdiag ] ) ) , each functional derivative @xmath192 will produce the term @xmath193 . in such a way @xmath194 can be written as @xmath195 e^ { - \mathcal{a}^{e=0}_\delta[\eta ] } \ , , \label{ltfunc}\end{aligned}\ ] ] where , strictly speaking , the functional @xmath196 $ ] is regularized in such a way that it depends on @xmath3 only for @xmath197 $ ] , and @xmath198 are sent to @xmath199 only at the end of the calculation . first , let us make the simple observation that the formula ( [ ltfunc ] ) reduces to ( [ offdiag ] ) for the choice @xmath190 = 1 $ ] . one of the most important mean values of a local - time functional , as evaluated with eq . ( [ ltfunc ] ) , is the mean of @xmath200 which gives the moment - generating functional . the local - time moment structure is particularly pertinent in various perturbative expansions , including low- and high - temperature expansions ( see section  [ sec : asymp ] ) . another important example , namely the case of a one - point distribution function will be discussed in the following subsection .    in passing we should note , that should we have started from ( [ replrepre ] ) and repeated the above procedure , an analog of eq . ( [ ltfunc ] ) for higher - dimensional spaces , @xmath99 , could be easily obtained . this would include the @xmath89-dimensional replica field @xmath107 in the @xmath201-dimensional euclidean configuration space . simple , though quite important consequence of eq . ( [ ltfunc ] ) is that it readily provides the @xmath116-point distribution functions of the local time . this is achieved when we set @xmath190 = \prod_{n=1}^n \delta ( l^{x_n } - l_n ) $ ] . in order to see what is involved let us now illustrate the calculation for @xmath202 ( with @xmath203 ) . our discussion will be greatly simplified by considering only a free particle ( i.e. , @xmath204 ) that starts and ends at the origin , i.e. , @xmath205 . this corresponds to a stochastic process known as brownian bridge . our goal is to derive the one - point distribution function , denoted @xmath206 , of the local time at @xmath207 . we define @xmath206 by eq . ( [ pifunc ] ) with @xmath190=\delta ( l^0 - l ) $ ] , and calculate it from the representation ( [ ltfunc ] ) as follows .    in the laplace picture , @xmath208 , the path integral ( [ ltfunc ] ) can be sliced at @xmath209 so that @xmath210 the @xmath211-integration can be done easily by realizing that @xmath212 . furthermore , the limits in @xmath111 and @xmath89 can be carried out with the help of formulas ( [ etalimminus ] ) and ( [ etalimplus ] ) from appendix  [ sec : app ] . consequently , we obtain @xmath213 , \ ] ] where , for the free - particle case , @xmath214/\sqrt{2me/\hbar^2}$ ] , and @xmath215/\sqrt{2me/\hbar^2}$ ] , as one can straightforwardly verify . in the limit @xmath216 , ( [ ex1 ] ) reduces to @xmath217 and its inverse - laplace transform yields @xmath218    we stress that @xmath206 thus obtained is , in fact , the ( unnormalized ) _ joint _ probability density for stochastic events  @xmath219 " and  @xmath220 " . by bayes theorem of the probability calculus , the desired conditional probability density @xmath221 $ ] is obtained from ( [ ltfree ] ) by dividing @xmath206 by the brownian - bridge probability density @xmath222 $ ] , which is ( omitting again normalization ) @xmath223 ( see , e.g. , ref . normalization factors mutually cancel in the fraction and we arrive at @xmath224 \ = \ \frac{\hbar^2 l \exp\left(- \frac{l^2 \hbar^2}{2\beta m } \right)}{{\beta m}}\ ,   , \end{aligned}\ ] ] which is clearly normalized to @xmath225 . one can proceed along the same lines also in more complicated higher - dimensional ( @xmath226 ) cases . our result agrees with the one found through other means in ref . a compelling feature of the local - time representation ( [ offdiag ] ) is that it naturally captures both small- and large-@xmath5 asymptotic regimes . this should be compared with the feynman  kac pi representation ( [ pix ] ) , which is typically suitable only for the small-@xmath5 ( i.e. , large - temperature ) analysis . the latter is epitomized either by wkb approximation  @xcite or wigner  kirkwood expansion  @xcite . in the large-@xmath5 ( small - temperature ) limit , the spectral representation of the gibbs operator , @xmath227 , reduces the bloch density matrix to the ground - state contribution @xmath228 that is not evident from the feynman  kac pi representation  @xcite . in connection with eq . ( [ spectralasymp ] ) it is useful to remind that in @xmath229 the discrete bound states can all be chosen to be real  @xcite , so that the bloch density matrix is real and symmetric and can be written in the form @xmath230    let us first comment on the small-@xmath5 regime of the local - time representation , assuming @xmath66 for simplicity . this case was discussed in detail in our previous article  @xcite . there one should first taylor - expand the potential @xmath50 around the point @xmath33 , and then expand the exponential part containing the structure @xmath231 , where @xmath232^m\ , .\end{aligned}\ ] ] after the term @xmath233 is factored out of the integral , the individual summands of the ensuing series are of the form ( [ ltfunc ] ) with the potential @xmath204 , and functional @xmath190 \propto \prod_n l^{x_n}/\hbar$ ] . the latter can be related to the power expansion in @xmath5 presented in @xcite through the equality of representations ( [ pifunc ] ) and ( [ ltfunc ] ) . the whole bloch density matrix ( containing also off - diagonal elements ) can be treated similarly in a full analogy with ref . @xcite .    the large-@xmath5 expansion of eq . ( [ offdiag ] ) can be conveniently studied after rescaling @xmath234 , in which case we can write @xmath235 & \times & \ \delta \ ! \left ( \int_{x_-}^{x_+ } \!\!\ ! \eta^2 dx - 1 \right ) \exp \left\ { -\int_{x_-}^{x_+ } \!\!\ ! dx \left[\frac{\beta \hbar^2}{2 m}{\eta'}^2 + \beta v(x)\eta^2 + \frac{m}{\beta \hbar^2}\frac{\delta(x)}{8 \eta^2 } \right ] \right\ } \nonumber \\[3 mm ] & = & \ \lim_{d \rightarrow 0 } \frac{2}{d^2 } \lim_{\eta_\pm \rightarrow 0 } ( \eta_- \eta_+)^{\frac{1-d}{2 } }   \frac{\beta \delta^2}{\delta j(x_a ) \delta j(x_b ) } \int_{c - i\infty}^{c+i\infty } \frac{d \kappa } { 2\pi i } \int_{\eta(x_-)=\eta_-}^{\eta(x_+)=\eta_+ } \mathcal{d}\eta(x)\nonumber\\[1 mm ] & \times & \ \exp\left\ { - \beta \left[\int_{x_-}^{x_+ } \!\!\ ! dx \left(\frac{\hbar^2}{2 m}{\eta'}^2 + v(x)\eta^2 - \kappa \eta^2\right )          + \kappa \right ] \right\}\nonumber\\[1 mm ]   & \times & \ \left . \exp\left\ { - \int_{x_-}^{x_+ } \!\!\ ! dx \left [ \frac{m}{\beta \hbar^2}\frac{\delta(x)}{8 \eta^2 } + j\eta   \right ] \right\}\right|_{j=0 } , \end{aligned}\ ] ] where @xmath236 is an arbitrary real number . with the _ method of images _  @xcite we can rewrite the radial pi involved as a superposition of two genuine one - dimensional pi s  @xcite @xmath237 - \langle j|\eta \rangle \right\ } \exp\left [   - \int_{x_-}^{x_+ }   dx\ \ ! \frac{m}{\beta \hbar^2}\frac{\delta(x)}{8 \eta^2}\right]_{\!r}\nonumber \\[3 mm ] & & \mbox{\hspace{10mm}}= \ \int_{\eta(x_-)=\eta_-}^{\eta(x_+)=\eta_+ }   \mathcal{d}\eta(x ) \exp\left\{-\beta\left[\langle \eta| \hat{h}| \eta \rangle -\kappa\left(\langle \eta|\eta \rangle -1 \right )   \right ] - \langle j||\eta| \rangle \right\}\nonumber \\[1 mm ] & & \mbox{\hspace{14 mm } } \times \ \exp\left [   - \int_{x_-}^{x_+ }   dx\ \ ! \frac{m}{\beta \hbar^2}\frac{\delta(x)}{8 \eta^2}\right]\nonumber \\[1 mm ] & & \mbox{\hspace{10 mm } } -\   \cos(\pi d/2 ) \int_{\eta(x_-)=-\eta_-}^{\eta(x_+)=\eta_+ }   \mathcal{d}\eta(x ) \exp\left\{-\beta\left[\langle \eta| \hat{h}| \eta \rangle -\kappa\left(\langle \eta|\eta \rangle -1 \right )   \right ] - \langle j||\eta| \rangle \right\ } \nonumber \\[1 mm ] & & \mbox{\hspace{14 mm } } \times \ \exp\left [   - \int_{x_-}^{x_+ }   dx\ \ ! \frac{m}{\beta \hbar^2}\frac{\delta(x)}{8 \eta^2}\right ] , \label{method - images-1aa}\end{aligned}\ ] ] where dirac s notation was employed . a few comments are in order about the right - hand - side of the above relation . first , it should be noticed the presence of the parity - even terms @xmath238 in pi s . second , pi s differ by their respective dirichlet boundary conditions . third , the super - index `` @xmath239 '' was used to stress the restricted nature of the fluctuations in the radial pi measure while the measure without `` @xmath239 '' represents a usual one - dimensional pi measure , i.e. @xmath240 here `` @xmath241 '' denotes de witt s `` equivalence '' symbol  @xcite . finally , the correct time - sliced form of the exponential with the centrifugal potential is ( cf . e.g. , refs . @xcite ) @xmath242_{\!r } \ \lim_{n\rightarrow \infty}\prod_{k=1}^n \sqrt{2\pi \frac{\beta \hbar^2}{m } \frac{\eta_k\eta_{k-1}}{\varepsilon } \tilde{\delta}_k } \ \ ! \exp\left(-\frac{\beta \hbar^2}{m } \frac{\eta_k\eta_{k-1}}{\varepsilon } \tilde{\delta}_k \right ) i_{\frac{d-2}{2}}\left ( \frac{\beta \hbar^2}{m } \frac{\eta_k\eta_{k-1}}{\varepsilon } \tilde{\delta}_k\right ) , \nonumber \\[2 mm ]   & & \exp\left [   - \int_{x_-}^{x_+ }   dx\ \ ! \frac{m}{\beta \hbar^2}\frac{\delta(x)}{8 \eta^2}\right ] \ \lim_{n\rightarrow \infty}\prod_{k=1}^n \psi_{\frac{d-2}{2}}\left(-\frac{\beta \hbar^2}{m } \frac{|\eta_k\eta_{k-1}|}{\varepsilon } \tilde{\delta}_k\right ) ,   \label{centrifugal - slicing-1a}\end{aligned}\ ] ] with @xmath243 \\    1 & { \rm otherwise } , \end{cases}\end{aligned}\ ] ] and ( see , e.g. ,  refs . @xcite ) @xmath244 \ = \ e^x \sqrt{\frac{2 x } { \pi } } k_p(x)\ , . \end{aligned}\ ] ] [ @xmath245 and @xmath246 are the modified bessel functions of the first and the second kind , respectively . ] in cases when @xmath247 meaning that @xmath248 ( e.g. , `` typical situation '' for very fine time slicings ) the asymptotic form of @xmath249 holds . with help of the preceding asymptotic behavior one obtains @xmath250\nonumber \\[1 mm ] & \sim & \ \prod_{k=1}^n \exp\left[- \frac{m}{\beta \hbar^2 } \frac{\delta_k}{8|\eta_k\eta_{k-1}| } \ \!\varepsilon + { \mathcal{o}}(\varepsilon^2 ) \right]\nonumber \\[1 mm ] & \sim & \   \exp\left [   - \int_{x_-}^{x_+ }   dx\ \ ! \frac{m}{\beta \hbar^2}\frac{\delta(x)}{8 \eta^2}\right ] . \label{48.aa}\end{aligned}\ ] ] potential singularities of the integral at @xmath251 can be regularized , e.g. , by a principal value prescription . unfortunately , the formula ( [ 48.aa ] ) can not be directly used in our case because the boundary values @xmath252 and @xmath253 are arbitrarily close to zero , and hence the assumed asymptotic behavior for @xmath254 is not fulfilled . this situation can be rectified by factorizing out the problematic boundary points as @xmath255   & & \sim \ \frac{1}{\pi}\left(\frac{\beta \hbar^2}{2 m \varepsilon}\right)^{d-1 } \left(|\eta_-\eta_{1}||\eta_+ \eta_{_{n-1}}|\right)^{(d-1)/2}\ \ ! \prod_{k=2}^{n-1 } \exp\left[- \frac{m}{\beta \hbar^2 } \frac{\delta_k}{8|\eta_k\eta_{k-1}| } \ \!\varepsilon + { \mathcal{o}}(\varepsilon^2 ) \right ] . \ ] ] here we have utilized the asymptotic form @xmath256 valid for @xmath257 and @xmath258 ( @xmath259 ) . the passage from the radial pi ( [ offdiagrescaled ] ) to the ordinary ( 1-dimensional ) pi brings about an important advantage , namely , one can perform the wkb approximation . in particular , one can use the laplace s formula of the asymptotic calculus  @xcite @xmath260   \ = \   \sqrt{\frac{2\pi}{\beta g''(t_0 ) } } \ \ ! f(t_0 ) \ \ ! \exp\left[-\beta g(t_0 )    \right ] \ + \ { \mathcal{o}}\left ( \frac{\exp\left[-\beta g(t_0 )    \right]}{\beta^{3/2 } } \right ) , \label{laplace}\end{aligned}\ ] ] with @xmath261 being a solution of @xmath262 ( provided @xmath263 has a smooth absolute minimum at the interior point @xmath264 ) . in case of need , the full asymptotic expansion can be found , e.g. , in ref . @xcite . in our case the customary pi substitution @xmath265^{-1/2}\ , , \label{vi.49a}\end{aligned}\ ] ] should be employed  @xcite . here @xmath266 is the ground state energy and the difference operators ( lattice derivatives ) are defined as  @xcite @xmath267\ , , \;\;\;\;\ ;    \overline{\nabla } \eta(x ) \ = \ \frac{1}{\varepsilon}\left[\eta(x )   \ - \   \eta(x-\varepsilon ) \right]\ , , \end{aligned}\ ] ] with @xmath268 being the hermitian operator on the space of `` time - sliced '' functions with vanishing end points , i.e. , @xmath269 . the ensuing determinant @xmath270 in ( [ vi.49a ] ) can be conveniently computed by means of the gelfand  yaglom method  @xcite ( or shifting method  @xcite ) . there the determinant @xmath271 where @xmath272 is finite in the @xmath273 limit and fulfills the differential equation @xmath274 d(x ) \ = \ 0\ , , \label{vi.51a}\end{aligned}\ ] ] with the cauchy conditions @xmath275 and @xmath276 . the solution @xmath272 can be written as a linear combination of two independent solutions of ( [ vi.51a ] ) . one is immediate , namely the wkb solution @xmath277 . a second solution can be constructed with the help of dalambert s formula directly from @xmath277 in the form ( cf . e.g. , ref . @xcite ) @xmath278 it is simple to verify that ( [ vi.52a ] ) is again a solution of ( [ vi.51a ] ) . the cauchy conditions in the gelfand  yaglom method uniquely fix the constants in front of the two independent solutions so that the full solution @xmath272 acquires the form @xmath279 as a result we obtain in the large @xmath116 ( i.e. , small @xmath115 ) limit @xmath280 \ \stackrel{\varepsilon \rightarrow 0}{\longrightarrow } \   \frac{\eta_+ \eta_-}{\varepsilon }   \int^{x_+}_{x_- } \frac{dx}{\eta_0 ^ 2(x)}\ , . \label{vi.55aa}\end{aligned}\ ] ] the great advantage of the gelfand  yaglom method is that it does not require any detailed knowledge of the spectrum of the operator whose determinant is being computed . it also shows that the determinant in ( [ vi.49a ] ) diverges linearly with @xmath116 in the large @xmath116 limit . we now substitute for @xmath281 in ( [ laplace ] ) the functional expression @xmath282 \ , , \end{aligned}\ ] ] where , @xmath277 comes from the path that minimizes the functional @xmath283 . according to the rayleigh  ritz variation principle ( see , e.g. , refs . @xcite ) , such a function @xmath177 is the ground - state wavefunction of the hamiltonian @xmath164 , i.e. , @xmath284 with @xmath285 . with a hindsight , we have used the latter value of @xmath286 already in eqs . ( [ vi.49a ] ) , ( [ vi.51a ] ) and ( [ vi.55aa ] ) . notice that the stationary point in @xmath286 is real but the integration contour in @xmath286 is parallel to the imaginary axis . both reality and positivity of @xmath177 pose no restriction in the rayleigh  ritz principle , because the ground state can always be chosen real and positive  @xcite . putting it all together , we get the leading large-@xmath5 behavior in the form @xmath287 as it should be ( cf . ( [ spectralasymp2 ] ) ) . we conclude the discussion of low - temperature expansion by noting that the rayleigh  ritz variation principle states that _ all _ eigenvalues and ( normalized ) eigenvectors of @xmath164 come from stationary solutions of @xmath288 , and conversely  @xcite . in the spirit of the wkb approximation one should sum over all path integrals evaluated about _ all _ stationary solutions . it is , however , only the ground state configuration @xmath289 that acquires the global minimum and which gives the larges contribution to the wkb approximation . this fact was implicitly used in our preceding reasonings . should we have included also other stationary solutions we would recover higher order terms in the spectral expansion of the bloch density matrix ( [ spectralasymp2 ] ) . in this paper , we have derived the local - time pi representation of the bloch density matrix . we have shown that the result obtained , apart from being of interest in pure mathematics ( stochastic theory , sturm  liouville theory , etc ) , can serve as a useful alternative to the traditional feynman  kac pi representation of green functions of fokker  planck equations . furthermore , by analytically continuing the result back to the real time via the inverse wick rotation , @xmath290 , one obtains the local - time pi representation of quantum - mechanical transition amplitudes , i.e. , matrix elements of the evolution operator @xmath291 . from the physicists point of view , perhaps the most important application of local - time pi s lies in statistical physics , and namely in the low and high temperature treatments of the bloch density matrix . this is because in conventional pi s only a very tiny subset of paths gives a relevant contribution in these asymptotic regimes . in particular , the high - temperature regime of the boltzmann density function @xmath292 is dominated by paths that spend a sizable amount of time in the vicinity of the point @xmath0 . similarly , the low - temperature regime is controlled by paths with a large local time near the global minimum of the potential . here we have exemplified the conceptual convenience of the local - time formulation by providing a generic analysis of the low - temperature behavior of the bloch density matrix . our formulation proved to be particularly instrumental in obtaining the correct asymptotic behavior ( known from the spectral theory ) which is otherwise notoriously difficult to obtain within the feynman  kac pi framework  @xcite . as a byproduct we have uncovered an interesting connection between a low - temperature pi expansion , the gelfand  yaglom formula and the rayleigh  ritz variational principle . in order to further reinforce our analysis , we formulated a local - time analog of the feynman  matthews  salam formula which is ( similarly as its qft counterpart ) expedient in number of statistical - physics contexts . the prescription obtained was substantiated by an explicit calculation of a one - point distribution function of the local time . it appears worthwhile to stress that our local - time representation ( with its build - in replica field trick ) is in its present form applicable only to one - dimensional quantum mechanical systems . with a hindsight we reflected this fact already in our choice of the incipient pi ( [ pix ] ) where we assumed @xmath293 and @xmath294 . though one may easily proceed up to eq . ( [ replrepre ] ) without any restriction on the value of @xmath201 ( in fact , eq . ( [ replrepre ] ) is valid for any @xmath99 with @xmath295 ) , a further progress in this direction is hindered by the fact that the replica fields depend on a @xmath201-dimensional argument @xmath0 , and thus the pi in ( [ replrepre ] ) can no longer be regarded as a quantum mechanical pi ( i.e. , pi over fluctuating paths ) . in effect , we can not use existing mathematical techniques of the pi calculus ( e.g. , transformation of pi s to polar coordinates ) , that we have employed to get the radial pi  ( [ diaglapl ] ) . the issue of the extension of our local - time pi representation to higher - dimensional configuration space is currently under active investigation . the authors acknowledge a.  alastuey , h.  kleinert , z.  haba and d.s . grebenkov for fruitful discussions . the work has been supported by the gar grant no . ga14 - 07983s . was also partially supported by the ctu in prague grant no . sgs13/217/ohk4/3t/14 and by the dfg grant : kl 256/54 - 1 . in this appendix we show that the representation ( [ laploffdiag ] ) reduces to the well established result ( [ sturmgreen ] ) of the sturm  liouville theory . since @xmath83 is symmetric in @xmath33 and @xmath34 , we will assume , without loss of generality , that @xmath296 .    the path integral in ( [ laploffdiag ] ) can be expressed via eq . ( [ radialhosol ] ) as @xmath297 the limits in eq . ( [ laploffdiag ] ) can be carried with the help of the asymptotic formulas @xmath135 , and @xmath136 , valid for @xmath298 . we obtain @xmath299 & & \xrightarrow{d \rightarrow 0 } \ \frac{1}{\sqrt{\eta_a } } \exp \left ( -\frac{\hbar^2}{2 m } \frac{f'_1(x_a)}{f_1(x_a ) } \eta_a^2 \right)\ ! , \end{aligned}\ ] ] where @xmath140 satisfies eq . ( [ diffeqfg ] ) with initial conditions @xmath300 and @xmath301 , and similarly , we find @xmath302 where @xmath303 satisfies eq . ( [ diffeqfg ] ) with initial conditions @xmath304 and @xmath305 . formula ( [ laploffdiag ] ) then reduces to @xmath306 & & \times \ \ !   \exp \left [ - \frac{\hbar^2}{2 m } \left ( \frac{w(g_2,f_1 ) \eta_a^2}{f_1(x_a ) g_2(x_a ) } + \frac{w(g_3,f_2 ) \eta_b^2}{f_2(x_b ) g_3(x_b ) } \right ) \right]\ , , \end{aligned}\ ] ] where @xmath307 and @xmath143 satisfy eq . ( [ diffeqfg ] ) with initial conditions @xmath308 and @xmath309 , and @xmath310 and @xmath311 , respectively . the wronskian @xmath127 is independent of @xmath0 , as discussed in section [ sec : radial ] , and antisymmetric , i.e. , @xmath312 . wronskians @xmath313 and @xmath314 assume a particularly simple form when evaluated at points @xmath34 and @xmath33 , respectively , due to the initial conditions satisfied by @xmath169 and @xmath315 . we find @xmath316 and @xmath317 . rescaling @xmath318 , @xmath319 , and using the relation @xmath320 we obtain @xmath321 the integrations are readily performed using the formula @xcite @xmath322 yielding @xmath323 to prove equality with ( [ sturmgreen ] ) , we only have to show that @xmath324 this is done by realizing that @xmath143 , being a solution of the second - order linear differential equation ( [ diffeqfg ] ) , can be uniquely composed as a linear combination of two other solutions @xmath140 and @xmath303 , @xmath325 indeed , thus defined @xmath169 satisfies the initial conditions @xmath310 and @xmath311 . the low temperature expansion is typically obtained indirectly via duality approaches . among these , a particularly powerful non - perturbative approximation scheme is the so - called _ variational _ or _ optimized _ perturbation theory  @xcite . f.  steiner , _ path integrals in polar coordinates from ev to gev _ , in _ path integrals from mev to mev , bielefeld , 1985 _ , edited by m.c . gutzwiller et al . ( world scientific , singapore , 1986 ) pp . 335 - 359 .
Abstract:

Answer: * abstract * +    we derive a local - time path - integral representation for a generic one - dimensional time - independent system . in particular , we show how to rephrase the matrix elements of the bloch density matrix as a path integral over @xmath0-dependent local - time profiles . the latter quantify the time that the sample paths @xmath1 in the feynman path integral spend in the vicinity of an arbitrary point @xmath0 . generalization of the local - time representation that includes arbitrary functionals of the local time is also provided . we argue that the results obtained represent a powerful alternative to the traditional feynman  kac formula , particularly in the high and low temperature regimes . to illustrate this point , we apply our local - time representation to analyze the asymptotic behavior of the bloch density matrix at low temperatures . further salient issues , such as connections with the sturm  liouville theory and the rayleigh ritz variational principle are also discussed .