🔭 astrophysics

Implementation of CR Energy SPectrum (CRESP) algorithm in PIERNIK MHD code. II. Propagation of Primary and Secondary nuclei in a magneto-hydrodynamical environment

यह शोध पत्र प्राथमिक और द्वितीयक ब्रह्मांडीय किरण नाभिकों के युग्मित उत्पादन और प्रसार को सिम्युलेट करने के लिए PIERNIK MHD कोड के CRESP मॉड्यूल का एक विस्तार प्रस्तुत करता है, जो यह प्रकट करता है कि चुंबकीय उत्प्लावकता प्रभावों और कम कण निवास समय के कारण विसरण गुणांकों में वृद्धि क्रमशः द्वितीयक-से-प्राथमिक फ्लक्स अनुपात को बढ़ाती है और सुपरनोवा दरों में वृद्धि इसे कम करती है।

मूल लेखक: Antoine Baldacchino-Jordan, Michał Hanasz, Mateusz Ogrodnik, Dominik Wóltański, Artur Gawryszczak, Andrew W. Strong, Philipp Girichidis

प्रकाशित 2026-03-16

📖 5 मिनट में पढ़ें🧠 गहराई से पढ़ें

CC BY 4.0

मूल लेखक: Antoine Baldacchino-Jordan, Michał Hanasz, Mateusz Ogrodnik, Dominik Wóltański, Artur Gawryszczak, Andrew W. Strong, Philipp Girichidis

मूल पेपर CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/) के तहत लाइसेंस किया गया है। ✨ नीचे दिए गए पेपर की यह व्याख्या AI से तैयार की गई है। इसे लेखकों ने न तो लिखा है, न इसका समर्थन किया है। तकनीकी सटीकता के लिए मूल पेपर देखें। पूरा डिस्क्लेमर पढ़ें

कल्पना कीजिए कि ब्रह्मांड एक शांत शून्य नहीं, बल्कि एक हलचल भरी, अराजक नगरी है जहाँ अदृश्य, उच्च गति वाले कण जिन्हें कॉस्मिक रेज़ (Cosmic Rays) कहा जाता है, अति-सक्रिय यात्रियों की तरह इधर-उधर दौड़ रहे हैं। ये केवल साधारण यात्री नहीं हैं; ये फटे हुए सितारों के "भूत" हैं, जो उस ऊर्जा को ले जाते हैं जो उस हवा (इंटरस्टेलर मीडियम) को आकार देती है जिससे होकर वे गुजरते हैं।

यह शोध पत्र इन कॉस्मिक यात्रियों के लिए एक अति-परिष्कृत ट्रैफिक सिम्युलेटर बनाने के बारे में है। लेखक, PIERNIK नामक एक शक्तिशाली कंप्यूटर कोड का उपयोग करके, अपने मॉडल को न केवल मुख्य यात्रियों (प्राइमरी पार्टिकल्स जैसे कार्बन और ऑक्सीजन) को ट्रैक करने के लिए, बल्कि उनके द्वारा रास्ते में बनाए गए "आकस्मिक" यात्रियों (सेकेंडरी पार्टिकल्स जैसे लिथियम, बेरिलियम और बोरॉन) को भी ट्रैक करने के लिए अपग्रेड कर रहे हैं।

यहाँ उनकी खोज का विवरण दिया गया, जिसे रोजमर्रा के उपमाओं के माध्यम से समझाया गया है:

1. सेटअप: एक बॉक्स में ब्रह्सीय नगरी

शोधकर्ताओं ने हमारे सूर्य के पास हमारी आकाशगंगा के एक हिस्से का प्रतिनिधित्व करने वाला एक डिजिटल "बॉक्स" बनाया। इस बॉक्स के भीतर, उन्होंने सिम्युलेट किया:

गैस: शहर की हवा।
चुंबकीय क्षेत्र (Magnetic Fields): अदृश्य रेल की पटरियाँ जो कणों का मार्गदर्शन करती हैं।
सुपरनोवा (Supernovae): यादृच्छिक विस्फोट (जैसे आतिशबाजी) जो सिस्टम में नए कॉस्मिक रेज़ इंजेक्ट करते हैं।

