Euclid preparation: The flat-sky approximation for the clustering of Euclid's photometric galaxies

1. بيان المشكلة

تهدف مهمة إقليدس (Euclid) إلى رسم خريطة للبنية واسعة النطاق للكون باستخدام مسح ضوئي للمجرات يضم حوالي $1.5 \times 10^9$ مجرة. ويتمثل أحد الأهداف الرئيسية في قياس أطياف القدرة الزاوية ($C_\ell$) لتعداد المجرات لاستنتاج المعلمات الكونية عبر تحليل سلسلة ماركوف مونت كارلو (MCMC).

• التحدي: حساب طيف القدرة الزاوي الكامل نظرياً (باستخدام أكواد مثل CLASS أو CAMB) مكلف حاسوبياً، مما يجعل من غير الممكن إجراء آلاف التكرارات المطلوبة لعملية MCMC.
• المعيار الحالي: يُستخدم تقريب ليمبر (Limber approximation) على نطاق واسع لتسريع الحسابات (بمقدار $\sim$100 ضعف) عن طريق استبدال التكامل على طول خط البصر بدالة دلتا ديراك، ومع ذلك فمن المعروف أنه يفشل في حالة نطاقات الإزاحة الحمراء الضيقة وعند المقاييس الكبيرة (قيم $\ell$ المنخفضة، $\ell \lesssim 100$)، مما يؤدي غالباً إلى أخطاء تصل إلى عامل 2 أو أكثر.
• السؤال: هل تقريب ليمبر دقيق بما يكفي لبيانات إصدارات إقليدس المحددة (DR1 و DR3)، أم أنه يتسبب في تحيزات تتجاوز الشكوك الرصدية؟ وإذا لم يكن كذلك، فهل يمكن لـ تقريب السماء المسطحة (flat-sky approximation) أن يوفر بديلاً أسرع وقابلاً للتطبيق؟

2. المنهجية

قارن المؤلفون بين ثلاث طرق لحساب أطياف القدرة الزاوية لتعداد المجرات:

1. الحساب الكامل: الحل الدقيق باستخدام حل بولتزمان CLASS، بما في ذلك تقلبات الكثافة، وتشوهات السرعة في الفضاء الأحمر (RSD)، والتأثيرات النسبوية (تكبير العدسة).
2. تقريب ليمبر: التقريب القياسي المستخدم في العديد من خطوط المعالجة (مثل CLOE).
3. تقريب السماء المسطحة: تقريب محسّن يعامل السماء كسطح مستوٍ ويحتفظ بالمكون الموازي لمتجه الموجة ($k_\parallel$)، بينما يقرب المكون العمودي. تستخدم هذه الطريقة نسخة "معاد معايرتها" حيث يتم استبدال $\ell$ بـ $\ell + 1/2$ لتحسين الدقة.

إعداد المحاكاة:

• مصدر البيانات: محاكاة Flagship، وهي محاكاة N-body عالية الدقة (صندوق بحجم $3600 h^{-1}$ ميجا بارسيك) تم ملؤها بالمجرات باستخدام توزيع إشغال الهالات (HOD) ومطابقة وفرة الهالات الفرعية.
• السيناريوهات:
  • DR1 (إصدار البيانات 1): سيناريو محافظ يتضمن 6 نطاقات متساوية التوزيع للإزاحة الحمراء مع قطع للمطال قدره $I_E < 23.5$.
  • DR3 (إصدار البيانات 3): سيناريو المهمة النهائي الذي يتضمن 13 نطاقاً متساوياً للتوزيع للإزاحة الحمراء مع قطع أعمق للمطال قدره $I_E < 24.5$.
• التقييم الإحصائي: قام المؤلفون بتوليد 1,000 طيف قدرة مشوش بناءً على مصفوفة التباين النظرية (التقريب الغاوسي). وقاموا بحساب إحصائية $\chi^2$ لتقريبي ليمبر والسماء المسطحة بالنسبة للنتيجة "الحقيقية" للسماء الكاملة لتحديد ما إذا كان بإمكان التقريبات التمييز إحصائياً عن الحقيقة بالنظر إلى ميزانية الخطأ في المسح.

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

• التحقق الكمي لإقليدس: توفر الورقة أول تقييم صارم قائم على المحاكاة لصلاحية تقريب ليمبر خصيصاً لتجمع المجرات الضوئية في إقليدس عبر تكوينات إصدارات البيانات المختلفة.
• أداء تقريب السماء المسطحة: تُظهر الورقة أن تقريب السماء المسطحة (تحديداً النسخة المعاد معايرتها) يستعيد طيف القدرة الزاوي الكامل بدقة أفضل من 5% عبر جميع المقاييس ونطاقات الإزاحة الحمراء، حتى في النطاقات الضيقة لـ DR3.
• تحليل المقايضة الحاسوبية: تسلط الدراسة الضوء على تحليل التكلفة مقابل الفائدة الحاسوبية. وبينما يعد تقريب السماء المسطحة أبطأ من ليمبر، إلا أنه يوفر تسريعاً كبيراً (بمقدار $\sim$2 ضعف) مقارنة بالحساب الكامل، في حين أن ليمبر أسرع بـ $\sim$200 ضعف.

