
The decay of gravitational potentials in the presence of dark energy leads to
an additional, latetime contribution to anisotropies in the cosmic microwave
background (CMB) at large angular scales. The imprint of this socalled
integrated SachsWolfe (ISW) effect to the CMB angular power spectrum has been
detected and studied in detail, but reconstructing its spatial contributions to
the CMB $\textit{map}$, which would offer the tantalizing possibility of
separating the early from the latetime contributions to CMB temperature
fluctuations, is more challenging. Here we study the technique for
reconstructing the ISW map based on information from galaxy surveys and focus
in particular on how its accuracy is impacted by the presence of photometric
calibration errors in input galaxy maps, which were previously found to be a
dominant contaminant for ISW signal estimation. We find that both including
tomographic information from a single survey and using data from multiple,
complementary galaxy surveys improve the reconstruction by effectively
selfcalibrating the estimator against spurious power contributions from
calibration errors. A highfidelity reconstruction further requires one to
account for the contribution of calibration errors to the observed galaxy power
spectrum in the model used to construct the ISW estimator. We find that if the
photometric calibration errors in galaxy surveys can be independently
controlled at the level required to obtain unbiased dark energy constraints,
then it is possible to reconstruct ISW maps with excellent accuracy using a
combination of maps from two galaxy surveys with properties similar to Euclid
and SPHEREx.

The discovery of the accelerating universe in the late 1990s was a watershed
moment in modern cosmology, as it indicated the presence of a fundamentally
new, dominant contribution to the energy budget of the universe. Evidence for
dark energy, the new component that causes the acceleration, has since become
extremely strong, owing to an impressive variety of increasingly precise
measurements of the expansion history and the growth of structure in the
universe. Still, one of the central challenges of modern cosmology is to shed
light on the physical mechanism behind the accelerating universe. In this
review, we briefly summarize the developments that led to the discovery of dark
energy. Next, we discuss the parametric descriptions of dark energy and the
cosmological tests that allow us to better understand its nature. We then
review the cosmological probes of dark energy. For each probe, we briefly
discuss the physics behind it and its prospects for measuring dark energy
properties. We end with a summary of the current status of dark energy
research.

The current $>3\sigma$ tension between the Hubble constant $H_0$ measured
from local distance indicators and from cosmic microwave background is one of
the most highly debated issues in cosmology, as it possibly indicates new
physics or unknown systematics. In this work, we explore whether this tension
can be alleviated by the sample variance in the local measurements, which use a
small fraction of the Hubble volume. We use a largevolume cosmological
$N$body simulation to model the local measurements and to quantify the
variance due to local density fluctuations and sample selection. We explicitly
take into account the inhomogeneous spatial distribution of type Ia supernovae.
Despite the faithful modelling of the observations, our results confirm
previous findings that sample variance in the local Hubble constant $(H_0^{\rm
loc})$ measurements is small; we find $\sigma(H_0^{\rm loc})=0.31\,{\rm km\
s^{1}Mpc^{1}}$, a nearly negligible fraction of the $\sim6\,{\rm km\
s^{1}Mpc^{1}}$ necessary to explain the difference between the local and the
global $H_0$ measurements. While the $H_0$ tension could in principle be
explained by our local neighbourhood being a underdense region of radius $\sim
150 \,\rm Mpc$ , the extreme required underdensity of such a void
$(\delta\simeq 0.8)$ makes it very unlikely in a $\Lambda$CDM universe, and it
also violates existing observational constraints. Therefore, sample variance in
a $\Lambda$CDM universe cannot appreciably alleviate the tension in $H_0$
measurements even after taking into account the inhomogeneous selection of type
Ia supernovae.

