A Lumped RC Equivalent Circuit Model of Head Tissues in sub-MHz Frequency Regimes

तकनीकी सारांश: सब-MHz आवृत्ति शासन (regime) में सिर के ऊतकों का एक लम्पड (Lumped) RC समतुल्य परिपथ मॉडल

समस्या विवरण
सब-मेगाहर्ट्ज़ (MHz) आवृत्ति सीमा में संचालित होने वाले न्यूरो-सेंसिंग (जैसे, EEG) और न्यूरो-स्टिमुलेशन (जैसे, tDCS, tACS, DBS) प्रणालियों के डिजाइन और मूल्यांकन के लिए सिर के ऊतकों के विद्युत विभव (electric potential) और धारा वितरण का सटीक मॉडलिंग आवश्यक है। जबकि परिमित तत्व विधि (FEM) और सीमा तत्व विधि (BEM) जैसी संख्यात्मक विधियाँ जटिल ज्यामिति को संभाल सकती हैं, उनकी उच्च कम्प्यूटेशनल लागत तीव्र प्रोटोटाइपिंग, वास्तविक समय के सिमुलेशन और सर्किट-स्तरीय डिजाइन वातावरण में सीधे एकीकरण की उपयोगिता को सीमित करती है। मौजूदा परिपथ-आधारित दृष्टिकोण, जैसे कि रेजिस्टर मेश मॉडल (RMM), अक्सर सूक्ष्म वॉल्यूमेट्रिक विविक्तिकरण (discretization) पर निर्भर करते हैं, जिसके परिणामस्वरूप गैर-नगण्य कम्प्यूटेशनल मांगें होती हैं, और वे अक्सर जैविक ऊतकों के कैपेसिटिव गुणों और आवृत्ति-निर्भर फैलाव (dispersion) को समझाने में विफल रहते हैं। एक कॉम्पैक्ट, लम्पड-पैरामीटर समतुल्य परिपथ की आवश्यकता है जो ऊतक फैलाव और विस्थापन धारा (displacement current) प्रभावों को शामिल करते हुए एक विस्तृत आवृत्ति स्पेक्ट्रम में सटीकता बनाए रखे।

कार्यप्रणाली
लेखक एक लम्पड RC समतुल्य परिपथ मॉडल प्रस्तावित करते हैं जिसे एक मानक तीन-परत गोलाकार सिर की ज्यामिति (स्कैल्प, खोपड़ी और मस्तिष्क) से प्राप्त किया गया है, जो एक अनंत वायु परत से घिरा हुआ है। यह मॉडल इलेक्ट्रो-क्वासी-स्टैटिक (EQS) शासन के भीतर कार्य करता है, जो कंडक्टिव और कैपेसिटिव दोनों प्रतिक्रियाओं को पकड़ने के लिए जटिल आवृत्ति-निर्भर चालकता ($\sigma_c(\omega) = \sigma(\omega) + i\omega\varepsilon(\omega)$) का उपयोग करता है।

• परिपथ टोपोलॉजी (Circuit Topology): मॉडल न्यूनतम प्रतिबाधा तत्वों (impedance elements) के सेट का उपयोग करता है। प्रत्येक ऊतक परत को दो युग्मित शाखाओं द्वारा दर्शाया गया है: एक रेडियल प्रतिबाधा ($Z^r_i$) जो परत के माध्यम से धारा प्रवाह को मॉडल करती है और एक स्पर्शरेखीय (tangential) प्रतिबाधा ($Z^\theta_i$) जो इंट्रा-लेयर रिटर्न करंट को मॉडल करती है। यह द्वि-शाखा संरचना इस तथ्य को ध्यान में रखती है कि सभी क्षेत्र रेखाएं (field lines) पूरे आयतन को पार नहीं करती हैं; कुछ धारा एक शेल के भीतर प्रसारित होकर वापस लौटती है।
• सीमा स्थितियाँ (Boundary Conditions): कुल धारा घनत्व के सामान्य घटक की निरंतरता को संतुष्ट करने के लिए इंटरफेशियल नोड्स पर किर्चॉफ के करंट लॉ (KCL) को लागू किया गया है, जो स्वाभाविक रूप से सतह आवेश संचय (surface charge accumulation) को पुनरुत्पादित करता है।
• पैरामीटराइजेशन (Parameterization):
  • ज्यामितीय परिवर्तनशीलता: रेडियल और स्पर्शरेखीय प्रतिबाधाओं के बीच का अनुपात एक आयामहीन ज्यामितीय गुणांक ($\gamma_i$) द्वारा नियंत्रित होता है, जो परतों की त्रिज्या और मोटाई का एक फलन है। बहुपद फिट (Polynomial fits) को यह अनुमति देने के लिए व्युत्पन्न किया गया था कि मॉडल खोपड़ी की मोटाई और मस्तिष्क-स्कैल्प अनुपात में विविधताओं के अनुकूल हो सके।
  • डिपोल सनक (Dipole Eccentricity): ऑफ-सेंटर न्यूरल स्रोतों को संभालने के लिए, मस्तिष्क की रेडियल प्रतिबाधा को दो असममित शाखाओं ($Z^{r_{up}}_{brain}$ और $Z^{r_{down}}_{brain}$) में विभाजित किया गया है, जिन्हें डिपोल की रेडियल स्थिति पर निर्भर असममित पैरामीटर $\alpha(\eta)$ द्वारा स्केल किया गया है।
  • आवृत्ति निर्भरता: मॉडल ऊतक फैलाव को मॉडल करने के लिए आवृत्ति-निर्भर चालकता और पारद्युतता (permittivity) डेटा (Wagner et al. के in vivo मापों पर आधारित) को शामिल करता है।
