🔭 astrophysics

Infall onto the Protoplanetary Disk during the Gravitational Collapse of a Molecular Cloud

यह शोध पत्र आणविक मेघ के पतन के दौरान पदार्थ के अंत:पतन (matter infall) के लिए व्यापक रूप से उपयोग किए जाने वाले नाकामोतो और नकागावा (1994) मॉडल में संशोधनों का प्रस्ताव और सत्यापन करता है, जो यह प्रदर्शित करता है कि ये संवर्द्धन युवा तारकीय पिंड (young stellar object) के विकास के अंतर्निहित चरण के लिए अनुमानित अवधि को कम करते हैं, विशेष रूप से उच्च प्रारंभिक घनत्व विक्षोभ वाले क्षेत्रों में।

मूल लेखक: E. R. Redkin, E. I. Vorobyov

प्रकाशित 2026-02-25

📖 6 मिनट में पढ़ें🧠 गहराई से पढ़ें

CC BY 4.0

मूल लेखक: E. R. Redkin, E. I. Vorobyov

मूल पेपर CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/) के तहत लाइसेंस किया गया है। ✨ नीचे दिए गए पेपर की यह व्याख्या AI से तैयार की गई है। इसे लेखकों ने न तो लिखा है, न इसका समर्थन किया है। तकनीकी सटीकता के लिए मूल पेपर देखें। पूरा डिस्क्लेमर पढ़ें

मुख्य चित्र: तारे और उनके ग्रहों के नर्सरी का जन्म कैसे होता है

कल्पना कीजिए कि अंतरिक्ष में गैस और धूल का एक विशाल, ठंडा बादल तैर रहा है। यह एक आणविक बादल (molecular cloud) है। कभी-कभी, कुछ ऐसा होता है (जैसे किसी पास के फटने वाले तारे से निकली शॉकवेव) जो इस बादल के एक छोटे से हिस्से को अपने स्वयं के गुरुत्वाकर्षण के कारण ढहने (collapse होने) के लिए मजबूर कर देता है।

जैसे-जैसे यह हिस्सा ढहता है, यह तेजी से घूमने लगता है (जैसे एक आइस स्केटर अपने हाथों को अंदर खींचकर घूमता है)। अधिकांश गैस सीधे केंद्र में गिर जाती है जिससे एक नन्हा तारा (protostar) बनता है। लेकिन उस घूमने के कारण, सारा गैस सीधे अंदर नहीं गिर सकता। इसका कुछ हिस्सा तारे के चारों ओर एक चपटे, घूमते हुए पैनकेक की तरह बाहर की ओर उछल जाता है। यह प्रोटोप्लेनेटरी डिस्क (protoplanetary disk) है, वह नर्सरी जहाँ अंततः ग्रह बनेंगे।

समस्या: पुराना नक्शा थोड़ा त्रुटिपूर्ण था

वैज्ञानिक इस प्रक्रिया को सिम्युलेट करने के लिए एक विशिष्ट "नक्शे" (नाकामोतो और नकागावा द्वारा 1994 में बनाया गया एक गणितीय मॉडल) का उपयोग कर रहे हैं कि कैसे यह गैस बादल से डिस्क पर गिरती है। यह एक लोकप्रिय नक्शा है क्योंकि इसे उपयोग करना सरल है।

हालाँकि, इस शोध पत्र के लेखक, रेडकिन और वोरोब्योव ने महसूस किया कि इस नक्शे में कुछ खामियां (glitches) हैं:

बादल अनंत था: पुराने मॉडल ने माना कि बादल अनंत रूप से बड़ा था। वास्तव में, बादलों के किनारे होते हैं।
घनत्व "परफेक्ट" था: पुराने मॉडल ने माना कि केंद्र में गैस गणितीय रूप से एकदम सटीक और एकल तरीके से भरी हुई थी। वास्तविक बादल थोड़े अधिक "ऊबड़-खाबड़" होते हैं और घने होने से पहले उनका केंद्र थोड़ा चपटा होता है।
बारिश कभी नहीं रुकी: पुराने मॉडल ने माना कि गैस तारे पर एक निरंतर दर से तब तक गिरती रहती है जब तक कि बादल खाली न हो जाए। वास्तव में, जैसे-जैसे बादल ढहता है, गैस की "बारिश" धीमी हो जानी चाहिए और अंततः रुक जानी चाहिए।

समाधान: सिमुलेशन को अपडेट करना

लेखकों ने इन खामियों को ठीक करने के लिए तीन मुख्य अपग्रेड प्रस्तावित किए, जिससे सिमुलेशन को कंप्यूटर के लिए बहुत जटिल बनाए बिना अधिक यथार्थवादी बनाया जा सके।

