Source author record

William Detmold

William Detmold appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

59works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

59 published item(s)

preprint2022arXiv

Finite-Volume Pionless Effective Field Theory for Few-Nucleon Systems with Differentiable Programming

Finite-volume pionless effective field theory provides an efficient framework for the extrapolation of nuclear spectra and matrix elements calculated at finite volume in lattice QCD to infinite volume, and to nuclei with larger atomic number. In this work, it is demonstrated how this framework may be implemented via a set of correlated Gaussian wavefunctions optimised using differentiable programming and via solution of a generalised eigenvalue problem. This approach is shown to be significantly more efficient than a stochastic implementation of the variational method based on the same form of correlated Gaussian wavefunctions, yielding comparably accurate representations of the ground-state wavefunctions with an order of magnitude fewer terms. The efficiency of representation allows such calculations to be extended to larger systems than in previous work. The method is demonstrated through calculations of the binding energies of nuclei with atomic number $A\in\{2,3,4\}$ in finite volume, matched to lattice QCD calculations at quark masses corresponding to $m_π=806$ MeV, and infinite-volume effective field theory calculations of $A\in\{2,3,4,5,6\}$ systems based on this matching.

preprint2022arXiv

Hadron Spectroscopy with Lattice QCD

The status and prospects for investigations of exotic and conventional hadrons with lattice QCD are discussed. The majority of hadrons decay strongly via one or multiple decay-channels, including most of the experimentally discovered exotic hadrons. Despite this difficult challenge, the properties of several hadronic resonances have been determined within lattice QCD. To further discern the spectroscopic properties of various hadrons and to help resolve their nature we present our suggestions for future analytic and lattice studies.

preprint2022arXiv

Novel Algorithms for Computing Correlation Functions of Nuclei

The computational cost required to calculate nuclear correlation functions grows factorially in the number of quarks, making the study of large nuclei inaccessible to ab initio study using lattice QCD at the present time. However, the tensor expressions corresponding to many of these correlation functions exhibit a high degree of permutation symmetry that can be exploited to reduce computational work. We present promising speed-ups for certain choices of interpolating operators using two new algorithms for computing nuclear correlation functions.

preprint2022arXiv

Nuclear Forces for Precision Nuclear Physics -- a collection of perspectives

This is a collection of perspective pieces contributed by the participants of the Institute of Nuclear Theory's Program on Nuclear Physics for Precision Nuclear Physics which was held virtually from April 19 to May 7, 2021. The collection represents the reflections of a vibrant and engaged community of researchers on the status of theoretical research in low-energy nuclear physics, the challenges ahead, and new ideas and strategies to make progress in nuclear structure and reaction physics, effective field theory, lattice QCD, quantum information, and quantum computing. The contributed pieces solely reflect the perspectives of the respective authors and do not represent the viewpoints of the Institute for Nuclear theory or the organizers of the program.

preprint2020arXiv

Lattice QCD Inputs for Nuclear Double Beta Decay

Second order beta-decay processes with and without neutrinos in the final state are key probes of nuclear physics and of the nature of neutrinos. Neutrinoful double-beta decay is the rarest Standard Model process that has been observed and provides a unique test of the understanding of weak nuclear interactions. Observation of neutrinoless double-beta decay would reveal that neutrinos are Majorana fermions and that lepton number conservation is violated in nature. While significant progress has been made in phenomenological approaches to understanding these processes, establishing a connection between these processes and the physics of the Standard Model and beyond is a critical task as it will provide input into the design and interpretation of future experiments. The strong-interaction contributions to double-beta decay processes are non-perturbative and can only be addressed systematically through a combination of lattice Quantum Chromoodynamics (LQCD) and nuclear many-body calculations. In this review, current efforts to establish the LQCD connection are discussed for both neutrinoful and neutrinoless double-beta decay. LQCD calculations of the hadronic contributions to the neutrinoful process $nn\to pp e^- e^- \barν_e\barν_e$ and to various neutrinoless pionic transitions are reviewed, and the connections of these calculations to the phenomenology of double-beta decay through the use of effective field theory (EFTs) is highlighted. At present, LQCD calculations are limited to small nuclear systems, and to pionic subsystems, and require matching to appropriate EFTs to have direct phenomenological impact. However, these calculations have already revealed qualitatively that there are terms in the EFTs that can only be constrained from double-beta decay processes themselves or using inputs from LQCD. Future prospects for direct calculations in larger nuclei are also discussed.

preprint2020arXiv

Path integral contour deformations for noisy observables

Monte Carlo studies of many quantum systems face exponentially severe signal-to-noise problems. We show that noise arising from complex phase fluctuations of observables can be reduced without introducing bias using path integral contour deformation techniques. A numerical study of contour deformations for correlation functions in Abelian gauge theory and complex scalar field theory demonstrates that variance can be reduced by orders of magnitude without modifying Monte Carlo sampling.

preprint2020arXiv

Quarkonium at non-zero isospin density

We calculate the energies of quarkonium bound states in the presence of a medium of nonzero isospin density using lattice QCD. The medium, created using a canonical (fixed isospin charge) approach, induces a reduction of the quarkonium energies. As the isospin density increases, the energy shifts first increase and then saturate. The saturation occurs at an isospin density close to that where previously a qualitative change in the behaviour of the energy density of the medium has been observed, which was conjectured to correspond to a transition from a pion gas to a Bose-Einstein condensed phase. The reduction of the quarkonium energies becomes more pronounced as the heavy-quark mass is decreased, similar to the behaviour seen in two-colour QCD at non-zero quark chemical potential. In the process of our analysis, the $η_b$-$π$ and $Υ$-$π$ scattering phase shifts are determined at low momentum. An interpolation of the scattering lengths to the physical pion mass gives $a_{η_b,π} = 0.0025(8)(6)$ fm and $a_{Υ,π} = 0.0030(9)(7)$ fm.

preprint2020arXiv

Two Nucleon Systems at $m_π\sim 450~{\rm MeV}$ from Lattice QCD

