The groundbreaking discoveries of gravitational waves from binary black-hole mergers and, most recently, coalescing neutron stars started a new era of Multi-Messenger Astrophysics and revolutionized our understanding of the Cosmos. Machine learning techniques such as artificial neural networks are already transforming many technological fields and have also proven successful in gravitational-wave astrophysics for detection and characterization of gravitational-wave signals from binary black holes. Here we use a deep-learning approach to rapidly identify transient gravitational-wave signals from binary neutron star mergers in noisy time series representative of typical gravitational-wave detector data. Specifically, we show that a deep convolution neural network trained on 100,000 data samples can rapidly identify binary neutron star gravitational-wave signals and distinguish them from noise and signals from merging black hole binaries. These results demonstrate the potential of artificial neural networks for real-time detection of gravitational-wave signals from binary neutron star mergers, which is critical for a prompt follow-up and detailed observation of the electromagnetic and astro-particle counterparts accompanying these important transients.

Determining the Equation of State (EOS) of dense neutron-rich nuclear matter is a shared goal of both nuclear physics and astrophysics. Except possible phase transitions, the density dependence of nuclear symmetry \esym is the most uncertain part of the EOS of neutron-rich nucleonic matter especially at supra-saturation densities. Much progresses have been made in recent years in predicting the symmetry energy and understanding why it is still very uncertain using various microscopic nuclear many-body theories and phenomenological models. Simultaneously, significant progresses have also been made in probing the symmetry energy in both terrestrial nuclear laboratories and astrophysical observatories. In light of the GW170817 event as well as ongoing or planned nuclear experiments and astrophysical observations probing the EOS of dense neutron-rich matter, we review recent progresses and identify new challenges to the best knowledge we have on several selected topics critical for understanding astrophysical effects of the nuclear symmetry energy.

Applying an equation of state (EOS) with its symmetric nuclear matter (SNM) contribution and low-density symmetry energy $E_{sym}(\rho)$ constrained by heavy-ion reaction data, we calculate the dimensionless tidal deformability $\Lambda$ of neutron stars in coalescing binary systems. Corresponding to the partially constrained EOS that previously predicted a radius of 11.5 km $\leq R_{1.4} \leq$ 13.6 km for canonical neutron-star configurations, $\Lambda$ is found to be in the range of 292 $\leq\Lambda_{1.4}\leq$ 680, consistent with the very recent observation of the GW170817 event. We investigate the effect of the high-density behavior of $E_{sym}(\rho)$ on the tidal properties of neutron stars and find that while $\Lambda$ depends strongly on the details of the symmetry energy, different trends of $E_{sym}(\rho)$ lead to very similar values of $\Lambda$. In particular, the transition from stiff/soft to soft/stiff $E_{sym}(\rho)$ could yield the same $\Lambda$. Thus, measuring $\Lambda$ alone may not determine completely the density dependence of the symmetry energy. Coherent analyses of the dense neutron-rich nuclear matter EOS underlying both nuclear laboratory experiments and astrophysical observations are therefore necessary to break this degeneracy and determine precisely the details of the $E_{sym}(\rho)$.

The recent direct detection of gravitational waves (GWs) from binary black hole mergers (2016, Phys. Rev. Lett. 116, no. 6, 061102; no. 24, 241103) opens up an entirely new non-electromagnetic window into the Universe making it possible to probe physics that has been hidden or dark to electromagnetic observations. In addition to cataclysmic events involving black holes, GWs can be triggered by physical processes and systems involving neutron stars. Properties of neutron stars are largely determined by the equation of state (EOS) of neutron-rich matter, which is the major ingredient in calculating the stellar structure and properties of related phenomena, such as gravitational wave emission from elliptically deformed pulsars and neutron star binaries. Although the EOS of neutron-rich matter is still rather uncertain mainly due to the poorly known density dependence of nuclear symmetry energy at high densities, significant progress has been made recently in constraining the symmetry energy using data from terrestrial nuclear laboratories. These constraints could provide useful information on the limits of GWs expected from neutron stars. Here after briefly reviewing our previous work on constraining gravitational radiation from elliptically deformed pulsars with terrestrial nuclear laboratory data in light of the recent gravitational wave detection, we estimate the maximum gravitational wave strain amplitude, using an optimistic value for the breaking strain of the neutron star crust, for 15 pulsars at distances 0.16 kpc to 0.91 kpc from Earth, and find it to be in the range of $\sim[0.2-31.1]\times 10^{-24}$, depending on the details of the EOS used to compute the neutron star properties. Implications are discussed.

