Benchmark Study of Core-Ionization Energies with the Generalized Active Space-Driven Similarity Renormalization Group.

X-ray photoelectron spectroscopy (XPS) is a powerful experimental technique for probing the electronic structure of molecules and materials; however, interpreting XPS data requires accurate computational methods to model core-ionized states. This work proposes and benchmarks a new approach based on the generalized active space-driven similarity renormalization group (GAS-DSRG) for calculating core-ionization energies and treating correlation effects at the perturbative and nonperturbative levels. We tested the GAS-DSRG across three data sets. First, the vertical core-ionization energies of small molecules containing first-row elements are evaluated. GAS-DSRG achieves mean absolute errors below 0.3 eV, which is comparable to high-level coupled cluster methods. Next, the accuracy of GAS-DSRG is evaluated for larger organic molecules using the CORE65 data set, with the DSRG-MRPT3 level yielding a mean absolute error of only 0.34 eV for 65 core-ionization transitions. Insights are provided into the treatment of static and dynamic correlation, the importance of high-order perturbation theory, and notable differences from density functional theory in the predicted energy ordering of core-ionized states for specific molecules. Finally, vibrationally resolved XPS spectra of diatomic molecules (CO, N2, and O2) are simulated, showing excellent agreement with experimental data.

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.

Toward Accurate Spin-Orbit Splittings from Relativistic Multireference Electronic Structure Theory.

Most nonrelativistic electron correlation methods can be adapted to account for relativistic effects, as long as the relativistic molecular spinor integrals are available, from either a four-, two-, or one-component mean-field calculation. However, relativistic multireference correlation methods remain a relatively unexplored area, with mixed evidence regarding the improvements brought by perturbative treatments. We report, for the first time, the implementation of state-averaged four-component relativistic multireference perturbation theories to second and third order based on the driven similarity renormalization group (DSRG). With our methods, named 4c-SA-DSRG-MRPT2 and 3, we find that the dynamical correlation included on top of 4c-CASSCF references can significantly improve the spin-orbit splittings in p-block elements and potential energy surfaces when compared to 4c-CASSCF and 4c-CASPT2 results. We further show that 4c-DSRG-MRPT2 and 3 are applicable to these systems over a wide range of the flow parameter, with systematic improvement from second to third order in terms of both improved error statistics and reduced sensitivity with respect to the flow parameter.

Frozen Natural Orbitals for the State-Averaged Driven Similarity Renormalization Group.

We present a reduced-cost implementation of the state-averaged driven similarity renormalization group (SA-DSRG) based on the frozen natural orbital (FNO) approach. The natural orbitals (NOs) are obtained by diagonalizing the one-body reduced density matrix from SA-DSRG second-order perturbation theory (SA-DSRG-PT2). We consider three criteria to truncate the virtual NOs for the subsequent electron correlation treatment beyond SA-DSRG-PT2. An additive second-order correction is applied to the SA-DSRG Hamiltonian to reintroduce correlation effects from the discarded orbitals. The FNO SA-DSRG method is benchmarked on 35 small organic molecules in the QUEST database. When keeping 98-99% of the cumulative occupation numbers, the mean absolute error in the vertical transition energies due to FNO is less than 0.01 eV. Using the same FNO threshold, we observe a speedup of 9 times compared to the conventional SA-DSRG implementation for nickel carbonyl with a quadruple-ζ basis set. The FNO approach enables nonperturbative SA-DSRG computations on chloroiron corrole [FeCl(C19H11N4)] with more than 1000 basis functions, surpassing the current limit of a conventional implementation.

Implementation of the Projective Quantum Eigensolver on a Quantum Computer.

