Quantum Calculations of Hydrogen Absorption and Diffusivity in Bulk CeO₂

CeO₂ (ceria) is an attractive material for heterogeneous catalysis applications involving hydrogen due to its favorable redox activity combined with its relative impermeability to hydrogen ions and molecules. However, to date, many bulk ceria/hydrogen properties remain unresolved in part due to a scarcity of experimental data combined with quantum calculation results that vary according to the approach used. In this regard, we have conducted a series of density functional theory (DFT) calculations utilizing generalized gradient (GGA), metaGGA, and hybrid functionals as well as several corrections for electronic correlations, applied to a number of properties regarding hydrogen in bulk stoichiometric CeO₂. Our calculations place reasonable bounds on the lattice constants, band gaps, hydrogen absorption energies, and O-H bond vibrational frequencies that can be determined by DFT. In addition, our results indicate that the activation energy barriers for hydrogen bulk diffusion are uniformly low (<0.15 eV) for the calculation parameters probed here and that, in general, the effect of hydrogen tunneling is small at ambient temperatures. Our study provides a recipe to determine fundamental physical chemical properties of Ce-O-H interactions while also determining realistic ranges for diffusion kinetics. This can facilitate the determination of future coarse-grained models that will be able to guide and elucidate experimental efforts in this area.