Improving Precision in Kinematic Weak Lensing with MIRoRS: Model-Independent Restoration of Reflection Symmetries

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

कमजोर गुरुत्वाकर्षण लेंसिंग (Weak gravitational lensing) ब्रह्मांड के द्रव्यमान वितरण का मानचित्रण करने के लिए एक प्राथमिक उपकरण है, लेकिन पारंपरिक विधियाँ "आकार शोर" (shape noise)—अर्थात बिना लेंस वाले आकाशगंगाओं की आंतरिक, यादृच्छिक दीर्घवृत्ताकारता (ellipticity) और अभिविन्यास (orientation) से ग्रस्त होती हैं। इससे निपटने के लिए, आमतौर पर हजारों से लाखों आकाशगंगाओं के सांख्यिकीय समूहों (statistical ensembles) की आवश्यकता होती है।

काइनेमैटिक लेंसिंग (Kinematic Lensing - KL) व्यक्तिगत आकाशगंगाओं के वेग क्षेत्रों (velocity fields) का उपयोग करके एक समाधान प्रदान करता है। एक बिना लेंस वाली डिस्क आकाशगंगा में विशिष्ट समरूपताएं होती हैं: इसके प्रमुख अक्ष (major axis) के सापेक्ष समरूपता (symmetry) और लघु अक्ष (minor axis) के सापेक्ष प्रति-समरूपता (antisymmetry)। कमजोर लेंसिंग शियर ($\gamma$) इन समरूपताओं को तोड़ देता है। जबकि पिछले KL तरीकों (जैसे, Donet & Wittman 2023) ने परिशुद्धता में सुधार किया है, वे आकाशगंगा के वेग प्रोफाइल के पैरामीट्रिक मॉडल (parametric models) पर निर्भर करते हैं। ये मॉडल अक्सर तब विफल हो जाते हैं जब वास्तविक आकाशगंगाएं आदर्श गोलाकारता या विशिष्ट वेग प्रोफाइलों से विचलित होती हैं, जिससे अवशिष्ट अनिश्चितताएं (residual uncertainties) (जैसे, $\sigma_{\gamma_\times} \approx 0.08$) उत्पन्न होती हैं जो स्वयं सिग्नल से भी अधिक होती हैं। इसके अलावा, शियर के क्रॉस-घटक ($\gamma_\times$) को फिट करना आकाशगंगा के आंतरिक अक्ष अनुपात (axis ratio) के बारे में धारणाओं के प्रति विशेष रूप से संवेदनशील है।

2. कार्यप्रणाली: MIRoRS

लेखक MIRoRS (Model-Independent Restoration of Reflection Symmetries) पेश करते हैं, जो एक ऐसी तकनीक है जो आकाशगंगा के वेग प्रोफाइल के लिए किसी विशिष्ट पैरामीट्रिक रूप को माने बिना $\gamma_\times$ को फिट करती है।

मुख्य अवधारणा (Core Concept):
यह विधि इस सिद्धांत पर कार्य करती है कि यदि सही शियर पैरामीटर लागू किया जाए, तो वह लेंसिंग प्रभाव को उलटने के लिए आकाशगंगा के वेग क्षेत्र की आंतरिक प्रतिबिंब समरूपताओं (reflection symmetries) को बहाल कर देगा।

एल्गोरिदम चरण (Algorithmic Steps):

1. डेटा रूपांतरण (Data Transformation): डिटेक्टर फ्रेम में देखे गए वेग और फोटोमेट्रिक डेटा से शुरू करते हुए, एल्गोरिदम डेटा को स्रोत प्ले (source plane) में वापस मैप करने के लिए रूपांतरणों की एक श्रृंखला लागू करता है:
   • स्थानांतरण (Translation): सेंट्रॉइड (centroid) को मूल बिंदु (origin) पर लाता है।
   • इनवर्स शियर (Inverse Shear): शियर विरूपण को हटाने के लिए $A^{-1}$ लागू करता है।
   • घूर्णन (Rotation): एक स्थिति कोण (position angle) पैरामीटर ($\phi$) का उपयोग करके काइनेमैटिक अक्षों को समन्वय अक्षों (coordinate axes) के साथ संरेखित करता है।
2. समरूपता बहाली (Symmetry Restoration): रूपांतरित क्षेत्र पर दो स्वतंत्र समरूपता संचालन किए जाते हैं:
   • प्रमुख अक्ष प्रतिबिंब (Major Axis Reflection): प्रमुख अक्ष के पार क्षेत्र को प्रतिबिंबित करता है।
   • लघु अक्ष व्युत्क्रमण (Minor Axis Inversion): लघु अक्ष के पार प्रतिबिंब बनाता है और चिह्न को उलट देता है (प्रति-समरूपता को ध्यान में रखते हुए)।
