Abstract:Electronic-structure calculations play a fundamental role in predicting important physical (optical, electrical, etc) properties of condensed matter. Based on density functional theory (DFT), Kohn-Sham (KS) equation replaces the interacting many-body (electrons) problem by an equivalent set of self-consistent single-particle equations. Different from well-adopted basis set approach, we present a real-space finite-difference method to discretize the KS equation where the Laplacian operator is represented by high-order differences. After spatial discretization, a derived nonlinear eigenvalue problem is solved by a self-consistent field (SCF) iteration scheme with the update of electron density and potentials. In spite of the simplicity and versatility of the finite-difference method, one has to perform large-scale calculations of both eigenvalue and electrostatic problems at each iteration step. For eigenvalue problems, we minimize the Rayleigh quotients by the conjugate gradient method; and find minimum eigenvalues of interest through subspace diagonalization without Gram-Schmidt procedure. For obtaining electrostatic potentials, we use Dirichlet boundary condition to truncate the computational domain after inserting electron density and Gaussian compensating charge density together into Poisson's equation. Then a preconditioned Krylov subspace solver is employed to solve the Poisson's equation. A three-dimensional quantum dot with few spin-polarized electrons (each electron has spin up) is simulated by the real-space DFT incorporating local spin density (LSD) approximation. Wave function, electron density, eigen-energy, and total energy are calculated and compared to published results. The real-space finite-difference method produces reliable simulation results with a high computational efficiency. This work is fundamentally important to quantum-mechanical ab-initio calculation.

The power of Hellmann-Feynman theorem: Kohn-Sham DFT energy derivatives with respect to the parameters of exchange-correlation functional at linear cost

Differentiable programming for Kohn-Sham DFT: energy derivatives with respect to the parameters of exchange-correlation functional at linear cost

Self-Consistent Solution Of Kohn-Sham Equation By Real-Space Finite-Difference Method

D4FT: A Deep Learning Approach to Kohn-Sham Density Functional Theory

By-passing the Kohn-Sham equations with machine learning

Efficient Calculation of Derivatives of Integrals in a Basis of Non-Separable Gaussians Through Exploitation of Sparsity

Accurate and Efficient Calculations of Hellmann-Feynman Forces for Quantum Computation.

Efficient HF exchange evaluation through Fourier convolution in Cartesian grid for orbital-dependent density functionals

Recent Progress in Density Functional Theory and Its Numerical Methods

Inverse molecular design and parameter optimization with Hückel theory using automatic differentiation

Linear Jacobi-Legendre expansion of the charge density for machine learning-accelerated electronic structure calculations

Towards Efficient and Unified Treatment of Static and Dynamic Correlations in Generalized Kohn-Sham Density Functional Theory

Even Faster Exact Exchange for Solids via Tensor Hypercontraction

Linear-Scaling Subspace-Iteration Algorithm with Optimally Localized Nonorthogonal Wave Functions for Kohn-Sham Density Functional Theory

Consistent analytical second derivatives of the Kohn–Sham DFT energy in the framework of the Conductor-like screening Model (COSMO) through Gaussian charge distributions

Efficient real space formalism for hybrid density functionals

Achieving Performance Portability in Gaussian Basis Set Density Functional Theory on Accelerator Based Architectures in NWChemEx

Analytic derivative couplings in time-dependent density functional theory: Quadratic response theory versus pseudo-wavefunction approach.

Accelerating finite-temperature Kohn-Sham density functional theory with deep neural networks

Approaching the Basis Set Limit for DFT Calculations Using an Environment-Adapted Minimal Basis with Perturbation Theory: Formulation, Proof of Concept, and a Pilot Implementation

Exploiting the hessian for a better convergence of the SCF RDMFT procedure