Using 324 numerically modelled galaxy clusters we investigate the radial and galaxy-halo alignment of dark matter subhaloes and satellite galaxies orbiting within and around them. We find that radial alignment depends on distance to the centre of the galaxy cluster but appears independent of the dynamical state of the central host cluster. Furthermore, we cannot find a relation between radial alignment of the halo or galaxy shape with its own mass. We report that backsplash galaxies, i.e. objects that have already passed through the cluster radius but are now located in the outskirts, show a stronger radial alignment than infalling objects. We further find that there exists a population of well radially aligned objects passing very close to the central cluster's centre which were found to be on highly radial orbit.

We introduce a new multi-power-law distribution for the Initial Mass Function (IMF) to explore its potential properties. It follows on prior work that introduced mechanisms accounting for mass accretion in star formation, developed within the framework of general evolution equations for the mass distribution of accreting and non-accreting (proto)stars. This paper uses the same fundamental framework to demonstrate that the interplay between a mass-dependent and a time-dependent step-like dropout rate from accretion leads to IMFs that exhibit multiple power laws for an exponential mass growth. While the mass-dependent accretion and its dropout is intrinsic to each star, the time-dependent dropout might be tied to a specific history such as the rapid consumption of nebular material by nearby stars or the sweeping away of some material by shock waves. The time-dependent dropout folded into the mass-dependent process of star formation is shown to have a significant influence on the IMFs.

The three-point correlation function (3PCF) of the 21cm brightness temperature from the Epoch of Reionization (EoR) probes complementary information to the commonly studied two-point correlation function (2PCF) about the morphology of ionized regions. We investigate the 21cm 2PCF and 3PCF in configuration space using semi-numerical simulations and test whether they can be described by the local quadratic bias model. We find that fits of bias model predictions for the 2PCF and 3PCF deviate from our measurements by $\sim 20\%$ at scales above the typical size of ionized regions ($\simeq 30$ Mpc) and at early times with global neutral fractions of $\langle x_{\rm HI} \rangle \gtrsim 0.7$. At later times and smaller scales these deviations increase strongly, indicating a break down of the bias model. The 2PCF and 3PCF fits of the linear bias parameter agree at the $10\%$ level for different EoR model configurations. This agreement holds, when adding redshift space distortions to the simulations. The relation between spatial fluctuations in the matter density and the 21cm signal, as predicted by the bias model, is consistent with direct measurements of this relation in simulations for large smoothing scales ($\gtrsim 30$ Mpc). From this latter test we conclude that negative amplitudes of the 21cm 3PCF result from negative bias parameters, which describe the anti-correlation between the matter over-densities and the 21cm signal during the EoR. However, a more detailed interpretation of the bias parameters may require a description of non-local contributions to the bias model.

We introduce a new dual power law (DPL) probability distribution function for the mass distribution of stellar and substellar objects at birth, otherwise known as the initial mass function (IMF). The model contains both deterministic and stochastic elements, and provides a unified framework within which to view the formation of brown dwarfs and stars resulting from an accretion process that starts from extremely low mass seeds. It does not depend upon a top down scenario of collapsing (Jeans) masses or an initial lognormal or otherwise IMF-like distribution of seed masses. Like the modified lognormal power law (MLP) distribution, the DPL distribution has a power law at the high mass end, as a result of exponential growth of mass coupled with equally likely stopping of accretion at any time interval. Unlike the MLP, a power law decay also appears at the low mass end of the IMF. This feature is closely connected to the accretion stopping probability rising from an initially low value up to a high value. This might be associated with physical effects of ejections sometimes (i.e., rarely) stopping accretion at early times followed by outflow driven accretion stopping at later times, with the transition happening at a critical time (therefore mass). Comparing the DPL to empirical data, the critical mass is close to the substellar mass limit, suggesting that the onset of nuclear fusion plays an important role in the subsequent accretion history of a young stellar object.

