The Missing Multipole Problem: Investigating biases from model starting frequency in gravitational-wave analyses

Dit onderzoek toont aan dat het verkeerd kiezen van de startfrequentie van golfvormmodellen bij zware, asymmetrische zwarte-gatbinaire systemen leidt tot vertekende schattingen van de bronparameters, waarbij de drempel voor deze bias afhangt van de totale massa, de massaverhouding en het signaal-ruisverhouding.

De kern

Het Verdwijnende Deel van de Melodie: Waarom we de start van een geluid moeten begrijpen

Stel je voor dat je naar een orkest luistert dat een prachtige symfonie speelt. Je wilt precies weten welke instrumenten er spelen en hoe hard ze spelen, zodat je de compositie kunt reconstrueren. Maar er is een probleem: je begint pas te luisteren op het moment dat de trompetten hun eerste noot blazen. De violen en de contrabassen, die veel lagere tonen spelen, hebben al een tijdje aan het werk zijn, maar jij hebt die lage, diepe klanken gemist.

Zou je op basis van alleen de trompetten kunnen concluderen wat het hele orkest deed? Waarschijnlijk niet. Je zou denken dat het orkest kleiner is dan het in werkelijkheid is, of dat de muziek anders klinkt dan hij doet.

Dit is precies wat er gebeurt in de wereld van zwaartekrachtgolven (gravitationele golven), en dit is waar dit nieuwe onderzoek over gaat.

1. Wat zijn zwaartekrachtgolven?

Wanneer twee zwarte gaten in het heelal tegen elkaar botsen, maken ze de ruimte zelf trillen. Dit zijn zwaartekrachtgolven. Het zijn als het ware "geluidsgolven" van het universum, maar dan gemaakt van kromme ruimte in plaats van lucht.

Wetenschappers gebruiken deze golven om te achterhalen hoe zwaar de zwarte gaten waren, hoe snel ze draaiden en hoe ver ze weg waren. Om dit te doen, vergelijken ze de gemeten trillingen met miljoenen theoretische modellen (zoals een enorme database van mogelijke melodieën).

2. Het Probleem: De "Verdwijnende" Multipoles

Elke botsing van zwarte gaten produceert een complex geluid. Het belangrijkste geluid komt van de twee hoofdgaten (de "bas" en de "melodie"). Maar er zijn ook subtiele, hogere tonen die meeklinken. In de wetenschap noemen we deze multipoles.

• De basis (2,2): Dit is het hoofdgeluid, de diepe bas.
• De hogere tonen (3,3 en 4,4): Dit zijn de hogere, subtiele harmonieën.

Het probleem is dat deze hogere tonen eerder beginnen dan de bas. Ze zijn als een fluit die al begint te spelen voordat de trompet zijn eerste noot haalt.

In de computermodellen die wetenschappers gebruiken, moeten ze beslissen: "Op welk moment beginnen we met het analyseren van het geluid?"
Meestal kiezen ze voor 20 Hertz (Hz). Dit is het laagste punt dat onze huidige detectors kunnen horen.

Maar hier zit de valstrik:

• Als je start bij 20 Hz, hoor je de bas (2,2) perfect.
• Maar de hogere tonen (3,3 en 4,4) beginnen al bij 10 Hz of 13 Hz.
• Door te wachten tot 20 Hz, snijdt de computer het begin van die hogere tonen gewoon af. Het is alsof je een filmpje start op seconde 10, terwijl de actie al bij seconde 1 begon.

Dit noemen de auteurs het "Missing Multipole Problem" (Het probleem van de verdwijnende multipoles).

3. Wat gebeurt er als je het begin mist?

De auteurs van dit artikel hebben gekeken naar zware zwarte gaten (zogenoemde "Intermediate Mass Black Holes"). Ze hebben simulaties gedaan en ontdekt dat als je het begin mist (door pas bij 20 Hz te starten), je de resultaten van je metingen verkeerd interpreteert.