Significant progress has been made in recent years in constraining the density dependence of nuclear symmetry energy using terrestrial nuclear laboratory data. Around and below the nuclear matter saturation density, the experimental constraints start to merge in a relatively narrow region. At supra-saturation densities, there are, however, still large uncertainties. After summarizing the latest experimental constraints on the density dependence of nuclear symmetry energy, we highlight a few recent studies examining imprints of nuclear symmetry energy on the binding energy, energy release during hadron-quark phase transitions as well as the $w$-mode frequency and damping time of gravitational wave emission of neutron stars.

Pulsars are among the most mysterious astrophysical objects in the Universe and are believed to be rotating neutron stars formed in supernova explosions. They are unique testing grounds of dense matter theories and gravitational physics and also provide links among nuclear physics, particle physics and General Relativity. Neutron stars may exhibit some of the most extreme and exotic characteristics that could not be found elsewhere in the Universe. Their properties are largely determined by the equation of state (EOS) of neutron-rich matter, which is the chief ingredient in calculating neutron star structure and properties of related phenomena, such as gravitational wave emission from deformed pulsars. Presently, the EOS of neutron-rich matter is still very uncertain mainly due to the poorly known density dependence of the nuclear symmetry energy especially at supra-saturation densities. Nevertheless, significant progress has been made recently in constraining the density dependence of the nuclear symmetry energy mostly at sub-saturation densities using terrestrial nuclear reactions. While there are still some uncertainties especially at supra-saturation densities, these constraints could provide useful information on the limits of the global properties of pulsars and the gravitational waves to be expected from them. Here we review our recent work on constraining properties of pulsars and gravitational radiation with data from terrestrial nuclear laboratories.

The eigen-frequencies of the axial w-modes of oscillating neutron stars are studied using the continued fraction method with an Equation of State (EOS) partially constrained by the recent terrestrial nuclear laboratory data. It is shown that the density dependence of the nuclear symmetry energy $E_{sym}(\rho)$ affects significantly both the frequencies and the damping times of these modes. Besides confirming the previously found universal behavior of the mass-scaled eigen-frequencies as functions of the compactness of neutron stars, we explored several alternative universal scaling functions. Moreover, the $w_{II}$-mode is found to exist only for neutron stars having a compactness of $M/R\geq 0.1078$ independent of the EOS used.

Heavy-ion reactions especially those induced by radioactive beams provide useful information about the density dependence of the nuclear symmetry energy, thus the Equation of State of neutron-rich nuclear matter, relevant for many astrophysical studies. The latest developments in constraining the symmetry energy at both sub- and supra-saturation densities from analyses of the isopsin diffusion and the $\pi^-/\pi^+$ ratio in heavy-ion collisions using the IBUU04 transport model are discussed. Astrophysical ramifications of the partially constrained symmetry energy on properties of neutron star crusts, gravitational waves emitted by deformed pulsars and the w-mode oscillations of neutron stars are presented briefly.

Gravitational waves are tiny disturbances in space-time and are a fundamental, although not yet directly confirmed, prediction of General Relativity. Rapidly rotating neutron stars are one of the possible sources of gravitational radiation dependent upon pulsar's rotational frequency, details of the equation of state of stellar matter, and distance to detector. Applying an equation of state with symmetry energy constrained by recent nuclear laboratory data, we set an upper limit on the strain-amplitude of gravitational waves emitted by rapidly rotating neutron stars.

Gravitational radiation is a fundamental prediction of General Relativity. Elliptically deformed pulsars are among the possible sources emitting gravitational waves (GWs) with a strain-amplitude dependent upon the star's quadrupole moment, rotational frequency, and distance from the detector. We show that the gravitational wave strain amplitude $h_0$ depends strongly on the equation of state of neutron-rich stellar matter. Applying an equation of state with symmetry energy constrained by recent nuclear laboratory data, we set an upper limit on the strain-amplitude of GWs produced by elliptically deformed pulsars. Depending on details of the EOS, for several millisecond pulsars at distances $0.18kpc$ to $0.35kpc$ from Earth, the {\it maximal} $h_0$ is found to be in the range of $\sim[0.4-1.5]\times 10^{-24}$. This prediction serves as the first {\it direct} nuclear constraint on the gravitational radiation. Its implications are discussed.

