Search

We introduce a class of efficient multiple right-hand side multigrid algorithm for domain wall fermions. The simultaneous solution for a modest number of right hand sides concurrently allows for a significant reduction in the time spent solving the coarse grid operator in a multigrid preconditioner. We introduce a preconditioned block conjuate gradient with a multigrid preconditioner, giving additional algorithmic benefit from the multiple right hand sides. There is also a very significant additional to computation rate benefit to multiple right hand sides. This both increases the arithmetic intensity in the coarse space and increases the amount of work being performed in each subroutine call, leading to excellent performance on modern GPU architectures. Further, the software implementation makes use of vendor linear algebra routines (batched GEMM) that can make use of high throughput tensor hardware on recent Nvidia, AMD and Intel GPUs. The cost of the coarse space is made sub-dominant in this algorithm, and benchmarks from the Frontier supercomputer system show up to a factor of twenty speed up over the standard red-black preconditioned conjugate gradient algorithm on a large system with physical quark masses., Comment: 33 pages

We present our study of the $\rho(770)$ and $K^*(892)$ resonances from lattice quantum chromodynamics (QCD) employing domain-wall fermions at physical quark masses. We determine the finite-volume energy spectrum in various momentum frames and obtain phase-shift parameterizations via the L\"uscher formalism, and as a final step the complex resonance poles of the $\pi \pi$ and $K \pi$ elastic scattering amplitudes via an analytical continuation of the models. By sampling a large number of representative sets of underlying energy-level fits, we also assign a systematic uncertainty to our final results. This is a significant extension to data-driven analysis methods that have been used in lattice QCD to date, due to the two-step nature of the formalism. Our final pole positions, $M+i\Gamma/2$, with all statistical and systematic errors exposed, are $M_{K^{*}} = 893(2)(8)(54)(2)~\mathrm{MeV}$ and $\Gamma_{K^{*}} = 51(2)(11)(3)(0)~\mathrm{MeV}$ for the $K^*(892)$ resonance and $M_{\rho} = 796(5)(15)(48)(2)~\mathrm{MeV}$ and $\Gamma_{\rho} = 192(10)(28)(12)(0)~\mathrm{MeV}$ for the $\rho(770)$ resonance. The four differently grouped sources of uncertainties are, in the order of occurrence: statistical, data-driven systematic, an estimation of systematic effects beyond our computation (dominated by the fact that we employ a single lattice spacing), and the error from the scale-setting uncertainty on our ensemble.

We present the first ab initio calculation at physical quark masses of scattering amplitudes describing the lightest pseudoscalar mesons interacting via the strong force in the vector channel. Using lattice quantum chromodynamics, we postdict the defining parameters for two short-lived resonances, the $\rho(770)$ and $K^*(892)$, which manifest as complex energy poles in $\pi \pi$ and $K \pi$ scattering amplitudes, respectively. The calculation proceeds by first computing the finite-volume energy spectrum of the two-hadron systems, and then determining the amplitudes from the energies using the L\"uscher formalism. The error budget includes a data-driven systematic error, obtained by scanning possible fit ranges and fit models to extract the spectrum from Euclidean correlators, as well as the scattering amplitudes from the latter. The final results, obtained by analytically continuing multiple parameterizations into the complex energy plane, are $M_\rho = 796(5)(50)~\mathrm{MeV}$, $\Gamma_\rho = 192(10)(31)~\mathrm{MeV}$, $M_{K^*} = 893(2)(54)~\mathrm{MeV}$ and $\Gamma_{K^*} = 51(2)(11)~\mathrm{MeV}$, where the subscript indicates the resonance and $M$ and $\Gamma$ stand for the mass and width, respectively, and where the first bracket indicates the statistical and the second bracket the systematic uncertainty.

We present a direct lattice QCD calculation of the $x$-dependence of the pion distribution amplitude (DA), which is performed using the quasi-DA in large momentum effective theory on a domain-wall fermion ensemble at physical quark masses and spacing $a\approx 0.084$ fm. The bare quais-DA matrix elements are renormalized in the hybrid scheme and matched to $\overline{\rm MS}$ with a subtraction of the leading renormalon in the Wilson-line mass. For the first time, we include threshold resummation in the perturbative matching onto the light-cone DA, which resums the large logarithms in the soft gluon limit at next-to-next-to-leading log. The resummed results show controlled scale-variation uncertainty within the range of momentum fraction $x\in[0.25,0.75]$ at the largest pion momentum $P_z\approx 1.85$~GeV. In addition, we apply the same analysis to quasi-DAs from a highly-improved-staggered-quark ensemble at physical pion mass and $a=0.076$ fm. By comparison we find with $2\sigma$ confidence level that the DA obtained from chiral fermions is flatter and lower near $x=0.5$., Comment: update to match the version published in journal

We present non-perturbative results for beyond the standard model kaon mixing matrix elements in the isospin symmetric limit ($m_u=m_d$) of QCD, including a complete estimate of all dominant sources of systematic error. Our results are obtained from numerical simulations of lattice QCD with $N_f = 2+1$ flavours of dynamical domain wall fermions. For the first time, these quantities are simulated directly at the physical pion mass $m_\pi$~$\sim$~$139\,\mathrm{MeV}$ for two different lattice spacings. We include data at three lattice spacings in the range $a = 0.11 $ - $ 0.07\,\mathrm{fm}$ and with pion masses ranging from the physical value up to 450$\,\mathrm{MeV}$. Compared to our earlier work, we have added both direct calculations at physical quark masses and a third lattice spacing making the removal of discretisation effects significantly more precise and eliminating the need for any significant mass extrapolation beyond the range of simulated data. We renormalise the lattice operators non-perturbatively using RI-SMOM off-shell schemes. These schemes eliminate the need to model and subtract non-perturbative pion poles that arises in the RI-MOM scheme and, since the calculations are performed with domain wall fermions, the unphysical mixing between chirality sectors is suppressed. Our results for the bag parameters in the $\overline{\mathrm{MS}}$ scheme at $3\,\mathrm{GeV}$ are $B_K~\equiv~\mathcal{B}_1 = 0.5240(17)(54)$, $\mathcal{B}_2 = 0.4794(25)(35)$, $\mathcal{B}_3 = 0.746(13)(17)$, $\mathcal{B}_4 = 0.897(02)(10)$ and $\mathcal{B}_5 = 0.6882(78)(94)$, where the first error is from lattice uncertainties and the second is the uncertainty due to the perturbative matching to $\overline{\mathrm{MS}}$., Comment: 31 pages, 18 figures, 21 tables, 2 ancillary files

I review recent research and advances in algorithms for solvers and gauge generation, with an emphasis on practical algorithms for four dimensional simulations. Particular consideration is given to advances in multigrid solvers, fourier acceleration and field transformation approaches to accelerating evolution dynamics, and to parallel tempering approaches to solving the topological tunneling problem. Particular consideration is given to the interaction between rapidly evolving computer architecture and optimal algorithms that exploit these. In this conference, nascent machine learning algorithms were separately reviewed, Comment: Lattice 2023