WavePID: Low-energy flavor identification using single-PMT time series in IceCube

Technische Samenvatting: WavePID: Laag-energetische flavorenidentificatie met behulp van single-PMT tijdreeksen in IceCube

Probleemstelling
De IceCube Neutrino Observatorie vertrouwt op de morfologie van gebeurtenissen om onderscheid te maken tussen neutrino-flavoren (tracks versus cascades). Echter, in het laag-energetische regime van 5–100 GeV, specifiek binnen de DeepCore sub-array, vertonen gebeurtenissen vaak een ijle fotondetectie (specifiek $10 \le n_{hit} < 20$). In dit regime biedt het lichtpatroon zwakke beperkingen aan de gebeurtenis-morfologie, wat de scheidingskracht van state-of-the-art graph neural networks (GNNs) zoals DynEdge beperkt. Hoewel GNNs gebruikmaken van ruwe per-puls informatie (positie, tijd, lading), normaliseren zij tijdstempels naar de microseconde-schaal van de gebeurtenis, wat potentieel de cruciale nanoseconde-schaal vroege-tijdstructuur kan comprimeren die muon-tracks onderscheidt van elektromagnetische cascades. Het artikel behandelt de vraag of nanoseconde-schaal timinginformatie op individuele Digital Optical Modules (DOMs) flavoren-identificatie-informatie bevat die complementair is aan bestaande morfologie-gebaseerde classificaties.

Methodologie
De auteurs introduceren WavePID, een compacte, op templates gebaseerde log-likelihood-ratio classifier ontworpen om nanoseconde-schaal timing op individuele DOMs uit te buiten. De methodologie verloopt in drie stadia:

1. Microfysica Motivatie (Geant4 Simulatie): Gebruikmakend van het Geant4 framework en het IceCube OMSim model, simuleerden de auteurs foton aankomsttijdverdelingen (PATDs) voor 5–100 GeV muon en elektron in ijs. Zij vonden dat muon-tracks een scherpe, directe Cherenkov-voortplantingsgolf produceren, wat resulteert in een snel stijgende PATD gedomineerd door primaire deeltjesemissie. Omgekeerd produceren elektromagnetische cascades (elektronen) licht via Bremsstrahlung secundaire deeltjes over meerdere nanoseconden, wat resulteert in een bredere PATD. Dit verschil is het meest uitgesproken in de eerste $\sim10$ ns na de aankomst van het eerste foton.
2. Observabele Constructie: WavePID werkt op gereconstrueerde ATWD-pulsen van individuele DOMs. Voor elke DOM $n$ construeert het een feature vector $\mathbf{x}_n$ bestaande uit drie natuurkundig gemotiveerde observabelen:
   • $r_n$: De afstand van de DOM tot de gereconstrueerde interactie-vertex.
   • $\xi^t_n$: De vroege-lading fractie (early-charge fraction), gedefinieerd als het fractie van de totale lading die binnen de eerste $t$ nanoseconden (specifiek $t=14$ ns) na de eerste puls arriveert.
   • $\Delta t_n$: De tijd van de eerste puls op de DOM relatief aan de vroegste puls op enige DOM in de gebeurtenis.
3. Classificatie Algoritme: De feature-ruimte wordt gediscretiseerd in 3D-bins. Het algoritme construeert vorm-alleen templates (waarschijnlijkheidsmassafuncties) voor "track" en "cascade" hypothesen met behulp van gelabelde Monte Carlo (MC) simulaties. De teststatistiek (TS) is een gebinde log-likelihood ratio:
   $$TS_e = \log \frac{\mathcal{L}_{cascade}}{\mathcal{L}_{track}} = \sum_{i,j,k} H^e_{ijk} \log \frac{\lambda^{cascade}_{ijk}}{\lambda^{track}_{ijk}}$$
   waarbij $H^e_{ijk}$ de geobserveerde bezetting is van DOMs in een specifieke bin voor gebeurtenis $e$, en $\lambda$ de template-waarschijnlijkheden representeert.

Belangrijkste Resultaten
De prestaties van WavePID werden geëvalueerd op IceCube DeepCore Monte Carlo ensembles die de dominante detector en ijs-model onzekerheden vertegenwoordigen (31 systematische variaties). De studie richtte zich op een sample verrijkt met cascades (geselecteerd door een DynEdge score $\le 0.4$) met hit-multipliciteiten tussen 10 en 20.

• Prestatieverbetering: In het ijle hit-regime ($10 \le n_{hit} < 20$) bereikte WavePID een hogere Area Under the Curve (AUC) van 0.571 vergeleken met 0.542 voor DynEdge alleen binnen dezelfde voorgeselecteerde sample. Dit demonstreert dat de vroege-tijd pulsstructuur flavoren-informatie codeert die niet wordt gevangen door de GNN.
• Zuiverheidswinst: Het toepassen van een drempelwaarde op de WavePID teststatistiek verhoogde de cascade-zuiverheid. Bij een werkpunt dat 20% van de cascades behoudt, steeg de zuiverheid van de baseline van 0.66 (DynEdge alleen) naar 0.75. Deze winst van $+7\% \pm 1\%$ bleef stabiel over alle 31 detector systematische ensembles, wat wijst op robuustheid tegen onzekerheden in DOM-efficiëntie, hole ice, en bulk ijs optische eigenschappen.
• Systematiek: De scheidingskracht bleef stabiel over variaties in detectorparameters, wat bevestigt dat het signaal wordt gedreven door fysieke timingverschillen in plaats van instrumentele artefacten.

Betekenis en Claims
Het artikel claimt dat WavePID de levensvatbaarheid aantoont van het gebruik van nanoseconde-schaal single-PMT timing als een op eerste principes gebaseerde observabele voor flavorenidentificatie in neutrino-telescopen. De auteurs benadrukken dat dit naar hun weten de eerste expliciete toepassing is van single-DOM nanoseconde timing voor deeltjesidentificatie in een neutrino-telescoop.

De betekenis ligt in de complementariteit van de methode: WavePID vervangt geen morfologie-gebaseerde classifiers, maar voegt een nieuwe dimensie van informatie toe. De discriminatie is geconcentreerd in het regime waar de primaire lepton leading Cherenkov-front wordt opgelost door per-DOM pulstiming, een kenmerk dat bestaande classifiers mogelijk missen door temporele normalisatie.

De auteurs concluderen dat deze resultaten motiveren om de exploitatie van fijnmazige per-module timing te gebruiken in toekomstige laag-energetische classifiers en in het ontwerp van volgende-generatie neutrino-telescopen (zoals IceCube-Gen2, KM3NeT/ORCA, en anderen) die gebruikmaken van multi-PMT of hybride optische modules met verbeterde timingresolutie. Het artikel merkt op dat de gepresenteerde configuratie een baseline implementatie is, en dat de optimale keuze van observabelen en binning afhangt van specifieke natuurkunde-doelstellingen. Toepassing op echte IceCube-data en kwantificering van de natuurkundige impact zijn voorbehouden aan toekomstige dedicated studies.