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.

Adaptive coordinate, real-space electronic structure calculations on parallel computers

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

Adaptive Finite Differencing in High Accuracy Electronic Structure Calculations

Parallel implementation of electronic structure eigensolver using a partitioned folded spectrum method

Parallel Mesh Refinement Of Higher Order Finite Elements For Electronic Structure Calculations

Adaptive Algorithm for Electronic Structure Calculations Using Reduction of Gaussian Mixtures

An auxilliary grid method for the calculation of electrostatic terms in density functional theory on a real-space grid.

An H-Adaptive Finite Element Solver for the Calculations of the Electronic Structures

Dgdft: A Massively Parallel Method for Large Scale Density Functional Theory Calculations

Two-Level Chebyshev Filter Based Complementary Subspace Method: Pushing the Envelope of Large-Scale Electronic Structure Calculations.

Low-rank approximations for accelerating plane-wave hybrid functional calculations in unrestricted and noncollinear spin density functional theory

Application of complex energy integration to selfconsistent electronic structure calculations

O(N) methods in electronic structure calculations

Auxiliary basis expansions for large-scale electronic structure calculations

A parallel orbital-updating based plane-wave basis method for electronic structure calculations

Adaptive Regularized Self-Consistent Field Iteration with Exact Hessian for Electronic Structure Calculation.

Efficient real space formalism for hybrid density functionals

Linear scaling Krylov subspace method for large scale {\it ab initio} electronic structure calculations of metals

2.5 Million-Atom Ab Initio Electronic-Structure Simulation of Complex Metallic Heterostructures with DGDFT

Solving the electronic structure problem for over 100,000 atoms in real-space