
Profiles of dark matterdominated halos at the group and cluster scales play
an important role in modern cosmology. Using results from two very large
cosmological $N$body simulations, which increase the available volume at their
mass resolution by roughly two orders of magnitude, we robustly determine the
halo concentrationmass $(cM)$ relation over a wide range of masses, employing
multiple methods of concentration measurement. We characterize individual halo
profiles, as well as stacked profiles, relevant for galaxygalaxy lensing and
nextgeneration cluster surveys; the redshift range covered is $0\leq z \leq
4$, with a minimum halo mass of $M_{200c}\sim2\times10^{11} M_\odot$. Despite
the complexity of a proper description of a halo (environmental effects, merger
history, nonsphericity, relaxation state), when the mass is scaled by the
nonlinear mass scale $M_\star(z)$, we find that a simple nonpowerlaw form for
the $cM/M_\star$ relation provides an excellent description of our simulation
results across eight decades in $M/M_{\star}$ and for $0\leq z \leq 4$. Over
the mass range covered, the $cM$ relation has two asymptotic forms: an
approximate power law below a mass threshold $M/M_\star\sim 5001000$,
transitioning to a constant value, $c_0\sim 3$ at higher masses. The relaxed
halo fraction decreases with mass, transitioning to a constant value of $\sim
0.5$ above the same mass threshold. We compare NavarroFrenkWhite (NFW) and
Einasto fits to stacked profiles in narrow mass bins at different redshifts; as
expected, the Einasto profile provides a better description of the simulation
results. At cluster scales at low redshift, however, both NFW and Einasto
profiles are in very good agreement with the simulation results, consistent
with recent weak lensing observations.

Report of the first workshop to identify approaches and techniques in the
domain of quantum sensing that can be utilized by future High Energy Physics
applications to further the scientific goals of High Energy Physics.

The use of highquality simulated sky catalogs is essential for the success
of cosmological surveys. The catalogs have diverse applications, such as
investigating signatures of fundamental physics in cosmological observables,
understanding the effect of systematic uncertainties on measured signals and
testing mitigation strategies for reducing these uncertainties, aiding analysis
pipeline development and testing, and survey strategy optimization. The list of
applications is growing with improvements in the quality of the catalogs and
the details that they can provide. Given the importance of simulated catalogs,
it is critical to provide rigorous validation protocols that enable both
catalog providers and users to assess the quality of the catalogs in a
straightforward and comprehensive way. For this purpose, we have developed the
DESCQA framework for the Large Synoptic Survey Telescope Dark Energy Science
Collaboration as well as for the broader community. The goal of DESCQA is to
enable the inspection, validation, and comparison of an inhomogeneous set of
synthetic catalogs via the provision of a common interface within an automated
framework. In this paper, we present the design concept and first
implementation of DESCQA. In order to establish and demonstrate its full
functionality we use a set of interim catalogs and validation tests. We
highlight several important aspects, both technical and scientific, that
require thoughtful consideration when designing a validation framework,
including validation metrics and how these metrics impose requirements on the
synthetic sky catalogs.

We present the clustering measurements of quasars in configuration space
based on the Data Release 14 (DR14) of the Sloan Digital Sky Survey IV extended
Baryon Oscillation Spectroscopic Survey. This dataset includes 148,659 quasars
spread over the redshift range $0.8\leq z \leq 2.2$ and spanning 2112.9 square
degrees. We use the Convolution Lagrangian Perturbation Theory (CLPT) approach
with a Gaussian Streaming (GS) model for the redshift space distortions of the
correlation function and demonstrate its applicability for dark matter halos
hosting eBOSS quasar tracers. At the effective redshift $z_{\rm eff} = 1.52$,
we measure the linear growth rate of structure $f\sigma_{8}(z_{\rm eff})= 0.426
\pm 0.077$, the expansion rate $H(z_{\rm eff})= 159^{+12}_{13}(r_{s}^{\rm
fid}/r_s){\rm km.s}^{1}.{\rm Mpc}^{1}$, and the angular diameter distance
$D_{A}(z_{\rm eff})=1850^{+90}_{115}\,(r_s/r_{s}^{\rm fid}){\rm Mpc}$, where
$r_{s}$ is the sound horizon at the end of the baryon drag epoch and
$r_{s}^{\rm fid}$ is its value in the fiducial cosmology. The quoted errors
include both systematic and statistical contributions. The results on the
evolution of distances are consistent with the predictions of flat
$\Lambda$Cold Dark Matter ($\Lambda$CDM) cosmology with Planck parameters,
and the measurement of $f\sigma_{8}$ extends the validity of General Relativity
(GR) to higher redshifts($z>1$) This paper is released with companion papers
using the same sample. The results on the cosmological parameters of the
studies are found to be in very good agreement, providing clear evidence of the
complementarity and of the robustness of the first fullshape clustering
measurements with the eBOSS DR14 quasar sample.

