Parallel Multi-Coordinate Descent Methods for Full Configuration Interaction

Yuejia Zhang,Weiguo Gao,Yingzhou Li

2024-11-12

Abstract:We develop a multi-threaded parallel coordinate descent full configuration interaction algorithm (mCDFCI), for the electronic structure ground-state calculation in the configuration interaction framework. The FCI problem is reformulated as an unconstrained minimization problem, and tackled by a modified block coordinate descent method with a deterministic compression strategy. mCDFCI is designed to prioritize determinants based on their importance, with block updates enabling efficient parallelization on shared-memory, multi-core computing infrastructure. We demonstrate the efficiency of the algorithm by computing an accurate benchmark energy for the chromium dimer in the Ahlrichs SV basis (48e, 42o), which explicitly includes $2.07 \times 10^9$ variational determinants. We also provide the binding curve of the nitrogen dimer under the cc-pVQZ basis set (14e, 110o). Benchmarks show up to $79.3\%$ parallel efficiency on 128 cores.

Chemical Physics,Computational Physics

What problem does this paper attempt to address?

### What problem does this paper attempt to solve?

This paper aims to solve the **computational efficiency and parallelization challenges faced by the Full Configuration Interaction (FCI) method in handling large - scale electronic structure calculations**. Specifically, the FCI method can provide numerically exact solutions, but in practical applications, its computational complexity grows exponentially with the number of orbitals and electrons, making it difficult to apply to large systems.

#### Main problems:
1. **High computational complexity**: When dealing with systems with a large number of orbitals and electrons, the traditional FCI method has a prohibitively high computational cost and is difficult to implement.
2. **Difficult parallelization**: Although the existing Coordinate Descent FCI (CDFCI) is effective, its inherent serial nature limits its parallelization potential, thus affecting the efficiency of large - scale calculations.

#### Solutions:
To overcome the above problems, the authors propose the **multi - threaded parallel coordinate descent full configuration interaction algorithm (mCDFCI)**. This method improves the traditional CDFCI method in the following ways:

- **Multi - coordinate update**: Update multiple coordinates simultaneously in each iteration instead of a single coordinate, thereby accelerating convergence.
- **Deterministic compression strategy**: Adopt a deterministic compression strategy, giving priority to important Slater determinants to reduce the computational burden.
- **Efficient parallelization**: Utilize shared - memory, multi - core computing infrastructure to achieve efficient parallelization, significantly improving computational efficiency.

#### Application examples:
The authors demonstrate the efficiency and accuracy of the mCDFCI algorithm by calculating the benchmark energies of the chromium dimer (Cr₂) and the nitrogen molecule (N₂). For example, when calculating the benchmark energy of Cr₂ in the Ahlrichs SV basis set, 2.07×10⁹ variational determinants are explicitly included, and a parallel efficiency of up to 79.3% is achieved.

### Summary:
The main objective of this paper is to develop an efficient parallelization algorithm to address the computational complexity and parallelization problems of the FCI method in handling large - scale electronic structure calculations. By introducing the multi - coordinate descent method and the deterministic compression strategy, the mCDFCI method not only improves computational efficiency but also significantly enhances its applicability to large systems.