
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.

Deepinelastic scattering, in the laboratory and on the lattice, is most
instructive for understanding how the nucleon is built from quarks and gluons.
The longterm 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, allencompassing 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+1flavor 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 finiterange regularization. We properly account for
omitted disconnected quark contractions using a partiallyquenched 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 GellMannOkubo 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_nM_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 finiterange 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 lowenergy 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
latticeinformed 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 FeynmanHellmann (FH) theorem. At our SU(3)flavour
symmetric point, we find that the connected quark spin fractions are
universally in the range 5570\% 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 quarkline 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 quarkbilinear operators are computed
nonperturbatively 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 290465 MeV, allowing a safe and stable chiral
extrapolation. Emphasis is given in the subtraction of the wellknown 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 nonperturbative 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 nonperturbative 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$)
oneloop lattice artifacts.

The renormalized nexttoleadingorder (NLO) chiral lowenergy constant,
$L_{10}^r$, is determined in a complete nexttonexttoleadingorder (NNLO)
analysis, using a combination of lattice and continuum data for the flavor $ud$
$VA$ correlator and results from a recent chiral sumrule analysis of the
flavorbreaking combination of $ud$ and $us$ $VA$ correlator differences. The
analysis also fixes two combinations of NNLO lowenergy constants, the
determination of which is crucial to the precision achieved for $L_{10}^r$.
Using the results of the flavorbreaking chiral $VA$ 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 twoflavor 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 byproduct we
obtain the lowenergy constant $\bar{l}_4=4.2(1)$.

Recent analyses of flavorbreaking hadronic$\tau$decaybased sum rules
produce values of $\vert V_{us}\vert$ $\sim 3\sigma$ low compared to 3family
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 flavorbreaking
correlators and show the fixedscale prescription is clearly preferred.
Preliminary updates of the conventional $\tau$based, and related mixed
$\tau$electroproductiondatabased, sum rule analyses incorporating Bfactory
results for lowmultiplicity 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 3family 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 grouptheoretical 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 lightquark VA
correlator, supplemented by results from a chiral sumrule analysis of the
flavorbreaking flavor $ud$$us$ VA correlator difference, is shown to make
possible a highprecision NNLO determination of the renormalized NLO chiral
lowenergy constant $L_{10}^r$. Key to this determination is the ability to
simultaneously fix the two combinations of NNLO lowenergy constants also
entering the analysis. With current versions of the strange hadronic $\tau$
branching fractions required as input to the flavorbreaking VA 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 flavorbreaking
combinations of hadronic tau decay data, are significantly lower than either
expectations based on 3family 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 flavorbreaking
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 crosssection 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 VA correlator. A comparison of the results for the correlator to
the corresponding mildlymodeldependent 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 flavourblindness 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
flavourconserving and flavourchanging 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 Kellyinspired parametrizations for the
Q^2dependence 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 nonanalytic 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 nonperturbative
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.