Because the baryon-to-photon ratio eta_{10} is in some doubt, we drop nucleosynthetic constraints on eta_{10} and fit the three cosmological parameters (h, Omega_M, eta_{10}) to four observational constraints: Hubble parameter h_o = 0.70+-0.15, age of the universe t_o = 14+7-2 Gyr, cluster gas fraction f_o \equiv f_G h^{3/2} = 0.060 +- 0.006, and effective shape parameter Gamma_o = 0.255 +- 0.017. Errors quoted are 1 sigma, and we assume Gaussian statistics. We experiment with a fifth constraint Omega_o = 0.2 +- 0.1 from clusters. We set the tilt parameter n = 1 and the gas enhancement factor Upsilon = 0.9. We consider CDM models (open and Omega_M = 1) and flat LambdaCDM models. We test goodness of fit and draw confidence regions by the Delta chi^2 method. CDM models with Omega_M = 1 (SCDM models) are accepted only because the large error on h_o allows h < 0.5. Baryonic matter plays a significant role in Gamma_o when Omega_M \sim 1. Open CDM models are accepted only for Omega_M \gtrsim 0.4. The combination of the four other constraints with Omega_o = 0.2 +- 0.1 is rejected in CDM models with 98% confidence, suggesting that light may not trace mass. LambdaCDM models give similar results. In all of these models, eta_{10} \gtrsim 6 is favored strongly over eta_{10} \lesssim 2. This suggests that reports of low deuterium abundances on QSO lines of sight may be correct and that observational determinations of primordial 4He may have systematic errors. Plausible variations on n and Upsilon in our models do not change the results much. If we drop or change the crucial Gamma_o constraint, lower values of Omega_M and eta_{10} are permitted. The constraint Gamma_o = 0.15 +- 0.04, derived recently from the IRAS redshift survey, favors Omega_M \approx 0.3 and eta_{10} \approx 5 but does not exclude eta_{10} \approx 2.

Since the baryon-to-photon ratio "eta" is in some doubt at present, we ignore the constraints on eta from big bang nucleosynthesis (BBN) and fit the three key cosmological parameters (h, Omega_M, eta) to four other observational constraints: Hubble parameter, age of the universe, cluster gas (baryon) fraction, and effective shape parameter "Gamma". We consider open and flat CDM models and flat "Lambda"-CDM models, testing goodness of fit and drawing confidence regions by the Delta-chi^2 method. CDM models with Omega_M = 1 (SCDM models) are accepted only because we allow a large error on h, permitting h < 0.5. Open CDM models are accepted only for Omega_M \gsim 0.4. Lambda-CDM models give similar results. In all of these models, large eta (\gsim 6) is favored strongly over small eta, supporting reports of low deuterium abundances on some QSO lines of sight, and suggesting that observational determinations of primordial 4He may be contaminated by systematic errors. Only if we drop the crucial Gamma constraint are much lower values of Omega_M and eta permitted.

It is shown that the metal abundance in the solar interior is constrained from the current solar neutrino experiment, even if one allows neutrino conversion in the sun due to neutrino oscillation. The result shows that the metal abundance in the interior should be within the range 0.4 - 1.4 times that in the surface, supporting the idea that the sun formed by a contraction of a gas cloud with an almost homogeneous composition.

We present an updated analysis of astrophysical solutions, two-flavor MSW solutions, and vacuum oscillation solutions to the solar neutrino anomaly. The recent results of each of the five solar neutrino experiments are incorporated, including both the zenith angle (day-night) and spectral information from the Kamiokande experiment, and the preliminary Super-Kamiokande results. New theoretical developments include the use of the most recent Bahcall-Pinsonneault flux predictions (and uncertainties) and density and production profiles, the radiative corrections to the neutrino-electron scattering cross section, and new constraints on the Ga absorption cross section inferred from the gallium source experiments. From a model independent analysis, arbitrary astrophysical solutions are excluded at more that 98% C.L. even if one ignores any one of the three classes of experiment, relaxes the luminosity constraint, or allows more suppression of the 7Be than 8B flux. The data is well described by large and small mixing angle two-flavor MSW conversions, MSW conversions into a sterile neutrino with small mixing, or vacuum oscillations. We also present MSW fits for nonstandard solar models parameterized by an arbitrary solar core temperature or arbitrary 8B flux.

