Massive hybrid stars within the extended three-flavor quark-meson diquark model

ملخص تقني: النجوم الهجينة الضخمة في نموذج كوارك-ميزون-دي كوارك الثلاثي الممتد

بيان المشكلة
لا يزال التركيب الداخلي للنجوم النيوترونية، وخاصة تلك التي تتجاوز كتلتها $2M_\odot$، يمثل سؤالاً مركزياً في الديناميكا الحرارية لكروموديناميكا الكم (QCD) الكثيفة. وبينما توفر النسبية العامة الإطار العملي (عبر معادلة تولمان-أوبنهايمر-فولكوف) لربط كتلة النجم ونصف قطره بمعادلة الحالة (EoS)، فإن معادلة الحالة نفسها تظل مجهولة عند الكثافات القصوى الموجودة في نوى النجوم. وتواجه نظرية الشبكة لـ QCD حالياً عدم قدرة على حساب الكميات الديناميكية الحرارية عند كثافة باريونية محدودة بسبب "مشكلة الإشارة". وبناءً على ذلك، يعتمد الباحثون على نظريات الحقل الفعال. ويكمن التحدي الكبير في سد الفجوة بين المادة النووية منخفضة الطاقة (التي تصفها نظرية الحقل الفعال الكيرالية $\chi$EFT) ومادة الكوارك عالية الكثافة (التي تصفها نظرية pQCD الاضطرابية). تعاني النماذج الحالية، مثل نموذج نامبو-جيونا-لاسنيو (NJL)، غالباً من عيوب التنظيم عند الكثافة المحدودة أو تفشل في توليد معادلة حالة صلبة بما يكفي لدعم النجوم النيوترونية الضخمة المرصودة دون انتهاك السببية أو شروط pQCD.

المنهجية
يستخدم المؤلفون نموذج كوارك-ميزون-دي كوارك الثلاثي الممتد (EQMD)، وهو نموذج فعال منخفض الطاقة قابل لإعادة التطبيع لـ QCD. تشمل درجات الحرية الفعالة الكواركات، والميزونات السلمية والمنوعة، والميزونات الدي كواركية، والأهم من ذلك، الميزونات المتجهة والمحورية.

1. بناء النموذج: تم بناء كثافة لاغرانج بحيث تحترم تناظرات QCD ثلاثية النكهات ($SU(3)_L \times SU(3)_R \times SU(3)_c \times U(1)_B \times U(1)_A$). قام المؤلفون بتوسيع نموذج الكوارك-ميزون-دي كوارك القياسي عبر إضافة حدود الكتلة والحدود الرباعية لحقول الميزونات المتجهة والمحورية الجديدة ($L_\mu$ و $R_\mu$) بشكل أدنى. ترتبط هذه الحقول الجديدة بالكواركات عبر ثابت اقتران $g_\omega$ وتعمل فعلياً على إزاحة الكيمياء الكيميائية للكواركات.
2. الديناميكا الحرارية: يتم حساب معادلة الحالة (EoS) في تقريب المجال المتوسط عند درجة حرارة صفر ($T=0$). يقوم النموذج بفرض الحياد الشحني المحلي للشحنات الكهربائية واللونية يدوياً، نظراً لافتقاره إلى حقول قياس ديناميكية. يتم إعادة تطبيع الجهد الديناميكي الحراري باستخدام مخطط "على الشاقول" (on-shell scheme)، حيث تُضبط كتل الميزونات وثوابت الاضمحلال لتناسب الملاحظات الفيزيائية.
3. بناء النجم الهجين: لنمذجة النجوم الهجينة، يقوم المؤلفون بربط معادلة الحالة عالية الكثافة لـ EQMD بمعادلة حالة نووية منخفضة الكثافة (تحديداً معادلة GMSR من مكتبة CompOSE). يستخدمون بناءً يتم فيه تحديد الانتقال بين المادة النووية ومادة الكوارك عبر مساواة كثافات الجسيمات ($n_{EQMD} = n_{nucl}$) عند جهد كيميائي انتقالي $\mu_t$، مما يضمن استمرارية كثافة الباريون. يتجنب هذا النهج الانقطاعات الناتجة عن إنشاء ماكسويل، لكنه يسمح بانتقال سلس بين الحالة النووية وحالة الكوارك.
4. فضاء المعلمات: يحتوي النموذج على عدة معلمات حرة، وأبرزها اقتران الكوارك-المتجه $g_\omega$. يختبر المؤلفون ثلاث مجموعات من المعلمات بقيم متفاوتة لـ $g_\omega$ لاستقصاء صلابة معادلة الحالة.

