Search

Wavefunction forms based on products of electron pairs are usually constructed as closed-shell singlets, which is insufficient when the molecular state has a nonzero spin or when the chemistry is determined by $d$- or $f-$electrons. A set of two-electron forms are considered as explicit couplings of second-quantized operators to open-shell singlets. Geminal wavefunctions are constructed and their structure is elaborated. Numerical results for small model systems clearly demonstrate improvement over closed-shell singlet pairs.

We extend ground-state density-functional theory to excited states and provide the theoretical formulation for the widely used $\Delta SCF$ method for calculating excited-state energies and densities. As the electron density alone is insufficient to characterize excited states, we formulate excited-state theory using the defining variables of a noninteracting reference system, namely (1) the excitation quantum number $n_{s}$ and the potential $w_{s}(\mathbf{r})$ (excited-state potential-functional theory, $n$PFT), (2) the noninteracting wavefunction $\Phi$ ($\Phi$-functional theory, $\Phi$FT), or (3) the noninteracting one-electron reduced density matrix $\gamma_{s}(\mathbf{r},\mathbf{r}')$ (density-matrix-functional theory, $\gamma_{s}$FT). We show the equivalence of these three sets of variables and their corresponding energy functionals. Importantly, the ground and excited-state exchange-correlation energy use the \textit{same} universal functional, regardless of whether $\left(n_{s},w_{s}(\boldsymbol{r})\right)$, $\Phi$, or $\gamma_{s}(\mathbf{r},\mathbf{r}')$ is selected as the fundamental descriptor of the system. We derive the excited-state (generalized) Kohn-Sham equations. The minimum of all three functionals is the ground-state energy and, for ground states, they are all equivalent to the Hohenberg-Kohn-Sham method. The other stationary points of the functionals provide the excited-state energies and electron densities, establishing the foundation for the $\Delta SCF$ method.

We propose a new Perturbation Theory framework that can be used to help with the projective solution of the Schr\"odinger equation for arbitrary wavefunctions. This Flexible Ansatz for N-body Perturbation Theory (FANPT) is based on our previously proposed Flexible Ansatz for 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 non-linear 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 (pCCD) ansatz. We present a new hierarchy of two-electron wavefunctions that extend 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.

We provide an inductive proof of Borchardt's theorem for calculating the permanent of a Cauchy matrix via the determinants of auxiliary matrices. This result has implications for antisymmetric products of interacting geminals (APIG), and suggests that the restriction of the APIG coefficients to Cauchy form (typically called APr2G) is special in its tractability.

Direct approaches to the quantum many-body problem suffer from the so-called "curse of dimensionality": the number of parameters needed to fully specify the exact wavefunction grows exponentially with increasing system size. This motivates the develop of accurate, but approximate, ways to parametrize the wavefunction, including methods like couple cluster theory and correlator product states (CPS). Recently, there has been interest in approaches based on machine learning both direct applications of neural network architecture and the combinations of conventional wavefunction parametrizations with various Boltzmann machines. While all these methods can be exact in principle, they are usually applied with only a polynomial number of parameters, limiting their applicability. This research's objective is to present a fresh approach to wavefunction parametrization that is size-consistent, rapidly convergent, and robust numerically. Specifically, we propose a hierarchical ansatz that converges rapidly (with respect to the number of least-squares optimization). The general utility of this approach is verified by applying it to uncorrelated, weakly-correlated, and strongly-correlated systems, including small molecules and the one-dimensional Hubbard model., Comment: 18 pages

We develop a bivariational principle for an antisymmetric product of nonorthogonal geminals. Special cases reduce to the antisymmetric product of strongly-orthogonal geminals (APSG), the generalized valence bond-perfect pairing (GVB-PP), and the antisymmetrized geminal power (AGP) wavefunctions. The presented method employs wavefunctions of the same type as Richardson-Gaudin (RG) states, but which are not eigenvectors of a model Hamiltonian which would allow for more freedom in the mean-field. The general idea is to work with the same state in a primal picture in terms of pairs, and in a dual picture in terms of pair-holes. This leads to an asymmetric energy expression which may be optimized bivariationally, and is strictly variational when the two representations are consistent. The general approach may be useful in other contexts, such as for computationally feasible variational coupled-cluster methods.

Fanpy is a free and open-source Python library for developing and testing multideterminant wavefunctions and related ab initio methods in electronic structure theory. The main use of Fanpy is to quickly prototype new methods by making it easier to transfer the mathematical conception of a new wavefunction ans\"{a}tze to a working implementation. Fanpy uses the framework of our recently introduced Flexible Ansatz for N-electron Configuration Interaction (FANCI), where multideterminant wavefunctions are represented by their overlaps with Slater determinants of orthonormal spin-orbitals. In the simplest case, a new wavefunction ansatz can be implemented by simply writing a function for evaluating its overlap with an arbitrary Slater determinant. Fanpy is modular in both implementation and theory: the wavefunction model, the system's Hamiltonian, and the choice of objective function are all independent modules. This modular structure makes it easy for users to mix and match different methods and for developers to quickly try new ideas. Fanpy is written purely in Python with standard dependencies, making it accessible for most operating systems; it adheres to principles of modern software development, including comprehensive documentation, extensive testing, and continuous integration and delivery protocols. This article is considered to be the official release notes for the Fanpy library.

