
Using Nbody simulations we study the impact of various systematic effects on
the bulk flow (BF) and the Cosmic Mach Number (CMN). We consider two types of
systematics: those related to survey properties and those induced by observer's
location in the Universe. In the former category we model sparse sampling,
velocity errors, and survey incompleteness. In the latter, we consider Local
Group (LG) analogue observers, placed in a specific location within the Cosmic
Web, satisfying various observational criteria. We differentiate such LG
observers from Copernican ones, who are at random locations. We report strong
systematic effects on the measured BF and CMN induced by sparse sampling,
velocity errors and radial incompleteness. For BF most of these effects exceed
10\% for scales $R\leq100h^{1}$Mpc. For CMN some of these systematics can be
catastrophically large ($>50\%$) also on bigger scales. Moreover, we find that
the position of the observer in the Cosmic Web significantly affects the
locally measured BF (CMN), with effects as large as $\sim20\%$ ($30\%)$ at
$R\leq50h^{1}$Mpc for a LGlike observer as compared to a random one. This
effect is comparable to the sample variance. To highlight the importance of
these systematics, we additionally study a model of modified gravity (MG) with
$\sim15\%$ enhanced growth rate. We found that the systematic effects can mimic
the modified gravity signal. The worstcase scenario is realized for a case of
a LGlike observer, when the effects induced by local structures are degenerate
with the enhanced growth rate fostered by MG. Our results indicate that
dedicated constrained simulations and realistic mock galaxy catalogs will be
absolutely necessary to fully benefit from the statistical power of the
forthcoming peculiar velocity data from surveys such as TAIPAN, WALLABY, Cosmic
Flows4 and SKA.

We study galaxy formation in sterile neutrino dark matter models that differ
signifi cantly from both cold and from `warm thermal relic' models. We use the
EAGLE code to carry out hydrodynamic simulations of the evolution of pairs of
galaxies chosen to resemble the Local Group, as part of the APOSTLE simulations
project. We compare cold dark matter (CDM) with two sterile neutrino models
with 7 keV mass: one, the warmest among all models of this mass (LA120) and the
other, a relatively cold case (LA10). We show that the lower concentration of
sterile neutrino subhalos compared to their CDM counterparts makes the inferred
inner dark matter content of galaxies like Fornax (or Magellanic Clouds) less
of an outlier in the sterile neutrino cosmologies. In terms of the galaxy
number counts the LA10 simulations are indistinguishable from CDM when one
takes into account halotohalo (or `simulationtosimulation') scatter. In
order for the LA120 model to match the number of Local Group dwarf galaxies, a
higher fraction of low mass haloes is required to form galaxies than is
predicted by the EAGLE simulations. As the census of the Local Group galaxies
nears completion, this population may provide a strong discriminant between
cold and warm dark matter models.

We elucidate the importance of the consistent treatment of gravitymodel
specific nonlinearities when estimating the growth of cosmological structures
from redshift space distortions (RSD). Within the context of standard
perturbation theory (SPT), we compare the predictions of two theoretical
templates with redshift space data from COLA (COmoving Lagrangian Acceleration)
simulations in the normal branch of DGP gravity (nDGP) and General Relativity
(GR). Using COLA for these comparisons is validated using a suite of full
Nbody simulations for the same theories. The two theoretical templates
correspond to the standard general relativistic perturbation equations and
those same equations modelled within nDGP. Gravitational clustering nonlinear
effects are accounted for by modelling the power spectrum up to one loop order
and redshift space clustering anisotropy is modelled using the Taruya,
Nishimichi and Saito (TNS) RSD model. Using this approach, we attempt to
recover the simulation's fiducial logarithmic growth parameter $f$. By
assigning the simulation data with errors representing an idealised survey with
a volume of $10\mbox{Gpc}^3/h^3$, we find the GR template is unable to recover
fiducial $f$ to within 1$\sigma$ at $z=1$ when we match the data up to $k_{\rm
max}=0.195h$/Mpc. On the other hand, the DGP template recovers the fiducial
value within $1\sigma$. Further, we conduct the same analysis for sets of mock
data generated for generalised models of modified gravity using SPT, where
again we analyse the GR template's ability to recover the fiducial value. We
find that for models with enhanced gravitational nonlinearity, the theoretical
bias of the GR template becomes significant for stage IV surveys. Thus, we show
that for the future large data volume galaxy surveys, the selfconsistent
modelling of nonGR gravity scenarios will be crucial in constraining theory
parameters.

