Black Hole Spin-down in Collapsars in 3D Neutrino Transport GRMHD Simulations

1. समस्या विवरण (Problem Statement)

लंबे समय तक चलने वाले गामा-रे बर्स्ट्स (LGRBs) के बारे में माना जाता है कि वे "कोलैप्सार्स" (collapsars) द्वारा संचालित होते हैं—ये वे विशाल तारे हैं जिनके कोर ढहकर ब्लैक होल (BH) में बदल जाते हैं, जिससे एक अभिवृद्धि डिस्क (accretion disk) बनती है जो ब्लैंडफोर्ड-ज़नाजेक (Blandford-Znajek - BZ) तंत्र के माध्यम से सापेक्षतावादी जेट (relativistic jets) लॉन्च करती है। इन जेट की शक्ति ($L_{BZ} \propto a^2 \Phi^2$) ब्लैक होल के स्पिन ($a$) और चुंबकीय फ्लक्स ($\Phi$) पर महत्वपूर्ण रूप से निर्भर करती है।

गैर-विकिरणी (non-radiative/adiabatic) डिस्क के पिछले जनरल रिलेटिविस्टिक मैग्नेटोहाइड्रोडायनामिक (GRMHD) सिमुलेशनों ने सुझाव दिया था कि मैग्नेटिकली अरेस्टेड डिस्क (MADs) ब्लैक होल को बहुत कम संतुलन स्पिन ($a_{eq} \approx 0.035 - 0.07$) तक तेजी से स्पिन-डाउन कर देती हैं। इतना कम स्पिन जेट को जनक तारे (progenitor star) से बाहर निकलने के लिए बहुत कमजोर बना देगा, जो शक्तिशाली LGRBs के प्रेक्षणों (observations) के विपरीत है। इसके अलावा, बाइनरी ब्लैक होल विलय के गुरुत्वाकर्षण तरंग (GW) अवलोकनों से कम स्पिन वितरण का संकेत मिलता है, लेकिन वास्तविक, शीतलन वातावरण (cooling environments) में कोलैप्सार-जनित ब्लैक होल के विशिष्ट स्पिन विकास के बारे में जानकारी अज्ञात थी।

मुख्य समस्या जिसे संबोधित किया गया है वह यह है: न्यूट्रिनो कूलिंग (neutrino cooling), जो घने, गर्म कोलैप्सार डिस्क में महत्वपूर्ण है, ब्लैक होल के संतुलन स्पिन और जेट पावर को कैसे प्रभावित करती है? पिछले अध्ययनों में गणनात्मक लागत (computational cost) के कारण न्यूट्रिनो ट्रांसपोर्ट का अभाव था, जिससे डिस्क के ऊष्मागतिकी (thermodynamics) और कोणीय संवेग परिवहन (angular momentum transport) का गलत चित्रण होने की संभावना थी।

2. कार्यप्रणाली (Methodology)

लेखकों ने $\nu$H-AMR नामक GPU-त्वरित कोड का उपयोग करके 3D GRMHD सिमुलेशन का एक समूह निष्पादित किया। मुख्य कार्यप्रणाली संबंधी विशेषताएं निम्नलिखित हैं:

• न्यूट्रिनो ट्रांसपोर्ट: तीन न्यूट्रिनो प्रजातियों ($\nu_e, \bar{\nu}_e, \nu_x$) के लिए एक टू-मोमेंट, ऊर्जा-एकीकृत ("ग्रे") M1 योजना लागू की गई। यह न्यूट्रिनो कूलिंग का स्व-सुसंगत (self-consistent) मॉडलिंग करने की अनुमति देता है, जो डिस्क की मोटाई को कम करता है और ऊष्मागतिकी को बदल देता है।
• प्रोजेनेटर मॉडल (Progenitor Models): द्रव्यमान अभिवृद्धि दर ($\dot{m}$) की एक सीमा को कवर करने के लिए दो अलग-अलग स्टेलर डेंसिटी प्रोफाइल का अनुकरण किया गया:
  • p# मॉडल: पावर-लॉ डेंसिटी प्रोफाइल ($\rho \propto r^{-1.5}$), जो कोलैप्सार के विशिष्ट हैं, जिससे $\dot{m} \sim 0.1 - 1 M_\odot/s$ प्राप्त होता है।
  • c# मॉडल: निरंतर-घनत्व कोर (constant-density cores) और उसके बाद तीव्र प्रोफाइल, जिससे उच्च प्रारंभिक-समय की अभिवृद्धि दर ($\dot{m} \gtrsim 10 M_\odot/s$) प्राप्त होती है।
