The late integrated Sachs-Wolfe (ISW) effect correlates the Cosmic Microwave Background (CMB) temperature anisotropies with foreground cosmic large-scale structures. As the correlation depends crucially on the growth history in the era of dark energy, it is a key observational probe for constraining the cosmological model. Here we present a detailed study based on full-sky and deep light cones generated from very large volume numerical N-body simulations, which allow us to avoid the use of standard replica techniques, while capturing the entirety of the late ISW effect on the large scales. We post-process the light cones using an accurate ray-tracing method and construct full-sky maps of the ISW temperature anisotropy for three different dark energy models. We quantify in detail the extent to which the ISW effect can be used to discriminate between different dark energy scenarios when cross-correlated with the matter distribution or the CMB lensing potential. We also investigate the onset of non-linearities, the so-called Rees-Sciama effect which provides a complementary probe of the dark sector. We find the signal of the lensing-lensing and ISW-lensing correlation of the three dark energy models to be consistent with measurements from the Planck satellite. Future surveys of the large-scale structures may provide cross-correlation measurements that are sufficiently precise to distinguish the signal of these models. Our methodology is very general and can be applied to any dark energy or modified gravity scenario as long as the metric seen by photons can still be characterized by a Weyl potential.

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.

The upcoming generation of galaxy surveys will probe the distribution of matter in the universe with unprecedented accuracy. Measurements of the matter power spectrum at different scales and redshifts will provide stringent constraints on the cosmological parameters. However, on non-linear scales this will require an accurate evaluation of the covariance matrix. Here, we compute the covariance matrix of the 3D matter density power spectrum for the concordance $\Lambda$CDM cosmology from an ensemble of N-body simulations of the Dark Energy Universe Simulation - Parallel Universe Runs (DEUS-PUR). This consists of 12288 realisations of a $(656\,h^{-1}\,\textrm{Mpc})^3$ simulation box with $256^3$ particles. We combine this set with an auxiliary sample of 96 simulations of the same volume with $1024^3$ particles. We find N-body mass resolution effect to be an important source of systematic errors on the covariance at high redshift and small intermediate scales. We correct for this effect by introducing an empirical statistical method which provide an accurate determination of the covariance matrix over a wide range of scales including the Baryon Oscillations interval. Contrary to previous studies that used smaller N-body ensembles, we find the power spectrum distribution to significantly deviate from expectations of a Gaussian random density field at $k\gtrsim 0.25\,h\,\textrm{Mpc}^{-1}$ and $z<0.5$. This suggests that in the case of finite volume surveys an unbiased estimate of the ensemble averaged band power at these scales and redshifts may require a careful assessment of non-Gaussian errors more than previously considered.

Observations of colliding galaxy clusters with high relative velocity probe the tail of the halo pairwise velocity distribution with the potential of providing a powerful test of cosmology. As an example it has been argued that the discovery of the Bullet Cluster challenges standard $\Lambda$CDM model predictions. Halo catalogs from N-body simulations have been used to estimate the probability of Bullet-like clusters. However, due to simulation volume effects previous studies had to rely on a Gaussian extrapolation of the pairwise velocity distribution to high velocities. Here, we perform a detail analysis using the halo catalogs from the Dark Energy Universe Simulation Full Universe Runs (DEUS-FUR), which enables us to resolve the high-velocity tail of the distribution and study its dependence on the halo mass definition, redshift and cosmology. Building upon these results we estimate the probability of Bullet-like systems in the framework of Extreme Value Statistics. We show that the tail of extreme pairwise velocities significantly deviates from that of a Gaussian, moreover it carries an imprint of the underlying cosmology. We find the Bullet Cluster probability to be two orders of magnitude larger than previous estimates, thus easing the tension with the $\Lambda$CDM model. Finally, the comparison of the inferred probabilities for the different DEUS-FUR cosmologies suggests that observations of extreme interacting clusters can provide constraints on dark energy models complementary to standard cosmological tests.

We investigate the non-linear evolution of Baryon Acoustic Oscillations (BAO) in the low-redshift matter power spectrum from the DEUS-FUR $\Lambda$CDM model simulation. This is the first cosmological N-body simulation encompassing the full observable cosmic volume, thus allowing cosmic variance limited predictions at BAO scales. We control the effect of numerical systematic errors using a series of large volume high-resolution simulations. The combined analysis allows us to measure the matter power spectrum between $z=0$ and $1$ to 1% over the entire BAO range, $0.03<k [\textrm{h Mpc}^{-1}]<0.3$, in bins of size $\Delta k/k\lesssim$ 1%. We define the BAO with respect to a non-linearly evolved wiggle-free spectrum and determine the characteristics of the BAO without recurring to extrapolation from global fitting functions. We quantify the effects of non-linearities on the position and amplitude of the BAO extrema, and the coupling to the broadband slope of the power spectrum. We use these estimates to test non-linear predictions from semi-analytical models. Quite remarkably from the analysis of the redshift evolution of BAO we find that the second dip and third peak remains unaltered by non-linear effects. Furthermore, we find that the square of the damping factor and the shift of the position of BAO extrema scale to good approximation as the square of the growth factor, in agreement with expectations from perturbation theory. This confirms the idea that, besides cosmic distances, an accurate measurement of BAO at different redshifts can directly probe the growth of cosmic structures.

