ExTraSS: a Domain Decomposed 3D NLTE Radiative Transfer spectral synthesis code for nebular phase transients

ملخص تقني لبرنامج ExTraSS: كود تخليق طيفي لنقل الإشعاع في حالة عدم التوازن المحلي الحراري (NLTE) ثلاثي الأبعاد يعتمد على تقسيم النطاقات

بيان المشكلة
يمثل نمذجة مرحلة السديم (nebular phase) للمستعمرات النجمية (SNe) تحديًا حاسوبياً كبيراً. في هذه المرحلة، تنخفض كثافات المقذوفات (ejecta) بدرجة كافية ليصبح الهيكل شفافاً بصرياً إلى حد كبير، مما يعني أن الحجم بأكله يساهم في الانبعاث. ومع ذلك، تتطلب النمذجة الدقيقة فيزياء عدم التوازن المحلي الحراري (NLTE) لمراعاة الميكروفيزياء المعقدة، بينما لا يمكن تجاهل تأثيرات نقل الإشعاع (RT)، لا سيما بالنسبة للمستعمرات النجمية الغنية بالهيدروجين والناتجة عن انهيار القلب. علاوة على ذلك، تستوجب المورفولوجيا غير المتماثلة للمستعمرات النجمية استخدام نماذج مدخلات ثلاثية الأبعاد.

تتمثل العقبة الرئيسية في عمليات حساب NLTE ثلاثية الأبعاد في توليد وتخزين معدلات الاستثارة الضوئية. فبالنسبة لشبكة تحتوي على حوالي $10^5$ خلية، قد يتطلب تخزين هذه المعدلات لملايين الانتقالات ما يزيد عن 1 تيرابايت من البيانات، وهو ما يتجاوز سعة الذاكرة للعقد الحسابية المنفردة في الحواسيب الفائقة الحديثة. كانت الجهود السابقة في مجال 3D NLTE إما محدودة بنماذج أحادية البعد، أو استخدمت تقريبات الشفافية البصرية، أو وظفت معالجات "في نفس الموقع" (on-the-spot) التي تهمل تأثيرات مجال الإشعاع العالمي. كانت هناك حاجة لكود قادر على إجراء نقل إشعاعي كامل ثلاثي الأبعاد بنظام NLTE مع اقتران عالمي بين المجال الإشعاعي وتجمعات المستويات الذرية.

المنهجية
يقدم المؤلفون ExTraSS (محاكي الطيف العابر الانفجاري)، وهو كود تم تطويره للقيام بنقل إشعاعي كامل ثلاثي الأبعاد بنظام NLTE. يعمل الكود في وضع "اللقطة" (snap-shot)، حيث يتبع انتقال الفوتونات بشكل مستقل عن الزمن، وهو ما يناسب المراحل المتأخرة حيث تكون الأعماق البصرية منخفضة.

• الشبكة والهندسة: يستخدم ExTraSS شبكة إحداثيات كروية، مما يسمح بدقة عالية في المقذوفات الداخلية الأبطل حركة وخلايا أكبر عند أنصاف الأقطار الخارجية. يتم تقسيم التباين الزاوي لمجال الإشعاع باستخدام $N_{\theta, VA} \times N_{\phi, VA}$ من المشاهدين (الافتراضي 20 $\times$ 20).
• نقل أشعة غاما ($\gamma$): يتتبع الكود الاضمحلال الإشعاعي (بشكل أساسي $^{56}$Ni و $^{56}$Co) باستخدام طريقة مونت كارلو مع تشتت كومبتون. ويُفترض أن البوزيترونات محتجزة محلياً. يقوم ترسب الطاقة من أشعة غاما بدفع التسخين غير الحراري، والتأين، والاستثارة عبر صيغة سبنسر-فانو (Spencer-Fano).
• حل NLTE: يحدد الحل بشكل تكراري درجة الحرارة المحلية وتجمعات المستويات. وهو يربط الفيزياء الحرارية وغير الحرارية، موازناً بين التسخين (من أشعة غاما والتأين الضوئي) والتبريد (الانتقالات بين المستويات المرتبطة، والارتباط الحر، وانبعاث الفوتون المزدوج). يستخدم الحل مخطط نيوتن-رافسون مع تخميد تكيفي للوصول إلى التقارب في تجمعات المستويات.
• نقل الإشعاع (RT): ينفذ الكود تقنية تتبع الأشعة الكاملة لحساب المجال الإشعاعي العالمي. تنبعث الأشعة من مراكز الخلايا باتجاه جميع زوايا الرؤية. يحسب الكود الأعماق البصرية للتأين الضوئي (باستخدام تقريبات هيدروجينية للحالات المثارة وبيانات Verner et al للحالات الأرضية) والاستثارة الضوئية (باستخدام تقريب سوبولف لتفاعلات الخطوط).
• تقسيم النطاقات (Domain Decomposition): للتغلب على قيود الذاكرة، يستخدم ExTraSS خوارزمية تقسيم النطاقات، حيث يتم تقسيم النطاق الحسابي إلى شرائح تعتمد على الزاوية السمتية ($\phi$)، حيث تتعامل كل عقدة مع شريحة محددة.
  • بنية الاتصال: تُعين كل عقدة نواة "مديرة" واحدة للتعامل مع الاتصال بين العقد، بينما تقوم بقية النوى "العاملة" بتتبع الأشعة والحسابات المحلية.
  • التحسين: يستخدم النظام اتصالات MPI غير حاصرة (non-blocking) لتقليل وقت خمول النوى العاملة. يتم نقل البيانات في حزم (أشعة) بدلاً من نقلها بشكل فردي. تقوم النوى المديرة بتنسيق نقل الأشعة بين النطاقات، مما يضمن قضاء النوى العاملة أقصى وقت في الحساب بدلاً من انتظار البيانات.

