Inferring lunar wake potentials from electron phase space densities

यहाँ JGR: Space Physics में प्रस्तुत पांडुलिपि "Inferring lunar wake potentials from electron phase space densities" का विस्तृत तकनीकी सारांश दिया गया है।

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

लूनर वेक (चंद्रमा की पूंछ) एक प्लाज्मा वॉइड (plasma void) है जो चंद्रमा के पीछे की ओर बनता है, जिसका कारण वैश्विक चुंबकीय क्षेत्र और वायुमंडल का अभाव है। इस वेक के भीतर विद्युत विभव संरचना ($\phi$) को समझना प्लाज्मा रिफिलिंग गतिकी (plasma refilling dynamics) के मॉडलिंग के लिए अत्यंत महत्वपूर्ण है। हालाँकि, इन विभवों का प्रत्यक्ष मापन सामान्यतः असंभव है क्योंकि:

• कमजोर क्षेत्र (Weak Fields): औसत समानांतर विद्युत क्षेत्र (parallel electric field) लगभग $0.06$ mV/m है, जो अंतरिक्ष यान के उपकरणों की संवेदनशीलता सीमा से काफी नीचे है।
• उपकरण सीमाएँ (Instrument Limitations): केवल इलेक्ट्रोस्टैटिक शॉक के पास स्थित स्थानीयकृत, तीव्र विभव वृद्धि को ही सीधे मापा जा सकता है।

परिणामस्वरूप, विभवों का अनुमान आमतौर पर इलेक्ट्रॉन फेज स्पेस डेंसिटी (PSD) मापन के आधार पर लगाया जाता है, जो अर्ध-स्थिर व्लासोव साम्यावस्था (quasi-static Vlasov equilibrium) ($f = f(H)$, जहाँ $H$ इलेक्ट्रॉन हैमिल्टोनियन है) की धारणा पर आधारित होता है। मौजूदा मानक, हैमिल्टोनियन शिफ्ट विधि (Hamiltonian shift method), लूनर वेक में दो महत्वपूर्ण सीमाओं का सामना करती है:

1. स्ट्राल-संचालित विषमता (Strahl-Driven Asymmetry): सौर पवन में सुप्राथर्मल इलेक्ट्रॉनों की एक क्षेत्र-संरेखित बीम (straal) होती है। यह वेक के सूर्य की ओर (sunward) और सूर्य से विपरीत (anti-sunward) पक्षों के बीच एक अंतर्निहित विषमता पैदा करती है, जिससे वैश्विक संबंध $f(H)$ बहु-मान वाली (multi-valued) हो जाती है और एकल-डोमेन व्युत्क्रमण (single-domain inversion) अमान्य हो जाता है।
2. शॉक-प्रेरित ट्रैपिंग (Shock-Induced Trapping): केंद्रीय वेक में, काउंटर-स्ट्रीमिंग आयन बीम आयनिक ध्वनिक शॉक (ion acoustic shocks) उत्पन्न करते हैं। ये शॉक नॉनलीन लैंडौ अनुनाद (nonlinear Landau resonance) के माध्यम से इलेक्ट्रॉनों को फँसाते (trap) हैं, जिससे "फ्लैट-टॉप" वेग वितरण बनते हैं जहाँ $f(H)$ अपरिभाषित हो जाता है। हैमिल्टोनियन शिफ्ट विधि यहाँ संघर्ष करती है क्योंकि यह उच्च-ऊर्जा वाले "पासिंग" इलेक्ट्रॉनों पर निर्भर करती है, जिनकी फेज स्पेस डेंसिटी अक्सर कम होती है और सांख्यिकी कमजोर होती है।

2. कार्यप्रणाली: हैमिल्टोनियन व्युत्क्रमण विधि (Methodology: The Hamiltonian Inversion Method)

लेखक एक हैमिल्टोनियन व्युत्क्रमण विधि प्रस्तावित करते हैं जो डोमेन-डिकंपोजिशन रणनीति के माध्यम से अर्ध-स्थिर स्थिति $f = f(H)$ का लाभ उठाकर पूर्ण स्थानिक विद्युत विभव प्रोफाइल $\phi(x)$ का अनुमान लगाती है। यह विधि मानती है कि चुंबकीय क्षेत्र की तीव्रता में चुंबकीय क्षेत्र रेखा के साथ नगण्य परिवर्तन होता है, जिससे समस्या 2D फेज स्पेस $(x, p_x)$ में बदल जाती है।

एल्गोरिदम पाँच चरणों में आगे बढ़ता है:

1. प्रारंभिक अनुमान और स्प्लिट पॉइंट (Initial Guess & Split Point):

   • एक प्रारंभिक विभव अनुमान $\tilde{\phi}(x)$ इलेक्ट्रॉन घनत्व से प्राप्त बोल्ट्ज़मैन संबंध से निकाला जाता है।
   • विभव न्यूनतम $x^*$ (जो घनत्व न्यूनतम के साथ मेल खाता है) को बाएं (sunward) और दाएं (anti-sunward) डोमेन के बीच प्रारंभिक विभाजन बिंदु के रूप में पहचाना जाता है।

2. डोमेन डिकंपोजिशन (Domain Decomposition):

   • बाएं/दाएं डोमेन (Left/Right Domains): उन क्षेत्रों के रूप में परिभाषित हैं जहाँ इलेक्ट्रॉन वितरण 'पासिंग' इलेक्ट्रॉनों द्वारा नियंत्रित होता है।
   • मध्य क्षेत्र (Middle Region): प्रभावी इलेक्ट्रॉन थर्मल वेलोसिटी $v_{th}(x)$ की गणना करके पहचाना जाता है। $x^*$ के पास $v_{th}(x)$ में एक तीव्र, स्थानीयकृत वृद्धि आयनिक ध्वनिक शॉक और फंसे हुए इलेक्ट्रॉनों की उपस्थिति का संकेत देती है। यह क्षेत्र $x_L$ और $x_R$ द्वारा सीमित है।

3. स्वतंत्र व्युत्क्रमण (बाएं और दाएं) (Independent Inversion - Left & Right):

   • प्रत्येक पक्ष के लिए, विधि यह मानती है कि $f = f(H)$ स्थानीय रूप से एकल-मान वाली है।
   • यह देखे गए फेज स्पेस पिक्सेल $(x, p_x)$ को उनके हैमिल्टोनियन मानों के आधार पर बिन (binning) करके एक कार्यात्मक संबंध $f_{interp}(\tilde{H})$ का निर्माण करती है।
   • एक त्रुटि फलन $L(\tilde{\phi})$ को L-BFGS-B एल्गोरिदम (एक क्वासी-न्यूटन विधि) का उपयोग करके न्यूनतम किया जाता है। इसका उद्देश्य वह विभव $\tilde{\phi}$ खोजना है जो देखे गए $\ln f$ और पुनर्निर्मित $\ln f_{interp}$ के बीच के अंतर को न्यूनतम करता है।
   • स्थानीय मिनिमा से बचने के लिए एक मल्टी-ग्रिड दृष्टिकोण (कोर्स से फाइन) का उपयोग किया जाता है।

4. प्रत्यक्ष व्युत्क्रमण (मध्य क्षेत्र) (Direct Inversion - Middle Region):

   • शॉक क्षेत्र में, $f(H)$ अपरिभाषित है। इसके बजाय, विधि फंसे हुए इलेक्ट्रॉनों के फ्लैट-टॉप स्ट्रक्चर का उपयोग करती है।
   • कटऑफ वेलोसिटी $v_{cut}(x)$ को उस वेग के रूप में पहचाना जाता है जहाँ PSD अपने प्लेटो (plateau) मान के 50% तक गिर जाता है।
   • विभव का अनुमान सीमा के सापेक्ष ऊर्जा संरक्षण के माध्यम से सीधे लगाया जाता है:
     $$\tilde{\phi}(x) = \tilde{\phi}(x_L) + \frac{m_e}{2e} [v_{cut}^2(x) - v_{cut}^2(x_L)]$$

5. स्टिचिंग (Stitching):

   • तीनों समाधानों को एक निरंतर प्रोफाइल में संयोजित किया जाता है। मध्य समाधान को बाएं डोमेन से जोड़ा जाता है, और दाएं डोमेन को $x_R$ पर निरंतरता सुनिश्चित करने के लिए एक स्थिर ऑफसेट द्वारा शिफ्ट किया जाता है।

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

• नवीन एल्गोरिदम: एक डोमेन-डिकंपोजिशन व्युत्क्रमण तकनीक का विकास जो स्ट्राल विषमता के कारण $f(H)$ की बहु-मान प्रकृति और शॉक ट्रैपिंग के कारण $f(H)$ की अपरिभाषित प्रकृति को स्पष्ट रूप से संभालती है।
• परावर्तित/फंसे हुए इलेक्ट्रॉनों का उपयोग: शिफ्ट विधि के विपरीत, जो परावर्तित इलेक्ट्रॉनों को त्याग देती है या शॉक में कम-सांख्यिकी वाले पासिंग इलेक्ट्रॉनों के साथ संघर्ष करती है, यह विधि विभव गिरावट का अनुमान लगाने के लिए फंसे हुए इलेक्ट्रॉनों के अच्छी तरह से मापे गए फ्लैट-टॉप वितरण का लाभ उठाती है।
• मजबूत सत्यापन: इस विधि को दो अलग-अलग विकास चरणों में वेक्टर पार्टिकल-इन-सेल (VPIC) सिमुलेशन के विरुद्ध कड़ाई से सत्यापित किया गया है और वास्तविक अंतरिक्ष यान डेटा पर लागू किया गया है।

