محاكاة عددية لاضطراب الكثافة وإنتاج موجات الجاذبية من انتقال الطور الكوني من الدرجة الأولى
Numerical simulations of density perturbation and gravitational wave production from cosmological first-order phase transition

شارك:
المجلة: Science China Physics Mechanics and Astronomy، المجلد: 69، العدد: 7
DOI: https://doi.org/10.1007/s11433-026-2971-y
تاريخ النشر: 2026-05-22
المؤلف: Jintao Zou وآخرون
الموضوع الرئيسي: علم الكون ونظريات الجاذبية

نظرة عامة

في هذه الدراسة، قمنا بإجراء محاكاة ثلاثية الأبعاد لشبكة لدراسة اضطرابات الكثافة والأمواج الجاذبية (GWs) المرتبطة بالانتقالات الطورية من الدرجة الأولى (FOPT). تشير نتائجنا إلى أنه بالنسبة لقوة الانتقال الطوري المميزة بـ $\alpha > 1$، فإن المصدر الرئيسي لاضطرابات الكثافة ينشأ من الحركة الأمامية لجدران الفقاعات. على العكس، عندما تكون $\alpha < 1$، فإن المساهمة السائدة تُنسب إلى التأخير في انهيار الفراغ. يكشف طيف الطاقة لاضطرابات الكثافة الناتجة خلال الانتقال الطوري عن ميل قدره $k^3$ عند أعداد موجية صغيرة و $k^{-1.5}$ عند أعداد موجية كبيرة. وبالمثل، تُظهر طيف الطاقة للأمواج الجاذبية المحسوبة ميلًا قدره $k^3$ عند أعداد موجية صغيرة و $k^{-2}$ عند أعداد موجية كبيرة. ومن الجدير بالذكر أن محاكياتنا العددية تشير إلى أن الانتقالات الطورية البطيئة يمكن أن تؤدي إلى تشكيل ثقوب سوداء بدائية (PBHs) وتوفر توقعات قيمة لطيف الطاقة للأمواج الجاذبية، مما يقدم دعمًا نظريًا لجهود الكشف عن الأمواج الجاذبية في المستقبل.

مقدمة

في هذا القسم، يقدم المؤلفون الثقوب السوداء البدائية (PBHs)، التي يُعتقد أنها تتشكل في الكون المبكر بسبب اضطرابات كثافة كبيرة. يتم مناقشة آليات مختلفة لتشكيل PBH، بما في ذلك الانتقالات الطورية من الدرجة الأولى (FOPT)، والانزلاق البطيء للانفلاتون، وانهيار الخيوط الكونية، بالإضافة إلى العمليات غير الاضطرابية التي تحفز الاضطرابات القياسية. تُعتبر PBHs بارزة لقدرتها على تغطية نطاق واسع من الكتلة، من المقاييس دون الكوكبية إلى المقاييس العملاقة، وتعتبر مرشحة للكتلة المظلمة ومصادر للأمواج الجاذبية القابلة للكشف من قبل المراصد مثل LIGO-Virgo والبعثات المستقبلية مثل LISA.

يؤكد المؤلفون على دور FOPT في توليد تقلبات كثافة كبيرة من خلال نواة وتوسع فقاعات الفراغ، والتي تساهم أيضًا في خلفية الأمواج الجاذبية العشوائية. تهدف الدراسة إلى محاكاة عددية لاضطرابات الكثافة الناشئة من FOPT، مع فحص كيفية تأثير معدلات وقوة الانتقالات الطورية المختلفة على كل من اضطرابات الكثافة وطيف الأمواج الجاذبية. تشير النتائج الأولية إلى أن عشوائية النفق الكمي خلال الانتقالات الطورية يمكن أن تنتج اضطرابات انحناء أكبر بكثير من تلك الناتجة عن التضخم، مع تفضيل القيم الأصغر من المعامل $\beta/H$ لتشكيل PBH. ومن الجدير بالذكر أن المحاكيات تشير إلى أن حالة مع $\beta/H = 6$ تؤدي إلى انهيار ناجح لـ PBH.

النتائج

في هذا القسم، يقدم المؤلفون نتائج عددية من محاكيات نموذج انتقال طوري من الدرجة الأولى (FOPT) في شبكة بحجم \( L^3 = (512 \, dx)^3 \) مع إجمالي \( 4^3 \) من أحجام هابل. تركز الدراسة على تحديد مناطق انهيار الفراغ الكاذب المتأخر، والتي تتميز بتغيرات كبيرة في كثافة الطاقة المحلية بسبب النفق الحراري. تم تصنيف آخر 16 حجم هابل شهدت نواة الفقاعات على أنها مؤجلة، مما يشير إلى مناطق تبقى في حالة الفراغ الكاذب لفترة أطول بسبب معدل نواة منخفض. تشمل المعلمات المستكشفة \( \alpha = 0.5, 1, 5, 10 \) و \( \beta/H = 6, 8, 10, 12 \)، مع اختيار معلمات محتملة محددة لتسهيل تحليل مستقل عن النموذج لاضطرابات الكثافة وتشكيل الثقوب السوداء البدائية (PBH).

