
We derive the evolution equation for the density matrix of a UV and IR
limited band of comoving momentum modes of the canonically normalized scalar
degree of freedom in two examples of nearly de Sitter universes. Including the
effects of a cubic interaction term from the gravitational action and tracing
out a set of longer wavelength modes, we find that the evolution of the system
is nonHamiltonian and nonMarkovian. We find linear dissipation terms for a
few modes with wavelength near the boundary between system and bath and
nonlinear dissipation terms for all modes. The nonHamiltonian terms persist to
late times when the scalar field dynamics is such that the curvature
perturbation continues to evolve on superHubble scales.

We study statistical anisotropies generated in the observed twopoint
function of the cosmic microwave background (CMB) fluctuations if the
primordial statistics are nonGaussian. Focusing on the dipole modulations of
the anisotropies, we find that the hemispherical power asymmetry observed in
the CMB temperature fluctuations can be modeled by a localtype trispectrum
with amplitude $\tau_{\rm NL}(k_p=0.05~{\rm Mpc}^{1}) \approx 2 \times 10^4$
and a large red tilt $n\approx 0.68$. We numerically evaluate the nonGaussian
covariance of the modulation estimators for both temperature and Emode
polarization fluctuations and discuss the prospects of constraining the model
using Planck satellite data. We then discuss other effects of the
scaledependent trispectrum that could be used to distinguish this scenario
from other explanations of the power asymmetry: higherorder modulations of the
twopoint function and the nonGaussian angular power spectrum covariance. As
an important consequence of the nonGaussian power spectrum covariance, we
discuss how the CMBinferred spectral index of primordial scalar fluctuations
can be significantly biased in the presence of a scaledependent localtype
trispectrum.

We explore the possible spectrum of binary mergers of subsolar mass black
holes formed out of dark matter particles interacting via a dark
electromagnetism. We estimate the properties of these dark black holes by
assuming that their formation process is parallel to PopulationIII star
formation; except that dark molecular cooling can yield smaller opacity limit.
We estimate the binary coalescence rates for the Advanced LIGO and Einstein
telescope, and find that scenarios compatible with all current constraints
could produce dark black holes at rates high enough for detection by Advanced
LIGO.

We elucidate the counting of the relevant small parameters in inflationary
perturbation theory. Doing this allows for an explicit delineation of the
domain of validity of the semiclassical approximation to gravity used in the
calculation of inflationary correlation functions. We derive an expression for
the dependence of correlation functions of inflationary perturbations on the
slowroll parameter $\epsilon = \dot{H}/H^2$, as well as on $H/M_p$, where $H$
is the Hubble parameter during inflation. Our analysis is valid for
singlefield models in which the inflaton can traverse a Plancksized range in
field values and where all slowroll parameters have approximately the same
magnitude. As an application, we use our expression to seek the boundaries of
the domain of validity of inflationary perturbation theory for regimes where
this is potentially problematic: models with small speed of sound and models
allowing eternal inflation.

We develop and study the positiondependent bispectrum. It is a
generalization of the recently proposed positiondependent power spectrum
method of measuring the squeezedlimit bispectrum. The positiondependent
bispectrum can similarly be used to measure the squeezedlimit trispectrum in
which one of the wavelengths is much longer than the other three. In this work,
we will mainly consider the case in which the three smaller wavelengths are
nearly the same (the equilateral configuration). We use the Fisher information
matrix to forecast constraints on bias parameters and the amplitude of
primordial trispectra from the positiondependent bispectrum method. We find
that the method can constrain the localtype $g_{\rm NL}$ at a level of
$\sigma(g_{\rm NL}^{\rm local}) \approx 3 \times 10^5$ for a large volume
SPHERExlike survey; improvements can be expected by including all the
triangular configurations of the bispectra rather than just the equilateral
configuration. However, the same measurement would also constrain a much larger
family of trispectra than local $g_{\rm NL}$ model. We discuss the implications
of the forecasted reach of future surveys in terms of super cosmic variance
uncertainties from primordial nonGaussianities.