1. बादल को एक सीमा देना ("बाड़" का उदाहरण)

कल्पना कीजिए कि बादल एक मैदान में लोगों की भीड़ की तरह है। पुराना मॉडल मानता था कि मैदान अनंत तक जाता है। नया मॉडल मैदान के चारों ओर एक बाड़ (fence) लगा देता है।

यह क्यों महत्वपूर्ण है: जैसे ही पतन (collapse) शुरू होता है, खालीपन की एक "लहर" केंद्र से बाहर की ओर बढ़ती है। जब यह लहर बाड़ से टकराती है, तो गैस की आपूर्ति कट जाती है। यह वैज्ञानिकों को यह गणना करने की अनुमति देता है कि बादल में कितनी ऊर्जा है (घूर्णन बनाम गुरुत्वाकर्षण), जो यह समझने के लिए महत्वपूर्ण है कि बनने वाली डिस्क कितनी बड़ी होगी।

2. केंद्र को सुचारू बनाना ("मफिन" का उदाहरण)

पुराने मॉडल ने कहा कि बादल का केंद्र अनंत रूप से घना था, जैसे कि एक गणितीय बिंदु। नया मॉडल कहता है कि केंद्र एक मफिन के ऊपरी हिस्से की तरह है—यह चपटा और घना है, लेकिन अनंत नहीं है।

यह क्यों महत्वपूर्ण है: यह वही है जो हम तारा निर्माण के वास्तविक 3D कंप्यूटर सिमुलेशन में देखते हैं। यह शुरुआती चरणों के दौरान गैस के गिरने की गति को बदल देता है।

3. "धीमा करने वाला बटन" ("लीकी बकेट" का उदाहरण)

पुराना मॉडल एक ऐसी बाल्टी की तरह था जिसमें छेद था और वह खाली होने तक एक स्थिर गति से लीक हो रही थी। नया मॉडल यह समझता है कि जैसे-जैसे बाल्टी खाली होती है, पानी का दबाव कम हो जाता है, और रिसाव धीमा हो जाता है।

यह क्यों महत्वपूर्ण है: उन्होंने एक नया फॉर्मूला बनाया जहाँ तारे पर गिरने वाली गैस की दर शुरुआत में उच्च होती है, कुछ समय तक स्थिर रहती है, और फिर जैसे-जैसे बादल में सामग्री खत्म होती है, यह एक्सपोनेंशियल रूप से घटती (fade out) जाती है। यह भौतिक रूप से बहुत अधिक सटीक है।

परिणाम: तारे हमारी सोच से कहीं अधिक तेजी से बड़े होते हैं

जब लेखकों ने अपने नए, बेहतर सिमुलेशन चलाए, तो उन्हें एक नन्हे तारे के जीवन के टाइमलाइन के बारे में एक आश्चर्यजनक बात पता चली।

खगोलशास्त्री नवजात तारों को इस आधार पर वर्गीकृत करते हैं कि उन पर कितनी गैस अभी भी गिर रही है:

Class 0: "नवजात" चरण। तारा गैस में गहराई से दबा हुआ है; बादल का 50% अभी भी गिर रहा है।
Class I: "टॉडलर" (छोटे बच्चे) चरण। अधिकांश गैस गिर चुकी है, लेकिन अभी भी एक मोटा आवरण मौजूद है।
Class II: "टीनएजर" (किशोर) चरण। गैस लगभग खत्म हो गई है, और हम डिस्क को स्पष्ट रूप से देख सकते हैं (जहाँ ग्रह बनते हैं)।

निष्कर्ष:
नए, अधिक यथार्थवादी मॉडल का उपयोग करते हुए, उन्होंने पाया कि Class 0 से Class I में संक्रमण पिछले मॉडलों के सुझाव की तुलना में बहुत अधिक तेजी से होता है।

उदाहरण: कल्पना कीजिए कि एक किशोर बड़ा हो रहा है। पुराने मॉडल कहते थे, "एक बच्चे से टॉडलर बनने में 500,000 साल लगते हैं।" नए मॉडल कहते हैं, "वास्तव में, इसमें केवल 100,000 साल लगते हैं।"

अंतर क्यों है?
नया मॉडल दिखाता है कि यदि प्रारंभिक बादल थोड़ा "ऊबड़-खाबड़" (उच्च घनत्व विचलन) है, तो गैस बहुत तेज़ी से गिरती है। तारा तेज़ी से बड़ा होता है।

"इसका क्या महत्व है?"

