Toward Real Chemical Accuracy on Current Quantum Hardware Through the Transcorrelated Method.

Quantum computing is emerging as a new computational paradigm with the potential to transform several research fields including quantum chemistry. However, current hardware limitations (including limited coherence times, gate infidelities, and connectivity) hamper the implementation of most quantum algorithms and call for more noise-resilient solutions. We propose an explicitly correlated Ansatz based on the transcorrelated (TC) approach to target these major roadblocks directly. This method transfers, without any approximation, correlations from the wave function directly into the Hamiltonian, thus reducing the resources needed to achieve accurate results with noisy quantum devices. We show that the TC approach allows for shallower circuits and improves the convergence toward the complete basis set limit, providing energies within chemical accuracy to experiment with smaller basis sets and, thus, fewer qubits. We demonstrate our method by computing bond lengths, dissociation energies, and vibrational frequencies close to experimental results for the hydrogen dimer and lithium hydride using two and four qubits, respectively. To demonstrate our approach's current and near-term potential, we perform hardware experiments, where our results confirm that the TC method paves the way toward accurate quantum chemistry calculations already on today's quantum hardware.

The electron density: a fidelity witness for quantum computation.

There is currently no combination of quantum hardware and algorithms that can provide an advantage over conventional calculations of molecules or materials. However, if or when such a point is reached, new strategies will be needed to verify predictions made using quantum devices. We propose that the electron density, obtained through experimental or computational means, can serve as a robust benchmark for validating the accuracy of quantum computation of chemistry. An initial exploration into topological features of electron densities, facilitated by quantum computation, is presented here as a proof of concept. Additionally, we examine the effects of constraining and symmetrizing measured one-particle reduced density matrices on noise-driven errors in the electron density distribution. We emphasize the potential benefits and future need for high-quality electron densities derived from diffraction experiments for validating classically intractable quantum computations of materials.

Optimizing Jastrow factors for the transcorrelated method.

We investigate the optimization of flexible tailored real-space Jastrow factors for use in the transcorrelated (TC) method in combination with highly accurate quantum chemistry methods, such as initiator full configuration interaction quantum Monte Carlo (FCIQMC). Jastrow factors obtained by minimizing the variance of the TC reference energy are found to yield better, more consistent results than those obtained by minimizing the variational energy. We compute all-electron atomization energies for the challenging first-row molecules C2, CN, N2, and O2 and find that the TC method yields chemically accurate results using only the cc-pVTZ basis set, roughly matching the accuracy of non-TC calculations with the much larger cc-pV5Z basis set. We also investigate an approximation in which pure three-body excitations are neglected from the TC-FCIQMC dynamics, saving storage and computational costs, and show that it affects relative energies negligibly. Our results demonstrate that the combination of tailored real-space Jastrow factors with the multi-configurational TC-FCIQMC method provides a route to obtaining chemical accuracy using modest basis sets, obviating the need for basis-set extrapolation and composite techniques.

The OpenMolcas Web: A Community-Driven Approach to Advancing Computational Chemistry.

The developments of the open-source OpenMolcas chemistry software environment since spring 2020 are described, with a focus on novel functionalities accessible in the stable branch of the package or via interfaces with other packages. These developments span a wide range of topics in computational chemistry and are presented in thematic sections: electronic structure theory, electronic spectroscopy simulations, analytic gradients and molecular structure optimizations, ab initio molecular dynamics, and other new features. This report offers an overview of the chemical phenomena and processes OpenMolcas can address, while showing that OpenMolcas is an attractive platform for state-of-the-art atomistic computer simulations.

Reference-State Error Mitigation: A Strategy for High Accuracy Quantum Computation of Chemistry.

Decoherence and gate errors severely limit the capabilities of state-of-the-art quantum computers. This work introduces a strategy for reference-state error mitigation (REM) of quantum chemistry that can be straightforwardly implemented on current and near-term devices. REM can be applied alongside existing mitigation procedures, while requiring minimal postprocessing and only one or no additional measurements. The approach is agnostic to the underlying quantum mechanical ansatz and is designed for the variational quantum eigensolver. Up to two orders-of-magnitude improvement in the computational accuracy of ground state energies of small molecules (H2, HeH+, and LiH) is demonstrated on superconducting quantum hardware. Simulations of noisy circuits with a depth exceeding 1000 two-qubit gates are used to demonstrate the scalability of the method.

