A possible universal origin of hadronic cosmic rays from ultrarelativistic ejecta of bipolar supernovae

Based on the “cannonball model” for gamma-ray bursts of Dar and De Rújula, it is proposed that masses of baryonic plasma (“cannonballs”), ejected in bipolar supernova explosions in our Galaxy are the sources of hadronic Galactic cosmic rays (CRs) at all energies. The propagation of the cannonballs in the Galactic disk and halo is studied. Two mechanisms for the acceleration of the observed CRs are proposed. The first is based on ultrarelativistic shocks in the interstellar medium and could accelerate the bulk of CRs up to the “knee” energy of 4×10¹⁵ eV. The second operates with second-order Fermi acceleration within the cannonball. If the total initial energy of the ejected plasmoids in a SN explosion is 10⁵³ erg or higher, this second mechanism may explain the CR spectrum above the knee up to the highest observed energies. It is shown that together with plausible assumptions about CR propagation in the Galactic confinement volume, the observed spectral indices of the CR spectrum can be theoretically understood to first order. The model allows a natural understanding of various basic CR observations like the absence of the Greisen–Zatsepin cutoff, the anisotropy of arrival directions as function of energy and the small Galactocentric gradient of the CR density.     

Non-thermal processes in large solar flares

We analyze particle acceleration processes in large solar flares, using observations of the August, 1972, series of large events. The energetic particle populations are estimated from the hard X-ray and γ-ray emission, and from direct interplanetary particle observations. The collisional energy losses of these particles are computed as a function of height, assuming that the particles are accelerated high in the solar atmosphere and then precipitate down into denser layers. We compare the computed energy input with the flare energy output in radiation, heating, and mass ejection, and find for large proton event flares that:

1. The ～10–10² keV electrons accelerated during the flash phase constitute the bulk of the total flare energy.
2. The flare can be divided into two regions depending on whether the electron energy input goes into radiation or explosive heating. The computed energy input to the radiative quasi-equilibrium region agrees with the observed flare energy output in optical, UV, and EUV radiation.
3. The electron energy input to the explosive heating region can produce evaporation of the upper chromosphere needed to form the soft X-ray flare plasma.
4. Very intense energetic electron fluxes can provide the energy and mass for interplanetary shock wave by heating the atmospheric gas to energies sufficient to escape the solar gravitational and magnetic fields. The threshold for shock formation appears to be ～10³¹ ergs total energy in >20 keV electrons, and all of the shock energy can be supplied by electrons if their spectrum extends down to 5–10 keV.
5. High energy protons are accelerated later than the 10–10² keV electrons and most of them escape to the interplanetary medium. The energetic protons are not a significant contributor to the energization of flare phenomena. The observations are consistent with shock-wave acceleration of the protons and other nuclei, and also of electrons to relativistic energies.
6. The flare white-light continuum emission is consistent with a model of free-bound transitions in a plasma with strong non-thermal ionization produced in the lower solar chromosphere by energetic electrons. The white-light continuum is inconsistent with models of photospheric heating by the energetic particles. A threshold energy of ～5×10³⁰ ergs in >20 keV electrons is required for detectable white-light emission.

The highly efficient electron energization required in these flares suggests that the flare mechanism consists of rapid dissipation of chromospheric and coronal field-aligned or sheet currents, due to the onset of current-driven Buneman anomalous resistivity. Large proton flares then result when the energy input from accelerated electrons is sufficient to form a shock wave.     

Multiple acceleration of protons on the Sun and their free propagation to the Earth on January 20, 2005

We analyze the observations of solar protons with energies >80 MeV near the Earth and the January 20, 2005, solar flare in various ranges of the electromagnetic spectrum. Within approximately the first 30 min after their escape into interplanetary space, the solar protons with energies above 80 MeV propagated without scattering to the Earth and their time profiles were determined only by the time profile of the source on the Sun and its energy spectrum. The 80–165 MeV proton injection function was nonzero beginning at 06:43:80 UT and can be represented as the product of the temporal part, the ACS (Anticoincidence System) SPI (Spectrometer on INTEGRAL) count rate, and the energy part, a power-law proton spectrum ～E ^?4.7±0.1. Protons with energies above 165 MeV and relativistic electrons were injected, respectively, 4 and 9 min later than this time. The close correlation between high-energy solar electromagnetic emission and solar proton fluxes near the Earth is evidence for prolonged and multiple proton acceleration in solar flares. The formation of a posteruptive loop system was most likely accompanied by successive energy releases and acceleration of charged particles with various energies. Our results are in conflict with the ideas of cosmic-ray acceleration in gradual solar particle events at the shock wave driven by a coronal mass ejection.     

