Comparative Analyses of Brookhaven National Laboratory Nuclear Decay Measurements and Super-Kamiokande Solar Neutrino Measurements: Neutrinos and Neutrino-Induced Beta-Decays as Probes of the Deep Solar Interior

An experiment carried out at the Brookhaven National Laboratory over a period of almost 8 years acquired 364 measurements of the beta-decay rates of a sample of \({}^{32}\mbox{Si}\) and, for comparison, of a sample of \({}^{36}\mbox{Cl}\). The experimenters reported finding “small periodic annual deviations of the data points from an exponential decay?…?of uncertain origin”. We find that power-spectrum and spectrogram analyses of these datasets show evidence not only of the annual oscillations, but also of transient oscillations with frequencies near 11 year^?1 and 12.5 year^?1. Similar analyses of 358 measurements of the solar neutrino flux acquired by the Super-Kamiokande neutrino observatory over a period of about 5 years yield evidence of an oscillation near 12.5 year^?1 and another near 9.5 year^?1. An oscillation near 12.5 year^?1 is compatible with the influence of rotation of the radiative zone. We suggest that an oscillation near 9.5 year^?1 may be indicative of rotation of the solar core, and that an oscillation near 11 year^?1 may have its origin in a tachocline between the core and the radiative zone. Modulation of the solar neutrino flux may be attributed to an influence of the Sun’s internal magnetic field by the Resonant Spin Flavor Precession (RSFP) mechanism, suggesting that neutrinos and neutrino-induced beta decays can provide information about the deep solar interior.     

Solar internal rotation rate and the latitudinal variation of the tachocline

A new set of accurately measured frequencies of solar oscillations is used to infer the rotation rate inside the Sun, as a function of radial distance as well as latitude. We have adopted a regularized least-squares technique with iterative refinement for both 1.5D inversion, using the splitting coefficients, and 2D inversion using individual m splittings. The inferred rotation rate agrees well with earlier estimates showing a shear layer just below the surface and another one around the base of the convection zone. The tachocline or the transition layer where the rotation rate changes from differential rotation in the convection zone to an almost latitudinally independent rotation rate in the radiative interior is studied in detail. No compelling evidence for any latitudinal variation in the position and width of the tachocline is found, although it appears that the tachocline probably shifts to a slightly larger radial distance at higher latitudes and possibly also becomes thicker. However, these variations are within the estimated errors and more accurate data would be needed to make a definitive statement about latitudinal variations.     

An analysis of apparent r-mode oscillations in solar activity,the solar diameter,the solar neutrino flux,and nuclear decay rates,with implications concerning the Sun’s internal structure and rotation,and neutrino processes

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–0.74, where the sidereal rotation rate is in the range 13.7–14.6 year⁻¹. 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 12.0 year⁻¹, 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 Resonant Spin Flavor Precession (RSFP) process, which leads to estimates of the neutrino magnetic moment and of the magnetic field strength in or near the solar core.     

Further Evidence Suggestive of a Solar Influence on Nuclear Decay Rates

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

Combined Analysis of Solar Neutrino and Solar Irradiance Data: Further Evidence for Variability of the Solar Neutrino Flux and Its Implications Concerning the Solar Core

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 year⁻¹, 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 a sidereal frequency of 12.85 year⁻¹. From Monte Carlo calculations, it is found 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 year⁻¹) 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.     

Intermediate degree p-mode frequency splittings near solar maximum

We report on observations of global solar Ca K-line intensity oscillations taken in May 1991 from Mees Solar Observatory, Hawaii. We measurep-mode frequency splittings for modes of spherical harmonic degrees between 20 and 129 averaged over the radial order of the modes. Our measurement of the antisymmetric component of the splittings is comparable with previous measurements and thus indicates a decrease in the latitudinal differential rotation with depth into the convection zone and the upper radiative zone. We find evidence for a 1% variation in the rotation rate of the upper convection zone roughly in phase with the solar activity cycle. Our measurement of the symmetric component of the splittings is of the same order as was reported from the previous solar maximum and is an order of magnitude larger than has been measured near solar minimum. From the degree dependence of the symmetric component of the splittings, we find evidence for an aspherical fractional sound speed perturbation located at a depth of 0.85 ± 0.05 solar radii. This perturbation has a magnitude ofc/c +9 × 10^–4 at the equator relative to the poles. Additionally, there is evidence for a near-surface aspherical sound speed perturbation of smaller magnitudec/c +4 × 10^–4 at the equator relative to the poles. If an intense global magnetic field were the dominant source of the observed symmetric component of the splittings, instead of latitudinal gradients in the sound speed, then global fields of order 10⁵ G would be required in the convection zone.     

