Analysis of the divide-and-conquer method for electronic structure calculations

We study the accuracy of the divide-and-conquer method for electronic structure calculations. The analysis is conducted for a prototypical subdomain problem in the method. We prove that the pointwise difference between electron densities of the global system and the subsystem decays exponentially as a function of the distance away from the boundary of the subsystem, under the gap assumption of both the global system and the subsystem. We show that gap assumption is crucial for the accuracy of the divide-and-conquer method by numerical examples. In particular, we show examples with the loss of accuracy when the gap assumption of the subsystem is invalid.

[1]  William W. Hager,et al.  Multilevel domain decomposition for electronic structure calculations , 2007, J. Comput. Phys..

[2]  R. Martin,et al.  Electronic Structure: Basic Theory and Practical Methods , 2004 .

[3]  P. Hohenberg,et al.  Inhomogeneous Electron Gas , 1964 .

[4]  Juan C. Meza,et al.  A Linear Scaling Three Dimensional Fragment Method for Large ScaleElectronic Structure Calculations , 2008 .

[5]  S. Goedecker.  Linear scaling electronic structure methods , 1999 .

[6]  W. Kohn,et al.  Self-Consistent Equations Including Exchange and Correlation Effects , 1965 .

[7]  D R Bowler,et al.  Calculations for millions of atoms with density functional theory: linear scaling shows its potential , 2009, Journal of physics. Condensed matter : an Institute of Physics journal.

[8]  Weitao Yang,et al.  A density‐matrix divide‐and‐conquer approach for electronic structure calculations of large molecules , 1995 .

[9]  Shmuel Agmon,et al.  On kernels, eigenvalues, and eigenfunctions of operators related to elliptic problems , 1965 .

[10]  Kohn,et al.  Density functional and density matrix method scaling linearly with the number of atoms. , 1996, Physical review letters.

[11]  Yang,et al.  Direct calculation of electron density in density-functional theory: Implementation for benzene and a tetrapeptide. , 1991, Physical review. A, Atomic, molecular, and optical physics.

[12]  R. Parr.  Density-functional theory of atoms and molecules , 1989 .

[13]  William W. Hager,et al.  Domain Decomposition and Electronic Structure Computations: A Promising Approach , 2008 .

[14]  A. Nakano,et al.  Divide-and-conquer density functional theory on hierarchical real-space grids: Parallel implementation and applications , 2008 .

[15]  Jianfeng Lu,et al.  The Electronic Structure of Smoothly Deformed Crystals: Wannier Functions and the Cauchy–Born Rule , 2011 .

[16]  Jacques Periaux,et al.  Numerical Analysis and Scientific Computing for PDEs and their Challenging Applications , 2005 .

[17]  Lin-wang Wang,et al.  A divide-and-conquer linear scaling three-dimensional fragment method for large scale electronic structure calculations , 2008 .

[18]  Michele Benzi,et al.  Decay Properties of Spectral Projectors with Applications to Electronic Structure , 2012, SIAM Rev..

[19]  Masato Kobayashi,et al.  Divide-and-conquer-based linear-scaling approach for traditional and renormalized coupled cluster methods with single, double, and noniterative triple excitations. , 2009, The Journal of chemical physics.

[20]  Rajiv K. Kalia,et al.  Large-scale atomistic simulations of nanostructured materials based on divide-and-conquer density functional theory , 2011 .

[21]  Ryo Kobayashi,et al.  Linear scaling algorithm of real-space density functional theory of electrons with correlated overlapping domains , 2012, Comput. Phys. Commun..

[22]  Yang,et al.  Direct calculation of electron density in density-functional theory. , 1991, Physical review letters.

[23]  W. Kohn,et al.  Nearsightedness of electronic matter. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[24]  Michael Aizenman,et al.  Moment analysis for localization in random Schrödinger operators , 2003, math-ph/0308023.

[25]  E Weinan,et al.  The Kohn-Sham Equation for Deformed Crystals , 2012 .

[26]  E Weinan,et al.  Linear-scaling subspace-iteration algorithm with optimally localized nonorthogonal wave functions for Kohn-Sham density functional theory , 2009 .