Cosmic thermal neutrino background: Can it be detected?

The detection of the Cosmic Thermal Neutrino Background (CTNB) would provide the “cleanest” evidence for the hot big bang model of the early Universe. I discuss some recent thoughts on the possibility of detecting the CTNB (especially if neutrinos have a small mass of ~ few eV) by looking for certain CTNB-induced features in the extremely high energy (E ≳ 10²⁰ eV) cosmic neutrino spectrum that may become measurable in the future by some of the large-area extensive air-shower detectors being built for detecting extremely high energy cosmic rays. NAS/NRC Senior Research Associate on sabbatical leave from Indian Institute of Astrophysics, Bangalore 560 034, India.     

Coherent Cherenkov radio pulses from hadronic showers up to EeV energies

The Cherenkov radio pulse emitted by hadronic showers of energies in the EeV range in ice is calculated for the first time using full three dimensional simulations of both shower development and the coherent radio pulse emitted as the excess charge develops in the shower. A Monte Carlo, ZHAireS, has been developed for this purpose combining the high energy hadronic interaction capabilities of AIRES, and the dense media propagation capabilities of TIERRAS, with the precise low energy tracking and specific algorithms developed to calculate the radio emission in ZHS. A thinning technique is implemented to allow the simulation of radio pulses induced by showers up to 10 EeV in ice. The code is validated comparing the results for electromagnetic and hadronic showers to those obtained with GEANT4 and ZHS codes. The contribution to the pulse of other shower particles in addition to electrons and positrons, mainly protons, pions and muons, is found to be below 3% for 10 PeV and above proton induced showers. The characteristics of hadronic showers and the corresponding Cherenkov frequency spectra are compared with those from purely electromagnetic showers. The dependence of the spectra on shower energy and high-energy hadronic model is addressed and parameterizations for the radio emission in hadronic showers in ice are given for practical applications.     

Monte Carlo simulations of radio pulses in atmospheric showers using ZHAireS

We present predictions for the radio pulses emitted by extensive air showers using ZHAireS, an AIRES-based Monte Carlo code that takes into account the full complexity of ultra-high energy cosmic-ray induced shower development in the atmosphere, and allows the calculation of the electric field in both the time and frequency domains. We do not presuppose any emission mechanism and our results are compatible with a superposition of geomagnetic and charge excess radio emission effects. We investigate the polarization of the electric field as well as the effects of the refractive index n and shower geometry on the radio pulses. We show that geometry, coupled to the relativistic effects that appear when using a realistic refractive index n > 1, play a prominent role on the radio emission of air showers.     

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.     

Thermo-acoustic sound generation in the interaction of pulsed proton and laser beams with a water target

The generation of hydrodynamic radiation in interactions of pulsed proton and laser beams with matter is explored. The beams were directed into a water target and the resulting acoustic signals were recorded with pressure sensitive sensors. Measurements were performed with varying pulse energies, sensor positions, beam diameters and temperatures. The obtained data are matched by simulation results based on the thermo-acoustic model with uncertainties at a level of 10%. The results imply that the primary mechanism for sound generation by the energy deposition of particles propagating in water is the local heating of the medium. The heating results in a fast expansion or contraction and a pressure pulse of bipolar shape is emitted into the surrounding medium. An interesting, widely discussed application of this effect could be the detection of ultra-high energetic cosmic neutrinos in future large-scale acoustic neutrino detectors. For this application a validation of the sound generation mechanism to high accuracy, as achieved with the experiments discussed in this article, is of high importance.     

Detecting filaments in the ultra-high energy cosmic ray distribution

We propose and test new statistical tools to study the distribution of cosmic rays based on the use of the minimal spanning tree. The method described is particularly sensitive to filamentary structures, as those expected to arise from strong sources of charged cosmic rays which get deflected by intervening magnetic fields. We also test the method with data available from the AGASA and SUGAR surface detector arrays.     

Galactic PeV neutrinos

The IceCube experiment has detected two neutrinos with energies between 1 and 10 PeV. They might have originated from Galactic or extragalactic sources of cosmic rays. In the present work we consider hadronic interactions of the diffuse very high energy cosmic rays with the interstellar matter within our Galaxy to explain the PeV neutrino events detected in IceCube. We also expect PeV gamma ray events along with the PeV neutrino events if the observed PeV neutrinos were produced within our Galaxy in hadronic interactions. PeV gamma rays are unlikely to reach us from sources outside our Galaxy due to pair production with cosmic background radiation fields. We suggest that in future with simultaneous detections of PeV gamma rays and neutrinos it would be possible to distinguish between Galactic and extragalactic origins of very high energy neutrinos.     

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.     