This book lays out the scientific goals to be addressed by the
nextgeneration groundbased cosmic microwave background experiment, CMBS4,
envisioned to consist of dedicated telescopes at the South Pole, the high
Chilean Atacama plateau and possibly a northern hemisphere site, all equipped
with new superconducting cameras. CMBS4 will dramatically advance cosmological
studies by crossing critical thresholds in the search for the Bmode
polarization signature of primordial gravitational waves, in the determination
of the number and masses of the neutrinos, in the search for evidence of new
light relics, in constraining the nature of dark energy, and in testing general
relativity on large scales.

We develop a mathematical construction of nonGaussian fields whose bispectra
satisfy the singleclock inflation consistency relation. At the same order that
our basis for bispectra recovers the two simplest single clock templates, we
also find a third orthogonal template which has the single clock squeezed
limit, peaks in folded configurations, and has very small coupling in the
equilateral limit. We explore the map between templates and operators in a very
general Lagrangian for singleclock fluctuations and find no significant
overlap between the new template and models in the literature. We comment on
the physical implications of this conclusion. Our findings add support for the
idea that both theory and data driven considerations will be best served if
next generation nonGaussianity constraints are made in a basis that uses the
degree of coupling between long and short wavelength modes as an organizing
principle.

We examine the squeezed limit of the bispectrum when a light scalar with
arbitrary nonderivative selfinteractions is coupled to the inflaton. We find
that when the hidden sector scalar is sufficiently light ($m\lesssim0.1\,H$),
the coupling between long and short wavelength modes from the series of higher
order correlation functions (from arbitrary order contact diagrams) causes the
statistics of the fluctuations to vary in subvolumes. This means that
observations of primordial nonGaussianity cannot be used to uniquely
reconstruct the potential of the hidden field. However, the local bispectrum
induced by modecoupling from these diagrams always has the same squeezed
limit, so the field's locally determined mass is not affected by this cosmic
variance.

We present a science forecast for the eBOSS survey, part of the SDSSIV
project, which is a spectroscopic survey using multiple tracers of largescale
structure, including luminous red galaxies (LRGs), emission line galaxies
(ELGs) and quasars (both as a direct probe of structure and through the
Ly$\alpha$ forest). Focusing on discrete tracers, we forecast the expected
accuracy of the baryonic acoustic oscillation (BAO), the redshiftspace
distortion (RSD) measurements, the $f_{\rm NL}$ parameter quantifying the
primordial nonGaussianity, the dark energy and modified gravity parameters. We
also use the lineofsight clustering in the Ly$\alpha$ forest to constrain
the total neutrino mass. We find that eBOSS LRGs ($0.6<z<1.0$) (combined with
the BOSS LRGs at $z>0.6$), ELGs ($0.6<z<1.2$) and Clustering Quasars (CQs)
($0.6<z<2.2$) can achieve a precision of 1%, 2.2% and 1.6% precisions,
respectively, for spherically averaged BAO distance measurements. Using the
same samples, the constraint on $f\sigma_8$ is expected to be 2.5%, 3.3% and
2.8% respectively. For primordial nonGaussianity, eBOSS alone can reach an
accuracy of $\sigma(f_{\rm NL})\sim1015$, depending on the external
measurement of the galaxy bias and our ability to model largescale systematic
errors. eBOSS can at most improve the dark energy Figure of Merit (FoM) by a
factor of $3$ for the ChevallierPolarskiLinder (CPL) parametrisation, and can
well constrain three eigenmodes for the general equationofstate parameter
(Abridged).

We derive a general expression for the probability of observing deviations
from statistical isotropy in the cosmic microwave background (CMB) if the
primordial fluctuations are nonGaussian and extend to superhorizon scales. The
primary motivation is to properly characterize the monopole and dipole
modulations of the primordial power spectrum that are generated by the coupling
between superhorizon and subhorizon perturbations. Unlike previous proposals
for generating the hemispherical power asymmetry, we do not assume that the
power asymmetry results from a single large superhorizon mode. Instead, we
extrapolate the observed power spectrum to superhorizon scales and compute the
power asymmetry that would result from a specific realization of nonGaussian
perturbations on scales larger than the observable universe. Our study
encompasses many of the scenarios that have been put forward as possible
explanations for the CMB hemispherical power asymmetry. We confirm our analytic
predictions for the probability of a given power asymmetry by comparing them to
numerical realizations of CMB maps. We find that nonlocal models of
nonGaussianity and scaledependent local nonGaussianity produce
scaledependent modulations of the power spectrum, thereby potentially
producing both a monopolar and a dipolar power modulation on large scales. We
then provide simple examples of finding the posterior distributions for the
parameters of the bispectrum from the observed monopole and dipole modulations.

The Extended Baryon Oscillation Spectroscopic Survey (eBOSS) will conduct
novel cosmological observations using the BOSS spectrograph at Apache Point
Observatory. Observations will be simultaneous with the Time Domain
Spectroscopic Survey (TDSS) designed for variability studies and the
Spectroscopic Identification of eROSITA Sources (SPIDERS) program designed for
studies of Xray sources. eBOSS will use four different tracers to measure the
distanceredshift relation with baryon acoustic oscillations (BAO). Using more
than 250,000 new, spectroscopically confirmed luminous red galaxies at a median
redshift z=0.72, we project that eBOSS will yield measurements of $d_A(z)$ to
an accuracy of 1.2% and measurements of H(z) to 2.1% when combined with the
z>0.6 sample of BOSS galaxies. With ~195,000 new emission line galaxy
redshifts, we expect BAO measurements of $d_A(z)$ to an accuracy of 3.1% and
H(z) to 4.7% at an effective redshift of z= 0.87. A sample of more than 500,000
spectroscopicallyconfirmed quasars will provide the first BAO distance
measurements over the redshift range 0.9<z<2.2, with expected precision of 2.8%
and 4.2% on $d_A(z)$ and H(z), respectively. Finally, with 60,000 new quasars
and reobservation of 60,000 quasars known from BOSS, we will obtain new
Lymanalpha forest measurements at redshifts z>2.1; these new data will enhance
the precision of $d_A(z)$ and H(z) by a factor of 1.44 relative to BOSS.
Furthermore, eBOSS will provide improved tests of General Relativity on
cosmological scales through redshiftspace distortion measurements, improved
tests for nonGaussianity in the primordial density field, and new constraints
on the summed mass of all neutrino species. Here, we provide an overview of the
cosmological goals, spectroscopic target sample, demonstration of spectral
quality from early data, and projected cosmological constraints from eBOSS.

We present examples of nonGaussian statistics that can induce bispectra
matching local and nonlocal (including equilateral) templates in biased
subvolumes. We find cases where the biasing from coupling to long wavelength
modes affects only the power spectrum, only the bispectrum or both. Our results
suggest that ruling out multifield scenarios is quite difficult: some
measurements of the scalar correlations, including the shape of the bispectrum,
can be consistent with singleclock predictions even when cosmic variance from
superhorizon modes is at work. Furthermore, if observations of the density
perturbations rule out singleclock inflation, we will face a serious cosmic
variance obstacle in drawing any further conclusions about the particle physics
origin of the scalar fluctuations.

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 perform cosmological Nbody simulations with nonGaussian initial
conditions generated from two independent fields. The dominant contribution to
the perturbations comes from a purely Gaussian field, but we allow the second
field to have local nonGaussianity that need not be weak. This scenario allows
us to adjust the relative importance of nonGaussian contributions beyond the
skewness, producing a scaling of the higher moments different from (and
stronger than) the scaling in the usual single field local ansatz. We compare
semianalytic prescriptions for the nonGaussian mass function, large scale
halo bias, and stochastic bias against the simulation results. We discuss
applications of this work to large scale structure measurements that can test a
wider range of models for the primordial fluctuations than is usually explored.

