Runaway OB stars within 1 kpc of the Sun

तकनीकी सारांश: सूर्य के 1 kpc के भीतर भागने वाले (Runaway) OB तारे

समस्या विवरण
भागने वाले तारे (Runaway stars) उच्च-वेग वाली वे वस्तुएं हैं जो अपने जन्म के वातावरण से बाहर फेंकी जाती हैं, जो तारकीय समूहों के गतिक इतिहास और तारकीय निष्कासन को संचालित करने वाले तंत्रों के बारे में महत्वपूर्ण अंतर्दृष्टि प्रदान करती हैं। हालांकि यह शब्द पारंपरिक रूप से OB तारों के लिए लागू होता है क्योंकि वे अवलोकन योग्य होते हैं, "रनवे अंश" (OB तारों का वह अनुपात जो रनवे तारे हैं) के अनुमानों ने ऐतिहासिक रूप से व्यापक प्रसार दिखाया है, जो ~1.5% से लेकर 50% से अधिक तक रहा है। यह विसंगति मुख्य रूप से अंतर्निहित तारकीय नमूनों (जो अक्सर परिमाण-सीमित या स्पेक्ट्रोस्कोपिक रूप से पक्षपाती होते हैं) और रनवे तारों को परिभाषित करने और चुनने के विविध तरीकों के कारण होती है। पिछले अध्ययन रनवे तारों का एक एकीकृत, वॉल्यूम-पूर्ण (volume-complete) गणना प्रदान करने में संघर्ष करते रहे हैं ताकि इन विसंगतियों को सुलझाया जा सके।

कार्यप्रणाली
यह अध्ययन सूर्य के 1 kpc के भीतर OB तारों के एक हाल ही में प्रकाशित, वॉल्यूम-पूर्ण नमूने का उपयोग करता है, जो क्विंटाना एट अल. (2025) के कैटलॉग से प्राप्त किया गया है। अंतिम नमूने में 40 O-प्रकार के तारे और 24,488 B-प्रकार के तारे शामिल हैं। विश्लेषण निम्नलिखित चरणों के माध्यम से आगे बढ़ता है:

1. गैलेक्टिक रोटेशन मॉडलिंग: लेखक MCMC दृष्टिकोण का उपयोग करके नमूने के प्रेक्षित 'प्रॉपर मोशन' (proper motions) पर एक गैलेक्टिक रोटेशन मॉडल फिट करते हैं। मॉडल गैलेक्टिक त्रिज्या और देशांतर के फलन के रूप में गैलेक्टोसेंट्रिक वेग को मानता है, जो त्रिज्या के सापेक्ष अज़ीमुथल वेग ($V_{\phi}$) और रेडियल वेग ($V_R$) तथा उनके प्रथम डेरिवेटिव को पैरामीटराइज़ करता है। मॉडल गैलेक्टिक सेंटर के सापेक्ष सूर्य की गति को भी ध्यान में रखता है।
2. अवशिष्ट वेग गणना (Residual Velocity Calculation): मॉडल द्वारा अनुमानित वेगों को प्रेक्षित वेगों से घटाकर, लेखक गैलेक्टिक प्लेन ($\ell$) और वर्टिकल ($b$) दिशाओं में इन तारों के विशिष्ट वेग (peculiar velocities) प्राप्त करते हैं।
3. सांख्यिकीय लक्षण वर्णन: इन अवशिष्ट वेगों के वितरण को गैर-रनवे आबादी की गतिक विखराव (kinematic dispersion) को चित्रित करने के लिए मैक्सवेल-बोल्ट्ज़मैन वितरण (Maxwell-Boltzmann distribution) पर फिट किया जाता है।
4. रनवे चयन मानदंड: रनवे तारों की पहचान करने के लिए दो अलग-अलग विधियां नियोजित की गई हैं:
   • विधि 1 (निश्चित थ्रेशोल्ड): एक निश्चित 2D विशिष्ट वेग थ्रेशोल्ड $23 \text{ km s}^{-1}$ लागू किया जाता है। यह मान प्रेक्षित संभाव्यता घनत्व और मैक्सवेल-बोल्ट्ज़मैन फिट के प्रतिच्छेदन पर आधारित है, जो उस वेग का प्रतिनिधित्व करता है जहाँ एक तारे के रनवे होने की संभावना अधिक होती है।
   • विधि 2 (सामान्यीकृत विखराव): तारों को तब रनवे के रूप में पहचाना जाता है यदि उनके अनिश्चितता-सामान्यीकृत विशिष्ट वेग एक महत्व थ्रेशोल्ड ($E > 1$) से अधिक होते हैं, जो प्रभावी रूप से उन तारों को चुनता है जिनके वेग उनके मापन त्रुटियों के सापेक्ष नमूना विखराव से काफी अधिक हैं।
