Exploring non-Poisson satellite occupation in HOD models and its impact on 2- and 3-point galaxy clustering

ملخص تقني: استكشاف التوزيع غير البواسوني لإشغال الأقمار الصناعية في نماذج HOD وتأثيره على التجميع المجري ثنائي وثلاثي النقاط

بيان المشكلة
يُعد الربط بين المجرات وهالات المادة المظلمة حجر الزاوية في علم الكونيات الحديث، ويتم نمذجته عادةً باستخدام إطار توزيع إشغال الهالة (HOD). يفترض افتراض قياسي في نماذج HOD أن عدد المجرات التابعة (الأقمار الصناعية) داخل هالة ذات كتلة ثابتة يتبع توزيعاً بواسونيًا (Poisson distribution). ومع ذلك، تشير المحاكاة العددية والتحليلات الرصدية إلى أن مجموعات الأقمار الصناعية قد تظهر انحرافات عن الإحصاء البواسوني الصارم نتيجة تأثيرات مثل انحياز تجميع الهالة (halo assembly bias)، والتقشير المديري (tidal stripping)، وتاريخ السقوط المترابط. وتتجلى هذه الانحرافات في شكل تشتت "دون البواسوني" (sub-Poissonian) أو "فوق البواسوني" (super-Poissonian) في أعداد الأقمار الصناعية، مما يغير العزم العاملي الثاني لعدد الإشغال. وبما أن هذا العزم يتحكم مباشرة في سعة حد الهالة الواحدة (one-halo term) في إحصاءات التجميع، فإن تجاهل السلوك غير البواسوني قد يؤدي إلى أخطاء نظامية في تحليلات التجميع على المقاييس الصغيرة أو القيود الكونية المستمدة من المسوح من المرحلة الرابعة (Stage-IV) مثل DESI وEuclid.

المنهجية
لمعالجة ذلك، قدم المؤلفون توزيع كونواي-ماكسويل-بواسون (CMP) كإمتداد مرن وأدنى لتوزيع بواسون القياسي. يُعرف توزيع CMP بدالة كتلة احتمالية $P(N_{sat}=k | \lambda, \nu) = \frac{\lambda^k}{(k!)^\nu Z(\lambda, \nu)}$، حيث $\nu$ هو معامل التشتت: $\nu=1$ يستعيد التوزيع البواسوني، $\nu < 1$ يعطي سلوكاً فوق بواسوني (التباين > المتوسط)، و $\nu > 1$ يعطي سلوكاً دون بواسوني (التباين < المتوسط).

تتضمن المنهجية ثلاثة مكونات تقنية رئيسية:

1. التقريبات التحليلية والمخططات العددية: نظرًا لأن ثابت التطبيع $Z(\lambda, \nu)$ ليس له شكل مغلق، فقد اشتق المؤلفون تقريبات تحليلية للمتوسط والتباين في نظامين:
   • عند $\lambda$ صغيرة ($\lambda \to 0$): تم اشتقاق توسع متسلسل حتى $O(\lambda^4)$ لتقريب المتوسط والتباين.
   • عند $\lambda$ كبيرة ($\lambda \to \infty$): تم استخدام التوسعات التقاربية القائمة على طرق نقطة السرج (saddle-point methods).
   • القياس العالمي (Global Scaling): تم إنشاء استكمال سلس باستخدام دالة الخطأ للربط بين هذه الأنظمة، مما يضمن خريطة مستمرة بين متوسط الإشغال المستهدف $\langle N_{sat}(M) \rangle$ ومعامل CMP $\lambda$ بدقة تقارب 5% لـ $0.5 < \nu < 2$.
2. توليد الكتالوجات الوهمية (Mock Catalogs): باستخدام حزمة HODDIES ومحاكاة AbacusSummit N-body (في كون $\Lambda$CDM لنموذج Planck 2018)، قام المؤلفون بتوليد 16 كتالوج مجري وهمي تمثل عينة من المجرات الحمراء اللامعة (LRG). تم تثبيت معاملات HOD القياسية، بينما تم تغيير معامل التشتت $\nu$ لـ CMP من 0.5 إلى 2.0.
3. تحليل إحصاءات التجميع: تم قياس تأثير تغيير $\nu$ باستخدام:
   • إحصاءات النقطتين: دالة الارتباط المسقطة ($w_p$)، والمنوال ($\xi_0$) والرباعي ($\xi_2$) لدالة الارتباط في الفضاء الأحمر.
   • إحصاءات النقاط الثلاث: البايسبيكتروم المجري (المنوال $B_{000}$ والرباعي $B_{202}$) المحسوب باستخدام مقدر Sugiyama في نطاق $k$ يتراوح بين $[0, 0.3] \, h/\text{Mpc}$.
   • اللحظات العليا: إحصاءات "العد في الأسطوانات" (Count-in-Cylinders - CiC) لاستقصاء البيئة المجرية المحلية واللحظات العليا للتوزيع.

