(Abridged) The direct precise measurements of spectra of cosmic rays (CR) species in the wide energy range form a basis for propagation models, for interpretation of gamma-ray and microwave observations, and for disentangling possible signatures of new phenomena. Composition and spectra of CR species are vital for studies of galactic nucleosynthesis and high-energy processes in variety of environments and on different scales. Since the beginning of its operation, AMS-02 has delivered outstanding quality measurements of the spectra of CR antiprotons, electrons, positrons, and nuclei: H-O, Ne, Mg, Si. However, the spectra of heavier nuclei are expected to be available only later in the mission, which is currently extended through the end of the International Space Station (ISS) in 2030. Meanwhile, a comparison of published AMS-02 results with data from earlier experiments, such as HEAO-3-C2, indicate that earlier data may be affected by significant undocumented systematic errors. This creates an awkward situation: the new data are still incomplete, while using the old data to make up for the lack of AMS-02 measurements could result in significant errors. In this paper we are using the existing AMS-02 data to test HEAO-3-C2 results. We show that a fraction of HEAO-3-C2 data match available AMS-02 measurements quite well and can be used together with Voyager 1 and ACE-CRIS data to make predictions for the local interstellar spectra (LIS) of nuclei that are not yet released by the AMS-02 collaboration. We are also updating our already published LIS to provide a complete set of LIS from H-Ni in the energy range from 1 MeV/nucleon to ~100-500 TeV/nucleon thus covering 8-9 orders of magnitude in energy. Our calculations employ the GalProp-HelMod framework that is proved to be a reliable tool in deriving the LIS of CR antiprotons, electrons, and nuclei from H to O.

Local interstellar spectra (LIS) of secondary cosmic ray (CR) nuclei, lithium, beryllium, boron, and partially secondary nitrogen, are derived in the rigidity range from 10 MV to ~200 TV using the most recent experimental results combined with the state-of-the-art models for CR propagation in the Galaxy and in the heliosphere. The lithium spectrum appears somewhat flatter at high energies compared to other secondary species that may imply a primary lithium component. Two propagation packages, GALPROP and HelMod, are combined to provide a single framework that is run to reproduce direct measurements of CR species at different modulation levels, and at both polarities of the solar magnetic field. An iterative maximum-likelihood method is developed that uses GALPROP-predicted LIS as input to HelMod, which provides the modulated spectra for specific time periods of the selected experiments for the model-data comparison. The proposed LIS accommodate the low-energy interstellar spectra measured by Voyager 1, HEAO-3, and ACE/CRIS as well as the high-energy observations by PAMELA, AMS-02, and earlier experiments that are made deep in the heliosphere. The interstellar and heliospheric propagation parameters derived in this study are consistent with our earlier results for propagation of CR protons, helium, carbon, oxygen, antiprotons, and electrons.

A measurement of the dipole anisotropy in galactic coordinates for different charged cosmic rays has been performed with the Alpha Magnetic Spectrometer (AMS) onboard the International Space Station (ISS). Results are presented for the first 7.5 years of data taking for protons, Helium, Carbon and Oxygen, and 6.5 years for positrons and electrons. All the species are found to be consistent with isotropy and upper limits to the dipole amplitude have been computed. In particular, for energies above 16 GeV a limit of delta < 1.9% and delta < 0.5% at the 95% C.I. is found for positrons and electrons respectively. For rigidities above 200 GV a limit of delta < 0.38%, delta < 0.36%, delta < 1.9% and delta < 1.7% is obtained for protons, Helium, Carbon and Oxygen.

HelMod is a Monte Carlo code developed to describe the transport of Galactic Cosmic Rays (GCRs) through the heliosphere from the interstellar space to the Earth. In the current HelMod version 4 the modulation process, based on Parker's equation, is applied to the propagation of GCRs in the inner and outer heliosphere, i.e., including the heliosheath. HelMod was proved to reproduce protons, nuclei and electrons cosmic rays spectra observed during solar cycles 23-24 by several detectors, for instance, PAMELA, BESS and AMS-02. In particular, the unprecedented accuracy of AMS-02 observations allowed one a better tuning of the description regarding the solar modulation mechanisms implemented in HelMod. In addition, HelMod demonstrated to be capable of reproducing the uxes observed by the Voyager probes in the inner and outer regions of heliosphere up to its border.