We study the performance of our previously proposed projective quantum eigensolver (PQE) on IBM's quantum hardware in conjunction with error mitigation techniques. For a single qubit model of H2, we find that we are able to obtain energies within 4 millihartree (2.5 kcal/mol) of the exact energy along the entire potential energy curve, with the accuracy limited by both the stochastic error and the inconsistent performance of the IBM devices. We find that an optimization algorithm using direct inversion of the iterative subspace can converge swiftly, even to excited states, but stochastic noise can prompt large parameter updates. For the 4-site transverse-field Ising model at its critical point, PQE with an appropriate application of qubit tapering can recover 99% of the correlation energy, even after discarding several parameters. The large number of CNOT gates needed for the additional parameters introduces a concomitant error that, on the IBM devices, results in a loss of accuracy despite the increased expressivity of the trial state. Error extrapolation techniques and tapering or postselection are recommended to mitigate errors in PQE hardware experiments.

Intruder-free cumulant-truncated driven similarity renormalization group second-order multireference perturbation theory.

Accurate multireference electronic structure calculations are important for constructing potential energy surfaces. Still, even in the case of low-scaling methods, their routine use is limited by the steep growth of the computational and storage costs as the active space grows. This is primarily due to the occurrence of three- and higher-body density matrices or, equivalently, their cumulants. This work examines the effect of various cumulant truncation schemes on the accuracy of the driven similarity renormalization group second-order multireference perturbation theory. We test four different levels of three-body reduced density cumulant truncations that set different classes of cumulant elements to zero. Our test cases include the singlet-triplet gap of CH2, the potential energy curves of the XΣg+1 and AΣu+3 states of N2, and the singlet-triplet splittings of oligoacenes. Our results show that both relative and absolute errors introduced by these cumulant truncations can be as small as 0.5 kcal mol-1 or less. At the same time, the amount of memory required is reduced from O(NA6) to O(NA5), where NA is the number of active orbitals. No additional regularization is needed to prevent the intruder state problem in the cumulant-truncated second-order driven similarity renormalization group multireference perturbation theory methods.

Unitary Coupled Cluster: Seizing the Quantum Moment.

Shallow, CNOT-efficient quantum circuits are crucial for performing accurate computational chemistry simulations on current noisy quantum hardware. Here, we explore the usefulness of noniterative energy corrections, based on the method of moments of coupled-cluster theory, for accelerating convergence toward full configuration interaction. Our preliminary numerical results relying on iteratively constructed ansätze suggest that chemically accurate energies can be obtained with substantially more compact circuits, implying enhanced resilience to gate and decoherence noise.

Linear-Scaling Quantum Circuits for Computational Chemistry.

We have recently constructed compact, CNOT-efficient, quantum circuits for Fermionic and qubit excitations of arbitrary many-body rank [Magoulas, I.; Evangelista, F. A. J. Chem. Theory Comput. 2023, 19, 822]. Here, we present approximations of these circuits that substantially reduce the CNOT counts even further. Our preliminary numerical data, using the selected projective quantum eigensolver approach, show up to a 4-fold reduction in CNOTs. At the same time, there is practically no loss of accuracy in the energies compared to the parent implementation, while the ensuing symmetry breaking is essentially negligible.

A study of core-excited states of organic molecules computed with the generalized active space driven similarity renormalization group.

This work examines the accuracy and precision of x-ray absorption spectra computed with a multireference approach that combines generalized active space (GAS) references with the driven similarity renormalization group (DSRG). We employ the x-ray absorption benchmark of organic molecule (XABOOM) set, consisting of 116 transitions from mostly organic molecules [Fransson et al., J. Chem. Theory Comput. 17, 1618 (2021)]. Several approximations to a full-valence active space are examined and benchmarked. Absolute excitation energies and intensities computed with the GAS-DSRG truncated to second-order in perturbation theory are found to systematically underestimate experimental and reference theoretical values. Third-order perturbative corrections significantly improve the accuracy of GAS-DSRG absolute excitation energies, bringing the mean absolute deviation from experimental values down to 0.32 eV. The ozone molecule and glyoxylic acid are particularly challenging for second-order perturbation theory and are examined in detail to assess the importance of active space truncation and intruder states.