Het is alsof je probeert het gewicht van een vrachtwagen te schatten door alleen naar de wielen te kijken, maar de wielen die je ziet zijn er net bijgekomen.

• Je denkt dat de zwarte gaten lichter zijn dan ze zijn.
• Je denkt dat ze ongelijk van gewicht zijn, terwijl ze misschien wel even zwaar zijn.
• Je schat de afstand verkeerd in.

Dit gebeurt vooral bij zware systemen en bij signalen die heel sterk zijn (hoge "signaal-ruisverhouding"). Als het signaal zwak is, is de ruis zo groot dat je de fouten niet eens merkt. Maar als het signaal helder is (zoals bij een zeer zware botsing), dan wordt de fout enorm.

4. De Oplossing: Start eerder!

De wetenschappers zeggen: "Stop met wachten tot 20 Hz!"

• Voor lichte signalen: Je kunt veilig wachten tot 20 Hz. De statistische ruis is dan nog zo groot dat de gemiste tonen niet uitmaken.
• Voor zware signalen (SNR > 20): Je moet de analyse starten bij 13 Hz. Dan hoor je de eerste extra toon (3,3) en krijg je een correct beeld.
• Voor zeer zware en asymmetrische signalen (SNR > 70): Je moet zelfs starten bij 10 Hz. Dan hoor je ook de tweede extra toon (4,4).

Het is alsof je zegt: "Als je een zware vrachtwagen wilt wegen, moet je de weegschaal al activeren voordat de eerste wielen erop staan, anders weeg je te licht."

5. Waarom is dit belangrijk?

Ons universum zit vol met zwarte gaten. Sommige botsen zijn zo zwaar dat ze de "gaten" in de massa-verdeling van sterren kunnen vullen (de zogenaamde pair-instability gap). Als we door dit "verdwijnende deel"-probleem de massa's verkeerd meten, kunnen we denken dat een zwart gat in een onmogelijk gebied zit, terwijl het dat eigenlijk niet is.

Dit onderzoek waarschuwt de wereld: We moeten onze modellen aanpassen. We moeten de "startknop" van onze analyses lager zetten (naar 10 of 13 Hz) om de volledige symfonie van het universum te horen, en niet alleen het stukje dat pas later begint.

Samenvatting in één zin

Als je luistert naar de botsing van zware zwarte gaten, mag je niet wachten tot het geluid "hard" genoeg is om te horen (20 Hz), want dan mis je de subtiele, vroege tonen die nodig zijn om te weten wat je eigenlijk ziet; je moet eerder beginnen met luisteren (bij 10 of 13 Hz) om de waarheid te horen.

Technische samenvatting

Titel: Het Vermiste Multipoolprobleem: Onderzoek naar biases door startfrequentie van modellen in gravitatiegolfanalyses

Auteurs: R. Ursell, C. Hoy, I. Harry, L. K. Nuttall (University of Portsmouth, UK)
Publicatie: MNRAS (Preprint 2025)

---

1. Het Probleem: Het "Vermiste Multipool"-Probleem

Gravitatiegolf (GW) signalen van samensmeltende zwarte gaten kunnen worden ontbonden in een som van spin-gewogen sferische harmonischen (multipoles). De quadrupool $(\ell=2, m=\pm2)$ is doorgaans de dominante bijdrage. Echter, hogere-orde multipoles (zoals $(3,3)$ en $(4,4)$) worden significant in systemen met ongelijke massa's, schuine banen en hoge totale massa's.

Het probleem ontstaat bij het gebruik van tijdsdomein-modellen (time-domain models) die worden omgezet naar het frequentiedomein voor analyse. Elke harmonische $(\ell, m)$ begint op een specifieke tijd $t_0$, wat correspondeert met een specifieke startfrequentie $f_{\ell m}$.

