Issue 8, 2018

Bimodal hole transport in bulk BiVO4 from computation


We report first principles and mesoscale kinetic simulations of intrinsic electron and hole transport in stoichiometric monoclinic BiVO4 (BVO). A main finding is that hole transport is bimodal: there are fast hole hops within VO4 tetrahedra that are not ‘transport-efficient’ and slower hole hops across VO4 tetrahedra that are ‘transport-efficient’. This bimodality in hole transport may explain the slow transport/high recombination characteristics of BVO. In this work, polaron hops were described via the two-state Marcus/Holstein model of VV–VIV and OI–OII interchanges for electron and hole polarons respectively. The relevant parameters in the theory were obtained using the density functional theory DFT+U method. Hopping activation energies were combined with Einstein diffusion theory to yield mobility as well as with kinetic Monte Carlo (KMC) modeling. All unique nearest neighbor electron and hole hopping processes were characterized. For V-to-V hops of electron polarons, we obtained activation energies of ∼0.37 eV that compare well with the experimental activation energy of ∼0.30 eV for drift mobility. Electron polarons remain localized on single sites along the hopping pathways before hopping via tunneling, suggestive of small non-adiabatic coupling and a diabatic character consistent with V-to-V electron transfer through space (no bridge linkers). These simulations yield a small electron transport anisotropy, with transport in the (a,b) plane being faster than along the c direction by a factor of ∼2. For hole polarons our calculations support the character of small polarons localized on O atoms as well. Some nearest neighbor hole hops have a diabatic character, while other nearest neighbor hops are adiabatic and phonon-assisted with an Arrhenius thermal dependence. Hops through O–V–O bridges and through space have low activation barriers (∼0.17 eV and ∼0.25 eV) while hops through O–Bi–O bridges have higher activation barriers (∼0.37 eV and higher). The mobility is gated by the slower transport across VO4 tetrahedra through O–Bi–O bridges owing to the underlying BVO lattice network. Hole mobility is determined to be about one order of magnitude larger than electron mobility consistent with THz spectroscopy measurements. This work sets the foundation to address facet selectivity of charge carriers in shape and size-tailored crystals.

Graphical abstract: Bimodal hole transport in bulk BiVO4 from computation

Supplementary files

Article information

Article type
19 Nov 2017
23 Jan 2018
First published
23 Jan 2018

J. Mater. Chem. A, 2018,6, 3714-3723

Bimodal hole transport in bulk BiVO4 from computation

T. Liu, V. Pasumarthi, C. LaPorte, Z. Feng, Q. Li, J. Yang, C. Li and M. Dupuis, J. Mater. Chem. A, 2018, 6, 3714 DOI: 10.1039/C7TA10170H

To request permission to reproduce material from this article, please go to the Copyright Clearance Center request page.

If you are an author contributing to an RSC publication, you do not need to request permission provided correct acknowledgement is given.

If you are the author of this article, you do not need to request permission to reproduce figures and diagrams provided correct acknowledgement is given. If you want to reproduce the whole article in a third-party publication (excluding your thesis/dissertation for which permission is not required) please go to the Copyright Clearance Center request page.

Read more about how to correctly acknowledge RSC content.

Social activity