General. All manipulations were performed under an inert atmosphere of argon by using Schlenk techniques or in a MBraun inert-gas glovebox. Hydrogen (99.999% purity) was purchased from Messer Austria and used as received. Solvents were purified according to standard procedures. Deuterated solvents were purchased from euriso-top and dried over 4 Å molecular sieves. All alkyne substrates were obtained from commercial sources and used as received. Complexes 1 and 2 were prepared according to literature procedures. 1,2 HRMS spectra were recorded on an Agilent 7200B GC/Q-TOF (EI) or a Bruker maXis UHR-TOF (ESI) spectrometer. 1 H, 13 C{ 1 H}, and 31 P{ 1 H} NMR spectra were recorded on Bruker AVANCE-250, AVANCE-400, DRX-500 and ASCEND-600 spectrometers. 1 H and 13 C{ 1 H} NMR spectra were referenced internally to residual protio-solvent, and solvent resonances, respectively, and are reported relative to tetramethylsilane (δ = 0 ppm). 31 P{ 1 H} NMR spectra were referenced externally to H3PO4 (85%) (δ = 0 ppm). Crystallization and X-ray structure determination of [Fe(PNP NMe-iPr)(H)(η 2-AlH4)]2 (1). Crystalls suitable for X-ray diffraction could be obtained by slow evaporation of a concentrated solution of 1 (0.04 mmol / 0.5 mL) diluted with 3.0 mL of n-pentane. X-ray diffraction data of 1 (CCDC 1951434) were collected at T = 100 K in a dry stream of nitrogen on a Bruker Kappa APEX II diffractometer system using graphite-monochromatized Mo-Kα radiation (λ = 0.71073 Å) and fine sliced φ-and ω-scans. Data were reduced to intensity values with SAINT and an absorption correction was applied with the multi-scan approach implemented in SADABS. 3 The structures were solved by the dual space method implemented in SHELXT 4 and refined against F 2 with SHELXL. 5 Non-hydrogen atoms were refined with anisotropic displacement parameters. The H atoms connected to C atoms were placed in calculated positions and thereafter refined as riding on the parent atoms. Hydride Hs were located in difference Fourier maps and refined freely. Molecular graphics were generated with the program MERCURY. 6 Computational Details. The computational results presented have been achieved in part using the Vienna Scientific Cluster (VSC). All calculations were performed using the GAUSSIAN 09 software package, 7 without symmetry constraints. The optimized geometries were obtained with the PBE0 functional. That functional uses a hybrid generalized gradient approximation (GGA), including 25 % mixture of Hartree-Fock 8 exchange with DFT 9 exchange-correlation, given by Perdew, Burke and Ernzerhof functional (PBE). 10 The basis set used for the geometry optimizations (basis b1) consisted of the Stuttgart/Dresden ECP (SDD) basis set 11 to describe the electrons of iron, and a standard 6-31G(d,p) basis set 12 for all other atoms. Transition state optimizations were performed with the Synchronous Transit-Guided Quasi-Newton Method (STQN) developed by Schlegel et al, 13 following extensive searches of the Potential Energy Surface. Frequency calculations were perfor...