Large primordial Black Hole (PBH) formation is enhanced if strongly coupled scalar and spinor fields ($\Phi$ and $\psi$) are a stable cosmic component since the primeval radiative expansion (SCDEW models). In particular, we show that PBH formation is easier at a specific time, i.e., when the asymptotic mass $m_H$, acquired by the $\psi$ field at the higgs scale, becomes dominant, so that the typical BH mass $M_{BH}$ depends on $m_H$ value. For instance, if $m_H \sim 100\,$ eV $(1$ keV$)$ and the coupling $\beta \sim 8.35 (37)$, PBH with $M_{BH} \simeq 10^7-10^8 $ $(10^3-10^4)\, M_\odot$ could form. The very mechanism enhancing PBH formation also causes technical difficulties to evaluate the transfer function of SCDEW models at high $k$. A tentative solution of this problem leaves only minor discrepancies from $\Lambda$CDM, also at these scales, gradually vanishing for greater $m_H$ values. We conclude that, for suitable parameter choices, SCDEW models could be the real physics underlying $\Lambda$ CDM, so overcoming its fine tuning and coincidence problems, with the extra bonus of yielding large BH seeds.

Strongly Coupled Dark Energy plus Warm dark matter (SCDEW) cosmologies are based on the finding of a conformally invariant (CI) attractor solution during the early radiative expansion, requiring then the stationary presence of $\sim 1\, \%$ of coupled-DM and DE, since inflationary reheating. In these models, coupled-DM fluctuations, even in the early radiative expansion, grow up to non-linearity, as shown in a previous associated paper. Such early non-linear stages are modelized here through the evolution of a top-hat density enhancement. As expected, its radius $R$ increases up to a maximum and then starts to decrease. Virial balance is reached when the coupled-DM density contrast is just 25-26 and DM density enhancement is $\cal O$$(10\, \%)$ of total density. Moreover, we find that this is not an equilibrium configuration as, afterwards, coupling causes DM particle velocities to increase, so that the fluctuation gradually dissolves. We estimate the duration of the whole process, from horizon crossing to dissolution, and find $z_{horizon}/z_{erasing} \sim 3 \times 10^4$. Therefore, only fluctuations entering the horizon at $z \lesssim 10^9$-$10^{10}$ are able to accrete WDM with mass $\sim 100\, $eV -as soon as it becomes non-relativistic- so avoiding full disruption. Accordingly, SCDEW cosmologies, whose WDM has mass $\sim 100\, $eV, can preserve primeval fluctuations down to stellar mass scale.

Strongly Coupled Dark Energy plus Warm dark matter (SCDEW) cosmologies admit the stationary presence of $\sim 1\, \%$ of coupled-DM and DE, since inflationary reheating. Coupled-DM fluctuations therefore grow up to non-linearity even in the early radiative expansion. Such early non-linear stages are modelized here through the evolution of a top-hat density enhancement, reaching an early virial balance when the coupled-DM density contrast is just 25-26 and DM density enhancement is $ \sim 10\, \%$ of total density. During the time needed to settle in virial equilibium, the virial balance conditions however continue to modify, so that "virialized" lumps undergo a complete evaporation. Here we outline that DM particles processed by overdentities preserve a fraction of their virial momentum. Although fully non-relativistic, the resulting velocities (moderately) affect the fluctuation dynamics over greater scales, entering the horizon later on.

Cosmologies including strongly Coupled (SC) Dark Energy (DE) and Warm dark matter (SCDEW) are based on a conformally invariant (CI) attractor solution modifying the early radiative expansion. Then, aside of radiation, a kinetic field $\Phi$ and a DM component account for a stationary fraction, $\sim 1\, \%$, of the total energy. Most SCDEW predictions are hardly distinguishable from LCDM, while SCDEW alleviates quite a few LCDM conceptual problems, as well as its difficulties to meet data below the average galaxy scale. The CI expansion begins at the inflation end, when $\Phi$ (future DE) possibly plays a role in reheating, and ends at the Higgs' scale. Afterwards, a number of viable options is open, allowing for the transition from the CI expansion to the present Universe. In this paper: (i) We show how the attractor is recovered when the spin degrees of freedom decreases. (ii) We perform a detailed comparison of CMB anisotropy and polarization spectra for SCDEW and LCDM, including tensor components, finding negligible discrepancies. (iii) Linear spectra exhibit a greater parameter dependence at large $k$'s, but are still consistent with data for suitable parameter choices. (iv) We also compare previous simulation results with fresh data on galaxy concentration. Finally, (v) we outline numerical difficulties at high $k$. This motivates a second related paper, where such problems are treated in a quantitative way.

In the early Universe, strongly interacting matter was a quark-gluon plasma. Both lattice computations and heavy ion collision experiments however tell us that, in the absence of chemical potentials, no plasma survives at $T <\sim 150$ MeV. The cosmological Quark-Hadron transition, however, seems to have been a crossover; cosmological consequences envisaged when it was believed to be a phase transition no longer hold. In this paper we discuss whether even a crossover transition can leave an imprint that cosmological observations can seek or, viceversa, there are questions cosmology should address to QCD specialists. In particular, we argue that it is still unclear how baryons (not hadrons) could form at the cosmological transition. A critical role should be played by diquark states, whose abundance in the early plasma needs to be accurately evaluated. We estimate that, if the number of quarks belonging to a diquark state, at the beginning of the cosmological transition, is $<\sim 1:10^6$, its dynamics could be modified by the process of B-transfer from plasma to hadrons. In turn, by assuming B-transfer to cause just mild perturbations and, in particular, no entropy input, we study the deviations from the tracking regime, in the frame of SCDEW models. We find that, in some cases, residual deviations could propagate down to primeval nuclesynthesis.

We discuss an extension of the Coyote emulator to predict non-linear matter power spectra of dark energy (DE) models with a scale factor dependent equation of state of the form w = w_0 + ( 1 - a )w_a . The extension is based on the mapping rule between non-linear spectra of DE models with constant equation of state and those with time varying one originally introduced in ref. [40]. Using a series of N-body simulations we show that the spectral equivalence is accurate to sub-percent level across the same range of modes and redshift covered by the Coyote suite. Thus, the extended emulator provides a very efficient and accurate tool to predict non-linear power spectra for DE models with w_0 - w_a parametrization. According to the same criteria we have developed a numerical code, and we have implemented in a dedicated module for the CAMB code, that can be used in combination with the Coyote Emulator in likelihood analyses of non-linear matter power spectrum measurements. All codes can be found at https://github.com/luciano-casarini/PKequal

We study the effects of sample variance in N--body simulations, as a function of the size of the simulation box, namely in connection with predictions on tomographic shear spectra. We make use of a set of 8 $\Lambda$CDM simulations in boxes of 128, 256, 512 $h^{-1}$Mpc aside, for a total of 24, differing just by the initial seeds. Among the simulations with 128 and 512 $h^{-1}$Mpc aside, we suitably select those closest and farthest from {\it average}. Numerical and linear spectra $P(k,z)$ are suitably connected at low $k$ so to evaluate the effects of sample variance on shear spectra $C_{ij}(\ell)$ for 5 or 10 tomographic bands. We find that shear spectra obtained by using 128 $h^{-1}$Mpc simulations can vary up to $\sim 25\, \%$, just because of the seed. Sample variance lowers to $\sim 3.3\, \%$, when using 512 $h^{-1}$Mpc. These very percentages could however slightly vary, if other sets of the same number of realizations were considered. Accordingly, in order to match the $\sim 1\, \%$ precision expected for data, if still using 8 boxes, we require a size $\sim 1300$ --$ 1700 \, h^{-1}$ Mpc for them.

