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 Semi-smooth Newton Method for Solving Semidefinite Programs in Electronic Structure Calculations

A SemiSmooth Newton Method for Semidefinite Programs and Its Applications in Electronic Structure Calculations.

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

A Newton-CG Augmented Lagrangian Method for Semidefinite Programming

A squared smoothing Newton method for semidefinite programming

SDPNAL+: A Matlab software for semidefinite programming with bound constraints (version 1.0)

A quadratically convergent semismooth Newton method for nonlinear semidefinite programming without generalized Jacobian regularity

Efficient Semidefinite Programming with approximate ADMM

A Parallel Reduced-Space Interior Point Method With Orthogonal Collocation for First-Swing Stability Constrained Emergency Control

Efficient numerical methods for computing the stationary states of phase field crystal models

PREDICTOR–CORRECTOR SMOOTHING NEWTON METHOD FOR SOLVING SEMIDEFINITE PROGRAMMING

An Efficient Inexact Symmetric Gauss-Seidel Based Majorized ADMM for High-Dimensional Convex Composite Conic Programming

Well-conditioned Primal-Dual Interior-point Method for Accurate Low-rank Semidefinite Programming

Enhancing Reduced Density Matrix Functional Theory Calculations by Coupling Orbital and Occupation Optimizations

Solving Low-Rank Semidefinite Programs via Manifold Optimization

Fast ADMM for Sum-of-Squares Programs Using Partial Orthogonality.

Towards Scalable Semidefinite Programming: Optimal Metric ADMM with A Worst-case Performance Guarantee

Inferring crystal electronic properties from experimental data sets through Semidefinite Programming

A Seamless Grid-Based Interface for Mean-Field QM/MM Coupled with Efficient Solvation Free Energy Calculations

Fast ADMM for Semidefinite Programs with Chordal Sparsity

A General Second Order Complete Active Space Self-Consistent-field Solver for Large-Scale Systems