Source author record

Evan Berkowitz

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

27works
15topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

27 published item(s)

preprint2022arXiv

Detailed analysis of excited state systematics in a lattice QCD calculation of $g_A$

Excited state contamination remains one of the most challenging sources of systematic uncertainty to control in lattice QCD calculations of nucleon matrix elements and form factors: early time separations are contaminated by excited states and late times suffer from an exponentially bad signal-to-noise problem. High-statistics calculations at large time separations $\gtrsim1$ fm are commonly used to combat these issues. In this work, focusing on $g_A$, we explore the alternative strategy of utilizing a large number of relatively low-statistics calculations at short to medium time separations (0.2--1 fm), combined with a multi-state analysis. On an ensemble with a pion mass of approximately 310 MeV and a lattice spacing of approximately 0.09 fm, we find this provides a more robust and economical method of quantifying and controlling the excited state systematic uncertainty. A quantitative separation of various types of excited states enables the identification of the transition matrix elements as the dominant contamination. The excited state contamination of the Feynman-Hellmann correlation function is found to reduce to the 1% level at approximately 1 fm while for the more standard three-point functions, this does not occur until after 2 fm. Critical to our findings is the use of a global minimization, rather than fixing the spectrum from the two-point functions and using them as input to the three-point analysis. We find that the ground state parameters determined in such a global analysis are stable against variations in the excited state model, the number of excited states, and the truncation of early-time or late-time numerical data.

preprint2021arXiv

Two-nucleon S-wave interactions at the $SU(3)$ flavor-symmetric point with $m_{ud}\simeq m_s^{\rm phys}$: a first lattice QCD calculation with the stochastic Laplacian Heaviside method

We report on the first application of the stochastic Laplacian Heaviside method for computing multi-particle interactions with lattice QCD to the two-nucleon system. Like the Laplacian Heaviside method, this method allows for the construction of interpolating operators which can be used to construct a positive definite set of two-nucleon correlation functions, unlike nearly all other applications of lattice QCD to two nucleons in the literature. It also allows for a variational analysis in which optimal linear combinations of the interpolating operators are formed that couple predominantly to the eigenstates of the system. Utilizing such methods has become of paramount importance in order to help resolve the discrepancy in the literature on whether two nucleons in either isospin channel form a bound state at pion masses heavier than physical, with the discrepancy persisting even in the $SU(3)$-flavor symmetric point with all quark masses near the physical strange quark mass. This is the first in a series of papers aimed at resolving this discrepancy. In the present work, we employ the stochastic Laplacian Heaviside method without a hexaquark operator in the basis at a lattice spacing of $a\sim0.086$~fm, lattice volume of $L=48a\simeq4.1$~fm and pion mass $m_π\simeq714$ MeV. With this setup, the observed spectrum of two-nucleon energy levels strongly disfavors the presence of a bound state in either the deuteron or dineutron channel.

preprint2020arXiv

$F_K / F_π$ from Möbius domain-wall fermions solved on gradient-flowed HISQ ensembles

We report the results of a lattice quantum chromodynamics calculation of $F_K/F_π$ using Möbius domain-wall fermions computed on gradient-flowed $N_f=2+1+1$ highly-improved staggered quark (HISQ) ensembles. The calculation is performed with five values of the pion mass ranging from $130 \lesssim m_π\lesssim 400$ MeV, four lattice spacings of $a\sim 0.15, 0.12, 0.09$ and $0.06$ fm and multiple values of the lattice volume. The interpolation/extrapolation to the physical pion and kaon mass point, the continuum, and infinite volume limits are performed with a variety of different extrapolation functions utilizing both the relevant mixed-action effective field theory expressions as well as discretization-enhanced continuum chiral perturbation theory formulas. We find that the $a\sim0.06$ fm ensemble is helpful, but not necessary to achieve a subpercent determination of $F_K/F_π$. We also include an estimate of the strong isospin breaking corrections and arrive at a final result of $F_{K^\pm}/F_{π^\pm} = 1.1942(45)$ with all sources of statistical and systematic uncertainty included. This is consistent with the Flavour Lattice Averaging Group average value, providing an important benchmark for our lattice action. Combining our result with experimental measurements of the pion and kaon leptonic decays leads to a determination of $|V_{us}|/|V_{ud}| = 0.2311(10)$.