المساهمات والنتائج الرئيسية

• تنفيذ CMP في HOD: توفر الورقة مخططاً عددياً قوياً لتنفيذ توزيع CMP ضمن أطر HOD، مما ينجح في فصل متوسط إشغال الأقمار الصناعية عن التباين عبر المعامل $\nu$ مع الحفاظ على المتوسط الفيزيائي $\langle N_{sat} \rangle$.
• التأثير على إحصاءات الإشغال: يغير تغيير $\nu$ بشكل كبير توزيع عدد الأقمار الصناعية. عندما يكون $\nu < 1$، يتسع التوزيع (فوق بواسوني)، مما يزيد من احتمالية وجود هالات ذات إشغال عالٍ؛ وعندما يكون $\nu > 1$، يضيق التوزيع (دون بواسوني). يتغير متوسط عدد الأقمار الصناعية بنسبة لا تتجاوز $\pm 5\%$ مقارنة بالحالة البواسونية، وهو انحراف يُعزى إلى عدم تناظر توزيع CMP عند قيم $\lambda$ المنخفضة.
• تأثير التجميع ثنائي النقاط: تؤثر الإحصاءات غير البواسونية بشكل أساسي على حد الهالة الواحدة (المقاييس الصغيرة، $r_p \lesssim 1 \, h^{-1}\text{Mpc}$).
  • يُظهر دالة الارتباط المسقطة $w_p$ انحرافات تصل إلى 10% مقارنة بالحالة البواسونية.
  • يُظهر المنوال والرباعي انحرافات تصل إلى 5%.
  • ترتبط هذه التغييرات بالتباين: التباين الأكبر ( $\nu$ أقل) يعزز قوة التجميع، بينما التباين الأصغر يثبطها.
• تأثير التجميع ثلاثي النقاط (البايسبيكتروم): وُجد أن البايسبيكتروم المجري عند مستوى الشجرة (في أساس Sugiyama) قوي (robust) ضد الافتراضات غير البواسونية في النطاق $k_{\text{max}} = 0.3$.
  • الانحرافات في منوال البايسبيكتروم أقل من 2%.
  • يستنتج المؤلفون أنه بالنسبة للتحليلات واسعة النطاق الحالية والمستقبلية القريبة ($k < 0.3$)، فمن غير المرجح أن تؤدي إحصاءات الأقمار الصناعية غير البواسونية إلى تحيز القيود الكونية المستمدة من البايسبيكتروم، حيث أن الإشارة أصغر من أخطاء البيانات المتوقعة.
• حساسية العد في الأسطوانات (CiC): تُظهر إحصاءات CiC حساسية عالية لتباين إشغال الأقمار الصناعية.
  • تصل الانحرافات عن الحالة البواسونية إلى 30%.
  • يؤكد هذا أن CiC هي أداة قوية لكسر التداخل في معاملات HOD وكشف السلوك غير البواسوني، وهي تتفوق على إحصاءات النقطتين القياسية في هذا الصدد.

الأهمية والادعاءات
يدعي البحث أن توزيع CMP يوفر "آلية دنيا ومنضبطة" لاختبار الفرضية البواسونية ضمن إطار HOD. ويؤكد المؤلفون أنه بينما تعتبر الإحصاءات غير البواسونية بالغة الأهمية لنمذجة دقيقة للتجميع على المقاييس الصغيرة (مما يؤثر على حد الهالة الواحدة بشكل كبير)، فإن تأثيرها على القيود الكونية المستمدة من قياسات الطيف القدرتي والبايسبيكترم عند المقاييس الكبيرة ($k < 0.3$) محدود.

يؤكد العمل أنه مع زيادة الدقة الرصدية، يجب أن يتطور فهم الاتصال بين المجرة والهالة إلى ما وراء الافتراض البواسوني للسيطرة على الأخطاء النظامية. ويقترح المؤلفون أنه بينما قد يكون البايسبيكترم قوياً، فإن الإحصاءات العليا مثل CiC والمتعددات المجرية (galaxy multiplets) ضرورية لتقييد معامل التشتت $\nu$. وتخلص الورقة إلى أن دمج مثل هذه الامتدادات ضروري لاستخلاص ميزانيات خطأ نظامية دقيقة للقياسات الكونية المستقبلية، رغم إشارتهم المتواضعة إلى أن التحليل التفصيلي لتأثير الاستدلال على المعلمات الكونية باستخدام البايسبيكترم يظل موضوعاً للعمل المستقبلي.