Correlating Particle Acceleration Rates with Plasma Conditions in Colliding Wind Binaries

तकनीकी सारांश: टकराने वाली पवन द्विआधारी प्रणालियों (Colliding Wind Binaries) में कण त्वरण दरों का प्लाज्मा स्थितियों के साथ सहसंबंध

समस्या विवरण
यद्यपि सुपरनोवा अवशेषों (SNe) को पारंपरिक रूप से $\sim$1 TeV तक के गैलेक्टिक कॉस्मिक किरणों के प्राथमिक त्वरक माना जाता रहा है, हाल ही में LHAASO प्रयोग द्वारा पहचाने गए तारा-निर्माण क्षेत्रों (SFRs) से प्राप्त अति-उच्च ऊर्जा (VHE) उत्सर्जन का संबंध सुपरनोवा से नहीं है। टकराने वाली पवन द्विआधारी प्रणालियों (CWBs) को इन कणों के संभावित स्रोतों के रूप में प्रस्तावित किया गया है, जो TeV और PeV ऊर्जाओं तक हैड्रोन को त्वरित करने में सक्षम हैं। हालांकि, CWBs के जटिल, विक्षुब्ध वातावरण में इस त्वरण को चलाने वाले विशिष्ट भौतिक तंत्र अभी भी अस्पष्ट हैं। पिछले अध्ययनों ने स्थापित किया है कि CWBs गैर-तापीय उत्सर्जन उत्पन्न करते हैं, लेकिन वे शास्त्रीय विसरित शॉक त्वरण (Diffusive Shock Acceleration - DSA) बनाम विक्षुब्ध, चुंबकीय रूप से संचालित प्रक्रियाओं (जैसे स्टोकेस्टिक त्वरण, चुंबकीय पुनर्संयोजन, और मिरर प्रभाव) के योगदान को पर्याप्त रूप से अलग करने में विफल रहे हैं। यह कार्य यह पहचानने का लक्ष्य रखता है कि कौन सी प्लाज्मा स्थितियाँ और भौतिक मात्राएँ कण त्वरण दरों के साथ सबसे दृढ़ता से सहसंबंध रखती हैं, ताकि CWBs में प्रमुख त्वरण तंत्रों को निर्धारित किया जा सके।

कार्यप्रणाली
लेखकों ने एक CWB प्रणाली के उच्च-रिज़ॉल्यूशन, तीन-आयामी मैग्नेटोहाइड्रोडायनामिक (MHD) सिमुलेशन किए, जिसे पेरिअस्ट्रॉन पासेज (periastron passage) के दौरान $\eta$ Carinae द्विआधारी के अनुरूप मॉडल किया गया था। सिमुलेशन सेटअप में शामिल था:

• MHD मॉडलिंग: सिस्टम को AMUN कोड का उपयोग करके हल किए गए आदर्श संपीड़ित MHD समीकरणों का उपयोग करके विकसित किया गया। द्विआधारी में एक प्राथमिक और एक माध्यमिक तारा शामिल था जिनके विशिष्ट पवन पैरामीटर (क्रमशः वेग 500 किमी/सेकंड और 2500 किमी/सेकंड) थे। सिमुलेशन ने प्रभावी पॉलीट्रोपिक इंडेक्स ($\Gamma_{\text{eff}}$) को एडियाबेटिक (5/3) से लेकर अत्यधिक कूलिंग (1.1) तक बदलकर विभिन्न प्लाज्मा संपीड़्यता शासन (compressibility regimes) की खोज की।
• रिज़ॉल्यूशन: शॉक, विक्षोभ (turbulence), और चुंबकीय पुनर्संयोजन संरचनाओं को सटीक रूप से पकड़ने के लिए एडेप्टिव मेश रिफाइनमेंट (AMR) का उपयोग किया गया, जिससे अधिकतम प्रभावी रिज़ॉल्यूशन $960 \times 1920 \times 1920$ सेल प्राप्त हुआ।
• टेस्ट पार्टिकल इंटीग्रेशन: MHD विकास के बाद, $10^5$ टेस्ट पार्टिकल्स (विभिन्न आयन प्रजातियों का प्रतिनिधित्व करते हुए, गणनात्मक दक्षता के लिए मुख्य रूप से Fe$^+$ का उपयोग किया गया, जिसके परिणाम H$^+$ पर सामान्यीकृत हैं) को विक्षुब्ध क्षेत्र में इंजेक्ट किया गया। उनके प्रक्षेप पथों (trajectories) को सापेक्षतावादी गति के समीकरण का उपयोग करके एकीकृत किया गया, जो उन्हें सिम्युलेटेड वेग और चुंबकीय क्षेत्रों से जोड़ता है।
• सहसंबंध विश्लेषण: अध्ययन ने स्थानीय कण त्वरण दर ($dE/dt$) को सिम्युलेशन स्नैपशॉट्स से निकाले गए चार विशिष्ट प्लाज्मा मात्राओं के साथ सहसंबंधित करने पर ध्यान केंद्रित किया: चुंबकीय क्षेत्र की शक्ति ($|B|$), धारा घनत्व ($|J|$), भंवरता (vorticity, $|\nabla \times \mathbf{u}|$), और वेग विचलन ($|\nabla \cdot \mathbf{u}|$)। सांख्यिकीय विश्लेषण 2D संभाव्यता वितरण कार्यों (2D-PDFs) और पियर्सन सहसंबंध गुणांकों का उपयोग करके किया गया, जिसमें केवल उन क्षेत्रों पर विचार किया गया जहाँ त्वरण दर 100 eV/s से अधिक थी।

मुख्य परिणाम

1. त्वरण तंत्र: विश्लेषण से पता चलता है कि कण त्वरण शास्त्रीय DSA के बजाय विक्षुब्ध, चुंबकीय रूप से संचालित प्रक्रियाओं द्वारा प्रधान है।

   • चुंबकीय क्षेत्र और भंवरता: त्वरण दर और चुंबकीय क्षेत्र की शक्ति ($|B|$) तथा भंवरता ($|\nabla \times \mathbf{u}|$) दोनों के बीच मजबूत सकारात्मक सहसंबंध पाए गए। पियर्सन सहसंबंध गुणांक इन चरों के लिए 0.46 से 0.55 के बीच थे। प्रतिगमन ढलान (regression slopes) सुपरलीनियर (1.79 तक) थे, जो इंगित करते हैं कि मजबूत चुंबकीय क्षेत्र और उच्च विक्षोभ ऊर्जा लाभ को महत्वपूर्ण रूप से बढ़ाते हैं।
   • वेग विचलन: इसके विपरीत, त्वरण दर और वेग विचलन (जो शॉक संपीड़न का एक प्रॉक्सी है) के बीच सहसंबंध कमजोर या नगण्य (पियर्सन गुणांक 0.18–0.33, ऋणात्मक या शून्य के करीब ढलान के साथ) था। यह सुझाव देता है कि सुस्पष्ट, समतलीय शॉक (planar shocks) दृश्य त्वरण में सीमित भूमिका निभाते हैं।
   • धारा घनत्व: धारा घनत्व ने मध्यम सहसंबंध (0.44–0.51) दिखाया लेकिन कम प्रतिगमन ढलान प्रदर्शित किए। लेखक इसे चुंबकीय क्षेत्र की जटिलता के एक अप्रत्यक्ष संकेतक के रूप में व्याख्या करते हैं, न कि एक प्रत्यक्ष चालक के रूप में, क्योंकि आदर्श MHD सीमा में प्रतिरोधक प्रभाव (resistive effects) नगण्य हैं।

