Source author record

Rainer Sommer

Rainer Sommer appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

65works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

65 published item(s)

preprint2022arXiv

A lattice QCD perspective on weak decays of b and c quarks Snowmass 2022 White Paper

Lattice quantum chromodynamics has proven to be an indispensable method to determine nonperturbative strong contributions to weak decay processes. In this white paper for the Snowmass community planning process we highlight achievements and future avenues of research for lattice calculations of weak $b$ and $c$ quark decays, and point out how these calculations will help to address the anomalies currently in the spotlight of the particle physics community. With future increases in computational resources and algorithmic improvements, percent level (and below) lattice determinations will play a central role in constraining the standard model or identifying new physics.

preprint2022arXiv

Investigation of the Perturbative Expansion of Moments of Heavy Quark Correlators for $N_f=0$

The QCD-coupling is a necessary input in the computation of many observables, and the parametric error on input parameters can be a dominant source of uncertainty. The coupling can be extracted by comparing high order perturbative computations and lattice evaluated moments of mesonic two-point functions with heavy quarks, which provide a high energy scale for perturbation theory. The truncation of the perturbative series is an important systematic uncertainty. We report on our attempt to study this issue by measuring pseudo-scalar two-point functions in volumes of $L=2\, \text{fm}$ with twisted-mass Wilson fermions in the quenched approximation. We use full twist, the non-perturbative clover term and lattice spacings down to $a=0.015\,\text{fm}$ to tame the sizable discretization effects. Our preliminary results indicate that either higher order perturbative corrections or the continuum limit are not under sufficient control despite our small lattice spacings and quark masses extending beyond $2\,m_{\text{charm}}$.

preprint2022arXiv

Logarithmic corrections to $\mathbf{a^2}$ scaling in lattice QCD with Wilson and Ginsparg-Wilson quarks

We analyse the leading logarithmic corrections to the $a^2$ scaling of lattice artefacts in QCD, following the seminal work of Balog, Niedermayer and Weisz in the O(n) non-linear sigma model. Limiting the discussion to contributions from the action, the leading logarithmic corrections can be determined by the anomalous dimensions of mass-dimension 6 operators. These operators form a minimal on-shell basis of the Symanzik Effective Theory. We present results for non-perturbatively O($a$) improved Wilson and Ginsparg-Wilson quarks.

preprint2019arXiv

Asymptotic behavior of cutoff effects in Yang-Mills theory and in Wilson's lattice QCD

Discretization effects of lattice QCD are described by Symanzik's effective theory when the lattice spacing, $a$, is small. Asymptotic freedom predicts that the leading asymptotic behavior is $\sim a^n [\bar g^2(a^{-1})]^{\hatγ_1} \sim a^n \left[\frac{1}{-\log(aΛ)}\right]^{\hatγ_1}$. For spectral quantities, $n=d$ is given in terms of the (lowest) canonical dimension, $d+4$, of the operators in the local effective Lagrangian and $\hatγ_1$ is proportional to the leading eigenvalue of their one-loop anomalous dimension matrix $γ^{(0)}$. We determine $γ^{(0)}$ for Yang-Mills theory ($n=2$) and discuss consequences in general and for perturbatively improved short distance observables. With the help of results from the literature, we also discuss the $n=1$ case of Wilson fermions with perturbative O$(a)$ improvement and the discretization effects specific to the flavor currents. In all cases known so far, the discretization effects are found to disappear faster than the naive $\sim a^n$ and the log-corrections are a rather weak modification -- in contrast to the two-dimensional O(3) sigma model.

preprint2019arXiv

Non-perturbative renormalization by decoupling

We propose a new strategy for the determination of the QCD coupling. It relies on a coupling computed in QCD with $N_{\rm f} \geq 3$ degenerate heavy quarks at a low energy scale $μ_{\rm dec}$, together with a non-perturbative determination of the ratio $Λ/μ_{\rm dec}$ in the pure gauge theory. We explore this idea using a finite volume renormalization scheme for the case of $N_{\rm f} = 3$ QCD, demonstrating that a precise value of the strong coupling $α_s$ can be obtained. The idea is quite general and can be applied to solve other renormalization problems, using finite or infinite volume intermediate renormalization schemes.

preprint2016arXiv

Continuum limit of the leading order HQET form factor in $B_s \to K\ellν$ decays

We discuss the computation of form factors for semi-leptonic decays of $\rm B$-, $\rm B_s$- mesons in lattice QCD. Considering in particular the example of the static $\rm B_s$ form factors we demonstrate that after non-perturbative renormalization the continuum limit can be taken with confidence. The resulting precision is of interest for extractions of $V_{\rm ub}$. The size of the corrections of order $1/m_{\rm b}$ is just estimated at present but it is expected that their inclusion does not pose significant difficulties.

