Testing the strong equivalence principle with multimessenger binary neutron star mergers

1. بيان المشكلة

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

تتناول هذه الورقة إمكانية وجود تغير زمني في ثابت الجاذبية ($\dot{G} \neq 0$). ومن شأن هذا التغير أن ينتهك مبدأ التكافؤ القوي ويظهر بطريقتين متمايزتين أثناء اندماج النجوم النيوترونية الثنائية (BNS):

1. تأثيرات المصدر: يؤدي $G$ المتغير زمنياً إلى تعديل البنية الداخلية (الطاقة الذاتية للجاذبية) والديناميكا المدارية للثنائيات المتراصة، مما يغير تطور طور الموجات الثقالية (GW).
2. تأثيرات الانتشار: يؤدي تغير $G$ عبر المسافات الكونية إلى تعديل سعة الموجات الثقالية أثناء انتقالها عبر كون فريدمان-لوميتر-روبرتسون-ووكر (FLRW) متمدد.

غالباً ما كانت الدراسات السابقة تعالج هذه التأثيرات بشكل منفصل، أو تعتمد على تقريبات مصفوفة فيشر (Fisher matrix)، أو تفتقر إلى قيود قوية من بيانات متعددة الرسل (multimessenger) حقيقية. يهدف هذا العمل إلى إجراء تحليل بايزي (Bayesian analysis) مشترك وصارم باستخدام بيانات حقيقية من اندماج النجوم النيوترونية الثنائية GW170817 ونظيره الكهرومغناطيسي GRB 170817A لتقييد $\dot{G}/G$.

2. المنهجية

أ. الإطار النظري

قام المؤلفون بتطوير نموذج موجة ثقالية متسق فيزيائياً يتضمن تغيراً بطيئاً في $G$ ضمن إطار ما بعد أينشتاين البارامتري (ppE):

• ديناميكا المصدر: اشتقوا التطور المداري لنظام ثنائي حيث تتطور الكتلة الكلية $M$ والزخم الزاوي $j$ بسبب $\dot{G}$. وهذا يستحدث تصحيحات لطور وسعة الموجات الثقالية تتميز بمعلمات ppE وهي $\alpha$ (السعة) و $\beta$ (الطور).
  • يتناسب تصحيح الطور مع $u^{-13}$ (حيث $u = (\pi M_c f)^{1/3}$).
  • يتناسب تصحيح السعة مع $u^{-8}$.
• تأثيرات الانتشار: أثبتوا أنه في كون FLRW مع متغير $G$، تكتسب سعة الموجات الثقالية عامل تصحيح $\sqrt{G_0/G_*}$، حيث $G_*$ هو ثابت الجاذبية عند المصدر و $G_0$ عند الكاشف.
• الحساسيات ($s$): تعتمد magnitude تأثيرات المصدر على "حساسية" الأجرام المتراصة ($s = -\partial \ln m / \partial \ln G$).
  • للنجوم النيوترونية (NS): تم حساب الحساسيات عددياً عن طريق حل معادلات تولمان-أوبنهايمر-فولكوف (TOV) لمختلف معادلات الحالة (EOS). أظهر المؤلفون أن التقريب الخطي التجريبي المستخدم غالباً في الأدبيات صالح فقط عند الكتل المنخفضة؛ وتتطلب الحساسيات الدقيقة حلولاً عددية تعتمد على معادلة الحالة (EOS).
  • للثقوب السوداء: قدموا برهاناً مستقلاً عن النموذج بأن الحساسية هي عالمياً $s=1/2$ للثقوب السوداء ذات التماثل الكروي في نظريات $G$ المتغيرة.

ب. بناء الموجة

قام المؤلفون ببناء قالب موجة ثقالية مخصص ($\tilde{h}_{VG}$) يجمع بين:

1. موجات النسبية العامة القياسية (المولدة عبر PyCBC باستخدام نماذج TaylorF2 و IMRphenomXAS_NRTidalv2).
2. حد تصحيح الطور $\exp(i \beta_{VG} u^{b_{VG}})$ المستمد من ديناميكا $G$ المتغيرة.
3. عامل تصحيح السعة $1/\sqrt{\kappa}$، حيث $\kappa \equiv G_*/G_0$.
4. القيود الفيزيائية: يفرض النموذج أن "كتلة معادلة الحالة" (كتلة معادلة إعادة قياس تعتمد على $\kappa$) يجب أن تقع ضمن النطاق المسموح به فيزيائياً للنجوم النيوترونية (0.08 إلى 2.01 $M_\odot$) بناءً على معادلة الحالة H4.

ج. الاستدلال البايزي