Signatures of diradicals in x-ray absorption spectroscopy.

Theoretical simulations are critical to analyze and interpret the x-ray absorption spectrum of transient open-shell species. In this work, we propose a model of the many-body core-excited states of symmetric diradicals. We apply this model to analyze the carbon K-edge transitions of o-, m-, and p-benzyne, three organic diradicals with diverse and unusual electronic structures. The predictions of our model are compared with high-level multireference computations of the K-edge spectrum of the benzynes obtained with the driven similarity renormalization group truncated to third order. Our model shows the importance of a many-body treatment of the core-excited states of the benzynes and provides a theoretical framework to understand which properties of the ground state of these diradicals can be extracted from their x-ray spectrum.

Ring-Polymer Instanton Tunneling Splittings of Tropolone and Isotopomers using a Δ-Machine Learned CCSD(T) Potential: Theory and Experiment Shake Hands.

Tropolone, a 15-atom cyclic molecule, has received much interest both experimentally and theoretically due to its H-transfer tunneling dynamics. An accurate theoretical description is challenging owing to the need to develop a high-level potential energy surface (PES) and then to simulate quantum-mechanical tunneling on this PES in full dimensionality. Here, we tackle both aspects of this challenge and make detailed comparisons with experiments for numerous isotopomers. The PES, of near CCSD(T)-quality, is obtained using a Δ-machine learning approach starting from a pre-existing low-level DFT PES and corrected by a small number of approximate CCSD(T) energies obtained using the fragmentation-based molecular tailoring approach. The resulting PES is benchmarked against DF-FNO-CCSD(T) and CCSD(T)-F12 calculations. Ring-polymer instanton calculations of the splittings, obtained with the Δ-corrected PES are in good agreement with previously reported experiments and a significant improvement over those obtained using the low-level DFT PES. The instanton path includes heavy-atom tunneling effects and cuts the corner, thereby avoiding passing through the conventional saddle-point transition state. This is in contradistinction with typical approaches based on the minimum-energy reaction path. Finally, the subtle changes in the splittings for some of the heavy-atom isotopomers seen experimentally are reproduced and explained.

CO Inversion on a NaCl(100) Surface: A Multireference Quantum Embedding Study.

We develop a multireference quantum embedding model to investigate a recent experimental observation of the isomerization of vibrationally excited CO molecules on a NaCl(100) surface [Science 2020, 367, 175-178]. To explore this mechanism, we built a reduced potential energy surface of CO interacting with NaCl(100) using a second-order multireference perturbation theory, modeling the adsorbate-surface interaction with our previously developed active space embedding theory (ASET). We considered an isolated CO molecule on NaCl(100) and a high-coverage CO monolayer (1/1), and for both we generated potential energy surfaces parametrized by the CO stretching, adsorption, and inversion coordinates. These surfaces are used to determine stationary points and adsorption energies and to perform a vibrational analysis of the states relevant to the inversion mechanism. We found that for near-equilibrium bond lengths, CO adsorbed in the C-down configuration is lower in energy than in the O-down configuration. Stretching of the C-O bond reverses the energetic order of these configurations, supporting the accepted isomerization mechanism. The vibrational constants obtained from these potential energy surfaces show a small (< 10 cm-1) blue- and red-shift for the C-down and O-down configurations, respectively, in agreement with experimental assignments and previous theoretical studies. Our vibrational analysis of the monolayer case suggests that the O-down configuration is energetically more stable than the C-down one beyond the 16th vibrational excited state of CO, a value slightly smaller than the one from quasi-classical trajectory simulations (22nd) and consistent with the experiment. Our analysis suggests that CO-CO interactions in the monolayer play an important role in stabilizing highly vibrationally excited states in the O-down configuration and reducing the barrier between the C-down and O-down geometries, therefore playing a crucial role in the inversion mechanism.

CNOT-Efficient Circuits for Arbitrary Rank Many-Body Fermionic and Qubit Excitations.