If long wavelength primordial tensor modes are coupled to short wavelength
scalar modes, the scalar curvature twopoint function will have an offdiagonal
component. This `fossil' remnant is a signature of a mode coupling that cannot
be achieved in single clock inflation. Any constraint on its presence allows a
cross check of the relationship between the dynamical generation of the
fluctuations and the evolution of the inflationary background. We use the
example of nonBunch Davies initial states for the tensor and scalar modes to
demonstrate that physically reasonable fossils, consistent with current data,
can be observable in the near future. We illustrate how the fossil offdiagonal
power spectrum is a complementary probe to the squeezed limit bispectra of the
scalar and tensor sectors individually. We also quantify the relation between
the observable signal and the squeezed limit bispectrum for a general
scalarscalarfossil coupling, and note the effect of superhorizon tensor modes
on the anisotropy in scalar modes.

Primordial nonGaussianity induces a scaledependent bias in large scale
structure (LSS) data, proportional to $f_{\rm NL}/k^2$ for the exact local
ansatz. Recent work has shown that models of inflation that predict a large
squeezed limit bispectrum, such as multifield models and single field
inflation with a modified initial state, typically give rise to a generalized
local ansatz, with the scaledependent bias now proportional to ${\cal A}_{\rm
NL}/k^\alpha$. We use photometric measurements of the angular power spectrum of
luminous red galaxies and quasars in the Sloan Digital Sky Survey Data Release
Eight (SDSS DR8) with the above parameterization to constrain the amplitude
${\cal A}_{\rm NL}$ and scaledependence $\alpha$. We find that the
marginalized upper limit on $\alpha$ is 2.0 at the $95\%$ confidence level,
consistent with the local ansatz. We also present Fisher forecasts for a survey
of the same size as DR8 to assess the role of systematics in current
photometric LSS data. Moreover, we present analytic results on the expected
mass dependence of ${\cal A}_{\rm NL}$ for different inflationary models, which
can be an important observable for future surveys, if primordial
nonGaussianity is nonzero.

We describe a general procedure for using number counts of any object to
constrain the probability distribution of the primordial fluctuations, allowing
for generic weak nonGaussianity. We apply this procedure to use limits on the
abundance of primordial black holes and dark matter ultracompact minihalos
(UCMHs) to characterize the allowed statistics of primordial fluctuations on
very small scales. We present constraints on the power spectrum and the
amplitude of the skewness for two different families of nonGaussian
distributions, distinguished by the relative importance of higher moments.
Although primordial black holes probe the smallest scales, ultracompact
minihalos provide significantly stronger constraints on the power spectrum and
so are more likely to eventually provide smallscale constraints on
nonGaussianity.

We demonstrate that local, scaledependent nonGaussianity can generate
cosmic variance uncertainty in the observed spectral index of primordial
curvature perturbations. In a universe much larger than our current Hubble
volume, locally unobservable long wavelength modes can induce a
scaledependence in the power spectrum of typical subvolumes, so that the
observed spectral index varies at a cosmologically significant level (\Delta
ns  ~ O(0.04)). Similarly, we show that the observed bispectrum can have an
induced scale dependence that varies about the global shape. If tensor modes
are coupled to long wavelength modes of a second field, the locally observed
tensor power and spectral index can also vary. All of these effects, which can
be introduced in models where the observed nonGaussianity is consistent with
bounds from the Planck satellite, loosen the constraints that observations
place on the parameters of theories of inflation with mode coupling. We suggest
observational constraints that future measurements could aim for to close this
window of cosmic variance uncertainty.

We study the behavior of npoint functions of the primordial curvature
perturbations, assuming our observed Universe is only a subset of a larger
space with statistically homogeneous and isotropic perturbations. If the larger
space has arbitrary npoint functions in a family of local type nonGaussian
statistics, sufficiently biased smaller volumes will have statistics from a
`natural' version of that family with moments that are weakly nonGaussian and
ordered, regardless of the statistics of the original field. We also describe
the effect of this bias on the shape of the bispectrum.

Thanks to the rapid progress in precision cosmology in the last few years, we
now have access to physical observables that may constrain the theory of
inflation through the nonGaussianity (NG) signatures in the cosmic microwave
background radiation and the distribution of largescale structure. Numerical
modeling of the NG signals from different inflation models is essential to
correctly interpret current and near future data from largescale structure
surveys. In this study, we use highresolution cosmological hydrodynamical
simulations to investigate the effects of primordial NG on the formation and
evolution of galaxies from the cosmic dawn to the present day. Focusing on the
local type primordial NG, we find that it may affect the formation history of
stars and black holes in galaxies, and their distribution. Compared to the
Gaussian case, large nonGaussian potential with $f_{NL} \gtrsim 10^3$ leads to
earlier collapse of the first structures, more massive galaxies especially at
high redshifts, stronger clustering of galaxies, and higher halo bias. However,
for smaller NG with $f_{NL} \lesssim 10^2$, the effect is significantly weaker.
Observations of the distribution and properties of highredshift, rare objects
such as the first galaxies and quasars may provide further constraints on the
primordial NG.

