Unsupervised Morphological Characterization of Gravitational-Wave Glitches in LIGO O4a Using Frozen DINOv2 Features

ملخص تقني: التوصيف المورفولوجي غير الخاضع للإشراف لخلل موجات الجاذبية في دورة LIGO O4a باستخدام ميزات DINOv2 المجمدة

بيان المشكلة
تتعرض كواشف موجات الجاذبية الأرضية، مثل Advanced LIGO، لظواهر عابرة من الضوضاء غير الغاوسية تُعرف باسم "الخلل" (glitches). هذه العيوب يمكن أن تحاكي الإشارات الفيزيائية الفلكية، وتقلل من حساسية الكاشف، وتصعب عملية تقدير المعلمات. وبينما وضع مشروع Gravity Spy معيارًا لتصنيف الخلل باستخدام التعلم الآلي الخاضع للإشراف والمشاركة المجتمعية، فإن هذه النهج محدودة بطبيعتها في التعرف على الأشكال المورفولوجية التي تم تدريبها عليها صراحةً؛ فهي لا تستطيع تحديد أنواع الخلل التي لم يسبق توصيفها. السؤال المركزي الذي تعالجه هذه الدراسة هو ما إذا كانت دورة الرصد O4a قد قدمت أشكالاً مورفولوجية للخلل غير موجودة في الدورات السابقة، وتحديداً تلك التي قد تكون غير مرئية للمصنفات الخاضعة للإشراف المدربة على البيانات التاريخية.

المنهجية
يقدم المؤلفون gravi-signal-ml، وهو مسار عمل (pipeline) كامل غير خاضع للإشراف وصفر-القدرة (zero-shot)، مصمم لتوصيف عابرات الضوضاء الآلية دون الحاجة إلى بيانات تدريب مصنفة، أو ضبط دقيق لنطاق محدد، أو تسريع عبر وحدة معالجة الرسومات (GPU). تسير المنهجية كما يلي:

1. المعالجة المسبقة للبيانات: يقوم مسار العمل بتحليل 1,277 ساعة من بيانات الإجهاد العامة من كاشفي هانفورد (H1) وليفينغستون (L1) عبر أربع نوافذ زمنية مستقلة في عام 2023. يتم تقسيم البيانات الخام إلى قطع مدة كل منها 32 ثانية، ثم تبييضها (whitening)، وتمريرها عبر مرشح تمرير نطاق (20–2000 هرتز)، وتحويلها إلى مخططات طيفية (spectrograms) بأبعاد 256×256 بكسل باستخدام لوحة ألوان cividis لضمان التوحيد الإدراكي عند التحويل إلى التدرج الرمادي.
2. استخراج الميزات: يتم ترميز المخططات الطيفية باستخدام نموذج DINOv2 Vision Transformer (ViT-S/14 مع رموز التسجيل - register tokens) مجمد، والمدرب مسبقاً على الصور الطبيعية. يُستخدم النموذج في وضع التقييم دون أي ضبط دقيق. يتم استخراج مخرج رمز CLS كتمثيل (embedding) ذي 384 بُعداً ويتم تطبيعه باستخدام معيار L2. الهدف من استخدام رموز التسجيل هو كبح الآثار الاصطناعية الشائعة عند معالجة الصور غير الطبيعية.
3. تقليل الأبعاد: يتم إسقاط التمثيلات من 384 إلى 50 بُعداً عبر تحليل المكونات الرئيسية (PCA)، مع الاحتفاظ بأكثر من 98% من التباين. بعد ذلك، يقوم خوارزم (UMAP) بتقليل البيانات إلى 10 أبعاد للتجميع، باستخدام مقياس جيب التمام (cosine metric) لمراعاة هندسة الكرة الوحدة للتمثيلات.
4. التجميع: يتم استخدام نموذج الخليط لعملية ديريكليه (DPMM) لتحديد عدد المجموعات المورفولوجية تلقائياً. وخلافاً للطرق القائمة على الكثافة (مثل HDBSCAN)، يتجنب DPMM الحساسية تجاه الحد الأدنى لحجم الكتلة عن طريق تعيين وزن ضئيل للمكونات الزائدة.
5. التحقق والمراجعة المتقاطعة:
   • المراجعة المورفولوجية المتقاطعة: تتم مقارنة الكتل الشاذة (تلك ذات درجات الاحتمالية اللوغاريتمية المنخفضة) مع مؤشر مرجعي داخل النطاق مبني من مجموعة بيانات Gravity Spy المصنفة لـ O3b، والتي تمت معالجتها عبر نفس المسار تماماً لإزالة الفجوات النطاقية. تُصنف العينات كـ معروفة، أو غامضة، أو جديدة بناءً على عتبات تشابه جيب التمام.
   • اختبارات المتانة: يخضع مسار العمل لدراسات الاستئصال (التحويل للرمادي، عكس البكسل، تغيير قياس الكثافة)، وتحليل الاستقرار (اضطراب المعلمات الفائقة)، واختبارات خلفية إزاحة الوقت لاستبعاد أي مصادفات ذات دلالة إحصائية بين H1 وL1.

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

