Generalized Non-linear Bayesian Pulsar Timing with Enterprise

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

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

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

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

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

بالنسبة للنباض PSR J1600–3053، وجد الفريق أن طريقتهم الجديدة وفرت قيوداً أكثر إحكاماً على كتلة النجم. كانت القياسات السابقة تشير إلى كتلة تبلغ حوالي 2.0 من الكتلة الشمسية، ولكن مع هامش خطأ كبير. وقد قلص التحليل الجديد هذا النطاق بشكل كبير، مما وفر تقديراً أكثر دقة. وبالمثل، بالنسبة لـ PSR J0740+6620، وهو نباض معروف بامتلاكه واحدة من أعلى الكتل التي تم قياسها على الإطلاق، كانت نتائج الفريق متسقة مع النتائج السابقة ولكنها قدمت قيوداً أكثر حدة. لقد أكدوا أن هذا النجم يمتلك كتلة تبلغ حوالي 2.06 ضعف كتلة شمسنا، مع هامش خطأ صغير جداً. هذه الدقة حيوية لأنها تساعد في استبعاد نظريات معينة حول ما يحدث للمادة داخل نجم نيوتروني. فإذا كانت الكتلة عالية جداً، فإن ذلك سيخالف قوانين الفيزياء كما نفهمها حالياً، مما يشير إلى ضرورة مراجعة نظرياتنا حول القوة النووية القوية.

كانت النتائج الأكثر دراماتيكية من النباض PSR J1640+2224. في الدراسات السابقة، أعطى هذا النباض تقديراً للكتلة مرتفعاً بشكل مفاجئ، حيث كان يوحي أحياناً بقيمة تزيد عن ثلاث مرات كتلة الشمس، وهو أمر مستحيل فيزيائياً لنجم نيوتروني. اشتبه الباحثون في أن الكتلة العالية كانت نتاجاً لأسلوب النمذجة وليس خاصية حقيقية للنجم. ولاختبار ذلك، فرضوا حداً فيزيائياً للكتلة في نموذجهم، حيث قيدوها بأن تكون أقل من ثلاث كتل شمسية. وعندما فعلوا ذلك، انخفضت الكتلة المقدرة إلى قيمة أكثر منطقية تبلغ حوالي 2.2 كتلة شمسية. يسلط هذا النتيجة الضوء على نقطة حرجة: إذا لم يكن النموذج الرياضي مرناً بما يكفي لاستيعاب التفاعل المعقد بين معلمات التوقيت والضجيج، فقد ينتج نتائج مستحيلة فيزيائياً. ومن خلال السماح للنموذج باستكشاف النطاق الكامل من الاحتمالات ثم تطبيق القيود الفيزيائية، تمكن الفريق من استعادة صورة أكثر واقعية للنجم.

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

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

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

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