Search for Annual Modulation in Combined ANAIS 112 and COSINE 100 Data using Bayesian Model Comparison

Met behulp van Bayesiaanse modelvergelijking op ongeveer drie jaar gecombineerde ANAIS-112- en COSINE-100-gegevens concludeert de studie dat er geen bewijs is voor een door donkere materie veroorzaakt jaarlijks modulatiesignaal, aangezien de Bayes-factor sterk de voorkeur geeft aan een constant achtergrondmodel boven een sinusvormig model.

De kern

Het Grote Mysterie: Wat is Donkere Materie?

Stel je voor dat ons sterrenstelsel een gigantische, onzichtbare oceaan is. We weten dat deze oceaan bestaat omdat de schepen (sterren) die erdoorheen varen zich zo gedragen dat het lijkt alsof ze worden weggeduwd door iets dat we niet kunnen zien. Wetenschappers noemen dit onzichtbare materiaal Donkere Materie.

Al decennia lang is de leidende theorie dat deze donkere materie bestaat uit tiny, spookachtige deeltjes die WIMPs worden genoemd (Weakly Interacting Massive Particles). Deze deeltjes zijn zo verlegen dat ze zelden tegen normaal materiaal botsen, maar als ze dat wel doen, kunnen ze een klein "krasje" of signaal achterlaten.

Het Geschil: Het Signaal van de "Jaarlijkse Modulatie"

In 1998 claimde een experiment genaamd DAMA/LIBRA deze krasjes te hebben gevonden. Ze zeiden dat ze een specifiek patroon zagen: het aantal "krasjes" ging eens per jaar omhoog en omlaag, zoals een hartslag.

De Analogie: Stel je voor dat je op een strand staat. Je merkt dat elk jaar, rond juni, de golven iets groter worden, en rond december worden ze kleiner. Je concludeert: "De oceaan ademt!" Dit is wat DAMA beweert: de Aarde beweegt door een "wind" van donkere materie, en naarmate we om de Zon draaien, raken we in juni meer ervan en in december minder.

De Skeptici: COSINE-100 en ANAIS-112

Andere wetenschappers waren sceptisch. Ze bouwden hun eigen detectoren, COSINE-100 (in Korea) en ANAIS-112 (in Spanje), met exact dezelfde materialen als DAMA, om te zien of ze hetzelfde "oceaanademen" konden horen.

Lange tijd zeiden deze nieuwe experimenten: "We horen niets." Ze keken naar de gegevens en zagen geen duidelijk omhoog-en-omlaag patroon.

De Nieuwe Studie: Een Vers Nieuwe Blik met een Nieuw Hulpmiddel

In dit artikel besloten twee onderzoekers (Om Godse en Shantanu Desai) om een frisse blik te werpen op de gecombineerde gegevens van COSINE-100 en ANAIS-112 (ongeveer drie jaar aan gegevens).

In plaats van alleen naar de cijfers te kijken, gebruikten ze een statistisch hulpmiddel genaamd Bayesiaanse Modelvergelijking.

De Analogie: Stel je voor dat je een detective bent die een zaak probeert op te lossen. Je hebt twee theorieën:

1. Theorie A (Het Cosine-model): De golven zijn ritmisch. Er is een echt patroon (het "oceaanademen").
2. Theorie B (Het Constante model): De golven zijn gewoon ruis. Er is geen patroon; het is gewoon een vlakke, saaie lijn.

Meestal kijk je misschien gewoon naar de gegevens en gok je welke theorie beter past. Maar dit artikel gebruikt een "Wiskundige Weegschaal" (de Bayes-factor) om het bewijs te wegen. Het vraagt: "Hoe veel waarschijnlijker is het dat de gegevens afkomstig zijn van een ritmisch patroon in vergelijking met willekeurige ruis?"

De Resultaten: De Weegschaal Zette Niet