Efficient quantum circuits are necessary for realizing quantum algorithms on noisy intermediate-scale quantum devices. Fermionic excitations entering unitary coupled-cluster (UCC) ansätze give rise to quantum circuits containing CNOT "staircases" whose number scales exponentially with the excitation rank. Recently, Yordanov et al. [Phys. Rev. A102, 062612 (2020); Commun. Phys.4, 228 (2021)] constructed CNOT-efficient quantum circuits for both Fermionic- (FEB) and qubit-excitation-based (QEB) singles and doubles and illustrated their usefulness in adaptive derivative-assembled pseudo-Trotterized variational quantum eigensolver (ADAPT-VQE) simulations. In this work, we extend these CNOT-efficient quantum circuits to arbitrary excitation ranks. To illustrate the benefits of these compact FEB and QEB quantum circuits, we perform numerical simulations using the recently developed selected projective quantum eigensolver (SPQE) approach, which relies on an adaptive UCC ansatz built from arbitrary-order particle-hole excitation operators. We show that both FEB- and QEB-SPQE decrease the number of CNOT gates compared to traditional SPQE by factors as large as 15. At the same time, QEB-SPQE requires, in general, more ansatz parameters than FEB-SPQE, in particular those corresponding to higher-than-double excitations, resulting in quantum circuits with larger CNOT counts. Although ADAPT-VQE generates quantum circuits with fewer CNOTs than SPQE, SPQE requires orders of magnitude less residual element evaluations than gradient element evaluations in ADAPT-VQE.

Automatic derivation of many-body theories based on general Fermi vacua.

This paper describes Wick&d, an implementation of the algebra of second-quantized operators normal ordered with respect to general correlated references and the corresponding Wick theorem [D. Mukherjee, Chem. Phys. Lett. 274, 561 (1997) and W. Kutzelnigg and D. Mukherjee, J. Chem. Phys. 107, 432 (1997)]. Wick&d employs a compact representation of operators and a backtracking algorithm to efficiently evaluate Wick contractions. Since Wick&d can handle both fully and partially contracted terms, it can be applied to both projective and Fock-space many-body formalisms. To demonstrate the usefulness of Wick&d, we use it to evaluate the single-reference coupled cluster equations up to octuple excitations and report an automated derivation and implementation of the second-order driven similarity renormalization group multireference perturbation theory.

Second-Order Active-Space Embedding Theory.

Quantum embedding schemes are a promising way to extend multireference computations to large molecules with strong correlation effects localized on a small number of atoms. This work introduces a second-order active-space embedding theory [ASET(2)] which improves upon mean-field frozen embedding by treating fragment-environment interactions via an approximate canonical transformation. The canonical transformation employed in ASET(2) is formulated using the driven similarity renormalization group. The ASET(2) scheme is benchmarked on the N═N bond dissociation in pentyldiazene, the S0 to S1 excitation in 1-octene, and the interaction energy of the O2-benzene complex. The ASET(2) explicit treatment of fragment-environment interactions beyond the mean-field level generally improves the accuracy of embedded computations, and it becomes necessary to achieve an accurate description of excitation energies of 1-octene and the singlet-triplet gap of the O2-benzene complex.

QForte: An Efficient State-Vector Emulator and Quantum Algorithms Library for Molecular Electronic Structure.

We introduce a novel open-source software package QForte, a comprehensive development tool for new quantum simulation algorithms. QForte incorporates functionality for handling molecular Hamiltonians, Fermionic encoding, ansatz construction, time evolution, and state-vector emulation, requiring only a classical electronic structure package as a dependency. QForte also contains black-box implementations of a wide variety of quantum algorithms, including variational and projective quantum eigensolvers, adaptive eigensolvers, quantum imaginary time evolution, and quantum Krylov methods. We highlight two features of QForte: (i) how the Python class structure of QForte enables the facile implementation of new algorithms, and (ii) how existing algorithms can be executed in just a few lines of code.

