Rapidly convergent quantum Monte Carlo using a Chebyshev projector

Zijun Zhao,Maria-Andreea Filip,Alexander James William Thom
DOI: https://doi.org/10.1039/d4fd00035h
2024-05-03
Faraday Discussions
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.
chemistry, physical
What problem does this paper attempt to address?
The problems that this paper attempts to solve mainly focus on improving the efficiency and applicability of the multi - reference coupled - cluster Monte Carlo (MR - CCMC) algorithm. Specifically, through introducing a series of algorithm improvements, the author aims to achieve the following goals: 1. **Accelerate convergence**: By introducing a wall - Chebyshev projector based on Chebyshev polynomial expansion, the number of Hamiltonian applications required to reach statistical convergence is reduced, thereby reducing the computational cost. This new projector can be applied to MR - CCMC and full - configuration - interaction quantum Monte Carlo (FCIQMC) methods, significantly improving their computational efficiency. 2. **Optimize multi - reference processing**: In order to enable the MR - CCMC algorithm to handle multi - reference problems more efficiently, the author proposes an efficient cluster acceptance algorithm and reference space compression technique. These improvements include: - **Efficient cluster acceptance algorithm**: Use the BK - tree (Burkhard - Keller tree) data structure to accelerate the process of checking whether the generated spawnees are within a specific excitation level, thereby significantly reducing the computation time. - **Reference space compression**: By selectively retaining important determinants in the reference space, the size of the reference space is reduced while maintaining the accuracy of the calculation results. This method can significantly reduce the demand for computational resources without significant loss of precision. 3. **Enhance the robustness and stability of the algorithm**: In response to the population control problems that may occur after using the Chebyshev projector, the author proposes an improved shift update procedure to better control the growth of the particle number and ensure the stable operation of the algorithm. Through these improvements, the author hopes that the MR - CCMC algorithm can be more efficient and practical when dealing with strongly correlated systems, especially in the calculation of large - scale systems, and can provide more accurate results while maintaining a low computational cost. These improvements not only improve the performance of the algorithm but also provide new tools and methods for future quantum - chemical calculations.