INTERFACE Force Field for Alumina with Validated Bulk Phases and a pH-Resolved Surface Model Database for Electrolyte and Organic Interfaces

TL;DR

This paper introduces a highly accurate, transferable force field and surface model database for alumina phases, enabling realistic simulations of electrolyte and organic interfaces across pH ranges with broad applications.

Contribution

The authors developed a unified, pH-resolved force field and surface model database for alumina, surpassing existing models in accuracy and transferability, and enabling simulations of complex interfaces.

Findings

Simulations reproduce experimental data with over 95% accuracy.

First transferable model for alumina surface ionization and charge regulation.

Predicted adsorption energies and kinetics match experimental trends.

Abstract

Alumina and aluminum oxyhydroxides underpin chemical-engineering technologies from heterogeneous catalysis, corrosion protection, functional coatings, energy-storage devices, to biomedical components. Yet molecular models that predictively connect phase structure, pH-dependent surface chemistry, electrolyte organization, and adsorption across operating conditions remain limited. Here we introduce a unified INTERFACE Force Field (IFF) parameterization together with a curated, ready-to-use pH-resolved surface model database that provides the most accurate and transferable atomistic description of major alumina phases to date. The framework covers a-Al2O3, g-Al2O3, boehmite, diaspore, and gibbsite using a single, physically interpretable parameter set that is directly compatible with CHARMM, AMBER, OPLS-AA, CVFF, and PCFF. Across structural, thermodynamic, mechanical, and interfacial benchmarks, simulations reproduce experimental reference data with more than 95 percent accuracy, exceeding existing force fields and the reliability of current density-functional approaches. A key advance is the first transferable treatment of surface ionization and charge regulation across alumina phases over a broad range of pH values, enabling simulations of realistic solid electrolyte interfaces without phase-specific reparameterization. Quantitative reliability is demonstrated by reproducing trends in zeta potentials and pH-dependent adsorption of a corrosion inhibitor at alumina-water interfaces. Predicted adsorption free energies and surface contact times correlate with experiments across more than an order of magnitude. Relative to ML-DFT workflows, the speed 100 to 1000 times faster, reaching system sizes and time scales inaccessible to quantum methods. The results establish a predictive computational platform to design alumina-containing functional materials under realistic process conditions.