The Large Scale Polarization Explorer (LSPE) is a cosmology program for the measurement of large scale curl-like features (B-modes) in the polarization of the Cosmic Microwave Background. Its goal is to constrain the background of inflationary gravity waves traveling through the universe at the time of matter-radiation decoupling. The two instruments of LSPE are meant to synergically operate by covering a large portion of the northern microwave sky. LSPE/STRIP is a coherent array of receivers planned to be operated from the Teide Observatory in Tenerife, for the control and characterization of the low-frequency polarized signals of galactic origin; LSPE/SWIPE is a balloon-borne bolometric polarimeter based on 330 large throughput multi-moded detectors, designed to measure the CMB polarization at 150 GHz and to monitor the polarized emission by galactic dust above 200 GHz. The combined performance and the expected level of systematics mitigation will allow LSPE to constrain primordial B-modes down to a tensor/scalar ratio of $10^{-2}$. We here report the status of the STRIP pre-commissioning phase and the progress in the characterization of the key subsystems of the SWIPE payload (namely the cryogenic polarization modulation unit and the multi-moded TES pixels) prior to receiver integration.

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.

In this paper we discuss the latest developments of the STRIP instrument of the "Large Scale Polarization Explorer" (LSPE) experiment. LSPE is a novel project that combines ground-based (STRIP) and balloon-borne (SWIPE) polarization measurements of the microwave sky on large angular scales to attempt a detection of the "B-modes" of the Cosmic Microwave Background polarization. STRIP will observe approximately 25% of the Northern sky from the "Observatorio del Teide" in Tenerife, using an array of forty-nine coherent polarimeters at 43 GHz, coupled to a 1.5 m fully rotating crossed-Dragone telescope. A second frequency channel with six-elements at 95 GHz will be exploited as an atmospheric monitor. At present, most of the hardware of the STRIP instrument has been developed and tested at sub-system level. System-level characterization, starting in July 2018, will lead STRIP to be shipped and installed at the observation site within the end of the year. The on-site verification and calibration of the whole instrument will prepare STRIP for a 2-years campaign for the observation of the CMB polarization.

We present an optimization scheme for the focal plane of a multi-frequency CMB B-modes experiment with a fixed number of detectors and apply it to the specific case of LSPE experiment. Optimal focal planes are identified on the ground of different figures of merit defined in terms of the forecasted uncertainty $\sigma_r$ on the tensor--to-scalar ratio $r$ and the expected map variance from foreground and instrumental noise residuals. We then perform a forecast analysis in order to assess the precision achievable in B-modes measurements.

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.

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.

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.

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.

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.

Halo masses and concentrations have been studied extensively, by means of N-body simulations as well as observationally, during the last decade. Nevertheless, the exact form of the mass-concentration relation is still widely debated. One of the most promising method to estimate masses and concentrations relies on gravitational lensing from massive halos. Here we investigate the impact of the mass-concentration relation on halo peak abundance in weak lensing shear maps relying on the aperture mass method for peak detections. After providing a prescription to take into account the concentration dispersion (always neglected in previous works) in peak number counts predictions, we assess their power to constrain the mass-concentration relation by means of Fisher matrix technique. We find that, when combined with different cosmological probes, peak statistics information from near-future weak lensing surveys provides an interesting and complementary alternative method to lessen the long standing controversy about the mass-concentration relation.

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.

Circular polarization of the Cosmic Microwave Background (CMB) offers the possibility of detecting rotations of the universe and magnetic fields in the primeval universe or in distant clusters of galaxies. We used the Milano Polarimeter (MIPOL) installed at the Testa Grigia Observatory, on the italian Alps, to improve the existing upper limits to the CMB circular polarization at large angular scales. We obtain 95% confidence level upper limits to the degree of the CMB circular polarization ranging between 5.0x10^{-4} and 0.7x10^{-4} at angular scales between 8 and 24 deg, improving by one order of magnitude preexisting upper limits at large angular scales. Our results are still far from the nK region where today expectations place the amplitude of the V Stokes parameter used to characterize circular polarization of the CMB but improve the preexisting limit at similar angular scales. Our observations offered also the opportunity of characterizing the atmospheric emission at 33 GHz at the Testa Grigia Observatory.

Weak gravitational lensing by galaxy clusters on faint higher redshift galaxies has been traditionally used to study the cluster mass distribution and as a tool to identify clusters as peaks in the shear maps. However, it becomes soon clear that peaks statistics can also be used as a way to constrain the underlying cosmological model due to its dependence on both the cosmic expansion rate and the growth rate of structures. This feature makes peak statistics particularly interesting from the point of view of discriminating between General Relativity and modified gravity. Here we consider a general class of $f(R)$ theories and compute the observable mass function based on the aperture mass statistics. We complement our theoretical analysis with a Fisher matrix forecast of the constraints that an Euclid\,-\,like survey can impose on the $f(R)$ model parameters. We show that peak statistics alone can in principle discriminate between General Relativity and $f(R)$ models and strongly constrain the $f(R)$ parameters that are sensitive to the non-linear growth of structure. However, further analysis is needed in order to include possible selection function in the peaks redshift determination.

