Spin-Flip Unitary Coupled Cluster Method: Toward Accurate Description of Strong Electron Correlation on Quantum Computers

Unitary Block-Correlated Coupled Cluster Ansatz Based on the Generalized Valence Bond Wave Function for Quantum Simulation

Strongly correlated (SC) systems present significant challenges for classical quantum chemistry methods. Quantum computing, particularly the variational quantum eigensolver (VQE), offers a promising framework to address these challenges by inherently supporting exponentially large configuration spaces. However, its application to SC systems remains limited due to the single-reference nature of the widely used ansatzes such as unitary coupled cluster (UCC). To address this challenge, we propose the generalized valence bond-based unitary block correlated coupled cluster (GVB-UBCCC) method. This novel ansatz incorporates the multiconfigurational nature of generalized valence bond (GVB) and the accuracy of block correlated coupled cluster (BCCC) methods, making it well-suited for SC systems. We have implemented the GVB-UBCCC method with up to two-block correlation (GVB-UBCCC2) and applied it to investigate ground-state energies for several SC systems, including H4, the water dimer, N2H2, and S6, at most described by 24 qubits. Our approach demonstrates that for these systems, GVB-UBCCC2 can achieve more accurate ground-state energies than UCCSD in most cases while requiring only O(N2) quantum gates and parameters, as opposed to the O(N4) scaling of UCCSD. The results highlight the effectiveness and potential advantages of GVB-UBCCC in SC systems.

Vibronic Structure of the UV/Visible Absorption Spectra of Phenol and Phenolate: A Hybrid Density Functional Theory─Doktorov's Quantum Algorithm Approach

The Doktorov's quantum algorithm has been enacted in combination with time-dependent density functional theory (TD-DFT) to simulate the vibronic structure of the UV/visible absorption spectra of the phenol and phenolate molecules. On the one hand, DFT and TD-DFT are employed with classical algorithms to calculate the ground and excited-state electronic structures as well as their vibrational frequencies and normal modes, whereas, on the other hand, quantum algorithms are employed for evaluating the vibrational transition intensities. In comparison to a previous study, J. Phys. Chem. A 2024, 128, 4369-4377, which demonstrated Doktorov's quantum algorithm as a proof of concept to predict the vibronic structure of ionization spectra, it is applied here to medium-size molecules with more than 30 vibrational normal modes, without accounting for Duschinsky rotations due to software limitations. This application to simulate the vibronic structures of the spectra of phenol and phenolate also enables assessing the impact of the differences in vibrational frequencies between the ground and excited electronic states.

Understanding and mitigating noise in molecular quantum linear response for spectroscopic properties on quantum computers

The promise of quantum computing to circumvent the exponential scaling of quantum chemistry has sparked a race to develop chemistry algorithms for quantum architecture. However, most works neglect the quantum-inherent shot noise, let alone the effect of current noisy devices. Here, we present a comprehensive study of quantum linear response (qLR) theory obtaining spectroscopic properties on simulated fault-tolerant quantum computers and present-day near-term quantum hardware. This work introduces novel metrics to analyze and predict the origins of noise in the quantum algorithm, proposes an Ansatz-based error mitigation technique, and reveals the significant impact of Pauli saving in reducing measurement costs and noise in subspace methods. Our hardware results using up to cc-pVTZ basis set serve as proof of principle for obtaining absorption spectra on quantum hardware in a general approach with the accuracy of classical multi-configurational methods. Importantly, our results exemplify that substantial improvements in hardware error rates and measurement speed are necessary to lift quantum computational chemistry from proof of concept to an actual impact in the field.

Divergences in classical and quantum linear response and equation of motion formulations

Calculating molecular properties using quantum devices can be performed through the quantum linear response (qLR) or, equivalently, the quantum equation of motion (qEOM) formulations. Different parameterizations of qLR and qEOM are available, namely naïve, projected, self-consistent, and state-transfer. In the naïve and projected parameterizations, the metric is not the identity, and we show that it depends on redundant orbital rotations. This dependency may lead to divergences in the excitation energies for certain choices of the redundant orbital rotation parameters in an idealized noiseless setting. Furthermore, this leads to a significant variance when calculations include statistical noise from finite quantum sampling.

