Source author record

Christoph Lehner

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

Catalog footprint

Research overview

26works
4topics
4close collaborators

Published work

26 publications

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 bb and cc 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

Discovering new physics in rare kaon decays

The decays and mixing of KK mesons are remarkably sensitive to the weak interactions of quarks and leptons at high energies. They provide important tests of the standard model at both first and second order in the Fermi constant GFG_F and offer a window into possible new phenomena at energies as high as 1,000 TeV. These possibilities become even more compelling as the growing capabilities of lattice QCD make high-precision standard model predictions possible. Here we discuss and attempt to forecast some of these capabilities.

preprint2022arXiv

Grid: OneCode and FourAPIs

We discuss a substantial update to the Grid software library for Lattice QCD, enabling it to port to multiple GPU architectures while retaining CPU vectorisation and SIMD execution within OpenMP threads. The GPU environments supported include vendor specific Nvidia CUDA and AMD HIP environments and a (mostly) standards based SYCL implementation. This is performed by an internal abstraction interface giving single source cross-platform performance portability across all number of planned Exascale architectures, and all those planned by the US Department of Energy.

preprint2020arXiv

Consistency of hadronic vacuum polarization between lattice QCD and the R-ratio

There are emerging tensions for theory results of the hadronic vacuum polarization contribution to the muon anomalous magnetic moment both within recent lattice QCD calculations and between some lattice QCD calculations and R-ratio results. In this paper we work towards scrutinizing critical aspects of these calculations. We focus in particular on a precise calculation of Euclidean position-space windows defined by RBC/UKQCD that are ideal quantities for cross-checks within the lattice community and with R-ratio results. We perform a lattice QCD calculation using physical up, down, strange, and charm sea quark gauge ensembles generated in the staggered formalism by the MILC collaboration. We study the continuum limit using inverse lattice spacings from a−1≈1.6a^{-1}\approx 1.6 GeV to 3.53.5 GeV, identical to recent studies by FNAL/HPQCD/MILC and Aubin et al. and similar to the recent study of BMW. Our calculation exhibits a tension for the particularly interesting window result of aμud,conn.,isospin,Wa_μ^{\rm ud, conn.,isospin, W} from 0.40.4 fm to 1.01.0 fm with previous results obtained with a different discretization of the vector current on the same gauge configurations. Our results may indicate a difficulty related to estimating uncertainties of the continuum extrapolation that deserves further attention. In this work we also provide results for aμud,conn.,isospina_μ^{\rm ud,conn.,isospin}, aμs,conn.,isospina_μ^{\rm s,conn.,isospin}, aμSIB,conn.a_μ^{\rm SIB,conn.} for the total contribution and a large set of windows. For the total contribution, we find aμHVP LO=714(27)(13)10−10a_μ^{\rm HVP~LO}=714(27)(13) 10^{-10}, aμud,conn.,isospin=657(26)(12)10−10a_μ^{\rm ud,conn.,isospin}=657(26)(12) 10^{-10}, aμs,conn.,isospin=52.83(22)(65)10−10a_μ^{\rm s,conn.,isospin}=52.83(22)(65) 10^{-10}, and aμSIB,conn.=9.0(0.8)(1.2)10−10a_μ^{\rm SIB,conn.}=9.0(0.8)(1.2) 10^{-10}, where the first uncertainty is statistical and the second systematic. We also comment on finite-volume corrections for the strong-isospin-breaking corrections.

preprint2020arXiv

The hadronic light-by-light scattering contribution to the muon anomalous magnetic moment from lattice QCD

We report the first result for the hadronic light-by-light scattering contribution to the muon anomalous magnetic moment with all errors systematically controlled. Several ensembles using 2+1 flavors of physical mass Möbius domain-wall fermions, generated by the RBC/UKQCD collaborations, are employed to take the continuum and infinite volume limits of finite volume lattice QED+QCD. We find aμHLbL=7.87(3.06)stat(1.77)sys×10−10a_μ^{\rm HLbL} = 7.87(3.06)_\text{stat}(1.77)_\text{sys}\times 10^{-10}. Our value is consistent with previous model results and leaves little room for this notoriously difficult hadronic contribution to explain the difference between the Standard Model and the BNL experiment.