We discuss the design and expected performance of STRIP (STRatospheric Italian Polarimeter), an array of coherent receivers designed to fly on board the LSPE (Large Scale Polarization Explorer) balloon experiment. The STRIP focal plane array comprises 49 elements in Q band and 7 elements in W-band using cryogenic HEMT low noise amplifiers and high performance waveguide components. In operation, the array will be cooled to 20 K and placed in the focal plane of a $\sim 0.6$ meter telescope providing an angular resolution of $\sim1.5$ degrees. The LSPE experiment aims at large scale, high sensitivity measurements of CMB polarization, with multi-frequency deep measurements to optimize component separation. The STRIP Q-band channel is crucial to accurately measure and remove the synchrotron polarized component, while the W-band channel, together with a bolometric channel at the same frequency, provides a crucial cross-check for systematic effects.

We provide an exhaustive analysis of the Integrated Sach-Wolfe effect (ISW) in the context of coupled Dark Energy cosmologies where a component of massive neutrinos is also present. We focus on the effects of both the coupling between Dark Matter and Dark Energy and of the neutrino mass on the cross-correlation between galaxy/quasar distributions and ISW effect. We provide a simple expression to appropriately rescale the galaxy bias when comparing different cosmologies. Theoretical predictions of the cross-correlation function are then compared with observational data. We find that, while it is not possible to distinguish among the models at low redshifts, discrepancies between coupled models and $\Lambda$CDM increase with $z$. In spite of this, current data alone seems not able to distinguish between coupled models and $\Lambda$CDM. However, we show that upcoming galaxy surveys will permit tomographic analysis which allow to better discriminate among the models. We discuss the effects on cross-correlation measurements of ignoring galaxy bias evolution, b(z), and magnification bias correction and provide fitting formulae for b(z) for the cosmologies considered. We compare three different tomographic schemes and investigate how the expected signal to noise ratio, snr, of the ISW-LSS cross-correlation changes when increasing the number of tomographic bins. The dependence of snr on the area of the survey and the survey shot noise is also discussed.

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.

We investigate the clustering properties of dynamical Dark Energy even in association of a possible coupling between Dark Energy and Dark Matter. We find that within matter inhomogeneities, Dark Energy migth form voids as well as overdensity depending on how its background energy density evolves. Consequently and contrarily to what expected, Dark Energy fluctuations are found to be slightly suppressed if a coupling with Dark Matter is permitted. When considering density contrasts and scales typical of superclusters, voids and supervoids, perturbations amplitudes range from $|\delta_\phi|\sim {\cal O} (10^{-6})$ to $|\delta_\phi|\sim {\cal O} (10^{-4})$ indicating an almost homogeneous Dark Energy component.

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.

The properties of dynamical Dark Energy (DE) and, in particular, the possibility that it can form or contribute to stable inhomogeneities, have been widely debated in recent literature, also in association to a possible coupling between DE and Dark Matter (DM). In order to clarify this issue, in this paper we present a general framework for the study of the nonlinear phases of structure formation, showing the equivalence between two possible descriptions of DE: a scalar field \phi self-interacting through a potential V(\phi) and a perfect fluid with an assigned negative equation of state w(a). This enables us to show that, in the presence of coupling, the mass of DE quanta may increase where large DM condensations are present, so that also DE may partake to the clustering process.

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.

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.

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.

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.

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..

The growth and virialization of spherical top-hat fluctuations, in coupled Dark Energy models, causes segregation between Dark Matter (DM) and baryons, as the gravitational infall into the potential well proceeds more slowly for the baryons than for DM. As a consequence, after attaining their turn-around and before full virialization, halos have outer layers rich of baryons. Accordingly, a natural ambiguity exists on the definition of the virial density contrast. In fact, when the outer baryon layers infall onto the DM-richer core, they carry with them DM materials outside the original fluctuation; hence, no time exists when all materials originally belonging to the fluctuation -and only them- have virialized. Baryon-DM segregation can have various astrophysical consequences on different length-scales. The smallest halos may loose up to 50% of the original baryonic contents and become hardly visible. Subhalos in cluster-size halos may loose much baryonic materials, which could then be observed as intra-cluster light. Isolated halos, in general, can be expected to have a baryon component richer than the cosmological proportions, due to the cosmic enrichement of baryons lost in small halo encounters.

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.

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.

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.

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 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.