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.

The collapse of a collisionless self-gravitating system, with the fast achievement of a quasi-stationary state, is driven by violent relaxation, with a typical particle interacting with the time-changing collective potential. It is traditionally assumed that this evolution is governed by the Vlasov-Poisson equation, in which case entropy must be conserved. We run N-body simulations of isolated self-gravitating systems, using three simulation codes: NBODY-6 (direct summation without softening), NBODY-2 (direct summation with softening) and GADGET-2 (tree code with softening), for different numbers of particles and initial conditions. At each snapshot, we estimate the Shannon entropy of the distribution function with three different techniques: Kernel, Nearest Neighbor and EnBiD. For all simulation codes and estimators, the entropy evolution converges to the same limit as N increases. During violent relaxation, the entropy has a fast increase followed by damping oscillations, indicating that violent relaxation must be described by a kinetic equation other than the Vlasov-Poisson, even for N as large as that of astronomical structures. This indicates that violent relaxation cannot be described by a time-reversible equation, shedding some light on the so-called "fundamental paradox of stellar dynamics". The long-term evolution is well described by the orbit-averaged Fokker-Planck model, with Coulomb logarithm values in the expected range 10-12. By means of NBODY-2, we also study the dependence of the 2-body relaxation time-scale on the softening length. The approach presented in the current work can potentially provide a general method for testing any kinetic equation intended to describe the macroscopic evolution of N-body systems.

In cosmology, phenomenologically motivated expressions for running vacuum are commonly parametrized as linear functions $\Lambda(H^2)$ or $\Lambda(R)$. Such kind of models assume an equation of state for vacuum given by $\,\overline P_\Lambda=-\,\overline\rho_\Lambda$, relating their background pressure $\,\overline P_\Lambda$ and mean energy density $\,\overline\rho_\Lambda\equiv\Lambda/8\pi G$. This equation of state requires that the dynamic for vacuum is due to the energy exchange with the material species. Most of the approaches to background level consider only the energy exchange between vacuum and the transient dominant material component of the universe. We extend such models assuming the running vacuum as the sum of independent contributions $\,\overline\rho_{\Lambda} =\sum_i\,\overline\rho_{\Lambda i}$, associated with (and interacting with) each of the $i$ material species. We derive the linear scalar perturbations for two running scenarios, modeling its cosmic evolution and identifying their different imprints on the cosmic microwave background anisotropies and the matter power spectrum. In the $\Lambda(H^2)$ scenario the running vacuum are coupled with all the material species in the universe, whereas the $\Lambda(R)$ description only leads to coupling between vacuum and the non-relativistic matter components; which produces different imprints of the two models on the matter power spectrum. A comparison with the Planck 2015 data was made in order to constrain the free parameters of the models. In the case of the $\Lambda(H^2)$ model, it was found that $\Omega_\Lambda=0.705\pm0.027$ and $H_0=69.6\pm2.9\, km\, Mpc^{-1}\, s^{-1}$, which diminish the tension with the low redshift expectations.

We propose a large class of nonsingular cosmologies of arbitrary spatial curvature whose cosmic history is determined by a primeval dynamical $\Lambda (t)$-term. For all values of the curvature, the models evolve between two extreme de Sitter phases driven by the relic time-varying vacuum energy density. The transition from inflation to the radiation phase is universal and points to a natural solution of the graceful exit problem regardless of the values of the curvature parameter. The flat case recovers the scenario recently discussed in the literature (Perico et al., Phys. Rev. D88, 063531, 2013). The early de Sitter phase is characterized by an arbitrary energy scale $H_I$ associated to the primeval vacuum energy density. If $H_I$ is fixed to be nearly the Planck scale, the ratio between the relic and the present observed vacuum energy density is $\rho_{vI}/\rho_{v0} \simeq 10^{123}$.

In the present mainstream cosmology, matter and spacetime emerged from a singularity and evolved through four distinct periods: early inflation, radiation, dark matter and late-time inflation (driven by dark energy). During the radiation and dark matter dominated stages, the universe is decelerating while the early and late-time inflations are accelerating stages. A possible connection between the accelerating periods remains unknown, and, even more intriguing, the best dark energy candidate powering the present accelerating stage (Lambda-vacuum) is plagued with the cosmological constant and coincidence puzzles. Here we propose an alternative solution for such problems based on a large class of time-dependent vacuum energy density models in the form of power series of the Hubble rate, Lambda=Lambda(H). The proposed class of Lambda(H)-decaying vacuum model provides: i) a new mechanism for inflation (different from the usual inflaton models), (ii) a natural mechanism for a graceful exit, which is universal for the whole class of models; iii) the currently accelerated expansion of the universe, iv) a mild dynamical dark energy at present; and v) a final de Sitter stage. Remarkably, the late-time cosmic expansion history of our class of models is very close to the concordance LambdaCDM model, but above all it furnishes the necessary smooth link between the initial and final de Sitter stages through the radiation- and matter-dominated epochs.

Modifications on the predictions about the matter power spectrum based on the hypothesis of a tiny contribution from a degenerate Fermi gas (DFG) test-fluid to some dominant cosmological scenario are investigated. Reporting about the systematic way of accounting for all the cosmological perturbations, through the Boltzmann equation we obtain the analytical results for density fluctuation, $\delta$, and fluid velocity divergence, $\theta$, of the DFG. Small contributions to the matter power spectrum are analytically obtained for the radiation-dominated background, through an ultra-relativistic approximation, and for the matter-dominated and $\Lambda$-dominated eras, through a non-relativistic approximation. The results can be numerically reproduced and compared with those of considering non-relativistic and ultra-relativistic neutrinos into the computation of the matter power spectrum. Lessons concerning the formation of large scale structures of a DFG are depicted, and consequent deviations from standard $\Lambda$CDM predictions for the matter power spectrum (with and without neutrinos) are quantified.

The hypothesis of a tiny fraction of the cosmic inventory evolving cosmologically as a degenerate Fermi gas test fluid at some dominant cosmological background is investigated. Our analytical results allow for performing preliminary computations to the evolution of perturbations for relativistic and non-relativistic test fluids. The density fluctuation, $\delta$, the fluid velocity divergence, $\theta$, and an explicit expression for the dynamics of the shear stress, $\sigma$, are obtained for a degenerate Fermi gas in the background regime of radiation. Extensions to the dominance of matter and to the $\Lambda$CDM cosmological background are also investigated and lessons concerning the formation of large structures of degenerate Fermi gas are depicted.