
We present a lattice determination of the $\Lambda$ parameter in threeflavor
QCD and the strong coupling at the Z pole mass. Computing the nonperturbative
running of the coupling in the range from $0.2\,$GeV to $70\,$GeV, and using
experimental input values for the masses and decay constants of the pion and
the kaon, we obtain $\Lambda_{\overline{\rm MS}}^{(3)}=341(12)\,$MeV. The
nonperturbative running up to very high energies guarantees that systematic
effects associated with perturbation theory are well under control. Using the
fourloop prediction for $\Lambda_{\overline{\rm
MS}}^{(5)}/\Lambda_{\overline{\rm MS}}^{(3)}$ yields
$\alpha^{(5)}_{\overline{\rm MS}}(m_{\rm Z}) = 0.11852(84)$.

We present results by the ALPHA collaboration for the $\Lambda$parameter in
3flavour QCD and the strong coupling constant at the electroweak scale,
$\alpha_s(m_Z)$, in terms of hadronic quantities computed on the CLS gauge
configurations. The first part of this proceedings contribution contains a
review of published material \cite{Brida:2016flw,DallaBrida:2016kgh} and yields
the $\Lambda$parameter in units of a low energy scale, $1/L_{\rm had}$. We
then discuss how to determine this scale in physical units from experimental
data for the pion and kaon decay constants. We obtain $\Lambda_{\overline{\rm
MS}}^{(3)} = 332(14)$ MeV which translates to $\alpha_s(M_Z)=0.1179(10)(2)$
using perturbation theory to match between 3, 4 and 5flavour QCD.

We apply the Symanzik improvement programme to the 4+1dimensional local
reformulation of the gradient flow in pure $SU(N)$ lattice gauge theories. We
show that the classical nature of the flow equation allows to eliminate all
cutoff effects at $\mathcal O(a^2)$ which originate either from the discretized
gradient flow equation or from the gradient flow observable. All the remaining
$\mathcal O(a^2)$ effects can be understood in terms of local counterterms at
the zero flow time boundary. We classify these counterterms and provide a
complete set as required for onshell improvement. Compared to the
4dimensional pure gauge theory only a single additional counterterm is
required, which corresponds to a modified initial condition for the flow
equation. A consistency test in perturbation theory is passed and allows to
determine all counterterm coefficients to lowest nontrivial order in the
coupling.

We review our new strategy and current status towards a high precision
computation of the Lambda parameter from threeflavour simulations in QCD. To
reach this goal we combine specific advantages of the Schr\"odinger functional
and gradient flow couplings.