We present a novel approach to constrain accelerating cosmologies with galaxy
cluster phase spaces. With the Fisher matrix formalism we forecast constraints
on the cosmological parameters that describe the cosmological expansion
history. We find that our probe has the potential of providing constraints
comparable to, or even stronger than, those from other cosmological probes.
More specifically, with 1000 (100) clusters uniformly distributed in the
redshift range $ 0 \leq z \leq 0.8$, after applying a conservative $80\%$ mass
scatter prior on each cluster and marginalizing over all other parameters, we
forecast $1\sigma$ constraints on the dark energy equation of state $w$ and
matter density parameter $\Omega_M$ of $\sigma_w = 0.138 (0.431)$ and
$\sigma_{\Omega_M} = 0.007 (0.025)$ in a flat universe. Assuming 40\% mass
scatter and adding a prior on the Hubble constant we can achieve a constraint
on the CPL parametrization of the dark energy equation of state parameters
$w_0$ and $w_a$ with 100 clusters in the same redshift range: $\sigma_{w_0} =
0.191 $ and $\sigma_{w_a} = 2.712$. Dropping the assumption of flatness and
assuming $w=1$ we also attain competitive constraints on the matter and dark
energy density parameters: $\sigma_{\Omega_M} = 0.101$ and
$\sigma_{\Omega_{\Lambda}} = 0.197$ for 100 clusters uniformly distributed in
the range $ 0 \leq z \leq 0.8$ after applying a prior on the Hubble constant.
We also discuss various observational strategies for tightening constraints in
both the near and far future.

Peculiar velocities of objects in the nearby universe are correlated due to
the gravitational pull of largescale structure. By measuring these velocities,
we have a unique opportunity to test the cosmological model at the lowest
redshifts. We perform this test, using current data to constrain the amplitude
of the "signal" covariance matrix describing the velocities and their
correlations. We consider a new, wellcalibrated "Supercal" set of lowredshift
SNe Ia as well as a set of distances derived from the fundamental plane
relation of 6dFGS galaxies. Analyzing the SN and galaxy data separately, both
results are consistent with the peculiar velocity signal of our fiducial
$\Lambda$CDM model, ruling out the noiseonly model with zero peculiar
velocities at greater than $7\sigma$ (SNe) and $8\sigma$ (galaxies). When the
two data sets are combined appropriately, the precision of the test increases
slightly, resulting in a constraint on the signal amplitude of $A =
1.05_{0.21}^{+0.25}$, where $A = 1$ corresponds to our fiducial model.
Equivalently, we report an 11% measurement of the product of the growth rate
and amplitude of mass fluctuations evaluated at $z_\text{eff} = 0.02$, $f
\sigma_8 = 0.428_{0.045}^{+0.048}$, valid for our fiducial $\Lambda$CDM model.
We explore the robustness of the results to a number of conceivable variations
in the analysis and find that individual variations shift the preferred signal
amplitude by less than ${\sim}0.5\sigma$. We briefly discuss our Supercal SN Ia
results in comparison with our previous results using the JLA compilation.

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.

The Integrated SachsWolfe (ISW) effect is a largeangle modulation of the
cosmic microwave background (CMB), generated when CMB photons traverse evolving
potential wells associated with large scale structure (LSS). Recent efforts
have been made to reconstruct maps of the ISW signal using information from
surveys of galaxies and other LSS tracers, but investigation into how survey
systematics affect their reliability has so far been limited. Using simulated
ISW and LSS maps, we study the impact of galaxy survey properties and
systematic errors on the accuracy of reconstructed ISW signal. We find that
systematics that affect the observed distribution of galaxies along the line of
sight, such as photoz and biasevolution related errors, have a relatively
minor impact on reconstruction quality. In contrast, however, we find that
directiondependent calibration errors can be very harmful. Specifically, we
find that in order to avoid significant degradation of our reconstruction
quality statistics, directiondependent number density fluctuations due to
systematics must be controlled so that their variance is smaller than $10^{6}$
(which corresponds to a 0.1% calibration). Additionally, we explore the
implications of our results for attempts to use reconstructed ISW maps to shed
light on the origin of largeangle CMB alignments. We find that there is only a
weak correlation between the true and reconstructed angular momentum
dispersion, which quantifies alignment, even for reconstructed ISW maps which
are fairly accurate overall.

