Investigating Twin Star Equation of States in Light of Recent Astrophysical Observations

ملخص تقني: استقصاء معادلات الحالة للنجوم التوأم في ضوء الأرصاد الفلكية الحديثة

بيان المشكلة
لا يزال التركيب الداخلي للنجوم النيوترونية (NSs) وطبيعة الانتقال الطوري (PT) من المادة الهادرونية (HM) إلى المادة الكواركية (QM) أمراً غير محل حل. وبينما تتنبأ الديناميكا اللونية الكمومية (QCD) بانتقال عند كثافات عالية، فإن الظروف المحددة — مثل كثافة البداية، وضغط الانتقال، وقوة الانقطاع — تظل مجهولة. يمكن أن يؤدي الانتقال الطوري من الدرجة الأولى القوي إلى وجود "النجوم التوأم" (TSs): وهي فرع ثالث من الأجرام المتراصة المستقرة في مخطط الكتلة-النصف قطر ($M-R$)، مفصول عن فرع النجم النيوتروني القياسي بمنطقة غير مستقرة. ستكون هذه النجوم التوأم ذات كتلة متساوية ولكن بنصف قطر أصغر من نظيراتها من النجوم النيوترونية. ويكمن التحدي في تحديد معادلات الحالة الهجينة (EoS) التي تستوفي كل من التوازن الهيدروستاتيكي النسبي العام والقيود المتزايدة الدقة من الفيزياء الفلكية متعددة الرسائل، بما في ذلك قياسات نصف قطر الكتلة من NICER وأرصاد موجات الجاذبية (GW170817).

المنهجية
يتبنى المؤلفون نهجاً غير مرتبط بنموذج محدد لبناء مجموعة شاملة من معادلات الحالة الهجينة القادرة على إنتاج حلول النجوم التوأم. تعتمد المنهجية على معادلات تولمان-أوبنهايمر-فولكوف (TOV) لحل التوازن الهيدروستاتيكي. يتضمن بناء معادلة الحالة أربعة مناطق متميزة:

1. منطقة القشرة: معادلة حالة "بوليتروبية" مجزأة وثابتة تحت $\sim 0.4\rho_{sat}$.
2. المنطقة الهادرونية: تُنمذج باستخدام "بوليتروب" بمعامل "أدياباتيكي" $\gamma$ حر، يتم تحديده بواسطة الشروط الحدية عند الانتقال.
3. منطقة الانتقال الطوري: تُنمذج عبر إنشاء "ماكسويل"، وتتميز بضغط انتقال ثابت ($P_{tr}$) وقفزة منفصلة في كثافة الطاقة ($\Delta e$).
4. المنطقة الكواركية: تُنمذج باستخدام بارامتريّة سرعة الصوت الثابتة ($c^2_{s,QM}$).

يتم تعريف فضاء المعلمات بواسطة أربعة متغيرات مستقلة: كثافة طاقة الانتقال ($e_{tr}$)، ضغط الانتقال ($P_{tr}$)، قوة انقطاع كثافة الطاقة ($\Delta e$)، وسرعة الصوت في المادة الكواركية ($c^2_{s,QM}$). تقوم الدراسة بتغيير هذه المعلمات بشكل منهجي ضمن حدود السببية لتحديد الظروف التي تنتج حلاً للفرع الثالث في تسلسل الكتلة-النصف قطر.

المساهمات الرئيسية

• تعريف منحنيات "قبعة الساحرة" (Witch-Hat Curves): يقدم البحث تمثيلاً هندسياً في فضاء المعلمات $(P_{tr}, \Delta e)$ لـ $e_{tr}$ و $c^2_{s,QM}$ ثابتين. تحدد هذه المنحنيات حدود قوة الانقطاع المسموح بها التي تسمح بتكوين النجوم التوأم. يحدد المؤلفون ثلاث مناطق متميزة داخل هذا الفضاء:
  • المنطقة الخالية من الحلول: لا تتكون نجوم توأم بغض النظر عن $\Delta e$.
  • المنطقة المحدودة: تتكون النجوم التوأم لنطاق محدد من $\Delta e$ تحت ضغط انتقال أقصى.
  • المنطقة غير المحدودة: تتكون النجوم التوأم لـ $\Delta e$ فوق حد أدنى، دون وجود حد أقصى لـ $\Delta e$ (ضمن الإعداد العددي).
• تحليل السببية والاكتمال: يكشف البحث أن منحنيات "قبعة الساحرة" يمكن أن تكون "مثقوبة" أو غير مكتملة. يحدث هذا عندما يتجاوز أقصى ضغط انتقال مسموح به ($P_{tr}^{max}$) مطلوب لتكوين نجم توأم حد السببية ($c_s = 1$) أو الضغط المركزي للتكوين الهادروني ذي الكتلة القصوى.
• دمج القيود الرصدية الحديثة: يطبق المؤلفون فترات الثقة $2\sigma$ (وكذلك $1\sigma/3\sigma$) من قياسات NICER الحديثة (PSR J0614–3329, PSR J0437–4715, PSR J0740+6620, PSR J0030+0451)، والنبض النابض الضخم PSR J0348+0432، والتشوه المدّي من GW170817 لتقييد فضاء المعلمات.

