Recent developments in the PySCF program package

Constructing Accurate and Efficient General-Purpose Atomistic Machine Learning Model with Transferable Accuracy for Quantum Chemistry

Density functional theory (DFT) has been a cornerstone in computational science, providing powerful insights into structure-property relationships for molecules and materials through first-principles quantum-mechanical (QM) calculations. However, the advent of atomistic machine learning (ML) is reshaping the landscape by enabling large-scale dynamics simulations and high-throughput screening at DFT-equivalent accuracy with drastically reduced computational cost. Yet, the development of general-purpose atomistic ML models as surrogates for QM calculations faces several challenges, particularly in terms of model capacity, data efficiency, and transferability across chemically diverse systems. This work introduces a novel extension of the polarizable atom interaction neural network (namely, XPaiNN) to address these challenges. Two distinct training strategies have been employed, one direct-learning and the other Δ-ML on top of a semiempirical QM method. These methodologies have been implemented within the same framework, allowing for a detailed comparison of their results. The XPaiNN models, in particular the one using Δ-ML, not only demonstrate competitive performance on standard benchmarks, but also demonstrate the effectiveness against other ML models and QM methods on comprehensive downstream tasks, including noncovalent interactions, reaction energetics, barrier heights, geometry optimization and reaction thermodynamics, etc. This work represents a significant step forward in the pursuit of accurate and efficient atomistic ML models of general-purpose, capable of handling complex chemical systems with transferable accuracy.

Fast and accurate nonadiabatic molecular dynamics enabled through variational interpolation of correlated electron wavefunctions

We build on the concept of eigenvector continuation to develop an efficient multi-state method for the rigorous and smooth interpolation of a small training set of many-body wavefunctions through chemical space at mean-field cost. The inferred states are represented as variationally optimal linear combinations of the training states transferred between the many-body bases of different nuclear geometries. We show that analytic multi-state forces and nonadiabatic couplings from the model enable application to nonadiabatic molecular dynamics, developing an active learning scheme to ensure a compact and systematically improvable training set. This culminates in application to the nonadiabatic molecular dynamics of a photoexcited 28-atom hydrogen chain, with surprising complexity in the resulting nuclear motion. With just 22 DMRG calculations of training states from the low-energy correlated electronic structure at different geometries, we infer the multi-state energies, forces and nonadiabatic coupling vectors at 12 000 geometries with provable convergence to high accuracy along an ensemble of molecular trajectories, which would not be feasible with a brute force approach. This opens up a route to bridge the timescales between accurate single-point correlated electronic structure methods and timescales of relevance for photo-induced molecular dynamics.

Quantum embedding for molecules using auxiliary particles - the ghost Gutzwiller Ansatz

Strong/static electronic correlation mediates the emergence of remarkable phases of matter, and underlies the exceptional reactivity properties in transition metal-based catalysts. Modeling strongly correlated molecules and solids calls for multi-reference Ansätze, which explicitly capture the competition of energy scales characteristic of such systems. With the efficient computational screening of correlated solids in mind, the ghost Gutzwiller (gGut) Ansatz has been recently developed. This is a variational Ansatz which can be formulated as a self-consistent embedding approach, describing the system within a non-interacting, quasiparticle model, yet providing accurate spectra in both low and high energy regimes. Crucially, small fragments of the system are identified as responsible for the strong correlation, and are therefore enhanced by adding a set of auxiliary orbitals, the ghosts. These capture many-body correlations through one-body fluctuations and subsequent out-projection when computing physical observables. gGut has been shown to accurately describe multi-orbital lattice models at modest computational cost. In this work, we extend the gGut framework to strongly correlated molecules, for which it holds special promise. Indeed, despite the asymmetric embedding treatment, the quasiparticle Hamiltonian effectively describes all major sources of correlation in the molecule: strong correlation through the ghosts in the fragment, and dynamical correlation through the quasiparticle description of its environment. To adapt the gGut Ansatz for molecules, we address the fact that, unlike in the lattice model previously considered, electronic interactions in molecules are not local. Hence, we explore a hierarchy of approximations of increasing accuracy capturing interactions between fragments and environment, and within the environment, and discuss how these affect the embedding description of correlations in the whole molecule. We will compare the accuracy of the gGut model with established methods to capture strong correlation within active space formulations, and assess the realistic use of this novel approximation to the theoretical description of correlated molecular clusters.

On the notion of strong correlation in electronic structure theory

Strong correlation has been said to have many faces, and appears to have many synonyms of questionable suitability. In this work we aim not to define the term once and for all, but to highlight one possibility that is both rigorously defined and physically transparent, and remains so in reference to molecules and quantum lattice models. We survey both molecular examples - hydrogen systems (Hn, n = 2, 4, 6), Be2, H-He-H, and benzene - and the half-filled Hubbard model over a range of correlation regimes. Various quantities are examined including the extent of spin symmetry breaking in correlated single-reference wave functions, energetic ratios inspired by the Hubbard model and the Virial theorem, and metrics derived from the one- and two-electron reduced density matrices (RDMs). The trace and the square norm of the cumulant of the two-electron reduced density matrix capture what may well be defined as strong correlation. Accordingly, strong correlation is understood as a statistical dependence between two electrons, and is distinct from the concepts of "correlation energy" and more general than entanglement quantities that require a partitioning of a quantum system into distinguishable subspaces. This work enables us to build a bridge between a rigorous and quantifiable regime of strong electron correlation and more familiar chemical concepts such as anti-aromaticity in the context of Baird's rule.

Challenges with relativistic GW calculations in solids and molecules

For molecules and solids containing heavy elements, accurate electronic-structure calculations require accounting not only for electronic correlations but also for relativistic effects. In molecules, relativity can lead to severe changes in the ground-state description. In solids, the interplay between both correlation and relativity can change the stability of phases or it can lead to an emergence of completely new phases. Traditionally, the simplest illustration of relativistic effects can be done either by including pseudopotentials in non-relativistic calculations or alternatively by employing large all-electron basis sets in relativistic methods. By analyzing different electronic properties (band structure, equilibrium lattice constant and bulk modulus) in semiconductors and insulators, we show that capturing the interplay of relativity and electron correlation can be rather challenging in Green's function methods. For molecular problems with heavy elements, we also observe that similar problems persist. We trace these challenges to three major problems: deficiencies in pseudopotential treatment as applied to Green's function methods, the scarcity of accurate and compact all-electron basis sets that can be converged with respect to the basis-set size, and linear dependencies arising in all-electron basis sets, particularly when employing Gaussian orbitals. Our analysis provides detailed insight into these problems and opens a discussion about potential approaches to mitigate them.

Restoring translational symmetry in periodic all-orbital dynamical mean-field theory simulations

