Search for Long-Transient Gravitational Waves from Supernova SN2023ixf using GFH-v2 Pipeline

ملخص تقني: البحث عن موجات الجاذبية طويلة العبور من SN 2023ixf باستخدام GFH-v2

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

يمثل SN 2023ixf، وهو مستعر أعظم من النوع الثاني (Type II) ناتج عن انهيار قلب تم اكتشافه في المجرة القريبة M101 (على مسافة $\approx$ 6.85 ميجابارسيك)، فرصة فريدة لمثل هذا البحث. إذ يتيح قربه ودقة توقيت رصده الكهرومغناطيسي (ضمن نافذة زمنية تبلغ حوالي ساعتين) إجراء بحث موجه. والأهم من ذلك، أن الحقبة الزمنية للحدث تتداخل مع فترة التشغيل الهندسي 15 (ER15) لمرصد LIGO، مما يوفر مجموعة بيانات متزامنة من كاشفي Hanford (H1) وLivingston (L1) تسبق مباشرة دورة الرصد الرابعة (O4).

المنهجية
يستخدم التحليل مسار المعالجة GFH-v2، وهو طريقة بحث شبه متماسكة هرمية تعتمد على تحويل "هوف" للتردد المعمم (Generalized Frequency Hough transform)، وهي مصممة خصيصاً للإشارات ذات التطور الترددي الذي يتبع قانون القوة.

• نموذج الإشارة: يستهدف البحث نجماً نيوترونياً صلباً غير مسبوق في الدوران، يتباطأ دورانه بشكل أساسي عبر انبعاث موجات الجاذبية. وهذا يتوافق مع مؤشر كبح $n=5$. يتطور تردد موجة الجاذبية $f_{gw}$ كالتالي: $\dot{f}_{gw} = -k f_{gw}^n$، وتتناقص سعة الشدة $h_0(t)$ كـ $t^{-1/2}$. يفترض البحث بارامترات نموذجية للنجم النيوتروني ($M=1.4 M_\odot$، $R=12$ كم) وموقعاً سماوياً ثابتاً بناءً على التحديد الموضعي الكهرومغناطيسي.
• البيانات: يستخدم التحليل قطاعات علمية متزامنة من H1 وL1 خلال فترة ER15، بإجمالي حوالي 48,619 ثانية (0.56 يوم) من البيانات ضمن نافذة المصدر المحافظة. وقد تم التحقق من جودة البيانات لتكون مقاربة لحساسية O4a في نطاق 600–2000 هرتز.
• إعدادات البحث: تغطي مساحة البارامترات الترددات الأولية ($f_0$) من 600 إلى 2000 هرتز (بخطوة 10 هرتز) واللاتمركزية ($\epsilon$) من $3 \times 10^{-4}$ إلى $3 \times 10^{-3}$ (مقسمة إلى خمس فترات فرعية). ينتج عن ذلك 700 تكوين بحث مستقل.
• خطوات المسار:
  1. تحضير البيانات: يتم تمرير بيانات قاعدة بيانات تحويل فوريه قصير المدى (SFDB) عبر مرشح تمرير حزمة (band-pass) وتحويلها إلى سلاسل زمنية تحليلية معقدة.
  2. توليد خريطة القمم (Peakmap): يتم اختيار القيم القصوى المحلية في طيف القدرة المتساوي لتشكيل خريطة قمم زمنية-ترددية.
  3. تحويل هوف (Hough Transform): يتم رسم خريطة للقمم من الإحداثيات $(t, f)$ إلى إحداثيات محولة $(t, x)$ حيث $x = 1/f^{n-1}$. في هذا الفضاء، يصبح التطور بقانون القوة خطياً. يقوم تحويل هوف برسم هذه الخريطة إلى فضاء البارامترات $(x_0, k)$ (التردد الأولي ومعامل التباطؤ).
  4. اختيار المرشحات: يتم تحديد المرشحات كفئات تحتوي على زيادات في العد بالنسبة للخلفية المحلية. ويتم حساب النسبة الحرجة (CR) باستخدام إحصاءات الضوضاء المحلية المقدرة (الوسيط والانتشار ضمن فترات 10 هرتز) لمراعاة الضوضاء غير المستقرة.
  5. التطابق (Coincidence): تتم مقارنة المرشحات من H1 وL1 في فضاء البارامترات $(x_0, k)$. يُعتبر الزوج متطابقاً إذا كانت مسافة التطابق $d_{coin} \le 2$.
  6. المتابعة (Follow-up): تخضع المرشحات المتطابقة التي تحقق $CR \ge 3.5$ لمتابعة باستخدام تقنية الهيتروداين (heterodyne). يتم إزالة تطور الطور المتوقع من البيانات، ثم تُحلل الإشارة مرة أخرى مع أوقات تماسك أطول ($3\times$ و $6\times$ مدة التماسك الأصلية) للتحقق من زيادة الأهمية، مما يشير إلى وجود إشارة فلكية حقيقية.

النتائج الرئيسية

• اختيار المرشحات: بعد تطبيق شرط التطابق وعتبة النسبة الحرجة ($CR \ge 3.5$)، تبقى 131 مرشحاً متطابقاً.
• نتيجة المتابعة: لم تظهر أي من المرشحات الـ 131 الزيادة المتوقعة في النسبة الحرجة أثناء مراحل متابعة الهيتروداين. وبناءً على ذلك، لم يتم ترقية أي مرشح إلى مرحلة المتابعة الثانية.
• الاكتشاف: لم يجد التحليل أي دليل على انبعاث موجات جاذبية طويلة العبور مرتبطة بـ SN 2023ixf في بيانات ER15.
• الحدود العليا: تم وضع حدود عليا لأقصى مسافة يمكن اكتشافها كدالة للتردد الأولي واللاتمركزية. بالنسبة لأعلى فئة لاتمركزية ($\epsilon \in [2.5 \times 10^{-3}, 3.0 \times 10^{-3}]$)، تتراوح الحدود العليا لكفاءة 90% بين 1 إلى 2.5 ميجابارسيك تقريباً عبر معظم نطاق التردد الذي تم تحليله. بالنسبة للقيم الأقل من اللاتمركزية، كانت الحدود إما غير قابلة للتحديد أو مقتصرة على مسافات أصغر بسبب محدودية مدة البيانات وانخفاض نسب الإشارة إلى الضوضاء.

الأهمية والادعاءات
يذكر المؤلفون أنه على الرغم من أن الحدود العليا المستمدة (1–2.5 ميجابارسيك) هي أقل من المسافة الفعلية لـ M101 ($\approx$ 6.85 ميجابارسيك)، إلا أن هذا العمل يمثل أول تطبيق لمسار GFH-v2 في بحث موجه لمستعر أعظم قريب ناتج عن انهيار قلب. لقد نجح البحث في توصيف أداء المسار على بيانات كاشف حقيقية من حقبة التكليف (ER15). كما يثبت أن طريقة GFH-v2 يمكن نشرها بفعالية للبحث عن الموجات الطويلة العبور (tCWs) من الأحداث المكونة للنجوم المغناطيسية، مما يؤسس إطاراً للبحوث المستقبلية باستخدام مجموعات بيانات أكثر حساسية من دورات الرصد القادمة (O5) وكواشف الجيل التالي. ويخلص البحث إلى أن الدراسة توفر تقديراً مباشراً للحساسية للبحث عن الموجات طويلة العبور في بيانات ER15 وتعمل كمعيار للمساعي المستقبلية.