Magnetic field and differential rotation in the radiative zone of the Sun

The problem of the interaction between magnetic fields and differential rotation in the radiative zone of the Sun is investigated. It is demonstrated that effects of magnetic buoyancy can be neglected in the analysis of this interaction. It is shown that hydromagnetic torsional waves propagating from the solar core cannot be responsible for the 22-year solar cycle. A possible geometry of the magnetic field that conforms with stationary differential rotation is considered. A verifying method for hypotheses on the structure of the magnetic field and torsional oscillations in the radiative zone of the Sun is proposed based on helioseismic data.     

Hydrodynamic Source of the 11-year Solar Variations

In the traditional axisymmetric models of the 11-year solar cycle, oscillations of the magnetic fields appear in the background of nonoscillating (over time scale considered) turbulent velocity fields and differential rotation. In this paper, an alternative approach is developed: The excitation of magnetic oscillations with the 22-year period is the consequence of hydrodynamic oscillations with the 11-year period. In the excitation of hydrodynamic oscillations, two processes taking place in high latitudes near the interface between the convective and radiative zones play a key role. One is forcing of the westerly zonal flow, the conditions for which are due to deformation of the interfacial surface. The other process is the excitation of a shear instability of zonal flow as a consequence of a strong radial gradient of angular velocity. The development of a shear instability at some stage brings about the disruption of the forcing of differential rotation. In the first (hydrodynamic) part of the paper, the dynamics of axisymmetric flows near the bottom of the convection zone is numerically simulated. Forcing of differential rotation having velocity shear in latitude and the existence of solutions in the form of torsional waves with the 11-year oscillation period are shown. In the second part the dynamics of the magnetic field is studied. The most pronounced peculiarities of the solutions are the existence of forced oscillations with the 22-year period and the drift of the toroidal magnetic field component from the mid latitudes to the equator. In high and low latitudes after cycle maximum, the toroidal component is of opposite sign in accordance with observations. In the third part, the transport of momentum from the bottom of the convection zone to the outer surface by virtue of diffusivity is considered. The existence of some sources of differential rotation in the convection zone is not implied. A qualitative correspondence of the differential rotation profile in the bulk of the convection zone and on its outer surface to experimental data is shown. The time correspondence between torsional and magnetic oscillations is also in accordance with observations.     

The slender solar tachocline: a magnetic model

A model for the sharp transition from differential rotation in the solar convection zone to rigid rotation in the radiative interior is presented. Differential rotation in the radiative zone is shown to be quenched efficiently by an internal magnetic field. The poloidal field amplitude, B₀, is the input parameter for our model which determines the transition layer thickness and the toroidal field strength. It is illustrated analytically and confirmed numerically that a rather small field, B₀ = 10⁻⁴ Gauss, suffices to satisfy the helioseismological restrictions on the depth of the differential rotation penetration below the convection zone. The transition layer thickness decreases further with increasing B₀. The toroidal field amplitude B ≃ 200 Gauss is almost independent of B₀.     

Helioseismology and the solar interior dynamics

The inversion of helioseismic modes leads to the sound velocity inside the Sun with a precision of about 0.1 per cent. Comparisoons of solar models with the “seismic sun” represent powerful tools to test the physics: depth of the convection zone, equation of state, opacities, element diffusion processes and mixing inside the radiative zone. We now have evidence that microscopic diffusion (element segregation) does occur below the convection zone, leading to a mild helium depletion in the solar outer layers. Meanwhile this process must be slowed down by some macroscopic effect, presumably rotation-induced mixing. The same mixing is also responsible for the observed lithium depletion. On the other hand, the observations of beryllium and helium 3 impose specific constraints on the depth of this mildly mixed zone. Helioseismology also gives information on the internal solar rotation: while differential rotation exists in the convection zone, solid rotation prevails in the radiative zone, and the transition layer (the so-called “tachocline”) is very small. These effects are discussed, together with the astrophysical constraints on the solar neutrino fluxes.     

The 22-year cycle of the sun

