Search

We present non-perturbative lattice calculations of the low-lying meson and baryon spectrum of the SU(4) gauge theory with fundamental fermion constituents. This theory is one instance of stealth dark matter, a class of strongly coupled theories, where the lowest mass stable baryon is the dark matter candidate. This work constitutes the first milestone in the program to study stealth dark matter self-interactions. Here, we focus on reducing excited state contamination in the single baryon channel by applying the Laplacian Heaviside method, as well as projecting our baryon operators onto the irreducible representations of the octahedral group. We compare our resulting spectrum to previous work involving Gaussian smeared non-projected operators and find good agreement with reduced statistical uncertainties. We also present the spectrum of the low-lying odd-parity baryons for the first time., Comment: 19 pages, 9 figures

Reducing memory traffic is critical to accelerate Lattice QCD computations on modern processors, given that such computations are memory-bandwidth bound. A commonly used strategy is mixed-precision solvers, however, these require careful treatment to ensure stable convergence. We give an overview of the strategies employed in QUDA to stabilize mixed-precision variants of Conjugate Gradient (CG), and its multi-shift brethren. Through the use of customized numerical storage formats we can significantly improve upon the precision achievable compared to IEEE numerical formats, increasing both the solver precision and stability achievable at fixed word size. We give examples using BiCGStab(l) and multi-shift CG solvers using the HISQ operator., Comment: 14 pages, 4 figures

Adaptive multi-grid methods have proven very successful in dealing with critical slow down for the Wilson-Dirac solver in lattice gauge theory. Multi-grid algorithms developed for Staggered fermions using the K\"ahler-Dirac preconditioning~\cite{Brower:2018ymy} have shown remarkable success. In this work, we discuss the performance of this staggered multi-grid algorithm in four dimensions. We also demonstrate that offloading some components of a multi-shift solve to a multi-grid solver leads to a significant performance improvement in an existing MILC spectrum workflow on the Summit and Selene supercomputers., Comment: Submission to Proceedings of Lattice 2022: the 39th International Symposium on Lattice Field Theory, Bonn, Germany

Contribution from the USQCD Collaboration to the Proceedings of the US Community Study on the Future of Particle Physics (Snowmass 2021)., Comment: 27 pp. main text, 4 pp. appendices, 29 pp. references, 1 p. index

The exascale race is at an end with the announcement of the Aurora and Frontier machines. This next generation of supercomputers utilize diverse hardware architectures to achieve their compute performance, providing an added onus on the performance portability of applications. An expanding fragmentation of programming models would provide a compounding optimization challenge were it not for the evolution of performance-portable frameworks, providing unified models for mapping abstract hierarchies of parallelism to diverse architectures. A solution to this challenge is the evolution of performance-portable frameworks, providing unified models for mapping abstract hierarchies of parallelism to diverse architectures. Kokkos is one such performance portable programming model for C++ applications, providing back-end implementations for each major HPC platform. Even with a performance portable framework, restructuring algorithms to expose higher degrees of parallelism is non-trivial. The Spectral Neighbor Analysis Potential (SNAP) is a machine-learned inter-atomic potential utilized in cutting-edge molecular dynamics simulations. Previous implementations of the SNAP calculation showed a downward trend in their performance relative to peak on newer-generation CPUs and low performance on GPUs. In this paper we describe the restructuring and optimization of SNAP as implemented in the Kokkos CUDA backend of the LAMMPS molecular dynamics package, benchmarked on NVIDIA GPUs. We identify novel patterns of hierarchical parallelism, facilitating a minimization of memory access overheads and pushing the implementation into a compute-saturated regime. Our implementation via Kokkos enables recompile-and-run efficiency on upcoming architectures. We find a $\sim$22x time-to-solution improvement relative to an existing implementation as measured on an NVIDIA Tesla V100-16GB for an important benchmark., Comment: Submitted to IPDPS 2021, October 19, 2020

The quantum extension of classical finite elements, referred to as quantum finite elements ({\bf QFE})~\cite{Brower:2018szu,Brower:2016vsl}, is applied to the radial quantization of 3d $\phi^4$ theory on a simplicial lattice for the $\mathbb R \times \mathbb S^2$ manifold. Explicit counter terms to cancel the one- and two-loop ultraviolet defects are implemented to reach the quantum continuum theory. Using the Brower-Tamayo~\cite{Brower:1989mt} cluster Monte Carlo algorithm, numerical results support the QFE ansatz that the critical conformal field theory (CFT) is reached in the continuum with the full isometries of $\mathbb R \times \mathbb S^2$ restored. The Ricci curvature term, while technically irrelevant in the quantum theory, is shown to dramatically improve the convergence opening, the way for high precision Monte Carlo simulation to determine the CFT data: operator dimensions, trilinear OPE couplings and the central charge., Comment: 8 pages, 7 figures

Critical slowing down for the Krylov Dirac solver presents a major obstacle to further advances in lattice field theory as it approaches the continuum solution. We propose a new multi-grid approach for chiral fermions, applicable to both the 5-d domain wall or 4-d Overlap operator. The central idea is to directly coarsen the 4-d Wilson kernel, giving an effective domain wall or overlap operator on each level. We provide here an explicit construction for the Shamir domain wall formulation with numerical tests for the 2-d Schwinger prototype, demonstrating near ideal multi-grid scaling. The framework is designed for a natural extension to 4-d lattice QCD chiral fermions, such as the M\"obius, Zolotarev or Borici domain wall discretizations or directly to a rational expansion of the 4-d Overlap operator. For the Shamir operator, the effective overlap operator is isolated by the use of a Pauli-Villars preconditioner in the spirit of the K\"ahler-Dirac spectral map used in a recent staggered MG algorithm [1]., Comment: 39 pages, 13 figures

This document is one of a series of whitepapers from the USQCD collaboration. Here, we discuss opportunities for lattice field theory research to make an impact on models of new physics beyond the Standard Model, including composite Higgs, composite dark matter, and supersymmetric theories., Comment: 24 pages

We present our lattice studies of SU(3) gauge theory with $N_f$ = 8 degenerate fermions in the fundamental representation. Using nHYP-smeared staggered fermions we study finite-temperature transitions on lattice volumes as large as $L^3 \times N_t = 48^3 \times 24$, and the zero-temperature composite spectrum on lattice volumes up to $64^3 \times 128$. The spectrum indirectly indicates spontaneous chiral symmetry breaking, but finite-temperature transitions with fixed $N_t \leq 24$ enter a strongly coupled lattice phase as the fermion mass decreases, which prevents a direct confirmation of spontaneous chiral symmetry breaking in the chiral limit. In addition to the connected spectrum we focus on the lightest flavor-singlet scalar particle. We find it to be degenerate with the pseudo-Goldstone states down to the lightest masses reached so far by non-perturbative lattice calculations. Using the same lattice approach, we study the behavior of the composite spectrum when the number of light fermions is changed from eight to four. A heavy flavor-singlet scalar in the 4-flavor theory affirms the contrast between QCD-like dynamics and the low-energy behavior of the 8-flavor theory., Comment: 31 pages, 36 figures, 8 tables. v2: update to published version

We present a method for defining a lattice realization of the $\phi^4$ quantum field theory on a simplicial complex in order to enable numerical computation on a general Riemann manifold. The procedure begins with adopting methods from traditional Regge Calculus (RC) and finite element methods (FEM) plus the addition of ultraviolet counter terms required to reach the renormalized field theory in the continuum limit. The construction is tested numerically for the two-dimensional $\phi^4$ scalar field theory on the Riemann two-sphere, $\mathbb{S}^2$, in comparison with the exact solutions to the two-dimensional Ising conformal field theory (CFT). Numerical results for the Binder cumulants (up to 12th order) and the two- and four-point correlation functions are in agreement with the exact $c = 1/2$ CFT solutions., Comment: 52 pages, 27 figures

Critical slowing down in Krylov methods for the Dirac operator presents a major obstacle to further advances in lattice field theory as it approaches the continuum solution. Here we formulate a multi-grid algorithm for the Kogut-Susskind (or staggered) fermion discretization which has proven difficult relative to Wilson multigrid due to its first-order anti-Hermitian structure. The solution is to introduce a novel spectral transformation by the K\"ahler-Dirac spin structure prior to the Galerkin projection. We present numerical results for the two-dimensional, two-flavor Schwinger model, however, the general formalism is agnostic to dimension and is directly applicable to four-dimensional lattice QCD., Comment: 48 pages, 37 figures

Lattice quantum chromodynamics simulations in nuclear physics have benefited from a tremendous number of algorithmic advances such as multigrid and eigenvector deflation. These improve the time to solution but do not alleviate the intrinsic memory-bandwidth constraints of the matrix-vector operation dominating iterative solvers. Batching this operation for multiple vectors and exploiting cache and register blocking can yield a super-linear speed up. Block-Krylov solvers can naturally take advantage of such batched matrix-vector operations, further reducing the iterations to solution by sharing the Krylov space between solves. However, practical implementations typically suffer from the quadratic scaling in the number of vector-vector operations. Using the QUDA library, we present an implementation of a block-CG solver on NVIDIA GPUs which reduces the memory-bandwidth complexity of vector-vector operations from quadratic to linear. We present results for the HISQ discretization, showing a 5x speedup compared to highly-optimized independent Krylov solves on NVIDIA's SaturnV cluster., Comment: 15 pages, 14 figures, in press

