Dynamo Simulations Confirm Predominantly Toroidal Fields in Near-Core Region of an Intermediate-Mass Star

तकनीकी सारांश: डायनेमो सिमुलेशन ने एक मध्यम-द्रव्यमान के तारे के निकट-कोर क्षेत्र में प्रमुख रूप से टोरोइडल क्षेत्रों की पुष्टि की

समस्या विवरण
हालिया एस्टरोसेइस्मिक विश्लेषण (Takata et al. 2025) ने मुख्य-अनुक्रम F-प्रकार के तारे KIC 9244992 पर आंतरिक चुंबकीय क्षेत्र की ज्यामिति के पहले प्रतिबंध प्रदान किए हैं, जिससे गहरे रेडिएटिव आंतरिक भाग में, विशेष रूप से ब्रंट-वैसाला आवृत्ति (Brunt–Väisälä frequency) के शिखर द्वारा परिभाषित निकट-कोर क्षेत्र में, एक मजबूत, प्रमुख रूप से टोरोइडल चुंबकीय क्षेत्र (अजीमुथल घटक $\sim 92$ kG) का संकेत मिलता है। यह निष्कर्ष एक सैद्धांतिक चुनौती प्रस्तुत करता है: वही तारा निकट-कोर और निकट-सतह क्षेत्रों के बीच उल्लेखनीय रूप से सपाट विभेदक घूर्णन (differential rotation) प्रदर्शित करता है (Saio et al. 2015), जिसमें केवल $\sim 4\%$ का अंतर है। मानक डायनेमो सिद्धांत बताता है कि मजबूत टोरोइडल क्षेत्रों के लिए महत्वपूर्ण शियर (shear) की आवश्यकता होती है, जिससे यह प्रश्न उठता है कि कमजोर वैश्विक विभेदक घूर्णन की उपस्थिति में इतने मजबूत क्षेत्रों को कैसे बनाए रखा जा सकता है। पिछले संख्यात्मक कार्यों (जैसे, Ratataingam et al. 2024) ने सुझाव दिया था कि स्थानीय शियर परतें टोरोइडल क्षेत्रों को प्रवर्धित कर सकती हैं, लेकिन KIC 9244992 के विशिष्ट द्रव्यमान और एस्टरोसेइस्मिक बाधाओं से मेल करने वाला एक सीधा सिमुलेशन उपलब्ध नहीं था।

कार्यप्रणाली
लेखकों ने एक घूमते हुए $2 M_\odot$ मुख्य-अनुक्रम तारे (जिसका केंद्रीय हाइड्रोजन द्रव्यमान अंश $X_c = 0.5$ है) को मॉडल करने के लिए SPIN कोड (Edelmann et al. 2019) का उपयोग करके त्रि-आयामी, एनैलिक (anelastic) मैग्नेटो-हाइड्रोडायनामिक (MHD) सिमुलेशन किए। सिमुलेशन 1500 रेडियल ग्रिड बिंदुओं और गोलाकार हार्मोनिक $\ell_{max} = 84$ के रिज़ॉल्यूशन के साथ एनैलिक सन्निकटन (approximation) में MHD समीकरणों को हल करते हैं।

मुख्य सेटअप पैरामीटर शामिल हैं:

• प्रारंभिक अवस्था: MESA द्वारा निर्मित एक 1D स्टेलर मॉडल, जिसमें एक "बूस्टेड" संवेक्टिव वेग (लगभग $10\times$ मिक्सिंग-लेंथ थ्योरी) का उपयोग किया गया है ताकि KIC 9244992 के तुलनीय रॉस्बी संख्या ($Ro \sim 0.1-0.3$) बनाए रखते हुए गणनात्मक रूप से व्यवहार्य रेनॉल्ड्स संख्या प्राप्त की जा सके।
• घूर्णन: तारे को सॉलिड-बॉडी रोटेशन ($\Omega \approx 1.0 \times 10^{-5}$ rad s$^{-1}$) के साथ इनिशियलाइज़ किया गया था, जो बूस्टेड संवहन की भरपाई के लिए KIC 9244992 की तुलना में लगभग 10 गुना तेज़ है।
• चुंबकीय क्षेत्र: एक डिपोल सीड फील्ड ($100$ kG बॉटम बाउंड्री पर) को एक हाइड्रोडायनामिक स्थिर अवस्था प्राप्त होने के बाद लागू किया गया था।
• पैरामीटर्स: प्राथमिक सिमुलेशन ने रेनॉल्ड्स संख्या ($Re$) 240, चुंबकीय रेनॉल्ड्स संख्या ($Rm$) 960, और चुंबकीय प्रैंड्टल संख्या ($Pm$)$\sim 4$ का उपयोग किया।
• विश्लेषण: डेटा का विश्लेषण भौतिक समय के 40 से 100 दिनों तक किया गया, जब चुंबकीय ऊर्जा एक सांख्यिकीय रूप से स्थिर अवस्था में पहुँच गई थी। अध्ययन मुख्य रूप से निकट-कोर क्षेत्र (संवेक्टिव-रेडिएटिव बाउंड्री, CRB के ठीक बाहर) पर केंद्रित था जहाँ ब्रंट-वैसाला आवृत्ति चरम पर होती है।

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