Probable spectral steepening of UHEγ-rays of nuclear origin and its consequences

Recent accelerator data based parameterization of the inclusive cross section (c −s) forπ0 production in hadronic collisions and an explicit incorporation of the finiteness of the relevant projectile hadron spectrum suggest a significant steepening in the spectrum (by as much as 0.4 in the spectral index) of the secondaryγ-ray towards the end of the spectrum. We emphasize here that this spectral steepening in conjunction with the possibility that in the bright X-ray binaries the maximum energy to which theγ-ray producing progenitor protons may be accelerated is only ∼ l0 PeV, may imply an effective efficiency forγ-ray production,ε, as reckoned by the PeV arrays, one or two orders smaller than the previous estimates. To explain the genesis of a given PeV photon flux from an X-ray binary, one, therefore, has to. accordingly consider a much higher value of the progenitor proton beam luminosity,L _p . This requirement may raise further questions regarding the actual genesis of PeVγ-rays in X-ray binaries, or alternatively, on the veracity of the high values of the PeV photon fluxes reported by earlier experiments.     

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.     

The sensitivity of the next generation of lunar Cherenkov observations to UHE neutrinos and cosmic rays

We present simulation results for the detection of ultra-high energy (UHE) cosmic ray (CR) and neutrino interactions in the Moon by radio-telescopes. We simulate the expected radio signal at Earth from such interactions, expanding on previous work to include interactions in the sub-regolith layer for single dish and multiple telescope systems. For previous experiments at Parkes, Goldstone (GLUE), and Kalyazin we recalculate the sensitivity to an isotropic flux of UHE neutrinos. We find the published sensitivity for the GLUE experiment to be too high (too optimistic) by an order of magnitude, and consequently the GLUE limit to be too low by an order of magnitude. Our predicted sensitivity for future experiments using the Australia Telescope Compact Array (ATCA) and the Australian SKA Pathfinder (ASKAP) indicate these instruments will be able to detect the more optimistic UHE neutrino flux predictions, while the square kilometre array (SKA) will also be sensitive to all bar one prediction of a diffuse ‘cosmogenic’, or ‘GZK’, neutrino flux.Outstanding theoretical uncertainties at both high-frequency and low-frequency limits currently prevent a reliable estimate of the sensitivity of the lunar Cherenkov technique for UHE cosmic ray (CR) astronomy. Here, we place limits on the effects of large-scale surface roughness on UHE CR detection, and find that when near-surface ‘formation-zone’ effects are ignored, the proposed SKA low-frequency aperture array could detect CR events above 56 EeV at a rate between 15 and 40 times that of the current Pierre Auger Observatory. Should further work indicate that formation-zone effects have little impact on UHE CR sensitivity, observations of the Moon with the SKA would allow directional analysis of UHE cosmic rays, and investigation of correlations with putative cosmic ray source populations, to be conducted with very high statistics.     

Are gamma-ray bursts the sources of ultra-high energy cosmic rays?

We reconsider the possibility that gamma-ray bursts (GRBs) are the sources of the ultra-high energy cosmic rays (UHECRs) within the internal shock model, assuming a pure proton composition of the UHECRs. For the first time, we combine the information from gamma-rays, cosmic rays, prompt neutrinos, and cosmogenic neutrinos quantitatively in a joint cosmic ray production and propagation model, and we show that the information on the cosmic energy budget can be obtained as a consequence. In addition to the neutron model, we consider alternative scenarios for the cosmic ray escape from the GRBs, i.e., that cosmic rays can leak from the sources. We find that the dip model, which describes the ankle in UHECR observations by the pair production dip, is strongly disfavored in combination with the internal shock model because (a) unrealistically high baryonic loadings (energy in protons versus energy in electrons/gamma-rays) are needed for the individual GRBs and (b) the prompt neutrino flux easily overshoots the corresponding neutrino bound. On the other hand, GRBs may account for the UHECRs in the ankle transition model if cosmic rays leak out from the source at the highest energies. In that case, we demonstrate that future neutrino observations can efficiently test most of the parameter space – unless the baryonic loading is much larger than previously anticipated.     

Radio pulses from ultra-high energy atmospheric showers as the superposition of Askaryan and geomagnetic mechanisms