A massive circumstellar envelope around the type-IIn supernova SN 1995G

We model the interaction of the supernova SN 1995G with a dense circumstellar (CS) gas in a thin-shell approximation. A model fit to the observed bolometric light curve combined with data on the supernova expansion velocity gives estimates for the density, mass (≈1 M_⊙), and age (≈8 yr) of the CS envelope. The determined CS-envelope density is shown to be virtually independent of the assumed mass of the supernova envelope because of the high CS-gas density at which the forward shock wave is essentially radiative. The derived CS-envelope density is consistent with the Hα luminosity and with the presence of distinct Thomson scattering in the red wing of this line. The mass of the CS envelope together with its expansion velocity and age indicate that the CS envelope was ejected by the presupernova eight years before the supernova explosion through violent energy release (～6×10⁴⁸ erg).     

Diagnostics of the early explosion phase of a classical nova using its X-ray emission: A model for the X-ray outburst of CI Camelopardalis in 1998

We have computed a spherically symmetric model for the interaction of matter ejected during the outburst of a classical nova with the stellar wind from its optical component. This model is used to describe the intense X-ray outburst (the peak 3–20 keV flux was ～2 Crab) of the binary system CI Camelopardalis in 1998. According to our model, the stellar wind from the optical component heated by a strong shock wave produced when matter is ejected from the white dwarf as the result of a thermonuclear explosion on its surface is the emission source in the standard X-ray band. Comparison of the calculated and observed time dependences of the mean radiation temperature and luminosity of the binary system during its outburst has yielded very important characteristics of the explosion. We have been able to measure the velocity of the ejected matter immediately after the onset of the explosion for the first time: it follows from our model that the ejected matter had a velocity of ～2700 km s^?1 even on 0.1–0.5 day after the outburst onset and it flew with such a velocity for the first 1–1.5 day under an external force, possibly, the radiation pressure from the white dwarf. Subsequently, the matter probably became transparent and began to decelerate. The time dependence of the mean radiation temperature at late expansion phases has allowed us to estimate the mass of the ejected matter, ～10^?7–10^?6 M _⊙. The mass loss rate in the stellar wind required to explain the observed peak luminosity of the binary system during its outburst has been estimated to be \(\dot M\) ～ (1 ? 2) × 10^?6 M _⊙ yr^?1.     

Ultrahigh energy cosmic ray acceleration in Seyfert galactic nuclei

We propose a model for the particle acceleration to energy E≈10²¹ eV in Seyfert galactic nuclei. The model is based on the theory of active galactic nuclei by Vilkoviskij et al. (1999). The acceleration takes place in hot spots of relativistic jets, which decay in a dense stellar kernel at a distance of 1–3 pc from the center. The maximum energy and chemical composition of the accelerated particles depend on the jet magnetic-field strength. Fe nuclei acquire the largest energy, E≈8×10²⁰ eV, if the jet field strength is B≈16 G. At a field strength B～5–40 G, the nuclei with Z≥10 acquire energy E≥2×10²⁰ eV; the lighter nuclei are accelerated to E≤10²⁰ eV. In a field B～1000 G, only the particles with Z≥23 gain energy E≤10²⁰ eV. The protons are accelerated to E<4×10¹⁹ eV, and they do not fall within the energy range concerned at any field strength B. Interactions with infrared photons do not affect the accelerated-particle escape from the sources if the galactic luminosity L≤10⁴⁶ erg s^?1 and if the angle between the normal to the galactic plane and the line of sight is sufficiently small, i.e., if the galactic-disk axial ratio is comparatively large. The particles do not lose their energy through magnetodrift radiation if their deflection from the jet axis does not exceed 0.03–0.04 pc at a distance R≈40–50 pc from the center. The synchrotron losses are small, because the magnetic field frozen in the galactic wind at R≤40–50 pc is directed (as in the jet) predominantly along the motion. If this model is correct, then the detected protons are nuclear fragments or are accelerated in other sources. The jet magnetic fields can be estimated by using the cosmic-ray energy spectrum and chemical composition.     