We explore the cosmological implications of anisotropic clustering
measurements of the quasar sample from Data Release 14 of the Sloan Digital Sky
Survey IV Extended Baryon Oscillation Spectroscopic Survey (eBOSS) in
configuration space. The $\sim 147,000$ quasar sample observed by eBOSS offers
a direct tracer of the density field and bridges the gap of previous BAO
measurements between redshift $0.8<z<2.2$. By analysing the twopoint
correlation function characterized by clustering wedges $\xi_{\rm w_i}(s)$ and
multipoles $\xi_{\ell}(s)$, we measure the angular diameter distance, Hubble
parameter and cosmic structure growth rate. We define a systematic error budget
for our measurements based on the analysis of $N$body simulations and mock
catalogues. Based on the DR14 large scale structure quasar sample at the
effective redshift $z_{\rm eff}=1.52$, we find the growth rate of cosmic
structure $f\sigma_8(z_{\rm eff})=0.396\pm 0.079$, and the geometric parameters
$D_{\rm V}(z)/r_{\rm d}=26.47\pm 1.23$, and $F_{\rm AP}(z)=2.53\pm 0.22$, where
the uncertainties include both statistical and systematic errors. These values
are in excellent agreement with the bestfitting standard ${\rm \Lambda CDM}$
model to the latest cosmic microwave background data from Planck.

