Theory of divalent main group H2 activation: electronics and quasiclassical trajectories

Density functional theory (DFT), absolutely localized molecular orbital (ALMO) analysis, and quasiclassical trajectories (QCTs) were used to study the structure, barrier heights, thermodynamics, electronic properties, and dynamics of dihydrogen (H2) activation by singlet divalent main group compounds (ER2; E = C, Si, Ge). ALMO energy and charge decomposition calculations reveal that in the transition state CR2 acts as an ambiphile toward H2 because of equal forward-bonding and back-bonding orbital stabilization while SiR2 and GeR2 act as nucleophiles with dominant orbital energy stabilization arising from ER2 to H2 donation. Frontier molecular orbital (FMO) energy gaps do not provide a reasonable estimate of energy stabilization gained between the ER2 and H2 in the transition state or an accurate description of the nucleophilic versus electrophilic character because of electron repulsion and orbital overlap influences that are neglected. In CR2 transition states, forward-bonding and back-bonding are maximized in the nonleast motion geometry. In contrast, SiR2/GeR2 transition states have side-on geometries to avoid electron-electron repulsion. Electron repulsion, rather than orbital interactions, also determines the relative barrier heights of CR2 versus SiR2/GeR2 reactions. Examination of barrier heights and reaction energies shows a clear kinetic-thermodynamic relationship for ER2 activation of H2. A computational survey of R groups on ER2 divalent atom centers was performed to explore the possibility for H2 activation to occur with a low barrier and thermodynamically reversible. QCTs show that dihydrogen approach and reaction with CR2 may involve geometries significantly different than the static transition-state structure. In contrast, trajectories for dihydrogen addition to SiR2 involve geometries close to the side-on approach suggested by the static transition-state structure. QCTs also demonstrate that addition of H2 to CR2 and SiR2 is dynamically concerted with the average time gap of bond formation between E-H bonds of approximately 11 and 21 fs, respectively.