Mass distribution of neutron stars in binary systems

तकनीकी सारांश: द्विआधारी प्रणालियों में न्यूट्रॉन सितारों का द्रव्यमान वितरण

समस्या विवरण
अवलोकन संकेत देते हैं कि ज्ञात द्विआधारी प्रणालियों में न्यूट्रॉन सितारों (NSs) का द्रव्यमान वितरण एक द्वि-मोडल (bimodal) पैटर्न प्रदर्शित करता है, जिसमें 1.4 M$_\odot$ और $\sim$1.8 M$_\odot$ के आसपास शिखर (peaks) हैं। जबकि 1.4 M$_\odot$ के शिखर को अक्सर मानक कोर-कोलैप्स सुपरनोवा या सामान्य आवरण विकास (common envelope evolution) के कारण माना जाता है, $\sim$1.8 M$_\odot$ के शिखर की उत्पत्ति अभी भी बहस का विषय है। अभिवृद्धि (accretion) के माध्यम से NS द्रव्यमान वृद्धि की व्याख्या करने का प्रयास करने वाले पिछले सैद्धांतिक मॉडल अक्सर एडिंगटन सीमा (Eddington limit) द्वारा सीमित रहे हैं या उन्होंने सुपर-एडिंगटन अभिवृद्धि शासन (super-Eddington accretion regimes) का पूर्णतः लेखा-जोखा नहीं लिया है। द्रव्यमान हस्तांतरण की दक्षता, प्रोपेलर प्रभाव (propeller effect), और चुंबकीय क्षेत्र के विकास के संबंध में अनिश्चितताएं, अभिवृत्त (accreting) NSs के अंतिम द्रव्यमान वितरण की सटीक भविष्यवाणियों में बाधा डालती रही हैं। इस अध्ययन का उद्देश्य यह जांचना है कि क्या एक सुपर-एडिंगटन अभिवृद्धि मॉडल द्विआधारी प्रणालियों में NSs के देखे गए द्वि-मोडल द्रव्यमान वितरण को सफलतापूर्वक पुनरुत्पादित कर सकता है।

कार्यप्रणाली
लेखकों ने दो अलग-अलग कोडों को संयोजित करते हुए एक हाइब्रिड जनसंख्या संश्लेषण (population synthesis) दृष्टिकोण का उपयोग किया है:

1. बाइनरी स्टार इवोल्यूशन (BSE): इसका उपयोग $10^7$ आदिम द्विआधारी प्रणालियों के निर्माण और प्रारंभिक विकास के अनुकरण के लिए किया गया। इस कोड ने मुख्य-अनुक्रम (main-sequence) साथियों वाली NS प्रणालियों की एक जनसंख्या उत्पन्न की, जिसमें प्रारंभिक द्रव्यमान फलन (initial mass functions), द्रव्यमान अनुपात, उत्केंद्रता (eccentricities), और जन्मजात किक वेग (natal kick velocities) शामिल थे। उन प्रणालियों को पहचानने के लिए BSE परिणामों को फ़िल्टर किया गया जो स्थिर द्रव्यमान हस्तांतरण में सक्षम थीं, जिसमें अस्थिर सामान्य आवरण विकास या अत्यधिक द्रव्यमान हस्तांतरण दरों ($>10^4 \dot{M}_{Edd}$) वाली प्रणालियों को बाहर रखा गया।
2. मॉड्यूल्स फॉर एक्सपेरिमेंट्स इन स्टेलर एस्ट्रोफिजिक्स (MESA): BSE आउटपुट से चयनित NS द्विआधारी प्रणालियों को MESA (संस्करण r24.08.1) का उपयोग करके विस्तृत विकासवादी गणनाओं के अधीन किया गया। यह चरण स्थिर रोश लोब ओवरफ्लो (Roche lobe overflow) स्थितियों के तहत अभिवृद्धि प्रक्रिया पर केंद्रित था।

प्रमुख भौतिक मॉडल और धारणाएं:

• सुपर-एडिंगटन अभिवृद्धि: एडिंगटन सीमा द्वारा सीमित पिछले मॉडलों के विपरीत, यह अध्ययन लिपुनोव (1987) द्वारा वर्णित सुपर-एडिंगटन अभिवृद्धि के भौतिकी को शामिल करता है। मॉडल अल्फ्वेन त्रिज्या ($R_a$) और स्फेरिज़ेशन त्रिज्या ($R_s$) के बीच संबंध के आधार पर उप-क्रांतिक (sub-critical) और अति-क्रांतिक (super-critical) शासन के बीच अंतर करता है।
• अभिवृद्धि दक्षता: हस्तांतरित द्रव्यमान का वह अंश जो NS द्वारा अभिवृत्त किया जाता है ($f_{mt,a}$), अभिवृद्धि अवस्था द्वारा निर्धारित होता है। अति-क्रांतिक शासन में, अतिरिक्त सामग्री NS के पास विकिरण दबाव द्वारा उड़ा दी जाती है, जो विशिष्ट कक्षीय कोणीय संवेग को ले जाती है। रिटेन (retained) द्रव्यमान के लिए 0.5 की अभिवृद्धि दक्षता मानी गई है।
• चुंबकीय क्षेत्र का विकास: NS चुंबकीय क्षेत्र अभिवृत्त द्रव्यमान ($\Delta M$) के फलन के रूप में एक घातांकीय नियम (exponential law) का पालन करते हुए क्षय होता है, जो $5 \times 10^{12}$ G के प्रारंभिक मान से घटकर न्यूनतम $5 \times 10^8$ G हो जाता है।
• स्पिन विकास: मॉडल तीन चरणों के माध्यम से NS स्पिन अवधि को ट्रैक करता है: रेडियो (कोई अभिवृद्धि नहीं), प्रोपेलर (पदार्थ का निष्कासन), और अभिवृद्धि (स्पिन-अप)। इन चरणों के बीच संक्रमण लाइट सिलेंडर, कोरोटेशन रेडियस और मैग्नेटोस्फेरिक रेडियस की सापेक्ष स्थितियों पर निर्भर करता है।
• कोणीय संवेग की हानि: सिमुलेशन गुरुत्वाकर्षण तरंग विकिरण, चुंबकीय ब्रेकिंग, और उत्सर्जित द्रव्यमान तथा परि-द्विआधारी डिस्क (circumbinary disk) के निर्माण के कारण कोणीय संवेग की हानि को ध्यान में रखता है।

