Analyses of neutrino measurements acquired by the Super-Kamiokande Neutrino Observatory (SK, for the time interval 1996 - 2001) and of radon decay measurements acquired by the Geological Survey of Israel (GSI, for the time interval 2007 - 2017) yield remarkably consistent detections of the same oscillation: frequency 9.43 +/- 0.04 year-1 (SK), 9.44 +/- 0.04 year-1 (GSI); amplitude 6.8 +/- 1.7 % (SK), 7.0 +/- 1.0 % (GSI); phase 124 +/- 15 deg. (SK), 124 +/- 9 deg. (GSI). We briefly discuss possible hypotheses that may be relevant to this experimental result.

Gamma radiation associated with radon decay exhibits variations in both time of year and time of day. The annual oscillation has its maximum value in June, suggestive of a galactic influence. Measurements made at midnight show strong evidence of an influence of solar rotation, but measurements made at noon do not. We find several pairs of oscillations with frequencies separated by 1 cycle per year that is suggestive of an influence of rotation that is oblique with respect to the normal to the ecliptic. We suggest that beta decays may be stimulated by neutrinos and that the decay products tend to travel in the same direction as the neutrinos. We estimate the relevant cross section.

An experiment carried out at the Brookhaven National Laboratory from February 1982 to December 1989 acquired 364 measurements of the beta-decay rates of a sample of 36Cl and of a sample of 32Si. The experimenters reported finding small periodic annual deviations of the data points from an exponential decay - of uncertain origin. We here analyze this dataset by power spectrum analysis and by forming spectrograms and phasegrams. We confirm the occurrence of annual oscillations but we also find evidence of oscillations in a band of frequencies appropriate for the internal rotation of the Sun. Both datasets show clear evidence of a transient oscillation with a frequency of 12.7 cycles per year that falls in the range of rotational frequencies for the solar radiative zone. We repeat these analyses for 358 neutrino measurements acquired by Super-Kamiokande over the interval May 1986 to August 2001. Spectrogram analysis yields a strong and steady oscillation at about 9.5 cycles per year and an intermittent oscillation with frequency in the range 12.5 - 12.7 cycles per year. We attribute the former to rotation of the solar core and the latter to rotation in the radiative zone. Since the flux of neutrinos (8B neutrinos) responsible for the Super-Kamiokande measurements is known, we are able to estimate the cross sections for the beta-decay oscillations at 12.7 cycles per year. These estimates are found to be 10-21.6 cm-2 for 36Cl and 10-18.4 cm-2 for 32Si. We suggest that the beta-decay process is influenced by neutrinos, and that the solar neutrino flux is modulated by magnetic field in the deep solar interior by Resonant Spin Flavor Precession.

Analysis of solar radius measurements acquired by the Michelson Doppler Imager on the SOHO spacecraft supports previously reported evidence of solar internal r-mode oscillations in Mt Wilson radius data and in nuclear-decay data acquired at the Lomonosov Moscow State University. The frequencies of these oscillations are compatible with oscillations in a putative inner tachocline that separates a slowly rotating core from the radiative envelope.

Claims by the DAMA (DArk MAtter) collaboration to have detected an annually varying signal consistent with models of dark matter appear to be at variance with results from other dark-matter searches. To further understand the DAMA results, we have carried out an independent analysis of DAMA data reconstructed from published figures. In addition to reexamining the Lomb-Scargle and chi-square analyses previously carried out by the DAMA collaboration, we carry out two new likelihood analyses and a new chi-square analysis, focusing attention on the treatment of experimental errors and binning. We confirm the existence of an annual oscillation, with a maximum in early June, but at a lower significance level than previously reported.

Power-spectrum analyses of the decay rates of certain nuclides reveal (at very high confidence levels) an annual oscillation and periodicities that may be attributed to solar rotation and to solar r-mode oscillations. A comparison of spectrograms (time-frequency displays) formed from decay data and from solar neutrino data reveals a common periodicity with frequency 12.5 year-1, which is indicative of the solar radiative zone. We propose that the neutrino flux is modulated by the solar magnetic field (via Resonant Spin Flavor Precession) in that region, and we estimate the force and the torque that could be exerted on a nuclide by the solar neutrino flux.

