A New PSF Deconvolution Algorithm: Simultaneous Spatial Resolution Enhancement and Point Source Removal for Morphological Analysis of AGN Host Galaxies

ملخص تقني: خوارزمية جديدة لفك التواء دالة انتشار النقطة (PSF) لتحليل المجرات المضيفة للنوى المجرة النشطة (AGN)

بيان المشكلة
غالبًا ما يصعب دراسة مورفولوجيا المجرات المضيفة للنوى المجرة النشطة (AGN) بسبب تلوث الضوء الصادر من المجرة المضيفة بواسطة المصدر النقطي المركزي الساطع (الـ AGN). وبينما تستطيع التلسكوبات الفضائية عالية الدقة مثل هابل (Hubble) وجيمس ويب (JWST) رصد هذه البنى، إلا أن مجال رؤيتها الضيق يحد من الدراسات الإحصائية لتشمل عينات صغيرة (أقل من $10^4$ جسم). توفر المسوحات الأرضية (مثل Subaru/HSC، وRubin، وEuclid، وRoman) عينات ضخمة، لكنها تعاني من دقة مكانية أقل وتشويه كبير ناتج عن دالة انتشار النقطة (PSF). وتعتبر الطرق التقليدية، مثل التفكيك ثنائي الأبعاد للصور (مثل GALFIT)، حساسة لعمليات الطرح الزائد للـ PSF والتركيبات المعقدة، كما أنها تواجه صعوبة في فصل المجرة المضيفة عن الـ AGN دون إدخال عدم يقين في المعلمات الهيكلية مثل مؤشر سيرسيك (Sérnic index). غالبًا ما تفشل خوارزميات فك التواء الـ PSF الموجودة في إخراج مكون ممتد منفصل، ومكون نقطي منفصل، وصورة معاد بناؤها مدمجة في آن واحد، وهي عناصر ضرورية للتحليل المورفولوجي والضوئي القوي للمجرات المضيفة للـ AGN في الأعماق البعيدة.

المنهجية
يقترح المؤلفون خوارزمية جديدة لفك التواء الـ PSF مصمحة لتعزيز الدقة المكانية وإزالة المصدر النقطي المركزي في آن واحد. تنمذج الخوارزمية الصورة الأصلية $x$ باعتبارها مجموع مكونين متمايزين: مكون ممتد $I_{sm}$ (المجرة المضيفة) ومكون نقطي $I_{sp}$ (الـ AGN). يفترض نموذج الرصد أن التدفق الملحوظ يتبع توزيعًا غاوسيًا بناءً على التواء هذين المكونين مع الـ PSF.

لإعادة بناء الصورة الأصلية، تقوم الخوارماية بتقليل دالة الهدف للوغاريتم الخلفي السالب (negative log-posterior) مع الخضوع لثلاثة قيود محددة:

1. قيد النعومة ($I_{sm}$): تنظيم يشبه "تيكونوف" (Tikhonov) يُطبق على المكون الممتد لفرض النعومة المكانية، مما ينمذج بنية المجرة المضيفة.
2. قيد التشتت ($I_{sp}$): تنظيم بمعيار $L_1$ (طريقة LASSO) يُطبق على المكون النقطي لتركيز تدفق الـ AGN في عدد قليل من البكسلات.
3. قيد توازن المصدر النقطي: قيد مبتكر يعتمد على حاصل الضرب البكسلي $I_{sm} \times I_{sp}$. يعمل هذا المصطلح على معاقبة الوجود المتزامن لقيم تدفق عالية في كل من المكون الممتد والمكون النقطي عند نفس موقع البكسل. وهذا يمنع "تسرب" تدفق المصدر النقطي إلى المجرة المضيفة (الطرح الزائد) أو كبت التدفق المركزي للمجرة المضيفة (الطرح الناقص).

يتم حل عملية التحسين باستخدام خوارزمية التقليص والعتبة التكرارية (ISTA) مع مخطط تراجع لثابت ليبشيتز (Lipschitz constant) لضمان تقارب مستقر ورتيب. تُخرج الخوارزمية ثلاث صور: المجرة المضيفة المنفصلة ($I_{sm}$)، والـ AGN المنفصل ($I_{sp}$)، وإعادة البناء المدمجة ($I_{sm} + I_{sp}$).

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

