DOI: https://doi.org/10.1088/1475-7516/2026/04/049
تاريخ النشر: 2026-04-01
المؤلف: Luis E. Padilla وآخرون
الموضوع الرئيسي: علم الكون ونظريات الجاذبية
نظرة عامة
في هذه الدراسة، نستكشف الانهيار الحرج الذي يؤدي إلى تشكيل الثقوب السوداء البدائية (PBH) داخل كون يتميز بحقل عددي ذاتي التفاعل مع إمكانات ربعوية، مقارنين إياه بالسيناريو التقليدي المهيمن على الإشعاع. باستخدام محاكاة عددية غير خطية كاملة النسبية في تناظر كروي، مستندة إلى صيغة ميسنر-شارب، نقوم بتحليل الديناميات التي تقترب من عتبة الانهيار ونراقب مقياس كتلة الثقب الأسود. تكشف نتائجنا أن كل من حالة الحقل العددي وحالة الإشعاع تظهر سلوكًا حرجًا من النوع الثاني، على الرغم من اختلافات طفيفة في الأسس الحرجة، تختلف بحوالي $2\sigma$.
هذا يشير إلى أنه بينما يمكن لحقل عددي ربعوي أن يحاكي بشكل فعال سائل الإشعاع في نظام الانهيار غير الخطي، تبقى هناك اختلافات طفيفة في الأساس الحرج. تدعم نتائجنا شبه عمومية الأساس الحرج في تشكيل الثقوب السوداء البدائية، مما يظهر اعتمادًا خفيفًا فقط على ما إذا كان الانهيار ناتجًا عن حقل عددي أو سائل مثالي. تسهم هذه الدراسة في فهم أعمق للشروط التي تتشكل فيها الثقوب السوداء البدائية والديناميات الأساسية لمثل هذه العمليات في سياقات كونية متنوعة.
مقدمة
تناقش مقدمة الورقة ظاهرة الانهيار الجاذبي الحرج، الذي يتوازى مع الانتقالات الطورية في الفيزياء الإحصائية. تسلط الضوء على وجود معلمة حرجة \( p_{\text{th}} \) تميز بين نتيجتين في الأنظمة المنهارة: بالنسبة لـ \( p < p_{\text{th}} \)، يبقى الزمكان منتظمًا، بينما بالنسبة لـ \( p > p_{\text{th}} \)، يتشكل ثقب أسود. تظهر كتلة الثقب الأسود الناتج سلوكًا من نوع القوة \( M \propto (p – p_{\text{th}})^\gamma \)، حيث \( \gamma \) هو أساس حرج. تصنف الورقة الانهيار الحرج إلى النوع الأول والنوع الثاني، اعتمادًا على سلوك الكتلة وطبيعة الحل الحرج، مع تحديد النوع بناءً على وجود مقاييس داخلية في معادلات الحقل.
يهدف المؤلفون إلى التحقيق في صحة تشبيه السائل المثالي في سياق تشكيل الثقوب السوداء البدائية (PBH)، وخاصة بالنسبة للحقول العددية ذات الإمكانات الرباعية. أشارت الأعمال السابقة إلى وجود اختلافات في سلوك الحقول العددية مقارنة بالسوائل المثالية، خاصة في الحالة التربيعية، حيث انحرفت ديناميات الانهيار عن سلوك الغبار. تستخدم هذه الدراسة محاكاة عددية لمقارنة الانهيار الحرج لحقل عددي ذاتي التفاعل ربعوي مع سائل إشعاعي، حيث وجدت أن كلاهما يظهر سلوكًا حرجًا مشابهًا وأسس حرجة متطابقة. تدعم هذه النتيجة استخدام نماذج السوائل المثالية في سيناريوهات الحقول العددية، بينما تؤكد أيضًا على الحاجة إلى التحقق الدقيق من مثل هذه التشبيهات، نظرًا لتداعياتها على طيف كتلة الثقوب السوداء البدائية وديناميات الانهيار.
طرق
في هذا القسم، يصف المؤلفون المنهجية العددية المستخدمة في محاكاتهم، والتي تعتمد على كود تم تأسيسه مسبقًا. يتم تنفيذ تطور الزمن للمتغيرات الديناميكية ($\psi$, $\chi$, $\Pi$, $R$, $\rho_{pf}$, $\tilde{U}$, $\lambda$) باستخدام مُدمج رانج-كوتا من الدرجة الرابعة على طول إحداثي زمني لوغاريتمي $\xi$. يتم حساب المشتقات المكانية باستخدام طرق الفرق المحدودة من الدرجة الرابعة، وتفرض شروط حدود عاكسة عند الأصل لضمان اختفاء المشتقات الأولى. يستخدم المؤلفون أيضًا قالب فرق محدود خارجي من سبع نقاط للمشتقات العددية عند حدود الاستئصال ويطبقون شروط حدود نيومان عند الحدود الخارجية.
لتحسين الاستقرار العددي وتقليل التذبذبات الزائفة، يتم تضمين مصطلح تشتت من الدرجة الثالثة لكريس-أوليجر. يتم تهيئة المحاكاة بعامل مقياس $a(t_0) = 1$ وزمن لوغاريتمي $\xi_0 = 0$، مع افتراض معادلة حالة سائل خلفية فعالة $w = 1/3$. يتم نمذجة الملف الأولي للزيادة الكثيفة كدالة غاوسية، مع كون المعلمة $p$ هي المتغير الحر الوحيد. يتناول المؤلفون الاعتماد الدائري في حساب الشرط الأولي لـ $\lambda$ من خلال استخدام إجراء تكراري للتقارب على قيمة متسقة لـ $\Gamma$. يتم مراقبة تشكيل الثقوب السوداء البدائية (PBHs) من خلال التحقق من شروط محددة، خاصة $\theta^+ = 0$، ويتم تتبع تطور الأفق الظاهر لتحديد متى تستقر كتلة PBH.
نقاش
في هذا القسم، يقدم المؤلفون صياغة رياضية مفصلة لدراسة تطور الاضطرابات في كل من الأكوان المهيمنة على الإشعاع والحقول العددية باستخدام صيغة ميسنر-شارب. يسمح هذا النهج بتحليل الأنظمة الكروية المتماثلة، مع تضمين سائل مثالي وحقل عددي. يستخرج المؤلفون مجموعة من المعادلات التفاضلية التي تحكم ديناميات هذه الحقول، مع التأكيد على أهمية الشروط الأولية وشروط الحدود لضمان الانتظام عند الأصل. يقدمون متغيرات بلا أبعاد لتعزيز الاستقرار والدقة العددية، مما يحول المعادلات إلى شكل أكثر قابلية للإدارة.
يناقش المؤلفون أيضًا المنهجية لتأسيس الشروط الأولية، خاصة للسيناريوهات المهيمنة على السوائل المثالية والحقول العددية. يسلطون الضوء على أهمية التوسعات الاضطرابية والعلاقة بين الاضطرابات الأولية وملفات الانحناء. علاوة على ذلك، يستكشفون الشروط اللازمة لتشكيل الثقوب السوداء البدائية (PBH)، مع تحديد الأسطح المحاصرة من خلال توسيع الجيوديسيات الصفرية. يتم تقديم دالة الانضغاط كأداة تشخيصية لتقييم الانهيار الجاذبي للاضطرابات، مع الإشارة إلى أن أقصى حد لها يدل على إمكانية تشكيل PBH. يختتم القسم بمناقشة السلوك الحرج الذي لوحظ في كل من الحالات المهيمنة على الإشعاع والحقول العددية، مؤكدًا على قوة ظواهر الانهيار الحرج وشبه عمومية الأساس الحرج عبر سيناريوهات مختلفة.
DOI: https://doi.org/10.1088/1475-7516/2026/04/049
Publication Date: 2026-04-01
Author(s): Luis E. Padilla et al.
Primary Topic: Cosmology and Gravitation Theories
Overview
In this study, we explore the critical collapse leading to primordial black hole (PBH) formation within a universe characterized by a self-interacting scalar field with a quartic potential, contrasting it with the conventional radiation-dominated scenario. Utilizing fully relativistic nonlinear numerical simulations in spherical symmetry, grounded in the Misner-Sharp formalism, we analyze the dynamics approaching the collapse threshold and monitor the scaling of the black hole mass. Our findings reveal that both the scalar field and radiation cases exhibit type II critical behavior, albeit with slightly different critical exponents, differing by approximately $2\sigma$.
This indicates that while a quartic scalar field can effectively simulate a radiation fluid in the nonlinear collapse regime, minor discrepancies in the critical exponent remain. Our results substantiate the near universality of the critical exponent in PBH formation, demonstrating only a mild dependence on whether the collapse is instigated by a scalar field or a perfect fluid. This work contributes to a deeper understanding of the conditions under which PBHs form and the underlying dynamics of such processes in varying cosmological contexts.
Introduction
The introduction of the paper discusses the phenomenon of critical gravitational collapse, which parallels phase transitions in statistical physics. It highlights the existence of a critical parameter \( p_{\text{th}} \) that distinguishes between two outcomes in collapsing systems: for \( p < p_{\text{th}} \), spacetime remains regular, while for \( p > p_{\text{th}} \), a black hole forms. The mass of the resulting black hole exhibits power-law scaling \( M \propto (p – p_{\text{th}})^\gamma \), where \( \gamma \) is a critical exponent. The paper categorizes critical collapse into type I and type II, depending on the mass behavior and the nature of the critical solution, with the type determined by the presence of intrinsic scales in the field equations.
The authors aim to investigate the validity of the perfect fluid analogy in the context of primordial black hole (PBH) formation, particularly for scalar fields with quartic potentials. Previous work indicated discrepancies in the behavior of scalar fields compared to perfect fluids, especially in the quadratic case, where the collapse dynamics deviated from dust-like behavior. This study employs numerical simulations to compare the critical collapse of a quartic self-interacting scalar field with that of a radiation fluid, finding that both exhibit similar critical behavior and matching critical exponents. This result supports the use of perfect fluid models in scalar field scenarios, while also emphasizing the need for careful verification of such analogies, given their implications for the PBH mass spectrum and collapse dynamics.
Methods
In this section, the authors describe the numerical methodology employed for their simulations, which builds upon a previously established code. The time evolution of the dynamical variables ($\psi$, $\chi$, $\Pi$, $R$, $\rho_{pf}$, $\tilde{U}$, $\lambda$) is executed using a fourth-order Runge-Kutta integrator along a logarithmic time coordinate $\xi$. Spatial derivatives are calculated with fourth-order finite difference methods, and reflective boundary conditions are imposed at the origin to ensure the first derivatives vanish. The authors also utilize a seven-point outward finite difference stencil for numerical derivatives at the excision boundary and apply Neumann boundary conditions at the outer boundary.
To enhance numerical stability and mitigate spurious oscillations, a third-order Kreiss-Oliger dissipation term is incorporated. The simulations are initialized with a scale factor $a(t_0) = 1$ and a logarithmic time $\xi_0 = 0$, assuming an effective background fluid equation of state $w = 1/3$. The initial profile for the overdensity is modeled as a Gaussian function, with the parameter $p$ being the only free variable. The authors address the circular dependency in calculating the initial condition for $\lambda$ by employing an iterative procedure to converge on a consistent value of $\Gamma$. The formation of primordial black holes (PBHs) is monitored by checking specific conditions, particularly $\theta^+ = 0$, and the evolution of the apparent horizon is tracked to determine when the PBH mass stabilizes.
Discussion
In this section, the authors present a detailed mathematical formulation to study the evolution of perturbations in both radiation- and scalar field-dominated universes using the Misner-Sharp formalism. This approach allows for the analysis of spherically symmetric systems, incorporating a perfect fluid and a scalar field. The authors derive a set of differential equations governing the dynamics of these fields, emphasizing the importance of initial conditions and boundary conditions to ensure regularity at the origin. They introduce dimensionless variables to enhance numerical stability and accuracy, transforming the equations into a more manageable form.
The authors also discuss the methodology for establishing initial conditions, particularly for scenarios dominated by perfect fluids and scalar fields. They highlight the significance of perturbative expansions and the relationship between initial perturbations and curvature profiles. Furthermore, they explore the conditions for primordial black hole (PBH) formation, identifying trapped surfaces through the expansion of null geodesics. The compaction function is introduced as a diagnostic tool to assess the gravitational collapse of perturbations, with its maximum indicating the potential for PBH formation. The section concludes by addressing the critical behavior observed in both radiation- and scalar field-dominated cases, confirming the robustness of critical collapse phenomena and the near universality of the critical exponent across different scenarios.