Local interstellar spectra (LIS) of primary cosmic ray (CR) nuclei, such as helium, oxygen, and mostly primary carbon are derived for the rigidity range from 10 MV to ~200 TV using the most recent experimental results combined with the state-of-the-art models for CR propagation in the Galaxy and in the heliosphere. Two propagation packages, GALPROP and HelMod, are combined into a single framework that is used to reproduce direct measurements of CR species at different modulation levels, and at both polarities of the solar magnetic field. The developed iterative maximum-likelihood method uses GALPROP-predicted LIS as input to HelMod, which provides the modulated spectra for specific time periods of the selected experiments for model-data comparison. The interstellar and heliospheric propagation parameters derived in this study are consistent with our prior analyses using the same methodology for propagation of CR protons, helium, antiprotons, and electrons. The resulting LIS accommodate a variety of measurements made in the local interstellar space (Voyager 1) and deep inside the heliosphere at low (ACE/CRIS, HEAO-3) and high energies (PAMELA, AMS-02).

The local interstellar spectrum (LIS) of cosmic-ray (CR) electrons for the energy range 1 MeV to 1 TeV is derived using the most recent experimental results combined with the state-of-the-art models for CR propagation in the Galaxy and in the heliosphere. Two propagation packages, GALPROP and HelMod, are combined to provide a single framework that is run to reproduce direct measurements of CR species at different modulation levels, and at both polarities of the solar magnetic field. An iterative maximum-likelihood method is developed that uses GALPROP-predicted LIS as input to HelMod, which provides the modulated spectra for specific time periods of the selected experiments for model-data comparison. The optimized HelMod parameters are then used to adjust GALPROP parameters to predict a refined LIS with the procedure repeated subject to a convergence criterion. The parameter optimization uses an extensive data set of proton spectra from 1997-2015. The proposed CR electron LIS accommodates both the low-energy interstellar spectra measured by Voyager 1 as well as the high-energy observations by PAMELA and AMS-02 that are made deep in the heliosphere; it also accounts for Ulysses counting rate features measured out of the ecliptic plane. The interstellar and heliospheric propagation parameters derived in this study agree well with our earlier results for CR protons, helium nuclei, and anti-protons propagation and LIS obtained in the same framework.

The heliospheric modulation model HelMod solve the transport-equation for the Galactic Cosmic Rays propagation through the heliosphere down to Earth. It is based on a 2-D Monte Carlo approach, that includes a general description of the symmetric and antisymmetric parts of the diffusion tensor, thus, properly treating the particle drift effects as well as convection within the solar wind and adiabatic energy loss. The model was tuned in order to fit 1) the data observed outside the ecliptic plane at several distances from the Earth and 2) the spectra observed near the Earth for both, high and low solar activity periods. Great importance was given to description of polar regions of the heliosphere. We present the flux for protons, antiprotons and helium nuclei computed along the solar cycle 23-24 in comparison with experimental observations and prediction for the full solar cycle 24.

Local interstellar spectra (LIS) for protons, helium and antiprotons are built using the most recent experimental results combined with the state-of-the-art models for propagation in the Galaxy and heliosphere. Two propagation packages, GALPROP and HelMod, are combined to provide a single framework that is run to reproduce direct measurements of cosmic ray (CR) species at different modulation levels and at both polarities of the solar magnetic field. To do so in a self-consistent way, an iterative procedure was developed, where the GALPROP LIS output is fed into HelMod that provides modulated spectra for specific time periods of selected experiments to compare with the data; the HelMod parameters optimization is performed at this stage and looped back to adjust the LIS using the new GALPROP run. The parameters were tuned with the maximum likelihood procedure using an extensive data set of proton spectra from 1997-2015. The proposed LIS accommodate both the low energy interstellar CR spectra measured by Voyager 1 and the high energy observations by BESS, Pamela, AMS-01, and AMS-02 made from the balloons and near-Earth payloads; it also accounts for Ulysses counting rate features measured out of the ecliptic plane. The found solution is in a good agreement with proton, helium, and antiproton data by AMS-02, BESS, and PAMELA in the whole energy range.

