FIRE‐GNN: Force‐Informed, Relaxed Equivariance Graph Neural Network for Rapid and Accurate Prediction of Surface Properties

The work function and cleavage energy of a surface are critical properties that determine the viability of materials in electronic emission applications, semiconductor devices, and heterogeneous catalysis. While first‐principles calculations are accurate in predicting these properties, their computational expense, combined with the vast search space of surfaces, makes a comprehensive screening approach with density functional theory (DFT) infeasible. Here, we introduce FIRE‐GNN (Force‐Informed, Relaxed Equivariance Graph Neural Network), which integrates surface‐normal symmetry breaking and machine learning interatomic potential (MLIP)‐derived force information, achieving a twofold reduction in mean absolute error (down to 0.065 eV) over the previous state‐of‐the‐art for work function prediction. We additionally benchmark recent invariant and equivariant architectures, analyze the impact of symmetry breaking, and evaluate out‐of‐distribution generalization, demonstrating that FIRE‐GNN consistently outperforms competing models for work function predictions. This model enables accurate and rapid predictions of the work function and cleavage energy across a vast chemical space and facilitates the discovery of materials with tuned surface properties.