Alloying transition metals, such as Mo, into BiVO4 has emerged as the primary mechanism for improving carrier transport in this photoanode for solar fuels production. The present work establishes the generality of improving photoelectrochemical performance through co-alloying with a transition metal electron donor and a structure-modulating rare earth. Further improvement for all such alloys is obtained by annealing the oxide materials in H2, ultimately producing photoanodes with above 3 mA cm-2 photocurrent density under AM 1.5G illumination, in the top tier of compact BiVO4 films.

Complex oxide β-Mn2V2O7 is identified as exhibiting near-optimal band energetics for solar fuel applications among known metal oxides. Experiments, corroborated by theory, indicate a bandgap near 1.8 eV. The calculations predict that β-Mn2V2O7 has well-aligned band edge energies for the hydrogen evolution reaction and oxygen evolution reaction. Photoelectrochemical measurements indicate appreciable photocurrent, corroborating the predictions.

Improving the efficiency of solar-powered oxygen evolution is both critical for development of solar fuels technologies and challenging due to the broad set of properties required of a solar fuels photoanode. Bismuth vanadate, in particular the monoclinic clinobisvanite phase, has received substantial attention and has exhibited the highest radiative efficiency among metal oxides with a band gap in the visible range. Efforts to further improve its photoelectrochemical performance have included alloying one or more metals onto the Bi and/or V sites, with progress on this frontier stymied by the difficulty in computational modelling of substitutional alloys and the high dimensionality of co-alloying composition spaces. Since substitutional alloying concurrently changes multiple materials properties, understanding the underlying cause for performance improvements is also challenging, motivating our application of combinatorial materials science techniques to map photoelectrochemical performance of 948 unique bismuth vanadate alloy compositions comprising 0 to 8% alloys of P, Ca, Mo, Eu, Gd, and W along with a variety of compositions from each pairwise combination of these elements. Upon identification of substantial improvements in the (Mo,Gd) co-alloying space, structural mapping was performed to reveal a remarkable correlation between performance enhancement and a lowered monoclinic distortion. First-principles density functional theory calculations indicate that the improvements are due to a lowered hole effective mass and hole polaron formation energy, and collectively, our results identify the monoclinic distortion as a critical parameter in the optimization and understanding of bismuth vanadate-based photoanodes.