David Muñoz Ramo

We develop a workflow within the ONIOM framework and demonstrate it on the hybrid computing system consisting of the supercomputer Fugaku and the Quantinuum Reimei trapped-ion quantum computer. This hybrid platform extends the layered approach for biomolecular chemical reactions to accurately treat the active site, such as a protein, and the large and often weakly correlated molecular environment. Our result marks a significant milestone in enabling scalable and accurate simulation of complex biomolecular reactions

Emergent quantum computing technologies are widely expected to provide novel approaches in the simulation of quantum chemistry. Despite rapid improvements in the scale and fidelity of quantum computers, high resource requirements make the execution of quantum chemistry experiments challenging. Typical experiments are limited in the number of qubits used, incur a substantial shot cost, or require complex architecture-specific optimization and error mitigation techniques. In this paper, we propose a conceptually simple benchmarking approach involving the use of multi-ancilla quantum phase estimation. Our approach is restricted to very small chemical systems, and does not scale favorably beyond molecular systems that can be described with $2$ qubits; however, this restriction allows us to generate circuits that scale quadratically in gate count with the number of qubits in the readout register. This enables the execution of quantum chemistry circuits that act on many qubits, while producing meaningful results with limited shot counts. We use this technique (with $200$ shots per experiment) to calculate the ground state energy of molecular hydrogen to $50$ bits of precision ($8.9 \times 10^{-16}$ hartree) on a $56$-qubit trapped-ion quantum computer, negating Trotter error. Including Trotter error, we obtain between $32$ and $36$ bits of precision ($1.5 * 10^{-10}$ and $6.0 * 10^{-11}$ hartree respectively), vastly exceeding chemical accuracy ($1.6 * 10^{-3}$ hartree) against Full Configuration Interaction. We consider application of the approach to deeper circuits, and discuss potential as a benchmark task for near-term quantum devices.

The linear combination of unitaries (LCU) method has proven to scale better than existing product formulas in simulating long time Hamiltonian dynamics. However, given the number of multi-control gate operations in the standard prepare-select-unprepare architecture of LCU, it is still resource-intensive to implement on the current quantum computers. In this work, we demonstrate LCU implementations on an ion trap quantum computer for calculating squared overlaps $|\langle ψ(t=0)|ψ(t>0)\rangle|^2$ of time-evolved states. This is achieved by an optimized LCU method, based on pre-selecting relevant unitaries, coupled with a compilation strategy which makes use of quantum multiplexor gates, leading to a significant reduction in the depth and number of two-qubit gates in circuits. For $L$ Pauli strings in a Taylor series expanded $n$-qubit-mapped time evolution operator, we find a two-qubit gate count of $2^{\lceil log_2(L)\rceil}(2n+1)-n-2$. We test this approach by simulating a Rabi-Hubbard Hamiltonian.

Quantifying correlation and entanglement between molecular orbitals can elucidate the role of quantum effects in strongly correlated reaction processes. However, accurately storing the wavefunction for a classical computation of those quantities can be prohibitive. Here we use the Quantinuum H1-1 trapped-ion quantum computer to calculate von Neumann entropies which quantify the orbital correlation and entanglement in a strongly correlated molecular system relevant to lithium-ion batteries (vinylene carbonate interacting with an O$_2$ molecule). As shown in previous works, fermionic superselection rules decrease correlations and reduce measurement overheads for constructing orbital reduced density matrices. Taking into account superselection rules we further reduce the number of measurements by finding commuting sets of Pauli operators. Using low overhead noise reduction techniques we calculate von Neumann entropies in excellent agreement with noiseless benchmarks, indicating that correlations and entanglement between molecular orbitals can be accurately estimated from a quantum computation. Our results show that the one-orbital entanglement vanishes unless opposite-spin open shell configurations are present in the wavefunction.

In this paper, we present a quantum computational method to calculate the many-body Green's function matrix in a spin orbital basis. We apply our approach to finite-sized fermionic Hubbard models and related impurity models within Dynamical Mean Field Theory, and demonstrate the calculation of Green's functions on Quantinuum's H1-1 trapped-ion quantum computer. Our approach involves a cumulant expansion of the Lanczos method, using Hamiltonian moments as measurable expectation values. This bypasses the need for a large overhead in the number of measurements due to repeated applications of the variational quantum eigensolver (VQE), and instead measures the expectation value of the moments with one set of measurement circuits. From the measured moments, the tridiagonalised Hamiltonian matrix can be computed, which in turn yields the Green's function via continued fractions. While we use a variational algorithm to prepare the ground state in this work, we note that the modularity of our implementation allows for other (non-variational) approaches to be used for the ground state.

Hydrogen has emerged as a promising energy source, holding the key to achieve low-carbon and sustainable mobility. However, its applications are still limited by modest conversion efficiency in the electrocatalytic oxygen reduction reaction (ORR) within fuel cells. Consequently, the development of novel catalysts and a profound understanding of the underlying reactions have become of paramount importance. The complex nature of the ORR potential energy landscape and the presence of strong electronic correlations present challenges to atomistic modelling using classical computers. This scenario opens new avenues for the implementation of novel quantum computing workflows to address these molecular systems. Here, we present a pioneering study that combines classical and quantum computational approaches to investigate the ORR on pure platinum and platinum/cobalt surfaces. Our research demonstrates, for the first time, the feasibility of implementing this workflow on the H1-series trapped-ion quantum computer and identify the challenges of the quantum chemistry modelling of this reaction. The results highlight the involvement of strongly correlated species in the cobalt-containing catalyst, suggesting their potential as ideal candidates for showcasing quantum advantage in future applications.

Quantum phase estimation (QPE) serves as a building block of many different quantum algorithms and finds important applications in computational chemistry problems. Despite the rapid development of quantum hardware, experimental demonstration of QPE for chemistry problems remains challenging due to its large circuit depth and the lack of quantum resources to protect the hardware from noise with fully fault-tolerant protocols. In the present work, we take a step towards fault-tolerant quantum computing by demonstrating a QPE algorithm on a Quantinuum trapped-ion computer. We employ a Bayesian approach to QPE and introduce a routine for optimal parameter selection, which we combine with a $[[ n+2,n,2 ]]$ quantum error detection code carefully tailored to the hardware capabilities. As a simple quantum chemistry example, we take a hydrogen molecule represented by a two-qubit Hamiltonian and estimate its ground state energy using our QPE protocol. In the experiment, we use the quantum circuits containing as many as 920 physical two-qubit gates to estimate the ground state energy within $6\times 10^{-3}$ hartree of the exact value.

We have demonstrated a prototypical hybrid classical and quantum computational workflow for the quantification of protein-ligand interactions. The workflow combines the Density Matrix Embedding Theory (DMET) embedding procedure with the Variational Quantum Eigensolver (VQE) approach for finding molecular electronic ground states. A series of $β$-secretase (BACE1) inhibitors is rank-ordered using binding energy differences calculated on the latest superconducting transmon (IBM) and trapped-ion (Honeywell) Noisy Intermediate Scale Quantum (NISQ) devices. This is the first application of real quantum computers to the calculation of protein-ligand binding energies. The results shed light on hardware and software requirements which would enable the application of NISQ algorithms in drug design.