The exchange energy, i.e. the splitting $\Delta E$ between gerade and ungerade states in the hydrogen molecule has proven very difficult in numerical calculation at large internuclear distances $R$, while known results are sparse and highly inaccurate. On the other hand, there are conflicting analytical results in the literature concerning its asymptotics... In this work we develop a flexible and efficient numerical approach using explicitly correlated exponential functions and demonstrate highly accurate exchange energies for internuclear distances as large as 57.5 au. This approach may find further applications in calculations of inter-atomic interactions. In particular, our results support the asymptotics form $\Delta E \sim R^{5/2}e^{-2R}$, but with the leading coefficient being $2\,\sigma$ away from the analytically derived value. (read more)