[论文解读] Accurate variational electronic structure calculations with the density matrix renormalization group

Sebastian Wouters|arXiv (Cornell University)|May 6, 2014

Quantum many-body systems参考文献 2被引用 3

一句话总结

本文提出了一种基于密度矩阵重整化群（DMRG）与矩阵乘积态（MPS）波函数的精确变分电子结构计算方法。通过利用SU(2)自旋和U(1)粒子数等对称性，并引入非阿贝尔对称性适应，该方法实现了高精度与高效率，能够精确处理氢链和碳二聚体等体系中的强关联效应。

ABSTRACT

During the past 15 years, the density matrix renormalization group (DMRG) has become increasingly important for ab initio quantum chemistry. The underlying matrix product state (MPS) ansatz is a low-rank decomposition of the full configuration interaction tensor. The virtual dimension of the MPS controls the size of the corner of the many-body Hilbert space that can be reached. Whereas the MPS ansatz will only yield an efficient description for noncritical one-dimensional systems, it can still be used as a variational ansatz for other finite-size systems. Rather large virtual dimensions are then required. The two most important aspects to reduce the corresponding computational cost are a proper choice and ordering of the active space orbitals, and the exploitation of the symmetry group of the Hamiltonian. By taking care of both aspects, DMRG becomes an efficient replacement for exact diagonalization in quantum chemistry. DMRG and Hartree-Fock theory have an analogous structure. The former can be interpreted as a self-consistent mean-field theory in the DMRG lattice sites, and the latter in the particles. It is possible to build upon this analogy to introduce post-DMRG methods. Based on an approximate MPS, these methods provide improved ansätze for the ground state, as well as for excitations. Exponentiation of the single-particle (single-site) excitations for a Slater determinant (an MPS with open boundary conditions) leads to the Thouless theorem for Hartree-Fock theory (DMRG), an explicit nonredundant parameterization of the entire manifold of Slater determinants (MPS wavefunctions). This gives rise to the configuration interaction expansion for DMRG. The Hubbard-Stratonovich transformation lies at the basis of auxiliary field quantum Monte Carlo for Slater determinants. An analogous transformation for spin-lattice Hamiltonians allows to formulate a promising variant for MPSs.

研究动机与目标

• 开发一种可系统性改进的、变分的电子结构方法，用于强关联体系，基于DMRG算法。
• 通过将精确对称性（SU(2)、U(1)、阿贝尔点群）引入MPS试探波函数，提升计算效率与精度。
• 实现对具有强静态关联的多参考体系（如碳二聚体和氢链）的精确处理。
• 将DMRG扩展至基态能量之外，引入后-DMRG方法以描述激发态并改进试探波函数。
• 探索在一维分子体系中，热力学极限下均匀MPS的可行性。

提出的方法

• 采用矩阵乘积态（MPS）波函数作为全配置相互作用（CI）波函数的低秩张量分解，其中虚维D控制关联的捕获能力。
• 通过构建属于哈密顿量对称群本征态的波函数，实现对称性自适应MPS，利用Clebsch-Gordan系数与约化张量。
• 应用Wigner-Eckart定理对MPS进行因子分解，引入块稀疏性，降低内存与计算成本。
• 在CheMPS2 DMRG代码中实现对SU(2)自旋、U(1)粒子数及阿贝尔点群对称性的支持。
• 通过单体激发的指数化，引入Thouless型参数化，实现对MPS流形的非冗余参数化。
• 探索基于Hubbard-Stratonovich变换的辅助场量子蒙特卡罗方法在自旋-晶格哈密顿量上的MPS应用。

实验结果

研究问题

• RQ1 如何在MPS试探波函数中利用对称性，在保持精度的同时降低DMRG计算的计算成本？
• RQ2 在分子体系的QC-DMRG中，实现快速收敛的最优轨道排序与选择是什么？
• RQ3 基于组态相互作用展开的后-DMRG方法能否改善DMRG基态与激发态的描述？
• RQ4 随着虚维D增加，DMRG能量中的系统性偏差如何表现？是否可被可靠地外推？
• RQ5 在一维从头算体系中，具有指数衰减关联的热力学极限下，均匀MPS试探波函数是否可行？

主要发现

• CheMPS2实现成功解析了具有多参考特征的碳二聚体的低激发态，得益于对称性子空间的利用。
• 对于氢链体系，QC-DMRG在较小虚维D下即实现数值收敛，归因于有效的电子屏蔽，从而可精确计算纵向响应性质。
• 对称性的使用通过块稀疏性与信息压缩显著节省了内存与时间开销，尤其在非阿贝尔群中效果显著。
• 基于MPS行走者的辅助场量子蒙特卡罗方法在后-DMRG中表现出前景，可用于改进能量与激发能估计。
• 随着D增加，DMRG能量中的系统性偏差可被估计并可能实现外推，支持可靠的收敛监测。
• 对于聚烯烃等一维体系，由于关联呈指数衰减，均匀MPS试探波函数是可行的，可实现热力学极限下的研究。