Dynamical mean-field theory (DMFT) and its cluster extensions provide an efficient Green's function formalism to simulate spectral properties of periodic systems at the quantum many-body level. However, traditional cluster DMFT breaks translational invariance in solid-state materials, and the best strategy to capture non-local correlation effects within cluster DMFT remains elusive. In this work, we investigate the use of overlapping atom-centered impurity fragments in recently-developed ab initio all-orbital DMFT, where all local orbitals within the impurity are treated with high-level quantum chemistry impurity solvers. We demonstrate how the translational symmetry of the lattice self-energy can be restored by designing symmetry-adapted embedding problems, which results in an improved description of spectral functions in two-dimensional boron nitride monolayers and graphene at the levels of many-body perturbation theory (GW) and coupled-cluster theory. Furthermore, we study the convergence of self-energy and density of states as the embedding size is systematically expanded in one-shot and self-consistent DMFT calculations.

Accelerated basis-set convergence of coupled-cluster excitation energies using the density-based basis-set correction method

We present the first application to real molecular systems of the recently proposed linear-response theory for the density-based basis-set correction method [J. Chem. Phys., 158, 234107 (2023)]. We apply this approach to accelerate the basis-set convergence of excitation energies in the equation-of-motion coupled-cluster singles and doubles (EOM-CCSD) method. We use an approximate linear-response framework that neglects the second-order derivative of the basis-set correction density functional and consists in simply adding to the usual Hamiltonian the one-electron potential generated by the first-order derivative of the functional. This additional basis-set correction potential is evaluated at the Hartree-Fock density, leading to a very computationally cheap basis-set correction. We tested this approach over a set of about 30 excitation energies computed for five small molecular systems and found that the excitation energies from the ground state to Rydberg states are the main source of basis-set error. These excitation energies systematically increase when the size of the basis set is increased, suggesting a biased description in favour of the excited state. Despite the simplicity of the present approach, the results obtained with the basis-set-corrected EOM-CCSD method are encouraging as they yield a mean absolute deviation of 0.02 eV for the aug-cc-pVTZ basis set, while it is 0.04 eV using the standard EOM-CCSD method. This might open the path to an alternative to explicitly correlated approaches to accelerate the basis-set convergence of excitation energies.

Towards efficient quantum computing for quantum chemistry: reducing circuit complexity with transcorrelated and adaptive ansatz techniques

The near-term utility of quantum computers is hindered by hardware constraints in the form of noise. One path to achieving noise resilience in hybrid quantum algorithms is to decrease the required circuit depth - the number of applied gates - to solve a given problem. This work demonstrates how to reduce circuit depth by combining the transcorrelated (TC) approach with adaptive quantum ansätze and their implementations in the context of variational quantum imaginary time evolution (AVQITE). The combined TC-AVQITE method is used to calculate ground state energies across the potential energy surfaces of H4, LiH, and H2O. In particular, H4 is a notoriously difficult case where unitary coupled cluster theory, including singles and doubles excitations, fails to provide accurate results. Adding TC yields energies close to the complete basis set (CBS) limit while reducing the number of necessary operators - and thus circuit depth - in the adaptive ansätze. The reduced circuit depth furthermore makes our algorithm more noise-resilient and accelerates convergence. Our study demonstrates that combining the TC method with adaptive ansätze yields compact, noise-resilient, and easy-to-optimize quantum circuits that yield accurate quantum chemistry results close to the CBS limit.

Adsorption and vibrational spectroscopy of CO on the surface of MgO from periodic local coupled-cluster theory

The adsorption of CO on the surface of MgO has long been a model problem in surface chemistry. Here, we report periodic Gaussian-based calculations for this problem using second-order perturbation theory (MP2) and coupled-cluster theory with single and double excitations (CCSD) and perturbative triple excitations [CCSD(T)], with the latter two performed using a recently developed extension of the local natural orbital approximation to problems with periodic boundary conditions. The low cost of periodic local correlation calculations allows us to calculate the full CCSD(T) binding curve of CO approaching the surface of MgO (and thus the adsorption energy) and the two-dimensional potential energy surface (PES) as a function of the distance from the surface and the CO stretching coordinate. From the PES, we obtain the fundamental vibrational frequency of CO on MgO, whose shift from the gas phase value is a common experimental probe of surface adsorption. We find that CCSD(T) correctly predicts a positive frequency shift upon adsorption of +14.7 cm-1, in excellent agreement with the experimental shift of +14.3 cm-1. We use our CCSD(T) results to assess the accuracy of MP2, CCSD, and several density functional theory (DFT) approximations, including exchange correlation functionals and dispersion corrections. We find that MP2 and CCSD yield reasonable binding energies and frequency shifts, whereas many DFT calculations overestimate the magnitude of the adsorption energy by 5-15 kJ mol-1 and predict a negative frequency shift of about -20 cm-1, which we attribute to self-interaction-induced delocalization errors that are mildly ameliorated with hybrid functionals. Our findings highlight the accuracy and computational efficiency of the periodic local correlation for the simulation of surface chemistry with accurate wavefunction methods.

Accurate and interpretable representation of correlated electronic structure via Tensor Product Selected CI

The task of computing wavefunctions that are accurate, yet simple enough mathematical objects to use for reasoning, has long been a challenge in quantum chemistry. The difficulty in drawing physical conclusions from a wavefunction is often related to the generally large number of configurations with similar weights. In Tensor Product Selected Configuration Interaction (TPSCI), we use a locally correlated tensor product state basis, which has the effect of concentrating the weight of a state onto a smaller number of physically interpretable degrees of freedom. In this paper, we apply TPSCI to a series of three molecular systems ranging in separability, one of which is the first application of TPSCI to an open-shell bimetallic system. For each of these systems, we obtain accurate solutions to large active spaces, and analyze the resulting wavefunctions through a series of different approaches including (i) direct inspection of the TPS basis coefficients, (ii) construction of Bloch effective Hamiltonians, and (iii) computation of cluster correlation functions.

Cavity-modified local and non-local electronic interactions in molecular ensembles under vibrational strong coupling

Resonant vibrational strong coupling (VSC) between molecular vibrations and quantized field modes of low-frequency optical cavities constitutes the conceptual cornerstone of vibro-polaritonic chemistry. In this work, we theoretically investigate the role of complementary nonresonant electron-photon interactions in the cavity Born-Oppenheimer (CBO) approximation. In particular, we study cavity-induced modifications of local and non-local electronic interactions in dipole-coupled molecular ensembles under VSC. Methodologically, we combine CBO perturbation theory (CBO-PT) [E. W. Fischer and P. Saalfrank, J. Chem. Theory Comput. 19, 7215 (2023)] with non-perturbative CBO Hartree-Fock (HF) and coupled cluster (CC) theories. In a first step, we derive up to second-order CBO-PT cavity potential energy surfaces, which reveal non-trivial intra- and inter-molecular corrections induced by the cavity. We then introduce the concept of a cavity reaction potential (CRP), minimizing the electronic energy in the cavity subspace to discuss vibro-polaritonic reaction mechanisms. We present reformulations of CBO-HF and CBO-CC approaches for CRPs and derive second-order approximate CRPs from CBO-PT for unimolecular and bimolecular scenarios. In the unimolecular case, we find small local modifications of molecular potential energy surfaces for selected isomerization reactions dominantly captured by the first-order dipole fluctuation correction. Excellent agreement between CBO-PT and non-perturbative wave function results indicates minor VSC-induced state relaxation effects in the single-molecule limit. In the bimolecular scenario, CBO-PT reveals an explicit coupling of interacting dimers to cavity modes besides cavity-polarization dependent dipole-induced dipole and van der Waals interactions with enhanced long-range character. An illustrative CBO-coupled cluster theory with singles and doubles-based numerical analysis of selected molecular dimer models provides a complementary non-perturbative perspective on cavity-modified intermolecular interactions under VSC.