Recent analyses of several isospin effects in heavy-ion reactions have allowed us to constrain the density dependence of nuclear symmetry energy at sub-saturation densities within a narrow range. Combined with constraints on the Equation of State (EOS) of symmetric nuclear matter obtained previously from analyzing the elliptic flow in relativistic heavy-ion collisions, the EOS of neutron-rich nuclear matter is thus partially constrained. Here we report effects of the partially constrained EOS of neutron-rich nuclear matter on the mass-radius correlation, moment of inertia, elliptical deformation and gravitational radiation of (rapidly) rotating neutron stars.

Properties and structure of neutron stars are determined by the equation of state (EOS) of neutron-rich stellar matter. While the collective flow and particle production in relativistic heavy-ion collisions have constrained tightly the EOS of symmetric nuclear matter up to about five times the normal nuclear matter density, the more recent experimental data on isospin-diffusion and isoscaling in heavy-ion collisions at intermediate energies have constrained considerably the density dependence of the nuclear symmetry energy at subsaturation densities. Although there are still many uncertainties and challenges to pin down completely the EOS of neutron-rich nuclear matter, the heavy-ion reaction experiments in terrestrial laboratories have limited the EOS of neutron-rich nuclear matter in a range much narrower than that spanned by various EOSs currently used in astrophysical studies in the literature. These nuclear physics constraints could thus provide more reliable information about properties of neutron stars. Within well established formalisms using the nuclear constrained EOSs we study the momenta of inertia of neutron stars. We put the special emphasis on the component A of the extremely relativistic double neutron star system PSR J0737-3039. Its moment of inertia is found to be between 1.30 and 1.63 $(\times10^{45}g$ $cm^2)$. Moreover, the transition density at the crust-core boundary is shown to be in the narrow range of $\rho_t=[0.091-0.093](fm^{-3})$.

Properties, structure, and thermal evolution of neutron stars are determined by the equation of state of stellar matter. Recent data on isospin-diffusion and isoscaling in heavy-ion collisions at intermediate energies as well as the size of neutron skin in $^{208}Pb$ have constrained considerably the density dependence of the nuclear symmetry energy and, in turn, the equation of state of neutron-rich nucleonic matter. These constraints could provide useful information about the global properties of rapidly rotating neutron stars. Models of rapidly rotating neutron stars are constructed applying several nucleonic equations of state. Particular emphasis is placed on configurations rotating rigidly at $716$ and $1122Hz$. The range of allowed hydrostatic equilibrium solutions is determined and tested for stability. The effect of rotation on the internal composition and thermal properties of neutron stars is also examined. At a given rotational frequency, each equation of state yields a range of possible neutron stars configurations restricted by the Keplerian (mass-shedding) limit, corresponding to the maximal circumferential radius, and the limit due to the onset of instabilities with respect to axial-symmetric perturbations, corresponding to the minimal equatorial radius of a stable neutron star model. We show that the mass of a neutron star rotating uniformly at $1122Hz$ is between $1.7$ and $2.1M_{\sun}$. Central stellar density and proton fraction decrease with increasing rotational frequency with respect to static models, and depending on the exact stellar mass and angular velocity, can drop below the Direct Urca threshold thus closing the fast cooling channel.

Testing the constancy of the gravitational constant G has been a longstanding fundamental question in natural science. As first suggested by Jofr\'{e}, Reisenegger and Fern\'{a}ndez [1], Dirac's hypothesis of a decreasing gravitational constant $G$ with time due to the expansion of the Universe would induce changes in the composition of neutron stars, causing dissipation and internal heating. Eventually, neutron stars reach their quasi-stationary states where cooling due to neutrino and photon emissions balances the internal heating. The correlation of surface temperatures and radii of some old neutron stars may thus carry useful information about the changing rate of G. Using the density dependence of the nuclear symmetry energy constrained by recent terrestrial laboratory data on isospin diffusion in heavy-ion reactions at intermediate energies and the size of neutron skin in $^{208}Pb$ within the gravitochemical heating formalism, we obtain an upper limit of the relative changing rate of $|\dot{G}/G|\le4\times 10^{-12}yr^{-1}$ consistent with the best available estimates in the literature.

Heavy-ion reactions provide a unique means to investigate the equation of state (EOS) of neutron-rich nuclear matter, especially the density dependence of the nuclear symmetry energy $E_{sym}(\rho)$. The latter plays an important role in understanding many key issues in both nuclear physics and astrophysics. Recent analyses of heavy-ion reactions have already put a stringent constraint on the $E_{sym}(\rho)$ around the saturation density. This subsequently allowed us to constrain significantly the radii and cooling mechanisms of neutron stars as well as the possible changing rate of the gravitational constant G.

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.

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

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.

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

