Non-LTE Ionization Modeling for Helium and Strontium in Neutron Star Merger Ejecta

यहाँ चिबा एवं अन्य (Chiba et al., ड्राफ्ट संस्करण 8 अप्रैल, 2026) के शोध पत्र "Non-LTE Ionization Modeling for Helium and Strontium in Neutron Star Merger Ejecta" का विस्तृत तकनीकी सारांश दिया गया है।

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

बाइनरी न्यूट्रॉन स्टार (BNS) विलय GW170817 से जुड़ी किलोनोवा AT2017gfo, अपने प्रारंभिक-चरण के स्पेक्ट्रा में 1 µm के आसपास एक प्रमुख अवशोषण विशेषता (absorption feature) प्रदर्शित करती है। इस विशेषता की व्याख्या करने के लिए दो प्राथमिक उम्मीदवारों का प्रस्ताव दिया गया है:

• स्ट्रोंटियम (Sr II): विशेष रूप से 1.03–1.09 µm के आसपास का संक्रमण।
• हीलियम (He I): विशेष रूप से 1.08 µm ($2^3S - 2^3P$) पर संक्रमण।

पिछले अध्ययन इन तत्वों के योगदान को समान रूप से निर्धारित करने में संघर्ष करते रहे हैं क्योंकि मुख्य रूप से दो सीमाएँ थीं:

1. LTE धारणा (LTE Assumption): प्रारंभिक-चरण के स्पेक्ट्रल मॉडलिंग में काफी हद तक लोकल थर्मोडायनामिक इक्विलिब्रियम (LTE) पर निर्भरता रही है। हालाँकि, किलोनोवा इजेक्टा (ejecta) का निम्न घनत्व और उच्च-ऊर्जा विकिरण वातावरण (जो रेडियोधर्मी क्षय द्वारा संचालित होता है) महत्वपूर्ण गैर-तापीय आयनीकरण (non-thermal ionization) उत्पन्न करता है, जिससे LTE से विचलन होता है।
2. समान मॉडलिंग का अभाव (Lack of Simultaneous Modeling): पिछले कार्यों में या तो He या Sr का अलग-अलग मॉडल बनाया गया, या वे देखे गए फीचर को पुनरुत्पादित करने के लिए आवश्यक द्रव्यमान अंशों (mass fractions) की पूरी श्रृंखला का व्यवस्थित रूप से अन्वेषण नहीं कर सके।

मूल समस्या एक स्व-संगत नॉन-एलटीई (Non-LTE) आयनीकरण मॉडल विकसित करके इजेक्टा में He और Sr की वास्तविक प्रचुरता (abundances) निर्धारित करना है, जो उच्च-ऊर्जा इलेक्ट्रॉनों द्वारा आयनीकरण को ध्यान में रखता है।

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

A. स्पेक्ट्रल फिटिंग और भौतिक स्थितियाँ (Spectral Fitting and Physical Conditions)

लेखकों ने पहले चार युगों (विलय के 1.43, 2.42, 3.41, और 4.40 दिन बाद) में VLT/X-shooter स्पेक्ट्रा का विश्लेषण किया। उन्होंने निरंतरता (continuum) और 1 µm फीचर को एक रिलेटिविस्टिक ब्लैकबॉडी और एक P-Cygni प्रोफाइल का उपयोग करके फिट किया।

• मुख्य निष्कर्ष: 1 µm फीचर की गहराई को पुनरुत्पादित करने के लिए, फोटोस्फीयर पर एक सोबोलेव ऑप्टिकल डेप्थ (Sobolev optical depth) $\tau \approx 1$ (विशेष रूप से $0.5 \le \tau \le 2$) की आवश्यकता होती है।
• पैरामीटर: फोटोस्फेरिक वेग ($v_{ph}$) $0.21c$ से $0.14c$ के बीच हैं, और तापमान ($T_0$) $4500$K से$2600$ K के बीच है।