preprint2016arXiv

Determination of the QCD $Λ$-parameter and the accuracy of perturbation theory at high energies

We discuss the determination of the strong coupling $α_\mathrm{\overline{MS}}^{}(m_\mathrm{Z})$ or equivalently the QCD $Λ$-parameter. Its determination requires the use of perturbation theory in $α_s(μ)$ in some scheme, $s$, and at some energy scale $μ$. The higher the scale $μ$ the more accurate perturbation theory becomes, owing to asymptotic freedom. As one step in our computation of the $Λ$-parameter in three-flavor QCD, we perform lattice computations in a scheme which allows us to non-perturbatively reach very high energies, corresponding to $α_s = 0.1$ and below. We find that (continuum) perturbation theory is very accurate there, yielding a three percent error in the $Λ$-parameter, while data around $α_s \approx 0.2$ is clearly insufficient to quote such a precision. It is important to realize that these findings are expected to be generic, as our scheme has advantageous properties regarding the applicability of perturbation theory.

preprint2016arXiv

The determination of $α_s$ by the ALPHA collaboration

We review the ALPHA collaboration strategy for obtaining the QCD coupling at high scale. In the three-flavor effective theory it avoids the use of perturbation theory at $α> 0.2$ and at the same time has the physical scales small compared to the cutoff $1/a$ in all stages of the computation. The result $Λ_\overline{MS}^{(3)}=332(14)$~MeV is translated to $α_\overline{MS}(m_Z)=0.1179(10)(2)$ by use of (high order) perturbative relations between the effective theory couplings at the charm and beauty quark "thresholds". The error of this perturbative step is discussed and estimated as $0.0002$.

preprint2015arXiv

A status update on the determination of $Λ_{\overline{\rm MS}}^{N_{\rm f}=3}$ by the ALPHA collaboration

The ALPHA collaboration aims to determine $α_s(m_Z)$ with a total error below the percent level. A further step towards this goal can be taken by combining results from the recent simulations of 2+1-flavour QCD by the CLS initiative with a number of tools developed over the years: renormalized couplings in finite volume schemes, recursive finite size techniques, two-loop renormalized perturbation theory and the (improved) gradient flow on the lattice. We sketch the strategy, which involves both the standard SF coupling in the high energy regime and a gradient flow coupling at low energies. This implies the need for matching both schemes at an intermediate switching scale, $L_{\rm swi}$, which we choose roughly in the range 2-4 GeV. In this contribution we present a preliminary result for this matching procedure, and we then focus on our almost final results for the scale evolution of the SF coupling from $L_{\rm swi}$ towards the perturbative regime, where we extract the $N_{\rm f} = 3$ $Λ$-parameter, $Λ_{\overline{\rm MS}}^{N_{\rm f}=3}$, in units of $L_{\rm swi}$ . Connecting $L_{\rm swi}$ and thus the $Λ$-parameter to a hadronic scale such as $F_K$ requires 2 further ingredients: first, the connection of $L_{\rm swi}$ to $L_{\rm max}$ using a few steps with the step-scaling function of the gradient flow coupling, and, second, the continuum extrapolation of $L_{\rm max} F_K$.

preprint2015arXiv

Determination of c_sw in N_f=3+1 Lattice QCD with massive Wilson fermions

We develop a strategy for the non-perturbative determination of the O(a)-improvement coefficient c_sw for Wilson fermions with massive sea quarks. The improvement condition is defined via the PCAC relation in the Schrödinger functional. It is imposed along a line of constant physics designed to be close to the correct mass of the charm quark. The numerical work uses the tree-level improved Lüscher-Weisz gauge action in N_f=3+1 Lattice QCD.

preprint2015arXiv

Non-perturbative computation of the strong coupling constant on the lattice

We review the long term project of the ALPHA collaboration to compute in QCD the running coupling constant and quark masses at high energy scales in terms of low energy hadronic quantities. The adapted techniques required to numerically carry out the required multiscale non-perturbative calculation with our special emphasis on the control of systematic errors are summarized. The complete results in the two dynamical flavor approximation are reviewed and an outlook is given on the ongoing three flavor extension of the programme with improved target precision.

preprint2015arXiv

Non-perturbative Heavy Quark Effective Theory: Introduction and Status

We give an introduction to Heavy Quark Effective Theory (HQET). Our emphasis is on its formulation non-perturbative in the strong coupling, including the non-perturbative determination of the parameters in the HQET Lagrangian. In a second part we review the present status of HQET on the lattice, largely based on work of the ALPHA collaboration in the last few years. We finally discuss opportunities and challenges.

preprint2015arXiv

Perturbative versus non-perturbative decoupling of heavy quarks

