Interpreting the Hubble tension with a cascade decaying dark matter sector

ملخص تقني: تفسير توتر هابل باستخدام قطاع مادة مظلمة متسلسل التحلل

بيان المشكلة
تتناول الورقة "توتر هابل"، وهو تباين قدره $\sim 5.8\sigma$ بين قياسات الكون المبكر لثابت هابل ($H_0 = 67.36 \pm 0.54$ كم ث$^{-1}$ ميجا بارسيك$^{-1}$) من بيانات إشعاع الخلفية الكونية الميكروي (Planck 2018)، والقياسات المباشرة للكون المتأخر ($H_0 = 73.01 \pm 0.92$ كم ث$^{-1}$ ميجا بارسيك$^{-1}$) من تعاون SH0ES. ضمن إطار نموذج $\Lambda$CDM القياسي، يُعتبر هذا التوتر على نطاق واسع مؤشراً على فيزياء جديدة. ركزت المحاولات السابقة لحله إما على تعديلات في الزمن المبكر (مثل الطاقة المظلمة المبكرة - EDE) أو تعديلات في الزمن المتأخر. ومع ذلك، يشير المؤلفون إلى أن التعديلات ذات التأثير الواحد تفشل عموماً في تقليل التوتر إلى مستوى الـ $3\sigma$ المقبول، مع كون الطاقة المظلمة المبكرة (EDE) استثناءً ملحوظاً. تبحث الورقة فيما إذا كان بإمكان نموذج موحد يدمج كلاً من تأثيرات الزمن المبكر والزمن المتأخر تحقيق قيمة أعلى لـ $H_0$ مع البقاء متوافقاً مع القيود الكونية الأخرى.

المنهجية
يقترح المؤلفون ويحللون نموذج المادة المظلمة متسلسلة التحلل (CDDM). في هذا الإطار، يتم استبدال المادة المظلمة الباردة (CDM) القياسية بقطاع يتضمن نوعين من الجسيمات: $\chi_M$ و $\chi_m$:

1. الإنتاج في الزمن المبكر: جسيم أب ثقيل $\chi_M$ (كتلته $M$) يتحلل إلى جسيم أخف $\chi_m$ (كتلته $m$) وحالة نهائية من النموذج المعياري (SM) وهي $X$ خلال الكون المبكر ($\tau_M \le 10^4$ ثانية). جسيمات $\chi_m$ الناتجة تكون نسبية في البداية، مما يساهم في عدد النيوترينوات الفعال ($N_{\text{eff}}$)، ثم تصبح غير نسبية لاحقاً.
2. التحلل في الزمن المتأخر: جسيمات $\chi_m$، التي تعمل كمادة مظلمة، تتحلل لاحقاً إلى نيوترينوهات من النموذج المعياري ($\chi_m \to \nu + \bar{\nu}$) في الأوقات المتأخرة ($\tau_m \ge 100$ مليار سنة).

يتم تمثيل النموذج بثلاثة بارامترات مستقلة هي: نسبة الكتلة $M/m$، وعمر الجسيم الأب $\tau_M$، ومعدل تحلل الجسيم الابن $\Gamma_m$ (المرتبط بـ $\tau_m$). يقوم المؤلفون بدمج معادلات الخلفية والاضطراب الخطي لهذا القطاع في برنامج الحل البولتزماني CLASS، ويجرون تحليل مونت كارلو لسلاسل ماركوف (MCMC) باستخدام Cobaya.

يتم مطابقة النموذج مع أحدث مجموعات البيانات:

• إشعاع الخلفية الكونية (CMB): طيف درجات الحرارة والاستقطاب لـ Planck 2018 عند الزوايا المنخفضة ($\ell$)، وأطياف TT/TE/EE عند الزوايا العالية، وإمكانات العدسة الجاذبية.
• التذبذبات الصوتية لبارونات الباريون (BAO): بيانات DESI DR2 (قيود المسافة المتناحية وغير المتناحية).
• المستعرات العظمى (Supernovae): تجميعة Pantheon+.
• الأولوية المحلية (Local Prior): قياس SH0ES لـ $H_0$ ($73.04 \pm 1.04$ كم ث$^{-1}$ ميجا بارسيك$^{-1}$).

تغير الدراسة أولويات البارامترات (priors) بشكل منهجي وتدرج في تضمين أو استبعاد أولوية $H_0$ المحلية لتقييم قوة النتائج. وأخيراً، يتم اختبار مناطق البارامترات المفضلة مقابل القيود التكميلية: تخليق العناصر (BBN)، حدود تدفق النيوترينو من التلسكوبات المختلفة، وتكون البنية (طول حرية الحركة).

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

1. قيم ثابت هابل ($H_0$):

   • عند مطابقة Planck 2018 + DESI BAO + Pantheon + SH0ES مع أولويات محددة (Priors III)، يعطي النموذج $H_0 = 69.05^{+0.31}_{-0.27}$ كم ث$^{-1}$ ميجا بارسيك$^{-1}$ (عند مستوى ثقة 68%).
   • بدون أولوية SH0ES المحلية، تنخفض القيمة إلى $H_0 = 68.76 \pm 0.35$ كم ث$^{-1}$ ميجا بارسيك$^{-1}$.
   • تقلل هذه القيم من التوتر إلى حوالي $3.8\sigma$.

2. تأثير الأولويات و $\Delta\chi^2$:

   • تكشف الدراسة عن اعتماد قوي على أولويات البارامترات. وبينما يمكن لتعديل الأولويات أن يؤدي إلى قيم أعلى قليلاً لـ $H_0$، إلا أن هذا يأتي على حساب زيادة كبيرة في $\Delta\chi^2$ (تدهور المطابقة الإحصائية).
   • بالنسبة لأفضل سيناريو مطابقة (Priors III)، فإن $\Delta\chi^2 = +16.0$ بالنسبة لنموذج $\Lambda$CDM، مما يشير إلى تفضيل إحصائي للنموذج القياسي على نموذج CDDM بالنظر إلى مجموعات البيانات المستخدمة.
   • يلاحظ المؤلفون اتجاهاً مفاده أن تقليل التوتر بشكل أكبر (أقل من $3\sigma$) يتطلب أولويات تؤدي إلى $\Delta\chi^2$ كبير للغاية.

3. القيود التكميلية:

   • تخليق العناصر (BBN): مناطق البارامترات التي تفضلها البيانات الكونية (تحديداً $\tau_M$ و $M/m$) تتوافق مع حدود BBN على وفرة العناصر الخفيفة (D/H و $^7$Li/H).
   • تدفق النيوترينو: التحلل المتأخر إلى نيوترينوهات يعني أن نطاق كتلة $\chi_m$ هو $m \sim 1-10$ MeV، وهو ما يظل متوافقاً مع قيود Borexino و KamLAND و Super-Kamiokande.
   • تكون البنية: تم حساب طول حرية الحركة ($\lambda_{fs}$) الناجم عن سرعة الركلة الناتجة عن التحلل المبكر ليكون $\sim 10^{-2}$ ميجا بارسيك، وهو ضمن الحد الرصدي $\lambda_{fs} < 0.1$ ميجا بارسيك.

الأهمية والادعاءات
تدعي الورقة أنها تصحح نتائج الأدبيات السابقة التي أشارت إلى إمكانية تقليل توتر هابل إلى ما دون مستوى $3\sigma$ ضمن سيناريوهات تعديل المادة المظلمة. يجادل المؤلفون بأن الدراسات السابقة ربما عانت من تحليلات MCMC غير دقيقة أو غير مكتملة.

الاستنتاج الأساسي هو أنه بينما ينجح نموذج CDDM في دمج كل من تأثيرات الزمن المبكر (النسبية) والزمن المتأخر (التحلل) لرفع قيمة $H_0$، فإنه لا يستطيع تقليل توتر هابل إلى ما دون مستوى $\sim 3\sigma$ دون تكبد عقوبة إحصائية كبيرة ( $\Delta\chi^2$ كبير). ويؤكد المؤلفون أن التوتر يظل تحدياً كبيراً، حيث إن النموذج لا يقدم حلاً "مجانيًا" يلبي جميع مجموعات البيانات والمعايير الإحصائية في آن واحد. يعمل هذا العمل كاختبار صارم لنموذج المادة المظلمة متسلسلة التحلل، موضحاً توافقه مع الحدود الحالية مع تسليط الضوء على صعوبة حل التوتر من خلال هذه الآلية وحدها.