William Detmold
Current Organizations / Institutes
The organizations/institutes from the scholar's latest article
Previous Organizations / Institutes
The organizations / institutes from the scholar's other articles in the past ten years
Research Subjects
Research Interests
William Detmold has a total of 547 co-authors. Here are the top co-authors by number of shared publications:
- Kostas Orginos
Physics Department, William & Mary
26 shared publications - André Walker-Loud
Lawrence Berkeley National Laboratory
21 shared publications - Assumpta Parreno
Departament de Física Quàntica i Astrofísica and Institut de Ciències del Cosmos, Universitat de Barcelona, Martí i Franquès 1, E08028, Spain
20 shared publications NPLQCD Collaboration
15 shared publicationsMartin J. Savage
Department of Physics, University of Washington, Seattle, Washington 98195-1560, USA
14 shared publications- Anthony W. Thomas
University of Adelaide
12 shared publications
Latest Publications
- Open Access
- 19 August 2026
- Published by Association for Computing Machinery (ACM) in ACM Transactions on AI for Science
Matrix-free Neural Preconditioner for the Dirac Operator in Lattice Gauge Theory
Linear systems arise in generating samples and in calculating observables in lattice quantum chromodynamics (QCD). Solving the Hermitian positive definite systems, which are sparse but ill-conditioned, involves using iterative methods, such as Conjugate Gradient (CG), which are time-consuming and computationally expensive. Preconditioners can effectively accelerate this process, with the state-of-the-art being multigrid preconditioners. However, constructing useful preconditioners can be challenging, adding additional computational overhead, especially in large linear systems. We propose a framework, leveraging operator learning techniques, to construct linear maps as effective preconditioners. The method in this work does not rely on explicit matrices from either the original linear systems or the produced preconditioners, allowing efficient model training and application in the CG solver. In the context of the Schwinger model ( \(U(1)\) gauge theory in 1+1 spacetime dimensions with two degenerate-mass fermions), this preconditioning scheme effectively decreases the condition number of the linear systems and approximately halves the number of iterations required for convergence in relevant parameter ranges. We further demonstrate the framework learns a general mapping dependent on the lattice structure which leads to zero-shot learning ability for the Dirac operators constructed from gauge field configurations of different sizes.- 19 August 2026
- Published by Association for Computing Machinery (ACM) in ACM Transactions on AI for Science
Open AccessMatrix-free Neural Preconditioner for the Dirac Operator in Lattice Gauge Theory
Linear systems arise in generating samples and in calculating observables in lattice quantum chromodynamics (QCD). Solving the Hermitian positive definite systems, which are sparse but ill-conditioned, involves using iterative methods, such as Conjugate Gradient (CG), which are time-consuming and computationally expensive. Preconditioners can effectively accelerate this process, with the state-of-the-art being multigrid preconditioners. However, constructing useful preconditioners can be challenging, adding additional computational overhead, especially in large linear systems. We propose a framework, leveraging operator learning techniques, to construct linear maps as effective preconditioners. The method in this work does not rely on explicit matrices from either the original linear systems or the produced preconditioners, allowing efficient model training and application in the CG solver. In the context of the Schwinger model ( \(U(1)\) gauge theory in 1+1 spacetime dimensions with two degenerate-mass fermions), this preconditioning scheme effectively decreases the condition number of the linear systems and approximately halves the number of iterations required for convergence in relevant parameter ranges. We further demonstrate the framework learns a general mapping dependent on the lattice structure which leads to zero-shot learning ability for the Dirac operators constructed from gauge field configurations of different sizes. - Open Access
- 20 April 2026
- Published by American Physical Society (APS) in Physical Review D
- Vol. 113 (7)
Empirical investigation of nuclear correlation function distributions in lattice QCD
Two-point correlation functions of systems with baryon number B ∈ { 1 , 2 , 3 , 4 } are investigated using lattice quantum chromodynamics (QCD). In particular, the empirical distributions of importance-sampling Monte-Carlo samples of these correlation functions are examined as a function of the spacetime separation between the two points and the baryon number. While the exact forms of these distributions are not known for QCD, recent work has determined asymptotic expressions for analogous correlation function distributions in simpler theories such as scalar field theory and the disordered phase of the O ( N ) model. The theoretical O ( N ) model distributions are found to provide an accurate description of the empirical QCD distributions at zero momentum over a wide range of temporal separations for each baryon number when assessed with a range of different statistical tests. In particular, the behavior of the baryon number B QCD correlation function at large temporal separation is well reproduced by the O ( N ∼ 2 / B ) model distribution.
- 20 April 2026
- Published by American Physical Society (APS) in Physical Review D
- Vol. 113 (7)
Open AccessEmpirical investigation of nuclear correlation function distributions in lattice QCD
Two-point correlation functions of systems with baryon number B ∈ { 1 , 2 , 3 , 4 } are investigated using lattice quantum chromodynamics (QCD). In particular, the empirical distributions of importance-sampling Monte-Carlo samples of these correlation functions are examined as a function of the spacetime separation between the two points and the baryon number. While the exact forms of these distributions are not known for QCD, recent work has determined asymptotic expressions for analogous correlation function distributions in simpler theories such as scalar field theory and the disordered phase of the O ( N ) model. The theoretical O ( N ) model distributions are found to provide an accurate description of the empirical QCD distributions at zero momentum over a wide range of temporal separations for each baryon number when assessed with a range of different statistical tests. In particular, the behavior of the baryon number B QCD correlation function at large temporal separation is well reproduced by the O ( N ∼ 2 / B ) model distribution.
- Open Access
- 30 January 2026
- Published by American Physical Society (APS) in Physical Review D
- Vol. 113 (1)
Parton physics from a heavy-quark operator product expansion: Lattice QCD calculation of the fourth moment of the pion distribution amplitude
The pion light-cone distribution amplitude (LCDA) is an essential nonperturbative input for a range of high-energy exclusive processes in quantum chromodynamics. Building on our previous work, the continuum limit of the fourth Mellin moment of the pion LCDA is determined in quenched QCD using quark masses which correspond to a pion mass of m π = 550 MeV . This calculation finds ⟨ ξ 2 ⟩ = 0.202 ( 8 ) ( 9 ) and ⟨ ξ 4 ⟩ = 0.039 ( 28 ) ( 11 ) where the first error indicates the combined statistical and systematic uncertainty from the analysis and the second indicates the uncertainty from working with Wilson coefficients computed to next-to-leading order. These results are presented in the MS ¯ scheme at a renormalization scale of μ = 2 GeV .
- 30 January 2026
- Published by American Physical Society (APS) in Physical Review D
- Vol. 113 (1)
Open AccessParton physics from a heavy-quark operator product expansion: Lattice QCD calculation of the fourth moment of the pion distribution amplitude
The pion light-cone distribution amplitude (LCDA) is an essential nonperturbative input for a range of high-energy exclusive processes in quantum chromodynamics. Building on our previous work, the continuum limit of the fourth Mellin moment of the pion LCDA is determined in quenched QCD using quark masses which correspond to a pion mass of m π = 550 MeV . This calculation finds ⟨ ξ 2 ⟩ = 0.202 ( 8 ) ( 9 ) and ⟨ ξ 4 ⟩ = 0.039 ( 28 ) ( 11 ) where the first error indicates the combined statistical and systematic uncertainty from the analysis and the second indicates the uncertainty from working with Wilson coefficients computed to next-to-leading order. These results are presented in the MS ¯ scheme at a renormalization scale of μ = 2 GeV .
- Open Access
Parton physics from a heavy-quark operator product expansion: Lattice QCD calculation of the fourth moment of the pion distribution amplitude
The pion light-cone distribution amplitude (LCDA) is an essential non-perturbative input for a range of high-energy exclusive processes in quantum chromodynamics. Building on our previous work, the continuum limit of the fourth Mellin moment of the pion LCDA is determined in quenched QCD using quark masses which correspond to a pion mass of $m_π= 550$ MeV. This calculation finds $\langleξ^2\rangle = 0.202(8)(9)$ and $\langle ξ^4 \rangle = 0.039(28)(11)$ where the first error indicates the combined statistical and systematic uncertainty from the analysis and the second indicates the uncertainty from working with Wilson coefficients computed to next-to-leading order. These results are presented in the $\overline{\text{MS}}$ scheme at a renormalization scale of $μ= 2$ GeV.Open AccessParton physics from a heavy-quark operator product expansion: Lattice QCD calculation of the fourth moment of the pion distribution amplitude
The pion light-cone distribution amplitude (LCDA) is an essential non-perturbative input for a range of high-energy exclusive processes in quantum chromodynamics. Building on our previous work, the continuum limit of the fourth Mellin moment of the pion LCDA is determined in quenched QCD using quark masses which correspond to a pion mass of $m_π= 550$ MeV. This calculation finds $\langleξ^2\rangle = 0.202(8)(9)$ and $\langle ξ^4 \rangle = 0.039(28)(11)$ where the first error indicates the combined statistical and systematic uncertainty from the analysis and the second indicates the uncertainty from working with Wilson coefficients computed to next-to-leading order. These results are presented in the $\overline{\text{MS}}$ scheme at a renormalization scale of $μ= 2$ GeV. - Open Access
Empirical investigation of nuclear correlation function distributions in lattice QCD
Two-point correlation functions of systems with baryon number $B \in \{1,2,3,4\}$ are investigated using lattice Quantum Chromodynamics (QCD). In particular, the empirical distributions of importance-sampling Monte-Carlo samples of these correlation functions are examined as a function of the spacetime separation between the two points and the baryon number. While the exact forms of these distributions are not known for QCD, recent work has determined asymptotic expressions for analogous correlation function distributions in simpler theories such as scalar field theory and the disordered phase of the $O(N)$ model. The theoretical O(N) model distributions are found to provide an accurate description of the empirical QCD distributions at zero momentum over a wide range of temporal separations for each baryon number when assessed with a range of different statistical tests. In particular, the behaviour of the baryon number $B$ QCD correlation function at large temporal separation is well-reproduced by the $O(N \sim 2/B)$ model distribution.Open AccessEmpirical investigation of nuclear correlation function distributions in lattice QCD
Two-point correlation functions of systems with baryon number $B \in \{1,2,3,4\}$ are investigated using lattice Quantum Chromodynamics (QCD). In particular, the empirical distributions of importance-sampling Monte-Carlo samples of these correlation functions are examined as a function of the spacetime separation between the two points and the baryon number. While the exact forms of these distributions are not known for QCD, recent work has determined asymptotic expressions for analogous correlation function distributions in simpler theories such as scalar field theory and the disordered phase of the $O(N)$ model. The theoretical O(N) model distributions are found to provide an accurate description of the empirical QCD distributions at zero momentum over a wide range of temporal separations for each baryon number when assessed with a range of different statistical tests. In particular, the behaviour of the baryon number $B$ QCD correlation function at large temporal separation is well-reproduced by the $O(N \sim 2/B)$ model distribution. - Open Access
- 17 June 2025
- Published by American Physical Society (APS) in Physical Review D
- Vol. 111 (11)
Aspects of propagator sparsening in lattice QCD
In lattice field theory, field sparsening aims to replace quantum fields, or objects constructed from them, with approximations that preserve the appropriate symmetries and maintain many aspects of the physics that the fields determine. For example, an effective sparsening of a quark propagator provides an efficient map from a quark propagator on a fine lattice geometry to a quark propagator defined on a coarser geometry in order to reduce storage and computational costs of subsequent calculational stages while maintaining long-distance correlations and corresponding low-energy physical information. Previous studies have focused on decimating lattice sites or randomly sampling lattice sites to reduce the size of the propagator and subsequent costs of Wick contractions. Here, we extend the study of sparsening to incorporate covariant averaging of spatial sites and examine the effects on two-point and three-point correlation functions involving various hadrons. We find that sparsening is most effective in reproducing the unsparsened versions of these correlation functions when weighted covariant-averaging is sequentially applied many times. Published by the American Physical Society 2025
- 17 June 2025
- Published by American Physical Society (APS) in Physical Review D
- Vol. 111 (11)
Open AccessAspects of propagator sparsening in lattice QCD
In lattice field theory, field sparsening aims to replace quantum fields, or objects constructed from them, with approximations that preserve the appropriate symmetries and maintain many aspects of the physics that the fields determine. For example, an effective sparsening of a quark propagator provides an efficient map from a quark propagator on a fine lattice geometry to a quark propagator defined on a coarser geometry in order to reduce storage and computational costs of subsequent calculational stages while maintaining long-distance correlations and corresponding low-energy physical information. Previous studies have focused on decimating lattice sites or randomly sampling lattice sites to reduce the size of the propagator and subsequent costs of Wick contractions. Here, we extend the study of sparsening to incorporate covariant averaging of spatial sites and examine the effects on two-point and three-point correlation functions involving various hadrons. We find that sparsening is most effective in reproducing the unsparsened versions of these correlation functions when weighted covariant-averaging is sequentially applied many times. Published by the American Physical Society 2025
- Open Access
- 2 June 2025
- Published by American Physical Society (APS) in Physical Review D
- Vol. 111 (11)
Constraints on the finite volume two-nucleon spectrum atThe low-energy, finite-volume spectrum of the two-nucleon system at a quark mass corresponding to a pion mass of m π ≈ 806 MeV is studied with lattice quantum chromodynamics (LQCD) using variational methods. The interpolating-operator sets used in [Variational study of two-nucleon systems with lattice QCD, .] are extended by including a complete basis of local hexaquark operators, as well as plane-wave dibaryon operators built from products of both positive- and negative-parity nucleon operators. Results are presented for the isosinglet and isotriplet two-nucleon channels. In both channels, noticeably weaker variational bounds on the lowest few energy eigenvalues are obtained from operator sets which contain only hexaquark operators or operators constructed from the product of two negative-parity nucleons, while other operator sets produce low-energy variational bounds which are consistent within statistical uncertainties. The consequences of these studies for the LQCD understanding of the two-nucleon spectrum are investigated. Published by the American Physical Society 2025
- 2 June 2025
- Published by American Physical Society (APS) in Physical Review D
- Vol. 111 (11)
Open AccessConstraints on the finite volume two-nucleon spectrum atThe low-energy, finite-volume spectrum of the two-nucleon system at a quark mass corresponding to a pion mass of m π ≈ 806 MeV is studied with lattice quantum chromodynamics (LQCD) using variational methods. The interpolating-operator sets used in [Variational study of two-nucleon systems with lattice QCD, .] are extended by including a complete basis of local hexaquark operators, as well as plane-wave dibaryon operators built from products of both positive- and negative-parity nucleon operators. Results are presented for the isosinglet and isotriplet two-nucleon channels. In both channels, noticeably weaker variational bounds on the lowest few energy eigenvalues are obtained from operator sets which contain only hexaquark operators or operators constructed from the product of two negative-parity nucleons, while other operator sets produce low-energy variational bounds which are consistent within statistical uncertainties. The consequences of these studies for the LQCD understanding of the two-nucleon spectrum are investigated. Published by the American Physical Society 2025
- Open Access
Extracting the distribution amplitude of pseudoscalar mesons using the HOPE method
The pseudoscalar meson light-cone distribution amplitudes (LCDAs) are essential non-perturbative inputs for a range of high-energy exclusive processes in quantum chromodynamics. In this proceedings, progress towards a determination of the low Mellin moments of the pion and kaon LCDAs by the HOPE Collaboration is reported.
Open AccessExtracting the distribution amplitude of pseudoscalar mesons using the HOPE method
The pseudoscalar meson light-cone distribution amplitudes (LCDAs) are essential non-perturbative inputs for a range of high-energy exclusive processes in quantum chromodynamics. In this proceedings, progress towards a determination of the low Mellin moments of the pion and kaon LCDAs by the HOPE Collaboration is reported.
- Open Access
- 21 February 2025
- Published by Sissa Medialab Srl
- Vol. 466, 457
Extracting the distribution amplitude of light pseudoscalar mesons using the HOPE method
The pseudoscalar meson light-cone distribution amplitudes (LCDAs) are essential non-perturbative inputs for a range of high-energy exclusive processes in quantum chromodynamics. In this proceedings, progress towards a determination of the low Mellin moments of the pion and kaon LCDAs by the HOPE Collaboration is reported.
- 21 February 2025
- Published by Sissa Medialab Srl
- Vol. 466, 457
Open AccessExtracting the distribution amplitude of light pseudoscalar mesons using the HOPE method
The pseudoscalar meson light-cone distribution amplitudes (LCDAs) are essential non-perturbative inputs for a range of high-energy exclusive processes in quantum chromodynamics. In this proceedings, progress towards a determination of the low Mellin moments of the pion and kaon LCDAs by the HOPE Collaboration is reported.
- Top Cited PapersOpen Access
- 6 January 2025
- Published by American Physical Society (APS) in Physical Review Letters
- Vol. 134 (1) , 011903
QCD Constraints on Isospin-Dense Matter and the Nuclear Equation of State
Understanding the behavior of dense hadronic matter is a central goal in nuclear physics as it governs the nature and dynamics of astrophysical objects such as supernovae and neutron stars. Because of the nonperturbative nature of quantum chromodynamics (QCD), little is known rigorously about hadronic matter in these extreme conditions. Here, lattice QCD calculations are used to compute thermodynamic quantities and the equation of state of QCD over a wide range of isospin chemical potentials with controlled systematic uncertainties. Agreement is seen with chiral perturbation theory when the chemical potential is small. Comparison to perturbative QCD at large chemical potential allows for an estimate of the gap in the superconducting phase, and this quantity is seen to agree with perturbative determinations. Since the partition function for an isospin chemical potential μ I bounds the partition function for a baryon chemical potential μ B = 3 μ I / 2 , these calculations also provide rigorous nonperturbative QCD bounds on the symmetric nuclear matter equation of state over a wide range of baryon densities for the first time. Published by the American Physical Society 2025
- 6 January 2025
- Published by American Physical Society (APS) in Physical Review Letters
- Vol. 134 (1) , 011903
Top Cited PapersOpen AccessQCD Constraints on Isospin-Dense Matter and the Nuclear Equation of State
Understanding the behavior of dense hadronic matter is a central goal in nuclear physics as it governs the nature and dynamics of astrophysical objects such as supernovae and neutron stars. Because of the nonperturbative nature of quantum chromodynamics (QCD), little is known rigorously about hadronic matter in these extreme conditions. Here, lattice QCD calculations are used to compute thermodynamic quantities and the equation of state of QCD over a wide range of isospin chemical potentials with controlled systematic uncertainties. Agreement is seen with chiral perturbation theory when the chemical potential is small. Comparison to perturbative QCD at large chemical potential allows for an estimate of the gap in the superconducting phase, and this quantity is seen to agree with perturbative determinations. Since the partition function for an isospin chemical potential μ I bounds the partition function for a baryon chemical potential μ B = 3 μ I / 2 , these calculations also provide rigorous nonperturbative QCD bounds on the symmetric nuclear matter equation of state over a wide range of baryon densities for the first time. Published by the American Physical Society 2025