Coupled cluster-inspired geminal wavefunctions

Time-reversal symmetry in RDMFT and pCCD with complex-valued orbitals

Reduced density matrix functional theory (RDMFT) and coupled cluster theory restricted to paired double excitations (pCCD) are emerging as efficient methodologies for accounting for the so-called non-dynamic electronic correlation effects. Up to now, molecular calculations have been performed with real-valued orbitals. However, before extending the applicability of these methodologies to extended systems, where Bloch states are employed, the subtleties of working with complex-valued orbitals and the consequences of imposing time-reversal symmetry must be carefully addressed. In this work, we describe the theoretical and practical implications of adopting time-reversal symmetry in RDMFT and pCCD when allowing for complex-valued orbital coefficients. The theoretical considerations primarily affect the optimization algorithms, while the practical implications raise fundamental questions about the stability of solutions. In particular, we find that complex solutions lower the energy when non-dynamic electronic correlation effects are pronounced. We present numerical examples to illustrate and discuss these instabilities and possible problems introduced by N-representability violations.

Graphical Approach to Interpreting and Efficiently Evaluating Geminal Wavefunctions

We consider wavefunctions built from antisymmetrized products of two-electron wavefunctions (geminals), which is arguably the simplest extension of the antisymmetrized product of one-electron wavefunctions (orbitals) (i.e., a Slater determinant). Extensive use of geminals in wavefunctions has been limited by their high cost stemming from the many combinations of the two-electron basis functions (orbital pairs) used to build the geminals. When evaluating the overlap of the APG wavefunction with an orthogonal Slater determinant, this cost can be interpreted as the cost of evaluating the permanent, resulting from the symmetry with respect to the interchange of orbital pairs, and the cost of assigning the occupied orbitals to the orbital pairs of the wavefunction. Focusing on the latter, we present a graphical interpretation of the Slater determinant and utilize the maximum weighted matching algorithm to estimate the combination of orbital pairs with the largest contribution to the overlap. Then, the cost due to partitioning the occupied orbitals in the overlap is reduced from 𝒪 ( ( N - 1 ) ! ! ) to 𝒪 ( N 3 log N ) . Computational results show that many of these combinations are not necessary to obtain an accurate solution to the wavefunction. Because the APG wavefunction is the most general of the geminal wavefunctions, this approach can be applied to any of the simpler geminal wavefunction ansätze. In fact, this approach may even be extended to generalized quasiparticle wavefunctions, opening the door to tractable wavefunctions built using components of arbitrary numbers of electrons, not just two electrons.

PyCI: A Python-scriptable library for arbitrary determinant CI

PyCI is a free and open-source Python library for setting up and running arbitrary determinant-driven configuration interaction (CI) computations, as well as their generalizations to cases where the coefficients of the determinant are nonlinear functions of optimizable parameters. PyCI also includes functionality for computing the residual correlation energy, along with the ability to compute spin-polarized one- and two-electron (transition) reduced density matrices. PyCI was originally intended to replace the ab initio quantum chemistry functionality in the HORTON library but emerged as a standalone research tool, primarily intended to aid in method development, while maintaining high performance so that it is suitable for practical calculations. To this end, PyCI is written in Python, adopting principles of modern software development, including comprehensive documentation, extensive testing, continuous integration/delivery protocols, and package management. Computationally intensive steps, notably operations related to generating Slater determinants and computing their expectation values, are delegated to low-level C++ code. This article marks the official release of the PyCI library, showcasing its functionality and scope.

ModelHamiltonian: A Python-scriptable library for generating 0-, 1-, and 2-electron integrals

ModelHamiltonian is a free, open source, and cross-platform Python library designed to express model Hamiltonians, including spin-based Hamiltonians (Heisenberg and Ising models) and occupation-based Hamiltonians (Pariser-Parr-Pople, Hubbard, and Hückel models) in terms of 1- and 2-electron integrals, so that these systems can be easily treated by traditional quantum chemistry software programs. ModelHamiltonian was originally intended to facilitate the testing of new electronic structure methods using HORTON but emerged as a stand-alone research tool that we recognize has wide utility, even in an educational context. ModelHamiltonian is written in Python and adheres to modern principles of software development, including comprehensive documentation, extensive testing, continuous integration/delivery protocols, and package management. While we anticipate that most users will use ModelHamiltonian as a Python library, we include a graphical user interface so that models can be built without programming, based on connectivity/parameters inferred from, for example, a SMILES string. We also include an interface to ChatGPT so that users can specify a Hamiltonian in plain language (without learning ModelHamiltonian's vocabulary and syntax). This article marks the official release of the ModelHamiltonian library, showcasing its functionality and scope.

Flexible Ansatz for N-Body Perturbation Theory

We propose a new perturbation theory framework that can be used to help with the projective solution of the Schrödinger equation for arbitrary wave functions. This Flexible Ansatz for N-body Perturbation Theory (FANPT) is based on our previously proposed Flexible Ansatz for the N-body Configuration Interaction (FANCI). We derive recursive FANPT expressions, including arbitrary orders in the perturbation hierarchy. We show that the FANPT equations are well-behaved across a wide range of conditions, including static correlation-dominated configurations and highly nonlinear wave functions.