
Reconstructing the expansion history of the Universe from type Ia supernovae
data, we fit the growth rate measurements and put modelindependent constraints
on some key cosmological parameters, namely, $\Omega_\mathrm{m},\gamma$, and
$\sigma_8$. The constraints are consistent with those from the concordance
model within the framework of general relativity, but the current quality of
the data is not sufficient to rule out modified gravity models. Adding the
condition that dark energy density should be positive at all redshifts,
independently of its equation of state, further constrains the parameters and
interestingly supports the concordance model.

The form of the primordial power spectrum (PPS) of cosmological scalar
(matter density) perturbations is not yet constrained satisfactorily in spite
of the tremendous amount of information from the Cosmic Microwave Background
(CMB) data. While a smooth powerlawlike form of the PPS is consistent with
the CMB data, some PPS with small nonsmooth features at large scales can also
fit the CMB temperature and polarization data with similar statistical
evidence. Future CMB surveys cannot help distinguish all such models due to the
cosmic variance at large angular scales. In this paper, we study how well we
can differentiate be tween such featured forms of the PPS not otherwise
distinguishable using CMB data. We ran 15 Nbody DESIlike simulations of these
models to explore this approach. Showing that statistics such as the halo mass
function and the twopoint correlation function are not able to distinguish
these models in a DESIlike survey, we advocate to avoid reducing the
dimensionality of the problem by demonstrating that the use of a simple
threedimensional countincell density field can be much more effective for
the purpose of model distinction.

Large deviation statistics is implemented to predict the statistics of cosmic
densities in cylinders applicable to photometric surveys. It yields few percent
accurate analytical predictions for the onepoint probability distribution
function (PDF) of densities in concentric or compensated cylinders; and also
captures the densitydependence of their angular clustering (cylinder bias).
All predictions are found to be in excellent agreement with the cosmological
simulation Horizon Run 4 in the quasilinear regime where standard perturbation
theory normally breaks down. These results are combined with a simple local
bias model that relates dark matter and tracer densities in cylinders and
validated on simulated halo catalogues. This formalism can be used to probe
cosmology with existing and upcoming photometric surveys like DES, Euclid or
WFIRST containing billions of galaxies.

A flat FriedmanRobersonWalker universe dominated by a cosmological constant
($\Lambda$) and cold dark matter (CDM) has been the working model preferred by
cosmologists since the discovery of cosmic acceleration. However, tensions of
various degrees of significance are known to be present among existing datasets
within the $\Lambda$CDM framework. In particular, the Lyman$\alpha$ forest
measurement of the Baryon Acoustic Oscillations (BAO) by the Baryon Oscillation
Spectroscopic Survey (BOSS) prefers a smaller value of the matter density
fraction $\Omega_{\rm M}$ compared to the value preferred by cosmic microwave
background (CMB). Also, the recently measured value of the Hubble constant,
$H_0=73.24\pm1.74 \ {\rm km}\ {\rm s}^{1} \ {\rm Mpc}^{1}$, is $3.4\sigma$
higher than $66.93\pm0.62 \ {\rm km}\ {\rm s}^{1} \ {\rm Mpc}^{1}$ inferred
from the Planck CMB data. In this work, we investigate if these tensions can be
interpreted as evidence for a nonconstant dynamical dark energy (DE). Using
the KullbackLeibler (KL) divergence to quantify the tension between datasets,
we find that the tensions are relieved by an evolving DE, with the dynamical DE
model preferred at a $3.5\sigma$ significance level based on the improvement in
the fit alone. While, at present, the Bayesian evidence for the dynamical DE is
insufficient to favour it over $\Lambda$CDM, we show that, if the current best
fit DE happened to be the true model, it would be decisively detected by the
upcoming DESI survey.

Starting from a very accurate model for densityincells statistics of dark
matter based on large deviation theory, a bias model for the tracer density in
spheres is formulated. It adopts a mean bias relation based on a quadratic bias
model to relate the logdensities of dark matter to those of massweighted dark
haloes in real and redshift space. The validity of the parametrised bias model
is established using a parametrisationindependent extraction of the bias
function. This average bias model is then combined with the dark matter PDF,
neglecting any scatter around it: it nevertheless yields an excellent model for
densitiesincells statistics of mass tracers that is parametrised in terms of
the underlying dark matter variance and three bias parameters. The procedure is
validated on measurements of both the one and two point statistics of subhalo
densities in the stateoftheart Horizon Run 4 simulation showing excellent
agreement for measured dark matter variance and bias parameters. Finally, it is
demonstrated that this formalism allows for a joint estimation of the nonlinear
dark matter variance and the bias parameters using solely the statistics of
subhaloes. Having verified that galaxy counts in hydrodynamical simulations
sampled on a scale of 10 Mpc/h closely resemble those of subhaloes, this work
provides important steps towards making theoretical predictions for
densityincells statistics applicable to upcoming galaxy surveys like Euclid
or WFIRST.

