Within the Halo Model of large scale structure, all matter is contained in dark matter halos. This simple yet powerful framework has been broadly applied to multiple data sets and enriched our comprehension of how matter is distributed in the Universe. In this work we extend this assumption by allowing for matter to rest not only inside halos but also within cosmic voids and in between halos and voids (which we call 'dust'). This assumption leads to additional contributions (1Void, 2Void, Halo-Void, etc.) to the predictions of correlation functions, spectra and profiles for both halos and voids. Whereas the Halo Model can only make predictions for halo quantities, the Halo Void Model extends those for void statistics and halo-void cross-correlations. We provide recipes for all new ingredients of the Halo Void (Dust) Model, such as the void abundance, linear bias and density profile and test their validity in a N-body simulation. Including voids and dust into the calculations improves the transition between the 1Halo and the 2Halo terms by up to $\sim 6\%$. It also eliminates the need to include low-mass structures on the normalization of large-scale terms, suggesting that halos and voids are complementary cosmic structures to effectively describe matter distribution on large scales of the Universe.

Modified gravity (MG) theories aim to reproduce the observed acceleration of the Universe by reducing the dark sector while simultaneously recovering General Relativity (GR) within dense environments. Void studies appear to be a suitable scenario to search for imprints of alternative gravity models on cosmological scales. Voids cover an interesting range of density scales where screening mechanisms fade out, which reaches from a density contrast $\delta \approx -1$ close to their centers to $\delta \approx 0$ close to their boundaries. We present an analysis of the level of distinction between GR and two modified gravity theories, the Hu-Sawicki $f(R)$ and the symmetron theory. This study relies on the abundance, linear bias, and density profile of voids detected in n-body cosmological simulations. We define voids as connected regions made up of the union of spheres with a {\it \textup{mean}} density given by $\overline\rho_v=0.2\,\overline\rho_m$, but disconnected from any other voids. We find that the height of void walls is considerably affected by the gravitational theory, such that it increases for stronger gravity modifications. Finally, we show that at the level of dark matter n-body simulations, our constraints allow us to distinguish between GR and MG models with $|f_{R0}| > 10^{-6}$ and $z_{SSB} > 1$. Differences of best-fit values for MG parameters that are derived independently from multiple void probes may indicate an incorrect MG model. This serves as an important consistency check.

We develop a new, simple, fast and parameter-free method to construct dark matter halo catalogues. This method requires as inputs only the linear matter power spectrum and the threshold density for halo formation in linear theory. It directly uses excursion set ideas and Lagrangian perturbation theory to produce halo catalogues with the correct abundance, large scale power spectrum, bispectrum and velocity field. These halo catalogues can be used for the fast construction of mock galaxy catalogues, allowing for the evaluation of covariance matrices for multiple observables. Because of its robustness and predictive nature, this method can be easily adapted to produce catalogues with e.g. primordial non-Gaussianities, modified theories of gravity and non-standard dark energy models, enabling detailed studies of these models in the context of next-generation surveys. We implement this method in a C code, and present numerical comparisons with theoretical predictions as well as full N-body cosmological simulations.

We investigate the relationship between the turnaround radius ($R_t$) and the virial mass ($M_v$) of cosmic structures in the context of $\Lambda$CDM model and in an $f(R)$ model of modified gravity -- namely, the Hu-Sawicki model. The $R_t$ is the distance from the center of the cosmic structure to the shell that is detaching from the Hubble flow at a given time, while the $M_v$ is defined, for this work, as the mass enclosed within the volume where the density is $200$ times the background density. We consider that gravitationally bound astrophysical systems follow a Navarro-Frenk-White density profile, while beyond the virial radius ($R_v$) the profile is approximated by the 2-halo term of the matter correlation function. By combining them together with the information drawn from solving the spherical collapse for the structures, we are able to connect two observables: the $R_t$ and the $M_v$. We show that, in $\Lambda$CDM, the turnaround mass ($M_t$) at $z=0$ is related to the $M_v$ of that same structure by $M_t \simeq 3.07 \, M_v$, while in terms of the radii we have that $R_t \simeq 3.7 \, R_v$ (for $M_v$ of $10^{13} \, h^{-1} \, M_\odot$). In the $f(R)$ model, we have $M_t \simeq 3.43 \, M_v$ and $R_t \simeq 4.1 \, R_v$, for $|f_{R0}|=10^{-6}$ and the same mass scale. Therefore, the difference between $\Lambda$CDM and $f(R)$ in terms of these observable relations is of order $\sim 10-20\%$ even for a relatively mild strength of the modification of gravity ($|f_{R0}|=10^{-6}$). For the $R_t$ itself we find a difference of $\sim 9\%$ between the weakly modification in gravity considered in this work ($|f_{R0}|=10^{-6}$) and $\Lambda$CDM for a mass of $10^{13} \, h^{-1} \, M_\odot$. Once observations allow precisions of this order or better in measurements $R_t$, as well as the $M_v$, these quantities will become powerful tests of modified gravity.

We investigate the turnaround radius in the spherical collapse model, both in General Relativity and in modified gravity, in particular $f(R)$ scenarios. The phases of spherical collapse are marked by the density contrast in the instant of turnaround $\delta_t$, and by the linear density contrast in the moment of collapse, $\delta_c$. We find that the effective mass of the extra scalar degree of freedom which arises in modified gravity models has an impact on $\delta_t$ of up to $\sim10\%$, and that $\delta_c$ can increase by $\sim1.0\%$. We also compute the turnaround radius, $R_t$, which in modified gravity models can increase by up to $\sim 6\%$ at $z \simeq 0$.

We use a spherical model and an extended excursion set formalism with drifting diffusive barriers to predict the abundance of cosmic voids in the context of general relativity as well as f(R) and symmetron models of modified gravity. We detect spherical voids from a suite of N-body simulations of these gravity theories and compare the measured void abundance to theory predictions. We find that our model correctly describes the abundance of both dark matter and galaxy voids, providing a better fit than previous proposals in the literature based on static barriers. We use the simulation abundance results to fit for the abundance model free parameters as a function of modified gravity parameters, and show that counts of dark matter voids can provide interesting constraints on modified gravity. For galaxy voids, more closely related to optical observations, we find that constraining modified gravity from void abundance alone may be significantly more challenging. In the context of current and upcoming galaxy surveys, the combination of void and halo statistics including their abundances, profiles and correlations should be effective in distinguishing modified gravity models that display different screening mechanisms.