preprint2020arXiv

Avoiding Ergodicity Problems in Lattice Discretizations of the Hubbard Model

The Hubbard model arises naturally when electron-electron interactions are added to the tight-binding descriptions of many condensed matter systems. For instance, the two-dimensional Hubbard model on the honeycomb lattice is central to the ab initio description of the electronic structure of carbon nanomaterials, such as graphene. Such low-dimensional Hubbard models are advantageously studied with Markov chain Monte Carlo methods, such as Hybrid Monte Carlo (HMC). HMC is the standard algorithm of the lattice gauge theory community, as it is well suited to theories of dynamical fermions. As HMC performs continuous, global updates of the lattice degrees of freedom, it provides superior scaling with system size relative to local updating methods. A potential drawback of HMC is its susceptibility to ergodicity problems due to so-called exceptional configurations, for which the fermion operator cannot be inverted. Recently, ergodicity problems were found in some formulations of HMC simulations of the Hubbard model. Here, we address this issue directly and clarify under what conditions ergodicity is maintained or violated in HMC simulations of the Hubbard model. We study different lattice formulations of the fermion operator and provide explicit, representative calculations for small systems, often comparing to exact results. We show that a fermion operator can be found which is both computationally convenient and free of ergodicity problems.

preprint2020arXiv

Renormalization of a Contact Interaction on a Lattice

Contact interactions can be used to describe a system of particles at unitarity, contribute to the leading part of nuclear interactions and are numerically non-trivial because they require a proper regularization and renormalization scheme. We explain how to tune the coefficient of a contact interaction between non-relativistic particles on a discretized space in 1, 2, and 3 spatial dimensions such that we can remove all discretization artifacts. By taking advantage of a latticized Lüscher zeta function, we can achieve a momentum-independent scattering amplitude at any finite lattice spacing.

preprint2020arXiv

The Semimetal-Mott Insulator Quantum Phase Transition of the Hubbard Model on the Honeycomb Lattice

We take advantage of recent improvements in the grand canonical Hybrid Monte Carlo algorithm, to perform a precision study of the single-particle gap in the hexagonal Hubbard model, with on-site electron-electron interactions. After carefully controlled analyses of the Trotter error, the thermodynamic limit, and finite-size scaling with inverse temperature, we find a critical coupling of $U_c/κ=3.834(14)$ and the critical exponent $zν=1.185(43)$. Under the assumption that this corresponds to the expected anti-ferromagnetic Mott transition, we are also able to provide a preliminary estimate $β=1.095(37)$ for the critical exponent of the order parameter. We consider our findings in view of the $SU(2)$ Gross-Neveu, or chiral Heisenberg, universality class. We also discuss the computational scaling of the Hybrid Monte Carlo algorithm, and possible extensions of our work to carbon nanotubes, fullerenes, and topological insulators.

preprint2019arXiv

Lattice QCD Determination of $g_A$

The nucleon axial coupling, $g_A$, is a fundamental property of protons and neutrons, dictating the strength with which the weak axial current of the Standard Model couples to nucleons, and hence, the lifetime of a free neutron. The prominence of $g_A$ in nuclear physics has made it a benchmark quantity with which to calibrate lattice QCD calculations of nucleon structure and more complex calculations of electroweak matrix elements in one and few nucleon systems. There were a number of significant challenges in determining $g_A$, notably the notorious exponentially-bad signal-to-noise problem and the requirement for hundreds of thousands of stochastic samples, that rendered this goal more difficult to obtain than originally thought. I will describe the use of an unconventional computation method, coupled with "ludicrously'" fast GPU code, access to publicly available lattice QCD configurations from MILC and access to leadership computing that have allowed these challenges to be overcome resulting in a determination of $g_A$ with 1% precision and all sources of systematic uncertainty controlled. I will discuss the implications of these results for the convergence of $SU(2)$ Chiral Perturbation theory for nucleons, as well as prospects for further improvements to $g_A$ (sub-percent precision, for which we have preliminary results) which is part of a more comprehensive application of lattice QCD to nuclear physics. This is particularly exciting in light of the new CORAL supercomputers coming online, Sierra and Summit, for which our lattice QCD codes achieve a machine-to-machine speed up over Titan of an order of magnitude.

preprint2018arXiv

Simulating the weak death of the neutron in a femtoscale universe with near-Exascale computing

