Determining the spectral structure of molecular Hamiltonians is one of the most important and challenging research fields in quantum chemistry. In the near future, such a problem may be solvable by means of quantum simulation. However, constrained by a handful of qubits and noisy quantum gates available in the near-term quantum devices, hybrid quantum-classical algorithms based on variational methods such as variational quantum eigensolver have become a preferable choice. Here we propose an interesting variational quantum algorithm for obtaining the excited states of many-body Hamiltonians. To obtain the kth excited states, the previous variational quantum deflation algorithm requires the use of k different variational circuits, which is now reduced to only two circuits with our algorithm. The feasibility of the algorithm was numerically tested on the hydrogen molecule with and without random noise, where the chemical accuracy is reached for both cases. This algorithm is expected to significantly reduce the resource requirements of quantum simulation.
Determining the spectral structure of molecular Hamiltonians is one of the most important and challenging research fields in quantum chemistry. In the near future, such a problem may be solvable by means of quantum simulation. However, constrained by a handful of qubits and noisy quantum gates available in the near-term quantum devices, hybrid quantum-classical algorithms based on variational methods such as variational quantum eigensolver have become a preferable choice. Here we propose an interesting variational quantum algorithm for obtaining the excited states of many-body Hamiltonians. To obtain the kth excited states, the previous variational quantum deflation algorithm requires the use of k different variational circuits, which is now reduced to only two circuits with our algorithm. The feasibility of the algorithm was numerically tested on the hydrogen molecule with and without random noise, where the chemical accuracy is reached for both cases. This algorithm is expected to significantly reduce the resource requirements of quantum simulation.
Utilizing quantum computer to investigate quantum chemistry is an important research field nowadays. In addition to the ground-state problems that have been widely studied, the determination of excited-states plays a crucial role in the prediction and modeling of chemical reactions and other physical processes. Here, we propose a non-variational full circuit-based quantum algorithm for obtaining the excited-state spectrum of a quantum chemistry Hamiltonian. Compared with previous classical-quantum hybrid variational algorithms, our method eliminates the classical optimization process, reduces the resource cost caused by the interaction between different systems, and achieves faster convergence rate and stronger robustness against noise without barren plateau. The parameter updating for determining the next energy-level is naturally dependent on the energy measurement outputs of the previous energy-level and can be realized by only modifying the state preparation process of ancillary system, introducing little additional resource overhead. Numerical simulations of the algorithm with hydrogen, LiH, H2O and NH3 molecules are presented. Furthermore, we offer an experimental demonstration of the algorithm on a superconducting quantum computing platform, and the results show a good agreement with theoretical expectations. The algorithm can be widely applied to various Hamiltonian spectrum determination problems on the fault-tolerant quantum computers.
Variational quantum algorithms on bosonic quantum processors are an emerging paradigm for quantum chemistry calculations, exploiting the natural alignment between molecular structure and harmonic oscillator-based hardware. We introduce the qumode-based variational quantum deflation framework (QumVQD) for finding both electronic and vibrational excited state energies on qumode-based architectures. We validate the approach through electronic structure calculations on H2 and linear H4, where we introduce Hamming-weight filtering of the Fock basis to enforce particle number conservation and eliminate spurious eigenstates by reducing the required Hilbert space, which reduces the required number of qumodes in turn. We achieve agreement with full configuration interaction (FCI) using the STO-3G basis set within the chemical accuracy threshold at most points along the potential energy surfaces. Extending to the vibrational structure, we combine QumVQD with an existing Hamiltonian fragmentation approach based on Cartan subalgebra, allowing us to compute the vibrational eigenenergies of CO2 and H2S to spectroscopic accuracy with per-fragment circuits that scale as O(N) in single-qumode gates and O(N2) in beam-splitter gates for N qumodes. For the case of CO2, we get total gate counts more than an order of magnitude smaller than those reported for qubit-based vibrational algorithms at this system size. These results demonstrate that bosonic quantum devices are a viable platform for excited-state quantum chemistry, particularly for vibrational problems where qubit-based methods incur substantial boson-to-qubit mapping overhead.
Variational quantum algorithms are emerging as promising candidates for near-term practical applications of quantum information processors, in the field of quantum chemistry. We implement the variational quantum eigensolver algorithm to calculate the molecular ground-state energy of the H2 molecule and experimentally demonstrated it on an NMR quantum processor. Further, we simulate the excited states of the H2 molecule using the variational quantum deflation algorithm and experimentally demonstrate it on the same NMR quantum processor. We also develop the first simulation of the energy calculation of the H2 molecule using only a single qubit, and verify the results on an NMR quantum computer. Our experimental results demonstrate that only a single NMR qubit suffices to calculate the molecular energies of the H2 molecule to the desired accuracy.
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.
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 electronic 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 comparison of three previously proposed methods to find excited states with numerical simulation of three molecules (lithium hydride, diazene, and azobenzene) and find that the VQD method exhibits the best performance among the three methods within the scope of noiseless simulation. Finally, we demonstrate the validity of our method by calculating the oscillator strength of lithium hydride, comparing results from numerical simulations and real-hardware experiments on the cloud enabled quantum computer IBMQ Rome. Our results illustrate the superiority of the VQD to find excited states and widen its applicability to various quantum systems.
Variational Quantum Algorithms (VQAs) are iterative algorithms suited to implementation on current-era quantum devices. VQAs employ classical optimization to minimize cost functions evaluated on quantum circuits. However, the extent to which VQAs manage noise is often insufficient for quantum chemistry applications. One method of improving VQAs is through accurate ansatz initialization. The CAFQA (Clifford Ansatz For Quantum Accuracy) protocol runs a discrete search through a classically simulatable subset of the entire state space to find a desirable initialization. Prior work has evaluated CAFQA applied to the Variational Quantum Eigensolver (VQE), a VQA that computes grounds states of a Hamiltonian. Motivated by CAFQA's success, we propose Excited-CAFQA initialization for Variational Quantum Deflation (VQD), a quantum algorithm that extends VQE by allowing the computation of excited states. VQD recursively computes excited states, by constraining the kth state to be orthogonal to the previous k-1 computed energy states via a penalty term appended to the standard VQE cost function. Just as with VQE, the VQD cost function can be efficiently computed classically for the states considered in the discrete CAFQA search, allowing for the discrete CAFQA optimizer to find good initial parameters for each energy level computation. Preliminary evaluation shows that Excited-CAFQA achieves 90 to 99+% accuracy across a variety of bond lengths and excited states for H2 and HeH+ molecular systems.