Seniority Number in Valence Bond Theory

Seniority and Hierarchy Configuration Interaction for Radicals and Excited States

Hierarchy configuration interaction (hCI) has recently been introduced as an alternative configuration interaction (CI) route combining excitation degree and seniority number and has been shown to efficiently recover both dynamic and static correlations for closed-shell molecular systems [ J. Phys. Chem. Lett. 2022, 13, 4342]. Here we generalize hCI for an arbitrary reference determinant, allowing calculations for radicals and excited states in a state-specific way. We gauge this route against excitation-based CI (eCI) and seniority-based CI (sCI) by evaluating how different ground-state properties of radicals converge to the full CI limit. We find that hCI outperforms or matches eCI, whereas sCI is far less accurate, in line with previous observations for closed-shell molecules. Employing second-order Epstein-Nesbet (EN2) perturbation theory as a correction significantly accelerates the convergence of hCI and eCI. We further explore various hCI and sCI models to calculate the excitation energies of closed- and open-shell systems. Our results underline that the choice of both the reference determinant and the set of orbitals drives the fine balance between correlation of ground and excited states. State-specific hCI2 and higher-order models perform similarly to their eCI counterparts, whereas lower orders of hCI deliver poor results unless supplemented by the EN2 correction, which substantially improves their accuracy. In turn, sCI1 produces decent excitation energies for radicals, encouraging the development of related seniority-based coupled-cluster methods.

Toward Complementary Characterization of the Chemical Bond

A precise discussion of a single bond requires consideration of two-particle wave function for the particles involved. Here we define and determine rigorously the intrinsic covalency and connected characteristics of the canonical example of the H₂ molecule. This is achieved by starting from an analytic form for the two-particle wave function for electrons forming the bond, in which we single out the atomic contribution (atomicity) in an unequivocal manner. The presence of the atomicity and ionicity factors complements the existing attributes of the bond. In this way, a gradual evolution of the molecular state to its two-atom correspondent is traced systematically with increasing interatomic distance. In effect, a direct relation to the onset of incipient Mott-Hubbard atomicity (Mottness) to the intrinsic covalency and ionicity is established. This goal is achieved formally by combining the single-particle wave function readjustment in the entangled state with a simultaneous determination of two-particle states in the particle (second quantization) representation.

Hierarchy Configuration Interaction: Combining Seniority Number and Excitation Degree

We propose a novel partitioning of the Hilbert space, hierarchy configuration interaction (hCI), where the excitation degree (with respect to a given reference determinant) and the seniority number (i.e., the number of unpaired electrons) are combined in a single hierarchy parameter. The key appealing feature of hCI is that each hierarchy level accounts for all classes of determinants whose number shares the same scaling with system size. By surveying the dissociation of multiple molecular systems, we found that the overall performance of hCI usually exceeds or, at least, parallels that of excitation-based CI. For higher orders of hCI and excitation-based CI, the additional computational burden related to orbital optimization usually does not compensate the marginal improvements compared with results obtained with Hartree-Fock orbitals. The exception is orbital-optimized CI with single excitations, a minimally correlated model displaying the qualitatively correct description of single bond breaking at a very modest computational cost.

N-Body Reduced Density Matrix-Based Valence Bond Theory and Its Applications in Diabatic Electronic-Structure Computations

