Status, performance, and first results of the IceTop array

We describe the design and performance of IceTop, the air shower array on top of the IceCube neutrino detector. After the 2008/09 antarctic summer season both detectors are deployed at almost 3/4 of their design size. With the current IceTop 59 stations we can start the study of showers of energy well above 10$^{17}$ eV. The paper also describes the first results from IceTop and our plans to study the cosmic ray composition using several different types of analysis.

💡 Research Summary

IceTop is the surface air‑shower array that sits atop the IceCube neutrino telescope at the South Pole. Designed to consist of 80 stations, each comprising two frozen water tanks instrumented with a pair of digital optical modules (DOMs) – one high‑gain and one low‑gain – the array provides a large dynamic range from a few vertical equivalent muons (VEM) up to several thousand VEM. By the end of the 2008/09 Antarctic summer, 59 stations (≈ ¾ of the planned footprint) were operational, giving IceTop a total area of about 1 km² and an atmospheric depth of 690 g cm⁻².

The paper first outlines the mechanical design: tanks of 0.9 m radius and 1 m height are made of white‑diffusive plastic, filled with ultra‑pure water that is slowly frozen to avoid air bubbles. Each tank houses two DOMs identical to those used in the deep IceCube array, equipped with 10‑inch Hamamatsu R7081‑2 photomultipliers, onboard high‑voltage supplies, flashers, and 300 MHz digitizers. The high‑gain DOM operates at a gain of 5 × 10⁶, the low‑gain at 10⁵, allowing seamless coverage of the full signal range.

Calibration proceeds in two stages. First, a muon calibration run without shower triggers uses vertical atmospheric muons to define the VEM scale. Small muon telescopes placed above selected tanks provide position‑dependent response, and the resulting charge distributions are fitted to obtain a conversion from photo‑electrons to VEM (≈ 170 PE for a typical DOM). Second, the high‑ and low‑gain channels are cross‑checked within each tank; their combined charge spectra are required to join smoothly, ensuring that signals above 10 VEM are correctly recorded by the low‑gain channel. A continuous muon calibration system is under development to maintain uniform response over the detector’s lifetime.

Triggering relies on a Simple Multiplicity Trigger (SMT) that fires when six DOMs fire within a 5 µs window. Upon an SMT, the full IceTop and InIce data are read out. Because of limited satellite bandwidth, only a subset of events is transmitted: at least three stations must participate for an STA‑3 trigger (≈ 12 Hz in 2008) and eight stations for STA‑8 (≈ 0.5 Hz). Coincidences with the deep IceCube SMT are also sent, with a 5 % transmission rate for InIce‑only triggers that have surface activity.

The first physics results were obtained with the 26‑station configuration (IT26) in 2007. The lateral distribution of the shower signal in VEM was fitted with a function S(r)=S_ref (r/R_ref)^{‑β‑κ log₁₀(r/R_ref)}; simulations indicated κ≈0.3. This form differs from the classic Greisen electron density because IceTop tanks also record gamma‑ray and muon components, which have broader lateral spreads. Using the signal at 125 m (S₁₂₅) as an energy estimator, a response matrix was built for proton and iron primaries. At 3 PeV the core position resolution was 9 m and the logarithmic energy error 0.05 for zenith angles below 30°. Assuming a mixed composition, the all‑particle spectrum shows a knee at 3.1 ± 0.3 (stat) ± 0.3 (syst) PeV, with spectral indices 2.71 ± 0.07 below and 3.11 ± 0.01 above the knee. The absolute normalization is slightly lower than most previous measurements, a discrepancy attributed to the early stage of the IceTop Monte‑Carlo code, which has since been substantially improved.

IceTop offers several composition‑sensitive observables: (1) the energy deposited by the muon bundle in the 1 km of ice (iron showers deposit ≈ 2.4 × more energy than protons at 10¹⁷ eV); (2) the shape of tank waveforms, which can be deconvolved to estimate the number of GeV muons far from the core; (3) the flatness of the lateral distribution (heavier primaries produce flatter LDFs); (4) the rise time of the shower front, related to X_max; (5) the zenith‑angle dependence of surface rates (heavy nuclei are absorbed more quickly at large angles); and (6) the angular distribution of in‑ice muon bundles in non‑coincident events (proton showers generate higher‑energy muons that survive to large zenith angles). By combining these complementary measurements, IceTop aims to reduce model‑dependent uncertainties in composition studies.

Future work will exploit the full 80‑station array, which will provide an acceptance of roughly 1/3 km² sr for coincident surface‑in‑ice events. Expected event rates are of order 10⁶ yr⁻¹ at 1 PeV, rising to ≈ 10⁴ yr⁻¹ at 100 PeV, and allowing the study of the transition from galactic to extragalactic cosmic rays around 10¹⁸ eV. The collaboration plans to refine energy calibration using the correlation between total in‑ice photo‑electron count (NPE) and the surface‑reconstructed shower energy, to develop robust muon‑bundle energy‑loss models, and to use high‑resolution waveform analysis for muon identification. Together with the deep IceCube detector, IceTop will deliver precise measurements of the cosmic‑ray spectrum and composition from 10¹⁵ eV up to the EeV range, providing critical input for astrophysical models of particle acceleration and propagation.