2025-11-18T09:58:13.321305

Walking on Archimedean Lattices: Insights from Bloch Band Theory

Joseph, Boettcher
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.
academic

Walking on Archimedean Lattices: Insights from Bloch Band Theory

基本信息

  • 论文ID: 2507.12662
  • 标题: Walking on Archimedean Lattices: Insights from Bloch Band Theory
  • 作者: Davidson Noby Joseph, Igor Boettcher (University of Alberta)
  • 分类: 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步后回到原点的路径序列。作者通过建立与Bloch能带理论的联系，推导出十一种二维阿基米德晶格上回归走步数量的解析表达式。研究结果通过计算大图邻接矩阵矩的替代方法进行验证。作为凝聚态物理应用，作者利用这些公式计算了阿基米德晶格上紧束缚模型的态密度，并解析确定了回归概率的渐近行为。

研究背景与动机

问题定义与重要性

  1. 核心问题: 确定晶格上长度为n的回归走步总数，这是一个典型的图论问题，在统计物理和凝聚态物理的晶格模型中有重要应用。
  2. 物理意义:
    • 回归走步数量与态密度(DOS)直接相关，后者是描述电子材料性质的核心物理量
    • 回归概率在统计模型的瞬态行为和无序固体中的Anderson局域化中起重要作用
    • 可用于定义晶格量子模型的离散路径积分
  3. 现有方法局限性:
    • 对于复杂晶格结构，简单的组合方法难以适用
    • 传统的直接计算方法在高维或大单元胞情况下计算复杂度过高
    • 缺乏统一的理论框架处理不同类型的周期性铺砌
  4. 研究动机:
    • 阿基米德晶格包含了最常研究的平面晶格(正方、三角、蜂窝、kagome晶格等)
    • 建立图论问题与物理能带理论的桥梁，为计算提供新的理论工具
    • 为更一般的二维或高维欧几里得晶格提供可推广的技术

核心贡献

  1. 建立了回归走步数量与Bloch能带理论的根本联系，推导出核心公式： Sn=1Nu∫kTr[A(k)n]S_n = \frac{1}{N_u}\int_k \text{Tr}[A(k)^n]
  2. 系统性地计算了所有十一种阿基米德晶格的回归走步数量，包括显式公式或生成函数
  3. 发展了构造大型有限图(flakes和clusters)的通用方法，支持开边界和周期边界条件
  4. 计算了七种阿基米德晶格的态密度解析表达式，其中部分为首次报道
  5. 确定了回归概率的渐近行为，得到pn∼α/np_n \sim \alpha/n的普遍形式

方法详解

任务定义

给定一个周期性晶格，定义回归走步数量Sn(i)=(An)iiS_n^{(i)} = (A^n)_{ii}，其中AA是邻接矩阵。对于顶点传递的阿基米德晶格，Sn=Sn(i)S_n = S_n^{(i)}与起始点无关。

核心理论框架

1. Bloch邻接矩阵构造

对于具有NuN_u个原子的单元胞的周期性铺砌，Bloch邻接矩阵A(k)A(k)是Nu×NuN_u \times N_u矩阵，其矩阵元为： A(k)iuju=∑v⟨iu∣A^T^v∣ju⟩e−iv⋅kA(k)_{i_u j_u} = \sum_v \langle i_u|\hat{A}\hat{T}_v|j_u\rangle e^{-iv \cdot k}

2. 关键恒等式推导

通过严格的数学推导，证明了： Sn(iu)=∫k[A(k)n]iuiuS_n^{(i_u)} = \int_k [A(k)^n]_{i_u i_u}

这个恒等式将无限晶格上的组合问题转化为有限维矩阵在布里渊区上的积分。

3. 生成函数方法

定义Bloch生成函数： G(z,k)=1NuTr(11−zA(k))G(z,k) = \frac{1}{N_u}\text{Tr}\left(\frac{1}{1-zA(k)}\right)

则有： G(z)=∫kG(z,k)=∑n≥0SnznG(z) = \int_k G(z,k) = \sum_{n \geq 0} S_n z^n

