Constrained Gaussian-process bridge prior for neutron-star equation-of-state inference

तकनीकी सारांश: न्यूट्रॉन-तारा समीकरण-अवस्था (Equation-of-State) अनुमान के लिए बाधित गॉसियन-प्रोसेस ब्रिज प्रायर (Constrained Gaussian-Process Bridge Prior)

समस्या विवरण
मैक्रोस्कोपिक अवलोकनीय (observables) से न्यूट्रॉन-तारा (NS) समीकरण-अवस्था (EoS) का अनुमान लगाने के लिए ऐसे प्रायर (priors) का निर्माण करना आवश्यक है जो मॉडल-स्वतंत्र हों लेकिन भौतिक रूप से सुसंगत भी हों। मौजूदा गैर-पैरामीट्रिक दृष्टिकोण, जैसे कि मानक गॉसियन प्रोसेस (GPs), वैश्विक भौतिक बाधाओं को एक साथ लागू करने में संघर्ष करते हैं: ऊष्मप्रवैगिकी सुसंगतता (thermodynamic consistency), यांत्रिक स्थिरता (mechanical stability), और कार्य-कारणता (causality, $c_s \leq c$)। जबकि पैरामीट्रिक मॉडल इन बाधाओं को संतुष्ट कर सकते हैं, वे अक्सर अनियंत्रित सहसंबंधों वाले पक्षपाती अनुमानक (biased estimators) के रूप में कार्य करते हैं। इसके विपरीत, विधियाँ जो उच्च-घनत्व की बाधाओं (perturbative QCD से) और निम्न-घनत्व की बाधाओं (chiral EFT से) को शामिल करती हैं, वे अक्सर "शूटिंग" प्रक्रियाओं, मध्यवर्ती लाइकलीहुड फंक्शनों, या विभिन्न श्रेणियों के बीच 'एड हॉक' स्विचिंग पर निर्भर करती हैं। ये दृष्टिकोण मिलान बिंदुओं (matching points) पर निर्भरता पैदा कर सकते हैं या अधिकांश उत्पन्न EoS को अस्वीकार कर सकते हैं क्योंकि बाधाएं प्रायर जनरेशन में सीधे तौरते शामिल नहीं होती हैं।

कार्यप्रणाली
लेखक एक नया गैर-पैरामीट्रिक प्रायर जनरेशन विधि प्रस्तावित करते हैं जो एक बाधित गॉसियन-प्रोसेस (GP) ब्रिज पर आधारित है। यह विधि उन EoS का निर्माण करती है जो निम्न-घनत्व सीमाओं (chiral EFT से) और उच्च-घनत्व सीमाओं (pQCD से) को जोड़ती हैं, जबकि भौतिक बाधाओं का निर्माण द्वारा ही कड़ाई से पालन करती हैं। यह प्रक्रिया दो मुख्य चरणों से मिलकर बनी है:

1. स्व-समान (फ्रैक्टल) परिशोधन (Self-Similar/Fractal Refinement):
   व्हाइट नॉइज़ के बजाय, यह विधि एक पुनरावर्ती, स्व-समान परिशोधन प्रक्रिया के माध्यम से "अधिकतम शोर" वाले EoS का स्थान निर्मित करती है। थर्मोडायनामिक ट्रिपलेट्स $\beta_L = (\mu_L, n_L, p_L)$ और $\beta_H = (\mu_H, n_H, p_H)$ द्वारा परिभाषित सीमा स्थितियों से शुरू करते हुए, एल्गोरिदम भौतिक रूप से अनुमत आयतन $V_{\beta_L, \beta_H}$ के भीतर मध्यवर्ती बिंदु $\beta_0$ को पुनरावर्ती रूप से पेश करता है। यह आयतन यांत्रिक स्थिरता (जहाँ $n(\mu)$ एकल-मान वाला है), कार्य-कारणता ($\mu/n \cdot dn/d\mu \geq 1$), और ऊष्मप्रवैगिकी सुसंगतता (सही दबाव एकीकरण) की बाधाओं द्वारा परिभाषित है। बिंदुओं को इस अनुमत आयतन के भीतर समान रूप से नमूनाबद्ध (sample) किया जाता है, जिससे सभी घनत्व पैमानों पर संरचना वाले लेकिन न्यूनतम अल्प-दूरी सहसंबंधों वाले EoS बनते हैं।

