We have evolved 10,000 solar models using 21 input parameters that are randomly drawn for each model from separate probability distributions for every parameter. We use the results of these models to determine the theoretical uncertainties in the predicted surface helium abundance, the profile of the sound speed versus radius, the profile of the density versus radius, the depth of the solar convective zone, the eight principal solar neutrino fluxes, and the fractions of nuclear reactions that occur in the CNO cycle or in the three branches of the p-p chains. We also determine the correlation coefficients of the neutrino fluxes for use in analysis of solar neutrino oscillations. Our calculations include the most accurate available input parameters, including radiative opacity, equation of state, and nuclear cross sections. We incorporate both the recently determined heavy element abundances recommended by Asplund, Grevesse & Sauval (2005) and the older (higher) heavy element abundances recommended by Grevesse & Sauval (1998). We present best-estimates of many characteristics of the standard solar model for both sets of recommended heavy element compositions.

We have evolved a series of thirteen complete solar models that utilize different assumed heavy element compositions. Models that are based upon the heavy element abundances recently determined by Asplund, Grevesse, and Sauval (2005) are inconsistent with helioseismological measurements. However, models in which the neon abundance is increased by 0.4-0.5 dex to log N(Ne) = 8.29 +- 0.05 (on the scale in which log N(H) = 12) are consistent with the helioseismological measurements even though the other heavy element abundances are in agreement with the determinations of Asplund et al. (2005). These results sharpen and strengthen an earlier study by Antia and Basu (2005). The predicted solar neutrino fluxes are affected by the uncertainties in the composition by less than their 1sigma theoretical uncertainties.

We show that uncertainties in the values of the surface heavy element abundances of the Sun are the largest source of the theoretical uncertainty in calculating the p-p, pep, 8B, 13N, 15O, and 17F solar neutrino fluxes. We evaluate for the first time the sensitivity (partial derivative) of each solar neutrino flux with respect to the surface abundance of each element. We then calculate the uncertainties in each neutrino flux using `conservative (preferred)' and `optimistic' estimates for the uncertainties in the element abundances. The total conservative (optimistic) composition uncertainty in the predicted 8B neutrino flux is 11.6% (5.0%) when sensitivities to individual element abundances are used. The traditional method that lumps all abundances into a single quantity (total heavy element to hydrogen ratio, Z/X) yields a larger uncertainty, 20%. The uncertainties in the carbon, oxygen, neon, silicon, sulphur, and iron abundances all make significant contributions to the uncertainties in calculating solar neutrino fluxes; the uncertainties of different elements are most important for different neutrino fluxes. The uncertainty in the iron abundance is the largest source of the estimated composition uncertainties of the important 7Be and 8B solar neutrinos. Carbon is the largest contributor to the uncertainty in the calculation of the p-p, 13N, and 15O neutrino fluxes. However, for all neutrino fluxes, several elements contribute comparable amounts to the total composition uncertainty.

We construct solar models with the newly calculated radiative opacities from the Opacity Project (OP) and recently determined (lower) heavy element abundances. We compare results from the new models with predictions of a series of models that use OPAL radiative opacities, older determinations of the surface heavy element abundances, and refinements of nuclear reaction rates. For all the variations we consider, solar models that are constructed with the newer and lower heavy element abundances advocated by Asplund et al. (2005) disagree by much more than the estimated measuring errors with helioseismological determinations of the depth of the solar convective zone, the surface helium composition, the internal sound speeds, and the density profile. Using the new OP radiative opacities, the ratio of the 8B neutrino flux calculated with the older and larger heavy element abundances (or with the newer and lower heavy element abundances) to the total neutrino flux measured by the Sudbury Neutrino Observatory is 1.09 (0.87) with a 9% experimental uncertainty and a 16% theoretical uncertainty, 1 sigma errors.

I provide a summary of the current theoretical knowledge of solar neutrino fluxes as derived from precise solar models.

This encyclopedia article addresses questions like the following. How does the Sun shine? Does the neutrino have a mass? Are there weak interactions beyond those described by the standard model of particle physics?

We show that standard solar models are in good agreement with the helioseismologically determined sound speed and density as a function of solar radius, the depth of the convective zone, and the surface helium abundance, as long as those models do not incorporate the most recent heavy element abundance determinations. However, sophisticated new analyses of the solar atmosphere infer lower abundances of the lighter metals (like C, N, O, Ne, and Ar) than the previously widely used surface abundances. We show that solar models that include the lower heavy element abundances disagree with the solar profiles of sound speed and density as well as the depth of the convective zone and the helium abundance. The disagreements for models with the new abundances range from factors of several to many times the quoted uncertainties in the helioseismological measurements. The disagreements are at temperatures below what is required for solar interior fusion reactions and therefore do not significantly affect solar neutrino emission. If errors in thecalculated OPAL opacities are solely responsible for the disagreements, then the corrections in the opacity must extend from 2 times 10^6 K (R = 0.7R_Sun)to 5 times 10^6 K (R = 0.4 R_Sun), with opacity increases of order 10%.

We compare, using a three neutrino analysis, the allowed neutrino oscillation parameters and solar neutrino fluxes determined by the experimental data available Before and After Neutrino 2004. New data available after Neutrino2004 include refined KamLAND and gallium measurements. We use six different approaches to analyzing the KamLAND data. We present detailed results using all the available neutrino and anti-neutrino data for Delta m^2_{12}, tan^2 theta_{12}, sin^2 theta_{13}, and sin^2 eta (sterile fraction). Using the same complete data sets, we also present Before and After determinations of all the solar neutrino fluxes, which are treated as free parameters, an upper limit to the luminosity fraction associated with CNO neutrinos, and the predicted rate for a 7Be solar neutrino experiment. The 1 sigma (3 sigma) allowed range of Delta m^2_{21} = (8.2 +- 0.3) (^+1.0_-0.8)times 10^{-5} eV^2 is decreased by a factor of 1.7 (5), but the allowed ranges of all other neutrino oscillation parameters and neutrino fluxes are not significantly changed. Maximal mixing is disfavored at 5.8 sigma and the bound on the mixing angle theta_{13} is slightly improved to sin^2 theta_{13}<0.048 at 3 sigma. The predicted rate in a 7Be neutrino-electron scattering experiment is (0.665 +-0.015) of the rate implied by the BP04 solar model in the absence of neutrino oscillations. The corresponding predictions for p-p and pep experiments are, respectively, 0.707 {+0.011}{-0.013} and 0.644 {+0.011}{-0.013}. We derive upper limits to CPT violation in the weak sector by comparing reactor anti-neutrino oscillation parameters with neutrino oscillation parameters. We also show that the recent data disfavor at 91 % CL a proposed non-standard interaction description of solar neutrino oscillations.

We evaluate the logarithmic derivative of the depth of the solar convective zone with respect to the logarithm of the radiative opacity. We use this expression to show that the radiative opacity near the base of the solar convective zone (CZ) must be known to an accuracy of +- 1% in order to calculate the CZ depth to the accuracy of the helioseismological measurement, R(CZ) = (0.713 +- 0.001)R(Sun). The radiative opacity near the base of the CZ that is obtained from OPAL tables must be increased by about 21% in the Bahcall-Pinsonneault (2004) solar model if one wants to invoke opacity errors in order to reconcile recent solar heavy abundance determinations with the helioseismological measurement of R(CZ). We show that the radiative opacity near the base of the convective zone depends sensitively upon the assumed heavy element mass fraction, Z. The uncertainty in the measured value of Z is currently the limiting factor in our ability to calculate the depth of the CZ. Different state-of-the-art interpolation schemes using the existing OPAL tables yield opacity values that differ by 4% . We describe the finer grid spacings that are necessary to interpolate the radiative opacity to 1%. Uncertainties due to the equation of state do not significantly affect the calculated depth of the convective zone.

The three years 2001 to 2003 were the golden years of solar neutrino research. In this period, scientists solved a mystery with which they had been struggling for four decades. The solution turned out to be important for both physics and for astronomy. In this article, I tell the story of those fabulous three years.

We assess how well next generation neutrinoless double beta decay and normal neutrino beta decay experiments can answer four fundamental questions. 1) If neutrinoless double beta decay searches do not detect a signal, and if the spectrum is known to be inverted hierarchy, can we conclude that neutrinos are Dirac particles? 2) If neutrinoless double beta decay searches are negative and a next generation ordinary beta decay experiment detects the neutrino mass scale, can we conclude that neutrinos are Dirac particles? 3) If neutrinoless double beta decay is observed with a large neutrino mass element, what is the total mass in neutrinos? 4) If neutrinoless double beta decay is observed but next generation beta decay searches for a neutrino mass only set a mass upper limit, can we establish whether the mass hierarchy is normal or inverted? We base our answers on the expected performance of next generation neutrinoless double beta decay experiments and on simulations of the accuracy of calculations of nuclear matrix elements.

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.

Solar model predictions of 8B and p-p neutrinos agree with the experimentally-determined fluxes (including oscillations): phi(pp)_{measured} = (1.02 +- 0.02 +- 0.01)phi(pp)_{theory}, and phi(8B)_{measured} =(0.88 +- .04 +- 0.23)phi(8B)_{theory}, 1 sigma experimental and theoretical uncertainties, respectively. We use improved input data for nuclear fusion reactions, the equation of state, and the chemical composition of the Sun. The solar composition is the dominant uncertainty in calculating the 8B and CNO neutrino fluxes; the cross section for the 3He(4He, gamma)7Be reaction is the largest uncertainty for the calculated 7Be neutrino flux.

