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.

The HI Content and Extent of Low Surface Brightness Galaxies - Could LSB Galaxies be Responsible for Damped Ly-alpha Absorption?

A Concise Approach to the Calculation of the Polaron Ground-State Energy

Symmetries of the Kac-Peterson Modular Matrices of Affine Algebras

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

A formula and numerical study on Ewald 1D summation

Rigorous Error Bounds for Ewald Summation of Electrostatics at Planar Interfaces

Lekner summations and Ewald summations for quasi-two dimensional systems

Pressure Calculation in Polar and Charged Systems using Ewald Summation: Results for the Extended Simple Point Charge Model of Water

Damped reaction field method and the accelerated convergence of the real space Ewald summation

Error estimate of the u-series method for molecular dynamics simulations

Non-Ewald methods for evaluating the electrostatic interactions of charge systems: similarity and difference

Application of Ewald's Method for Efficient Summation of Dyon Long-Range Potentials

Kawamata-Viehweg vanishing and quint-canonical map of a Complex threefold

Unified analysis of finite-size error for periodic Hartree-Fock and second order Møller-Plesset perturbation theory

Ewald method for polytropic potentials in arbitrary dimensionality

Ewald summation for ferroelectric perovskites with charges and dipoles

Optimizing working parameters of the twin-range cutoff method in terms of accuracy and efficiency

Numerical Analysis of Finite Dimensional Approximations of Kohn-Sham Models

Coulomb-hole summations and energies for GW calculations with limited number of empty orbitals: a modified static remainder approach

The ENUF Method -- Ewald Summation based on Non-Uniform Fast Fourier Transform: Implementation, Parallelization, and Application

A note on accurate pressure calculations of Coulomb systems with periodic boundary conditions