技术创新点

  1. 理论突破: 首次系统性地建立了图论中的回归走步问题与凝聚态物理中Bloch能带理论的深层联系
  2. 计算方法创新:
    • 发展了纯代数技术评估动量积分，使用常数部分提取方法
    • 对于复杂晶格，通过计算矩阵逆的对角元避免了特征值的显式计算
  3. 构造方法: 提出了通用的主公式(56)来构造任意大小的有限图： Ap,q=(1p⊗1q)⊗Au+∑d[(Rp(d)⊗Rq(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.}]

实验设置

研究对象

十一种阿基米德晶格：

  • 双分图: 正方(Square)、蜂窝(Honeycomb)、CaVO、SHD
  • 非双分图: 三角(Triangular)、Kagome、Trellis、Star、SrCuBO、Ruby、Maple-Leaf

验证方法

  1. 理论计算: 使用Bloch邻接矩阵和生成函数方法
  2. 数值验证: 构造大型有限图并计算邻接矩阵的幂
  3. 收敛性检验: 通过增加p,qp, q参数检验SnS_n的收敛性

计算工具

  • 使用Mathematica等计算代数软件计算矩阵逆
  • 通过Cayley-Hamilton定理从前几阶矩计算生成函数
  • 采用2500个均匀采样点进行布里渊区积分

实验结果

主要结果

1. 回归走步数量

成功计算了所有十一种晶格的回归走步序列，例如：

  • 正方晶格: S2n=(2nn)2S_{2n} = \binom{2n}{n}^2
  • 蜂窝晶格: S2n=∑l=0n(2ll)(nl)2S_{2n} = \sum_{l=0}^n \binom{2l}{l}\binom{n}{l}^2
  • 三角晶格: Sn=∑l=0n(nl)(−3)n−lS2l(H)S_n = \sum_{l=0}^n \binom{n}{l}(-3)^{n-l}S_{2l}^{(H)}

2. 态密度计算

获得了七种晶格的态密度解析表达式，包括：

  • 正方晶格: D□(E)=12π2K(1−E216)D_\square(E) = \frac{1}{2\pi^2}K(1-\frac{E^2}{16})
  • 蜂窝晶格: DH(E)=2∣E∣π2Z0(E)K(Z1(E)Z0(E))D_H(E) = \frac{2|E|}{\pi^2\sqrt{Z_0(E)}}K(\frac{Z_1(E)}{Z_0(E)})

3. 渐近行为

确定了回归概率的渐近系数α\alpha：

  • 双分图晶格: S2n∼αq2n2nS_{2n} \sim \alpha \frac{q^{2n}}{2n}
  • 非双分图晶格: Sn∼αqnnS_n \sim \alpha \frac{q^n}{n}

具体数值如正方晶格α=2π\alpha = \frac{2}{\pi}，三角晶格α=32π\alpha = \frac{\sqrt{3}}{2\pi}等。

验证结果

  1. 数值一致性: 理论计算结果与大型有限图的数值计算完全一致
  2. 已知结果重现: 成功重现了正方、蜂窝、三角、kagome晶格的已知序列
  3. 新序列发现: 为其余七种晶格首次计算了完整的回归走步序列

相关工作

传统方法

  1. 组合方法: 仅适用于简单晶格如正方晶格
  2. 邻接矩阵方法: 直接计算AnA^n，但对大系统计算困难
  3. 连分数方法: 通过Lanczos算法构造，但存在数值舍入误差

理论发展

  1. 谱矩定理: 本文结果可视为普遍谱矩定理的特殊情况
  2. 路径积分: 与离散路径积分理论相关
  3. 超几何函数: 与椭圆积分和超几何函数理论密切相关

本文优势

  1. 统一框架: 提供了处理所有阿基米德晶格的统一方法
  2. 解析结果: 获得了精确的解析表达式而非数值近似
  3. 物理洞察: 揭示了图论与能带理论的深层联系

结论与讨论

主要结论

  1. 理论贡献: 建立了回归走步问题与Bloch能带理论的根本联系
  2. 计算成果: 完整计算了十一种阿基米德晶格的回归走步数量
  3. 应用价值: 为态密度计算和渐近分析提供了新工具

局限性

  1. 适用范围: 主要针对symmorphic晶体学铺砌，对非周期结构不适用
  2. 计算复杂性: 对于大单元胞(如SHD的12个原子)，解析计算仍然复杂
  3. 维度限制: 虽然理论上可推广到高维，但具体计算复杂度会显著增加

未来方向

  1. 推广应用: 扩展到Laves晶格和其他周期性铺砌
  2. 非欧几何: 应用到双曲晶格等非欧几里得铺砌
  3. 相关问题: 扩展到自回避走步、带面积约束的走步等相关组合问题
  4. 物理应用: 应用到二维Ising模型的临界温度和自由能计算

深度评价

优点

  1. 理论创新性: 首次系统性建立图论与能带理论的联系，具有重要理论价值
  2. 计算完整性: 覆盖了所有阿基米德晶格，提供了完整的数据集
  3. 方法普适性: 发展的技术可推广到其他周期性铺砌
  4. 物理意义: 结果在态密度计算等凝聚态物理问题中有直接应用
  5. 数学严谨性: 推导严格，包含详细的数学证明

不足

  1. 新颖性限制: 核心恒等式(68)在某种程度上可视为已知谱矩定理的应用
  2. 计算效率: 对于某些复杂晶格，计算仍需依赖符号计算软件
  3. 实验验证: 缺乏与实际物理实验的对比验证
  4. 应用深度: 虽然提供了工具，但在具体物理问题中的深入应用有限

影响力

  1. 学术价值: 为图论、数学物理和凝聚态物理的交叉研究提供了新范例
  2. 实用价值: 为材料科学中的态密度计算提供了新工具
  3. 可复现性: 提供了详细的计算方法和结果，便于验证和应用
  4. 启发意义: 为相关组合问题的研究开辟了新思路

适用场景

  1. 理论研究: 晶格模型、统计物理、图论研究
  2. 材料计算: 新型二维材料的电子结构计算
  3. 算法开发: 大型稀疏矩阵特征值问题的新算法
  4. 教学应用: 作为数学物理交叉领域的典型案例

参考文献

论文引用了75篇相关文献，涵盖了从经典的Pólya随机走步理论到最新的双曲晶格研究，体现了该领域的历史发展和当前前沿。重要参考文献包括Wallace的石墨能带理论、Kitaev的量子自旋液体模型、以及近期的双曲晶格实验工作等。