Search

The generation of ensembles of gauge configurations is a considerable expense. The preservation and curation of these ensembles constitutes a valuable shared resource for the lattice field theory community. The organizers of Lattice 2022 dedicated a parallel session to the presentation of gauge ensembles and their generation, plans for ensemble publication and data management/storage activities of different collaborations. A summary of the twelve contributions is presented here., Comment: Contribution to the 39th International Symposium on Lattice Field Theory (LATTICE2022), 8-13 August, 2022, Bonn, Germany, 23 pages

We show that using the multi-splitting algorithm as a preconditioner for the domain wall Dirac linear operator, arising in lattice QCD, effectively reduces the inter-node communication cost, at the expense of performing more on-node floating point and memory operations. Correctly including the boundary \textit{snake} terms, the preconditioner is implemented in the QUDA framework, where it is found that utilizing kernel fusion and the tensor cores on NVIDIA GPUs is necessary to achieve a sufficiently performant preconditioner. A reduced-dimension (reduced-$L_s$) strategy is also proposed and tested for the preconditioner. We find the method achieves lower time to solution than regular CG at high node count despite the additional local computational requirements from the preconditioner. This method could be useful for supercomputers with more on-node flops and memory bandwidth than inter-node communication bandwidth., Comment: Add DOI

We present a lattice QCD calculation of the $\Delta I=1/2$, $K\to\pi\pi$ decay amplitude $A_0$ and $\varepsilon'$, the measure of direct CP-violation in $K\to\pi\pi$ decay, improving our 2015 calculation of these quantities. Both calculations were performed with physical kinematics on a $32^3\times 64$ lattice with an inverse lattice spacing of $a^{-1}=1.3784(68)$ GeV. However, the current calculation includes nearly four times the statistics and numerous technical improvements allowing us to more reliably isolate the $\pi\pi$ ground-state and more accurately relate the lattice operators to those defined in the Standard Model. We find ${\rm Re}(A_0)=2.99(0.32)(0.59)\times 10^{-7}$ GeV and ${\rm Im}(A_0)=-6.98(0.62)(1.44)\times 10^{-11}$ GeV, where the errors are statistical and systematic, respectively. The former agrees well with the experimental result ${\rm Re}(A_0)=3.3201(18)\times 10^{-7}$ GeV. These results for $A_0$ can be combined with our earlier lattice calculation of $A_2$ to obtain ${\rm Re}(\varepsilon'/\varepsilon)=21.7(2.6)(6.2)(5.0) \times 10^{-4}$, where the third error represents omitted isospin breaking effects, and Re$(A_0)$/Re$(A_2) = 19.9(2.3)(4.4)$. The first agrees well with the experimental result of ${\rm Re}(\varepsilon'/\varepsilon)=16.6(2.3)\times 10^{-4}$. A comparison of the second with the observed ratio Re$(A_0)/$Re$(A_2) = 22.45(6)$, demonstrates the Standard Model origin of this "$\Delta I = 1/2$ rule" enhancement., Comment: Updated to published version. 95 pages, 12 figures

The RBC and UKQCD Collaborations continue to produce 2+1 flavor domain wall fermion ensembles, currently focusing on an ensemble with a $96^3 \times 192$ volume on SUMMIT at ORNL with $1/a \approx 2.8$ GeV, and smaller ensembles at stronger couplings. The $1/a \approx 2.8$ GeV ensemble uses the Exact One Flavor Algorithm for the strange quark, along with the Multisplitting Preconditioned Conjugate Gradient for solving the Dirac equation. We report on our progress and experience to date with the evolution of this ensemble., Comment: 7 pages

We show that using the multisplitting algorithm as a preconditioner for conjugate gradient inversion of the domain wall Dirac operator could effectively reduce the inter- node communication cost, at the expense of performing more on-node floating point operations. This method could be useful for supercomputers with far more on-node flops than inter-node communication bandwidth.

We present the first calculation of the kaon semileptonic form factor with sea and valence quark masses tuned to their physical values in the continuum limit of 2+1 flavour domain wall lattice QCD. We analyse a comprehensive set of simulations at the phenomenologically convenient point of zero momentum transfer in large physical volumes and for two different values of the lattice spacing. Our prediction for the form factor is f+(0)=0.9685(34)(14) where the first error is statistical and the second error systematic. This result can be combined with experimental measurements of K->pi decays for a determination of the CKM-matrix element for which we predict |Vus|=0.2233(5)(9) where the first error is from experiment and the second error from the lattice computation., Comment: 21 pages, 7 figures, 6 tables

We study the autocorrelations of observables constructed from the topological charge density, such as the topological charge on a time slice or in a subvolume, using a series of hybrid Monte Carlo simulations of pure SU(3) gauge theory with both periodic and open boundary conditions. We show that the autocorrelation functions of these observables obey a simple diffusion equation and we measure the diffusion coefficient, finding that it scales like the square of the lattice spacing. We use this result and measurements of the rate of tunneling between topological charge sectors to calculate the scaling behavior of the autocorrelation times of these observables on periodic and open lattices. There is a characteristic lattice spacing at which open boundary conditions become worthwhile for reducing autocorrelations and we show how this lattice spacing is related to the diffusion coefficient, the tunneling rate, and the lattice Euclidean time extent., Comment: 36 pages, 14 figures

We present preliminary results for the calculation of the Kaon Bag parameter $B_K$ in $N_f=2+1$ lattice QCD, using M\"obius Domain Wall Fermion ensembles generated by the RBC-UKQCD collaboration. This computation is done directly at physical meson masses, so that we do not have to rely on chiral perturbation theory or any other mass extrapolation. In parallel, the four-quark operator is renormalised through the Rome-Southampton technique. Finally, we compare our value with previous results and draw some conclusions about the remaining dominant contributions in our error budget., Comment: presented at the 31st International Symposium on Lattice Field Theory (Lattice 2013), 29 July - 3 August 2013, Mainz, Germany

We present results from a 2+1 flavor DWF calculation at 1/a = 3 GeV and discuss strategies for similar calculations at finer lattice spacings which will target charm physics. At weak coupling the autocorrelation time of the global topological charge becomes very long because the HMC algorithm has trouble moving between topological sectors. We report the results of simulations that test two ideas for reducing the autocorrelation time of topological charge. In weak coupling quenched simulations we find that the open boundary conditions suggested by L\"uscher and Schaefer do not prevent the appearance of extremely long autocorrelation times for topological observables. We discuss the idea of a "dislocation-enhancing determinant" and show that it can produce an increase in topological tunneling., Comment: 7 pages, 4 figures. Presented at the 31st International Symposium on Lattice Field Theory (Lattice 2013), 29 July - 3 August 2013, Mainz, Germany

With a substantial amount of simulations, we have explored the system across a wide range of lattice scales. We have located a lattice artifact, first order bulk transition, have studied its properties, and found that the flavor-singlet scalar meson mass vanishes at the critical endpoint. We will discuss the lattice phase diagrams and the continuum limits for both a spontaneous chiral symmetry breaking phase and an infrared conformal phase, and compare results with other groups., Comment: 7 pages, 7 figures; Contribution to SCGT12 "KMI-GCOE Workshop on Strong Coupling Gauge Theories in the LHC Perspective", 4-7 Dec. 2012, Nagoya University

We report on new data from additional zero temperature simulations of QCD with 12 flavors. This is a continuation of previous studies using the DBW2 gauge action and naive staggered fermions. With the use of the force gradient integrator and a multiple-quark-mass preconditioned HMC, we have done simulations with input quark masses from $m_q=0.003$ to $m_q=0.008$. We have observed a metastable, first order, bulk transition that occurs at small input quark masses. As the quark mass increases, this first order bulk transition ends at a second order critical point, and, for still heavier quark masses, becomes the cross-over we have previously reported. We present measurements of hadron masses, decay constants and other low energy observables in the small quark mass region on the weak coupling side of the bulk transition. Our results show that the behavior of the system is still consistent with spontaneously broken chiral symmetry. We also discuss a preliminary investigation into the behavior of the bulk transition itself. We have found that, as the system approaches the second order critical end point, the scalar singlet meson becomes lighter. Thus it appears that the critical endpoint corresponds to a continuum limit theory only involving scalars and, following known triviality arguments, this is likely a free field theory. The presence of this critical endpoint could influence scaling of lattice observables in the conventional continuum limit., Comment: 7 pages, 13 figures, 2 tables. Contribution to the XXIX International Symposium on Lattice Field Theory, July 10-16, 2011, Squaw Valley, Lake Tahoe, California

We have implemented a variant of the force gradient integrator proposed by Kennedy et.al. and are using it in our production 2+1 flavor DWF simulations with pion masses of 180 MeV in (4.5fm)3 volumes. We find modest speed-ups (\sim 20%) from using the force gradient integrator, compared to our previously used Omelyan integrator. On other ensembles, primarily finite temperature 2+1 flavor DWF QCD, we have extensively tuned the Hasenbusch preconditioning masses and achieved speed-ups of 2-3x. Here we have also switched to the force gradient integrator, but this change has not had any impact on the speed. We also report on an improved solver for DWF, which uses M\"obius fermions, with a smaller fifth dimension than the original DWF fermions, as an intermediate step in the generation of solutions of the Dirac equation. This approach cuts the number of effective Dirac applications by approximately a factor of 2 when the conjugate gradient iteration count is large., Comment: 7 pages, 1 figure, Lattice 2011

As part of our ongoing investigations of QCD with many flavors of quarks, here we report on studies of the finite temperature phase transition for eight-flavor QCD with the DBW2 gauge action and na\"ive staggered fermions. We find a clear first order phase transition between the chirally asymmetric phase at zero temperature and the chirally symmetric phase at finite temperature, signaled by a two-state signal for $\langle\bar{\psi}\psi\rangle$ at a non-zero temperature. We see this signal at a gauge coupling of $\beta=0.54$, where, to set the scale, the zero temperature value for $f_\pi$, in the chiral limit, is 0.06661(92). This strong, first-order signal is seen for two different values of the quark mass, $m_q=0.007$ and 0.0195, at $N_\tau=8$ and 6 respectively. Using $f_\pi(m_q)$ as the scale, the critical temperature is measured to be $T_c/f_\pi=1.638(93)$ at $m_\pi/f_\pi=3.329(30)$ for $m_q=0.007$, and $T_c/f_\pi=1.779(27)$ at $m_\pi/f_\pi=4.093(15)$ for $m_q=0.0195$. At a weaker coupling $\beta=0.56$, where at zero temperature and in the chiral limit we find $f_\pi=0.0312(10)$, the first order signal becomes numerically invisible to us for the $N_\tau \leq 14$ lattices we have investigated so far., Comment: 7 pages, 8 figures, 1 table, Talk presented at The XXVIII International Symposium on Lattice Field Theory, LATTICE 2010 - Villasimius, Sardinia Italy / June 14-19, 2010

We compare extensive simulations of QCD with 8 and 12 flavors of degenerate quarks, using the DBW2 gauge action, naive staggered fermions, and the rational hybrid Monte Carlo algorithm. A variety of values of the coupling constant, quark mass, and lattice size have been used. Our data suggests that, as the bare coupling is decreased, a rapid cross-over, which dramatically changes the lattice scale, exists with both 8 and 12 flavors. The scale change across this cross-over is much larger with 12 flavors than it is with 8 flavors. All of the observables we have measured, in both the zero and finite temperature systems, are consistent with a chiral symmetry breaking phase for the zero temperature theory on the weak coupling side of the rapid cross-over., Comment: 7 pages, 8 figures, Talk presented at The XXVII International Symposium on Lattice Field Theory, LATTICE 2009 - Peking University, Beijing, China / July 26 - 31, 2009

We study the use of NLO and NNLO formulae from SU(2) chiral perturbation theory to fit results from the 2+1 flavor DWF QCD ensembles that have been generated by the RBC and UKQCD collaborations. These ensembles are at two different lattice spacings, contain multiple dynamical light quark masses, and include a variety of partially quenched valence quark masses. Both NLO and complete NNLO fits well represent our data, which has m_pi in the range 220 to 420 MeV. With our data, the NNLO fits have NLO and NNLO contributions of similar size, making the series not convergent and the extrapolation to physical light quark masses imprecise. Thus, we use NLO fit results for our predictions of f_pi, f_K and the light quark masses., Comment: Proceedings contribution for Lattice 2009

At stronger gauge-field couplings, the domain wall fermion (DWF) residual mass, a measure of chiral symmetry breaking, grows rapidly. This measure is largely due to near zero fermion eigenmodes of logarithm of the 4D transfer matrix along the fifth dimension, and these eigenmodes increase rapidly at strong coupling. To suppress these eigenmodes, we have added to the DWF path integral a multiplicative weighting factor consisting of a ratio of determinants of Wilson-Dirac fermions having a chirally twisted mass with a large negative real component and a small imaginary chiral component. Numerical results show that this weighting factor with an appropriate choice of twisted masses significantly suppresses the residual mass while allowing adequate topological tunneling., Comment: PDFLaTeX, 8 pages, 4 Pdf figures, uses JHEP.bst, POS.cls, and POSlogo.pdf

We report on simulations of QCD with many flavors of degenerate quarks, the DBW2 gauge action and naive staggered fermions, using the rational hybrid Monte Carlo algorithm. We primarily focus on eight degenerate quark flavors where a variety of values of the coupling constant and quark mass have been used in the simulations. The scaling behavior of the hadron spectrum and the string tension of the heavy quark potential is studied, to probe whether the zero temperature, continuum limit of the theory breaks chiral symmetry., Comment: 7 pages, 6 figures, Talk presented at The XXVI International Symposium on Lattice Field Theory, LATTICE 2008 - Williamsburg, Virginia, USA / July 14 - 19, 2008

We present results for the dependence of the residual mass of domain wall fermions (DWF) on the size of the fifth dimension and its relation to the density and localization properties of low-lying eigenvectors of the corresponding hermitian Wilson Dirac operator relevant to simulations of 2+1 flavor domain wall QCD. Using the DBW2 and Iwasaki gauge actions, we generate ensembles of configurations with a $16^3\times 32$ space-time volume and an extent of 8 in the fifth dimension for the sea quarks. We demonstrate the existence of a regime where the degree of locality, the size of chiral symmetry breaking and the rate of topology change can be acceptable for inverse lattice spacings $a^{-1} \ge 1.6$ GeV., Comment: 59 Pages, 23 figures, 1 MPG linked

Numerical simulations of full QCD on anisotropic lattices provide a convenient way to study QCD thermodynamics with fixed physics scales and reduced lattice spacing errors. We report results from calculations with two flavors of dynamical staggered fermions, where all bare parameters and the renormalized anisotropy are kept constant and the temperature is changed in small steps by varying only the number of time slices. Including results from zero-temperature scale setting simulations, which determine the Karsch coefficients, allows for the calculation of the equation of state at finite temperatures., Comment: 27 pages, 11 figures, 8 tables

We are using domain wall fermions to study $K \to \pi \pi$ matrix elements by measuring $K \to \pi$ and $K \to 0$ matrix elements on the lattice and employing chiral perturbation theory to relate these to the desired physical result. The residual chiral symmetry breaking of domain wall fermions with a finite extent in the fifth dimension impacts these measurements. Using the Ward-Takahashi identities, we investigate residual chiral symmetry breaking effects for divergent quantities and study pathologies of the quenched approximation for small quark mass. We then discuss the $\Delta S = 1$ operator $O_2$, where chiral symmetry is vital for the subtraction of unphysical effects., Comment: 4 pages, 3 figures, Lattice 2000 (Hadronic Matrix Elements), RBC Collaboration, corrected equations 2 and 7

The QCDSP computer (Quantum Chromodynamics on Digital Signal Processors) is an inexpensive, massively parallel computer intended primarily for simulations in lattice gauge theory. Currently, two large QCDSP machines are in full-time use: an 8,192 processor, 0.4 Teraflops machine at Columbia University and an 12,288 processor, 0.6 Teraflops machine at the RIKEN-BNL Research Center at Brookhaven National Laboratory. We describe the design process, architecture, software and current physics projects of these computers., Comment: 19 pages, 3 figures

One of the major frontiers of lattice field theory is the inclusion of light fermions in simulations, particularly in pursuit of accurate, first principles predictions from lattice QCD. With dedicated Teraflops-scale computers currently simulating QCD, another step towards precision full QCD simulations is underway. In addition to ongoing staggered and Wilson fermion simulations, first results from full QCD with domain wall fermions are available. After some discussion of work toward better algorithms, simulations completed to date will be discussed., Comment: LATTICE99 (plenary) - 10 pages, 11 figures

A serious difficulty in conventional lattice field theory calculations is the coupling between the chiral and continuum limits. With both staggered and Wilson fermions, the chiral limit cannot be realized without first taking the limit of vanishing lattice spacing. In this talk, we report on extensive studies of the domain wall formulation of lattice fermions, which avoids this difficulty at the expense of requiring that fermion propagators be computed in five dimensions. A variety of results will be described for quenched and dynamical simulations at both zero and finite temperature. Conclusions about the benefits of this new method and some new physical results will be presented. These results were obtained on the QCDSP machines recently put into operation at Columbia and the RIKEN Brookhaven Research Center., Comment: Talk presented at ICHEP'98, July 23-29, 1998, Vancouver, Canada. 7 pages, 9 figures Omission from electronic author list corrected

We examine the chiral limit of domain wall fermions in quenched QCD. One expects that in a quenched simulation, exact fermion zero modes will give a divergent, 1/m behavior in the chiral condensate for sufficiently small valence quark masses. Unlike other fermion formulations, domain wall fermions clearly demonstrate this behavior., Comment: LATTICE98(spectrum), G. R. Fleming presented talk, 5 pages, 3 figures, corrected typos in printed version

We describe the completed 8,192-node, 0.4Tflops machine at Columbia as well as the 12,288-node, 0.6Tflops machine assembled at the RIKEN Brookhaven Research Center. Present performance as well as our experience in commissioning these large machines is presented. We outline our on-going physics program and explain how the configuration of the machine is varied to support a wide range of lattice QCD problems, requiring a variety of machine sizes. Finally a brief discussion is given of future prospects for large-scale lattice QCD machines., Comment: LATTICE98(machines), 3 pages, 1 picture, 1 figure

We study the anomalous breaking of U_A(1) symmetry just above the QCD phase transition for zero and two flavors of quarks, using a staggered fermion, lattice discretization. The properties of the QCD phase transition are expected to depend on the degree of U_A(1) symmetry breaking in the transition region. For the physical case of two flavors, we carry out extensive simulations on a 16^3 x 4 lattice, measuring a difference in susceptibilities which is sensitive to U_A(1) symmetry and which avoids many of the staggered fermion discretization difficulties. The results suggest that anomalous effects are at or below the 15% level., Comment: 10 pages including 2 figures and 1 table

The QCDSP machine at Columbia University has grown to 2,048 nodes achieving a peak speed of 100 Gigaflops. Software for quenched and Hybrid Monte Carlo (HMC) evolution schemes has been developed for staggered fermions, with support for Wilson and clover fermions under development. We provide an overview of the runtime environment, the current status of the QCDSP construction program and preliminary results not presented elsewhere in these proceedings., Comment: 3 pages, 1 figure, LaTeX with espcrc2.sty, talk presented at LATTICE97(Poster) by G. R. Fleming

We have done simulations of lattice QCD with different numbers of light dynamical quarks. Since we cannot reach the continuum limit with our current computers, we have done comparisons with 0 (quenched), 2 and 4 light quark flavors with the physical volume and lattice spacing constant, when these are determined from the mass of the rho. We find a 7% (2 standard deviation) difference in the nucleon to rho mass ratio for 2 and 4 quarks. More importantly, the effects of chiral symmetry breaking are dramatically decreased for the case of 4 light quarks., Comment: 15 pages, 6 figures. Presented at RHIC Summer Study '96

We present a summary of the progress on QCDSP in the last year. QCDSP, Quantum Chromodynamics on Digital Signal Processors, is an inexpensive computer being built at Columbia that can achieve 0.8 teraflops for three million dollars., Comment: 4 pages, 1 figure

We have measured the hadron spectrum in lattice QCD, using staggered fermions, for 0 (the quenched approximation), 2 and 4 light degenerate dynamical quarks. In addition to earlier results involving extrapolations in valence quark masses for fixed dynamical mass, we also report results where we extrapolate in the dynamical mass for 4 flavors. We see a marked difference in the hadron spectrum for 2 and 4 flavors; the hadron spectrum is nearly parity doubled for 4 flavors, indicating smaller effects of chiral symmetry breaking. This pronounced effect in the hadron spectrum cannot be removed by a simple change in scale as the number of light quark flavors is changed. Further simulations at larger volume are needed to rule out finite volume effects., Comment: 5 pages, 4 figures. Presented at Lattice QCD on Parallel Computers, Tsukuba, Japan, march 1997

We have studied the hadron spectrum while varying the number of light dynamical quarks when the physical lattice spacing and volume are held fixed relative to the rho mass. For two and zero flavors of staggered fermions, we find the nucleon to rho mass ratios (extrapolated to zero valence quark mass) are very similar. However, for four flavors the ratio is 7% (2 standard deviations) above the two flavor result., Comment: 3 pages, 3 figures

We have measured the valence pion mass and the valence chiral condensate on lattice configurations generated with and without dynamical fermions. We find that our data and that of others is well represented by a linear relationship between $m_{\pi}^2$ and the valence quark mass, with a non-zero intercept. For our data, we relate the intercept to finite volume effects visible in the valence chiral condensate. We see no evidence for the singular behavior expected from quenched chiral logarithms., Comment: 5 pages, Latex with included macro. 4 encapsulated postscript figures included. Contribution to Lattice '95

We discuss a class of saddle-point configurations in SU(2) lattice gauge theory in three Euclidean dimensions. These configurations are smooth on the scale of the lattice and have an action density exhibiting localized peaks, as has been seen in cooled and extremized Monte Carlo generated lattices. Large Wilson loops centered on the action peaks show a unit of Z(2) flux. We discuss the generation of these configurations and measurements of the Creutz ratios on them., Comment: Contribution to Confinement '95, March 1995, Osaka, Japan. 5 page uuencoded postscript file

The current status of United States projects pursuing Teraflops-scale computing resources for lattice field theory is discussed. Two projects are in existence at this time: the Multidisciplinary Teraflops Project, incorporating the physicists of the QCD Teraflops Collaboration, and a smaller project, centered at Columbia, involving the design and construction of a 0.8 Teraflops computer primarily for QCD., Comment: Contribution to Lattice 94. 7 pages. Latex source followed by compressed, uuenocded postscript file of the complete paper. Individual figures available from rdm@phys.columbia.edu

QCD with eight flavors is studied on $16^3\times N_t$ lattices with $N_t=4$, 6, 8, 16 and 32, a dynamical quark mass $ma=0.015$ and lattice coupling $\beta=6/g^2$ between 4.5 and 5.0. For $N_t=16$ and 32, hadron masses and screening lengths are computed for a variety of valence quark masses. The previously observed, strong, first-order transition for $N_t=4$, 6 and 8 is seen, for $N_t=16$, to become a $\beta$-independent, zero-temperature transition characterized by a factor of $\approx 3$ change in lattice scale. This strong, first-order transition restores chiral symmetry, at least for $N_t=4$, 6 and 8, producing a chirally symmetric, weak-coupling phase. However, as $N_t$ increases to 16, the chiral symmetry properties of the weak-coupling side of the zero-temperature transition are unclear and offer a hint of a normal, finite-temperature, chiral symmetry breaking transition in the weak-coupling phase., Comment: 36 pages of text plus 17 pages of figures