Abstract
We combined both density functional theory and classical molecular dynamics simulations to investigate the molecular mechanisms governing hydrogen solvation in a total of 12 ionic liquids. Overall, the analysis of the structural properties under high temperature and pressure conditions revealed weak interactions between hydrogen and the ionic liquids, with a slight preference of this gas to be placed at the apolar domains. Interestingly, those ionic liquids comprising nitrate anions allow the accommodation of hydrogen molecules also in the polar areas. The study of the hydrogen velocity autocorrelation functions supports this observation. In addition, the structure of all of the tested ionic liquids was almost insensitive to the addition of hydrogen, so the available free volume and cavity formation are presumably the most important factors affecting solubility.