We analyze the latest published solar neutrino data assuming an arbitrary neutrino oscillation/conversion mechanism suppresses the electron neutrino flux from the Sun independent of energy. For oscillations/transitions into active (sterile) neutrinos such mechanisms are ruled out at 99.96 (99.9997) % C.L. assuming the standard solar model by Bahcall and Pinnsoneault '95 correctly predicts all solar neutrino fluxes within their estimated uncertainties. Even if one allows for $^8{\rm B}$ and $^7{\rm Be}$ solar neutrino fluxes that are vastly different from the ones in contemporary standard solar models these mechanisms are strongly disfavored by the data.

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.

On the basis of present solar neutrino observations and relaxing the constraints from solar models it is possible that most (or nearly all) of the flux of electron neutrinos detected comes from electron capture in $^7{Be}$. These solutions arise from neutrino oscillations in which $\nu_e$ $-$ $\nu_\tau$ mixing suppresses high energy $\nu_e$ and $\nu_e - \nu_\mu$ mixing suppresses low energy $\nu_e$ as qualitatively suggested from some SO(10) grand unified models. The importance of future observations is emphasized.

The most popular solutions of the solar neutrino problem which assume massive neutrinos predict varying with time signals in the solar neutrino detectors. The variations represent a distinctive feature of these solutions and do not depend on the theoretically predicted total flux of neutrinos from the Sun. The sensitivity of SuperKamiokande and SNO to day-night variations due to regeneration in the Earth is estimated.

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 study the stability of the two--neutrino MSW solution of the solar neutrino problem, corresponding to solar $\nu_e$ transitions into sterile neutrino, $\nu_e \rightarrow \nu_s$, with respect to changes of the total fluxes of $^{8}$B and $^{7}$Be neutrinos, $\Phi_{{\rm B}}$ and $\Phi_{{\rm Be}}$. For any value of $\Phi_{{\rm Be}}$ from the interval $0.7\Phi^{{\rm BP}}_{{\rm Be}}\leq \Phi_{{\rm Be}} \leq 1.3\Phi^{{\rm BP}}_{{\rm Be}}$ (for $\Phi_{{\rm Be}} = 0.7\Phi^{{\rm BP}}_{{\rm Be}}$) the $\nu_e \rightarrow \nu_{s}$ MSW transitions provide at 95\% C.L. a description of the existing solar neutrino data for $0.40~(0.39)~\Phi^{{\rm BP}}_{{\rm B}} \ltap \Phi_{{\rm B}} \ltap 36~(40)\Phi^{{\rm BP}}_{{\rm B}}$, $\Phi^{{\rm BP}}_{{\rm B}}$ and $\Phi^{{\rm BP}}_{{\rm Be}}$ being the fluxes in the solar model of Bahcall--Pinsonneault from 1992. The corresponding allowed regions of values of the parameters \dm2 and \sn2, characterizing the solar neutrino transitions, are derived. The physical implications of the found MSW $\nu_e \rightarrow \nu_s$ solutions for the future solar neutrino experiments are considered as well.

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.

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.

There are large uncertainties in the predictions of the boron neutrino flux from the Sun which cannot be considered as being of purely statistical origin. We treat the magnitude of this flux, $\Phi_B$, as a parameter to be found from experiment. The properties of the MSW solution to the solar neutrino problem for different values of $\Phi_B$ are studied. Present data give the bounds: $0.38 < \Phi_B/ \Phi_B^0 < 3.1$ (2$\sigma$), where $\Phi_B^0 \equiv 5.7 \cdot 10^6$ cm$^{-2}$s$^{-1}$ is the flux in the reference SSM. The variations of the flux in this interval enlarge the allowed region of mixing angles: $\sin^2 2\theta = 8\cdot 10^{-4} \div 2\cdot 10^{-2}$ (small mixing solutions) and $\sin^2 2\theta = 0.2 \div 0.85$ (large mixing solution). If the value of the original boron neutrino flux is about that measured by Kamiokande, a consistent description of the data is achieved for $\sin^2 2\theta \sim (0.8 \div 2)\cdot 10^{-3}$ (``very small mixing solution"). The solution is characterized by a strong suppression of the beryllium neutrino line, a weak distortion of the high energy part of the boron neutrino spectrum and a value of the double ratio $(CC/NC)^{exp}/(CC/NC)^{SSM}$ at $E > 5$ MeV close to 1. We comment on the possibility to measure the neutrino parameters and the original boron neutrino flux in future experiments.