The heliospheric modulation model \helmod{} is a two dimensional treatment dealing with the helio-colatitude and radial distance from Sun and is employed to solve the transport-equation for the GCR propagation through the heliosphere down to Earth. This work presents the current version 3 of the \helmod{} model and reviews how main processes involved in GCR propagation were implemented.~The treatment includes the so-called particle drift effects --~e.g., those resulting, for instance, from the extension of the neutral current sheet inside the heliosphere and from the curvature and gradient of the IMF --, which affect the transport of particles entering the solar cavity as a function of their charge sign.~The \helmod{} model is capable to provide modulated spectra which well agree within the experimental errors with those measured by AMS-01, BESS, PAMELA and AMS-02 during the solar cycles 23 and 24.~Furthermore, the counting rate measured by Ulysses at +/- 80degree of solar latitude and 1 to 5\,AU was also found in agreement with that expected by \helmod{} code version 3.

Propagation of cosmic rays (CRs) from their sources to the observer is described mainly as plain diffusion at high energies, while at lower energies there are other physical processes involved, both in the interstellar space and in the heliosphere. The latter was a subject of considerable uncertainty until recently. New data obtained by several CR missions can be used to and the local interstellar spectra (LIS) of CR species that would significantly reduce the uncertainties associated with the heliospheric propagation. In this paper we present the LIS of CR protons and helium outside the heliospheric boundary. The proposed LIS are tuned to accommodate both, the low energy CR spectra measured by Voyager 1, and the high energy observations publicly released by BESS, Pamela, AMS-01 and AMS-02. The proton and helium LIS are derived by combining CR propagation in the Galaxy, as described by GALPROP, with the heliospheric modulation computed using the HelMod Monte Carlo Tool. The proposed LIS are tuned to reproduce the modulated spectra for both, high and low, levels of solar activity.

The search for cosmic ray anisotropy has been performed using particles collected by the Alpha Magnetic Spectrometer after the first five years of operation. The positron to electron ratio is consistent with isotropy at all energies and angular scales. The limit on the dipole anisotropy at the 95% confidence level for energies above 16GeV is $\delta<0.020$ for positrons and $\delta<0.006$ for electrons. In the case of high rigidity protons the upper limit in the rigidity range from 80 to 1800GV is $\delta<0.003$. No indication of seasonal excess is observed for all particle species within the present statistics.

HelMod is a code evaluating the transport of Galactic cosmic rays through the inner heliosphere down to Earth. It is based on a 2-D Monte Carlo approach and includes a general description of the symmetric and antisymmetric parts of the diffusion tensor, thus, properly treating the particle drift effects. The model has been tuned in order to fit the data observed outside the ecliptic plane at several distances from the Earth and the spectra observed near the Earth for both, high and low solar activity levels. A stand-alone python module, fully compatible with GalProp, was developed for a comprehensive calculation of solar modulation effects, resulting in a newly suggested set of local interstellar spectra.

Euclid is a European Space Agency medium class mission selected for launch in 2020 within the Cosmic Vision 2015 2025 program. The main goal of Euclid is to understand the origin of the accelerated expansion of the universe. Euclid will explore the expansion history of the universe and the evolution of cosmic structures by measuring shapes and redshifts of galaxies as well as the distribution of clusters of galaxies over a large fraction of the sky. Although the main driver for Euclid is the nature of dark energy, Euclid science covers a vast range of topics, from cosmology to galaxy evolution to planetary research. In this review we focus on cosmology and fundamental physics, with a strong emphasis on science beyond the current standard models. We discuss five broad topics: dark energy and modified gravity, dark matter, initial conditions, basic assumptions and questions of methodology in the data analysis. This review has been planned and carried out within Euclid's Theory Working Group and is meant to provide a guide to the scientific themes that will underlie the activity of the group during the preparation of the Euclid mission.

