A high-order accurate moving mesh finite element method for the radial Kohn--Sham equation

Zheming Luo,Yang Kuang
2024-11-07
Abstract:In this paper, we introduce a highly accurate and efficient numerical solver for the radial Kohn--Sham equation. The equation is discretized using a high-order finite element method, with its performance further improved by incorporating a parameter-free moving mesh technique. This approach greatly reduces the number of elements required to achieve the desired precision. In practice, the mesh redistribution involves no more than three steps, ensuring the algorithm remains computationally efficient. Remarkably, with a maximum of $13$ elements, we successfully reproduce the NIST database results for elements with atomic numbers ranging from $1$ to $92$.
Numerical Analysis,Computational Physics
What problem does this paper attempt to address?
### What problems does this paper attempt to solve? This paper aims to solve the problem of efficient and high - precision numerical solution of the radial Kohn - Sham equation. Specifically, the author proposes a new algorithm based on the moving - grid finite - element method to improve the computational efficiency and accuracy of solving the radial Kohn - Sham equation. The following are the specific problems that the paper attempts to solve: 1. **Numerical challenges**: - Kohn - Sham density functional theory (KSDFT) is an electron - structure method widely used in condensed - matter physics and quantum chemistry. However, there are numerical challenges in solving the Kohn - Sham equation. In particular, the singular potential near the atomic nucleus causes the wave function to change sharply in the core region and smoothly far from the atomic nucleus. - Traditional methods such as the uniform - grid method require a very high grid density in the core region to maintain accuracy, while the radial grid, although more flexible, requires the user to define multiple parameters, which increases the complexity of use. 2. **Improving computational efficiency and accuracy**: - The author introduces the high - order finite - element method, which has an exponential convergence characteristic with respect to the polynomial order and high grid flexibility. - To further improve efficiency, the author proposes a parameter - free moving - grid technique, which significantly reduces the required number of grids by adaptively concentrating grid points in areas requiring higher resolution. 3. **Verifying the effectiveness of the method**: - The author verifies the accuracy and efficiency of the proposed moving - grid finite - element method in handling elements from hydrogen (H) to uranium (U) by comparing with the results in the NIST database. The results show that the results in the NIST database can be reproduced using no more than 13 grid points and a 10 - order polynomial, demonstrating the high efficiency and accuracy of this method. ### Method overview - **High - order finite - element method**: The numerical solution accuracy is improved by approximating the wave function with high - order polynomials. - **Moving - grid technique**: Based on the equidistribution principle, the grid - point distribution is adaptively adjusted so that the error is evenly distributed in each sub - interval. - **Pre - conditioned LOBCPG method**: It is used to solve the generalized eigenvalue problem, significantly reducing the number of iterations and computational time. In conclusion, the method proposed in this paper not only solves the numerical challenges in traditional methods but also greatly improves computational efficiency and accuracy, providing a new and effective way to solve the radial Kohn - Sham equation.