Joint searches by FACT, H.E.S.S., MAGIC and VERITAS for VHE gamma-ray emission associated with neutrinos detected by IceCube

The sources of the astrophysical flux of high-energy neutrinos detected by IceCube are still largely unknown, but searches for temporal and spatial correlation between neutrinos and electromagnetic radiation are a promising approach in this endeavor. All major imaging atmospheric Cherenkov telescopes (IACTs) - FACT, H.E.S.S., MAGIC, and VERITAS - operate an active follow-up program of target-of-opportunity observations of neutrino alerts issued by IceCube. These programs use several complementary neutrino alert streams. A publicly distributed alert stream is formed by individual high-energy neutrino candidate events of potentially astrophysical origin, such as IceCube-170922A (which could be linked to the flaring blazar TXS 0506+056). A privately distributed alert stream is formed by clusters of neutrino events in time and space around either pre-selected gamma-ray sources or anywhere in the sky. Here, we present joint searches for multi-wavelength emission associated with a set of IceCube alerts, both private and public, received through mid-January 2021. We will give an overview of the programs of the participating IACTs. We will showcase the various follow-up and data analysis strategies employed in response to the different alert types and various possible counterpart scenarios. Finally, we will present results from a combined analysis of the VHE gamma-ray observations obtained across all involved instruments, as well as relevant multi-wavelength data

Searching for VHE gamma-ray emission associated with IceCube neutrino alerts using FACT, H.E.S.S., MAGIC, and VERITAS

The realtime follow-up of neutrino events is a promising approach to searchfor astrophysical neutrino sources. It has so far provided compelling evidencefor a neutrino point source: the flaring gamma-ray blazar TXS 0506+056 observedin coincidence with the high-energy neutrino IceCube-170922A detected byIceCube. The detection of very-high-energy gamma rays (VHE, $\mathrm{E} >100\,\mathrm{GeV}$) from this source helped establish the coincidence andconstrained the modeling of the blazar emission at the time of the IceCubeevent. The four major imaging atmospheric Cherenkov telescope arrays (IACTs) -FACT, H.E.S.S., MAGIC, and VERITAS - operate an active follow-up program oftarget-of-opportunity observations of neutrino alerts sent by IceCube. Thisprogram has two main components. One are the observations of known gamma-raysources around which a cluster of candidate neutrino events has been identifiedby IceCube (Gamma-ray Follow-Up, GFU). Second one is the follow-up of singlehigh-energy neutrino candidate events of potential astrophysical origin such asIceCube-170922A. GFU has been recently upgraded by IceCube in collaborationwith the IACT groups. We present here recent results from the IACT follow-upprograms of IceCube neutrino alerts and a description of the upgraded IceCubeGFU system.<br

A model-independent analysis of neutrino flares detected in IceCube from X-ray selected blazars

Blazars are among the most powerful steady sources in the Universe. Multi-messenger searches for blazars have traditionally focused on their gamma-ray emission, which can be produced simultaneously with neutrinos in photohadronic interactions. However, X-ray data can be equally vital to constrain the SED of these sources, since the hadronically co-produced gamma-rays could get absorbed by the ambient photon fields and cascade down to X-ray energies before escaping. In this work, we present the outline for an untriggered, time-dependent analysis of neutrino flares from the direction of X-ray selected blazars using 10 years of IceCube data. A binomial test will be performed on the population to reveal if a subcategory of sources has statistically significant emission. The sources are selected from RomaBZCat, and the p-values and best-fit flare parameters are obtained for each source using the method of unbinned likelihood maximisation

A next-generation optical sensor for IceCube-Gen2

