Total and differential cross sections for mutual neutralization in He + and H − collisions at low to intermediate (0.001 eV to 100 eV) are calculated ab initio and fully quantum mechanically. Atomic final state distributions and isotope effects are investigated. The theoretical model includes dynamics on eleven coupled states of 2 Σ + symmetry, where autoionization is incorporated. The potential energy curves, autoionization widths and non-adiabatic couplings of electronic resonant states of HeH are computed by combining structure calculations with electron scattering calculations. The nuclear dynamics is studied using a strict diabatic representation of the resonant states. Effects of rotational couplings between 2 Σ + and 2 Π electronic states are investigated in the pure precession approximation.