تشير النتائج إلى أن وفرة PBHs تزداد مع \( \alpha \) وتتناقص مع \( \beta/H \)، مع اعتماد أقوى على \( \beta/H \). ومن الجدير بالذكر أن تشكيل PBH يُفضل خلال الانتقالات الطورية البطيئة مع \( \beta/H < 6 \). كما تحسب الدراسة اضطراب الكثافة \( \sigma_\delta \) عند مقياس هابل، كاشفة أن قيمته القصوى تتأثر بالتحولات الطاقية خلال الانتقال الطوري. يظهر طيف الطاقة لاضطرابات الكثافة ميولًا مميزة في نطاقات مختلفة، مع ميل نهائي قدره \( k^{-1.5} \) يُنسب إلى تأثيرات التنعيم للطاقة الحركية وطاقة التدرج. بالإضافة إلى ذلك، يتم تحليل أطياف الأمواج الجاذبية (GW)، مما يظهر حساسية للمعلمات \( \alpha \) و \( \beta/H \)، مع إمكانية الكشف عنها من قبل مراصد مختلفة مثل LIGO وSKA.

المناقشة

في هذا القسم، يستكشف المؤلفون ديناميات الانتقالات الطورية (PTs) التي تشمل حقلًا قياسيًا $\phi$ وتشكيل الثقوب السوداء البدائية (PBHs) من خلال محاكيات عددية. يتم تعريف إمكانات الحقل القياسي على أنها \( V(\phi) = \frac{1}{2} M^2 \phi^2 + \frac{1}{3} \kappa \phi^3 + \frac{1}{4} \lambda \phi^4 \). تحدث نواة فقاعات الفراغ بمعدل أسي، وتتحكم معادلة كلاين-غوردون في تطور الحقل في كون متوسع. تجد الدراسة أن المناطق التي تعاني من تأخير في انهيار الفراغ يمكن أن تطور كثافات محلية زائدة، وهو أمر حاسم لتشكيل PBH. يُظهر وفرة PBHs أنها تزداد مع قوة معامل الانتقال الطوري $\alpha$ ولكن تتناقص مع معامل المعدل $\beta/H$. ومن الجدير بالذكر أن حساسية وفرة PBH للتغيرات في $\beta/H$ أقوى من تلك المتعلقة بـ $\alpha$، خاصة عندما تكون $\alpha \gtrsim 10$.

كما يحلل المؤلفون طيف الطاقة لاضطرابات الكثافة، كاشفين عن قمتين متميزتين: واحدة بسبب انهيار الفراغ المتأخر والأخرى من حركة جدران الفقاعات. تقترح الدراسة أنه بالنسبة لـ $\alpha < 1$، تهيمن القمة الأولى، بينما بالنسبة لـ $\alpha > 1$، تصبح القمة الثانية أكثر أهمية. علاوة على ذلك، يتم فحص طيف الطاقة للأمواج الجاذبية (GW) الناتجة خلال FOPTs، مما يظهر أنها تعتمد بشكل ضعيف على معامل حالة المادة $\omega$. تشير النتائج إلى أن الانتقالات الطورية الأسرع (أكبر $\beta/H$) تقمع تشكيل PBH بسبب توزيع الطاقة الأكثر تجانسًا، بينما يبقى طيف GW متسقًا عبر ظروف مختلفة. بشكل عام، توفر هذه الدراسة رؤى حول التفاعل المعقد بين الانتقالات الطورية، واضطرابات الكثافة، وتشكيل PBHs في الكون المبكر.

Journal: Science China Physics Mechanics and Astronomy, Volume: 69, Issue: 7
DOI: https://doi.org/10.1007/s11433-026-2971-y
Publication Date: 2026-05-22
Author(s): Jintao Zou et al.
Primary Topic: Cosmology and Gravitation Theories

Overview

