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.

WKB Methods for Finite Difference Schrodinger Equations

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

High-order WKB-based Method For The 1D Stationary Schrödinger Equation In The Semi-classical Limit

WKB-based third order method for the highly oscillatory 1D stationary Schrödinger equation

Mathematical and Computational Methods for Semiclassical Schrödinger Equations

Multi-phase Computations of the Semiclassical Limit of the Schrödinger Equation and Related Problems: Whitham Vs Wigner

Study on a Quantization Condition and the Solvability of Schrödinger-type Equations

Numerical Methods for the Wigner Equation with Unbounded Potential

A High Order Efficient Numerical Method for 4-D Wigner Equation of Quantum Double-Slit Interferences

Explicit summation of the constituent WKB series and new approximate wave functions

On nonlinear effects in multiphase WKB analysis for the nonlinear Schrodinger equation

Structured Backward Error for the WKB method

An energy-based summation-by-parts finite difference method for the wave equation in second order form

On nonlinear effects in multiphase WKB analysis for the nonlinear Schrödinger equation *

Modified Semi-Classical Methods for Nonlinear Quantum Oscillations Problems

Optimally truncated WKB approximation for the highly oscillatory stationary 1D Schrödinger equation

Stationary Schrödinger equation in the semi-classical limit: numerical coupling of oscillatory and evanescent regions

Optimally truncated WKB approximation for the 1D stationary Schrödinger equation in the highly oscillatory regime

On numerical solutions of the time-dependent Schrödinger equation

PERTURBATION-WKB APPROXIMATION : QUARTIC ANHARMONIC OSCILLATOR

Exact WKB in all sectors I: Potentials with degenerate saddles