We simulate a theory with $N_f=2$ heavy quarks of mass $M$. At energies much smaller than $M$ the heavy quarks decouple and the theory can be described by an effective theory which is a pure gauge theory to leading order in $1/M$. We present results for the mass dependence of ratios such as $t_0(M)/t_0(0)$. We compute these ratios from simulations and compare them to the perturbative prediction. The latter relies on a factorisation formula for the ratios which is valid to leading order in $1/M$.

preprint2014arXiv

Chiral symmetry breaking in QCD Lite

A distinctive feature of the presence of spontaneous chiral symmetry breaking in QCD is the condensation of low modes of the Dirac operator near the origin. The rate of condensation must be equal to the slope of (Mpi^2 Fpi^2)/2 with respect to the quark mass m in the chiral limit, where Mpi and Fpi are the mass and the decay constant of the Nambu-Goldstone bosons. We compute the spectral density of the (Hermitian) Dirac operator, the quark mass, the pseudoscalar meson mass and decay constant by numerical simulations of lattice QCD with two light degenerate Wilson quarks. We use CLS lattices at three values of the lattice spacing in the range 0.05-0.08 fm, and for several quark masses corresponding to pseudoscalar mesons masses down to 190 MeV. Thanks to this coverage of parameters space, we can extrapolate all quantities to the chiral and continuum limits with confidence. The results show that the low quark modes do condense in the continuum as expected by the Banks-Casher mechanism, and the rate of condensation agrees with the Gell-Mann-Oakes-Renner (GMOR) relation. For the renormalisation-group-invariant ratios we obtain [Σ^RGI]^(1/3)/F =2.77(2)(4) and Lambda^MSbar/F = 3.6(2), which correspond to [Σ^\MSbar(2 GeV)]^(1/3) =263(3)(4) MeV and F=85.8(7)(20) MeV if FK is used to set the scale by supplementing the theory with a quenched strange quark.

preprint2014arXiv

Form factors for $\mathrm B_\mathrm s \to \mathrm K \ell ν$ decays in Lattice QCD

We present the current status of the computation of the form factor $f_+ (q^2)$ for the semi-leptonic decay $\mathrm B_\mathrm s \to \mathrm K \ell ν$ by the ALPHA collaboration. We use gauge configurations which were generated as part of the Coordinated Lattice Simulations (CLS) effort. They have $N_\mathrm f=2$ non-perturbatively $O(a)$ improved Wilson fermions, and pion masses down to $\approx 250 \,\mathrm {MeV}$ with $m_πL \geq 4$. The heavy quark is treated in non-perturbative Heavy Quark Effective Theory (HQET). We discuss how to extract the form factors from the correlation functions and present first results for the form factor at $q^2 = 21.23\,\mathrm{GeV}^2$ extrapolated to the continuum. Next-to-leading order terms in HQET and the chiral extrapolation still need to be included in the analysis.

preprint2014arXiv

On the $N_f$-dependence of gluonic observables

We compute $t_0$, $w_0$ and the topological susceptibility, defined at finite gradient flow time for two-flavour QCD. The use of three lattice spacings and pion masses between 192 and 500 MeV together with a careful error analysis allow to approach the continuum limit of the two-flavour theory despite significant auto-correlations. A comparison to $N_f=0$ results shows the size of sea quark effects in $t_0^2χ$, with $χ$ the topological susceptibility, and low energy observables such as $t_0/w_0^2$ and $t_0/r_0^2$.

preprint2014arXiv

On the effects of heavy sea quarks at low energies

We present a factorisation formula for the dependence of light hadron masses and low energy hadronic scales on the mass $M$ of a heavy quark: apart from an overall factor $Q$, ratios such as $r_0(M)/r_0(0)$ are computable in perturbation theory at large $M$. The mass-independent factor $Q$ is obtained from the theory in the limit $M\to0$ and the decoupled theory with the heavy quark removed. The perturbation theory part is stable concerning different loop orders and our non-perturbative results match on quantitatively to the perturbative prediction. Upon taking ratios of different hadronic scales at the same mass, the perturbative function drops out and the ratios are given by the decoupled theory up to $M^{-2}$ corrections. Our present numerical results are obtained in a model calculation where there are no light quarks and a heavy doublet of quarks is decoupled. They are limited to masses a factor two below the charm. This is not large enough to see the $M^{-2}$ scaling predicted by the theory, but it is sufficient to verify - in the continuum limit - that the sea quark effects of quarks with masses around the charm mass are very small.

preprint2014arXiv

Physical and cut-off effects of heavy sea quarks

