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 General Second Order Complete Active Space Self-Consistent-field Solver for Large-Scale Systems

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

An Efficient DFT Solver for Nanoscale Simulations and Beyond

Self-Consistent Field Approach for the Variational Quantum Eigensolver: Orbital Optimization Goes Adaptive

Efficient single-grid and multi-grid solvers for real-space orbital-free density functional theory

Efficient and Robust Ab Initio Self-Consistent Field Acceleration Algorithm Based on a Semiempirical Model Hamiltonian

Automatic Selection of Active Orbitals from Generalized Valence Bond Orbitals.

Accelerating the Convergence of Auxiliary-Field Quantum Monte Carlo inSolids with Optimized Gaussian Basis Sets

Two-Level Iterative Solver For Linear Response Time-Dependent Density Functional Theory With Plane Wave Basis Set

MBE-CASSCF Approach for the Accurate Treatment of Large Active Spaces

A black-box, general purpose quadratic self-consistent field code with and without Cholesky Decomposition of the two-electron integrals

Self-consistent Field for Fragmented Quantum Mechanical Model of Large Molecular Systems.

Economical Quasi-Newton Self Consistent Field Solver

Ab initio Calculations in Atoms, Molecules, and Solids, Treating Spin-Orbit Coupling and Electron Interaction on Equal Footing

Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory Based $Ab$ $Initio$ Molecular Dynamics I: Theory, Algorithm, and Performance

Efficient real space formalism for hybrid density functionals

Two-Level Chebyshev Filter Based Complementary Subspace Method: Pushing the Envelope of Large-Scale Electronic Structure Calculations.

High-order finite element method for atomic structure calculations

A Semi-smooth Newton Method for Solving Semidefinite Programs in Electronic Structure Calculations

Adaptive Regularized Self-Consistent Field Iteration with Exact Hessian for Electronic Structure Calculation.