Kick Velocities and Mass Function of Free-Floating Planets from Dynamical Ejection in Hierarchical Three-Body Systems

तकनीकी सारांश: पदानुक्रमित त्रि-पिंड प्रणालियों (Hierarchical Three-Body Systems) में गतिक निष्कासन (Dynamical Ejection) से मुक्त-तैरते ग्रहों की किक वेग (Kick Velocities) और द्रव्यमान फलन (Mass Function)

समस्या विवरण
मुक्त-तैरते ग्रह (FFPs), या आवारा ग्रह (rogue planets), वे उप-तारा पिंड हैं जो किसी मेजबान तारे से बंधे नहीं होते। उनका वेग और द्रव्यमान वितरण उनके निर्माण और निष्कासन इतिहास के गतिज जीवाश्म (kinematic fossils) के रूप में कार्य करते हैं। जबकि गुरुत्वाकर्षण माइक्रोलेंसिंग ने पृथ्वी-द्रव्यमान से लेकर बृहस्पति-द्रव्यमान तक के FFPs की एक आबादी का पता लगाया है, उनके निष्कासन वेगों और परिणामी द्रव्यमान फलन को नियंत्रित करने वाले गतिक तंत्र अभी भी अपूर्ण रूप से लक्षणित हैं। विशेष रूप से, पदानुक्रमित त्रि-पिंड प्रणालियों (एक मेजबान तारा, एक विशाल विचलित करने वाला पिंड, और एक हल्का आंतरिक ग्रह) के निष्कासन परिणामों को FFPs के अवलोकन योग्य गुणों, जैसे कि उनके अनंत वेग ($v_\infty$) और निष्कासन के समय ($t_{ejec}$) के साथ जोड़ने की आवश्यकता है।

कार्यप्रणाली
लेखक उच्च सटीकता सुनिश्चित करने के लिए (विशेष रूप से निकट संपर्क के दौरान) IAS15 इंटीग्रेटर (15वीं-क्रम का अनुकूली टाइमस्टेपिंग) के साथ rebound इंटीग्रेशन पैकेज का उपयोग करके प्रत्यक्ष $N$-बॉडी सिमुलेशन का एक समूह प्रस्तुत करते हैं। यह अध्ययन प्रतिबंधित त्रि-पिंड समस्या (restricted three-body problem) पर केंद्रित है, जिसमें एक सौर-द्रव्यमान वाले मेजबान तारे ($M_\star = 1 M_\odot$) और एक विशाल विचलित करने वाले पिंड के सापेक्ष आंतरिक ग्रह को एक परीक्षण कण (test particle) के रूप में माना गया है।

सिमुलेशन ग्रिड नियंत्रित और व्यवस्थित है:

• निश्चित पैरामीटर: मेजबान द्रव्यमान ($1 M_\odot$), विशाल ग्रह का अर्ध-दीर्घ अक्ष (semi-major axis, $a_J = 5.2$ AU, एक बृहस्पति का प्रतिरूप), और प्रारंभिक झुकाव ($0^\circ$–$3^\circ$)।
• परिवर्तित पैरामीटर:
  • निष्कासित ग्रह द्रव्यमान ($m_p$): $10^{-5}$ से $10^{-1} M_{Jupiter}$।
  • विशाल ग्रह द्रव्यमान ($M_J$): $1$से$75 M_{Jupiter}$।
  • निष्कासित ग्रह का अर्ध-दीर्घ अक्ष ($a_p$): $0.7$से$1.5 a_J$।
  • दोनों ग्रहों की उत्केंद्रता (eccentricities) ($e_p$ और $e_J$): $0.0$से$0.9$।
• साक्षात्कार (Realizations): प्रत्येक विन्यास को पैरामीटर प्रभावों को अलग करने के लिए साझा रैंडम सीड के साथ $N=10,000$ स्वतंत्र साक्षात्कारों के साथ सिम्युलेट किया गया है।
• मेट्रिक्स: अध्ययन अनंत निष्कासन वेग ($v_\infty$), निष्कासन पुष्टि के क्षण पर परिमित-त्रिज्या किक वेग ($v_{kick}$), और निष्कासन समय ($t_{ejec}$) को मापता है। निष्कासन को विशिष्ट ऊर्जा धनात्मकता, प्रणाली के अधिकतम अर्ध-दीर्घ अक्ष से दूरी, और दस कक्षीय अवधियों तक बाहर की ओर रेडियल वेग के बने रहने द्वारा परिभाषित किया गया है।

विश्लेषणात्मक ढांचा
संख्यात्मक परिणामों की व्याख्या निम्नलिखित के संयोजन वाले एक अर्ध-विश्लेषणात्मक ढांचे के माध्यम से की जाती है:

1. ऊर्जा बजट: निष्कासन के लिए आवश्यक ऊर्जा वृद्धि को बाध्य कक्षीय ऊर्जा और अनंत गतिज ऊर्जा के बीच के अंतर से प्राप्त किया जाता है।
2. गुरुत्वाकर्षण स्लिंगशॉट (Gravitational Slingshot): अधिकतम किक वेग विशाल ग्रह के विश्राम फ्रेम में स्लिंगशॉट तंत्र द्वारा सीमित है, जो मुठभेड़ की गति $u$ और विशाल ग्रह के कक्षीय वेग पर निर्भर करता है।
3. हिल-स्केल स्कैटरिंग (Hill-Scale Scattering): औसत निष्कासन गति विशाल ग्रह के हिल क्षेत्र (Hill sphere) के भीतर विशिष्ट मुठभेड़ों द्वारा नियंत्रित होती है, जो $v_H \propto M_J^{1/3}$ के रूप में स्केल करती है।
4. टिसरैंड पैरामीटर (Tisserand Parameter): कक्षीय तत्वों को मुठभेड़ की गति से जोड़ने के लिए एक संरक्षित नैदानिक (diagnostic) के रूप में उपयोग किया जाता है।
5. सेकुलर विकास (Secular Evolution): अस्थिरता के लिए समय (timescales) को सेकुलर पंपिंग और मीन-मोशन रेजोनेंस (MMR) ओवरलैप से जोड़ा जाता है।

मुख्य परिणाम

1. टेस्ट-पार्टिकल सीमा: चार दशकों की सीमा ($10^{-5}$ से $10^{-1} M_{Jupiter}$) में निष्कासन सांख्यिकी पर निष्कासित ग्रह के द्रव्यमान ($m_p$) का नगण्य प्रभाव पड़ता है। निष्कासन अंश 2.3% से कम भिन्न होता है, और औसत वेग स्थिर रहता है। यह पुष्टि करता है कि $m_p/M_J \lesssim 10^{-2}$ के लिए, गतिकी विशाल विचलित करने वाले पिंड द्वारा हावी होती है, और निष्कासित पिंड एक परीक्षण कण की तरह व्यवहार करता है।
2. विशाल ग्रह द्रव्यमान स्केलिंग: विशाल ग्रह का द्रव्यमान ($M_J$) निष्कासन पैमाने और समय को निर्धारित करता है।
   • औसत निष्कासन वेग $\langle v_\infty \rangle \propto M_J^{0.42}$ के रूप में स्केल करता है।
   • निष्कासन समय $\langle t_{ejec} \rangle \propto M_J^{-1.09}$ के रूप में स्केल करता है, जो सेकुलर भविष्यवाणी $t_{sec} \propto M_J^{-1}$ के अनुरूप है।
   • निष्कासन अंश गैर-मोनोटोनिक (non-monotonic) है, जो $10 M_{Jupiter}$ के पास चरम पर पहुँचता है और बहुत कम तथा बहुत उच्च द्रव्यमान पर घट जाता है (सेकुलर समय सीमाओं और हिल क्षेत्र कैप्चर के कारण)।
3. अर्ध-दीर्घ अक्ष निर्भरता: रेजोनेंस-ओवरलैप ज़ोन ($|a_p - a_J|/a_J \lesssim 0.40$ के लिए $10 M_{Jupiter}$ विचलित करने वाले पिंड हेतु) के भीतर विन्यासों के लिए, निष्कासन समय फ्लैट ($\sim 1.9 \times 10^4$ वर्ष) रहता है, जो सेकुलर विकास के बजाय तेज़ अराजक ओवरलैप समय द्वारा शासित होता है। औसत वेग $a_p = 0.9 a_J$ के पास चरम पर होता है।
4. उत्केंद्रता प्रभाव (दो-गति प्रतिक्रिया):
   • आंतरिक ग्रह की उत्केंद्रता ($e_p$): मुख्य रूप से उच्च-वेग वाली पूंछ (high-velocity tail) को प्रभावित करती है। $e_p$ को 0 से 0.9 तक बढ़ाने से अधिकतम किक वेग लगभग 3 गुना ($\sim 7.8$ से $\sim 23$ किमी/सेकेंड) बढ़ जाता है, लेकिन केवल औसत वेग को लगभग 1.7 गुना ही बढ़ाता है।
   • विशाल ग्रह की उत्केंद्रता ($e_J$): पूंछ पर काफी अधिक प्रभाव डालती है। $e_J$ को 0 से 0.9 तक बढ़ाने से अधिकतम किक वेग लगभग 10 गुना ($\sim 7.8$ से $\sim 80$ किमी/सेकेंड) बढ़ जाता है। यह विशाल ग्रह के पेरीसेंटर वेग ($v_{peri} \propto \sqrt{(1+e_J)/(1-e_J)}$) के गैर-रेखीय विकास द्वारा संचालित है, जो दुर्लभ, उच्च-ऊर्जा स्लिंगशॉट मुठभेड़ों की अनुमति देता है।
