Population synthesis of Galactic middle-aged pulsar wind nebulae I. Detection prospects for current and future instruments

समस्या विवरण
पल्सर विंड नेबुला (PWNe) गैलेक्टिक बहुत-उच्च-ऊर्जा (VHE; $E > 100$ GeV) $\gamma$-रे स्रोतों की सबसे बड़ी आबादी का प्रतिनिधित्व करते हैं और कण त्वरण (particle acceleration) तथा पल्सर-सुपरनोवा अवशेष (SNR) अंतःक्रियाओं के अध्ययन के लिए महत्वपूर्ण प्रयोगशाला के रूप में कार्य करते हैं। जबकि H.E.S.S., MAGIC, और VERITAS जैसे वर्तमान उपकरण कई PWNe की स्पष्ट पहचान कर चुके हैं, यथार्थवादी जनसंख्या-स्तरीय भविष्यवाणियों में ऐतिहासिक रूप से "रिवर्बरेशन चरण" (reverberation phase) के विस्तृत उपचार का अभाव रहा है। यह चरण तब होता है जब SNR का रिवर्स शॉक विस्फोट केंद्र की ओर वापस यात्रा करता है, जिससे PWN से टकराता है और उसे संकुचित (compress) करता है। यह घटना नेबुला की गतिशीलता को महत्वपूर्ण रूप से बदल देती है, चुंबकीय क्षेत्रों को प्रवर्धित करती है, कण आबादी को पुनर्जीवित करती है, और रेडिएटिव नुकसान को अत्यधिक बढ़ा देती है। पिछले मॉडलिंग दृष्टिकोणों में एक द्वैत (dichotomy) देखा गया था: सरल एक-ज़ोन थिन-शेल मॉडल गणनात्मक रूप से कुशल हैं लेकिन रिवर्बरेशन चरण की जटिल, गैर-रेखीय फीडबैक को पकड़ने में विफल रहते हैं, जबकि बहुआयामी हाइड्रोडायनामिक (HD) और मैग्नेटोहाइड्रोडायनामिक (MHD) सिमुलेशन सटीक हैं लेकिन $10^5$ वर्षों में हजारों स्रोतों के सांख्यिकीय अध्ययन के लिए अत्यधिक महंगे हैं। फलस्वरूप, गैलेक्टिक TeV PWN आबादी, विशेष रूप से आगामी वेधशालाओं के लिए, यदि इस चरण की उपेक्षा की जाती है या इसे सरल बनाया जाता है, तो भविष्यवाणियां पक्षपाती होने का जोखिम उठाती हैं।

कार्यप्रणाली
इन सीमाओं को संबोधित करने के लिए, लेखक हाइब्रिड TIDE+L फ्रेमवर्क का उपयोग करते हैं, जो मुक्त-विस्तार चरण (free-expansion phase) के लिए एक थिन-शेल डायनेमिकल मॉडल (TIDE) की गणनात्मक दक्षता को रिवर्बरेशन के दौरान SNR संरचना के लैग्रेंजियन (Lagrangian) उपचार के साथ जोड़ता है। यह दृष्टिकोण रिवर्बरेशन के जटिल संपीड़न प्रकरणों सहित सभी चरणों में PWN-SNR प्रणाली के स्व-संगत विकास की अनुमति देता है, जबकि रेडिएटिव और एडियाबेटिक नुकसान के अधीन कण स्पेक्ट्रम को ट्रैक करता है।

अध्ययन लगभग 2,400 गैलेक्टिक PWNe की एक सिंथेटिक आबादी का निर्माण करता है, जिसमें से सांख्यिकीय सुदृढ़ता सुनिश्चित करने के लिए 1,000 वास्तविकताओं (realizations) में 1,600 स्रोतों को यादृच्छिक रूप से चुना गया है। आबादी को निम्नलिखित वितरणों से लिया गया है:

• पल्सर गुण: प्रारंभिक स्पिन अवधि (गॉसियन वितरण, 10 ms पर ट्रंकेटेड), प्रारंभिक चुंबकीय क्षेत्र (लॉग-नॉर्मल), और ब्रेकिंग इंडेक्स ($n=2.33$)।
• SNR और पर्यावरणीय गुण: विस्फोट ऊर्जा (ट्रंकेटेड गॉसियन), इजेक्टा द्रव्यमान (GENEVA स्टेलर इवोल्यूशन मॉडल के माध्यम से ZAMS द्रव्यमान से व्युत्पन्न), और अंतरतारकीय माध्यम (ISM) घनत्व (यूनिफॉर्म वितरण)।
• स्पेक्ट्रल पैरामीटर: कण इंजेक्शन स्पेक्ट्रा, कंटेनमेंट कारक, और चुंबकीय अंश।