We report a reanalysis of data on the measured decay rate ratio $^{22}$Na/$^{44}$Ti which were originally published by Norman et al., and interpreted as supporting the conventional hypothesis that nuclear decay rates are constant and not affected by outside influences. We find upon a more detailed analysis of both the amplitude and the phase of the Norman data that they actually favor the presence of an annual variation in $^{22}$Na/$^{44}$Ti, albeit weakly. Moreover, this conclusion holds for a broad range of parameters describing the amplitude and phase of an annual sinusoidal variation in these data. The results from this and related analyses underscore the growing importance of phase considerations in understanding the possible influence of the Sun on nuclear decays. Our conclusions with respect to the phase of the Norman data are consistent with independent analyses of solar neutrino data obtained at Super-Kamiokande-I and the Sudbury Neutrino Observatory (SNO).

This article presents a comparative analysis of solar activity data, Mt Wilson diameter data, Super-Kamiokande solar neutrino data, and nuclear decay data acquired at the Lomonosov Moscow State University (LMSU). We propose that salient periodicities in all of these datasets may be attributed to r-mode oscillations. Periodicities in the solar activity data and in Super-Kamiokande solar neutrino data may be attributed to r-mode oscillations in the known tachocline, with normalized radius in the range 0.66 to 0.74, where the sidereal rotation rate is in the range 13.7 to 14.6 year-1. We propose that periodicities in the Mt Wilson and LMSU data may be attributed to similar r-mode oscillations where the sidereal rotation rate is approximately which we attribute to a hypothetical inner tachocline separating a slowly rotating core from the radiative zone. We also discuss the possible role of the RSFP (Resonant Spin Flavor Precession) process, which leads to estimates of the neutrino magnetic moment and of the magnetic field strength in or near the solar core.

Additional experimental evidence is presented in support of the recent hypothesis that a possible solar influence could explain fluctuations observed in the measured decay rates of some isotopes. These data were obtained during routine weekly calibrations of an instrument used for radiological safety at The Ohio State University Research Reactor using Cl-36. The detector system used was based on a Geiger-Mueller gas detector, which is a robust detector system with very low susceptibility to environmental changes. A clear annual variation is evident in the data, with a maximum relative count rate observed in January/February, and a minimum relative count rate observed in July/August, for seven successive years from July 2005 to June 2011. This annual variation is not likely to have arisen from changes in the detector surroundings, as we show here.

Experiments showing a seasonal variation of the nuclear decay rates of a number of different nuclei, and decay anomalies apparently related to solar flares and solar rotation, have suggested that the Sun may somehow be influencing nuclear decay processes. Recently, Cooper searched for such an effect in $^{238}$Pu nuclei contained in the radioisotope thermoelectric generators (RTGs) on board the Cassini spacecraft. In this paper we modify and extend Cooper's analysis to obtain constraints on anomalous decays of $^{238}$Pu over a wider range of models, but these limits cannot be applied to other nuclei if the anomaly is composition-dependent. We also show that it may require very high sensitivity for terrestrial experiments to discriminate among some models if such a decay anomaly exists, motivating the consideration of future spacecraft experiments which would require less precision.

This article presents an analysis of about 29,000 measurements of gamma radiation associated with the decay of radon in a sealed container at the Geological Survey of Israel (GSI) Laboratory in Jerusalem between 28 January 2007 and 10 May 2010. These measurements exhibit strong variations in time of year and time of day, which may be due in part to environmental influences. However, time-series analysis reveals a number of periodicities, including two at approximately 11.2 year$^{-1}$ and 12.5 year$^{-1}$. We have previously found these oscillations in nuclear-decay data acquired at the Brookhaven National Laboratory (BNL) and at the Physikalisch-Technische Bundesanstalt (PTB), and we have suggested that these oscillations are attributable to some form of solar radiation that has its origin in the deep solar interior. A curious property of the GSI data is that the annual oscillation is much stronger in daytime data than in nighttime data, but the opposite is true for all other oscillations. This may be a systematic effect but, if it is not, this property should help narrow the theoretical options for the mechanism responsible for decay-rate variability.

