Spectral Trends of Solar Bursts at Sub-THz Frequencies

Previous sub-THz studies were derived from single-event observations. We here analyze for the first time spectral trends for a larger collection of sub-THz bursts. The collection consists of a set of 16 moderate to small impulsive solar radio bursts observed at 0.2 and 0.4 THz by the Solar Submillimeter-wave Telescope (SST) in 2012?–?2014 at El Leoncito, in the Argentinean Andes. The peak burst spectra included data from new solar patrol radio telescopes (45 and 90 GHz), and were completed with microwave data obtained by the Radio Solar Telescope Network, when available. We critically evaluate errors and uncertainties in sub-THz flux estimates caused by calibration techniques and the corrections for atmospheric transmission, and introduce a new method to obtain a uniform flux scale criterion for all events. The sub-THz bursts were searched during reported GOES soft X-ray events of class C or larger, for periods common to SST observations. Seven out of 16 events exhibit spectral maxima in the range 5?–?40 GHz with fluxes decaying at sub-THz frequencies (three of them associated to GOES class X, and four to class M). Nine out of 16 events exhibited the sub-THz spectral component. In five of these events, the sub-THz emission fluxes increased with a separate frequency from that of the microwave spectral component (two classified as X and three as M), and four events have only been detected at sub-THz frequencies (three classified as M and one as C). The results suggest that the THz component might be present throughout, with the minimum turnover frequency increasing as a function of the energy of the emitting electrons. The peculiar nature of many sub-THz burst events requires further investigations of bursts that are examined from SST observations alone to better understand these phenomena.     

A 3-dimensional calculation of the atmospheric neutrino fluxes

An extensive 3-dimensional Monte Carlo calculation of the atmospheric neutrino flux is in progress with the FLUKA Monte Carlo code. The results are compared to those obtained under the 1-dimensional approximation, where secondary particles and decay products are assumed to be collinear to the primary cosmic ray, as usually done in most of the already existing flux calculations. It is shown that the collinear approximation gives rise to a wrong angular distribution of neutrinos, essentially in the Sub-GeV region. However, the angular smearing introduced by the experimental inability of detecting recoils in neutrino interactions with nuclei is large enough to wash out, in practice, most of the differences between 3-dimensional and 1-dimensional flux calculations. Therefore, the use of the collinear approximation should not have introduced a significant bias in the determination of the flavor oscillation parameters in current experiments.     

Neutrino Event Rates from Gamma-Ray Bursts

We recalculate the diffuse flux of high-energy neutrinos produced by gamma-ray bursts in the relativistic fireball model. Although we confirm that the average single burst produces only approximately 10-2 high-energy neutrino events in a detector with a 1 km2 effective area, i.e., about 10 events yr-1, we show that the observed rate is dominated by burst-to-burst fluctuations that are very large. We find event rates that are expected to be larger by 1 order of magnitude, likely more, which are dominated by a few very bright bursts. This greatly simplifies their detection.     

Field-aligned suprathermal electron fluxes below 270 km in the auroral zone

Observations are reported of field aligned etectron fluxes in the energy range 50–500 eV at altitudes below 270 km from two rocket flights in the auroral zone. The regions of field aligned suprathermal electrons occurred in bursts of a few seconds duration, and in some instances the energy of the peak field aligned flux was in the range 100–500 eV. Theoretical calculations of the pitch angle distribution were made using the Monte Carlo technique for two model atmospheres having exospheric temperatures of 750 and 1500 K bracketing the expected auroral zone exospheric temperature. The calculations were made for the case of incident field aligned suprathermal fluxes with no local parallel electric field and also for the case of a local constant parallel electric field. Comparison of theoretical and experimental pitch angle distributions showed that in one case at 270 km a parallel electric field of 1–2 mV/m fitted the data whereas another burst at 210 km required a parallel electric field of about 10 mV/m to produce a field aligned distribution of 230 eV electrons as pronounced as was observed. Furthermore in this latter case the lack of strong field alignment at 500 eV pointed to localisation of the parallel electric field to an altitude range of 20–30 km about the rocket altitude.     

Research on neutrino events of SN 1987a

We carried out a Monte Carlo simulation calculation on the observed data of neutrino events of the supernova SN1987A. Our results show that the best value for the rest-mass of neutrino is 3–4 eV. We also calculated the light curve of the neutrino bursts and the energy spectrum.Our results impose constraints on the rest-mass of neutrino and the mass of the precursor of SN 1987A.     