प्रत्येक सिस्टम को अधिकतम आयु $10^5$ वर्ष तक विकसित किया जाता है। इसके बाद प्राप्त $\gamma$-रे फ्लक्स की तुलना वर्तमान और भविष्य के उपकरणों: H.E.S.S., HAWC, चेरेनकोव टेलीस्कोप एरे ऑब्जर्वेटरी (CTAO), LHAASO, और सदर्न वाइड-फील्ड गामा-रे ऑब्जर्वेटरी (SWGO) की संवेदनशीलता सीमाओं और आकाश कवरेज के विरुद्ध की जाती है। विश्लेषण विस्तारित स्रोतों के लिए संवेदनशीलता को कम करने वाले स्रोत विस्तार प्रभावों (source extension effects) को भी ध्यान में रखता है।

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

1. यथार्थवादी रिवर्बरेशन मॉडलिंग: यह शोध पत्र एक सांख्यिकीय रूप से सुदृढ़ जनसंख्या संश्लेषण प्रस्तुत करता है जो रिवर्बरेशन के लैग्रेंजियन उपचार को स्पष्ट रूप से शामिल करता है, जो पिछले थिन-शेल मॉडलों के सरलीकृत विश्लेषणात्मक नुस्खों से आगे बढ़ता है।
2. हाइब्रिड फ्रेमवर्क अनुप्रयोग: यह एक बड़े पैमाने के जनसंख्या अध्ययन में TIDE+L फ्रेमवर्क के अनुप्रयोग को प्रदर्शित करता है, जो विस्तृत हाइड्रोडायनामिक भौतिकी और जनसंख्या संश्लेषण की गणनात्मक आवश्यकताओं के बीच के अंतर को पाटता है।
3. डिटेक्शन संभावनाएं: यह अध्ययन अगली पीढ़ी के उपकरणों के लिए, विशेष रूप से CTAO, SWGO, और LHAASO की प्राप्ति (yield) को मात्रात्मक रूप से निर्धारित करते हुए, विस्तृत गैलेक्टिक PWN आबादी के अनुमान प्रदान करने वाला पहला विवरण है।
4. विकासवादी चरण: यह जनसंख्या को मुक्त-विस्तार, संपीड़न (compressing), और पोस्ट-कंप्रेशन चरणों में वर्गीकृत करता है, और प्रत्येक के पता लगाने की क्षमता (detectability) और स्पेक्ट्रल विशेषताओं का विश्लेषण करता है।

परिणाम

