
Galaxy spectra are essential to probe the spatial distribution of galaxies in
our Universe. To better interpret current and future spectroscopic galaxy
redshift surveys, it is important to be able to simulate these data sets. We
describe Uspec, a forward modeling tool to generate galaxy spectra taking into
account intrinsic galaxy properties as well as instrumental responses of a
given telescope. The model for the intrinsic properties of the galaxy
population was developed in an earlier work for broadband imaging surveys [1].
We apply Uspec to the SDSS/CMASS sample of Luminous Red Galaxies (LRGs). We
construct selection cuts that match those used to build this LRG sample, which
we then apply to data and simulations in the same way. The resulting real and
simulated average spectra show a very good agreement overall, with the
simulated one showing a slightly bluer galaxy population. For a quantitative
comparison, we perform Principal Component Analysis (PCA) of the sets of
spectra. By comparing the PCs constructed from simulations and data, we find
very good agreement for the first four components, and moderate for the fifth.
The distributions of the eigencoefficients also show an appreciable overlap. We
are therefore able to properly simulate the LRG sample taking into account the
SDSS/BOSS instrumental responses. The small residual differences between the
two samples can be ascribed to the intrinsic properties of the simulated galaxy
population, which can be reduced by adjusting the model parameters in the
future. This provides good prospects for the forward modeling of upcoming large
spectroscopic surveys.

Upcoming weak lensing surveys will probe large fractions of the sky with
unprecedented accuracy. To infer cosmological constraints, a large ensemble of
survey simulations are required to accurately model cosmological observables
and their covariances. We develop a parallelized multilensplane pipeline
called UFalcon, designed to generate fullsky weak lensing maps from lightcones
within a minimal runtime. It makes use of LPICOLA, an approximate numerical
code, which provides a fast and accurate alternative to cosmological $N$Body
simulations. The UFalcon maps are constructed by nesting 2 simulations covering
a redshiftrange from $z=0.1$ to $1.5$ without replicating the simulation
volume. We compute the convergence and projected overdensity maps for LPICOLA
in the lightcone or snapshot mode. The generation of such a map, including the
LPICOLA simulation, takes about 3 hours walltime on 220 cores. We use the maps
to calculate the spherical harmonic power spectra, which we compare to
theoretical predictions and to UFalcon results generated using the full
$N$Body code GADGET2. We then compute the covariance matrix of the fullsky
spherical harmonic power spectra using 150 UFalcon maps based on LPICOLA in
lightcone mode. We consider the PDF, the higherorder moments and the variance
of the smoothed field variance to quantify the accuracy of the covariance
matrix, which we find to be a few percent for scales $\ell \sim 10^2$ to
$10^3$. We test the impact of this level of accuracy on cosmological
constraints using an optimistic survey configuration, and find that the final
results are robust to this level of uncertainty. The speed and accuracy of our
developed pipeline provides a basis to also include further important features
such as masking, varying noise and will allow us to compute covariance matrices
for models beyond $\Lambda$CDM. [abridged]

Assessing the consistency of parameter constraints derived from different
cosmological probes is an important way to test the validity of the underlying
cosmological model. In an earlier work [Nicola et al., 2017], we computed
constraints on cosmological parameters for $\Lambda$CDM from an integrated
analysis of CMB temperature anisotropies and CMB lensing from Planck, galaxy
clustering and weak lensing from SDSS, weak lensing from DES SV as well as Type
Ia supernovae and Hubble parameter measurements. In this work, we extend this
analysis and quantify the concordance between the derived constraints and those
derived by the Planck Collaboration as well as WMAP9, SPT and ACT. As a measure
for consistency, we use the Surprise statistic [Seehars et al., 2014], which is
based on the relative entropy. In the framework of a flat $\Lambda$CDM
cosmological model, we find all data sets to be consistent with one another at
a level of less than 1$\sigma$. We highlight that the relative entropy is
sensitive to inconsistencies in the models that are used in different parts of
the analysis. In particular, inconsistent assumptions for the neutrino mass
break its invariance on the parameter choice. When consistent model assumptions
are used, the data sets considered in this work all agree with each other and
$\Lambda$CDM, without evidence for tensions.

