In the past decade, the properties of annihilating dark matter models were examined by various kinds of data, including the data of gamma rays, radio waves, X-ray, positrons, electrons, antiprotons and neutrinos. In particular, most of the studies focus on the data of our Galaxy, nearby galaxies (e.g. M31 galaxy) or nearby galaxy clusters (e.g. Fornax cluster). In this article, we examine the archival radio continuum spectral data of a relatively high-redshift galaxy cluster (A697 cluster) to constrain the properties of annihilating dark matter. We find that leptophilic annihilation channels ($e^+e^-$, $\mu^+\mu^-$ and $\tau^+\tau^-$) can give very good fits to the radio continuum spectrum of the A697 cluster.

Flavor transitions in supernova neutrinos are yet to be determined. We present a method to probe whether or not the Mikheyev-Smirnov-Wolfenstein effects occur as SN neutrinos propagate outward from the SN core by investigating time evolutions of neutrino event rates for different flavors in different kinds of detectors. As the MSW effect occurs, the $\nu_e$ flux swaps with the $\nu_x$ flux, which represents any one of $\nu_{\mu}$, $\nu_{\tau}$, $\bar{\nu}_{\mu}$, and $\bar{\nu}_{\tau}$ flux, either fully or partially depending on the neutrino mass hierarchy. During the neutronization burst, the $\nu_e$ emission evolves in a much different shape from the emissions of $\bar{\nu}_e$ and $\nu_x$ while the latter two evolve in a similar pattern. Meanwhile, the luminosity of the the $\nu_e$ emission is much larger than those of the $\bar{\nu}_e$ and $\nu_x$ emissions while the latter two are roughly equal. As a consequence, the time-evolution pattern of the $\nu_e{\rm Ar}$ event rates in the absence of the MSW effect will be much different from that in the occurrence of the MSW effect, in either mass hierarchy. With the simulated SN neutrino emissions, the $\nu_e{\rm Ar}$ and inverse beta decay event rates are evaluated. The ratios of the two cumulative event rates are calculated for different progenitor masses up to $100~{\rm ms}$. We show that the time evolutions of this cumulative ratios can effectively determine whether MSW effects really occur for SN neutrinos or not.

Fast radio bursts (FRBs) are bright, millisecond-duration radio transients originating from extragalactic distances. Their origin is unknown. Some FRB sources emit repeat bursts, ruling out cataclysmic origins for those events. Despite searches for periodicity in repeat burst arrival times on time scales from milliseconds to many days, these bursts have hitherto been observed to appear sporadically, and though clustered, without a regular pattern. Here we report the detection of a $16.35\pm0.15$ day periodicity (or possibly a higher-frequency alias of that periodicity) from a repeating FRB 180916.J0158+65 detected by the Canadian Hydrogen Intensity Mapping Experiment Fast Radio Burst Project (CHIME/FRB). In 38 bursts recorded from September 16th, 2018 through February 4th, 2020, we find that all bursts arrive in a 5-day phase window, and 50% of the bursts arrive in a 0.6-day phase window. Our results suggest a mechanism for periodic modulation either of the burst emission itself, or through external amplification or absorption, and disfavour models invoking purely sporadic processes.

Magnetars are highly magnetized young neutron stars that occasionally produce enormous bursts and flares of X-rays and gamma-rays. Of the approximately thirty magnetars currently known in our Galaxy and Magellanic Clouds, five have exhibited transient radio pulsations. Fast radio bursts (FRBs) are millisecond-duration bursts of radio waves arriving from cosmological distances. Some have been seen to repeat. A leading model for repeating FRBs is that they are extragalactic magnetars, powered by their intense magnetic fields. However, a challenge to this model has been that FRBs must have radio luminosities many orders of magnitude larger than those seen from known Galactic magnetars. Here we report the detection of an extremely intense radio burst from the Galactic magnetar SGR 1935+2154 using the Canadian Hydrogen Intensity Mapping Experiment (CHIME) FRB project. The fluence of this two-component bright radio burst and the estimated distance to SGR 1935+2154 together imply a 400-800 MHz burst energy of $\sim 3 \times 10^{34}$ erg, which is three orders of magnitude brighter than those of any radio-emitting magnetar detected thus far. Such a burst coming from a nearby galaxy would be indistinguishable from a typical FRB. This event thus bridges a large fraction of the radio energy gap between the population of Galactic magnetars and FRBs, strongly supporting the notion that magnetars are the origin of at least some FRBs.