preprint2016arXiv

Emerging lattice approach to the K-Unitarity Triangle

It has been clear for past many years that in low energy observables new physics can only appear as a perturbation. Therefore precise theoretical predictions and precise experimental measurements have become mandatory. Here we draw attention to the significant advances that have been made on the lattice in recent years in K→ππK\to ππ, ΔMKΔM_K, the long-distance part of ε\varepsilon and rare K-decays. Thus, in conjunction with experiments, the construction of a unitarity triangle purely from Kaon physics should soon become feasible. Along with the B-unitarity triangle, this should allow for more stringent tests of the Standard Model and tighter constraints on new physics.

preprint2016arXiv

Form factors for semi-leptonic BB decays

Semi-leptonic BB decays provide promising channels to test the Standard Model, search for signs of new physics, or determine fundamental parameters like CKM matrix elements. We present an update on our calculation of short distance contributions to GIM suppressed rare BB decays focusing in particular on Bs→ϕℓ+ℓ−B_s\to ϕ\ell^+ \ell^- decays. Furthermore we show first results for our calculation of B(s)→D(s)(∗)ℓνB_{(s)}\to D_{(s)}^{(*)}\ellν semi-leptonic decays involving b→cb\to c transitions. Our calculations are based on RBC-UKQCD's 2+1 flavor domain-wall fermion and Iwasaki gauge field configurations featuring three lattice spacings in the range 1.731.73 GeV ≤a−1≤2.77\le a^{-1} \le 2.77 GeV and pion masses down to the physical value. We calculate the form factors by simulating bb-quarks using the relativistic heavy quark action, create light u/du/d and ss quarks with standard domain-wall kernel, and use optimised Möbius domain-wall fermions for charm quarks.

preprint2016arXiv

The connected and leading disconnected diagrams of the hadronic light-by-light contribution to muon g−2g - 2

We report our recent lattice calculation of hadronic light-by-light contribution to muon g−2g-2 using our recently developed moment method. The connected diagrams and the leading disconnected diagrams are included. The calculation is performed on a 483×9648^3 \times 96 lattice with physical pion mass and 5.5 fm box size. We expect sizable finite volume and finite lattice spacing corrections to the results of these calculations which will be estimated in calculations to be carried out over the next 1-2 years.

preprint2015arXiv

Λb→pℓ−νˉℓΛ_b \to p \ell^- \barν_\ell and Λb→Λcℓ−νˉℓΛ_b \to Λ_c \ell^- \barν_\ell form factors from lattice QCD with relativistic heavy quarks

Measurements of the Λb→pℓ−νˉℓΛ_b \to p \ell^- \barν_\ell and Λb→Λcℓ−νˉℓΛ_b \to Λ_c \ell^- \barν_\ell decay rates can be used to determine the magnitudes of the CKM matrix elements VubV_{ub} and VcbV_{cb}, provided that the relevant hadronic form factors are known. Here we present a precise calculation of these form factors using lattice QCD with 2+1 flavors of dynamical domain-wall fermions. The bb and cc quarks are implemented with relativistic heavy-quark actions, allowing us to work directly at the physical heavy-quark masses. The lattice computation is performed for six different pion masses and two different lattice spacings, using gauge-field configurations generated by the RBC and UKQCD collaborations. The b→ub \to u and b→cb \to c currents are renormalized with a mostly nonperturbative method. We extrapolate the form factor results to the physical pion mass and the continuum limit, parametrizing the q2q^2-dependence using zz-expansions. The form factors are presented in such a way as to enable the correlated propagation of both statistical and systematic uncertainties into derived quantities such as differential decay rates and asymmetries. Using these form factors, we present predictions for the Λb→pℓ−νˉℓΛ_b \to p \ell^- \barν_\ell and Λb→Λcℓ−νˉℓΛ_b \to Λ_c \ell^- \barν_\ell differential and integrated decay rates. Combined with experimental data, our results enable determinations of ∣Vub∣|V_{ub}|, ∣Vcb∣|V_{cb}|, and ∣Vub/Vcb∣|V_{ub}/V_{cb}| with theory uncertainties of 4.4%, 2.2%, and 4.9%, respectively.

