Fitting an ab initio HF–HF potential surface