論文要約：異方性双ドメイン演算子の固有値問題に対する境界積分アプローチ

1. 問題の背景と目的

本論文は、心臓の電気生理現象を記述する双ドメインモデル（Bidomain Model）において、左心室、中隔（セプタム）、右心室の 3 つの領域からなる異方性不均質領域における固有値問題を研究対象としています。

従来の双ドメインモデルの数値解析は、有限要素法などの体積離散化に依存しており、複雑な心臓幾何学と組織の異方性を正確に捉えるために計算コストが極めて高いという課題がありました。また、中隔が機能性バリアとして振る舞い、心室間での波動伝播が阻害される現象の数学的メカニズム（特に伝導率テンソルの不連続性）を、境界条件として明確に定式化した研究は不足していました。

本研究の目的は、ポテンシャル理論を用いてこの境界値問題を境界積分方程式系に変換し、効率的な固有値計算手法を構築することにあります。

2. 数学的定式化とモデル設定

• 領域設定: 2 次元空間 $\Omega$ を、左心室 ($\Omega_I$)、中隔 ($\Omega_S$)、右心室 ($\Omega_D$) の 3 つの矩形部分領域に分割します。
• 支配方程式: 細胞内・細胞外電位 $u_i, u_e$ に対する異方性ヘルムホルツ方程式を解きます。
  $$\nabla \cdot (\sigma(x) \nabla u) + \lambda^2 u = 0$$
  ここで、$\sigma(x)$ は異方性伝導率テンソル、$\lambda$ は固有値です。
• 異方性構造: 各領域において、繊維方向（長手方向）とそれに直交する方向（横方向）の伝導率 $\sigma_l, \sigma_t$ が定義されます。特に、中隔における繊維方向が心室とは異なる方向を持つことを考慮し、テンソルが領域境界で不連続になるように設定しています。
• 境界・接触条件:
  • 外部境界では、ノイマン型境界条件（法線方向の電流フラックスがゼロ）を課します。
  • 領域間境界（$\Gamma_{IS}, \Gamma_{SD}$）では、**完全接触条件（Perfect Contact Conditions）**を課します。すなわち、電位と電流フラックスの両方が連続であることを要求します。これは、伝導率テンソルが不連続であっても電流が連続に流れるという生理学的仮定に基づいています。

3. 手法：ポテンシャル理論と境界積分方程式への帰着

本研究の核心的な手法は、ポテンシャル理論を用いて偏微分方程式を境界積分方程式に変換することです。

1. 基本解の導出:
   各領域における異方性ヘルムホルツ演算子に対する基本解を明示的に導出しました。変数変換により等方性のヘルムホルツ方程式に変換し、その解を第 2 種ベッセル関数 $N_0$（対数特異性を持つ）を用いて表現します。
2. 直接定式化（Direct Formulation）:
   解を単層ポテンシャルと双層ポテンシャルの線形結合として表現します。
   $$u(x) = \int_{\partial \Omega} u(y) p^*(x,y) ds(y) - \int_{\partial \Omega} p(y) u^*(x,y) ds(y)$$
   ここで、$u$ は電位、$p$ は電流フラックス、$u^*, p^*$ はそれぞれ基本解とその法線微分です。
3. 境界積分方程式系（BIEs）の構築:
   境界上の極限値に関する跳躍関係（Proposition 9）を用いることで、元の偏微分方程式問題を、境界上の未知関数（電位と電流フラックス）に関するフレドホルム型の第二種境界積分方程式系に帰着させました。
   • 本問題では、3 つの領域と 2 つの界面、外部境界を考慮し、合計7 つの方程式からなる連立方程式系を導出しました。
   • 境界の角点における特異性を補正するための係数 $\kappa(z)$ を導入し、積分の収束性を確保しています。

4. 数値解法と結果

• 離散化: 境界と界面を均等に分割し、特異積分（弱特異性を持つ対数特異点を含む）を扱うための特殊な数値積分公式（クアドラチュア）を適用して、線形代数方程式系に変換しました。
• 核関数の評価: 核関数はベッセル関数 $N_0, N_1$ で表され、級数展開の最初の数項を用いて近似計算を行いました。
• 固有値の決定: 得られた線形方程式系が自明でない解を持つようなパラメータ $\lambda$ を探索することで、双ドメイン演算子の固有値を近似しました。行列のランク欠損（最小特異値のゼロへの収束、または行列式の零点）を検出することで固有値を同定しました。
• 数値結果:
  • 計算された電位分布は、左心室、中隔、右心室の各領域で連続的かつ整合的な挙動を示しました。
  • 完全接触条件を課しているにもかかわらず、伝導率の不連続性により電位が界面で急激に変化する（ジャンプする）挙動が確認されました。これは、完全な双ドメイン問題（拡散項を含む場合）において、波動伝播速度の変化を引き起こす可能性を示唆しています。
  • 計算された電位の絶対値は非常に小さく（$10^{-46}$ 程度）、外部刺激がない受動状態におけるモデルの特性を反映しています。

5. 主な貢献と意義

1. 次元削減と計算効率: 従来の体積離散化（FEM など）に代わり、境界積分法を採用することで、問題の次元を 1 つ削減しました。これにより、複雑な心臓幾何学における異方性問題の解析効率が向上します。
2. 理論的枠組みの確立: 異方性ヘルムホルツ演算子に対する基本解をベッセル関数で明示的に導出し、それを基にした境界積分定式化を完成させました。
3. 中隔の機能性バリアの数学的記述: 伝導率テンソルの不連続性を境界条件として取り込むことで、中隔が心室間の電気的伝播に与える影響を数学的にモデル化し、そのメカニズムを解明する新たなアプローチを提供しました。
4. 数値アルゴリズムの提案: 特異積分を含む境界積分方程式系を効率的に離散化し、固有値を計算する具体的な数値スキームを提案しました。

6. 結論

本論文は、心臓の異方性双ドメインモデルにおける固有値問題を、ポテンシャル理論に基づいた境界積分方程式法によって効率的に解くための包括的な枠組みを提示しました。このアプローチは、心臓電気生理学的な複雑な現象（特に界面での波動伝播の挙動）を、従来の体積法よりも低コストかつ高精度に解析する可能性を開くものであり、将来の心臓シミュレーションや病態解析への応用が期待されます。