Returning walks on a lattice are sequences of moves that start at a given lattice site and return to the same site after $n$ steps. Determining the total number of returning walks of a given length $n$ is a typical graph-theoretical problem with connections to lattice models in statistical and condensed matter physics. We derive analytical expressions for the returning walk numbers on the eleven two-dimensional Archimedean lattices by developing a connection to the theory of Bloch energy bands. We benchmark our results through an alternative method that relies on computing the moments of adjacency matrices of large graphs, whose construction we explain explicitly. As condensed matter physics applications, we use our formulas to compute the density of states of tight-binding models on the Archimedean lattices and analytically determine the asymptotics of the return probability. While the Archimedean lattices provide a sufficiently rich structure and are chosen here for concreteness, our techniques can be generalized straightforwardly to other two- or higher-dimensional Euclidean lattices.
論文ID : 2507.12662タイトル : Walking on Archimedean Lattices: Insights from Bloch Band Theory著者 : Davidson Noby Joseph、Igor Boettcher(アルバータ大学)分類 : cond-mat.stat-mech、cond-mat.mes-hall、cond-mat.str-el、math-ph、math.MP発表時期 : 2025年1月(arXiv プレプリント)論文リンク : https://arxiv.org/abs/2507.12662 本論文は格子上の帰還歩行問題を研究しており、与えられた格子点から出発してn歩後に原点に戻る経路の数列を扱っている。著者らはブロッホ帯理論との関連性を確立することで、11種類の2次元アルキメデス格子上の帰還歩行数の解析表現式を導出した。研究結果は大規模グラフの隣接行列の累乗を計算する代替方法により検証されている。凝聚系物理学への応用として、著者らはこれらの公式を用いてアルキメデス格子上の強束縛模型の状態密度を計算し、帰還確率の漸近挙動を解析的に決定した。
中心的問題 : 格子上の長さnの帰還歩行の総数を決定することは、統計物理学および凝聚系物理学の格子模型において重要な応用を持つ典型的なグラフ理論問題である。物理的意義 :帰還歩行数は状態密度(DOS)と直接関連しており、後者は電子材料の物性を記述する中核的な物理量である 帰還確率は統計模型の過渡的挙動およびランダム固体中のアンダーソン局在において重要な役割を果たす 格子量子模型の離散経路積分を定義するために使用できる 既存方法の限界 :複雑な格子構造に対して、単純な組合せ論的方法は適用困難である 従来の直接計算方法は高次元または大きな単位胞の場合、計算複雑性が過度に高い 異なる種類の周期的タイリングを処理する統一的な理論枠組みが欠けている 研究動機 :アルキメデス格子は最も研究されている平面格子(正方、三角、ハニカム、カゴメ格子など)を含んでいる グラフ理論問題と物理的帯理論の橋渡しを確立し、計算のための新しい理論的ツールを提供する より一般的な2次元または高次元ユークリッド格子に対して推広可能な技術を提供する 帰還歩行数とブロッホ帯理論の根本的な関連性を確立 し、中核公式を導出:
S n = 1 N u ∫ k Tr [ A ( k ) n ] S_n = \frac{1}{N_u}\int_k \text{Tr}[A(k)^n] S n = N u 1 ∫ k Tr [ A ( k ) n ] 11種類すべてのアルキメデス格子の帰還歩行数を体系的に計算 し、明示的な公式または生成関数を提供大規模有限グラフ(フレークおよびクラスタ)を構築するための汎用方法を開発 し、開放境界条件および周期的境界条件に対応7種類のアルキメデス格子の状態密度の解析表現式を計算 し、その一部は初めて報告される帰還確率の漸近挙動を決定 し、p n ∼ α / n p_n \sim \alpha/n p n ∼ α / n の普遍的形式を得た周期的格子が与えられたとき、帰還歩行数をS n ( i ) = ( A n ) i i S_n^{(i)} = (A^n)_{ii} S n ( i ) = ( A n ) ii と定義する。ここでA A A は隣接行列である。頂点推移的なアルキメデス格子に対して、S n = S n ( i ) S_n = S_n^{(i)} S n = S n ( i ) は開始点に依存しない。
N u N_u N u 個の原子を持つ単位胞を有する周期的タイリングに対して、ブロッホ隣接行列A ( k ) A(k) A ( k ) はN u × N u N_u \times N_u N u × N u 行列であり、その行列要素は以下の通りである:
A ( k ) i u j u = ∑ v ⟨ i u ∣ A ^ T ^ v ∣ j u ⟩ e − i v ⋅ k A(k)_{i_u j_u} = \sum_v \langle i_u|\hat{A}\hat{T}_v|j_u\rangle e^{-iv \cdot k} A ( k ) i u j u = ∑ v ⟨ i u ∣ A ^ T ^ v ∣ j u ⟩ e − i v ⋅ k
厳密な数学的導出を通じて、以下が証明される:
S n ( i u ) = ∫ k [ A ( k ) n ] i u i u S_n^{(i_u)} = \int_k [A(k)^n]_{i_u i_u} S n ( i u ) = ∫ k [ A ( k ) n ] i u i u
この恒等式は無限格子上の組合せ論的問題を、ブリルアン領域上の有限次元行列の積分に変換する。
ブロッホ生成関数を定義する:
G ( z , k ) = 1 N u Tr ( 1 1 − z A ( k ) ) G(z,k) = \frac{1}{N_u}\text{Tr}\left(\frac{1}{1-zA(k)}\right) G ( z , k ) = N u 1 Tr ( 1 − z A ( k ) 1 )
すると以下が成立する:
G ( z ) = ∫ k G ( z , k ) = ∑ n ≥ 0 S n z n G(z) = \int_k G(z,k) = \sum_{n \geq 0} S_n z^n G ( z ) = ∫ k G ( z , k ) = ∑ n ≥ 0 S n z n
理論的突破 : グラフ理論における帰還歩行問題と凝聚系物理学におけるブロッホ帯理論の深層的な関連性を初めて体系的に確立計算方法の革新 :運動量積分を評価するための純代数的技術を開発し、定数部分抽出法を使用 複雑な格子に対して、固有値の明示的な計算を回避するため、行列逆の対角要素を計算 構成方法 : 任意の大きさの有限グラフを構築するための汎用主公式(56)を提案:
A p , q = ( 1 p ⊗ 1 q ) ⊗ A u + ∑ d [ ( R p ( d ) ⊗ R q ( d ) ) ⊗ Γ ( d ) + h.c. ] A_{p,q} = (1_p \otimes 1_q) \otimes A_u + \sum_d [(R_p^{(d)} \otimes R_q^{(d)}) \otimes \Gamma^{(d)} + \text{h.c.}] A p , q = ( 1 p ⊗ 1 q ) ⊗ A u + ∑ d [( R p ( d ) ⊗ R q ( d ) ) ⊗ Γ ( d ) + h.c. ] 11種類のアルキメデス格子 :
二部格子 : 正方(Square)、ハニカム(Honeycomb)、CaVO、SHD非二部格子 : 三角(Triangular)、カゴメ(Kagome)、トレリス(Trellis)、スター(Star)、SrCuBO、ルビー(Ruby)、メープルリーフ(Maple-Leaf)理論計算 : ブロッホ隣接行列および生成関数法を使用数値検証 : 大規模有限グラフを構築し、隣接行列の累乗を計算収束性検査 : p , q p, q p , q パラメータを増加させることでS n S_n S n の収束性を検査Mathematicaなどの計算代数ソフトウェアを使用して行列逆を計算 ケーリー・ハミルトン定理を使用して前数項の矩から生成関数を計算 ブリルアン領域積分に2500個の均一サンプリング点を採用 11種類すべての格子の帰還歩行数列の計算に成功。例えば:
正方格子 : S 2 n = ( 2 n n ) 2 S_{2n} = \binom{2n}{n}^2 S 2 n = ( n 2 n ) 2 ハニカム格子 : S 2 n = ∑ l = 0 n ( 2 l l ) ( n l ) 2 S_{2n} = \sum_{l=0}^n \binom{2l}{l}\binom{n}{l}^2 S 2 n = ∑ l = 0 n ( l 2 l ) ( l n ) 2 三角格子 : S n = ∑ l = 0 n ( n l ) ( − 3 ) n − l S 2 l ( H ) S_n = \sum_{l=0}^n \binom{n}{l}(-3)^{n-l}S_{2l}^{(H)} S n = ∑ l = 0 n ( l n ) ( − 3 ) n − l S 2 l ( H ) 7種類の格子の状態密度解析表現式を取得。以下を含む:
正方格子 : D □ ( E ) = 1 2 π 2 K ( 1 − E 2 16 ) D_\square(E) = \frac{1}{2\pi^2}K(1-\frac{E^2}{16}) D □ ( E ) = 2 π 2 1 K ( 1 − 16 E 2 ) ハニカム格子 : D H ( E ) = 2 ∣ E ∣ π 2 Z 0 ( E ) K ( Z 1 ( E ) Z 0 ( E ) ) D_H(E) = \frac{2|E|}{\pi^2\sqrt{Z_0(E)}}K(\frac{Z_1(E)}{Z_0(E)}) D H ( E ) = π 2 Z 0 ( E ) 2∣ E ∣ K ( Z 0 ( E ) Z 1 ( E ) ) 帰還確率の漸近係数α \alpha α を決定:
二部格子 : S 2 n ∼ α q 2 n 2 n S_{2n} \sim \alpha \frac{q^{2n}}{2n} S 2 n ∼ α 2 n q 2 n 非二部格子 : S n ∼ α q n n S_n \sim \alpha \frac{q^n}{n} S n ∼ α n q n 具体的な数値として、正方格子ではα = 2 π \alpha = \frac{2}{\pi} α = π 2 、三角格子ではα = 3 2 π \alpha = \frac{\sqrt{3}}{2\pi} α = 2 π 3 など。
数値的一貫性 : 理論計算結果は大規模有限グラフの数値計算と完全に一致既知結果の再現 : 正方、ハニカム、三角、カゴメ格子の既知数列を成功裏に再現新規数列の発見 : 残りの7種類の格子に対して、完全な帰還歩行数列を初めて計算組合せ論的方法 : 正方格子などの単純な格子にのみ適用可能隣接行列法 : A n A^n A n を直接計算するが、大規模システムでは計算困難連分数法 : ランツォス法により構築するが、数値丸め誤差が存在スペクトル矩定理 : 本論文の結果は普遍的スペクトル矩定理の特殊ケースと見なせる経路積分 : 離散経路積分理論と関連超幾何関数 : 楕円積分および超幾何関数理論と密接に関連統一的枠組み : すべてのアルキメデス格子を処理するための統一的方法を提供解析的結果 : 数値近似ではなく、正確な解析表現式を取得物理的洞察 : グラフ理論と帯理論の深層的な関連性を明らかにする理論的貢献 : 帰還歩行問題とブロッホ帯理論の根本的な関連性を確立計算成果 : 11種類のアルキメデス格子の帰還歩行数を完全に計算応用価値 : 状態密度計算および漸近解析のための新しいツールを提供適用範囲 : 主にシンモルフィック結晶学的タイリングを対象とし、非周期構造には不適用計算複雑性 : 大きな単位胞(例えばSHDの12個の原子)に対して、解析計算は依然として複雑次元制限 : 理論的には高次元への推広が可能だが、具体的な計算複雑性は著しく増加推広応用 : ラベス格子およびその他の周期的タイリングへの拡張非ユークリッド幾何学 : 双曲格子などの非ユークリッド的タイリングへの応用関連問題 : 自己回避歩行、面積制約付き歩行などの関連組合せ論的問題への拡張物理的応用 : 2次元イジング模型の臨界温度および自由エネルギー計算への応用理論的革新性 : グラフ理論と帯理論の関連性を初めて体系的に確立し、重要な理論的価値を有する計算の完全性 : すべてのアルキメデス格子をカバーし、完全なデータセットを提供方法の普遍性 : 開発された技術は他の周期的タイリングに推広可能物理的意義 : 結果は状態密度計算などの凝聚系物理学の問題に直接応用可能数学的厳密性 : 導出は厳密であり、詳細な数学的証明を含む新規性の制限 : 中核恒等式(68)はある程度、既知のスペクトル矩定理の応用と見なせる計算効率 : 複雑な格子に対して、計算は依然として記号計算ソフトウェアに依存実験的検証 : 実際の物理実験との比較検証が欠けている応用の深さ : ツールを提供する一方で、具体的な物理問題への深い応用は限定的学術的価値 : グラフ理論、数学物理学、凝聚系物理学の学際的研究に新しい範例を提供実用的価値 : 材料科学における状態密度計算のための新しいツールを提供再現可能性 : 詳細な計算方法と結果を提供し、検証と応用を容易にする啓発的意義 : 関連する組合せ論的問題の研究に新しい思考方向を開く理論研究 : 格子模型、統計物理学、グラフ理論研究材料計算 : 新規2次元材料の電子構造計算アルゴリズム開発 : 大規模疎行列固有値問題の新しいアルゴリズム教育応用 : 数学物理学学際領域の典型的なケーススタディとして論文は古典的なポーリャ確率歩行理論から最新の双曲格子研究まで、75篇の関連文献を引用しており、この分野の歴史的発展と現在の最前線を反映している。重要な参考文献にはウォレスのグラフェン帯理論、キタエフの量子スピン液体模型、および最近の双曲格子実験研究などが含まれる。