The multimessenger view of Pulsar Timing Array black holes with the Horizon-AGN simulation

بيان المشكلة
لا تزال أصول وتطور الثقوب السوداء فائقة الكتلة (MBHs) وتطورها غير مؤكدة، لا سيما فيما يتعلق بدمجها في الكون المحلي. وبينما رصدت مصفوفات التوقيت النبضي (PTAs) مؤخراً إشارة خلفية موجات جاذبية (GWB) تتوافق مع تراكب عشوائي لثنائيات الثقما السوداء الضخمة (MBHBs)، إلا أن هناك توتراً مستمراً بين السعة المرصودة والتوقعات النظرية، والتي غالباً ما تفضل سعات أقل. كما أن تفسير طيف خلفية موجات الجاذبية (GWB) يزداد تعقيداً بسبب المدة المحدودة لمجموعات البيانات واحتمال وجود مصادر موجات مستمرة (CW) قابلة للحل بشكل فردي وتبرز فوق الخلفية. تعتمد الطرق الحالية لمقارنة النماذج النظرية مع بيانات مصفوفات التوقيت النبضي (PTA) على افتراضات مبسطة قد تؤدي إلى تحيزات. تهدف هذه الورقة إلى معالجة هذه القضايا من خلال استخدام محاكاة "Horizon-AGN" الكونية لدراسة خصائص ثنائيات الثقما السوداء الضخمة (MBHBs) المساهمة في نطاق مصفوفات التوقيت النبضي، وتوليد مجموعات واقعية لتوقعات خلفية موجات الجاذبية (GWB) والموجات المستمرة (CW)، وتقييم منهجيات استنتاج الخصائص الطيفية لخلفية موجات الجاذبية من بقايا التوقيت بشكل نقدي.

المنهجية
تستخدم الدراسة محاكاة "Horizon-AGN" الهيدروديناميكية الكونية، والتي تنمذج حجماً قدره $142 \, \text{cMpc}^3$ ضمن كون $\Lambda$CDM القياسي. قام المؤلفون بإنشاء فهرس لاندماجات الثقما السوداء الضخمة عبر تتبع معرفات (IDs) هذه الثقوب عبر المحاكاة وتحديد عمليات الاندماج الرقمي. وللأخذ في الاعتبار التأخير الفيزيائي بين الاندماج الرقمي (عند مسافة انفصال تبلغ $\sim 4 \, \text{kpc}$) والالتحام الفعلي، طبقوا خط معالجة لاحقة يتضمن:

1. الاحتكاك الديناميكي: تم حسابه بناءً على الاحتكاك الديناميكي النجمي ضمن تقريب الكرة المتساوية الحرارة.
2. تصلب الثنائي: تم نمذجته عبر التصلب النجمي، وعزوم الدوران اللزوجية في الأقراص المحيطة بالثنائي، وانبعاثات موجات الجاذبية (GW)، باتباع وصفات راسخة (مثل Sesana & Khan 2015؛ Dotti et al. 2015).
3. تطور اللامركزية: تم تتبعه عبر تشتت فقدان المخروط، والسحب اللزوجي، وانبعاثات موجات الجاذبية، مع تعيين قيم أولية للامركزية بناءً على نتائج عددية سابقة.

لتقدير خلفية موجات الجاذبية (GWB)، استخدم المؤلفون نهجين متميزين:

• نهج المجموعة الغاوسية (GE): حساب حتمي لطيف شدة الخصائص بافتراض توزيع مستمر للمصادر، وهو نهج يُستخدم غالباً في الأدبيات.
• نهج المجموعات السكانية: طريقة عشوائية تولد 2,000 "تحقيق للكون" عن طريق سحب ثنائيات في حالة اندماج من توزيع بواسون بناءً على فهرس الاندماج. وهذا يأخذ في الاعتبار التباين الكوني والطبيعة المنفصلة للمصادر.

قام المؤلفون بمحاكاة بقايا توقيت مصفوفات التوقيع النبضي (PTA) لهذه المجموعات باستخدام طريقتين:

1. المجموعة الغاوسية: تحويل شدة الخصائص إلى كثافة الطيف القدرتي (PSD) لبقايا التوقيت عبر استجابة نبض متوسطة.
2. نهج المجموعات السكانية: جمع مساهمات شدة موجات الجاذبية الفردية لـ "طرف الأرض" من جميع الثنائيات في تحقيق معين، مع نمذجة صريحة لمدة المراقبة ($T_{\text{obs}}$) وآثار التسرب الطيفي عبر دالات sinc.

