

%%%  useful macros for econf proceedings articles:

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%  basic formatting macros:
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%  single-line equations:

\def\beq{\begin{equation}}
\def\eeq#1{\label{#1}\end{equation}}
\def\eeqn{\end{equation}}

%%  multiple-line equations  (use \CR as the carriage return):

\newenvironment{Eqnarray}%
   {\arraycolsep 0.14em\begin{eqnarray}}{\end{eqnarray}}
\def\beqa{\begin{Eqnarray}}
\def\eeqa#1{\label{#1}\end{Eqnarray}}
\def\eeqan{\end{Eqnarray}}
\def\CR{\nonumber \\ }

%%  reference to an equation number:

\def\leqn#1{(\ref{#1})}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%  bibliographic entries   (use this or the SPIRES LaTeX output)

%%   Journal or bibliographic formatting macros are obsolete!
%%   The SLAC/SPIRES database supplies properly formatted citations;
%%      click on  LaTeX(US) or LaTeX(EU)   

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%  sizing and bars

\def\st{\scriptstyle}
\def\sst{\scriptscriptstyle}

\def\overbar#1{\overline{#1}}
\let\littlebar=\bar
\let\bar=\overbar

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%  text-mode macros:

\def\etal{{\it et al.}}
\def\ie{{\it i.e.}}
\def\eg{{\it e.g.}}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%  expectation values:

\def\VEV#1{\left\langle{ #1} \right\rangle}
\def\bra#1{\left\langle{ #1} \right|}
\def\ket#1{\left| {#1} \right\rangle}
\def\vev#1{\langle #1 \rangle}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%  caligraphic letters (for matrix elements, luminosity, etc.)

\def\D{{\cal D}}
\def\L{{\cal L}}
\def\M{{\cal M}}
\def\O{{\cal O}}
\def\W{{\cal W}}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%  matrix operations and fractions:

\def\One{{\bf 1}}
\def\hc{{\mbox{\rm h.c.}}}
\def\tr{{\mbox{\rm tr}}}
\def\half{\frac{1}{2}}
\def\thalf{\frac{3}{2}}
\def\third{\frac{1}{3}}
\def\tthird{\frac{2}{3}}

\def\Dslash{\not{\hbox{\kern-4pt $D$}}}
\def\dslash{\not{\hbox{\kern-2pt $\del$}}}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%  high-energy physics terminology:

\def\Pl{{\mbox{\scriptsize Pl}}}
\def\eff{{\mbox{\scriptsize eff}}}
\def\CM{{\mbox{\scriptsize CM}}}
\def\BR{\mbox{\rm BR}}
\def\ee{e^+e^-}
\def\sstw{\sin^2\theta_w}
\def\cstw{\cos^2\theta_w}
\def\mz{m_Z}
\def\gz{\Gamma_Z}
\def\mw{m_W}
\def\mt{m_t}
\def\gt{\Gamma_t}
\def\mh{m_h}
\def\gmu{G_\mu}
\def\GF{G_F}
\def\alphas{\alpha_s}
\def\msb{{\bar{\ssstyle M \kern -1pt S}}}
\def\lmsb{\Lambda_{\msb}}
\def\ELER{e^-_Le^+_R}
\def\EREL{e^-_Re^+_L}
\def\ELEL{e^-_Le^+_L}
\def\ERER{e^-_Re^+_R}
\def\eps{\epsilon}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%  supersymmetry:

\def\ch#1{\widetilde\chi^+_{#1}}
\def\chm#1{\widetilde\chi^-_{#1}}
\def\neu#1{\widetilde\chi^0_{#1}}
\def\s#1{\widetilde{#1}}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5

%% miscellaneous

\def\lsim{\mathrel{\raise.3ex\hbox{$<$\kern-.75em\lower1ex\hbox{$\sim$}}}}
\def\gsim{\mathrel{\raise.3ex\hbox{$>$\kern-.75em\lower1ex\hbox{$\sim$}}}}
%  sample eprint article in LaTeX           --- M. Peskin, 9/7/00
%  enhancements added by H.E. Haber, 11/16/00

\documentclass[12pt]{article}
\usepackage{epsfig}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% basic data for the eprint:
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\textwidth=6.0in  \textheight=8.25in

%%  Adjust these for your printer:
\leftmargin=-0.3in   \topmargin=-0.20in
\hoffset=-.25in
%% preprint number data:
%% Change the preprint number below to your institutional preprint number
\newcommand\pubnumber{}
\newcommand\pubdate{}
%% If you will submit to hep-ex, change hep-ph to hep-ex below
%% Do not ; this will be automatically converted
%% to the Los Alamos hep number during the submission process
\newcommand\hepnumber

%%  address and funding acknowledgement data:
\def\infn{INFN and SISSA\\
Via Beirut 4, I-34013 Trieste, Italy}
\def\support{\footnote{...}}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%   document style macros
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\def\Title#1{\begin{center} {\Large\bf #1 } \end{center}}
\def\Author#1{\begin{center}{ \sc #1} \end{center}}
\def\Address#1{\begin{center}{ \it #1} \end{center}}
\def\andauth{\begin{center}{and} \end{center}}
\def\submit#1{\begin{center}Submitted to {\sl #1} \end{center}}
\newcommand\pubblock{\rightline{\begin{tabular}{l} \pubnumber\\
         \pubdate  \end{tabular}}}
\newenvironment{Abstract}{\begin{quotation}  }{\end{quotation}}
\newenvironment{Presented}{\begin{quotation} \begin{center} 
             Presented at the \end{center}\bigskip 
      \begin{center}\begin{large}}{\end{large}\end{center} \end{quotation}}
\def\Acknowledgements{\bigskip  \bigskip \begin{center} \begin{large}
             \bf Acknowledgments \end{large}\end{center}}

\makeatletter
\def\section{\@startsection{section}{0}{\z@}{5.5ex plus .5ex minus
 1.5ex}{2.3ex plus .2ex}{\large\bf}}
\def\subsection{\@startsection{subsection}{1}{\z@}{3.5ex plus .5ex minus
 1.5ex}{1.3ex plus .2ex}{\normalsize\bf}}
