BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Kostas Orginos

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

ResearcherUnclaimed source record
hep-lathep-phnucl-thhep-exhep-thnucl-exphysics.comp-phDistributed, Parallel, and Cluster Computingmath.NANumerical Analysis

Catalog footprint

What is connected

52works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

52 published item(s)

preprint2026arXiv

The Spectrum and Scale Setting on 2+1-flavor NME Lattices

This paper describes the thirteen ensembles, named NME, generated with 2+1-flavor Wilson-clover fermions by the JLab/W\&M/LANL/MIT/Marseille collaborations, and presents an analysis of the meson and baryon spectrum, decay constants $f_π$ and $f_K$, flow scales $t_0$ and $w_0$, and time histories of the $Θ$ and Weinberg operators under gradient flow. Using these quantities, the physical point values of the two flow scales, ${t_0^{\rm Phy}}$ and ${w_0^{\rm Phy}}$, and the ratio $\mathop{f_K / f_π}^{\rm Phy}$ are determined. The masses of the octet and decuplet baryons are analyzed using both the next-to-leading order (NLO) and the next-next-to-leading order (NNLO) ansatz from heavy baryon chiral perturbation theory (HB$χ$PT). The NNLO fit to the octet baryons, $M_N$, $M_Σ$, $M_Λ$ and $M_Ξ$, is preferred while the corresponding fits to the decuplet Omega mass, $M_Ω$, are not distinguished. We also present a study of the autocorrelations in the data and show that there is no evidence, even at large flow time, of the freezing of the topological charge or the Weinberg three-gluon operator.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Lattice QCD Calculations of Parton Physics

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The transversity parton distribution function of the nucleon using the pseudo-distribution approach

We present a determination of the non-singlet transversity parton distribution function (PDF) of the nucleon, normalized with respect to the tensor charge at $μ^2=2$ GeV$^2$ from lattice quantum chromodynamics. We apply the pseudo-distribution approach, using a gauge ensemble with a lattice spacing of 0.094 fm and the light quark mass tuned to a pion mass of 358 MeV. We extract the transversity PDF from the analysis of the short-distance behavior of the Ioffe-time pseudo-distribution using the leading-twist next-to-leading order (NLO) matching coefficients calculated for transversity. We reconstruct the $x$-dependence of the transversity PDF through an expansion in a basis of Jacobi polynomials in order to reduce the PDF ansatz dependence. Within the limitations imposed by a heavier-than-physical pion mass and a fixed lattice spacing, we present a comparison of our estimate for the valence transversity PDF with the recent global fit results based on single transverse spin asymmetry. We find the intrinsic nucleon sea to be isospin symmetric with respect to transversity.

0 citations0 reviewsNo code linkedOpen access
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)$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Distillation at High-Momentum

Extraction of hadronic observables at finite-momenta from Lattice QCD (LQCD) is constrained by the well-known signal-to-noise problems afflicting all such LQCD calculations. Traditional quark smearing algorithms are commonly used tools to improve the statistical quality of hadronic $n$-point functions, provided operator momenta are small. The momentum smearing algorithm of Bali et al. extends the range of momenta that are cleanly accessible, and has facilitated countless novel lattice calculations. Momentum smearing has, however, not been explicitly demonstrated within the framework of distillation. In this work we extend the momentum-smearing idea, by exploring a few modifications to the distillation framework. Together with enhanced time slice sampling and expanded operator bases engendered by distillation, we find ground-state nucleon energies can be extracted reliably for $\left|\vec{p}\right|\lesssim3\text{ GeV}$ and matrix elements featuring a large momentum dependence can be resolved.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Parton Distribution Functions from Ioffe time pseudo-distributions

In this paper, we present a detailed study of the unpolarized nucleon parton distribution function (PDF) employing the approach of parton pseudo-distribution functions. We perform a systematic analysis using three lattice ensembles at two volumes, with lattice spacings $a=$ 0.127 fm and $a=$ 0.094 fm, for a pion mass of roughly 400 MeV. With two lattice spacings and two volumes, both continuum limit and infinite volume extrapolation systematic errors of the PDF are estimated. In addition to the $x$ dependence of the PDF, we compute their first two moments and compare them with the pertinent phenomenological determinations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Pion valence quark distribution from current-current correlation in lattice QCD