Nucleon-nucleon systems are studied with lattice quantum chromodynamics at a pion mass of $m_π\sim 450~{\rm MeV}$ in three spatial volumes using $n_f=2+1$ flavors of light quarks. At the quark masses employed in this work, the deuteron binding energy is calculated to be $B_d = 14.4^{+3.2}_{-2.6} ~{\rm MeV}$, while the dineutron is bound by $B_{nn} = 12.5^{+3.0}_{-5.0}~{\rm MeV}$. Over the range of energies that are studied, the S-wave scattering phase shifts calculated in the 1S0 and 3S1-3D1 channels are found to be similar to those in nature, and indicate repulsive short-range components of the interactions, consistent with phenomenological nucleon-nucleon interactions. In both channels, the phase shifts are determined at three energies that lie within the radius of convergence of the effective range expansion, allowing for constraints to be placed on the inverse scattering lengths and effective ranges. The extracted phase shifts allow for matching to nuclear effective field theories, from which low energy counterterms are extracted and issues of convergence are investigated. As part of the analysis, a detailed investigation of the single hadron sector is performed, enabling a precise determination of the violation of the Gell-Mann--Okubo mass relation. [An Erratum to the published version is included as an appendix. It details the impact of an error discovered in 2020 and corrects typographical errors.]

preprint2019arXiv

Hadrons and Nuclei

This document is one of a series of whitepapers from the USQCD collaboration. Here, we discuss opportunities for lattice QCD calculations related to the structure and spectroscopy of hadrons and nuclei. An overview of recent lattice calculations of the structure of the proton and other hadrons is presented along with prospects for future extensions. Progress and prospects of hadronic spectroscopy and the study of resonances in the light, strange and heavy quark sectors is summarized. Finally, recent advances in the study of light nuclei from lattice QCD are addressed, and the scope of future investigations that are currently envisioned is outlined.

preprint2019arXiv

Lattice QCD and Neutrino-Nucleus Scattering

This document is one of a series of whitepapers from the USQCD collaboration. Here, we discuss opportunities for lattice QCD in neutrino-oscillation physics, which inevitably entails nucleon and nuclear structure. In addition to discussing pertinent lattice-QCD calculations of nucleon and nuclear matrix elements, the interplay with models of nuclei is discussed. This program of lattice- QCD calculations is relevant to current and upcoming neutrino experiments, becoming increasingly important on the timescale of LBNF/DUNE and HyperK.

preprint2016arXiv

Axial-Current Matrix Elements in Light Nuclei from Lattice QCD

I present results from the first lattice QCD calculations of axial-current matrix elements in light nuclei, performed by the NPLQCD collaboration. Precision calculations of these matrix elements, and the subsequent extraction of multi-nucleon axial-current operators, are essential in refining theoretical predictions of the proton-proton fusion cross section, neutrino-nucleus cross sections and $ββ$-decay rates of nuclei. In addition, they are expected to shed light on the phenomenological quenching of $g_A$ that is required in nuclear many-body calculations.

preprint2016arXiv

Low energy scattering phase shifts for meson-baryon systems

In this work, we calculate meson-baryon scattering phase shifts in four channels using lattice QCD methods. From a set of calculations at four volumes, corresponding to spatial sizes of 2, 2.5, 3, and 4 fm, and a pion mass of m_pi ~ 390 MeV, we determine the scattering lengths and effective ranges for these systems at the corresponding quark masses. We also perform the calculation at a lighter quark mass, m_pi ~ 230 MeV, on the largest volume. Using these determinations, along with those in previous work, we perform a chiral extrapolation of the scattering lengths to the physical point after correcting for the effective range contributions using the multi-volume calculations performed at m_pi ~ 390 MeV.

preprint2016arXiv

Multiscale Monte Carlo equilibration: Two-color QCD with two fermion flavors

We demonstrate the applicability of a recently proposed multiscale thermalization algorithm to two-color quantum chromodynamics (QCD) with two mass-degenerate fermion flavors. The algorithm involves refining an ensemble of gauge configurations that had been generated using a renormalization group (RG) matched coarse action, thereby producing a fine ensemble that is close to the thermalized distribution of a target fine action; the refined ensemble is subsequently rethermalized using conventional algorithms. Although the generalization of this algorithm from pure Yang-Mills theory to QCD with dynamical fermions is straightforward, we find that in the latter case, the method is susceptible to numerical instabilities during the initial stages of rethermalization when using the hybrid Monte Carlo algorithm. We find that these instabilities arise from large fermion forces in the evolution, which are attributed to an accumulation of spurious near-zero modes of the Dirac operator. We propose a simple strategy for curing this problem, and demonstrate that rapid thermalization--as probed by a variety of gluonic and fermionic operators--is possible with the use of this solution. In addition, we study the sensitivity of rethermalization rates to the RG matching of the coarse and fine actions, and identify effective matching conditions based on a variety of measured scales.

preprint2016arXiv

Unitary Limit of Two-Nucleon Interactions in Strong Magnetic Fields

Two-nucleon systems are shown to exhibit large scattering lengths in strong magnetic fields at unphysical quark masses, and the trends toward the physical values indicate that such features may exist in nature. Lattice QCD calculations of the energies of one and two nucleons systems are performed at pion masses of $m_π\sim 450$ and 806 MeV in uniform, time-independent magnetic fields of strength {\bf B}| \sim 10^{19}$-$10^{20}$ Gauss to determine the response of these hadronic systems to large magnetic fields. Fields of this strength may exist inside magnetars and in peripheral relativistic heavy ion collisions, and the unitary behavior at large scattering lengths may have important consequences for these systems.

preprint2015arXiv

$Λ_b \to p \ell^- \barν_\ell$ and $Λ_b \to Λ_c \ell^- \barν_\ell$ form factors from lattice QCD with relativistic heavy quarks

