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

Rajan Gupta

Rajan Gupta 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-thquant-phhep-thcond-mat.str-elnucl-excond-mathep-exPopulations and Evolution

Catalog footprint

What is connected

64works
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

64 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
preprint2024arXiv

Update on flavor diagonal nucleon charges from clover fermions

We present a summary of the full calculation of the axial, scalar and tensor flavor diagonal charges of the nucleon carried out using Wilson-clover fermions on eight ensembles generated using 2+1+1-flavors of highly improved staggered quarks (HISQ) by the MILC collaboration. We also give results for the $3\times 3$ matrix of renormalization factors between the RI-sMOM and $\overline{\rm MS}$ scheme for the 2+1 flavor theory that include flavor mixing. Preliminary results for $g_{A,S,T}^{u,d,s}$ are presented in the $\overline{\rm MS}$ scheme at scale 2 GeV.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

All About the Neutron from Lattice QCD

I describe how simulations of lattice QCD using the path integral formulation provide the two basic quantum mechanical properties of QCD, its ground state in which correlation functions are calculated, and Fock state wavefunctions between which matrix elements of operators are calculated. Both constructs are stochastic, so unfortunately one gets no intuitive picture or even a qualitative understanding of what they look like, nevertheless they contain and display all the subtleties of the quantum field theory. Today, these simulations provide many quantities that are impacting phenomenology and experiments. I illustrate the methods and the steps in the analysis using, as examples, three observables: the isovector charges of the nucleon, the contribution of the quark's intrinsic spin to the nucleon spin, and the pion-nucleon sigma term.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

nEDM from the theta-term and chromoEDM operators

In a previous work, we showed that unresolved excited state contaminations provide a major source of systematic uncertainty in the calculation of the nucleon electric dipole moment due to the QCD topological term theta. Here we extend the calculation to the quark chromo-electric dipole moment operator (qcEDM). We also show quantitatively the impact of the mixing of the qcEDM with lower-dimensional operators on the lattice. Finally, we present preliminary results from a unitary clover-on-clover calculation for the QCD topological term.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Nucleon form factors and the pion-nucleon sigma term

This talk summarizes the progress made since Lattice 2021 in understanding and controlling the contributions of towers of multihadron excited states with mass gaps starting lower than of radial excitations, and in increasing our confidence in the extraction of ground state nucleon matrix elements. The most clear evidence for multihadron excited state contributions (ESC) is in axial/pseudoscalar form factors that are required to satisfy the PCAC relation between them. The talk examines the broader question--which and how many of the theoretically allowed positive parity states $N(\textbf p)π(-\textbf p)$, $N(\textbf 0)π(\textbf 0)π(\textbf 0)$, $N(\textbf p)π(\textbf 0)$, $N(\textbf 0)π(\textbf p),\ \ldots$ make significant contributions to a given nucleon matrix element? New data for the axial, electric and magnetic form factors are presented. They continue to show trends observed in Ref[1]. The N${}^2$LO $χ$PT analysis of the ESC to the pion-nucleon sigma term, $σ_{πN}$, has been extended to include the $Δ$ as an explicit degree of freedom [2]. The conclusion reached in Ref [3] that $N π$ and $N ππ$ states each contribute about 10 MeV to $σ_{πN}$, and the consistency between the lattice result with $N π$ state included and the phenomenological estimate is not changed by this improvement.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Phases of 2d massless QCD with qubit regularization

We investigate the possibility of reproducing the continuum physics of 2d SU(N) gauge theory coupled to a single flavor of massless Dirac fermions using qubit regularization. The continuum theory is described by N free fermions in the ultraviolet (UV) and a coset Wess-Zumino-Witten (WZW) model in the infrared (IR). In this work, we explore how well these features can be reproduced using the Kogut-Susskind Hamiltonian with a finite-dimensional link Hilbert space and a generalized Hubbard coupling. Using strong coupling expansions, we show that our model exhibits a gapped dimer phase and another phase described by a spin-chain. Furthermore, for N=2, using tensor network methods, we show that there is a second-order phase transition between these two phases. The critical theory at the transition can be understood as an SU(2)_1 WZW model, using which we determine the phase diagram of our model quantitatively. Using the confinement properties of the model we argue how the UV physics of free fermions could also emerge, but may require further modifications to our model.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Synergy between NP and HEP research goals and efforts in fundamental symmetries and interactions

The aim of this white paper is to highlight several areas for which the Department of Energy's Office of Nuclear Physics has primary stewardship or significant investment and expertise, and for which there is also significant interest and expertise within the HEP community. These areas of overlap offer exciting opportunities for collaboration.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Update on Flavor Diagonal Nucleon Charges

This talk provides an update on the calculation of matrix elements of flavor diagonal axial, scalar and tensor quark bilinear operators between the nucleon ground state. The simulations are done using Wilson-clover fermions on a sea of eight 2+1+1-flavor HISQ ensembles generated by the MILC collaboration. We discuss the signal in the sum of the connected and disconnected contributions for the up, down and strange quarks, control over fits to remove excited state contamination, and the simultaneous chiral-continuum fit used to extract the charges.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

AI and Theoretical Particle Physics

Theoretical particle physicists continue to push the envelope in both high performance computing and in managing and analyzing large data sets. For example, the goals of sub-percent accuracy in predictions of quantum chromodynamics (QCD) using large scale simulations of lattice QCD and in finding signals of rare events and new physics in exabytes of data produced by experiments at the high luminosity large hadron collider (LHC) require new tools beyond just developments in hardware. Machine learning and artificial intelligence offer the promise of dramatically reducing the computational cost and time. This chapter reviews selected areas where AI/ML tools could have a major impact, provides an overview of the challenges, and discusses how new ideas such as normalizing flows can speed up the generation of gauge configurations needed in lattice QCD calculations; the growth of ML in surrogate models and pattern matching to reduce the cost of event generators and in the analysis of experimental data; and in the search for viable vacua in the landscape of string theories. While such approaches transform aspects of particle theory into computational problems, and thus black boxes, we argue that physics-aware development of these tools combined with algorithms that ensure that the results are bias free will continue to require a deep understanding of the physics. We see this broader transformation as akin to formulating and extracting observables from simulations of lattice QCD, a numerical integration of the path integral formulation of QCD that nevertheless requires a deep understanding of the underlying quantum field theory, the standard model of particle physics and effective field theory methods.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Calculation of neutron electric dipole moment due to the QCD topological term, Weinberg three-gluon operator and the quark chromoelectric moment

