Bell Instability and Cosmic-Ray Acceleration in Active Galactic Nuclei Ultrafast Outflow Shocks

तकनीकी सारांश: सक्रिय गैलेक्टिक नाभिकों (AGN) के अल्ट्राफास्ट आउटफ्लो शॉक में बेल अस्थिरता (Bell Instability) और कॉस्मिक-रे त्वरण

समस्या विवरण
"नी" (knee) ऊर्जा ($\sim 10^{15}$ eV) के स्तर पर और उससे ऊपर कॉस्मिक किरणों (CRs) की उत्पत्ति का मुद्दा अभी भी अनसुलझा है, जिसमें सक्रिय गैलेक्टिक नाभिकों (AGN) को प्राथमिक उम्मीदवारों के रूप में प्रस्तावित किया गया है। विशेष रूप से, अल्ट्राफास्ट आउटफ्लो (UFOs)—जो कि AGN से निकलने वाली मध्यम सापेक्षतावादी पवनें ($\sim 0.1c-0.4c$) हैं—PeV–EeV ऊर्जाओं तक विसरित शॉक त्वरण (DSA) के लिए आशाजनक स्थल माने जाते हैं। पिछले अध्ययन अक्सर घटनात्मक मापदंडों (जैसे, $\epsilon_B$) पर निर्भर रहे हैं या चुंबकीय क्षेत्र प्रवर्धन के भौतिकी को स्व-सुसंगत रूप से मॉडल किए बिना मजबूत विक्षोभ सीमाओं को माना है। इसके अलावा, जबकि पार्टिकल-इन-सेल (PIC) सिमुलेशन ने कुशल चुंबकीय प्रवर्धन के माध्यम से गैर-अनुनादी हाइब्रिड (NRH या बेल) अस्थिरता का प्रदर्शन किया है, वे आमतौर पर कम अधिकतम ऊर्जाओं तक सीमित होते हैं जहाँ CR करंट उच्च होता है। उच्च ऊर्जाओं पर NRH अस्थिरता की दक्षता, जहाँ पलायन करने वाले CR करंट को $E_{\text{max}}$ के पास के कणों द्वारा ले जाया जाता है, अनिश्चित बनी हुई है। यह अध्ययन इस बात की जांच करता है कि क्या NRH अस्थिरता वास्तविक मापदंडों के तहत AGN UFO रिवर्स शॉक में PeV–EeV त्वरण को सक्षम करने के लिए चुंबकीय क्षेत्रों को स्व-सुसंगत रूप से प्रवर्धित कर सकती है।

कार्यप्रणाली
लेखक एक एक-आयामी मैग्नेटोहाइड्रोडायनामिक्स–कॉस्मिक रे (MHD–CR) संख्यात्मक ढांचे का उपयोग करते हैं, जो T. Inoue (2019) द्वारा विकसित कोड पर आधारित है। यह दृष्टिकोण CR वितरण फलन के लिए युग्मित टेलीग्राफ-प्रकार के डिफ्यूजन-कन्वेक्शन समीकरणों के साथ MHD समीकरणों को हल करता है। सेटअप की प्रमुख विशेषताओं में शामिल हैं:

• भौतिक मॉडल: सिमुलेशन एक UFO के अंतरतारकीय माध्यम (ISM) से टकराने वाले रिवर्स शॉक पर केंद्रित है। CR वितरण को समदैशिक ($f_0$) और अनिसोट्रोपिक ($f_1$) घटकों में विस्तारित किया गया है, जिसमें अनिसोट्रोपिक भाग उस CR करंट को संचालित करता है जो NRH अस्थिरता को उत्तेजित करता है।
• समीकरण: सिस्टम में फ्लुइड के लिए निरंतरता, संवेग, ऊर्जा और प्रेरण (induction) समीकरण शामिल हैं, जो इंजेक्शन, एडवेक्शन, डिफ्यूजन और $p\gamma$ ऊर्जा हानि (फोटोमेसॉन इंटरैक्शन) को ध्यान में रखते हुए CR के विकास के साथ जुड़े हुए हैं।
• मापदंड: अध्ययन व्यवस्थित रूप से बैकग्राउंड चुंबकीय क्षेत्र ($B_0$ $10^{-5}$ से $10^{-2}$ G), प्रारंभिक चुंबकीय विक्षोभ आयाम ($\xi_{B,\text{ini}}$), CR इंजेक्शन दर ($\eta$), और ISM घनत्व ($n_{\text{ISM}}$) को परिवर्तित करता है।
• रिज़ॉल्यूशन: NRH अस्थिरता की अधिकतम विकास तरंग दैर्ध्य और निम्न-ऊर्जा CR के डिफ्यूजन लेंथ को हल करने के लिए उच्च स्थानिक रिज़ॉल्यूशन बनाए रखा गया है, जबकि मोमेंटम स्पेस को $10^{19}$ eV तक लॉगरिदमिक रूप से विविक्त (discretized) किया गया है।