यह इसलिए महत्वपूर्ण है क्योंकि खगोलशास्त्री आकाश में वास्तविक तारों को देखते हैं और उनके आसपास कितनी गैस है इसके आधार पर उनकी उम्र का अनुमान लगाने की कोशिश करते हैं। यदि हमारे कंप्यूटर मॉडल कहते हैं कि "यह प्रक्रिया 500,000 वर्षों में होती है," लेकिन वास्तविक तारे इसे 100,000 वर्षों में कर रहे हैं, तो हम उनकी उम्र का गलत अनुमान लगा सकते हैं।

लेखक सुझाव देते हैं कि शायद वास्तविक तारे हमारी सोच से कहीं अधिक तेज़ी से बड़े हो रहे हैं, या शायद जिस तरह से हम उन्हें वर्गीकृत करते हैं वह पेचीदा है (जैसे कि केवल ऊंचाई देखकर यह बताना कि कोई किशोर 13 साल का है या 15 साल का)।

सारांश

यह शोध पत्र तारा जन्म के सिमुलेशन के लिए एक पुराने, लोकप्रिय उपकरण के छेद को भरने (patching the holes) का काम करता है। उन्होंने एक सीमा जोड़ी, केंद्र को सुचारू बनाया, और गैस के प्रवाह को स्वाभाविक रूप से धीमा किया। परिणाम? तारों और उनके ग्रहीय डिस्क के जन्म की एक अधिक यथार्थवादी तस्वीर, जो बताती है कि "एम्बेडेड" चरण (जहाँ तारा गैस में छिपा होता है) हमारी पिछली धारणा से छोटा है।

यहाँ रेडकिन और वोरोब्योव के शोध पत्र "Infall onto the Protoplanetary Disk during the Gravitational Collapse of a Molecular Cloud" का विस्तृत तकनीकी सारांश दिया गया है।

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

प्रोटोप्लेनेटरी डिस्क का निर्माण इस बात पर निर्भर करता है कि एक ढहते हुए आणविक बादल (molecular cloud) से केंद्रीय प्रोटोस्टार और उसके आसपास की डिस्क पर पदार्थ कैसे गिरता है, इसके सटीक मॉडल क्या हैं। व्यापक रूप से उपयोग किया जाने वाला नाकामोतो और नकागावा (NN) मॉडल (कासेन और मूस के कार्य पर आधारित) इस इनफॉल दर (infall rate) की गणना करने के लिए एक सरलीकृत दृष्टिकोण प्रदान करता है। हालाँकि, लेखक मूल NN मॉडल में तीन महत्वपूर्ण सीमाओं की पहचान करते हैं जो संख्यात्मक सिमुलेशन (numerical simulations) की भौतिक वास्तविकता को प्रभावित करती हैं:

आदर्शित घनत्व वितरण (Idealized Density Distribution): मॉडल एक सिंगुलर आइसोथर्मल स्फीयर (SIS) घनत्व प्रोफाइल ( $\rho \propto r^{-2}$ ) मानता है, जो केंद्र में अनंत घनत्व का संकेत देता है और इसमें एक भौतिक बाहरी सीमा का अभाव है।
असीमित क्लाउड आकार (Unlimited Cloud Size): मूल मॉडल में एक परिमित क्लाउड त्रिज्या (finite cloud radius) को परिभाषित नहीं किया गया है, जिससे महत्वपूर्ण भौतिक मापदंडों जैसे कि घूर्णन और गुरुत्वाकर्षण ऊर्जा के अनुपात ( $\beta$ ) की गणना करना या सीमित द्रव्यमान भंडार (mass reservoir) को ध्यान में रखना असंभव हो जाता है।
स्थिर इनफॉल दर (Constant Infall Rate): मॉडल एक स्थिर द्रव्यमान इनफॉल दर ( $\dot{M}_c$ ) मानता है जब तक कि क्लाउड समाप्त न हो जाए, जो इस भौतिक वास्तविकता की अनदेखी करता है कि क्लाउड की बाहरी सीमा से एक रिफ्रेक्शन वेव (rarefaction wave) केंद्र तक पहुँचने के बाद इनफॉल दर कम होनी चाहिए।

ये सीमाएँ युवा तारकीय पिंडों (YSOs) के "एम्बेडेड चरण" (Class 0 और Class I अवस्थाओं) की अवधि के संबंध में सिमुलेशन परिणामों और अवलोकन संबंधी डेटा के बीच विसंगतियों का कारण बनती हैं।

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