We compare reduced three-point correlations $Q$ of matter, haloes (as proxies for galaxies) and their cross correlations, measured in a total simulated volume of $\sim100 \ (h^{-1} \text{Gpc})^{3}$, to predictions from leading order perturbation theory on a large range of scales in configuration space. Predictions for haloes are based on the non-local bias model, employing linear ($b_1$) and non-linear ($c_2$, $g_2$) bias parameters, which have been constrained previously from the bispectrum in Fourier space. We also study predictions from two other bias models, one local ($g_2=0$) and one in which $c_2$ and $g_2$ are determined by $b_1$ via an approximately universal relation. Overall, measurements and predictions agree when $Q$ is derived for triangles with $(r_1r_2r_3)^{1/3} \gtrsim 60 h^{-1}\text{Mpc}$, where $r_{1-3}$ are the sizes of the triangle legs. Predictions for $Q_{matter}$, based on the linear power spectrum, show significant deviations from the measurements at the BAO scale (given our small measurement errors), which strongly decrease when adding a damping term or using the non-linear power spectrum, as expected. Predictions for $Q_{halo}$ agree best with measurements at large scales when considering non-local contributions. The universal bias model works well for haloes and might therefore be also useful for tightening constraints on $b_1$ from $Q$ in galaxy surveys. Such constraints are independent of the amplitude of matter density fluctuation ($\sigma_8$) and hence break the degeneracy between $b_1$ and $\sigma_8$, present in galaxy two-point correlations.

We study the linear and non-linear bias parameters which determine the mapping between the distributions of galaxies and the full matter density fields, comparing different measurements and predictions. Associating galaxies with dark matter haloes in the MICE Grand Challenge N-body simulation we directly measure the bias parameters by comparing the smoothed density fluctuations of haloes and matter in the same region at different positions as a function of smoothing scale. Alternatively we measure the bias parameters by matching the probability distributions of halo and matter density fluctuations, which can be applied to observations. These direct bias measurements are compared to corresponding measurements from two-point and different third-order correlations, as well as predictions from the peak-background model, which we presented in previous articles using the same data. We find an overall variation of the linear bias measurements and predictions of $\sim 5 \%$ with respect to results from two-point correlations for different halo samples with masses between $\sim 10^{12} - 10^{15}$ $h^{-1}M_\odot$ at the redshifts $z=0.0$ and $0.5$. Variations between the second- and third-order bias parameters from the different methods show larger variations, but with consistent trends in mass and redshift. The various bias measurements reveal a tight relation between the linear and the quadratic bias parameters, which is consistent with results from the literature based on simulations with different cosmologies. Such a universal relation might improve constraints on cosmological models, derived from second-order clustering statistics at small scales or higher-order clustering statistics.

Using dark matter simulations we show how halo bias is determined by local density and not by halo mass. This is not totally surprising, as according to the peak-background split model, local density is the property that constraints bias at large scales. Massive haloes have a high clustering because they reside in high density regions. Small haloes can be found in a wide range of environments which determine their clustering amplitudes differently. This contradicts the assumption of standard Halo Occupation Distribution (HOD) models that the bias and occupation of haloes is determined solely by their mass. We show that the bias of central galaxies from semi-analytic models of galaxy formation as a function of luminosity and colour is not correctly predicted by the standard HOD model. Using local density instead of halo mass the HOD model correctly predicts galaxy bias. These results indicate the need to include information about local density and not only mass in order to correctly apply HOD analysis in these galaxy samples. This new model can be readily applied to observations and has the advantage that the galaxy density can be directly observed, in contrast with the dark matter halo mass.

We study halo clustering bias with second- and third-order statistics of halo and matter density fields in the MICE Grand Challenge simulation. We verify that two-point correlations deliver reliable estimates of the linear bias parameters at large scales, while estimations from the variance can be significantly affected by non-linear and possibly non-local contributions to the bias function. Combining three-point auto- and cross-correlations we find, for the first time in configuration space, evidence for the presence of such non-local contributions. These contributions are consistent with predicted second-order non-local effects on the bias functions originating from the dark matter tidal field. Samples of massive haloes show indications of bias (local or non-local) beyond second order. Ignoring non-local bias causes $20-30$\% and $5-10$\% overestimation of the linear bias from three-point auto- and cross-correlations respectively. We study two third-order bias estimators which are not affected by second-order non-local contributions. One is a combination of three-point auto- and cross- correlation. The other is a combination of third-order one- and two-point cumulants. Both methods deliver accurate bias estimations of the linear bias. Furthermore their estimations of second-order bias agree mutually. Ignoring non-local bias causes higher values of the second-order bias from three-point correlations. Our results demonstrate that third-order statistics can be employed for breaking the growth-bias degeneracy.

