The "Intensity" Countoscope: Measuring particle dynamics in real space from microscopy images

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

يعد قياس ديناميكيات الجسيمات (مثل معاملات الانتشار والسلوك الجماعي) في الأنظمة المجهرية تحديًا جوهريًا في فيزياء المواد اللينة والبيولوجيا. تندرج أحدث الطرق الحالية عمومًا تحت ثلاث فئات:

1. الارتباط القائم على الشدة: تحليل تقلبات تشتت الضوء (مثل FCS وDLS وICS).
2. القائم على التجزئة (Segmentation): تحديد سمات جسيمات محددة (المواقع، الاتجاهات).
3. القائم على المسارات (Trajectory-based): ربط مواقع الجسيمات عبر الإطارات لحساب متوسط مربع الإزاحة (MSD).

محددات الطرق الحالية:

• الطرق القائمة على المسارات تتطلب تصويرًا عالي الدقة حيث يمكن تمييز الجسيمات الفردية بوضوح وتتبعها. وهي تفشل في الأنظمة الكثيفة أو عند الدقة المنخفضة حيث تتداخل الجسيمات أو تصبح غير قابلة للتمييز.
• طرق فضاء فوريه (مثل DDM أو k-space ICS) قوية للتحليل متعدد المقاييس ولكنها تعمل في فضاء التردد، مما يجعل الربط المباشر بالآليات الفيزيائية في الفضاء الحقيقي (مثل دخول/خروج الجسيمات من منطقة ما) أقل بديهية.
• الفجوة: هناك نقص في وجود طريقة قوية في الفضاء الحقيقي تستخدم تقلبات الشدة داخل صناديق ملاحظة افتراضية ذات أحجام عشوائية لاستخراج معاملات الانتشار، لا سيالما في الأنظمة التي لا يمكن فيها تمييز الجرسمات الفردية.

2. المنهجية

يقترح المؤلفون تقنية تحليل مبتكرة في الفضاء الحقيقي تسمى "Intensity Countoscope" (مجهر عد الشدة).

• النظام التجريبي: معلق مخفف من كولودات البوليستيرين الفلورية (قطر 4 ميكرومتر) محصورة في خلية مائعية ثنائية الأبعاد، يتم تصويرها عبر مجهر مسح ليزري بؤري (CLSM).
• صناديق الملاحظة الافتراضية: بدلًا من تحليل مجال الرؤية الكامل أو تتبع الجسيمات الفردية، تقسم الطريقة صورة المجهر إلى صناديق مربعة افتراضية بحجم $L \times L$.
• استخراج الإشارة: لكل صندوق، يتم حساب إجمالي الشدة $I(t)$ عن طريق جمع شدة البكسلات. تحلل الطريقة متوسط تغير مربع الشدة بمرور الوقت:
  $$\langle \Delta I^2(t) \rangle = \langle (I(t) - I(0))^2 \rangle = 2\langle I^2 \rangle - 2\langle I(t)I(0) \rangle$$
  هذا المقياس يزيل تلقائيًا ضوضاء الخلفية الثابتة والجسيمات "العالقة".
• الإطار النظري:
  • اشتق المؤلفون نظرية عامة تربط تقلبات الشدة بتوزيع كثافة الجسيمات $\hat{\rho}$ ودالة انتشار النقطة (PSF) للمجهر.
  • باستخدام تحويلات فوريه، قاموا بربط ارتباط الشدة بدالة التشتت المتوسطة.
  • قاموا بنمذجة الـ PSF كدالة غاوسية (نموذج شائع في المجهر البؤري) وصندوق الملاحظة كدالة مؤشر مربعة.
  • التعبير التحليلي الناتج عن $\langle \Delta I^2(t) \rangle$ يعتمد على معامل الانتشار $D$، ومتوسط شدة الجسيم $I_0$، وعرض الجسيم الفعال $\sigma$، وحجم الصندوق $L$.

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

أ. اكتشاف أنظمة زمنية متميزة

يكشف التحليل أن قياس تقلبات الشدة يعتمد بشكل حاسم على نسبة حجم صندوق الملاحظة ($L$) إلى عرض الجسيم الفعال ($\sigma$):

1. الصناديق الصغيرة / الأطر الزمنية القصيرة ($\sigma \gg L$ أو وقت مبكر $t$):
   • يكون الصندوق أصغر من PSF الخاص بالجسيم.
   • يتناسب $\langle \Delta I^2(t) \rangle$ خطيًا مع متوسط مربع الإزاحة (MSD)، أي أن $\langle \Delta I^2(t) \rangle \propto \langle \Delta x^2(t) \rangle \propto t$.
   • التفسير الفيزيائي: الصندوق "ينتشر" فعليًا عبر ملف شدة الجسيم.