Cosmic ray acceleration in hypernova remnants expanding in wind bubbles of Wolf-Rayet progenitor stars

The determination of the origin of cosmic rays with observed energies in excess of 10¹⁷ eV that exceed the expected energies of cosmic rays accelerated by supernova remnants in the galaxy is a pressing problem in modern astrophysics. Hypernova remnants are one of the possible galactic sources of cosmic rays with energies of up to 10¹⁹ eV. Hypernovae constitute a class of extremely powerful supernova explosions, whose supposed progenitors are massive Wolf-Rayet stars. We analyze the special aspects of acceleration of cosmic rays in hypernova remnants that expand in wind bubbles of Wolf-Rayet progenitor stars. We show that these cosmic rays may attain maximum energies of 10¹⁸ eV even with a relatively conservative choice of acceleration parameters and account for tens of percent of the total cosmic ray flux observed in the vicinity of the earth in the energy range of 10¹⁶–10¹⁸ eV if the galactic hypernova explosion rate in the modern epoch reaches ? _S ～ 10^?4 year^?1.     

WD + He star systems as the progenitors of Type Ia supernovae and their surviving companion stars

Employing Eggleton’s stellar evolution code with the optically thick wind assumption, we have systematically studied the WD + He star channel of Type Ia supernovae (SNe Ia), in which a carbon–oxygen WD accretes material from a He main-sequence star or a He subgiant to increase its mass to the Chandrasekhar mass. We mapped out the parameter spaces for producing SNe Ia. According to a detailed binary population synthesis approach, we find that the Galactic SN Ia birthrate from this channel is ～0.3×10^?3 yr^?1, and that this channel can produce SNe Ia with short delay times (～45–140 Myr). We also find that the surviving companion stars in this channel have a high spatial velocity (>400 km/s) after the SN explosion, which could be an alternative origin for hypervelocity stars (HVSs), especially for HVSs such as US 708.     

Energetic proton irradiation history of the howardite parent body regolith and implications for ancient solar activity

Abstract— Previous studies have shown that the Kapoeta howardite, as well as several other meteorites, contains excess concentrations of cosmogenic Ne in the darkened, solar-irradiated phase compared to the light, non-irradiated phase. The two explanations offered for the nuclear production of these Ne excesses in the parent body regolith are either from galactic cosmic-ray proton (GCR) irradiation or from a greatly enhanced flux of energetic solar “cosmic-ray” protons (SCR), as compared to the recent solar flux. Combining new isotopic data we obtained on acid-etched, separated feldspar from Kapoeta light and dark phases with literature data, we show that the cosmogenic ²¹Ne/²²Ne ratio of light phase feldspar (0.80) is consistent with only GCR irradiation in space for ～3 Ma. However, the ²¹Ne/²²Ne ratio (0.68) derived for irradiation of dark phase feldspar in the Kapoeta regolith indicates that cosmogenic Ne was produced in roughly equal proportions from galactic and solar protons. Considering a simple model of an immature Kapoeta parent body regolith, the duration of this early galactic exposure was only ～3–6 Ma, which would be an upper limit to the solar exposure time of individual grains. Concentrations of cosmogenic ²¹Ne in pyroxene separates and of cosmogenic ¹²⁶Xe in both feldspar and pyroxene are consistent with this interpretation. The near-surface irradiation time of individual grains in the Kapoeta regolith probably varied considerably due to regolith mixing to an average GCR irradiation depth of ～10 cm. Because of the very different depth scales for production of solar ～Fe tracks, SCR Ne, and GCR Ne, the actual regolith exposure times for average grains probably differed correspondingly. However, both the SCR ²¹Ne and solar track ages appear to be longer because of enhanced production by early solar activity. The SCR/GCR production ratio of ²¹Ne inferred from the Kapoeta data is larger by a at least a factor of 10 and possibly as much as a factor of ～50 compared to recent solar particle fluxes. Thus, this study indicates that our early Sun was much more active and emitted a substantially higher flux of energetic (>10 MeV/nucleon) protons.     

