From cosmological simulations to binary black hole mergers: The impact of using analytical star formation history models on gravitational-wave source populations

तकनीकी सारांश: कॉस्मोलॉजिकल सिमुलेशन से बाइनरी ब्लैक होल विलय तक

समस्या विवरण
बाइनरी ब्लैक होल (BBH) विलय के अवलोकन ब्रह्मांडीय समय के दौरान विशाल तारा विकास के बारे में महत्वपूर्ण अंतर्दृष्टि प्रदान करते हैं। हालाँकि, उनके पूर्वज तारों (progenitor stars) को जोड़ने के लिए सटीक मॉडलों की आवश्यकता होती है। वर्तमान में, यह क्षेत्र कॉस्मिक स्टार फॉर्मेशन हिस्ट्री (SFRH) और रासायनिक संवर्धन के लिए विश्लेणात्मक फिट (analytical fits) पर निर्भर करता है, जो जटिल अंतर्निहित इतिहास के प्रॉक्सी के रूप में कार्य करते हैं जिन्हें कॉस्मोलॉजिकल सिमुलेशन या अवलोकनों से प्राप्त किया गया है। ये विश्लेणात्मक प्रिस्क्रिप्शन अक्सर अनियंत्रित होते हैं और उनके प्रभाव BBH आबादी के अनुमानित गुणों पर अपर्याप्त रूप से समझे गए हैं। इसके अलावा, व्युत्पन्न $S(Z, z)$ और उसके बाद के BBH विलय दरों पर कॉस्मोलॉजिकल सिमुलेशन रिज़ॉल्यूशन और बॉक्स आकार के प्रभावों को व्यवस्थित रूप से मापा नहीं गया है।

कार्यप्रणाली
लेखक $S(Z, z)$ को सीधे गैस सेल डेटा से प्राप्त करने के लिए IllustrisTNG सूट (विशेष रूप से TNG50-1, TNG100-1, और TNG300-1) का उपयोग करते हैं, जो रिज़ॉल्यूशन और बॉक्स आकार में भिन्न होते हैं। वे अपने प्रत्यक्ष सिमुलेशन परिणामों की तुलना एक व्यापक रूप से उपयोग किए जाने वाले विश्लेणात्मक फिट मॉडल (L. A. C. van Son et al. 2023 पर आधारित) से करते हैं, जो कुल SFR घनत्व $S(z)$ के लिए एक Madau & Dickinson (2014) कार्यात्मक रूप को स्क्यूड लॉग-नॉर्मल मेटालिसिटी वितरण के साथ जोड़ता है।

BBH आबादी पर इन विभिन्न $S(Z, z)$ इनपुट के प्रभाव का आकलन करने के लिए, लेखक दोनों सिमुलेशन-व्युत्पन्न और विश्लेणात्मक $S(Z, z)$ मॉडलों को COMPAS (Compact Object Mergers: Population Astrophysics and Statistics) रैपिड बाइनरी पॉपुलेशन सिंथेसिस कोड के साथ जोड़ते हैं। वे रेडशिफ्ट के माध्यम से विलय दर और द्रव्यमान वितरण की गणना करने के लिए एक "कॉस्मिक इंटीग्रेशन" ढांचे का उपयोग करते हैं। अध्ययन विशेष रूप से निम्नलिखित की जांच करता है:

1. विभिन्न TNG रिज़ॉल्यूशन से $S(Z, z)$ के आकार और नॉर्मलाइजेशन को पुनरुत्पादित करने में विश्लेणात्मक फिट की निष्ठा (fidelity)।
2. पूर्ण सिमुलेशन डेटा बनाम विश्लेणात्मक फिट का उपयोग करने पर BBH विलय दर और प्राथमिक द्रव्यमान वितरण ($M_{BH,1}$) में होने वाले विचलन।
3. BBH परिणामों पर सिमुलेशन रिज़ॉल्यूशन (TNG50-1 के उच्च-रिज़ॉल्यूशन से लेकर TNG300-1 के निम्न-रिज़ॉल्यूशन तक) का प्रभाव।

मुख्य परिणाम
अध्ययन पाता है कि $S(Z, z)$ के लिए विश्लेणात्मक प्रिस्क्रिप्शन BBH आबादी के मॉडलिंग में महत्वपूर्ण व्यवस्थित त्रुटियाँ (systematic errors) उत्पन्न करते हैं:

• $S(Z, z)$ में विसंगतियाँ: विश्लेणात्मक फिट सिमुलेशन-आधारित मेटालिसिटी वितरण की प्रमुख संरचनात्मक विशेषताओं को पकड़ने में विफल रहता है। विशेष रूप से, यह उच्च-मेटालिसिटी शासन ($Z \gtrsim 3 Z_\odot$) में एक माध्यमिक "बम्प" (bump) और सिमुलेशन में मौजूद एक सपाट निम्न-मेटालिसिटी टेल (tail) को मिस करता है। फिट उच्च रेडशिफ्ट ($z \gtrsim 6$) और उच्च मेटालिसिटी पर स्टार फॉर्मेशन रेट को बढ़ा-चढ़ाकर दिखाता है, जबकि निम्न मेटालिसिटी और बहुत उच्च रेडशिफ्ट पर योगदान को कम करके दिखाता है, जो सिमुलेशन रिज़ॉल्यूशन पर निर्भर करता है।
• विलय दरों पर प्रभाव: उच्च रेडशिफ्ट ($z \gtrsim 6$) पर, विश्लेणात्मक फिट सिमुलेशन-आधारित दरों की तुलना में 10 से $10^4$ के कारकों तक BBH विलय दर को बढ़ाकर दिखाता है, जिसमें त्रुटि का परिमाण सिमुलेशन रिज़ॉल्यूशन घटने के साथ बढ़ता जाता है। निम्न रेडशिफ्ट ($z \approx 0$) पर, फिट से प्राप्त स्थानीय विलय दरें आम तौर पर सिमुलेशन-आधारित दरों के साथ $\sim 1.5$ के कारक के भीतर सहमत होती हैं, हालांकि TNG50-1 $\sim 15$ Gpc$^{-3}$ yr$^{-1}$ का ओवरएस्टिमेशन दिखाता है।
• द्रव्यमान वितरण पर प्रभाव: विश्लेणात्मक फिट का उपयोग BBH प्राथमिक द्रव्यमान वितरण में स्प्यूरियस (spurious) विशेषताएं पेश करता है। सबसे उल्लेखनीय रूप से, $z \lesssim 2$ पर, फिट प्राथमिक द्रव्यमान वितरण में $8 M_\odot$ के पास एक कृत्रिम अधिकता (bump) पैदा करता है, जो पूर्ण सिमुलेशन-आधारित $S(Z, z)$ का उपयोग करने पर अनुपस्थित होता है। इस विशेषता का श्रेय स्टेबल मास ट्रांसफर चैनल को दिया जाता है और यह उच्च-मेटालिसिटी स्टार फॉर्मेशन पीक को सही ढंग से मॉडल करने में फिट की अक्षमता के कारण होता है। इसके विपरीत, फिट $M_{BH,1} \gtrsim 35 M_\odot$ पर उच्च-द्रव्यमान विशेषताओं को 3–5 के कारकों से दबा देता है।
• रिज़ॉल्यूशन निर्भरता: सिमुलेशन रिज़ॉल्यूशन $S(z)$ के नॉर्मलाइजेशन और पीक लोकेशन को महत्वपूर्ण रूप से प्रभावित करता है। उच्च-रिज़ॉल्यूशन सिमुलेशन (TNG50-1) निम्न-रिज़ॉल्यूशन वाले रन ($z \approx 2.5$) की तुलना में उच्च स्टार फॉर्मेशन दरें प्रदान करते हैं और थोड़े कम रेडशिफ्ट ($z \approx 2$) पर पीक करते हैं। विश्लेणात्मक फिट की सटीकता निम्न सिमुलेशन रिज़ॉल्यूशन के साथ घटती जाती है, विशेष रूप से उच्च-रेडशिफ्ट ढलान (slope) को पुनरुत्पादित करने में।

महत्व और दावे
लेखक निष्कर्ष निकालते हैं कि जबकि विश्लेणात्मक फिट कम्प्यूटेशनल सुविधा प्रदान करते हैं, वे BBH आबादी की व्याख्या करते समय जटिल, अंतर्निहित $S(Z, z)$ के लिए अविश्वसनीय प्रॉक्सी हैं। अध्ययन प्रदर्शित करता है कि सरल विश्लेणात्मक मॉडल उन महत्वपूर्ण भौतिक विशेषताओं को पकड़ने में विफल रहते हैं—जैसे उच्च-मेटालिसिटी बम्प्स और विशिष्ट निम्न-मेटालिसिटी टेल्स—जो सीधे अनुमानित BBH विलय दरों और द्रव्यमान वितरण को प्रभावित करते हैं।

यह अध्ययन इस बात पर जोर देता है कि विसंगतियाँ इसलिए उत्पन्न होती हैं क्योंकि वर्तमान विश्लेणात्मक फिट के कार्यात्मक रूप सिमुलेशन में पाए जाने वाले मेटालिसिटी वितरण की मल्टी-मोडल प्रकृति को समायोजित नहीं कर सकते। फलस्वरूप, इन फिट्स पर निर्भर रहने से BBH के निर्माण चैनलों और विकासवादी इतिहास के बारे में गलत निष्कर्ष निकाले जा सकते हैं। लेखक तर्क देते हैं कि ग्रेविटेशनल-वेव डेटा की व्याख्या में स्प्यूरियस फीचर्स पेश करने से बचने के लिए कॉस्मोलॉजिकल सिमुलेशन, विश्लेणात्मक फिट और पॉपुलेशन सिंथेसिस के सावधानीपूर्वक एकीकरण की आवश्यकता है। वे यह दावा नहीं करते कि विश्लेणात्मक फिट पूरी तरह से बेकार हैं, बल्कि वे उनकी सीमाओं और उच्च-फिडेलिटी सिमुलेशन डेटा के विरुद्ध उन्हें मान्य करने की आवश्यकता को रेखांकित करते हैं, विशेष रूप से जैसे-जैसे अगली पीढ़ी के डिटेक्टर उच्च रेडशिफ्ट की जांच करेंगे।