%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% 2/5/2002
%% Proceedings Trim Size 8.5"x6" has changed into Trim Size 9"x6". 
%% This Latex2e file rewritten from various sources for use in the
%% preparation of the (smaller [9"x6"]) single-column proceedings 
%% Volume, latest version by R. Sankaran with acknowledgements to 
%% Susan Hezlet and Lukas Nellen.			        
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\documentclass{ws-procs9x6}
%\usepackage{showkeys}
\usepackage{graphicx}
\usepackage{latexsym}
\usepackage{amssymb}

\newcommand{\dsp}{\displaystyle}
\newcommand{\ba}{\begin{eqnarray}}
\newcommand{\ea}{\end{eqnarray}}
\newcommand{\be}{\begin{equation}}
\newcommand{\ee}{\end{equation}}
\newcommand{\re}{\mbox{Re}}
\newcommand{\im}{\mbox{Im}}
\newcommand{\kob}{{\overline{K^0}}}
\newcommand{\epspeps}{\ensuremath{\varepsilon^\prime/\varepsilon}}
\newcommand{\epspepsfrac}{\ensuremath{\frac{\varepsilon^\prime}{\varepsilon}}}
\newcommand{\hepph}[1]{{hep-ph/#1}}



\newcommand{\citebk}[1]{[\hspace{0.mm}\raisebox{-1.85mm}[0mm][0mm]
  {\Large\cite{#1}}\hspace{-0.1mm}]}

\newcommand{\citebkcap}[1]{[\hspace{0.mm}\raisebox{-1.5mm}[0mm][0mm]
  {\large\cite{#1}}\hspace{-0.2mm}]}

\begin{document}
%\begin{titlepage}
\thispagestyle{empty}
\begin{flushright}
LU TP 02-27\\
\\
July 2002
\end{flushright}
\vfill
\begin{center}
{\huge{\bf  Penguins 2002: Penguins in $K\to\pi\pi$ Decays}$^a$}
\end{center}

\vfill

\begin{center}
{\large\bf Johan~Bijnens}\\[1cm]
{\large Department of Theoretical Physics, Lund University,\\[0.5cm]
S\"olvegatan 14A, S22362 Lund, Sweden}
\end{center}

\vfill
\begin{center}
{\bf Abstract}
\end{center}
{This talk contains a short overview of the history of the interplay of
the weak and the strong interaction and $CP$-violation. It describes
the phenomenology and the basic physics mechanisms involved in the
Standard Model calculations of $K\to\pi\pi$ decays with an emphasis
on the evaluation of Penguin operator matrix-elements.}

\vfill

\noindent $^a$
{\small
Presented at the
Symposium and Workshop "Continuous Advances in QCD 2002/Arkadyfest,"
Minneapolis, USA, 17-23 May 2002, to be published in the proceedings.}

\setcounter{page}{0}
\newpage
%\end{titlepage}


\title{Penguins 2002: Penguins in $K\to\pi\pi$ Decays}

\author{Johan Bijnens}

\address{Department of Theoretical Physics 2, Lund University,\\
S\"olvegatan 14A, S22362 Lund, Sweden\\
E-mail: bijnens@thep.lu.se}


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% You may repeat \author \address as often as necessary      %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\maketitle

\abstracts{
This talk contains a short overview of the history of the interplay of
the weak and the strong interaction and $CP$-violation. It describes
the phenomenology and the basic physics mechanisms involved in the
Standard Model calculations of $K\to\pi\pi$ decays with an emphasis
on the evaluation of Penguin operator matrix-elements.}

\section{Introduction}

In this conference in honour of Arkady Vainshtein's 60th birthday, a
discussion of the present state of the art in analytical calculations
of relevance for Kaon decays is very appropriate given Arkady's large
contributions to the field. He has summarised his own contributions
on the occasion of accepting the Sakurai Prize.\cite{Penguin3}
This also contains the story of how Penguins, at least the diagram
variety, got their name. In Fig.~\ref{fig1} I show what a real
(Linux) Penguin looks like and the diagram in a ``Penguinized'' version.
\begin{figure}
\begin{center}
\includegraphics[width=0.455\textwidth]{penguin.eps}
\nolinebreak
\includegraphics[width=0.379\textwidth]{penguindiagram.eps}
\end{center}
\caption{\label{fig1} The comparison between a (Linux) Penguin and
the Penguin diagram. (Linux Penguin from Neal Tucker,
({\tt http://www.isc.tamu.edu/$\sim$lewing/linux/}).)}
\end{figure}
This talk could easily have had other titles, examples are
``QCD and Weak Interactions of Light Quarks'' or
``Penguins and Other Graphs.'' In fact I have left out many
manifestations of Penguin diagrams. In particular I do not cover the
importance in $B$ decays where Penguins were first experimentally verified
via $B\to K^{(*)}\gamma$, but are at present more considered a nuisance
and often referred to as ``Penguin Pollution.'' Penguins also play
a major role in other Kaon decays, reviews of rare decays where
pointers to the literature can be found are
Refs~\citebk{Isidori,Littenberg,Buras2}.

In Sect.~\ref{history} I give a very short historical overview.
Sects~\ref{deltaIhalf} and \ref{epspeps} discuss the main physics
issues and present the relevant phenomenology of the $\Delta I=1/2$ rule
and the $CP$-violating quantities $\varepsilon$ and $\epspeps$.
The underlying Standard Model diagrams responsible for $CP$-violation are 
shown there as well. The more challenging part is to actually evaluate
these diagrams in the presence of the strong interaction. We can
distinguish several regimes of momenta which have to be treated using
different methods. An overview is given in Sect.~\ref{EFT}
where also the short-distance part is discussed. The more difficult
long-distance part has a long history and some approaches are mentioned,
but only my favourite method, the $X$-boson or fictitious Higgs exchange,
is described in more detail in Sect.~\ref{Xboson} where I also present results
for the main quantities. For two particular matrix-elements, 
those of the electroweak
Penguins $Q_7$ and $Q_8$, a dispersive analysis allows to evaluate these
in the chiral limit from experimental data. This is discussed in
Sect.~\ref{dispersive}. We summarise our results in the conclusions and compare
with the original hopes from Arkady and his collaborators.
This talk is to a large extent a shorter version of the
review~\citebk{QCDweak}.


\section{A short historical overview}
\label{history}

The weak interaction was discovered in 1896 by Becquerel when he discovered
spontaneous radioactivity. The next step towards a more fundamental
study of the weak interaction was taken in the 1930s when the neutron was
discovered and its $\beta$-decay studied in detail. 
The fact that the proton and electron energies did not add up
to the total energy corresponding to the mass of the neutron,
made Pauli suggest the neutrino as a solution.
Fermi then incorporated it in the first full fledged theory of the
weak interaction, the famous Fermi four-fermion~\cite{Fermi} interaction. 

\be
\label{Fermi1}
{\mathcal L}_{\mbox{Fermi}} =
\frac{G_F}{\sqrt{2}}\; [\overline p \gamma_\mu\left(1-\gamma_5\right) n]\;
[\overline e \gamma^\mu\left(1-\gamma_5\right)\nu]\,.
\ee

The first fully nonhadronic weak interaction came after world-war two
with the muon discovery and the study of its $\beta$-decay.
The analogous Lagrangian to Eq.~(\ref{Fermi1}) was soon written down.
At that point T.D.~Lee and C.N.~Yang~\cite{LeeYang} realized that there was no
evidence that parity was conserved in the weak interaction. This quickly
led to a search for parity violation both in nuclear decays~\cite{CSWu}
and in the decay chain
$\pi^+\to\mu^+\nu_\mu\to e^+\nu_e\bar\nu_\mu\nu_\mu$.\cite{FriedmanTelegdi}
Parity violation
was duly observed in both cases. These experiments and others
led to the final form of the
Fermi Lagrangian given in
Eq.~(\ref{Fermi1}).\cite{SudarshanMarshak,FeynmanGell-Mann}

During the 1950s steadily more particles were discovered
providing many puzzles.
These were solved by the introduction of
strangeness,\cite{Pais,Gell-Mann}
of what is now known as the $K_L$ and the $K_S$~\cite{Gell-MannPais} and
the ``eightfold way'' of classifying the hadrons
into symmetry-multiplets.\cite{eightfold}

Subsequently Cabibbo realized that the weak interactions of
the strange particles were very similar to those of the
nonstrange particles.\cite{Cabibbo} He proposed that the weak interactions
of hadrons occurred through a current which was a mixture of the
strange and non-strange currents with a mixing angle now universally
known as the Cabibbo angle. The hadron symmetry group
led to the introduction of quarks~\cite{quarks} as a means of
organising which $SU(3)_V$ multiplets were present in the
spectrum.

In the same time period the Kaons provided another surprise.
Measurements at
Brookhaven~\cite{CCFT} indicated that the long-lived state, the $K_L$,
did occasionally decay to two pions in the final state as well,
showing that $CP$ was violated.
Since the $CP$-violation was small, explanations could be sought at
many scales, an early phenomenological analysis can be found
in Ref. ~\citebk{WuYang}, but as 
the socalled superweak model~\cite{Wolfenstein} showed,
the scale of the interaction involved in $CP$-violation could be much higher.

The standard model for the weak and electromagnetic interactions of leptons
was introduced in the same period.
The Fermi theory is nonrenormalisable.
Alternatives based on Yang-Mills~\cite{Yang-Mills}
theories had been proposed by Glashow~\cite{Glashow} but struggled with
the problem of having massless gauge bosons. This was solved by 
the introduction of the
Higgs mechanism by Weinberg and Salam. The model
could be extended to include the weak interactions of hadrons by adding
quarks in doublets, similar to the way the leptons were included.
One problem this produced was that loop-diagrams provided a much too
high probability for the decay $K_L\to\mu^+\mu^-$ compared to
the experimental limits. These socalled flavour changing neutral currents
(FCNC) needed to be suppressed. The solution was found in
the Glashow-Iliopoulos-Maiani mechanism.\cite{GIM} A fourth quark, the charm
quark, was introduced beyond the up, down and strange quarks.
If all the
quark masses were equal, the dangerous loop contributions to FCNC processes
cancel, the socalled GIM mechanism.
This allowed a prediction of the charm quark
mass,~\cite{GaillardLee1} soon confirmed with
the discovery of the $J/\psi$.

In the mean time, 
QCD was formulated.\cite{QCD} The property of
asymptotic freedom~\cite{Asympt} was established which explained why
quarks at short distances could behave as free particles and at the
same time at large distances be confined inside hadrons.

The study of Kaon decays still went on, and an already old problem,
the $\Delta I=1/2$ rule saw the first signs of a solution. It was
shown~\cite{GaillardLee2,AltarelliMaiani} that the short-distance QCD part
of the nonleptonic weak decays provided already an enhancement of
the $\Delta I=1/2$ weak $\Delta S=1$ transition over the $\Delta I=3/2$
one. The ITEP group extended first the Gaillard-Lee analysis
for the charm mass,\cite{Vainshtein1} but then realized that in
addition to the effects that were included in
Refs~\citebk{GaillardLee2,AltarelliMaiani},
there was a new class of diagrams that only contributed to
the $\Delta I=1/2$ transition.\cite{Penguin1,Penguin2} While, as we
will discuss in more detail later, the general class of these
contributions, the socalled Penguin-diagrams, is the most likely main
cause of the $\Delta I=1/2$ rule, the short-distance part of them
provide only a small enhancement contrary to the original hope.
A description of the early history of Penguin diagrams, including the
origin of the name, can be found in the 1999 Sakurai Prize lecture  of
Vainshtein.\cite{Penguin3}

Penguin diagrams at short distances provide nevertheless a large amount
of physics. The origin of $CP$-violation was (and partly is)
still a mystery. The superweak model explained it, but introduced new
physics that had no other predictions. 
Kobayashi and Maskawa~\cite{KobayashiMaskawa} realized that
the framework established by Ref.~\citebk{GIM} could be extended
to three generations.
The really new aspect this brings in is that $CP$-violation
could easily be produced at the weak scale and not at the much
higher superweak scale. In this Cabibbo-Kobayashi-Maskawa (CKM)
scenario, $CP$-violation comes from
the mixed quark-Higgs sector, the Yukawa sector,
and is linked with the masses and mixings of the quarks.
Other mechanisms at the weak scale
also exist, as e.g. an extended Higgs sector.\cite{WeinbergCP}

The inclusion of the CKM mechanism into the
calculations for weak decays was done by Gilman and
Wise~\cite{GilmanWise1,GilmanWise2} which provided the prediction
that $\epspeps$ should be nonzero and of the order of $10^{-3}$.
Guberina and Peccei~\cite{GuberinaPeccei} confirmed this.
This prediction spurred on the experimentalists and after
two generations of major experiments, NA48 at CERN and KTeV at Fermilab
have now determined this quantity and the qualitative prediction
that $CP$-violation at the weak scale exists is now confirmed.
Much stricter tests of this picture will happen at other Kaon
experiments as well as in $B$ meson studies.

The $K^0$-$\kob$ mixing has QCD corrections and $CP$-violating contributions
as well. The calculations of these required a proper treatment
of box diagrams and inclusions of the effects of the $\Delta S=1$
interaction squared. This was accomplished at one-loop
by Gilman and Wise a few years later.\cite{GilmanWise3,GilmanWise4}

That Penguins had more surprises in store was shown some years later
when it was realized that the enhancement originally expected on chiral grounds
for the Penguin diagrams~\cite{Penguin1,Penguin2} was present,
not for the Penguin diagrams with gluonic intermediate states, but
for those with a photon.\cite{BijnensWise} This contribution
was also enhanced in its effects by the $\Delta I=1/2$ rule. 
This lowered the
expectation for $\epspeps$, but it became significant
after it was found that the top quark had a very large mass.
Flynn and Randall~\cite{RandallFlynn} reanalysed the electromagnetic
Penguin with a large top quark mass and included also $Z^0$ exchange.
The final effect was that
the now rebaptized electroweak Penguins could have a very large contribution
that could even cancel the contribution to $\epspeps$ from
gluonic Penguins. This story still continues at present and the
cancellation, though not complete, is one of the major impediments
to accurate theoretical predictions of $\epspeps$.

The first calculation of two-loop effects in the short-distance part
was done in Rome~\cite{twoloopold} in 1981.
The value of $\Lambda_{\mbox{QCD}}$ has risen
from values of about 100~MeV to more than 300~MeV. A full calculation
of all operators at two loops thus became necessary, taking into account
all complexities of higher order QCD. This program was finally
accomplished by two independent groups. One in Munich around A.~Buras
and one in Rome around G.~Martinelli. 

\section{$K\to\pi\pi$ and the $\Delta I=1/2$ rule}
\label{deltaIhalf}

The underlying qualitative difference we want to understand is the
$\Delta I=1/2$ rule.
We can try to calculate $K\to\pi\pi$
decays by simple $W^+$ exchange.
For $K^+\to\pi^+\pi^0$ we can draw the two Feynman diagrams
of Fig.~\ref{figkpipiplus}(a).
\begin{figure}[t]
\begin{minipage}{0.66\textwidth}
\centerline{\includegraphics[width=0.99\textwidth]{figKPOdecay.ps}}
\centerline{(a)}
\end{minipage}
\begin{minipage}{0.33\textwidth}
\centerline{\includegraphics[width=0.99\textwidth]{figKOOdecay.ps}}
\centerline{(b)}
\end{minipage}
\caption{(a) The two naive $W^+$-exchange diagrams for
$K^+\longrightarrow \pi^+\pi^0$.
(b) No simple $W^+$-exchange diagram is possible for
$K^0\longrightarrow \pi^0\pi^0$.}
\label{figkpipiplus}
\end{figure}
The $W^+$-hadron couplings are known from semi-leptonic
decays. This approximation agrees with the
measured decay within a factor of two.


A much worse result appears when we try the same for
$K^0\to\pi^0\pi^0$. As shown in Fig. \ref{figkpipiplus}(b) there is no
possibility to draw diagrams similar to those in Fig. \ref{figkpipiplus}(a).
The needed vertices always violate charge-conservation.
So we expect that the neutral decay should be small compared
with the ones with charged pions. Well, if we look at the experimental
results, we see
\ba
\Gamma(K^0\longrightarrow\pi^0\pi^0)=&
\dsp\frac{1}{2}\Gamma(K_S\longrightarrow\pi^0\pi^0)&=2.3\cdot10^{-12}\mbox{ MeV}
\nonumber\\
\Gamma(K^+\longrightarrow\pi^+\pi^0)=&\dsp
1.1\cdot10^{-14}\mbox{ MeV}&
\ea
So the expected zero one is by far the largest !!!

The same conundrum can be expressed in terms of the
 isospin amplitudes:\,\footnote{The  sign
convention is the one used in the work by J.~Prades and myself.}
\ba
A[K^0 \to \pi^0 \pi^0]  &\equiv& \sqrt{{1}/{3}} { A_0}
-\sqrt{{2}/{3}} \, { A_2} \nonumber\\ 
A[K^0 \to \pi^+ \pi^-]  &\equiv& \sqrt{{1}/{3}} { A_0}
+\sqrt{{1}/{6}} \, { A_2} \nonumber\\ 
A[K^+ \to \pi^+ \pi^0]  &\equiv& ({\sqrt{3}}/{2}) { A_2}\,.
\ea
The above quoted experimental results can now be rewritten as
\be 
\label{dIhalf1}
\left|{A_0}/{A_2}\right|_{\mbox{exp}} = 22
\ee
while the naive $W^+$-exchange discussed would give
\be
\label{dInaive}
\left|{A_0}/{A_2}\right|_{\mbox{naive}} = \sqrt{2}\,.
\ee
This discrepancy is known as the problem of the $\Delta I=1/2$ rule.

Some enhancement comes from final state $\pi\pi$-rescattering.
Removing these and higher order effects in the light quark masses
one obtains~\cite{KMW,BDP}
\be
\label{dIhalf2}
\dsp\left|{A_0}/{A_2}\right|_{\chi} = 17.8\,.
\ee
This changes the discrepancy somewhat but is still different by an order
of magnitude from the naive result (\ref{dInaive}). The difference
will have to be explained by pure strong interaction effects
and it is a {\em qualitative} change, not just a quantitative one.

We also use amplitudes without the final state interaction
phase:
\be
A_I=-ia_I e^{i\delta_I}
\ee
for $I=0,2$. $\delta_I$ is the angular momentum zero, isospin I scattering
phase at the Kaon mass.

\section{$K\to\pi\pi$, $\varepsilon$, \epspeps}
\label{epspeps}


The $K^0$, $\kob$\ states have $\bar sd$,  $\bar ds$
quark content. $CP$ acts on these states as
\be
CP |K^0\rangle = -|\kob\rangle\,.
\ee
{}We can construct eigenstates with a definite $CP$ transformation:
\be
K^0_{1(2)} = \frac{1}{\sqrt{2}}
\left(K^0 -(+) \kob\right),\quad\quad
 CP|K_{1(2)} = +(-)|K_{1(2)}\,.
\ee
The main decay mode of $K^0$-like states is $\pi\pi$. A two pion state
with charge zero in spin zero is always CP even.
Therefore the decay $K_1\to\pi\pi$ is possible
but $K_2\to\pi\pi$ is {\em impossible}; $K_2\to\pi\pi\pi$ is possible.
Phase-space for the $\pi\pi$ decay is much larger than for the
three-pion final state. Therefore if we start out with a pure $K^0$ or $\kob$\
state, the
$K_2$ component in its wave-function lives much longer than the $K_1$
component such that after a long time only the $K_2$ component
survives.

In the early sixties, as you see it pays off to do precise experiments,
one actually measured \cite{CCFT}
\be
\frac{\Gamma(K_L\to\pi^+\pi^-)}{\Gamma(K_L\to\mbox{all})}=
(2\pm0.4)\cdot10^{-3}\ne 0\,,
\ee
showing that {\em $CP$ is violated}\,.
This leaves us with the questions:
\begin{itemize}
\item[{\bf ???}] Does $K_1$ turn in to $K_2$
(mixing or indirect $CP$ violation)?
\item[{\bf ???}] Does $K_2$ decay directly into $\pi\pi$ 
(direct $CP$ violation)?
\end{itemize}
In fact, the 
answer to both is {\em YES}\, and is major qualitative test
of the standard model Higgs-fermion sector and the $CKM$-picture
of $CP$-violation.

The presence of $CP$-violation means that $K_1$ and $K_2$ are not the mass
eigenstates, these are
\be
K_{S(L)} = \frac{1}{\sqrt{1+|\tilde\varepsilon|^2}}
\left(K_{1(2)}+\tilde\varepsilon K_{2(1)}\right)\,.
\ee
They are not orthogonal since the Hamiltonian is not hermitian.


We define the observables
\ba
\varepsilon &\equiv& \frac{A(K_L\to(\pi\pi)_{I=0})}{A(K_S\to(\pi\pi)_{I=0})}
\nonumber\\
\varepsilon^\prime &=& \frac{1}{\sqrt{2}}\left(
\frac{A(K_L\to(\pi\pi)_{I=2})}{A(K_S\to(\pi\pi)_{I=0})}-
\varepsilon
\frac{A(K_S\to(\pi\pi)_{I=2})}{A(K_S\to(\pi\pi)_{I=0})}\right)\,.
\ea
The latter has been specifically constructed to remove the 
$K^0$-$\kob$\ transition.
$|\varepsilon|$ is a directly measurable as ratios of decay rates.

We now make a series of experimentally valid approximations,
\be
 |\im a_0|,|\im a_2| << |\re a_2| << |\re a_0|,\quad
 |\varepsilon|,|\tilde \varepsilon| << 1,
 \quad |\varepsilon^\prime| << | \varepsilon|\,,
\ee
to obtain the usually quoted expression
\be
\label{epsp}
\varepsilon^\prime = \frac{i}{\sqrt{2}}e^{i(\delta_2-\delta_0)}
\frac{\re a_2}{\re a_0}\left(\frac{\im a_2}{\re a_2}-\frac{\im a_0}{\re a_0}
\right)\,.
\ee


Experimentally,\cite{PDG}
\be
|\varepsilon| = (2.271\pm0.017)\cdot 10^{-3}\,.
\ee

The set of diagrams, depicted schematically in Fig.~\ref{figbox}(a),
responsible for $K^0\kob$\ mixing are known as box diagrams.
It is the presence of the virtual intermediate quark lines
of up, charm and top quarks that produces the $CP$-violation. 
\begin{figure}[t]
\begin{minipage}{0.45\textwidth}
\centerline{\includegraphics[width=0.99\textwidth]{figbox.ps}}
\centerline{(a)}
\end{minipage}
\begin{minipage}{0.54\textwidth}
\centerline{\includegraphics[width=0.95\textwidth]{figpenguin.ps}}
\centerline{(b)}
\end{minipage}
\caption{\label{figbox}
(a) The box diagram contribution to $K^0\kob$\ mixing.
Crossed versions and diagrams with extra gluons etc. are not shown.
(b) The Penguin diagram contribution
to $K\to\pi\pi$.
Extra gluons and crossed versions etc. are not shown.}
\end{figure}

The experimental situation on $\epspeps$ was unclear for a long time.
Two large experiments, NA31 at CERN and E731 at FNAL,
obtained conflicting results
in the mid 1980's. Both groups have since gone on and build improved versions
of their detectors, NA48 at CERN and KTeV at FNAL.
$\epspeps$ is
measured via the double ratio
\be
\re\left(\frac{\varepsilon^\prime}{\varepsilon}\right)
= \frac{1}{6}\left\{
1-\frac{\dsp \Gamma(K_L\to\pi^+\pi^-)/\Gamma(K_S\to\pi^+\pi^-)}
{\dsp \Gamma(K_L\to\pi^0\pi^0)/\Gamma(K_S\to\pi^0\pi^0)}
\right\}\,.
\ee
The two main experiments follow a somewhat different strategy in
measuring this double ratio, mainly in the way the relative normalisation
of $K_L$ and $K_S$ components is treated.
After some initial disagreement with the first results, KTeV has
reanalysed their systematic errors and the
situation for $\epspeps$ is now quite clear.
We show the recent results in Table \ref{tabepspeps}.
The data are taken from Ref.~\citebk{epspeps} and the recent
reviews in the Lepton-Photon conference.\cite{Kessler,NA4801} 
\begin{table}[t]
\tbl{Recent results on $\varepsilon^\prime/\varepsilon$. The years
refer to the data sets.}
{\footnotesize
%\begin{center}
{\begin{tabular}{|cc|}
\hline
NA31 & $(23.0\pm6.5)\times 10^{-4}$\\
E731 & $(7.4\pm5.9)\times 10^{-4}$\\
\hline
KTeV 96 & $(23.2\pm4.4)\times10^{-4}$\\
KTeV 97 & $(19.8\pm2.9)\times10^{-4}$\\
NA48 97 & $(18.5\pm7.3)\times10^{-4}$\\
NA48 98+99 & $(15.0\pm2.7)\times10^{-4}$\\
\hline
ALL & $(17.2\pm1.8)\times10^{-4}$\\
\hline
\end{tabular}}}
%\end{center}
\label{tabepspeps}
\end{table}


The Penguin diagram shown in Fig. \ref{figbox}(b) contributes to the direct
$CP$-violation as given by $\varepsilon^\prime$.
Again, $W$-couplings to all three generations show up so $CP$-violation
is possible in $K\to\pi\pi$. This is a qualitative prediction of the
standard model and borne out by experiment.
The main problem is now to embed these diagrams and the simple $W$-exchange
in the full strong interaction. The $\Delta I=1/2$ rule shows that there
will have to be large corrections to the naive picture.

\section{From Quarks to Mesons: a Chain of Effective Field Theories}
\label{EFT}

The full calculation in the presence of the strong interaction is quite
difficult. Even at short distances, due to the presence of logarithms
of large ratios of scales, a simple one-loop calculation gives
very large effects. These need to be resummed which fortunately can be done
using renormalisation group methods. 

The three steps of the
full calculation are depicted in 
Fig. \ref{figsteps}.
\begin{figure}
\begin{tabular}{ccc}
ENERGY SCALE & FIELDS & Effective Theory\\
\hline
\\[2mm]
$M_W$ & \framebox{\parbox{5cm}{\begin{center}
$W,Z,\gamma,g$;\\$\tau,\mu,e,\nu_\ell$;\\
$t,b,c,s,u,d$\end{center}
}} & Standard Model\\[2mm]
 & {\em $\Downarrow$ using OPE} & \\[2mm]
$\lesssim m_c$ & \framebox{\parbox{5cm}{\begin{center}
$\gamma,g$; $\mu,e,\nu_\ell$;\\ $s,d,u$\end{center}
}} & \parbox{2cm}{QCD,QED,\\${\mathcal H}_{\mbox{eff}}^{|\Delta S|=1,2}$}\\[2mm]
  & {\em$\Downarrow$ ???} & \\[2mm]
$M_K$  & \framebox{\parbox{5cm}{\begin{center}$\gamma$; $\mu,e,\nu_\ell$;\\
$\pi$, $K$, $\eta$\end{center}}} & CHPT\\ \\\hline
\end{tabular}
\caption{A schematic exposition of the various steps in the calculation
of nonleptonic matrix-elements.}
\label{figsteps}
\end{figure}
First we integrated out the heaviest particles step by step using Operator
Product Expansion methods. The steps OPE we describe in the next subsections
while step {\em ???} we will split up in more subparts later.


\subsection{Step I: from SM to OPE}

The first step concerns the
standard model diagrams of Fig. \ref{figSM}(a).
\begin{figure}[t]
\begin{minipage}{0.495\textwidth}
\centerline{\includegraphics[width=0.999\textwidth]{figSM.ps}}
\centerline{(a)}
\end{minipage}
\begin{minipage}{0.495\textwidth}
\centerline{\includegraphics[width=0.999\textwidth]{figQ2.ps}}
\centerline{(b)}
\end{minipage}
\caption{(a) The standard model diagrams to be calculated at a high scale.
(b) The diagrams needed for the
matrix-elements calculated at a scale $\mu_H\approx m_W$
using the effective Hamiltonian.}
\label{figSM}
\end{figure}
We replace their effect
with a contribution of an effective Hamiltonian given by
\be
\label{defHeff}
 {\mathcal H}_{\mbox{eff}} = \sum_i C_i(\mu) Q_i(\mu)
= \frac{G_F}{\sqrt{2}} V_{ud}V_{us}^*
\sum_i\left(z_i-y_i\frac{V_{td}V_{ts}^*}{V_{ud}V_{us}^*}\right)Q_i\,.
\ee
In the last part we have real coefficients $z_i$ and $y_i$ and
the CKM-matrix-elements occurring are shown explicitly.
The four-quark operators $Q_i$ can be found in e.g. Ref.~\citebk{BPscheme}.

We calculate now matrix-elements between quarks and gluons in the
standard model using the diagrams of Fig. \ref{figSM}(a)
and equate those to the same matrix-elements calculated using the
effective Hamiltonian of Eq. (\ref{defHeff}) and the diagrams
of Fig. \ref{figSM}(b). This determines the value of the $z_i$ and $y_i$.
The top quark and the $W$ and $Z$ bosons are integrated out all at the
same time. There should be no large logarithms present due to that.
The scale $\mu = \mu_H$ in the diagrams of Fig.~\ref{figSM}(b)
of the OPE expansion diagrams
should be chosen of the order of the $W$ mass.
The scale $\mu_W$ in the Standard Model diagrams of Fig.~\ref{figSM}(a)
should be chosen of the same order.

{\underline{Notes}:}\\
$\bullet$ In the Penguin diagrams
$CP$-violation shows up since all 3 generations
are present.\\
$\bullet$ The equivalence is done by calculating matrix-elements between
{\em Quarks and Gluons}\\
$\bullet$ The SM part is $\mu_W$-independent to $\alpha_S^2(\mu_W)$.\\
$\bullet$ OPE part: The $\mu_H$ dependence of
$C_i(\mu_H)$ cancels the $\mu_H$ dependence of the diagrams
to order $\alpha_S^2(\mu_H)$.

This procedure gives at $\mu_W = \mu_H=M_W$
in the NDR-scheme~\footnote{The precise definition
of the four-quark operators $Q_i$ comes in here as well. See the lectures by
Buras \cite{Buras1} for a more extensive description of that.}
the numerical values given in Table \ref{tabzimw}.
\begin{table}[t]
\tbl{The Wilson coefficients and their main source at
the scale $\mu_H=m_W$ in the NDR-scheme.}
{\footnotesize
\begin{tabular}{|ccc|ccc|}
\hline
$z_1$   & 0.053     &  { $g,\gamma$-box}               &
$y_6$	&$-$0.0019  & {$g$-Penguin}\\
$z_2$	& 0.981	    &{ $W^+$-exchange  $g,\gamma$-box} &
$y_7$	& 0.0009    & {$\gamma,Z$-Penguin}\\
$y_3$	& 0.0014    & {$g,Z$-Penguin  $WW$-box}        &
$y_8$	& 0.	    & \\
$y_4$	&$-$0.0019  & {$g$-Penguin}                    &
$y_9$	& $-$0.0074 & {$\gamma,Z$-Penguin  $WW$-box}\\
$y_5$	& 0.0006    & {$g$-Penguin}                    &
$y_{10}$& 0.        & \\
\hline
\end{tabular}
}
\label{tabzimw}
\end{table}
In the same table I have
given the main source of these numbers. Pure tree-level $W$-exchange
would have only given $z_2=1$ and all others zero.
Note that the coefficients from $\gamma,Z$ exchange are similar to the gluon
exchange ones since $\alpha_S$ at this scale is not very big.


\subsection{Step II}
\label{shortdistance}

Now comes the main advantage of the OPE formalism. Using the
renormalisation group equations we can  calculate
the change with $\mu$ of the $C_i$, thus resumming the
$\log\left(m_W^2/\mu^2\right)$ effects.
The renormalisation group equations (RGEs) for the
strong coupling and the Wilson coefficients are
\be
\label{RGE}
\mu\frac{d}{d\mu}g_S(\mu) = \beta(g_S(\mu)),\quad
\mu\frac{d}{d\mu}C_i(\mu) = \gamma_{ji}(g_S(\mu),\alpha) C_j(\mu)\,.
\ee
$\beta$ is the QCD beta function for the running coupling.
The coefficients
$\gamma_{ij}$ are the elements
of the anomalous dimension matrix $\hat\gamma$. They can be derived from
the infinite parts of loop diagrams and this has been done to
one~\cite{one-loop} and two loops.\cite{two-loop}
The series in $\alpha$ and $\alpha_S$ is known to
\be
\hat\gamma = \hat\gamma^0_S \frac{\alpha_S}{4\pi} 
            + \hat\gamma^1_S \left(\frac{\alpha_S}{4\pi}\right)^2
            + \hat\gamma_e \frac{\alpha}{4\pi}
            + \hat\gamma_{se}\frac{\alpha_S}{4\pi} \frac{\alpha}{4\pi}
            +\cdots 
\ee
Many subtleties are involved in this calculation.\cite{Buras1,two-loop}
They all are related to the fact that everything at higher loop orders need
to be specified correctly, and many things which are equal at tree
level are no longer so in $d\ne4$ and at higher loops, see the
lectures~\citebk{Buras1} or the review~\citebk{twoloopreview}.
The numbers below are obtained by numerically
integrating Eq. (\ref{RGE}).\cite{BPdIhalf,BPeps}

We perform the following steps to get down to a scale $\mu_{OPE}$
 around 1~GeV. Starting from the $z_i$ and $y_i$
at the scale {\boldmath$\mu_H$}:\\
(1)
solve Eqs. (\ref{RGE}); run from {\boldmath${\mu_H}$}
to {\boldmath${\mu_b}$}.\\
(2)
At {\boldmath$\mu_b$}$\approx m_b$ remove $b$-quark and match to the
theory without $b$ by calculating matrix-elements of the
effective Hamiltonian in the five and in the four-quark picture
and putting them equal.\\
(3)
Run from {\boldmath$\mu_b$} to {\boldmath$\mu_c$}$\approx m_c$.\\
(4)
At {\boldmath$\mu_c$} remove the $c$-quark and match to
the theory without $c$.\\
(5)
Run from {\boldmath$\mu_c$} to {\boldmath$\mu_{\mbox{OPE}}$}.\\
Then {\em all}\, large logarithms including $m_W$, $m_Z$, $m_t$,
$m_b$ and $m_c$, are summed.

With the inputs
$m_t(m_t)=166~GeV$, $\alpha=1/137.0$, $\alpha_S(m_Z) = 0.1186$
which led to the initial conditions shown in Table~\ref{tabzimw}, we can
perform the above procedure down to $\mu_{OPE}$.
Results for 900~MeV are shown in columns two and three of Table \ref{tabzimu}.
\begin{table}[t]
\tbl{The Wilson coefficients
$z_i$ and $y_i$ at a scale $\mu_{\mbox{OPE}}=$ 900~MeV
in the NDR scheme and in the $X$-boson scheme at $\mu_X =$ 900~MeV.}
{\footnotesize
\begin{tabular}{|c|cc|cc|}
\hline
i        &\hspace{6mm} $z_i$\hspace{6mm} &\hspace{6mm} $y_i$\hspace{6mm} &\hspace{6mm} $z_i$\hspace{6mm} &\hspace{6mm} $y_i$\hspace{6mm} \\
\hspace{1cm}         &\multicolumn{2}{c}{$\mu_{OPE} = 0.9$ GeV}\vline
&\multicolumn{2}{c}{$\mu_X =0.9$ GeV}\vline\\
\hline
$z_{1}$ &$-$0.490           & 0.                & $-$0.788        & 0.\\
$z_{2}$	&1.266 		    & 0.                & 1.457           & 0. \\
$z_{3}$	&0.0092		    &  0.0287           & 0.0086          & 0.0399\\
$z_{4}$	&$-$0.0265	    & $-$0.0532         & $-$0.0101       & $-$0.0572\\
$z_{5}$	& 0.0065	    & 0.0018            & 0.0029          & 0.0112\\
$z_{6}$	&$-$0.0270	    & $-$0.0995         &  $-$0.0149      & $-$0.1223\\
$z_{7}$	& 2.6$~10^{-5}$	    &$-$0.9$~10^{-5}$   & 0.0002          &$-$0.00016\\
$z_{8}$	& 5.3$~10^{-5}$	    & 0.0013            &  6.8$~10^{-5}$  & 0.0018\\
$z_{9}$	& 5.3$~10^{-5}$	    & $-$0.0105         &  0.0003         & $-$0.0121\\
$z_{10}$&$-$3.6$~10^{-5}$   &0.0041             & $-$8.7$~10^{-5}$& 0.0065\\
\hline
\end{tabular}
}
\label{tabzimu}
\end{table}
$z_1$ and $z_2$ have changed much from $0$ and $1$.
This is the short-distance contribution to the $\Delta I=1/2$ rule.
We also see a large enhancement of $y_6$ and $y_8$, which will
lead to our value of $\varepsilon^\prime$.


\subsection{Step III: Matrix-elements}
\label{longdistance1}

Now
remember that the $C_i$ depend on $\mu_{OPE}$ (scale dependence)
and on the definition of the $Q_i$ (scheme dependence)
and the numerical change in the coefficients
due to the various choices for the $Q_i$ possible is not negligible.
It is therefore important both from the phenomenological and fundamental
point of view that this dependence is correctly accounted for in
the evaluation of the matrix-elements.
We can solve this in various ways.\\
$\bullet$ {\bf Stay in QCD} $\Rightarrow$
Lattice calculations.\cite{Sachrajda}\\
$\bullet$ {\bf ITEP Sum Rules} or QCD sum rules.~\cite{SVZ}\\
$\bullet$ {\bf Give up} $\Rightarrow$ Naive factorisation.\\
$\bullet$ {\bf Improved factorisation}\\
$\bullet$ {\bf  $X$-boson method} (or fictitious Higgs method)\\
$\bullet$ {\bf  Large $N_c$} (in combination with something like
the $X$-boson method.) Here the difference is mainly in the
treatment of the low-energy hadronic physics. Three main approaches
exist of increasing sophistication.\footnote{Which of course means that
calculations exist only for simpler matrix-elements for the more sophisticated
approaches.}
\begin{flushright}
\begin{minipage}{0.97\textwidth}
$\circledast$ CHPT: As originally proposed by Bardeen-Buras-G\'erard~\cite{BBG}
and now pursued mainly by Hambye and collaborators.\cite{Hambye}\\
$\circledast$ ENJL (or extended Nambu-Jona-Lasinio model~\cite{ENJL}):
As mainly done by
myself and J.~Prades.\cite{BPBK,BPscheme,kptokpp,BPdIhalf,BPeps}\\
$\circledast$ LMD or lowest meson dominance approach.\cite{LMD} 
These papers stay
with dimensional regularisation throughout. The $X$-boson corrections
discussed below, show up here as part of the QCD corrections.
\end{minipage}
\end{flushright}
$\bullet$ {\bf Dispersive methods} Some matrix-elements can in principle
be deduced
from experimental spectral functions.

Notice that there other approaches as well, e.g.~the chiral quark
model.\cite{PichdeRafael} 
These have no underlying arguments why the $\mu$-dependence
should cancel, but the importance of several effects
was first discussed in this context. I will also not treat the
calculations done using bag models and potential models which
similarly do not address the $\mu$-dependence issue.

\section{The $X$-boson Method and Results using ENJL for the Long Distance}
\label{Xboson}

We want to have a consistent calculational scheme that takes the scale
and scheme dependence into account correctly. Let us therefore have a closer
look at how we calculate the matrix-elements using naive factorisation.
We start from the four-quark operator:\\
$\Rightarrow$ See it as a product of currents or densities.\\
$\Rightarrow$ Evaluate current matrix-elements in low energy
theory or model or from experiment.\\
$\Rightarrow$ Neglect extra momentum transfer between the current matrix
elements.\\
The main lesson here is that {\em currents and densities are easier
to deal with.} We also need to go beyond the approximation in the last step.
 To obtain well defined currents, we replace the four-quark operators
by exchanges of fictitious massive $X$-bosons
coupling to two-quark currents or densities.
\be
{\mathcal H} = \sum_i C_i(\mu_{OPE}) Q_i
\quad
\Longrightarrow
\quad
\sum_i g_i X_i J_i\,.
\ee
$\bullet$ This is a well defined scheme of nonlocal operators.\\
$\bullet$ The matching to obtain the coupling constants
$g_i$ from the $C_i$ is done with
matrix-elements of {\em quarks and gluons}.

A simple example is the one needed for the $B_K$ parameter. The
four-quark operator is replaced by the exchange of one $X$-boson $X_B$:
\be
C(\mu) \bar d\gamma_\mu(1-\gamma_5)s\,\bar d\gamma_\mu(1-\gamma_5)s
\quad
\Longrightarrow
\quad
g_B X_B^\mu \bar d\gamma_\mu(1-\gamma_5)s\,.
\ee
Taking a matrix-element between quarks at
next-to-leading order in $\alpha_S$
gives
\be
\label{Xmatching}
C(\mu_{OPE}) \left(1+ r\alpha_S(\mu_{OPE})\right) =
(g_B^2/M_{X_B}^2) \left(1+r'\alpha_S+r''\log\left(M_{X_B}^2/\mu^2\right)\right)
\,.
\ee
The coefficients $r$ and $r'$ take care of the scheme dependence.
The l.h.s. is scale independent to the required order in $\alpha_S$.
The effect of these coefficients surprisingly always went in the
direction to improve agreement with experiment~\cite{BPscheme,BPeps}
as can be seen from columns 4 and 5 in Table~\ref{tabzimu}.

The final step is the matrix-element of $X_B$-boson exchange.
For this we split the integral over the  $X_B$ momentum $q_X$ in two parts
\be
\int dq_X^2 \Longrightarrow
\int_0^{\mu^2} dq_X^2 +\int_{\mu^2}^\infty dq_X^2\,.
\ee
\begin{figure}[t]
\centerline{\includegraphics[width=0.8\textwidth]{figX.eps}}
\caption{(a) The leading in $1/N_c$ contribution from $X_B$ exchange.
(b) The large momentum part of the $X_B$ exchange matrix-element.}
\label{figX}
\end{figure}
The leading in $N_c$ contribution is depicted in Fig.~\ref{figX}(a)
and corresponds to the large $N_c$ factorisation.
The large momentum regime is evaluated by the diagram in Fig.~\ref{figX}(b),
since the large momentum must flow back through quarks and gluons.
Hadronic exchanges are power suppressed because of the form factors involved.
The $\alpha_S$ present already suppresses by $N_c$ so the matrix-element
of this part can be evaluated using factorisation.
This part cancels the $r''\log(M_{X_B}^2/\mu^2)$ present
in (\ref{Xmatching}). The final part with small $q_X$ momentum in the
integral then needs to be evaluated nonperturbatively.
Here one can use Chiral Perturbation Theory, the ENJL model
or meson exchange approximations with various short-distance constraints.
 
Let me now show some results from Refs~\citebk{BPscheme,BPeps}.
The chiral limit coupling $G_8$ responsible for the octet contribution,
it is 1 in the naive approximation and about 6 when fitted to
experiment,\cite{KMW,BDP} is shown in Fig.~\ref{figG8}.
\begin{figure}[t]
\centerline{\includegraphics[angle=270,width=0.495\textwidth]{G8.real.ps}
\includegraphics[angle=270,width=0.495\textwidth]{G8.imag.ps}
}
\caption{Left: the results for the real part of the octet $\re G_8$
as a function of $\mu$. Right: the same for the imaginary part.}
\label{figG8}
\end{figure}
As can be seen, the matching between the short and long distance is reasonable
both for the real and imaginary parts.
The value of $\re G_8$ is dominated by the matrix-element of $Q_1$ and $Q_2$
but about 30-60\% comes from the long distance Penguin part of $Q_2$.
The result for the $\im G_8$ corresponds to a value of the $Q_6$ 
matrix-element much larger than usually assumed. We obtained $B_6\approx2$-2.5
while it is usually {\em assumed} to be less than 1.5.

Putting our results in (\ref{epsp}) we obtain a chiral limit value
for $(\epspeps)_\chi$ of about $6\cdot 10^{-3}$.

We now add the main isospin breaking component $\Omega$~\cite{EPetal}
and the effect of final state interaction (FSI).\cite{PP}
The latter in our case has mainly effect on the forefactor $\re a_2/\re a_0$
in (\ref{epsp}) since the ratios of imaginary parts have been evaluated
to the same order in $p^2$ in CHPT and thus receive no FSI corrections.
The final result is~\cite{BPeps} $\epspeps\approx 1.5\cdot 10^{-3}$
with an error $\gtrapprox$ 50\%.



\section{Dispersive Estimates for $\langle Q_7\rangle$ and
$\langle Q_8\rangle$}
\label{dispersive}

Some of the matrix-elements we want can be extracted from
experimental information in a different way. The canonical
example is the mass difference between the charged and the neutral pion
in the chiral limit which can be extracted from a dispersive integral
over the difference of the vector and axial vector spectral
functions.\cite{Dasetal}

This idea has been pursued in the context of weak decay in a series of
papers by Donoghue, Golowich and collaborators.\cite{DG1}
The matrix-element of $Q_7$ could be extracted directly
from these data. To get at the matrix-element of $Q_8$ is somewhat more
difficult. Ref. \citebk{DG1} extracted it first
by requiring $\mu$-independence, this corresponds
to extracting  the matrix
element of $Q_8$ from the spectral functions
via the coefficient of the dimension 6 term in the operator
product expansion of the underlying Green's function.
The most recent papers using this method are
Refs.~\citebk{KnechtQ7Q8,NAR01,DG2} and \citebk{BPQ7Q8}.
In the last two papers also some QCD corrections were included which had a
substantial impact on the numerical results.

The results are given in Table.~\ref{tableQ7Q8}. The operator
$O_6^{(1)}$ is related by a chiral transformation to $Q_7$
and $O_6^{(2)}$ to $Q_8$.
The numbers are valid in the chiral limit.
\begin{table}[t]
\tbl{The values of the VEVs in the NDR scheme at $\mu_R=2$ GeV.
The most recent dispersive results are line 3 to 5.
The other results are shown for comparison. Errors are those quoted in the
papers.
Adapted from Ref.~[71].
}
{\footnotesize
\begin{tabular}{|c|c|c|}
\hline
Reference&$\langle 0| O_6^{(1)}|0 \rangle^{NDR}_\chi$&
$\langle 0| O_6^{(2)}|0 \rangle^{NDR}_\chi$ 
 \\ 
\hline
$B_7=B_8=1$                      & $-(5.4\pm2.2)\cdot10^{-5}$ GeV$^6$ & $(1.0 \pm 0.4)\cdot10^{-3}$ GeV$^6$   \\
\hline				   
Bijnens et al. \citebkcap{BPQ7Q8}     &$-(4.0\pm0.5)\cdot10^{-5}$ GeV$^6$& $(1.2  \pm 0.5)\cdot10^{-3}$ GeV$^6$  \\  
Knecht et al. \citebkcap{KnechtQ7Q8}  &$-(1.9\pm0.6)\cdot10^{-5}$ GeV$^6$& $(2.3 \pm 0.7)\cdot10^{-3}$ GeV$^6$   \\   
Cirigliano et al. \citebkcap{DG2}     &$-(2.7\pm1.7)\cdot10^{-5}$ GeV$^6$& $(2.2 \pm 0.7)\cdot10^{-3}$ GeV$^6$   \\  
\hline				   
Donoghue et al.\citebkcap{DG1}        &$-(4.3\pm0.9)\cdot10^{-5}$ GeV$^6$& $(1.5 \pm 0.4)\cdot10^{-3}$ GeV$^6$      \\  
Narison        \citebkcap{NAR01}      &$-(3.5\pm1.0)\cdot10^{-5}$ GeV$^6$& $(1.5 \pm 0.3)\cdot10^{-3}$ GeV$^6$      \\  
lattice   \citebkcap{Q7Q8lattice}     &$-(2.6\pm0.7)\cdot10^{-5}$ GeV$^6$& $(0.74 \pm 0.15)\cdot10^{-3}$ GeV$^6$    \\
ENJL  \citebkcap{BPeps}               &$-(4.3\pm0.5)\cdot10^{-5}$ GeV$^6$& $(1.3\pm0.2)\cdot10^{-3}$ GeV$^6$\\
\hline
\end{tabular}
}
\label{tableQ7Q8}
\end{table}
The various results for the matrix-element of $O_6^{(1)}$ are in reasonable
agreement with each other. The underlying spectral integral,
evaluated directly from data in Refs.~\citebk{NAR01},\citebk{DG2} and
 \citebk{BPQ7Q8},
or via the minimal hadronic ansatz~\citebk{KnechtQ7Q8} are in
better agreement. 
The largest source of the differences is the way the different
results for the underlying evaluation of $O_6^{(2)}$ come back into
$O_6^{(1)}$.

The results for $O_6^{(2)}$ are also in
reasonable agreement. Ref. \citebk{BPQ7Q8}
uses two approaches. First, the matrix-element for $O_6^{(2)}$ can be extracted
via a similar dispersive integral over the scalar and pseudoscalar spectral
functions. The requirements of short-distance matching for this spectral
function combined with a saturation with a few states imposes that
the nonfactorisable part is suppressed and the number and error quoted follows
from this. Extracting the coefficient of the dimension 6 operator
in the expansion of the vector and axial-vector spectral functions yields a
result comparable but with a larger error of about 0.9.
Ref. \citebk{KnechtQ7Q8} uses a derivation
based on a single resonance plus continuum
ansatz for the spectral functions and assumes a typical large $N_c$
error of 30\%. This ansatz worked well for lower moments of the spectral
functions which can be tested experimentally. Adding more resonances allows
for a broader range of results.\cite{BPQ7Q8}
Ref.~\citebk{DG2} chose to enforce all the known constraints on the
vector and axial-vector spectral functions to obtain a result. This resulted
in rather large cancellations between the various contributions making an error
analysis more difficult. A reasonable estimate lead to the value quoted.

The reason why the central value based on the same data can be so different is
that the quantity in question is sensitive to the energy regime
above 1.3~GeV where the accuracy of the data is rather low.

\section{Conclusions}
\label{conclusions}

Penguins are alive and well, they provide a sizable part of the $\Delta I=1/2$
enhancement though mainly through long distance Penguin like
topologies in the evaluation of the matrix-element of $Q_2$.
They have found a much richer use in the $CP$ violation phenomenology.
For the electroweak Penguins, calculations are in qualitative agreement
but more work is still needed to get the errors down. For the strong Penguins,
the work I have presented here shows a strong enhancement over factorisation
with $B_6$ significantly larger than one. The latter conclusion is similar
to the one derived from the older more phenomenological arguments where
the coefficients were taken at a low scale and the matrix-elements
for $Q_6$ taken from the the value of the $\Delta I=1/2$ rule.
This also indicated a rather large enhancement of the matrix-element of $Q_6$
over the naive factorisation.

\section*{Acknowledgements}
This work has been partially supported  by the Swedish Research Council
and by the European Union TMR Network
EURODAPHNE (Contract No. ERBFMX-CT98-0169).
I thank the organisers for a nice and well organised meeting and Arkady for
many discussions and providing a good reason to organise a meeting.
\begin{thebibliography}{99}

\bibitem{Penguin3}
A.~I.~Vainshtein,
the 1999 Sakurai Prize Lecture,
%``How penguins started to fly,''
{\it Int.\ J.\ Mod.\ Phys.\ }{\bf A14}, 4705 (1999)
.
%%CITATION = ;%%

\bibitem{Isidori}
G.~Isidori,
%``Standard model vs new physics in rare kaon decays,''
 talk KAON99;
%%CITATION = ;%%
G.~Isidori,
%``Rare decays: Theory vs. experiments,''
;
%%CITATION = ;%%
G.~Buchalla,
%``Rare kaon decays: Overview,''
.
%%CITATION = ;%%

\bibitem{Littenberg}
L.~Littenberg and G.~Valencia,
%``Rare and radiative kaon decays,''
{\it Ann.\ Rev.\ Nucl.\ Part.\ Sci.\  }{\bf 43}, 729 (1993)
.
%%CITATION = ;%%

\bibitem{Buras2}
A.~J.~Buras,
``CP violation and rare decays of K and B mesons,''
 Lake Louise lectures.
%%CITATION = ;%%

\bibitem{QCDweak}
J.~Bijnens,
%``QCD and weak interactions of light quarks,''

%%CITATION = ;%%
to be published in 'At the Frontier of Particle Physics/ Handbook of QCD', 
edited by M. Shifman, Volume 4. 

\bibitem{Fermi}
E.~Fermi,
%``Trends To A Theory Of Beta Radiation. (In Italian),''
{\it Nuovo Cim.\  }{\bf 11}, 1 (1934);
%%CITATION = NUCIA,11,1;%%
%E.~Fermi,
%``An Attempt Of A Theory Of Beta Radiation. 1,''
{\it Z.\ Phys.\  }{\bf 88}, 161 (1934).
%%CITATION = ZEPYA,88,161;%%

\bibitem{LeeYang}
T.~D.~Lee and C.~N.~Yang,
%``Question Of Parity Conservation In Weak Interactions,''
{\it Phys.\ Rev.\  }{\bf 104}, 254 (1956).
%%CITATION = PHRVA,104,254;%%

\bibitem{CSWu}
C.~S.~Wu {\it et al.}, 
%E.~Ambler, R.~W.~Hayward, D.~D.~Hoppes and R.~P.~Hudson,
%``Experimental Test Of Parity Conservation In Beta Decay,''
{\it Phys.\ Rev.\  }{\bf 105}, 1413 (1957).
%%CITATION = PHRVA,105,1413;%%

\bibitem{FriedmanTelegdi}
J.~I.~Friedman and V.~L.~Telegdi,
%``Nuclear Emulsion Evidence For Parity Nonconservation In The Decay Chain Pi+ Mu+ E+,''
{\it Phys.\ Rev.\  }{\bf 105}, 1681 (1957).
%%CITATION = PHRVA,105,1681;%%


\bibitem{SudarshanMarshak}
E.~C.~Sudarshan and R.~E.~Marshak,
%``Chirality Invariance And The Universal Fermi Interaction,''
{\it Phys.\ Rev.\  }{\bf 109}, 1860 (1958).
%%CITATION = PHRVA,109,1860;%%

\bibitem{FeynmanGell-Mann}
R.~P.~Feynman and M.~Gell-Mann,
%``Theory Of The Fermi Interaction,''
{\it Phys.\ Rev.\  }{\bf 109}, 193 (1958).
%%CITATION = PHRVA,109,193;%%

\bibitem{Pais}
A.~Pais,{\it  Phys. Rev. }{\bf 86}, 663 (1952).
%%CITATION = PHRVA,86,663;%%


\bibitem{Gell-Mann}
M. Gell-Mann,
{\it Phys. Rev. }{\bf 92}, 833 (1953).
%%CITATION = PHRVA,92,833;%%

\bibitem{Gell-MannPais}
M. Gell-Mann and A. Pais,{\it  Phys. Rev. }{\bf 97}, 1387 (1955).
%%CITATION = PHRVA,97,1387;%%


\bibitem{eightfold}
Y. Ne'eman, {\it Nucl. Phys. }{\bf 26}, 222 (1961);
%%CITATION = NUPHA,26,222;%%
M.~Gell-Mann, {\it Phys. Rev. }{\bf 125}, 1067 (1962).
%%CITATION = PHRVA,125,1067;%%

\bibitem{Cabibbo}
N.~Cabibbo,
%``Unitary Symmetry And Leptonic Decays,''
{\it Phys.\ Rev.\ Lett.\  }{\bf 10}, 531 (1963).
%%CITATION = PRLTA,10,531;%%

\bibitem{quarks}
M.~Gell-Mann,
%``A Schematic Model Of Baryons And Mesons,''
{\it Phys.\ Lett.\  }{\bf 8}, 214 (1964);
%%CITATION = PHLTA,8,214;%%
G.~Zweig, CERN report, unpublished.

\bibitem{CCFT}
J.~H.~Christenson {\it et al},
%, J.~W.~Cronin, V.~L.~Fitch and R.~Turlay,
%``Evidence For The 2 Pi Decay Of The K(2)0 Meson,''
{\it Phys.\ Rev.\ Lett.\  }{\bf 13}, 138 (1964).
%%CITATION = PRLTA,13,138;%%

\bibitem{WuYang}
T.~T.~Wu and C.~N.~Yang,
%``Phenomenological Analysis Of Violation Of CP Invariance In Decay Of K0 And Anti-K0,''
{\it Phys.\ Rev.\ Lett.\  }{\bf 13}, 380 (1964).
%%CITATION = PRLTA,13,380;%%

\bibitem{Wolfenstein}
L.~Wolfenstein,
%``Violation Of CP Invariance And The Possibility Of Very Weak Interactions,''
{\it Phys.\ Rev.\ Lett.\  }{\bf 13}, 562 (1964).
%%CITATION = PRLTA,13,562;%%

\bibitem{Yang-Mills}
C.~N.~Yang and R.~L.~Mills,
%``Conservation Of Isotopic Spin And Isotopic Gauge Invariance,''
{\it Phys.\ Rev.\  }{\bf 96}, 191 (1954).
%%CITATION = PHRVA,96,191;%%

\bibitem{Glashow}
S.~L.~Glashow,
%``Partial Symmetries Of Weak Interactions,''
{\it Nucl.\ Phys.\  }{\bf 22}, 579 (1961).
%%CITATION = NUPHA,22,579;%%

\bibitem{GIM}
S.~L.~Glashow, J.~Iliopoulos and L.~Maiani,
%``Weak Interactions With Lepton - Hadron Symmetry,''
{\it Phys.\ Rev.\  }{\bf D2}, 1285 (1970).
%%CITATION = PHRVA,D2,1285;%%

\bibitem{GaillardLee1}
M.~K.~Gaillard and B.~W.~Lee,
%``Rare Decay Modes Of The K - Mesons In Gauge Theories,''
{\it Phys.\ Rev.\  }{\bf D10}, 897 (1974).
%%CITATION = PHRVA,D10,897;%%

\bibitem{QCD}
H.~Fritzsch, M.~Gell-Mann and H.~Leutwyler,
%``Advantages Of The Color Octet Gluon Picture,''
{\it Phys.\ Lett.\  }{\bf B47}, 365 (1973).
%%CITATION = PHLTA,B47,365;%%

\bibitem{Asympt}
D.~J.~Gross and F.~Wilczek,
%``Ultraviolet Behavior Of Nonabelian Gauge Theories,''
{\it Phys.\ Rev.\ Lett.\  }{\bf 30}, 1343 (1973);
%%CITATION = PRLTA,30,1343;%%
H.~D.~Politzer,
%``Reliable Perturbative Results For Strong Interactions?,''
{\it Phys.\ Rev.\ Lett.\  }{\bf 30}, 1346 (1973).
%%CITATION = PRLTA,30,1346;%%

\bibitem{GaillardLee2}
M.~K.~Gaillard and B.~W.~Lee,
%``Delta I = 1/2 Rule For Nonleptonic Decays In Asymptotically Free Field Theories,''
{\it Phys.\ Rev.\ Lett.\  }{\bf 33}, 108 (1974).
%%CITATION = PRLTA,33,108;%%

\bibitem{AltarelliMaiani}
G.~Altarelli and L.~Maiani,
%``Octet Enhancement Of Nonleptonic Weak Interactions In Asymptotically Free Gauge Theories,''
{\it Phys.\ Lett.\  }{\bf B52}, 351 (1974).
%%CITATION = PHLTA,B52,351;%%

\bibitem{Vainshtein1}
A.~I.~Vainshtein, V.~I.~Zakharov, V.~A.~Novikov and M.~A.~Shifman,
%``Bounds On The Charmed Particles Masses From The K(L) $\to$ 2 
%Mu Decay Width And K(L) - K(S) Mass Difference In An Asymptotically 
%Free Theory Of Strong Interactions,''
{\it Sov.\ J.\ Nucl.\ Phys.\  }{\bf 23}, 540 (1976)
[{\it Yad.\ Fiz.\  }{\bf 23}, 1024 (1976)].
%%CITATION = SJNCA,23,540;%%

\bibitem{Penguin1}
A.~I.~Vainshtein, V.~I.~Zakharov and M.~A.~Shifman,
%``A Possible Mechanism For The Delta T = 1/2 Rule In 
%Nonleptonic Decays Of Strange Particles,''
{\it JETP Lett.\ } {\bf 22}, 55 (1975)
[{\it Pisma Zh.\ Eksp.\ Teor.\ Fiz.\ } {\bf 22}, 123 (1975)];
%%CITATION = JTPLA,22,55;%%
M.~A.~Shifman, A.~I.~Vainshtein and V.~I.~Zakharov,
%``Light Quarks And The Origin Of The Delta I = 1/2 Rule In
% The Nonleptonic Decays Of Strange Particles,''
{\it Nucl.\ Phys.\  }{\bf B120}, 316 (1977).
%%CITATION = NUPHA,B120,316;%%

\bibitem{Penguin2}
M.~A.~Shifman, A.~I.~Vainshtein and V.~I.~Zakharov,
%``Nonleptonic Decays Of K Mesons And Hyperons,''
{\it Sov.\ Phys.\ JETP }{\bf 45}, 670 (1977)
[{\it Zh.\ Eksp.\ Teor.\ Fiz.\  }{\bf 72}, 1275 (1977)].
%%CITATION = SPHJA,45,670;%%

\bibitem{KobayashiMaskawa}
M.~Kobayashi and T.~Maskawa,
%``CP Violation In The Renormalizable Theory Of Weak Interaction,''
{\it Prog.\ Theor.\ Phys.\  }{\bf 49}, 652 (1973).
%%CITATION = PTPKA,49,652;%%

\bibitem{WeinbergCP}
S.~Weinberg,
%``Gauge Theory Of CP Violation,''
{\it Phys.\ Rev.\ Lett.\  }{\bf 37}, 657 (1976).
%%CITATION = PRLTA,37,657;%%

\bibitem{GilmanWise1}
F.~J.~Gilman and M.~B.~Wise,
%``The Delta I = 1/2 Rule And Violation Of CP In The Six Quark Model,''
{\it Phys.\ Lett.\  }{\bf B83}, 83 (1979).
%%CITATION = PHLTA,B83,83;%%

\bibitem{GilmanWise2}
F.~J.~Gilman and M.~B.~Wise,
%``Effective Hamiltonian For Delta S = 1 Weak Nonleptonic Decays In The Six Quark Model,''
{\it Phys.\ Rev.\  }{\bf D20}, 2392 (1979).
%%CITATION = PHRVA,D20,2392;%%

\bibitem{GuberinaPeccei}
B.~Guberina and R.~D.~Peccei,
%``Quantum Chromodynamic Effects And CP Violation In The Kobayashi-Maskawa Model,''
{\it Nucl.\ Phys.\  }{\bf B163}, 289 (1980).
%%CITATION = NUPHA,B163,289;%%

\bibitem{GilmanWise3}
F.~J.~Gilman and M.~B.~Wise,
%``Strong Interaction Corrections To K0 Anti-K0 Mixing In The Six Quark Model,''
{\it Phys.\ Lett.\  }{\bf B93}, 129 (1980).
%%CITATION = PHLTA,B93,129;%%

\bibitem{GilmanWise4}
F.~J.~Gilman and M.~B.~Wise,
%``K0 Anti-K0 Mixing In The Six Quark Model,''
{\it Phys.\ Rev.\  }{\bf D27}, 1128 (1983).
%%CITATION = PHRVA,D27,1128;%%

\bibitem{BijnensWise}
J.~Bijnens and M.~B.~Wise,
%``Electromagnetic Contribution To Epsilon-Prime / Epsilon,''
{\it Phys.\ Lett.\  }{\bf B137}, 245 (1984).
%%CITATION = PHLTA,B137,245;%%

\bibitem{RandallFlynn}
J.~M.~Flynn and L.~Randall,
%``The Electromagnetic Penguin Contribution To Epsilon-Prime / Epsilon For Large Top Quark Mass,''
{\it Phys.\ Lett.\  }{\bf B224}, 221 (1989)
[{\it Phys.\ Lett.\  }{\bf B235}, 412 (1989)].
%%CITATION = PHLTA,B224,221;%%

\bibitem{twoloopold}
G.~Altarelli {\it et al.},
%, G.~Curci, G.~Martinelli and S.~Petrarca,
%``Weak Nonleptonic Decays Beyond Leading Logarithms In QCD,''
{\it Phys.\ Lett.\  }{\bf B99}, 141 (1981);
%%CITATION = PHLTA,B99,141;%%
%G.~Altarelli, G.~Curci, G.~Martinelli and S.~Petrarca,
%``QCD Nonleading Corrections To Weak Decays As An Application Of Regularization By Dimensional Reduction,''
{\it Nucl.\ Phys.\  }{\bf B187}, 461 (1981).
%%CITATION = NUPHA,B187,461;%%

\bibitem{KMW}
J.~Kambor, J.~Missimer and D.~Wyler,
%``K $\to$ 2 pi and K $\to$ 3 pi decays in next-to-leading order
%chiral perturbation theory,''
{\it Phys.\ Lett.\  }{\bf B261}, 496 (1991).
%%CITATION = PHLTA,B261,496;%%

\bibitem{BDP}
J.~Bijnens, P. Dhonte and F. Persson,
%``K $\to$ 3pi decays in chiral perturbation theory,''
.
%%CITATION = ;%%

\bibitem{PDG}
D.~E.~Groom {\it et al.},
%``Review of particle physics,''
{\it Eur.\ Phys.\ J.\  }{\bf C15}, 1 (2000).
%%CITATION = EPHJA,C15,1;%%

\bibitem{epspeps}
A. Alavi-Harati et al. (KTeV), {\it Phys. Rev. Lett.}{\bf 83}, 22 (1999);
%%CITATION = PRLTA,83,22;%%
V.~Fanti~et~al. (NA48),~{\it Phys. Lett. }{\bf B465}, 335 (1999);
%%CITATION = PHLTA,465,335;%%
H. Burkhart et al. (NA31), {\it Phys. Lett. }{\bf B206}, 169 (1988);
%%CITATION = PHLTA,206,169;%%
G.D. Barr et al. (NA31), {\it Phys. Lett. }{\bf B317}, 233 (1993);
%%CITATION = PHLTA,317,233;%%
L.K. Gibbons et al. (E731), {\it Phys. Rev. Lett. }{\bf 70}, 1203 (1993)
%%CITATION = PRLTA,70,1203;%%

\bibitem{Kessler}
R. Kessler, ``Recent KTeV results'',
 .
 
\bibitem{NA4801}
L.  Iconomidou-Fayard, NA48 Collaboration,
%``Results on CP Violation from the NA48 Experiment at CERN,''
%Presented at 20th International Symposium on Lepton and Photon Interactions
%at High Energies (Lepton Photon 01), Rome, Italy, 23-28 Jul 2001, 
.
%%CITATION = ;%%

\bibitem{BPscheme}
J.~Bijnens and J.~Prades,
%``Scheme dependence of weak matrix elements in the 1/N(c) expansion,''
{\it J. High Energy Phys. }{\bf 0001}, 002 (2000)
[\hepph{9911392}].
%%CITATION = ;%%

\bibitem{Buras1}
A.~J.~Buras,
%``Weak Hamiltonian, CP violation and rare decays,''
 Les Houches lectures.
%%CITATION = ;%%

\bibitem{one-loop}
Refs. \citebk{GaillardLee2,AltarelliMaiani,Penguin1,GilmanWise1,GilmanWise2,GuberinaPeccei,BijnensWise,RandallFlynn} and
M.~Lusignoli,
%``Electromagnetic Corrections To The Effective Hamiltonian
%For Strangeness Changing Decays And Epsilon-Prime / Epsilon,''
{\it Nucl.\ Phys.\  }{\bf B325}, 33 (1989);
%%CITATION = NUPHA,B325,33;%%

\bibitem{two-loop}
A.~J.~Buras and P.~H.~Weisz,
%``QCD Nonleading Corrections To Weak Decays In Dimensional
%Regularization And 'T Hooft-Veltman Schemes,''
{\it Nucl.\ Phys.\  }{\bf B333}, 66 (1990);
%%CITATION = NUPHA,B333,66;%%
A.~J.~Buras {\it et al.},
%, M.~Jamin, E.~Lautenbacher and P.~H.~Weisz,
%``Effective Hamiltonians for Delta S = 1 and Delta B = 1 nonleptonic
% decays beyond the leading logarithmic approximation,''
{\it Nucl.\ Phys.\  }{\bf B370}, 69 (1992), 
Addendum-ibid. {\bf B375}, 501 (1992);
%%CITATION = NUPHA,B370,69;%%
A.~J.~Buras, M.~Jamin and M.~E.~Lautenbacher,
%``Two loop anomalous dimension matrix for Delta S = 1 weak
% nonleptonic decays. 2. O(alpha-alpha-s),''
{\it Nucl.\ Phys.\  }{\bf B400}, 75 (1993)
[\hepph{9211321}];
%%CITATION = ;%%
A.~J.~Buras {\it et al.},
%, M.~Jamin, M.~E.~Lautenbacher and P.~H.~Weisz,
%``Two loop anomalous dimension matrix for Delta S = 1 weak
% nonleptonic decays. 1. O(alpha-s**2),''
{\it Nucl.\ Phys.\  }{\bf B400}, 37 (1993)
[\hepph{9211304}];\\
%%CITATION = ;%%
M.~Ciuchini {\it et al.},
%, E.~Franco, G.~Martinelli and L.~Reina,
%``The Delta S = 1 effective Hamiltonian including
% next-to-leading order QCD and QED corrections,''
{\it Nucl.\ Phys.\  }{\bf B415}, 403 (1994)
[\hepph{9304257}].
%%CITATION = ;%%

\bibitem{twoloopreview}
G.~Buchalla, A.~J.~Buras and M.~E.~Lautenbacher,
%``Weak Decays Beyond Leading Logarithms,''
{\it Rev.\ Mod.\ Phys.\  }{\bf 68}, 1125 (1996)
.
%%CITATION = ;%%

\bibitem{BPdIhalf}
J.~Bijnens and J.~Prades,
%``The Delta(I) = 1/2 rule in the chiral limit,''
{\it J. High Energy Phys. }{\bf 9901}, 023 (1999)
[\hepph{9811472}].
%%CITATION = ;%%

\bibitem{BPeps}
J.~Bijnens and J.~Prades,
%``epsilon'(K)/epsilon (K) in the chiral limit,''
{\it J. High Energy Phys. }{\bf 0006}, 035 (2000)
.
%%CITATION = ;%%

\bibitem{Sachrajda}
C.~T.~Sachrajda,
%``Phenomenology from lattice QCD,''
%Presented at 20th International Symposium on Lepton and Photon
%Interactions at High Energies (Lepton Photon 01), Rome, Italy, 23-28 Jul 2001,
;
%%CITATION = ;%%
G.~Martinelli,
%``Kaon physics from lattice QCD,''
.
%%CITATION = ;%%

\bibitem{SVZ}
M.~A.~Shifman, A.~I.~Vainshtein and V.~I.~Zakharov,
%``QCD And Resonance Physics. Sum Rules,''
{\it Nucl.\ Phys.\  }{\bf B147}, 385 (1979);
%%CITATION = NUPHA,B147,385;%%
%M.~A.~Shifman, A.~I.~Vainshtein and V.~I.~Zakharov,
%``QCD And Resonance Physics: Applications,''
{\it Nucl.\ Phys.\  }{\bf B147}, 448 (1979).
%%CITATION = NUPHA,B147,448;%%

\bibitem{BBG}
W.~A.~Bardeen, A.~J.~Buras and J.-M.~G\'erard,
%``A Consistent Analysis Of The Delta I = 1/2 Rule For K Decays,''
{\it Phys.\ Lett.\  }{\bf B192}, 138 (1987);
%%CITATION = PHLTA,B192,138;%%
%W.~A.~Bardeen, A.~J.~Buras and J.-M.~G\'erard,
%``The K $\to$ Pi Pi Decays In The Large N Limit: Quark Evolution,''
{\it Nucl.\ Phys.\  }{\bf B293}, 787 (1987).
%%CITATION = NUPHA,B293,787;%%

\bibitem{Hambye}
T.~Hambye {\it et al.},
%, G.~O.~K\"ohler, E.A.~Paschos, P.H.~Soldan and W.A.~Bardeen,
%``1/N(c) corrections to the hadronic matrix elements of Q(6) and Q(8) in  K --> pi pi decays,''
{\it Phys.\ Rev.\  }{\bf D58}, 014017 (1998)
[\hepph{9802300}];
%%CITATION = ;%%
T.~Hambye, G.~O.~K\"ohler and P.~H.~Soldan,
%``New analysis of the Delta(I) = 1/2 rule in kaon decays and the  B-hat(K) parameter,''
{\it Eur.\ Phys.\ J.\  }{\bf C10}, 271 (1999)
[\hepph{9902334}];
%%CITATION = ;%%
T.~Hambye {\it et al.},
%, G.~O.~K\"ohler, E.~A.~Paschos and P.~H.~Soldan,
%``Analysis of epsilon'/epsilon in the 1/N(c) expansion,''
{\it Nucl.\ Phys.\  }{\bf B564}, 391 (2000)
[\hepph{9906434}].
%%CITATION = ;%%

\bibitem{ENJL}
J.~Bijnens, C.~Bruno and E.~de Rafael,
%``Nambu-Jona-Lasinio like models and the low-energy effective action of QCD,''
{\it Nucl.\ Phys.\  }{\bf B390}, 501 (1993)
[\hepph{9206236}];
%%CITATION = ;%%
J.~Bijnens,
%``Chiral Lagrangians and Nambu-Jona-Lasinio - like models,''
{\it Phys.\ Rept.\  }{\bf 265}, 369 (1996)
[\hepph{9502335}] and references therein.
%%CITATION = ;%%

\bibitem{BPBK}
J.~Bijnens and J.~Prades,
%``B(K) and explicit chiral symmetry breaking,''
{\it Phys.\ Lett.\  }{\bf B342}, 331 (1995)
[\hepph{9409255}];
%%CITATION = ;%%
%J.~Bijnens and J.~Prades,
%``The B(K) parameter in the 1/N(c) expansion,''
{\it Nucl.\ Phys.\  }{\bf B444}, 523 (1995)
[\hepph{9502363}].
%%CITATION = ;%%

\bibitem{kptokpp}
J.~Bijnens, E.~Pallante and J.~Prades,
%``Obtaining K --> pi pi from off-shell K --> pi amplitudes,''
{\it Nucl.\ Phys.\  }{\bf B521}, 305 (1998)
[\hepph{9801326}].
%%CITATION = ;%%

\bibitem{LMD}
M.~Knecht, S.~Peris and E.~de Rafael,
%``Matrix elements of electroweak penguin operators in the 1/N(c)  expansion,''
{\it Phys.\ Lett.\  }{\bf B457}, 227 (1999)
[\hepph{9812471}];
%%CITATION = ;%%
S.~Peris and E.~de Rafael,
%``K0 anti-K0 mixing in the 1/N(c) expansion,''
{\it Phys.\ Lett.\  }{\bf B490}, 213 (2000)
[hep-ph/\.
%%CITATION = ;%%
E.~de Rafael,
%``Large N(c) methods in kaon physics,''
arXiv:.
%%CITATION = ;%%

\bibitem{PichdeRafael}
A.~Pich and E.~de Rafael,
%``Four quark operators and nonleptonic weak transitions,''
{\it Nucl.\ Phys.\  }{\bf B358}, 311 (1991);
%%CITATION = NUPHA,B358,311;%%
V.~Antonelli {\it et al.},
%, S.~Bertolini, J.~O.~Eeg, M.~Fabbrichesi and E.~I.~Lashin,
%``The $\Delta S = 1$ Weak Chiral Lagrangian as the Effective Theory of the Chiral Quark Model,''
{\it Nucl.\ Phys.\  }{\bf B469}, 143 (1996)
[\hepph{9511255}];
%%CITATION = ;%%
V.~Antonelli {\it et al.},
%, S.~Bertolini, M.~Fabbrichesi and E.~I.~Lashin,
%``The $\Delta I = 1/2$ Selection Rule,''
{\it Nucl.\ Phys.\  }{\bf B469}, 181 (1996)
[\hepph{9511341}];
%%CITATION = ;%%
S.~Bertolini, J.~O.~Eeg and M.~Fabbrichesi,
%``A New Estimate of $\varepsilon '/\varepsilon$,''
{\it Nucl.\ Phys.\  }{\bf B476}, 225 (1996)
[\hepph{9512356}];
%%CITATION = ;%%
S.~Bertolini {\it et al.},
%, J.~O.~Eeg, M.~Fabbrichesi and E.~I.~Lashin,
%``epsilon'/epsilon at O(p**4) in the chiral expansion,''
{\it Nucl.\ Phys.\  }{\bf B514}, 93 (1998)
[\hepph{9706260}];
%%CITATION = ;%%
S.~Bertolini, J.~O.~Eeg and M.~Fabbrichesi,
%``An updated analysis of epsilon'/epsilon in the standard model with  hadronic matrix elements from the chiral quark model,''
{\it Phys.\ Rev.\ }{\bf D63}, 056009 (2001)
;
%%CITATION = ;%%
M.~Franz, H.~C.~Kim and K.~Goeke,
%``Effective weak chiral Lagrangian to O(p**4) in the chiral quark model,''
{\it Nucl.\ Phys.\  }{\bf B562}, 213 (1999)
.
%%CITATION = ;%%

\bibitem{EPetal}
G.~Ecker {\it et al.},
% G.~Muller, H.~Neufeld and A.~Pich,
%``pi0 eta mixing and CP violation,''
{\it Phys.\ Lett.\  }{\bf B477}, 88 (2000)
.
%%CITATION = ;%%

\bibitem{PP}
E.~Pallante and A.~Pich,
%``Strong enhancement of epsilon'/epsilon through final state  interactions,''
{\it Phys.\ Rev.\ Lett.\  }{\bf 84}, 2568 (2000)
;
%%CITATION = ;%%
%E.~Pallante and A.~Pich,
%``Final state interactions in kaon decays,''
{\it Nucl.\ Phys.\ }{\bf B592}, 294 (2001)
.
%%CITATION = ;%%

\bibitem{Dasetal}
T.~Das, G.~S.~Guralnik, V.~S.~Mathur, F.~E.~Low and J.~E.~Young,
%``Electromagnetic Mass Difference Of Pions,''
{\it Phys.\ Rev.\ Lett.\  }{\bf 18}, 759 (1967).
%%CITATION = PRLTA,18,759;%%

\bibitem{DG1}
J.~F.~Donoghue and E.~Golowich,
%``Anatomy Of A Weak Matrix Element,''
{\it Phys.\ Lett.\  }{\bf B315}, 406 (1993)
;
%%CITATION = ;%%
%J.~F.~Donoghue and E.~Golowich,
%``Dispersive calculation of B7(3/2) and B8(3/2) in the chiral limit,''
{\it Phys.\ Lett.\  }{\bf B478}, 172 (2000)
;\\
%%CITATION = ;%%
V.~Cirigliano and E.~Golowich,
%``Analysis of O(p**2) corrections to ,''
{\it Phys.\ Lett.\  }{\bf B475}, 351 (2000)
;\\
%%CITATION = ;%%
V.~Cirigliano, J.~F.~Donoghue and E.~Golowich,
%``Dimension-eight operators in the weak OPE,''
{\it J. High Energy Phys. }{\bf 0010}, 048 (2000)
;\\
%%CITATION = ;%%
V.~Cirigliano and E.~Golowich,
%``Comment on 'Analysis of O(p**2) corrections to ',''
{\it Phys.\ Rev.\  }{\bf D65}, 054014 (2002)
.
%%CITATION = ;%%

\bibitem{KnechtQ7Q8}
M.~Knecht, S.~Peris and E.~de Rafael,
%``A critical reassessment of Q(7) and Q(8) matrix elements,''
{\it Phys.\ Lett.\  }{\bf B508}, 117 (2001)
 and private communication.
%%CITATION = ;%%

\bibitem{NAR01} S.~Narison,
%``New QCD-estimate of the kaon penguin matrix elements and 
% epsilon'/epsilon,''
 {\it Nucl.\ Phys.\  }{\bf B593}, 3 (2001) ;
%%CITATION = ;%%
%S.~Narison,
%``V-A hadronic tau decays: A QCD laboratory,''
 {\it Nucl.\ Phys.\ Proc.\ Suppl.\  }{\bf 96}, 364 (2001)
 .
%%CITATION = ;%%

\bibitem{DG2}
V.~Cirigliano, J.~F.~Donoghue, E.~Golowich and K.~Maltman,
%``Determination of <(pi pi)(I=2)$|$Q(7,8)$|$K0> in the chiral limit,''
{\it Phys.\ Lett.\  }{\bf B522}, 245 (2001)
.
%%CITATION = ;%%

\bibitem{BPQ7Q8}\label{BPQ7Q8b}
J.~Bijnens, E.~Gamiz and J.~Prades,
%``Matching the electroweak penguins Q(7), Q(8) and spectral correlators,''
{\it J. High Energy Phys. }{\bf 0110}, 009 (2001)
.
%%CITATION = ;%%

\bibitem{Q7Q8lattice}  A. Donini {\it et al.},
%, V. Gim\'enez, L. Giusti, and G. Martinelli,
%''Renormalization Group Invariant Matrix Elements of DS = 2 
%and DI = 3/2Four-Fermion Operators without Quark Mass,''
 {\it Phys. Lett.  }{\bf B470}, 233 (1999)
%%CITATION = ;%%

\end{thebibliography}

\end{document}