We present results for the neutron electric dipole moment due to the dimension 4 and dimension 6 gluonic CP violation, and the isovector quark chromoelectric dipole moment using clover valence quarks on HISQ dynamical ensembles generated by the MILC Collaboration. For the gluonic operators, we use the gradient flow scheme to obtain divergence-free continuum results. For the chromoelectric dipole moment operator, we use the unflowed local operator but discuss how the quadratically divergent mixing with the pseudoscalar operator can be controlled nonperturbatively.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Excited states and precision results for nucleon charges and form factors

The exponentially falling signal-to-noise ratio in all nucleon correlation functions, and the presence of towers of multihadron excited states with relatively small mass gaps makes extraction of matrix elements of various operators within the ground state nucleon challenging. Theoretically, the allowed positive parity states with the smallest mass gaps are the $N(\bm p)π(-\bm p)$, $N(\bm 0)π(\bm 0)π(\bm 0)$, $N(\bm p)π(\bm 0)$, $N(\bm 0)π(\bm p),\ \ldots$, states. A priori, the contribution of these states arises at one loop in chiral perturbation theory ($χ$PT), however, in many cases the contributions are enhanced. In this talk, I will review four such cases: the correlation functions from which the axial form factors, electric and magnetic form factors, the $Θ$-term contribution to neutron electric dipole moment (nEDM), and the pion-nucleon sigma term are extracted. Including appropriate multihadron states in the analysis can lead to significantly different results compared to standard analyses with the mass gaps taken from fits to 2-point functions. The $χ$PT case for $N π$ states is the most clear in the axial/pseudoscalar form factors which need to satisfy the PCAC relation between them. Our analyses, supported by $χ$PT, suggests similarly large effects in the calculations of the $Θ$-term and the pion-nucleon sigma term that have significant phenomenological implications.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Flavor Diagonal Nucleon Charges

This talk provides an update on the calculation of matrix elements of flavor diagonal axial, scalar and tensor quark bilinear operators between the nucleon ground state. The simulations are done using Wilson-clover fermions on a sea of eight 2+1+1-flavor HISQ ensembles generated by the MILC collaboration. We discuss the signal in the connected and disconnected contributions, calculation of the renormalization constants and mixing in the RI-sMOM scheme, and control over the simultaneous chiral-continuum-finite-volume fit used to extract the final charges.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Precision Nucleon Charges and Form Factors Using 2+1-flavor Lattice QCD

We present high statistics results for the isovector nucleon charges and form factors using seven ensembles of 2+1-flavor Wilson-clover fermions. The axial and pseudoscalar form factors obtained on each ensemble satisfy the PCAC relation once the lowest energy $Nπ$ excited state is included in the spectral decomposition of the correlation functions used for extracting the ground state matrix elements. Similarly, we find evidence that the $Nππ$ excited state contributes to the correlation functions with the vector current, consistent with the vector meson dominance model. The resulting form factors are consistent with the Kelly parameterization of the experimental electric and magnetic data. Our final estimates for the isovector charges are $g_{A}^{u-d} = 1.31(06)(05)_{sys}$, $g_{S}^{u-d} = 1.06(10)(06)_{sys}$, and $g_{T}^{u-d} = 0.95(05)(02)_{sys}$, where the first error is the overall analysis uncertainty and the second is an additional combined systematic uncertainty. The form factors yield: (i) the axial charge radius squared, ${\langle r_A^2 \rangle}^{u-d}=0.428(53)(30)_{sys}\ {\rm fm}^2$, (ii) the induced pseudoscalar charge, $g_P^\ast=7.9(7)(9)_{sys}$, (iii) the pion-nucleon coupling $g_{π{\rm NN}} = 12.4(1.2)$, (iv) the electric charge radius squared, ${\langle r_E^2 \rangle}^{u-d} = 0.85(12)(19)_{sys} \ {\rm fm}^2$, (v) the magnetic charge radius squared, ${\langle r_M^2 \rangle}^{u-d} = 0.71(19)(23)_{\rm sys} \ {\rm fm}^2$, and (vi) the magnetic moment $μ^{u-d} = 4.15(22)(10)_{\rm sys}$. All our results are consistent with phenomenological/experimental values but with larger errors. Lastly, we present a Padé parameterization of the axial, electric and magnetic form factors over the range $0.04< Q^2 <1$ GeV${}^2$ for phenomenological studies.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The pion-nucleon sigma term from Lattice QCD

We summarize recent evidence, both from lattice QCD and chiral perturbation theory, that suggests that larger-than-expected excited-state contamination could be the reason for the tension between phenomenological determinations and previous direct lattice-QCD calculations of the pion--nucleon sigma term $σ_{πN}$. In addition, we extend the $χ$PT analysis by calculating the corrections due to including the $Δ(1232)$ resonance as an explicit degree of freedom. This correction is found to be small, thereby corroborating the excited-state effects found in the $Δ$-less calculation and the result for $σ_{πN}$.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Moments of nucleon isovector structure functions in $2+1+1$-flavor QCD

We present results on the isovector momentum fraction, $\langle x \rangle_{u-d}$, helicity moment, $\langle x \rangle_{Δu-Δd}$, and the transversity moment, $\langle x \rangle_{δu-δd}$, of the nucleon obtained using nine ensembles of gauge configurations generated by the MILC collaboration using $2+1+1$-flavors of dynamical highly improved staggered quarks (HISQ). The correlation functions are calculated using the Wilson-Clover action and the renormalization of the three operators is carried out nonperturbatively on the lattice in the RI${}^\prime$-MOM scheme. The data have been collected at lattice spacings $a \approx 0.15,\ 0.12,\ 0.09,$ and 0.06 fm and $M_π\approx 310,\ 220$ and 135 MeV, which are used to obtain the physical values using a simultaneous chiral-continuum-finite-volume fit. The final results, in the $\overline{MS}$ scheme at 2 GeV, are $\langle x \rangle_{u-d} = 0.173(14)(07)$, $\langle x \rangle_{Δu-Δd} = 0.213(15)(22)$ and $\langle x \rangle_{δu-δd} = 0.208(19)(24)$, where the first error is the overall analysis uncertainty and the second is an additional systematic uncertainty due to possible residual excited-state contributions. These results are consistent with other recent lattice calculations and phenomenological global fit values.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Nuclear two point correlation functions on a quantum-computer

