DOI: https://doi.org/10.1093/mnras/stag646
تاريخ النشر: 2026-04-09
المؤلف: Swapnil Shankar وآخرون
الموضوع الرئيسي: انفجارات أشعة غاما والسوبرنوفا
نظرة عامة
تقدم هذه الدراسة مجموعة شاملة من محاكاة الديناميكا الزمنية العامة ثلاثية الأبعاد (GRMHD) لانفجارات السوبرنوفا الناتجة عن انهيار النواة (CCSNe) لنجم أصلي كتلة 25 $M_\odot$. تم محاكاة ما مجموعه 12 نموذجًا على مدى 190-260 مللي ثانية، مع دراسة تأثير معدلات الدوران ($\Omega_0$) وقوة الحقول المغناطيسية ($B_0$) على تشكيل النفاثات من خلال آلية الدوران المغناطيسي. استخدمت المحاكاة كود GRaM-X المعزز بوحدات معالجة الرسوميات، وكشفت أن النماذج التي لديها $B_0 = 10^{11}$ G لم تحقق الانفجار، بينما تلك التي لديها $B_0 = 10^{12}$ G أظهرت أشكال نفاثات متنوعة ونتائج انفجارية تعتمد على $\Omega_0$. ومن الجدير بالذكر أن النماذج التي لديها $\Omega_0 \geq 2.0$ راديان/ثانية أنتجت نفاثات بسرعات تتجاوز 15,000 كم/ثانية، مما يجعلها مرشحة محتملة لنجم سوبرنوفا من النوع Ic ذو الخطوط العريضة.
تؤكد النتائج على الأدوار الحاسمة للدوران والحقول المغناطيسية في آلية CCSNe، مع تسليط الضوء على قيود البيانات الرصدية الحالية ونماذج تطور النجوم أحادية الأبعاد في تقييد ملفات الدوران والحقول المغناطيسية للنجوم الضخمة. تؤكد الدراسة على ضرورة إجراء تحقيقات منهجية في الظروف اللازمة للانفجار والمعلمات الناتجة، مدعومة بالكفاءة الحاسوبية لمحاكاة معززة بوحدات معالجة الرسوميات. تشير النتائج إلى أنه بينما تختلف طاقات الانفجار وكتل المواد المنبعثة بشكل كبير مع $\Omega_0$ و$B_0$، فإن النماذج التي لديها $B_0 = 10^{12}$ G و$\Omega_0 \geq 0.5$ راديان/ثانية تظهر زيادة في طاقات الانفجار وكتل المواد المنبعثة بمرور الوقت، مما يشير إلى تفاعل معقد بين هذه العوامل في تحديد ديناميات CCSNe. ستستكشف الأعمال المستقبلية كتل النجوم الأصلية المختلفة وكثافتها لتوضيح هذه العلاقات بشكل أكبر.
مقدمة
في هذه المقدمة، يناقش المؤلفون عمليات نهاية حياة النجوم الضخمة (كتلة أكبر من 10 $M_\odot$)، والتي تت culminate في انهيار جاذبي لنوى الحديد الخاصة بها، مما يؤدي إلى تشكيل نجم نيوتروني أولي (PNS) أو ثقب أسود (BH). يطلق هذا الانهيار حوالي $10^{53}$ إرج من الطاقة، بشكل رئيسي على شكل نيوترينوات، والتي يمكن أن تدفع انفجار سوبرنوفا ناتج عن انهيار النواة (CCSN). بينما تفسر الآلية المدفوعة بالنيوترينوات معظم CCSNe، فإن مجموعة فرعية من السوبرنوفا، المعروفة باسم سوبرنوفا من النوع Ic ذو الخطوط العريضة (SNe Ic-bl) أو هايبرنوفا، تظهر طاقات حركية حوالي $10^{52}$ إرج وتدفقات نسبية، مما يتطلب تفسيرات بديلة مثل الآلية المدفوعة بالنفاثات أو الدوران المغناطيسي.
يسلط المؤلفون الضوء على أهمية الدوران السريع والحقول المغناطيسية القوية في الآلية المدفوعة بالدوران المغناطيسي، والتي يمكن أن تطلق تدفقات شبيهة بالنفاثات قادرة على تفسير الطاقات الحركية العالية التي لوحظت في هايبرنوفا. يشيرون إلى أن عدم استقرار الدوران المغناطيسي (MRI) يمكن أن يعزز الحقول المغناطيسية الضعيفة إلى قوة مغناطيسية، لكن محاكاة CCSN النموذجية تفتقر إلى الدقة المكانية اللازمة لحل هذا عدم الاستقرار بالكامل. تقدم الورقة دراسة شاملة تشمل 12 محاكاة ثلاثية الأبعاد (3D) للديناميكا المائية المغناطيسية العامة (GRMHD) لنجم أصلي كتلة 25 $M_\odot$، مع تغيير منهجي لمعدلات الدوران الأولية والحقول المغناطيسية للتحقيق في الظروف اللازمة لتشكيل النفاثات بنجاح. يمثل هذا العمل، الذي يستخدم الكود المعزز بوحدات معالجة الرسوميات GRaM-X المطور حديثًا، أكبر مجموعة من محاكاة GRMHD ثلاثية الأبعاد لنجم أصلي واحد، مما يظهر إمكانية إجراء دراسات شاملة للمعلمات في CCSNe المدفوعة بالدوران المغناطيسي.
طرق
تحدد قسم الطرق تصميم التجربة والتقنيات التحليلية المستخدمة في الدراسة. استخدم الباحثون نهجًا كميًا، حيث استخدموا طرقًا إحصائية لتحليل البيانات المجمعة من تجارب مختلفة. تم تحديد وقياس المتغيرات الرئيسية، مما يضمن أن عملية جمع البيانات تتماشى مع معايير صارمة لتعزيز الموثوقية والصلاحية.
شمل التحليل تطبيق نماذج إحصائية متقدمة، بما في ذلك تحليل الانحدار واختبار الفرضيات، لتفسير النتائج بشكل فعال. كما استخدم الباحثون أدوات برمجية لتصور البيانات، مما سهل فهمًا أوضح للعلاقات بين المتغيرات المدروسة. بشكل عام، تم تصميم الإطار المنهجي لضمان نتائج قوية تساهم في الجسم المعرفي القائم في هذا المجال.
نتائج
في هذا القسم، يتم تقديم نتائج محاكاة انفجارات السوبرنوفا الناتجة عن انهيار النواة، مع التركيز على الديناميات بعد ارتداد النواة. تبدأ المحاكاة بنجم أصلي، مما يؤدي إلى زيادة سريعة في الكثافة المركزية من حوالي $10^9 \, \text{g/cm}^3$ إلى الكثافة النووية ($\sim 10^{14} \, \text{g/cm}^3$) وانخفاض في الانخفاض المركزي $\alpha$ من حوالي 1 إلى 0.84. يحدث ارتداد النواة عند حوالي $t_b \sim 164 – 167 \, \text{ms}$، متأثرًا بالحقول المغناطيسية الأولية ومعلمات الدوران. بعد الارتداد، يتم توليد صدمة ديناميكية مائية كروية ولكنها تتوقف عند حوالي 10 مللي ثانية بسبب فقدان الطاقة من التفاعلات مع المواد المتراكمة وتفكك مجموعة الحديد. تصل مكون الحقل المغناطيسي التوربيني، $B_{\text{tor}}$، إلى قيم قوة مغناطيسية، متغيرة مع معدلات الدوران، وتظهر النماذج ذات معدلات الدوران الأعلى حقول مغناطيسية أقوى.
تظل تطورات الإنتروبيا المحددة متسقة عبر النماذج حتى حوالي 80 مللي ثانية بعد الارتداد، مع زيادة الحد الأقصى للإنتروبيا من حوالي $6 k_b \, \text{baryon}^{-1}$ عند الارتداد إلى $12 k_b \, \text{baryon}^{-1}$ عند 80 مللي ثانية. يتم استخدام معامل البلازما-$\beta$، $\beta = P_{\text{fluid}} / P_{\text{magnetic}}$، لتقييم تأثير الحقول المغناطيسية على الديناميات. تشير المناطق التي لديها $\beta \lesssim 1$ إلى الهيمنة المغناطيسية، وهو أمر حاسم لتشكيل النفاثات المحتملة. من الجدير بالذكر أن النماذج ذات معدلات الدوران الأعلى (B12) تظهر انتقالًا إلى الهيمنة المغناطيسية وتطوير تدفقات ثنائية القطب غير متناظرة، بينما لا تحقق نماذج B11 تشكيل النفاثات ضمن إطار زمن المحاكاة. تشير النتائج إلى أنه بينما من غير المحتمل أن تنتج نماذج B11 انفجارات مدفوعة بالنفاثات على الفور، فإن نماذج B12 تظهر مجموعة من النتائج، حيث تنجح النماذج ذات الدوران الأسرع في تشكيل نفاثات تنتشر على طول محور الدوران، على الرغم من وجود عدم تماثل كبير بسبب عدم الاستقرار.
مناقشة
في هذا القسم، يناقش المؤلفون إعداد ونتائج محاكاة السوبرنوفا الناتجة عن انهيار النواة (CCSN) ثلاثية الأبعاد باستخدام نموذج نجم أصلي كتلة 25 $M_\odot$ في بداية عمره. يبنون على محاكاة سابقة أحادية البعد تشير إلى أن معامل كثافة نواة النجم الأصلي ($\xi_{2.5} = 0.47$) يشير إلى سيناريو انفجار صعب. يقدم المؤلفون مجموعة متنوعة من ملفات الدوران وتكوينات الحقول المغناطيسية لاستكشاف تأثير هذه المعلمات على ديناميات الانفجار. يتم نمذجة الدوران التفاضلي الأولي باستخدام ملف تعريف سرعة زاوية محددة، ويتم تحديد الحقل المغناطيسي لفحص تأثيراته على آلية CCSN.
تكشف المحاكاة أنه بينما لا تحقق النماذج ذات الحقول المغناطيسية الضعيفة (B11) انفجارات، فإن تلك ذات الحقول المغناطيسية الأعلى (B12) تظهر مجموعة من طاقات الانفجار وكتل المواد المنبعثة، متأثرة بشكل خاص بمعدلات الدوران الأولية. تظهر تطورات الصدمة اختلافات كبيرة بين النماذج، حيث تؤدي معدلات الدوران الأعلى إلى انتشار أسرع للصدمة وأحجام أكبر للصدمة. يشير المؤلفون إلى أن طاقات الانفجار لنماذج B12، على الرغم من زيادتها المستمرة، تظل أقل من القيمة الكانونية المتوقعة لسوبرنوفا من النوع Ic-bl، على الأرجح بسبب الكثافة العالية للنجم الأصلي. تؤكد الدراسة على ضرورة إجراء محاكاة ثلاثية الأبعاد منهجية لفهم أفضل لتفاعل الدوران والحقول المغناطيسية في CCSNe، مما يمهد الطريق للتحقيقات المستقبلية في الظروف التي تؤدي إلى انفجارات ناجحة.
DOI: https://doi.org/10.1093/mnras/stag646
Publication Date: 2026-04-09
Author(s): Swapnil Shankar et al.
Primary Topic: Gamma-ray bursts and supernovae
Overview
This research presents a comprehensive set of three-dimensional, dynamical-spacetime general relativistic magnetohydrodynamics (GRMHD) simulations of core-collapse supernovae (CCSNe) for a 25 $M_\odot$ progenitor. A total of 12 models were simulated over 190-260 ms, examining the influence of rotation rates ($\Omega_0$) and magnetic field strengths ($B_0$) on jet formation through the magnetorotational mechanism. The simulations utilized the GPU-accelerated GRaM-X code, revealing that models with $B_0 = 10^{11}$ G did not achieve explosion, while those with $B_0 = 10^{12}$ G exhibited diverse jet morphologies and explosive outcomes contingent on $\Omega_0$. Notably, models with $\Omega_0 \geq 2.0$ rad/s produced jets with velocities exceeding 15,000 km/s, marking them as potential progenitors for broad-lined type Ic supernovae.
The findings underscore the critical roles of rotation and magnetic fields in the CCSNe mechanism, while highlighting the limitations of current observational data and one-dimensional stellar evolution models in constraining the rotation profiles and magnetic fields of massive stars. The study emphasizes the necessity for systematic investigations into the conditions for explosion and the resultant parameters, facilitated by the computational efficiency of GPU-accelerated simulations. The results indicate that while the explosion energies and ejecta masses vary significantly with $\Omega_0$ and $B_0$, models with $B_0 = 10^{12}$ G and $\Omega_0 \geq 0.5$ rad/s show increasing explosion energies and ejecta masses over time, suggesting a complex interplay between these factors in determining the dynamics of CCSNe. Future work will explore different progenitor masses and compactness to further elucidate these relationships.
Introduction
In this introduction, the authors discuss the end-of-life processes of massive stars (mass greater than 10 $M_\odot$), which culminate in the gravitational collapse of their iron cores, leading to the formation of a protoneutron star (PNS) or black hole (BH). This collapse releases approximately $10^{53}$ erg of energy, predominantly in the form of neutrinos, which can drive a core-collapse supernova (CCSN) explosion. While the neutrino-driven mechanism accounts for most CCSNe, a subset of supernovae, known as broad-lined type-Ic supernovae (SNe Ic-bl) or hypernovae, exhibit kinetic energies around $10^{52}$ erg and relativistic outflows, necessitating alternative explanations such as the jet-driven or magnetorotational mechanism.
The authors highlight the importance of rapid rotation and strong magnetic fields in the magnetorotational mechanism, which can launch jet-like outflows capable of explaining the high kinetic energies observed in hypernovae. They note that the magnetorotational instability (MRI) can amplify weak magnetic fields to magnetar strength, but typical CCSN simulations lack the spatial resolution needed to fully resolve this instability. The paper presents a comprehensive study involving 12 three-dimensional (3D) general relativistic magnetohydrodynamic (GRMHD) simulations of a 25 $M_\odot$ progenitor, systematically varying initial rotation rates and magnetic fields to investigate conditions for successful jet formation. This work, utilizing the newly developed GPU-accelerated code GRaM-X, represents the largest set of 3D GRMHD simulations for a single progenitor, demonstrating the feasibility of extensive parameter studies in magnetorotationally driven CCSNe.
Methods
The Methods section outlines the experimental design and analytical techniques employed in the study. The researchers utilized a quantitative approach, employing statistical methods to analyze the data collected from various experiments. Key variables were identified and measured, ensuring that the data collection process adhered to rigorous standards to enhance reliability and validity.
The analysis involved the application of advanced statistical models, including regression analysis and hypothesis testing, to interpret the results effectively. The researchers also employed software tools for data visualization, which facilitated a clearer understanding of the relationships between the variables studied. Overall, the methodological framework was designed to ensure robust findings that contribute to the existing body of knowledge in the field.
Results
In this section, the results of simulations of core collapse supernovae are presented, focusing on the dynamics following core bounce. The simulations begin with a progenitor star, leading to a rapid increase in central density from approximately $10^9 \, \text{g/cm}^3$ to nuclear density ($\sim 10^{14} \, \text{g/cm}^3$) and a decrease in the central lapse $\alpha$ from about 1 to 0.84. Core bounce occurs at approximately $t_b \sim 164 – 167 \, \text{ms}$, influenced by initial magnetic field and rotation parameters. Postbounce, a spherical magnetohydrodynamic shock is generated but stalls at around 10 ms due to energy loss from interactions with accreting material and iron group dissociation. The magnetic field’s toroidal component, $B_{\text{tor}}$, reaches magnetar-strength values, varying with rotation rates, and models with higher rotation rates exhibit stronger magnetic fields.
The evolution of specific entropy remains consistent across models until about 80 ms postbounce, with maximum entropy increasing from approximately $6 k_b \, \text{baryon}^{-1}$ at bounce to $12 k_b \, \text{baryon}^{-1}$ at 80 ms. The plasma-$\beta$ parameter, $\beta = P_{\text{fluid}} / P_{\text{magnetic}}$, is utilized to assess the influence of magnetic fields on dynamics. Regions with $\beta \lesssim 1$ indicate magnetic dominance, which is crucial for potential jet formation. Notably, models with higher rotation rates (B12) show a transition to magnetic dominance and the development of asymmetric bipolar outflows, while B11 models do not achieve jet formation within the simulation timeframe. The findings suggest that while B11 models are unlikely to produce prompt jet-driven explosions, B12 models exhibit a range of outcomes, with the most rapidly rotating models successfully forming jets that propagate along the rotation axis, albeit with significant asymmetries due to instabilities.
Discussion
In this section, the authors discuss the setup and findings of their three-dimensional magnetorotational core-collapse supernova (CCSN) simulations using a 25 $M_\odot$ zero-age main sequence progenitor model. They build upon previous one-dimensional simulations that indicated the progenitor’s core compactness parameter ($\xi_{2.5} = 0.47$) suggests a challenging explosion scenario. The authors introduce various rotation profiles and magnetic field configurations to explore the impact of these parameters on the explosion dynamics. The initial differential rotation is modeled using a specific angular velocity profile, and the magnetic field is parameterized to examine its effects on the CCSN mechanism.
The simulations reveal that while models with lower magnetic field strengths (B11) do not achieve explosions, those with higher magnetic fields (B12) exhibit a range of explosion energies and ejecta masses, particularly influenced by the initial rotation rates. The shock evolution shows significant differences between models, with higher rotation rates leading to more rapid shock propagation and larger shock radii. The authors note that the explosion energies for the B12 models, although increasing steadily, remain below the canonical value expected for type Ic-bl CCSNe, likely due to the progenitor’s high compactness. The study emphasizes the necessity of systematic 3D simulations to better understand the interplay of rotation and magnetic fields in CCSNe, paving the way for future investigations into the conditions that lead to successful explosions.