
Strong lensing is a sensitive probe of the smallscale density fluctuations
in the Universe. We implement a novel approach to modeling strongly lensed
systems using probabilistic cataloging, which is a transdimensional,
hierarchical, and Bayesian framework to sample from a metamodel (union of
models with different dimensionality) consistent with observed photon count
maps. Probabilistic cataloging allows us to robustly characterize modeling
covariances within and across lens models with different numbers of subhalos.
Unlike traditional cataloging of subhalos, it does not require model subhalos
to improve the goodness of fit above the detection threshold. Instead, it
allows the exploitation of all information contained in the photon count maps,
for instance, when constraining the subhalo mass function. We further show
that, by not including these small subhalos in the lens model,
fixeddimensional inference methods can significantly mismodel the data. Using
a simulated Hubble Space Telescope (HST) dataset, we show that the subhalo mass
function can be probed even when many subhalos in the sample catalogs are
individually below the detection threshold and would be absent in a traditional
catalog. With the planned Wide Field Infrared Space Telescope (WFIRST),
simultaneous probabilistic cataloging of dark subhalos in highresolution, deep
strong lens images has the potential to constrain the subhalo mass function at
even lower masses.

We present updated constraints on the freestreaming nature of cosmological
neutrinos from cosmic microwave background (CMB) power spectra, baryonic
acoustic oscillation data, and local measurements of the Hubble constant.
Specifically, we consider a Fermilike fourfermion interaction between
massless neutrinos, characterized by an effective coupling constant $ G_{\rm
eff}$, and resulting in a neutrino opacity $\dot{\tau}_\nu\propto G_{\rm eff}^2
T_\nu^5$. Using a conservative prior on the parameter $\log_{10}\left(G_{\rm
eff} {\rm MeV}^2\right)$, we find a bimodal posterior distribution. The first
of these modes is consistent with the standard $\Lambda$CDM cosmology and
corresponds to neutrinos decoupling at redshift $z_{\nu,{\rm dec}} >
1.3\times10^5$. The other mode of the posterior, dubbed the "interacting
neutrino mode", corresponds to neutrino decoupling occurring within a narrow
redshift window centered around $z_{\nu,{\rm dec}}\sim8300$. This mode is
characterized by a high value of the effective neutrino coupling constant,
together with a lower value of the scalar spectral index and amplitude of
fluctuations, and a higher value of the Hubble parameter. Using both a maximum
likelihood analysis and the ratio of the two mode's Bayesian evidence, we find
the interacting neutrino mode to be statistically disfavored compared to the
standard $\Lambda$CDM cosmology. Interestingly, the addition of CMB
polarization and direct Hubble constant measurements significantly raises the
statistical significance of this secondary mode, indicating that new physics in
the neutrino sector could help explain the difference between local
measurements of $H_0$, and those inferred from CMB data. A robust consequence
of our results is that neutrinos must be free streaming long before the epoch
of matterradiation equality.

Studying the smallest selfbound dark matter structure in our Universe can
yield important clues about the fundamental particle nature of dark matter.
Galaxyscale strong gravitational lensing provides a unique way to detect and
characterize dark matter substructures at cosmological distances from the Milky
Way. Within the cold dark matter (CDM) paradigm, the number of lowmass
subhalos within lens galaxies is expected to be large, implying that their
contribution to the lensing convergence field is approximately Gaussian and
could thus be described by their power spectrum. We develop here a general
formalism to compute from first principles the substructure convergence power
spectrum for different populations of dark matter subhalos. As an example, we
apply our framework to two distinct subhalo populations: a truncated
NavarroFrenkWhite subhalo population motivated by standard CDM, and a
truncated cored subhalo population motivated by selfinteracting dark matter
(SIDM). We study in detail how the subhalo abundance, mass function, internal
density profile, and concentration affect the amplitude and shape of the
substructure power spectrum. We determine that the power spectrum is mostly
sensitive to a specific combination of the subhalo abundance and moments of the
mass function, as well as to the average tidal truncation scale of the largest
subhalos included in the analysis. Interestingly, we show that the asymptotic
slope of the substructure power spectrum at large wave number reflects the
internal density profile of the subhalos. In particular, the SIDM power
spectrum exhibits a characteristic steepening at large wave number absent in
the CDM power spectrum, opening the possibility of using this observable, if at
all measurable, to discern between these two scenarios.