The calculation of dynamic response functions is expected to be an early application benefiting from rapidly developing quantum hardware resources. The ability to calculate real-time quantities of strongly-correlated quantum systems is one of the most exciting applications that can easily reach beyond the capabilities of traditional classical hardware. Response functions of fermionic systems at moderate momenta and energies corresponding roughly to the Fermi energy of the system are a potential early application because the relevant operators are nearly local and the energies can be resolved in moderately short real time, reducing the spatial resolution and gate depth required. This is particularly the case in quasielastic electron and neutrino scattering from nuclei, a topic of great interest in the nuclear and particle physics communities and directly related to experiments designed to probe neutrino properties. In this work we use current quantum hardware and error mitigation protocols to calculate response functions for a highly simplified nuclear model through calculations of a 2-point real time correlation function for a modified Fermi-Hubbard model in two dimensions with three distinguishable nucleons on four lattice sites.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Nucleon isovector momentum fraction, helicity and transversity moment using Lattice QCD

We present our recent high precision calculations (Phys. Rev. D102 (2020) no.5, 054512 and JHEP 04 (2021) 044, JHEP 21 (2020) 004) of the first moment of nucleon isovector polarized, unpolarized and transversity distributions, i.e., momentum fraction, helicity and transversity moment, respectively. We use the standard method for the calculation of these moments (via matrix elements of twist two operators), and carry out a detailed analysis of the sources of systematic uncertainty, in particular of excited state contributions. Our calculations have been performed using two different lattice setups (Clover-on-HISQ and Clover-on-Clover), each with several ensembles. They give consistent results that are in agreement with global fit analyses.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Space-time symmetric qubit regularization of the asymptotically free two-dimensional O(4) model

We explore if space-time symmetric lattice field theory models with a finite Hilbert space per lattice site can reproduce asymptotic freedom in the two-dimensional $O(4)$ model. We focus on a simple class of such models with a five dimensional local Hilbert space. We demonstrate how even the simplest model reproduces asymptotic freedom within the D-theory formalism but at the cost of increasing the size of the Hilbert space through coupling several layers of a two-dimensional lattice. We then argue that qubit regularization can be viewed as an effective field theory (EFT) even if the continuum limit cannot be reached, as long as we can tune the model close enough to the continuum limit where perturbation theory, or other analytical techniques, become viable. We construct a simple lattice model on a single layer with a four dimensional local Hilbert space that acts like an excellent EFT of the original theory.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Neutron Electric Dipole Moments with Clover Fermions

We present preliminary results for the contributions to the neutron EDM arising from the QCD $θ$-term, the Weinberg three-gluon and the quark chromo-EDM operators from our ongoing lattice calculations using clover valence quarks on the MILC HISQ lattices. We use the gradient-flow technique to smooth the lattices and renormalize the gluonic operators, and use the Schwinger source method to incorporate the quark chromo-EDM interactions in the quark propagator. For the QCD $θ$-term and the Weinberg three-gluon operator, we report results in the gradient-flow scheme from 8 ensembles at four lattice spacings and three pion masses, including 2 physical pion mass ensembles. For the quark chromo-EDM, unrenormalized results are presented at two lattice spacings, $a=0.12$ and $0.09$ fm, and two pion masses, $M_π= 310$ MeV and $220$ MeV.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Nucleon Axial Form Factors from Clover Fermion on 2+1+1-flavor HISQ Lattice

The nucleon axial form factors -- axial $G_A$, induced pseudoscalar $\widetilde{G}_P$ and pseudoscalar $G_P$ -- have displayed large systematics in lattice QCD calculations. The major symptoms were the violation of the partially conserved axial current (PCAC) relation between the three form factors, and the underestimation of the induced pseudoscalar coupling $g_P^\ast$ and the axial charge radius $r_A$ compared to phenomenological estimates. The small $g_P^\ast$ was a consequence of the failure of the pion-pole dominance (PPD) hypothesis, especially at low $M_π^2$. The small charge radius $r_A$ and the underestimate of $g_A$ were related. The dominant systematic responsible is the lack of inclusion of low-energy ($N π$) states that are not manifest in the multiexponential fit to the nucleon two-point correlator. We show that this low-energy state can be determined from the three-point correlator $\langle N A_4 N \rangle $ with the insertion of the temporal component of the axial current $A_4$ within the nucleon state, ie, the strategy labeled $S_{A4}$ [1]. Including this low-energy state in fits to control excited-state contamination (ESC) gives results for $g_A$, $r_A$, and $g_P^\ast$ that are consistent with experimental/phenomenological values. However, the systematic uncertainties, especially in data at small $Q^2$, are now much larger.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Nucleon charges and form factors using clover and HISQ ensembles

We present high statistics ($\mathcal{O}(2\times 10^5)$ measurements) preliminary results on (i) the isovector charges, $g^{u-d}_{A,S,T}$, and form factors, $G^{u-d}_E(Q^2)$, $G^{u-d}_M(Q^2)$, $G^{u-d}_A(Q^2)$, $\widetilde G^{u-d}_P(Q^2)$, $G^{u-d}_P(Q^2)$, on six 2+1-flavor Wilson-clover ensembles generated by the JLab/W&M/LANL/MIT collaboration with lattice parameters given in Table 1. Examples of the impact of using different estimates of the excited state spectra are given for the clover-on-clover data, and as discussed in [1], the biggest difference on including the lower energy (close to $Nπ$ and $Nππ$) states is in the axial channel. (ii) Flavor diagonal axial, tensor and scalar charges, $g^{u,d,s}_{A,S,T}$, are calculated with the clover-on-HISQ formulation using nine 2+1+1-flavor HISQ ensembles generated by the MILC collaboration [2] with lattice parameters given in Table 2. Once finished, the calculations of $g^{u,d,s}_{A,T}$ will update the results given in Refs.[3,4]. The estimates for $g^{u,d,s}_{S}$ and $σ_{Nπ}$ are new. Overall, a large part of the focus is on understanding the excited state contamination (ESC), and the results discussed provide a partial status report on developing defensible analyses strategies that include contributions of possible low-lying excited states to individual nucleon matrix elements.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Recent results of nucleon structure & matrix element calculations

