Significance: Changes in the structure of tissue occur in many disease processes, such as the boundaries of cancerous tumors and burn injuries. Spectroscopic and polarimetric alterations of terahertz light caused by Mie scattering patterns has the potential to be a diagnostic marker.
Aim: We present an analysis of Monte Carlo simulation of Mie scattering of polarized terahertz light from cancerous tumor budding, compare the simulation to experimental results obtained in phantom models, and present an analysis of a polarization-sensitive terahertz scan of an ex vivo porcine burn injury.
Approach: Using a Monte Carlo simulation, we modeled the changes in diffuse intensity and degree of polarization of broadband off-specular terahertz light due to scattering particles in highly attenuating tissue. We extracted the Mueller matrix of the tissue using this model and analyzed the Lu-Chipman product decomposition matrices. We compared this model to experimental data from four phantoms consisting of polypropylene particles of varying sizes embedded in gelatin. Finally, we induced a full-thickness burn injury in ex vivo porcine skin samples and compared experimental data from burned and healthy regions of the tissue.
Results: Simulation revealed contrast in the Stokes vectors and Mueller Matrix elements for varying scattering particle sizes. Experimental phantom results showed contrast between different sizes of scattering particles in degree of polarization and diffuse intensity in agreement with Monte Carlo simulation results. Finally, we demonstrated a similar diffused imaging signal contrast between burned and healthy regions of ex vivo porcine skin.
Conclusion: Polarimetric terahertz imaging has the potential to detect structural changes due to biological disease processes.
- معرّف الورقة: 2501.00630
- العنوان: Mie scattering due to tissue structures in the terahertz regime: Experimental and Monte Carlo verification using diffused polarimetric imaging in highly attenuating tissue phantoms
- المؤلفون: إيريكا هيلر، كوانغيي شو، زاكري هاريس، م. حسن أربب (جامعة ستوني بروك)
- التصنيف: physics.med-ph physics.bio-ph physics.optics
- تاريخ الإرسال: 31 ديسمبر 2024
- رابط الورقة: https://arxiv.org/abs/2501.00630
تستكشف هذه الدراسة ظاهرة تشتت مي الناجمة عن هياكل الأنسجة في نطاق التيراهيرتز، مع التركيز على تطبيقاتها المحتملة في التشخيص الطبي. تحلل الدراسة خصائص التشتت للضوء المستقطب في التيراهيرتز في الأنسجة عالية التوهين من خلال محاكاة مونت كارلو والتحقق التجريبي. كشفت النتائج عن تباينات واضحة بين جزيئات التشتت ذات الأحجام المختلفة من حيث درجة الاستقطاب (DOP) وشدة التشتت المنتشر، وتم التحقق من الإمكانات التشخيصية للتقنية في عينات حروق الجلد الخنزيري السابقة للحيوية.
التغييرات في هياكل الأنسجة تمثل سمات مهمة في العديد من عمليات المرض، مثل براعم الأورام (tumor budding) عند حدود الأورام السرطانية وإصابات الحروق. تركز التطبيقات الطبية الحيوية التقليدية للتيراهيرتز بشكل أساسي على الخصائص العازلة للأنسجة، خاصة محتوى الماء، مع الاستفادة القليلة من خصائص التشتت للضوء في التيراهيرتز كوسيلة تشخيصية.
- احتياجات التشخيص الطبي: براعم الأورام والعناقيد الخلوية سيئة التمايز تمثل عوامل تنبؤية مستقلة لانتقال العقد اللمفاوية وبقاء المريض
- القيود التكنولوجية: الطرق الحالية مثل صبغة H&E والكيمياء النسيجية المناعية تتطلب شرائح نسيجية، مما يحول دون التشخيص في الوقت الفعلي
- آليات التشتت غير المفهومة: لم تتم دراسة آليات الاستجابة الاستقطابية في نطاق التيراهيرتز بشكل كافٍ
مع تطور مقاييس الاستقطاب في التيراهيرتز وأنظمة القياس الاستقطابية، هناك حاجة ملحة لنماذج حسابية لتفسير الظواهر المرصودة، مما يوفر دعماً نظرياً لتطبيقات التيراهيرتز في الفوتونيات الحيوية.
- إنشاء نموذج مونت كارلو لتشتت مي في التيراهيرتز: بناءً على النموذج الذي طوره شو وأربب، مع تحسينات لجزيئات التشتت في الوسائط الماصة
- التحقق من الإمكانات التشخيصية للقياسات الاستقطابية: إثبات التباين بين جزيئات التشتت ذات الأحجام المختلفة من خلال تجارب الأنسجة الشبحية
- توفير تحليل تحلل مصفوفة مويلر: استخدام طريقة تحلل لو-تشيبمان لتحليل الخصائص الاستقطابية للأنسجة
- عرض التطبيقات العملية: التحقق من جدوى التقنية في عينات حروق الجلد الخنزيري السابقة للحيوية
- تبسيط متطلبات القياس: إثبات أنه بالنسبة للجزيئات الكروية، يكفي قياس حالة استقطاب واحدة فقط لبناء مصفوفة مويلر كاملة
نظام التيراهيرتز في المجال الزمني:
- طرف الإرسال: استخدام ليزر فيمتوثانية بطول موجة 1560 نانومتر لتحفيز هوائي موصل الضوء (PCA)
- طرف الاستقبال: ذراع كشف قابلة للدوران من 0-160 درجة، تتضمن مستقطب شبكة سلكية (WGP)
- نطاق التردد: 0.2-1.5 تيراهيرتز
تكوين قياس التشتت:
- اختيار زاوية 140 درجة لقياس التشتت المنتشر (أوضح تباين بين القوة المتماسكة وغير المتماسكة)
- مسح الشبكة: فاصل بكسل 2 ملم، منطقة 2 سم × 2 سم
- تسجيل 100 قياس لكل حالة استقطاب
اختيار المواد:
- المصفوفة: محلول جيلاتين بقري بنسبة 15% (محاكاة خلفية الأنسجة)
- جزيئات التشتت: حبيبات البولي بروبيلين (محاكاة الهياكل منخفضة معامل الانكسار مثل الجريبات والغدد العرقية)
- نطاق حجم الجزيئات: 115-280 ميكرومتر (يقابل معاملات تشتت مي من 0.2-20)
عملية التحضير:
- تسخين محلول الجيلاتين إلى 40°C
- خلط سريع مع حبيبات البولي بروبيلين
- صب في طبق الثقافة للتصلب (سمك حوالي 5 ملم)
- تحليل توزيع حجم الجزيئات والتركيز باستخدام المجهر البصري وبرنامج ImageJ
الخوارزمية الأساسية:
بناءً على رمز مونت كارلو للضوء المستقطب في المستوى الزوالي الذي طوره رامييلا-رومان وآخرون، مع تعديلات لتشتت مي في الوسائط الماصة.
حسابات المعاملات الرئيسية:
- معامل التشتت: μs=∣nm∣2k022πρ∑j=1∞(2j+1)(∣aj∣2+∣bj∣2)
- معامل الامتصاص: μa=ℜ(nm)k022πρℑ[(∑j=1∞(2j+1)(∣cj∣2ψj(z)ψj′∗(z)−∣dj∣2ψj′(z)ψj∗(z)))/np]
بناء مصفوفة مويلر:
حساب معاملات ستوكس المكانية باستخدام أربع حالات استقطاب مختلفة للإدخال (استقطاب خطي أفقي، استقطاب خطي عمودي، استقطاب خطي بزاوية 45 درجة، استقطاب دائري يميني)، ثم بناء مصفوفة مويلر.
تم تحضير أربعة أشباح أنسجة مختلفة:
| الشبح | حجم الجزيئات (ميكرومتر) | الكثافة السطحية (جزيء/ملم²) | الكثافة الحجمية (ملم⁻³) |
|---|
| A | 115 | 3 | 10 |
| B | 130 | 2 | 7 |
| C | 180 | 3.5 | 8.5 |
| D | 280 | 6.5 | 10.5 |
الخصائص العازلة:
- البولي بروبيلين: معامل انكسار ~1.51، معامل امتصاص <2 سم⁻¹
- الجيلاتين: معامل انكسار ومعامل امتصاص بين الجلد والماء
- العينة: جلد الخنزير السابق للحيوية
- تحضير الحروق: قضيب نحاسي بدرجة حرارة 180°C لمدة 5 دقائق، مما يشكل حرقاً كاملاً بحجم 1 بوصة
- منطقة القياس: 2 × 3.5 سم، حجم بكسل 1 ملم
- جمع البيانات: 10 قياسات لكل بكسل، 10 مرات متوسط زمني
كفاءة التشتت ومعامل عدم التناسق:
مع زيادة حجم الجزيئات، تزداد كفاءة التشتت ومعامل عدم التناسق في النطاق منخفض التردد (<1 تيراهيرتز)، مما يشير إلى أن الجزيئات الأكبر تنتج تشتتاً أمامياً أقوى.
تحلل مصفوفة مويلر:
تحلل مصفوفة مويلر للشبح A و D باستخدام تحلل لو-تشيبمان:
- الشبح A: MΔ,A≈diag[1,0.921,0.925,0.849]
- الشبح D: MΔ,D≈diag[1,0.586,0.586,0.207]
تشير النتائج إلى أن كلا الشبحين يمكن تقريبهما كمزيج من مخفف موحد وإزالة استقطاب غير موحدة.
شدة التشتت المنتشر:
- أظهرت المحاكاة والتجربة: كلما زاد حجم الجزيئات، زادت شدة التشتت
- الاعتماد على التردد: الشدة منخفضة التردد عالية، تتحلل بشكل أسي مع التردد
- اتساق الاتجاه: تتطابق اتجاهات المحاكاة والنتائج التجريبية
تباين درجة الاستقطاب (DOP):
- ينخفض DOP للشبح D إلى حوالي 0.6 بالقرب من 0.6 تيراهيرتز
- ينخفض DOP للشبح C بدرجة متوسطة بين A و B و D
- تحديد التردد: يرتبط موقع التردد لانخفاض DOP بحجم الجزيئات
التباين المكاني:
في تردد 0.6 تيراهيرتز، تظهر منطقة الحروق والأنسجة السليمة تباينات واضحة في DOP وشدة التشتت المنتشر.
الخصائص الطيفية:
- ينخفض DOP للأنسجة السليمة بشكل أسرع في نطاق 0.2-0.8 تيراهيرتز
- تكون شدة التشتت المنتشر لمنطقة الحروق أعلى بكثير من الأنسجة السليمة
- قد يعزى هذا إلى تدمير الهياكل الجلدية الكبيرة (الجريبات والغدد العرقية) بسبب الحروق وفقدان الرطوبة
تركز الأبحاث الموجودة على:
- التشخيص السرطاني: الدراسات الطيفية لسرطان الثدي وسرطان الخلايا الحرشفية الفموي السابقة للحيوية
- تصنيف حروق الجلد: تصنيف الحروق في الجسم الحي والتنبؤ بالشفاء
- تحليل الخصائص العازلة: التركيز الأساسي على تغييرات محتوى الماء في الأنسجة
تطبيقات واسعة في نطاقات أخرى:
- التصوير بالمنظار الداخلي: كشف الاورام الحميدة في القولون وآفات الصفاق
- التشخيص السرطاني: سرطان الجلد والمعدة والثدي وأنواع سرطانية متعددة أخرى
- نمذجة مونت كارلو: تطبيق نماذج التشتت الضوئي في الفوتونيات الحيوية
تتأخر أبحاث آليات الاستجابة الاستقطابية في نطاق التيراهيرتز، مع نقص في النماذج الحسابية الناضجة لتفسير الظواهر المرصودة.
- نجاح التحقق من النموذج: تتطابق نتائج محاكاة مونت كارلو مع البيانات التجريبية بشكل كبير من حيث الاتجاه
- تأكيد الإمكانات التشخيصية: يمكن لـ DOP وشدة التشتت المنتشر التمييز بين جزيئات التشتت ذات الأحجام المختلفة
- تبسيط القياس: بالنسبة للجزيئات الكروية، يمكن الحصول على معلومات مصفوفة مويلر الكاملة من قياس حالة استقطاب واحدة
- الجدوى السريرية: تم إثبات القدرة على كشف الحالات المرضية بنجاح في الأنسجة السابقة للحيوية
- قيود نطاق النظام: يحد النطاق الحالي 0.2-1.5 تيراهيرتز من القدرة على تمييز الفروقات الصغيرة في حجم الجزيئات
- قيود نسبة الإشارة إلى الضوضاء: تنخفض موثوقية قياسات DOP فوق 0.8 تيراهيرتز
- تأثير خشونة السطح: يفترض أن سطح العينة أملس ومسطح، مما يتطلب في التطبيقات العملية الأخذ في الاعتبار تأثيرات إزالة الاستقطاب السطحي
- متطلبات الفروقات في حجم الجزيئات: يتطلب فروقات كبيرة في حجم الجزيئات لإنتاج تباين إشارة واضح
- توسيع نطاق النظام: استخدام بلورات كهروضوئية بدلاً من PCA لتحقيق نطاق 8 تيراهيرتز
- التحقق من أنسجة السرطان السابقة للحيوية: توسيع دراسة الشبح إلى أنسجة سرطانية حقيقية
- تحسين معالجة الإشارات: تطوير تقنيات تحلل موجة الإشارة وتحليل البقع الإحصائي
- فصل مصادر إزالة الاستقطاب المتعددة: التمييز بين المساهمات النسبية لخشونة السطح والهياكل الداخلية
- دمج النظرية والتجربة: دمج محاكاة مونت كارلو مع التحقق التجريبي يوفر إطار عمل بحثي شامل
- الابتكار المنهجي: أول دراسة منهجية لخصائص الاستقطاب لتشتت مي في نطاق التيراهيرتز
- قيمة عملية عالية: تبسيط متطلبات قياس مصفوفة مويلر يحسن جدوى التطبيق السريري
- التحقق الكافي: استراتيجية التحقق التدريجي من الأشباح المثالية إلى الأنسجة الحية الحقيقية معقولة
- التحليل العميق: يوفر تحلل لو-تشيبمان فهماً أعمق للآليات الفيزيائية
- تنوع العينات محدود: استخدام أشباح البولي بروبيلين والجيلاتين فقط قد تختلف خصائص المواد عن الأنسجة الحقيقية
- التحليل الإحصائي غير كافٍ: نقص اختبارات الدلالة الإحصائية التفصيلية وتحليل الأخطاء
- التحقق السريري محدود: تم التحقق فقط على عينة حروق واحدة، مما يتطلب دراسات سريرية أكبر نطاقاً
- عدم كفاية اعتبار التشتت المعقد: التركيز الأساسي على الجزيئات الكروية، بينما هياكل الأنسجة الحقيقية أكثر تعقيداً
- القيمة الأكاديمية: توفير أساس نظري وطرق تجريبية مهمة لفوتونيات التيراهيرتز الحيوية
- التقدم التكنولوجي: تعزيز تطبيق تقنية التصوير الاستقطابي في التيراهيرتز في التشخيص الطبي
- الأهمية متعددة التخصصات: ربط نظرية التشتت الكهرومغناطيسي ونمذجة مونت كارلو والتصوير الطبي الحيوي
- الآفاق الصناعية: توفير مسار تكنولوجي لتطوير أجهزة التيراهيرتز الطبية
- تشخيص أمراض الجلد: تصنيف الحروق وكشف سرطان الجلد
- تحديد حدود الأورام: الملاحة الجراحية وتحليل الأمراض
- توصيف هياكل الأنسجة: تحليل ألياف الكولاجين والعناقيد الخلوية
- مراقبة فعالية الأدوية: تتبع التغييرات في هياكل الأنسجة أثناء العلاج
تستشهد الورقة بـ 41 مرجعاً ذا صلة، تغطي أعمالاً مهمة في مجالات متعددة بما في ذلك تقنية التيراهيرتز والبصريات الطبية الحيوية ونمذجة مونت كارلو، مما يوفر أساساً نظرياً متيناً للبحث.
التقييم الشامل: هذه ورقة بحثية عالية الجودة متعددة التخصصات ذات قيمة نظرية وعملية مهمة في مجال التصوير الطبي الحيوي للتيراهيرتز. تتسم طرق البحث بالدقة والنتائج موثوقة، مما يضع أساساً جيداً للتطور المستقبلي في هذا المجال.