1. क्षेत्र ज्यामिति और शक्ति: सिमुलेशन ने स्वतः उत्पन्न (self-consistently) चुंबकीय क्षेत्र बनाए जहाँ निकट-कोर क्षेत्र में अजीमुथल (टोरोइडल) घटक, रेडियल (पोलॉइडल) घटक पर हावी रहता है। ब्रंट-वैसाला स्पाइक के शिखर पर, रूट-मीन-स्क्वायर अजीमुथल और रेडियल क्षेत्र का अनुपात ($B_\phi/B_r$) अधिकतम 20 तक पहुँच गया।

   • CRB पर: $B_r \approx 33.8$ kG, $B_\phi \approx 120.8$ kG।
   • ब्रंट-वैसाला स्पाइक के शिखर पर: $B_r \approx 0.29$ kG, $B_\phi \approx 7.65$ kG।
   • ये मान Takata et al. (2025) द्वारा KIC 9244992 के लिए अनुमानित न्यूनतम क्षेत्र शक्ति के कुछ गुना के भीतर हैं।

2. विभेदक घूर्णन और शियर: जबकि गोलाकार-औसत घूर्णन प्रोफाइल ने निकट-कोर क्षेत्र और सतह के बीच केवल एक छोटा अंतर ($\sim 3.7\%$) दिखाया—जो एस्टरोसेइस्मिक निष्कर्षों के अनुरूप है—सिमुलेशन ने महत्वपूर्ण अक्षांशीय (latitudinal) विभेदक घूर्णन का खुलासा किया। घूर्णन दर अक्षांशों केAcross (विशेष रूप से भूमध्य रेखा के पास) $20\%$ तक भिन्न हुई। यह स्थानीयकृत शियर, न कि वैश्विक रेडियल शियर, $\Omega$-इफेक्ट को संचालित करता है, जो पोलॉइडल घटक से टोरोइडल क्षेत्र को प्रवर्धित करता है।

3. अक्षांशीय संरचना: टोरोइडल क्षेत्र की शक्ति ने अक्षांश में डबल-पीक संरचना प्रदर्शित की, जिसमें मध्य-अक्षांशों पर अधिकतम और ध्रुवों पर न्यूनतम मान देखे गए। यह ज्यामिति Takata et al. (2025) द्वारा रिपोर्ट किए गए लेजेंड्रे पॉलीनोमियल गुणांकों ($a < 0, b > 0$) के चिह्न के साथ मेल खाती है।

4. मजबूती (Robustness): लेखकों ने सत्यापित किया कि परिणाम प्रारंभिक क्षेत्र शक्ति के आर्टिफैक्ट (artifacts) नहीं हैं (एक रन जिसमें प्रारंभिक क्षेत्र $1/5$ था, उसी स्थिर अवस्था तक पहुँचा) या विशिष्ट पैरामीटर विकल्पों के कारण नहीं हैं, बशर्ते कि रेनॉल्ड्स संख्या पर्याप्त रूप से उच्च ($Re > \sim 100$) हो और रॉस्बी संख्या उपयुक्त हो। संवेदनशीलता परीक्षणों ने दिखाया कि $Re$को$\sim 30$तक कम करने से काफी कमजोर क्षेत्र ($1-2$ kG) प्राप्त हुए, जो डायनेमो संतृप्ति (saturation) के लिए उच्च $Re$ की आवश्यकता को रेखांकित करता है।

महत्व और दावे
पेपर का दावा है कि निकट-कोर क्षेत्र में एक प्रभावी टोरोइडल चुंबकीय क्षेत्र की उपस्थिति कोर-कन्वक्टिंग तारों की एक सामान्य विशेषता है, जो उपयुक्त आयामहीन संख्याओं ($Re, Ro$) का उपयोग करने पर वैश्विक डायनेमो सिमुलेशन से स्वाभाविक रूप से उत्पन्न होती है।

लेखक इस बात पर जोर देते हैं कि उनके परिणाम KIC 9244992 के एस्टरोसेइस्मिक व्याख्या के लिए स्वतंत्र संख्यात्मक समर्थन प्रदान करते हैं। महत्वपूर्ण रूप से, यह अध्ययन दर्शाता है कि मजबूत टोरोइडल क्षेत्र कमजोर गोलाकार-औसत विभेदक घूर्णन के साथ सह-अस्तित्व में रह सकते हैं क्योंकि क्षेत्र का उत्पादन स्थानीयकृत अक्षांशीय शियर द्वारा संचालित होता है। यह निष्कर्ष बताता है कि भविष्य के एस्टरोसेइस्मिक अनुमानों को पूरे तारे में मुख्य रूप से रेडियल चुंबकीय क्षेत्र ज्यामिति का अनुमान नहीं लगाना चाहिए, क्योंकि निकट-कोर क्षेत्र संभवतः त्रि-आयामी शियर प्रवाह द्वारा उत्पन्न टोरोइडल क्षेत्रों द्वारा प्रभावी होता है। यह कार्य अवलोकन संबंधी बाधाओं और संख्यात्मक डायनेमो सिद्धांत के बीच के अंतर को पाटता है, यह पुष्टि करता है कि KIC 9244992 के लिए अनुमानित विशिष्ट क्षेत्र ज्यामिति स्टेलर मैग्नेटोहाइड्रोडायनामिक्स के ढांचे के भीतर भौतिक रूप से प्रशंसनीय है।