The recent discovery of PeV electrons from the Crab nebula, produced on rapid time scales of one day or less with a sharply peaked gamma-ray spectrum without hard X-rays, challenges traditional models of diffusive shock acceleration followed by synchrotron radiation. Here we outline an accleration model involving a DC electric field parallel to the magnetic field in a twisted toroidal field around the pulsar. Sudden developments of resistivity in localized regions of the twisted field are thought to drive the particle acceleration, up to PeV energies, resulting in flares. This model can reproduce the observed time scales of $T \approx 1$ day, the peak photon energies of $U_{\Phi,rr} \approx 1$ MeV, maximum electron energies of $U_{e,rr} \approx 1$ PeV, and luminosities of $L \approx 10^{36}$ erg s$^{-1}$.

This article presents a power-spectrum analysis of 2,350 measurements of the $^{90}$Sr/$^{90}$Y decay process acquired over the interval 4 August 2002 to 6 February 2009 at the Lomonosov Moscow State University (LMSU). As we have found for other long sequences of decay measurements, the power spectrum is dominated by a very strong annual oscillation. However, we also find a set of low-frequency peaks, ranging from 0.26 year$^{-1}$ to 3.98 year$^{-1}$, which are very similar to an array of peaks in a power spectrum formed from Mt Wilson solar diameter measurements. The Mt Wilson measurements have been interpreted in terms of r-mode oscillations in a region where the sidereal rotation frequency is 12.08 year$^{-1}$. We find that the LMSU measurements may also be attributed to the same type of r-mode oscillations in a solar region with the same sidereal rotation frequency. We propose that these oscillations occur in an inner tachocline that separates the radiative zone from a more slowly rotating solar core.

Recent analyses of datasets acquired at the Brookhaven National Laboratory and at the Physikalisch-Technische Bundesanstalt both show evidence of pronounced annual variations, suggestive of a solar influence. However, the phases of decay-rate maxima do not correspond precisely to the phase of minimum Sun-Earth distance, as might then be expected. We here examine the hypothesis that decay rates are influenced by an unknown solar radiation, but that the intensity of the radiation is influenced not only by the variation in Sun-Earth distance, but also by a possible North-South asymmetry in the solar emission mechanism. We find that this can lead to phases of decay-rate maxima in the range 0 to 0.183 or 0.683 to 1 (September 6 to March 8) but that, according to this hypothesis, phases in the range 0.183 to 0.683 (March 8 to September 6) are "forbidden." We find that phases of the three datasets here analyzed fall in the allowed range.

Since the 1930s, and with very few exceptions, it has been assumed that the process of radioactive decay is a random process, unaffected by the environment in which the decaying nucleus resides. There have been instances within the past few decades, however, where changes in the chemical environment or physical environment brought about small changes in the decay rates. But even in light of these instances, decaying nuclei that were undisturbed or un-"pressured" were thought to behave in the expected random way, subject to the normal decay probabilities which are specific to each nuclide. Moreover, any "non-random" behavior was assumed automatically to be the fault of the detection systems, the environment surrounding the detectors, or changes in the background radiation to which the detector was exposed. Recently, however, evidence has emerged from a variety of sources, including measurements taken by independent groups at Brookhaven National Laboratory, Physikalisch-Technische Bundesanstalt, and Purdue University, that indicate there may in fact be an influence that is altering nuclear decay rates, albeit at levels on the order of $10^{-3}$. In this paper, we will discuss some of these results, and examine the evidence pointing to the conclusion that the intrinsic decay process is being affected by a solar influence.

Unexplained annual variations in nuclear decay rates have been reported in recent years by a number of groups. We show that data from these experiments exhibit not only variations in time related to Earth-Sun distance, but also periodicities attributable to solar rotation. Additionally, anomalous decay rates coincident in time with a series of solar flares in December 2006 also point to a solar influence on nuclear decay rates. This influence could arise from some flavor of solar neutrinos, or through some other objects we call "neutrellos" which behave in some ways like neutrinos. The indication that neutrinos or neutrellos must interact weakly in the Sun implies that we may be able to use data on time-varying nuclear decay rates to probe the interior of the Sun, a technique which we may call "helioradiology".

