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.

Novel linear algebraic theory and one-hundred-million-atom quantum material simulations on the K computer

Ten-million-atom electronic structure calculations on the K computer with a massively parallel order-N theory

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

One-hundred-nm-scale electronic structure and transport calculations of organic polymers on the K computer

Efficient and Accurate Linear Algebraic Methods for Large-scale Electronic Structure Calculations with Non-orthogonal Atomic Orbitals

Redesigning OpenKMC for Multi-Component Trillion-Atom Simulations on the New Sunway Supercomputer

2.5 Million-Atom Ab Initio Electronic-Structure Simulation of Complex Metallic Heterostructures with DGDFT

Quantum Krylov-Subspace Method Based Linear Solver

Large-Scale Simulation of Quantum Computational Chemistry on a New Sunway Supercomputer

10-Million Atoms Simulation of First-Principle Package LS3DF

152K-computer-node parallel scalable implicit solver for dynamic nonlinear earthquake simulation

Atomistic simulation of low-dimensional nanostructures toward extreme-scale supercomputing

Parallel Implementation of Large-Scale Linear Scaling Density Functional Theory Calculations with Numerical Atomic Orbitals in HONPAS.

A robust quantum nonlinear solver based on the asymptotic numerical method

ELSI: A Unified Software Interface for Kohn-Sham Electronic Structure Solvers

Real-Time Krylov Theory for Quantum Computing Algorithms

Electronic structure calculations and molecular dynamics simulations with linear system-size scaling

Scalable Program Implementation and Simulation of the Large-Scale Quantum Algorithm: $1024\times 1024$ Quantum Linear Solver and Beyond

Multiscale quantum algorithms for quantum chemistry

Development of simulation package 'ELSES' for extra-large-scale electronic-structure calculation

The Improved Krylov Subspace Methods for Large and Sparse Linear Systems on Bulk Synchronous Parallel Architectures