In this first paper we discuss the linear theory and the background evolution of a new class of models we dub SCDEW: Strongly Coupled DE, plus WDM. In these models, WDM dominates today's matter density; like baryons, WDM is uncoupled. Dark Energy is a scalar field $\Phi$; its coupling to ancillary CDM, whose today's density is $\ll 1\, \%$, is an essential model feature. Such coupling, in fact, allows the formation of cosmic structures, in spite of very low WDM particle masses ($\sim 100$ eV). SCDEW models yields Cosmic Microwave Background and linear Large Scale features substantially undistinguishable from $\Lambda$CDM, but thanks to the very low WDM masses they strongly alleviate $\Lambda$CDM issues on small scales, as confirmed via numerical simulations in the II associated paper. Moreover SCDEW cosmologies significantly ease the coincidence and fine tuning problems of $\Lambda$CDM and, by using a field theory approach, we also outline possible links with inflationary models. We also discuss a possible fading of the coupling at low redshifts which prevents non linearities on the CDM component to cause computational problems. The (possible) low-$z$ coupling suppression, its mechanism, and its consequences are however still open questions -not necessarily problems- for SCDEW models. The coupling intensity and the WDM particle mass, although being extra parameters in respect to $\Lambda$CDM, are found to be substantially constrained a priori so that, if SCDEW is the underlying cosmology, we expect most data to fit also $\Lambda$CDM predictions.

In this second paper we present the first Nbody cosmological simulations of strongly coupled Dark Energy models (SCDEW), a class of models that alleviates theoretical issues related to the nature of dark energy. SCDEW models assume a strong coupling between Dark Energy (DE) and an ancillary Cold Dark Matter (CDM) component together with the presence of an uncoupled Warm Dark Matter component. The strong coupling between CDM and DE allows us to preserve small scale fluctuations even if the warm particle is quite light ($\approx 100$ eV). Our large scale simulations show that, for $10^{11}<M/M_\odot<10^{14}$, SCDEW haloes exhibit a number density and distribution similar to a standard Lambda Cold Dark Matter (LCDM) model, even though they have lower concentration parameters. High resolution simulation of a galactic halo ($M\sim 10^{12} M_{\odot} $) shows $\sim 60\%$ less substructures than its LCDM counterpart, but the same cuspy density profile. On the scale of galactic satellites ($M\sim 10^{9} M_{\odot}$) SCDEW haloes dramatically differ from LCDM. Due to the high thermal velocities of the WDM component they are almost devoid of any substructures and present strongly cored dark matter density profiles. These density cores extend for several hundreds of parsecs, in very good agreement with Milky Way satellites observations. Strongly coupled models, thanks to their ability to match observations on both large and small scales might represent a valid alternative to a simple LCDM model.

Models including an energy transfer from CDM to DE are widely considered in the literature, namely to allow DE a significant high-z density. Strongly Coupled cosmologies assume a much larger coupling between DE and CDM, together with the presence of an uncoupled warm DM component, as the role of CDM is mostly restricted to radiative eras. This allows us to preserve small scale fluctuations even if the warm particle, possibly a sterile neutrino, is quite light, O(100 eV). Linear theory and numerical simulations show that these cosmologies agree with LCDM on supergalactic scales; e.g., CMB spectra are substantially identical. Simultaneously, simulations show that they significantly ease problems related to the properties of MW satellites and cores in dwarfs. SC cosmologies also open new perspectives on early black hole formation, and possibly lead towards unificating DE and inflationary scalar fields.

In the early Universe, a dual component made of coupled CDM and a scalar field $\Phi$, if their coupling $\beta > \sqrt{3}/2$, owns an attractor solution, making them a stationary fraction of cosmic energy during the radiation dominated era. Along the attractor, both such components expand $\propto a^{-4}$ and have early density parameters $\Omega_{d} = 1/ (4\beta^2)$ and $\Omega_c= 2, \Omega_d$ (field and CDM, respectively). In a previous paper it was shown that, if a further component, expanding $\propto a^{-3}$, breaks such stationary expansion at $z \sim 3$--$5 \times 10^3$, cosmic components gradually acquire densities consistent with observations. This paper, first of all, considers the case that this component is warm. However, its main topic is the analysis of fluctuation evolution: out of horizon modes are then determined; their entry into horizon is numerically evaluated as well as the dependence of Meszaros effect on the coupling $\beta$; finally, we compute: (i) transfer function and linear spectral function; (ii) CMB $C_l$ spectra. Both are close to standard $\Lambda$CDM models; in particular, the former one can be so down to a scale smaller than Milky Way, in spite of its main DM component being made of particles of mass $<1$ keV. The previously coupled CDM component, whose present density parameter is $\cal O$$(10^{-3})$, exhibits wider fluctuations $\delta \rho/\rho$, but approximately $\beta$-independent $\delta \rho$ values. We discuss how lower scale features of these cosmologies might ease quite a few problems that $\Lambda$CDM does not easily solve.

We assume the DE state equations w(a) = w_0+w_a(a_p-a), and study the dependence of the constraints on w_0 and w_a coefficients on the pivoting redshift 1+z_p=1/a_p. Coefficients are fitted to data including WMAP7, SNIa (Union 2.1), BAO's (including WiggleZ and SDSS results) and H_0 constraints. The fitting algorithm is CosmoMC. We find specific differences between the cases when neutrino mass is allowed or disregarded. More in detail: i) The z_p value yielding uncorrelated constraints on w_0 and w_a is different in the two cases, holding ~0.25 and ~0.35, respectively. (ii) If we consider the intervals allowed to w_0, we find that they shift when z_p increases, in opposite directions for vanishing or allowed neutrino mass. This leads to no overlap between 1sigma intervals already at z_p >~0.4. (iii) The known effect that a more negative state parameter is required to allow for neutrino mass displays its effects on w_a, rather than on w_0. (iv) The w_0-w_a constraints found by using any pivot z_p can be translated into constraints holding at a specific z_p value (0 or the z_p where errors are uncorrelated). When we do so, error ellipses exhibit a satisfactory overlap.

Forthcoming large angle surveys are planned to obtain high precision tomographic shear data. In principle, they will allow us to recover the spectra of matter density fluctuation, at various redshift, through the inversion of the expressions yielding shear from fluctuation spectra. This was discussed in previous work, where SVD techniques for matrix inversion were also shown to be the optimal tool to this aim. Here we show the significant improvements obtainable by using a 7 bin tomography, as allowed by future Euclid data, as well as the question of error propagation from shear to fluctuation spectra. We find that the technique is a promising tool, namely for the analysis of baryon physics throug high-l shear spectra and to test the consistency between expansion rate and fluctuation growth.

A dual component made of non-relativistic particles and a scalar field, exchanging energy, naturally falls onto an attractor solution, making them a (sub)dominant part of the cosmic energy during the radiation dominated era, provided that the constant \beta, measuring the coupling, is strong enough. The density parameters of both components are then constant, as they expand as a^{-4}. If the field energy is then prevalently kinetic, as is expected, its energy is exactly half of the pressureless component; the dual component as a whole, then, has a density parameter \Omega_{cd} = 3/4\beta^2 (e.g., for \beta~2.5, \Omega_{cd}~0.1, in accordance with Dark Radiation expectations). The stationary evolution can only be broken by the rising of other component(s), expanding as a^{-3}. In a realistic scenario, this happens when z~3-5x10^3. When such extra component(s) become(s) dominant, the densities of the dual components also rise above radiation. The scalar field behavior can be easily tuned to fit Dark Energy data, while the coupled DM density parameter becomes O(10^{-3}). This model however requires that, at present, two different DM components exist. The one responsible for the break of the stationary regime could be made, e.g., by thermally distributed particles with mass even >>1$-2 keV (or non-thermal particles with analogous average speed) so accounting for the size of observed galactic cores; in fact, a fair amount of small scale objects is however produced by fluctuation re-generated by the coupled DM component, in spite of its small density parameter, after the warm component has become non-relativistic.

