Molybdenum defect complexes in bismuth vanadate.

Monoclinic bismuth vanadate (BiVO\textsubscript{4}) is a promising \textit{n}-type semiconductor for applications in sunlight-driven water splitting. Several studies have shown that its photocatalytic activity is greatly enhanced by high concentrations of Mo and W dopants. In the present work, we performed \textit{ab initio} calculations to assess the most favorable relative position between Mo-related pairs in BiVO\textsubscript{4}. Surprisingly, we verify that the lowest energy configuration for Mo\textsubscript{V} pairwise defects in BiVO\textsubscript{4} occurs on nearest-neighbor sites, despite the higher electrostatic repulsion and larger strain on the crystal lattice. Similar results were observed for W\textsubscript{V} defect pairs in W-doped BiVO\textsubscript{4}. We show that the origin of this effect lies in a favorable hybridization between the atomic orbitals of the impurities that is only verified when they are closest to each other, resulting in an enthalpy gain that overcomes the repulsive components of the pair formation energy. As a consequence, Mo and/or W doped BiVO\textsubscript{4} are likely to present donor-donor defect complexes, which is an outcome that can be applied in experimental approaches for improving the photocatalytic activity of these metal oxides.