We present a Flexible Ansatz for N-body Configuration Interaction (FANCI) that includes any multideterminant wavefunction. This ansatz is a generalization of the Configuration Interaction (CI) wavefunction, where the coefficients are replaced by a specified function of certain parameters. By making an appropriate choice for this function, we can reproduce popular wavefunction structures like CI, Coupled-Cluster, Tensor Network States, and geminal-product wavefunctions. The universality of this framework suggests a programming structure that allows for the easy construction and optimization of arbitrary wavefunctions. Here, we will discuss the structures of the FANCI framework and its implications for wavefunction properties, particularly accuracy, cost, and size-consistency. We demonstrate the flexibility of this framework by reconstructing popular wavefunction ans\"{a}tze and modifying them to construct novel wavefunction forms. FANCI provides a powerful framework for exploring, developing, and testing new wavefunction forms.

Hydrogen cyanide (HCN) and formaldehyde (H2CO) are key precursors to biomolecules such as nucleobases and amino acids in planetary atmospheres; However, many reactions which produce and destroy these species in atmospheres containing CO2 and H2O are still missing from the literature. We use a quantum chemistry approach to find these missing reactions and calculate their rate coefficients using canonical variational transition state theory and Rice-Ramsperger-Kassel-Marcus/master equation theory at the BHandHLYP/aug-cc-pVDZ level of theory. We calculate the rate coefficients for 126 total reactions, and validate our calculations by comparing with experimental data in the 39% of available cases. Our calculated rate coefficients are most frequently within an factor of 2 of experimental values, and generally always within an order of magnitude of these values. We discover 45 previously unknown reactions, and identify 6 from this list that are most likely to dominate H2CO and HCN production and destruction in planetary atmospheres. We highlight $^1$O + CH3 $\rightarrow$ H2CO + H as a new key source, and H2CO + $^1$O $\rightarrow$ HCO + OH as a new key sink, for H2CO in upper planetary atmospheres. In this effort, we develop an oxygen extension to our consistent reduced atmospheric hybrid chemical network (CRAHCN-O), building off our previously developed network for HCN production in N2-, CH4-, and H2-dominated atmospheres (CRAHCN). This extension can be used to simulate both HCN and H2CO production in atmospheres dominated by any of CO2, N2, H2O, CH4, and H2., Comment: Accepted for publication in J Phys Chem A, 21 pages, 4 tables, 35 supplementary pages

HCN is a key ingredient for synthesizing biomolecules such as nucleobases and amino acids. We calculate 42 reaction rate coefficients directly involved with or in competition with the production of HCN in the early Earth or Titan atmospheres. These reactions are driven by methane and nitrogen radicals produced via UV photodissociation or lightning. For every reaction in this network, we calculate rate coefficients at 298 K using canonical variational transition state theory (CVT) paired with computational quantum chemistry simulations at the BHandHLYP/augcc-pVDZ level of theory. We also calculate the temperature dependence of the rate coefficients for the reactions that have barriers from 50 to 400 K. We present 15 new reaction rate coefficients with no previously known value; 93% of our calculated coefficients are within an order of magnitude of the nearest experimental or recommended values. Above 320 K, the rate coefficient for the new reaction H2CN -> HCN + H dominates. Contrary to experiments, we find the HCN reaction pathway, N + CH3 -> HCN + H2, to be inefficient and suggest that the experimental rate coefficient actually corresponds to an indirect pathway, through the H2CN intermediate. We present CVT using energies computed with density functional theory as a feasible and accurate method for calculating a large network of rate coefficients of small-molecule reactions., Comment: 34 pages, 8 figures, 14 tables, accepted for publication in J Phys Chem A

Background: The nuclear many-body system is a strongly correlated quantum system, posing serious challenges for perturbative approaches starting from uncorrelated reference states. The last decade has witnessed considerable progress in the accurate treatment of pairing correlations, one of the major components in medium-sized nuclei, reaching accuracies below the 1% level of the correlation energy. Purpose: Development of a quantum many-body method for pairing correlations that is (a) competitive in the 1% error range, and (b) can be systematically improved with a fast (exponential) convergence rate. Method: The present paper capitalizes upon ideas from Richardson-Gaudin integrability. The proposed method is a two-step approach. The first step consists of the optimization of a Richardson-Gaudin ground state as variational trial state. At the second step, the complete set of excited states on top of this Richardson-Gaudin ground state is used as an optimal basis for a Configuration Interaction method in an increasingly large effective Hilbert space. Results: The performance of the variational Richardson-Gaudin (varRG) and Richardson-Gaudin Configuration Interaction (RGCI) method is benchmarked against exact results using an effective $G$-matrix interaction for the Sn region. The varRG already reaches accuracies around the 1% level of the correlation energies, and the RGCI step sees an additional improvement scaling exponentially with the size of the effective Hilbert space. Conclusions: The Richardson-Gaudin models of integrability provide an optimized complete basis set for pairing correlations.

