Radiation hydrodynamic simulations for the origin of quasi-periodic oscillations for accretion onto supermassive black holes

यहाँ Lin et al. (2026) के शोध पत्र "Radiation hydrodynamic simulations for the origin of quasi-periodic oscillations for accretion onto supermassive black holes" का विस्तृत तकनीकी सारांश दिया गया है।

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

विभिन्न सुपरमैसिव ब्लैक होल (SMBH) प्रणालियों, जिनमें सक्रिय गैलेक्टिक नाभिक (AGN) और टाइडल डिसरप्शन इवेंट्स (TDEs) शामिल हैं, में अर्ध-आवधिक दोलन (Quasi-periodic oscillations - QPOs) देखे गए हैं। दशकों के सैद्धांतिक कार्य और संख्यात्मक सिमुलेशन (जो स्यूडो-न्यूटनियन हाइड्रोडायनामिक्स से लेकर जनरल रिलेटिविस्टिक मैग्नेटोहाइड्रोडायनामिक्स (GRMHD) तक फैले हुए हैं) के बावजूद, SMBH प्रणालियों में इन QPOs का भौतिक मूल अभी भी अस्पष्ट है।

• साहित्य में अंतराल (Gap in Literature): पिछले अध्ययन मुख्य रूप से स्टेलर-मास ब्लैक होल्स पर केंद्रित रहे हैं या विकिरण (radiation) और गुरुत्वाकर्षण के सरलीकृत उपचारों का उपयोग करते रहे हैं। बहुत कम सिमुलेशनों ने स्व-सुसंगत रेडिएशन हाइड्रोडायनामिक्स (RHD) का उपयोग करके SMBHs को संबोधित किया है जो रेडिएटिव ट्रांसफर के साथ एक्रेशन फ्लो डायनेमिक्स को जोड़ता है।
• उद्देश्य: लेखक उन्नत RHD सिमुलेशन करने और परिणामों की विश्लेतिक सिद्धांतों एवं अवलोकन संबंधी डेटा के साथ तुलना करके, SMBH एक्रेशन फ्लो में QPOs को संचालित करने वाले भौतिक तंत्र को निर्धारित करने का लक्ष्य रखते हैं।

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

लेखकों ने Athena++ कोड का उपयोग करके दो-आयामी अक्षीय सममित विस्कस रेडिएशन हाइड्रोडायनामिक (RHD) सिमुलेशन करने के लिए किया।

• भौतिक सेटअप (Physical Setup):
  • ब्लैक होल द्रव्यमान ($M_{BH}$): श्वासरशिल्ड ब्लैक होल के आसपास GR प्रभावों की नकल करने के लिए पैसिनी-विटा (Paczyński-Wiita) के स्यूडो-न्यूटनियन (PN) पोटेंशियल का उपयोग करते हुए इसे $10^7 M_\odot$ पर स्थिर रखा गया है: $\phi = -GM_{BH}/(r - R_S)$।
  • विस्कोसिटी (Viscosity): $\alpha = 0.1$ के साथ एक $\alpha$-विस्कोसिटी मॉडल अपनाया गया।
  • विकिरण (Radiation): रेडिएटिव ट्रांसफर के लिए सबसे उन्नत उपचार (Jiang 2021, 2022) का उपयोग किया गया, जो हाइड्रोडायनामिक्स को रेडिएटिव ट्रांसफर के साथ स्व-सुसंगत रूप से जोड़ता है।
• सिमुलेशन पैरामीटर्स:
  • द्रव्यमान इंजेक्शन (Mass Injection): गैस को एक निश्चित त्रिज्या $R_{inject} = 10 R_S$ (श्वार्ज़स्चाइल्ड त्रिज्या) पर स्थानीय कीपलरियन वेग के बराबर घूर्णी वेग के साथ इंजेक्ट किया गया था।
  • इंजेक्शन दरें (Injection Rates): तीन अलग-अलग द्रव्यमान एक्रेशन दरों का परीक्षण किया गया:
    1. $\dot{M}_{inject} = 0.01 \dot{M}_{Edd}$ (ज्यामितीय रूप से पतला, ऑप्टिकली थिक डिस्क)।
    2. $\dot{M}_{inject} = 0.1 \dot{M}_{Edd}$ (ज्यामितीय रूप से पतला, ऑप्टिकली थिक डिस्क)।
    3. $\dot{M}_{inject} = 2 \dot{M}_{Edd}$ (ज्यामितीय रूप से मोटा, ऑप्टिकली थिक "स्लिम" डिस्क)।
  • डोमेन (Domain): रेडियल रेंज $2 R_S \leq r \leq 10^5 R_S$ जिसमें $768 \times 256$ रेजोल्यूशन का उपयोग किया गया। सिमुलेशन 400 दिनों (लगभग 1394 ऑर्बिटल पीरियड, $10 R_S$ पर) के लिए चलाए गए।
• विश्लेषण (Analysis):
  • विभिन्न त्रिज्याओं पर द्रव्यमान प्रवाह दर ($\dot{M}_{in}$) की गणना की गई।
  • QPO आवृत्तियों की पहचान करने के लिए पावर स्पेक्ट्रा और डायनेमिकल पावर स्पेक्ट्रा उत्पन्न करने के लिए फास्ट फूरियर ट्रांसफॉर्म (FFT) लागू किया गया।
  • सिमुलेशन-व्युत्पन्न आवृत्तियों की तुलना सैद्धांतिक रेडियल एपिकिलिक आवृत्तियों ($\nu_r$) से की गई।

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

A. एक्रेशन फ्लो का विकास और QPO का उद्भव

• कम/मध्यम एक्रेशन दरें ($0.01$और$0.1 \dot{M}_{Edd}$): प्रवाह एक स्थिर, ज्यामितीय रूप से पतली डिस्क में स्थापित हो जाता है। प्रारंभिक क्षणिक चरण (जब इंजेक्शन त्रिज्या और BH के बीच का अंतर भर जाता है) के बाद QPO सिग्नल उभरते हैं।
• उच्च एक्रेशन दर ($2 \dot{M}_{Edd}$): प्रवाह रेडिएशन प्रेशर अस्थिरता द्वारा संचालित उच्च और निम्न अवस्थाओं के बीच महत्वपूर्ण स्टेट ट्रांजिशन (state transitions) प्रदर्शित करता है।
  • उच्च अवस्था (High State): डिस्क ज्यामितीय रूप से मोटी हो जाती है और इसमें मजबूत आउटफ्लो होते हैं; रेडियल दोलन दब जाते हैं, और QPOs कमजोर या अनुपस्थित होते हैं।
  • निम्न अवस्था (Low State): प्रवाह भूमध्यरेखीय तल (equatorial plane) तक सीमित रहता है; रेडियल दोलन पुनः प्राप्त होते हैं, और मजबूत QPOs दिखाई देते हैं।
• निष्कर्ष: QPOs द्रव्यमान प्रवाह दर में मजबूती से पाए जाते हैं, विशेष रूप से सुपर-एडिंगटन एक्रेशन की निम्न अवस्था में और स्थिर पतली डिस्क शासन (regimes) में।

B. आवृत्ति विश्लेषण और सैद्धांतिक निरंतरता

• फ्रीक्वेंसी मैचिंग: $r > 3.8 R_S$ पर सिमुलेशन से प्राप्त QPO आवृत्तियाँ ($\nu_{QPO}$) विश्लेतिक रूप से गणना की गई रेडियल एपिकिलिक फ्रीक्वेंसी ($\nu_r$) के साथ अच्छी तरह से सुसंगत हैं।
• क्रिटिकल रेडियस: सिमुलेशन एक क्रिटिकल रेडियस $r \approx 3.8 R_S$ (विशेष रूप से $(2+\sqrt{3})R_S$) प्रकट करता है, जहाँ रेडियल एपिकिलिक फ्रीक्वेंसी अपने अधिकतम मान ($\nu_{r,max}$) पर पहुँचती है।
• प्रसार (Propagation): $r < 3.8 R_S$ के लिए, QPO फ्रीक्वेंसी लगभग स्थिर रहती है (जो $\nu_{r,max}$ से मेल खाती है), लेकिन इसकी शक्ति (power) काफी कम हो जाती है। लेखक प्रस्तावित करते हैं कि $3.8 R_S$ पर अधिकतम एपिकिलिक फ्रीक्वेंसी आंतरिक सीमा की ओर प्रसारित होती है, जो देखे गए QPOs के चालक के रूप में कार्य करती है।

C. अवलोकनों के साथ स्केलिंग (Scaling to Observations)

• लेखकों ने स्केलिंग संबंध $\nu_{r,max} \propto M_{BH}^{-1}$ का परीक्षण किया।
• AGNs (RE J1034+396, 1ES 1927+654) और TDEs (Swift J164449.3+573451, ASASSN-14li) के नमूनों के लिए देखी गई QPO आवृत्तियों को ब्लैक होल द्रव्यमान के विरुद्ध प्लॉट करके, उन्होंने पाया कि $\nu_{r,max}$ का सैद्धांतिक वक्र अवलोकनात्मक डेटा के साथ अच्छी तरह मेल खाता है।
• अपवाद: ASASSN-14li उच्च फ्रीक्वेंसी दिखाता है (जो मध्यम ब्लैक होल स्पिन का सुझाव देता है), और 2XMM J123103.2+110648 कम फ्रीक्वेंसी दिखाता है (संभवतः थर्मल अस्थिरता के कारण), लेकिन सामान्य रुझान $\nu_{r,max}$ मॉडल का समर्थन करता है।

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