Measurements of the $Λ_b \to p \ell^- \barν_\ell$ and $Λ_b \to Λ_c \ell^- \barν_\ell$ decay rates can be used to determine the magnitudes of the CKM matrix elements $V_{ub}$ and $V_{cb}$, provided that the relevant hadronic form factors are known. Here we present a precise calculation of these form factors using lattice QCD with 2+1 flavors of dynamical domain-wall fermions. The $b$ and $c$ quarks are implemented with relativistic heavy-quark actions, allowing us to work directly at the physical heavy-quark masses. The lattice computation is performed for six different pion masses and two different lattice spacings, using gauge-field configurations generated by the RBC and UKQCD collaborations. The $b \to u$ and $b \to c$ currents are renormalized with a mostly nonperturbative method. We extrapolate the form factor results to the physical pion mass and the continuum limit, parametrizing the $q^2$-dependence using $z$-expansions. The form factors are presented in such a way as to enable the correlated propagation of both statistical and systematic uncertainties into derived quantities such as differential decay rates and asymmetries. Using these form factors, we present predictions for the $Λ_b \to p \ell^- \barν_\ell$ and $Λ_b \to Λ_c \ell^- \barν_\ell$ differential and integrated decay rates. Combined with experimental data, our results enable determinations of $|V_{ub}|$, $|V_{cb}|$, and $|V_{ub}/V_{cb}|$ with theory uncertainties of 4.4%, 2.2%, and 4.9%, respectively.

preprint2015arXiv

Ab initio calculation of the $np \to d γ$ radiative capture process

Lattice QCD calculations of two-nucleon systems are used to isolate the short-distance two-body electromagnetic contributions to the radiative capture process $np \to dγ$, and the photo-disintegration processes $γ^{(\ast)} d \to np$. In nuclear potential models, such contributions are described by phenomenological meson-exchange currents, while in the present work, they are determined directly from the quark and gluon interactions of QCD. Calculations of neutron-proton energy levels in multiple background magnetic fields are performed at two values of the quark masses, corresponding to pion masses of $m_π\sim 450$ and 806 MeV, and are combined with pionless nuclear effective field theory to determine these low-energy inelastic processes. Extrapolating to the physical pion mass, a cross section of $σ^{lqcd}(np\to dγ)=332.4({\tiny \begin{array}{l}+5.4 \\ - 4.7\end{array}})\ mb$ is obtained at an incident neutron speed of $v=2,200\ m/s$, consistent with the experimental value of $σ^{expt}(np \to dγ) = 334.2(0.5)\ mb$.

preprint2015arXiv

Composite Vector Particles in External Electromagnetic Fields

Lattice quantum chromodynamics (QCD) studies of electromagnetic properties of hadrons and light nuclei, such as magnetic moments and polarizabilities, have proven successful with the use of background field methods. With an implementation of nonuniform background electromagnetic fields, properties such as charge radii and higher electromagnetic multipole moments (for states of higher spin) can be additionally obtained. This can be achieved by matching lattice QCD calculations to a corresponding low-energy effective theory that describes the static and quasi-static response of hadrons and nuclei to weak external fields. With particular interest in the case of vector mesons and spin-1 nuclei such as the deuteron, we present an effective field theory of spin-1 particles coupled to external electromagnetic fields. To constrain the charge radius and the electric quadrupole moment of the composite spin-1 field, the single-particle Green's functions in a linearly varying electric field in space are obtained within the effective theory, providing explicit expressions that can be used to match directly onto lattice QCD correlation functions. The viability of an extraction of the charge radius and the electric quadrupole moment of the deuteron from the upcoming lattice QCD calculations of this nucleus is discussed.

preprint2015arXiv

Multiscale Monte Carlo equilibration: Pure Yang-Mills theory

We present a multiscale thermalization algorithm for lattice gauge theory, which enables efficient parallel generation of uncorrelated gauge field configurations. The algorithm combines standard Monte Carlo techniques with ideas drawn from real space renormalization group and multigrid methods. We demonstrate the viability of the algorithm for pure Yang-Mills gauge theory for both heat bath and hybrid Monte Carlo evolution, and show that it ameliorates the problem of topological freezing up to controllable lattice spacing artifacts.

preprint2015arXiv

On the Implementation of General Background Electromagnetic Fields on a Periodic Hypercubic Lattice

Nonuniform background electromagnetic fields, once implemented in lattice quantum chromodynamics calculations of hadronic systems, provide a means to constrain a large class of electromagnetic properties of hadrons and nuclei, from their higher electromagnetic moments and charge radii to their electromagnetic form factors. We show how nonuniform fields can be constructed on a periodic hypercubic lattice under certain conditions and determine the precise form of the background U(1) gauge links that must be imposed on the quantum chromodynamics gauge-field configurations to maintain periodicity. Once supplemented by a set of quantization conditions on the background-field parameters, this construction guarantees that no nonuniformity occurs in the hadronic correlation functions across the boundary of the lattice. The special cases of uniform electric and magnetic fields, a nonuniform electric field that varies linearly in one spatial coordinate (relevant to the determination of quadruple moment and charge radii), nonuniform electric and magnetic fields with given temporal and spatial dependences (relevant to the determination of nucleon spin polarizabilities) and plane-wave electromagnetic fields (relevant to the determination of electromagnetic form factors) are discussed explicitly.

preprint2015arXiv

The Magnetic Structure of Light Nuclei from Lattice QCD

