• We review the domain wall charm physics program of the RBC and UKQCD collaborations based on simulations including ensembles with physical pion mass. We summarise our current set-up and present a status update on the decay constants $f_D$, $f_{D_s}$, the charm quark mass, heavy-light and heavy-strange bag parameters and the ratio $\xi$.
  • We report on a calculation of the B*Bpi coupling in lattice QCD. The strong matrix element for a B* to Bpi transition is directly related to the leading order low-energy constant in heavy meson chiral perturbation theory (HMChPT) for B mesons. We carry out our calculation directly at the b-quark mass using a non-perturbatively tuned clover action that controls discretization effects of order pa and (ma)^n for all n. Our analysis is performed on RBC/UKQCD gauge configurations using domain-wall fermions and the Iwasaki gauge action at two lattice spacings of ainverse = 1.729(25) GeV, ainverse = 2.281(28) GeV, and unitary pion masses down to 290 MeV. We achieve good statistical precision and control all systematic uncertainties, giving a final result for the HMChPT coupling g_b = 0.56(3)stat(7)sys in the continuum and at the physical light-quark masses. This is the first calculation performed directly at the physical b-quark mass and lies in the region one would expect from carrying out an interpolation between previous results at the charm mass and at the static point.
  • We propose to construct a chirally broken model based on the infrared fixed point of a conformal system by raising the mass of some flavors while keeping the others massless. In the infrared limit the massive fermions decouple and the massless fermions break chiral symmetry. The running coupling of this system "walks" and the energy range of walking can be tuned by the mass of the heavy flavors. Renormalization group considerations predict that the spectrum of such a system shows hyperscaling. We have studied a model with four light and eight heavy flavors coupled to SU(3) gauge fields and verified the above expectations. We determined the mass of several hadronic states and found that some of them are in the 2-3 TeV range if the scale is set by the pseudoscalar decay constant $F_\pi \approx 250$ GeV. The $0^{++}$ scalar state behaves very differently from the other hadronic states. In most of our simulations it is nearly degenerate with the pion and we estimate its mass to be less than half of the vector resonance mass.
  • The dynamical origin of electroweak symmetry breaking is an open question with many possible theoretical explanations. Strongly coupled systems predicting the Higgs boson as a bound state of a new gauge-fermion interaction form one class of candidate models. Due to increased statistics, LHC run II will further constrain the phenomenologically viable models in the near future. In the meanwhile it is important to understand the general properties and specific features of the different competing models. In this work we discuss many-flavor gauge-fermion systems that contain both massless (light) and massive fermions. The former provide Goldstone bosons and trigger electroweak symmetry breaking, while the latter indirectly influence the infrared dynamics. Numerical results reveal that such systems can exhibit a light $0^{++}$ isosinglet scalar, well separated from the rest of the spectrum. Further, when we set the scale via the $vev$ of electroweak symmetry breaking, we predict a 2 TeV vector resonance which could be a generic feature of SU(3) gauge theories.
  • We calculate the $B \to\pi\ell\nu$ and $B_s \to K \ell\nu$ form factors in dynamical lattice QCD. We use the (2+1)-flavor RBC-UKQCD gauge-field ensembles generated with the domain-wall fermion and Iwasaki gauge actions. For the $b$ quarks we use the anisotropic clover action with a relativistic heavy-quark interpretation. We analyze two lattice spacings $a \approx 0.11, 0.086$ fm and unitary pion masses as light as $M_\pi \approx 290$ MeV. We simultaneously extrapolate our numerical results to the physical light-quark masses and to the continuum and interpolate in the pion/kaon energy using SU(2) "hard-pion" chiral perturbation theory. We provide complete error budgets for the form factors $f_+(q^2)$ and $f_0(q^2)$ at three momenta that span the $q^2$ range accessible in our numerical simulations. We extrapolate these results to $q^2 = 0$ using a model-independent $z$-parametrization and present our final form factors as the $z$-coefficients and the matrix of correlations between them. Our results agree with other lattice determinations using staggered light quarks and provide important independent cross-checks. Both $B \to\pi\ell\nu$ and $B_s \to K \ell\nu$ decays enable a determination of the CKM matrix element $|V_{ub}|$. To illustrate this, we perform a combined $z$-fit of our numerical $B\to\pi\ell\nu$ form-factor data with the experimental branching-fraction measurements leaving the relative normalization as a free parameter; we obtain $|V_{ub}| = 3.61(32) \times 10^{-3}$, where the error includes statistical and systematic uncertainties. This approach can be applied to $B_s\to K \ell\nu$ decay to determine $|V_{ub}|$ once the process has been measured experimentally. Finally, in anticipation of future measurements, we make predictions for $B \to \pi\ell\nu$ and $B_s\to K \ell\nu$ Standard-Model differential branching fractions and forward-backward asymmetries.
  • We present a novel technique for the determination of the topological susceptibility (related to the variance of the distribution of global topological charge) from lattice gauge theory simulations, based on maximum-likelihood analysis of the Markov-chain Monte Carlo time series. This technique is expected to be particularly useful in situations where relatively few tunneling events are observed. Restriction to a lattice subvolume on which topological charge is not quantized is explored, and may lead to further improvement when the global topology is poorly sampled. We test our proposed method on a set of lattice data, and compare it to traditional methods.
  • We report on a calculation of the $B^*B\pi$ coupling in lattice QCD. The strong matrix element $\langle B \pi | B^*\rangle$ is directly related to the leading order low-energy constant in heavy meson chiral perturbation theory (HM$\chi$PT) for $B$-mesons. We carry out our calculation directly at the $b$-quark mass using a non-perturbatively tuned clover action that controls discretisation effects of order $|\vec{p}a|$ and $(ma)^n$ for all $n$. Our analysis is performed on RBC/UKQCD gauge configurations using domain wall fermions and the Iwasaki gauge action at two lattice spacings of $a^{-1}=1.73(3)$ GeV, $a^{-1}=2.28(3)$ GeV, and unitary pion masses down to 290 MeV. We achieve good statistical precision and control all systematic uncertainties, giving a final result for the HM$\chi$PT coupling $g_b = 0.569(48)_{stat}(59)_{sys}$ in the continuum and at the physical light-quark masses. This is the first calculation performed directly at the physical $b$-quark mass and lies in the region one would expect from carrying out an interpolation between previous results at the charm mass and at the static point.
  • We discuss some large effects of dynamical fermions. One is a cutoff effect, others concern the contribution of multi-pion states to correlation functions and are expected to survive the continuum limit. We then turn to the preparation for simulations at small lattice spacings which we are planning down to around a=0.04fm in order to understand the size of O(a^2)-effects of the standard O(a)-improved theory. The dependence of the lattice spacing on the bare coupling is determined through the Schr"odinger functional renormalized coupling.
  • Under the constraints of HQET the equations of motion of heavy meson distribution amplitudes of definite geometric twist, using the knowledge of their off-cone structure, are reformulated as a set of algebraic equations. Together with equations due to various Dirac structures various relations between the (sets of) independent two- and three particle distribution amplitudes of definite geometric twist are derived and presented using both the notion of (double) Mellin moments and re-summed non-local distribution amplitudes. Resolving these relations for the independent two-particle moments in terms of three-particle double moments we confirmed the representation of $\Phi_{\pm|n}$ by Kawamura et al. (Phys. Lett. B \bf{523} (2001) 111).