Forbidden Formation Histories: The Binary Black Hole Merger Rate Disfavors Long Delay Times

Technische Samenvatting: Verboden vormingsgeschiedenissen: De samensmeltingsfrequentie van dubbele zwarte gaten sluit lange vertragingstijden uit

Probleemstelling
De samensmeltingsfrequentie van dubbele zwarte gaten (BBH), $R_{\text{merge}}(t)$, is de convolutie van de vormingsfrequentie van de voorlopers, $R_{\text{form}}(t)$, en de verdeling van de vertragingstijden tussen vorming en samensmelting, $p(\tau)$. Traditionele benaderingen om de oorsprong van BBH's te begrijpen, vertrouwen vaak op forward modeling: het aannemen van specifieke functionele vormen voor $R_{\text{form}}$ (bijvoorbeeld gekoppeld aan de globale stervormingsfrequentie) en $p(\tau)$ (afgeleid van populatiesynthese), en vervolgens het vergelijken van de voorspelde samensmeltingsfrequentie met data van zwaartekrachtgolven (GW). Deze aanpak lijdt aan inherente degeneraties; verschillende combinaties van vormingsgeschiedenissen en vertragingstijden kunnen identieke samensmeltingsfrequenties produceren, en rigide parametrische aannames kunnen datakenmerken verbergen of fysiek inconsistente resultaten afdwingen. Bovendien kunnen sterke modellen voor $R_{\text{form}}$ en $p(\tau)$ fundamenteel in strijd zijn met de waargenomen evolutie van de samensmeltingsfrequentie met de roodverschuiving.

Methodologie
De auteurs stellen een "backward-modeling"-kader voor om de waargenomen BBH-samensmeltingsfrequentie te deconvolueren in zijn samenstellende componenten, zonder een specifieke functionele vorm voor de vormingsfrequentie van de voorlopers aan te nemen.

1. Wiskundig Kader: Met gebruikmaking van de convolutiestelling transformeren de auteurs de integraalvergelijking $R_{\text{merge}}(t) = \int R_{\text{form}}(t-\tau)p(\tau)d\tau$ naar het frequentiedomein via de Fast Fourier Transform (FFT). Dit maakt directe reconstructie van $R_{\text{form}}$ mogelijk, gegeven een waargenomen $R_{\text{merge}}$ en een bekende $p(\tau)$:
   $$R_{\text{form}} = \mathcal{F}^{-1} \left\{ \frac{\mathcal{F}\{R_{\text{merge}}\}}{\mathcal{F}\{p(\tau)\}} \right\}$$
2. Numerieke Implementatie: Om om te gaan met het niet-periodieke karakter van de data en de eindige observatiehorizon (GWTC-4.0-data tot $z \approx 1.5$), gebruiken de auteurs zero-padding om wrap-around-artefacten te voorkomen en passen ze cosinustapering toe om Gibbs-fenomeenruis aan de randen te mitigeren. Ze veronderstellen een conservatieve randvoorwaarde waarbij $R_{\text{merge}}(z > 1.5) = 0$, hoewel ze asymptotische waarden tot 500 Gpc$^{-3}$yr$^{-1}$ testen.
3. Regularisatie: Om numerieke stabiliteit te waarborgen bij het delen door de Fourier-transformatie van de vertragingstijdverdeling (die nul kan benaderen), maken ze gebruik van Wiener-deconvolutie, waarbij een regularisatieterm $\lambda$ aan de noemer wordt toegevoegd.
4. Data-invoer: De methode wordt toegepast met behulp van de niet-parametrische B-spline-aanpassing aan de BBH-samensmeltingsfrequentie uit GWTC-4.0 en een reeks vertragingstijdverdelingen gegenereerd door de COMPAS-populatiesynthesecode, die parameters variëren gerelateerd aan massatransfer, supernova's en restmassa's.

Belangrijkste bijdragen en resultaten