Euclid is a European Space Agency medium class mission selected for launch in 2019 within the Cosmic Vision 2015-2025 programme. The main goal of Euclid is to understand the origin of the accelerated expansion of the Universe. Euclid will explore the expansion history of the Universe and the evolution of cosmic structures by measuring shapes and redshifts of galaxies as well as the distribution of clusters of galaxies over a large fraction of the sky. Although the main driver for Euclid is the nature of dark energy, Euclid science covers a vast range of topics, from cosmology to galaxy evolution to planetary research. In this review we focus on cosmology and fundamental physics, with a strong emphasis on science beyond the current standard models. We discuss five broad topics: dark energy and modified gravity, dark matter, initial conditions, basic assumptions and questions of methodology in the data analysis. This review has been planned and carried out within Euclid's Theory Working Group and is meant to provide a guide to the scientific themes that will underlie the activity of the group during the preparation of the Euclid mission.

The cosmic rays modulation inside the heliosphere is well described by a transport equation introduced by Parker in 1965. To solve this equation several approaches were followed in the past. Recently the Monte Carlo approach becomes widely used in force of his advantages with respect to other numerical methods. In the Monte Carlo approach, the transport equation is associated to a fully equivalent set of Stochastic Differential Equations. This set is used to describe the stochastic path of a quasi-particle from a source, e.g., the interstellar medium, to a specific target, e.g., a detector at Earth. In this work, we present both the Forward-in-Time and Backward-in-Time Monte Carlo solutions. We present an implementation of both algorithms in the framework of HelMod Code showing that the difference between the two approach is below 5\% that can be quoted as the systematic uncertain of the Method itself.

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.

We aim to constrain mass-varying neutrino models using large scale structure observations and produce forecast for the Euclid survey. We investigate two models with different scalar field potential and both positive and negative coupling parameters \beta. These parameters correspond to growing or decreasing neutrino mass, respectively. We explore couplings up to |\beta|<5. In the case of the exponential potential, we find an upper limit on $\Omega_\nu h^2$<0.004 at 2-$\sigma$ level. In the case of the inverse power law potential the null coupling can be excluded with more than 2-\sigma significance; the limits on the coupling are \beta>3 for the growing neutrino mass and \beta<-1.5 for the decreasing mass case. This is a clear sign for a preference of higher couplings. When including a prior on the present neutrino mass the upper limit on the coupling becomes |\beta|<3 at 2-$\sigma$ level for the exponential potential. Finally, we present a Fisher forecast using the tomographic weak lensing from an Euclid-like experiment and we also consider the combination with the cosmic microwave background (CMB) temperature and polarisation spectra from a Planck-like mission. If considered alone, lensing data is more efficient in constraining $\Omega_\nu$ with respect to CMB data alone. There is, however, a strong degeneracy in the \beta-$\Omega_\nu h^2$ plane. When the two data sets are combined, the latter degeneracy remains, but the errors are reduced by a factor ~2 for both parameters.

Energetic (suprathermal) solar particles, accelerated in the interplanetary medium, contribute to the solar wind pressure, in particular during high solar activity periods. We estimated the effect of the increase of solar wind pressure due to suprathermal particles on magnetospheric transmissivity of galactic cosmic rays in the case of one recent solar event.

We used a backtracing code to reconstruct particle trajectory inside the Earth Magnetosphere during the last solar active period (2011 and 2012) when very high Solar Wind pressure values were measured. We compared our results on AMS-02 proton and electron data with 2 different External Field models, namely Tsyganenko 1996 (T96) and 2005 (T05), both for quiet (defined as the periods when the solar wind pressure is below the average value, set at 2nPa) and active periods. Although T05 has been specifically designed for storm events, at high energy the particle trajectory is similar for the two models. For instance at rigidities larger than 50 GV, the RMS of angular difference between reconstructed asymptotic direction outside the Magnetosphere is of the order of the millirad, while it increases at intermediate energies. We also confirmed, as a function of the pointing direction, the well known East-West effect on the trajectory of primary particles and on the access solid angle on the AMS detector.

