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.

A first encounter with the Hartree-Fock self-consistent-field method

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

Convergence criteria for the Hartree iterative solution of the general self-consistent field equations

A iterative finite method approach to Kohn-Sham equations of light atoms: Density Functional Theory tutorial for undergraduate students

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

Extended Hartree-Fock method based on pair density functional theory

Efficient Hartree-Fock Exchange Algorithm with Coulomb Range Separation and Long-Range Density Fitting

Efficient Four-Component Dirac-Coulomb-Gaunt Hartree--Fock in Pauli Spinor Representation

Efficient and Robust Ab Initio Self-Consistent Field Acceleration Algorithm Based on a Semiempirical Model Hamiltonian

Resolution-of-identity approach to Hartree–Fock, hybrid density functionals, RPA, MP2 and GW with numeric atom-centered orbital basis functions

Solving the relativistic Hartree-Bogoliubov equation with the finite-difference method

An assessment of initial guesses for self-consistent field calculations. Superposition of Atomic Potentials: simple yet efficient

Self-consistent Density Functional Calculations without Diagonalization

Efficient Four-Component Dirac-Coulomb-Gaunt Hartree-Fock in the Pauli Spinor Representation.

Hartree-Fock Corrections in a Mean-field Limit for Fermions

Assessment of Basis Set Projection and Many-Body Expansion for Initial Guesses of Self-Consistent Field Calculations

Efficient periodic resolution-of-the-identity Hartree–Fock exchange method with k -point sampling and Gaussian basis sets

Numerically Stable Resonating Hartree-Fock

Robust formulation of Wick's theorem for computing matrix elements between Hartree-Fock-Bogoliubov wavefunctions

Implementation of time-dependent Hartree Fock in real space

Fortran 90 implementation of the Hartree-Fock approach within the CNDO/2 and INDO models