In methods like geminal-based approaches or coupled cluster that are solved using the projected Schr\"odinger equation, direct computation of the 2-electron reduced density matrix (2-RDM) is impractical and one falls back to a 2-RDM based on response theory. However, the 2-RDMs from response theory are not $N$-representable. That is, the response 2-RDM does not correspond to an actual physical $N$-electron wave function. We present a new algorithm for making these non-$N$-representable 2-RDMs approximately $N$-representable, i.e. it has the right symmetry and normalization and it fulfills the $P$-, $Q$- and $G$-conditions. Next to an algorithm which can be applied to any 2-RDM, we have also developed a 2-RDM optimization procedure specifically for seniority-zero 2-RDMs. We aim to find the 2-RDM with the right properties that is the closest (in the sense of the Frobenius norm) to the non-N-representable 2-RDM by minimizing the square norm of the difference between the initial 2-RDM and the targeted 2-RDM under the constraint that the trace is normalized and the 2-RDM, $Q$- and $G$-matrices are positive semidefinite, i.e. their eigenvalues are non-negative. Our method is suitable for fixing non-N-respresentable 2-RDMs which are close to being N-representable. Through the N-representability optimization algorithm we add a small correction to the initial 2-RDM such that it fulfills the most important N-representability conditions., Comment: 13 pages, 8 figures

Smooth model potentials with parameters selected to reproduce the spectrum of one-electron atoms are used to approximate the singular Coulomb potential. Even when the potentials do not mimic the Coulomb singularity, much of the spectrum is reproduced within the chemical accuracy. For the Hydrogen atom, the smooth approximations to the Coulomb potential are more accurate for higher angular momentum states. The transferability of the model potentials from an attractive interaction (Hydrogen atom) to a repulsive one (Harmonium and the uniform electron gas) is discussed.

Atomic partial charges appear in the Coulomb term of many force-field models and can be derived from electronic structure calculations with a myriad of atoms-in-molecules (AIM) methods. More advanced models have also been proposed, using the distributed nature of the electron cloud and atomic multipoles. In this work, an electrostatic force field is defined through a concise approximation of the electron density, for which the Coulomb interaction is trivially evaluated. This approximate "pro-density" is expanded in a minimal basis of atom-centered s-type Slater density functions, whose parameters are optimized by minimizing the Kullback-Leibler divergence of the pro-density from a reference electron density, e.g. obtained from an electronic structure calculation. The proposed method, Minimal Basis Iterative Stockholder (MBIS), is a variant of the Hirshfeld AIM method but it can also be used as a density-fitting technique. An iterative algorithm to refine the pro-density is easily implemented with a linear-scaling computational cost, enabling applications to supramolecular systems. The benefits of the MBIS method are demonstrated with systematic applications to molecular databases and extended models of condensed phases. A comparison to 14 other AIM methods shows its effectiveness when modeling electrostatic interactions. MBIS is also suitable for rescaling atomic polarizabilities in the Tkatchenko-Sheffler scheme for dispersion interactions., Comment: 61 pages, 12 figures, 2 tables

Actinide-containing complexes present formidable challenges for electronic structure methods due to the large number of degenerate or quasi-degenerate electronic states arising from partially occupied 5f and 6d shells. Conventional multi-reference methods can treat active spaces that are often at the upper limit of what is required for a proper treatment of species with complex electronic structures, leaving no room for verifying their suitability. In this work we address the issue of properly defining the active spaces in such calculations, and introduce a protocol to determine optimal active spaces based on the use of the Density Matrix Renormalization Group algorithm and concepts of quantum information theory. We apply the protocol to elucidate the electronic structure and bonding mechanism of volatile plutonium oxides (PuO$_3$ and PuO$_2$(OH)$_2$), species associated with nuclear safety issues for which little is known about the electronic structure and energetics. We show how, within a scalar relativistic framework, orbital-pair correlations can be used to guide the definition of optimal active spaces which provide an accurate description of static/non-dynamic electron correlation, as well as to analyse the chemical bonding beyond a simple orbital model. From this bonding analysis we are able to show that the addition of oxo- or hydroxo-groups to the plutonium dioxide species considerably changes the pi-bonding mechanism with respect to the bare triatomics, resulting in bent structures with considerable multi-reference character., Comment: 12 pages, 5 figures

We present a Linearized Coupled Cluster (LCC) correction based on an Antisymmetric Product of 1 reference orbital Geminals (AP1roG) reference state. In our LCC ansatz, the cluster operator is restricted to double and to single and double excitations as in standard single-reference CC theory. The performance of the AP1roG-LCC models is tested for the dissociation of diatomic molecules (C$_2$ and F$_2$), spectroscopic constants of the uranyl cation (UO$_2^{2+}$), and the symmetric dissociation of the H$_{50}$ hydrogen chain. Our study indicates that an LCC correction based on an AP1roG reference function is more robust and reliable than corrections based on perturbation theory, yielding spectroscopic constants that are in very good agreement with theoretical reference data., Comment: 9 pages, 4 figures