Constraining r-process nucleosynthesis with multi-objective Galactic chemical evolution models

तकनीकी सारांश: बहु-उद्देश्यीय गैलेक्टिक रासायनिक विकास मॉडलों के साथ r-प्रक्रिया न्यूक्लियोसिंथेसिस को सीमित करना

समस्या विवरण
तीव्र न्यूट्रॉन-कैप्चर प्रक्रिया (r-प्रक्रिया) का खगोलीय उद्गम अनिश्चित बना हुआ है, जिसमें न्यूट्रॉन स्टार मर्जर्स (NSMs) और मैग्नेटो-रोटेशनल सुपरनोवा (MR-SNe) जैसे प्रतिस्पर्धी परिदृश्य शामिल हैं। जबकि मल्टी-मैसेंजर अवलोकनों ने NSMs को एक व्यवहार्य स्थल के रूप में पुष्ट किया है, उनकी विलंबित उपस्थिति और कम दरें शुरुआती, निम्न-धातुता वाले तारों में r-प्रक्रिया संवर्धन की व्याख्या करने की उनकी क्षमता को चुनौती देती हैं। इसके विपरीत, prompt स्रोत जैसे MR-SNe सैद्धांतिक रूप से प्रशंसनीय हैं लेकिन उनकी आवृत्ति और यील्ड (yields) पर प्रत्यक्ष अवलोकन संबंधी बाधाओं का अभाव है। गैलेक्टिक केमिकल इवोल्यूशन (GCE) मॉडलिंग में एक महत्वपूर्ण मुद्दा "r-प्रक्रिया सार्वभौमिकता" (r-process universality) की धारणा है—वह परिकल्पना कि यूरोपियम (Eu) से अन्य तत्वों तक r-प्रक्रिया के सापेक्ष प्रचुरता का पैटर्न सौर r-प्रक्रिया वितरण के आधार पर समान रूप से स्केल होता है। निम्न-धातुता वाले तारों के अवलोकन सुझाव देते हैं कि यह सार्वभौमिकता भारी तत्वों ($Z \ge 56$) के लिए तो बनी रह सकती है, लेकिन हल्के न्यूट्रॉन-कैप्चर तत्वों ($38 < Z < 50$) के लिए विफल हो जाती है, जो अक्सर r-प्रक्रिया-विहीन तारों में Eu के सापेक्ष उच्च अनुपात दिखाते हैं। वर्तमान मॉडल अक्सर एकल-उद्देश्यीय अनुकूलन (केवल Eu को फिट करना) पर निर्भर करते हैं, जो कई तत्वों के रासायनिक विकास को एक साथ दोहराने के लिए आवश्यक ट्रेड-ऑफ (trade-offs) को छिपा सकते हैं।

कार्यप्रणाली
लेखक लगभग $1.5 \times 10^5$ वन-इनफॉल, होमोजेनियस गैलेक्टिक केमिकल इवोल्यूशन मॉडलों के ग्रिड का उपयोग करके एक व्यवस्थित, डेटा-संचालित दृष्टिकोण अपनाते हैं। r-प्रक्रिया को किसी विशिष्ट खगोलीय स्थल (जैसे NSMs या MR-SNe) के लिए निश्चित करने के बजाय, वे एक लचीले, पैरामीट्रिक ढांचे को अपनाते हैं जिसमें चार प्रमुख मापदंडों को परिवर्तित किया जाता है:

1. प्रति घटना Eu यील्ड ($Y_{Eu,r}$): $[10^{-7}, 10^{-5}] M_\odot$ में लॉगरिदमिक रूप से सैंपल्ड।
2. घटना दर अंश ($\alpha$): उन पूर्वज तारों का अंश जो r-प्रक्रिया सामग्री उत्पन्न करते हैं, $[0, 1]$ में लीनियर रूप से सैंपल्ड।
3. विलंब समय ($\tau$): तारकीय मृत्यु और संवर्धन के बीच का समय, $[0, 14]$ Gyr में क्वाजी-लॉगारिदमिक रूप से सैंपल्ड।
4. पूर्वज द्रव्यमान सीमा ($[M_l, M_u]$): $[8, 10)$ से $[95, 100] M_\odot$ तक $5 M_\odot$ के चरणों में परिवर्तित।

अन्य न्यूट्रॉन-कैप्चर तत्वों की यील्ड को सौर r-प्रक्रिया पैटर्न (Simmerer et al. 2004) के अनुसार Eu यील्ड को स्केल करके प्राप्त किया जाता है। मॉडलों की तुलना $-2.7 \le [\text{Fe/H}] \le \text{solar}$ धातुता सीमा के विस्तृत रेंज वाले अवलोकन संबंधी डेटासेट के विरुद्ध की जाती है, जो Gaia-ESO सर्वे, MINCE, और SAGA डेटाबेस को मिलाकर निर्मित किए गए हैं।