We simulate a theory with two dynamical O($a$) improved Wilson quarks whose mass $M$ ranges from a factor eight up to a factor two below the charm quark mass and at three values of the lattice spacing ranging from 0.066 to 0.034 fm. This theory is a prototype to study the decoupling of heavy quarks. We measure the mass and cut-off dependence of ratios of gluonic observables defined from the Wilson flow or the static potential. The size of the 1/$M$ corrections can be determined and disentangled from the lattice artifacts. The difference with the pure gauge theory is at the percent level when two quarks with a mass of the charm quark are present.

preprint2014arXiv

Precision lattice QCD computation of the $B^*Bπ$ coupling

The static $B^{*}Bπ$ coupling, $\hat{g}_χ$, a low energy constant in the leading order heavy meson chiral Lagrangian, is determined using $N_\mathrm{f} = 2$ lattice QCD. We use CLS ensembles with lattice spacings and pion masses down to $a = 0.05\mathrm{fm}$ and $m_π=270\mathrm{MeV}$, and perform combined continuum and chiral extrapolations of our results which have a much better accuracy than previous numbers in the literature. As a by-product, we determine the coupling between the first radial excitations in the $B$ and $B^{*}$ channels ($\hat{g}_{22}$). Accounting for all uncertainties, which are dominated by the chiral extrapolation, we obtain $\hat{g}_χ= 0.492(29)$, while $\hat{g}_{22}$ is somewhat smaller. The comparison to a precise quenched computation suggests that there is little influence by the sea quarks and $\hat{g}_χ$ will not change much when a dynamical strange quark is included.

preprint2014arXiv

Scale setting in lattice QCD

The principles of scale setting in lattice QCD as well as the advantages and disadvantages of various commonly used scales are discussed. After listing criteria for good scales, I concentrate on the main presently used ones with an emphasis on scales derived from the Yang-Mills gradient flow. For these I discuss discretisation errors, statistical precision and mass effects. A short review on numerical results also brings me to an unpleasant disagreement which remains to be explained.

preprint2014arXiv

Spectral density of the Dirac operator in two-flavour QCD

We compute the spectral density of the (Hermitean) Dirac operator in Quantum Chromodynamics with two light degenerate quarks near the origin. We use CLS/ALPHA lattices generated with two flavours of O(a)-improved Wilson fermions corresponding to pseudoscalar meson masses down to 190 MeV, and with spacings in the range 0.05-0.08 fm. Thanks to the coverage of parameter space, we can extrapolate our data to the chiral and continuum limits with confidence. The results show that the spectral density at the origin is non-zero because the low modes of the Dirac operator do condense as expected in the Banks-Casher mechanism. Within errors, the spectral density turns out to be a constant function up to eigenvalues of approximately 80 MeV. Its value agrees with the one extracted from the Gell-Mann-Oakes-Renner relation.

preprint2014arXiv

Topological susceptibility and the sampling of field space in $N_f=2$ lattice QCD simulations

We present a measurement of the topological susceptibility in two flavor QCD. In this observable, large autocorrelations are present and also sizable cutoff effects have to be faced in the continuum extrapolation. Within the statistical accuracy of the computation, the result agrees with the expectation from leading order chiral perturbation theory.

preprint2013arXiv

A one-loop study of matching conditions for static-light flavor currents

Heavy Quark Effective Theory (HQET) computations of semi-leptonic decays, e.g. B -> pi l nu, require the knowledge of the parameters in the effective theory for all components of the heavy-light flavor currents. So far non-perturbative matching conditions have been employed only for the time component of the axial current. Here we perform a check of matching conditions for the time component of the vector current and the spatial component of the axial vector current up to one-loop order of perturbation theory and to lowest order of the 1/m-expansion. We find that the proposed observables have small higher order terms in the 1/m-series and are thus excellent candidates for a non-perturbative matching procedure.

preprint2013arXiv

Chiral condensate from the Banks-Casher relation

We report on our ongoing project of determining the chiral condensate of two-flavor QCD from the Banks-Casher relation. We compute the mode number of the O(a)-improved Wilson-Dirac operator for several values of Λ, and we discuss different fitting strategies to extract the chiral condensate from its mass and Λdependence. Our preliminary results haven been obtained at two different lattice spacings by using CLS-configurations.

preprint2012arXiv

B $\to$ $π$ form factor with 2 flavours of $O(a)$ improved Wilson quarks

The determinations of $|V_{\rm ub}|$ from the exclusive branching ratios of $B\to τν$ and $B \to πl ν$ tend to show a tension at the level of $3σ$ \cite{Beringer:1900zz}. On the theoretical side they depend on the lattice computation of the hadronic matrix elements $f_{\rm B}$ and the $B\to π$ form factor $f_+(q^2)$. To understand the tension, improved precision and a careful analysis of the systematics involved are necessary. Working towards this goal, we present preliminary lattice results of the ALPHA collaboration for the $B\to π$ form factor $f_+(q^2)$ with $N_{\rm f}=2$ flavours of $O(a)$-improved Wilson fermions. Our computation uses HQET in the static limit, pion masses ranging down to $\sim250$ MeV, large volumes with $m_πL >4$, three lattice spacings, and non-perturbative renormalization. We describe the techniques adopted to reduce the statistical noise (stochastic all-to-all with full time dilution) and the contamination from excited states (smearing for the B and the pion). We estimate the size of the chiral and continuum extrapolations. We discuss the impact our result could have to clarify the above mentioned discrepancy in the determination of $|V_{\rm ub}|$.