अपने कोड के पिछले संस्करणों में, वे केवल "मुख्य" कणों को ट्रैक कर सकते थे। इस नए संस्करण में, उन्होंने यह ट्रैक करने की क्षमता भी जोड़ दी है कि क्या होता है जब एक तेज़ गति वाला प्राइमरी पार्टिकल गैस के परमाणु से टकराता है। यह एक बिलियर्ड बॉल के क्लस्टर में टकराने जैसा है; मूल गेंद धीमी हो जाती है, और नई, छोटी गेंदें अलग-अलग दिशाओं में बिखर जाती हैं।

2. बड़ी हैरानी: "लीकी बॉक्स" बनाम "जीवंत शहर"

द दशकों से, वैज्ञानिक एक सरल मॉडल का उपयोग करते आए हैं जिसे "लीकी बॉक्स" (Leaky Box) कहा जाता है। कल्पना कीजिए कि एक बाल्टी है जिसके नीचे एक छेद है।

पुरानी तर्कशक्ति: यदि आप छेद को बड़ा करते हैं (तेज़ डिफ्यूजन), तो कण तेज़ी से बाहर निकल जाते हैं। इसलिए, आपके पास अंदर कम कण होते हैं, और कम टकराव होते हैं, जिसका अर्थ है कम "आकस्मिक" सेकेंडरी कण।
नई खोज: लेखकों ने पाया कि उनके जीवंत, सांस लेते सिम्युलेशन में, इसके विपरीत होता है!
- उपमा: कल्पना कीजिए कि कॉस्मिक रेज़ गर्म हवा के गुब्बारे हैं। जब वे बहुत तेज़ चलते हैं (उच्च डिफ्यूजन), तो वे केवल बाहर नहीं निकलते; वे वास्तव में अपने आस-पास की हवा को धक्का देते हैं। यह एक "उत्प्लावकता" (buoyancy) प्रभाव पैदा करता है जो चुंबकीय "रेल की पटरियों" के आकार को बदल देता है।
- परिणाम: जब डिफ्यूजन अधिक होता है, तो चुंबकीय पटरियाँ अधिक क्षैज़ (horizontal) और कम ऊर्ध्वाधर (vertical) हो जाती हैं। यह प्राइमरी कणों को "जमीन" (गैलेक्टिक डिस्क) के करीब अधिक समय तक फंसा कर रखता है क्योंकि वे आसानी से ऊपर नहीं चढ़ सकते।
- परिणामस्वरूप: क्योंकि प्राइमरी कण भीड़भाड़ वाली डिस्क में अधिक समय तक रहते हैं, वे गैस परमाणुओं से अधिक बार टकराते हैं, जिससे अधिक सेकेंडरी कण बनते हैं। इसलिए, तेज़ डिफ्यूजन = अधिक सेकेंडरी कण। यह पुरानी थ्योरी को पूरी तरह उलट देता है!

3. सुपरनोवा कारक: "हवा" का प्रभाव

टीम ने यह भी परीक्षण किया कि विस्फोटों (सुपरनोवा) की संख्या बदलने पर क्या होता है।

अधिक विस्फोट: यह आकाशगंगा से ऊपर की ओर बहने वाली एक मजबूत "हवा" पैदा करता है। यह एक विशाल पंखे के चालू होने जैसा है।
प्रभाव: यह हवा कॉस्मिक रेज़ को डिस्क से बहुत जल्दी बाहर धकेल देती है। वे गैस परमाणुओं से टकराकर सेकेंडरी कण बनाने के लिए पर्याप्त समय तक नहीं रुक पाते।
परिणाम: अधिक विस्फोट वास्तव में कम सेकेंडरी कणों (एक कम बोरॉन-टू-कार्बन अनुपात) की ओर ले जाते हैं क्योंकि "यात्री" शहर से इतनी जल्दी बाहर निकाल दिए जाते हैं कि वे दुर्घटनाएं करने का समय ही नहीं पाते।

4. यह क्यों महत्वपूर्ण है

वैज्ञानिकों ने अपने सिम्युलेशन परिणामों की तुलना अंतर्राष्ट्रीय अंतरिक्ष स्टेशन पर AMS-02 उपकरण द्वारा एकत्र किए गए वास्तविक डेटा से की।