Constraints on big bang nucleosynthesis (BBN) and on cosmological parameters from conflicting deuterium observations in different high red-shift QSO systems are discussed. The high deuterium observations by Carswell et al., Songaila et al., and Rugers & Hogan is consistent with 4He and 7Li observations and Standard BBN (N_nu =3) and allows N_nu \leq 3.6 at 95% C.L., but is inconsistent with local observations of D and 3He in the context of conventional theories of stellar and Galactic evolution. In contrast, the low deuterium observations by Tytler, Fan & Burles and Burles & Tytler are consistent with the constraints from local Galactic observations, but require N_nu = 1.9 +/- 0.3 at 68% C.L., excluding Standard BBN at 99.9% C.L., unless the systematic uncertainties in the 4He observations have been underestimated by a large amount. The high and low primordial deuterium abundances imply, respectively, Omega_B h^2 = 0.005 - 0.01 and Omega_B h^2 = 0.02 - 0.03 at 95% C.L. When combined with the high baryon fraction inferred from x-ray observations of rich clusters, the corresponding total mass densities (for 50 \le H_0 \le 90) are Omega_M = 0.05 - 0.20 and Omega_M = 0.2 - 0.7, respectively (95\% C.L.) The range of Omega_M corresponding to high D is in conflict with dynamical constraints (Omega_M > 0.2 - 0.3) and with the shape parameter constraint (Gamma = Omega_M h = 0.25 +/- 0.05) from large scale structure formation in CDM and Lambda-CDM models.

A new evaluation of the constraint on the number of light neutrino species (N_nu) from big bang nucleosynthesis suggests a discrepancy between the predicted light element abundances and those inferred from observations, unless the inferred primordial 4He abundance has been underestimated by 0.014 +/- 0.004 (1 sigma) or less than 10% (95%C.L.) of 3He survives stellar processing. With the quoted systematic errors in the observed abundances and a conservative chemical evolution parameterization, the best fit to the combined data is N_nu = 2.1 +/- 0.3 (1 sigma) and the upper limit is N_nu < 2.6 (95% C.L.). The data are inconsistent with the Standard Model (N_nu = 3) at the 98.6% C.L.

We present new upper and lower bounds to the primordial abundances of deuterium and helium-3 based on observational data from the solar system and the interstellar medium. Independent of any model for the primordial production of the elements we find (at the 95\% C.L.): $1.5 \times 10^{-5} \le (D/H)_P \le 10.0 \times 10^{-5}$ and $(^3He/H)_P \le 2.6\times 10^{-5}$. When combined with the predictions of standard big bang nucleosynthesis, these constraints lead to a 95\% C.L. bound on the primordial abundance of deuterium: $(D/H)_{best} = (3.5^{+2.7}_{-1.8})\times 10^{-5}$. Measurements of deuterium absorption in the spectra of high redshift QSOs will directly test this prediction. The implications of this prediction for the primordial abundances of helium-4 and lithium-7 are discussed, as well as those for the universal density of baryons.

We argue that, prior to the recent GALLEX $^{51}$Cr source experiment, the excited state contributions to the $^{71}$Ga capture cross section for $^{51}$Cr and $^7$Be neutrinos were poorly constrained, despite forward-angle (p,n) measurements. We describe the origin of the uncertainties and estimate their extent. We explore the implications of the source experiment for solar neutrino capture in light of these uncertainties. A reanalysis of the $^7$Be and $^8$B flux constraints and MSW solutions of the solar neutrino puzzle is presented.

Besides the opportunity for discovering new neutrino physics, solar neutrino measurements provide a sensitive probe of the solar interior, and thus a rigorous test of solar model predictions. We present model independent determinations of the neutrino spectrum by using relevant flux components as free parameters subject only to the luminosity constraint. (1) Without the Mikheyev-Smirnov-Wolfenstein (MSW) effect, the best fit for the combined data is poor. Furthermore, the data indicate a severe suppression of the $^7$Be flux relative to the $^8$B, contradicting both standard and nonstandard solar models in general; the $pp$ flux takes its maximum value allowed by the luminosity constraint. This pathology consistently appears even if we ignore any one of the three data. (2) In the presence of the two-flavor MSW effect, the current constraint on the initial $^8$B flux is weak, but consistent with the SSM and sufficient to exclude nonstandard models with small $^8$B fluxes. No meaningful constraint is obtained for the other fluxes. In the future, even allowing MSW, the $^8$B and $^7$Be fluxes can be determined at the $\pm$(15 -- 20)\% level, making competing solar models distinguishable. We emphasize that the neutral current sensitivity for $^7$Be neutrinos in BOREXINO, HELLAZ, and HERON is essential for determining the initial fluxes. The constraints on the MSW parameters in the model independent analysis are also discussed.

The present status of the solar neutrino problem is reviewed. The strongest motivation for new neutrino physics comes from the complete phenomenological failure of astrophysical solutions: (1) The standard solar model is excluded by each of the solar neutrino experiments. (2) The combined results of Homestake and Kamiokande are incompatible with any astrophysical solution. (3) Even if the Homestake results were ignored entirely, there is no realistic solar model so far that describes the Kamiokande and gallium data simultaneously within their experimental uncertainties. On the other hand, the MSW solutions provide a complete description of the data, strongly hinting at the existence of mass and mixing of neutrinos. The MSW predictions for the future experiments are given in detail. It is stressed that solar model uncertainties are significant in determining the oscillation parameters and in making predictions for the future experiments. Especially, a $^8$B flux (or $S_{17}$) smaller by 20\% than that of the Bahcall-Pinsonneault model can make the interpretation of the charged to neutral current ratio measurement in SNO ambiguous.

