Development of a Multicenter Density Functional Tight Binding Model for Plutonium Surface Hydriding

We detail the creation of a multicenter density functional tight binding (DFTB) model for hydrogen on 3 plutonium, using a framework of new Slater-Koster interaction parameters and a repulsive energy based on the Chebyshev Interaction Model for Efficient Simulation (ChIMES), where two- and three-center atomic interactions are represented by linear combinations of Chebyshev polynomials. We find that our DFTB/ChIMES model yields a total electron density of states for bulk delta-Pu that compares well to that from Density Functional Theory, as well as to a grid of energy calculations representing approximate H-2 dissociation paths on the delta-Pu (100) surface. We then perform molecular dynamics simulations and minimum energy pathway calculations to determine the energetics of surface dissociation and subsurface diffusion on the (100) and (111) surfaces. Our approach allows for the efficient creation of multicenter repulsive energies with a relatively small investment in initial DFT calculations. Our efforts are particularly pertinent to studies that rely on quantum calculations for interpretation and validation, such as experimental determination of chemical reactivity both on surfaces and in condensed phases.