[论文解读] Tensor Product Variational Formulation for Quantum Systems

Yukinobu Nishio, Nobuya Maeshima+2 more|arXiv (Cornell University)|Jan 8, 2004

Quantum many-body systems被引用 46

一句话总结

本文提出了一种张量积变分方法（TPVA），用于近似二维量子自旋系统的基态，通过在二维晶格的单胞上使用局部权重的均匀乘积来表示试探波函数。通过将能量期望值映射为二维经典统计模型，该方法利用CTMRG实现精确的能量最小化，对于正方晶格上的S=1/2海森堡模型和XY模型，其变分能量与蒙特卡罗基准值的偏差在2.3%以内。

ABSTRACT

We consider a variational problem for the two-dimensional (2D) Heisenberg and XY models, using a trial state which is constructed as a 2D product of local weights. Variational energy is calculated by use of the the corner transfer matrix renormalization group (CTMRG) method, and its upper bound is surveyed. The variational approach is a way of applying the density matrix renormalization group method (DMRG) to infinite size 2D quantum systems.

研究动机与目标

• 开发一种适用于无限尺寸二维量子系统的变分方法，以应对标准DMRG在密度矩阵本征值衰减缓慢时失效的问题。
• 通过由局部权重构建的张量积态，将DMRG的变分结构推广至高维系统。
• 评估该张量积变分公式在正方晶格上S=1/2海森堡模型和XY模型中的数值效率与精度。
• 将变分能量结果与已建立的蒙特卡罗基准值及先前的变分方法进行比较。

提出的方法

• 将试探波函数表示为定义在二维晶格单胞上的局部权重W(σ₁,σ₂,σ₃,σ₄)的均匀乘积，其中包含三个独立参数a、b、c。
• 通过利用平移不变性与乘积结构，将能量期望值映射为单个键上的局部能量计算。
• 将内积⟨Ψ|Ψ⟩和能量分子表示为各向同性围面相互作用（IRF）模型的配分函数，从而可使用CTMRG进行数值计算。
• 利用CTMRG高精度计算配分函数，实现对参数a、b、c的系统性扫描，以最小化变分能量。
• 将该方法应用于S=1/2 XXZ模型（α=1为海森堡模型，α=0为XY模型），通过变换哈密顿量简化变分公式。
• 通过网格搜索优化参数a、b、c，以最小化局部能量期望值，从而获得变分基态估计。

实验结果

研究问题

• RQ1 张量积变分方法能否为如S=1/2海森堡模型和XY模型等二维量子自旋系统提供精确的基态能量估计？
• RQ2 与已建立的蒙特卡罗基准值及早期变分方法相比，TPVA的性能如何？
• RQ3 通过TPVA获得的变分态在海森堡模型中是否表现出长程序或呈现无序相？
• RQ4 通过增加变分参数数量或引入辅助变量，该方法是否可系统性地得到改进？

主要发现

• 对于S=1/2海森堡模型，TPVA得到的变分能量为-1.3089，与蒙特卡罗基准值-1.33888的偏差在2.3%以内。
• 对于XY模型，变分能量为-1.0848，与蒙特卡罗结果-1.09764相比偏差为1.2%。
• 海森堡模型的最优参数（a≈0.834，b≈0.748，c≈0.508）表明基态为无序态，未表现出反铁磁序的迹象。
• 非对角自旋关联SˣSˣ小于对角关联SᶻSᶻ，表明即使在海森堡点，变分态也表现出各向异性。
• 尽管能量分子结构复杂，该方法成功避免了符号问题，并通过CTMRG实现了高精度能量计算。
• 结果表明TPVA是无限二维量子系统的一种可行替代方法，但进一步改进可能需要引入更多变分参数或辅助自由度。