We study the imprint of Dark Energy on the density profile of Dark Matter halos using a set of high-resolution large volume cosmological N-body simulations from the Dark Energy Universe Simulation Series (DEUSS). We first focus on the analysis of the goodness-of-fit of the Navarro-Frenk-White (NFW) profile which we find to vary with halo mass and redshift. We also find that the fraction of halos ill-fitted by NFW varies with cosmology, thus indicating that the mass assembly of halos with perturbed density profiles carries a characteristic signature of Dark Energy. To access this information independently of any parametric profile, we introduce a new observable quantity: the halo sparsity $s_\Delta$. This is defined as the mass ratio $M_{200}/M_\Delta$, i.e. the ratio of mass inside a sphere of radius $r_{200}$ to that contained within a radius $r_\Delta$, enclosing 200 and $\Delta$ times the mean matter density respectively. We find the average sparsity to be nearly independent of the total halo mass, while its value can be inferred to better than a few percent from the ratio of the integrated halo mass functions at overdensities $\Delta$ and 200 respectively. This provides a consistency relation that can validate observational measurements of the halo sparsity. Most importantly, the sparsity significantly varies with the underlying Dark Energy model, thus providing an alternative cosmological probe.

We have performed the first-ever numerical N- body simulation of the full observable universe (DEUS "Dark Energy Universe Simulation" FUR "Full Universe Run"). This has evolved 550 billion particles on an Adaptive Mesh Refinement grid with more than two trillion computing points along the entire evolutionary history of the universe and across 6 order of magnitudes length scales, from the size of the Milky Way to that of the whole observable universe. To date, this is the largest and most advanced cosmological simulation ever run. It provides unique information on the formation and evolution of the largest structure in the universe and an exceptional support to future observational programs dedicated to mapping the distribution of matter and galaxies in the universe. The simulation has run on 4752 (of 5040) thin nodes of BULL supercomputer CURIE, using more than 300 TB of memory for 10 million hours of computing time. About 50 PBytes of data were generated throughout the run. Using an advanced and innovative reduction workflow the amount of useful stored data has been reduced to 500 TBytes.

The universality of the halo mass function is investigated in the context of dark energy cosmologies. This widely used approximation assumes that the mass function can be expressed as a function of the matter density omega_m and the rms linear density fluctuation sigma only, with no explicit dependence on the properties of dark energy or redshift. In order to test this hypothesis we run a series of 15 high-resolution N-body simulations for different cosmological models. These consists of three LCDM cosmologies best fitting WMAP-1, 3 and 5 years data, and three toy-models characterized by a Ratra-Peebles quintessence potential with different slopes and amounts of dark energy density. These toy models have very different evolutionary histories at the background and linear level, but share the same sigma8 value. For each of these models we measure the mass function from catalogues of halos identified in the simulations using the Friend-of-Friend (FoF) algorithm. We find redshift dependent deviations from a universal behaviour, well above numerical uncertainties and of non-stochastic origin, which are correlated with the linear growth factor of the investigated cosmologies. Using the spherical collapse as guidance, we show that such deviations are caused by the cosmology dependence of the non-linear collapse and virialization process. For practical applications, we provide a fitting formula of the mass function accurate to 5 percents over the all range of investigated cosmologies. We also derive an empirical relation between the FoF linking parameter and the virial overdensity which can account for most of the deviations from an exact universal behavior. Overall these results suggest that the halo mass function contains unique cosmological information since it carries a fossil record of the past cosmic evolution.

In this "Invisible Universe" proceedings, we introduce the Dark Energy Universe Simulation Series (DEUSS) which aim at investigating the imprints of realistic dark energy models on cosmic structure formation. It represents the largest dynamical dark energy simulation suite to date in term of spatial dynamics. We first present the 3 realistic dark energy models (calibrated on latest SNIa and CMB data): LambdaCDM, quintessence with Ratra-Peebles potential, and quintessence with Sugra potential. We then isolate various contributions for non-linear matter power spectra from a series of pre-DEUSS high-resolution simulations (130 million particles). Finally, we introduce DEUSS which consist in 9 Grand Challenge runs with 1 billion particles each thus probing scales from 4 Gpc down to 3 kpc at z=0. Our goal is to make these simulations available to the community through the "Dark Energy Universe Virtual Observatory" (DEUVO), and the "Dark Energy Universe Simulations" (DEUS) consortium.

We build a mechanism of gravitational symmetry breaking (GSB) of a global U(1) symmetry based on the relaxation of the equivalence principle due to the mass variation of pseudo Nambu-Goldstone dark matter (DM) particles. This GSB process is described by the modified cosmological convergence mechanism of the Abnormally Weighting Energy (AWE) Hypothesis previously introduced by the authors. Several remarkable constraints from the Hubble diagram of far-away supernovae are derived, notably on the explicit and gravitational symmetry breaking energy scales of the model. We then briefly present some consequences on neutrino masses when this mechanism is applied to the particular case of the breaking of lepton number symmetry.

