Random Green's function method for large-scale electronic structure calculation

Mingfa Tang,Chang Liu,Aixia Zhang,Qingyun Zhang,Shengjun Yuan,Youqi Ke

2024-03-03

Abstract:We report a linear-scaling random Green's function (rGF) method for large-scale electronic structure calculation. In this method, the rGF is defined on a set of random states to stochastically express the density matrix, and rGF is calculated with the linear-scaling computational cost. We show the rGF method is generally applicable to the nonorthogonal localized basis, and circumvent the large Chebyshev expansion for the density matrix. As a demonstration, we implement rGF with density-functional Tight-Binding method and apply it to self-consistently calculate water clusters up 9984 H2Os. We find the rGF method combining with a simple fragment correction can reach an error of ~1meV per H2O in total energy, compared to the deterministic calculations, due to the self-average. The development of rGF method advances the stochastic electronic structure theory to a new stage of the efficiency and applicability.

Mesoscale and Nanoscale Physics,Materials Science,Computational Physics

What problem does this paper attempt to address?

This paper attempts to solve the problem of computational complexity in large - scale electronic structure calculations, especially for systems containing hundreds or thousands of atoms. The traditional Kohn - Sham density functional theory (KSDFT) faces the problem of cubic - level computational complexity when dealing with these large - scale systems, resulting in huge consumption of computational resources and low efficiency.

To solve this problem, the author proposes a linear - scaling stochastic Green's function (rGF) method based on random states. This method effectively calculates the Fermi - Dirac operator by projecting random states onto the Krylov subspace. Compared with the traditional polynomial expansion method, the rGF method can significantly reduce the amount of calculation while maintaining high precision, and is suitable for narrow - band - gap materials at low temperatures.

Specifically, the main contributions of this method include:
1. **Linear - scaling calculation**: Through Krylov subspace projection, a linear - scaling computational complexity is achieved, thereby greatly improving the computational efficiency.
2. **Efficient handling of low - temperature systems**: For narrow - band - gap materials at low temperatures, the rGF method only requires a relatively small Krylov subspace dimension to achieve high precision, while traditional methods such as Chebyshev polynomial expansion require a very high order.
3. **Combination with density - functional tight - binding method (DFTB)**: The author combines the rGF method with DFTB and verifies its effectiveness and accuracy in practical applications, especially in self - consistent calculations on water molecule clusters and silicon hydride clusters.

In conclusion, this paper aims to develop an efficient, linear - scaling electronic structure calculation method to meet the computational challenges in large - scale material simulations.