Improved stability of low Fourier modes in inverse problems for potentials

技術要約：ポテンシャルに関する逆問題における低次フーリエモードの安定性の向上

問題設定
本研究は、楕円型方程式の逆問題、特にディリクレ・ネーマン（DN）写像 $\Lambda_q$ から未知のポテンシャル $q$ を復元することに焦必している。このような逆問題における中心的な課題は、固有の不良設定性にあり、通常は測定ノイズに対して非常に敏感な対数的な安定性評価（logarithmic stability estimates）をもたらす。著者らは、ポテンシャルの特定の構造的特徴、すなわちその低周波フーリエモードが、ポテンシャルの差 $Q = q - q'$ に対する特定の正則性仮定の下で、改善された安定性（リプシッツ、サブ・ヘルダー、またはヘルダー）をもって復元できるかどうかを調査している。

本研究では、$Q$ の3つの異なるシナリオを検討している：

1. バンドリミテッド（帯域制限）な差： $Q$ は $M$-バンドリミテッドである（$|k| > M$ でフーリエモードが消失する）。
2. 実解析的な差： $Q$ は実解析的であり、かつ周期関数である。
3. 超指数関数的に減衰する差： $Q$ のフーリエ係数が $e^{-c|k|^{n/2}}$ の割合で減衰する。

支配方程式は、高周波 $\lambda$ におけるヘルムホルツ方程式 $(-\Delta - \lambda^2 + q)u = 0$ および $\lambda=0$ におけるシュレーディンガー方程式 $(-\Delta + q)u = 0$ である。

手法
証明は、DN写像の差と、前向き問題の解との積分を関連付けるアレッサンドリーニの恒等式（Alessandrini identity）に基づいている：
$$\langle (\Lambda_q - \Lambda_{q'})u, u' \rangle_{\partial \Omega} = \int_\Omega (q - q') u \bar{u}' \, dx$$
フーリエ係数 $Q$ の情報を抽出するために、著者らは特定のフーリエモードを分離するために高周波で振動する特定のテスト関数（解 $u, u'$）を構成する。

• 実幾何光学（RGO）： バンドリミテッドの場合（定理 1.1）において、著者らは $u = e^{i\eta \cdot x}(1+r)$ （ここで $|\eta| = \lambda$）という形の、実幾何光学解を用いる。これらの解は、高周波におけるポテンシャルの差を探索するために用いられる。解析には、剰余項 $r$ の推定、およびフーリエ係数とDN写像の差との関係を反転させるための、ランダウ–ポラック–スレピアン型の作用素の正定値性の利用が含まれる。
• 複素幾何光学（CGO）： 実解析および超指数減衰の場合（定理 1.2 および 1.3）において、著者らは複素位相 $\rho, \rho' \in \mathbb{C}^n$ を持つCGO解を用いる。これにより、周期的な設定 $\Omega = (-\pi, \pi)^n$ においてフーリエモードを分離することが可能となる。
• フーリエ減衰の推定： 手法における重要な要素は、フーリエ級数の「テイル（裾）」（モード $|k| > M$）を抑え込むことである。
  • 実解析関数については、著者らは複素近傍への正則拡張に対するコーシーの評価を利用して、フーリエ係数の指数関数的減衰を導出する。
  • 超指数減衰については、$|Q_k| \leq C e^{-c|k|^{n/2}}$ という減衰を仮定する。
  • 補題 3.2 は、高周波シェルにおける減衰する指数関数のテイルの和に関する推定を提供しており、これは安定性評価の誤差項を制御するために不可欠である。

主な貢献と結果

1. バンドリミテッドな差に対するリプシッツ安定性（定理 1.1）：
   差 $Q = q - q'$ が $M$-バンドリミテッドである場合、周波数 $\lambda$ が帯域幅に対して十分に大きい（$\lambda \gtrsim M^{n/2}$）ならば、そのすべてのフーリエモードはリプシッツ安定性をもって復元できることを著者らは証明している。

   • 意義： 次元に応じて指数関数的に増大する一般的な有限次元のリプシッツ結果とは異なり、ここでのリプシッツ定数は $M$ に対して一様である（$\Omega, n, E$ にのみ依存する）。
   • 仮定： バンドリミテッドである必要があるのは「差」 $Q$ だけであり、個々のポテンシャル $q, q'$ は $L^\infty$ で有界であるだけでよい。

2. 実解析的な差に対するサブ・ヘルダー安定性（定理 1.2）：
   $Q$ が実解析的かつ周期的である場合、著者らはサブ・ヘルダー安定性を確立している。安定性係数は $\exp(-c(\log(1/\delta))^{2/n})$ のように振る舞い、これは対数的よりも良く、任意のヘルダー係数 $\delta^\alpha$ よりも劣る。

   • 解決策： 復元可能なフーリエモードの数 $M(\delta)$ は、測定誤差 $\delta$ が減少するにつれて増大する。具体的には、$M(\delta) \sim (\log(1/\delta))^{2/n}$ である。
   • 仮定： $Q$ は半径 $\epsilon$ の複素近傍に正則に拡張できる。

3. 超指数減衰に対するヘルダー安定性（定理 1.3）：
   もし $Q$ のフーリエ係数が超指数的な速度 $e^{-c|k|^{n/2}}$ で減衰する場合、安定性はヘルダー安定性（$\delta^\alpha$）へと向上する。

   • 指数： ヘルダー指数は $\alpha = c/(\gamma + c)$ である（ここで $c$ は減衰率、$\gamma$ は幾何学的定数）。$c \to \infty$ となるにつれ、$\alpha \to 1$ となり、リプシッツ安定性に近づく。
   • 解決策： 解析的な場合と同様に、復元可能なモードの数 $M(\delta)$ は $\delta \to 0$ のとき無限大に発散する。具体的には、$M(\delta) \sim (\log(1/\delta))^{2/n}$ である。

意義と主張
本論文は、「増加する安定性（increasing stability）」および「増加する分解能（increasing resolution）」という非線形逆問題における理解を深めるものであると主張している。

• 仮定の緩和： 主要な貢献の一つは、正則性の仮定（バンドリミテッド、解析的、または超指数減衰）が、個々のポテンシャルではなく、ポテンシャルの差 $Q$ に対して適用される点である。ポテンシャル自体は $L^\infty$ で有界であることのみが要求される。これにより、未知数が既知の有限次元部分空間に属さなければならないという制限的な仮定を回避している。
• 一様な定数： バンドリミテッドなケースにおいて、安定性定数は未知数の次元（モードの数 $M$）に対して指数関数的に増大しない。これは、標準的な有限次元の逆問題の結果とは対照的である。
• 定量的な分解能： 本結果は、測定精度（$\delta$）と分解能（復元可能なフーリエモードの数）との間の定量的な関係を提供している。測定がより精密になるにつれて、より多くの低次モードが安定して復元可能となり、差の正則性に応じて、対数的からヘルダー／リプシッツ安定性の領域へと移行する。
• 非線形性： 本研究は、同様の増加する安定性が観察されていた線形化された設定（線形散乱における特異値推定など）を超えて、非線形逆問題を直接扱っている。

著者らは、推定における定数は復元されるモード数に対して一様であり、結果は一般的な有界領域（定理 2.1）またはフーリエ解析が最も自然なトーラス／立方体の設定に対して成立することを述べている。本論文は実験的な応用を提案するものではなく、数学的な逆問題に対する理論的な安定性評価に焦点を当てている。