Type-Ia supernovae in dense circumstellar gas

We propose a model for the bolometric light curve of a type-Ia supernova (SN Ia) that explodes in a dense circumstellar (CS) envelope. Our modeling of the light curves for SN 2002ic and SN 1997cy shows that the densities of the CS envelopes around both supernovae at a radius of ～7×10¹⁵ cm are similar, while the characteristic ejection time for this envelope around SN 1997cy does not exceed 600 yr. We analyze two possible evolutionary scenarios that could lead to the explosion of a SN Ia inside a dense C S hydrogen envelope: accretion onto a CO white dwarf in a symbiotic binary and the evolution of a single star with an initial mass of about 8M_⊙. If the hypothesis of a SN Ia explosion in a dense CS envelope is correct for SN 2002ic and SN 1997cy, then we must assume that the the rapid loss of the red-supergiant envelope in several hundred years and the subsequent explosion of the CO white dwarf are synchronized by some physical mechanism. This mechanism may be related to the contraction of the white dwarf as it approaches the Chandrasekhar limit. We show that the formation of a (super-)Chandrasekhar mass due to the merger of a CO white dwarf and the CO core of a red supergiant followed by a supernovae explosion is unlikely, since this mechanism does not provide the required synchronization of the rapid mass loss and the explosion.     

4U 0115+63 from RXTE and INTEGRAL data: Pulse profile and cyclotron line energy

We analyze the observations of the transient X-ray pulsar 4U 0115+63 with the RXTE and INTEGRAL observatories in a wide X-ray (3–100 keV) energy band during its intense outbursts in 1999 and 2004. The energy of the fundamental harmonic of the cyclotron resonance absorption line near the maximum of the X-ray flux from the source (luminosity range 5 × 10³⁷–2 × 10³⁸ erg s^?1) is ～11 keV. When the pulsar luminosity falls below ～5 × 10³⁷ erg s^?1, the energy of the fundamental harmonic is displaced sharply toward the high energies, up to ～16 keV. Under the assumption of a dipole magnetic field configuration, this change in cyclotron harmonic energy corresponds to a decrease in the height of the emitting region by ～2 km, while other spectral parameters, in particular, the cutoff energy, remain essentially constant. At a luminosity ～7 × 10³⁷ erg s^?1, four almost equidistant cyclotron line harmonics are recorded in the spectrum. This suggests that either the region where the emission originates is compact or the emergent spectrum from different (in height) segments of the accretion column is uniform. We have found significant pulse profile variations with energy, luminosity, and time. In particular, we show that the profile variations from pulse to pulse are not reduced to a simple modulation of the accretion rate specified by external conditions.     

Estimation of 56Ni Needed in the Explosion Process of the Super Luminous Supernova ASASSN-15lh

ASASSN-15lh is a super luminous supernova, whose light curve is similar to that of the type Ia supernova (SN Ia). Since the luminosity of SN Ia is directly related to the decay of ⁵⁶Ni, in this paper, we consider the de-excitation energy of the new nuclei, and calculate the energy generated by the decay of ⁵⁶Ni in the explosive environment of ASASSN-15lh. The calculated mass of ⁵⁶Ni needed by the ASASSN-15lh explosion is 31.32 M_⊙. This result agrees with the estimation of the mass of ⁵⁶Ni ≥ 30 M_⊙ derived from the observed light curve of ASASSN-15lh. No agreement has reached for the explosion mechanism of supernova ASASSN-15lh so far. The calculation in this paper provides a reference for the further study on the progenitor and explosion mechanism of the supernova ASASSN-15lh.     

On timing and spectral characteristics of the X-ray pulsar 4U 0115+63: Evolution of the pulsation period and the cyclotron line energy