Determining the redshift distribution $n(z)$ of galaxy samples is essential
for several cosmological probes including weak lensing. For imaging surveys,
this is usually done using photometric redshifts estimated on an
objectbyobject basis. We present a new approach for directly measuring the
global $n(z)$ of cosmological galaxy samples, including uncertainties, using
forward modeling. Our method relies on image simulations produced using UFig
(Ultra Fast Image Generator) and on ABC (Approximate Bayesian Computation)
within the $MCCL$ (MonteCarlo Control Loops) framework. The galaxy population
is modeled using parametric forms for the luminosity functions, spectral energy
distributions, sizes and radial profiles of both blue and red galaxies. We
apply exactly the same analysis to the real data and to the simulated images,
which also include instrumental and observational effects. By adjusting the
parameters of the simulations, we derive a set of acceptable models that are
statistically consistent with the data. We then apply the same cuts to the
simulations that were used to construct the target galaxy sample in the real
data. The redshifts of the galaxies in the resulting simulated samples yield a
set of $n(z)$ distributions for the acceptable models. We demonstrate the
method by determining $n(z)$ for a cosmic shear like galaxy sample from the
4band Subaru SuprimeCam data in the COSMOS field. We also complement this
imaging data with a spectroscopic calibration sample from the VVDS survey. We
compare our resulting posterior $n(z)$ distributions to the one derived from
photometric redshifts estimated using 36 photometric bands in COSMOS and find
good agreement. This offers good prospects for applying our approach to current
and future large imaging surveys.

Weak Gravitational Lensing is a powerful probe of the dark sector of the
Universe. One of the main challenges for this technique is the treatment of
systematics in the measurement of cosmic shear from galaxy shapes. In an
earlier work, Refregier & Amara (2014) have proposed the Monte Carlo Control
Loops (MCCL) to overcome these effects using a forward modeling approach. We
focus here on one of the control loops in this method, the task of which is the
calibration of the shear measurement. For this purpose, we first consider the
requirements on the shear systematics for a given survey and propagate them to
different systematics terms. We use two onepoint statistics to calibrate the
shear measurement and six further onepoint statistics as diagnostics. We also
propagate the systematics levels that we estimate from the onepoint functions
to the twopoint functions for the different systematic error sources. This
allows us to assess the consistency between the systematics levels measured in
different ways. To test the method, we construct synthetic sky surveys with an
area of 1,700 deg$^2$. With some simplifying assumptions, we are able to meet
the requirements on the shear calibration for this survey configuration.
Furthermore, we account for the total residual shear systematics in terms of
the contributing sources. We discuss how this MCCL framework can be applied to
current and future weak lensing surveys.

Recent progress in cosmology has relied on combining different cosmological
probes. In earlier work, we implemented an integrated approach to cosmology
where the probes are combined into a common framework at the map level. This
has the advantage of taking full account of the correlations between the
different probes, to provide a stringent test of systematics and of the
validity of the cosmological model. We extend this analysis to include not only
CMB temperature, galaxy clustering, weak lensing from SDSS but also CMB
lensing, weak lensing from the DES SV survey, Type Ia SNe and $H_{0}$
measurements. This yields 12 auto and crosspower spectra as well as
background probes. Furthermore, we extend the treatment of systematic
uncertainties. For $\Lambda$CDM, we find results that are consistent with our
earlier work. Given our enlarged data set and systematics treatment, this
confirms the robustness of our analysis and results. Furthermore, we find that
our bestfit cosmological model gives a good fit to the data we consider with
no signs of tensions within our analysis. We also find our constraints to be
consistent with those found by WMAP9, SPT and ACT and the KiDS weak lensing
survey. Comparing with the Planck Collaboration results, we see a broad
agreement, but there are indications of a tension from the marginalized
constraints in most pairs of cosmological parameters. Since our analysis
includes CMB temperature Planck data at $10 < \ell < 610$, the tension appears
to arise between the Planck high$\ell$ and the other measurements.
Furthermore, we find the constraints on the probe calibration parameters to be
in agreement with expectations, showing that the data sets are mutually
consistent. In particular, this yields a confirmation of the amplitude
calibration of the weak lensing measurements from SDSS, DES SV and Planck CMB
lensing from our integrated analysis. [abridged]

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.

