Neutrinos of astrophysical origin could be detected through the electromagnetic radiation of the particle showers induced in the atmosphere by their interaction in the Earth. This applies in particular for tau neutrinos of energies $E>10^{16}$eV following Earth-skimming trajectories. The $\sim1^{\circ}$ beaming of the radio emission in the forward direction however implies that the radio signal will likely fly above a detector deployed over a flat site and would therefore not be detected. We study here how a non-flat detector topography can improve the detection probability of these neutrino-induced air showers. We do this by computing with three distinct tools the neutrino detection efficiency for a radio array deployed over a toy-model mountainous terrain, also taking into account experimental and topographic constraints. We show in particular that ground topographies inclined by few degrees only induce detection efficiencies typically three times larger than those obtained for flat areas. We conclude that the topography of the area where the detector is deployed will be a key factor for an experiment like GRAND.

Binary systems composed of a recycled millisecond pulsar and a stellar companion in close orbit could be excellent sites to diagnose pulsar winds. In such systems, the pulsar outflow irradiates and heats up the companion atmosphere, which can lead to the observation of strong day/night modulations in temperature. We demonstrate with particle shower simulations that the particle energy of the wind affects the heating depth in the atmosphere: the wind heat can be deposited above or below the photosphere, leading to different signatures in the observed spectra. We apply our method to four specific systems: We find that systems with cool night side companions showing strong temperature variations can give interesting lower limits on the particle energy in the winds. In particular, if the companion night side of PSR B1957+20 were to be suddenly irradiated, deep heating would only take place if particles with energy > 100 TeV were present. Observational evidence of deep heating in this system thus suggests that i) such particles exist in the pulsar wind and/or ii) binary evolution non-trivially takes the companion to the observed temperature asymmetry. Besides, the observed temperature difference can be maintained only with particle energies of the order of 100 MeV.

Over the last years, radio detection has matured to become a competitive method for the detection of air showers. Arrays of thousands of antennas are now envisioned for the detection of cosmic rays of ultra high energy or neutrinos of astrophysical origin. The data exploitation of such detectors requires to run massive air-shower simulations to evaluate the radio signal at each antenna position. In order to reduce the associated computational cost, we have developed a semi-analytical method for the computation of the emitted radio signal called Radio Morphing. The method consists in computing the radio signal of any air-shower at any location from the simulation of one single reference shower at given positions by i) a scaling of the electric-field amplitude of this reference shower, ii) an isometry on the simulated positions and iii) an interpolation of the radio pulse at the desired position. This technique enables one to compute electric field time traces with characteristics very similar to those obtained with standard computation methods, but with computation times reduced by several orders of magnitude. In this paper, we present this novel tool, explain its methodology, and discuss its limitations. Furthermore, we validate the method on a typical event set for the future GRAND experiment showing that the calculated peak amplitudes are consistent with the results from ZHAireS simulations with a mean offset of +8.5% and a standard deviation of 27.2% in this specific case. This overestimation of the signal strength by Radio Morphing arises mainly from the choice of the underlying reference shower.

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 detection of air-shower events via radio signals requires to develop a trigger algorithm for a clean discrimination between signal and background events in order to reduce the data stream coming from false triggers. In this contribution we will describe an approach to trigger air-shower events on a single-antenna level as well as performing an online reconstruction of the shower parameters using neural networks.

Over the last decades, radio detection of air showers has been established as a promising detection technique for ultrahigh-energy cosmic rays and neutrinos. Very large or dense antenna arrays are necessary to be proficient at collecting information about these particles and understanding their properties accurately. The exploitation of such arrays requires to run massive air-shower simulations to evaluate the radio signal at each antenna position, taking into account features such as the ground topology. In order to reduce computational costs, we have developed a fast computation of the emitted radio signal on the basis of generic shower simulations, called Radio Morphing. The method consists in the calculation of the radio signal of any air-shower by i) a scaling of the electric-field amplitude of a reference air shower to the target shower, ii) an isometry on the simulated positions and iii) an interpolation of the radio pulse at the desired position. This technique enables one to gain many orders of magnitude in CPU time compared to a standard computation. In this contribution, we present this novel tool and explain its methodology. In particular, Radio Morphing will be a key element for the simulation chain of the Giant Radio Array for Neutrino Detection (GRAND) project, that aims at detecting ultra-high-energy neutrinos with an array of 200 000 radio antennas in mountainous regions.os with an array of 200 000 radio antennas in mountainous regions.

