論文要約：3 次元線形化カルデロン問題における非有擾乱の再構成

1. 問題設定 (Problem)

本論文は、3 次元単位球 $B \subset \mathbb{R}^3$ 内における線形化カルデロン問題（Linearised Calderón Problem）の解法を扱っています。

• 背景: 導電率 $\gamma$ の分布を境界電流密度から内部を特定する非侵襲的画像化技術（電気インピーダンストモグラフィなど）の基礎となる逆問題です。
• 線形化: 導電率が基準値 $\gamma \equiv 1$ から摂動 $\eta$ だけずれている場合（$\gamma = 1 + \eta$）、非線形な写像 $\Lambda(\gamma)$ を線形近似します。この線形化写像 $F = D\Lambda(1; \eta)$ は、摂動 $\eta$ と境界データ（Neumann-to-Dirichlet 写像の変化量）を結びつけます。
• 従来の課題: 既存の研究（特に 2 次元）では、摂動 $\eta$ が $L^2$ 空間に属することを前提としていました。しかし、より一般的な物理的状況や特異な摂動を扱うためには、$\eta$ が有界でない（unbounded）場合、すなわち $\eta \in L^3(B)$ であっても再構成可能であることが望まれます。
• 本研究の目的: 3 次元空間において、$L^3(B)$ に属する任意の摂動 $\eta$ を、線形化データ $F\eta$ から厳密に直接再構成するアルゴリズムを構築することです。

2. 手法 (Methodology)

本研究の核心は、3 次元 Zernike 基底（3D Zernike basis）を用いた展開と、それに基づく前方代入法（forward substitution）の採用にあります。

2.1 3 次元 Zernike 基底の導入

摂動 $\eta$ を以下の 3 次元 Zernike 基底関数 $\psi_{\ell}^{k,m}$ の級数として展開します。
$$\eta = \sum_{k, \ell} \sum_{m} c_{\ell}^{k,m} \psi_{\ell}^{k,m}$$
ここで、基底関数は以下の通り定義されます。

• $\psi_{\ell}^{k,m}(r, \theta, \phi) = R_{\ell}^k(r) Y_{\ell}^m(\theta, \phi)$
• $Y_{\ell}^m$: 球面調和関数（Spherical Harmonics）
• $R_{\ell}^k$: 3 次元径向 Zernike 多項式
• $k, \ell \in \mathbb{N}_0, m \in \{-\ell, \dots, \ell\}$

この基底は $L^2(B)$ における正規直交基底を形成し、$L^3(B)$ の関数もこの空間に含まれるため展開可能です。

2.2 前方代入による係数再構成

線形化データ $F\eta$ と特定の境界測定の組み合わせを用いることで、展開係数 $c_{\ell}^{k,m}$ を再帰的（inductive）に計算する公式を導出しました。

• 測定の選択: 境界電流密度として、球面調和関数 $Y_{k+1}^0$ のみを使用します（$m$ 指標を変化させません）。これは 2 次元法と比較して必要なデータ量が大幅に少ないことを意味します。
• 再構成のロジック:
  1. $k=0$ の係数（調和関数成分）は、データ $(F\eta)Y_1^0$ から直接計算されます。
  2. $k=1$ の係数は、データ $(F\eta)Y_2^0$ と既知の $k=0$ の係数から計算されます。
  3. 一般に、$k$ 次の係数は、データ $(F\eta)Y_{k+1}^0$ と、それ以前の $k' < k$ の係数を用いて計算されます。
• 数式: 定理 1.1 に示されるように、係数 $c_{\ell}^{k,m}$ は以下の形（前方代入）で得られます。
  $$c_{\ell}^{k,m} = (Q_{\ell,0}^{k,m,k})^{-1} \left( \langle (F\eta)Y_{k+1}^0, Y_{\ell+k+1}^m \rangle - \sum_{q=0}^{k-1} \sum_{s=0}^{k-q} Q_{\ell,s}^{k,m,q} c_{\ell+2s}^{q,m} \right)$$
  ここで $Q$ は $\eta$ に依存しない定数（Gaunt 係数やポッハハマー記号を含む）です。

2.3 数学的拡張

付録 A では、摂動空間を $L^\infty$ から $L^d$（$d$ は空間次元）へ拡張する理論的根拠を示しています。Neumann 境界条件が $L^2$ にあるという仮定が、この拡張を可能にする鍵となります。

3. 主要な貢献と結果 (Key Contributions & Results)

1. 非有界摂動の再構成:
   3 次元において、$L^3(B)$ に属する（有界でない可能性のある）摂動 $\eta$ を、線形化データから厳密に再構成する最初のアルゴリズムを提供しました。

2. 効率的な前方代入アルゴリズム:
   2 次元の手法 [Garde & Hyvönen, 2024] を 3 次元に一般化し、三角構造（triangular structure）を持つ連立方程式を解くことで、反復計算を不要とした直接再構成法を実現しました。これにより数値計算が非常に効率的になります。

3. データ量の削減:
   2 次元ではすべてのフーリエ基底関数を電流密度として使用していましたが、3 次元法では球面調和関数の $m$ 指標（方位角成分）を固定した $Y_{k+1}^0$ のみを使用すれば十分であることが示されました。これは完全な $L^2$ 基底に比べて境界測定データを大幅に削減できることを意味します。

4. 数値的検証:

   • 高精度シミュレーション: 数値計算ツール（Mathematica）を用いた高精度なデータで、切断次数を制限してもほぼ完全な再構成が可能であることを確認しました。
   • 現実的なノイズ: 粗い有限要素法（FEM）メッシュとノイズを含むデータを用いたシミュレーションでも、安定した近似再構成が得られることを示しました。

4. 意義と考察 (Significance)

• 理論的進展: 線形化カルデロン問題における一意性（injectivity）が $L^3$ 空間でも成り立つことを示し、その具体的な再構成アルゴリズムを構築しました。これは、より一般的な物理モデルへの適用可能性を開くものです。
• 実用性: 前方代入法であるため、大規模な逆問題ソルバー（反復法など）に比べて計算コストが低く、リアルタイムに近い処理や、限られた測定データからの再構成に有望です。
• 数値的安定性への課題: 3 次元では、2 次元には存在しなかった「追加の和（sum）」項（既知の係数の中で $\ell$ 指標が大きいものが必要になる）が存在します。これは球面調和関数の積が複数の次数の線形結合になることに起因し、数値的に不安定になる可能性がありますが、実験的には良好な結果が得られています。
• 将来展望: この手法は、不完全データ（partial data）や、より複雑な幾何学形状への拡張、あるいは非線形問題への適用への道筋を示唆しています。

総じて、本論文は 3 次元電気インピーダンストモグラフィの理論的基盤を強化し、非有界な導電率変化を効率的に画像化する新しい計算手法を確立した重要な研究です。