DOI: https://doi.org/10.1103/z5xb-fxdf
تاريخ النشر: 2026-01-02
المؤلف: Krinio Marouda وآخرون
الموضوع الرئيسي: الثقوب السوداء والفيزياء النظرية
نظرة عامة
في هذه الدراسة، نستكشف عتبة الانهيار لحقل قياسي مركب بلا كتلة في الزمان والمكان المحوريين، معتمدين على الإطار الذي وضعه شوبتويك وآخرون وتركزنا على المعامل الزاوي \( m \). من خلال تنفيذ تقنية جديدة لتقليل التماثل تُسمى “m-cartoon” واستخدام كود بامبس، نقوم بتحليل بيانات أولية قريبة من الحرجة لوحدات الزاوية \( m = 1 \) و\( m = 2 \). تكشف نتائجنا أنه بالنسبة لـ \( m = 1 \)، يظهر النظام تماثلًا ذاتيًا متقطع (DSS) مع فترة صدى \( \Delta \approx 0.42 \) ومؤشر قياس القوة \( \gamma \approx 0.11 \). بالمقابل، بالنسبة لـ \( m = 2 \)، نلاحظ حلاً مميزًا لـ DSS يتميز بمعاملات حرجة أصغر، \( \Delta \approx 0.09 \) و\( \gamma \approx 0.035 \)، مما يشير إلى أنه بينما تستمر العالمية ضمن عائلات \( m \) الثابتة، فإن المؤشرات الحرجة تعتمد على وضع الزاوية.
تظهر التحليلات الإضافية للعلاقة بين الزخم الزاوي \( J_{AH} \) والكتلة \( M_{AH} \) عند لحظة تشكيل الأفق الظاهر أن الدوران غير البعدي \( \chi_{AH} = \frac{J_{AH}}{M_{AH}^2} \) يقترب من الصفر، مما يشير إلى أن الزخم الزاوي له تأثير ضئيل عند عتبة الانهيار. يتم تأكيد ذلك من خلال مؤشر ليابونو السلبي الفرعي \( \lambda_1 \)، والذي يشير إلى أن الزخم الزاوي يتناقص بسرعة أكبر من الكتلة مع اقتراب العتبة. كما تسلط نتائجنا الضوء على أنه، على عكس النتائج السابقة في السيناريوهات غير المتشابكة، لا يوجد تنافس بين عتبات القياس والفراغ، وتحافظ الحلول الحرجة على خصائص عالمية على الرغم من التغيرات في الزخم الزاوي. تؤكد هذه الدراسة على أهمية وحدات الزاوية في ظواهر الانهيار الحرجة وتقترح طرقًا للبحث المستقبلي في حلول الفراغ مع التواء.
مقدمة
تتناول مقدمة هذه الورقة البحثية تطور فهم الانهيار الجاذبي في النسبية العامة (GR) منذ العمل الرائد لشوبتويك، الذي حدد عتبة حرجة بين الانهيار والتشتت في الزمان والمكان المميزين بحقول قياسية بلا كتلة. تحدد هذه العتبة الزمان والمكان تحت الحرجة، التي تتشتت، من تلك فوق الحرجة التي تنهار إلى ثقوب سوداء. تسلط الورقة الضوء على القياس العالمي للقوة الملاحظ في كلا النظامين، مع كتلة الثقوب السوداء على الجانب فوق الحرجة وثوابت الانحناء على الجانب تحت الحرجة، مما يشير إلى نظام حرجي حيث يوجد تماثل ذاتي متقطع (DSS). يشير المؤلفون إلى أن هذه الظواهر قوية عبر عائلات بيانات أولية مختلفة وصيغ من GR، مما يجعلها أرض اختبار قيمة لرموز النسبية العددية.
توسع الدراسة العمل السابق من خلال فحص حالة الحقول القياسية المركبة مع الزخم الزاوي في إعدادات محورية، باستخدام فرضية محددة للحقل القياسي. يؤكد المؤلفون أن إدخال الزخم الزاوي يعقد الديناميات، حيث لا يوجد حد كروي في هذا السياق. يهدفون إلى استكشاف السلوك الحرجي في هذا النظام غير الكروي، مع التركيز بشكل خاص على إمكانية الثقوب السوداء كير القصوى عند عتبة الانهيار. توضح الورقة الأساليب العددية المستخدمة، بما في ذلك التعديلات لتقليل التماثل واكتشاف الأفق، وتضع الأساس لفحص مفصل للعتبة لتشكيل الثقوب السوداء كير، وهو جانب مهم من التحقيق المستمر في الظواهر الحرجة في الانهيار الجاذبي.
طرق
في هذا القسم، يصف المؤلفون الصياغة والأساليب العددية المستخدمة لتنفيذ نظام في التماثل المحوري باستخدام كود بامبس، وهو أداة طيفية مصممة للتطور الزمني للأنظمة الهيدروستاتيكية المتناظرة من الدرجة الأولى. يقدمون طريقة m-cartoon، وهي تقنية جديدة لتقليل التماثل تعمم الطريقة التقليدية لتناسب الزمان والمكان المميزين بالزخم الزاوي الصافي في الحقول القياسية المركبة.
يقدم المؤلفون نظرة عامة موجزة عن المعادلات المستمرة المعالجة، وتكوين الشبكة، وحل بيانات البداية، واستراتيجية تحسين الشبكة hp، بالإضافة إلى تحسينات على معلمات تخفيف القيود. يشيرون إلى وثائق فنية إضافية للحصول على تفاصيل شاملة حول هذه المكونات. علاوة على ذلك، ينتهي القسم بمناقشة التعديلات التي تم إجراؤها على مكتشف الأفق الظاهر لدمج الالتواء، مما يبرز التقدم المنهجي في نهجهم.
نتائج
في هذا القسم، يقدم المؤلفون نتائج عددية من دراستهم، بدءًا من الحالة التي يكون فيها المعامل \( m = 1 \) ثم تحليل الحالة التي يكون فيها \( m = 2 \). تشمل النتائج فحصًا مفصلًا للزخم الزاوي ومحتوى الموجات الجاذبية، مع التركيز بشكل خاص على سلوك هذه الكميات بالقرب من عتبة الظواهر المدروسة. يسمح هذا النهج المنظم بفهم شامل للديناميات المعنية في كلا السيناريوهين، مع تسليط الضوء على الاختلافات والتشابهات الرئيسية في تطور النظام مع تغير المعامل \( m \).
مناقشة
في هذا القسم، يستكشف المؤلفون ديناميات حقل قياسي مركب بلا كتلة مع الزخم الزاوي ضمن إطار النسبية العامة (GR). يتم صياغة العمل باستخدام وحدات هندسية، مما يؤدي إلى معادلات حقل أينشتاين المرتبطة بمعادلة كلاين-غوردون بلا كتلة. تستخدم الدراسة فرضية محددة تقدم وضعًا زاويًا واحدًا للحقل القياسي، مما يحافظ على التماثل المحوري للزمان والمكان بينما يسمح بمكونات مقياس غير صفرية تساهم في الزخم الزاوي كومار الصافي. يتم تقديم متجه الالتواء 4 لوصف فشل المتجه القاتل المرتبط بالتماثل الزاوي ليكون عموديًا على السطح، مما يشير إلى وجود درجات حرية دورانية في الزمان والمكان.
يستخدم المؤلفون مقياس هارموني عام (GHG) لصياغتهم لـ GR، مما يضمن نظامًا هيدروستاتيكيًا متناظرًا من الدرجة الأولى. يقومون بتنفيذ نظام تخفيف للقيود الهارمونية ويستخدمون اختيار مقياس محدد يتكيف مع التماثل خلال مرحلة الاستقرار الديناميكي للتطور. يتم تقليل المعادلات التي تحكم الحقل القياسي المركب إلى نظام من الدرجة الأولى، مع تصميم بيانات البداية لضمان الانتظام في المركز، وهو أمر حاسم لدقة أساليبهم العددية الطيفية. تتضمن بيانات البداية الزخم الزاوي وتم بناؤها لتقليل الانتهاكات المتبقية للقيود الهاميلتونية والزخم. كما يوضح المؤلفون إعدادهم الحاسوبي، بما في ذلك استراتيجية تحسين الشبكة التكيفية للتعامل مع الدقة العالية المطلوبة في منطقة الحقل القوي، ويقدمون طريقة “m-cartoon” لتطوير الزمان والمكان المحوري بكفاءة مع الأخذ في الاعتبار الاعتماد الزاوي للحقل القياسي.
بشكل عام، يوفر هذا العمل إطارًا شاملاً لدراسة التفاعل بين الحقول القياسية والزخم الزاوي في سياق نسبي عام، مع آثار لفهم تشكيل الثقوب السوداء والديناميات في سيناريوهات الانهيار الجاذبي.
DOI: https://doi.org/10.1103/z5xb-fxdf
Publication Date: 2026-01-02
Author(s): Krinio Marouda et al.
Primary Topic: Black Holes and Theoretical Physics
Overview
In this study, we explore the collapse threshold of a massless complex scalar field in axisymmetric spacetimes, employing the framework established by Choptuik et al. and focusing on the azimuthal parameter \( m \). By implementing a new symmetry reduction technique termed “m-cartoon” and utilizing the pseudospectral code bamps, we analyze near-critical initial data for angular modes \( m = 1 \) and \( m = 2 \). Our findings reveal that for \( m = 1 \), the system exhibits discrete self-similarity (DSS) with an echoing period \( \Delta \approx 0.42 \) and a power-law scaling exponent \( \gamma \approx 0.11 \). In contrast, for \( m = 2 \), we observe a distinct DSS solution characterized by smaller critical parameters, \( \Delta \approx 0.09 \) and \( \gamma \approx 0.035 \), indicating that while universality persists within fixed \( m \) families, the critical exponents are dependent on the angular mode.
Further analysis of the relationship between angular momentum \( J_{AH} \) and mass \( M_{AH} \) at the moment of apparent horizon formation shows that the dimensionless spin \( \chi_{AH} = \frac{J_{AH}}{M_{AH}^2} \) approaches zero, suggesting that angular momentum has a negligible effect at the threshold of collapse. This is corroborated by a negative subdominant Lyapunov exponent \( \lambda_1 \), which indicates that angular momentum diminishes more rapidly than mass as the threshold is approached. Our results also highlight that, unlike previous findings in non-twisting scenarios, there is no competition between scalar and vacuum thresholds, and the critical solutions maintain universal characteristics despite variations in angular momentum. This study underscores the significance of angular modes in critical collapse phenomena and suggests avenues for future research into vacuum solutions with twist.
Introduction
The introduction of this research paper discusses the evolution of understanding gravitational collapse in general relativity (GR) since Choptuik’s seminal work, which identified a critical threshold between collapse and dispersion in spacetimes characterized by massless scalar fields. This threshold delineates subcritical spacetimes, which disperse, from supercritical ones that collapse into black holes. The paper highlights the universal power-law scaling observed in both regimes, with the mass of black holes on the supercritical side and curvature invariants on the subcritical side, indicating a critical regime where discrete self-similarity (DSS) is present. The authors note that these phenomena are robust across various initial data families and formulations of GR, making them a valuable testing ground for numerical relativity codes.
The study extends previous work by examining the case of complex scalar fields with angular momentum in axisymmetric settings, utilizing a specific ansatz for the scalar field. The authors emphasize that introducing angular momentum complicates the dynamics, as no spherical limit exists in this context. They aim to explore critical behavior in this non-spherical regime, particularly focusing on the potential for extremal Kerr black holes at the threshold of collapse. The paper outlines the numerical methods employed, including adaptations for symmetry reduction and horizon detection, and sets the stage for a detailed examination of the threshold for Kerr black hole formation, which is a significant aspect of the ongoing investigation into critical phenomena in gravitational collapse.
Methods
In this section, the authors describe the formulation and numerical methods employed to implement a system in axisymmetry using the bamps code, a pseudospectral tool designed for the time evolution of first-order symmetric hyperbolic systems. They introduce the m-cartoon method, a novel symmetry reduction technique that generalizes the traditional cartoon method to accommodate spacetimes characterized by net angular momentum in complex scalar fields.
The authors provide a concise overview of the continuum equations addressed, the grid configuration, the initial data solver, and the hp-mesh-refinement strategy, along with enhancements to the constraint damping parameters. They reference additional technical documentation for comprehensive details on these components. Furthermore, the section concludes with a discussion of modifications made to the apparent horizon finder to incorporate twist, highlighting the methodological advancements in their approach.
Results
In this section, the authors present numerical results from their study, beginning with the case where the parameter \( m = 1 \) and subsequently analyzing the case where \( m = 2 \). The findings include a detailed examination of the angular momentum and the gravitational wave content, particularly focusing on the behavior of these quantities near the threshold of the studied phenomena. This structured approach allows for a comprehensive understanding of the dynamics involved in both scenarios, highlighting key differences and similarities in the evolution of the system as the parameter \( m \) varies.
Discussion
In this section, the authors explore the dynamics of a massless complex scalar field with angular momentum within the framework of General Relativity (GR). The action is formulated using geometric units, leading to Einstein’s field equations coupled with the massless Klein-Gordon equation. The study employs a specific ansatz that introduces a single angular mode to the scalar field, preserving the axisymmetry of the spacetime while allowing for non-zero metric components that contribute to the net Komar angular momentum. The twist 4-vector is introduced to characterize the failure of the Killing vector associated with azimuthal symmetry to be hypersurface orthogonal, indicating the presence of rotational degrees of freedom in the spacetime.
The authors utilize a generalized harmonic gauge (GHG) for their formulation of GR, ensuring a first-order symmetric hyperbolic system. They implement a damping scheme for harmonic constraints and employ a specific gauge choice that adapts to the symmetry during the dynamical stability phase of the evolution. The equations governing the complex scalar field are reduced to a first-order system, with the initial data designed to ensure regularity at the center, which is crucial for the accuracy of their pseudospectral numerical methods. The initial data incorporates angular momentum and is constructed to minimize residual violations of the Hamiltonian and momentum constraints. The authors also detail their computational setup, including an adaptive mesh refinement strategy to handle the high resolution required in the strong-field region, and introduce the “m-cartoon” method to efficiently evolve the axisymmetric spacetime while accounting for the azimuthal dependence of the scalar field.
Overall, this work provides a comprehensive framework for studying the interplay between scalar fields and angular momentum in a general relativistic context, with implications for understanding black hole formation and dynamics in gravitational collapse scenarios.