Forward-modelling Milky Way Cepheids: selection effects and physical priors in the Gaia-HST calibration

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

यह शोध पत्र हबल तनाव (Hubble tension) को संबोधित करता है, जो स्थानीय दूरी-सीढ़ी (local distance-ladder) द्वारा निर्धारित हबल स्थिरांक ($H_0$) और कॉस्मिक माइक्रोवेव बैकग्राउंड (CMB) से प्राप्त प्रारंभिक ब्रह्मांड के निष्कर्षों के बीच $\sim5\sigma$ का अंतर है। $H_0$ पर सबसे सटीक स्थानीय प्रतिबंध, मिल्की वे (MW) के सेफिड वेरिएबल्स (Cepheid variables) के अंशांकन (calibration) पर निर्भर करता है, जो Gaia EDR3 पैरलैक्स (parallaxes) और हबल स्पेस टेलीस्कोप (HST) फोटोमेट्री का उपयोग करता है।

Högas & Mörtsell (2026, HM26) के हालिया कार्य ने दावा किया कि मिल्की वे सेफिड्स के लिए एक यूनिफॉर्म-इन-वॉल्यूम डिस्टेंस प्रायर ($\pi(d) \propto d^2$) अपनाने से हबल तनाव कम हो जाता है। हालांकि, इस दृष्टिकोण ने सैंपल सिलेक्शन फंक्शन (sample selection function) (उन विशिष्ट मानदंडों की अनदेखी की जिनका उपयोग अवलोकन के लिए सितारों को चुनने में किया गया था) की उपेक्षा की। लेखक तर्क देते हैं कि चयन प्रभावों (selection effects) को निर्दिष्ट करते समय यदि सिलेक्शन फंक्शन को नजरअंदाज किया जाता है, तो यह पक्षपाती निष्कर्ष (biased inferences) उत्पन्न करता है, जिससे ऐसे मॉक पॉपुलेशन बनते हैं जो प्रेक्षित डेटा के साथ सुसंगत नहीं होते। यह शोध पत्र एक पूर्णतः फॉरवर्ड-मॉडलिंग बेयसियन फ्रेमवर्क (forward-modeled Bayesian framework) का निर्माण करके, जो गैलेक्टिक ज्यामिति, चयन प्रभावों और भौतिक प्रायर (physical priors) को स्वतः समाहित करता है, इन धारणाओं का कठोरता से परीक्षण करने का लक्ष्य रखता है।

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

लेखकों ने एक बेयसियन फॉरवर्ड मॉडल विकसित किया है जो पीरियड-ल्यूमिनोसिटी (PL) संबंध के मापदंडों, Gaia पैरलैक्स ज़ीरो-पॉइंट ऑफसेट ($\delta\varpi$), और व्यक्तिगत सितारों की दूरियों का संयुक्त रूप से अनुमान लगाता है।

• डेटा (Data): विश्लेषण में दो अलग-अलग HST अभियानों (Cycle 22 और Cycle 27) के 66 मिल्की वे सेफिड्स का उपयोग किया गया है, जिनके साथ Gaia EDR3 पैरलैक्स भी उपलब्ध हैं, और इसमें लार्ज मैगेलैनिक क्लाउड (LMC) तथा गैलेक्सी N4258 से प्राप्त ज्यामितीय एंकरों (geometric anchors) को भी शामिल किया गया है।
  • C22: $\sim6.6$ kpc तक के लंबे आवर्त वाले सेफिड्स ($P > 8$ दिन) को लक्षित करता है।
  • C27: निकटवर्ती, उच्च-पैरलैक्स वाले सेफिड्स ($\varpi_{phot} > 0.8$ mas) को लक्षित करता है ताकि विसंगतियों (degeneracies) को दूर किया जा सके।
• भौतिक प्रायर (दूरी के लिए) (Physical Priors - Distance): यूनिफॉर्म-इन-वॉल्यूम प्रायर के बजाय, मॉडल मिल्की वे सेफिड्स के लिए एक थिन-डिस्क ज्योमेट्री प्रायर (thin-disk geometry prior) अपनाता है:
  $$\pi(d, \ell, b) \propto d^2 \exp\left(-\frac{R_{GC}}{R_d}\right) \exp\left(-\frac{|z|}{z_d}\right)$$
  जहाँ $R_{GC}$ गैलेक्टोसेंट्रिक रेडियस है, $z$ प्ले के ऊपर की ऊंचाई है, और $R_d, z_d$ स्केल लेंथ और स्केल हाइट हैं। यह आकाशगंगा में सेफिड्स के वास्तविक वितरण को दर्शाता है।
• सिलेक्शन मॉडलिंग (Selection Modelling): मुख्य नवाचार सिलेक्शन फंक्शन का स्पष्ट विश्लेत्मक उपचार है। लेखक सभी लेटेंट वेरिएबल्स (दूरी, आवर्त, धात्विकता) और प्रेक्षित मात्राओं पर समाकलन (integration) करके डिटेक्शन प्रोबेबिलिटी $p(S=1|\Lambda)$ को व्युत्पन्न करते हैं, जिसमें निम्नलिखित कट्स (cuts) को ध्यान में रखा जाता है:
  • आवर्त ($P$)
  • आभासी परिमाण ($m_W$)
  • पैरलैक्स ($\varpi$)
  • एक्सटिंक्शन ($A_H$)
    नीचे दिए गए मॉडल में चयन फंक्शन को एक हार्ड कट के बजाय एक स्मूथ ट्रांजिशन (संचयी सामान्य वितरण $\Phi$ का उपयोग करके) के रूप में मॉडल किया गया है, जिससे उच्च-आयामी इंटीग्रल को दूरी और आकाश की स्थिति पर एक सुलभ 2D इंटीग्रल में बदला जा सका।
• इन्फरेंस इंजन (Inference Engine): मॉडल को JAX में लागू किया गया है और numpyro के माध्यम से No-U-Turn Sampler (NUTS) का उपयोग करके सैंपल किया गया है। इसमें 91 सैंपल्ड पैरामीटर्स शामिल हैं, जिनमें वैश्विक PL पैरामीटर्स, पैरलैक्स ऑफसेट, इनट्रिन्सिक स्कैटर और प्रति-तारा दूरियां शामिल हैं।

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

1. सिलेक्शन के साथ प्रथम फॉरवर्ड मॉडल: यह मिल्की वे सेफिड आबादी का पहला बेयसियन फॉरवर्ड मॉडल है जो चयन फंक्शन और गैलेक्टिक डिस्क ज्योमेट्री को एक साथ स्पष्ट रूप से समाहित करता है।
2. विश्लेत्मक चयन उपचार (Analytic Selection Treatment): लेखकों ने जटिल, बहु-आयामी चयन कट्स (परिमाण, पैरलैक्स, आवर्त, एक्सटिंक्शन) के लिए डिटेक्शन प्रोबेबिलिटी की गणना करने के लिए एक विश्लेत्मक विधि विकसित की है, जिससे 9-आयामी इंटीग्रल को मल्टीवेरिएट नॉर्मल CDF वाले संख्यात्मक रूप से सुलभ रूप में बदला गया।
3. "कम हुए तनाव" के दावे का खंडन (Debunking the "Reduced Tension" Claim): यह शोध पत्र सिद्ध करता है कि HM26 द्वारा रिपोर्ट किया गया कम हुआ हबल तनाव, चयन प्रभावों की अनदेखी का एक आर्टिफैक्ट (artefact) है। लेखक दिखाते हैं कि बिना सिलेक्शन मॉडलिंग के यूनिफॉर्म-इन-वॉल्यूम प्रायर का उपयोग करने से परिणाम पक्षपाती हो जाते हैं, जबकि ज्योमेट्री और सिलेक्शन दोनों का सुसंगत उपचार करने पर परिणाम मानक SH0ES विश्लेषण के अनुरूप रहते हैं।
4. वैधीकरण (Validation): मॉडल पोस्टीरियर प्रेडिक्टिव चेक्स (PPC) को सफलतापूर्वक पास करता है, जो पैरलैक्स, मैग्नीट्यूड और पीरियड के प्रेक्षित वितरणों को सफलतापूर्वक पुनरुत्पादित करता है, जबकि सिलेक्शन को अनदेखा करने वाले मॉडल इन चेक्स में विफल रहते हैं।

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

