DOI: https://doi.org/10.1103/x9fk-ykm1
تاريخ النشر: 2026-01-20
المؤلف: Zhenyun Du وآخرون
الموضوع الرئيسي: الدراسات النجمية والكوكبية والمجرية
نظرة عامة
تقدم هذه الدراسة مجموعة شاملة من جداول معادلات الحالة (EOS) لمادة النجوم الكثيفة الباردة، مع التركيز بشكل خاص على كسور البروتون \( Y_p = 0.01 – 0.65 \) وكثافات عدد الباريونات \( n_b = 10^{-8} – 2 \, \text{fm}^{-3} \). تستخدم الدراسة 13 دالة كثافة نسبية مختلفة (RDFs) لاشتقاق هذه المعادلات، والتي تلبي شرط الاستقرار β للنجوم النيوترونية بكثافات \( n_b = 7.6 \times 10^{-11} – 2 \, \text{fm}^{-3} \). تشير النتائج إلى أنه عند الكثافات المنخفضة (\( n_b \lesssim 0.015 \, \text{fm}^{-3} \))، تظهر مادة النجوم بشكل رئيسي في مرحلة القطرات، بينما تظهر هياكل أكثر تعقيدًا مثل القضبان والألواح والأنابيب والفقاعات مع زيادة الكثافة. تم تحديد كسور البروتون الحرجة لتسرب النيوترونات لتكون تقريبًا \( Y_{\text{drip}}^p \approx 0.26 – 0.31 \)، مع عدم حساسية المعادلات لاختيار RDFs عند الكثافات الصغيرة أو كسور البروتون العالية.
تستكشف الدراسة أيضًا العلاقة بين الكثافة وكسور البروتون والهياكل المجهرية الناتجة، كاشفة أن كثافات الانتقال بين المراحل المختلفة تتأثر بكل من \( Y_p \) وRDFs المختارة. من الجدير بالذكر أن الضغط \( P \) والطاقة لكل باريون \( E/A \) تزداد مع الكثافة، مع تغير غير أحادي لـ \( E/A \) بالنسبة لـ \( Y_p \). تسلط الدراسة الضوء على عدم اليقين الكبير في المعادلات والهياكل، خاصة حول بداية تسرب النيوترونات وعند الكثافات الأعلى (\( n_b \gtrsim 0.1 \, \text{fm}^{-3} \))، حيث تنتقل المادة إلى حالة موحدة وتظهر الميونات. تمتد تداعيات هذه النتائج إلى فهم اهتزازات g-mode في النجوم النيوترونية، مما يشير إلى أن الاختلافات في سرعات الصوت قد توفر رؤى حول ميل طاقة التماثل \( L \) من خلال الدراسات الرصدية المستقبلية.
مقدمة
تؤكد مقدمة الورقة على الدور الحاسم لمعادلة الحالة (EOS) في فهم سلوك مادة النجوم الكثيفة، وهو أمر محوري لمجموعة متنوعة من الأحداث الفلكية مثل انفجارات السوبرنوفا الناتجة عن انهيار النواة، وتطور النجوم النيوترونية، واندماجات النجوم النيوترونية الثنائية. يبرز المؤلفون التحديات التي تطرحها الكروموديناميكا الكمومية الشبكية (QCD) في تحديد المعادلة بدقة عبر نطاقات الكثافة ودرجة الحرارة الواسعة ذات الصلة بهذه الظواهر. على الرغم من هذه التحديات، فإن خصائص المادة النووية عند كثافة التشبع ($n_0 \approx 0.16 \, \text{fm}^{-3}$) تعتبر راسخة نسبيًا، مع معلمات رئيسية تشمل طاقة ربط قدرها $B \approx -16 \, \text{MeV}$، وعدم الانضغاطية $K = 240 \pm 20 \, \text{MeV}$، وطاقة التماثل $S = 31.7 \pm 3.2 \, \text{MeV}$، وميلها $L = 58.7 \pm 28.1 \, \text{MeV}$.
من المتوقع أن تؤدي الجهود التجريبية الأخيرة، وخاصة قياسات PREX-II وCREX لسمك قشرة النيوترون، إلى تحسين ميل طاقة التماثل. ومع ذلك، فإن هذه القياسات تظهر حاليًا تباينات كبيرة، مما يستلزم مزيدًا من التحقيق. بالإضافة إلى ذلك، يشير المؤلفون إلى أنه يمكن اشتقاق قيود إضافية على المادة النووية عند الكثافات الأعلى (حتى حوالي $4.5n_0$) من تصادمات الأيونات الثقيلة النسبية، والتي توفر رؤى حول معلمات مثل معامل الانحراف $J = -390 \pm 60$.
النتائج
في هذا القسم، يستكشف المؤلفون خصائص وهياكل مادة النجوم الكثيفة المميزة بكثافات باريونية محددة ($n_b$) وكسور بروتون ($Y_p$). يوضحون الظهور المتسلسل لمراحل مختلفة—قطرة، قضيب، لوح، أنبوب، وفقاعة—داخل خلايا ويجنر-سايتز عند كسر بروتون ثابت قدره $Y_p = 0.3$، كما هو موضح في الشكل 1. تشير ملفات الكثافة إلى أنه في مراحل القطرات والقضبان، تكون كثافات البروتونات والنيوترونات أعلى في المركز ($r = 0$) مقارنة بالحافة ($r = R_W$)، بينما العكس صحيح بالنسبة لمراحل الأنابيب والفقاعات. يتم تحديد حجم القطرات ($R_d$) باستخدام معادلة محددة، مع تمثيل النتائج بصريًا بواسطة خطوط برتقالية متقطعة عمودية في الشكل.
تستخدم الدراسة 13 دالة كثافة نسبية مختلفة (RDFs)، بما في ذلك كل من الاقترانات الذاتية غير الخطية والاقترانات المعتمدة على الكثافة، لتحليل طاقة الربط لكل نوكليون ($\epsilon(n_b, Y_p)$) لمادة نووية موحدة. يتم التعبير عن طاقة الربط كدالة لكثافة التشبع النووي ($n_0$) ومعلمات مختلفة مثل طاقة الربط ($B$)، وعدم الانضغاطية ($K$)، وطاقة التماثل ($S$). يقدم المؤلفون تحليلًا مفصلًا لنطاقات وأحجام الخطوات للمعلمات المستخدمة في معادلاتهم، كما هو ملخص في الجدول I، ويبلغون أن الخصائص المشبعة التي تتنبأ بها RDFs تتماشى مع نطاقات عدم اليقين الحالية في الأدبيات.
مناقشة
في هذا القسم، تركز المناقشة على خصائص ومعادلات الحالة (EOS) لمادة النجوم الكثيفة الباردة، خاصة في سياق النجوم النيوترونية وتشكيلها. تسلط الدراسة الضوء على أهمية علم الفلك متعدد الرسائل في تقييد خصائص مادة النجوم الكثيفة، باستخدام قياسات من النباضات وأحداث موجات الجاذبية. من الجدير بالذكر أن الدراسة تؤكد على أهمية دوال الكثافة النسبية المختلفة (RDFs) في نمذجة المعادلة، مع استخدام 13 دالة مختلفة لالتقاط عدم اليقين في دوال كثافة الطاقة النووية. تشير النتائج إلى أن الضغط والطاقة لكل باريون لمادة النجوم الكثيفة تزداد مع كثافة عدد الباريونات ($n_b$)، مع ظهور مراحل هيكلية متميزة عند كثافات وكسور بروتون مختلفة ($Y_p$).
يكشف التحليل أنه عند الكثافات المنخفضة ($n_b \lesssim 10^{-3} \text{fm}^{-3}$)، يتأثر الضغط بشكل رئيسي بغاز الإلكترونات المتدهورة، بينما عند الكثافات الأعلى ($n_b \gtrsim 0.2 \text{fm}^{-3}$)، تهيمن المادة النووية. تحدد الدراسة أيضًا انتقالًا من الدرجة الأولى من مرحلة السائل إلى الغاز عند الكثافات دون التشبع، مما يؤدي إلى تشكيل هياكل غير متجانسة، مثل المعكرونة النووية، التي تلعب دورًا حاسمًا في خصائص النقل والمرونة للنجوم النيوترونية. تختتم الدراسة بأن المعادلة والهياكل المجهرية لمادة النجوم الكثيفة ضرورية لفهم ظواهر فلكية متنوعة، بما في ذلك انفجارات السوبرنوفا الناتجة عن انهيار النواة واندماجات النجوم النيوترونية، وأن الملاحظات المستقبلية ستعمل على تحسين هذه النماذج بشكل أكبر.
DOI: https://doi.org/10.1103/x9fk-ykm1
Publication Date: 2026-01-20
Author(s): Zhenyun Du et al.
Primary Topic: Stellar, planetary, and galactic studies
Overview
This research presents a comprehensive set of equations of state (EOS) tables for cold dense stellar matter, specifically focusing on proton fractions \( Y_p = 0.01 – 0.65 \) and baryon number densities \( n_b = 10^{-8} – 2 \, \text{fm}^{-3} \). The study employs 13 different relativistic density functionals (RDFs) to derive these EOSs, which fulfill the β-stability condition for neutron stars with densities \( n_b = 7.6 \times 10^{-11} – 2 \, \text{fm}^{-3} \). The findings indicate that at lower densities (\( n_b \lesssim 0.015 \, \text{fm}^{-3} \)), the stellar matter predominantly exhibits a droplet phase, while more complex structures such as rods, slabs, tubes, and bubbles emerge as density increases. The critical proton fractions for neutron drip are identified to be approximately \( Y_{\text{drip}}^p \approx 0.26 – 0.31 \), with the EOSs showing insensitivity to the choice of RDFs at small densities or high proton fractions.
The study further explores the relationship between density, proton fraction, and the resulting microscopic structures, revealing that the transition densities between different phases are influenced by both \( Y_p \) and the selected RDFs. Notably, the pressure \( P \) and energy per baryon \( E/A \) increase with density, with a nonmonotonic variation of \( E/A \) concerning \( Y_p \). The research highlights significant uncertainties in the EOSs and structures, particularly around the onset of neutron drip and at higher densities (\( n_b \gtrsim 0.1 \, \text{fm}^{-3} \)), where the matter transitions to a uniform state and muons appear. The implications of these findings extend to the understanding of g-mode oscillations in neutron stars, suggesting that the differences in sound velocities could provide insights into the symmetry energy slope \( L \) through future observational studies.
Introduction
The introduction of the paper emphasizes the critical role of the equation of state (EOS) in understanding the behavior of dense stellar matter, which is pivotal for various astronomical events such as core-collapse supernovae, neutron star evolution, and binary neutron star mergers. The authors highlight the challenges posed by lattice Quantum Chromodynamics (QCD) in accurately determining the EOS across the extensive density and temperature ranges relevant to these phenomena. Despite these challenges, the properties of nuclear matter at saturation density ($n_0 \approx 0.16 \, \text{fm}^{-3}$) are relatively well established, with key parameters including a binding energy of $B \approx -16 \, \text{MeV}$, incompressibility $K = 240 \pm 20 \, \text{MeV}$, symmetry energy $S = 31.7 \pm 3.2 \, \text{MeV}$, and its slope $L = 58.7 \pm 28.1 \, \text{MeV}$.
Recent experimental efforts, particularly the PREX-II and CREX measurements of neutron skin thicknesses, are anticipated to refine the slope of the symmetry energy. However, these measurements currently exhibit significant discrepancies, necessitating further investigation. Additionally, the authors note that further constraints on nuclear matter at higher densities (up to approximately $4.5n_0$) can be derived from relativistic heavy-ion collisions, which provide insights into parameters such as the skewness coefficient $J = -390 \pm 60$.
Results
In this section, the authors investigate the properties and microscopic structures of dense stellar matter characterized by specific baryon densities ($n_b$) and proton fractions ($Y_p$). They illustrate the sequential emergence of various phases—droplet, rod, slab, tube, and bubble—within Wigner-Seitz cells at a fixed proton fraction of $Y_p = 0.3$, as shown in Figure 1. The density profiles indicate that in the droplet and rod phases, proton and neutron densities are higher at the center ($r = 0$) compared to the edge ($r = R_W$), while the opposite is true for the tube and bubble phases. The droplet size ($R_d$) is determined using a specified equation, with the results visually represented by vertical orange dotted lines in the figure.
The study employs 13 different relativistic density functionals (RDFs), including both nonlinear self-couplings and density-dependent couplings, to analyze the binding energy per nucleon ($\epsilon(n_b, Y_p)$) of uniform nuclear matter. The binding energy is expressed as a function of the nuclear saturation density ($n_0$) and various parameters such as binding energy ($B$), incompressibility ($K$), and symmetry energy ($S$). The authors provide a detailed breakdown of the ranges and step sizes for the parameters used in their equations, as summarized in Table I, and report that the saturation properties predicted by the RDFs align with current uncertainty bands in the literature.
Discussion
In this section, the discussion focuses on the properties and equations of state (EOS) of cold dense stellar matter, particularly in the context of neutron stars and their formation. The research highlights the significance of multimessenger astronomy in constraining the characteristics of dense stellar matter, utilizing measurements from pulsars and gravitational wave events. Notably, the study emphasizes the importance of various relativistic density functionals (RDFs) in modeling the EOS, with 13 different RDFs employed to capture the uncertainties in nuclear energy density functionals. The findings indicate that the pressure and energy per baryon of dense stellar matter increase with baryon number density ($n_b$), with distinct structural phases emerging at varying densities and proton fractions ($Y_p$).
The analysis reveals that at low densities ($n_b \lesssim 10^{-3} \text{fm}^{-3}$), the pressure is primarily influenced by the degenerate electron gas, while at higher densities ($n_b \gtrsim 0.2 \text{fm}^{-3}$), nuclear matter dominates. The study also identifies a first-order liquid-gas phase transition at subsaturation densities, leading to the formation of nonuniform structures, such as nuclear pasta, which play a critical role in the transport and elastic properties of neutron stars. The research concludes that the EOS and microscopic structures of dense stellar matter are crucial for understanding various astrophysical phenomena, including core-collapse supernovae and neutron star mergers, and that future observations will further refine these models.