Building three-dimensional giant stellar models for common envelope simulations

तकनीकी सारांश: कॉमन एनवेलप सिम्युलेशन के लिए त्रि-आयामी विशाल तारकीय मॉडल का निर्माण

समस्या विवरण
कॉमन एनवेलप इवोल्यूशन (CEE) सिमुलेशन महत्वपूर्ण चुनौतियों का सामना करते हैं क्योंकि एक विकसित प्राथमिक तारे (RGB, AGB, या RSG) के फूले हुए आवरण (envelope) और सघन द्वितीयक साथी (secondary companion) दोनों को हल करने के लिए एक विशाल गतिशील रेंज (dynamical range) की आवश्यकता होती है। मौजूदा संख्यात्मक कोड अक्सर बड़े आवरण, संवहन (convection), अभिवृद्धि प्रक्रियाओं (accretion processes) और संभावित जेट लॉन्चिंग को एक साथ हल करने में संघर्ष करते हैं। वर्तमान पद्धतियों में एक महत्वपूर्ण सीमा प्रारंभिक तारकीय मॉडल का उपचार है। कई अध्ययन एक 1D तारकीय मॉडल को 3D ग्रिड पर ले जाते हैं और फिर मॉडल को "रिलैक्स" (वेगों को कम करना) करते हैं ताकि सिमुलेशन को शून्य आवरण वेग के साथ शुरू किया जा सके। यह दृष्टिकोण विशाल तारों में निहित प्राकृतिक स्पंदनों (pulsations) और संवहन को दबा देता है, जो द्रव्यमान हानि, अभिवृद्धि दर और CEE की शुरुआत को प्रभावित करने के लिए जाने जाते हैं। इसके अलावा, बिना ऊर्जा इनपुट के अनपरर्ब्तित (unperturbed) 3D मॉडलों को सिम्युलेट करने के पिछले प्रयासों में केवल दो अवधियों के बाद दोलनों का तेजी से क्षय हुआ, जिससे वास्तविक विशाल तारे की गतिशील अवस्था को बनाए रखने में विफलता मिली।

कार्यप्रणाली
लेखकों ने CEE सिमुलेशन के लिए 3D रेड सुपरजाइंट (RSG) मॉडल के निर्माण के लिए एक नया दृष्टिकोण विकसित किया है, जिसमें FLASH हाइड्रोडायनामिकल कोड का उपयोग किया गया है। कार्यप्रणाली में एक 1D RSG मॉडल (जो MESA से प्राप्त है, जिसका प्रारंभिक द्रव्यमान $15 M_\odot$ है और जो $881 R_\odot$ की त्रिज्या के साथ $12.5 M_\odot$ तक विकसित हुआ है) को 3D कार्टेशियन ग्रिड पर स्थानांतरित करना शामिल है। मुख्य संख्यात्मक घटक निम्नलिखित हैं:

• कोई रिलैक्सेशन नहीं: पिछले अध्ययनों के विपरीत, मॉडल को रिलैक्स नहीं किया जाता है; इसे अपने प्राकृतिक वेगों के साथ विकसित होने दिया जाता है, जिससे स्पंदन (pulsation) संभव होता है।
• इनर्ट कोर (Inert Core): गणना का समय बचाने के लिए, एक आंतरिक जड़ (inert) गोला ($R_{inert} = 0.2 R_{RSG}$) स्थिर रखा जाता है, जबकि आवरण को गतिशील रूप से विकसित किया जाता है।
• फोटोस्फेरिक कूलिंग की नकल करना: विकिरण ऊर्जा हानि को सिम्युलेट करने के लिए, एक निर्धारित "संख्यात्मक फोटोस्फेरिक घनत्व" ($\rho_{ph} = 2 \times 10^{-9}$ g cm$^{-3}$) से कम घनत्व वाले संख्यात्मक सेल को "संख्यात्मक क्रोमोस्फेरिक तापमान" ($T_{ch} = 1000$ K) तक ठंडा किया जाता है।
• परमाणु ऊर्जा की नकल करना: कोर परमाणु शक्ति को सिम्युलेट करने के लिए, ऊर्जा को इनर्ट कोर के ठीक ऊपर स्थित एक आंतरिक शेल ($201 R_\odot < r < 230 R_\odot$) में 1D तारकीय चमक ($L \approx 7 \times 10^4 L_\odot$) के बराबर दर से जमा किया जाता है।
• सीमा स्थितियाँ (Boundary Conditions): ग्रिड सीमा प्रभावों को रोकने के लिए $1.63 R_{*,0}$ पर एक गोलाकार संख्यात्मक सीमा लगाई गई है, जिसके परे एक धीमी रेडियल आउटफ्लो वेलोसिटी सेट की गई है ताकि पवन (wind) की नकल की जा सके (हालांकि विस्तृत विंड मॉडलिंग भविष्य के कार्य के लिए आरक्षित है)।

अध्ययन दो प्राथमिक परिदृश्यों की तुलना करता है: एक जिसमें केवल फोटोस्फेरिक कूलिंग है और दूसरा जिसमें फोटोस्फेरिक कूलिंग और ऊर्जा जमाव दोनों शामिल हैं।

मुख्य परिणाम

1. निरंतर स्पंदन (Sustained Pulsation): ऊर्जा जमाव के बिना किए गए सिमुलेशन में, दोलनों का क्षय लगभग दो अवधियों (लगभग 3 वर्ष) के बाद हो गया। हालांकि, जब फोटोस्फेरिक कूलिंग और ऊर्जा जमाव दोनों को शामिल किया गया, तो दोलन क्षय नहीं हुए; इसके बजाय, उनका आयाम समय के साथ धीरे-धीरे बढ़ता गया।
2. स्पंदन अवधि: मुख्य स्पंदन अवधि लगभग 1 वर्ष पाई गई, जो तारकीय डायनेमिकल समय ($\tau_d \approx 0.76$ वर्ष) के तुल्य है, जो एक मौलिक रेडियल पल्स मोड (fundamental radial pulsation mode) का सुझाव देती है।
3. नॉन-रेडियल मोड: तारकीय मॉडल ने गैर-गोलाकार संरचनाएं और औसत तारकीय त्रिज्या में तीव्र, निम्न-आयाम वाले अस्थायी परिवर्तन प्रदर्शित किए। गतिज ऊर्जा ग्राफ के विश्लेषण ने शिखर (peaks) और गर्त (troughs) के विभाजन को प्रकट किया, जिसे लेखक मौलिक रेडियल मोड पर सुपरइम्पोज़्ड नॉन-रेडियल दोलन मोड की उपस्थिति के रूप में देखते हैं।
4. संवहन (Convection): मॉडलों ने जोरदार संवहन को सफलतापूर्वक पुनरुत्पादित किया, जो RSG तारों की विशेषता है। ट्रेसर मैप्स ने संकेत दिया कि संवहन ने लगभग 3 वर्षों (लगभग चार डायनेमिकल समय) के समय पैमाने पर आवरण गैस को मिश्रित किया, जिससे प्रारंभ में अलग-थलग शेल अविभेद्य हो गए।
5. वेग संरचना: वेलोसिटी मैप्स ने $30-50$ km s$^{-1}$ का परिमाण दिखाया, जो तारकीय सतह पर केप्लरियन वेगों के अनुरूप है, जो जोरदार संवहन के विकास की पुष्टि करता है।

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