Direct energy functional minimization under orthogonality constraints

Abstract

The direct energy functional minimization problem in electronic structure theory, where the single-particle orbitals are optimized under the constraint of orthogonality, is explored. We present an orbital transformation based on an efficient expansion of the inverse factorization of the overlap matrix that keeps orbitals orthonormal. The orbital transformation maps the orthogonality constrained energy functional to an approximate unconstrained functional, which is correct to some order in a neighborhood of an orthogonal but approximate solution. A conjugate gradient scheme can then be used to find the ground state orbitals from the minimization of a sequence of transformed unconstrained electronic energy functionals. The technique provides an efficient, robust, and numerically stable approach to direct total energy minimization in first principles electronic structure theory based on tight-binding, Hartree-Fock, or density functional theory. For sparse problems, where both the orbitals and the effective single-particle Hamiltonians have sparse matrix representations, the effort scales linearly with the number of basis functions N in each iteration. For problems where only the overlap and Hamiltonian matrices are sparse the computational cost scales as O(M-2 N), where M is the number of occupied orbitals. We report a single point density functional energy calculation of a DNA decamer hydrated with 4003 water molecules under periodic boundary conditions. The DNA fragment containing a cis-syn thymine dimer is composed of 634 atoms and the whole system contains a total of 12 661 atoms and 103 333 spherical Gaussian basis functions. The direct energy functional minimization problem in electronic structure theory, where the single-particle orbitals are optimized under the constraint of orthogonality, is explored. We present an orbital transformation based on an efficient expansion of the inverse factorization of the overlap matrix that keeps orbitals orthonormal. The orbital transformation maps the orthogonality constrained energy functional to an approximate unconstrained functional, which is correct to some order in a neighborhood of an orthogonal but approximate solution. A conjugate gradient scheme can then be used to find the ground state orbitals from the minimization of a sequence of transformed unconstrained electronic energy functionals. The technique provides an efficient, robust, and numerically stable approach to direct total energy minimization in first principles electronic structure theory based on tight-binding, Hartree-Fock, or density functional theory. For sparse problems, where both the orbitals and the effective single-particle Hamiltonians have sparse matrix representations, the effort scales linearly with the number of basis functions N in each iteration. For problems where only the overlap and Hamiltonian matrices are sparse the computational cost scales as O͑M 2 N͒, where M is the number of occupied orbitals. We report a single point density functional energy calculation of a DNA decam...

Efficient Construction of Nonorthogonal Localized Molecular Orbitals in Large Systems†

“…[1][2][3][4][5][6][7][8][9][10][11][12][13][14][15][16][17][18] Development of linear-scaling techniques, in particular for excited states, remains an active area of electronic structure theory. With sparse matrix techniques, linear-scaling calculations can be achieved by means of DAC-style fragment methods, localized molecular orbitals (LMOs), and density matrices.…”

Fast, accurate evaluation of exact exchange: The occ-RI-K algorithm

“…The basic idea of computing just K iν or, for that matter, just F iν is of course extremely simple and undoubtedly has been thought of before. It has been employed in the iterative updates of plane-wave DFT calculations, 60 where the full Fock matrix is enormous and thus, evaluating a narrow rectangular matrix is far preferable to an enormous square matrix. Indeed, the idea of constructing this economized version of K was mentioned briefly by Aquilante et al, 61 but we have found no prior or subsequent discussion in the literature.…”

Efficient Linear-Scaling Density Functional Theory for Molecular Systems

“…To demonstrate the computational efficiency of the ALMO methods their performance is compared to that of the orbital transformation (OT) approach 47,48 . OT is a highly efficient and well-optimized SCF algorithm that performs PCG optimization of the occupied MOs.…”