To enable rational design of alloy nanoparticle catalysts,
we develop
an approach to generate catalytic activity maps of alloy nanoparticles
on a grid of particle size and composition. The catalytic activity
maps are created by using a quaternary cluster expansion to explicitly
predict adsorbate binding energies on alloy nanoparticles of varying
shape, size, and atomic order while accounting for interactions among
the adsorbates. This cluster expansion is used in kinetic Monte Carlo
simulations to predict activated nanoparticle structures and turnover
frequencies on all surface sites. We demonstrate our approach on Pt–Ni
octahedral nanoparticle catalysts for the oxygen reduction reaction
(ORR), revealing that the specific activity is predicted to be optimized
at an edge length of larger than 5.5 nm and a composition of about
Pt0.85Ni0.15 and the mass activity is predicted
to be optimized at an edge length of 3.3–3.8 nm and a composition
of about Pt0.8Ni0.2