preprint2012arXiv

Parameters of Heavy Quark Effective Theory from Nf=2 lattice QCD

We report on a non-perturbative determination of the parameters of the lattice Heavy Quark Effective Theory (HQET) Lagrangian and of the time component of the heavy-light axial-vector current with Nf=2 flavors of massless dynamical quarks. The effective theory is considered at the 1/mh order, and the heavy mass mh covers a range from slightly above the charm to beyond the beauty region. These HQET parameters are needed to compute, for example, the b-quark mass, the heavy-light spectrum and decay constants in the static approximation and to order 1/mh in HQET. The determination of the parameters is done non-perturbatively. The computation reported in this paper uses the plaquette gauge action and two different static actions for the heavy quark described by HQET. For the light-quark action we choose non-perturbatively O(a)-improved Wilson fermions.

preprint2012arXiv

The strange quark mass and Lambda parameter of two flavor QCD

We complete the non-perturbative calculations of the strange quark mass and the Lambda parameter in two flavor QCD by the ALPHA collaboration. The missing lattice scale is determined via the kaon decay constant, for whose chiral extrapolation complementary strategies are compared. We also give a value for the scale r_0 in physical units as well as an improved determination of the renormalization constant Z_A.

preprint2011arXiv

Determination of the Static Potential with Dynamical Fermions

We present in detail a technique to extract the potential between a static quark and anti-quark pair from Wilson loops measured on dynamical configurations. This technique is based on HYP smearing and leads to an exponential improvement of the noise-to-signal ratio of Wilson loops. We explain why the correct continuum potential is obtained and show numerical evidence that the cut-off effects are small. We present precise results for the non-perturbative potential. As applications, we determine the scale r_0/a and study the shape of the static potential in the range of distances around r_0, where it can be compared with phenomenological potential models.

preprint2011arXiv

Determination of the Wilson ChPT low energy constant c_2

Following a suggestion by Aoki and Bär, c2 can be extracted by analyzing volume effects in 2-pion states. To this end we consider renormalized ratios of four point to two point correlation functions. We present the results from various CLS lattices, with pion masses ranging from 280 to 450 MeV and lattice spacings of 0.07 fm and 0.08 fm. This low energy constant is useful to understand discretization effects in chiral extrapolations with Wilson fermions, especially for quantities which vanish in the chiral limit, like the pion mass. Since our procedure is computationally cheap and straightforward, it is recommended as a routine study for any Wilson-type simulation, as a check of discretization effects.

preprint2011arXiv

M_b and f_B from non-perturbatively renormalized HQET with Nf=2 light quarks

We present an updated analysis of the non-perturbatively renormalized b-quark mass and B meson decay constant based on CLS lattices with two dynamical non-perturbatively improved Wilson quarks. This update incorporates additional light quark masses and lattice spacings in large physical volume to improve chiral extrapolations and to reach the continuum limit. We use Heavy Quark Effective Theory (HQET) including 1/m_b terms with non-perturbative coefficients based on the matching of QCD and HQET developed by the ALPHA collaboration during the past years.

preprint2011arXiv

On the computation of hadron-to-hadron transition matrix elements in lattice QCD

We discuss the accurate determination of matrix elements < f| h_w | i > where neither |i> nor |f> is the vacuum state and h_w is some operator. Using solutions of the Generalized Eigenvalue Problem (GEVP) we construct estimators for matrix elements which converge rapidly as a function of the Euclidean time separations involved. |i> and |f> may be either the ground state in a given hadron channel or an excited state. Apart from a model calculation, the estimators are demonstrated to work well for the computation of the B*B pi-coupling in the quenched approximation. They are also compared to a standard ratio as well as to the "summed ratio method" of [1,2,3]. In the model, we also illustrate the ordinary use of the GEVP for energy levels.

preprint2011arXiv

One-loop lattice artifacts of a dynamical charm quark

For a few observables in O(a) improved lattice QCD, we compute discretization effects arising from the vacuum polarization of a heavy quark at one-loop order. In particular, the force between static quarks, the running coupling in the Schrodinger functional and a related quantity, {\bar \upsilon}, are considered. Results show that the cutoff effects of a dynamical charm quark are typically smaller than those present in the pure gauge theory. This perturbative result is a good indication that dynamical charm quarks are feasible already now.