We extract the pion valence quark distribution $q^π_{\rm v}(x)$ from lattice QCD (LQCD) calculated matrix elements of spacelike correlations of one vector and one axial vector current analyzed in terms of QCD collinear factorization, using a new short-distance matching coefficient calculated to one-loop accuracy. We derive the Ioffe time distribution of the two-current correlations in the physical limit by investigating the finite lattice spacing, volume, quark mass, and higher-twist dependencies in a simultaneous fit of matrix elements computed on four gauge ensembles. We find remarkable consistency between our extracted $q^π_{\rm v}(x)$ and that obtained from experimental data across the entire $x$-range. Further, we demonstrate that the one-loop matching coefficient relating the LQCD matrix computed in position space to the $q_{\rm v}^π(x)$ in momentum space has well-controlled behavior with Ioffe time. This justifies that LQCD calculated current-current correlations are good observables for extracting partonic structures by using QCD factorization, which complements to the global effort to extract partonic structure from experimental data.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Report from the A.I. For Nuclear Physics Workshop

This report is an outcome of the workshop "AI for Nuclear Physics" held at Thomas Jefferson National Accelerator Facility on March 4-6, 2020. The workshop brought together 184 scientists to explore opportunities for Nuclear Physics in the area of Artificial Intelligence. The workshop consisted of plenary talks, as well as six working groups. The report includes the workshop deliberations and additional contributions to describe prospects for using AI across Nuclear Physics research.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Hadrons and Nuclei

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

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Multigrid deflation for Lattice QCD

Computing the trace of the inverse of large matrices is typically addressed through statistical methods. Deflating out the lowest eigenvectors or singular vectors of the matrix reduces the variance of the trace estimator. This work summarizes our efforts to reduce the computational cost of computing the deflation space while achieving the desired variance reduction for Lattice QCD applications. Previous efforts computed the lower part of the singular spectrum of the Dirac operator by using an eigensolver preconditioned with a multigrid linear system solver. Despite the improvement in performance in those applications, as the problem size grows the runtime and storage demands of this approach will eventually dominate the stochastic estimation part of the computation. In this work, we propose to compute the deflation space in one of the following two ways. First, by using an inexact eigensolver on the Hermitian, but maximally indefinite, operator $A γ_5$. Second, by exploiting the fact that the multigrid prolongator for this operator is rich in components toward the lower part of the singular spectrum. We show experimentally that the inexact eigensolver can approximate the lower part of the spectrum even for ill-conditioned operators. Also, the deflation based on the multigrid prolongator is more efficient to compute and apply, and, despite its limited ability to approximate the fine level spectrum, it obtains similar variance reduction on the trace estimator as deflating with approximate eigenvectors from the fine level operator.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Algorithms for Disconnected Diagrams in Lattice QCD

Computing disconnected diagrams in Lattice QCD (operator insertion in a quark loop) entails the computationally demanding problem of taking the trace of the all to all quark propagator. We first outline the basic algorithm used to compute a quark loop as well as improvements to this method. Then, we motivate and introduce an algorithm based on the synergy between hierarchical probing and singular value deflation. We present results for the chiral condensate using a 2+1-flavor clover ensemble and compare estimates of the nucleon charges with the basic algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Axial-Current Matrix Elements in Light Nuclei from Lattice QCD

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Matrix elements from moments of correlation functions

Momentum-space derivatives of matrix elements can be related to their coordinate-space moments through the Fourier transform. We derive these expressions as a function of momentum transfer $Q^2$ for asymptotic in/out states consisting of a single hadron. We calculate corrections to the finite volume moments by studying the spatial dependence of the lattice correlation functions. This method permits the computation of not only the values of matrix elements at momenta accessible on the lattice, but also the momentum-space derivatives, providing {\it a priori} information about the $Q^2$ dependence of form factors. As a specific application we use the method, at a single lattice spacing and with unphysically heavy quarks, to directly obtain the slope of the isovector form factor at various $Q^2$, whence the isovector charge radius. The method has potential application in the calculation of any hadronic matrix element with momentum transfer, including those relevant to hadronic weak decays.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Unitary Limit of Two-Nucleon Interactions in Strong Magnetic Fields

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

High-precision calculation of the strange nucleon electromagnetic form factors

