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.

A direct calculation of eigenvalues of any arbitrary range of energy within local-subspace

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

Local Self-Energy Approach For Electronic Structure Calculations

Computing Generalized Eigenfunctions in Rigged Hilbert Spaces

Efficient Spectral and Spectral Element Methods for Eigenvalue Problems of Schrödinger Equations with an Inverse Square Potential

Properties of Definite Bethe-Salpeter Eigenvalue Problems

Alternative Approach to the Optimized Effective Potential Method

Modified-operator method for the calculation of band diagrams of crystalline materials

An Efficient Step-Forward Way To Solve The Schrodinger Eigenvalue Equation In Self-Consistent Calculations

Grid method for computation of generalized spheroidal wave functions based on discrete variable representation.

Energy spectrum and applications of Eckart plus Hellmann potential in hyperspherical coordinates

Computing eigenfunctions and eigenvalues of boundary value problems with the orthogonal spectral renormalization method

Strict Relationship: Potential - energy levels

Extended Eckart Theorem and New Variation Method for Excited States of Atoms

Subspace recursive Fermi-operator expansion strategies for large-scale DFT eigenvalue problems on HPC architectures

Direct energy functional minimization under orthogonality constraints

Efficient Hartree-Fock Exchange Algorithm with Coulomb Range Separation and Long-Range Density Fitting

Theory and numerics of subspace approximation of eigenvalue problems

Efficient HF exchange evaluation through Fourier convolution in Cartesian grid for orbital-dependent density functionals

FEAST nonlinear eigenvalue algorithm for $GW$ quasiparticle equations

Analytic Expression for Exact Ground State Energy Based on an Operator Method for a Class of Anharmonic Potentials