प्रमुख योगदान और परिणाम
पेपर बैकग्राउंड चुंबकीय क्षेत्र ($B_0$) की शक्ति के आधार पर त्वरण दक्षता और चुंबकीय व्यवहार में एक विशिष्ट संक्रमण की पहचान करता है:

1. कमजोर क्षेत्र शासन ($B_0 \lesssim 10^{-4}$ G):

   • इस शासन में, NRH अस्थिरता कुशलतापूर्वक कार्य करती है। पलायन करने वाला CR करंट अपस्ट्रीम चुंबकीय विक्षोभ के महत्वपूर्ण प्रवर्धन को संचालित करता है।
   • सिस्टम एक स्व-नियमित अवस्था तक पहुँच जाता है जहाँ अधिकतम CR ऊर्जा ($E_{\text{max}}$) और चुंबकीय ऊर्जा अंश ($\epsilon_B$) प्रारंभिक विक्षोभ आयाम ($\xi_{B,\text{ini}}$) से काफी हद तक स्वतंत्र होकर अभिसरित होते हैं।
   • $B_0 = 10^{-4}$ G और इंजेक्शन दर $\eta = 10^{-4}$ के लिए, $E_{\text{max}}$ लगभग $\sim 10^{16}$ eV तक पहुँच जाता है। $\eta$ को $10^{-3}$ तक बढ़ाने से $E_{\text{max}}$ में और वृद्धि होती है।
   • चुंबकीय प्रवर्धन CR करंट द्वारा संचालित होता है, जिसमें अस्थिरता का e-folding नंबर शॉक के अपस्ट्रीम में एक परिमित क्षेत्र में चरम पर होता है।

2. मजबूत क्षेत्र शासन ($B_0 \gtrsim 10^{-3}$ G):

   • जैसे-जैसे $B_0$ बढ़ता है, CR की जाइरो-त्रिज्या (gyroradius) कम होती जाती है, जिससे कणों के लिए बहुत दूर अपस्ट्रीम पलायन करना कठिन हो जाता है। परिणामस्वरूप, CR करंट घनत्व ($j_{\text{CR}}$) काफी गिर जाता है।
   • कम हुए करंट के कारण, NRH अस्थिरता के लिए विकास समय ($t_{\text{NRH}}$), एडवेक्शन समय ($t_{\text{adv}}$) से अधिक हो जाता है, जिससे अस्थिरता अप्रभावी हो जाती है।
   • इस शासन में, $E_{\text{max}}$ अब स्व-नियमित नहीं है बल्कि स्पष्ट रूप से प्रारंभिक स्थितियों ($B_0$ और $\xi_{B,\text{ini}}$) पर निर्भर है।
   • इसके अतिरिक्त, उच्च प्रारंभिक विक्षोभ ($\xi_{B,\text{ini}} = 0.1$) वाले मजबूत क्षेत्रों में, लघु-तरंगदैर्घ्य चुंबकीय उतार-चढ़ाव पैरामीट्रिक अस्थिरताओं (विशेष रूप से स्टिमुलेटेड ब्रिलुइन स्कैटरिंग) के माध्यम से तेजी से क्षय होते हैं, जिससे ऊर्जा ध्वनि तरंगों और बैकवर्ड-प्रोपगेटिंग अल्फवेन तरंगों में स्थानांतरित हो जाती है। यह डैम्पिंग (damping) स्कैटरिंग के लिए उपलब्ध विक्षोभ के आयाम को कम कर देती है, जो संभावित रूप से त्वरण दक्षता को कम कर सकती है जब तक कि लंबी-तरंगदैर्घ्य उतार-चढ़ाव हावी न हों।