We report a direct lattice QCD calculation of the strange nucleon electromagnetic form factors $G_E^s$ and $G_M^s$ in the kinematic range $0 \leq Q^2 \lesssim 1.2\: {\rm GeV}^2$. For the first time, both $G_E^s$ and $G_M^s$ are shown to be nonzero with high significance. This work uses closer-to-physical lattice parameters than previous calculations, and achieves an unprecedented statistical precision by implementing a recently proposed variance reduction technique called hierarchical probing. We perform model-independent fits of the form factor shapes using the $z$-expansion and determine the strange electric and magnetic radii and magnetic moment. We compare our results to parity-violating electron-proton scattering data and to other theoretical studies.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Improved Sampling Algorithms in Lattice QCD

Reverse Monte Carlo (RMC) is an algorithm that incorporates stochastic modification of the action as part of the process that updates the fields in a Monte Carlo simulation. Such update moves have the potential of lowering or eliminating potential barriers that may cause inefficiencies in exploring the field configuration space. The highly successful Cluster algorithms for spin systems can be derived from the RMC framework. In this work we apply RMC ideas to pure gauge theory, aiming to alleviate the critical slowing down observed in the topological charge evolution as well as other long distance observables. We present various formulations of the basic idea and report on our numerical experiments with these algorithms.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Locally smeared operator product expansions in scalar field theory

We propose a new locally smeared operator product expansion to decompose nonlocal operators in terms of a basis of smeared operators. The smeared operator product expansion formally connects nonperturbative matrix elements determined numerically using lattice field theory to matrix elements of nonlocal operators in the continuum. These nonperturbative matrix elements do not suffer from power-divergent mixing on the lattice, which significantly complicates calculations of quantities such as the moments of parton distribution functions, provided the smearing scale is kept fixed in the continuum limit. The presence of this smearing scale complicates the connection to the Wilson coefficients of the standard operator product expansion and requires the construction of a suitable formalism. We demonstrate the feasibility of our approach with examples in real scalar field theory.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Multiscale Monte Carlo equilibration: Pure Yang-Mills theory

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Magnetic Structure of Light Nuclei from Lattice QCD

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Charmed bottom baryon spectroscopy from lattice QCD

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Finite volume renormalization scheme for fermionic operators

We propose a new finite volume renormalization scheme. Our scheme is based on the Gradient Flow applied to both fermion and gauge fields and, much like the Schrödinger functional method, allows for a nonperturbative determination of the scale dependence of operators using a step-scaling approach. We give some preliminary results for the pseudo-scalar density in the quenched approximation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Locally-smeared operator product expansions

We propose a "locally-smeared Operator Product Expansion" (sOPE) to decompose non-local operators in terms of a basis of locally-smeared operators. The sOPE formally connects nonperturbative matrix elements of smeared degrees of freedom, determined numerically using the gradient flow, to non-local operators in the continuum. The nonperturbative matrix elements do not suffer from power-divergent mixing on the lattice, provided the smearing scale is kept fixed in the continuum limit. The presence of this smearing scale prevents a simple connection to the standard operator product expansion and therefore requires the construction of a two-scale formalism. We demonstrate the feasibility of our approach using the example of real scalar field theory.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Möbius Domain Wall Fermion Algorithm

We present a review of the properties of generalized domain wall Fermions, based on a (real) Möbius transformation on the Wilson overlap kernel, discussing their algorithmic efficiency, the degree of explicit chiral violations measured by the residual mass ($m_{res}$) and the Ward-Takahashi identities. The Möbius class interpolates between Shamir's domain wall operator and Boriçi's domain wall implementation of Neuberger's overlap operator without increasing the number of Dirac applications per conjugate gradient iteration. A new scaling parameter ($α$) reduces chiral violations at finite fifth dimension ($L_s$) but yields exactly the same overlap action in the limit $L_s \rightarrow \infty$. Through the use of 4d Red/Black preconditioning and optimal tuning for the scaling $α(L_s)$, we show that chiral symmetry violations are typically reduced by an order of magnitude at fixed $L_s$. At large $L_s$ we argue that the observed scaling for $m_{res} = O(1/L_s)$ for Shamir is replaced by $m_{res} = O(1/L_s^2)$ for the properly tuned Möbius algorithm with $α= O(L_s)$

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Uncertainty Quantification in Lattice QCD Calculations for Nuclear Physics

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Hierarchical probing for estimating the trace of the matrix inverse on toroidal lattices