Over the last decades, radio detection of air showers has been established as a detection technique for ultra-high-energy cosmic-rays impinging on the Earth's atmosphere with energies far beyond LHC energies. Today's second-generation of digital radio-detection experiments, as e.g. AERA or LOFAR, are becoming competitive in comparison to already standard techniques e.g. fluorescence light detection. Thanks to a detailed understanding of the physics of the radio emission in extensive air showers, simulations of the radio signal are already successfully tested and applied in the reconstruction of cosmic rays. However, the limits of the computational power resources are easily reached when it comes to computing electric fields at the numerous positions requested by large or dense antenna arrays. In the case of mountainous areas as e.g. for the GRAND array, where 3D shower simulations are necessary, the problem arises with even stronger acuity. Therefore we developed a full parametrisation of the emitted radio signal on the basis of generic shower simulations which will reduce the simulation time by orders of magnitudes. In this talk we will present this concept after a short introduction to the concept of the radio detection of air-shower induced by cosmic rays.

The Giant Radio Array for Neutrino Detection (GRAND) aims at detecting ultra-high energy extraterrestrial neutrinos via the extensive air showers induced by the decay of tau leptons created in the interaction of neutrinos under the Earth's surface. Consisting of an array of $\sim10^5$ radio antennas deployed over $\sim 2\times10^5\,\rm {km}^2$, GRAND plans to reach, for the first time, an all-flavor sensitivity of $\sim1.5\times10^{-10} \,\rm GeV\, cm^{-2} \,s^{-1}\, sr^{-1}$ above $5\times10^{17}$ eV and a sub-degree angular resolution, beyond the reach of other planned detectors. We describe here preliminary designs and simulation results, plans for the ongoing, staged approach to the construction of GRAND, and the rich research program made possible by GRAND's design sensitivity and angular resolution.

As of 2023, the low-frequency part of the Square Kilometre Array will go online in Australia. It will constitute the largest and most powerful low-frequency radio-astronomical observatory to date, and will facilitate a rich science programme in astronomy and astrophysics. With modest engineering changes, it will also be able to measure cosmic rays via the radio emission from extensive air showers. The extreme antenna density and the homogeneous coverage provided by more than 60,000 antennas within an area of one km$^2$ will push radio detection of cosmic rays in the energy range around 10$^{17}$ eV to ultimate precision, with superior capabilities in the reconstruction of arrival direction, energy, and an expected depth-of-shower-maximum resolution of 6~g/cm${^2}$.

The SLAC T-510 experiment was designed to compare controlled laboratory measurements of radio emission of particle showers to predictions using particle-level simulations, which are relied upon in ultra-high-energy cosmic-ray air shower detection. Established formalisms for the simulation of radio emission physics, the "endpoint" formalism and the "ZHS" formalism, lead to results which can be explained by a superposition of magnetically induced transverse current radiation and charge-excess radiation due to the Askaryan effect. Here, we present the results of Geant4 simulations for the SLAC T-510 experiment, taking into account the details of the experimental setup (beam energy, target geometry and material, magnetic field configuration, and refraction effects) and their comparison to measured data with respect to e.g. signal polarisation, linearity with magnetic field, and angular distribution. We find that the microscopic calculations reproduce the measurements within uncertainties and describe the data well.

As LOFAR has shown, using a dense array of radio antennas for detecting extensive air showers initiated by cosmic rays in the Earth's atmosphere makes it possible to measure the depth of shower maximum for individual showers with a statistical uncertainty less than $20\,g/cm^2$. This allows detailed studies of the mass composition in the energy region around $10^{17}\,eV$ where the transition from a Galactic to an Extragalactic origin could occur. Since SKA1-low will provide a much denser and very homogeneous antenna array with a large bandwidth of $50-350\,MHz$ it is expected to reach an uncertainty on the $X_{\max}$ reconstruction of less than $10\,g/cm^2$. We present first results of a simulation study with focus on the potential to reconstruct the depth of shower maximum for individual showers to be measured with SKA1-low. In addition, possible influences of various parameters such as the numbers of antennas included in the analysis or the considered frequency bandwidth will be discussed.

The SLAC T-510 experiment was designed to reproduce the physics of radio emission from air showers caused by ultra-high energy cosmic rays in a controlled lab experiment with the goal to test established formalisms for simulation of radio emission physics: the "end-point" formalism and the "ZHS" formalism. Simulation results derived with these formalisms can be explained by a superposition of magnetically induced transverse current radiation and the Askaryan (charge-excess) effect. Here, we present results of Geant4 simulations of the experiment with both formalisms, taking into account the details of the experimental setup (beam energy, target geometry and material, magnetic field configuration, and refraction effects) to test this hypothesis