The cosmic web is one of the most striking features of the distribution of
galaxies and dark matter on the largest scales in the Universe. It is composed
of dense regions packed full of galaxies, long filamentary bridges, flattened
sheets and vast low density voids. The study of the cosmic web has focused
primarily on the identification of such features, and on understanding the
environmental effects on galaxy formation and halo assembly. As such, a variety
of different methods have been devised to classify the cosmic web  depending
on the data at hand, be it numerical simulations, large sky surveys or other.
In this paper we bring twelve of these methods together and apply them to the
same data set in order to understand how they compare. In general these cosmic
web classifiers have been designed with different cosmological goals in mind,
and to study different questions. Therefore one would not {\it a priori} expect
agreement between different techniques however, many of these methods do
converge on the identification of specific features. In this paper we study the
agreements and disparities of the different methods. For example, each method
finds that knots inhabit higher density regions than filaments, etc. and that
voids have the lowest densities. For a given web environment, we find
substantial overlap in the density range assigned by each web classification
scheme. We also compare classifications on a halobyhalo basis; for example,
we find that 9 of 12 methods classify around a third of groupmass haloes (i.e.
$M_{\rm halo}\sim10^{13.5}h^{1}M_{\odot}$) as being in filaments. Lastly, so
that any future cosmic web classification scheme can be compared to the 12
methods used here, we have made all the data used in this paper public.

We use a pair of high resolution Nbody simulations implementing two dark
matter models, namely the standard cold dark matter (CDM) cosmogony and a warm
dark matter (WDM) alternative where the dark matter particle is a 1.5keV
thermal relic. We combine these simulations with the GALFORM semianalytical
galaxy formation model in order to explore differences between the resulting
galaxy populations. We use GALFORM model variants for CDM and WDM that result
in the same z=0 galaxy stellar mass function by construction. We find that most
of the studied galaxy properties have the same values in these two models,
indicating that both dark matter scenarios match current observational data
equally well. Even in underdense regions, where discrepancies in structure
formation between CDM and WDM are expected to be most pronounced, the galaxy
properties are only slightly different. The only significant difference in the
local universe we find is in the galaxy populations of "Local Volumes", regions
of radius 1 to 8Mpc around simulated Milky Way analogues. In such regions our
WDM model provides a better match to observed local galaxy number counts and is
five times more likely than the CDM model to predict subregions within them
that are as empty as the observed Local Void. Thus, a highly complete census of
the Local Volume and future surveys of void regions could provide constraints
on the nature of dark matter.

We describe and demonstrate the potential of a new and very efficient method
for simulating certain classes of modified gravity theories, such as the widely
studied $f(R)$ gravity models. High resolution simulations for such models are
currently very slow due to the highly nonlinear partial differential equation
that needs to be solved exactly to predict the modified gravitational force.
This nonlinearity is partly inherent, but is also exacerbated by the specific
numerical algorithm used, which employs a variable redefinition to prevent
numerical instabilities. The standard NewtonGaussSeidel iterative method used
to tackle this problem has a poor convergence rate. Our new method not only
avoids this, but also allows the discretised equation to be written in a form
that is analytically solvable. We show that this new method greatly improves
the performance and efficiency of $f(R)$ simulations. For example, a test
simulation with $512^3$ particles in a box of size $512 \, \mathrm{Mpc}/h$ is
now 5 times faster than before, while a Millenniumresolution simulation for
$f(R)$ gravity is estimated to be more than 20 times faster than with the old
method. Our new implementation will be particularly useful for running very
high resolution, largesized simulations which, to date, are only possible for
the standard model, and also makes it feasible to run large numbers of lower
resolution simulations for covariance analyses. We hope that the method will
bring us to a new era for precision cosmological tests of gravity.

