Journal of Physical Chemistry B, Vol.122, No.42, 9763-9774, 2018
Molecular Origins of the Apparent Ideal CO2 Solubilities in Binary Ionic Liquid Mixtures
Molecular dynamics (MD) simulations were conducted to investigate the variation of Henry's constant of CO2 in two binary ionic liquid mixtures. One of the mixtures is formed by pairing the cation 1-n-butyl-3-methylimidazolium [C(4)mim](+) with chloride Cl- and methylsulfate [MeSO4](-), whereas the other binary ionic liquid mixture contains [C(4)mim](+) in combination with the anions Cl- and bis(trifluoromethanesulfonyl)imide [NTf2](-). In order to provide a microscopic understanding of the behavior of the Henry's constant with the anion composition, MD simulations of ionic liquid mixtures with and without CO2 saturation were performed at 353 K and 10 bar. Our calculations indicate that the Henry's constant for CO2 follows a highly nonlinear, although expected based on ideal solubility, trend with respect to the increasing concentration of Cl- in [C(4)mim]Cl-x[NTf2](1-x), whereas the Henry's constant is almost independent of the anion composition in the [C(4)mim]Cl-x[MeSO4](1-x) system. Structural analyses presented in terms of radial, spatial, and angular distribution functions point to significant structural reorganization of the anions around cations in the [C(4)mim]Cl-x[NTf2](1-x) system. Because of the weakly coordinating ability of the [NTf2](-) anion with the cation, the [NTf2](-) anion is displaced from the equatorial plane of the imidazolium ring and occupies positions above and below the ring, enabling enhanced CO2-[NTf2](-) association. The rearrangement also weakens the cation pi-pi interactions, resulting in the formation of increased local free volume aiding CO2 accommodation. On the contrary, such structural transitions are absent in the [C(4)mim]Cl-x[MeSO4](1-x) mixture system.