2. ऊर्जा वितरण:

   • कण कुछ घंटों से लेकर दिनों के भीतर थर्मल वितरण से सापेक्षवादी पावर-लॉ (power-laws) में विकसित हुए, सैकड़ों TeV की ऊर्जा तक पहुँच गए।
   • अधिकतम ऊर्जा और ऊर्जा स्पेक्ट्रम में "नी" (knee) की स्थिति सतह चुंबकीय क्षेत्र की शक्ति ($B_*$) पर निर्भर करती है, जो हिलस कन्फाइनमेंट मानदंड (Hillas confinement criterion) के अनुरूप है।
   • पॉलीट्रोपिक इंडेक्स ($\Gamma_{\text{eff}}$) ने स्पेक्ट्रल आकार को महत्वपूर्ण रूप से प्रभावित किया। अधिक कूलिंग (कम $\Gamma_{\text{eff}}$) वाले मॉडलों ने कठोर स्पेक्ट्रा (harder spectra) और उच्च त्वरण दक्षता उत्पन्न की क्योंकि उनमें विक्षोभ और चुंबकीय क्षेत्र की जटिलता अधिक थी। एडियाबेटिक मॉडल ($\Gamma_{\text{eff}} = 5/3$) ने तीव्र ढलान (steeper slopes) और कम कुशल त्वरण दिखाया।

3. स्थानिक वितरण: उच्चतम त्वरण वाले क्षेत्र सख्ती से सबसे मजबूत शॉक फ्रंट्स (वेग विचलन) के साथ मेल नहीं खाते थे, बल्कि विक्षुब्ध डाउनस्ट्रीम क्षेत्रों के साथ मेल खाते थे जहाँ चुंबकीय क्षेत्र प्रवर्धित और जटिल होते हैं।

महत्व और दावे
यह शोध पत्र दावा करता है कि यह CWBs में कण त्वरण दक्षता और स्थानीय प्लाज्मा स्थितियों के बीच एक निर्णायक लिंक प्रदान करता है, जो स्थिर-अवस्था (steady-state) धारणाओं से आगे जाता है। इसकी प्राथमिक महत्ता यह प्रदर्शित करने में है कि विक्षोभ और चुंबकीय क्षेत्र की जटिलता कण त्वरण के प्रमुख चालक हैं, जो शास्त्रीय विसरित शॉक त्वरण (DSA) को पीछे छोड़ देते हैं।

लेखक निष्कर्ष निकालते हैं कि:

• CWBs सैकड़ों TeV या यहाँ तक कि PeV ऊर्जाओं तक हैड्रोनिक कणों को त्वरित कर सकते हैं, जिसमें 1% से अधिक कण VHE रेंज तक पहुँचते हैं, बशर्ते कि मध्यम चुंबकन (moderate magnetization, >1 G) हो।
• त्वरण प्रक्रिया स्वाभाविक रूप से स्टोकेस्टिक और चुंबकीय रूप से संचालित है, जो सरल शॉक क्रॉसिंग के बजाय विक्षुब्ध संरचनाओं (Alfvén waves, magnetic mirrors, reconnection) के साथ कणों की अंतःक्रिया पर निर्भर करती है।
• इन निष्कर्षों के भविष्य के उच्च-संवेदनशीलता वाले $\gamma$-रे अवलोकनों (जैसे, LACT और CTAO द्वारा) की व्याख्या करने के लिए प्रत्यक्ष निहितार्थ हैं, जो यह सुझाव देते हैं कि CWBs की कॉस्मिक किरण स्रोतों के रूप में पहचान क्षमता अत्यधिक विक्षोभ और चुंबकीय क्षेत्र की जटिलता पर निर्भर करती है, जो कूलिंग दक्षता और कक्षीय चरण (orbital phase) द्वारा नियंत्रित होती है।

यह कार्य नए प्रयोगात्मक सेटअप का प्रस्ताव नहीं करता है, बल्कि मौजूदा और भविष्य के अवलोकन संबंधी डेटा की व्याख्या करने के लिए एक सैद्धांतिक ढांचा प्रदान करता है, इस बात पर जोर देता है कि कुछ CWBs में मजबूत गामा-रे डिटेक्शन की अनुपस्थिति त्वरण क्षमता की कमी के बजाय कम विक्षोभ या कम चुंबकन के कारण हो सकती है।