We use a set of Nbody simulations employing a modified gravity (MG) model
with Vainshtein screening to study matter and halo hierarchical clustering. As
testcase scenarios we consider two normal branch DvaliGabadadzePorrati
(nDGP) gravity models with mild and strong growth rate enhancement. We study
higherorder correlation functions $\xi_n(R)$ up to $n=9$ and associated
hierarchical amplitudes $S_n(R)\equiv\xi_n(R)/\sigma(R)^{2n2}$. We find that
the matter PDFs are strongly affected by the fifthforce on scales up to
$50h^{1}$Mpc, and the deviations from GR are maximised at $z=0$. For reduced
cumulants $S_n$, we find that at small scales $R\leq10h^{1}$Mpc the MG is
characterised by lower values, with the deviation growing from $7\%$ in the
reduced skewness up to even $40\%$ in $S_5$. To study the halo clustering we
use a simple abundance matching and divide haloes into thee fixed number
density samples. The halo twopoint functions are weakly affected, with a
relative boost of the order of a few percent appearing only at the smallest
pair separations ($r\leq 5h^{1}$Mpc). In contrast, we find a strong MG signal
in $S_n(R)$'s, which are enhanced compared to GR. The strong model exhibits a
$>3\sigma$ level signal at various scales for all halo samples and in all
cumulants. In this context, we find that the reduced kurtosis to be an
especially promising cosmological probe of MG. Even the mild nDGP model leaves
a $3\sigma$ imprint at small scales $R\leq3h^{1}$Mpc, while the stronger model
deviates from a GRsignature at nearly all scales with a significance of
$>5\sigma$. Since the signal is persistent in all halo samples and over a range
of scales, we advocate that the reduced kurtosis estimated from galaxy
catalogues can potentially constitute a strong MGmodel discriminatory as well
as GR selfconsistency test.

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.

We assess the effect of the local large scale structure on the estimation of
twopoint statistics of the observed radial peculiar velocities of galaxies. A
large Nbody simulation is used to examine these statistics from the
perspective of random observers as well as "Local Group (LG)like" observers
conditioned to reside in an environment resembling the observed universe within
20 Mpc. The local environment systematically distorts the shape and amplitude
of velocity statistics with respect to ensembleaveraged measurements made by a
Copernican (random) observer. The Virgo cluster has the most significant
impact, introducing large systematic deviations in all the statistics. For a
simple "tophat" selection function, an idealized survey extending to $\sim
160h^{1}\,{\rm Mpc}$ or deeper is needed to completely mitigate the effects of
the local environment. Using shallower catalogues leads to systematic
deviations of the order of $50$ to $200\%$ depending on the scale considered.
For a flat redshift distribution similar to the one of the CosmicFlows3
survey, the deviations are even more prominent in both the shape and amplitude
at all separations considered $({\stackrel{<}{{}_\sim}} 100h^{1}\,{\rm Mpc})$.
Conclusions based on statistics calculated without taking into account the
impact of the local environment should be revisited.

We use the EAGLE galaxy formation simulation to study the effects of baryons
on the power spectrum of the total matter and dark matter distributions and on
the velocity fields of dark matter and galaxies. On scales
$k{\stackrel{>}{{}_\sim}} 4{h\,{\rm Mpc}^{1}}$ the effect of baryons on the
amplitude of the total matter power spectrum is greater than $1\%$. The
backreaction of baryons affects the density field of the dark matter at the
level of $\sim3\%$ on scales of $1\leq k/({h\,{\rm Mpc}^{1}})\leq 5$. The dark
matter velocity divergence power spectrum at
$k{\stackrel{<}{{}_\sim}}0.5{h\,{\rm Mpc}^{1}}$ is changed by less than $1\%$.
The 2D redshiftspace power spectrum is affected at the level of $\sim6\%$ at
$\vec{k}{\stackrel{>}{{}_\sim}} 1{h\,{\rm Mpc}^{1}}$ (for $\mu>0.5$), but
for $\vec{k}\leq 0.4{h\,{\rm Mpc}^{1}}$ it differs by less than $1\%$. We
report vanishingly small baryonic velocity bias for haloes: the peculiar
velocities of haloes with $M_{200}>3\times10^{11}{{\rm M}_{\odot}}$ (hosting
galaxies with $M_{*}>10^9{{\rm M}_{\odot}}$) are affected at the level of at
most $1~$km/s, which is negligible for $1\%$precision cosmology. We caution
that since EAGLE overestimates cluster gas fractions it may also underestimate
the impact of baryons, particularly for the total matter power spectrum.
Nevertheless, our findings suggest that for theoretical modelling of redshift
space distortions and galaxy velocitybased statistics, baryons and their
backreaction can be safely ignored at the current level of observational
accuracy. However, we confirm that the modelling of the total matter power
spectrum in weak lensing studies needs to include realistic galaxy formation
physics in order to achieve the accuracy required in the precision cosmology
era.