1. مسار عمل صفر-القدرة (Zero-Shot): أول تطبيق لميزات DINOv2 المجمدة في توصيف مورفولوجيا خلل موجات الجاذبية، مما لا يتطلب بيانات مصنفة أو موارد وحدة معالجة رسومات.
2. تحقق منهجي: إطار عمل للتحقق متعدد الطبقات مطبق باستمرار عبر أربع نوافذ مستقلة لـ O4a، بما في ذلك الاستئصال، والاستقرار، وتقدير خلفية إزاحة الوقت.
3. خط أساس قابل لإعادة الإنتاج: توفير مسار عمل موثق بالكامل ومفتوح المصدر ونتائج لتوصيف الخلل بنمط صفر-القدرة.
4. ملاحظة تجريبية: تحديد اختلاف منهجي في متانة التمثيلات بين كاشفي H1 وL1.

النتائج
عالج مسار العمل أكثر من 188,000 مخطط طيفي عبر الجلسات الأربع. أسفرت النتائج عن نتيجة سلبية (null result) فيما يتعلق بالأشكال المورفولوجية الجديدة:

• لا وجود لخلل مورفولوجي جديد: لم يتم تحديد أي مرشح لخلل مورفولوجي جديد. جميع الكتل الشاذة طابقت فئات Gravity Spy المعروفة بمتوسط تشابه جيب تمام (top-1) أكبر من 0.98.
• استقرار الكتل: حدد نموذج DPMM ما بين 10 إلى 16 مجموعة مورفولوجية لكل كاشف لكل جلسة. أكد تحليل الاستقرار قابلية إعادة الإنتاج، حيث تجاوز متوسط قيم مؤشر راند المعدل (ARI) 0.83 عبر التشغيلات المضطربة.
• عدم التماثل بين الكواشف: لوحظ عدم تماثل ملحو-ظ في متانة التمثيلات تجاه التحويل للرمادي. أظهرت تمثيلات L1 متانة عالية (ARI للرمادي > 0.90 في معظم الجلسات)، بينما أظهرت تمثيلات H1 متانة أقل وأكثر تباينًا (تراوح ARI للرمادي بين 0.62 و0.90). يشير هذا إلى اختلاف هيكلي في فضاء ضوضاء H1 تحت استخراج ميزات DINOv2، وهو ما قد يرتبط بعدم الاستقرار المعروف في أرضية ضوضاء H1.
• استبعاد الخلفية: أظهرت اختبارات إزاحة الوقت عدم وجود مصادفات ذات دلالة إحصائية بين H1 وL1 (p ≥ 0.1)، مما يستبعد الآثار الآلية المترابطة كسبب لهذه المجموعات.
• المقارنة المرجعية: عند مقارنتها بمرجع Gravity Spy، حقق نموذج DPMM مؤشر ARI قدره 0.133. يشير المؤلفون إلى أن هذه الدرجة المتوسطة تعكس التباعد الجوهري بين الاصطلاحات المورفولوجية البشرية والتشابه البصري الذاتي في فضاء التمثيل، وليس فشلاً في الطريقة.

الأهمية والادعاءات
يزعم البحث أن نتيجته السلبية ذات معنى علمي. فهي تثبت أنه عند مستوى الحساسية والدقة الزمنية لهذا المسار المحدد (نوافذ 32 ثانية، هيكل Vi-S/14)، فإن مجتمع خلل O4a متسق مورفولوجيًا مع كتالوج Gravity Spy O3b. وهذا يتماشى مع النتائج السابقة التي تشير إلى أن خصائص O4a رفعت معدلات فئات الخلل المعروفة (مثل Low Frequency Lines، 1400Ripples) بدلاً من تقديم مورفولوجيا جديدة جذريًا.

يوضح المؤلفون بتواضع أن هذه النتيجة لا تثبت عدم وجود خلل جديد في O4a؛ بل تحدد انتشارها بما يقل عن عتبة الكشف لهذا النهج. تشمل القيود المعترف بها النافذة الثابتة (32 ثانية، والتي قد تفقد الأحداث التي تقل عن ثانية واحدة)، واستخدام هيكل ViT أصغر، وحقيقة أن النوافذ التي تم تحليلها تغطي حوالي 5% فقط من مجموعة بيانات O4a الكاملة. يعمل هذا العمل كخط أساس قابل لإعادة الإنتاج لجهود التوصيف المستقبلي صفر-القدرة؛ وبشكل خاص، فإن مسار العمل مفتوح المصدر 'gravi-signal-ml' متاح على GitHub على الرابط https://github.com/lucacirfeta/dante-gravi-signal-ml، مع إصدار مؤرشف برقم DOI 10.5281/zenodo.20121860.

تعالج هذه الدراسة الحاجة المتزايدة لأدوات قابلة للتوسع وغير خاضعة للإشراف في عصر توسع فهارس موجات الجاذبية؛ وللسياق، في مايو 2026، أصدر تعاون LIGO-Virgo-KAGRA كتالوج GWTC-5.0، الذي أبلغ عن 161 حدثًا جديدًا لموجات الجاذبية، مما يعكس حجم البيانات المتزايد الذي يحفز مثل هذه النهج في التوصيف الآلي.