Investigation of surface and near-surface effects on hydrogen desorption kinetics of MgH2

Desorption of hydrogen atoms from the (110) surface of rutile magnesium hydride (MgH₂) was investigated using density functional theory (DFT) and pseudopotential method. System was represented by (110) (2×2) slab MgH₂ supercell with 12 atomic layers along the z-axis. The H-desorption was modeled by the successive release of the four two-fold bonded H atoms from the (110) surface of MgH₂. Dependence of the H-desorption energy on number and configuration of remaining surface hydrogen atoms has been determined. The features of the H atoms diffusion from the bulk towards the surface have been investigated, too. The results suggest that decrease in number of surface H atoms generally lowers the H-desorption energy in each desorption step and that both the H–H and the Mg–H interatomic interactions strongly influence the H-desorption process. The hydrogen vacancy formation energy in the first three sub-surface layers also exhibits a pronounced dependence on concentration. These findings lead to the conclusion that tendency of the MgH₂ (110) surface to preserve a maximum possible surface H concentration in its most stable configuration is the limiting factor for the H-desorption kinetics. In principle, the obtained results allow us to determine preferred paths of surface and sub-surface H-diffusion for a wide range of H concentrations and the principle features of the MgH₂ dehydrogenation process, at least for the H-rich region. Being rather comprehensive, the approach is applicable for other metal hydrides, as well.