Alchemical insights into approximately quadratic energies of iso-electronic atoms

Accurate quantum mechanics based predictions of property trends are so important for material design and discovery that even inexpensive approximate methods are valuable. We use the alchemical integral transform to study multi-electron atoms and to gain a better understanding of the approximately quadratic behavior of energy differences between iso-electronic atoms in their nuclear charges. Based on this, we arrive at the following simple analytical estimate of energy differences between any two iso-electronic atoms, ΔE≈-(1+2γNe-1)ΔZZ̄. Here, γ ≈ 0.3766 ± 0.0020 Ha corresponds to an empirical constant, and Ne, ΔZ, and Z̄, respectively, to electron number, nuclear charge difference, and average. We compare the formula's predictive accuracy using experimental numbers and non-relativistic, numerical results obtained via density functional theory (pbe0) for the entire periodic table up to Radon. A detailed discussion of the atomic helium-series is included.

A static quantum embedding scheme based on coupled cluster theory

We develop a static quantum embedding scheme that utilizes different levels of approximations to coupled cluster (CC) theory for an active fragment region and its environment. To reduce the computational cost, we solve the local fragment problem using a high-level CC method and address the environment problem with a lower-level Møller-Plesset (MP) perturbative method. This embedding approach inherits many conceptual developments from the hybrid second-order Møller-Plesset (MP2) and CC works by Nooijen [J. Chem. Phys. 111, 10815 (1999)] and Bochevarov and Sherrill [J. Chem. Phys. 122, 234110 (2005)]. We go beyond those works here by primarily targeting a specific localized fragment of a molecule and also introducing an alternative mechanism to relax the environment within this framework. We will call this approach MP-CC. We demonstrate the effectiveness of MP-CC on several potential energy curves and a set of thermochemical reaction energies, using CC with singles and doubles as the fragment solver, and MP2-like treatments of the environment. The results are substantially improved by the inclusion of orbital relaxation in the environment. Using localized bonds as the active fragment, we also report results for N=N bond breaking in azomethane and for the central C-C bond torsion in butadiene. We find that when the fragment Hilbert space size remains fixed (e.g., when determined by an intrinsic atomic orbital approach), the method achieves comparable accuracy with both a small and a large basis set. Additionally, our results indicate that increasing the fragment Hilbert space size systematically enhances the accuracy of observables, approaching the precision of the full CC solver.

Unified Implementation of Relativistic Wave Function Methods: 4C-iCIPT2 as a Showcase

In parallel to the unified construction of relativistic Hamiltonians based solely on physical arguments (J. Chem. Phys. 2024, 160, 084111), a unified implementation of relativistic wave function methods is achieved here via programming techniques (e.g., template metaprogramming and polymorphism in C++). That is, once the code for constructing the Hamiltonian matrix is made ready, all the rest can be generated automatically from existing templates used for the nonrelativistic counterparts. This is facilitated by decomposing a second-quantized relativistic Hamiltonian into diagrams that are topologically the same as those required for computing the basic coupling coefficients between spin-free configuration state functions (CSF). Moreover, both time reversal and binary double point group symmetries can readily be incorporated into molecular integrals and Hamiltonian matrix elements. The latter can first be evaluated in the space of (randomly selected) spin-dependent determinants and then transformed to that of spin-dependent CSFs, thanks to simple relations in between. As a showcase, we consider here the no-pair four-component relativistic iterative configuration interaction with selection and perturbation correction (4C-iCIPT2), which is a natural extension of the spin-free iCIPT2 (J. Chem. Theory Comput. 2021, 17, 949), and can provide near-exact numerical results within the manifold of positive energy states (PES), as demonstrated by numerical examples.

Periodic Local Coupled-Cluster Theory for Insulators and Metals

We describe the implementation details of periodic local coupled-cluster theory with single and double excitations (CCSD) and perturbative triple excitations [CCSD(T)] using local natural orbitals (LNOs) and k-point symmetry. We discuss and compare several choices for orbital localization, fragmentation, and LNO construction. By studying diamond and lithium, we demonstrate that periodic LNO-CC theory can be applied with equal success to both insulators and metals, achieving speedups of 2 to 3 orders of magnitude even for moderately sized k-point meshes. Our final predictions of the equilibrium cohesive energy, lattice constant, and bulk modulus for diamond and lithium are in good agreement with previous theoretical predictions and experimental results.

Semiclassical Nonadiabatic Molecular Dynamics Using Linearized Pair-Density Functional Theory

Nonadiabatic molecular dynamics is an effective method for modeling nonradiative decay in electronically excited molecules. Its accuracy depends strongly on the quality of the potential energy surfaces, and its affordability for long direct-dynamic simulations with adequate ensemble averaging depends strongly on the cost of the required electronic structure calculations. Linearized pair-density functional theory (L-PDFT) is a recently developed post-self-consistent-field multireference method that can model potential energy surfaces with an accuracy similar to expensive multireference perturbation theories but at a computational cost similar to the underlying multiconfiguration self-consistent field method. Here, we integrate the SHARC dynamics and PySCF electronic structure code to utilize L-PDFT for electronically nonadiabatic calculations and use the combined programs to study the photoisomerization reaction of cis-azomethane. We show that L-PDFT is able to successfully simulate the photoisomerization without crashes, and it yields results similar to the more expensive extended multistate complete active space second-order perturbation theory. This shows that L-PDFT can model internal conversion, and it demonstrates its promise for broader photodynamics applications.

Restricted Open-Shell Cluster Mean-Field theory for Strongly Correlated Systems

The cluster-based Mean Field method (cMF) and it is second order perturbative correction was introduced by Jiménez-Hoyos and Scuseria to reduce the cost of modeling strongly correlated systems by dividing an active space up into small clusters, which are individually solved in the mean-field presence of each other. In that work, clusters with unpaired electrons are treated by allowing the α and β orbitals to spin polarize. While that provided significant energetic stabilization, the resulting cMF wave function was spin-contaminated, making it difficult to use as a reference state for spin-pure post-cMF methods. In this work, we propose the Restricted Open-shell cMF (RO-cMF) method, extending the cMF approach to systems with open-shell clusters, while not permitting spin-polarization. While the resulting RO-cMF energies are necessarily higher in energy than the unrestricted orbital cMF, the new RO-cMF provides a simple reference state for post-cMF methods that recover the missing intercluster correlations. We provide a detailed explanation of the method, and report demonstrative calculations of exchange coupling constants for three systems: a di-iron complex, a dichromium complex, and a dimerized organic radical. We also report the first perturbatively corrected RO-cMF-PT2 results as well.

Heron: Visualizing and Controlling Chemical Reaction Explorations and Networks

