Adsorption and evolution of hydrogen molecules on hexagonal boron nitride monolayer: a combined DFT and kinetic monte-carlo simulations study

Density functional theory (DFT) and kinetic Monte-Carlo (kMC) simulation code has been combinedly used to study the adsorption and evolution dynamics of hydrogen molecules over a hexagonal boron nitride (h-BN) monolayer. Maximum adsorption energy from van der Waals curve is predicted to be around 60 to 70 meV using two different DFT functionals. Repulsive lateral interaction between two hydrogen molecules plays a key role in determining the maximum number of adsorptions inside one unit cell of h-BN. Bader charge analysis, electron localization function (ELF), total and partial density of states (DOS) plots have been included to understand the weak interaction going on between the adsorbent and substrate. The input energy parameters from the DFT calculation has been used to perform the kMC simulation for describing the adsorption, desorption and the diffusion pattern of hydrogen molecules with a given time of exposure to an empty h-BN substrate along with the overall surface coverage.