Symbolic regression (SR) is an emerging method for building analytical formulas to find models that best fit data sets. Here, SR was used to guide the design of new oxide perovskite catalysts with improved oxygen evolution reaction (OER) activities. An unprecedentedly simple descriptor, {mu}/t, where {mu} and t are the octahedral and tolerance factors, respectively, was identified, which accelerated the discovery of a series of new oxide perovskite catalysts with improved OER activity. We successfully synthesized five new oxide perovskites and characterized their OER activities. Remarkably, four of them, Cs0.4La0.6Mn0.25Co0.75O3, Cs0.3La0.7NiO3, SrNi0.75Co0.25O3, and Sr0.25Ba0.75NiO3, outperform the current state-of-the-art oxide perovskite catalyst, Ba0.5Sr0.5Co0.8Fe0.2O3 (BSCF). Our results demonstrate the potential of SR for accelerating data-driven design and discovery of new materials with improved properties.
In order to estimate the reactivity of a large number of potentially complex heterogeneous catalysts while searching for novel and more efficient materials, physical as well as data-centric models have been developed for a faster evaluation of adsorption energies compared to first-principles calculations. However, global models designed to describe as many materials as possible might overlook the very few compounds that have the appropriate adsorption properties to be suitable for a given catalytic process. Here, the subgroup-discovery (SGD) local artificial-intelligence approach is used to identify the key descriptive parameters and constrains on their values, the so-called SG rules, which particularly describe transition-metal surfaces with outstanding adsorption properties for the oxygen reduction and evolution reactions. We start from a data set of 95 oxygen adsorption energy values evaluated by density-functional-theory calculations for several monometallic surfaces along with 16 atomic, bulk and surface properties as candidate descriptive parameters. From this data set, SGD identifies constraints on the most relevant parameters describing materials and adsorption sites that (i) result in O adsorption energies within the Sabatier-optimal range required for the oxygen reduction reaction and (ii) present the largest deviations from the linear scaling relations between O and OH adsorption energies, which limit the performance in the oxygen evolution reaction. The SG rules not only reflect the local underlying physicochemical phenomena that result in the desired adsorption properties but also guide the challenging design of alloy catalysts.
NiFe oxyhydroxide is one of the most promising oxygen evolution reaction (OER) catalysts for renewable hydrogen production, and deciphering the identity and reactivity of the oxygen intermediates on its surface is a key challenge but is critical to understanding the OER mechanism as well as designing water-splitting catalysts with higher efficiencies. Here, we screened and utilized in situ reactive probes that can selectively target specific oxygen intermediates with high rates to investigate the OER intermediates and pathway on NiFe oxyhydroxide. Most importantly, the oxygen atom transfer (OAT) probes (e.g. 4-(Diphenylphosphino) benzoic acid) could efficiently inhibit the OER kinetics by scavenging the OER intermediates, exhibiting lower OER currents, larger Tafel slopes and larger kinetic isotope effect values, while probes with other reactivities demonstrated much smaller effects. Combining the OAT reactivity with electrochemical kinetic and operando Raman spectroscopic techniques, we identified a resting Fe=O intermediate in the Ni-O scaffold and a rate-limiting O-O chemical coupling step between a Fe=O moiety and a vicinal bridging O. DFT calculation further revealed a longer Fe=O bond formed on the surface and a large kinetic energy barrier of the O-O chemical step, corroborating the experimental results. These results point to a new direction of liberating lattice O and expediting O-O coupling for optimizing NiFe-based OER electrocatalyst.
Considering the recent breakthroughs in the synthesis of novel two-dimensional (2D) materials from layered bulk structures, ternary layered transition metal borides, known as MAB phases, have come under scrutiny as a means of obtaining novel 2D transition metal borides, so-called MBene. Here, based on a set of phonon calculations, we show the dynamic stability of many Al-containing MAB phases, MAlB (M = Ti, Hf, V, Nb, Ta, Cr, Mo, W, Mn, Tc), M$_2$AlB$_2$ (Sc, Ti, Zr, Hf, V, Cr, Mo, W, Mn, Tc, Fe, Rh, Ni), M$_3$Al$_2$B$_2$ (M = Sc, T, Zr, Hf, Cr, Mn, Tc, Fe, Ru, Ni), M$_3$AlB$_4$ (M = Sc, Ti, Zr, Hf, V, Nb, Ta, Cr, Mo, W, Mn, Fe), and M$_4$AlB$_6$ (M = Sc, Ti, Zr, Hf, V, Nb, Ta, Cr, Mo). By comparing the formation energies of these MAB phases with those of their available competing binary M$-$B and M$-$Al, and ternary M$-$Al$-$B phases, we find that some of the Sc-, Ti-, V-, Cr-, Mo-, W-, Mn-, Tc-, and Fe-based MAB phases could be favorably synthesized in an appropriate experimental condition. In addition, by examining the strengths of various bonds in MAB phases via crystal orbital Hamilton population and spring constant calculations, we find that the B$-$B and then M$-$B bonds are stiffer than the M$-$Al and Al$-$B bonds. The different strength between these bonds implies the etching possibility of Al atoms from MAB phases, consequently forming various 2D MB, M$_2$B$_3$, and M$_3$B$_4$ MBenes. Furthermore, we employ the nudged elastic band method to investigate the possibility of the structural phase transformation of the 2D MB MBenes into graphene-like boron sheets sandwiched between transition metals and find that the energy barrier of the transformation is less than $0.4$ eV/atom.