A Reproducible Workflow for Extracting Quantum Hamiltonians from Surface-Adsorbate Models
Manisha Malhotra, Neil Ramarapu, F. Rinaldi
Abstract
Abstract Bridging realistic surface-catalysis models from density functional theory (DFT) and the many-body Hamiltonian inputs required by quantum algorithms remains a practical obstacle: surface DFT yields total energies and densities, whereas algorithms such as the variational quantum eigensolver require an explicit fermionic or qubit operator. We present an open-source, reproducible workflow that takes a DFT-optimised slab, extracts a hydrogen-capped cluster around the active site, selects a frontier active space automatically, and constructs the corresponding fermionic and qubit Hamiltonians. Heavy-element cores are treated with an effective core potential (ECP) at the Hartree–Fock stage, an established molecular-cluster strategy for quantum-computing catalysis workflows. WedemonstratethepipelineonH*/OH* adsorption on rutile IrO2(110) at two adsorption sites, producing 14-qubit Hamiltonians from (8e, 7o) active spaces. The pipeline is validated end to end: the Hartree–Fock reference reconstructed from the qubit Hamiltonian reproduces the self-consistent-field energy to better than 10 −6 Hartree, and exact diagonalisation yields physically sensible static correlation energies at the milli-Hartree scale. The entire active-space pipeline runs on commodity hardware without high-performance computing. We additionally document a methodological caution: the def2-SVP basis is defined for Iridium only together with its ECP, so an all-electron treatment in this basis is ill-defined rather than merely approximate.