Showing 1–50 of 207 results for author: Shanahan, P

QCD constraints on isospin-dense matter and the nuclear equation of state

Authors: Ryan Abbott, William Detmold, Marc Illa, Assumpta Parreño, Robert J. Perry, Fernando Romero-López, Phiala E. Shanahan, Michael L. Wagman

Abstract: 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 non-perturbative nature of quantum chromodynamics (QCD), little is known rigorously about hadronic matter in these extreme conditions. Here, lattice QCD calculations are used to compute thermodyn… ▽ More 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 non-perturbative 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. Agreement is seen with chiral perturbation theory predictions when the chemical potential is small. Comparison to perturbative QCD calculations 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/2μ_I$, these calculations also provide rigorous non-perturbative QCD bounds on the nuclear equation of state over a wide range of baryon densities for the first time. △ Less

Submitted 13 June, 2024; originally announced June 2024.

Comments: 5 pages, 4 figures

Report number: MIT-CTP/5729

Characterization and Novel Application of Power Over Fiber for Electronics in a Harsh Environment

Authors: M. A. Arroyave, B. Behera, F. Cavanna, A. Feld, F. Guo, A. Heindel, C. K. Jung, K. Koch, D. Leon Silverio, D. A. Martinez Caicedo, C. McGrew, A. Paudel, W. Pellico, R. Rivera, J. Rodríguez Rondon, S. Sacerdoti, P. Shanahan, W. Shi, D. Torres Muñoz, D. Totani, C. Uy, C. Vermeulen, H. Vieira de Souza

Abstract: Power-over-Fiber (PoF) technology has been used extensively in settings where high voltages require isolation from ground. In a novel application of PoF, power is provided to photon detector modules located on a surface at $\sim$ 300 kV with respect to ground in the planned DUNE experiment. In cryogenic environments, PoF offers a reliable means of power transmission, leveraging optical fibers to t… ▽ More Power-over-Fiber (PoF) technology has been used extensively in settings where high voltages require isolation from ground. In a novel application of PoF, power is provided to photon detector modules located on a surface at $\sim$ 300 kV with respect to ground in the planned DUNE experiment. In cryogenic environments, PoF offers a reliable means of power transmission, leveraging optical fibers to transfer power with minimal system degradation. PoF technology excels in maintaining low noise levels when delivering power to sensitive electronic systems operating in extreme temperatures and high voltage environments. This paper presents the R$\&$D effort of PoF in extreme conditions and underscores its capacity to revolutionize power delivery and management in critical applications, offering a dependable solution with low noise, optimal efficiency, and superior isolation. △ Less

Submitted 27 May, 2024; originally announced May 2024.

Comments: 29 pages, 32 figures

Report number: FERMILAB-PUB-24-0265-AD-ETD-LBNF

Constraints on the finite volume two-nucleon spectrum at $m_π\approx 806$ MeV

Authors: William Detmold, Marc Illa, William I. Jay, Assumpta Parreño, Robert J. Perry, Phiala E. Shanahan, Michael L. Wagman

Abstract: The low-energy finite-volume spectrum of the two-nucleon system at a pion mass of $m_π\approx 806$ MeV is studied with lattice quantum chromodynamics (LQCD) using variational methods. The interpolating-operator sets used in [Phys.Rev.D 107 (2023) 9, 094508] are extended by including a complete basis of local hexaquark operators, as well as plane-wave dibaryon operators built from products of both… ▽ More The low-energy finite-volume spectrum of the two-nucleon system at a pion mass of $m_π\approx 806$ MeV is studied with lattice quantum chromodynamics (LQCD) using variational methods. The interpolating-operator sets used in [Phys.Rev.D 107 (2023) 9, 094508] 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, these results provide compelling evidence for the presence of an additional state in the low-energy spectrum not present in a non-interacting two-nucleon system. Within the space of operators that are considered in this work, the strongest overlap of this state is onto hexaquark operators, and variational analyses based on operator sets containing only dibaryon operators are insensitive to this state over the temporal extent studied here. The consequences of these studies for the LQCD understanding of the two-nucleon spectrum are investigated. △ Less

Submitted 18 April, 2024; originally announced April 2024.

Comments: 49 pages, 19 Figures

Report number: FERMILAB-PUB-24-0126-T, MIT-CTP/5700

Practical applications of machine-learned flows on gauge fields

Authors: Ryan Abbott, Michael S. Albergo, Denis Boyda, Daniel C. Hackett, Gurtej Kanwar, Fernando Romero-López, Phiala E. Shanahan, Julian M. Urban