Exact Ansatz of Fermion-Boson Systems for a Quantum Device

We present an exact Ansatz for the eigenstate problem of mixed fermion-boson systems that can be implemented on quantum devices. Based on a generalization of the electronic contracted Schrödinger equation (CSE), our approach guides a trial wave function to the ground state of any arbitrary mixed Hamiltonian by directly measuring residuals of the mixed CSE on a quantum device. Unlike density functional and coupled cluster theories applied to electron-phonon or electron-photon systems, the accuracy of our approach is not limited by the unknown exchange-correlation functional or the uncontrolled form of the exponential Ansatz. To test the performance of the method, we study the Tavis-Cummings model, commonly used in polaritonic quantum chemistry. Our results demonstrate that the CSE is a powerful tool in the development of quantum algorithms for solving general fermion-boson many-body problems.

Toward Accurate Calculation of Excitation Energies on Quantum Computers with ΔADAPT-VQE: A Case Study of BODIPY Derivatives

Quantum chemistry simulations offer a cost-effective way to computationally design BODIPY photosensitizers. However, accurate predictions of excitation energies pose a challenge for time-dependent density functional theory and equation-of-motion coupled-cluster singles and doubles methods. By contrast, reliable predictions can be achieved by multireference quantum chemistry methods; unfortunately, their computational cost increases exponentially with the number of electrons. Alternatively, quantum computing holds potential for an exact simulation of the photophysical properties in a computationally more efficient way. Herein, we introduce the state-specific ΔUCCSD-VQE (unitary coupled-cluster singles and doubles-variational quantum eigensolver) and ΔADAPT-VQE methods in which the electronically excited state is calculated via a non-Aufbau configuration. We show for six BODIPY derivatives that the proposed methods predict accurate excitation energies that are in good agreement with those from experiments. Due to its performance and simplicity, we believe that ΔADAPT will become a useful approach for the simulation of BODIPY photosensitizers on near-term quantum devices.

Simulating Vibronic Spectra by Direct Application of Doktorov Formulas on a Superconducting Quantum Simulator

In this work, a direct quantum implementation of the Doktorov formulas for calculating the vibronic spectrum of molecules under the harmonic approximation is presented. It is applied to the three-atom molecules H2O, SO2, ClO2, HS2, and ZnOH. The method solves the classically hard problem of estimating the Franck-Condon (FC) factors by using the Duschinsky matrices as the only input via the Doktorov quantum circuit. This has the advantage of avoiding basis changes, artificial squeezing parameters, and symmetry dependencies. In other words, it is a general method for three-atom molecules that can easily be generalized to bigger molecules. The results are compared with other quantum algorithms and classical anharmonic algorithms. Furthermore, the circuit requirements are studied in order to estimate its applicability on real superconducting quantum hardware.

Quantum Equation of Motion with Orbital Optimization for Computing Molecular Properties in Near-Term Quantum Computing

Determining the properties of molecules and materials is one of the premier applications of quantum computing. A major question in the field is how to use imperfect near-term quantum computers to solve problems of practical value. Inspired by the recently developed variants of the quantum counterpart of the equation-of-motion (qEOM) approach and the orbital-optimized variational quantum eigensolver (oo-VQE), we present a quantum algorithm (oo-VQE-qEOM) for the calculation of molecular properties by computing expectation values on a quantum computer. We perform noise-free quantum simulations of BeH2 in the series of STO-3G/6-31G/6-31G* basis sets and of H4 and H2O in 6-31G using an active space of four electrons and four spatial orbitals (8 qubits) to evaluate excitation energies, electronic absorption, and, for twisted H4, circular dichroism spectra. We demonstrate that the proposed algorithm can reproduce the results of conventional classical CASSCF calculations for these molecular systems.

Simulating Polaritonic Ground States on Noisy Quantum Devices