B. नॉन-एलटीई आयनीकरण ढांचा (Non-LTE Ionization Framework)

अध्ययन दोनों He और Sr के लिए एक दर-समीकरण सॉल्वर (rate-equation solver) विकसित करता है, जिसमें रेडियोधर्मी क्षय ($\beta$-decay और $\gamma$-rays) द्वारा उत्पन्न उच्च-ऊर्जा इलेक्ट्रॉनों के कारण होने वाले गैर-तापीय आयनीकरण को शामिल किया गया है।

• हीलियम (He) मॉडल:

  • 21 अवस्थाओं वाला एक विस्तृत परमाणु मॉडल (He I की 19 बाउंड अवस्थाएं $n=4$ तक, साथ ही He II और He III की ग्राउंड स्टेट्स)।
  • इसमें उत्तेजना (excitation), डी-एक्साइटेशन, आयनीकरण और पुनर्संयोजन (recombination) के लिए रेडिएटिव और कोलिजनल (इलेक्ट्रॉन-इम्पैक्ट) प्रक्रियाओं को शामिल किया गया है।
  • महत्वपूर्ण तंत्र: 1.08 µm संक्रमण मेटास्टेबल $2^3S$ अवस्था से उत्पन्न होता है। विशिष्ट किलोनोवा तापमान ($<5000$ K) में, तापीय उत्तेजना इस अवस्था को कुशलतापूर्वक भरने में सक्षम नहीं होती है। मॉडल गैर-तापीय आयनीकरण और उसके बाद पुनर्संयोजन पर निर्भर करता है ताकि $2^3S$ अवस्था को भरा जा सके।

• स्ट्रोंटियम (Sr) मॉडल:

  • भारी तत्वों की जटिलता के कारण, एक सरलीकृत आयनीकरण मॉडल का उपयोग किया गया है, जो Sr XX तक के आयनीकरण अवस्थाओं को ट्रैक करता है (हालाँकि केवल Sr VII तक ही महत्वपूर्ण हैं)।
  • Sr II की उत्तेजित अवस्थाओं के लिए बोल्ट्ज़मैन वितरण (Boltzmann distribution) मानता है।
  • आयनीकरण संतुलन थर्मल फोटोआयनीकरण (diluted blackbody) और गैर-तापीय आयनीकरण ($\Gamma_{nt} = \dot{q}/w_i$) के बीच प्रतिस्पर्धा द्वारा नियंत्रित होता है।
  • मुख्य तंत्र: उच्च-ऊर्जा इलेक्ट्रॉन ओवर-आयनीकरण (over-ionization) का कारण बनते हैं, जिससे Sr II से उच्च आयनीकरण अवस्थाओं (Sr III, IV, आदि) में इलेक्ट्रॉन निकल जाते हैं, जिससे Sr II की जनसंख्या कम हो जाती है।

C. प्रचुरता प्रतिबंध (Abundance Constraints)

लेखक प्रत्येक युग में लक्षित ऑप्टिकल डेप्थ ($\tau \approx 1$) प्राप्त करने के लिए दर समीकरणों को पुनरावृत्ति (iteratively) से हल करते हैं, जिसमें आयनीकरण अवस्थाओं से प्राप्त स्व-संगत इलेक्ट्रॉन घनत्व को भी ध्यान में रखा जाता है।

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

1. एक स्व-संगत नॉन-एलटीई मॉडल का विकास: यह दोनों He और Sr के आयनीकरण को एक साथ मॉडल करने वाला पहला अध्ययन है, जो उच्च-ऊर्जा इलेक्ट्रॉनों द्वारा गैर-तापीय आयनीकरण को स्पष्ट रूप से शामिल करता है।
2. नॉन-एलटीई प्रभावों का परिमाणीकरण:
   • Sr के लिए: ओवर-आयनीकरण के कारण LTE भविष्यवाणियों की तुलना में Sr II की जनसंख्या को नॉन-एलटीई प्रभाव महत्वपूर्ण रूप से दबा देते हैं (कारक $10^2$–$10^4$ द्वारा)। फलस्वरूप, फीचर को पुनरुत्पादित करने के लिए आवश्यक Sr का द्रव्यमान अंश LTE की तुलना में बहुत अधिक है।
   • He के लिए: मेटास्टेबल $2^3S$ अवस्था की जनसंख्या तापमान और गैर-तापीय आयनीकरण दरों के प्रति अत्यधिक संवेदनशील है।
