Controlling artificial surface heating in neutron star simulations: Application to hybrid equations of state

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

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

في دراسة حديثة، بحث باحث في جامعة غوته فرانكفورت في كيفية إيقاف هذا التسخين الاصطناعي. ركز العمل على نوع معين من الأكواد الحاسوبية المستخدمة لمحاكاة النجوم النيوترونية، والمعروف باسم كود "BAM". يُستخدم هذا الكود على نطاق واسع لنمذجة الاصطدامات العنيفة للنجوم النيوترونية الثنائية، حيث يلتف نجمان نحو بعضهما البعض ويندمجان في النهاية. اختبر الباحث استراتيجية رياضية جديدة تسمى "تحديد التدفق القائم على الإنتروبيا" (entropy-based flux-limiting). وببساطة، تعمل هذه الطريقة كمرشح ذكي لحسابات الحاسوب؛ فهي تراقب تدفق الطاقة والمادة عند سطح النجم وتعدل طريقة الحساب تلقائيًا عندما تكتشف أنواع التغيرات الحادة التي تسبب الأخطاء عادةً. ومن خلال القيام بذلك، تمنع الحاسوب من توليد تلك الحرارة الزائفة في المقام الأول.

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

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

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