The fundamental particle theory called Quantum Chromodynamics (QCD) dictates everything about protons and neutrons, from their intrinsic properties to interactions that bind them into atomic nuclei. Quantities that cannot be fully resolved through experiment, such as the neutron lifetime (whose precise value is important for the existence of light-atomic elements that make the sun shine and life possible), may be understood through numerical solutions to QCD. We directly solve QCD using Lattice Gauge Theory and calculate nuclear observables such as neutron lifetime. We have developed an improved algorithm that exponentially decreases the time-to solution and applied it on the new CORAL supercomputers, Sierra and Summit. We use run-time autotuning to distribute GPU resources, achieving 20% performance at low node count. We also developed optimal application mapping through a job manager, which allows CPU and GPU jobs to be interleaved, yielding 15% of peak performance when deployed across large fractions of CORAL.

preprint2016arXiv

A microscopic description of black hole evaporation via holography

We propose a description of how a large, cold black hole (black zero-brane) in type IIA superstring theory evaporates into freely propagating D0-branes, by solving the dual gauge theory quantitatively. The energy spectrum of emitted D0-branes is parametrically close to thermal when the black hole is large. The black hole, while initially cold, gradually becomes an extremely hot and stringy object as it evaporates. As it emits D0-branes, its emission rate speeds up and it evaporates completely without leaving any remnant. Hence this system provides us with a concrete holographic description of black hole evaporation without information loss.

preprint2016arXiv

Neutrinoless double beta decay from lattice QCD

While the discovery of non-zero neutrino masses is one of the most important accomplishments by physicists in the past century, it is still unknown how and in what form these masses arise. Lepton number-violating neutrinoless double beta decay is a natural consequence of Majorana neutrinos and many BSM theories, and many experimental efforts are involved in the search for these processes. Understanding how neutrinoless double beta decay would manifest in nuclear environments is key for understanding any observed signals. In these proceedings we present an overview of a set of one- and two-body matrix elements relevant for experimental searches for neutrinoless double beta decay, describe the role of lattice QCD calculations, and present preliminary lattice QCD results.

preprint2016arXiv

Precision lattice test of the gauge/gravity duality at large-$N$

We pioneer a systematic, large-scale lattice simulation of D0-brane quantum mechanics. The large-$N$ and continuum limits of the gauge theory are taken for the first time at various temperatures $0.4 \leq T \leq 1.0$. As a way to directly test the gauge/gravity duality conjecture we compute the internal energy of the black hole directly from the gauge theory and reproduce the coefficient of the supergravity result $E/N^2=7.41T^{14/5}$. This is the first confirmation of the supergravity prediction for the internal energy of a black hole at finite temperature coming directly from the dual gauge theory. We also constrain stringy corrections to the internal energy.

preprint2016arXiv

Supergravity from D0-brane Quantum Mechanics

The gauge/gravity duality conjecture claims the equivalence between gauge theory and superstring/M-theory. In particular, the one-dimensional gauge theory of D0-branes and type IIA string theory should agree on properties of hot black holes. Type IIA superstring theory predicts the leading $N^2$ behavior of the black hole internal energy to be $E/N^2=a_0T^{14/5}+ a_1T^{23/5}+a_2T^{29/5}+\cdots$ with the supergravity prediction $a_0=7.41$ and unknown coefficients $a_1$, $a_2$, $\ldots$ associated with stringy corrections. In order to test this duality we perform a lattice study of the gauge theory and extract a continuum, large-$N$ value of $a_0=7.4\pm 0.5$---the first direct confirmation of the supergravity prediction at finite temperature---and constrain the stringy corrections ($a_1=-9.7\pm2.2$ and $a_2=5.6\pm1.8$). We also study the sub-leading $1/N^2$ corrections to the internal energy.

preprint2016arXiv

Supergravity from Gauge Theory

Gauge/gravity duality is the conjecture that string theories have dual descriptions as gauge theories. Weakly-coupled gravity is dual to strongly-coupled gauge theories, ideal for lattice calculations. I will show precision lattice calculations that confirm large-N continuum D0-brane quantum mechanics correctly reproduces the leading-order supergravity prediction for a black hole's internal energy---the first leading-order test of the duality---and constrains stringy corrections.

preprint2016arXiv

Two-Nucleon Higher Partial-Wave Scattering from Lattice QCD

We present a determination of nucleon-nucleon scattering phase shifts for l >= 0. The S, P, D and F phase shifts for both the spin-triplet and spin-singlet channels are computed with lattice Quantum ChromoDynamics. For l > 0, this is the first lattice QCD calculation using the Luscher finite-volume formalism. This required the design and implementation of novel lattice methods involving displaced sources and momentum-space cubic sinks. To demonstrate the utility of our approach, the calculations were performed in the SU(3)-flavor limit where the light quark masses have been tuned to the physical strange quark mass, corresponding to m_pi = m_K ~ 800 MeV. In this work, we have assumed that only the lowest partial waves contribute to each channel, ignoring the unphysical partial wave mixing that arises within the finite-volume formalism. This assumption is only valid for sufficiently low energies; we present evidence that it holds for our study using two different channels. Two spatial volumes of V ~ (3.5 fm)^3 and V ~ (4.6 fm)^3 were used. The finite-volume spectrum is extracted from the exponential falloff of the correlation functions. Said spectrum is mapped onto the infinite volume phase shifts using the generalization of the Luscher formalism for two-nucleon systems.

preprint2015arXiv

Direct Detection of Stealth Dark Matter through Electromagnetic Polarizability

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.

preprint2015arXiv

Lattice QCD and Axion Cosmology

The Strong CP Problem can be resolved by introducing an additional global symmetry known as Peccei-Quinn symmetry. Once PQ symmetry is broken the associated particle, the QCD axion, is a plausible dark matter candidate. Calculating the cosmological energy density of the axion requires nonperturbative QCD input---the high-temperature topological susceptibility. I will show results from a pure-glue calculation and examine the implications for the axion mass and coupling.

preprint2015arXiv

Lattice QCD input for axion cosmology

One intriguing BSM particle is the QCD axion, which could simultaneously provide a solution to the Strong CP problem and account for some, if not all, of the dark matter density in the universe. This particle is a pNGB of the conjectured Peccei-Quinn (PQ) symmetry of the Standard Model. Its mass and interactions are suppressed by a heavy symmetry breaking scale, $f_a$, whose value is roughly greater than $10^{9}$ GeV (or, conversely, the axion mass, $m_a$, is roughly less than $10^4\ μ\text{eV}$). The density of axions in the universe, which cannot exceed the relic dark matter density and is a quantity of great interest in axion experiments like ADMX, is a result of the early-universe interplay between cosmological evolution and the axion mass as a function of temperature. The latter quantity is proportional to the second derivative of the QCD free energy with respect to the CP-violating phase, $θ$. However, this quantity is generically non-perturbative and previous calculations have only employed instanton models at the high temperatures of interest (roughly 1 GeV). In this and future works, we aim to calculate the temperature-dependent axion mass at small $θ$ from first-principle lattice calculations, with controlled statistical and systematic errors. Once calculated, this temperature-dependent axion mass is input for the classical evolution equations of the axion density of the universe. Due to a variety of lattice systematic effects at the very high temperatures required, we perform a calculation of the leading small-$θ$ cumulant of the theta vacua on large volume lattices for SU(3) Yang-Mills with high statistics as a first proof of concept, before attempting a full QCD calculation in the future. From these pure glue results, the misalignment mechanism yields the axion mass bound $m_a \geq (14.6\pm0.1) \ μ\text{eV}$ when PQ-breaking occurs after inflation.

preprint2015arXiv

Nuclear Parity Violation from Lattice QCD

The electroweak interaction at the level of quarks and gluons are well understood from precision measurements in high energy collider experiments. Relating these fundamental parameters to Hadronic Parity Violation in nuclei however remains an outstanding theoretical challenge. One of the most interesting observables in this respect is the parity violating hadronic neutral current: it is hard to measure in collider experiments and is thus the least constrained observable of the Standard Model. Precision measurements of parity violating transitions in nuclei can help to improve these constraints. In these systems however, the weak interaction is masked by effects of the seven orders of magnitude stronger non-perturbative strong interaction. Therefore, in order to relate experimental measurements of the parity violating pion-nucleon couplings to the fundamental Lagrangian of the SM, these non-perturbative effects have to be well understood. In this paper, we are going to present a Lattice QCD approach for computing the $ΔI{=}2$ parity violating matrix element in proton proton scattering. This process does not involve disconnected diagrams in the isospin symmetric limit and is thus a perfect testbed for studying the feasibility of the more involved calculation of the parity violating pion-nucleon coupling.