Strong lensing timedelay systems constrain cosmological parameters via the
socalled timedelay distance and the angular diameter distance to the lens. In
previous studies, only the former information was used. In this paper, we show
that the cosmological constraints improve significantly when the latter
information is also included. Specifically, the angular diameter distance plays
a crucial role in breaking the degeneracy between the curvature of the Universe
and the timevarying equation of state of dark energy. Using a mock sample of
55 bright quadruple lens systems based on expectations for ongoing/future
imaging surveys, we find that adding the angular diameter distance information
to the timedelay distance information and the cosmic microwave background data
of Planck improves the constraint on the constant equation of state by 30%, on
the time variation in the equation of state by a factor of two, and on the
Hubble constant in the flat $\Lambda$CDM model by a factor of two. Therefore,
previous forecasts for the statistical power of timedelay systems were
significantly underestimated, i.e., timedelay systems are more powerful than
previously appreciated.

We revisit the effect of peculiar velocities on lowredshift type Ia
supernovae. Velocities introduce an additional guaranteed source of
correlations between supernova magnitudes that should be considered in all
analyses of nearby supernova samples but has largely been neglected in the
past. Applying a likelihood analysis to the latest compilation of nearby
supernovae, we find no evidence for the presence of these correlations,
although, given the significant noise, the data is also consistent with the
correlations predicted for the standard $\Lambda$CDM model. We then consider
the dipolar component of the velocity correlations  the frequently studied
"bulk velocity"  and explicitly demonstrate that including the velocity
correlations in the data covariance matrix is crucial for drawing correct and
unambiguous conclusions about the bulk flow. In particular, current supernova
data is consistent with no excess bulk flow on top of what is expected in
$\Lambda$CDM and effectively captured by the covariance. We further clarify the
nature of the apparent bulk flow that is inferred when the velocity covariance
is ignored. We show that a significant fraction of this quantity is expected to
be noise bias due to uncertainties in supernova magnitudes and not any physical
peculiar motion.

Upcoming or future deep galaxy samples with wide sky coverage can provide
independent measurement of the kinematic dipole  our motion relative to the
rest frame defined by the largescale structure. Such a measurement would
present an important test of the standard cosmological model, as the standard
model predicts the galaxy measurement should precisely agree with the existing
precise measurements made using the CMB. However, the required statistical
precision to measure the kinematic dipole typically makes the measurement
susceptible to bias from the presence of the localstructureinduced dipole
contamination. In order to minimize the latter, a sufficiently deep survey is
required. We forecast both the statistical error and the systematic bias in the
kinematic dipole measurements. We find that a survey covering $\sim 75\%$ of
the sky in both hemispheres and having $\sim 30$ million galaxies can detect
the kinematic dipole at $5\sigma$, while its median redshift should be at least
$z_{med} \sim 0.75$ for negligible bias from the local structure.

Several unexpected features have been observed in the microwave sky at large
angular scales, both by WMAP an by Planck. Among those features is a lack of
both variance and correlation on the largest angular scales, alignment of the
lowest multipole moments with one another and with the motion and geometry of
the Solar System, a hemispherical power asymmetry or dipolar power modulation,
a preference for odd parity modes and an unexpectedly large cold spot in the
Southern hemisphere. The individual pvalues of the significance of these
features are in the per mille to per cent level, when compared to the
expectations of the bestfit inflationary $\Lambda$CDM model. Some pairs of
those features are demonstrably uncorrelated, increasing their combined
statistical significance and indicating a significant detection of CMB features
at angular scales larger than a few degrees on top of the standard model.
Despite numerous detailed investigations, we still lack a clear understanding
of these largescale features, which seem to imply a violation of statistical
isotropy and scale invariance of inflationary perturbations. In this
contribution we present a critical analysis of our current understanding and
discuss several ideas of how to make further progress.