3. व्यवस्थित प्रचुरता प्रतिबंध: शोध पत्र केवल एकल प्रतिनिधि मानों के बजाय, इजेक्टा की वेग सीमा ($v \sim 0.15c$) में द्रव्यमान अंशों ($X_{He}$ और $X_{Sr}$) पर प्रतिबंध प्रदान करता है।

4. परिणाम (Results)

A. आवश्यक द्रव्यमान अंश (Required Mass Fractions)

AT2017gfo में $v \sim 0.15c$ पर 1 µm फीचर को पुनरुत्पादित करने के लिए:

• हीलियम परिदृश्य (Helium Scenario): इसके लिए $X_{He} \approx 1\%$ के द्रव्यमान अंश की आवश्यकता है।
• स्ट्रोंटियम परिदृश्य (Strontium Scenario): इसके लिए $X_{Sr} \approx 1\% - 10\%$ के द्रव्यमान अंश की आवश्यकता है।
  • नोट: नॉन-एलटीई मामले में आवश्यक Sr द्रव्यमान अंश, LTE मॉडलों की तुलना में लगभग दो परिमाण (two orders of magnitude) अधिक है।
  • प्राप्त Sr द्रव्यमान अंश (1-10%) सौर r-प्रक्रिया प्रचुरता (solar r-process abundance) पैटर्न के साथ उल्लेखनीय रूप से मेल खाता है।

B. न्यूक्लियोसिंथेसिस निहितार्थ (Nucleosynthesis Implications)

पैरामेट्रिक न्यूक्लियोसिंथेसिस गणनाओं (S. Wanajo 2018) के साथ अपने प्राप्त प्रचुरता की तुलना करके, लेखक इजेक्टा की भौतिक स्थितियों को सीमित करते हैं:

• इलेक्ट्रॉन अंश (Electron Fraction - $Y_e$): परिणाम निम्न इलेक्ट्रॉन अंश ($Y_e \lesssim 0.35$) का समर्थन करते हैं। उच्च $Y_e$ ($>0.35$) $\alpha$-रिच फ्रीज-आउट की ओर ले जाता है, जिससे He और Sr दोनों का अत्यधिक उत्पादन होता है, जो एक ऐसा अवशोषण फीचर बनाएगा जो प्रेक्षणों से बहुत अधिक मजबूत होगा।
• एन्ट्रॉपी (Entropy - $s$): परिणाम निम्न एन्ट्रॉपी ($s \lesssim 30 k_B/\text{nucleon}$) का समर्थन करते हैं। उच्च एन्ट्रॉपी भी $\alpha$-रिच फ्रीज-आउट को बढ़ा देती है।
• 1 µm फीचर के दो संभावित मूल:
  1. मध्यवर्ती $Y_e$ ($0.15 \lesssim Y_e \lesssim 0.35$): यह फीचर Sr द्वारा प्रधान है। यह इस परिकल्पना का समर्थन करता है कि बाइनरी न्यूट्रॉन स्टार (BNS) ब्रह्मांड में r-प्रक्रिया तत्वों के प्राथमिक स्रोत हैं।
  2. अत्यंत निम्न $Y_e$ ($Y_e \lesssim 0.15$): यह फीचर He द्वारा प्रधान है। इस शासन में, He का उत्पादन trans-Pb नाभिकों के $\alpha$-क्षय के माध्यम से होता है। यदि इसकी पुष्टि होती है, तो यह तीसरे r-प्रक्रिया शिखर (third r-process peak) से परे तत्वों के उत्पादन का एक अप्रत्यक्ष संकेत प्रदान करेगा।