We report on the detection of seven bursts from the periodically active, repeating fast radio burst (FRB) source FRB 180916.J0158+65 in the 300-400-MHz frequency range with the Green Bank Telescope (GBT). Emission in multiple bursts is visible down to the bottom of the GBT band, suggesting that the cutoff frequency (if it exists) for FRB emission is lower than 300 MHz. Observations were conducted during predicted periods of activity of the source, and had simultaneous coverage with the Low Frequency Array (LOFAR) and the FRB backend on the Canadian Hydrogen Intensity Mapping Experiment (CHIME) telescope. We find that one of the GBT-detected bursts has potentially associated emission in the CHIME band (400-800 MHz) but we detect no bursts in the LOFAR band (110-190 MHz), placing a limit of $\alpha > -1.0$ on the spectral index of broadband emission from the source. We also find that emission from the source is severely band-limited with burst bandwidths as low as $\sim$40 MHz. In addition, we place the strictest constraint on observable scattering of the source, $<$ 1.7 ms, at 350 MHz, suggesting that the circumburst environment does not have strong scattering properties. Additionally, knowing that the circumburst environment is optically thin to free-free absorption at 300 MHz, we find evidence against the association of a hyper-compact HII region or a young supernova remnant (age $<$ 50 yr) with the source.

We report on simultaneous radio and X-ray observations of the repeating fast radio burst source FRB 180916.J0158+65 using the Canadian Hydrogen Intensity Mapping Experiment (CHIME), Effelsberg, and Deep Space Network (DSS-14 and DSS-63) radio telescopes and the Chandra X-ray Observatory. During 33 ks of Chandra observations, we detect no radio bursts in overlapping Effelsberg or Deep Space Network observations and a single radio burst during CHIME/FRB source transits. We detect no X-ray events in excess of the background during the Chandra observations. These non-detections imply a 5-$\sigma$ limit of $<5\times10^{-10}$ erg cm$^{-2}$ for the 0.5--10 keV fluence of prompt emission at the time of the radio burst and $1.3\times10^{-9}$ erg cm$^{-2}$ at any time during the Chandra observations at the position of FRB 180916.J0158+65. Given the host-galaxy redshift of FRB 180916.J0158+65 ($z\sim0.034$), these correspond to energy limits of $<1.6\times10^{45}$ erg and $<4\times10^{45}$ erg, respectively. We also place a 5-$\sigma$ limit of $<8\times10^{-15}$ erg s$^{-1}$ cm$^{-2}$ on the 0.5--10\,keV absorbed flux of a persistent source at the location of FRB 180916.J0158+65. This corresponds to a luminosity limit of $<2\times10^{40}$ erg s$^{-1}$. Using Fermi/GBM data we search for prompt gamma-ray emission at the time of radio bursts from FRB 180916.J0158+65 and find no significant bursts, placing a limit of $4\times10^{-9}$ erg cm$^{-2}$ on the 10--100 keV fluence. We also search Fermi/LAT data for periodic modulation of the gamma-ray brightness at the 16.35-day period of radio-burst activity and detect no significant modulation. We compare these deep limits to the predictions of various fast radio burst models, but conclude that similar X-ray constraints on a closer fast radio burst source would be needed to strongly constrain theory.

We report on the discovery and analysis of bursts from nine new repeating fast radio burst (FRB) sources found using the Canadian Hydrogen Intensity Mapping Experiment (CHIME) telescope. These sources span a dispersion measure (DM) range of 195 to 1380 pc cm$^{-3}$. We detect two bursts from three of the new sources, three bursts from four of the new sources, four bursts from one new source, and five bursts from one new source. We determine sky coordinates of all sources with uncertainties of $\sim$10$^\prime$. We detect Faraday rotation measures for two sources, with values $-20(1)$ and $-499.8(7)$ rad m$^{-2}$, that are substantially lower than the RM derived from bursts emitted by FRB 121102. We find that the DM distribution of our events, combined with the nine other repeaters discovered by CHIME/FRB, is indistinguishable from that of thus far non-repeating CHIME/FRB events. However, as previously reported, the burst widths appear statistically significantly larger than the thus far non-repeating CHIME/FRB events, further supporting the notion of inherently different emission mechanisms and/or local environments. These results are consistent with previous work, though are now derived from 18 repeating sources discovered by CHIME/FRB during its first year of operation. We identify candidate galaxies that may contain FRB 190303.J1353+48 (DM = 222.4 pc cm$^{-3}$).