Performance of a one-parameter correlation factor for transcorrelation: Study on a series of second row atomic and molecular systems.

In this work, we investigate the performance of a recently proposed transcorrelated (TC) approach based on a single-parameter correlation factor [E. Giner, J. Chem. Phys. 154, 084119 (2021)] for systems involving more than two electrons. The benefit of such an approach relies on its simplicity as efficient numerical-analytical schemes can be set up to compute the two- and three-body integrals occurring in the effective TC Hamiltonian. To obtain accurate ground state energies within a given basis set, the present TC scheme is coupled to the recently proposed TC-full configuration interaction quantum Monte Carlo method [Cohen et al., J. Chem. Phys. 151, 061101 (2019)]. We report ground state total energies on the Li-Ne series, together with their first cations, computed with increasingly large basis sets and compare to more elaborate correlation factors involving electron-electron-nucleus coordinates. Numerical results on the Li-Ne ionization potentials show that the use of the single-parameter correlation factor brings on average only a slightly lower accuracy (1.2 mH) in a triple-zeta quality basis set with respect to a more sophisticated correlation factor. However, already using a quadruple-zeta quality basis set yields results within chemical accuracy to complete basis set limit results when using this novel single-parameter correlation factor. Calculations on the H2O, CH2, and FH molecules show that a similar precision can be obtained within a triple-zeta quality basis set for the atomization energies of molecular systems.

Spin-Pure Stochastic-CASSCF via GUGA-FCIQMC Applied to Iron-Sulfur Clusters.

In this work, we demonstrate how to efficiently compute the one- and two-body reduced density matrices within the spin-adapted full configuration interaction quantum Monte Carlo (FCIQMC) method, which is based on the graphical unitary group approach (GUGA). This allows us to use GUGA-FCIQMC as a spin-pure configuration interaction (CI) eigensolver within the complete active space self-consistent field (CASSCF) procedure and hence to stochastically treat active spaces far larger than conventional CI solvers while variationally relaxing orbitals for specific spin-pure states. We apply the method to investigate the spin ladder in iron-sulfur dimer and tetramer model systems. We demonstrate the importance of the orbital relaxation by comparing the Heisenberg model magnetic coupling parameters from the CASSCF procedure to those from a CI-only (CASCI) procedure based on restricted open-shell Hartree-Fock orbitals. We show that the orbital relaxation differentially stabilizes the lower-spin states, thus enlarging the coupling parameters with respect to the values predicted by ignoring orbital relaxation effects. Moreover, we find that, while CASCI results are well fit by a simple bilinear Heisenberg Hamiltonian, the CASSCF eigenvalues exhibit deviations that necessitate the inclusion of biquadratic terms in the model Hamiltonian.

Resolution of Low-Energy States in Spin-Exchange Transition-Metal Clusters: Case Study of Singlet States in [Fe(III)4S4] Cubanes.

Polynuclear transition-metal (PNTM) clusters owe their catalytic activity to numerous energetically low-lying spin states and stable oxidation states. The characterization of their electronic structure represents one of the greatest challenges of modern chemistry. We propose a theoretical framework that enables the resolution of targeted electronic states with ease and apply it to two [Fe(III)4S4] cubanes. Through direct access to their many-body wave functions, we identify important correlation mechanisms and their interplay with the geometrical distortions observed in these clusters, which are core properties in understanding their catalytic activity. The simulated magnetic coupling constants predicted by our strategy allow us to make qualitative connections between spin interactions and geometrical distortions, demonstrating its predictive power. Moreover, despite its simplicity, the strategy provides magnetic coupling constants in good agreement with the available experimental ones. The complexes are intrinsically frustrated anti-ferromagnets, and the obtained spin structures together with the geometrical distortions represent two possible ways to release spin frustration (spin-driven Jahn-Teller distortion). Our paradigm provides a simple, yet rigorous, route to uncover the electronic structure of PNTM clusters and may be applied to a wide variety of such clusters.