Automated and high-throughput quantum chemical investigations into chemical processes have become feasible in great detail and broad scope. This results in an increase in complexity of the tasks and in the amount of generated data. An efficient and intuitive way for an operator to interact with these data and to steer virtual experiments is required. Here, we introduce Heron, a graphical user interface that allows for advanced human-machine interactions with quantum chemical exploration campaigns into molecular structure and reactivity. Heron offers access to interactive and automated explorations of chemical reactions with standard electronic structure modules, haptic force feedback, microkinetic modeling, and refinement of data by automated correlated calculations including black-box complete active space calculations. It is tailored to the exploration and analysis of vast chemical reaction networks. We show how interoperable modules enable advanced workflows and pave the way for routine low-entrance-barrier access to advanced modeling techniques.

Variational quantum imaginary time evolution for matrix product state Ansatz with tests on transcorrelated Hamiltonians

The matrix product state (MPS) Ansatz offers a promising approach for finding the ground state of molecular Hamiltonians and solving quantum chemistry problems. Building on this concept, the proposed technique of quantum circuit MPS (QCMPS) enables the simulation of chemical systems using a relatively small number of qubits. In this study, we enhance the optimization performance of the QCMPS Ansatz by employing the variational quantum imaginary time evolution (VarQITE) approach. Guided by McLachlan's variational principle, the VarQITE method provides analytical metrics and gradients, resulting in improved convergence efficiency and robustness of the QCMPS. We validate these improvements numerically through simulations of H2, H4, and LiH molecules. In addition, given that VarQITE is applicable to non-Hermitian Hamiltonians, we evaluate its effectiveness in preparing the ground state of transcorrelated Hamiltonians. This approach yields energy estimates comparable to the complete basis set (CBS) limit while using even fewer qubits. In particular, we perform simulations of the beryllium atom and LiH molecule using only three qubits, maintaining high fidelity with the CBS ground state energy of these systems. This qubit reduction is achieved through the combined advantages of both the QCMPS Ansatz and transcorrelation. Our findings demonstrate the potential practicality of this quantum chemistry algorithm on near-term quantum devices.

Vertical Ionization Energies, Generalized Kohn-Sham Orbital Energies, and the Curious Case of the Copper Oxide Anions

Are the vertical ionization energies from a bound electronic system, initially in its ground state, equal to minus the corresponding exact Kohn-Sham orbital energies of density functional theory (DFT)? This is known to be true for the first or lowest vertical ionization energy. We show that the correction from time-dependent DFT arises from the continuum and need not vanish. Recent work compared the experimental photoemission thresholds of the molecules Cu2O-, CuO-, CuO2-, and CuO3- with minus the corresponding orbital energies from a generalized gradient approximation (GGA) and its global and range-separated hybrids with exact exchange, finding striking differences which were attributed to self-interaction error, strong correlation, or both. Here, we extend that work to include the local spin density approximation (LSDA), its Perdew-Zunger self-interaction correction with Fermi-Löwdin localized orbitals (LSDA-SIC), a quasi-self-consistent locally scaled-down version of LSDA-SIC (QLSIC), and the Quantum Theory Project QTP02 range-separated hybrid functional, all but LSDA implemented in a generalized Kohn-Sham approach. QTP02 impressively yields a near equality for many sp-bonded molecules. However, for the copper oxide anions studied here, none of the tested methods reproduces the experimental photoemission thresholds.

Complete Active Space Iterative Coupled Cluster Theory

In this work, we investigate the possibility of improving multireference-driven coupled cluster (CC) approaches with an algorithm that iteratively combines complete active space (CAS) calculations with tailored CC and externally corrected CC. This is accomplished by establishing a feedback loop between the CC and CAS parts of a calculation through a similarity transformation of the Hamiltonian with those CC amplitudes that are not encompassed by the active space. We denote this approach as the complete active space iterative coupled cluster (CASiCC) ansatz. We investigate its efficiency and accuracy in the singles and doubles approximation by studying the prototypical molecules H4, H8, H2O, and N2. Our results demonstrate that CASiCC systematically improves on the single-reference CCSD and the externally corrected CCSD methods across entire potential energy curves while retaining modest computational costs. However, the tailored coupled cluster method shows superior performance in the strong correlation regime, suggesting that its accuracy is based on error compensation. We find that the iterative versions of externally corrected and tailored coupled cluster methods converge to the same results.

A Critical Evaluation of the Hybrid KS DFT Functionals Based on the KS Exchange-Correlation Potentials

We have developed a critical methodology for the evaluation of the quality of hybrid exchange-correlation (XC) density functional approximations (DFAs) based on very fundamental quantities, i.e., Kohn-Sham (KS) XC potentials, self-consistent electron densities, first ionization potentials (IPs), and total energies. Since the XC potentials, the primary objects in the current study, are not directly accessible for the hybrids, we calculate them by inverting the KS electron densities. Utilizing this methodology, we tested 155 hybrid DFAs available in the LIBXC library using FCI and CCSD(T) methods as a reference. We have found that a group of functionals produces very decent XC potentials, mainly those with a large mixture of Hartree-Fock exchange. Moreover, the value of IP strongly depends on the XC potential quality. On the other hand, we show that the XC energy is dominated by functional-driven error, which in some cases leads to substantial errors in electronic densities. The study shows new directions for constructing more accurate XC functionals within the KS-DFT framework.

Quantum Algorithms for the Variational Optimization of Correlated Electronic States with Stochastic Reconfiguration and the Linear Method

Solving the electronic Schrodinger equation for strongly correlated ground states is a long-standing challenge. We present quantum algorithms for the variational optimization of wave functions correlated by products of unitary operators, such as Local Unitary Cluster Jastrow (LUCJ) ansatzes, using stochastic reconfiguration (SR) and the linear method (LM). While an implementation on classical computing hardware would require exponentially growing compute cost, the cost (number of circuits and shots) of our quantum algorithms is polynomial in system size. We find that classical simulations of optimization with the linear method consistently find lower energy solutions than with the L-BFGS-B optimizer across the dissociation curves of the notoriously difficult N2 and C2 dimers; LUCJ predictions of the ground-state energies deviate from exact diagonalization by 1 kcal/mol or less at all points on the potential energy curve. While we do characterize the effect of shot noise on the LM optimization, these noiseless results highlight the critical but often overlooked role that optimization techniques must play in attacking the electronic structure problem (on both classical and quantum hardware), for which even mean-field optimization is formally NP hard. We also discuss the challenge of obtaining smooth curves in these strongly correlated regimes, and propose a number of quantum-friendly solutions ranging from symmetry-projected ansatz forms to a symmetry-constrained optimization algorithm.

Reconciling Accuracy and Feasibility for Barrierless Reaction Steps by the PCS/DDCI/MC-PDFT Protocol: Methane and Ethylene Dissociations as Case Studies

