Comparative Analysis of Super-Kamiokande and SNO Solar-Neutrino Data and the Photospheric Magnetic Field

We analyze Super-Kamiokande, SNO, and photospheric magnetic-field data for the common time interval, namely the SNO D2O phase. Concerning rotational modulation, the magnetic-field power spectrum shows the strongest peaks at the second and sixth harmonics of the solar synodic rotation frequency [3 nu(rot) and 7 nu(rot)]. The restricted Super-Kamiokande dataset shows strong modulation at the second harmonic. The SNO D2O dataset shows weak modulation at that frequency, but strong modulation in the sixth-harmonic frequency band. We estimate the significance level of the correspondence of the Super-Kamiokande second-harmonic peak with the corresponding magnetic-field peak to be 0.0004, and the significance level of the correspondence of the SNO D2O sixth-harmonic peak with the corresponding magnetic-field peak to be 0.009. By estimating the amplitude of the modulation of the solar neutrino flux at the second harmonic from the restricted Super-Kamiokande dataset, we find that the weak power at that frequency in the SNO D2O power spectrum is not particularly surprising. Concerning 9.43 yr-1, we find no peak at this frequency in the power spectrum formed from the restricted Super-Kamiokande dataset, so it is no surprise that this peak does not show up in the SNO D2O dataset, either.Comment: 32 pages, 8 tables, 16 figure

MaxEnt power spectrum estimation using the Fourier transform for irregularly sampled data applied to a record of stellar luminosity

The principle of maximum entropy is applied to the spectral analysis of a data signal with general variance matrix and containing gaps in the record. The role of the entropic regularizer is to prevent one from overestimating structure in the spectrum when faced with imperfect data. Several arguments are presented suggesting that the arbitrary prefactor should not be introduced to the entropy term. The introduction of that factor is not required when a continuous Poisson distribution is used for the amplitude coefficients. We compare the formalism for when the variance of the data is known explicitly to that for when the variance is known only to lie in some finite range. The result of including the entropic measure factor is to suggest a spectrum consistent with the variance of the data which has less structure than that given by the forward transform. An application of the methodology to example data is demonstrated.Comment: 15 pages, 13 figures, 1 table, major revision, final version, Accepted for publication in Astrophysics & Space Scienc

Comparative analysis of Gallex and GNO solar neutrino data

Since the GALLEX and GNO datasets were derived from closely related experiments, there is a natural tendency to merge them. This is perhaps appropriate for any analysis based on the hypothesis that the solar neutrino flux is constant, but it is not necessarily appropriate for an analysis that allows for possible variability, since the GALLEX and GNO experiments belong to different solar cycles. Moreover, we find significant differences between the GALLEX and GNO datasets. It appears, from inspection of the time series and histograms, that GNO measurements are compatible with the assumption that the solar neutrino flux is constant, but GALLEX measurements are not. Furthermore, power-spectrum analysis yields evidence of rotational modulation in GALLEX data but not in GNO data. We compare our results with those of Pandola, who claims that GALLEX-GNO data show no evidence for variability.Comment: 20 pages plus 6 tables plus 11 figure

Multidimensional relativistic MHD simulations of Pulsar Wind Nebulae: dynamics and emission

Pulsar Wind Nebulae, and the Crab nebula in particular, are the best cosmic laboratories to investigate the dynamics of magnetized relativistic outflows and particle acceleration up to PeV energies. Multidimensional MHD modeling by means of numerical simulations has been very successful at reproducing, to the very finest details, the innermost structure of these synchrotron emitting nebulae, as observed in the X-rays. Therefore, the comparison between the simulated source and observations can be used as a powerful diagnostic tool to probe the physical conditions in pulsar winds, like their composition, magnetization, and degree of anisotropy. However, in spite of the wealth of observations and of the accuracy of current MHD models, the precise mechanisms for magnetic field dissipation and for the acceleration of the non-thermal emitting particles are mysteries still puzzling theorists to date. Here we review the methodologies of the computational approach to the modeling of Pulsar Wind Nebulae, discussing the most relevant results and the recent progresses achieved in this fascinating field of high-energy astrophysics.Comment: 29 pages review, preliminary version. To appear in the book "Modelling Nebulae" edited by D. Torres for Springer, based on the invited contributions to the workshop held in Sant Cugat (Barcelona), June 14-17, 201

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 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.Comment: 22 pages, 9 tables, 10 figure

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-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}.Comment: 12 pages, 6 figures, v2 has a color corrected Fig 6, a corrected reference, and a corrected typ

