تتم دراسة النجوم النيوترونية المختلطة بالمادة المظلمة عبر نموذج ثنائي السائل باستخدام معادلات حالة هادرونية ومادة مظلمة غير متحيزة. تُنمذَج المادة المظلمة كغاز فيرمي يُعرَّف فقط بمعادلة حالته عند الكثافات المنخفضة وكتلته. تُرسى المادة الهادرونية عند الكثافات المنخفضة بنظرية المجال الفعال الكيرالية (chiral EFT) وعند الكثافات العالية بالكروموديناميكا الكمية الاضطرابية (pQCD)؛ ويُستخدم بارامتر لسرعة الصوت لتغطية الكثافات المتوسطة للهادرونات والكثافات العالية للمادة المظلمة، مما يضمن الاتساق الترموديناميكي دون تحيز. تشكل المادة المظلمة الخفيفة هالات ممتدة تزيد من القابلية للتشوه المدّي، بينما تشكل المادة المظلمة الثقيلة أنوية مدمجة تُقلل منها. وهكذا تنتقل القيود الرصدية من التشوه المدّي عبر الموجات الثقالية للمادة المظلمة الخفيفة إلى بيانات الكتلة ونصف القطر من NICER للمادة المظلمة الثقيلة. عند مستوى ثقة 1σ، يتم تقييد نسبة المادة المظلمة إلى f_DM ≲ 0.11 للمادة المظلمة الخفيفة. هذا الإطار غير المتحيز يُعطي حدودًا محافظة وقابلة للتطبيق على نطاق واسع، وقد تكون النجوم النيوترونية ذات الكتل المتشابهة لكن بقيم تشوه مدّي مختلفة جدًا توقيعًا دامغًا على وجود المادة المظلمة.
النجوم النيوترونية — بقايا المستعرات العظمى — من بين أكثر الأجسام كثافة في الكون. كثافتها المركزية تتجاوز الكثافة النووية بأضعاف، وجاذبيتها لا تفوقها سوى جاذبية الثقوب السوداء. ولهذا تعد النجوم النيوترونية مختبراً مثالياً لاختبار الفيزياء الأساسية: من سلوك المادة عند كثافات فوق نووية إلى خصائص المادة المظلمة. المادة المظلمة، التي تنبأت بها فيرا روبين ونشأت بعد الانفجار العظيم، تشكل ربع طاقة الكون، لكن لم تُرصد مباشرة حتى الآن. يمكن للنجوم النيوترونية أن تأسر جسيمات المادة المظلمة، مما يغير بارامتراتها المرصودة — الكتلة، ونصف القطر، وقابلية التشوه المدّي التي تقاس عبر الموجات الثقالية. البيانات الحديثة من النباضات، التي اكتشفتها جوسلين بيل بورنيل وتُدرس باستخدام التحليل الطيفي بالأشعة السينية (مهمة NICER)، ومن هوائيات الموجات الثقالية (LIGO، راينر فايس وآخرون) تتيح فرض قيود على نسبة المادة المظلمة في هذه الأجسام، دون الالتزام بنماذج محددة خارج النموذج المعياري.
لتجنب التحيز، طبق العلماء نهجاً لا أدرياً على معادلات الحالة لكل من المادة النووية والمادة المظلمة. يُرسى الجزء النووي على حسابات نظرية الحقل الفعالة الكيرالية عند كثافات منخفضة، والكروموديناميكا الكمومية الاضطرابية عند كثافات عالية، في حين تُبارمتر المنطقة الوسيطة بمقاطع عشوائية لسرعة الصوت (c²_s)، لا تقيدها إلا السببية (لا يمكن لسرعة الصوت أن تتجاوز سرعة الضوء) والاستقرار الديناميكي الحراري. بالنسبة للمادة المظلمة، استُخدمت بارامترية مماثلة، لكن دون مرساة عالية الكثافة: عند كثافات أقل من 0.1 من الكثافة النووية، وُصفت كغاز فيرمي منحل بكتلة جسيم m_D، وعند كثافات أعلى — باستيفاء سرعة الصوت. عُومل النجم كنظام ثنائي السوائل: يتفاعل المكونان عبر الجاذبية فقط، وتُحل لكل منهما معادلة توازن هيدروستاتيكي خاصة (معادلة تولمان-أوبنهايمر-فولكوف المعدلة). وُلدت مجموعات من نحو 10⁵ تسلسل للكتلة-نصف القطر لقيم مختلفة من m_D (0.2–1.1 جيجا إلكترون فولت) ونسب المادة المظلمة f_DM (0.01–0.15)، ثم قورنت بالبيانات الرصدية.
أظهر التحليل أن تأثير المادة المظلمة يعتمد جذرياً على كتلة جسيماتها. المادة المظلمة الخفيفة (m_D ~ 0.2–0.3 جيجا إلكترون فولت) تشكل هالة متخلخلة حول النجم، مما يزيد نصف قطره الفعال. يؤدي هذا إلى زيادة قابلية التشوه المدّي Λ (المتناسبة مع R⁵)، وهذه التكوينات مقيدة بشدة ببيانات حدث الموجات الثقالية GW170817. على العكس، المادة المظلمة الثقيلة (m_D > 0.7 جيجا إلكترون فولت) تتركز في نواة مدمجة، مما يصغر النجم ويخفض Λ. وهنا تكون قياسات الكتلة ونصف القطر من NICER هي المقيِّد الرئيسي. عند كتل الجسيمات في المنطقة الوسيطة (0.5–0.7 جيجا إلكترون فولت)، يكون كلا الخيارين ممكنين، إضافة إلى بنى انتقالية من نوع "نواة-هالة". أتاح التحليل الإحصائي باستخدام جميع البيانات المتاحة (NICER، GW170817، واشتراط كتلة قصوى ≥ 2.01 كتلة شمسية) تحديد أن النسبة المسموح بها من المادة المظلمة الخفيفة لا تتجاوز 11%. بالنسبة للمادة المظلمة الثقيلة، تكون القيود أشد عند النسب الكبيرة، لكن عموماً تبقى منطقة البارامترات أوسع بكثير. ومن التنبؤات المثيرة أن اكتشاف نجمين نيوترونيين بكتلتين متشابهتين، لكن بقيمتي قابلية تشوه مدّي مختلفتين بشدة، قد يكون مؤشراً قطعياً على وجود المادة المظلمة في أحدهما.
النتائج المتحصل عليها ذات طبيعة مستقلة عن النموذج، وقابلة للتطبيق على فئة واسعة من نظريات المادة المظلمة. وهي تظهر أن الأرصاد الحديثة متعددة الموجات قادرة على تقييد محتوى المادة المظلمة في النجوم النيوترونية بشكل ملموس، حتى دون معرفة طبيعتها المجهرية. إن دور القيود نفسه — قابلية التشوه المدّي أو الكتلة-نصف القطر — يعمل كأداة تشخيصية: هيمنة قيود الموجات الثقالية تشير إلى بنى هالية، بينما تشير هيمنة القيود السينية إلى تكوينات نووية للمادة المظلمة.
من المتوقع مستقبلاً أن تتيح قياسات جديدة من NICER لنجوم نابضة أخرى، إضافة إلى بيانات من كاشفات الموجات الثقالية للجيل القادم (بما في ذلك الفضائية)، تنقيح هذه الحدود بشكل كبير. قد يشمل التطوير اللاحق للطريقة اللا أدرية إدراج المادة المظلمة البوزونية ومراعاة دوران النجوم، مما يجعل النماذج أقرب إلى الواقع. كما أن دراسة اعتماد النتائج على اختيار كثافة الربط للقطاع المظلم تبدو واعدة.
ستؤثر النتائج على فيزياء المادة المظلمة، والفيزياء الفلكية للنجوم النيوترونية، وعلم فلك الموجات الثقالية، موفرةً أداة جديدة للبحث عن مرشحات الجسيمات الخفيفة للمادة المظلمة.
تشمل الخطوات التالية تحليل بيانات ذات دقة محسّنة من مهمات مستقبلية مثل STROBE-X أو Athena، وإعادة الدراسة للمادة المظلمة البوزونية.
يرتبط البحث مباشرة بالمشكلة غير المحلولة لطبيعة المادة المظلمة، وبمسألة تركيب النجوم النيوترونية عند كثافات قصوى، حيث يمكن أن تظهر مادة غريبة.
🎯 ملعقة شاي من مادة النجم النيوتروني تزن حوالي مليار طن — وهو نفس وزن جبل إيفرست تقريباً.