• خوارزمية مبتكرة: تطوير طريقة لفك التواء الـ PSF ذات مخرجات ثلاثية مصممة خصيصًا لصور الـ AGN، حيث تفصل بوضوح بين المجرة المضيفة والـ AGN مع إعادة بناء الصورة الكلية.
• قيد توازن المصدر النقطي: تقديم حد تنظيم جديد ($I_{sm} \times I_{sp}$) يعمل بفعالية على فك الارتباط بين المكونات الممتدة والنقطية، مما يعالج مشكلة الاعتماد المتبادل بين المعلمات الفائقة التي تعيب الطرق السابقة.
  টি تعزيز الدقة من الأرض إلى الفضاء: إثبات أن الصور الأرضية (Subaru/HSC) يمكن فك التواءها لتحقيق دقة مكانية تضاهي صور تلسكوب هابل (HST) مع إزالة تلوث الـ AGN.

النتائج
تم اختبار الخوارزمية على كل من صور الـ AGN الاصطناعية وبيانات حقيقية من مجال Subaru/HSC COSMOS (الأجرام J0959 عند $z=0.34$ وJ1000 عند $z=1.38$)، باستخدام صور HST/ACS كمرجع حقيقي.

• التحليل النوعي: تُظهر صور المجرة المضيفة المعاد بناؤها ($I_{sm}$) تعزيزًا ملحوظًا في الدقة المكانية، مما يكشف عن بنى دقيقة مثل الأذرع الحلزونية والتركيبات المتكتلة بتباين عالٍ. كما تمت إزالة المصدر النقطي المركزي بفعالية دون حدوث طرح زائد أو ناقص ملموس. وتتشابه الصور المعاد بناؤها المدمجة ($I_{sm} + I_{sp}$) بشكل وثيق مع الصور المرجعية من هابل.
• التحليل الكمي:
  • التقارب: انخفضت مسافة هيلينجر (Hellinger distance) بين التكرارات بسرعة، واستقرت بعد حوالي 100 تكرار.
  • الدقة: كان جذر متوسط مربع الخطأ (RMSE) بين الصورة المعاد بناؤها (بعد التوائها مع PSF الخاص بـ HSC) والملاحظة الأصلية لـ HSC منخفضًا للغاية ($\sim 10^{-6}$ إلى $10^{-5}$). كما كان RMSE بين الصورة المعاد بناؤها ومرجع هابل صغيرًا أيضًا ($\sim 10^{-2}$).
  • تحسين الدقة: تحسن عرض نصف القمة الكامل (FWHM) لبنى المجرة المضيفة المعاد بناؤها بمعامل يتراوح بين $2-3$مقارنة بصور HSC الأصلية، ليصل إلى مستويات تضاهي هابل (على سبيل المثال، من$1.06''$إلى$0.57''$ لـ J0959).
  • فعالية القيود: ثبت أن إدراج قيد توازن المصدر النقطي ($\alpha_{bl}$) كان أمرًا حاسمًا للـ AGNs الحقيقية، حيث أدى إلى تقليل قيم FWHM بشكل أكبر وتحسين إزالة المصدر النقطي مقارنة بالتشغيل بدون هذا القيد.

الأهمية
يدعي البحث أن هذه الخوارزمية تمكن من إجراء دراسات مورفولوجية إحصائية للمجرات المضيفة للـ AGN البعيدة باستخدام بيانات المسوحات واسعة المجال من المرافق الحالية والمستقبلية، بما في ذلك مرصد فيرا روبين (Vera C. Rubin)، وتلسكوب إقليدس (Euclid)، وتلسكوب رومان (Roman). ومن خلال الإزالة الفعالة لتلوث الـ AGN وتعزيز الدقة، تسمح هذه الطريقة بتطبيق المؤشرات المورفولوجية غير البارامترية (مثل التركيز، وعدم التماثل، ومعامل جيني) على عينات ضخمة (أكثر من $10^4$ جسم) دون التحيزات التي كانت تسببها طاقة الـ AGN سابقًا. ومن المتوقع أن تحل هذه القدرة الجدل القائم حول العلاقة بين نشاط الـ AGN وعمليات اندماج/تفاعل المجرات. بالإضافة إلى ذلك، تقلل هذه الطريقة من الاعتماد على التصنيف البصري ويمكن توسيعها لتشمل دراسات علاقات الثقب الأسود الفائق الكتلة والانتفاخ (bulge)، واستخراج البنى الأخرى مثل الكتل النجمية العملاقة.