We generalize tensor-scalar theories of gravitation by the introduction of an abnormally weighting type of energy. This theory of tensor-scalar anomalous gravity is based on a relaxation of the weak equivalence principle that is now restricted to ordinary visible matter only. As a consequence, the convergence mechanism toward general relativity is modified and produces naturally cosmic acceleration as an inescapable gravitational feedback induced by the mass-variation of some invisible sector. The cosmological implications of this new theoretical framework are studied. This glimpses at an enticing new symmetry between the visible and invisible sectors, namely that the scalar charges of visible and invisible matter are exactly opposite.

We present N-body cosmological numerical simulations including a primordial non-Gaussianity in the density fluctuation field quantified by the non-linear parameter $f_{NL}$. We have used MPGRAFIC code to produce initial conditions and the Adaptive Mesh Refinement (AMR) code RAMSES to evolve the large scale structure formation. We estimated the higher order momenta of the initial distribution of density fluctuations, investigated the redshift evolution of the non-linear power spectrum and estimated the discrepancy introduced by the primordial non-Gaussianity in the non-linear power spectrum.

Dark energy as a quintessence component causes a typical modification of the background cosmic expansion, which in addition to its clustering properties, can leave a potentially distinctive signature on large scale structures. Many previous studies have investigated this topic, particularly in relation to the non-linear regime of structure formation. However, no careful pre-selection of viable quintessence models with high precision cosmological data was performed. Here we show that this has led to a misinterpretation (and underestimation) of the imprint of quintessence on the distribution of large scale structures. To this purpose we perform a likelihood analysis of the combined Supernova Ia UNION dataset and WMAP5-years data to identify realistic quintessence models. Differences from the vanilla LambdaCDM are especially manifest in the predicted amplitude and shape of the linear matter power spectrum, though these remain within the uncertainties of the SDSS data. We use these models as benchmark for studying the clustering properties of dark matter halos by performing a series of high resolution N-body simulations. We find that realistic quintessence models allow for relevant differences of the dark matter distribution with the respect to the LambdaCDM scenario well into the non-linear regime, with deviations up to 40% in the non-linear power spectrum. Such differences are shown to depend on the nature of DE, as well as the scale and epoch considered. At small scales (k~1-5 h Mpc^{-1}, depending on the redshift) the structure formation process is about 20% more efficient than in LambdaCDM. We show that these imprints are a specific record of the cosmic structure formation history in DE cosmologies and therefore cannot be accounted in standard fitting functions of the non-linear matter power spectrum.

In order to quantitatively test the ability of averaged inhomogeneous cosmologies to correctly describe observations of the large scale properties of the Universe, we introduce a smoothed template metric corresponding to a constant spatial curvature model at any time, but with an evolving curvature parameter. This metric is used to compute quantities along an approximate effective lightcone of the averaged model of the Universe. As opposed to the standard Friedmann model, we parameterize this template metric by exact scaling properties of an averaged inhomogeneous cosmology, and we also motivate this form of the metric by results on a geometrical smoothing of inhomogeneous cosmological hypersurfaces. We test our hypothesis for the template metric against supernova data and the position of the CMB peaks, and infer the goodness--of--fit and parameter uncertainties. We find that averaged inhomogeneous models can reproduce the observations without requiring an additional Dark Energy component (though a volume acceleration is still needed), and that current data do not disfavour our main assumption on the effective lightcone structure. We also show that the experimental uncertainties on the angular diameter distance and the Hubble parameter from Baryon Acoustic Oscillations measurements -- forseen in future surveys like the proposed EUCLID satellite project -- are sufficiently small to distinguish between a FLRW template geometry and the template geometry with consistently evolving curvature. (Abridged)

We generalize tensor-scalar theories of gravitation by the introduction of an abnormally weighting type of energy. This theory of tensor-scalar anomalous gravity is based on a relaxation of the weak equivalence principle that is now restricted to ordinary visible matter only. As a consequence, the convergence mechanism toward general relativity is modified and produces naturally cosmic acceleration as an inescapable gravitational feedback induced by the mass-variation of some invisible sector. The cosmological implications of this new theoretical framework are studied. From the Hubble diagram cosmological test \textit{alone}, this theory provides an estimation of the amount of baryons and dark matter in the Universe that is consistent with the independent cosmological tests of Cosmic Microwave Background (CMB) and Big Bang Nucleosynthesis (BBN). Cosmic coincidence is naturally achieved from a equally natural assumption on the amplitude of the scalar coupling strength. Finally, from the adequacy to supernovae data, we derive a new intriguing relation between the space-time dependences of the gravitational coupling and the dark matter mass, providing an example of crucial constraint on microphysics from cosmology. This glimpses at an enticing new symmetry between the visible and invisible sectors, namely that the scalar charges of visible and invisible matter are exactly opposite.

