How X-rays heat the IGM in different 21-cm simulation codes: a comparison between Licorice and Beorn

تخيل الكون كمحيط عملاق غير مرئي من الغاز، يتكون معظمه من الهيدروجين، ويمتد في كل الاتجاهات. خلال بضع مئات الملايين من السنين الأولى بعد الانفجار العظيم، كان هذا المحيط مظلماً وبارداً وصامتاً. ثم اشتعلت أول النجوم والمجرات مثل المنارات في الضباب. بدأت إشعاعاتها المكثفة في تسخين الغاز وتمزيق الإلكترونات بعيداً عن الذرات، وهي عملية تسمى "إعادة التأين". يتوق علماء الفلك لرؤية هذا العصر وهو يتكشف لأن هذا العصر يحمل أسرار كيفية نمو كوننا. وللقيام بذلك، يقومون ببناء تلسكوب راديوي ضخم يسمى مصفوفة الكيلومتر المربع (SKA)، والذي يعمل كأذن عملاقة تستمع إلى همس خافت: "إشارة الـ 21 سم". هذه الإشارة هي موجة راديوية محددة تنبعث من الهيدروجين المتعادل، وتتغير قوتها بناءً على مدى سخونة الغاز. إذا كان الغاز بارداً، فستبدو الإشارة مختلفة عما لو كان دافئاً. ومن خلال قياس هذه الإشارة، يأمل العلماء في رسم خريطة لتاريخ تسخين الكون.

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

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

ولجعل المواجهة عادلة، لم يكتفِ العلماء بجعل الكودين يعملان بإعداداتهما العشوائية الخاصة. بدلاً من ذلك، أخذوا "المكونات" (النجوم وإنتاجها الضوئي) من محاكاة Licorice التفصيلية وغذوا بها Beorn. هذا ضمن أن كلا الكودين يحاولان طهي نفس الكون، ولكن باستخدام طرق مختلفة. ثم قارنوا النتائج: درجة حرارة الغاز، سطوع إشارة الـ 21 سم، والأنماط الإحصائية لكيفية تقلب تلك الإشارة عبر السماء.

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

تبين أن هذه الاختلافات الصغيرة في "نسيج" الحرارة مهمة بشكل مفاجئ. فعندما حسب الفريق "طيف القدرة" (power spectrum) — وهي طريقة معقدة لقياس مدى تقلب الإشارة عند مستويات مختلفة — وجدوا فرقاً بنحو 30% بين الكودين. ولوضع ذلك في المنظور، إذا كنت تحاول تخمين وصفة كعكة بناءً على قضمة واحدة، فإن فرق 30% في الطعم سيعني أنك ستخمن المكونات بشكل خاطئ تماماً.

ثم تعمق الباحثون لمعرفة لماذا اختلف الكودان. واكتشفوا أن الطرق المختصرة المستخدمة في الكود الأسرع هي الجاني. وتحديداً، افترض الكود أحادي البعد أن كثافة الغاز هي نفسها في كل مكان (مثل افتراض أن المحيط مسطح تماماً) واستخدم صيغة مبسطة لكيفية ترسيب الأشعة السينية للحرارة فعلياً. في الواقع، الغاز متكتل، ويعتمد التسخين على مكان وجود الغاز بالضبط. وعندما أجبر الفريق محاكاة Licorice التفصيلية ثلاثية الأبعاد على استخدام نفس هذه الاختصارات المبسطة، اقتربت النتائج من كود Beorn، مما أكد أن هذه التقريبات كانت مصدر الخطأ.

أخيراً، طرح الفريق السؤال النهائي: "إذا استخدمنا هذه المحاكيات المختلفة لتحليل البيانات من تلسكوب SKA المستقبلي، فهل سنحصل على إجابات خاطئة؟". قاموا بمحاكاة سيناريو حيث سيقوم التلسكوب بالمراقبة لمدة 100 ساعة. ووجدوا أن الاختلافات بين الكودين ستؤدي إلى انحياز كبير في الخصائص المقدرة للنجوم الأولى. من الناحية الإحصائية، كان الخطأ كبيراً بما يكفي ليكون أكبر من 1 سيجما (مقياس للثقة)، وبالنسبة لبعض المعايير، كان حتى أعلى من ذلك. وهذا يعني أنه إذا استخدم علماء الفلك الكود الأسرع والأبسط لتفسير البيانات الحقيقية، فقد يستنتجون بثقة أن النجوم الأولى كانت مختلفة عما هي عليه في الواقع.

تخلص الدراسة إلى أنه بينما تعد الأكواد الأسرع مفيدة وتصيب في سرد "القصة المتوسطة"، إلا أنها ليست دقيقة بما يكفي لعصر SKA. إن فرق الـ 30% في طيف القدرة وما نتج عنه من انحياز في البيانات يشير إلى أنه لا يمكننا ببساطة الوثوق بالطرق المختصرة عندما نحاول قياس الكون بدقة متناهية. ويشير المؤلفون إلى أنه بينما تعد الأكود متشابهة بما يكفي للتلسكوبات الحالية الأقل حساسية، إلا أن الكثير من العمل مطلوب لتحسين الأكود الأسرع أو تطوير طرق جديدة لتصحيح أخطائها قبل أن يبدأ SKA في الاستماع. إنه تذكير بأنه في السعي لفهم فجر الكون، تكمن التفاصيل في أدق التفاصيل.