Stochastic density functional theory

Linear‐scaling implementations of density functional theory (DFT) reach their intended efficiency regime only when applied to systems having a physical size larger than the range of their Kohn–Sham density matrix (DM). This causes a problem since many types of large systems of interest have a rather broad DM range and are therefore not amenable to analysis using DFT methods. For this reason, the recently proposed stochastic DFT (sDFT), avoiding exhaustive DM evaluations, is emerging as an attractive alternative linear‐scaling approach. This review develops a general formulation of sDFT in terms of a (non)orthogonal basis representation and offers an analysis of the statistical errors (SEs) involved in the calculation. Using a new Gaussian‐type basis‐set implementation of sDFT, applied to water clusters and silicon nanocrystals, it demonstrates and explains how the standard deviation and the bias depend on the sampling rate and the system size in various types of calculations. We also develop a basis‐set embedded‐fragments theory, demonstrating its utility for reducing the SEs for energy, density of states and nuclear force calculations. Finally, we discuss the algorithmic complexity of sDFT, showing it has CPU wall‐time linear‐scaling. The method parallelizes well over distributed processors with good scalability and therefore may find use in the upcoming exascale computing architectures.

[1]  M. J. Gillan,et al.  Order-N first-principles calculations with the conquest code , 2007, Comput. Phys. Commun..

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

[3]  R. Baer,et al.  Equilibrium configurations of large nanostructures using the embedded saturated-fragments stochastic density functional theory. , 2017, The Journal of chemical physics.

[4]  D. Sánchez-Portal,et al.  The SIESTA method for ab initio order-N materials simulation , 2001, cond-mat/0111138.

[5]  K. Barros,et al.  Gradient-based stochastic estimation of the density matrix , 2017, 1711.10570.

[6]  D. Bowler,et al.  O(N) methods in electronic structure calculations. , 2011, Reports on progress in physics. Physical Society.

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

[8]  Leonard Kleinman,et al.  Efficacious Form for Model Pseudopotentials , 1982 .

[9]  R. Baer,et al.  Overlapped embedded fragment stochastic density functional theory for covalently-bonded materials. , 2018, The Journal of chemical physics.

[10]  Eran Rabani,et al.  Communication: Embedded fragment stochastic density functional theory. , 2014, The Journal of chemical physics.

[11]  Yousef Saad,et al.  Approximating Spectral Densities of Large Matrices , 2013, SIAM Rev..

[12]  Martin Head-Gordon,et al.  Sparsity of the Density Matrix in Kohn-Sham Density Functional Theory and an Assessment of Linear System-Size Scaling Methods , 1997 .

[13]  M. W. van der Kamp,et al.  Combined quantum mechanics/molecular mechanics (QM/MM) methods in computational enzymology. , 2013, Biochemistry.

[14]  Wolfram Koch,et al.  A Chemist's Guide to Density Functional Theory , 2000 .

[15]  D. Truhlar,et al.  QM/MM: what have we learned, where are we, and where do we go from here? , 2007 .

[16]  Giulia Galli,et al.  Large‐Scale Electronic Structure Calculations Using Linear Scaling Methods , 2000 .

[17]  R. Baer,et al.  Stochastic Density Functional Theory at Finite Temperatures , 2018, 1801.02163.

[18]  J. Kussmann,et al.  CHAPTER 1 Linear-Scaling Methods in Quantum Chemistry , 2006 .

[19]  J. Kussmann,et al.  Linear‐Scaling Methods in Quantum Chemistry , 2007 .

[20]  R. Baer,et al.  Self-averaging stochastic Kohn-Sham density-functional theory. , 2013, Physical review letters.

[21]  Matthias Scheffler,et al.  Efficient O(N) integration for all-electron electronic structure calculation using numeric basis functions , 2009, J. Comput. Phys..

[22]  Chen Huang,et al.  Potential-functional embedding theory for molecules and materials. , 2011, The Journal of chemical physics.

[23]  G. G. Hall The molecular orbital theory of chemical valency VIII. A method of calculating ionization potentials , 1951, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

[24]  Soler,et al.  Self-consistent order-N density-functional calculations for very large systems. , 1996, Physical review. B, Condensed matter.

[25]  Victor Guallar,et al.  Importance of accurate charges in molecular docking: Quantum mechanical/molecular mechanical (QM/MM) approach , 2005, J. Comput. Chem..

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

[27]  Sankey,et al.  Projected random vectors and the recursion method in the electronic-structure problem. , 1994, Physical review. B, Condensed matter.

[28]  Stefan Goedecker,et al.  Low Complexity Algorithms for Electronic Structure Calculations , 1995 .

[29]  Yihan Shao,et al.  An improved J matrix engine for density functional theory calculations , 2000 .

[30]  Dominika Zgid,et al.  Stochastic Self-Consistent Second-Order Green's Function Method for Correlation Energies of Large Electronic Systems. , 2017, Journal of chemical theory and computation.

[31]  David R. Bowler,et al.  A Comparison of Linear Scaling Tight Binding Methods , 1997 .

[32]  Kieron Burke,et al.  Partition Density Functional Theory , 2009, 0901.0942.

[33]  G. M. Stocks,et al.  Order-N multiple scattering approach to electronic structure calculations. , 1995, Physical review letters.

[34]  Colombo,et al.  Efficient linear scaling algorithm for tight-binding molecular dynamics. , 1994, Physical review letters.

[35]  Johannes M. Dieterich,et al.  Orbital-free density functional theory for materials research , 2018 .

[36]  常田 貴夫,et al.  Density functional theory in quantum chemistry , 2014 .