\def\subsubsection{\@startsection{subsubsection}{2}{\z@}{-3.5ex plus
-1ex minus  -.2ex}{2.3ex plus .2ex}{\normalsize\sl}}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% small size table and figure captions %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\renewcommand{\@makecaption}[2]{%
   \vskip 10pt
   \setbox\@tempboxa\hbox{\small #1: #2}
   \ifdim \wd\@tempboxa >\hsize     % IF longer than one line:
       \small #1: #2\par          %   THEN set as ordinary paragraph.
     \else                        %   ELSE  center.
       \hbox to\hsize{\hfil\box\@tempboxa\hfil}
   \fi}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% macros to collapse citation numbers to ranges %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% \citenum emits the plain citation number without ornament
% \citea puts it's argument into the ornamentation for citations
% thus \cite{foo} is equivalent to \citea{\citenum{foo}}
 \def\citenum#1{{\def\@cite##1##2{##1}\cite{#1}}}
%\def\citea#1{\@cite{#1}{}}
 
% Collapse citation numbers to ranges.  Non-numeric and undefined labels
% are handled.  No sorting is done.  E.g., 1,3,2,3,4,5,foo,1,2,3,?,4,5
% gives 1,3,2-5,foo,1-3,?,4,5
\newcount\@tempcntc
\def\@citex[#1]#2{\if@filesw\immediate\write\@auxout{\string\citation{#2}}\fi
  \@tempcnta\z@\@tempcntb\m@ne\def\@citea{}\@cite{\@for\@citeb:=#2\do
    {\@ifundefined
       {b@\@citeb}{\@citeo\@tempcntb\m@ne\@citea\def\@citea{,}{\bf ?}\@warning
       {Citation `\@citeb' on page \thepage \space undefined}}%
    {\setbox\z@\hbox{\global\@tempcntc0\csname b@\@citeb\endcsname\relax}%
     \ifnum\@tempcntc=\z@ \@citeo\@tempcntb\m@ne
       \@citea\def\@citea{,}\hbox{\csname b@\@citeb\endcsname}%
     \else
      \advance\@tempcntb\@ne
      \ifnum\@tempcntb=\@tempcntc
      \else\advance\@tempcntb\m@ne\@citeo
      \@tempcnta\@tempcntc\@tempcntb\@tempcntc\fi\fi}}\@citeo}{#1}}
\def\@citeo{\ifnum\@tempcnta>\@tempcntb\else\@citea\def\@citea{,}%
  \ifnum\@tempcnta=\@tempcntb\the\@tempcnta\else
  {\advance\@tempcnta\@ne\ifnum\@tempcnta=\@tempcntb \else\def\@citea{--}\fi
    \advance\@tempcnta\m@ne\the\@tempcnta\@citea\the\@tempcntb}\fi\fi}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\makeatother


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% corrects bad spacing in LaTeX multiple line equations
% replace eqnarray with Eanarray (or use abbreviations
% in econfmacros2.tex
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%\newenvironment{Eqnarray}%
%   {\arraycolsep 0.14em\begin{eqnarray}}{\end{eqnarray}}
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%  personal abbreviations and macros
%    the following package contains macros used in this document:
\input econfmacros2.tex
%
\def\vev#1{\left\langle #1\right\rangle}
\def\Im{\mathop{\mbox{Im}}}
\def\Re{\mathop{\mbox{Re}}}
\def\Tr{\mathop{\mbox{Tr}}\,}

\def\qq{$\vev{\bar qq}$~}
\def\GG{$\vev{\alpha_s GG/\pi}$~}
\def\eoe{$\varepsilon'/\varepsilon$~}
\def\ratio{\eoe}
\def\eps{$\varepsilon$~}
%
\def\be{\begin{equation}}
\def\ee{\end{equation}}
\def\bea{\begin{eqnarray}}
\def\eea{\end{eqnarray}}
% A useful Journal macro
\def\Journal#1#2#3#4{{\sl #1} {\bf #2} (#4) #3}
% Some useful journal names
\def\NCA{Nuovo Cimento}
\def\NPB{Nucl. Phys. B}
\def\PLB{Phys. Lett. B}
\def\PRD{Phys. Rev. D}
\def\PRL{Phys. Rev. Lett.}
\def\RMP{Rev. Mod. Phys.}
\def\ZPC{Z. Phys. C}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\begin{document}
\begin{titlepage}
\pubblock

\vfill
\Title{Theory of \eoe}
\vfill
\Author{Stefano Bertolini}
\Address{\infn}
\vfill
\begin{Abstract}
I shortly review the present status of the
theoretical estimates of \eoe.
I consider a few aspects of the theoretical calculations
which may be relevant in understanding the present experimental
results. I discuss the role of higher order chiral corrections 
and in general of non-factorizable contributions for the
explanation of the $\Delta I = 1/2$ selection rule in kaon decays
and \eoe. Lacking reliable lattice calculations, the $1/N$ expansion
and phenomenological approaches may help in understanding 
correlations among theoretical effects and the experimental data.
The same dynamics which underlies the CP conserving selection rule
appears to drive \eoe in the range of the recent experimental measurements.
\end{Abstract}
\vfill
\begin{Presented}
5th International Symposium on Radiative Corrections \\ 
(RADCOR--2000) \\[4pt]
Carmel CA, USA, 11--15 September, 2000
\end{Presented}
\vfill
\end{titlepage}
\def\thefootnote{\fnsymbol{footnote}}
\setcounter{footnote}{0}
%

\section{Introduction}
The results announced during the last year by the KTeV~\cite{KTeV} 
and NA48~\cite{NA48} collaborations
have marked a great experimental achievement,
establishing 35 years after the discovery of CP violation
in the neutral kaon system~\cite{Christenson}
the existence of a much smaller violation acting directly in the
decays. 

While the Standard Model (SM) of strong and electroweak interactions
provides an economical and elegant understanding
of indirect~($\varepsilon$) and direct~($\varepsilon'$) 
CP violation in term of a single phase,
the detailed calculation of the size of these effects 
implies mastering strong interactions at a scale 
where perturbative methods break down. In addition, 
CP violation in $K\to\pi\pi$ decays
is the result of a destructive interference between
two sets of contributions,
which may inflate up to an order of magnitude the uncertainties
on the hadronic matrix elements of the effective 
four-quark operators.
All that makes predicting \eoe a complex and subtle theoretical
challenge~\cite{review}.

\begin{figure}[b!]
\begin{center}
\epsfig{file=fig_pre.eps,height=7.0cm}
\caption[0]{
The 1-$\sigma$ results of the
NA31 and E731 Collaborations (early 90's)
are shown by the gray horizontal bands.
The old M\"unchen, Roma and Trieste theoretical predictions for \eoe are
depicted by the vertical bars with their central values.
For comparison, the VSA estimate is shown using two renormalization schemes.
\label{fig:pre}}
\end{center}
\end{figure}

The status of the theoretical 
predictions and experimental data available before 
the KTeV announcement in February 1999
is summarized in Fig.~\ref{fig:pre}. 

The gray horizontal bands show
the one-sigma average of the old (early 90's)
NA31 (CERN) and E731 (Fermilab) results. 
The vertical lines show the ranges of the
published theoretical $predictions$ (before February 1999), 
identified with the cities where most of the group members reside.
The range of the naive Vacuum Saturation Approximation (VSA) is
shown for comparison. 

The experimental and theoretical scenarios have changed substantially
after the first KTeV data and the subsequent NA48 results.   
Fig.~\ref{fig:post} shows the present experimental world average
for \eoe compared with the revised or new theoretical calculations
that appeared during the last year. 

Notwithstanding the complexity of the problem,
all theoretical calculations 
show a remarkable overall agreement,
most of them pointing to a non-vanishing positive effect in the SM
(which is by itself far from trivial).

On the other hand, if we focus our attention on
the central values,
many of the predictions
prefer the $10^{-4}$ regime, whereas only a few of them stand 
above $10^{-3}$.
Is this just a ``noise'' in the theoretical calculations?

The answer is no.
Without entering the details of the various estimates, 
it is possible to explain most of
the abovementioned difference in terms
of a single effect: the different size of the 
hadronic matrix element of the gluonic penguin $Q_6$
obtained in the various approaches. 

\begin{figure}[t]
\begin{center}
\epsfig{file=expvsth.eps,height=7.0cm}
\caption[0]{
The latest theoretical calculations of \eoe are compared with
the combined 1-$\sigma$ average of the
NA31, E731, KTeV and NA48 results (\eoe = $19.2\pm 4.6\times 10^{-4}$),
depicted by the gray horizontal band (the error is inflated
according to the Particle Data Group procedure when averaging over
data with substantially different central values).  
\label{fig:post}}
\end{center}
\end{figure}

While some of the calculations, as the early M\"unchen and Rome predictions,
assume for $\vev{\pi\pi|Q_6|K}$ values in the neighboroud of the 
leading $1/N$ result (naive factorization), other approaches,
first of which the Trieste and Dortmund calculations and more recently the 
Lund and Valencia analyses,
find a substantial enhancement 
of this matrix element with respect to the simple factorization result.
The bulk of such an effect
is actually a global enhancement of the
$I=0$ components of the $K\to\pi\pi$ amplitudes, 
which affects $both$ current-current $and$ penguin operators, and
it can be at least partially understood in terms of 
chiral dynamics (final-state interactions).


\section{Final State Interactions}
As a matter of fact, one should in general expect an enhancement of
\eoe with respect to the naive VSA due to final-state interactions (FSI). 
As Fermi first argued~\cite{Fermi}, in potential scattering 
the isospin $I=0$ two-body~states feel an 
attractive interaction, of a sign opposite to that of the $I=2$
components thus affecting the size of the corresponding amplitudes. 
This feature is at the root of the enhancement of the
$I=0$ amplitude over the $I=2$ one and of the corresponding enhancement
of \eoe beyond factorization.

The question is how to make of a qualitative statement 
a quantitative one.
A dispersive analysis of the $K\to\pi\pi$ amplitudes has been recently
presented in ref.~\cite{dispersive}. 
The Omn\`es-Mushkelishvili representation~\cite{Omnes}
\be 
M(s+i\epsilon) = P(s)\ {\rm exp}\left(\frac{1}{\pi}
\int_{4 m_\pi^2}^\infty\frac{\delta(s')}{s'-s-i\epsilon} ds'\right)
\label{omnes}
\ee
is used in order to resum FSI effects from the knowledge of the
$\pi\pi$ rescattering phase $\delta(s)$ in the elastic 
regime ($s < 1$ GeV$^2$).
$P(s)$ is a polinomial function of $s$ which is related
to the factorized amplitude.
A solution of the above dispersive relation for the $A_{0,2}$
amplitudes can be written as
\be
A_I(s) = A'\ (s-m_\pi^2)\ {\cal R}_I(s)\ e^{i\delta_I(s)}\ ,
\label{omnes-sol}
\ee
where $A'$ is the derivative of the amplitude at the subtraction
point $s=m_\pi^2$. The coefficient $\cal R$ represents the 
rescaling effect related to the FSI. By replacing $A'_I$ with
the value given by LO chiral perturbation theory,
Pich and Pallante found ${\cal R}(m_k^2)_{0,2} \simeq  1.4,\ 0.9$
thus confirming via the resummation of the leading chiral logs related
to FSI the enhancement
of the $I=0$ amplitudes, together with a mild depletion of
the $I=2$ components. 

The numerical significance of these results has been 
questioned~\cite{burasetal} on the basis that
the precise size of the effect 
depends on boundary conditions of the factorized amplitude which are
not unambiguously known, due to 
higher order chiral corrections.
As a matter of fact,
while the choice of a low subtraction scale minimizes
the effect of momentum dependent chiral corrections 
the result of ref.~\cite{dispersive} cannot account
for polinomial corrections due to 
contact interaction terms whose size is unknown (and renormalization-scheme
dependent).

The analysis of ref.~\cite{dispersive} shows non-perturbatively
the presence of a potentially large departure from factorization 
which affects the $I=0$ $K\to\pi\pi$ matrix elements.
Nevertheless, the question whether the FSI rescaling of the factorized 
isospin amplitudes 
leads by itself to a satisfactory calculation of \ratio remains open.


\section{CP conserving versus CP violating amplitudes}

Given the possibility that common systematic uncertainties may 
affect the calculation of \eoe and the $\Delta I = 1/2$ rule
(see for instance the present difficulties in calculating on the lattice
the ``penguin contractions'' for CP violating as well as for
CP conserving amplitudes~\cite{Martinelli})
a convincing calculation of \eoe must involve at the same time
a reliable explanation of the $\Delta I = 1/2$ selection rule,
which is still missing. FSI effects alone are {\em not} enough to account
for the large ratio of the $I=0,2$ CP conserving amplitudes.
Other sources of large non-factorizable corrections are needed, which
may affect the determination of \eoe as well. 

The $\Delta I = 1/2$ selection rule in $K\to\pi\pi$ decays is known since 
45 years~\cite{Pais-Gell-Mann} and it states the experimental
evidence that kaons
are 400 times more likely to decay in the $I=0$ two-pion state
than in the $I=2$ component. This rule is not justified by any
general symmetry consideration and, although it is common understanding 
that its explanation must be rooted in the dynamics of strong interactions,
there is no up to date derivation of this effect from first principle QCD. 

The ratio of $I=2$ over $I=0$  amplitudes appears 
directly in the definition of \eoe:
\bea
\frac{\varepsilon'}{\varepsilon} &=& \frac{1}{\sqrt{2}} \left\{
\frac{\langle ( \pi \pi )_{I=2} | {\cal H}_W | {K_L} \rangle}
{\langle ( \pi \pi )_{I=0} | {\cal H}_W | {K_L} \rangle} \right.
 - \left.
\frac{\langle ( \pi \pi )_{I=2} | {\cal H}_W | {K_S} \rangle}
{\langle ( \pi \pi )_{I=0} | {\cal H}_W | {K_S} \rangle} \right\}\ .
\label{eedi}
\eea
As a consequence, a self-consistent calculation of \eoe must also address 
the determination of the CP conserving amplitudes. 

The way we approach the calculation of the hadronic $K\to\pi\pi$ transitions
in gauge theories is provided by the Operator Product Expansion 
which allows us to write
the relevant amplitudes in terms of the hadronic matrix elements of
effective $\Delta S = 1$ four quark operators (at a scale $\mu$)
and of the corresponding Wilson coefficients,
which encode the information about the dynamical degrees
of freedom heavier than the chosen renormalization scale:   
\be
{\cal H}_{\Delta S = 1} =
\frac{G_F}{\sqrt{2}} V_{ud}\,V^*_{us} \sum_i \Bigl[{z_i}(\mu) 
+ {\tau}\ {y_i}(\mu) \Bigr] {Q_i} (\mu) \ .
\label{Leff}
\ee
The entries $V_{ij}$ of the $3\times 3$ Cabibbo-Kobayashi-Maskawa
matrix describe the flavour mixing in the SM
and ${\tau} = - {V_{td}V_{ts}^{*}}/V_{ud}V_{us}^{*}$.
For $\mu < m_c$ ($q=u,d,s$), the relevant quark operators are:
\be
\begin{array}{lcl}
\left.
\begin{array}{lcl}
{Q_{1}} & = & \left( \overline{s}_{\alpha} u_{\beta}  \right)_{\rm V-A}
            \left( \overline{u}_{\beta}  d_{\alpha} \right)_{\rm V-A}
\\[1ex]
{Q_{2}} & = & \left( \overline{s} u \right)_{\rm V-A}
            \left( \overline{u} d \right)_{\rm V-A}
\end{array}
\right\} &&\hspace{-0em} \mbox{Current-Current}
\\[4ex]
\left.
\begin{array}{lcl}
{Q_{3,5}} & = & \left( \overline{s} d \right)_{\rm V-A}
   \sum_{q} \left( \overline{q} q \right)_{\rm V\mp A}
\\[1ex]
{Q_{4,6}} & = & \left( \overline{s}_{\alpha} d_{\beta}  \right)_{\rm V-A}
   \sum_{q} ( \overline{q}_{\beta}  q_{\alpha} )_{\rm V\mp A}
\end{array}
\right\} &&\hspace{-0em} \mbox{Gluon ``penguins''}
\\[4ex]
\left.
\begin{array}{lcl}
{Q_{7,9}} & = & \frac{3}{2} \left( \overline{s} d \right)_{\rm V-A}
         \sum_{q} \hat{e}_q \left( \overline{q} q \right)_{\rm V\pm A}
\\[1ex]
{Q_{8,10}} & = & \frac{3}{2} \left( \overline{s}_{\alpha} 
                                                 d_{\beta} \right)_{\rm V-A}
     \sum_{q} \hat{e}_q ( \overline{q}_{\beta}  q_{\alpha})_{\rm V\pm A}
\end{array}
\right\} &&\hspace{-0em} \mbox{Electroweak ``penguins''}
\end{array} 
\label{quarkeff}
\ee
Current-current operators are induced by tree-level W-exchange whereas
the so-called penguin (and ``box'') diagrams are generated via an 
electroweak loop.
Only the latter ``feel'' all three quark families via the virtual quark
exchange and are therefore sensitive to the weak CP phase.
Current-current operators control instead the CP conserving
transitions. This fact suggests already that the connection 
between \eoe and the
$\Delta I = 1/2$ rule is by no means a straightforward one.

Using the effective $\Delta S=1$ quark Hamiltonian we can write \eoe as
\be
\frac{{\varepsilon'}}{\varepsilon} =  
e^{i \phi} \frac{G_F \omega}{2|\epsilon|\Re{A_0}} \:
{\mbox{Im}\, \lambda_t} \: \:
 \left[ {\Pi_0} - \frac{1}{\omega} \: {\Pi_2} \right]
\label{main}
\ee
where
\be
\begin{array}{lcl}
 {\Pi_0} & = & \frac{1}{{\cos\delta_0}} 
\sum_i {y_i} \, 
\Re\langle  Q_i  \rangle _0\ (1 - {\Omega_{\eta +\eta'}}) 
\\[1ex]
 {\Pi_2} & = & \frac{1}{{\cos\delta_2}} \sum_i {y_i} \, 
\Re\langle Q_i \rangle_2 \quad ,
\end{array} 
\label{PI02}
\ee
and $\langle Q_i \rangle \equiv \langle \pi\pi | Q_i | K \rangle$.
The rescattering phases $\delta_{0,2}$ can be extracted from elastic $\pi-\pi$
scattering data~\cite{FSIphases} and are such that $\cos\delta_0 \simeq 0.8$
and $\cos\delta_2 \simeq 1$. Given that the phase of $\varepsilon$,
$\theta_\varepsilon$, is approximately $\pi/4$, 
as well as $\delta_0-\delta_2$, 
$\phi = \frac{\pi}{2} + {\delta_2} - {\delta_0} - \theta_\varepsilon$ 
turns out to be consistent with zero.

Two key ingredients appear in eq. \ref{main}:
\begin{itemize}
\item[1.]
The isospin breaking $\pi^0-\eta-\eta'$ mixing, parametrized by
$\Omega_{\eta+\eta'}$, which is estimated to give
a $positive$ correction to the $A_2$ amplitude of about 15-35\%.
The complete inclusion of NLO chiral corrections to the $\pi^0-\eta-\eta'$ 
mixing~\cite{GardnerValencia,etapiNLO} and of additional
isospin breaking effects 
($\Delta I =5/2$~\cite{omegaNLOstr,omegaNLOmodel,omegaNLOem}) 
in the extraction of the isospin amplitudes may sizeably affect
the determination of \ratio. Although a partial cancellation
of the new terms in \ratio
reduces their numerical impact, we must await for further analyses
in order to confidently assess their relevance.

\item[2.]
The combination of CKM elements  
{$\Im \lambda_t \equiv \Im (V_{ts}^*V_{td})$}, which
affects directly the size of \eoe and
the range of the uncertainty.
The determination of $\Im \lambda_t$ depends on B-physics
constraints and on $\varepsilon$~\cite{Ciuchinietal}.
In turn, the fit of $\varepsilon$ depends on the 
theoretical determination of $B_K$, the $\bar K^0-K^0$ hadronic parameter,
which should be consistently determined within every analysis.
The theoretical uncertainty on $B_K$ affects subtantially the
final uncertainty on $\Im \lambda_t$.
A better determination of the unitarity triangle
is expected from the B-factories and the hadronic colliders.
In K-physics, the decay {$K_L\to\pi^0\nu\bar\nu$} gives the cleanest 
``theoretical''
determination of $\Im\lambda_t$, albeit representing a great experimental
challenge.
\end{itemize}


\section{Summary of theory results}
A satisfactory approach to the
calculation of \eoe should comply with the following requirements: 
\begin{itemize}
\item[A:]
A consistent definition of renormalized operators leading to
the correct scheme and scale matching with the short-distance
perturbative analysis.\\[-3ex]
\item[B:]
A self-contained calculation of $all$ relevant hadronic matrix elements
(including $B_K$).\\[-3ex]
\item[C:]
A simultaneous explanation of the $\Delta I = 1/2$ selection rule
and \eoe.
\end{itemize}

None of the existing calculations satisfies all previous 
requirements. I summarize very briefly the various attempts
to calculate \eoe which have appeared so far leading to the
estimates shown in Figs. \ref{fig:pre}
and \ref{fig:post}.
\begin{itemize}
\item VSA:
A simple and naive approach to the problem is
the {VSA}, which is based on two drastic assumptions:
the factorization of the four quark operators
in products of currents or densities and the saturation of 
the completeness of the intermediate states by the vacuum. 
As an example: 
\bea
\langle \pi^+ \pi^-|Q_6| K^0 \rangle & = & 
 2\  \langle \pi^-|\overline{u}\gamma_5 d|0 \rangle
\langle \pi^+|\overline{s} u |K^0 \rangle 
- 2\  \langle \pi^+ \pi^-|\overline{d} d|0 \rangle  
\langle 0|\overline{s} \gamma_5 d |K^0 \rangle
\nonumber \\
& & +\ 2  \left[\langle 0|\overline{s} s|0 \rangle - 
\langle 0|\overline{d}d|0 \rangle\right] 
\langle \pi^+ \pi^-|\overline{s}\gamma_5 d |K^0 \rangle
\eea
The VSA does not allow
for a consistent matching of the scale and scheme dependence
of the Wilson coefficients (the HV and NDR results are shown in
Fig. \ref{fig:pre}) and it carries potentially
large systematic uncertainties~\cite{review}. 
It is best used for LO estimates. 

\item Taipei's:
Generalized factorization represents an attempt to 
parametrize the hadronic matrix elements
in the framework of factorization 
without a-priori assumptions~\cite{neubertstech}. 
Phenomenological parameters are introduced to account for
non-factorizable effects. 
Experimental data are used in order to extract
as much information as possible on the non-factorizable parameters.
This approach has been applied to the $K\to\pi\pi$ amplitudes
in ref.~\cite{Cheng}. 
The effective Wilson coefficients,
which include the perturbative QCD running
of the quark operators,
are matched to
the factorized matrix elements at the scale {$\mu_F$}
which is arbitrarily chosen in the perturbative
regime. A residual scale dependence remains in the penguin matrix 
elements via the quark masses. The analysis shows that in order to
reproduce the $\Delta I = 1/2$ rule and
\ratio\ sizable non-factorizable
contributions are required both in the current-current
and the penguin matrix elements. However, some assumptions on the 
phenomenological parameters and ad hoc subtractions
of scheme-dependent terms in the Wilson coefficients make the
numerical results questionable.
The quoted error does not include any short-distance uncertainty.  

\item M\"unchen's:
In the M\"unchen approach 
(phenomenological $1/N$) some of the matrix elements are obtained by
fitting the $\Delta I =1/2$ rule at $\mu=m_c=1.3$ GeV.
On the other hand, the relevant gluonic and electroweak penguin 
$\vev{Q_6}$ and $\vev{Q_8}_2$ remain undetermined and 
are taken around their leading $1/N$ values 
(which implies a scheme dependent result).
In Fig. \ref{fig:post} the HV (left) and NDR (right) results
are shown~\cite{Buras}. The dark range represents the result of gaussian
treatment of the input parameters compared to flat scannning 
(complete range). 

\item Dortmund's:
In the recent years the Dortmund group has revived and improved
the approach of Bardeen, Buras and Gerard~\cite{BBG} based on the $1/N$
expansion. 
Chiral loops are regularized via a cutoff 
and the amplitudes are arranged in a $p^{2n}/N$ expansion. 
A particular attention has been given
to the matching procedure between the scale dependence of the chiral loops
and that arising from the short-distance analysis~\cite{Dortmund}.  
The renormalization-scheme dependence remains and it is included in the 
final uncertainty. The $\Delta I = 1/2$ rule is reproduced, but the
presence of the quadratic cutoff induces a matching scale instability
(which is very large for $B_K$).
The NLO corrections to $\vev{Q_6}$ induce a substantial enhancement
of the matrix element (right bar in Fig. \ref{fig:post}) 
compared to the leading order result (left bar).
The darker ranges correspond to central values of
$m_s$, $\Omega_{\eta+\eta'}$, $\Im \lambda_t$ and $\Lambda_{QCD}$. 

\item Dubna's:
In the Nambu, Jona-Lasinio (NJL) modelling of QCD~\cite{NJL} 
the Dubna group~\cite{Belkov}
has calculated \eoe including
chiral loops up to $O(p^6)$ and the effects of 
scalar, vector and axial-vector resonances.
Chiral loops are regularized via
the heat-kernel method, which leaves unsolved the
problem of the renormalization-scheme dependence.
A phenomenological fit of the $\Delta I = 1/2$ rule implies deviations up to 
a factor two on the calculated $\vev{Q_6}$.  
The reduced (dark) range in Fig. \ref{fig:post} corresponds to taking the   
central values of the NLO chiral couplings and varying the short-distance 
parameters.

\item Trieste's:
In the approach of the Trieste group, based on the
Chiral Quark Model ($\chi$QM)~\cite{chiQM}, 
all hadronic matrix elements are computed up to $O(p^4)$ in the chiral 
expansion in terms of the three model parameters:
the constituent quark mass, the quark condensate
and the gluon condensate.
These parameters are
phenomenologically fixed by fitting the $\Delta I =1/2$ rule~\cite{ts98a}.
This step is crucial in order to make the model predictive, 
since there is no a-priori argument
for the consistency of the matching procedure
(dimensional regularization and minimal subtraction are used in the
effective chiral theory). 
As a matter of fact, all computed observables turn
out to be very weakly dependent on the
scale (and the renormalization scheme)
in a few hundred MeV range around the 
matching scale, which is taken to be $0.8$ GeV as a
compromise between the ranges of validity 
of perturbation theory and chiral expansion.
The $I=0$ matrix elements are strongly enhanced by non-factorizable
chiral corrections and drive \eoe in the $10^{-3}$ regime.
The dark (light) ranges in Fig. \ref{fig:post} correspond 
to Gaussian (flat) scan of the input parameters. 
The bar on the left represents the result
of ref.~\cite{ts00a} which updates the 1997 calculation~\cite{ts98b}.
That on the right is a new estimate~\cite{ts00b}, similarly 
based on the $\chi$QM hadronic matrix elements, in which however
\ratio\ is computed
by including the  explicit computation of \eps\ in the ratio as
opposed to  the usual procedure of taking its value from the experiments. 
This approach has the advantage of
being independent from the determination of the 
CKM parameters $\Im \lambda_t$ and of showing
more directly the dependence on the long-distance parameter $\hat B_K$
as determined within the model.

\item Roma's:
Lattice regularization of QCD
is $the$ consistent approach to the problem. 
On the other hand, there are presently important numerical and 
theoretical limitations, like the quenching approximation
and the implementation of chiral symmetry, which 
may substantially affect the calculation of the weak matrix elements. 
In addition, chiral perturbation
theory is needed in order to obtain $K\to\pi\pi$ amplitudes  
from the computed $K\to\pi$ transitions.
As summarized in ref.~\cite{Roma}
lattice cannot provide us at present with reliable calculations
of the $I=0$ penguin operators relevant to \eoe, as well as of the 
$I=0$ components of the hadronic matrix elements of the
current-current operators (penguin contractions), which are relevant to the
$\Delta I = 1/2$ rule. 
This is due to large renormalization uncertainties,
partly related to the breaking of chiral symmetry on the lattice.
In the recent Roma re-evaluation of \eoe 
$\vev{Q_6}$ is taken at the VSA value with a 100\% uncertainty~\cite{Roma}.
The result is therefore scheme dependent (the HV and NDR results
are shown in Fig. \ref{fig:post}). The dark (light) ranges correspond 
to Gaussian (flat) scan of the input parameters. 

\item Montpellier's:
The analysis in ref.~\cite{Narison} is based on QCD Sum Rules and uses
recent data on the $\tau$ hadronic total decay rates. The value of
the $Q_8$ matrix element thus found is substantially larger
than the leading $1/N$ result. At the same time, the matrix element
of the $Q_6$ gluonic penguin, computed assuming scalar meson dominance, 
is found in agreement with leading order $1/N$.
The combined effect is a strong cancellation between electroweak 
and gluonic penguins which leads to a vanishingly small \ratio.
Various sources of uncertainties in the calculation and the comparison
with other analyses are discussed in~\cite{Narison}. 

\item Lund's:
The $\Delta I = 1/2$ rule and $B_K$ have been 
studied in the NJL framework and $1/N$ expansion
by Bijnens and Prades~\cite{BijnensPrades} showing an impressive scale 
stability when including vector and axial-vector resonances.
The same authors have
recently produced a calculation of \ratio\ at the 
NLO in $1/N$~\cite{lund}. The calculation is 
done in the chiral limit and it is eventually corrected 
by estimating the largest $SU(3)$ breaking effects. 
Particular attention is devoted to the matching between
long- and short-distance components by use of the $X$-boson 
method~\cite{X-boson0,X-boson1}. The couplings of the $X$-bosons are computed
within the ENJL model which improves the high-energy behavior.
The $\Delta I = 1/2$ rule is reproduced and the computed amplitudes show a 
satisfactory renormalization -scale and -scheme stability. 
A sizeable enhancement of the
$Q_6$ matrix element is found which brings the central value of \ratio
at the level of $3\times 10^{-3}$.

\item Valencia's:
The standard model estimate given by Pallante and Pich is obtained
by applying the FSI correction factors obtained using a dispersive
analysis \`a la Omn\`es-Mushkelishvili~\cite{Omnes} to the leading (factorized)
$1/N$ amplitudes. The detailed numerical outcome has been questioned on the
basis of ambiguities related to the choice of the subtraction
point at which the factorized amplitude is taken~\cite{burasetal}.
Large corrections may also be induced by unknown local terms
which are unaccounted for by the dispersive resummation of the
leading chiral logs. 
Nevertheless, the analysis of ref.~\cite{dispersive} confirms the 
crucial role of higher order chiral corrections for \ratio, 
even though FSI effects alone leave
the problem of reproducing the $\Delta I = 1/2$ selection rule open. 
\end{itemize}

Other attempts to reproduce the measured \eoe using the linear 
$\sigma$-model, which include the effect of a scalar resonance 
with $m_\sigma \simeq 900$ MeV, obtain the needed enhancement
of $\vev{Q_6}$~\cite{sigmamodel}. However, it is not possible
to reproduce simultaneously the experimental values of \eoe and of
the CP conserving $K\to\pi\pi$ amplitudes. 

Studies on the matching between long- and short- distances in large
$N$ QCD, with the calculation of the $Q_7$ penguin matrix element
and of $\hat B_K$ at the NLO in the $1/N$ expansion
have been presented in ref.~\cite{PerisEdR}. However,
a complete calculation of the $K\to\pi\pi$ matrix elements
relevant to \eoe is not available yet.

It has been recently emphasized~\cite{dim8} 
that cut-off based approaches should 
pay attention to higher-dimension
operators which become relevant for matching scales below 2 GeV
and may represent one of the largest sources of uncertainty in present
calculations.
On the other hand, the calculations based on dimensional regularization (for
instance the Trieste one) 
may be safe in this respect if phenomenological input is used 
in order to encode 
in the hadronic matrix elements the physics at all scales. 

Lattice, as a regularization of QCD, is {\em the} first-principle
approach to the problem. 
Presently, very promising developments are being undertaken
to circumvemt the technical and conceptual shortcomings related
to the calculation of weak matrix elements (for a recent survey see
ref.~\cite{Lellouch}).
Among these are the Domain Wall Fermion
approach~\cite{DWF,Soni} which allows us to decouple the chiral symmetry 
from the continuum limit, and the possibility to circumvent
the Maiani-Testa theorem~\cite{MaianiTesta} using the fact that
lattice calculations are performed in finite volume~\cite{LellouchLuscher}.
All these developments need a tremendous effort in machine power
and in devising faster algorithms. Preliminary results 
for the calculations of both \ratio 
and the $\Delta I = 1/2$ selection rule
are expected during the next year. 
 
  
%%%%%%%%%

\section{The $\Delta I = 1/2$ selection rule}
Without entering into the details of the various calculations
I wish to illustrate with a simple exercise the crucial
role of higher order chiral corrections (in general of 
non-factorizable contributions) for \eoe and the $\Delta I = 1/2$ 
selection rule.
In order to do that I focus on two semi-phenomenological approaches.

\begin{figure}[t]
\begin{center}
\epsfig{file=istobuqm.eps,height=7.0cm}
\caption[0]{
Predicting $\varepsilon'/\varepsilon$:
a (Penguin) Comparative Anatomy of the M\"unchen (dark gray) and 
Trieste (light gray) results in Fig. \ref{fig:pre} (in units of $10^{-3}$).
\label{istobuqm}}
\end{center}
\end{figure}

A commonly used way of comparing the estimates of hadronic matrix elements
in different approaches is via the so-called $B$ factors which
represent the ratio of the model matrix elements to the corresponding VSA
values. However, care must be taken in the comparison of
different models due to the scale
dependence of the $B$'s and the values used by different groups
for the parameters that enter the VSA expressions. 
An alternative pictorial and synthetic
way of analyzing different outcomes for \eoe
is shown in Fig.~\ref{istobuqm}, where a ``comparative
anatomy'' of the early Trieste and M\"unchen predictions
is presented. 

From the inspection of the various contributions it is apparent that
the different outcome on the central value of \eoe is almost entirely
due to the difference in the size of the $Q_6$ contribution. 

In the M\"unich approach~\cite{Buras} the $\Delta I = 1/2$ rule
is used in order to determine phenomenologically the  
matrix elements of $Q_{1,2}$ and, 
via operatorial relations, some of the matrix elements of the
left-handed penguins. The approach does not allow
for a phenomenological determination of the matrix elements of the penguin
operators which are most relevant for \eoe, namely the gluonic penguin $Q_6$
and the electroweak penguin $Q_8$. These matrix elements are taken
around their leading $1/N$ values (factorization). 

In the semi-phenomenological approach of the Trieste group 
the size of the effects on the $I=0,2$ amplitudes is controlled
by the phenomenological embedding of the $\Delta I= 1/2$ selection rule
which determines the ranges of the model paremeters: 
the constituent quark mass, the quark and the gluon condensates.
In terms of these parameters all matrix elements are computed.

\begin{figure}[t]
\begin{center}
\epsfig{file=road.eps,height=7.0cm}
\caption[0]{
Anatomy of the $\Delta I = 1/2$ rule in the 
$\chi$QM. See the text for explanations.
The cross-hairs indicate the experimental point. 
\label{road}}
\end{center}
\end{figure}

Fig.~\ref{road} shows an anatomy of the $\chi$QM contributions
which lead to the experimental value
of the $\Delta I = 1/2$ selection rule for central values of the 
input parameters.
Point (1) represents the result obtained by neglecting QCD
and taking the factorized matrix element for the
tree-level operator $Q_2$, which is the leading electroweak 
contribution. 
The ratio $A_0/A_2$ is thus found equal
to $\sqrt{2}$: by far off the experimental point (8).
Step (2) includes the effects of perturbative QCD renormalization
on the operators $Q_{1,2}$~\cite{Gaillard-etc}. 
Step (3) shows the effect of including the gluonic 
penguin operators~\cite{VSZ-GWise-CFGeorgi}. 
Electroweak penguins~\cite{Lusignoli-etc} are numerically
negligeable in the CP conserving amplitudes and 
are responsible for the very small shift in the $A_2$ direction.
Therefore, perturbative QCD and factorization lead us from (1) to (4):
a factor five away from the experimental ratio.

Non-factorizable gluon-condensate corrections,
a crucial model dependent effect entering at the leading order
in the chiral expansion, produce a substantial
reduction of the $A_2$ amplitude (5), 
as it was first observed by Pich and de Rafael~\cite{Pich-deRafael}. 
Moving the analysis to $O(p^4)$,
the chiral loop corrections, computed on the LO chiral
lagrangian via dimensional regularization and minimal subtraction, 
lead us from (5) to (6), while 
the finite parts of the NLO counterterms
calculated in the $\chi$QM approach lead to the point (7).
Finally, step (8) represents the inclusion of $\pi$-$\eta$-$\eta'$ 
isospin breaking effects~\cite{Ometapeta}. 

This model dependent anatomy 
shows the relevance of non-factorizable contributions
and higher-order chiral corrections. The suggestion that
chiral dynamics may be relevant to
the understanding of the $\Delta I = 1/2$ selection rule goes back to the
work of Bardeen, Buras and Gerard~\cite{BBG} in the $1/N$
framework with a cutoff regularization. 
A pattern similar to that shown
in Fig.~\ref{road} for the chiral loop corrections to $A_0$ and $A_2$
was previously obtained, using dimensional regularization in a NLO chiral
lagrangian analysis, by Missimer, Kambor and Wyler~\cite{MKWyler}.
The Trieste group has extended their calculation to include
the NLO contributions
to the electroweak penguin matrix elements~\cite{ts96}.

\begin{figure}[t]
\begin{center}
\epsfig{file=charteps.eps,height=7.0cm}
\caption[0]{
Anatomy of \eoe (in units of $10^{-3}$) in the $\chi$QM approach. 
In black the LO results (which includes the non-factorizable gluonic
corrections), in half-tone the effect of the inclusion of chiral-loop 
corrections and in light gray the complete $O(p^4)$ estimate.
\label{charteps}}
\end{center}
\end{figure}

Fig.~\ref{charteps} shows the contributions to \eoe
of the various penguin operators, providing us
with a finer anatomy of the NLO chiral corrections.
It is clear that chiral-loop dynamics plays a subleading role in the 
electroweak penguin sector ($Q_{8-10}$) while enhancing by 60\% the gluonic
penguin ($I=0$) matrix elements. 
The NLO enhancement of the $Q_6$ matrix element is what drives \eoe 
in the $\chi$QM to the $10^{-3}$ ballpark.

As a consequence, the $\chi$QM analysis shows that the same dynamics 
that is relevant to the
reproduction of the CP conserving $A_0$ amplitude 
(Fig.~\ref{road}) is at work
in the CP violating sector, albeit with a reduced strenght. 

In order to ascertain
whether these model features represent real QCD effects we must
wait for future improvements in lattice calculations~\cite{Lellouch}.
On the other hand, indications for such a dynamics follow also from
the analitic properties
of the $K\to\pi\pi$ amplitudes~\cite{dispersive}.
It is important to notice however that the size of
the effect so derived is generally
not enough to fully account for the $\Delta I = 1/2$ rule.
Other non-factorizable contributions are needed to further
enhance the CP conserving $I=0$ amplitude, and to
reduce the large $I=2$ amplitude obtained from perturbative QCD
and factorization.
In the $\chi$QM approach, for instance, the fit of the  $\Delta I = 1/2$ rule
is due to the interplay of NLO chiral corrections 
and non-factorizable soft gluonic
contributions (at LO in the chiral expansion). 


\section{Conclusions}
In summary,
those semi-phenomenological approaches which reproduce the
$\Delta I = 1/2$ selection rule in $K\to\pi\pi$ decays,
generally agree in the pattern and size    
of the $I = 2$ hadronic matrix elements with the existing
lattice calculations. 
On the other hand,
the $\Delta I = 1/2$ rule forces upon us large deviations from
the naive factorization for the $I=0$ amplitudes: 
$B-$factors of $O(10)$ are required for $\vev{Q_{1,2}}_0$.
Here is were 
lattice calculations presently suffer from large sistematic uncertainties. 

In the Trieste and Dortmund calculations, which reproduce the CP conserving
$K\to\pi\pi$ amplitudes, non-factorizable effects
(mainly due to final-state interactions) enhance the hadronic matrix 
elements of the gluonic penguins, and give $B_6/B_8^{(2)} \approx 2$. 
Similar indications stem from recent $1/N$~\cite{lund} 
and dispersive~\cite{dispersive} approaches. 
The direct calculation of $K\to\pi\pi$ amplitudes and unquenching
are needed in the lattice calculations in order to account
for final state interactions. Promising and exiciting work in this direction
is in progress.


\Acknowledgements
I wish to thank the organizers for their invitation to this interesting
meeting hosted in a most enjoyable environment.


\begin{thebibliography}{99}

%%
%%  bibliographic items can be constructed using the LaTeX format in SPIRES:
%%    see    http://www.slac.stanford.edu/spires/hep/latex.html
%%  SPIRES will also supply the CITATION line information; please include it.
%%

\bibitem{KTeV} 
A. Alavi-Harati {\ et al.}, \Journal{\PRL}{83}{22}{1999}; 
A. Barker, these Proceedings.

\bibitem{NA48} 
V. Fanti {\ et al.}, \Journal{\PLB}{465}{335}{1999};  
J. Ocariz, {\sl JHEP Proc.}, 
Third Latin American Symposium on High Energy
Physics, Cartagena de Indias, Colombia, 2-8 April 2000
[PREHEP-silafae-III/027].

\bibitem{Christenson} 
J.H. Christenson {\ et al.}, \Journal{\PRL}{13}{138}{1964}.

\bibitem{review} 
S. Bertolini, J.O. Eeg and M. Fabbrichesi, 
\Journal{Rev. Mod. Phys.}{72}{65}{2000}.

\bibitem{Fermi} E. Fermi, \Journal{Suppl. Nuovo Cim.}{2}{17}{1955}.

\bibitem{dispersive}
E. Pallante and A. Pich, \Journal{\PRL}{84}{2568}{2000};
\Journal{\NPB}{592}{294}{2000};
A. Pich, these Proceedings. 
%J.F. Donoghue and E. Golowich, \Journal{\PLB }{478}{172}{2000};
%E. A. Paschos, .

\bibitem{Omnes} 
N.I. Mushkelishvili, {\sl Singular Integral Equations}
(Noordhoff, Gronigen, 1953), p. 204; 
R. Omn\`es, \Journal{Nuovo Cim.}{8}{316}{1958}. 

\bibitem{burasetal} {A.~J. Buras et al.}, \Journal{\PLB}{480}{80}{2000}. 

\bibitem{Martinelli} M. Ciuchini and G. Martinelli, .

\bibitem{Pais-Gell-Mann} 
M. Gell-Mann and A. Pais, \Journal{Proc. Glasgow Conf.}{}{342}{1955}.

\bibitem{FSIphases} 
J. Gasser and U.G. Meissner, \Journal{\PLB}{258}{219}{1991}; 
E. Chell and M.G. Olsson, \Journal{\PRD}{48}{4076}{1993}.

\bibitem{GardnerValencia} S. Gardner and G. Valencia, 
\Journal{\PLB }{466}{355}{1999}.

\bibitem{etapiNLO}{ G. Ecker et al.}, \Journal{\PLB}{467}{88}{2000}.

\bibitem{omegaNLOstr}{ S. Gardner and G. Valencia}, 
\Journal{\PRD}{62}{094024}{2000}.

\bibitem{omegaNLOmodel}{ K. Maltman and C. Wolfe},
\Journal{\PLB}{482}{77}{2000}; \Journal{\PRD}{63}{014008}{2001}.

\bibitem{omegaNLOem}{ V. Cirigliano, J.F. Donoghue and E. Golowich},
\Journal{\PLB}{450}{241}{1999}; \Journal{\PRD}{61}{093001}{2000};
\Journal{\PRD}{61}{093002}{2000}; \Journal{Eur. Phys. J. C}{18}{83}{2000}.

\bibitem{Ciuchinietal} M. Ciuchini {et al.}, .

\bibitem{neubertstech} {M. Neubert and B. Stech}, in "Heavy Flavours", 
edited by A.J. Buras and M. Lindner (World Scientific, Singapore),
.

\bibitem{Cheng} H.Y. Cheng, \Journal{Chin. J. Phys.}{38}{1044}{2000}.

\bibitem{Buras}
A.J. Buras, Proceedings of the Kaon 99 Conference, Chicago, USA,
. 

\bibitem{BBG} W.A. Bardeen, A.J. Buras and J.M. Gerard, 
\Journal{\NPB }{293}{787}{1987}.

\bibitem{Dortmund} 
T. Hambye {\ et al.}, \Journal{\PRD}{58}{014017}{1998}; 
\Journal{\NPB }{564}{391}{2000}.

\bibitem{NJL} Y. Nambu and G. Jona-Lasinio, 
\Journal{Phys. Rev.}{122}{345}{1961}; J. Bijnens,  
\Journal{Phys. Rept.}{265}{369}{1996}.

\bibitem{Belkov} A.A. Bel'kov {\ et al.}, ; 

\bibitem{chiQM} 
S. Weinberg, \Journal{Physica} A {96}{327}{1979};
A. Manhoar and H. Georgi, \Journal{\NPB }{234}{189}{1984}.

\bibitem{ts98a} 
S. Bertolini {\ et al.}, \Journal{\NPB }{514}{63}{1998}. 

\bibitem{ts98b} 
S. Bertolini {\ et al.}, \Journal{\NPB }{514}{93}{1998}. 

\bibitem{ts00a} 
S. Bertolini, J.O. Eeg and M. Fabbrichesi,  to appear
in \Journal{\PRD}{63}{0560xx}{2001}. 

\bibitem{ts00b} 
M. Fabbrichesi, \Journal{\PRD}{62}{097902}{2000}.

\bibitem{Roma} G. Martinelli, Proceedings of the Kaon 99 Conference, 
Chicago, USA, .

\bibitem{Narison} S. Narison, \Journal{\NPB}{593}{3}{2001}. 

\bibitem{BijnensPrades} J. Bijnens and J. Prades, 
\Journal{JHEP}{9901}{023}{1999}.

\bibitem{lund} J. Bijnens and J. Prades, \Journal{JHEP}{0006}{035}{2000}.

\bibitem{X-boson0}{W.~A.~Bardeen, J.~Bijnens and J.-M.~G\'erard},
\Journal{\PRL}{62}{1343}{1989};
{J.~Bijnens, J.-M.~G\'erard and G.~Klein},
\Journal{\PLB}{257}{191}{1991};
{J.~P.~Fatelo and J.-M.~G\'erard},
\Journal{\PLB}{347}{136}{1995}.

\bibitem{X-boson1}{J.~Bijnens and J.~Prades},
\Journal{\PLB}{342}{331}{1995};
\Journal{\NPB}{444}{523}{1995};
\Journal{JHEP}{9901}{023}{1999};
\Journal{JHEP}{0001}{002}{2000}.

\bibitem{sigmamodel} Y.Y. Keum, U. Nierste and A.I. Sanda,
\Journal{\PLB }{457}{157}{1999}; 
M. Harada {\ et al.}, \Journal{\PRD}{62}{014002}{2000};
J.C.R. Bloch {\ et al.}, \Journal{Phys. Rev. C}{62}{025206}{2000}.

\bibitem{PerisEdR} 
S. Peris, M. Perrottet and E. de Rafael, \Journal{JHEP}{9805}{011}{1998}; 
M. Knecht, S. Peris and E. de Rafael, \Journal{\PLB}{443}{255}{1998},
\Journal{\PLB}{457}{227}{1999};
S. Peris and E. de Rafael, \Journal{\PLB}{490}{213}{2000}.

\bibitem{dim8}{ V. Cirigliano, J.F. Donoghue and E. Golowich},
\Journal{JHEP}{0010}{048}{2000}. 

\bibitem{Lellouch} L. Lellouch, Proceedings of Lattice 2000, Bangalore,
India, .

\bibitem{DWF} D. Kaplan, \Journal{\PLB}{288}{342}{1992}.

\bibitem{Soni} A. Soni, Proceedings of the $B_{CP}$ 99 Conference, 
Taipei, Taiwan, . 

\bibitem{MaianiTesta} L. Maiani and M. Testa, \Journal{\PLB}{245}{585}{1990}.

\bibitem{LellouchLuscher} L. Lellouch and M. L\"uscher, .

%%%%
\bibitem{Gaillard-etc} 
M.K. Gaillard and B. Lee, \Journal{\PRL}{33}{108}{1974}; 
G. Altarelli and L. Maiani, \Journal{\PLB }{52}{413}{1974}. 
 
\bibitem{VSZ-GWise-CFGeorgi} M.A. Shifman, A.I. Vainshtein and V.I. Zacharov,
\Journal{JETP Lett.}{22}{55}{1975}; \Journal{\NPB }{120}{316}{1977}; 
\Journal{JETP}{45}{670}{1977};
J. Ellis {\ et al.},
\Journal{\NPB }{131}{285}{1977}, (E) {132}{541}{1978};
F.J. Gilman and M.B. Wise, \Journal{\PLB }{83}{83}{1979};
\Journal{\PRD }{20}{2392}{1979}; 
R.S. Chivukula, J.M. Flynn and H. Georgi, \Journal{\PLB }{171}{453}{1986}.

\bibitem{Lusignoli-etc} 
J. Bijnens and M.B. Wise, \Journal{\PLB }{137}{245}{1984}; 
A.J. Buras and J.M. Gerard, \Journal{\PLB }{192}{156}{1987}; 
M. Lusignoli, \Journal{\NPB }{325}{33}{1989}; 
J.M. Flynn and L. Randall, \Journal{\PLB }{224}{221}{1989}.   

\bibitem{Pich-deRafael} 
A. Pich and E. de Rafael, \Journal{\NPB }{358}{311}{1991}.

\bibitem{Ometapeta} B. Holstein, \Journal{\PRD }{20}{1187}{1979}; 
J.F. Donoghue, E. Golowich, B.R. Holstein and J. Trampetic, 
\Journal{\PLB }{179}{361}{1986}; 
A.J. Buras and J.M. Gerard, in ref.~\cite{Lusignoli-etc}.

\bibitem{MKWyler} J. Kambor, J. Missimer and D. Wyler, 
\Journal{\NPB }{346}{17}{1990}; \Journal{\PLB }{261}{496}{1991}. 

\bibitem{ts96} 
V. Antonelli {\ et al.}, \Journal{\NPB }{469}{143}{1996};
\Journal{\NPB }{469}{181}{1996}; 
S. Bertolini, J.O. Eeg and M. Fabbrichesi, \Journal{\NPB }{476}{225}{1996}; 

\end{thebibliography}

 
\end{document}