Abstract: Normalizing flows are machine-learned maps between different lattice theories which can be used as components in exact sampling and inference schemes. Ongoing work yields increasingly expressive flows on gauge fields, but it remains an open question how flows can improve lattice QCD at state-of-the-art scales. We discuss and demonstrate two applications of flows in replica exchange (parallel tempe… ▽ More Normalizing flows are machine-learned maps between different lattice theories which can be used as components in exact sampling and inference schemes. Ongoing work yields increasingly expressive flows on gauge fields, but it remains an open question how flows can improve lattice QCD at state-of-the-art scales. We discuss and demonstrate two applications of flows in replica exchange (parallel tempering) sampling, aimed at improving topological mixing, which are viable with iterative improvements upon presently available flows. △ Less

Submitted 17 April, 2024; originally announced April 2024.

Comments: 9 pages, 5 figures, proceedings of the 40th International Symposium on Lattice Field Theory (Lattice 2023)

Report number: FERMILAB-CONF-24-0007-T, MIT-CTP/5669

Multiscale Normalizing Flows for Gauge Theories

Authors: Ryan Abbott, Michael S. Albergo, Denis Boyda, Daniel C. Hackett, Gurtej Kanwar, Fernando Romero-López, Phiala E. Shanahan, Julian M. Urban

Abstract: Scale separation is an important physical principle that has previously enabled algorithmic advances such as multigrid solvers. Previous work on normalizing flows has been able to utilize scale separation in the context of scalar field theories, but the principle has been largely unexploited in the context of gauge theories. This work gives an overview of a new method for generating gauge fields u… ▽ More Scale separation is an important physical principle that has previously enabled algorithmic advances such as multigrid solvers. Previous work on normalizing flows has been able to utilize scale separation in the context of scalar field theories, but the principle has been largely unexploited in the context of gauge theories. This work gives an overview of a new method for generating gauge fields using hierarchical normalizing flow models. This method builds gauge fields from the outside in, allowing different parts of the model to focus on different scales of the problem. Numerical results are presented for $U(1)$ and $SU(3)$ gauge theories in 2, 3, and 4 spacetime dimensions. △ Less

Submitted 16 April, 2024; originally announced April 2024.

Comments: Submitted as a proceedings to the 40th International Symposium on Lattice Field Theory (Lattice 2023)

Report number: MIT-CTP/5666,FERMILAB-CONF-24-0008-T

Search for CP-violating Neutrino Non-Standard Interactions with the NOvA Experiment

Authors: NOvA Collaboration, M. A. Acero, B. Acharya, P. Adamson, L. Aliaga, N. Anfimov, A. Antoshkin, E. Arrieta-Diaz, L. Asquith, A. Aurisano, A. Back, N. Balashov, P. Baldi, B. A. Bambah, A. Bat, K. Bays, R. Bernstein, T. J. C. Bezerra, V. Bhatnagar, D. Bhattarai, B. Bhuyan, J. Bian, A. C. Booth, R. Bowles, B. Brahma , et al. (182 additional authors not shown)

Abstract: This Letter reports a search for charge-parity (CP) symmetry violating non-standard interactions (NSI) of neutrinos with matter using the NOvA Experiment, and examines their effects on the determination of the standard oscillation parameters. Data from $ν_μ(\barν_μ)\rightarrowν_μ(\barν_μ)$ and $ν_μ(\barν_μ)\rightarrowν_{e}(\barν_{e})$ oscillation channels are used to measure the effect of the NSI… ▽ More This Letter reports a search for charge-parity (CP) symmetry violating non-standard interactions (NSI) of neutrinos with matter using the NOvA Experiment, and examines their effects on the determination of the standard oscillation parameters. Data from $ν_μ(\barν_μ)\rightarrowν_μ(\barν_μ)$ and $ν_μ(\barν_μ)\rightarrowν_{e}(\barν_{e})$ oscillation channels are used to measure the effect of the NSI parameters $\varepsilon_{eμ}$ and $\varepsilon_{eτ}$. With 90% C.L. the magnitudes of the NSI couplings are constrained to be $|\varepsilon_{eμ}| \, \lesssim 0.3$ and $|\varepsilon_{eτ}| \, \lesssim 0.4$. A degeneracy at $|\varepsilon_{eτ}| \, \approx 1.8$ is reported, and we observe that the presence of NSI limits sensitivity to the standard CP phase $δ_{\tiny\text{CP}}$. △ Less

Submitted 11 March, 2024; originally announced March 2024.

Report number: FERMILAB-PUB-24-0108-PPD

Performance of a modular ton-scale pixel-readout liquid argon time projection chamber

Authors: DUNE Collaboration, A. Abed Abud, B. Abi, R. Acciarri, M. A. Acero, M. R. Adames, G. Adamov, M. Adamowski, D. Adams, M. Adinolfi, C. Adriano, A. Aduszkiewicz, J. Aguilar, B. Aimard, F. Akbar, K. Allison, S. Alonso Monsalve, M. Alrashed, A. Alton, R. Alvarez, T. Alves, H. Amar, P. Amedo, J. Anderson, D. A. Andrade , et al. (1340 additional authors not shown)

Abstract: The Module-0 Demonstrator is a single-phase 600 kg liquid argon time projection chamber operated as a prototype for the DUNE liquid argon near detector. Based on the ArgonCube design concept, Module-0 features a novel 80k-channel pixelated charge readout and advanced high-coverage photon detection system. In this paper, we present an analysis of an eight-day data set consisting of 25 million cosmi… ▽ More The Module-0 Demonstrator is a single-phase 600 kg liquid argon time projection chamber operated as a prototype for the DUNE liquid argon near detector. Based on the ArgonCube design concept, Module-0 features a novel 80k-channel pixelated charge readout and advanced high-coverage photon detection system. In this paper, we present an analysis of an eight-day data set consisting of 25 million cosmic ray events collected in the spring of 2021. We use this sample to demonstrate the imaging performance of the charge and light readout systems as well as the signal correlations between the two. We also report argon purity and detector uniformity measurements, and provide comparisons to detector simulations. △ Less

Submitted 5 March, 2024; originally announced March 2024.

Comments: 47 pages, 41 figures

Report number: FERMILAB-PUB-24-0073-LBNF

Multi-particle interpolating operators in quantum field theories with cubic symmetry

Authors: William Detmold, William I. Jay, Gurtej Kanwar, Phiala E. Shanahan, Michael L. Wagman

Abstract: Numerical studies of lattice quantum field theories are conducted in finite spatial volumes, typically with cubic symmetry in the spatial coordinates. Motivated by these studies, this work presents a general algorithm to construct multi-particle interpolating operators for quantum field theories with cubic symmetry. The algorithm automates the block diagonalization required to combine multiple ope… ▽ More Numerical studies of lattice quantum field theories are conducted in finite spatial volumes, typically with cubic symmetry in the spatial coordinates. Motivated by these studies, this work presents a general algorithm to construct multi-particle interpolating operators for quantum field theories with cubic symmetry. The algorithm automates the block diagonalization required to combine multiple operators of definite linear momentum into irreducible representations of the appropriate little group. Examples are given for distinguishable and indistinguishable particles including cases with both zero and non-zero spin. An implementation of the algorithm is publicly available at https://github.com/latticeqcdtools/mhi. △ Less

Submitted 1 March, 2024; originally announced March 2024.

Comments: 27 pages. An implementation of the algorithm is publicly available at https://github.com/latticeqcdtools/mhi

Report number: MIT-CTP/5685, FERMILAB-PUB-24-0101-T

Long-Distance Nuclear Matrix Elements for Neutrinoless Double-Beta Decay from Lattice QCD

Authors: Zohreh Davoudi, William Detmold, Zhenghao Fu, Anthony V. Grebe, William Jay, David Murphy, Patrick Oare, Phiala E. Shanahan, Michael L. Wagman

Abstract: Neutrinoless double-beta ($0νββ$) decay is a heretofore unobserved process which, if observed, would imply that neutrinos are Majorana particles. Interpretations of the stringent experimental constraints on $0νββ$-decay half-lives require calculations of nuclear matrix elements. This work presents the first lattice quantum-chromodynamics (LQCD) calculation of the matrix element for $0νββ$ decay in… ▽ More Neutrinoless double-beta ($0νββ$) decay is a heretofore unobserved process which, if observed, would imply that neutrinos are Majorana particles. Interpretations of the stringent experimental constraints on $0νββ$-decay half-lives require calculations of nuclear matrix elements. This work presents the first lattice quantum-chromodynamics (LQCD) calculation of the matrix element for $0νββ$ decay in a multi-nucleon system, specifically the $nn \rightarrow pp ee$ transition, mediated by a light left-handed Majorana neutrino propagating over nuclear-scale distances. This calculation is performed with quark masses corresponding to a pion mass of $m_π= 806$ MeV at a single lattice spacing and volume. The statistically cleaner $Σ^- \rightarrow Σ^+ ee$ transition is also computed in order to investigate various systematic uncertainties. The prospects for matching the results of LQCD calculations onto a nuclear effective field theory to determine a leading-order low-energy constant relevant for $0νββ$ decay with a light Majorana neutrino are investigated. This work, therefore, sets the stage for future calculations at physical values of the quark masses that, combined with effective field theory and nuclear many-body studies, will provide controlled theoretical inputs to experimental searches of $0νββ$ decay. △ Less

Submitted 14 February, 2024; originally announced February 2024.

Comments: 16 pages, 7 figures

Report number: FERMILAB-PUB-24-0067-T, MIT-CTP/5682, UMD-PP-024-03

Determination of the Collins-Soper kernel from Lattice QCD

Authors: Artur Avkhadiev, Phiala E. Shanahan, Michael L. Wagman, Yong Zhao

Abstract: This work presents a determination of the quark Collins-Soper kernel, which relates transverse-momentum-dependent parton distributions (TMDs) at different rapidity scales, using lattice quantum chromodynamics (QCD). This is the first such determination with systematic control of quark mass, operator mixing, and discretization effects. Next-to-next-to-leading logarithmic matching is used to match l… ▽ More This work presents a determination of the quark Collins-Soper kernel, which relates transverse-momentum-dependent parton distributions (TMDs) at different rapidity scales, using lattice quantum chromodynamics (QCD). This is the first such determination with systematic control of quark mass, operator mixing, and discretization effects. Next-to-next-to-leading logarithmic matching is used to match lattice-calculable distributions to the corresponding TMDs. The continuum-extrapolated lattice QCD results are consistent with several recent phenomenological parameterizations of the Collins-Soper kernel and are precise enough to disfavor other parameterizations. △ Less

Submitted 9 February, 2024; originally announced February 2024.

Report number: FERMILAB-PUB-24-0037-T, MIT-CTP/5677

Real-time Dynamics of the Schwinger Model as an Open Quantum System with Neural Density Operators

Authors: Joshua Lin, Di Luo, Xiaojun Yao, Phiala E. Shanahan

Abstract: Ab-initio simulations of multiple heavy quarks propagating in a Quark-Gluon Plasma are computationally difficult to perform due to the large dimension of the space of density matrices. This work develops machine learning algorithms to overcome this difficulty by approximating exact quantum states with neural network parametrisations, specifically Neural Density Operators. As a proof of principle d… ▽ More Ab-initio simulations of multiple heavy quarks propagating in a Quark-Gluon Plasma are computationally difficult to perform due to the large dimension of the space of density matrices. This work develops machine learning algorithms to overcome this difficulty by approximating exact quantum states with neural network parametrisations, specifically Neural Density Operators. As a proof of principle demonstration in a QCD-like theory, the approach is applied to solve the Lindblad master equation in the 1+1d lattice Schwinger Model as an open quantum system. Neural Density Operators enable the study of in-medium dynamics on large lattice volumes, where multiple-string interactions and their effects on string-breaking and recombination phenomena can be studied. Thermal properties of the system at equilibrium can also be probed with these methods by variationally constructing the steady state of the Lindblad master equation. Scaling of this approach with system size is studied, and numerical demonstrations on up to 32 spatial lattice sites and with up to 3 interacting strings are performed. △ Less

Submitted 9 February, 2024; originally announced February 2024.

Comments: 23 pages, 11 figures

Report number: MIT-CTP/5679,IQuS@UW-21-073

Doping Liquid Argon with Xenon in ProtoDUNE Single-Phase: Effects on Scintillation Light

Authors: DUNE Collaboration, A. Abed Abud, B. Abi, R. Acciarri, M. A. Acero, M. R. Adames, G. Adamov, M. Adamowski, D. Adams, M. Adinolfi, C. Adriano, A. Aduszkiewicz, J. Aguilar, B. Aimard, F. Akbar, K. Allison, S. Alonso Monsalve, M. Alrashed, A. Alton, R. Alvarez, H. Amar Es-sghir, P. Amedo, J. Anderson, D. A. Andrade, C. Andreopoulos , et al. (1300 additional authors not shown)

Abstract: Doping of liquid argon TPCs (LArTPCs) with a small concentration of xenon is a technique for light-shifting and facilitates the detection of the liquid argon scintillation light. In this paper, we present the results of the first doping test ever performed in a kiloton-scale LArTPC. From February to May 2020, we carried out this special run in the single-phase DUNE Far Detector prototype (ProtoDUN… ▽ More Doping of liquid argon TPCs (LArTPCs) with a small concentration of xenon is a technique for light-shifting and facilitates the detection of the liquid argon scintillation light. In this paper, we present the results of the first doping test ever performed in a kiloton-scale LArTPC. From February to May 2020, we carried out this special run in the single-phase DUNE Far Detector prototype (ProtoDUNE-SP) at CERN, featuring 770 t of total liquid argon mass with 410 t of fiducial mass. The goal of the run was to measure the light and charge response of the detector to the addition of xenon, up to a concentration of 18.8 ppm. The main purpose was to test the possibility for reduction of non-uniformities in light collection, caused by deployment of photon detectors only within the anode planes. Light collection was analysed as a function of the xenon concentration, by using the pre-existing photon detection system (PDS) of ProtoDUNE-SP and an additional smaller set-up installed specifically for this run. In this paper we first summarize our current understanding of the argon-xenon energy transfer process and the impact of the presence of nitrogen in argon with and without xenon dopant. We then describe the key elements of ProtoDUNE-SP and the injection method deployed. Two dedicated photon detectors were able to collect the light produced by xenon and the total light. The ratio of these components was measured to be about 0.65 as 18.8 ppm of xenon were injected. We performed studies of the collection efficiency as a function of the distance between tracks and light detectors, demonstrating enhanced uniformity of response for the anode-mounted PDS. We also show that xenon doping can substantially recover light losses due to contamination of the liquid argon by nitrogen. △ Less

Submitted 9 February, 2024; v1 submitted 2 February, 2024; originally announced February 2024.

Comments: 35 pages, 20 figures

Report number: CERN-EP-2024-024; FERMILAB-PUB-23-0819-LBNF

Applications of flow models to the generation of correlated lattice QCD ensembles

Authors: Ryan Abbott, Aleksandar Botev, Denis Boyda, Daniel C. Hackett, Gurtej Kanwar, Sébastien Racanière, Danilo J. Rezende, Fernando Romero-López, Phiala E. Shanahan, Julian M. Urban

Abstract: Machine-learned normalizing flows can be used in the context of lattice quantum field theory to generate statistically correlated ensembles of lattice gauge fields at different action parameters. This work demonstrates how these correlations can be exploited for variance reduction in the computation of observables. Three different proof-of-concept applications are demonstrated using a novel residu… ▽ More Machine-learned normalizing flows can be used in the context of lattice quantum field theory to generate statistically correlated ensembles of lattice gauge fields at different action parameters. This work demonstrates how these correlations can be exploited for variance reduction in the computation of observables. Three different proof-of-concept applications are demonstrated using a novel residual flow architecture: continuum limits of gauge theories, the mass dependence of QCD observables, and hadronic matrix elements based on the Feynman-Hellmann approach. In all three cases, it is shown that statistical uncertainties are significantly reduced when machine-learned flows are incorporated as compared with the same calculations performed with uncorrelated ensembles or direct reweighting. △ Less

Submitted 28 May, 2024; v1 submitted 19 January, 2024; originally announced January 2024.

Comments: 12 pages, 2 tables, 5 figures. v2: accepted for publication

Report number: MIT-CTP/5658, FERMILAB-PUB-24-0014-T

The DUNE Far Detector Vertical Drift Technology, Technical Design Report

Authors: DUNE Collaboration, A. Abed Abud, B. Abi, R. Acciarri, M. A. Acero, M. R. Adames, G. Adamov, M. Adamowski, D. Adams, M. Adinolfi, C. Adriano, A. Aduszkiewicz, J. Aguilar, B. Aimard, F. Akbar, K. Allison, S. Alonso Monsalve, M. Alrashed, A. Alton, R. Alvarez, H. Amar, P. Amedo, J. Anderson, D. A. Andrade, C. Andreopoulos , et al. (1304 additional authors not shown)

Abstract: DUNE is an international experiment dedicated to addressing some of the questions at the forefront of particle physics and astrophysics, including the mystifying preponderance of matter over antimatter in the early universe. The dual-site experiment will employ an intense neutrino beam focused on a near and a far detector as it aims to determine the neutrino mass hierarchy and to make high-precisi… ▽ More DUNE is an international experiment dedicated to addressing some of the questions at the forefront of particle physics and astrophysics, including the mystifying preponderance of matter over antimatter in the early universe. The dual-site experiment will employ an intense neutrino beam focused on a near and a far detector as it aims to determine the neutrino mass hierarchy and to make high-precision measurements of the PMNS matrix parameters, including the CP-violating phase. It will also stand ready to observe supernova neutrino bursts, and seeks to observe nucleon decay as a signature of a grand unified theory underlying the standard model. The DUNE far detector implements liquid argon time-projection chamber (LArTPC) technology, and combines the many tens-of-kiloton fiducial mass necessary for rare event searches with the sub-centimeter spatial resolution required to image those events with high precision. The addition of a photon detection system enhances physics capabilities for all DUNE physics drivers and opens prospects for further physics explorations. Given its size, the far detector will be implemented as a set of modules, with LArTPC designs that differ from one another as newer technologies arise. In the vertical drift LArTPC design, a horizontal cathode bisects the detector, creating two stacked drift volumes in which ionization charges drift towards anodes at either the top or bottom. The anodes are composed of perforated PCB layers with conductive strips, enabling reconstruction in 3D. Light-trap-style photon detection modules are placed both on the cryostat's side walls and on the central cathode where they are optically powered. This Technical Design Report describes in detail the technical implementations of each subsystem of this LArTPC that, together with the other far detector modules and the near detector, will enable DUNE to achieve its physics goals. △ Less

Submitted 5 December, 2023; originally announced December 2023.

Comments: 425 pages; 281 figures Central editing team: A. Heavey, S. Kettell, A. Marchionni, S. Palestini, S. Rajogopalan, R. J. Wilson

Report number: Fermilab Report no: TM-2813-LBNF

Expanding neutrino oscillation parameter measurements in NOvA using a Bayesian approach

Authors: NOvA Collaboration, M. A. Acero, B. Acharya, P. Adamson, N. Anfimov, A. Antoshkin, E. Arrieta-Diaz, L. Asquith, A. Aurisano, A. Back, N. Balashov, P. Baldi, B. A. Bambah, A. Bat, K. Bays, R. Bernstein, T. J. C. Bezerra, V. Bhatnagar, D. Bhattarai, B. Bhuyan, J. Bian, A. C. Booth, R. Bowles, B. Brahma, C. Bromberg , et al. (174 additional authors not shown)

Abstract: NOvA is a long-baseline neutrino oscillation experiment that measures oscillations in charged-current $ν_μ \rightarrow ν_μ$ (disappearance) and $ν_μ \rightarrow ν_{e}$ (appearance) channels, and their antineutrino counterparts, using neutrinos of energies around 2 GeV over a distance of 810 km. In this work we reanalyze the dataset first examined in our previous paper [Phys. Rev. D 106, 032004 (20… ▽ More NOvA is a long-baseline neutrino oscillation experiment that measures oscillations in charged-current $ν_μ \rightarrow ν_μ$ (disappearance) and $ν_μ \rightarrow ν_{e}$ (appearance) channels, and their antineutrino counterparts, using neutrinos of energies around 2 GeV over a distance of 810 km. In this work we reanalyze the dataset first examined in our previous paper [Phys. Rev. D 106, 032004 (2022)] using an alternative statistical approach based on Bayesian Markov Chain Monte Carlo. We measure oscillation parameters consistent with the previous results. We also extend our inferences to include the first NOvA measurements of the reactor mixing angle $θ_{13}$ and the Jarlskog invariant. We use these results to quantify the strength of our inferences about CP violation, as well as to examine the effects of constraints from short-baseline measurements of $θ_{13}$ using antineutrinos from nuclear reactors when making NOvA measurements of $θ_{23}$. Our long-baseline measurement of $θ_{13}$ is also shown to be consistent with the reactor measurements, supporting the general applicability and robustness of the PMNS framework for neutrino oscillations. △ Less

Submitted 27 May, 2024; v1 submitted 13 November, 2023; originally announced November 2023.

Comments: 20 pages, 17 figures; version accepted by Phys. Rev. D. Data associated with this paper is available at https://doi.org/10.15484/2349444

Report number: FERMILAB-PUB-23-667-AD-CSAID-ND

Journal ref: Phys.Rev.D 110 (2024) 1, 012005

Gravitational form factors of the proton from lattice QCD

Authors: Daniel C. Hackett, Dimitra A. Pefkou, Phiala E. Shanahan

Abstract: The gravitational form factors (GFFs) of a hadron encode fundamental aspects of its structure, including its shape and size as defined from e.g., its energy density. This work presents a determination of the flavor decomposition of the GFFs of the proton from lattice QCD, in the kinematic region $0\leq -t\leq 2~\text{GeV}^2$. The decomposition into up-, down-, strange-quark, and gluon contribution… ▽ More The gravitational form factors (GFFs) of a hadron encode fundamental aspects of its structure, including its shape and size as defined from e.g., its energy density. This work presents a determination of the flavor decomposition of the GFFs of the proton from lattice QCD, in the kinematic region $0\leq -t\leq 2~\text{GeV}^2$. The decomposition into up-, down-, strange-quark, and gluon contributions provides first-principles constraints on the role of each constituent in generating key proton structure observables, such as its mechanical radius, mass radius, and $D$-term. △ Less

Submitted 23 June, 2024; v1 submitted 12 October, 2023; originally announced October 2023.

Comments: Version published in PRL

Report number: MIT-CTP/5630,FERMILAB-PUB-23-592-T

Journal ref: Physical Review Letters, Vol. 132, Iss. 25 (2024)

Advances in machine-learning-based sampling motivated by lattice quantum chromodynamics

Authors: Kyle Cranmer, Gurtej Kanwar, Sébastien Racanière, Danilo J. Rezende, Phiala E. Shanahan

Abstract: Sampling from known probability distributions is a ubiquitous task in computational science, underlying calculations in domains from linguistics to biology and physics. Generative machine-learning (ML) models have emerged as a promising tool in this space, building on the success of this approach in applications such as image, text, and audio generation. Often, however, generative tasks in scienti… ▽ More Sampling from known probability distributions is a ubiquitous task in computational science, underlying calculations in domains from linguistics to biology and physics. Generative machine-learning (ML) models have emerged as a promising tool in this space, building on the success of this approach in applications such as image, text, and audio generation. Often, however, generative tasks in scientific domains have unique structures and features -- such as complex symmetries and the requirement of exactness guarantees -- that present both challenges and opportunities for ML. This Perspective outlines the advances in ML-based sampling motivated by lattice quantum field theory, in particular for the theory of quantum chromodynamics. Enabling calculations of the structure and interactions of matter from our most fundamental understanding of particle physics, lattice quantum chromodynamics is one of the main consumers of open-science supercomputing worldwide. The design of ML algorithms for this application faces profound challenges, including the necessity of scaling custom ML architectures to the largest supercomputers, but also promises immense benefits, and is spurring a wave of development in ML-based sampling more broadly. In lattice field theory, if this approach can realize its early promise it will be a transformative step towards first-principles physics calculations in particle, nuclear and condensed matter physics that are intractable with traditional approaches. △ Less

Submitted 3 September, 2023; originally announced September 2023.

Comments: 11 pages, 5 figures

Journal ref: Nature Reviews Physics 5, 526-535 (2023)

Signal-to-noise improvement through neural network contour deformations for 3D $SU(2)$ lattice gauge theory

Authors: William Detmold, Gurtej Kanwar, Yin Lin, Phiala E. Shanahan, Michael L. Wagman

Abstract: Complex contour deformations of the path integral have been demonstrated to significantly improve the signal-to-noise ratio of observables in previous studies of two-dimensional gauge theories with open boundary conditions. In this work, new developments based on gauge fixing and a neural network definition of the deformation are introduced, which enable an effective application to theories in hig… ▽ More Complex contour deformations of the path integral have been demonstrated to significantly improve the signal-to-noise ratio of observables in previous studies of two-dimensional gauge theories with open boundary conditions. In this work, new developments based on gauge fixing and a neural network definition of the deformation are introduced, which enable an effective application to theories in higher dimensions and with generic boundary conditions. Improvements of the signal-to-noise ratio by up to three orders of magnitude for Wilson loop measurements are shown in $SU(2)$ lattice gauge theory in three spacetime dimensions. △ Less

Submitted 1 September, 2023; originally announced September 2023.

Comments: 9 pages, 3 figures. Proceedings for the 40th Lattice conference at Fermilab from July 31 to August 4, 2023

Lattice quantum chromodynamics at large isospin density: 6144 pions in a box

Authors: Ryan Abbott, William Detmold, Fernando Romero-López, Zohreh Davoudi, Marc Illa, Assumpta Parreño, Robert J. Perry, Phiala E. Shanahan, Michael L. Wagman

Abstract: We present an algorithm to compute correlation functions for systems with the quantum numbers of many identical mesons from lattice quantum chromodynamics (QCD). The algorithm is numerically stable and allows for the computation of $n$-pion correlation functions for $n \in \{ 1, \dots, N\}$ using a single $N \times N$ matrix decomposition, improving on previous algorithms. We apply the algorithm t… ▽ More We present an algorithm to compute correlation functions for systems with the quantum numbers of many identical mesons from lattice quantum chromodynamics (QCD). The algorithm is numerically stable and allows for the computation of $n$-pion correlation functions for $n \in \{ 1, \dots, N\}$ using a single $N \times N$ matrix decomposition, improving on previous algorithms. We apply the algorithm to calculations of correlation functions with up to 6144 $π^+$s using two ensembles of gauge field configurations generated with quark masses corresponding to a pion mass $m_π= 170$ MeV and spacetime volumes of $(4.4^3\times 8.8)\ {\rm fm}^4$ and $(5.8^3\times 11.6)\ {\rm fm}^4$. We also discuss statistical techniques for the analysis of such systems, in which the correlation functions vary over many orders of magnitude. In particular, we observe that the many-pion correlation functions are well approximated by log-normal distributions, allowing the extraction of the energies of these systems. Using these energies, the large-isospin-density, zero-baryon-density region of the QCD phase diagram is explored. A peak is observed in the energy density at an isospin chemical potential $μ_I\sim 1.5 m_π$, signalling the transition into a Bose-Einstein condensed phase. The isentropic speed of sound in the medium is seen to exceed the ideal-gas (conformal) limit ($c_s^2\leq 1/3$) over a wide range of chemical potential before falling towards the asymptotic expectation at $μ_I\sim 15 m_π$. These, and other thermodynamic observables, indicate that the isospin chemical potential must be large for the system to be well described by an ideal gas or perturbative QCD. △ Less

Submitted 27 July, 2023; originally announced July 2023.

Comments: 16 pages, 18 figures, 1 table

Report number: MIT-CTP/5560,UMD-PP-023-03,FERMILAB-PUB-23-382-T

Collins-Soper kernel from lattice QCD at the physical pion mass

Authors: Artur Avkhadiev, Phiala Shanahan, Michael Wagman, Yong Zhao

Abstract: This work presents a determination of the quark Collins-Soper kernel, which relates transverse-momentum-dependent parton distributions (TMDs) at different rapidity scales, using lattice Quantum Chromodynamics (QCD). This is the first lattice QCD calculation of the kernel at quark masses corresponding to a close-to-physical value of the pion mass, with next-to-next-leading logarithmic matching to T… ▽ More This work presents a determination of the quark Collins-Soper kernel, which relates transverse-momentum-dependent parton distributions (TMDs) at different rapidity scales, using lattice Quantum Chromodynamics (QCD). This is the first lattice QCD calculation of the kernel at quark masses corresponding to a close-to-physical value of the pion mass, with next-to-next-leading logarithmic matching to TMDs from the corresponding lattice-calculable distributions, and includes a complete analysis of systematic uncertainties arising from operator mixing. The kernel is extracted at transverse momentum scales $240\,\text{MeV}\lesssim q_{T}\lesssim 1.6\,\text{GeV}$ with a precision sufficient to begin to discriminate between different phenomenological models in the non-perturbative region. △ Less

Submitted 11 January, 2024; v1 submitted 23 July, 2023; originally announced July 2023.

Comments: 52 pages, 47 figures, 3 tables

Report number: MIT-CTP/5587, FERMILAB-PUB-23-375-T

Journal ref: Phys. Rev. D 108, 114505 (2023)

Gravitational form factors of the pion from lattice QCD

Authors: Daniel C. Hackett, Patrick R. Oare, Dimitra A. Pefkou, Phiala E. Shanahan

Abstract: The two gravitational form factors of the pion, $A^π(t)$ and $D^π(t)$, are computed as functions of the momentum transfer squared $t$ in the kinematic region $0\leq -t< 2~\text{GeV}^2$ on a lattice QCD ensemble with quark masses corresponding to a close-to-physical pion mass $m_π\approx 170~\text{MeV}$ and $N_f=2+1$ quark flavors. The flavor decomposition of these form factors into gluon, up/down… ▽ More The two gravitational form factors of the pion, $A^π(t)$ and $D^π(t)$, are computed as functions of the momentum transfer squared $t$ in the kinematic region $0\leq -t< 2~\text{GeV}^2$ on a lattice QCD ensemble with quark masses corresponding to a close-to-physical pion mass $m_π\approx 170~\text{MeV}$ and $N_f=2+1$ quark flavors. The flavor decomposition of these form factors into gluon, up/down light-quark, and strange-quark contributions is presented in the $\overline{\text{MS}}$ scheme at energy scale $μ=2~\text{GeV}$, with renormalization factors computed nonperturbatively via the RI-MOM scheme. Using monopole and $z$-expansion fits to the gravitational form factors, we obtain estimates for the pion momentum fraction and $D$-term that are consistent with the momentum fraction sum rule and the next-to-leading order chiral perturbation theory prediction for $D^π(0)$. △ Less

Submitted 11 December, 2023; v1 submitted 21 July, 2023; originally announced July 2023.

Comments: 30 pages, 18 figures, 7 tables; version accepted for publication at PRD

Report number: MIT-CTP/5585

The case for an EIC Theory Alliance: Theoretical Challenges of the EIC

Authors: Raktim Abir, Igor Akushevich, Tolga Altinoluk, Daniele Paolo Anderle, Fatma P. Aslan, Alessandro Bacchetta, Baha Balantekin, Joao Barata, Marco Battaglieri, Carlos A. Bertulani, Guillaume Beuf, Chiara Bissolotti, Daniël Boer, M. Boglione, Radja Boughezal, Eric Braaten, Nora Brambilla, Vladimir Braun, Duane Byer, Francesco Giovanni Celiberto, Yang-Ting Chien, Ian C. Cloët, Martha Constantinou, Wim Cosyn, Aurore Courtoy , et al. (146 additional authors not shown)

Abstract: We outline the physics opportunities provided by the Electron Ion Collider (EIC). These include the study of the parton structure of the nucleon and nuclei, the onset of gluon saturation, the production of jets and heavy flavor, hadron spectroscopy and tests of fundamental symmetries. We review the present status and future challenges in EIC theory that have to be addressed in order to realize thi… ▽ More We outline the physics opportunities provided by the Electron Ion Collider (EIC). These include the study of the parton structure of the nucleon and nuclei, the onset of gluon saturation, the production of jets and heavy flavor, hadron spectroscopy and tests of fundamental symmetries. We review the present status and future challenges in EIC theory that have to be addressed in order to realize this ambitious and impactful physics program, including how to engage a diverse and inclusive workforce. In order to address these many-fold challenges, we propose a coordinated effort involving theory groups with differing expertise is needed. We discuss the scientific goals and scope of such an EIC Theory Alliance. △ Less

Submitted 23 May, 2023; originally announced May 2023.

Comments: 44 pages, ReVTeX, White Paper on EIC Theory Alliance

Constraint of pionless EFT using two-nucleon spectra from lattice QCD

Authors: William Detmold, Fernando Romero-López, Phiala E. Shanahan

Abstract: Finite-volume pionless effective field theory (FVEFT$_{ π\!/ }$) at next-to-leading order (NLO) is used to analyze the two-nucleon lattice QCD spectrum of Ref.~\cite{Amarasinghe:2021lqa}, performed at quark masses corresponding to a pion mass of approximately $800 $ MeV. Specifically, the effective theory is formulated in finite volume, and variational sets of wave functions are optimized using di… ▽ More Finite-volume pionless effective field theory (FVEFT$_{ π\!/ }$) at next-to-leading order (NLO) is used to analyze the two-nucleon lattice QCD spectrum of Ref.~\cite{Amarasinghe:2021lqa}, performed at quark masses corresponding to a pion mass of approximately $800 $ MeV. Specifically, the effective theory is formulated in finite volume, and variational sets of wave functions are optimized using differential programming. Using these wave functions projected to the appropriate finite-volume symmetry group, variational bounds from FVEFT$_{π\!/ }$ are obtained for the ground state, as well as excited states. By comparison with the lattice QCD GEVP spectrum, different low energy constants (LECs) are constrained. Relativistic corrections are incorporated, allowing for the extractions of NLO LECs, as well as the leading $s$-$d$-wave mixing term in the deuteron channel. △ Less

Submitted 10 May, 2023; originally announced May 2023.

Comments: 17 pages, 8 figs, 2 tables

Report number: MIT-CTP/5538

Normalizing flows for lattice gauge theory in arbitrary space-time dimension

Authors: Ryan Abbott, Michael S. Albergo, Aleksandar Botev, Denis Boyda, Kyle Cranmer, Daniel C. Hackett, Gurtej Kanwar, Alexander G. D. G. Matthews, Sébastien Racanière, Ali Razavi, Danilo J. Rezende, Fernando Romero-López, Phiala E. Shanahan, Julian M. Urban

Abstract: Applications of normalizing flows to the sampling of field configurations in lattice gauge theory have so far been explored almost exclusively in two space-time dimensions. We report new algorithmic developments of gauge-equivariant flow architectures facilitating the generalization to higher-dimensional lattice geometries. Specifically, we discuss masked autoregressive transformations with tracta… ▽ More Applications of normalizing flows to the sampling of field configurations in lattice gauge theory have so far been explored almost exclusively in two space-time dimensions. We report new algorithmic developments of gauge-equivariant flow architectures facilitating the generalization to higher-dimensional lattice geometries. Specifically, we discuss masked autoregressive transformations with tractable and unbiased Jacobian determinants, a key ingredient for scalable and asymptotically exact flow-based sampling algorithms. For concreteness, results from a proof-of-principle application to SU(3) lattice gauge theory in four space-time dimensions are reported. △ Less

Submitted 3 May, 2023; originally announced May 2023.

TMD Handbook

Authors: Renaud Boussarie, Matthias Burkardt, Martha Constantinou, William Detmold, Markus Ebert, Michael Engelhardt, Sean Fleming, Leonard Gamberg, Xiangdong Ji, Zhong-Bo Kang, Christopher Lee, Keh-Fei Liu, Simonetta Liuti, Thomas Mehen, Andreas Metz, John Negele, Daniel Pitonyak, Alexei Prokudin, Jian-Wei Qiu, Abha Rajan, Marc Schlegel, Phiala Shanahan, Peter Schweitzer, Iain W. Stewart, Andrey Tarasov , et al. (4 additional authors not shown)

Abstract: This handbook provides a comprehensive review of transverse-momentum-dependent parton distribution functions and fragmentation functions, commonly referred to as transverse momentum distributions (TMDs). TMDs describe the distribution of partons inside the proton and other hadrons with respect to both their longitudinal and transverse momenta. They provide unique insight into the internal momentum… ▽ More This handbook provides a comprehensive review of transverse-momentum-dependent parton distribution functions and fragmentation functions, commonly referred to as transverse momentum distributions (TMDs). TMDs describe the distribution of partons inside the proton and other hadrons with respect to both their longitudinal and transverse momenta. They provide unique insight into the internal momentum and spin structure of hadrons, and are a key ingredient in the description of many collider physics cross sections. Understanding TMDs requires a combination of theoretical techniques from quantum field theory, nonperturbative calculations using lattice QCD, and phenomenological analysis of experimental data. The handbook covers a wide range of topics, from theoretical foundations to experimental analyses, as well as recent developments and future directions. It is intended to provide an essential reference for researchers and graduate students interested in understanding the structure of hadrons and the dynamics of partons in high energy collisions. △ Less

Submitted 6 April, 2023; originally announced April 2023.

Comments: 471 pages, many figures

Report number: JLAB-THY-23-3780, LA-UR-21-20798, MIT-CTP/5386

Impact of cross-section uncertainties on supernova neutrino spectral parameter fitting in the Deep Underground Neutrino Experiment

Authors: DUNE Collaboration, A. Abed Abud, B. Abi, R. Acciarri, M. A. Acero, M. R. Adames, G. Adamov, M. Adamowski, D. Adams, M. Adinolfi, C. Adriano, A. Aduszkiewicz, J. Aguilar, Z. Ahmad, J. Ahmed, B. Aimard, F. Akbar, K. Allison, S. Alonso Monsalve, M. Alrashed, A. Alton, R. Alvarez, P. Amedo, J. Anderson, D. A. Andrade , et al. (1294 additional authors not shown)

Abstract: A primary goal of the upcoming Deep Underground Neutrino Experiment (DUNE) is to measure the $\mathcal{O}(10)$ MeV neutrinos produced by a Galactic core-collapse supernova if one should occur during the lifetime of the experiment. The liquid-argon-based detectors planned for DUNE are expected to be uniquely sensitive to the $ν_e$ component of the supernova flux, enabling a wide variety of physics… ▽ More A primary goal of the upcoming Deep Underground Neutrino Experiment (DUNE) is to measure the $\mathcal{O}(10)$ MeV neutrinos produced by a Galactic core-collapse supernova if one should occur during the lifetime of the experiment. The liquid-argon-based detectors planned for DUNE are expected to be uniquely sensitive to the $ν_e$ component of the supernova flux, enabling a wide variety of physics and astrophysics measurements. A key requirement for a correct interpretation of these measurements is a good understanding of the energy-dependent total cross section $σ(E_ν)$ for charged-current $ν_e$ absorption on argon. In the context of a simulated extraction of supernova $ν_e$ spectral parameters from a toy analysis, we investigate the impact of $σ(E_ν)$ modeling uncertainties on DUNE's supernova neutrino physics sensitivity for the first time. We find that the currently large theoretical uncertainties on $σ(E_ν)$ must be substantially reduced before the $ν_e$ flux parameters can be extracted reliably: in the absence of external constraints, a measurement of the integrated neutrino luminosity with less than 10\% bias with DUNE requires $σ(E_ν)$ to be known to about 5%. The neutrino spectral shape parameters can be known to better than 10% for a 20% uncertainty on the cross-section scale, although they will be sensitive to uncertainties on the shape of $σ(E_ν)$. A direct measurement of low-energy $ν_e$-argon scattering would be invaluable for improving the theoretical precision to the needed level. △ Less

Submitted 7 July, 2023; v1 submitted 29 March, 2023; originally announced March 2023.

Comments: 25 pages, 21 figures

Report number: FERMILAB-PUB-23-132-CSAID-LBNF-ND-T

Journal ref: Phys. Rev. D 107, 112012 (2023)

Colloquium: Gravitational Form Factors of the Proton

Authors: V. D. Burkert, L. Elouadrhiri, F. X. Girod, C. Lorce, P. Schweitzer, P. E. Shanahan

Abstract: The physics of the gravitational form factors of the proton, and their understanding within quantum chromodynamics, has advanced significantly in the past two decades through both theory and experiment. This Colloquium provides an overview of this progress, highlights the physical insights unveiled by studies of gravitational form factors, and reviews their interpretation in terms of the mechanica… ▽ More The physics of the gravitational form factors of the proton, and their understanding within quantum chromodynamics, has advanced significantly in the past two decades through both theory and experiment. This Colloquium provides an overview of this progress, highlights the physical insights unveiled by studies of gravitational form factors, and reviews their interpretation in terms of the mechanical properties of the proton. △ Less

Submitted 16 January, 2024; v1 submitted 14 March, 2023; originally announced March 2023.

Comments: 22 pages, 15 figures

The Present and Future of QCD

Authors: P. Achenbach, D. Adhikari, A. Afanasev, F. Afzal, C. A. Aidala, A. Al-bataineh, D. K. Almaalol, M. Amaryan, D. Androić, W. R. Armstrong, M. Arratia, J. Arrington, A. Asaturyan, E. C. Aschenauer, H. Atac, H. Avakian, T. Averett, C. Ayerbe Gayoso, X. Bai, K. N. Barish, N. Barnea, G. Basar, M. Battaglieri, A. A. Baty, I. Bautista , et al. (378 additional authors not shown)

Abstract: This White Paper presents the community inputs and scientific conclusions from the Hot and Cold QCD Town Meeting that took place September 23-25, 2022 at MIT, as part of the Nuclear Science Advisory Committee (NSAC) 2023 Long Range Planning process. A total of 424 physicists registered for the meeting. The meeting highlighted progress in Quantum Chromodynamics (QCD) nuclear physics since the 2015… ▽ More This White Paper presents the community inputs and scientific conclusions from the Hot and Cold QCD Town Meeting that took place September 23-25, 2022 at MIT, as part of the Nuclear Science Advisory Committee (NSAC) 2023 Long Range Planning process. A total of 424 physicists registered for the meeting. The meeting highlighted progress in Quantum Chromodynamics (QCD) nuclear physics since the 2015 LRP (LRP15) and identified key questions and plausible paths to obtaining answers to those questions, defining priorities for our research over the coming decade. In defining the priority of outstanding physics opportunities for the future, both prospects for the short (~ 5 years) and longer term (5-10 years and beyond) are identified together with the facilities, personnel and other resources needed to maximize the discovery potential and maintain United States leadership in QCD physics worldwide. This White Paper is organized as follows: In the Executive Summary, we detail the Recommendations and Initiatives that were presented and discussed at the Town Meeting, and their supporting rationales. Section 2 highlights major progress and accomplishments of the past seven years. It is followed, in Section 3, by an overview of the physics opportunities for the immediate future, and in relation with the next QCD frontier: the EIC. Section 4 provides an overview of the physics motivations and goals associated with the EIC. Section 5 is devoted to the workforce development and support of diversity, equity and inclusion. This is followed by a dedicated section on computing in Section 6. Section 7 describes the national need for nuclear data science and the relevance to QCD research. △ Less

Submitted 4 March, 2023; originally announced March 2023.

Comments: QCD Town Meeting White Paper, as submitted to 2023 NSAC LRP committee on Feb. 28, 2023

Journal ref: Nucl.Phys.A 1047 (2024) 122874

Highly-parallelized simulation of a pixelated LArTPC on a GPU

Authors: DUNE Collaboration, A. Abed Abud, B. Abi, R. Acciarri, M. A. Acero, M. R. Adames, G. Adamov, M. Adamowski, D. Adams, M. Adinolfi, C. Adriano, A. Aduszkiewicz, J. Aguilar, Z. Ahmad, J. Ahmed, B. Aimard, F. Akbar, K. Allison, S. Alonso Monsalve, M. Alrashed, C. Alt, A. Alton, R. Alvarez, P. Amedo, J. Anderson , et al. (1282 additional authors not shown)

Abstract: The rapid development of general-purpose computing on graphics processing units (GPGPU) is allowing the implementation of highly-parallelized Monte Carlo simulation chains for particle physics experiments. This technique is particularly suitable for the simulation of a pixelated charge readout for time projection chambers, given the large number of channels that this technology employs. Here we pr… ▽ More The rapid development of general-purpose computing on graphics processing units (GPGPU) is allowing the implementation of highly-parallelized Monte Carlo simulation chains for particle physics experiments. This technique is particularly suitable for the simulation of a pixelated charge readout for time projection chambers, given the large number of channels that this technology employs. Here we present the first implementation of a full microphysical simulator of a liquid argon time projection chamber (LArTPC) equipped with light readout and pixelated charge readout, developed for the DUNE Near Detector. The software is implemented with an end-to-end set of GPU-optimized algorithms. The algorithms have been written in Python and translated into CUDA kernels using Numba, a just-in-time compiler for a subset of Python and NumPy instructions. The GPU implementation achieves a speed up of four orders of magnitude compared with the equivalent CPU version. The simulation of the current induced on $10^3$ pixels takes around 1 ms on the GPU, compared with approximately 10 s on the CPU. The results of the simulation are compared against data from a pixel-readout LArTPC prototype. △ Less

Submitted 28 February, 2023; v1 submitted 19 December, 2022; originally announced December 2022.

Comments: 26 pages, 15 figures

Report number: FERMILAB-PUB-22-926-LBNF

Aspects of scaling and scalability for flow-based sampling of lattice QCD

Authors: Ryan Abbott, Michael S. Albergo, Aleksandar Botev, Denis Boyda, Kyle Cranmer, Daniel C. Hackett, Alexander G. D. G. Matthews, Sébastien Racanière, Ali Razavi, Danilo J. Rezende, Fernando Romero-López, Phiala E. Shanahan, Julian M. Urban

Abstract: Recent applications of machine-learned normalizing flows to sampling in lattice field theory suggest that such methods may be able to mitigate critical slowing down and topological freezing. However, these demonstrations have been at the scale of toy models, and it remains to be determined whether they can be applied to state-of-the-art lattice quantum chromodynamics calculations. Assessing the vi… ▽ More Recent applications of machine-learned normalizing flows to sampling in lattice field theory suggest that such methods may be able to mitigate critical slowing down and topological freezing. However, these demonstrations have been at the scale of toy models, and it remains to be determined whether they can be applied to state-of-the-art lattice quantum chromodynamics calculations. Assessing the viability of sampling algorithms for lattice field theory at scale has traditionally been accomplished using simple cost scaling laws, but as we discuss in this work, their utility is limited for flow-based approaches. We conclude that flow-based approaches to sampling are better thought of as a broad family of algorithms with different scaling properties, and that scalability must be assessed experimentally. △ Less

Submitted 14 November, 2022; originally announced November 2022.

Comments: 22 pages, 8 figures

Report number: MIT-CTP/5496

Identification and reconstruction of low-energy electrons in the ProtoDUNE-SP detector

Authors: DUNE Collaboration, A. Abed Abud, B. Abi, R. Acciarri, M. A. Acero, M. R. Adames, G. Adamov, M. Adamowski, D. Adams, M. Adinolfi, C. Adriano, A. Aduszkiewicz, J. Aguilar, Z. Ahmad, J. Ahmed, B. Aimard, F. Akbar, K. Allison, S. Alonso Monsalve, M. Alrashed, C. Alt, A. Alton, R. Alvarez, P. Amedo, J. Anderson , et al. (1235 additional authors not shown)

Abstract: Measurements of electrons from $ν_e$ interactions are crucial for the Deep Underground Neutrino Experiment (DUNE) neutrino oscillation program, as well as searches for physics beyond the standard model, supernova neutrino detection, and solar neutrino measurements. This article describes the selection and reconstruction of low-energy (Michel) electrons in the ProtoDUNE-SP detector. ProtoDUNE-SP is… ▽ More Measurements of electrons from $ν_e$ interactions are crucial for the Deep Underground Neutrino Experiment (DUNE) neutrino oscillation program, as well as searches for physics beyond the standard model, supernova neutrino detection, and solar neutrino measurements. This article describes the selection and reconstruction of low-energy (Michel) electrons in the ProtoDUNE-SP detector. ProtoDUNE-SP is one of the prototypes for the DUNE far detector, built and operated at CERN as a charged particle test beam experiment. A sample of low-energy electrons produced by the decay of cosmic muons is selected with a purity of 95%. This sample is used to calibrate the low-energy electron energy scale with two techniques. An electron energy calibration based on a cosmic ray muon sample uses calibration constants derived from measured and simulated cosmic ray muon events. Another calibration technique makes use of the theoretically well-understood Michel electron energy spectrum to convert reconstructed charge to electron energy. In addition, the effects of detector response to low-energy electron energy scale and its resolution including readout electronics threshold effects are quantified. Finally, the relation between the theoretical and reconstructed low-energy electron energy spectrum is derived and the energy resolution is characterized. The low-energy electron selection presented here accounts for about 75% of the total electron deposited energy. After the addition of lost energy using a Monte Carlo simulation, the energy resolution improves from about 40% to 25% at 50~MeV. These results are used to validate the expected capabilities of the DUNE far detector to reconstruct low-energy electrons. △ Less

Submitted 31 May, 2023; v1 submitted 2 November, 2022; originally announced November 2022.

Comments: 19 pages, 10 figures

Report number: FERMILAB-PUB-22-784, CERN-EP-DRAFT-MISC-2022-008

Journal ref: Phys. Rev. D 107, 092012 (2023)

The Future of High Energy Physics Software and Computing

Authors: V. Daniel Elvira, Steven Gottlieb, Oliver Gutsche, Benjamin Nachman, S. Bailey, W. Bhimji, P. Boyle, G. Cerati, M. Carrasco Kind, K. Cranmer, G. Davies, V. D. Elvira, R. Gardner, K. Heitmann, M. Hildreth, W. Hopkins, T. Humble, M. Lin, P. Onyisi, J. Qiang, K. Pedro, G. Perdue, A. Roberts, M. Savage, P. Shanahan , et al. (3 additional authors not shown)

Abstract: Software and Computing (S&C) are essential to all High Energy Physics (HEP) experiments and many theoretical studies. The size and complexity of S&C are now commensurate with that of experimental instruments, playing a critical role in experimental design, data acquisition/instrumental control, reconstruction, and analysis. Furthermore, S&C often plays a leading role in driving the precision of th… ▽ More Software and Computing (S&C) are essential to all High Energy Physics (HEP) experiments and many theoretical studies. The size and complexity of S&C are now commensurate with that of experimental instruments, playing a critical role in experimental design, data acquisition/instrumental control, reconstruction, and analysis. Furthermore, S&C often plays a leading role in driving the precision of theoretical calculations and simulations. Within this central role in HEP, S&C has been immensely successful over the last decade. This report looks forward to the next decade and beyond, in the context of the 2021 Particle Physics Community Planning Exercise ("Snowmass") organized by the Division of Particles and Fields (DPF) of the American Physical Society. △ Less

Submitted 8 November, 2022; v1 submitted 11 October, 2022; originally announced October 2022.

Comments: Computational Frontier Report Contribution to Snowmass 2021; 41 pages, 1 figure. v2: missing ref and added missing topical group conveners. v3: fixed typos

Report of the Snowmass 2021 Topical Group on Lattice Gauge Theory

Authors: Zohreh Davoudi, Ethan T. Neil, Christian W. Bauer, Tanmoy Bhattacharya, Thomas Blum, Peter Boyle, Richard C. Brower, Simon Catterall, Norman H. Christ, Vincenzo Cirigliano, Gilberto Colangelo, Carleton DeTar, William Detmold, Robert G. Edwards, Aida X. El-Khadra, Steven Gottlieb, Rajan Gupta, Daniel C. Hackett, Anna Hasenfratz, Taku Izubuchi, William I. Jay, Luchang Jin, Christopher Kelly, Andreas S. Kronfeld, Christoph Lehner , et al. (13 additional authors not shown)

Abstract: Lattice gauge theory continues to be a powerful theoretical and computational approach to simulating strongly interacting quantum field theories, whose applications permeate almost all disciplines of modern-day research in High-Energy Physics. Whether it is to enable precision quark- and lepton-flavor physics, to uncover signals of new physics in nucleons and nuclei, to elucidate hadron structure… ▽ More Lattice gauge theory continues to be a powerful theoretical and computational approach to simulating strongly interacting quantum field theories, whose applications permeate almost all disciplines of modern-day research in High-Energy Physics. Whether it is to enable precision quark- and lepton-flavor physics, to uncover signals of new physics in nucleons and nuclei, to elucidate hadron structure and spectrum, to serve as a numerical laboratory to reach beyond the Standard Model, or to invent and improve state-of-the-art computational paradigms, the lattice-gauge-theory program is in a prime position to impact the course of developments and enhance discovery potential of a vibrant experimental program in High-Energy Physics over the coming decade. This projection is based on abundant successful results that have emerged using lattice gauge theory over the years: on continued improvement in theoretical frameworks and algorithmic suits; on the forthcoming transition into the exascale era of high-performance computing; and on a skillful, dedicated, and organized community of lattice gauge theorists in the U.S. and worldwide. The prospects of this effort in pushing the frontiers of research in High-Energy Physics have recently been studied within the U.S. decadal Particle Physics Planning Exercise (Snowmass 2021), and the conclusions are summarized in this Topical Report. △ Less

Submitted 21 September, 2022; originally announced September 2022.

Comments: 57 pages, 1 figure. Submitted to the Proceedings of the US Community Study on the Future of Particle Physics (Snowmass 2021). Topical Group Report for TF05 - Lattice Gauge Theory

Report number: UMD-PP-022-08, LA-UR-22-29361, FERMILAB-CONF-22-703-T

Snowmass 2021 Computational Frontier CompF03 Topical Group Report: Machine Learning

Authors: Phiala Shanahan, Kazuhiro Terao, Daniel Whiteson

Abstract: The rapidly-developing intersection of machine learning (ML) with high-energy physics (HEP) presents both opportunities and challenges to our community. Far beyond applications of standard ML tools to HEP problems, genuinely new and potentially revolutionary approaches are being developed by a generation of talent literate in both fields. There is an urgent need to support the needs of the interdi… ▽ More The rapidly-developing intersection of machine learning (ML) with high-energy physics (HEP) presents both opportunities and challenges to our community. Far beyond applications of standard ML tools to HEP problems, genuinely new and potentially revolutionary approaches are being developed by a generation of talent literate in both fields. There is an urgent need to support the needs of the interdisciplinary community driving these developments, including funding dedicated research at the intersection of the two fields, investing in high-performance computing at universities and tailoring allocation policies to support this work, developing of community tools and standards, and providing education and career paths for young researchers attracted by the intellectual vitality of machine learning for high energy physics. △ Less

Submitted 15 September, 2022; originally announced September 2022.

Comments: Contribution to Snowmass 2021

Strategies for quantum-optimized construction of interpolating operators in classical simulations of lattice quantum field theories

Authors: A. Avkhadiev, P. E. Shanahan, R. D. Young

Abstract: It has recently been argued that noisy intermediate-scale quantum computers may be used to optimize interpolating operator constructions for lattice quantum field theory (LQFT) calculations on classical computers. Here, two concrete realizations of the method are developed and implemented. The first approach is to maximize the overlap, or fidelity, of the state created by an interpolating operator… ▽ More It has recently been argued that noisy intermediate-scale quantum computers may be used to optimize interpolating operator constructions for lattice quantum field theory (LQFT) calculations on classical computers. Here, two concrete realizations of the method are developed and implemented. The first approach is to maximize the overlap, or fidelity, of the state created by an interpolating operator acting on the vacuum state to the target eigenstate. The second is to instead minimize the energy expectation value of the interpolated state. These approaches are implemented in a proof-of-concept calculation in (1+1)-dimensions for a single-flavor massive Schwinger model to obtain quantum-optimized interpolating operator constructions for a vector meson state in the theory. Although fidelity maximization is preferable in the absence of noise due to quantum gate errors, it is found that energy minimization is more robust to these effects in the proof-of-concept calculation. This work serves as a concrete demonstration of how quantum computers in the intermediate term might be used to accelerate classical LQFT calculations. △ Less

Submitted 2 September, 2022; originally announced September 2022.

Comments: 13 pages, 5 figures

Report number: MIT-CTP/5459; ADP-22-26/T1197

On the nonorientable 4-genus of double twist knots

Authors: Jim Hoste, Patrick D. Shanahan, Cornelia A. Van Cott

Abstract: We investigate the nonorientable 4-genus $γ_4$ of a special family of 2-bridge knots, the twist knots and double twist knots $C(m,n)$. Because the nonorientable 4-genus is bounded by the nonorientable 3-genus, it is known that $γ_4(C(m,n)) \le 3$. By using explicit constructions to obtain upper bounds on $γ_4$ and known obstructions derived from Donaldson's diagonalization theorem to obtain lower… ▽ More We investigate the nonorientable 4-genus $γ_4$ of a special family of 2-bridge knots, the twist knots and double twist knots $C(m,n)$. Because the nonorientable 4-genus is bounded by the nonorientable 3-genus, it is known that $γ_4(C(m,n)) \le 3$. By using explicit constructions to obtain upper bounds on $γ_4$ and known obstructions derived from Donaldson's diagonalization theorem to obtain lower bounds on $γ_4$, we produce infinite subfamilies of $C(m,n)$ where $γ_4=0,1,2,$ and $3$, respectively. However, there remain infinitely many double twist knots where our work only shows that $γ_4$ lies in one of the sets $\{1,2\}, \{2,3\}$, or $\{1,2,3\}$. We tabulate our results for all $C(m,n)$ with $|m|$ and $|n|$ up to 50. We also provide an infinite number of examples which answer a conjecture of Murakami and Yasuhara. △ Less

Submitted 28 March, 2023; v1 submitted 16 August, 2022; originally announced August 2022.

Comments: Some exposition is revised, a figure is added, and typos are corrected, following comments from the referee

MSC Class: 57K10

Neutrinoless Double Beta Decay from Lattice QCD: The Short-Distance $π^-\rightarrowπ^+ e^- e^-$ Amplitude

Authors: William Detmold, William I. Jay, David J. Murphy, Patrick R. Oare, Phiala E. Shanahan

Abstract: This work presents a determination of potential short-distance contributions to the unphysical $π^-\rightarrowπ^+ e^- e^-$ decay through lattice QCD calculations. The hadronic contributions to the transition amplitude are described by the pion matrix elements of five Standard Model Effective Field Theory operators, which are computed on five ensembles of domain-wall fermions with $N_f = 2 + 1$ qua… ▽ More This work presents a determination of potential short-distance contributions to the unphysical $π^-\rightarrowπ^+ e^- e^-$ decay through lattice QCD calculations. The hadronic contributions to the transition amplitude are described by the pion matrix elements of five Standard Model Effective Field Theory operators, which are computed on five ensembles of domain-wall fermions with $N_f = 2 + 1$ quark flavors with a range of heavier-than-physical values of the light quark masses. The matrix elements are extrapolated to the continuum, physical light-quark mass, and infinite volume limit using a functional form derived in chiral Effective Field Theory ($χ\mathrm{EFT}$). This extrapolation also yields the relevant low-energy constants of $χ\mathrm{EFT}$, which are necessary input for $χ\mathrm{EFT}$ calculations of neutrinoless double beta decay of nuclei. △ Less

Submitted 15 December, 2022; v1 submitted 10 August, 2022; originally announced August 2022.

Comments: 24 pages, 17 figures

Report number: MIT-CTP/5414

Sampling QCD field configurations with gauge-equivariant flow models

Authors: Ryan Abbott, Michael S. Albergo, Aleksandar Botev, Denis Boyda, Kyle Cranmer, Daniel C. Hackett, Gurtej Kanwar, Alexander G. D. G. Matthews, Sébastien Racanière, Ali Razavi, Danilo J. Rezende, Fernando Romero-López, Phiala E. Shanahan, Julian M. Urban

Abstract: Machine learning methods based on normalizing flows have been shown to address important challenges, such as critical slowing-down and topological freezing, in the sampling of gauge field configurations in simple lattice field theories. A critical question is whether this success will translate to studies of QCD. This Proceedings presents a status update on advances in this area. In particular, it… ▽ More Machine learning methods based on normalizing flows have been shown to address important challenges, such as critical slowing-down and topological freezing, in the sampling of gauge field configurations in simple lattice field theories. A critical question is whether this success will translate to studies of QCD. This Proceedings presents a status update on advances in this area. In particular, it is illustrated how recently developed algorithmic components may be combined to construct flow-based sampling algorithms for QCD in four dimensions. The prospects and challenges for future use of this approach in at-scale applications are summarized. △ Less

Submitted 20 August, 2022; v1 submitted 7 August, 2022; originally announced August 2022.

Comments: Submitted as a proceedings to the 39th International Symposium on Lattice Field Theory (Lattice 2022)

Neural-network preconditioners for solving the Dirac equation in lattice gauge theory

Authors: Salvatore Calì, Daniel C. Hackett, Yin Lin, Phiala E. Shanahan, Brian Xiao

Abstract: This work develops neural-network--based preconditioners to accelerate solution of the Wilson-Dirac normal equation in lattice quantum field theories. The approach is implemented for the two-flavor lattice Schwinger model near the critical point. In this system, neural-network preconditioners are found to accelerate the convergence of the conjugate gradient solver compared with the solution of unp… ▽ More This work develops neural-network--based preconditioners to accelerate solution of the Wilson-Dirac normal equation in lattice quantum field theories. The approach is implemented for the two-flavor lattice Schwinger model near the critical point. In this system, neural-network preconditioners are found to accelerate the convergence of the conjugate gradient solver compared with the solution of unpreconditioned systems or those preconditioned with conventional approaches based on even-odd or incomplete Cholesky decompositions, as measured by reductions in the number of iterations and/or complex operations required for convergence. It is also shown that a preconditioner trained on ensembles with small lattice volumes can be used to construct preconditioners for ensembles with many times larger lattice volumes, with minimal degradation of performance. This volume-transferring technique amortizes the training cost and presents a pathway towards scaling such preconditioners to lattice field theory calculations with larger lattice volumes and in four dimensions. △ Less

Submitted 4 August, 2022; originally announced August 2022.

Comments: 11 pages, 6 figures, and 2 tables

The Profiled Feldman-Cousins technique for confidence interval construction in the presence of nuisance parameters

Authors: M. A. Acero, B. Acharya, P. Adamson, L. Aliaga, N. Anfimov, A. Antoshkin, E. Arrieta-Diaz, L. Asquith, A. Aurisano, A. Back, C. Backhouse, M. Baird, N. Balashov, P. Baldi, B. A. Bambah, S. Bashar, A. Bat, K. Bays, R. Bernstein, V. Bhatnagar, D. Bhattarai, B. Bhuyan, J. Bian, A. C. Booth, R. Bowles , et al. (196 additional authors not shown)

Abstract: Measuring observables to constrain models using maximum-likelihood estimation is fundamental to many physics experiments. The Profiled Feldman-Cousins method described here is a potential solution to common challenges faced in constructing accurate confidence intervals: small datasets, bounded parameters, and the need to properly handle nuisance parameters. This method achieves more accurate frequ… ▽ More Measuring observables to constrain models using maximum-likelihood estimation is fundamental to many physics experiments. The Profiled Feldman-Cousins method described here is a potential solution to common challenges faced in constructing accurate confidence intervals: small datasets, bounded parameters, and the need to properly handle nuisance parameters. This method achieves more accurate frequentist coverage than other methods in use, and is generally applicable to the problem of parameter estimation in neutrino oscillations and similar measurements. We describe an implementation of this method in the context of the NOvA experiment. △ Less

Submitted 1 August, 2022; v1 submitted 28 July, 2022; originally announced July 2022.

Comments: 19 pages, 12 figures

Report number: FERMILAB-PUB-22-476-ND

Gauge-equivariant flow models for sampling in lattice field theories with pseudofermions

Authors: Ryan Abbott, Michael S. Albergo, Denis Boyda, Kyle Cranmer, Daniel C. Hackett, Gurtej Kanwar, Sébastien Racanière, Danilo J. Rezende, Fernando Romero-López, Phiala E. Shanahan, Betsy Tian, Julian M. Urban

Abstract: This work presents gauge-equivariant architectures for flow-based sampling in fermionic lattice field theories using pseudofermions as stochastic estimators for the fermionic determinant. This is the default approach in state-of-the-art lattice field theory calculations, making this development critical to the practical application of flow models to theories such as QCD. Methods by which flow-base… ▽ More This work presents gauge-equivariant architectures for flow-based sampling in fermionic lattice field theories using pseudofermions as stochastic estimators for the fermionic determinant. This is the default approach in state-of-the-art lattice field theory calculations, making this development critical to the practical application of flow models to theories such as QCD. Methods by which flow-based sampling approaches can be improved via standard techniques such as even/odd preconditioning and the Hasenbusch factorization are also outlined. Numerical demonstrations in two-dimensional U(1) and SU(3) gauge theories with $N_f=2$ flavors of fermions are provided. △ Less

Submitted 16 October, 2022; v1 submitted 18 July, 2022; originally announced July 2022.

Comments: 15 pages, 7 figures. v3: accepted version for publication. New appendix C

Report number: MIT-CTP/5446, INT-PUB-22-017

Journal ref: Phys.Rev.D 106 (2022) 7, 074506

Lattice QCD and Particle Physics

Authors: Andreas S. Kronfeld, Tanmoy Bhattacharya, Thomas Blum, Norman H. Christ, Carleton DeTar, William Detmold, Robert Edwards, Anna Hasenfratz, Huey-Wen Lin, Swagato Mukherjee, Konstantinos Orginos, Richard Brower, Vincenzo Cirigliano, Zohreh Davoudi, Bálint Jóo, Chulwoo Jung, Christoph Lehner, Stefan Meinel, Ethan T. Neil, Peter Petreczky, David G. Richards, Alexei Bazavov, Simon Catterall, Jozef J. Dudek, Aida X. El-Khadra , et al. (57 additional authors not shown)

Abstract: Contribution from the USQCD Collaboration to the Proceedings of the US Community Study on the Future of Particle Physics (Snowmass 2021). Contribution from the USQCD Collaboration to the Proceedings of the US Community Study on the Future of Particle Physics (Snowmass 2021). △ Less

Submitted 2 October, 2022; v1 submitted 15 July, 2022; originally announced July 2022.

Comments: 27 pp. main text, 4 pp. appendices, 29 pp. references, 1 p. index

Report number: FERMILAB-CONF-22-531-T

Towards Precise and Accurate Calculations of Neutrinoless Double-Beta Decay: Project Scoping Workshop Report

Authors: V. Cirigliano, Z. Davoudi, J. Engel, R. J. Furnstahl, G. Hagen, U. Heinz, H. Hergert, M. Horoi, C. W. Johnson, A. Lovato, E. Mereghetti, W. Nazarewicz, A. Nicholson, T. Papenbrock, S. Pastore, M. Plumlee, D. R. Phillips, P. E. Shanahan, S. R. Stroberg, F. Viens, A. Walker-Loud, K. A. Wendt, S. M. Wild

Abstract: We present the results of a National Science Foundation (NSF) Project Scoping Workshop, the purpose of which was to assess the current status of calculations for the nuclear matrix elements governing neutrinoless double-beta decay and determine if more work on them is required. After reviewing important recent progress in the application of effective field theory, lattice quantum chromodynamics, a… ▽ More We present the results of a National Science Foundation (NSF) Project Scoping Workshop, the purpose of which was to assess the current status of calculations for the nuclear matrix elements governing neutrinoless double-beta decay and determine if more work on them is required. After reviewing important recent progress in the application of effective field theory, lattice quantum chromodynamics, and ab initio nuclear-structure theory to double-beta decay, we discuss the state of the art in nuclear-physics uncertainty quantification and then construct a road map for work in all these areas to fully complement the increasingly sensitive experiments in operation and under development. The road map contains specific projects in theoretical and computational physics as well as an uncertainty-quantification plan that employs Bayesian Model Mixing and an analysis of correlations between double-beta-decay rates and other observables. The goal of this program is a set of accurate and precise matrix elements, in all nuclei of interest to experimentalists, delivered together with carefully assessed uncertainties. Such calculations will allow crisp conclusions from the observation or non-observation of neutrinoless double-beta decay, no matter what new physics is at play. △ Less

Submitted 3 July, 2022; originally announced July 2022.

Comments: This Project Scoping Workshop report is focused on the US context for the theory of neutrinloess double beta decay. Its authors plan to produce a journal article that addresses similar issues, but is more inclusive as regards non-US efforts on this problem. We would be happy to receive further input that will help us refine our text before it is submitted to the journal

Report number: INT-PUB-22-018

Journal ref: J. Phys. G: Nucl. Part. Phys. 49, 120502 (2022)

Measurement of the $ν_e-$Nucleus Charged-Current Double-Differential Cross Section at $\left< E_ν \right> = $ 2.4 GeV using NOvA

Authors: M. A. Acero, P. Adamson, L. Aliaga, N. Anfimov, A. Antoshkin, E. Arrieta-Diaz, L. Asquith, A. Aurisano, A. Back, C. Backhouse, M. Baird, N. Balashov, P. Baldi, B. A. Bambah, S. Bashar, K. Bays, R. Bernstein, V. Bhatnagar, D. Bhattarai, B. Bhuyan, J. Bian, A. C. Booth, R. Bowles, B. Brahma, C. Bromberg , et al. (190 additional authors not shown)

Abstract: The inclusive electron neutrino charged-current cross section is measured in the NOvA near detector using $8.02\times10^{20}$ protons-on-target (POT) in the NuMI beam. The sample of GeV electron neutrino interactions is the largest analyzed to date and is limited by $\simeq$ 17\% systematic rather than the $\simeq$ 7.4\% statistical uncertainties. The double-differential cross section in final-sta… ▽ More The inclusive electron neutrino charged-current cross section is measured in the NOvA near detector using $8.02\times10^{20}$ protons-on-target (POT) in the NuMI beam. The sample of GeV electron neutrino interactions is the largest analyzed to date and is limited by $\simeq$ 17\% systematic rather than the $\simeq$ 7.4\% statistical uncertainties. The double-differential cross section in final-state electron energy and angle is presented for the first time, together with the single-differential dependence on $Q^{2}$ (squared four-momentum transfer) and energy, in the range 1 GeV $ \leq E_ν < $6 GeV. Detailed comparisons are made to the predictions of the GENIE, GiBUU, NEUT, and NuWro neutrino event generators. The data do not strongly favor a model over the others consistently across all three cross sections measured, though some models have especially good or poor agreement in the single differential cross section vs. $Q^{2}$. △ Less

Submitted 21 June, 2022; originally announced June 2022.

Report number: FERMILAB-PUB-22-446-ND

Theoretical tools for neutrino scattering: interplay between lattice QCD, EFTs, nuclear physics, phenomenology, and neutrino event generators

Authors: L. Alvarez Ruso, A. M. Ankowski, S. Bacca, A. B. Balantekin, J. Carlson, S. Gardiner, R. Gonzalez-Jimenez, R. Gupta, T. J. Hobbs, M. Hoferichter, J. Isaacson, N. Jachowicz, W. I. Jay, T. Katori, F. Kling, A. S. Kronfeld, S. W. Li, H. -W. Lin, K. -F. Liu, A. Lovato, K. Mahn, J. Menendez, A. S. Meyer, J. Morfin, S. Pastore , et al. (36 additional authors not shown)

Abstract: Maximizing the discovery potential of increasingly precise neutrino experiments will require an improved theoretical understanding of neutrino-nucleus cross sections over a wide range of energies. Low-energy interactions are needed to reconstruct the energies of astrophysical neutrinos from supernovae bursts and search for new physics using increasingly precise measurement of coherent elastic neut… ▽ More Maximizing the discovery potential of increasingly precise neutrino experiments will require an improved theoretical understanding of neutrino-nucleus cross sections over a wide range of energies. Low-energy interactions are needed to reconstruct the energies of astrophysical neutrinos from supernovae bursts and search for new physics using increasingly precise measurement of coherent elastic neutrino scattering. Higher-energy interactions involve a variety of reaction mechanisms including quasi-elastic scattering, resonance production, and deep inelastic scattering that must all be included to reliably predict cross sections for energies relevant to DUNE and other accelerator neutrino experiments. This white paper discusses the theoretical status, challenges, required resources, and path forward for achieving precise predictions of neutrino-nucleus scattering and emphasizes the need for a coordinated theoretical effort involved lattice QCD, nuclear effective theories, phenomenological models of the transition region, and event generators. △ Less

Submitted 20 April, 2022; v1 submitted 16 March, 2022; originally announced March 2022.

Comments: 81 pages, contribution to Snowmass 2021

Report number: DESY-22-05, FERMILAB-FN-1161-T, MITP-22-027

Symmetry Group Equivariant Architectures for Physics

Authors: Alexander Bogatskiy, Sanmay Ganguly, Thomas Kipf, Risi Kondor, David W. Miller, Daniel Murnane, Jan T. Offermann, Mariel Pettee, Phiala Shanahan, Chase Shimmin, Savannah Thais

Abstract: Physical theories grounded in mathematical symmetries are an essential component of our understanding of a wide range of properties of the universe. Similarly, in the domain of machine learning, an awareness of symmetries such as rotation or permutation invariance has driven impressive performance breakthroughs in computer vision, natural language processing, and other important applications. In t… ▽ More Physical theories grounded in mathematical symmetries are an essential component of our understanding of a wide range of properties of the universe. Similarly, in the domain of machine learning, an awareness of symmetries such as rotation or permutation invariance has driven impressive performance breakthroughs in computer vision, natural language processing, and other important applications. In this report, we argue that both the physics community and the broader machine learning community have much to understand and potentially to gain from a deeper investment in research concerning symmetry group equivariant machine learning architectures. For some applications, the introduction of symmetries into the fundamental structural design can yield models that are more economical (i.e. contain fewer, but more expressive, learned parameters), interpretable (i.e. more explainable or directly mappable to physical quantities), and/or trainable (i.e. more efficient in both data and computational requirements). We discuss various figures of merit for evaluating these models as well as some potential benefits and limitations of these methods for a variety of physics applications. Research and investment into these approaches will lay the foundation for future architectures that are potentially more robust under new computational paradigms and will provide a richer description of the physical systems to which they are applied. △ Less

Submitted 11 March, 2022; originally announced March 2022.

Comments: Contribution to Snowmass 2021

Low-Energy Physics in Neutrino LArTPCs

Authors: D. Caratelli, W. Foreman, A. Friedland, S. Gardiner, I. Gil-Botella, G. Karagiorgi, M. Kirby, G. Lehmann Miotto, B. R. Littlejohn, M. Mooney, J. Reichenbacher, A. Sousa, K. Scholberg, J. Yu, T. Yang, S. Andringa, J. Asaadi, T. J. C. Bezerra, F. Capozzi, F. Cavanna, E. Church, A. Himmel, T. Junk, J. Klein, I. Lepetic , et al. (264 additional authors not shown)

Abstract: In this white paper, we outline some of the scientific opportunities and challenges related to detection and reconstruction of low-energy (less than 100 MeV) signatures in liquid argon time-projection chamber (LArTPC) detectors. Key takeaways are summarized as follows. 1) LArTPCs have unique sensitivity to a range of physics and astrophysics signatures via detection of event features at and below… ▽ More In this white paper, we outline some of the scientific opportunities and challenges related to detection and reconstruction of low-energy (less than 100 MeV) signatures in liquid argon time-projection chamber (LArTPC) detectors. Key takeaways are summarized as follows. 1) LArTPCs have unique sensitivity to a range of physics and astrophysics signatures via detection of event features at and below the few tens of MeV range. 2) Low-energy signatures are an integral part of GeV-scale accelerator neutrino interaction final states, and their reconstruction can enhance the oscillation physics sensitivities of LArTPC experiments. 3) BSM signals from accelerator and natural sources also generate diverse signatures in the low-energy range, and reconstruction of these signatures can increase the breadth of BSM scenarios accessible in LArTPC-based searches. 4) Neutrino interaction cross sections and other nuclear physics processes in argon relevant to sub-hundred-MeV LArTPC signatures are poorly understood. Improved theory and experimental measurements are needed. Pion decay-at-rest sources and charged particle and neutron test beams are ideal facilities for experimentally improving this understanding. 5) There are specific calibration needs in the low-energy range, as well as specific needs for control and understanding of radiological and cosmogenic backgrounds. 6) Novel ideas for future LArTPC technology that enhance low-energy capabilities should be explored. These include novel charge enhancement and readout systems, enhanced photon detection, low radioactivity argon, and xenon doping. 7) Low-energy signatures, whether steady-state or part of a supernova burst or larger GeV-scale event topology, have specific triggering, DAQ and reconstruction requirements that must be addressed outside the scope of conventional GeV-scale data collection and analysis pathways. △ Less