उन्होंने पाया कि उनका जटिल, स्व-सुसंगत मॉडल (जहाँ कण और गैस एक-दूसरे से बात करते हैं) पुराने, स्थिर मॉडलों की तुलना में अलग व्यवहार करता है।
उन्होंने पाया कि वास्तविक ब्रह्मांड से मेल खाने के लिए, आप केवल एक संख्या (जैसे डिफ्यूजन दर) को बदलकर काम नहीं चला सकते। आपको यह समझना होगा कि कैसे हवा, चुंबकीय क्षेत्र और विस्फोट एक साथ नृत्य करते हैं।

निष्कर्ष (The Takeaway)

यह शोध पत्र एक बड़ी प्रगति है क्योंकि यह आकाशगंगा को एक स्थिर कमरे के रूप में मानना बंद करता है और इसे एक गतिशील पारिस्थितिकी तंत्र (dynamic ecosystem) के रूप में मानना शुरू करता है।

पुराना दृष्टिकोण: कॉस्मिक रेज़ केवल हवा में उड़ती धूल हैं।
नया दृष्टिकोण: कॉस्मिक रेज़ स्वयं हवा हैं। वे गैस को धकेलते हैं, चुंबकीय क्षेत्रों को मोड़ते हैं, और ऐसा करके अपना रास्ता खुद बदलते हैं।

इस जटिल नृत्य को समझकर, वैज्ञानिक अंततः यह पता लगा सकते हैं कि ये कण हमारे आकाशगंगा में कितने समय तक रहते हैं और वे कहाँ से आते हैं, जिसमें वे "आकस्मिक" कणों (जैसे बोरॉन) का उपयोग रहस्य को सुलझाने के लिए सुराग के रूप में करते हैं।

यहाँ PIERNIK MHD कोड में CR ऊर्जा स्पेक्ट्रम (CRESP) एल्गोरिदम के कार्यान्वयन पर आधारित शोध पत्र का विस्तृत तकनीकी सारांश दिया गया है: "Implementation of CR Energy SPectrum (CRESP) algorithm in PIERNIK MHD code. II. Propagation of Primary and Secondary nuclei in a magneto-hydrodynamical environment."

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

कॉस्मिक रे (CR) प्रसार को पारंपरिक रूप से "लीकी-बॉक्स" (leaky-box) मॉडलों या स्थिर विसरण-संवहन समीकरणों (जैसे GALPROP) का उपयोग करके मॉडल किया जाता है जो एक स्थिर अंतरतारकीय माध्यम (ISM) मान लेते हैं। हालाँकि, ये मॉडल अक्सर ISM पर CR के गतिशील फीडबैक, जैसे कि CR-संचालित पवनें (winds) और चुंबकीय क्षेत्र का पुनर्गठन, की उपेक्षा करते हैं। जबकि स्व-सुसंगत मैग्नेटो-हाइड्रोडायनामिकल (MHD) मॉडल मौजूद हैं, वे अक्सर CR को एक स्पेक्ट्रली अनरिज़ॉल्व्ड "ग्रे" (grey) फ्लूइड के रूप में देखते हैं या केवल इलेक्ट्रॉनों पर ध्यान केंद्रित करते हैं।

एक स्व-सुसंगत, स्पेक्ट्रली रिज़ॉल्व्ड मॉडल की आवश्यकता है जो एक साथ कई प्राथमिक और द्वितीयक CR नाभिक प्रजातियों के प्रसार को ट्रैक कर सके और उन्हें विकसित होते तापीय गैस, चुंबकीय क्षेत्रों और सुपरनोवा (SN) फीडबैक के साथ गतिशील रूप से जोड़ सके। विशेष रूप से, लेखक प्राथमिक नाभिकों (C, N, O) के स्पलेशन (spallation) से उत्पन्न होने वाले द्वितीयक नाभिकों (Li, Be, B) के उत्पादन और प्रसार को एक वास्तविक, गतिशील गैलेक्टिक वातावरण के भीतर मान्य करने का लक्ष्य रखते हैं।

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

संख्यात्मक ढांचा: PIERNIK MHD कोड

