We introduce the orbital-optimized multistate contracted variational quantum eigensolver (oo-MC-VQE) method with spin-adapted operators for the computation of ground and excited states, as well as state-specific and transition properties. The use of spin-adapted operators ensures that the spin symmetry of the reference states is conserved throughout the VQE optimization. In multistate variational approaches, achieving a balanced description of an increasing number of electronic states places growing demands on the expressibility of the underlying ansatz, thereby introducing a fundamental trade-off between accuracy and circuit complexity. We consider the effects of this trade-off explicitly and find that the number of circuit parameters required to obtain accurate results is reported to scale approximately linearly in the number of states. We further present an explicit quantum-circuit implementation of the oo-MC-VQE method and demonstrate its integration with quantum error mitigation techniques. Finally, we execute the method on real quantum devices to compute absorption spectra for two benchmark molecular systems.
High Resolution Image Download MS PowerPoint Slide The recent developments of quantum computing present novel potential pathways for quantum chemistry as the scaling of the computational power of quantum computers could be harnessed to naturally encode and solve electronic structure problems. Theoretically exact quantum algorithms for chemistry have been proposed (e.g., quantum phase estimation), but the limited capabilities of current noisy intermediate-scale quantum devices motivated the development of less demanding hybrid algorithms. In this context, the variational quantum eigensolver (VQE) algorithm was successfully introduced as an effective method to compute the ground-state energies of small molecules. This study investigates the folded spectrum (FS) method as an extension of the VQE algorithm for the computation of molecular excited states. It provides the possibility of directly computing excited states around a selected target energy using the same or quantum circuit as for the ground-state calculation. Inspired by the variance-based methods from the quantum Monte Carlo literature, the FS method minimizes the energy variance, thus, in principle, requiring a computationally expensive squared Hamiltonian to be applied. We alleviate this potentially poor scaling by employing a Pauli grouping procedure to identify sets of commuting Pauli strings that can be evaluated simultaneously. This allows for a significant reduction in the computational cost. We applied the FS-VQE method to small molecules (H 2, LiH), obtaining all electronic excited states with chemical accuracy on ideal quantum simulators. Furthermore, we explore the application of quantum error mitigation techniques, demonstrating improved energy accuracy on noisy simulators compared with simulations without mitigation.
Calculations of molecular spectral properties, like photodissociation rates and absorption bands, rely on knowledge of the excited state energies of the molecule of interest. Protocols based on the variational quantum eigensolver (VQE) are promising candidates to calculate such energies on emerging noisy intermediate-scale quantum (NISQ) computers. The successful implementation of these protocols on NISQ computers, relies on ans\"atze that can accurately approximate the molecular states and that can be implemented by shallow quantum circuits. We introduce the excited qubit-excitation-based adaptive (e-QEB-ADAPT)-VQE protocol to calculate molecular-excited-state energies. The e-QEB-ADAPT-VQE constructs efficient problem-tailored ans\"atze by iteratively appending evolutions of qubit excitation operators. The e-QEB-ADAPT-VQE also improves on previous ADAPT-VQE protocol in that it is designed to be independent on the choice of initial reference state. We perform classical numerical simulations for LiH and
BeH2 to benchmark the performance of the e-QEB-ADAPT-VQE. We demonstrate that the e-QEB-ADAPT-VQE can construct highly accurate ans\"atze that require at least an order of magnitude fewer cnot gates than standard fixed unitary coupled-cluster ans\"atze, such as the UCCSD and the GUCCSD. We also show that the e-QEB-ADAPT-VQE is more successful in constructing ans\"atze for excited molecular states than other ADAPT-VQE protocols.
Variational quantum eigensolver (VQE) is an appealing candidate for the application of near-term quantum computers. A technique introduced in [Higgot et al., Quantum 3, 156 (2019)], which is named variational quantum deflation (VQD), has extended the ability of the VQE framework for finding excited states of a Hamiltonian. However, no method to evaluate transition amplitudes between the eigenstates found by the VQD without using any costly Hadamard-test-like circuit has been proposed despite its importance for computing properties of the system such as oscillator strengths of molecules. Here we propose a method to evaluate transition amplitudes between the eigenstates obtained by the VQD avoiding any Hadamard-test-like circuit. Our method relies only on the ability to estimate overlap between two states, so it does not restrict to the VQD eigenstates and applies for general situations. To support the significance of our method, we provide a comprehensive comparison of three previously proposed methods to find excited states with numerical simulation of three molecules (lithium hydride, diazene, and azobenzene) in a noiseless situation and find that the VQD method exhibits the best performance among the three methods. Finally, we demonstrate the validity of our method by calculating the oscillator strength of lithium hydride in numerical simulations with shot noise. Our results illustrate the superiority of the VQD to find excited states and widen its applicability to various quantum systems.
Within the evolving domain of quantum computational chemistry, the Variational Quantum Eigensolver (VQE) has been developed to explore not only the ground state but also the excited states of molecules. In this study, we compare the performance of Variational Quantum Deflation (VQD) and Subspace-Search Variational Quantum Eigensolver (SSVQE) methods in determining the low-lying excited states of
LiH. Our investigation reveals that while VQD exhibits a slight advantage in accuracy, SSVQE stands out for its efficiency, allowing the determination of all low-lying excited states through a single parameter optimization procedure. We further evaluate the effectiveness of optimizers, including Gradient Descent (GD), Quantum Natural Gradient (QNG), and Adam optimizer, in obtaining
LiH's first excited state, with the Adam optimizer demonstrating superior efficiency in requiring the fewest iterations. Moreover, we propose a novel approach combining Folded Spectrum VQE (FS-VQE) with either VQD or SSVQE, enabling the exploration of highly excited states. We test the new approaches for finding all three
H4's excited states. Folded Spectrum SSVQE (FS-SSVQE) can find all three highly excited states near
−1.0 Ha with only one optimizing procedure, but the procedure converges slowly. In contrast, although Folded spectrum VQD (FS-VQD) gets highly excited states with individual optimizing procedures, the optimizing procedure converges faster.
Calculations of molecular spectral properties, like photodissociation rates and absorption bands, rely on knowledge of the excited state energies of the molecule of interest. Protocols based on the variational quantum eigensolver (VQE) are promising candidates to calculate such energies on emerging noisy intermediate scale quantum (NISQ) computers. The successful implementation of these protocols on NISQ computers, relies on ansätze that can accurately approximate the molecular states and that can be implemented by shallow quantum circuits. In this paper, we introduce the excited qubit-excitation-based adaptive (e-QEB-ADAPT)-VQE protocol to calculate molecular excited state energies. The e-QEB-ADAPT-VQE constructs efficient problem-tailored ansätze by iteratively appending evolutions of qubit excitation operators. The e-QEB-ADAPT-VQE is an adaptation of the QEB-ADAPT-VQE protocol, which is designed to be independent on the choice of an initial reference state. We perform classical numerical simulations for LiH and BeH
2 to benchmark the performance of the e-QEB-ADAPT-VQE. We demonstrate that the e-QEB-ADAPT-VQE can construct highly accurate ansätze that require at least an order of magnitude fewer
CNOTs than standard fixed UCC ansätze, such as the UCCSD and the GUCCSD.