We model the abundance of haloes in the $\sim(3 \ \text{Gpc}/h)^3$ volume of the MICE Grand Challenge simulation by fitting the universal mass function with an improved Jack-Knife error covariance estimator that matches theory predictions. We present unifying relations between different fitting models and new predictions for linear ($b_1$) and non-linear ($c_2$ and $c_3$) halo clustering bias. Different mass function fits show strong variations in their performance when including the low mass range ($M_h \lesssim 3 \ 10^{12} \ M_{\odot}/h$) in the analysis. Together with fits from the literature we find an overall variation in the amplitudes of around $10$% in the low mass and up to $50$% in the high mass (galaxy cluster) range ($M_h > 10^{14} \ M_{\odot}/h$). These variations propagate into a $10$% change in $b_1$ predictions and a $50$% change in $c_2$ or $c_3$. Despite these strong variations we find universal relations between $b_1$ and $c_2$ or $c_3$ for which we provide simple fits. Excluding low mass haloes, different models fitted with reasonable goodness in this analysis, show percent level agreement in their $b_1$ predictions, but are systematically $5-10$% lower than the bias directly measured with two-point halo-mass clustering. This result confirms previous findings derived from smaller volumes (and smaller masses). Inaccuracies in the bias predictions lead to $5-10$% errors in growth measurements. They also affect any HOD fitting or (cluster) mass calibration from clustering measurements.

Measurements of the linear growth factor $D$ at different redshifts $z$ are key to distinguish among cosmological models. One can estimate the derivative $dD(z)/d\ln(1+z)$ from redshift space measurements of the 3D anisotropic galaxy two-point correlation $\xi(z)$, but the degeneracy of its transverse (or projected) component with galaxy bias $b$, i.e. $\xi_{\perp}(z) \propto\ D^2(z) b^2(z)$, introduces large errors in the growth measurement. Here we present a comparison between two methods which break this degeneracy by combining second- and third-order statistics. One uses the shape of the reduced three-point correlation and the other a combination of third-order one- and two-point cumulants. These methods use the fact that, for Gaussian initial conditions and scales larger than $20$ $h^{-1}$Mpc, the reduced third-order matter correlations are independent of redshift (and therefore of the growth factor) while the third-order galaxy correlations depend on $b$. We use matter and halo catalogs from the MICE-GC simulation to test how well we can recover $b(z)$ and therefore $D(z)$ with these methods in 3D real space. We also present a new approach, which enables us to measure $D$ directly from the redshift evolution of second- and third-order galaxy correlations without the need of modelling matter correlations. For haloes with masses lower than $10^{14}$ $h^{-1}$M$_\odot$, we find $10%$ deviations between the different estimates of $D$, which are comparable to current observational errors. At higher masses we find larger differences that can probably be attributed to the breakdown of the bias model and non-Poissonian shot noise.

We study the shapes of subhalo distributions from four dark-matter-only simulations of Milky Way type haloes. Comparing the shapes derived from the subhalo distributions at high resolution to those of the underlying dark matter fields we find the former to be more triaxial if theanalysis is restricted to massive subhaloes. For three of the four analysed haloes the increased triaxiality of the distributions of massive subhaloes can be explained by a systematic effect caused by the low number of objects. Subhaloes of the fourth halo show indications for anisotropic accretion via their strong triaxial distribution and orbit alignment with respect to the dark matter field. These results are independent of the employed subhalo finder. Comparing the shape of the observed Milky Way satellite distribution to those of high-resolution subhalo samples from simulations, we find an agreement for samples of bright satellites, but significant deviations if faint satellites are included in the analysis. These deviations might result from observational incompleteness.

A quantitative study of the clustering properties of the cosmic web as a function of absolute magnitude and colour is presented using the SDSS Data Release 7 galaxy survey. Mark correlations are included in the analysis. We compare our results with mock galaxy samples obtained with four different semi-analytical models of galaxy formation imposed on the merger trees of the Millenium simulation. The clustering of both red and blue galaxies is studied separately.