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

ملخص تقني: نموذج دائرة مكافئة (RC) مجمعة لأنسجة الرأس في نطاقات التردد دون الميغاهرتز

بيان المشكلة
يعد النمذجة الدقيقة للجهد الكهربائي وتوزيع التيار في أنسجة الرأس أمراً ضرورياً لتصميم وتقييم أنظمة الاستشعار العصبي (مثل تخطيط كهربية الدماغ EEG) وأنظمة التحفيز العصبي (مثل tDCS، وtACS، وDBS) التي تعمل في نطاق الترددات دون الميغاهرتز. وبينما يمكن للطرق العددية مثل طريقة العناصر المحدودة (FEM) وطريقة العناصر الحدودية (BEM) التعامل مع الهندسات المعقدة، فإن تكلفتها الحسابية العالية تحد من فائدتها في عمليات النمذجة السريعة، والمحاكاة في الوقت الفعلي، والتكامل المباشر في بيئات تصميم الدوائر على مستوى الدائرة. تعتمد النهج القائمة على الدوائر، مثل نماذج شبكة المقاومات (RMM)، غالباً على تقسيم حجمي دقيق، مما يؤدي إلى متطلبات حسابية غير مستهان بها، وكثيراً ما تفشل في مراعاة الخصائص السعوية والتشتت المعتمد على التردد للأنسجة البيولوجية. هناك حاجة إلى دائرة مكافئة ذات معاملات مجمعة (lumped-parameter) تحافظ على الدقة عبر طيف ترددي واسع مع دمج تشتت الأنسجة وتأثيرات تيار الإزاحة.

المنهجية
يقترح المؤلفون نموذج دائرة مكافئة (RC) مجمعة مشتق من هندسة رأس كروية نموذجية مكونة من ثلاث طبقات (فروة الرأس، الجمجمة، والدماغ) محاطة بطبقة هواء لانهائية. يعمل النموذج ضمن نظام الكهروسكوني (EQS)، باستخدام موصلية معقدة تعتمد على التردد ($\sigma_c(\omega) = \sigma(\omega) + i\omega\varepsilon(\omega)$) لالتقاط الاستجابات الموصلة والسعوية معاً.

• طوبولوجيا الدائرة: يستخدم النموذج مجموعة دنيا من عناصر المعاوقة. تمثل كل طبقة نسيجية عبر فرعين مقترنين: معاوقة شعاعية ($Z^r_i$) لنمذجة تدفق التيار عبر الطبقة، ومعاوقة مماسية ($Z^\theta_i$) لنمذجة تيارات العودة داخل الطبقة. يراعي هذا الهيكل ثنائي الفروع حقيقة أن بعض خطوط المجال لا تعبر الحجم بالكامل؛ بل يعود بعض التيار للانتشار داخل غلاف قبل العودة.
• الشروط الحدودية: يتم فرض قانون كيرشوف للتيار (KCL) عند العقد البينية لتحقيق استمرارية المكون العمودي لكثافة التيار الكلية، مما يعيد إنتاج تراكم الشحنة السطحية بسبب عدم استمرارية السماحية بشكل طبيعي.
• البارامترات (المعلمات):
  • التباين الهندسي: يتم التحكم في نسبة المعاوقة الشعاعية إلى المماسية عبر معامل هندسي لا بُعدي ($\gamma_i$)، وهو دالة لنصاف أقطار الطبقات وسماكاتها. تم اشتقاق ملاءمات متعددة الحدود للسماح للنموذج بالتكيف مع التباينات في سمك الجمجمة ونسب الدماغ إلى فروة الرأس.
  • لامركزية ثنائي القطب: للتعامل مع المصادر العصبية غير المركزية، يتم تقسيم المعاوقة الشعاعية للدماغ إلى فرعين غير متماثلين ($Z^{r_{up}}_{brain}$ و $Z^{r_{down}}_{brain}$)، ويتم قياسهما بمعامل عدم تماثل $\alpha(\eta)$ يعتمد على الموضع الشعاعي لثنائي القطب.
  • الاعتماد على التردد: يدمج النموذج بيانات الموصلية والسماحية المعتمدة على التردد (بناءً على قياسات in vivo بواسطة Wagner وآخرون) لنمذجة تشتت الأنسجة.