Time-dependent Stochastic Modeling of Solar Active Region Energy

A time-dependent model for the energy of a flaring solar active region is presented based on a stochastic jump-transition model (Wheatland and Glukhov 1998; Wheatland 2008; Wheatland 2009). The magnetic free energy of the model active region varies in time due to a prescribed (deterministic) rate of energy input and prescribed (random) flare jumps downwards in energy. The model has been shown to reproduce observed flare statistics, for specific time-independent choices for the energy input and flare transition rates. However, many solar active regions exhibit time variation in flare productivity, as exemplified by NOAA active region AR 11029 (Wheatland 2010). In this case a time-dependent model is needed. Time variation is incorporated for two cases: 1. a step change in the rates of flare jumps; and 2. a step change in the rate of energy supply to the system. Analytic arguments are presented describing the qualitative behavior of the system in the two cases. In each case the system adjusts by shifting to a new stationary state over a relaxation time which is estimated analytically. The new model retains flare-like event statistics. In each case the frequency-energy distribution is a power law for flare energies less than a time-dependent rollover set by the largest energy the system is likely to attain at a given time. For Case 1, the model exhibits a double exponential waiting-time distribution, corresponding to flaring at a constant mean rate during two intervals (before and after the step change), if the average energy of the system is large. For Case 2 the waiting-time distribution is a simple exponential, again provided the average energy of the system is large. Monte Carlo simulations of Case~1 are presented which confirm the analytic estimates. The simulation results provide a qualitative model for observed flare statistics in active region AR 11029.Comment: 25 pages, 9 figure

Observations of quasi-periodic solar X-ray emission as a result of MHD oscillations in a system of multiple flare loops

We investigate the solar flare of 20 October 2002. The flare was accompanied by quasi-periodic pulsations (QPP) of both thermal and nonthermal hard X-ray emissions (HXR) observed by RHESSI in the 3-50 keV energy range. Analysis of the HXR time profiles in different energy channels made with the Lomb periodogram indicates two statistically significant time periods of about 16 and 36 seconds. The 36-second QPP were observed only in the nonthermal HXR emission in the impulsive phase of the flare. The 16-second QPP were more pronounced in the thermal HXR emission and were observed both in the impulsive and in the decay phases of the flare. Imaging analysis of the flare region, the determined time periods of the QPP and the estimated physical parameters of magnetic loops in the flare region allow us to interpret the observations as follows. 1) In the impulsive phase energy was released and electrons were accelerated by successive acts with the average time period of about 36 seconds in different parts of two spatially separated, but interacting loop systems of the flare region. 2) The 36-second periodicity of energy release could be caused by the action of fast MHD oscillations in the loops connecting these flaring sites. 3) During the first explosive acts of energy release the MHD oscillations (most probably the sausage mode) with time period of 16 seconds were excited in one system of the flare loops. 4) These oscillations were maintained by the subsequent explosive acts of energy release in the impulsive phase and were completely damped in the decay phase of the flare.Comment: 14 pages, 4 figure

Nonlinear time-series analysis of Hyperion's lightcurves

Hyperion is a satellite of Saturn that was predicted to remain in a chaotic rotational state. This was confirmed to some extent by Voyager 2 and Cassini series of images and some ground-based photometric observations. The aim of this aticle is to explore conditions for potential observations to meet in order to estimate a maximal Lyapunov Exponent (mLE), which being positive is an indicator of chaos and allows to characterise it quantitatively. Lightcurves existing in literature as well as numerical simulations are examined using standard tools of theory of chaos. It is found that existing datasets are too short and undersampled to detect a positive mLE, although its presence is not rejected. Analysis of simulated lightcurves leads to an assertion that observations from one site should be performed over a year-long period to detect a positive mLE, if present, in a reliable way. Another approach would be to use 2---3 telescopes spread over the world to have observations distributed more uniformly. This may be achieved without disrupting other observational projects being conducted. The necessity of time-series to be stationary is highly stressed.Comment: 34 pages, 12 figures, 4 tables; v2 after referee report; matches the version accepted in Astrophysics and Space Scienc

Power Spectrum Analysis of Physikalisch-Technische Bundesanstalt Decay-Rate Data: Evidence for Solar Rotational Modulation

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.Comment: 15 pages, 13 figure