The cold dark matter (CDM) cosmological model unambigously predicts that a
large number of haloes should survive as subhaloes when they are accreted into
a larger halo. The CDM model would be ruled out if such substructures were
shown not to exist. By contrast, if the dark matter consists of warm particles
(WDM), then below a threshold mass that depends on the particle mass far fewer
substructures would be present. Finding subhaloes below a certain mass would
then rule out warm particle masses below some value. Strong gravitational
lensing provides a clean method to measure the subhalo mass function through
distortions in the structure of Einstein rings and giant arcs.Using mock
lensing observations constructed from highresolution Nbody simulations, we
show that measurements of approximately 100 strong lens systems with a
detection limit of $10^7 h^{1} M_{\odot}$ would clearly distinguish CDM from
WDM in the case where this consists of 7 keV sterile neutrinos such as those
that might be responsible for the 3.5 keV Xray emission line recently detected
in galaxies and clusters.

We introduce Copernicus Complexio (COCO), a highresolution cosmological
Nbody simulation of structure formation in the $\Lambda{\rm CDM}{}$ model.
COCO follows an approximately spherical region of radius $\sim 17.4h^{1}\,{\rm
Mpc}$ embedded in a much larger periodic cube that is followed at lower
resolution. The high resolution volume has a particle mass of
$1.135\times10^5h^{1}{\rm M}_{\odot}$ (60 times higher than the MillenniumII
simulation). COCO gives the dark matter halo mass function over eight orders of
magnitude in halo mass; it forms $\sim 60$ haloes of galactic size, each
resolved with about 10 million particles. We confirm the powerlaw character of
the subhalo mass function, $\bar{N}(>\mu)\propto\mu^{s}$, down to a reduced
subhalo mass $M_{sub}/M_{200}\equiv\mu=10^{6}$, with a bestfit powerlaw
index, $s=0.94$, for hosts of mass $\langle M_{200}\rangle=10^{12}h^{1}{\rm
M}_{\odot}$. The concentrationmass relation of COCO haloes deviates from a
single power law for masses $M_{200}<\textrm{a few}\times 10^{8}h^{1}{\rm
M}_{\odot}$, where it flattens, in agreement with results by SanchezConde et
al. The host mass invariance of the reduced maximum circular velocity function
of subhaloes, $\nu\equiv V_{max}/V_{200}$, hinted at in previous simulations,
is clearly demonstrated over five orders of magnitude in host mass. Similarly,
we find that the average, normalised radial distribution of subhaloes is
approximately universal (i.e. independent of subhalo mass), as previously
suggested by the Aquarius simulations of individual haloes. Finally, we find
that at fixed physical subhalo size, subhaloes in lower mass hosts typically
have lower central densities than those in higher mass hosts.

Selfconsistent ${\it N}$body simulations of modified gravity models are a
key ingredient to obtain rigorous constraints on deviations from General
Relativity using largescale structure observations. This paper provides the
first detailed comparison of the results of different ${\it N}$body codes for
the $f(R)$, DGP, and Symmetron models, starting from the same initial
conditions. We find that the fractional deviation of the matter power spectrum
from $\Lambda$CDM agrees to better than $1\%$ up to $k \sim 510~h/{\rm Mpc}$
between the different codes. These codes are thus able to meet the stringent
accuracy requirements of upcoming observational surveys. All codes are also in
good agreement in their results for the velocity divergence power spectrum,
halo abundances and halo profiles. We also test the quasistatic limit, which
is employed in most modified gravity ${\it N}$body codes, for the Symmetron
model for which the most significant nonstatic effects among the models
considered are expected. We conclude that this limit is a very good
approximation for all of the observables considered here.

We present results of analysis of the dark matter (DM) pairwise velocity
statistics in different Cosmic Web environments. We use the DM velocity and
density field from the Millennium 2 simulation together with the NEXUS+
algorithm to segment the simulation volume into voxels uniquely identifying one
of the four possible environments: nodes, filaments, walls or cosmic voids. We
show that the PDFs of the mean infall velocities $v_{12}$ as well as its
spatial dependence together with the perpendicular and parallel velocity
dispersions bear a significant signal of the largescale structure environment
in which DM particle pairs are embedded. The pairwise flows are notably colder
and have smaller mean magnitude in wall and voids, when compared to much denser
environments of filaments and nodes. We discuss on our results, indicating that
they are consistent with a simple theoretical predictions for pairwise motions
as induced by gravitational instability mechanism. Our results indicate that
the Cosmic Web elements are coherent dynamical entities rather than just
temporal geometrical associations. In addition it should be possible to
observationally test various Cosmic Web finding algorithms by segmenting
available peculiar velocity data and studying resulting pairwise velocity
statistics

