Rapidly convergent quantum Monte Carlo using a Chebyshev projector

Zijun Zhao,Maria-Andreea Filip,Alex J W Thom
DOI: https://doi.org/10.1039/D4FD00035H
2024-05-02
Abstract:The multi-reference 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 multi-reference 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 at 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.
Strongly Correlated Electrons,Chemical Physics
What problem does this paper attempt to address?
The main problem this paper attempts to address is improving the convergence speed and computational efficiency of the Multi-Reference Coupled Cluster Monte Carlo (MR-CCMC) algorithm, making it more effective in handling large systems with significant strong correlation. Specifically, the authors introduce the following improvements: 1. **Efficient Search Algorithm**: - Using a logarithmic time complexity excitation acceptance algorithm based on BK trees, replacing the original linear scan algorithm, thereby significantly speeding up the process of finding determinants connected to the reference state in the multi-reference space. - Introducing a symmetric screening procedure to reduce the number of determinants in the active reference space with almost no loss of accuracy. 2. **Chebyshev Projector**: - Proposing an approximate wall projector based on Chebyshev polynomials to accelerate any projector-based quantum Monte Carlo algorithm. - This projector is implemented by truncating the expansion of Chebyshev polynomials, significantly reducing the number of Hamiltonian applications required to achieve the same statistical convergence, thereby lowering computational costs. 3. **Reference Space Compression**: - Developing a method for compressing the reference space, generating a compressed set of reference determinants that captures the most important part of the determinants in the CAS, thus significantly reducing computational costs while maintaining high accuracy. Through these improvements, the authors demonstrate that the new method can significantly enhance computational efficiency and convergence speed when dealing with systems such as carbon dimers and beryllium dimers, particularly excelling in handling challenging multi-reference problems. These improvements make the MR-CCMC algorithm more practical and closer to algorithmic maturity.