Analytical model for the photomultiplier single photoelectron response including the electron back-scattering contribution

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

لا تزال أنابيب المضاعف الضوئي (PMTs) هي المعيار المعتمد للكشف عن مستويات الضوء المنخفضة جدًا عبر مساحات واسعة. وتعد استجابة الفوتوإلكترون المنفرد (Single Photoelectron - SPE) خاصية بالغة الأهمية لمعايرتها ومحاكاتها.

• القيود الحالية: تصف النماذج الحالية بشكل جيد ذروة الفوتوإلكترون المنفرد المضخمة بالكامل (باستخدام توزيعات بواسون أو غاما الموزونة)، لكنها تفشل في نمذجة المنطقة الطيفية الواقعة بين قاعدة الإلكترونيات (pedestal) وذروة الـ SPE الرئيسية بدقة.
• فجوة "الضجيج": غالبًا ما تُعزى هذه المنطقة المتوسطة إلى "ضجيج" عام، ويتم نمذجتها باستخدام دوال عشوائية (مثل التوزيعات الغاوسية المعدلة أسيًا) دون أساس فيزيائي صارم.
• الآلية الفيزيائية: جزء كبير من الفوتوإلكترونات الأولية التي تصطدم بالدينود الأول يخضع لعملية تشتت خلفي غير مرن (inelastic back-scattering). فبدلاً من إيداع كامل طاقتها الحركية لتوليد إلكترونات ثانوية، تتشتت وتفقد جزءًا فقط من طاقتها. يؤدي هذا إلى توليد إشارات "مضخمة جزئيًا" تخلق توزيعًا متميزًا غير غاوسي في منطقة الشحنة المنخفضة من الطيف.

2. المنهجية

قام المؤلفون بتطوير نموذج تحليلي شامل يعتمد على الوصف الفيزيائي للتشتت الخلفي للإلكترونات كما ورد في كتاب (The Photomultiplier Handbook لـ A. G. Wright). يفصل النموذج عملية التضخيم إلى مرحلتين: الدينود الأول وسلسلة الـ $N-1$ دينود المتبقية.

أ. الإطار النظري

يبني النموذج دالة استجابة الـ SPE الكلية، $\mathcal{P}_{SPE}(x)$، كمجموع موزون لأربعة مكونات فيزيائية متميزة:

1. الإلكترونات المضخمة بالكامل (FA):

   • الإلكترونات الأولية التي تودع كامل طاقتها عند الدينود الأول.
   • يتم نمذجتها كـ توزيع بواسون مستمر موزون مدمج (convolved) مع دقة التوزيع الغاوسي لسلسلة الدينودات اللاحقة.
   • الصيغة: $\mathcal{P}_{FA}(x) = \rho e^{-\rho \mu} (\rho \mu)^{\rho x} / \Gamma(1 + \rho x)$.

2. الإلكترونات المضخمة جزئيًا (PA) (التشتت الخلفي):

   • الإلكترونات الأولية التي تتشتت خلفيًا بشكل غير مرن، وتودع فقط كسرًا $(1-\kappa)$ من طاقتها.
   • يتم اشتقاق الاحتمال المنفصل لتوليد $n$ من الإلكترونات الثانوية عن طريق تكامل توزيع بواسون فوق توزيع فقدان طاقة منتظم.
   • بعد الدمج مع سلسلة الدينودات اللاحقة، يشكل هذا المكون شكل "صندوق" مستدير موصوف بدوال الخطأ ($\text{erf}$).
   • هذا المكون هو الابتكار الجوهري، حيث يستبدل دوال الضجيج العشوائية باشتقاق قائم على الفيزياء.

3. النبضات المسبقة ($\gamma_{1Dy}$):

   • الفوتونات التي تتجاوز الكاثود الضوئي وتصطدم مباشرة بالدينود الأول.
   • تولد هذه الفوتونات إشارات ذات كسب (gain) أصغر بنسبة تقارب $1/G_1$ من الذروة الرئيسية.
   • يتم نمذجتها كتوزيع بواسون موزون بمعلمات كسب ودقة معدلة.

4. أحداث الشحنة المنخفضة الأسية ($\mathcal{P}_{Exp}$):

   • الإلكترونات ذات المسارات السيئة التي تخطئ الدينودات أو غيرها من حالات الشحنة المنخفضة.
   • يتم نمذجتها كـ توزيع غاوسي معدل أسيًا متمركز عند الأصل.

ب. التحقق التجريبي

تم اختبار النموذج باستخدام نموذجين مختلفين من أنابيب المضاعف الضوئي:

• Hamamatsu R5912-100 (بقطر 8 بوصات، مستخدم في تجارب XENONnT و LZ).
• Hamamatsu 6233 (بقط diameter 3 بوصات).

الإعداد:

• مصدر الضوء: ليزرات بيكو-ثانية (405 نانومتر و 507 نانومتر) بنبضات أقل من النانو-ثانية.
• القراءة: مضخم منخفض الضجيج (RFBAY LNA-150) ومحول رقمي عالي السرعة (2 جيجا عينة/ثانية).
• الظروف: أُجريت القياسات عند إعدادات جهد عالٍ (HV) مختلفة، ومواضع إضاءة مختلفة، وأطوال موجية مختلفة لاختبار استقرار المعلمات.
• معالجة البيانات: تم طرح العدات المظلمة (dark counts)، ودمج الموجات لتوليد أطياف الشحنة.

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

• أول اشتقاق تحليلي لإشارات الـ PA: تقدم الورقة أول دالة تحليلية صارمة لطيف الإلكترونات المتشتتة خلفيًا والمضخمة جزئيًا، مشتقة مباشرة من فيزياء التشتت الخلفي، وليس من مجرد ملاءمة تجريبية.
• إطار عمل موحد: يدمج النموذج ذرى التضخيم الكامل، وذيول التشتت الخلفي، والنبضات المسبقة، وذيول الشحنة المنخفضة في إطار رياضي واحد متسق مع عدد محدود من المعلمات الجوهرية.
• تقليل المعلمات: يقلل النموذج وصف طيف الـ SPE المعقد إلى عدد قليل من المعلمات ذات المعنى الفيزيائي:
  • $G_1$: كسب الدينود الأول.
  • $\eta$: كسر الإلكترونات المتشتتة خلفيًا.
  • $R$: معلمة الدقة لسلسلة الدينودات.
  • $\zeta$: نسبة الكسب الكامل مقابل متوسط الكسب عند الدينود الثاني (للنبضات المسبقة).
  • $\alpha$: ثابت الاضمحلال لأحداث الشحنة المنخفضة الأسية.

4. النتائج

تمت ملاءمة النموذج للبيانات التجريبية بدقة عالية:

• جودة الملاءمة: حقق النموذج توافقًا ممتازًا مع البيانات في جميع الظروف المختبرة، مع قيم مربع كاي المختزل ($\chi^2/\text{ndf}$) تراوحت بين 1.01 و 1.36.
• Hamamatsu R5912-100:
  • وُجد أن كسر التشتت الخلفي ($\eta$) يبلغ حوالي 27%، وهو ما يتوافق مع التوقعات النظرية للمعادن عالية العدد الذري (high-Z).
  • كان مكون الشحنة المنخفضة الأسي ($A_{Exp}$) مهمًا (~4%)، بينما كانت النبضات المسبقة مهملة.
  • أظهرت المعلمات $G_1$ وإجمالي الكسب اعتمادًا خطيًا متوقعًا على الجهد العالي (HV).
  • وُجد أن الكسر $\eta$ مستقل عن الـ HV ولكنه يعتمد على موضع الإضاءة على الكاثود الضوئي.
• Hamamatsu 6233:
  • أظهر كسر تشتت خلفي أقل ($\eta \approx 17\%$).
  • لوحظ وجود مساهمة نبضات مسبقة كبيرة ($A_{\gamma 1Dy} \approx 7\%$)، تم حلها بنجاح بواسطة النموذج.
  • تم استخراج معلمة الدقة $R$ بنجاح كمعلمة حرة نظرًا للفصل الواضح لذروة النبضة المسبقة.

5. الأهمية والتأثير

• المعايرة القائمة على الفيزياء: من خلال استبدال دوال "الضجيج" التجريبية بنموذج مشتق فيزيائيًا، يسمح هذا العمل بتحديد أكثر دقة لكسب الـ PMT، والكفاءة الكمية، والقبول الفعال في التجارب.
• تحسين المحاكاة: يمكن تطبيق الدوال التحليلية مباشرة في محاكاة مونت كارلو (مثل تجارب المادة المظلمة أو النيوترينو مثل XENONnT و LZ) لنمذجة الخلفية منخفضة الطاقة واستعادة الإشارة بشكل أفضل.
• أداة تشخيصية: تتيح القدرة على استخراج كسور محددة ($\eta$، $A_{Exp}$، $A_{\gamma 1Dy}$) للباحثين تشخيص سلوكيات محددة في الـ PMT، مثل خصائص سطح الدينود أو مشكلات المحاذاة البصرية، دون الاعتماد على الملاءمة "الصندوق الأسود".
• المتانة: تثبت صلاحية النموذج عبر أحجام PMT ومواد وظروف تشغيل مختلفة مدى عالميته لتطبيقات عد الفوتونات المنفردة.

في الختام، نجحت هذه الورقة في سد الفجوة بين فيزياء الإلكترونات النظرية وتوصيف الـ PMT التجريبي، مما يوفر أداة قوية وفعالة المعلمات للجيل القادم من تجارب الكشف عن الضوء المنخفض.