RotCurves: A PYTHON package for efficient modelling and fitting of galactic rotation curves at high-z

ملخص تقني: RotCurves – حزمة برمجية بلغة بايثون للنمذجة والملائمة الفعالة لمنحنيات الدوران المجري عند الإنزياح الأحمر العالي

بيان المشكلة
تعد منحنيات الدوران أساسية لفهم بنية المجرات، وديناميكيتها، وهالات المادة المظلمة. ومع ذلك، فإن تحليل هذه المنحنيات عند الإنزياح الأحمر العالي ($z$) يفرض تحديات كبيرة بسبب خفوت السطوع السطحي ومحدودية الدقة الزاوية. تؤدي هذه القيود إلى "تلطيخ الحزمة" (beam smearing)، وهو تنعيم اصطناعي لتدرجات السرعة ناتج عن تراكب الضوء من مساحات فيزيائية كبيرة فوق المصدر. تصحيح هذا التأثير أمر بالغ الأهمية ولكنه مكلف حاسوبياً. تقوم أدوات النمذجة الأمامية البارامترية الموجودة، مثل dysmalpy، بإنشاء مكعبات بيانات ثلاثية أو رباعية الأبعاد كاملة (تشمل الأبعاد المكانية والطيفية) لمحاكاة الملاحظات. وبينما تعد هذه الطرق دقيقة، إلا أنها مكثفة حاسوبياً، حيث تتطلب غالباً ثوانٍ أو دقائق لكل عملية محاكاة، مما يعيق الاستكشاف السريع للمعلمات، واختبار فرضيات نماذج الكتلة المتعددة، ومعالجة عينات كبيرة من مسوحات وحدة المجال البؤري (IFU) القادمة. هناك حاجة لأداة توازن بين دقة تصحيح تلطيخ الحزمة والكفاءة الحسابية.

المنهجية
يقدم المؤلفون RotCurves، وهي أداة نمذجة أمامية بارامترية مصممة لملائمة منحنيات الدوران المجري بكفاءة من خلال تقريب عملية تلطيخ الحزمة في بُعدين بدلاً من ثلاثة أبعاد.

• بناء النموذج: يقوم الكود ببناء نموذج مجري مثالي متماثل المحور في حالة توازن طرد مركزي. ويقوم بحساب السرعة الدائرية الجوهرية في المستوى المجري المتوسط عن طريق جمع مساهمات مختلف مكونات الكتلة:
  • الباريونات: تُنمذج كأقراص رقيقة أو سميكة للغاية وتكتلات (bulges) باستخدام ملفات تعريف سيرسيك (Sérsic profiles) (بما في ذلك ملفات التعريف الأسية وde Vaucouleurs).
  • المادة المظلمة: تُنمذج كهالات متماثلة كروياً باستخدام ملفات تعريف مثل NFW، و$\alpha$-NFW، وBurkert، وEinasto، وDekel-Zhao. يتضمن الكود خياراً للانكماش الأديباتي (adiabatic contraction) للهالة نتيجة لتراكم الباريونات.
  • الدعم بالضغط: يقوم النموذج بتصحيح سرعة الدوران المرصودة مقابل تدرجات الضغط في الوسط بين النجمي (ISM) المضطرب باستخدام افتراضات التوازن الهيدروستاتيكي، مما يحول السرعة الدائرية الجوهرية إلى سرعة دوران مرصودة بناءً على تشتت السرعة ($\sigma_0$).
• تقريب تلطيخ الحزمة: على عكس الأدوات التي تقوم بتدوير مجرة ثلاثية الأبعاد لتتطابق مع السماء وتراكبها مع دالة انتشار النقطة (PSF)، يقوم RotCurves بإسقاط الحزمة على المستوى المجري المتوسط. بالنسبة لزاوية ميل معينة $i$، يتم تحويل دالة انتشار نقطة غاوسية دائرية إلى دالة غاوسية ثنائية المتغيرات ذات نسبة محور $1/\cos(i)$. ويتم حساب السرعة والتشتت في السرعة المرصودة عند نصف قطر معين كمتوسطات موزونة فوق هذا "الشعاع الفعال"، مع مراعاة ملف الانبعاث (الوزن الضوئي) وإسقاط خط البصر.
• إجراء الملائمة: تتضمن الحزمة محرك استدلال بايزي مدمج باستخدام عينات MCMC من نوع emcee. وهي تقوم بملائمة المنحنيات المرصودة لمنحنيات الدوران عن طريق تقليل دالة الاحتمالية اللوغاريتمية التي تتكون من احتمالية المربعات الصغرى والاعتبارات المسبقة (priors) المحددة من قبل المستخدم (منتظمة أو غاوسية) على المعلمات مثل الكتلة الباريونية، ونصف القطر الفعال، ونسبة الانتفاخ إلى الإجمالي، وتشتت السرعة الجوهري، وكتلة هالة المادة المظلمة.

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