We consider the Schr\"odinger functional with staggered onecomponent
fermions on a fine lattice of size $(L/a)^3 \times (T/a)$ where $T/a$ must be
an odd number. In order to reconstruct the fourcomponent spinors, two
different setups are proposed, corresponding to the coarse lattice having size
$(L/2a)^3 \times (T'/2a)$, with $T' = T \pm a$. The continuum limit is then
defined at fixed $T'/L$. Both cases have previously been investigated in the
pure gauge theory. Here we define fermionic correlation functions and study
their approach to the continuum limit at treelevel of perturbation theory.

We simulate two variants of quenched twisted mass QCD (tmQCD), with
degenerate Wilson quarks of masses equal to or heavier than half the strange
quark mass. We use Ward identities in order to measure the twist angles of the
theory and thus check the quality of the tuning of mass parameters to a physics
condition which stays constant as the lattice spacing is varied. Flavour
symmetry breaking in tmQCD is studied in a framework of two fully twisted and
two standard Wilson quark flavours, tuned to be degenerate in the continuum.
Comparing pseudoscalar masses, obtained from connected quark diagrams made of
tmQCD and/or standard Wilson quark propagators, we confirm that flavour
symmetry breaking effects, which are at most 5%, decrease as we approach the
continuum limit. We also compute the pseudoscalar decay constant in the
continuum limit, with reduced systematics. As a consequence of improved tuning
of the mass parameters at $\beta = 6.1$, we reanalyse our previous $B_K$
results. Our main phenomenological findings are $r_0 f_K = 0.421(7)$ and $\hat
B_K = 0.735(71)$.

We apply finitesize recursion techniques based on the Schrodinger functional
formalism to determine the renormalization group running of fourfermion
operators which appear in the Delta S=2 effective weak Hamiltonian of the
Standard Model. Our calculations are done using O(a) improved Wilson fermions
with N_f=2 dynamical flavours. Preliminary results are presented for the
fourfermion operator which determines the B_K parameter in tmQCD.

The $B_K$ parameter is computed in quenched lattice QCD with Wilson twisted
mass fermions. Two variants of tmQCD are used; in both of them the relevant
$\Delta S = 2$ fourfermion operator is renormalised multiplicatively. The
renormalisation adopted is nonperturbative, with a Schroedinger functional
renormalisation condition. Renormalisation group running is also
nonperturbative, up to very high energy scales. In one of the two tmQCD
frameworks the computations have been performed at the physical $K$meson mass,
thus eliminating the need of mass extrapolations. Simulations have been
performed at several lattice spacings and the continuum limit was reached by
combining results from both tmQCD regularisations. Finite volume effects have
been partially checked and turned out to be small. Exploratory studies have
also been performed with nondegenerate valence flavours. The final result for
the RGI bag parameter, with all sources of uncertainty (except quenching) under
control, is $\hat B_K =0.789 \pm 0.046$.

We define a family of Schroedinger Functional renormalization schemes for the
fourquark multiplicatively renormalizable operators of the $\Delta F = 1$ and
$\Delta F = 2$ effective weak Hamiltonians. Using the lattice regularization
with quenched Wilson quarks, we compute nonperturbatively the renormalization
group running of these operators in the continuum limit in a large range of
renormalization scales. Continuum limit extrapolations are well controlled
thanks to the implementation of two fermionic actions (Wilson and Clover). The
ratio of the renormalization group invariant operator to its renormalized
counterpart at a low energy scale, as well as the renormalization constant at
this scale, is obtained for all schemes.

Starting from the QCD Schroedinger functional (SF), we define a family of
renormalization schemes for two fourquark operators, which are, in the chiral
limit, protected against mixing with other operators. With the appropriate
flavour assignments these operators can be interpreted as part of either the
$\Delta F=1$ or $\Delta F=2$ effective weak Hamiltonians. In view of lattice
QCD with Wilsontype quarks, we focus on the parity odd components of the
operators, since these are multiplicatively renormalized both on the lattice
and in continuum schemes. We consider 9 different SF schemes and relate them to
commonly used continuum schemes at oneloop order of perturbation theory. In
this way the twoloop anomalous dimensions in the SF schemes can be inferred.
As a byproduct of our calculation we also obtain the oneloop cutoff effects
in the stepscaling functions of the respective renormalization constants, for
both O(a) improved and unimproved Wilson quarks. Our results will be needed in
a separate study of the nonperturbative scale evolution of these operators.

Twisted mass lattice QCD (tmQCD), generalised to four Wilson quark flavours,
can be used for the computation of some weak matrix elements related to $\Delta
I=1/2$ transitions. Besides eliminating unphysical zero modes, tmQCD may
alleviate longstanding renormalisation problems of the fourquark operators
which contribute to CPconserving $K\to\pi$ transitions. With an active charm
quark, the renormalisation of the $K\to\pi$ matrix elements requires at most
the subtraction of a linearly divergent counterterm. Furthermore, in the
(partially) quenched approximation the twist angles can be chosen so that only
a finite counterterm needs to be subtracted.

We present the results of a precision computation of B_K with Wilson
fermions. Simulations are performed at different lattice spacings, enabling
continuum limit extrapolations. Two different twisted mass QCD (tmQCD)
regularisations are considered for the computation of bare matrix elements. In
both cases the relevant fourfermion operator renormalises multiplicatively. In
one regularisation it is possible to perform the computation directly at the
physical kaon mass value, thus avoiding extrapolations in the mass.
Nonperturbative renormalisation is carried out using available Schroedinger
Functional results.

We present some preliminary results for $B_K$ at $\beta=6.0$, using the
twisted mass QCD formalism for the computation of bare matrix elements of the
$\Delta S=2$ operator. The main advantage of the method is that mixing with
other $d=6$ operators under renormalisation is avoided. Moreover the operator
renormalisation is performed in the Schr\"odinger functional (SF) framework,
using earlier results of our collaboration for the corresponding step scaling
function.

We show that the application of twisted mass QCD (tmQCD) with four (Wilson)
quark flavours to the computation of lattice weak matrix elements relevant to
$\Delta I=1/2$ transitions has important advantages: the renormalisation of $K
\to \pi$ matrix elements does not require the subtraction of other dimension
six operators, the divergence arising from the subtraction of lower dimensional
operators is softened by one power of the lattice spacing and quenched
simulations do not suffer from exceptional configurations at small pion mass.
This last feature is also retained in the tmQCD computation of $K \to \pi\pi$
matrix elements, which, as far as renormalisation and power subtractions are
concerned, has properties analogous to the standard Wilson case.

We apply the Schroedinger Functional (SF) formalism to determine the
renormalisation group running of fourfermion operators which appear in the
effective weak Hamiltonian of the Standard Model. Our calculations are done
using Wilson fermions and the parityodd components of the operators.
Preliminary results are presented for the operator $O_{VA}=(\bar s \gamma_\mu
d)(\bar s \gamma_\mu \gamma_5 d)$.

We present a scaling study of lattice QCD with O(a) improved Wilson fermions
and a chirally twisted mass term. In order to get precise results with a
moderate computational effort, we have considered a system of physical size of
0.75^3 * 1.5 fm^4 with Schroedinger functional boundary conditions in the
quenched approximation. Looking at meson observables in the pseudoscalar and
vector channels, we find that O(a) improvement is effective and residual cutoff
effects are fairly small.

We show that the renormalization factor relating the renormalization group
invariant quark masses to the bare quark masses computed in lattice QCD can be
determined nonperturbatively. The calculation is based on an extension of a
finitesize technique previously employed to compute the running coupling in
quenched QCD. As a byproduct we obtain the $\Lambda$parameter in this theory
with completely controlled errors.

It is shown how onshell O(a) improvement can be implemented
nonperturbatively in lattice QCD with Wilson quarks. Improvement conditions
are obtained by requiring the PCAC relation to hold exactly in certain matrix
elements. These are derived from the QCD Schr\"odinger functional which enables
us to simulate directly at vanishing quark masses. In the quenched
approximation and for bare couplings in the range $0\leq g_0\leq 1$, we
determine the improved action, the improved axial current, the additive
renormalization of the quark mass and the isospin current normalization
constants Z_A and Z_V.