Depending on the software package and the BaTiO 3 x H y electronic state either the generalized gradient approximation exchange-correlation functional of Perdew, Burke, and Ernzerhof (PBE) , PBE with the Hubbard correction or PBE0 were used for the calculations. A detailed description of the computational methods is provided in the supplementary material.