• Reconstructie van vormingsgeschiedenissen: Door GWTC-4.0-samensmeltingsfrequenties te deconvolueren met COMPAS-vertragingstijdverdelingen, reconstrueren de auteurs de geïmpliceerde vormingsfrequentie van de voorlopers. Ze constateren dat de afgeleide $R_{\text{form}}$ steiler afneemt naar lage roodverschuivingen dan de globale stervormingsfrequentie (SFR), wat consistent is met een metaliciteitsafhankelijke vormingsefficiëntie van BBH's.
• Identificatie van "verboden" geschiedenissen: De studie toont aan dat bepaalde combinaties van samensmeltingsfrequenties en vertragingstijdverdelingen fysiek onmogelijk zijn, omdat ze vereisen dat $R_{\text{form}}$ op latere tijden negatief wordt. Dit treedt op wanneer een vertragingstijdverdeling met een lange staart (een ondiepe machtsverdelingsindex) meer samensmeltingen voorspelt uit oude vorming bij hoge roodverschuiving dan daadwerkelijk wordt waargenomen bij lage roodverschuiving.
• Beperkingen aan vertragingstijdverdelingen: Met gebruikmaking van de eis dat $R_{\text{form}} > 0$, beperken de auteurs de machtsverdelingsindex $\alpha$ van de vertragingstijdverdeling (waarbij $p(\tau) \propto \tau^\alpha$).
  • Vertragingstijdverdelingen met ondiepe indices ($\alpha \gtrsim -0.7$), kenmerkend voor kanalen met stabiele massatransfer, worden sterk afgekeurd.
  • Specifiek, voor de GWTC-4.0-data, wordt $\alpha < -1.08$ 19 keer meer geprefereerd dan $\alpha > -0.67$.
  • COMPAS-simulaties met steile verdelingen (die korte vertragingen bevorderen, $\alpha < -1.2$) zijn consistent met bijna alle posterior-trekkingen van de samensmeltingsfrequentie van GWTC-4.0, terwijl ondiepe verdelingen compatibel zijn met slechts ongeveer 30% van de trekkingen.
• Beperkingen aan de fysica van binaire evolutie: Als bewijs van concept beperken de auteurs de parameter voor massatransfer-efficiëntie $\beta$ in COMPAS. Lagere $\beta$-waarden (niet-conservatieve massatransfer, leidend tot orbitale inkrimping en korte vertragingen) tonen een hogere compatibiliteit met GW-data. Hoge $\beta$-waarden (conservatieve massatransfer, lange vertragingen) worden met een factor van ongeveer 3.7 afgekeurd.
• Impact op inferentie van samensmeltingsfrequentie: Het opleggen van de fysieke beperking van een positieve vormingsfrequentie als prior op het GWTC-4.0 B-spline-model verlaagt de mediaan afgeleide samensmeltingsfrequentie bij $z=1.5$ met ongeveer 10% (van 334 naar 299 Gpc$^{-3}$yr$^{-1}$), wat een vlakkere evolutie van de samensmeltingsfrequentie bevoordeelt dan de standaard datagedreven inferentie.
• Mapping van de parameter ruimte: De auteurs mappen het verboden gebied in de $\alpha$–$\kappa$-parameter ruimte (waarbij $\kappa$ de evolutie van de samensmeltingsfrequentie beschrijft, $R_{\text{merge}} \propto (1+z)^\kappa$). Ze tonen aan dat voor een gegeven $\kappa$ er een kritieke $\alpha_{\text{crit}}$ bestaat, waarboven de vormingsfrequentie onfysisch wordt. Huidige data sluit machtsverdelingsvertragingstijdverdelingen uit met $\alpha \gtrsim -0.7$.

Betekenis
Het artikel presenteert een modelonafhankelijk kader om gebieden van binaire evolutie en parameter ruimten van samensmeltingsfrequenties uit te sluiten. Door direct de compatibiliteit van parameters van populatiesynthese met waarnemingen van zwaartekrachtgolven te evalueren door de lens van causaliteit (d.w.z. dat de vormingsfrequentie niet negatief kan zijn), omzeilt de methode de noodzaak van hoge-dimensionale priors voor stervorming en evolutie van metaliciteit. De auteurs betogen dat deze aanpak "verboden" vormingsgeschiedenissen blootlegt die forward modeling zou kunnen missen, en zo een robuuste consistentiecheck biedt voor theorieën over de evolutie van sterrenbinaire systemen. De resultaten suggereren dat vormingskanalen die lange vertragingstijden voorspellen (zoals stabiele massatransfer) in spanning staan met de waargenomen steile afname van de BBH-samensmeltingsfrequentie bij lage roodverschuivingen, en dat kanalen die kortere vertragingen produceren worden bevoordeeld.