We provide software with a graphical user interface to calculate the
phenomenology of a wide class of dark energy models featuring multiple scalar
fields. The user chooses a subclass of models and, if desired, initial
conditions, or else a range of initial parameters for Monte Carlo. The code
calculates the energy density of components in the universe, the equation of
state of dark energy, and the linear growth of density perturbations, all as a
function of redshift and scale factor. The output also includes an approximate
conversion into the average equation of state, as well as the common $(w_0,
w_a)$ parametrization. The code is available here:
http://github.com/kahinton/DarkEnergyUIandMC

We revisit the alignments of the largest structures observed in the cosmic
microwave background (CMB) using the seven and nineyear WMAP and firstyear
Planck data releases. The observed alignments  the quadrupole with the
octopole and their joint alignment with the direction of our motion with
respect to the CMB (the dipole direction) and the geometry of the Solar System
(defined by the Ecliptic plane)  are generally in good agreement with results
from the previous WMAP data releases. However, a closer look at fullsky data
on the largest scales reveals discrepancies between the earlier WMAP data
releases (three to sevenyear) and the final nineyear release. There are also
discrepancies between all the WMAP data releases and the firstyear Planck
release. Nevertheless, both the WMAP and Planck data confirm the alignments of
the largest observable CMB modes in the Universe. In particular, the pvalues
for the mutual alignment between the quadrupole and octopole, and the alignment
of the plane defined by the two with the dipole direction, are both at the
greater than 3sigma level for all three Planck maps studied. We also calculate
conditional statistics on the various alignments and find that it is currently
difficult to unambiguously identify a leading anomaly that causes the others or
even to distinguish correlation from causation.