We analyze all available solar and related reactor neutrino experiments, as well as simulated future 7Be, p-p, pep, and ^8B solar neutrino experiments. We treat all solar neutrino fluxes as free parameters subject to the condition that the total luminosity represented by the neutrinos equals the observed solar luminosity (the `luminosity constraint'). Existing experiments show that the p-p solar neutrino flux is 1.02 +- 0.02 (1 sigma) times the flux predicted by the BP00 standard solar model; the 7Be neutrino flux is 0.93^{+0.25}_{-0.63} the predicted flux; and the ^8B flux is 1.01 +- 0.04 the predicted flux. The neutrino oscillation parameters are: Delta m^2 = 7.3^{+0.4}_{-0.6}\times 10^{-5} eV^2 and tan^2 theta_{12} = 0.41 +- 0.04. We evaluate how accurate future experiments must be to determine more precisely neutrino oscillation parameters and solar neutrino fluxes, and to elucidate the transition from vacuum-dominated to matter-dominated oscillations at low energies. A future 7Be nu-e scattering experiment accurate to +- 10 % can reduce the uncertainty in the experimentally determined 7Be neutrino flux by a factor of four and the uncertainty in the p-p neutrino flux by a factor of 2.5 (to +- 0.8 %). A future p-p experiment must be accurate to better than +- 3 % to shrink the uncertainty in tan^2 theta_{12} by more than 15 %. The idea that the Sun shines because of nuclear fusion reactions can be tested accurately by comparing the observed photon luminosity of the Sun with the luminosity inferred from measurements of solar neutrino fluxes. Based upon quantitative analyses of present and simulated future experiments, we answer the question: Why perform low-energy solar neutrino experiments?

I summarize 40 years of development of the standard solar model that is used to predict solar neutrino fluxes and then describe the current uncertainties in the predictions. I will also attempt to explain why it took so long, about three and a half decades, to reach a consensus view that new physics is being learned from solar neutrino experiments.

We use the strong nebular lines of O III, 5007 A and 4959 A, to set a robust upper limit on the time dependence of the fine structure constant. We find |alpha^{-1} d alpha(t)/dt | < 2*10^{-13} yr, corresponding to Delta alpha/alpha(0) = (0.7 +- 1.4) * 10^{-4} for quasars with 0.16 < z < 0.80 obtained from the SDSS Early Data Release. Using a blind analysis, we show that the upper limit given here is invariant with respect to 17 different ways of selecting the sample and analyzing the data. As a by-product, we show that the ratio of transition probabilities for the 5007 A and the 4959 A lines is2.99 +- 0.02, in good agreement with (but more accurate than) theoretical estimates. We compare and contrast the O III emission line method used here with the Many-Multiplet method that has been used recently to suggest evidence for a time-dependent alpha. In an Appendix, we have analyzed the larger quasar sample from the SDSS Data Release One Sample and find Delta alpha/alpha(0) = (1.2 +- 0.7) * 10^{-4}.

We derive the disk I-band luminosity function from the Zheng et al. sample of \~1400 disk M dwarfs observed with the Hubble Space Telescope. We adopt a Galactic-height-dependent color-magnitude relation to account for the metallicity gradient above the Galactic plane. The resultant I-band luminosity function peaks at M_I~9.5 and drops sharply toward M_I~10.5.

The energy spectra of ultra high energy cosmic rays reported by the AGASA, Fly's Eye, Haverah Park, HiRes, and Yakutsk experiments are all shown to be in agreement with each other for energies below 10^{20} eV (after small adjustments, within the known uncertainties, of the absolute energy scales). The data from HiRes, Fly's Eye, and Yakutsk are consistent with the expected flux suppression above 5\times 10^{19} eV due to interactions of cosmic rays with the cosmic microwave background, the Greisen-Zatsepin-Kuzmin (GZK) "supression," and are inconsistent with a smooth extrapolation of the observed cosmic ray energy spectrum to energies > 5\times 10^{19}$ eV. AGASA data show an excess of events above 10^{20} eV, compared to the predicted GZK suppression and to the flux measured by the other experiments.

We use the recently reported KamLAND measurements on oscillations of reactor anti-neutrinos, together with the data of previously reported solar neutrino experiments, to show that: (1) the total 8B neutrino flux emitted by the Sun is 1.00(1.0 \pm 0.06) of the standard solar model (BP00) predicted flux, (2) the KamLAND measurements reduce the area of the globally allowed oscillation regions that must be explored in model fitting by six orders of magnitude in the Delta m^2-tan^2 theta plane, (3) LMA is now the unique oscillation solution to a CL of 4.7sigma, (4) maximal mixing is disfavored at 3.1 sigma, (5) active-sterile admixtures are constrained to sin^2 eta<0.13 at 1 sigma, (6) the observed ^8B flux that is in the form of sterile neutrinos is 0.00^{+0.09}_{-0.00} (1 sigma), of the standard solar model (BP00) predicted flux, and (7) non-standard solar models that were invented to completely avoid solar neutrino oscillations are excluded by KamLAND plus solar at 7.9 sigma . We also refine quantitative predictions for future 7Be and p-p solar neutrino experiments.

We show that solar neutrino experiments set an upper limit of 7.8% (7.3% including the recent KamLAND measurements) to the fraction of energy that the Sun produces via the CNO fusion cycle, which is an order of magnitude improvement upon the previous limit. New experiments are required to detect CNO neutrinos corresponding to the 1.5% of the solar luminosity that the standard solar model predicts is generated by the CNO cycle.

This article describes the unconventional and occassionally stormy birth and early history of the HST Snapshot Program.

I summarize in four slides the 40 years of development of the standard solar model that is used to predict solar neutrino fluxes and then describe the current uncertainties in the predictions. I next dispel the misconception that the p-p neutrino flux is determined by the solar luminosity and present a related formula that gives, in terms of the p-p and 7Be neutrino fluxes, the ratio of the rates of the two primary ways of terminating the p-p fusion chain. I will also attempt to explain why it took so long, about three and a half decades, to reach a consensus view that new physics is being learned from solar neutrino experiments. Finally, I close with a personal confession.

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.

The 8B solar neutrino flux inferred from a global analysis of solar neutrino experiments is within 11% (1 sigma) of the predicted standard solar model value if only active neutrinos exist, but could be as large as 1.7 times the standard prediction if sterile neutrinos exist. We show that the total 8B neutrino flux (active plus sterile neutrinos) can be determined experimentally to about 10% (1 sigma) by combining charged current measurements made with the KamLAND reactor experiment and with the SNO CC solar neutrino experiment, provided the LMA neutrino oscillation solution is correct and the simulated performance of KamLAND is valid. Including also SNO NC data, the sterile component of the 8B neutrino flux can be measured by this method to an accuracy of about 12% (1 sigma) of the standard solar model flux. Combining Super-Kamiokande and KamLAND measurements and assuming the oscillations occur only among active neutrinos, the 8B neutrino flux can be measured to 6% (1 sigma); the total flux can be measured to an accuracy of about 9%. The total 7Be solar neutrino flux can be determined to an accuracy of about 28% (1 sigma) by combining measurements made with the KamLAND, SNO, and gallium neutrino experiments. One can determine the total 7Be neutrino flux to a one sigma accuracy of about 11% or better by comparing data from the KamLAND experiment and the BOREXINO solar neutrino experiment provided both detectors work as expected. The pp neutrino flux can be determined to about 15% using data from the gallium, KamLAND, BOREXINO, and SNO experiments.

With the goal of identifying signatures that select specific neutrino oscillation parameters, we test the robustness of global oscillation solutions that fit all the available solar and reactor experimental data. We use three global analysis strategies previously applied by different authors and also determine the sensitivity of the oscillation solutions to the critical nuclear fusion cross section, S_{17}(0), for the production of 8B. The favored solutions are LMA, LOW, and VAC in order of g.o.f. The neutral current to charged current ratio for SNO is predicted to be 3.5 +- 0.6 (1 sigma), which is separated from the no-oscillation value of 1.0 by much more than the expected experimental error. The predicted range of the day-night difference in charged current rates is (8.2 +- 5.2)% and is strongly correlated with the day-night effect for neutrino-electron scattering. A measurement by SNO of either a NC to CC ratio > 3.3 or a day-night difference > 10%, would favor a small region of the currently allowed LMA neutrino parameter space. The global oscillation solutions predict a 7Be neutrino-electron scattering rate in BOREXINO and KamLAND in the range 0.66 +- 0.04 of the BP00 standard solar model rate, a prediction which can be used to test both the solar model and the neutrino oscillation theory. Only the LOW solution predicts a large day-night effect(< 42%) in BOREXINO and KamLAND. For the KamLAND reactor experiment, the LMA solution predicts 0.44 of the standard model rate; we evaluate 1 sigma and 3 sigma uncertainties and the first and second moments of the energy spectrum.

We show that the combined data from solar neutrino experiments and from the KamLAND reactor neutrino experiment can establish an upper limit on, or detect, potential CPT violation in the neutrino sector of order 10^{-20} GeV to 10^{-21} GeV.

A specific linear combination of the total solar neutrino fluxes must equal the measured solar photon luminosity if nuclear fusion reactions among light elements are responsible for solar energy generation. This luminosity constraint, previously used in a limited form in testing the no neutrino oscillation hypothesis, is derived in a generality that includes all of the relevant solar neutrino fluxes and which is suitable for analyzing the results of many different solar neutrino experiments. With or without allowing for neutrino oscillations, the generalized luminosity constraint can be used in future analyses of solar neutrino data. Accurate numerical values for the linear coefficients are provided.

The minimal standard electroweak model can be tested by allowing all the solar neutrino fluxes, with undistorted energy spectra, to be free parameters in fitting the measured solar neutrino event rates, subject only to the condition that the total observed luminosity of the sun is produced by nuclear fusion. The rates of the five experiments prior to SNO (chlorine, Kamiokande, SAGE, GALLEX, Super-Kamiokande) cannot be fit by an arbitrary choice of undistorted neutrino fluxes at the level of 2.5 sigma (formally 99% C.L.). Considering just SNO and Super-Kamiokande, the discrepancy is at the 3.3 sigma level(10^{-3} C.L.). If all six experiments are fit simultaneously, the formal discrepancy increases to 4 sigma (7*10^{-5} C.L.). If the relative scaling in temperature of the nuclear reactions that produce 7Be and 8B neutrinos is taken into account, the formal discrepancy is at the 7.4 sigma level.

For active and sterile neutrinos, we present the globally allowed solutions for two neutrino oscillations. We include the SNO CC measurement and all other relevant solar neutrino and reactor data. Five active neutrino oscillation solutions (LMA, LOW, SMA, VAC, and Just So2) are currently allowed at 3 sigma; three sterile neutrino solutions (Just So2, SMA, and VAC) are allowed at 3 sigma. The goodness of fit is satisfactory for all eight solutions. We also investigate the robustness of the allowed solutions by carrying out global analyses with and without: 1) imposing solar model constraints on the 8B neutrino flux, 2) including the Super-Kamiokande spectral energy distribution and day-night data, 3) including a continuous mixture of active and sterile neutrinos, 4) using an enhanced CC cross section for deuterium (due to radiative corrections), and 5) a optimistic, hypothetical reduction by a factor of three of the error of the SNO CC rate. For every analysis strategy used in this paper, the most favored solutions all involve large mixing angles: LMA, LOW, or VAC. The favored solutions are robust, but the presence at 3 sigma of individual sterile solutions and the active Just So2 solution is sensitive to the analysis assumptions.

I discuss why we need solar neutrino experiments below 1 MeV. I also express my prejudices about the desired number and types of such experiments, emphasizing the importance of p-p solar neutrino experiments

We present a global analysis of all the available solar neutrino data treating consistently the 8B and hep neutrino fluxes as free parameters. The analysis reveals at 99.7% C.L. eight currently-allowed discrete regions in two-neutrino oscillation space, five regions corresponding to active neutrinos and three corresponding to sterile neutrinos. Most of the allowed solutions are robust with respect to changes in the analysis procedure, but the traditional vacuum solution is fragile. The globally-permitted range of the 8B neutrino flux, 0.45 to 1.95 in units of the BP2000 flux, is comparable to the 3 sigma range allowed by the standard solar model. We discuss the implications for SNO of a low mass, Delta m^2 ~ 6 times 10^{-12} eV^2, vacuum oscillation solution, previously found by Raghavan, and by Krastev and Petcov, but absent in recent analyses that included Super-Kamiokande data. For the SNO experiment, we present refined predictions for the charged-current rate and the ratio of the neutral-current rate to charged-current rate. The predicted charged-current rate can be clearly distinguished from the no-oscillation rate only for the LMA solution. The predicted ratio of the neutral-current rate to charged-current rate is distinguishable from the no-oscillation ratio for the LMA, SMA, LOW, and VAC solutions for active neutrinos.

We calculate accurate solar models and report the detailed time dependences of important solar quantities. We use helioseismology to constrain the luminosity evolution of the sun and report the discovery of semi-convection in evolved solar models that include diffusion. In addition, we compare the computed sound speeds with the results of p-mode observations by BiSON, GOLF, GONG, LOWL, and MDI instruments. We contrast the neutrino predictions from a set of eight standard-like solar models and four deviant (or deficient) solar models with the results of solar neutrino experiments. For solar neutrino and for helioseismological applications, we present present-epoch numerical tabulations of characteristics of the standard solar model as a function of solar radius, including the principal physical and composition variables, sound speeds, neutrino fluxes, and functions needed for calculating solar neutrino oscillations.

We study a sample of about 1400 disk M dwarfs that are found in 148 fields observed with the Wide Field Camera 2 (WFC2) on the Hubble Space Telescope and 162 fields observed with pre-repair Planetary Camera 1 (PC1), of which 95 of the WFC2 fields are newly analyzed. The method of maximum likelihood is applied to derive the luminosity function and the Galactic disk parameters. At first, we use a local color-magnitude relation and a locally determined mass-luminosity relation in our analysis. The results are consistent with those of previous work but with considerably reduced statistical errors. These small statistical errors motivate us to investigate the systematic uncertainties. Considering the metallicity gradient above the Galactic plane, we introduce a modified color-magnitude relation that is a function of Galactic height. The resultant M dwarf luminosity function has a shape similar to that derived using the local color-magnitude relation but with a higher peak value. The peak occurs at $M_V \sim 12$ and the luminosity function drops sharply toward $M_V \sim 14$. We then apply a height-dependent mass-luminosity function interpolated from theoretical models with different metallicities to calculate the mass function. Unlike the mass function obtained using local relations, which has a power-law index $\alpha = 0.47$, the one derived from the height-dependent relations tends to be flat ($\alpha = -0.10$). The resultant local surface density of disk M dwarfs (12.2 +/- 1.6 M_sun pc^{-2}) is somewhat smaller than the one obtained using local relations (14.3 +/- 1.3 M_sun pc^{-2}). Our measurement favors a short disk scale length, H = 2.75 +/- 0.16 (statistical) +/- 0.25 (systematic) kpc.

This historical essay describes the struggle to understand how the sun shines. The story begins with the disagreement between Kelvin and Darwin over the age of the sun in the middle of the 19th century and continues through the revelations of solar neutrino experiments at the conclusion of the 20th century.

We elucidate the physical basis for the upper bound on high energy neutrino fluxes implied by the observed cosmic ray flux. We stress that the bound is valid for neutrinos produced either by p,gamma reactions or by p-p(n) reactions in sources which are optically thin for high energy protons to photo-meson and nucleon-meson interactions. We show that the upper bound is robust and conservative. The Waxman-Bahcall bound overestimates the most likely neutrino flux by a factor ~ 5/tau, for small optical depths tau. The upper limit cannot be plausibly evaded by invoking magnetic fields, optically thick AGNs, or large hidden fluxes of extragalactic protons. We describe the implications of the bound for future experiments including the AMANDA, ANTARES, Auger, ICECUBE, NESTOR, and OWL/AIRWATCH detectors.

Neutrino oscillation scenarios predict correlations, and zones of avoidance, among measurable quantities such as spectral energy distortions, total fluxes, time dependences, and flavor content. The comparison of observed and predicted correlations will enhance the diagnostic power of solar neutrino experiments. A general test of all presently-allowed (two neutrino) oscillation solutions is that future measurements must yield values outside the predicted zones of avoidance. To illustrate the discriminatory power of the simultaneous analysis of multiple observables, we map currently allowed regions of neutrino masses and mixing angles onto planes of quantities measurable with the Sudbury Neutrino Observatory (SNO). We calculate the correlations that are predicted by vacuum and MSW (active and sterile) neutrino oscillation solutions that are globally consistent with all available neutrino data. We derive approximate analytic expressions for the dependence of individual observables and specific correlations upon neutrino oscillations parameters. We also discuss the prospects for identifying the correct oscillation solution using multiple SNO observables.

We review five different derivations that demonstrate that the Salpeter formula for the plasma corrections to fusion rates is valid at the center of the sun with insignificant errors (~ percent). We point out errors in several recent papers that have obtained a variety of answers, some even with the wrong sign or the wrong functional dependence.

I summarize the first four decades of solar neutrino research and suggest what may be possible to learn with extragalactic neutrinos and with solar neutrinos in the next decade.

A gamma-ray burst fireball is likely to contain an admixture of neutrons, in addition to protons, in essentially all progenitor scenarios. Inelastic collisions between differentially streaming protons and neutrons in the fireball produce muon neutrinos (antineutrinos) of ~ 10 GeV as well as electron neutrinos (antineutrinos) of ~ 5 GeV, which could produce ~ 7 events/year in kilometer cube detectors, if the neutron abundance is comparable to that of protons. Photons of ~ 10 GeV from pi-zero decay and ~ 100 MeV electron antineutrinos from neutron decay are also produced, but will be difficult to detect. Photons with energies < 1 MeV from shocks following neutron decay produce a characteristic signal which may be distinguishable from the proton-related MeV photons.

We calculate the range of predicted values for 10 quantities that will be measured by the Sudbury Neutrino Observatory (SNO). We use neutrino oscillation solutions (vacuum and MSW; active and sterile neutrinos) that are globally consistent with all available neutrino data and estimate realistic theoretical and experimental uncertainties. The neutral current to charged current double ratio is predicted to be more than 9 sigma from the no-oscillation solution for all of the currently favored neutrino oscillation solutions. The best-fit oscillation solutions predict a CC day-night rate difference between -0.1% and +12.5% and a NC day-night difference < 0.01%. We present also the predicted range for the first and the second moments of the charged current electron recoil energy spectrum, the charged current, the neutral current, and the neutrino-electron scattering rates, the seasonal dependence of the charged current rate, and the double ratio of neutrino-electron scattering rate to charged current rate.

Clustering of cosmic-ray sources affects the flux observed beyond the cutoff imposed by the cosmic microwave background and may be important in interpreting the AGASA, Fly's Eye, and HiRes data. The standard deviation, sigma, in the predicted number, N, of events above 10^{20} eV is sigma/N = 0.9(r_0/10 Mpc)^{0.9}, where r_0 is the unknown scale length of the correlation function (r_0 = 10 Mpc for field galaxies). Future experiments will allow the determination of r_0 through the detection of anisotropies in arrival directions of ~ 10^{20} eV cosmic-rays over angular scales of Theta ~ r_0/30 Mpc.

We show that a significant fraction of the energy of a gamma-ray burst(GRB) is probably converted to a burst of 10^{17}-10^{19} eV neutrinos and multiple GeV gammas that follow the GRB by > 10 s . If, as previously suggested, GRB's accelerate protons to ~10^{20} eV, then both the neutrinos and the gammas may be detectable.

How did neutrino astronomy evolve? Are there any useful lessons for astronomers and physicists embarking on new observational ventures today? We answer the first question; the reader can can decide if there are any useful parallels for other fields.

The ratio of the measured to the predicted standard model CC event rates in SNO will be 0.47 if no oscillations occur. The best-fit active oscillation predictions for the CC ratio are: 0.35-39 (MSW) and 0.38-42 (vacuum) (all for a 5 MeV energy threshold), typically about 20% less than the no-oscillation expectation. We calculate the predicted ratios for six active and sterile neutrino oscillation solutions allowed at 99% CL and determine the dependence of the ratios on energy threshold. If the high-energy anomaly observed by SuperKamiokande is due to an enhanced hep flux, MSW active solutions predict that out of a total of 5000 CC events above 5 MeV in SNO between 49 and 54 events will be observed above 13 MeV whereas only 19 events are expected for no-oscillations and a nominal standard hep flux.

I summarize the current state of solar neutrino research and then give my answer to the question: What should we do next?

Recent results on solar neutrinos provide hints that the LMA MSW solution could be correct. We perform accurate calculations for potential `smoking-gun' effects of the LMA solution in the SuperKamiokande solar neutrino experiment, including: (1) an almost constant reduction of the standard recoil electron energy spectrum (with a weak, < 10%, relative increase below 6.5 MeV); (2) an integrated difference in day-night rates (2% to 14%); (3) an approximately constant zenith-angle dependence of the nighttime event rate; (4) a new test for the difference in the shape of the equally-normalized day-night energy spectra (~ 1%); and (5) annual variations of the signal due to the regeneration effect (~ 6 times smaller than the integrated day-night effect). We also establish a relation between the integrated day-night asymmetry and the seasonal asymmetry due to LMA regeneration. As a cautionary example, we simulate the effect of an absolute energy calibration error on the shape (distortion) of the recoil energy spectrum. We compare LMA predictions with SuperKamiokande data and discuss the possibilities for distinguishing experimentally between LMA and vacuum oscillations.

We investigate systematic uncertainties in determining the profiles of the solar sound speed, density, and adiabatic index by helioseismological techniques. We find that rms uncertainties-averaged over the sun of ~ 0.2%-0.4% are contributed to the sound speed profile by each of three sources: 1)the choice of assumed reference model, 2) the width of the inversion kernel, and 3) the measurements errors. The density profile is about an order of magnitude less well determined by the helioseismological measurements. The profile of the adiabatic index is determined to an accuracy of about 0.2% . We find that even relatively crude reference models yield reasonably accurate solar parameters.

