Search for point-like sources of cosmic rays with energies above 1018.5 eV in the HiRes-I monocular data set

We report the results of a search for point-like deviations from isotropy in the arrival directions of ultra-high energy cosmic rays in the northern hemisphere. In the monocular data set collected by the High-Resolution Fly’s Eye, consisting of 1525 events with energy exceeding 10^18.5 eV, we find no evidence for point-like excesses. We place a 90% c.l. upper limit of 0.8 hadronic cosmic rays/km² yr on the flux from such sources for the northern hemisphere and place tighter limits as a function of position in the sky.     

Cosmic ray acceleration to very high energy through the non-linear amplification by cosmic rays of the seed magnetic field

The maximum energy for cosmic ray acceleration at supernova shock fronts is usually thought to be limited to around 10¹⁴–10¹⁵ eV by the size of the shock and the time for which it propagates at high velocity. We show that the magnetic field can be amplified non-linearly by the cosmic rays to many times the pre-shock value, thus increasing the acceleration rate and facilitating acceleration to energies well above 10¹⁵ eV. A supernova remnant expanding into a uniform circumstellar medium may accelerate protons to 10¹⁷ eV and heavy ions, with charge Ze , to Z ×10¹⁷ eV. Expansion into a pre-existing stellar wind may increase the maximum cosmic ray energy by a further factor of 10.     

Production of neutrons, neutrinos and gamma-rays by a very fast pulsar in the Galactic Centre region

We consider the possibility that the excess of cosmic rays near ∼10¹⁸ eV, reported by the AGASA and SUGAR groups from the direction of the Galactic Centre, is caused by a young, very fast pulsar in the high-density medium. The pulsar accelerates iron nuclei to energies ∼10²⁰ eV, as postulated by the Galactic models for the origin of the highest-energy cosmic rays. The iron nuclei, about 1 yr after pulsar formation, leave the supernova envelope without energy losses and diffuse through the dense central region of the Galaxy. Some of them collide with the background matter creating neutrons (from disintegration of Fe), neutrinos and gamma-rays (in inelastic collisions). We suggest that neutrons produced at a specific time after the pulsar formation are responsible for the observed excess of cosmic rays at ∼10¹⁸ eV. From normalization of the calculated neutron flux to the one observed in the cosmic ray excess, we predict the neutrino and gamma-ray fluxes. It has been found that the 1 km² neutrino detector of the IceCube type should detect from a few up to several events per year from the Galactic Centre, depending on the parameters of the considered model. Moreover, future systems of Cherenkov telescopes (CANGAROO III, HESS, VERITAS) should be able to observe  1–10 TeV  gamma-rays from the Galactic Centre if the pulsar was created inside a huge molecular cloud about  3–10×10³ yr  ago.     

Estimation prospects of the source number density of ultra-high-energy cosmic rays

We discuss the possibility of accurately estimating the source number density of ultra-high-energy cosmic rays (UHECRs) using small-scale anisotropy in their arrival distribution. The arrival distribution has information on their source and source distribution. We calculate the propagation of UHE protons in a structured extragalactic magnetic field (EGMF) and simulate their arrival distribution at the Earth using our previously developed method. The source number density that can best reproduce observational results by Akeno Giant Air Shower Array is estimated at about 10⁻⁵ Mpc⁻³ in a simple source model. Despite having large uncertainties of about one order of magnitude, due to small number of observed events in current status, we find that more detection of UHECRs in the Auger era can sufficiently decrease this so that the source number density can be more robustly estimated. Two hundred event observation above 4 × 10¹⁹ eV in a hemisphere can discriminate between 10⁻⁵ and 10⁻⁶ Mpc⁻³. Number of events to discriminate between 10⁻⁴ and 10⁻⁵ Mpc⁻³ is dependent on EGMF strength. We also discuss the same in another source model in this paper.     

Cosmic rays from remnants of quasars?

Considerations of the collision losses for protons traversing the 2.7-K blackbody microwave radiation field have led to the conclusion that the highest energy cosmic rays, those observed at ≥10²⁰ eV, must come from sources within the present epoch. In light of this constraint, it is here suggested that these particles may be accelerated near the event horizons of spinning supermassive black holes associated with presently inactive quasar remnants. The required emf is generated by the black hole induced rotation of externally supplied magnetic field lines threading the horizon. Producing the observed flux of the highest energy cosmic rays would constitute a negligible drain on the black hole dynamo. Observations with upcoming air shower arrays and space missions may lead to the identification of candidate dormant galaxies that harbour such black holes. Although the highest energy events observed so far are accounted for within the context of this scenario, a spectral upper bound at ∼10²¹ eV is expected since the acceleration to higher energies appears to be precluded, on general grounds.     

Extragalactic UHE proton spectrum and prediction for iron-nuclei flux at 10–10 GeV

We investigate the problem of transition from galactic cosmic rays to extragalactic ultra-high energy cosmic rays. Using the model for extragalactic ultra-high energy cosmic rays and observed all-particle cosmic ray spectrum, we calculate the galactic spectrum of iron nuclei in the energy range 10⁸–10⁹ GeV. The flux and spectrum predicted at lower energies agree well with the KASCADE data. The transition from galactic to extragalactic cosmic rays is distinctly seen in spectra of protons and iron nuclei, when they are measured separately. The shape of the predicted iron spectrum agrees with the Hall diffusion.     

Implications of particle acceleration in active galactic nuclei for cosmic rays and high energy neutrino astronomy

We consider the production of high energy neutrinos and cosmic rays in radio-quiet active galactic nuclei (AGN) or in the central regions of radio-loud AGN. We use a model in which acceleration of protons takes place at a shock in an accretion flow onto a supermassive black hole, and follow the cascade that results from interactions of the accelerated protons in the AGN environment. We use our results to estimate the diffuse high energy neutrino intensity and cosmic ray intensity due to AGN. We discuss our results in the context of high energy neutrino telescopes under construction, and measurements of the cosmic ray composition in the region of the “knee” in the energy spectrum at 10⁷ GeV.     

Non-linear momentum diffusion of heliospheric cosmic rays

Besides parallel and perpendicular spatial diffusion, momentum diffusion can be seen as the third important process of cosmic ray transport. In this paper, the recently derived weakly non-linear theory is applied for a simple quasi-magnetostatic composite model to determine the momentum diffusion coefficient. It is demonstrated that non-linear effects are essential and cannot be neglected. Therefore, the weakly non-linear transport theory has to be preferred over the traditional quasi-linear approach. Within this improved theory, we find for the rigidity dependence of the momentum diffusion coefficient   A ∼ R ^1.4  for relativistic and   A ∼ R ^0.4  for non-relativistic cosmic rays.     

Cosmic rays at ultra high energies (Neutrinos!)

Resonant photo-pion production with the cosmic microwave background predicts a suppression of extragalactic protons above the famous Greisen–Zatsepin–Kuzmin cutoff at about E_GZK ≈ 5 × 10¹⁰ GeV. Current cosmic ray data measured by the AGASA and HiRes Collaborations do not unambiguously confirm the GZK cutoff and leave a window for speculations about the origin and chemical composition of the highest energy cosmic rays. In this work we analyze the possibility of strongly interacting neutrino primaries and derive model-independent quantitative requirements on the neutrino–nucleon inelastic cross section for a viable explanation of the cosmic ray data. Search results on weakly interacting cosmic particles from the AGASA and RICE experiments are taken into account simultaneously. Using a flexible parameterization of the inelastic neutrino–nucleon cross section we find that a combined fit of the data does not favor the Standard Model neutrino–nucleon inelastic cross section, but requires, at 90% confidence level, a steep increase within one energy decade around E_GZK by four orders of magnitude. We illustrate such an enhancement within some extensions of the Standard Model. The impact of new cosmic ray data or cosmic neutrino search results on this scenario, notably from the Pierre Auger Observatory soon, can be immediately evaluated within our approach.     

Cosmic rays and the primordial gas

One of the most-outstanding problems in the gravitational collapse scenario of early structure formation is the cooling of primordial gas to allow for small-mass objects to form. As the neutral primordial gas is a poor radiator at temperatures   T ≤ 10⁴ K  , molecular hydrogen is needed for further cooling down to temperatures   T ∼ 100 K  . The formation of molecular hydrogen is catalyzed by the presence of free electrons, which could be provided by the ionization due to an early population of cosmic rays (CRs). In order to investigate this possibility, we developed a code to study the effects of ionizing CRs on the thermal and chemical evolution of primordial gas. We found that CRs can provide enough free electrons needed for the formation of molecular hydrogen, and therefore can increase the cooling ability of such primordial gas under following conditions. A dissociating photon flux with   F < 10⁻¹⁸ erg cm⁻² Hz⁻¹ s⁻¹  , initial temperature of the gas  ∼10³ K  , total gas number densities   n ≥ 1 cm⁻³  , and cosmic-ray sources with     .     

Ultra high energy cosmic rays and the Galactic halo

An analysis has been made of the fraction of ultra high energy cosmic rays (above 10¹⁸ eV) which could be due to processes involved in two possible ‘Models’. The first is the Giant Magnetic Halo Model and the second is the Dark Matter Halo Model. We find that the former, in which heavy nuclei are trapped in a giant halo, fails for energies above about 3 × 10¹⁹ eV. For the Dark Matter Halo Model, in which relic particles follow the “conventional” dark matter and whose decays give ultra high energy cosmic rays, the predicted anisotropies are much higher than those observed. The lack of observation of a finite flux from the Andromeda Galaxy means that the conclusion is insensitive to the spatial scale size of the assumed halo distribution. It is concluded that less than 10% of the ultra high energy cosmic rays come from relic particles in the Galactic halo.     

Cosmic and Galactic neutrino backgrounds from thermonuclear sources

We estimate energy spectra and fluxes at the Earth’s surface of the cosmic and Galactic neutrino backgrounds produced by thermonuclear reactions in stars. The extra-galactic component is obtained by combining the most recent estimates of the cosmic star formation history and the stellar initial mass function with accurate theoretical predictions of the neutrino yields all over the thermonuclear lifetime of stars of different masses. Models of the structure and evolution of the Milky Way are used to derive maps of the expected flux generated by Galactic sources as a function of sky direction. The predicted neutrino backgrounds depend only slightly on model parameters. In the relevant 50 keV–10 MeV window, the total flux of cosmic neutrinos ranges between 20 and 65 cm⁻² s⁻¹. Neutrinos reaching the Earth today have been typically emitted at redshift z2. Their energy spectrum peaks at E0.1–0.3 MeV. The energy and entropy densities of the cosmic background are negligible with respect to the thermal contribution of relic neutrinos originated in the early universe. In every sky direction, the cosmic background is outnumbered by the Galactic one, whose integrated flux amounts to 300–1000 cm⁻² s⁻¹. The emission from stars in the Galactic disk contributes more than 95% of the signal.     

Stability of C60 and C70 fullerenes toward corpuscular and γ radiation

The stability of C₆₀ and C₇₀ fullerenes in the interstellar medium deposited on dust surface or embedded in meteorites and comets has been simulated with γ irradiation and with He⁺ ion bombardment. It is shown by vibrational spectroscopy that a γ radiation dose of 2.6 MGy (1 Gy = 1 joule absorbed energy per kilogram) causes partial oligomerization of both C₆₀ and C₇₀ fullerenes. Oligomers are made by fullerene cages chemically connected each other which can yield back free fullerenes by a thermal treatment. The amount of irreversibly polymerized fullerenes caused by 2.6 MGy as deduced as the toluene insoluble fraction has been determined as 1.7 and 15 per cent by weight, respectively, for C₆₀ and C₇₀ fullerene. The radiation dose generated by radionuclides decay and expected to be delivered to fullerenes buried at a depth of more than 20 m in comets and meteorites is about 3 MGy per 10⁹ yr. Since fullerenes are by far resistant to such radiation dose they can survive for at least some billion years inside comets and meteorites and in fact have been detected inside certain carbonaceous chondrites. On the other hand, the direct exposure of fullerenes to cosmic rays for instance when they are adsorbed or deposited on the surface of carbon dust corresponds to the delivery of a radiation dose comprised between 30 and 65 MGy per 10⁹ yr. Experimental bombardment of both C₆₀ and C₇₀ fullerenes for instance with He⁺ ions has shown that the complete amorphization occurs at about 250 MGy. Thus in ∼4 Gyr exposure to cosmic rays it is expected a complete amorphization.     

Contribution of nuclei accelerated by gamma-ray pulsars to cosmic rays in the Galaxy

We consider the galactic population of gamma-ray pulsars as possible sources of cosmic rays at and just above the “knee” in the observed cosmic ray spectrum at 10¹⁵–10¹⁶ eV. We suggest that iron nuclei may be accelerated in the outer gaps of pulsars, and then suffer partial photo-disintegration in the non-thermal radiation fields of the outer gaps. As a result, protons, neutrons, and surviving heavier nuclei are injected into the expanding supernova remnant. We compute the spectra of nuclei escaping from supernova remnants into the interstellar medium, taking into account the observed population of radio pulsars.

Our calculations, which include a realistic model for acceleration and propagation of nuclei in pulsar magnetospheres and supernova remnants, predict that heavy nuclei accelerated directly by gamma-ray pulsars could contribute about 20% of the observed cosmic rays in the knee region. Such a contribution of heavy nuclei to the cosmic ray spectrum at the knee can significantly increase the average value of lnA with increasing energy as is suggested by recent observations.     

On the interaction between cosmic rays and dark matter molecular clouds in the Milky Way – I. Basic considerations

We explore some basic observational consequences of assuming that the dark matter in the Milky Way consists mainly of molecular clouds, and that cosmic rays can penetrate these clouds. In a favoured model of the clouds, this penetration would have the following consequences, all of which agree with observation.
(i) Cosmic ray nuclei would be fragmented when they enter a cloud, giving them a lifetime in the Galaxy of ∼10¹⁵ s (for relativistic nuclei).
(ii) Pionic γ -rays emitted by the clouds, after proton–proton (pp) collisions, would have a diffuse flux in the Galactic plane comparable to the flux from known sources for photon energies ≳1 GeV .
(iii) The heat input into the clouds from cosmic rays would be re-radiated mainly in the far-infrared. The resulting radiation background agrees, in both intensity and spectrum in different directions, with a known excess in the far‐infrared background of the galaxy over emission by warm dust.     

Decay of the vacuum energy into cosmic microwave background photons

We examine the possibility of the decay of the vacuum energy into a homogeneous distribution of a thermalized cosmic microwave background (CMB), which is characteristic of an adiabatic vacuum energy decay into photons. It is shown that observations of the primordial density fluctuation spectrum, obtained from CMB and galaxy distribution data, restrict the possible decay rate. When photon creation due to an adiabatic vacuum energy decay takes place, the standard linear temperature dependence   T ( z ) = T ₀(1 + z )  is modified, where T ₀ is the present CMB temperature, and can be parametrized by a modified CMB temperature dependence     . From the observed CMB and galaxy distribution data, a strong limit on the maximum value of the decay rate is obtained by placing a maximum value  β_max≃ 3.4 × 10⁻³  on the β parameter.     

Gamma ray signatures of ultra high energy cosmic ray accelerators: electromagnetic cascade versus synchrotron radiation of secondary electrons

We discuss the possibility of observing ultra high energy cosmic ray sources in high energy gamma rays. Protons propagating away from their accelerators produce secondary electrons during interactions with cosmic microwave background photons. These electrons start an electromagnetic cascade that results in a broad band gamma ray emission. We show that in a magnetized Universe (B≳10⁻¹² G) such emission is likely to be too extended to be detected above the diffuse background. A more promising possibility comes from the detection of synchrotron photons from the extremely energetic secondary electrons. Although this emission is produced in a rather extended region of size ∼10 Mpc, it is expected to be point-like and detectable at GeV energies if the intergalactic magnetic field is at the nanogauss level.      

Direct observation of cosmic strings via their strong gravitational lensing effect – I. Predictions for high-resolution imaging surveys

We use present theoretical estimates for the density of long cosmic strings to predict the number of strong gravitational lensing events in astronomical imaging surveys as a function of the angular resolution and survey area. We show that angular resolution is the single most important factor, and that interesting limits on the dimensionless string tension   G μ/ c ²  can be obtained by existing and planned surveys. At the resolution of the Hubble Space Telescope ( HST ) (0.14 arcsec), it is sufficient to survey of the order of a few square degrees – well within reach of the current HST archive – to probe the regime   G μ/ c ²∼ 10⁻⁷  . If lensing by cosmic strings is not detected, such a survey would improve the limit on the string tension by a factor of two over that available from the cosmic microwave background. Future high resolution imaging surveys, covering a few hundred square degrees or more, either from space in the optical or from large-format radio telescopes on the ground, would be able to further lower this limit to   G μ/ c ²∼ 10⁻⁸  . These limits will not be improved significantly by increasing the solid angle of the survey.     

Description of Solar Cosmic Ray Propagation in the Interplanetary Medium on the Basis of the Kinetic Equation

The propagation of solar cosmic rays in interplanetary space is considered based on the kinetic equation. The expression for cosmic ray density under instantaneous particle injection by a point-like source is obtained. The set of a differential equation system for harmonics of cosmic ray distribution function is obtained starting from the kinetic equation. The cosmic ray transport equation, taking into account the presence of the second harmonic of particle angular distribution, is derived and the solution of this equation is obtained.     

Statistical distribution of gravitational-lensing excursion angles: winding ways to us from the deep Universe

We investigate statistical distributions of differences in gravitational-lensing deflections between two light rays, the so-called lensing excursion angles. A probability distribution function of the lensing excursion angles, which plays a key role in estimates of lensing effects on angular clustering of objects (such as galaxies, quasi-stellar objects and also the cosmic microwave background temperature map), is known to consist of two components: a Gaussian core and an exponential tail. We use numerical gravitational-lensing experiments in a ΛCDM cosmology for quantifying these two components. We especially focus on the physical processes responsible for generating these two components. We develop a simple empirical model for the exponential tail which allows us to explore its origin. We find that the tail is generated by the coherent lensing scatter by massive haloes with   M > 10¹⁴  h ⁻¹ M_⊙  at   z < 1  and that its exponential shape arises due to the exponential cut-off of the halo mass function at that mass range. On scales larger than 1 arcmin, the tail does not have a practical influence on the lensing effects on the angular clustering. Our model predicts that the coherent scatter may have non-negligible effects on angular clustering at subarcminute scales.