Lattice QCD with background magnetic fields is used to calculate the magnetic moments and magnetic polarizabilities of the nucleons and of light nuclei with $A\le4$, along with the cross-section for the $M1$ transition $np\rightarrow dγ$, at the flavor SU(3)-symmetric point where the pion mass is $m_π\sim 806$ MeV. These magnetic properties are extracted from nucleon and nuclear energies in six uniform magnetic fields of varying strengths. The magnetic moments are presented in a recent Letter. For the charged states, the extraction of the polarizability requires careful treatment of Landau levels, which enter non-trivially in the method that is employed. The nucleon polarizabilities are found to be of similar magnitude to their physical values, with $β_p=5.22(+0.66/-0.45)(0.23) \times 10^{-4}$ fm$^3$ and $β_n=1.253(+0.056/-0.067)(0.055) \times 10^{-4}$ fm$^3$, exhibiting a significant isovector component. The dineutron is bound at these heavy quark masses and its magnetic polarizability, $β_{nn}=1.872(+0.121/-0.113)(0.082) \times 10^{-4}$ fm$^3$ differs significantly from twice that of the neutron. A linear combination of deuteron scalar and tensor polarizabilities is determined by the energies of the $j_z=\pm 1$ deuteron states, and is found to be $β_{d,\pm 1}=4.4(+1.6/-1.5)(0.2) \times 10^{-4}$ fm$^3$. The magnetic polarizabilities of the three-nucleon and four-nucleon systems are found to be positive and similar in size to those of the proton, $β_{^{3}\rm He}=5.4(+2.2/-2.1)(0.2) \times 10^{-4}$ fm$^3$, $β_{^{3}\rm H}=2.6(1.7)(0.1) \times 10^{-4}$ fm$^3$, $β_{^{4}\rm He}=3.4(+2.0/-1.9)(0.2) \times 10^{-4}$ fm$^3$. Mixing between the $j_z=0$ deuteron state and the spin-singlet $np$ state induced by the background magnetic field is used to extract the short-distance two-nucleon counterterm, ${\bar L}_1$, of the pionless effective theory for $NN$ systems (equivalent to the meson-exchange current contribution in nuclear potential models), that dictates the cross-section for the $np\to dγ$ process near threshold. Combined with previous determinations of NN scattering parameters, this enables an ab initio determination of the threshold cross-section at these unphysical masses.

preprint2014arXiv

Charmed bottom baryon spectroscopy from lattice QCD

We calculate the masses of baryons containing one, two, or three heavy quarks using lattice QCD. We consider all possible combinations of charm and bottom quarks, and compute a total of 36 different states with $J^P = \frac12^+$ and $J^P = \frac32^+$. We use domain-wall fermions for the up, down, and strange quarks, a relativistic heavy-quark action for the charm quarks, and nonrelativistic QCD for the bottom quarks. Our analysis includes results from two different lattice spacings and seven different pion masses. We perform extrapolations of the baryon masses to the continuum limit and to the physical pion mass using $SU(4|2)$ heavy-hadron chiral perturbation theory including $1/m_Q$ and finite-volume effects. For the 14 singly heavy baryons that have already been observed, our results agree with the experimental values within the uncertainties. We compare our predictions for the hitherto unobserved states with other lattice calculations and quark-model studies.

preprint2014arXiv

Dark Nuclei I: Cosmology and Indirect Detection

In a companion paper (to be presented), lattice field theory methods are used to show that in two-color, two-flavor QCD there are stable nuclear states in the spectrum. As a commonly studied theory of composite dark matter, this motivates the consideration of possible nuclear physics in this and other composite dark sectors. In this work, early Universe cosmology and indirect detection signatures are explored for both symmetric and asymmetric dark matter, highlighting the unique features that arise from considerations of dark nuclei and associated dark nuclear processes. The present day dark matter abundance may be composed of dark nucleons and/or dark nuclei, where the latter are generated through it dark nucleosynthesis. For symmetric dark matter, indirect detection signatures are possible from annihilation, dark nucleosynthesis, and dark nuclear capture and we present a novel explanation of the galactic center gamma ray excess based on the latter. For asymmetric dark matter, dark nucleosynthesis may alter the capture of dark matter in stars, allowing for captured particles to be processed into nuclei and ejected from the star through dark nucleosynthesis in the core. Notably, dark nucleosynthesis realizes a novel mechanism for indirect detection signals of asymmetric dark matter from regions such as the galactic center, without having to rely on a symmetric dark matter component.

preprint2014arXiv

Dark Nuclei II: Nuclear Spectroscopy in Two-Colour QCD

We consider two-colour QCD with two flavours of quarks as a possible theory of composite dark matter and use lattice field theory methods to investigate nuclear spectroscopy in the spin $J=0$ and $J=1$ multi-baryon sectors. We find compelling evidence that $J=1$ systems with baryon number $B=2,3$ (and their mixed meson-baryon counterparts) are bound states - the analogues of nuclei in this theory. In addition, we estimate the $σ$-terms of the $J=0$ and $J=1$ single baryon states which are important for the coupling of the theory to scalar currents that may mediate interactions with the visible sector.

preprint2014arXiv

Fast Evaluation of Multi-Hadron Correlation Functions

Calculating the values of nuclear correlation functions is computationally intensive due to the fact that the number of terms in a nuclear wave function scales exponentially with atomic number. To speed up this computation, we represent a correlation function as a sum of the determinants of many small matrices, and exploit similarities between the matrices to speed up the calculations of those determinants.

preprint2014arXiv

QCD inequalities for hadron interactions

We derive generalisations of the Weingarten--Witten QCD mass inequalities for particular multi-hadron systems. For systems of any number of identical pseudo-scalar mesons of maximal isospin, these inequalities prove that interactions between the constituent mesons must be repulsive and that no bound states can form in these channels. Similar constraints in less symmetric systems are also extracted. These results are compatible with experimental results (where known) and recent lattice QCD calculations, and also lead to a more stringent bound on the nucleon mass than previously derived, $m_N \ge 3/2 m_π$.

preprint2014arXiv

Signal/noise enhancement strategies for stochastically estimated correlation functions

We develop strategies for enhancing the signal/noise ratio for stochastically sampled correlation functions. The techniques are general and offer a wide range of applicability. We demonstrate the potential of the approach with a generic two-state system, and then explore the practical applicability of the method for single hadron correlators in lattice quantum chromodynamics. In the latter case, we determine the ground state energies of the pion, proton, and delta baryon, as well as the ground and first excited state energy of the rho meson using matrices of correlators computed on an exemplary ensemble of anisotropic gauge configurations. In the majority of cases, we find a modest reduction in the statistical uncertainties on extracted energies compared to conventional variational techniques. However, in the case of the delta baryon, we achieve a factor of three reduction in statistical uncertainties. The variety of outcomes achieved for single hadron correlators illustrates an inherent dependence of the method on the properties of the system under consideration and the operator basis from which the correlators are constructed.

preprint2014arXiv

Signal/noise optimization strategies for stochastically estimated correlation functions

