Analytische kernen voor efficiënte constant-Q-transformaties in zoektochten naar donkere materie met LIGO

Naar donkere materie luisteren op een nieuwe manier

Zwaartekrachtsgolfobservatoria zoals LIGO behoren tot de gevoeligste instrumenten ooit gebouwd en kunnen ook krachtige antennes voor donkere materie zijn. Het volledig benutten van hun data werd echter beperkt door een fundamenteel rekenprobleem: de beste manier om naar bepaalde signalen van donkere materie te zoeken is zo kostbaar dat ze onpraktisch wordt voor echte datasets. Dit artikel presenteert een nieuwe signaalverwerkingsmethode die de volledige gevoeligheid van de ideale benadering behoudt terwijl de rekenkosten drastisch dalen, waardoor grondiger zoektochten naar donkere materie in huidige en toekomstige detectoren mogelijk worden.

Waarom donkere materie een smal muzikaal spoor achterlaat

De auteurs richten zich op een populaire klasse ideeën waarbij donkere materie zich gedraagt als een zacht oscillerend veld dat de ruimte vult. In dit plaatje veroorzaakt donkere materie, in plaats van zeldzame deeltjesbotsingen, kleine bijna continue rimpels die natuurconstanten of optische componenten in zwaartekrachtsgolfdetectoren verstoren. Deze rimpels verschijnen als extreem smalle pieken in het frequentiespectrum van de detectoruitgang. De beweging van de aarde door onze Melkweg vervaagt echter elke piek iets, zodat de beste observeertijd afhangt van de frequentie: lage tonen blijven coherent gedurende vele uren, terwijl hoge tonen binnen enkele minuten veranderen. Elke succesvolle zoekmethode moet zich aanpassen aan deze wisselende coherentie-tijd over een zeer breed frequentiebereik.

De uitdaging van inzoomen over veel toonhoogtes

Standaardtools zoals de snelle Fouriertransformatie knippen data in even grote stukken en werken goed wanneer voor alle frequenties dezelfde tijdspanne passend is. Voor ultralichte donkere materie faalt die aanname. Een geschiktere tool is de “constant Q”-transformatie, oftewel logaritmische vermogensspectraliteitsdichtheid, die het tijdvenster voor elke frequentieband aanpast zodat elk deel van het spectrum optimaal wordt behandeld. Helaas schaalt een eenvoudige implementatie met het kwadraat van de datalengte, waardoor ze duizenden tot miljoenen keren langzamer wordt dan snelle algoritmes en feitelijk onbruikbaar voor lange LIGO-datasets. Vorige zoektochten naar donkere materie vertrouwden daarom op slimme benaderingen, waarbij frequenties in banden werden gegroepeerd met vaste vensters en kleine gevoeligheidsverliezen en extra naverwerkingsstappen werden geaccepteerd.

Een snelkoppeling geïnspireerd door digitale muziek

Geput uit technieken uit de computermuziekanalyse herschrijven de auteurs de constant-Q-transformatie zodat het zware werk in het frequentiedomein plaatsvindt in plaats van in het tijdsdomein. Zij scheiden de berekening in de feitelijke detectordata en een wiskundige kernel die codificeert hoe elke frequentiebak wordt gewogen. Terwijl deze kernel in het tijdsdomein breed en kostbaar is om te behandelen, is het tegenbeeld in frequentie scherp gepiekt: slechts een handvol waarden is van belang en de rest kan praktisch als nul worden behandeld. Door deze sparsity te benutten ontwerpen ze een “nul-onderdrukte” versie van de transformatie die het exacte resultaat behoudt maar vrijwel alle onnodige bewerkingen vermijdt. Een sleutelvooruitgang is dat zij een analytische vorm voor de kernel afleiden, zodat deze nooit vooraf berekend of opgeslagen hoeft te worden voor miljoenen frequentiebakken.

Snelheid omzetten in sterkere limieten

Met dit nieuwe raamwerk is één enkele snelle Fouriertransformatie van de data voldoende om alle frequentiebakken van het logaritmische spectrum van gegevens te voorzien, waarna alleen lichte, zeer selectieve bewerkingen nodig zijn. Het team past de methode toe op LIGO’s derde observeerperiode en herziet data die eerder met een benaderende methode waren bestudeerd. Ze vinden dat de nieuwe methode de signaal-tegen-ruisverhouding tot de theoretische maximum versterkt terwijl de rekenkosten dalen, en bereikt ongeveer een orde van grootte sneller dan de eerdere snelle-Fourier-gebaseerde benadering en presteert veel beter dan een brute-forceberekening. Met gedetailleerde modellen van de detectorachtergrond gebaseerd op flexibele splinefits en scheve statistische verdelingen zoeken zij naar overtollig vermogen dat zou wijzen op scalaire-veld donkere materie en verkrijgen in plaats daarvan aangescherpte bovengrenzen voor de mogelijke koppelingen daarvan.

Wat dit betekent voor toekomstige observatoria

Hoewel in deze studie geen signaal van donkere materie wordt gevonden, is de methode op zichzelf een krachtig nieuw instrument. Elk experiment dat een logaritmisch spectrum nodig heeft afgestemd op wisselende coherentie-tijden — van grondgebaseerde interferometers zoals LIGO en GEO600 tot geplande ruimtemissies zoals LISA — kan nu volledig optimale analyses uitvoeren zonder onoverkomelijke rekenkosten. Door het meest gevoelige type zoekopdracht op schaal haalbaar te maken, vergroot dit werk de kans dat toekomstige zwaartekrachtsgolfdetectoren niet alleen verre kosmische botsingen horen maar ook het zwakke, constante gezoem van donkere materie zelf opvangen.

Bronvermelding: Göttel, A.S., Raymond, V. Analytical kernels for efficient constant Q transforms in dark matter searches with LIGO. Sci Rep 16, 15364 (2026). https://doi.org/10.1038/s41598-025-33428-2

Trefwoorden: LIGO, donkere materie, zwaartekrachtsgolven, signaalverwerking, constant Q-transformatie