وقارنوا هذه التوقعات النظرية مع مجموعات بيانات "لعبة" (toy datasets) تمثل حساسية EPTA الحالية (25 نبضاً، 10.33 سنة) وحساسية SKA المستقبلية (80 نبضاً، 10.33 سنة). كما حددوا مرشحي الموجات المستمرة (CW) المحتملين (ثنائيات فردية ذات نسبة إشارة إلى ضجيج $S/N > 3$) وحللوا خصائصهم الكهرومغناطيسية (البصرية، والأشعة تحت الحمراء القريبة، والراديو، والأشعة السينية) من خلال نمذجة انبعاث النوى المجريّة النشطة (AGN) وخصائص المجرات المضيفة.

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

• استنتاج طيف GWB والتحيز: تسلط الدراسة الضوء على تحفظ منهجي هام عند مقارنة النماذج النظرية ببيانات PTA. النهج القياسي للمجموعة الغاوسية، الذي يفترض مدة مراقبة لانهائية، يفشل في مراعاة التسرب الطيفي الناتج عن مدة المراقبة المحدودة ($T_{\text{obs}}$). عندما يُستخدم نهج المجموعات السكانية مع تحليل واقعي لـ PTA (بما في ذلك تهميش نموذج التوقيت)، فإن طيف GWB المستنتج ينحاز نحو منحدرات أكثر ضحالة وسعات أعلى عند التردد $f = 1/\text{yr}$ مقارنة بالطيف الحقيقي الأساسي. هذا التحيز، المدفوع بالتسرب الطيفي من المصادر منخفضة التردد، قد يفسر جزئياً التوتر بين التوقعات النظرية وسعات PTA المرصودة.
• العشوائية وأعداد المصادر: عند الترددات المنخفضة ($\sim 3 \, \text{nHz}$)، تهيمن مئات أو آلاف الثنائيات على GWB، مما يؤكد صحة تقريب الخلفية العشوائية. ومع ذلك، عند الترددات الأعلى ($\sim 10\text{--}50 \, \text{nHz}$)، ينخفض عدد المصادر المساهمة بشكل كبير (إلى بضع عشرات)، مما يؤدي إلى فقدان العشوائية وزيادة التباين الكوني. تحافظ المجموعات ذات اللامركزية على العشوائية إلى ترددات أعلى من المجموعات الدائرية بسبب الانبعاث في التوافقيات المدارية الأعلى.
• خصائص المصادر (الخلفية مقابل CW):
  • الخلفية: تهيمن على GWB ثنائيات بكتلة تشير (chirp masses) تبلغ $M_c \sim 10^{8.5} - 10^{9.5} \, M_\odot$ عند انزياحات حمراء $z \in 0.05 - 1$. هذه الثنائيات تسكن أساساً في مجرات ضخمة ($M_\star > 10^{11} \, M_\odot$) ضمن هالات كتلتها $\sim 10^{13} \, M_\odot$.
  • مرشحو CW: تميل المصادر الفردية القابلة للحل (مرشحو CW) إلى أن تكون أكثر كتلة وتقع في انزياحات حمراء أقل ($z < 0.4$) مقارنة بالخلفية العامة. وهي تسكن غالباً في المجرات المركزية للمجموعات والعناقيد الضخمة ($M_{\text{halo}} > 10^{13} \, M_\odot$).
• النظائر الكهرومغناطيسية:
  • الطبيعة: تظهر غالبية مرشحي CW كنوى مجرية نشطة (AGN) بدلاً من كوازارات، وتتميز بنسب إدنجتون منخفضة ($f_{\text{Edd}} < 10^{-2}$).
  • القابلية للرصد: في نطاقات الضوء المرئي والأشعة تحت الحمراء القريبة (NIR)، تطغى المجرات المضيفة على معظم مرشحي CW، مما يجعل عمليات البحث عن تقلبات منحنى الضوء صعبة ومنحازة للأنظمة المضيئة والنادرة.
  • أفضل الأطوال الموجية: توفر نطاقات الراديو والأشعة السينية أفضل الآفاق لتحديد الهوية. بالنسبة للمرشحين الذين تبلغ نسبة $f_{\text{Edd}} > 10^{-3}$، يمكن لانبعاث الـ AGN أن يتفوق على المجرة المضيفة في هذه النطاقات. حتى بالنسبة للأنظمة ذات التراكم المنخفض، تشير التقديرات المتفائلة إلى أن ما يصل إلى 80% (راديو) و30% (أشعة سينية) يمكن اكتشافها فوق مستوى المجرة المضيفة.
• احتمالات الكشف: بافتراض خط أساس مدته 10 سنوات، تُقدر احتمالية اكتشاف مرشح واحد على الأقل من نوع CW بنسبة $S/N > 3$ بنسبة 4% لحساسية شبيهة بـ EPTA و20% لحساسية شبيهة بـ SKA.

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