Numerical studies of quantum field theories usually rely upon an accurate determination of stochastically estimated correlation functions in order to extract information about the spectrum of the theory and matrix elements of operators. The reliable determination of such correlators is often hampered by an exponential degradation of signal/noise at late time separations. We demonstrate that it is sometimes possible to achieve significant enhancements of signal/noise by appropriately optimizing correlators with respect to the source and sink interpolating operators, and highlight the large range of possibilities that are available for this task. The ideas are discussed for both a toy model, and single hadron correlators in the context of quantum chromodynamics.

preprint2014arXiv

Uncertainty Quantification in Lattice QCD Calculations for Nuclear Physics

The numerical technique of Lattice QCD holds the promise of connecting the nuclear forces, nuclei, the spectrum and structure of hadrons, and the properties of matter under extreme conditions with the underlying theory of the strong interactions, quantum chromodynamics. A distinguishing, and thus far unique, feature of this formulation is that all of the associated uncertainties, both statistical and systematic can, in principle, be systematically reduced to any desired precision with sufficient computational and human resources. We review the sources of uncertainty inherent in Lattice QCD calculations for nuclear physics, and discuss how each is quantified in current efforts.

preprint2013arXiv

$Λ_b \to p l^- \barν$ form factors from lattice QCD with static b quarks

We present a lattice QCD calculation of form factors for the decay $Λ_b \to p μ^- \barν$, which is a promising channel for determining the CKM matrix element $|V_{ub}|$ at the Large Hadron Collider. In this initial study we work in the limit of static b quarks, where the number of independent form factors reduces to two. We use dynamical domain-wall fermions for the light quarks, and perform the calculation at two different lattice spacings and at multiple values of the light-quark masses in a single large volume. Using our form factor results, we calculate the $Λ_b \to p μ^- \barν$ differential decay rate in the range $14 GeV^2 \leq q^2 \leq q^2_{max}$, and obtain the integral $\int_{14 GeV^2}^{q^2_{max}} [dΓ/dq^2] dq^2 / |V_{ub}|^2 = 15.3 \pm 4.2 ps^{-1}$. Combined with future experimental data, this will give a novel determination of $|V_{ub}|$ with about 15\% theoretical uncertainty. The uncertainty is dominated by the use of the static approximation for the b quark, and can be reduced further by performing the lattice calculation with a more sophisticated heavy-quark action.

preprint2013arXiv

Baryon masses at nonzero isospin/kaon density

We present a lattice QCD calculation of the ground-state energy shifts of various baryons in a medium of pions or kaons at a single value of the quark mass corresponding to a pion mass of m_π~390 MeV and a kaon mass of m_K~540 MeV, and in a spatial volume V (4fm)^3. All systems are created using a canonical formalism in which quark propagators are contracted into correlation functions of fixed isospin/kaon density. We study four different systems, Σ^+(π^+)^n, Ξ^0(π^+)^n, p(K^+)^n, and n(K^+)^n, for up to n=11 mesons. From the ground-state energy shifts we extract two- and three-body scattering parameters, as well as linear combinations of low-energy constants appearing in tree-level chiral perturbation theory.

preprint2013arXiv

Baryon properties in meson mediums from lattice QCD

We present results for the ground-state mass shifts of octet baryons due to the presence of a medium of pions or kaons from a lattice QCD calculation performed at a single value of the quark mass, corresponding to a pion mass of $m_π$ ~ 390 MeV, and a spatial volume V ~ (4fm)^3. We use a canonical approach in which correlators are formed using a single baryon propagator and a fixed number of meson propagators, up to n=9. From the ground-state energies we calculate two- and three-body interaction parameters. We also extract combinations of low-energy constants by comparing our results to tree level chiral perturbation theory at non-zero isospin/kaon chemical potential.

preprint2012arXiv

Axial couplings and strong decay widths of heavy hadrons

We calculate the axial couplings of mesons and baryons containing a heavy quark in the static limit using lattice QCD. These couplings determine the leading interactions in heavy hadron chiral perturbation theory and are central quantities in heavy quark physics, as they control strong decay widths and the light-quark mass dependence of heavy hadron observables. Our analysis makes use of lattice data at six different pion masses, 227 MeV < m_π< 352 MeV, two lattice spacings, a=0.085, 0.112 fm, and a volume of (2.7 fm)^3. Our results for the axial couplings are g_1=0.449(51), g_2=0.84(20), and g_3=0.71(13), where g_1 governs the interaction between heavy-light mesons and pions and g_{2,3} are similar couplings between heavy-light baryons and pions. Using our lattice result for g_3, and constraining 1/m_Q corrections in the strong decay widths with experimental data for Σ_c^{(*)} decays, we obtain Γ[Σ_b^{(*)} \to Λ_b π^\pm] = 4.2(1.0), 4.8(1.1), 7.3(1.6), 7.8(1.8) MeV for the Σ_b^+, Σ_b^-, Σ_b^{*+}, Σ_b^{*-} initial states, respectively. We also derive upper bounds on the widths of the Ξ_b^{\prime(*)} baryons.

preprint2012arXiv

Axial couplings of heavy hadrons from domain-wall lattice QCD

We calculate matrix elements of the axial current for static-light mesons and baryons in lattice QCD with dynamical domain wall fermions. We use partially quenched heavy hadron chiral perturbation theory in a finite volume to extract the axial couplings g_1, g_2, and g_3 from the data. These axial couplings allow the prediction of strong decay rates and enter chiral extrapolations of most lattice results in the b sector. Our calculations are performed with two lattice spacings and with pion masses down to 227 MeV.

preprint2012arXiv

Calculation of the heavy-hadron axial couplings g_1, g_2, and g_3 using lattice QCD

