-
The spin-independent and transversity generalised form factors (GFFs) of the
$\phi$ meson are studied using lattice QCD calculations with light quark masses
corresponding to a pion mass $m_\pi\sim450(5)$ MeV. One transversity and three
spin-independent GFFs related to the lowest moments of leading-twist
spin-independent and transversity gluon distributions are obtained at six
non-zero values of the momentum transfer up to 1.2 GeV$^2$. These quantities
are compared with the analogous spin-independent quark GFFs and the
electromagnetic form factors determined on the same lattice ensemble. The
results show quantitative distinction between the spatial distribution of
transversely polarised gluons, unpolarised gluons, and quarks, and point the
way towards further investigations of the gluon structure of nucleons and
nuclei.
-
We renew our intent to submit a proposal to perform a search for a non-zero
value of the unmeasured hadronic double helicity flip structure function
$\Delta(x,Q^2)$, predicted to be sensitive to gluons in the nucleus. This would
be performed with an unpolarized electron beam and transversely polarized,
spin-1, nuclear target. This structure function was first identified by Jaffe
and Manohar in 1989 as "a clear signature for exotic gluonic components in the
target," and a recent lattice QCD result by our collaborators has prompted
renewed interest in the topic. An inclusive search with deep inelastic
scattering, below $x$ of 0.3, via single spin tensor asymmetries may be
feasible using the CEBAF 12 GeV electron beam and JLab/UVa solid polarized
target, and would represent the first experimental exploration of this
quantity.
-
We present an exploratory study of the gluonic structure of the $\phi$ meson
using lattice QCD (LQCD). This includes the first investigation of gluonic
transversity via the leading moment of the twist-two double-helicity-flip
gluonic structure function $\Delta(x,Q^2)$. This structure function only exists
for targets of spin $J\ge1$ and does not mix with quark distributions at
leading twist, thereby providing a particularly clean probe of gluonic degrees
of freedom. We also explore the gluonic analogue of the Soffer bound which
relates the helicity flip and non-flip gluonic distributions, finding it to be
saturated at the level of 80%. This work sets the stage for more complex LQCD
studies of gluonic structure in the nucleon and in light nuclei where
$\Delta(x,Q^2)$ is an 'exotic glue' observable probing gluons in a nucleus not
associated with individual nucleons.
-
Quarkonium-nucleus systems are composed of two interacting hadronic states
without common valence quarks, which interact primarily through multi-gluon
exchanges, realizing a color van der Waals force. We present lattice QCD
calculations of the interactions of strange and charm quarkonia with light
nuclei. Both the strangeonium-nucleus and charmonium-nucleus systems are found
to be relatively deeply bound when the masses of the three light quarks are set
equal to that of the physical strange quark. Extrapolation of these results to
the physical light-quark masses suggests that the binding energy of charmonium
to nuclear matter is B < 40 MeV.
-
It has been argued that the leading scalar-isoscalar WIMP-nucleus
interactions receive parametrically enhanced contributions in the context of
nuclear effective field theories. These contributions arise from meson-exchange
currents (MECs) and potentially modify the impulse approximation estimates of
these interactions by 10--60%. We point out that these MECs also contribute to
the quark mass dependence of nuclear binding energies, that is, nuclear
\sigma-terms. In this work, we use recent lattice QCD calculations of the
binding energies of the deuteron, He-3 and He-4 at pion masses near 500 MeV and
800 MeV, combined with the experimentally determined binding energies at the
physical point, to provide approximate determinations of the \sigma-terms for
these light nuclei. For each nucleus, we find that the deviation of the
corresponding nuclear \sigma-term from the single-nucleon estimate is at the
few percent level, in conflict with the conjectured enhancement. As a
consequence, lattice QCD calculations currently indicate that the cross
sections for scalar-isoscalar WIMP-nucleus interactions arising from
fundamental WIMP interactions with quarks do not suffer from significant
uncertainties due to enhanced meson-exchange currents.
-
We present the results of lattice QCD calculations of the magnetic moments of
the lightest nuclei, the deuteron, the triton and ${}^3$He, along with those of
the neutron and proton. These calculations, performed at quark masses
corresponding to $m_\pi \sim 800$ MeV, reveal that the structure of these
nuclei at unphysically heavy quark masses closely resembles that at the
physical quark masses. In particular, we find that the magnetic moment of
${}^3$He differs only slightly from that of a free neutron, as is the case in
nature, indicating that the shell-model configuration of two spin-paired
protons and a valence neutron captures its dominant structure. Similarly a
shell-model-like moment is found for the triton, $\mu_{{}^3{\rm H}} \sim
\mu_p$. The deuteron magnetic moment is found to be equal to the nucleon
isoscalar moment within the uncertainties of the calculations.
-
The scattering lengths and effective ranges that describe low-energy
nucleon-nucleon scattering are calculated in the limit of SU(3)-flavor symmetry
at the physical strange-quark mass with Lattice Quantum Chromodynamics. The
calculations are performed with an isotropic clover discretization of the quark
action in three volumes with spatial extents of L \sim 3.4 fm, 4.5fm and 6.7
fm, and with a lattice spacing of b \sim 0.145 fm. With determinations of the
energies of the two-nucleon systems (both of which contain bound states at
these up and down quark masses) at rest and moving in the lattice volume,
Luscher's method is used to determine the low-energy phase shifts in each
channel, from which the scattering length and effective range are obtained. The
scattering parameters, in the 1S0 channel are found to be m_pi a^(1S0) =
9.50^{+0.78}_{-0.69}^{+1.10}_{-0.80} and m_pi r^(1S0) =
{4.61^{+0.29}_{-0.31}^{+0.24}_{-0.26}, and in the 3S1 channel are m_pi a^(3S1)
= 7.45^{+0.57}_{-0.53}^{+0.71}_{-0.49} and m_pi r^(3S1) =
3.71^{+0.28}_{-0.31}^{+0.28}_{-0.35}. These values are consistent with the
two-nucleon system exhibiting Wigner's supermultiplet symmetry, which becomes
exact in the limit of large-N_c. In both spin channels, the phase shifts change
sign at higher momentum, near the start of the t-channel cut, indicating that
the nuclear interactions have a repulsive core even at the SU(3)-symmetric
point.
-
The binding energies of a range of nuclei and hypernuclei with atomic number
A <= 4 and strangeness |s| <= 2, including the deuteron, di-neutron,
H-dibaryon, 3He, Lambda 3He, Lambda 4He, and Lambda Lambda 4He, are calculated
in the limit of flavor-SU(3) symmetry at the physical strange quark mass with
quantum chromodynamics (without electromagnetic interactions). The nuclear
states are extracted from Lattice QCD calculations performed with n_f=3
dynamical light quarks using an isotropic clover discretization of the
quark-action in three lattice volumes of spatial extent L ~ 3.4 fm, 4.5 fm and
6.7 fm, and with a single lattice spacing b ~ 0.145 fm.
-
An analysis of the pion mass and pion decay constant is performed using
mixed-action Lattice QCD calculations with domain-wall valence quarks on
ensembles of rooted, staggered n_f = 2+1 MILC configurations. Calculations were
performed at two lattice spacings of b~0.125 fm and b~0.09 fm, at two strange
quark masses, multiple light quark masses, and a number of lattice volumes. The
ratios of light quark to strange quark masses are in the range 0.1 <= m_l / m_s
<= 0.6, while pion masses are in the range 235 < m_\pi < 680 MeV. A two-flavor
chiral perturbation theory analysis of the Lattice QCD calculations constrains
the Gasser-Leutwyler coefficients bar{l}_3 and bar{l}_4 to be bar{l}_3 =
4.04(40)(+73-55) and bar{l}_4 = 4.30(51)(+84-60). All systematic effects in the
calculations are explored, including those from the finite lattice space-time
volume, the finite lattice spacing, and the finite fifth dimension in the
domain-wall quark action. A consistency is demonstrated between a chiral
perturbation theory analysis at fixed lattice spacing combined with a leading
order continuum extrapolation, and the mixed-action chiral perturbation theory
analysis which explicitly includes the leading order discretization effects.
Chiral corrections to the pion decay constant are found to give f_\pi / f =
1.062(26)(+42-40) where f is the decay constant in the chiral limit. The most
recent scale setting by the MILC Collaboration yields a postdiction of f_\pi =
128.2(3.6)(+4.4-6.0)(+1.2-3.3) MeV at the physical pion mass.
-
The low-energy neutron-Sigma^- interactions determine, in part, the role of
the strange quark in dense matter, such as that found in astrophysical
environments. The scattering phase shifts for this system are obtained from a
numerical evaluation of the QCD path integral using the technique of Lattice
QCD. Our calculations, performed at a pion mass of m_pi ~ 389 MeV in two large
lattice volumes, and at one lattice spacing, are extrapolated to the physical
pion mass using effective field theory. The interactions determined from QCD
are consistent with those extracted from hyperon-nucleon experimental data
within uncertainties, and strengthen theoretical arguments that the strange
quark is a crucial component of dense nuclear matter.
-
The pi+pi+ s-wave scattering phase-shift is determined below the inelastic
threshold using Lattice QCD. Calculations were performed at a pion mass of
m_pi~390 MeV with an anisotropic n_f=2+1 clover fermion discretization in four
lattice volumes, with spatial extent L~2.0, 2.5, 3.0 and 3.9 fm, and with a
lattice spacing of b_s~0.123 fm in the spatial direction and b_t b_s/3.5 in the
time direction. The phase-shift is determined from the energy-eigenvalues of
pi+pi+ systems with both zero and non-zero total momentum in the lattice volume
using Luscher's method. Our calculations are precise enough to allow for a
determination of the threshold scattering parameters, the scattering length a,
the effective range r, and the shape-parameter P, in this channel and to
examine the prediction of two-flavor chiral perturbation theory: m_pi^2 a r =
3+O(m_pi^2/Lambda_chi^2). Chiral perturbation theory is used, with the Lattice
QCD results as input, to predict the scattering phase-shift (and threshold
parameters) at the physical pion mass. Our results are consistent with
determinations from the Roy equations and with the existing experimental phase
shift data.
-
Results of a high-statistics, multi-volume Lattice QCD exploration of the
deuteron, the di-neutron, the H-dibaryon, and the Xi-Xi- system at a pion mass
of m ~ 390 MeV are presented. Calculations were performed with an anisotropic
n_f = 2+1 Clover discretization in four lattice volumes of spatial extent L ~
2.0, 2.5, 3.0 and 4.0 fm, with a lattice spacing of b_s ~ 0.123 fm in the
spatial-direction, and b_t ~ b_s/3.5 in the time-direction. The Xi-Xi- is found
to be bound by B_{Xi-Xi-} = 14.0(1.4)(6.7) MeV, consistent with expectations
based upon phenomenological models and low-energy effective field theories
constrained by nucleon-nucleon and hyperon-nucleon scattering data at the
physical light-quark masses. We find weak evidence that both the deuteron and
the di-neutron are bound at this pion mass, with binding energies of B_d =
11(05)(12) MeV and B_{nn} = 7.1(5.2)(7.3) MeV, respectively. With an increased
number of measurements and a refined analysis, the binding energy of the
H-dibaryon is B_H = 13.2(1.8)(4.0) MeV at this pion mass, updating our previous
result.
-
The current constraints from lattice QCD on the existence of the H-dibaryon
are discussed. With only two significant lattice QCD calculations of the
H-dibaryon binding energy at approximately the same lattice spacing, the forms
of the chiral and continuum extrapolations to the physical point are not
determined. In this brief report, we consider the constraints on the H-dibaryon
imposed by two simple chiral extrapolations. In both instances, the
extrapolation to the physical pion mass allows for a bound H-dibaryon or a
near-threshold scattering state. Further lattice QCD calculations are required
to clarify this situation.
-
Chiral perturbation theory makes definitive predictions for the extrinsic
behavior of hadrons in external electric and magnetic fields. Near the chiral
limit, the electric and magnetic polarizabilities of pions, kaons, and nucleons
are determined in terms of a few well-known parameters. In this limit, hadrons
become quantum mechanically diffuse as polarizabilities scale with the inverse
square-root of the quark mass. In some cases, however, such predictions from
chiral perturbation theory have not compared well with experimental data.
Ultimately we must turn to first principles numerical simulations of QCD to
determine properties of hadrons, and confront the predictions of chiral
perturbation theory. To address the electromagnetic polarizabilities, we
utilize the background field technique. Restricting our attention to
calculations in background electric fields, we demonstrate new techniques to
determine electric polarizabilities and baryon magnetic moments for both
charged and neutral states. As we can study the quark mass dependence of
observables with lattice QCD, the lattice will provide a crucial test of our
understanding of low-energy QCD, which will be timely in light of ongoing
experiments, such as at COMPASS and HI\gamma S.
-
We present evidence for the existence of a bound H-dibaryon, an I=0, J=0,
s=-2 state with valence quark structure uuddss, at a pion mass of m_pi ~ 389
MeV. Using the results of Lattice QCD calculations performed on four ensembles
of anisotropic clover gauge-field configurations, with spatial extents of L ~
2.0, 2.5, 3.0 and 3.9 fm at a spatial lattice spacing of b ~ 0.123 fm, we find
an H-dibaryon bound by B = 16.6 +- 2.1 +- 4.6 MeV at a pion mass of m_pi ~ 389
MeV.
-
The volume dependence of the octet baryon masses and relations among them are
explored with Lattice QCD. Calculations are performed with n_f=2+1 clover
fermion discretization in four lattice volumes, with spatial extent L ~ 2.0,
2.5, 3.0 and 3.9 fm, with an anisotropic lattice spacing of b_s ~ 0.123 fm in
the spatial direction, and b_t = b_s/3.5 in the time direction, and at a pion
mass of m_pi ~ 390 MeV. The typical precision of the ground-state baryon mass
determination is ~0.2%, enabling a precise exploration of the volume dependence
of the masses, the Gell-Mann--Okubo mass relation, and of other mass
combinations. A comparison of the volume dependence with the predictions of
heavy baryon chiral perturbation theory is performed in both the SU(2)_L X
SU(2)_R and SU(3)_L X SU(3)_R expansions. Predictions of the three-flavor
expansion for the hadron masses are found to describe the observed volume
dependences reasonably well. Further, the Delta-N-pi axial coupling constant is
extracted from the volume dependence of the nucleon mass in the two-flavor
expansion, with only small modifications in the three-flavor expansion from the
inclusion of kaons and etas. At a given value of m_pi L, the finite-volume
contributions to the nucleon mass are predicted to be significantly smaller at
m_pi ~ 140 MeV than at m_pi ~ 390 MeV due to a coefficient that scales as ~
m_pi^3. This is relevant for the design of future ensembles of lattice
gauge-field configurations. Finally, the volume dependence of the pion and kaon
masses are analyzed with two-flavor and three-flavor chiral perturbation
theory.
-
We review recent progress toward establishing lattice Quantum Chromodynamics
as a predictive calculational framework for nuclear physics. A survey of the
current techniques that are used to extract low-energy hadronic scattering
amplitudes and interactions is followed by a review of recent two-body and
few-body calculations by the NPLQCD collaboration and others. An outline of the
nuclear physics that is expected to be accomplished with Lattice QCD in the
next decade, along with estimates of the required computational resources, is
presented.
-
We discuss the extraction of the electromagnetic and spin polarisabilities of
nucleons from lattice QCD. We show that the external field method can be used
to measure all the electromagnetic and spin polarisabilities including those of
charged particles. We then turn to the extrapolations required to connect such
calculations to experiment in the context of chiral perturbation theory,
finding a strong dependence on the lattice volume and quark masses.
-
We discuss the extraction of the electromagnetic and spin polarisabilities of
nucleons from lattice QCD. We show that the external field method can be used
to measure all the electromagnetic and spin polarisabilities including those of
charged particles. We then turn to the extrapolations required to connect such
calculations to experiment in the context of finite volume chiral perturbation
theory. We derive results relevant for lattice simulations of QCD,
partially-quenched QCD and quenched QCD. Our results for the polarisabilities
show a strong dependence on the lattice volume and quark masses, typically
differing from the infinite volume limit by ~10% for current lattice volumes
and quark masses.
-
We study the influence of the strong nuclear vector potential, treated using
the mean-field approximation, in deep inelastic scattering. A consistent
treatment of the electromagnetic current operator, combined with the use of the
operator product expansion is presented and discussed.
-
We show that hadronic matrix elements can be extracted from lattice
simulations with background fields that arise from operator exponentiation.
Importantly, flavour-singlet matrix elements can be evaluated without requiring
the computation of disconnected diagrams, thus facilitating a calculation of
the quark contribution to the spin of the proton and the singlet axial
coupling, $g_A^0$. In the two nucleon sector, a background field approach will
allow calculation of the magnetic and quadrupole moments of the deuteron and an
investigation of the EMC effect directly from lattice QCD. Matrix elements
between states of differing momenta are also analysed in the presence of
background fields.
-
The results of different non-perturbative studies agree on a power law as the
infrared behavior of the Landau gauge gluon propagator. This propagator
violates positivity and thus indicates the absence of the transverse gluons
from the physical spectrum, i.e. gluon confinement. A simple analytic structure
for the gluon propagator is proposed capturing all of its features. We comment
also on related investigations for the Landau gauge quark propagator.
-
We derive double distributions for the proton in a simple model that contains
scalar as well as axial-vector diquark correlations. The model parameters are
tuned so that the experimentally measured electromagnetic form factors are
reproduced for small momentum transfer. Resulting generalized parton
distributions satisfy known constraints, including the positivity bounds.
-
We explore the analytic structure of the gluon and quark propagators of
Landau gauge QCD from numerical solutions of the coupled system of renormalized
Dyson--Schwinger equations and from fits to lattice data. We find sizable
negative norm contributions in the transverse gluon propagator indicating the
absence of the transverse gluon from the physical spectrum. A simple analytic
structure for the gluon propagator is proposed. For the quark propagator we
find evidence for a mass-like singularity on the real timelike momentum axis,
with a mass of 350 to 500 MeV. Within the employed Green's functions approach
we identify a crucial term in the quark-gluon vertex that leads to a positive
definite Schwinger function for the quark propagator.
-
In Landau gauge QCD the infrared behavior of the propagator of transverse
gluons can be analytically determined to be a power law from Dyson-Schwinger
equations. This propagator clearly shows positivity violation, indicating the
absence of the transverse gluons from the physical spectrum, i.e. gluon
confinement. A simple analytic structure for the gluon propagator is proposed
capturing all important features. We provide arguments that the Landau gauge
quark propagator possesses a singularity on the real timelike axis. For this
propagator we find a positive definite Schwinger function.