The recent advent of quantum algorithms for noisy quantum devices offers a new route toward simulating strong light-matter interactions of molecules in optical cavities for polaritonic chemistry. In this work, we introduce a general framework for simulating electron-photon-coupled systems on small, noisy quantum devices. This method is based on the variational quantum eigensolver (VQE) with the polaritonic unitary coupled cluster (PUCC) ansatz. To achieve chemical accuracy, we exploit various symmetries in qubit reduction methods, such as electron-photon parity, and use recently developed error mitigation schemes, such as the reference zero-noise extrapolation method. We explore the robustness of the VQE-PUCC approach across a diverse set of regimes for the bond length, cavity frequency, and coupling strength of the H2 molecule in an optical cavity. To quantify the performance, we measure two properties: ground-state energy, fundamentally relevant to chemical reactivity, and photon number, an experimentally accessible general indicator of electron-photon correlation.

Variational Quantum Eigensolver Boosted by Adiabatic Connection

In this work, we integrate the variational quantum eigensolver (VQE) with the adiabatic connection (AC) method for efficient simulations of chemical problems on near-term quantum computers. Orbital-optimized VQE methods are employed to capture the strong correlation within an active space, and classical AC corrections recover the dynamical correlation effects comprising electrons outside of the active space. On two challenging strongly correlated problems, namely, the dissociation of N2 and the electronic structure of the tetramethyleneethane biradical, we show that the combined VQE-AC approach enhances the performance of VQE dramatically. Moreover, since the AC corrections do not bring any additional requirements on quantum resources or measurements, they can actually boost the VQE algorithms. Our work paves the way toward quantum simulations of real-life problems on near-term quantum computers.

Toward Accurate Post-Born-Oppenheimer Molecular Simulations on Quantum Computers: An Adaptive Variational Eigensolver with Nuclear-Electronic Frozen Natural Orbitals

Nuclear quantum effects such as zero-point energy and hydrogen tunneling play a central role in many biological and chemical processes. The nuclear-electronic orbital (NEO) approach captures these effects by treating selected nuclei quantum mechanically on the same footing as electrons. On classical computers, the resources required for an exact solution of NEO-based models grow exponentially with system size. By contrast, quantum computers offer a means of solving this problem with polynomial scaling. However, due to the limitations of current quantum devices, NEO simulations are confined to the smallest systems described by minimal basis sets, whereas realistic simulations beyond the Born-Oppenheimer approximation require more sophisticated basis sets. For this purpose, we herein extend a hardware-efficient ADAPT-VQE method to the NEO framework in the frozen natural orbital (FNO) basis. We demonstrate on H2 and D2 molecules that the NEO-FNO-ADAPT-VQE method reduces the CNOT count by several orders of magnitude relative to the NEO unitary coupled cluster method with singles and doubles while maintaining the desired accuracy. This extreme reduction in the CNOT gate count is sufficient to permit practical computations employing the NEO method─an important step toward accurate simulations involving nonclassical nuclei and non-Born-Oppenheimer effects on near-term quantum devices. We further show that the method can capture isotope effects, and we demonstrate that inclusion of correlation energy systematically improves the prediction of difference in the zero-point energy (ΔZPE) between isotopes.

Unitary coupled-cluster for quantum computation of molecular properties in a strong magnetic field

In truncated coupled-cluster (CC) theories, non-variational and/or generally complex ground-state energies can occur. This is due to the non-Hermitian nature of the similarity transformed Hamiltonian matrix in combination with CC truncation. For chemical problems that deal with real-valued Hamiltonian matrices, complex CC energies rarely occur. However, for complex-valued Hamiltonian matrices, such as those that arise in the presence of strong magnetic fields, complex CC energies can be regularly observed unless certain symmetry conditions are fulfilled. Therefore, in the presence of magnetic fields, it is desirable to pursue CC methods that are guaranteed to give upper-bound, real-valued energies. In this work, we present the first application of unitary CC to chemical systems in a strong magnetic field. This is achieved utilizing the variational quantum eigensolver algorithm applied to the unitary coupled-cluster singles and doubles (UCCSD) method. We benchmark the method on the H2 molecule in a strong magnetic field and then calculate UCCSD energies for the H4 molecule as a function of both geometry and field angle. We show that while standard CCSD can yield generally complex energies that are not an upper-bound to the true energy, UCCSD always results in variational and real-valued energies. We also show that the imaginary components of the CCSD energy are largest in the strongly correlated region. Last, the UCCSD calculations capture a large percentage of the correlation energy.