We formulate an effective theory of structure formation (ETHOS) that enables
cosmological structure formation to be computed in almost any microphysical
model of dark matter physics. This framework maps the detailed microphysical
theories of particle dark matter interactions into the physical effective
parameters that shape the linear matter power spectrum and the selfinteraction
transfer cross section of nonrelativistic dark matter. These are the input to
structure formation simulations, which follow the evolution of the cosmological
and galactic dark matter distributions. Models with similar effective
parameters in ETHOS but with different dark particle physics would nevertheless
result in similar dark matter distributions. We present a general method to map
an ultraviolet complete or effective field theory of low energy dark matter
physics into parameters that affect the linear matter power spectrum and carry
out this mapping for several representative particle models. We further propose
a simple but useful choice for characterizing the dark matter selfinteraction
transfer cross section that parametrizes selfscattering in structure formation
simulations. Taken together, these effective parameters in ETHOS allow the
classification of dark matter theories according to their structure formation
properties rather than their intrinsic particle properties, paving the way for
future simulations to span the space of viable dark matter physics relevant for
structure formation.

The amplitude of largescale matter fluctuations inferred from the observed
SunyaevZeldovich (SZ) cluster mass function and from weak gravitational
lensing studies, when taken at face value, is in tension with measurements of
the cosmic microwave background (CMB) and baryon acoustic oscillation (BAO). In
this work, we revisit whether this possible discrepancy can be attributed to
new interactions in the dark matter sector. Focusing on a cosmological model
where dark matter interacts with a dark radiation species until the epoch of
matterradiation equality, we find that measurements of the Lymanalpha flux
power spectrum from the Sloan Digital Sky Survey provides no support to the
hypothesis that new dark matter interactions can resolve the possible tension
between CMB and largescale structure (LSS). Indeed, while the addition of dark
matterdark radiation interactions leads to an improvement of
$2\Delta\ln\mathcal{L}=12$ with respect to the standard $\Lambda$ cold dark
matter ($\Lambda$CDM) model when only CMB, BAO, and LSS data are considered,
the inclusion of Lymanalpha data reduces the improvement of the fit to
$2\Delta\ln\mathcal{L}=6$ relative to $\Lambda$CDM. We thus conclude that the
statistical evidence for new dark matter interactions (largely driven by the
Planck SZ dataset) is marginal at best, and likely caused by systematics in the
data. We also perform a Fisher forecast analysis for the reach of a future
dataset composed of a CMBS4 experiment combined with the Large Synoptic Survey
Telescope galaxy survey. We find that the constraint on the effective number of
fluidlike dark radiation species, $\Delta N_{\rm fluid}$, will be improved by
an order of magnitude compared to current bounds.

Recently it was shown that dark matter with mass of order the weak scale can
be charged under a new longrange force, decoupled from the Standard Model,
with only weak constraints from early Universe cosmology. Here we consider the
implications of an additional charged particle $C$ that is light enough to lead
to significant dissipative dynamics on galactic times scales. We highlight
several novel features of this model, which can be relevant even when the $C$
particle constitutes only a small fraction of the number density (and energy
density). We assume a small asymmetric abundance of the $C$ particle whose
charge is compensated by a heavy $X$ particle so that the relic abundance of
dark matter consists mostly of symmetric $X$ and $\bar{X}$, with a small
asymmetric component made up of $X$ and $C$. As the universe cools, it
undergoes asymmetric recombination binding the free $C$s into $(XC)$ dark atoms
efficiently. Even with a tiny asymmetric component, the presence of $C$
particles catalyzes tight coupling between the heavy dark matter $X$ and the
dark photon plasma that can lead to a significant suppression of the matter
power spectrum on small scales and lead to some of the strongest bounds on such
dark matter theories. We find a viable parameter space where structure
formation constraints are satisfied and significant dissipative dynamics can
occur in galactic haloes but show a large region is excluded. Our model shows
that subdominant components in the dark sector can dramatically affect
structure formation.

We revisit constraints on dark matter that is charged under a $U(1)$ gauge
group in the dark sector, decoupled from Standard Model forces. We find that
the strongest constraints in the literature are subject to a number of
mitigating factors. For instance, the naive dark matter thermalization
timescale in halos is corrected by saturation effects that slow down
isotropization for modest ellipticities. The weakened bounds uncover
interesting parameter space, making models with weakscale charged dark matter
viable, even with electromagnetic strength interaction. This also leads to the
intriguing possibility that dark matter selfinteractions within small dwarf
galaxies are extremely large, a relatively unexplored regime in current
simulations. Such strong interactions suppress heat transfer over scales larger
than the dark matter mean free path, inducing a dynamical cutoff length scale
above which the system appears to have only feeble interactions. These effects
must be taken into account to assess the viability of darklycharged dark
matter. Future analyses and measurements should probe a promising region of
parameter space for this model.

