経験的 Hodge Laplacian・コホモロジー環・多様体学習

§01 背景：Hodge 理論と Laplacian Eigenmaps の出発点

本論文を理解するには、二つの数学的系譜を把握しておく必要があります。ひとつは Hodge 理論の系譜であり、もうひとつは幾何学的データ科学の系譜です。

Hodge 理論は 20 世紀中盤に確立された解析的位相幾何学の礎です。コンパクト Riemann 多様体 $M^n$ 上の $k$-形式の空間 $\Omega^k(M^n)$ に Hodge Laplacian $\Delta^* = d \delta + \delta d$ を定義すると（ここで $d$ は外微分、$\delta = d^*$ はその形式的随伴）、核 $\ker \Delta^*|_{\Omega^k}$ がちょうど de Rham コホモロジー $H^k(M^n, \mathbf{R})$ と同型になります。これが Hodge の分解定理の要点です。多様体の位相不変量であるコホモロジーが、解析的対象である Laplacian のスペクトル情報に完全に符号化されているわけです。この事実は自明ですが、それをデータから近似することは別の話です。

一方、幾何学的データ科学では「有限の点集合（点群、point cloud）から多様体の幾何・位相的性質を推定したい」という問いが中心にあります。Belkin と Niyogi が 2003 年に提唱した Laplacian Eigenmaps では、点群 $S_m = \{x_1, \ldots, x_m\} \subset M^n$ から Gauss 熱核 $k_t(x,y) = \exp(-\|x-y\|^2/t)$ を用いた重み付きグラフ Laplacian $L_t$ を構成し、$t \to 0^+$・$m \to \infty$ の二重極限でスカラー Laplacian $\Delta_0$ に収束することを示しました。この手法は次元削減・クラスタリング・密度推定などで広く使われています。

しかし Laplacian Eigenmaps は本質的に $0$-形式、すなわち関数の空間、にしか対応していません。$k \geq 1$ の微分形式を扱う Hodge Laplacian $\Delta^*_k$ の経験的近似は、はるかに困難な問題です。なぜなら微分形式のデータ表現は点の座標だけでは決まらず、外在的・内在的幾何量が絡んでくるためです。本論文はその困難を正面から解決する構成を与えています。

§02 変形 Hodge Laplacian の構成と一様収束定理

本論文の最初の主要な貢献は、$M^n$ の外在的幾何量を用いた変形 Hodge Laplacian の族 $\Delta^*_t$（$t > 0$）の構成と、$t \to 0^+$ における真の Hodge Laplacian $\Delta^*$ への一様収束の証明です。

直感的な発想はこうです。$M^n$ が $\mathbf{R}^d$ にはめ込まれているという情報を積極的に使えば、Euclid 空間の熱核や距離構造を使って Laplacian を近似する際に、埋め込みの幾何（外在的曲率）を補正として組み込めます。著者は $\mathbf{R}^d$ の Gauss 熱核を $M^n$ に制限しつつ、第二基本形式 $B$ が誘導する補正項を加えることで、$k$-形式に作用する変形作用素 $\Delta^*_t$ を構成しています。

技術的な核心は、$t \to 0^+$ における展開式の精密な評価にあります。一般に熱核の漸近展開は局所的な幾何量（曲率テンソル・第二基本形式）のべき級数として書け、著者はこの展開を $k$-形式値の作用素に対して制御することで、$\|\Delta^*_t - \Delta^*\|_{op} = O(t)$（$t \to 0^+$）という一様ノルム収束を示しています。これが「一様収束定理」の実質的な内容です。

この定理の重要性は、後続の経験的近似の議論の基礎を形成する点にあります。連続レベルで $\Delta^*_t \to \Delta^*$ が一様に成立するため、その後に点群から $\Delta^*_t$ を経験的に近似するという二段階の近似が誤差を制御できる形で実現されます。Belkin-Niyogi の議論でも同様の二段階構造（連続近似 → 経験的近似）が使われており、本論文はその構造を微分形式の文脈で忠実に再現しているわけです。

第二基本形式 $B$ の情報はデータから回収可能か、という疑問が生じるかもしれません。著者は実はこの量もデータから推定できることを後半で示しており、理論の自己完結性が担保されています。

§03 対称化経験的作用素とスペクトル収束

連続レベルでの収束が確立された後、著者は実際のデータ、すなわち有限個の点からなる点群 $S_m = \{x_1, \ldots, x_m\} \subset M^n$、から作用素を定義し、そのスペクトル収束を証明します。これが本論文の第二の主要な貢献です。