Observations with the Wide-Field/Planetary Camera-2 of the Hubble Space Telescope (HST) are presented for three radio-loud quasars: 3C 48 (z=0.367), B2 1425+267 (z=0.366), and 3C 345 (z=0.594). All three quasars have luminous (~4 L^*) galaxies as hosts, which are either elliptical (B2 1425+267 and 3C 345) or interacting (3C 48), and all hosts are 0.5 - 1.0 mag bluer in (V-I) than other galaxies with the same overall morphology at similar redshifts to the quasars. The host of 3C 48 has many H II regions and a very extended tidal tail. All nine of the radio-loud quasars studied here and in Bahcall et al. (1997) either have bright elliptical hosts or occur in interacting systems. There is a robust correlation between the radio emission of the quasar and the luminosity of host galaxy; the radio-loud quasars reside in galaxies that are on average about 1 mag brighter than hosts of the radio-quiet quasars.

We describe a program of star counts in the range 19 < I < 26 made with the WFPC cameras aboard the Hubble Space Telescope. Red (V-I > 1.0) stars at these magnitudes are primarily disk and spheroid M dwarfs. The stars are found both on dedicated images as part of the parallel program and by using appropriate archive data. We measure the faint end of the luminosity functions of the disk and spheroid (i.e. stellar halo). We measure the low mass end of the mass function and show that M dwarfs do not dominate the total disk or spheroid mass. We place strong I-band constraints on the amount of halo dark matter in the form of low mass stars (such as M dwarfs or cool white dwarfs). The disk and spheroid contribute only a minor amount of optical depth toward the Magellanic clouds.