An S-brane solution with two non-composite electric branes and a set of l scalar fields is considered. The intersection rule for branes corresponds to the Lie algebra A_2. The solution contains five factor spaces with the fifth one interpreted as ``our'' 3-dimensional space. It is shown that there exists a time interval where accelerating expansion of ``our'' 3-dimensional space is compatible with small enough value of effective gravitational ``constant'' variation.

The Abnormally Weighting Energy (AWE) hypothesis consists of assuming that the dark sector of cosmology violates the weak equivalence principle (WEP) on cosmological scales, which implies a violation of the strong equivalence principle for ordinary matter. In this paper, dark energy (DE) is shown to result from the violation of WEP by pressureless (dark) matter. This allows us to build a new cosmological framework in which general relativity (GR) is satisfied at low scales, as WEP violation depends on the ratio of the ordinary matter over dark matter densities, but at large scales, we obtain a general relativity-like theory with a different value of the gravitational coupling. This explanation is formulated in terms of a tensor-scalar theory of gravitation without WEP for which there exists a revisited convergence mechanism toward GR. The consequent DE mechanism build upon the anomalous gravity of dark matter (i) does not require any violation of the strong energy condition $p<-\rho c^2/3$, (ii) offers a natural way-out of the coincidence problem thanks to the non-minimal couplings to gravitation, (iii) accounts fairly for supernovae data from various very simple couplings and with density parameters very close to the ones of the concordance model $\Lambda CDM$, therefore suggesting an explanation to its remarkable adequacy. Finally, (iv) this mechanism ends up in the future with an Einstein-de Sitter expansion regime once the attractor is reached.

We introduce the Abnormally Weighting Energy (AWE) hypothesis in which dark energy (DE) is presented as a consequence of the violation of the weak equivalence principle (WEP) at cosmological scales by some dark sector. Indeed, this implies a violation of the strong equivalence principle (SEP) for ordinary matter and consequent cosmic acceleration in the observable frame as well as variation of the gravitational constant. The consequent DE mechanism build upon the AWE hypothesis (i) does not require a violation of the strong energy condition $p<-\rho c^2/3$, (ii) assumes rather non-negligible direct couplings to the gravitational scalar field (iii) offers a natural convergence mechanism toward general relativity (iv) accounts fairly for supernovae data from various couplings and equations of state of the dark sector as well as density parameters very close to the ones of the concordance model $\Lambda CDM$. Finally (v), this AWE mechanism typically ends up with an Einstein-de Sitter expansion regime once the attractor is reached.

In the context of averaged cosmologies, the effective equations can be written in the form of "regional" Friedmannian equations with additional sources arising from the so-called backreaction of inhomogeneities. We propose a mean field description of this backreaction in terms of a regionally homogeneous scalar field: this provides a physical motivation to the phenomenological scalar fields generically called quintessence fields. We explicitly reconstruct the potential of the scalar field for a one-parameter family of scaling solutions to the backreaction problem, showing that it entails most of the standard scalar fields including e.g. standard and phantom quintessence scenarii.

Combined with other CMB experiments, the WMAP survey provides an accurate estimate of the baryon density of the Universe. In the framework of the standard Big Bang Nucleosynthesis (BBN), such a baryon density leads to predictions for the primordial abundances of $^{4}$He and D in good agreement with observations. However, it also leads to a significant discrepancy between the predicted and observed primordial abundance of $^{7}$Li. Such a discrepancy is often termed as 'the lithium problem'. In this paper, we analyze this problem in the framework of scalar-tensor theories of gravity. It is shown that an expansion of the Universe slightly slower than in General Relativity before BBN, but faster during BBN, solves the lithium problem and leads to $^4$He and D primordial abundances consistent with the observational constraints. This kind of behavior is obtained in numerous scalar-tensor models, both with and without a self-interaction potential for the scalar field. In models with a self-interacting scalar field, the convergence towards General Relativity is ensured without any condition, thanks to an attraction mechanism which starts to work during the radiation-dominated epoch.

Spatially averaged inhomogeneous cosmologies in classical general relativity can be written in the form of effective Friedmann equations with sources that include backreaction terms. In this paper we propose to describe these backreaction terms with the help of a homogeneous scalar field evolving in a potential; we call it the `morphon field'. This new field links classical inhomogeneous cosmologies to scalar field cosmologies, allowing to reinterpret, e.g., quintessence scenarios by routing the physical origin of the scalar field source to inhomogeneities in the Universe. We investigate a one-parameter family of scaling solutions to the backreaction problem. Subcases of these solutions (all without an assumed cosmological constant) include scale-dependent models with Friedmannian kinematics that can mimic the presence of a cosmological constant or a time-dependent cosmological term. We explicitly reconstruct the scalar field potential for the scaling solutions, and discuss those cases that provide a solution to the Dark Energy and coincidence problems. In this approach, Dark Energy emerges from morphon fields, a mechanism that can be understood through the proposed correspondence: the averaged cosmology is characterized by a weak decay (quintessence) or growth (phantom quintessence) of kinematical fluctuations, fed by `curvature energy' that is stored in the averaged 3-Ricci curvature. We find that the late-time trajectories of those models approach attractors that lie in the future of a state that is predicted by observational constraints.