preprint2015arXiv

B-meson decay constants from 2+1-flavor lattice QCD with domain-wall light quarks and relativistic heavy quarks

We calculate the B-meson decay constants f_B, f_Bs, and their ratio in unquenched lattice QCD using domain-wall light quarks and relativistic b-quarks. We use gauge-field ensembles generated by the RBC and UKQCD collaborations using the domain-wall fermion action and Iwasaki gauge action with three flavors of light dynamical quarks. We analyze data at two lattice spacings of a ~ 0.11, 0.086 fm with unitary pion masses as light as M_pi ~ 290 MeV; this enables us to control the extrapolation to the physical light-quark masses and continuum. For the b-quarks we use the anisotropic clover action with the relativistic heavy-quark interpretation, such that discretization errors from the heavy-quark action are of the same size as from the light-quark sector. We renormalize the lattice heavy-light axial-vector current using a mostly nonperturbative method in which we compute the bulk of the matching factor nonperturbatively, with a small correction, that is close to unity, in lattice perturbation theory. We also improve the lattice heavy-light current through O(alpha_s a). We extrapolate our results to the physical light-quark masses and continuum using SU(2) heavy-meson chiral perturbation theory, and provide a complete systematic error budget. We obtain f_B0 = 199.5(12.6) MeV, f_B+ = 195.6(14.9) MeV, f_Bs = 235.4(12.2) MeV, f_Bs/f_B0 = 1.197(50), and f_Bs/f_B+ = 1.223(71), where the errors are statistical and total systematic added in quadrature. These results are in good agreement with other published results and provide an important independent cross check of other three-flavor determinations of BB-meson decay constants using staggered light quarks.

preprint2015arXiv

Covariant approximation averaging

We present a new class of statistical error reduction techniques for Monte-Carlo simulations. Using covariant symmetries, we show that correlation functions can be constructed from inexpensive approximations without introducing any systematic bias in the final result. We introduce a new class of covariant approximation averaging techniques, known as all-mode averaging (AMA), in which the approximation takes account of contributions of all eigenmodes through the inverse of the Dirac operator computed from the conjugate gradient method with a relaxed stopping condition. In this paper we compare the performance and computational cost of our new method with traditional methods using correlation functions and masses of the pion, nucleon, and vector meson in Nf=2+1N_f=2+1 lattice QCD using domain-wall fermions. This comparison indicates that AMA significantly reduces statistical errors in Monte-Carlo calculations over conventional methods for the same cost.

preprint2015arXiv

Hadronic Light by Light Contributions to the Muon Anomalous Magnetic Moment With Physical Pions

The current measurement of muonic g−2g - 2 disagrees with the theoretical calculation by about 3 standard deviations. Hadronic vacuum polarization (HVP) and hadronic light by light (HLbL) are the two types of processes that contribute most to the theoretical uncertainty. The current value for HLbL is still given by models. I will describe results from a first-principles lattice calculation with a 139 MeV pion in a box of 5.5 fm extent. Our current numerical strategies, including noise reduction techniques, evaluating the HLbL amplitude at zero external momentum transfer, and important remaining challenges, in particular those associated with finite volume effects, will be discussed.

preprint2015arXiv

Lattice Calculation of Hadronic Light-by-Light Contribution to the Muon Anomalous Magnetic Moment