We model the linear and nonlinear growth of large scale structure in the
Cubic Galileon gravity model, by running a suite of Nbody cosmological
simulations using the {\tt ECOSMOG} code. Our simulations include the
Vainshtein screening effect, which reconciles the Cubic Galileon model with
local tests of gravity. In the linear regime, the amplitude of the matter power
spectrum increases by $\sim 25%$ with respect to the standard $\Lambda$CDM
model today. The modified expansion rate accounts for $\sim 20%$ of this
enhancement, while the fifth force is responsible for only $\sim 5%$. This is
because the effective unscreened gravitational strength deviates from standard
gravity only at late times, even though it can be twice as large today. In the
nonlinear regime ($k \gtrsim 0.1 h\rm{Mpc}^{1}$), the fifth force leads to
only a modest increase ($\lesssim 8%$) in the clustering power on all scales
due to the very efficient operation of the Vainshtein mechanism. Such a strong
effect is typically not seen in other models with the same screening mechanism.
The screening also results in the fifth force increasing the number density of
halos by less than 10%, on all mass scales. Our results show that the screening
does not ruin the validity of linear theory on large scales which anticipates
very strong constraints from galaxy clustering data. We also show that, whilst
the model gives an excellent match to CMB data on small angular scales ($l
\gtrsim 50$), the predicted integrated SachsWolf effect is in tension with
Planck/WMAP results.

Numerical simulations in modified gravity have commonly been performed under
the quasistatic approximation  that is, by neglecting the effect of time
derivatives in the equation of motion of the scalar field that governs the
fifth force in a given modified gravity theory. To test the validity of this
approximation, we analyse the case of $f(R)$ gravity beyond this quasistatic
limit, by considering effects, if any, these terms have in the matter and
velocity divergence cosmic fields. To this end, we use the adaptive mesh
refinement code ECOSMOG to study three variants ($f_{R}= 10^{4}[$F4$],
10^{5}[$F5$]$ and $10^{6}[$F6$]$) of the HuSawicki $f(R)$ gravity model,
each of which refers to a different magnitude for the scalar field that
generates the fifth force. We find that for F4 and F5, which show stronger
deviations from standard gravity, a lowresolution simulation is enough to
conclude that time derivatives make a negligible contribution to the matter
distribution. The F6 model shows a larger deviation from the quasistatic
approximation, but one that diminishes when resimulated at higher resolution.
We therefore come to the conclusion that the quasistatic approximation is
valid for the most practical applications in $f(R)$ cosmologies.

We study the nonlinear growth of structure in nonlocal gravity models with
the aid of Nbody simulation and the spherical collapse and halo models. We
focus on a model in which the inversesquared of the d'Alembertian operator
acts on the Ricci scalar in the action. For fixed cosmological parameters, this
model differs from $\Lambda{\rm CDM}$ by having a lower latetime expansion
rate and an enhanced and timedependent gravitational strength ($\sim 6\%$
larger today). Compared to $\Lambda{\rm CDM}$ today, in the nonlocal model,
massive haloes are slightly more abundant (by $\sim 10\%$ at $M \sim 10^{14}
M_{\odot}/h$) and concentrated ($\approx 8\%$ enhancement over a range of mass
scales), but their linear bias remains almost unchanged. We find that the
ShethTormen formalism describes the mass function and halo bias very well,
with little need for recalibration of free parameters. The fitting of the halo
concentrations is however essential to ensure the good performance of the halo
model on small scales. For $k \gtrsim 1 h/{\rm Mpc}$, the amplitude of the
nonlinear matter and velocity divergence power spectra exhibits a modest
enhancement of $\sim 12\%$ to $15\%$, compared to $\Lambda{\rm CDM}$ today.
This suggests that this model might only be distinguishable from $\Lambda{\rm
CDM}$ by future observational missions. We point out that the absence of a
screening mechanism may lead to tensions with Solar System tests due to local
time variations of the gravitational strength, although this is subject to
assumptions about the local time evolution of background averaged quantities.