Regionally averaged relativistic cosmologies have recently been considered as a possible explanation for the apparent late time acceleration of the Universe. This contribution reports on a mean field description of the backreaction in terms of a minimally coupled regionally homogeneous scalar field evolving in a potential, then giving a physical origin to the various phenomenological scalar fields generically called quintessence fields. As an example, the correspondence is then applied to scaling solutions.

We present a new interpretation of dark energy in terms of an \textit{Abnormally Weighting Energy} (AWE). This means that dark energy does not couple to gravitation in the same way as ordinary matter, yielding a violation of the weak and strong equivalence principles on cosmological scales. The resulting cosmological mechanism accounts for the Hubble diagram of type Ia supernovae in terms of both cosmic acceleration and variation of the gravitational constant while still accounting for the present tests of general relativity. This explanation allows to build dark energy models (i) without violation of the strong energy condition $p<-\rho c^2/3$ (ii) with non-negligible direct couplings to gravitation and (iii) natural convergence mechanism toward general relativity.

We investigate the possibility that dark energy does not couple to gravitation in the same way than ordinary matter, yielding a violation of the weak and strong equivalence principles on cosmological scales. We build a transient mechanism in which gravitation is pushed away from general relativity by a Born-Infeld gauge interaction acting as an "Abnormally Weighting" (dark) Energy. This mechanism accounts for the Hubble diagram of far-away supernovae by cosmic acceleration and time variation of the gravitational constant while accounting naturally for the present tests on general relativity.

The non-abelian Einstein-Born-Infeld-Dilaton theory, which rules the dynamics of tensor-scalar gravitation coupled to a $su(2)$-valued gauge field ruled by Born-Infeld lagrangian, is studied in a cosmological framework. The microscopic energy exchange between the gauge field and the dilaton which results from a non-universality of the coupling to gravity modifies the usual behaviour of tensor-scalar theories coupled to matter fluids. General cosmological evolutions are derived for different couplings to gravitation and a comparison to universal coupling is highlighted. Evidences of cosmic acceleration are presented when the evolution is interpreted in the Jordan physical frame of a matter respecting the weak equivalence principle. The importance for the mechanism of cosmic acceleration of the dynamics of the Born-Infeld gauge field, the attraction role of the matter fluid and the non-universality of the gravitational couplings is briefly outlined.

The problem of galaxy formation and its dependence on thermodynamic properties is addressed by using Eulerian hydrodynamic numerical simulations of large scale structure formation. Global galaxy properties are explored in simulations including gravitation, shock heating and cooling processes, and following self-consistently the chemical evolution of a primordial composition hydrogen-helium plasma without assuming collisional ionization equilibrium. The galaxy formation model is mainly based on the identification of converging dense cold gas regions. We show that the evolution at low redshift of the observed cosmic star formation rate density is reproduced, and that the galaxy-like object mass function is dominated by low-mass objects. The galaxy mass functions are well described by a two power-law Schechter function whose parameters are in good agreement with observational fits of the galaxy luminosity function. The high-mass end of the galaxy mass function includes objects formed at early epochs and residing in high-mass dark matter halos whereas the low-mass end includes galaxies formed at later epochs and active in their ``stellar'' mass formation. Finally, the influence of two other physical processes, photoionization and non-equipartition processes between electrons, ions and neutrals of the cosmological plasma is discussed and the modifications on galaxy formation are examined.

Using N-body/hydrodynamic simulations, the influence of non-equipartition processes on the thermal and dynamical properties of cosmological baryonic gas is investigated. We focus on a possible departure from equilibrium between electrons, ions and neutral atoms in low temperature (10^4-10^6 K) and weakly ionized regions of the intergalactic medium. The simulations compute the energy exchanges between ions, neutrals and electrons, without assuming thermal equilibrium. They include gravitation, shock heating and cooling processes, and follow self-consistently the chemical evolution of a primordial composition hydrogen-helium plasma without assuming collisional ionization equilibrium. At high redshift, a significant fraction of the intergalactic medium is found to be warmer and weakly ionized in simulations with non-equipartition processes than in simulations in which the cosmological plasma is considered to be in thermodynamic equilibrium. With a semi-analytical study of the out of equilibrium regions we show that, during the formation of cosmic structures, departure from equilibrium in accreted plasma results from the competition between the atomic cooling processes and the elastic processes between heavy particles and electrons. Our numerical results are in agreement with this semi-analytical model. Therefore, since baryonic matter with temperatures around 10^4 K is a reservoir for galaxy formation, non-equipartition processes are expected to modify the properties of the objects formed.

