Probing globular clusters parameters through gravitational wave lensing with stellar-mass black hole binaries

Technische Samenvatting: Het onderzoeken van parameters van bolvormige sterrenhopen door middel van gravitatiegolf-lenzing met binaire systemen van stellaire zwarte gaten

Probleemstelling
Bolvormige sterrenhopen (Globular Clusters, GC's) zijn dichte stellaire systemen die dienen als laboratoria voor het bestuderen van stellaire dynamica en vroege galactische formatie. Een cruciale parameter voor het begrijpen van hun interne dynamica, formatie en evolutie is hun massa, die doorgaans wordt afgeleid van de centrale stellaire snelheidsdispersie ($\sigma_v$). Huidige methoden voor het bepalen van de massa van GC's vertrouwen op elektromagnetische (EM) waarnemingen, zoals metingen van de snelheidsdispersie en oppervlaktehelderheidprofielen. Deze methoden kampen echter met aanzienlijke beperkingen, waaronder projectie-effecten langs de gezichtslijn, stellaire crowding in dichte centrale regio's, shot noise van heldere reuzensternen en onzekerheden bij het identificeren van het werkelijke centrum van de cluster. Deze uitdagingen worden verergerd in clusters met een kerninstorting (core-collapsed) of systemen die mogelijk een intermediair zwart gat (IMBH) huisvesten. Bov�elijk is het zo dat toekomstige breedveld-surveys (bijv. LSST, Euclid) de catalogi van GC's zullen uitbreiden met precieze posities en afstanden, maar het meten van interne kinematica zoals snelheidsdispersie blijft moeilijk en vereist doorgaans gerichte spectroscopische follow-up.

De auteurs stellen een complementair onderzoek voor: het gebruik van gravitatiegolven (GW's) van compacte binaire coalescenties (CBC's) die worden gelinst door GC's. In tegenstelling tot sterke lenzing door sterrenstelsels, die meerdere beelden produceert, werkt lenzing door GC's in het golfoptica-regime. In dit regime creëert de lens geen meerdere signalen, maar laat het een frequentieafhankelijke amplificatie achter op het GW-signaal. Het artikel onderzoekt of deze golfoptica-signaturen kunnen worden gebruikt om de intrinsieke eigenschappen van de lenseffectieve GC te herstellen, specifiek haar centrale snelheidsdispersie.

Methodologie
De studie maakt gebruik van een simulatie- en Bayesiaans parameter-schattingsframework om de herwinbaarheid van GC-parameters uit gelinst GW-signalen te testen.

1. Lensmodellering: GC's worden gemodelleerd als Singular Isothermal Spheres (SIS). De massaverdeling wordt benaderd door een dichtheidsprofiel $\rho(r) = \sigma_v^2 / (2\pi G r^2)$, waarbij de snelheidsdispersie $\sigma_v$ als constant wordt aangenomen. De effectieve lensingmassa ($M_{SIS}$) wordt afgeleid als een functie van $\sigma_v$ en de afstand tot de lens ($D_l$).
2. Signaal Simulatie: De auteurs simuleren GW-signalen die vergelijkbaar zijn met GW150914 (een binaire fusie van stellaire zwarte gaten). Deze signalen worden gelinst met de SIS-amplificatiefactor afgeleid van het Fresnel-Kirchhoff diffractie-integraal. De amplificatie is afhankelijk van de dimensieloze frequentie ($w$) en de impactparameter ($y$).
3. Injectiestudie: De gelinst signalen worden geïnjecteerd in Gaussische ruis gegenereerd uit de O4 power spectral densities (PSDs) van een netwerk van drie detectoren (LIGO Hanford, LIGO Livingston en Virgo).
4. Parameter-schatting: De geïnjecteerde signalen worden herwonnen met behulp van Bayesiaanse inferentie (via de Bilby-package en dynesty nested sampling) met waveform-templates die zowel bronparameters (15 binaire parameters met de IMRPhenomXPHM approximant) als lensparameters (SIS lensing massa $M_{SIS}$ en impactparameter $y$) bevatten.
5. Identificatie van Kandidaten: Vanuit de herwonnen posterior distributies wordt een 90% betrouwbaarheidsregio voor de hemellokalisatie geconstrueerd. Deze regio wordt gekruist met een catalogus van bekende Melkweg-GC's (uit [37]) om een potentiële lens te identificeren.
6. Herwinning van Snelheidsdispersie: Zodra een kandidaat GC is geïdentificeerd, wordt de bekende afstand ($D_l$) van deze cluster uit de catalogus gebruikt om de de werkelijke $\sigma_v$ te schatten door de herwonnen $M_{SIS}$ posterior te inverteren met behulp van de relatie $M_{SIS} \propto \sigma_v^4 D_l$.

Drie specifieke GC's uit de Melkweg-catalogus werden geselecteerd voor de studie op basis van hun onderscheidende centrale snelheidsdispersies: Liller-1 ($\sigma_v \approx 19,6$ km s$^{-1}$), NGC 6093 ($\sigma_v \approx 11,1$ km s$^{-1}$) en NGC 2419 ($\sigma_v \approx 5,1$ km s$^{-1}$). Voor elke cluster werden 20 injecties uitgevoerd met variërende impactparameters ($y \in [0,001, 5]$).

Belangrijkste Resultaten

• Herwinning van Lensingmassa en Snelheidsdispersie: Voor gunstige bron-lens uitlijningen (kleine impactparameters) werden de geïnjecteerde snelheidsdispersiewaarden succesvol herwonnen binnen de 68% betrouwbaarheidsintervallen. Voor Liller-1 was de mediaan herwonnen waarde 19,7 km s$^{-1}$, wat nauw overeenkomt met de geïnjecteerde 19,6 km s$^{-1}$.
• Impact van de Impactparameter: De nauwkeurigheid van de herwinning is sterk afhankelijk van de impactparameter $y$. Naarmate $y$ toeneemt, neemt de amplificatie af, waardoor de lensing-afdruk op het signaal verzwakt. Dit leidt tot een slechtere bemonstering van de lensingmassa en een grotere onzekerheid in de afgeleide snelheidsdispersie. Bij grotere $y$ kunnen de herwonnen mediane waarden verschuiven van de werkelijke waarde door ruisrealisaties.
• Hemellokalisatie en Catalogus-kruisverwijzing: De geïnjecteerde cluster werd correct geïdentificeerd binnen de 90% betrouwbaarheidsregio voor de hemellokalisatie in alle geteste gevallen. De methode slaagde erin de lensing-cluster te onderscheiden van andere in de catalogus, mits de bron zich achter de specifieke doelcluster bevond.
• Parameterdegeneratie: De studie observeerde een anti-correlatie tussen de herwonnen lensingmassa ($M_{SIS}$) en de impactparameter ($y$) voor $y \sim 1$, wat consistent is met eerdere bevindingen voor puntmassa-lenzen.

Betekenis en Claims
Het artikel demonstreert dat gelinst GW's kunnen dienen als een complementaire sonde voor de dynamica van GC's, specif gezien de mogelijkheid om de snelheidsdispersie en daarmee de massa van de cluster te schatten. Dit biedt een nieuwe methode om de beperkingen te overwinnen die inherent zijn aan EM-gebaseerde massabepaling, met name in dichte of met kerninstorting te maken hebbende regio's.

De auteurs merken op dat de waarschijnlijkheid dat een Melkweg-GC als lens fungeert voor een GW-gebeurtenis laag is ($\tau \sim 10^{-8}$), wat de potentie van volgende generaties detectoren om dergelijke gebeurtenissen te ontdekken benadrukt. De methodologie wordt gepresenteerd als een instrument om:

1. GW-gebeurtenissen met vertrouwen te identificeren als gelinst signalen, indien een Melkweg-GC de lens is.
2. De snelheidsdispersie en massa van de lenseffectieve cluster te beperken.
3. Potentieel de aanwezigheid van IMBH's in de kernen van clusters aan te geven.

De studie is bescheiden in haar reikwijdte, maakt gebruik van het SIS-model en richt zich op Melkweg-GC's waar de afstanden goed bekend zijn. De auteurs erkennen dat er voor extragalactische GC's (waar catalogi incompleet kunnen zijn) een degeneratie bestaat tussen de afstand van de lens en de snelheidsdispersie, wat uitgebreidere parameter-schatting vereist. Ze merken ook op dat hoewel het huidige werk het SIS-model gebruikt, de methodologie in de toekomst kan worden uitgebreid naar meer realistische modellen (bijv. cored isothermal spheres, King-modellen). Het artikel concludeert dat een globaal netwerk van ten minste drie GW-detectoren met vergelijkbare gevoeligheid essentieel is om de vereiste nauwkeurigheid van de hemellokalisatie te verbeteren zodat deze techniek effectief kan zijn.