when a quantum field theory is constructed on a curved background spacetime , the gravitationally induced vacuum polarization typically induces a non - zero vacuum expectation value for the stress - energy tensor  @xcite .
it should be pointed out that fully self consistent calculations are horrendously difficult ; consequently most calculations in the literature are carried out in the test - field limit .
( wherein the vacuum polarization is _ not _ permitted to back - react on the geometry via the einstein field equations . ) even in the test - field limit , most known results are obtained by numerical ( rather than analytic ) computations .    _ perturbative _ self - consistent solutions around flat spacetime have recently been investigated by flanagan and wald  @xcite , building on earlier work of horowitz  @xcite , and horowitz and wald  @xcite .
another type of energy condition , based on the `` quantum inequalities '' , is investigated for a schwarzschild background in recent papers by ford and roman  @xcite .
ford and roman have also discussed the anec and awec in ( 1 + 1)dimensional and ( 3 + 1)dimensional evaporating black holes ( unruh vacuum )  @xcite , and for the ( 1 + 1)dimensional boulware vacuum .
in contrast , the present paper works strictly in ( 3 + 1 ) dimensions and treats the equilibrium hartle  hawking vacuum state .
now we know , on rather general grounds ( the existence of the hawking radiation effect and concomitant violation of the _ classical _ area increase theorem of black hole dynamics ) that at least some of the energy conditions must be violated at or near the event horizon of a black hole . in this paper
i shall explore these issues in a little more detail .
i shall continue to work in the test - field limit : i discuss the effect of the gravitational vacuum polarization on various energy conditions .
in particular , i discuss the point - wise null , weak , strong , and dominant energy conditions ( nec , wec , sec , and dec ) , the averaged null energy condition ( anec ) , and shall furthermore introduce and discuss several new , potentially interesting , energy conditions ( pnec , scri
nec , and scri  anec ) .    for the purposes of this paper ,
i restrict attention to the most well - studied curved - space quantum field theory : the conformally - coupled massless scalar field on schwarzschild spacetime , in the hartle  hawking vacuum
. for this geometry and vacuum state one has both ( 1 ) a useful analytic approximation to the gravitational polarization
page s approximation  @xcite , and ( 2 ) numerical estimates of the vacuum polarization  see the numerical calculations of howard  @xcite , howard and candelas  @xcite , and anderson , hiscock , and samuel  @xcite . )
having established the basic machinery in this paper , a study of the boulware vacuum will be presented in a companion paper , while further studies of the unruh vacua should be straightforward , although tedious .
( for the current state of affairs in the unruh vacuum , see  @xcite . ) similarly , extensions to non - conformal couplings , massive scalars , and reissner
nordstrm geometries should be straightforward .
to set the stage , recall some basic definitions  + * definition [ nec ] * * definition [ anec ] * here @xmath0 is a generalized affine parameter ( see  @xcite ) for the null curve , the tangent vector being denoted by @xmath1 .
note that in almost all applications it suffices to consider null geodesics , rather than generic null curves .
however , isolated discussions of the anec on non  geodesic null curves do occur in the literature  @xcite  the extended definition presented here allows one to compare and contrast the present results with other calculations .    as a practical matter
i shall often replace the word `` inextendible '' by the phrase `` inextendible past the event horizon '' .
this is a purely pragmatic decision based on a lack of numerical data inside the event horizon , coupled with the general feeling in the community that page s analytic approximation will become progressively worse as one approaches the central singularity .
( thus my studies of the anec can more precisely be referred to as studies of a `` truncated  anec '' . )    in an earlier publication  @xcite , ( see also  @xcite ) , i derived a no - go theorem for the anec : i showed that the anec is violated ( in the test - field limit ) whenever the background spacetime has a non - zero _ scale anomaly_. specifically , i proved the + * anec no - go theorem * this no - go theorem , since it depends on a special case of the quantum field theoretic ( conformal ) anomaly , is known to be independent of the choice of vacuum state .
( one can easily generalize this argument to include non - conformally coupled fields , see  @xcite . )    unfortunately , a brief calculation suffices to show that the tensor @xmath2 vanishes on schwarzschild spacetime .
indeed it vanishes on any ricci - flat spacetime ( or indeed on any spacetime conformal to an einstein spacetime ) .
thus schwarzschild spacetime avoids the no  go theorem enunciated above , and it becomes a matter of explicit checking to see whether or not the anec is satisfied or violated on this geometry . note that in doing this explicit checking we will have to be cognizant of our particular choice of vacuum state , and indeed our results will also depend on the particular class of null geodesics under consideration .
in performing the explicit checks alluded to above , i have found it useful to define several new energy conditions that are in some sense intermediate between the nec and anec :    * definition [ pnec ] * now , if one does not place any further constraints on the set @xmath3 the resulting definition , while general , is not particularly useful .
( it is best used as a diagnostic tool to probe interesting regions of the tangent bundle . ) a suitable restriction is :    * definition [ scri  nec ] * ( this definition basically says that one should not worry too much if violations of the nec occur only when one looks along null geodesics that never make it out to null infinity . )    we shall soon see that in the hartle  hawking vacuum the gravitational vacuum polarization of a conformally - coupled massless scalar field residing on the schwarzschild background geometry everywhere satisfies this scri  nec energy condition , though it does not everywhere satisfy the nec . we shall also use the pnec to investigate the regions in phase space ( the tangent bundle ) where suitable positivity conditions hold .    finally , for completeness , i enunciate + * definition [ scri  anec ] * ( this definition basically says that one should not worry too much if violations of the anec occur only along null geodesics that never make it out to null infinity .
after all , null infinity is where you want to work to establish results such as the penrose
woolgar positive mass theorem  @xcite or the friedman  schleich
witt topological censorship theorem  @xcite . )    from the previous comment regarding scri
nec it automatically follows that , in the hartle  hawking vacuum , the gravitational vacuum polarization of a conformally coupled scalar field residing on the schwarzschild background geometry satisfies this scri  anec .
we shall also see that not all null geodesics satisfy the anec itself .
in the particular case of conformally coupled scalar fields residing on schwarzschild spacetime one has a lot of information about the vacuum polarization . by spherical symmetry
one knows that @xmath4.\ ] ] where @xmath5 , @xmath6 and @xmath7 are functions of @xmath8 , @xmath9 and @xmath10 .
note that i set @xmath11 , and choose to work in a local - lorentz basis attached to the fiducial static observers ( fidos ) .
page s analytic approximation gives  @xcite a polynomial approximation to the stress  energy tensor : @xmath12 , \\ % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % %              2        3        4        5        6 %     6 r    9 r     12 r     15 r     18 r     99 r % 3 + --- + ---- + ----- + ----- + ----- - ----- %      r       2       3        4        5        6 %            r       r        r        r        r % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % %    \tau(r ) & = & - p_\infty   \left [   1 + 2\left({2m\over r}\right ) + 3\left({2m\over r}\right)^2   + 4\left({2m\over r}\right)^3 + 5\left({2m\over r}\right)^4   + 6\left({2m\over r}\right)^5 + 15\left({2m\over r}\right)^6 \right ] , \\ % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % %              2       3       4       5        6 %     2 r    3 r     4 r     5 r     6 r     15 r % 1 + --- + ---- + ---- + ---- + ---- + ----- %      r       2       3       4       5       6 %            r       r       r       r       r % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % p(r )     & = & + p_\infty   \left [   1 + 2\left({2m\over r}\right ) + 3\left({2m\over r}\right)^2   + 4\left({2m\over r}\right)^3 + 5\left({2m\over r}\right)^4   + 6\left({2m\over r}\right)^5 - 9\left({2m\over r}\right)^6 \right ] , \\
% % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % %              2       3       4       5       6 %     2 r    3 r     4 r     5 r     6 r     9 r % 1 + --- + ---- + ---- + ---- + ---- - ---- %      r       2       3       4       5       6 %            r       r       r       r       r % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % %         & = & + p_\infty   \left [   1 + 2\left({2m\over r}\right ) + 3\left({2m\over r}\right)^2   \right ] \left [ 1 +   4\left({2m\over r}\right)^3 - 3\left({2m\over r}\right)^4 \right ] .
\end{aligned}\ ] ] here i have defined a constant ( to be interpreted as the pressure at spatial infinity ) by @xmath13 note that i have explicitly expanded the functions given in references  @xcite as polynomials in @xmath14 to exhibit the fact that all components of the stress  energy are well behaved at the horizon .
it is perhaps somewhat surprising that the rather messy rational _ quotients _ of polynomials exhibited in  @xcite reduce to such relatively nice compact expressions .
these formulae have also been checked for consistency with the brown  ottewill extensions of the original page approximation  ( * ? ? ?
* see esp .
see also elster  @xcite .
if one prefers , an alternative form for the page approximation is @xmath15 + q(r )       \left [ \matrix{+3&0&0&0\cr                 0&+1&0&0\cr                 0&0&0&0\cr             0&0&0&0\cr } \right ] \right\}.\ ] ] the polynomials @xmath16 and @xmath17 being @xmath18 and @xmath19    the trace of the stress  energy tensor is given by @xmath20 this result , because it is simply a restatement of the conformal anomaly , is known to be exact .    at the event horizon ,
page s analytic approximation gives  @xcite @xmath21.\ ] ]    the explicit numerical calculations of howard  @xcite , howard and candelas  @xcite , and anderson , hiscock , and samuel  @xcite show that page s analytic approximation is reasonably good  the worst deviations occur in the immediate vicinity of the event horizon where the anderson  hiscock
samuel numerical analysis gives @xmath22.\ ] ] a convenient table of numerical results , covering the range @xmath23 to @xmath24 may be found on page 2541 of reference  @xcite .
( near the event horizon , the howard
candelas data is believed to be accurate to two significant digits .
the accuracy is expected to improve with increasing @xmath8 , as the numerical results converge upon the page approximation . )    additionally , i have obtained via private communication the more intensive numerical data of anderson  hiscock
samuel  @xcite , which covers the range @xmath23 to @xmath25 at much higher resolution .
( near the event horizon , the anderson  hiscock
samuel data is believed to be accurate to three significant digits . )
i have collated the data and used it to construct suitable interpolating functions to be used for the numerical aspects of the following discussion .    to construct the interpolating functions i use the anderson  hiscock
samuel data in the range @xmath23 to @xmath25 , and the howard
candelas data in the range @xmath25 to @xmath24 , augmented by the fact that we know the exact result at @xmath26 .
using mathematica , i have fit the data using third - order interpolating polynomials in the variable @xmath27 .
( note that @xmath28 conveniently runs from @xmath29 at the horizon to @xmath30 at spatial infinity ; @xmath28 is also _ the _ natural variable to use when numerically analyzing the page approximation . )    of course , there are may other papers extant which calculate the stress - energy tensor in somewhat different configurations .
see , for example , .
more recently , in a very interesting development , matyjasek  @xcite has developed a curve - fitting analysis that fits the numerical data to high accuracy . in this paper
i prefer to work directly with the numerical data . in this paper
i am only presenting the minimum requirements for the particular job at hand and it is clear that with additional effort more information can be extracted by mining this additional vein of data .
outside the event horizon , the nec reduces to the pair of constraints @xmath31 page s approximation yields @xmath32.\ ] ] ( note the factorization ! ) consequently @xmath33 is always explicitly positive outside the horizon .
furthermore @xmath34.\ ] ] numerically finding the roots of this sixth order polynomial shows that @xmath35 is negative from @xmath23 to @xmath36 .
thus the page approximation suggests that the null energy condition is violated in the range @xmath37 $ ] .    turning to the numerical data , one easily verifies that that outside the horizon @xmath38 , while @xmath39 at the horizon .
this is in complete agreement with page s analytic approximation .    on the other hand , inspection of the numerical data
indicates that @xmath40 for @xmath41 .
( this number was obtained by fitting third - order interpolating polynomials to the numeric data ; and then numerically finding the root . )
this is qualitatively , though not quantitatively in agreement with page s analytic approximation .
thus the numerical data indicate that the null energy condition is violated in the range @xmath42 $ ] .      inside the event horizon ,
the radial coordinate becomes timelike , and the roles played by @xmath43 and @xmath44 are interchanged .
the nec reduces to the pair of constraints @xmath45    unfortunately , inside the horizon i am reduced to reliance upon the page approximation by a total lack of numerical data .
now there are some subtleties associated with extrapolating the page approximation inside the event horizon .
the existence of a killing vector that points in the @xmath46 direction is absolutely critical , ( even if the @xmath46 direction is no longer timelike ) .
thus if you wish to extrapolate the page approximation inside the event horizon it is necessary to limit attention to an _ eternal _ black hole ( that is , the maximally extended kruskal
szekeres manifold ) .
i should remind the reader that in the hartle  hawking vacuum the stress - energy is at least known to be regular at the horizon , so the page approximation should not be too far wrong just inside the event horizon of an eternal black hole . on the other hand , one does expect the page approximation to get progressively worse as one moves further in toward the singularity , so one should probably not take these results too seriously once one is far inside the horizon .
i should especially warn the reader against extrapolating the page approximation to the interior of an astrophysical black hole : we really do not expect the interior of an astrophysical black hole to have the same killing vectors as the exterior , but instead expect the interior to be fully dynamical .
( in addition we do not expect an astrophysical black hole to even to be in the hartle  hawking quantum state . )
the present calculations inside the horizon are strictly limited to eternal black holes and are exhibited merely because they are do - able , and because they may be suggestive of the actual situation .
inspection of page s analytic approximation indicates that @xmath47 for all @xmath48 , so this part of the nec is satisfied inside the event horizon . on the other hand @xmath49 which is everywhere negative , ( both inside and outside the horizon ) .
thus the page approximation suggests that the null energy condition is violated throughout the interior of the black hole .
this is enough to inform us that , ( assuming the reliability of the page approximation in this matter ) , the weak , strong , and dominant energy conditions are also violated throughout the entire interior of the black hole .
outside the horizon , the weak energy condition is equivalent to the three constraints @xmath50 the last two constraints have already been discussed
they are simply the null energy condition . in the page approximation , finding the roots of the relevant sixth - order polynomial shows that the density constraint is violated for @xmath51 .
therefore the weak energy condition is violated in the range @xmath52 $ ] .
in contrast , the numerical data indicate that violations of the weak energy condition extend to the region @xmath53 $ ] .      inside the horizon ,
the weak energy condition is equivalent to the three constraints @xmath54 the last two constraints have already been discussed ( when dealing with the nec ) .
the violation of the nec is already sufficient to tell us that the wec is violated everywhere inside the horizon , but we can add in passing that the page approximation additionally implies that @xmath44 is everywhere negative inside ( and outside ) the horizon .
the page approximation suggests that the weak energy condition is violated throughout the interior of the black hole .
outside the horizon , the strong energy condition is equivalent to the three constraints @xmath55 we have already looked at the first two constraints when discussing the null energy condition .
the third condition is always satisfied outside the horizon ( since both @xmath56 , and @xmath57 in this region ) .
this is true both in the page approximation , and by appeal to the numerical data .
thus the strong energy condition is violated in the same region as the null energy condition .      inside the horizon ,
the strong energy condition is equivalent to the three constraints @xmath58 we have already looked at the first two constraints when discussing the null energy condition , and thereby know that ( the page approximation suggests that ) the strong energy condition is violated throughout the entire interior of the black hole .    for completeness
i point out that the first condition is satisfied throughout the interior , the second of these conditions is violated throughout the interior , while the third condition also fails throughout the interior .
outside the horizon , the dominant energy condition is equivalent to the three constraints @xmath59 ?
\nonumber\\ & & p(r)\in [ -\rho(r),+\rho(r)]?\end{aligned}\ ] ] we can rephrase this as @xmath60    using the page approximation , and restricting attention to the range @xmath61 $ ] , one has : @xmath62 ;    \nonumber\\ & \rho(r)-\tau(r )
< 0   & \qquad      r\in\emptyset ;       \nonumber\\ & \rho(r)+\tau(r ) < 0 & \qquad      r\in[2m,2.77256 m ] ;   \nonumber\\ & \rho(r)-p(r ) < 0     & \qquad      r\in[2m,2.58512 m ] ;   \nonumber\\ & \rho(r)+p(r ) < 0     & \qquad      r\in[2m,2.18994m].\end{aligned}\ ] ] pulling this all together , the page approximation suggests that dominant energy condition fails in the region @xmath63 $ ] .
the numerical data implies quantitative though not qualitative modifications . in the range
@xmath61 $ ] one finds : @xmath64 ;    \nonumber\\ & \rho(r)-\tau(r ) < 0   & \qquad      r\in\emptyset ;      \nonumber\\ & \rho(r)+\tau(r ) < 0 & \qquad      r\in[2m,2.992 m ] ;    \nonumber\\ & \rho(r)-p(r ) < 0     & \qquad      r\in[2m,2.628 m ] ;    \nonumber\\ & \rho(r)+p(r ) < 0     & \qquad      r\in[2m,2.298m].\end{aligned}\ ] ] pulling this all together , the numerical data indicate that the dominant energy condition fails in the region @xmath65 $ ] .
note that this is suspiciously close to @xmath66the unstable circular photon orbit  and _ might _ be trying to tell us something .    if one actually calculates @xmath67 at @xmath66 one gets @xmath68given the expected three - significant - digit numerical reliability of the data this is ( questionably ) compatible with zero .
if one takes this suggestion seriously , it would imply a hidden ( accidental ? ) symmetry in the stress tensor at @xmath66 : @xmath69?\ ] ] it is quite possible however , that this is purely a numerical accident .    finally , i point out that this conjectured property certainly does not survive the introduction of non - conformal coupling , nor is there any particular reason to expect it to survive the introduction of non - zero rest mass .
this conjecture also most definitely does not hold in the boulware or unruh vacuum states  @xcite .      inside the horizon ,
the dominant energy condition is equivalent to the three constraints @xmath70 ?   \nonumber\\   & & p(r)\in [ -\tau(r),+\tau(r)]?\end{aligned}\ ] ] we can rephrase this as @xmath71 we already know that @xmath72 inside the horizon ( in fact everywhere ) .
so , provided the page approximation is not misleading in this regard , the dominant energy condition is violated throughout the interior of the black hole .
no additional information comes from the other conditions .
to analyse the partial null energy condition introduced earlier in this paper , consider a generic null vector inclined at an angle @xmath73 away from the radial direction . then without loss of generality , in an orthonormal frame attached to the @xmath74 coordinate system , @xmath75 ignoring the ( presently irrelevant ) overall normalization of the null vector , one has @xmath76
+ [ \tau + p ] \sin^2\psi).\end{aligned}\ ] ] we have already seen that @xmath77 is positive outside the event horizon . on the other hand
the page approximation gives @xmath78 which is everywhere negative .
a glance at howard s numerical data confirms that the numerical data also satisfies @xmath79 .
thus the partial null energy condition _ fails _ for those radii @xmath8 , and those angles @xmath73 , such that @xmath80 \equiv    \sqrt { \rho(r)-\tau(r ) \over |\tau(r ) + p(r)| } .\ ] ] that is , the partial null energy condition _ fails _ for @xmath81 $ ] . defining @xmath27
this critical angle is ( in the page approximation ) given by @xmath82.\ ] ] for @xmath83 ( @xmath84 ) there are no real solutions to this equation . at @xmath85 one
has @xmath86 , while @xmath87 moves monotonically to zero as @xmath88 .
( see figure  [ fig1 ] where @xmath89 is plotted as a function of @xmath28 . )
thus , sufficiently near the horizon , almost all directions violate the partial null energy condition .
as one moves further away from the horizon the violations of the partial energy condition are confined to null vectors that are progressively more and more transverse , finally at @xmath85 the violations of the partial null energy condition disappear completely .
using the numerical data will modify the precise location where this behaviour manifests itself , but will not qualitatively modify this picture .
the violations of the partial null energy condition vanish outside @xmath90 ( @xmath91 ) , and the numerically determined values of @xmath89 are superimposed on figure  [ fig1 ] .      inside the event horizon , one should consider a generic null vector inclined at an angle @xmath92 away from the @xmath46 direction ( which is now spacelike ) . then without loss of generality , in an orthonormal frame attached to the @xmath74 coordinate system , @xmath93 one should now consider the quantity @xmath94   + [ \rho + p ] \sin^2\tilde\psi).\end{aligned}\ ] ] note that inside the horizon @xmath95 is strictly positive , while @xmath96 is strictly negative .
in fact , the page approximation yields @xmath97.\ ] ] one deduces that violations of the partial null energy condition occur whenever @xmath98 \equiv   \sqrt { \tau(r)-\rho(r ) \over |\rho(r)+p(r)| } \ ] ] consequently , violations of the partial null energy condition are confined to the range @xmath99 $ ] . again defining @xmath27 this critical angle is ( in the page approximation ) given by @xmath100.\ ] ] one has , at @xmath101 , @xmath102 , while @xmath103 falls to zero as @xmath8 approaches @xmath104 . qualitatively : near ( but inside ) the horizon almost all null directions suffer violations of the partial energy condition , while near the central singularity somewhat fewer directions violate the partial null energy condition .
again , i remind the reader to not take the page approximation too seriously as one approaches the singularity .
i now turn attention to the asymptotic null energy condition ( scri  nec ) that i introduced earlier in this paper . discussing
this energy condition will require some standard results concerning the null geodesics of the schwarzschild spacetime .
( see , for instance , such standard textbooks as misner  thorne
wheeler  @xcite , wald  @xcite , or weinberg  @xcite . )
there are two conserved quantities for null geodesic motion in schwarzschild spacetime ; the energy and angular momentum . using these conservation laws
the affine parameter can be chosen in such a way that @xmath105 the parameter @xmath106 is the angular momentum per unit energy . if the null geodesic reaches asymptotic spatial infinity then this parameter can also be interpreted as the `` impact parameter '' .
however , there is a large class of null geodesics that never reaches spatial infinity , for these null geodesics the notion of `` impact parameter '' is at best an abuse of language .
the angle between the null geodesic and the radial direction is given by @xmath107 we are ultimately interested in the quantity
@xmath108   + [ \tau + p ] \sin^2\psi ) \nonumber\\ & = &          \{[\rho -\tau ]   + [ \tau + p ] ( 1 - 2m / r)a^2/r^2\},\end{aligned}\ ] ] but will need to go through a few preliminaries .
the radial motion of null geodesics is governed by @xmath109.\ ] ] so the turning points , @xmath110 , are given by the cubic @xmath111 if the impact parameter is small , @xmath112 , then it is a standard result that there are no turning points : the null geodesic either plunges into the future singularity , or emerges from the past singularity of the maximally extended schwarzschild spacetime .
( for @xmath113 the geodesic may make a large number of `` orbits '' before crossing the event horizon . )
if the impact parameter is marginal , @xmath114 , then it is a standard result that @xmath66 , corresponding to the ( unstable ) circular photon orbit .    if the impact parameter is large , @xmath115 , then it is a standard result that there are _ two _ turning points at physical values of @xmath8 .
one of these turning points lies in the range @xmath116 , while the other lies in the range @xmath117 .
note that if @xmath118 then the three mathematical roots of the cubic are approximately @xmath119 , and @xmath120 $ ] . the two physical roots are approximately @xmath121 and @xmath122 $ ] .
the first of these turning points , @xmath116 , corresponds to the obvious class of null geodesics with high impact parameter
those geodesics that come in from spatial infinity , bounce off the angular momentum barrier , and return to spatial infinity .
( for @xmath123 the geodesic may make a large number of `` orbits '' before returning to spatial infinity . )
the second of these turning points , @xmath117 , corresponds to a completely separate class of null geodesics with high `` impact parameter''these geodesics emerge from the event horizon at @xmath124 , with high angular momentum , make a large number of `` orbits '' before reaching their maximum height above the event horizon , and then make an equally large number of `` orbits '' before returning to re - cross the event horizon at @xmath125 . for these geodesics
the use of the phrase `` impact parameter '' to describe the parameter @xmath106 is most definitely an abuse of language .    to analyze the asymptotic null energy condition ,
start by first considering the quantity @xmath126 which i shall refer to as the nec density .
( remember that @xmath127 is an explicitly known function of @xmath8 and @xmath106 . )
those null geodesics that come in from infinity and return to infinity never get closer to the origin than @xmath66 .
see , for instance ,  @xcite .
inspection of either the page approximation or of the numerical data indicates that @xmath5 , @xmath128 , and @xmath7 are all positive for @xmath129 .
thus the asymptotic null energy condition is satisfied along all null geodesics that come from , and return to , infinity .    for other classes of null geodesics
it proves convenient to rewrite the quantity @xmath130 as : @xmath131   + [ \tau + p ] ( 1 - 2m / r)a^2/r^2\}.\ ] ] now , for incoming null geodesics with smaller than critical impact parameter , the null geodesic may circle the black hole a large number of times , but is guaranteed to ultimately plunge into the event horizon  @xcite .
this makes the analysis a little more subtle .
inspection of either the page approximation or the numerical data shows that ( outside the horizon ) @xmath77 is always positive , while @xmath132 is always negative .
now use the fact that for an infalling null geodesic @xmath112 . since @xmath132 is negative , this implies ( for this class of null geodesics ) a lower bound @xmath133
+ [ \tau + p ] ( 1 - 2m / r)27m^2/r^2\}.\ ] ] inspection of page s approximation indicates that this lower bound is strictly positive ( zero at the horizon ) .
in fact @xmath134.\ ] ] ( note the factorization ! ) thus the page approximation suggests that the asymptotic null energy condition holds on all infalling null geodesics .
( at least until one crosses the horizon ! )
this lower bound @xmath135 is plotted in figure  [ fig3 ] .    by time reversal
, this suggests that the asymptotic null energy condition also holds on all outgoing null geodesics that reach infinity .    with respect to the numerical data
, one can easily evaluate @xmath135 numerically .
the results are superimposed on figure  [ fig3 ] .
while there are marked differences between the analytic approximation and the numerical data , both curves are seen to be strictly positive outside the horizon .
we conclude that scri
nec is satisfied . for _ all _ null geodesics that reach spatial infinity
the nec density @xmath130 is strictly positive everywhere outside the event horizon .    for completeness , and future reference ,
i point out that the page approximation yields @xmath136 .
\label{eq_pi_r_a}\ ] ] again , note the factorization !
now turn attention to the trapped null geodesics .
( these are null geodesics with high `` impact parameter '' @xmath137 , trapped in the region @xmath138 $ ] . ) while these trapped null geodesics are , by definition , not directly relevant to the scri
nec , the tools developed above permit us to gain additional insight into the pnec on trapped null geodesics .
pick some value of @xmath8 in the range @xmath139 .
then @xmath130 is guaranteed to be negative if one chooses @xmath140 ( setting @xmath27 is again a convenient choice of variables for both analytic and numeric work . in this section
we will only be interested in the region @xmath141 $ ] . )
note that what we are doing is guaranteeing that with these definitions the pnec is violated for the region of the @xmath142 plane above the curve @xmath143 .
using page s approximation , this critical impact parameter is given by @xmath144 the critical impact parameter implied by the numerical data was also determined and both curves are plotted on figure  [ fig4 ] .
it should be noted that , for a given value of @xmath117 , one can not choose @xmath106 arbitrarily : in order for a trapped null geodesic with impact parameter @xmath106 to ever reach radius @xmath8 it is necessary that @xmath106 be small enough .
indeed one must have @xmath145 this may be thought of as a kinematic bound on trapped null geodesics : the region of the @xmath142 plane above the curve @xmath146 is kinematically inaccessible to trapped null geodesics .
( note that the region below @xmath114 is also kinematically inaccessible . )
the relevant curve is plotted in figures  [ fig5 ] and  [ fig6 ] , where it is overlain with @xmath147 as obtained from figure  [ fig4 ] . by looking at where these curves cross one another
we can draw some general conclusions .    in the page approximation , now restricting attention to the region outside the event horizon
:    1 .   some trapped null geodesics ( those with @xmath148 = [ 5.19615 m , 6.1101 m]$ ] )
_ never _ experience nec density violations .
that is : @xmath149 is _ satisfied _ along the entire portion of the null geodesic that lies outside the event horizon .
( and therefore the [ truncated ] anec must be satisfied on this entire class of null geodesics . ) 2 .   some trapped null geodesics ( those with @xmath150 ) _ always _ experience nec density violations .
that is : @xmath151 along the entire [ truncated ] null geodesic .
( and therefore the [ truncated ] anec must be violated on this entire class of null geodesics . )
all other trapped null geodesics ( those with @xmath152 $ ] ) will experience some nec violations when they get sufficiently close to the event horizon .
( and therefore investigating anec violations on this class of null geodesics requires more work . ) 4 .
note that @xmath153 corresponds to a null geodesic that reaches a maximum radius @xmath154 , a number that we have seen before ( when discussing the nec and pnec ) .
use of the numerical data implies quantitative though not qualitative changes .
we observe that    1 .
the set of trapped null geodesics which _ never _ experience nec density violations is much smaller  those with @xmath155 = [ 5.19615 m , 5.276 m]$ ] .
( the [ truncated ] anec must be satisfied on this entire class of null geodesics . )
2 .   some trapped null geodesics ( those with @xmath156 ) _ always _ experience nec density violations .
( therefore the [ truncated ] anec must be violated on this entire class of null geodesics . )
all other trapped null geodesics ( those with @xmath157 $ ] ) will experience some nec violations when they get sufficiently close to the event horizon .
( and therefore investigating anec violations on this class of null geodesics requires more work . ) 4 .   note that @xmath158 corresponds to null geodesic that reaches a maximum radius @xmath159 , a number that we have seen before ( when discussing the pnec ) .
the message to be extracted is this : the numeric data and analytic approximations are in good qualitative agreement with each other .
null geodesics that reach null infinity are well behaved ( scri  nec is satisfied ) , but _ some _ of the trapped null geodesics encounter nec violations ( and anec violations ) .
we are finally in a position to pin down precisely violations of the anec itself .
although significant information can already be extracted by using the point - wise energy conditions already discussed , beyond a certain stage explicit use of the affine parameterization must be invoked .
write the anec integral as  @xcite @xmath160   + [ \tau + p ] \sin^2\psi ) \ ; \xi^2 \ ; d\lambda , \\ & = & \int_\gamma \pi(r , a ) \ ; \xi^2 \ ; d\lambda.\end{aligned}\ ] ]    note that the _ integrand _ is proportional to the quantity @xmath130 which has already been extensively discussed in the context of the scri  nec .
the arguments previously adduced for the scri  nec can be carried over wholesale to the scri  anec . in particular :    * for null geodesics that come in from spatial infinity and return to spatial infinity @xmath161 the _ integrand _ is everywhere positive and so the anec holds on this entire class of geodesics .
( note that satisfaction of the anec does not arise from the trivial observation that the heat bath contributes an asymptotically constant and positive energy density far from the black hole  rather one has the stronger statement that the integrand itself is positive along the entire geodesic . ) * the fact that the integrand is positive implies that anec will also be satisfied for curves with transverse smearing  @xcite , or curves with arbitrary positive weighting functions  @xcite . * for infalling or outfalling null geodesics that reach asymptotic spatial infinity @xmath162
the integrand is everywhere positive and anec is satisfied .
[ caveat : i stop the anec integral once it reaches the horizon . ]
* this may be interpreted as follows : because scri
nec is satisfied , it automatically follows that scri
anec is satisfied .
* for the unstable circular photon orbit @xmath163 the integrand is everywhere positive and anec is satisfied .
* however , for trapped null geodesics @xmath164)$ ] the integrand is no longer necessarily positive , and this is the only case for which we will explicitly need to look at the `` weighting function '' @xmath165 appearing in the anec integral .      observe that the `` weighting function '' appearing in the anec satisfies  ( * ? ? ?
* , equations ( 12.60)(12.63 ) ) @xmath166 although this calculation was carried out for the schwarzschild geometry the result remains true for an arbitrary static spacetime
the anec integral is simply the time - average along a null geodesic of the local - lorentz nec integrand where the time average is to be taken with respect to the natural static time coordinate .    for actual calculations it is much more practical to re - express this as an integral with respect to the radial variable @xmath8 by using @xmath167 so that @xmath168 now this observation appears to weight the region near the event horizon very heavily  because of the explicit pole at @xmath23 .
however , the integrand @xmath130 has a zero at the event horizon  in the page approximation one discovers an explicit factor of @xmath169 [ _ cf .
_ equation [ eq_pi_r_a ] ] , while the numerical data also exhibits a first - order zero in @xmath130 at the horizon .
thus for trapped null geodesics one may write @xmath170 and have some confidence that the integral actually converges at the lower bound @xmath23 .
as a penultimate step , recall that calculating @xmath171 involves solving a cubic .
it is more convenient to parameterize the trapped geodesic by calculating the impact parameter in terms of the maximum height attained by the null geodesic : @xmath172 so for these trapped null geodesics one finds the anec integral is @xmath173/[(1 - 2m / r_{\rm max})r^2]}},\end{aligned}\ ] ] to actually evaluate this integral numerically it is useful to change variables to @xmath27 , and @xmath174 , with the result that @xmath175/[z_0 ^ 2(1-z_0)]}},\end{aligned}\ ] ] with a little work , the square root can be seen to factorize explicitly @xmath176 this integral , though singular at the lower limit ( @xmath177 , corresponding to @xmath178 ) , is now certainly finite .
while mathematica has resources to deal with @xmath179 singularities at the endpoints of the integration range , it reacts badly to such singularities when the location is chosen dynamically . that is :
@xmath180 is handled badly .
for this reason the change of variables @xmath181 , while being a formal mathematical identity , leads to much better numerical behaviour for the integral : @xmath182    this integral was determined numerically both for the page approximation and for the anderson  hiscock
samuel numeric data .
the results are graphed in figure  [ fig7 ] .
remember that it is meaningless to force this particular integral out of the range @xmath138 $ ] , corresponding to @xmath183 $ ]    thus we have ( finally ! ) managed to characterize the precise class of null geodesics on which the [ truncated ] anec is violated .
the critical values of @xmath184 , @xmath185 , and @xmath106 are , for the page approximation : @xmath186 for the numerical data , one obtains : @xmath187 note that these numbers are qualitatively reasonable and in agreement with the violations of the point - wise energy conditions , and the numerical investigation of the pnec .
furthermore the anec violations in the numerical data are seen to extend out to larger distances than the anec violations in the page approximation , in agreement with the general trend .    though i have quoted four - significant - digit accuracy for the numerical data one should probably not take anything past the second significant digit too seriously . in fact , given the vagaries of numerical analysis on singular integrals ( and the inherent uncertainties in the anderson  hiscock
samuel data ) it is conceivable that the exact critical impact parameters might be @xmath188 and @xmath189 respectively .      if one considers arbitrary non - geodesic curves the anec loses much ( though not quite all ) of its power : by a judicious choice of non - geodesic null curve one could try to remain in a region where the nec is violated and thereby `` trivially '' violate the anec  @xcite .
for instance , along any circular null curve at fixed @xmath8 ( _ not _ a geodesic except in the case of @xmath66 ) the anec integral is still proportional to @xmath96 .
inspection of the numeric data indicates that @xmath40 for @xmath41 .
( page s analytic approximation giving a slightly different result , @xmath190 for @xmath191 . )
thus the anec is _ not _
satisfied for this particular class of non  geodesic null curves .
investigation of the properties of the averaged null energy condition is of considerable interest to diverse applications in semiclassical quantum gravity .
it is now abundantly clear that , in the test - field limit , semiclassical quantum fields do _ not _ generally satisfy the anec : indeed anec violations are related to the existence of a non - zero scale anomaly  @xcite . even if the scale anomaly vanishes , this does not necessarily imply that the anec is satisfied : one has to do a case by case analysis . as an example
, this paper investigates the situation in schwarzschild spacetime .
the analysis presented herein is somewhat of an attempt to crack a walnut with a sledgehammer in the sense that it is a collection of rather general techniques applied to a rather particular problem
it is certainly true that this type of analysis can now be carried forward to other geometries , other quantum fields , and other vacuum states by by straightforward but tedious computation .
* point - wise energy conditions * * inside the event horizon , with suitable caveats regarding the applicability of the page approximation , almost any energy condition you can think of will be violated .
* * between the event horizon @xmath192 and the unstable photon orbit @xmath193 many of the energy conditions are violated , in a series of onion - like layers . * * outside the unstable photon orbit @xmath193 all energy conditions are satisfied . *
* all null geodesics that reach asymptotic infinity are well - behaved : if you look _ along _ the null geodesic you never see nec violations . *
averaged null energy condition * * if your null geodesic ever reaches infinity , the [ truncated ] anec is definitely satisfied , and is satisfied for a non - trivial reason : the _ integrand _ is strictly positive all the way from the event horizon to null infinity . * * between the event horizon @xmath192 and
the unstable photon orbit @xmath193 some of the trapped null geodesics violate the [ truncated ] anec . * * if you are willing to look at non - geodesic null curves even more violations of the anec can be found .      1
.   it would be very nice to have an analytic understanding of the precise role played by the unstable photon orbit  the numerical evidence is suggestive but not definitive .
2 .   for that matter it would be nice to know if the relevance of the unstable photon orbit generalizes to other geometries .
3 .   it would be nice to go beyond the numerics ; to develop some exact analytic arguments that go beyond the page approximation .
generalizations to the boulware vacuum will be presented in a companion paper .
generalizations to the unruh vacuum , other quantum fields , non - conformal couplings , particle masses , and the reissner  nordstrm geometry will be straightforward if tedious .
the new energy conditions i introduce , scri  nec and scri  anec , are interesting in that they focus attention on null infinity .
and null infinity is where all the interesting details arise in the friedman
witt  @xcite topological censorship theorem , and the penrose  sorkin
woolgar version of the positive mass theorem  @xcite .
i suspect that with a little more work suitable generalizations of these theorems can be constructed in terms of the scri  anec .
finally i should point out that even though the various energy conditions are violated in many regions , this does not give one a completely free hand to design spacetime geometries to taste : it seems quite likely that the `` quantum inequality '' approach of ford and roman  @xcite will allow us to place constraints on spacetime geometries even if all the usual types of energy condition fail .