In this paper, we analyze the conditions for convergence toward General Relativity of scalar-tensor gravity theories defined by an arbitrary coupling function $\alpha$ (in the Einstein frame). We show that, in general, the evolution of the scalar field $(\phi)$ is governed by two opposite mechanisms: an attraction mechanism which tends to drive scalar-tensor models toward Einstein's theory, and a repulsion mechanism which has the contrary effect. The attraction mechanism dominates the recent epochs of the universe evolution if, and only if, the scalar field and its derivative satisfy certain boundary conditions. Since these conditions for convergence toward general relativity depend on the particular scalar-tensor theory used to describe the universe evolution, the nucleosynthesis bounds on the present value of the coupling function, $\alpha_0$, strongly differ from some theories to others. For example, in theories defined by $\alpha \propto \mid\phi\mid$ analytical estimates lead to very stringent nucleosynthesis bounds on $\alpha_0$ ($\lesssim 10^{-19}$). By contrast, in scalar-tensor theories defined by $\alpha \propto \phi$ much larger limits on $\alpha_0$ ($\lesssim 10^{-7}$) are found.

We present a method which enables exact solutions to be found for at homogeneous and isotropic scalar-tensor cosmologies with an arbitrary $\omega(\Phi)$ function and satisfying the general perfect fluid state equation $P=(\gamma-1)\rho c^2$. This method has been used to analyze a wide range of asymptotic analytical solutions at early and late times for different epochs in the cosmic history: false vacuum inflationary models, vacuum and radiation-dominated models, and matter-dominated models. We also describe the qualitative behavior of models at intermediary times and give exact solutions at any time for some particular scalar-tensor theories.

We study the stability of rotating collisionless self-gravitating spherical systems by using high resolution N-body experiments on a Connection Machine CM-5. We added rotation to Ossipkov-Merritt (hereafter OM) anisotropic spherical systems by using two methods. The first method conserves the anisotropy of the distribution function defined in the OM algorithm. The second method distorts the systems in velocity space. We then show that the stability of systems depends both on their anisotropy and on the value of the ratio between the total kinetic energy and the rotational kinetic energy. We also test the relevance of the stability parameters introduced by Perez et al. (1996) for the case of rotating systems.

We investigate the influence of non equilibrium thermodynamics on cosmological structure formation. In this paper, we consider the collapse of planar perturbations usually called "Zel'dovich pancakes". We have developed for that purpose a new two fluids (gas and dark matter) hydrodynamical code, with three different thermodynamical species: electrons, ions and neutral particles (T_e\ne T_i \ne T_n). We describe in details the complex structure of accretion shock waves. We include several relevant processes for a low density, high temperature, collisional plasma such as non-equilibrium chemical reactions, cooling, shock heating, thermal energy equipartition between electrons, ions and neutral particles and electronic conduction. We find two different regions in the pancake structure: a thermal precursor ahead of the compression front and an equipartition wave after the compression front where electrons and ions temperatures differ significantly. This complex structure may have two interesting consequences: pre-heating of unshocked regions in the vicinity of massive X-ray clusters and ions and electrons temperatures differences in the outer regions of X-rays clusters.

We compare the mass functions obtained analytically, in the framework of an extended Press & Schechter (PS) formalism, in a previous paper to the numerical mass functions obtained in N -- body simulations, using different algorithms to define objects in the density field. After discussing the properties of the algorithms, we show that the mass function obtained using the friend -- of -- friend algorithm reproduces best the scaling behaviors predicted in the extended PS formalism. Following this statistical analysis, we show that it is possible in the framework of our extended PS formalism to reproduce the mass function but also, and for the first time, the initial statistical properties of structures and their collapse time. This allow to present a ``coherent'' picture of structure formation which can account for the initial, final and dynamical properties of structures.

We investigate the role of radiative cooling within the core of large X-ray clusters using multi-phase hydrodynamics. We developed for that purpose a spherically symmetric hydrodynamical code, coupled to a "fluid model" that describes accurately the dark matter component. Cooling is included using a self-consistent multi-phase approach, leading to "cooled gas" mass deposition throughout the flow. We simulate the collapse and the subsequent evolution of a Coma-like X-ray cluster, avoiding the well-known "cooling catastrophe". The total mass profile of our simulated cluster is very similar to the "universal" profile proposed by Navarro, Frenk & White (1995). More interestingly, we also obtain a quasi-isothermal temperature profile, which is a direct consequence of multi-phase cooling within such a potential well.

We simulate the formation of a large X-ray cluster using a fully 3D hydrodynamical code coupled to a Particle-Mesh scheme which models the dark matter component. We focus on a possible decoupling between electrons and ions temperatures. We then solve the energy transfer equations between electrons, ions and neutrals without assuming thermal equilibrium between the three gases (T_e <> T_i <> T_n). We solve self-consistently the chemical equations for an hydrogen/helium primordial plasma without assuming ionization-recombination equilibrium. We find that the electron temperature differs from the true dynamical temperature by 20% at the Virial radius of our simulated cluster. This could lead marginally to an underestimate of the total mass in the outer regions of large X-ray clusters.

