Vertex reconstruction for a search for neutron-antineutron conversions with HIBEAM
تقيم هذه الورقة وتقارن بين الطرق الكلاسيكية والطرق القائمة على الشبكات العصبية الرسومية لإعادة بناء رؤوس الفناء في تجربة HIBEAM، وتجد أنه بينما يحقق كلا النهجين دقة مماثلة لإحداثيات الرأس، فإن التعلم الآلي يقدم مزايا فريدة في استخراج معلومات شكل الحدث للتحليل اللاحق.
المؤلفون الأصليون:Alexander Burgman, Sze Chun Yiu, Yamna Shaikh, David Milstead, Eily Merjy, Lucas Åstrand, Kenneth Österberg, Fredrik Oljemark, Matthias Holl, Valentina Santoro, André Nepomuceno, Joshua L. Barrow
المؤلفون الأصليون: Alexander Burgman, Sze Chun Yiu, Yamna Shaikh, David Milstead, Eily Merjy, Lucas Åstrand, Kenneth Österberg, Fredrik Oljemark, Matthias Holl, Valentina Santoro, André Nepomuceno, Joshua L. Barrow
الكون مليء بالمادة، من النجوم في الأعالي إلى الذرات في أجسادنا، ومع ذلك يظل لغزاً عميقاً سبب وجود هذا القدر الهائل منها ومن القليل جداً من نقيضها، المادة المضادة. وفقاً لقوانيع الفيزياء كما نفهمها حالياً، كان ينبغي للانفجار العظيم أن يخلق كميات متساوية من كليهما، مما يؤدي إلى فنائهما فوراً، ولا يترك خلفهما شيئاً. ولتفسير سبب وجودنا، يبحث العلماء عن انتهاكات نادرة لقاعدة أساسية تسمى "حفظ العدد الباريوني". تنص هذه القاعدة على أن العدد الإجمالي للجسيمات الثقيلة، مثل البروتونات والنيوترونات، يجب أن يظل ثابتاً. فإذا استطاع نيوترون أن يتحول تلقائياً إلى نيوترون مضاد، فإن ذلك سيخرق هذه القاعدة ويقدم دليلاً على عدم التوازن الكوني. إن مثل هذا التحول يصعب رصده للغاية لأنه لم يُرَ من قبل، كما أن الإشارات التي قد ينتجها يسهل طمسها بسبب الضجيج الخلفي.
وللبحث عن هذا الحدث المراوغ، يجري التخطيط لتجربة من جيل جديد في "مصدر السبليشن الأوروبي" (European Spallation Source) في السويد. وتتضمن الخطة إرسال حزمة من النيوترونات الباردة عبر أنبوب محمي وإلى رقاقة معدنية رقيقة. إذا تحول نيوترون إلى نيوترون مضاد، فسيصطدم فوراً بالرقاقة وينفجر إلى وابل من الجسيمات المشحونة. والمفتاح لتحديد هذا الانفجار هو تحديد مكان حدوثه بدقة على الرقاقة. تستخدم التجربة كاشفاً كبيراً مملوءاً بالغاز، يسمى "غرفة إسقاط زمن التصوير" (time projection chamber)، والذي يسجل مسارات الجسيمات أثناء تطايرها. ومع ذلك، يرى الكاشف أيضاً جسيمات أخرى، مثل الإلكترونات المتشتتة بفعل الإشعاع الشارد، والتي يمكن أن تبدو مشابهة جداً للإشارة الأصلية. ويتمثل التحدي في فرز الانفجار الحقيقي عن الضجيج وتحديد مركز الحدث بدقة متناهية.
قام فريق من الباحثين مؤخراً باختبار مدى كفاءة طرق حاسوبية مختلفة في حل هذا اللغز باستخدام عمليات محاكاة تفصيلية لكاشف (HIBEAM). لقد أنشأوا الملايين من الأحداث الافتراضية حيث يفني نيوترون مضاد نفسه على الرقاقة، منتجاً بضعة بيونات مشحونة تسافر إلى الكاشف. ولجعل الاختبار واقعياً، قاموا بحقن مسارات خلفية إضافية، لمحاكاة الإلكترونات العشوائية التي سيراها الكاشف في تشغيل حقيقي. ثم قام الباحثون بتشغيل أربع استراتيجيات إعادة بناء مختلفة على هذه البيانات لمعرفة أي منها يمكنه الأفضل في حساب موقع الانفجار. استخدمت إحدى الطرق نهجاً إحصائياً كلاسيكياً استُخدم في الفيزياء لعقود، بينما اعتمدت طريقة أخرى على خدعة هندسية بحتة تتجاهل مفهوم المسارات تماماً. واستخدمت طريقتان أخريان تعلم الآلة الحديث، وتحديداً نوعاً من الذكاء الاصطناعي يسمى "الشبكة العصبية الرسومية" (graph neural network)، والتي تتعامل مع ضربات الكاشف كشبكة معلومات مترابطة لإيجاد الأنماط.
أظهرت النتائج أنه بالنسبة للهندسة البسيطة لهذا الكاشف المحدد، أدت الطرق الإحصائية التقليدية أداءً يضاهي أدوات تعلم الآلة الجديدة عندما يتعلق الأمر بإيجاد الإحداثيات الدقيقة للحدث. وفي غياب الضجيج الخلفي، تمكنت جميع الطرق من تحديد مركز الانفجار بدقة تصل إلى بضعة مليمترات. ومع ذلك، تصرفت الطرق بشكل مختلف عندما زاد الضجيج الخلفي؛ فقد كانت الطريقة الهندسية البحتة ونموذج تعلم الآلة الذي ينظر إلى الحدث ككل مستقرين للغاية، حيث لم تتغير دقتها إلا قليلاً حتى عند مضاعفة أو أربعة أضعاف عدد الإلكترونات الخلفية. وفي المقابل، بدأت الطريقة الإحصائية الكلاسيكية وسلسلة تعلم الآلة الهجينة الأكثر تعقيداً في المعاناة مع نمو الضجيج، مما أدى إلى إنتاج بعض التخمينات السيئة للغاية التي وسعت نطاق الخطأ بشكل كبير.
وعلى الرغم من الأداء المتشابه في إيجاد المركز، قدمت طرق تعلم الآلة شيئاً لم تستطع الطرق الأخرى تقديمه؛ فقد وفرت معلومات إضافية حول شكل الحدث، مثل عدد الجسيمات التي يُحتمل وجودها ودرجة الثقة لكل تخمين. هذه البيانات الإضافية قيمة لأنها تساعد العلماء لاحقاً في تقرير ما إذا كان الحدث اكتشافاً حقيقياً أم مجرد صدفة عابرة. وخلصت الدراسة إلى أنه بينما تعتبر الطرق الكلاسيكية كافية للمهمة الأساسية المتمثلة في إيجاد الرأس (vertex)، فإن مناهج تعلم الآلة تقدم صورة أكثر ثراءً للحدث. كما وجد الباحثون أنه إذا خففوا القواعد لقبول الأحداث التي تحتوي على مسار مرئي واحد فقط، فيمكنهم رصد كل حدث تقريباً، ولكن موقع الانفجار سيصبح أقل يقيناً بكثير. وفي النهاية، يشير العمل إلى أن المسار الأفضل لهذا البحث تحديداً قد يكون الجمع بين موثوقية الطرق التقليدية والقدرة الوصفية لتعلم الآلة لضمان أنه إذا تحول نيوترون إلى نيوترون مضاد، فإن التجربة لن تجده فحسب، بل ستفهمه تماماً.
ملخص تقني: إعادة بناء الرأس (Vertex Reconstruction) لتحولات النيوترون-أنتينيوترون باستخدام HIBEAM
بيان المشكلة يهدف برنامج HIBEAM/NNBAR في مصدر السبوندات الأوروبي (ESS) إلى البحث عن انتهاك عدد الباريون عبر تحولات النيوترون-أنتينيوترون (n→nˉ). وتعد إعادة بناء الرأس (vertex) حيث يفنى الأنتينيوترون على نوية داخل رقاقة هدف كربونية رقيقة، معلماً مرصوداً حاسماً لتحديد حدث الإشارة. ينتج عن هذا الفناء حالة نهائية من 2 إلى 4 بيونات مشحونة تدخل إلى غرفة إسقاط زمن (TPC). ويكمن التحدي في إعادة بناء إحداثيات الرأس المستعرضة (x,y) بدقة على مستوى رقاقة الهدف (z=0) من قطع مسارات قصيرة ومستقيمة تقريباً. يجب أن تظل عملية إعادة البناء قوية ضد المسارات الخلفية، وتحديداً إلكترونات كومبتون المشتتة منخفضة الطاقة التي تدخل من المنطقة الداخلية لـ TPC، والتي يمكن أن تحاكي طوبولوجيا الإشارة أو تعقد عملية تجميع المسارات (clustering). حققت التجارب السابقة (مثل ILL) دقة شعاعية مستعرضة بلغت حوالي 4 سم؛ وتسعى هذه الدراسة إلى قياس أداء تقنيات إعادة البناء الحديثة لتحديد ما إذا كان بإمكان HIBEAM تحسين هذه الدقة بشكل كبير.
المنهجية تستخدم الدراسة محاكاة للكاشف تعتمد على هندسة HIBEAM، وتتكون من رقاقة هدف كربونية بنصف قطر 20 سم محاطة بـ TPC مملوء بغاز Ar/CO2. تم استخدام نماذج محاكاة لأحداث الإشارة باستخدام نموذج فناء نووي محيطي، مما أنتج 106 حدثاً موزعة بانتظام فوق الرقاقة. ولاختبار الخوارزميات تحت الضغط، تم إنشاء عينات خلفية عن طريق حقن 0، 1، 2، 4، أو 8 مسارات من إلكترونات كومبتون لكل حدث.
تم تقييم أربع خوارزميات متميزة لإعادة بناء الرأس:
Kalman: نهج كلاسيكي يستخدم DBSCAN لتجميع النقاط (hits)، ومرشح كالمان (Kalman filter) لملاءمة المسار (مع أربعة بارامترات للحالة)، وملاءمة رأس تكيفية مقيدة بمستوى الرقاقة. تستخدم هذه الطريقة التخمير الحتمي (deterministic annealing) وتزن المسارات بناءً على امتدادها إلى الرقاقة.
TrackLess: طريقة هندسية بحتة تتجنب ملاءمة المسار الصريحة. تقوم بإسقاط أزواج النقاط غير المرتبة على مستوى الرقاقة، ثم تجمع الإسقاطات باستخدام DBSCAN، وتعرف الرأس كمركز ثقل لأكبر تجمع. وتتضمن قطعاً للاتساق الزاوي لرفض أحداث المسار الواحد.
GraphNeT: تنفيذ لإطار عمل الشبكات العصبية الرسومية (GNN) مفتوح المص المصدر، والذي تم تطويره في الأصل لـ IceCube. تتعامل هذه الطريقة مع نقاط الكاشف كعُقد في رسم بياني متصلة عبر أقرب k من الجيران. وباستخدام بنية DynEdge، تتعلم الشبكة الارتباطات المحلية والعالمية عبر تمرير الرسائل (message passing) لتقدير إحداثيات الرأس وعدد الجسيمات المشحونة.
IterGNN: خط معالجة هجين يجمع بين الخطوات الهندسية الحتمية ومرحلتين من GNN. يبدأ بتجميع النقاط وملاءمة المسارات المستقيمة، يليه مرحلة TrackGNN لتصنيف المسارات كإشارة أو خلفية. بعد ذلك، يتم حساب بذرة (seed) لسطح الرقاقة موزونة باستخدام هذه التصنيفات، والتي يتم صقلها بواسطة مرحلة VertexGNN التي تتنبأ بالتصحيحات (Δx,Δy) وتقديرات عدم اليقين لكل حدث.
المساهمات الرئيسية
التحليل المقارن: تقدم الورقة مقارنة منهجية بين الطرق الإحصائية الكلاسيكية، والتقنيات الهندسية البحتة، ونهج التعلم الآلي (ML) لطوبولوجيا ذات تعددية منخفضة (HIBEAM TPC).
البنية الهجينة: يوضح تسلسل IterGNN طريقة تستفيد من القيود الهندسية الصريحة (مستوى الرقاقة، اتجاه المسار) مع استخدام GNN للجوانب التوافقية لتصنيف المسارات وصقل الرأس.
اختبار الضغط بالخلفية: تقيم الدراسة بدقة أداء الخوارزميات مقابل زيادة كثافة خلفيات إلكترونات كومبتون المحقونة، وهو عامل حاسم في عمليات البحث القائمة على انعدام الخلفية.
المعلومات الإضافية: تسلط الورقة الضوء على أن طرق التعلم الآلي (GraphNeT و IterGNN) توفر مخرجات إضافية لا تتوفر في الطرق الكلاسيكية، مثل احتمالات مسار الإشارة، وعدد المسارات لكل حدث، وتقديرات عدم اليقين للرأس لكل حدث.
النتائج
الدقة: في غياب الخلفية (0 إلكترون كومبتون)، تحقق جميع الطرق الأربع خطأ نصف قطر متوسط (Δr) في نطاق بضعة مليمترات (4.3 ملم إلى 6.1 ملم).
تُظهر طريقتي TrackLess و GraphNeT أدق توزيعات للخطأ وأصغر أنصاف أقطار احتواء بنسبة 90% (12–15 ملم)، مع تدهور ضئيل مع زيادة تعدد الخلفية.
تُظهر طريقة Kalman ذيلاً كبيراً غير غاوسي ينمو بشكل ملحوظ مع الخلفية (يصل نصف قطر الاحتواء بنسبة 90% إلى 244 ملم عند وجود 8 إلكترونات كومبتون).
تؤدي IterGNN أداءً متوسطاً، بمتوسط خطأ قدره 5.6 ملم (عند 0 إلكترون كومبتون) واحتواء بنسبة 90% يصل إلى 28 ملم، والذي ينمو إلى 116 ملم عند 8 إلكترونات كومبتون.
الكفاءة: تحقق الطرق الكلاسيكية والهجينة (Kalman, IterGNN) عموماً كفاءات إعادة بناء أعلى (∼80–94%) مقار بـ TrackLess (∼68–70%) في وجود الخلفية. وتزديد كفاءة IterGNN و Kalman مع تعدد الخلفية، ويرجع ذلك على الأرجح إلى تضمين مسارات الخلفية في عدد المسارات.
تقدير عدم اليقين: نجحت IterGNN في التنبؤ بعدم اليقين المستعرض لكل حدث (σx,σy)، والتي ترتبط بعدد النقاط وكثافة الخلفية.
أحداث المسار الواحد: ترفض جميع الطرق مبدئياً الأحداث التي تحتوي على مسار واحد فقط لضمان موثوقية الرأس. يؤدي إزالة هذا القيد إلى رفع الكفاءة إلى ∼100% لجميع الطرق، ولكنه يضعف الدقة بشكل كبير، خاصة في الذيول.
الأهمية والادعاءات تخلص الورقة إلى أنه بالنسبة للأحداث بسيطة الهندسة التي تميز HIBEAM TPC، فإن طرق إعادة البناء الكلاسيكية تؤدي أداءً يضاهي طرق التعلم الآلي فيما يتعلق بدقة إعادة بناء إحداثيات الرأس. وتحديداً، في الظروف الخالية من الخلفية، يكون الأداء المتوسط لجميع الطرق متقارباً، مما يشير إلى أن القيد الأساسي هو التشتت المتعدد لجسيمات الإشارة في مادة أنبوب الشعاع وليس خوارزمية إعادة البناء نفسها.
ومع ذلك، يزعم المؤلفون أن طرق التعلم الآلي المخصصة توفر مزايا متميزة تتجاوز مجرد إعادة بناء الإحداثيات. فهي توفر معلومات حول شكل الحدث (تصنيف المسار، عد المسارات، وتقديرات عدم اليقين) وهي معلومات لا يمكن الوصول إليها عبر الطرق الهندسية أو الكلاسيكية المنفذة. وتعتبر هذه الكميات الإضافية مهمة للتحليلات اللاحقة. تشير الدراسة إلى أن العمل المستقبلي يجب أن يركز على الجمع بين هذه الطرق لتحقيق كل من دقة إحداثيات الرأس العالية وإعادة بناء شكل الحدث الموثوق، بدلاً من الاعتماد على نهج واحد. لا تدعي الورقة أن طرق التعلم الآلي تحل مشكلة إعادة بناء الرأس بشكل أفضل بطبيعتها في هذه الطوبولوجيا المحددة، لكنها تقدم معلومات تكميلية قيمة لاستراتيجية البحث الأوسع.