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.

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

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

Multiscale analysis and numerical algorithm for the Schrödinger equations in heterogeneous media

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

Combined Semiclassical And Effective-Mass Schrodinger Approach For Multiscale Analysis Of Semiconductor Nanostructures

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

A nonlinear poisson-schrodinger solver under cylindrical coordinate for quantum effect in nanowire MOSFET

A Multiscale Method For Optical Responses Of Nanostructures

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

A new analytical method for modeling a 2D electrostatic potential in MOS devices, applicable to compact modeling

A finite volume method for the multi subband Boltzmann equation with realistic 2D scattering in double gate MOSFETs

A Multiscale Model Reduction Method for Partial Differential Equations

A pseudo-spectral multiscale method: Interfacial conditions and coarse grid equations

Two-scale Analysis for Multiscale Landau-Lifshitz-Gilbert Equation: Theory and Numerical Methods

Low order nonconforming finite element method for time-dependent nonlinear Schrödinger equation

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

A Time Two-Mesh Compact Difference Method for the One-Dimensional Nonlinear Schrödinger Equation

A Finite Element Method with Energy-Adaptive Grids for the Coupled Schrödinger-Poisson-drift-diffusion Model

A New Coordinate Conversion Method for Spatial Transmission of Scattering Electrons

An Adaptive Grid Algorithm for Self-Consistent K·p Schrodinger and Poisson Equations in UTB InSb-based Pmosfets

Solving Dirac Equations on A 3d Lattice with Inverse Hamiltonian and Spectral Methods