We analyse the clustering of the Sloan Digital Sky Survey IV extended Baryon
Oscillation Spectroscopic Survey Data Release 14 quasar sample (DR14Q). We
measure the redshift space distortions using the power spectrum monopole,
quadrupole and hexadecapole inferred from 148,659 quasars between redshifts 0.8
and 2.2 covering a total sky footprint of 2112.9 deg$^2$. We constrain the
logarithmic growth of structure times the amplitude of dark matter density
fluctuations, $f\sigma_8$, and the AlcockPaczynski dilation scales which allow
constraints to be placed on the angular diameter distance $D_A(z)$ and the
Hubble $H(z)$ parameter. At the effective redshift of $z_{\rm eff}=1.52$,
$f\sigma_8(z_{\rm eff})=0.420\pm0.076$, $H(z_{\rm eff})=[162\pm 12]\, (r_s^{\rm
fid}/r_s)\,{\rm km\, s}^{1}{\rm Mpc}^{1}$, and $D_A(z_{\rm eff})=[1.85\pm
0.11]\times10^3\,(r_s/r_s^{\rm fid})\,{\rm Mpc}$, where $r_s$ is the comoving
sound horizon at the baryon drag epoch and the superscript `fid' stands for its
fiducial value. The errors take into account the full error budget, including
systematics and statistical contributions. These results are in full agreement
with the current $\Lambda$Cold Dark Matter ($\Lambda$CDM) cosmological model
inferred from Planck measurements.
Finally, we compare our measurements with other eBOSS companion papers and
find excellent agreement, demonstrating the consistency and complementarity of
the different methods used for analysing the data.

Spatially coupled lowdensity paritycheck (SCLDPC) codes are sparse graph
codes that have recently become of interest due to their capacityapproaching
performance on memoryless binary input channels. In this paper, we unify all
existing SCLDPC code construction methods under a new generalized description
of SCLDPC codes based on algebraic lifts of graphs. We present an improved
lowcomplexity counting method for the special case of $(3,3)$absorbing sets
for arraybased SCLDPC codes, which we then use to optimize permutation
assignments in SCLDPC code construction. We show that codes constructed in
this way are able to outperform previously published constructions, in terms of
the number of dominant absorbing sets and with respect to both standard and
windowed decoding.

We introduce a new cosmic emulator for the matter power spectrum covering
eight cosmological parameters. Targeted at optical surveys, the emulator
provides accurate predictions out to a wavenumber k~5/Mpc and redshift z<=2.
Besides covering the standard set of LCDM parameters, massive neutrinos and a
dynamical dark energy of state are included. The emulator is built on a sample
set of 36 cosmological models, carefully chosen to provide accurate predictions
over the wide and large parameter space. For each model, we have performed a
highresolution simulation, augmented with sixteen mediumresolution
simulations and TimeRG perturbation theory results to provide accurate coverage
of a wide krange; the dataset generated as part of this project is more than
1.2Pbyte. With the current set of simulated models, we achieve an accuracy of
approximately 4%. Because the sampling approach used here has established
convergence and errorcontrol properties, followon results with more than a
hundred cosmological models will soon achieve ~1% accuracy. We compare our
approach with other prediction schemes that are based on halo model ideas and
remapping approaches. The new emulator code is publicly available.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV groundbased dark
energy experiment that will study baryon acoustic oscillations (BAO) and the
growth of structure through redshiftspace distortions with a widearea galaxy
and quasar redshift survey. To trace the underlying dark matter distribution,
spectroscopic targets will be selected in four classes from imaging data. We
will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to
even higher redshift, DESI will target bright [O II] emission line galaxies up
to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying
dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the
Ly$\alpha$ forest absorption features in their spectra, which will be used to
trace the distribution of neutral hydrogen. When moonlight prevents efficient
observations of the faint targets of the baseline survey, DESI will conduct a
magnitudelimited Bright Galaxy Survey comprising approximately 10 million
galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy
and quasar redshifts will be obtained to measure the BAO feature and determine
the matter power spectrum, including redshift space distortions.

DESI (Dark Energy Spectropic Instrument) is a Stage IV groundbased dark
energy experiment that will study baryon acoustic oscillations and the growth
of structure through redshiftspace distortions with a widearea galaxy and
quasar redshift survey. The DESI instrument is a roboticallyactuated,
fiberfed spectrograph capable of taking up to 5,000 simultaneous spectra over
a wavelength range from 360 nm to 980 nm. The fibers feed ten threearm
spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500,
depending on wavelength. The DESI instrument will be used to conduct a
fiveyear survey designed to cover 14,000 deg$^2$. This powerful instrument
will be installed at prime focus on the 4m Mayall telescope in Kitt Peak,
Arizona, along with a new optical corrector, which will provide a threedegree
diameter field of view. The DESI collaboration will also deliver a
spectroscopic pipeline and data management system to reduce and archive all
data for eventual public use.

Gravitational lensing has become one of the most powerful tools available for
investigating the 'dark side' of the universe. Cosmological strong
gravitational lensing, in particular, probes the properties of the dense cores
of dark matter halos over decades in mass and offers the opportunity to study
the distant universe at flux levels and spatial resolutions otherwise
unavailable. Studies of stronglylensed variable sources offer yet further
scientific opportunities. One of the challenges in realizing the potential of
strong lensing is to understand the statistical context of both the individual
systems that receive extensive followup study, as well as that of the larger
samples of strong lenses that are now emerging from survey efforts. Motivated
by these challenges, we have developed an imagesimulation pipeline, PICS
(Pipeline for Images of Cosmological Strong lensing) to generate realistic
strong gravitational lensing signals from group and cluster scale lenses. PICS
uses a lownoise and unbiased density estimator based on (resampled) Delaunay
Tessellations to calculate the density field; lensed images are produced by
raytracing images of actual galaxies from deep Hubble Space Telescope
observations. Other galaxies, similarly sampled, are added to fill in the light
cone. The pipeline further adds clustermember galaxies and foreground stars
into the lensed images. The entire image ensemble is then observed using a
realistic point spread function which includes appropriate detector artifacts
for bright stars. Noise is further added, including such nonGaussian elements
as noise windowpaning from mosaiced observations, residual bad pixels, and
cosmic rays. The aim is to produced simulated images that appear identical  to
the eye (expert or otherwise)  to real observations in various imaging
surveys.

The pairwise kinematic SunyaevZel'dovich (kSZ) signal from galaxy clusters
is a probe of their lineofsight momenta, and thus a potentially valuable
source of cosmological information. In addition to the momenta, the amplitude
of the measured signal depends on the properties of the intracluster gas and
observational limitations such as errors in determining cluster centers and
redshifts. In this work we simulate the pairwise kSZ signal of clusters at z<1,
using the output from a cosmological Nbody simulation and including the
properties of the intracluster gas via a model that can be varied in
postprocessing. We find that modifications to the gas profile due to star
formation and feedback reduce the pairwise kSZ amplitude of clusters by ~50%,
relative to the naive 'gas traces mass' assumption. We demonstrate that
miscentering can reduce the overall amplitude of the pairwise kSZ signal by up
to 10%, while redshift errors can lead to an almost complete suppression of the
signal at small separations. We confirm that a highsignificance detection is
expected from the combination of data from currentgeneration, highresolution
CMB experiments, such as the South Pole Telescope, and cluster samples from
optical photometric surveys, such as the Dark Energy Survey. Furthermore, we
forecast that future experiments such as Advanced ACTPol in conjunction with
data from the Dark Energy Spectroscopic Instrument will yield detection
significances of at least 20{\sigma}, and up to 57{\sigma} in an optimistic
scenario. Our simulated maps are publicly available at:
http://www.hep.anl.gov/cosmology/ksz.html

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
leadershipclass 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) Longrange 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 longterm 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 nextgeneration 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
nextgeneration systems.

Largescale structure surveys in the coming years will measure the
redshiftspace power spectrum to unprecedented accuracy, allowing for powerful
new tests of the LambdaCDM picture as well as measurements of particle physics
parameters such as the neutrino masses. We extend the TimeRG perturbative
framework to redshift space, computing the power spectrum P_s(k,mu) in massive
neutrino cosmologies with timedependent dark energy equations of state w(z).
TimeRG is uniquely capable of incorporating scaledependent growth into the
P_s(k,mu) computation, which is important for massive neutrinos as well as
modified gravity models. Although changes to w(z) and the neutrino mass
fraction both affect the latetime scaledependence of the nonlinear power
spectrum, we find that the two effects depend differently on the lineofsight
angle mu. Finally, we use the HACC Nbody code to quantify errors in the
perturbative calculations. For a LambdaCDM model at redshift z=1, our procedure
predicts the monopole~(quadrupole) to 1% accuracy up to a wave number 0.19h/Mpc
(0.28h/Mpc), compared to 0.08h/Mpc (0.07h/Mpc) for the Kaiser approximation and
0.19h/Mpc (0.16h/Mpc) for the current stateoftheart perturbation scheme. Our
calculation agrees with the simulated redshiftspace power spectrum even for
neutrino masses above the current bound, and for rapidlyevolving dark energy
equations of state, dw/dz ~ 1. Along with this article, we make our
redshiftspace TimeRG implementation publicly available as the code redTime.

Precision measurements of the large scale structure of the Universe require
large numbers of high fidelity mock catalogs to accurately assess, and account
for, the presence of systematic effects. We introduce and test a scheme for
generating mock catalogs rapidly using suitably derated Nbody simulations. Our
aim is to reproduce the large scale structure and the gross properties of dark
matter halos with high accuracy, while sacrificing the details of the internal
structure of the halos. By adjusting global and local timesteps in an Nbody
code, we demonstrate that we recover halo masses to better than 0.5% and the
power spectrum to better than 1% both in real and redshift space for k =
1h/Mpc, while requiring a factor of 4 less CPU time. We also calibrate the
redshift spacing of outputs required to generate simulated light cones. We find
that outputs separated by every z = 0.05 allow us to interpolate particle
positions and velocities to reproduce the real and redshift space power spectra
to better than 1% (out to k = 1h/Mpc). We apply these ideas to generate a suite
of simulations spanning a range of cosmologies, motivated by the Baryon
Oscillation Spectroscopic Survey (BOSS) but broadly applicable to future large
scale structure surveys including eBOSS and DESI. As an initial demonstration
of the utility of such simulations, we calibrate the shift in the baryonic
acoustic oscillation peak position as a function of galaxy bias with higher
precision than has been possible so far. This paper also serves to document the
simulations, which we make publicly available.

Computing plays an essential role in all aspects of high energy physics. As
computational technology evolves rapidly in new directions, and data throughput
and volume continue to follow a steep trendline, it is important for the HEP
community to develop an effective response to a series of expected challenges.
In order to help shape the desired response, the HEP Forum for Computational
Excellence (HEPFCE) initiated a roadmap planning activity with two key
overlapping drivers  1) software effectiveness, and 2) infrastructure and
expertise advancement. The HEPFCE formed three working groups, 1) Applications
Software, 2) Software Libraries and Tools, and 3) Systems (including systems
software), to provide an overview of the current status of HEP computing and to
present findings and opportunities for the desired HEP computational roadmap.
The final versions of the reports are combined in this document, and are
presented along with introductory material.

The halo occupation distribution (HOD) approach has proven to be an effective
method for modeling galaxy clustering and bias. In this approach, galaxies of a
given type are probabilistically assigned to individual halos in Nbody
simulations. In this paper, we present a fast emulator for predicting the fully
nonlinear galaxy power spectrum over a range of freely specifiable HOD modeling
parameters. The emulator is constructed using results from 100 HOD models run
on a large LCDM Nbody simulation, with Gaussian Process interpolation applied
to a PCAbased representation of the galaxy power spectrum. The total error is
currently ~3% (~2% in the simulation and ~1% in the emulation process) from z=1
to z=0, over the considered parameter range. We use the emulator to investigate
parametric dependencies in the HOD model, as well as the behavior of galaxy
bias as a function of HOD parameters. The emulator is publicly available at
http://www.hep.anl.gov/cosmology/CosmicEmu/emu.html.

Ground and spacebased sky surveys enable powerful cosmological probes based
on measurements of galaxy properties and the distribution of galaxies in the
Universe. These probes include weak lensing, baryon acoustic oscillations,
abundance of galaxy clusters, and redshift space distortions; they are
essential to improving our knowledge of the nature of dark energy. On the
theory and modeling front, largescale simulations of cosmic structure
formation play an important role in interpreting the observations and in the
challenging task of extracting cosmological physics at the needed precision.
These simulations must cover a parameter range beyond the standard six
cosmological parameters and need to be run at high mass and force resolution.
One key simulationbased task is the generation of accurate theoretical
predictions for observables, via the method of emulation. Using a new sampling
technique, we explore an 8dimensional parameter space including massive
neutrinos and a variable dark energy equation of state. We construct trial
emulators using two surrogate models (the linear power spectrum and an
approximate halo mass function). The new sampling method allows us to build
precision emulators from just 26 cosmological models and to increase the
emulator accuracy by adding new sets of simulations in a prescribed way. This
allows emulator fidelity to be systematically improved as new observational
data becomes available and higher accuracy is required. Finally, using one LCDM
cosmology as an example, we study the demands imposed on a simulation campaign
to achieve the required statistics and accuracy when building emulators for
dark energy investigations.

SPHEREx (SpectroPhotometer for the History of the Universe, Epoch of
Reionization, and Ices Explorer) ( http://spherex.caltech.edu ) is a proposed
allsky spectroscopic survey satellite designed to address all three science
goals in NASA's Astrophysics Division: probe the origin and destiny of our
Universe; explore whether planets around other stars could harbor life; and
explore the origin and evolution of galaxies. SPHEREx will scan a series of
Linear Variable Filters systematically across the entire sky. The SPHEREx data
set will contain R=40 spectra fir 0.75$<\lambda<$4.1$\mu$m and R=150 spectra
for 4.1$<\lambda<$4.8$\mu$m for every 6.2 arc second pixel over the entiresky.
In this paper, we detail the extragalactic and cosmological studies SPHEREx
will enable and present detailed systematic effect evaluations. We also outline
the Ice and Galaxy Evolution Investigations.

Modeling largescale sky survey observations is a key driver for the
continuing development of high resolution, largevolume, cosmological
simulations. We report the first results from the 'Q Continuum' cosmological
Nbody simulation run carried out on the GPUaccelerated supercomputer Titan.
The simulation encompasses a volume of (1300 Mpc)^3 and evolves more than half
a trillion particles, leading to a particle mass resolution of ~1.5 X 10^8
M_sun. At this mass resolution, the Q Continuum run is currently the largest
cosmology simulation available. It enables the construction of detailed
synthetic sky catalogs, encompassing different modeling methodologies,
including semianalytic modeling and subhalo abundance matching in a large,
cosmological volume. Here we describe the simulation and outputs in detail and
present first results for a range of cosmological statistics, such as mass
power spectra, halo mass functions, and halo massconcentration relations for
different epochs. We also provide details on challenges connected to running a
simulation on almost 90% of Titan, one of the fastest supercomputers in the
world, including our usage of Titan's GPU accelerators.

Current and future surveys of largescale cosmic structure are associated
with a massive and complex datastream to study, characterize, and ultimately
understand the physics behind the two major components of the 'Dark Universe',
dark energy and dark matter. In addition, the surveys also probe primordial
perturbations and carry out fundamental measurements, such as determining the
sum of neutrino masses. Largescale simulations of structure formation in the
Universe play a critical role in the interpretation of the data and extraction
of the physics of interest. Just as survey instruments continue to grow in size
and complexity, so do the supercomputers that enable these simulations. Here we
report on HACC (Hardware/Hybrid Accelerated Cosmology Code), a recently
developed and evolving cosmology Nbody code framework, designed to run
efficiently on diverse computing architectures and to scale to millions of
cores and beyond. HACC can run on all current supercomputer architectures and
supports a variety of programming models and algorithms. It has been
demonstrated at scale on Cell and GPUaccelerated systems, standard multicore
node clusters, and Blue Gene systems. HACC's design allows for ease of
portability, and at the same time, high levels of sustained performance on the
fastest supercomputers available. We present a description of the design
philosophy of HACC, the underlying algorithms and code structure, and outline
implementation details for several specific architectures. We show selected
accuracy and performance results from some of the largest high resolution
cosmological simulations so far performed, including benchmarks evolving more
than 3.6 trillion particles.

Over the next decade, cosmological measurements of the largescale structure
of the Universe will be sensitive to the combined effects of dynamical dark
energy and massive neutrinos. The matter power spectrum is a key repository of
this information. We extend higherorder perturbative methods for computing the
power spectrum to investigate these effects over quasilinear scales. Through
comparison with Nbody simulations we establish the regime of validity of a
TimeRenormalization Group (TimeRG) perturbative treatment that includes
dynamical dark energy and massive neutrinos. We also quantify the accuracy of
Standard (SPT), Renormalized (RPT) and Lagrangian Resummation (LPT)
perturbation theories without massive neutrinos. We find that an approximation
that neglects neutrino clustering as a source for nonlinear matter clustering
predicts the Baryon Acoustic Oscillation (BAO) peak position to 0.25% accuracy
for redshifts 1 < z < 3, justifying the use of LPT for BAO reconstruction in
upcoming surveys. We release a modified version of the public Copter code which
includes the additional physics discussed in the paper.

Modern sky surveys are returning precision measurements of cosmological
statistics such as weak lensing shear correlations, the distribution of
galaxies, and cluster abundance. To fully exploit these observations, theorists
must provide predictions that are at least as accurate as the measurements, as
well as robust estimates of systematic errors that are inherent to the modeling
process. In the nonlinear regime of structure formation, this challenge can
only be overcome by developing a largescale, multiphysics simulation
capability covering a range of cosmological models and astrophysical processes.
As a first step to achieving this goal, we have recently developed a prediction
scheme for the matter power spectrum (a socalled emulator), accurate at the 1%
level out to k~1/Mpc and z=1 for wCDM cosmologies based on a set of
highaccuracy Nbody simulations. It is highly desirable to increase the range
in both redshift and wavenumber and to extend the reach in cosmological
parameter space. To make progress in this direction, while minimizing
computational cost, we present a strategy that maximally reuses the original
simulations. We demonstrate improvement over the original spatial dynamic range
by an order of magnitude, reaching k~10 h/Mpc, a fourfold increase in redshift
coverage, to z=4, and now include the Hubble parameter as a new independent
variable. To further the range in k and z, a new set of nested simulations run
at modest cost is added to the original set. The extension in h is performed by
including perturbation theory results within a multiscale procedure for
building the emulator. This economical methodology still gives excellent error
control, ~5% near the edges of the domain of applicability of the emulator. A
public domain code for the new emulator is released as part of the work
presented in this paper.

This document presents (offline) computing requrements and challenges for
Cosmic Frontier science, covering the areas of data management, analysis, and
simulations. We invite contributions to extend the range of covered topics and
to enhance the current descriptions.

Many of the most exciting questions in astrophysics and cosmology, including
the majority of observational probes of dark energy, rely on an understanding
of the nonlinear regime of structure formation. In order to fully exploit the
information available from this regime and to extract cosmological constraints,
accurate theoretical predictions are needed. Currently such predictions can
only be obtained from costly, precision numerical simulations. This paper is
the third in a series aimed at constructing an accurate calibration of the
nonlinear mass power spectrum on Mpc scales for a wide range of currently
viable cosmological models, including dark energy. The first two papers
addressed the numerical challenges, and the scheme by which an interpolator was
built from a carefully chosen set of cosmological models. In this paper we
introduce the "Coyote Univers"' simulation suite which comprises nearly 1,000
Nbody simulations at different force and mass resolutions, spanning 38 wCDM
cosmologies. This large simulation suite enables us to construct a prediction
scheme, or emulator, for the nonlinear matter power spectrum accurate at the
percent level out to k~1 h/Mpc. We describe the construction of the emulator,
explain the tests performed to ensure its accuracy, and discuss how the central
ideas may be extended to a wider range of cosmological models and applications.
A power spectrum emulator code is released publicly as part of this paper.