Bayesian Optimization of High-Entropy Alloy Compositions for Electrocatalytic Oxygen Reduction*.

Active, selective and stable catalysts are imperative for sustainable energy conversion, and engineering materials with such properties are highly desired. High-entropy alloys (HEAs) offer a vast compositional space for tuning such properties. Too vast, however, to traverse without the proper tools. Here, we report the use of Bayesian optimization on a model based on density functional theory (DFT) to predict the most active compositions for the electrochemical oxygen reduction reaction (ORR) with the least possible number of sampled compositions for the two HEAs Ag-Ir-Pd-Pt-Ru and Ir-Pd-Pt-Rh-Ru. The discovered optima are then scrutinized with DFT and subjected to experimental validation where optimal catalytic activities are verified for Ag-Pd, Ir-Pt, and Pd-Ru binary alloys. This study offers insight into the number of experiments needed for optimizing the vast compositional space of multimetallic alloys which has been determined to be on the order of 50 for ORR on these HEAs.

What Atomic Positions Determines Reactivity of a Surface? Long-Range, Directional Ligand Effects in Metallic Alloys.

Ligand and strain effects can tune the adsorption energy of key reaction intermediates on a catalyst surface to speed up rate-limiting steps of the reaction. As novel fields like high-entropy alloys emerge, understanding these effects on the atomic structure level is paramount: What atoms near the binding site determine the reactivity of the alloy surface? By statistical analysis of 2000 density functional theory calculations and subsequent host/guest calculations, it is shown that three atomic positions in the third layer of an fcc(111) metallic structure fourth-nearest to the adsorption site display significantly increased influence on reactivity over any second or third nearest atomic positions. Subsequently observed in multiple facets and host metals, the effect cannot be explained simply through the d-band model or a valence configuration model but rather by favorable directions of interaction determined by lattice geometry and the valence difference between host and guest elements. These results advance the general understanding of how the electronic interaction of different elements affect adsorbate-surface interactions and will contribute to design principles for rational catalyst discovery of better, more stable and energy efficient catalysts to be employed in energy conversion, fuel cell technologies, and industrial processes.

Complex-Solid-Solution Electrocatalyst Discovery by Computational Prediction and High-Throughput Experimentation*.

Complex solid solutions ("high entropy alloys"), comprising five or more principal elements, promise a paradigm change in electrocatalysis due to the availability of millions of different active sites with unique arrangements of multiple elements directly neighbouring a binding site. Thus, strong electronic and geometric effects are induced, which are known as effective tools to tune activity. With the example of the oxygen reduction reaction, we show that by utilising a data-driven discovery cycle, the multidimensionality challenge raised by this catalyst class can be mastered. Iteratively refined computational models predict activity trends around which continuous composition-spread thin-film libraries are synthesised. High-throughput characterisation datasets are then used as input for refinement of the model. The refined model correctly predicts activity maxima of the exemplary model system Ag-Ir-Pd-Pt-Ru. The method can identify optimal complex-solid-solution materials for electrocatalytic reactions in an unprecedented manner.