In a recent letter [Phys. Rev. Lett. 108, 172003 (2012), arXiv:1109.2480] we have reported on a lattice QCD calculation of the heavy-hadron axial couplings $g_1$, $g_2$, and $g_3$. These quantities are low-energy constants of heavy-hadron chiral perturbation theory (HH$χ$PT) and are related to the $B^*Bπ$, $Σ_b^*Σ_bπ$, and $Σ_b^{(*)}Λ_bπ$ couplings. In the following, we discuss important details of the calculation and give further results. To determine the axial couplings, we explicitly match the matrix elements of the axial current in QCD with the corresponding matrix elements in HH$χ$PT. We construct the ratios of correlation functions used to calculate the matrix elements in lattice QCD, and study the contributions from excited states. We present the complete numerical results and discuss the data analysis in depth. In particular, we demonstrate the convergence of $SU(4|2)$ HH$χ$PT for the axial-current matrix elements at pion masses up to about 400 MeV and show the impact of the nonanalytic loop contributions. Finally, we present additional predictions for strong and radiative decay widths of charm and bottom baryons.

preprint2012arXiv

Form factors for Lambda_b -> Lambda transitions from lattice QCD

The rare baryonic decays $Λ_b \to Λμ^+ μ^-$ and $Λ_b \to Λγ$ can complement rare $B$ meson decays in constraining models of new physics. In this work, we calculate the relevant $Λ_b \to Λ$ transition form factors at leading order in heavy-quark effective theory using lattice QCD. Our analysis is based on RBC/UKQCD gauge field ensembles with 2+1 flavors of domain-wall fermions, and with lattice spacings of $a\approx 0.11$ fm and $a\approx 0.08$ fm. We compute appropriate ratios of three-point and two-point correlation functions for a wide range of source-sink separations, and extrapolate to infinite separation in order to eliminate excited-state contamination. We then extrapolate the form factors to the continuum limit and to the physical values of the light-quark masses.

preprint2012arXiv

Lambda_b -> Lambda l+ l- form factors and differential branching fraction from lattice QCD

We present the first lattice QCD determination of the $Λ_b \to Λ$ transition form factors that govern the rare baryonic decays $Λ_b \to Λl^+ l^-$ at leading order in heavy-quark effective theory. Our calculations are performed with 2+1 flavors of domain-wall fermions, at two lattice spacings and with pion masses down to 227 MeV. Three-point functions with a wide range of source-sink separations are used to extract the ground-state contributions. The form factors are extrapolated to the physical values of the light-quark masses and to the continuum limit. We use our results to calculate the differential branching fractions for $Λ_b \to Λl^+ l^-$ with $l=e,μ,τ$ within the standard model. We find agreement with a recent CDF measurement of the $Λ_b \to Λμ^+ μ^-$ differential branching fraction.

preprint2012arXiv

Lattice QCD at non-zero isospin chemical potential

Quantum chromodynamics (QCD) at non-zero isospin chemical potential is studied in a canonical approach by analyzing systems of fixed isospin number density. To construct these systems, we develop a range of new algorithms for performing the factorially large numbers of Wick contractions required in multi-hadron systems. We then use these methods to study systems with the quantum numbers of up to 72 $π^+$'s on three ensembles of gauge configurations with spatial extents $L\sim$ 2.0, 2.5 and 3.0 fm, and light quark masses corresponding to a pion mass of {390 MeV}. The ground state energies of these systems are extracted and the volume dependence of these energies is utilized to determine the two- and three- body interactions amongst $π^+$'s. The systems studied correspond to isospin densities of up to $ρ_I\sim 9\ {\rm fm}^{-3}$ and probe isospin chemical potentials, $μ_I$, in the range $m_π \lsim μ_I\ \lsim 4.5\ m_π$, allowing us to investigate aspects of the QCD phase diagram at low temperature and for varying isospin chemical potential. By studying the energy density of the system, we provide numerical evidence for the conjectured transition of the system to a Bose-Einstein condensed phase at $μ_I\ \gsim m_π$.

preprint2012arXiv

Nuclear correlation functions in lattice QCD

We consider the problem of calculating the large number of Wick contractions necessary to compute states with the quantum numbers of many baryons in lattice QCD. We consider a constructive approach and a determinant-based approach and show that these methods allow the required contractions to be performed in computationally manageable amount of time for certain choices of interpolating operators. Examples of correlation functions computed using these techniques are shown for the quantum numbers of the light nuclei, He, Be, C, O and Si.

preprint2011arXiv

Application of Quadrature Methods for Re-Weighting in Lattice QCD

Re-weighting is a useful tool that has been employed in Lattice QCD in different contexts including, tuning the strange quark mass, approaching the light quark mass regime, and simulating electromagnetic fields on top of QCD gauge configurations. In case of re-weighting the sea quark mass, the re-weighting factor is given by the ratio of the determinants of two Dirac operators $D_a$ and $D_b$. A popular approach for computing this ratio is to use a pseudofermion representation of the determinant of the composite operator $Ω=D_a(D_b^\dagger D_b)^{-1} D_a^\dagger$. Here, we study using quadrature methods together with noise vectors to compute the ratio of determinants. We show that, with quadrature methods each determinant can be computed separately using the operators $Ω_a=D_a^\dagger D_a$ and $Ω_b=D_b^\dagger D_b$. We also discuss using bootstrap re-sampling to remove the bias from the determinant estimator.

preprint2011arXiv

Axial couplings in heavy hadron chiral perturbation theory at the next-to-leading order

We present calculations of axial-current matrix elements between various heavy-meson and heavy-baryon states to the next-to-leading order in heavy hadron chiral perturbation theory in the p-regime. When compared with data from lattice computations or experiments, these results can be used to determine the axial couplings in the chiral Lagrangian. Our calculation is performed in partially-quenched chiral perturbation theory for both SU(4|2) and SU(6|3). We incorporate finite-size effects arising from a single Goldstone meson wrapping around the spatial volume. Results for full QCD with two and three flavours can be obtained straightforwardly by taking the sea-quark masses to be equal to the valence-quark masses. To illustrate the impact of our chiral perturbation theory calculation on lattice computations, we analyse the SU(2) full QCD results in detail. We also study one-loop contributions relevant to the heavy hadron strong-decay amplitudes involving final-state Goldstone bosons, and demonstrate that the quark-mass dependence of these amplitudes can be significantly different from that of the axial current matrix elements containing only single hadron external states.