Vacuum oscillations are considered for the combined solar neutrino observations, including the Kamiokande II spectrum data and incorporating theoretical uncertainties and their correlations. Despite the conceptual difficulty of the fine tuning between the neutrino parameters and the Sun-Earth distance, 2-flavor vacuum oscillations provide phenomenologically acceptable solutions. There are allowed regions at 99\% C.L. for $\Delta m^2 = (0.45 - 1.2) \times 10^{-10} \; \mbox{eV}\,^2$ and $\sin^22\theta = 0.6 - 1$; the best fit solution is $\chi^2 / \mbox{d.f.} = 19.2 / 16$, which is acceptable at 16\% C.L. Oscillations for sterile neutrinos are, however, excluded by the averaged data at 99.4\% C.L. The vacuum oscillation solutions predict characteristic energy spectrum distortions and seasonal variations in Sudbury Neutrino Observatory, Super-Kamiokande, and BOREXINO. Those predictions are given in detail, emphasizing that the vacuum solutions are distinguishable from the MSW solutions.

Various theoretical uncertainties in the standard solar model and in the Mikheyev-Smirnov-Wolfenstein (MSW) analysis are discussed. It is shown that two methods of estimating the solar neutrino flux uncertainties are equivalent: (a) a simple parametrization of the uncertainties using the core temperature and the nuclear production cross sections; (b) the Monte Carlo method of Bahcall and Ulrich. In the MSW analysis, we emphasize proper treatments of correlation of theoretical uncertainties between flux components and between different detectors, the Earth effect, and multiple solutions in a combined $\chi^2$ procedure. The MSW solutions for various standard and nonstandard solar models are also shown. The MSW predictions of the global solutions for the future solar neutrino experiments are given, emphasizing the measurement of the energy spectrum and the day-night effect in Sudbury Neutrino Observatory and Super-Kamiokande to distinguish the two solutions.

We update the analysis of the MSW and general astrophysical solutions to the combined solar neutrino observations by including the GALLEX II result. We also show that our parametrized flux uncertainties are equivalent to the Monte-Carlo results of Bahcall and Ulrich.

We consider the most general solar model, using the neutrino fluxes as free parameters constrained only by the solar luminosity, and show that the combined solar neutrino data exclude any astrophysical solution at 98\% C.L.\ Our best fit to the $^7$Be and $^8$B fluxes is respectively $<$7\% and 37$\pm$4\% of the standard solar model prediction, but only with a large $\chi^2$ (5.6 for 1 d.f.). This best fit to the fluxes contradicts explicit nonstandard solar models, which generally reduce the $^8$B flux more than the $^7$Be. Those models are well parameterized by a single parameter, the central temperature.

We consider the Earth effect in the MSW analysis of the Homestake, Kamiokande, GALLEX, and SAGE solar neutrino experiments. Using the time-averaged data and assuming two-flavor oscillations, the large-angle region of the combined fit extends to much smaller angles (to $\sin^22\theta \simeq 0.1$) than when the Earth effect is ignored. However, the additional constraint from the Kamiokande II day-night data excludes most of the parameter space sensitive to the Earth effect independent of astrophysical uncertainties, and leaves only a small large-angle region close to maximal mixing at 90\% C.L. The nonadiabatic solution remains unaffected by the Earth effect and is still preferred. Both theoretical and experimental uncertainties are included in the analysis.

Constraints on the core temperature (T_c) of the Sun and on neutrino- oscillation parameters are obtained from the existing solar neutrino data, including the recent GALLEX and Kamiokande III results. (1) A purely astrophysical solution to the solar neutrino problem is strongly disfavored by the data: the best fit in a cooler Sun model requires an 8% reduction in T_c, but the chi-sqaured test rejects this hypothesis at 99.99% C.L., suggesting new neutrino physics. (2) Assuming the Standard Solar Model (SSM) and MSW oscillations, the MSW parameters are constrained to two small regions: one in the non-adiabatic region and the other in the large-mixing region. The non-adiabatic solution gives a considerably better fit. For nu_e oscillations into sterile neutrinos, the allowed region (90%) is constrained to non- adiabatic oscillations. As long as the SSM is assumed, the neutrino mixing angles are at least four times larger, or considerably smaller, than the corresponding quark mixing angles. (3) Allowing both MSW oscillations and a non-standard core temperature, a) the experiments determine the core temperature at the 5% level, yielding a value consistent with the SSM prediction. b) When T_c is used as a free parameter, the allowed MSW region is broadened: a cooler Sun (T_c=0.95) allows mass and mixing implied by the SUSY SO(10) GUT, while a warmer Sun (T_c=1.05) allows parameter space suggested by intermediate-scale SO(10) GUTs. Superstring-inspired models are consistent with all solutions. (4) From the narrowed parameter space, we predict the neutrino spectral shape which should be observed in SNO. Throughout the calculation we use the latest Bahcall-Pinsonneault SSM, and include nuclear and astrophysical uncertainties in a simplified, but physically transparent way.