Search of ultra high energy cosmic rays’ sources. FRI-radiogalaxy Centaurus A

The propagation of ultra high energy cosmic rays in Galactic and extragalactic magnetic fields is investigated in the present paper. The motion of charged particles of different energies and chemical composition is simulated using different Galactic magnetic field models. Positions for the real sources of events registered at the Auger observatory are calculated taking into account the influence of both Galactic and extragalactic turbulent fields. The possibility of their correlation with the Centaurus A radio galaxy is analyzed.     

Is sky distribution of gamma-ray bursts random?

We investigate the full randomness of the angular distribution of gamma-ray bursts (GRBs) detected by the BATSE (the Burst and Transient Source Experiment). We divided the BATSE sample into 5 subsamples (short1, short2, intermediate, long1, long2 based on their durations and peak fluxes and studied the angular distributions separately. We used three methods, the Voronoi tesselation, minimal spanning tree and multifractal spectra to search for non-randomness in the subsamples. To investigate the eventual non-randomness in the subsamples we defined 13 test-variables (9 from the Voronoi tesselation, 3 from the minimal spanning tree and one from the multifractal spectrum). We performed Monte Carlo simulations taking into account the BATSE sky-exposure function. We tested the randomness by introducing squared Euclidean distances in the parameter space of the test-variables. We concluded that the short1, short2 groups deviate significantly (99.90%, 99.98%) from the full randomness in the distribution of the squared Euclidean distances, however, it is not the case with the long subsamples. In the intermediate group, the squared Euclidean distances also yield a significant deviation (98.51%).     

A new component of cosmic rays?

Recent direct measurements of the energy spectra of the major mass components of cosmic rays have indicated the presence of a ‘kink’ in the region of 200 GeV per nucleon. The kink, which varies in magnitude from one element to another, is much sharper than predicted by our cosmic ray origin model in which supernova remnants are responsible for cosmic ray acceleration and it appears as though a new, steeper component is responsible.The component amounts to about 20 percent of the total at 30 GeV/nucleon for protons and helium nuclei and its magnitude varies with nuclear charge; the unweighted fraction for all cosmic rays being 36%.The origin of the new component is subject to doubt but the contenders include O, B, A, supergiant and Wolf-Rayet stars, by way of their intense stellar winds. Another explanation is also in terms of these particles as the sources but then being trapped, and even further accelerated, in the Local Bubble.     

On the variability of afterglows from cosmic gamma-ray bursts—the possibility of a transition to a nonrelativistic motion

Variability on time scales δt < t is observed on numerous occasions in the afterglows of cosmic gamma-ray bursts (GRBs). It is well known that the radiation originating in an external shock produced by the interaction of an ultrarelativistic jet with the ambient interstellar medium should not contain such variability within the framework of simple models. The corresponding constraints were established by Ioka et al. (2005) and, in some instances, are inconsistent with observations. On the other hand, if the motion is not relativistic, then the rapid afterglow variability can be explained much more easily. Various estimates of the transition time to a nonrelativistic motion in a GRB source are discussed in this connection. It has been shown that this transition should occur on an observed time scale of ～10 days. In the case of a higher density of the surrounding material, ～10²?10⁴ cm^?3, or a stellar wind with ? ～ 10^?5?10^?4 M _⊙ yr^?1, the transition to a nonrelativistic motion can occur on a time scale of ～1 day. Such densities may well be expected in star-forming regions and around massive Wolf-Rayet stars.     

On the spectrum of ultrahigh energy cosmic rays and the γ-ray burst origin hypothesis

It has been suggested that cosmological γ-ray bursts (GRBs) can produce the observed flux of cosmic rays at the highest energies. However, recent studies of GRBs indicate that their redshift distribution likely follows the average star formation rate of the universe and that GRBs were more numerous at high redshifts. As a consequence, we show that photomeson production energy losses suffered by ultrahigh energy cosmic rays coming from GRBs would produce too sharp a spectral energy cutoff to be consistent with the air shower data. Futhermore, we show that cosmolgical GRBs fail to supply the energy input required to account for the cosmic ray flux above 10¹⁹ eV by a factor of 100–1000.     

Gamma-ray emitting supernova remnants as the origin of Galactic cosmic rays?

The origin of cosmic rays is one of the long-standing mysteries in physics and astrophysics. Simple arguments suggest that a scenario of supernova remnants (SNRs) in the Milky Way as the dominant sources for the cosmic ray population below the knee could work: a generic calculation indicates that these objects can provide the energy budget necessary to explain the observed flux of cosmic rays. However, this argument is based on the assumption that all sources behave in the same way, i.e. they all have the same energy budget, spectral behavior and maximum energy. In this paper, we investigate if a realistic population of SNRs is capable of producing the cosmic ray flux as it is observed below the knee. We use 21 SNRs that are well-studied from radio wavelengths up to gamma-ray energies and derive cosmic ray spectra under the assumption of hadronic emission. The cosmic ray spectra show a large variety in their energy budget, spectral behavior and maximum energy. These sources are assumed to be representative for the total class of SNRs, where we assume that about 100–200 cosmic ray emitting SNRs should be present today. Finally, we use these source spectra to simulate the cosmic ray transport from individual SNRs in the Galaxy with the GALPROP code for cosmic ray propagation. We find that the cosmic ray budget can be matched well for these sources. We conclude that gamma-ray emitting SNRs can be a representative sample of cosmic ray emitting sources. In the future, experiments like CTA and HAWC will help to distinguish hadronic from leptonic sources and to further constrain the maximum energy of the sources and contribute to producing a fully representative sample in order to further investigate the possibility of SNRs being the dominant sources of cosmic rays up to the knee.     

A beam of particles in ultrahigh-energy cosmic rays?

Three particles with energies of 36, 35, and 58 EeV arrived from one sky region were recorded by two EAS arrays during a day. The events are assumed to have been produced by the beam of particles that resulted from the interaction of cosmic rays with a relativistic shock front.     

Optical afterglows of gamma-ray bursts: a bimodal distribution?

The luminosities of the optical afterglows of gamma-ray bursts, 12 h (rest-frame time) after the trigger, show a surprising clustering, with a minority of events being at a significantly smaller luminosity. If real, this dichotomy would be a crucial clue to understand the nature of optically dark afterglows, i.e. bursts that are detected in the X-ray band, but not in the optical. We investigate this issue by studying bursts of the pre- Swift era, both detected and undetected in the optical. The limiting magnitudes of the undetected ones are used to construct the probability that a generic burst is observed down to a given magnitude limit. Then, by simulating a large number of bursts with pre-assigned characteristics, we can compare the properties of the observed optical luminosity distribution with the simulated one. Our results suggest that the hints of bimodality present in the observed distribution reflect a real bimodality: either the optical luminosity distribution of bursts is intrinsically bimodal, or there exists a population of bursts with a quite significant grey absorption, i.e. wavelength-independent extinction. This population of intrinsically weak or grey-absorbed events can be associated with dark bursts.     

Classification of BeppoSAX’s gamma-ray bursts

The BeppoSAX Catalog has been very recently published. In this paper we analyze—using the Maximum Likelihood (ML) method—the duration distribution of the 1003 GRBs listed in the catalog with duration. The ML method can identify the long and the intermediate duration groups. The short population of the bursts is identified only at a 96% significance level. MC simulation has been also applied and gives a similar significance level; 95%. However, the existence of the short bursts is not a scientific question after the Compton Gamma-Ray Observatory’s observation. Our minor result is this well-known fact that in the BeppoSAX data the short bursts are under-represented, mainly caused by the different triggering system. Our major result is the identification of the intermediate group in the BeppoSAX data. Therefore, four different satellites (CGRO, Swift, RHESSI and BeppoSAX) observed the intermediate type Gamma-Ray Burst.     

Spectrum and anisotropy of cosmic rays at TeV–PeV-energies and contribution of nearby sources

The role of nearby galactic sources, the supernova remnants, in formation of observed energy spectrum and large-scale anisotropy of high-energy cosmic rays is studied. The list of these sources is made up based on radio, X-ray and gamma-ray catalogues. The distant sources are treated statistically as ensemble of sources with random positions and ages. The source spectra are defined based on the modern theory of cosmic ray acceleration in supernova remnants while the propagation of cosmic rays in the interstellar medium is described in the frameworks of galactic diffusion model. Calculations of dipole component of anisotropy are made to reproduce the experimental procedure of “two-dimensional” anisotropy measurements. The energy dependence of particle escape time in the process of acceleration in supernova remnants and the arm structure of sources defining the significant features of anisotropy are also taken into account. The essential new trait of the model is a decreasing number of core collapse SNRs being able to accelerate cosmic rays up to the given energy, that leads to steeper total cosmic ray source spectrum in comparison with the individual source spectrum. We explained simultaneously the new cosmic ray data on the fine structure of all particle spectrum around the knee and the amplitude and direction of the dipole component of anisotropy in the wide energy range 1 TeV–1 EeV. Suggested assumptions do not look exotic, and they confirm the modern understanding of cosmic ray origin.     

A catalog of cosmic gamma-ray bursts detected by the PHEBUS instrument on the Granat observatory: October 1994–December 1996

We present the final part of the catalog of cosmic gamma-ray bursts (GRBs) observed in the PHEBUS experiment on the Granat orbiting astrophysical observatory. The first three parts of the catalog were published by Terekhov et al. (1994, 1995a) and Tkachenko et al. (1998). The fourth part contains information on 32 events recorded from October 1994 until December 1996. We give burst light curves in the energy range 100 keV to 1.6 MeV, integrated energy spectra, and information on the fluence and energy flux at the luminosity peak for energies above 100 keV. Over the entire period of its operation, the PHEBUS instrument detected 206 cosmic GRBs. The mean ?V/V_max? was 0.336±0.007. The mean hardness corresponding to the ratio of count numbers in the energy ranges 400–1000 and 100–400 keV is 0.428±0.018 for events with a duration shorter than 2 s and 0.231±0.004 for events with a duration longer than 2 s.     

Soft gamma-ray bursts from the superluminal source GRS 1915+105?

The power required to eject relativistic plasma clouds in the hard X-ray transient GRS 1915+105 is at least 100 times the luminosity of the soft-ray bursts (SGRs) that were observed by BATSE from the same region of the sky in the year 1992. We show that there are spatial, time, and spectral coincidences between GRS 1915+105 and the SGRs observed by BATSE which suggest that they are one and the same source. However, the position of the SGRs is rather uncertain and until better positions are obtained, the question on the association of GRS 1915+105 with the SGRs must remain open.     

A search for discrete gamma-ray sources of energy greater than 2×1012 eV

The results of the observations to search gamma-ray sources with the energy greater than 2×10¹² eV, which were made in Crimean Astrophysical Observatory during the years 1969–73 are presented. A technique of the detection of the EAS Cerenkov flashes was used.The quality of the data obtained is analysed. The criteria for the selection of the data free from meteorological variations are considered.It was shown that two objects, namely, Cyg X-3 and Cas -1, may be the sources of high-energy gamma quanta. It is probable that the object with the coordinates =05^h15^m, =+1° is the source of gamma-rays as well. An unidentified object Cas -1 is variable: gamma-ray flux was observed twice — in Sepember–October 1971 and in December 1972. It is possible that the flux from Cyg X-3 has a period of 4.8 hr.

---

I I , I I , - >2.10¹² . I . I , I I, I ., - -1 Cyg -3- -I . , =0515 ·=+1° -.I -1 I: I J I- - 1971 1972 . Cyg -3, , - T=4.8 .

     

Are short-time variations of the solar S-component emission identical with microwave bursts?

Extended time series (time resolution about 2–3 min) of spatially resolved observations (≫ 17 arcsec) in one dimension of solar S-component sources obtained at the Siberian Solar Radio Telescope (SSRT) at 5.2 cm wavelength allow the detection of evolutional features of the growth and decay of active regions in the solar corona. Characteristic slow flux variations with timescales of about 1–2 hours occurring during the decay phase of the radio emission in the low corona above plages and sunspots are compared with recently detected step-like flux increases on timescales of about 10–20 min followed by quasi-constant periods appearing in the initial phase of the development of active regions. Superimposed on this basic behaviour, also fluctuations at shorter timescales (or even periodic oscillations) have been observed. As it is well known from emission-model calculations, the variations of the S-component radiation can be due to variations of the magnetic field and/or changes of the energy of the radiating particles, which is basically the same emission mechanism as for microwave bursts. Since the “S-component” is originally defined by its long timescale behaviour derived from whole-Sun flux density measurements, the presently detected small-timescale features in S-component sources require either a revised definition of S-component emission or must be considered as “burst-like”.     

Can the bump in the composite spectrum of GRB 910503 be an emission line feature of gamma-ray bursts？

Appearing in the composite spectral data of BATSE, EGRET and COMPTEL for GRB 910503, there is a bump at around 1600 keV. We perform a statistical analysis on the spectral data, trying to find out if the bump could be accounted for by a blue-shifted and significantly broadened rest frame line due to the Doppler effect of an expanding fireball surface. We made an F-test and adopted previously proposed criteria. The study reveals that the criteria are well satisfied and the feature can be interpreted as the blue shifted 6.4 keV line. Prom the fit with this line taken into account, we find the Lorentz factor of this source to be P = 116-9+9 (at the 68% confident level,△x2 = 1) and the rest frame spectral peak energy to be E0,p- 2.96-0.18+0.24 keV. Although the existence of the emission line feature requires other independent tests to confirm, the analysis suggests that it is feasible to detect emission line features in the high energy range of GRB spectra when taking into account the Doppler effect of fireball expansion.     

Heliospheric modulation of cosmic rays in solar cycles 19–23

The previously developed basic theory of the heliospheric modulation of high-energy cosmic rays is generalized to lower energies. Comparison of the theory with the results of long-term observations of cosmic rays in the stratosphere carried out by the group from the Lebedev Institute of Physics inMoscow andMurmansk shows satisfactory agreement. The cosmic rays are shown to behave quite differently when even and odd solar cycles alternate. Possible causes of the anomalously high cosmic-ray intensity recorded during the last solar activity minimum are discussed.     

Diurnal anisotropy of cosmic rays during intensive solar activity for the period 2001–2014

The diurnal variation of cosmic ray intensity, based on the records of two neutron monitor stations at Athens (Greece) and Oulu (Finland) for the time period 2001 to 2014, is studied. This period covers the maximum and the descending phase of the solar cycle 23, the minimum of the solar cycles 23/24 and the ascending phase of the solar cycle 24.These two stations differ in their geographic latitude and magnetic threshold rigidity. The amplitude and phase of the diurnal anisotropy vectors have been calculated on annual and monthly basis.From our analysis it is resulted that there is a different behaviour in the characteristics of the diurnal anisotropy during the different phases of the solar cycle, depended on the solar magnetic field polarity, but also during extreme events of solar activity, such as Ground Level Enhancements and cosmic ray events, such as Forbush decreases and magnetospheric events. These results may be useful to Space Weather forecasting and especially to Biomagnetic studies.     

Solar magnetic cycle dependence in corotating modulation of galactic cosmic rays

We study the temporal evolution of cosmic ray intensity during ～27-day Carrington rotation periods applying the method of superposed epoch analysis. We discuss about the average oscillations in the galactic cosmic ray intensity, as observed by ground based neutron monitors, during the course of Carrington rotation in low solar activity conditions and in different polarity states of the heliosphere (A<0 and A>0). During minimum and decreasing phases in low solar activity conditions, we compare the oscillation in one polarity state with that observed in other polarity state in similar phases of solar activity. We find difference in the evolution and amplitude of ～27-day variation during A<0 and A>0 epoch. We also compare the average variations in cosmic ray intensity with the simultaneous variations of solar wind parameters such as solar wind speed and interplanetary magnetic field strength. From the correlation analysis between the cosmic ray intensity and the solar wind speed during the course of Carrington rotation, we find that the correlation is stronger for A>0 than A<0.     

Solar modulation of galactic cosmic rays during 19–23 solar cycles

We make a detailed analysis of cross-correlation and time-lag between monthly data of galactic cosmic rays (GCRs) intensity and different solar activity indices (e.g., sunspot number, sunspot area, green coronal Fe line and 10.7 cm solar radio flux) during 19–23 solar cycles. GCRs time-series data from Kiel neutron monitor station and solar data from the last 50 years period, covering five solar cycles (19–23), and alternating solar polarity states (i.e., five A < 0 and four A > 0) have been investigated. We find a clear asymmetry in the cross-correlation between GCRs and solar activity indicators for both odd and even-numbered solar cycles. The time-lags between GCRs and solar parameters are found different in different solar cycles as well as in the opposite polarity states (A < 0 and A > 0) within the same solar cycle. Possible explanations of the observed results are discussed in light of modulation models, including drift effects.     

Are ‘negative bursts’ due to absorption?

Optical activity near the time of two microwave negative bursts is distinctive enough to demonstrate a real relation, but does not occur simultaneously with the radio events. Such a loose association between microwave flux decreases and H prominence or filament activity is typical of a sample of 23 negative bursts. A model of a microwave occulting cloud differs from that of an H absorber: the H dark flocculus is relatively low, dense, and compact; the microwave occulter is higher in the corona, larger, and more completely ionized. The two types of absorption are not expected to be seen simultaneously, although they could be separate phases of an ejection of chromospheric material into the corona. The association of microwave decreases with H activity in some cases, and their usual non-simultaneity, is consistent with interpretation of the decrease as absorption, but it does not rule out alternative interpretations such as an intrinsic change in the emission of the microwave source.The location of decrease-associated activity has a suggestive, though not statistically significant, bias toward east limb that leads one toward an absorption interpretation, with asymmetry introduced by a tendency for ejected material to move from the leading part of an active region toward the following part, agreeing with Liszka's observation of asymmetry in the line-of-sight velocities of prominences.