Valence bond (VB) theory, as a helpful complement to the more popular molecular orbital theory, is a fundamental electronic-structure theory that aims at interpreting molecular structure and chemical reactions in a lucid way. Both theoretical and experimental chemists have shown great interest in VB theory because of its capability of providing intuitive insight into the nature of chemical bonding and the mechanism of chemical reaction in a clear and comprehensible language rooted in Lewis structure. Therefore, there is a great call for the renaissance of VB theory. Nevertheless, this is possible only after a series of methods and algorithms were developed and efficiently implemented in user-friendly programs so as to serve computational chemists for general applications. In the past three decades, we have devoted a great amount of scientific enthusiasm toward this goal. In this Account, we will concisely summarize and briefly but insightfully discuss recent developments in ab initio VB theory, especially the N-body reduced density matrices (RDM)-based approach and its applications in diabatic electronic-structure computations, which is very useful for the vivid interpretation of many fundamental chemical processes such as electron and energy transfers. Furthermore, because of the fundamentally important role that the diabatic state plays in electron and energy transfers, which are two frontier research topics in both molecular and biochemical sciences, there are a broad range of applications that VB theory can handle.We start by briefly reviewing the general feature of ab initio VB wave functions. In particular, we focus on the multistructural ab initio VB theory that uses strictly localized orbitals, including the fundamental VB self-consistent field (VBSCF) and two post-SCF methods, VBCI and VBPT2, that use the VBSCF wave function as reference. We then allot a section to describing the recent developments of the RDM-based VB approach in the second quantization language. In this section, the enhanced Wick theorem is first outlined, followed by a brief discussion of its applications in evaluating VBSCF energy gradients and a Hessian with respect to the orbital expansion coefficients, together with a short review of the implementation of an automatic formula and code generator (AFCG) designed for many-body methods with nonorthogonal orbitals. Then, we introduce the application of the RDM-based approach in implementing the post-SCF method that addresses dynamic electronic correlation via perturbation theory, viz., the icVBPT2 method that adopts an internal contraction technique naturally. We finish this section by incorporating VB theory with the concept of seniority number, in which the tensor analysis technique is carefully exploited with the RDM-based approach, resulting in significant improvements in both the number of the active electrons/orbitals and in the speedup of the computational efficiency, thus pushing VB theory to its new limit. With these achievements available, we present the applications of VB theory in diabatic electronic-structure computations by using the intuitive insight rendered by VB theory. Therefore, we believe that there is a bright future in VB theory with true opportunities and new challenges coexisting both for theoretical developments and computational applications.

Ab initio valence bond theory: A brief history, recent developments, and near future

This Perspective presents a survey of several issues in ab initio valence bond (VB) theory with a primary focus on recent advances made by the Xiamen VB group, including a brief review of the earlier history of the ab initio VB methods, in-depth discussion of algorithms for nonorthogonal orbital optimization in the VB self-consistent field method and VB methods incorporating dynamic electron correlation, along with a concise overview of VB methods for complex systems and VB models for chemical bonding and reactivity, and an outlook of opportunities and challenges for the near future of the VB theory.

Fast transformations between configuration state function and Slater determinant bases for direct configuration interaction

A hybrid configuration state function (CSF) and Slater determinant (SD) basis full configuration interaction (CI) program was developed to simultaneously take advantage of fast SD basis algorithms for σ = Hc formation and the smaller CI vector length and more robust convergence offered by a CSF basis. Graphical processing unit acceleration of the direct CSF-SD and SD-CSF basis transformation algorithms ensures that the combined transformation time per iteration relative to σ formation is small (∼15%). In addition to the obvious benefits of reducing the memory footprint of the CI vector, additional computational savings are demonstrated that rely directly on the size of the CI basis, in one particular case reducing the CI time-to-solution of a HF-CAS-(16,16)-CI/6-31G calculation of ethylene from 1954.79 s to 956 s by using a CSF basis, a 2.0× speedup.

Quantum orbital-optimized unitary coupled cluster methods in the strongly correlated regime: Can quantum algorithms outperform their classical equivalents?

The Coupled Cluster (CC) method is used to compute the electronic correlation energy in atoms and molecules and often leads to highly accurate results. However, due to its single-reference nature, standard CC in its projected form fails to describe quantum states characterized by strong electronic correlations and multi-reference projective methods become necessary. On the other hand, quantum algorithms for the solution of many-electron problems have also emerged recently. The quantum unitary variant of CC (UCC) with singles and doubles (q-UCCSD) is a popular wavefunction Ansatz for the variational quantum eigensolver algorithm. The variational nature of this approach can lead to significant advantages compared to its classical equivalent in the projected form, in particular, for the description of strong electronic correlation. However, due to the large number of gate operations required in q-UCCSD, approximations need to be introduced in order to make this approach implementable in a state-of-the-art quantum computer. In this work, we evaluate several variants of the standard q-UCCSD Ansatz in which only a subset of excitations is included. In particular, we investigate the singlet and pair q-UCCD approaches combined with orbital optimization. We show that these approaches can capture the dissociation/distortion profiles of challenging systems, such as H₄, H₂O, and N₂ molecules, as well as the one-dimensional periodic Fermi-Hubbard chain. These results promote the future use of q-UCC methods for the solution of challenging electronic structure problems in quantum chemistry.