3. मॉडल निर्माण (Model Construction): दोनों रूपांतरित क्षेत्रों को एक एकल गैर-पैरामीट्रिक मॉडल बनाने के लिए औसत निकाला जाता है। इस मॉडल को फिर डिटेक्टर फ्रेम में वापस रूपांतरित किया जाता है।
4. लाइकलीहुड अनुकूलन (Likelihood Optimization): मॉडल की तुलना मूल डेटा से की जाती है ताकि एक अवशेष (residual) और लॉग-लाइकलीहुड की गणना की जा सके। अनुकूलन (dual annealing algorithm का उपयोग करके) $\gamma_\times$, स्थिति कोण $\phi$, सेंट्रॉइड, और केंद्रीय वेग जैसे मापदंडों को बदलकर लाइकलीहुड को अधिकतम करता है।
5. फोटोमेट्रिक बाधा (Photometric Constraint): अपभ्रष्टता (degeneracies) को तोड़ने और परिशुद्धता में सुधार करने के लिए, यह विधि फोटोमेट्रिक डेटा को शामिल करती है। यह देखे गए फोटोमेट्रिक स्थिति कोण और शियर रूपांतरण (DiGiorgio etл. 2021 से व्युत्पन्न) द्वारा अनुमानित काइनेमैटिक स्थिति कोण के बीच के अंतर की गणना करता है। इस पद को लाइकलीहुड फंक्शन में जोड़ा जाता है।
6. आउटलियर रिजेक्शन (Outlier Rejection): अवशेषों (residuals) पर एक ग्लोबल मोरन'स I टेस्ट (Global Moran's I test) लागू किया जाता है। अवशेषों में स्थानिक स्व-सहसंबंध (spatial autocorrelation) उन अनमॉडल भौतिक संरचनाओं (जैसे, वार्प्स, बार्स) को दर्शाता है जो समरूपता धारणाओं का उल्लंघन करते हैं। महत्वपूर्ण स्थानिक सहसंबंध वाले हेलो (halos) को फ़िल्टर कर दिया जाता है।

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

• मॉडल-स्वतंत्रता (Model-Independence): पिछले KL दृष्टिकोणों के विपरीत, MIRoRS आकाशगंगा के वेग प्रोफाइल (जैसे, टली-फिशर या आर्कटैंजेंट प्रोफाइल) के लिए किसी विशिष्ट पैरामीट्रिक धारणा का उपयोग नहीं करता है, जिससे यह आकाशगंगा की आकृति विज्ञान (morphology) में विचलन के प्रति मजबूत बनता है।
• हाइब्रिड डेटा एकीकरण (Hybrid Data Integration): यह एक ही लाइकलीहुड ढांचे के भीतर काइमैटिक (वेग क्षेत्र) और फोटोमेट्रिक (इसोफोथल आकार) डेटा को अनूठे रूप से जोड़ता है, जिससे काइमैटिक बनाम फोटोमेट्रिक अक्षों को प्रभावित करने वाले शियर के विभिन्न तरीकों का लाभ उठाया जाता है।
• सांख्यिकीय आउटलियर रिजेक्शन (Statistical Outlier Rejection): काइनेमैटिक लेंसिंग के लिए एक मानक फिल्टर के रूप में मोरन'स I टेस्ट का परिचय, जो शियर माप को पक्षपाती बनाने वाली जटिल आंतरिक संरचनाओं वाली आकाशगंगाओं की वस्तुनिष्ठ पहचान और निष्कासन की अनुमति देता है।
• **यथार्थवादी सिमुलेशन पर सत्यापन (Validation on Realistic Simulations):)_{} यह विधि आदर्श मॉक डेटा पर कठोरता से परीक्षण की गई और Illustris TNG सिमुलेशन के 358 हेलो के एक बड़े डेटासेट पर लागू की गई, जिससे पिछले पैरामीट्रिक अध्ययनों के साथ सीधा तुलना संभव हुई।

4. परिणाम (Results)

लेख तीन अलग-अलग विश्लेषणों के परिणाम प्रस्तुत करता है:

A. आदर्श मॉक डेटा (Idealized Mock Data):

