Population synthesis of active galactic nuclei based on the radiation-regulated unification model

तकनीकी सारांश: विकिरण-विनियमित एकीकरण मॉडल पर आधारित सक्रिय गैलेक्टिक नाभिक (AGN) का जनसंख्या संश्लेषण

समस्या विवरण
सक्रिय गैलेक्टिक नाभिक (AGNs) सुपरमैसिव ब्लैक होल (SMBHs) द्वारा संचालित होते हैं जिनका विकास मेजबान आकाशगंगा के गुणों से जुड़ा हुआ है। जबकि एक्स-रे सर्वेक्षण एक्स-रे उत्सर्जन, अवरोध (obscuration) और अभिवृद्धि (accretion) दरों पर महत्वपूर्ण बाधाएं प्रदान करते हैं, पिछले जनसंख्या संश्लेषण मॉडल प्रतिबिंब (reflection) और अवशोषण (absorption) को स्वतंत्र कारकों के रूप में मानते थे, जो अक्सर स्व-संगत भौतिक मॉडलिंग के बजाय घटनात्मक टेम्प्लेट (phenomenological templates) पर निर्भर करते थे। इसके अलावा, कॉम्पटन-थिक (CTK) AGNs ($N_H \ge 10^{24} \text{ cm}^{-2}$) का आंतरिक अंश विवादास्पद बना हुआ है, जिसके मॉडल भविष्यवाणियां 10% से 50% के बीच रहती हैं। एक प्रमुख चुनौती अवरोध और चमक (luminosity) के बीच प्रति-सहसंबंध (anti-correlation) को "एकीकरण मॉडल" (unification model) के साथ सामंजस्य बिठाना है, जो यह प्रस्तावित करता है कि केवल अभिविन्यास (orientation) ही AGN वर्गीकरण को निर्धारित करता है। हालिया प्रस्ताव सुझाव देते हैं कि एक विकिरण-विनियमित एकीकरण मॉडल जहाँ अवरोधक सामग्री का कवरेज फैक्टर (covering factor) एडिंगटन अनुपात ($\lambda_{\text{Edd}}$) द्वारा संचालित होता है, न कि चमक द्वारा, क्योंकि विकिरण दबाव परिनाभिकीय वातावरण (circumnuclear environment) को साफ कर देता है। यह कार्य एक व्यापक एक्स-रे अवलोकन सेट के विरुद्ध इस विकिरण-विनियमित मॉडल की निरंतरता का परीक्षण करने का लक्ष्य रखता है और स्व-संगत रूप से डस्टी टोरस (dusty torus) की ज्यामिति और स्थानीय AGN जनसंख्या के विकास को नियंत्रित करता है।

कार्यप्रणाली
लेखक AGN जनसंख्या के ज्यामितीय और भौतिक गुणों को नियंत्रित करने के लिए एक सिमुलेशन-आधारित अनुमान (simulation-based inference - SBI) दृष्टिकोण का उपयोग करते हैं। कार्यप्रणाली में निम्नलिखित चरण शामिल हैं:

1. सिंथेटिक जनसंख्या निर्माण: एक्स-रे ल्यूमिनोसिटी फंक्शन (XLF) से नमूना लेने के बजाय, सिंथेटिक जनसंख्या को $z < 0.3$ के लिए अनना एट अल. (2022, A22) द्वारा व्युत्पन्न आंतरिक स्थानीय सक्रिय ब्लैक होल मास फंक्शन (BHMF) और एडिंगटन अनुपात वितरण फंक्शन (ERDF) से नमूना लेकर बनाया जाता है। एक सरल विकास नुस्खा (evolutionary prescription) लागू किया जाता है जहाँ सक्रिय SMBHs की संख्या घनत्व रेडशिफ्ट ($e(z)$) के साथ विकसित होता है, जबकि BHMF और ERDF के आकार स्थिर रहते हैं।
2. विकिरण-विनियमित ज्यामिति: RefleX रे-ट्रेसिंग कोड का उपयोग करके एक जटिल ज्यामितीय मॉडल बनाया गया है। इस मॉडल में शामिल हैं:
   • एक "लैंप-पोस्ट" एक्स-रे स्रोत (10 $r_g$ पर डबल पॉइंट सोर्स)।
   • एक अभिवृद्धि डिस्क (एक समतल वलय/annulus)।
   • एक ब्रॉड-लाइन रीजन (BLR) जिसे 0.2 के निश्चित कवरेज फैक्टर के साथ एक ठंडे, धूल-रहित वेज (wedge) के रूप में मॉडल किया गया है।
   • एक डस्टी टोरस जिसका आंतरिक किनारा BLR की बाहरी सीमा ($R_{\text{BLR,out}}$) पर सेट है। महत्वपूर्ण रूप से, टोरस का कवरेज फैक्टर ($C_F$) $\lambda_{\text{Edd}}$ के सिग्मॉइड फंक्शन (sigmoid function) के रूप में मॉडल किया गया है, जो अभिवृद्धि दर बढ़ने के साथ घटता है, जो विकिरण-विनियमित मॉडल के अनुरूप है।
   • उत्सर्जन स्पेक्ट्रा (हाई-एनर्जी कटऑफ के साथ पावर-लॉ) को $\lambda_{\text{Edd}}$ के साथ अनुभवजन्य संबंधों के आधार पर आवंटित किया गया है।
