[論文レビュー] The Dirac operator of a graph

Oliver Knill|arXiv (Cornell University)|Jun 10, 2013

Matrix Theory and Algorithms参考文献 11被引用数 26

ひとこと要約

本稿では、元のグラフのすべてのクリーク（完全部分グラフ）からなる単体グラフに基づいて、符号付き隣接行列を構築することで、有限単純グラフに対する離散ディラック作用素 D を導入する。向き付けられたクリークと離散外微分 d を用いて、作用素 D = d + d* はグラフ微積分（勾配、回転、発散）を符号化し、D² はラプラス＝ベルトラミ作用素 L を与える。主な貢献は、基本的な行列演算とグラフデータ構造のみを用いて、グラフコホモロジーとスペクトル解析を計算的に効率的に行うための線形代数ベースのフレームワークを提供することである。

ABSTRACT

We discuss some linear algebra related to the Dirac matrix D of a finite simple graph G=(V,E).

研究の動機と目的

• 有限単純グラフに対する微分幾何学のディラック作用素の離散的類似を構築すること。
• 向き付けられたクリークと離散外微分 d を用いて、グラフ微積分（勾配、回転、発散）を形式化すること。
• ディラック作用素 D = d + d* が異なる向き付けに対してユニタリ同値であること、したがってスペクトル特性が保存されることを示すこと。
• 標準的な線形代数とグラフデータ構造のみを用いて、コホモロジーとラプラシアンを実装可能な実用的手法を提供すること。
• D のスペクトルと L = D² のブロックが、向き付けの選択に依存せず、物理におけるゲージ不変性に類似していることを示すこと。

提案手法

• グラフ G のすべての K_{k+1} 部分グラフとして、クリークの集合 G_k を定義し、これにより単体グラフ G を形成する。
• 各クリークに向き付けを割り当て、次元 k ごとに符号付きインシデント行列 d_k を定義する。
• D = d + d* としてディラック作用素 D を構築する。ここで d は下三角行列で d² = 0 を満たし、クリーク上の関数空間の直和に作用する。
• 行列 |D_ij| を隣接行列として用いて単体グラフを定義し、符号は向きの整合性によって決定する。
• クリークの列挙と符号付きインシデント行列を用いて、エッジと頂点のリストから D を計算する Mathematica 関数 Dirac[s] を実装する。
• ラプラス＝ベルトラミ作用素を L = D² として計算し、これは k-形式（k-クリーク上の関数）に作用するブロック L_k に分解される。

実験結果

リサーチクエスチョン

• RQ1 組合せ的および線形代数的道具のみを用いて、有限単純グラフ上でディラック作用素をどのように離散化できるか。
• RQ2 クリークの向き付けとディラック作用素 D のスペクトル特性との関係は何か。
• RQ3 純粋に組合せ的設定において、可定向性を要件とせず、離散外微分 d が d² = 0 を満たすことは可能か。
• RQ4 ラプラス＝ベルトラミ作用素 D² は、スカラーラプラシアン L₀ = B - A などの既知のグラフラプラシアンとどのように関係するか。
• RQ5 グラフコホモロジーとスペクトル不変量は、基本的な行列演算とグラフデータ構造のみを用いて、どの程度効率的に計算可能か。

主な発見

• ディラック作用素 D は、異なるクリークの向き付けに対してユニタリ同値であり、対角成分が ±1 の行列による共役で与えられるため、D のスペクトルは向き付けに依存しない。
• 例題グラフにおける D の特性多項式は p_D(x) = x¹⁸ - 24x¹⁶ + 242x¹⁴ - 1334x¹² + 4377x¹⁰ - 8706x⁸ + 10187x⁶ - 6370x⁴ + 1624x² であり、正の固有値は 0.92, 1.05, 1.41, 1.69, 1.78, 2.00, 2.15, 2.38 である。
• スカラーラプラシアン L₀ の核は定数ベクトル [1,1,1,1,1,1,1]ᵀ で張られ、1次元サイクルにおける自明な0次コホモロジーを示している。
• L₁ の核は [1, -1, -3, 2, -5, 8, -8, 0, 8]ᵀ で張られ、グラフのサイクル構造に起因する非自明な1次コホモロジーを示している。
• ラプラス＝ベルトラミ作用素 L = D² は3つのブロック L₀（7×7）、L₁（9×9）、L₂（2×2）に分解され、L₂ = [[3,1],[1,3]] は三角形関数に作用する。
• L₂ の対角成分は 3 であり、p=2 で四面体が存在しない条件下で deg_p(x) = L_p(x,x) - (p+1) と整合しており、すべての三角形について deg(x) = 0 である。