Forthcoming experiments will enable us to determine tomographic shear spectra at a high precision level. Most predictions about them have until now been biased on algorithms yielding the expected linear and non-linear spectrum of density fluctuations. Even when simulations have been used, so-called Halofit (Smith et al 2003) predictions on fairly large scales have been needed. We wish to go beyond this limitation. We perform N-body and hydrodynamical simulations within a sufficiently large cosmological volume to allow a direct connection between simulations and linear spectra. While covering large length-scales, the simulation resolution is good enough to allow us to explore the high-l harmonics of the cosmic shear (up to l ~ 50000), well into the domain where baryon physics becomes important. We then compare shear spectra in the absence and in presence of various kinds of baryon physics, such as radiative cooling, star formation, and supernova feedback in the form of galactic winds. We distinguish several typical properties of matter fluctuation spectra in the different simulations and test their impact on shear spectra. We compare our outputs with those obtainable using approximate expressions for non--linear spectra, and identify substantial discrepancies even between our results and those of purely N-body results. Our simulations and the treatment of their outputs however enable us, for the first time, to obtain shear results taht are fully independent of any approximate expression, also in the high-l range, where we need to incorporate a non-linear power spectrum of density perturbations, and the effects of baryon physics. This will allow us to fully exploit the cosmological information contained in future high--sensitivity cosmic shear surveys, exploring the physics of cosmic shears via weak lensing measurements.

Forthcoming experiments will enable us to determine high precision tomographic shear spectra. Matter density fluctuation spectra, at various $z$, should then be worked out of them, in order to constrain the model and determine the DE state equation. Available analytical expressions, however, do the opposite, enabling us to derive shear spectra from fluctuation spectra. Here we find the inverse expression, yielding density fluctuation spectra from observational tomographic shear spectra. The procedure involves SVD techniques for matrix inversion. We show in detail how the approach works and provide a few examples.

At variance from 2H, 3He and 4He abundances, 7Li abundance data yield an extra constraint to cosmological parameters, in top of those deriving from CMB, BAO, SNIa or H0 data. This constraint, often disregarded, would favor smaller $\Omega_b h^2$ values, also indicating a preference for Dark Energy state equations well in the phantom regime, simultaneously softening the upper limit on the sum of neutrino masses, up to ~1.6eV.

We investigate the impact of non-linear corrections on dark energy parameter estimation from weak lensing probes. We find that using halofit expressions, suited to LCDM models, implies substantial discrepancies with respect to results directly obtained from N-body simulations, when w(z)\neq-1. Discrepancies appear strong when using models with w'(z=0)>0, as fiducial models; they are however significant even in the neighborhood of LCDM, where neglecting the degrees of freedom associated with the DE state equation can lead to a misestimate of the matter density parameter \Omega_m.

Mildly mixed coupled models include massive neutrinos and CDM--DE coupling. We present new tests of their likelihood vs. recent data including WMAP7, confirming it to exceed LCDM, although at ~2\sigma's. We then show the impact on the physics of the dark components of neutrino mass detection in tritium beta decay or neutrinoless double beta decay experiments.

The next generation mass probes will obtain information on non--linear power spectra P(k,z) and their evolution, allowing us to investigate the nature of Dark Energy. To exploit such data we need high precision simulations, extending at least up to scales of k 10 h/Mpc, where the effects of baryons can no longer be neglected. In this paper, we present a series of large scale hydrodynamical simulations for LCDM and dynamical Dark Energy (dDE) models, in which the equation of state parameter is z-dependent. The simulations include gas cooling, star formation and Supernovae feedback. They closely approximate the observed star formation rate and the observationally derived star/Dark Matter mass ratio in collapsed systems. Baryon dynamics cause spectral shifts exceeding 1% at k > 2-3 h/Mpc compared to pure n-body simulations in the LCDM simulations. This agrees with previous studies, although we find a smaller effect (~50%) on the power spectrum amplitude at higher k's. dDE exhibits similar behavior, even though the dDE simulations produce ~20% less stars than the analogous LCDM cosmologies. Finally, we show that the technique introduced in Casarini et al. to obtain spectra for any $w(z)$ cosmology from constant-w models at any redshift still holds when gas physics is taken into account. While this relieves the need to explore the entire functional space of dark energy state equations, we illustrate a severe risk that future data analysis could lead to misinterpretation of the DE state equation.

The next generation mass probes will investigate DE nature by measuring non-linear power spectra at various z, and comparing them with high precision simulations. Producing a complete set of them, taking into account baryon physics and for any DE state equation w(z), would really be numerically expensive. Regularities reducing such duty are essential. This paper presents further n-body tests of a relation we found, linking models with DE state parameter w(z) to const.-w models, and also tests the relation in hydro simulations.

We consider cosmological models with dynamical dark energy (dDE) coupled to cold dark matter (CDM), while simultaneously allowing neutrinos to be massive. Using a MCMC approach, we compare these models with a wide range of cosmological data sets. We find a strong correlation between this coupling strength and the neutrino mass. This correlation persists when BAO data are included in the analysis. We add then priors on $\nu$ mass from particle experiments. The claimed detection of $\nu$ mass from the Heidelberg-Moscow neutrinoless double--$\beta$ decay experiment would imply a 7--$8 \sigma$ detection of CDM-DE coupling. Similarly, the detection of $\nu$ mass from coming KATRIN tritium $\beta$ decay experiment will imply a safe detection of a coupling in the dark sector. Previous attempts to accommodate cosmic phenomenology with such possible $\nu$ mass data made recourse to a $w < -1$ eoS. We compare such an option with the coupling option and find that the latter allows a drastic improvement.

We fit WMAP5 and related data by allowing for a CDM--DE coupling and non--zero neutrino masses, simultaneously. We find a significant correlation between these parameters, so that simultaneous higher coupling and \nu--masses are allowed. Furthermore, models with a significant coupling and \nu--mass are statistically favoured in respect to a cosmology with no coupling and negligible neutrino mass (our best fits are: C ~ 1/2m_p, m_\nu ~ 0.12eV per flavor). We use a standard Monte Carlo Markov Chain approach, by assuming DE to be a scalar field self--interacting through Ratra--Peebles or SUGRA potentials.

We allow simultaneously for a CDM--DE coupling and non--zero neutrino masses and find that significant coupling and neutrino mass are (slightly) statistically favoured in respect to a cosmology with no coupling and negligible neutrino mass (our best fits are: C~1/2m_p, m_\nu~0.12eV each flavor). We assume DE to be a self--interacting scalar field and use a standard Monte Carlo Markov Chain approach.

Accurate predictions on non--linear power spectra, at various redshift z, will be a basic tool to interpret cosmological data from next generation mass probes, so obtaining key information on Dark Energy nature. This calls for high precision simulations, covering the whole functional space of w(z) state equations and taking also into account the admitted ranges of other cosmological parameters; surely a difficult task. A procedure was however suggested, able to match the spectra at z=0, up to k~3, hMpc^{-1}, in cosmologies with an (almost) arbitrary w(z), by making recourse to the results of N-body simulations with w = const. In this paper we extend such procedure to high redshift and test our approach through a series of N-body gravitational simulations of various models, including a model closely fitting WMAP5 and complementary data. Our approach detects w= const. models, whose spectra meet the requirement within 1% at z=0 and perform even better at higher redshift, where they are close to a permil precision. Available Halofit expressions, extended to (constant) w \neq -1 are unfortunately unsuitable to fit the spectra of the physical models considered here. Their extension to cover the desired range should be however feasible, and this will enable us to match spectra from any DE state equation.

Within the frame of cosmologies where Dark Energy (DE) is a self--interacting scalar field, we allow for a CDM--DE coupling and non--zero neutrino masses, simultaneously. In their 0--0 version, i.e. in the absence of coupling and neutrino mass, these cosmologies provide an excellent fit to WMAP, SNIa and deep galaxy sample spectra, at least as good as \LambdaCDM. When the new degrees of freedom are open, we find that CDM--DE coupling and significant neutrino masses (~0.1eV per \nu species) are at least as likely as the 0--0 option and, in some cases, even statistically favoured. Results are obtained by using a Monte Carlo Markov Chain approach.

Cosmological limits on neutrino masses are softened, by more than a factor 2, if Cold Dark Matter (CDM) and Dark Energy (DE) are coupled. In turn, a neutrino mass yielding $\Omega_\nu$ up to $\sim0.20$ allows coupling levels $\beta \simeq 0.15, $ or more, already easing the coincidence problem. The coupling, in fact, displaces both $P(k)$ and $C_l$ spectra in a fashion opposite to neutrino mass. Estimates are obtained through a Fisher--matrix technique.