3. संख्यात्मक सिमुलेशन: RefleX का उपयोग करते हुए, लेखक एक्स-रे फोटॉन (0.1 keV – 1 MeV) को परिनाभिकीय सामग्री के माध्यम से प्रसारित करने के लिए मोंटे कार्लो सिमुलेशन करते हैं, जिसमें कॉम्पटन स्कैटरिंग, फोटोइलेक्ट्रिक अवशोषण और फ्लोरेसेंस (fluorescence) को ध्यान में रखा जाता है। $6 \times 10^6$ स्पेक्ट्रा ज्यामितीय और स्पेक्ट्रल मापदंडों के ग्रिड पर उत्पन्न किए जाते हैं।
4. अनुमान (Inference): मॉडल के मापदंडों को सीक्वेंशियल न्यूरल पोस्टीरियर एस्टीमेशन (SNPE) विधि का उपयोग करके अवलोकन संबंधी बाधाओं के विरुद्ध तुलना करके निकाला जाता है। बाधाओं में शामिल हैं:
   • कॉस्मिक एक्स-रे बैकग्राउंड (CXB) स्पेक्ट्रम (2–100 keV)।
   • तीन ऊर्जा बैंडों (14–195 keV, 8–24 keV, 2–10 keV) में विभेदक AGN संख्या गणना (number counts)।
   • स्थानीय AGNs का प्रेक्षित $N_H$ वितरण।
   • $\lambda_{\text{Edd}}$ के फलन के रूप में अवरुद्ध (obscured) AGNs का अंश।
   • सीमित फ्लक्स के फलन के रूप में प्रेक्षित CTK AGNs का अंश।

प्रमुख योगदान

• स्व-संगत मॉडलिंग: यह पहला जनसंख्या संश्लेषण अध्ययन है जो सीधे आंतरिक BHMF और ERDF से नमूना लेता है और एक एकीकृत ढांचे के भीतर एडिंगटन-अनुपात-निर्भर ज्यामिति और उत्सर्जन स्पेक्ट्रम को लागू करता है जहाँ अवरोध (obscuration) और प्रतिबिंब (reflection) स्थानिक वितरण के माध्यम से स्वाभाविक रूप से जुड़े हुए हैं।
• विकिरण-विनियमित सत्यापन: यह अध्ययन विकिरण-विनियमित एकीकरण मॉडल का कड़ाई से परीक्षण करता है, यह प्रदर्शित करते हुए कि एक ज्यामिति जहाँ टोरस का कवरेज फैक्टर $\lambda_{\text{Edd}}$ पर निर्भर करता है, वह CXB, संख्या गणना और अवशोषण गुणों को एक साथ पुनरुत्पादित कर सकता है।
• आंतरिक CTK अंश: लेखक $40 \pm 3\%$ का एक मजबूत आंतरिक CTK अंश प्राप्त करते हैं, जो पिछले अनुमानों से काफी अधिक है, जिसका श्रेय BLR अवरोध को दिया जाता है जो $\sim 20\%$ का एक बेसलाइन CTK अंश निर्धारित करता है।
• मापदंड नियंत्रण: अध्ययन प्रमुख भौतिक मापदंडों के लिए पोस्टीरियर वितरण प्रदान करता है, जिसमें टोरस कॉलम घनत्व वितरण ($\mu, \sigma$), $\lambda_{\text{Edd}}$ के साथ टोरस कवरेज फैक्टर की निर्भरता ($k$), और सक्रिय SMBH अंश का विकास ($p_1, z_c$) शामिल है।