यह अध्ययन PIERNIK का उपयोग करता है, जो एक ग्रिड-आधारित, मल्टीफ्लुइड MHD कोड है। लेखक मौजूदा CRESP (Cosmic Ray Energy SPectrum) एल्गोरिदम का विस्तार करते हैं, जिसका उपयोग पहले इलेक्ट्रॉनों के लिए किया गया था, ताकि इसमें हैड्रोनिक CR प्रजातियों को शामिल किया जा सके।

टू-मोमेंट दृष्टिकोण (Two-Moment Approach): कोड मोमेंट बिन्स के भीतर संख्या घनत्व ( $n_a$ ) और गतिज ऊर्जा घनत्व ( $e_a$ ) का उपयोग करके वितरण फलनों (distribution functions) के लिए फॉकर-प्लैंक समीकरण को हल करता है।
पीस-वाइज पावर-लॉ सन्निकटन (Piece-wise Power-Law Approximation): प्रत्येक मोमेंट बिन के भीतर, वितरण फलन को एक पावर लॉ के रूप में अनुमानित किया जाता है। यह प्रत्येक चरण में पूर्ण वितरण फलन को पूरी तरह से हल किए बिना ऊर्जा लाभ/हानि पदों के कुशल संख्यात्मक एकीकरण की अनुमति देता है।
स्पेक्ट्रल रिज़ॉल्यूशन: मॉडल गैर-सापेक्षिक (non-relativistic) से अति-सापेक्षिक (ultra-relativistic) ऊर्जाओं ( $5 \times 10^{-1}$ से $10^5$ प्रोटॉन द्रव्यमान) तक 16 मोमेंट बिन्स को ट्रैक करता है।

कार्यान्वित भौतिक प्रक्रियाएं

प्राथमिक नाभिक: $^{12}\text{C}$ , $^{14}\text{N}$ , और $^{16}\text{O}$ को सुपरनोवा अवशेषों (SNRs) के माध्यम से इंजेक्ट किया जाता है। प्रोटॉन को ISM की गतिशीलता को संचालित करने के लिए एक "ग्रे" (स्पेक्ट्रली अनरिज़ॉल्वड) रिलेटिविस्टिक फ्लूइड के रूप में माना जाता है।
द्वितीयक नाभिक: द्वितीयक प्रजातियों ( $^{7}\text{Li}$ , $^{9}\text{Be}$ , $^{10}\text{Be}$ , $^{10}\text{B}$ , $^{11}\text{B}$ ) को प्राथमिक नाभिकों के तापीय ISM प्रोटॉनों के विरुद्ध स्पलेशन (हैड्रोनिक टकराव) के माध्यम से उत्पन्न किया जाता है। क्रॉस-सेक्शन GALPROP से लिए गए हैं।
रेडियोधर्मी क्षय (Radioactive Decay): अस्थिर आइसोटोप (विशेष रूप से $^{10}\text{Be}$ ) के क्षय को शामिल किया गया है।
ऊर्जा हानि: मॉडल में एडियाबेटिक विस्तार हानि और कूलिंग (Bethe-Block फॉर्मूला का उपयोग करके) शामिल है। स्ट्रीमिंग प्रभावों को वर्तमान में छोड़ दिया गया है।
परिवहन तंत्र: CR तापीय गैस के साथ एडवेक्शन (advection) और चुंबकीय क्षेत्र रेखाओं के अनुदिश एनिसोट्रोपिक, रिजिडिटी-डिपेंडेंट डिफ्यूजन ( $D_{\parallel} \propto R^{\delta}$ ) का अनुभव करते हैं।

सिमुलेशन सेटअप

डोमेन: एक 3D गुरुत्वाकर्षण-स्तरीकृत बॉक्स ( $0.5 \times 1 \times 8$ kpc) जो सौर मंडल के पास के स्थानीय ISM का प्रतिनिधित्व करता है।
प्रारंभिक स्थितियाँ: $n_{\text{ISM}} \approx 1 \text{ cm}^{-3}$ और $T \approx 7000$ K के साथ हाइड्रोस्टैटिक संतुलन।
फीडबैक: यादृच्छिक रूप से वितरित सुपरनोवा (SN) प्रोटॉन इंजेक्ट करते हैं (प्रत्येक SN से $10^{51}$ erg का 10%) ताकि सिस्टम को संचालित किया जा सके।
पैरामीटर अध्ययन: लेखकों ने तीन प्रमुख मापदंडों को बदलकर सिमुलेशन किया:
1. SN दर: 20, 60, और 80 kpc $^{-2}$ Myr $^{-1}$ ।
2. डिफ्यूजन गुणांक ( $D_{\parallel}$ ): $3 \times 10^{28}$ और $6 \times 10^{28}$ cm $^2$ s $^{-1}$ ।
3. रिजिडिटी इंडेक्स ( $\delta$ ): 0.3 और 0.5।

