The Peng-Robinson equation of state (PR EOS) is used to model the vapor-liquid equilibria (VLE) of binary systems of hydrogen with five heavy-n-alkanes: n-decane (n-C10H22), n-hexadecane (n-C16H34), n-octacosane (n-C28H58), n-hexatriacontane (n-C36H74), and n-hexatetracontane (n-C46H94). Using literature experimental data for these systems, binary interaction parameters (BIPs) were calculated using the PR EOS coupled with three different alpha () functions (Soave, Twu, and Gasem). The calculated BIPs have been fitted to a generalized correlation that can be used to estimate the BIPs and model the VLE within the temperature range of 283.2-449.6K, pressure range of 1.151-15.970MPa, and hydrogen solubility range of 0.016-0.257mole fraction. It is found that PR EOS combined with one specific form of function is capable of reproducing the experimental VLE data with an overall %AARD of 1.1%.