Perturber-Driven Dynamics of Supermassive Black Hole Binaries in Galaxy Merger

तकनीकी सारांश: गैलेक्सी विलय में सुपरमैसिव ब्लैक होल बाइनरीज की विचलित-प्रेरित गतिशीलता (Perturber-Driven Dynamics)

समस्या विवरण
सुपरमैसिव ब्लैक होल बाइनरीज (MBHBs) के निर्माण के समय उनकी कक्षीय उत्केंद्रता (orbital eccentricity) एक महत्वपूर्ण पैरामीटर है, जो पल्सर टाइमिंग एरेज़ (PTAs) द्वारा पता योग्य गुरुत्वाकर्षण तरंग पृष्ठभूमि (GWB) को निर्धारित करती है। पिछले उच्च-रिज़ॉल्यूशन $N$-बॉडी सिमुलेशन ने प्रदर्शित किया है कि निर्माण उत्केंद्रता (formation eccentricity) में रन-टू-रन महत्वपूर्ण बिखराव (scatter) होता है, विशेष रूप से लगभग रेडियल विलय कक्षाओं के लिए। जबकि यह दिखाया गया है कि यह बिखराव रिज़ॉल्यूशन के साथ $\sim 1/\sqrt{N}$ के रूप में कम होता है, जो पॉइसन सैंपलिंग शोर (Poisson sampling noise) के अनुरूप है, यह स्पष्ट नहीं है कि क्या भौतिक उप-संरचना—विशेष रूप से भारी विचलितकर्ता (massive perturbers) जैसे कि ग्लोबुलर क्लस्टर, स्ट्रिप्ड न्यूक्लिआई, या इंटरमीडिएट-मास ब्लैक होल्स (IMBHs)—एक वास्तविक खगोलीय स्टोकेस्टिसिटी (astrophysical stochasticity) पेश करते हैं जो संख्यात्मक शोर के स्तर (numerical noise floor) से अधिक हो। मुख्य प्रश्न यह है कि क्या कई कम द्रव्यमान वाले विचलितकर्ताओं का संचयी प्रभाव (डिफ्यूसिव रिजीम) या दुर्लभ, भारी मुठभेड़ (इम्पल्सिव रिजीम) बड़े दीर्घवृत्ताकार (elliptical) आकाशगंगाओं में MBHBs की निर्माण उत्केंद्रता को महत्वपूर्ण रूप से यादृच्छिक (randomize) कर सकते हैं।

कार्यप्रणाली
लेखकों ने $Griffin$ $N$-body कोड का उपयोग करके एक प्रमुख गैलेक्सी विलय (IllustrisTNG100-1 ब्रह्मांडीय सिमुलेशन से Merger 6) के उच्च-रिज़ॉल्यूशन पुन: सिमुलेशन आयोजित किए। सिमुलेशन सेटअप में शामिल था:

• प्रारंभिक स्थितियाँ: दो आकाशगंगाओं के बीच एक लगभग रेडियल विलय ($e_{orb} \approx 0.987$), जिन्हें हर्नक्विस्ट प्रोफाइल (Hernquist profiles) का पालन करने वाले कोलिजनलेस सिस्टम (तारकीय बल्ज, डार्क मैटर हेलो और केंद्रीय MBHs) के रूप में मॉडल किया गया था। प्राथमिक आकाशगंगा ने एक छोटा गैस अंश ($f_{gas} \approx 0.19$) बनाए रखा, लेकिन सिमुलेशन ने लंबे समय तक चलने वाले तारकीय सिस्टमों पर ध्यान केंद्रित करने के लिए सिस्टम को कोलिजनलेस माना।
• विचलितकर्ता परिदृश्य (Perturber Scenarios): तीन सूट के सिमुलेशन की तुलना की गई, जिनमें से प्रत्येक में चार स्वतंत्र वास्तविक स्वरूप (realisations) थे:
  1. कंट्रोल (Control): कोई भारी विचलितकर्ता नहीं; चिकना तारकीय बल्ज बैकग्राउंड के रूप में कार्य करता है।
  2. $10^7 M_\odot$ सूट: $3,162$समान-द्रव्यमान विचलितकर्ता ($M_p = 10^7 M_\odot$) जो प्राथमिक बल्ज द्रव्यमान का$10%$पुनर्वित्रित करते हैं। यह एक विचलितकर्ता-से-बाइनरी द्रव्यमान अनुपात$\mu_p \approx 3.2 \times 10^{-3}$ के अनुरूप है।
  3. $10^8 M_\odot$ सूट: $316$समान-द्रव्यमान विचलितकर्ता ($M_p = 10^8 M_\odot$) जो समान कुल द्रव्यमान अंश को पुनर्वित्रित करते हैं। यह$\mu_p \approx 3.2 \times 10^{-2}$ के अनुरूप है।
