Using the method of molecular dynamics (MD), we have estimated the solubility of hydrogen in heavy hydrocarbons for a range of temperatures and pressures. Heavy hydrocarbon systems are known to be challenging not only for experimental measurements but also for reliable estimations using traditional equations of state (EOS). The simulation system used was designed with semi-permeable membranes to mimic actual experimental studies of gas solubility. Binary interaction parameters between the solute gas and the solvent (heavy hydrocarbons) components were adjusted when necessary to improve agreement with experimental results and then used in subsequent multi-component studies. Temperature and pressure ranges studied included higher temperatures and pressures which are especially difficult to investigate experimentally. Simulation results were finally used to adjust the binary interaction parameters (BIP) in simulation packages (e.g. Aspen) to enable quick and reliable predictions. (C) 2010 Elsevier B.V. All rights reserved.