Numerical simulations of black hole-neutron star mergers with equal and near-equal mass ratios

तकनीकी सारांश: समान और लगभग समान द्रव्यमान अनुपात वाले ब्लैक होल-न्यूट्रॉन स्टार विलय के संख्यात्मक सिमुलेशन (Numerical Simulations)

समस्या और प्रेरणा
गुरुत्वाकर्षण तरंग घटना GW230529 का पता लगाना, जिसे द्रव्यमान अनुपात $q \approx 1/2.6$ के साथ एक ब्लैक होल-न्यूट्रॉन स्टार (BHNS) विलय के रूप में व्याख्यायित किया गया है, पहले की अपेक्षा अधिक सममित (symmetric) BHNS प्रणालियों के अस्तित्व का सुझाव देता है। जबकि मानक तारकीय विकास मॉडल द्रव्यमान अनुपात $q \lesssim 1/3.5$ की भविष्यवाणी करते हैं, सघन समूहों (dense clusters) में गतिशील निर्माण चैनल या दूसरी पीढ़ी के विलय $q \sim 1/2$ या यहाँ तक कि $q \approx 1$ वाले सिस्टम बना सकते हैं। ये सममित प्रणालियाँ सैद्धांतिक रूप से विशाल अभिवृद्धि डिस्क (accretion disks) उत्पन्न करने की उम्मीद करती हैं, जो चमकीले विद्युत चुम्बकीय क्षणिक (electromagnetic transients) जैसे कि किलोनोवा (kilonovae) को संचालित कर सकती हैं। हालाँकि, समान और लगभग समान द्रव्यमान अनुपातों के लिए पैरामीटर स्पेस संख्यात्मक-सापेक्षता (numerical-relativity - NR) सिमुलेशन द्वारा विरल रूप से भरा गया है। डेटा की यह कमी गुरुत्वाकर्षण तरंग मॉडल और अवशेष गुणों (remnant properties) के विश्लेषणात्मक फिटिंग सूत्रों की सटीकता को सीमित करती है, जिससे वास्तविक घटनाओं का गलत वर्गीकरण और विद्युत चुम्बकीय समकक्षों (counterparts) को पहचानने के अवसरों को खोने का जोखिम बढ़ जाता है।

कार्यप्रणाली (Methodology)
लेखक BAM कोड का उपयोग करके गैर-स्पिनिंग BHNS विलय के बारह NR सिमुलेशन करते हैं। सिमुलेशन द्रव्यमान अनुपात $q \in \{1, 1/2, 1/3\}$ और न्यूट्रॉन स्टार (NS) द्रव्यमान $m_{NS} \in \{0.8, 1.2, 1.4, 2.2\} M_\odot$ को कवर करते हैं। दो समीकरण अवस्थाओं (Equations of State - EOS) का उपयोग किया गया है: नरम SLy और कठोर DD2। सिमुलेशन स्पेसटाइम के लिए कॉन्फ़ॉर्मल Z4 फॉर्मूलेशन और हाइड्रोडायनामिक्स के लिए वैलेंसिया (Valencia) फॉर्मूलेशन का उपयोग करते हुए लगभग नौ कक्षाओं (orbits) के लिए विकसित होते हैं। प्रारंभिक डेटा FUKA कोड के साथ निर्मित किया गया है, जिसमें सनकीवद्धता (eccentricity) को $10^{-3}$ से नीचे कम किया गया है।

अध्ययन में शामिल हैं:

1. गुरुत्वाकर्षण तरंग विश्लेषण (Gravitational Waveform Analysis): कई त्रिज्याओं पर तरंगों (waveforms) को निकालना, द्वितीय-क्रम एक्सट्रपलेशन (second-order extrapolation) करना, और केंद्र-द्रव्यमान ड्रिफ्ट (center-of-mass drift) के लिए सुधार करना। इनकी तुलना IMRPhenomNSBH और SEOBNRv4_ROM_NRTidalv2_NSBH तरंग मॉडलों से की जाती है।
2. अवशेष गुण सत्यापन (Remnant Property Verification): गतिशील इजेक्टा (dynamical ejecta) के द्रव्यमान और अभिवृद्धि डिस्क के द्रव्यमान, तथा विलय के 10 ms बाद अवशेष ब्लैक होल (BH) के द्रव्यमान और स्पिन की गणना करना। इनकी तुलना संदर्भ [28, 31, 40, 41] से प्राप्त विश्लेषणात्मक फिटिंग सूत्रों से की जाती है।
3. डिस्क गतिशीलता (Disk Dynamics): महत्वपूर्ण डिस्क द्रव्यमान वाले विन्यासों के लिए अवशेष डिस्क की आकृति विज्ञान (morphology), रोटेशन नियम और वैश्विक घनत्व दोलनों (global density oscillations) का विश्लेषण करना।
4. अभिवृक्क मॉडलिंग (Accretion Modeling): अभिवृक्क व्यवहारों को वर्गीकृत करना और डिस्क दोलनों द्वारा होने वाले मॉड्यूलेशन की पहचान करना।
5. किलोनोवा मॉडलिंग (Kilonova Modeling): गतिशील इजेक्टा प्रोफाइल और एक अतिरिक्त गोलाकार विंड घटक (डिस्क द्रव्यमान का 30%) के आधार पर POSSIS रेडिएटिव ट्रांसफर कोड का उपयोग करके किलोनोवा प्रकाश वक्रों (lightcurves) को मॉडल करना।

मुख्य परिणाम

• गुरुत्वाकर्षण तरंगें (Gravitational Waveforms): वर्तमान BHNS तरंग मॉडल (IMRPhenomNSBH और SEOBNRv4) सममित द्रव्यमान अनुपात शासन में NR तरंगों के साथ खराब सामंजस्य दिखाते हैं। अधिकांश मामलों में विलय के समय लगभग 1 रेडियन का डीफेजिंग (dephasing) देखा जाता है। मॉडल अक्सर NR डेटा की तुलना में पहले विलय समय और विस्तारित रिंगडाउन चरणों की भविष्यवाणी करते हैं। विशेष रूप से, सब-सोलर-मास NS मामले ($m_{NS}=0.8 M_\odot$) के लिए, टाइडल डिफॉर्मेबिलिटी मौजूदा मॉडलों के वैध पैरामीटर स्पेस से बाहर है; एक्सट्रपलेशन करने पर, मॉडल महत्वपूर्ण चरण (phase) और आयाम विसंगतियां दिखाते हैं।
• अवशेष गुण (Remnant Properties):
  • इजेक्टा और डिस्क द्रव्यमान: गतिशील इजेक्टा द्रव्यमान और डिस्क द्रव्यमान के लिए विश्लेषणात्मक फिटिंग सूत्र सिमुलेशन डेटा के साथ अच्छा सामंजस्य दिखाते हैं।
  • अवशेष BH द्रव्यमान: विश्लेषणात्मक सूत्र अवशेष BH द्रव्यमान की अत्यधिक सटीकता के साथ भविष्यवाणी करते हैं (त्रुटियां आम तौर पर $< 1\%$ हैं)।
  • अवशेष BH स्पिन: अपडेट किए गए फिटिंग सूत्र अवशेष स्पिन की भविष्यवाणी करने के लिए कम सटीक हैं, विशेष रूप से उच्च टाइडल डिफॉर्मेबिलिटी (DD2 EOS) वाले सिस्टम के लिए। सब-सोलर NS मामले के लिए, स्पिन भविष्यवाणी लगभग 50% विचलित होती है। लेखक सुझाव देते हैं कि यह कम टाइडल डिफॉर्मेबिलिटी और अधिक विषम द्रव्यमान अनुपात वाले सिस्टम पर ओवरफिटिंग के कारण हो सकता है।
