[论文解读] A Renormalisation-Group Algorithm for Eigenvalue Density Functions of Interacting Quantum Systems

Tobias J. Osborne|arXiv (Cornell University)|May 8, 2006

Quantum many-body systems被引用 3

一句话总结

该论文提出了一种基于矩阵乘积算符（MPO）的可验证重整化群算法，以多项式资源缩放的方式计算一维相互作用量子自旋系统中 smeared 特征值密度函数。该方法提供了数值稳定、可系统性改进的态密度近似，包括在指定置信区间内对基态能隙的精确估计，适用于无序系统和临界系统。

ABSTRACT

We present a certifiable algorithm to calculate the eigenvalue density function -- the number of eigenvalues within an infinitesimal interval -- for an arbitrary 1D interacting quantum spin system. Our method provides an arbitrarily accurate numerical representation for the smeared eigenvalue density function, which is the convolution of the eigenvalue density function with a gaussian of prespecified width. In addition, with our algorithm it is possible to investigate the density of states near the ground state. This can be used to numerically determine the size of the ground-state energy gap for the system to within a prespecified confidence interval. Our method exploits a finitely correlated state/matrix product state representation of the propagator and applies equally to disordered and critical interacting 1D quantum spin systems. We illustrate our method by calculating an approximation to the eigenvalue density function for a random antiferromagnetic Heisenberg model.

研究动机与目标

• 开发一种用于近似一维相互作用量子自旋系统特征值密度函数的数值可验证方法。
• 克服标准 DMRG 及相关方法中缺乏验证的问题，这些方法可能收敛至局部极小值。
• 通过分析基态附近的态密度，实现对基态能隙的精确估计。
• 将方法推广至非平移不变及非临界系统，包括无序和临界量子自旋链。
• 通过基于 MPO 的传播子表示，提供对 smeared 特征值密度函数的多项式时间、系统性近似。

提出的方法

• 该方法使用矩阵乘积算符（MPO）表示时间演化算符 $ U(t) = e^{-iHt} $，利用 Lie-Trotter 分裂将哈密顿量分解为偶 sites 和奇 sites 之间的相互作用。
• 通过迭代组合短时间演化算符 $ W = e^{iA\frac{t}{m}}e^{iB\frac{t}{m}} $ 构建 MPO 表示，并通过希尔伯特-施密特范数最小化实现截断，以控制辅助维度。
• 利用迹规则 $ \operatorname{tr}(\sigma^\alpha) = \delta_{\alpha,0} $ 高效计算 MPO 近似 $ \widetilde{V}(t) $ 的迹，得到离散时间序列 $ g_k \sim \operatorname{tr}(W^k) $。
• 对时间序列应用高斯窗 $ e^{-k^2/(2s^2)} $ 以抑制高频噪声，生成对 smeared 特征值密度的平滑近似。
• 对加窗后的时间序列应用离散傅里叶变换，得到 smeared 密度 $ \widetilde{\mu}_H(x) $ 的频域表示，采样受奈奎斯特准则控制。
• 该方法确保最终输出中所有特征值均被表示，通过 smearing 以可控不确定性实现对精确特征值定位的权衡。

实验结果

研究问题

• RQ1 能否开发一种可验证算法，以受控误差界计算一维相互作用量子自旋系统中特征值密度函数？
• RQ2 如何利用系统性近似方法，在指定置信区间内对基态能隙进行数值估计？
• RQ3 该方法是否可应用于无序和临界量子自旋链，而无需依赖平移不变性或非临界性？
• RQ4 该算法在系统尺寸和期望精度方面的缩放行为如何？
• RQ5 基于 MPO 的截断策略如何在近似时间演化算符的计算成本与保真度之间实现平衡？

主要发现

• 该算法提供了一种多项式可扩展的方法来计算 smeared 特征值密度函数 $ \widetilde{\mu}_H(x) $，误差可通过高斯 smearing 宽度 $ s $ 控制。
• 该方法保证输出中包含所有特征值，避免了标准 DMRG 等非验证方法中可能遗漏态的风险。
• 基态能隙可在一个预设置信区间内进行估计，从而实现对量子自旋链中能隙的可靠数值分析。
• 对于 20 个自旋的随机反铁磁海森堡模型，该方法成功以高分辨率捕捉了态密度，证明了其实际适用性。
• MPO 截断过程确保辅助维度仅随系统尺寸多项式增长，使该方法可扩展至大规模系统。
• 通过改进的 MPO 技术，该方法可推广至高维系统，从而将适用范围扩展至一维系统之外。