• पीरियड-ल्यूमिनोसिटी ज़ीरो-पॉइंट: फिडुशियल मॉडल (डिस्क प्रायर + सिलेक्शन) का ज़ीरो-पॉइंट $M_W^H,1 = -5.909 \pm 0.022$ mag है। यह $0.7\sigma$ के भीतर आधारभूत SH0ES मान ($M_W^H,1 \approx -5.894$) के साथ सहमत है।
• Gaia पैरलैक्स ऑफसेट: अनुमानित ज़ीरो-पॉइंट ऑफसेट $\delta\varpi = -12.4 \pm 5.3$ $\mu$as है, जो स्वतंत्र सत्यापन और पिछले SH0ES विश्लेषणों के अनुरूप है।
• सिलेक्शन मॉडलिंग का प्रभाव:
  • सिलेक्शन की अनदेखी करना: जब सिलेक्शन को अनदेखा किया जाता है और यूनिफॉर्म-इन-वॉल्यूम प्रायर का उपयोग किया जाता है (HM26 की नकल करते हुए), तो ज़ीरो-पॉइंट अधिक चमकदार मानों ($\sim -5.95$ से $-5.98$ mag) की ओर शिफ्ट हो जाता है और $\delta\varpi$ अधिक नकारात्मक ($\sim -26$ से $-36$ $\mu$as) हो जाता है। यह शिफ्ट कम $H_0$ के अनुरूप है, जो कृत्रिम रूप से हबल तनाव को कम करता है।
  • प्रायर की भूमिका: एक बार जब सिलेक्शन को सही ढंग से मॉडल कर लिया जाता है, तो विशिष्ट दूरी प्रायर (डिस्क बनाम यूनिफॉर्म-इन-वॉल्यूम) का प्रभाव नगण्य होता है ($<0.002$ mag शिफ्ट)। पक्षपात मुख्य रूप से सिलेक्शन मॉडल की अनुपस्थिति से उत्पन्न होता है, न कि स्वयं प्रायर से।
• इनट्रिन्सिक स्कैटर (Intrinsic Scatter): मॉडल दोनों अभियानों के लिए अलग-अलग इनट्रिन्सिक स्कैटर का अनुमान लगाता है: $\sigma_{C22} \approx 0.071$ mag और $\sigma_{C27} \approx 0.040$ mag, जो संभवतः अधिक दूरस्थ C22 सैंपल में अवशिष्ट एक्सटिंक्शन अनिश्चितताओं के कारण है।
• हबल स्थिरांक निहितार्थ: SH0ES Gaia बाधा को लेखकों के अनुमानित ज़ीरो-पॉइंट से बदलने पर $H_0 = 73.12 \pm 1.08$ km s$^{-1}$ Mpc$^{-1}$ प्राप्त होता है। यह बदलाव नगण्य ($\sim0.1\sigma$) है, जो पुष्टि करता है कि हबल तनाव $\sim5\sigma$ पर बना हुआ है।

5. महत्व (Significance)

• सांख्यिकीय कठोरता (Statistical Rigor): यह शोध पत्र स्थापित करता है कि कॉस्मिक डिस्टेंस लैडर के प्रतिशत-स्तर की सटीकता के लिए सिलेक्शन प्रभावों का एक सिद्धांतबद्ध उपचार आवश्यक है। सिलेक्शन को अनदेखा करने से पक्षपाती जेनेरेटिव मॉडल बनते हैं जो प्रेक्षित डेटा के साथ असंगत होते हैं।
• विसंगतियों का समाधान: यह स्पष्ट करता है कि HM26 द्वारा दावा किया गया हबल तनाव का "न्यूनीकरण" नए भौतिक अंतर्दृष्टि के कारण नहीं, बल्कि एक सांख्यिकीय त्रुटि (मिसिंग सिलेक्शन मॉडलिंग) के कारण था।
• भविष्य का ढांचा: यह कार्य दूरी-सीढ़ी विश्लेषणों के लिए एक नया मानक स्थापित करता है। यह अन्य दूरी संकेतकों (जैसे टाइप Ia सुपरनोवा, TRGB) के लिए भविष्य के बेयसियन फॉरवर्ड मॉडलिंग के लिए एक टेम्पलेट प्रदान करता है और सुझाव देता है कि आगामी उच्च-परिशुद्धता ब्रह्मांडीय सर्वेक्षणों के लिए सिमुलेशन-आधारित अनुमान (simulation-based inference) या कठोर फॉरवर्ड मॉडलिंग आवश्यक है।
• SH0ES का वैधीकरण: इसके परिणाम SH0ES डिस्टेंस लैडर की मजबूती और हबल तनाव की निरंतरता को पुष्ट करते हैं, जिससे पता चलता है कि तनाव संभवतः भौतिक है न कि स्थानीय सेफिड कैलिब्रेशन के सांख्यिकीय प्रबंधन का परिणाम।