Wide parameter-space O3 search for continuous gravitational waves from unknown neutron stars in binary systems

ملخص تقني: البحث في فضاء بارامترات واسع لبيانات O3 عن موجات جاذبية مستمرة من نجوم نيوترونية مجهولة في أنظمة ثنائية

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

المنهجية
أجرى المؤلفون بحثاً شبه متماسك (semi-coherent search) باستخدام بيانات عامة من دورة الرصد الثالثة (O3) لكواشف لايف (LIGO) المتقدمة (H1 و L1). تم تنظيف البيانات من الضوضاء غير الغاوسية عند التوافقيات الـ 60 وخطوط المعايرة، كما تمت إزالة النبضات (glitches) قصيرة المدى عبر عملية التغطية (gating).

• نموذج الإشارة: استخدم البحث نموذج إشارة CW الخاص بـ Jaranowski et al. (1998)، والذي يتضمن بارامترات السعة الذاتية ($h_0, \cos \iota, \psi, \phi_0$) وبارامترات تطور الطور (التردد $f_0$، التباطؤ في الدوران $f_1$، الموقع في السماء، والبارامترات المدارية الثنائية: $P_{orb}, a_p, t_{asc}, \omega, e$).
• فضاء البارامترات: غطى البحث ترددات موجات الجاذبية من 50 هرتز إلى 1,000 هرتز. ولإدارة التكاليف الحوسبية، تم تقسيم فضاء البارامترات المدارية $\{P_{orb}, a_p\}$ إلى خمس مناطق متميزة (من A إلى E). ومن الجدير بالذكر أن المنطقة E ($P_{orb} \in [0.2, 0.4]$ يوماً، $a_p \in [0.03, 0.12]$ ثانية ضوئية) لم تُستكشف تماماً قبل هذا العمل. افترض البحث مدارات ذات انحرافات صغيرة ($e \leq 0.1$) وحدود تباطؤ في الدوران تتوافق مع مجتمع النجوم النيوترونية الثنائية المرصودة.
• مسار البحث (Pipeline): اعتمد التحليل على مسار البحث شبه المتماسك BinarySkyHouF.
  • المرحلة الأولية: تم تقسيم البيانات إلى $N_{seg} = 31,675$ جزءاً، طول كل منها $T_{seg} = 900$ ثانية. تم حساب إحصائية الكشف ($2\hat{F}_{AB}$) باستخدام مزيج خطي من الإحصائيات من القوالب ذات البارامترات المدارية الصفرية، والموزونة لتتبع الأنماط الزمنية-الترددية بما في ذلك البارامترات المدارية غير الصفرية. تم تحسين تباعد الشبكة للحفاظ على نسبة عدم التطابق (mismatch) دون 2%.
  • المتابعة (Follow-up): تم اختيار المرشحين بناءً على الأهمية وخضعوا لمتابعة متماسكة (coherent follow-up) بمدة $T_{seg} = 2,700$ ثانية. استُخدمت خوارزمية أخذ العينات المتداخلة (nested sampling) لحساب إحصائية الكشف وإحصائية المتانة ($\log_{10} \hat{B}_{1;S/GLtL}$) للتمييز بين الإشارات الفيزيائية الفلكية والخطوط الأدواتية والضوضاء.
• التحقق: تم حقن إشارات اختبار محاكية في البيانات لمعايرة عتبات الكشف وتقدير الحدود العليا. تم فحص المرشحين مقابل الحقن العتادي (hardware injections) والاضطرابات غير الغاوسية.

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

