A Reproducible Workflow for Extracting Quantum Hamiltonians from Surface-Adsorbate Models
Manisha Malhotra, F. Rinaldi, Neil Ramarapu
Abstract
Abstract Bridging the gap between realistic surface catalysis models from density functional theory (DFT) and the input requirements of quantum algorithms is a non-trivial task. Surface DFT calculations typically model extended slabs and yield total energies or reaction energetics, whereas quantum algorithms (such as variational quantum eigensolvers) require a many-body Hamiltonian (e.g. a list of fermionic or qubit operators) as input. We present an open-source, reproducible software workflow that takes an optimized DFT slab model of a catalyst surface and produces a few-body Hamiltonian suitable for quantum simulation. Starting from a DFT-optimized slab geometry, the workflow extracts a localized cluster around the active site, identifies an appropriate set of active orbitals, and constructs the corresponding fermionic and qubit Hamiltonians. Orbital selection is guided by manual cross-referencing of semiempirical orbital contributions across neutral and anionic cluster states, using localisation and energetic proximity to the HOMO as selection criteria. This end-to-end pipeline enables researchers to bridge surface science and quantum computing in a transparent and automated manner. We demonstrate the workflow on a representative electrochemical system: a hydrogen-based adsorbate (H* or OH*) on an IrO$_2$ catalyst surface. For this system, our software produces a 14-qubit Hamiltonian (derived from a 7-orbital active space) that can be directly used in quantum algorithms such as the Variational Quantum Eigensolver (VQE \cite{VQE}). This work provides a reproducible, end-to-end pipeline bridging periodic surface DFT models and quantum-ready Hamiltonians for heterogeneous catalytic systems, enabling direct integration of realistic surface chemistry with near-term quantum algorithms.