• सत्यापन (Validation): मॉडल को QSPICE में कार्यान्वित किया गया और स्केलर गोलाकार हार्मोनिक्स (SSH) विस्तार पर आधारित एक अर्ध-विश्लेलेतिक संदर्भ समाधान के विरुद्ध सत्यापित किया गया। सत्यापन ने 10 Hz से 50 kHz तक की आवृत्ति सीमा, विभिन्न डिपोल सनक और विभिन्न खोपड़ी की मोटाई को कवर किया।

प्रमुख योगदान

1. कॉम्पैक्ट सर्किट प्रतिनिधित्व: यह पेपर एक न्यूनतम घटकों वाले लम्पड-एलिमेंट सर्किट को पेश करता है जो एक बहु-परत गोलाकार सिर के विद्युत व्यवहार को पुनरुत्पादित करता है, जिससे मानक सर्किट सिम्युलेटर (SPICE) में एकीकरण संभव होता है।
2. फैलाव और कैपेसिटिव प्रभावों का समावेश: पिछले सरलीकृत मॉडलों के विपरीत, यह दृष्टिकोण जटिल प्रतिबाधा तत्वों के माध्यम से आवृत्ति-निर्भर चालकता और पारद्युतता, साथ ही विस्थापन धाराओं को स्पष्ट रूप से शामिल करता है।
3. ज्यामितीय और स्रोत अनुकूलन क्षमता: मॉडल के पैरामीटर शारीरिक विविधताओं (खोपड़ी की मोटता, परत की त्रिज्या) और स्रोत सनक के साथ स्केल करने के लिए डिज़ाइन किए गए हैं, जो उन्हें विभिन्न शारीरिक विन्यासों के लिए मजबूत बनाते हैं।
4. सरलीकरण का मात्रात्मक मूल्यांकन: अध्ययन विस्थापन धाराओं और आवृत्ति फैलाव को अनदेखा करने से होने वाली त्रुटियों का मात्रात्मक मूल्यांकन करता है, जो मॉडल की निष्ठा (fidelity) के लिए उनकी महत्वपूर्ण भूमिका को प्रदर्शित करता है।

परिणाम

• सटीकता: प्रस्तावित लम्पड सर्किट मॉडल ने 10 Hz–50 kHz की सीमा में अर्ध-विश्लेलेतिक SSH संदर्भ समाधान के साथ घनिष्ठ सहमति दिखाई। मीन रिलेटिव फ्रीक्वेंसी एरर (MRFE) कम रहा, जिसमें मुख्य रूप से उच्च डिपोल सनक पर मामूली विचलन देखे गए जहाँ स्रोत मस्तिष्क-खोपड़ी इंटरफेस के करीब पहुँच जाता है।
• फैलाव और विस्थापन धाराओं का प्रभाव: संख्यात्मक विश्लेषण से पता चला कि इन प्रभावों को अनदेखा करने से महत्वपूर्ण त्रुटियां होती हैं।
  • एक शुद्ध ओमिक मॉडल (स्थिर प्रतिरोध, कोई कैपेसिटेंस नहीं) के परिणामस्वरूप 50 kHz पर लगभग 160% की सापेक्ष त्रुटियां हुईं।
  • केवल आवृत्ति-फैलाव वाली चालकता को शामिल करने से त्रुटियां काफी कम हो गईं, लेकिन फिर भी 600 Hz से 10% से अधिक और 1.7 kHz पर 15% तक का अंतर बना रहा।
  • केवल आवृत्ति-फैलाव वाले प्रतिरोध और कैपेसिटेंस (विस्थापन धारा) दोनों को शामिल करने से भौतिक रूप से सुसंगत और सटीक भविष्यवाणियां सुनिश्चित हुईं।
• ज्यामितीय संवेदनशीलता: स्कैल्प विभव (potential) के अनुमान में त्रुटि, खोपड़ी की मोटाई में भिन्नता की तुलना में डिपोल सनक के प्रति अधिक संवेदनशील पाई गई, हालांकि दोनों कारकों ने मॉडल के प्रदर्शन को प्रभावित किया।

महत्व और दावे
लेखक दावा करते हैं कि यह कार्य विशिष्ट न्यूरल सेंसिंग और स्टिमुलेशन अनुप्रयोगों के लिए एक व्यावहारिक, कम्प्यूटेशनल रूप से कुशल विकल्प प्रदान करता है। एक कॉम्पैक्ट इलेक्ट्रिकल प्रतिनिधित्व प्रदान करके जो सर्किट सिमुलेशन वातावरण के अनुकूल है, यह मॉडल तीव्र डिजाइन, सिग्नल-टू-नॉइज़ रेशियो मूल्यांकन और स्टिमुलेशन मापदंडों के वास्तविक समय के अनुकूलन (adaptive tuning) को सुगम बनाता है। अध्ययन इस बात पर जोर देता है कि सब-MHz शासन में सटीक मॉडलिंग के लिए स्कैल्प विभव के अत्यधिक आकलन से बचने के लिए कैपेसिटिव प्रभावों और ऊतक फैलाव दोनों को शामिल करना आवश्यक है। हालांकि वर्तमान सत्यापन रेडियल रूप से उन्मुख डिपोल और संकेंद्रित, आइसोट्रोपिक गोलाकार परतों तक सीमित है, यह ढांचा अधिक जटिल ज्यामिति और अनिसोट्रोपिक ऊतक गुणों के लिए भविष्य के विस्तार के लिए एक आधार स्थापित करता है।