In the past few years, some studies claimed that annihilating dark matter with mass $\sim 10-100$ GeV can explain the GeV gamma-ray excess in our Galaxy. However, recent analyses of the Fermi-LAT and radio observational data rule out the possibility of the thermal relic annihilating dark matter with mass $m \le 100$ GeV for some popular annihilation channels. By using the new observed radio data of the Andromeda galaxy, we rule out the existence of $\sim 100-300$ GeV thermal relic annihilating dark matter for ten annihilation channels. The lower limits of annihilating dark matter mass are improved to larger than 330 GeV for the most conservative case, which is a few times larger than the current best constraints. Moreover, these limits strongly disfavor the benchmark model of weakly interacting massive particle (WIMP) produced through the thermal freeze-out mechanism.

In the past few decades, observations indicated that an unexplained high production rate of positrons (the strong 511 keV line) exists in the Milky Way center. By using the fact that a large amount of high density gas used to exist near the Milky Way center million years ago, we model the rate of positrons produced due to dark matter annihilation. We consider the effect of adiabatic contraction of dark matter density due to the supermassive black hole at the Milky Way center and perform a detailed calculation to constrain the possible annihilation channel and dark matter mass range. We find that only three annihilation channels ($\mu^+\mu^-$, $4e$ and $4 \mu$) can provide the required positron production rate and satisfy the stringent constraint of gamma-ray observations. In particular, the constrained mass range for the $\mu^+\mu^-$ channel is $m \approx 80-100$ GeV, which is close to the mass range obtained for the dark matter interpretation of the GeV gamma-ray and positron excess. In other words, the proposed scenario can simultaneously provide the required positron production rate to explain the 511 keV emission, the positron excess and the GeV gamma-ray excess in our Milky Way, and it is compatible with the density spike due to adiabatic growth model of the supermassive black hole.

We study the equations of motion of the massive and massless particles in the Schwarzschild geometry of general relativity by using the Laplace-Adomian Decomposition Method, which proved to be extremely successful in obtaining series solutions to a wide range of strongly nonlinear differential and integral equations. After introducing a general formalism for the derivation of the equations of motion in arbitrary spherically symmetric static geometries, and of the general mathematical formalism of the Laplace-Adomian Decomposition Method, we obtain the series solution of the geodesics equation in the Schwarzschild geometry. The truncated series solution, containing only five terms, can reproduce the exact numerical solution with a high precision. In the first order of approximation we reobtain the standard expression for the perihelion precession. We study in detail the bending angle of light by compact objects in several orders of approximation. The extension of this approach to more general geometries than the Schwarzschild one is also briefly discussed.

Photons from distant astronomical sources can be used as a classical source of randomness to improve fundamental tests of quantum nonlocality, wave-particle duality, and local realism through Bell's inequality and delayed-choice quantum eraser tests inspired by Wheeler's cosmic-scale Mach-Zehnder interferometer gedankenexperiment. Such sources of random numbers may also be useful for information-theoretic applications such as key distribution for quantum cryptography. Building on the design of an "astronomical random-number generator" developed for the recent "cosmic Bell" experiment [Handsteiner et al., Phys. Rev. Lett. 118, 060401 (2017)], in this paper we report on the design and characterization of a device that, with 20-nanosecond latency, outputs a bit based on whether the wavelength of an incoming photon is greater than or less than 700 nm. Using the one-meter telescope at the Jet Propulsion Laboratory (JPL) Table Mountain Observatory, we generated random bits from astronomical photons in both color channels from 50 stars of varying color and magnitude, and from 12 quasars with redshifts up to $z = 3.9$. With stars, we achieved bit rates of $\sim 1 \times 10^6$ Hz / m$^2$, limited by saturation for our single-photon detectors, and with quasars of magnitudes between 12.9 and 16, we achieved rates between $\sim 10^2$ and $2 \times 10^3$ Hz /m$^2$. For bright quasars, the resulting bitstreams exhibit sufficiently low amounts of statistical predictability as quantified by the mutual information. In addition, a sufficiently high fraction of bits generated are of true astronomical origin in order to address both the locality and freedom-of-choice loopholes when used to set the measurement settings in a test of the Bell-CHSH inequality.

