Radiative GRMHD simulations of puffy accretion discs: Numerical versus analytical models of sub-Eddington accretion

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

एक्स-रे बाइनरी सिस्टम्स (XRBs) में चमकदार सॉफ्ट स्पेक्ट्रल स्टेट के दौरान अभिव्र्द्धि प्रवाह (accretion flows) के लिए मानक सैद्धांतिक ढांचा ज्यामितीय रूप से पतला और ऑप्टिकली थिक शाकुरा-सुनाएव डिस्क (नोविकोव-थॉर्न विस्तार) है। हालांकि, यह विश्लेषणात्मक मॉडल दो महत्वपूर्ण समस्याओं का सामना करता है:

1. तापीय अस्थिरता (Thermal Instability): मॉडल भविष्यवाणी करता है कि विकिरण दबाव (radiation pressure) द्वारा नियंत्रित डिस्क (जो उच्च-दीप्ति वाले सॉफ्ट स्टेट में सामान्य है) तापीय रूप से अस्थिर होती है, फिर भी अवलोकन दर्शाते हैं कि ये सिस्टम लंबे समय तक स्थिर रहते हैं।
2. अवलोकन संबंधी विसंगतियां (Observational Discrepancies): XRBs के सॉफ्ट स्टेट के हालिया पोलारिमेट्रिक अवलोकनों (जैसे, IXIXE द्वारा) ने उच्च ध्रुवीकरण डिग्री (polarization degrees) प्रकट की है जिसे मानक पतली डिस्क मॉडल स्पष्ट नहीं कर सकते, भले ही स्पेक्ट्रल फिटिंग सटीक हो। इसके अतिरिक्त, ब्लैक होल (BH) स्पिन के अनुमान के लिए 'इनरमोस्ट स्टेबल सर्कुलर ऑर्बिट' (ISCO) पर आधारित विधियां गलत हो सकती हैं यदि डिस्क की संरचना विश्लेषणात्मक धारणाओं से विचलित होती है।

जबकि चुंबकीय क्षेत्रों को एक स्थिर करने वाले तंत्र के रूप में सुझाया गया है, विश्लेषणात्मक मॉडल अक्सर विकिरण, चुंबकीय क्षेत्र और ऊर्ध्वाधर स्तरीकरण (vertical stratification) के बीच जटिल परस्पर क्रिया को पकड़ने में विफल रहते हैं। सब-एडिंगटन अभिवृद्धि प्रवाहों की वास्तविक संरचना को समझने के लिए एक स्व-संगत (self-consistent), 3D संख्यात्मक सिमुलेशन की आवश्यकता है।

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

लेखकों ने KORAL कोड का उपयोग करके जनरल रिलेटिविस्टिक रेडिएटिव मैग्नेटोहाइड्रोडायनामिक (GRRMHD) सिमुलेशन का उपयोग किया।

• सिमुलेशन सेटअप:

  • सिस्टम: एक गैर-घूर्णन (श्वार्ज़चाइल्ड) स्टेलर-मास ब्लैक होल ($M = 10 M_\odot$)।
  • प्रारंभिक स्थितियां: एक मोटा, सब-केपलरियन टॉरस जिसमें क्वाड्रुपोल चुंबकीय क्षेत्र टोपोलॉजी (प्रेशर अनुपात $p_g/p_m = 20$) है, जो BH से दूर ($r_{in} = 35 r_g$) स्थित है।
  • भौतिकी: कोड M1 क्लोजर स्कीम का उपयोग करके रेडिएशन ट्रांसपोर्ट के साथ युग्मित रिलेटिविस्टिक आइडियल MHD समीकरणों को हल करता है। इसमें फ्री-फ्री अवशोषण, थॉमसन स्कैटरिंग, और सिनक्रोट्रॉन/ब्रेम्सस्ट्रालुंग ओपेसिटी शामिल है।
  • पैरामीटर्स: चार सिमुलेशन अलग-अलग मास एक्रेशन रेट ($\dot{m} = \dot{M}/\dot{M}_{Edd}$) के साथ किए गए: 0.45, 0.58, 0.89, और 1.51। मुख्य ध्यान सिमुलेशन S06 ($\dot{m} = 0.58$) पर है।
  • अवधि: सिमुलेशन को एक अर्ध-स्थिर अवस्था (quasi-steady state) तक पहुँचने के लिए दसियों हजार गुरुत्वाकर्षण समय इकाइयों ($t_g$) तक चलाया गया, जिसमें डेटा को समय और अज़ीमुथ (azimuth) के ऊपर औसत निकाला गया।

• तुलनात्मक विश्लेषण:

  • सिमुलेशन परिणामों की तुलना मानक विश्लेषणात्मक मॉडलों के साथ की गई: थिन डिस्क (नोविकोव-थॉर्न) और स्लिम डिस्क (अब्रामोविच एट अल.), विशेष रूप से ऊर्ध्वाधर संरचना, सतह घनत्व (surface density), तापमान और प्रभावी श्यानता (effective viscosity) का परीक्षण किया गया।

3. प्रमुख योगदान और परिणाम (Key Contributions and Results)

A. "पफी डिस्क" संरचना की खोज (Discovery of the "Puffy Disc" Structure)

सिमुलेशन एक स्थिर अभिवृद्धि प्रवाह प्रकट करते हैं जो ज्यामितीय रूप से मोटा ($h_\tau/R \sim 1$) है, जो पतली डिस्क की धारणा के विपरीत है। संरचना तीन विशिष्ट क्षेत्रों में ऊर्ध्वाधर रूप से स्तरित है:

1. कोर (C): भूमध्यरेखीय तल के पास एक सघन, ऑप्टिकली थिक, ज्यामितीय रूप से पतला क्षेत्र ($h_\rho/R \approx 0.1$)।
2. पफी क्षेत्र (P): कोर के चारों ओर एक ज्यामितीय और ऑप्टिकली थिक "वार्म कोरोना", जो फोटोस्फीयर तक फैला हुआ है।
3. फनल (F): घूर्णन अक्ष के साथ एक ऑप्टिकली थिन क्षेत्र।

• स्थिरीकरण: डिस्क को चुंबकीय दबाव (magnetic pressure) द्वारा तापीय अस्थिरता के विरुद्ध स्थिर किया जाता है। चुंबकीय दबाव ($p_m$), गैस दबाव ($p_g$) पर हावी है और विकिरण दबाव ($p_r$) के तुल्य है। महत्वपूर्ण रूप से, स्थिरता के लिए एक क्वाड्रुपोल चुंबकीय क्षेत्र टोपोलॉजी (जो एक शुद्ध ऊर्ध्वाधर/त्रिज्यीय फ्लक्स प्रदान करती है) आवश्यक है; द्विध्रुवीय (dipolar) क्षेत्र प्रवाह को स्थिर करने में विफल रहते हैं।

B. विश्लेषणात्मक मॉडलों से विचलन (Deviations from Analytic Models)

पफी डिस्क मानक विश्लेषणात्मक भविष्यवाणियों से महत्वपूर्ण अंतर प्रदर्शित करती है:

• सतह घनत्व (Surface Density): सिम्युलेटेड सतह घनत्व ($\Sigma$) विश्लेषणात्मक मॉडलों द्वारा अनुमानित मान से काफी कम है, जिसका अर्थ है उच्च अभिवृद्धि वेग।
• तापमान: केंद्रीय गैस तापमान ($T_c$) आंतरिक क्षेत्र में उतना नहीं गिरता जितना विश्लेषणात्मक मॉडल भविष्यवाणी करते हैं; इसके बजाय, यह उच्च बना रहता है क्योंकि डिस्क इवेंट होराइजन तक ऑप्टिकली थिक रहती है।
• आंतरिक किनारे का स्थान (Inner Edge Location):
  • विश्लेषणात्मक मॉडल मानते हैं कि आंतरिक किनारा ISCO (जहाँ स्ट्रेस शून्य हो जाता है) या सोनिक पॉइंट पर है।
  • पफी डिस्क में, प्रवाह ISCO के काफी अंदर तक सबसोनिक और ऑप्टिकली थिक रहता है। मैग्नेटोसोनिक पॉइंट (आंतरिक किनारा) ISCO ($6 r_g$) की तुलना में इवेंट होराइजन के काफी करीब ($R_{sp} \approx 3.9 r_g$ S06 के लिए) स्थित है।
  • इसका तात्पर्य है कि BH स्पिन अनुमान के लिए मानक स्पेक्ट्रल फिटिंग विधियां, जो मानती हैं कि डिस्क ISCO पर समाप्त हो जाती है, गलत (अति-अनुमानित) स्पिन मान दे सकती हैं।