Several enhancements have been introduced into state-of-the-art computational protocols for the treatment of barrierless reaction steps in the framework of variable reaction coordinate variational transition state theory. The first step is the synergistic integration of the Iterative Difference Dedicated Configuration Interaction (I-DDCI) and Pisa Composite Scheme, which defines a reduced cost, yet very accurate, computational workflow. This approach provides a near black box tool for obtaining 1D reference potentials. Then, a general strategy has been devised for tuning the level of theory used in Monte Carlo (MC) sampling, employing Multiconfiguration Pair Density Functional Theory (MC-PDFT) with dynamically adjusted Hartree-Fock exchange. Concurrently, partial geometry optimizations during the MC simulations account for the coupling between the reaction coordinates and conserved modes. The protocol closely approaches full size consistency and yields highly accurate results, with several test computations suggesting rapid convergence of the I-DDCI correction with the basis set dimensions. The capabilities of the new platform are illustrated by two case studies (the hydrogen dissociation from CH4 and C2H4), which highlight its flexibility in handling different carbon hybridizations (sp3 and sp2). The remarkable accuracy of the computed rate constants confirms the robustness of the proposed method. Together with their intrinsic interest, these results pave the way for systematic investigations of complex gas-phase reactions through a reliable, user-friendly tool accessible to specialists and nonspecialists alike.

Configuration Interaction in Frontier Molecular Orbital Basis for Screening the Spin-Correlated, Spatially Separated Triplet Pair State 1(T···T) Formation

In the theoretical screening of Singlet Fission rates in molecular aggregates, often the frontier molecular orbital model for dimers is employed. However, the dimer approach fails to account for recent experimental findings that suggest singlet fission progresses through a further intermediate state featuring two spatially separated, spin-correlated triplets, specifically a 1(T···T) state. We address this limitation by generalizing the often used frontier molecular orbital model for singlet fission by incorporation of both separated Charge Transfer (C···T) and 1(T···T) states as well as mixed triplet-charge transfer states, delivering analytic expressions for the diabatic matrix elements. Applying the methodology to the perylene diimide trimer, we examine the packing dependence of competing formation pathways of the 1(T···T) state by evaluation of diabatic matrix elements.

Exploring Ground and Excited States Via Single Reference Coupled-Cluster Theory and Algebraic Geometry

The exploration of the root structure of coupled cluster (CC) equations holds both foundational and practical significance for computational quantum chemistry. This study provides insight into the intricate root structures of these nonlinear equations at both the CCD and CCSD level of theory. We utilize computational techniques from algebraic geometry, specifically the monodromy and parametric homotopy continuation methods, to calculate the full solution set. We compare the computed CC roots against various established theoretical upper bounds, shedding light on the accuracy and efficiency of these bounds. We hereby focus on the dissociation processes of four-electron systems such as (H2)2 in both D2h and D∞h configurations, H4 symmetrically distorted on a circle, and lithium hydride. We moreover investigate the ability of single-reference CC solutions to approximate excited state energies. We find that multiple CC roots describe energies of excited states with high accuracy. Our investigations reveal that for systems like lithium hydride, CC not only provides high-accuracy approximations to several excited state energies but also to the states themselves.

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.

Quantum Information Driven Ansatz (QIDA): Shallow-Depth Empirical Quantum Circuits from Quantum Chemistry

Hardware-efficient empirical variational ansätze for Variational Quantum Eigensolver (VQE) simulations of quantum chemistry often lack a direct connection to classical quantum chemistry methods. In this work, we propose a method to bridge this gap by introducing a novel approach to constructing a starting point for variational quantum circuits, leveraging quantum mutual information from classical quantum chemistry states to design simple yet effective heuristic ansätze with a topology reflecting the molecular system's correlations. As a first step, we make use of quantum chemistry calculations, such as Mo̷ller-Plesset (MP2) perturbation theory, to initially provide approximate Natural Orbitals, which have been shown to be the best candidate one-electron basis for developing compact empirical wave functions.1 Second, we evaluate the quantum mutual information matrix, which provides insights about the main correlations between qubits of the quantum circuit, and enables a direct design of entangling blocks for the circuit. The resulting ansatz is then used with a VQE to obtain a short-depth variational ground state of the electronic Hamiltonian. To validate our approach, we perform a comprehensive statistical analysis through simulations of various molecular systems (H2, LiH, H2O) and apply it to the more complex NH3 molecule. The reported results demonstrate that the proposed methodology gives rise to highly effective ansätze, outperforming the standard empirical ladder-entangler ansatz. Overall, our approach can be used as an effective state preparation, providing a promising route for designing efficient variational quantum circuits for large molecular systems.

Wannier Function Localization Using Bloch Intrinsic Atomic Orbitals

We extend the intrinsic atomic orbital (IAO) method for the localization of molecular orbitals to calculate well-localized generalized Wannier functions in crystals in the spirit of the Pipek-Mezey method. We furthermore present a one-shot diabatic Wannierization procedure that aligns the phases of the Bloch functions, providing immediate Wannier localization, which serves as an excellent initial guess for optimization. We test our Wannier localization implementation on a number of solid-state systems, highlighting the effectiveness of the diabatic preparation, especially for localizing core bands. Partial charges of Wannier functions generated using Bloch IAOs align well with chemical intuition, which we demonstrate through the example of the adsorption of CO on a MgO surface.

Simulation of a Diels-Alder reaction on a quantum computer

The simulation of chemical reactions is an anticipated application of quantum computers. Using a Diels-Alder reaction as a test case, in this study we explore the potential applications of quantum algorithms and hardware in investigating chemical reactions. Our specific goal is to calculate the activation barrier of a reaction between ethylene and cyclopentadiene forming a transition state. To achieve this goal, we use quantum algorithms for near-term quantum hardware (entanglement forging and quantum subspace expansion) and classical post-processing (many-body perturbation theory) in concert. We conduct simulations on IBM quantum hardware using up to 8 qubits, and compute accurate activation barrier in the reaction between cyclopentadiene and ethylene by accounting for both static and dynamic electronic correlation. This work illustrates a hybrid quantum-classical computational workflow to study chemical reactions on near-term quantum devices, showcasing the potential for performing quantum chemistry simulations on quantum hardware to predict activation barriers in agreement with those predicted by CASCI.

Triple Electron Attachments with a New Intermediate-Hamiltonian Fock-Space Coupled-Cluster Method

The implementation of a new intermediate-Hamiltonian Fock-space coupled-cluster (IHFSCC) scheme for the (3,0) sector of the Fock space is reported. In this IHFSCC approach, the three-body contributions in the cluster operator S(3,0) corresponding to the (3,0) sector of the Fock space are considered, while S(1,0) and S(2,0) at the (1,0) and (2,0) level only include one- and two-body operators. By introducing a suitable partition of the wave operator, the intermediate Hamiltonian, which only depends on the two-body operator of S(1,0), is obtained. S(2,0) and S(3,0) are not required within this new IHFSCC scheme, and a large reference space can be possibly employed. The performance of this (3,0) IHFSCC method in calculating triple ionization potentials and excitation energies for atoms and cations with a ground p3 configuration as well as adiabatic excitation energies for some molecules is investigated. The effect of the number of active virtual orbitals and three different types of orbitals, i.e., reference orbitals, restricted open-shell Hartree-Fock orbitals (ROHF) of the target systems, and canonicalized ROHF orbitals, on IHFSCC results, is also studied. Our calculations indicate that reasonable results can be achieved with this (3,0) IHFSCC method when a minimal reference space is employed. Further increasing the number of active orbitals does not necessarily improve the results. In addition, the IHFSCC results using canonicalized ROHF orbitals generally agree well with reference values, and they are not very sensitive to the number of active orbitals compared with results using the reference orbitals. The new (3,0) IHFSCC method can be applied to open-shell systems with three unpaired electrons with reasonable accuracy at a relatively low computational cost.