For fifty years, cosmic-ray air showers have been detected by their radio emission. We present the first laboratory measurements that validate electrodynamics simulations used in air shower modeling. An experiment at SLAC provides a beam test of radio-frequency (RF) radiation from charged particle cascades in the presence of a magnetic field, a model system of a cosmic-ray air shower. This experiment provides a suite of controlled laboratory measurements to compare to particle-level simulations of RF emission, which are relied upon in ultra-high-energy cosmic-ray air shower detection. We compare simulations to data for intensity, linearity with magnetic field, angular distribution, polarization, and spectral content. In particular, we confirm modern predictions that the magnetically induced emission in a dielectric forms a cone that peaks at the Cherenkov angle and show that the simulations reproduce the data within systematic uncertainties.

We have performed a search for muon neutrinos from dark matter annihilation in the center of the Sun with the 79-string configuration of the IceCube neutrino telescope. For the first time, the DeepCore sub-array is included in the analysis, lowering the energy threshold and extending the search to the austral summer. The 317 days of data collected between June 2010 and May 2011 are consistent with the expected background from atmospheric muons and neutrinos. Upper limits are therefore set on the dark matter annihilation rate, with conversions to limits on spin-dependent and spin-independent WIMP-proton cross-sections for WIMP masses in the range 20 - 5000 GeV. These are the most stringent spin-dependent WIMP-proton cross-sections limits to date above 35 GeV.

As of 2023, the Square Kilometre Array will constitute the world's largest radio telescope, offering unprecedented capabilities for a diverse science programme in radio astronomy. At the same time, the SKA will be ideally suited to detect extensive air showers initiated by cosmic rays in the Earth's atmosphere via their radio emission. With its very dense and uniform antenna spacing in a fiducial area of one km$^2$ and its large bandwidth of 50-350 MHz, the low-frequency part of the SKA will provide very precise measurements of individual cosmic ray air showers. These precision measurements will allow detailed studies of the mass composition of cosmic rays in the energy region of transition from a Galactic to an extragalactic origin. Also, the SKA will facilitate three-dimensional "tomography" of the electromagnetic cascades of air showers, allowing the study of particle interactions at energies beyond the reach of the LHC. Finally, studies of possible connections between air showers and lightning initiation can be taken to a new level with the SKA. We discuss the science potential of air shower detection with the SKA and report on the technical requirements and project status.

We report on the observation of two neutrino-induced events which have an estimated deposited energy in the IceCube detector of 1.04 $\pm$ 0.16 and 1.14 $\pm$ 0.17 PeV, respectively, the highest neutrino energies observed so far. These events are consistent with fully contained particle showers induced by neutral-current $\nu_{e,\mu,\tau}$ ($\bar\nu_{e,\mu,\tau}$) or charged-current $\nu_{e}$ ($\bar\nu_{e}$) interactions within the IceCube detector. The events were discovered in a search for ultra-high energy neutrinos using data corresponding to 615.9 days effective livetime. The expected number of atmospheric background is $0.082 \pm 0.004 \text{(stat)}^{+0.041}_{-0.057} \text{(syst)}$. The probability to observe two or more candidate events under the atmospheric background-only hypothesis is $2.9\times10^{-3}$ ($2.8\sigma$) taking into account the uncertainty on the expected number of background events. These two events could be a first indication of an astrophysical neutrino flux, the moderate significance, however, does not permit a definitive conclusion at this time.

A search for muon neutrinos from dark matter annihilations in the Galactic Center region has been performed with the 40-string configuration of the IceCube Neutrino Observatory using data collected in 367 days of live-time starting in April 2008. The observed fluxes were consistent with the atmospheric background expectations. Upper limits on the self-annihilation cross-section are obtained for dark matter particle masses ranging from 100 GeV to 10 TeV. In the case of decaying dark matter, lower limits on the lifetime have been determined for masses between 200 GeV and 20 TeV.

We report the first measurement of the atmospheric electron neutrino flux in the energy range between approximately 80 GeV and 6 TeV, using data recorded during the first year of operation of IceCube's DeepCore low energy extension. Techniques to identify neutrinos interacting within the DeepCore volume and veto muons originating outside the detector are demonstrated. A sample of 1029 events is observed in 281 days of data, of which 496 $\pm$ 66(stat.) $\pm$ 88(syst.) are estimated to be cascade events, including both electron neutrino and neutral current events. The rest of the sample includes residual backgrounds due to atmospheric muons and charged current interactions of atmospheric muon neutrinos. The flux of the atmospheric electron neutrinos is consistent with models of atmospheric neutrinos in this energy range. This constitutes the first observation of electron neutrinos and neutral current interactions in a very large volume neutrino telescope optimized for the TeV energy range.