• البيانات: بيانات الإجهاد الثقالي من LIGO Hanford و Livingston (GW170817) مع قطع ترددي عند 20 هرتز.
• قيود تعدد الرسل: لكسر التحلل في المعلمات (parameter degeneracies) المتأصل في التحليلات التي تعتمد على الموجات الثقالية فقط (تحديداً بين $\kappa$، ومسافة المسافة السطوعية $D_L$، والميل $\iota$)، أدرج المؤلفون قيوداً كهرومغناطيسية مستقلة:
  • مسافة السطوع ($D_L$): مثبتة عند 40.7 ميجابارسيك (من تقلبات سطوع السطح لتلسكوب هابل).
  • الموقع في السماء: مثبت بناءً على ملاحظات كاميرا الطاقة المظلمة (Dark Energy Camera).
  • الميل ($\iota$): مقيد بالفترة الموثوقة 90% $[2.70, 2.81]$ راديان باستخدام تداخل خط القاعدة الطويل جداً (VLBI) للمتوهج اللاحق.
• الاستدلال: تم إجراء تحليل بايزي مشترك لاستنتاج التوزيعات البعدية للمعلمات القياسية للنسبية العامة والمعلمة غير النسبية $\kappa$. ويتم اشتقاق المشتق الزمني الكسري كـ $\dot{G}/G = (1-\kappa)/T$، حيث $T$ هو زمن النظر إلى الخلف (lookback time) وقدره $1.31 \times 10^8$ سنة.

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

1. إطار موحد: تطوير نموذج موجة ثقالية شامل يأخذ في الاعتبار في وقت واحد تأثيرات توليد المصدر (الديناميكا المدارية، الحساسيات) وتأثيرات الانتشار (تخميد السعة) لثابت جاذبية متغير.
2. حساسيات تعتمد على معادلة الحالة (EOS): تجاوزوا التقريبات التجريبية عبر حساب حساسيات النجوم النيوترونية عددياً لعدة معادلات حالة واقعية (AkmalPR, AP4, SLy4, إلخ)، مما يضمن الاتساق الفيزيائي لحدود الكتلة المستخدمة في التحليل.
3. كسر التحلل عبر تعدد الرسل: أظهروا القدرة الحاسمة لدمج بيانات الموجات الثقالية مع النظائر الكهرومغناطيسية المستقلة (تحديداً الميل والمسافة) لفك الارتباط بين تأثيرات $G$ المتغير والمعلمات الفيزية الفلكية القياسية.
4. استدلال بايزي صارم: إجراء تقدير كامل للمعلمات البايزية باستخدام بيانات حقيقية (GW170817) بدلاً من الاعتماد على إسقاطات مصفوفة فيشر أو محاكاة مونت كارلو، مع استخدام نموذجين مختلفين للموجات لضمان المتانة.

4. النتائج

• لا يوجد دليل على التغير: لم يجد التحليل أي دليل ذي دلالة إحصائية على تغير زمني في ثابت الجاذبية. التوزيعات البعدية لمعلمة $G$ المتغيرة $\kappa$ تتوافق مع تنبؤ النسبية العامة ($\kappa = 1$).
• قيود المعلمات:
  • باستخدام نموذج TaylorF2: $\kappa = 1.17^{+0.27}_{-0.19}$.
  • باستخدام نموذج IMRphenomXAS_NRTidalv2: $\kappa = 1.11^{+0.25}_{-0.18}$.
• الحد النهائي لـ $\dot{G}/G$:
  من خلال تحويل القيود على $\kappa$ باستخدام زمن النظر إلى الخلف، وضع المؤلفون حداً تحفظياً عند مستوى 3$\sigma$:
  $$\frac{\dot{G}}{G} \in [-3.36 \times 10^{-9}, 5.34 \times 10^{-10}] \, \text{yr}^{-1}$$
  ويمثل هذا أشد قيود تم استخلاصها حتى الآن من ملاحظات الموجات الثقالية الحقيقية.

5. الأهمية

• اختبار دقيق لمبدأ التكافؤ القوي (SEP): يوفر هذا العمل أضيق اختبار للحقل القوي لمبدأ التكافؤ القوي باستخدام اندماجات الثنائيات المتراصة، مما يحسن بشكل كبير من الحدود السابقة القائمة على الموجات الثقالية (والتي كانت أضعف بعدة مراتب أو تعتمد على بيانات متوقعة).
• التحقق من النسبية العامة: تعزز النتائج صحة النسبية العامة في نظام الحقل القوي والديناميكي، مؤكدة أن $G$ ثابت ضمن نطاق $\sim 10^{-9}$ في السنة عبر المقاييس الكونية.
• التآزر متعدد الرسل: تسلط الدراسة الضوء على أن الإمكانات الكاملة لعلم الفلك متعدد الرسل (الجمع بين الموجات الثقالية والنظائر الكهرومغناطيسية) ضرورية للفيزياء الأساسية الدقيقة، حيث تسمح بكسر حالات التحلل التي قد تجعل مثل هذه الاختبارات مستحيلة باستخدام بيانات الموجات الثقالية وحدها.
• معيار منهجي: يعمل الإطار المطور وحسابات الحساسية العددية كنموذج متين للاختبارات المستقبلية للفيزياء الأساسية باستخدام كواشف الموجات الثقالية القادمة (مثل Einstein Telescope أو Cosmic Explorer) أو أحداث متعددة الرسل أكثر دقة.