A DFT/MRCI Hamiltonian parameterized using only ab initio data. II. Core-excited states

A newly parameterized combined density functional theory and multi-reference configuration interaction (DFT/MRCI) Hamiltonian, termed core-valence separation (CVS)-QE12, is defined for the computation of K-shell core-excitation and core-ionization energies. This CVS counterpart to the recently reported QE8 Hamiltonian [Costain et al., J. Chem. Phys, 160, 224106 (2024)] is parameterized by fitting to benchmark quality ab initio data. The definition of the CVS-QE12 and QE8 Hamiltonians differ from previous CVS-DFT/MRCI parameterizations in three primary ways: (i) the replacement of the BHLYP exchange-correlation functional with QTP17 to yield a balanced description of both core and valence excitation energies, (ii) the adoption of a new, three-parameter damping function, and (iii) the introduction of separate scaling of the core-valence and valence-valence Coulombic interactions. Crucially, the parameters of the CVS-QE12 Hamiltonian are obtained via fitting exclusively to highly accurate ab initio vertical core-excitation and ionization energies computed at the CVS-EOM-CCSDT level of theory. The CVS-QE12 Hamiltonian is validated against further benchmark computations and is found to furnish K-edge core vertical excitation and ionization energies exhibiting absolute errors ≤0.5 eV at low computational cost.

Modulating molecular plasmons in naphthalene via intermolecular interactions and strong light-matter coupling

We conducted a theoretical investigation on the modulation of plasmon-like resonances in naphthalene - the so-called molecular plasmons - through intermolecular interactions and strong light-matter coupling. The configuration interaction with single excitations (CIS) approach and its quantum electrodynamics extension (QED-CIS-1) are used to describe the molecular plasmon states under these interactions. We detail the effects of changing intermolecular distances of the naphthalene dimer and incorporating the naphthalene molecule into optical cavities, both allowing for precise control of naphthalene's plasmonic responses. Our results show significant shifts of the plasmon peak in the absorption spectra of naphthalene, depending on the spatial configuration of the dimer and cavity parameters such as polarization, frequency, and coupling strength. Further investigation of the naphthalene dimer in a cavity reveals a synergistic effect on the plasmon peak when the two types of interactions are combined. This research provides insights into the plasmonic behavior of simple polyacenes like naphthalene and opens up possibilities for plasmon modulation in more complex systems.

Numerical Investigation of the Quantum Inverse Algorithm on Small Molecules

We evaluate the accuracy of the quantum inverse (Q-Inv) algorithm, in which the multiplication of Ĥ-k to the reference wave function is replaced by the Fourier transformed multiplication of e-iλĤ, as a function of the integration parameters and the iteration power k for various systems, including H2, LiH, BeH2 and the notorious H4 molecule at square geometry. We further consider the possibility of employing the Gaussian-quadrature rule as an alternate integration method and compared it to the results employing trapezoidal integration. The Q-Inv algorithm is compared to the inverse iteration method using the Ĥ-1 inverse (I-Iter) and the exact inverse by lower-upper decomposition. Energy values are evaluated as the expectation values of the Hamiltonian. Results suggest that the Q-Inv method provides lower energy results than the I-Iter method up to a certain k, after which the energy increases due to errors in the numerical integration that are dependent on the integration interval. A combined Gaussian-quadrature and trapezoidal integration method proved to be more effective at reaching convergence while decreasing the number of operations. For systems like H4, in which the Q-Inv cannot reach the expected error threshold, we propose a combination of Q-Inv and I-Iter methods to further decrease the error with k at lower computational cost. Finally, we summarize the recommended procedure when treating unknown systems.

The Localized Active Space Method with Unitary Selective Coupled Cluster

We introduce a hybrid quantum-classical algorithm, the localized active space unitary selective coupled cluster singles and doubles (LAS-USCCSD) method. Derived from the localized active space unitary coupled cluster (LAS-UCCSD) method, LAS-USCCSD first performs a classical LASSCF calculation, then selectively identifies the most important parameters (cluster amplitudes used to build the multireference UCC ansatz) for restoring interfragment interaction energy using this reduced set of parameters with the variational quantum eigensolver method. We benchmark LAS-USCCSD against LAS-UCCSD by calculating the total energies of (H2)2, (H2)4, and trans-butadiene, and the magnetic coupling constant for a bimetallic compound [Cr2(OH)3(NH3)6]3+. For these systems, we find that LAS-USCCSD reduces the number of required parameters and thus the circuit depth by at least 1 order of magnitude, an aspect which is important for the practical implementation of multireference hybrid quantum-classical algorithms like LAS-UCCSD on near-term quantum computers.

Training Machine-Learned Density Functionals on Band Gaps

The systematic underestimation of band gaps is one of the most fundamental challenges in semilocal density functional theory (DFT). In addition to hindering the application of DFT to predicting electronic properties, the band gap problem is intimately related to self-interaction and delocalization errors, which make the study of charge transfer mechanisms with DFT difficult. To expand the range of available tools for addressing the band gap problem, we design an approach for machine learning density functionals based on Gaussian processes to explicitly fit single-particle energy levels. We also introduce nonlocal features of the density matrix that are expressive enough to fit these single-particle levels. Combining these developments, we train a machine-learned functional for the exact exchange energy that predicts molecular energy gaps and reaction energies of a wide range of molecules in excellent agreement with reference hybrid DFT calculations. In addition, while being trained solely on molecular data, our model predicts reasonable formation energies of polarons in solids, showcasing its transferability and robustness. We discuss how this approach can be generalized to full exchange-correlation functionals, thus paving the way to the design of state-of-the-art functionals for the prediction of electronic properties of molecules and materials.

Multi-GPU RI-HF Energies and Analytic Gradients─Toward High-Throughput Ab Initio Molecular Dynamics

This article presents an optimized algorithm and implementation for calculating resolution-of-the-identity Hartree-Fock (RI-HF) energies and analytic gradients using multiple graphics processing units (GPUs). The algorithm is especially designed for high throughput ab initio molecular dynamics simulations of small and medium size molecules (10-100 atoms). Key innovations of this work include the exploitation of multi-GPU parallelism and a workload balancing scheme that efficiently distributes computational tasks among GPUs. Our implementation also employs techniques for symmetry utilization, integral screening, and leveraging sparsity to optimize memory usage. Computational results show that the implementation achieves significant performance improvements, including over 3 × speedups in single GPU AIMD throughput compared to previous GPU-accelerated RI-HF and traditional HF methods. Furthermore, utilizing multiple GPUs can provide superlinear speedup when the additional aggregate GPU memory allows for the storage of decompressed three-center integrals.