In a previous paper (Serna & Alimi 1996), we have pointed out the existence of some particular scalar-tensor gravity theories able to relax the nucleosynthesis constraint on the cosmic baryonic density. In this paper, we present an exhaustive study of primordial nucleosynthesis in the framework of such theories taking into account the currently adopted observational constraints. We show that a wide class of them allows for a baryonic density very close to that needed for the universe closure. This class of theories converges soon enough towards General Relativity and, hence, is compatible with all solar-system and binary pulsar gravitational tests. In other words, we show that primordial nucleosynthesis does not always impose a very stringent bound on the baryon contribution to the density parameter.

We reconsider the collapse of cosmic structures in an Einstein-de Sitter Universe, using the self similar initial conditions of Fillmore & Goldreich (1984). We first derive a new approximation to describe the dark matter dynamics in spherical geometry, that we refer to the "fluid approach". This method enables us to recover the self-similarity solutions of Fillmore & Goldreich for dark matter. We derive also new self-similarity solutions for the gas. We thus compare directly gas and dark matter dynamics, focusing on the differences due to their different dimensionalities in velocity space. This work may have interesting consequences for gas and dark matter distributions in large galaxy clusters, allowing to explain why the total mass profile is always steeper than the X-ray gas profile. We discuss also the shape of the dark matter density profile found in N-body simulations in terms of a change of dimensionality in the dark matter velocity space. The stable clustering hypothesis has been finally considered in the light of this analytical approach.

We investigate the formation of spherical cosmological structures following both dark matter and gas components. We focus on the dynamical aspect of the collapse assuming an adiabatic, $\gamma = 5/3$, fully ionized primordial plasma. We use for that purpose a fully Lagrangian hydrodynamical code designed to describe highly compressible flows in spherical geometry. We investigate also a "fluid approach" to describe the mean physical quantities of the dark matter flow. We test its validity for a wide range of initial density contrast. We show that an homogeneous isentropic core forms in the gas distribution, surrounded by a self-similar hydrostatic halo, with much higher entropy generated by shock dissipation. We derive analytical expressions for the size, density and temperature of the core, as well as for the surrounding halo. We show that, unless very efficient heating processes occur in the intergalactic medium, we are unable to reproduce within adiabatic models the typical core sizes in X-ray clusters. We also show that, for dynamical reasons only, the gas distribution is naturally antibiased relative to the total mass distribution, without invoking any reheating processes. This could explain why the gas fraction increases with radius in very large X-ray clusters. As a preparation for the next study devoted to the thermodynamical aspect of the collapse, we investigate the initial entropy level required to solve the core problem in X-ray clusters.

We investigate some modifications to the Press & Schechter (1974) (PS) prescription resulting from shear and tidal effects. These modifications rely on more realistic treatments of the collapse process than the standard approach based on the spherical model. First, we show that the mass function resulting from a new approximate Lagrangian dynamic (Audit & Alimi 96), contains more objects at high mass, than the classical PS mass function and is well fitted by a PS-like function with a threshold density of $\delta_c \simeq 1.4$. However, such a Lagrangian description can underestimate the epoch of structure formation since it defines it as the collapse of the first principal axis. We therefore suggest some analytical prescriptions, for computing the collapse time along the second and third principal axes, and we deduce the corresponding mass functions. The collapse along the third axis is delayed by the shear and the number of objects of high mass then decreases. Finally, we show that the shear also strongly affects the formation of low-mass halos. This dynamical effect implies a modification of the low-mass slope of the mass function and allows the reproduction of the observed luminosity function of field galaxies. In a companion paper, we present results of numerical simulations which complete this work.

In this paper we present a new local Lagrangian approximation to the gravitational dynamics of cold matter. We describe the dynamics of a Lagrangian fluid element through only one quantity, the deformation tensor. We show that this tensor is clearly suited to the study of gravitational dynamics and, moreover, that knowing its evolution is enough to completely describe a fluid element. Our approximate dynamics treat exactly the conservation of mass, of the velocity divergence and of the shear and is exact in the case of planar, cylindrical and spherical collapses. It also reproduces very accurately the evolution of all dynamical quantities for a very wide class of initial conditions as illustrated by a detailed comparison with the homogeneous ellipsoid model. Beside, we highlight, for the first time, the important dynamical role played by the Newtonian counterpart of the magnetic part of the Weyl tensor.

We present and test a general-purpose code, called PPASPH, for evolving self-gravitating fluids in astrophysics, both with and without a collisionless component. In PPASPH, hydrodynamical properties are computed by using the SPH (Smoothed Particle Hydrodynamics) method while, unlike most previous implementations of SPH, gravitational forces are computed by a PP (Particle-Particle) approach. another important feature of this code is that hydrodynamics equations optionally include the correction terms appearing when $h(t,\br)$ is not constant. Our code has been implemented by using the data parallel programming model on CM5(CM). PPASPH has been applied to study the importance of adaptative smoothing correction terms on the entropy conservation. We confirm Hernquist's (1993) interpretation of the entropy violation observed in previous SPH simulations as a result of having neglected these terms. An improvement on the entropy conservation is not found by merely considering larger numbers of particles or different $N_S$ choices. The correct continuum description is only obtained if the $\bn h$ correction terms are included. Otherwise, the entropy conservation is always rather poor as compared to that found for the total energy.