For the in-ice component of the next generation neutrino observatory at the South Pole, IceCubeGen2, a new sensor module is being developed, which is an evolution of the D-Egg and mDOM sensors developed for the IceCube Upgrade. The sensor design features up to 18 4-inch PMTs distributed homogeneously in a borosilicate glass pressure vessel. Challenges arise for the mechanical design from the tight constraints on the bore hole diameter (which will be 2 inches smaller than for IceCube Upgrade) and from the close packing of the PMTs. The electronics design must meet the space constraints posed by the mechanical design as well as the power consumption and cost considerations from over 10,000 optical modules being deployed. This contribution presents forward-looking solutions to these design considerations. Prototype modules will be installed and integrated in the IceCube Upgrade

A Combined Fit of the Diffuse Neutrino Spectrum using IceCube Muon Tracks and Cascades

The IceCube Neutrino Observatory first observed a diffuse flux of high energy astrophysical neutrinos in 2013. Since then, this observation has been confirmed in multiple detection channels such as high energy starting events, cascades, and through-going muon tracks. Combining these event selections into a high statistics global fit of 10 years of IceCube’s neutrino data could strongly improve the understanding of the diffuse astrophysical neutrino flux: challenging or confirming the simple unbroken power-law flux model as well as the astrophysical neutrino flux composition. One key component of such a combined analysis is the consistent modelling of systematic uncertainties of different event selections. This can be achieved using the novel SnowStorm Monte Carlo method which allows constraints to be placed on multiple systematic parameters from a single simulation set. We will report on the status of a new combined analysis of through-going muon tracks and cascades. It is based on a consistent all flavor neutrino signal and background simulation using, for the first time, the SnowStorm method to analyze IceCube’s high-energy neutrino data. Estimated sensitivities for the energy spectrum of the diffuse astrophysical neutrino flux will be shown

A novel microstructure-based model to explain the IceCube ice anisotropy

The IceCube Neutrino Observatory instruments about 1 km3 of deep, glacial ice at the geographic South Pole using 5160 photomultipliers to detect Cherenkov light of charged relativistic particles. Most of IceCube’s science goals rely heavily on an ever more precise understanding of the optical properties of the instrumented ice. A curious light propagation effect observed by the experiment is an anisotropic attenuation, which is aligned with the local flow of the ice. Having recently identified curved photon trajectories resulting from asymmetric light diffusion in the birefringent polycrystalline microstructure of the ice as the most likely underlying cause of this effect, work is now ongoing to optimize the model parameters (effectively deducing the average crystal size and shape in the detector). We present the parametrization of the birefringence effect in our photon propagation simulation, the fitting procedures and results as well as the impact of the new ice model on data-MC agreement

A Posterior Analysis on IceCube Double Pulse Tau Neutrino Candidates

The IceCube Neutrino Observatory at the South Pole detects Cherenkov light emitted by charged secondary particles created by primary neutrino interactions. Double pulse waveforms can arise from charged current interactions of astrophysical tau neutrinos with nucleons in the ice and the subsequent decay of tau leptons. The previous 8-year tau double pulse analysis found three tau neutrino candidate events. Among them, the most promising one observed in 2014 is located very near the dust layer in the middle of the detector. A posterior analysis on this event will be presented in this paper, using a new ice model treatment with continuously varying nuisance parameters to do the targeted Monte Carlo re-simulation for tau and other background neutrino ensembles. The impact of different ice models on the expected signal and background statistics will also be discussed

Update on the Combined Analysis of Muon Measurements from Nine Air Shower Experiments

Over the last two decades, various experiments have measured muon densities in extensive air showers over several orders of magnitude in primary energy. While some experiments observed differences in the muon densities between simulated and experimentally measured air showers, others reported no discrepancies. We will present an update of the meta-analysis of muon measurements from nine air shower experiments, covering shower energies between a few PeV and tens of EeV and muon threshold energies from a few 100 MeV to about 10GeV. In order to compare measurements from different experiments, their energy scale was cross-calibrated and the experimental data has been compared using a universal reference scale based on air shower simulations. Above 10 PeV, we find a muon excess with respect to simulations for all hadronic interaction models, which is increasing with shower energy. For EPOS-LHC and QGSJet-II.04 the significance of the slope of the increase is analyzed in detail under different assumptions of the individual experimental uncertainties