Sunspot simulations with MURaM -- I. Parameter study using potential field initial conditions

यहाँ Schmassmann et al. (2026) के शोध पत्र "Sunspot simulations with MURaM I. Parameter study using potential field initial conditions" का विस्तृत तकनीकी सारांश दिया गया है।

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

सूर्यकलंक (sunspots) के पिछले मैग्नेटोहाइड्रोडायनामिक (MHD) सिमुलेशन सटीक चुंबकीय क्षेत्र वितरण को पुनरुत्पादित करने में संघर्ष करते रहे हैं। मौजूदा मॉडलों (जैसे Rempel 2012) में एक प्राथमिक सीमा कृत्रिम ऊपरी सीमा स्थितियों (artificial top boundary conditions) पर निर्भरता है, जो चुंबकीय क्षेत्र को फोटोसफेरिक स्तर पर अधिक क्षैतिज होने के लिए मजबूर करती हैं।

• कृत्रिम सीमाओं के परिणाम: ये स्थितियाँ ऐसे सूर्यकलंक उत्पन्न करती हैं जिनका उम्ब्रा-पेनुम्ब्रा सीमा झुकाव $\approx 60^\circ$ होता है (जो प्रेक्षित $\approx 45^\circ$ की तुलना में काफी अधिक क्षैतिज है) और उम्ब्रा सीमा पर ऊर्ध्वाधर चुंबकीय क्षेत्र घटक ($B_z$) प्रेक्षित मानों ($\approx 1.8$ kG) से कम होता है।
• लक्ष्य: लेखक कृत्रिम ऊपरी सीमा बल के बिना, बेहतर सूर्यकलंक चुंबकीय और गतिशील गुणों को पुनरुत्पादित करने के लिए वैकल्पिक प्रारंभिक स्थितियों और पैरामीटर विविधताओं का पता लगाने का लक्ष्य रखते हैं, विशेष रूप से पेनुम्ब्रा के निर्माण और एवर्सहेड प्रवाह (Evershed flow) पर ध्यान केंद्रित करते हुए।

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

यह अध्ययन MURaM (MPS/University of Chicago Radiative Magneto-hydrodynamics) कोड का उपयोग करता है। लेखकों ने स्व-समान (self-similar) या फ्लक्स-उद्भव (flux-emergence) सेटअप के बजाय पोटेंशियल फील्ड प्रारंभिक स्थितियों (Nordlund 2015 द्वारा प्रस्तावित) का उपयोग करके एक व्यापक पैरामीटर अध्ययन किया।

सिमुलेशन सेटअप

• प्रारंभिक स्थितियाँ: बॉक्स के निचले भाग में एक गॉसियन वितरण ($B_{z,0}$) से एक पोटेंशियल चुंबकीय क्षेत्र निकाला गया था। स्पॉट के बाहर लुप्त चुंबकीय क्षेत्र की नकल करने के लिए, एक समान विपरीत ऊर्ध्वाधर क्षेत्र ($B_{opp}$) को घटाया गया।
• सीमा स्थितियाँ (Boundary Conditions):
  • ऊपरी (Top): इनफ्लो (inflows) के लिए खुला, आउटफ्लो (outflows) के लिए बंद, जिसमें चुंबकीय क्षेत्र को पोटेंशियल फील्ड के रूप में निकाला गया था (कोई कृत्रिम झुकाव बल नहीं)।
  • निचली (Bottom): आम तौर पर द्रव्यमान प्रवाह (mass flows) के लिए खुला, जिसमें चुंबकीय क्षेत्र को सममित रूप से निकाला गया था। एक विशिष्ट रन ('160_opp000b') ने चुंबकीय ट्रंक के भीतर एक बंद निचली सीमा का परीक्षण किया।
• परिवर्तित पैरामीटर:
  • प्रारंभिक चुंबकीय क्षेत्र शक्ति ($B_0$): 20, 40, 80, और 160 kG।
  • चुंबकीय फ्लक्स: गॉसियन फ्लक्स ($F_{Gauss}$) और कुल फ्लक्स ($F_{tot}$) को $10^{22}$ Mx के आसपास बदला गया।
  • विपरीत क्षेत्र ($B_{opp}$): 0 से 300 G।
  • बॉक्स की चौड़ाई: 49.152 Mm (मानक) बनाम 98.304 Mm (चौड़ा)।
  • रेज़ोल्यूशन (Resolution): निम्न रेज़ोल्यूशन (96 km क्षैतिज / 32 km ऊर्ध्वाधर) बनाम उच्च रेज़ोल्यूशन (32 km / 16 km)।
  • रेडिएटिव ट्रांसफर: ग्रे (निम्न रेज़ोल्यूशन) बनाम 4 ओपेसिटी बिन्स के साथ नॉन-ग्रे (उच्च रेज़ोल्यूशन)।

डेटा प्रोसेसिंग

