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.

Linear Scaling Electronic Structure Calculations With Numerical Atomic Basis Set

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

Electronic structure calculations and molecular dynamics simulations with linear system-size scaling

A Sub-Linear Scaling Algorithm for Computing the Electronic Structure of Materials

O(N) methods in electronic structure calculations

Linear-scaling ab-initio calculations for large and complex systems

Parallel Implementation of Large-Scale Linear Scaling Density Functional Theory Calculations with Numerical Atomic Orbitals in HONPAS.

HONPAS: A Linear Scaling Open‐source Solution for Large System Simulations

Linear Scaling Local Correlation Approach for Solving the Coupled Cluster Equations of Large Systems.

Linear scaling calculation of band edge states and doped semiconductors

Linear Scaling Calculation of Maximally Localized Wannier Functions with Atomic Basis Set.

Low scaling quantum chemistry software

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

Linear scaling Krylov subspace method for large scale {\it ab initio} electronic structure calculations of metals

Linear scaling computation of the Fock matrix VII. Periodic Density Functional Theory at the $Γ$-point

Linear scaling DFT calculations for large Tungsten systems using an optimized local basis

A review on non-relativistic fully numerical electronic structure calculations on atoms and diatomic molecules

Linear Scaling Discontinuous Galerkin Density Matrix Minimization Method with Local Orbital Enriched Finite Element Basis: 1-D Lattice Model System

Self-consistent field theory based molecular dynamics with linear system-size scaling

A Parallel Orbital-Updating Approach for Electronic Structure Calculations

Toward linear scaling auxiliary field quantum Monte Carlo with local natural orbitals