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.

Tuning the range separation parameter in periodic systems

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

Study on tunable periodic graphene split-ring resonator (GSRR) structures using FDTD method

Band gaps of crystalline solids from Wannier-localization–based optimal tuning of a screened range-separated hybrid functional

Range-separated double-hybrid density-functional theory applied to periodic systems

Hybrid functionals for periodic systems in the density functional tight-binding method

Quasiparticle spectra from a non-empirical optimally-tuned range-separated hybrid density functional

All-electron self-consistent GW in the Matsubara-time domain: implementation and benchmarks of semiconductors and insulators

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

Exact Exchange with Range-separated Algorithm for Thermodynamic Limit of Periodic Hartree-Fock Theory

Realizing Effective Cubic-Scaling Coulomb Hole Plus Screened Exchange Approximation in Periodic Systems Via Interpolative Separable Density Fitting with a Plane-Wave Basis Set.

Improved Band Gaps and Structural Properties from Wannier-Fermi-Löwdin Self-Interaction Corrections for Periodic Systems

Range-separated density-functional theory with random phase approximation applied to noncovalent intermolecular interactions

Precision benchmarks for solids: G0W0 calculations with different basis sets

Understanding Band Gaps of Solids in Generalized Kohn-Sham Theory.

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

Tuning Range-Separated Density Functional Theory for Photocatalytic Water Splitting Systems

Nonempirical Range-Separated Hybrid Functional with Spatially Dependent Screened Exchange

Ionization Potentials and Fundamental Gaps in Atomic Systems from the Ensemble-DFT Approach

Range Separation Method for Density Functional Theory Based on Two-Electron Infinite-Order Two-Component Hamiltonian

Approaching Periodic Systems in Ensemble Density Functional Theory via Finite One-Dimensional Models