Short-wave admittance correction for a time-domain cochlear transmission line model

ملخص تقني: تصحيح الممانعة للموجات القصيرة لنموذج خط النقل القوقعي في المجال الزمني

بيان المشكلة
تعد نماذج خط النقل (TL) في المجال الزمني أدوات فعالة لمحاكاة إزاحة الغشاء القاعدي (BM) استجابةً للأصوات العابرة. ومع ذلك، فإن طبيعتها أحادية البعد (1-D) تحد من دمج تأثيرات الهيدروديناميكا متعددة الأبعاد، وتحديداً تركيز الضغط والتخميد اللزوجي المستعرض، وهي تأثيرات بارزة في منطقة الموجات القصيرة في القوقعة. في نموذج V-1D القياسي (Verhulst et al., 2018)، يرتبط الكسب والاصطفاء الترددي ارتباطاً وثيقاً: فالحصول على كسب عالٍ يتطلب اصطفاءً حاداً، بينما يؤدي الاصطفاء الأوسع إلى كسب أقل. يشكل هذا الارتباط تحدياً كبيراً عند تكييف النموذج مع الثدييات الصغيرة مثل الجربوع (gerbil). تشير البيانات التجريبية إلى أن مرشحات قوقعة الجربوع تظهر قممًا عالية مع اصطفاء واسع نسبياً عند مستويات الصوت المنخفضة، ومع ذلك تحافظ على نطاق واسع من اللاخطية الانضغاطية. يفشل نموذج V-1D القياسي في إعادة إنتاج ذلك، حيث يظهر انضغاطاً غير كافٍ (أقل بنحو 10 ديسيبل) لأن صياغته أحادية البعد لا يمكنها فك الارتباط بين الكسب والاصطفاء الترددي دون المساس بالواقعية الفسيولوجية.

المنهجية
يقترح المؤلفون نموذجاً معدلاً، يُسمى نموذج $V^*$، والذي يدمج تصحيحاً عددياً للممانعة للموجات القصيرة ضمن إطار عمل V-1D الحالي في المجال الزمني. تتضمن المنهجية ثلاث مراحل رئيسية:

1. النمذجة المرجعية (S-2D): يستخدم المؤلفون نموذج S-2D (Sisto et al., 2021)، وهو نموذج خط نقل في المجال الترددي يتضمن تأثيرات السوائل ثنائية الأبعاد (2-D) (تركيز الضغط والتخميد اللزوجي). يوفر هذا النموذج "هدفاً" لعامل تعزيز الكسب ($\alpha$) الناتج عن تركيز الضغط. تم تعريف عامل تصحيح معياري، $\beta = \alpha / \alpha_0$، بالنسبة لمرجع قوة العملية النشطة ($G_{REF}$).
2. الانحدار وتصميم المرشح: بما أن $\beta$ يتغير مع التردد، والتردد المميز للغشاء القاعدي ($\omega_{BM}$)، وقوة العملية النشطة ($G$)، ولا يمكن تنفيذه مباشرة كحد بسيط سببي في المجال الزمني، فقد قام المؤلفون بتقريب $\beta$ باستخدام مرشحات ذاتية الانحدار (all-pole filters).
   • الخطوة 1: بالنسبة لأزواج ثابتة من $(\omega_{BM}, G)$، تم تحديد 32 معامل انحدار ذاتي عبر انحدار المربعات المعاقب لمطابقة الاستجابة الترددية لـ $\beta$ المستهدف. تضمن حدود العقوبة استقرار المعاملات وتمنع عدم الاستقرار في الترددات العالية.
   • الخطوة 2: تم تدريب شبكة عصبية ذات وظيفة أساس شعاعية (RBF-NN) لربط المدخلات المستمرة $(\omega_{BM}, G)$ بمجموعة من معاملات المرشح الـ 32، مما يسمح بتحديثات ديناميكية أثناء المحاكاة.
3. التكامل في المجال الزمني: يدمج نموذج $V^*$ هذه المرشحات التكرارية في معادلات ممانعة الغشاء القاعدي (BM). يتم تحديث معاملات المرشح ديناميكياً بناءً على متوسط قصير المدى لشدة المخرجات، والذي يعدل قوة العملية النشطة ($G$). ولتحقيق التوازن بين التكلفة الحسابية والدقة، استخدم المؤلفون جدول بحث (lookup table) للمعاملات المحسوبة مسبقاً، واختاروا معدل تحديث قدره 0.03 مللي ثانية (أبطأ بـ 6 مرات من خطوة زمن المحاكاة)، وهو ما وُجد أنه يمثل حلاً وسطاً كافياً.

النتائج الرئيسية

• فك الارتباط بين الكسب والاصطفاء: نجح نموذج $V^*$ في إدخال درجة من فك الارتباط بين الاصطفاء الترددي والكسب. وبينما لم يستطع نموذج V-1D تحقيق كل من الاصطفاء الواسع والكسب العالي الملاحظين في الجرابيع في آن واحد، استطاع نموذج $V^*$ توسيع نطاق النظام الانضغاطي بمقدار 10 ديسيبل.
• تعزيز الكسب: حقق النموذج كسباً إضافياً قدره 5 ديسيبل عند ذروة الموجة مقارنة بنموذج V-1D غير المصحح. سمح هذا التحسين للنموذج بمطابقة دوال النمو التجريبية لسرعة الغشاء القاعدي في الجرابيع (He et al., 2022) بشكل أفضل دون الحاجة إلى اصطفاء حاد بشكل غير واقعي عند مستويات الصوت المنخفضة.
• الاستقرار والدقة: ظلت المرشحات ذاتية الانحدار مستقرة، مع عدم ملاحظة أي عدم استقرار في نطاق z. نجحت طريقة الانحدار في التقاط الاتجاه العام لعامل تركيز الضغط، وإن كانت قد قللت قليلاً من تعزيز الكسب وبالغت في تقدير خفض الكسب، مما يشير إلى أن الطريقة هي الأنسب للتصحيحات المتوسطة.
• الكفاءة الحسابية: على الرغم من أن التحديث الديناميكي لمعاملات المرشح زاد وقت الحساب في البداية بمعامل 20، إلا أن استخدام جدول بحث محوسب قلل أوقات المحاكة لتصبح مقاربة لنموذج V-1D الأصلي.

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

1. التكامل المنهجي: يوضح نهجاً هجيناً يجمع بين ديناميكا السوائل التحليلية (S-2D) وتصميم المرشحات القائم على الانحدار (RBF-NN والمرشحات ذاتية الانحدار) لترجمة تأثيرات المجال الترددي ثنائية الأبعاد (2-D) إلى إطار عمل أحادي البعد (1-D) في المجال الزمني.
2. الجمع بين اللاخطية: يجمع النموذج بشكل فريد بين اللاخطية اللحظية (المعتمدة على سرعة الغشاء القاعدي اللحظية) واللاخطية غير اللحظية (المعتمدة على غلاف الإشارة عبر حلقة التغذية الراجعة لتصحيح الكسب). وهذا يسمح للنموذج بمحاكاة ديناميكيات هجينة معقدة حيث يعمل تركيز الضغط كظاهرة سلبية يتم تعديلها بواسطة العملية النشطة.

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