• التحقق من الصحة: تم تنفيذ النموذج في برنامج QSPICE والتحقق من صحته مقابل حل مرجعي شبه تحليلي يعتمد على توسيع التوافقيات الكروية القياسية (SSH). غطت عملية التحقق نطاقاً ترددياً من 10 هرتز إلى 50 كيلوهرتز، مع مختلف درجات لامركزية ثنائي القطب، ومختلف سماكات الجمجمة.

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

1. تمثيل دائرة مدمج: يقدم البحث دائرة ذات عناصر مجمعة تعيد إنتاج السلوك الكهربائي لرأس كروي متعدد الطبقات باستخدام الحد الأدنى من المكونات، مما يسمح بالتكامل في برامج محاكاة الدووت القياسية (SPole).
2. إدراج تأثيرات التشتت والسعة: على عكس النماذج المبسطة السابقة، يتضمن هذا النهج صراحةً الموصلية والسماحية المعتمدة على التردد، بالإضافة إلى تيارات الإزاحة، عبر عناصر المعاوقة المعقدة.
3. القابلية للتكيف الهندسي والمصدري: تم تصميم بارامترات النموذج لتتدرج مع التباينات التشريحية (سمك الجمجمة، أنصاف أقطار الطبقات) ولامركزية المصدر، مما يجعله قوياً لمختلف التكوينات الفسيولوجية.
4. التقييم الكمي للتبسيطات: تقيم الدراسة كمياً الأخطاء الناتجة عن إهمال تيارات الإزاحة وتشتت التردد، مما يثبت دورها الحاسم في دقة النموذج.

النتائج

• الدقة: أظهر نموذج الدائرة المجمعة المقترح توافقاً وثيقاً مع الحل المرجعي شبه التحليلي (SSH) عبر نطاق 10 هرتز–50 كيلوهرتز. ظل متوسط خطأ التردد النسبي (MRFE) منخفضاً، مع ملاحظة انحرافات طفيفة خاصة عند حالات لامركزية ثنائي القطب العالية حيث يقترب المصدر من الواجهة بين الدماغ والجمجمة.
• تأثير التشتت وتيارات الإزاحة: كشف التحليل العددي أن إهمال هذه التأثيرات يؤدي إلى أخطاء كبيرة.
  • أدى النموذج الأومي الصرف (مقاومة ثابتة، لا سعة) إلى وصول الأخطاء النسبية إلى حوالي 160% عند 50 كيلوهرتز.
  • أدى تضمين الموصلية المشتتة للتردد وحده إلى تقليل الأخطاء بشكل كبير، لكنه ظل ينتج تباينات تتجاوز 10% من 600 هرتز وتصل إلى 15% عند 1.7 كيلوهرتز.
  • فقط من خلال تضمين كل من المقاومة المشتتة للتردد والسعة (تيارات الإزاحة) تم ضمان تنبؤات دقيقة ومتسقة فيزيائياً.
• الحساسية الهندسية: وُجد أن الخطأ في تقدير جهد فروة الرأس أكثر حساسية للامركزية ثنائية القطب منه لتباينات سمك الجمجمة، رغم أن كلا العاملين يؤثران على أداء النموذج.

الأهمية والادعاءات
يزعم المؤلفون أن هذا العمل يوفر بديلاً عملياً وفعالاً حسابياً للمحللات العددية الحجمية الكاملة لتطبيقات محددة في الاستشعار والتحفيز العصبي. ومن خلال تقديم تمثيل كهربائي مدمج متوافق مع بيئات محاكاة الدوائر، يسهل النموذج عمليات التصميم السريع، وتقييم نسبة الإشارة إلى الضوضاء، والضبط التكيفي في الوقت الفعلي لمعلمات التحفيز. وتؤكد الدراسة أن النمذجة الدقيقة في نطاق ما دون الميغاهرتز تتطلب إدراج كل من التأثيرات السعوية وتشتت الأنسجة لتجنب المبالغة الكبيرة في تقدير جهود فروة الرأس. وبينما اقتصر التحقق الحالي على ثنائيات القطب الموجهة شعاعياً في طبقات كروية متحدة المركز ومتناحية، فإن هذا الإطار يضع أساساً للتوسعات المستقبلية نحو هندسات أكثر تعقيداً وخصائص أنسجة متباينة المناحي.