The exchange energy, i.e. the splitting Δ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 ΔE∼R5/2e−2R, but with
the leading coefficient being 2σ away from the analytically derived
value