When CMB data are used to derive cosmological parameters, their very choice does matter: some parameter values can be biased if the parameter space does not cover the "true" model. This is a problem, because of the difficulty to parametrize Dark Energy (DE) physics. We test this risk through numerical experiments. We create artificial data for dynamical or coupled DE models and then use MCMC techniques to recover model parameters, by assuming a constant DE state parameter w and no DM--DE coupling. For the DE potential considered, no serious bias arises when coupling is absent. On the contrary, \omega_{o,c}, and thence H_o and \Omega_{o,m}, suffer a serious bias when the "true" cosmology includes even just a mild DM--DE coupling. Until the dark components keep an unknown nature, therefore, it can be important to allow for a degree of freedom accounting for DM--DE coupling, even more than increasing the number of parameters accounting for the w(a) behavior.

A phaenomenological DM-DE coupling could indicate their common origin. Various constraint however exist to such coupling; here we outline that it can suppress Meszaros' effect, yielding transfered spectra with a softer bending above k_{hor,eq}. It could be therefore hard to reconcile these models with both CMB and deep sample data, using a constant spectral index.

The dual axion model (DAM), yielding bot DM and DE form a PQ-like scalar field solving the strong CP problem, is known to allow a fair fit of CMB data. Recently, however, it was shown that its transfer function exhibits significant anomalies, causing difficulties to fit deep galaxy sample data. Here we show how DAM can be modified to agree with the latter data set. The modification follows the pattern suggested to reconcile any PQ-like approach with gravity. Modified DAM allows precise predictions which can be testable against future CMB and/or deep sample data.

DM--DE coupling can be a phenomenological indication of a common origin of the dark cosmic components. In this work we outline a new constraint to coupled--DE models: the coupling can partially or totally suppress the Meszaros effect, yielding transfered spectra with quite a soft bending above $k_{hor,eq}$. Models affected by this anomaly do not show major variation in the CMB anisotropy spectrum and it is herefore hard to reconcile them with both CMB and deep sample data, through the same value of the primeval spectral index.

We compare a large set of cosmologies with WMAP data, performing a fit based on a MCMC algorithm. Besides of LCDM models, we take dynamical DE models, where DE and DM are uncoupled or coupled, both in the case of constant coupling and in the case when coupling varies with suitable laws. DE however arises from a scalar field self-interacting through a SUGRA potential. We find that the best fitting model is SUGRA dynamical DE, almost indipendently from the exponent alpha in the self-interaction potential. The main target of this work are however coupled DE models, for which we find limits on the DE-DM coupling strength. In the case of variable coupling, we also find that greater values of the Hubble constant are preferred.

Although the cosmic concordance cosmology is quite successful in fitting data, fine tuning and coincidence problems apparently weaken it. We review several possibilities to ease its problems, by considering various kinds of dynamical Dark Energy and possibly its coupling to Dark Matter, trying to set observational limits on Dark Energy state equation and coupling.

The strong CP problem was solved by Peccei & Quinn by introducing axions, a viable candidate for Dark Matter (DM). Here the PQ approach is modified so to yield also Dark Energy (DE). DM and DE arise, in fai proportions, from a single scalar field, without tuning any extra parameter. In the present epoch, they are weakly coupled. Fluctuations have a fair evolution. The model is also fitted to the WMAP1 release, using a Markov Chain Monte Carlo technique, and performs as well as $\Lambda$CDM, coupled or uncoupled DE. Best--fit cosmological parameters for different models are mostly within 2--$\sigma$ level. Here, the main peculiarity of the model is to favor high values of the Hubble parameter.

We evaluate the mass function of virialized halos, by using Press & Schechter (PS) and/or Steth & Tormen (ST) expressions, for cosmologies where Dark Energy (DE) is due to a scalar self-interacting field, coupled with Dark Matter (DM). We keep to coupled DE (cDE) models known to fit linear observables. To implement the PS-ST approach, we start from reviewing and extending the results of a previous work on the growth of a spherical top-hat fluctuation in cDE models, confirming their most intriguing astrophysical feature, i.e. a significant baryon-DM segregation, occurring well before the onset of any hydrodynamical effect. Accordingly, the predicted mass function depends on how halo masses are measured. For any option, however, the coupling causes a distortion of the mass function, still at z=0. Furthermore, the z-dependence of cDE mass functions is mostly displaced, in respect to LambdaCDM, in the opposite way of uncoupled dynamical DE. This is an aspect of the basic underlying result, that even a little DM-DE coupling induces relevant modifications in the non-linear evolution. Therefore, without causing great shifts in linear astrophysical observables, the DM-baryon segregation induced by the coupling can have an impact on a number of cosmological problems, it e.g., galaxy satellite abundance, spiral disk formation, apparent baryon shortage, entropy input in clusters, etc..

Dynamical Dark Energy (DE) is a viable alternative to the cosmological constant. Yet, constructing tests to discriminate between Lambda and dynamical DE models is difficult because the differences are not large. In this paper we explore tests based on the galaxy mass function, the void probability function (VPF), and the number of galaxy clusters. At high z the number density of clusters shows large differences between DE models, but geometrical factors reduce the differences substantially. We find that detecting a model dependence in the cluster redshift distribution is a hard challenge. We show that the galaxy redshift distribution is potentially a more sensitive characteristics. We do so by populating dark matter halos in Nbody simulations with galaxies using well-tested Halo Occupation Distribution (HOD). We also estimate the Void Probability Function and find that, in samples with the same angular surface density of galaxies in different models, the VPF is almost model independent and cannot be used as a test for DE. Once again, geometry and cosmic evolution compensate each other. By comparing VPF's for samples with fixed galaxy mass limits, we find measurable differences.

Axions are likely to be the Dark Matter (DM) that cosmological data require. They arise in the Peccei Quinn solution of the strong CP problem. In a previous work we showed that their model has a simple and natural generalization which yields also Dark Energy (DE), in fair proportions, without tuning any parameter: DM and DE arise from a single scalar field and are weakly coupled in the present era. In this paper we extend the analysis of this dual-axion cosmology and fit it to WMAP data, by using a Markov chain technique. We find that LCDM, dynamical DE with a SUGRA potential, DE with a SUGRA potential and a constant DE-DM coupling, as well as the dual-axion model with a SUGRA potential, fit data with a similar accuracy. The best-fit parameters are however fairly different, although consistency is mostly recovered at the 2-sigma level. A peculiarity of the dual-axion model with SUGRA potential is to cause more stringent constraints on most parameters and to favor high values of the Hubble parameter.

Neutral hydrogen around high-z QSO and an optical depth tau ~ 0.17 can be reconciled if reionization is more complex than a single transition at z ~ 6-8. Tracing its details could shed a new light on the first sources of radiation. Here we discuss how far such details can be inspected through planned experiments on CMB large-scale anisotropy and polarization, by simulating an actual data analysis. By considering a set of double reionization histories of Cen (2003) type, a relevant class of models not yet considered by previous works, we confirm that large angle experiments rival high resolution ones in reconstructing the reionization history. We also confirm that reionization histories, studied with the prior of a single and sharp reionization, yield a biased tau, showing that this bias is generic. We further find a monotonic trend in the bias for the models that we consider, and propose an explanation of the trend, as well as the overall bias. We also show that in long-lived experiments such a trend can be used to discriminate between single and double reionization patterns.

We show that using the redshift dependence of the deep cluster abundance to detect the nature of Dark Energy is a serious challenge. We compare the expected differences between flat LCDM models, with different Omega_mo, with the difference between LCDM and dynamical DE models. In the former case, cluster abundances in comoving volume and geometrical factors act in the same direction, yielding a significant difference between the expected angular densities. On the contrary, when we keep a constant Omega_mo and change the DE nature,abundances in comoving volume and geometrical factors act in the opposite direction, so that the expected differences in angular densities reduce to small factors.