In the context of the ESA M5 (medium mission) call we proposed a new
satellite mission, Theia, based on relative astrometry and extreme precision to
study the motion of very faint objects in the Universe. Theia is primarily
designed to study the local dark matter properties, the existence of Earthlike
exoplanets in our nearest star systems and the physics of compact objects.
Furthermore, about 15 $\%$ of the mission time was dedicated to an open
observatory for the wider community to propose complementary science cases.
With its unique metrology system and "point and stare" strategy, Theia's
precision would have reached the sub microarcsecond level. This is about 1000
times better than ESA/Gaia's accuracy for the brightest objects and represents
a factor 1030 improvement for the faintest stars (depending on the exact
observational program). In the version submitted to ESA, we proposed an optical
(3501000nm) onaxis TMA telescope. Due to ESA Technology readiness level, the
camera's focal plane would have been made of CCD detectors but we anticipated
an upgrade with CMOS detectors. Photometric measurements would have been
performed during slew time and stabilisation phases needed for reaching the
required astrometric precision.

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.

In the standard structure formation scenario based on the cold dark matter
paradigm, galactic halos are predicted to contain a large population of dark
matter subhalos. While the most massive members of the subhalo population can
appear as luminous satellites and be detected in optical surveys, establishing
the existence of the low mass and mostly dark subhalos has proven to be a
daunting task. Galaxyscale strong gravitational lenses have been successfully
used to study mass substructures lying close to lensed images of bright
background sources. However, in typical galaxyscale lenses, the strong lensing
region only covers a small projected area of the lens's dark matter halo,
implying that the vast majority of subhalos cannot be directly detected in
lensing observations. In this paper, we point out that this large population of
dark satellites can collectively affect gravitational lensing observables,
hence possibly allowing their statistical detection. Focusing on the region of
the galactic halo outside the strong lensing area, we compute from first
principles the statistical properties of perturbations to the gravitational
time delay and position of lensed images in the presence of a mass substructure
population. We find that in the standard cosmological scenario, the statistics
of these lensing observables are well approximated by Gaussian distributions.
The formalism developed as part of this calculation is very general and can be
applied to any halo geometry and choice of subhalo mass function. Our results
significantly reduce the computational cost of including a large substructure
population in lens models and enable the use of Bayesian inference techniques
to detect and characterize the distributed satellite population of distant lens
galaxies.

We perform a detailed study of the weak interactions of standard model
neutrinos with the primordial plasma and their effect on the resonant
production of sterile neutrino dark matter. Motivated by issues in cosmological
structure formation on small scales, and reported Xray signals that could be
due to sterile neutrino decay, we consider $7$ keVscale sterile neutrinos.
Oscillationdriven production of such sterile neutrinos occurs at temperatures
$T \gtrsim 100$ MeV, where we study two significant effects of weakly charged
species in the primordial plasma: (1) the redistribution of an input lepton
asymmetry; (2) the opacity for active neutrinos. We calculate the
redistribution analytically above and below the quarkhadron transition, and
match with lattice QCD calculations through the transition. We estimate
opacities due to tree level processes involving leptons and quarks above the
quarkhadron transition, and the most important mesons below the transition. We
report final sterile neutrino dark matter phase space densities that are
significantly influenced by these effects, and yet relatively robust to
remaining uncertainties in the nature of the quarkhadron transition. We also
provide transfer functions for cosmological density fluctuations with cutoffs
at $k \simeq 10 \ h \ {\rm Mpc}^{1}$, that are relevant to galactic structure
formation.

