Very accurate interaction energies and equilibrium geometries are obtained for the C2v Cd(1S, 3P, 1P)–H2 van der Waals (VDW) complexes using extensive ab initio multireference configuration interaction plus second order multireference Möller–Plesset calculations. Gaussian basis sets are of triple-ζ plus polarization quality and a relativistic effective core potential for cadmium has been used. For the interaction energies the Boys and Bernardi counterpoise (CP) correction has been used. The ground state Cd(1S)–H2 van der Waals complex was found to be very weakly bound with a 108 cm−1 well at 8.4 a.u., as compared with the estimated experimental values of De ≂50 cm−1 at 8.3 a.u. The computed (De=1580 cm−1, Re=4.9 a.u.) constants for the 1B1 Cd(1 P)–H2 exciplex are in good agreement with the recent experimentally derived values (De=1713 cm−1, Re=4.7 a.u.) obtained by Wallace et al. [J. Chem. Phys. 97, 3135 (1992)]. The well depth and equilibrium geometry for the 3B1 Cd (3P) –H2 exciplex were computed as 420 cm−1 and 5.7 a.u. while those of the 3B2 exciplex are 650 cm−1 and 4.75 a.u. The 1B2 surface did not lead to a VDW complex but instead it produced a strongly bound (De≂10 kcal/mol) bent H–Cd–H molecule responsible for the experimentally observed CdH+H and Cd+H+H reaction channels. The 3,1A1 states arising from the Cd(3,1P)+H2 asymptotes were found to be totally repulsive. In all cases the CP correction was a non-negligible fraction of the interaction energy. Finally, the energetic position and the geometry of the 1B1/3A1 surface crossing, crucial to explain the observed Cd(1B1,v′)–H2 predissociation lifetimes, are accurately determined.