النتائج
يؤدي تطبيق البيانات الرصدية إلى تقييد فضاء المعلمات القابل للتحقق بشكل كبير:

• القيود على كثافة الانتقال: تشير القياسات الحديثة لـ PSR J0614–3329 و PSR J0437–4715 إلى معادلة حالة أكثر ليونة حول الكتلة النموذجية، مما يقيد بشدة كثافة طاقة الانتقال. تحدد الدراسة حداً أقصى لكثافة طاقة الانتقال عند $e_{tr}^{max} \approx 1.40 \times 10^{15}$ جم/سم$^3$ ($\approx 4.03\rho_{sat}$) لأصلب معادلة حالة كواركية ($c^2_{s,QM}=1.0$).
• القيود على قوة الانقطاع: وُجد أن أقصى انقطاع مسموح به في كثافة الكتلة الراكنة هو $\Delta\rho_{max} \approx 7.76\rho_{sat}$ ($\Delta e_{max} \approx 2.94 \times 10^{15}$ جم/سم$^3$). هذا الحد محكوم بمعادلات الحالة التي تنتج نجوم توأم من الفئة الثانية (حيث يدعم الفرع الثاني $>2.0 M_\odot$ بينما لا يدعم الفرع الثالث ذلك).
• الكتلة القصوى للنجوم التوأم: تم حساب الكتلة القصوى التي يمكن أن تصل إليها النجوم التوأم لتكون $M_{TS}^{max} \approx 2.05 M_\odot$ (عند $c^2_{s,QM}=1.0$ وثقة $2\sigma$). يضع هذا الحد نجوم توأم الفئة الأولى (حيث يدعم كلا الفرعين $>2.0 M_\odot$).
• تأثير مستويات الثقة: النتائج حساسة لمستوى الثقة للاستدلالات الرصدية. عند مستوى ثقة صارم قدره $1\sigma$، يتم استبعاد نجوم توأم الفئة الأولى والثانية تماماً، مما يترك فقط حلول الفئة الثالثة والرابعة ذات الكتل القصوى الأقل بكثير ($\approx 1.37 M_\_\odot$). وعلى العكس من ذلك، عند $3\sigma$، تسترخي القيود، مما يسمح بنطاق أوسع من المعلمات، بما في ذلك الحالات التي يقع فيها ضغط الانتقال في المنطقة "غير المحدودة".
• الاعتماد على الفئة: تصنف الدراسة النجوم التوأم بناءً على الكتل القصوى للفرعين الثاني ($M_{TOV,2}$) والثالث ($M_{TOV,3}$). وتجد أن نجوم توأم الفئة الأولى (الضخمة في كلا الفرعين) مسموح بها فقط عند كثافات انتقال عالية جداً ($\approx 4\rho_{sat}$)، بينما تفضل كثافات الانتقال الأقل الفئات الثانية والثالثة والرابعة.

الأهمية
يؤكد البحث أن رصد النجوم التوأم، وخاصة الضخمة منها، سيكون بمثابة "الدليل القاطع" (smoking gun) على وجود انتقال طوري قوي من الدرجة الأولى من المادة الهادرونية إلى المادة الكواركية في نوى النجوم النيوترونية. ومن خلال رسم منحنيات "قبعة الساحرة" وتطبيق القيود الرصدية الصارمة، توفر الدراسة حدوداً ثرموديناميكية عليا محددة لكثافة الانتقال ($e_{tr}$)، وقوة الانقطاع ($\Delta e$)، والكتلة القصوى للنجوم التوأم. توفر هذه الحدود إطاراً نقدياً لتفسير البيانات الفلكية المستقبلية وتضييق نطاق الطبيعة المحتملة للمادة الكثيفة. ويشير المؤلفون إلى أنه على الرغم من أن تحليلهم يعتمد على التوازن الهيدروستاتيكي، إلا أن وجود نجوم توأم من الفئة الثانية هش للغاية؛ إذ تشير الدراسات الديناميكية السابقة إلى أنها قد تكون غير مستقرة ضد الاضطرابات، مما يؤدي إلى الهجرة إلى الفروع الهادرونية أو الانهيار إلى ثقوب سوداء. وبالتالي، فإن تحديد معادلات الحالة القابلة للنجوم التوأم هو خطوة ضرورية، وإن كانت قد لا تكون كافية، نحو تأكيد وجود المادة الكواركية في النجوم النيوترونية.