The Giant Radio Array for Neutrino Detection (GRAND) is a planned large-scale observatory of ultra-high-energy (UHE) cosmic particles, with energies exceeding 10^8 GeV. Its goal is to solve the long-standing mystery of the origin of UHE cosmic rays. To do this, GRAND will detect an unprecedented number of UHE cosmic rays and search for the undiscovered UHE neutrinos and gamma rays associated to them with unmatched sensitivity. GRAND will use large arrays of antennas to detect the radio emission coming from extensive air showers initiated by UHE particles in the atmosphere. Its design is modular: 20 separate, independent sub-arrays, each of 10 000 radio antennas deployed over 10 000 km^2. A staged construction plan will validate key detection techniques while achieving important science goals early. Here we present the science goals, detection strategy, preliminary design, performance goals, and construction plans for GRAND.

The aim of the Air Microwave Yield (AMY) experiment is to investigate the Molecular Bremsstrahlung Radiation (MBR) emitted from an electron beam induced air-shower. The measurements have been performed with a 510 MeV electron beam at the Beam Test Facility (BTF) of Frascati INFN National Laboratories in a wide frequency range between 1 and 20 GHz. We present the experimental apparatus and the first results of the measurements. Contrary to what have been reported in a previous similar experiment~\cite{Gorham-SLAC}, we have found that the intensity of the emission is strongly influenced by the particular time structure of the accelerator beam. This makes very difficult the interpretation of the emission process and a realistic extrapolation of the emission yield to the plasma generated during the development of an atmospheric shower.

The lunar Askaryan technique is a method to study the highest-energy cosmic rays, and their predicted counterparts, the ultra-high-energy neutrinos. By observing the Moon with a radio telescope, and searching for the characteristic nanosecond-scale Askaryan pulses emitted when a high-energy particle interacts in the outer layers of the Moon, the visible lunar surface can be used as a detection area. Several previous experiments, at Parkes, Goldstone, Kalyazin, Westerbork, the ATCA, Lovell, LOFAR, and the VLA, have developed the necessary techniques to search for these pulses, but existing instruments have lacked the necessary sensitivity to detect the known flux of cosmic rays from such a distance. This will change with the advent of the SKA. The Square Kilometre Array (SKA) will be the world's most powerful radio telescope. To be built in southern Africa, Australia and New Zealand during the next decade, it will have an unsurpassed sensitivity over the key 100 MHz to few-GHZ band. We introduce a planned experiment to use the SKA to observe the highest-energy cosmic rays and, potentially, neutrinos. The estimated event rate will be presented, along with the predicted energy and directional resolution. Prospects for directional studies with phase 1 of the SKA will be discussed, as will the major technical challenges to be overcome to make full use of this powerful instrument. Finally, we show how phase 2 of the SKA could provide a vast increase in the number of detected cosmic rays at the highest energies, and thus to provide new insight into their spectrum and origin.

The lunar Askaryan technique, which involves searching for Askaryan radio pulses from particle cascades in the outer layers of the Moon, is a method for using the lunar surface as an extremely large detector of ultra-high-energy particles. The high time resolution required to detect these pulses, which have a duration of around a nanosecond, puts this technique in a regime quite different from other forms of radio astronomy, with a unique set of associated technical challenges which have been addressed in a series of experiments by various groups. Implementing the methods and techniques developed by these groups for detecting lunar Askaryan pulses will be important for a future experiment with the Square Kilometre Array (SKA), which is expected to have sufficient sensitivity to allow the first positive detection using this technique. Key issues include correction for ionospheric dispersion, beamforming, efficient triggering, and the exclusion of spurious events from radio-frequency interference. We review the progress in each of these areas, and consider the further progress expected for future application with the SKA.

The Pierre Auger Observatory is the world's largest cosmic ray observatory. Our current exposure reaches nearly 40,000 km$^2$ str and provides us with an unprecedented quality data set. The performance and stability of the detectors and their enhancements are described. Data analyses have led to a number of major breakthroughs. Among these we discuss the energy spectrum and the searches for large-scale anisotropies. We present analyses of our X$_{max}$ data and show how it can be interpreted in terms of mass composition. We also describe some new analyses that extract mass sensitive parameters from the 100% duty cycle SD data. A coherent interpretation of all these recent results opens new directions. The consequences regarding the cosmic ray composition and the properties of UHECR sources are briefly discussed.

Joint contributions of the Pierre Auger and Telescope Array Collaborations to the 33rd International Cosmic Ray Conference, Rio de Janeiro, Brazil, July 2013: cross-calibration of the fluorescence telescopes, large scale anisotropies and mass composition.

The current status of searches for ultra-high energy neutrinos and photons using air showers is reviewed. Regarding both physics and observational aspects, possible future research directions are indicated.

The Cherenkov radio pulse emitted by hadronic showers in ice is calculated for showers of energies in the EeV range. This is obtained with 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 and optimized 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 pions and muons, is found to be below 1%. 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.

