Modeling the refractive index profile n(z) of polar ice for ultra-high energy neutrino experiments

Titel: Modellierung des Brechungsindexprofils $n(z)$ von polarem Eis für Experimente mit ultrahohen Neutrinoenergien

Autoren: S. Ali et al. (ARA-Kollaboration)
Datum: 2. April 2026 (vorgelegt am 1. April 2024)

---

1. Problemstellung

Experimente zur Detektion von Ultrahoch-Energie-Neutrinos (UHEN), wie das Askaryan Radio Array (ARA) am Südpol, die Radio Neutrino Observatory in Greenland (RNO-G) und das geplante IceCube-Gen2, nutzen Radiowellen, um Neutrino-Wechselwirkungen im Eis zu detektieren. Da Radiosignale im Eis über Kilometer weit propagieren können, ist dies eine kosteneffiziente Methode, um große Detektionsvolumina zu erschließen.

Ein kritischer Faktor für die Sensitivität dieser Detektoren ist das Brechungsindexprofil $n(z)$ des Eises in Abhängigkeit von der Tiefe. Da der Druck mit der Tiefe zunimmt, verdichtet sich das Eis, was den Brechungsindex verändert. Dies führt gemäß dem Fermatschen Prinzip zu einer Krümmung der Radiowellenstrahlen (Refraktion).

• Herausforderung: Eine ungenaue Modellierung von $n(z)$ führt zu Fehlern bei der Rekonstruktion von Neutrino-Ereignissen und zur Unterschätzung oder Überschätzung des „schattierten Volumens" (Shadow Zone), in dem Signale den Detektor nicht erreichen können.
• Aktueller Ansatz: Bisherige Simulationen basierten oft auf einem einfachen einstufigen exponentiellen Modell (ähnlich einem selbstgravitierenden Fluid).
• Physikalische Realität: Glaziologische Studien zeigen jedoch, dass die Verdichtung von Schnee zu Eis in mehreren Phasen erfolgt (Schnee $\to$ Firn $\to$ blasiges Eis $\to$ festes Eis), was durch unterschiedliche Dichtegrenzen und Verdichtungsraten gekennzeichnet ist.

2. Methodik

Die Autoren entwickelten ein in-situ Profil des Brechungsindex, indem sie die Laufzeitunterschiede von Radiosignalen analysierten, die von englischen Sendern zu Empfängern gesendet wurden.

• Datenquelle:
  • IDP (IceCube Deep Pulser): Ein stationärer Sender in 1400 m Tiefe.
  • SDP (SPICE Deep Pulser): Ein Sender, der in den SPICE-Eisbohrkern (bis 1751 m Tiefe) abgesenkt wurde.
  • Empfänger: Die Antennen der ARA-Stationen (A2–A5), die in Tiefen von 170–200 m installiert sind.
• Messgröße: Der Zeitunterschied $\Delta t(D,R)$ zwischen dem direkten Signalpfad ($D$) und dem refraktierten/reflectierten Signalpfad ($R$).
  • Der direkte Pfad durchläuft hauptsächlich das tiefe, blasige Eis.
  • Der refraktierte Pfad durchläuft auch die oberen Schichten (Firn und Schnee), wodurch das Profil $n(z)$ in verschiedenen Tiefen getestet wird.
• Simulationswerkzeug: Der Ray-Tracer RadioPropa (integriert in NuRadioMC) wurde verwendet, um die Strahlengänge basierend auf verschiedenen $n(z)$-Modellen zu simulieren und mit den gemessenen $\Delta t(D,R)$-Werten zu vergleichen.
• Modellvergleich:
  1. Einstufiges Modell: Einfache exponentielle Abhängigkeit ($n(z) = n_{\infty} - B_0 e^{Cz}$).
  2. Dreistufiges Modell: Ein stückweise definiertes Modell, das drei Phasen der Verdichtung abbildet (Schnee, Firn, blasiges Eis) mit unterschiedlichen Verdichtungsraten ($c_1, c_2, c_3$) und Dichtegrenzen ($\rho_1, \rho_2$).

3. Wichtige Beiträge und Ergebnisse

A. Validierung des dreistufigen Modells

Durch Anpassung der Simulationsparameter an die gemessenen $\Delta t(D,R)$-Daten der Station A2 (als Trainingsdatensatz) und Validierung an den Stationen A3, A4 und A5 ergab sich:

• Das dreistufige Modell passt die gemessenen Zeitdifferenzen signifikant besser an als das ein- oder zweistufige exponentielle Modell.
• Statistische Kennzahlen ($\chi^2$ pro Freiheitsgrad):
  • 1-Stufen-Modell: 406,5
  • 2-Stufen-Modell: 295,8
  • 3-Stufen-Modell: 9,3
• Die Diskrepanz bei den ein- und zweistufigen Modellen war besonders bei flacheren Sendertiefen und in der Firn-Region (ca. 20–100 m Tiefe) offensichtlich.

B. Korrelation mit Dichtedaten (SPICE-Kern)

Die Autoren wandelten Dichtedaten aus dem SPICE-Eisbohrkern unter Verwendung einer linearen Beziehung ($n = 1 + A\rho$) in Brechungsindizes um.

• Das dreistufige Modell stimmt in der Firn-Region (zwischen 20,5 m und 96,6 m Tiefe) sehr gut mit den umgerechneten Dichtedaten überein.
• In der Region des blasigen Eises (>96,6 m) zeigt das dreistufige Modell eine leichte Abweichung (Überschreitung) gegenüber den umgerechneten Dichtedaten, was darauf hindeuten könnte, dass die lineare Beziehung in großen Tiefen nicht mehr perfekt gilt oder zusätzliche Daten benötigt werden.

C. Einfluss auf den Schattenbereich (Shadow Zone)

Die Krümmung der Strahlen bestimmt, wo ein Sender für einen Empfänger „unsichtbar" wird (Schattenzone).

• Das dreistufige Modell verändert die Ausdehnung der Schattenzone im Vergleich zum einstufigen Modell.
• Für Empfänger in 100 m Tiefe (typisch für RNO-G) wird die Schattenzone kleiner.
• Für Empfänger in 170 m Tiefe (typisch für ARA) wird die Schattenzone leicht größer, was die Reichweite für tiefere Quellen beeinflusst.

D. Berücksichtigung von Schneeanhäufung

Da die ARA-Antennen zwischen 2011 und 2013 installiert wurden, die SDP-Daten aber 2018–2019 aufgenommen wurden, wurde der Effekt der jährlichen Schneeanhäufung (ca. 27 cm/Jahr an der Oberfläche, aber weniger an der Tiefe durch Verdichtung) analysiert.

• Die Analyse der Zeitverschiebungen über mehrere Jahre ergab eine effektive Tiefenänderung von ca. 0,24 m pro Jahr.
• Diese Korrektur war notwendig, um die Sender- und Empfängertiefen für den Modellvergleich präzise zu bestimmen.

4. Signifikanz und Auswirkungen

• Steigerung der Empfindlichkeit: Simulationen des effektiven Volumens für einen ARA-ähnlichen Detektor (170 m Tiefe) zeigen, dass die Verwendung des dreistufigen Modells im Vergleich zum einstufigen exponentiellen Modell eine Steigerung der Neutrino-Empfindlichkeit um 14 % im Energiebereich von $10^{18}$ bis $10^{21}$ eV bewirkt.
• Kalibrierung und Rekonstruktion: Ein präziseres $n(z)$-Modell verbessert die Genauigkeit der Rekonstruktion von Neutrino-Ereignissen (Richtung und Energie) und reduziert systematische Unsicherheiten bei der Kalibrierung zukünftiger Experimente (z. B. IceCube-Gen2 Radio).
• Glaziologische Erkenntnis: Die Studie liefert eine unabhängige, radio-basierte Bestätigung für die Existenz von drei Verdichtungsphasen im polaren Eis, die mit glaziologischen Modellen übereinstimmen.

Fazit

Die Arbeit demonstriert, dass die Annahme eines einfachen exponentiellen Brechungsindexprofils für polare Eisdetektoren unzureichend ist. Ein glaziologisch motiviertes, dreistufiges Modell, das durch in-situ Radiomessungen validiert wurde, liefert nicht nur eine bessere Übereinstimmung mit den gemessenen Laufzeiten, sondern erhöht auch die nachweisbare Neutrinoflussrate erheblich. Dies ist entscheidend für die Optimierung zukünftiger Radio-Neutrino-Observatorien.