The Peccei Quinn (PQ) solution of the strong CP problem requires the existence of axions, which are a viable candidate for Dark Matter. Here we show that, if the Nambu Goldstone potential of the PQ model is replaced by a potential V(|\Phi|) admitting a tracker solution, the scalar field |\Phi| can account for Dark Energy, while the phase of \Phi yields axion Dark Matter. Such Dark Matter and Dark Energy turn out to be weakly coupled. If V is a SUGRA potential, the model essentially depends on a single parameter, the energy scale \Lambda. Once we set \Lambda \simeq 10^{10} GeV, at the quark--hadron transition, |\Phi| naturally passes through values suitable to solve the strong CP problem, later growing to values providing fair amounts of Dark Matter and Dark Energy. In this model, the linear growth factor, from recombination to now, is quite close to \LambdaCDM. The selected \Lambda value can be an indication of the scale where the soft breaking of SUSY occurred.

SPOrt is an ASI-funded experiment specifically designed to measure the sky polarization at 22, 32 and 90 GHz, which was selected in 1997 by ESA to be flown on the International Space Station. Starting in 2006 and for at least 18 months, it will be taking direct and simultaneous measurements of the Stokes parameters Q and U at 660 sky pixels, with FWHM=7 degrees. Due to development efforts over the past few years, the design specifications have been significantly improved with respect to the first proposal. Here we present an up-to-date description of the instrument, which now warrants a pixel sensitivity of 1.7 microK for the polarization of the cosmic background radiation, assuming two years of observations. We discuss SPOrt scientific goals in the light of WMAP results, in particular in connection with the emerging double-reionization cosmological scenario.

We present the results of a series of N-body simulations in cosmologies where dark matter (DM) is coupled to dark energy (DE), so easing the cosmic coincidence problem. The dark-dark coupling introduces two novel effects in N-body dynamics: (i) DM particle masses vary with time; (ii) gravity between DM particles is ruled by a constant $G^{*}$, greater than Newton's constant $G$, holding in other 2-body interactions. As a consequence, baryons and DM particle distributions develop a large scale bias. Here we investigate DE models with Ratra-Peebles (RP) potentials; the dark-dark coupling is set in a parametric range compatible with observations, for as concern background and linear perturbation properties. We study the halo mass function, the halo density profile and the behavior of the non-linear bias. We find that non-linear dynamics puts additional constraints to the coupling parameter. They mostly arise from density profiles, that we find to yield higher concentrations, in coupled RP models, with respect to (uncoupled) dynamical DE cosmologies. Such enhancement, although being a strong effect in some coupling parameter range, is just a minor change for smaller but significant values of the coupling parameter. With these further restrictions, coupled DE models with RP potential are consistent with non-linear observables.

Polarization measures, on wide angular scales, together with anisotropy data, can fix DE parameters. Here we discuss the sensitivity needed to provide significant limits. Our analysis puts in evidence that a class of models predicts low correlation or anticorrelation between polarization and anisotropy at low l. This class includes open models and models with DE due to a Ratra--Peebles (RP) potential. Results on this point, given in a previous paper of ours, are updated and partially corrected. We outline that, with the sensitivity of experiments like SPOrt or WMAP, high values of Lambda (energy scale in the RP potential) can be excluded. With the sensitivity expected for PLANCK, the selection will extend to much lower Lambda's.

We perform N-body simulations for models with a DE component. Besides of DE with constant negative state parameter w, we consider DE due to scalar fields, self-interacting through RP or SUGRA potentials. According to our post-linear analysis, at z=0, DM power spectra and halo mass functions do not depend on DE nature. This is welcome, as LCDM fits observations. Halo profiles, instead, are denser than LCDM. For example, the density at 15 kpc of a DE halo with M=10^13 exceeds LCDM by ~45%. Differences, therefore, are small but, however, DE does not ease the problem with cuspy DM profiles. We study also subhalos and find that, at $z=0$, the number of satellites coincides in all DE models. At higher z, DE models show increasing differences from LCDM and among themselves (i.e. in the mass function evolution); this is the obvious pattern to distinguish between different DE state equations.

High sensitivity polarization measures, on wide angular scales, together with data on anisotropy, can be used to fix DE parameters. In this paper, first of all, we aim to determine the sensitivity needed to provide significant limits. Our analysis puts in evidence that there is a class of DE models that polarization measures can possibly exclude soon. This class includes models with DE due to a Ratra-Peebles (RP) potential. Using a likelihood analysis, we show that it is possible to distinguish RP models from LCDM and other dynamical DE models, already with the sensitivity of experiments like SPOrt or WMAP, thanks to their negative TE correlation at low-l, when the optical depth tau is sufficiently large. On the contrary, fixing the energy scale Lambda for RP potentials or distinguishing between LCDM and other DE potentials requires a much lower pixel noise, that no planned polarization experiment will achieve. While reviewing this paper after the referee report, the first-year WMAP data were released. WMAP finds large positive anisotropy-polarization correlations at low l; this apparently excludes DE models with RP potentials.

Cosmological models with different types of Dark Energy are becoming viable alternatives for standard models with the cosmological constant. Yet, such models are more difficult to analyze and to simulate. We present analytical approximations and discuss ways of making simulations for two families of models, which cover a wide range of possibilities and include models with both slow and fast changing ratio w=p\rho. More specifically, we give analytical expressions for the evolution of the matter density parameter Omega_m(z) and the virial density contrast Delta_c at any redshift z. The latter is used to identify halos and to find their virial masses. We also provide an approximation for the linear growth factor of linear fluctuations between redshift z=40 and z=0. This is needed to set the normalization of the spectrum of fluctuations. Finally, we discuss the expected behavior of the halo mass function and its time evolution.

We study properties of dark matter halos in a variety of models which include Dark Energy (DE). We consider both DE due to a scalar field self--interacting through Ratra-Peebles or SUGRA potentials, and DE with constant negative w=p\rho >-1. We find that at redshift zero the nonlinear power spectrum of the dark matter, and the mass function of halos, practically do not depend on DE state equation and are almost indistinguishable from predictions of the LCDM model. This is consistent with the nonlinear analysis presented in the accompanying paper. It is also a welcome feature because LCDM models fit a large variety of data. On the other hand, at high redshifts DE models show substantial differences from LCDM and substantial differences among themselves. Halo profiles differ even at z=0. DE halos are denser than LCDM in their central parts because the DE halos collapse earlier. Nevertheless, differences between the models are not so large. For example, the density at 10 kpc of a DE \~10^{13}Msun halo deviates from LCDM by not more than 50%. This, however, means that DE is not a way to ease the problem with cuspy dark matter profiles. Addressing another cosmological problem - abundance of subhalos -- we find that the number of satellites of halos in various DE models does not change relative to the LCDM, when normalized to the same circular velocity of the parent halo. To summarize, the best way to find which DE model fits the observed Universe is to look for evolution of halo properties. For example, the abundance of galaxy groups with mass larger than 10^{13}\Msun at z> 2 can be used to discriminate between the models, and, thus, to constrain the nature of DE.

We study the possibility of mixed dark matter obtained through the phase space distribution of a single particle. An example is offered in the context of SUSY models with a Peccei-Quinn symmetry. Axinos in the 100 keV range can naturally have both thermal and non-thermal components. The latter one arises from the lightest neutralino decays and derelativizes at z ~ 10^4.

We show that DM halos, in n--body simulations, have a boundary layer (BL), separating bound from unbound mass. Let T(r) and W(r) be the kinetic and potential energies in shells of halos. We find that, in almost all halos: (i) The virial ratio R=-2T/W has at least one persistent (i.e. resolution independent) minimum r_c, such that, close to it, (ii) the function w=-d log W/d log r has a maximum and (iii) the relation R(r_c)=w(r_c) is almost exactly fulfilled. The radius r_c is the position of BL's, in halos found in simulations of TCDM and LCDM models, run using ART and GADGET codes at various resolutions. We find that 97% of the ~300 halos (per model) with M>4.2 10^14 M_s h^-1 owns a BL. Those with no BL are undergoing major mergings. The mass M_c enclosed in r_c almost coincides with the mass evaluated from velocities, according to virial theorem. Particles at r>r_c are not in virial equilibrium. Using r_c we have a density contrast \Delta_c for each halo. For each mass scale, \Delta_v=178 Omega_m^0.45 is within the range of \Delta_c's, but the spread in \Delta_c is wide and the average \Delta_c is ~25% smaller than the corresponding \Delta_v. The matching of properties derived under the assumption of spherical symmetry is a consequence of violent relaxation destroying features related to previous ellipsoidal non-linear growth. In turn, the spread of final \Delta_c's is an imprint of the different 3-D geometries and of the variable environment during collapses, as suggested by a comparison with Sheth & Tormen analysis.

