Ab initio investigations at the MP2, CCSD(T), and MRCISD levels of theory with augmented triple-basis sets have identified and characterized various stationary points on the Be/(H 2 ) n , n ) 1-3, hypersurfaces. The van der Waals complexes, Be(H 2 ) n , are very weakly bound (D e ) 0.08-0.32 kcal/mol with respect to H 2 loss) with H 2 /H 2 interactions playing an important role in determining equilibrium structures which can be understood in terms of the various relevant long-range potentials. The covalent molecule, BeH 2 , is found to have a linear, centrosymmetric structure and to be strongly bound with respect to Be + H 2 , in agreement with previous calculations. BeH 2 interacts weakly with additional H 2 molecules (D e < 0.75 kcal/mol) which are positioned parallel to the near-linear BeH 2 moiety in the equilibrium structures of the BeH 2 (H 2 ) n-1 complexes. Of particular interest is the dramatic change in the nature of the transition state for BeH 2 production depending on the number of H 2 molecules present. For n ) 1, the reaction proceeds stepwise: first breaking the H 2 bond and forming one BeH bond followed by forming the second BeH bond. This process has an activation energy of about 56 kcal/mol. For n ) 2, the reaction proceeds via a pericyclic mechanism through a planar cyclic transition state where two H 2 bonds are broken while simultaneously two BeH bonds and one new H 2 bond are formed. The activation energy for this process decreases from the n ) 1 value to about 38 kcal/mol. For n ) 3, the reaction proceeds through a true insertion mechanism with the addition of the third H 2 molecule, decreasing the activation energy to about 33 kcal/mol. The results are discussed in comparison to the isoelectronic B + /nH 2 systems where significant σ bond activation through a cooperative interaction mechanism has been identified.