A review of recent lattice calculations of nucleon structure and matrix elements of operators in nucleons is presented. It primarily covers developments in the calculation of the matrix elements of the scalar, tensor, pseudo-scalar, axial-vector and vector operators in the ground state of neutrons and protons in the isospin symmetric limit. Some preliminary calculations of the electric dipole moment, the gravitational moments and stress-energy distribution, and the magnetic polarizability are briefly described.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

SEIR and Regression Model based COVID-19 outbreak predictions in India

COVID-19 pandemic has become a major threat to the country. Till date, well tested medication or antidote is not available to cure this disease. According to WHO reports, COVID-19 is a severe acute respiratory syndrome which is transmitted through respiratory droplets and contact routes. Analysis of this disease requires major attention by the Government to take necessary steps in reducing the effect of this global pandemic. In this study, outbreak of this disease has been analysed for India till 30th March 2020 and predictions have been made for the number of cases for the next 2 weeks. SEIR model and Regression model have been used for predictions based on the data collected from John Hopkins University repository in the time period of 30th January 2020 to 30th March 2020. The performance of the models was evaluated using RMSLE and achieved 1.52 for SEIR model and 1.75 for the regression model. The RMSLE error rate between SEIR model and Regression model was found to be 2.01. Also, the value of R0 which is the spread of the disease was calculated to be 2.02. Expected cases may rise between 5000-6000 in the next two weeks of time. This study will help the Government and doctors in preparing their plans for the next two weeks. Based on the predictions for short-term interval, these models can be tuned for forecasting in long-term intervals.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

State preparation and measurement in a quantum simulation of the O(3) sigma model

Recently, Singh and Chandrasekharan showed that fixed points of the non-linear O(3) sigma model can be reproduced near a quantum phase transition of a spin model with just two qubits per lattice site. In a paper by the NuQS collaboration, the proposal is made to simulate such field theories on a quantum computer using the universal properties of a similar model. In this paper, following that direction, we demonstrate how to prepare the ground state of the model from and measure a dynamical quantity of interest, the O(3) Noether charge, on a quantum computer. In particular, we apply Trotter methods to obtain results for the complexity of adiabatic ground state preparation in both the weak-coupling and quantum-critical regimes and use shadow tomography to measure the dynamics of local observables. We then present and analyze a quantum algorithm based on non-unitary randomized simulation methods that may yield an approach suitable for intermediate-term noisy quantum devices.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Axial Vector Form Factors from Lattice QCD that Satisfy the PCAC Relation

Previous lattice QCD calculations of axial vector and pseudoscalar form factors show significant deviation from the partially conserved axial current (PCAC) relation between them. Since the original correlation functions satisfy PCAC, the observed deviations from the operator identity cast doubt on whether all the systematics in the extraction of form factors from the correlation functions are under control. We identify the problematic systematic as a missed excited state, whose energy as a function of the momentum transfer squared, $Q^2$, is determined from the analysis of the 3-point functions themselves. Its mass is much smaller than those of the excited states previously considered and including it impacts the extraction of all the ground state matrix elements. The form factors extracted using these mass/energy gaps satisfy PCAC and other consistency conditions, and validate the pion-pole dominance hypothesis. We also show that the extraction of the axial charge $g_A$ is very sensitive to the value of the mass gaps of the excited states used and current lattice data do not provide an unambiguous determination of these, unlike the $Q^2 \neq 0$ case. To highlight the differences and improvement between the conventional versus the new analysis strategy, we present a comparison of results obtained on a physical pion mass ensemble at $a\approx 0.0871\,\mathrm{fm}$. With the new strategy, we find $g_A = 1.30(6)$. A very significant improvement over previous lattice results is found for the axial charge radius $r_A = 0.74(6)\,\mathrm{fm}$, extracted using the $z$-expansion to parameterize the $Q^2$ behavior of $G_A(Q^2)$, and $g_P^\ast = 8.06(44)$ obtained using the pion pole-dominance ansatz to fit the $Q^2$ behavior of the induced pseudoscalar form factor $\widetilde{G}_P(Q^2)$.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Lattice QCD and Neutrino-Nucleus Scattering

This document is one of a series of whitepapers from the USQCD collaboration. Here, we discuss opportunities for lattice QCD in neutrino-oscillation physics, which inevitably entails nucleon and nuclear structure. In addition to discussing pertinent lattice-QCD calculations of nucleon and nuclear matrix elements, the interplay with models of nuclei is discussed. This program of lattice- QCD calculations is relevant to current and upcoming neutrino experiments, becoming increasingly important on the timescale of LBNF/DUNE and HyperK.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Nucleon Electromagnetic Form Factors in the Continuum Limit from 2+1+1-flavor Lattice QCD

Results are presented for the nucleon isovector electromagnetic form factors using 11 ensembles generated by the MILC collaboration using the 2+1+1-flavors HISQ action. They span 4 lattice spacings $a \sim$ 0.06, 0.09, 0.12 and 0.15~fm and 3 values of $M_π\sim 135, 225$ and 315 MeV. High-statistics estimates are used to perform a simultaneous extrapolation in the lattice spacing, lattice volume and light-quark masses. The $Q^2$ dependence over the range 0.05-1.4 ${\rm GeV}^2$ is investigated using both the $z$-expansion and the dipole form. Final $z$-expansion estimates for the isovector r.m.s. radius are $r_E = 0.769(27)(30)$ fm $r_M = 0.671(48)(76)$ fm and $μ^{p-n} = 3.939(86)(138)$ Bohr magneton. The first error is the combined uncertainty from the leading-order analysis, and the second is an estimate of the additional uncertainty due to using the leading order chiral-continuum-finite-volume fits. The dipole estimates, $r_E = 0.765(11)(8)$ fm, $r_M = 0.704(21)(29)$ fm and $μ^{p-n} = 3.975(84)(125)$, are consistent with those from the $z$-expansion but with smaller errors. Our analysis highlights three points. First, all data from the eleven ensembles and existing lattice data on, or close to, physical mass ensembles from other collaborations collapses more clearly onto a single curve when plotted versus $Q^2/M_N^2$ as compared to $Q^2$ with the scale set by quantities other than $M_N$. The difference between these two analyses is indicative of discretization errors, some of which presumably cancel when the data are plotted versus $Q^2/M_N^2$. Second, the size of the remaining deviation of this common curve from the Kelly curve is small and can be accounted for by statistical and possible systematic uncertainties. Third, to improve lattice estimates, high statistics data for $Q^2 < 0.1$ ${\rm GeV}^2$ are needed.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Quantum Computing for Neutrino-nucleus Scattering