We report on Hubble Space Telescope (HST) observations of the host galaxy of GRB 971214, about four months after the burst. The redshift of the proposed host galaxy at z=3.418, combined with optical and radio observations of the burst afterglow, implies the extremely large isotropic energy release from the burst in gamma-rays of approximately 3x10^53 ergs, some two orders of magnitude higher than the previously commonly assumed numbers. The positional offset between the optical transient observed at the Keck telescope and the centroid of the proposed host galaxy in the HST image is 0.14 -+0.07 arcsec. We find no evidence in our deep HST image for a chance foreground galaxy superposed along the line of sight to the proposed host at z=3.418. The morphology and photometric properties of this galaxy, such as the total flux, morphology, radial surface profile and scale length, are typical as compared to other, spectroscopically confirmed z>3 galaxies.

We show that cosmic-ray observations set a model-independent upper bound to the flux of high-energy, > 10^14 eV, neutrinos produced by photo-meson (or p-p) interactions in sources of size not much larger than the proton photo-meson (or pp) mean-free-path. The bound applies, in particular, to neutrino production by either AGN jets or GRBs. This upper limit is two orders of magnitude below the flux predicted in some popular AGN jet models, but is consistent with our predictions from GRB models. We discuss the implications of these results for future km^2 high-energy neutrino detectors.

I review recent developments that affect standard solar model predictions of solar neutrino fluxes.

We present a catalogue of absorption lines obtained from the analysis of the ultra-violet spectra of 66 quasars. The data were acquired with the Faint Object Spectrograph of the HST as part of the Quasar Absorption Line Survey, a Key Project for the first four cycles of HST observations. This is the third of a series of catalogues of absorption lines produced from the survey and increases the number of quasars whose higher resolution (R=1300) spectra we have published from 17 to 83. The general properties and execution of the survey are reviewed, including descriptions of the final sample of observed objects and the algorithmic processes used to construct the catalogue. The detection of a single damped Ly-a system in a path length of Delta_z=49 yields an observed number of damped systems per unit redshift of (dN/dz)_{damp}(z=0.58)=0.020 with 95% confidence boundaries of 0.001 to 0.096 systems per unit redshift. We include notes on our analysis of each of the observed quasars and the absorption systems detected in each spectrum. Some especially interesting systems include low redshift Ly-a absorbers suitable for extensive follow-up observations (in the spectra of TON28 and PG1216+069), possibly physically associated pairs of extensive metal line absorption systems (e.g., in the spectrum of PG0117+213), and systems known to be associated with galaxies (e.g., in the spectrum of 3C232).