An overview of the results of observations for the transient X-ray pulsar 4U 0115+63, amember of a binary system with a Be star, since its discovery to the present day (～40 years) based on data from more than dozen observatories and instruments is presented. An overall light curve and the history of change in the spin frequency of the neutron star over the entire history of its observations, which also includes the results of recent measurements made by the INTEGRAL observatory during the 2004, 2008, and 2011 outbursts, are provided. The source’s energy spectra have also been constructed from the INTEGRAL data obtained during the 2011 outburst for a dynamic range of its luminosities 10³⁷?7 × 10³⁷ erg s^?1. We show that apart from the fundamental harmonic of the cyclotron absorption line at energy～11 keV, its four higher harmonics at energies ?24, 35.6, 48.8, and 60.7 keV are detected in the spectrum. We have performed a detailed analysis of the source’s spectra in the 4–28 keV energy band based on all of the available RXTE archival data obtained during bright outbursts in 1995–2011. We have confirmed that modifying the source’s continuum model can lead to the disappearance of the observed anticorrelation between the energy of the fundamental harmonic of the cyclotron absorption line and the source’s luminosity. Thus, the question about the evolution of the cyclotron absorption line energy with the luminosity of the X-ray pulsar 4U 0115+63 remains open and a physically justified radiation model for X-ray pulsars is needed to answer it.     

Broadband spectrum of the total X-ray emission from the galaxy M31

We present the results of measurements of the total X-ray flux from the Andromeda galaxy (M31) in the 3-100 keV band based on data from the RXTE/PCA, INTEGRAL/ISGRI, and SWIFT/BAT space experiments. We show that the total emission from the galaxy has a multicomponent spectrum whose main characteristics are specified by binaries emitting in the optically thick and optically thin regimes. The galaxy’s luminosity at energies 20–100 keV gives about 6% of its total luminosity in the 3–100 keV band. The emissivity of the stellar population in M31 is L _{2–20 keV} ～ 1.1 × 10²⁹ erg s^?1 M _⊙ ^?1 in the 2–20 keV band and L _{20–100 keV} ～ 8 × 10²⁷ erg s^?1 M _⊙ ^?1 in the 20–100 keV band. Since low-mass X-ray binaries at high luminosities pass into a soft state with a small fraction of hard X-ray emission, the detection of individual hard X-ray sources in M31 requires a sensitivity that is tens of times better (up to 10^?13 erg s^?1 cm^?2) than is needed to detect the total hard X-ray emission from the entire galaxy. Allowance for the contribution from the hard spectral component of the galaxy changes the galaxy’s effective Compton temperature approximately by a factor of 2, from ～1.1 to ～2.1 keV.     

Gamma-ray emission expected from Kepler’s SNR

Nonlinear kinetic theory of cosmic ray (CR) acceleration in supernova remnants (SNRs) is used to investigate the properties of Kepler’s SNR and, in particular, to predict the γ-eay spectrum expected from this SNR. Observations of the nonthermal radio and X-ray emission spectra as well as theoretical constraints for the total supernova (SN) explosion energy E _sn are used to constrain the astronomical and particle acceleration parameters of the system. Under the assumption that Kepler’s SN is a type Ia SN we determine for any given explosion energy E _sn and source distance d the mass density of the ambient interstellar medium (ISM) from a fit to the observed SNR size and expansion speed. This makes it possible to make predictions for the expected γ-eay flux. Exploring the expected distance range we find that for a typical explosion energy E _sn=10⁵¹ erg the expected energy flux of TeV γ-rays varies from 2×10⁻¹¹ to 10⁻¹³ erg/(cm² s) when the distance changes from d=3.4 kpc to 7 kpc. In all cases the γ-eay emission is dominated by π ⁰-decay γ-rays due to nuclear CRs. Therefore Kepler’s SNR represents a very promising target for instruments like H.E.S.S., CANGAROO and GLAST. A non-detection of γ-rays would mean that the actual source distance is larger than 7 kpc.     

Pulsars and the origin of cosmic rays

