Ab initio simulations are combined with in situ infrared spectroscopy to unveil the molecular transport of H2, CO2, and H2O in the metal organic framework MOF-74-Mg. Our study uncovers-at the atomistic level-the major factors governing the transport mechanism of these small molecules. In particular, we identify four key diffusion mechanisms and calculate the corresponding diffusion barriers, which are nicely confirmed by time-resolved infrared experiments. We also answer a longstanding question about the existence of secondary adsorption sites for the guest molecules, and we show how those sites affect the macroscopic diffusion properties. Our findings are important to gain a fundamental understanding of the diffusion processes in these nano-porous materials, with direct implications for the usability of MOFs in gas sequestration and storage applications. [5,6], MOFs are also interesting for carbon-capture applications. These qualities, combined with their low production cost, have made MOFs the target of many studies, focusing mostly on their adsorption properties [2, 3,[7][8][9]. However, their performance for practical storage and capture applications also critically depends on the diffusion and interaction of gases (e.g. with water) in the MOF environment, which are currently poorly understood [6,10,11].To address this issue and elucidate the diffusion process of small molecules in the nano-pores of MOFs, we use a combination of ab initio simulations and in situ infrared (IR) spectroscopy. Specifically, we study the diffusion of H 2 , CO 2 , and H 2 O in MOF-74-Mg and investigate the role of water on diffusion of small molecules. We focus on this particular MOF, as it has attracted a lot of attention due to its enhanced reactivity with small molecules, caused by the exposure of open metal sites in its nano-pores. Britt et al. [6] have shown that MOF-74-Mg is extremely efficient in capturing CO 2 , compared to iso-structural MOFs with other metal sites such as Zn, Mn, Fe, Co, Ni, and Cu.To model the molecular diffusion in the MOF structure, we use climbing-image nudged elastic band [12] (NEB) simulations coupled with density functional theory (DFT), as implemented in QuantumEspresso [13]. To correctly capture the weak van der Waals forceswhich are critical in this study-we use the truly nonlocal functional vdW-DF [14][15][16], which has already been successfully applied to a number of related studies [17][18][19][20][21]. We use ultrasoft pseudopotentials with wavefunction and density cutoffs of 480 eV and 4800 eV, respectively. Tests show that Γ-point calculations are sufficient and yield total energies converged to within 5 meV with respect to denser k -point meshes; however, energy differences-important for our diffusion barriersare converged to within less than 1 meV. The selfconsistency tolerance was set to 1.4 × 10 −10 eV and during optimizations the total forces were relaxed to less than 2.6 × 10 −4 eV/Å. We started from the experimental MOF-74-Mg structure, optimizing the internal parameters and keep...