We determine the neutrino parameters for MSW and vacuum oscillations (active and sterile neutrinos) that are allowed by the separate, and collective, imposition of the constraints from total event rates in the chlorine, GALLEX, SAGE, and SuperKamiokande experiments (504 days), the SuperKamiokande electron energy spectrum, and the SuperKamiokande zenith-angle dependence. The small mixing angle MSW solution is acceptable at the 7% C.L. (8% for sterile nu's) and the vacuum solution is acceptable at the 6% C.L. . The best-fit global MSW solution for active neutrinos is: Delta m^2 = 5 x 10^-6 eV^2, sin^2 (2 theta) = 5.5 x 10^{-3} (and for sterile neutrinos: Delta m^2 = 4 x 10^-6 eV^2, sin^2 (2 theta) = 7 x 10^-3). For vacuum oscillations, the best-fit solution is: Delta m^2 = 6.5 x 10^-11 eV^2, sin^2 (2 theta) = 0.75 . An arbitrary combination of undistorted (no oscillations) pp, 7Be, 8B, and CNO neutrino fluxes is inconsistent with the combined data sets at the 3.5 sigma C.L., independent of astrophysical considerations. We use improved calculations of solar model fluxes, neutrino absorption cross sections and energy spectra, and a detailed evaluation of regeneration effects.

If the low energy cross section for 3He + p goes to 4He + e + nu_e, the `hep' reaction, is > 20 times larger than the best (but uncertain) theoretical estimates, then this reaction could significantly influence the electron energy spectrum produced by solar neutrino interactions and measured in the SuperKamiokande, SNO, and ICARUS experiments. We compare predicted energy spectra for different assumed hep fluxes and different neutrino oscillation scenarios with the observed SuperKamiokande spectrum. The spectra with enhanced hep contributions provide better fits to the SuperKamiokande data.

We present an analysis of the evolution of the Ly alpha absorption lines in the redshift interval 0.0 to \sim 1.5 based upon a sample of 987 lines identified in the spectra of 63 QSOs obtained with the Faint Object Spectrograph on the Hubble Space Telescope. Fits to the number of absorbers per unit redshift dN/dz = A \times (1+z)^\gamma yield values of \gamma in the range 0.1 to 0.3, decidedly flatter than results from groundbased data pertaining to the redshift range z > 1.7. The combination of the space and groundbased data suggest a marked transition in the rate of evolution of the Ly alpha lines at a redshift of about 1.7. The distribution in redshift is well represented by a power law and the distribution of equivalent widths above a rest equivalent width of 0.1 \AA is fit well by an exponential. The sample of Ly alpha systems which has ions from heavy elements has a much larger value of \gamma than a sample of lines of comparable strength without heavy elements. We argue that this is not likely due to either line misidentification or incomplete spectral coverage and suggest that it is due to rapid evolution of the very highest column density systems, rather than differences in metallicity. We also find evidence that the rate of evolution increases with increasing equivalent width. We compare our results for the variation of line density with redshift to recent numerical simulations of Ly alpha absorbers, and find fairly good agreement. Our results for the stronger Ly alpha lines in our sample are compatible with the recent analysis of the absorber--galaxy correlation by Chen et al. 1998 but for the weaker lines our results suggest that the association between absorbers and galaxies is different from that for the stronger lines.

We present F555W (V), F439W (B), and F336W (U) photometry of 9507 stars in the central 2' of the dense, post core collapse cluster M30 (NGC 7099) derived from HST WFPC2 images. These data are used to study the mix of stellar populations in the central region of the cluster. Forty eight blue straggler stars are identified and are found to be strongly concentrated towards the center. The specific frequency of blue stragglers, F_BSS\equiv N(BSS)/N(V<V_HB+2), is 0.25\pm 0.05 in the inner region of M30 (r<20''), significantly higher than the frequency found in other clusters. The shape of M30's BSS luminosity function resembles the prediction of the collisional formation model and is inconsistent with the binary merger model, of Bailyn and Pinsonneault. An unusually blue star (B=18.6,B-V=-0.97), possibly a cataclysmic variable based on its color, is found about 1.2'' from the center. Bright RG stars appear to be depleted by a factor of 2--3 in the inner r<10'' relative to fainter giants, subgiants, and MSTO stars (95% significance). We confirm that there is a radial gradient in the color of the cluster light, going from B-V=0.82 at r=1' to B-V=0.45 in the central 10''. The central depletion of the bright RG is responsible for about half of the observed color gradient; the rest of the gradient is caused by the relative underabundance of faint red MS stars near the center (presumably a result of mass segregation). The luminosity function of M30's evolved stars does not match the luminosity function shape derived from standard stellar evolutionary models: the ratio of the number of bright giants to the number of turnoff stars in the cluster is 30% higher than predicted by the model (3.8\sigma effect), roughly independent of RG brightness over the range M_V=-2 to +2. (abridged)

Solar neutrino fluxes and sound speeds are calculated using a systematic reevaluation of nuclear fusion rates. The largest uncertainties are identified and their effects on the solar neutrino fluxes are estimated.

We review and analyze the available information for nuclear fusion cross sections that are most important for solar energy generation and solar neutrino production. We provide best values for the low-energy cross-section factors and, wherever possible, estimates of the uncertainties. We also describe the most important experiments and calculations that are required in order to improve our knowledge of solar fusion rates.

We analyze 166 spheroid subdwarfs $(6.5<M_V<14.5)$ found in 53 fields observed with the Wide Field Camera on the {\it Hubble Space Telescope}. The fields cover 221 square arcmin over a wide range of directions. The spheroid luminosity function (LF) is inconsistent at about the 3 sigma level with the local spheroid LF of Dahn et al. even when the normalization of the latter is corrected to take account of the latest data on spheroid kinematics. The difference may reflect systematic errors in one of the two studies or features of the spheroid spatial distribution that are not included in the simplest models. The mass function, which shows no obvious structure, can be represented by a power law, d N/d ln M ~ M^\alpha, with \alpha=0.25 +/- 0.32 over the mass range 0.71 Msun > M >0.09 Msun. The spheroid therefore does not contribute significantly to microlensing unless the mass function changes slope dramatically in the substellar range. The total local mass density of spheroid stars (including remnants and unseen binary companions) is \rho ~ 6.4E-5 Msun/pc^3, with an uncertainty of about 50%. The power-law indices \alpha=0.25 for the spheroid and \alpha=0.44 for the disk (both uncorrected for binaries) are similar to those of globular clusters of moderate to high metallicity.

We test the Salpeter formalism for the electron screening of the solar proton-proton fusion reaction by solving numerically the relevant Schrodinger equation. We evaluate exactly the square of the overlap integral of the two-proton wave function and the deuteron wave function and compare with the usual analytic approximation. The usual WKB solution agrees with the numerical result to $O(10^{-4})$. The WKB approximation should be even more precise for the other nuclear fusion reactions in the $pp$ chain and CNO cycles.

We evaluate the effect of electrostatic screening by ions and electrons on low-Z thermonuclear reactions in the sun. We use a mean field formalism and calculate the electron density of the screening cloud using the appropriate density matrix equation of quantum statistical mechanics. Because of well understood physical effects that are included for the first time in our treatment, the calculated enhancement of reaction rates does not agree with the frequently used interpolation formulae. Our result does agree, within small uncertainties, with Salpeter's weak screening formula. If weak screening is used instead of the commonly employed screening prescription of Graboske et al., the predicted $^8$B neutrino flux is increased by 7% and the predicted chlorine rate is increased by 0.4 SNU.

The first lecture describes the conflicts between the combined standard model (minimal electroweak and solar) predictions and the results of solar neutrino experiments. The second lecture discusses the possibilities for detecting with solar neutrinos ``smoking gun'' indications of departures from minimal electroweak theory. Viewgraphs and related data are available at http://www.sns.ias.edu/~jnb .

Neutrino absorption cross sections for 71Ga are calculated for all solar neutrino sources with standard energy spectra, and for laboratory sources of 51Cr and 37Ar; the calculations include, where appropriate, the thermal energy of fusing solar ions and use improved nuclear and atomic data. The ratio, R, of measured (in GALLEX and SAGE) to calculated 51Cr capture rate is R = 0.95 +/- 0.07 (exp)} + ^{+0.04}_{-0.03} (theory). Cross sections are also calculated for specific neutrino energies chosen so that a spline fit determines accurately the event rates in a gallium detector even if new physics changes the energy spectrum of solar neutrinos. Theoretical uncertainties are estimated for cross sections at specific energies and for standard neutrino energy spectra. Standard energy spectra are presented for pp and CNO neutrino sources in the appendices. Neutrino fluxes predicted by standard solar models, corrected for diffusion, have been in the range 120 SNU to 141 SNU since 1968.

We calculate accurately the number of solar neutrino events expected as a function of solar zenith angle, with and without neutrino oscillations, for detectors at the locations of Super-Kamiokande, SNO, and the Gran Sasso National Laboratory. Using different earth models to estimate geophysical uncertainties, and different solar models to estimate solar uncertainties, we evaluate distortions predicted by the MSW effect in the zenith angle distributions of solar neutrino events. The distortions are caused by oscillations and by $\nu-e$ interactions in the earth that regenerate $\nu_e$ from $\nu_\mu$ or $\nu_\tau$. We show that the first two moments of the zenith-angle distribution are more sensitive to the small mixing angle MSW solution than the conventionally studied day-night asymmetry. We present iso-sigma contours that illustrate the potential of Super-Kamiokande, SNO, BOREXINO, ICARUS and HERON/HELLAZ for detecting the earth regeneration effect at their actual locations (and at the equator). MSW solutions favored by the four pioneering solar neutrino experiments predict characteristic distortions for Super-Kamiokande, SNO, BOREXINO, and ICARUS that range from being unmeasurably small to $> 5\sigma$ (stat) after only a few years of observations.

In a series of recent papers, Ivanov et al. and Oberhummer et al. have calculated the rate for the $p + p \to d + e^+ + \nu_e$ reaction with a zero-range four-fermion effective interaction and find a result 2.9 times higher than the standard value calculated from non-relativistic potential theory. Their procedure is shown to give a wrong answer because their assumed interaction disagrees with low-energy $pp$ scattering data.

Localized differences between the real sun and standard solar models are shown to be small. The sound speeds of the real and the standard model suns typically differ by less than 0.3% in the solar core.

This article summarizes the process by which astronomers and astrophysicists reached a consensus view as to what are the highest priority scientific projects for the decade of the 1990s. I also review the remarkable success of these recommendations.

For solar conditions, we numerically integrate the density matrix equation for a thermal electron in the field of a Be-7 ion and other plasma ions and smeared-out electrons. Our results are in agreement with previous calculations that are based on a different physical picture, a picture which postulates the existence of distinct continuum and bound state orbits for electrons. The density matrix calculation of the electron capture rate is independent of the nature of electron states in the solar plasma. To within a 1% accuracy, the effects of screening can be described at high temperatures by a Salpeter-like factor which can be derived from the density matrix equation. The total theoretical uncertainty in the electron capture rate is about 2%.

This talk compares standard model predictions for solar neutrino experiments with the results of actual observations. Here `standard model' means the combined standard model of minimal electroweak theory plus a standard solar model. I emphasize the importance of recent analyses in which the neutrino fluxes are treated as free parameters, independent of any constraints from solar models, and the stunning agreement between the predictions of standard solar models and helioseismological measurements.

Observations suggest that $\gamma$-ray bursts (GRBs) are produced by the dissipation of the kinetic energy of a relativistic fireball. We show that a large fraction, $\ge 10%$, of the fireball energy is expected to be converted by photo-meson production to a burst of $\sim10^{14} eV$ neutrinos. A km^2 neutrino detector would observe at least several tens of events per year correlated with GRBs, and test for neutrino properties (e.g. flavor oscillations, for which upward moving $\tau$'s would be a unique signature, and coupling to gravity) with an accuracy many orders of magnitude better than is currently possible.

The sound speeds of solar models that include element diffusion agree with helioseismological measurements to a rms discrepancy of better than 0.2% throughout almost the entire sun. Models that do not include diffusion, or in which the interior of the sun is assumed to be significantly mixed, are effectively ruled out by helioseismology. Standard solar models predict the measured properties of the sun more accurately than is required for applications involving solar neutrinos.

Observations with the Wide-Field Camera of the Hubble Space Telescope are presented for a representative sample of 20 intrinsically luminous quasars with redshifts smaller than 0.30. These observations show that luminous quasars occur in diverse environments that include ellipticals as bright as the brightest cluster galaxies (2), apparently normal ellipticals (10), apparently normal spirals with H II regions (3), complex systems of gravitationally interacting components (3), and faint surrounding nebulosity (2). The quasar host galaxies are centered on the quasar to the accuracy of our measurements, 0.4 kpc. There are more radio quiet quasars in galaxies that appear to be ellipticals (7) than in spiral hosts (3), contrary to expectations. The average luminosity of the hosts is 2.2 magnitudes brighter than expected for a field galaxy luminosity function. We detect eight companion galaxies within projected distances of 10 kpc from quasar nuclei. Best figs. at http://www.sns.ias.edu/~jnb/Html/preprints.html

We analyze the disk M dwarfs found in 31 new fields observed with the Wide Field Camera (WFC2) on the Hubble Space Telescope (HST) together with the sample previously analyzed from 22 WFC2 fields and 162 prerepair Planetary Camera (PC1) fields. The new observations, which include the 28 high-latitude fields comprising the Large Area Multi-Color Survey (``Groth Strip''), increase the total sample to 337 stars, and more than double the number of late M dwarfs (M_V>13.5) from 23 to 47. The mass function changes slope at M~0.6 Msun, from a near-Salpeter power-law index of \alpha=-1.21 to \alpha=0.44. In both regimes the mass function at the Galactic plane is given by ${d^3 N / d\log M d M_V d V} = 8.1\times 10^{-2}\pc^{-3} ({M / 0.59 M_\odot})^{\alpha}$. The correction for secondaries in binaries changes the low-mass index from \alpha=0.44 to \alpha\sim 0.1. If the Salpeter slope continued to the hydrogen-burning limit, we would expect 500 stars in the last four bins (14.5<M_V<18.5), instead of the 25 actually detected. The explanation of the observed microlensing rate towards the Galactic bulge requires either a substantial population of bulge brown dwarfs or that the disk and bulge mass functions are very different for stars with M~< 0.5 Msun.

To probe the nature of stellar evolution in dense environments, we study the dense core of the globular cluster Messier 13 using pre-refurbishment Planetary Camera-I images obtained with the Hubble Space Telescope. We find 15 blue straggler star candidates and 10 other possible blue stragglers in this region of M13. Their specific frequency is in the range 0.04-0.07, comparable to what is observed near the centers of other dense clusters. A comparison between M13's observed V band stellar luminosity function and a theoretical model (Bergbusch & Vandenberg 1992) for the luminosity function of an old, metal-poor cluster shows that the model predicts too few of the brightest red giants (V~12.5-15) by a factor of two relative to subgiants/turnoff stars (> 6 sigma effect). The radial distributions of the red giants, blue stragglers, and subgiants are consistent with one another. Stars in the blue horizontal branch of M13, however, appear to be centrally depleted relative to other stellar types. We combine data from three dense `King model clusters', M13, M3, and 47 Tuc, and two post core collapse clusters, M30 and M15, and compare the distributions of various stellar types as a function of (r/r_half light) and (r/r_core). The horizontal branch stars in the combined sample appear to be centrally depleted relative to the giants (97% significance) - this depletion is only marginal (< about 1.5 sigma) in each of the clusters taken individually. The blue stragglers in the combined sample are centrally concentrated relative to the giants.

The neutrino fluxes calculated from 14 standard solar models published recently in refereed journals are inconsistent with the results of the 4 pioneering solar neutrino experiments if nothing happens to the neutrinos after they are created in the solar interior. The sound speeds calculated from standard solar models are in excellent agreement with helioseismological measurements of sound speeds. Some statements made by Dar at Neutrino 96 are answered here.

We show that the effects of neutrino oscillations on 8B solar neutrinos are described well by the first two moments (the average and the variance) of the energy distribution of scattered or recoil electrons. For the SuperKamiokande and the Sudbury Neutrino Observatory experiments, the differences between the moments calculated with oscillations and the standard, no-oscillation moments are greater than 3 standard deviations for a significant fraction of the neutrino mass-mixing (Delta m^2, sin^2 2 theta) parameter space.

We report the discovery of a high-ionization broad absorption line system at a redshift of z_abs=0.695 in the spectrum of the z_em=1.052 radio-quiet quasar PG 2302+039. Broad absorption with FWHM from 3,000 to 5,000 km/s is detected from C IV, N V, and O VI in HST Faint Object Spectrograph spectra of the quasar. A narrow line system (FWHM approximately 250 km/s) at z_abs=0.7016 is resolved from the broad blend and includes absorption by Ly-alpha and the C IV, N V, and O VI doublets. No absorption by low-ionization metal species (e.g. Si II and Mg II) is detected in the HST or ground-based spectra for either the broad or the narrow system. The centroids of the broad system lines are displaced by approximately 56,000 km/s to the blue of the quasar's broad emission lines. The reddest extent of the broad line absorption is more than 50,000 km/s from the quasar. The properties of this system are unprecedented whether it is an intervening or an ejected system.

Measurements of the interaction rates of the solar neutrino lines of $^7$Be and $pep$ can be used, independent of solar models, to test whether electron flavor is conserved, to determine survival probabilities of electron-type neutrinos at specific energies, and to test for the existence of sterile neutrinos. We present analytic descriptions of these tests. We also illustrate by numerical simulations, assuming matter-enhanced and vacuum neutrino oscillations, what measurements of solar neutrino lines can teach us about neutrino masses and mixing angles.

We analyze tests of electron flavor conservation that can be performed at the Sudbury Neutrino Observatory (SNO). These tests, which utilize $^8$B solar neutrinos interacting with deuterium, measure: 1) the shape of the recoil electron spectrum in charged-current (CC) interactions (the CC spectrum shape); and 2) the ratio of the number of charged current to neutral current (NC) events (the CC/NC ratio). We determine standard model predictions for the CC spectral shape and for the CC/NC ratio, together with realistic estimates of their errors and the correlations between errors. We consider systematic uncertainties in the standard neutrino spectrum and in the charged-current and neutral current cross-sections, the SNO energy resolution and absolute energy scale, and the SNO detection efficiencies. Assuming that either matter-enhanced or vacuum neutrino oscillations solve the solar neutrino problems, we calculate the confidence levels with which electron flavor non-conservation can be detected using either the CC spectrum shape or the CC/NC ratio, or both. If the SNO detector works as expected, the neutrino oscillation solutions that best-fit the results of the four operating solar neutrino experiments can be distinguished unambiguously from the standard predictions of electron flavor conservation.

The neutrino fluxes calculated from the 14 standard solar models published recently in refereed journals are inconsistent with the results of the 4 pioneering solar neutrino experiments if nothing happens to the neutrinos after they are created in the solar interior. The calculated fluxes and the experimental results are in good agreement if neutrino oscillations occur.

