We compute a new generation of standard solar models (SSMs) that includes recent updates on some important nuclear reaction rates and a more consistent treatment of the equation of state. Models also include a novel and flexible treatment of opacity uncertainties based on opacity kernels, required in the light of recent theoretical and experimental works on radiative opacity. Two large sets of SSMs, each based on a different canonical set of solar abundances with high and low metallicity (Z), are computed to determine model uncertainties and correlations among different observables. We present detailed comparisons of high- and low-Z models against different ensembles of solar observables including solar neutrinos, surface helium abundance, depth of convective envelope and sound speed profile. A global comparison, including all observables, yields a p-value of 2.7$\sigma$ for the high-Z model and 4.7$\sigma$ for the low-Z one. When the sound-speed differences in the narrow region of $0.65 < r/R_{sun} < 0.70$ are excluded from the analysis, results are 0.9$\sigma$ and 3.0$\sigma$ for high- and low-Z models respectively. These results show that: high-Z models agree well with solar data but have a systematic problem right below the bottom of the convective envelope linked to steepness of molecular weight and temperature gradients, and that low-Z models lead to a much more general disagreement with solar data. We also show that, while simple parametrizations of opacity uncertainties can strongly alleviate the solar abundance problem, they are insufficient to substantially improve the agreement of SSMs with helioseismic data beyond that obtained for high-Z models due to the intrinsic correlations of theoretical predictions. All the detailed information of the new solar models, can be found online: http://www.ice.cat/personal/aldos/Solar_Data.html

We investigate the non-linear evolution of the relic cosmic neutrino background by running large box-size, high resolution N-body simulations. Our set of simulations explore the properties of neutrinos in a reference $\Lambda$CDM model with total neutrino masses between 0.05-0.60 eV in cold dark matter haloes of mass $10^{11}-10^{15}$ $h^{-1}$M$_{\odot}$, over a redshift range $z=0-2$. We compute the halo mass function and show that it is reasonably well fitted by the Sheth-Tormen formula. More importantly, we focus on the CDM and neutrino properties of the density and peculiar velocity fields in the cosmological volume, inside and in the outskirts of virialized haloes. The dynamical state of the neutrino particles depends strongly on their momentum: whereas neutrinos in the low velocity tail behave similarly to CDM particles, neutrinos in the high velocity tail are not affected by the clustering of the underlying CDM component. We find that the neutrino (linear) unperturbed momentum distribution is modified and mass and redshift dependent deviations from the expected Fermi-Dirac distribution are in place both in the cosmological volume and inside haloes. The neutrino density profiles around virialized haloes have been carefully investigated and a simple fitting formula is provided. The neutrino profile, unlike the cold dark matter one, is found to be cored with core size and central density that depend on the neutrino mass, redshift and mass of the halo, for halos of masses larger than $\sim 10^{13.5}h^{-1}$M$_{\odot}$. For lower masses the neutrino profile is best fitted by a simple power-law relation in the range probed by the simulations. Our findings are particularly important in view of upcoming large-scale structure surveys, like Euclid, that are expected to probe the non-linear regime at the percent level with lensing and clustering observations.

We study whether non-linear gravitational effects of relic neutrinos on the development of clustering and large-scale structure may be observable by weak gravitational lensing. We compute the density profile of relic massive neutrinos in a spherical model of a cluster of galaxies, for several neutrino mass schemes and cluster masses. Relic neutrinos add a small perturbation to the mass profile, making it more extended in the outer parts. In principle, this non-linear neutrino perturbation is detectable in an all-sky weak lensing survey such as EUCLID by averaging the shear profile of a large fraction of the visible massive clusters in the universe, or from its signature in the general weak lensing power spectrum or its cross-spectrum with galaxies. However, correctly modeling the distribution of mass in baryons and cold dark matter and suppressing any systematic errors to the accuracy required for detecting this neutrino perturbation is severely challenging.

We provide a summary of the current knowledge, theoretical and experimental, of solar neutrino fluxes and of the masses and mixing angles that characterize solar neutrino oscillations. We also summarize the principal reasons for doing new solar neutrino experiments and what we think may be learned from the future measurements.

