• We present a lattice determination of the $\Lambda$ parameter in three-flavor 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 four-loop 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 3-flavour 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 5-flavour QCD.
  • We apply the Symanzik improvement programme to the 4+1-dimensional local re-formulation 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 on-shell improvement. Compared to the 4-dimensional 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 non-trivial order in the coupling.
  • We review our new strategy and current status towards a high precision computation of the Lambda parameter from three-flavour 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 one-component 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 four-component spinors, two different set-ups 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 tree-level 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 finite-size recursion techniques based on the Schrodinger functional formalism to determine the renormalization group running of four-fermion 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 four-fermion 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$ four-fermion operator is renormalised multiplicatively. The renormalisation adopted is non-perturbative, with a Schroedinger functional renormalisation condition. Renormalisation group running is also non-perturbative, 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 non-degenerate 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 four-quark 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 non-perturbatively 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 four-quark 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 Wilson-type 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 one-loop order of perturbation theory. In this way the two-loop anomalous dimensions in the SF schemes can be inferred. As a by-product of our calculation we also obtain the one-loop cutoff effects in the step-scaling 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 non-perturbative 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 long-standing renormalisation problems of the four-quark operators which contribute to CP-conserving $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 four-fermion 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 four-fermion operators which appear in the effective weak Hamiltonian of the Standard Model. Our calculations are done using Wilson fermions and the parity-odd 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 non-perturbatively. The calculation is based on an extension of a finite-size technique previously employed to compute the running coupling in quenched QCD. As a by-product we obtain the $\Lambda$--parameter in this theory with completely controlled errors.
  • It is shown how on-shell O(a) improvement can be implemented non-perturbatively 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.