Search

This draft report summarizes and details the findings, results, and recommendations derived from the ASCR/HEP Exascale Requirements Review meeting held in June, 2015. The main conclusions are as follows. 1) Larger, more capable computing and data facilities are needed to support HEP science goals in all three frontiers: Energy, Intensity, and Cosmic. The expected scale of the demand at the 2025 timescale is at least two orders of magnitude -- and in some cases greater -- than that available currently. 2) The growth rate of data produced by simulations is overwhelming the current ability, of both facilities and researchers, to store and analyze it. Additional resources and new techniques for data analysis are urgently needed. 3) Data rates and volumes from HEP experimental facilities are also straining the ability to store and analyze large and complex data volumes. Appropriately configured leadership-class facilities can play a transformational role in enabling scientific discovery from these datasets. 4) A close integration of HPC simulation and data analysis will aid greatly in interpreting results from HEP experiments. Such an integration will minimize data movement and facilitate interdependent workflows. 5) Long-range planning between HEP and ASCR will be required to meet HEP's research needs. To best use ASCR HPC resources the experimental HEP program needs a) an established long-term plan for access to ASCR computational and data resources, b) an ability to map workflows onto HPC resources, c) the ability for ASCR facilities to accommodate workflows run by collaborations that can have thousands of individual members, d) to transition codes to the next-generation HPC platforms that will be available at ASCR facilities, e) to build up and train a workforce capable of developing and using simulations and analysis to support HEP scientific research on next-generation systems., Comment: 77 pages, 13 Figures; draft report, subject to further revision

These reports present the results of the 2013 Community Summer Study of the APS Division of Particles and Fields ("Snowmass 2013") on the future program of particle physics in the U.S. Chapter 9, on Computing, discusses the computing challenges for future experiments in the Energy, Intensity, and Cosmic Frontiers, for accelerator science, and for particle theory, as well as structural issues in supporting the intense uses of computing required in all areas of particle physics., Comment: 36 pages

These reports present the results of the 2013 Community Summer Study of the APS Division of Particles and Fields ("Snowmass 2013") on the future program of particle physics in the U.S. Chapter 1 contains the Executive Summary and the summaries of the reports of the nine working groups., Comment: 51 pages

This report summarizes the findings of the DPF Theory Panel which was formed with a goal of understanding the scientific problems and opportunities of the next decade, as well as the challenges involved in sustaining a first-class program in theoretical particle physics research in the United States., Comment: 23 pages, LaTeX, Chapter 5 of the report "Planning the Future of U.S. Particle Physics (Snowmass 2013)"