We present "Before and After" global oscillation solutions, as well as predicted "Before and After" values and ranges for ten future solar neutrino observables (for BOREXINO, KamLAND, SNO, and a generic p-p neutrino detector). We have performed global analyses using the full SNO day-night energy spectrum and, alternatively, just the SNO CC and NC rates and the day-night asymmetry.The ``Before'' case includes all solar neutrino data (and some theoretical improvements) available prior to April 20, 2002 and the ``After'' case includes, in addition, the new SNO data on the CC, NC, and day-night asymmetry. The LMA solution is the only currently allowed MSW oscillation solution at ~99% CL. The LOW solution is allowed only at more than 2.5 sigma, SMA at 3.7 sigma or 4.7 sigma (depending upon analysis strategy), and pure sterile oscillations at 4.7sigma. Small mixing angles are ``out''(pure sterile is ``way out''); MSW with large mixing angles is definitely ``in''. Vacuum oscillations are allowed at 3 sigma, but not a 2 sigma. Precise maximal mixing is excluded at 3.2 sigma for MSW solutions and at 2.8 sigma for vacuum solutions. Most of the predicted values for future observables for the BOREXINO, KamLAND, and future SNO measurements are changed only by minor amounts by the inclusion of the recent SNO data. In order to test the robustness of the allowed neutrino oscillation regions and the predictions for future observables, we have carried out calculations using a variety of strategies to analyze the SNO and other solar neutrino data.

We discuss the potential of SNO to constraint the four-neutrino mixing schemes favoured by the results of all neutrino oscillations experiments. These schemes allow simultaneous transitions of solar $\nu_e's$ into active $\nu_\mu$'s, $\nu_\tau$'s and sterile $\nu_s$ controlled by the additional parameter $\cos^2(\vartheta_{23}) \cos^2(\vartheta_{24})$ and they contain as limiting cases the pure $\nu_e$-active and $\nu_e$-sterile neutrino oscillations. We first obtain the solutions allowed by the existing data in the framework of the BP00 standard solar model and quantify the corresponding predictions for the CC and the NC/CC event ratios at SNO for the different allowed regions as a function of the active-sterile admixture. Our results show that some information on the value of $\cos^2(\vartheta_{23}) \cos^2(\vartheta_{24})$ can be obtained by the first SNO measurement of the CC ratio, while considerable improvement on the knowledge of this mixing will be achievable after the measurement of the NC/CC ratio.

We have performed a detailed study of the zenith angle dependence of the regeneration factor and distributions of events at SNO and SK for different solutions of the solar neutrino problem. In particular, we discuss oscillatory behaviour and the synchronization effect in the distribution for the LMA solution, the parametric peak for the LOW solution, etc.. Physical interpretation of the effects is given. We suggest a new binning of events which emphasizes distinctive features of zenith angle distributions for the different solutions. We also find the correlations between the integrated day-night asymmetry and the rates of events in different zenith angle bins. Study of these correlations strengthens the identification power of the analysis.

We study the phenomenological consequences of maximal and near-maximal mixing of the electron neutrino with other ($x$=tau and/or muon) neutrinos. We describe the deviations from maximal mixing in terms of a parameter $\epsilon\equiv1-2\sin^2\theta_{ex}$ and quantify the present experimental status for $|\epsilon|<0.3$. We find that the global analysis of solar neutrino data allows maximal mixing with confidence level better than 99% for $10^{-8}$ eV$^2\lsim\Delta m^2\lsim2\times10^{-7}$ eV$^2$. In the mass ranges $\Delta m^2\gsim 1.5\times10^{-5}$ eV$^2$ and $4\times10^{-10}$ eV$^2\lsim\Delta m^2\lsim2\times10^{-7}$ eV$^2$ the full interval $|\epsilon|<0.3$ is allowed within 4$\sigma$(99.995 % CL). We suggest ways to measure $\epsilon$ in future experiments. The observable that is most sensitive to $\epsilon$ is the rate [NC]/[CC] in combination with the Day-Night asymmetry in the SNO detector. With theoretical and statistical uncertainties, the expected accuracy after 5 years is $\Delta \epsilon\sim 0.07$. We also discuss the effects of maximal and near-maximal $\nu_e$-mixing in atmospheric neutrinos, supernova neutrinos, and neutrinoless double beta decay.

We study the effect of random matter density fluctuations in the sun on resonant neutrino conversion in matter by solving numerically the evolution equation for the neutrino system including the full effect of the random matter density fluctuations of given amplitude and correlation length. In order to establish the possible effect on the MSW solutions to the solar neutrino problem we perform a global analysis of all the existing observables including the measured total rates as well as the Super-Kamiokande measurement on the time dependence of the event rates during the day and night and the recoil electron energy spectrum. We find the effects of random noise to be larger for small mixing angles and they are mostly important for correlation lengths in the range few 100 km $\lesssim L_0\lesssim$ few 1000 km. They can be understood as due to a parametric resonance occuring when the phase acquired by the oscillating neutrino state on one fluctuation length L_0 is a multiple of 2$\pi$. We find that this resonant parametric condition is mainly achieved for low energy neutrinos such as the pp-neutrinos and therefore its effect is mostly seen on the total event rates while the other Super-Kamiokande observables are very marginally sensitive to the presence of noise due to the higher energy threshold.

