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.
- पेपर ID: 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
- लेखक: Erica Heller, Kuangyi Xu, Zachery Harris, M. Hassan Arbab (Stony Brook University)
- वर्गीकरण: physics.med-ph physics.bio-ph physics.optics
- प्रस्तुति समय: 31 दिसंबर 2024
- पेपर लिंक: https://arxiv.org/abs/2501.00630
यह अध्ययन टेराहर्ट्ज़ आवृत्ति में ऊतक संरचनाओं के कारण Mie बिखराव घटना की जांच करता है, जिसमें रोग निदान में संभावित अनुप्रयोग पर ध्यान केंद्रित किया गया है। अनुसंधान Monte Carlo सिमुलेशन और प्रायोगिक सत्यापन के माध्यम से, उच्च क्षीणकारी ऊतक में ध्रुवीकृत टेराहर्ट्ज़ प्रकाश की बिखराव विशेषताओं का विश्लेषण किया गया है। अध्ययन में पाया गया कि विभिन्न आकार के बिखराव कणों ने ध्रुवीकरण की डिग्री (DOP) और विसरित बिखराव तीव्रता में स्पष्ट विपरीतता प्रदर्शित की, और पूर्व में सूअर की त्वचा के जलने के नमूनों में इस तकनीक की निदान क्षमता को सत्यापित किया गया।
ऊतक संरचना में परिवर्तन कई रोग प्रक्रियाओं की महत्वपूर्ण विशेषता है, जैसे कि कैंसर ट्यूमर सीमा में ट्यूमर कलियों (tumor budding) और जलने की चोट। पारंपरिक टेराहर्ट्ज़ जैव चिकित्सा अनुप्रयोग मुख्य रूप से ऊतक के परावैद्युत गुणों पर ध्यान केंद्रित करते हैं, विशेष रूप से जल सामग्री, और टेराहर्ट्ज़ प्रकाश की बिखराव विशेषताओं को निदान साधन के रूप में शायद ही कभी उपयोग करते हैं।
- रोग निदान की आवश्यकता: ट्यूमर कलियां और विभेदित कोशिकाओं के समूह लिम्फ नोड मेटास्टेसिस और रोगी जीवन के स्वतंत्र पूर्वानुमानात्मक कारक हैं
- तकनीकी सीमाएं: H&E धुंधलापन और इम्यूनोहिस्टोकेमिस्ट्री जैसी मौजूदा विधियों को ऊतक स्लाइस की आवश्यकता होती है, वास्तविक समय निदान को सक्षम नहीं कर सकते
- बिखराव तंत्र अस्पष्ट: टेराहर्ट्ज़ आवृत्ति में ध्रुवीकरण प्रतिक्रिया अंतर तंत्र को पर्याप्त रूप से अध्ययन नहीं किया गया है
टेराहर्ट्ज़ एलिप्सोमेट्री और ध्रुवीकरण माप प्रणाली के विकास के साथ, टेराहर्ट्ज़ बायोफोटोनिक्स अनुप्रयोगों के लिए सैद्धांतिक समर्थन प्रदान करने के लिए देखी गई घटनाओं को समझाने के लिए कम्प्यूटेशनल मॉडल की तत्काल आवश्यकता है।
- टेराहर्ट्ज़ Mie बिखराव का Monte Carlo मॉडल स्थापित किया: Xu और Arbab द्वारा विकसित मॉडल के आधार पर, अवशोषक माध्यम में बिखराव कणों के लिए सुधार किया गया
- ध्रुवीकरण माप की निदान क्षमता को सत्यापित किया: ऊतक फैंटम प्रयोगों के माध्यम से विभिन्न आकार के बिखराव कणों की विपरीतता को प्रदर्शित किया
- Mueller मैट्रिक्स अपघटन विश्लेषण प्रदान किया: Lu-Chipman अपघटन विधि का उपयोग करके ऊतक की ध्रुवीकरण विशेषताओं का विश्लेषण किया
- व्यावहारिक अनुप्रयोग प्रदर्शित किया: पूर्व में सूअर की त्वचा के जलने के नमूनों में तकनीक की व्यवहार्यता को सत्यापित किया
- माप आवश्यकताओं को सरल बनाया: गोलाकार कणों के लिए, केवल एक ध्रुवीकरण स्थिति माप पूर्ण Mueller मैट्रिक्स का निर्माण कर सकता है
टेराहर्ट्ज़ समय-डोमेन स्पेक्ट्रोस्कोपी प्रणाली:
- संचरण अंत: 1560 nm फेमटोसेकंड लेजर द्वारा उत्तेजित फोटोकंडक्टिव एंटीना (PCA) का उपयोग
- रिसीवर अंत: 0-160 डिग्री घुमाने योग्य पहचान भुजा, जिसमें तार ग्रिड ध्रुवीकरण (WGP) शामिल है
- आवृत्ति श्रेणी: 0.2-1.5 THz
बिखराव माप विन्यास:
- 140 डिग्री कोण पर विसरित बिखराव माप के लिए चयन (सुसंगत शक्ति और गैर-सुसंगत शक्ति विपरीतता सबसे स्पष्ट)
- ग्रिड स्कैन: 2mm पिक्सेल रिक्ति, 2cm × 2cm क्षेत्र
- प्रत्येक ध्रुवीकरण स्थिति के लिए 100 माप दर्ज किए गए
सामग्री का चयन:
- मैट्रिक्स: 15% गोजातीय त्वचा जिलेटिन समाधान (ऊतक पृष्ठभूमि का अनुकरण)
- बिखराव कण: पॉलीप्रोपाइलीन कणिकाएं (बालों के रोम, पसीने की ग्रंथियों आदि कम अपवर्तक सूचकांक संरचनाओं का अनुकरण)
- कण आकार श्रेणी: 115-280 μm (Mie बिखराव पैरामीटर 0.2-20 के अनुरूप)
तैयारी प्रक्रिया:
- जिलेटिन समाधान को 40°C तक गर्म किया
- पॉलीप्रोपाइलीन कणिकाओं को तेजी से मिलाया
- संस्कृति पकवान में डाला और ठोस किया (मोटाई लगभग 5mm)
- ऑप्टिकल माइक्रोस्कोप और ImageJ सॉफ्टवेयर का उपयोग करके कण आकार वितरण और एकाग्रता का विश्लेषण किया
मुख्य एल्गोरिथ्म:
Ramella-Roman आदि द्वारा मेरिडियनल समतल ध्रुवीकृत प्रकाश Monte Carlo कोड के आधार पर, अवशोषक माध्यम में Mie बिखराव के लिए संशोधन किया गया।
मुख्य पैरामीटर गणना:
- बिखराव गुणांक: μ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]
Mueller मैट्रिक्स निर्माण:
चार विभिन्न आपतन ध्रुवीकरण अवस्थाओं (क्षैतिज रैखिक ध्रुवीकरण, ऊर्ध्वाधर रैखिक ध्रुवीकरण, 45° रैखिक ध्रुवीकरण, दाहिने हाथ की गोलाकार ध्रुवीकरण) का उपयोग करके स्थानिक Stokes पैरामीटर की गणना की गई, और फिर Mueller मैट्रिक्स का निर्माण किया गया।
चार विभिन्न ऊतक फैंटम तैयार किए गए:
| फैंटम | कण आकार (μm) | सतह घनत्व (कण/mm²) | आयतन घनत्व (mm⁻³) |
|---|
| A | 115 | 3 | 10 |
| B | 130 | 2 | 7 |
| C | 180 | 3.5 | 8.5 |
| D | 280 | 6.5 | 10.5 |
परावैद्युत गुण:
- पॉलीप्रोपाइलीन: अपवर्तक सूचकांक ~1.51, अवशोषण गुणांक <2 cm⁻¹
- जिलेटिन: अपवर्तक सूचकांक और अवशोषण गुणांक त्वचा और पानी के बीच
- नमूना: पूर्व में सूअर की त्वचा
- जलने की तैयारी: 180°C पीतल की छड़ 5 मिनट के लिए, 1 इंच पूर्ण परत जलने का निर्माण
- माप क्षेत्र: 2×3.5 cm, 1mm पिक्सेल आकार
- डेटा संग्रह: प्रति पिक्सेल 10 माप, 10 बार समय औसत
बिखराव दक्षता और असमरूपता कारक:
कण आकार बढ़ने के साथ, बिखराव दक्षता और असमरूपता कारक कम आवृत्ति खंड (<1 THz) में दोनों बढ़ते हैं, जो दर्शाता है कि बड़े कण अधिक मजबूत अग्रगामी बिखराव उत्पन्न करते हैं।
Mueller मैट्रिक्स अपघटन:
फैंटम A और D के Mueller मैट्रिक्स के लिए Lu-Chipman अपघटन:
- फैंटम A: MΔ,A≈diag[1,0.921,0.925,0.849]
- फैंटम D: MΔ,D≈diag[1,0.586,0.586,0.207]
परिणाम दर्शाते हैं कि दोनों फैंटम को समान क्षीणकारी और गैर-समान depolarizer के संयोजन के रूप में अनुमानित किया जा सकता है।
विसरित बिखराव तीव्रता:
- सिमुलेशन और प्रयोग दोनों दिखाते हैं: कण आकार जितना बड़ा, बिखराव तीव्रता उतनी अधिक
- आवृत्ति निर्भरता: कम आवृत्ति तीव्रता अधिक, आवृत्ति के साथ घातांकीय क्षय
- प्रवृत्ति संगति: सिमुलेशन और प्रायोगिक परिणाम प्रवृत्ति में सुसंगत
ध्रुवीकरण की डिग्री (DOP) विपरीतता:
- फैंटम D 0.6 THz के पास DOP को लगभग 0.6 तक गिराता है
- फैंटम C की DOP में गिरावट A, B और D के बीच मध्यवर्ती है
- आवृत्ति स्थिति: DOP गिरावट की आवृत्ति स्थिति कण आकार से संबंधित है
स्थानिक विपरीतता:
0.6 THz आवृत्ति पर, जलने वाले क्षेत्र और स्वस्थ ऊतक DOP और विसरित बिखराव तीव्रता दोनों में स्पष्ट विपरीतता दिखाते हैं।
वर्णक्रमीय विशेषताएं:
- स्वस्थ ऊतक का DOP 0.2-0.8 THz श्रेणी में तेजी से गिरता है
- जलने वाले क्षेत्र की विसरित बिखराव तीव्रता स्वस्थ ऊतक से काफी अधिक है
- यह संभवतः जलने के कारण बड़ी त्वचा संरचनाओं (बालों के रोम, पसीने की ग्रंथियों) के विनाश और जल हानि के कारण है
मौजूदा अनुसंधान मुख्य रूप से निम्नलिखित पर केंद्रित है:
- कैंसर निदान: पूर्व में स्तन कैंसर, मुखीय स्क्वामस सेल कार्सिनोमा आदि का स्पेक्ट्रोस्कोपी अध्ययन
- त्वचा जलने का वर्गीकरण: शरीर में जलने की ग्रेडिंग और उपचार पूर्वानुमान
- परावैद्युत गुण विश्लेषण: मुख्य रूप से ऊतक जल सामग्री परिवर्तन पर ध्यान केंद्रित
अन्य तरंग दैर्ध्य में व्यापक अनुप्रयोग:
- एंडोस्कोपिक इमेजिंग: कोलन पॉलिप्स और पेरिटोनियल रोग का पता लगाना
- कैंसर निदान: त्वचा कैंसर, पेट कैंसर, स्तन कैंसर आदि कई कैंसर उप-प्रकार
- Monte Carlo मॉडलिंग: बायोफोटोनिक्स में प्रकाश बिखराव मॉडल का अनुप्रयोग
टेराहर्ट्ज़ आवृत्ति में ध्रुवीकरण प्रतिक्रिया तंत्र अनुसंधान अपेक्षाकृत पिछड़ा है, देखी गई घटनाओं को समझाने के लिए परिपक्व कम्प्यूटेशनल मॉडल की कमी है।
- मॉडल सत्यापन सफल: Monte Carlo सिमुलेशन परिणाम प्रायोगिक डेटा के साथ प्रवृत्ति में अत्यधिक सुसंगत
- निदान क्षमता की पुष्टि: DOP और विसरित बिखराव तीव्रता विभिन्न आकार के बिखराव कणों को अलग कर सकते हैं
- माप सरलीकरण: गोलाकार कणों के लिए, एकल ध्रुवीकरण स्थिति माप पूर्ण Mueller मैट्रिक्स जानकारी प्राप्त कर सकता है
- नैदानिक व्यवहार्यता: पूर्व में ऊतक में रोग संबंधी स्थिति का पता लगाने की क्षमता को सफलतापूर्वक प्रदर्शित किया
- प्रणाली बैंडविड्थ सीमा: वर्तमान 0.2-1.5 THz बैंडविड्थ छोटे कण आकार अंतर के विभेदन क्षमता को सीमित करता है
- संकेत-से-शोर अनुपात बाधा: 0.8 THz से ऊपर DOP माप विश्वसनीयता में कमी
- सतह खुरदरापन प्रभाव: नमूना सतह को चिकना और समतल मानता है, वास्तविक अनुप्रयोग में सतह depolarization प्रभाव पर विचार करने की आवश्यकता है
- कण आकार अंतर आवश्यकता: स्पष्ट संकेत विपरीतता उत्पन्न करने के लिए बड़े कण आकार अंतर की आवश्यकता है
- प्रणाली बैंडविड्थ विस्तार: इलेक्ट्रो-ऑप्टिक क्रिस्टल का उपयोग करके PCA को प्रतिस्थापित करना, 8 THz बैंडविड्थ प्राप्त करना
- पूर्व में कैंसर ऊतक सत्यापन: फैंटम अध्ययन को वास्तविक कैंसर ऊतक तक विस्तारित करना
- संकेत प्रसंस्करण अनुकूलन: तरंग संकेत अपघटन और सांख्यिकीय speckle विश्लेषण तकनीकें विकसित करना
- बहु-स्रोत depolarization अलगाव: सतह खुरदरापन और आंतरिक संरचना के सापेक्ष योगदान को अलग करना
- सिद्धांत और प्रयोग का संयोजन: Monte Carlo सिमुलेशन और प्रायोगिक सत्यापन का संयोजन, एक पूर्ण अनुसंधान ढांचा प्रदान करता है
- पद्धति नवाचार: टेराहर्ट्ज़ आवृत्ति में Mie बिखराव की ध्रुवीकरण विशेषताओं का पहला व्यवस्थित अध्ययन
- उच्च व्यावहारिक मूल्य: Mueller मैट्रिक्स माप आवश्यकताओं को सरल बनाता है, नैदानिक अनुप्रयोग की व्यवहार्यता में सुधार करता है
- पर्याप्त सत्यापन: आदर्श फैंटम से वास्तविक जैविक ऊतक तक क्रमिक सत्यापन रणनीति तर्कसंगत है
- गहन विश्लेषण: Lu-Chipman अपघटन भौतिक तंत्र की गहन समझ प्रदान करता है
- नमूना विविधता सीमित: केवल पॉलीप्रोपाइलीन-जिलेटिन फैंटम का उपयोग किया गया, सामग्री विशेषताएं वास्तविक ऊतक से भिन्न हो सकती हैं
- सांख्यिकीय विश्लेषण अपर्याप्त: विस्तृत सांख्यिकीय महत्व परीक्षण और त्रुटि विश्लेषण की कमी
- सीमित नैदानिक सत्यापन: केवल एक जलने के नमूने पर सत्यापन, बड़े पैमाने पर नैदानिक अनुसंधान की आवश्यकता है
- जटिल बिखराव विचार अपर्याप्त: मुख्य रूप से गोलाकार कणों पर विचार करता है, वास्तविक ऊतक संरचना अधिक जटिल है
- शैक्षणिक मूल्य: टेराहर्ट्ज़ बायोफोटोनिक्स के लिए महत्वपूर्ण सैद्धांतिक आधार और प्रायोगिक विधि प्रदान करता है
- तकनीकी प्रगति: चिकित्सा निदान में टेराहर्ट्ज़ ध्रुवीकरण इमेजिंग तकनीक के अनुप्रयोग को आगे बढ़ाता है
- अंतः-विषय महत्व: विद्युत चुंबकीय बिखराव सिद्धांत, Monte Carlo मॉडलिंग और जैव चिकित्सा इमेजिंग को जोड़ता है
- औद्योगिक संभावनाएं: टेराहर्ट्ज़ चिकित्सा उपकरणों के विकास के लिए तकनीकी पथ प्रदान करता है
- त्वचा रोग निदान: जलने की ग्रेडिंग, त्वचा कैंसर का पता लगाना
- ट्यूमर सीमा पहचान: सर्जिकल नेविगेशन, रोग संबंधी विश्लेषण
- ऊतक संरचना विशेषता: कोलेजन फाइबर, कोशिका समूह विश्लेषण
- दवा चिकित्सा प्रभाव निगरानी: उपचार प्रक्रिया में ऊतक संरचना परिवर्तन ट्रैकिंग
पेपर में 41 संबंधित संदर्भों का हवाला दिया गया है, जिसमें टेराहर्ट्ज़ तकनीक, जैव चिकित्सा प्रकाशिकी, Monte Carlo मॉडलिंग आदि कई क्षेत्रों के महत्वपूर्ण कार्य शामिल हैं, जो अनुसंधान के लिए एक मजबूत सैद्धांतिक आधार प्रदान करते हैं।
समग्र मूल्यांकन: यह टेराहर्ट्ज़ जैव चिकित्सा इमेजिंग क्षेत्र में एक उच्च गुणवत्ता वाला अंतः-विषय अनुसंधान पत्र है, जिसमें महत्वपूर्ण सैद्धांतिक और व्यावहारिक मूल्य है। अनुसंधान विधि कठोर है, परिणाम विश्वसनीय हैं, और इस क्षेत्र के आगे विकास के लिए एक अच्छा आधार प्रदान करते हैं।