Recent analyses of nuclear decay data show evidence of variations suggestive of a solar influence. Analyses of datasets acquired at the Brookhaven National Laboratory (BNL) and at the Physikalisch-Technische Bundesanstalt (PTB) both show evidence of an annual periodicity and of periodicities with sidereal frequencies in the neighborhood of 12.25 year^{-1} (at a significance level that we have estimated to be 10^{-17}). It is notable that this implied rotation rate is lower than that attributed to the solar radiative zone, suggestive of a slowly rotating solar core. This leads us to hypothesize that there may be an "inner tachocline" separating the core from the radiative zone, analogous to the "outer tachocline" that separates the radiative zone from the convection zone. The Rieger periodicity (which has a period of about 154 days, corresponding to a frequency of 2.37 year^{-1}) may be attributed to an r-mode oscillation with spherical-harmonic indices l=3, m=1, located in the outer tachocline. This suggests that we may test the hypothesis of a solar influence on nuclear decay rates by searching BNL and PTB data for evidence of a "Rieger-like" r-mode oscillation, with l=3, m=1, in the inner tachocline. The appropriate search band for such an oscillation is estimated to be 2.00-2.28 year^{-1}. We find, in both datasets, strong evidence of a periodicity at 2.11 year^{-1}. We estimate that the probability of obtaining these results by chance is 10^{-12}.

Evidence for an anomalous annual periodicity in certain nuclear decay data has led to speculation concerning a possible solar influence on nuclear processes. We have recently analyzed data concerning the decay rates of Cl-36 and Si-32, acquired at the Brookhaven National Laboratory (BNL), to search for evidence that might be indicative of a process involving solar rotation. Smoothing of the power spectrum by weighted-running-mean analysis leads to a significant peak at frequency 11.18/yr, which is lower than the equatorial synodic rotation rates of the convection and radiative zones. This article concerns measurements of the decay rates of Ra-226 acquired at the Physikalisch-Technische Bundesanstalt (PTB) in Germany. We find that a similar (but not identical) analysis yields a significant peak in the PTB dataset at frequency 11.21/yr, and a peak in the BNL dataset at 11.25/yr. The change in the BNL result is not significant since the uncertainties in the BNL and PTB analyses are estimated to be 0.13/yr and 0.07/yr, respectively. Combining the two running means by forming the joint power statistic leads to a highly significant peak at frequency 11.23/yr. We comment briefly on the possible implications of these results for solar physics and for particle physics.

This article presents a power-spectrum analysis of 39,024 measurements of the solar diameter made at the Mount Wilson Observatory from 1968.670 to 1997.965. This power spectrum contains a number of very strong peaks. We find that eight of these peaks agree closely with the frequencies of r-mode oscillations for a region of the Sun where the sidereal rotation frequency is 12.08 year$^{-1}$. We estimate that there is less than one chance in ten to the sixth power of finding this pattern by chance.

Recent reports of periodic fluctuations in nuclear decay data of certain isotopes have led to the suggestion that nuclear decay rates are being influenced by the Sun, perhaps via neutrinos. Here we present evidence for the existence of an additional periodicity that appears to be related to the Rieger periodicity well known in solar physics.

Evidence for an anomalous annual periodicity in certain nuclear decay data has led to speculation concerning a possible solar influence on nuclear processes. As a test of this hypothesis, we here search for evidence in decay data that might be indicative of a process involving solar rotation, focusing on data for 32Si and 36Cl decay rates acquired at the Brookhaven National Laboratory. Examination of the power spectrum over a range of frequencies (10 - 15 year^-1) appropriate for solar synodic rotation rates reveals several periodicities, the most prominent being one at 11.18 year^-1 with power 20.76. We evaluate the significance of this peak in terms of the false-alarm probability, by means of the shuffle test, and also by means of a new test (the "shake" test) that involves small random time displacements. The last two tests indicate that the peak at 11.18 year^-1 would arise by chance only once out of about 10^7 trials. Since there are several peaks in the search band, we also investigate the running mean of the power spectrum, and identify a major peak at 11.93 year^-1 with peak running-mean power 4.08. Application of the shuffle test and the shake test indicates that there is less than one chance in 10^11, and one chance in 10^15, respectively, finding by chance a value as large as 4.08.

The term "false-alarm probability" denotes the probability that at least one out of M independent power values in a prescribed search band of a power spectrum computed from a white-noise time series is expected to be as large as or larger than a given value. The usual formula is based on the assumption that powers are distributed exponentially, as one expects for power measurements of normally distributed random noise. However, in practice one typically examines peaks in an over-sampled power spectrum. It is therefore more appropriate to compare the strength of a particular peak with the distribution of peaks in over-sampled power spectra derived from normally distributed random noise. We show that this leads to a formula for the false-alarm probability that is more conservative than the familiar formula. We also show how to combine these results with a Bayesian method for estimating the probability of the null hypothesis (that there is no oscillation in the time series), and we discuss as an example the application of these procedures to Super-Kamiokande solar neutrino data.

