Failed ejection and oscillations of a current-carrying filament balanced by gravity

समस्या
सौर फिलामेंट्स (प्रॉमिनेंस) कोरोना में निलंबित घने, ठंडे प्लाज्मा संरचनाएं हैं, जो अक्सर सौर विस्फोटों और कोरोनल मास इजेक्शन (CMEs) से जुड़े होते हैं। जबकि कई मॉडल फिलामेंट स्थिरता और विस्फोट का वर्णन करते हैं, मौजूदा संख्यात्मक सिमुलेशन में एक महत्वपूर्ण सीमा शून्य प्लाज्मा बीटा का अनुमान है, जो गुरुत्वाकर्षण के प्रभावों की उपेक्षा करता है। यथार्थवादी कोरोनल स्थितियाँ परिमित प्लाज्मा दबाव और गुरुत्वाकर्षण स्तरीकरण (gravitational stratification) शामिल करती हैं, जो फिलामेंट स्थिरता, द्रव्यमान लोडिंग और तापीय संरचना को समझने के लिए अत्यंत महत्वपूर्ण हैं। इसके अतिरिक्त, फिलामेंट विस्फोटों पर प्राप्त अवलोकन संबंधी डेटा अक्सर "विफल" घटनाओं को दर्शाते हैं जहाँ निष्कासन प्रक्रिया रुक जाती है और फिलामेंट वापस गिर जाता है, एक ऐसी घटना जिसके लिए अंतर्निहित भौतिकी को पूरी तरह से पकड़ने हेतु गुरुत्वाकर्षण बलों और यथार्थवादी वर्तमान विन्यास वाले मॉडलों की आवश्यकता होती है।

कार्यप्रणाली
लेखक सोलोवेव (2010) द्वारा प्रस्तावित एक गुरुत्वाकर्षण-संतुलित मॉडल का उपयोग करके एक फिलामेंट के अस्थिरता-पश्चात विकास की जांच करते हैं। अध्ययन एक दोहरी पद्धति का उपयोग करता है:

1. विश्लेषणात्मक जांच: लेखक सोलोवेव मॉडल के आधार पर विद्युत धारा घनत्व, प्लाज्मा घनत्व और उनके स्थानिक वितरण के लिए विश्लेषणात्मक व्यंजक (expressions) प्राप्त करते हैं। इस मॉडल में, एक विद्युत धारा ले जाने वाले घने फिलामेंट को उसके उप-फोटोस्फेरिक इमेज करंट और फिलामेंट करंट के बीच प्रतिकर्षक लोरेन्त्ज़ बल (Lorentz force) द्वारा गुरुत्वाकर्षण के विरुद्ध सहारा दिया जाता है।
2. संख्यात्मक सिमुलेशन: Lare2d संख्यात्मक कोड का उपयोग करके एक द्वि-आयामी मैग्नेटोहाइड्रोडायनामिक (MHD) मॉडल बनाया गया है। मॉडल एक कार्टेशियन समन्वय प्रणाली में रेसिस्टिव, कंप्रेसिबल MHD समीकरणों को हल करता है, जिसमें गुरुत्वाकर्षण और प्रतिरोधकता (resistivity) को स्पष्ट रूप से शामिल किया गया है। सिमुलेशन डोमेन ($x \in [-10, 10]$ Mm, $y \in [0, 40]$ Mm) को $2000 \times 4000$ ग्रिड द्वारा हल किया गया है। प्रारंभिक साम्यावस्था (equilibrium state) को एक चुंबकीय क्षेत्र $B_0 = 10^{-3}$ T, द्रव्यमान घनत्व $\rho_0 \sim 1.3 \times 10^{-9}$ kg m$^{-3}$, और तापमान $T \sim 13000$ K द्वारा अभिलक्षित किया गया है। फोटोस्फीयर और कोरोना के बीच संक्रमण क्षेत्र को सिम्युलेट करने के लिए एक संशोधित तापमान प्रोफाइल का उपयोग किया गया है।

प्रणाली को दो अलग-अलग परिदृश्यों के माध्यम से अस्थिर किया जाता है: (a) फिलामेंट में विद्युत धारा में वृद्धि (चुंबकीय क्षेत्र की शक्ति $B_0$ बढ़ाकर) और (b) फिलामेंट के द्रव्यमान घनत्व को कम करना। निष्कासन, रुकने और वापस गिरने की गतिशीलता को देखने के लिए सिस्टम के विकास को ट्रैक किया जाता है।

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

