
The kinetic Sunyaev Zel'dovich (kSZ) and polarized Sunyaev Zel'dovich (pSZ)
effects are temperature and polarization anisotropies induced by the scattering
of CMB photons from structure in the postreionization Universe. In the case of
the kSZ effect, small angular scale anisotropies in the optical depth are
modulated by the cosmic microwave background (CMB) dipole field, i.e. the CMB
dipole observed at each spacetime point, which is sourced by the primordial
dipole and especially the local peculiar velocity. In the case of the pSZ
effect, similar smallscale anisotropies are modulated by the CMB quadrupole
field, which receives contributions from both scalar and tensor modes.
Statistical anisotropies in the cross correlations of CMB temperature and
polarization with tracers of the inhomogeneous distribution of electrons
provide a means of isolating and reconstructing the dipole and quadrupole
fields. In this paper, we present a set of unbiased minimum variance quadratic
estimators for the reconstruction of the dipole and quadrupole fields, and
forecast the ability of future CMB experiments and large scale structure
surveys to perform this reconstruction. Consistent with previous work, we find
that a high fidelity reconstruction of the dipole and quadrupole fields over a
variety of scales is indeed possible, and demonstrate the sensitivity of the
pSZ effect to primordial tensor modes. Using a principle component analysis, we
estimate how many independent modes could be accessed in such a reconstruction.
We also comment on a few first applications of a detection of the dipole and
quadrupole fields, including a reconstruction of the primordial contribution to
our locally observed CMB dipole, a test of statistical homogeneity on large
scales from the first modes of the quadrupole field, and a reconstruction
technique for the primordial potential on the largest scales.

Secondary CMB polarization is induced by the latetime scattering of CMB
photons by free electrons on our past light cone. This polarized Sunyaev
Zel'dovich (pSZ) effect is sensitive to the electrons' locally observed CMB
quadrupole, which is sourced primarily by long wavelength inhomogeneities. By
combining the remote quadrupoles measured by free electrons throughout the
Universe after reionization, the pSZ effect allows us to obtain additional
information about large scale modes beyond what can be learned from our own
last scattering surface. Here we determine the power of pSZ tomography, in
which the pSZ effect is crosscorrelated with the density field binned at
several redshifts, to provide information about the long wavelength Universe.
The signal we explore here is a power asymmetry in the crosscorrelation
between $E$ or $B$ mode CMB polarization and the density field. We compare this
to the cosmic variance limited noise: the random chance to get a power
asymmetry in the absence of a large scale quadrupole field. By computing the
necessary transfer functions and crosscorrelations, we compute the
signaltonoise ratio attainable by idealized next generation CMB experiments
and galaxy surveys. We find that a signaltonoise ratio of $\sim 110$ is in
principle attainable over a significant range of power multipoles, with the
strongest signal coming from the first multipoles in the lowest redshift bins.
These results prompt further assessment of realistically measuring the pSZ
signal and the potential impact for constraining cosmology on large scales.

It is generally expected that heavy fields are present during inflation,
which can leave their imprint in latetime cosmological observables. The main
signature of these fields is a small amount of distinctly shaped
nonGaussianity, which if detected, would provide a wealth of information about
the particle spectrum of the inflationary Universe. Here we investigate to what
extent these signatures can be detected or constrained using futuristic 21cm
surveys. We construct modelindependent templates that extract the
squeezedlimit behavior of the bispectrum, and examine their overlap with
standard inflationary shapes and secondary nonGaussianities. We then use these
templates to forecast detection thresholds for different masses and couplings
using a 3D reconstruction of modes during the dark ages ($z\sim 30100$). We
consider interactions of several broad classes of models and quantify their
detectability as a function of the baseline of a dark ages interferometer. Our
analysis shows that there exists the tantalizing possibility of discovering new
particles with different masses and interactions with future 21cm surveys.