The quark-connected part of the hadronic light-by-light scattering contribution to the muon's anomalous magnetic moment is computed using lattice QCD with chiral fermions. We report several significant algorithmic improvements and demonstrate their effectiveness through specific calculations which show a reduction in statistical errors by more than an order of magnitude. The most realistic of these calculations is performed with a near-physical, 171171 MeV pion mass on a (4.6  fm)3(4.6\;\mathrm{fm})^3 spatial volume using the 323×6432^3\times 64 Iwasaki+DSDR gauge ensemble of the RBC/UKQCD Collaboration.

preprint2015arXiv

Lattice Calculation of the Connected Hadronic Light-by-Light Contribution to the Muon Anomalous Magnetic Moment

The anomalous magnetic moment of muon, g−2g-2, is a very precisely measured quantity. However, the current measurement disagrees with standard model by about 3 standard deviations. Hadronic vacuum polarization and hadronic light by light are the two types of processes that contribute most to the theoretical uncertainty. I will describe how lattice methods are well-suited to provide a first-principle's result for the hadronic light by light contribution, the various numerical strategies that are presently being used to evaluate it, our current results and the important remaining challenges which must be overcome.

preprint2015arXiv

Neutral BB meson mixings and BB meson decay constants with static heavy and domain-wall light quarks

Neutral BB meson mixing matrix elements and BB meson decay constants are calculated. Static approximation is used for bb quark and domain-wall fermion formalism is employed for light quarks. The calculations are carried out on 2+12+1 flavor dynamical ensembles generated by RBC/UKQCD Collaborations with lattice spacings 0.0860.086fm (a−1∼2.3a^{-1}\sim 2.3GeV) and 0.110.11fm (1.71.7GeV), and a fixed physical spatial volume of about (2.7fm)3(2.7{\rm fm})^3. In the static quark action, link-smearings are used to improve the signal-to-noise ratio. We employ two kinds of link-smearings, HYP1 and HYP2, and their results are combined in taking the continuum limit. For the matching between the lattice and the continuum theory, one-loop perturbative O(a)O(a) improvements are made to reduce discretization errors. As the most important quantity of this work, we obtain SU(3) breaking ratio ξ=1.208(60)ξ=1.208(60), where the error includes statistical and systematic one. (Uncertainty from infinite bb quark mass is not included.) We also find other neutral BB meson mixing quantities fBB^B=240(22)f_B\sqrt{\hat{B}_B}=240(22)MeV, fBsB^Bs=290(22)f_{B_s}\sqrt{\hat{B}_{B_s}}=290(22)MeV, B^B=1.17(22)\hat{B}_B=1.17(22), B^Bs=1.22(13)\hat{B}_{B_s}=1.22(13) and BBs/BB=1.028(74)B_{B_s}/B_B=1.028(74), BB meson decay constants fB=219(17)f_B=219(17)MeV, fBs=264(19)f_{B_s}=264(19)MeV and fBs/fB=1.193(41)f_{B_s}/f_B=1.193(41), in the static limit of bb quark, which do not include infinite bb quark mass uncertainty.

preprint2015arXiv

On calculating disconnected-type hadronic light-by-light scattering diagrams from lattice QCD

For reliable comparison of the standard model prediction to the muon g-2 with its experimental value, the hadronic light-by-light scattering (HLbL) contribution must be calculated by lattice QCD simulation. HLbL contribution has many types of disconnected-type diagrams. Here, we start with recalling the point that must be taken care of in every method to calculate them by lattice QCD, and present one concrete method called nonperturbative QED method.

preprint2013arXiv

Neutral B meson mixing with static heavy and domain-wall light quarks

Neutral B meson mixing matrix elements and B meson decay constants are calculated. Static approximation is used for b quark and domain-wall fermion formalism is employed for light quarks. The calculations are done on 2+1 flavor dynamical ensembles, whose lattice spacings are 0.086 fm and 0.11 fm with a fixed physical spatial volume of about (2.7 fm)^3. In the static quark action, link-smearings are used to improve the signal-to-noise ratio. We employ two kinds of link-smearings and their results are combined in taking a continuum limit. For the matching between the lattice and the continuum theory, one-loop perturbative calculations are used including O(a) improvements to reduce discretization errors. We obtain SU(3) braking ratio ξ=1.222(60) in the static limit of b quark.