The capabilities and limitations of pulsars as sources of cosmic rays are reviewed in the light of experimental observations. Pulsars can supply the cosmic ray power if they have rotational velocities in excess of 700 rad s^?1 at birth. Though this is theoretically possible, there is no experimental proof for the same. Pulsars can accelerate particles to the highest energies of 10²⁰ eV, but in general, the spectra on simple considerations, turn out to be flatter than the observed cosmic ray spectrum. At the highest energies, absorption processes due to fragmentation and photodisintegration dominate for heavy nuclei. The existence of a steady flux of cosmic rays of energy greater than 10¹⁷ eV demands acceleration of particles to last over fifty years, the time interval between supernovae outbursts, whereas the expected period of activity is less than a few years. Finally, the problem of anisotropy with relevance to pulsars as sources and the possibility of observing pulsar accelerated particles from galactic clusters is considered.     

Diffuse cosmic gamma-ray background in the 28 ke V-4.1 MeV range from Kosmos 461 observations

Diffuse cosmic background and atmospheric gamma-radiation in the range 28 keV-4.1 MeV were studied with a scintillation spectrometer on board of the Kosmos 461 satellite. Separation of the cosmic and atmospheric components was made possible through a reliable determination of the geomagnetic dependences of albedo gamma-radiation: The spectrum of diffuse background in the energy range covered cannot be fitted with a common law. At energies below 400 keV the spectrum follows a power-law $$I = (5.6 \pm 0.5) \times 10^{ - 3} E^{ - (2.80 \pm 0.05)} cm^{ - 2} s^{ - 1} sr^{ - 1} MeV^{ - 1} .$$ Starting from 400 keV, this power-law breaks down; the spectrum revealing a clearly pronounced shoulder. Extrapolation of the power-law spectrum to higher energies shows that the gamma-ray component responsible for the change in the shape of the spectrum is quite strong, becoming predominant in the diffuse background in the range 1–100 MeV. The intensity of excess radiation is maximum in the region of 700–800 keV reaching ～1.8×10^?2 cm^?2s^?1sr^?1 MeV^?1. The shape of the high energy component spectrum of the diffuse background constructed using the data of Kosmos 461 and SAS-2 is in agreement with the hypotheses of the cosmological origin of the radiation.     

Efficient cosmic-ray production in supernova remnant G347.3–0.5

Broad-band models of the bright NW limb of G347.3–0.5 give convincing evidence that the forward shock of this supernova remnant is accelerating cosmic rays efficiently, placing >25% of the shock kinetic energy flux into relativistic ions. Despite this high efficiency, the maximum electron and proton energies are well below the observed ‘knee’ at ∼10¹⁵ eV in the Galactic cosmic-ray spectrum.     

On time variations of the intensity of galactic cosmic rays for the recent billion years from the data on exposure ages of iron meteorites

To ascertain probable variations of the intensity of galactic cosmic rays (GCR) for the recent billion years, the distribution of exposure ages T of iron meteorites has been analyzed. We considered all ~80 values of ages from the data by Voshage and Feldmann (1979), Voshage et al. (1983), and Voshage (1984), as well as a set of values obtained from the correction for eliminating the meteorites formed in a single collision. To correct the data, the Akaike information criterion was used. For the distributions of the phase values Ph = T/t–int(T/t), the dependence of the criterion χ ² on the presumable period t in the exposure age variations was analyzed. For t ~ 400–500 Myr and, partly, for t ~ 150 Myr, the significant deviations of this criterion from the corresponding mean values were found. To clear up the influence of the GCR intensity variations on the age distribution, the numerical models were calculated with an account of the set of ages randomly distributed in the interval of 0–1000 Myr with the presumptive mean lifetime of iron meteorites in outer space τ = 700 Myr. It has been ascertained that, for variations with a period of t = 450 Myr, the distribution of exposure ages of the model set is similar to that found for iron meteorites. The obtained data suggest that the GCR intensity variations with a period of approximately 400–500 Myr have probably existed during the recent billion years. These variations may be caused by periodic passages of the Solar System through spiral arms of the Galaxy. It has been shown that the earlier discussed changes in the GCR intensity with a period of ~150 Myr (Shaviv, 2002; 2003; Scherer et al., 2006) are less defined.