
The axion is a hypothetical elementary particle postulated by the
PecceiQuinn theory to resolve the strong CP problem in QCD. If axions exist
and have low mass, they are a candidate for dark matter as well. So far our
knowledge of the properties of the QCD axion rests on semiclassical arguments
and effective theory. In this Letter we perform a fully dynamical investigation
of the PecceiQuinn theory, focussing on the impact of QCD on key axion
parameters, by simulating the PecceiQuinnWeinbergWilczek action on the
lattice. The results of the simulation, including the axion mass and effective
potential, are found to be in contradiction with current axion phenomenology
and question the validity and use of the PecceiQuinn theory.

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.

The generalization of the FeynmanHellmann theorem for resonance states in
quantum field theory is derived. On the basis of this theorem, a criterion is
proposed to study the possible exotic nature of certain hadronic states
emerging in QCD. It is shown that this proposal is supported by explicit
calculations in Chiral Perturbation Theory and by large$N_c$ arguments.
Analyzing recent lattice data on the quark mass dependence in the pseudoscalar,
vector meson, baryon octet and baryon decuplet sectors, we conclude that, as
expected, these are predominately quarkmodel states, albeit the corrections
are nonnegligible.

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.

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 FeynmanHellmann 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.

We apply a formalism inspired by heavy baryon chiral perturbation theory with
finiterange 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 finitevolume
corrections have been applied, at six fixed values of $Q^2$ in the range
0.21.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.

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.

We investigate implications of the use of the pointsplit 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.

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 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 [heplat]. 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
selfenergies of the neutral mesons are zero, and discuss how the selfenergies
change when we transform to a scheme such as $\bar{MS}$, in which Dashen's
theorem for neutral mesons is violated.

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.

The FeynmanHellmann (FH) relation offers an alternative way of accessing
hadronic matrix elements through artificial modifications to the QCD
Lagrangian. In particular, a FHmotivated method provides a new approach to
calculations of disconnected contributions to matrix elements and highmomentum
nucleon and pion form factors. Here we present results for the total nucleon
axial charge, including a statistically significant nonnegative total
disconnected quark contribution of around $5\%$ at an unphysically heavy pion
mass. Extending the FH relation to finitemomentum transfers, we also present
calculations of the pion and nucleon electromagnetic form factors up to
momentum transfers of around 78 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.

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
pointsplit or onelink axial vector current whose divergence exactly satisfies
a lattice Ward identity, involving the pseudoscalar density and a number of
irrelevant operators. We check in oneloop 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.

The spin decomposition of the proton is a longstanding 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 strangequark 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 FeynmanHellmann 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 nonzero result for the
disconnected quark spin contribution to the proton of about 5% at a pion mass
of 470 MeV.

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.

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}.

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.

Experimental tests of QCD through its predictions for the strangequark
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 electronproton 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.

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.

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
FeynmanHellmann relation, and involves computing twopoint 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 threepoint function techniques is found.