Neutrino-nucleus cross section uncertainties are expected to be a dominant systematic in future accelerator neutrino experiments. The cross sections are determined by the linear response of the nucleus to the weak interactions of the neutrino, and are dominated by energy and distance scales of the order of the separation between nucleons in the nucleus. These response functions are potentially an important early physics application of quantum computers. Here we present an analysis of the resources required and their expected scaling for scattering cross section calculations. We also examine simple small-scale neutrino-nucleus models on modern quantum hardware. In this paper, we use variational methods to obtain the ground state of a three nucleon system (the triton) and then implement the relevant time evolution. In order to tame the errors in present-day NISQ devices, we explore the use of different error-mitigation techniques to increase the fidelity of the calculations.

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

Iso-vector and Iso-scalar Tensor Charges of the Nucleon from Lattice QCD

We present results for the iso-vector and flavor diagonal tensor charges $g^{u-d}_T$, $g^{u}_T$, $g^{d}_T$, and $g^{s}_T$ needed to probe novel tensor interactions at the TeV scale in neutron and nuclear $β$-decays and the contribution of the quark electric dipole moment (EDM) to the neutron EDM. The lattice QCD calculations were done using nine ensembles of gauge configurations generated by the MILC collaboration using the HISQ action with 2+1+1 dynamical flavors. These ensembles span three lattice spacings $a \approx 0.06, 0.09$ and $0.12 $ fm and three quark masses corresponding to the pion masses $M_π\approx 130, 220$ and $310 $ MeV. Using estimates from these ensembles, we quantify all systematic uncertainties and perform a simultaneous extrapolation in the lattice spacing, volume and light quark masses for the connected contributions. The final estimates of the connected nucleon (proton) tensor charge for the iso-vector combination is $g_T^{u-d} = 1.020(76) $ in the $\bar{\text{MS}}$ scheme at $2$ GeV. The additional disconnected quark loop contributions needed for the flavor-diagonal matrix elements are calculated using a stochastic estimator employing the truncated solver method with the all-mode-averaging technique. We find that the size of the disconnected contribution is smaller than the statistical error in the connected contribution. This allows us to bound the disconnected contribution and include it as an additional uncertainty in the flavor-diagonal charges. After a continuum extrapolation, we find $g_T^{u} = 0.774(66) $, $g_T^{d} = -0.233(28) $ and $g_T^{u+d} = 0.541(67) $. The strangeness tensor charge, that can make a significant contribution to the neutron EDM due to the large ratio $m_s/m_{u,d}$, is $g_T^{s}=0.008(9)$ in the continuum limit.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Neutron Electric Dipole Moment and Tensor Charges from Lattice QCD

We present Lattice QCD results on the neutron tensor charges including, for the first time, a simultaneous extrapolation in the lattice spacing, volume, and light quark masses to the physical point in the continuum limit. We find that the "disconnected" contribution is smaller than the statistical error in the "connected" contribution. Our estimates in the $\bar{\text{MS}}$ scheme at $2$ GeV, including all systematics, are $g_T^{d-u}=1.020(76)$, $g_T^d = 0.774(66)$, $g_T^u = - 0.233(28)$, and $g_T^s = 0.008(9)$. The flavor diagonal charges determine the size of the neutron electric dipole moment (EDM) induced by quark EDMs that are generated in many new scenarios of CP-violation beyond the Standard Model (BSM). We use our results to derive model-independent bounds on the EDMs of light quarks and update the EDM phenomenology in split Supersymmetry with gaugino mass unification, finding a stringent upper bound of $d_n < 4 \times 10^{-28} \, e$ cm for the neutron EDM in this scenario.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Neutron Electric Dipole Moment from quark Chromoelectric Dipole Moment

The connection between a regularization-independent symmetric momentum substraction (RI-$\tilde{\rm S}$MOM) and the $\overline{\rm MS}$ scheme for the quark chromo EDM operators is discussed. A method for evaluating the neutron EDM from quark chromoEDM is described. A preliminary study of the signal in the matrix element using clover quarks on a highly improved staggered quark (HISQ) ensemble is shown.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Nucleon Charges, Form-factors and Neutron EDM

We present an update of our analysis of statistical and systematic errors in the calculation of iso-vector scalar, axial and tensor charges of the nucleon. The calculations are done using $N_f=2+1+1$ flavor HISQ ensembles generated by the MILC Collaboration at three values of the lattice spacing ($a=0.12,\ 0.09,$ and $0.06$ fm) and three values of the quark mass ($M_π\approx 310,\ 220$ and $130$ MeV); and clover fermions for calculating the correlation functions, i.e., we use a clover-on-HISQ lattice formulation. The all-mode-averaging method allows us to increase the statistics by a factor of eight for the same computational cost leading to a better understanding of and control over excited state contamination. Our current results, after extrapolation to the continuum limit and physical pion mass are $g_A^{u-d} = 1.21(3)$, $g_T^{u-d} = 1.005(59)$ and $g_S^{u-d} = 0.95(12) $. Further checks of control over all systematic errors, especially in $g_A^{u-d}$, are still being performed. Using results for the flavor-diagonal charges, $g_T^{u} $, $g_T^{d} $ and $g_T^{s}$, we analyze contributions of the quark electric dipole moment to the neutron EDM and the consequences for split SUSY model.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Quark Chromoelectric Dipole Moment Contribution to the Neutron Electric Dipole Moment

