Ground state eigenvectors of the reduced Bardeen-Cooper-Schrieffer Hamiltonian are employed as a wavefunction ansatz to model strong electron correlation in quantum chemistry. This wavefunction is a product of weakly-interacting pairs of electrons.While other geminal wavefunctions may only be employed in a projected Schrödinger equation, the present approach may be solved variationally with polynomial cost. The resulting wavefunctions are used to compute expectation values of CoulombHamiltionans and we present results for atoms and dissociation curves which are in agreement with doubly-occupied configuration interaction (DOCI) data. The present approach will serve as the starting point for a many-body theory of pairs, much as Hartree-Fock is the starting point for weakly-correlated electrons.

Eigenvectors of the reduced Bardeen–Cooper–Schrieffer Hamiltonian have recently been employed as a variational wavefunction ansatz in quantum chemistry. This wavefunction is a mean-field of pairs of electrons (geminals). In this contribution, we report optimal expressions for their reduced density matrices in both the original physical basis and the basis of the Richardson–Gaudin pairs. Physical basis expressions were originally reported by Gorohovsky and Bettelheim [Phys. Rev. B 84, 224503 (2011)]. In each case, the expressions scale like O(N4), with the most expensive step being the solution of linear equations. Analytic gradients are also reported in the physical basis. These expressions are an important step toward practical mean-field methods to treat strongly correlated electrons.

In this work, we have developed a simple synthetic approach using Et3N·3HF as an alternative to the DAST reagent. We controlled the stereochemistry of the nucleophilic fluorination at C4 of 1,6-anhydro-2,3-dideoxy-2,3-difluoro-4-O-triflate-β-ᴅ-talopyranose using Et3N·3HF or in situ generated Et3N·1HF. The influence of the fluorine atom at C2 on reactivity at C4 could contribute to a new fluorine effect in nucleophilic substitution. Finally, with the continuous objective of synthesizing novel multi-vicinal fluorosugars, we prepared one difluorinated and one trifluorinated alditol analogue.

Scalar products and density matrix elements of closed-shell pair geminal wavefunctions are evaluated directly in terms of the pair amplitudes, resulting in an analog of Wick’s theorem for fermions or bosons. This expression is, in general, intractable, but it is shown how it becomes feasible in three distinct ways for Richardson–Gaudin (RG) states, the antisymmetrized geminal power, and the antisymmetrized product of strongly orthogonal geminals. Dissociation curves for hydrogen chains are computed with off-shell RG states and the antisymmetrized product of interacting geminals. Both are near exact, suggesting that the incorrect results observed with ground state RG states (a local maximum rather than smooth dissociation) may be fixable using a different RG state.

Eigenvectors of the reduced Bardeen–Cooper–Schrieffer (BCS) Hamiltonian, Richardson–Gaudin (RG) states, are used as a variational wavefunction ansatz for strongly correlated electronic systems. These states are geminal products whose coefficients are solutions of non-linear equations. Previous results showed an un-physical apparent avoided crossing in ground state dissociation curves for hydrogen chains. In this paper, it is shown that each seniority-zero state of the molecular Coulomb Hamiltonian corresponds directly to an RG state. However, the seniority-zero ground state does not correspond to the ground state of a reduced BCS Hamiltonian. The difficulty is in choosing the correct RG state. The systems studied showed a clear choice, and we expect that it should always be possible to reason physically which state to choose.