• ### Charmed states and flavour symmetry breaking(1711.02485)

Nov. 7, 2017 hep-lat
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.
• ### Nucleon structure functions from lattice operator product expansion(1703.01153)

May 18, 2017 hep-ph, hep-ex, nucl-th, hep-lat
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.
• ### Electromagnetic form factors at large momenta from lattice QCD(1702.01513)

Feb. 6, 2017 hep-lat
Accessing hadronic form factors at large momentum transfers has traditionally presented a challenge for lattice QCD simulations. Here we demonstrate how a novel implementation of the Feynman-Hellmann method can be employed to calculate hadronic form factors in lattice QCD at momenta much higher than previously accessible. Our simulations are performed on a single set of gauge configurations with three flavours of degenerate mass quarks corresponding to $m_\pi \approx 470 \text{ MeV}$. We are able to determine the electromagnetic form factors of the pion and nucleon up to approximately $6 \text{ GeV}^2$, with results for $G_E/G_M$ in the proton agreeing well with experimental results.
• ### Electric form factors of the octet baryons from lattice QCD and chiral extrapolation(1403.1965)

Jan. 7, 2017 hep-ph, nucl-th, hep-lat
We apply a formalism inspired by heavy baryon chiral perturbation theory with finite-range regularization to dynamical $2+1-$flavor CSSM/QCDSF/UKQCD Collaboration lattice QCD simulation results for the electric form factors of the octet baryons. The electric form factor of each octet baryon is extrapolated to the physical pseudoscalar masses, after finite-volume corrections have been applied, at six fixed values of $Q^2$ in the range 0.2-1.3 GeV$^2$. The extrapolated lattice results accurately reproduce the experimental form factors of the nucleon at the physical point, indicating that omitted disconnected quark loop contributions are small. Furthermore, using the results of a recent lattice study of the magnetic form factors, we determine the ratio $\mu_p G^p_E/G^p_M$. This quantity decreases with $Q^2$ in a way qualitatively consistent with recent experimental results.
• ### Magnetic form factors of the octet baryons from lattice QCD and chiral extrapolation(1401.5862)

Jan. 7, 2017 hep-ph, nucl-th, hep-lat
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.
• ### Partially conserved axial vector current and applications(1612.04992)

Dec. 15, 2016 hep-lat
We investigate implications of the use of the point-split axial vector current derived from a Wilson like fermionic action. We compute the corresponding renormalization factor nonperturbatively for one beta value. The axial charge gA calculated from this nonlocal current is found to be nearer to the physical value than computed with the local axial vector current -- computed both on the same lattice with the same action.
• ### Flavour breaking effects in the pseudoscalar meson decay constants(1612.04798)

Dec. 14, 2016 hep-ph, hep-lat
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.
• ### Isospin splittings of meson and baryon masses from three-flavor lattice QCD + QED(1508.06401)

July 15, 2016 hep-ph, nucl-th, hep-lat
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.
• ### Review of lattice results concerning low-energy particle physics(1607.00299)

July 1, 2016 hep-ph, hep-lat
We review lattice results related to pion, kaon, D- and B-meson physics with the aim of making them easily accessible to the particle physics community. More specifically, we report on the determination of the light-quark masses, the form factor f+(0), arising in the semileptonic K -> pi transition at zero momentum transfer, as well as the decay constant ratio fK/fpi and its consequences for the CKM matrix elements Vus and Vud. Furthermore, we describe the results obtained on the lattice for some of the low-energy constants of SU(2)LxSU(2)R and SU(3)LxSU(3)R Chiral Perturbation Theory. We review the determination of the BK parameter of neutral kaon mixing as well as the additional four B parameters that arise in theories of physics beyond the Standard Model. The latter quantities are an addition compared to the previous review. For the heavy-quark sector, we provide results for mc and mb (also new compared to the previous review), as well as those for D- and B-meson decay constants, form factors, and mixing parameters. These are the heavy-quark quantities most relevant for the determination of CKM matrix elements and the global CKM unitarity-triangle fit. Finally, we review the status of lattice determinations of the strong coupling constant alpha_s.
• ### QED effects in the pseudoscalar meson sector(1509.00799)

Feb. 25, 2016 hep-ph, nucl-th, hep-lat
We present results on the pseudoscalar meson masses from a fully dynamical simulation of QCD+QED. We concentrate particularly on violations of isospin symmetry. We calculate the $\pi^+$-$\pi^0$ splitting and also look at other isospin violating mass differences. We have presented results for these isospin splittings in arXiv:1508.06401 [hep-lat]. In this paper we give more details of the techniques employed, discussing in particular the question of how much of the symmetry violation is due to QCD, arising from the different masses of the $u$ and $d$ quarks, and how much is due to QED, arising from the different charges of the quarks. This decomposition is not unique, it depends on the renormalisation scheme and scale. We suggest a renormalisation scheme in which Dashen's theorem for neutral mesons holds, so that the electromagnetic self-energies of the neutral mesons are zero, and discuss how the self-energies change when we transform to a scheme such as $\bar{MS}$, in which Dashen's theorem for neutral mesons is violated.
• ### Determining the scale in Lattice QCD(1512.05745)

Dec. 17, 2015 hep-lat
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.
• ### Applications of the Feynman-Hellmann theorem in hadron structure(1511.07090)

Nov. 23, 2015 hep-lat
The Feynman-Hellmann (FH) relation offers an alternative way of accessing hadronic matrix elements through artificial modifications to the QCD Lagrangian. In particular, a FH-motivated method provides a new approach to calculations of disconnected contributions to matrix elements and high-momentum nucleon and pion form factors. Here we present results for the total nucleon axial charge, including a statistically significant non-negative total disconnected quark contribution of around $-5\%$ at an unphysically heavy pion mass. Extending the FH relation to finite-momentum transfers, we also present calculations of the pion and nucleon electromagnetic form factors up to momentum transfers of around 7-8 GeV$^2$. Results for the nucleon are not able to confirm the existence of a sign change for the ratio $\frac{G_E}{G_M}$, but suggest that future calculations at lighter pion masses will provide fascinating insight into this behaviour at large momentum transfers.
• ### Improving the lattice axial vector current(1511.05304)

Nov. 17, 2015 hep-lat
For Wilson and clover fermions traditional formulations of the axial vector current do not respect the continuum Ward identity which relates the divergence of that current to the pseudoscalar density. Here we propose to use a point-split or one-link axial vector current whose divergence exactly satisfies a lattice Ward identity, involving the pseudoscalar density and a number of irrelevant operators. We check in one-loop lattice perturbation theory with SLiNC fermion and gauge plaquette action that this is indeed the case including order $O(a)$ effects. Including these operators the axial Ward identity remains renormalisation invariant. First preliminary results of a nonperturbative check of the Ward identity are also presented.
• ### Disconnected contributions to the spin of the nucleon(1508.06856)

Aug. 27, 2015 hep-lat
The spin decomposition of the proton is a long-standing topic of much interest in hadronic physics. Lattice QCD has had much success in calculating the connected contributions to the quark spin. However, complete calculations, which necessarily involve gluonic and strange-quark contributions, still present some challenges. These "disconnected" contributions typically involve small signals hidden against large statistical backgrounds and rely on computationally intensive stochastic techniques. In this work we demonstrate how a Feynman-Hellmann approach may be used to calculate such quantities, by measuring shifts in the proton energy arising from artificial modifications to the QCD action. We find a statistically significant non-zero result for the disconnected quark spin contribution to the proton of about -5% at a pion mass of 470 MeV.
• ### SU(3) breaking in hyperon transition vector form factors(1508.06923)

Aug. 27, 2015 hep-th, nucl-th, hep-lat
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.
• ### Wilson flow and scale setting from lattice QCD(1508.05916)

Aug. 24, 2015 hep-lat
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.
• ### Reply to "Comment on `Lattice determination of Sigma - Lambda mixing' "(1507.07825)

July 28, 2015 hep-ph, nucl-ex, nucl-th, hep-lat
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.
• ### The electric dipole moment of the neutron from 2+1 flavor lattice QCD(1502.02295)

July 7, 2015 hep-ph, nucl-th, hep-lat
We compute the electric dipole moment d_n of the neutron from a fully dynamical simulation of lattice QCD with 2+1 flavors of clover fermions and nonvanishing theta term. The latter is rotated into the pseudoscalar density in the fermionic action using the axial anomaly. To make the action real, the vacuum angle theta is taken to be purely imaginary. The physical value of d_n is obtained by analytic continuation. We find d_n = -3.8(2)(9) x 10^{-16} [theta e cm], which, when combined with the experimental limit on d_n, leads to the upper bound theta < 7.6 x 10^{-11}.
• ### A lattice determination of Sigma - Lambda mixing(1411.7665)

April 16, 2015 hep-ph, hep-lat
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.
• ### Charge Symmetry Violation in the Electromagnetic Form Factors of the Proton(1503.01142)

March 3, 2015 hep-ph, nucl-th, hep-lat
Experimental tests of QCD through its predictions for the strange-quark content of the proton have been drastically restricted by our lack of knowledge of the violation of charge symmetry (CSV). We find unexpectedly tiny CSV in the proton's electromagnetic form factors by performing the first extraction of these quantities based on an analysis of lattice QCD data. The resulting values are an order of magnitude smaller than current bounds on proton strangeness from parity violating electron-proton scattering experiments. This result paves the way for a new generation of experimental measurements of the proton's strange form factors to challenge the predictions of QCD.
• ### Determination of the strange nucleon form factors(1403.6537)

Feb. 16, 2015 hep-ph, nucl-th, hep-lat
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.
• ### Connected and disconnected quark contributions to hadron spin(1412.6569)

Dec. 20, 2014 hep-lat
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.
• ### Determining Sigma - Lambda mixing(1412.0970)

Dec. 2, 2014 hep-lat
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.
• ### A novel approach to nonperturbative renormalization of singlet and nonsinglet lattice operators(1410.3078)

Oct. 12, 2014 hep-lat
A novel method for nonperturbative renormalization of lattice operators is introduced, which lends itself to the calculation of renormalization factors for nonsinglet as well as singlet operators. The method is based on the Feynman-Hellmann relation, and involves computing two-point correlators in the presence of generalized background fields arising from introducing additional operators into the action. As a first application, and test of the method, we compute the renormalization factors of the axial vector current $A_\mu$ and the scalar density $S$ for both nonsinglet and singlet operators for $N_f=3$ flavors of SLiNC fermions. For nonsinglet operators, where a meaningful comparison is possible, perfect agreement with recent calculations using standard three-point function techniques is found.
• ### Renormalization of local quark-bilinear operators for Nf=3 flavors of SLiNC fermions(1408.6047)

Aug. 26, 2014 hep-lat
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.