We have performed a series of high resolution N-body experiments on a Connection Machine CM-5 in order to study the stability of collisionless self-gravitating spherical systems. We interpret our results in the framework of symplectic mechanics, which provides the definition of a new class of particular perturbations: The preserving perturbations, which are a generalization of the radial ones. Using models defined by the Ossipkov-Merritt algorithm, we show that the stability of a spherical anisotropic system is directly related to the preserving or non-preserving nature of the perturbations acting on the system. We then generalize our results to all spherical systems. Since the ``isotropic component'' of the linear variation of the distribution function cannot be used to predict the stability or instability of a spherical system, we propose a more useful stability parameter which is derived from the ``anisotropic'' component of the linear variation.

We analyze the qualitative behaviors of scalar-tensor cosmologies with an arbitrary monotonic $\omega(\Phi)$ function. In particular, we are interested on scalar-tensor theories distinguishable at early epochs from General Relativity (GR) but leading to predictions compatible with solar-system experiments. After extending the method developed by Lorentz-Petzold and Barrow, we establish the conditions required for convergence towards GR at $t\rightarrow\infty$. Then, we obtain all the asymptotic analytical solutions at early times which are possible in the framework of these theories. The subsequent qualitative evolution, from these asymptotic solutions until their later convergence towards GR, has been then analyzed by means of numerical computations. From this analysis, we have been able to establish a classification of the different qualitative behaviors of scalar-tensor cosmological models with an arbitrary monotonic $\omega(\Phi)$ function.

We present a detailed calculation of the light element production in the framework of Scalar-Tensor (ST) theories of gravitation. The coupling function \omega has been described by an appropriate form which reproduces all the possible asymptotic behaviors at early times of viable ST cosmological models with a monotonic \omega(\Phi). This form gives an exact representation for most of the particular theories proposed in the literature, but also a first-order approximation to many other theories. In most of ST theories, the nucleosynthesis bounds lead to cosmological models which do not significantly differ from the standard FRW ones. We have found however a particular class of ST theories where the expansion rate of the universe during nucleosynthesis can be very different from that found in GR, while the present value of \omega is high enough to ensure compatibility with solar-system experiments. In the framework of this class of theories, right primordial abundances can be obtained for a baryon density range much wider (2.8 \lesssim \eta_{10} \lesssim 58.7) than in GR. Consequently, the usual constraint on the baryon contribution to the density parameter of the universe can be drastically relaxed (0.01 \lesssim \Omega_{b0} \lesssim 1.38) by considering these gravity theories. This is the first time that a ST theory is found to be compatible both with primordial nucleosynthesis and solar-system experiments while implying cosmological models significantly different from the FRW ones.

We propose a morphological multi-scale analysis of large-scale structures obtained by computer simulations and by observations. Structures are obtained at different scales by applying a wavelet transform on the observed and simulated data. Application of a segmentation algorithm allows a quantitative morphological description of the structures at each scale. The morphological parameter which we propose represents the deviation of a structure from sphericity. The dependence on scale of this simple parameter is shown to characterize in simulations the underlying cosmological model. We compare the HDM, CDM and mixed models with the CfA catalogue. This comparison favours a mixed model containing $65\%$ of CDM, $30\%$ of HDM, and $5\%$ baryonic matter.

We present the first results of a galaxy redshift survey, ESO Slice Project (ESP), we are accomplishing as an ESO Key-Project over about 30 square degrees in a region near the South Galactic Pole. The limiting magnitude is b_J = 19.4. Observations have been almost completed and about 90% of the data obtained so far has been reduced providing about 3000 galaxy redshifts. We present some preliminary results concerning the large scale galaxy distribution and their luminosity function.

We present the first results of a galaxy redshift survey, ESO Slice Project (ESP), we are accomplishing as an ESO Key--Project over about 40 square degrees in a region near the South Galactic Pole. The limiting magnitude is $b_J = 19.4$. Up to now $\sim 85\%$ of the observations has been completed and $\sim 65\%$ of the data has been reduced providing $\sim 2000$ galaxy redshifts. We present some preliminary results concerning the large scale galaxy distribution, the luminosity function and the properties of galaxies with spectra showing prominent emission lines.

We have performed a series of N-body experiments on Connection Machine-5 in order to simulate the formation of galaxy clusters gravitationally dominated by a massive dark background. In accordance with previous authors we find an extremely inhomogeneous evolution where subcondensations are continually formed and merged. The final distribution of galaxies is more centrally condensed than that of dark matter particles. We have analyzed the origin of this galaxy-dark matter segregation and also the origin of subclustering leading to this effect. We have then analysed, analytically and numerically, the dependence of final segregation on the physical parameters characterizing the model of protocluster. We also find that such a segregation persists even when inelastic encounters of dark halos around galaxies are taken into account. We conclude that this effect cannot be in general avoided in any hierarchical clustering scenario.