preprint2011arXiv

Strange quark mass and Lambda parameter by the ALPHA collaboration

We determine f_K for lattice QCD in the two flavor approximation with non-perturbatively improved Wilson fermions. The result is used to set the scale for dimensionful quantities in CLS/ALPHA simulations. To control its dependence on the light quark mass, two different strategies for the chiral extrapolation are applied. Combining f_K and the bare strange quark mass with non-perturbative renormalization factors and step scaling functions computed in the Schroedinger Functional, we determine the RGI strange quark mass and the Lambda parameter in units of f_K.

preprint2010arXiv

B meson spectrum and decay constant from Nf=2 simulations

We report on the status of an ALPHA Collaboration project to extract quantities for B physics phenomenology from Nf=2 lattice simulations. The framework is Heavy Quark Effective Theory (HQET) expanded up to the first order of the inverse b-quark mass. The couplings of the effective theory are determined by imposing matching conditions of observables computed in HQET with their counterpart computed in QCD. That program, based on Nf=2 simulations in a small physical volume with SF boundary conditions, is now almost finished. On the other side the analysis of configurations selected from the CLS ensembles, in order to measure HQET hadronic matrix elements, has just started recently so that only results obtained at a single lattice spacing, a=0.07 fm, will be discussed. We give our first results for the b-quark mass and for the B meson decay constant.

preprint2010arXiv

Critical slowing down and error analysis in lattice QCD simulations

We study the critical slowing down towards the continuum limit of lattice QCD simulations with Hybrid Monte Carlo type algorithms. In particular for the squared topological charge we find it to be very severe with an effective dynamical critical exponent of about 5 in pure gauge theory. We also consider Wilson loops which we can demonstrate to decouple from the modes which slow down the topological charge. Quenched observables are studied and a comparison to simulations of full QCD is made. In order to deal with the slow modes in the simulation, we propose a method to incorporate the information from slow observables into the error analysis of physical observables and arrive at safer error estimates.

preprint2010arXiv

HQET at order $1/m$: II. Spectroscopy in the quenched approximation

Using Heavy Quark Effective Theory with non-perturbatively determined parameters in a quenched lattice calculation, we evaluate the splittings between the ground state and the first two radially excited states of the $B_s$ system at static order. We also determine the splitting between first excited and ground state, and between the $B_s^*$ and $B_s$ ground states to order $1/m_b$. The Generalized Eigenvalue Problem and the use of all-to-all propagators are important ingredients of our approach.

preprint2010arXiv

HQET at order 1/m: I. Non-perturbative parameters in the quenched approximation

We determine non-perturbatively the parameters of the lattice HQET Lagrangian and those of the time component of the heavy-light axial-vector current in the quenched approximation. The HQET expansion includes terms of order 1/m. Our results allow to compute, for example, the heavy-light spectrum and B-meson decay constants in the static approximation and to order 1/m in HQET. The determination of the parameters is separated into universal and non-universal parts. The universal results can be used to determine the parameters for various discretizations. The computation reported in this paper uses the plaquette gauge action and the "HYP1/2" action for the b-quark described by HQET. The parameters of the current also depend on the light-quark action, for which we choose non-perturbatively O(a)-improved Wilson fermions.

preprint2010arXiv

HQET at order 1/m: III. Decay constants in the quenched approximation

