Reversible-jump MCMC reveals binary black hole subpopulations with distinct redshift evolution

तकनीकी सारांश: रिवर्सिबल-जंप MCMC ने विशिष्ट रेडशिफ्ट विकास के साथ द्विआधारी ब्लैक होल उप-जनसंख्याओं का अनावरण किया

समस्या विवरण
बाइनरी ब्लैक होल (BBH) के निर्माण तंत्र अभी भी खगोल भौतिकी में एक केंद्रीय खुला प्रश्न बना हुआ है। विभिन्न चैनल—जैसे कि पृथक बाइनरी विकास (isolated binary evolution), घने क्लस्टरों में गतिशील संयोजन (dynamical assembly), रासायनिक रूप से सजातीय विकास (CHE), और सक्रिय गैलेक्टिक न्यूक्लियस (AGN) डिस्क में पदानुक्रमित विलय (hierarchical mergers)—द्रव्यमान, स्पिन और रेडशिफ्ट के बहुआयामी पैरामीटर स्पेस में विशिष्ट संकेतों को उत्पन्न करने की भविष्यवाणी करते हैं। जबकि लाइगो-विर्गो-कैगरा (LVK) सहयोग ने GWTC-4 कैटलॉग जारी किया है, जिससे देखे गए आयोजनों की संख्या दोगुनी हो गई है, उप-जनसंख्याओं की पहचान करने वाली मौजूदा विधियों को एक समझौते (trade-off) का सामना करना पड़ता है। दृढ़ता से-मॉडल किए गए दृष्टिकोण (जैसे कि विशिष्ट पैरामेट्रिक रूपों का मिश्रण) व्याख्यात्मकता प्रदान करते हैं लेकिन डेटा में अनुपस्थित संरचनाओं को थोपने या अप्रत्याशित विशेषताओं को चूक जाने का जोखिम उठाते है। इसके विपरीत, दुर्बल रूप से-मॉडल किए गए दृष्टिकोण (जैसे कि गैर-पैरामीट्रिक स्मूथिंग) डेटा को आकार निर्धारित करने देते हैं लेकिन अक्सर खगोल भौतिकीय व्याख्यात्मकता की कमी रखते हैं और उन्हें उप-जनसंख्याओं की पहचान करने के लिए अनुमानी (heuristic) पोस्ट-प्रोसेसिंग की आवश्यकता होती है।

कार्यप्रणाली
लेखक BBH उप-जनसंख्याओं की खोज के लिए रिवर्सिबल-जंप मार्कोव चेन मोंटे कार्लो (RJ MCMC) का उपयोग करते हुए एक नवीन ढांचा प्रस्तुत करते हैं। यह विधि पैरामीटर अनुमान के साथ-साथ विभिन्न जटिलताओं वाले मॉडलों (विभिन्न उप-जनसंख्याओं की संख्या, $n$) के बीच बेयसियन मॉडल तुलना करती है।

• ढांचा: BBH जनसंख्या को 4-आयामी पैरामीटर स्पेस $\theta = [m_1, q, \chi_{\text{eff}}, z]$ पर $n$ घटकों के मिश्रण के रूप में मॉडल किया गया है। विभेदक विलय दर स्थानीय दरों और प्रत्येक उप-जनसंख्या के वितरण का योग है।
• मॉडल वेरिएंट्स: मजबूती सुनिश्चित करने के लिए, लेखक कई जनसंख्या मॉडलों का परीक्षण करते हैं:
  1. skewt: प्राथमिक द्रव्यमान वितरण को जोन्स और फैडी स्क्यू-टी (skew-t) वितरण के रूप में मॉडल किया गया है (जो लचीला है और पावर-लॉ एवं गॉसियन व्यवहारों को कैप्चर करता है)।
  2. NPLNP: पावर-लॉ और गॉसियनों का एक मिश्रण, जो मानक घटनात्मक (phenomenological) मॉडलों के तुलनीय है, जहाँ $n$ एक स्वतंत्र पैरामीटर है।
  3. skewt ID: एक ऐसा मॉडल जहाँ दोनों ब्लैक होल एक ही अंतर्निहित द्रव्यमान वितरण से लिए जाते हैं और एक सहसंबंध गुणांक $\rho$ द्वारा जुड़े होते हैं, जो द्रव्यमान युग्मन (mass pairing) तंत्र की जांच करता है।
  4. skewt z: एक वेरिएंट जो प्रत्येक उप-जनसंख्या के लिए स्वतंत्र रेडशिफ्ट विकास ($\gamma_k$) की अनुमति देता है।
