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

技術要約：時間領域コークレア伝搬線モデルにおける短波長アドミッタンス補正

問題提起
時間領域の伝搬線（TL）モデルは、過渡的な音に対する基底膜（BM）変位をシミュレートするための効率的なツールである。しかし、その1次元（1-D）的な性質により、短波長領域のコークレアにおいて顕著となる高次元の流体力学的効果（具体的には圧力集束および横方向の粘性減衰）を統合することには限界がある。標準的なV-1Dモデル（Verulst et al., 2018）では、利得（ゲイン）と周波数選択性が密接に結合している。すなわち、高い利得を得るには鋭いチューニングが必要であり、一方でチューニングを広げると利得が低下するという関係にある。この結合は、ジャービル（オコジョモグラット）のような小型哺乳類にモデルを適応させる際に大きな課題となる。実験データによれば、ジャービルのコークレア・フィルターは、低音圧レベルにおいて比較的広いチューニングを持ちつつも高いピークを示すことが示されている。一方で、広範な圧縮非線形性も維持している。標準的なV-1Dモデルはこれを再現できず、生理学的なリアリズムを損なうことなく利得と周波数選択性を分離できないため、不十分な圧縮（約10 dB不足）を示す。

手法
著者らは、既存の時間領域V-1Dフレームワークに短波長アドミッタンスの数値補正を統合した修正モデルである$V^*$モデルを提案している。その手法は、主に以下の3つのフェーズで構成される。

1. 参照モデリング（S-2D）： 著者らは、2次元の流体効果（圧力集束および粘性減衰）を組み込んだ周波数領域のTLモデルであるS-2Dモデル（Sisto et al., 2021）を利用している。このモデルは、圧力集束に起因する利得増幅係数（$\alpha$）の「ターゲット」を提供する。基準となる能動過程の強さ（$G_{REF}$）に対して、正規化された補正係数$\beta = \alpha / \alpha_0$が定義される。
2. 回帰およびフィルター設計： $\beta$は周波数、基底膜特性周波数（$\omega_{BM}$）、および能動過程の強さ（$G$）によって変化し、単純な因果的（causal）な時間領域項として直接実装することはできないため、著者らは自己回帰（オールポール）フィルターを用いて$\beta$を近似している。
   • ステップ1： 固定された$(\omega_{BM}, G)$のペアに対して、ターゲットとなる$\beta$の周波数応答に一致するように、ペナルティ付き最小二乗回帰法を用いて32個の自己回帰係数を決定する。ペナルティ項は、係数の安定性を確保し、高周波での不安定性を防止する。
   • ステップ2： 連続的な入力$(\omega_{BM}, G)$から32個のフィルター係数のセットへのマッピングを行うために、放射基底関数ニューラルネットワーク（RBF-NN）を学習させ、シミュレーション中の動的な更新を可能にする。
3. 時間領域への統合： $V^*$モデルは、これらの再帰的フィルターをBMアドミッタンス方程式に組み込んでいる。フィルター係数は、出力強度（intensity）の短期平均に基づいて動的に更新され、これが能動過程の強さ（$G$）を変調する。計算コストと精度のバランスをとるため、著者らは事前計算された係数のルックアップテーブルを使用し、更新レートを0.03 ms（シミュレーションのタイムステップより6倍遅い設定）とした。これは十分なトレードオフであると判断された。

主な結果

• 利得と選択性のデカップリング： $V^*$モデルは、周波数選択性と利得の間のある程度のデカップリング（分離）に成功した。V-1Dモデルでは、ジャービルで見られるような広いチューニングと高い圧縮を同時に達成できなかったが、$V^*$モデルは圧縮領域を10 dB拡張した。
• 利得増幅： 本モデルは、補正なしのV-1Dモデルと比較して、波のピークにおいて5 dBの追加利得を達成した。この改善により、低音圧レベルにおいて非現実的に鋭いチューニングを必要とすることなく、ジャービルのBM速度の成長関数（He et al., 2022）にモデルをより良く適合させることが可能となった。
• 安定性と精度： 自己回帰フィルターは安定しており、z領域における不安定性は観察されなかった。回帰法は圧力集束因子の一般的な傾向を捉えたが、利得の増幅を過小評価し、利得の減少を過大評価する傾向があった。これは、本手法が中程度の補正に適していることを示唆している。
• 計算効率： フィルター係数の動的更新により、当初は計算時間が20倍増加したが、事前計算されたルックアップテーブルを使用することで、シミュレーション時間を元のV-1Dモデルと同等のレベルまで短縮できた。

意義と主張
本論文は、$V^*$モデルが、時間領域TLモデルにおいて波長依存の効果を実装するための実用的な解決策を提供すると主張している。これは特に、小型哺乳類のモデルにおける圧縮不足の問題に対処するものである。本研究の主な貢献は以下の2点である。

1. 手法の統合： 解析的な流体力学（S-2D）と、回帰ベースのフィルター設計（RBF-NNおよび自己回帰フィルター）を組み合わせることで、周波数領域の2次元効果を時間領域の1次元フレームワークへと変換するハイブリッドなアプローチを実証したこと。
2. 非線形性の組み合わせ： 本モデルは、瞬時的な非線形性（瞬時のBM速度に依存）と、非瞬時的な非線形性（利得補正のフィードバックループを介した信号エンベロープに依存）を独自に組み合わせている。これにより、圧力集束が能動過程によって変調される受動現象として機能する、複雑なハイブリッド力学をシミュレートすることが可能となる。

著者らは、控えめながらも、$V^*$モデルはジャービルの基底回転における実験データとの一致を改善するものの、「チップ・トゥ・テイル（先端から尾部にかけて）」の比率の不一致（絶対的な利得が依然として実験的な観察値よりも低いこと）を完全には解決していないと述べている。これは、高利得かつ因果的なフィルターを時間領域モデルに実装することの根本的な困難さに起因すると考えている。そこでは、局所的な因果性の制約が、2次元の物理モデルに見られるグローバルな波の因果性や瞬時的な圧力集束を再現することを制限するためである。結論として、本アプローチは、豊富な生理学的データが存在しながらも標準的な1次元の仮定ではモデリングが困難な、小型実験動物を用いたトランスレーショナル研究において特に価値がある。