• बाइनरी निर्माण मानदंड: बाइनरी निर्माण समय ($t_b$) को स्थायी गुरुत्वाकर्षण बंधन के बाद उत्केंद्रता के पहले स्थानीय न्यूनतम (local minimum) के रूप में परिभाषित किया गया था। निर्माण उत्केंद्रता ($e_b$) और उनके विभिन्न वास्तविक स्वरूपों में मानक विचलन ($\sigma_e$) प्राथमिक मेट्रिक्स थे।
• गतिक निदान (Dynamical Diagnostics): अध्ययन ने $t_b$ के साथ संरेखित, मैच किए गए कंट्रोल रन के सापेक्ष कक्षीय ऊर्जा ($E$), कोणीय संवेग ($J_z$), टॉर्क ($\tau_z$), और बाइनरी सेंटर ऑफ मास (COM) के विचलन (wander) में अवशेषों (residuals) का विश्लेषण किया।

मुख्य योगदान और परिणाम
यह शोध पेपर विचलितकर्ता-से-बाइनरी द्रव्यमान अनुपात ($\mu_p$) के आधार पर डिफ्यूसिव और निकट-इम्पल्सिव स्कैटरिंग रिजीम के बीच एक संक्रमण स्थापित करता है:

1. संख्यात्मक शोर स्तर (Numerical Noise Floor): नो-पर्टर्बर कंट्रोल सूट ने एक निर्माण उत्केंद्रता बिखराव $\sigma_e \approx 0.11$ दिया, जो अपनाए गए रिज़ॉल्यूशन ($N_{bulge} \approx 2 \times 10^6$) के लिए अपेक्षित पॉइसन शोर स्तर के अनुरूप है।
2. डिफ्यूसिव रिजीम ($10^7 M_\odot$): $10^7 M_\odot$ सूट ने $\sigma_e \approx 0.115$ उत्पन्न किया, जो सांख्यिकीय रूप से कंट्रोल से अभिन्न (indistinguishable) है। ऊर्जा और कोणीय संवेग में गतिक अवशेष सुचारू थे, और कोई महत्वपूर्ण टॉर्क स्पाइक्स नहीं देखे गए। यह पुष्टि करता है कि $\mu_p \approx 3.2 \times 10^{-3}$ के लिए, सिस्टम डिफ्यूसिव रिजीम में रहता है जहाँ व्यक्तिगत मुठभेड़ इतनी कमजोर होती है कि वे निर्माण उत्केंद्रता पर कोई दृश्य छाप नहीं छोड़ पाती हैं।
3. निकट-इम्पल्सिव रिजीम ($10^8 M_\odot$): $10^8 M_\odot$ सूट ने काफी अधिक बिखराव प्रदर्शित किया, $\sigma_e \approx 0.26$ (कंट्रोल से $\sim 2.4$ गुना अधिक)। इस अतिरिक्त बिखराव के साथ निम्नलिखित देखे गए:
   • अर्ध-दीप्त अक्ष (semi-major axis) और कक्षीय ऊर्जा में विविक्त उछाल (discrete jumps)।
   • $t_b$ के पास कोणीय संमेद में बड़े, इवेंट-अलाइन्ड अवशेष और टॉर्क स्पाइक्स।
   • बाइनरी सेंटर ऑफ मास का बढ़ा हुआ विचलन (कंट्रोल की तुलना में $\sim 1.25$–$1.4\times$ अधिक)।
     ये संकेत दर्शाते हैं कि इस रिजीम में व्यक्तिगत करीबी मुठभेड़ इतने बड़े इम्पल्सिव किक (impulsive kicks) देती है जो हार्डनिंग चरण शुरू होने से पहले बाइनरी के कोणीय संवेग और उत्केंद्रता को बदल देते हैं।
