An analytical potential energy surface is developed from high quality ab initio calculations for the electrostatic region of the Li−+H2 interaction. The Li−(H2) electrostatic complex is found to have a linear minimum energy structure with a De of 64.44 cm−1. A numerical determination of the bound levels supported by this potential indicates a D0 of only about 7 cm−1 for Li−(para H2) and a considerably larger D0 of about 22 cm−1 for Li−(ortho H2). Altogether, the Li−(para H2) interaction is predicted to support 11 bound levels: ν3=0, J=0–6; and ν3=1, J=0–3, whereas the Li−(ortho H2) interaction is predicted to support 28 bound levels: ν3=0, J=0–10; ν3=1, J=0–8; ν3=2, J=0–5; and ν2=1, J=1–2. Analogous results for the D2 and HD isotopolouges are reported.