Proton Acceleration by Collisionless Shocks in Supermassive Black Hole Coronae: Implications for High-Energy Neutrinos

तकनीकी सारांश: सुपरमैसिव ब्लैक होल कोरोना में टकराव रहित झटकों (Collisionless Shocks) द्वारा प्रोटॉन त्वरण

समस्या विवरण
आइसक्यूब न्यूट्रिनो ऑब्जर्वेटरी (IceCube Neutrino Observatory) के हालिया अवलोकनों ने पास के रेडियो-शांत सक्रिय गैलेक्टिक नाभिकों (AGNs), विशेष रूप से सिफर्ट गैलेक्सी NGC 1068 से उच्च-ऊर्जा न्यूट्रिनो की एक महत्वपूर्ण अधिकता का पता लगाया है। ये अवलोकन एक विसंगति प्रस्तुत करते हैं: जबकि उच्च-ऊर्जा न्यूट्रिनो का पता लगाया जाता है, उच्च-ऊर्जा गामा-किरणों का संगत फ्लक्स मानक हैड्रोनिक इंटरैक्शन मॉडल द्वारा अनुमानित स्तर से काफी कम है। यह सुझाव देता है कि न्यूट्रिनो एक ऐसे क्षेत्र में उत्पन्न होते हैं जो गामा-किरणों के लिए अपारदर्शी (opaque) है, जो संभवतः केंद्रीय सुपरमैसिव ब्लैक होल (SMBH) के चारों ओर का गर्म कोरोना है, जहाँ एक्स-रे फोटॉन के साथ $\gamma\gamma$ युग्म उत्पादन (pair production) के माध्यम से गामा-किरणों का क्षीणन (attenuation) होता है।

अवलोकित न्यूट्रिनो स्पेक्ट्रम की व्याख्या करने के लिए, जनक प्रोटॉन को कोरोना के भीतर $\sim 100$ TeV की ऊर्जा तक त्वरित होना चाहिए। हालांकि, SMBH कोरोना में प्रमुख त्वरण तंत्र अभी भी अनिश्चित है। जबकि डिफ्यूसिव शॉक एक्सीलरेशन (DSA) एक अग्रणी उम्मीदवार है, कोरोनल स्थितियों के तहत इसकी व्यवहार्यता अनिश्चित है। विशेष रूप से, कोरोनल शॉक्स में निम्न सोनिक मच संख्या ($M_s \sim 2-10$) और ट्रांस-रिलेटिविस्टिक वेग होने की उम्मीद है, जो पारंपरिक रूप से अक्षम त्वरण वाले शासन (regimes) माने जाते हैं। इसके अलावा, दो-तापमान वाले प्लाज्मा ($T_i \gg T_e$) का व्यवहार और इन वातावरणों में प्रोटॉन बनाम इलेक्ट्रॉन त्वरण की विशिष्ट दक्षता के लिए प्रथम-सिद्धांतों (first-principles) के बाधाओं का अभाव है।

कार्यप्रणाली
इन अनिश्चितताओं को दूर करने के लिए, लेखक SMILEI कोड का उपयोग करके प्रथम-सिद्धांतों पर आधारित, एक-आयामी पार्टिकल-इन-सेल (PIC) सिमुलेशन का उपयोग करते हैं। अध्ययन का केंद्र क्वैसी-पैरेलल कोलिजनलेस शॉक्स ($\theta = 30^\circ$) है, जिन्हें SMBH कोरोना के भौतिक मापदंडों के अनुरूप ट्यून किया गया है।

सिमुलेशन सर्वेक्षण एक व्यवस्थित रूप से व्यापक पैरामीटर स्पेस का अन्वेषण करता है जो सिफर्ट गैलेक्सी के लिए प्रेक्षित बाधाओं से प्रेरित है:

• सोनिक मच संख्या ($M_s$): ट्रांस-सोनिक ($M_s \approx 2$) से लेकर उच्च मान ($M_s \approx 30$) तक।
• अल्वेनिक मच संख्या ($M_A$): $4$से$46$ तक।
• शॉक वेग ($v_{sh}$): ट्रांस-रिलेटिविस्टिक गति ($v_{sh}/c \approx 0.17 - 0.66$)।
• आयन-टू-इलेक्ट्रॉन तापमान अनुपात ($T_i/T_e$): $0.01$से$50$ तक, जो विच्छेदित (decoupled) थर्मल अवस्थाओं को दर्शाता है जो कोलिजनलेस कोरोना में अपेक्षित है।

सिमुलेशन में एक कम किया गया आयन-टू-इलेक्ट्रॉन द्रव्यमान अनुपात ($m_i/m_e = 100$) का उपयोग किया गया है ताकि मुख्य आयामहीन मापदंडों को सुरक्षित रखते हुए गणनात्मक व्यवहार्यता सुनिश्चित की जा सके। रन को पर्याप्त समय ($t \gtrsim 400 \Omega_{ci}^{-1}$) के लिए विकसित किया गया है ताकि इलेक्ट्रोमैग्नेटिक अस्थिरताओं का संतृप्ति (saturation) और नॉनथर्मल पावर-लॉ टेल्स का विकास हो सके।

मुख्य परिणाम

1. मजबूत प्रोटॉन त्वरण: अध्ययन पाता है कि DSA, SMBH कोरोना में प्रोटॉन त्वरण के लिए एक मजबूत और कुशल तंत्र है। निम्न सोनिक मच संख्याओं के बावजूद, शॉक्स लगातार लगभग 10% शॉक की गतिज ऊर्जा को नॉनथर्मल आयनों में परिवर्तित करते हैं। यह दक्षता पूरे पैरामीटर सर्वेक्षण में स्थिर रहती है, जिसमें $M_s \approx 2$ वाले मामले भी शामिल हैं, जो इस धारणा को चुनौती देती है कि निम्न-$M_s$ वाले शॉक्स अप्रभावी त्वरक होते हैं।
2. त्वरण गतिकी (Acceleration Dynamics): त्वरित प्रोटॉन की अधिकतम ऊर्जा शॉक संरचना के क्वैसी-स्टेडी स्टेट तक पहुँचने के बाद समय के साथ रैखिक रूप से बढ़ती है। त्वरण दर DSA सिद्धांत द्वारा अनुमानित स्केलिंग, $t_{acc} \propto \gamma \beta_{sh}^{-2} \Omega_{ci}^{-1}$ का अनुसरण करती है। लेखक एक सामान्यीकृत त्वरण गुणांक, $C^{-1} \equiv d(\gamma \beta_{sh}^{-2})/d(t\Omega_{ci})$, प्राप्त करते हैं, जो एक सार्वभौमिक मान $C^{-1} \approx 0.1$ प्रदर्शित करता है, जो सभी परीक्षण किए गए मापदंडों में दोगुने के कारक से भी कम भिन्न होता है।
3. इलेक्ट्रॉन त्वरण की अक्षमता: प्रोटॉन के विपरीत, इलेक्ट्रॉन त्वरण अत्यधिक अक्षम है, जहाँ नॉनथर्मल इलेक्ट्रॉन ऊर्जा अंश 1% से नीचे (आमतौर पर कोरोनल-जैसे $T_i/T_e$ अनुपातों के लिए $\sim 10^{-3}$) रहता है। इलेक्ट्रॉन त्वरण की दक्षता $T_i/T_e$ अनुपात बढ़ने के साथ घटती जाती है। विशिष्ट ट्रांस-रिलेटिविस्टिक मामलों में, जहाँ चुंबकीय क्षेत्र प्रवर्धन (amplification) मजबूत है, इलेक्ट्रॉन पावर-लॉ टेल उच्च ऊर्जा तक विस्तारित नहीं हो पाती है, जो यह सुझाव देता है कि जब प्रवर्धित क्षेत्र अंतराल पर शॉक को इलेक्ट्रॉनों के लिए प्रभावी रूप से सुपरल्यूमिनल बना देते हैं, तो निरंतर DSA बाधित हो जाता है।
4. अस्थिरता तंत्र (Instability Mechanisms): स्कैटरिंग के लिए आवश्यक अपस्ट्रीम चुंबकीय विक्षोभ (turbulence) अधिकांश पैरामीटर स्पेस में मुख्य रूप से नॉन-रेजोनेंट इंस्टेबिलिटी (NRI, या बेल की अस्थिरता) द्वारा संचालित होता है। रेजोनेंट इंस्टेबिलिटी (RI) केवल बहुत कम अल्वेनिक मच संख्या ($M_A \lesssim 8$) पर प्रभावी होती है।