We present a systematic evaluation of the shape of the neutrino energy spectrum produced by beta-decay of $^8$B. We place special emphasis on determining the range of uncertainties permitted by existing laboratory data and theoretical ingredients (such as forbidden and radiative corrections). We review and compare the available experimental data on the $^8$B$(\beta^+){}^8$Be$(2\alpha)$ decay chain. We analyze the theoretical and experimental uncertainties quantitatively. We give a numerical representation of the best-fit (standard-model) neutrino spectrum, as well as two extreme deviations from the standard spectrum that represent the total (experimental and theoretical) effective $\pm3\sigma$ deviations. Solar neutrino experiments that are currently being developed will be able to measure the shape of the $^8$B neutrino spectrum above about 5 MeV. An observed distortion of the $^8$B solar neutrino spectrum outside the range given in the present work could be considered as evidence, at an effective significance level greater than three standard deviations, for physics beyond the standard electroweak model. We use the most recent available experimental data on the Gamow--Teller strengths in the $A=37$ system to calculate the $^8$B neutrino absorption cross section on chlorine: $\sigma_{\rm Cl}=(1.14\pm0.11)\times10^{-42}$~cm$^2$ ($\pm3\sigma$ errors). The chlorine cross section is also given as a function of the neutrino energy. The $^8$B neutrino absorption cross section in gallium is $\sigma_{\rm Ga}=(2.46^{+2.1}_{-1.1})\times10^{-42}$ cm$^2$ ($\pm3\sigma$ errors).

We use a new technique to search for faint red stars in the Hubble Deep Field (HDF) imaged by the Wide Field Camera (WFC2) on the {\it Hubble Space Telescope}. We construct a densely sampled stellar light profile from a set of undersampled images of bright stars in a low-latitude field. Comparison of this stellar light profile to densely sampled profiles of individual objects in the HDF constructed from multiple undersampled dithers allows us to distinguish unambiguously between stars and galaxies to $I=26.3$. We find no stars with $V-I>1.8$ in the outer 90\% of the volume probed. This result places strong and general constraints on the $I$ band luminosity of the constituents of the Galactic dark halo: $$M_I > 15.9 + {5\over 3}\log\biggl(f{0.5 M_\odot\over M}\biggr)\qquad (V-I>1.8),$$ where $M$ is the mass of the objects and $f$ is their density as a fraction of the local halo density, taken to be $\rho_0=9\times 10^{-3}\,M_\odot\,\pc^{-3}$. If the halo is made of white dwarfs, this limit implies that these objects have $M_V\gsim 18.4$ and $V-I\gsim 2.5$. That is, they are $\gsim 2$ magnitudes fainter than the end of the disk white dwarf sequence. Faint red dwarfs account for $<1\%$ of the Galactic dark halo for $M_I<14$, and $<6\%$ for $M_I<15$ at the 95\% confidence level. The density of inter-galactic Local Group stars is at least a factor 3000 smaller than the density of local Population II stars.

We examine the clustering of galaxies around a sample of 20 luminous low redshift (z<0.30) quasars observed with the Wide Field Camera-2 on the Hubble Space Telescope. The HST resolution makes possible galaxy identification brighter than V=23.5 and as close as 2'' to the quasar. We find a significant enhancement of galaxies within a projected separation of < 100 kpc/h of the quasars. If we model the qso/galaxy correlation function as a power law with a slope given by the galaxy/galaxy correlation function, we find that the ratio of the qso/galaxy to galaxy/galaxy correlation functions is $3.8\pm 0.8$. The galaxy counts within r<15 kpc/h of the quasars are too high for the density profile to have an appreciable core radius ( > 100 kpc). Our results reinforce the idea that low redshift quasars are located preferentially in groups of 10-20 galaxies rather than in rich clusters. We see no significant difference in the clustering amplitudes derived from radio-loud and radio-quiet subsamples.

Assuming that MSW neutrino oscillations occur and ignoring all solar physics except for the constraint that nuclear fusion produces the solar luminosity, we show that new solar neutrino experiments are required to rule out empirically the hypothesis that the sun shines via the CNO cycle.

By comparing neutrino fluxes and central temperatures calculated from 1000 detailed numerical solar models, we derive improved scaling laws which show how each of the neutrino fluxes depends upon the central temperature (flux $\propto T^m$); we also estimate uncertainties for the temperature exponents. With the aid of a one-zone model of the sun, we derive expressions for the temperature exponents of the neutrino fluxes. For the most important neutrino fluxes, the exponents calculated with the one-zone model agree to within 20\% or better with the exponents extracted from the detailed numerical models. The one-zone model provides a physical understanding of the temperature dependence of the neutrino fluxes. For the $pp$ neutrino flux, the one-zone model explains the (initially-surprising) dependence of the flux upon a negative power of the temperature and suggests a new functional dependence. This new function makes explicit the strong anti-correlation between the $^7$Be and $pp$ neutrino fluxes. The one-zone model also predicts successfully the average linear relations between neutrino fluxes, but cannot predict the appreciable scatter in a $\Delta \phi_i/\phi_i$ versus $\Delta \phi_j/\phi_j$ diagram.

We argue that the solar g-modes are unlikely to have caused the discrete peaks in the power spectrum of the solar wind flux observed by Thomson et al. (1995). The lower limit to the energy of individual g-modes, using the amplitudes given by Thomson et al., is estimated to be at least 10$^{36}$ erg for low order g-modes; the resulting surface velocity amplitude is at least 50 cm s$^{-1}$, larger than the observational upper limit (5 cm s$^{-1}$). We suggest that the most likely source for the excitation of solar g-modes is turbulent stresses in the convection zone. The surface velocity amplitude of low degree and low order g-modes resulting from this process is estimated to be of order 10$^{-2}$ cm s$^{-1}$. This amplitude is interestingly close to the detection threshold of the SOHO satellite. The long lifetime of g-modes ($\sim 10^6$ years for low order modes) should be helpful in detecting these small amplitude pulsations.

This talk answers a series of questions. Why study solar neutrinos? What does the combined standard model (solar plus electroweak) predict for solar neutrinos? Why are the calculations of neutrino fluxes robust? What are the three solar neutrino problems? What have we learned in the first thirty years of solar neutrino research? For the next decade, what are the most important solvable problems in the physics of solar neutrinos? What are the most important problems in the astrophysics of solar neutrinos?

We describe images of the center of the dense globular cluster M15 (NGC 7078) obtained with the Hubble Space Telescope Wide Field and Planetary Camera 2 (WFPC2). Data taken in the F336W, F439W, and F555W filters (approximately U, B, and V) are used to study the surface density distribution of the $\sim3\times10^4$~stars detected in a 5~arcmin$^2$ region within $r<2'$ (6.7$\,$pc) of the cluster center. The surface density of stars in M15 (after correction for the effects of incompleteness and photometric bias/scatter) is well represented by a power law in radius: $N(r)\sim{r}^{-0.82\pm0.12}$, over the radial range $0.3''$ (0.017$\,$pc) to $6''$. The observed power law is remarkably similar to what is expected if the center of the cluster harbors a massive black hole. Non-parametric estimates of the density profile show a monotonic rise with decreasing radius all the way in to $r=0.3''$, the smallest radius at which the density can be reliably measured; there is no indication that the profile flattens at smaller radii. Any flat core of radius larger than $2''$ (0.11$\,$pc) in the stellar distribution is ruled out at the $>95$\% significance level. The star count profile is consistent with that expected from core-collapse models or with the predicted distribution around a massive (few times $10^3\,M_\odot$) black hole. The close triplet of bright stars, AC~214, is within $0.5''$ ($1.5\sigma$) of the cluster centroid position.

HST images (with WFPC2) of PHL~909\ ($z = 0.171$) and PG~0052$+$251\ ($z = 0.155$) show that these luminous radio-quiet quasars each occur in an apparently normal host galaxy. The host galaxy of PHL~909 is an elliptical galaxy ($\sim$ E4) and the host of PG~0052$+$251 is a spiral ($\sim$~Sb). Both host galaxies are several tenths of a magnitude brighter than $L^*$, the characteristic Schechter luminosity of field galaxies. The images of PHL~909 and PG~0052$+$251, when compared with HST images of other objects in our sample of 20 luminous, small-redshift ($z \leq 0.30$) quasars, show that luminous quasars occur in a variety of environments. The local environments of the luminous quasars range from luminous ellipticals, to apparently normal host galaxies, to complex systems of interacting components, to faint (and as yet undetected) hosts. The bright HII regions of the host galaxy of PG~0052$+$251 provide an opportunity to measure directly the metallicity of the host of a luminous quasar, to establish an upper limit to the mass of the nuclear AGN (i.e., the putative black hole source), and to test stringently the cosmological hypothesisthat the galaxy and the quasar are both at the distance indicated by the quasar redshift.

We present red and blue images of the jet of the quasar 3C273 obtained with the WFPC2 on the Hubble Space Telescope as well as a new radio map made with the MERLIN array. The images are of significantly better quality than that of previous data. The two maps are aligned to an accuracy of 0.020$''$; this accuracy is achieved because both the quasar and the jet are contained in both the radio and optical images. The start of the optical jet is marked by an elongated knot which appears identical at radio and optical wavelengths. Other knots in the optical jet correspond to narrow oblique features within the radio outline. The total width of the smooth emission in the optical jet is~$0.7''$; the FWHM of the optical knots is~$0.3''$. The knots may trace the current location of a narrow, perhaps helical jet lying within the outlines of the older radio cocoon.

We present evidence that clumps of \lya lines are physically associated with about half of the extensive metal-line systems (absorption systems with four or more observed metal-line species) found in this paper, demonstrate that all four \lylimit systems discussed here correspond to extensive metal-line absorption systems, and present an extraordinary pair of extensive metal-line absorption systems within 2000~km/s of each other at $z~=~0.95$ that are probably an early manifestation of large scale structure. These results are obtained using ultraviolet spectra, taken with the higher-resolution gratings of the Faint Object Spectrograph of the Hubble Space Telescope, for four quasars with emission-line redshifts between 1.0 and 1.3. We also determine the evolution of \lya absorption lines at redshifts less than 1.3 by combining the results for 13 smaller redshift quasars discussed in Paper I of this series with the 4 moderate redshift quasars analyzed in the present paper. For $z_{\rm abs}~\leq~1.3$, the density of \lya lines with equivalent widths greater than 0.24~\AA\ is adequately fit by $\left(dN/dz\right) =\left(dN/dz\right)_0 \cdot(1 + z)^\gamma$ with $(dN/dz)_0 = 24.3 \pm 6.6$ \lya lines per unit redshift, and $\gamma = 0.58 \pm 0.50$ (1-$\sigma$ uncertainties). This rate of evolution at low redshifts is less than the evolutionary rate inferred from several different ground-based data samples that pertain to high redshifts. The observed gaseous structures at redshifts of $0.5$ to $1.0$ with velocity dispersions of $6 \times 10^2$ \kms to $1.4 \times 10^3$ \kms (or velocity spans of $1.2 \times 10^3$ \kms to $3 \times 10^3$ \kms) constitute a constraint on cosmological models of structure formation.