Radio emission in atmospheric showers is currently interpreted in terms of radiation due to the deviation of the charged particles in the magnetic field of the Earth and to the charge excess (Askaryan effect). Each of these mechanisms has a distinctive polarization. The complex signal patterns can be qualitatively explained as the interference (superposition) of the fields induced by each mechanism. In this work we explicitly and quantitatively test a simple phenomenological model based on this idea. The model is constructed by isolating each of the two components at the simulation level and by making use of approximate symmetries for each of the contributions separately. The results of the model are then checked against full ZHAireS Monte Carlo simulations of the electric field calculated from first principles. We show that the simple model describes radio emission at a few percent level in a wide range of shower-observer geometries and on a shower-by-shower basis. As a consequence, this approach provides a simple method to reduce the computing time needed to accurately predict the electric field of radio pulses emitted from air showers, with many practical applications in experimental situations of interest.     

The Antarctic Impulsive Transient Antenna ultra-high energy neutrino detector: Design, performance, and sensitivity for the 2006–2007 balloon flight

We present a comprehensive report on the experimental details of the Antarctic Impulsive Transient Antenna (ANITA) long-duration balloon payload, including the design philosophy and realization, physics simulations, performance of the instrument during its first Antarctic flight completed in January of 2007, and expectations for the limiting neutrino detection sensitivity.     

Prediction of the diffuse neutrino flux from cosmic ray interactions near supernova remnants

In this paper, we present high-energy neutrino spectra from 21 Galactic supernova remnants (SNRs), derived from gamma-ray measurements in the GeV–TeV range. We find that only the strongest sources, i.e. G40.5-0.5 in the north and Vela Junior in the south could be detected as single point sources by IceCube or KM3NeT, respectively. For the first time, it is also possible to derive a diffuse signal by applying the observed correlation between gamma-ray emission and radio signal. Radio data from 234 supernova remnants listed in Green’s catalog are used to show that the total diffuse neutrino flux is approximately a factor of 2.5 higher compared to the sources that are resolved so far. We show that the signal at above 10 TeV energies can actually become comparable to the diffuse neutrino flux component from interactions in the interstellar medium. Recently, the IceCube collaboration announced the detection of a first diffuse signal of astrophysical high-energy neutrinos. Directional information cannot unambiguously reveal the nature of the sources at this point due to low statistics. A number of events come from close to the Galactic center and one of the main questions is whether at least a part of the signal can be of Galactic nature. In this paper, we show that the diffuse flux from well-resolved SNRs is at least a factor of 20 below the observed flux.     

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.     

Thinned simulations of extremely energetic showers in dense media for radio applications

This paper presents our results for full 3-D simulations of very-high to ultra-high energy electromagnetic cascades – and the associated coherent Cherenkov radiation – as might be produced by high-energy neutrino interactions in dense media. Using “thinning” techniques, we develop an algorithm based on the existing “ZHS” code, and demonstrate that the new “ZHS-thinned” code can produce fast and accurate results for showers up to . Using ZHS-thinned, we develop new parameterisations for the radiation from showers in ice, salt, and the lunar regolith, with a separate treatment of the megaregolith (deep regolith). Our parameterisations include for the first time a method to simulate fluctuations in shower length induced by the LPM effect. Our results, which avoid the pit-falls of scaling simulations from lower energies, allow improved calculations of the detection probability for experiments searching for high-energy neutrinos using the radio technique.     

Primary cosmic-ray nucleon spectrum in the energy range 40–400 TeV reconstructed from the secondary gamma-ray spectrum in the stratosphere

We present the spectrum of all primary cosmic-ray nucleons in the energy range 40–400 TeV reconstructed from the spectra of gamma rays recorded along with charged particles at an altitude of 30 km in the stratosphere in the RUNJOB balloon experiment. Gamma rays are produced by the interaction of primary cosmic-ray particles with atomic nuclei of the residual atmosphere and are a component whose spectrum follows the nucleon spectrum at the atmospheric boundary in the approximation of quasi-scaling models. Comparison of the nucleon spectrum with the proton spectrum indicates that the contribution of nucleons from nuclei increases toward the 100-TeV energy region compared to the 1-TeV region. This result confirms the conclusions reached in the JACEE, SOKOL, ATIC-2, and KASCADE experiments: the fraction of helium nuclei and, possibly, heavier nuclei in the primary cosmic-ray flux increases with energy.