Over the past five years, graphics processing units (GPUs) have had a transformational effect on numerical lattice quantum chromodynamics (LQCD) calculations in nuclear and particle physics. While GPUs have been applied with great success to the post-Monte Carlo "analysis" phase which accounts for a substantial fraction of the workload in a typical LQCD calculation, the initial Monte Carlo "gauge field generation" phase requires capability-level supercomputing, corresponding to O(100) GPUs or more. Such strong scaling has not been previously achieved. In this contribution, we demonstrate that using a multi-dimensional parallelization strategy and a domain-decomposed preconditioner allows us to scale into this regime. We present results for two popular discretizations of the Dirac operator, Wilson-clover and improved staggered, employing up to 256 GPUs on the Edge cluster at Lawrence Livermore National Laboratory., Comment: 11 pages, 10 figures, to appear in the proceedings of the 2011 ACM/IEEE International Conference for High Performance Computing, Networking, Storage and Analysis (SC'11)

The separation of a heavy quark and antiquark pair leads to the formation of a tube of flux, or "string", which should break in the presence of light quark-antiquark pairs. This expected zero-temperature phenomenon has proven elusive in simulations of lattice QCD. We study mixing between the string state and the two-meson decay channel in QCD with two flavors of dynamical sea quarks. We confirm that mixing is weak and find that it decreases at level crossing. While our study does not show direct effects of internal quark loops, our results, combined with unitarity, give clear confirmation of string breaking., Comment: 20 pages, 7 figures. With small clarifications and two additions to references. Submitted to Phys. Rev. D

We present results for pseudoscalar decay constants of heavy-light mesons using both quenched and N_f = 2 dynamical fermion configurations. A variety of fermion actions is investigated : Wilson, nonperturbative clover, and fat-link clover. For heavy quarks the Fermilab formalism is applied. In the quenched approximation, results with the nonperturbatively improved clover action of the Alpha collaboration allow us to study the systematic error of the continuum extrapolation from the Wilson action. In addition, we use quenched configurations to explore the effects of fattening. The lessons from the quenched analyses are then applied to data with dynamical fermions, where both Wilson and fat-link clover actions have been used. This allows us to attempt a continuum extrapolation of the dynamical results., Comment: Lattice 2000 (Heavy Quark Physics). 4 pages, 2 figures

We study the effects of dynamical quarks on the static quark potential at distances shorter than those where string breaking is expected. Quenched calculations and calculations with three flavors of dynamical quarks are done on sets of lattices with the lattice spacings matched within about one percent. The effect of the sea quarks on the shape of the potential is clearly visible. We investigate the consequences of these effects in a very crude model, namely solving Schroedinger's equation in the resulting potential., Comment: 20 pages RevTeX, color postscript figures included Table heading fixed, Fig. 4 improved, reference added, lines tightened

Improved lattice actions for Kogut-Susskind quarks have been shown to improve rotational symmetry and flavor symmetry. In this work we find improved scaling behavior of the rho and nucleon masses expressed in units of a length scale obtained from the static quark potential, and better behavior of the Dirac operator in instanton backgrounds., Comment: 4 pages, 4 figures, Revtex

We are studying a variety of semileptonic decays of heavy-light mesons in an effort to improve the determination of the heavy-quark Standard-Model CKM matrix elements. Our fermion action is a novel, improved ``fat'' clover action that promises to reduce problems with exceptional configurations. Dynamical sea quarks are included in a mixed approach, i.e. we use staggered sea quarks and fat-clover valence quarks. Here we report preliminary results., Comment: LATTICE99(heavyqk) - 3p, 4 Postscript figs

Quantum chromodynamics with two zero-mass flavors is expected to exhibit a phase transition with O(4) critical behavior. Fixing the universality class is important for phenomenology and for facilitating the extrapolation of simulation data to physical quark mass values. Other groups have reported results from lattice QCD simulations with dynamical staggered quarks at Nt = 4, which suggest a departure from the expected critical behavior. We have pushed simulations to the largest volumes and smallest quark mass to date. Strong discrepancies in critical exponents and the scaling equation of state persist., Comment: 23 pp RevTeX, 10 Postscript figures, submitted to Phys. Rev

We present an update on the MILC Collaboration's light hadron spectrum calculation with two flavors of dynamical, staggered quarks. We present extrapolations of the nucleon to rho mass ratio to the continuum limit for fixed values of the pi to rho mass ratio including the physical one., Comment: 3 pages, LaTex with espcrc2 and epsf, 5 postscript figures included, Lattice '97 Proceedings

Quantum chromodynamics with two zero mass flavors is expected to exhibit a phase transition with O(4) critical behavior. Fixing the universality class is important for phenomenology and for facilitating the extrapolation of simulation data to physical quark mass values. At Lattice '96 the Tsukuba and Bielefeld groups reported results from new simulations with dynamical staggered quarks at $N_t = 4$, which suggested a departure from the expected critical behavior. We report observations of similar deviations and discuss efforts in progress to understand this phenomenon., Comment: 3 pp, LaTeX with 6 encapsulated Postscript figures. Lattice '97 proceedings

We calculate the two flavor equation of state for QCD on lattices with lattice spacing a=(6T)^{-1} and find that cutoff effects are substantially reduced compared to an earlier study using a=(4T)^{-1}. However, it is likely that significant cutoff effects remain. We fit the lattice data to expected forms of the free energy density for a second order phase transition at zero-quark-mass, which allows us to extrapolate the equation of state to m_q=0 and to extract the speed of sound. We find that the equation of state depends weakly on the quark mass for small quark mass., Comment: 24 pages, latex, 11 postscipt figures

We have extended our earlier simulations of the high temperature behaviour of lattice QCD with two light flavours of staggered quarks on a $16^3 \times 8$ lattice to lower quark mass (m_q=0.00625). The transition from hadronic matter to a quark-gluon plasma is observed at $6/g^2=5.49(2)$ corresponding to a temperature of $T_c \approx 140$MeV. We present measurements of observables which probe the nature of the quark-gluon plasma and serve to distinguish it from hadronic matter. Although the transition is quite abrupt, we have seen no indications that it is first order., Comment: 23 pages, RevteX, 6 encapsulated postscript figures

We have carried out spectrum calculations with two flavors of dynamical Kogut-Susskind quarks on four lattice sizes from $8^3\times 24$ to $16^3\times24$ at couplings that correspond to chiral symmetry restoration for a lattice with 6 time slices. We estimate that the linear spatial sizes of the lattices range from 1.8 to 3.6 fm. We find significant finite size effects for all particles between the smallest and largest volume with the larger quark mass that we study, $am_q=0.025$, where $a$ is the lattice spacing. The nucleon experiences the largest effect of about 6 percent. We also study a lighter quark mass, $am_q=0.0125$, on the two largest lattices. Effects of the dynamical and valence quark masses on the hadron spectrum are studied both directly, by comparing the two simulations, and by extracting mass derivatives from the correlation functions. We do not find much improvement in the nucleon to rho mass ratio as we decrease the quark mass at this lattice spacing. Finally, we report on an unsuccessful attempt to see effects of the $\rho \rightarrow 2\pi$ decay on the $\rho$ mass, and on studies of Wilson and Kogut-Susskind hadron masses with large valence quark masses. (The paper is also available via anonymous ftp in a compressed PostScript file: fuji.physics.indiana.edu:/pub/milc/paper.ps.Z .), Comment: IUHET-232, 54 pages. Posted here is the PostScript file including 15 embedded figures

We have measured some simple matrix elements for pseudoscalar and vector mesons made of Wilson valence quarks and staggered sea quarks at $\beta=5.6$ at sea quark masses $am_q=0.01$ and 0.025. Our measurements include the decay constants of pseudoscalars (including $f_D$), the wave function at the origin (or decay constant) of vector mesons, and the calculation of quark masses from current algebra. The effects of sea quarks on the simulations are small. We make comparisons to quenched simulations at similar values of the lattice spacing ($1/a \simeq 2$ GeV).

As part of an ongoing effort to characterize the high temperature phase of QCD, we measure the quark baryon density in the vicinity of a fixed test quark and compare it with similar measurements at low temperature and at the crossover temperature. Such an observable has also been studied by the Vienna group. We find an extremely weak correlation at high temperature, suggesting that small color singlet clusters are unimportant in the thermal ensemble. We also find that at T = 0.75 T_c the induced quark number shows a surprisingly large component attributable to baryonic screening. A simulation of a simple flux tube model produces results that suggest a plausible scenario: As the crossover temperature is approached from below, baryonic states proliferate. Above the crossover temperature the mean size of color singlet clusters grows explosively, resulting in an effective electrostatic deconfinement., Comment: LaTeX document. Requires Elsevier espcrc2 style macro. 4 pages with 4 Postscript figures bundled in a Unix shar file. Lattice '92 talk by C. DeTar

We present an analysis of hadronic spectroscopy for Wilson valence quarks with dynamical staggered fermions at lattice coupling $6/g^2 = \beta=5.6$ at sea quark mass $am_q=0.01$ and 0.025, and of Wilson valence quarks in quenched approximation at $\beta=5.85$ and 5.95, both on $16^3 \times 32$ lattices. We make comparisons with our previous results with dynamical staggered fermions at the same parameter values but on $16^4$ lattices doubled in the temporal direction., Comment: 32 pages

Lattice QCD with 2 light staggered quark flavours is being simulated on a $16^3\times8$ lattice to study the transition from hadronic matter to a quark gluon plasma. We have completed runs at $m_q=0.0125$ and are extending this to $m_q=0.00625$. We also examine the addition of a non-dynamical "strange" quark. Thermodynamic order parameters are being measured across the transition and further into the plasma phase, as are various screening lengths. No evidence for a first order transition is seen, and we estimate the transition temperature to be $TY_c=143(7) MeV$.