Abstract:Quantum mechanical calculations for material modelling using Kohn-Sham density functional theory (DFT) involve the solution of a nonlinear eigenvalue problem for $N$ smallest eigenvector-eigenvalue pairs with $N$ proportional to the number of electrons in the material system. These calculations are computationally demanding and have asymptotic cubic scaling complexity with the number of electrons. Large-scale matrix eigenvalue problems arising from the discretization of the Kohn-Sham DFT equations employing a systematically convergent basis traditionally rely on iterative orthogonal projection methods, which are shown to be computationally efficient and scalable on massively parallel computing architectures. However, as the size of the material system increases, these methods are known to incur dominant computational costs through the Rayleigh-Ritz projection step of the discretized Kohn-Sham Hamiltonian matrix and the subsequent subspace diagonalization of the projected matrix. This work explores the potential of polynomial expansion approaches based on recursive Fermi-operator expansion as an alternative to the subspace diagonalization of the projected Hamiltonian matrix to reduce the computational cost. Subsequently, we perform a detailed comparison of various recursive polynomial expansion approaches to the traditional approach of explicit diagonalization on both multi-node CPU and GPU architectures and assess their relative performance in terms of accuracy, computational efficiency, scaling behaviour and energy efficiency.

Calculation of the eigenfunctions and eigenvalues of Schrödinger type equations by asymptotic Taylor expansion method (ATEM)

Remarks on the Solution of the Position Dependent Mass (PDM) Schrödinger Equation

An Innovative Pseudo-Spectral Galerkin Algorithm for the Time-Fractional Tricomi-Type Equation

A new approach in handling one-dimensional time-fractional Schrödinger equations

Eigenvalue superposition for Toeplitz matrix-sequences with matrix order dependent symbols

A novel analytical iterative approach to time-fractional vibration equation using Aboodh transform

An Asymptotic Green's Function Method for Time-Dependent Schrödinger Equations with Application to Kohn-Sham Equations

Fast Non-Hermitian Toeplitz Eigenvalue Computations, Joining Matrixless Algorithms and FDE Approximation Matrices

Exact Calculation of the Mean-Square Error in the Method of Expansion of Iterated Ito Stochastic integrals Based on Generalized Multiple Fourier Series

Taylor series for the Askey-Wilson operator and classical summation formulas

ALGORITHMIC CALCULATION OF TWO-LOOP FEYNMAN DIAGRAMS

The Extended Weierstrass Transformation Method for the Biswas–Arshed Equation with Beta Time Derivative

The power-law and the logarithmic potentials

Stability analysis, lump and exact solutions to Sharma–Tasso–Olver–Burgers equation

Accurate analytical evaluation of the generalized logarithmic and double Fermi–Dirac and Bose–Einstein functions

On the convergence properties of generalized Szász–Kantorovich type operators involving Frobenious–Euler–Simsek-type polynomials

New optical soliton solutions of the popularized anti-cubic nonlinear Schrödinger equation versus its numerical treatment

Computational Solution of the Time-Fractional Schrödinger Equation by Using Trigonometric B-Spline Collocation Method

Subspace recursive Fermi-operator expansion strategies for large-scale DFT eigenvalue problems on HPC architectures

Numerical computations of coupled fractional resonant Schrödinger equations arising in quantum mechanics under conformable fractional derivative sense

Novel analytical STFT expressions for nonlinear power engineering problem solving