• प्रदर्शन: संयुक्त काइमैटिक+फोटोमेट्रिक मॉडल ने विभिन्न शियर मानों, झुकाव (inclinations) और शोर स्तरों पर दोनों काइमैटिक-केवल और पारंपरिक पैरामीट्रिक मॉडलों को पछाड़ दिया।
• झुकाव संवेदनशीलता (Inclination Sensitivity): संयुक्त मॉडल ने मध्यम झुकाव ($40^\circ$–$60^\circ$) पर बेहतर प्रदर्शन दिखाया।
• स्थिति कोण संवेदनशीलता (Position Angle Sensitivity): काइमैटिक-केवल मॉडल पूरी तरह विफल रहा जब आकाशगंगा को समन्वय अक्षों के सापेक्ष $45^\circ$ घुमाया गया (क्योंकि $\gamma_\times$, $\gamma_+$ की नकल करता है उस फ्रेम में)। संयुक्त मॉडल मजबूत बना रहा क्योंकि फोटोमेट्रिक पद ने फिट को नियंत्रित किया।
• व्यवस्थित त्रुटियां (Systematics): लेखकों ने एक छोटा पूर्वाग्रह (bias) पहचाना जब $\gamma_+$ (समानांतर शियर) गैर-शून्य होता है, विशेष रूप से उच्च झुकाव पर, जहाँ यह फिट में द्वि-मोडल (bimodal) वितरण का कारण बन सकता है। हालाँकि, अनिश्चितता पैरामीट्रिक विधियों की तुलना में कम रही।

B. इलस्ट्रिस TNG सिमुलेशन विश्लेषण (Illustris TNG Simulation Analysis):

• डेटासेट: ज्ञात वास्तविक शियर $\gamma_\times = 0$ के साथ $z=0$ पर 358 हेलो।
• अनिश्चितता में कमी:
  • पैरामीट्रिक मॉडल (DW23): समग्र अनिश्चितता $\sigma \approx 0.078$।
  • MIRO RS (संयुक्त): समग्र अनिश्चितता घटकर $\sigma = 0.039$ हो गई।
  • कोर अनिश्चितता (Core Uncertainty): वितरण का संकीर्ण घटक (जो सुव्यवस्थित आकाशगंगाओं का प्रतिनिधित्व करता है) में सुधार होकर $\sigma_{core} = 0.016$ हो गया।
• पूर्वाग्रह (Bias): एक मामूली नकारात्मक पूर्वाग्रह ($\approx -0.007$) देखा गया, जिसका कारण सिमुलेटेड आकाशगंगाओं का विशिष्ट क्लॉकवाइज रोटेशन (जो वार्प्ड डिस्क की नकल करता है) है, जिसकी अपेक्षा है कि यह एक वास्तविक, यादृच्छिक उन्मुख समूह में औसत रूप से समाप्त हो जाएगा।

C. आउटलियर रिजेक्शन (Moran's I):

• मोरन'स I टेस्ट ($p \leq 0.001$) लागू करने से लगभग $44\%$ नमूने (1074 अवलोकनों में से 472) को फ़िल्टर कर दिया गया।
• अंतिम परिशुद्धता (Final Precision): फ़िल्टर करने के बाद, समग्र अनिश्चितता और गिरकर $\sigma = 0.028$ हो गई, जिसमें कोर अनिश्चितता $\sigma_{core} = 0.014$ रही।
• सत्यापन: अस्वीकृत हेलो के दृश्य निरीक्षण ने पुष्टि की कि उनमें उप-संरचनाएं (बार्स, वार्प्स) थीं जो समरूपता धारणाओं का उल्लंघन करती हैं, जबकि स्वीकृत हेलो ने स्वच्छ अवशेष (clean residuals) दिखाए।

5. महत्व (Significance)

यह कार्य सिंगल-सोर्स वीक लेंसिंग में एक महत्वपूर्ण छलांग है। पैरामीट्रिक धारणाओं से दूर हटकर और फोटोमेट्रिक बाधाओं को एकीकृत करके, MIRoRS उसी डेटासेट पर पिछले अत्याधुनिक तरीकों की तुलना में शियर ($\gamma_\times$) को मापने की अनिश्चितता को लगभग 2 से 3 गुना कम कर देता है।

एक ही दृष्टि रेखा (line of sight) पर $\sigma_{\gamma_\times} \approx 0.028$ (और आउटलियर रिमूवल के साथ संभावित रूप से इससे भी कम) प्राप्त करने की क्षमता यह सुझाव देती है कि काइनेमैटिक लेंसिंग अंततः विशाल सांख्यिकीय समूहों की आवश्यकता के बिना व्यक्तिगत हेलो के लिए सटीक द्रव्यमान माप प्रदान कर सकती है। यह विशिष्ट गैलेक्सी क्लस्टर्स के द्रव्यमान वितरण और छोटे पैमाने पर डार्क मैटर की प्रकृति का अध्ययन करने के लिए नए मार्ग खोलता है। यह कार्य भविष्य के उच्च-रिज़ॉल्यूशन स्पेक्ट्रोस्कोपिक सर्वेक्षणों (जैसे, ELT या SKA के साथ) में "आकार शोर" को संभालने के लिए एक मजबूत ढांचा भी प्रदान करता है।