On a re-examination of neutron star cooling in transient sources -- No shallow heating required?

तकनीकी सारांश: क्षणिक स्रोतों (transient sources) में न्यूट्रॉन तारे के शीतलन (cooling) के पुनरीक्षण पर

समस्या विवरण
एक अभिवृद्धि विस्फोट (accretion outburst) के बाद लो-मास एक्स-रे बाइनरीज़ (LMXBs) में न्यूट्रॉन सितारों (NSs) के तापीय विकास (thermal evolution) के लिए पारंपरिक रूप से "उथले हीटिंग" (shallow heating) नामक एक अतिरिक्त, एड-हॉक (ad hoc) हीटिंग स्रोत का आह्वान करने की आवश्यकता होती है, जिसका भौतिक मूल अज्ञात है। इस तंत्र को आमतौर पर द्रव्यमान अभिवृद्धि दर ($\dot{M}$) के समानुपाती माना जाता है और यह विस्फोट के कुछ सौ दिनों के बाद देखी जाने वाली उच्च प्रभावी तापमान ($T_{\infty}^{\text{eff}}$) को दर्शाने के लिए आवश्यक होता है। मानक मॉडलिंग यह मानती है कि न्यूट्रॉन स्टार का लिफाफा (envelope) $\dot{M}$ और निरंतर थर्मोन्यूक्लियर दहन (thermonuclear burning) से स्वतंत्र एक सीमा स्थिति (boundary condition) के रूप में कार्य करता है, जो लिफाफे में उत्पन्न ऊष्मा के क्रस्ट (crust) में रिसाव की संभावना की उपेक्षा करता है।

कार्यप्रणाली
लेखक न्यूट्रॉन स्टार के क्रस्ट और कोर के तापीय विकास का अनुकरण करने के लिए पूर्णतः सापेक्षतावादी (fully relativistic) 1D कोड NSCool का उपयोग करते हैं। यह अध्ययन लिफाफे/क्रस्ट संक्रमण (सीमा घनत्व $\rho_b \approx 10^7$ g cm$^{-3}$) पर एक नवीन सतह सीमा स्थिति (BC) पेश करता है। यह "New BC", जो पिछले मॉडलों के विपरीत है जो एक स्थिर $T_b - T_{\text{eff}}$ संबंध का उपयोग करते हैं, $\dot{M}$ के प्रति संवेदनशील है और लिफाफे से क्रस्ट में होने वाले थर्मोन्यूक्लियर हीटिंग लीकेज को ध्यान में रखता है, जो "वॉटरशेड मॉडल" (Hanawa & Fujimoto 1984; Fujimoto et al. 1987) के अनुरूप है।

अध्ययन सात विशिष्ट स्रोतों से आठ विस्फोट प्रकरणों (outburst episodes) का विश्लेषण करता है: MXB 1659-29 (दो विस्फोट), XTE J1701-462, EXO 0748-676, IGR J17480-2446 (Ter5 X-2), Swift J174805.3-244637 (Ter5 X-3), 1RXS J180408.9-342058, और KS 1731-260। इन स्रोतों को उनकी विविध विस्फोट अवधि ($t_{\text{acc}} \sim 0.1$ से 24 वर्ष) और अभिवृद्धि दरों ($10^{16}$ से $10^{18}$ g s$^{-1}$) के आधार पर चुना गया था।

दो फिटिंग परिदृश्यों (fitting scenarios) का मोंटे कार्लो मार्कोव चेन (MCMC) विधियों द्वारा परीक्षण किया गया:

• परिदृश्य A (Scenario A): उथले हीटिंग के बिना शीतलन मॉडल। मुक्त पैरामीटर $\dot{M}$, हल्के तत्वों का कॉलम डेप्थ ($y_L$), रेड-शिफ्टेड कोर तापमान ($\tilde{T}_{\text{core}}$), और दो अशुद्धि कारक (impurity factors) ($Q_{\text{imp},1}$ बाहरी क्रस्ट के लिए, $Q_{\text{imp},2}$ आंतरिक क्रस्ट के लिए) हैं।
• परिदृश्य B (Scenario B): उथले हीटिंग सहित एक नल टेस्ट (null test)। यह तीन पैरामीटर जोड़ता है: हीटिंग स्ट्रेंथ ($Q_{\text{sh}}$) और वह घनत्व सीमा $[\rho_{\text{sh,min}}, \rho_{\text{sh,max}}]$ जहाँ यह कार्य करता है।

मुख्य परिणाम