The velocity field of dark matter and galaxies reflects the continued action
of gravity throughout cosmic history. We show that the loworder moments of the
pairwise velocity distribution, $v_{12}$, are a powerful diagnostic of the laws
of gravity on cosmological scales. In particular, the projected lineofsight
galaxy pairwise velocity dispersion, $\sigma_{12}(r)$, is very sensitive to the
presence of modified gravity. Using a set of highresolution Nbody simulations
we compute the pairwise velocity distribution and its projected lineofsight
dispersion for a class of modified gravity theories: the chameleon \fR gravity
and Galileon gravity (cubic and quartic). The velocities of dark matter halos
with a wide range of masses would exhibit deviations from General Relativity at
the $(510)\sigma$ level. We examine strategies for detecting these deviations
in galaxy redshift and peculiar velocity surveys. If detected, this signature
would be a "smoking gun" for modified gravity.

We investigate the performance of semianalytical modelling of largescale
structure in Galileon gravity cosmologies using results from Nbody
simulations. We focus on the Cubic and Quartic Galileon models that provide a
reasonable fit to CMB, SNIa and BAO data. We demonstrate that the ShethTormen
mass function and linear halo bias can be calibrated to provide a very good fit
to our simulation results. We also find that the halo concentrationmass
relation is well fitted by a power law. The nonlinear matter power spectrum
computed in the halo model approach is found to be inaccurate in the mildly
nonlinear regime, but captures reasonably well the effects of the Vainshtein
screening mechanism on small scales. In the Cubic model, the screening
mechanism hides essentially all of the effects of the fifth force inside
haloes. In the case of the Quartic model, the screening mechanism leaves behind
residual modifications to gravity, which make the effective gravitational
strength timevarying and smaller than the standard value. Compared to normal
gravity, this causes a deficiency of massive haloes and leads to a weaker
matter clustering on small scales. For both models, we show that there are
realistic halo occupation distributions of Luminous Red Galaxies that can match
both the observed largescale clustering amplitude and the number density of
these galaxies.

We investigate how longrange scalar interactions affect the properties of
dark matter haloes. For doing so we employ the ReBEL model which implements an
additional interaction between dark matter particles. On the phenomenological
level this is equivalent to a modification of gravity. We analyse the
differences between five ReBEL models and $\Lambda$CDM using a series of high
resolution cosmological simulations. Emphasis is placed on investigating how
halo properties change in the presence of a fifth force. We report that the
density profile of ReBEL haloes is well described by the NFW profile but with
mean concentrations from $5\%$ to a few times higher than the standard
$\Lambda$CDM value. We also find a slight increase of the halo spin for haloes
more massive than $5\times10^{11}\M_{\odot}$, reflecting a higher rotational
support of those haloes due to scalar forces. In addition, the dark matter
haloes in our models are more spherical than their counterparts in
$\Lambda$CDM. The ReBEL haloes are also more virialised, with a large
difference from $\Lambda$CDM for strong fifth forces and a much smaller change
for weak scalar interactions.

We develop a numerical algorithm to solve the highorder nonlinear
derivativecoupling equation associated with the quartic Galileon model, and
implement it in a modified version of the RAMSES Nbody code to study the
effect of the Galileon field on the largescale matter clustering. The
algorithm is tested for several matter field configurations with different
symmetries, and works very well. This enables us to perform the first
simulations for a quartic Galileon model which provides a good fit to the
cosmic microwave background (CMB) anisotropy, supernovae and baryonic acoustic
oscillations (BAO) data. Our result shows that the Vainshtein mechanism in this
model is very efficient in suppressing the spatial variations of the scalar
field. However, the time variation of the effective Newtonian constant caused
by the curvature coupling of the Galileon field cannot be suppressed by the
Vainshtein mechanism. This leads to a significant weakening of the strength of
gravity in highdensity regions at late times, and therefore a weaker matter
clustering on small scales. We also find that without the Vainshtein mechanism
the model would have behaved in a completely different way, which shows the
crucial role played by nonlinearities in modified gravity theories and the
importance of performing selfconsistent Nbody simulations for these theories.