• डिटेबिलिटी (पता लगाने की क्षमता): मॉडल भविष्यवाणी करता है कि CTAO वर्तमान में TeV रेंज में दृढ़ता से पहचाने गए PWNe की तुलना में एक क्रम अधिक PWNe का पता लगाएगा। विशेष रूप से, CTAO लगभग 172 स्रोतों (जनसंख्या का 10.7%) का पता लगाने की उम्मीद करता है, जो एक अपग्रेड किए गए "CTAO+" कॉन्फ़िगरेशन के साथ बढ़कर 226 (14.1%) हो जाता है। यह H.E.S.S. द्वारा अनुमानित ~52 स्रोतों की तुलना में काफी अधिक है।
• रिवर्बरेशन का प्रभाव: TIDE+L फ्रेमवर्क और सरलीकृत TIDE मॉडल के बीच तुलना से पता चलता है कि यथार्थवादी रिवर्बरेशन गतिशीलता की उपेक्षा करने से पता लगाने योग्य स्रोतों का कम अनुमान लगाया जाता है। सरल मॉडल संपीड़न और चुंबकीय क्षेत्र प्रवर्धन को अत्यधिक बढ़ाकर दिखाता है, जिसके परिणामस्वरूप अत्यधिक सिनक्रोट्रॉन नुकसान और एक क्षीण उच्च-ऊर्जा कण आबादी होती है, जिससे $\gamma$-रे फ्लक्स कम हो जाता है।
• विकासवादी वितरण: सिंथेटिक आबादी का अधिकांश हिस्सा (87.5%) पोस्ट-कंप्रेशन चरण में है। हालांकि, इस चरण के स्रोत आमतौर पर पुराने, बड़े और स्वाभाविक रूप से मंद होते हैं, जिसके परिणामस्वरूप मुक्त-विस्तार या संपीड़न चरणों की तुलना में कम डिटेक्टेबिलिटी अंश प्राप्त होता है।
• स्पेक्ट्रल विकास: व्यक्तिगत केस स्टडीज दर्शाती हैं कि कंप्रेशन फैक्टर (CF) रिवर्बरेशन के दौरान चमक प्रवर्धन और स्पेक्ट्रल हार्डनिंग/सॉफ्टनिंग की शक्ति और अवधि को नियंत्रित करता है। उच्च-CF वाले स्रोत स्पेक्ट्रल एनर्जी डिस्ट्रीब्यूशन (SED) में तीव्र, तीव्र परिवर्तन प्रदर्शित करते हैं, जबकि निम्न-CF वाले स्रोत मंद, विलंबित विकास दिखाते हैं।
• सहसंबंध (Correlations): अध्ययन एकीकृत 1–10 TeV ल्यूमिनोसिटी और पल्सर स्पिन-डाउन पावर के बीच एक सहसंबंध ($L_{1-10 \text{ TeV}} \propto \dot{E}^{1.15}$) पाता है। लेखक नोट करते हैं कि HGPS डेटा में देखे गए उथले ढलान ($L \propto \dot{E}^{0.59}$) को पुनः प्राप्त करने के लिए विशिष्ट चयन कट्स (जैसे, कम स्पिन-डाउन पावर या पुराने युग के स्रोतों को बाहर करना) लागू करने की आवश्यकता है, जो यह सुझाव देता है कि सर्वेक्षण पूर्वाग्रह और यंत्र संबंधी सीमाएं देखे गए सहसंबंधों को प्रभावित करती हैं।
• स्रोत विस्तार: सिम्युलेटेड कोणीय विस्तार आम तौर पर HGPS डेटा के अनुरूप हैं, हालांकि मॉडल HGPS में न्यूनतम देखे गए रेडियस से छोटे रेडियस वाले स्रोतों का उच्च अंश भविष्यवाणी करता है। यह गोलाकार समरूपता (spherical symmetry) धारणा और मिक्सिंग एवं पल्सर प्रॉपर मोशन की उपेक्षा के कारण है।

महत्व
यह शोध पत्र दावा करता है कि गैलेक्टिक TeV PWN आबादी की सटीक भविष्यवाणी करने के लिए रिवर्बरेशन चरण का यथार्थवादी मॉडलिंग आवश्यक है। SNR रिवर्स शॉक इंटरैक्शन की जटिल गतिशीलता को शामिल करके, यह अध्ययन वर्तमान स्रोत कैटलॉग की व्याख्या करने और आगामी सर्वेक्षणों की प्राप्ति का अनुमान लगाने के लिए एक अधिक विश्वसनीय सैद्धांतिक आधार प्रदान करता है। परिणाम पुष्टि करते हैं कि CTAO आगामी TeV जनसंख्या गणना में प्रमुख भूमिका निभाएगा, जो संभवतः HGPS में उच्च-फ्लक्स अज्ञात स्रोतों की प्रकृति को सुलझा सकता है। इसके अलावा, यह कार्य मध्य-आयु वाले PWNe और विस्तारित TeV हेलो के उद्भव के बीच एक प्राकृतिक विकासवादी लिंक का सुझाव देता है, यह देखते हुए कि भविष्य के उपकरणों द्वारा पहचाने गए पुराने सिस्टम शास्त्रीय कॉम्पैक्ट नेबुला के बजाय हेलो जैसी रूपात्मकता (morphologies) से मेल खाएंगे। यह अध्ययन एक श्रृंखला के पहले भाग के रूप में कार्य करता है, जो दृश्यता (visibility) के लिए सांख्यिकीय आधार स्थापित करता है, जिसके बाद एक आगामी पेपर रिवर्बरेशन के "सुपर-एफिशिएंसी" परिणामों पर केंद्रित होगा।