NGC 1068 और न्यूट्रिनो उत्पादन के लिए निहितार्थ
NGC 1068 के लिए इन परिणामों को लागू करते हुए, लेखक प्रदर्शित करते हैं कि:

• समय-सीमा (Timescales): प्रोटॉन प्रासंगिक कूलिंग (बेथे-हीटले, फोटोमेसन, $pp$) और गतिशील (फ्री-फॉल) समय-सीमाओं से कम समय में$\sim 100$ TeV तक त्वरित हो सकते हैं।
• ऊर्जा बजट: देखे गए न्यूट्रिनो चमक (luminosity) की व्याख्या करने के लिए आवश्यक नॉनथर्मल प्रोटॉन शक्ति कुल शॉक शक्ति का $\sim 3\%$ अनुमानित है। यह सिमुलेशन में मापी गई $\sim 10\%$ त्वरण दक्षता के भीतर है।
• गामा-रे घाटा: इलेक्ट्रॉन त्वरण की निम्न दक्षता और गामा-किरणों के प्रति कोरोना की उच्च अपारदर्शिता स्वाभाविक रूप से देखी गई उच्च-ऊर्जा गामा-किरणों की कमी की व्याख्या करती है, जबकि न्यूट्रिनो को बाहर निकलने की अनुमति देती है।

महत्ता और दावे
यह शोध पत्र दावा करता है कि यह SMBH कोरोना के विशिष्ट संदर्भ में DSA पर पहला प्रथम-सिद्धांत आधारित, स्व-सुसंगत (self-consistent) प्रतिबंध प्रदान करता है। इसकी प्राथमिक महत्ता यह स्थापित करने में है कि कोलिजनलेस शॉक्स, निम्न सोनिक मच संख्या और उच्च आयन-टू-इलेक्ट्रॉन तापमान अनुपात जैसी चुनौतीपूर्ण स्थितियों के तहत भी, उच्च-ऊर्जा प्रोटॉन उत्पन्न करने के लिए एक व्यवहार्य और कुशल तंत्र हैं।

लेखक स्पष्ट रूप से कहते हैं कि जबकि उनके परिणाम न्यूट्रिनो उत्पादन के लिए हैड्रोनिक मॉडल का पुरजोर समर्थन करते हैं, एक जटिल खगोलीय वस्तु जैसे कि AGN कोरोना पर स्थानीय गतिज सिमुलेशन परिणामों को लागू करने के लिए सावधानीपूर्वक विचार की आवश्यकता है। ऊर्जा बजट विश्लेषण को सूक्ष्म भौतिकी (microphysics) पर आधारित एक ऑर्डर-ऑफ-मैग्निट्यूड मूल्यांकन के रूप में प्रस्तुत किया गया है, जो सूक्ष्म भौतिकी और स्थूल उच्च-ऊर्जा घटनाओं के बीच के अंतर को पाटता है। अध्ययन यह दावा नहीं करता है कि यह अधिकतम प्रोटॉन ऊर्जा या 3D में अस्थिरताओं के विस्तृत संक्रमण से जुड़ी सभी अनिश्चितताओं को हल कर देता है, यह नोट करते हुए कि इनके लिए आगे के अन्वेषण की आवश्यकता है।