DOI: https://doi.org/10.3847/1538-4357/ae731d
تاريخ النشر: 2026-06-18
المؤلف: Meng-Yuan Ho وآخرون
الموضوع الرئيسي: علم الفلك ودراسات تكوين النجوم
نظرة عامة
في هذه الدراسة، نقدم محاكاة عالية الدقة لتشكيل أول سحب تشكيل النجوم داخل 15 هالة صغيرة، بأحجام تتراوح من حوالي $10^5$ إلى $10^7 \, M_\odot$ عند انزياحات حمراء $z \sim 17-20$. باستخدام كود GIZMO، تتضمن محاكاة لدينا فيزياء الغاز البدائي بالتفصيل وتستفيد من الظروف الأولية من محاكاة إيلستريسTNG الكونية. من خلال استخدام تقنية تقسيم الجسيمات، نقوم بتحسين الدقة بمقدار حوالي $10^5$، مما يحقق كتل جسيمات الغاز والمادة المظلمة بمقدار $0.2 \, M_\odot$ و $80 \, M_\odot$، على التوالي. هذه الدقة العالية تسمح لنا بحل عمليات تراكم الغاز بفعالية ومراقبة ظهور اضطرابات قوية فوق الصوتية، تتميز بأرقام ماخ تتراوح من 1.8 إلى 4.2، والتي تزداد مع كتلة الهالة. تؤدي الاضطرابات إلى تجزئة سحابة الغاز المركزية إلى كتل كثيفة متعددة، بأحجام تتراوح بين $2.6 \, M_\odot$ و $66.5 \, M_\odot$، متجاوزة كتل جانز المقابلة لها وانهار لتشكيل النجوم الأولى.
تشير نتائجنا إلى أن الاضطرابات فوق الصوتية هي سمة شائعة في الهالات الصغيرة وتلعب دورًا حاسمًا في تشكيل سحب تشكيل النجوم المتكتلة، مما له آثار كبيرة على دالة الكتلة الأولية للنجوم الأولى. تكشف المحاكاة أن سرعات الغاز داخل الهالات الصغيرة يمكن أن تصل إلى مستويات فوق صوتية، مع أرقام ماخ النموذجية تتراوح بين 2 و 5، ترتبط إيجابيًا مع كتلة الهالة. تعطل الاضطرابات الانهيار الكروي للغاز، مما يؤدي إلى هياكل خيطية متكتلة، متوافقة مع التنبؤات السابقة. توفر الخصائص الفيزيائية للاضطرابات التي تم تحديدها في محاكاة لدينا رؤى أساسية حول تشكيل نجوم المجموعة الثالثة وآليات ردود الفعل الخاصة بها، والتي تعتبر حيوية في تشكيل خصائص المجرات الأولى. هذه الآثار ذات صلة خاصة بالملاحظات المستقبلية للمجرات المبكرة باستخدام تلسكوب جيمس ويب الفضائي (JWST).
مقدمة
تناقش مقدمة الورقة تشكيل وخصائص نجوم المجموعة الثالثة (Pop III)، التي نشأت من الغاز البدائي المكون أساسًا من الهيدروجين والهيليوم في هالات المادة المظلمة بأحجام تتراوح بين $10^5$ و $10^6 \, M_\odot$. تشكلت هذه النجوم المبكرة، التي تفتقر إلى تبريد خطوط المعادن، بأحجام جانز أعلى بكثير، مما أدى إلى كتل نجمية تتراوح بين 40 إلى 500 $M_\odot$، على عكس الكتل الأقل التي تُلاحظ عادة في الكون الحالي. تشير المحاكاة الحديثة إلى توزيع كتلة أوسع لنجوم Pop III، مع نتائج تقترح دوال الكتلة الأولية (IMFs) تمتد من حوالي $10^{-3}$ إلى $100 \, M_\odot$، متأثرة بعوامل مثل الاضطرابات وآليات ردود الفعل.
تهدف الدراسة إلى التحقيق في الخصائص الفيزيائية لسحب تشكيل نجوم Pop III داخل الهالات الصغيرة، مع معالجة التباينات بين التنبؤات النظرية والبيانات الرصدية، لا سيما فيما يتعلق بكتلة النجوم من الجيل الثاني التي تم إثراؤها بواسطة انفجارات نجوم Pop III. يبرز المؤلفون قيود المحاكاة السابقة في التقاط الاضطرابات أثناء تجميع الهالات الصغيرة ويقترحون توسيع النتائج الحديثة التي توصل إليها تشين وآخرون (2025) حول الاضطرابات فوق الصوتية في الهالات الصغيرة. من خلال استخدام مجموعة من الهالات الصغيرة من محاكاة إيلستريسTNG وتعزيز الدقة من خلال تقنية تقسيم الجسيمات، تسعى الدراسة إلى وصف خصائص الغاز المضطرب وآثارها على الظروف الأولية لتشكيل نجوم Pop III. تم هيكلة الورقة لتفصيل الأساليب، والنتائج حول تشكيل الاضطرابات، والآثار الفلكية لهذه النتائج.
طرق
في هذا القسم، يوضح المؤلفون منهجية المحاكاة الخاصة بهم، والتي تستند إلى الظروف الأولية المأخوذة من مشروع إيلستريسTNG، وهو محاكاة كونية شاملة تقوم بنمذجة تطور الهيكل والمادة في الكون. يركزون بشكل خاص على 15 هالة صغيرة بأحجام تتراوح من $10^5$ إلى $10^7 \, M_\odot$، تم تحديدها من لقطة إيلستريسTNG عند انزياح أحمر $z \approx 20$، والتي تمثل بيئات ملائمة لتشكيل نجوم المجموعة الثالثة.
لزيادة دقة هذه الهالات الصغيرة، يستخدم المؤلفون كود GIZMO في وضع الكتلة النهائية بدون شبكة (MFM)، والذي يتميز بقدرته على محاكاة ديناميات الغاز في بيئات فلكية معقدة. نظرًا لأن الدقة الأصلية لإيلستريسTNG غير كافية لالتقاط ديناميات الغاز على النطاقات الصغيرة، فإنهم ينفذون خوارزمية تقسيم الجسيمات التي تزيد من الدقة بمقدار $10^5$. يسمح هذا التحسين بفحص أكثر تفصيلًا لعمليات تراكم الغاز المضطرب والتجزئة أثناء تجميع الهالات الصغيرة. يستخدم GIZMO، الذي يبني على إطار GADGET-2، خوارزمية شجرة الجسيمات-الشبكة (TreePM) للحسابات الجاذبية، مما يسهل تقييمات القوة الدقيقة والفعالة من خلال التحلل الهرمي لحجم المحاكاة.
مناقشة
إطار محاكاة GIZMO، الذي يدمج نهج MFM/MFV، يعزز الديناميكا الهيدروليكية التقليدية المعتمدة على الجسيمات (SPH) من خلال استخدام صيغة موزعة الحجم تجمع بين مزايا SPH والأساليب الإيولرية المعتمدة على الشبكة. يحسن هذا النهج الهجين بشكل كبير من دقة المحاكاة في التقاط الظواهر المعقدة مثل الصدمات والتدفقات المضطربة مع ضمان خصائص الحفاظ القوية والثبات العددي. تم التحقق من صحة الإطار عبر سيناريوهات فلكية متنوعة، مما يؤكد فعاليته في المحاكاة عالية الدقة. يتم نمذجة عمليات التبريد والعمليات الكيميائية للغاز البدائي باستخدام مكتبة GRACKLE، التي تتعقب شبكة من اثني عشر نوعًا بدائيًا وتدمج آليات تبريد مختلفة، مما يمكّن من فهم تفصيلي للتطور الحراري والكيميائي للغاز البدائي.
في سياق مشروع إيلستريسTNG، تركز الدراسة على الهالات الصغيرة، باستخدام مجموعات بيانات عالية الدقة لتحديد مرشحي الهالة بناءً على معايير محددة للكتلة والعزل. يتم استخدام تقنية تقسيم الجسيمات لتعزيز دقة الكتلة، مما يسمح بحل فعال لهياكل تدفق الغاز أثناء تشكيل الهالات الصغيرة. تكشف المحاكاة عن ديناميات غاز معقدة تتميز باضطرابات فوق صوتية وهياكل خيطية مدفوعة بالتفاعلات الجاذبية. تظهر طيف طاقة الحركة للغاز انحدارًا على شكل قانون القوة يتماشى مع طيف كولموغوروف، مما يشير إلى وجود اضطرابات عبر مقاييس مختلفة. توضح التحليلات وجود ارتباط إيجابي بين رقم ماخ المضطرب وكتلة الهالة، مما يشير إلى أن الهالات الصغيرة الأكثر ضخامة تولد تدفقات مضطربة أقوى بسبب آفاق جاذبية أعمق. بشكل عام، توفر النتائج رؤى حول تشكيل وتطور الغاز البدائي في الهالات الصغيرة، مما يبرز التفاعل المعقد بين الديناميات الجاذبية وفيزياء الغاز في الكون المبكر.
DOI: https://doi.org/10.3847/1538-4357/ae731d
Publication Date: 2026-06-18
Author(s): Meng-Yuan Ho et al.
Primary Topic: Astrophysics and Star Formation Studies
Overview
In this study, we present high-resolution simulations of the formation of the first star-forming clouds within 15 minihalos, with masses ranging from approximately $10^5$ to $10^7 \, M_\odot$ at redshifts $z \sim 17-20$. Utilizing the GIZMO code, our simulations incorporate detailed primordial gas physics and leverage initial conditions from the IllustrisTNG cosmological simulations. By employing a particle-splitting technique, we enhance the resolution by a factor of about $10^5$, achieving gas and dark matter particle masses of $0.2 \, M_\odot$ and $80 \, M_\odot$, respectively. This high resolution allows us to effectively resolve gas accretion processes and observe the emergence of strong supersonic turbulence, characterized by Mach numbers ranging from 1.8 to 4.2, which increases with halo mass. The turbulence leads to the fragmentation of the central gas cloud into multiple dense clumps, with masses between $2.6 \, M_\odot$ and $66.5 \, M_\odot$, surpassing their corresponding Jeans masses and collapsing to form the first stars.
Our findings indicate that supersonic turbulence is a prevalent feature in minihalos and plays a crucial role in the formation of clumpy star-forming clouds, which has significant implications for the initial mass function of the first stars. The simulations reveal that gas velocities within the minihalos can reach supersonic levels, with characteristic Mach numbers between 2 and 5, positively correlating with halo mass. The turbulence disrupts the spherical collapse of gas, resulting in clumpy filamentary structures, consistent with previous predictions. The physical properties of the turbulence identified in our simulations provide essential insights into the formation of Population III stars and their feedback mechanisms, which are pivotal in shaping the characteristics of the first galaxies. These implications are particularly relevant for future observations of early galaxies using the James Webb Space Telescope (JWST).
Introduction
The introduction of the paper discusses the formation and characteristics of Population III (Pop III) stars, which originated from primordial gas primarily composed of hydrogen and helium in dark matter halos with masses between $10^5$ and $10^6 \, M_\odot$. These early stars, lacking metal-line cooling, formed with significantly higher Jeans masses, resulting in stellar masses ranging from 40 to 500 $M_\odot$, as opposed to the lower masses typically observed in the current universe. Recent simulations indicate a broader mass distribution for Pop III stars, with findings suggesting initial mass functions (IMFs) spanning from approximately $10^{-3}$ to $100 \, M_\odot$, influenced by factors such as turbulence and feedback processes.
The study aims to investigate the physical properties of Pop III star-forming clouds within minihalos, addressing the discrepancies between theoretical predictions and observational data, particularly regarding the mass of second-generation stars enriched by Pop III supernovae. The authors highlight the limitations of previous simulations in capturing turbulence during minihalo assembly and propose to extend recent findings by Chen et al. (2025) on supersonic turbulence in minihalos. By utilizing a suite of minihalos from the IllustrisTNG simulation and enhancing resolution through a particle-splitting technique, the study seeks to characterize the turbulent gas properties and their implications for the initial conditions of Pop III star formation. The paper is structured to detail the methods, findings on turbulence formation, and the astrophysical implications of these results.
Methods
In this section, the authors detail their simulation methodology, which is based on initial conditions sourced from the IllustrisTNG project, a comprehensive cosmological simulation that models the evolution of structure and matter in the universe. They specifically focus on 15 minihalos with masses ranging from $10^5$ to $10^7 \, M_\odot$, identified from an IllustrisTNG snapshot at redshift $z \approx 20$, which are representative of environments conducive to Population III star formation.
To enhance the resolution of these minihalos, the authors utilize the GIZMO code in meshless finite mass (MFM) mode, which is adept at simulating gas dynamics in complex astrophysical settings. Given that the original resolution of IllustrisTNG is inadequate for capturing small-scale gas dynamics, they implement a particle-splitting algorithm that increases the resolution by a factor of $10^5$. This refinement allows for a more detailed examination of turbulent gas accretion and fragmentation processes during minihalo assembly. GIZMO, which builds on the GADGET-2 framework, employs a Tree Particle-Mesh (TreePM) algorithm for gravitational calculations, facilitating precise and efficient force evaluations through hierarchical decomposition of the simulation volume.
Discussion
The GIZMO simulation framework, which integrates the MFM/MFV approach, enhances traditional smoothed-particle hydrodynamics (SPH) by employing a volume-discretized formulation that merges the advantages of SPH and grid-based Eulerian methods. This hybrid approach significantly improves the accuracy of simulations in capturing complex phenomena such as shocks and turbulent flows while ensuring robust conservation properties and numerical stability. The framework has been validated across various astrophysical scenarios, confirming its efficacy for high-resolution simulations. The cooling and chemical processes of primordial gas are modeled using the GRACKLE library, which tracks a network of twelve primordial species and incorporates various cooling mechanisms, thereby enabling a detailed understanding of the thermal and chemical evolution of primordial gas.
In the context of the IllustrisTNG project, the study focuses on minihalos, utilizing high-resolution datasets to identify halo candidates based on specific mass and isolation criteria. The particle-splitting technique is employed to enhance mass resolution, allowing for the effective resolution of gas flow structures during minihalo formation. The simulations reveal complex gas dynamics characterized by supersonic turbulence and filamentary structures driven by gravitational interactions. The kinetic energy power spectra of the gas exhibit a power-law decay consistent with the Kolmogorov spectrum, indicating the presence of turbulence across various scales. The analysis demonstrates a positive correlation between the turbulent Mach number and halo mass, suggesting that more massive minihalos generate stronger turbulent flows due to deeper gravitational potentials. Overall, the findings provide insights into the formation and evolution of primordial gas in minihalos, highlighting the intricate interplay between gravitational dynamics and gas physics in the early universe.
