We derive chiral Ward identities for lattice QCD with Wilson quarks and $$N_{\mathrm{f}}\ge 3$$Nf≥3 flavours, on small lattices with Schrödinger functional boundary conditions and vanishingly small quark masses. These identities relate the axial variation of the non-singlet pseudoscalar density to the scalar one, thus enabling the non-perturbative determination of the scale-independent ratio $$Z_{\mathrm {S}}/Z_{\mathrm {P}}$$ZS/ZP of the renormalisation parameters of these operators. We obtain results for $$N_{\mathrm{f}}=3$$Nf=3 QCD with tree-level Symanzik-improved gluons and Wilson-Clover quarks, for bare gauge couplings which cover the typical range of large-volume $$N_{\mathrm{f}}= 2+1$$Nf=2+1 simulations with Wilson fermions at lattice spacings below $$0.1\,$$0.1fm. The precision of our results varies from 0.3 to 1%, except for the coarsest lattice, where it is 2%. We discuss how the $$Z_{\mathrm {S}}/Z_{\mathrm {P}}$$ZS/ZP ratio can be used in the non-perturbative calculations of $${\mathrm {O}}(a)$$O(a) improved renormalised quark masses.

We present the results of a non-perturbative determination of the improvement coefficient $$c_\mathrm{V}$$ c V and the renormalisation factor $$Z_\mathrm{V}$$ Z V , which define the renormalised vector current in three-flavour $$\mathrm{O}(a)$$ O ( a ) improved lattice QCD with Wilson quarks and tree-level Symanzik-improved gauge action. In case of the improvement coefficient, we consider both lattice descriptions of the vector current, the local as well as the conserved (i.e., point-split) one. Our improvement and normalisation conditions are based on massive chiral Ward identities and numerically evaluated in the Schrödinger functional setup, which allows to eliminate finite quark mass effects in a controlled way. In order to ensure a smooth dependence of the renormalisation constant and improvement coefficients on the bare gauge coupling, our computation proceeds along a line of constant physics, covering the typical range of lattice spacings $$0.04\,\mathrm{fm}\lesssim a\lesssim 0.1\,\mathrm{fm}$$ 0.04 fm ≲ a ≲ 0.1 fm that is useful for phenomenological applications. Especially for the improvement coefficient of the local vector current, we report significant differences between the one-loop perturbative estimates and our non-perturbative results.

We determine non-perturbatively the normalisation factor $$r_{\mathrm{m}}\equiv Z_{\mathrm{S}}/Z_{\mathrm{S}}^{0}$$ r m ≡ Z S / Z S 0 , where $$Z_{\mathrm{S}}$$ Z S and $$Z_{\mathrm{S}}^{0}$$ Z S 0 are the renormalisation parameters of the flavour non-singlet and singlet scalar densities, respectively. This quantity is required in the computation of quark masses with Wilson fermions and for instance the renormalisation of nucleon matrix elements of scalar densities. Our calculation involves simulations of finite-volume lattice QCD with the tree-level Symanzik-improved gauge action, $$N_{\mathrm{f}}= 3$$ N f = 3 mass-degenerate $${\mathrm{O}}(a)$$ O ( a ) improved Wilson fermions and Schrödinger functional boundary conditions. The slope of the current quark mass, as a function of the subtracted Wilson quark mass is extracted both in a unitary setup (where nearly chiral valence and sea quark masses are degenerate) and in a non-unitary setup (where all valence flavours are chiral and the sea quark masses are small). These slopes are then combined with $$Z \equiv Z_{\mathrm{P}}/(Z_{\mathrm{S}}Z_{\mathrm{A}})$$ Z ≡ Z P / ( Z S Z A ) in order to obtain $$r_{\mathrm{m}}$$ r m . A novel chiral Ward identity is employed for the calculation of the normalisation factor Z. Our results cover the range of gauge couplings corresponding to lattice spacings below $$0.1\,$$ 0.1 fm, for which $$N_{\mathrm{f}}= 2+1$$ N f = 2 + 1 QCD simulations in large volumes with the same lattice action are typically performed.

We report on non-perturbative computations of the improvement coefficient cV and the renormalization factor ZV of the vector current in three-flavour O(a) improved lattice QCD with Wilson quarks and tree-level Symanzik improved gauge action. To reduce finite quark mass effects, our improvement and normalization conditions exploit massive chiral Ward identities formulated in the Schrödinger functional setup, which also allow deriving a new method to extract the ratio ZS/ZP of scalar to pseudoscalar renormalization constants. We present preliminary results of a numerical evaluation of ZV and cV along a line of constant physics with gauge couplings corresponding to lattice spacings of about 0:09 fm and below, relevant for phenomenological applications.

We present our progress in the non-perturbative O(a) improvement and renormalization of tensor currents in three-flavor lattice QCD with Wilson-clover fermions and tree-level Symanzik improved gauge action. The mass-independent O(a) improvement factor of tensor currents is determined via a Ward identity approach, and their renormalization group running is calculated via recursive finite-size scaling techniques, both implemented within the Schrödinger functional framework. We also address the matching factor between bare and renormalization group invariant currents for a range of lattice spacings < 0.1 fm, relevant for phenomenological large-volume lattice QCD applications.