DOCSLIB.ORG

Stochastic Density Functional Theory

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Stochastic Density Functional Theory Stochastic density functional theory Marcel David Fabian∗ Fritz Haber Center for Molecular Dynamics and Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 9190401, Israel Ben Shpiro∗ Fritz Haber Center for Molecular Dynamics and Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 9190401, Israel Eran Rabaniy Department of Chemistry, University of California, Berkeley, California 94720, USA, and Materials Science Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA, and The Sackler Center for Computational Molecular and Materials Science, Tel Aviv University, Tel Aviv 69978, Israel Daniel Neuhauserz Department of Chemistry and Biochemistry, University of California, Los Angeles California 90095, USA Roi Baerx Fritz Haber Center for Molecular Dynamics and Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 9190401, Israel 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 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. I. INTRODUCTION free DFT [36, 37] approaches using density-dependent kinetic energy functionals. The first two types of ap- proaches mentioned above are designed to answer ques- Density functional theory (DFT) is emerging as a tions typically asked about molecules, while for materials usefully-accurate general-purpose computational plat- and other large scale systems, we are more interested in form for predicting from first principles the ground-state coarse-grained properties. For example, with molecules, structure and properties of systems spanning a wide one is interested in bond orders, bond lengths spectral range of length scales, from single atoms and gas-phase lines; while for large systems we are more interested molecules, through macromolecules, proteins, nanocrys- in atomic densities, pair-correlation distributions (mea- tals, nanosheets, nanoribbons, surfaces, interfaces up to sured using neutron scattering) as well as charge/spin periodic or amorphous homogeneous or heterogeneous densities, polarizabilities and optical and electrical con- materials [1{5]. Significant efforts have been diverted to- ductivity. In molecules, we strive to understand each oc- wards the development of numerical and computational cupied/unoccupied Kohn-Sham eigenstate while in large methods enabling the use of DFT for studying exten- systems we are concerned with the density of hole and sive molecular systems. Several routes have been sug- electron states. gested: linear-scaling approaches [6{32], relying on the sparsity of the density matrix [33], DFT-based tight- Of course, detailed \molecular type" questions can also binding (DFTB) methods [34{36] which reduce the nu- arise in large systems, primarily when the processes of in- merical scaling using model Hamiltonians. Moreover, terest occur in small pockets or localized regions | for significant efforts have gone towards developing orbital- example, biochemical processes in proteins, localized cat- 2 alytic events on a surface, impurities in solids, etc. Here, A. Traditional basis-set formulation of Kohn-Sham a combination of methods, where the small subsystem equations with cubic scaling can be embedded in the larger environment is required. The Kohn-Sham (KS) density functional theory (KS- In this advanced review, we will focus on the stochas- DFT) is a molecular orbitals (MOs) approach which can tic DFT (sDFT) approach, developed using grids and be applied to a molecular system of N electrons using a plane-waves in recent years [38{42] but also based on e basis-set of atom-centered orbitals φ (r), α = 1; : : : ; K: ideas taken from works starting in the early 1990's, α The basis functions were developed to describe the elec- mainly within the tight-binding electronic structure tronic structure of the parent atom, and for molecules framework [43{50]. We make the point that the effi- they are the building blocks from which the orthonormal ciency of sDFT results from its adherence to answering MOs are built as superpositions: the coarse-grained \large system questions" mentioned above, rather than those asked for molecules. K X The new viewpoint taken here is that of stochastic (r) = φ (r) C ; n = 1; : : : ; K: (1) DFT using non-orthogonal localized basis-sets. The pri- n α αn α=1 mary motivation behind choosing local basis-sets is that they are considerably more compact than plane-waves In the simplest \population" model, each MO can either and therefore may enable studying significantly larger \occupy" two electrons (of opposing spin) or be empty. systems. Deterministic calculations using local basis-sets The occupied MOs (indexed as the first Nocc = Ne=2 are more readily applicable to large systems, and thus MOs) are used to form the total electron density: can generate useful benchmarks with which the statisti- Nocc cal errors and other properties characterizing sDFT can X 2 be studied in detail. n (r) = 2 × j n (r)j : (2) n The review includes three additional sections, further divided into subsections, to be described later. Sec- The coefficient matrix C in Eq. (1) can be obtained tion II reviews the theory and techniques used for non- from the variational principle applied to the Schr¨odinger orthogonal sDFT and studies in detail the statistical er- equation, leading to the Roothaan-Hall generalized eigen- rors and their dependence on sampling and system size. value equations [51, 52] (we follow the notations in In section III we explain the use of embedded fragments refs. [53{55]): and show their efficacy in reducing the stochastic errors FC = SCE: (3) of sDFT. Section IV summarizes and discusses the find- ings. Here, F = T + V en + J [n] + V xc [n] is the K × K KS Fock matrix, Sαα0 = hφαjφα0 i is the overlap matrix of the AO's and E is a diagonal matrix containing the MO energies, "1;:::;"K . The Fock matrix Fαα0 includes 1 2 the kinetic energy integrals, Tαα0 = φα − 2 r φα0 , en II. THEORY AND METHODS the nuclear attraction integrals Vαα0 = hφα jv^enj φα0 i, wherev ^en is the electron-nuclear interaction opera- tor, the Coulomb integrals Jαα0 = hφα jvH [n](r^)j φα0 i, 0 In this section, we discuss three formulations of KS- R n(r ) 3 0 where vH [n](r) = jr−r0j d r is the Hartree potential, DFT represented in non-orthogonal basis-sets. Since the xc and finally, the exchange-correlation integrals, V 0 = issue of algorithmic scaling is at the heart of develop- αα hφα jvxc [n](r^)j φα0 i where vxc [n](r^) is the exchange cor- ing DFT methods for large systems, we emphasize for relation potential. each formulation the associated algorithmic complexity In KS theory, the Fock matrix F and the electron den- (so-called system-size scaling). We start with the tra- sity n (r) are mutually dependent on each other and must ditional basis-set formulation of the Kohn-Sham equa- be obtained self-consistently. This is usually achieved by tions leading to standard cubic-scaling (subsection II A). converging an iterative procedure, Then, showing how, by focusing on observables and ex- ploiting the sparsity of the matrices, a quadratic-scaling · · · −! n (r) −! fvH [n](r) ; vxc [n](r)g −! (4) approach can be developed with no essential loss of rigor O(K3) or accuracy (subsection II B). Most of the discussion will −! F −−−−!fC; Eg −! n (r) −! :::; revolve around the third and final approach, stochastic DFT, which estimates expectation values using stochas- where in each iteration, a previous density iterate n (r) tic sampling methods, as described in subsection II C. is used to generate the Hartree vH [n](r) and exchange- This latter approach leads, to linear-scaling complexity. correlation vxc [n](r) potentials from which we construct 3 the Fock matrix F . Then, by solving Eq. (3) the coef- Here, S−1F is \plugged" in place of " into the function ficient matrix C is obtained from which a new density f of Eq. (7) [58]. Just like P is an operator, our method iterate n (r) is generated via Eqs. (1)-(2). The iterations also views S−1 as an operator which is applied to any continue until convergence (density stops changing with vector u with linear-scaling cost using a preconditioned a predetermined threshold), and a self-consistent field so- conjugate gradient method [59,60]). The operator P , ap- lution is thus obtained. plied to an arbitrary vector u, uses a Chebyshev expan- PNC l This implementation of the basis-set based approach sion [9,17,44,61] of length NC : P u = l=0 al (T; µ) u 0 −1 becomes computationally expensive for very large sys- where al are the expansion coefficients and u = S u, tems due to the cubic scaling of solving the algebraic u1 = Hu0 and then ul+1 = 2Hul − ul−1, l = 2; 3; :::.
Load more

Recommended publications
  • Arxiv:1805.00538V1 [Physics.Atm-Clus] 1 May 2018 Conduction Band Electrons Excited by an Oscillating Electric field, Typically a Photon
    First principles absorption spectra of Au nanoparticles: from quantum to classical Samuel Hernandeza∗, Yantao Xiab∗, Vojtˇech Vlˇceka, Robert Boutellea, Roi Baerc, Eran Rabanid,e, Daniel Neuhausera a Department of Chemistry and Biochemistry, University of California, Los Angeles, California 90095, USA.; b Department of Chemical and Biomolecular Engineering, University of California, Los Angeles, California 90095, USA.; c Fritz Haber Center for Molecular Dynamics and Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 9190401, Israel; d Department of Chemistry, University of California, Berkeley, California 94720, USA. e The Raymond and Beverly Sackler Center for Computational Molecular and Materials Science, Tel Aviv University, Tel Aviv 69978, Israel. ARTICLE HISTORY Compiled May 3, 2018 Abstract Absorption cross-section spectra for gold nanoparticles were calculated using fully quantum Stochastic Density Functional Theory and a classical Finite-Difference Time Domain (FDTD) Maxwell solver. Spectral shifts were monitored as a function of size (1.3{3.1 nm) and shape (octahedron, cubeoctahedron, and truncated cube). Even though the classical approach is forced to fit the quantum TDDFT at 3.1nm, at smaller sizes there is a significant deviation as the classical theory is unable to account for peak splitting and spectral blue shifts even after quantum spectral corrections. We attribute the failure of classical methods at predicting these features to quantum effects and low density of states in small nanoparticles. Classically, plasmon resonances are modeled as collective conduction electron excitations, but at small nanoparticle size these excitations transition to few or even individual conductive electron excitations, as indicated by our results. Introduction The unique physical and chemical properties of nanoparticles have generated intense academic and industrial interest, in hope that these properties, once well-understood, could be used for technological advances.
    [Show full text]
  • Stochastic Resolution of Identity for Real-Time Second-Order Green's
    Article Cite This: J. Chem. Theory Comput. XXXX, XXX, XXX−XXX pubs.acs.org/JCTC Stochastic Resolution of Identity for Real-Time Second-Order Green’s Function: Ionization Potential and Quasi-Particle Spectrum † ‡ † § ∥ ⊥ Wenjie Dou,*, Tyler Y. Takeshita,*, Ming Chen,*, , Roi Baer,*, Daniel Neuhauser,*, † § # and Eran Rabani*, , , † Department of Chemistry, University of California Berkeley, Berkeley, California 94720, United States ‡ Mercedes-Benz Research and Development North America, Sunnyvale, California 94085, United States § Materials Sciences Devision, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States ∥ Fritz Haber Research Center for Molecular Dynamics, Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 9190401, Israel ⊥ Department of Chemistry and Biochemistry, University of California Los Angeles, Los Angeles, California 90095, United States # The Raymond and Beverly Sackler Center of Computational Molecular and Materials Science, Tel Aviv University, Tel Aviv 69978, Israel ABSTRACT: We develop a stochastic resolution of identity approach to the real-time second-order Green’s function (real-time sRI-GF2) theory, extending our recent work for imaginary-time Matsubara Green’s function [Takeshita et al. J. Chem. Phys. 2019, 151, 044114]. The approach provides a framework to obtain the quasi-particle spectra across a wide range of frequencies and predicts ionization potentials and electron affinities. To assess the accuracy of the real-time sRI-GF2, we study a series of molecules and compare our results to experiments as well as to a many-body perturbation approach based on the GW approximation, where we find that the real-time sRI-GF2 is as accurate as self-consistent GW. The stochastic formulation reduces the formal computatinal O N 5 O N 3 N scaling from ( e ) down to ( e ) where e is the number of electrons.
    [Show full text]
  • Forces from Stochastic Density Functional Theory Under Nonorthogonal Atom-Centered Basis Sets
    Forces from stochastic density functional theory under nonorthogonal atom-centered basis sets Ben Shpiro Fritz Haber Center for Molecular Dynamics and Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 9190401, Israel Marcel David Fabian Fritz Haber Center for Molecular Dynamics and Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 9190401, Israel Eran Rabani∗ Department of Chemistry, University of California, Berkeley, California 94720, USA Materials Science Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA and The Sackler Center for Computational Molecular and Materials Science, Tel Aviv University, Tel Aviv 69978, Israel Roi Baer† Fritz Haber Center for Molecular Dynamics and Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 9190401, Israel We develop a formalism for calculating forces on the nuclei within the linear-scaling stochastic density functional theory (sDFT) in a nonorthogonal atom-centered basis-set representation (Fabian et al. WIREs Comput Mol Sci. 2019;e1412. https://doi.org/10.1002/wcms.1412) and apply it to Tryptophan Zipper 2 (Trp-zip2) peptide solvated in water. We use an embedded-fragment approach to reduce the statistical errors (fluctuation and systematic bias), where the entire peptide is the main fragment and the remaining 425 water molecules are grouped into small fragments. We analyze the magnitude of the statistical errors in the forces and find that the systematic bias is of the order of 3 0.065 eV/Å( 1.2 10− Eh/a0) when 120 stochastic orbitals are used, independently of systems ∼ × size. This magnitude of bias is sufficiently small to ensure that the bond lengths estimated by stochastic DFT (within a Langevin molecular dynamics simulation) will deviate by less than 1% from those predicted by a deterministic calculation.
    [Show full text]
  • Program Sunday 17/2 Nonadiabatic Dynamics Chair: Marco Garavelli 11:00 11:10 Sanford Ruhman Opening Remarks 11:10 11:50 Eberhard K
    Dead Sea Workshop on Excitons in Natural and Manmade Materials Ein Gedi, February 17-20 2013 Program Sunday 17/2 Nonadiabatic Dynamics Chair: Marco Garavelli 11:00 11:10 Sanford Ruhman Opening Remarks 11:10 11:50 Eberhard K. U. Gross MPI for Microstructure Physics, What's the correct classical force on the nuclei: A fresh look at Halle, Germany potential energy surfaces and Berry phases in the time domain 11:50 12:30 Siva Umapathy Indian Institute of Science, Frequency and time-resolved resonance Raman approaches to Bangalore, India understanding conical intersections: examples and potential experiments 12:30 13:10 Giulio Cerullo Politecnico di Milano, Italy Real-time wavepacket dynamics through a conical intersection: the primary event of vision 13:10 14:40 Lunch 14:40 15:20 Wolfgang Domcke Technische Universität Computational studies of the photoinduced homolytic München, Germany dissociation of water in hydrogen-bonded chromophore-water complexes 15:20 16:00 Bern Kohler Montana State University, USA Unraveling excited state localization and decay in single- and double-stranded nucleic acids 16:00 16:40 Alan Aspuru Guzik Harvard University, Towards the happy marriage of correlated electronic structure Massachusetts, USA and open quantum systems theory: correlated polaron theory 16:40 17:00 Coffee Magnetoreception and Avian Navigation Chair: Ksenia Bravaya 17:00 17:40 Abraham Nitzan Tel Aviv University, Israel Excitonic effects in molecular electronics and molecular plasmonics 17:40 18:20 Kiminori Maeda University of Oxford, UK Magnetically
    [Show full text]