Finite-Temperature Equations of State for Neutron-Star Mergers within the Quark Mean-Field Model

في أعماق قلب نجم نيوتروني، تُسحق المادة إلى كثافات تتحدى التجربة اليومية. هذه البقايا النجمية هي النوى المنهارة لنجوم ضخمة، وهي متراصة بإحكام شديد لدرجة أن ملعقة صغيرة واحدة من مادتها ستزن مليارات الأطنان على الأرض. لعقود من الزمن، درس العلماء هذه الأجرام من خلال مراقبة حالتها الباردة والهادئة، محاولين فهم كيف توازن الضغط الناتج عن المادة الكثيفة داخلها الجاذبية الهائلة. ومع ذلك، نادراً ما يكون الكون هادئاً. فعندما يصطدم نجمان نيوتونيان، فإنهما يخلقان حدثاً عنيفاً وفوضوويًا ترتفع فيه درجات الحرارة إلى مستويات من شأنها أن تبخر أي مادة معروفة على الأرض. في هذه اللحظات من الحرارة والضغط الشديدين، تتغير القواعد التي تحكم المادة الباردة، وتتغير الطريقة التي تقاوم بها النجوم الجاذبية بطرق بدأنا للتو في فهمها.

للتنبؤ بما يحدث أثناء هذا الاصطدام الكوني، يعتمد علماء الفلك على مجموعة من القواعد تسمى "معادلة الحالة". فكر في هذا كأنه وصفة تخبر محاكاة حاسوبية بالضبط مقدار الضغط الذي ستولده كمية معينة من المادة الساخنة والكثيفة. لفترة طويلة، استخدم العلماء نسخة مبسطة من هذه الوصفة، بافتراض أن الضغط الإضافي الناتج عن الحرارة يتصرف بطريقة بسيطة وثابتة. لكن الملاحظات الأخيرة لموجات الجاذبية — وهي تموجات في الزمكان ناتجة عن هذه الاصطدامات — تشير إلى أن الواقع أكثر تعقيداً بكفاءة. فالحرارة المتولدة أثناء الاندماج لا تضيف مجرد دفعة إضافية بسيطة؛ بل تتفاعل مع المادة بطرق معقدة تعتمد على مدى كثافة المادة ونوع الجسيمات الموجودة. إذا استخدمنا وصفة خاطئة، فإن عمليات المحاكاة لهذه الاصطدامات ستكون غير دقيقة، مما يؤدي بنا إلى سوء فهم مصير النجم المتبقي والإشارات التي يرسلها إلينا.

لقد اتخذ فريق من الباحثين خطوة كبيرة نحو إصلاح ذلك من خلال تطوير وصف أكثر تفصيلاً وواقعية للمادة الساخنة والكثيفة. لقد ركزوا على إطار نظري محدد يُعرف باسم "نموذج ميدان الكوارك المتوسط" (quark mean-field model). في هذه الرؤية، لا يتم التعامل مع البروتونات والنيوترونات التي تشكل النجم ككرات صلبة وغير قابلة للانقسام. بدلاً من ذلك، يقر النموذج بأن هذه الجسيمات تتكون في الواقع من قطع أصغر وأساسية تسمى "الكواركات". ومع زيادة الضغط، يتغير سلوك هذه الكواركات الداخلية، مما يغير بدوره طريقة استجابة البروتونات والنيوترونات لمحيطها. وقد وسع الباحثون هذا النموذج ليشمل تأثيرات درجات الحرارة العالية، مما خلق مجموعة جديدة من القواعد التي تصف كيفية سلوك هذه المادة عندما تُسخن إلى عشرات الملايين من الدرجات.

ولضمان موثوقية نموذجهم الجديد، قارن العلماء بينه وبين نهج تقليدي أكثر يتعامل مع البروتونات والنيوترونات كجسيمات بسيطة وبلا بنية. لقد قاموا بمعايرة كلا النموذجين بعناية ليتطابقا مع الخصائص المعروفة للمادة النووية عند الكثافات العادية، مما يضمن انطلاقهما من نفس الأرضية الصلبة. ثم دفعوا كلا النموذجين إلى الظروف القصوى الموجودة في اندماج النجوم النيوترونية، حيث اختبراهما عند كثافات تصل إلى ستة أضعاف كثافة النواة الذرية العادية ودرجات حرارة تتراوح من 50 مليون إلى أكثر من مليار درجة. كما فحصوا كيف تتغير طاقة المادة وضغطها وإنتروبيا (عشوائيتها) تحت هذه الظروف، مع إيلاء اهتمام وثيق لكيفية تغير الكتلة الفعالة للجسيمات مع ازدياد الازدحام والحرارة في البيئة المحيطة.

وكشفت النتائج أنه بينما اتفق كلا النموذجين على السلوك العام للمادة عند الكثافات المنخفضة، فقد بدآ في التباعد بشكل كبير مع زيادة الكثافة. في النموذج التقليدي، أصبحت الجسيمات أخف وزناً وأكثر خمولاً مع ارتفاع الكثافة. أما في النموذج الجديد القائم على الكواركات، فقد احتفظت الجسيمات بمزيد من كتلتها واستجابت للحرارة بشكل مختلف. وقد تبين أن هذا الاختلاف كان حاسماً؛ إذ وجد الباحثون أن الطريقة التي تولد بها الحرارة ضغطاً في نموذجهم الجديد ليست قيمة ثابتة وبسيطة كما يُفترض غالباً في عمليات المحاكاة. بدلاً من ذلك، فإن "المؤشر الحراري" — وهو مقياس لمدى كفاءة تحول الحرارة إلى ضغط — يتغير اعتماداً على الكثافة ودرجة الحرارة ومزيج الجسيمات المحدد. وفي بعض المناطق عالية الكثافة، توقع النموذج الجديد أن الضغط الحراري سيكون أقل مما كان يُعتقد سابقاً، مما قد يعني أن النجم المتبقي الساخن من الاندماج أقل قدرة على دعم نفسه ضد الجاذبية مما اقترحته عمليات المحاكاة السابقة.

ولجعل هذه النتائج المعقدة مفيدة للعلماء الآخرين الذين يجرون عمليات محاكاة حاسوبية ضخمة للنجوم المتصادمة، قام الفريق بترجمة حساباتهم التفصيلية إلى صيغة رياضية موجزة وسهلة الاستخدام. تسمح هذه الصيغة الجديدة للباحثين بحساب الضغط الحراري بسرعة لأي مزيج من الكثافة ودرجة الحرارة وتكوين الجسيمات ضمن نطاق الظروف المتوقعة في الاندماج. وتتميز هذه الصيغة بدقتها التي تصل إلى ضمن اثنين بالمائة من الحسابات الكاملة والمعقدة، مما يجعلها أداة عملية للجيل القادم من عمليات المحاكاة الفيزيائية الفلكية. ومن خلال استخدام هذا الوصف الأكثر واقعية للمادة الساخنة، قد تتمكن عمليات المحاكاة المستقبلية من التنبؤ باللحظة الدقيقة لانهيار النجم المندمج إلى ثقب أسود أو بقائه كجسم فائق الكتلة، وقد تفسر بشكل أفضل إشارات موجات الجاذبية المحددة التي ترصدها المراصد على الأرض. إن هذا العمل لا يكتفي فقط بتدقيق رقم ما، بل يوفر نافذة أوضح على الفيزياء العنيفة وعالية الطاقة التي تحكم أكثر الاصطدامات تطرفاً في الكون.