• सिमुलेशन का विश्लेषण $\approx 4$ घंटों के संदर्भ समय (या उच्च-रेज़ोल्यूशन रन के लिए विशिष्ट बाद के समय) पर किया गया।
• तीव्रता, चुंबकीय क्षेत्र वेक्टर्स ($B_z, B_r, |B|$), झुकाव, और वेगों ($v_r, v_z$) के अज़ीमुथल औसत (azimuthal averages) की गणना की गई।
• उम्ब्रा और स्पॉट की सीमाओं को तीव्रता कंटूर ($I = 0.5 I_{qs}$ और $I = 0.9 I_{qs}$) द्वारा परिभाषित किया गया।

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

1. व्यवस्थित पैरामीटर अध्ययन: यह MURaM कोड में विभिन्न $B_0$, रेज़ोल्यूशन और बॉक्स आकार के साथ पोटेंशियल फील्ड प्रारंभिक स्थितियों का उपयोग करके सूर्यकलंक निर्माण का पहला व्यापक पैरामीटर अध्ययन है।
2. रेज़ोल्यूशन निर्भरता: यह अध्ययन पेनुम्ब्रा की सूक्ष्म संरचना और प्रवाह पैटर्न को पुनरुत्पादित करने में संख्यात्मक रेज़ोल्यूशन की महत्वपूर्ण भूमिका को रेखांकित करता है।
3. पोटेंशियल प्रारंभिक स्थितियों का सत्यापन: यह प्रदर्शित करता है कि पोटेंशियल फील्ड प्रारंभिक स्थितियाँ, विशिष्ट फ्लक्स और क्षेत्र शक्ति मापदंडों के साथ मिलकर, कृत्रिम ऊपरी सीमा बल के बिना यथार्थवादी सूर्यकलंक आकारिकी (morphologies) प्रदान कर सकती हैं।
4. फ्लक्स उद्भव गतिशीलता (Flux Emergence Dynamics): पेपर देखे गए द्वि-दिशीय प्रवाहों (bi-directional flows) को चल रहे फ्लक्स उद्भव से जोड़ता है, जो प्रारंभिक पेनुम्ब्रा निर्माण के लिए एक भौतिक तंत्र प्रदान करता है।

4. मुख्य परिणाम (Key Results)

A. आकारिकी और आकार (Morphology and Size)

• क्षेत्र शक्ति ($B_0$): $B_0$ को 20 kG से बढ़ाकर 160 kG करने से स्पॉट, उम्ब्रा और पेनुम्ब्रा का आकार व्यवस्थित रूप से बढ़ता है।
  • $B_0 \le 40$ kG: कोई पेनुम्ब्रा नहीं बनता; केवल इनफ्लो के साथ लंबे उत्तोलित संवेदी सेल (elongated convective cells) बनते हैं।
  • $B_0 = 80$ kG: छोटे पेनुम्ब्रा फिलामेंट्स के साथ मध्यवर्ती अवस्था।
  • $B_0 = 160$ kG: लंबे, पतले पेनुम्ब्रा फिलामेंट्स बनते हैं।
• पेनुम्ब्रा अंश (Penumbra Fraction): यहाँ तक कि सबसे अच्छे रन ($B_0=160$ kG) भी देखे गए स्थिर सूर्यकलंक ($\approx 60\%$) की तुलना में कम पेनुम्ब्रा अंश ($r_{pu}/r_s \approx 47-52\%$) उत्पन्न करते हैं। यह सुझाव देता है कि सिमुलेशन सूर्यकलंक निर्माण के प्रारंभिक चरणों का प्रतिनिधित्व करते हैं न कि परिपक्व, स्थिर स्पॉट का।
• फ्लक्स सामग्री: $F_{Gauss} > 10^{22}$ Mx वाले रन उम्ब्रा में अवास्तविक रूप से मजबूत ऊर्ध्वाधर क्षेत्र ($B_z > 5$ kG) उत्पन्न करते हैं। सबसे यथार्थवादी परिणाम $F_{Gauss} \approx 10^{22}$ Mx से प्राप्त होते हैं।

B. चुंबकीय क्षेत्र के गुण (Magnetic Field Properties)

• ऊर्ध्वाधर क्षेत्र ($B_z$): सबसे अच्छे रन ($160$ kG, $10^{22}$ Mx) एक $B_z$ प्रोफाइल उत्पन्न करते हैं जो त्रिज्या के साथ यथार्थवादी रूप से घटता है, हालांकि अधिकतम केंद्रीय क्षेत्र ($\approx 4.2-4.5$ kG) सामान्य अवलोकनों की तुलना में अधिक रहता है।
• झुकाव (Inclination): झुकाव कोण त्रिज्या के साथ बढ़ता है, जो $B_0=160$ kG के लिए स्पॉट सीमा पर $\approx 80^\circ$ तक पहुँच जाता है, जो कृत्रिम-सीमा मॉडलों की तुलना में अवलोकनों के अधिक करीब है।
• विपरीत क्षेत्र और बॉक्स का आकार: $B_{opp}$ को बदलना या बॉक्स की चौड़ाई बढ़ाना स्पॉट के बाहर के बैकग्राउंड क्षेत्र को प्रभावित करता है लेकिन आंतरिक स्पॉट की गतिशीलता या संरचना पर बहुत कम प्रभाव डालता है।