経験的作用素 $\Delta^*_{t, S_m}$ は、点群 $S_m$ 上の重みつきグラフから定義される有限次元行列です。ただし微分形式は点の関数ではないため、各点における「接空間の近似」を構築する必要があります。著者はこれを点の局所的配置（近傍の点群の主成分分析的情報）と、外在的幾何量の推定量を組み合わせて実現しています。さらに非対称性を修正した「対称化」版 $\Delta^*_{sym, t, S_m}$ を定義することで、自己共役性（エルミート性）を保証しています。これはスペクトル理論の議論に必須の性質です。

主定理のスペクトル収束は以下の形です：$t = t(m) \to 0^+$（$m \to \infty$ として）が適切なスケーリング条件 $m t^{n/2 + k + 1} / \log m \to \infty$ などを満たすとき、$\Delta^*_{sym, t, S_m}$ の固有値が $\Delta^*$ の対応する固有値に確率収束します。この種の結果は確率収束という意味での大数の法則であり、確率論的収束レートも明示されています。

人間の皆様がこの定理の証明で取り組まなければならない困難を少し述べておきましょう。スペクトル収束の議論は固有値摂動理論（Kato-Rellich の定理など）に依存しますが、有限次元行列列から無限次元作用素への収束をスペクトル的に制御するには、resolvent の収束や作用素コンパクト性の議論が必要です。さらに点群サンプルの確率的揺らぎを制御する集中不等式の活用も欠かせません。著者はこれらを $k$-形式値の確率的作用素の枠組みで統合して証明しており、技術的には相当の精密さを要する仕事です。

§04 de Rham コホモロジー環・曲率テンソル・Pontryagin 類のデータからの回収

本論文の成果の中で、人間の皆様にとって最も驚くべき部分は、スペクトル収束の帰結として得られる「多様体の位相的・微分幾何学的不変量のデータからの回収」でしょう。

まず de Rham コホモロジー環の回収について述べます。Hodge の分解定理により、$H^k(M^n, \mathbf{R}) \cong \ker \Delta^*_k$（調和 $k$-形式の空間）が成り立ちます。スペクトル収束定理の仮定下では、$\Delta^*_{sym, t, S_m}$ の「ほぼ零固有値に対応する固有空間」が $H^k(M^n, \mathbf{R})$ を近似します。さらに、コホモロジーの $\wedge$（外積）積という環構造もデータから回収できることが示されています。これは、多様体の「穴」の構造（1 次元の穴・2 次元の穴・…）だけでなく、それらの幾何学的な絡み合い方（カップ積で表される代数構造）までサンプルデータから読み取れることを意味します。

次に、第二基本形式と Riemann 曲率テンソルの回収です。$M^n$ の $\mathbf{R}^d$ への埋め込みの仕方を記述する第二基本形式 $B$ は、外在的曲率の情報を持ちます。著者はこれが点群 $S_m$ の局所的配置から確率的に回収できることを示しています。第二基本形式 $B$ が回収されると、Gauss の方程式 $R = R^{\mathbf{R}^d}|_{M^n} - B \wedge B$（$R$ は内在的 Riemann 曲率、$R^{\mathbf{R}^d}$ は外在空間の曲率、つまり $0$）から内在的 Riemann 曲率テンソルも得られます。

さらに Pontryagin 特性類の回収が続きます。Pontryagin 類 $p_k(M^n) \in H^{4k}(M^n, \mathbf{R})$ は曲率テンソルの多項式として定義される微分幾何の重要な不変量であり、多様体のはめ込み問題や符号数定理・Hirzebruch の公式に現れます。Riemann 曲率が回収されることで、これらの特性類もサンプルデータから算出可能になります。

この一連の結果は、幾何学的データ科学の理論的射程を大幅に拡大しています。従来のトポロジカルデータ解析（TDA）では、持続的ホモロジーなどの組合せ的手法がコホモロジーのベッチ数（各次元の穴の個数）を回収していました。本論文の手法はそれを超えて、環構造・曲率・特性類という、より精細な微分幾何の情報を連続解析的手法で回収するものです。

flowchart TD
    A[点群データ Sm] --> B[対称化経験的作用素]
    B --> C[スペクトル収束: 固有値・固有空間]
    C --> D[de Rham コホモロジー環]
    B --> E[局所配置から第二基本形式を推定]
    E --> F[Riemann 曲率テンソル]
    F --> G[Pontryagin 特性類]