6533b823fe1ef96bd127eb67

RESEARCH PRODUCT

Nucleon matrix elements from lattice QCD with all-mode-averaging and a domain-decomposed solver: An exploratory study

Eigo ShintaniHartmut WittigThomas RaeGeorg M. Von Hippel

subject

QuarkPhysicsNuclear and High Energy PhysicsParticle physics010308 nuclear & particles physicsHigh Energy Physics - Lattice (hep-lat)Scalar (physics)FOS: Physical sciencesCharge (physics)Lattice QCDComputational Physics (physics.comp-ph)Solver01 natural sciencesMatrix (mathematics)High Energy Physics - Lattice13. Climate actionQuantum electrodynamics0103 physical scienceslcsh:QC770-798ddc:530lcsh:Nuclear and particle physics. Atomic energy. RadioactivityTensor010306 general physicsNucleonPhysics - Computational Physics

description

We study the performance of all-mode-averaging (AMA) when used in conjunction with a locally deflated SAP-preconditioned solver, determining how to optimize the local block sizes and number of deflation fields in order to minimize the computational cost for a given level of overall statistical accuracy. We find that AMA enables a reduction of the statistical error on nucleon charges by a factor of around two at the same cost when compared to the standard method. As a demonstration, we compute the axial, scalar and tensor charges of the nucleon in $N_f=2$ lattice QCD with non-perturbatively O(a)-improved Wilson quarks, using O(10,000) measurements to pursue the signal out to source-sink separations of $t_s\sim 1.5$ fm. Our results suggest that the axial charge is suffering from a significant amount (5-10%) of excited-state contamination at source-sink separations of up to $t_s\sim 1.2$ fm, whereas the excited-state contamination in the scalar and tensor charges seems to be small.

10.1016/j.nuclphysb.2016.11.003http://www.sciencedirect.com/science/article/pii/S0550321316303534