The histogram formed from published capture-rate measurements for the GALLEX solar neutrino experiment is bimodal, showing two distinct peaks. On the other hand, the histogram formed from published measurements derived from the similar GNO experiment is unimodal, showing only one peak. However, the two experiments differ in run durations: GALLEX runs are either three weeks or four weeks (approximately) in duration, whereas GNO runs are all about four weeks in duration. When we form 3-week and 4-week subsets of the GALLEX data, we find that the relevant histograms are unimodal. The upper peak arises mainly from the 3-week runs, and the lower peak from the 4-week runs. The 4-week subset of the GALLEX dataset is found to be similar to the GNO dataset. A recent re-analysis of GALLEX data leads to a unimodal histogram.

The usual procedure for estimating the significance of a peak in a power spectrum is to calculate the probability of obtaining that value or a larger value by chance (known as the "p-value"), on the assumption that the time series contains only noise - typically that the measurements are derived from random samplings of a Gaussian distribution. We really need to know the probability that the time series is - or is not - compatible with the null hypothesis that the measurements are derived from noise. This probability can be calculated by Bayesian analysis, but this requires one to specify and evaluate a second hypothesis, that the time series does contain a contribution other than noise. We approach the problem of identifying this function in two ways. We first propose three simple conditions that it seems reasonable to impose on this function, and show that these conditions may be satisfied by a simple function with one free parameter. We then define two different ways of combining information derived from two independent power estimates. We find that this consistency condition may be satisfied, to good approximation, by a special case of the previously proposed likelihood function. We find that the resulting significance estimates are considerably more conservative than those usually associated with the p-values. As two examples, we apply the new procedure to two recent analyses of solar neutrino data: (a) power spectrum analysis of Super-Kamiokande data, and (b) the combined analysis of radiochemical neutrino data and irradiance data.

Recent coordinated power-spectrum analyses of radiochemical solar neutrino data and the solar irradiance have revealed a highly significant, high-Q common modulation at 11.85 yr-1. Since the stability of this frequency points to an explanation in terms of rotation, this result may be attributable to non-spherically-symmetric nuclear burning in a solar core with sidereal rotation frequency 12.85 yr-1. The variability of the amplitude (on a timescale of years) suggests that the relevant nuclear burning is variable as well as asymmetric. Recent analysis of Super-Kamiokande solar neutrino data has revealed r-mode-type modulations with frequencies corresponding to a region with sidereal rotation frequency 13.97 yr-1. If this modulation is attributed to the RSFP (Resonant Spin Flavor Precession) process, it provides a measurement of the rotation rate deep in the radiative zone. These two results suggest that the core rotates significantly more slowly than the radiative zone. If one accepts an upper limit of 7 MG for the Sun's internal magnetic field, an RSFP interpretation of the Super-Kamiokande results leads to a lower limit of 10-12 Bohr magnetons for the neutrino transition magnetic moment.

The usual procedure for estimating the significance of a peak in a power spectrum is to calculate the probability of obtaining that value or a larger value by chance, on the assumption that the time series contains only noise (e.g. that the measurements were derived from random samplings of a Gaussian distribution). However, it is known that one should regard this P-Value approach with caution. As an alternative, we here examine a Bayesian approach to estimating the significance of a peak in a power spectrum. This approach requires that we consider explicitly the hypothesis that the time series contains a periodic signal as well as noise. The challenge is to identify a probability distribution function for the power that is appropriate for this hypothesis. We propose what seem to be reasonable conditions to require of this function, and then propose a simple function that meets these requirements. We also propose a consistency condition, and check to see that our function satisfies this condition. We find that the Bayesian significance estimates are considerably more conservative than the conventional estimates. We apply this procedure to three recent analyses of solar neutrino data: (a) bimodality of GALLEX data; (b) power spectrum analysis of Super-Kamiokande data; and (c) the combined analysis of radiochemical neutrino data and irradiance data.