3. ISM घनत्व और कूलिंग का प्रभाव:

   • उच्च ISM घनत्व ($n_{\text{ISM}}$) रिवर्स शॉक वेग ($v_{\text{sh}}$) को बढ़ाता है, जो CR करंट और NRH विकास दर को बढ़ाता है, जिससे उच्च $E_{\text{max}}$ प्राप्त होता है।
   • $B_0 \gtrsim 10^{-3}$ G के लिए $p\gamma$ कूलिंग एक सीमित कारक बन जाती है। $B_0 = 10^{-2}$ G के लिए, कूलिंग $E_{\text{max}}$ को $\sim 10^{18}$ eV से घटाकर सब-EeV शासन ($\sim 3.5 \times 10^{17}$ eV) तक ले आती है।

महत्व और दावे
लेखक दावा करते हैं कि EeV-स्केल त्वरण प्राप्त करना सरल एक्सट्रपलेशन या निम्न-ऊर्जा PIC सिमुलेशन द्वारा सुझाए गए अनुमानों की तुलना में अधिक बाधित है। अध्ययन दर्शाता है कि:

• NRH अस्थिरता सार्वभौमिक रूप से कुशल नहीं है; इसकी प्रभावशीलता उच्च चुंबकीय क्षेत्रों में बाधित होती है जिनकी आवश्यकता EeV त्वरण के लिए होती है क्योंकि आवश्यक पलायन करने वाला CR करंट अस्थिरता को चलाने के लिए बहुत कमजोर होता है।
• कमजोर-क्षेत्र शासन में, NRH अस्थिरता $E_{\text{max}}$ निर्धारित करने के लिए एक स्व-सुसंगत तंत्र प्रदान करती है, लेकिन परिणामी ऊर्जा ($\sim 10^{15}-10^{16}$ eV) EeV स्केल से कम रह जाती है।
• $\sim 10^{18}$ eV तक पहुँचने के लिए, सिस्टम को मजबूत प्रारंभिक स्थितियों ($B_0 \ge 10^{-2}$ G और $\xi_{B,\text{ini}} \ge 0.1$) पर निर्भर रहना चाहिए जहाँ NRH अस्थिरता अप्रभावी है, और त्वरण प्रारंभिक विक्षोभ द्वारा शासित है। हालांकि, इसके लिए यह आवश्यक है कि लघु-तरंगदैर्घ्य उतार-चढ़ाव पैरामीट्रिक अस्थिरताओं के माध्यम से तेजी से क्षय न हों और $p\gamma$ कूलिंग अक्षम बनी रहे।
• परिणाम बताते हैं कि यथार्थवादी सैद्धांतिक मॉडलों को NRH प्रवर्धन को एक निश्चित प्रवर्धन कारक मानने के बजाय $E_{\text{max}}$ के एक फलन (function) के रूप में मानना चाहिए, क्योंकि अस्थिरता को चलाने वाला करंट अधिकतम ऊर्जा बढ़ने के साथ घटता जाता है।

पेपर निष्कर्ष निकालता है कि हालांकि UFO रिवर्स शॉक्स आशाजनक उम्मीदवार हैं, लेकिन EeV त्वरण (मजबूत क्षेत्र, लंबी-तरंगदैर्घ्य विक्षोभ, कमजोर फोटॉन क्षेत्र) के लिए आवश्यक स्थितियों की एक साथ संतुष्टि करना गैर-तुच्छ (non-trivial) है और इसके लिए आगे के अवलोकन संबंधी और सैद्धांतिक अन्वेषण की आवश्यकता है, विशेष रूप से UFO की क्लंपी प्रकृति और उच्च द्रव्यमान बहिर्वाह दरों के संबंध में।