Search for Annual Modulation in Combined ANAIS 112 and COSINE 100 Data using Bayesian Model Comparison

ملخص تقني: مقارنة النماذج البايزية لبيانات ANAIS-112 و COSINE-100 المجمعة

بيان المشكلة
على مدار أكثر من عقدين، أبلغت تجربة DAMA/LIBRA عن إشارة تعديل سنوي ذات دلالة إحصائية في معدلات الأحداث المتبقية لكواشف يوديد الصوديوم (NaI)، وهي إشارة تتوافق مع تفاعلات المادة المظلمة من نوع الـ WIMPs. ومع ذلك، لا تزال هذه المزاعم غير مؤكدة من قبل تجارب الكشف المباشر الأخرى التي تستخدم نفس المادة المستهدفة، وتحديداً COSINE-100 (كوريا) و ANAIS-112 (إسبانيا). وبينما استخدمت التحليلات المشتركة السابقة لأول ثلاث سنوات من البيانات المجمعة من COSINE-100 و ANAIS-112 (المشار إليها بـ C25) الانحدار التكراري (Frequentist regression) لإيجاد عدم وجود دليل على التعديل، تسعى هذه الدراسة إلى استكمال تلك النتائج من خلال إطار عمل بايزي (Bayesian framework). وتتمثل المشكلة الجوهرية التي تعالجها الدراسة في ما إذا كانت المتبقيات المجمعة من هاتين التجربتين توفر دليلاً إحصائياً يرجح نموذج التعديل الجيبي (المتوقع من المادة المظلمة) على نموذج الخلفية الثابتة، بشكل مستقل عن تقنيات الانحدار المحددة المستخدمة في العمل السابق.

المنهجية
يقوم المؤلفون بإجراء مقارنة نماذج بايزية بين فرضيتين:

• النموذج $M_2$ (إشارة جيبية): تتبع المتبقيات في البيانات دالة جيبية $R(t) = S_m \cos(\omega(t - t_0))$، حيث $S_m$ هو السعة، و $\omega$ هو التردد الزاوي، و $t_0$ هو الطور.
• النموذج $M_1$ (خلفية ثابتة): توصف المتبقيات بقيمة ثابتة $A$.

يستخدم التحليل ما يقرب من ثلاث سنوات من البيانات المجمعة من COSINE-100 و ANAIS-112، والتي تغطي تعرضاً إجمالياً قدره 485 كجم-سنة. تم تقسيم البيانات إلى فترات زمنية مدتها 15 يوماً عبر نطاقين من الطاقة: 1–6 كيلو إلكترون فولت و 2–6 كيلو إلكترون فولت.

الخطوات المنهجية الرئيسية تشمل:

1. بناء دالة الإمكان (Likelihood Construction): يُفترض وجود دالة إمكان غاوسية للمتبقيات المرصودة، مع معاملة الأخطاء كأخطاء غاوسية.
2. اختيار التوزيعات المسبقة (Prior Selection): لضمان المتانة، يختبر المؤلفون أربع مجموعات متميزة من التوزيعات المسبقة للمعلمات الحرة ($S_m$، $\omega$، $t_0$):
   • السعة ($S_m$): توزيعات موحدة تتراوح من الحد الأدنى إلى الحد الأقصى للمتبقيات المرصودة، وبشكل منفصل من الصفر إلى الحد الأقصى للمتبقي المطلق.
   • التردد ($\omega$) والطور ($t_0$): كلاهما توزيعات موحدة (تغطي كامل مدة مجموعة البيانات و 0–365 يوماً على التوالي)، وتوزيعات طبيعية متمركزة حول أفضل القيم الملائمة التي أبلغت عنها DAMA/LIBRA ($T \approx 1$ سنة، $t_0 \approx 145$ يوماً).
3. حساب الدليل (Evidence Calculation): يتم حساب الدليل البايزي (الاحتمال الهامشي) لكلا النموذجين باستخدام خوارزمية العينات المتداخلة (Nested Sampling) المنفذة في Dynesty. معيار المقارنة هو اللوغاريتم الطبيعي لعامل بايز، $\ln(B_{21})$، حيث $B_{21}$ هو نسبة دليل النموذج الجيبي إلى النموذج الثابت.

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

• التحقق المستقل: يوفر هذا العمل تحققاً مستقلاً لنتائج C25 باستخدام نهج مقارنة النماذج البايزية بدلاً من انحدار المربعات الصغرى وسلاسل ماركوف مونت كارلو (MCMC) المستخدم في التحليل المشترك الأصلي.
• متانة التوزيعات المسبقة: من خلال اختبار تكوينات متعددة للتوزيعات المسبقة — بما في ذلك تلك المستندة إلى القيم المثلى لـ DAMA/LIBRA — تُظهر الدراسة أن الاستنتاج المتعلق بغياب الإشارة متين ضد التغيرات في افتراضات التوزيعات المسبقة.
• العلم المفتوح: جعل المؤلفون أكواد التحليل الخاصة بهم متاحة علناً لضمان شفافية وقابلية إعادة إنتاج حساب الدليل البايزي.

النتائج
يُسفر التحليل عن اللوغاريتم الطبيعي لعامل بايز ($\ln(B_{21})$) لجميع مجموعات التوزيعات المسبقة الأربع عبر كلا الفاصلين الطاقيين. وتُلخص النتائج كما يلي:

• بالنسبة لجميع مجموعات التوزيعات المسبقة المختبرة، تكون قيمة $\ln(B_{21})$ أقل من 1.15.
• وفقاً لمقياس جيفريز لتفسير عوامل بايز، فإن قيمة $\ln(B_{21}) < 1.15$ (ما يقابل $B_{21} \leq \sqrt{10}$) تشير إلى أن الدليل لصالح النموذج الجيبي على النموذج الثابت "بالكاد يستحق الذكر".
• على وجه التحديد، في الفترات 1–6 كيلو إلكترون فولت و 2–6 كيلو إلكترون فولت، تتراوح القيم المحسوبة تقريباً بين -0.13 و 1.36، مع تمركز معظمها بالقرب من أو تحت عتبة 1.15.
• تُظهر مخططات التوزيع اللاحق (Posterior contours) المجمعة للسعة والتردد والطور عدم وجود تفضيل قوي للمعلمات المفضلة من DAMA على فرضية العدم المتمثلة في المعدل الثابت.

الأهمية والادعاءات
تخلص الورقة إلى أن مقارنة النماذج البايزية لا تقدم أي دليل لدعم فرضية أن البيانات المجمعة لـ ANAIS-112 و COSINE-100 تحتوي على إشارة تعديل سنوي متوافقة مع تفاعلات المادة المظلمة. والنتائج متوافقة مع فرضية الخلفية الثابتة.

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