Based on first principles calculations, we theoretically predict the new two-dimensional (2D) MgH 2 . The thermodynamic stability, partial density of states, electron localization function, and Bader charge of pure and the transition metal (Ti, V, and Mn) doped 2D MgH 2 are investigated. The results show that all the systems are dynamically stable, and the dehydrogenation properties indicate that the decomposition temperature can be reduced by introducing the transition metal, and the Mn doped system exhibits good performance for better hydrogen storage and dehydrogenation kinetics. 482 2 of 10 phonon spectra and heat of formation. The calculated heat of formation for pure and Ti/V/Mn doped 2D MgH 2 are −37.57, −25.67, −18.14, and −23.90 kJ/mol·H 2 , respectively, which are significantly lower than that of −75.99 kJ/mol·H 2 of bulk MgH 2 . The electronic structure and hydrogen desorption kinetics results show that the predicted two-dimensional magnesium hydride are promising candidates for hydrogen storage.
Computational DetailsThe structural optimization and electronic property calculations were performed using the projector augmented plane-wave method (PAW) based on the density functional theory (DFT) in the Vienna ab initio simulation package (VASP) [18,19]. The exchange-correlation potential was approximated by generalized gradient approximation (GGA) in the Perdew-Burke-Ernzerhof (PBE) form [20,21]. To avoid the interlayer effects of the c-axis, the vacuum region around 15 Å was set in all the systems. The energy cutoff of 600 eV and the 9 × 9 × 1 Γ-centered Monkhorst-Pack k-points [22] were employed for all calculations. The atomic positions were fully relaxed and the force tolerance between each atom was less than 0.01 eV/Å for the structural optimization. The convergence criteria of 10 −6 eV per atom was applied to be self-consistent. Meanwhile, for calculation of electronic structures, we also applied the local density approximation (LDA) [23] and HSE06 [24] was functional. The kinetic stability was discussed using the phonon spectra calculations in PHONOPY code coupled with VASP using the density functional perturbation theory (DFPT) method [25][26][27].
Results and Discussion