\documentstyle[sprocl]{article}

\addtolength{\textwidth}{1in}
\addtolength{\textheight}{1in}
\hoffset -.5in
\voffset -.5in

\input epsf


\arraycolsep1.5pt


 \newcommand{\zr}[1]{\mbox{\hspace*{#1em}}}
 \newcommand{\ID}{\mbox{{\sf 1}\zr{-0.14}\rule{0.04em}{1.55ex}\zr{0.1}}}
 \newcommand{\NN}{\mbox{\zr{0.1}\rule{0.04em}{1.6ex}\zr{-0.05}{\sf N}}}
 \newcommand{\RR}{\mbox{\zr{0.1}\rule{0.04em}{1.6ex}\zr{-0.05}{\sf R}}}
 \newcommand{\CC}{\mbox{\zr{0.1}\rule{0.04em}{1.6ex}\zr{-0.30}{\sf C}}}
 \newcommand{\ZZ}{\mbox{\sf Z\zr{-0.45}Z}}

\begin{document}

{\noindent\small FAU-TP3-00/10 \hfill  UNITU-THEP-00/15 
\hfill \vspace{5mm} 

\title{What the Infrared Behavior of QCD Green Functions can tell
       us about Confinement in the Covariant Gauge\footnote{Invited talk
       by L. v. Smekal at ``Quark Confinement and the Hadron
       Spectrum IV'', Vienna, July 4-8, 2000.} 
}

\author{Lorenz von Smekal}
\address{Universit\"at Erlangen-N\"urnberg, 
	 Institut f\"ur Theoretische Physik III,\\
	 Staudtstr.~7, 91058 Erlangen, Germany\\
	 E-mail: smekal@theorie3.physik.uni-erlangen.de}

\author{Reinhard Alkofer}
\address{Universit\"at T\"ubingen, Institut f\"ur Theoretische Physik,\\
         Auf der Morgenstelle 14, 72076 T\"ubingen, Germany\\
	 E-mail: Reinhard.Alkofer@uni-tuebingen.de}


\maketitle	 
	 
%PACS: 02.30.Rz 11.10.Gh 12.38.Aw 14.70.Dj

\begin{abstract}
We review aspects of confinement in the covariant and local 
description of QCD and discuss to what extend our present knowledge of the
infrared behavior of QCD Green functions can support this description.
In particular, we emphasize: the 
positivity violations of transverse gluon and quark states, 
the Kugo-Ojima confinement criterion, and the conditions necessary to avoid 
the decomposition property for colored clusters.
We summarize how these issues relate to the infrared behavior of the 
propagators in Landau gauge QCD as extracted from solutions to truncated
Dyson-Schwinger equations and lattice simulations.  
\end{abstract}

\section{The Covariant Description of Confinement}
Covariant quantum theories of gauge fields require indefinite metric
spaces. Abandoning the positivity of the representation space
already implies some modifications to the standard
(axiomatic) framework of quantum field theory. 
Maintaining the much stronger principle of locality,
great emphasis has been put on the idea of relating confinement to the
violation of positivity in QCD.   
Just as in QED, where the Gupta-Bleuler prescription is to
enforce the Lorentz condition on physical states, a semi-definite {\em
physical subspace} can be defined as the kernel of an operator.  
The physical states then correspond to equivalence classes of states 
in this subspace differing by zero norm components.  
Besides transverse photons covariance implies the existence of 
longitudinal and scalar photons in QED. The latter two form metric partners
in the indefinite space. The Lorentz condition eliminates
half of these leaving unpaired states of zero norm which do not contribute to
observables. Since the Lorentz condition commutes with the 
$S$-Matrix, physical states scatter into physical ones exclusively. 

Confinement in QCD can be ascribed to an analogous
mechanism:
Within the framework of BRS algebra, in the simplest version for
the BRS-charge $Q_B$ and the ghost number $Q_c$ (both hermitean with respect
to an indefinite inner product) given by,
\begin{equation} 
           Q_B^2 = 0 \; , \quad \left[ iQ_c , Q_B \right] = Q_B \; ,
\end{equation}
completeness of the nilpotent BRS-charge $Q_B$ in a state space $\mathcal{V}$
of indefinite metric is assumed. This charge generates the BRS 
transformations by ghost number graded commutators $\{ \, , \}$,
{\it i.e.}, by commutators or anticommutators for fields with even
or odd ghost number, respectively.
The semi-definite {\em physical} subspace 
${\mathcal{V}}_{\mbox{\tiny phys}}  = \mbox{Ker}\, Q_B  $
is defined on the basis of this algebra by those states which are annihilated
by the BRS charge $Q_B$. 
Since $Q_B^2 =0 $, this subspace contains the space $ \mbox{Im}\, Q_B $
of so-called daughter
states which are images of others, their parent states in $\mathcal{V}$.
A physical Hilbert space is then obtained as (the completion of) the 
covariant space of equivalence classes, the BRS-cohomology of states in the
kernel modulo those in the image of $Q_B$,
\begin{equation}
     {\mathcal{H}}(Q_B,{\mathcal{V}}) = {\mbox{Ker}\, Q_B}/{\mbox{Im}\, Q_B} 
       \simeq  {\mathcal{V}}_s \; , 
\end{equation}
which is isomorphic to the space ${\mathcal{V}}_s$ of BRS singlets. 
It is easy to see that the image is furthermore contained in the orthogonal
complement of the kernel (given completeness they are identical).
It follows that states in $\mbox{Im}\, Q_B$
do not contribute to the inner product in $\mathcal{V_{\mbox{\tiny phys}}}$.   
Completeness is thereby important in the proof of positivity for physical
states,\cite{Kug79,Nak90} because it assures the absence of metric
partners of BRS-singlets, so-called ``singlet pairs''. 

With completeness, all states in $\mathcal{V}$ are either BRS
singlets in ${\mathcal{V}}_s$ or belong to so-called quartets which are 
metric-partner pairs of BRS-doublets (of parent with daughter states).
This exhausts all possibilities. The generalization of the
Gupta-Bleuler condition on physical states, $Q_B |\psi\rangle = 0$ in
$\mathcal{V}_{\mbox{\tiny phys}}$, eliminates half of these metric partners
leaving unpaired states of zero norm 
which do not contribute to any observable.
This essentially is the quartet mechanism: 
Just as in QED, one such quartet, the elementary quartet, is formed by
the massless asymptotic states of longitudinal and timelike gluons together 
with ghosts and antighosts which are thus all unobservable. 
In contrast to QED, however, one expects the quartet mechanism also 
to apply to transverse gluon and quark states, as far as they exist
asymptotically. A violation of positivity for such states then entails
that they have to be unobservable also. 
The combined evidence for
this collected below provides strong indication 
in favor of such a violation for possible transverse gluon states. 

The members of quartets are frequently said to be confined
kinematically. This is no comprehensive explanation of confinement but it is
one aspect in its present description.   
In particular, asymptotic transverse gluon and quark states
may or may not exist in the indefinite metric space $\mathcal{V}$. If either 
of them do exist and the Kugo-Ojima criterion is realized (see below), they
belong to unobservable quartets. In that case, the BRS-transformations of their
asymptotic fields entail that they form these quartets together with
ghost-gluon and/or ghost-quark bound states, respectively.\cite{Nak90}
It is furthermore  crucial for confinement, however, to have a mass gap in
transverse gluon correlations, {\it i.e.}, the massless transverse gluon
states of perturbation theory have to dissappear even though they should
belong to quartets due to color antiscreening and superconvergence in QCD for
less than 10 quark flavors.\cite{Oeh80,Nis94,Alk00} 

The quantum mechanical interpretation 
in terms of transition probabilities 
and measurements holds between physical states  and 
for expectation values in $\mathcal{V}_{\mbox{\tiny 
phys}}$ of those operators for which all contributions from zero norm
states vanish. For a (smeared local) operator $A$ to be
observable in this sense it is necessary to be BRS-closed, $\{ iQ_B , A \} \,
= 0 $, which coincides with the requirement of its local gauge invariance.
It then follows that for all states generated from the vacuum 
$|\Omega\rangle$ by any such observable one has
     $ \langle \Omega | A^\dagger A |\Omega \rangle \ge 0 $ .

If this construction is shown to apply to a QCD description of hadrons as the
genuine physical particles of $\mathcal{H}$, one concludes from physical
$S$-matrix unitarity (with respect to the indefinite inner product) that
absorptive thresholds in hadronic amplitudes can only be due to intermediate
hadronic states. The $S$-matrix commutes with the BRS-charge, {\it i.e.}, 
it is an observable in the above sense which transforms physical states into
physical ones exclusively and without measurable effects from
possible zero norm components.\cite{Nak90} One shows likewise in this
description that anomalous thresholds arise only from the
substructure of a given hadron as being composed of other hadrons. 
The argument to establish this employs standard analyticity 
properties for hadronic amplitudes and crossing to relate them to
absorptive singularities of other hadronic amplitudes which can in turn 
be due to intermediate hadronic states only.\cite{Oeh95}
 
Confinement depends on the 
realization of the unfixed global gauge symmetries in this formulation.
The identification of the 
BRS-singlets in the physical Hilbert space $\mathcal{H}$ with
color singlets is possible only if the charge of global gauge transformations
is BRS-exact {\em and} unbroken, {\it i.e.}, well-defined in the whole of the
indefinite metric space $\mathcal{V}$. The sufficent conditions for this are
provided by the Kugo-Ojima criterion: Considering the 
globally conserved current     
\begin{equation} 
    J^a_\mu = \partial_\nu F_{\mu\nu}^a  + \{ Q_B , D_{\mu}^{ab} \bar c^b \} 
    \qquad (\mbox{with} \; \partial_\mu J^a_\mu = 0 \,) \; ,
       \label{globG}
\end{equation}
one realizes that the first of its two terms corresponds to a coboundary 
with respect to the space-time exterior derivative while the second term 
is a BRS-coboundary with charges denoted by $G^a$ and $N^a$, respectively, 

\begin{equation} 
      Q^a =  \int d^3x \,  \partial_i F_{0 i}^a \,  +\,  \int d^3x \, 
             \{ Q_B , D_{0}^{ab} \bar c^b \} \, = \, G^a \, + \, N^a \; .
        \label{globC}
\end{equation}
For the first term herein there are only two options, it is either ill-defined
due to massless states in the spectrum of $\partial_\nu F_{\mu\nu}^a $, or else
it vanishes. 

In QED massless photon states contribute to the analogues of both currents
in~(\ref{globG}), and both charges on the r.h.s. in (\ref{globC}) are
separately ill-defined. One can employ an arbitrariness in the 
definition of the generator of the global gauge transformations
(\ref{globC}), however, to multiply the first term by a suitable
constant so chosen that these massless contributions cancel.
This way one obtains a well-defined and unbroken global
gauge charge which replaces the naive definition in (\ref{globC})
above.\cite{Kug95} Roughly speaking, there are two independent structures in
the globally conserved gauge currents in QED which both contain massless
photon contributions. These can be combined 
to yield one well-defined charge as the generator of global gauge
transformations leaving any other combination spontaneously broken,
such as the displacement symmetry which lead to the identification 
the photon with its massless Goldstone boson.\cite{Nak90,Fer71} 

If $\partial_\nu F_{\mu\nu}^a $ contains no massless
discrete spectrum on the other hand, {\it i.e.}, if there is no massless
particle pole in the Fourier transform of transverse gluon correlations, then
$G^a \equiv 0$.
In particular, this is the case for channels containing massive vector fields
in theories with Higgs mechanism, and it is expected to be also the case in
any color channel for QCD with confinement for which it actually represents one
of the two conditions formulated by Kugo and Ojima. 
In both these situations one first has equally, however, 
\begin{equation}
                       Q^a \, = \, N^a \, = \, \Big\{   Q_B \, , 
       \int d^3x \,   D_{0}^{ab} \bar c^b \Big\} \; ,
\end{equation}
which is BRS-exact. The second of the two conditions for confinement
is that this charge be well-defined in the whole of the indefinite metric space
$\mathcal{V}$. Together these conditions 
are sufficient to establish that all BRS-singlet physical
states in $\mathcal{H}$ are also color singlets, and that all colored states
are thus subject to the quartet mechanism. The 
second condition thereby provides the essential 
difference between Higgs mechanism and confinement. 
The operator $D_\mu^{ab}\bar c^b$ determining the charge $N^a$ will in
general contain a  {\em massless} contribution from the elementary
quartet due to the asymptotic field $\bar\gamma^a(x)$ in the  
antighost field,  $\bar c^a\, \stackrel{x_0 \to \pm\infty}{\longrightarrow}
\, \bar\gamma^a + \cdots $ (in the weak asymptotic limit), 
\begin{equation}
          D_\mu^{ab}\bar c^b \; \stackrel{x_0 \to \pm\infty}{\longrightarrow}
              \;   ( \delta^{ab} + u^{ab} )\,   \partial_\mu \bar\gamma^b(x) +
                 \cdots  \;  .
\end{equation}
Here, the dynamical parameters $ u^{ab} $ determine the contribution 
of the massless asymptotic state to the composite field $g f^{abc} A^c_\mu
\bar c^b  \, \stackrel{x_0 \to \pm\infty}{\longrightarrow}  \,
u^{ab} \partial_\mu \bar\gamma^b + \cdots $. These parameters can be obtained
in the limit $p^2\to 0$ from the Euclidean correlation functions of this
composite field, {\it e.g.},
\begin{equation}
\int d^4x \; e^{ip(x-y)} \,
\langle  \; D^{ae}_\mu c^e(x) \; gf^{bcd}A_\nu^d(y) \bar c^c (y) \; \rangle
\; =: \; \Big(\delta_{\mu \nu} -{p_\mu p_\nu \over p^2} \Big) \, u^{ab}(p^2)
\; .  \label{Corru}
\end{equation}
The theorem by Kugo and Ojima asserts that all $Q^a = N^a$ are
well-defined in the whole of  $\mathcal{V}$ (and do not suffer from
spontaneous breakdown), if and only if
\begin{eqnarray}
                 u^{ab} \equiv u^{ab}(0)  \stackrel{!}{=} - \delta^{ab} \; .
\label{KO1}
\end{eqnarray}
Then the massless states from the elementary quartet do not contribute to 
the spectrum of the current in $N^a$, and the equivalence between physical
BRS-singlet states and color singlets is established.\cite{Kug79,Nak90,Kug95}

In contrast, if $\mbox{det}(  \ID + u ) \not=0$, the global
gauge symmetry generated by the charges $Q^a$ in eq.~(\ref{globC}) is
spontaneuosly broken in each channel in which the gauge potential 
contains an asymptotic massive vector field.\cite{Kug79,Nak90} 
While this massive vector state 
results to be a BRS-singlet, the massless Goldstone boson states which 
usually occur in some components of the Higgs field, replace the 
third component of the vector field in the elementary
quartet and are thus unphysical. 
Since the broken charges
are BRS-exact, this {\em hidden} 
symmetry breaking is not directly observable in the Hilbert space of physical
states $\mathcal{H}$.  

It is nevertheless instructive to classify the different
szenarios according to the realization of the global gauge symmetry on the
whole of the indefinite metric space $\mathcal{V}$ of covariant gauge
theories. If it is unbroken, {\it i.e.}, as for QED and QCD, 
the first condition is crucial for confinement.
Namely, it is then necessary  to
have a mass gap in transverse gluon correlations ({\it i.e.}, in
$\partial_\nu F_{\mu\nu}^a $), since otherwise one could in principle have
{\em non-local} physical (BRS-singlet and thus gauge invariant) states which
are no color singlets, just as one has non-local gauge invariant
charged states in QED ({\it e.g.}, the state of one electron alone in the
world with its long-range Coulomb tail). 
Indeed, with unbroken global gauge invariance QED and QCD have in common
that any gauge invariant localized state must be
chargeless/colorless.\cite{Nak90} The
question is the extension to non-local states as approximated by local ones.
In QED this leads to the so-called charge superselection sectors,\cite{Haa96}
and non-local physical states with charge arise.    
If in QCD, with  unbroken global gauge symmetry {\em and}  mass gap, {\em
every} gauge-invariant state can be approximated by gauge-invariant localized
ones (which are colorless), one concludes that {\em every} gauge invariant
(BRS-singlet) state in $\mathcal{H}$ must be a color singlet.  

The (2nd condition in the) Kugo-Ojima confinement criterion,
$u = -\ID$ leading to well-defined charges $N^a$, can in Landau gauge be
shown by standard arguments employing Dyson-Schwinger
equations and Slavnov-Taylor identities to be 
equivalent to an infrared enhanced ghost propagator.\cite{Kug95}
In momentum space the nonperturbative ghost propagator of Landau gauge QCD  
is related to the form factor occuring in the correlations of
eq.~(\ref{Corru}) as follows,
\begin{equation}
    D_G(p) = \frac{-1}{p^2}      \, \frac{1}{ 1 + u(p^2) } \, , \;\;
                 \mbox{with}  \; \;   
                 u^{ab}(p^2)  = \delta^{ab}  u(p^2) \, . \label{DGdef}
\end{equation}
The Kugo-Ojima criterion, $u(0) = -1 $, thus entails that the Landau gauge
ghost propagator should be more singular than a massless particle pole in the
infrared. Indeed, we will present quite compelling evidence for this exact
infrared enhancement of ghosts in Landau gauge. 

The remaining dynamical aspect of confinement in this
formulation resides in the cluster decomposition property of local quantum 
field theory.\cite{Haa96} 
Including the indefinite metric spaces of covariant gauge
theories it can roughly be summarized as
follows: For the vacuum expectation values of
two (smeared local) operators $A$ and $B$, translated to a large spacelike
separaration $R$ of each other one obtains the following bounds depending on
the existence of a finite gap $M$ in the spectrum of the 
mass operator in $\mathcal{V}$,\cite{Nak90}   
\begin{eqnarray} 
        \Big|  \langle  \Omega | A(x) B(0) |\Omega \rangle  &-&         
 \langle  \Omega | A(x) |\Omega \rangle  \,  \langle  \Omega
             |  B(0) |\Omega \rangle  \Big|  \\
   && \hskip -.2cm \le  \;  \bigg\{ \begin{array}{ll} 
   \mbox{\small Const.} \, \times \, R^{-3/2 + 2N} \, e^{-MR} \!\!, \quad 
                      & \mbox{mass gap } M \; , \\
   \mbox{\small Const.} \, \times \, R^{-2 + 2N} \,, \;\; 
                      & \mbox{no mass gap} \; ,  \end{array}  \nonumber  
\end{eqnarray}
for $R^2 = - x^2 \to \infty $. Herein, positivity entails that $N = 0$, but a
positive integer $N$ is possible for the indefinite inner product structure in
$\mathcal{V}$. Therefore, in order to avoid the decomposition property
for products of unobservable operators $A$ and $B$ which 
together with the Kugo-Ojima criterion is equivalent to avoiding the
decomposition property for colored clusters, there should 
be no mass gap in the indefinite space $\mathcal{V}$. 
Of course, this implies nothing on the physical spectrum of the mass operator
in $\mathcal{H}$ which certainly should have a mass gap. 
In fact, if the cluster decomposition property holds for a product $A(x) B(0)$ 
forming a (smeared local) observable, it can be shown that both 
$A$ and $B$ are observables themselves. 
This then eliminates the possibility of scattering a physical state into
color singlet states consisting of widely separated colored clusters (the
``behind-the-moon'' problem).\cite{Nak90} 

The necessity for the absence of the massless particle pole in $\partial_\nu
F^a_{\mu\nu} $ in the Kugo-Ojima criterion shows that the (unphysical)
massless correlations to avoid the cluser decomposition property are {\em
not} the transverse gluon correlations. An infrared suppressed propagator for
the transverse gluons in Landau gauge confirms this condition. This holds
equally well for the infrared vanishing propagator obtained from
DSEs,\cite{Sme98} and conjectured in the studies of the implications of the
Gribov horizon,\cite{Gri78,Zwa92} 
as for the infrared suppressed but possibly finite ones extraced from
improved lattice actions for quite large volumes.\cite{Bon00} 

\section{The Infrared Behavior of Gluon and Ghost Propagators}
In Landau gauge the two invariant functions $Z(k^2)$ and $G(k^2)$ 
in (Euclidean) momentum space parametrize the structure of the gluon and ghost
propagators, respectively, as follows (with $G(k^2) = 1/(1+u(k^2))$, {\it
c.f.}, eq.~(\ref{DGdef})),
\begin{eqnarray} 
        D_{\mu\nu}(k) = \frac{Z(k^2)}{k^2} \, \left( \delta_{\mu\nu} -
        \frac{k_\mu k_\nu}{k^2} \right)  \; ,\quad
        D_G(k)  &=& - \frac{G(k^2)}{k^2}          \;.
\end{eqnarray}
One approach that proved quite useful 
to study the non-perturbative infrared behavior of these 
functions is provided by solutions to truncated 
Dyson-\-Schwin\-ger equations(DSEs) for the propagators, {\it i.e.},
their equations of motion.\cite{Alk00,Rob00}
The known structures in the 3-point vertex functions, most importantly from
their Slavnov-Taylor identities and exchange symmtries, are thereby employed
to establish closed systems of non-linear integral equations that are
complete on the level of the gluon, ghost and quark propagators in Landau
gauge QCD. This is possible with systematically neglecting
contributions from explicit  4-point vertices to the propagator 
DSEs as well as non-trivial 4-point scattering kernels 
in the constructions of the 3-point vertices.\cite{Sme98,Alk00}  
For the pure gauge theory this leads to the propagators DSEs 
diagrammatically represented in Fig.~\ref{GlGh} with
the 3-gluon and ghost-gluon vertices (the open circles) 
expressed in terms of the two functions $Z$ and $G$. 
Employing a one-dimensional approximation one obtains the numerical 
solutions sketched in Fig.~\ref{ZG}.\cite{Sme98,Hau98} 

\begin{figure}[t]
 \centerline{\epsfxsize=0.9\linewidth \epsfbox{GluonGhost.eps}}
\caption{Diagrammatic representation of the truncated system of 
gluon and ghost DSEs.}
\label{GlGh}
\end{figure}

Asymptotic expansions of the solutions in the infrared are known 
analytically. The leading infrared behavior is thereby uniquely determined 
by one exponent $\kappa = (61 - \sqrt{1897})/19 \approx 0.92 $,
\begin{equation} 
   Z(k^2) \, \stackrel{k^2\to 0}{\sim}   
      \,  \left(\frac{k^2}{\sigma}\right)^{2\kappa}  \quad \mbox{and} \quad
          G(k^2) \, \stackrel{k^2\to 0}{\sim}    
         \, \left(\frac{\sigma}{k^2}\right)^{\kappa} \; , \label{IRB}
\end{equation}
for which the bounds $0 < \kappa < 1$ can be established.\cite{Sme98} 
The renormalization group invariant momentum scale $\sigma $ represents the 
free parameter at this point. Its relation to the scale $\Lambda$ of
perturbative QCD is rather indirect, but a rough estimate of its value may be
obtained from the corresponding running coupling shown in Fig.~\ref{alpha}
(see below). From $\alpha_S(\mu) = 0.118$ at $\mu = M_Z = 91.2$GeV, the mass
of the $Z$-boson, one then obtains $ \sigma \cong (350{\rm MeV})^2$ which, as
a qualitative test, with $M_Z/M_\tau \cong 51.5$ yields  $\alpha_S(M_\tau) =
0.38$ at the $\tau$-mass.\cite{Alk00,Sme98} 

The infrared behavior in eqs.~(\ref{IRB})
was later confirmed qualitatively by studies of further truncated
DSEs.\cite{Atk97} Neither does it thus seem to depend on the particular
3-point vertices nor on the one-dimensional approximation employed in the
original solutions. 
While the gluon propagator is found to vanish
for small spacelike momenta in this way, an apparent contradiction with
earlier DSE studies that implied its infrared enhancement\cite{Hau96} 
can be understood
from the observation that the previously neglected ghost propagator now
assumes just this: An infrared enhancement of ghosts corresponding to $u(0)
= - 1$ which alongside with the absence of massless asymptotic 
transverse gluon states for $Z(0) = 0$ is predicted by the
Kugo-Ojima confinement criterion.\cite{Kug95}  

There are also recent lattice simulations which test this criterion
directly.\cite{Nak99} Instead of $u^{ab} = -\delta^{ab}$ they obtain
numerical values of around $u = -0.7$  for the 
unrenormalized diagonal parts and zero (within  statistical errors) for the
off-diagonal parts. After renormalization, diagonal parts very close to
$-1$ result. Taking into account the finite size effects on the
lattices employed in the simulations, these preliminary results
might be perfectly reconciled with the Kugo-Ojima confinement criterion
(\ref{KO1}).



\begin{figure*}[t]
\parbox{.49\linewidth}{\epsfxsize=\linewidth\epsfbox{ZG.eps}
\refstepcounter{figure} \label{ZG}
\centerline{\footnotesize Fig. \thefigure: DSE solutions for $Z(x)$ and
$G(x)$.\protect\cite{Sme98}}}  
\hskip -.2cm
\parbox{.49\linewidth}{\vskip -.5cm
\epsfxsize=1.12\linewidth\epsfbox{alpha.epsc}
\refstepcounter{figure} \label{alpha}
\vskip -.25cm 
\centerline{\hskip 1cm\footnotesize Fig. \thefigure: 
{\small $\alpha_S$} from the solutions in Fig.~\protect\ref{ZG}.}} 
\end{figure*}
 

Positivity violations of transverse gluon states are manifest in 
the spectral representation of (the relevant part of) 
the gluon propagator,\footnote{This expresses the fact that 
any 2-point function is analytic in the cut complex $p^2$-plane with
singularities along the time-like real axis only which holds
due to the spectrum condition in the local description also with indefinite
metric.}  
\begin{equation} 
       D(p^2) := \frac{Z(p^2)}{p^2}  = \int_0^\infty dm^2
         \, \frac{\rho(m^2)}{p^2 +  m^2} \; .
\end{equation}
From color antiscreening and unbroken global gauge symmetry
in QCD the spectral density herein asymptotically is negative and {\em
superconvergent},\cite{Oeh80,Nis94,Alk00}  
\begin{equation} 
  \rho(k^2) \stackrel{k^2\to\infty}{\sim}  - \frac{\gamma g^2}{k^2}
  \Big(g^2 \ln\frac{k^2}{\Lambda^2}\Big)^{-\gamma-1} \hskip -.3cm ,
  \quad\mbox{and}\;
  \int_0^\infty \!\! dm^2  \rho(m^2) = \left(
\frac{g_0^2}{g^2} \right)^\gamma  \to 0 \; , 
\end{equation}
since $ \gamma > 0  $ for $ N_f\le 9 $ in Landau gauge.
This imlies that it contains contributions from quartet states 
(and does therefore not need to be gauge invariant unlike in QED).
Here, we consider the one-dimensional Fourier transform 
\begin{eqnarray}
  D(t,\hbox{\bf p}^2) =  \int \frac{dp_0}{2\pi}
    \frac{Z(p_0^2 + \hbox{\bf p}^2)}{p_0^2 + \hbox{\bf p}^2} \; e^{i p_0 t}
               \, = \, 
\int_{\sqrt{\mbox{\footnotesize\bf p}^2}}^\infty  d\omega \, \rho(\omega^2\! -\!
\mbox{\bf p}^2) \,  e^{-\omega t}\; , \label{eq:gluon_FT}
\end{eqnarray}
which for $\rho \ge 0$ had to be positive definite (and one had
$\frac{d^2}{dt^2} \ln D(t,\hbox{\bf p}) \ge 0$). 
This is clearly not the case for the DSE solution shown in
Fig.~\ref{gluon_ft} which violates reflection positivity.\cite{Alk00,Sme98} 
Even though no negative $D(t, \hbox{\bf p}^2)$ have been reported in 
lattice calculations yet, the available results,\cite{Man99} {\it e.g.}, see
Fig.~\ref{Nak95_Fig2}, agree in indicating that ln $D(t, \hbox{\bf p}^2)$ is
not the convex function of the Euclidean time it should be for positive
$\rho$.\cite{Man87,Nak95} These are non-perturbative verifications of the
positivity violation for transverse gluon states which already occur in
perturbation theory. More significant for
confinement is the fact that no massless single transverse gluon 
contribution to $\rho$ exists for $Z(0) = 0$.

\begin{figure*}[t]
\parbox{.46\linewidth}{\hskip .5cm\epsfxsize=0.83\linewidth
\epsfbox{gluon_ft.eps}}
%\hskip .3cm
\hfill
\parbox{.49\linewidth}{\vskip -.7cm\hskip -.1cm
\epsfxsize=1.05\linewidth\epsfbox{Nak95_Fig2.ps}} 

\hskip .8cm\parbox{0.4\linewidth}{\refstepcounter{figure}  \label{gluon_ft}
{\footnotesize Fig. \thefigure: $D(t, \hbox{\bf p}^2)$ from DSEs for
the gluon\\[-3pt]
renormalization function $Z$ in Fig.~\ref{ZG}.}}  
\hfill
\parbox{0.4\linewidth}{\refstepcounter{figure}  \label{Nak95_Fig2}
{\footnotesize Fig. \thefigure: $D(t,\mbox{\bf p})$ at
  $(a\mbox{\bf p}) = (2\pi/48,0,0)$ for\\[-3pt] 
$\beta = 6.8$ on a $48^3   \times 64$ lattice and the free\\[-3pt] 
periodic one $\propto  \cosh(2\pi(t - N_t/2)/48)$.\cite{Nak95}}}
\vskip -.2cm
\end{figure*}


In Landau gauge, the non-renormalization 
of the gluon-ghost vertex offers a convenient
possibility to define a nonperturbative running coupling.\cite{Alk00,Sme98}  
An infrared fixed point, as obtained from the DSE solutions for 
this coupling shown in Fig.~\ref{alpha},
\begin{equation}
      \alpha_{S}(\mu) =  \frac{g^2}{4\pi\beta_0} 
      Z(\mu^2) G^2(\mu^2) \stackrel{\mu \to 0}{\longrightarrow}  
 \frac{16\pi}{3N_c} \left(\frac{1}{\kappa} - \frac{1}{2} \right)^{-1} \!\!\! 
         \simeq  9.5 \; ,
\end{equation}
which determines the 2-point interactions of
color-octet quark currents, thereby 
implied the existence of the unphysical massless
states necessary to escape the cluster decomposition of colored clusters.
The corresponding massless single particle poles should occur in colored
composite operators and by virtue of the Kugo-Ojima criterion belong to
unobservable quartets. An independent verification of this result would thus
be desirable. Some implications towards an infrared finite running coupling 
as analogously being extracted from the 3-gluon vertex might be seen in lattice
calculations.\cite{Bou98} A potential problem thereby arises in presence of
the infrared enhanced ghost propagator as entailed by the Kugo-Ojima
criterion for the Landau gauge, if asymmetric momentum subtraction
schemes are employed.\cite{Alk00,Alk99} An ideal comparison would therefore 
be possible from a lattice study of the ghost-gluon vertex in a symmetric
momentum subtraction scheme.\cite{Nak99}  

Confirmation of the important result that the gluon renormalization function
vanishes in the infrared and no massless asymptotic transverse gluon states
occur, {\it i.e.}, $Z(0) =0$, is seen in Fig.~\ref{Der98}, where 
the DSE solution of Fig.~\ref{ZG} is compared to lattice data,\cite{Lei98} 
and it was further verified recently with improved lattice actions for large
volumes.\cite{Bon00} This infrared suppression as seen in lattice
calculations thereby seems to be weaker than the DSE result, 
apparently giving rise to an infrared finite gluon
propagator $D(k) \sim Z(k^2)/k^2 $ (corresponding to an exponent 
$\kappa = 1/2 $ in~(\ref{IRB})), but a vanishing one does not seem to be 
ruled out for the infinite volume limit.\cite{Cuc97}
Similar results with finite $D(0)$ are also 
reported from the Laplacian gauge which practically 
avoids Gribov copies.\cite{Ale00} 

The infrared enhanced DSE solution for ghost propagator is compared
to lattice data in Fig.~\ref{Sum95col}. 
One observes quite compelling agreement, the numerical DSE solution fits the
lattice data at low momenta ($x \le 1$) significantly better than the fit to an
infrared singular form with integer exponents, $D_G(k^2) = c/k^2  + d/k^4$.
Though low momenta ($x<2$) were excluded in this fit, the authors concluded
that no reasonable fit of such a form was otherwise possible.\cite{Sum96} 
Therefore, apart from the question about a possible maximum at the very
lowest momenta, the lattice calculation seems to confirm the infrared
enhanced ghost propagator with a non-integer exponent $0 < \kappa < 1$. 
The same qualitative conclusion has in fact been obtained in a 
lattice calculation of the ghost propagator in $SU(2)$,\cite{Cuc97} where its
infrared dominant part was fitted best by $D_G \sim 1/(k^2)^{1+\kappa}$ for
an exponent $\kappa $ of roughly $ 0.3$ (for $\beta = 2.7$). 

To summarize, the qualitative infrared behavior in eqs.~(\ref{IRB}), an
infrared suppression of the gluon propagator together with an infrared
enhanced ghost propagator as predicted by the Kugo-Ojima criterion for the
Landau gauge, is confirmed by the presently availabe lattice
calculations. The exponents obtained therein ($0 < \kappa \le 0.5 $) both
seem to be consistently smaller than the one obtained in solving their DSEs.
Whether the lattice data for 
the infrared behavior of both propagators  can thereby also be determined 
from one unique exponent $0<\kappa <1 $, has not yet been investigated
to our knowledge. 

L.v.S. thanks the organizers for this stimulating conference and for the
hospitality extended to him during his stay at Vienna. This work was 
supported in part by the DFG (Al 279/3-3).

\begin{figure*}[t]
\parbox{.49\linewidth}{\hskip -.5cm\epsfxsize=1.07\linewidth
\epsfbox{Der98.eps}}
\hfill
\parbox{.48\linewidth}{\hskip .2cm
\epsfxsize=0.97\linewidth\epsfbox{Sum95col.eps}}

\hskip 0.6cm\parbox{0.4\linewidth}{
\refstepcounter{figure} \label{Der98} 
{\footnotesize Fig. \thefigure: The gluon renormalization function\\[-3pt] 
from the DSE solutions of ref.\protect\cite{Sme98} (solid line)\\[-3pt] 
and from the lattice data of ref.\protect\cite{Lei98}.}}  
\hskip 1.6cm
\parbox{0.44\linewidth}{\refstepcounter{figure}  \label{Sum95col} 
{\footnotesize Fig. \thefigure: The ghost propagator from
DSEs in ref.\protect\cite{Sme98}\\[-3pt] 
(solid line) compared to data and fit (dashed with\\[-3pt]  
$ca^2 =0.744 , \; da^4=0.256 $ for $x \ge 2$) from ref.\protect\cite{Sum96}.}}
\end{figure*}

\vspace{.5cm}
\hrule width 5cm
\vspace{.1cm}

\begin{thebibliography}{99}
\setlength{\itemsep}{0pt}
\bibitem{Kug79}
T.~Kugo and I.\ Ojima, Prog. Theor. Phys. Supl. {\bf 66} (1979) 1.
\bibitem{Nak90}  
N.~Nakanishi and I.~Ojima,
{\em Covariant Operator Formalism of Gauge Theories and Quantum
Gravity}, World Scientific Lecture Notes in Physics, Vol. 27, 1990; and
references therein.
\bibitem{Oeh80}
R.~Oehme and W.~Zimmermann,
Phys.~Rev. {\bf D21} (1980) 471; {\it ibid.}, 1661;\\
R.~Oehme, Phys.~Rev. {\bf D42} (1990) 4209, Phys.~Lett. {\bf B252} (1990) 641.
\bibitem{Nis94}
K.~Nishijima, Int.~J.~Mod.~Phys. {\bf A9} (1994) 3799,
{\it ibid.}, {\bf A10} (1995) 3155; see also,
Czech.~J.~Phys. {\bf 46} (1996) 1.
\bibitem{Alk00} 
R. Alkofer and L. v. Smekal, {\em The Infrared Behavior of QCD Green's
Func\-tions}, submitted to Phys. Rep., ; and references
therein.   
\bibitem{Oeh95}
R.~Oehme, Int.~J.~Mod.~Phys. {\bf A10} (1995) 1995.
\bibitem{Kug95}
T.~Kugo, at {\em Int. Symp. BRS Sym.}, Kyoto, Sep. 1995,
.
\bibitem{Fer71}
F.~Lenz et al., Ann.~Phys. {\bf 233} (1994) 17; {\it ibid.}, 51.
\bibitem{Haa96}
R.~Haag, {\em Local Quantum Physics},
Springer, 1996; and references therein.
\bibitem{Sme98} L.~v.~Smekal, A.~Hauck and R.~Alkofer, Phys.~Rev.~Lett.\ 
{\bf 79} (1997) 3591;\\ Ann.~Phys. {\bf 267} (1998) 1.
\bibitem{Gri78}
V.~N. Gribov, Nucl.~Phys. {\bf B139} (1978) 1.
\bibitem{Zwa92}
D.~Zwanziger, Nucl.~Phys. {\bf B378} (1992) 525; {\it ibid.},
{\bf B412} (1994) 657.
\bibitem{Bon00}
F.~Bonnet et al., ; see A.~G. Williams, these proceedings.
\bibitem{Rob00}
For applications of DSEs to QCD at finite temperature and density, see 
C. D. Roberts and S. M. Schmidt, Prog. Part. Nucl. Phys. {\bf 45} (2000).
\bibitem{Hau98}
A. Hauck, L. v. Smekal and R. Alkofer, Comput. Phys. Commun. {\bf 112} (1998) 166.
\bibitem{Atk97} 
D. Atkinson and J. C. R. Bloch,
Mod.~Phys.~Lett. {\bf A13} (1998) 1055;
Phys. Rev. {\bf D58} (1998) 094036.
\bibitem{Hau96}
R. Alkofer et al., in {\em Quark Confinement and the Hadron
Spectrum II}, World Scientific (1997), 258, ;  
A. Hauck, L. v. Smekal and R.~Alkofer, Comput. Phys. Commun. {\bf 112} (1998)
149; and references therein. 
\bibitem{Nak99}
H.~Nakajima and S.~Furui,
Nucl. Phys. B (Proc. Suppl.) {\bf 83-84} (2000) 521;\\
; see S.~Furui, these proceedings.
\bibitem{Man99} J.~E. Mandula,
Phys. Rep. {\bf 315} (1999) 273; and references therein.
 \bibitem{Man87}
J.~E. Mandula and M.~Ogilvie,
Phys.~Lett. {\bf B185} (1987) 127.
\bibitem{Nak95}
A.~Nakamura et al., in {\em RCNP Confinement 1995}, ;\\ 
H.~Aiso et al., Nucl.~Phys.~(Proc. Suppl.) {\bf B53} (1997) 570.
\bibitem{Alk99}
R. Alkofer, S. Ahlig and L. v. Smekal, Fizika {\bf B8} (1999), 277;
in {\em Understanding Deconfinement in QCD}, World Scientific
(1999), 196, ; R. Alkofer and\\ 
L. v. Smekal, .
\bibitem{Bou98}
B.~All\'es et al., Nucl.~Phys. {\bf B502} (1997) 325;\\
P.~Boucaud et al., JHEP  {\bf 10} (1998) 017; {\it ibid.} {\bf 12} (1998) 004. 
\bibitem{Lei98}
D. Leinweber, J. I. Skullerud, and A. G. Williams,
Phys. Rev. {\bf D58} (1998) 031501.
\bibitem{Ale00}
C. Alexandrou, Ph. de Forcrand and E. Follana, ;
; see C. Alexandrou, these proceedings.
\bibitem{Sum96}
H.~Suman and K.~Schilling, Phys.~Lett. {\bf B373} (1996) 314.
\bibitem{Cuc97}
A.~Cucchieri, Nucl.~Phys. {\bf B508} (1997) 353; in {\em Understanding
Deconfinement in QCD}, World Scientific (1999), .
\end{thebibliography}
\end{document}










%%% ====================================================================
%%%   This file is freely redistributable and placed into the
%%%   public domain by Tomas Rokicki.
%%%  @TeX-file{
%%%     author          = "Tom Rokicki",
%%%     version         = "2.7",
%%%     date            = "25 October 1996",
%%%     time            = "10:00:05 MDT",
%%%     filename        = "epsf.tex",
%%%     address         = "Tom Rokicki
%%%                        Box 2081
%%%                        Stanford, CA 94309
%%%                        USA",
%%%     telephone       = "+1 415 855 9989",
%%%     email           = "rokicki@cs.stanford.edu (Internet)",
%%%     codetable       = "ISO/ASCII",
%%%     keywords        = "PostScript, TeX",
%%%     supported       = "yes",
%%%     abstract        = "This file contains macros to support the inclusion
%%%                        of Encapsulated PostScript files in TeX documents.",
%%%     docstring       = "This file contains TeX macros to include an
%%%                        Encapsulated PostScript graphic.  It works
%%%                        by finding the bounding box comment,
%%%                        calculating the correct scale values, and
%%%                        inserting a vbox of the appropriate size at
%%%                        the current position in the TeX document.
%%%
%%%                        To use, simply say
%%%
%%%                        \input epsf % somewhere early on in your TeX file
%%%
%%%                        % then where you want to insert a vbox for a figure:
%%%                        \epsfbox{filename.ps}
%%%
%%%                        Alternatively, you can supply your own
%%%                        bounding box by
%%%
%%%                        \epsfbox[0 0 30 50]{filename.ps}
%%%
%%%                        This will not read in the file, and will
%%%                        instead use the bounding box you specify.
%%%
%%%                        The effect will be to typeset the figure as
%%%                        a TeX box, at the point of your \epsfbox
%%%                        command. By default, the graphic will have
%%%                        its `natural' width (namely the width of
%%%                        its bounding box, as described in
%%%                        filename.ps). The TeX box will have depth
%%%                        zero.
%%%
%%%                        You can enlarge or reduce the figure by
%%%                        saying
%%%
%%%                          \epsfxsize=<dimen> \epsfbox{filename.ps}
%%%                        or
%%%                          \epsfysize=<dimen> \epsfbox{filename.ps}
%%%
%%%                        instead. Then the width of the TeX box will
%%%                        be \epsfxsize and its height will be scaled
%%%                        proportionately (or the height will be
%%%                        \epsfysize and its width will be scaled
%%%                        proportionately).
%%%
%%%                        The width (and height) is restored to zero
%%%                        after each use, so \epsfxsize or \epsfysize
%%%                        must be specified before EACH use of
%%%                        \epsfbox.
%%%
%%%                        A more general facility for sizing is
%%%                        available by defining the \epsfsize macro.
%%%                        Normally you can redefine this macro to do
%%%                        almost anything.  The first parameter is
%%%                        the natural x size of the PostScript
%%%                        graphic, the second parameter is the
%%%                        natural y size of the PostScript graphic.
%%%                        It must return the xsize to use, or 0 if
%%%                        natural scaling is to be used.  Common uses
%%%                        include:
%%%
%%%                           \epsfxsize  % just leave the old value alone
%%%                           0pt         % use the natural sizes
%%%                           #1          % use the natural sizes
%%%                           \hsize      % scale to full width
%%%                           0.5#1       % scale to 50% of natural size
%%%                           \ifnum #1>\hsize\hsize\else#1\fi
%%%                                       % smaller of natural, hsize
%%%
%%%                        If you want TeX to report the size of the
%%%                        figure (as a message on your terminal when
%%%                        it processes each figure), say
%%%                        `\epsfverbosetrue'.
%%%
%%%                        If you only want to get the bounding box
%%%                        extents, without producing any output boxes
%%%                        or \special{}, then say
%%%                        \epsfgetbb{filename}.  The extents will be
%%%                        saved in the macros \epsfllx \epsflly
%%%                        \epsfurx \epsfury in PostScript units of
%%%                        big points.
%%%
%%%                        Revision history:
%%%
%%%                        ---------------------------------------------
%%%                        epsf.tex macro file:
%%%                        Originally written by Tomas Rokicki of
%%%                        Radical Eye Software, 29 Mar 1989.
%%%
%%%                        ---------------------------------------------
%%%                        Revised by Don Knuth, 3 Jan 1990.
%%%
%%%                        ---------------------------------------------
%%%                        Revised by Tomas Rokicki, 18 Jul 1990.
%%%                        Accept bounding boxes with no space after
%%%                        the colon.
%%%
%%%                        ---------------------------------------------
%%%                        Revised by Nelson H. F. Beebe
%%%                        <beebe@math.utah.edu>, 03 Dec 1991 [2.0].
%%%                        Add version number and date typeout.
%%%
%%%                        Use \immediate\write16 instead of \message
%%%                        to ensure output on new line.
%%%
%%%                        Handle nested EPS files.
%%%
%%%                        Handle %%BoundingBox: (atend) lines.
%%%
%%%                        Do not quit when blank lines are found.
%%%
%%%                        Add a few percents to remove generation of
%%%                        spurious blank space.
%%%
%%%                        Move \special output to
%%%                        \epsfspecial{filename} so that other macro
%%%                        packages can input this one, then change
%%%                        the definition of \epsfspecial to match
%%%                        another DVI driver.
%%%
%%%                        Move size computation to \epsfsetsize which
%%%                        can be called by the user; the verbose
%%%                        output of the bounding box and scaled width
%%%                        and height happens here.
%%%
%%%                        ---------------------------------------------
%%%                        Revised by Nelson H. F. Beebe
%%%                        <beebe@math.utah.edu>, 05 May 1992 [2.1].
%%%                        Wrap \leavevmode\hbox{} around \vbox{} with
%%%                        the \special so that \epsffile{} can be
%%%                        used inside \begin{center}...\end{center}
%%%
%%%                        ---------------------------------------------
%%%                        Revised by Nelson H. F. Beebe
%%%                        <beebe@math.utah.edu>, 09 Dec 1992 [2.2].
%%%                        Introduce \epsfshow{true,false} and
%%%                        \epsfframe{true,false} macros; the latter
%%%                        suppresses the insertion of the PostScript,
%%%                        and instead just creates an empty box,
%%%                        which may be handy for rapid prototyping.
%%%
%%%                        ---------------------------------------------
%%%                        Revised by Nelson H. F. Beebe
%%%                        <beebe@math.utah.edu>, 14 Dec 1992 [2.3].
%%%                        Add \epsfshowfilename{true,false}.  When
%%%                        true, and \epsfshowfalse is specified, the
%%%                        PostScript file name will be displayed
%%%                        centered in the figure box.
%%%
%%%                        ---------------------------------------------
%%%                        Revised by Nelson H. F. Beebe
%%%                        <beebe@math.utah.edu>, 20 June 1993 [2.4].
%%%                        Remove non-zero debug setting of \epsfframemargin,
%%%                        and change margin handling to preserve EPS image
%%%                        size and aspect ratio, so that the actual
%%%                        box is \epsfxsize+\epsfframemargin wide by
%%%                        \epsfysize+\epsfframemargin high.
%%%                        Reduce output of \epsfshowfilenametrue to
%%%                        just the bare file name.
%%%
%%%                        ---------------------------------------------
%%%                        Revised by Nelson H. F. Beebe
%%%                        <beebe@math.utah.edu>, 13 July 1993 [2.5].
%%%                        Add \epsfframethickness for control of
%%%                        \epsfframe frame lines.
%%%
%%%                        ---------------------------------------------
%%%                        Revised by Nelson H. F. Beebe
%%%                        <beebe@math.utah.edu>, 02 July 1996 [2.6]
%%%                        Add missing initialization \epsfatendfalse;
%%%                        the lack of this resulted in the wrong
%%%                        BoundingBox being picked up, mea culpa, sigh...
%%%                        ---------------------------------------------
%%%
%%%                        ---------------------------------------------
%%%                        Revised by Nelson H. F. Beebe
%%%                        <beebe@math.utah.edu>, 25 October 1996 [2.7]
%%%                        Update to match changes in from dvips 5-600
%%%                        distribution: new user-accessible macros:
%%%                        \epsfclipon, \epsfclipoff, \epsfdrafton,
%%%                        \epsfdraftoff, change \empty to \epsfempty.
%%%                        ---------------------------------------------
%%%  }
%%% ====================================================================

 \immediate\write16{This is `epsf.tex' v2.7 <25 October 1996>}%
%
\newread\epsffilein    % file to \read
\newif\ifepsfatend     % need to scan to LAST %%BoundingBox comment?
\newif\ifepsfbbfound   % success?
\newif\ifepsfdraft     % use draft mode?
\newif\ifepsffileok    % continue looking for the bounding box?
\newif\ifepsfframe     % frame the bounding box?
\newif\ifepsfshow      % show PostScript file, or just bounding box?
\epsfshowtrue          % default is to display PostScript file
\newif\ifepsfshowfilename % show the file name if \epsfshowfalse specified?
\newif\ifepsfverbose   % report what you're making?
\newdimen\epsfframemargin % margin between box and frame
\newdimen\epsfframethickness % thickness of frame rules
\newdimen\epsfrsize    % vertical size before scaling
\newdimen\epsftmp      % register for arithmetic manipulation
\newdimen\epsftsize    % horizontal size before scaling
\newdimen\epsfxsize    % horizontal size after scaling
\newdimen\epsfysize    % vertical size after scaling
\newdimen\pspoints     % conversion factor
%
\pspoints = 1bp        % Adobe points are `big'
\epsfxsize = 0pt       % default value, means `use natural size'
\epsfysize = 0pt       % ditto
\epsfframemargin = 0pt % default value: frame box flush around picture
\epsfframethickness = 0.4pt % TeX's default rule thickness
%
\def\epsfbox#1{\global\def\epsfllx{72}\global\def\epsflly{72}%
   \global\def\epsfurx{540}\global\def\epsfury{720}%
   \def\lbracket{[}\def\testit{#1}\ifx\testit\lbracket
   \let\next=\epsfgetlitbb\else\let\next=\epsfnormal\fi\next{#1}}%
%
% We use \epsfgetlitbb if the user specified an explicit bounding box,
% and \epsfnormal otherwise.  Because \epsfgetbb can be called
% separately to retrieve the bounding box, we move the verbose
% printing the bounding box extents and size on the terminal to
% \epsfstatus.  Therefore, when the user provided the bounding box,
% \epsfgetbb will not be called, so we must call \epsfsetsize and
% \epsfstatus ourselves.
%
\def\epsfgetlitbb#1#2 #3 #4 #5]#6{%
   \epsfgrab #2 #3 #4 #5 .\\%
   \epsfsetsize
   \epsfstatus{#6}%
   \epsfsetgraph{#6}%
}%
%
\def\epsfnormal#1{%
    \epsfgetbb{#1}%
    \epsfsetgraph{#1}%
}%
%
\def\epsfgetbb#1{%
%
%   The first thing we need to do is to open the
%   PostScript file, if possible.
%
    \openin\epsffilein=#1
    \ifeof\epsffilein
        \errmessage{Could not open file #1, ignoring it}%
    \else                       %process the file
        {%                      %start a group to contain catcode changes
            % Make all special characters, except space, to be of type
            % `other' so we process the file in almost verbatim mode
            % (TeXbook, p. 344).
            \chardef\other=12
            \def\do##1{\catcode`##1=\other}%
            \dospecials
            \catcode`\ =10
            \epsffileoktrue         %true while we are looping
            \epsfatendfalse     %[02-Jul-1996]: add forgotten initialization
            \loop               %reading lines from the EPS file
                \read\epsffilein to \epsffileline
                \ifeof\epsffilein %then no more input
                \epsffileokfalse %so set completion flag
            \else                %otherwise process one line
                \expandafter\epsfaux\epsffileline:. \\%
            \fi
            \ifepsffileok
            \repeat
            \ifepsfbbfound
            \else
                \ifepsfverbose
                    \immediate\write16{No BoundingBox comment found in %
                                    file #1; using defaults}%
                \fi
            \fi
        }%                      %end catcode changes
        \closein\epsffilein
    \fi                         %end of file processing
    \epsfsetsize                %compute size parameters
    \epsfstatus{#1}%
}%
%
% Clipping control:
\def\epsfclipon{\def\epsfclipstring{ clip}}%
\def\epsfclipoff{\def\epsfclipstring{\ifepsfdraft\space clip\fi}}%
\epsfclipoff % default for dvips is OFF
%
% The special that is emitted by \epsfsetgraph comes from this macro.
% It is defined separately to allow easy customization by other
% packages that first \input epsf.tex, then redefine \epsfspecial.
% This macro is invoked in the lower-left corner of a box of the
% width and height determined from the arguments to \epsffile, or
% from the %%BoundingBox in the EPS file itself.
%
% This version is for dvips:
\def\epsfspecial#1{%
     \epsftmp=10\epsfxsize
     \divide\epsftmp\pspoints
     \ifnum\epsfrsize=0\relax
       \special{PSfile=\ifepsfdraft psdraft.ps\else#1\fi\space
                llx=\epsfllx\space
                lly=\epsflly\space
                urx=\epsfurx\space
                ury=\epsfury\space
                rwi=\number\epsftmp
                \epsfclipstring
               }%
     \else
       \epsfrsize=10\epsfysize
       \divide\epsfrsize\pspoints
       \special{PSfile=\ifepsfdraft psdraft.ps\else#1\fi\space
                llx=\epsfllx\space
                lly=\epsflly\space
                urx=\epsfurx\space
                ury=\epsfury\space
                rwi=\number\epsftmp
                rhi=\number\epsfrsize
                \epsfclipstring
               }%
     \fi
}%
%
% \epsfframe macro adapted from the TeXbook, exercise 21.3, p. 223, 331.
% but modified to set the box width to the natural width, rather
% than the line width, and to include space for margins and rules
\def\epsfframe#1%
{%
  \leavevmode                   % so we can put this inside
                                % a centered environment
  \setbox0 = \hbox{#1}%
  \dimen0 = \wd0                                % natural width of argument
  \advance \dimen0 by 2\epsfframemargin         % plus width of 2 margins
  \advance \dimen0 by 2\epsfframethickness      % plus width of 2 rule lines
  \vbox
  {%
    \hrule height \epsfframethickness depth 0pt
    \hbox to \dimen0
    {%
      \hss
      \vrule width \epsfframethickness
      \kern \epsfframemargin
      \vbox {\kern \epsfframemargin \box0 \kern \epsfframemargin }%
      \kern \epsfframemargin
      \vrule width \epsfframethickness
      \hss
    }% end hbox
    \hrule height 0pt depth \epsfframethickness
  }% end vbox
}%
%
\def\epsfsetgraph#1%
{%
   %
   % Make the vbox and stick in a \special that the DVI driver can
   % parse.  \vfil and \hfil are used to place the \special origin at
   % the lower-left corner of the vbox.  \epsfspecial can be redefined
   % to produce alternate \special syntaxes.
   %
   \leavevmode
   \hbox{% so we can put this in \begin{center}...\end{center}
     \ifepsfframe\expandafter\epsfframe\fi
     {\vbox to\epsfysize
     {%
        \ifepsfshow
            % output \special{} at lower-left corner of figure box
            \vfil
            \hbox to \epsfxsize{\epsfspecial{#1}\hfil}%
        \else
            \vfil
            \hbox to\epsfxsize{%
               \hss
               \ifepsfshowfilename
               {%
                  \epsfframemargin=3pt % local change of margin
                  \epsfframe{{\tt #1}}%
               }%
               \fi
               \hss
            }%
            \vfil
        \fi
     }%
   }}%
   %
   % Reset \epsfxsize and \epsfysize, as documented above.
   %
   \global\epsfxsize=0pt
   \global\epsfysize=0pt
}%
%
%   Now we have to calculate the scale and offset values to use.
%   First we compute the natural sizes.
%
\def\epsfsetsize
{%
   \epsfrsize=\epsfury\pspoints
   \advance\epsfrsize by-\epsflly\pspoints
   \epsftsize=\epsfurx\pspoints
   \advance\epsftsize by-\epsfllx\pspoints
%
%   If `epsfxsize' is 0, we default to the natural size of the picture.
%   Otherwise we scale the graph to be \epsfxsize wide.
%
   \epsfxsize=\epsfsize{\epsftsize}{\epsfrsize}%
   \ifnum \epsfxsize=0
      \ifnum \epsfysize=0
        \epsfxsize=\epsftsize
        \epsfysize=\epsfrsize
        \epsfrsize=0pt
%
%   We have a sticky problem here:  TeX doesn't do floating point arithmetic!
%   Our goal is to compute y = rx/t. The following loop does this reasonably
%   fast, with an error of at most about 16 sp (about 1/4000 pt).
%
      \else
        \epsftmp=\epsftsize \divide\epsftmp\epsfrsize
        \epsfxsize=\epsfysize \multiply\epsfxsize\epsftmp
        \multiply\epsftmp\epsfrsize \advance\epsftsize-\epsftmp
        \epsftmp=\epsfysize
        \loop \advance\epsftsize\epsftsize \divide\epsftmp 2
        \ifnum \epsftmp>0
           \ifnum \epsftsize<\epsfrsize
           \else
              \advance\epsftsize-\epsfrsize \advance\epsfxsize\epsftmp
           \fi
        \repeat
        \epsfrsize=0pt
      \fi
   \else
     \ifnum \epsfysize=0
       \epsftmp=\epsfrsize \divide\epsftmp\epsftsize
       \epsfysize=\epsfxsize \multiply\epsfysize\epsftmp
       \multiply\epsftmp\epsftsize \advance\epsfrsize-\epsftmp
       \epsftmp=\epsfxsize
       \loop \advance\epsfrsize\epsfrsize \divide\epsftmp 2
       \ifnum \epsftmp>0
          \ifnum \epsfrsize<\epsftsize
          \else
             \advance\epsfrsize-\epsftsize \advance\epsfysize\epsftmp
          \fi
       \repeat
       \epsfrsize=0pt
     \else
       \epsfrsize=\epsfysize
     \fi
   \fi
}%
%
% Issue some status messages if the user requested them
%
\def\epsfstatus#1{% arg = filename
   \ifepsfverbose
     \immediate\write16{#1: BoundingBox:
                  llx = \epsfllx\space lly = \epsflly\space
                  urx = \epsfurx\space ury = \epsfury\space}%
     \immediate\write16{#1: scaled width = \the\epsfxsize\space
                  scaled height = \the\epsfysize}%
   \fi
}%
%
%   We still need to define the tricky \epsfaux macro. This requires
%   a couple of magic constants for comparison purposes.
%
{\catcode`\%=12 \global\let\epsfpercent=%\global\def\epsfbblit{%BoundingBox}}%
\global\def\epsfatend{(atend)}%
%
%   So we're ready to check for `%BoundingBox:' and to grab the
%   values if they are found.
%
%   If we find a line
%
%   %%BoundingBox: (atend)
%
%   then we ignore it, but set a flag to force parsing all of the
%   file, so the last %%BoundingBox parsed will be the one used.  This
%   is necessary, because EPS files can themselves contain other EPS
%   files with their own %%BoundingBox comments.
%
%   If we find a line
%
%   %%BoundingBox: llx lly urx ury
%
%   then we save the 4 values in \epsfllx, \epsflly, \epsfurx, \epsfury.
%   Then, if we have not previously parsed an (atend), we flag completion
%   and can stop reading the file.  Otherwise, we must keep on reading
%   to end of file so that we find the values on the LAST %%BoundingBox.
\long\def\epsfaux#1#2:#3\\%
{%
   \def\testit{#2}%             % save second character up to just before colon
   \ifx#1\epsfpercent           % then first char is percent (quick test)
       \ifx\testit\epsfbblit    % then (slow test) we have %%BoundingBox
            \epsfgrab #3 . . . \\%
            \ifx\epsfllx\epsfatend % then ignore %%BoundingBox: (atend)
                \global\epsfatendtrue
            \else               % else found %%BoundingBox: llx lly urx ury
                \ifepsfatend    % then keep parsing ALL %%BoundingBox lines
                \else           % else stop after first one parsed
                    \epsffileokfalse
                \fi
                \global\epsfbbfoundtrue
            \fi
       \fi
   \fi
}%
%
%   Here we grab the values and stuff them in the appropriate definitions.
%
\def\epsfempty{}%
\def\epsfgrab #1 #2 #3 #4 #5\\{%
   \global\def\epsfllx{#1}\ifx\epsfllx\epsfempty
      \epsfgrab #2 #3 #4 #5 .\\\else
   \global\def\epsflly{#2}%
   \global\def\epsfurx{#3}\global\def\epsfury{#4}\fi
}%
%
%   We default the epsfsize macro.
%
\def\epsfsize#1#2{\epsfxsize}%
%
%   Finally, another definition for compatibility with older macros.
%
\let\epsffile=\epsfbox
\endinput