We present an analysis of the MSW neutrino oscillation solutions of the solar neutrino problem in the framework of two-neutrino mixing in the enlarged parameter space $(\Delta m^2, \tan^2\theta)$ with $\theta \in (0,\frac{\pi}{2})$. Recently, it was pointed out that the allowed region of parameters from a fit to the measured total rates can extend to values $\theta \geq \frac{\pi}{4}$ (the so called ``dark side'') when higher confidence levels are allowed. The purpose of this letter is to reanalize the problem including all the solar neutrino data available, to discuss the dependence on the statistical criteria in the determination of the CL of the ``dark side'' and to extract the corresponding limits on the largest mixing allowed by the data. Our results show that when the Super-Kamiokande data on the zenith angle distribution of events and the spectrum information is included, the regions extend more into the dark side.

We present an analysis of the neutrino oscillation solutions of the solar neutrino problem in the framework of four-neutrino mixing where a sterile neutrino is added to the three standard ones. We perform a fit to the full data set corresponding to the 825-day Super-Kamiokande data sample as well as to Chlorine, GALLEX and SAGE and Kamiokande experiments. In our analysis we use all measured total event rates as well as all Super-Kamiokande data on the zenith angle dependence and the recoil electron energy spectrum. We consider both transitions via the Mikheyev-Smirnov-Wolfenstein (MSW) mechanism as well as oscillations in vacuum (just-so) and find the allowed solutions for different values of the additional mixing angles. This framework permits transitions into active or sterile neutrinos controlled by the additional parameter $\cos^2(\vartheta_{23}) \cos^2(\vartheta_{24})$ . We discuss the maximum allowed values of this additional mixing parameter for the different solutions.

We present an updated global analysis of two-flavor MSW solutions to the solar neutrino problem in terms of conversions of nu_e into active or sterile neutrinos. We perform a fit to the full data set corresponding to the 825-day Super-Kamiokande data sample as well as to Chlorine, GALLEX and SAGE experiments. We use all measured total event rates as well as Super-Kamiokande data on the zenith angle dependence, energy spectrum and seasonal variation of the events. For conversions into active neutrinos we find that, although the data on the total event rates favours the Small Mixing Angle (SMA) solution, once the full data set is included both SMA and Large Mixing Angle (LMA) solutions give an equivalent fit to the data. The best-fit points for the combined analysis are Delta m^2=3.6 10^{-5} eV^2 sin^2(2theta)=0.79 with chi^2_min=35.4/30 d.o.f and Delta m^2=5.1~ 10^{-6} eV^2 and sin^2(2theta)=5.5 10^{-3} with chi^2_min=37.4/30 d.o.f. In contrast with the earlier 504-day study of Bahcall-Krastev-Smirnov our results indicate that the LMA solution is not only allowed, but slightly preferred. We also show that seasonal effects, although small, may still reach 8.5 % at the best-fit LMA solution without conflict with the hints of a day-night variation (6% is due to the eccentricity of the Earth's orbit). For conversions into sterile neutrinos only the SMA solution is allowed with best-fit point Delta m^2=5.0 10^{-6} eV^2 and sin^2(2theta)=3. 10^{-3} and chi^2_min=40.2/30 d.o.f. We also consider departures of the Standard Solar Model of Bahcall and Pinsonneault 1998 by allowing arbitrary 8^B and hep fluxes. Best fit is obtained for ^8B/^8B_SSM=0.61 (1.37) and hep/hep_SSM=12 (38) for the SMA (LMA) solution.

MSW solutions of the solar neutrino problem predict a seasonal dependence of the zenith angle distribution of the event rates, due to the non-zero latitude at the Super-Kamiokande site. We calculate this seasonal dependence and compare it with the expectations in the no-oscillation case as well as just-so scenario, in the light of the latest Super-Kamiokande 708-day data. The seasonal dependence can be sizeable in the large mixing angle MSW solution and would be correlated with the day-night effect. This may be used to discriminate between MSW and just-so scenarios and should be taken into account in refined fits of the data.