The standard approach for computing the trace of the inverse of a very large, sparse matrix $A$ is to view the trace as the mean value of matrix quadratures, and use the Monte Carlo algorithm to estimate it. This approach is heavily used in our motivating application of Lattice QCD. Often, the elements of $A^{-1}$ display certain decay properties away from the non zero structure of $A$, but random vectors cannot exploit this induced structure of $A^{-1}$. Probing is a technique that, given a sparsity pattern of $A$, discovers elements of $A$ through matrix-vector multiplications with specially designed vectors. In the case of $A^{-1}$, the pattern is obtained by distance-$k$ coloring of the graph of $A$. For sufficiently large $k$, the method produces accurate trace estimates but the cost of producing the colorings becomes prohibitively expensive. More importantly, it is difficult to search for an optimal $k$ value, since none of the work for prior choices of $k$ can be reused.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Lattice QCD at non-zero isospin chemical potential

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Nuclear correlation functions in lattice QCD

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Tetraquark bound states in the heavy-light heavy-light system

A calculation of the interaction potential of two heavy-light mesons in lattice QCD is used to study the existence of tetraquark bound states. The interaction potential of the tetraquark system is calculated on the lattice with 2+1 flavours of dynamical fermions with lattice interpolating fields constructed using colorwave propagators. These propagators provide a new method for constructing all-to-all spatially smeared the interpolating fields, a technique which allows for a better overlap with the ground state wavefunction as well as reduced statistical noise. Potentials are extracted for 24 distinct channels, and are fit with a phenomenological non-relativistic quark model potential, from which a determination of the existence of bound states is made via numerical solution of the two body radial Schrödinger equation.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Application of Quadrature Methods for Re-Weighting in Lattice QCD

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Heavy-Baryon Spectroscopy from Lattice QCD

We use a four-dimensional lattice calculation of the full-QCD (quantum chromodynamics, the non-abliean gauge theory of the strong interactions of quarks and gluons) path integrals needed to determine the masses of the charmed and bottom baryons. In the charm sector, our results are in good agreement with experiment within our systematics, except for the spin-1/2 $Ξ_{cc}$, for which we found the isospin-averaged mass to be $Ξ_{cc}$ to be $3665\pm17\pm14^{+0}_{-78}$ MeV. We predict the mass of the (isospin-averaged) spin-1/2 $Ω_{cc}$ to be $3763\pm19\pm26^{+13}_{-79}$ {MeV}. In the bottom sector, our results are also in agreement with experimental observations and other lattice calculations within our statistical and systematic errors. In particular, we find the mass of the $Ω_b$ to be consistent with the recent CDF measurement. We also predict the mass for the as yet unobserved $Ξ^\prime_b$ to be 5955(27) MeV.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Lattice Calculations of Nucleon Electromagnetic Form Factors at Large Momentum Transfer

In this work, we report a novel technique in lattice QCD for studying the high momentum-transfer region of nucleon form factors. These calculations could give important theoretical input to experiments, such as those of JLab's 12-GeV program and studies of nucleon deformation. There is an extensive history of form-factor calculations on the lattice, primarily with ground states for both the initial and final state. However, determining form factors at large momentum transfer ($Q^2$) has been difficult due to large statistical and systematic errors in this regime. We study the nucleon form factors using three pion masses with both quenched and 2+1-flavor anisotropic lattice configurations with $Q^2$ as large as $6 GeV^2$. These form factors are further processed to obtain transverse charge and magnetization densities across 2-dimensional impact-parameter space. Our approach can be applied to isotropic lattices and lattices with smaller lattice spacing to calculate even larger-$Q^2$ form factors.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Nucleon structure from mixed action calculations using 2+1 flavors of asqtad sea and domain wall valence fermions

We present high statistics results for the structure of the nucleon from a mixed-action calculation using 2+1 flavors of asqtad sea and domain wall valence fermions. We perform extrapolations of our data based on different chiral effective field theory schemes and compare our results with available information from phenomenology. We discuss vector and axial form factors of the nucleon, moments of generalized parton distributions, including moments of forward parton distributions, and implications for the decomposition of the nucleon spin.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Extending the eigCG algorithm to non-symmetric linear systems with multiple right-hand sides