Force-Free Identification of Minimum-Energy Pathways and Transition States for Stochastic Electronic Structure Theories

The accurate mapping of potential energy surfaces (PESs) is crucial to our understanding of the numerous physical and chemical processes mediated by atomic rearrangements, such as conformational changes and chemical reactions, and the thermodynamic and kinetic feasibility of these processes. Stochastic electronic structure theories, e.g., Quantum Monte Carlo (QMC) methods, enable highly accurate total energy calculations that in principle can be used to construct the PES. However, their stochastic nature poses a challenge to the computation and use of forces and Hessians, which are typically required in algorithms for minimum-energy pathway (MEP) and transition state (TS) identification, such as the nudged elastic band (NEB) algorithm and its climbing image formulation. Here, we present strategies that utilize the surrogate Hessian line-search method, previously developed for QMC structural optimization, to efficiently identify MEP and TS structures without requiring force calculations at the level of the stochastic electronic structure theory. By modifying the surrogate Hessian algorithm to operate in path-orthogonal subspaces and at saddle points, we show that it is possible to identify MEPs and TSs by using a force-free QMC approach. We demonstrate these strategies via two examples, the inversion of the ammonia (NH3) molecule and the nucleophilic substitution (SN2) reaction F- + CH3F → FCH3 + F-. We validate our results using Density Functional Theory (DFT)- and Coupled Cluster (CCSD, CCSD(T))-based NEB calculations. We then introduce a hybrid DFT-QMC approach to compute thermodynamic and kinetic quantities, free energy differences, rate constants, and equilibrium constants that incorporates stochastically optimized structures and their energies, and show that this scheme improves upon DFT accuracy. Our methods generalize straightforwardly to other systems and other high-accuracy theories that similarly face challenges computing energy gradients, paving the way for highly accurate PES mapping, transition state determination, and thermodynamic and kinetic calculations at significantly reduced computational expense.

Discriminating Clockwise and Counterclockwise Photoisomerization Paths in Achiral Photoswitches by Excited-State Electronic Circular Dichroism

Despite the numerous investigations of photoisomerization reactions from both the computational and experimental points of view, even in complex environments, to date there is no direct demonstration of the direction of rotation of the retinal chromophore, initiating the vision process in several organisms, occurring upon light irradiation. In the literature, many proposals have been formulated to shed light on the details of this process, most of which are extracted from semiclassical simulations. Although high hopes are held in the development of time-resolved X-ray spectroscopy, I argue in this work that simpler but less known techniques can be used to unravel the details of this fascinating photochemical process. In fact, chiroptical spectroscopy would unambiguously prove the direction of the rotatory motion of the chromophore during the photoisomerization process by probing excited state chirality, a piece of information that, so far, has been exclusively extracted from atomistic simulations. I demonstrate this statement by computing the expected chiroptical response along photoisomerization pathways for several models of the retinal chromophores that are found in nature bound to rhodopsins, including nuclear ensemble spectra from semiclassical dynamics simulations, that can be compared with time-resolved experiments.

Range-separated density functional theory using multiresolution analysis and quantum computing

Quantum computers are expected to outperform classical computers for specific problems in quantum chemistry. Such calculations remain expensive, but costs can be lowered through the partition of the molecular system. In the present study, partition was achieved with range-separated density functional theory (RS-DFT). The use of RS-DFT reduces both the basis set size and the active space size dependence of the ground state energy in comparison with the use of wave function theory (WFT) alone. The utilization of pair natural orbitals (PNOs) in place of canonical molecular orbitals (MOs) results in more compact qubit Hamiltonians. To test this strategy, a basis-set independent framework, known as multiresolution analysis (MRA), was employed to generate PNOs. Tests were conducted with the variational quantum eigensolver for a number of molecules. The results show that the proposed approach reduces the number of qubits needed to reach a target energy accuracy.

An Improved Penalty-Based Excited-State Variational Monte Carlo Approach with Deep-Learning Ansatzes

We introduce several improvements to the penalty-based variational quantum Monte Carlo (VMC) algorithm for computing electronic excited states of Entwistle et al. [Nat. Commun. 14, 274 (2023)] and demonstrate that the accuracy of the updated method is competitive with other available excited-state VMC approaches. A theoretical comparison of the computational aspects of these algorithms is presented, where several benefits of the penalty-based method are identified. Our main contributions include an automatic mechanism for tuning the scale of the penalty terms, an updated form of the overlap penalty with proven convergence properties, and a new term that penalizes the spin of the wave function, enabling the selective computation of states with a given spin. With these improvements, along with the use of the latest self-attention-based ansatz, the penalty-based method achieves a mean absolute error below 1 kcal/mol for the vertical excitation energies of a set of 26 atoms and molecules, without relying on variance matching schemes. Considering excited states along the dissociation of the carbon dimer, the accuracy of the penalty-based method is on par with that of natural-excited-state (NES) VMC, while also providing results for additional sections of the potential energy surface, which were inaccessible with the NES method. Additionally, the accuracy of the penalty-based method is improved for a conical intersection of ethylene, with the predicted angle of the intersection agreeing well with both NES-VMC and multireference configuration interaction.

Generalized Spin in the Variance-Based Wave Function Optimization Method within the Doubly Occupied Configuration Interaction Framework

In this work, we implement a generalized spin formulation of the doubly occupied configuration interaction methodology using the energy variance of the N-electron Hamiltonian. We perform the optimization of the N-electron wave functions and calculate their corresponding energies, using a unified variational treatment for ground and excited states based on the energy variance, which allows us to describe the entire energy spectra on an equal footing. We analyze the effects produced by the breakdown of the Ŝ2 and Ŝz symmetries in the spectra of model hydrogenic clusters in terms of energies and spin-related quantities, arising from the restricted, unrestricted, and generalized spin methods. The results are compared with other related methods as well as full configuration interaction.

Exchange-correlation potential built on the derivative discontinuity of electron density

Electronic structures are fully determined by the exchange-correlation (XC) potential. In this work, we develop a new method to construct reliable XC potentials by properly mixing the exact exchange and the local density approximation potentials in real space. The spatially dependent mixing parameter is derived based on the derivative discontinuity of electron density and is first-principle. We derived the equations for solving the mixing parameter and proposed an approximation to simplify these equations. Based on this approximation, this new method gives reasonable predictions for the ionization energies, fundamental gaps, and singlet-triplet energy differences for various molecular systems. The impact of the approximation on the constructed XC potentials is examined, and it is found that the quality of the XC potentials can be further improved by removing the approximation. This work demonstrates that the derivative discontinuity of electron density is a promising constraint for constructing high-quality XC potentials.

Ab initio extended Hubbard model of short polyenes for efficient quantum computing