2. विसरण स्मूथिंग (बाधित GP ब्रिज):
   स्थानीय सहसंबंधों को पेश करने और अवकलनीयता (differentiability) सुनिश्चित करने के लिए, लेखक एक फ्लो समय $\tau$ के फलन के रूप में रासायनिक क्षमता $\mu(n)$ पर एक हीट-कर्नेल स्मूथिंग (डिफ्यूजन) लागू करते हैं। इसका विकास डिफ्यूजन समीकरण द्वारा नियंत्रित होता है:
   $$\partial_\tau \mu(n, \tau) = \partial_n [D(n) \partial_n \mu(n, \tau)]$$
   महत्वपूर्ण रूप से, डिफ्यूजन गुणांक $D(n)$ को भौतिक बाधाओं को संरक्षित करने के लिए चुना गया है:

• यांत्रिक स्थिरता: यदि $D(n) > 0$ हो, तो यह संरक्षित रहती है।
• कार्य-कारणता (Causality): यदि $D'' - D'/n \leq 0$ हो, तो यह संरक्षित रहती है।
• ऊष्मप्रवैगिकी सुसंगतता: यह प्रक्रिया डोमेन के भीतर ऊर्जा घनत्व को संरक्षित करती है, जिसमें सीमाओं पर केवल नगण्य उल्लंघन होते हैं जिन्हें $D(n)$ को कम करके या ऊर्जा को पुनः इंजेक्ट करके कम किया जाता है।

$D(n) \propto n^2$ चुनकर, यह विधि एक सहसंबंध लंबाई $\sigma$ प्राप्त करती है जो घनत्व का एक निश्चित अंश ($\sigma/n = c$) है, जिससे निम्न घनत्वों पर सूक्ष्म संरचना (fine structure) और pQCD घनत्वों तक विस्तृत विस्तार संभव होता है। इसके परिणामस्वरूप एक संशोधित GP-ब्रिज निर्माण होता है जहाँ स्थानीय सहसंबंध संरचना गॉसियन है, लेकिन वैश्विक संरचना भौतिक बाधाओं द्वारा आवश्यक गैर-तुच्छ (non-trivial) एंटी-कोरिलेशन बनाए रखती है।

मुख्य योगदान

• एकीकृत ढांचा: यह विधि बिना किसी मध्यवर्ती लाइकलीहुड या शूटिंग प्रक्रिया के, निम्न-घनत्व (chiral EFT) और उच्च-घनत्व (pQCD) बाधाओं को एक एकल, एकीकृत प्रायर में एकीकृत करती है।
• बाधा संतुष्टि: यह निर्माण गारंटी देता है कि प्रत्येक उत्पन्न EoS डिजाइन द्वारा ऊष्मप्रवैगिकी रूप से सुसंगत, यांत्रिक रूप से स्थिर और कार्य-कारणता (causal) का पालन करने वाला है, जिससे अमान्य नमूनों को अस्वीकार करने की आवश्यकता समाप्त हो जाती है।
• लचीला सहसंबंध नियंत्रण: डिफ्यूजन दृष्टिकोण ध्वनि-वेग (sound-speed) के सहसंबंध लंबाई को ट्यून करने की अनुमति देता है, जो रूढ़िवादी (कम सहसंबंध) और सिद्धांत-सूचित (लंबा सहसंबंध) प्रायरों के बीच के अंतर को पाटता है।
• वैश्विक एंटी-कोरिलेशन: यह विधि ध्वनि-वेग में दीर्घ-रेंज एंटी-कोरिलेशन स्वाभाविक रूप से प्रेरित करती है: जो EoS मध्यम घनत्वों पर सख्त (stiffen) होते हैं, उन्हें कार्य-कारणता और pQCD सीमाओं को संतुष्ट करने के लिए उच्च घनत्वों पर नरम (soften) होना आवश्यक है।