For Hermitian positive definite linear systems and eigenvalue problems, the eigCG algorithm is a memory efficient algorithm that solves the linear system and simultaneously computes some of its eigenvalues. The algorithm is based on the Conjugate-Gradient (CG) algorithm, however, it uses only a window of the vectors generated by the CG algorithm to compute approximate eigenvalues. The number and accuracy of the eigenvectors can be increased by solving more right-hand sides. For Hermitian systems with multiple right-hand sides, the computed eigenvectors can be used to speed up the solution of subsequent systems. The algorithm was tested on Lattice QCD problems by solving the normal equations and was shown to give large speed up factors and to remove the critical slowing down as we approach light quark masses. Here, an extension to the non-symmetric case based on the two-sided Lanczos algorithm is given. The new algorithm is tested on Lattice QCD problems and is shown to give promising results. We also study the removal of the critical slowing down and compare results with those of the eigCG algorithm. We also discuss the case when the system is gamma5-Hermitian.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Meson-Baryon Scattering Lengths from Mixed-Action Lattice QCD

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

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Singly and Doubly Charmed J=1/2 Baryon Spectrum from Lattice QCD

We compute the masses of the singly and doubly charmed baryons in full QCD using the relativistic Fermilab action for the charm quark. For the light quarks we use domain-wall fermions in the valence sector and improved Kogut-Susskind sea quarks. We use the low-lying charmonium spectrum to tune our heavy-quark action and as a guide to understanding the discretization errors associated with the heavy quark. Our results are in good agreement with experiment within our systematics. For the Ξ_{cc}, we find the isospin-averaged mass to be M_{Ξ_{cc}} = 3665 +/- 17 +/- 14 +0 -78 MeV; the three given uncertainties are statistical, systematic and an estimate of lattice discretization errors, respectively. In addition, we predict the mass splitting of the (isospin-averaged) spin-1/2 Ω_{cc} with the Ξ_{cc} to be M_{Ω_{cc}} - M_{Ξ_{cc}} = 98 +/- 9 +/- 22 +/- 13 MeV (in this mass splitting, the leading discretization errors are also suppressed by SU(3) symmetry). Combining this splitting with our determination of M_{Ξ_{cc}} leads to our prediction of the spin-1/2 Ω_{cc} mass, M_{Ω_{cc}} = 3763 +/- 19 +/- 26 +13 -79 MeV.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Strange Baryon Electromagnetic Form Factors and SU(3) Flavor Symmetry Breaking

We study the nucleon, Sigma and cascade octet baryon electromagnetic form factors and the effects of SU(3) flavor symmetry breaking from 2+1-flavor lattice calculations. We find that electric and magnetic radii are similar; the maximum discrepancy is about 10%. In the pion-mass region we explore, both the quark-component and full-baryon moments have small SU(3) symmetry breaking. We extrapolate the charge radii and the magnetic moments using three-flavor heavy-baryon chiral perturbation theory (HBXPT). The systematic errors due to chiral and continuum extrapolations remain significant, giving rise to charge radii for $p$ and $Σ^-$ that are 3--4 standard deviations away from the known experimental ones. Within these systematics the predicted $Σ^+$ and $Ξ^-$ radii are 0.67(5) and 0.306(15) fm$^2$ respectively. When the next-to-next-to-leading order of HBXPT is included, the extrapolated magnetic moments are less than 3 standard deviations away from PDG values, and the discrepancy is possibly due to remaining chiral and continuum extrapolation errors.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Computing and deflating eigenvalues while solving multiple right hand side linear systems in Quantum Chromodynamics

We present a new algorithm that computes eigenvalues and eigenvectors of a Hermitian positive definite matrix while solving a linear system of equations with Conjugate Gradient (CG). Traditionally, all the CG iteration vectors could be saved and recombined through the eigenvectors of the tridiagonal projection matrix, which is equivalent theoretically to unrestarted Lanczos. Our algorithm capitalizes on the iteration vectors produced by CG to update only a small window of vectors that approximate the eigenvectors. While this window is restarted in a locally optimal way, the CG algorithm for the linear system is unaffected. Yet, in all our experiments, this small window converges to the required eigenvectors at a rate identical to unrestarted Lanczos. After the solution of the linear system, eigenvectors that have not accurately converged can be improved in an incremental fashion by solving additional linear systems. In this case, eigenvectors identified in earlier systems can be used to deflate, and thus accelerate, the convergence of subsequent systems. We have used this algorithm with excellent results in lattice QCD applications, where hundreds of right hand sides may be needed. Specifically, about 70 eigenvectors are obtained to full accuracy after solving 24 right hand sides. Deflating these from the large number of subsequent right hand sides removes the dreaded critical slowdown, where the conditioning of the matrix increases as the quark mass reaches a critical value. Our experiments show almost a constant number of iterations for our method, regardless of quark mass, and speedups of 8 over original CG for light quark masses.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Hadron structure with light dynamical quarks

