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.

Waveform Krylov Subspace Methods for Tightly Coupled Systems

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

An Adaptive Wavelet Method for Nonlinear Circuit Simulation

Nonlinear circuit simulation based on adaptive wavelet method

On convergence of waveform relaxation for nonlinear systems of ordinary differential equations

Convergence Analysis of a Krylov Subspace Spectral Method for the 1-D Wave Equation in an Inhomogeneous Medium

Adaptive wavelet-based method for simulation of electronic circuits

Convergence of Waveform Relaxation Methods for Hermitian Positive Definite Linear Systems

A Fast Wavelet Collocation Method for High-Speed Circuit Simulation

An Efficient Sylvester Equation Solver for Time Domain Circuit Simulation by Wavelet Collocation Method.

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

Subspace-Based Approach for Circuit Simulation

An efficient off-step spline algorithm for wave simulation of nonlinear Kuramoto–Sivashinsky and Korteweg-de Vries equations

A mixed finite element method for solving coupled wave equation of Kirchhoff type with nonlinear boundary damping and memory term

Two linear energy-preserving compact finite difference schemes for coupled nonlinear wave equations

A Wavelet-Balance Approach for Steady-State Analysis of Nonlinear Circuits

Waveform Transmission Method, a New Waveform-relaxation Based Algorithm to Solve Ordinary Differential Equations in Parallel

A wavelet method for solving a class of nonlinear boundary value problems

Adaptive Multi-rate Wavelet Method for Circuit Simulation

The Time Domain Steady State Simulation for Nonlinear Dynamic Circuit Based on FWCM