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.

Linear Scaling Solution of the Time-Dependent Self-Consistent-Field Equations

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

Self-consistent field theory based molecular dynamics with linear system-size scaling

Sublinear scaling in non-Markovian open quantum systems simulations

Linear-Scaling Density Matrix Perturbation Treatment Of Electric Fields In Solids

Time-Dependent Quantum Transport: An Efficient Method Based On Liouville-Von-Neumann Equation For Single-Electron Density Matrix

Scale Transformation Method in Self-Consistent Solution of Schrödinger and Poisson Equations

Linear Scaling Approach for Optical Excitations Using Maximally Localized Wannier Functions

Virtual states in the coupled-channel problems with an improved complex scaling method

Linear scaling DFT calculations for large Tungsten systems using an optimized local basis

Linear Scaling Local Correlation Approach for Solving the Coupled Cluster Equations of Large Systems.

Stochastic Multiconfigurational Self-Consistent Field Theory.

Linear scaling calculation of band edge states and doped semiconductors

Linear scaling computation of the Fock matrix VII. Periodic Density Functional Theory at the $Γ$-point

A Sub-Linear Scaling Algorithm for Computing the Electronic Structure of Materials

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

Dynamics and Scaling of 2D Polymers in a Dilute Solution

Linear-Scaling Selected Inversion based on Hierarchical Interpolative Factorization for Self Green's Function for Modified Poisson-Boltzmann Equation in Two Dimensions

Density Matrix Based Time-Dependent Density Functional Theory and the Solution of Its Linear Response in Real Time Domain.