Comparison of the FLUKA calculations with CAPRICE94 data on muons in atmosphere

In order to benchmark the three-dimensional calculation of the atmospheric neutrino flux based on the FLUKA Monte Carlo code, muon fluxes in the atmosphere have been computed and compared with data taken by the CAPRICE94 experiment at ground level and at different altitudes in the atmosphere. For this purpose only two additions have been introduced with respect to the neutrino flux calculation: the specific solar modulation corresponding to the period of data taking and the bending of charged particles in the atmosphere. Results are in good agreement with experimental data, although improvements in the model are possible. At this level, however, it is not possible to disentangle the interplay between the primary flux and the interaction model.     

Shock-associated kilometric radio emission and solar metric type II bursts

We present statistics relating shock-associated (SA) kilometric bursts (Cane et al., 1981) to solar metric type II bursts. An SA burst is defined here to be any 1980 kHz emission temporally associated with a reported metric type II burst and not temporally associated with a reported metric type III burst. In this way we extend to lower flux densities and shorter durations the original SA concept of Cane et al. About one quarter of 316 metric type II bursts were not accompanied by any 1980 kHz emission, another quarter were accompanied by emission attributable to preceding or simultaneous type III bursts, and nearly half were associated with SA bursts. We have compared the time profiles of 32 SA bursts with Culgoora Observatory dynamic spectral records of metric type II bursts and find that the SA emission is associated with the most intense and structured part of the metric type II burst. On the other hand, the generally poor correlation found between SA burst profiles and Sagamore Hill Observatory 606 and 2695 MHz flux density profiles suggests that most SA emission is not due to energetic electrons escaping from the microwave emission region. These results support the interpretation that SA bursts are the long wavelength extension of type II burst herringbone emission, which is presumed due to the shock acceleration of electrons.Also: Department of Physics and Astronomy, University of Maryland, College Park, MD 20742, U.S.A.     

Discovery of type I X-ray bursts from the low-mass X-ray binary 4U 1708 – 40

We report the discovery of type I X-ray bursts from the low-mass X-ray binary  4U 1708 − 40  during the 100-ks observation performed by BeppoSAX on 1999 August 15–16. Six X-ray bursts have been observed. The unabsorbed 2–10 keV fluxes of the bursts range from ∼3 to  9 × 10⁻¹⁰ erg cm⁻² s⁻¹  . A correlation between peak flux and fluence of the bursts is found, in agreement with the behaviour observed in other similar sources. There is a trend of the burst flux to increase with the time interval from the previous burst. From the value of the persistent flux we infer a mass accretion rate     , which may correspond to the mixed hydrogen/helium burning regime triggered by thermally unstable hydrogen. We have also analysed a BeppoSAX observation performed on 2001 August 22 and previous RXTE observations of  4U 1708 − 40  , where no bursts have been observed; we find persistent fluxes of more than a factor of 7 higher than the persistent flux observed during the BeppoSAX observation showing X-ray bursts.     

Evidence the Galactic contribution to the IceCube astrophysical neutrino flux

We show that the Galactic latitude distribution of IceCube astrophysical neutrino events with energies above 100 TeV is inconsistent with the isotropic model of the astrophysical neutrino flux. Namely, the Galactic latitude distribution of the events shows an excess at low latitudes |b| < 10° and a deficit at high Galactic latitude |b| ≳50°. We use Monte–Carlo simulations to show that the inconsistency of the isotropic signal model with the data is at ≳3σ level, after the account of trial factors related to the choice of the low-energy threshold and Galactic latitude binning in our analysis.     

Development of flare morphology in X-rays,and the flare scenario

We define the impulsive phase of a flare as its first phase, characterized by: X-ray bursts of short (seconds to tens of seconds) duration, a patchy X-ray morphology, and injection of energy. It lasts some five to ten minutes. The gradual or diffuse phase starts virtually at the same time as the impulsive one and is characterized by a gradually varying X-ray flux from a larger, diffuse, area situated higher than the sources of the impulsive X-ray bursts. The diffuse cloud is initially (during the first five minutes) hotter by a few million degrees than the sources of the impulsive phase bursts and is assumed to be caused by convective motions with upward velocities of a few hundred km s^?1. It contains about the same number of energetic electrons as the impulsive burst patches contained initially. It cools gradually down by radiative and conductive losses, a process that may last for about an hour.     

Solar neutrino flux,cosmic rays,and the solar activity cycle