المساهمات والنتائج الرئيسية

• صلابة معادلة الحالة (EoS): يُعتبر إدراج الميزونات المتجهة والمحورية أمراً جوهرياً. تتيح درجات الحرية هذه للنموذج توليد معادلة حالة صلبة بما يكفي عند الكثافات المتوسطة لدعم النجوم الضخمة، وتلين عند الكثافات الأعلى لتظل متوافقة مع pQCD.
• الاتساق مع الملاحظات: تحقق النجوم الهجينة الناتجة القيود الفيزيائية الفلكية الحالية، بما في ذلك:
  • قيود الكتلة القصوى للنباضات PSR J1614-2230 ($1.97 M_\odot$)، و PSR J0348+0432 ($2.01 M_\odot$)، و PSR J2215+5135 ($2.27 M_\odot$)، و PSR J0952-0607 ($2.35 M_\odot$).
  • قيود الكتلة ونصف القطر من ملاحظات NICER للنباضين PSR J0030+0451 و PSR J0740+6620.
  • تتوافق مجموعة معلمات واحدة (المجموعة 1) مع جميع الملاحظات، بينما تتوافق المجموعتان 2 و 3 مع كل شيء باستثناء النباض الأكثر ضخامة PSR J0952-0607.
• الكثافات المركزية والنوى الكواركية: يجد المؤلفون أن النجوم ذات الكتل $M \gtrsim 2M_\odot$ تمتلك نواة كواركية بكثافة باريونية مركزية $n_B \geq 3.9 n_{sat}$ (حيث $n_{sat} \approx 0.165 \, \text{fm}^{-3}$). بالنسبة لتكوينات الكتلة القصوى، تُقدر أنصاف أقطار النوى الكواركية بحوالي 3.8 كم إلى 4.6 كم.
• بنية سرعة الصوت: يظهر مربع سرعة الصوت ($c_s^2$) بنية ذروة مزدوجة:
  • تنشأ الذروة الأولى من الانتقال من معادلة الحالة النووية إلى معادلة حالة مادة الكوارك.
  • يتبع ذلك انخفاض حاد، ناتج عن الذوبان السريع لتكاثف الكوارك الغريب مع زيادة الجهد الكيميائي.
  • تحدث ذروة ثانية عند كثافات أعلى، حيث تتجاوز $c_s^2$ الحد المطابق ($1/3$)، قبل أن تقترب من هذا الحد من الأعلى عند الكثافات النهائية.
  • يتناقض هذا السلوك مع نموذج NJL (حيث $c_s \to 1$) ومع التوسعات ذات الاقتران الضعيف (التي تقترح الاقتراب من الحد من الأسفل)، ولكنه يتوافق مع الديناميكيات المحددة لنموذج EQMD، حيث تقل كتلة الكوارك الغريب مع زيادة $\mu_B$.
• تحديد ثابت الحقيبة (Bag Constant): على عكس نماذج الحقيبة التقليدية حيث يكون ثابت الحقيبة $B$ معلماً حراً، يتم هنا تحديد $B$ بشكل متسق ذاتياً من خلال شرط المطابقة عند كثافة الانتقال. القيم المستمدة ($B^{1/4} \approx 234\text{--}242$ MeV) تتوافق مع تقديرات نماذج NJL والحسابات القائمة على QCD.

الأهمية والادعاءات
يزعم هذا العمل أن نموذج EQMD، وتحديداً مع إضافة الميزونات المتجهة، يقدم إطاراً قابلاً لإعادة التطبيع لوصف طور المادة منزوعة الحجز في النجوم النيوترونية. تكمن الأهمية الأساسية في إثبات إمكانية بناء نجوم هجينة تحقق في آن واحد:

1. الصلابة المطلوبة لدعم نجوم بكتلة $2M_\odot$.
2. اللين الذي تطلبه قيود pQCD عند الكثافات العالية.
3. قيود $\chi$EFT عند الكثافات المنخفضة.

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