परिणाम

• CXB और संख्या गणना: मॉडल CXB स्पेक्ट्रम की समग्र तीव्रता और आकार को सफलतापूर्वक पुनरुत्पादित करता है, हालांकि यह 60 keV से ऊपर के फ्लक्स का थोड़ा अधिक अनुमान लगाता है। विभेदक संख्या गणनाओं को अनिश्चितताओं के भीतर पुनरुत्पादित किया जाता है, जिसमें सॉफ्ट बैंड के कमजोर छोर पर और 14–195 keV बैंड के उज्ज्वल छोर पर मामूली अति-अनुमान (overpredictions) देखे गए हैं।
• अवशोषण गुण: मॉडल द्वि-मोडल (bimodal) $N_H$ वितरण और $\lambda_{\text{Edd}}$ के साथ अवरुद्ध अंश के प्रति-सहसंबंध को पुनरुत्पादित करता है। हालांकि, पता लगाए गए अवरुद्ध और अनब्लॉक (unobscured) AGNs की पूर्ण संख्या लगभग 15–25% अधिक अनुमानित है, जो संभवतः सिंथेटिक जनसंख्या में पूरी तरह से कैप्चर न किए गए चयन प्रभावों (selection effects) के कारण है।
• टोरस ज्यामिति: अनुमानित टोरस का आंतरिक किनारा $-3 < \log \lambda_{\text{Edd}} < 0.5$ के लिए $\sim 0.5 - 8$ pc तक फैला हुआ है, जो मिड-इन्फ्रारेड इंटरफेरोमेट्रिक आकार और विकिरण-विनियमित परिदृश्य के अनुरूप है। कवरेज फैक्टर $\lambda_{\text{Edd}}$ के बढ़ने के साथ घटता है, जो स्थानीय AGNs में देखे गए रुझान को पुनः प्राप्त करता है।
• प्रतिबिंब-अवरोध सहसंबंध: मॉडल गुणात्मक रूप से उस रुझान को पुनरुत्पादित करता है जहाँ अवरुद्ध स्रोतों के लिए प्रतिबिंब $N_H$ के साथ बढ़ता है, और अनब्लॉक AGNs के लिए प्रतिबिंब स्तरों की भविष्यवाणी करता है जो पिछले मॉडलों की तुलना में अवलोकनों के साथ बेहतर ढंग से मेल खाती है।
• विकास: सक्रिय SMBH अंश के अनुमानित विकास में $z \approx 1$ तक तीव्र वृद्धि और उसके बाद एक पठार (plateau) दिखाई देता है, जो XLF अध्ययनों में उपयोग किए जाने वाले प्योर डेंसिटी इवोल्यूशन मॉडल्स के अनुरूप है।

महत्व और दावे
यह शोध पत्र दावा करता है कि विकिरण-विनियमित एकीकरण मॉडल एक बहुत ही व्यापक सेट के एक्स-रे अवलोकनों के साथ सुसंगत है। घटनात्मक से भौतिक मॉडलिंग की ओर बढ़ते हुए, लेखक प्रदर्शित करते हैं कि सक्रिय SMBHs के सक्रिय अंश को नियंत्रित करने वाला एक सरल विकास नुस्खा, स्थानीय रूप से नमूना लिए गए BHMF और ERDF के साथ मिलकर CXB को पुनरुत्पादित करने के लिए पर्याप्त है। प्राप्त $40 \pm 3\%$ का आंतरिक CTK अंश स्थानीय ब्रह्मांड के हालिया NuSTAR अवलोकनों और पिछले जनसंख्या संश्लेषण भविष्यवाणियों के साथ अच्छी तरह से मेल खाता है। अध्ययन इस निष्कर्ष पर निकलता है कि अवरुद्ध और अनब्लॉक AGNs के बीच BHMF और ERDF के अंतर को एक एकल पैरेंट डिस्ट्रीब्यूशन से नमूना लेकर पुनरुत्पादित किया जा सकता है, जिसमें अंतर केवल अवरोध को नियंत्रित करने वाले भौतिक और ज्यामितीय गुणों से उत्पन्न होता है। लेखक सीमाओं को नोट करते हैं, जैसे कि सजातीय घटकों (homogeneous components) की धारणा और अन्य संभावित अवरोधकों (जैसे ध्रुवीय धूल, आउटफ्लो) की कमी, जो भविष्य के परिशोधन के लिए क्षेत्रों के रूप में सुझाव देते हैं।