[论文解读] Hybrid Monte-Carlo simulation of interacting tight-binding model of graphene

Dominik Smith, Lorenz von Smekal|arXiv (Cornell University)|Nov 5, 2013

Graphene research and applications参考文献 12被引用 3

一句话总结

本研究采用哈伯德-斯特拉托诺维奇变换，对石墨烯的相互作用紧束缚模型进行混合蒙特卡洛模拟，以模拟长程库仑相互作用。研究识别出在有效精细结构常数 α_C ≈ 4 时发生导体-绝缘体相变，与先前工作一致，并表明尽管存在高阶校正，紧致哈伯德场与二阶离散化方案仍能产生相同的连续极限。

ABSTRACT

In this work, results are presented of Hybrid-Monte-Carlo simulations of the tight-binding Hamiltonian of graphene, coupled to an instantaneous long-range two-body potential which is modeled by a Hubbard-Stratonovich auxiliary field. We present an investigation of the spontaneous breaking of the sublattice symmetry, which corresponds to a phase transition from a conducting to an insulating phase and which occurs when the effective fine-structure constant $α$ of the system crosses above a certain threshold $α_C$. Qualitative comparisons to earlier works on the subject (which used larger system sizes and higher statistics) are made and it is established that $α_C$ is of a plausible magnitude in our simulations. Also, we discuss differences between simulations using compact and non-compact variants of the Hubbard field and present a quantitative comparison of distinct discretization schemes of the Euclidean time-like dimension in the Fermion operator.

研究动机与目标

• 通过非微扰格点模拟研究石墨烯相互作用紧束缚模型中的自发亚晶格对称性自发破缺。
• 确定有效精细结构常数 α_C 的临界值，该值对应于导体-绝缘体相变的发生。
• 评估不同离散化方案以及哈伯德场变体（紧致与非紧致）对模拟精度与收敛性的影响。
• 在一致的势参数与边界条件下，将模拟框架与先前结果进行对比验证，并为未来采用改进边界条件的大规模研究做好准备。

提出的方法

• 推导紧束缚哈密顿量的路径积分表述，通过哈伯德-斯特拉托诺维奇辅助场引入长程两体相互作用势。
• 利用 cRPA 值构建非局部、分段定义的势能，涵盖局域、最近邻、次近邻及六边形交叉项，并保留未屏蔽的库仑尾部。
• 在 6×6 矩形石墨烯晶格上，采用 GPU 加速的混合蒙特卡洛方法与欧几里得时间离散化进行模拟。
• 通过测量多个配置下的序参量 ⟨Δ_N⟩ 来探测自发对称性破缺，并外推至手征极限（m→0）。
• 通过将 ⟨Δ_N⟩ 拟合至 1/N_t 与 1/N_t² 项，量化离散化误差，分别针对标准与改进的二阶费米子作用量。
• 由于收敛性更优且能检测到非零序参量，优先采用紧致哈伯德场而非非紧致哈伯德场。

实验结果

研究问题

• RQ1 石墨烯的相互作用紧束缚模型是否在物理上合理的 α_C 值下表现出从导体到绝缘体的相变？
• RQ2 紧致与非紧致哈伯德场的表述方式如何影响亚晶格对称性破缺的模拟？
• RQ3 一阶与二阶离散化方案对费米子作用量连续极限的定量影响为何？
• RQ4 在使用一致势参数与边界条件时，该模拟框架能否重现先前关于 α_C 的结果？

主要发现

• 模拟识别出在 α_C ≈ 4 时发生导体-绝缘体相变，与先前研究报道的 α_C ≈ 3.12 在定性上一致。
• 非紧致哈伯德场无法产生非零序参量，表明存在数值不稳定性与收敛性差。
• 一阶与二阶离散化方案均收敛至相同的连续极限，拟合系数 c1 与 c2 的差异可忽略不计。
• 二阶方案在本设置中未表现出相对于标准方案的可测量改进，表明高阶项在此情形下无实际益处。
• 结果验证了该模拟框架作为未来大规模研究可靠平台的可行性，尤其适用于改进边界处理与势能匹配。
• 采用基于 FFT 的势能计算方法，使未来在更大体系规模下进行模拟并保持一致边界条件成为可能。