The lack of largeangle correlations in the observed microwave background
temperature fluctuations persists in the finalyear maps from WMAP and the
first cosmological data release from Planck. We find a statistically robust and
significant result: pvalues for the missing correlations lying below 0.24 per
cent (i.e. evidence at more than 3 sigma) for foreground cleaned maps, in
complete agreement with previous analyses based upon earlier WMAP data. A
conservative cutsky analysis of the Planck HFI 100 GHz frequency band, the
`cleanest CMB channel', returns a pvalue as small as 0.03 per cent, based on
the conservative mask defined by WMAP. These findings are in stark contrast to
expectations from the inflationary Lambda cold dark matter model and still lack
a convincing explanation. If this lack of largeangle correlations is a true
feature of our Universe, and not just a statistical fluke, then the
cosmological dipole must be considerably smaller than that predicted in the
bestfitting model.

We perform parametric tests of the consistency of the standard $w$CDM model
in the framework of general relativity by carefully separating information
between the geometry and growth of structure. We replace each lateuniverse
parameter that describes the behavior of dark energy with two parameters: one
describing geometrical information in cosmological probes, and the other
controlling the growth of structure. We use data from all principal
cosmological probes: of these, Type Ia supernovae, baryon acoustic
oscillations, and the peak locations in the cosmic microwave background angular
power spectrum constrain the geometry, while the redshift space distortions,
weak gravitational lensing and the abundance of galaxy clusters constrain both
geometry and growth. Both geometry and growth separately favor the $\Lambda$CDM
cosmology with the matter density relative to critical $\Omega_M\simeq 0.3$.
When the equation of state is allowed to vary separately for probes of growth
and geometry, we find again a good agreement with the $\Lambda$CDM value
($w\simeq 1$), with the major exception of redshiftspace distortions which
favor less growth than in $\Lambda$CDM at 3$\sigma$ confidence, favoring the
equation of state $w^{\rm grow}\simeq 0.8$. The anomalous growth favored by
redshift space distortions has been noted earlier, and is common to all
redshift space distortion data sets, but may well be caused by systematics, or
be explained by the sum of the neutrino masses higher than that expected from
the simplest mass hierarchies, $m_\nu \simeq 0.45$ eV. On the whole, the
constraints are tight even in the new, larger parameter space due to impressive
complementarity of different cosmological probes.

Two of the most commonly used tools to constrain the primordial
nonGaussianity are the bispectrum and the Minkowski functionals of CMB
temperature anisotropies. These two measures of nonGaussianity in principle
provide distinct (though correlated) information, but in the past constraints
from them have only been loosely compared, and not statistically combined. In
this work we evaluate, for the first time, the covariance matrix between the
local nonGaussianity coefficient fnl estimated through the bispectrum and
Minkowski functionals. We find that the estimators are positively correlated,
with corerlation coefficient r ~ 0.3. Using the WMAP7 data to combine the two
measures and accounting for the pointsource systematics, we find the combined
constraint fnl=37+/28, which has a ~20% smaller error than either of the
individual constraints.

We develop a general method to "selfcalibrate" observations of galaxy
clustering with respect to systematics associated with photometric calibration
errors. We first point out the danger posed by the multiplicative effect of
calibration errors, where largeangle error propagates to small scales and may
be significant even if the largescale information is cleaned or not used in
the cosmological analysis. We then propose a method to measure the arbitrary
largescale calibration errors and use these measurements to correct the
smallscale (highmultipole) power which is most useful for constraining the
majority of cosmological parameters. We demonstrate the effectiveness of our
approach on synthetic examples and briefly discuss how it may be applied to
real data.

Continuum and HI surveys with the Square Kilometre Array (SKA) will allow us
to probe some of the most fundamental assumptions of modern cosmology,
including the Cosmological Principle. SKA allsky surveys will map an enormous
slice of spacetime and reveal cosmology at superhorizon scales and redshifts
of order unity. We illustrate the potential of these surveys and discuss the
prospects to measure the cosmic radio dipole at high fidelity. We outline
several potentially transformational tests of cosmology to be carried out by
means of SKA allsky surveys.

The statistics of primordial curvature fluctuations are our window into the
period of inflation, where these fluctuations were generated. To date, the
cosmic microwave background has been the dominant source of information about
these perturbations. Large scale structure is however from where drastic
improvements should originate. In this paper, we explain the theoretical
motivations for pursuing such measurements and the challenges that lie ahead.
In particular, we discuss and identify theoretical targets regarding the
measurement of primordial nonGaussianity. We argue that when quantified in
terms of the local (equilateral) template amplitude $f_{\rm NL}^{\rm loc}$
($f_{\rm NL}^{\rm eq}$), natural target levels of sensitivity are $\Delta
f_{\rm NL}^{\rm loc, eq.} \simeq 1$. We highlight that such levels are within
reach of future surveys by measuring 2, 3 and 4point statistics of the
galaxy spatial distribution. This paper summarizes a workshop held at CITA
(University of Toronto) on October 2324, 2014.

We investigate how the use of higher wavenumbers (smaller scales) in the
galaxy clustering power spectrum influences cosmological constraints. We take
into account uncertainties from nonlinear density fluctuations, (scale
dependent) galaxy bias, and baryonic effects. Allowing for substantially model
independent uncertainties through separate fit parameters in each wavenumber
bin that also allow for the redshift evolution, we quantify strong gains in
dark energy and neutrino mass leverage with increasing maximum wavenumber,
despite marginalizing over numerous (up to 125) extra fit parameters. The
leverage is due to not only an increased number of modes but, more
significantly, breaking of degeneracies beyond the linear regime.

We test the statistical isotropy of the universe by analyzing the
distribution of WISE extragalactic sources that were also observed by 2MASS. We
pay particular attention to color cuts and foreground marginalization in order
to cull a uniform sample of extragalactic objects and avoid stars. We detect a
dipole gradient in the numbercounts with an amplitude of $\sim$0.05, somewhat
larger than expectations based on local structures corresponding to the depth
and (independently measured) bias of our WISE2MASS sources. The direction of
the dipole, $(l,b)\simeq (310\,^{\circ}, 15\,^{\circ})$, is in reasonably good
agreement with that found previously in the (shallower) 2MASS Extended Source
Catalog alone. Interestingly, the dipole direction is not far from the
direction of the dipolar modulation in the CMB found by Planck, and also fairly
closely matches largescalestructure bulkflow directions found by various
groups using galaxies and type Ia supernovae. It is difficult, however, to draw
specific conclusions from the nearagreement of these directions.

Some recent observations provide $> 2\sigma$ evidence for phantom dark energy
 a value of the dark energy equation of state less than the
cosmologicalconstant value of $1$. We focus on constraining the equation of
state by combining current data from the most mature geometrical probes of dark
energy: Type Ia supernovae (SNe Ia) from the Supernova Legacy Survey (SNLS3),
the Supernova Cosmology Project (Union2.1), and the PanSTARRS1 survey (PS1);
cosmic microwave background measurements from Planck and WMAP9; and a
combination of measurements of baryon acoustic oscillations. The combined data
are consistent with $w = 1$ for the Union2.1 sample, though they present
moderate ($\sim 1.9\sigma$) evidence for a phantom value when either the SNLS3
or PS1 sample is used instead. We study the dependence of the constraints on
the redshift, stretch, color, and host galaxy stellar mass of SNe, but we find
no unusual trends. In contrast, the constraints strongly depend on any external
$H_0$ prior: a higher adopted value for the direct measurement of the Hubble
constant ($H_0 \gtrsim 71~\text{km/s/Mpc}$) leads to $\gtrsim 2\sigma$ evidence
for phantom dark energy. Given Planck data, we can therefore make the following
statement at $2\sigma$ confidence: either the SNLS3 and PS1 data have
systematics that remain unaccounted for or the Hubble constant is below 71
km/s/Mpc; else the dark energy equation of state is indeed phantom.

The massive dark matter halos that host groups and clusters of galaxies have
observable properties that appear to be lognormally distributed about
powerlaw mean scaling relations in halo mass. Coupling this assumption with
either quadratic or cubic approximations to the mass function in log space, we
derive closedform expressions for the space density of halos as a function of
multiple observables as well as forms for the loworder moments of properties
of observableselected samples. Using a Tinker mass function in a {\Lambda}CDM
cosmology, we show that the cubic analytic model reproduces results obtained
from direct, numerical convolution at the 10 percent level or better over
nearly the full range of observables covered by current observations and for
redshifts extending to z = 1.5. The model provides an efficient framework for
estimating effects arising from selection and covariance among observable
properties in survey samples.

The quantity and quality of cosmic structure observations have greatly
accelerated in recent years. Further leaps forward will be facilitated by
imminent projects, which will enable us to map the evolution of dark and
baryonic matter density fluctuations over cosmic history. The way that these
fluctuations vary over space and time is sensitive to the nature of dark matter
and dark energy. Dark energy and gravity both affect how rapidly structure
grows; the greater the acceleration, the more suppressed the growth of
structure, while the greater the gravity, the more enhanced the growth. While
distance measurements also constrain dark energy, the comparison of growth and
distance data tests whether General Relativity describes the laws of physics
accurately on large scales. Modified gravity models are able to reproduce the
distance measurements but at the cost of altering the growth of structure
(these signatures are described in more detail in the accompanying paper on
Novel Probes of Gravity and Dark Energy). Upcoming surveys will exploit these
differences to determine whether the acceleration of the Universe is due to
dark energy or to modified gravity. To realize this potential, both wide field
imaging and spectroscopic redshift surveys play crucial roles. Projects
including DES, eBOSS, DESI, PFS, LSST, Euclid, and WFIRST are in line to map
more than a 1000 cubicbillionlightyear volume of the Universe. These will
map the cosmic structure growth rate to 1% in the redshift range 0<z<2, over
the last 3/4 of the age of the Universe.