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.

Real-time adaptive finite element solution of time-dependent Kohn–Sham equation

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

Numerical Solution of the Kohn-Sham Equation by Finite Element Methods with an Adaptive Mesh Redistribution Technique

Finite Element Method for Solving Kohn-Sham Equations Based on Self-Adaptive Tetrahedral Mesh

A Nonnested Augmented Subspace Method for Kohn-Sham Equation

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

Adaptive Finite Element Approximations for Kohn-Sham Models

An Arbitrary Lagrangian–Eulerian Method Based on the Adaptive Riemann Solvers for General Equations of State

On accelerating a multilevel correction adaptive finite element method for Kohn-Sham equation

An Adaptive High-Order Transient Algorithm to Solve Large-Scale Anisotropic Maxwell’s Equations

An adaptive multiresolution ultra-weak discontinuous Galerkin method for nonlinear Schrodinger equations

An Adaptive Time-Stepping Procedure Based On The Scaled Boundary Finite Element Method For Elastodynamics

An Asymptotic Green's Function Method for Time-Dependent Schrödinger Equations with Application to Kohn-Sham Equations

AN EFFICIENT H-ADAPTIVE SCALED BOUNDARY ELEMENT METHOD FOR TRANSIENT ELASTODYNAMICS

A high-order accurate moving mesh finite element method for the radial Kohn--Sham equation

Adaptive residual refinement in an RBF finite difference scheme for 2D time-dependent problems

Towards Translational Invariance of Total Energy with Finite Element Methods for Kohn-Sham Equation

Adaptive finite element approximation for steady-state Poisson-Nernst-Planck equations

Robust a Posteriori Error Estimates of Time-Dependent Poisson–Nernst–Planck Equations

A fast energy conserving finite element method for the nonlinear fractional Schrödinger equation with wave operator

3D Adaptive Finite Element Method for a Phase Field Model for the Moving Contact Line Problems