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.

Finite Volume Discretizations for Eigenvalue Problems with Applications to Electronic Structure Calculations.

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

Finite Element Approximations for Schrodinger Equations with Applications to Electronic Structure Computations

High-order Symplectic FDTD Scheme for Solving a Time-Dependent Schrödinger Equation

Local and Parallel Finite Element Discretizations for Eigenvalue Problems

Numerical Approximations of a Nonlinear Eigenvalue Problem and Applications to a Density Functional Model

Numerical Analysis of Finite Dimensional Approximations of Kohn-Sham Models

Convergence of Adaptive Finite Element Approximations for Nonlinear Eigenvalue Problems

Non-Disjoint Hexahedral Space Discretization For The Finite-Volume Technique

A Time-Dependent Random State Approach for Large-Scale Density Functional Calculations

Numerical Analysis and Computation of the Finite Volume Element Method for the Nonlinear Coupled Time-Fractional Schrödinger Equations

Theory and numerics of subspace approximation of eigenvalue problems

A Two-Scale Discretization Scheme for Mixed Variational Formulation of Eigenvalue Problems

The Partition of Unity Finite Element Method for the Schrödinger Equation

Finite-Difference Approximations and Local Algorithm for the Poisson and Poisson-Boltzmann Electrostatics

Finite Element Approximation to the Extremal Eigenvalue Problem for Inhomogenous Materials

Multiscale Finite Element Discretizations Based on Local Defect Correction for the Biharmonic Eigenvalue Problem of Plate Buckling

Multigrid Discretization and Iterative Algorithm for Mixed Variational Formulation of the Eigenvalue Problem of Electric Field

Discrete mechanics Based on Finite Element Methods

An Efficient Adaptive Finite Element Method for Eigenvalue Problems

Adaptive Finite Element Approximations for Kohn-Sham Models