Coupled cluster-inspired geminal wavefunctions

Electron pairs have an illustrious history in chemistry, from powerful concepts to understanding structural stability and reactive changes to the promise of serving as building blocks of quantitative descriptions of the electronic structure of complex molecules and materials. However, traditionally, two-electron wavefunctions (geminals) have not enjoyed the popularity and widespread use of the more standard single-particle methods. This has changed recently, with a renewed interest in the development of geminal wavefunctions as an alternative to describing strongly correlated phenomena. Hence, there is a need to find geminal methods that are accurate, computationally tractable, and do not demand significant input from the user (particularly via cumbersome and often ill-behaved orbital optimization steps). Here, we propose new families of geminal wavefunctions inspired by the pair coupled cluster doubles ansatz. We present a new hierarchy of two-electron wavefunctions that extends the one-reference orbital idea to other geminals. Moreover, we show how to incorporate single-like excitations in this framework without leaving the quasiparticle picture. We explore the role of imposing seniority restrictions on these wavefunctions and benchmark these new methods on model strongly correlated systems.

Communication: chemical bonding in carbon dimer isovalent series from the natural orbital functional theory perspective

The natural orbital functional theory admits two unique representations in the orbital space. On the one hand, we have the natural orbitals themselves that minimize the energy functional, and which afford for a diagonal one-particle reduced density matrix but not for a diagonal Lagrangian orbital energy multipliers matrix. On the other hand, since it is possible to reverse the situation but only once the energy minimization has been achieved, we have the so-called canonical representation, where the Lagrangian orbital energy multipliers matrix is diagonal but the one-particle reduced density matrix is not. Here it is shown that the former representation, the natural orbital representation, accounts nicely for the quadrupole bond character of the ground states of C2, BN, CB(-), and CN(+), and for the double bond order character of the isovalent (1)Σg (+) state of Si2. Similarly, the canonical orbital representation accounts correctly for the ionization spectra of all these species.