NECI: N-Electron Configuration Interaction with an emphasis on state-of-the-art stochastic methods.

We present NECI, a state-of-the-art implementation of the Full Configuration Interaction Quantum Monte Carlo (FCIQMC) algorithm, a method based on a stochastic application of the Hamiltonian matrix on a sparse sampling of the wave function. The program utilizes a very powerful parallelization and scales efficiently to more than 24 000 central processing unit cores. In this paper, we describe the core functionalities of NECI and its recent developments. This includes the capabilities to calculate ground and excited state energies, properties via the one- and two-body reduced density matrices, as well as spectral and Green's functions for ab initio and model systems. A number of enhancements of the bare FCIQMC algorithm are available within NECI, allowing us to use a partially deterministic formulation of the algorithm, working in a spin-adapted basis or supporting transcorrelated Hamiltonians. NECI supports the FCIDUMP file format for integrals, supplying a convenient interface to numerous quantum chemistry programs, and it is licensed under GPL-3.0.

Compression of Spin-Adapted Multiconfigurational Wave Functions in Exchange-Coupled Polynuclear Spin Systems.

We present a protocol based on unitary transformations of molecular orbitals to reduce the number of nonvanishing coefficients of spin-adapted configuration interaction expansions. Methods that exploit the sparsity of the Hamiltonian matrix and compactness of its eigensolutions, such as the full configuration interaction quantum Monte Carlo (FCIQMC) algorithm in its spin-adapted implementation, are well suited to this protocol. The wave function compression resulting from this approach is particularly attractive for antiferromagnetically coupled polynuclear spin systems, such as transition-metal cubanes in biocatalysis, and Mott and charge-transfer insulators in solid-state physics. Active space configuration interaction calculations on N2 and CN- at various bond lengths, the stretched square N4 compounds, the chromium dimer, and a [Fe2S2]2- model system are presented as a proof-of-concept. For the Cr2 case, large and intermediate bond distances are discussed, showing that the approach is effective in cases where static and dynamic correlations are equally important. The [Fe2S2]2- case shows the general applicability of the method.

Efficient formulation of full configuration interaction quantum Monte Carlo in a spin eigenbasis via the graphical unitary group approach.

We provide a spin-adapted formulation of the Full Configuration Interaction Quantum Monte Carlo (FCIQMC) algorithm, based on the Graphical Unitary Group Approach (GUGA), which enables the exploitation of SU(2) symmetry within this stochastic framework. Random excitation generation and matrix element calculation on the Shavitt graph of GUGA can be efficiently implemented via a biasing procedure on the branching diagram. The use of a spin-pure basis explicitly resolves the different spin-sectors and ensures that the stochastically sampled wavefunction is an eigenfunction of the total spin operator S^2. The method allows for the calculation of states with low or intermediate spin in systems dominated by Hund's first rule, which are otherwise generally inaccessible. Furthermore, in systems with small spin gaps, the new methodology enables much more rapid convergence with respect to walker number and simulation time. Some illustrative applications of the GUGA-FCIQMC method are provided: computation of the 2F - 4F spin gap of the cobalt atom in large basis sets, achieving chemical accuracy to experiment, and the Σg+1, Σg+3, Σg+5, and Σg+7 spin-gaps of the stretched N2 molecule, an archetypal strongly correlated system.

Time Propagation and Spectroscopy of Fermionic Systems Using a Stochastic Technique.

We present a stochastic method for solving the time-dependent Schrödinger equation, generalizing a ground state full configuration interaction quantum Monte Carlo method. By performing the time integration in the complex plane close to the real-time axis, the numerical effort is kept manageable and the analytic continuation to real frequencies is efficient. This allows us to perform ab initio calculation of electron spectra for strongly correlated systems. The method can be used as a cluster solver for embedding schemes.