Submitted 1 March, 2022; originally announced March 2022.

Comments: Contribution to Snowmass 2021

Flow-based sampling in the lattice Schwinger model at criticality

Authors: Michael S. Albergo, Denis Boyda, Kyle Cranmer, Daniel C. Hackett, Gurtej Kanwar, Sébastien Racanière, Danilo J. Rezende, Fernando Romero-López, Phiala E. Shanahan, Julian M. Urban

Abstract: Recent results suggest that flow-based algorithms may provide efficient sampling of field distributions for lattice field theory applications, such as studies of quantum chromodynamics and the Schwinger model. In this work, we provide a numerical demonstration of robust flow-based sampling in the Schwinger model at the critical value of the fermion mass. In contrast, at the same parameters, conven… ▽ More Recent results suggest that flow-based algorithms may provide efficient sampling of field distributions for lattice field theory applications, such as studies of quantum chromodynamics and the Schwinger model. In this work, we provide a numerical demonstration of robust flow-based sampling in the Schwinger model at the critical value of the fermion mass. In contrast, at the same parameters, conventional methods fail to sample all parts of configuration space, leading to severely underestimated uncertainties. △ Less

Submitted 23 February, 2022; originally announced February 2022.

Comments: 5 pages main text, 3 pages supplementary material. 4 figures

