The nature of the dark matter can affect the collapse time of dark matter haloes, and can therefore be imprinted in observables such as the stellar population ages and star formation histories of dwarf galaxies. In this paper we use high resolution hydrodynamical simulations of Local Group-analogue (LG) volumes in cold dark matter (CDM), sterile neutrino warm dark matter (WDM) and self-interacting dark matter (SIDM) models with the EAGLE galaxy formation code to study how galaxy formation times change with dark matter model. We are able to identify the same haloes in different simulations, since they share the same initial density field phases. We show that the stellar mass varies systematically with resolution by over a factor of two, in a manner that depends on the final stellar mass. The evolution of the stellar populations in SIDM is largely identical to that of CDM, but in WDM early star formation is instead suppressed. The time at which LG haloes can begin to form stars through atomic cooling is delayed by $\sim$200~Myr in WDM models compared to CDM. 70~per~cent of WDM haloes of mass $>10^{8}M_{\odot}$ collapse early enough to form stars before $z=6$, compared to 90~per~cent of CDM and SIDM galaxies. It will be necessary to measure stellar ages of old populations to a precision of better than 100~Myr, and to address degeneracies with the redshift of reionization, in order to use these observables to distinguish between dark matter models.

We study the reduced skewness, $S_{3,g}\equiv\bar{\xi}_{3,g}/\bar{\xi}_{2,g}^2$ of galaxy distribution at $z=0.5$ in two families of modfied gravity models: the Hu-Sawicki $f(R)$-gravity and normal-branch of Dvali-Gabadadze-Porrati (nDGP) models. We use a set of mock galaxy catalogues specifally designed to match CMASS spectroscopic galaxy sample. For the first time we investigate the third reduced moment of such galaxy distributions both in the redshift space. Our analysis confirms that the signal previously indicated only for dark matter halo catalogues persists also in realistic mock galaxy samples. This result offers a possibility to extract a potential modified gravity signal in $S_3$ from spectroscopic galaxy data without a need for a very precise and self-consistent RSD models constructed for each and every modified gravity scenario separately. We show that the relative deviations from $\Lambda$CDM~ $S_{3,g}$ of various modified gravity models can vary from $7$ down to $\sim 2-3\%$ effects. Albeit, the effect looks small, we show that for considered models it can foster a $2-3\sigma$ falsification. Finally we argue that galaxy sample of a significantly higher number density should provide even stronger constraints by limiting shot-noise effects affecting the $S_{3,g}$ estimates at small comoving separations.

The nearby Perseus galaxy cluster is a key target for indirect detection searches for decaying dark matter. We use the C-EAGLE simulations of galaxy clusters to predict the flux, width and shape of a dark matter decay line, paying particular attention to the unexplained 3.55keV line detected in the spectra of some galaxies and clusters, and the upcoming XRISM X-ray observatory mission. We show that the line width in C-EAGLE clusters similar to Perseus is typically [600-800]$\mathrm{kms^{-1}}$, and therefore narrower than the amplitude of the velocity dispersion of galaxies in the cluster. Haloes that are significantly disturbed can, however, exhibit galaxy velocity dispersions higher than $1000\mathrm{kms^{-1}}$, and in this case will show a large difference between the line profiles of on- and off-center observations. We show that the line profile is likely to be slightly asymmetric, but still well approximated by a Gaussian at the 10% level, and that the halo asymmetry can lead to fluxes that vary by a factor of two. In summary, we predict that, if the previously reported 3.55keV line detections do originate from dark matter decay, the XRISM mission will detect a line with a roughly Gaussian profile at a rest frame energy of 3.55keV, with a width $>600\mathrm{kms^{-1}}$ and flux approximately in the range $[4-9]\times10^{-8}\mathrm{counts/sec/cm^{2}}$.

Dark matter particles may decay, emitting photons. Drawing on the EAGLE family of hydrodynamic simulations of galaxy formation -- including the APOSTLE and C-EAGLE simulations -- we assess the systematic uncertainties and scatter on the decay flux from different galaxy classes, from Milky Way satellites to galaxy clusters, and compare our results to studies of the 3.55~keV line. We demonstrate that previous detections and non-detections of this line are consistent with a dark matter interpretation. For example, in our simulations the width of the the dark matter decay line for Perseus-analogue galaxy clusters lies in the range 1300-1700~\kms. Therefore, the non-detection of the 3.55~keV line in the centre of the Perseus cluster by the {\it Hitomi} collaboration is consistent with detections by other instruments. We also consider trends with stellar and halo mass and evaluate the scatter in the expected fluxes arising from the anisotropic halo mass distribution and from object-to-object variations. We provide specific predictions for observations with {\it XMM-Newton} and with the planned X-ray telescopes {\it XRISM} and {\it ATHENA}. If future detections of unexplained X-ray lines match our predictions, including line widths, we will have strong evidence that we have discovered the dark matter.