1. نطاق ترددي ممتد: هذا هو أول بحث يغطي ترددات موجات الجاذبية فوق 520 هرتز، وصولاً إلى 1,000 هرتز. وهذا أمر مهم بسبب سيناريو "إعادة التدوير" (recycling) حيث تكتسب النجوم النيوترونية زخماً زاوياً عبر التراكم، مما قد يؤدي للوصول إلى ترددات دوران عالية.
2. فترات مدارية قصيرة: هذا هو أول بحث يستخدم الكواشف المتقدمة لاستكشاف فترات مدارية أقل من 3 أيام، مستهدفاً تحديداً المنطقة E ($P_{orb} \approx 0.2-0.4$ يوماً).
3. اتساع فضاء البارامترات: يغطي البحث فضاء بارامترات أكبر بنحو أربع مراتب عشرية من أي تحقيق سابق باستخدام بيانات Advanced LIGO.
4. الاستراتيجية الحوسبية: استخدم المؤلفون نهج "البحث أولاً بالعرض" (breadth-first)، من خلال استخدام شبكات قوالب أخشن وفترات أساس متماسكة أقصر مقارنة بالبحوث المحسنة من حيث الحساسية، مما سمح بتغطية مناطق شاسعة وغير مستكشفة سابقاً ضمن ميزانية حوسبية ثابتة.

النتائج

• الكشف: لم يتم اكتشاف أي إشارات موجات جاذبية مستمرة فيزيائية فلكية.
• المرشحون: من أصل حوالي 127,500 مرشح، صمد 15 فقط أمام عمليات الفيتو (vetoes) للمتابعة. ومع ذلك، تم تحديد 12 منها كإشارات مرتبطة بالحقن العتادي (إشارات نجوم نيوترونية معزولة عند ترددات محددة)، ونُسبت الثلاثة المتبقية إلى اضطرابات غير غاوسية. لم تُعتبر أي منها إشارات CW يمكن الدفاع عنها.
• الحدود العليا: وضع المؤلفون أشد القيود صرامة حتى الآن على سعة مثل هذه الإشارات.
  • عمق الحساسية: تم تحديد متوسط عمق الحساسية بـ $\langle D_{95\%}^{low} \rangle = 18.1 \pm 1.1$ هرتز$^{-1/2}$ (50–500 هرتز) و $\langle D_{95\%}^{high} \rangle = 17.1 \pm 0.9$ هرتز$^{-1/2}$ (500–1000 هرتز).
  • قيود الشكل البيضاوي (Ellipticity): بالنسبة للنجوم النيوترونية التي تبعد 100 فرسخ فلكي وتدور بسرعة تزيد عن ~495 هرتز، يُستبعد وجود أشكال بيضاوية أعلى من $5.2 \times 10^{-8}$ بثقة 95% (بافتراض $I_{zz} = 10^{38}$ كجم م$^2$).
  • قيود نمط r-mode: بالنسبة لنفس المسافة ومعدلات الدوران (> ~740 هرتز)، يُستبعد وجود سعة نمط r-mode أعلى من $1.5 \times 10^{-6}$.

الأهمية
تدعي هذه الورقة أن هذا العمل يمثل البحث الأكثر شمولاً حتى الآن عن الموجات المستمرة للجاذبية من النجوم النيوترونية المجهولة في الأنظمة الثنائية. ومن خلال دفع حد التردد إلى 1,000 هرتز واستكشاف فترات مدارية تصل إلى 0.2 يوم، يعالج البحث فجوات حرجة في فضاء البارامترات حيث يُتوقع وجود النجوم النيوترونية المعاد تدويرها والمتراكمة. ورغم عدم العثور على إشارات، إلا أن النتائج تقدم أضيق القيود حتى الآن على الشكل البيضاوي وسعة نمط r لهذه المصادر. يضع المؤلفون هذا البحث كـ "مخطط" (blueprint) لاستقصاء مجموعات البيانات الجديدة الأكثر حساسية، موضحين قدرات التقنيات المتطورة في نظام الترددات العالية الصعب رغم التوسع الحوسبي ($\propto f_0^5$). ويقر العمل بأن مناطق شاسعة من فضاء البارמترات لا تزال غير مستكشفة، لكنه يضع إطاراً للبحوث المستقبلية الأكثر حساسية.