• Extending the SU(3) flavour symmetry breaking expansion from up, down and strange sea quark masses to partially quenched valence quark masses allows an extrapolation to the charm quark mass. This approach leads to a determination of charmed quark hadron masses and decay constants. We describe our recent progress and give preliminary results in particular with regard to the recently discovered doubly charmed baryon by the LHCb Collaboration.
  • Deep-inelastic scattering, in the laboratory and on the lattice, is most instructive for understanding how the nucleon is built from quarks and gluons. The long-term goal is to compute the associated structure functions from first principles. So far this has been limited to model calculations. In this Letter we propose a new method to compute the structure functions directly from the virtual, all-encompassing Compton amplitude, utilizing the operator product expansion. This overcomes issues of renormalization and operator mixing, which so far have hindered lattice calculations of power corrections and higher moments.
  • We present a 2+1--flavor lattice QCD calculation of the electromagnetic Dirac and Pauli form factors of the octet baryons. The magnetic Sachs form factor is extrapolated at six fixed values of $Q^2$ to the physical pseudoscalar masses and infinite volume using a formulation based on heavy baryon chiral perturbation theory with finite-range regularization. We properly account for omitted disconnected quark contractions using a partially-quenched effective field theory formalism. The results compare well with the experimental form factors of the nucleon and the magnetic moments of the octet baryons.
  • The SU(3) flavour symmetry breaking expansion in up, down and strange quark masses is extended from hadron masses to meson decay constants. This allows a determination of the ratio of kaon to pion decay constants in QCD. Furthermore when using partially quenched valence quarks the expansion is such that SU(2) isospin breaking effects can also be determined. It is found that the lowest order SU(3) flavour symmetry breaking expansion (or Gell-Mann-Okubo expansion) works very well. Simulations are performed for 2+1 flavours of clover fermions at four lattice spacings.
  • Lattice QCD simulations are now reaching a precision where isospin breaking effects become important. Previously, we have developed a program to systematically investigate the pattern of flavor symmetry beaking within QCD and successfully applied it to meson and baryon masses involving up, down and strange quarks. In this Letter we extend the calculations to QCD + QED and present our first results on isospin splittings in the pseudoscalar meson and baryon octets. In particular, we obtain the nucleon mass difference of $M_n-M_p=1.35(18)(8)\,\mbox{MeV}$ and the electromagnetic contribution to the pion splitting $M_{\pi^+}-M_{\pi^0}=4.60(20)\,\mbox{MeV}$. Further we report first determination of the separation between strong and electromagnetic contributions in the $\bar{MS}$ scheme.
  • We discuss scale setting in the context of 2+1 dynamical fermion simulations where we approach the physical point in the quark mass plane keeping the average quark mass constant. We have simulations at four beta values, and after determining the paths and lattice spacings, we give an estimation of the phenomenological values of various Wilson flow scales.
  • We present a calculation of the SU(3)-breaking corrections to the hyperon transition vector form factors to $\mathcal{O}(p^4)$ in heavy baryon chiral perturbation theory with finite-range regularisation. Both octet and decuplet degrees of freedom are included. We formulate a chiral expansion at the kinematic point $Q^2=-(M_{B_1}-M_{B_2})^2$, which can be conveniently accessed in lattice QCD. The two unknown low-energy constants at this point are constrained by lattice QCD simulation results for the $\Sigma^-\rightarrow n$ and $\Xi^0\rightarrow \Sigma^+$ transition form factors. Hence we determine lattice-informed values of $f_1$ at the physical point. This work constitutes progress towards the precise determination of $|V_{us}|$ from hyperon semileptonic decays.
  • We give a determination of the phenomenological value of the Wilson (or gradient) flow scales t0 and w0 for 2+1 flavours of dynamical quarks. The simulations are performed keeping the average quark mass constant, which allows the approach to the physical point to be made in a controlled manner. O(a) improved clover fermions are used and together with four lattice spacings this allows the continuum extrapolation to be taken.
  • In this Reply, we respond to the above Comment. Our computation [Phys. Rev. D 91 (2015) 074512] only took into account pure QCD effects, arising from quark mass differences, so it is not surprising that there are discrepancies in isospin splittings and in the Sigma - Lambda mixing angle. We expect that these discrepancies will be smaller in a full calculation incorporating QED effects.
  • Isospin breaking effects in baryon octet (and decuplet) masses are due to a combination of up and down quark mass differences and electromagnetic effects and lead to small mass splittings. Between the Sigma and Lambda this mass splitting is much larger, this being mostly due to their different wavefunctions. However when isospin is broken, there is a mixing between between these states. We describe the formalism necessary to determine the QCD mixing matrix and hence find the mixing angle and mass splitting between the Sigma and Lambda particles due to QCD effects.
  • The strange contribution to the electric and magnetic form factors of the nucleon is determined at a range of discrete values of $Q^2$ up to $1.4$ GeV$^2$. This is done by combining recent lattice QCD results for the electromagnetic form factors of the octet baryons with experimental determinations of those quantities. The most precise result is a small negative value for the strange magnetic moment: $G_M^s(Q^2=0) = -0.07\pm0.03\,\mu_N$. At larger values of $Q^2$ both the electric and magnetic form factors are consistent with zero to within $2$-sigma.
  • By introducing an external spin operator to the fermion action, the quark spin fractions of hadrons are determined from the linear response of the hadron energies using the Feynman-Hellmann (FH) theorem. At our SU(3)-flavour symmetric point, we find that the connected quark spin fractions are universally in the range 55-70\% for vector mesons and octet and decuplet baryons. There is an indication that the amount of spin suppression is quite sensitive to the strength of SU(3) breaking. We also present first preliminary results applying the FH technique to calculations of quark-line disconnected contributions to hadronic matrix elements of axial and tensor operators. At the SU(3)-flavour symmetric point we find a small negative contribution to the nucleon spin from disconnected quark diagrams, while the corresponding tensor matrix elements are consistent with zero.
  • SU2 isospin breaking effects in baryon octet (and decuplet) masses are due to a combination of up and down quark mass differences and electromagnetic effects. These mass differences are small. Between the Sigma and Lambda the mass splitting is much larger, but this is mostly due to their different wavefunctions. However there is now also mixing between these states. We determine the QCD mixing matrix and hence find the mixing angle and mass splitting.
  • The renormalization factors of local quark-bilinear operators are computed non-perturbatively for $N_f=3$ flavors of SLiNC fermions, with emphasis on the various procedures for the chiral and continuum extrapolations. The simulations are performed at a lattice spacing $a=0.074$ fm, and for five values of the pion mass in the range of 290-465 MeV, allowing a safe and stable chiral extrapolation. Emphasis is given in the subtraction of the well-known pion pole which affects the renormalization factor of the pseudoscalar current. We also compute the inverse propagator and the Green's functions of the local bilinears to one loop in perturbation theory. We investigate lattice artifacts by computing them perturbatively to second order as well as to all orders in the lattice spacing. The renormalization conditions are defined in the RI$'$-MOM scheme, for both the perturbative and non-perturbative results. The renormalization factors, obtained at different values of the renormalization scale, are translated to the ${\bar{\rm MS}}$ scheme and are evolved perturbatively to 2 GeV. Any residual dependence on the initial renormalization scale is eliminated by an extrapolation to the continuum limit. We also study the various sources of systematic errors. Particular care is taken in correcting the non-perturbative estimates by subtracting lattice artifacts computed to one loop perturbation theory using the same action. We test two different methods, by subtracting either the ${\cal O}(g^2\,a^2)$ contributions, or the complete (all orders in $a$) one-loop lattice artifacts.
  • The renormalized next-to-leading-order (NLO) chiral low-energy constant, $L_{10}^r$, is determined in a complete next-to-next-to-leading-order (NNLO) analysis, using a combination of lattice and continuum data for the flavor $ud$ $V-A$ correlator and results from a recent chiral sum-rule analysis of the flavor-breaking combination of $ud$ and $us$ $V-A$ correlator differences. The analysis also fixes two combinations of NNLO low-energy constants, the determination of which is crucial to the precision achieved for $L_{10}^r$. Using the results of the flavor-breaking chiral $V-A$ sum rule obtained with current versions of the strange hadronic $\tau$ branching fractions as input, we find $L_{10}^r(m_\rho )\, =\, -0.00346(32)$. This result represents the first NNLO determination of $L_{10}^r$ having all inputs under full theoretical and/or experimental control, and the best current precision for this quantity.
  • The axial charge of the nucleon $g_A$ and the pion decay constant $f_\pi$ are computed in two-flavor lattice QCD. The simulations are carried out on lattices of various volumes and lattice spacings. Results are reported for pion masses as low as $m_\pi=130\,\mbox{MeV}$. Both quantities, $g_A$ and $f_\pi$, suffer from large finite size effects, which to leading order ChEFT and ChPT turn out to be identical. By considering the naturally renormalized ratio $g_A/f_\pi$, we observe a universal behavior as a function of decreasing quark mass. From extrapolating the ratio to the physical point, we find $g_A^R=1.29(5)(3)$, using the physical value of $f_\pi$ as input and $r_0=0.50(1)$ to set the scale. In a subsequent calculation we attempt to extrapolate $g_A$ and $f_\pi$ separately to the infinite volume. Both volume and quark mass dependencies of $g_A$ and $f_\pi$ are found to be well decribed by ChEFT and ChPT. We find at the physical point $g_A^R=1.24(4)$ and $f_\pi^R=89.6(1.1)(1.8)\,\mbox{MeV}$. Both sets of results are in good agreement with experiment. As a by-product we obtain the low-energy constant $\bar{l}_4=4.2(1)$.
  • Recent analyses of flavor-breaking hadronic-$\tau$-decay-based sum rules produce values of $\vert V_{us}\vert$ $\sim 3\sigma$ low compared to 3-family unitarity expectations. An unresolved systematic issue is the significant variation in $\vert V_{us}\vert$ produced by different prescriptions for treating the slowly converging $D=2$ OPE series. We investigate the reliability of these prescriptions using lattice data for various flavor-breaking correlators and show the fixed-scale prescription is clearly preferred. Preliminary updates of the conventional $\tau$-based, and related mixed $\tau$-electroproduction-data-based, sum rule analyses incorporating B-factory results for low-multiplicity strange $\tau$ decay mode distributions are then performed. Use of the preferred FOPT $D=2$ OPE prescription is shown to significantly reduce the discrepancy between 3-family unitarity expectations and the sum rule results.
  • By extending the SU(3) flavour symmetry breaking expansion from up, down and strange sea quark masses to partially quenched valence quark masses we propose a method to determine charmed quark hadron masses including possible QCD isospin breaking effects. Initial results for some open charmed pseudoscalar meson states and singly and doubly charmed baryon states are encouraging and demonstrate the potential of the procedure. Essential for the method is the determination of the scale using singlet quantities, and to this end we also give here a preliminary estimation of the recently introduced Wilson flow scales.
  • Lattice QCD simulations are now reaching a precision where electromagnetic corrections from QED become important. In investigating the effects of SU(3) breaking due to quark mass differences within QCD, a group-theoretical analysis of the mass dependence greatly helped us organize our results. We now do the same with electromagnetic charge effects by extending the calculations to dynamical 1+1+1 flavor QCD + QED.
  • A combination of lattice and continuum data for the light-quark V-A correlator, supplemented by results from a chiral sum-rule analysis of the flavor-breaking flavor $ud$-$us$ V-A correlator difference, is shown to make possible a high-precision NNLO determination of the renormalized NLO chiral low-energy constant $L_{10}^r$. Key to this determination is the ability to simultaneously fix the two combinations of NNLO low-energy constants also entering the analysis. With current versions of the strange hadronic $\tau$ branching fractions required as input to the flavor-breaking V-A sum rule, we find $L_{10}^r(m_\rho ) = -0.00346(29)$. This represents both the best current precision for $L_{10}^r$, and the first NNLO determination having all errors under full control.
  • Recent sum rule determinations of |V_us|, employing flavor-breaking combinations of hadronic tau decay data, are significantly lower than either expectations based on 3-family unitarity or determinations from K_ell3 and Gamma[K_mu2]/Gamma[pi_mu2]. We use lattice data to investigate the accuracy/reliability of the OPE representation of the flavor-breaking correlator combination entering the tau decay analyses. The behavior of an alternate correlator combination, constructed to reduce problems associated with the slow convergence of the D = 2 OPE series, and entering an alternate sum rule requiring both electroproduction cross-section and hadronic tau decay data, is also investigated. Preliminary updates of both analyses, with the lessons learned from the lattice data in mind, are also presented.
  • We present preliminary results on extractions of the chiral LECs L_10 and C_87 and constraints on the excited pseudoscalar state pi(1300) and pi(1800) decay constants obtained from an analysis of lattice data for the flavor ud light quark V-A correlator. A comparison of the results for the correlator to the corresponding mildly-model-dependent continuum results (based primarily on experimental hadronic tau decay data) is also given
  • By considering a flavour expansion about the SU(3)-flavour symmetric point, we investigate how flavour-blindness constrains octet baryon matrix elements after SU(3) is broken by the mass difference between the strange and light quarks. We find the expansions to be highly constrained along a mass trajectory where the singlet quark mass is held constant, which proves beneficial for extrapolations of 2+1 flavour lattice data to the physical point. We investigate these effects numerically via a lattice calculation of the flavour-conserving and flavour-changing matrix elements of the vector and axial operators between octet baryon states.
  • We present a comprehensive analysis of the electromagnetic form factors of the nucleon from a lattice simulation with two flavors of dynamical O(a)-improved Wilson fermions. A key feature of our calculation is that we make use of an extensive ensemble of lattice gauge field configurations with four different lattice spacings, multiple volumes, and pion masses down to m_\pi ~ 180 MeV. We find that by employing Kelly-inspired parametrizations for the Q^2-dependence of the form factors, we are able to obtain stable fits over our complete ensemble. Dirac and Pauli radii and the anomalous magnetic moments of the nucleon are extracted and results at light quark masses provide evidence for chiral non-analytic behavior in these fundamental observables.
  • By determining the quark momentum fractions of the octet baryons from N_f=2+1 lattice simulations, we are able to predict the degree of charge symmetry violation in the parton distribution functions of the nucleon. This is of importance, not only as a probe of our understanding of the non-perturbative structure of the proton but also because such a violation constrains the accuracy of global fits to parton distribution functions and hence the accuracy with which, for example, cross sections at the LHC can be predicted. A violation of charge symmetry may also be critical in cases where symmetries are used to guide the search for physics beyond the Standard Model.