Quantum computation of reactions on surfaces using local embedding
Tanvi P. Gujarati, Mário Motta, Triet Friedhoff, Julia E. Rice, Nam Nguyen, Panagiotis Kl. Barkoutsos, Richard J. Thompson, Tyler Smith, Marna M. Kagele, Mark Brei, Barbara Jones, Kristen Williams
Abstract Modeling electronic systems is an important application for quantum computers. In the context of materials science, an important open problem is the computational description of chemical reactions on surfaces. In this work, we outline a workflow to model the adsorption and reaction of molecules on surfaces using quantum computing algorithms. We develop and compare two local embedding methods for the systematic determination of active spaces. These methods are automated and based on the physics of molecule-surface interactions and yield systematically improvable active spaces. Furthermore, to reduce the quantum resources required for the simulation of the selected active spaces using quantum algorithms, we introduce a technique for exact and automated circuit simplification. This technique is applicable to a broad class of quantum circuits and critical to enable demonstration on near-term quantum devices. We apply the proposed combination of active-space selection and circuit simplification to the dissociation of water on a magnesium surface using classical simulators and quantum hardware. Our study identifies reactions of molecules on surfaces, in conjunction with the proposed algorithmic workflow, as a promising research direction in the field of quantum computing applied to materials science.
Abstract Modeling electronic systems is an important application for quantum computers. In the context of materials science, an important open problem is the computational description of chemical reactions on surfaces. In this work, we outline a workflow to model the adsorption and reaction of molecules on surfaces using quantum computing algorithms. We develop and compare two local embedding methods for the systematic determination of active spaces. These methods are automated and based on the physics of molecule-surface interactions and yield systematically improvable active spaces. Furthermore, to reduce the quantum resources required for the simulation of the selected active spaces using quantum algorithms, we introduce a technique for exact and automated circuit simplification. This technique is applicable to a broad class of quantum circuits and critical to enable demonstration on near-term quantum devices. We apply the proposed combination of active-space selection and circuit simplification to the dissociation of water on a magnesium surface using classical simulators and quantum hardware. Our study identifies reactions of molecules on surfaces, in conjunction with the proposed algorithmic workflow, as a promising research direction in the field of quantum computing applied to materials science.
Quantum-chemical simulations are essential for predicting energies of chemical reactions. Accurately solving the many-body Schrödinger equation for reagent and product states of most relevant chemical processes is, however, unfeasible. Quantum computing offers a pathway for predicting energies of correlated electronic systems with localized interactions. Here, we apply a quantum embedding approach for investigating oxygen reduction reactions at the electrode surface in Lithium batteries, a representative example of energetic analysis in localized chemical reactions. We employ an active space selection method based on density difference analysis for identifying the orbitals involved in the reaction. Leveraging the Local Unitary Cluster Jastrow ansatz for state preparation, the active-space orbitals are then processed on a quantum computer. As quantum algorithms, we use Sample-based Quantum Diagonalization, SQD, and its extended version, Ext-SQD, which integrates electronic excitations into the quantum-selected electronic configuration subspace. The largest configurations are represented by quantum circuits mapped onto 80 qubits of an IBM Heron R2 quantum processing unit. For up to 12 orbitals, we are able to benchmark the quantum-computed reaction energies against results obtained with Complete Active Space Configuration Interaction. For benchmarking results in active spaces as large as 32 orbitals, we resort to Heat-Bath Configuration Interaction and Coupled Cluster Singles and Doubles calculations, respectively. At 27 orbitals, the Ext-SQD results exhibit prediction accuracy improvements with regard to the standard, quantum-chemical reference methods that remain computationally feasible at that scale. The results indicate the potential of sample-based quantum diagonalization for performing high-accuracy reaction modeling in chemistry and materials science.
Predictive simulations of catalytic interfaces require correlated electronic-structure treatments that describe localized chemical transformations while retaining the influence of the extended metallic environment. We introduce QC-DFET, a quantum-computing density-functional embedding framework that maps surface-reaction active spaces to compact, environment-aware qubit Hamiltonians. A reaction-consistent active-space protocol preserves orbital continuity along reaction coordinates, while quantum-selected configuration interaction based on measurements from the Zuchongzhi superconducting quantum processor and strongly contracted perturbation theory capture static and dynamic correlation. On Cu(111), QC-DFET treats active spaces up to 28 qubits and is validated through a hierarchy of experimentally constrained surface-chemistry challenges. H2 dissociation/desorption tests balanced bond breaking and recombination barriers, CO adsorption tests site selectivity and metal-adsorbate bonding, and formate hydrogenation tests competing hydrogenation branches with different kinetic and thermodynamic signatures. Across these cases, QC-DFET reproduces bidirectional H2 barriers, recovers the observed top-site preference and adsorption strength of CO, and reconciles the experimentally benchmarked H2COO* reverse barrier with the lower forward barrier to HCOOH*. These results establish embedded quantum computing as a practical route to correlated surface-reaction energetics.
methods which account for static and dynamic electron correlation, while arbitrarily accurate for small systems, are generally too expensive to be applied to modelling of catalytic cycles, due to prohibitive time and space computational complexity with respect to the size of the active space. Recent advances in quantum computing give hope for enabling access to accurate ab initio methods at scale. Herein, we present a prototype hybrid quantum-classical workflow for modeling chemical reactions on surfaces, applied to proof-of-concept models of activation and dissociation of nitrogen on small Fe clusters and a single-layer (221) iron surface. First, we determined the structures of species present in the catalytic cycle at DFT level and studied their electronic structure using CASSCF. We show that it is possible to decouple the half-filled Fe-3d band from the Fe-N and N-N bond orbitals, thereby reducing the active space significantly. Subsequently, we translated the CASSCF wavefunctions into corresponding qubit quantum states, using the Adaptive Variational Quantum Eigensolver, and estimated their energies using a state vector simulator, H1-1E quantum emulator and (for selected systems) H1-1 quantum computer. We demonstrated that if a sufficiently small active orbital space is chosen, ground state energies obtained with classical methods and with the quantum computer are in reasonable agreement. We argue that once quantum computing methods are scaled up so that larger active spaces are accessible, they can offer a tremendous practical advantage to the computational catalysis community.
Using quantum computers to tackle problems in chemistry is one of the first and most promising fields of application for near-term available quantum computers in the noisy intermediate-scale quantum computer (NISQ) era. To present applications, we start from small examples well known in classical computing and move to open questions we—all of humanity—face in solving issues like reducing the amount of energy needed to produce fertilizers and optimizing the underlying chemical reaction. We begin with the general process of simulating chemical molecules with quantum computers. Our first goal is to use potential energy surface calculation to determine the bonding energy of H2. We then adapt the example to CO2 and the metal—organic frameworks used in direct air capture. We show the adaptions and restrictions for current quantum computer hardware with a limited number of qubits during the NISQ era. We present the Haber— Bosch process and its importance for all humanity and investigate how and in which areas quantum computing can be applied in future research. We conclude the chapter by outlining what still needs to be done in research to constructively apply quantum computing in chemistry.
The predictive simulation of molecules and materials has had a broad and significant impact. It nevertheless remains constrained by the cost of accurately treating electronic correlation, excited states, and complex energy landscapes. Quantum computing offers a fundamentally different computational paradigm in which quantum states are encoded and manipulated directly rather than approximated on classical hardware. Here we discuss where this approach may provide a genuine scientific advantage in chemistry, materials science, and biochemistry. Promising directions include the high-accuracy treatment of correlated active spaces, improved excited-state simulations, and accelerated exploration of combinatorial structure spaces. The central challenge is therefore not qubit scaling alone, but demonstrably chemically meaningful gains in predictive reliability. We argue that near-term value is most likely to come from disciplined workflow integration rather than wholesale replacement of classical methods. Noisy physical devices, error-mitigated utility experiments, early fault-tolerant devices, and fully fault-tolerant quantum computers offer different scientific prospects, and claims of usefulness must be tied to the specific regime being discussed. Quantum computing will become scientifically valuable when it demonstrably reduces uncertainty in computed energies, rates, spectra, or materials stability after the full costs of state preparation, measurement, error handling, and coupling to classical simulation are included.
We present a multilevel embedding framework for quantum chemistry calculations on a quantum computer. In our framework, a quantum algorithm treats the strongly correlated active space, while a high-level wave-function method such as coupled cluster theory or multireference perturbation theory recovers the remaining correlation in the surrounding region. A sampling-based quantum algorithm, quantum-selected configuration interaction, bridges the quantum and classical treatments. The entire calculation is embedded in a low-cost density functional theory description of the surrounding environment using Manby's projection technique. We apply the framework to organic, metal-organic, and metallic systems, computing bond dissociation energies, adsorption energies, and reaction barriers using only the subset of qubits of a 144-qubit superconducting quantum computer at the University of Osaka and achieving $\sim$1 kcal/mol agreement with classical references for a Menshutkin SN2 reaction inside a carbon nanotube. Our results may open the way to quantitatively reliable quantum-classical hybrid calculations for large-scale chemical systems.