C. द्रव्यमान उत्सर्जन तंत्र (Mass Ejection Mechanism)

अपने प्रतिबंधों की तुलना दीर्घकालिक विस्कस हाइड्रोडायनामिक्स सिमुलेशन (विशेष रूप से DD2-135-135 मॉडल, जो एक लंबे समय तक जीवित रहने वाले न्यूट्रॉन स्टार अवशेष को छोड़ता है) से करके:

• मानक लंबे समय तक जीवित रहने वाले अवशेष मॉडल He और Sr दोनों का अत्यधिक उत्पादन (overproduce) करते हैं क्योंकि उनके पोस्ट-मर्जर विंड्स में उच्च $Y_e$ और उच्च एन्ट्रॉपी होती है।
• लेखक सुझाव देते हैं कि GW170817 में महत्वपूर्ण द्रव्यमान उत्सर्जन वर्तमान लंबे समय तक जीवित रहने वाले मॉडलों की तुलना में कम $Y_e$ और कम एन्ट्रॉपी वाली स्थितियों में हुआ होगा।
• यह संकेत देता है कि द्रव्यमान उत्सर्जन एक छोटे समय अंतराल (shorter timescale) पर होना चाहिए (इससे पहले कि $Y_e$ बहुत अधिक बढ़ जाए), जो हालिया मैग्नेटोहाइड्रोडायनामिक्स (MHD) सिमुलेशन के अनुरूप है।

5. महत्व (Significance)

• 1 µm फीचर का समाधान: यह अध्ययन स्पष्ट करता है कि 1 µm फीचर को या तो Sr या He द्वारा समझाया जा सकता है, लेकिन आवश्यक प्रचुरता आयनीकरण भौतिकी के आधार पर नाटकीय रूप से भिन्न होती है। नॉन-एलटीई उपचार अनिवार्य है; LTE मॉडल आवश्यक Sr प्रचुरता का बहुत कम अनुमान लगाते हैं।
• कॉस्मोकेमिकल प्रतिबंध: प्राप्त Sr प्रचुरता और सौर r-प्रक्रिया पैटर्न के बीच का सामंजूल्य इस बात को मजबूत करता है कि BNS विलय भारी तत्वों के उत्पादन के प्रमुख स्थल हैं।
• इक्वेशन ऑफ स्टेट (EoS) की जांच: He की प्रचुरता अवशेष न्यूट्रॉन स्टार के जीवनकाल के लिए एक संभावित प्रोब (probe) के रूप में कार्य करती है। यदि यह फीचर He-प्रधान है, तो यह बहुत कम $Y_e$ वातावरण का संकेत देता है, जो संभावित रूप से परमाणु EoS और अवशेष की स्थिरता को सीमित करता है।
• भविष्य की दिशाएं: यह शोध पत्र He उत्सर्जन रेखाओं (जैसे 2.06 µm) के अधिक विस्तृत मॉडलिंग और He तथा Sr परिदृश्यों के बीच के द्वंद्व (degeneracy) को तोड़ने के लिए सुसंगत लाइट कर्व मॉडलिंग की आवश्यकता पर बल देता है।

निष्कर्षतः, यह कार्य स्थापित करता है कि किलोनोवा स्पेक्ट्रा की व्याख्या के लिए नॉन-एलटीई आयनीकरण अत्यंत महत्वपूर्ण है। यह GW170817 के इजेक्टा की संरचना पर मजबूत प्रतिबंध प्रदान करता है, जो एक निम्न-एन्ट्रॉपी, निम्न-इलेक्ट्रॉन-अंश वाले वातावरण का सुझाव देता है जो वर्तमान लंबे समय तक जीवित रहने वाले अवशेष मॉडलों से भिन्न हो सकता है।