Ab initio study on the hydrogen desorption from $\rm {MH\text{–}NH}_3$MH–NH3 (M = Li, Na, K) hydrogen storage systems

The hydrogen storage system LiH + \documentclass[12pt]{minimal}\begin{document}$\rm {NH}_3$\end{document} NH 3 ↔ \documentclass[12pt]{minimal}\begin{document}$\rm {LiNH}_2$\end{document} LiNH 2 + \documentclass[12pt]{minimal}\begin{document}$\rm {H}_2$\end{document}H2 is one of the most promising hydrogen storage systems, where the reaction yield can be increased by replacing Li in LiH with other alkali metals (Na or K) in order of Li < Na < K. In this paper, we have studied the alkali metal M (M = Li, Na, K) dependence of the reactivity of MH with \documentclass[12pt]{minimal}\begin{document}$\rm {NH}_3$\end{document} NH 3 by calculating the potential barrier of the \documentclass[12pt]{minimal}\begin{document}$\rm {H}_2$\end{document}H2 desorption process from the reaction of an \documentclass[12pt]{minimal}\begin{document}$\rm {M}_2$\end{document}M2\documentclass[12pt]{minimal}\begin{document}$\rm {H}_2$\end{document}H2 cluster with an \documentclass[12pt]{minimal}\begin{document}$\rm {NH}_3$\end{document} NH 3 molecule based on the ab initio structure optimization method. We have shown that the height of the potential barrier becomes lower in order of Li, Na, and K, where the difference of the potential barrier in Li and Na is relatively smaller than that in Na and K, and this tendency is consistent with the recent experimental results. We have also shown that the H–H distance of the \documentclass[12pt]{minimal}\begin{document}$\rm {H}_2$\end{document}H2 dimer at the transition state takes larger distance and the change of the potential energy around the transition state becomes softer in order of Li, Na, and K. There are almost no M dependence in the charge of the H atom in \documentclass[12pt]{minimal}\begin{document}$\rm {NH}_3$\end{document} NH 3 before the reaction, while that of the H atom in \documentclass[12pt]{minimal}\begin{document}$\rm {M}_2$\end{document}M2\documentclass[12pt]{minimal}\begin{document}$\rm {H}_2$\end{document}H2 takes larger negative value in order of Li, Na, and K. We have also performed molecular dynamics simulations on the \documentclass[12pt]{minimal}\begin{document}$\rm {M}_2$\end{document}M2\documentclass[12pt]{minimal}\begin{document}$\rm {H}_2$\end{document}H2–\documentclass[12pt]{minimal}\begin{document}$\rm {NH}_3$\end{document} NH 3 system and succeeded to reproduce the \documentclass[12pt]{minimal}\begin{document}$\rm {H}_2$\end{document}H2 desorption from the reaction of \documentclass[12pt]{minimal}\begin{document}$\rm {Na}_2$\end{document} Na 2\documentclass[12pt]{minimal}\begin{document}$\rm {H}_2$\end{document}H2 with \documentclass[12pt]{minimal}\begin{document}$\rm {NH}_3$\end{document} NH 3.