परिणाम
लेखकों ने द्रव्यमान मापन (PSR J0348+0432 और PSR J1614−2230), NICER (PSR J0740+6620, PSR J0614−3329, PSR J0437−4715) के संयुक्त द्रव्यमान-त्रिज्या प्रतिबंधों, और GW170817 से टाइडल डिफॉर्मेबिलिटी सहित अपडेटेड ऑब्जर्वेशनल डेटा का उपयोग करके EoS का अनुमान लगाने के लिए इस ढांचे को लागू किया।

• EoS व्यवहार: पोस्टीरियर वितरण लगातार $2M_\odot$ सितारों को सहारा देने के लिए chiral-EFT शासन के तुरंत बाद EoS के सख्त होने (stiffening) को दिखाते हैं, जिसके बाद उच्च घनत्वों ($n \gtrsim 1 \text{ fm}^{-3}$) पर एक व्यवस्थित मृदुकरण (softening) होता है।
• सहसंबंध लंबाई निर्भरता: हालांकि EoS के स्थानीय गुण (जैसे, पीक साउंड स्पीड) चुनी गई सहसंबंध लंबाई पर निर्भर करते हैं, द्रव्यमान-त्रिज्या संबंध काफी हद तक मजबूत (robust) रहता है। हालांकि, लंबी सहसंबंध लंबाई वाले EoS को विस्तारित घनत्व श्रेणियों में सख्त रहने की अनुमति देती है, जिससे अधिकतम त्रिज्या थोड़ी बढ़ जाती है।
• फेज स्ट्रक्चर संकेतक: पोस्टीरियर नमूने उच्च घनत्वों पर शून्य के करीब पहुँचते हुए एक ट्रेस एनोमली $\Delta$ और उच्च घनत्वों पर घटते हुए एक कॉन्फॉर्मल डिस्टेंस $d_c$ दिखाते हैं, जो क्वार्क मैटर के आगमन के अनुरूप है। हालांकि, एक मजबूत प्रथम-क्रम चरण संक्रमण (दबाव के व्युत्पन्न में विच्छिन्नता) को खारिज नहीं किया जा सकता है।
• अवलोकन संबंधी बाधाएं: वैश्विक ऊष्मप्रवैगिकी सहसंबंध व्यक्तिगत अवलोकन (जैसे, PSR J0614−3329) को देखे गए सितारों के केंद्रीय घनत्वों से बहुत आगे के घनत्वों पर EoS को प्रतिबंधित करने की अनुमति देते हैं। विश्लेषण PSR J0740+6620 के त्रिज्या प्रतिबंधों और अन्य मापों के बीच तनाव को उजागर करता है।

महत्व और दावे
यह शोध दावा करता है कि बाधित GP-ब्रिज प्रायर, मौजूदा विधियों के लिए एक भौतिक रूप से प्रेरित, गणनात्मक रूप से कुशल विकल्प प्रदान करता है। पिछले दृष्टिकोणों के विपरीत, जो भौतिक रूप से स्वीकार्य व्यवहारों को बाहर कर सकते हैं या मनमाने मिलान बिंदुओं पर निर्भर कर सकते हैं, यह विधि बिना सीमा से बाहर जाए पूरे अनुमत आयतन को भरती है। लेखक इस बात पर जोर देते हैं कि उच्च घनत्वों पर देखा गया मृदुकरण (softening) वैश्विक बाधाओं (कार्य-कारणता और pcd) के कारण उत्पन्न एक स्वाभाविक परिणाम है, न कि स्थानीय सहसंबंधों का आर्टिफैक्ट। वे यह भी स्पष्ट करते हैं कि यह ढांचा ध्वनि-वेग सहसंबंध लंबाई को ट्यून करने का एक पारदर्शी तरीका प्रदान करता है और मानक GP प्रायरों के साथ सीधे तुलना की सुविधा देता है, जबकि उनके आवश्यक गैर-तुच्छ वैश्विक संरचना को बनाए रखता है। वे नोट करते हैं कि इस ढांचे की प्रतिबंध शक्ति बढ़ने की उम्मीद है जैसे-जैसे pQCD गणनाएं उच्च क्रमों (जैसे, N3LO) तक पहुँचती हैं और अनिश्चितताएं कम होती हैं।