ELLIPSECT: A surface brightness analysis tool for GALFIT 3

1. بيان المشكلة

بينما يُعد GALFIT 3 الأداة القياسية لنمذجة سطوع السطح (SB) ثنائي الأبعاد للمجرات، فإنه يفتقر إلى واجهة رسومية مخصصة وسهلة الاستخدام لتصور مخرجاته وتقييمها كمياً.

• فجوة التصور: يعتمد المستخدمون تقليدياً على أدوات خارجية (مثل ellipse في نظام IRAF) لمقارنة صور المدخلات مع نماذج GALFIT. ومع ذلك، فإن نظام IRAF قد أصبح مهجوراً، كما أن تحويل تنسيقات مخرجات GALFIT المعقدة (مكعبات FITS وملفات المعلمات) لاستخدامها مع أكواد أخرى هو عملية تستغرق وقتاً طويلاً وعرضة للخطأ.
• قصور المنهجيات أحادية البعد: تعاني عمليات مطابقة الإيزوفوت (isophotal fitting) التقليدية أحادية البعد (مثل ellipse في IRAF) في الحقول المزدحمة (العناقيد، المجموعات المتراصة) حيث تتداخل المجرات. وغالباً ما تتطلب عملية "القناع" (masking)، والتي تفشل في مراعاة الضوء الممتد للمجرات المجاورة. علاوة على ذلك، تفترض الطرق أحادية البعد تغيراً في الشكل الإهليلجي (ellipticity) وزوايا الموضع مع نصف القطر، وهو أمر قد يكون إشكالياً عند تحليل عمليات المطابقة متعددة المكونات المتزامنة حيث تكون المكونات مستقلة.
• المقاييس المفقودة: يوفر GALFIT 3 معلمات للمكونات الفردية (مثل نصف قطر الفعالية $r_e$ لانتفاخ المجرة - bulge) ولكنه لا يحسب المقاييس العالمية غير المعلمية (non-parametric) للنظام متعدد المكونات بالكامل، مثل نصف قطر الفعالية الكلي ($r^t_e$)، أو نصف قطر بيتروسيان (Petrosian radius)، أو إحصائيات جودة المطابقة المحددة مثل AIC و BIC المصممة خصيصاً للنموذج.

2. المنهجية

EllipSect هي أداة قائمة على لغة بايثون (Python) صُممت كأداة تحليل لما بعد المطابقة، حيث تستوعب ملفات مخرجات GALFIT 3 مباشرة دون الحاجة إلى تحويل يدوي للتنسيقات.

• بيانات المدخلات: تقرأ ملف FITS cube الناتج عن GALFIT 3 (الصورة الأصلية، صورة النموذج، وصورة البواقي) وملف المعلمات (galfit.nn).
• القياس الضوئي القائم على القطاعات: بدلاً من مطابقة الأشكال الإهليلجية مع الإيزوفوتات، يقوم EllipSect بتقسيم المجرة ونموذجها إلى قطاعات زاوية حول المركز باستخدام روتين sectors_photometry من مكتبة mgefit.
  • يستخدم المعلمات الهندسية (زاوية الموضع، الشكل الإهليلجي) من آخر مكون في نموذج GALFIT (عادةً ما يكون القرص الخارجي أو المكون الأسي) كقيم افتراضية ثابتة للتحليل بأكمله، رغم أنها قابلة للتعديل من قبل المستخدم.
  • يحسب سطوع السطح (SB) عند أنصاف أقطار متعددة داخل هذه القطاعات.
• تقدير خلفية السماء: يقدم EllipSect طريقتين لحساب خلفية السماء، وهو ما يفتقر إليه GALFIT 3:
  1. طريقة تدرج السماء (Sky Gradient Method): تبني حلقات إهليلجية متمركزة وتحلل تدرج قيم متوسط السماء لإيجاد المنطقة المستقرة (اتباعاً لـ Barden et al., 2012).
  2. طريقة الصندوق العشوائي (Random Box Method): تضع صناديق عشوائية في الحقل خارج قناع المجرة لحساب وسيط قيمة السماء (اتباعاً لـ Gao &-Ho, 2017).
• الحسابات غير المعلمية: يحسب مقاييس لا يوفرها GALFIT 3، بما في ذلك:
  • نصف القطر الفعال الكلي ($r^t_e$): يُحسب عن طريق جمع تدفق جميع المكونات وإيجاد نصف القطر الذي يشمل 50% من الضوء الكلي عبر طريقة التنصيف (bisection method).
  • أنصاف أقطار بيتروسيان وكرون (Petrosian and Kron Radii).
  • جودة المطابقة (Goodness-of-Fit): معيار معلومات أكايكي (AIC) ومعيار معلومات بايز (BIC)، بما في ذلك BIC بديل يعتمد على عناصر الدقة.
  • المعلمات الهيكلية: التعرج (Bumpiness/Irregularity)، معلمة المد والجزر (Tidal parameter/Deviation from smooth model)، ونسبة الإشارة إلى الضوضاء (S/N).
