Ab initio linear scaling response theory: Electric polarizability by perturbed projection

Abstract

A linear scaling method for calculation of the static {\em ab inito} response within self-consistent field theory is developed and applied to calculation of the static electric polarizability. The method is based on density matrix perturbation theory [Niklasson and Challacombe, cond-mat/0311591], obtaining response functions directly via a perturbative approach to spectral projection. The accuracy and efficiency of the linear scaling method is demonstrated for a series of three-dimensional water clusters at the RHF/6-31G** level of theory. Locality of the response under a global electric field perturbation is numerically demonstrated by approximate exponential decay of derivative density matrix elements.

arXiv comment4.25 pages in PRL format, 2 figures

References in corpus (3)

• First principles calculation of vibrational Raman spectra in large systems: signature of small rings in crystalline SiO2
• Expansion algorithm for the density matrix
• Density Matrix Perturbation Theory

Cited by in corpus (12)

• O(N) methods in electronic structure calculations
• Density Matrix Perturbation Theory
• Graph-based linear scaling electronic structure theory
• Next generation extended Lagrangian first principles molecular dynamics
• Molecular-orbital-free algorithm for excited states in time-dependent perturbation theory
• Canonical density matrix perturbation theory
• Linear scaling computation of the Fock matrix VII. Periodic Density Functional Theory at the $\Gamma$-point
• Linear scaling calculation of band edge states and doped semiconductors
• Graph-based Quantum Response Theory and Shadow Born-Oppenheimer Molecular Dynamics
• Efficient Parallel Linear Scaling Method to get the Response Density Matrix in All-Electron Real-Space Density-Functional Perturbation Theory
• Notes on density matrix perturbation theory
• Linear scaling calculation of a n-type GaAs quantum dot