1. الكفاءة الحسابية: يحقق RotCurves زمن تشغيل يبلغ حوالي $10-20$ ميلي ثانية لكل عملية محاكاة لشبكة نموذجية، وهو ما يعادل تقريباً 250 ضعف سرعة dysmalpy (التي تستغرق حوالي $3-4$ ثوانٍ لنفس المهمة). تم تحقيق هذه السرعة من خلال تجنب بناء مكعبات فائقة رباعية الأبعاد كاملة.
2. الإسقاط ثنائي الأبعاد للحزمة: قدمت الورقة البحثية وصادقت على طريقة إسقاط ثنائية الأبعاد لتلطيخ الحزمة تأخذ في الاعتبار زاوية الميل دون الحاجة إلى تكامل حجمي ثلاثي الأبعاد كامل.
3. القياس المرجعي: قام المؤلفون بإجراء قياس مرجعي صارم لـ RotCurves مقابل dysmalpy باستخدام ملاحظات اصطناعية لنموذج مجري نموذجي عند إزاحة حمراء $z=2$.

النتائج

• الدقة مقابل dysmalpy: بالنسبة للأنظمة جيدة الدقة (حيث FWHM لـ PSF $< R_{\text{eff}}$)، يعيد RotCurves إنتاج منحنيات الدوران والتشتت المرصودة المستخرجة من dysmalpy ضمن نطاق 5% حتى $3R_{\text{eff}}$. وفي الأنظمة ذات الدقة الحدية (حيث FWHM لـ PSF $\gtrsim 1.5 R_{\text{eff}}$)، تزدัง الفروقات لتصل إلى 15%.
  • الاختلافات المنهجية: يميل RotCurves إلى التقليل قليلاً من تشتت السرعة في المناطق الداخلية (بنسبة $\sim 10\%$) لأنه يأخذ في الاعتبار الحركات في المستوى المتوسط فقط، بينما تقوم dysmalpy بدمج الجسيمات خارج المستوى المتوسط والتي تؤدي إلى توسيع الخط.
• استعادة المعلمات: باستخدام ملائمة MCMC على البيانات الاصطناعية، نجح RotCurves في استعادة المعلمات الجوهرية للنموذج (الكتلة، نصف القطر، التشتت، كسر المادة المظلمة) عبر نطاق واسع من خصائص المجرات.
  • الانحيازات: تظهر انحيازات منهجية لنصف القطر الفعال ($R_{\text{eff}}$)، والذي غالباً ما يتم تقديره بأكثر من الحقيقة بمعامل $\sim 1.5$، وللكتل القرصية المنخفضة ($M_{\text{disk}} \lesssim 3 \times 10^9 M_\odot$)، حيث يكون كتلة القرص غير محددة بدقة.
  • نسبة الإشارة إلى الضوضاء (S/N): استعادة المعلمات ممتازة عند مستويات S/N عالية ($\gtrsim 25$). عند المستويات المنخفضة، تزداد الانحرافات، خاصة لكتلة القرص وكسور المادة المظلمة. ومع ذلك، تظل كتلة الانتفاخ (bulge mass) وتشتت السرعة الجوهري ثابتين نسبياً حتى عند مستويات S/N المنخفضة.
• حدود زاوية الميل: الأداة هي الأنسب لزوايا الميل بين $10^\circ$ و $80^\circ$.

الأهمية والادعاءات
تضع الورقة البحثية RotCurves بشكل أساسي كـ أداة استكشافية للاختبار السريع لفرضيات نماذج الكتلة ودراسات المعلمات، وليس كبديل للنمذجة ثلاثية الأبعاد الكاملة في جميع السياقات. تكمن أهميتها في قدرتها على:

• معالجة عينات كبيرة من البيانات الرصدية بكفاءة من المسوحات الحالية والمستقبلية لوحدات المجال البؤري الكبيرة (مثل VLT/MAVIS، وJWST/NIRSpec، وELT/MOSAIC).
• تسهيل اختبار الاعتبارات المسبقة المختلفة، وأشكال الملفات التعريفية، والآثار المنهجية دون العبء الحسابي للمحاكاة ثلاثية الأبعاد الكاملة.
• توفير تقديرات سريعة ومعقولة للمعلمات للمجرات عند الإزاحة الحمراء العالية، مما يعمل كنقطة انطلاق لتحليل أكثر تفصيلاً.

يؤكد المؤلفون أنه بينما يعد RotCurves فعالاً ودقيقاً للغاية للأنظمة جيدة الدقة، يجب على المستخدمين توخي الحذر عند تفسير النتائج للأنظمة ضعيفة الدقة ومنخفضة الـ S/N، حيث تكون آثار تلطيخ الحزمة والتحيزات المنهجية (خاصة في $R_{\text{eff}}$) في ذروتها. الكود متاح للجمهور على GitHub.