Zhao Zijun, Filip Maria-Andreea, Thom Alex J W
Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, UK.
Faraday Discuss. 2024 Nov 6;254(0):429-450. doi: 10.1039/d4fd00035h.
The multireference coupled-cluster Monte Carlo (MR-CCMC) algorithm is a determinant-based quantum Monte Carlo (QMC) algorithm that is conceptually similar to Full Configuration Interaction QMC (FCIQMC). It has been shown to offer a balanced treatment of both static and dynamic correlation while retaining polynomial scaling, although application to large systems with significant strong correlation remained impractical. In this paper, we document recent algorithmic advances that enable rapid convergence and a more black-box approach to the multireference problem. These include a logarithmically scaling metric-tree-based excitation acceptance algorithm to search for determinants connected to the reference space at the desired excitation level and a symmetry-screening procedure for the reference space. We show that, for moderately sized reference spaces, the new search algorithm brings about an approximately 8-fold acceleration of one MR-CCMC iteration, while the symmetry screening procedure reduces the number of active reference space determinants with essentially no loss of accuracy. We also introduce a stochastic implementation of an approximate wall projector, which is the infinite imaginary time limit of the exponential projector, using a truncated expansion of the wall function in Chebyshev polynomials. Notably, this wall-Chebyshev projector can be used to accelerate any projector-based QMC algorithm. We show that it requires significantly fewer applications of the Hamiltonian to achieve the same statistical convergence. We benchmark these acceleration methods on the beryllium and carbon dimers, using initiator FCIQMC and MR-CCMC with basis sets up to cc-pVQZ quality.
多参考耦合簇蒙特卡罗(MR - CCMC)算法是一种基于行列式的量子蒙特卡罗(QMC)算法,在概念上与全组态相互作用QMC(FCIQMC)相似。尽管将其应用于具有显著强关联的大型系统仍然不切实际,但已证明它能在保持多项式缩放比例的同时,对静态和动态关联进行平衡处理。在本文中,我们记录了近期的算法进展,这些进展能够实现快速收敛,并为多参考问题提供更黑箱式的方法。这些进展包括一种基于对数缩放度量树的激发接受算法,用于在期望的激发水平下搜索与参考空间相连的行列式,以及一种针对参考空间的对称筛选程序。我们表明,对于中等规模的参考空间,新的搜索算法使一次MR - CCMC迭代加速约8倍,而对称筛选程序减少了活跃参考空间行列式的数量,且基本不损失精度。我们还引入了近似壁投影算符的随机实现,它是指数投影算符的无限虚时极限,使用切比雪夫多项式对壁函数进行截断展开。值得注意的是,这种壁 - 切比雪夫投影算符可用于加速任何基于投影算符的QMC算法。我们表明,实现相同的统计收敛所需的哈密顿量应用次数显著减少。我们使用初始FCIQMC和MR - CCMC,并采用高达cc - pVQZ质量的基组,在铍和碳二聚体上对这些加速方法进行了基准测试。