A note on the spectral distribution of non-Hermitian block matrices with Toeplitz blocks

この論文「A note on the spectral distribution of non-Hermitian block matrices with Toeplitz blocks（Toeplitz ブロックを持つ非エルミットブロック行列のスペクトル分布に関する注）」の技術的な要約を以下に記します。

1. 研究の背景と問題設定

• 対象: 数値解析において頻繁に現れる、行列列 $\{A_n\}_n$ のスペクトル分布（固有値の漸近的分布）の解析。
• 行列構造: $k \times k$ ($k \ge 2$) のブロック行列であり、各ブロックが Toeplitz 行列 $T_n(f_{i,j})$ で構成されている。
  $$A_n = \begin{bmatrix} T_n(f_{1,1}) & \cdots & T_n(f_{1,k}) \\ \vdots & \ddots & \vdots \\ T_n(f_{k,1}) & \cdots & T_n(f_{k,k}) \end{bmatrix}$$
  ここで、$f_{i,j} \in L^\infty([-\pi, \pi])$ は生成関数である。
• 課題:
  • 行列列がエルミットである場合、既知の理論（GLT 理論や特異値分布の理論）を用いて完全なスペクトル解析が可能である。
  • しかし、非エルミット（Non-Hermitian）であり、かつブロック間の次元比が特定の有理数・無理数に収束しない一般的なケースにおいて、固有値分布を記述する符号（Symbol）を明示的に求める手法は確立されていなかった。
  • 従来の GLT 理論は特異値分布には有効だが、非エルミット行列の固有値分布については直接的な記述を提供しない。

2. 手法と理論的枠組み

本研究は、以下の 2 つの主要な数学的ツールを組み合わせてアプローチしている。

1. 一般化局所 Toeplitz (GLT) 列の理論:
   • 行列列の漸近的なスペクトル分布を記述する関数（GLT 符号）の枠組み。
   • Toeplitz 行列、対角サンプリング行列、零分布列などの基本クラスが GLT 代数を構成し、その演算（和、積、逆行列など）が符号の対応関係を保つ性質を利用する。
2. 行列の幾何学的平均 (Geometric Mean):
   • 2 つの正定値エルミット行列 $A, B$ に対する幾何学的平均 $G(A, B) = A^{1/2}(A^{-1/2}BA^{-1/2})^{1/2}A^{1/2}$ の概念を導入。
   • 非エルミットブロック行列の固有値分布を、ブロック間の幾何学的平均を用いて変形・近似することで解析する。
   • 具体的には、各 Toeplitz ブロックを巡回行列（Circulant matrix）で近似し、その幾何学的平均の性質を利用することで、元の行列列と相似な行列列を構成する。

3. 主要な貢献と定理

論文の中心的な成果は、以下の定理 1.1（および定理 3.8）の証明である。

• 仮定:

  • 生成関数 $f_{i,j}$ は実数値であり、対角ブロック $f_{i,i}$ は実数値。
  • 隣接する非対角ブロック $f_{i, i+1}$ と $f_{i+1, i}$ について、その本質的下限（essential infimum）が正である（$\text{ess inf} > 0$）。

• 結論:
  行列列 $\{A_n\}_n$ の固有値分布は、以下の行列値関数 $F(\theta)$ （またはその相似変換 $\tilde{F}(\theta)$）によって記述される。
  $$\{A_n\}_n \sim_\lambda F(\theta) = \begin{bmatrix} f_{1,1}(\theta) & f_{1,2}(\theta) & \cdots \\ f_{2,1}(\theta) & f_{2,2}(\theta) & \cdots \\ \vdots & \vdots & \ddots \end{bmatrix}$$
  または、幾何学的平均を用いて対称化された形式：
  $$\tilde{F}(\theta) = \begin{bmatrix} f_{1,1} & \sqrt{f_{1,2}f_{2,1}} & \cdots \\ \sqrt{f_{1,2}f_{2,1}} & f_{2,2} & \cdots \\ \vdots & \vdots & \ddots \end{bmatrix}$$
  注：非対角要素 $f_{i,j}$ と $f_{j,i}$ の積の平方根が現れる。これは、非エルミット行列の固有値分布が、実質的には「対称化された」符号の固有値分布と一致することを示唆している。

• 証明の鍵:

  • フロベニウスノルムにおける Toeplitz 行列の巡回行列による近似。
  • 幾何学的平均の連続性と、巡回行列の可換性を利用した誤差評価（$o(\sqrt{n})$ のオーダー）。
  • 行列 $A_n$ を相似変換し、対角ブロックと非対角ブロック（幾何学的平均）からなる行列 $\hat{A}_n$ に近づけることで、GLT 理論を適用可能にする。

4. 数値実験と結果

理論的導出を検証するための数値実験が行われた。

• ケース $k=2, 3$: 様々な生成関数（三角関数、多項式など）を用いた実験。
• 結果:
  • 理論的に予測された符号 $F(\theta)$ の固有値と、実際に計算された行列 $A_n$ の固有値の実部が、$n$ が大きくなるにつれて極めて良く一致することを確認した。
  • 虚部は数値的に無視できるほど小さい（$O(10^{-14})$）場合が多いが、仮定を緩和したケース（$\text{ess inf} = 0$ となる場合）では虚部が現れるものの、$n$ の増加とともに減少する傾向（$O(1/n)$）が観測された。
  • 仮定 $\text{ess inf} > 0$ を満たさない「限界ケース」においても、同様の分布が観測されることから、仮定は技術的なものであり、本質的には緩和可能である可能性が示唆された。

5. 意義と今後の展望

• 理論的意義:
  • 非エルミットブロック Toeplitz 行列の固有値分布を、生成関数の幾何学的平均を用いて記述する初めての体系的なアプローチを提供した。
  • GLT 理論と行列の幾何学的平均を結びつけることで、非対称な構造を持つ行列列の解析手法を拡張した。
• 実用的意義:
  • 偏微分方程式の離散化などで生じる非対称な大規模行列のスペクトル特性を予測可能にする。
  • クリロフ部分空間法（Krylov methods）などの反復解法の収束速度予測や、前処理（Preconditioning）技術の設計に寄与する。
• 今後の課題（著者の提言）:
  • 仮定 $\text{ess inf} > 0$ を「ブロックが正定値である」というより自然な条件に緩和できるか（予想 5.1）。
  • 多変数符号や、ブロックサイズが異なる場合への拡張。
  • $L^1$ 関数への拡張や、具体的な線形方程式ソルバーへの応用。

総じて、この論文は非エルミットなブロック Toeplitz 行列のスペクトル分布を、幾何学的平均の概念を通じて明確に記述する画期的な結果であり、数値線形代数と関数解析の接点において重要な貢献を果たしている。