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 Strategy to Improve the Accuracy of Digital Simulation for Electroanalytical Chemistry

Direct Simulation for CAD Models Undergoing Parametric Modifications

A Multilevel Method With Novel Correction Strategy for Parallel Finite-Element Analysis of Electromagnetic Problems

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

Internal-growth Strategies for Dynamic Process Optimization with Differential-algebraic Equations

Digital simulation of chronoamperometric current at microdisk electrode and its verification

Time-frequency Analysis of the Spatial Distribution for Li+concentration in the Electrolyte Throughout the Full Battery

Computer simulations of electrodiffusion problems based on Nernst–Planck and Poisson equations

Single Alternating Group Explicit (sage) Method for Electrochemical Finite Difference Digital Simulation

Towards chemical accuracy using a multi-mesh adaptive finite element method in all-electron density functional theory

Numerical Solution of 3D Poisson-Nernst-Planck Equations Coupled with Classical Density Functional Theory for Modeling Ion and Electron Transport in a Confined Environment

Adaptive Finite Differencing in High Accuracy Electronic Structure Calculations

Recent Advances in Electrochemical Kinetics Simulations and Their Applications in Pt-based Fuel Cells

Accurate numerical simulation of electrodiffusion and water movement in brain tissue

Time filtering methods for electrohydrodynamics models

Interpolation Based Local Postprocessing for Adaptive Finite Element Approximations in Electronic Structure Calculations

A High-Accuracy Discretization Scheme for Numerical Device Simulations

A Simple and Efficient Lattice Summation Method for Metallic Electrodes in Constant Potential Molecular Dynamics Simulation

Accuracy Improvement of Immersed Boundary-Lattice Boltzmann and Finite Element Method by Iterative Velocity Correction

Establishment of a finite element model for extracting chemical reaction kinetics in a micro-flow injection system with high throughput sampling

Residual Based Error Estimator for Chemical-Mechanically Coupled Battery Active Particles