1. परिदृश्य A की सफलता: नया सीमा प्रतिबंध (boundary condition) बिना किसी उथले हीटिंग की आवश्यकता के सात में से आठ स्रोतों के शीतलन वक्रों (cooling curves) को सफलतापूर्वक पुनरुत्पादित करता है। यह मॉडल विस्फोट के तुरंत बाद उच्च $T_{\infty}^{\text{eff}}$ को स्वाभाविक रूप से समझाता है।
2. पैरामीटरों की निरंतरता: परिदृश्य A में, फिट किए गए $\dot{M}$ मान आम तौर पर प्रेक्षण योग्य औसत अभिवृद्धि दरों के साथ अच्छी तरह मेल खाते हैं (Ter5 X-2 को छोड़कर, जहाँ फिट के लिए दर 3-4 गुना अधिक आवश्यक है)। अशुद्धि कारक (impurity factors) लगातार गहराई के साथ घटते रुझान का समर्थन करते हैं ($Q_{\text{imp},1} > Q_{\text{imp,2}}$), जिनका परिमाण बाहरी क्रस्ट में $\sim 10$ और आंतरिक क्रस्ट में $\sim 1$ है। हल्के तत्व का कॉलम डेप्थ लगातार $y_L \sim 10^8$ g cm$^{-2}$ के आसपास पाया गया है।
3. परिदृश्य B की विफलता (भौतिक तर्कसंगतता): हालांकि परिदृश्य B डेटा को फिट कर सकता है, लेकिन इसके लिए अवास्तविक (unphysical) पैरामीटरों की आवश्यकता होती है। विशेष रूप से, यह फिट किए गए $\dot{M}$ को प्रेक्षण योग्य अनुमानों से एक क्रम (order of magnitude) कम करने के लिए मजबूर करता है, जिसकी भरपाई अत्यधिक उच्च उथले हीटिंग स्ट्रेंथ ($Q_{\text{sh}} \gtrsim 5$ MeV baryon$^{-1}$) द्वारा की जाती है।
4. मॉडल चयन: सूचना मानदंड (AICc और BIC) 7 में से 5 स्रोतों के लिए परिदृश्य A का पक्ष लेते हैं। EXO 0748-676 के लिए, परिदृश्य B सांख्यिकीय रूप से पसंदीदा है, लेकिन आवश्यक उथला हीटिंग स्ट्रेंथ अभी भी अवास्तविक रूप से बड़ा है। EXO 0748-676 का शीतलन वक्र एकमात्र ऐसा मामला है जहाँ परिदृश्य A संतोषजनक फिट प्रदान करने में विफल रहता है, जो संभवतः इस स्रोत के उच्च द्रव्यमान ($2M_{\odot}$) और लंबी विस्फोट अवधि के कारण है।

महत्व और दावे
यह शोध पत्र तर्क देता है कि "उथला हीटिंग" (shallow heating), जिसे पहले न्यूट्रॉन स्टार (NS) शीतलन मॉडलों में बुलाया जाता था, संभवतः सरलीकृत सतह सीमा स्थितियों का एक आर्टिफैक्ट (artifact) है जो अभिवृद्धि दर, लिफाफे में थर्मोन्यूक्लियर दहन और अंतर्निहित क्रस्ट के बीच जुड़ाव की उपेक्षा करती है। एक ऐसी सीमा स्थिति को लागू करके जो लिफाफे से क्रस्ट में ऊर्जा के प्रवाह की अनुमति देती है (वॉटरशेड मॉडल), लेखक यह प्रदर्शित करते हैं कि क्षणिक स्रोतों के देखे गए शीतलन व्यवहार को बिना किसी एड-हॉक हीटिंग तंत्र के मानक भौतिकी द्वारा समझाया जा सकता है।

लेखक निष्कर्ष निकालते हैं कि न्यूट्रॉन स्टार के लिफाफे में थर्मोन्यूक्लियर हीटिंग, शीतलन वक्रों को आकार देने में पहले की तुलना में अधिक महत्वपूर्ण भूमिका निभाता है। हालांकि वे स्वीकार करते हैं कि सुधार संभव हैं—जैसे कि समय-निर्भर अभिवृद्धि दर, परिवर्तनशील NS द्रव्यमान/त्रिज्या, और गहराई के फलन के रूप में अशुद्धियों का मॉडलिंग करना—उनके परिणाम दृढ़ता से सुझाव देते हैं कि कई, यदि नहीं तो सभी, क्षणिक स्रोतों के लिए उथले हीटिंग की आवश्यकता नहीं हो सकती है। EXO 0748-676 का अपवाद यह सुझाव देता है कि विशिष्ट उच्च-द्रव्यमान या लंबी अवधि वाले सिस्टम के लिए अतिरिक्त भौतिकी की अभी भी आवश्यकता हो सकती है।