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.

KS-pies: Kohn-Sham inversion toolkit

KS-pies: Kohn-Sham Inversion Toolkit

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

Kohn-Sham Potentials From Electron Densities Using A Matrix Representation Within Finite Atomic Orbital Basis Sets

Inverse Kohn-Sham Density Functional Theory: Progress and Challenges

Kohn-Sham inversion with mathematical guarantees

A Story of Three Levels of Sophistication in SCF/KS-DFT Orbital Optimization Procedures

On the Analysis of the Discretized Kohn-Sham Density Functional Theory

Inverting the Kohn-Sham equations with physics-informed machine learning

KSSOLV—a MATLAB Toolbox for Solving the Kohn-Sham Equations

KSSOLV 2.0: an Efficient MATLAB Toolbox for Solving the Kohn-Sham Equations with Plane-Wave Basis Set.

Robust Determination Of The Chemical Potential In The Pole Expansion And Selected Inversion Method For Solving Kohn-Sham Density Functional Theory

Measuring density-driven errors using Kohn-Sham inversion

A Critical Evaluation of the Hybrid KS DFT Functionals Based on the KS Exchange-Correlation Potentials

Accelerating inverse Kohn-Sham calculations using reduced density matrices

Accelerating Atomic Orbital-based Electronic Structure Calculation via Pole Expansion and Selected Inversion

Fractionalized Kohn-Sham ansatz for strongly-correlated electrons

A Simple and Accurate Method for Computing Optimized Effective Potentials for Exact Exchange Energy

On the Convergence of the Self-Consistent Field Iteration in Kohn--Sham Density Functional Theory

Computing the Self-Consistent Field in Kohn-Sham Density Functional Theory