5. विश्लेषणात्मक ऊपरी लिफाफा (Analytic Upper Envelope): लेखक स्लिंगशॉट तंत्र और टिसरैंड पैरामीटर के आधार पर अधिकतम निष्कासन गति के लिए एक विश्लेषणात्मक सूत्र व्युत्पन्न करते हैं। यह लिफाफा सिम्युलेटेड वेग वितरण की ऊपरी सीमा का सटीक अनुमान लगाता है, विशेष रूप से $e_p \gtrsim 0.5$ के लिए।

महत्व और निहितार्थ
लेखक दावा करते हैं कि ये परिणाम ग्रहों की व्यवस्था की वास्तुकला और गैलेक्टिक FFP आबादी के अवलोकन योग्य गुणों के बीच एक सीधा संबंध प्रदान करते हैं:

• द्रव्यमान फलन (Mass Function): चूंकि निष्कासन अंश ग्रह के द्रव्यमान से स्वतंत्र है, इसलिए वर्तमान FFP द्रव्यमान फलन सीधे अंतर्निहित ग्रह उपस्थिति दर ($dN_{FFP}/dm_p \approx \bar{f}_{eject} \cdot dN_{planet}/dm_p$) को दर्शाता है। इसका अर्थ है कि यदि निम्न द्रव्यमान की ओर ग्रह उपस्थिति दर बढ़ती है (जैसा कि केप्लर डेटा द्वारा सुझाया गया है), तो उप-पृथ्वी (sub-Earth) FFP आबादी संख्यात्मक रूप से प्रमुख होगी।
• माइक्रोलेंसिंग संकेत: रूपांतरित माइक्रोलेंसिंग समय-पैमाने का वितरण ($dN/d\log t_E$) $\log t_E$ में फ्लैट होगा यदि ग्रह उपस्थिति दर फ्लैट है, या छोटे समय-पैमानों पर प्रबल होगा यदि उपस्थिति दर निम्न द्रव्यमान की ओर बढ़ती है। यह रोमन वाइड फील्ड सर्वे और यूक्लिड (Euclid) के लिए एक परीक्षण योग्य भविष्यवाणी प्रदान करता है।
• गतिकी (Kinematics): अधिकांश निष्कासित ग्रहों का अनंत वेग केवल कुछ किमी/सेकेंड होता है, जिससे वे अपने मूल तारकीय जनसंख्या से गतिक रूप से अविभेद्य हो जाते हैं (3D वेग फैलाव को $\sim 1\%$ बदलते हैं)। हालांकि, अत्यधिक उत्केंद्रित विशाल विचलित करने वाले पिंडों ($e_J \gtrsim 0.5$) वाले सिस्टम एक उच्च-वेग वाली पूंछ ($v_\infty \sim 80$ किमी/सेकेंड) उत्पन्न कर सकते हैं, जो संभवतः गाया (Gaia) या रोमन द्वारा पता लगाने योग्य प्रॉपर-मोशन आउटलेर्स बना सकते हैं।
• निष्कासन समय-सीमा: विशाल या उत्केंद्रित प्रणालियों से अधिकांश निष्कासन तेजी से ($10^3$–$10^4$ वर्ष) होते हैं, अक्सर प्री-मेन-सीक्वेंस चरण के दौरान। फलस्वरूप, सामान्य FFP आबादी का स्थानिक वितरण फील्ड (field) से प्रभावी रूप से अविभेद्य है, केवल युवा, हाल ही में निष्कासित पिंड ही अपने जन्म के वातावरण की स्मृति बनाए रखते हैं।

सीमाएं
लेखक नोट करते हैं कि केवल तीन पिंडों तक सीमित रहना प्राथमिक सीमा है, क्योंकि वास्तविक प्रणालियों में कई ग्रह शामिल हो सकते हैं जो कैस्केडिंग अस्थिरता (cascading instabilities) का कारण बनते हैं। इसके अतिरिक्त, सिमुलेशन में प्रोटोप्लेनेटरी डिस्क, तारकीय द्रव्यमान हानि, या पूर्ण रेंज के झुकाव (विशेष रूप से कोज़ाई-लिडव (Kozai-Lidov) तंत्र) को शामिल नहीं किया गया है, जो निष्कासन सांख्यिकी को बदल सकते हैं। भविष्य के कार्य को इन कारकों को पूर्ण जनसंख्या संश्लेषण मॉडल में शामिल करने के लिए आवश्यक माना गया है।