The 22-yr solar cycle is being explained as due to a spin-aligned magnetic quadrupole frozen into the (radiative) core. Differential rotation of the poorly conducting convection zone gives rise to flux winding and to an alternating dominance, near the surface, of one of the two constituent dipoles. Equatorial superrotation and polar subrotation are stabilized by the same magnetic flux that transfers angular momentum from the (radiative) core to the (escaping) solar wind. The time- and radius-dependent magnetic torque is also responsible for the rigid and twisting oscillations of a thin surface layer, and for the relative velocities of different tracers.     

The Dynamics of the Solar Radiative Zone

The dynamics of the solar radiative interior are still poorly constrained by comparison to the convective zone. This disparity is even more marked when we attempt to derive meaningful temporal variations. Many data sets contain a small number of modes that are sensitive to the inner layers of the Sun, but we found that the estimates of their uncertainties are often inaccurate. As a result, these data sets allow us to obtain, at best, a low-resolution estimate of the solar-core rotation rate down to approximately 0.2R _⊙. We present inferences based on mode determination resulting from an alternate peak-fitting methodology aimed at increasing the amount of observed modes that are sensitive to the radiative zone, while special care was taken in the determination of their uncertainties. This methodology has been applied to MDI and GONG data, for the whole Solar Cycle 23, and to the newly available HMI data. The numerical inversions of all these data sets result in the best inferences to date of the rotation in the radiative region. These results and the method used to obtain them are discussed. The resulting profiles are shown and analyzed, and the significance of the detected changes is discussed.     

Solar rotation and the giant cells

Departures from the mean solar differential rotation rate as a function of latitude, longitude, and epoch of the solar cycle, together with variations in the rotation rate as determined by spectroscopic and tracer measurements are reviewed. It is shown that, if giant convection cells do exist as predicted, real variations in the subsurface rotation rate should occur and that this may be responsible for the observed surface anomalies.In terms of this hypothesis, a simple account is given for the anomalous rotation rates of sunspots. Furthermore, the torsional oscillations are identified as a modulation of the differential rotation produced by a system of toroidal convective rolls generated near the poles and propagating towards the equator. It is suggested that, as these rolls progress through lower latitudes, they break up into a system of cells which are the long sought for giant cells of the convection zone. Thus the torsional oscillations are identified as direct surface evidence for the existence of these cells.Solar Cycle Workshop Paper.     

Differential rotation and meridional flow in the solar convection zone and beneath

The influence of the basic rotation on anisotropic and inhomogeneous turbulence is discussed in the context of differential rotation theory. An improved representation for the original turbulence leads to a Λ‐effect which complies with the results of 3D numerical simulations. The resulting rotation law and meridional flow agree well with both the surface observations (∂Ω/∂r < 0 and meridional flow towards the poles) and with the findings of helioseismology. The computed equatorward flow at the bottom of convection zone has an amplitude of about 10 m/s and may be significant for the solar dynamo. The depth of the meridional flow penetration into the radiative zone is proportional to ν^0.5_core, where ν_core is the viscosity beneath the convection zone. The penetration is very small if the tachocline is laminar. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Differential rotation and meridional flow of Arcturus

