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 .


when a commensurate one - dimensional charge density wave ( cdw ) is doped , the resulting charges may be viewed as defects in the cdw . if the interactions between them are sufficiently strong compared to their effective kinetic energy , then they will form an ordered lattice ( which can be stabilized by a crossover to 3d order if there are many @xmath2 chains coupled together ) of defects and the cdw will become incommensurate with the lattice .
this can be more energetically favorable than simply doping holes into the rigid band formed in the presence of cdw order at a fixed wavevector .
the reason is that , in the latter case , we must pay the single - particle energy gap of the cdw . in a two - dimensional @xmath0-density wave ( ddw )
state , however , there are nodes in the order parameter , so it is relatively painless for doped holes to simply go to the nodes , which expand into fermi pockets ( in quasi-@xmath2 systems , however , this might not be possible , so doping must lead to incommensurability @xcite ) .
however , this can not continue indefinitely if the ddw state is stabilized by approximate nesting , since the fermi surface ( fs ) would eventually move away from the nesting wavevector .
thus , incommensurate ddw order is a strong possibility , at least over some range of hole dopings . in this paper
, we explore this possibility .
near half - filling , the hubbard model , the @xmath3 model , and generalizations of these appear to have many phases which are nearby in energy @xcite .
thus , small perturbations can strongly influence the competition between them .
consequently , it should not be too surprising that experiments on the cuprate superconductors have also uncovered evidence for a cornucopia of phases , particularly on the underdoped side of the phase diagram , which appear in particular materials and for certain values of the doping level , temperature , magnetic field , etc . @xcite .
depending on the interfacial energies between these phases , one way in which their competition can be resolved is through the formation of stripes @xcite or other inhomogeneous patterns of microscale phase - separation .
neutron scattering and scanning tunneling microscopy experiments point towards this possibility@xcite .
furthermore , as the cuprates are doped , their fermi surfaces evolve away from the commensurate nesting wavevector @xmath4 , so we would expect translational symmetry - breaking order parameters to occur at incommensurate wavevectors .
hence , it would seem important for any theory of the pseudogap to incorporate the tendency towards incommenuration , which we take as further impetus to study incommensurate ( ic ) order parameters related to ddw order .
in this paper we consider various ways in which @xmath0-density wave ( ddw ) order can go incommensurate . before turning to this discussion
, we review the reasons for expecting ddw order to occur in the pseudogap regime of the cuprates . it has been suggested@xcite that competition between ddw order and @xmath5-wave superconductivity ( dsc ) can explain many of the anomalous properties of the pseudogap@xcite .
the ddw order parameter is given by @xmath6 with @xmath7 and a _ commensurate _ ( c ) ordering wavevector @xmath8 .
the particle - hole condensate in breaks translation by one lattice spacing , parity , and time - reversal but is invariant under the combination of any two of the above symmetries . in real space
it is represented by currents which alternate from one plaquette to the next .
one motivation for considering ddw order is to explain @xcite the abrupt depletion of the superfluid density below the critical @xmath9 in bi2212 and similar dopings in other materials @xcite .
because the order has a @xmath5-wave momentum dependence , its onset would explain arpes measurements above @xmath10 in the pseudogap @xcite .
furthermore , since ddw is a spin - singlet condensate , it would also explain the suppression of the spin susceptibility in nmr measurements and opening of a spin gap in inelastic neutron scattering @xcite .
also , the depletion of the c - axis conductivity @xcite , whose contribution mostly comes from the antinodal regions , would also be naturally explained by the development of ddw order , which gaps these regions .
direct attempts to measure ddw order through neutron scattering have neither ruled it out @xcite nor unambiguously verified its presence , but there are intriguing suggestions that it may be present@xcite .
thus , it is natural to ask whether ddw order is compatible with incommensurability .
( for a discussion of the experimental signatures of ddw order , the reader is referred to the literature @xcite .
the case of ddw in @xmath11 is discussed in@xcite . for numerical evidence of the existence of ddw in other systems
, one is referred to@xcite . )
the plan of our paper is as follows . in section [ susc ]
we develop some preliminary intuition about the energetics of ic density wave order by analyzing the ddw susceptibility of a fermi liquid ( both the nested and non - nested cases ) . in section
[ c - ic ] we review some of the theory of commensurate ( c ) ddw order and show various ways to make it incommensurate .
we also suggest the kind of microscopic physics that could lead to c and ic ddw , and in addition discuss possible experimental signatures of ic ddw .
finally , in section [ num ] , we compute the mean - field phase diagrams of an extended hubbard model .
our main results are that ( a ) ic ddw develops for a wide range of parameters , but ( b ) the ic wavevectors can be very small .
their relevance to the cuprates may , at most , be for a rather narrow set of dopings .
we also offer a physical explanation of our results , which sheds some light on the nature of c - ic ddw transition in extended hubbard models .
we note that similar incommensurate order has been found in refs .
the easiest way to see that ic ddw phases could arise in effective models ( which , in this paper , will , in turn , emerge from microscopic extended hubbard models with correlated hopping terms@xcite ; see section [ micro ] for details ) is to consider the bare susceptibility as a function of the incommensurability wavevector @xmath12 , where @xmath13 : @xmath14 where @xmath15 is defined after eq .
[ cop ] . here , for simplicity
, ddw is considered without the presence of any other order . in eq .
[ suscep ] , @xmath16 is the fermi occupation number , and the energy dispersion contains nearest , @xmath17 , and next - nearest neighbor , @xmath18 , hopping parameters , @xmath19 .
the ddw equivalent of the stoner criterion is satisfied when the right - hand - side of this equation approaches the inverse of the pertinent coupling @xmath20 .
this observation allows us to interpret the positions of the peaks of @xmath21 as the wavevectors at which ddw order is likely to occur as we lower the temperature .
a similar analysis was done by schulz@xcite for the case of af and ic sdw order in the hubbard model close to half - filling , where perfect nesting at half - filling was asumed .
in fact ,     is at nonzero @xmath12 . here , @xmath22 is measured along the zone diagonal in units of @xmath23 , where @xmath24 is the lattice spacing.,height=326 ]    his results can be recovered by letting @xmath25 in if we set @xmath26 .
the support of @xmath27 in the integrand in is at @xmath28 and symmetry - related antinodal points of the fs since this is where the ddw order parameter is large .
however , since these points are also the location of van hove singularities , the regions around them also dominate the af susceptibility .
thus , at least in the @xmath26 case , one expects that whenever an ic sdw instability occurs , an ic instability for ddw should also occur .
indeed , for both orders , the peak of @xmath21 occurs at @xmath29 for low enough temperatures and sufficiently away from half - filling , e.g. if @xmath30 , as shown in fig .
[ suscfig ] .
note that the same analysis repeated with any other density wave or spin - density wave order would yield an ic instability as long as the dominant channels of the modulation coincide with fermi surface points which give large contributions to the susceptibility . for instance , a @xmath31-density wave@xcite with wavevector @xmath32 has a peak in @xmath33 for @xmath34 away from half filling at @xmath35 , in the case @xmath26 .
a natural question to ask is : to what extent does the ic ddw instability mimic the ic sdw one ? to answer this , one needs to consider specific models as we do in [ micro ] .
however , two differences can be anticipated even beforehand .
first , the ddw order parameter has nodes , while the af order parameter does not .
second , in a more realistic model , one must consider competition with dsc , which also could change features of the transition to ic state .
a better understanding of ic ddw order parameter is necessary to better understand the the situation away from half - filling .
given that the susceptibility points towards an ic instability , the next issue is what ic ddw order would look like . before addressing this question ,
let us first summarize some features of the c - case.@xcite in the c case , a general singlet particle - hole condensate can be written as ( neglecting the @xmath36 contribution ) : @xmath37 here , @xmath38 is the magnitude of the order parameter .
the wavevector @xmath39 is a reciprocal lattice vector .
the spin index @xmath40 included above will be omitted in the future for simplicity unless clarity requires otherwise .
in , the form factor @xmath15 , which can transform non - trivially under the point group of the 2-d lattice , is related to the angular momentum of the particle - hole condensate . for ddw order , @xmath41 which corresponds to @xmath42 and @xmath0 symmetry .
a @xmath31-density - wave with @xmath43 has @xmath44 , which corresponds to @xmath45 .    in real space , becomes @xmath46 where the vertices @xmath47 of the orders discussed above are given by @xmath48 note that the factors of @xmath49 in @xmath50 and @xmath51 signify that the corresponding phases break time - reversal , e.g. current flows along the bonds .
positive / negative signs in represent current going in / out of the vertex .
translation invariance , by default , is broken by any density wave order since there is a preferred vector @xmath52 , but in the c case , invariance under a _ subset _ of the lattice group is retained .
ddw order , for example , is invariant under translation by linear combinations of even reciprocal lattice vectors along each bond direction .
the current along the bonds is given by : @xmath53 here @xmath54 is along either the @xmath55- or @xmath56-direction . because of the equivalence of @xmath52 and @xmath57 in the c case , the current is simply @xmath58 .
in the ddw case this results in an alternating plaquette current / checkerboard magnetic field patterns.@xcite      to obtain an ic version of ddw order , one first notes that ddw is chiral and hence breaks the ising symmetry which reverses all the clockwise plaquette currents into counter - clockwise and vice versa .
therefore one can construct anti - phase domain walls as in fig.[fig : domains ] . a simple caricature of a bond - oriented domain wall is given in fig.[fig : domains]a .
its real space representation is @xmath59 where @xmath60 is the step function , so @xmath61 with @xmath62 . at the midpoint of the domain wall ,
the currents in the @xmath63-direction vanish while the currents in the @xmath64-direction are halved .
the diagonal domain wall depicted in figure [ fig : domains]b has a similar representation .
note that because of the halving of the @xmath64-direction currents at the midpoint of the bond - domain wall in , we need to include the term : @xmath65 where the new vertices @xmath66 are given by : @xmath67 such vertices seem necessary for ic transitions of current density waves .
one could imagine situations in which the transition to a disordered state is driven by thermal proliferation of such vertices , similar to earlier analysis of the 6-vertex model@xcite .
although the analysis of such a transition is an interesting topic on its own , we will not pursue it in this paper .
we are more interested in the energetics of the ic phase , so a long - wavelength mean - field treatment is sufficient .    to analyze the c - ic transition ,
one then considers an array of alternating domain walls with average separation @xmath68 .
then , in the single harmonic approximation one retains only the @xmath69 term in and identifies @xmath70 with the incommensurability wavevector .
the new short - range vertex contributions are also neglected in the lowest harmonic approximation .
the neglect of higher harmonics should not change the transition temperature@xcite or the value of the critical doping for the onset of ic order since @xmath38 is small and the corrections are of higher order in @xmath38 .
an incommensurate instability in a clean system occurs when @xmath68 becomes finite .
( in a disordered system , it occurs when when @xmath71 , where @xmath72 is a characteristic length scale for the disorder . )    within the single harmonic approximation , an array of anti - phase domain walls with separation @xmath73 , looks simpler in momentum space : @xmath74 note that in the limit @xmath75 , the above reduces to .
one could have guessed the result by writing @xmath76 in the commensurate case , and then let @xmath52 go incommensurate .
this suggests other ways to get ic order .
in particular the oddness of ddw order with respect to rotations by @xmath77 and transpositions about , say , the @xmath4 direction lead to the following example of a ddw _ checkerboard pattern _ : @xmath78 where the operation @xmath79 can be either a transposition , or a rotaion by @xmath77 . the checkerboard pattern in corresponds to simply superimposing two domain walls rotated by @xmath77 with respect to each other .    in a similar spirit
, one could also investigate the possibility of _ non - topological _
domain walls@xcite , e.g. ic domains as the above , superimposed with an uniform order , so the overall order parameter does not change chirality . however , as we will see in sec .
[ num ] , much of the physics of the c - ic transition will already be transparent from the simpler choices and of order parameter .
finally , note that within the single harmonic approximation , and are current - conserving only to lowest order in @xmath22 .
for the purpose of clarifying the energetics of the ic order close to the c - ic transition , or in the cases where the magnitude of the incommensurability is small , this lowest order approximation is good enough .
our numerical results , sec .
[ num ] , _ a posteriori _ show that we are in such regime
.      we will consider the above order parameters in the context of the following hamiltonian @xcite : @xmath80 in this formula , @xmath81 is hopping with @xmath82 for nearest neighbors , @xmath83 for next nearest neighbors and @xmath84 otherwise .
@xmath85 is a correlated hopping term which simultaneously hops an electron from site @xmath86 to site @xmath49 and hops an electron from @xmath87 into the vacated site @xmath86 .
the on - site and nearest - neighbor repulsions are , respectively , @xmath88 and @xmath89 . the indices @xmath90 signify lattice sites and @xmath91 the spin .
we will consider this hamiltonian at the mean - field level . at this level ,
the energetics of ddw order is the same as for a ` bcs - reduced ' hamiltonian for ddw order of the form : @xmath92 with @xmath93 .
a similar reduced hamiltonian for @xmath5-wave superconductivity has @xmath94 .
note that the processes responsible for both dsc and ddw are essentially kinetic .
a hamiltonian of the form ( [ micro - ham ] ) has one important virtue : it has a @xmath5-wave superconducting ground state over a range of dopings ( and an antiferromagnetic ground state at half - filling ) .
since this is the _ sine qua non _ for any description of the cuprates , we believe that this is a good starting point . one may object that correlated hopping terms have to be put in by hand in a microscopic hamiltonian .
however , they naturally arise even from the one - band hubbard model away from half - filling , where their contribution can be computed in @xmath95 perturbation theory : @xmath96 at dopings @xmath97 , for @xmath98 , we would have @xmath99 . thus ,
not far away from half - filling , correlated hopping terms are certainly not negligible .
let us generalize the mean field description of c ddw to the ic case .
the path we will take is similar to that in sec .
[ icop ] , where we wrote down the _ commensurate _ order parameter without making use of the @xmath100 symmetry , and as a result the generalization to incommensurate @xmath52 was straightforward .
we thus obtain a functional @xmath101 , equation , whose minimization yields self - consistently both the order parameters @xmath102 and the deviation from commensurability @xmath103 .
the only subtlety , as we discuss at the end of this section , is that although @xmath104 reduces to the mean - field free energy in the c phase , in the ic phase it is _ not _ the free energy even within the single harmonic approximation . instead , minimization of @xmath104 is equivalent to a variational search of the ground state .
our derivations generalize the treatment of nayak and pivovarov @xcite of competition of c ddw with dsc and possibly af .
the c ddw mean field hamiltonian we start with is : @xmath105 where @xmath106 .
note that in the @xmath107 symmetry for _ incommensurate _
@xmath52 is not equivalent to the hermiticity of the hamiltonian due to the @xmath108-dependance of @xmath109 ( e.g. @xmath110 for ic * q * ) .
this is unlike the case of af order , where the interaction is simply @xmath111 and we have @xmath112 .    in order to check with earlier results on the c - ic transition of af order in the _ weak coupling _
hubbard model@xcite ( see sec . [ susc ] ) , we will also include a term like in our full reduced hamiltonian .
of course , in the strong coupling , @xmath113 , limit , one is not justified in treating af as a fermi surface instability , especially close to half filling .
a @xmath3 description is more appropriate in this case .
however , as mentioned in sec .
[ micro ] , both ddw and dsc are essentially kinetic in origin within the extended hubbard model , so even in the strong onsite repulsion limit , there is a good justification to use mean field theory to describe their competition away from half filling .    in the c case ,
the standard dsc term in our reduced hamiltonian is : @xmath114 where @xmath115 . using a nambu basis @xmath116
, the full hamiltonian @xmath117 can be rewritten as @xmath118 where : @xmath119 with @xmath120 .
the action derived from the above hamiltonian yields , upon integrating the fermions out , the free energy in the ordered state : @xmath121 with @xmath122 the energy eigenvalues of , and @xmath123 if we assume @xmath100 symmetry , @xmath124 would be given by : @xmath125 ^ 2}\ ] ] where @xmath126 and @xmath127 .
however , in the ic case , @xmath100 need not hold , in which case the above would generalize to :    @xmath128\delta^2_\text{dsc}+ f^4_-({\bf k})\delta^4_\text{dsc}. \end{split}\ ] ]    the effect of a nonzero incommensurability @xmath12 can be seen not only in the appearance of the terms @xmath129 and @xmath130 , but also in the dsc - and doping - dependent function @xmath131 .
our final step before letting @xmath52 be incommensurate is to pass the integration in from rbz to bz : @xmath132 once we convert the rbz integrals to bz integrals as in , we numerically minimize ( see sec .
[ num ] ) the resulting ` free energy ' with respect to both the order parameters @xmath102 and the deviation from commensurability @xmath12 .
ev and @xmath133ev are large ( compare to figs .
[ fig : ddwpd3 ] and [ fig : ddwpd4 ] ) @xmath134ev .
the next nearest neighbor hopping , @xmath135ev , is small compared to typical values fitted to arpes@xcite , but ic order was nevertheless robust away from half - filling .
inset : ic wavevector @xmath22 at @xmath136ev as a function of doping .
wavevectors are in units of @xmath137 , where @xmath24 is the lattice spacing.,height=288 ]    it must be emphasized that the expression to be minimized in sec .
[ num ] is _ not _ the true free energy of the system if @xmath52 is incommensurate .
one way to see this is to note that the eigenvalues are not invariant under @xmath138 if @xmath139 and therefore they can not be interpreted as band energies . in a c - ic transition , the excitations split into a generally uncountably many bands even within the single harmonic approximation@xcite .
therefore , the problem of finding the mean - field free energy is notoriously difficult even in the single harmonic limit .
the calculation we are doing , on the other hand , is essentially a mean - field variational search for the ground state , and not for the excitations .
a similar calculation was done by schulz@xcite for the case of the c - ic transition of the weak - coupling hubbard model close to half - filling . by continuity of the true free energy@xcite
, one expects that the resulting variational estimates for @xmath140 would be accurate only close to the c - ic transition , where @xmath22 is small .
ev , @xmath141ev are ` large ' ( note temperature scale ) similar to fig . [
fig : ddwpd1 ] , with @xmath142ev .
however this time we taka a more realistic @xmath143ev . again
ic order is robust at dopings relevant for the pseudogap .
inset : ic wavevector @xmath144 at @xmath136ev as a function of doping .
[ fig : qdw4 - 2tp012 ] for a more detailed view of @xmath144 ( and also fig .
[ fig : qdw4 - 2tp006 ] for the case of @xmath145ev ) .,height=288 ]    ev , @xmath146ev , @xmath147ev .
the next nearest neighbor hopping , @xmath148ev .
note the appearance of dsc proximate to ` magic ' dopings .
inset : the ic wavevector @xmath149 at @xmath150ev .
see fig.[fig : qdw25 - 1tp012 ] for comparison between @xmath149 and @xmath151 .
, height=288 ]    to obtain the phase diagrams for the models in , for various fixed dopings @xmath63 and temperatures @xmath152 , we minimize @xmath153 with respect to both the order parameters @xmath102 ( as was done previously in the c - case@xcite ) and the ic wavevector magnitude @xmath22 .
each choice of ic order parameter ( e.g. bond or diagonal ) corresponds to a separate minimization procedure . to distinguish which ic order dominates at a particular point in @xmath154 space , we compared the free energies of the orders considered , and chose the ic phase of minimal energy .
note that , as given in , the free energy is a function of @xmath155 , so to obtain the dependance on the doping @xmath63 , we first perform a legendre transformation and minimize instead @xmath156 where the chemical potential is obtained from @xmath157 with @xmath158 .
our approach is slightly different than that previously used in the c - case@xcite .
there , the competition between af , ddw and dsc was analyzed , and some sensitivity on the details of the microscopics was found . phase diagrams resembling that of the cuprates were obtained .
the general trend was that any phase diagram which included @xmath5-wave superconductivity also had ddw order competing with it . by changing the bare coupling constants ,
rather different phase diagrams could be obtained and these generally did not have superconductivity .    in our treatment here , we do not consider the competition with af order , which is not a fermi surface instability anyway in the context of the cuprate superconductors .
in , we still include af , but only in order to check with previous results@xcite of ic af order away from half filling in the _ weak coupling _ @xmath159 limit . on the other hand ,
the af order appearing at half filling in the cuprates is a strong coupling phenomenon .
a reduced hamiltonian is the wrong starting point at half filling in the limit @xmath113 .
therefore , we only consider competition between ddw and dsc and look for ic order in this case .
all the phase diagrams constructed this way are robust in the sense that ddw and dsc are found to coexist at mean field level for almost all reasonable values of the coupling constants .
sample phase diagrams are shown in figs .
[ fig : ddwpd1]-[fig : ddwpd4 ] . in all of them , ic order is present at dopings @xmath160 where ddw and dsc compete . in the first two of these figures ,
we vary @xmath18 at relatively high @xmath161 and @xmath162 , while in the last two we explore the phase diagram for more realistic ddw and dsc couplings ( compare the temperature scales to those in the cuprates ) . note the complexity of figs .
[ fig : ddwpd3 ] , figs . [ fig : ddwpd4 ] compared to figs . [
fig : ddwpd1 ] , figs .
[ fig : ddwpd2 ] . at realistic couplings
, the ic wavevector profile is complicated , and furthermore dsc order appears in more than one doping region within the ddw dome .
in fact , the appearance of incommensurability seems correlated with the appearance of dsc . a more systematic analysis to be discussed below reveals that generally the types of c - ic transitions that were produced in our phase diagrams fall into two categories : ( a ) those which correspond to motion of the bare fermi surface away from commensurate nesting ( shaded regions ) ;     except @xmath163ev .
there is a distinct ic response ( inset ) at @xmath164 and @xmath165 .
[ fig : qdw2 - 1tp012 ] for comparison between @xmath151 and @xmath149.,height=288 ]     where ddw and dsc compete .
note that nonzero @xmath22 is found above as well as below the dsc critical temperature.,height=288 ]
\(b ) those which appear near certain ` magic ' filling fractions . case ( a ) is similar to the familiar c - ic mechanism discussed by schulz@xcite , while case ( b ) appears to be a consequence of negative effective stiffness at special fractions , the exact nature of which at present is not clear . in both cases
@xmath22 is independent of the onset of dsc , fig .
[ fig : qoft ] , but because incommensurability relieves the ddw - dsc competition , the onset of dsc correlates to the presence of ic ddw .
let us first discuss ( a ) .
it is well known@xcite that in the weak - coupling hubbard model at @xmath26 , the af instability at half - filling is driven ic for any non - zero doping ( see sec .
[ suscep ] ) .
the source of the c - ic transition can be traced to the approximate nesting of the bare fermi surface ( fs ) at ic wavevectors away from half - filling .
close to half - filling , this picture would produce a linear deviation from commensurability , @xmath166 , fig .
[ fig : afpd1 ] .
the incommensurability ( deviation from @xmath4 ) is horizontal or vertical , not diagonal , as we would expect because the bare susceptibility is more divergent due to a greater overlap between nested portions of the fermi surface@xcite .
a similar scenario happens in the case of ddw order .
however , two notable differences from the af case are present .
first , ddw is less sensitive to @xmath167 , since it is driven mainly by antinodal regions of high density of states ( d.o.s . ) .
therefore , while there is no af instability without nesting ( at half - filling af fermi surface instability is only possible if @xmath168 ) , ddw is still possible in this regime .
this suggests that whenever the antinodal regions of the bare fs evolve away from the anti - nodal regons ( which happen to also be regions with a high density of states as a result of the presence of van hove singularities ) , a c - ic transition is likely to occur .    .
the results of schulz@xcite were recovered . here , we used @xmath169ev with @xmath134ev .
bond - oriented ic order was found to be energetically favorable , as a direct computation@xcite in the weak repulsion regime would yield .
_ inset _ : as expected , the ic wavevector at @xmath170 is linear in @xmath63 .
, height=288 ]    the second difference from the af case is that ddw has nodes .
this implies that the chemical potential in the ordered state need not be higher than the maximal ddw gap . as a result
, ddw order can stay commensurate for wider ranges of doping .
this is indeed what we find .
for example , for all @xmath171 , we obtain @xmath172 near zero temperature . in general
, we find that the doping values @xmath173 for which the above type of c - ic ddw is possible are unrealistically high in the context of the cuprates , if one considers more realistic bare fs , @xmath174 .
furthermore , in the ic ddw case , there need not be a simple relation @xmath166 between the doping and doping and deviation from commensurability .
let us now discuss the c - ic transitions that are present in our phase diagrams proximate to ` magic ' filling fractions @xmath175 .
although our understanding as to why they arise is incomplete , we would like to point out several interesting features about our results .
first , to our knowledge , there is nothing special about the bare fs close to these filling fractions , likewise for the fs of the commensurate ordered state .
the c - ic transition is not driven by nesting , since upon variation of @xmath18 we see similar signatures proximate to the ` magic ' dopings .
second , the magnitude of the ic wavevector is generally small ( corresponding to periodicity of more than @xmath176 lattice spacings ) and shows no clear dependance on @xmath177 , figs .
[ fig : qdw4 - 2tp006 ] -[fig : qdw2 - 1tp012 ] .
third , we checked that similar behavior is present deep into the af ordered state in the total free energy .
[ fig : qaf004t0tp004 ] . unlike the ddw case , however , the ic lattice spacing in the af case depends on @xmath177 : @xmath22 follows an increasing staircase pattern at ` magic ' dopings .
note that our results for af order seem to contradict schulz s results@xcite which predict a linear dependence of @xmath149 .
the naive contradiction is easily resolved by realizing that schulz s computation@xcite of @xmath149 relies on the stoner criterion , which is strictly valid at @xmath10 .
indeed , while deep into the ordered state the ic wavevector is stepwise increasing at ` magic ' fractions , at higher temperatures we regain the linearity of @xmath178 , fig .
[ fig : qaf004t035tp004 ] .
moreover , the fact that both the bare and the gapped fs in the c case are generally not nested at the the ` magic ' fractions suggests that the transition is not driven by a divergent susceptibility , but is first order .
indeed , at low temperatures , the chemical potential near the special dopings is non - monotonic for both the af and the ddw+dsc diagrams .
similar behavior of @xmath151 was found by nayak and pivovarov in the context of c ddw competing with dsc@xcite .
a simple view of the energetics is given by the following landau - ginsburg free energy of two nodeless competing orders ( the generalization to order parameters with nodes is straightforward ) : @xmath179 in this free energy , @xmath180 plays the role of the ddw order parameter and @xmath181 the role of superconductivity . when @xmath182 becomes negative , @xmath180 orders at non - zero wavevector @xmath183 .
this occurs at the dopings indicated in figures [ fig : ddwpd3 ] and [ fig : ddwpd4 ] .
when @xmath181 orders , this further favors incommensurability if @xmath184 since the competition between the two orders is partially alleviated .
this alleviation is a possible explanation of why the onset of dsc at low @xmath152 is correlated to the presence of ic - ddw .     and @xmath151 ( inset ) at @xmath136 for @xmath185ev , @xmath186ev .
note that @xmath187ev which is smaller than the value taken in fig .
[ fig : ddwpd2].,height=288 ]     except this time @xmath143ev as in fig .
[ fig : ddwpd2].,height=288 ]     and @xmath149 ( inset ) at @xmath150ev corresponding to the parameters in fig .
[ fig : ddwpd3 ] .
although here the ic response is not at ` magic ' fractions , note that in fig .
[ fig : ddwpd3 ] there is dsc both at @xmath164 and @xmath165 . ,
height=288 ]     ( and at the same temperature ) except this time there is a clear c - ic transition at ` magic ' fractions ( inset ) . at the same time
@xmath151 is non - monotonic at these dopings .
the coupling constants @xmath188ev and @xmath189ev are the same as in fig . [
fig : ddwpd4].,height=288 ]     ev ) staircase ic spectrum of af around ` magic ' filling fractions .
the parameters we have taken here are @xmath190 ev , @xmath142 ev .
inset : @xmath151 is non - monotonic at both jumps of @xmath149 , but is smoothly decreasing whenever @xmath22 is pinned.,height=288 ]    the same phenomenon can be viewed from the opposite perspective : when ddw order becomes incommensurate , there are regions in which it is suppressed ( as in the domain walls of fig .
2 ) . superconductivity can occur more easily there .
hence , incommensurability alleviates the competition between dsc and ddw order at the mean - field level , thereby enhancing superconductivity .
of course , ic order appears to suppress superconductivity in the cuprates , but this could be due to effects beyond a simple mean - field theory .    when the superconducting order parameter becomes strong enough , it may become energetically favorable for phase separation to occur .
we can guess when this occurs from the usual maxwell construction whenever the mean field chemical potential @xmath151 does not vary monotonically .
this appears to occur at a few isolated ` magic ' doping values . however , in a real system , proximity to a second order critical point@xcite , coulomb repulsion@xcite , or disorder could drive the system in a state of mesoscopic phase separation .
moreover , the lattice pins any ic density wave order , so once the system is ic in order to minimize the free energy , the deviation from commensurability @xmath22 would be pinned at a finite value . finally , it is important to remember that disorder could have important effects .
it would not only disorder the array of domain walls , but could lead to droplets of one phase in regions of another
.     regains its linearity around @xmath191 as predicted by schulz@xcite . here
@xmath192ev and all other parameters are as in fig .
[ fig : qaf004t0tp004 ] .
inset : @xmath151 is monotonic throughout , suggesting that at higher temperatures , the c - ic transition is driven purely by ic nesting of the fs away from half - filling.,height=288 ]    we stress that while the general trend towards ic order may be clear , many of the quantitative aspects are not clear , in particular the precise copings at which ddw order becomes incommensurate .
it is not clear to us why such c - ic transitions happen near ` magic ' filling fractions .
the fact that similar behavior occurs in the case of af order suggests that the tendency of phase separation at ` magic ' filling fractions is a general trend for a variety of density wave order . because the af only case was explored upon taking the @xmath193 limit in the full free energy , it is not clear yet whether the tendency of phase separation at ` magic ' filling fractions is purely a feature of density waves competing with other order parameters , or of the lattice
. it might be worth reconsidering the af case in the presence of other order parameters ( such as ddw , dsc ) .
what makes our results surprising is that unlike in the frustrated phase separation picture for stripe formation , there is no coulomb repulsion present .
hence , there is no obvious connection between phase separation and incommensurability .
( please note that while ic ddw order can induce charge density wave order at twice the wavevector , we have not included this effect , which we expect to be small , in our calculation . ) therefore , while coulomb repulsion may play a role in determining the cdw spacing at @xmath164 in @xmath194 , in our calculation , the phase separation at ` magic ' fractions can only arise as interplay between the lattice , competing order physics and energetics .
what exactly is the combination of the three deserves some future research as the analog of the fulde - ferrel - larkin - ovchinnikov ( fflo ) state of dsc in a magnetic field parallel to the plane@xcite . in the absence of next - nearest neighbor hopping
, the two mean fields are related via a particle - hole transformation in the same way _ commensurate _ ddw is related to dsc at half - filling . however , while in the latter case one can rotate c ddw to dsc@xcite , in the case of finite @xmath155 , the particle - hole transformation can not interchange fflo and ic ddw .
while the mean - field hamiltonians of ic ddw and the fflo are related by a particle - hole transformation , the spectra are not  in the ic ddw state there is a splitting to a generally uncountably many bands ( if @xmath195 is irrational as we discussed in sec .
[ num ] ) since the order parameter connects @xmath108 and @xmath196 while in the case of dsc there is no such splitting .
we are indebted to the authors of @xcite for pointing out their work to us . ] .
according to the discussion following it seems likely that within the dsc state a competing ddw order parameter is ic as it develops .
hence , we briefly describe the experimental signatures that would distinguish between c and ic ddw order .
originally , c ddw was also named ` hidden order ' in the context of the pseudogap , because most of the experimental signatures associated with it are indirect .
the magnetic field , created by the alternating currents around neighboring plaquettes , is too small to be detected by neutron scattering ( @xmath197 g for @xmath198 mev ) . at the same time , because of the symmetry of the c ddw magnetic field , there is no signature in nmr .
the charge of the c ddw state is uniform so stm would not show signatures , unless they are associated with interference of nodal quasiparticles due to scattering off impurities . it has been argued recently@xcite that such interference disperses weakly enough to account for the pseudogap d.o.s .
incommensurate modulations observed by stm . for a further discussion of the signatures of c ddw , the reader is referred to the literature @xcite .
an ic ddw is more easily detectable directly .
the local electronic density of states would be modulated due to the ic order parameter , and hence could be observable by stm .
secondly , the ic modulations of the staggered currents would produce nmr line splitting of both the cu and o atoms when the applied magnetic fields are _ perpendicular _ to the a - b planes .
this is in sharp contrast with ic sdw , where there is an nmr line - splitting due to spatially modulated local magnetization , which is directed _ along _ the a - b plane .
let us estimate the modulated magnetic field strength arising from ic ddw ( the reader is referred to@xcite where similar considerations are applied to the case of orbital antiferromagnetism proposed to exist in @xmath199 ) .
suppose we have an array of domain walls , fig .
[ fig : domains ] , along the y - axis .
deep between the domains the field is zero as in the case of c ddw . along a domain at @xmath200 , where the field is maximal , again most of the fields cancel except those created by the vertical alternating currents of strength @xmath201 at @xmath202 and those of strength @xmath203 at @xmath204 .
equivalently this would be the field of a single vertical array of alternating currents of strength @xmath203 at @xmath202 : @xmath205 where @xmath206/[\sqrt{((i-1)a + y_0)^2 + x^2_0}]$ ] and @xmath207 for o , cu atoms respectively . clearly at @xmath208 , cu ,
atoms the field is zero by symmetry . for o atoms
the sum can be performed yielding @xmath209 . in the dilute domain limit , this field is independent of the distance between the domains .    for a sinusoidal current modulation , @xmath210 would depend on the modulation length @xmath211 . if @xmath212 then a good estimate of @xmath210 is given by with @xmath201 replaced by @xmath213 where @xmath177 is the periodicity of the modulation ( in units of the lattice spacing ) .
therefore , for large @xmath177 we obtain @xmath214 .
incommensurate order is an important possibility when competing orders are considered . at a first - order transition
the system is phase separated and hence necessarily inhomogeneous .
if the competition between the orders is weak , coexistence is possible , but modulations of the competing orders are likely to be induced .
the effects of the lattice , disorder , lomg - range coulomb repulsion@xcite and proximity to criticality@xcite is to generally stabilize incommensurability on mesoscopic length scales .
this is unlike the case of , say , the liquid to vapor phase transition where the competition results in phase separation on macroscopic length scales .    in the case of ddw and dsc ,
the competition is certainly weak within the framework of the extended hubbard models we have considered . because the processes stabilizing both orders are of similar nature and are at the same time smaller than all of onsite repulsion , exchange and hopping , it does not cost much energy to convert ddw order into dsc and vice versa .
a simple resolution of the competition is for ddw order to become incommensurate .
we find that this occurs for some doping levels , but the incommensuration tends to be relatively small .
such ic ddw could , in principle , be observable at low temperatures by nmr since it creates an inhomogeneous magnetic field at o atoms _ perpendicular _ to the ab - plane .
most of nmr experiments test for much stronger magnetic fields _ parallel _ to the ab - plane on the cu atoms , and at higher temperatures such ic ddw phase would be smeared by thermal fluctuations , which may exlpain why such ic order , if it exists , it hasnt been seen by experiment .
we would like to thank s. brown , s. chakravarty , r. jamei , s. kivelson , and e. pivovarov for helpful discussions .
special thanks to e. pivovarov for letting i. d. generalize some of his computer code .
this work has been supported by the nsf under grant no .
dmr-0411800 .
99 u. schollwck , s. chakravarty , j. o. fjaerestad , j. b. marston , and m. troyer , phys .
lett . * 90 * , 186401 ( 2003 ) ; e. dagotto , rev .
phys . * 66 * , 763 ( 1994 ) ; c. s. hellberg and e. manousakis , phys .
lett . * 83 * , 132 ( 1999 ) ; s. sorella , g.b .
martins , f. becca , c. gazza , l. capriotti , a. parola , e. dagotto , phys .
lett . * 88 * , 117002 ( 2002 ) .
s.  r. white and d.  j. scalapino , phys .
b * 61 * , 6320 ( 2000 ) ; d. f. b. ten haaf , p. w. brouwer , p. j. h. denteneer , and j. m. j. w. o. putikka , m. u. luchini , and r. r. p. singh , phys .
lett . * 81 * , 2966 ( 1998 ) ; l.  p. pryadko , s.  a. kivelson , o. zachar , phys . rev . lett .
* 92 * , 067002 ( 2004 ) ; e. dagotto and t. m. rice , science 271 , 618 - 623 ( 1996 ) ; g. roux , s. r. white , s. capponi , a. laeuchli , d. poilblanc , phys .
b 72 , 014523 ( 2005 ) .
experiments in the cuprates have been interpreted in terms of competing orders discussed in refs . ; and also in c.  m. varma , phys .
b 55 , 14554 ( 1997 ) ; s.  c. zhang , science 275 , 1089 ( 1997 ) .
e. w. carlson , v. j. emery , s. a. kivelson , and d. orgad , in _ the physics of conventional and unconventional superconductors _ ,
k. h. bennemann and j. b. ketterson ( springer - verlag 2003 ) , cond - mat/0206217 ; a.  h. castro neto and c. morais smith , in _ strong interactions in low dimensions _ , ed .
d. baeriswyl and l. degiorgi ( kluwer acad . pub .
2004 ) , cond - mat/0304094 ; v.  j. emery , s.  a. kivelson , and j.  m. tranquada , proc .
usa * 96 * , 8814 ( 1999 ) ; m. vojta and s. sachdev , phys .
lett . * 83 * , 3916 - 3919 ( 1999 ) ; j. zaanen , nature 404 , 714 ( 2000 ) ; and references therein . j.  m. tranquada , physica b 241 - 243 , 745 ( 1998 ) ; in _ neutron scattering in layered copper - oxide superconductors _ , ed.a . furrer ( kluwer , dordrecht , the netherlands ) , p. 225 . j. e. hoffman , k. mcelroy , d .- h .
lee , k. m lang , h. eisaki , s. uchida , and j. c. davis , science * 297 * , 1148 ( 2002 ) ; k. mcelroy , r. w. simmonds , j. e. hoffman , d .- h .
lee , j. orenstein , h. eisaki , s. uchida , j. c. davis , nature * 422 * , 592 ( 2003 ) ; c. howald , h. eisaki , n. kaneko , m. greven , and a. kapitulnik , phys .
b * 67 * , 014533 ( 2003 ) .
s. chakravarty , r. b. laughlin , d. k. morr , and c. nayak , phys . rev . *
b 63 * , 094503 ( 2001 ) . c. nayak , phys
b 62 * , 4880 ( 2000 ) . t. timusk and b. statt , rep .
progr . phys .
* 62 * , 61 - 122 ( 1999 ) .
s. tewari _
et al . _ , , 224516 ( 2001 ) .
j.l . tallon and j.w .
loram , physica c * 349 * , 53 ( 2001 ) ; j.w.loram _ et al .
_ , j. phys .
solids @xmath215,2091 ( 1998 ) ; j.l.tallon _ et al .
b @xmath216 , 531 ( 1999 ) . y. j. uemura _
et al . _ ,
lett , @xmath217 , 2317 ( 1989 ) ; c. panagopoulos _ et al .
b * 60 * , 14617 ( 1999 ) ; c. panagopoulos _ et al .
_ , physica c , * 341 - 348 * , 843 ( 2000 ) .
s. chakravarty , c. nayak , s. tewari , , 100504/1 - 4 ( 2003 ) .
homes _ et al .
_ , , 1645 ( 1993 ) . c. stock _ et al
_ , , 024505 ( 2002 )
. h. a. mook _
et al . _ , , 144513 ( 2002 ) .
s. chakravarty , h. y. kee , and c. nayak , _ intl . j.
phys . _ * 15 * , 2901 ( 2001 ) .
s. chakravarty , c. nayak , s. tewari , x. yang , , 277003/1 - 4 ( 2002 ) .
et al . _ ,
nature * 424 * , 912 ( 2003 ) .
s. tewari _
et al . _ , , 014514 ( 2004 ) .
_ , , 137002 ( 2004 ) . c. wu , j. zaanen , s - c
zhang , , 247007 ( 2005 ) .
u. schollwoeck _
et al . _ , 186401 ; j. o. fjaerestad , j. b. marston , u. schollwoeck ann .
( n.y . ) * 321 * , 894 ( 2006 ) . ;
s. capponi , c. wu , s - c .
zhang , 220505(r ) ( 2004 ) .
e. cappelluti and r. zeyher , phys . rev . *
b 59 * , 6475 - 6486 ( 1999 ) .
kee and y. b. kim , phys . rev . * b 66 * , 012505 ( 2002 ) .
c. nayak , e. pivovarov , , 064508 ( 2002 ) .
h. j. schulz , j. de physique * 50 * , 2833 ( 1989 ) .
s. chakravarty , phys .
rev.*b 66 * , 224505 ( 2002 ) .
a. kotani , j. phys .
jpn , * 42 * , 408 ( 1977 ) .
l. pryadko , _ et . al .
_ , phys . rev . *
b 60 * , 7541 ( 1999 ) .
a. damascelli , z. hussain , z - x shen , rev .
phys . * 75 * , 473 ( 2003 ) .
s. kivelson , _ et .
al . _ , rev .
phys . * 75 * , 1201 ( 2003 ) .
r. jamei , s. kivelson , b. spivak , cond - mat/0408066 . c. bena _ et .
_ , phys . rev . *
b 69 * , 134517 ( 2004 ) ; c. bena _ et .
_ cond - mat/0311299 ; a. ghosal , a. kopp , s. chakravarty , cond - mat/0412241 .
v. tripathi , p. chandra , p. coleman , j. phys . : condens .
matter 17 , 5285 ( 2005 ) .