Molecular simulation is used to compute the solubility of small gases in nitrile bu-
tadiene rubber (NBR) with a Widom particle-insertion technique biased by local free
volume. The convergence of the method is examined as a function of the number
of snapshots upon which the insertions are performed and the number of insertions
per snapshot, and is compared to the convergence of the unbiased Widom insertion
technique. The effect of varying the definition of the local free volume is also investi-
gated. The acrylonitrile content of the polymer is altered to examine its influence on
the solubility of helium, CO
2
and H
2
O, and the solubilities of polar gases are found to
be enhanced relative to nonpolar gases, in qualitative agreement with experiment. To
probe this phenomenon further, the solubilities are decomposed into contributions from
neighbourhoods of different atoms, using a Voronoi cell construction, and a strong bias
is found for CO
2
and H
2
O in particular to be situated near nitrogen sites in the elas-
tomer. Temperature is shown to suppress the solubility of CO
2
and H
2
O, but increase
that of helium. Increasing pressure is found to suppress the solubility of all gases, but
at different rates according to a balance between their molecular size and electrostatic
interaction with the polymer. These results are relevant to the use of NBR seals at
elevated temperatures and pressures, such as in oil and gas wells.