प्रमुख परिणाम

1. द्रव्यमान वृद्धि तंत्र: अध्ययन प्रदर्शित करता है कि NS द्रव्यमान वृद्धि द्विआधारी कक्षीय अवधि और डोनर (दाता) तारे के द्रव्यमान पर अत्यधिक निर्भर है। सुपर-एडिंगटन मॉडल में, NS दीर्घ समय तक अभिवृद्धि संतुलन अवस्था में रह सकते हैं, जिससे महत्वपूर्ण द्रव्यमान संचय (1.4 M$_\odot$ जन्म द्रव्यमान से $\sim$1.0 M$_\odot$ तक की वृद्धि) संभव होता है। यह ली एट अल. (2021) के मॉडल (एडिंगटन-सीमित) के विपरीत है, जहाँ NS अक्सर एक लंबे प्रोपेलर चरण में प्रवेश करते हैं जो द्रव्यमान वृद्धि को सीमित करता है।
2. द्वि-मोडालिटी का पुनरुत्पादन: जनसंख्या संश्लेषण के परिणाम देखे गए द्वि-मोडल वितरण को सफलतापूर्वक पुनरुत्पादित करते हैं:
   • $\sim$1.8 M$_\odot$ का शिखर: यह शिखर मुख्य रूप से उन NS द्विआधारी प्रणालियों से उत्पन्न होता है जिनमें डोनर द्रव्यमान $M_d < \sim$1.6 M$_\odot$ और कक्षीय अवधि $P_{orb} < 20$ दिन है। इन प्रणालियों में द्रव्यमान हस्तांतरण दर लंबी अवधि में कुशल, अर्ध-संतुलन उप-एडिंगटन अभिवृद्धि की अनुमति देती है।
   • 1.4 M$_\text{odot}$ का शिखर: यह शिखर उन NS प्रणालियों के कारण है जो सामान्य आवरण विकास (अक्सर अधिक द्रव्यमान वाले डोनर्स या विशिष्ट कक्षीय विन्यास शामिल) से गुजरते हैं, जहाँ द्रव्यमान वृद्धि न्यूनतम होती है।
   • उच्च-द्रव्यमान पूँछ (High-Mass Tail): मध्यवर्ती द्रव्यमान वाले डोनर्स ($1.6 \le M_d \le 3.0$ M$_\odot$) और लघु कक्षीय अवधि वाली प्रणालियाँ 2.1 M$_\odot$ से अधिक द्रव्यमान अभिवृत्त कर सकती हैं।
3. स्पिन अवधि वितरण: मॉडल लो-मास एक्स-रे बाइनरीज (LMXBs) के लिए 1 ms और 15 ms के बीच केंद्रित NS स्पिन अवधि की भविष्यवाणी करता है। परिणाम कुछ प्रेक्षित डेटा की तुलना में थोड़े बड़े हैं, जिसे लेखक $5 \times 10^8$ G के न्यूनतम चुंबकीय क्षेत्र की धारणा के कारण बताते हैं, जो अभिवृद्धि डिस्क को NS की सतह के बहुत करीब जाने से रोकता है।

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

सीमाएं और अनिश्चितताएं
लेखक अपने मॉडल में कई अनिश्चितताओं को स्वीकार करते हैं:

• अभिवृद्धि दक्षता: माना गया 0.5 का दक्षता स्तर और NS के पास द्रव्यमान हानि का उपचार सरलीकृत है।
• चुंबकीय क्षेत्र का क्षय: क्षय दर और न्यूनतम क्षेत्र शक्ति वे पैरामीटर हैं जो स्पिन विकास और अभिवृद्धि दक्षता को महत्वपूर्ण रूप से प्रभावित करते हैं।
• इक्वेशन ऑफ स्टेट (EoS): जड़त्व आघूर्ण (moment of inertia) को $10^{45}$ g cm$^2$ पर स्थिर रखा गया है, जो निचले स्तर के करीब है। NS EoS में परिवर्तन जड़त्व आघूर्ण को बदल सकता है, जिससे स्पिन-अप समय और प्रोपेलर शासन के संक्रमण पर प्रभाव पड़ेगा।
• बाहर रखे गए चैनल: अध्ययन यह मात्रात्मक निर्धारण नहीं करता है कि 'एक्रीशन-इंड्यूस्ड कोलैप्स' (AIC) के माध्यम से बने NS या द्वि-मोडल द्रव्यमान वितरण के साथ पैदा हुए NS का कितना योगदान है, बल्कि यह मुख्य रूप से मानक कोर-कोलैप्स और उसके बाद के द्विआधारी विकास के माध्यम से बने NS पर केंद्रित है।