Theoretical studies
on molecule–metal surface reactions
have so far been limited to small surface unit cells due to computational
costs. Here, for the first time molecular dynamics simulations on
very large surface unit cells at the level of density functional theory
are performed, allowing a direct comparison to experiments performed
on a curved crystal. Specifically, the reaction of D
2
on
a curved Pt crystal is investigated with a neural network potential
(NNP). The developed NNP is also accurate for surface unit cells considerably
larger than those that have been included in the training data, allowing
dynamical simulations on very large surface unit cells that otherwise
would have been intractable. Important and complex aspects of the
reaction mechanism are discovered such as diffusion and a shadow effect
of the step. Furthermore, conclusions from simulations on smaller
surface unit cells cannot always be transfered to larger surface unit
cells, limiting the applicability of theoretical studies of smaller
surface unit cells to heterogeneous catalysts with small defect densities.