De onderzoekers testten deze weegschaal met verschillende "regels" (priors) om te zien of het antwoord veranderde afhankelijk van hoe ze de wiskunde instelden. Ze keken naar de gegevens in twee verschillende energiebereiken (alsof je de oceaan bekijkt met twee verschillende gekleurde filters).

Het Oordeel:
Het resultaat was duidelijk: de weegschaal bewoog nauwelijks.

• De "Wiskundige Weegschaal" gaf een score van minder dan 1,15.
• In de wereld van de statistiek is een zo lage score als het ware een zwak gefluister in een lawaaierige kamer. Het is "nauwelijks vermeldenswaard".

Wat dit betekent: De gegevens van COSINE-100 en ANAIS-112 ondersteunen niet het idee dat er een ritmisch "oceaanademen" is (jaarlijkse modulatie). De gegevens lijken veel meer op willekeurige ruis (een constante waarde) dan op een signaal van donkere materie.

De Conclusie

Het artikel concludeert dat wanneer je de gegevens van deze twee experimenten combineert en ze door deze strenge wiskundige test haalt, er geen bewijs is voor het signaal van de jaarlijkse modulatie dat DAMA beweert te zien.

De onderzoekers hebben ook hun "recept" (de computercode) met het publiek gedeeld, zodat iedereen hun werk kan controleren. Dit bevestigt de eerdere bevindingen dat het "oceaanademen" misschien gewoon een lichtbedrog is en geen echt fenomeen.

Technische samenvatting

Technische Samenvatting: Bayesiaanse Modelvergelijking van Gecombineerde ANAIS-112- en COSINE-100-Data

Probleemstelling
Al meer dan twee decennia rapporteert het DAMA/LIBRA-experiment een statistisch significant jaarlijks modulerend signaal in de residuele gebeurtenisfrequenties van natriumjodide (NaI)-detectoren, een kenmerk dat consistent is met interacties van zwakker interacterende zware deeltjes (WIMP) als donkere materie. Deze claim blijft echter onbevestigd door andere directe detectie-experimenten die hetzelfde doelwitmateriaal gebruiken, namelijk COSINE-100 (Korea) en ANAIS-112 (Spanje). Waar eerdere gezamenlijke analyses van de eerste drie jaar aan gecombineerde data van COSINE-100 en ANAIS-112 (verwezen als C25) frequentistische regressie gebruikten om geen bewijs voor modulatie te vinden, tracht deze studie die bevindingen aan te vullen via een Bayesiaans raamwerk. Het kernprobleem dat wordt aangepakt, is of de gecombineerde residuen van deze twee experimenten statistisch bewijs leveren dat een sinusvormig modulerend model (verwacht van donkere materie) de voorkeur verdient boven een constant achtergrondmodel, onafhankelijk van de specifieke regressietechnieken die in eerdere werken werden gebruikt.

Methodologie
De auteurs voeren een Bayesiaanse modelvergelijking uit tussen twee hypothesen:

• Model $M_2$ (Cosinesignaal): De residuale data volgt een sinusfunctie $R(t) = S_m \cos(\omega(t - t_0))$, waarbij $S_m$ de amplitude is, $\omega$ de hoekfrequentie en $t_0$ de fase.
• Model $M_1$ (Constante Achtergrond): De residuale data wordt beschreven door een constante waarde $A$.

De analyse maakt gebruik van ongeveer drie jaar aan gecombineerde data van COSINE-100 en ANAIS-112, met een totale blootstelling van 485 kg-jaar. De data is ingedeeld in 15-daagse intervallen over twee energiebereiken: 1–6 keV en 2–6 keV.

Belangrijke methodologische stappen omvatten:

1. Constructie van de Likelihood: Er wordt een Gaussische likelihood verondersteld voor de waargenomen residuen, waarbij de fouten als Gaussisch worden behandeld.
2. Selectie van Priors: Om robuustheid te waarborgen, testen de auteurs vier verschillende sets priors voor de vrije parameters ($S_m$, $\omega$, $t_0$):
   • Amplitude ($S_m$): Uniforme priors variërend van de minimum tot de maximum waargenomen residuen, en afzonderlijk van nul tot de maximum absolute residu.
   • Frequentie ($\omega$) en Fase ($t_0$): Zowel uniforme priors (die de volledige duur van de dataset respectievelijk 0–365 dagen bestrijken) als normale priors gecentreerd op de best-fit waarden gerapporteerd door DAMA/LIBRA ($T \approx 1$ jaar, $t_0 \approx 145$ dagen).
3. Berekening van het Bewijs: Het Bayesiaanse bewijs (marginale likelihood) voor beide modellen wordt berekend met behulp van het geneste sampling-algoritme geïmplementeerd in Dynesty. De vergelijkingsmetriek is de natuurlijke logaritme van de Bayes-factor, $\ln(B_{21})$, waarbij $B_{21}$ de verhouding is van het bewijs voor het cosinusmodel tot het constante model.

Belangrijkste Bijdragen

• Onafhankelijke Verificatie: Dit werk biedt een onafhankelijke validatie van de C25-resultaten met behulp van een Bayesiaanse modelvergelijkingsbenadering in plaats van de kleinste-kwadraten- en Markov Chain Monte Carlo (MCMC)-regressie die in de oorspronkelijke gezamenlijke analyse werd gebruikt.
• Robuustheid van Priors: Door meerdere prior-configuraties te testen, waaronder die geïnformeerd door de best-fit waarden van DAMA/LIBRA, toont de studie aan dat de conclusie omtrent het ontbreken van een signaal robuust is tegen variaties in prior-aannames.
• Open Wetenschap: De auteurs hebben hun analysecodes publiek beschikbaar gesteld om transparantie en reproduceerbaarheid van de berekening van het Bayesiaanse bewijs te waarborgen.

Resultaten
De analyse levert de natuurlijke logaritme van de Bayes-factor ($\ln(B_{21})$) op voor alle vier prior-sets over beide energie-intervallen. De resultaten worden als volgt samengevat:

• Voor alle geteste prior-combinaties is de waarde van $\ln(B_{21})$ kleiner dan 1,15.
• Volgens de schaal van Jeffreys voor het interpreteren van Bayes-factoren, geeft een waarde van $\ln(B_{21}) < 1,15$ (corresponderend met $B_{21} \leq \sqrt{10}$) aan dat het bewijs voor het cosinusmodel ten opzichte van het constante model "amper vermeldenswaardig" is.
• Specifiek, in de 1–6 keV- en 2–6 keV-intervallen, variëren de berekende waarden van ongeveer -0,13 tot 1,36, waarbij de meerderheid zich nabij of onder de drempel van 1,15 concentreert.
• Gemarginaliseerde postieure contouren voor de amplitude, frequentie en fase tonen geen sterke voorkeur voor de door DAMA-bevoordeelde parameters boven de nulhypothese van een constante snelheid.

Betekenis en Claims
De paper concludeert dat de Bayesiaanse modelvergelijking geen bewijs levert om de hypothese te ondersteunen dat de gecombineerde ANAIS-112- en COSINE-100-data een jaarlijks modulerend signaal bevat dat consistent is met interacties van donkere materie. De resultaten zijn consistent met de hypothese van een constante achtergrond.

De auteurs formuleren hun bijdrage bescheiden als een aanvulling op de C25-analyse. Zij stellen dat hun bevindingen de conclusie versterken dat de gecombineerde data uit de eerste drie jaar van deze experimenten de claim van DAMA over jaarlijkse modulatie niet ondersteunt. De studie claimt niet om interacties van donkere materie in het algemeen definitief uit te sluiten, maar vindt specifiek geen statistische voorkeur voor een sinusvormig signaal boven een constante achtergrond in deze specifieke dataset, ongeacht de prior-aannames met betrekking tot de fase en frequentie van het signaal.