المساهمات الرئيسية

1. أول نظام 3D NLTE RT للمستعمرات النجمية: يمثل هذا العمل، حسب علم المؤلفين، أول تنفيذ لتقسيم النطاقات خصيصاً لكود نقل إشعاعي للمستعمرات النجمية، مما يتيح حسابات NLTE ثلاثية الأبعاد كاملة كانت غير ممكنة سابقاً بسبب قيود الذاكرة.
2. تحسين الذاكرة: تفصل الورقة تنفيذ الذاكرة المشتركة، و"المصفوفات المسننة" (تخصيص الذاكرة فقط للأنواع/المستويات الموجودة)، ومخطط تقسيم النطاقات، الذي يجعل استهلاك الذاكرة يتناسب عكسياً مع عدد العقد المستخدمة.
3. التحسينات الخوارزمية: قام المؤلفون بتحسين وتيرة استدعاء روتين سبنسر-فانو المكلف حاسوبياً (باستخدام مخطط 4+10: يُستدعى في التكرارات 1، 4، ثم كل 10 تكرارات بعد ذلك) وحسنوا مخططات التخميد لمحلل تجمعات المستويات بنيوتن-رافسون لتحقيق التوازن بين سرعة التقارب والاستقرار.

النتائج والتحقق
تحقق المؤلفون من صحة ExTraSS من خلال عدة اختبارات مرجعية:

• التقارب والكفاءة: أظهرت الاختبارات على حاسوب Dardel (باستاستخدام عقد بـ 128 نواة) أن الكود يحقق كفاءة حاسوبية تبلغ $\sim 0.8$ لوحدة نقل الإشعاع عند تحسين توازن الأحمال. يتم الوصول إلى التقارب العالمي (نسبة اللمعان البولتومتري إلى ترسب أشعة غاما) في غضون $\sim 17$ تكراراً، وتستقر عادةً حول نسبة 0.96–1.0.
• المقارنة مع SUMO (أحادي البعد): تمت مقارنة ExTraSS مع الكود أحادي البعد SUMO باستخدام نموذج "أحادي البعد المحول لثلاثي الأبعاد" (J18).
  • الظروف الفيزيائية: اتفقت درجات الحرارة والكسور الإلكترونية بشكل جيد، رغم أن ExTraSS أظهر درجات حرارة أعلى قليلاً في المقذوفات الداخلية بسبب انخفاض الكسور المتعادلة.
  • المعدلات: كانت معدلات التأين الضوئي والاستثارة الضوئية متوافقة عموماً، وإن أظهر ExTraSS معدلات أعلى قليلاً في المناطق الأكثر داخلية.
  • الأطياف: أظهرت الأطياف الاصطناعية توافقاً جيداً، مع وجود اختلافات طفيفة في انبعاث Fe II وقوة خط $H\alpha$، والتي عُزيت إلى الاختلافات في عدد المستويات الذرية المدرجة (ExTraSS يضع حداً أقصى بـ 100 مستوى لكل أيون مقابل مئات في SUMO).
• المقارنة مع اختبار مرجعي متعدد الأكواد (Blondin et al. 2022): تم اختبار ExTraSS مقابل نموذج 'toy06' جنباً إلى جنب مع ARTIS-nebular وCMFGEN وSUMO.
  • الظروف الفيزيائية: حسب ExTraSS درجات حرارة أقل قليلاً في المقذوفات الداخلية وكسور تأين كوبالت أقل مقارنة ببعض الأكواد، لكن النتائج كانت متسقة عموماً مع المجموعة.
  • الأطياف: اتفقت المخرجات الطيفية مع SUMO والأكواد الأخرى، خاصة في نطاقات 5100–6400 Å و >8800 Å. لوحظت اختلافات في ذرى انبعاث الحديد (Fe) حول 4600–5000 Å، حيث توقع ExTraSS لمعاناً أقل، وهو ما قد يعود إلى بتر تجمعات المستويات أو الاختلافات في التأين.

الأهمية
تزعم الورقة أن Ex-TraSS أصبح الآن كوداً وظيفياً كاملاً لتخليق الطيف بنظام نقل إشعاعي 3D NLTE يعتمد على تقسيم النطاقات. تكمن أهميته في قدرته على التعامل مع متطلبات الذاكرة الهائلة لفيزياء 3D NLTE، مما يفتح الطريق لنمذجة دقيقة للغاية للمستعمرات النجمية غير المتماثلة في مرحلة السديم. وبينما يُطبق حالياً على نماذج المستعمرات النجمية الناتجة عن انهيار القلب (تحديداً النوع الثاني)، يذكر المؤلفون أن الكود مصمم ليكون قابلاً للتطبيق على مجموعة واسعة من نماذج المدخلات والظواهر العابرة الانفجارية، بشرط دمج مكونات فيزيائية إضافية (مثل الاعتماد على الزمن، أو الغبار، أو الجزيئات) في التحديثات المستقبلية. يضع هذا العمل أساساً لفهم أفضل لخصائص الانفجار، وكتل المقذوفات، والتركيبات من خلال العلامات الطيفية التفصيلية.