We present the relation between the genus in cosmology and the Betti numbers
for excursion sets of three and twodimensional smooth Gaussian random fields,
and numerically investigate the Betti numbers as a function of threshold level.
Betti numbers are topological invariants of figures that can be used to
distinguish topological spaces. In the case of the excursion sets of a
threedimensional field there are three possibly nonzero Betti numbers;
$\beta_0$ is the number of connected regions, $\beta_1$ is the number of
circular holes, and $\beta_2$ is the number of threedimensional voids. Their
sum with alternating signs is the genus of the surface of excursion regions. It
is found that each Betti number has a dominant contribution to the genus in a
specific threshold range. $\beta_0$ dominates the highthreshold part of the
genus curve measuring the abundance of high density regions (clusters).
$\beta_1$ dominates the genus near the median thresholds which measures the
topology of negatively curved isodensity surfaces, and $\beta_2$ corresponds
to the lowthreshold part measuring the void abundance. We average the Betti
number curves (the Betti numbers as a function of the threshold level) over
many realizations of Gaussian fields and find that both the amplitude and shape
of the Betti number curves depend on the slope of the power spectrum $n$ in
such a way that their shape becomes broader and their amplitude drops less
steeply than the genus as $n$ decreases. This behaviour contrasts with the fact
that the shape of the genus curve is fixed for all Gaussian fields regardless
of the power spectrum. Even though the Gaussian Betti number curves should be
calculated for each given power spectrum, we propose to use the Betti numbers
for better specification of the topology of large scale structures in the
universe.

We study the topology of the Megaparsec Cosmic Web in terms of the
scaledependent Betti numbers, which formalize the topological information
content of the cosmic mass distribution. While the Betti numbers do not fully
quantify topology, they extend the information beyond conventional cosmological
studies of topology in terms of genus and Euler characteristic. The richer
information content of Betti numbers goes along the availability of fast
algorithms to compute them.
For continuous density fields, we determine the scaledependence of Betti
numbers by invoking the cosmologically familiar filtration of sublevel or
superlevel sets defined by density thresholds. For the discrete galaxy
distribution, however, the analysis is based on the alpha shapes of the
particles. These simplicial complexes constitute an ordered sequence of nested
subsets of the Delaunay tessellation, a filtration defined by the scale
parameter, $\alpha$. As they are homotopy equivalent to the sublevel sets of
the distance field, they are an excellent tool for assessing the topological
structure of a discrete point distribution. In order to develop an intuitive
understanding for the behavior of Betti numbers as a function of $\alpha$, and
their relation to the morphological patterns in the Cosmic Web, we first study
them within the context of simple heuristic Voronoi clustering models.
Subsequently, we address the topology of structures emerging in the standard
LCDM scenario and in cosmological scenarios with alternative dark energy
content. The evolution and scaledependence of the Betti numbers is shown to
reflect the hierarchical evolution of the Cosmic Web and yields a promising
measure of cosmological parameters. We also discuss the expected Betti numbers
as a function of the density threshold for superlevel sets of a Gaussian random
field.

We use a suite of high resolution stateoftheart Nbody Dark Matter
simulations of chameleon $f(R)$ gravity to study the higher order volume
averaged correlation functions $\bar{\xi_n}$ together with the hierarchical
$n$th order correlation amplitudes $S_n=\bar{\xi}_n/\bar{\xi}_2^{n1}$ and
density distribution functions (PDF). We show that under the nonlinear
modifications of gravity the hierarchical scaling of the reduced cumulants is
preserved. This is however characterised by significant changes of both the
$\bar{\xi_n}$ and $S_n$'s values and their scale dependence with respect to
General Relativity gravity (GR). In addition, we measure a significant increase
of the non linear $\sigma_8$ parameter reaching $14, 5$ and 0.5% in excess of
the GR value for the three flavours of our $f(R)$ models. We further note that
the values of the reduced cumulants up to order $n=9$ are significantly
increased in $f(R)$ gravity for all our models at small scales $R\simlt
30\hmpc$. In contrast the values of the hierarchical amplitudes, $S_n$'s, are
smaller in $f(R)$ indicating that the modified gravity density distribution
functions are deviating from the GR case. Furthermore we find that the redshift
evolution of relative deviations of the $f(R)$ hierarchical correlation
amplitudes is fastest at high and moderate redshifts $1\leq z \leq4$. The
growth of these deviations significantly slows down in the low redshift
universe. We also compute the PDFs and show that for scales below $\sim
20\hmpc$ they are significantly shifted in $f(R)$ gravity towards the low
densities. Finally we discuss the implications of our theoretical predictions
for measurements of the hierarchical clustering in galaxy redshifts surveys,
including the important problems of the galaxy biasing and redshifts space
distortions.