Exploiting Locality in Quantum Computation for Quantum Chemistry

Folded Spectrum VQE: A Quantum Computing Method for the Calculation of Molecular Excited States

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 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 (H2, 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.

Online quantum time series processing with random oscillator networks

Reservoir computing is a powerful machine learning paradigm for online time series processing. It has reached state-of-the-art performance in tasks such as chaotic time series prediction and continuous speech recognition thanks to its unique combination of high computational power and low training cost which sets it aside from alternatives such as traditionally trained recurrent neural networks, and furthermore is amenable to implementations in dedicated hardware, potentially leading to extremely compact and efficient reservoir computers. Recently the use of random quantum systems has been proposed, leveraging the complexity of quantum dynamics for classical time series processing. Extracting the output from a quantum system without disturbing its state too much is problematic however, and can be expected to become a bottleneck in such approaches. Here we propose a reservoir computing inspired approach to online processing of time series consisting of quantum information, sidestepping the measurement problem. We illustrate its power by generalizing two paradigmatic benchmark tasks from classical reservoir computing to quantum information and introducing a task without a classical analogue where a random system is trained to both create and distribute entanglement between systems that never directly interact. Finally, we discuss partial generalizations where only the input or only the output time series is quantum.

Evaluating the evidence for exponential quantum advantage in ground-state quantum chemistry

Due to intense interest in the potential applications of quantum computing, it is critical to understand the basis for potential exponential quantum advantage in quantum chemistry. Here we gather the evidence for this case in the most common task in quantum chemistry, namely, ground-state energy estimation, for generic chemical problems where heuristic quantum state preparation might be assumed to be efficient. The availability of exponential quantum advantage then centers on whether features of the physical problem that enable efficient heuristic quantum state preparation also enable efficient solution by classical heuristics. Through numerical studies of quantum state preparation and empirical complexity analysis (including the error scaling) of classical heuristics, in both ab initio and model Hamiltonian settings, we conclude that evidence for such an exponential advantage across chemical space has yet to be found. While quantum computers may still prove useful for ground-state quantum chemistry through polynomial speedups, it may be prudent to assume exponential speedups are not generically available for this problem.

Sub-system self-consistency in coupled cluster theory

In this article, we provide numerical evidence indicating that the single-reference coupled-cluster (CC) energies can be calculated alternatively to their copybook definition. We demonstrate that the CC energy can be reconstructed by diagonalizing the effective Hamiltonians describing correlated sub-systems of the many-body system. In the extreme case, we provide numerical evidence that the CC energy can be reproduced through the diagonalization of the effective Hamiltonian describing sub-system composed of a single electron. These properties of the CC formalism can be exploited to design protocols to define effective interactions in sub-systems used as probes to calculate the energy of the entire system and introduce a new type of self-consistency for approximate CC approaches.

Romero S, Oregi I, Osaba E, Villar-Rodriguez E, Ban Y. Digital Quantum Simulation and Circuit Learning for the Generation of Coherent States

Coherent states, known as displaced vacuum states, play an important role in quantum information processing, quantum machine learning, and quantum optics. In this article, two ways to digitally prepare coherent states in quantum circuits are introduced. First, we construct the displacement operator by decomposing it into Pauli matrices via ladder operators, i.e., creation and annihilation operators. The high fidelity of the digitally generated coherent states is verified compared with the Poissonian distribution in Fock space. Secondly, by using Variational Quantum Algorithms, we choose different ansatzes to generate coherent states. The quantum resources-such as numbers of quantum gates, layers and iterations-are analyzed for quantum circuit learning. The simulation results show that quantum circuit learning can provide high fidelity on learning coherent states by choosing appropriate ansatzes.

Amplitude Reordering Accelerates the Adaptive Variational Quantum Eigensolver Algorithms

The variational quantum eigensolver (VQE) algorithm can simulate the chemical systems such as molecules in the noisy-intermediate-scale quantum devices and shows promising applications in quantum chemistry simulations. The accuracy and computational cost of the VQE simulations are determined by the underlying ansatz. Therefore, the most important issue is to generate a compact and accurate ansatz, which requires a shallower parametric quantum circuit and can achieve an acceptable accuracy. The newly developed adaptive algorithms (AAs) such as the adaptive derivative-assembled pseudo-Trotter VQE (ADAPT-VQE) can solve this issue via generating compact and accurate ansatzes. However, these AAs show very low computational efficiency because they require a large number of additional measurements. Here we propose an amplitude reordering (AR) strategy to accelerate the promising but expensive AAs by adding operators in a "batched" fashion in a way that their order is still quasi-optimal. We first introduce the AR method into ADAPT-VQE and build the AR-ADAPT-VQE algorithm. We then endow the energy-sorting VQE (ES-VQE) algorithm with the adaptive feature and introduce the AR into AES-VQE to form the AR-AES-VQE algorithm. To demonstrate the performance of these algorithms, we calculate the dissociation curves of three small molecules, LiH, linear BeH₂, and linear H₆, by using (AR-)ADAPT-VQE and (AR-)AES-VQE algorithms. It is found that all of the AR-equipped AAs (AR-AAs) can significantly reduce the number of iterations and subsequently accelerate the calculations with a speedup of up to more than ten times without the obvious loss of accuracy. The final ansatz generated by the AR-AAs not only avoids extra circuit depth but also maintains the computational accuracy; sometimes the AR-AAs even outperforms their original counterparts.

Quantum Algorithm of the Divide-and-Conquer Unitary Coupled Cluster Method with a Variational Quantum Eigensolver

The variational quantum eigensolver (VQE) with shallow or constant-depth quantum circuits is one of the most pursued approaches in the noisy intermediate-scale quantum (NISQ) devices with incoherent errors. In this study, the divide-and-conquer (DC) linear scaling technique, which divides the entire system into several fragments, is applied to the VQE algorithm based on the unitary coupled cluster (UCC) method, denoted as DC-qUCC/VQE, to reduce the number of required qubits. The unitarity of the UCC ansatz that enables the evaluation of the total energy as well as various molecular properties as expectation values can be easily implemented on quantum devices because the quantum gates are unitary operators themselves. Based on this feature, the present DC-qUCC/VQE algorithm is designed to conserve the total number of electrons in the entire system using the density matrix evaluated on a quantum computer. Numerical assessments clarified that the energy errors of the DC-qUCC/VQE calculations decrease by using the constraint of the total number of electrons. Furthermore, the DC-qUCC/VQE algorithm could reduce the number of quantum gates and shows the possibility of decreasing incoherent errors.

Adiabatic state preparation of correlated wave functions with nonlinear scheduling functions and broken-symmetry wave functions

Adiabatic state preparation (ASP) can generate the correlated wave function by simulating the time evolution of wave function under the time-dependent Hamiltonian that interpolates the Fock operator and the full electronic Hamiltonian. However, ASP is inherently unsuitable for studying strongly correlated systems, and furthermore practical computational conditions for ASP are unknown. In quest for the suitable computational conditions for practical applications of ASP, we performed numerical simulations of ASP in the potential energy curves of N₂, BeH₂, and in the C_2v quasi-reaction pathway of the Be atom insertion to the H₂ molecule, examining the effect of nonlinear scheduling functions and the ASP with broken-symmetry wave functions with the S² operator as the penalty term, contributing to practical applications of quantum computing to quantum chemistry. Eventually, computational guidelines to generate the correlated wave functions having the square overlap with the complete-active space self-consistent field wave function close to unity are discussed.

Probing the Action of Screened Anticancer Triazole–Tetrazole Derivatives Against COVID-19 Using Molecular Docking and DFT Investigations

Drugs are continuously being evaluated for novel therapeutic uses. The purpose of this work was to screen anticancer triazole/tetrazole derivatives for effectiveness against the SARS-CoV-2 main protease (Mpro). First, the chemical structures’ activity was derived from conceptual quantum chemical calculations. According to molecular docking analysis, the compounds scored good interactions against SAR-COV-2's Mpro, with binding energies extending from −8.21 to −8.97 kcal/mol. The docked complexes included various bindings with His41 and Cys145, both catalytic residues responsible for cleavage of the SARS-CoV-2 Mpro. Among the 4 studied compounds, TD3 exhibited the highest affinity by achieving the most stable binding energy and lowest value for the inhibition constant. Most striking was that TD3 not only formed strong bonds with the catalytic residues His41 and Cys145, but also captured the residues of the catalytic loop (Cys44 to Pro52), which flank the catalytic dyads in Mpro's active site. As a result, the studied triazole/tetrazole derivatives, notably TD3, must be reviewed as potent drugs that could be repurposed for SARS-CoV-2 treatment.

What the foundations of quantum computer science teach us about chemistry

With the rapid development of quantum technology, one of the leading applications that has been identified is the simulation of chemistry. Interestingly, even before full scale quantum computers are available, quantum computer science has exhibited a remarkable string of results that directly impact what is possible in a chemical simulation with any computer. Some of these results even impact our understanding of chemistry in the real world. In this Perspective, we take the position that direct chemical simulation is best understood as a digital experiment. While on the one hand, this clarifies the power of quantum computers to extend our reach, it also shows us the limitations of taking such an approach too directly. Leveraging results that quantum computers cannot outpace the physical world, we build to the controversial stance that some chemical problems are best viewed as problems for which no algorithm can deliver their solution, in general, known in computer science as undecidable problems. This has implications for the predictive power of thermodynamic models and topics such as the ergodic hypothesis. However, we argue that this Perspective is not defeatist but rather helps shed light on the success of existing chemical models such as transition state theory, molecular orbital theory, and thermodynamics as models that benefit from data. We contextualize recent results, showing that data-augmented models are a more powerful rote simulation. These results help us appreciate the success of traditional chemical theory and anticipate new models learned from experimental data. Not only can quantum computers provide data for such models, but they can also extend the class and power of models that utilize data in fundamental ways. These discussions culminate in speculation on new ways for quantum computing and chemistry to interact and our perspective on the eventual roles of quantum computers in the future of chemistry.

Fermionic Partial Tomography via Classical Shadows

We propose a tomographic protocol for estimating any k-body reduced density matrix (k-RDM) of an n-mode fermionic state, a ubiquitous step in near-term quantum algorithms for simulating many-body physics, chemistry, and materials. Our approach extends the framework of classical shadows, a randomized approach to learning a collection of quantum-state properties, to the fermionic setting. Our sampling protocol uses randomized measurement settings generated by a discrete group of fermionic Gaussian unitaries, implementable with linear-depth circuits. We prove that estimating all k-RDM elements to additive precision ϵ requires on the order of (n/k)k^{3/2}log(n)/ϵ^{2} repeated state preparations, which is optimal up to the logarithmic factor. Furthermore, numerical calculations show that our protocol offers a substantial improvement in constant overheads for k≥2, as compared to prior deterministic strategies. We also adapt our method to particle-number symmetry, wherein the additional circuit depth may be halved at the cost of roughly 2-5 times more repetitions.

Quantum algorithm for alchemical optimization in material design

The development of tailored materials for specific applications is an active field of research in chemistry, material science and drug discovery. The number of possible molecules obtainable from a set of atomic species grow exponentially with the size of the system, limiting the efficiency of classical sampling algorithms. On the other hand, quantum computers can provide an efficient solution to the sampling of the chemical compound space for the optimization of a given molecular property. In this work, we propose a quantum algorithm for addressing the material design problem with a favourable scaling. The core of this approach is the representation of the space of candidate structures as a linear superposition of all possible atomic compositions. The corresponding ‘alchemical’ Hamiltonian drives the optimization in both the atomic and electronic spaces leading to the selection of the best fitting molecule, which optimizes a given property of the system, e.g., the interaction with an external potential as in drug design. The quantum advantage resides in the efficient calculation of the electronic structure properties together with the sampling of the exponentially large chemical compound space. We demonstrate both in simulations and with IBM Quantum hardware the efficiency of our scheme and highlight the results in a few test cases. This preliminary study can serve as a basis for the development of further material design quantum algorithms for near-term quantum computers.

‘Alchemical’ quantum algorithm for the simultaneous optimisation of chemical composition and electronic structure for material design. By exploiting quantum mechanical principles this approach will boost drug discovery in the near future.

A Full Quantum Eigensolver for Quantum Chemistry Simulations

Quantum simulation of quantum chemistry is one of the most compelling applications of quantum computing. It is of particular importance in areas ranging from materials science, biochemistry, and condensed matter physics. Here, we propose a full quantum eigensolver (FQE) algorithm to calculate the molecular ground energies and electronic structures using quantum gradient descent. Compared to existing classical-quantum hybrid methods such as variational quantum eigensolver (VQE), our method removes the classical optimizer and performs all the calculations on a quantum computer with faster convergence. The gradient descent iteration depth has a favorable complexity that is logarithmically dependent on the system size and inverse of the precision. Moreover, the FQE can be further simplified by exploiting a perturbation theory for the calculations of intermediate matrix elements and obtaining results with a precision that satisfies the requirement of chemistry application. The full quantum eigensolver can be implemented on a near-term quantum computer. With the rapid development of quantum computing hardware, the FQE provides an efficient and powerful tool to solve quantum chemistry problems.

Quantum Chemistry in the Age of Quantum Computing

Practical challenges in simulating quantum systems on classical computers have been widely recognized in the quantum physics and quantum chemistry communities over the past century. Although many approximation methods have been introduced, the complexity of quantum mechanics remains hard to appease. The advent of quantum computation brings new pathways to navigate this challenging and complex landscape. By manipulating quantum states of matter and taking advantage of their unique features such as superposition and entanglement, quantum computers promise to efficiently deliver accurate results for many important problems in quantum chemistry, such as the electronic structure of molecules. In the past two decades, significant advances have been made in developing algorithms and physical hardware for quantum computing, heralding a revolution in simulation of quantum systems. This Review provides an overview of the algorithms and results that are relevant for quantum chemistry. The intended audience is both quantum chemists who seek to learn more about quantum computing and quantum computing researchers who would like to explore applications in quantum chemistry.

Digital quantum simulation of molecular vibrations

Molecular vibrations underpin important phenomena such as spectral properties, energy transfer, and molecular bonding. However, obtaining a detailed understanding of the vibrational structure of even small molecules is computationally expensive. While several algorithms exist for efficiently solving the electronic structure problem on a quantum computer, there has been comparatively little attention devoted to solving the vibrational structure problem with quantum hardware. In this work, we discuss the use of quantum algorithms for investigating both the static and dynamic vibrational properties of molecules. We introduce a physically motivated unitary vibrational coupled cluster ansatz, which also makes our method accessible to noisy, near-term quantum hardware. We numerically test our proposals for the water and sulfur dioxide molecules.

Quantum Chemistry on Quantum Computers: A Method for Preparation of Multiconfigurational Wave Functions on Quantum Computers without Performing Post-Hartree-Fock Calculations

The full configuration interaction (full-CI) method is capable of providing the numerically best wave functions and energies of atoms and molecules within basis sets being used, although it is intractable for classical computers. Quantum computers can perform full-CI calculations in polynomial time against the system size by adopting a quantum phase estimation algorithm (QPEA). In the QPEA, the preparation of initial guess wave functions having sufficiently large overlap with the exact wave function is recommended. The Hartree-Fock (HF) wave function is a good initial guess only for closed shell singlet molecules and high-spin molecules carrying no spin-β unpaired electrons, around their equilibrium geometry, and thus, the construction of multiconfigurational wave functions without performing post-HF calculations on classical computers is highly desired for applying the method to a wide variety of chemistries and physics. In this work, we propose a method to construct multiconfigurational initial guess wave functions suitable for QPEA-based full-CI calculations on quantum computers, by utilizing diradical characters computed from spin-projected UHF wave functions. The proposed approach drastically improves the wave function overlap, particularly in molecules with intermediate diradical characters.

Barren plateaus in quantum neural network training landscapes

Many experimental proposals for noisy intermediate scale quantum devices involve training a parameterized quantum circuit with a classical optimization loop. Such hybrid quantum-classical algorithms are popular for applications in quantum simulation, optimization, and machine learning. Due to its simplicity and hardware efficiency, random circuits are often proposed as initial guesses for exploring the space of quantum states. We show that the exponential dimension of Hilbert space and the gradient estimation complexity make this choice unsuitable for hybrid quantum-classical algorithms run on more than a few qubits. Specifically, we show that for a wide class of reasonable parameterized quantum circuits, the probability that the gradient along any reasonable direction is non-zero to some fixed precision is exponentially small as a function of the number of qubits. We argue that this is related to the 2-design characteristic of random circuits, and that solutions to this problem must be studied.

A Comparison of the Bravyi-Kitaev and Jordan-Wigner Transformations for the Quantum Simulation of Quantum Chemistry

The ability to perform classically intractable electronic structure calculations is often cited as one of the principal applications of quantum computing. A great deal of theoretical algorithmic development has been performed in support of this goal. Most techniques require a scheme for mapping electronic states and operations to states of and operations upon qubits. The two most commonly used techniques for this are the Jordan-Wigner transformation and the Bravyi-Kitaev transformation. However, comparisons of these schemes have previously been limited to individual small molecules. In this paper, we discuss resource implications for the use of the Bravyi-Kitaev mapping scheme, specifically with regard to the number of quantum gates required for simulation. We consider both small systems, which may be simulatable on near-future quantum devices, and systems sufficiently large for classical simulation to be intractable. We use 86 molecular systems to demonstrate that the use of the Bravyi-Kitaev transformation is typically at least approximately as efficient as the canonical Jordan-Wigner transformation and results in substantially reduced gate count estimates when performing limited circuit optimizations.

Quantum information processing with superconducting circuits: a review

During the last ten years, superconducting circuits have passed from being interesting physical devices to becoming contenders for near-future useful and scalable quantum information processing (QIP). Advanced quantum simulation experiments have been shown with up to nine qubits, while a demonstration of quantum supremacy with fifty qubits is anticipated in just a few years. Quantum supremacy means that the quantum system can no longer be simulated by the most powerful classical supercomputers. Integrated classical-quantum computing systems are already emerging that can be used for software development and experimentation, even via web interfaces. Therefore, the time is ripe for describing some of the recent development of superconducting devices, systems and applications. As such, the discussion of superconducting qubits and circuits is limited to devices that are proven useful for current or near future applications. Consequently, the centre of interest is the practical applications of QIP, such as computation and simulation in Physics and Chemistry.

Elucidating reaction mechanisms on quantum computers

With rapid recent advances in quantum technology, we are close to the threshold of quantum devices whose computational powers can exceed those of classical supercomputers. Here, we show that a quantum computer can be used to elucidate reaction mechanisms in complex chemical systems, using the open problem of biological nitrogen fixation in nitrogenase as an example. We discuss how quantum computers can augment classical computer simulations used to probe these reaction mechanisms, to significantly increase their accuracy and enable hitherto intractable simulations. Our resource estimates show that, even when taking into account the substantial overhead of quantum error correction, and the need to compile into discrete gate sets, the necessary computations can be performed in reasonable time on small quantum computers. Our results demonstrate that quantum computers will be able to tackle important problems in chemistry without requiring exorbitant resources.

Efficient Bayesian Phase Estimation

We introduce a new method called rejection filtering that we use to perform adaptive Bayesian phase estimation. Our approach has several advantages: it is classically efficient, easy to implement, achieves Heisenberg limited scaling, resists depolarizing noise, tracks time-dependent eigenstates, recovers from failures, and can be run on a field programmable gate array. It also outperforms existing iterative phase estimation algorithms such as Kitaev's method.

Error Sensitivity to Environmental Noise in Quantum Circuits for Chemical State Preparation

Calculating molecular energies is likely to be one of the first useful applications to achieve quantum supremacy, performing faster on a quantum than a classical computer. However, if future quantum devices are to produce accurate calculations, errors due to environmental noise and algorithmic approximations need to be characterized and reduced. In this study, we use the high performance qHiPSTER software to investigate the effects of environmental noise on the preparation of quantum chemistry states. We simulated 18 16-qubit quantum circuits under environmental noise, each corresponding to a unitary coupled cluster state preparation of a different molecule or molecular configuration. Additionally, we analyze the nature of simple gate errors in noise-free circuits of up to 40 qubits. We find that, in most cases, the Jordan-Wigner (JW) encoding produces smaller errors under a noisy environment as compared to the Bravyi-Kitaev (BK) encoding. For the JW encoding, pure dephasing noise is shown to produce substantially smaller errors than pure relaxation noise of the same magnitude. We report error trends in both molecular energy and electron particle number within a unitary coupled cluster state preparation scheme, against changes in nuclear charge, bond length, number of electrons, noise types, and noise magnitude. These trends may prove to be useful in making algorithmic and hardware-related choices for quantum simulation of molecular energies.

Adiabatic quantum simulation of quantum chemistry

We show how to apply the quantum adiabatic algorithm directly to the quantum computation of molecular properties. We describe a procedure to map electronic structure Hamiltonians to 2-body qubit Hamiltonians with a small set of physically realizable couplings. By combining the Bravyi-Kitaev construction to map fermions to qubits with perturbative gadgets to reduce the Hamiltonian to 2-body, we obtain precision requirements on the coupling strengths and a number of ancilla qubits that scale polynomially in the problem size. Hence our mapping is efficient. The required set of controllable interactions includes only two types of interaction beyond the Ising interactions required to apply the quantum adiabatic algorithm to combinatorial optimization problems. Our mapping may also be of interest to chemists directly as it defines a dictionary from electronic structure to spin Hamiltonians with physical interactions.