• विश्लेषणात्मक व्यंजकों का व्युत्पन्न: शोध पत्र कुल विद्युत धारा और फिलामेंट के भार सहित मॉडल मापदंडों के फलन (function) के रूप में विद्युत धारा घनत्व, प्लाज्मा घनत्व और उनके स्थानिक वितरण के लिए विश्लेषणात्मक सूत्र प्रदान करता है।
• गुरुत्वाकर्षण-संतुलित MHD मॉडल: अध्ययन एक 2D MHD मॉडल लागू करता है जो गुरुत्वाकर्षण प्रभावों की उपेक्षा करने वाले शून्य-बीटा मॉडलों की सीमाओं को संबोधित करते हुए, प्रतिकर्षक चुंबकीय बल को गुरुत्वाकर्षण के विरुद्ध स्पष्ट रूप से संतुलित करता है।
• विफल निष्कासन का लक्षण वर्णन: लेखक "विफल निष्कासन" (failed ejection) घटना का व्यवस्थित रूप से सिमुलेशन और लक्षण वर्णन करते हैं, जहाँ एक अस्थिर फिलामेंट को उत्सर्जित किया जाता है, एक विशिष्ट ऊंचाई पर रोका जाता है, और उसके बाद वापस गिरता है, जिससे एक दोलन चक्र (oscillatory cycle) में प्रवेश करता है।
• मात्रात्मक निर्भरता: अधिकतम निष्कासन ऊंचाइयों को अस्थिर करने वाली धाराओं के परिमाण और फिलामेंट प्लाज्मा के तनुकरण (dilution) की डिग्री, दोनों के फलन के रूप में निर्धारित किया गया है।

परिणाम

• दोलन व्यवहार: दोनों अस्थिरता व्यवस्थाओं (धारा वृद्धि और घनत्व कमी) में, फिलामेंट ऊपर की ओर उत्सर्जित होता है, एक निश्चित ऊंचाई पर रुकता है, और वापस गिरता है। यह चक्र लगभग 600 सेकंड की अवधि के साथ दोहराया जाता है।
• निष्कासन वेग: अस्थिरता तंत्र के आधार पर अधिकतम निष्कासन वेग भिन्न होते हैं: धारा-वृद्धि परिदृश्य के लिए लगभग 80 km s$^{-1}$ और घनत्व-कमी परिदृश्य के लिए लगभग 40 km s$^{-1}$।
• करंट शीट का निर्माण: उत्सर्जित फिलामेंट के नीचे, एक करंट शीट बनती है जहाँ चुंबकीय पुनर्संयोजन (magnetic reconnection) होता है। यह पुनर्संयोजन ऊपरी और निचले धाराओं को जोड़ने वाली चुंबकीय रेखाओं को बाधित करता है।
• मापदंड निर्भरता:
  • धारा (चुंबकीय क्षेत्र की शक्ति) में वृद्धि से उच्च निष्कासन ऊंचाई प्राप्त होती है।
  • फिलामेंट घनत्व को कम करने से भी उच्च निष्कासन ऊंचाई प्राप्त होती है, जिसमें कम घनत्व अनुपात ($\rho/\rho_0$) अधिक ऊंचाई की ओर ले जाते हैं।
  • घनत्व-कमी परिदृश्य में फिलामेंट अपनी अधिकतम ऊंचाई तक जल्दी पहुँचता है और वहां अधिक समय तक (200 सेकंड) रहता है, जबकि धारा-वृद्धि परिदृश्य (100 सेकंड) की तुलना में, हालांकि दोलन अवधि समान रहती है।
• असममिति (Asymmetry): धारा-वृद्धि के मामले में, फिलामेंट के नीचे का घनत्व संरचना अंततः विषमता विकसित करती है, जिससे फिलामेंट का ऊर्ध्वाधर अक्ष से विस्थापन होता है। घनत्व-कमी के मामले में, फिलामेंट ऊर्ध्वाधर अक्ष के साथ एक स्थिर अवस्था तक दोलन करता है।
• सत्यापन: एक विस्तारित ऊर्ध्वाधर डोमेन के साथ तुलनात्मक सिमुलेशन ने पुष्टि की कि मूल संख्यात्मक डोमेन का आकार चुंबकीय क्षेत्र विन्यास या प्राप्त अधिकतम ऊंचाई को प्रभावित नहीं करता है।

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