लेखक मूल NN मॉडल में किए गए सुधारों का एक प्रस्ताव देते हैं ताकि एक-आयामी श्यान विकास कोड (one-dimensional viscous evolution codes) के लिए इसकी कम्प्यूटेशनल सरलता को बनाए रखते हुए पहचानी गई कमियों को दूर किया जा सके।

A. एक परिमित क्लाउड सीमा का परिचय (Introduction of a Finite Cloud Boundary)

लेखक कुल क्लाउड द्रव्यमान ( $M_c$ ) के आधार पर आणविक बादल की बाहरी त्रिज्या ( $r_{out}$ ) की गणना का परिचय देते हैं। यह घूर्णन ऊर्जा और गुरुत्वाकर्षण ऊर्जा के अनुपात ( $\beta$ ) की गणना करने की अनुमति देता है, जो डिस्क निर्माण के लिए एक महत्वपूर्ण पैरामीटर है।

परिणाम: एक मानक सौर-द्रव्यमान वाले क्लाउड के लिए, यह $r_{out} \approx 8337$ au और $\beta \approx 0.19\%$ प्राप्त करता है, जो प्रेक्षित प्रीस्टेलर कोर (prestellar cores) के अनुरूप है।

B. घनत्व वितरण में संशोधन (Modification of Density Distribution)

सिंगुलर आइसोथर्मल स्फीयर (SIS) प्रोफाइल को एक बोनर-एबर्ट स्फीयर-जैसे प्रोफाइल (विशेष रूप से लार्सन-पेन्स्टन प्रोफाइल) से बदल दिया गया है, जो एक केंद्रीय पठार (central plateau) द्वारा अभिलक्षित है:
$\rho(r) = \frac{\rho_0}{1 + (r/r_0)^2}$

औचित्य: यह प्रोफाइल ढहते हुए बादलों के केंद्रीय क्षेत्रों में भौतिक स्थितियों का बेहतर प्रतिनिधित्व करता है (जो $r=0$ पर विलक्षणता/singularity से बचता है) और SIS प्रोफाइल की तुलना में 3D हाइड्रोडायनामिक सिमुलेशन डेटा के अधिक करीब है।
प्रभाव: यह परिवर्तन समय और सेंट्रीफ्यूगल रेडियस ( $R_{cf}$ ) के बीच के संबंध को बदल देता है, जिसके लिए $R_{cf}(t)$ के नए व्युत्पत्ति की आवश्यकता होती है।

C. समय-निर्भर इनफॉल दर (Time-Dependent Infall Rate)

स्थिर इनफॉल दर के अनुमान को क्लाउड के परिमित द्रव्यमान को ध्यान में रखते हुए एक समय-निर्भर फलन (time-dependent function) से बदल दिया गया है।

तंत्र: इनफॉल दर स्थिर ( $\dot{M}_c = m_0 a^3/G$ ) रहती है जब तक कि $t_d = r_{out}/2a$ समय पर बाहरी सीमा से यात्रा करने वाली रिफ्रेक्शन वेव केंद्र तक नहीं पहुँच जाती। $t_d$ के बाद, दर तेजी से (exponentially) घटती है।
परिशोधन: लेखक एक विशिष्ट क्षय फलन (decay function) (संशोधन M2) का प्रस्ताव करते हैं जो यह सुनिश्चित करता है कि समय के साथ एकीकृत कुल द्रव्यमान प्रारंभिक क्लाउड द्रव्यमान ( $M_c$ ) के बराबर हो, जिससे क्लाउड में "लॉक" (locked) द्रव्यमान की समस्या को रोका जा सके।
गुणांक समायोजन: शोध पत्र घनत्व विचलन आयाम (density perturbation amplitude) $A$ के फलन के रूप में गुणांक $m_0$ के लिए एक बहुपद सन्निकटन (polynomial approximation) प्रदान करता है, जो $2 < A < 5$ के लिए पिछले साहित्य में पाई गई अति-अनुमानों को ठीक करता है।

D. डिस्क-मध्यस्थ अभिवृद्धि (Disk-Mediated Accretion)

लेखक इनफॉल समीकरण को परिष्कृत करते हैं ताकि सीधे प्रोटोस्टार पर गिरने वाले पदार्थ और डिस्क पर गिरने वाले पदार्थ के बीच अंतर किया जा सके। एक सुधार कारक $\xi(t)$ पेश किया गया है ताकि यह सुनिश्चित किया जा सके कि द्रव्यमान एकीकरण (mass integration) डिस्क की आंतरिक सीमा ( $R_{in}$ ) से शुरू होता है, न कि शून्य से, एक बार जब सेंट्रीफ्यूगल रेडियस $R_{in}$ से अधिक हो जाता है।

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

भौतिक यथार्थवाद (Physical Realism): संशोधित मॉडल (NN-1 से NN-7) परिमित क्लाउड सीमाओं और यथार्थवादी घनत्व प्रोफाइल को शामिल करता है, जिससे $\beta$ की गणना और अधिक सटीक प्रारंभिक स्थितियाँ संभव होती हैं।
विश्लेषणात्मक क्षय नियम (Analytical Decay Law): क्षय होने वाली इनफॉल दर के लिए एक विश्लेषणात्मक अभिव्यक्ति का व्युत्पत्ति जो सख्ती से कुल क्लाउड द्रव्यमान का संरक्षण करती है, जो पिछले मॉडलों में "अनंत" या "लॉक" द्रव्यमान की समस्या को हल करती है।
सामान्यीकृत इनफॉल समीकरण (Generalized Infall Equation): सतह घनत्व विकास के लिए एक सामान्यीकृत समीकरण का प्रावधान जो सब-केपलरियन इनफॉल वेग और डिस्क के विशिष्ट आंतरिक सीमा को ध्यान में रखता है।
पैरामीटर संवेदनशीलता विश्लेषण (Parameter Sensitivity Analysis): इस बात का व्यवस्थित अध्ययन कि घनत्व विचलन आयाम ( $A$ ) और क्लाउड तापमान ( $T_c$ ) एम्बेडेड चरण की अवधि को कैसे प्रभावित करते हैं।

4. परिणाम (Results)

लेखकों ने विभिन्न मापदंडों ( $A$ , $T_c$ , $\omega_c$ ) के साथ संशोधित मॉडलों (NN-1 से NN-7) का उपयोग करके सिमुलेशन चलाए और परिणामों की तुलना Class 0 और Class I YSOs की अवधि के अवलोकन संबंधी अनुमानों से की।

एम्बेडेड चरण की अवधि: सिमुलेशन आम तौर पर अवलोकन संबंधी अनुमानों (0.15–0.48 Myr) की तुलना में एम्बेडेड चरण की छोटी अवधि की भविष्यवाणी करते हैं।
- उच्च घनत्व विचलन आयाम ( $A > 3$ ) और मानक तापमान ( $T_c = 10-15$ K) वाले मॉडल बहुत तीव्र पतन दिखाते हैं, जो 0.2 Myr से भी कम समय में Class 0 से Class II में परिवर्तित हो जाते हैं।
- केवल कम तापमान ( $T_c \approx 7$ K) या बहुत कम विचलन आयाम ( $A \approx 2.2$ ) वाले मॉडल ही अवलोकनों के तुलनीय संक्रमण समय उत्पन्न करते हैं (जैसे, मॉडल NN-7: $t_{C0} \approx 231$ kyr, $t_{CI} \approx 706$ kyr)।
द्रव्यमान वितरण: सभी मॉडलों में, अधिकांश द्रव्यमान (लगभग 80%) सीधे तारे के बजाय डिस्क पर गिरता है, जो "डिस्क-मध्यस्थ अभिवृद्धि" परिदृश्य की पुष्टि करता है।
विसंगति विश्लेषण: लेखक सुझाव देते हैं कि सिमुलेशन और अवलोकन के बीच बेमेल होने का कारण निम्नलिखित हो सकता है:
- अवलोकन संबंधी पूर्वाग्रह (Observational Bias): धूल के अवरोध के कारण एज-ऑन (edge-on) सिस्टम Class 0 के रूप में दिखाई दे सकते हैं, जिससे प्रारंभिक चरणों की देखी गई अवधि कृत्रिम रूप से बढ़ सकती है।
- मॉडल की सीमाएं: एक अलग (isolated) क्लाउड का अनुमान आस-पास के फिलामेंट्स से द्रव्यमान अभिवृद्धि को अनदेखा करता है, जो वास्तव में इनफॉल चरण को लंबा कर सकता है।

5. महत्व (Significance)

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

निष्कर्षों से पता चलता है कि प्रारंभिक घनत्व विचलन आयाम तारा निर्माण की समयसीमा निर्धारित करने में एक महत्वपूर्ण कारक है। यह शोध पत्र भविष्य के अनुसंधान के लिए आधार तैयार करता है जो इन उन्नत इनफॉल मॉडलों का उपयोग करके प्रोटोप्लेनेटरी डिस्क के श्यान विकास और ग्रह निर्माण का अनुकरण करेगा, जिससे आणविक क्लाउड पतन और ग्रहीय प्रणाली संरचना के बीच एक अधिक सटीक सेतु प्रदान होगा।