Recent measurements suggest that extensive air showers initiated by ultra-high energy cosmic rays (UHECR) emit signals in the microwave band of the electromagnetic spectrum caused by the collisions of the free-electrons with the atmospheric neutral molecules in the plasma produced by the passage of the shower. Such emission is isotropic and could allow the detection of air showers with 100% duty cycle and a calorimetric-like energy measurement, a significant improvement over current detection techniques. We have built MIDAS (MIcrowave Detection of Air Showers), a prototype of microwave detector, which consists of a 4.5 m diameter antenna with a cluster of 53 feed-horns in the 4 GHz range. The details of the prototype and first results will be presented.

By means of Monte Carlo simulations of extensive air showers (EAS), we have performed a comprehensive study of the shower to shower fluctuations affecting the longitudinal and lateral development of EAS. We split the fluctuations into physical fluctuations and those induced by the thinning procedure customarily applied to simulate showers at EeV energies and above. We study the influence of thinning on the calculation of the shower to shower fluctuations in the simulations. For thinning levels larger than 10^(-5) - 10^(-6), the determination of the shower to shower fluctuations is hampered by the artificial fluctuations induced by the thinning procedure. However, we show that shower to shower fluctuations can still be approximately estimated, and we provide expressions to calculate them. The influence of fluctuations of the depth of first interaction on the determination of shower to shower fluctuations is also addressed.

The electric field of the Cherenkov radio pulse produced by a single charged particle track in a dielectric medium is derived from first principles. An algorithm is developed to obtain the pulse in the time domain for numerical calculations. The algorithm is implemented in a Monte Carlo simulation of electromagnetic showers in dense media (specifically designed for coherent radio emission applications) as might be induced by interactions of ultra-high energy neutrinos. The coherent Cherenkov radio emission produced by such showers is obtained simultaneously both in the time and frequency domains. A consistency check performed by Fourier-transforming the pulse in time and comparing it to the frequency spectrum obtained directly in the simulations yields, as expected, fully consistent results. The reversal of the time structure inside the Cherenkov cone and the signs of the corresponding pulses are addressed in detail. The results, besides testing algorithms used for reference calculations in the frequency domain, shed new light into the properties of the radio pulse in the time domain. The shape of the pulse in the time domain is directly related to the depth development of the excess charge in the shower and its width to the observation angle with respect to the Cherenkov direction. This information can be of great practical importance for interpreting actual data.

Inclined air showers - those arriving at ground with zenith angle with respect to the vertical theta > 60 deg - are characterised by the dominance of the muonic component at ground which is accompanied by an electromagnetic halo produced mainly by muon decay and muon interactions. By means of Monte Carlo simulations we give a full characterisation of the particle densities at ground in ultra-high energy inclined showers as a function of primary energy and mass composition, as well as for different hadronic models assumed in the simulations. We also investigate the effect of intrinsic shower-to-shower fluctuations in the particle densities.

LUNASKA (Lunar UHE Neutrino Astrophysics with the Square Kilometre Array) is a theoretical and experimental project developing the lunar Cherenkov technique for the next generation of giant radio-telescope arrays. Here we report on a series of observations with ATCA (the Australia Telescope Compact Array). Our current observations use three of the six 22m ATCA antennas with a 600 MHz bandwidth at 1.2-1.8 GHz, analogue dedispersion filters to correct for the typical night-time ionospheric dispersion, and state-of-the-art 2 GHz FPGA-based digital pulse detection hardware. We have observed so as to maximise the UHE neutrino sensitivity in the region surrounding the galactic centre and to Centaurus A, to which current limits on the highest-energy neutrinos are relatively weak.