BaR-SPOrt (Balloon-borne Radiometers for Sky Polarisation Observations) is an experiment to measure the linearly polarized emission of sky patches at 32 and 90 GHz with sub-degree angular resolution. It is equipped with high sensitivity correlation polarimeters for simultaneous detection of both the U and Q stokes parameters of the incident radiation. On-axis telescope is used to observe angular scales where the expected polarization of the Cosmic Microwave Background (CMBP) peaks. This project shares most of the know-how and sophisticated technology developed for the SPOrt experiment onboard the International Space Station. The payload is designed to flight onboard long duration stratospheric balloons both in the Northern and Southern hemispheres where low foreground emission sky patches are accessible. Due to the weakness of the expected CMBP signal (in the range of microK), much care has been spent to optimize the instrument design with respect to the systematics generation, observing time efficiency and long term stability. In this contribution we present the instrument design, and first tests on some components of the 32 GHz radiometer.

SPOrt (Sky Polarization Observatory) is a space experiment to be flown on the International Space Station during Early Utilization Phase aimed at measuring the microwave polarized emission with FWHM = 7deg, in the frequency range 22-90 GHz. The Galactic polarized emission can be observed at the lower frequencies and the polarization of Cosmic Microwave Background (CMB) at 90 GHz, where contaminants are expected to be less important. The extremely low level of the CMB Polarization signal (< 1 uK) calls for intrinsically stable radiometers. The SPOrt instrument is expressly devoted to CMB polarization measurements and the whole design has been optimized for minimizing instrumental polarization effects. In this contribution we present the receiver architecture based on correlation techniques, the analysis showing its intrinsic stability and the custom hardware development carried out to detect such a low signal.

We study the cluster mass function and its evolution in different models with Dark Energy arising from a self--interacting scalar field, with Ratra-Peebles and SUGRA potentials. Computations are based on a Press & Schechter approximation. The mass functions we obtain are compared with results holding for open models or models with Dark Energy due to a cosmological constant. Evolution results, in the Dark Energy cases, closely approach open models.

Anisotropy data analysis leaves a significant degeneracy between primeval spectral index (n_s) and cosmic opacity to CMB photons (\tau). Low--l polarization measures, in principle, can remove it. We perform a likelihood analysis to see how cosmic variance possibly affects such a problem. We find that, for a sufficiently low noise level (\sigma_{pix}) and if \tau is not negligibly low, the degeneracy is greatly reduced, while the residual impact of cosmic variance on n_s and \tau determinations is under control. On the contrary, if \sigma_{pix} is too high, cosmic variance effects appear to be magnified. We apply general results to specific experiments and find that, if favorable conditions occur, it is possible that a 2--\sigma detection of a lower limit on \tau is provided by the SPOrt experiment. Furthermore, if the PLANCK experiment will measure polarization with the expected precision, the error on low--l harmonics is adequate to determine \tau, without significant magnification of the cosmic variance. This however indicates that high sensitivity might be more important than high resolution in \tau determinations. We also outline that a determination of \tau is critical to perform detailed analyses on the nature of dark energy and/or on the presence of primeval gravitational waves.

We show that dark matter haloes, in n--body simulations, have a boundary layer (BL) with precise features. In particular, it encloses all dynamically stable mass while, outside it, dynamical stability is lost soon. Particles can pass through such BL, which however acts as a confinement barrier for dynamical properties. BL is set by evaluating kinetic and potential energies (T(r) and W(r)) and calculating R=-2T/W. Then, on BL, R has a minimum which closely approaches a maximum of w= -dlog W/dlog r. Such $Rw$ ``requirement'' is consistent with virial equilibrium, but implies further regularities. We test the presence of a BL around haloes in spatially flat CDM simulations, with or without cosmological constant. We find that the mass M_c, enclosed within the radius r_c, where the $Rw$ requirement is fulfilled, closely approaches the mass M_{dyn}, evaluated from the velocities of all particles within r_c, according to the virial theorem. Using r_c we can then determine an individual density contrast Delta_c for each virialized halo, which can be compared with the "virial" density contrast $\Delta_v ~178 \Omega_m^{0.45}$ (Omega_m: matter density parameter) obtained assuming a spherically symmetric and unperturbed fluctuation growth. The spread in Delta_c is wide, and cannot be neglected when global physical quantities related to the clusters are calculated, while the average Delta_c is ~25 % smaller than the corresponding Delta_v; moreover if $M_{dyn}$ is defined from the radius linked to Delta_v, we have a much worse fit with particle mass then starting from {\it Rw} requirement.

Using COBE-DMR 4-year data, we find a general expression yielding the likelihood distribution in the 3-dimensional parameter space spanned by the spectral index $n$, the spectral amplitude $a_{10}$ and the false-vacuum density parameter $\Omega_\Lambda$. Using such simple expression, the range of possible normalizations, within a given likelihood interval from top-likelihood normalization, is readily found, with fair approximation, for any model with total density parameter $\Omega_o = 1$ and assigned $n$ and $\Omega_\Lambda$.

We examine how the statistics of the quadrupoles of (projected) cluster masses can discriminate between flat cold dark matter (CDM) universes with or without a cosmological constant term. Even in the era of high precision cosmology that cosmic microwave background experiments should open soon, it is important to devise self consistency tests of cosmogonic theories tuned at the matter radiation decoupling epoch using data from the non--linear evolved universe. We build cluster catalogs from two large volume simulations of a ``tilted'' CDM model and a $\Lambda$CDM model with cosmic density parameter $\Omega_m=0.35$ and cosmological constant contribution $\Omega_\Lambda=0.65$. From the projected mass distribution of the clusters we work out the quadrupoles Q and examine their dependence on cluster mass and the cosmological model. We find that TCDM clusters have systematically larger quadrupoles than their \lcdm counterpart. The effect is mass dependent: massive clusters ($M\gsim 10^{15}\Msunh$) have quadrupoles differing by more than 30% in the two models, while for $M\lsim 4\times10^{14}\Msunh$ the difference rapidly drops to $\sim 1$%. Performing a K-S test of the Q distributions, we estimate that using just the 15 most massive clusters in the simulation volume ($360\Mpch$ a side) we can discriminate between TCDM and \lcdm at a confidence level better than 99.9%. In the volume probed by exhisting observations, there are potentially several hundred clusters with masses above the threshold for which the differences in the quadrupoles become relevant. Should weak lensing data become available for this whole set, a quadrupole analysis may be expected to discriminate among different values of $\Lambda$.

We evaluate the dependence of the cluster correlation length r_c on the mean intercluster separation D_c, for three models with critical matter density, vanishing vacuum energy (Lambda = 0) and COBE normalized: a tilted CDM (tCDM) model (n=0.8) and two blue mixed models with two light massive neutrinos yielding Omega_h = 0.26 and 0.14 (MDM1 and MDM2, respectively). All models approach the observational value of sigma_8 (and, henceforth, the observed cluster abundance) and are consistent with the observed abundance of Damped Lyman_alpha systems. Mixed models have a motivation in recent results of neutrino physics; they also agree with the observed value of the ratio sigma_8/sigma_25, yielding the spectral slope parameter Gamma, and nicely fit LCRS reconstructed spectra. We use parallel AP3M simulations, performed in a wide box (side 360/h Mpc) and with high mass and distance resolution, enabling us to build artificial samples of clusters, whose total number and mass range allow to cover the same D_c interval inspected through APM and Abell cluster clustering data. We find that the tCDM model performs substantially better than n=1 critical density CDM models. Our main finding, however, is that mixed models provide a surprisingly good fit of cluster clustering data.

