  We develop meanfield approximation and numerical quadrature schemes for the
evaluation of Angular-Dependent interatomic Potentials (ADPs) for magnesium and
magnesium hydrides at finite temperature (thermalization) and arbitrary atomic
molar fractions (mixing) within a non-equilibrium statistical mechanical
framework and derive local equilibrium relations. We numerically verify and
experimentally validate the accuracy and fidelity of the resulting
thermalized/mixed ADPs (TADPs) by means of selected numerical tests including
free entropy, heat capacity, thermal expansion, molar volumes, equation of
state and elastic constants. We show that the local equilibrium properties
predicted by TADPs agree closely with those computed directly from ADP by means
of Molecular Dynamics (MD).
