Periodic Reporting for period 1 - HEAVYDWF (Charm and Bottom flavour physics with chiral lattice fermions)
Berichtszeitraum: 2015-10-01 bis 2017-09-30
Out of the six known quarks, the b-quarks play a special role: due to their large mass of 4.18 GeV [PDG] they allow to explore processes at otherwise not accessible energies and couple more strongly to the Higgs boson than the lighter quarks. Of particular interest are weak decays of b-quarks to another lighter quark where in addition a pair of leptons (electron, muon, tau or neutrino) is emitted. If the b-quark is bound as a quark-antiquark pair, the particle is named B-meson and such decays are classified as semileptonic. Lately, such semileptonic B decays have received much attention because several experiments (BaBar, Belle, and LHCb) see indications for violations of lepton flavor universality i.e. results exhibit an unexpected dependence whether the lepton is e.g. an electron or a tau. When combined these results amount to the strongest indications of new physics seen so far. However, improvements of the SM predictions are highly desirable which include predictions for so called form factors based on lattice field theory calculations. These form factors parametrize contributions from the strong force and the calculation faces the challenge to simulate b-quarks carrying a mass larger than the lattice cut-off currently possible to simulate.
The emphasis of this project is to improve simulations involving such heavy quarks and calculate quantities relevant to the experimental research program. First we explored using so called domain-wall fermions for simulations of heavy quarks and found them suitable for quark masses in the range of c-quarks (1.28 GeV [PDG]). With a modification (link smearing), up to half the mass of b-quarks can be simulated and the results extrapolated to the physical b-quark mass. In addition concepts to nonperturbatively renormalize heavy-light operators have been investigated. While such extrapolations are promising to determine simple quantities (decay constants, neutral meson mixing parameters), the application to a form factor calculations is more difficult. To explore improvements for the most timely form factor calculations involving the decay of a b-quark to a c-quark, we therefore chose to study B_(s) -> D_(s)^(*) l nu semileptonic decays using a different, effective action to simulate b-quarks at their physical mass but use the newly established concept of heavy domain-wall fermions to simulate c-quarks. We have completed the massive scale numerical simulations and implemented all parts of the analysis required for the different projects. Preliminary results have been presented at workshops and conferences and work is in progress to prepare a manuscript for submission to a peer-reviewed journal.
To test domain-wall fermions for simulating heavy quarks, several nonperturbative quantities have been calculated. Each calculation requires to implement required matrix elements describing the processes of interest. Next nonperturbative simulations using super-computers are carried out. These simulations are based on using Monte Carlo methods to statistically integrate the phase space using Feynman's path integral. Subsequently a statistical analysis is performed leading to quantities at fixed quark masses and values of the lattice spacing. After repeating the simulations for different quark masses and values of the lattice spacing, the results have to be renormalized before all data are analyzed together by carrying out a fit to extrapolate the quantities to the physical values of the quarks masses and the continuum limit. At this point systematic effects originating from the simulation or the analysis need to be estimated in order to obtain results with full error budget and suitable for comparison or combination with experimental results.
We have completed such simulations for quantities including decay constants for D, D_s, B, and B_s mesons, parameters for neutral meson mixing for D, B, and Bs_mesons as well as their ratio xi. Furthermore, form factors for semileptonic decays of B and Bs-mesons have been calculated. We have established the different parts of the analysis and are currently finalizing nonperturbative determinations of the renormalization factors as well as estimates for the systematic uncertainties. The project will be completed by preparing manuscripts for submission to peer-reviewed journals. In case of the semileptonic form factor calculation we are focusing at the determination of Bs -> Ds l nu and Bs -> K l nu decays because these will enable LHCb to obtain an independent determination of the ratio of CKM matrix elements |V_cb|/|V_ub|. Afterwards we will analyze B->D^(*) l nu decays and obtain theoretical predictions for the ratios R_D and R_D^* to address LFUV.
Preliminary results have been presented at the annual Lattice conferences in Southampton (2016) and Granada (2017), as well as at the major phenomenological conferences ICHEP 2016 in Chicago, and the EPS-HEP 2017 in Venice with summaries written up as proceedings. In addition OW was invited to report on this work at the LHCb UK meeting in Oxford in January 2017 as well as the LHCb workshop "Implications of LHCb measurements and future prospects" at CERN, Geneva, Switzerland in November 2017. Moreover, he has been invited to present his work at several seminars and other workshops.
PB and OW organized very successfully 9th International workshop on QCD and numerical algorithms workshop (QCDNA) in Edinburgh in 2016. Given the success, OW was invited to help organize the follow-up workshop QCDNA X in Coimbra (2017).