The spectroscopic variability of Arcturus hints at cyclic activity cycle and differential rotation. This could provide a test of current theoretical models of solar and stellar dynamos. To examine the applicability of current models of the flux transport dynamo to Arcturus, we compute a mean‐field model for its internal rotation, meridional flow, and convective heat transport in the convective envelope. We then compare the conditions for dynamo action with those on the Sun. We find solar‐type surface rotation with about 1/10th of the shear found on the solar surface. The rotation rate increases monotonically with depth at all latitudes throughout the whole convection zone. In the lower part of the convection zone the horizontal shear vanishes and there is a strong radial gradient. The surface meridional flow has maximum speed of 170 m/s and is directed towards the equator at high and towards the poles at low latitudes. Turbulent magnetic diffusivity is of the order 10¹⁵–10¹⁶ cm²/s. The conditions on Arcturus are not favorable for a circulation‐dominated dynamo (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

X-ray variability of SS 433: Evidence for supercritical accretion

We study the X-ray variability of SS 433 based on data from the ASCA observatory and the MAXI and RXTE/ASM monitoring missions. Based on the ASCA data, we have constructed the power spectrum of SS 433 in the frequency range from 10^?6 to 0.1 Hz, which confirms the presence of a flat portion in the spectrum at frequencies 3 × 10^?5?10^?3 Hz. The periodic variability (precession, nutation, eclipses) begins to dominate significantly over the stochastic variability at lower frequencies, which does not allow the stochastic variability to be studied reliably. The model in which the flat portion extends to 9.5 × 10^?6 Hz, while a power-law rise with an index of 2.6 occurs below provides the best agreement with the observations. The nutational oscillations of the jets with a period of about three days suggests that the time for the passage of material through the disk is less than this value. At frequencies below 4 × 10^?6 Hz, the shape of the power spectrum probably does not reflect the disk structure but is determined by external factors, for example, by a change in the amount of material supplied by the donor. The flat portion can arise from a rapid decrease in the viscous time in the supercritical or radiative disk zones. The flat spectrum is associated with the variability of the X-ray jets that are formed in the supercritical disk region.     

Present problems of the solar interior

The standard model of solar evolution is reviewed and a number of problems highlighted. A fundamental question is whether there is any mixing of matter in the central regions, since such mixing could radically alter the model of the present Sun and modify our understanding of the evolution of other stars. Standard models of solar evolution become unstable to ³He driven global oscillations at an age of 3 × 10⁸ years and this may drive some mixing, even if this is not the case the finite amptitude limit of these oscillations is likely to produce modifications in the standard model. Convective overshooting at the bottom of the outer convective zone leads to an increased depth of this zone and small changes in the interior. It is pointed out that the young Sun had a ¹²C driven convective core whose extent and duration depends on the extent of overshooting. Such a core is likely to produce a magnetic field which will affect the internal dynamics. The internal rotation of the Sun remains an enigma and absence of knowledge of any internal magnetic field makes it difficult to study the problem. Rotationally driven instabilities are ineffective in the central chemically inhomogenous regions but may contribute to the inward diffusion of lithium from the convective zone. These and other problems are considered, but few solutions are proposed.     

Evidence for solar neutrino flux variability and its implications

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

Evidence for solar neutrino flux variability and its implications

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

Effect of Gravity Waves on Transport Process in Stellar Radiative Region

The aim of this work is to present a transport process which is likely to have a great importance for the internal constitution of the stars. In order to set the problem, we first give a short presentation of the physical properties of the Sun and stars, described usually under the names of `Standard Solar Model' or `Standard Stellar Models' (SSM). Next we show that an important question about SSM is that they do not explain the age dependance of lithium deficiency of stars of known age: stars of galactic clusters and the Sun. It has been suggested a long time ago to assume the presence of a macrosocpic diffusion process in the radiative zone, below the surface convective zone of solar like stars. It is then possible for the lithium present in the convective zone to be carried to the thermonuclear burning level below the convective zone. The first assumption was that differential rotation generates turbulence and therefore that a turbulent diffusion process takes place. However, this model predicts a lithium abundance which is strongly rotation dependant, contrary to the observations. Furthermore, the diffusion coefficient being large all over the radiative zone, it prevents the possibility of gravitational separation by diffusion and consequently leads to an impossibility of explaining the difference of helium abundance between the surface and the center of the Sun. The consequence is obviously that we need to take into account another physical process. Stars having a mass M < 1.3 M _⊙ have a convective zone which begins close to the stellar surface and extends down to a depth which is an appreciable fraction of stellar radius. In the convective zone, strong stochastic motions take care, at least partially, of heat transfer. These motions do not vanish at the lower boundary and generate internal waves into the radiative zone. These random internal waves are at the origin of a diffusion process which can be considered as responsible of the diffusive transport of lithium down to the lithium burning level. This is certainly not the only physical process responsible of lithium deficiency in main sequence stars, but its properties open the way to a completely consistent analysis of lithium deficiency. The model of generation of gravity waves is based on a model of heat transport in the convective zone by diving plumes. The horizontal component of the turbulent motion at the boundary of the convective zone is supposed to generate the horizontal motion of internal waves. The result is a large horizontal component of the diffusion coefficient, which produces in a short time an horizontal uniform chemical composition. It is known that gravity waves, in the absence of any dissipative process, cannot generate vertical mixing. Therefore, the vertical component of the diffusion coefficient is entirely dependant of radiative damping. It decreases quickly in the radiative zone, but is large enough to be responsible of lithium burning. Due to the radial dependance of velocity amplitude, the diffusion coeficient increases when approaching the stellar center. However, very close to the center, non-linear dissipative and radiative damping of internal waves become large and the diffusion coefficient vanishes at the very center. The development of this abstract can be found in E. Schatzman (1996, J. Fluid Mech. 322, 355). This revised version was published online in July 2006 with corrections to the Cover Date.