$$E(\z_1,\y_1;\z_2,\y_2) 
 \equiv \frac{\langle \Psi | 
\left[ \frac{g}{1 - \lambda}\right]_{\z_1,\y_1}
  \left[ \frac{g}{1 - \lambda}\right]_{\y_2,\z_2}|\Psi  \rangle }
{ \langle \Psi|\Psi \rangle }. \label{ee}$$ The color factors leading to $C_F$ can be extracted from the expectation value in [Eq. (\[ee\])]{} (the ground state expectation value of the inverse of two FP operators is an identity in the adjoint representation). In the Abelian limit, $E(\z_1 \cdots \y_2) \to \delta(\y_1 -
\z_1)\delta(\y_2 - \z_2)$ and Eq. (\[qq\]) gives the dipole field distribution, $\langle \E^2 \rangle_{QED}$. One should note that [Eq. (\[qq\])]{} contains the two self energies. These self energies are necessary to produce the correct asymptotic behavior at $x >>
 R$ for charge-neutral systems, (in QED and QCD) [*i.e*]{} $\E^2$ has to fall-off at least as $1/\x^4$ at large distances from the sources.

The infrared, $|\x| \sim |\R| >> 1/\mu$ enhancement in QCD arises from the expectation value of the inverse of the FP operator. If $\langle \E^2(\x,\R)  \rangle$ is integrated over $\x$ one obtains the expectation value of the Coulomb energy of the $Q{\bar Q}$ source. The mutual interaction energy is given by, $$\begin{aligned}
V_C(\R) &= {1\over 2}\int d\x \langle \E^2(\x,\R) \rangle \nonumber \\ 
&\hspace{-5mm}=  -C_F {{\langle \Psi | 
\left[ {g\over {1 - \lambda}} 