C. प्रभावी श्यानता ($\alpha$) (Effective Viscosity)

• गैर-स्थिर श्यानता: विश्लेषणात्मक मॉडलों में एक स्थिर $\alpha$ पैरामीटर की धारणा के विपरीत, सिमुलेशन दिखाते हैं कि $\alpha$ स्थिर नहीं है।
• त्रिज्यीय प्रोफाइल: $\alpha$ ब्लैक होल के निकटतम आंतरिक क्षेत्रों में तेजी से बढ़ता है।
• तंत्र: श्यानता मुख्य रूप से मैक्सवेल स्ट्रेस (चुंबकीय विक्षोभ/turbulence) द्वारा संचालित होती है, न कि रेनॉल्ड्स स्ट्रेस द्वारा। आंतरिक भाग में $\alpha$ का बढ़ना प्रवाह के पतला होने और चुंबकीय तनाव के तीव्र होने से जुड़ा है।

D. एडवेक्शन और फोटॉन ट्रैपिंग (Advection and Photon Trapping)

• फोटॉन ट्रैपिंग: उच्च ऑप्टिकल गहराई और ऊर्ध्वाधर मोटाई के कारण, फोटॉन ट्रैप हो जाते हैं और $\sim 20 r_g$ तक प्रवाह के साथ अंदर की ओर एडवेक्ट (advect) होते हैं।
• स्टैग्नेशन सतह (Stagnation Surface): फनल क्षेत्र में एक गैस स्टैग्नेशन सतह मौजूद है जहाँ त्रिज्यीय वेग अपना चिह्न बदल देता है, जो इनफ्लो (inflow) को आउटफ्लो (outflow) से अलग करता है।

4. महत्व और निहितार्थ (Significance and Implications)

• अस्थिरता का समाधान: यह अध्ययन प्रदर्शित करता है कि चुंबकीय क्षेत्र विकिरण-दबाव-प्रधान अभिवृद्धि प्रवाह को स्थिर कर सकते हैं, जो XRB सॉफ्ट स्टेट्स में देखे गए दीर्घकालिक स्थिरता की व्याख्या करता है।
• अवलोकनों की पुनर्व्याख्या:
  • ध्रुवीकरण (Polarization): ज्यामितीय रूप से मोटा, स्तरित संरचना और फोटोस्फीयर के पास मजबूत ऊर्ध्वाधर इनफ्लो वेग, IXE द्वारा देखी गई उच्च ध्रुवीकरण डिग्री के लिए एक प्राकृतिक स्पष्टीकरण प्रदान करता है, जिसे पतली डिस्क मॉडल पुनरुत्पादित नहीं कर सकते।
  • स्पिन माप: यह निष्कर्ष कि आंतरिक किनारा ISCO के काफी अंदर स्थित है, यह सुझाव देता है कि थर्मल कॉन्टिन्यूम फिटिंग के माध्यम से BH स्पिन मापने की वर्तमान विधियां व्यवस्थित रूप से पक्षपाती (biased) हो सकती हैं।
• मॉडल की सीमाएं: पेपर रेखांकित करता है कि मानक विश्लेषणात्मक और अर्ध-विश्लेषणात्मक मॉडल (थिन और स्लिम डिस्क) चमकदार, सब-एडिंगटन अभिवृद्धि की जटिल भौतिकी का वर्णन करने के लिए अपर्याप्त हैं। यद्यपि वे गणनात्मक रूप से महंगे हैं, वास्तविक ऊर्ध्वाधर संरचना और अवलोकन संबंधी संकेतों को पकड़ने के लिए 3D GRRMHD सिमुलेशन आवश्यक हैं।
• अन्य घटनाओं से संबंध: "पफी" संरचना डिस्क-कोरोना मॉडलों में देखे गए "वार्म कोरोना" के समान है और सुपर-एडिंगटन प्रवाह (ULXs) और 'पोलिश डोनट्स' में देखे गए बीमिंग प्रभावों के साथ समानता साझा करती है, जो विभिन्न अभिवृद्धि व्यवस्थाओं के लिए एक एकीकृत भौतिक चित्र का सुझाव देती है।

निष्कर्ष (Conclusion)

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