preprint2011arXiv

Investigations of QCD at non-zero isospin density

We investigate the QCD phase diagram as a function of isospin chemical potential at a fixed temperature by directly putting large numbers of π^+s into the system. Correlation functions of N π^+s systems involves N!N! contractions, and become extremely expensive when N is large. In order to alleviate this problem, a recursion relation of correlation functions has been derived in Ref. [1] that substantially reduces the number of independent contractions needed and makes the study of many pions systems be possible. In this proceeding this method is investigated numerically. We have also constructed a new method that is even more efficient, enabling us to study systems of up to 72 π^+s.

preprint2011arXiv

Lattice QCD study of mixed systems of pions and kaons

The different ground state energies of N-pion and M-kaon systems for N+M <=12 are studied in lattice QCD. These energies are then used to extract the various two- and three- body interactions that occur in these systems. Particular attention is paid to additional thermal states present in the spectrum because of the finite temporal extent. These calculations are performed using one ensemble of 2+1 flavor anisotropic lattices with a spatial lattice spacing a_s ~ 0.125 fm, an anisotropy factor ξ=a_s/a_t=3.5, and a spatial volume L^3 ~ (2.5 fm)^3. The quark masses used correspond to pion and kaon masses of m_π ~ 383 MeV and m_K ~ 537 MeV, respectively. The isospin and strangeness chemical potentials of these systems are found to be in the region where chiral perturbation theory and hadronic models predict a phase transition between a pion condensed phase and a kaon condensed phase.

preprint2011arXiv

Medium Effects in Parton Distributions

Understanding the effects of a background hadronic medium on hadronic observables is important in the context of hadron structure. Many experiments probing nucleon structure make use of nuclear targets and unraveling the modifications that ensue is a complex task. Using lattice QCD, we investigate the ab initio computation of hadron structure in a medium, focusing on the structure of the pion in a Bose-condensed gas of pions.

preprint2010arXiv

A method to study complex systems of mesons in Lattice QCD

Finite density systems can be explored with Lattice QCD through the calculation of multi-hadron correlation functions. Recently, systems with up to 12 $π^+$'s or $K^+$'s have been studied to determine the 3-$π^+$ and 3-$K^+$ interactions, and the corresponding chemical potentials have been determined as a function of density. We derive recursion relations between correlation functions that allow this work to be extended to systems of arbitrary numbers of mesons and to systems containing many different types of mesons, such as $π^+$'s, $K^+$'s, $\bar{D}^0$'s and $B^+$'s. These relations allow for the study of finite-density systems in arbitrary volumes, and for the study of high-density systems.

preprint2010arXiv

Extracting Nucleon Magnetic Moments and Electric Polarizabilities from Lattice QCD in Background Electric Fields

Nucleon properties are investigated in background electric fields. As the magnetic moments of baryons affect their relativistic propagation in constant electric fields, electric polarizabilities cannot be determined without knowledge of magnetic moments. This is analogous to the experimental situation, for which determination of polarizabilities from the Compton amplitude requires subtraction of Born terms. With the background field method, we devise combinations of nucleon correlation functions in constant electric fields that isolate magnetic moments and electric polarizabilities. Using an ensemble of anisotropic gauge configurations with dynamical clover fermions, we demonstrate how both observables can be determined from lattice QCD simulations in background electric fields. We obtain results for the neutron and proton, however, our study is currently limited to electrically neutral sea quarks. The value we extract for the nucleon isovector magnetic moment is comparable to those obtained from measuring lattice three-point functions at similar pion masses.

preprint2010arXiv

Nucleon Magnetic Moments and Electric Polarizabilities

Electromagnetic properties of the nucleon are explored with lattice QCD using a novel technique. Focusing on background electric fields, we show how the electric polarizability can be extracted from nucleon correlation functions. A crucial step concerns addressing contributions from the magnetic moment, which affects the relativistic propagation of nucleons in electric fields. By properly handing these contributions, we can determine both magnetic moments and electric polarizabilities. Lattice results from anisotropic clover lattices are presented. Our method is not limited to the neutron; we show results for the proton as well.

preprint2009arXiv

Bottom hadrons from lattice QCD with domain wall and NRQCD fermions

Dynamical 2+1 flavor lattice QCD is used to calculate the masses of bottom hadrons, including B mesons, singly and doubly bottom baryons, and for the first time also the triply-bottom baryon Omega_bbb. The domain wall action is used for the up-, down-, and strange quarks (both valence and sea), while the bottom quark is implemented with non-relativistic QCD. A calculation of the bottomonium spectrum is also presented.

preprint2009arXiv

Extracting Electric Polarizabilities from Lattice QCD

Charged and neutral, pion and kaon electric polarizabilities are extracted from lattice QCD using an ensemble of anisotropic gauge configurations with dynamical clover fermions. We utilize classical background fields to access the polarizabilities from two-point correlation functions. Uniform background fields are achieved by quantizing the electric field strength with the proper treatment of boundary flux. These external fields, however, are implemented only in the valence quark sector. A novel method to extract charge particle polarizabilities is successfully demonstrated for the first time.

preprint2009arXiv

High Statistics Analysis using Anisotropic Clover Lattices: (I) Single Hadron Correlation Functions

We present the results of high-statistics calculations of correlation functions generated with single-baryon interpolating operators on an ensemble of dynamical anisotropic gauge-field configurations generated by the Hadron Spectrum Collaboration using a tadpole-improved clover fermion action and Symanzik-improved gauge action. A total of 292,500 sets of measurements are made using 1194 gauge configurations of size 20^3 x 128 with an anisotropy parameter ξ= b_s/b_t = 3.5, a spatial lattice spacing of b_s=0.1227\pm 0.0008 fm, and pion mass of m_π~ 390 MeV. Ground state baryon masses are extracted with fully quantified uncertainties that are at or below the ~0.2%-level in lattice units. The lowest-lying negative-parity states are also extracted albeit with a somewhat lower level of precision. In the case of the nucleon, this negative-parity state is above the Nπthreshold and, therefore, the isospin-1/2 πN s-wave scattering phase-shift can be extracted using Luescher's method. The disconnected contributions to this process are included indirectly in the gauge-field configurations and do not require additional calculations. The signal-to-noise ratio in the various correlation functions is explored and is found to degrade exponentially faster than naive expectations on many time-slices. This is due to backward propagating states arising from the anti-periodic boundary conditions imposed on the quark-propagators in the time-direction. We explore how best to distribute computational resources between configuration generation and propagator measurements in order to optimize the extraction of single baryon observables.