preprint2015arXiv

Two-nucleon scattering in multiple partial waves

We determine scattering phase shifts for S,P,D, and F partial wave channels in two-nucleon systems using lattice QCD methods. We use a generalization of Luscher's finite volume method to determine infinite volume phase shifts from a set of finite volume ground- and excited-state energy levels on two volumes, V=(3.4 fm)^3 and V=(4.5 fm)^3. The calculations are performed in the SU(3)-flavor limit, corresponding to a pion mass of approximately 800 MeV. From the energy dependence of the phase shifts we are able to extract scattering parameters corresponding to an effective range expansion.

preprint2014arXiv

Composite bosonic baryon dark matter on the lattice: SU(4) baryon spectrum and the effective Higgs interaction

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.

preprint2013arXiv

Multi-channel S-matrices from energy levels in finite boxes

We show that for a generic quantum mechanical system with more than one open scattering channel, it is not possible to fully reconstruct the theory's S-matrix from spectral information obtained in large finite volumes with periodic boundary conditions. Physically distinct S-matrices can have identical finite-volume spectra for large finite boxes of arbitrary sizes. If the theory is not time-reversal symmetric, there exists an uncountably infinite set of distinct S-matrices with the same spectra. If the theory respects time-reversal symmetry there exists a discrete set of S-matrices with identical energy levels for finite boxes. We illustrate the issue for simple quantum mechanical systems in 1+1 dimensions.

preprint2012arXiv

Neutrino Emission from Helium White Dwarfs with Condensed Cores

The possibility that ions in a helium white dwarf star are in a Bose-Einstein condensed state has been explored recently. In particular, it has been argued that the resulting novel quantum liquid has a new kind of quasiparticle excitation with a phonon-like dispersion relation. We investigate the neutrino emission rate due to this gapless state and the resulting impact on the total luminosity of helium white dwarf stars, as a possible observable way of detecting this exotic phase. If the condensation temperature for the quantum liquid state, which is currently not known very precisely, turns out to be high enough, our calculations indicate that neutrino emission due to the gapless mode would make a large contribution to the total luminosity of the helium white dwarf stars.

preprint2012arXiv

Thermodynamics of nuclear condensates and phase transitions in white dwarfs

We study the thermodynamics of helium at densities relevant for white dwarf physics. We find evidence that, as the temperature is increased, there is first a first order transition between two superconducting phases followed by a second order transition to the normal state. These transitions occur, for realistic densities, at temperatures below the crystallization temperature and the crystalline state is likely to remain as the true ground state of the system. The calculations are performed with a screening but non-dynamical electron background and we comment on the impact of this and other approximations to our result.

preprint2011arXiv

Nuclear condensate and helium white dwarfs

We consider a high density region of the helium phase diagram, where the nuclei form a Bose-Einstein condensate rather than a classical plasma or a crystal. Helium in this phase may be present in helium-core white dwarfs. We show that in this regime there is a new gapless quasiparticle not previously noticed, arising when the constraints imposed by gauge symmetry are taken into account. The contribution of this quasiparticle to the specific heat of a white dwarf core turns out to be comparable in a range of temperatures to the contribution from the particle-hole excitations of the degenerate electrons. The specific heat in the condensed phase is two orders of magnitude smaller than in the uncondensed plasma phase, which is the ground state at higher temperatures, and four orders of magnitude smaller than the specific heat that an ion lattice would provide, if formed. Since the specific heat of the core is an important input for setting the rate of cooling of a white dwarf star, it may turn out that such a change in the thermal properties of the cores of helium white dwarfs has observable implications.

preprint2011arXiv

Vortons in dense quark matter

At large baryon number density, it is likely that the ground state of QCD is a color-flavor-locked phase with a K0 condensate. The CFL+K0 phase is known to support superconducting vortex strings, and it has been previously suggested that it may also support vortons, which are superconducting vortex rings. We reexamine the question of the stability of vortons, taking into account electromagnetic effects, which make leading-order contributions to vorton dynamics but were not investigated in previous work. We find that current-carrying and electrically charged vortons can be stabilized either by their angular momentum, by Coulomb repulsion, or by a combination of both effects.