4. النتائج

تكوين DR1 (6 نطاقات)

• أداء ليمبر: تقريب ليمبر دقيق بما يكفي. الانحرافات عن الحساب الكامل تقع ضمن أشرطة الخطأ الرصدية لـ $1\sigma$.
• الاختبار الإحصائي: توزيع $\Delta\chi^2$ يبلغ ذروته عند 12 ($p \approx 0.038$)، وهو مرفوض بشكل هامشي عند مستوى ثقة 95%، لكن تقريب السماء المسطحة يعمل بشكل مطابق تقريباً للحساب الكامل ($\Delta\chi^2 \approx 0.2$، $p \approx 0.39$).
• الاستنتاج: بالنسبة للنطاقات الأوسع في DR1، فإن السرعة الحاسوبية لليمبر تجعله خياراً قابلاً للتطبيق، حيث أن التحيز ليس ذا دلالة إحصائية بالنسبة لضجيج المسح.

تكوين DR3 (13 نطاقاً)

• أداء ليمبر: تقريب ليمبر يفشل. النطاقات الضيقة للإزاحة الحمراء وعمق المسح الأكبر يجعلان التقريب غير صالح.
  • تتجاوز الأخطاء النسبية 5% وغالباً ما تتجاوز أشرطة الخطأ لـ $1\sigma$ عند قيم $\ell$ المنخفضة ($\ell \lesssim 100$).
  • يعود الفشل إلى عدم قدرة ليمبر على نمذجة حدود RSD و اعتماد $k_\parallel$ بشكل صحيح في النطاقات الضيقة.
• الاختبار الإحصائي: توزيع $\Delta\chi^2$ لليمبر يبلغ ذروته عند 410 ($p < 0.001$)، مما يرفض بقوة فرضية العدم بأن النتائج متسقة مع الحساب الكامل.
• أداء السماء المسطحة: يظل تقريب السماء المسطحة دقيقاً للغاية ($\Delta\chi^2 \approx 5$، $p \approx 0.131$)، حيث يفشل في رفض فرضية العدم، فهو يستعيد الطيف ضمن الشكوك الرصدية.
• الاستنتاج: بالنسبة لإصدار DR3 النهائي، يؤدي تقريب ليمبر إلى نتيجة منحازة لا يمكن استخدامها لعلم الكونيات الدقيق. لذا، فإن تقريب السماء المسطحة مطلوب.

الرؤى الفيزيائية

• تكون أخطاء ليمبر أكثر وضوحاً في الارتباطات الذاتية لنطاقات الإزاحة الحمراء العالية والارتباطات المتبادلة بين النطاقات المتباعدة.
• يُعزى الفشل بشكل كبير إلى تشوهات السرعة في الفضاء الأحمر (RSD) وحدود الكثافة، والتي يقرّبها ليمبر بشكل سيء في النطاقات الضيقة. أما حدود تكبير العدسة (L) فتظل مقربة بشكل جيد نسبياً بواسطة ليمبر حتى في DR3.

5. الأهمية

• توصيات خطوط المعالجة: يوصي المؤلفون بأن خطوط معالجة بيانات إقليدس (مثل CLOE) يجب أن تنتقل إلى تقريب السماء المسطحة لإصدار البيانات النهائي (DR3) لتجنب القيود الكونية المنحازة.
• النهج الهجين: يمكن أن يكون النهج الهجين استراتيجية للمستقبل: استخدام تقريب ليمبر للنطاقات الواسعة (DR1) أو للحدود التي يهيمن عليها التعدس (L-L, D-L) حيث يظل دقيقاً، واستخدام تقريب السماء المسطحة للحدود الحرجة للكثافة و RSD في النطاقات الضيقة.
• التحسين الحاسوبي: على الرغم من أن تقريب السماء المسطحة أبطأ حالياً بنحو $\sim$2 ضعف من الحساب الكامل (بسبب الحاجة إلى تكاملات ثلاثية)، إلا أن الورقة تشير إلى أن المزيد من التحسينات (مثل الحساب المسبق لتكامل $k$ أو استخدام تقنيات شبيهة بـ BLAST) يمكن أن تجعله معياراً عملياً لتحليل MCMC.
• الآثار على المسح الطيفي: تشير النتائج أيضاً إلى أن تحليل 6x2pt للمسح الطيفي (الذي يجمع بين التجمع الطيفي مع التعدس الضعيف) يتطلب تقريب السماء المسطحة إذا استُخدمت نطاقات إزاحة حمراء ضيقة، حيث سيكون ليمبر غير كافٍ.

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