Using N-body 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 LG-like 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 worst-case scenario is realized for a case of a LG-like 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 Flows-4 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 halo-to-halo (or `simulation-to-simulation') 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 gravity-model specific non-linearities 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 N-body 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 non-linear 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 non-linearity, 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 self-consistent modelling of non-GR 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 halo-by-halo basis; for example, we find that 9 of 12 methods classify around a third of group-mass 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 N-body 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 semi-analytical 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 under-dense 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 sub-regions 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 Newton-Gauss-Seidel 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 Millennium-resolution 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, large-sized 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 N-body simulations employing a modified gravity (MG) model with Vainshtein screening to study matter and halo hierarchical clustering. As test-case scenarios we consider two normal branch Dvali-Gabadadze-Porrati (nDGP) gravity models with mild and strong growth rate enhancement. We study higher-order correlation functions $\xi_n(R)$ up to $n=9$ and associated hierarchical amplitudes $S_n(R)\equiv\xi_n(R)/\sigma(R)^{2n-2}$. We find that the matter PDFs are strongly affected by the fifth-force 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 two-point 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 GR-signature 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 MG-model discriminatory as well as GR self-consistency test.

We present two-point correlation function statistics of the mass and the halos in the chameleon $f(R)$ modified gravity scenario using a series of large volume N-body simulations. Three distinct variations of $f(R)$ are considered (F4, F5 and F6) and compared to a fiducial $\Lambda$CDM model in the redshift range $z \in [0,1]$. We find that the matter clustering is indistinguishable for all models except for F4, which shows a significantly steeper slope. The ratio of the redshift- to real-space correlation function at scales $> 20 h^{-1} \mathrm{Mpc}$ agrees with the linear General Relativity (GR) Kaiser formula for the viable $f(R)$ models considered. We consider three halo populations characterized by spatial abundances comparable to that of luminous red galaxies (LRGs) and galaxy clusters. The redshift-space halo correlation functions of F4 and F5 deviate significantly from $\Lambda$CDM at intermediate and high redshift, as the $f(R)$ halo bias is smaller or equal to that of the $\Lambda$CDM case. Finally we introduce a new model independent clustering statistic to distinguish $f(R)$ from GR: the relative halo clustering ratio -- $\mathcal{R}$. The sampling required to adequately reduce the scatter in $\mathcal{R}$ will be available with the advent of the next generation galaxy redshift surveys. This will foster a prospective avenue to obtain largely model-independent cosmological constraints on this class of modified gravity models.

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree 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.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area 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 magnitude-limited 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.

We use the Copernicus Complexio (COCO) high resolution $N$-body simulations to investigate differences in the properties of small-scale structures in the standard cold dark matter (CDM) model and in a model with a cutoff in the initial power spectrum of density fluctuations consistent with both a thermally produced warm dark matter (WDM) particle or a sterile neutrino with mass 7 keV and leptogenesis parameter $L_6=8.7$. The latter corresponds to the "coldest" model with this sterile neutrino mass compatible with the identification of the recently detected 3.5 keV X-ray line as resulting from particle decay. CDM and WDM predict very different number densities of subhaloes with mass $\leq 10^9\,h^{-1}\,M_\odot$ although they predict similar, nearly universal, normalised subhalo radial density distributions. Haloes and subhaloes in both models have cuspy NFW profiles, but WDM subhaloes below the cutoff scale in the power spectrum (corresponding to maximum circular velocities $V_{\mathrm{max}}^{z=0} \leq50~\mathrm{kms}^{-1}$) are less concentrated than their CDM counterparts. We make predictions for observable properties using the GALFORM semi-analytic model of galaxy formation. Both models predict Milky Way satellite luminosity functions consistent with observations, although the WDM model predicts fewer very faint satellites. This model, however, predicts slightly more UV bright galaxies at redshift $z>7$ than CDM, but both are consistent with observations. Gravitational lensing offers the best prospect of distinguishing between the models.

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.

We assess the effect of the local large scale structure on the estimation of two-point statistics of the observed radial peculiar velocities of galaxies. A large N-body 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 ensemble-averaged 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 "top-hat" 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 CosmicFlows-3 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 a suite of cosmological simulations to study the mass-concentration-redshift relation, $c({\rm M},z)$, of dark matter halos. Our simulations include standard $\Lambda$-cold dark matter (CDM) models, and additional runs with truncated power spectra, consistent with a thermal warm dark matter (WDM) scenario. We find that the mass profiles of CDM and WDM halos are self-similar and well approximated by the Einasto profile. The $c({\rm M},z)$ relation of CDM halos is monotonic: concentrations decrease with increasing virial mass at fixed redshift, and decrease with increasing redshift at fixed mass. The mass accretion histories (MAHs) of CDM halos are also scale-free, and can be used to infer concentrations directly. These results do not apply to WDM halos: their MAHs are not scale-free because of the characteristic scale imposed by the power-spectrum suppression. Further, the WDM $c({\rm M},z)$ relation is non-monotonic: concentrations peak at a mass scale dictated by the truncation scale, and decrease at higher and lower masses. We show that the assembly history of a halo can still be used to infer its concentration, provided that the total mass of its progenitors is considered (the "collapsed mass history"; CMH), rather than just that of its main ancestor. This exploits the scale-free nature of CMHs to derive a simple scaling that reproduces the mass-concentration-redshift relation of both CDM and WDM halos over a vast range of halo masses and redshifts. Our model therefore provides a robust account of the mass, redshift, cosmology and power spectrum dependence of dark matter halo concentrations.

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 back-reaction 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 redshift-space 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 velocity-based statistics, baryons and their back-reaction 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 high-resolution N-body 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 X-ray emission line recently detected in galaxies and clusters.

We introduce Copernicus Complexio (COCO), a high-resolution cosmological N-body 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 Millennium-II 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 power-law 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 best-fit power-law index, $s=0.94$, for hosts of mass $\langle M_{200}\rangle=10^{12}h^{-1}{\rm M}_{\odot}$. The concentration-mass 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 Sanchez-Conde 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.

Self-consistent ${\it N}$-body simulations of modified gravity models are a key ingredient to obtain rigorous constraints on deviations from General Relativity using large-scale 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 5-10~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 quasi-static limit, which is employed in most modified gravity ${\it N}$-body codes, for the Symmetron model for which the most significant non-static effects among the models considered are expected. We conclude that this limit is a very good approximation for all of the observables considered here.

The detection of planar structures within the satellite systems of both the Milky Way (MW) and Andromeda (M31) has been reported as being in stark contradiction to the predictions of the standard cosmological model ($\Lambda$CDM). Given the ambiguity in defining a planar configuration, it is unclear how to interpret the low incidence of the MW and M31 planes in $\Lambda$CDM. We investigate the prevalence of satellite planes around galactic mass haloes identified in high resolution cosmological simulations. We find that planar structures are very common, and that ~10% of $\Lambda$CDM haloes have even more prominent planes than those present in the Local Group. While ubiquitous, the planes of satellite galaxies show a large diversity in their properties. This precludes using one or two systems as small scale probes of cosmology, since a large sample of satellite systems is needed to obtain a good measure of the object-to-object variation. This very diversity has been misinterpreted as a discrepancy between the satellite planes observed in the Local Group and $\Lambda$CDM predictions. In fact, ~10% of $\Lambda$CDM galactic haloes have planes of satellites that are as infrequent as the MW and M31 planes. The look-elsewhere effect plays an important role in assessing the detection significance of satellite planes and accounting for it leads to overestimating the significance level by a factor of 30 and 100 for the MW and M31 systems, respectively.

The recent detection of a 3.5 keV X-ray line from the centres of galaxies and clusters by Bulbul et al. (2014a) and Boyarsky et al. (2014a) has been interpreted as emission from the decay of 7 keV sterile neutrinos which could make up the (warm) dark matter (WDM). As part of the COpernicus COmplexio (COCO) programme, we investigate the properties of dark matter haloes formed in a high-resolution cosmological $N$-body simulation from initial conditions similar to those expected in a universe in which the dark matter consists of 7 keV sterile neutrinos. This simulation and its cold dark matter (CDM) counterpart have $\sim13.4$bn particles, each of mass $\sim 10^5\, h^{-1} M_\odot$, providing detailed information about halo structure and evolution down to dwarf galaxy mass scales. Non-linear structure formation on small scales ($M_{200}\, \leq\, 2 \times 10^9\,h^{-1}\,M_\odot$) begins slightly later in COCO-Warm than in COCO-Cold. The halo mass function at the present day in the WDM model begins to drop below its CDM counterpart at a mass $\sim 2 \times 10^{9}\,h^{-1}\,M_\odot$ and declines very rapidly towards lower masses so that there are five times fewer haloes of mass $M_{200}= 10^{8}\,h^{-1}\,M_\odot$ in COCO-Warm than in COCO-Cold. Halo concentrations on dwarf galaxy scales are correspondingly smaller in COCO-Warm, and we provide a simple functional form that describes its evolution with redshift. The shapes of haloes are similar in the two cases, but the smallest haloes in COCO-Warm rotate slightly more slowly than their CDM counterparts.

We investigate the properties of dark matter haloes and subhaloes in an $f(R)$ gravity model with $|f_{R0}|=10^{-6}$, using a very high-resolution N-body simulation. The model is a borderline between being cosmologically interesting and yet still consistent with current data. We find that the halo mass function in this model has a maximum 20% enhancement compared with the $\Lambda$CDM predictions between $z=1$ and $z=0$. Because of the chameleon mechanism which screens the deviation from standard gravity in dense environments, haloes more massive than $10^{13}h^{-1}M_\odot$ in this $f(R)$ model have very similar properties to haloes of similar mass in $\Lambda$CDM, while less massive haloes, such as that of the Milky Way, can have steeper inner density profiles and higher velocity dispersions due to their weaker screening. The halo concentration is remarkably enhanced for low-mass haloes in this model due to a deepening of the total gravitational potential. Contrary to the naive expectation, the halo formation time $z_f$ is later for low-mass haloes in this model, a consequence of these haloes growing faster than their counterparts in $\Lambda$CDM at late times and the definition of $z_f$. Subhaloes, especially those less massive than $10^{11}h^{-1}M_\odot$, are substantially more abundant in this $f(R)$ model for host haloes less massive than $10^{13}h^{-1}M_\odot$. We discuss the implications of these results for the Milky Way satellite abundance problem. Although the overall halo and subhalo properties in this borderline $f(R)$ model are close to their $\Lambda$CDM predictions, our results suggest that studies of the Local Group and astrophysical systems, aided by high-resolution simulations, can be valuable for further tests of it.

We apply a semi-analytic galaxy formation model to two high resolution cosmological N-body simulations to investigate analogues of the Milky Way system. We select these according to observed properties of the Milky Way rather than by halo mass as in most previous work. For disk-dominated central galaxies with stellar mass (5--7) x 10d10Msun, the median host halo mass is 1.4 x 10d12Msun, with 1 sigma dispersion in the range [0.86, 3.1] x 10d12Msun, consistent with dynamical measurements of the Milky Way halo mass. For any given halo mass, the probability of hosting a Milky Way system is low, with a maximum of ~20% in haloes of mass ~10d12Msun. The model reproduces the V-band luminosity function and radial profile of the bright (MV < -9) Milky Way satellites. Galaxy formation in low mass haloes is found to be highly stochastic, resulting in an extremely large scatter in the relation between MV (or stellar mass) for satellites and the depth of the subhalo potential well in which they live, as measured by the maximum of the rotation curve, Vmax. We conclude that the "too big to fail" problem is an artifact of selecting satellites in N-body simulations according to subhalo properties: in 10% of cases we find that three or fewer of the brightest (or most massive) satellites have Vmax > 30 km/s. Our model predicts that around half of the dark matter subhaloes with Vmax > 20 km/s host satellites fainter than MV = -9 and so may be missing from existing surveys.

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 large-scale 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 N-body 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 Sachs-Wolf effect is in tension with Planck/WMAP results.

Numerical simulations in modified gravity have commonly been performed under the quasi-static 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 quasi-static 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 Hu-Sawicki $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 low-resolution 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 quasi-static approximation, but one that diminishes when re-simulated at higher resolution. We therefore come to the conclusion that the quasi-static approximation is valid for the most practical applications in $f(R)$ cosmologies.

We use the distribution of maximum circular velocities, $V_{max}$, of satellites in the Milky Way (MW) to constrain the virial mass, $M_{200}$, of the Galactic halo under an assumed prior of a $\Lambda$CDM universe. This is done by analysing the subhalo populations of a large sample of halos found in the Millennium II cosmological simulation. The observation that the MW has at most three subhalos with $V_{max}\ge30 km/s$ requires a halo mass $M_{200}\le1.4\times10^{12} M_\odot$, while the existence of the Magellanic Clouds (assumed to have $V_{max}\ge60 km/s$) requires $M_{200}\ge1.0\times10^{12} M_\odot$. The first of these conditions is necessary to avoid the "too-big-to-fail" problem highlighted by Boylan-Kolchin et al., while the second stems from the observation that massive satellites like the Magellanic Clouds are rare. When combining both requirements, we find that the MW halo mass must lie in the range $0.25 \le M_{200}/(10^{12} M_\odot) \le 1.4$ at $90\%$ confidence. The gap in the abundance of Galactic satellites between $30 km/s\le V_{max} \le 60 km/s$ places our galaxy in the tail of the expected satellite distribution.

We study the substructure population of Milky Way (MW)-mass halos in the $\Lambda$CDM cosmology using a novel procedure to extrapolate subhalo number statistics beyond the resolution limit of N-body simulations. The technique recovers the mean and the variance of the subhalo abundance, but not its spatial distribution. It extends the dynamic range over which precise statistical predictions can be made by the equivalent of performing a simulation with 50 times higher resolution, at no additional computational cost. We apply this technique to MW-mass halos, but it can easily be applied to halos of any mass. We find up to $20\%$ more substructures in MW-mass halos than found in previous studies. Our analysis lowers the mass of the MW halo required to accommodate the observation that the MW has only three satellites with a maximum circular velocity $V_{max}\ge30 km/s$ in the $\Lambda$CDM cosmology. The probability of having a subhalo population similar to that in the MW is $20\%$ for a virial mass, $M_{200}=1\times10^{12} M_\odot$ and practically zero for halos more massive than $M_{200}=2\times10^{12} M_\odot$.

We study the nonlinear growth of structure in nonlocal gravity models with the aid of N-body simulation and the spherical collapse and halo models. We focus on a model in which the inverse-squared 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 late-time expansion rate and an enhanced and time-dependent 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 Sheth-Tormen 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 low-order moments of the pairwise velocity distribution, $v_{12}$, are a powerful diagnostic of the laws of gravity on cosmological scales. In particular, the projected line-of-sight galaxy pairwise velocity dispersion, $\sigma_{12}(r)$, is very sensitive to the presence of modified gravity. Using a set of high-resolution N-body simulations we compute the pairwise velocity distribution and its projected line-of-sight 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 $(5-10)\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 semi-analytical modelling of large-scale structure in Galileon gravity cosmologies using results from N-body 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 Sheth-Tormen 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 concentration-mass 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 time-varying 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 large-scale clustering amplitude and the number density of these galaxies.

We investigate how long-range 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.

Comparison of peculiar velocities of galaxies with their gravitational accelerations (induced by the density field) is one of the methods to constrain the redshift distortion parameter \beta=(\Omega_m^0.55)/b, where \Omega_m is the non-relativistic matter density parameter and b is the linear bias. In particular, one can use the motion of the Local Group (LG) for that purpose. Its peculiar velocity is known from the dipole component of the cosmic microwave background, whereas its acceleration can be estimated with the use of an all-sky galaxy catalog, from the so-called clustering dipole. At the moment, the biggest dataset of that kind is the Two Micron All Sky Survey Extended Source Catalog (2MASS XSC) containing almost 1 million galaxies and complete up to ~300 Mpc/h. We applied 2MASS data to measure LG acceleration and used two methods to estimate the beta parameter. Both of them yield \beta~0.4 with an error of several per cent, which is the most precise determination of this parameter from the clustering dipole to date.

We develop a numerical algorithm to solve the high-order nonlinear derivative-coupling equation associated with the quartic Galileon model, and implement it in a modified version of the RAMSES N-body code to study the effect of the Galileon field on the large-scale 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 high-density 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 self-consistent N-body simulations for these theories.

We present the relation between the genus in cosmology and the Betti numbers for excursion sets of three- and two-dimensional 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 three-dimensional field there are three possibly non-zero 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 three-dimensional 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 high-threshold 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 iso-density surfaces, and $\beta_2$ corresponds to the low-threshold 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 scale-dependent 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 scale-dependence 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 scale-dependence 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 state-of-the-art N-body 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^{n-1}$ and density distribution functions (PDF). We show that under the non-linear 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.

In this contribution we present the preliminary results regarding the non-linear BAO signal in higher-order statistics of the cosmic density field. We use ensembles of N-body simulations to show that the non-linear evolution changes the amplitudes of the BAO signal, but has a negligible effect on the scale of the BAO feature. The latter observation accompanied by the fact that the BAO feature amplitude roughly doubles as one moves to higher orders, suggests that the higher-order correlation amplitudes can be used as probe of the BAO signal.

On 28th and 29th of January 2013 we held an international meeting in Zielona G\'ora (Poland) honouring the first anniversary of premature passing away^(1) of Professor Roman Juszkiewicz. We have celebrated an opening of a new seminar room at the University of Zielona G\'ora commemorated to the memory of Roman Juszkiewicz and we have shared our anecdotes and memories of this great scientist and friend. Here we want to present a limited and short memoir and obituary for Roman Juszkiewicz. 1: Roman Juszkiewicz passed away on 28th of January 2012

In this study we show that the skewness S_3 of the cosmic density field contains a significant and potentially detectable and clean imprint of Baryonic Acoustic Oscillations. Although the BAO signal in the skewness has a lower amplitude than second order measures like the two-point correlation function and power spectrum, it has the advantage of a considerably lower sensitivity to systematic influences. Because it lacks a direct dependence on bias if this concerns simple linear bias, skewness will be considerably less beset by uncertainties due to galaxy bias. Also, it has a weaker sensitivity to redshift distortion effects. We use perturbation theory to evaluate the magnitude of the effect on the volume-average skewness, for various cosmological models. One important finding of our analysis is that the skewness BAO signal occurs at smaller scales than that in second order statistics. For an LCDM spectrum with WMAP7 normalization, the BAO feature has a maximum wiggle amplitude of ~3% and appears at a scale of ~82Mpc/h. We conclude that the detection of BAO wiggles in future extensive galaxy surveys via the skewness of the observed galaxy distribution may provide us with a useful, and potentially advantageous, measure of the nature of Dark Energy.

One of the important unknowns of current cosmology concerns the effects of the large scale distribution of matter on the formation and evolution of dark matter haloes and galaxies. One main difficulty in answering this question lies in the absence of a robust and natural way of identifying the large scale environments and their characteristics. This work summarizes the NEXUS+ formalism which extends and improves our multiscale scale-space MMF method. The new algorithm is very successful in tracing the Cosmic Web components, mainly due to its novel filtering of the density in logarithmic space. The method, due to its multiscale and hierarchical character, has the advantage of detecting all the cosmic structures, either prominent or tenuous, without preference for a certain size or shape. The resulting filamentary and wall networks can easily be characterized by their direction, thickness, mass density and density profile. These additional environmental properties allows to us to investigate not only the effect of environment on haloes, but also how it correlates with the environment characteristics.

We study the matter and velocity divergence power spectra in a f(R) gravity theory and their time evolution measured from several large-volume N-body simulations with varying box sizes and resolution. We find that accurate prediction of the matter power spectrum in f(R) gravity places stronger requirements on the simulation than is the case with LCDM, because of the nonlinear nature of the fifth force. Linear perturbation theory is shown to be a poor approximation for the f(R) models, except when the chameleon effect is very weak. We show that the relative differences from the fiducial LCDM model are much more pronounced in the nonlinear tail of the velocity divergence power spectrum than in the matter power spectrum, which suggests that future surveys which target the collection of peculiar velocity data will open new opportunities to constrain modified gravity theories. A close investigation of the time evolution of the power spectra shows that there is a pattern in the evolution history, which can be explained by the properties of the chameleon-type fifth force in f(R) gravity. Varying the model parameter |f_R0|, which quantifies the strength of the departure from standard gravity, mainly varies the epoch marking the onset of the fifth force, as a result of which the different f(R) models are in different stages of the same evolutionary path at any given time

We introduce a new descriptor of the weblike pattern in the distribution of galaxies and matter: the scale dependent 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 used in earlier analyses of cosmological models. The richer information content of Betti numbers goes along with the availability of fast algorithms to compute them. When measured as a function of scale they provide a "Betti signature" for a point distribution that is a sensitive yet robust discriminator of structure. The signature is highly effective in revealing differences in structure arising in different cosmological models, and is exploited towards distinguishing between different dark energy models and may likewise be used to trace primordial non-Gaussianities. In this study we demonstrate the potential of Betti numbers by studying their behaviour in simulations of cosmologies differing in the nature of their dark energy.

In this work we investigate higher order statistics for the $\lcdm$ and ReBEL scalar-interacting dark matter models by analyzing $180\hmpc$ dark matter N-body simulation ensembles. The N-point correlation functions and the related hierarchical amplitudes, such as skewness and kurtosis, are computed using the Count-In-Cells method. Our studies demonstrate that the hierarchical amplitudes $S_n$ of the scalar-interacting dark matter model significantly deviate from the values in the $\lcdm$ cosmology on scales comparable and smaller then the screening length $r_s$ of a given scalar-interacting model. The corresponding additional forces that enhance the total attractive force exerted on dark matter particles at galaxy scales lowers the values of the hierarchical amplitudes $S_n$. We conclude that hypothetical additional exotic interactions in the dark matter sector should leave detectable markers in the higher-order correlation statistics of the density field. We focussed in detail on the redshift evolution of the dark matter field's skewness and kurtosis. From this investigation we find that the deviations from the canonical $\lcdm$ model introduced by the presence of the ``fifth'' force attain a maximum value at redshifts $0.5<z<2$. We therefore conclude that moderate redshift data are better suited for setting observational constraints on the investigated ReBEL models.

We investigate the effect of long-range scalar interactions in dark matter (DM) models of cosmic structure formation with a particular focus on the formation times of haloes. Utilising $N$-body simulations with $512^3$ DM particles we show that in our models dark matter haloes form substantially earlier: tracing objects up to redshift $z\sim6$ we find that the formation time, as characterised by the redshift $z_{1/2}$ at which the halo has assembled half of its final mass, is gradually shifted from $z_{1/2}\approx 1.83$ in the fiducial LCDM model to $z_{1/2}\approx 2.54$ in the most extreme self-interaction model. This is accompanied by a shift of the redshift that marks the transition between merger and steady accretion epochs from $z_{*}\approx 4.32$ in the \lcdm\ halos to $z_{*}\approx 6.39$ in our strongest interaction model. In other words, the scalar-interacting model employed in this work produces more structures at high redshifts, prolonging at the same time the steady accretion phases. These effects taken together can help the LCDM model to account for a high redshift reionisation as indicated by the WMAP data and can alleviate issues related to the survival of the thin-disk dominated galaxies at low redshifts.

Based on a set of cosmological N-body simulations we analyze properties of the dark matter haloes (DM) in a galaxy mass range ($10^{11} - 10^{13} h^{-1}M_{\odot}$) in modified $\lcdm$ cosmology with additional dynamically screened scalar interactions in DM sector. Our simulations show that scalar interactions support picture of the Island Universe. Rapid structure formation processes are shifted into higher redshifts resulting in a much smaller accretion and merging rates for galactic haloes at low redshifts. Finally, we present how this "fifth" force affects halo properties, like density profile, triaxiality, ellipticities and the spin parameter.

We explore the possibility of improving the $\Lambda$CDM model at megaparsec scales by introducing a scalar interaction which increases the mutual gravitational attraction of dark matter particles. Using N-body simulations, we study the spatial distribution of dark matter particles and halos. We measure the effect of modifications in the Newton's gravity on properties of the two-point correlation function, the dark matter power spectrum, the cumulative halo mass function and density probability distribution functions. The results look promising: the scalar interactions produce desirable features at megaparsec scales without spoiling the $\Lambda$CDM successes at larger scales.

In this paper we present preliminary results from cosmological simulations of modified gravity in the dark matter sector. Our results show improvements over standard cold dark matter cosmology. The abundance of low-mass haloes in the modified gravity model fit observations better than the conventional theory, while the differences of the modified density fluctuation power spectrum differs from the standard, $\Lambda$CDM power spectrum are small enough to make these two models observationally indistinguishable at large scales.