Generalized parton distributions encompass a wealth of information concerning the three-dimensional quark and gluon structure of the nucleon, and thus provide an ideal focus for the study of hadron structure using lattice QCD. The special limits corresponding to form factors and parton distributions are well explored experimentally, providing clear tests of lattice calculations, and the lack of experimental data for more general cases provides opportunities for genuine predictions and for guiding experiment. We present results from hybrid calculations with improved staggered (Asqtad) sea quarks and domain wall valence quarks at pion masses down to 350 MeV.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Nucleon structure functions with domain wall fermions

We present a quenched lattice QCD calculation of the first few moments of the polarized and un-polarized structure functions of the nucleon. Our calculations are done using domain wall fermions and the DBW2 gauge action with inverse lattice spacing ~1.3GeV, physical volume approximatelly (2.4 fm)^3, and light quark masses down to about 1/4 the strange quark mass. Values of the individual moments are found to be significantly larger than experiment, as in past lattice calculations, but interestingly the chiral symmetry of domain wall fermions allows for a precise determination of the ratio of the flavor non-singlet momentum fraction to the helicity distribution, which is in very good agreement with experiment. We discuss the implications of this result. Next, we show that the chiral symmetry of domain wall fermions is useful in eliminating mixing of power divergent lower dimensional operators with twist-3 operators. Finally, we find the isovector tensor charge at renormalization scale 2 GeV in the MS bar scheme to be 1.192(30), where the error is the statistical error only.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Deconfining Phase Transition as a Matrix Model of Renormalized Polyakov Loops

We discuss how to extract renormalized from bare Polyakov loops in SU(N) lattice gauge theories at nonzero temperature in four spacetime dimensions. Single loops in an irreducible representation are multiplicatively renormalized without mixing, through a renormalization constant which depends upon both representation and temperature. The values of renormalized loops in the four lowest representations of SU(3) were measured numerically on small, coarse lattices. We find that in magnitude, condensates for the sextet and octet loops are approximately the square of the triplet loop. This agrees with a large $N$ expansion, where factorization implies that the expectation values of loops in adjoint and higher representations are just powers of fundamental and anti-fundamental loops. For three colors, numerically the corrections to the large $N$ relations are greatest for the sextet loop, $\leq 25%$; these represent corrections of $\sim 1/N$ for N=3. The values of the renormalized triplet loop can be described by an SU(3) matrix model, with an effective action dominated by the triplet loop. In several ways, the deconfining phase transition for N=3 appears to be like that in the $N=\infty$ matrix model of Gross and Witten.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Mobius Fermions: Improved Domain Wall Chiral Fermions

A new class of domain wall fermions is defined that interpolates between Shamir's and Boriçi's form without increasing the number of Dirac applications per CG iteration. This class represents a full (real) Möbius transformation of the Wilson kernel. Simulations on quenched Wilson lattices with $β= 6.0$ show that the number of lattice sites ($L_s$) in the fifth dimension can be reduced by a factor of 2 or more at fixed value of chiral symmetry violations measured by the residual mass ($m_{res}$).

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Variants of fattening and flavor symmetry restoration

We study the effects of different "fat link" actions for Kogut-Susskind quarks on flavor symmetry breaking. Our method is mostly empirical - we compute the pion spectrum with different valence quark actions on common sets of sample lattices. Different actions are compared, as best we can, at equivalent physical points. We find significant reductions in flavor symmetry breaking relative to the conventional or to the "link plus staple" actions, with a reasonable cost in computer time. We also develop and test a scheme for approximate unitarization of the fat links. While our tests have concentrated on the valence quark action, our results will be useful in designing simulations with dynamical quarks.

0 citations0 reviewsNo code linkedOpen access