The IceCube Neutrino Observatory, approximately 1 km^3 in size, is now complete with 86 strings deployed in the Antarctic ice. IceCube detects the Cherenkov radiation emitted by charged particles passing through or created in the ice. To realize the full potential of the detector, the properties of light propagation in the ice in and around the detector must be well understood. This report presents a new method of fitting the model of light propagation in the ice to a data set of in-situ light source events collected with IceCube. The resulting set of derived parameters, namely the measured values of scattering and absorption coefficients vs. depth, is presented and a comparison of IceCube data with simulations based on the new model is shown.

We report on searches for neutrino sources at energies above 200 GeV in the Northern sky of the galactic plane, using the data collected by the South Pole neutrino telescopes IceCube and AMANDA. The galactic region considered here includes the Local Arm towards the Cygnus region and our closest approach to the Perseus Arm. The data have been collected between 2007 and 2009 when AMANDA was an integrated part of IceCube, which was still under construction and operated with 22-strings (2007-8) and 40-strings (2008-9) of optical modules deployed in the ice. By combining the larger IceCube detector with the lower energy threshold of the more compact AMANDA detector, we obtain an improved sensitivity at energies below $\sim$10 TeV with respect to previous searches. The analyses presented here are: a scan for point sources within the galactic plane; a search optimized for multiple and extended sources in the Cygnus region, which might be below the sensitivity of the point source scan; and studies of seven pre-selected neutrino source candidates. For one of them, Cygnus X-3, a time-dependent search for neutrinos in coincidence with observed radio and X-ray flares has been performed. No evidence of a signal is found, and upper limits are reported for each of the searches. We investigate neutrino spectra proportional to E$^{-2}$ and E$^{-3}$ to cover the entire range of possible spectra. The soft E$^{-3}$ spectrum results in an energy distribution similar to a source with cut-off below $\sim$50 TeV. For the considered region of the galactic plane, the 90% confidence level muon neutrino flux upper limits are in the range E$^3$dN/dE$\sim 5.4 - 19.5 \times 10^{-11} \rm{TeV^{2} cm^{-2} s^{-1}}$ for point-like neutrino sources in the energy region [180.0 GeV - 20.5 TeV]. These represent the most stringent upper limits for soft-spectra neutrino sources within the Galaxy reported to date.

We report on the observation of anisotropy in the arrival direction distribution of cosmic rays at PeV energies. The analysis is based on data taken between 2009 and 2012 with the IceTop air shower array at the South Pole. IceTop, an integral part of the IceCube detector, is sensitive to cosmic rays between 100 TeV and 1 EeV. With the current size of the IceTop data set, searches for anisotropy at the 10^-3 level can, for the first time, be extended to PeV energies. We divide the data set into two parts with median energies of 400 TeV and 2 PeV, respectively. In the low energy band, we observe a strong deficit with an angular size of about 30 degrees and an amplitude of (-1.58 +/- 0.46 (stat) +/- 0.52 (sys)) x 10^(-3) at a location consistent with previous observations of cosmic rays with the IceCube neutrino detector. The study of the high energy band shows that the anisotropy persists to PeV energies and increases in amplitude to (-3.11 +/- 0.38 (stat) +/- 0.96 (sys)) x 10^(-3).

In cosmic ray air showers, the muon lateral separation from the center of the shower is a measure of the transverse momentum that the muon parent acquired in the cosmic ray interaction. IceCube has observed cosmic ray interactions that produce muons laterally separated by up to 400 m from the shower core, a factor of 6 larger distance than previous measurements. These muons originate in high pT (> 2 GeV/c) interactions from the incident cosmic ray, or high-energy secondary interactions. The separation distribution shows a transition to a power law at large values, indicating the presence of a hard pT component that can be described by perturbative quantum chromodynamics. However, the rates and the zenith angle distributions of these events are not well reproduced with the cosmic ray models tested here, even those that include charm interactions. This discrepancy may be explained by a larger fraction of kaons and charmed particles than is currently incorporated in the simulations.

We present the first results in the search for relativistic magnetic monopoles with the IceCube detector, a subsurface neutrino telescope located in the South Polar ice cap containing a volume of 1 km$^{3}$. This analysis searches data taken on the partially completed detector during 2007 when roughly 0.2 km$^{3}$ of ice was instrumented. The lack of candidate events leads to an upper limit on the flux of relativistic magnetic monopoles of $\Phi_{\mathrm{90%C.L.}}\sim 3\e{-18}\fluxunits$ for $\beta\geq0.8$. This is a factor of 4 improvement over the previous best experimental flux limits up to a Lorentz boost $\gamma$ below $10^{7}$. This result is then interpreted for a wide range of mass and kinetic energy values.