5. मोंटे कार्लो सिमुलेशन: मापन अनिश्चितताओं और मॉडल पैरामीटर विविधताओं को ध्यान में रखने के लिए, 500 सिमुलेशन का एक मोंटे कार्लो प्रयोग चलाया जाता है। एक तारे के रनवे होने की संभावना इस आधार पर गिनी जाती है कि वह सिमुलेशन के दौरान कितनी बार चयन मानदंडों को पार करता है।
6. 3D अंशांकन (Calibration): 3D विशिष्ट वेगों की गणना करने और 3D मैक्सवेल-बोल्ट्ज़मैन फिट के विरुद्ध 2D थ्रेशोल्ड को सत्यापित करने के लिए, उपलब्ध रेडियल वेगों (APOGEE और LAMOST से) वाले 3,884 तारों के उपसमूह का उपयोग करके 2D परिणामों को अंशांकित किया जाता है।

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

• रनवे अंश (Runaway Fractions): निश्चित वेग थ्रेशोल्ड (विधि 1) का उपयोग करते हुए, अध्ययन O-प्रकार के तारों के लिए 17.5% ($^{+0.1}_{-2.5}$) और B-प्रकार के तारों के लिए 7.0% ($\pm 0.1$) का रनवे अंश पाता है। विधि 2 निम्न अंश प्रदान करती है: O-प्रकार के लिए 10.0% और B-प्रकार के लिए 2.5%।
• वेग वितरण: नमूने के अवशिष्ट वेगों को मैक्सवेल-बोल्ट्ज़मैन वितरण द्वारा अच्छी तरह से वर्णित किया गया है, जिसमें $b \approx 11 \text{ km s}^{-1}$ (2D) और $b \approx 14 \text{ km s}^{-1}$ (3D) है। 3D विश्लेषण पुष्टि करता है कि 2D थ्रेशोल्ड $23 \text{ km s}^{-1}$, पूर्ण 3D उपचार के लिए एक उचित सन्निकटन है, जो लगभग $28 \text{ km s}^{-1}$ के 3D थ्रेशोल्ड के अनुरूप है।
• व्यक्तिगत उम्मीदवार: यह अध्ययन प्रॉपर मोशन के आधार पर कई प्रसिद्ध उम्मीदवारों (जैसे, $\mu$ Col, AE Aur, $\zeta$ Pup, 68 Cyg, $\lambda$ Cep) की रनवे स्थिति की पुष्टि करता है। यह उन तारों की भी पहचान करता है जिन्हें केवल रेडियल वेग के आधार पर रनवे माना गया था (जैसे, $\xi$ Per, X Per) और उन्हें इन मानदंडों के तहत रनवे के रूप में वर्गीकृत करने के लिए अपर्याप्त 2D प्रॉपर मोशन के रूप में चिह्नित करता है, जो विशिष्ट वस्तुओं के लिए 2D-मात्र विश्लेषण की सीमाओं को उजागर करता है।
• स्थानिक वितरण: रनवे तारे, विशेष रूप से B-प्रकार के, सामान्य OB आबादी की तुलना में अधिक विसरित (diffuse) स्थानिक वितरण दिखाते हैं। O-प्रकार के रनवे तारे गैलेक्टिक प्लेन के करीब रहते हैं, जो संभवतः उनकी कम आयु के कारण है। रनवे तारों का ऊर्ध्वाधर वितरण सामान्य OB आबादी की तुलना में "थिन डिस्क" घटक के लिए कम भार दिखाता है, जो यह संकेत देता है कि रनवे तारे उच्च गैलेक्टिक अक्षांशों पर पाए जाने की अधिक संभावना रखते हैं।

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

परिणाम दर्शाते हैं कि O-प्रकार के तारों का रनवे अंश B-प्रतकार के तारों की तुलना में काफी अधिक है, एक ऐसा रुझान जो पिछले साहित्य में संकेत दिया गया था लेकिन अब एक पूर्ण स्थानीय नमूने के साथ मात्रात्मक रूप से सिद्ध किया गया है। लेखक सुझाव देते हैं कि O-प्रकार के तारों का उच्च रनवे अंश उनके निष्कासन के लिए प्रमुख तंत्र के रूप में डायनामिकल इजेक्शन परिदृश्य (DES) का समर्थन करता है, जबकि B-प्रकार के तारों के लिए कम अंश बाइनरी सुपरनोवा परिदृश्य (BSS) के बड़े योगदान को इंगित कर सकता है। शोध निष्कर्ष निकालता है कि हालांकि प्राप्त अंश पिछले कुछ अध्ययनों के अनुरूप हैं, साहित्य में देखी गई अंतर मुख्य रूप से अंतर्निहित नमूनों और चयन पद्धतियों के कारण है, न कि अंतर्निगत खगोलीय विविधताओं के कारण। यह अध्ययन इन तारों को उनके जन्म क्लस्टरों तक ट्रैक करने के लिए भविष्य के कार्यों हेतु एक मजबूत आधार प्रदान करता है।