4. सैद्धांतिक संक्रमण: परिणाम बाइनरी-सिंगल स्कैटरिंग थ्योरी के अनुरूप हैं, जहाँ संक्रमण $\mu_p$ द्वारा नियंत्रित होता है। $10^8 M_\odot$ का मामला निकट-इम्पल्सिव सीमा ($\Delta e_{single} \sim \mu_p$) के करीब पहुँच जाता है, जबकि $10^7 M_\odot$ का मामला डिफ्यूसिव बना रहता है।

महत्व और दावे
यह पेपर दावा करता है कि हालांकि भारी विचलितकर्ता (massive perturbers) वास्तव में MBHB निर्माण उत्केंद्रता में भौतिक स्टोकेस्टिसिटी उत्पन्न कर सकते हैं, यह प्रभाव विचलितकर्ताओं के द्रव्यमान पैमाने पर अत्यधिक निर्भर है जो बाइनरी के सापेक्ष है।

• खगोलीय संदर्भ: विशाल दीर्घवृत्ताकार आकाशगंगाओं (जो PTA-प्रासंगिक MBHBs के प्राथमिक होस्ट हैं) में, प्रमुख विचलितकर्ता आबादी ग्लोबुलर क्लस्टर ($10^5$–$10^6 M_\odot$) और स्ट्रिप्ड न्यूक्लिआई है। PTAs के लिए प्रासंगिक बाइनरी द्रव्यमानों ($10^8$–$10^{10} M_\odot$) के लिए, ये वस्तुएं $\mu_p \lesssim 10^{-4}$ के अनुरूप हैं, जो उन्हें दृढ़ता से डिफ्यूसिव रिजीम में रखती हैं।
• निष्कर्ष: फलस्वरूप, भारी विचलितकर्ता द्वारा संचालित उत्केंद्रता यादृच्छिककरण (randomization), GWB डिटेक्शन के लिए प्रासंगिक अधिकांश MBHB विलयों के लिए एक महत्वपूर्ण मार्ग होने की संभावना नहीं है। वर्तमान सिमुलेशन में देखा गया स्टोकेस्टिसिटी, वास्तविक विचलितकर्ता आबादी के लिए भौतिक उप-संरचना के बजाय संख्यात्मक शोर द्वारा हावी है।
• सीमाएँ (Caveats): लेखक नोट करते हैं कि उनका निष्कर्ष केवल एकल $10^8 M_\odot$ वृद्धि की सांख्यिकीय सार्थकता के बजाय द्रव्यमान-अनुपात मार्जिन पर निर्भर करता है (जो केवल चार वास्तविक स्वरूपों के साथ मामूली है)। वे यह भी स्वीकार करते हैं कि गैस-समृद्ध विलय (जो यहाँ पूरी तरह से मॉडल नहीं किए गए हैं) विशाल आणविक बादलों या सर्कमबाइनरी डिस्क के माध्यम से अलग गतिशीलता पेश कर सकते हैं, और विशिष्ट कक्षीय विन्यास (लगभग रेडियल) सटीक सीमा को प्रभावित कर सकता है। हालाँकि, प्रमुख खगोलीय निष्कर्ष यह है कि विशाल दीर्घवृत्ताकार आकाशगंगाओं में वास्तविक विचलितकर्ता आबादी, निर्माण उत्केंद्रता को महत्वपूर्ण रूप से यादृच्छिक करने के लिए आवश्यक निकट-इम्पल्सिव थ्रेशोल्ड से नीचे है।