Heavy element and helium diffusion are both included for the first time in this series of papers on precise solar models, along with improvements in the input data for nuclear reaction rates, the solar luminosity, the solar age, heavy element abundances, radiative opacities, helium and metal diffusion rates, and neutrino interaction cross sections. The calculated neutrino event rates, including all of the improvements, are $9.3_{-1.4}^{+1.2}$ SNU for the \hbox{$^{37}$Cl} experimentand $137_{-7}^{+8} $ SNU for the \hbox{$^{71}$Ga} experiments. The calculated flux of \hbox{$^{7}$Be} neutrinos is $5.1(1.00_{-0.07}^{+0.06}) \times 10^{9}~{\rm cm^{-2}s^{-1}}$ and the flux of ${\rm ^8B}$ neutrinos is $6.5 (1.00_{-0.17}^{+0.14})\times 10^{6}~{\rm cm^{-2}s^{-1}}$. The present-day surface helium abundance of the model is $Y_s = 0.247$, in agreement with the helioseismological measurement of $Y_s = 0.242 \pm 0.003$ determined by Hernandez and Christensen-Dalsgaard (1994). The computed depth of the convective zone is \hbox{$R = 0.712~ R_{\odot}$} in agreement with the observed value determined from $p$-mode oscillation data of \hbox{$R = 0.713 \pm 0.003~ R_{\odot}$} found by Christensen-Dalsgaard {\it et al.}~(1991). Solar models that do not include diffusion disagree with the helioseismology measurements. Although the present results increase the predicted event rate in the four operating solar neutrino experiments by almost $1\sigma$ (theoretical uncertainty), they only slightly increase the difficulty of explaining the existing experiments with standard physics (i.e., by assuming that nothing happens to the neutrinos after they are created in the center of the sun).