Recent observational progress has led to the establishment of the standard
$\Lambda$CDM model for cosmology. This development is based on different
cosmological probes that are usually combined through their likelihoods at the
latest stage in the analysis. We implement here an integrated scheme for
cosmological probes, which are combined in a common framework starting at the
map level. This treatment is necessary as the probes are generally derived from
overlapping maps and are thus not independent. It also allows for a thorough
test of the cosmological model and of systematics through the consistency of
different physical tracers. As a first application, we combine current
measurements of the Cosmic Microwave Background (CMB) from the Planck
satellite, and galaxy clustering and weak lensing from SDSS. We consider the
spherical harmonic power spectra of these probes including all six auto and
crosscorrelations along with the associated full Gaussian covariance matrix.
This provides an integrated treatment of different analyses usually performed
separately including CMB anisotropies, cosmic shear, galaxy clustering,
galaxygalaxy lensing and the Integrated SachsWolfe (ISW) effect with galaxy
and shear tracers. We derive constraints on $\Lambda$CDM parameters that are
compatible with existing constraints and highlight tensions between data sets,
which become apparent in this integrated treatment. We discuss how this
approach provides a complete and powerful integrated framework for probe
combination and how it can be extended to include other tracers in the context
of current and future wide field cosmological surveys.

We use a suite a cooling halo simulations to study a new mechanism for rapid
accretion of hot halo gas onto starforming galaxies. Correlated supernovae
events create converging 'superbubbles' in the halo gas. Where these collide,
the density increases, driving cooling filaments of low metallicity gas that
feed the disc. At our current numerical resolution (20 pc) we are only able to
resolve the most dramatic events; these could be responsible for the buildup
of galaxy discs after the most massive gasrich mergers have completed (z < 1).
As we increase the numerical resolution, we find that the filaments persist for
longer, driving continued latetime star formation. This suggests that
SNedriven accretion could act as an efficient mechanism for extracting cold
gas from the hot halo, driving latetime star formation in disc galaxies. We
show that such filament feeding leads to a peak star formation rate (SFR) of
$\sim 3$ M$_{\rm sun}$ yr$^{1}$, consistent with estimates for the Milky Way.
By contrast, direct cooling from the hot halo ('hotmode' accretion, not
present in the simulations that show filament feeding) falls short of the
SNedriven SFR by a factor of 34, and is sustained over a shorter time period.
The filaments we resolve extend to $\sim$ 50 kpc, reaching column densities of
$\sim 10^{18}$ cm$^{2}$. We show that such structures can plausibly explain
the broad dispersion in Mg II absorption seen along sight lines to quasars. Our
results suggest a dual role for stellar feedback in galaxy formation,
suppressing hotmode accretion while promoting coldmode accretion along
filaments. This ultimately leads to more star formation, suggesting that the
positive feedback effect outweighs the negative. Finally, since the filamentary
gas has higher angular momentum than that coming from hotmode accretion, we
show that this leads to the formation of substantially larger gas discs.

Spectroscopic redshift surveys offer great prospects for constraining the
dark sector in cosmology. Future surveys will however be both deep and wide and
will thus require an analysis in 3dimensional spherical geometry. We review
and compare several methods which have been proposed in the literature for this
purpose, focusing in particular on implementations of the spherical harmonic
tomography (SHT) power spectrum $C^{i j}_{l}$ and the spherical Fourier Bessel
(SFB) power spectrum $C_{l} (k, k')$. Using a Fisher analysis, we compare the
forecasted constraints on cosmological parameters using these statistics. These
constraints typically rely on approximations such as the Limber approximation
and make specific choices in the numerical implementation of each statistic.
Using a series of toy models, we explore the applicability of these
approximations and study the sensitivity of the SHT and SFB statistics to the
details of their implementation. In particular, we show that overlapping
redshift bins may improve cosmological constraints using the SHT statistic when
the number of bins is small, and that the SFB constraints are quite robust to
changes in the assumed distanceredshift relation. We also find that the SHT
can be tailored to be more sensitive to modes at redshifts close to the survey
boundary, while the SFB appears better suited to capture information beyond the
smooth shape of the power spectrum. In this context, we discuss the pros and
cons of the different techniques and their impact on the design and analysis of
future wide field spectroscopic surveys.