We calculate the neutrino emission from individual gamma-ray bursts observed by the BATSE detector on the Compton Gamma-Ray Observatory. Neutrinos are produced by photoproduction of pions when protons interact with photons in the region where the kinetic energy of the relativistic fireball is dissipated allowing the acceleration of electrons and protons. We also consider models where neutrinos are predominantly produced on the radiation surrounding the newly formed black hole. From the observed redshift and photon flux of each individual burst, we compute the neutrino flux in a variety of models based on the assumption that equal kinetic energy is dissipated into electrons and protons. Where not measured, the redshift is estimated by other methods. Unlike previous calculations of the universal diffuse neutrino flux produced by all gamma-ray bursts, the individual fluxes (compiled at http://www.arcetri.astro.it/~dafne/grb/) can be directly compared with coincident observations by the AMANDA telescope at the South Pole. Because of its large statistics, our predictions are likely to be representative for future observations with larger neutrino telescopes.

We explore the constant intensity cut method that is widely used for the derivation of the cosmic ray energy spectrum, for comparisons of data obtained at different atmospheric depths, for measuring average shower profiles, and for estimates of the proton-air cross section from extensive air shower data. The constant intensity cut method is based on the selection of air showers by charged particle or muon size and therefore is subject to intrinsic shower fluctuations. We demostrate that, depending on the selection method, shower fluctuations can strongly influence the characteristics of the selected showers. Furthermore, a mixture of different primaries in the cosmic ray flux complicates the interpretation of measurements based on the method of constant intensity cuts. As an example we consider data published by the Akeno Collaboration. The interpretation of the Akeno measurements suggests that more than 60-70% of cosmic ray primaries in the energy range 10^{16}-19^{17} eV are heavy nuclei. Our conclusions depend only weakly on the hadronic interaction model chosen to perform the simulations, namely SIBYLL and QGSjet.

We compare in detail the results of simulations of electromagnetic showers in ice in the GeV-TeV energy range, using both the GEANT package and the ZHS Monte Carlo, a code specifically designed to calculate coherent Cerenkov radio pulses from electromagnetic showers in dense media. The longitudinal and lateral profiles as well as the tracklengths, and excess tracklengths are shown to agree at the 10% level. We briefly comment on the negligible influence of the Landau-Pomeranchuk-Migdal effect on the total shower tracklength. Our results are relevant for experiments exploiting the radio Cerenkov technique in dense media as well as for other detectors that rely on the Cerenkov effect in water or ice.

If spacetime has more than four dimensions, ultra-high energy cosmic rays may create microscopic black holes. Black holes created by cosmic neutrinos in the Earth will evaporate, and the resulting hadronic showers, muons, and taus may be detected in neutrino telescopes below the Earth's surface. We simulate such events in detail and consider black hole cross sections with and without an exponential suppression factor. We find observable rates in both cases: for conservative cosmogenic neutrino fluxes, several black hole events per year are observable at the IceCube detector; for fluxes at the Waxman-Bahcall bound, tens of events per year are possible. We also present zenith angle and energy distributions for all three channels. The ability of neutrino telescopes to differentiate hadrons, muons, and possibly taus, and to measure these distributions provides a unique opportunity to identify black holes, to experimentally constrain the form of black hole production cross sections, and to study Hawking evaporation.

We show that a kilometer-scale neutrino observatory, though optimized for detecting neutrinos of TeV to PeV energy, can reveal the science associated with the enigmatic super-EeV radiation in the Universe. Speculations regarding its origin include heavy relics from the early Universe, particle interactions associated with the Greisen cutoff, and topological defects which are remnant cosmic structures associated with phase transitions in grand unified gauge theories. We show that it is a misconception that new instruments optimized to EeV energy can exclusively do this important science. Because kilometer-scale neutrino telescopes such as IceCube can reject the atmospheric neutrino background by identifying the very high energy of the signal events, they have sensitivity over the full solid angle, including the horizon where most of the signal is concentrated. This is critical because upgoing neutrino-induced muons, considered in previous calculations, are absorbed by the Earth. Previous calculations have underestimated the event rates of IceCube for EeV signals by over one order of magnitude.

We show that a kilometer-scale neutrino observatory, though optimized for TeV to PeV energy, is sensitive to the neutrinos associated with super-EeV sources. These include super-heavy relics, neutrinos associated with the Greisen cutoff, and topological defects which are remnant cosmic structures associated with phase transitions in grand unified gauge theories. It is a misconception that new instruments optimized to EeV energy are required to do this important science. Because kilometer-scale neutrino telescopes can reject atmospheric backgrounds by establishing the very high energy of the signal events, they have sensitivity over the full solid angle, including the horizon where most of the signal is concentrated. This is important because up-going neutrino-induced muons, routinely considered in previous calculations, are absorbed by the Earth.

Following Waxman and Bahcall we calculate the event rate, energy and zenith angle dependence of neutrinos produced in the fireball model of gamma ray bursts (GRB). We emphasize the primary importance of i) burst-to-burst fluctuations and ii) absorption of the neutrinos in the Earth. From the astronomical point of view, we draw attention to the sensitivity of neutrino measurements to the boost Lorentz factor of the fireball $\Gamma$, which is central to the fireball model, and only indirectly determined by follow-up observations. Fluctuations result in single bursts emitting multiple neutrinos, making it possible to determine the flavor composition of a beam observed after a baseline of thousands of Megaparsecs.

Because of the limited size of the satellite-borne instruments, it has not been possible to observe the flux of gamma ray bursts (GRB) beyond GeV energy. We here show that it is possible to detect the GRB radiation of TeV energy and above, by detecting the muon secondaries produced when the gamma rays shower in the Earth's atmosphere. Observation is made possible by the recent commissioning of underground detectors (AMANDA, the Lake Baikal detector and MILAGRO) which combine a low muon threshold of a few hundred GeV or less, with a large effective area of 10^3 m^2 or more. Observations will not only provide new insights in the origin and characteristics of GRB, they also provide quantitative information on the diffuse infrared background.