Abstract:Fermi operator expansion (FOE) methods are powerful alternatives to diagonalization type methods for solving Kohn-Sham density functional theory (KSDFT). One example is the pole expansion and selected inversion (PEXSI) method, which approximates the Fermi operator by rational matrix functions and reduces the computational complexity to at most quadratic scaling for solving KSDFT. Unlike diagonalization type methods, the chemical potential often cannot be directly read off from the result of a single step of evaluation of the Fermi operator. Hence multiple evaluations are needed to be sequentially performed to compute the chemical potential to ensure the correct number of electrons within a given tolerance. This hinders the performance of FOE methods in practice. In this paper, we develop an efficient and robust strategy to determine the chemical potential in the context of the PEXSI method. The main idea of the new method is not to find the exact chemical potential at each self-consistent-field (SCF) iteration but to dynamically and rigorously update the upper and lower bounds for the true chemical potential, so that the chemical potential reaches its convergence along the SCF iteration. Instead of evaluating the Fermi operator for multiple times sequentially, our method uses a two-level strategy that evaluates the Fermi operators in parallel. In the regime of full parallelization, the wall clock time of each SCF iteration is always close to the time for one single evaluation of the Fermi operator, even when the initial guess is far away from the converged solution. We demonstrate the effectiveness of the new method using examples with metallic and insulating characters, as well as results from ab initio molecular dynamics. Published by AIP Publishing.

Closed-form solution for inverse problems of Fermi systems (vol 48, pg 1558, 1993)

Unified Expression for Fermi and Bose Distributions

Fermi-Dirac Integrals in Degenerate Regimes: A Novel Asymptotic Expansion

Numerical Methods for the Inverse Problem of Density Functional Theory

Finite-Temperature Evaluation of the Fermi Density Operator

Quantitative relaxation towards equilibrium for solutions to the Boltzmann-Fermi-Dirac equation with cutoff hard potentials

Inverse problems on Fermi systems and ionic crystals

Solution of the Dirichlet problem by a finite difference analog of the boundary integral equation

Spectral functions of the strongly interacting 3D Fermi gas

Highly Efficient and Accurate Sum-over-poles Expansion of Fermi and Bose Functions at Near Zero Temperatures: Fano Spectrum Decomposition Scheme.

Closed-form solutions of the nonlinear Schrödinger equation with arbitrary dispersion and potential

Solving Dirac Equations on A 3d Lattice with Inverse Hamiltonian and Spectral Methods

A Symmetry of Boundary Functions Method for Solving the Backward Time-Fractional Diffusion Problems

Robust Determination Of The Chemical Potential In The Pole Expansion And Selected Inversion Method For Solving Kohn-Sham Density Functional Theory

Inverse problem in energy-dependent potentials using semiclassical methods

Analytical form of the density of states of a Fermi gas in low dimensions

A Geometric Approach to the Landauer-Büttiker Formula

A New Method for Inverse Black Body Radiation Problem

Exact Solution of the Problem of a Two-Dimensional Dirac Particle with Flux 1/2.

Global well-posedness of Boltzmann-Fermi-Dirac equation for hard potential

Population of Closed-Channel Molecules in Trapped Fermi Gases with Broad Feshbach Resonances