In this study, we performed three-dimensional lattice simulations to investigate density perturbations and gravitational waves (GWs) associated with first-order phase transitions (FOPT). Our findings indicate that for phase transition strengths characterized by $\alpha > 1$, the primary source of density perturbations arises from the forward motion of bubble walls. Conversely, when $\alpha < 1$, the predominant contribution is attributed to the delay in vacuum decay. The power spectrum of density perturbations generated during the phase transition reveals a slope of $k^3$ at small wavenumbers and $k^{-1.5}$ at large wavenumbers. Similarly, the calculated GW power spectra show a slope of $k^3$ at small wavenumbers and $k^{-2}$ at large wavenumbers. Notably, our numerical simulations suggest that slow phase transitions can lead to the formation of primordial black holes (PBHs) and provide valuable predictions for the GW power spectrum, thereby offering theoretical support for future GW detection efforts.

Introduction

In this section, the authors introduce primordial black holes (PBHs), which are theorized to form in the early universe due to significant density perturbations. Various mechanisms for PBH formation are discussed, including first-order phase transitions (FOPT), inflaton ultraslow roll, and the collapse of cosmic strings, as well as non-perturbative processes that induce scalar perturbations. PBHs are notable for their potential to span a wide mass range, from subplanetary to supermassive scales, and are considered candidates for dark matter and sources of gravitational waves detectable by observatories like LIGO-Virgo and future missions such as LISA.

The authors emphasize the role of FOPT in generating large density fluctuations through the nucleation and expansion of vacuum bubbles, which also contribute to a stochastic gravitational wave background. The study aims to numerically simulate density perturbations arising from FOPT, examining how different rates and strengths of phase transitions affect both density perturbations and gravitational wave spectra. Preliminary findings suggest that quantum tunneling randomness during phase transitions can produce curvature perturbations significantly larger than those from inflation, with smaller values of the parameter $\beta/H$ favoring PBH formation. Notably, simulations indicate that a case with $\beta/H = 6$ leads to successful PBH collapse.

Results

In this section, the authors present numerical results from simulations modeling a first-order phase transition (FOPT) in a lattice of size \( L^3 = (512 \, dx)^3 \) with a total of \( 4^3 \) Hubble volumes. The study focuses on identifying regions of delayed false vacuum decay, which are characterized by significant local energy density variations due to thermal tunneling. The last 16 Hubble volumes to undergo bubble nucleation are classified as postponed, indicating areas that remain in the false vacuum state longer due to a low nucleation rate. The parameters explored include \( \alpha = 0.5, 1, 5, 10 \) and \( \beta/H = 6, 8, 10, 12 \), with specific potential parameters chosen to facilitate a model-independent analysis of density perturbations and primordial black hole (PBH) formation.

The results indicate that the abundance of PBHs increases with \( \alpha \) and decreases with \( \beta/H \), with a stronger dependence on \( \beta/H \). Notably, PBH formation is favored during slow phase transitions with \( \beta/H < 6 \). The study also computes the density perturbation \( \sigma_\delta \) at the Hubble scale, revealing that its maximum value is influenced by the energy transformations during the phase transition. The power spectrum of density perturbations exhibits distinct slopes in different ranges, with a final slope of \( k^{-1.5} \) attributed to the smoothing effects of kinetic and gradient energies. Additionally, gravitational wave (GW) spectra are analyzed, showing a sensitivity to the parameters \( \alpha \) and \( \beta/H \), with potential detectability by various observatories such as LIGO and SKA.

Discussion

In this section, the authors investigate the dynamics of phase transitions (PTs) involving a scalar field $\phi$ and the formation of primordial black holes (PBHs) through numerical simulations. The potential for the scalar field is defined as \( V(\phi) = \frac{1}{2} M^2 \phi^2 + \frac{1}{3} \kappa \phi^3 + \frac{1}{4} \lambda \phi^4 \). The nucleation of vacuum bubbles occurs at an exponential rate, and the evolution of the field is governed by the Klein-Gordon equation in an expanding universe. The study finds that regions with delayed vacuum decay can develop local overdensities, which are crucial for PBH formation. The abundance of PBHs is shown to increase with the strength of the phase transition parameter $\alpha$ but decrease with the rate parameter $\beta/H$. Notably, the sensitivity of PBH abundance to variations in $\beta/H$ is stronger than that to $\alpha$, particularly when $\alpha \gtrsim 10$.

The authors also analyze the power spectrum of density perturbations, revealing two distinct peaks: one due to delayed vacuum decay and another from the motion of bubble walls. The study suggests that for $\alpha < 1$, the first peak dominates, while for $\alpha > 1$, the second peak becomes more significant. Furthermore, the gravitational wave (GW) power spectrum generated during FOPTs is examined, showing that it is weakly dependent on the equation of state parameter $\omega$. The findings indicate that faster phase transitions (larger $\beta/H$) suppress PBH formation due to a more homogeneous energy distribution, while the GW spectrum remains consistent across different conditions. Overall, this work provides insights into the complex interplay between phase transitions, density perturbations, and the formation of PBHs in the early universe.

شارك: