BLAST: Beyond Limber Angular power Spectra Toolkit. A fast and efficient algorithm for 3x2 pt analysis

ملخص تقني: Blast.jl – خوارزمية سريعة وفعالة لتحليل 3 × 2 pt

بيان المشكلة
يعد بزوغ فجر المسوحات الضوئية والطيفية من الجيل التالي (مثل Euclid، وLSST، وRoman، وDESI) بتقديم قياسات عالية الدقة للهيكل واسع النطاق (LSS). وتعتمد هذه المسوحات بشكل كبير على إحصائيات "3 × 2 pt"، والتي تجمع بين تكتل المجرات، وارتباطات قص (shear-shear) للعدسة الضعيفة، وارتباطات مجرة-مجرة للعدسة (galaxy-galaxy lensing). يتطلب النمذجة النظرية لهذه الإحصائيات حساب أطياف القدرة الزاوية ($C_\ell$) عبر تكامل حاصل ضرب طيف القدرة المادي ثلاثي الأبعاد ودوال بسل الكروية شديدة التذبذب عبر ثلاثة أبعاد.

النهج القياسي لتقليل التكلفة الحسابية هو تقريب "ليمبر" (Limber approximation)، الذي يبسط التكامل ثلاثي الأبعاد إلى تكامل أحادي البعد بافتراض أن الأنماط عند حقب زمنية مختلفة لا تساهم في التكامل على طول خط البصر. ومع ذلك، فإن هذا التقريب يؤدي إلى أخطاء كبيرة، لا سيما في المقاييس الزاوية الكبيرة (low $\ell$). ومن أجل استغلال القدرة التقييدية الكاملة للمسوحات المستقبلية والتحكم في الأخطاء المنهجية، تعد الحسابات غير الليمبرية (non-Limber) الدقيقة ضرورية. ومع ذلك، فإن إجراء هذه التكاملات ثلاثية الأبعاد الكاملة يعد أمراً مستحيلاً من الناحية الحسابية لطرق الاستدلال البارامتري مثل "ماركوف تشين مونتي كارلو" (MCMC)، والتي تتطلب ملايين التقييمات. على سبيل المثال، تطلبت التحليلات الأخيرة التي استخدمت تقريب ليمبر 0.5 مليون ساعة من المعالجة المركزية (CPU) للوصول إلى التقارب.

المنهجية
يقدم المؤلفون Blast.jl، وهي خوارزمية مبتكرة مصممة لحساب أطياف القدرة الزاوية الكاملة غير الليمبرية بكفاءة دون افتراض معدل نمو يعتمد على المقياس. تعتمد المنهجية الأساسية على تفكيك طيف القدرة المادي باستخدام متعددات حدود تشيبيشيف (Chebyshev polynomials).

1. تفكيك تشيبيشيف: تقوم الخوارزمية بتقريب طيف القدرة المادي السلس $P(k, \chi_1, \chi_2)$ على الفترة $[k_{min}, k_{max}]$ باستخدام توسيع سلسلة تشيبيشيف. يسمح هذا بفصل التكامل إلى مكونات هندسية ومكونات كونية.
2. الحساب المسبق للمصطلحات الهندسية: تمت إعادة هيكلة التكامل بحيث يصبح الجزء الأكثر استهلاكاً للوقت حسابياً — وهو تكامل حاصل ضرب دالتي بسل كرويتين ومتعددات حدود تشيبيشيف — مستقلاً عن المعلمات الكونية. يتم حساب هذه المصطلحات ($\tilde{T}^{AB}_{n; \ell}$) وتخزينها مرة واحدة.
3. التقييم الفعال: لأي نموذج كوني جديد، تحتاج الخوارزمية فقط إلى إعادة حساب معاملات تشيبيشيف لطيف القدرة، وهي عملية غير مكلفة حسابياً. أما التكاملات الخارجية المتبقية على المسافة المادية ($\chi$) والنسبة $R = \chi_2/\chi_1$ فتتم باستخدام قواعد تربيع قياسية (قاعدة سمبسون لـ $\chi$ وقاعدة كلينشو-كورتيس لـ $R$).
4. تحويل الإحداثيات: تم استخدام تغيير في المتغيرات من $(\chi_1, \chi_2)$ إلى $(\chi, R)$ لتحسين أخذ العينات في المنطقة حيث $\chi_1 \approx \chi_2$، والتي تساهم بشكل كبير في التكامل.
5. التصحيح غير الخطي: تعالج الطريقة أطياف القدرة الخطية وغير الخطية بشكل منفصل. حيث يتم تصحيح المساهمة غير الخطية — التي تكون كبيرة فقط في المقاييس الصغيرة حيث يكون تقريب ليمبر صالحاً — باستخدام نهج قائم على ليمبر، بينما يُطبق التكامل غير الليمبري الكامل على المكون الخطي.

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