• التصحيحات: يسترجع تلقائياً معاملات المسافة والامتصاص المجري من NED لحساب القدر المطلق واللمعان. كما يطبق تصحيحات لبهتان سطوع السطح، لكنه يستبعد تصحيحات $k$-corrections (التي تتطلب بيانات لونية).

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

• سير عمل متكامل: يلغي الحاجة إلى التحويل اليدوي للبيانات بين GALFIT 3 وأدوات الرسم، مما يوفر رسوماً توضيحية بجودة النشر (مثل ملفات تعريف SB، والبواقي، وتفكيك المكونات) مباشرة من مخرجات GALFIT.
• التحليل متعدد المكونات: يتيح تصور وتقييم كمي للمكونات الفردية (مثل الانتفاخ، القرص، القضيب) ضمن مطابقة واحدة، مما يسمح للمستخدمين بتحديد المكونات الزائدة أو غير ذات الصلة.
• المقاييس العالمية للأنظمة المعقدة: يقدم حساب نصف القطر الفعال الكلي ($r^t_e$) للمطابقات متعددة المكونات، وهو معلمة حاسمة لمقارنة المجرات مع المحاكاة العددية أو تقييم تدرجات المعدنية، والتي لا يمكن اشتقاقها ببساال عبر قيم $r_e$ للمكونات الفردية.
• تحديد قوي لخلفية السماء: يوفر طرقاً مستقلة لتقدير خلفية السماء للتحقق من أو تحسين قيمة السماء المستخدمة في مطابقة GALFIT، مما يقلل من تدهور المعلمات (parameter degeneracy).
• جاهز للأتمتة: صُمم ليكون قابلاً للبرمجة ومتوافقاً مع خطوط معالجة المطابقة الآلية (مثل كود DGCG لمجرات العناقيد).

4. النتائج والأمثلة

يستعرض البحث قدرات EllipSect باستخدام ثلاث مجموعات بيانات متميزة:

1. Holm 15A (Abell 85 BCG): مجرة من نوع cD تم مطابقتها باستخدام نموذج 7-Gaussian Multi-Gaussian Expansion (MGE). نجح EllipSect في رسم ملف تعريف SB، ومكونات Gaussian الفردية، والبواقي، مما أظهر قدرة النموذج على التقاط النواة الممتدة للمجرة.
2. M61 (2MASS Ks-band): مجرة حلزونية تم مطابقتها بمكونات Sérsic وExponential. أنتجت الأداة تصورات متعددة الرسوم توضح تغيرات SB عند زوايا سمتية مختلفة، مما يسلط الضوء على المناطق التي ينحرف فيها النموذج عن البيانات.
3. NGC 4261 (HST): مطابقة Sérsic واحدة لصورة من تلسكوب هابل (HST). أنتجت الأداة رسوماً عالية الدقة للمجرة، والنموذج، والبواقي، مما أظهر فائدتها مع البيانات الفضائية.

في جميع الحالات، أنتج EllipSect رسوماً توضيحية جاهزة للنشر وصدر ملفات بيانات تحتوي على المتغيرات الفوتومترية (القدر، الأنصاف الأقطار، AIC/BIC) لمزيد من التحليل.

5. الأهمية

• تعزيز الموثوقية: من خلال توفير رابط مباشر بين النمذجة ثنائية الأبعاد والتقييم الكمي، يساعد EllipSect علماء الفلك على التحقق من أن نماذج GALFIT مناسبة فيزيائياً للبيانات، خاصة في البيئات المعقدة والمزدحمة حيث تفشل الطرق أحادية البعد.
• التقييس: يعمل على توحيد حساب المعلمات العالمية (مثل $r^t_e$) للمجرات متعددة المكونات، مما يسهل المقارنات الأكثر دقة بين البيانات الرصدية والنماذج النظرية.
• سهولة الوصول: كأداة مفتوحة المصدر قائمة على بايثون، يقلل EllipSect من حواجز الدخول لتحليل مخرجات GALFIT، مستبدلاً سير عمل IRAF المهجور بأكواد حديثة وقابلة لإعادة الإنتاج.
• قابلية التوسع المستقبلي: صُممت الأداة للتعامل مع مجموعات بيانات كبيرة (مئات المجرات في العناقيد)، مما يجعلها مكوناً حيوياً للمسوحات واسعة النطاق القادمة ودراسات مورفولوجيا المجرات المؤتمتة.

باختصار، يسد EllipSect الفجوة بين النمذجة المعلمية ثنائية الأبعاد والتحليل الفوتومتري التفصيلي، حيث يقدم مجموعة شاملة من الأدوات للتحقق من هياكل المجرات وتصورها وقياسها المستمد من GALFIT 3.