We study the properties of dark matter haloes in a wide range of modified
gravity models, namely, $f(R)$, DGP, and interacting dark energy models. We
study the effects of modified gravity and dark energy on the internal
properties of haloes, such as the spin and the structural parameters. We find
that $f(R)$ gravity enhance the median value of the Bullock spin parameter, but
could not detect such effects for DGP and coupled dark energy. $f(R)$ also
yields a lower median sphericity and oblateness, while coupled dark energy has
the opposite effect. However, these effects are very small. We then study the
interaction rate of haloes in different gravity, and find that only strongly
coupled dark energy models enhance the interaction rate. We then quantify the
enhancement of the alignment of the spins of interacting halo pairs by modified
gravity. Finally, we study the alignment of the major axes of haloes with the
largescale structures. The alignment of the spins of interacting pairs of
haloes in DGP and coupled dark energy models show no discrepancy with GR, while
$f(R)$ shows a weaker alignment. Strongly coupled dark energy shows a stronger
alignment of the halo shape with the largescale structures.

Using measurements of $H(z)$ and $d_\mathrm{A}(z)$ from the Baryon
Oscillation Spectroscopic Survey (BOSS) DR12 and luminosity distances from the
Joint Lightcurve Analysis (JLA) compilation of supernovae (SN), we measure $H_0
r_\mathrm{d}$ without any model assumption. Our measurement of $H_0
r_\mathrm{d} = $($10033.20^{+333.10}_{371.81}$ (SN) $\pm$ 128.19 (BAO)) km
s$^{1}$ is consistent with Planck constrains for the flat {\Lambda}CDM model.
We also report that higher expansion history rates $h(z)$ (among the
possibilities) as well as lowerbound values of $H_0 r\mathrm{d}$ result in
better internal consistency among the independent data (H(z)rd and
$d_\mathrm{A}(z)/r_\mathrm{d}$ from BAO at $z=0.32$ and z=0.57 and
$d_\mathrm{L}$ from JLA) we used in this work. This can be interpreted as an
interesting and independent support of Planck cosmology without using any
cosmic microwave background data. We then combine these observables to test the
FriedmannLema\^itreRobertsonWalker (FLRW) metric and the flatness of the
Universe in a modelindependent way at two redshifts, namely 0.32 and 0.57, by
introducing a new diagnostic for flatFLRW, $\Theta(z)$, which only depends on
observables of BAO and SN data. Our results are consistent with a flatFLRW
Universe within $2{\sigma}$.

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.

Interactions such as mergers and flybys play a fundamental role in shaping
galaxy morphology. Using the Horizon Run 4 cosmological Nbody simulation, we
studied the frequency and type of halo interactions, and their redshift
evolution as a function of the environment defined by the largescale density,
pair separation, mass ratio, and target halo mass. Most interactions happen at
largescale density contrast $\delta \approx 20$, regardless of the redshift,
corresponding to groups and relatively dense part of filaments. However, the
fraction of interacting target is maximum at $\delta \approx 1000$. We provide
a new empirical fitting form for the interaction rate as a function of the halo
mass, largescale density, and redshift. We also report the existence of two
modes of interactions from the distributions of mass ratio and relative
distance, implying two different physical origins of the interaction. Satellite
targets lose their mass as they proceed deeper into the host halo. The relative
importance of these two trends strongly depends on the largescale density,
target mass, and redshift.

Galaxies accrete their mass by means of both smooth accretion from the cosmic
web, and the mergers of smaller entities. We wish to quantify the respective
role of these two modes of accretion, which could determine the morphological
types of galaxies observed today. Multizoom cosmological simulations are used
to estimate as a function of time the evolution of mass in bound systems, for
dark matter as well as baryons. The baryonic contents of dark matter haloes are
studied. Merger histories are followed as a function of external density, and
the different ways in which mass is assembled in galaxies and the stellar
component accumulated are quantified. We find that most galaxies assemble their
mass through smooth accretion, and only the most massive galaxies also grow
significantly through mergers. The mean fraction of mass assembled by accretion
is 77 %, and by mergers 23 %. We present typical accretion histories of
hundreds of galaxies: masses of the most massive galaxies increase
monotonically in time, mainly through accretion, many intermediatemass objects
also experience massloss events such as tidal stripping and evaporation.
However, our simulations suffer from the overcooling of massive galaxies caused
by the neglect of active galaxy nuclei (AGN) feedback. The time by which half
of the galay mass has assembled, both in dark matter and baryons, is a
decreasing function of mass, which is compatible with the observations of a
socalled downsizing. At every epoch in the universe, there are lowmass
galaxies actively forming stars, while more massive galaxies form their stars
over a shorter period of time within half the age of the universe.

According to the hierarchical model, small galaxies form first and merge
together to form bigger objects. In parallel, galaxies assemble their mass
through accretion from cosmic filaments. Recently, the increased spatial
resolution of the cosmological simulations have emphasised that a large
fraction of cold gas can be accreted by galaxies. In order to compare the role
of both phenomena and the corresponding star formation history, one has to
detect the structures in the numerical simulations and to follow them in time,
by building a merger tree.