Structural properties and enthalpy of formation of magnesium hydride from quantum Monte Carlo calculations

We have used diffusion Monte Carlo (DMC) calculations to study the structural properties of magnesium hydride (MgH₂), including the pressure-volume equation of state, the cohesive energy, and the enthalpy of formation from magnesium bulk and hydrogen gas. The calculations employ pseudopotentials and B-spline basis sets to expand the single particle orbitals used to construct the trial wave functions. Extensive tests on system size, time step, and other sources of errors, performed on periodically repeated systems of up to 1050 atoms, show that all these errors together can be reduced to below 100.3em0exmeV∕f.u.. We find excellent agreement with the experiments for the equilibrium volume of both the Mg and the MgH₂ crystals. The cohesive energy of the Mg crystal is found to be 1.51(1)0.3em0exeV and agrees perfectly with the experimental value of 1.510.3em0exeV. The enthalpy of formation of MgH₂ from Mg bulk and H₂ gas is found to be 0.85±0.010.3em0exeV∕f.u., or 82±10.3em0exkJ∕mole, which is off the experimental one of 76.1±10.3em0exkJ∕mole only by 60.3em0exkJ∕mole. This shows that DMC can almost achieve chemical accuracy (10.3em0exkcal∕mole) on this system. Density functional theory errors are shown to be much larger and depend strongly on the functional employed.