In the standard model neutrinos are assumed to have streamed across the
Universe since they last scattered at the weak decoupling epoch when the
temperature of the standardmodel plasma was ~MeV. The shear stress of
freestreaming neutrinos imprints itself gravitationally on the Cosmic
Microwave Background (CMB) and makes the CMB a sensitive probe of neutrino
scattering. Yet, the presence of nonstandard physics in the neutrino sector may
alter this standard chronology and delay neutrino freestreaming until a much
later epoch. We use observations of the CMB to constrain the strength of
neutrino selfinteractions G_eff and put limits on new physics in the neutrino
sector from the early Universe. Recent measurements of the CMB at large
multipoles made by the Planck satellite and highl experiments are critical for
probing this physics. Within the context of conventional LambdaCDM parameters
cosmological data are compatible with G_eff < 1/(56 MeV)^2 and neutrino
freestreaming might be delayed until their temperature has cooled to as low as
~25 eV. Intriguingly, we also find an alternative cosmology compatible with
cosmological data in which neutrinos scatter off each other until z~10^4 with a
preferred interaction strength in a narrow region around $G_{\rm eff} \simeq
1/({\rm 10 \, MeV})^{2} \simeq 8.6\times10^8 G_{\rm F}$, where $G_{\rm F}$ is
the Fermi constant. This distinct selfinteracting neutrino cosmology is
characterized by somewhat lower values of both the scalar spectral index and
the amplitude of primordial fluctuations. While we phrase our discussion here
in terms of a specific scenario in which a late onset of neutrino
freestreaming could occur, our constraints on the neutrino visibility function
are very general.

We demonstrate that selfinteracting dark matter models with interactions
mediated by light particles can have significant deviations in the matter
powerspectrum and detailed structure of galactic halos when compared to a
standard cold dark matter scenario. While these deviations can take the form of
suppression of small scale structure that are in some ways similar to that of
warm dark matter, the selfinteracting models have a much wider range of
possible phenomenology. A longrange force in the dark matter can introduce
multiple scales to the initial power spectrum, in the form of dark acoustic
oscillations and an exponential cutoff in the power spectrum. Using
simulations we show that the impact of these scales can remain observationally
relevant up to the present day. Furthermore, the selfinteraction can continue
to modify the smallscale structure of the dark matter halos, reducing their
central densities and creating a dark matter core. The resulting phenomenology
is unique to this type of models.

If all or a fraction of the dark matter (DM) were coupled to a bath of dark
radiation (DR) in the early Universe we expect the combined DMDR system to
give rise to acoustic oscillations of the dark matter until it decouples from
the DR. Much like the standard baryon acoustic oscillations, these dark
acoustic oscillations (DAO) imprint a characteristic scale, the sound horizon
of dark matter, on the matter power spectrum. We compute in detail how the
microphysics of the DMDR interaction affects the clustering of matter in the
Universe and show that the DAO physics also gives rise to unique signatures in
the temperature and polarization spectra of the cosmic microwave background
(CMB). We use cosmological data from the CMB, baryon acoustic oscillations
(BAO), and largescale structure to constrain the possible fraction of
interacting DM as well as the strength of its interaction with DR. Like nearly
all knowledge we have gleaned about dark matter since inferring its existence
this constraint rests on the betrayal by gravity of the location of otherwise
invisible dark matter. Although our results can be straightforwardly applied to
a broad class of models that couple dark matter particles to various light
relativistic species, in order to make quantitative predictions, we model the
interacting component as dark atoms coupled to a bath of dark photons. We find
that linear cosmological data and CMB lensing put strong constraints on
existence of DAO features in the CMB and the largescale structure of the
Universe. Interestingly, we find that at most ~5% of all DM can be very
strongly interacting with DR. We show that our results are surprisingly
constraining for the recently proposed Doubledisk DM model, a novel example of
how largescale precision cosmological data can be used to constrain galactic
physics and subgalactic structure.

While, to ensure successful cosmology, dark matter (DM) must kinematically
decouple from the standard model plasma very early in the history of the
Universe, it can remain coupled to a bath of "dark radiation" until a
relatively late epoch. One minimal theory that realizes such a scenario is the
Atomic Dark Matter model, in which two fermions oppositely charged under a new
U(1) dark force are initially coupled to a thermal bath of "dark photons" but
eventually recombine into neutral atomlike bound states and begin forming
gravitationallybound structures. As dark atoms have (dark) atomsized
geometric cross sections, this model also provides an example of
selfinteracting DM with a velocitydependent cross section. Delayed kinetic
decoupling in this scenario predicts novel DM properties on small scales but
retains the success of cold DM on larger scales. We calculate the atomic
physics necessary to capture the thermal history of this dark sector and show
significant improvements over the standard atomic hydrogen calculation are
needed. We solve the Boltzmann equations that govern the evolution of
cosmological fluctuations in this model and find in detail the impact of the
atomic DM scenario on the matter power spectrum and the cosmic microwave
background (CMB). This scenario imprints a new length scale, the
DarkAcousticOscillation (DAO) scale, on the matter density field. This DAO
scale shapes the smallscale matter power spectrum and determines the minimal
DM halo mass at late times which may be many orders of magnitude larger than in
a typical WIMP scenario. This model necessarily includes an extra dark
radiation component, which may be favoured by current CMB experiments, and we
quantify CMB signatures that distinguish an atomic DM scenario from a standard
$\Lambda$CDM model containing extra freestreaming particles. [Abridged]