A large set of TREESPH simulations is used to test the global morphology of galaxy clusters and its evolution against X-ray data. A powerful method to investigate substructures in galaxy clusters are the power ratios introduced by Buote & Tsai. We consider three flat cosmological models: CDM, LCDM($\Omega_\Lambda = 0.7$) and CHDM($\Omega_h = 0.2$, 1 massive $\nu$), all normalized so to fit the observed number of clusters. For each model we built 40 clusters, using a TREESPH code, and performed a statistical comparison with a data sample including nearby clusters observed with ROSAT PSPC instrument. The comparison disfavors the LCDM model, as clusters appear too relaxed, while CDM and CHDM clusters, in which a higher degree of complexity occurs, seem to be closer to observations. A better fit of data can be expected for some different DM mix. If DM distributions are used instead of baryons, we find substructures more pronounced than in gas and models have a different score. Using hydrodynamical simulations is therefore essential to our aims.

The passage from angular to spatial correlations, in the case of spatial clustering length depending on the average distance between nearby objects is studied. We show that, in a number of cases, the scaling law of angular correlation amplitudes, which holds for constant spatial clustering length, is still true also for a luminosity dependent spatial correlation. If the Limber equation is then naively used to obtain `the' spatial clustering length from the angular function amplitude, a quantity close to the average object separation is obtained. The case of cluster clustering is explicitly considered.

We study the cluster mass function in mixed dark matter (MDM) models, using two COBE normalized simulations with Omega_h = 0.26 and n=1.2, and Omega_h = 0.14 and n = 1.05, both with 2 massive nu's (MDM1 and MDM2, respectively). For the sake of comparison, we also simulate a CDM model with spectral index n=0.8 (TCDM), also COBE normalized. We argue that, in our non--hydro simulations, where CDM particles describe both actual CDM and baryons, the galaxy distribution traces CDM particles. Therefore, we use them to define clusters and their velocities to work out cluster masses. As CDM particles are more clustered than HDM and therefore have, in average, greater velocities, this leads to significant differences from PS predictions. Such predictions agree with simulations if both HDM and CDM are used to define clusters. Clusters defined through CDM in MDM models, instead, are less numerous than PS estimates, by a factor ~0.3, at the low mass end; the factor becomes \~0.6-0.8, depending on the mix, on intermediate mass scales (~4-5, h^-1 10^14 Msun) and almost vanishes on the high mass end. Therefore: (i) MDM models expected to overproduce clusters over intermediate scales are viable; (ii) the greater reduction factor at small scales agrees with the observational data dependence on the cluster mass M (which, however, may be partially due to sample incompleteness); (iii) the higher spectral normalization is felt at large scales, where MDM models produce more objects (hence, large clusters) than CDM. MDM1 even exceeds Donahue et al. (1998) findings, while MDM2 is consistent with them. (abridged)

The global structure of galaxy clusters and its evolution are tested within a large set of TREESPH simulations, so to allow a fair statistical comparison with available X-ray data. Structure tests are based on the "power ratios", introduced by Buote & Tsai. Cosmological models considered are CDM, LCDM (Omega_L=0.7) and CHDM (1 mass.neu., Omega_h = 0.2). All models are normalized to provide a fair number density of clusters. For each model we run a P3M simulation in a large box, where we select the most massive 40 clusters. Going back to the initial redshift we run a hydro-TREESPH simulation for each of them. In this way we perform a statistical comparison of the global morphology of clusters, for each cosmological model, with ROSAT data, using Student t-test, F-test and K-S test. The last test and its generalization to 2--D distributions are also used to compare the joint distributions of 2 or 3 power ratios. We find that, using DM distribution, instead of gas, as done by some authors, leads to biased results, as baryons are distributed in a less structured way than DM. We also find that the cosmological models considered have different behaviours in these tests: LCDM has the worst performance. CDM and our CHDM have similar scores. The general trend of power ratio distributions is already fit by these models, but a further improvement is expected either from a different DM mix or a non-flat CDM model.

We compute CBR anisotropies in mixed models with different hot components, including neutrinos or volatile HDM arising from the decay of heavier particles. The CBR power spectra of these models exhibit a higher doppler peak than CDM, and the discrepancy is even stronger in volatile models when the decay gives rise also to a neutral scalar. CBR experiments, together with Large Scale Structure (LSS) data, are then used to constrain the space parameter of mixed models, when values of the primeval spectral index n > 1 are also considered. Even if n>1 is allowed, however, LSS alone prescribes that \Omega_h \mincir 0.30. LSS can be fitted by taking simultaneously a low derelativization redshift z_der (down to \simeq 600) and a high n, while CBR data from baloon--borne experiment cause a severe selection on this part of the parameter space. In fact, while late derelativization and $n>1$ have opposite effects on the fluctuation spectrum P(k), they sum their action on the angular spectrum C_l. Henceforth n \magcir 1.3 seems excluded by baloon--borne experiment outputs, while a good fit of almost all CBR and LSS data is found for Omega_h values between 0.11 and 0.16, n \sim 1.1 and z_der \sim 2000--5000. A smaller n is allowed, but z_der should never be smaller than \simeq 1200.

Recent data on CBR anisotropies show a Doppler peak higher than expected in CDM cosmological models, if the spectral index $n=1$. However, CDM and LCDM models with n>1 can hardly be consistent with LSS data. Mixed models, instead, whose transfer function is naturally steeper because of free--streaming in the hot component, may become consistent with data if n>1, when Omega_h is large. This is confirmed by our detailed analysis, extended both to models with a hot component whose momentum space distribution had a thermal origin (like massive neutrinos), and to models with a non--cold component arising from heavier particle decay. In this work we systematically search models which fulfill all constraints which can be implemented at the linear level. We find that a stringent linear constraint arises from fitting the extra-power parameter Gamma. Other significant constraints arise comparing the expected abundances of galaxy clusters and high-z systems with observational data. Keeping to models with Gamma \geq 0.13, a suitable part of the space parameter still allows up to \sim 30% of hot component (it is worth outlining that our stringent criteria allow only models with 0.10 \mincir Omega_h \mincir 0.16, if n \leq 1). We also outline that models with such large non--cold component would ease the solution of the so--called baryon catastrophe in galaxy clusters.

We compare the evolution of morphology in X-ray cluster data and in clusters obtained in a simulation of flat CDM, normalized to the observed cluster abundance. We find that the evolution rate in model clusters is significantly higher than in data clusters. The results is stricking, as our cluster data are all just for $z \simlt 0.2$, but is not contrary to expectations, as is known that the CDM spectrum is too flat, around the cluster scale, and therefore induces the presence of too many substructures. The cluster simulations were run using a TREESPH code and therefore include hydrodynamical effects. The test we performed, which turns out to be very sensitive, starts from the so--called power ratios introduced by Buote and Tsai, and makes use of the 2--dimensional generalization of the Kolmogorov--Smirnov test. The discrepancy between data and model is however directly visible through linear regression techniques and using simpler statistical tests.

We apply several statistical estimators to high-resolution N-body simulations of two currently viable cosmological models: a mixed dark matter model, having $\Omega_\nu=0.2$ contributed by two massive neutrinos (C+2\nuDM), and a Cold Dark Matter model with Cosmological Constant (\LambdaCDM) with $\Omega_0=0.3$ and h=0.7. Our aim is to compare simulated galaxy samples with the Perseus-Pisces redshift survey (PPS). We consider the n-point correlation functions (n=2-4), the N-count probability functions P_N, including the void probability function P_0, and the underdensity probability function U_\epsilon (where \epsilon fixes the underdensity threshold in percentage of the average). We find that P_0 (for which PPS and CfA2 data agree) and P_1 distinguish efficiently between the models, while U_\epsilon is only marginally discriminatory. On the contrary, the reduced skewness and kurtosis are, respectively, S_3\simeq 2.2 and S_4\simeq 6-7 in all cases, quite independent of the scale, in agreement with hierarchical scaling predictions and estimates based on redshift surveys. Among our results, we emphasize the remarkable agreement between PPS data and C+2\nuDM in all the tests performed. In contrast, the above \LambdaCDM model has serious difficulties in reproducing observational data if galaxies and matter overdensities are related in a simple way.

We investigate the clustering properties of loose groups in the Perseus--Pisces redshift Survey (PPS). Previous analyses based on CfA and SSRS surveys led to apparently contradictory results. We investigate the source of such discrepancies, finding satisfactory explanations for them. Furthermore, we find a definite signal of group clustering, whose amplitude $A_G$ exceeds the amplitude $A_g$ of galaxy clustering ($A_G=14.5^{+3.8}_{-3.0}$, $A_g=7.42^{+0.20}_{-0.19}$ for the most significant case; distances are measured in $\hMpc$). Groups are identified with the adaptive Friends--Of--Friends (FOF) algorithms HG (Huchra \& Geller 1982) and NW (Nolthenius \& White 1987), systematically varying all search parameters. Correlation strenght is especially sensitive to the sky--link $D_L$ (increasing for stricter normalization $D_0$), and to the (depth $\mlim$ of the) galaxy data. It is only moderately dependent on the galaxy luminosity function $\phi(L)$, while it is almost insensitive to the redshift--link $V_L$ (both to the normalization $V_0$ and to the scaling recipes HG or NW).

We use the void probability function to compare the redshift--space galaxy distribution in the Perseus--Pisces survey with artificial samples from N-body simulations of standard cold dark matter (CDM) and broken scale invariance (BSI) models. Galaxies are identified as residing in peaks of the evolved density field in such a way as to reproduce both the observed luminosity and two--point correlation functions. Using a similar approach, it was recently shown that the VPF can discriminate between CDM and a cold+hot dark matter (CHDM) model with cold/hot/baryon fractions = 0.6/0.3/0.1. Our main result is that both CDM (as expected from a previous analysis) and BSI fit observational data. The robustness of the result is checked against changing the observer's position in the simulations and the galaxy identification in the evolved density field. Therefore, while the void statistics is sensitive to the passage from CDM to CHDM (different spectrum and different nature of dark matter), it is not to the passage from CDM to BSI (different spectrum but same dark matter). On such a basis, we conjecture that the distribution of voids could be directly sensitive to the nature of dark matter, but scarcely sensitive to the shape of the transfer function.

We investigate the properties of large--scale structure predicted in a class of mixed dark matter models in which the volatile component (made of particles with high rms velocity) derives from the decay of a heavier particle. Such models based on cold+volatile dark matter (CVDM) differ from the standard mixture of CDM and massive neutrinos, usually known as CHDM, in that they involve a component which has a non--thermal phase space distribution function. As a consequence, and differently from CHDM models, the value of the redshift at which volatile particles become non relativistic, $z_{nr}$, can be varied almost independently of the volatile fraction, $\Omega_X$. We compute transfer functions for a selection of such models, having $0.1\le \Omega_X\le 0.5$ and different values of $z_{nr}$. Using linear theory and assuming a scale--free primordial spectrum, we compare such models with observational constraints on large--scale galaxy clustering and bulk flows, as well as on the abundance of galaxy clusters and high--redshift damped Ly$\alpha$ systems. We find that these constraints enable us to discriminate between different $\Omega_X$ and $z_{nr}$; within the range of the models inspected, those which can be most easily accommodated by the data correspond to the parameter choice $\Omega_X\simeq 0.2$ and $z_{nr} \simeq 2\times 10^4\Omega_X$.

We discuss the phase--space distribution of $\mu$ neutrinos if $\tau$ neutrinos are unstable and decay into $\nu_\mu + scalar$. If this scalar is a familon or a Majoron, in the generic case the $\nu_\mu$ background is NOT the straightforward overlap of neutrinos of thermal and decay origins. A delay in $\nu_\tau$ decay, due to the Pauli exclusion principle, can modify it in a significant way. We provide the equations to calculate the $\nu_\mu$ distribution and show that, in some cases, there exists a good approximate solution to them. However, even when such solution is not admitted, the equations can be numerically solved following a precise pattern. We give such a solution for a number of typical cases. If $\nu_\mu$ has a mass $\sim 2$ eV and the see--saw argument holds, $\nu_\tau$ must be unstable and the decay into $\nu_\mu + scalar$ is a reasonable possibility. The picture leads to a delayed equivalence redshift, which could allow to reconcile COBE data with a bias parameter $b\ge 1$.

We calculate the linear evolution for a class of mixed dark matter models, where the hot component derives from the decay of a heavier particle (CntHDM models). These models differ from ordinary mixed models based on massive neutrinos (CHDM models), for which the hot component has a phase space distribution which derives from a thermal one. In CntHDM models the density of the hot component and the derelativisation redshift of its quanta are indipendent parameters. In this work we provide the spectra for a number of CntHDM models, and compare them with CDM and CHDM spectra. If PQ symmetry and SUSY simultaneously hold, the lightest standard neutralino can be expected to decay into axino and photon. We briefly summarise the features of this particle model which gives rise to the cosmological framework discussed here, as a fairly generic consequence. Other frameworks which lead to similar models are also briefly discussed.

We use the void probability statistics to study the redshift-space galaxy distribution as described by a volume-limited subsample of the Perseus-Pisces survey. We compare the results with the same analysis realized on artificial samples, extracted from high-resolution N-body simulations by reproducing the observational biases of the real data set. Simulations are run for the Cold+HotDM model (CHDM) and for unbiased and biased (b=1.5) CDM models in a 50 Mpc/h box. We identify galaxies as residing in peaks of the evolved density field. We fragment overmerged structures into individual galaxies so as to reproduce both the correct luminosity function (after assuming M/ L values for the resulting galaxy groups) and the two-point correlation function. Our main result is that a void-probability function (VPF) from the standard CHDM model with fractions 60% cold, 30% hot, 10% barions, exceeds the observational VPF with a high confidence level. CDM models produce smaller VPF independent of the biasing parameter. We verify the robustness of this result against changing the observer position in the simulations and the galaxy identification in the evolved density field.

We discuss the effects of redshift--space distortions on the estimate of high order galaxy correlation functions. We use both the Perseus--Pisces redshift survey and the results of high--resolution N--body simulations to explore the consequences of working in redshift space on the detection of deviations from hierarchical clustering. Both for real and simulated data, significant deviations from hierarchical clustering seem to be present in real space. Their behaviour is coherent with that expected from an initially biased galaxy field, once the displacement from the sites where galaxies formed, due to the nonlinear gravitational evolution, is taken into account. The passage to redshift space has the net effect of filtering out the higher powers required to fit the distribution in real space. This magnifies the distortions operated by nonlinear evolution on the initial distribution, erasing the residual narrow scale range (2--5~$h^{-1}$Mpc) where deviations from hierarchical clustering are still detectable in real space. We conclude that such deviations can hardly be estimated using data in redshift space.

We use high resolution PM N-body simulations to follow the development of non-linear clustering in a flat Universe, dominated by Cold + Hot Dark Matter (CHDM) with 60% of CDM, 30% of HDM and 10% of baryons; a simulation box of 100 Mpc a side ($h=0.5$) is used. We analyze two CHDM simulations with $b =1.5$ (COBE normalization). We also compare them with CDM simulations with $b=1.5$ and $b=1$. We evaluate high-order correlation functions and the void-probability-function (VPF). Correlation functions are obtained both from counts in cells and counts of neighbors. The analysis is made for DM particles and for galaxies, identified as massive halos in the evolved density field. We also check the effects of dynamical evolution and redshift space distortions. We find that clustering of DM particles exhibits deviations from the hierarchical scaling, which decrease somewhat in redshift space. Galaxies follow hierarchical scaling far more closely, with coefficients S_3=2.5 and S_4=7.5, in general agreement with observational results, just marginally affected by redshift distortions and on the choices for the initial spectra. The hierarchical scaling of galaxy clustering is confirmed by the VPF analysis and in all the cases considered a good agreement with observational results is obtained.