Henry constant simulations employed the Monte Carlo (MC) method in the uniform ensemble, in which single guest molecules are inserted at random, and the total energies obtained for a sufficiently large ensemble of configurations are used to compute the Henry constant KH. The selectivity towards a binary mixture of two species A and B in the limit of infinite dilution (zero coverage) can be calculated from the Henry constants as: