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.

Accurate Calculation of Equilibrium Reduced Density Matrix for the System-Bath Model: A Multilayer Multiconfiguration Time-Dependent Hartree Approach and Its Comparison to a Multi-Electronic-state Path Integral Molecular Dynamics Approach

Numerical computation of the equilibrium-reduced density matrix for strongly coupled open quantum systems

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

A Time-Dependent Random State Approach for Large-Scale Density Functional Calculations

Time-Dependent Quantum Transport: An Efficient Method Based On Liouville-Von-Neumann Equation For Single-Electron Density Matrix

Matrix-Product-State Algorithm for Finite Fractional Quantum Hall Systems

Self-consistent Density Functional Calculations without Diagonalization

Reduced Quantum Dynamics With Arbitrary Bath Spectral Densities: Hierarchical Equations Of Motion Based On Several Different Bath Decomposition Schemes

Ab-Initio Approach to Many-Body Quantum Spin Dynamics

A Unified Stochastic Formulation of Dissipative Quantum Dynamics. II. Beyond Linear Response of Spin Baths.

A multi-layer multi-configurational time-dependent Hartree approach to lattice models beyond one dimension

A New Method to Account for the Difference Between Classical and Quantum Baths in Quantum Dissipative Dynamics

Exact Quantum Master Equation Via the Calculus on Path Integrals

A new method for directly computing reduced density matrices

Continuous-time mean-variance efficiency: the 80% rule

Density Matrix Equation for a Bathed Small System and its Application to Molecular Magnets

Calculation Of State-To-State Cross Sections For Triatomic Reaction By The Multi-Configuration Time-Dependent Hartree Method

Quantum Density Mechanics: Accurate, purely density-based \textit{ab initio} implementation of many-electron quantum mechanics

Density-functional, density-functional tight-binding, and wave function calculations on biomolecular systems

Path Integral Molecular Dynamics for Exact Quantum Statistics of Multi-Electronic-State Systems

Solving the relativistic Hartree-Bogoliubov equation with the finite-difference method