The inclusion of heavy quarks, such as the charm quark, poses difficulties in Lattice QCD simulations. These need to be overcome in order to make predictions of hadronic observables including a charm quark such as the decay constants <i>f<sub>D</sub></i> and <i>f<sub>Ds</sub></i>. We first establish the region of validity of simulating heavy quarks with domain wall fermions (DWFs) in a feasibility study. We then carry out dynamical simulations with RBC/UKQCD’s <i>N<sub>f</sub></i><sub> </sub>= 2+1<i>f</i> ensembles at three lattice spacings (a <sup>−1</sup> = 1.73 − 2.77 GeV) including two physical pion mass ensembles. From this we make a prediction for the decay constants <i>f<sub>D</sub></i> and <i>f<sub>Ds</sub></i> and, using experimental input, the corresponding CKM matrix elements. In the first part of this work we investigate the suitability of domain wall fermions (DWF) as a lattice regularisation for heavy quarks. We generate four quenched QCD gauge ensembles with the tree-level improved Symanzik gauge action and inverse lattice spacings in the range 2.0 − 5.7 GeV. On these we carry out an exploratory study to identify a region in the DWF parameter space that displays minimal cut-off effects. We find this region for the domain wall height M<sub>5</sub> = 1.6 allowing for bare heavy quark masses satisfying <i>am<sub>h</sub></i> <<sub>∼ </sub>0.4, independent of the lattice spacing. Below this limit we maintain desirable features of DWF such as approximate chiral symmetry and O(<i>a</i>)-improvement. Based on this, we carry out a detailed scaling study of the decay constants of heavy-strange pseudoscalar mesons and the dispersion relation of heavy-heavy and heavy-strange pseudoscalar mesons. We find mild <i>a</i><sup>2</sup> discretisation effects for the heavy-strange decay constants and the heavy-strange dispersion relation. The cut-off effects for the heavy-heavy pseudoscalar dispersion relation are somewhat more pronounced. On our ensembles we find a 4 effects for heavy masses beyond the charm quark mass and large momenta (|<i>p</i>| ∼ 1.6 GeV). The findings of the pilot study establish the basis for simulations of charm quarks with domain wall fermions. In the second part of the presented work we simulate heavy-light and heavy-strange pseudoscalar mesons on RBC/UKQCD’s 2 + 1f ensembles with pion masses as low as 139 MeV and inverse lattice spacings ranging from 1.73 GeV to 2.77 GeV. From this data we extract the decay constants <i>f<sub>D</sub></i> and <i>f</i><sub>D</sub><i><sub>s</sub></i>. We devise two different fit ansatze to process this data and carry out a full systematic error analysis. We find <i>f<sub>D</sub></i> = 208.7(2.8)<sub>stat</sub> (<sup>+2.1</sup><sub>−1.8</sub>)<sub>sys</sub> MeV, <i>f<sub>Ds</sub></i> = 246.4(1.9)<sub>stat</sub> (<sup>+1.3</sup>−1.9)<sub>sys</sub> MeV and <i>f</i><sub style="font-style: italic;">Ds</sub>/<i>f</i><i><sub>D</sub></i> = 1.1667(77)<sub>stat</sub> (<sup>+60</sup>−46)<sub>sys</sub>. Finally, using experimental input we extract the corresponding CKM matrix elements. We find |V|<i><sub>cd</sub></i> = 0.2185(50)<sub>exp.</sub> (<sup>+35</sup>−37)<sub>lat</sub> and |V |<i><sub>cs</sub></i> = 1.011(16)<sub>exp.</sub> (<sup>+11</sup><sub>− 9</sub>)<sub>lat</sub> where the first error comes from the experimental input and the second from our determination of the decay constants. We compare our results with the existing literature and find good agreement between our central values and competitive errors. This work constitutes the basis of RBC/UKQCD’s wider charm physics program.
| Identifer | oai:union.ndltd.org:bl.uk/oai:ethos.bl.uk:749736 |
| Date | January 2017 |
| Creators | Tsang, Justus Tobias |
| Contributors | Juttner, Andreas |
| Publisher | University of Southampton |
| Source Sets | Ethos UK |
| Detected Language | English |
| Type | Electronic Thesis or Dissertation |
| Source | https://eprints.soton.ac.uk/405419/ |
Page generated in 0.362 seconds