GPU-accelerated spectrum reweighting for new-physics searches in solar neutrino--electron scattering

Technische Samenvatting: GPU-versnelde Spectrum-herweging voor Zonne-neutrino–elektronverstrooiing

Probleemstelling
Precieze metingen van de elastische verstrooiing van zonne-neutrino’s op elektronen dienen als gevoelige sondes voor fysica buiten het Standaardmodel (BSM), waaronder niet-standaard interacties (NSI's) en afwijkende neutrino-magnetische momenten. Deze nieuwe fysica-effecten modificeren het differentiële dwarsdoorsnedeproces, waardoor zowel de gebeurtenisratio als de spectrale vorm van de terugstuwende elektronen worden gewijzigd. In conventionele likelihood-analyses vereist het testen van een nieuw fysica-hypothese het propageren van een gemodificeerd terugstuwend spectrum door een detectorresponsmodel dat is afgeleid van volledige Monte Carlo (MC)-simulaties (bijv. Geant4). Dit proces is computationeel prohibitief wanneer het wordt uitgevoerd voor elk punt in een meerdimensionale parameterruimte, aangezien het regenereren van detector-niveau MC-monsters voor elke variatie in parameters de kosten van fits, scans en Markov-keten Monte Carlo-studies domineert.

Methodologie
De auteurs stellen een GPU-versneld spectrum-herwegingsframework voor dat de variatie van de fysica-parameters ontkoppelt van de detectorsimulatie. De kern van de methodologie berust op de volgende stappen:

1. Factorisatie-aanname: Het framework gaat ervan uit dat de detectorrespons voor een terugstuwend elektron met een bepaalde kinetische energie $T$ onafhankelijk is van de onderliggende interactie (Standaardmodel, NSI of magnetisch moment). Bijgevolg volstaat een enkele baseline MC-steekproef, gegenereerd onder de Standaardmodel (SM) hypothese, om de mapping tussen de werkelijke terugstuwende energie en gereconstrueerde observabelen (bijv. foto-elektronen) te definiëren.
2. Bin-naar-bin Herweging: In plaats van MC-monsters te regenereren, past de methode multiplicatieve gewichten toe op de bins van het terugstuwende energiespectrum. Voor een set nieuwe fysica-parameters $\varepsilon$ wordt het terugstuwende spectrum berekend als $f_T(T; \varepsilon) = r(T; \varepsilon) \cdot \hat{f}^{SM}_T(T)$, waarbij $r(T; \varepsilon)$ de herwegingsfunctie is die is afgeleid van de gewijzigde differentiële dwarsdoorsneden.
3. Vaste Respons-propagatie: Het hergewogen spectrum wordt gepropageerd door een vooraf berekende, tweedimensionale detectorresponsmatrix $R_{ij}$ (die de mapping vormt tussen de werkelijke terugstuwende bin $i$ en de gereconstrueerde observabele bin $j$). Deze matrix, afgeleid van Geant4, bevat effecten zoals scintillatie-quenching, energie-resolutie en reconstructie-efficiëntie.
4. GPU-implementatie: De berekening is geïmplementeerd met behulp van NVIDIA's Thrust template library, waarbij gebruik wordt gemaakt van parallelle primitieven (transform, reduce, transform-reduce). De code is geschreven in een enkele bron die compileert voor zowel CUDA (GPU) als OpenMP (CPU) backends. De pijplijn omvat:
   • Fysica-update: Het berekenen van oscillatiekansen en differentiële dwarsdoorsneden op een gediscretiseerd kinematisch rooster.
   • MC Herweging: Het toepassen van de herwegingsfunctie op het terugstuwende spectrum en het convolueren ervan met de vaste responsmatrix.
   • Likelihood Evaluatie: Het accumuleren van het voorspelde spectrum en het evalueren van de Poisson log-likelihood.

Belangrijkste Bijdragen

• Algoritmische Efficiëntie: De primaire versnelling wordt bereikt door het vermijden van de regeneratie van detector-niveau MC-monsters voor elk parameterpunt. Dit reduceert de computationele kosten van een parameter-update van uren (volledige simulatie) naar milliseconden (herweging).
• Hardware-agnosticisme: Het framework is ontworpen om te draaien op zowel gangbare consumenten-GPU's (bijv. NVIDIA RTX 3080 Ti) als datacenter-accelerators (bijv. NVIDIA A30X), evenals op multi-core CPU's, met behulp van een uniforme codebase.
• Validatie: De methode wordt gevalideerd tegen een analytische referentie waarbij de exacte oplossing bekend is. De herwegingsvoorspelling reproduceert de analytische waarheid met een integrale afwijking onder de $10^{-6}$ voor een 2000-bin terugstuwend rooster, wat bevestigt dat de bin-niveau benadering een verwaarloosbare fout introduceert ten opzichte van de discretisatieschaal.

Resultaten

• Prestatiebenchmarks:
  • Op een datacenter NVIDIA A30X GPU duurt een enkele likelihood-evaluatie ongeveer 52 ms. Dit vertegenwoordigt een 58× versnelling ten opzichte van een enkele CPU-thread en een 2,5× versnelling ten opzichte van een volledig belaste 64-thread CPU.
  • Op een consumenten NVIDIA RTX 3080 Ti duurt de evaluatie ongeveer 87 ms, wat interactieve parameter-scans op een enkele werkstation mogelijk maakt.
• Nauwkeurigheid: Het framework slaagt erin geïnjecteerde parameters te herstellen in end-to-end "injectie-reconstructie" tests voor zowel flavor-diagonale NSI's als neutrino-magnetische momenten. De steile $1/T$ afhankelijkheid van het magnetische moment dwarsdoorsnedeproces, die een strenge test vormt voor de nauwkeurigheid van de binning, wordt met hoge precisie afgehandeld.
• Schaalbaarheid: De GPU-implementatie schaalt efficiënt; de computationele bottleneck is de herhaalde toepassing van de respons-kernel, die van nature geparallelliseerd is over de gereconstrueerde energie-bins.

Betekenis en Claims
Het artikel beweert dat dit framework interactieve parameter-scans mogelijk maakt voor analyses van zonne-neutrino–elektronverstrooiing op hardware die toegankelijk is voor individuele onderzoekers. De auteurs benadrukken dat de belangrijkste versnelling voortkomt uit de algoritmische strategie van het vermijden van MC-regeneratie in plaats van enkel door GPU-hardwareversnelling. Hoewel de GPU een verdere versnelling biedt ten opzichte van multi-core CPU's, ligt de waarde van de methode in het mogelijk maken van hoogwaardige detector-modellering binnen iteratieve fitting-loops die voorheen computationeel onhaalbaar waren.

De reikwijdte is expliciet beperkt tot scenario's waarin nieuwe fysica kan worden gerepresenteerd als een herweging van een bestaand terugstuwend spectrum (bijv. flavor-diagonale NSI's en magnetische momenten), terwijl de detectorrespons constant blijft. De auteurs merken op dat scenario's met kwalitatief andere eindtoestanden-deeltjes of flavor-veranderende structuren een aanvullende formalisering vereisen. Het werk dient als een methodologische studie, gevalideerd op synthetische analytische responsen en Geant4-afgeleide matrices, zonder gebruik te maken van specifieke experimentele datasets.