We propose introducing an extended Hubbard Hamiltonian derived via the ab initio downfolding method, which was originally formulated for periodic materials, toward efficient quantum computing of molecular electronic structure calculations. By utilizing this method, the first-principles Hamiltonian of chemical systems can be coarse-grained by eliminating the electronic degrees of freedom in higher energy space and reducing the number of terms of electron repulsion integral from O(N4) to O(N2). Our approach is validated numerically on the vertical excitation energies and excitation characters of ethylene, butadiene, and hexatriene. The dynamical electron correlation is incorporated within the framework of the constrained random phase approximation in advance of quantum computations, and the constructed models capture the trend of experimental and high-level quantum chemical calculation results. As expected, the L1-norm of the fermion-to-qubit mapped model Hamiltonians is significantly lower than that of conventional ab initio Hamiltonians, suggesting improved scalability of quantum computing. Those numerical outcomes and the results of the simulation of excited-state sampling demonstrate that the ab initio extended Hubbard Hamiltonian may hold significant potential for quantum chemical calculations using quantum computers.

Ensemble Generalization of the Perdew-Zunger Self-Interaction Correction: A Way Out of Multiple Minima and Symmetry Breaking

The Perdew-Zunger (PZ) self-interaction correction (SIC) is an established tool to correct unphysical behavior in density functional approximations. Yet, the PZ-SIC is well-known to sometimes break molecular symmetries. An example of this is the benzene molecule, for which the PZ-SIC predicts a symmetry-broken electron density and molecular geometry, since the method does not describe the two possible Kekulé structures on an even footing, leading to local minima [Lehtola et al. J. Chem. Theory Comput. 2016, 12, 3195]. The PZ-SIC is often implemented with Fermi-Löwdin orbitals (FLOs), yielding the FLO-SIC method, which likewise has issues with symmetry breaking and local minima [Trepte et al. J. Chem. Phys. 2021, 155, 224109]. In this work, we propose a generalization of the PZ-SIC─the ensemble PZ-SIC (E-PZ-SIC) method─which shares the asymptotic computational scaling of the PZ-SIC (albeit with an additional prefactor). The E-PZ-SIC is straightforwardly applicable to various molecules, merely requiring one to average the self-interaction correction over all possible Kekulé structures, in line with chemical intuition. We showcase the implementation of the E-PZ-SIC with FLOs, as the resulting E-FLO-SIC method is easy to realize on top of an existing implementation of the FLO-SIC. We show that the E-FLO-SIC indeed eliminates symmetry breaking, reproducing a symmetric electron density and molecular geometry for benzene. The ensemble approach suggested herein could also be employed within approximate or locally scaled variants of the PZ-SIC and its FLO-SIC versions.

Estimating Molecular Thermal Averages with the Quantum Equation of Motion and Informationally Complete Measurements

By leveraging the Variational Quantum Eigensolver (VQE), the "quantum equation of motion" (qEOM) method established itself as a promising tool for quantum chemistry on near-term quantum computers and has been used extensively to estimate molecular excited states. Here, we explore a novel application of this method, employing it to compute thermal averages of quantum systems, specifically molecules like ethylene and butadiene. A drawback of qEOM is that it requires measuring the expectation values of a large number of observables on the ground state of the system, and the number of necessary measurements can become a bottleneck of the method. In this work, we focus on measurements through informationally complete positive operator-valued measures (IC-POVMs) to achieve a reduction in the measurement overheads by estimating different observables of interest through the measurement of a single set of POVMs. We show with numerical simulations that the qEOM combined with IC-POVM measurements ensures satisfactory accuracy in the reconstruction of the thermal state with a reasonable number of shots.

MP2-Based Composite Extrapolation Schemes Can Predict Core-Ionization Energies for First-Row Elements with Coupled-Cluster Level Accuracy

X-ray photoelectron spectroscopy (XPS) measures core-electron binding energies (CEBEs) to reveal element-specific insights into the chemical environment and bonding. Accurate theoretical CEBE prediction aids XPS interpretation but requires proper modeling of orbital relaxation and electron correlation upon core-ionization. This work systematically investigates basis set selection for extrapolation to the complete basis set limit of CEBEs from ΔMP2 and ΔCC energies across 94 K-edges in diverse organic molecules. We demonstrate that an alternative composite scheme using ΔMP2 in a large basis corrected by ΔCC-ΔMP2 difference in a small basis can quantitatively recover optimally extrapolated ΔCC CEBEs within 0.02 eV. Unlike ΔCC, MP2 calculations do not suffer from convergence issues and are computationally cheaper, and thus, the composite ΔMP2/ΔCC scheme balances accuracy and cost, overcoming limitations of solely using either method. We conclude by providing a comprehensive analysis of the choice of small and large basis sets for the composite schemes and provide practical recommendations for highly accurate (within 0.10-0.15 eV MAE) ab initio prediction of XPS data.

Renormalized Internally Contracted Multireference Coupled Cluster with Perturbative Triples

In this work, we combine the many-body formulation of the internally contracted multireference coupled cluster (ic-MRCC) method with Evangelista's multireference formulation of the driven similarity renormalization group (DSRG). The DSRG method can be viewed as a unitary multireference coupled cluster theory, which renormalizes the amplitudes based on a flow equation approach to eliminate numerical instabilities. We extend this approach by demonstrating that the unitary flow equation approach can be adapted for nonunitary transformations, rationalizing the renormalization of ic-MRCC amplitudes. We denote the new approach, the renormalized ic-MRCC (ric-MRCC) method. To achieve high accuracy with a reasonable computational cost, we introduce a new approximation to the Baker-Campbell-Hausdorff expansion. We fully consider the linear commutator while approximating the quadratic commutator, for which we neglect specific contractions involving amplitudes with active indices. Moreover, we introduce approximate perturbative triples to obtain the ric-MRCCSD[T] method. We demonstrate the accuracy of our approaches in comparison to advanced multireference methods for the potential energy curves of H8, F2, H2O, N2, and Cr2. Additionally, we show that ric-MRCCSD and ric-MRCSSD[T] match the accuracy of CCSD(T) for evaluating spectroscopic constants and of full configuration interaction energies for a set of small molecules.

Error Control and Automatic Detection of Reference Active Spaces in Many-Body Expanded Full Configuration Interaction

We present a wide-reaching revamp of the generalized many-body expanded full configuration interaction (MBE-FCI) method. First, we outline how to automatize the selection of reference active spaces, whereby the inherent bias introduced through a manual identification is reduced, also within the context of traditional complete active space methods. Second, we allow for the use of compact orbital clusters as expansion objects, which works to circumvent the unfavorable scaling with the number of orbitals included in the space complementary to the reference orbitals. Finally, we present a new algorithm for ensuring that many-body expansions can be efficiently terminated while conservatively accounting for resulting errors. These developments are all tested on a variety of molecular systems and different orbital representations to illustrate the abilities of our algorithm to produce correlation energies within predetermined error bounds, significantly broadening the overall applicability of the MBE-FCI method.

Forte: A suite of advanced multireference quantum chemistry methods

Forte is an open-source library specialized in multireference electronic structure theories for molecular systems and the rapid prototyping of new methods. This paper gives an overview of the capabilities of Forte, its software architecture, and examples of applications enabled by the methods it implements.