preprint2012arXiv

Automated lattice perturbation theory and relativistic heavy quarks in the Columbia formulation

We introduce a new computer algebra system optimized for use in lattice perturbation theory as well as continuum perturbation theory and a new framework to perform automated perturbative calculations on top of said computer algebra system. The new framework is used to tune the relativistic heavy quark action in the Columbia formulation at one loop in meanfield-improved perturbation theory. Preliminary results for the matching and O(a)-improvement of heavy-light axial vector currents with light domain-wall quarks are also presented.

preprint2012arXiv

Nonperturbative tuning of an improved relativistic heavy-quark action with application to bottom spectroscopy

We calculate the masses of bottom mesons using an improved relativistic action for the b-quarks and the RBC/UKQCD Iwasaki gauge configurations with 2+1 flavors of dynamical domain-wall light quarks. We analyze configurations with two lattice spacings: a^{-1} = 1.729 GeV (a ~ 0.11 fm) and a^{-1} = 2.281 GeV (a ~ 0.086 fm). We use an anisotropic, clover-improved Wilson action for the b-quark, and tune the three parameters of the action nonperturbatively such that they reproduce the experimental values of the B_s and B_s* heavy-light meson states. The masses and mass-splittings of the low-lying bottomonium states (such as the eta_b and Upsilon) can then be computed with no additional inputs, and comparison between these predictions and experiment provides a test of the validity of our method. We obtain bottomonium masses with total uncertainties of ~0.5-0.6% and fine-structure splittings with uncertainties of ~35-45%; for all cases we find good agreement with experiment. The parameters of the relativistic heavy-quark action tuned for b-quarks presented in this work can be used for precise calculations of weak matrix elements such as B-meson decay constants and mixing parameters with lattice discretization errors that are of the same size as in light pseudoscalar meson quantities. This general method can also be used for charmed meson masses and matrix elements if the parameters of the heavy-quark action are appropriately tuned.

preprint2011arXiv

Geometry dependence of RMT-based methods to extract the low-energy constants Sigma and F

The lowest-order low-energy constants ΣΣ and FF of chiral pertubation theory can be extracted from lattice data using methods based on the equivalence of random matrix theory (RMT) and QCD in the epsilon regime. We discuss how the choice of the lattice geometry affects such methods. In particular, we show how to minimize systematic deviations from RMT by an optimal choice of the lattice geometry in the case of two light quark flavors. We illustrate our findings by determining ΣΣ and FF from lattice configurations with two dynamical overlap fermions generated by JLQCD, using two different lattice geometries.

preprint2011arXiv

Matching factors for Delta S=1 four-quark operators in RI/SMOM schemes

The non-perturbative renormalization of four-quark operators plays a significant role in lattice studies of flavor physics. For this purpose, we define regularization-independent symmetric momentum-subtraction (RI/SMOM) schemes for Delta S=1 flavor-changing four-quark operators and provide one-loop matching factors to the MS-bar scheme in naive dimensional regularization. The mixing of two-quark operators is discussed in terms of two different classes of schemes. We provide a compact expression for the finite one-loop amplitudes which allows for a straightforward definition of further RI/SMOM schemes.

preprint2010arXiv

The epsilon expansion at next-to-next-to-leading order with small imaginary chemical potential

We discuss chiral perturbation theory for two and three quark flavors in the epsilon expansion at next-to-next-to-leading order (NNLO) including a small imaginary chemical potential. We calculate finite-volume corrections to the low-energy constants ΣΣ and FF and determine the non-universal modifications of the theory, i.e., modifications that cannot be mapped to random matrix theory (RMT). In the special case of two quark flavors in an asymmetric box we discuss how to minimize the finite-volume corrections and non-universal modifications by an optimal choice of the lattice geometry. Furthermore we provide a detailed calculation of a special version of the massless sunset diagram at finite volume.

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.