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.

Towards the solution of the many-electron problem in real materials: equation of state of the hydrogen chain with state-of-the-art many-body methods

Publisher’s Note: Towards the Solution of the Many-Electron Problem in Real Materials: Equation of State of the Hydrogen Chain with State-of-the-Art Many-Body Methods [phys. Rev. X 7, 031059 (2017)]

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

Unveiling the Finite Temperature Physics of Hydrogen Chains via Auxiliary Field Quantum Monte Carlo

Large Scale and Quantum Accurate Molecular Dynamics Simulation:Liquid Iron under Extreme Condition

First principles simulations of dense hydrogen

Structural and electronic properties of solid molecular hydrogen from many-electron theories

Towards the same line of liquid-liquid phase transition of dense hydrogen from various theoretical predictions

A Comparative Study Using State-of-the-art Electronic Structure Theories on Solid Hydrogen Phases under High Pressures

Theoretical and Experimental Constraints for the Equation of State of Dense and Hot Matter

Solutions of the Two Dimensional Hubbard Model: Benchmarks and Results from a Wide Range of Numerical Algorithms

Describing a Strongly Correlated Model System with Density Functional Theory.

Genuine converging solution of self-consistent field equations for extended many-electron systems

Solution to the key problem of statistical physics -- calculations of partition function of many-body systems

Stochastic density functional theory combined with Langevin dynamics for warm dense matter

Equation of State for Hard-Sphere Chain Molecules

Untangling Strongly and Weakly Interacting Configurations in Many-electron Wave Functions

Stable liquid Hydrogen at high pressure by a novel ab-initio molecular dynamics

Bound state breaking and the importance of thermal exchange-correlation effects in warm dense hydrogen

Ab initio path integral Monte Carlo simulations of warm dense two-component systems without fixed nodes: Structural properties