one of the most active areas of condensed matter theory is the development of new algorithms to simulate and predict the behavior of materials exhibiting strong correlations .
recent developments of the powerful many - body approach , the dynamical mean field theory ( dmft ) @xcite and its cluster extensions @xcite , hold a great promise in being able to accurately predict physical properties of this challenging class of materials . in recent years , the dmft method has substantially advanced our understanding of the physics of the mott transition and demonstrated its power to explain such problems as the structural phase diagrams of actinides @xcite , phonon response @xcite , optical conductivity @xcite , valence and x - ray absorption @xcite and transport @xcite of some of the archetype materials with strong correlations .
the success of the dynamical mean field theory in the context of the electronic structure revitalized the search for fast and flexible impurity solvers which could treat hund s coupling and spin - orbit coupling of the parent atomic constituents in the crystal environment of the lattice .
further , the newly developed cluster extensions of dmft @xcite require faster impurity solver which could access low temperature strong correlation limit .
many impurity solvers were developed over the last few decades , including hirsh - fye quanum monte carlo method @xcite , exact diagonalization @xcite , non - crossing approximation @xcite , and its extensions like one - crossing approximation @xcite or sunca @xcite , iterative perturbation theory @xcite , wilson s numerical renormalization group @xcite expansion around the atomic limit @xcite , and many others .
each method has some advantages and disadvantages , but at the present time , there is no method that works efficiently and produces accurate solutions for the green s function in all regimes of parameters . concentrating only on the most often employed method , numerically exact hirsch - fye quanum monte carlo ,
the following weaknesses limit its usefulness in many realistic materials and clusters of strongly correlated models : i ) it can not treat realistic multiplet structure which is very important in strongly correlated @xmath0 materials .
ii ) the discretization of the imaginary time leads to considerable systematic error @xcite and requires extrapolation to infinitesimally small time slices .
iii ) the low temperatures regime in the strong correlation limit of large @xmath1 is computationally very expensive since it requires many time slices and infinite @xmath1 models like t - j model are unaccessible .
all of the above mentioned shortcomings of the conventional quantum monte carlo algorithm are eliminated by the novel continuous time quantum monte carlo method @xcite .
in addition , the new method is even much faster for most of applications we tested , including cluster dmft for the hubbard model and application of lda+dmft to the actinides .
the basic idea and its implementation for the hubbard model was recently published in ref .  .
further extensions to more general impurity model was implemented in ref .  .
first demonstration of the power of our implementation was presented in ref .   by detail
low temperature study of sodium doped cobaltates and in ref .   to study plutonium valence .    here
we want to describe the powerful new implementation of the method for applications to realistic materials with complicated multiplet structure , including all interaction terms : hubbard @xmath1 , hund s coupling and spin - orbit coupling .
the extension to clusters of few sites and superconducting state within cluster dmft will be addressed and its power will be demonstrated by studying the low temperature coherence of both the hubbard and the t - j model . the kondo volume collapse transition in elemental cerium will be reexamined showing the valence histogram of the alpha and gamma phase of the material .
more technical details of the implementation are given in the appendix [ appa ] .
in this section , we will explain the main steps of the recently developed quantum monte carlo method with emphasis on the generalization to clusters and multiplets of real materials .
the cluster dynamical mean field approach can be conveniently expressed by the functional of the local green s function @xcite @xmath2 = { \mathrm{tr}}\log(g_0^{-1}-\sigma)-{\mathrm{tr}}[\sigma g ] + \phi[g_{loc}].\end{aligned}\ ] ] here @xmath3 stands for the green s function of the cluster or single site under consideration and is in general a matrix of the size equal to the number of sites times the number of correlated orbitals per site .
further , @xmath4 is the non - interacting green s function @xmath5 and contains all quadratic terms of the hamiltonian including periodic ionic potential of the crystal ( @xmath6 ) . the interacting part of the functional @xmath7 $ ] contains _ all _ two particle irreducible skeleton diagrams inside the chosen cluster , i.e. , @xmath8/\delta g$ ] .    within the dmft method , the summation of all the diagrams
is achieved by solving the corresponding quantum impurity problem @xmath9 e^{-s_{c } -     \int_0^\beta d\tau\int_0^\beta d\tau ' \sum_{{\alpha}{\alpha } ' }     \psi_{\alpha}^\dagger(\tau)\delta_{{\alpha}{\alpha}'}(\tau,\tau')\psi_{{\alpha}'}(\tau ' ) }   \label{zimp}\end{aligned}\ ] ] where the anderson impurity model hybridization @xmath10 term plays the role of the generalized weiss field that needs to be added to the cluster effective action @xmath11 to obtain the local green s function . the self - consistency condition which determines this weiss field ( hybridization @xmath10 ) is @xmath12    in the general impurity problem defined in eq .
( [ zimp ] ) the electrons in the cluster @xmath11 hybridize with a matrix of weiss fields @xmath13 . in the case of single - site dmft for an @xmath0 material ,
this hybridization is a @xmath14 matrix while in cluster dmft for plaquette , it is @xmath15 matrix . in some cases ,
hybridization can be block diagonalized as we will show in the appendix  [ appa ] below on the example of cellular - dmft for the normal state of the hubbard model .
however , in the superconducting state of the same model , the hybridization acquires off - diagonal components due to anomalous components of self - energy
.    the cluster part of the action @xmath11 can be very complicated and the power of this method is that it can treat arbitrary interaction within the cluster . for real materials study ,
the most important on - site terms are the hund s couplings which are usually expressed by the slater integrals @xcite ( @xmath16 , @xmath17 , @xmath18 in case of @xmath0 electrons ) .
in addition , there is spin - orbit coupling @xmath19 and crystal field splittings as well as various hoppings and non - local interactions within the cluster .
the continuous time impurity solvers are based on the diagrammatic expansion of the partition function and stochastic sampling of the diagrams .
two expansions were recently implemented : the expansion around the band limit @xcite and the expansion in the hybridization strength @xcite .
the later seems to be superior in the strongly correlated regime due to substantial reduction of the size of matrices that need to be manipulated @xcite and , most importantly , the empirical finding is that the minus sign problem in this approach is severely reduced or maybe even eliminated @xcite .
the expansion in hybridization strength has a long history starting with the famous non - crossing approximation @xcite and various extensions of it such as oca @xcite , ctma @xcite , sunca @xcite .
all these approximations can be viewed as partial summation of the same type of diagrams . with stochastic sampling ,
the summation of essentially all the diagrams is now possible .
the only weakness of the approach is that it works exclusively with the imaginary time green s functions and analytic continuation to real axis and access to real axis self - energy , for example , is still hard to achieved .
however , we believe that the substantially enhanced precision of the method , as compared to hirsh - fye qmc , will make this step easier .
the idea of expanding partition function in terms of the hybridization with the conduction band , dates back to work of g. yuval and p.w .
anderson @xcite . in this work ,
mapping the hybridization expansion to coulomb gas model lead to one of the first breakthroughs in the area of the kondo problem .
generalization to asymmetric anderson model was published by h.d.m .
haldane in ref .  .
similar approach was later used in exploring the physics of the generalized hubbard model in the context of dmft @xcite . in this work ,
the hybridization expansion of the partition function was analyzed by renormalization - group analysis technique .
an early implementation of the related idea to sum up the terms of the partition function by monte carlo sampling was implemented in ref . to solve the two two - channel anderson impurity model .
_ sampled partition function : _ in the new quantum monte carlo method , the partition function is expanded in terms of hybridization strength @xmath10 and the resulting diagrams are summed up by stochastic metropolis sampling .
taylor expansion of the impurity partition function eq .
( [ zimp ] ) gives    @xmath20 e^{-s_{c}}\sum_k \frac{1}{k ! }    \left [      \sum_{{\alpha}{\alpha } ' }      \int_0^\beta d\tau\int_0^\beta d\tau '      \psi_{{\alpha}'}(\tau')\psi^\dagger_{\alpha}(\tau)\delta_{{\alpha}{\alpha}'}(\tau,\tau ' )    \right]^k\end{aligned}\ ] ]    by separating the cluster contribution from the bath contribution the partition function can be cast into the form    @xmath21 e^{-s_{c}}\sum_k \frac{1}{k ! }    \int_0^{\beta}\prod_{i=1}^k d\tau_i \int_0^{\beta}\prod_{i=1}^k d\tau_i '    \sum_{{\alpha}{{\alpha}}'}\prod_{i=1}^k\left[\psi_{{{\alpha}}'_i}(\tau_i')\psi_{{\alpha}_i}^\dagger    ( \tau_i)\right ]    \times \prod_{i=1}^k \delta_{{\alpha}_i{{\alpha}}'_i}(\tau_i,\tau_i ' ) .
\label{z_f}\end{aligned}\ ] ]    it becomes clear that partition function is a product of two terms : the average over the cluster states @xmath22 and average over the bath degrees of freedom @xmath10 .
it was pointed out in ref .
that naive sampling of the above diagrams would run into a very bad minus sign problem .
the reason is that crossing diagrams ( vertex corrections to famous non - crossing approximation to anderson impurity model ) can have either sign and thus weights that correspond to crossing diagrams could be negative .
the ingenious idea proposed in ref .   is to combine all the diagrams of the same order , crossing and non - crossing , into a determinant .
mathematically , this can be expressed by    @xmath23    where @xmath24 e^{-s_{c}}$ ] and average of the operator is @xmath25 e^{-s_{c } } o$ ] .
this is the central equation of the continuous time monte carlo sampling around the atomic limit . to derive eq .
( [ z_sample ] ) from eq .
( [ z_f ] ) one needs to permute time integration variables @xmath26 in all possible ways .
permutation of fermions gives minus sign in an odd permutation .
this minus sign can be absorbed in the minus sign of the product of hybridizations , resulting in the determinant of hybridizations .
_ simulation : _ the set of diagrams , which are associated with the set of imaginary times @xmath27 and corresponding band indexes @xmath28 are visited by monte carlo ( metropolis ) algorithm with the weights given by eq .
( [ z_sample ] ) .
the effect of the hybridization @xmath29 is to create a kink in the time evolution of the cluster , i.e. , to destroy one electron at time @xmath30 on the cluster and create another electron at some other time @xmath26 on the cluster .
the number of kinks is always even due to particle number conservation .
two monte carlo steps which need to be implemented are : i ) insertion of two kinks at random times @xmath31 and @xmath32 ( chosen uniformly @xmath33 ) , corresponding to a random baths @xmath34 , @xmath35 .
ii ) removal of a two kinks by removing one creation operator and one annhilation operator .
many other steps can greatly reduce the sampling time , for example displacing a randomly chosen operator ( either @xmath22 or @xmath36 ) to a new location chosen uniformly @xmath33 .
the double step of inserting or removing two kinks is also possible and is relevant when off - diagonal components of @xmath10 are dominant .
the detailed balance condition requires that the probability to insert two kinks at random times @xmath26 , @xmath30 , being chosen uniformly in the interval @xmath33 , is @xmath37 \label{padd}\ ] ] where @xmath38 is the number of baths , @xmath39 is the current perturbation order ( number of kinks/2 ) , @xmath40 is the cluster matrix element @xmath41 and @xmath42 is the ratio between the new and the old determinant of baths @xmath10 and can be evaluated using usual shermann - morrison formulas .
the factors of @xmath43 enter because of the increase of the phase space when adding a kinks ( increase of entropy ) while the factor @xmath44 comes from factorials in eq .
[ z_sample ] .
similarly , the probability to remove two kinks , chosen randomly between @xmath45 $ ] is @xmath46 .
\label{premove}\ ] ]    an important simplification occurs if the hybridization is block diagonal . since @xmath47 vanishes for some combination of @xmath48 the determinant in eq .
( [ z_sample ] ) can be written as a product of smaller determinants , one for each block of hybridization .
a specially simple case occurs when all the blocks are of size one , the partition function becomes a product of @xmath38 terms , where @xmath38 is number of all baths    @xmath49    the probability to add two new kinks eq .
( [ padd ] ) now depends on the number of kinks of this particular type @xmath50 ( rather than total perturbation order @xmath39 ) and the dimension of the bath subspace @xmath38 becomes unity .
in general the probability to add two kinks is @xmath51 , where @xmath52 is the number of bands which form an off - diagonal sub - block in hybridization @xmath10 and need to be treated simultaneously in one determinant @xmath53 in eq .
( [ z_sample ] ) .
the size of hybridization determinants can thus be greatly reduced if hybridization @xmath10 is block diagonal .
the cluster term @xmath54 however can not be broken into separate contributions for each bath , rather the full trace needs to be computed numerically .
it is therefore essential to find a fast way to compute the cluster average @xmath55 in eq .
( [ z_sample ] ) .    _ exact diagonalization of the cluster : _ it is crucial to evaluate the cluster trace in eigenbase of the cluster hamiltonian and take into account the conservation of various quantum numbers , such as the particle number , the total spin , and the total momentum of the cluster states .
typical contribution to the cluster part of the trace that needs to be evaluated at each monte carlo step , takes the form    @xmath56    where the matrix elements are @xmath57 and @xmath58 are eigenvalues of the cluster . the actual order of operators in eq .
( [ evolution ] ) depends on their time arguments and creation operator is not necessary followed by annhilation operator .
the bottleneck of the approach is that the number of cluster states @xmath59 is very large ( for example , single site dmft for the @xmath0 shell requires 2@xmath60 states ) .
consequently , the matrix elements @xmath61 are in general very large matrices and the typical diagram order is inversely proportional to temperature ( see fig .
[ fig1 ] ) therefore one typically needs to multiply few hundred large matrices at each monte carlo step .
it is clear that this is very impractical and the progress can be achieved only by various trick which we have implemented    * most of matrix elements vanish .
a fast algorithm is needed to determine which matrix elements are nonzero .
* symmetries of the problem can be taken into account to reduce the size of the @xmath61 matrix . *
the number of trial steps is usually much bigger ( 100 times ) than the number of accepted steps and the insertion or removal of a kink is very local in time operation .
it is convenient to store the product eq .
[ evolution ] ( from both sides , left and right ) and when trying to insert a new kink , recompute the trace only between the inserted times @xmath31 and @xmath32 . * during simulation , the probability for visiting any cluster state can be recorded and can be used in the next step to eliminate the irrelevant atomic states from the trace in eq .
[ evolution ] .
the cluster base can hence be adjusted dynamically to describe the particular regime studied by the minimum number of relevant cluster states .    to illustrate the method ,
let us consider a concrete example of the cluster of one band model ( hubbard or t - j ) in the normal state .
the bath index @xmath34 runs over cluster momenta @xmath62 and spin orientation @xmath63 .
the eigenstates of the cluster can be written in a form @xmath64 , where @xmath65 is total number of electrons in the state , @xmath66 and @xmath67 are the total spin and its @xmath68 components , while @xmath69 is the total momentum of the cluster state and @xmath70 stands for the rest of the quantum numbers .
_ concept of superstates : _ in this base , the matrix elements of the creation operator are greatly simplified @xmath71 and @xmath72 .
it is therefore convenient to group together states with the same @xmath73 and treat the rest of the quantum numbers as internal degrees of freedom of a cluster superstate @xmath74 .
the superstate @xmath75 is multidimensional state with internal quantum numbers @xmath76\rangle\equiv|\{s,\gamma\rangle\}$ ] .
it is then clear that creation operator acting on a state @xmath75 gives a unique state @xmath77 and it is enough to store a single index array @xmath78 to figure out how the hilbert subspaces are visited under application of a sequence of creation and annhilation operators such as in eq .
( [ evolution ] ) : @xmath79 .
this sequence is very often truncated in few steps only , because most of the sequences contain either multiple application of the same creation or annhilation operator ( pauli principle ) or because they lead to a state outside the base ( for example @xmath80 or @xmath81 ) .
once the nonzero matrix elements are found by this simple index lookup , the value of the matrix element needs to be computed by matrix multiplication . by this breakup of the hilbert space and
introducing superstates @xmath75 , the largest matrix which needs to be treated for the t - j model on a plaquette is @xmath82 and for the hubbard model it is @xmath83 , thus substantially smaller than the original @xmath84 and @xmath85 dimensional hilbert space . to compute the trace in eq .
( [ evolution ] ) we start with unity matrix in each subspace of a superstate @xmath75 and apply both the time evolution operator @xmath86 ( by multiplying each row of a matrix with its time evolution ) and the kink ( by multiplication with the matrix @xmath87 or @xmath88 ) . the operation of @xmath61 brings us to the next superstate @xmath89 where we repeat both the time evolution and the kink application .
after @xmath39 steps , the trace of a matrix gives the desired matrix element .
_ storing the time evolution : _ the number of kinks is proportional to inverse temperature @xmath90 and kinetic energy of the system @xmath91 ( see eq .
( [ avek ] ) ) .
it thus becomes large at low temperatures .
however , an insertion of a kink with large time difference ( @xmath92 with large @xmath93 ) has a very low probability .
the reason is that pauli principle forbids insertion of the pair @xmath94 if another kink of the same species @xmath95 is between the two times @xmath96 . at the same time
, @xmath97 is like @xmath98 peaked at small times @xmath93 and falls of at large times making the long time intervals rare .
the insertion of a kink is thus fairly local in time operation , therefore it is convenient to store a whole chain of products that appear in eq .
( [ evolution ] ) from both sides , left and right , to make trial step very cheap .
it takes only few matrix multiplications ( almost independent of temperature ) to compute the trace in eq .
( [ evolution ] ) .
when the move is accepted , the trace needs to be updated which takes somewhat more time .
however , the acceptance rate is typically small and on average does not require much computational time .
_ adjusting the cluster base : _ finally , the ultimate speedup can be achieved by dynamically adjusting `` the best '' cluster base .
the probability for a cluster state @xmath59 can be computed during simulation . for a given diagram with particular configuration of kinks ,
the probability for a cluster state @xmath59 is proportional to its matrix element defined in eq .
( [ evolution ] ) , i.e. , @xmath99 the sampled average of this quantity gives probability for cluster eigenstate @xmath59 .
a large number of cluster states have very small probability and can be eliminated in simulation to ultimately speed up the simulation .
it is important that the probability for any cluster state depends on the particular problem at hand and the program adjusts the base dynamically after a few million of monte carlo steps .
_ green s function evaluation : _ like in other impurity solvers which are based on the expansion of the hybridization ( nca , oca , sunca ) , the green s function is computed from the bath electron t matrix
. using equation of motion , it is easy to see that the bath green s function @xmath100 is connected to the local green s function @xmath3 through the following identity @xmath101 or summing over momenta @xmath102 this green s function @xmath103 is equal to the ratio between the determinant of @xmath104 s @xmath105 where one row and one column is added to the bath electron determinant .
the reason for this simple form is that the conduction band is non - interacting and thus obeys wicks theorem .
here we used the definition @xmath106 using shermann morrison formulas to express enlarged determinant by the original determinant of @xmath104 , @xmath103 becomes @xmath107 finally , comparing eq .
( [ gl1 ] ) and ( [ gl2 ] ) we see that @xmath108 and in imaginary frequency @xmath109 this equation is the central equation of the approach since it relates local green s function with the quantities computed in qmc importance sampling .
the equation shows that only matrix @xmath110 needs to be stored and manipulated in simulation .    from the above consideration
, it is clear that one could also add two rows and two columns to matrix @xmath104 and compute the two particle vertex function in similar way without much overhead .
_ fast updates : _ the green s function can be updated in linear time ( rather than quadratic ) .
when adding a construction and annhilation operator at @xmath26 and @xmath30 , adding a column at @xmath111 and row at @xmath112 to matrix @xmath113 , leads to the following relation between the matrix @xmath114 and @xmath115 @xmath116 here one row and one column of zeros is added to @xmath117 to match the size of @xmath118 .
the arrays @xmath119 and @xmath120 are given by @xmath121 where @xmath122 and @xmath123 is @xmath124 finally , the local green s function becomes @xmath125 it is clear from eq .
( [ dgn ] ) that only linear amount of time is needed to update the local green s function .    when removing two kinks of construction and annhilation operators at @xmath111 and @xmath112 , old and new matrix @xmath126 are related by @xmath127 green s function therefore becomes @xmath128    finally , the exponential factors @xmath129 do not need to be recomputed at each monte carlo step since all `` old '' times can be stored and the exponents need to be computed only for the new pair of times and only at each accepted move . in the present implementation , the algorithm to sample directly @xmath130 is sufficiently fast that it does not introduce any performance costs .
since it does not introduce systematic error in binning @xmath98 we believe it is superior than the alternative implementations which sample @xmath98 .
_ large frequencies and moments : _ similarly to the hirsch - fye qmc , the low frequency points of green s function converge very fast to the exact value while the high - frequency points , when sampled directly , contain a lot of noise .
it is therefore not very useful to sample large frequencies in the above described way .
usually we sample 200 - 300 matsubara points while the rest are replaced by the high frequency moments of the self - energy computed analytically .
the high frequency moments of the self - energy are computed from the green s function moments , which , in general , take the following form @xmath131\cdots ] ] , \psi^\dagger_{\beta}\right\}\right\rangle\ ] ] to compute the moments within the present approach , few operators need to be sampled in simulation . in the one - band model , only density is required , but in more complicated situation , higher order density - density and exchange terms enter . in general , green s function moments can be expressed by the average of a few equal time operator @xmath132 , which take particular simple form in the local eigenbase @xmath133 the lowest order self - energy moments can then be computed in the following way @xmath134    _ sampling of other quantities : _ equal time operator like density or those in eq .
( [ gmoments1 ] ) and ( [ gmoments2 ] ) can be straightforwardly computed in simulation . to improve the precision of simulation ,
operators are averaged over all times , i.e. , @xmath135\ ] ] the contribution to @xmath136 of each particular diagram is @xmath137 \label{o1}\ ] ] where @xmath138 is the cluster part of the trace eq .
( [ evolution ] ) , @xmath139 are the operators in eq .
( [ evolution ] ) which appear before @xmath140 and @xmath141 contains operator with time arguments greater than @xmath140 .
@xmath39 is the number of kinks in the diagram @xmath142 .
as we mentioned before , these time evolution operators ( @xmath139 , @xmath143 ) are stored and regularly updated during simulation .
further simplification is possible , if operator @xmath136 commutes with @xmath144 @xmath145.\ ] ] if operator @xmath146 is one of the good quantum numbers of a superstate ( one of @xmath73 in our example above ) , matrix @xmath147 is proportional to identity matrix in each superstate @xmath148,m_2[i]}=\delta_{m_1,m_2 } o_i$ ] .
starting with a superstate @xmath149 , kinks of a diagram @xmath142 generate a sequence of superstates @xmath150 , as explained above .
the equal time average of cluster conserved quantity therefore simplifies to @xmath151 where @xmath152 is probability for a superstate @xmath75 defined in eq .
( [ pm ] ) and @xmath153 is the value of the conserved quantity in each superstate @xmath154 in the interval @xmath155 .
this time difference is the time each superstate `` lives '' in the particular diagram .    in this way
, the electron density @xmath156 and average magnetization @xmath67 is computed to very high precision .
similar simplification in evaluating the total spin susceptibility , defined by @xmath157 where @xmath67 is the total spin of the cluster ( or atom ) leads to @xmath158 where @xmath159 is the total spin of the cluster in cluster eigenstate @xmath160 .
if @xmath146 does not commute with @xmath161 , the integral ( [ o1 ] ) needs to be evaluated . in cluster eigenbase
it takes the form @xmath162    _ sampling the total energy : _ the average of the potential energy @xmath163 can be computed from the average of the local energy since @xmath164\end{aligned}\ ] ] here we concentrate on the case of general single site dmft rather than the cluster extensions .
the reason is that kinetic energy in cluster extensions depends on the periodization scheme and we will not go into this details here .
kinetic energy of the general single site dmft @xmath165 $ ] can be computed by @xmath166 + { \mathrm{tr}}[(\mu+e_{imp})n]\end{aligned}\ ] ] the total energy is therefore given by @xmath167     + \mu n\end{aligned}\ ] ] the first term @xmath168 can be computed very precisely in simulation .
the sampled quantity @xmath136 is just the energy of an atomic state and can be simply obtained from the probabilities for atomic states @xmath169 .
computing kinetic energy from the green s function gives in general worse accuracy because the high - frequency behavior of the green s function can not be directly sampled and augmentation with analytically computed tails is necessary .
however , it is simple to show that the average value of the perturbation order is related to the average of the kinetic energy @xmath170 \label{avek}\ ] ] where @xmath171 is the average perturbation order and @xmath172 is temperature . the later quantity is directly sampled in the present algorithm and it is just the center of gravity of the histogram ( presented in fig .
[ fig1 ] ) .
finally , the total energy @xmath173 is given by @xmath174 all quantities in this equation can be computed to very high accuracy and since low temperatures can be reachable in this method , the entropy can be obtained by integrating the specific heat .    _ superconductivity : _ the power of the method can be further demonstrated by studying the superconducting state of the strongly correlated systems at low temperature with essentially no performance cost . by employing the nambu formalism , the translationally invariant cluster methods ( for details on cellular - dmft on a plaquette see appendix [ appa ] ) results in @xmath175 two dimensional baths , where @xmath175 is the number of cluster momenta .
namely , the baths @xmath176 and @xmath177 are coupled through the anomalous component of hybridization and require simultaneous treatment in @xmath178 in eq .
( [ z1b1 ] ) .
the determinants are on average twice as large as in normal state , however , the cluster part of the trace in eq .
( [ z1b1 ] ) remains unchanged . since most of the time is
usually spend in evaluating the local part of @xmath179 , the performance is not noticeably degraded in superconducting state . in typical run presented below , the histogram is peaked around @xmath180 which is equal to the order of a typical diagram . in the translational invariant representation employed here , the size of a typical determinant in eq .
( [ z1b1 ] ) is only @xmath181 and using fast - update scheme presented above , the trace over the bath states is not expensive part of the algorithm .
_ hund s coupling and spin - orbit coupling : _ in materials with open @xmath0 orbitals , the multiplet effects are very strong and @xmath182 approximation is not adequate .
simultaneous inclusion of hund s coupling and spin - orbit coupling in dmft method is crucial for quantitative description of actinides @xcite .
minimal local hamiltonian for lanthanide and actinide materials is @xmath183 here @xmath184 is the impurity level without the spin - orbit coupling since later is included explicitely .
the hund s coupling and spin - orbit coupling take the following form    @xmath185 c^{j m_j}_{lm,\sigma}c^{j m_j}_{l m',\sigma ' } f_{lm\sigma}^\dagger f_{lm'\sigma'}\end{aligned}\ ] ]    here @xmath186 are spheric harmonics , @xmath187 are clebsch - gordan coefficients and @xmath188 are slater integrals .
the @xmath189 , usual hubbard @xmath1 , is commonly computed by constrained lda , while the rest of the slater integrals ( @xmath190 , @xmath191 and @xmath192 ) are computed using atomic physics program @xcite . finally , the spin - orbit strength @xmath193 is computed within lda program and needs to be updated during self - consistent lda+dmft calculation .
exact diagonalization of the atomic hamiltonian leads to eigenstates with conserved number of particles @xmath65 , total angular momentum @xmath194 and its @xmath68 component @xmath195 .
for the one electron base , it is more convenient to chose the electron total angular momentum @xmath196 rather than @xmath197 since the construction and annhilation operators conserve @xmath68 component of total spin @xmath198 superstates @xmath75 are chosen in a way to simplify the matrix elements of the creation operator @xmath36 . in this case ,
a convenient choice is @xmath199 in the absence of crystal field splitting , the hybridization @xmath10 is a diagonal matrix ( all @xmath200 baths are one - dimensional ) .
however , strong crystal fields splitting generates off - diagonal components of @xmath201 and in general , the determinant @xmath178 in eq .
( [ z_sample ] ) can not be broken into small determinants as in eq .
( [ z1b1 ] ) .
we implemented cellular dmft for both the hubbard and the t - j model on a plaquette .
the energy scales are given in units of @xmath202 , hubbard @xmath1 is fixed at @xmath203 while @xmath194 of the t - j model is set to @xmath204 .
we simulated @xmath205 monte carlo steps per processor and results were averaged over 64 processors .
one dmft step for the hubbard model takes approximately 45 minutes on 1.7@xmath206ghz pc processor therefore each monte carlo step requires on average one million clock - ticks ( 0.9 mflops ) .    .
]    fig .
[ fig1 ] shows histograms ( probability distribution for the perturbation order @xmath39 ) for few dopings of the hubbard model @xmath207 at @xmath208 .
the average perturbation order is increasing with doping since the electrons are getting more delocalized ( absolute value of the kinetic energy is increasing ) and the creation of kinks becomes less expensive .
it is of the order of 450 in overdoped regime , but the typical size of the determinant is only 112 in superconducting state or 56 in normal state .
we recently addressed the problem of coherence scale in the t - j model @xcite and we found , using nca as the impurity solver , that the imaginary part of the cluster self - energy @xmath209 , which plays the crucial role in the optical conductivity and transport , becomes very large at optimal doping and consequently the coherence scale vanishes around optimal doping . here
we extend this study to the hubbard model using much lower temperature .
we will show that the system become strongly incoherent at optimal doping and the maximum of @xmath210 tracks the maximum scattering rate in both the hubbard and the t - j model .
[ fig2 ] shows the imaginary part of the self - energy @xmath211 at few different dopings and temperature @xmath212 which is around the superconducting critical temperature of this approach .
the system is still in normal state .
it is clear that the self - energy at large frequencies is a monotonic function of doping and is largest at the mott transition @xmath213 .
however , the low frequency region is distinctly different and the crossing of self - energies is observed around @xmath214 .
low doping as well as large doping self - energies can be extrapolated to zero , while at optimal doping self - energy remains of the order of unity even around the critical temperature .     on imaginary axis .
temperature @xmath212 is around the critical temperature but still in normal state . at low frequency
, the most incoherent self - energy corresponds to optimal doped and not to the underdoped system . ]
self - energy on the imaginary axis for both the t - j ( left ) and hubbard ( right ) model showing the large scattering rate at optimal doping .
the lower panel shows the anomalous self - energy for the same models and doping levels and can be used to locate the optimal doped regime . ]
the upper panels of figure [ fig3 ] show the same self - energy at small frequency for both the t - j ( left ) and the hubbard ( right ) model . in this figure , temperature is @xmath215 and is far below @xmath210 therefore all curves extrapolate to zero and the system becomes coherent in superconducting state .
however , the reminiscent of the strong incoherence at optimal doping is the large slope of imaginary part of @xmath216 which induces very small quasiparticle residue in this regime .
the precise position of the optimal doping , characterized by the largest low frequency anomalous self - energy , is different in the t - j and in the hubbard model ( see lower panel of fig .
[ fig3 ] ) .
it is around @xmath217 in the former and around @xmath218 in the later .
however , the large scattering rate is ultimately connected with the largest anomalous self - energy and hence largest critical temperature in both the t - j and the hubbard model . in our view
, the most important advantage of the new quantum monte carlo method is that it can treat the realistic multiplet structure of the atom .
it was shown in ref .   that hubbard term only ( @xmath219 ) ,
leads to sever underestimation of the interaction strength in actinides and for realistic values of @xmath219 , dmft predicts heavy fermion state in curium rather than magnetic state .
negligence of the multiplet structure of plutonium misses the fine structure of the quasiparticles ( two peaks around 0.5 and 0.85 ev ) and more importantly , predicts only weakly correlated metallic state in delta phase of plutonium .    to demonstrate the advantage of the new method , we recomputed the localization - delocalization transition in the archetype material exhibiting kondo collapse , namely the @xmath220 transition of elemental cerium .
since the number of electrons in cerium fluctuates between states with zero , one and two electrons in the @xmath0-shell , the number of atomic states that needs to be kept is relatively small hence solving the impurity problem requires very little computational power in this case . in fig .
[ fig4 ] we show the `` valence histogram '' @xcite of the two phases of cerium , i.e. , the projection of the density matrix to the eigenstates of the atom .
the plot shows the probability to find an @xmath0 electron of cerium in any of the atomic eigenstates and demonstrates how strongly is the atom fluctuating between atomic states .
the typical fluctuating time is inversely proportional to the kondo temperature of the phase , being around 2000@xmath221 in the alpha phase and around 80@xmath206k in the gamma phase .
the itinerant alpha phase histogram is peaked for many atomic states , including the spin - orbit split @xmath222 and @xmath223 singly occupied states as well as the empty state . on the other hand ,
the local - moment gamma phase is peaked only at the ground state of the singly occupied sector with @xmath222 spin showing that the dmft ground state closely resembles the atomic @xmath224 ground state .
electron in the solid spends appreciable time in a few atomic configurations .
the height of the peak corresponds to the fraction of the time the @xmath0 electron of the solid spends in one of the eigenstates of the atom , denoted by the total spin @xmath194 of the atom .
we summed up the probabilities for the atomic states which differ only in the z component of the total spin @xmath225 .
the @xmath226 axis indicates the energy of atomic eigenstates in the following way : @xmath227 and @xmath228 , where n is between 0 and 2 . ]
we generalized the recently developed continuous time quantum monte carlo expansion around the atomic limit @xcite to clusters treated within cellular dmft or dynamical cluster approximation as well as to materials that require realistic hund s coupling and multiplet splitting of the atomic state .
we explained the steps necessary for the efficient implementation of the method .
our low temperature data in the strongly correlated regime of the hubbard and t - j model show the efficiency of the method and demonstrate its superiority compared to the conventional hirsh - fye quantum monte carlo method .
the long - standing problem of adequate treatment of the multiplet splitting within dmft is resolved .
this splitting is crucial in actinides and its omission can lead to wrong prediction of the magnetic nature of the dmft ground state solution .
we showed that the optimal doped regime in both the t - j and the hubbard model is characterized by the largest scattering rate and the system becomes more coherent in the underdoped and overdoped regime .
the precise position of the optimal doping , as determined from the maximum in anomalous self - energy , is different in the two models .
however , the strong incoherence is always found at the doping corresponding to maximum @xmath229 .
we computed the valence histogram across the cerium alpha to gamma transition with emphasis on the multiplet splitting of the atomic states .
we showed that the atom fluctuates between many atomic states in itinerant alpha phase and both the @xmath222 and @xmath223 spin - orbit split states have large probability in the ground state of the system .
the empty state and the doubly occupied states , which are substantially split by hund s coupling , acquire a finite probability .
the gamma phase , on the other hand , shows well defined valence @xmath230 and charge fluctuations become rare .
the author wish to thank g. kotliar for very stimulating discussion and careful reading of the manuscript .
for the present continuous time monte carlo method , it is convenient to pick the base such that the local quantities are diagonal .
cellular dmft uses a generalized open boundary conditions and therefore , in general , momentum is not a good quantum number .
however , for small clusters , such as plaquette , all sites are equivalent and translational invariance is still obeyed .
the local green s function and hybridization in momentum base take the diagonal form @xmath231 in superconducting state @xmath232 is a @xmath233 matrix in nambu notation .
although the local quantities are diagonal , non - interacting hamiltonian at general @xmath234-points is not .
the self - consistency condition in superconducting state takes the following form      where we defined @xmath236 and assumed @xmath237 here @xmath238 is the normal on - site self - energy , @xmath239 is the nearest - neighbor and @xmath240 is the next - nearest neighbor self - energy and @xmath241 s are the anomalous components of the self - energy . for d - wave symmetry ,
@xmath242 and @xmath243 vanish and @xmath244 .
the advantage of this formulation of the cellular dmft is that the hybridization becomes block diagonal and hence the determinants ( @xmath245 ) which enter eq .
( [ z_sample ] ) can be broken up into separate contribution for each momentum @xmath69 point like in eq .
( [ z1b1 ] ) .
a. georges , g. kotliar , w. krauth and m.j .
rozenberg , rev .
phys . * 68 * , 13 ( 1996 ) .
g. kotliar , and d. vollhardt , physics today * 57 * , 53 ( 2004 ) .
g. kotliar , s.y .
savrasov , k. haule , v.s .
oudovenko , o. parcollet , c.a .
marianetti , rev .
mod . phys . * 78 * , 865 ( 2006 ) .
savrasov , g. kotliar , e. abrahams , nature * 410 * , 793 ( 2001 ) .
savrasov , k. haule , and g. kotliar , phys .
lett . * 96 * , 036404 ( 2006 ) .
x. dai , s.y .
savrasov , g. kotliar , a. migliori , h. ledbetter , e. abrahams , science * 300 * , 5621 , 953 - 955 ( 2003 ) . k. haule , v. oudovenko , s.y .
savrasov , and g. kotliar , phys .
. lett . * 94 * , 036401 ( 2005 ) . j.h .
shim , k. haule , and g. kotliar , cond - mat/0611760 .
oudovenko , g. palsson , k. haule , g. kotliar , and s.y .
savrasov , phys .
b * 73 * , 035120 ( 2006 ) .
h. kajueter and g. kotliar , phys .
lett . * 77 * , 131 ( 1996 ) .
r. bulla , phys .
lett . * 83 * , 136 ( 1999 ) ; r. bulla , t. a. costi , and d. vollhardt , phys .
b * 64 * , 045103 ( 2001 ) .
x. dai , k. haule , and g. kotliar , phys . rev .
b * 72 * , 045111 ( 2005 ) .