• प्रारंभिक स्थितियाँ (Initial Conditions): सिमुलेशन विभिन्न प्रारंभिक स्पिन ($a = 0.0$ से $0.8$) और MAD गठन के अनुकूल चुंबकीय क्षेत्र विन्यासों वाले ब्लैक होल्स के साथ शुरू हुए।
• विश्लेषण ढांचा (Analysis Framework): लेखकों ने स्पिन-अप पैरामीटर ($s$) की गणना की और एक अर्ध-विश्लेषणात्मक स्पिन विकास मॉडल प्राप्त करने के लिए फ्लक्स को हाइड्रोडायनामिक (HD), इलेक्ट्रोमैग्नेटिक (EM), और न्यूट्रिनो (NU) घटकों में विभाजित किया। उन्होंने अपने परिणामों को विशिष्ट LGRB इंजन जीवनकाल ($t_{eng} \sim 40$ s) पर विस्तारित (extrapolate) किया।

3. प्रमुख योगदान (Key Contributions)

• प्रथम न्यूट्रिनो-कूल्ड MAD सिमुलेशन: यह न्यूट्रिनो ट्रांसपोर्ट के साथ पूर्ण 3D GRMHD को मिलाकर कोलैप्सार में ब्लैक होल स्पिन विकास की जांच करने वाला पहला अध्ययन है।
• संशोधित संतुलन स्पिन (Revised Equilibrium Spin): यह अध्ययन स्थापित करता है कि न्यूट्रिनो कूलिंग संतुलन स्पिन को महत्वपूर्ण रूप से बदल देती है, जिससे $a_{eq} \approx 0.13$ का मान प्राप्त होता है। यह गैर-विकिरणी MAD सिमुलेशन ($a_{eq} \approx 0.035 - 0.07$) में पाए गए मानों की तुलना में 2–4 गुना अधिक है।
• स्पिन-डाउन तंत्र का परिमाणीकरण: लेखक मात्रात्मक रूप से प्रदर्शित करते हैं कि हालांकि न्यूट्रिनो फ्लक्स स्वयं ब्लैक होल पर नगण्य टॉर्क (torque) लगाते हैं, लेकिन न्यूट्रिनो कूलिंग अत्यंत महत्वपूर्ण है क्योंकि यह डिस्क को पतला करती है। पतली डिस्क चुंबकीय क्षेत्र द्वारा ब्लैक होल पर लगाए जाने वाले स्पिन-डाउन टॉर्क की दक्षता को कम कर देती है, जिससे ब्लैक होल उच्च स्पिन बनाए रखने में सक्षम होता है।
• जेट स्थिरता विश्लेषण: पेपर कम ब्लैक होल स्पिन को जेट अस्थिरता (kink instability) से जोड़ता है, यह समझाते हुए कि कम-स्पिन वाले मॉडल अस्थिर या विफल जेट क्यों उत्पन्न करते हैं, विशेष रूप से उथले घनत्व प्रोफाइल वाले प्रोजेनेटर्स में।

4. मुख्य परिणाम (Key Results)

• संतुलन स्पिन ($a_{eq}$):
  • MAD अवस्था में, ब्लैक होल स्पिन तेजी से एक ऐसे संतुलन की ओर विकसित होता है जहाँ स्पिन-अप पैरामीटर $s(a_{eq}) = 0$ होता है।
  • प्राप्त संतुलन स्पिन $a_{eq} \approx 0.13$ है।
  • यह मान LIGO/Virgo/Kagra द्वारा अनुमानित बाइनरी BH विलय के स्पिन वितरण (जो $a \lesssim 0.3$ के आसपास केंद्रित है) के अनुरूप है।
• जेट पावर में वृद्धि:
  • चूंकि जेट पावर $a^2$ के साथ स्केल करती है, इसलिए उच्च संतुलन स्पिन ($0.13$बनाम$0.07$) गैर-विकिरणी MAD मॉडलों की तुलना में 4–16 गुना अधिक शक्तिशाली जेट का परिणाम देता है।
  • ये उन्नत जेट जनक तारे से सफलतापूर्वक बाहर निकलने (drill out) में सक्षम हैं, जो पिछले कम-स्पिन मॉडलों की "विफल जेट" समस्या को हल करता है।
• स्पिन विकास समय-सीमा (Spin Evolution Timescales):
  • विशिष्ट LGRB अवधियों ($t \gtrsim 30$ s) के लिए, उच्च- और निम्न-अभिवृद्धि दर वाले दोनों प्रकार के प्रोजेनेटर्स में ब्लैक होल के पास $a_{eq} \approx 0.13$ तक पहुँचने के लिए पर्याप्त द्रव्यमान एकत्र करने का समय होता है।
  • छोटे या कम-अभिवृद्धि वाले इवेंट थोड़े उच्च स्पिन ($a > 0.13$) के साथ ब्लैक होल छोड़ सकते हैं।