We report on the computation of the $B_s$ meson decay constant in Heavy Quark Effective Theory on the lattice. The next to leading order corrections in the HQET expansion are included non-perturbatively. We estimate higher order contributions to be very small. The results are extrapolated to the continuum limit, the main systematic error affecting the computation is therefore the quenched approximation used here. The Generalized Eigenvalue Problem and the use of all-to-all propagators are important technical ingredients of our approach that allow to keep statistical and systematic errors under control. We also report on the decay constant $f_{B'_s}$ of the first radially excited state in the $B_s$ sector, computed in the static limit.

preprint2010arXiv

Running of the SF-coupling with four massless flavours

We discuss the status of different determinations of alpha_s, motivating a precise and reliable computation from lattice QCD. In order to suppress perturbative errors, the non-perturbative computation has to reach high energy scales mu. Such results already exist in the SF-scheme for N_f=0,2 and N_f=3. We recently added the running with four massless flavours in a range of alpha from about 0.07 to 0.3 . It is based on our recent determination of the Sheikholeslami Wohlert coefficient in the four-flavour theory.

preprint2010arXiv

The B*Bpi Coupling in the Static Limit

We study an accurate method for the lattice calculation of the B*Bpi Coupling in the static limit, paying particular attention to excited state contamination. As this coupling is a parameter of the heavy meson chiral Lagrangian, it is useful for constraining the chiral behaviour of various observables in B-physics. We present a precise study of the continuum limit in the quenched approximation and preliminary results with 2 flavours of improved Wilson quarks (using CLS lattices) for pion masses down to around 250 MeV. With dynamical quarks both the lattice spacing and the light quark mass dependences are found to be very weak. We can quote g = 0.51(2) for the continuum value in the chiral limit, where the error will be reduced when a more complete analysis has been performed.

preprint2009arXiv

Spectroscopy and Decay Constants from Nonperturbative HQET at Order 1/m

We carry out a thorough analysis with the GEVP method to obtain ground-state and first-excited-state masses and decay constants of bottom-strange (pseudo-scalar and vector) mesons. This computation is done for quenched, nonperturbatively renormalized HQET, including order $1/m_b$ terms. The continuum limit is obtained using three lattice spacings and two static actions.

preprint2007arXiv

Precision for B-meson matrix elements

We demonstrate how HQET and the Step Scaling Method for B-physics, pioneered by the Tor Vergata group, can be combined to reach a further improved precision. The observables considered are the mass of the b-quark and the B_s-meson decay constant. The demonstration is carried out in quenched lattice QCD. We start from a small volume, where one can use a standard O(a)-improved relativistic action for the b-quark, and compute two step scaling functions which relate the observables to the large volume ones. In all steps we extrapolate to the continuum limit, separately in HQET and in QCD for masses below m_b. The physical point m_b is then reached by an interpolation of the continuum results in 1/m. The essential, expected and verified, feature is that the step scaling fuctions have a weak mass-dependence resulting in an easy interpolation to the physical point. With r_0=0.5fm and the experimental B_s and K masses as input, we find F_Bs=191(6)MeV and the renormalization group invariant mass M_b=6.88(10)GeV, translating into m_b(m_b)=4.42(6)GeV in the MSbar scheme. This approach seems very promising for full QCD.

preprint2006arXiv

m_b and f_Bs from a combination of HQET and QCD

We compute the mass of the b-quark and the B_s meson decay constant in quenched lattice QCD using a combination of HQET and the standard relativistic QCD Lagrangian. We start from a small volume, where one can directly deal with the b-quark, and compute the evolution to a big volume, where the finite size effects are negligible through step scaling functions which give the change of the observables when L is changed to 2L. In all steps we extrapolate to the continuum limit, separately in HQET and in QCD for masses below m_b. The point m_b is then reached by an interpolation of the continuum results. With r_0=0.5 fm and the experimental B_s and K masses we find f_Bs=191(6) MeV and the renormalization group invariant mass M_b=6.89(11) GeV, translating into m_b(m_b)=4.42(7) GeV in the MS-bar scheme.

preprint2004arXiv

Computation of the strong coupling in QCD with two dynamical flavours

We present a non-perturbative computation of the running of the coupling alpha_s in QCD with two flavours of dynamical fermions in the Schroedinger functional scheme. We improve our previous results by a reliable continuum extrapolation. The Lambda-parameter characterizing the high-energy running is related to the value of the coupling at low energy in the continuum limit. An estimate of Lambda*r_0 is given using large-volume data with lattice spacings a from 0.07 fm to 0.1 fm. It translates into Lambda_{MSbar}^{(2)}=245(16)(16) MeV [assuming r_0=0.5 fm]. The last step still has to be improved to reduce the uncertainty.

preprint2004arXiv

Non-perturbative Heavy Quark Effective Theory

We explain how to perform non-perturbative computations in HQET on the lattice. In particular the problem of the subtraction of power-law divergences is solved by a non-perturbative matching of HQET and QCD. As examples, we present a full calculation of the mass of the b-quark in the combined static and quenched approximation and outline an alternative way to obtain the B-meson decay constant at lowest order. Since no excessively large lattices are required, our strategy can also be applied including dynamical fermions.

preprint2002arXiv

Non-perturbative renormalization of HQET and QCD

We discuss the necessity of non-perturbative renormalization in QCD and HQET and explain the general strategy for solving this problem. A few selected topics are discussed in some detail, namely the importance of off-shell improvement in the MOM-scheme on the lattice, recent progress in the implementation of finite volume schemes and then particular emphasis is put on the recent idea to carry out a non-perturbative renormalization of the Heavy Quark Effective Theory.

preprint2001arXiv

A strategy to compute the b-quark mass with non-perturbative accuracy

We describe a strategy for a non-perturbative computation of the b-quark mass to leading order in 1/m in the Heavy Quark Effective Theory (HQET). The approach avoids the perturbative subtraction of power law divergencies, and the continuum limit may be taken. First numerical results in the quenched approximation demonstrate the potential of the method with a preliminary result m_b(4GeV)=4.56(2)(7) GeV. In principle, the idea may also be applied to the matching of composite operators or the computation of 1/m corrections in HQET.

preprint2001arXiv

Heavy Quark Effective Theory at one-loop order: An explicit example

We consider correlation functions containing the axial current of one light and one heavy quark in the static approximation as well as in full QCD, using the lattice regularization. Up to one-loop order of perturbation theory, we study the difference between the full and the effective theory in the continuum limit. In the full theory we find a term non-analytic in 1/m, revealing the asymptotic character of the 1/m-expansion. In general, deviations from the m-to-infinity limit turn out to be small and are well described by the first non-trivial terms when m is a factor 2-3 above the external scale. We also investigate the mass dependence of discretization errors, and find that the behaviour of the correlation functions at finite lattice spacing differs significantly from that in the continuum limit when the quark mass is large.

preprint2001arXiv

The Nf=0 heavy quark potential from short to intermediate distances

We study the potential of a static quark anti-quark pair in the range 0.05fm \leq r \leq 0.8fm, employing a sequence of lattices up to 64^4. Lattice artifacts in potential and force are investigated theoretically as well as numerically and continuum quantities are obtained by extrapolation of the results at finite lattice spacing. Consistency of the numerical results with the form of scaling violations predicted by an analysis `a la Symanzik is found. The scale r_0/a is determined for the Wilson action up to beta=6.92.

preprint2000arXiv

Effective Chiral Lagrangians and Lattice QCD

We propose a general method to obtain accurate estimates for some of the "low-energy constants" in the one-loop effective chiral Lagrangian by means of simulating lattice QCD. In particular, the method is sensitive to those constants whose values are required to test the hypothesis of a massless up-quark. Initial tests performed in the quenched approximation confirm that good statistical precision can be achieved. As a byproduct we obtain an accurate estimate for the ratio of pseudoscalar decay constants, F_K/F_pi, in the quenched approximation, which lies 10% below the experimental result. The quantities that serve to extract the low-energy constants also allow a test of the scaling behaviour of different discretizations of QCD and a search for the effects of dynamical quarks.

preprint2000arXiv

Renormalization and O(a)-improvement of the static-light axial current

A systematic treatment of O(a)-improvement in lattice theories with static quarks is presented. The Schrödinger functional is discussed and a renormalization condition for the static axial current in the SF-scheme is introduced. Its relation to other schemes is computed to 1-loop order and the 2-loop anomalous dimension is derived. In finite volume renormalization schemes such as the SF-scheme, the renormalization scale dependence of the renormalized quantities is described by the step scaling function which can be computed by MC- simulations. We evaluate its lattice spacing effects in perturbation theory.

preprint1999arXiv

Hadron masses and matrix elements from the QCD Schr"odinger functional

We explain how masses and matrix elements can be computed in lattice QCD using Schr"odinger functional boundary conditions. Numerical results in the quenched approximation demonstrate that good precision can be achieved. For a statistical sample of the same size, our hadron masses have a precision similar to what is achieved with standard methods, but for the computation of matrix elements such as the pseudoscalar decay constant the Schr"odinger functional technique turns out to be much more efficient than the known alternatives.

preprint1996arXiv

Non-perturbative determination of the axial current normalization constant in O(a) improved lattice QCD

A finite-size technique is employed to compute the normalization constant $Z_A$ of the isovector axial current in lattice QCD. The calculation is carried out in the quenched approximation for values of the bare gauge coupling $g_0$ ranging from 0 to 1. In the lattice action and the lattice expression for the axial current we include the counterterms required for O(a) improvement, with non-perturbatively determined coefficients. With little additional work the normalization constant $Z_V$ of the improved isospin current is also obtained.

preprint1996arXiv

Non-perturbative O(a) improvement of lattice QCD

The coefficients multiplying the counterterms required for O($a$) improvement of the action and the isovector axial current in lattice QCD are computed non-perturbatively, in the quenched approximation and for bare gauge couplings $g_0$ in the range $0 \leq g_0 \leq 1$. A finite-size method based on the Schrödinger functional is employed, which enables us to perform all calculations at zero or nearly zero quark mass. As a by-product the critical hopping parameter $κ_c$ is obtained at all couplings considered.

preprint1992arXiv

Computation of the Running Coupling in the SU(2) Yang-Mills Theory

A finite-size scaling technique is applied to the SU(2) gauge theory (without matter fields) to compute a non-perturbatively defined running coupling alpha(q) for a range of momenta q given in units of the string tension K. We find that already at rather low q, the evolution of alpha(q) is well described by the 2-loop approximation to the Callan-Symanzik beta-function. At the highest momentum reached, q=20 sqrt(K), we obtain alpha_MSbar(q)=0.187 +/- 0.005 +/- 0.009 for the running coupling in the MSbar scheme of dimensional regularization.