We study a sample of 257 Galactic disk M dwarfs (8<M_V<18.5) found in images obtained using HST. These include 192 stars in 22 fields imaged with the repaired WFC2 with mean limiting mag I=23.7 and 65 stars in 162 fields imaged with the pre-repair Planetary Camera with mean limiting mag V=21.3. We find that the disk luminosity function (LF) drops sharply for M_V>12 (M<0.25 \ms), decreasing by a factor \gsim 3 by M_V~14 (M~0.14\ms). This decrease in the LF is in good agreement with the ground-based photometric study of nearby stars by Stobie et al. (1989), and in mild conflict with the most recent LF measurements based on local parallax stars by Reid et al. (1995). The local LF of the faint Galactic disk stars can be transformed into a local mass function using an empirical mass-M_V relation. The mass function can be represented analytically over the mass range 0.1\ms<M<1.6\ms by \log(\phi)=-1.35-1.34\log(M/\ms)-1.85 [\log(M/\ms)]^2 where \phi is the number density per logarithmic unit of mass. The total column density of M stars is only \Sigma_M=11.8\pm 1.8\ms\pc^{-2}, implying a total `observed' disk column density of \Sigma_\obs~=39\ms\pc^{-2}, lower than previously believed, and also lower than all estimates with which we are familiar of the dynamically inferred mass of the disk. The measured scale length for the M-star disk is 3.0\pm 0.4 kpc. The optical depth to microlensing toward the LMC by the observed stars in the Milky Way disk is \tau~1x10^{-8}, compared to the observed optical depth found in ongoing experiments \tau_\obs~ 10^{-7}. The M-stars show evidence for a population with characteristics intermediate between thin disk and spheroid populations. Approximating what may be a continuum of populations by two separate component, we find characteristic exponential scale heights of ~210 pc and ~740 pc.

HST images (WFC2) of \pks show that this luminous nearby quasar is interacting with diffuse (presumably galactic) material. Two thin wisps that have a total extent of about 20 kpc (for $H_0 = 100$ \kms ${\rm Mpc}^{-1}$ and $\Omega_0 = 1.0$) are observed to approximately surround the quasar. One of the wisps appears to pass through a companion galaxy that is located at a projected distance of three~kpc from the center of the quasar light. The companion galaxy, if located at the distance of \pks, has an intrinsic size and luminosity similar to the Large Magellanic Cloud. A faint extended nebulosity, which is detected over a region of 35 kpc by 50 kpc and is centered about 5 kpc from the quasar nucleus, overlaps the wisps. The immediate environment of \pks is different from the enivronments of the other eight luminous quasars that we have studied previously with HST. If the multiple light components of the HST images are fit to a single de~Vaucouleurs profile, as was done in previous analyses of ground--based data, then the result obtained for the total luminosity of the model galaxy is in agreement with the earlier ground-based studies.

The neutrino deficits observed in four solar neutrino experiments, relative to the theoretical predictions, have led to fresh insights into neutrino and solar physics. Neutrino emission from distant, energetic astronomical systems may form the basis for a new field of astronomy. Additional experiments are needed to test explanations of the solar neutrino deficits and to allow the detection of more distant sources.

Radiative corrections to the electron recoil-energy spectra and to total cross sections are computed for neutrino-electron scattering by solar neutrinos. Radiative corrections change monotonically the electron recoil spectrum for incident \b8 neutrinos, with the relative probability of observing recoil electrons being reduced by about 4 \% at the highest electron energies. For $p-p$ and $\be7$ neutrinos, the recoil spectra are not affected significantly. Total cross sections for solar neutrino-electron scattering are reduced by about 2 \% compared to previously computed values. We also calculate the recoil spectra from $^{13}$N and $^{15}$O neutrinos including radiative corrections.

HST observations of eight intrinsically luminous quasars with redshifts between $0.16$ and~$0.29$ are presented. Seven companion galaxies brighter than $M_V = -16.5$ lie within a projected distance of 25 kpc of the quasars; three of the companions are located closer than $3''$ (6 kpc projected distance) from the quasars, well within the volume that would be enclosed by a typical $L^*$ host galaxy. The observed association of quasars and companion galaxies is statistically significant and may be an important element in the luminous-quasar phenomenon. Evidence for candidate host galaxies is presented for the three most promising cases: PG 1116+215, 3C 273, and PG 1444+407, but additional observations are required before the characteristics of the candidate hosts can be regarded as established. Upper limits are placed on the visual-band brightnesses of representative galactic hosts for all of the quasars. These limits are established by placing galaxy images obtained with HST underneath the quasars and measuring at what faintness level the known galaxies are detected. On average, the HST spirals would have been detected if they were one magnitude fainter than $L^*$ or brighter and the early-type galaxies could have been detected down to a brightness level of about $L^*$, where $L^*$ is the Schechter characteristic luminosity of field galaxies.

We study the proposed solution of the solar neutrino problem which requires a flavor nondiagonal coupling of neutrinos to gravity. We adopt a phenomenological point of view and investigate the consequences of the hypothesis that the neutrino weak interaction eigenstates are linear combinations of the gravitational eigenstates which have slightly different couplings to gravity, $f_1G$ and $f_2G$, $|f_1-f_2| << 1$, corresponding to a difference in red-shift between electron and muon neutrinos, $\Delta z/(1+z) \sim |f_1 - f_2|$. We perform a $\chi^2$ analysis of the latest available solar neutrino data and obtain the allowed regions in the space of the relevant parameters. The existing data rule out most of the parameter space which can be probed in solar neutrino experiments, allowing only $|f_1 - f_2| \sim 3 \times 10^{-14}$ for small values of the mixing angle ($2 \times 10^{-3} \le \sin^2(2\theta_G) \le 10^{-2}$) and $10^{-16} \stackrel{<}{\sim} |f_1 - f_2| \stackrel{<}{\sim}10^{-15}$ for large mixing ($0.6 \le \sin^2(2\theta_G) \le 0.9$). Measurements of the $^8{\rm B}$-neutrino energy spectrum in the SNO and Super-Kamiokande experiments will provide stronger constraints independent of all considerations related to solar models. We show that these measurements will be able to exclude part of the allowed region as well as to distinguish between conventional oscillations and oscillations due to the violation of the equivalence principle.

Strong upper limits are placed on the visual-band brightnesses of galactic hosts for four luminous, radio-quiet quasars with redshifts between $0.16$ and~$0.24$ that were studied with the HST's Wide Field/Planetary Camera-2. Typical upper limits on the luminosities of galactic hosts are about $1.4$ mag fainter than $L^*$ for spirals and about $0.5$ mag fainter than $L^*$ for ellipticals. The galactic hosts of the quasars are more than a magnitude and a half fainter than the median integrated absolute magnitude of Seyfert galaxies. If the detection limits are determined using featureless simulated galaxies instead of observed galaxy images, then the detection limits for spirals are 0.5--1.0 mag less stringent. These results are consistent with the hypothesis that the quasar phenomenon corresponds to the early stages of galaxy formation, before extensive star formation occurs.

Low-ionization nuclear emission-line regions (LINERs), which exist in a large fraction of galaxies, may be the least luminous manifestation of quasar activity. The nature of LINERs has, however, remained controversial because an AGN-like nonstellar continuum source has not been directly observed in them. We report the detection of bright, unresolved (FWHM $\ltorder 0.1''$) point sources of UV ($\sim 2300$ \AA) emission in the nuclei of nine nearby galaxies from a complete sample of 110 nearby galaxies imaged with {\it HST}. Ground-based optical spectroscopy reveals that five of the nuclei are LINERs, three are starburst nuclei, and one is a Seyfert nucleus. The observed UV flux in each of the five LINERs implies an ionizing flux that is sufficient to account for the observed emission lines through photoionization. The detection of a strong UV continuum in the LINERs argues against shock excitation as the source of the observed emission lines, and supports the idea that photoionization excites the lines in at least some objects of this class. Among the Northern-hemisphere galaxies in the sample, 26 are LINERs, of which only the above five LINERs have a detected nuclear UV source. There are no obvious differences in the optical line intensity ratios between the UV-bright and UV-dark LINERs. If all LINERs are photoionized, then the continuum source is unobscured along our line of sight in $5/26\approx 20\%$ of LINERs. Alternatively, spectrally-similar LINERs may be produced by various excitation mechanisms, with photoionization responsible in only about 20\% of the cases. The high angular resolution allows us to set upper limits, typically several parsecs, on the physical size of the compact star-cluster or AGN-type continuum source that is emitting the UV light in these objects.

Ten recently-published solar models give $\7be$ neutrino fluxes that lie within a range of $\pm 10$\% of the average value, a convergence that is independent of uncertainties in the measured laboratory rate of the $\7be(p,\gamma)\8b$ reaction. If nothing happens to solar neutrinos after they are created ({\it a la} standard electroweak theory) and the operating solar neutrino experiments are correct, then the $\7be$ solar neutrino flux must be less than 50\% of the solar model value. At least three of the four existing solar neutrino experiments must be wrong {\it if}: (1) standard electroweak theory is correct, and (2) the true $\7be$ neutrino flux lies within the range predicted by standard solar models.

Predicted rates for solar neutrino experiments that are obtained with a modified stellar evolution code originally developed to study the advanced stages of stellar evolution are shown to be in agreement with other recently-calculated precise solar models to about 2% (0.2 SNU for the chlorine experiment). Different scenarios for pre-main sequence evolution are explored and are found to change the predicted rates for solar neutrino experiments by less than or of order 1%. The influence of the depth of the solar convective zone on the predicted solar neutrino fluxes is established by direct calculation. It is shown that a change in the calculated depth of the convective zone that is five times larger than the quoted helioseismological measurement uncertainty determined by Christensen-Dalsgaard, Gough, & Thompson causes a change in the predicted $^8$B neutrino flux of less than 7% and a change in the $^7$Be neutrino flux of less than 4%. Additionally, it is shown that the radiative opacities near the depth of the convective zone cannot differ from the standard OPAL values by more than about 12% without causing the calculated and measured depths of the convective zone to differ by more than five times the helioseismological measurement uncertainty.

We show that faint red stars do not contribute significantly to the mass budget of the Galaxy or to microlensing statistics. Our results are obtained by analyzing two long exposures of a high-latitude field taken with the Wide Field Camera on the newly repaired {\it Hubble Space Telescope}. Stars are easily distinguished from galaxies essentially to the limiting magnitudes of the images. We find five stars with $2.0<V-I<3.0$ and $I<25.3$ and no stars with $V-I>3.0$. Therefore, main-sequence stars with $M_I>10$ that are above the hydrogen-burning limit in the dark halo or the spheroid contribute $<6%$ of the unseen matter. Faint red disk stars, M-dwarfs, contribute at most $ 15%$ to the mass of the disk. We parameterize the faint end of the cumulative distribution of stars, $\Phi$, as a function of luminosity $L_V$, $d\Phi/d \ln L_V \propto L_V^{-\gamma}$. For spheroid stars, $\gamma<0.32$ over the range $6<M_V<17$, with 98% confidence. The disk luminosity function falls, $\gamma<0$, for $15\lsim M_V\lsim 19$. Faint red stars in the disk or thick disk, and stars with $M_V<16$ in the spheroid contribute $\tau< 10^{-8}$ to the optical depth to microlensing toward the Large Magellanic Cloud.

There are four important facts about solar neutrinos. They are listed in order of importance in this abstract and discussed more in the text of the talk. First, solar neutrinos have been detected in four experiments with approximately the energies and fluxes predicted by the standard solar model, confirming the hypothesis that the energy source for the solar luminosity is the fusion of light elements. Second, the measured event rates are significantl less than the event rates predicted by the combined standard solar and electroweak models in all four experiments. Third, a comparison of the event rates measured in the chlorine experiment (threshold 0.8 Mev) and the neutrino-electron scattering experiment (Kamiokande II, threshold 7.5 MeV) indicates that the deficiency of electron-type neutrinos at the earth is energy dependent, if the rates and the uncertainties in both experiments have been correctly understood. The inference that the deficiency is energy-dependent conflicts with the simplest version of standard electroweak theory. Fourth, experiments are being constructed that have the capabilities to determine conclusively if new neutrino physics is required.

The claim by Dar and Shaviv that they have found a standard model solution to the solar neutrino problem is based upon an incorrect assumption made in extrapolating nuclear cross sections and the selective use of a small fraction of the nuclear physics and of the neutrino data. In addition, five different solar model codes show that the rate obtained for the chlorine experiment using the Dar-Shaviv stated parameters differs by at least $14 \sigma$ from the observed rate.

A precise test of the theory of stellar evolution can be performed by measuring the average difference in energy between the neutrino line produced by ${\rm ^7Be}$ electron capture in the solar interior and the corresponding neutrino line produced in a terrestrial laboratory. This energy shift is calculated to be 1.29~keV (to an accuracy of a few percent) for the dominant ground-state to ground-state transition. The energy shift is approximately equal to the average temperature of the solar core, computed by integrating the temperature over the solar interior with a weighting factor equal to the locally-produced $^7$Be neutrino emission. The characteristic modulation of the ${\rm ^7Be}$ line shape that would be caused by either vacuum neutrino oscillations or by matter-enhanced (MSW) neutrino oscillations is shown to be small. Other frequently-discussed weak interaction solutions to the solar neutrino problem are also not expected to change significantly the line profile. Therefore, a measurement of the energy shift is a measurement of the central temperature distribution of the sun.

The four operating solar neutrino experiments confirm the hypothesis that the energy source for solar luminosity is hydrogen fusion. However, the measured rate for each of the four solar neutrino experiments differs significantly (by factors of 2.0 to 3.5) from the corresponding theoretical prediction that is based upon the standard solar model and the simplest version of the standard electroweak theory. If standard electroweak theory is correct, the energy spectrum for \b8 neutrinos created in the solar interior must be the same (to one part in $10^5$) as the known laboratory \b8 neutrino energy spectrum. Direct comparison of the chlorine and the Kamiokande experiments, both sensitive to \b8 neutrinos, suggests that the discrepancy between theory and observations depends upon neutrino energy, in conflict with standard expectations. Monte Carlo studies with 1000 implementations of the standard solar model confirm that the chlorine and the Kamiokande experiments cannot be reconciled unless new weak interaction physics changes the shape of the \b8 neutrino energy spectrum. The results of the two gallium solar neutrino experiments strengthen the conclusion that new physics is required and help determine a relatively small allowed region for the MSW neutrino parameters.

The vacuum-polarization corrections to rates for nuclear-fusion reactions in the $pp$ chain and in the CNO cycle are calculated. For the reactions of particular importance to the solar-neutrino problem, the $\Hethree$, $\Hefour$, $\Be$, and $\Ni$ reactions, we find the magnitude of the effect to be less than 2\%. The effect of VP on all the other reaction rates is expected to be of a similar order of magnitude. We discuss how these results affect the predicted fluxes of solar neutrinos.

A precise test of the theory of stellar evolution can be performed by measuring the difference in average energy between the neutrino line produced by ${\rm ^7Be}$ electron capture in the solar interior and the corresponding neutrino line produced in a terrestrial laboratory. The high temperatures in the center of the sun broaden the line asymmetrically, FWHM = 1.6~keV, and cause an average energy shift of 1.3~keV. The width of the $^7$Be neutrino line should be taken into account in calculations of vacuum neutrino oscillations.

Physics beyond the simplest version of the standard electroweak model is required to reconcile the results of the chlorine and the Kamiokande solar neutrino experiments. None of the 1000 solar models in a full Monte Carlo simulation is consistent with the results of the chlorine or the Kamiokande experiments. Even if the solar models are forced articficially to have a ${}^8 B$ neutrino flux in agreement with the Kamiokande experiment, none of the fudged models agrees with the chlorine observations. This comparison shows that consistency of the chlorine and Kamiokande experiments requires some physical process that changes the shape of the ${}^8 B$ neutrino energy spectrum. The GALLEX and SAGE experiments, which currently have large statistical uncertainties, differ from the predictions of the standard solar model by $2 \sigma$ and $3 \sigma$, respectively. The possibility that the neutrino experiments are incorrect is briefly discussed.

We re-evaluate the matrix element for the proton-proton reaction which is important for stellar-evolution calculations and for the solar-neutrino problem. We self-consistently determine the effect of vacuum polarization on the matrix element by first correcting the low-energy scattering data to account for vacuum polarization. We then calculate the proton-proton wave function by integrating the Schrodinger equation with vacuum polarization included. We use improved data for proton-proton scattering and for the deuteron wave function. We evaluate the uncertainties that are due to experimental error and estimate those that are due to theoretical inadequacies. We estimate the theoretical uncertainty by using six different deuteron potentials and five different proton-proton potentials. Vacuum polarization decreases the calculated value by $0.6_{-0.4}^{+0.1}$\%. The complete result is $\Lambda^2=6.92\times(1\pm0.002^{+0.014}_{-0.009})$ where the first uncertainty is due to experimental errors and the second uncertainty is due to theoretical uncertainties. Our value of $\Lambda^2$ is 2\% smaller than the value obtained in 1969 by Bahcall and May. The improved calculations of the rate of the $pp$ reaction described here increase slightly the predicted event rates for the chlorine and the Kamiokande solar-neutrino experiments.

We analyze the ultraviolet (UV) emission line and continuum properties of five low-redshift active galactic nuclei (four luminous quasars: PKS~0405$-$123, H1821+643, PG~0953+414, and 3C273, and one bright Seyfert 1 galaxy: Mrk~205). The HST spectra have higher signal-to-noise ratios (typically $\sim 60$ per resolution element) and spectral resolution ($R = 1300$) than all previously- published UV spectra used to study the emission characteristics of active galactic nuclei. We include in the analysis ground-based optical spectra covering \hb\ and the narrow [O III] $\lambda\lambda$4959,5007 doublet. The following new results are obtained: \lyb/\lya=0.03$-$0.12 for the four quasars, which is the first accurate measurement of the long-predicted \lyb\ intensity in QSOs. The cores of \lya\ and C~IV are symmetric to an accuracy of better than 2.5% within about 2000 km s$^{-1}$ of the line peak. This high degree of symmetry of \lya\ argues against models in which the broad line cloud velocity field has a significant radial component. The observed smoothness of the \lya\ and C~IV line profiles requires at least $\sim 10^4$ individual clouds if bulk velocity is the only line-broadening mechanism. The overall similarity of the \lya\ and C IV $\lambda$1549 profiles rules out models for the broad line region (BLR) with a radial distribution of virialized....

We study the diffusion of helium and other heavy elements in the solar interior by solving exactly the set of flow equations developed by Burgers for a multi-component fluid, including the residual heat-flow terms. No approximation is made concerning the relative concentrations and no restriction is placed on the number of elements considered. We give improved diffusion velocities for hydrogen, helium, oxygen and iron, in the analytic form derived previously by Bahcall and Loeb. These expressions for the diffusion velocities are simple to program in stellar evolution codes and are expected to be accurate to $\sim 15\%$. Our complete treatment of element diffusion can be directly incorporated in a standard stellar evolution code by means of an exportable subroutine, but, for convenience, we also give simple analytical fits to our numerical results.

Contents: 1. The Conference 2. Solar Neutrinos -A. Operating Experiments -B. New Funded Experiments -C. Are Solar Neutrino Fluxes Time Dependent? -D. Have $pp$ Neutrinos Been Detected? -E. Is Directionality Required for Astronomical Observations? -F. Is There a Solar Neutrino Problem? -G. What Have We Learned? 3. Supernova Neutrinos 4. Atmospheric Neutrinos 5. High-Energy Neutrinos 6. Our Field Is Flourishing

None of the 1000 solar models in a full Monte Carlo simulation is consistent with the results of the chlorine or the Kamiokande experiments. Even if the solar models are forced artifically to have a \b8 neutrino flux in agreeement with the Kamiokande experiment, none of the fudged models agrees with the chlorine observations. The GALLEX and SAGE experiments, which currently have large statistical uncertainties, differ from the predictions of the standard solar model by $2 \sigma$ and $3 \sigma$, respectively.

If the published event rates of the chlorine and Kamiokande solar neutrino experiments are correct, then the energy spectrum of neutrinos produced by the decay of $^8$B in the sun must be different from the energy spectrum determined from laboratory nuclear physics measurements. This change in the energy spectrum requires physics beyond the standard electroweak model. In addition, the GALLEX and SAGE experiments, which currently have large statistical uncertainties, differ from the predictions of the standard solar model by $2 \sigma$ and $3 \sigma$, respectively.

Recent results on standard solar models are reviewed. I shall summarize briefly three of the themes that I stressed at the Neutrino '92 Conference: 1) Different solar model codes give the same answers when the same input data are used; 2) Improved calculations of standard solar models include helium diffusion, the Livermore radiative opacity, the meteoritic iron abundance, and a variety of other corrections; and 3) There are a few basic rules that should be followed in using standard solar models. At the Neutrino '92 Conference, I reviewed in more detail the recent work on standard solar models by Marc Pinsonneault and myself. This work has by now appeared in print (Rev. Mod. Phys. 64, 885, 1992, hereafter Paper I, and ApJ Letters, 69, 717, 1992, Paper II). Therefore, there is no need for me to repeat the details here.