Report number: MIT-CTP/5409

Lattice QCD Calculations of Parton Physics

Authors: Martha Constantinou, Luigi Del Debbio, Xiangdong Ji, Huey-Wen Lin, Keh-Fei Liu, Christopher Monahan, Kostas Orginos, Peter Petreczky, Jian-Wei Qiu, David Richards, Nobuo Sato, Phiala Shanahan, C. -P. Yuan, Jian-Hui Zhang, Yong Zhao

Abstract: In this document, we summarize the status and challenges of calculating parton physics in lattice QCD for the US Particle Physics Community Planning Exercise (a.k.a. "Snowmass"). While PDF-moments calculations have been very successful and been continuously improved, new methods have been developed to calculate distributions directly in $x$-space. Many recent lattice studies have been focused on c… ▽ More In this document, we summarize the status and challenges of calculating parton physics in lattice QCD for the US Particle Physics Community Planning Exercise (a.k.a. "Snowmass"). While PDF-moments calculations have been very successful and been continuously improved, new methods have been developed to calculate distributions directly in $x$-space. Many recent lattice studies have been focused on calculating isovector PDFs of the pion and nucleon, learning to control systematics associated with excited-state contamination, renormalization and continuum extrapolations, pion-mass and finite-volume effects, etc. Although in some cases, the lattice results are already competitive with experimental data, to reach the level of precision in a wide range of $x$ for unpolarized nucleon PDFs impactful for future collider physics remains a challenge, and may require exascale supercomputing power. The new theoretical methods open the door for calculating other partonic observables which will be the focus of the experimental program in nuclear physics, including generalized parton distributions and transverse-momentum dependent PDFs. A fruitful interplay between experimental data and lattice-QCD calculations will usher in a new era for parton physics and hadron structure. △ Less

Submitted 15 February, 2022; originally announced February 2022.

Comments: 14 pages

Report number: JLAB-THY-22-3564,MIT-CTP/5408,MSUHEP-22-004

Applications of Machine Learning to Lattice Quantum Field Theory

Authors: Denis Boyda, Salvatore Calì, Sam Foreman, Lena Funcke, Daniel C. Hackett, Yin Lin, Gert Aarts, Andrei Alexandru, Xiao-Yong Jin, Biagio Lucini, Phiala E. Shanahan

Abstract: There is great potential to apply machine learning in the area of numerical lattice quantum field theory, but full exploitation of that potential will require new strategies. In this white paper for the Snowmass community planning process, we discuss the unique requirements of machine learning for lattice quantum field theory research and outline what is needed to enable exploration and deployment… ▽ More There is great potential to apply machine learning in the area of numerical lattice quantum field theory, but full exploitation of that potential will require new strategies. In this white paper for the Snowmass community planning process, we discuss the unique requirements of machine learning for lattice quantum field theory research and outline what is needed to enable exploration and deployment of this approach in the future. △ Less

Submitted 10 February, 2022; originally announced February 2022.

Comments: 10 pages, contribution to Snowmass 2022

Report number: MIT-CTP/5405