Detecting a deviation from a featureless primordial power spectrum of
fluctuations would give profound insight into the physics of the primordial
Universe. Depending on their nature, primordial features can either provide
direct evidence for the inflation scenario or pin down details of the inflation
model. Thus far, using the cosmic microwave background (CMB) we have only been
able to put stringent constraints on the amplitude of features, but no
significant evidence has been found for such signals. Here we explore the limit
of the experimental reach in constraining such features using 21 cm tomography
at high redshift. A measurement of the 21 cm power spectrum from the Dark Ages
is generally considered as the ideal experiment for early Universe physics,
with potentially access to a large number of modes. We consider three different
categories of theoretically motivated models: the sharp feature models,
resonance models, and standard clock models. We study the improvements on
bounds on features as a function of the total number of observed modes and
identify parameter degeneracies. The detectability depends critically on the
amplitude, frequency and scalelocation of the features, as well as the angular
and redshift resolution of the experiment. We quantify these effects by
considering different fiducial models. Our forecast shows that a cosmic
variance limited 21 cm experiment measuring fluctuations in the redshift range
$30\leq z \leq 100$ with a 0.01MHz bandwidth and subarcminute angular
resolution could potentially improve bounds by several orders of magnitude for
most features compared to current Planck bounds. At the same time, 21 cm
tomography also opens up a unique window into features that are located on very
small scales.

We develop the tools necessary to assess the statistical significance of
resonant features in the CMB correlation functions, combining power spectrum
and bispectrum measurements. This significance is typically addressed by
running a large number of simulations to derive the probability density
function (PDF) of the featureamplitude in the Gaussian case. Although these
simulations are tractable for the power spectrum, for the bispectrum they
require significant computational resources. We show that, by assuming that the
PDF is given by a multivariate Gaussian where the covariance is determined by
the Fisher matrix of the sine and cosine terms, we can efficiently produce
spectra that are statistically close to those derived from full simulations. By
drawing a large number of spectra from this PDF, both for the power spectrum
and the bispectrum, we can quickly determine the statistical significance of
candidate signatures in the CMB, considering both single frequency and
multifrequency estimators. We show that for resonance models, cosmology and
foreground parameters have little influence on the estimated amplitude, which
allows to simplify the analysis considerably. A more precise likelihood
treatment can then be applied to candidate signatures only. We also discuss a
modal expansion approach for the power spectrum, aimed at quickly scanning
through large families of oscillating models.

We propose optimal estimators for bispectra from excited states. Two common
properties of such bispectra are the enhancement in the collinear limit, and
the prediction of oscillating features. We review the physics behind excited
states and some of the choices made in the literature. We show that the
enfolded template is a good template in the collinear limit, but does poorly
elsewhere, establishing a strong case for an improved estimator. Although the
detailed scale dependence of the bispectra differs depending on various
assumptions, generally the predicted bispectra are either effectively 1 or
2dimensional and a simple Fourier basis suffices for accurate reconstruction.
For an optimal CMB data analysis, combining all npoint functions, the choice
for the excited state needs to be the same when computing power spectrum,
bispectrum and higher order correlation functions. This has not always been the
case, which could lead to wrong conclusions. We calculate the bispectrum for
different choices previously discussed for the power spectrum, setting up a
consistent framework to search for evidence of excited states in the CMB data.

We propose an (optimal) estimator for a CMB bispectrum containing
logarithmically spaced oscillations. There is tremendous theoretical interest
in such bispectra, and they are predicted by a plethora of models, including
axion monodromy models of inflation and initial state modifications. The number
of resolved logarithmical oscillations in the bispectrum is limited due to the
discrete resolution of the multipole bispectrum. We derive a simple relation
between the maximum number of resolved oscillations and the frequency. We
investigate several ways to factorize the primordial bispectrum, and conclude
that a one dimensional expansion in the sum of the momenta $\sum k_i = k_t$ is
the most efficient and flexible approach. We compare the expansion to the exact
result in multipole space and show for $\omega_{\rm eff}=100$ that
$\mathcal{O}(10^3)$ modes are sufficient for an accurate reconstruction. We
compute the expected $\sigma_{f_{\rm NL}}$ and find that within an effective
field theory (EFT) the overall signal to noise scales as $S/N\propto
\omega^{3/2}$. Using only the temperature data we find
$S/N\sim\mathcal{O}(110^2)$ for the frequency domain set by the EFT.

Contributions of the Pierre Auger Collaboration to the 33rd International
Cosmic Ray Conference, Rio de Janeiro, Brazil, July 2013