• Voor de dominante quadrupool geldt vaak $f_{22} = f_{start}$.
• Voor hogere multipoles geldt $f_{\ell m} \approx \frac{m}{2} f_{22}$.

Als een analyse begint bij de standaard detectorfrequentie van 20 Hz ($f_{start}=20$ Hz), wordt de $(2,2)$-component volledig vastgelegd. Echter, de $(3,3)$-component begint pas bij $\sim30$ Hz en de $(4,4)$-component pas bij $\sim40$ Hz. Dit betekent dat er niet-negligibel vermogen (power) van hogere-orde multipoles in het lage frequentiebereik (20–30 Hz en 20–40 Hz) wordt gemist. Dit fenomeen, het "missing multipole problem", kan leiden tot systematische biases in de geschatte bronparameters (zoals massa's en spins), vooral bij kortdurende signalen met weinig inspiratie.

2. Methodologie

De auteurs voerden een uitgebreide systematische studie uit om de impact van de startfrequentie van templates op Bayesiaanse parameter-schatting te kwantificeren.

• Modellen: Ze gebruikten het NRSur7dq4-model, een numeriek-relativistisch surrogaatmodel dat tot $\ell_{max}=4$ nauwkeurig is en spin-precessie bevat.
• Injectie/Recovery: Ze injecteerden synthetische GW-signalen in "zero-noise" (om systematische fouten te isoleren) met bekende parameters ($\theta_{inj}$).
• Variatie: Ze varieerden de startfrequentie van de template ($f_{22}$) in drie scenario's:
  1. 10 Hz: Alle multipoles tot $\ell=4$ zijn aanwezig vanaf de analyseband (20 Hz). Dit dient als de "onbevooroordeelde" baseline.
  2. 13 Hz: De $(3,3)$-multipool is aanwezig, maar $(4,4)$ mist het lage frequentiebereik.
  3. 20 Hz: Alleen de $(2,2)$-multipool is aanwezig bij 20 Hz; $(3,3)$ en $(4,4)$ missen significant vermogen.
• Parameterbereik: De studie focuste op lichte Intermediaire Zware Black Holes (IMBH) met totale massa's $M \in [200, 450] M_\odot$, massaverhoudingen $q \in [0.25, 1]$, en signaal-ruisverhoudingen (SNR) $\rho \in [20, 100]$.
• Analyse: Bayesiaanse inferentie werd uitgevoerd met de nested sampling algoritme dynesty (via bilby).
• Beoordelingsmetriek: Om biases te kwantificeren, gebruikten ze de Mahalanobis Recovery Score ($r_M$). Deze score meet het percentage van de posterior-verdelingen waarin de ware waarde binnen het 90% betrouwbaarheidsinterval valt. Een lagere score duidt op een grotere bias. Ook werden Bayes-factoren berekend om te bepalen welk model de data beter ondersteunt.

3. Belangrijkste Resultaten

De resultaten tonen aan dat de keuze van de startfrequentie kritiek is, afhankelijk van de SNR en de fysieke eigenschappen van het systeem:

• Invloed van Totale Massa:

  • Voor $M \lesssim 250 M_\odot$ zijn de biases verwaarloosbaar; alle startfrequenties geven vergelijkbare resultaten.
  • Voor $M \gtrsim 300 M_\odot$ ontstaan significante biases als $f_{22}=20$ Hz wordt gebruikt. De geschatte totale massa wordt onderschat (bij $300 M_\odot$) of overschat (bij $>350 M_\odot$), en de massaverhouding wordt vertekend.
  • Bij $M \gtrsim 450 M_\odot$ wordt zelfs de $(4,4)$-multipool essentieel; hierbij divergeert zelfs de $f_{22}=13$ Hz analyse van de baseline.

• Invloed van Signaal-ruisverhouding (SNR):

  • $\rho < 20$: Statistische onzekerheden domineren. Het missen van lage frequenties heeft geen significant effect; $f_{22}=20$ Hz is acceptabel.
  • $20 \lesssim \rho \lesssim 70$: Systematische biases worden zichtbaar. Het missen van de $(3,3)$-multipool (door starten bij 20 Hz) leidt tot fouten. Een startfrequentie van 13 Hz is vereist om de $(3,3)$-multipool volledig te vangen.
  • $\rho \gtrsim 70$: Zelfs het missen van de $(4,4)$-multipole wordt kritiek. Voor hoge SNR's en zware systemen is een startfrequentie van 10 Hz noodzakelijk voor nauwkeurige parameter-schatting.

• Invloed van Massaverhouding en Spin:

  • Bij gelijke massa's ($q=1$) zijn de biases het grootst bij $f_{22}=20$ Hz omdat de inspiratie korter is en hogere multipoles minder tijd hebben om de band binnen te komen.
  • Bij sterk ongelijke massa's en hoge spins neemt de impact van het missen van multipoles toe, hoewel de trends complexer zijn.
  • De inclinatiehoek ($\theta_{JN}$) wordt over het algemeen robuust gerecupereerd, maar bij hoge SNR's en specifieke hoeken (bijv. $\pi/3$) wordt de $(4,4)$-multipool belangrijk.

• Toepassing op Real Data (GW231123):

  • De auteurs pasten hun analyse toe op het reële signaal GW231123 (SNR $\approx 22.6$).
  • Bij gebruik van $f_{22}=20$ Hz werd de totale massa bimodaal en vertoonde meer steun voor ongelijke massa's.
  • Analyses met $f_{22}=10$ Hz en $f_{22}=13$ Hz gaven vergelijkbare resultaten, wat aangeeft dat de $(3,3)$-multipool bij dit signaal belangrijk is, maar de $(4,4)$-multipool minder kritiek (gezien de SNR).
  • De Bayes-factor gaf de voorkeur aan de $f_{22}=10$ Hz analyse ten opzichte van de $f_{22}=20$ Hz analyse ($\log_{10} B = 2.2$).

4. Aanbevelingen en Conclusies

De auteurs formuleren de volgende richtlijnen voor het kiezen van de startfrequentie ($f_{22}$) in GW-analyses:

1. Voor $\rho < 20$: Start bij 20 Hz. Statistische onzekerheden zijn groot genoeg om biases door het missen van multipoles te maskeren.
2. Voor $20 \lesssim \rho \lesssim 70$: Start bij 13 Hz (of lager). Dit is noodzakelijk om de $(3,3)$-multipool volledig te includeren en biases in massa en spin te voorkomen.
3. Voor $\rho \gtrsim 70$ en/of $M \gtrsim 300 M_\odot$: Start bij 10 Hz. Dit is vereist om ook de $(4,4)$-multipool volledig te vangen, vooral voor zware en asymmetrische systemen.

5. Significatie

Dit onderzoek heeft grote implicaties voor de interpretatie van huidige en toekomstige GW-observaties:

• Astrofysische conclusies: Systematische biases kunnen leiden tot verkeerde conclusies over de vorming van IMBH's en of deze zich in het "pair-instability mass gap" bevinden.
• Modelgebruik: Het benadrukt dat het simpelweg starten van tijdsdomein-modellen bij de standaard detectorfrequentie (20 Hz) onvoldoende is voor nauwkeurige inferentie bij hoge SNR's.
• Toekomstige detectoren: Aangezien toekomstige detectoren (zoals Einstein Telescope of Cosmic Explorer) een betere gevoeligheid bij lage frequenties hebben, zal het vermogen om hogere multipoles te detecteren toenemen, waardoor het correct hanteren van startfrequenties nog kritischer wordt.

Samenvattend waarschuwt de paper dat het negeren van het "missing multipole problem" leidt tot onbetrouwbare astrofysische parameters, en biedt het een kwantitatief kader om de startfrequentie van templates te optimaliseren op basis van de SNR en de massa van het bronstelsel.