The lattice Dirac equation is formulated on a simplicial complex which approximates a smooth Riemann manifold by introducing a lattice vierbein on each site and a lattice spin connection on each link. Care is taken so the construction applies to any smooth D-dimensional Riemannian manifold that permits a spin connection. It is tested numerically in 2D for the projective sphere ${\mathbb S}^2$ in the limit of an increasingly refined sequence of triangles. The eigenspectrum and eigenvectors are shown to converge rapidly to the exact result in the continuum limit. In addition comparison is made with the continuum Ising conformal field theory on ${\mathbb S}^2$. Convergence is tested for the two point, $\langle \epsilon(x_1) \epsilon(x_2) \rangle$, and the four point, $\langle \sigma(x_1) \epsilon(x_2) \epsilon(x_3 )\sigma(x_4) \rangle $, correlators for the energy, $\epsilon(x) = i \bar \psi(x)\psi(x)$, and twist operators, $\sigma(x)$, respectively., Comment: 53 pages, 29 figures

We present results for the spectrum of a strongly interacting SU(3) gauge theory with $N_f = 8$ light fermions in the fundamental representation. Carrying out non-perturbative lattice calculations at the lightest masses and largest volumes considered to date, we confirm the existence of a remarkably light singlet scalar particle. We explore the rich resonance spectrum of the 8-flavor theory in the context of the search for new physics beyond the standard model at the Large Hadron Collider (LHC). Connecting our results to models of dynamical electroweak symmetry breaking, we estimate the vector resonance mass to be about 2 TeV with a width of roughly 450 GeV, and predict additional resonances with masses below ~3 TeV., Comment: 6 pages, 6 figures. Added report number. Version submitted to journal

Viable non-perturbative methods for lattice quantum field theories on curved manifolds are difficult. By adapting features from the traditional finite element methods (FEM) and Regge Calculus, a new simplicial lattice Quantum Finite Element (QFE) Lagrangian is constructed for fields on a smooth Riemann manifold. To reach the continuum limit additional counter terms must be constructed to cancel the ultraviolet distortions. This is tested by the comparison of phi 4-th theory at the Wilson-Fisher fixed point with the exact Ising (c =1/2) CFT on a 2D Riemann sphere. The Dirac equation is also constructed on a simplicial lattice approximation to a Riemann manifold by introducing a lattice vierbein and spin connection on each link. Convergence of the QFE Dirac equation is tested against the exact solution for the 2D Riemann sphere. Future directions and applications to Conformal Field Theories are suggested., Comment: 21 pages, 17 figures, The 33rd International Symposium on Lattice Field Theory

We calculate the spin-independent scattering cross section for direct detection that results from the electromagnetic polarizability of a composite scalar baryon dark matter candidate -- "Stealth Dark Matter", that is based on a dark SU(4) confining gauge theory. In the nonrelativistic limit, electromagnetic polarizability proceeds through a dimension-7 interaction leading to a very small scattering cross section for dark matter with weak scale masses. This represents a lower bound on the scattering cross section for composite dark matter theories with electromagnetically charged constituents. We carry out lattice calculations of the polarizability for the lightest baryons in SU(3) and SU(4) gauge theories using the background field method on quenched configurations. We find the polarizabilities of SU(3) and SU(4) to be comparable (within about 50%) normalized to the baryon mass, which is suggestive for extensions to larger SU(N) groups. The resulting scattering cross sections with a xenon target are shown to be potentially detectable in the dark matter mass range of about 200-700 GeV, where the lower bound is from the existing LUX constraint while the upper bound is the coherent neutrino background. Significant uncertainties in the cross section remain due to the more complicated interaction of the polarizablity operator with nuclear structure, however the steep dependence on the dark matter mass, $1/m_B^6$, suggests the observable dark matter mass range is not appreciably modified. We briefly highlight collider searches for the mesons in the theory as well as the indirect astrophysical effects that may also provide excellent probes of stealth dark matter., Comment: 6 pages, 2 figures, citations added, typos fixed, minor clarifications