If dark matter is composed of massive bosons, a Bose-Einstein Condensation process must have occurred during the cosmological evolution. Therefore galactic dark matter may be in a form of a condensate, characterized by a strong self-interaction. We consider the effects of rotation on the Bose-Einstein Condensate dark matter halos, and we investigate how rotation might influence their astrophysical properties. In order to describe the condensate we use the Gross-Pitaevskii equation, and the Thomas-Fermi approximation, which predicts a polytropic equation of state with polytropic index $n=1$. By assuming a rigid body rotation for the halo, with the use of the hydrodynamic representation of the Gross-Pitaevskii equation we obtain the basic equation describing the density distribution of the rotating condensate. We obtain the general solutions for the condensed dark matter density, and we derive the general representations for the mass distribution, boundary (radius), potential energy, velocity dispersion, tangential velocity and for the logarithmic density and velocity slopes, respectively. Explicit expressions for the radius, mass, and tangential velocity are obtained in the first order of approximation, under the assumption of slow rotation. In order to compare our results with the observations we fit the theoretical expressions of the tangential velocity of massive test particles moving in rotating Bose-Einstein Condensate dark halos with the data of 12 dwarf galaxies, and the Milky Way, respectively.

The Weak Equivalence Principle states that the geodesics of a test particle in a gravitational field are independent of the particle's constitution. To constrain violations of the Weak Equivalence Principle, we use the one-meter telescope at Table Mountain Observatory near Los Angeles to monitor the relative arrival times of pulses from the Crab Pulsar in the optical ($\lambda \approx 585$ nm) and near-infrared ($\lambda \approx 814$ nm) using an instrument which detects single photons with nanosecond-timing resolution in those two bands. The infrared pulse arrives slightly before the visible pulse. Our three analysis methods give delays with statistical errors of $\Delta t_{obs} = 7.41 \pm 0.58$, $0.4 \pm 3.6$, and $7.35 \pm 4.48$ microseconds (at most 1/4000 of the pulsar period). We attribute this discrepancy to systematic error from the fact that the visible and infrared pulses have slightly different shapes. Whether this delay emerges from the pulsar, is caused by passing through wavelength-dependent media, or is caused by a violation of the equivalence principle, unless there is a fine-tuned cancellation among these, we set the first upper limit on the differential post-Newtonian parameter at these wavelengths of $\Delta \gamma < 1.07 \times 10^{-10}~(3\sigma)$. This result falls in an unexplored region of parameter space and complements existing limits on equivalence-principle violation from fast radio bursts, gamma ray bursts, as well as previous limits from the Crab.

Recently, some very tight constraints of annihilating dark matter have been obtained from gamma-ray data of the Milky Way and Milky Way dwarf spheroidal satellite galaxies. In this article, we report that there are two excellent galaxy clusters (A2877 and Fornax) which can provide interesting constraints for annihilating dark matter. The lower limits of the dark matter mass for the thermal relic annihilation cross section are 25 GeV, 6 GeV, 130 GeV and 100 GeV respectively for the $e^+e^-$, $\mu^+\mu^-$, $\tau^+\tau^-$ and $b\bar{b}$ channels. For some configuration of our working assumptions, our results improve the Fermi-LAT upper limits of annihilation cross sections by a factor of $1.3-1.8$ for wide ranges of dark matter mass for $e^+e^-$, $\mu^+\mu^-$ and $b\bar{b}$ channels, and a factor of $1.2-1.8$ for $\tau^+\tau^-$ channel with dark matter mass $\le 100$ GeV. These limits basically rule out most of the existing popular dark matter interpretation of the GeV excess in the Milky Way.

