cloelib: A Flexible Python Library for Computing Cosmological Observables in the Euclid Era
تقدم الورقة البحثية **cloelib**، وهي مكتبة بايثون مرنة وقابلة للاشتقاق مصممة لحساب الملاحظات الكونية لمهمة إقليدس (Euclid) عبر دمج بروتوكولات نمطية مع التفاضل التلقائي القائم على JAX لسد الفجوة بين علم الكونيات العددي التقليدي وتقنيات الاستدلال الحديثة للتعلم الآلي.
المؤلفون الأصليون:Marco Bonici (cloe-org maintainers), Guadalupe Cañas-Herrera (cloe-org maintainers), Pedro Carrilho (cloe-org maintainers), Santiago Casas (cloe-org maintainers), Chiara Moretti (cloe-org maintainers)Marco Bonici (cloe-org maintainers), Guadalupe Cañas-Herrera (cloe-org maintainers), Pedro Carrilho (cloe-org maintainers), Santiago Casas (cloe-org maintainers), Chiara Moretti (cloe-org maintainers), Andrea Pezzotta (cloe-org maintainers), Michel Aguena (cloe-org contributors), Giovanni Aricó (cloe-org contributors), Zahra Baghkhani (cloe-org contributors), Matteo Baratto (cloe-org contributors), Emilio Bellini (cloe-org contributors), Jip de Buck (cloe-org contributors), Klara Bertmann (cloe-org contributors), Ben Bose (cloe-org contributors), Jeger C. Broxterman (cloe-org contributors), Pierre Burger (cloe-org contributors), Carmelita Carbone (cloe-org contributors), Chaitanya Chawak (cloe-org contributors), Jose Coloma-Nadal (cloe-org contributors), Martin Crocce (cloe-org contributors), Stefano Davini (cloe-org contributors), Christopher A. J. Duncan (cloe-org contributors), Samuel Farrens (cloe-org contributors), Lisa Goh (cloe-org contributors), Nastassia Grimm (cloe-org contributors), Alex Hall (cloe-org contributors), Ryusei R. Kano (cloe-org contributors), Martin Kärcher (cloe-org contributors), Felicitas Keil (cloe-org contributors), Raphaël Kou (cloe-org contributors), Laila Linke (cloe-org contributors), Arthur Loureiro (cloe-org contributors), Katarina Markovic (cloe-org contributors), David Navarro-Gironés (cloe-org contributors), Filippo Oppizzi (cloe-org contributors), Gabriele Parimbelli (cloe-org contributors), Anna Porredon (cloe-org contributors), Robert Reischke (cloe-org contributors), Jaime Ruiz Zapatero (cloe-org contributors), Iñigo Sáez-Casares (cloe-org contributors), Ziad Sakr (cloe-org contributors), Neel Shah (cloe-org contributors), Davide Sciotti (cloe-org contributors), Matthieu Schaller (cloe-org contributors), Ivan Sladoljev (cloe-org contributors), Arghavan Souki (cloe-org contributors), Sankarshana Srinivasan (cloe-org contributors), Konstantinos Tanidis (cloe-org contributors), Peter L. Taylor (cloe-org contributors), Nicolas Tessore (cloe-org contributors), Linus Thummel (cloe-org contributors), Maria Tsedrik (cloe-org contributors), Isaac Tutusaus (cloe-org contributors), Casper Vedder (cloe-org contributors), Angus H. Wright (cloe-org contributors), Miguel Zumalacárregui (cloe-org contributors), Joe Zuntz (cloe-org contributors)
تخيل الكون كأنه أحجية (بازل) عملاقة ومعقدة يحاول علماء الفلك حلها. لسنوات، امتلك العلماء قطعاً من هذه الأحجية (بيانات من تلسكوبات مثل Euclid وDESI وغيرها)، لكنهم كانوا بحاجة إلى طريقة أفضل لتجميع هذه القطع معاً لرؤية الصورة الكاملة.
تقدم هذه الورقة البحثية cloelib، وهو "صندوق أدوات" رقمي جديد صُمم خصيصاً لهذه المهمة. إليك شرح مبسط لماهية هذا الصندوق ولماذا يكتسب أهميته، باستخدام بعض التشبيهات من الحياة اليومية.
١. المشكلة: "السكين السويسري" الذي كان ثقيلاً للغاية
في الماضي، كان لدى علماء الفلك برمجيات يمكنها حساب كيفية توسع الكون أو كيفية تجمع المجرات. ومع ذلك، كانت هذه الأدوات تشبه صندوق أدوات صلباً وثقيلاً. إذا أردت إضافة أداة جديدة (مثل طريقة جديدة لقياس أشكال المجرات) أو استبدال جزء ما (مثل استخدام حاسبة مختلفة لتاريخ الكون)، فغالباً ما كان عليك تفكيك صندوق الأدوات بالكامل وإعادة بنائه. جعل هذا الأمر العملية بطيئة، وصعبة التحديث، ويصعب فيها دمج الأفكار المختلفة.
٢. الحل: نظام "ليجو" (Lego)
cloelib يشبه مجموعة من قطع الليجو القابلة للتبديل.
التصميم المجزأ (Modular Design): بدلاً من كونه كتلة واحدة ضخمة غير قابلة للتغيير، بُني البرنامج في طبقات. لديك طبقة أساسية (تاريخ الكون)، وطبقة وسطى (كيف تتحرك المادة)، وطبقة عليا (ما نراه فعلياً عبر التلسكوبات).
التوصيل والتشغيل (Plug-and-Play): نظرًا لأن القطع تتناسب مع بعضها البعض بشكل مثالي (باستخدام "بروتوكولات"، وهي بمثابة موصلات قياسية)، يمكنك استبدال حاسبة بأخرى دون كسر النظام بأكره. إذا أردت استخدام "محاكي" (emulator) فائق السرعة (وهو بمثابة حاسبة مختصرة) بدلاً من واحدة مفصلة وبطيئة، يمكنك ببساطة تركيبه. لن يحتاج باقي عملك إلى أي تغيير.
٣. القوة الخارقة: الرياضيات "المفكرة"
أحد أكثر الميزات تميزاً في cloelib هو أنه قابل للتفاضل (differentiable).
التشبيه: تخيل أنك تحاول إيجاد الوصفة المثالية لكعكة. في الأيام الخوالي، كنت ستخبز الكعكة، ثم تتذوقها، ثم تخمن ما يجب تغييره، ثم تخبزها مرة أخرى. هذا أمر بطيء.
طريقة cloelib: بفضل ارتباطه بتقنية تسمى JAX، يمكن لهذا الصندوق أن يخبرك بالضبط كيف سيغير تغيير مكون واحد (مثل إضافة القليل من السكر) طعم الكعكة، بشكل فوري. إنه يحسب "الميل" (gradient) أو درجة التغير تلقائياً. وهذا يسمح للعلماء باستخدام تقنيات حاسوبية متقدمة (مثل تعلم الآلة) لإيجاد أفضل ملاءمة لبيانات الكون بشكل أسرع بكثير من ذي قبل.
٤. ماذا يفعل في الواقع؟
يساعد صندوق الأدوات علماء الفلك في حساب نوعين رئيسيين من الملاحظات:
الصور الضوئية (Photometric - التقاط الصور): يساعد في تحليل الصور الضبابية للمجرات لقياس مدى تمددها بفعل الجاذبية (القص الكوني - cosmic shear) وكيفية تجمعها معاً. كما يعالج "الأخطاء المنهجية"، والتي تشبه البقع على عدسة الكاميرا التي تحتاج إلى تنظيف رياضياً.
التحليل الطيفي (Spectroscopic - أخذ البصمات): يحلل "ألوان" الضوء (الأطياف) المنبعثة من المجرات لقياس سرعتها ومسافتها بدقة عالية.
٥. لماذا يهم هذا في "عصر إقليدس" (Euclid Era)؟
مهمة وكالة الفضاء الأوروبية Euclid على وشك إرسال كمية هائلة من البيانات عالية الدقة حول الكون. هذه البيانات ضخمة ودقيقة جداً لدرجة أن البرمجيات القديمة ستعاني لمواكبتها.
السرعة: تزعم الورقة البحثية أن cloelib سريع للغاية. على سبيل المثال، حساب خرائط معقدة للكون كان يستغرق وقتاً، أصبح الآن يستغرق أجزاءً من الثانية.
المرونة: يتيح للعلماء اختبار نظريات مختلفة حول الكون (مثل كيفية عمل الطاقة المظلمة) دون الحاجة لإعادة كتابة الكود الخاص بهم في كل مرة.
معيار مجتمعي: تم تصميمه ليكون "الأداة المرجعية الرسمية" لمهمة Euclid، مما يضمن أن الجميع في المجتمع العلمي يستخدمون نفس القواعد لتفسير البيانات.
الملخص
فكر في cloelib كأنه مترجم عالمي وحاسبة فائقة السرعة لعلماء الكونيات. إنه يأخذ الرياضيات المعقدة والفوضوية للكون ويحولها إلى نظام مرن، سريع، وسهل الاستخدام. إنه يسمح للعلماء باستبدال النظريات المختلفة مثل تغيير العدسات في الكاميرا، مما يضمن أنه عندما تصل البيانات الجديدة من تلسكوب Euclid، سيكونون مستعدين لحل أعظم أسرار الكون بكفاءة.
ملخص تقني لمكتبة cloelib: مكتبة بايثون مرنة لحساب الملاحظات الكونية في عصر إقليدس
بيان المشكلة يدخل مجال علم الكونيات الرصدي عصر الدقة مدفوعاً بمسوحات المجرات من المرحلة الرابعة، بما في ذلك مهمة "إقليدس" (Euclid) التابعة لوكالة الفضاء الأوروبية، وDESI، ومرصد فيرا روبين، وتلسكوب نانسي غريس رومان الفضائي. تولد هذه المهام مجموعات بيانات هائلة تتطلب تنبؤات نظرية دقيقة وفعالة حاسوبياً لتقييد النماذج الكونية. وبينما توفر الأطر الحالية (مثل CCL وCosmoSIS وCAMB وCLASS) أدوات قوية، إلا أنها غالباً ما تفتقر إلى المرونة اللازمة لدمج المكونات المتنوعة الموجودة مسبقاً بسلاسة، أو لاستكشاف نماذج نظرية واسعة إلى جانب التأثيرات المنهجية الشاملة. علاوة على ذلك، عانت برمجية CLOE الأصلية (التي كانت مركزية لتحليلات إقليدس) من قيود هيكلية؛ حيث كان توسيعها ليشمل ملاحظات جديدة أو خيارات نمذجة يتطلب تعديلات جذرية، مما أدى إلى تراكم الديون التقنية وانخفاض القدرة على الصيانة. فضلاً عن ذلك، لم يكن الإطار مصمماً في الأصل لدعم التطوير المستقل والدمج للنماذج النظرية الجديدة أو التأثيرات المنهجية ضمن مسار موحد، كما لم يدعم بشكل أصيل تقنيات الاستدلال الحديثة القائمة على التدرج (gradient-based inference).
المنهجية وتصميم البرمجيات تُعد cloelib مكتبة تعتمد لغة بايثون بالكامل، صُممت لتشغيل المرونة والكفاءة المطلوبة لعلم الكونيات الحديث. وهي تمثل تطوراً هيكلياً لبرمجية CLOE، حيث قدمت بنية معيارية تعتمد على بروتوكولات بايثون (PEP 544) لفرض نهج "التوصيل والتشغيل" (plug-and-play).
تنظم المكتبة الوظائف في أربع طبقات متميزة:
الخلفيات الكونية (Cosmological Backgrounds): التعامل مع تاريخ التوسع والمسافات.
نظرية الاضطراب (Perturbation Theory): إدارة أطياف القدرة للمادة الخطية وغير الخطية.
واجهات قائمة على البروتوكولات (Protocol-Based Interfaces): تتفاعل المكونات عبر واجهات محددة جيداً (مثل Background وPerturbations وTracer وSpectroPower) بدلاً من هرميات الوراثة الجامدة. وهذا يسمح للباحثين بالتبديل بين خلفيات نظرية مختلفة (مثل حلول بولتزمان مثل CAMB وCLASS، أو المماثلات/Emulators مثل CosmoPower وBACCOemu وEuclidEmulator2) دون تعديل مسار التحليل.
التكامل مع JAX: الميزة الجوهرية هي التكامل الأصيل مع JAX، مما يتيح التفاضل التلقائي (automatic differentiation)، والحساب القائم على التدرج الفعال، والتجميع في الوقت المناسب (JIT compilation). يحول هذا المسارات التقليدية إلى برامج قابلة للتفاضل بالكامل، مما يسهل طرق الاستدلال المتقدمة مثل مونت كارلو لهاملتون (Hamiltonian Monte Carlo) والنمذجة القائمة على الشبكات العصبية.
نمذجة التأثيرات المنهجية (Systematic Modeling): تتضمن المكتبة تنفيذات للتأثيرات المنهجية الخاصة بالمسح، مثل تصحيحات ألكوك-باتشينسكي (Alcock–Paczynski)، ومعايرة القص والإزاحة الحمراء الضوئية، والنقاء الطيفي.
تحسين الأداء: تستخدم الأجزاء الحرجة التجميع في الوقت المناسب (JIT compilation)، ويقوم نظام التخزين المؤقت (caching) بتجنب التقييمات غير الضرورية. تدعم المكتبة نظريات الاضطراب الخطية وغير الخطية، وتتصل مع المماثلات (emulators) المتطورة (مثل HMCode2020Emu وFlamingoBaryonResponseEmulator) والأكواد الطيفية (PBJ وcomet-emu).
المساهمات والقدرات الرئيسية
واجهة موحدة: تُعد cloelib أول كود للهياكل الكونية واسعة النطاق ينفذ واجهة موحدة لعدة أكواد كونية باستخدام بروتوكولات بايثون، مما يتيح مقارنة سلسة بين مختلف التنفيذات النظرية ضمن إعداد متسق.
علم الكونيات القابل للتفاضل: من خلال الاستفادة من JAX، تتيح المكتبة تدريب المماثلات (emulators) للشبكات العصبية واستخدام العينات القائمة على التدرج، والتي غالباً ما يصعب تنفيذها بكفاءة في الأطر التقليدية. ،* التحليل الشامل (End-to-End): متكاملة مع وحدة الاحتمالية المصاحبة cloelike، وتدعم المكتبة السلسلة الكاملة من حساب الملاحظات إلى تقييم الاحتمالية وأخذ عينات التوزيع البعدي للمعلمات.
القابلية للتوسع المعيارية: تسمح استخدامات "الخلط" (mixins) والبروتوكولات في بايثون للمستخدمين بتوسيع الوظائف (مثل إضافة التأثيرات الباريونية) دون تغيير التنفيذات الأساسية.
النتائج والأداء يسجل البحث اختبارات أداء حاسوبي على شريحة Apple M4:
التهيئة: تستغرق تهيئة فئات تتبع القص والموضع حوالي 0.4 ثانية؛ وتصبح الاستدعاءات اللاحقة فورية تقريباً بسبب التخزين المؤقت.
المجسات الضوئية (Photometric Probes): يستغرق حساب أطياف القدرة الزاوية لتحليل 3x2-pt (3000 متعدد قطب، 512 قيمة k، 500 إزاحة حمراء، 6 حاويات) باستخدام تقريب ليمبر (Limber approximation) حوالي 0.06 ثانية لكل استدعاء بعد مرحلة تهيئة مدتها ثانية واحدة. بالنسبة لتكوين شبيه بـ Euclid مع 13 حاوية (351 طيفاً)، يبلغ الوقت حوالي 0.1 ثانية لكل استدعاء.
المجسات الطيفية (Spectroscopic Probes): يستغرق حساب متعددات أقطاب ليجاندر الطيفية كاملة الشكل لحاوية إزاحة حمراء واحدة حوالي 3 مللي ثانية.
المقارنة: تمثل هذه النتائج تسريعاً كبيراً مقارمة ببرمجية CLOE السابقة، مع الحفاظ على دقة تصل إلى مستوى النسب المئوية من خلال استخدام أحدث المماثلات (emulators).
الأهمية يضع هذا البحث cloelib كبنية تحتية مرجعية للحسابات النظرية لمهمة إقليدس في إصدارها الكوني الأول. تكمن أهميتها في سد الفجوة بين علم الكونيات الرقمي التقليدي وتقنيات التحسين الحديثة ونهج التعلم الآلي الناشئة. ومن خلال توفير منصة مرنة وقابلة للتوسع وقابلة للتفاضل، تعالج cloelib القيود الهيكلية للأدوات السابقة، مما يتيح تحليلات إحصائية أسرع وأكثر سلاسة. لقد صُممت لتلبية متطلبات الجيل القادم من تجارب علم الكونيات الدقيق، مما يسهل الدراسات المنهجية، ومقارنات النماذج، والتحقق المتقاطع القوي للنتائج الكونية ضمن كونسورتيوم إقليدس والمجتمع العلمي الأوسع.