The quark chromo-electric dipole moment operator and the pseudo-scalar fermion bilinear with which it mixes under renormalization can both be included in a calculation of the electromagnetic form factor of the nucleon using the Schwinger source method. A preliminary calculation of these operators using clover quarks on HISQ lattices generated by MILC collaboration will be presented showing the quality of the signal in the correlators necessary for calculating the neutron electric dipole moment.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Tuning the hopping parameter in the Oktay-Kronfeld action for charm and bottom quarks on a MILC HISQ ensemble

The first step in the calculation of semi-leptonic form factors in the decay of heavy mesons is the tuning of the hopping parameter $κ$ for the charm and bottom quark masses. Results for the Oktay-Kronfeld (OK) action are presented for one $N_f=2+1+1$ HISQ ensemble generated by the MILC collaboration at $a\approx 0.12\,\mathrm{fm}$ and $M_π\approx 310$ MeV. Estimates of hyperfine splitting of heavy-light and heavy-heavy mesons are presented and the inconsistency parameter is evaluated.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Disconnected Quark Loop Contributions to Nucleon Structure

We calculate the disconnected contribution to isoscalar nucleon charges for scalar, axial and tensor channels of light and strange quarks. The calculation has been done with the Clover valence quarks on the MILC $N_f=2+1+1$ HISQ lattices whose light quark masses corresponding to the pion masses of 305 MeV and 217 MeV at $a \approx 0.12$ fm and 312 MeV at $a \approx 0.09$ fm. All-mode-averaging technique is used for the evaluation two-point functions. Disconnected quark loops are estimated by using the truncated solver method with Gaussian random noise sources. Contamination from the excited states is removed by fitting the results of various source-sink separations and operator insertions to the formula including up to the first excited state, simultaneously.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Precision calculations of nucleon charges $g_A$, $g_S$, $g_T$

We present a detailed analysis of statistical and systematic errors in the calculation of matrix elements of iso-vector scalar, axial and tensor charges between a neutron and a proton state. These analyses are being done on dynamical $N_f=2+1+1$ HISQ configurations generated by the MILC Collaboration using valence clover fermions. Using ensembles at three values of the lattice spacing ($a=0.12,\ 0.09,$ and $0.06$ fm) and three values of the quark mass ($M_π\approx 310,\ 220$ and $130$ MeV) we find that the estimates of the tensor charge are stable and it can be extracted with $5\%$ precision with O(10,000) measurements. We also find that higher statistics are needed to resolve the various uncertainties in the calculation of $g_A$ and improve the signal in $g_S$, which with present data has large errors. A brief status report on the mixing and renormalization of novel operators contributing to nEDM is also given.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Neutron Electric Dipole Moments from Beyond the Standard Model Physics

Neutron Electric Dipole Moment (nEDM), a generic feature of CP-violation, is predicted to be very small in the Standard Model, but can be much larger in most extensions of the model. In this talk, I will discuss the classification of the CP violating operators up to dimension 6 that can give rise to nEDM, and then describe the mixing and renormalization structure of the operators of dimension 5 and lower in both dimensional and cutoff regularizations in general terms. Finally I will describe how to connect the dimension 5 operators, in particular, the Chromoelectric Dipole Moment of the quarks, between MSbar scheme and a Regularization Independent prescription in the chiral limit.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Nucleon Charges and Electromagnetic Form Factors from 2+1+1-Flavor Lattice QCD

We present lattice-QCD results on the nucleon isovector axial, scalar and tensor charges, the isovector electromagnetic Dirac and Pauli form factors, and the connected parts of the isoscalar charges. The calculations have been done using two ensembles of HISQ lattices generated by the MILC Collaboration with 2+1+1 dynamical flavors at a lattice spacing of 0.12 fm and with light-quark masses corresponding to pions with masses 310 and 220 MeV. We perform a systematic study including excited-state degrees of freedom and examine the dependence of the extracted nucleon matrix elements on source-sink separation. This study demonstrates with high-statistics data that including excited-state contributions and generating data at multiple separations is necessary to remove contamination that would otherwise lead to systematic error. We also determine the renormalization constants of the associated quark bilinear operators in the RI-sMOM scheme and make comparisons of our renormalized results with previous dynamical-lattice calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Probing TeV scale physics in precision UCN decays

We present the calculation of matrix elements of iso-vector scalar, axial and tensor charges between a neutron and a proton state on dynamical $N_f=2+1+1$ HISQ configurations generated by the MILC Collaboration using valence clover fermions. These matrix elements are needed to probe novel scalar and tensor interactions in neutron beta-decay that can arise in extensions to the Standard Model at the TeV scale. Results are presented at one value of the lattice spacing, $a=0.12$ fm, and two values of light quarks corresponding to $M_π=310$ and $220$ MeV. We discuss two sources of systematic errors, contribution of excited states to these matrix elements and the renormalization constants, and the efficacy of methods used to control them.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The equation of state in (2+1)-flavor QCD

We present results for the equation of state in (2+1)-flavor QCD using the highly improved staggered quark action and lattices with temporal extent $N_τ=6,~8,~10$, and $12$. We show that these data can be reliably extrapolated to the continuum limit and obtain a number of thermodynamic quantities and the speed of sound in the temperature range $(130-400)$ MeV. We compare our results with previous calculations, and provide an analytic parameterization of the pressure, from which other thermodynamic quantities can be calculated, for use in phenomenology. We show that the energy density in the crossover region, $145~ {\rm MeV} \leq T \leq 163$ MeV, defined by the chiral transition, is $ε_c=(0.18-0.5)~{\rm GeV}/{\rm fm}^3$, $i.e.$, $(1.2-3.1)\ ε_{\rm nuclear}$. At high temperatures, we compare our results with resummed and dimensionally reduced perturbation theory calculations. As a byproduct of our analyses, we obtain the values of the scale parameters $r_0$ from the static quark potential and $w_0$ from the gradient flow.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The QCD phase transition with physical-mass, chiral quarks

