Abstract

We propose and implement an algorithm to calculate the norm and reduced density matrices of the antisymmetrized geminal power (AGP) of any rank with polynomial cost. Our method scales quadratically per element of the reduced density matrices. Numerical tests indicate that our method is very fast and capable of treating systems with a few thousand orbitals and hundreds of electrons reliably in double-precision. In addition, we present reconstruction formulae that allows one to decompose higher order reduced density matrices in terms of linear combinations of lower order ones, thereby reducing the computational cost significantly.

Explore this paper’s authors, topics and related work.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Reviews 0

Write a reviewWrite

No reviews yet.

Discussion 0

Add a commentWrite

No comments yet.