The absence of a physically motivated model for large scale profiles of cosmic voids limits our ability to extract valuable cosmological information from their study. In this paper, we address this problem by introducing the spherically compensated cosmic regions, named CoSpheres. Such cosmic regions are identified around local extrema in the density field and admit a unique compensation radius $R_1$ where the internal spherical mass is exactly compensated. Their origin is studied by extending the standard peak model and implementing the compensation condition. Since the compensation radius evolves as the Universe itself, $R_1(t)\propto a(t)$, CoSpheres behave as bubble Universes with fixed comoving volume. Using the spherical collapse model, we reconstruct their profiles with a very high accuracy until $z=0$ in N-body simulations. CoSpheres are symmetrically defined and reconstructed for both central maximum (seeding haloes and galaxies) and minimum (identified with cosmic voids). We show that the full non linear dynamics can be solved analytically around this particular compensation radius, providing useful predictions for cosmology. This formalism highlights original correlations between local extremum and their large scale cosmic environment. The statistical properties of these spherically compensated cosmic regions and the possibilities to constrain efficiently both cosmology and gravity will be investigated in companion papers.

The statistical properties of cosmic structures are well known to be strong probes for cosmology. In particular, several studies tried to use the cosmic void counting number to obtain tight constrains on Dark Energy. In this paper we address this question by using the CoSphere model as introduced in de Fromont & Alimi (2017a). We derive their exact statistics in both primordial and non linearly evolved Universe for the standard $\Lambda$CDM model. We first compute the full joint Gaussian probability distribution for the various parameters describing these profiles in the Gaussian Random Field. We recover the results of Bardeen et al. (1986) only in the limit where the compensation radius becomes very large, i.e. when the central extremum decouples from its cosmic environment. We derive the probability distribution of the compensation size in this primordial field. We show that this distribution is redshift independent and can be used to model cosmic void size distribution. Interestingly, it can be used for central maximum such as DM haloes. We compute analytically the statistical distribution of the compensation density in both primordial and evolved Universe. We also derive the statistical distribution of the peak parameters already introduced by Bardeen et al. (1986) and discuss their correlation with the cosmic environment. We thus show that small central extrema with low density are associated with narrow compensation regions with a small $R_1$ and a deep compensation density $\delta_1$ while higher central extrema are located in larger but smoother over/under massive regions.

We re-explore the Bi- and Multi-Galileon models with trivial asymptotic conditions at infinity and show that propagation of superluminal fluctuations is a common and unavoidable feature of these theories, unlike previously claimed in the literature. We show that all Multi-Galileon theories containing a Cubic Galileon term exhibit superluminalities at large distances from a point source, and that even if the Cubic Galileon is not present one can always find sensible matter distributions in which there are superluminal modes at large distances. In the Bi-Galileon case we explicitly show that there are always superluminal modes around a point source even if the Cubic Galileon is not present. Finally, we briefly comment on the possibility of avoiding superluminalities by modifying the asymptotic conditions at infinity.