3. मुख्य योगदान (Key Contributions)

हैड्रोनिक प्रजातियों के लिए CRESP का विस्तार: यह हैड्रोनिक प्रजातियों को एक साथ संभालने के लिए PIERNIK में CRESP एल्गोरिदम का पहला कार्यान्वयन है, जिसमें स्पलेशन सोर्स टर्म और रेडियोधर्मी क्षय शामिल है।
स्व-सुसंगत CR-ISM कपलिंग: स्थिर मॉडलों के विपरीत, यह अध्ययन प्रदर्शित करता है कि कैसे CR गतिशील रूप से ISM को आकार देते हैं (बॉयन्सी और विंड्स के माध्यम से) और कैसे वह बदला हुआ वातावरण बदले में CR परिवहन को प्रभावित करता है।
B/C अनुपात पर नवीन खोज: यह शोध पत्र चुंबकीय क्षेत्र के पुनर्गठन द्वारा संचालित डिफ्यूजन गुणांक और बोरोन-टू-कार्बन (B/C) अनुपात के बीच एक प्रति-सहज (counter-intuitive) संबंध की पहचान करता है।

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

गतिक विकास (Dynamical Evolution)

CR-संचालित आउटफ्लो: CR इंजेक्शन ऊर्ध्वाधर गैस आउटफ्लो को प्रेरित करता है और पार्कर अस्थिरता संरचनाएं (चुंबकीय लूप) बनाता है। उच्च SN दर अधिक तीव्र विंड्स और तेज़ ऊर्ध्वाधर गैस वेगों की ओर ले जाती है।
चुंबकीय क्षेत्र का पुनर्गठन: CR बॉयन्सी ऊर्ध्वाधर चुंबकीय घटक बनाती है। इन ऊर्ध्वाधर क्षेत्रों की ताकत डिफ्यूजन दक्षता पर निर्भर करती है।

स्पेक्ट्रल गुण

स्लोप विचलन (Slope Deviations): विकसित CR स्पेक्ट्रा प्रारंभिक इंजेक्शन स्लोप ( $q \approx 4.1$ $q \approx 4.1$ ) से विचलन दिखाते हैं।
- रिलेटिविस्टिक ऊर्जाओं पर, प्राथमिक और द्वितीयक के बीच स्पेक्ट्रल स्लोप का अंतर ( $\Delta a_R$ ) डिफ्यूजन पावर इंडेक्स $\delta$ के बराबर होता है।
- नॉन-रिलेटिविस्टिक ऊर्जाओं पर, कूलिंग और स्पलेशन के कारण स्लोप अधिक तीव्र (steeper) होता है, जिसमें डिस्क हेलो की तुलना में हार्डर स्पेक्ट्रा दिखाई देता है।
ऊर्जा घनत्व शिखर (Energy Density Peaks): स्पेक्ट्रा इंजेक्शन, कूलिंग और डिफ्यूजन के बीच परस्पर क्रिया के कारण 10–30 GeV/nucleon के आसपास शिखर प्रदर्शित करते हैं।

अवलोकन योग्य: B/C और $^{10}\text{Be}/^{9}\text{Be}$ अनुपात