Bell's theorem states that some predictions of quantum mechanics cannot be reproduced by a local-realist theory. That conflict is expressed by Bell's inequality, which is usually derived under the assumption that there are no statistical correlations between the choices of measurement settings and anything else that can causally affect the measurement outcomes. In previous experiments, this "freedom of choice" was addressed by ensuring that selection of measurement settings via conventional "quantum random number generators" was space-like separated from the entangled particle creation. This, however, left open the possibility that an unknown cause affected both the setting choices and measurement outcomes as recently as mere microseconds before each experimental trial. Here we report on a new experimental test of Bell's inequality that, for the first time, uses distant astronomical sources as "cosmic setting generators." In our tests with polarization-entangled photons, measurement settings were chosen using real-time observations of Milky Way stars while simultaneously ensuring locality. Assuming fair sampling for all detected photons, and that each stellar photon's color was set at emission, we observe statistically significant $\gtrsim 7.31 \sigma$ and $\gtrsim 11.93 \sigma$ violations of Bell's inequality with estimated $p$-values of $ \lesssim 1.8 \times 10^{-13}$ and $\lesssim 4.0 \times 10^{-33}$, respectively, thereby pushing back by $\sim$600 years the most recent time by which any local-realist influences could have engineered the observed Bell violation.

Saturn's moon Enceladus offers a unique opportunity in the search for life and habitable environments beyond Earth, a key theme of the National Research Council's 2013-2022 Decadal Survey. A plume of water vapor and ice spews from Enceladus's south polar region. Cassini data suggest that this plume, sourced by a liquid reservoir beneath the moon's icy crust, contain organics, salts, and water-rock interaction derivatives. Thus, the ingredients for life as we know it-- liquid water, chemistry, and energy sources-- are available in Enceladus's subsurface ocean. We have only to sample the plumes to investigate this hidden ocean environment. We present a New Frontiers class, solar-powered Enceladus orbiter that would take advantage of this opportunity, Testing the Habitability of Enceladus's Ocean (THEO). Developed by the 2015 Jet Propulsion Laboratory Planetary Science Summer School student participants under the guidance of TeamX, this mission concept includes remote sensing and in situ analyses with a mass spectrometer, a sub-mm radiometer-spectrometer, a camera, and two magnetometers. These instruments were selected to address four key questions for ascertaining the habitability of Enceladus's ocean within the context of the moon's geological activity: (1) How are the plumes and ocean connected? (2) Are the abiotic conditions of the ocean suitable for habitability? (3) How stable is the ocean environment? (4) Is there evidence of biological processes? By taking advantage of the opportunity Enceladus's plumes offer, THEO represents a viable, solar-powered option for exploring a potentially habitable ocean world of the outer solar system.

We consider a description of the stochastic oscillations of the general relativistic accretion disks around compact astrophysical objects interacting with their external medium based on a generalized Langevin equation with colored noise, which accounts for the general memory and retarded effects of the frictional force, and on the fluctuation-dissipation theorem. The presence of the memory effects influences the response of the disk to external random interactions, and modifies the dynamical behavior of the disk, as well as the energy dissipation processes. The generalized Langevin equation of the motion of the disk in the vertical direction is studied numerically, and the vertical displacements, velocities and luminosities of the stochastically perturbed disks are explicitly obtained for both the Schwarzschild and the Kerr cases. The Power Spectral Distribution (PSD) of the disk luminosity is also obtained. As a possible astrophysical application of the formalism we investigate the possibility that the Intra Day Variability (IDV) of the Active Galactic Nuclei (AGN) may be due to the stochastic disk instabilities. The perturbations due to colored/nontrivially correlated noise induce a complicated disk dynamics, which could explain some astrophysical observational features related to disk variability.

We consider a description of the stochastic oscillations of the general relativistic accretion disks around compact astrophysical objects based on the generalized Langevin equation, which accounts for the general retarded effects of the frictional force, and on the fluctuation-dissipation theorems. The vertical displacements, velocities and luminosities of the stochastically perturbed disks are explicitly obtained for both the Schwarzschild and the Kerr cases. The Power Spectral Distribution of the luminosity it is also obtained, and it is shown that it has non-standard values. The theoretical predictions of the model are compared with the observational data for the luminosity time variation of the BL Lac S5 0716+714 object.

