Abstract
The computational efficiency of local correlation methods is strongly dependent on the size of the domain of functions used to expand local correlating orbitals such as orbital specific or pair natural orbitals. Here, we define a principal domain of order m as the subset of m one-particle functions that provides the best support for a given n-electron wave function by maximizing the partial trace of the one-body reduced density matrix. Principal domains maximize the overlap between the wave function and its approximant for two-electron systems and are the domain selection equivalent of Löwdin's natural orbitals. We present an efficient linear scaling greedy algorithm for obtaining principal domains of projected atomic orbitals and demonstrate its utility in the context of the pair natural orbital local correlation theory. We numerically determine thresholds such that the projected atomic orbital domain error is an order of magnitude smaller than the pair natural orbital truncation error.
Sign-in/Register to access full text options 
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a callDisclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.
