the qcd sum rule ( qcd - sr ) approach to qcd continues to be a highly active field .
the influence of the field is reflected in over 1500 references to the seminal paper of shifman , vainshtein and zakharov ( svz ) @xcite , with over 380 of these references following the beginning of 1992 @xcite .
while the approach has been applied to a variety of hadronic observables in both vacuum and finite density nuclear matter , a comprehensive review of the systematic errors associated with the approach has not been considered .
one of the key assumptions of the approach is the use of the so called `` continuum model '' to remove excited state contaminations from the hadron correlator under investigation .
recently , the continuum model for the nucleon was tested by using correlation functions calculated via lattice regularized qcd @xcite .
the results suggest that the continuum model adequately removes the excited state contaminations and allows the isolation of the ground state .
the success of the continuum model warrants a more careful investigation of other aspects of the qcd - sr approach so that uncertainties in the predictions may be reliably and quantitatively determined .
the focus of this paper is to establish a rigorous procedure for extracting quantities of phenomenological interest from qcd sum rules .
the main conclusion of this investigation is that qcd sum rules work and they are predictive when the analysis is done rigorously .
the ground state masses of the @xmath0-meson and the nucleon are used to introduce the concepts .
however , the procedures introduced here are more generally applicable to quantities of current experimental interest .
the goal is to shed light on the stability of qcd sum rule analyses and thus the accuracy to which observables may be reliably determined .
implementation of the procedure introduced here does not lead to slight adjustments in the extracted parameters , but rather leads to qualitatively different conclusions .
it will become apparent that most existing qcd - sr investigations of nucleon properties are in fact unreliable .
previous calculations are based on sum rules in which neither the operator product expansion ( ope ) nor the phenomenological description are under acceptable control .
the qcd - sr method is presently the best fundamentally based approach for investigating the properties of hadrons in nuclear matter , as the lattice approach is challenged by a number of formidable obstacles @xcite .
issues surrounding the finite width of the rho - meson make reliable lattice calculations of @xmath2 mixing very difficult .
these are two of the many topics of current interest to nuclear physicists where the best fundamentally based approach is that of qcd - srs .
moreover , the approach has minimal model dependence .
hence it is imperative to probe the predictive ability of the approach and the ideas presented here are currently being applied to these topics mentioned above @xcite .
the most significant source of uncertainty in the qcd sum rule approach is an imprecise knowledge of the vacuum condensates appearing in the operator product expansion ( ope ) .
this is particularly problematic for the higher dimension operators , where one usually invokes the vacuum saturation hypothesis ( factorization ) to replace higher - dimension - operator vacuum - expectation values by products of lower - dimension operators .
a list of key sources of uncertainty in the qcd - sr approach should include :    * the unfactorized condensate values , * factorization of higher dimensional operators , * small but neglected @xmath3 corrections , * the truncation of the operator product expansion ( ope ) , * the selection of the regime for matching the qcd and phenomenological sides of the sum rules , * uncertainties associated with the summation of perturbation theory in large orders , and * the possibility of significant direct instanton contributions .    with the exception of direct instanton contributions ,
the effects of these uncertainties will be estimated via a monte - carlo error analysis .
gaussian distributions for the condensate values are generated via monte carlo .
these distributions are selected to reflect the spread of values assumed in previously published qcd - sr analyses , and the uncertainties such as operator factorization listed above .
these distributions provide a distribution for the ope and thus uncertainty estimates for the ope which will be used in a @xmath4 fit .
in turn , the ope distribution provides distributions for the phenomenological fit parameters , thus establishing the predictive ability of the qcd - sr approach .      in the field of qcd inspired models
, the focus is often on whether or not the model is able to encompass the known experimental data , as opposed to assessing the predictive potential of the method .
a review of the qcd - sr literature leaves no doubt as to whether or not the qcd - sr method is able to encompass the established data . the poorly determined condensate values and arbitrary aspect of the matching borel regime provides tremendous freedom in devising ways to achieve the desired result
in fact , it is difficult to find qcd - sr - based predictions of hadronic observables which fail to agree with experiment .
this is somewhat surprising , since the physics responsible for some observables may be poorly represented in a truncated ope commonly limited to distances of less than 0.3 fm .
most qcd - sr calculations of nucleon properties in the literature fall short in the analysis stage where the phenomenological parameters are determined by matching the qcd and phenomenological sides of the sum rules .
typical analysis short - comings include :    * selecting a single value for the borel parameter which gives `` nice '' results .
a borel _ regime _ should be selected in order to evaluate the stability and reliability of the results . *
a selection of the borel regime without careful regard to ope convergence nor maintaining ground state dominance of the phenomenological side of the sum rules . * the fixing of search parameters ( such as the continuum threshold ) to preferred values .
it will become apparent that this introduces a strong bias to the remaining fit parameters which may not reflect the properties of qcd .
* claiming an accuracy for qcd - sr predictions without supporting calculations .
occasionally a `` stability analysis '' @xcite is considered , in which fit parameters are monitored as a single condensate value is varied .
however , such analyses explore a relatively small corner of the condensate parameter space .    in general , former procedures adopted for matching the two sides of the sum rules lack the level of rigor which will be presented here .
this general lack of rigor can be traced back to the original paper of svz @xcite . in their conclusions , they comment `` we prefer not to deepen into the debris of computer calculations , sophisticated fit programs , arranging error bars here and there@xmath5 '' while this may be appropriate for a seminal paper on the qcd - sr approach , continuation of this philosophy has lead to strong negative criticism of the qcd - sr field .
it is not uncommon to hear remarks such as `` qcd sum rules are useful only if you know the answer . '' or `` you can get anything you want from qcd sum rules . '' and `` 10% here , 10% there  pretty soon you re talking about real numbers ! ''
these remarks are unfortunate , as a rigorous and proficuous analysis is possible .      in this paper ,
a method for the quantitative determination of phenomenological quantities with uncertainties will be presented . since this paper
is directed to both qcd - sum - rule practitioners and skeptics alike , we begin by briefly outlining the qcd - sr approach in section [ qcdsr ] .
a more pedagogical review of the approach may be found in the recent review article of ref .
this section also introduces the `` ioffe formula '' @xcite , argued by many to qualitatively encapsulate a description of the nucleon mass in terms of the chiral symmetry breaking order parameter @xmath1 .
uncertainties associated with the evaluation of the ope will be addressed via a new monte - carlo based procedure .
estimates for these uncertainties are presented in section [ uncert ] .
section [ analysis ] details the monte - carlo uncertainty analysis and discusses the optimization algorithm and cautious convergence criteria used in the following .
section [ rho ] presents an analysis of @xmath0-meson sum rules . here
the criteria for the selection of the regime for matching the qcd and phenomenological sides of the sum rules is reviewed .
reasonable alternatives are considered via monte carlo .
these considerations lead to the selection of the optimal interpolating field for nucleon sum rules .
new insights into unconventional nucleon interpolators obtained from the lattice qcd investigation of ref .
@xcite play a significant role here .
section [ nucleon11 ] reviews the selection of the optimal nucleon interpolating field for correlators obtained from spin-1/2 interpolating fields .
the effects of choosing non - optimal interpolators is demonstrated .
section [ nucleon13 ] explores additional sum rules obtained from the overlap of the generalized spin-1/2 and a spin-3/2 interpolator .    in section [ correl ] we explore the correlations between condensate values , fit parameters and sum rule consistency .
these correlations serve to identify the roles of the leading terms of the ope in hadronic physics .
the validity of the `` ioffe formula '' is evaluated here , and these results may be of interest to those modeling the qcd vacuum .
the contingency table analysis reveals that the intimate relationship between the quark condensate and the nucleon mass suggested in the `` ioffe formula '' is invalid .
the necessity of direct instanton contributions to the qcd sum rules is determined in section [ instantons ] .
we demonstrate that there is no evidence indicating that direct instanton contributions are required to maintain sum rule consistency .
the implications of the results of this analysis on coordinate space correlators are discussed in section [ coordimpl ] .
finally , section [ conclusions ] summarizes the findings of this in - depth investigation .
the underlying principle of field theoretic approaches to hadron phenomenology is the ansatz of duality .
that is , it is possible to simultaneously describe a hadron as quarks propagating in the qcd vacuum , and as a phenomenological field with the appropriate quantum numbers .
hadron masses are extracted from an analysis of the two - point function @xmath6 the interpolating field @xmath7 is typically constructed from quark field operators combined to give the quantum numbers of the hadron under investigation .
the @xmath0-meson interpolator related to the decay constant is @xmath8 to maintain maximal overlap with the ground state relative to excited states , only interpolators without derivatives are considered . for the nucleon we begin by considering the most general spin-1/2 nucleon interpolator @xmath9 where    @xmath10    and @xmath11
@xmath12 is the interpolator typically used in lattice qcd analyses .
@xmath13 vanishes in the nonrelativistic limit .
however in a theory with light relativistic current quarks , there is no reason to exclude such an interpolating field _ a priori _ @xcite . with the use of the fierz relations
, the combination of the above two interpolating fields for @xmath14 may be written    @xmath15    [ chisr ]    giving the proton interpolating field advocated by ioffe @xcite and often found in qcd - sr calculations .      at the phenomenological level one
proceeds by inserting a complete set of eigenstates with the quantum numbers of the interpolator .
the ability of the interpolator to annihilate a positive parity baryon to the qcd vacuum is described by the parameter @xmath16 in @xmath17 for state @xmath18 of momentum @xmath19 and spin @xmath20 .
@xmath21 is a dirac spinor .
negative parity states require a @xmath22 on the right - hand side of ( [ lambdadef ] ) .
for the nucleon , the two point function has the form @xmath23 where @xmath24 corresponds to positive / negative parity baryons . while there is some suppression of excited state contributions to the two - point function ,
there is little hope of isolating ground state properties from such a function .
hence one uses the celebrated borel transform @xmath25 = \lim_{{\scriptstyle -p^2 \to \infty \atop                          \scriptstyle     n \to \infty } \atop                          \scriptstyle -p^2/n = m^2 } { 1\over n ! } ( -p^2)^{n+1 } \left ( { d \over d p^2 } \right ) ^n f ( p^2 ) \ , , \ ] ] at each dirac-@xmath26 structure to obtain exponential suppression of excited states . the parameter @xmath27 is commonly referred to as the borel mass .
the borel transform is easily applied to a dispersion relation for @xmath28 at each dirac-@xmath26 structure @xmath29 and the polynomial subtraction terms are eliminated by the borel transform . at the structure @xmath30 one has    @xmath31    where we have introduced the spectral density @xmath32 , @xmath33[gp ]    at the structure @xmath34    @xmath35    where @xmath36    here , @xmath37 and @xmath38 account for excited state contributions , and @xmath24 in ( [ rho1 ] ) corresponds to positive / negative parity states .
for the vector meson one has    @xmath39    @xmath40    [ rhov ]    at the structure @xmath41 .      at the quark level ,
one exploits the operator product expansion ( ope ) to describe the short distance behavior of the two - point function    @xmath42    [ ope ]    where @xmath43 is the normalization point at which the coefficient functions and the operators are defined .
here we have explicitly included all operators up to dimension eight , to leading order in the quark mass @xmath44 , having the quantum numbers of the vacuum .
the first digit of the subscript of the wilson coefficients @xmath45 gives the energy dimension of the operator .
@xmath46 may take any of the 16 independent dirac-@xmath26 matrices .    upon fourier transforming to momentum space , ( [ ope ] )
becomes an expansion in @xmath47 , valid for large momentum transfers .
thus one may use perturbation theory to calculate the wilson coefficients .
the contribution of the gluon operator term 4.1 in the expansion for the proton is suppressed due a factor of @xmath48 which appears in the denominator of the wilson coefficient arising from integration over the quark - loop momentum which introduces the factor @xmath49 .
terms involving more than two gluon field strength tensors @xmath50 and products of @xmath44 and @xmath51 are estimated to be small and are typically neglected .
the standard treatment of the ope proceeds via inserting the explicit forms of the interpolating fields into the two - point function of ( [ twopt ] ) and contracting out pairs of time - ordered quark - field operators which are the fully interacting quark propagators of qcd .
wick s theorem for the time - ordered product of two fermion fields provides @xmath52 taylor expansion of the normal ordered piece leads to @xmath53 the well known trick is to select the coordinate gauge @xmath54 such that @xmath55 and the partial derivatives of ( [ taylor ] ) may be replaced by covariant derivatives @xmath56 generally , one works with ( [ aexp ] ) to leading order .
the following relations allow an easy determination of the covariant derivatives    @xmath57    @xmath58 = i g_c g^{\mu \nu } \ , , \ ] ]    @xmath59    in summary , the quark correlators employed in calculating the ope under the usual assumption of vacuum saturation of the intermediate states of composite operators include @xmath60 \bigl ( g_c g_{\alpha \beta}^n \bigr ) \sigma^{\alpha \beta } { \lambda^n_{ab } \over 2 } \ , .
\label{quarkcorr } \\ & & \nonumber
\\ \lefteqn{\langle \omega \bigm | t \left \ { q^a ( x ) \ , g_c   g_{\alpha \beta}^n \ , \overline q^b ( 0 ) \right \ } \bigm | \omega \rangle = } \hspace{36pt } \nonumber \\ & & - { 1 \over 2 ^ 6 \ , 3 } \langle \overline q g_c \sigma \cdot g q \rangle \sigma_{\alpha \beta } { \lambda^n_{ab } \over 2 }   \nonumber \\ & & + { i \over 2 ^ 8 \ , 3 } m_q \langle \overline q g_c \sigma \cdot g q \rangle   \bigl (   \gamma \cdot x \ , \sigma_{\alpha \beta } + \sigma_{\alpha \beta }   \ , \gamma \cdot
x   \bigr ) { \lambda^n_{ab } \over 2 } \nonumber \\ & & + { x^2 \over 2^{10 } \ , 3 ^ 2 } \langle g_c^2 g^2 \rangle   \langle \overline q q \rangle \sigma_{\alpha \beta } { \lambda^n_{ab } \over 2 } \ , .   \label{quarkgluoncorr}\end{aligned}\ ] ]    finally the momentum space correlator is borel transformed
. here one encounters the factorial suppression of higher - dimension operators as indicated in @xmath61 = ( -1)^k \ , { 1 \over ( k-1 ) ! } \ , { 1 \over ( m^2)^{k-1 } } \ , .\ ] ] it is worth noting that the factorial suppression does not really set in until one reaches the term @xmath62 .
for the @xmath0 meson this term corresponds to an operator of dimension eight , and for the nucleon the dimensions are twelve and thirteen for the two nucleon sum rules obtained from spin-1/2 interpolating fields .
since the ope for the @xmath0-meson is traditionally truncated at dimension eight @xcite , the correlator should be more reliable than for the nucleon sum rules which are truncated at dimension nine .
it is well known that the standard treatment of the ope encounters difficulties at large orders in the perturbative expansion of the wilson coefficients .
perturbative expansions are divergent at large orders and are asymptotic at best @xcite . in the standard treatment , these divergences limit the accuracy to which the perturbative coefficients may be determined .
renormalons @xcite , a particular set of perturbative graphs , provide a simple illustration of the factorial growth in the contributions of perturbative graphs at large orders @xcite .
the factorial growth in large - order contributions associated with soft infrared ( ir ) virtual momenta is problematic , as the series is not borel summable . as such , the wilson coefficients , by themselves ,
are ill defined in the standard treatment . the uncertainty resulting from the restriction to finite orders of the asymptotic expansion
is commonly referred to as the ir renormalon ambiguity . of course the use of perturbative propagators to describe soft virtual momenta is incorrect .
instead , one must utilize fully dressed nonperturbative propagators . hence
, it is not surprising one encounters difficulties in the ir regime , as one is simply encountering the landau pole of the perturbative coupling constant .
a solution to the ir renormalon ambiguity problem was provided by the itep group @xcite through the introduction of the normalization / separation scale @xmath43 of ( [ opemu ] ) .
the ope becomes a separation of scales , with virtual momenta lying above @xmath43 represented in the wilson coefficients , and momenta below @xmath43 in the vacuum expectation values ( vevs ) of higher - dimension operators beginning at dimension four .
the introduction of the scale @xmath43 eliminates the ir renormalon ambiguity of the perturbative expansion .
however , in principle , both the wilson coefficients and the vevs of operators contain perturbative as well as nonperturbative physics .    in principle , the ope defined in ( [ opemu ] ) is free of ir renormalon ambiguities .
however , in practice it is difficult to implement the formal procedure of truncating virtual momenta below @xmath43 when calculating higher order @xmath3 corrections to the wilson coefficients . while these contributions are small for leading order corrections , their inclusion is a source of error which can be accommodated in the monte - carlo uncertainty analysis .    in the ultraviolet ( uv ) virtual momenta regime , the renormalon series is sign alternating and borel summable .
summation of the uv series can give rise to a power correction proportional to @xmath63 which closely resembles the squared quark mass contribution to the ope @xcite .
for the vector correlator , the uv renormalon gives rise to an effective mass of roughly @xmath64 , which is much larger than the usual current quark mass the order of 5 mev @xcite .
hence there is some concern that important power corrections associated with uv renormalons has been overlooked in the standard qcd - sr treatment .
moreover , it is the power corrections that provide the link to nonperturbative phenomena . to the extent that the power correction is scheme dependent and
not reliably known , we will refer to these quark - mass like power corrections as uv renormalon uncertainties .
additional evidence for @xmath65-type power corrections independent of renormalons has been obtained through a consideration of constraints imposed by asymptotic freedom and analyticity on the large - order behavior of perturbation theory @xcite . in this case
the perturbative series is not borel summable .
given the absence of a local gauge - invariant field - operator product of dimension two , some might argue that the coefficient of the singularity in the borel plane governing the large - order behavior of perturbative series must have a zero .
however , perturbative estimates of the coefficient do not provide evidence of such a zero @xcite and the issue is unsettled .
hence the issue of possible dimension - two power corrections is of current interest .
the techniques presented in this manuscript might be used to provide some phenomenological insight into the debate . however
, such considerations would take us too far afield and will be deferred to a subsequent analysis @xcite . on the other hand ,
it is interesting to look for discrepancies which might be rooted in such power corrections , and this is done in the following analysis .      to summarize the points of this section , we present the traditional qcd sum rules for the nucleon @xcite obtained from the consideration of the interpolating field @xmath66 of ( [ chisr ] ) to operator dimension 8 .
@xmath67    [ tradit ]    the dots allow for excited state contributions . )
has had an elusive history .
early calculations @xcite appear to have calculated a subset of the diagrams contributing to the dimension 7 operator @xcite .
however a more persistent error lies in a factor of 3 correction in the term @xmath68
\left < \overline q q \right > \left < g_c^2 g^2 \right > \sigma_{\alpha \beta } \lambda_{ab}^n/2 $ ] of ( [ quarkgluoncorr ] ) from that used in previous analyses
. ]    the ioffe formula for the nucleon mass is obtained from ( [ tradit ] ) by keeping the leading terms of the opes and assuming dominance of the ground state on the right - hand side of the equations ,    @xmath69    [ bigassumption ]    the ioffe formula follows from the ratio of these two equations ,    @xmath70    if one prefers to eliminate the borel mass by selecting @xmath71 , @xmath72    here the quark condensate is estimated from the partially conserved axial current ( pcac ) relation and pion decay to give @xmath73 which corresponds to @xmath74 gev , @xmath75 gev and the average light quark mass is 7 mev .
the qualitative features of chiral symmetry breaking coupled with the quantitative success of these relations has lead to enormous faith being placed into the validity of the ioffe formula .
most practitioners expect this relation to be valid , at least in a qualitative manner of proportionality .
the operator product expansion provides knowledge of the two - point function ( [ twopt ] ) in the near perturbative regime .
as such the correlator represents contributions from both the ground state hadron and excited states . to maintain some predictive ability in the qcd - sr approach , a `` continuum model '' for excited state contributions to the correlator is introduced
guided by the principle of duality , the leading terms of the ope surviving in the limit @xmath76 are used to model the form of the spectral density accounting for excited state contributions . consider , as an example , the identity operator for the nucleon interpolator @xmath77 of ( [ chisr ] ) . explicitly including the excited state term of ( [ gp ] ) in equation ( [ ioffe1a ] )
leads to @xmath78 for large borel mass .
the laplace transform readily provides @xmath79 the continuum model contribution is obtained by introducing a threshold , @xmath80 , at which strength in the excited states becomes significant .
such a model is motivated well by the experimental cross section for @xmath81 hadrons .
in addition , this approach accommodates the finite widths and the multi - particle continuum of qcd .
hence the excited state contributions associated with the identity operator are @xmath82 continuum model contributions associated with @xmath83 and @xmath84 terms are calculated in a similar manner . usually the continuum contributions are placed on the ope side of the sum rules .
hence , terms surviving in the limit @xmath76 have the following factors associated with them :    @xmath85 \ , , \\ m^4 \ , : & & \left [ 1 - e^{-w^2/m^2 } \
, \left ( { w^2 \over m^2 } + 1 \right ) \right ] \ , , \\ m^2 \ , : & &   \left [ 1 - e^{-w^2/m^2 } \right ] \ , .\end{aligned}\ ] ]    [ continuum ]    the contributions of this model relative to the ground state , whose properties one is really trying to determine , are not small .
they are typically 10 to 50% @xcite .
the validity of this somewhat crude model is relied upon to cleanly remove the excited state contaminations .    in the lattice qcd investigation of the continuum model formulated in euclidean space @xcite , the short time regime of point - to - point lattice correlation functions
is described well by the qcd - sum - rule - inspired continuum model . there , the laplace transform of the spectral density appears to be sufficient to render any structure in the spectral density insignificant in the short euclidean time regime of point - to - point correlators .
similar conclusions are expected to hold for the borel transformed sum rules under investigation here .    it should be noted that the spectral density of the physical world is different than that encountered on the lattice .
for example , in the nucleon channel one expects strength in the correlator above the ground state to start at the @xmath86 threshold @xcite , whereas in the lattice calculations such contributions are altered by the heavier quark mass and the quenched approximation .
while the lattice results are encouraging , one may find that in some cases the continuum model is not sufficiently detailed when the continuum model contributions are large .
an example of this will be presented in the discussion surrounding the @xmath0-meson sum rules , where discrepancies may be due to an over simplified continuum model .
the origin of the uncertainty estimate for the ope lies in the uncertainties assigned to the qcd parameters appearing in the truncated ope . in this section
we assess the uncertainties on these qcd parameters .    in assigning these uncertainties one would like to select values conservatively enough such that the qcd - sr approach can test the validity of our present understanding of qcd .
qcd - sr predictions that fail to agree with experiment at the @xmath87 level should indicate a possibly interesting discrepancy . at the same time
one would also like to select uncertainties in accord with the generally accepted values in the literature .
in addition , the error assessments need to reflect uncertainties associated with the factorization of higher dimensional operators and small but never - the - less neglected @xmath3 corrections .
these considerations have been balanced in the following . while some may argue that some values are known better than that indicated here
, others will no doubt argue that the errors are underestimated .    in any event
, one will learn for the first time how the uncertainties in the qcd parameters are mapped into uncertainties in the phenomenological fit parameters . in the lattice qcd investigation of ref .
@xcite it was noted that interplay between the pole position and continuum threshold lead to rather large uncertainties in the fit parameters .
it will be interesting to see how reasonable estimates of the qcd uncertainties are revealed in the fit parameters .    for the quark condensate , @xmath88
we note that there are two values commonly found in qcd - sr analyses .
pcac and pion decay considerations coupled with quark masses of @xmath89 mev lead to a value of @xmath90 gev@xmath91 .
qcd - sr considerations of @xmath92 and octet baryon magnetic moments @xcite prefer a smaller magnitude at @xmath93 gev@xmath91 .
expecting the value to lie some where in this regime , we take the average of these values with an uncertainty of half the difference at an 80% confidence level . introducing the standard qcd - sr notation
we define @xmath94 this value and uncertainty is also in accord with @xmath95 obtained from more recent quark mass estimates @xcite .
early estimates of the gluon condensate @xmath96 from charmonium sum rules @xcite place the value at @xmath97 gev@xmath98 , which is generally referred to as the `` standard value '' @xmath99 however , a number of more recent investigations place the gluon condensate at much larger values @xcite . for example , a finite - energy sum rule ( fesr ) analysis of the @xmath0-meson channel places the gluon condensate at 2 to 5 times the standard value @xcite . in this analysis the spectral density is taken from a fit of the isospin - one channel of @xmath100 scattering data .
first results from the aleph and cleo ii experiments @xcite also suggest a larger value at @xmath101 gev@xmath98 .
a recent laplace based sum rule analyses @xcite suggests a value of @xmath102 gev@xmath98 , more in accord with the standard value .
however , we are somewhat skeptical of the refinements of this result as the analysis is based on an approach which is criticized in the following . the most reliable sum rules of sec . [ nucleon13 ] also prefer a larger gluon condensate in order to maintain consistency between the sum rules . there
it is found that the nucleon mass may be reproduced with a value of @xmath103 gev@xmath98 which agrees with experimental estimate .
since the gluon condensate plays a critical role in the nucleon sum rules as part of a factorized dimension seven operator we adopt an uncertainty of 50% and take @xmath104    the mixed condensate is parameterized as @xmath105 . early analyses of baryons @xcite as well as more recent analyses of heavy - light quark systems @xcite place @xmath106 gev@xmath107 , while an analysis of @xmath108 qcd - sr sum rules demand @xmath109 to be smaller @xcite at @xmath110 gev@xmath107 .
hence we take the average and half the difference as the uncertainty . @xmath111 since the parameter @xmath109 is selected via monte carlo , the value of the mixed condensate is obtained by multiplying @xmath109 by the central value of the quark condensate , as opposed to a randomly selected value of the quark condensate .
this is in contrast to factorized operators where both condensates appearing in the operator are selected via monte carlo .    relatively little
is known about the magnitude of the dimension - six , four - quark operators .
early arguments placed the values of these condensates within 10% of the factorized values @xcite . however other analyses claimed significant violation of factorization @xcite for these operators in both nucleon and vector meson sum rules . parameterizing the four quark - operators as @xmath112 , estimates of factorization violation place @xmath113 or more . for the nucleon we consider @xmath114
if one adopts the standard value for the gluon condensate , then @xmath115 is sufficient to reproduce the @xmath0-meson mass
. however , selection of the more recent and reliable estimates of @xmath116 demands @xmath117 to reproduce the @xmath0-meson mass .
this value is in accord with the observations of refs .
@xcite , however the uncertainty on this parameter is large .
hence we adopt as opposed to 2 , even a small suppression of the chirally even four - quark condensate at finite density can give rise to substantial suppression of the @xmath0-meson mass and width . further investigation of these issues is warranted . ]
@xmath118 figure [ distrib ] displays histograms drawn from 1000 qcd parameter sets for four of the qcd parameters discussed to this point .
=    variation of the qcd scale parameter @xmath119 has little effect on the results .
however , we restrict the values to the conventional range adopted in qcd - srs @xmath120 variation of @xmath119 above 200 mev is only sensible if the separation scale @xmath43 of the ope , taken to be 500 mev , is increased to maintain plausible convergence of the power corrections of the ope .
a change in @xmath43 redefines the condensate values .
a systematic determination of the new condensate values is beyond the scope of this investigation and we choose to simply restrict the value as above .
it is perhaps interesting to note that if @xmath119 is very large , say near the upper limit of the current world average @xcite , the ope separation scale would necessarily become large .
consequently , the borel mass would be restricted to larger values , and access to ground state properties may not be possible .
the correlator may be dominated by strength from excited states , similar to the present situation for the pseudoscalar pion correlator .
fortunately , the lower bound on the current world average for @xmath121 bodes well for present qcd sum rule analyses .
however , if the current world average for @xmath119 persists , a comprehensive reformulation of the approach should be pursued .
the @xmath0-meson sum rule requires knowledge of @xmath122 at the scale of 1 gev@xmath107 .
since this scale is rather low for two - loop perturbation theory , we adopt a 10% uncertainty @xmath123    finally , the @xmath0-meson sum rule requires knowledge of the average light quark mass @xmath44 .
since it appears multiplied by the quark condensate , we will use the relation of ( [ gor ] ) and equate @xmath124    monte - carlo techniques might also be used to probe ope truncation errors by randomly selecting the size of the next term to follow the truncated series .
the problem lies in determining a sensible and realistic estimate without actually calculating the term .
instead , we will emphasize the importance of ope convergence and select the borel analysis regime such that truncated terms are very unlikely to be important .
gaussian distributions are generated using established algorithms for the conversion of uniform distributions on the interval ( 0,1 ) @xcite .
a particularly fast algorithm used in this investigation is as follows .
if @xmath125 and @xmath126 are uniform on ( 0,1 ) , then construct @xmath127 and @xmath128 .
calculate @xmath129 . if @xmath130 start over .
otherwise , calculate @xmath131 @xmath132 and @xmath133 are independent and normally distributed with mean 0 and variance 1 .
a set of gaussianly - distributed randomly - selected condensate values is generated , from which an ope in borel space , @xmath134 , is constructed . by repeating this procedure , an uncertainty for the ope may be determined .
selecting @xmath135 evenly distributed points in the borel parameter space @xmath136 , the standard deviation in the ope at the @xmath137th borel mass is given by @xmath138 ^ 2   \ , , \label{sigmaope}\ ] ] where @xmath18 denotes the @xmath18th set of @xmath139 qcd parameter sets and @xmath140 in practice , @xmath141 was used in estimating the ope uncertainty .
however , @xmath142 was sufficient to get a stable uncertainty estimate .
the uncertainties in the ope are not uniform throughout the borel regime .
they are larger at the lower nonperturbative end of the borel region where uncertainties in the higher - dimensional vacuum condensates dominate .
hence , it is crucial that the appropriate weight is used in the calculation of @xmath4 .
having determined the ope uncertainty , a @xmath4 measure is easily constructed . for the ope obtained from the @xmath143th set of qcd parameters ,
the @xmath4 per degree of freedom is @xmath144 ^ 2       \over      \sigma_{{\rm ope}}^2(m_j ) } \ , , \label{chisq}\ ] ] where @xmath145 is the number of phenomenological search parameters , and @xmath146 denotes the phenomenological spectral representation of the qcd - sr . in practice ,
@xmath147 points were used along the borel axis .
@xmath4 is minimized by adjusting the phenomenological parameters including the pole residue @xmath148 , the hadron mass @xmath149 , and the continuum threshold @xmath150 .
distributions for the phenomenological parameters are obtained by minimizing @xmath4 for many qcd parameter sets .
provided the resulting distributions are gaussian , an estimate of the uncertainty in the phenomenological results such as the hadron mass is obtained from @xmath151 in the event that the resultant distribution is not gaussian , we will report the median and asymmetric standard deviations from the median .
we note that ( [ sigmam ] ) is the standard deviation of the distribution and not the standard error obtained by further dividing by @xmath152 , the number of qcd parameter sets in the sample .
the standard deviation is roughly independent of the number of qcd parameter sets used in the analysis , and directly reflects the input parameter uncertainties . in the following , we generally select @xmath153 . while 100 qcd parameter configurations are sufficient for obtaining reliable uncertainty estimates , we also wish to explore correlations among the qcd parameters and the phenomenological fit parameters .
the extra sets aid in resolving more subtle correlations .
some care must be taken in the interpretation of the @xmath154 . while we are free to choose any number of points @xmath135 along the borel axis for fitting the two sides of the sum rule , it is important to recognize that the actual number of degrees of freedom is determined by the number of terms in the ope .
for the @xmath154 to have true statistical meaning , a correlated @xmath4 calculation is required .
@xmath155 \times \nonumber \\ & & \quad c_{j\,j'}^{-1 }   \left [ \pi^{\rm ope}_k(m_{j ' } ) - \pi^{\rm sr}(m_{j ' } ; \lambda_k , m_k , w_k )   \right ]   \ , , \label{corrchisq}\end{aligned}\ ] ] where @xmath156 is the inverse covariance matrix .
the covariance matrix may be estimated by @xmath157    \left [ \pi^{\rm ope}_i(m_{j ' } ) - \overline{\pi^{\rm ope}}(m_{j ' } ) \right ] \label{covar}\ ] ] unfortunately , for most of the fits considered here the covariance matrix is ill - conditioned . inverting the covariance matrix
via singular value decomposition leads to well known pathological problems @xcite such as best fit lines which do not pass through the data points etc . while one could experiment with different approaches to stepping around this problem
, we prefer to use the robust method of determining best fit parameters by minimizing the uncorrelated @xmath4 of ( [ chisq ] ) .
even if some correlation among the data is suspected , it is still acceptable to use the uncorrelated @xmath4 fit provided the errors on the parameters are estimated using a monte - carlo based procedure like that discussed here rather than from the dependence of @xmath4 on the parameters @xcite .
when simultaneously fitting multiple sum rules which may have some fit parameters in common , ( [ chisq ] ) is simply modified to @xmath158 ^ 2       \over      \sigma_{{\rm ope}}^2{}_i(m_{ij } ) } \ , , \label{chisqmulti}\ ] ] where @xmath18 runs over the considered sum rules .
( in this example , @xmath149 is a common parameter to all @xmath159 sum rules . )
the subscript @xmath18 in @xmath160 indicates that each sum rule has a unique set of @xmath135 borel masses @xmath161 spanning the region of interest in borel space .
the subscript @xmath143 in ( [ chisq ] ) , denoting a particular set of qcd parameters , has been suppressed . in the following discussion
it will become apparent that it is advantageous to weight the contributions in the sum over the various sum rules by the size of their valid borel regimes .
if @xmath162 is the size of the valid borel regime for the @xmath18th sum rule , then ( [ chisqmulti ] ) is modified to @xmath163 ^ 2       \over      \sigma_{{\rm ope}}^2{}_i(m_{ij } ) } \ , .
\label{chisqweight}\ ] ]      the optimization method utilized in this investigation is an updated version of a direction - set routine by powell originally published in ref .
this algorithm finds the minimum value of a function by iterative variation of the function parameters which need not be independent .
modifications of the original routine improve the selection of conjugate directions in the search and avoid the possible generation of linearly dependent search directions @xcite .
powell s method is among the best methods for finding minima in a multidimensional parameter space when derivatives are not readily available @xcite .
part of the reason for the robust nature of powell s routine is the cautious criterion for ultimate convergence .
once a minimum is found , the optimal vector , @xmath164 , in the parameter space is displaced by ten times the requested accuracy .
the search proceeds from this new point until a minimum is found at vector @xmath165 .
the minimum on the line joining @xmath164 and @xmath165 is found at @xmath166 .
convergence is reached if the components of the vectors @xmath167 and @xmath168 are all less than 10% of the required accuracy .
otherwise the direction @xmath167 is utilized in a new search .    the displacement of the first minimum by ten times the requested accuracy plays a crucial role in removing the sensitivity of the ultimate minimum from the initial parameter estimates . while there is sufficient flexibility in the algorithm parameters to search in the locality of the initial parameter estimates ,
the opposite criteria has been utilized in this investigation .
one of the main goals of this investigation is to eliminate the possibility of fine tuning the initial parameter estimates which otherwise can bias the results .
instead , we are interested in the predictions of the qcd sum rules .
the fundamental qcd - sr for the @xmath0 meson follows from the consideration of ( [ twopt ] ) , ( [ chirho ] ) , ( [ rhov ] ) and ( [ continuum ] ) .
we consider the correlator , @xmath169 , where the continuum model is subtracted from the phenomenological side and placed on the ope side . @xmath170 + c_1 + { c_2\over m^2 } + { c_3\over 2!\ , m^4 } + \cdots \nonumber \\ & & + { c_m \over ( m-1 ) ! \ , ( m^2)^{m-1 } } + \cdots \ , .
\label{rhosr}\end{aligned}\ ] ] the wilson coefficients are given by @xcite @xmath171    the phenomenology of the qcd - sr is described by the vector - meson pole of interest plus the continuum model accounting for the contributions of all excited states . by working in a region where the pole dominates the phenomenology
, one can minimize sensitivity to the model and have assurance that it is the spectral parameters of the ground state of interest that are being determined by matching the sum rules . in practice
, these considerations effectively set an upper limit in the borel parameter space , beyond which the model for excited states dominates the phenomenological side .    at the same time
, the truncated ope must be sufficiently convergent as to accurately describe the true ope . since the ope is an expansion in the inverse squared borel mass
, this consideration sets a lower limit in borel parameter space , beyond which higher order terms not present in the truncated ope are significant and important . monitoring ope convergence is absolutely crucial to recovering nonperturbative phenomena in the sum - rule approach , as it is the lower end of the borel region where the nonperturbative information of the ope is most significant .
as we shall see , this information must also be accurate .    in short
, one should not expect to extract information on the ground state spectral properties unless the ground state dominates the contributions on the phenomenological side and the ope is sufficiently convergent . in this investigation
, we will analyze each individual sum rule with regard to the above criteria . a sum rule with an upper limit in borel space lower than the lower limit
is considered invalid . as a measure of the relative reliability of various sum rules we consider the size of the regime in borel space where both sides of the sum rules are valid .
this quantity is used as a weight in the calculation of the @xmath154 as indicated in ( [ chisqweight ] ) .
in addition , the size of continuum contributions throughout the borel region can also serve as a measure of reliability , with small continuum model contributions being more reliable .
the fiducial borel region is chosen @xcite such that the highest - dimensional operator(s ) ( hdo ) contribute no more than @xmath172 to the qcd side while the continuum contribution is less than @xmath173 of the total phenomenological side ( i.e. , the sum of the pole and the continuum contribution ) .
the former sets a criterion for the convergence of the ope while the latter controls the continuum contribution . while the selection of @xmath174 is obvious for pole dominance , the selection of @xmath175 is a reasonably conservative criterion that has not failed in practice . and
@xmath174 criteria are automatically explored in the monte - carlo error analysis , as the condensate values and the continuum threshold change in each sample . ]
see ref .
@xcite for further specific examples .
figure [ sr1hdo ] displays the valid borel window for the @xmath0-meson sum rule of ( [ rhosr ] ) .
the hdo contributions limited to @xmath175 of the ope and the continuum model contributions limited to 50% of the phenomenology are illustrated .
the continuum contributions are the order of @xmath176 of the phenomenological contributions at the lower end of the borel regime . here
, the pole of interest truly dominates the qcd - sr .
figure [ sr1fit ] displays the corresponding fit obtained by minimizing @xmath4 in a three parameter search including @xmath177 , @xmath178 and @xmath179 .
figure [ binrhomass ] displays the distribution of @xmath0-meson masses obtained from the 1000 sum rule fits .
the distribution corresponds to @xmath180 the uncertainty is somewhat larger than the commonly assumed 10% .
figure [ binrhof ] displays the distribution of the square of @xmath177 which corresponds to @xmath181 which agrees with the experimental value of @xmath182 gev .
the continuum threshold is @xmath183 gev which compares favorably with @xmath184 gev ( the @xmath185 mass less the half width ) , given the approximate nature of the continuum model .
= 11.7truecm =    = 11.7truecm =    = 11.7truecm =    = 11.7truecm =      in addition to the fundamental sum rule of ( [ rhosr ] ) one also finds the use of sum rules obtained by differentiating ( [ rhosr ] ) with respect to @xmath186 . for the @xmath187th derivative @xmath188 \ , \label{derivsr}\end{aligned}\ ] ] while these sum rules are obviously unnecessary for the extraction of @xmath0-meson properties , the hope is that these additional sum rules will provide new useful information that might be used to further constrain the phenomenological properties .
it is possible to isolate the @xmath0-meson mass as a function of the borel mass , by taking a ratio of the first derivative sum rule with the fundamental sum rule .
this approach has some aesthetic appeal , and has widely become the method of choice for analyzing qcd - srs .
the continuum threshold is selected to make the ratio of the two sum rules as flat as possible as a function of borel mass .
finally , the @xmath0-meson mass is selected from the point at which the ratio is most flat or stable .
unfortunately , this method has some very unpleasant shortcomings .
the ratio method does not check the validity of each individual sum rule .
it is possible to have individual sum rules that are not valid while their ratio is flat as function of borel mass . in addition , the ratio method can not account for the fact that sum rules do not work equally well .
the borel region where a sum rule is valid can vary from one sum rule to another .
moreover , large uncertainties in one sum rule can spoil the predictions of another sum rule .
finally , the continuum contributions to the sum rules are not monitored in the ratio method .
if the continuum contribution is dominant in a sum rule , one should not expect to get reliable information about the lowest resonance .
the fact that sum rules do not work equally well deserves further discussion . in ref .
@xcite a detailed comparison of fundamental and derivative sum rules is carried out for the @xmath0 meson .
these arguments are not broadly known and are generally applicable to any derivative borel - improved qcd - sr .
hence the important observations are highlighted here .
if one could carry out the ope to arbitrary accuracy and use a spectral density independent of models for excited state contributions , the predictions based on ( [ rhosr ] ) and those based on the derivative sum rules of ( [ derivsr ] ) should be the same . in practical calculations , however , one has to truncate the ope and use a simple phenomenological ansatz for the spectral density .
thus it is unrealistic to expect the sum rules to work equally well .
compare the @xmath187th derivative sum rule of ( [ derivsr ] ) with the direct sum rule of ( [ rhosr ] ) .    1 .
the perturbative contribution in the derivative sum rule has an extra factor @xmath189 relative to the corresponding term in ( [ rhosr ] ) , implying that the perturbative contribution is more important in the derivative sum rules than in the fundamental sum rule , and becomes increasingly important as @xmath187 increases .
since the perturbative term mainly contributes to the continuum of the spectral density , maintaining dominance of the lowest resonance pole in the sum rule will become increasingly difficult as @xmath187 increases .
2 .   in ( [ rhosr ] ) ,
the term proportional to @xmath190 is suppressed by a factor of @xmath191 , while it is only suppressed by @xmath192 in the derivative sum rule ( @xmath193 ) .
this implies that the convergence of the ope is much slower in the derivative sum rule than in the fundamental sum rule .
consequently , the poorly known high - order power corrections are more important in the derivative sum rule than in eq .
( [ rhosr ] ) , and become more and more important as @xmath187 increases . if one would like to restrict the size of the last term of the ope to maintain some promise of ope convergence , the size of the borel region in which the sum rules are believed to be valid is restricted .
3 .   the power corrections proportional to @xmath194 do _ not _ contribute to the @xmath195 derivative sum rule but do contribute to ( [ rhosr ] ) .
if one truncates the ope , part or all of the nonperturbative information will be lost in the derivative sum rules .
it is also worth noting that the leading power corrections are the most desirable terms to have
. they do not give rise to a term in the continuum model and they are not the last term in the ope , whose relative contribution should be restricted to maintain ope convergence .    in practice ,
the predictions based on the fundamental sum rule of  ( [ rhosr ] ) are more reliable than those from the derivative sum rules , which become less and less reliable as @xmath187 increases .
this can also be demonstrated by analyzing the sum rules numerically .
figure [ sr12hdo ] displays the valid borel windows for a simultaneous fit of the @xmath0-meson sum rules of ( [ rhosr ] ) and the first derivative @xmath196 of ( [ derivsr ] ) .
the reliability of the fundamental sum rule over the derivative sum rule is obvious .
the valid borel regime for the derivative sum rule is nearly nonexistent .
all three points of the discussion surrounding derivative sum rules are displayed in this figure .
item 1 indicates the continuum model contributions will be much larger in the derivative sum rule for a given borel mass .
indeed the continuum model contribution never drops below 45% in the valid regime .
the relative contribution of the last term in the truncated ope is larger in the derivative sum rule for a given borel mass , in agreement with the expectations of items 2 and 3 .
figure [ sr12fit ] displays the corresponding fit obtained by minimizing the weighted @xmath154 of ( [ chisqweight ] ) in a three parameter search including @xmath177 , @xmath178 and @xmath179 .
this time the logarithms of the left- and right - hand sides of equations ( [ rhosr ] ) and ( [ derivsr ] ) are plotted as a function of the inverse squared borel mass .
this method has the advantage of revealing the possible breakdown of the ope signified by deviations of the dashed curves from linearity .
uncertainties are displayed only at the ends of the borel region for clarity .
= 11.7truecm =    = 11.7truecm =    while the central values of the best fit parameters @xmath197 gev , @xmath198 gev@xmath107 , and @xmath199 gev , remain essentially unchanged , the reduction of the uncertainties comes from rather unreliable information as discussed above .
these arguments are further supported by the failure of the two sides of ( [ derivsr ] ) for @xmath200 to agree as illustrated in figure [ sr12fit ] . to further probe the reasons for this failure ,
the analysis of these sum rules is repeated with a spectral representation that accounts for the finite width of the @xmath0 meson .    to account for the finite width of the @xmath0-meson , the @xmath201 function in ( [ rhov ] ) is replaced by a normalized breit - wigner form @xmath202 and the integral over the @xmath0-meson ground state @xmath203
is handled numerically @xcite .
using the experimental width of 151.5 mev , a three parameter fit of the fundamental and derivative sum rules with the breit - wigner form leads to essentially the same results reported in figure [ sr12fit ] .
the mass , coupling , and continuum threshold all agree within uncertainties , and have similar uncertainties .
hence , these spectral properties may be reliably estimated without precise information on the @xmath0 width .
however , @xmath2 mixing provides an interesting example where finite mesonic widths can play a central role in arriving at the correct physics @xcite .
the origin of the discrepancy for the derivative sum rule is more likely to lie in short - comings of continuum model . for this sum rule
the continuum model contributions never drop below 45% , which may be too large for the simple model typically used in qcd - sr calculations .
alternatively , poorly known higher order terms , not present in the truncated ope may become important in the derivative sum rule . in the derivative sum rule , the next term in the series
is suppressed by merely @xmath204 .
the importance of carefully monitoring the validity of qcd - srs is nicely displayed in the recent debate over the behavior of vector - meson masses in nuclear matter .
there an analysis involving invalid sum rules leads to the conclusion that vector - meson masses increase in finite density @xcite , in contradiction to a previous analysis @xcite also based on the ratio method .
the debate is resolved in ref .
@xcite where valid sum rules indicate that vector - meson masses decrease for increasing nuclear matter density .
having firmly established the analysis procedure in the well understood @xmath0-meson channel , we are ready to apply the same approach to the nucleon sum rules .
we will find that many results contradict the conventional wisdom .
here we focus on the borel improved qcd sum rules for the generalized interpolator of ( [ chio ] ) .
the sum rules proportional to the correlator @xmath205 where continuum model contributions are subtracted from the phenomenological side and placed on the ope side are    @xmath206 \nonumber \\ & & \quad + { 5 + 2 \beta + 5 \beta^2 \over 256 } \ , b \ , m^2
\ , l^{-4/9 }   \left [ 1 - e^{-w_1 ^ 2/m^2 } \right ]   \nonumber \\ & & \quad + { 7 - 2 \beta - 5 \beta^2 \over 24 } \ , a^2 \ , l^{4/9 }           -{13 - 2
\beta - 11 \beta^2 \over 96 } \ ,            { m_0 ^ 2 \ , a^2 \over m^2 } \ , l^{-2/27 }   \nonumber \\ & & \qquad = \widetilde \lambda_{\cal o}^2 \ , e^{-m_n^2/m^2 } \ , , \label{nucl11atp } \end{aligned}\ ] ]    and @xmath207 \nonumber \\ & & \quad -{3 ( 1 - \beta^2 ) \over 16 } \ , m_0 ^ 2 \ , a \ , m^2
\ , l^{-14/27 }   \left [ 1 - e^{-w_2 ^ 2/m^2 } \right ] \nonumber \\ & & \quad + { 19 + 10 \beta - 29 \beta^2 \over 2 ^ 7 \ , 3 ^ 2 } \ , a \ , b    \nonumber \\ & & \qquad = \widetilde \lambda_{\cal o}^2 \ ,
m_n \ , e^{-m_n^2/m^2 } \
, , \label{nucl11at1}\]][nucl11 ]    where the dirac-@xmath26 structures are indicated on the left , and @xmath208 .
the factors of @xmath209 $ ] account for the anomalous scaling of the operators and interpolators in the leading logarithmic approximation @xcite .
let us begin with the traditional or standard nucleon sum rule analysis where @xmath3 corrections are ignored .
the more common interpolator of ( [ chisr ] ) is utilized by setting @xmath14 in ( [ nucl11 ] ) .
in addition , the traditional condensate parameter set with @xmath210 is selected in order to evaluate the validity of previously published results . for a preliminary analysis of the validities of the sum rules ,
the nucleon mass is fixed to 1 gev .
figure [ sr34k2hdo ] illustrates the valid borel regimes for these two sum rules found by implementing the procedure discussed for the @xmath0-meson .
= 11.7truecm =    a comparison of figure [ sr34k2hdo ] with the analogous figure [ sr12hdo ] for the @xmath0-meson suggests the nucleon sum rules will not be as reliable . in the nucleon sum rules
, the continuum model contributes 20% or more to the phenomenological side of the sum rules .
in particular , the continuum model contribution in the sum rule of ( [ nucl11atp ] ) exceeds that for the derivative @xmath0-meson sum rule which has already been shown to be incompatible with the fundamental sum rule .
it should not be surprising to find a similar disaster for the nucleon case .
as in the @xmath0-meson case , there is an analogous ratio method for obtaining the nucleon mass that has also become popular . here
the nucleon mass is isolated as a function of the borel mass by taking the ratio of equations ( [ nucl11at1])/([nucl11atp ] ) .
the mass is usually extracted from a regime in which the ratio is flat .
again ope convergence and ground - state dominance criteria are not checked in this approach , and too often the mass is extracted from a regime in which the continuum model dominates .
figure [ sr34k2fit ] illustrates the three parameter fit of the residue @xmath148 , and the two continuum model thresholds of @xmath211 and @xmath212 . the small borel regime and
the relatively huge ope uncertainties of the sum rule of ( [ nucl11atp ] ) indicate that this sum rule plays a negligible role in determining the spectral properties of the nucleon .
this is a perfect example of a case where large uncertainties in the sum rule of ( [ nucl11atp ] ) would spoil the predictive ability of ( [ nucl11at1 ] ) if the ratio method was used .
the origin of the large uncertainties lies in the poorly known value of the four - quark condensate .
the resultant fit parameters for this sum rule are also unphysical as the continuum threshold lies at zero below the nucleon pole .
in other words , the present form of the sum rule at @xmath213 is unable to provide any information on the ground state that is consistent with the more reliable sum rule of ( [ nucl11at1 ] ) at the structure 1 .
= 11.7truecm =    unfortunately , the community has incorrectly identified ( [ nucl11atp ] ) as the more reliable sum rule .
this belief is based on the simple fact that the leading term of the ope ( the identity operator ) has a well known vacuum expectation value . however , the @xmath3 corrections for the identity operator @xcite are uncomfortably large .
the first order corrections are independent of @xmath214 and are @xmath215 where @xmath216 .
this correction is @xmath217% and suggests perturbative corrections are completely out of control for this sum rule . since the magnitude of the identity operator contribution
is closely tied to the size of continuum model contributions , it is not surprising to find the valid borel regime for this sum rule disappears when @xmath3 corrections are added .
another factor in the traditional preference of ( [ nucl11atp ] ) over ( [ nucl11at1 ] ) is that the second term of ( [ nucl11at1 ] ) vanishes for ioffe s choice of nucleon interpolator @xmath218 .
this is really nothing more than a good reason to exclude the interpolator @xmath77 of ( [ chisr ] ) as an optimal interpolator .
when one considers the monstrous leading - order @xmath3 corrections , the possibility of a significant dimension - two power correction arising from a summation of the perturbative series @xcite , and the large uncertainties associated with the value of the four - quark condensate , one can not help but reach the conclusion that the qcd side of the sum rule of ( [ nucl11atp ] ) is really unknown .      since the sum rule of ( [ nucl11atp ] ) has traditionally been the nucleon sum rule of choice , it is interesting to consider the predictions of this sum rule alone .
applying the same techniques
as used for the @xmath0-meson to ( [ nucl11atp ] ) , a three parameter fit of @xmath148 , @xmath211 and @xmath219 yields @xmath220 while setting @xmath221 slightly increases the nucleon mass to 0.47 gev , adding leading order @xmath3 corrections has a negligible effect on the fit parameters .
consideration of the more recent condensate parameter set discussed in sec .
[ uncert ] leads to similar results @xmath222 this may be the first qcd sum rule to predict spectral properties that significantly contradict experiment .
as such , this may be the most interesting of all qcd sum rules , and we shall return to this sum rule after examining the selection of the most optimal interpolator for nucleon sum rule analyses .
the predictions of ( [ nucl11atp ] ) summarized in ( [ failure ] ) are extremely different from that found in the literature .
hence it is important to understand how these new results come about . a critical examination of the nucleon ground state results reported in ref
.  @xcite by the present author sheds considerable light on the importance of the rigorous qcd - sr analysis presented here .
@xcite was the first analysis to give careful regard to the issues of ope convergence and ground - state dominance of the phenomenological side .
moreover , all phenomenological parameters were included as fit parameters to be determined by the sum rules
. however , there are three important differences between the analysis presented here and that of ref .  @xcite .
1 .   the standard logarithmic measure @xcite of sum rule agreement was used to fit the sum rules .
this approach fails to apply an appropriate weight to sum rule discrepancies , unlike the @xmath4 measure of ( [ chisq ] ) .
2 .   the algorithm used to fit the sum rules did not employ the cautious convergence criteria discussed in section [ algorithm ] .
implementation of the new criteria aids in avoiding insignificant local minima .
3 .   only the standard condensate parameter set was considered in determining the nucleon mass . here
1000 sets are considered .
the optimal parameter set found in ref .
@xcite corresponds to a local minimum which is sensitive to the measure of sum rule agreement , the initial parameter estimates and the condensate parameter set .
other authors have managed to keep the nucleon mass closer to the experimental result by extracting the result from a regime in borel space that provides the correct answer or fixing the continuum threshold to values suggested by the particle data tables . as we shall see in section [ correl ] , fixing the continuum threshold to phenomenologically motivated values largely determines the mass .
hence you _ can _ get anything you want provided you are clever enough to fix the threshold to the _ right _ value . and since the excited states are modeled there is tremendous freedom in the threshold parameter space .
of course , such input may not be in agreement with the content of the ope . in any event , such results must not be regarded as predictions of qcd sum rules .
having established the nucleon sum rule of ( [ nucl11at1 ] ) at the structure 1 to be the potentially more reliable sum rule , it is important to determine the optimal mixing of spin-1/2 interpolators for this sum rule .
it has long been established that there are two independent interpolating fields with no derivatives having the quantum numbers of spin 1/2 and isospin 1/2 .
both are expected to excite the ground state nucleon from the vacuum . as such
, there is a history of argument over the optimal nucleon interpolating field to be used in analyses @xcite and this issue is recognized to be of paramount importance @xcite .    while some advocate an interpolating field for which the leading terms of the ope are stationary with respect to the interpolating field mixing parameter @xcite
, others argue that a balance between ope convergence and qcd continuum contributions must be maintained @xcite .
ideally , one would like to simply calculate with alternate interpolating fields and confirm that the nucleon properties remain unchanged .
however , previous implementations of the qcd - sr approach have hindered any firm conclusions .
the rigorous analysis introduced here will shed considerable light on these issues .
fortunately , the lattice approach to qcd is not plagued with the same limitations and the debate over the optimal interpolating field was resolved in ref .
there is was found that , to a good approximation , @xmath223 of ( [ chin2 ] ) excites pure qcd continuum .
since @xmath223 has negligible overlap with the ground state nucleon , it is tempting to simply conclude that the optimal interpolating field is @xmath224 .
while this is certainly the case for lattice qcd investigations , it is not the case for qcd - sr analyses . by including @xmath223 components in an interpolating field
, one can reduce the continuum contributions excited by @xmath224 and allow a broader borel analysis window in which better contact with the ground state of interest is obtained .
one of the most difficult things to monitor in the qcd - sr approach is whether the ope is sufficiently convergent for a particular value of borel mass .
the lattice results of ref .
@xcite indicate the @xmath225 correlator has the fastest converging ope , as its overlap with the nucleon ground state is negligible .
similarly , the combination @xmath226 produces an ope with the slowest convergence , as this correlation function is dominated by the ground state nucleon for small borel masses .
hence , errors made in truncating the ope are dominated by errors in the @xmath226 component of the general correlator . the relative error in the ope truncation
can be reduced by adding @xmath223 components to the correlator .
however , the @xmath223 components in the ope are simply subtracted off again by the continuum model terms .
hence the relevant error is the absolute error . for @xmath227 , this error
is dominated by @xmath226 components of the correlator . as a result ,
ope truncation errors are approximately independent of @xmath214 .
this crucial point has been neglected in previous arguments regarding the optimal nucleon interpolating field .
since the errors made in truncating the ope are dominated by errors in the @xmath226 component of the general correlator the criteria of a 10% limit on the hdo contribution to the ope should be applied to sum rules when @xmath228 .
once this lower limit is established , it may be used for all @xmath227 .
we will refer to this approach for determining the valid borel regime as the `` refined '' method in the following .
the optimal nucleon interpolator must involve @xmath224 , as this interpolator is required to maintain overlap with the ground state .
the task is to determine the optimal mixing of @xmath223 .
since @xmath223 has negligible overlap with the nucleon , the ground state contribution to a sum rule is also independent of @xmath214 .
hence , the size of the continuum model contribution is the predominant factor in determining the optimal interpolator .
figure [ optcon ] illustrates the contributions of the continuum model terms in ( [ nucl11at1 ] ) for @xmath229 gev and @xmath230 gev .
the following discussion is not dependent on the precise values of these parameters .
the first point to be made is that contributions from the continuum model are largest for @xmath231 .
this selection of mixing is the worst possible choice for extracting information on the ground state nucleon .
= 11.7truecm =    figure [ optcon ] also indicates it is possible to have vanishing continuum model contributions at @xmath232 or @xmath233 .
however , we are relying on the continuum model to account for strength in the correlator that does not have its origin in the ground state . without a continuum model
, one would need to include additional poles on the right - hand - side of ( [ nucl11at1 ] ) to account for positive and negative parity excitation strength . for @xmath234 or @xmath235
the correlator is negative indicating the sum rule is saturated by a negative parity state .
thus the optimal interpolator is @xmath236 or @xmath237 . to discriminate between these two regimes
, we turn to the higher - dimension - operator contributions which do not contribute to the continuum model .
it is these terms that provide crucial information on whether the strength in the correlator lies in the ground state or the excited states .
if these terms are absent , the optimal fit of the correlator is obtained when @xmath238 and @xmath239 . in this case the continuum model and the ope become equivalent via the laplace transform of ( [ laplacetrans ] ) and the fit is perfect . hence the higher - dimension - operator contributions should be large in magnitude .
a change in sign from the leading terms of the ope will also assist in distinguishing ground state strength from excited state strength as the change in the curvature of the correlator will be more prominent .
the last term of ( [ nucl11at1 ] ) is independent of the continuum model , and its value is plotted as a function of beta in figure [ optcon ] .
the contributions are larger for @xmath236 than for @xmath240 .
in addition , the sign of the contribution is opposite that of the continuum model contributions .
hence the preferred regime is @xmath236 .
these arguments are supported by the numerical analysis described below .    in summary ,
the lattice results indicate ope convergence , and ground state pole contributions are approximately independent of @xmath214 .
considerations of the size of continuum model contributions and the sign and magnitude of operators independent of the continuum model leads to the preferred value of @xmath241 a more precise determination of @xmath214 will depend on the details of limits for continuum model and hdo contributions , condensate values , and other parameters of the sum rules .    at @xmath242 , where the leading term of the ope is stationary with respect to @xmath214 @xcite ,
the continuum contributions are maximal .
the positive value and small magnitude of the hdo contribution indicates the stability of the leading terms of the ope will not be realized as stability in the ground state mass , coupling , nor in the continuum threshold .
another point of view argues that the optimal interpolating field is one that eliminates direct instanton contributions from the correlator @xcite .
however , such arguments are based on the assumption that direct instantons necessarily contribute to the nucleon correlator in a significant manner .
section [ instantons ] will demonstrate that direct instantons are not necessary to maintain sum rule self consistency for a wide range of interpolating fields .
the @xmath214 dependence of the direct instanton contributions to ( [ nucl11at1 ] ) is proportional to @xmath243 .
minimization of instanton contributions would place @xmath244 , very near @xmath245 where the hdo term of ( [ nucl11at1 ] ) vanishes and where continuum contributions are near their peak .
isolation of ground state contributions is nearly impossible in this regime .
therefore minimization of direct instanton contributions does not lead to an optimal interpolating field .    while it is important to establish the optimal mixing of interpolating fields for qcd - sr analyses
, one should not overlook the fact that there is a range of values for @xmath214 where the sum rules are expected to work .
moreover , the ground state contribution to all these sum rules is equivalent to the 1% level @xcite . in other words ,
the right - hand side of ( [ nucl11at1 ] ) for a single pole plus continuum model is independent of @xmath214 .
after the first sum rule is written down , additional sum rules may be introduced with merely one new fit parameter ( the continuum threshold ) per sum rule . since direct instanton contributions to the sum rules
are not independent of @xmath214 @xcite , one has an excellent opportunity to see if direct instanton contributions really are necessary to maintain sum - rule consistency .
this issue is addressed in section [ instantons ] below .
having discovered that @xmath223 has negligible overlap with the nucleon ground state , it might be interesting to return to the `` ioffe formula '' this time generalized for arbitrary @xmath214 @xmath246 \ ,               { 4 \ , \left ( 2 \pi \right ) ^2 \ ,                \langle \overline q q \rangle               \over                m^2 } \ , .\ ] ] for @xmath228 , @xmath247 gev when @xmath248 gev , or @xmath249 gev when @xmath71 gev .
the explicit dependence on the choice of interpolating field and the borel mass limits the significance of the `` ioffe formula '' to a qualitative role at best .
here we focus on the more reliable sum rule of ( [ nucl11at1 ] ) at the optimal interpolator mixing of @xmath250 .
fortunately , the @xmath3 corrections for the leading order operator of the ope are much smaller than for the identity operator at the order of 10% .
the @xmath214-dependent coefficient of the quark condensate in ( [ nucl11at1 ] ) becomes @xcite @xmath251 these corrections are small relative to uncertainties in the other terms of the correlator and can be absorbed by a small shift in the condensate itself . since our qcd parameter set
is extracted from analyses with neglected @xmath3 corrections , we shall do the same here .
armed with knowledge of the optimal interpolator , and a refined understanding of ope convergence tests , we now present the qcd - sr predictions for the nucleon spectral properties from spin-1/2 interpolators .
figure [ sr4hdo ] displays the valid borel regime for this sum rule .
the hdo contributions are less than 10% of the ope contributions for this fit when @xmath228 , as indicated by the solid curve .
as described above , contributions from correlator components known to be more convergent increase the hdo operator contributions making resolution of the pole from the continuum more probable .
= 11.7truecm =    = 11.7truecm =    figure [ sr4fit ] displays the fit for this sum rule .
the near perfect linearity of the ope@xmath252continuum side of @xmath205 supports the arguments presented above . despite 40% contributions of the hdo term to this sum rule , the ope@xmath252continuum side displays the anticipated behavior of reasonable ope convergence .
the optimal fit parameters obtained from consideration of 1000 qcd parameter sets are @xmath253 the uncertainty of the nucleon mass is more than double the traditional expectation at 260 mev .
figures [ sr4binm ] , [ sr4binw ] , and [ sr4binlambda ] display histograms for the nucleon mass , threshold and residue respectively .
these histograms will aid in the discussion of the following section .
= 11.7truecm =    = 11.7truecm =    table [ comparesr ] compares the prediction for @xmath254 with other estimates of @xmath255 .
this residue corresponds to ioffe s interpolator of ( [ chisr ] ) and is expected to be independent of @xmath227 to the order of 1% @xcite .    0.0(0 )    .comparison of predictions for @xmath256 for various approaches to qcd . [ cols="<,<,^ " , ]     what this indicates is that the intimate relationship between the quark condensate and the nucleon mass suggested in the `` ioffe formula '' is invalid .
considering the assumptions made in deriving the `` ioffe formula '' one might expect the quark condensate value to be tied more tightly to the continuum model than the ground state pole , and this is borne out in table [ nucleoncorr ] . there
the quark condensate plays the dominant role in determining the continuum model threshold .
of course , all the condensate values themselves are correlated by the theory of qcd , and this correlation has not been included in this analysis .
it may very well be that the nucleon mass is positively correlated with the quark condensate when the higher dimension condensates are also correlated with the quark condensate . the important point to recognize
is that the correlation does not come about through some direct connection such as the `` ioffe formula '' .
moreover , we are extracting nucleon properties in a borel regime in which these leading order condensates just happen to be the dominant ones .
one could also extract the mass from much deeper in the non - perturbative regime , perhaps by a lattice qcd calculation .
there , it is much higher dimension operators which account for the same nucleon mass @xcite .
therefore , it is really impossible to identify a particular condensate for giving rise to ground state hadron properties .
the only true statement that may be made about the relationship between condensates and spectral properties is what we have known all along : the condensates reflecting chiral symmetry breaking play a significant role in determining hadron properties .
correlations among the phenomenological fit parameters themselves are qualitatively different from the correlations with the qcd parameters .
figures [ rhowlambda ] , [ rhomw ] and [ rhomlambda ] display the close correspondence between all three of the fit parameters , @xmath179 , @xmath177 , and @xmath178 obtained from 1000 fits of ( [ rhosr ] ) .
this close association between the @xmath0-meson spectral parameters suggests an underlying constraint on the fit parameters .
indeed , local duality necessarily leads to correlations among the spectral parameters . to examine the correlations of local duality in an analytical but qualitative fashion ,
consider the finite energy sum rules ( fesr ) @xcite obtained from @xmath257 here , @xmath258 is given by ( [ rhovb ] ) and @xmath259 is the qcd evaluation of ( [ twopt ] ) at the structure @xmath260 .
the weighting function @xmath261 for fesr is taken to be @xmath262 . these functions weight excited state contributions in favor of ground state contributions and , as such , fesr are inappropriate for a quantitative determination of ground state properties .
however , the first moment and , to a lesser extent , the second moment of the fesrs can provide some general insight into the relationships among the spectral parameters .
taking @xmath263 and @xmath264 yields the following two relationships    @xmath265    [ fesrrho ]    in a qualitative sense , the contributions from the ope on the right - hand side of ( [ fesrs ] ) are negligible . in this case
( [ fesrrho ] ) implies    @xmath266    [ fesrrel ]    these relationships , capturing the features of local duality in a qualitative sense , are independent of the precise values of the qcd parameters .
the relationships remain prominent in the more quantitative borel sum rule analysis as indicated in figures [ rhowlambda ] , [ rhomw ] and [ rhomlambda ] . of course
, the borel analysis is done independent of these qualitative constraints .
the first relation , the most reliable fesr relation , is satisfied remarkably well .
= 11.7truecm =    = 11.7truecm =    = 11.7truecm =    = 11.7truecm =    = 11.7truecm =    an even more striking correlation is illustrated in figure [ sr4corrmw ] for the nucleon sum rule of ( [ nucl11at1 ] ) with the optimal @xmath250 .
the association between the continuum threshold and the nucleon mass is nearly one - to - one despite the underlying scatter of qcd parameter values .
the explanation for this phenomena is that there is a single term in ( [ nucl11at1 ] ) which is capable of resolving the pole from the continuum .
the two leading terms of the ope are used in deriving the continuum model .
hence if the third term is vanishingly small , the best fit is most certainly pure continuum @xmath239 as the sum rule is satisfied as a perfect laplace transform .
the few points lying above the curve in figure [ sr4corrmw ] where the mass equals the threshold are cases where the third term in the sum rule has failed to resolve the pole from the continuum .
figure [ sr4corrmlambda ] displays the corresponding association for the mass and the residue of the pole .    in the @xmath0-meson case
there are in fact two terms independent of the continuum model in ( [ rhosr ] ) .
both terms act together to resolve the pole and continuum and hence give rise to some scatter in the mass  threshold association .
scatter in the residue of the pole reflects the fact that all terms affect the strength of the pole .
the strong associations among nucleon spectral parameters reflect the underlying constraints of local duality .
relationships analogous to the fesr for the @xmath0 meson may be derived for the nucleon spectral parameters @xcite .
however , the larger continuum model contributions to the nucleon correlator make the use of fesr less desirable .
the approximate nature of the continuum model may be too crude for a weighting function that fails to suppress excited state strength in the spectral density .
the resultant fesr relations are relatively crude and do not capture the essence of local duality as well as in the @xmath0-meson case .
a very different picture develops for the correlations among the nucleon fit parameters obtained from the analysis of ( [ nucl11at1 ] ) for the unfavorable value of @xmath267 .
figure [ sr4at0.8corrmw ] displays the correlation between the mass and the continuum threshold , @xmath212 .
the corresponding correlation between the mass and the residue is illustrated in figure [ sr4at0.8corrmlambda ] .
= 11.7truecm =    = 11.7truecm =    the nature of the three regimes discussed earlier in section [ secnucl11at1.8 ] is resolved here .
the near - pure continuum fits are revealed in the lower - left corner of the mass - threshold scatter plot of figure [ sr4at0.8corrmw ] .
the collision of the continuum threshold and the nucleon mass gives rise to the intermediate regime where the gap between ground and excited states has not been resolved .
the smaller residue compensating for strength from the continuum model in these fits is illustrated in figure [ sr4at0.8corrmlambda ] .
the third pole - plus - continuum regime displays the same near one - to - one correspondence between the mass and the threshold as seen in figure [ sr4corrmw ] for the optimal @xmath250 .
the correlations among the spectral parameters are interesting because they reveal the manner in which the sum rules work .
correlations among the fit parameters from the simultaneous analysis of ( [ nucl13atgg ] ) and ( [ nucl13atgqq ] ) at @xmath228 are illustrated in figures [ sr467corrmw2 ] and [ sr467corrmlambda2 ] .
these distributions are similar to those for the celebrated @xmath0-meson sum rule , and reflect the presence of two terms separating pole from continuum .
= 11.7truecm =    = 11.7truecm =
arguments favoring the need for direct instanton contributions to the qcd side of qcd - srs @xcite are based on the assumption that the average instanton size @xcite , @xmath268 fm @xmath269 fm , the separation scale of the ope .
even though @xmath270 , proponents of direct instanton contributions still argue that the effects of instantons are not sufficiently accounted for in qcd condensates alone .
the analysis of direct instantons in nucleon sum rules are based on many additional assumptions . here
we focus on the most recent investigation of ref .
@xcite . there ,
@xmath3 corrections to the identity operator are neglected . as discussed in section [ nucleon11 ]
these leading order corrections are 50% and suggest that @xmath3 corrections are completely out of control .
moreover the possibility of a significant dimension - two power correction arising from a summation of the perturbative series , uncertainties surrounding the factorization of the four - quark operator , coupled with an absence of a valid borel regime for the sum rule of ( [ nucl11atp ] ) make any conclusions dubious .
in addition , the ope is truncated at dimension six . as we saw in section [ nucleon11 ] the dimension six , seven , and eight
operators are crucial to separating the ground state pole from the continuum . without the dimension seven term ,
the best fit of ( [ nucl11at1 ] ) is pure continuum .
the fit is perfect when @xmath239 and @xmath271 .
hence , the nucleon sum rule at the structure 1 in ref .
@xcite ( equation 12 ) is unstable by design .
any further discussion of stabilizing the sum rule by adding instanton contributions is uninteresting .
finally , issues surrounding the problems of double counting were handled crudely by simply eliminating ope contributions beyond dimension five . hence one
should not be prejudiced by existing analyses suggesting nucleon sum rules require direct instanton contributions .
the issue is really not settled .    for the spin-1/2 interpolator of ( [ chin1 ] ) , direct instanton contributions to the sum rule of ( [ nucl11at1 ] )
are governed by the factor @xcite @xmath272 the direct instanton contributions were argued to be of vital importance in ref .
@xcite in the vicinity of @xmath273 .    in the lattice qcd investigation of @xcite
it was discovered that the overlap of the interpolator @xmath223 with the nucleon ground state is one hundred times smaller than that for @xmath224 . while the nucleon mass is certainly independent of the interpolating field mixing parameter @xmath214 , one also has the residue of the pole , @xmath148 , independent of @xmath214 at the 1% level .
this anticipated independence of both ground state properties provides new opportunities for evaluating the necessity of direct instanton effects .    as a first test of sum rule consistency
, we utilize the sum rule of ( [ nucl11at1 ] ) to extract the mass as a function of the interpolating field mixing parameter @xmath214 .
as @xmath274 the terms contributing to the continuum model become very small and additional poles must be added to the phenomenological side of the sum rule to account for strength in the correlator lying above the ground state mass . as a result
we restrict @xmath275 .
valid borel regimes were found over the range @xmath276 .
the coefficient @xmath277 changes by 50% over this range . if direct instanton contributions are indeed important at @xmath278 , the absence of direct instanton contributions in this investigation might be revealed by a significant @xmath214 dependence in the nucleon mass .
figure [ sr4massevol ] displays the ratio of the nucleon mass extracted at @xmath214 to the value obtained at the optimal @xmath279 .
uncertainties are obtained from two uncorrelated samples of 200 qcd parameter sets .
as such the central value at @xmath250 need not be 1 .
the large uncertainties associated with this sum rule prevent a sensitive test of the necessity of direct instantons .
it is clear that all ratios accommodate the ratio of 1 .
= 11.7truecm =    a more demanding test may be performed using the sum rules arising from the consideration of the overlap of spin-1/2 and spin-3/2 interpolators .
the smaller uncertainties allow a more interesting examination of these issues .
in addition , the continuum model terms are independent of @xmath214 and provide a wider range of @xmath214 in which a valid borel regime exists .
figure [ massevol ] addresses the dependence of the nucleon mass obtained from the sum rules of ( [ nucl13atgg ] ) and ( [ nucl13atgqq ] ) on the interpolating field mixing parameter @xmath214 . here
, ratios of the mass obtained at finite @xmath214 to @xmath228 are displayed .
sum rule consistency demands that all ratios agree with one , and this is satisfied in figure [ massevol ] without direct instanton contributions .
the uncertainty estimates are obtained from two uncorrelated samples of 200 qcd parameter sets .
we note that a valid regime does not exist for @xmath280 .
hence the use of the ioffe interpolator ( @xmath14 ) is a particularly unfavorable choice for these sum rules . the large error bar at @xmath281 reflects the small valid borel regime used in the fit .
= 11.7truecm =    the corresponding dependence of @xmath148 on @xmath214 is plagued by huge uncertainties .
small changes in the nucleon mass which appears squared in the exponential of the phenomenological side of the sum rules leads to large fluctuations in @xmath148 .
hence it is difficult to gain any information on the possible size of direct instanton contributions .
however , they are certainly unnecessary as all ratios encompass 1 .    in order to proceed
, we consider fixing the nucleon mass to the value determined at @xmath228 in an attempt to reduce the uncertainty in @xmath148 .
figure [ lambdaevolmfixed ] displays these ratios .
the corresponding plot for the continuum threshold is given in figure [ contevol ] .
= 11.7truecm =    = 11.7truecm =    for @xmath282 the continuum threshold tends to rise and the distribution of the residue remains largely independent of @xmath214 without direct instanton contributions . for @xmath283
the continuum threshold drops to account for strength that would otherwise be absorbed by the rising nucleon mass as illustrated in figure [ massevol ] . with the threshold near the pole position , strength which might otherwise rest in the pole has been suppressed as discussed in section [ phenomcorr ] .
hence the drop in the residue is an artifact of our fixing the nucleon mass in an attempt to reduce the error bar .    in summary ,
the independence of the mass and residue from @xmath227 is satisfied within uncertainties in the absence of direct instanton contributions .
as the sum rules are improved in the future and qcd input parameters are better determined , it will be important to repeat this analysis to again evaluate the necessity of direct instanton contributions or other refinements of the qcd - sr approach .    at present ,
direct instanton contributions are not required .
this is good news indeed as the amount of modeling required for instanton estimates is relatively unpalatable when compared to that for the standard qcd - srs themselves .
in fact a preferable approach would be to simply raise the ope separation scale further into the perturbative regime to ensure that direct instanton effects are more fully accounted .
however , this does not yet appear to be necessary .
ground state properties have been found to be independent of the interpolating field mixing parameter @xmath214 . on the other hand ,
the leading terms of the ope for ( [ nucl11 ] ) are very dependent on @xmath214 .
this situation has interesting ramifications for previous analyses of coordinate space point - to - point correlation functions @xcite where there has been some confusion surrounding the physics represented in point - to - point correlation functions .    in ref .
@xcite some conclusions are drawn from an incorrect interpretation of the physics represented in two - point correlators .
these authors suggest that the behavior of a ratio of the interacting to free correlator can reveal information on clustering in the scalar diquark channel .    it should be emphasized that the phenomenological description of the two - point function involves the baryon mass , the continuum threshold describing the effective onset of excited states and the residue of the pole which indicates the ability of the particular interpolator to excite the baryon from the qcd vacuum .
information on quark configurations in the nucleon wave function is absent in point - to - point correlation functions
.    the coordinate - space correlator may be obtained from ( [ rhop ] ) or ( [ rho1 ] ) via @xmath284 the important point here is that for large @xmath285 , @xmath286 is dominated by the contribution from the ground state pole . as such @xmath286 is independent of @xmath227 at the 1% level , or
alternatively @xmath286 is independent of @xmath287 at the 10% level @xcite .    the coordinate space analyses focused on a ratio of @xmath288 where @xmath289 is the free coordinate - space correlator . for the correlator at the structure
@xmath290 studied in @xcite at @xmath14 , the generalized free correlator is given by @xmath291 since the denominator is @xmath214 dependent and the numerator is independent of @xmath214 for large @xmath285 , the ratio is dependent on @xmath214 .
the rate at which the curve turns up as the ground state dominates and the free correlator continues to drop off for increasing @xmath285 is dependent on the interpolator mixing parameter @xmath214 .
the rate at which the ratio increases above 1 describes the relative weightings of ground versus excited states generated by the interpolating field .
it has nothing to do with the possibility of scalar diquark degrees of freedom being predominant in the nucleon wave function .
moreover , a comparison with a @xmath292 correlation function for example has no meaning .
a rigorous procedure for extracting quantities of phenomenological interest from qcd sum rules has been established .
we have performed a comprehensive analysis of ground state @xmath0-meson and nucleon qcd sum rules .
the inclusion of uncertainty estimates for the phenomenological parameters has allowed us to evaluate the predictive ability and self consistency of the qcd - srs . when the analysis is done rigorously ,
qcd sum rules are predictive .
some have argued that fitting sum rules appears to be more of an art than a science .
indeed , many of the findings here contradict the conventional wisdom of both practitioners and skeptics alike @xcite and suggests that past arguments have been based more on rhetoric than sound scientific arguments .      *
the nucleon sum rule of ( [ nucl11atp ] ) , traditionally taken to be the most reliable sum rule , fails to have a borel regime where both the truncated ope is reasonably convergent and the ground state dominates the phenomenology . * the nucleon sum rule of ( [ nucl11at1 ] ) , also obtained from spin-1/2 interpolating fields , resolves the nucleon mass to @xmath293 mev .
this uncertainty is more than double the value commonly assumed by most practitioners . *
the nucleon sum rules of ( [ nucl13atgg ] ) and ( [ nucl13atgqq ] ) obtained from the consideration of a spin-3/2 isospin-1/2 interpolator projected by the generalized spin-1/2 isospin-1/2 interpolator are much more reliable than the traditional sum rules .
they have greater overlap with the ground state pole relative to the continuum model , have broader valid borel regimes and provide phenomenological estimates of nucleon spectral properties which are more stable than those obtained from the conventional sum rules . *
the nucleon mass obtained from a simultaneous analysis of ( [ nucl13atgg ] ) and ( [ nucl13atgqq ] ) at the optimal @xmath228 is resolved to @xmath294 mev using qcd parameter estimates from the literature . *
the dependence of the spectral parameters on the interpolating field mixing parameter is in accord with expectations without resort to direct instanton contributions .
this is compelling evidence that instanton physics is adequately accounted for in the nonperturbative vacuum expectation values .
the latter point is truly exciting for the qcd sum rule approach .
the amount of modeling required for instanton estimates is relatively unpalatable when compared to that for the standard qcd - srs themselves .      *
the leading order @xmath3 corrections are the order of 50% . as such , perturbative corrections to this sum rule are very likely out of control .
* summation of the perturbative series may give rise to an important dimension - two contribution @xcite .
the size of such a correction is presently unknown . *
contributions from the four - quark condensate are not well known and appear to be large . unlike the @xmath0-meson sum rule , the four - quark condensate is not suppressed by a factor of @xmath3 .
* a regime in which the ope is reasonably convergent while the phenomenological description is dominated by the ground state pole is non - existent .
if one proceeds with this sum rule , the results are very sensitive to the model for the continuum , or inaccuracies in the ope lead to incorrect results as seen here .
the `` ioffe formula '' directly relating the nucleon mass to the local chiral condensate has been found to be misleading at best . while there is a strong connection between chiral symmetry breaking and hadron properties , such detailed or direct relationships are impossible to make and it is inappropriate to be more specific on the origin of the nucleon mass .
the fundamental and derivative sum rules for the @xmath0 meson have also been examined . once again , the popular sum rule analysis has been found to be undesirable .
the use of derivative sum rules for the determination of @xmath0-meson spectral properties is unfavorable . increased continuum model contributions , poorer ope convergence , and loss of vital information in the derivative sum rule results in discrepancies between the fundamental and derivative sum rules .
the associations among the phenomenological fit parameters are particularly interesting as they reveal how the sum rules resolve the spectral properties .
we have found that the separation of the ground state from the excited state continuum model is completely determined by the high dimension operators of the ope which are not used in formulating the continuum model .
the practice of fixing the continuum threshold to a preferred value has been exposed as strongly determining the outcome of the sum rule analysis .
qcd sum rule determinations of phenomenological spectral properties are necessarily dependent on the treatment of high dimension operators where a factorization assumption is usually invoked .
moreover , careful attention to ope convergence is crucial to extracting meaningful results .
these issues are somewhat unsettling , as we are working with expansions which are asymptotic .
the emphasis here has been on exploring the qcd parameter space via monte carlo .
further refinements of the parameter space might be made by adjusting the valid borel regime for each qcd parameter set such that the 10%50% criteria are satisfied .
this was not done here in order to explore other criteria and reduce the sensitivity to a reasonable yet somewhat arbitrary criteria .
in addition , the hlder inequalities might be applied as a further check of the validity and consistency of the sum rules @xcite .
we are now in an excellent position to attack phenomena which is of current interest to the nuclear physics community .
the techniques are currently being applied to an in - medium study of vector meson properties @xcite . there
the debates over the behavior of the vector - meson masses and the sum rules to be used in extracting vector meson properties in nuclear matter are resolved . using standard in - medium condensate estimates
, the analysis leaves no doubt that vector - meson masses decrease with increasing density .
these techniques are also being applied to in - medium properties of the nucleon @xcite , by generalizing the vacuum sum rules analyzed here to finite nuclear matter density .
the monte - carlo techniques are resolving the predictive ability of the in - medium sum rules .
another investigation utilizing these monte - carlo techniques studies the qcd - sr determination of the @xmath295 meson mixing parameter @xmath16 @xcite .
there , correlations between the ratio of @xmath296 and @xmath297 quark masses are particularly interesting .
the relatively large uncertainty in the residue of the pole does not bode well for qcd sum rule analysis beyond two - point functions .
for example , the residue of the pole appears as a pre - factor to any nucleon matrix element such as @xmath92 or magnetic form factors @xcite .
it will not be surprising to find uncertainties in these observables as large as the 35% relative error in the squared residue itself . as the uncertainties in the additional vacuum matrix elements relevant to the three - point functions are folded in the uncertainties are likely to exceed 50% and may in fact approach 100% .
analyses of these issues are currently underway @xcite .    when compared to uncertainty estimates previously published in the literature , our uncertainty estimates are somewhat larger .
this is not due to any shortcoming of the approach presented here . instead ,
previous authors have failed to provide a reliable method for determining realistic uncertainties .
often , uncertainty estimates have been conjectured without any supporting quantitative analysis .
this in - depth examination of qcd sum rule self consistency paints a favorable picture for further quantitative refinements of the approach .
of course such progress will be very difficult , as there are many aspects of the approach which must be refined for any particular improvement to be relevant .
the techniques introduced here , such as maintaining independence of ground state properties from the interpolating field mixing parameter , might be used to further narrow the uncertainties of the approach .
it is extremely important to refine the qcd - sr approach such that uncertainties in the qcd input parameters may be reduced while maintaining @xmath298 .
research in this direction is required before meaningful qcd - sr - based predictions for nucleon matrix elements and moments of structure functions are possible .
i am indebted to kai hencken for his careful checking of the wilson coefficients and his discovery of new corrections in the dimension seven operators .
his contribution is key to resolving the discrepancy between the sum rules presented here and modern estimates of the vacuum condensates .
thanks also to xuemin jin for numerous beneficial discussions .
dick furnstahl and xuemin jin have provided a critical examination of this analysis and i thank them for their insightful comments .
thanks also to javed iqbal for his contributions to the evaluation of the finite width integrals for the @xmath0-meson correlator .
many years ago , jimmy law provided the optimization routine used in this investigation , and i thank him for providing this useful program .
this research was supported by the natural sciences and engineering research council of canada and the u.s .
department of energy under grant de - fg06 - 88er40427 .
d.  b. leinweber , in _ lattice 93 , proceedings of the international symposium _ , edited by t.  draper , s.  gottlieb , a.  soni , and d.  toussaint , * 34 * , ( 1994 ) 407 , dallas , tx , 1993 , nucl .
b ( proc .  suppl . ) .
f.  x. lee , d.  b. leinweber , and x.  jin , new qcd sum rules for nucleon axial vector coupling constants , in _
panic 96 , proceedings of the xiv international conference on particles and nucleii _ , river edge , nj , ( 1996 ) , tjnaf , world scientific publishing co. , u. washington pp doe / er/40427 - 10-n96 .