C. सतही प्रवाह (Surface Flows - Dynamics)

• निम्न रेज़ोल्यूशन ($96/32$ km):
  • $B_0 \le 40$ kG: शुद्ध इनफ्लो और डाउनफ्लो (कोई एवर्सहेड प्रवाह नहीं)।
  • $B_0 \ge 80$ kG: द्वि-दिशीय प्रवाह (Bi-directional flows) दिखाई देते हैं। आंतरिक पेनुम्ब्रा में मजबूत इनफ्लो और डाउनफ्लो दिखाई देते हैं, जबकि बाहरी पेनुम्ब्रा में आउटफ्लो दिखाई देता है। यह अवलोकनों (काउंटर-एवर्सहेड प्रवाह) की नकल करता है जो प्रारंभिक पेनुम्ब्रा निर्माण के लक्षण हैं।
• उच्च रेज़ोल्यूशन ($32/16$ km):
  • सिमुलेशन द्वि-दिशीय प्रवाह वाले फिलामेंट्स और शुद्ध एवर्सहेड प्रवाह (पूरे फिलामेंट की लंबाई के साथ रेडियल आउटफ्लो) प्रदर्शित करने वाले फिलामेंट्स का मिश्रण उत्पन्न करते हैं।
  • आउटफ्लो की गति उच्च-रेज़ोल्यूशन रन में अधिक होती है, जो Rempel (2012) के पिछले निष्कर्षों के अनुरूप है।
• फ्लक्स उद्भव (Flux Emergence): अधिकांश सफल रन निरंतर फ्लक्स उद्भव ($\approx 8-13\%$ प्रति घंटा) दिखाते हैं। अध्ययन बताता है कि स्पॉट के बाहरी किनारों पर फ्लक्स उद्भव द्वि-दिशीय प्रवाह को संचालित करता है और यह स्थिर पेनुम्ब्रा गठन का पूर्वगामी (precursor) है।

D. विफल विन्यास (Failed Configurations)

• बंद निचली सीमा (Closed Bottom Boundary): चुंबकीय ट्रंक के नीचे द्रव्यमान प्रवाह को रोकने से यथार्थवादी पेनुम्ब्रा और आउटफ्लो के निर्माण में बाधा आती है।
• अत्यधिक फ्लक्स: $10^{22}$ Mx से अधिक फ्लक्स अवास्तविक रूप से मजबूत चुंबकीय क्षेत्र उत्पन्न करते हैं।

5. महत्व और निष्कर्ष (Significance and Conclusions)

• सर्वश्रेष्ठ विन्यास: सबसे यथार्थवादी सूर्यकलंक सिमुलेशन एक पोटेंशियल फील्ड, $\approx 10^{22}$ Mx के कुल फ्लक्स, 160 kG की मजबूत प्रारंभिक क्षेत्र शक्ति, और एक खुली निचली सीमा के साथ प्राप्त किए जाते हैं।
• भौतिक अंतर्दृष्टि: सिमुलेशन बताते हैं कि प्रारंभिक पेनुम्ब्रा चरणों में देखे गए द्वि-दिशीय प्रवाह फ्लक्स उद्भव द्वारा संचालित होते हैं। स्थिर, शुद्ध एवर्सहेड प्रवाह में संक्रमण के लिए उच्च संख्यात्मक रेज़ोल्यूशन और/या प्राकृतिक चुंबकीय कैनोपी (canopy) के निर्माण की आवश्यकता हो सकती है (जो इन रन में पूरी तरह से प्राप्त नहीं हुआ था)।
• भविष्य की दिशाएँ: लेखक अनुमान लगाते हैं कि सूर्यकलंक को बहुत बड़े सिमुलेशन डोमेन में स्थापित करने से, जिससे प्राकृतिक लघु-स्तरीय डायनेमो क्रिया और कैनोपी निर्माण संभव हो सके, निम्न रेज़ोल्यूशन पर भी पूर्ण विकसित एवर्सहेड प्रवाह का विकास किया जा सके।
• कुल प्रभाव: यह कार्य पोटेंशियल फील्ड प्रारंभिक स्थिति दृष्टिकोण को मान्य करता है और यह पहचानता है कि संख्यात्मक रेज़ोल्यूशन MHD सिमुलेशन में पूरी तरह विकसित सूर्यकलंक पेनुम्ब्रा को पुनरुत्पादित करने में एक महत्वपूर्ण बाधा है।