4. परिणाम (Results)

A. सिमुलेशन सत्यापन (Simulation Validation)

• प्रारंभिक चरण (कोई शॉक नहीं, $r/R_L \approx 1.65$):
  • इसमें मजबूत स्ट्राल विषमता है लेकिन कोई केंद्रीय शॉक नहीं है।
  • परिणाम: हैमिल्टोनियन व्युत्क्रमण विधि ने सटीक रूप से विषम (asymmetric) विभव प्रोफाइल को पुनः प्राप्त किया। हैमिल्टोनियन शिफ्ट विधि वेक केंद्र के पास कम गिनती सांख्यिकी (counting statistics) के कारण विफल रही।
• अंतिम चरण (शॉक के साथ, $r/R_L \approx 5.5$):
  • इसमें आयनिक ध्वनिक शॉक और फ्लैट-टॉप इलेक्ट्रॉन वितरण मौजूद हैं।
  • परिणाम: डोमेन डिकंपोजिशन ने मध्य क्षेत्र की सफलतापूर्वक पहचान की। $v_{cut}$-आधारित प्रत्यक्ष व्युत्क्रमण ने शॉक से जुड़े विभव संवर्द्धन को सटीक रूप से पकड़ा, एक ऐसी विशेषता जिसे शिफ्ट विधि केवल उच्च अनिश्चितता के साथ ही अनुमानित कर सकती थी।

B. ARTEMIS अंतरिक्ष यान अनुप्रयोग (ARTEMIS Spacecraft Application)

इस विधि को दो ARTEMIS लूनर वेक क्रॉसिंग पर लागू किया गया:

1. प्रारंभिक चरण क्रॉसिंग (Early Stage Crossing):
   • अनुमानित सामान्यीकृत विभव गिरावट: $e\Delta\phi/T_e \sim 15$ ($\Delta\phi \approx 800$ V)।
   • स्ट्राल द्वारा संचालित मजबूत विषमता को कैप्चर किया।
   • परिणाम खुले क्षेत्रों में शिफ्ट विधि के साथ सहमत थे लेकिन केंद्र में बेहतर रिज़ॉल्यूशन प्रदान करते थे।
2. अंतिम चरण क्रॉसिंग (Late Stage Crossing):
   • अनुमानित सामान्यीकृत विभव गिरावट: $e\Delta\phi/T_e \sim 5$ ($\Delta\phi \approx 200$ V)।
   • केंद्रीय वेक में शॉक-संबंधित विभव संवर्द्धन को कैप्चर किया।
   • विषमता कम स्पष्ट थी, जो सिमुलेशन रुझानों के अनुरूप है।
   • शिफ्ट विधि के साथ विसंगतियां शॉक सीमाओं के पास देखी गईं, जिसका कारण फ्लैट-टॉप किनारे को सटीक रूप से परिभाषित करने में कठिनाई को माना गया।

5. महत्व (Significance)

• वैज्ञानिक प्रगति: यह विधि लूनर वेक के पूर्ण विद्युत विभव प्रोफाइल को पुनर्गठित करने के लिए पहला मजबूत, स्व-सुसंगत ढांचा प्रदान करती है, जिसमें जटिल केंद्रीय शॉक क्षेत्र भी शामिल हैं जिन्हें पहले चरित्रित करना कठिन था।
• व्यापक प्रयोज्यता: यह तकनीक केवल चंद्रमा तक सीमित नहीं है। यह किसी भी गैर-चुंबकीयकृत या कमजोर चुंबकीय पिंड (जैसे क्षुद्रग्रह, धूमकेतु, मंगल के चंद्रमा) के लिए लागू करने योग्य है जहाँ प्लाज्मा रिफिलिंग समान वेक संरचनाएं बनाती है और इलेक्ट्रॉन अर्ध-स्थिर साम्यावस्था में होते हैं।
• भविष्य की उपयोगिता: लेखक इस विधि को ARTEMIS डेटा के सांख्यिकीय नमूने पर लागू करने की योजना बना रहे हैं ताकि रेडियल दूरी और सौर पवन स्थितियों के रूप में लूनर वेक विभव के विकास को मैप किया जा सके, जो मौलिक प्लाज्मा प्रक्रियाओं जैसे कि इलेक्ट्रोस्टैटिक शॉक गठन और कण ट्रैपिंग के बारे में नई अंतर्दृष्टि प्रदान करेगा।