Chebyshev interpolation in Einstein-Boltzmann codes

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

لقد وجد باحث في جامعة أوسلو، يُدعى هيرمان سلتيموين، طريقة لجعل لعبة التخمين هذه أكثر دقة وسرعة بشكل ملحوظ. في دراسة جديدة، يوضح سلتيموين أنه من خلال تغيير التقنية الرياضية المستخدمة لملء تلك الفجوات، يمكن للحاسوب أن يحقق دقة أعلى بكثير بعدد أقل بكثير من الحسابات. فالطريقة التقليدية، التي تستخدم "الشرائح التكعيبية" (cubic splines) — وهي تقنية تربط النقاط بمنحنيات صغيرة وسلسة — تعمل بشكل جيد ولكنها تتطلب عددًا كبيرًا من نقاط البيانات للوصول إلى مستويات عالية من الدقة. وقد استبدل سلتيموين ذلك بنهج مختلف يعتمد على "متعددات حدود تشيبيشيف" (Chebyshev polynomials)، وهي أداة رياضية بارعة للغاية في ملاءمة المنحنيات السلسة عندما يتم وضع نقاط البيانات في فترات محددة وغير منتظمة. ومن خلال حل معادلات الفيزياء الصعبة فقط عند هذه النقاط الخاصة، ثم استخدام الطريقة الجديدة لملء الباقي، يمكن للحاسوب إنتاج نتائج تقل فيها أخطاء الاستكمال الرياضي بمقدار عشرة آلاف مرة عن الطريقة القديمة لنفس مقدار العمل.

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

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

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

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