We present a new model of "Stealth Dark Matter": a composite baryonic scalar of an $SU(N_D)$ strongly-coupled theory with even $N_D \geq 4$. All mass scales are technically natural, and dark matter stability is automatic without imposing an additional discrete or global symmetry. Constituent fermions transform in vector-like representations of the electroweak group that permit both electroweak-breaking and electroweak-preserving mass terms. This gives a tunable coupling of stealth dark matter to the Higgs boson independent of the dark matter mass itself. We specialize to $SU(4)$, and investigate the constraints on the model from dark meson decay, electroweak precision measurements, basic collider limits, and spin-independent direct detection scattering through Higgs exchange. We exploit our earlier lattice simulations that determined the composite spectrum as well as the effective Higgs coupling of stealth dark matter in order to place bounds from direct detection, excluding constituent fermions with dominantly electroweak-breaking masses. A lower bound on the dark baryon mass $m_B \gtrsim 300$ GeV is obtained from the indirect requirement that the lightest dark meson not be observable at LEP II. We briefly survey some intriguing properties of stealth dark matter that are worthy of future study, including: collider studies of dark meson production and decay; indirect detection signals from annihilation; relic abundance estimates for both symmetric and asymmetric mechanisms; and direct detection through electromagnetic polarizability, a detailed study of which will appear in a companion paper., Comment: 15 pages, 3 figures, citations added, typos fixed, minor clarifications

The light Higgs boson of the Standard Model could arise as the consequence of the weakly broken conformal symmetry in a strongly interacting gauge theory. Here we present a novel idea to study the transition from conformal to confining behavior using an SU(3) gauge theory with four light and eight heavy flavors. This system interpolates between the 12-flavor conformal and the 4 flavor chirally broken theory as the mass of the heavy flavors are varied. We show first results on our determination of the iso-singlet 0++ state., Comment: 6 pages, 4 figures, contributed talk to the XXVI IUPAP Conference on Computational Physics (CCP2014), August 11-14 2014, Boston, Massachusetts, USA

We study the transition between spontaneous chiral symmetry breaking and conformal behavior in the SU(3) theory with multiple fermion flavors. Instead of the traditional approach of changing the number of flavors, we keep the number of fermions fixed but lift the mass of a subset, keeping the remaining fermions near to the massless chiral limit. This way we can interpolate continuously between the conformal and chirally broken dynamics. In particular, we consider four light and eight heavy flavors and investigate the running/walking gauge coupling and the low energy meson spectrum, including the 0++ iso-singlet scalar state in this system. Our preliminary data reveal an iso-singlet scalar that is considerably lighter than the pion at large fermion mass but becomes heavier at smaller masses. This behavior is of particular phenomenological interest., Comment: 14 pages, 8 figures, combined contributions to the 32nd International Symposium on Lattice Field Theory (Lattice 2014), June 23-28 2014, New York, USA

We investigate the transition between spontaneous chiral symmetry breaking and conformal behavior in the SU(3) theory with multiple fermion flavors. We propose a new strategy for studying this transition. Instead of changing the number of flavors, we lift the mass of a subset of the fermions, keeping the rest of the fermions near to the massless chiral limit in order to probe the transition. Dedicated to the 60th birthday of Academician Valery Rubakov., Comment: 11 pages, 5 figures, fixed typos, added references

We present the spectrum of baryons in a new SU(4) gauge theory with fundamental fermion constituents. The spectrum of these bosonic baryons is of significant interest for composite dark matter theories. Here, we compare the spectrum and properties of SU(3) and SU(4) baryons, and then compute the dark-matter direct detection cross section via Higgs boson exchange for TeV-scale composite dark matter arising from a confining SU(4) gauge sector. Comparison with the latest LUX results leads to tight bounds on the fraction of the constituent-fermion mass that may arise from electroweak symmetry breaking. Lattice calculations of the dark matter mass spectrum and the Higgs-dark matter coupling are performed on quenched $16^{3} \times 32$, $32^{3} \times 64$, $48^{3} \times 96$, and $64^{3} \times128$ lattices with three different lattice spacings, using Wilson fermions with moderate to heavy pseudoscalar meson masses. Our results lay a foundation for future analytic and numerical study of composite baryonic dark matter., Comment: 18 pages, 18 figures

We compare observables to the finite size scaling hypothesis in SU(2) lattice gauge theory with two Dirac fermions in the adjoint representation. The fits that we obtain yield an estimate of the anomalous mass dimension that is consistent with four loop perturbation theory: $\gamma = 0.51 \pm 0.16$, with the error due to systematic uncertainties in the finite size scaling analysis. The result is somewhat larger than one Schr\"odinger functional study (by 1.3$\sigma$) but consistent with another., Comment: 5 pages, 4 figures

We extract the running coupling based on Creutz ratios in SU(2) lattice gauge theory with two Dirac fermions in the adjoint representation. Depending on how the extrapolation to zero fermion mass is performed, either backward running or an absence of running is observed at strong bare coupling. This behavior is consistent with other findings which indicate that this theory has an infrared fixed point., Comment: 12 pages, 3 figures; v2 fifth mass point added, three extrapolations now performed, conclusions softened