A bond-order LiFH potential energy surface for 3D quantum-mechanical calculations