A search for any particular feature in any single solar neutrino dataset is unlikely to establish variability of the solar neutrino flux since the count rates are very low. It helps to combine datasets, and in this article we examine data from both the Homestake and GALLEX experiments. These show evidence of modulation with a frequency of 11.85 yr-1, which could be indicative of rotational modulation originating in the solar core. We find that precisely the same frequency is prominent in power spectrum analyses of the ACRIM irradiance data for both the Homestake and GALLEX time intervals. These results suggest that the solar core is inhomogeneous and rotates with sidereal frequency 12.85 yr-1. We find, by Monte Carlo calculations, that the probability that the neutrino data would by chance match the irradiance data in this way is only 2 parts in 10,000. This rotation rate is significantly lower than that of the inner radiative zone (13.97 yr-1) as recently inferred from analysis of Super-Kamiokande data, suggesting that there may be a second, inner tachocline separating the core from the radiative zone. This opens up the possibility that there may be an inner dynamo that could produce a strong internal magnetic field and a second solar cycle.

There has for some time been evidence of variability in radiochemical solar neutrino measurements, but this evidence has seemed suspect since the cerenkov experiments have not shown similar evidence of variability. The present re-analysis of Super-Kamiokande data shows strong evidence of r-mode oscillations. The frequencies of these oscillations correspond to a region with a sidereal rotation rate of 13.97 year-1, not far from the value 13.88 year-1 found some time ago in Homestake data, or 13.87 year-1 found recently in GALLEX data. These estimates are incompatible with the rotation rate in the convection zone, but fully compatible with current estimates of the rotation rate in the radiative zone, including the solar core. These results are suggestive of variability originating in the core from fluctuating and asymmetric nuclear burning.

Time-frequency analysis of data from the GALLEX and GNO solar neutrino experiments shows that some features in power-spectrum analyses of those datasets are due to aliasing. We normalize data for each of the four experiments in the GALLEX series, concatenate the resulting normalized data, and then form time-frequency displays. We also form time-frequency displays of the two principal modulations found in the data, at 11.87 yr-1 and at 13.63 yr-1, and of the datasets formed by subtracting these modulations from the actual (normalized) data. The results suggest that the true modulation is that at 11.87 yr-1, consistent with the results of our earlier analysis of Homestake solar neutrino data. Comparison with helioseismology data indicates that modulation is occurring either in the radiative zone just below the tachocline, presumably by the RSFP (resonant spin-flavor precession) process, or in the core, presumably due to inhomogeneities and fluctuations in the nuclear-burning process.

A histogram display of the solar neutrino capture-rate measurements made by the GALLEX experiment appears to be bimodal, but that of the follow-on GNO experiment does not. To assess the significance of these results, we introduce a "bimodality index" based on the probability-transform procedure. This confirms that the GALLEX measurements are indeed bimodal (at the 99.98 percent confidence level), and that the GNO measurements are not. Tracking the bimodality index as a function of time shows that the strongest contribution to bimodality comes from runs 42 to 62, i.e. from the time interval 1995.1 to 1996.9. The bimodality index for the first half (runs 1 through 33) is 2.56, whereas that for the second half (runs 33 through 65) is 7.05. Power-spectrum analysis shows a similar distinction: the peaks in the power spectrum formed from the second half are stronger than those in the power spectrum formed from the first half, suggesting that bimodality and rotational modulation are related.

According to current practice, the results of each run of a radiochemical solar neutrino experiment comprise an estimate of the flux and upper and lower error estimates. These estimates are derived by a maximum-likelihood procedure from the times of decay events in the analysis chamber. This procedure has the following shortcomings: (a) Published results sometimes include negative flux estimates. (b) Even if the flux estimate is non-negative, the probability distribution function implied by the flux and error estimates will extend into negative territory; and (c) The overall flux estimate derived from the results of a sequence of runs may differ substantially from an estimate made by a global analysis of all of the timing data taken together. These defects indicate that the usual packaging of data in radiochemical solar neutrino experiments provides an inadequate summary of the data, which implies a loss of information. This article reviews this problem from a Bayesian perspective, and suggests an alternative scheme for the packaging of radiochemical solar neutrino data, which is we believe free from the above objections.