1. स्व-सुसंगत RHD सिमुलेशन: यह उन पहले अध्ययनों में से एक है जो SMBH एक्रेशन के लिए पूरी तरह से युग्मित रेडिएशन हाइड्रोडायनामिक्स का उपयोग करके QPOs का अध्ययन करने के लिए किया गया है, जो रेडिएशन और डायनेमिक्स के अलग-अलग उपचारों से आगे बढ़ता है।
2. उत्पत्ति की पहचान: यह अध्ययन इस बात के पुख्ता संख्यात्मक प्रमाण प्रदान करता है कि रेडियल एपिकिलिक ऑसिलेशन्स (radial epicyclic oscillations) ही SMBH प्रणालियों में QPOs का भौतिक मूल हैं।
3. प्रॉक्सी प्रस्ताव: लेखक प्रस्तावित करते हैं कि $r \approx 3.8 R_S$ पर अधिकतम रेडियल एपिकिलिक फ्रीक्वेंसी ($\nu_{r,max}$), देखी गई QPO फ्रीक्वेंसी ($\nu_{QPO}$) के लिए एक मजबूत प्रॉक्सी के रूप में कार्य करती है।
4. स्टेट-डिपेंडेंट QPOs: सुपर-एडिंगटन प्रवाह में QPOs का अवस्था-निर्भर होना (निम्न अवस्थाओं में मौजूद, वर्टिकल डिस्क विस्तार के कारण उच्च अवस्थाओं में दमित) एक नया स्पष्टीकरण प्रदान करता है, जो TDEs की परिवर्तनशीलता (variability) को समझाता है।

तुलनात्मक महत्व और चर्चा (Significance and Discussion)

• भौतिक तंत्र: परिणाम स्टेलर-मास और सुपरमैसिव ब्लैक होल्स दोनों में QPOs की समझ को एकीकृत करते हैं, जो रेडियल एपिकिलिक रेजोनेंस पर आधारित एक सार्वभौमिक तंत्र का सुझाव देते हैं।
• ब्लैक होल गुणों पर प्रतिबंध: चूंकि $\nu_{r,max}$ द्रव्यमान ($M_{BH}$) और स्पिन पैरामीटर ($a$) पर निर्भर करता है, इसलिए यह मॉडल एक संभावित विधि प्रदान करता है जिससे यदि द्रव्यमान स्वतंत्र रूप से ज्ञात हो, तो SMBH स्पिन को सीमित किया जा सकता है।
• सीमाएँ और भविष्य का कार्य:
  • गुरुत्वाकर्षण: इस अध्ययन में PN पोटेंशियल का उपयोग किया गया है। हालांकि यह इस उद्देश्य के लिए GR का अच्छा अनुमान लगाता है, भविष्य के कार्य के लिए झुकी हुई डिस्क (tilted disks) में फ्रेम-ड्रैगिंग (Lense-Thirring precession) प्रभावों को पकड़ने के लिए पूर्ण GRMHD की आवश्यकता होगी।
  • चुंबकीय क्षेत्र: चुंबकीय क्षेत्रों की उपेक्षा की गई है। QPO दृश्यता पर चुंबकीय हवाओं (magnetic winds) के प्रभाव का आकलन करने के लिए भविष्य में RMHD सिमुलेशन की आवश्यकता है।
  • अवलोकन योग्य (Observables): वर्तमान अध्ययन द्रव्यमान प्रवाह दर का विश्लेषण करता है। भविष्य के कार्यों को व्यूइंग एंगल प्रभावों और आउटफ्लो अवरोध (outflow obscuration) को ध्यान में रखने के लिए, जो QPO संकेतों को कम कर सकते हैं, वैश्विक उत्सर्जित स्पेक्ट्रा (global emergent spectra) और लाइट कर्व्स की गणना करनी चाहिए।

संक्षेप में, Lin et al. (2026) यह प्रदर्शित करते हैं कि रेडिएशन हाइड्रोडायनामिक सिमुलेशन स्वाभाविक रूप से रेडियल एपिकिलिक ऑसिलेशन्स द्वारा संचालित QPOs उत्पन्न करते हैं, जो एक सम्मोहक सैद्धांतिक ढांचा प्रदान करते हैं जो विभिन्न SMBH द्रव्यमान और एक्रेशन अवस्थाओं के माध्यम से अवलोकनात्मक डेटा के साथ संरेखित होता है।