SN दर का प्रभाव:
- उच्च SN दर $\rightarrow$ निम्न B/C। तीव्र विंड्स एडवेक्शन को बढ़ाती हैं, जिससे प्राथमिक CR का निवास समय (residence time) कम हो जाता है जहाँ स्पलेशन होता है।
- उच्च SN दरें तेजी से पलायन और क्षय के कारण निम्न $^{10}\text{Be}/^{9}\text{Be}$ अनुपात की ओर भी ले जाती हैं।
डिफ्यूजन गुणांक ( $D_{\parallel}$ ) का प्रभाव:
- नवीन खोज: मानक लीकी-बॉक्स भविष्यवाणियों के विपरीत (जहाँ उच्च $D$ का अर्थ निम्न B/C है), यह मॉडल दिखाता है कि उच्च $D_{\parallel}$ से उच्च B/C होता है।
- तंत्र: उच्च डिफ्यूजन CR को क्षैतिज चुंबकीय रेखाओं के साथ अधिक समान रूप से वितरित करता है, जिससे स्थानीय CR दबाव प्रवणता (pressure gradients) कम हो जाती है जो चुंबकीय बॉयन्सी को संचालित करती है। इसके परिणामस्वरूप कमजोर ऊर्ध्वाधर चुंबकीय क्षेत्र ( $\langle |B_z|/|\vec{B}| \rangle$ घटता है) बनते हैं। कमजोर ऊर्ध्वाधर क्षेत्र प्राथमिक CR को डिस्क में अधिक समय तक रोक कर रखते हैं, जिससे द्वितीयक (बोरोन) का उत्पादन बढ़ जाता है।
रिजिडिटी इंडेक्स ( $\delta$ ) का प्रभाव:
- $\delta$ को 0.3 से बढ़ाकर 0.5 करने से रिलेटिविस्टिक ऊर्जाओं पर B/C अनुपात काफी तीव्र हो जाता है, जो सिमुलेशन परिणामों को AMS-02 के अवलोकन संबंधी डेटा ( $\delta \approx 1/3$ का समर्थन करता है) के साथ बेहतर तालमेल में लाता है।

5. महत्व और सीमाएँ (Significance and Limitations)

महत्व:

यह अध्ययन प्रमाणित करता है कि CR अवलोकनों (observables) की सटीक व्याख्या करने के लिए स्व-सुसंगत MHD सिमुलेशन आवश्यक हैं। CR द्वारा अनुभव किया जाने वाला "प्रभावी" डिफ्यूजन केवल इनपुट पैरामीटर $D_{\parallel}$ का फलन नहीं है, बल्कि चुंबकीय क्षेत्र की ज्यामिति द्वारा नियंत्रित होता है, जो स्वयं CR फीडबैक का उत्पाद है।
यह जटिल CR प्रसार मॉडलों की उच्च-परिशुद्धता डेटा (जैसे AMS-02 और Voyager) के साथ तुलना करने के लिए एक मजबूत संख्यात्मक ढांचा प्रदान करता है।

सीमाएँ और भविष्य का कार्य:

रिज़ॉल्यूशन: वर्तमान ग्रिड रिज़ॉल्यूशन (62.5 pc) अत्याधुनिक सिमुलेशन (1–4 pc) की तुलना में मोटा है, जिससे छोटे पैमाने की विक्षोभ (turbulence) छूट सकती है।
सरलीकृत भौतिकी: मॉडल में CR स्ट्रीमिंग, बहु-चरणीय (multi-phase) ISM तापीय संरचना, और विस्तृत SN तापीय/गतिज फीडबैक को छोड़ दिया गया है।
ज्यामिति: सिमुलेशन एक स्थानीय पैच का उपयोग करता है जिसमें आवधिक सीमाएं (periodic boundaries) हैं, जो पूर्ण गैलेक्टिक हेलो मॉडल की तुलना में पूर्वाग्रह (bias) पेश कर सकता है।

निष्कर्ष:
लेखकों ने सफलतापूर्वक प्रदर्शित किया है कि PIERNIK MHD कोड, जिसे CRESP एल्गोरिदम के साथ विस्तारित किया गया है, स्व-सुसंगत रूप से CR नाभिक के उत्पादन और प्रसार को मॉडल कर सकता है। परिणाम इस बात पर प्रकाश डालते हैं कि CR परिवहन संबंधी अवलोकन (जैसे B/C) गतिशील रूप से CR, गैस और चुंबकीय क्षेत्रों के बीच कपलिंग के प्रति संवेदनशील होते हैं, जो स्थिर, फेनोमेनोलॉजिकल मॉडलों की धारणाओं को चुनौती देते हैं।