The measurement of muon energy is critical for many analyses in large Cherenkov detectors, particularly those that involve separating extraterrestrial neutrinos from the atmospheric neutrino background. Muon energy has traditionally been determined by measuring the specific energy loss (dE/dx) along the muon's path and relating the dE/dx to the muon energy. Because high-energy muons (E_mu > 1 TeV) lose energy randomly, the spread in dE/dx values is quite large, leading to a typical energy resolution of 0.29 in log10(E_mu) for a muon observed over a 1 km path length in the IceCube detector. In this paper, we present an improved method that uses a truncated mean and other techniques to determine the muon energy. The muon track is divided into separate segments with individual dE/dx values. The elimination of segments with the highest dE/dx results in an overall dE/dx that is more closely correlated to the muon energy. This method results in an energy resolution of 0.22 in log10(E_mu), which gives a 26% improvement. This technique is applicable to any large water or ice detector and potentially to large scintillator or liquid argon detectors.

The mass composition of high energy cosmic rays depends on their production, acceleration, and propagation. The study of cosmic ray composition can therefore reveal hints of the origin of these particles. At the South Pole, the IceCube Neutrino Observatory is capable of measuring two components of cosmic ray air showers in coincidence: the electromagnetic component at high altitude (2835 m) using the IceTop surface array, and the muonic component above ~1 TeV using the IceCube array. This unique detector arrangement provides an opportunity for precision measurements of the cosmic ray energy spectrum and composition in the region of the knee and beyond. We present the results of a neural network analysis technique to study the cosmic ray composition and the energy spectrum from 1 PeV to 30 PeV using data recorded using the 40-string/40-station configuration of the IceCube Neutrino Observatory.

Gamma-ray induced air showers are notable for their lack of muons, compared to hadronic showers. Hence, air shower arrays with large underground muon detectors can select a sample greatly enriched in photon showers by rejecting showers containing muons. IceCube is sensitive to muons with energies above ~500 GeV at the surface, which provides an efficient veto system for hadronic air showers with energies above 1 PeV. One year of data from the 40-string IceCube configuration was used to perform a search for point sources and a Galactic diffuse signal. No sources were found, resulting in a 90% C.L. upper limit on the ratio of gamma rays to cosmic rays of 1.2 x 10^(-3)for the flux coming from the Galactic Plane region (-80 deg < l < -30 deg; -10 deg < b < 5 deg) in the energy range 1.2 - 6.0 PeV. In the same energy range, point source fluxes with E^(-2) spectra have been excluded at a level of (E/TeV)^2 d\Phi/dE ~ 10^(-12)-10^(-11) cm^2/s/TeV depending on source declination. The complete IceCube detector will have a better sensitivity, due to the larger detector size, improved reconstruction and vetoing techniques. Preliminary data from the nearly-final IceCube detector configuration has been used to estimate the 5 year sensitivity of the full detector. It is found to be more than an order of magnitude better, allowing the search for PeV extensions of known TeV gamma-ray emitters.

We present a fast likelihood method for including event-level neutrino telescope data in parameter explorations of theories for new physics, and announce its public release as part of DarkSUSY 5.0.6. Our construction includes both angular and spectral information about neutrino events, as well as their total number. We also present a corresponding measure for simple model exclusion, which can be used for single models without reference to the rest of a parameter space. We perform a number of supersymmetric parameter scans with IceCube data to illustrate the utility of the method: example global fits and a signal recovery in the constrained minimal supersymmetric standard model (CMSSM), and a model exclusion exercise in a 7-parameter phenomenological version of the MSSM. The final IceCube detector configuration will probe almost the entire focus-point region of the CMSSM, as well as a number of MSSM-7 models that will not otherwise be accessible to e.g. direct detection. Our method accurately recovers the mock signal, and provides tight constraints on model parameters and derived quantities. We show that the inclusion of spectral information significantly improves the accuracy of the recovery, providing motivation for its use in future IceCube analyses.

IceTop, the surface component of the IceCube Neutrino Observatory at the South Pole, is an air shower array with an area of 1 km2. The detector allows a detailed exploration of the mass composition of primary cosmic rays in the energy range from about 100 TeV to 1 EeV by exploiting the correlation between the shower energy measured in IceTop and the energy deposited by muons in the deep ice. In this paper we report on the technical design, construction and installation, the trigger and data acquisition systems as well as the software framework for calibration, reconstruction and simulation. Finally the first experience from commissioning and operating the detector and the performance as an air shower detector will be discussed.