We report on the first lattice calculation of the QCD phase transition using chiral fermions at physical values of the quark masses. This calculation uses 2+1 quark flavors, spatial volumes between (4 fm$)^3$ and (11 fm$)^3$ and temperatures between 139 and 196 MeV . Each temperature was calculated using a single lattice spacing corresponding to a temporal Euclidean extent of $N_t=8$. The disconnected chiral susceptibility, $χ_{\rm disc}$ shows a pronounced peak whose position and height depend sensitively on the quark mass. We find no metastability in the region of the peak and a peak height which does not change when a 5 fm spatial extent is increased to 10 fm. Each result is strong evidence that the QCD ``phase transition'' is not first order but a continuous cross-over for $m_π=135$ MeV. The peak location determines a pseudo-critical temperature $T_c = 155(1)(8)$ MeV. Chiral $SU(2)_L\times SU(2)_R$ symmetry is fully restored above 164 MeV, but anomalous $U(1)_A$ symmetry breaking is non-zero above $T_c$ and vanishes as $T$ is increased to 196 MeV.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Fluctuations and Correlations of net baryon number, electric charge, and strangeness: A comparison of lattice QCD results with the hadron resonance gas model

We calculate the quadratic fluctuations of net baryon number, electric charge and strangeness as well as correlations among these conserved charges in (2+1)-flavor lattice QCD at zero chemical potential. Results are obtained using calculations with tree level improved gauge and the highly improved staggered quark (HISQ) actions with almost physical light and strange quark masses at three different values of the lattice cut-off. Our choice of parameters corresponds to a value of 160 MeV for the lightest pseudo scalar Goldstone mass and a physical value of the kaon mass. The three diagonal charge susceptibilities and the correlations among conserved charges have been extrapolated to the continuum limit in the temperature interval 150 MeV <T < 250 MeV. We compare our results with the hadron resonance gas (HRG) model calculations and find agreement with HRG model results only for temperatures T<= 150 MeV. We observe significant deviations in the temperature range 160 MeV < T < 170 MeV and qualitative differences in the behavior of the three conserved charge sectors. At T < 160 MeV quadratic net baryon number fluctuations in QCD agree with HRG model calculations while, the net electric charge fluctuations in QCD are about 10% smaller and net strangeness fluctuations are about 20% larger. These findings are relevant to the discussion of freeze-out conditions in relativistic heavy ion collisions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Probing novel TeV physics through precision calculations of scalar and tensor charges of the nucleon

We present an update on the calculation of matrix elements of iso-vector scalar, axial and tensor charges between a neutron and a proton state. These matrix elements are needed to probe novel scalar and tensor interactions in neutron beta-decay that can arise in extensions of the Standard Model at the TeV scale. Our calculations are being done using valence clover fermions on dynamical N_f=2+1+1 HISQ configurations generated by the MILC Collaboration. We provide preliminary estimates of the dependence of these matrix elements on the light quark masses, lattice spacing, and the time separation between the source and sink of the nucleons. We also find that the renormalization constants calculated using the RI-sMOM scheme are close to unity for the HYP smeared HISQ lattices.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Probing TeV scale physics via ultra cold neutron decays and calculating non-standard baryon matrix elements

We motivate undertaking precision analyses of neutron decays to look for signatures of new scalar and tensor interactions that can arise in extensions of the Standard Model at the TeV scale. The key ingrediant needed to connect experimental data with theoretical analysis are high-precision calculations of matrix elements of isovector bilinear operators between the decaying neutron and final state proton. We describe the status of our Lattice QCD program of using valence clover fermions on dynamical N_f=2+1+1 HISQ configurations generated by the MILC Collaboration. On the theoretical side we use the effective field theory method and provide both model independent and dependent analyses to obtain bounds on possible scalar and tensor interactions, both from low energy experiments and LHC data.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The chiral transition and U(1)_A symmetry restoration from lattice QCD using Domain Wall Fermions

We present results on both the restoration of the spontaneously broken chiral symmetry and the effective restoration of the anomalously broken U(1)_A symmetry in finite temperature QCD at zero chemical potential using lattice QCD. We employ domain wall fermions on lattices with fixed temporal extent N_τ= 8 and spatial extent N_σ= 16 in a temperature range of T = 139 - 195 MeV, corresponding to lattice spacings of a \approx 0.12 - 0.18 fm. In these calculations, we include two degenerate light quarks and a strange quark at fixed pion mass m_π= 200 MeV. The strange quark mass is set near its physical value. We also present results from a second set of finite temperature gauge configurations at the same volume and temporal extent with slightly heavier pion mass. To study chiral symmetry restoration, we calculate the chiral condensate, the disconnected chiral susceptibility, and susceptibilities in several meson channels of different quantum numbers. To study U(1)_A restoration, we calculate spatial correlators in the scalar and pseudo-scalar channels, as well as the corresponding susceptibilities. Furthermore, we also show results for the eigenvalue spectrum of the Dirac operator as a function of temperature, which can be connected to both U(1)_A and chiral symmetry restoration via Banks-Casher relations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Theoretical Bounds on New Four-Fermion Interactions and TeV Scale Physics

The standard model weak interactions can be described by four-fermion V-A operators at low energies. New physics at the TeV scale can, however, generate the other Lorentz structures. In this talk, we review the constraints on such interactions from nuclear and hadronic decays, as well as from collider searches. Currently the most stringent bounds come from the analysis of the 0+ to 0+ nuclear and the pi to e nu gamma radiative pion decays. In the near future, the ultracold neutron beta decay experiments and the direct LHC measurements will compete in setting the most stringent bounds, provided, however, that the neutron-to-proton non-perturbative transition matrix elements can be calculated to a level of 10-20% accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Equation of State from Lattice QCD Calculations

We provide a status report on the calculation of the Equation of State (EoS) of QCD at finite temperature using lattice QCD. Most of the discussion will focus on comparison of recent results obtained by the HotQCD and Wuppertal-Budapest (W-B) collaborations. We will show that very significant progress has been made towards obtaining high precision results over the temperature range of T=150-700 MeV. The various sources of systematic uncertainties will be discussed and the differences between the two calculations highlighted. Our final conclusion is that the lattice results of EoS are getting precise enough to justify being used in the phenomenological analysis of heavy ion experiments at RHIC and LHC.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Probing Novel Scalar and Tensor Interactions from (Ultra)Cold Neutrons to the LHC

