Advancing Surface Chemistry with Large-Scale Ab-Initio Quantum Many-Body Simulations

Predictive simulation of surface chemistry is of paramount importance for progress in fields from catalysis to electrochemistry and clean energy generation. Ab-initio quantum many-body methods should be offering deep insights into these systems at the electronic level, but are limited in their efficacy by their steep computational cost. In this work, we build upon state-of-the-art correlated wavefunctions to reliably converge to the `gold standard' accuracy in quantum chemistry for application to extended surface chemistry. Efficiently harnessing graphics processing unit acceleration along with systematically improvable multiscale resolution techniques, we achieve linear computational scaling up to 392 atoms in size. These large-scale simulations demonstrate the importance of converging to these extended system sizes, achieving a validating handshake between simulations with different boundary conditions for the interaction of water on a graphene surface. We provide a new benchmark for this water-graphene interaction that clarifies the preference for water orientations at the graphene interface. This is extended to the adsorption of carbonaceous molecules on chemically complex surfaces, including metal oxides and metal-organic frameworks, where we consistently achieve chemical accuracy compared to experimental references, and well inside the scatter of traditional density functional material modeling approaches. This pushes the state of the art for simulation of molecular adsorption on surfaces, and marks progress into a post-density functional era for more reliable and improvable approaches to first-principles modeling of surface problems at an unprecedented scale and accuracy using ab-initio quantum many-body methods.