Althogh KamLAND apparently rules out Resonant-Spin-Flavor-Precession (RSFP) as an explanation of the solar neutrino deficit, the solar neutrino fluxes in the Cl and Ga experiments appear to vary with solar rotation. Added to this evidence, summarized here, a power spectrum analysis of the Super-Kamiokande data reveals significant variation in the flux matching a dominant rotation rate observed in the solar magnetic field in the same time period. Three frequency peaks, all related to this rotation rate, can be explained quantitatively. A Super-Kamiokande paper reported no time variation of the flux, but showed the same peaks, there interpreted as statistically insignificant, due to an inappropriate analysis. This modulation is small (7%) in the Super-Kamiokande energy region (and below the sensitivity of the Super-Kamiokande analysis) and is consistent with RSFP as a subdominant neutrino process in the convection zone. The data display effects that correspond to solar-cycle changes in the magnetic field, typical of the convection zone. This subdominant process requires new physics: a large neutrino transition magnetic moment and a light sterile neutrino, since an effect of this amplitude occurring in the convection zone cannot be achieved with the three known neutrinos. It does, however, resolve current problems in providing fits to all experimental estimates of the mean neutrino flux, and is compatible with the extensive evidence for solar neutrino flux variability.

In solar physics, especially in exploratory stages of research, it is often necessary to compare the power spectra of two or more time series. One may, for instance, wish to estimate what the power spectrum of the combined data sets might have been, or one may wish to estimate the significance of a particular peak that shows up in two or more power spectra. One may also on occasion need to search for a complex of peaks in a single power spectrum, such as a fundamental and one or more harmonics, or a fundamental plus sidebands, etc. Visual inspection can be revealing, but it can also be misleading. This leads one to look for one or more ways of forming statistics, which readily lend themselves to significance estimation, from two or more power spectra. We derive formulas for statistics formed from the sum, the minimum, and the product of two or more power spectra. A distinguishing feature of our formulae is that, if each power spectrum has an exponential distribution, each statistic also has an exponential distribution.

It is often necessary to compare the power spectra of two or more time series: one may, for instance, wish to estimate what the power spectrum of the combined data sets might have been, or one may wish to estimate the significance of a particular peak that shows up in two or more power spectra. Also, one may occasionally need to search for a complex of peaks in a single power spectrum, such as a fundamental and one or more harmonics, or a fundamental plus sidebands, etc. Visual inspection can be revealing, but it can also be misleading. This leads one to look for one or more ways of forming statistics, which readily lend themselves to significance estimation, from two or more power spectra. The familiar chi-square statistic provides a convenient mechanism for combining variables drawn from normal distributions, and one may generalize the chi-square statistic to be any function of any number of variables with arbitrary distributions. In dealing with power spectra, we are interested mainly in exponential distributions. One well-known statistic, formed from the sum of two or more variables with exponential distributions, satisfies the gamma distribution. We show that a transformation of this statistic has the convenient property that it has an exponential distribution. We introduce two additional statistics formed from two or more variables with exponential distributions. For certain investigations, we may wish to study the minimum power (as a function of frequency) drawn from two or more power spectra. In other investigations, it may be helpful to study the product of the powers. We give numerical examples and an example drawn from our solar-neutrino research.

While KamLAND apparently rules out Resonant-Spin-Flavor-Precession (RSFP) as an explanation of the solar neutrino deficit, the solar neutrino fluxes in the Cl and Ga experiments vary with solar rotation rates. Added to this evidence, summarized here, a power spectrum analysis of the Super-Kamiokande data reveals (99.9% CL) an oscillation in the band of twice the equatorial rotation frequencies of the solar interior. An m=2 magnetic structure and RSFP, perhaps as a subdominant process, would give this effect. Solar cycle data changes are seen, as expected for convection zone modulations.

It is often necessary to compare the power spectra of two or more time series. One may, for instance, wish to estimate what the power spectrum of the combined data sets might have been. One might also wish to estimate the significance of a particular peak that shows up in two or more power spectra. Visual comparison can be revealing, but it can also be misleading. This leads one to look for one or more ways of forming statistics, which lend themselves to significance estimation, from two or more power spectra. We here propose two such statistics, one that is most useful for the combined study of two or more similar time series, and another that is more useful for the study of dissimilar time series.