• कार्यान्वयन: यह विश्लेषण GWTC-4 से 153 BBH आयोजनों पर लागू eryn पैकेज का उपयोग करके RJ MCMC का उपयोग करता है। उप-जनसंख्या दरों में डैजेनेरेसी (degeneracy) को संभालने के लिए कस्टम गिब्स प्रपोजल स्टेप्स लागू किए गए हैं। यह विधि 'ऑकेम के रेज़र' (Occam's razor) के माध्यम से ओवरफिटिंग को स्वतः दंडित करती है, और सरल मॉडलों को प्राथमिकता देती है जब तक कि अतिरिक्त घटक महत्वपूर्ण रूप से लाइकलीहुड (likelihood) में सुधार न करें।
• पोस्ट-प्रोसेसिंग: विभिन्न पोस्टीरियर सैंपल्स के बीच घटकों को क्लस्टर करने के लिए एक "व्हाइटनिंग" फीचर स्पेस मैपिंग का उपयोग किया जाता है, जिससे यह सुनिश्चित होता है कि विशिष्ट उप-जनसंख्याओं की पहचान विशिष्ट मॉडल रियलाइजेशन के बावजूद सुसंगत रूप से की जा सके।

मुख्य योगदान और परिणाम
GWTC-4 डेटा पर इस ढांचे को लागू करते हुए, लेखक तीन अलग-अलग उप-जनसंख्याओं की पहचान मजबूती से करते हैं जो विभिन्न जनसंख्या मॉडल वेरिएंट्स में सुसंगत हैं:

1. $10\,M_\odot$ पीक:

   • विशेषताएं: $m_1 \approx 10\,M_\odot$ पर केंद्रित एक संकीर्ण उप-जनसंख्या, जिसमें एक छोटा, सकारात्मक औसत प्रभावी स्पिन ($\chi_{\text{eff}} \approx 0.05$) और असमान द्रव्यमान अनुपात ($q < 1$) के प्रति प्राथमिकता है।
   • व्याख्या: यह पृथक क्षेत्र बाइनरी विकास (isolated field binary evolution), विशेष रूप से स्थिर द्रव्यमान हस्तांतरण (SMT) परिदृश्यों के अनुरूप है जहाँ टाइडल इंटरैक्शन स्पिन को संरेखित करते हैं और द्रव्यमान हस्तांतरण असमान युग्मों को प्राथमिकता देता है।
   • रेडशिफ्ट विकास: यह उप-जनसंख्या $30\,M_\odot$ पीक की तुलना में काफी तेज़ रेडशिफ्ट विकास ($\gamma \approx 5.7$) प्रदर्शित करती है। यह एक लघु विलंब-समय वितरण (delay-time distribution) और निम्न मेटालिसिटी कटऑफ ($Z_{\text{max}} \lesssim 0.1\,Z_\odot$) के ऊपर अक्षम BBH निर्माण का संकेत देता है।

2. $30\,M_\odot$ पीक:

   • विशेषताएं: $m_1 \sim 30\,M_\odot$ के आसपास केंद्रित एक विस्तृत वितरण, जिसमें $\chi_{\text{eff}}$ वितरण शून्य पर केंद्रित है और समान द्रव्यमान अनुपात ($q \approx 1$) के प्रति मजबूत प्राथमिकता है।
   • व्याख्या: यह घने क्लस्टरों में गतिशील संयोजन (dynamical assembly) के अनुरूप है, जहाँ द्रव्यमान पृथक्करण (mass segregation) और आइसोट्रोपिक स्पिन ओरिएंटेशन अपेक्षित होते हैं।

3. हाई-स्पिन कॉन्टिनम:

   • विशेषताएं: एक निरंतरता (continuum) में फैले द्रव्यमानों (उच्च-द्रव्यमान पूंछ सहित) वाला एक उप-जनसंख्या, जिसमें एक विस्तृत, सकारात्मक-औसत $\chi_{\text{eff}}$ वितरण ($\mu_{\chi} \approx 0.23$) है। इसमें कैटलॉग के कई असाधारण आयोजन शामिल हैं (जैसे, GW190412, GW231123)।
   • व्याख्या: लेखक इसे एक "कैच-ऑल" घटक के रूप में व्याख्यायित करते हैं जो कई उप-प्रमुख चैनलों के संगम का प्रतिनिधित्व करता है। सकारात्मक औसत स्पिन, पदानुक्रमित विलय के लिए शुद्ध क्लस्टर-आधारित अपेक्षाओं (जो शून्य के चारों ओर सममित $\chi_{\text{eff}}$ की भविष्यवाणी करते हैं) का खंडन करता है। इसके बजाय, यह CHE या AGN-डिस्क पदानुक्रमित विलय के योगदान का सुझाव देता है। डेटा इस विशेषता के लिए 2-घटक मॉडल की तुलना में 3-घटक मॉडल को मामूली रूप से प्राथमिकता देता है, हालांकि साक्ष्य अत्यधिक नहीं है।

महत्व
यह शोध दावा करता है कि यह एक नवीन, डेटा-संचालित ढांचे की नींव रखता है जो दृढ़ रूप से-मॉडल किए गए और दुर्बल रूप से-मॉडल किए गए दृष्टिकोणों के बीच के अंतर को पाटता है। RJ MCMC का उपयोग करके, लेखक यह प्रदर्शित करते हैं कि उप-जनसंख्याओं की संख्या या उनके विशिष्ट कार्यात्मक रूपों पर कड़े पूर्व अनुमान (prior assumptions) लगाए बिना, BBH जनसंख्या में मजबूत, खगोल भौतिकीय रूप से व्याख्यात्मक विशेषताओं को हल करना संभव है।

• अज्ञेयवादी खोज (Agnostic Discovery): ढांचा बिना हाई-स्पिन कॉन्टिनम के $\chi_{\text{eff}}$ वितरण पर समरूपता लागू किए, तीन उप-जनसंख्याओं को सफलतापूर्वक रिकवर करता है, जो कि पदानुक्रमित विलय के लिए लक्षित खोजों में मौजूद एक बाधा थी।
• रेडशिफ्ट विकास: यह कार्य पहला डेटा-संचालित प्रमाण प्रदान करता है कि $10\,M_\odot$ पीक, $30\,M_\odot$ पीक की तुलना में रेडशिफ्ट के साथ अलग तरह से विकसित होती है, जो विलंब-समय वितरण और मेटालिसिटी निर्भरता पर नए प्रतिबंध प्रदान करता है।
• भविष्य की दृष्टि: लेखक नोट करते हैं कि जबकि वर्तमान विश्लेषण अधिकांश आयोजनों को विशिष्ट उप-जनसंख्याओं को आवंटित करता है, "हाई-स्पिन कॉन्टिनम" संभावित रूप से अलग-अलग चैनलों का एक मिश्रण बना हुआ है। वे तर्क देते हैं कि $\mathcal{O}(10^3)$ आयोजनों और बेहतर स्पिन मापों ( $\chi_{\text{eff}}$ से परे) के साथ भविष्य के कैटलॉग इस कॉन्टिनम को इसके घटक निर्माण चैनलों में विभाजित करने में सक्षम होंगे, जिससे BBH मूलों का प्रश्न एक अनुभवजन्य (empirical) प्रश्न में बदल जाएगा।