preprint2009arXiv

High Statistics Analysis using Anisotropic Clover Lattices: (III) Baryon-Baryon Interactions

Low-energy baryon-baryon interactions are calculated in a high-statistics lattice QCD study on a single ensemble of anisotropic clover gauge-field configurations at a pion mass of m_π~ 390 MeV, a spatial volume of L^3 ~ (2.5 fm)^3, and a spatial lattice spacing of b~0.123 fm. Lüscher's method is used to extract nucleon-nucleon, hyperon-nucleon and hyperon-hyperon scattering phase shifts at one momentum from the one- and two-baryon ground-state energies in the lattice volume. The isospin-3/2 NΣinteractions are found to be highly spin-dependent, and the interaction in the ^3S_1 channel is found to be strong. In contrast, the NΛinteractions are found to be spin-independent, within the uncertainties of the calculation, consistent with the absence of one-pion-exchange. The only channel for which a negative energy-shift is found is ΛΛ, indicating that the ΛΛinteraction is attractive, as anticipated from model-dependent discussions regarding the H-dibaryon. The NN scattering lengths are found to be small, clearly indicating the absence of any fine-tuning in the NN-sector at this pion mass. This is consistent with our previous Lattice QCD calculation of NN interactions. The behavior of the signal-to-noise ratio in the baryon-baryon correlation functions, and in the ratio of correlation functions that yields the ground-state energy splitting is explored. In particular, focus is placed on the window of time slices for which the signal-to-noise ratio does not degrade exponentially, as this provides the opportunity to extract quantitative information about multi-baryon systems.

preprint2009arXiv

Meson-Baryon Scattering Lengths from Mixed-Action Lattice QCD

The $π^+Σ^+$, $π^+Ξ^0$, $K^+p$, $K^+n$, and $\overline{K}{}^0 Ξ^0$ scattering lengths are calculated in mixed-action Lattice QCD with domain-wall valence quarks on the asqtad-improved coarse MILC configurations at four light-quark masses, and at two light-quark masses on the fine MILC configurations. Heavy Baryon Chiral Perturbation Theory with two and three flavors of light quarks is used to perform the chiral extrapolations. To the order we work in the three-flavor chiral expansion, the kaon-baryon processes that we investigate show no signs of convergence. Using the two-flavor chiral expansion for extrapolation, the pion-hyperon scattering lengths are found to be $a_{π^+Σ^+}=-0.197\pm0.017$ fm, and $a_{π^+Ξ^0}=-0.098\pm0.017$ fm, where the comprehensive error includes statistical and systematic uncertainties.

preprint2008arXiv

Bottom hadron mass splittings in the static limit from 2+1 flavour lattice QCD

Dynamical 2+1 flavour lattice QCD is used to calculate the splittings between the masses of mesons and baryons containing a single static heavy quark and domain-wall light and strange quarks. Our calculations are based on the dynamical domain-wall gauge field configurations generated by the RBC and UKQCD collaborations at a spatial volume of (2.7 fm)^3 and a range of quark masses with a lightest value corresponding to a (partially-quenched) pion mass of 275 MeV. When extrapolated to the physical values of the light quark masses, the results of our calculations are generally in good agreement with experimental determinations in the bottom sector. However, the static limit splittings between the Omega_b^- baryon and other bottom hadrons tend to slightly underestimate those obtained using the recent D-zero measurement of the Omega_b^-.

preprint2008arXiv

Electric Polarizabilities from Lattice QCD

The response of hadrons to electromagnetic probes is highly constrained by chiral dynamics; but, in some cases, predictions have not compared well with experimental data. The lattice can be used to test the chiral electromagnetism of hadrons and ultimately confront experiment. We use background field techniques to study the electromagnetic polarizabilities of hadrons. Focusing on simulations in background electric fields, we present preliminary results for both charged and neutral particle polarizabilities. The former are extracted using a novel method.

preprint2008arXiv

Multi-meson States in Lattice QCD

In this contribution, I summarise the studies of the properties of Bose-Einstein condensed systems composed of up to twelve pions or kaons carried out by the NPLQCD collaboration. These investigations have provided precise determination the I=2 pi-pi and I=1 K-K scattering lengths and the first determination of three-hadron interactions from QCD, finding a repulsive three-pion interaction of size consistent with naive dimensional analysis and a three kaon interaction consistent with zero. We have also determined the isospin (strangeness) density dependence of the isospin (strangeness) chemical potential, finding results in surprisingly good agreement with chiral perturbation theory.

preprint2007arXiv

Chiral behaviour of Delta B = 2 and Delta C =2 matrix elements

We investigate the light-quark mass and spatial volume dependence of matrix elements of Delta B = 2 and Delta C = 2 four-fermion operators. These operators are relevant for Neutral B(s) and D mixing in and beyond the Standard Model. An important conclusion of this work is that the chiral extrapolations for matrix elements of heavy-light meson mixing beyond the Standard Model are more complicated than that for the Standard Model mixing matrix elements.

preprint2007arXiv

Multi-Pion Systems in Lattice QCD and the Three-Pion Interaction

The ground-state energies of 2, 3, 4 and 5 π^+'s in a spatial volume V (2.5 fm)^3 are computed with lattice QCD. By eliminating the leading contribution from three-π^+ interactions, particular combinations of these n-π^+ ground-state energies provide precise extractions of the π^+π^+ scattering length in agreement with that obtained from calculations involving only two π^+'s. The three-π^+ interaction can be isolated by forming other combinations of the n-π^+ ground-state energies. We find a result that is consistent with a repulsive three-π^+ interaction for m_π< 352 MeV.