[論文レビュー] On Helmholtz free energy for finite abstract simplicial complexes

Oliver Knill|arXiv (Cornell University)|Mar 20, 2017

Topological and Geometric Data Analysis参考文献 18被引用数 4

ひとこと要約

本稿は、有限抽象単体複体の接続ラプラシアンのグリーン関数を介して、ヘルムホルツ自由エネルギーとオイラー乗数の間の位相的同一性を確立する。全ポテンシャルエネルギー ∑ₓ,ᵧg(x,y) が χ(G) に等しいことを示し、ヘルムホルツ自由エネルギー F(p) = βU(p) − TS(p) の臨界点が分岐と不連続性を示すことを明らかにした。これは、K₃ などの小さな複体でさえも、相転移に類似した現象を示すことを示している。

ABSTRACT

We prove a Gauss-Bonnet formula X(G) = sum_x K(x), where K(x)=(-1)^dim(x) (1-X(S(x))) is a curvature of a vertex x with unit sphere S(x) in the Barycentric refinement G1 of a simplicial complex G. K(x) is dual to (-1)^dim(x) for which Gauss-Bonnet is the definition of Euler characteristic X. Because the connection Laplacian L'=1+A' of G is unimodular, where A' is the adjacency matrix of of the connection graph G', the Green function values g(x,y) = (1+A')^-1_xy are integers and 1-X(S(x))=g(x,x). Gauss-Bonnet for K^+ reads therefore as str(g)=X(G), where str is the super trace. As g is a time-discrete heat kernel, this is a cousin to McKean-Singer str(exp(-Lt)) = X(G) for the Hodge Laplacian L=dd^* +d^*d which lives on the same Hilbert space than L'. Both formulas hold for an arbitrary finite abstract simplicial complex G. Writing V_x(y)= g(x,y) for the Newtonian potential of the connection Laplacian, we prove sum_y V_x(y) = K(x), so that by the new Gauss-Bonnet formula, the Euler characteristic of G agrees with the total potential theoretic energy sum_x,y g(x,y)=X(G) of G. The curvature K now relates to the probability measure p minimizing the internal energy U(p)=sum_x,y g(x,y) p(x) p(y) of the complex. Since both the internal energy (here linked to topology) and Shannon entropy are natural and unique in classes of functionals, we then look at critical points p the Helmholtz free energy F(p)=(1-T) U(p)-T S(p) which combines the energy functional U and the entropy functional S(p)=-sum_x p(x) log(p(x)). As the temperature T changes, we observe bifurcation phenomena. Already for G=K_3 both a saddle node bifurcation and a pitchfork bifurcation occurs. The saddle node bifurcation leads to a catastrophe: the function T -> F(p(T),T) is discontinuous if p(T) is a free energy minimizer.

研究の動機と目的

• 有限抽象単体複体における接続ラプラシアンのグリーン関数成分の位相的解釈を確立すること。
• 新しいガウス-ボンネット型の公式を介して、全ポテンシャルエネルギー ∑ₓ,ᵧg(x,y) をオイラー乗数 χ(G) に結びつけること。
• 内部エネルギー U(p) とシャノンエントロピー S(p) を用いたヘルムホルツ自由エネルギー F(p) = βU(p) − TS(p) の臨界点を調査すること。
• 逆温度 β の関数としての自由エネルギー関数の分岐現象と不連続性を探索すること。
• オイラー乗数とシャノンエントロピーが、それぞれの定義域において正規化され乗法的関数として一意に特徴づけられることを示すこと。

提案手法

• 頂点 x に対する曲率 K⁺(x) = 1 − χ(S(x)) が、Barycentric refinement G₁ においてガウス-ボンネット公式 ∑ₓ K⁺(x) = χ(G) を満たすことを証明する。ここで S(x) は x の単位球面を表す。
• 接続ラプラシアン L′ = 1 + A′ がユニモジュラーであることを利用し、そのグリーン関数成分 g(x,y) = (1 + A′)⁻¹ₓᵧ が整数であることを示す。特に対角成分 g(x,x) = 1 − χ(S(x)) である。
• グリーン関数のスレッドが満たす str(g) = ∑ₓ g(x,x) = χ(G) を確立し、トレースとオイラー乗数の間の関係を結ぶ。
• 全ポテンシャルエネルギー ∑ₓ,ᵧ g(x,y) = χ(G) を示し、エネルギーと位相の間の直接的関係を確立する。
• ヘルムホルツ自由エネルギー F(p) = βU(p) − TS(p) を分析し、U(p) = ∑ₓ,ᵧ g(x,y)p(x)p(y) および S(p) = −∑ₓ p(x) log p(x) であることを用い、β に関する臨界点を検討する。
• G = K₃ の場合、自由エネルギー関数が逆温度 β の変化に伴い、サドルノード分岐とピッチフォーク分岐の両方を示すことを証明する。

実験結果

リサーチクエスチョン

• RQ1 任意の有限抽象単体複体に対して、接続ラプラシアンのグリーン関数の全ポテンシャルエネルギー ∑ₓ,ᵧ g(x,y) はオイラー乗数 χ(G) に等しいか？
• RQ2 曲率 K⁺(x) = 1 − χ(S(x)) はグリーン関数のスーパitraceとして解釈可能であり、ガウス-ボンネット公式を満たすか？
• RQ3 逆温度 β の変化に伴い、ヘルムホルツ自由エネルギー F(p) = βU(p) − TS(p) の臨界点にどのような分岐現象が現れるか？
• RQ4 複体の位相的不変量、たとえば χ(G) は、エネルギーとエントロピーを含む熱力学的関数とどのように関係するか？
• RQ5 オイラー乗数とシャノンエントロピーは、それぞれ単体複体と確率測度の定義域において、正規化され乗法的関数として一意に特徴づけられるか？

主な発見

• 接続ラプラシアンのグリーン関数の全ポテンシャルエネルギー ∑ₓ,ᵧ g(x,y) はオイラー乗数 χ(G) に等しく、エネルギーと位相の間の直接的同一性を確立する。
• 曲率 K⁺(x) = 1 − χ(S(x)) はガウス-ボンネット公式 ∑ₓ K⁺(x) = χ(G) を満たし、K⁺(x) はグリーン関数のスーパitraceとして解釈可能である。
• G = K₃ の場合、ヘルムホルツ自由エネルギー関数は逆温度 β の変化に伴い、サドルノード分岐とピッチフォーク分岐の両方を示す。
• 自由エネルギー F(p(β)) は臨界 β 値で不連続となり、系の熱力学的挙動における崩壊を示している。
• 内部エネルギー関数 U(p) = ∑ₓ,ᵧ g(x,y)p(x)p(y) とシャノンエントロピー関数 S(p) は、それぞれの定義域において正規化され乗法的関数として一意に特徴づけられる。
• グリーン関数の対角成分は g(x,x) = 1 − χ(S(x)) を満たし、非対角成分はエネルギー-位相同一性における役割を通じて位相的意義を持つことが示された。