[论文解读] Simulating the Generalized Gibbs Ensemble (GGE): a Hilbert space Monte Carlo approach

Vincenzo Alba|arXiv (Cornell University)|Jul 24, 2015

Quantum many-body systems被引用 11

一句话总结

该论文提出了一种希尔伯特空间蒙特卡罗方法，用于在可积量子自旋链中模拟截断广义吉布斯系综（TGGE），利用贝特 ansatz 状态和基于 GGE 概率测度的重要度采样。该方法可精确计算热力学极限可观测量，并直接提取贝特-高金-田中根分布，有限尺寸效应呈指数衰减，结果与热力学贝特 ansatz 和量子转移矩阵基准一致。

ABSTRACT

By combining classical Monte Carlo and Bethe ansatz techniques we devise a numerical method to construct the Truncated Generalized Gibbs Ensemble (TGGE) for the spin-1/2 isotropic Heisenberg ($XXX$) chain. The key idea is to sample the Hilbert space of the model with the appropriate GGE probability measure. The method can be extended to other integrable systems, such as the Lieb-Liniger model. We benchmark the approach focusing on GGE expectation values of several local observables. As finite-size effects decay exponentially with system size, moderately large chains are sufficient to extract thermodynamic quantities. The Monte Carlo results are in agreement with both the Thermodynamic Bethe Ansatz (TBA) and the Quantum Transfer Matrix approach (QTM). Remarkably, it is possible to extract in a simple way the steady-state Bethe-Gaudin-Takahashi (BGT) roots distributions, which encode complete information about the GGE expectation values in the thermodynamic limit. Finally, it is straightforward to simulate extensions of the GGE, in which, besides the local integral of motion (local charges), one includes arbitrary functions of the BGT roots. As an example, we include in the GGE the first non-trivial quasi-local integral of motion.

研究动机与目标

• 开发一种数值方法，用于在有限尺寸可积量子多体系统中模拟截断广义吉布斯系综（TGGE）。
• 实现在热力学极限下对局域可观测量的 GGE 期望值进行高精度计算。
• 通过蒙特卡罗采样直接提取稳态贝特-高金-田中（BGT）根分布，其编码了 GGE 的全部信息。
• 将该方法扩展至包含准局域荷和 BGT 根的任意函数的 GGE。
• 与已建立的方法（如热力学贝特 ansatz（TBA）和量子转移矩阵（QTM））进行基准测试。

提出的方法

• 该方法在希尔伯特空间中使用基于局域和准局域守恒荷构造的 GGE 概率测度进行重要度采样。
• 通过贝特-高金-田中（BGT）方程生成本征态，从而精确表示多体谱。
• 蒙特卡罗采样通过基于 GGE 哈密顿量指数的权重实现，拉格朗日乘子由初始态期望值固定。
• 代表性密度矩阵 ρ^sp 作为采样根分布的鞍点极限获得，收敛至真实的热力学分布。
• 应用有限尺寸标度以获得热力学极限，观察到有限尺寸修正呈指数衰减。
• 该方法被推广以包含 BGT 根的任意函数，如第一个非平凡的准局域荷。

实验结果

研究问题

• RQ1 能否在有限尺寸可积系统中，通过希尔伯特空间上的蒙特卡罗采样高效模拟截断广义吉布斯系综（TGGE）？
• RQ2 在 XXX 自旋链中，GGE 期望值的有限尺寸效应在多大程度上随系统尺寸衰减？
• RQ3 能否通过蒙特卡罗采样直接提取稳态贝特-高金-田中（BGT）根分布，且其是否收敛至热力学极限？
• RQ4 该方法在多大程度上能准确再现热力学贝特 ansatz（TBA）和量子转移矩阵（QTM）方法的已知结果？
• RQ5 该方法能否扩展至包含 GGE 中的准局域荷和 BGT 根的任意函数？

主要发现

• GGE 期望值中的有限尺寸修正随系统尺寸呈指数衰减，即使在中等大小的自旋链（如 L=30）中，也能实现对热力学极限的高精度外推。
• 局域可观测量（包括能量、自旋磁化率和更高阶荷）的蒙特卡罗结果与 TBA 和 QTM 基准结果高度一致。
• 该方法成功提取了稳态 BGT 根分布，数值直方图随系统尺寸增大而收敛至解析的 TBA 结果。
• 对于无限温度吉布斯系综，根分布呈现 ∼1/x⁴ 的尾部行为，与解析预测一致。
• 在 GGE 中引入第一个非平凡的准局域荷后，ρ₁(x) 分布在小 x 处表现出双峰结构，表明存在非平凡关联。
• 使用单个 CPU 核对 30 链 sites 系统的模拟耗时约两天，表明该方法在标准硬件上具有可行性。