• डिस्क संरचना और गतिशीलता (Disk Structure and Dynamics):
  • आकृति विज्ञान (Morphology): समान-द्रव्यमान ($q=1$) डिस्क लगभग तुरंत ($\sim 3$ ms विलय के बाद) गोलाकार हो जाती है, जबकि $q=1/2$ डिस्क स्पाइरल डेंसिटी वेव्स और फॉलबैक मैटर से लंबे समय तक विक्षुब्ध रहती है, और $\sim 8$ ms के बाद ही अक्षीय समरूपता (axisymmetry) में स्थिर होती है।
  • रोटेशन नियम (Rotation Law): विशिष्ट कोणीय संवेग (specific angular momentum) एक पावर लॉ $j \propto \Omega^{-n}$ का पालन करता है। घातांक $n$ रेडियल रूप से भिन्न होता है, जो ISCO के अंदर स्थिर कोणीय संवेग ($n=0$) से बाहरी डिस्क में उप-केपलरियन (sub-Keplerian) मानों में परिवर्तित होता है। यह दीर्घकालिक अक्षीय सममित सिमुलेशन में सामान्य निश्चित रोटेशन नियमों की धारणाओं का खंडन करता है।
  • दोलन (Oscillations): $10^{-2} M_\odot$ से अधिक द्रव्यमान वाली डिस्क में वैश्विक ट्रैप्ड $g$-mode घनत्व दोलनों की पहचान की गई है। उनकी आवृत्तियाँ अवशेष BH की अधिकतम प्रोग्रेड एपिसाइकिल फ्रीक्वेंसी के अनुरूप हैं।
  • अभिवृक्क (Accretion): तीन प्रकार के अभिवृक्क को अलग किया गया है। टाइप B ($q=1/2$) वैश्विक $g$-mode दोलनों द्वारा सीधे मॉड्यूलेटेड अभिवक्क दर दिखाता है। टाइप A ($q=1/2$) और टाइप C (डायरेक्ट प्लंज) अलग व्यवहार दिखाते हैं। कुछ सिस्टम के लिए $g$-mode आवृत्तियाँ गामा-रे बर्स्ट में देखे गए हाई-फ्रीक्वेंसी क्यू-पी-ओ (QPOs) के अनुकूल हैं।
• किलोनोवा उत्सर्जन (Kilonova Emission): मॉडल किए गए किलोनोवा मुख्य रूप से डिस्क विंड इजेक्टा द्वारा संचालित होते हैं। कठोर DD2 EOS वाले सिस्टम नरम SLy EOS वाले सिस्टम की तुलना में अधिक चमकीले किलोनोवा उत्पन्न करते हैं क्योंकि उनमें अधिक विशाल डिस्क होती है। 200 Mpc की दूरी के लिए, अधिकांश सिस्टम वेरा सी रुबिन ऑब्जर्वेटरी द्वारा चार दिनों के भीतर और DECam द्वारा दो दिनों के भीतर $r$ और $i$ बैंड में पता लगाने योग्य होंगे।

महत्व और दावे
यह शोध पत्र सममित BHNS विलय के NR सिमुलेशन पैरामीटर स्पेस में एक महत्वपूर्ण अंतर को भरने का दावा करता है। प्राथमिक महत्व इस बात में निहित है कि वर्तमान गुरुत्वाकर्षण तरंग मॉडल समान और लगभग समान द्रव्यमान अनुपातों के लिए अपर्याप्त रूप से सटीक हैं, जिससे GW230529 जैसी भविष्य की खोजों के लिए पैरामीटर अनुमान (parameter estimation) में पूर्वाग्रह पैदा हो सकता है। अध्ययन अवशेष द्रव्यमान के लिए विश्लेषणात्मक फिटिंग सूत्रों को मान्य करता है लेकिन उच्च-टाइडल-डिफॉर्मेबिलिटी सिस्टम के लिए अवशेष स्पिन की भविष्यवाणी करने में उनकी सीमाओं को भी रेखांकित करता है। इसके अलावा, डिस्क गतिशीलता का विस्तृत विश्लेषण जटिल रोटेशन नियमों और दोलन-संचालित अभिवक्क मॉड्यूलेशन को प्रकट करता है, जो यह सुझाव देता है कि $g$-modes गामा-रे बर्स्ट लाइटकर्व पर दृश्य छाप छोड़ सकते हैं। अंततः, किलोनोवा मॉडलिंग इंगित करती है कि सममित BHNS विलय विद्युत चुम्बकीय फॉलो-अप के लिए आशाजनक लक्ष्य हैं, जो पिछले अवलोकनों में हावी रहे विषम सिस्टम की तुलना में उच्च पहचान संभावना प्रदान करते हैं।