We implemented a website to deal with main effects on Cosmic Ray access to the Earth, i.e. the Solar Modulation and the Geomagnetic Field effect. In helmod.org the end user can easily access a web interface to results catalog of the HelMod Monte Carlo Code. This Model uses a Monte Carlo Approach to solves the Parker Transport Equation, obtaining a modulated proton flux for a period (monthly average) between January 1990 and december 2007. geomagsphere.org is instead based on GeoMag Backtracing Code, that solves the Lorentz equation with a Runge-Kutta method of 6th order, and, reversing charge sign and velocity, reconstruct particle trajectories in the Earth Magnetosphere back in time. We use last models of internal (IGRF-11) and external (Tsyganenko 1996 -T96- and 2005 -T05-) field components valid up to 2015. Particles are backtraced to the outer (magnetopause) or inner boundary to separate Primary (allowed trajectory) from Secondary (forbidden) Cosmic Rays. This code has been used both for reproducing known effects as East-West effect and rigidity cutoff calculations. In geomagsphere.org the user can choose the external field model from Tsyganenko (T96 or T05) and obtain for a fixed position and date from 1st Jan. 1968 (T96) and 1st Jan. 1995 (T05) respectively till 31$^{st}$ Dec 2012, the vertical rigidity cutoff estimation obtained with the backtracing technique with a rigidity step of 0.1 GV. For a more precise calculation (0.01 GV), requiring more CPU time, results are sent to the user by email (mail model)

The Cosmic Rays propagation was studied in details using the HelMod-2D Monte Carlo code, that includes a general description of the diffusion tensor, and polar magnetic-field. The Numerical Approach used in this work is based on a set of Stochastic Differential Equations fully equivalent to the well know Parker Equation for the transport of Cosmic Rays. In our approach the Diffusion tensor in the frame of the magnetic field turbolence does not depends explicitly by Solar Latitude but varies with time using a diffusion parameter obtained by Neutron Monitors. The parameters of the Model were tuned using data during the solar Cycle 23 and Ulysses latitudinal Fast Scan in 1995. The actual parametrization is able to well reproduce the observed latitudinal gradient of protons and the southward shift of the minimum of latitudinal intensity. The description of the model is also available online at website www.helmod.org. The model was then applied on Pamela/Ulysses proton intensity from 2006 up to 2009. The model during this 4-year continous period agree well with both PAMELA (at 1 AU) and Ulysses data (at various solar distance and solar latitude). The agreement improves when considering the ratio between this data. Studies done also with particles with different charge (e.g. electrons) allow us to explain the presence (or not) of protons and electrons latitudinal gradients observed by Ulysses during the Latitudinal Fast Scan in 1995 and 2007.

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.

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.

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 use WMAP5 and other cosmological data to constrain model parameters in quintessence cosmologies, focusing also on their shift when we allow for non-vanishing neutrino masses. The Ratra-Peebles (RP) and SUGRA potentials are used here, as examples of slowly or fastly varying state parameter w(a). Both potentials depend on an energy scale \Lambda. Here we confirm the results of previous analysis with WMAP3 data on the upper limits on \Lambda, which turn out to be rather small (down to ~10^{-9} in RP cosmologies and ~10^{-5} for SUGRA). Our constraints on \Lambda are not heavily affected by the inclusion of neutrino mass as a free parameter. On the contrary, when the neutrino mass degree of freedom is opened, significant shifts in the best-fit values of other parameters occur.

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.

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.

Upcoming Weak Lensing (WL) surveys can be used to constrain Dark Energy (DE) properties, namely if tomographic techniques are used to improve their sensitivity. In this work, we use a Fisher matrix technique to compare the power of CMB anisotropy and polarization data with tomographic WL data, in constraining DE parameters. Adding WL data to available CMB data improves the detection of all cosmological parameters, but the impact is really strong when DE--DM coupling is considered, as WL tomography can then succeed to reduce the errors on some parameters by factors >10.