• الابتكار الخوارزمي: تقدم Blast.jl نهجاً قائماً على التفكيك باستخدام متعددات حدود تشيبيشيف، مما يتجنب قيود التباعد الشبكي اللوغاريتمي المستخدمة في طرق FFTLog (مثل FKEM وmatter) والحل التفاضلي التكراري لطريقة ليفين (Levin).
• استراتيجية الحساب المسبق: من خلال عزل التكاملات الهندسية المستقلة عن الكونيات، تقلل الخوارژية بشكل كبير من الوقت المطلوب للتقييمات المتكررة أثناء الاستدلال البارامتري.
• المرونة: على عكس بعض الطرق المنافسة، لا تفترض Blast.jl معدل نمو مستقل عن المقياس، مما يجعلها قابلاً للتطبيق على نماذج الكون التي تحتوي على نيوترينوات ضخمة أو جاذبية معدلة.
• المصدر المفتوح: الكود متاح علناً على GitHub، مرفقاً بمستودع يفصل كيفية استخدامه.

النتائج
تم تقييم الخوارزمية مقابل تحدي "N5K" (Non-local No-Nonsense Non-Limber Numerical Knock-out)، والذي نظمه تعاون علم طاقة الظلام لـ LSST لتقييم طرق التكامل غير الليمبرية. أُجريت الاختبارات على عقدة حوسبة عالية الأداء (Narval) باستخدام 64 مساراً (threads)، مع مقارنة Blast.jl مع الفائزين في التحدي: FKEM (CosmoLike)، وmatter، وLevin.

• السرعة: عند متطلبات الدقة النموذجية ($\Delta\chi^2 < 0.2$ لـ $\ell < 200$)، حققت Blast.jl وقت تشغيل قدره 0.00512 ثانية باستخدام 64 نواة. وهذا أسرع بنحو 10 مرات من FKEM (أسرع طريقة سابقة، 0.05871 ثانية) وأسرع بنحو 15 مرة عند تحسينها لتعمل بـ 32 نواة (0.00360 ثانية).
• التوسع (Scaling): تُظهر Blast.jl توسعاً ممتازاً فيما يتعلق بعدد صناديق الإزاحة الحمراء التوماتية (tomographic redshift bins). وبينما تُظهر الطرق الأخرى زيادة في وقت التشغيل مع زيادة عدد الصناديق، يظل وقت تشغيل Blast.jl مستقلاً فعلياً عن عدد الصناديق (يتراوح من 2 إلى 15).
• الدقة: تلبي الطريقة متطلبات الدقة الخاصة بالتحدي. ويشير المؤلفون إلى أنه يمكن ضبط الدقة عن طريق تعديل عدد متعددات حدود تشيبيشيف ونقاط أخذ العينات دون التأثير بشكل كبير على وقت التشغيل.
• المفاضلة بين الدقة والسرعة: أدى تقليل دقة الفاصلة العائمة من Float64 إلى Float32 إلى تحسين وقت التشغيل بمعامل قدره 2 مع تأثير ضئيل على الدقة ($\Delta\chi^2$ ارتفعت من 0.1054 إلى 0.1064).

الأهمية
يزعم البحث أن Blast.jl تمثل أسرع خوارزمية تكامل غير ليمبرية معروفة متاحة حالياً، حيث تتفوق على الطرق المتطورة بأكثر من رتبة مقدار واحدة في السيناريو النموذجي. تكمن أهميتها الأساسية في جعل التحليلات الكونية باستخدام إحصائيات 3 × 2 pt ممكنة حسابياً للمسوحات من الجيل القادم دون الاعتماد على تقريب ليمبر المثير للجدل.

يؤكد المؤلفون أن هذه السرعة أمر بالغ الأهمية لتمكين الاستدلال البارامتري القائم على MCMC مع مجموعات البيانات الضخمة المتوقعة من المسوحات المستقبلية. علاوة على ذلك، يسمح تصميم الخوارزمية بتوسعات مستقبلية، مثل تضمين تشوهات الفضاء الأحمر (redshift space distortions)، وانحياز التكبير (magnification bias)، والتصحيحات النسبوية، بالإضافة إلى دمج عدسة الخلفية الميكروية للكون (CMB lensing). كما يسلط المؤلفون الضوء على الخطط لإدراج الاشتقاق التلقائي (automatic differentiation)، مما سيمكن من استخدام طرق الاستدلال القائمة على التدرج، وبالتالي تسريع عمليات تحليل الكونيات.