Scalar and tensor interactions were once competitors to the now well-established V-A structure of the Standard Model weak interactions. We revisit these interactions and survey constraints from low-energy probes (neutron, nuclear, and pion decays) as well as collider searches. Currently, the most stringent limit on scalar and tensor interactions arise from 0+ -> 0+ nuclear decays and the radiative pion decay pi -> e nu gamma, respectively. For the future, we find that upcoming neutron beta decay and LHC measurements will compete in setting the most stringent bounds. For neutron beta decay, we demonstrate the importance of lattice computations of the neutron-to-proton matrix elements to setting limits on these interactions, and provide the first lattice estimate of the scalar charge and a new average of existing results for the tensor charge. Data taken at the LHC is currently probing these interactions at the 10^-2 level (relative to the standard weak interactions), with the potential to reach the 10^-3 level. We show that, with some theoretical assumptions, the discovery of a charged spin-0 resonance decaying to an electron and missing energy implies a lower limit on the strength of scalar interactions probed at low energy.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

The QCD EoS from simulations on BlueGene L Supercomputers at LLNL and NYBlue

We present results for the QCD Equation of State (EoS) obtained using simulations of lattice QCD at zero chemical potential. Our high statistics results compare improved asqtad and p4fat3 staggered quarks on lattices with a temporal extent N_tau = 6 and 8 and light quark masses approximately one fifth and one tenth the strange quark mass. We find that the two actions give consistent results and estimate that the trace anomaly (epsilon - 3p)/T^4 obtained on N_tau = 8 lattices represents the continuum value to better than 20% uncertainty over the temperature range 140-700 MeV. The precision in the estimates of energy density and pressure are better, therefore, we conclude that lattice estimates of the energy density and pressure should be used in the phenomenological analysis of RHIC and LHC data. We also find a consistent picture of the crossover temperature from all observables studied, with the best estimated range to be 185-195 MeV. These calculations are being carried out on the IBM BlueGene/L supercomputer at Lawrence Livermore National Laboratory and at the New York Center for Computational Science (NYBlue).

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Testing improved staggered fermions with $m_s$ and $B_K$

We study the improvement of staggered fermions using hypercubically smeared (HYP) links. We calculate the strange quark mass and the kaon B-parameter, $B_K$, in quenched QCD on a $16^3 \times 64$ lattice at $β=6.0$. We find $m_s(\bar{\rm MS},2 {\rm GeV})=101.2\pm1.3\pm4 $MeV and $B_K(\bar{\rm MS},2 {\rm GeV}) = 0.578 \pm 0.018\pm 0.042$, where the first error is from statistics and fitting, and the second from using one-loop matching factors. The scale ($1/a=1.95$GeV) is set by $M_ρ$, and $m_s$ is determined using the kaon mass. Comparing to quenched results obtained using unimproved staggered fermions and other discretizations, we argue that the size of discretization errors in $B_K$ is substantially reduced by improvement.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Scaling behavior of improvement and renormalization constants

This talk summarizes results for all the scale independent renormalization constants for bilinear currents ($Z_A$, $Z_V$, and $Z_S/Z_P$), the improvement constants ($c_A$, $c_V$, and $c_T$), the quark mass dependence of $Z_O$, and the coefficients of the equation of motion operators for $O(a)$ improved lattice QCD. Using data at $β=6.0$, 6.2 and 6.4 we study the scaling behavior of these quantities and quantify residual discretization errors.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Non-perturbative Renormalization Constants using Ward Identities

We extend the application of axial Ward identities to calculate $b_A, b_P$ and $b_T$, coefficients that give the mass dependence of the renormalization constants of the corresponding bilinear operators in the quenched theory. The extension relies on using operators with non-degenerate quark masses. It allows a complete determination of the O(a) improvement coefficients for bilinears in the quenched approximation using Ward Identities alone. Only the scale dependent normalization constants $Z_P^0$ (or $Z_S^0$) and $Z_T$ are undetermined. We present results of a pilot numerical study using hadronic correlators.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

Critical Exponents of the 3-D Ising Model

We present a status report on the ongoing analysis of the 3D Ising model with nearest-neighbor interactions using the Monte Carlo Renormalization Group (MCRG) and finite size scaling (FSS) methods on $64^3$, $128^3$, and $256^3$ simple cubic lattices. Our MCRG estimates are $K_{nn}^c=0.221655(1)(1)$ and $ν=0.625(1)$. The FSS results for $K^c$ are consistent with those from MCRG but the value of $ν$ is not. Our best estimate $η= 0.025(6)$ covers the spread in the MCRG and FSS values. A surprise of our calculation is the estimate $ω\approx 0.7$ for the correction-to-scaling exponent. We also present results for the renormalized coupling $g_R$ along the MCRG flow and argue that the data support the validity of hyperscaling for the 3D Ising model.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

Prospects of Solving Grand Challenge Problems

The recent woes of the supercomputer industry and changes in federal funding have caused some scientists to re-evaluate the means by which they hope to solve Grand Challenge problems. I evaluate the potential of Massively Parallel Processors (MPP) within this context and the state of today's MPP. I stress that for solving large-scale problems MPP are crucial and that it is essential to seek a balance between CPU performance, memory access time, inter-node communications, and I/O. To achieve this it is important to preserve certain characteristics of the hardware while selecting the hottest processor to design the machine around. I emphasize that for long term stability and growth of parallel computing priority should be given to standardizing software so that the same code can run on different platforms and on machines ranging from clusters of workstations to MPP.

0 citations0 reviewsNo code linkedOpen access
preprint1994arXiv

A pot-pourri of results in QCD from large lattice simulations on the CM5

We present a status report on simulations being done on $32^3 \times 64$ lattices at $β= 6.0$ using quenched Wilson fermions. Results for the spectrum, decay constants, the kaon B-parameter $B_K$, and semi-leptonic form factors are given based on the current statistical sample of 28 configurations. These ``Grand Challenge'' calculations are being done on the CM5 at Los Alamos in collaboration with G. Kilcup, S. Sharpe and P. Tamayo. We end with a brief statement of code performance.

0 citations0 reviewsNo code linkedOpen access