We carry out a comparative analysis of the GALLEX-GNO solar neutrino data and estimates of the solar internal rotation rate derived from the MDI helioseismology experiment on the SOHO spacecraft. We introduce a statistic, which we evaluate as a function of radius and latitude, that is a measure of the degree of "resonance" of oscillations in the neutrino flux and the synodic solar rotation rate at that radius and latitude. A map of this statistic indicates that the probable location is deep in the convection zone near the equator. We also examine the integral of this statistic over the equatorial section of the convection zone. This provides a measure of the likelihood that the variability of the solar neutrino flux, as measured by the GALLEX-GNO data, has its origin in the equatorial section of the convection zone. We apply the shuffle test, randomly reassigning measurements among runs, to estimate the significance of the value of the statistic computed from the actual data. This test implies that the result is significant at the 0.2% level. When, for comparison, we repeat this analysis for the radiative zone, we find that the integral resonance statistic is not significant. These results support earlier evidence for rotational modulation of the solar neutrino flux.

If the solar neutrino flux were constant, as is widely assumed, the histogram of flux measurements would be unimodal. On the other hand, sinusoidal or square-wave modulation (either periodic or stochastic) may lead to a bimodal histogram. We here present evidence that the neutrino flux histogram is in fact bimodal. We analyze all available data from gallium experiments, coordinating results from the GALLEX and GNO experiments into one data set, and adopting results from the SAGE experiment as another data set. The two histograms, from the two data sets, are consistent in showing peaks in the range 45-75 SNU and 90-120 SNU, with a valley in between. By combining the data into one data set, we may form more detailed histograms; these strengthen the case that the flux is bimodal. A preliminary statistical analysis indicates that the bimodal character of the solar neutrino flux is highly significant. Since the upper peak is close to the expected flux (120-140 SNU), we may infer that the neutrino deficit is due to time-varying attenuation of the flux produced in the core. We estimate the time scale of this variation to be in the range 10-60 days. Attenuation that varies on such a time scale is suggestive of the influence of solar rotation, and points towards a process involving the solar magnetic field in conjunction with a nonzero neutrino magnetic moment.

Recent analysis of the Homestake data indicates that the solar neutrino flux contains a periodic variation that may be attributed to rotational modulation occurring deep in the solar interior, either in the tachocline or in the radiative zone. This paper presents an analysis of GALLEX data that yields supporting evidence of this rotational modulation at the 0.1% significance level. The depth of modulation inferred from the rotational signature is large enough to explain the neutrino deficit. The Rieger 157-day periodicity, first discovered in solar gamma-ray flares, is present also in Homestake data. A related oscillation with period 52 days is found in the GALLEX data. The relationship of these periods to the rotational period inferred from neutrino data suggests that they are due to r-mode oscillations.

In order to explore a recent proposal that the solar core may contain a component that varies periodically with a period in the range 21.0 - 22.4 days, due either to rotation or to some form of oscillation, we have examined the time series formed from measurements of the solar neutrino flux by means of the GALLEX, Homestake and Kamiokande experiments. Direct Fourier transform analysis of the Homestake data shows that the most prominent peak in the entire spectrum (examined down to 5 days period) is found at a frequency of approximately 17.2 y-1 corresponding to a period of approximately 21.3 days. According to the "shuffle test," the probability of finding this large a peak in the prescribed search band is about 0.03%, if it is assumed that there is no correlation between count rate and time. The GALLEX and Kamiokande data are examined in a way that searches for similarity in the shapes of the two spectra in sliding windows in frequency. We find that the "spectral correlation measure" peaks at 17.2 y-1, and the shuffle test indicates that the probability of finding this large a peak at a specified frequency is 2%, if it is again assumed that for each time series there is no correlation between count rate and time. The combined significance estimate is of order 1 part in 105 that the results are due to chance, on the assumption that there is no real structure to the count-rate time series.

It has been established that the Rieger periodicity of approximately 153 days is part of a complex of periodicities, all multiples of a basic period of approximately 25.5 days. However, it has not been clear why the sixth subharmonic of this periodicity should be preferentially manifested. We here note that if the Sun contains two rotating elements, with different periods and different axes, a special situation will arise if the two periods have a lowest common multiple, for in this case the relative configuration of the two rotators would repeat exactly at that (LCM) period. This chain of thought leads us to suspect that the Sun contains a second rotating element with rotational period in the range 21 - 22.5 days.