Localtype primordial nonGaussianity couples statistics of the curvature
perturbation \zeta on vastly different physical scales. Because of this
coupling, statistics (i.e. the polyspectra) of \zeta in our Hubble volume may
not be representative of those in the larger universe  that is, they may be
biased. The bias depends on the local background value of \zeta, which includes
contributions from all modes with wavelength k ~< H_0 and is therefore enhanced
if the entire postinflationary patch is large compared with our Hubble volume.
We study the bias to locallymeasured statistics for general localtype
nonGaussianity. We consider three examples in detail: (i) the usual fNL, gNL
model, (ii) a strongly nonGaussian model with \zeta ~ \zeta_G^p, and (iii)
twofield nonGaussian initial conditions. In each scenario one may generate
statistics in a Hubblesize patch that are weakly Gaussian and consistent with
observations despite the fact that the statistics in the larger,
postinflationary patch look very different.

We report constraints on primordial nonGaussianity from the abundance of
Xray detected clusters. Our analytic prescription for adding nonGaussianity
to the cluster mass function takes into account moments beyond the skewness,
and we demonstrate that those moments should not be ignored in most analyses of
cluster data. We constrain the amplitude of the skewness for two scenarios that
have different overall levels of nonGaussianity, characterized by how
amplitudes of higher cumulants scale with the skewness. We find that current
data can constrain these oneparameter nonGaussian models at a useful level,
but are not sensitive to adding further details of the corresponding inflation
scenarios. Combining cluster data with Cosmic Microwave Background constraints
on the cosmology and power spectrum amplitude, we find the dimensionless
skewness to be 1000*M3=1+2428 for one of our scaling scenarios, and
1000*M3=4+/7 for the other. These are the first constraints on
nonGaussianity from Large Scale Structure that can be usefully applied to any
model of primordial nonGaussianity. When applied to the standard local ansatz,
the constraint corresponds to fNL=3+7891. When applied to a model with a
localshape bispectrum but higher cumulants that are relatively more important,
the amplitude of the localshape bispectrum is constrained to be fNL=14+2221.
This sensitivity implies that cluster counts could be used to distinguish
qualitatively different models for the primordial fluctuations that have
identical bispectra. (Abstract Abridged.)

We calculate Large Scale Structure observables for nonGaussianity arising
from nonBunchDavies initial states in single field inflation. These scenarios
can have substantial primordial nonGaussianity from squeezed (but observable)
momentum configurations. They generate a term in the halo bias that may be more
strongly scaledependent than the contribution from the local ansatz. We also
discuss theoretical considerations required to generate an observable
signature.

Primordial nonGaussianity is generated by interactions of the inflaton
field, either selfinteractions or couplings to other sectors. These two
physically different mechanisms can lead to nearly indistinguishable bispectra
of the equilateral type, but generate distinct patterns in the relative scaling
of higher order moments. We illustrate these classes in a simple effective
field theory framework where the flatness of the inflaton potential is
protected by a softly broken shift symmetry. Since the distinctive difference
between the two classes of interactions is the scaling of the moments, we
investigate the implications for observables that depend on the series of
moments. We obtain analytic expressions for the Minkowski functionals and the
halo mass function for an arbitrary structure of moments, and use these to
demonstrate how different classes of interactions might be distinguished
observationally. Our analysis casts light on a number of theoretical issues, in
particular we clarify the difference between the physics that keeps the
distribution of fluctuations nearly Gaussian, and the physics that keeps the
calculation under control.

In pursuing the intriguing resemblance of the Einstein equations to
thermodynamic equations, most sharply seen in systems possessing horizons, we
suggest that eternal inflation of the stochastic type may be a fruitful
phenomenon to explore. We develop a thermodynamic first law for quaside Sitter
space, valid on the horizon of a single observer's Hubble patch and explore
consistancy with previous proposals for horizons of various types in dynamic
and static situations. We use this framework to demonstrate that for the local
observer fluctuations of the type necessary for stochastic eternal inflation
fall within the regime where the thermodynamic approach is believed to apply.
This scenario is interesting because of suggestive parallels to black hole
evaporation.