• न्यूट्रिनो की भूमिका:
  • ब्लैक होल पर प्रत्यक्ष न्यूट्रिनो टॉर्क नगण्य हैं।
  • न्यूट्रिनो की प्राथमिक भूमिका कूलिंग (शीतलन) है, जो डिस्क के पहलू अनुपात (aspect ratio, $h/r$) को कम करती है। पतली डिस्क चुंबकीय टॉर्क के माध्यम से ब्लैक होल से कोणीय संवेग को दूर ले जाने में कम प्रभावी होती है, जिसके परिणामस्वरूप उच्च संतुलन स्पिन प्राप्त होता है।
• परिवर्तनशीलता और अस्थिरता (Variability and Instability):
  • कम प्रारंभिक स्पिन ($a=0.1$) और उथले घनत्व प्रोफाइल (p#) वाले मॉडल चुंबकीय फ्लक्स और जेट दक्षता में महत्वपूर्ण परिवर्तनशीलता प्रदर्शित करते हैं।
  • यह ग्लोबल किंक इंस्टेबिलिटी (global kink instability) के कारण होता है जो कमजोर जेट में देखी जाती है, जिससे डिस्क का मिसअलाइनमेंट (misalignment) और MAD अवस्था से अस्थायी निष्कासन होता है। यह तंत्र कम-ल्यूमिनोसिटी या विफल GRBs की व्याख्या कर सकता है।

5. महत्व (Significance)

• "लो स्पिन" विरोधाभास का समाधान: यह अध्ययन उच्च जेट पावर (जिसके लिए उच्च स्पिन आवश्यक है) और GW अवलोकनों के कम-स्पिन बाइनरी BH के बीच के सामंजस्य को सुलझाता है। यह सुझाव देता है कि कोलैप्सार ब्लैक होल स्वाभाविक रूप से $a \approx 0.13$ पर स्थिर होते हैं, जो GW डेटा से मेल खाने के लिए पर्याप्त कम है लेकिन सफल LGRBs को शक्ति प्रदान करने के लिए पर्याप्त उच्च है।
• GW प्रोजेनेटर्स के लिए निहितार्थ: यदि LGRBs बाइनरी BH विलय के प्रोजेनेटर हैं, तो अनुमानित स्पिन वितरण ($a \sim 0.1 - 0.2$) वर्तमान GW जनसंख्या अध्ययनों के साथ पूरी तरह से मेल खाता है, जो बाइनरी ब्लैक होल (BBH) विलय के एक हिस्से के लिए कोलैप्सार मॉडल का समर्थन करता है।
• GRB ल्यूमिनोसिटी पर ऊपरी सीमाएँ: इंजन गतिविधि के पहले कुछ सेकंड के भीतर उच्च-स्पिन ब्लैक होल्स ($a \sim 1$) का $a \approx 0.13$ तक तेजी से स्पिन-डाउन होना, देखे गए GRB ल्यूमिनोसिटी के ऊपरी कटऑफ के लिए एक भौतिक तंत्र प्रदान करता है। भले ही कोई तारा अधिकतम स्पिनिंग कोर के साथ ढह जाए, फिर भी जेट पावर को जेट के तारे से बाहर निकलने से पहले तेजी से स्पिन-डाउन द्वारा सीमित कर दिया जाता है।
• भविष्य का मॉडलिंग: यह कार्य कोलैप्सार सिमुलेशन में जेट ऊर्जा और r-प्रक्रिया न्यूक्लियोसिंथेसिस (r-process nucleosynthesis) स्थितियों की सटीक भविष्यवाणी करने के लिए न्यूट्रिनो ट्रांसपोर्ट को शामिल करने की आवश्यकता पर बल देता है।

निष्कर्षतः, यह शोध पत्र दर्शाता है कि न्यूट्रिनो-कूल्ड MADs को कोलैप्सार ब्लैक होल्स को पहले के अनुमानों की तुलना में उच्च संतुलन स्पिन ($a_{eq} \approx 0.13$) बनाए रखने की अनुमति मिलती है, जिससे सफल गामा-रे बर्स्ट्स का उत्पादन संभव होता है और साथ ही यह बाइनरी ब्लैक होल स्पिन के गुरुत्वाकर्षण तरंग अवलोकनों के अनुरूप भी रहता है।