The different roles of Pu-oxide overlayers in the hydrogenation of Pu-metal: An ab initio molecular dynamics study based on vdW-DFT+U
arXiv:1403.1087 · doi:10.1063/1.4873418
Abstract
Based on the van der Waals density functional theory (vdW-DFT)+U scheme, we carry out the ab initio molecular dynamics (AIMD) study of the interaction dynamics for H impingement against the stoichiometric PuO(111), the reduced PuO(111), and the stoichiometric -PuO(111) surfaces. The hydrogen molecular physisorption states, which can not be captured by pure DFT+\textit{U} method, are obtained by employing the vdW-DFT+\textit{U} scheme. We show that except for the weak physisorption, PuO(111) surfaces are so difficult of access that almost all of the H molecules will bounce back to the vacuum when their initial kinetic energies are not sufficient. Although the dissociative adsorption of H on PuO(111) surfaces is found to be very exothermic, the collision-induced dissociation barriers of H are calculated to be as high as eV and eV for stoichiometric and reduced PuO surfaces, respectively. Unlike PuO, our AIMD study directly reveals that the hydrogen molecules can penetrate into -PuO(111) surface and diffuse easily due to the \ native O vacancies located along the 111 diagonals of -PuO matrix. By examining the temperature effect and the internal vibrational excitations of H, we provide a detailed insight into the interaction dynamics of H in -PuO. The optimum pathways for hydrogen penetration and diffusion, the corresponding energy barriers ( eV and eV, respectively) and rate constants are systematically calculated. Overall, our study fairly reveals the different interaction mechanisms between H and Pu-oxide surfaces, which have strong implications to the interpretation of experimental observations.
29 pages, 8 figures