Many models of inflation predict oscillatory features in the bispectrum of
primordial fluctuations. Since it has been shown that primordial
nonGaussianity can lead to a scaledependent halo bias, we investigate the
effect of oscillations in the threepoint function on the clustering of
darkmatter halos. Interestingly, we find that features in the inflaton
potential such as oscillations or sharp steps get imprinted in the mass
dependence of the nonGaussian halo bias. In this paper, we focus on models
displaying a sharp feature in the inflaton potential as well as Resonant
nonGaussianity. In both cases, we find a strong scale dependence for the
nonGaussian halo bias with a slope similar to that of the local model. In the
resonant case, we find that the nonGaussian bias oscillates with halo mass, a
novel feature that is unique to this type of models. In the case of a sharp
feature in the inflaton potential, we find that the clustering of halos is
enhanced at the mass scale corresponding to the Fourier mode that exited the
horizon when the inflaton was crossing the feature in the potential. Both of
these are new effects that open the possibility of characterizing the
inflationary potential with largescalestructure surveys. We briefly discuss
the prospects for detecting these nonGaussian effects.

Prior to recombination photons, electrons, and atomic nuclei rapidly
scattered and behaved, almost, like a single tightlycoupled photonbaryon
plasma. We investigate here the accuracy of the tightcoupling approximation
commonly used to numerically evolve the baryon and photon perturbation
equations at early times. By solving the exact perturbations equations with a
stiff solver starting deep in the radiationdominated epoch we find the level
of inaccuracy introduced by resorting to the standard firstorder
tightcoupling approximation. We develop a new secondorder approximation in
the inverse Thomson opacity expansion and show that it closely tracks the full
solution, at essentially no extra numerical cost. We find the bias on estimates
of cosmological parameters introduced by the firstorder approximation is, for
most parameters, negligible. Finally, we show that our secondorder
approximation can be used to reduce the time needed to compute cosmic microwave
background angular spectra by as much as ~17%.

Reheating is an important part of inflationary cosmology. It describes the
production of Standard Matter particles after the phase of accelerated
expansion. We give a review of the reheating process, focusing on an indepth
discussion of the preheating stage which is characterized by exponential
particle production due to a parametric resonance or tachyonic instability. We
give a brief overview of the thermalization process after preheating and end
with a survey of some applications to supersymmetric theories and to other
issues in cosmology such as baryogenesis, dark matter and metric preheating.

It has recently been argued that the PAMELA, ATIC and PPBBETS data showing
an anomalous excess of highenergy cosmic ray positrons and electrons might be
explained by dark matter annihilating in the Galactic halo with a cross section
resonantly enhanced compared to its value in the primeval plasma. We find that
with a very large annihilation cross section the flash of energetic photons and
electronpositron pairs expected from darkmatter annihilation in the first
protohalos that form at redshift z~40 is likely substantial and observable. As
a consequence, bounds on the allowed energy injection into the primordial gas
and the energy density of the diffuse gammaray background give rise to limits
on the lowvelocity dark matter cross section that can be difficult to
reconcile with this interpretation of the PAMELA, ATIC and PPBBETS results.

We study how the resonant decay of moduli fields arising in the Minimal
Supersymmetric Standard Model (MSSM) could affect large scale curvature
perturbations in the early universe. It has been known for some time that the
presence of entropy perturbations in a multicomponent system can act as seeds
for the curvature perturbations on all scales. These entropy perturbations
could be amplified exponentially if one of the moduli decays via stochastic
resonance, affecting the curvature power spectrum in the process. By imposing
the COBE normalization on this power spectrum, one could put constraints on the
masses and couplings of the underlying particle physics model without having to
rely on collider experiments. We discuss in detail the case of the MSSM but
this method could be applied to other theories beyond the Standard Model.