Fractals are a basic tool to phenomenologically describe natural objects having a high degree of temporal or spatial variability. From a physical point of view the fractal properties of natural systems can also be interpreted by using the standard formalism of thermodynamical fluctuations. In the present paper we introduce and analyze the fractal dimension of the Intra - Day- Variability (IDV) light-curves of the BL Lac objects, in the optical, radio and X-ray bands, respectively. A general description of the fluctuation spectrum of these systems based on general thermodynamical principles is also proposed. Based on the general fractal properties of a given physical system, we also introduce the predictability index for the IDV light curves. We have explicitly determined the fractal dimension for the R-band observations of five blazars, as well as for the radio band observations of the compact extragalactic radio source J 1128+5925, and of several X-ray sources. Our results show that the fractal dimension of the optical and X-ray observations indicates an almost pure "Brownian noise" (random walk) spectrum, with a very low predictability index, while in the radio band the predictability index is much higher. We have also studied the spectral properties of the IDV light curves, and we have shown that their spectral index is very closely correlated with the corresponding fractal dimension.

In this paper, we introduce a simplified model for explaining the observations of the optical intraday variability (IDV) of the BL Lac Objects. We assume that the source of the IDV are the stochastic oscillations of an accretion disk around a supermassive black hole. The Stochastic Fluctuations on the vertical direction of the accretion disk are described by using a Langevin type equation with a damping term and a random, white noise type force. Furthermore, the preliminary numerical simulation results are presented, which are based on the numerical analysis of the Langevin stochastic differential equation.

It is shown that, in the context of split supersymmetry, a simple model with a single complex scalar field can produce chaotic inflation and generate the observed amount of baryon asymmetry via the Affleck-Dine mechanism. While the inflaton quantum fluctuations give rise to curvature perturbation, we show that quantum fluctuations of the phase of the scalar field can produce baryonic isocurvature perturbation. Combining with constraints from WMAP data, all parameters in the model can be determined to within a narrow range.

It was recently conjectured by Shabad and Usov that there exists in QED a maximum magnetic field of 10^{42} G, above which the magnetized vacuum becomes unstable. Using a nonperturbative analysis that consistently incorporates the effective electron mass and the screening effect in a strong magnetic field, we show that the conjectured phenomenon of positronium collapse never takes place. Thus, there does not exist a maximum magnetic field in QED and the magnetized vacuum is stable for all values of the magnetic field.

The gauge independence of the dynamical fermion mass generated through chiral symmetry breaking in QED in a strong, constant external magnetic field is critically examined. We show that the bare vertex approximation, in which the vertex corrections are ignored, is a consistent truncation of the Schwinger-Dyson equations in the lowest Landau level approximation. The dynamical fermion mass, obtained as the solution of the truncated Schwinger-Dyson equations evaluated on the fermion mass shell, is shown to be manifestly gauge independent. By establishing a direct correspondence between the truncated Schwinger-Dyson equations and the 2PI (two-particle-irreducible) effective action truncated at the lowest nontrivial order in the loop expansion as well as in the 1/N_f expansion (N_f is the number of fermion flavors), we argue that in a strong magnetic field the dynamical fermion mass can be reliably calculated in the bare vertex approximation.

The gauge independence of the dynamical fermion mass generated through chiral symmetry breaking in QED in a strong, constant external magnetic field is critically examined. We present a (first, to the best of our knowledge) consistent truncation of the Schwinger-Dyson equations in the lowest Landau level approximation. We demonstrate that the dynamical fermion mass, obtained as the solution of the truncated Schwinger-Dyson equations evaluated on the fermion mass shell, is manifestly gauge independent.

Possible effects due to nonequilibrium dynamics in the Affleck-Dine mechanism of baryogenesis are examined. Using the closed-time-path formalism, the quantum fluctuation and the back reaction of the Affleck-Dine scalar field are incorporated self-consistently into the dynamical equations of the system. It is found that such nonequilibrium effects can significantly affect the amount of baryon asymmetry that can be generated. In particular, it is possible to generate the observed baryon asymmetry with suitable initial conditions. The methodology described in this paper as well as some of the results obtained are quite general, and can be applied to any complex scalar field in a cosmological background.

We consider the effect of a long range, flavor changing tensor interaction of possible gravitational origin. Neutrino mixing experiments provide the most sensitive probe to date for such forces---testing the equivalence principle at levels below $10^{-20}$. Here we justify and generalize a formalism for describing such effects. The constraints from neutrino mixing experiments on gravitationally induced mixing are calculated. Our detailed analysis of the atmospheric neutrino data confirms a remarkable result: the atmospheric neutrino data implies the same size force as does the solar neutrino data. Additional tests of this suggestive result are discussed.

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.