मॉडल प्रदर्शन का मूल्यांकन करने के लिए, लेखक पारेटो फ्रंट विश्लेषण (Pareto front analysis) पर आधारित एक बहु-उद्देश्यीय अनुकूलन रणनीति का उपयोग करते हैं। एक एकल "सर्वश्रेष्ठ-फिट" मॉडल खोजने के बजाय, वे पारेटो-इष्टतम समाधानों (Pareto-optimal solutions) के सेट की पहचान करते हैं—ऐसे मॉडल जहाँ एक तत्व के फिट में सुधार करने से अनिवार्य रूप से दूसरे तत्व के फिट में गिरावट आती है। यह विभिन्न तत्वों के बीच ट्रेड-ऑफ (जैसे, Eu बनाम Zr, या हल्का समूह Sr-Y-Zr बनाम भारी समूह Ba-La-Ce) को मैप करने की अनुमति देता है।

प्रमुख योगदान और परिणाम

• पैरामीटर बाधाएं: सर्वश्रेष्ठ प्रदर्शन करने वाले मॉडल लगातार लघु विलंब समय ($\tau \le 30$ Myr), कम-द्रव्यमान वाले पूर्वजों (लगभग $20-25 M_\odot$ के आसपास केंद्रित), और प्रति घटना $\sim 2 \times 10^{-7} M_\odot$ की प्रभावी Eu इंजेक्शन दर का पक्ष लेते हैं। $\sim 80 M_\odot$ से अधिक द्रव्यमान वाले पूर्वज इतने दुर्लभ पाए जाते हैं कि वे उच्च यील्ड के बावजूद अकेले r-प्रक्रिया संवर्धन की व्याख्या करने के लिए पर्याप्त नहीं हैं।
• हल्के तत्वों के लिए सौर स्केलिंग की विफलता: जब सर्वश्रेष्ठ-Eu मॉडलों को सौर r-प्रक्रिया पैटर्न के अनुसार स्केल किया जाता है, तो वे मध्यवर्ती से निम्न धातुता पर हल्के तत्वों (Sr, Zr) की प्रचुरता को व्यवस्थित रूप से कम आंकते हैं, जबकि भारी तत्वों (Ba, La, Ce) को अधिक सफलतापूर्वक दोहराते हैं।
• पारेटो ट्रेड-ऑफ:
  • Eu बनाम Zr: द्वि-आयामी अनुकूलन में, कोई भी एकल मॉडल दोनों Eu और Zr रुझानों को एक साथ नहीं दोहरा सकता। जो मॉडल Zr को अच्छी तरह से फिट करता है, वह व्यवस्थित रूप से Eu को बढ़ाकर दिखाता है, और इसके विपरीत भी।
  • हल्का बनाम भारी समूह: हल्के समूह (Sr, Y, Zr) के लिए त्रि-आयामी अनुकूलन से पता चलता है कि इसमें महत्वपूर्ण संघर्ष हैं, जहाँ पारेटो फ्रंट इष्टतम समाधान से बड़ी दूरियां दिखाता है। इसके विपरीत, भारी समूह (Ba, La, Ce) को बहुत कम ट्रेड-ऑफ और इष्टतम बिंदु से छोटी दूरियों के साथ Eu के साथ दोहराया जा सकता है।
  • वैश्विक तुलना: जब Eu को हल्के समूह (Sr-Y-Zr) के साथ जोड़ा जाता है, तो ऑब्जेक्टिव स्पेस में इष्टतम बिंदु से दूरी काफी अधिक ($D \approx 200$) होती है, तुलनात्मक रूप से भारी समूह (Ba, La, Ce) के साथ संयोजन ($D \approx 50$) के, जो पुष्टि करता कि सौर-स्केल्ड संबंध हल्के न्यूट्रॉन-कैप्चर तत्वों की ओर बढ़ते हुए प्रगतित्मक रूप से कम वैध होता जाता है।

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

इसके बजाय, परिणाम कम से कम दो अलग-अलग स्केलिंग घटकों की आवश्यकता की ओर संकेत करते हैं:

1. एक मुख्य r-प्रक्रिया घटक (H-घटक), जो सौर पैटर्न के समान है और r-प्रक्रिया-समृद्ध तारों में देखा जाता है, जो भारी तत्वों ($Z \ge 56$) के उत्पादन में प्रमुख है।
2. एक कमजोर r-प्रक्रिया घटक (L-घटक), जिसकी विशेषता लैंथेनाइड्स के सापेक्ष हल्के तत्वों (Sr, Y, Zr) का उन्नत उत्पादन है, जो r-प्रक्रिया-विहीन तारों में देखे गए पैटर्न के समान है।

लेखक तर्क देते हैं कि ये घटक संभवतः विभिन्न खगोलीय स्थितियों या घटनाओं से उत्पन्न होते हैं, हालांकि वे विशिष्ट भौतिक स्थलों की पहचान नहीं करते हैं। यह अध्ययन सौर-स्केल्ड यील्ड धारणा के उपयोग की सीमाओं को उजागर करता है और जटिल न्यूक्लियोसिंथेटिक योगदानों को अलग करने में बहु-उद्देश्यीय अनुकूलन की उपयोगिता को प्रदर्शित करता है। भविष्य के कार्य के रूप में GCE मॉडलों के भीतर या सिम्युलेटेड एबंडेंस पैटर्न के माध्यम से कमजोर घटक को सीधे बाधित करने का सुझाव दिया गया है।