It is suggested that the experimental data on the solar neutrino flux as measured by Davis and his collaborators from 1970 to 1982 vary with the solar activity cycle to a very high level of statistical significance for all the available tests of the hypothesis (e.g., (t-test, ²-test, run test, Wilcoxon-Mann-Whitney test) when the solar neutrino flux data are computed from the weighted moving averages of order 5. The above tests have also been applied to the data that have been generated by the Monte Carlo simulation with production rate and background rate parameters that are typical of those in the actual experiment. It is shown that the Monte Carlo simulated data do not indicate a variation within the solar cycle. Thus the moving average data strongly favours the variation within the solar activity cycle.     

Spectral distributions of microwave bursts

The lack of open literature publication of the distributional properties of the cm-λ spectra of solar microwave bursts has lead to some erroneous concepts of the typical characteristics of these spectra. To provide more accurate information, this paper sets forth various distributions of the peak flux density spectra of large numbers of bursts, based on observations of the Sagamore Hill Radio Observatory at nine discrete frequencies between 245 and 35000 MHz over the years 1968–1971. As a foundation for the distribution studies, the basic spectral classification system is outlined. The majority of burst spectra were found to contain a cm-λ component having a single spectral maximum in the 1400 to 35000 MHz range; such spectra are designated C type. A study of the correlation of the spectral maximum frequency f _max of the cm component and the photospheric magnetic field strength of the associated region shows a tendency for greater correlation at higher f _max for stronger magnetic sssfields. A study of the correlation for C type spectra between f _max and the quasi-cutoff frequency f _qc on the low-frequency side shows that for bursts of moderate peak flux density (50–500 sfu) f _qc is well correlated with f _max; a good fit to the relation f _max=A f _qc is found with A =3.4. The possible attenuating mechanisms responsible for the spectral shaping of the cm component are discussed.     

Characteristics of soft solar X-ray bursts

The burst component of the solar X-ray flux in the soft wavelength range 2 < < 12 Å observed from Explorer 33 and Explorer 35 from July 1966 to September 1968 was analyzed. In this period 4028 burst peaks were identified.The differential distributions of the temporal and intensity parameters of the bursts revealed no separation into more than one class of bursts. The most frequently observed value for rise time was 4 min and for decay time was 12 min. The distribution of the ratio of rise to decay time can be represented by an exponential with exponent -2.31 from a ratio of 0.3 to 2.7; the maximum in this distribution occurred at a ratio of 0.3. The values of the total observed flux, divided by the background flux at burst maximum, can be represented by a power law with exponent -2.62 for ratios between 1.5 and 32. The distribution of peak burst fluxes can be represented by a power law with exponent - 1.75 over the range 1–100 milli-erg (cm² sec)^–1. The flux time integral values are given by a power law with exponent -1.44 over the range 1–50 erg cm^–2.The distribution of peak burst flux as a function of H importance revealed a general tendency for larger peak X-ray fluxes to occur with both larger H flare areas and with brighter H flares. There is no significant dependence of X-ray burst occurrence on heliographic longitude; the emission thus lacks directivity.The theory of free-free emission by a thermal electron distribution was applied to a composite quantitative discussion of hard X-ray fluxes (data from Arnoldy et al., 1968; Kane and Winckler, 1969; and Hudson et al., 1969) and soft X-ray fluxes during solar X-ray bursts. Using bursts yielding measured X-ray intensities in three different energy intervals, covering a total range of 1–50 keV, temperatures and emission measures were derived. The emission measure was found to vary from event to event. The peak time of hard X-ray events was found to occur an average of 3 min before the peak time of the corresponding soft X-ray bursts. Thus a changing emission measure during the event is also required. A free-free emission process with temperatures of 12–39 × 10⁶K and with an emission measure in the range 3.6 × 10⁴⁷ to 2.1 × 10⁵⁰ cm^–3 which varies both from event to event and within an individual event is required by the data examined.Now at Department of Astrophysical Sciences, Princeton University, Princeton, New Jersey.     

GRBs on probation: Testing the UHE CR paradigm with IceCube

Gamma ray burst (GRB) fireballs provide one of very few astrophysical environments where one can contemplate the acceleration of cosmic rays to energies that exceed 10²⁰ eV. The assumption that GRBs are the sources of the observed cosmic rays generates a calculable flux of neutrinos produced when the protons interact with fireball photons. With data taken during construction IceCube has already reached a sensitivity to observe neutrinos produced in temporal coincidence with individual GRBs provided that they are the sources of the observed extra-galactic cosmic rays. We here point out that the GRB origin of cosmic rays is also challenged by the IceCube upper limit on a possible diffuse flux of cosmic neutrinos which should not be exceeded by the flux produced by all GRB over Hubble time. Our alternative approach has the advantage of directly relating the diffuse flux produced by all GRBs to measurements of the cosmic ray flux. It also generates both the neutrino flux produced by the sources and the associated cosmogenic neutrino flux in a synergetic way.     

Role of non-thermal electrons on the generation of X-ray and EUV lines in solar flares

The energy and angular distributions of electrons have been studied by combining small angle scatterings using analytical treatment with large angle collisions using Monte Caroo calculations as a function of column density for initially power-law electron distributions and incidence angles of 0, 30, and 60°. Using these distributions the X-ray and EUV line flux as a function of column density has been computed. The flux increases with increase in column density. At the initial column densities the contribution of non-thermal electrons for the production of line flux is negligible. However, it becomes significant at intermediate column densities at which the electron energy and angular distributions have non-Maxwellian nature. X-ray and EUV flux have also been calculated as a function of electron spectral index at a fixed column density. It falls steeply with increase in spectral index. The calculated flux is compared with the observations.     

The influence of magnetic fields on neutrino-dominated accretion disc

We consider the influence of magnetic fields on the model of neutrino-dominated accretion flows (NDAFs) for gamma-ray bursts (GRBs) via the assumption that the accretion rate of the disc is totally caused by the torque of the Lorentz force, i.e. the magnetic braking of large-scale magnetic fields and magnetic viscosity of small-scale magnetic fields. We calculate the structure, composition, luminosity of neutrino emission and the Poynting flux, and the rate of mass loss driven by neutrino heating or launched centrifugally by large-scale magnetic fields, based on the physical condition of the magnetized NDAFs. It is shown that the magnetized disc is favourable to interpret the diverse prompt emissions as well as the X-ray flares observed in the early afterglow of GRBs.     

An index of impulsiveness for 2800 MHz impulsive solar noise bursts

The initial growth of intensity, I, of the impulsive solar noise burst observed at 2800 MHz has been fitted with a parabolic curve of the form I=ct ²and the quantity c taken as the index of impulsiveness. Two groups of bursts comprising 85% of all impulsive bursts observed in 1962–63 and 1966 were selected for study. A good fit has been obtained for bursts having peak flux density up to 20 flux units, while for more intense bursts, the average observed growth is more rapid than the parabolic rate. The distribution of the index in the range 0.1 to more than 100 shows two peaks, one for c values 1–10 and another apparent peak for those with c greater than 100. The index is independent of the peak intensity of the burst and its position on the solar disc, while there is a small trend, indicating that shorter bursts are more impulsive than longer duration events. The more easily derived linear rate of rise, b = Peak Flux/Interval from start to peak is related to the parabolic impulsive index by b = 1.86 c ^0.57. The non-linear rate of expansion of a flaring volume suggested by Pneuman when applied to explain the parabolic rise of microwave bursts indicates that the impulsiveness of bursts is inversely related to the contained magnetic field.     

Searching for neutrinos from WIMP annihilations in the Galactic stellar disc

Weakly interacting massive particles (WIMPs) are a viable candidate for the relic abundance of dark matter (DM) produced in the early universe. So far, WIMPs have eluded direct detection through interactions with baryonic matter. Neutrino emission from accumulated WIMP annihilations in the solar core has been proposed as a signature of DM, but has not yet been detected. These null results may be due to small-scale DM density fluctuations in the halo with the density of our local region being lower than the average  (∼0.3 GeV cm⁻³)  . However, the accumulated neutrino signal from WIMP annihilations in the Galactic stellar disc would be insensitive to local density variations. Inside the disc, DM can be captured by stars causing an enhanced annihilation rate and therefore a potentially higher neutrino flux than what would be observed from elsewhere in the halo. We estimate a neutrino flux from the WIMP annihilations in the stellar disc to be enhanced by more than an order of magnitude compared to the neutrino fluxes from the halo. We offer a conservative estimate for this enhanced flux, based on the WIMP–nucleon cross-sections obtained from direct-detection experiments by assuming a density of  ∼0.3 GeV cm⁻³  for the local DM. We also compare the detectability of these fluxes with a signal of diffuse high-energy neutrinos produced in the Milky Way by the interaction of cosmic rays with the interstellar medium. These comparative signals should be observable by large neutrino detectors.