[37]  M. Hutchinson A stochastic estimator of the trace of the influence matrix for laplacian smoothing splines , 1989 .

[38]  Vanderbilt,et al.  Generalization of the density-matrix method to a nonorthogonal basis. , 1994, Physical review. B, Condensed matter.

[39]  Michael W. Schmidt,et al.  Noncovalent interactions in extended systems described by the effective fragment potential method: theory and application to nucleobase oligomers. , 2010, The journal of physical chemistry. A.

[40]  J. Arriaga Physica Status Solidi (B) Basic Research , 2005 .

[41]  Benny G. Johnson,et al.  THE CONTINUOUS FAST MULTIPOLE METHOD , 1994 .

[42]  Wang Calculating the density of states and optical-absorption spectra of large quantum systems by the plane-wave moments method. , 1994, Physical review. B, Condensed matter.

[43]  David E. Manolopoulos,et al.  Canonical purification of the density matrix in electronic-structure theory , 1998 .

[44]  David Alan Drabold,et al.  Maximum entropy approach for linear scaling in the electronic structure problem. , 1993, Physical review letters.

[45]  Thierry Deutsch,et al.  Challenges in large scale quantum mechanical calculations , 2016, 1609.00252.

[46]  R. Baer,et al.  Stochastic Optimally Tuned Range-Separated Hybrid Density Functional Theory. , 2015, The journal of physical chemistry. A.

[47]  Alán Aspuru-Guzik,et al.  Advances in molecular quantum chemistry contained in the Q-Chem 4 program package , 2014, Molecular Physics.

[48]  Sándor Suhai,et al.  Self-consistent-charge density-functional tight-binding method for simulations of complex materials properties , 1998 .

[49]  Martins,et al.  Efficient pseudopotentials for plane-wave calculations. , 1991, Physical review. B, Condensed matter.

[50]  E Weinan,et al.  Pole-Based approximation of the Fermi-Dirac function , 2009, 0906.1319.

[51]  Susanne Hertz,et al.  Advances in Quantum Chemistry , 2019, Quantum Systems in Physics, Chemistry and Biology - Theory, Interpretation, and Results.

[52]  Hernández,et al.  Self-consistent first-principles technique with linear scaling. , 1995, Physical review. B, Condensed matter.

[53]  Michele Parrinello,et al.  Quickstep: Fast and accurate density functional calculations using a mixed Gaussian and plane waves approach , 2005, Comput. Phys. Commun..

[54]  A. Szabó,et al.  Modern quantum chemistry : introduction to advanced electronic structure theory , 1982 .

[55]  R. Dreizler,et al.  Density Functional Theory: An Advanced Course , 2011 .

[56]  R. Abrol,et al.  An optimal adiabatic-to-diabatic transformation of the 1 2A′ and 2 2A′ states of H3 , 2002 .

[57]  Wei Wu,et al.  Density embedded VB/MM: a hybrid ab initio VB/MM with electrostatic embedding. , 2008, The journal of physical chemistry. A.

[58]  Martin,et al.  Unconstrained minimization approach for electronic computations that scales linearly with system size. , 1993, Physical review. B, Condensed matter.

[59]  J. Moussa,et al.  Minimax rational approximation of the Fermi-Dirac distribution. , 2016, The Journal of chemical physics.

[60]  R. Silver,et al.  Kernel polynomial method for a nonorthogonal electronic-structure calculation of amorphous diamond , 1997 .

[61]  G. Wellein,et al.  The kernel polynomial method , 2005, cond-mat/0504627.

[62]  Li,et al.  Density-matrix electronic-structure method with linear system-size scaling. , 1993, Physical review. B, Condensed matter.

[63]  Mark E. Tuckerman,et al.  A reciprocal space based method for treating long range interactions in ab initio and force-field-based calculations in clusters , 1999 .

[64]  R. Baer,et al.  Augmented Lagrangian method for order-N electronic structure , 2001 .

[65]  Clemens C. J. Roothaan,et al.  New Developments in Molecular Orbital Theory , 1951 .

[66]  上海科学技術文献出版社 数学年刊 = Chinese annals of mathematics , 1982 .

[67]  V. Barone,et al.  Chapter 2 - Extending the Range of Computational Spectroscopy by QM/MM Approaches: Time-Dependent and Time-Independent Routes , 2010 .

[68]  Taisuke Ozaki,et al.  Efficient low-order scaling method for large-scale electronic structure calculations with localized basis functions , 2010 .

[69]  M. Desjarlais,et al.  Frontiers and challenges in warm dense matter , 2014 .

[70]  S. Trickey,et al.  Frank Discussion of the Status of Ground-State Orbital-Free DFT , 2015, 1502.04017.

[71]  Alexei A Kananenka,et al.  Communication: Towards ab initio self-energy embedding theory in quantum chemistry. , 2015, The Journal of chemical physics.

[72]  Y. Saad,et al.  An estimator for the diagonal of a matrix , 2007 .

[73]  Martin Head-Gordon,et al.  Chebyshev expansion methods for electronic structure calculations on large molecular systems , 1997 .

[74]  Gustavo E. Scuseria,et al.  Linear Scaling Density Functional Calculations with Gaussian Orbitals , 1999 .

[75]  Chris-Kriton Skylaris,et al.  Introducing ONETEP: linear-scaling density functional simulations on parallel computers. , 2005, The Journal of chemical physics.

[76]  John R. Sabin,et al.  Combining quantum mechanics and molecular mechanics : some recent progresses in QM/MM methods , 2010 .

[77]  T. Frauenheim,et al.  DFTB+, a sparse matrix-based implementation of the DFTB method. , 2007, The journal of physical chemistry. A.