Theoretical Calculation of Core-Excited States along Dissociative Pathways beyond Second-Order Perturbation Theory.

We extend the multireference driven similarity renormalization (MR-DSRG) method to compute core-excited states by combining it with a GASSCF treatment of orbital relaxation and static electron correlation effects. We consider MR-DSRG treatments of dynamical correlation truncated at the level of perturbation theory (DSRG-MRPT2/3) and iterative linearized approximations with one- and two-body operators [MR-LDSRG(2)] in combination with a spin-free exact-two-component (X2C) one-electron treatment of scalar relativistic effects. This approach is calibrated and tested on a series of 16 core-excited states of five closed- and open-shell diatomic molecules containing first-row elements (C, N, and O). All GASSCF-MR-DSRG theories show excellent agreement with experimental adiabatic transitions energies, with mean absolute errors ranging between 0.17 and 0.35 eV, even for the challenging partially doubly excited states of the N2+ molecule. The vibrational structure of all these transitions, obtained from using a full potential energy scan, shows a mean absolute error as low as 25 meV for DSRG-MRPT2 and 12/13 meV for DSRG-MRPT3 and MR-LDSRG(2). We generally find that a treatment of dynamical correlation that goes beyond the second-order level in perturbation theory improves the accuracy of the potential energy surface, especially in the bond-dissociation region.

Analytic Energy Gradients for the Driven Similarity Renormalization Group Multireference Second-Order Perturbation Theory.

We derive analytic energy gradients of the driven similarity renormalization group (DSRG) multireference second-order perturbation theory (MRPT2) using the method of Lagrange multipliers. In the Lagrangian, we impose constraints for a complete-active-space self-consistent-field reference wave function and the semicanonical orthonormal molecular orbitals. Solving for the associated Lagrange multipliers is found to share the same asymptotic scaling of a single DSRG-MRPT2 energy computation. A pilot implementation of the DSRG-MRPT2 analytic gradients is used to optimize the geometry of the singlet and triplet states of p-benzyne. The equilibrium bond lengths and angles are similar to those computed via other MRPT2s and Mukherjee's multireference coupled cluster theory. An approximate DSRG-MRPT2 method that neglects the contributions of the three-body density cumulant is found to introduce negligible errors in the geometry of p-benzyne, lending itself to a promising low-cost approach for molecular geometry optimizations using large active spaces.

Spin-free formulation of the multireference driven similarity renormalization group: A benchmark study of first-row diatomic molecules and spin-crossover energetics.

We report a spin-free formulation of the multireference (MR) driven similarity renormalization group (DSRG) based on the ensemble normal ordering of Mukherjee and Kutzelnigg [J. Chem. Phys. 107, 432 (1997)]. This ensemble averages over all microstates of a given total spin quantum number, and therefore, it is invariant with respect to SU(2) transformations. As such, all equations may be reformulated in terms of spin-free quantities and they closely resemble those of spin-adapted closed-shell coupled cluster (CC) theory. The current implementation is used to assess the accuracy of various truncated MR-DSRG methods (perturbation theory up to third order and iterative methods with single and double excitations) in computing the constants of 33 first-row diatomic molecules. The accuracy trends for these first-row diatomics are consistent with our previous benchmark on a small subset of closed-shell diatomic molecules. We then present the first MR-DSRG application on transition-metal complexes by computing the spin splittings of the [Fe(H2O)6]2+ and [Fe(NH3)6]2+ molecules. A focal point analysis (FPA) shows that third-order perturbative corrections are essential to achieve reasonably converged energetics. The FPA based on the linearized MR-DSRG theory with one- and two-body operators and up to a quintuple-ζ basis set predicts the spin splittings of [Fe(H2O)6]2+ and [Fe(NH3)6]2+ to be -35.7 and -17.1 kcal mol-1, respectively, showing good agreement with the results of local CC theory with singles, doubles, and perturbative triples.