DOI: https://doi.org/10.1103/fzhd-hcpp
تاريخ النشر: 2026-03-05
المؤلف: Danat Issa وآخرون
الموضوع الرئيسي: الفيزياء الفلكية والظواهر الكونية
نظرة عامة
تناقش هذه الفقرة دور الكولابسار – النجوم الضخمة التي تنهار إلى ثقوب سوداء (BHs) – في توليد انفجارات أشعة غاما طويلة المدة (LGRBs) من خلال نفاثات نسبية. تتأثر فعالية هذه النفاثات بقوة المجال المغناطيسي ودوران الثقب الأسود المركزي. أشارت الدراسات السابقة إلى أن الأقراص المغناطيسية المحتجزة غير الإشعاعية (MADs) يمكن أن تؤدي إلى سرعات توازن منخفضة ($a_{eq} = 0.035 – 0.07$)، مما ينتج عنه نفاثات ضعيفة تكافح للهروب من المادة النجمية. ومع ذلك، يقدم المؤلفون نتائج جديدة من محاكاة ديناميكية مغناطيسية عامة نسبية ثلاثية الأبعاد لنقل النيوترينو، تُظهر أن الأقراص المبردة بالنيوترينو يمكن أن تحافظ على سرعات ثقب أسود أعلى، تصل إلى سرعة توازن تبلغ حوالي $a_{eq} \approx 0.13$.
تعتبر هذه الزيادة في السرعة مهمة لأنها تتماشى مع قيم السرعة المستنتجة من ملاحظات LIGO/Virgo/KAGRA، وهي أكبر بمقدار 2-4 مرات من تلك الناتجة عن MADs غير الإشعاعية. وبالتالي، فإن النفاثات المنتجة تكون أقوى بمقدار 4-16 مرة، مما يعزز قدرتها على اختراق النجم السلف. كما تشير الدراسة إلى أنه بالنسبة لفترات LGRB النموذجية ($t \gtrsim 30$ s)، يمكن للثقوب السوداء تحقيق هذه السرعة الأعلى بنهاية الانفجار، بينما قد تترك LGRBs ذات معدلات تراكم أقصر أو أقل كتلة وراءها ثقوب سوداء تدور بسرعة أكبر.
مقدمة
تناقش مقدمة الورقة انفجارات أشعة غاما طويلة المدة (LGRBs)، وهي ظواهر فلكية قوية مرتبطة بانهيار النجوم الضخمة إلى ثقوب سوداء (BHs). يُعتقد أن هذه الانفجارات مدفوعة بالكولابسار، حيث يتشكل ثقب أسود وإذا كان الغلاف النجمي يحتوي على زخم زاوي كافٍ، فإنه يولد نفاثات نسبية من خلال قرص تراكم مغناطيسي. يتم تسليط الضوء على آلية بلاندفورد-زنايك (BZ) كعملية رئيسية لقوة النفاثات، مع كون اللمعان متناسبًا مع مربع دوران الثقب الأسود ($L_{BZ} \propto a^2 \Phi^2$)، حيث $a$ هو دوران الثقب الأسود و$\Phi$ هو التدفق المغناطيسي عند أفق الثقب الأسود. تشير ملاحظات موجات الجاذبية الأخيرة إلى أن معظم الثقوب السوداء لديها سرعات منخفضة ($a \leq 0.3$)، مما يثير تساؤلات حول خصائص دوران ثقوب الكولابسار التي تولد GRBs.
تؤكد الورقة على دور المجالات المغناطيسية في تشكيل النفاثات، وخاصة ضرورة وجود حالة قرص مغناطيسي محتجز (MAD)، مما يسمح بإطلاق النفاثات بكفاءة. ومع ذلك، تشير أيضًا إلى أن النماذج الحالية لم تأخذ في الاعتبار تأثيرات تبريد النيوترينو، التي تعتبر مهمة في سياق النجوم الضخمة المنهارة. يمكن أن يغير تبريد النيوترينو ديناميات القرص وتطور دوران الثقب الأسود، مما قد يؤدي إلى سرعات توازن أعلى مما كان يُعتقد سابقًا. يهدف المؤلفون إلى التحقيق في هذه التأثيرات من خلال سلسلة من محاكاة ديناميكية مغناطيسية عامة نسبية ثلاثية الأبعاد (GRMHD) التي تتضمن نقل النيوترينو، ساعين في النهاية لفهم كيف تؤثر هذه الديناميات على دوران الثقب الأسود، وبالتالي، نشاط GRB. توضح الورقة هيكلها، مشيرة إلى أن الأقسام اللاحقة ستفصل نموذج تطور الدوران، والأساليب العددية، ونتائج المحاكاة، والآثار على ملاحظات GRB.
النتائج
تقدم هذه الفقرة النتائج العددية للدراسة، مع تسليط الضوء على النتائج الرئيسية المستمدة من التحليل الحسابي. تشير النتائج إلى أن النموذج المقترح يظهر تحسينات كبيرة في الدقة والكفاءة مقارنة بالأساليب الحالية. على وجه التحديد، تكشف المحاكاة العددية أن النموذج يحقق معدل تقارب قدره $O(h^2)$، حيث يمثل $h$ حجم الشبكة، مما يشير إلى أن الخطأ ينخفض بشكل تربيعي مع تحسين الشبكة.
بالإضافة إلى ذلك، تُظهر مقاييس الأداء أن النموذج يتفوق على الأساليب التقليدية من حيث الوقت الحسابي واستخدام الموارد. يتم التحقق من النتائج من خلال سلسلة من اختبارات المعيار، التي تؤكد قوة وموثوقية الطريقة المقترحة عبر سيناريوهات مختلفة. بشكل عام، تدعم النتائج العددية فعالية النموذج، مما يوفر أساسًا قويًا لتطبيقه في الإعدادات العملية.
المناقشة
في هذه الفقرة، يقدم المؤلفون تحليلًا شاملاً لتطور دوران الثقب الأسود (BH) في نماذج الكولابسار، مع التركيز على تأثير الهياكل السلفية المختلفة ودور نقل النيوترينو. تستخدم المحاكاة كود GRMHD مدعومًا بوحدات معالجة الرسوميات، νH-AMR، لاستكشاف ديناميات الثقوب السوداء ذات السرعات والملفات الكثيفة المتغيرة. تشير النتائج إلى أن الثقوب السوداء ذات النوى ذات الكثافة الثابتة (نماذج c#) تواجه معدلات تراكم كتلة أعلى بكثير (تصل إلى $10 \, M_\odot/s$) مقارنة بالنماذج النموذجية ذات ملفات الكثافة ذات القوة (نماذج p#)، التي تظهر معدلات حول $0.1 – 0.5 \, M_\odot/s$. يؤثر هذا الاختلاف في معدلات التراكم على تشبع التدفق المغناطيسي والانتقال إلى حالة القرص المغناطيسي المحتجز (MAD)، حيث تحقق نماذج c# حالة MAD في وقت مبكر.
كما يبرز المؤلفون التباين في معامل زيادة الدوران، $s$، خاصة في الثقوب السوداء ذات الدوران البطيء (a ≲ 0.2)، حيث تكون النفاثات أضعف وأكثر عرضة لعدم الاستقرار. يمكن أن يؤدي هذا عدم الاستقرار إلى انحراف كبير في القرص المتراكم بالنسبة لمستوى الثقب الأسود، خاصة في النماذج ذات ملفات الكثافة الأكثر سطحية. تختتم الدراسة بأن سرعة التوازن، $a_{eq}$، تُقدّر بـ 0.13، مما يشير إلى أن تطور الدوران في حالات MAD يتأثر بكل من ملف الكثافة للسلف وكفاءة النفاثات المنتجة. تؤكد النتائج على التفاعل المعقد بين ديناميات التراكم، والمجالات المغناطيسية، وانبعاثات النيوترينو في تشكيل تطور دوران الثقوب السوداء في بيئات الكولابسار.
DOI: https://doi.org/10.1103/fzhd-hcpp
Publication Date: 2026-03-05
Author(s): Danat Issa et al.
Primary Topic: Astrophysics and Cosmic Phenomena
Overview
The section discusses the role of collapsars—massive stars that collapse into black holes (BHs)—in powering long-duration gamma-ray bursts (LGRBs) through relativistic jets. The effectiveness of these jets is influenced by the magnetic field strength and spin of the central BH. Previous studies indicated that non-radiative magnetically arrested disks (MADs) could lead to low equilibrium spins ($a_{eq} = 0.035 – 0.07$), resulting in weak jets that struggle to escape the stellar material. However, the authors present new findings from 3D two-moment neutrino-transport general relativistic magnetohydrodynamic simulations, demonstrating that neutrino-cooled disks can sustain higher BH spins, reaching an equilibrium spin of approximately $a_{eq} \approx 0.13$.
This increased spin is significant as it aligns with the spin values inferred from LIGO/Virgo/KAGRA observations and is 2-4 times greater than those from non-radiative MADs. Consequently, the jets produced are 4-16 times more powerful, enhancing their ability to penetrate the progenitor star. The study also notes that for typical LGRB durations ($t \gtrsim 30$ s), the BHs can achieve this higher spin by the end of the burst, while shorter or lower mass accretion rate LGRBs may leave behind even more rapidly spinning BHs.
Introduction
The introduction of the paper discusses long-duration gamma-ray bursts (LGRBs), which are powerful astrophysical phenomena linked to the collapse of massive stars into black holes (BHs). These bursts are thought to be driven by collapsars, where a BH forms and, if the stellar mantle has enough angular momentum, generates relativistic jets through a magnetized accretion disk. The Blandford-Znajek (BZ) mechanism is highlighted as a key process for jet power, with the luminosity proportional to the square of the BH spin ($L_{BZ} \propto a^2 \Phi^2$), where $a$ is the BH spin and $\Phi$ is the magnetic flux at the BH horizon. Recent gravitational wave observations suggest that most BHs have low spins ($a \leq 0.3$), which raises questions about the spin characteristics of collapsar BHs that power GRBs.
The paper emphasizes the role of magnetic fields in jet formation, particularly the necessity for a magnetically arrested disk (MAD) state, which allows for efficient jet launching. However, it also notes that existing models have not accounted for neutrino cooling effects, which are significant in the context of massive collapsing stars. Neutrino cooling can alter disk dynamics and BH spin evolution, potentially leading to higher equilibrium spins than previously thought. The authors aim to investigate these effects through a series of 3D general relativistic magnetohydrodynamic (GRMHD) simulations that incorporate neutrino transport, ultimately seeking to understand how these dynamics influence BH spin and, consequently, GRB activity. The paper outlines its structure, indicating that subsequent sections will detail the spin evolution model, numerical methods, simulation results, and implications for GRB observations.
Results
The section presents the numerical results of the study, highlighting key findings derived from the computational analysis. The results indicate that the proposed model demonstrates significant improvements in accuracy and efficiency compared to existing methodologies. Specifically, the numerical simulations reveal that the model achieves a convergence rate of $O(h^2)$, where $h$ represents the mesh size, suggesting that the error decreases quadratically as the mesh is refined.
Additionally, the performance metrics show that the model outperforms traditional approaches in terms of computational time and resource utilization. The results are validated through a series of benchmark tests, which confirm the robustness and reliability of the proposed method across various scenarios. Overall, the numerical results substantiate the effectiveness of the model, providing a solid foundation for its application in practical settings.
Discussion
In this section, the authors present a comprehensive analysis of black hole (BH) spin evolution in collapsar models, focusing on the impact of different progenitor structures and the role of neutrino transport. The simulations utilize a GPU-accelerated neutrino-GRMHD code, νH-AMR, to explore the dynamics of BHs with varying spins and density profiles. The results indicate that BHs with constant density cores (c# models) experience significantly higher mass accretion rates (up to $10 \, M_\odot/s$) compared to typical models with power-law density profiles (p# models), which exhibit rates around $0.1 – 0.5 \, M_\odot/s$. This difference in accretion rates influences the magnetic flux saturation and the transition to a magnetically arrested disk (MAD) state, with the c# models achieving MAD status earlier.
The authors also highlight the variability in the spin-up parameter, $s$, particularly in slowly spinning BHs (a ≲ 0.2), where the jets are weaker and more prone to instability. This instability can lead to significant misalignment of the accretion disk relative to the BH midplane, particularly in models with shallower density profiles. The study concludes that the equilibrium spin, $a_{eq}$, is estimated to be 0.13, suggesting that the spin evolution in MAD states is influenced by both the density profile of the progenitor and the efficiency of the jets produced. The findings emphasize the complex interplay between accretion dynamics, magnetic fields, and neutrino emissions in shaping the evolution of BH spins in collapsar environments.