2. الصناديق الكبيرة / الأطر الزمنية الطويلة ($\sigma \ll L$ أو وقت متأخر $t$):
   • يكون الصندوق أكبر من الجسيم.
   • يتناسب $\langle \Delta I^2(t) \rangle$ مع الجذر التربيعي لـ MSD، أي أن $\langle \Delta I^2(t) \rangle \propto \sqrt{\langle \Delta x^2(t) \rangle} \propto \sqrt{t}$.
   • التفسير الفيزيائي: تهيمن التقلبات الناتجة عن الدخول والخروج العشوائي لجسيمات كاملة من وإلى حجم الملاحظة.

ب. استخراج المعلمات بمتانة

تسمح الطريقة باستخراج المعلمات الفيزيائية في آن واحد دون الحاجة لتمييز الجسيمات الفردية:

• معامل الانتشار ($D$): يتم استخراجه عن طريق مطابقة المنحنى المعتمد على الزمن.
• عرض الجسيم ($\sigma$) والشدة ($I_0$): يتم استخراجهما من خلال تحليل قيمة الهضبة (التباين) لتقلبات الشدة عند الأوقات الطويلة.

ج. الاستقلال عن الدقة

أظهر المؤلفون أن الطريقة تظل قوية حتى لو تم تقليل دقة الصورة اصطناعيًا (تقليل العينات). حتى عندما تصبح الجسيمات الفردية غير قابلة للتمييز (حجم البكسل > حجم الجسيم)، يمكن بعد ذلك تحديد معامل الانتشار بدقة، بشرط أن تأخذ عملية المطابقة فقدان نظام $\sigma \gg L$ في الاعتبار.

4. النتائج

• التحقق من الصحة: يتطابق النموذج النظري (المعادلة 6 في الورقة) تمامًا مع البيانات التجريبية عبر جميع الأطر الزمنية وأحجام الصناديق لمعلق كولودات مخفف.
• الدقة الكمية:
  • معامل الانتشار المستخرج عبر الـ Intensity Countoscope كان $D = 0.076 \, \mu m^2/s$.
  • وهذا يطابق القيمة التي تم الحصول عليها من تتبع الجسيمات القياسي ($D = 0.078 \, \mu m^2/s$) ويتوافق مع النظرية الهيدروديناميكية للجسيمات بالقرب من جدار (تنبؤ ستوكس-أينشتاين كان $0.107 \, \mu m^2/s$ في الكتلة، وانخفض بسبب احتكاك الجدار).
  • عرض الجسيم الذي تم استخراجه ($\sigma \approx 1.05 \, \mu m$) يتوافق مع نصف قطر الجسيم الفيزيائي ($2 \, \mu m$) بالنظر إلى PSF البصري.
• الأداء عند الدقة المنخفضة: عندما تم تقليل دقة الصورة بعوامل 16 و64 (مما جعل الجسيمات غير قابلة للتمييز)، ظلت الطريقة تعطي معاملات انتشار ضمن نطاق 10-15% من القيمة الحقيقية. تحسنت الدقة عندما تم تثبيت عرض الجسيم $\sigma$ بناءً على معرفة مسبقة.

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

• سد الفجوة: تسد هذه الطريقة الفجوة بين تحليل الارتباط القائم على الشدة (مثل FCS) وتتبع الجسيمات في الفضاء الحقيقي. وهي توسع مفهوم "Countoscope" (الذي كان مقتصرًا سابقًا على عد الجسيمات) ليشمل إشارات الشدة.
• القابلية للتطبيق في الأنظمة غير القابلة للتمييز: تتيح دراسة الديناميكيات في الأنظمة التي لا يمكن تجزئة الجسيمات فيها (مثل المعلقات الكثيفة، أو الأنسجة البيولوجية، أو التصوير منخفض الدقة)، مما يوفر أداة لدراسة كل من السلوك الفردي والجماعي.
• الكفاءة الحسابية: بخلاف الطرق القائمة على فوريه (DDM)، تعمل هذه الطريقة بالكامل في الفضاء الحقيقي، مما يتجنب تحويلات فوريه المكلفة حسابيًا ويجعل من السهل تفسير الآليات الفيزيائية مباشرة.
• تعدد الاستخدامات: الإطار النظري عام ويمكن تكييفه مع أشكال صناديق مختلفة (مثل التوزيع الغاوسي لـ FCS) وكذلك مع نماذج PSF مختلفة (مثل أقراص Airy للمجهر ذي المجال المضيء)، مما يجعلها قابلة للتطبيق على مجموعة واسعة من تقنيات المجهر (TIRF، المجال المضيء، إلخ).

باختصار، يوفر "Intensity Countoscope" أداة قوية ومتينة وفعالة حسابيًا لتحديد ديناميكيات الجسيمات مباشرة من شدة الصور الخام، متجاوزًا قيود الدقة التي تعيق تحليل المسارات التقليدي.