Representation of the virtual space in extended systems – a correlation energy convergence study
We present an investigation of the convergence behaviour of the local second-order Møller-Plesset perturbation theory (MP2) correlation energy toward the canonical result for three insulating crystals with either projected atomic orbitals (PAOs) or various orthonormal representations of the virtual orbital space. Echoing recent results for finite molecular systems, we find that significantly fewer PAOs than localised orthonormal virtual orbitals are required to reproduce the canonical correlation energy. We find no clear-cut correlation between conventional measures of orbital locality and the ability of the representation to span the excitation space of local domains. We show that the PAOs of the reference unit cell span parts of the excitation space that can only be reached with distant local orthonormal virtual orbitals.
A page showing Hydejack-specific markdown content.
Hydejack offers a few additional features to markup your markdown. Don’t worry, these are merely CSS classes added with kramdown’s {:...} syntax, so that your content remains compatible with other Jekyll themes.
A page showing Hydejack-specific markdown content.
Hydejack offers a few additional features to markup your markdown. Don’t worry, these are merely CSS classes added with kramdown’s {:...} syntax, so that your content remains compatible with other Jekyll themes.
Divide–Expand–Consolidate Second-Order Møller–Plesset Theory with Periodic Boundary Conditions
We present a generalization of the divide–expand–consolidate (DEC) framework for local coupled-cluster calculations to periodic systems and test it at the second-order Møller–Plesset (MP2) level of theory. For simple model systems with periodicity in one, two, and three dimensions, comparisons with extrapolated molecular calculations and the local MP2 implementation in the Cryscor program show that the correlation energy errors of the extended DEC (X-DEC) algorithm can be controlled through a single parameter, the fragment optimization threshold. Two computational bottlenecks are identified: the size of the virtual orbital spaces and the number of pair fragments required to achieve a given accuracy of the correlation energy. For the latter, we propose an affordable algorithm based on cubic splines interpolation of a limited number of pair-fragment interaction energies to determine a pair cutoff distance in accordance with the specified fragment optimization threshold.
Diagrammatic derivation of the Coupled-Cluster equations
The Coupled-Cluster (CC) equations at the triples level and beyond is generally cumbersome to work with due to the number of terms and complexity of the contractions. Both the book-keeping and the actual implementation of the contractions are prone to errors.
