Stable isotope ratios of C1–C5 alkanes, the major constituents of subsurface gaseous hydrocarbons, can provide valuable insights on their origins, transport, and fates. Equilibrium isotope effects are fundamental to interpreting stable isotope signatures, as recognition of them in natural materials indicates reversible processes and constrains the temperatures of equilibrated systems. Hydrogen isotope equilibrium of C1–C5 alkanes is of particular interest because evidence shows that alkyl H can undergo isotopic exchange with coexisting compounds under subsurface conditions. We present the results of a combined experimental and theoretical effort to determine equilibrium hydrogen isotope distributions in mixtures of these hydrocarbon compounds. We created two mixtures: one with C1, C2 and C3 (where C1 indicates methane, C2 ethane, etc., hereinafter called ‘C1–C3 mixture’) and another one with C2, C3, iC4, nC4, iC5 and nC5 (hereinafter called ‘C2–C5 mixture’); in both cases, the mixtures were created to start out of hydrogen isotope equilibrium. We tested the performance of several metal catalysts as aids to H isotopic exchange by exposing these mixtures to different metal catalysts at 100 or 200 °C and analyzing the compound-specific hydrogen isotope ratios of the product gases. The C1–C3 mixture exchanged hydrogen isotopes among the starting gas components rapidly in the presence of Ru/Al2O3 at 200 °C. The isotope ratios reach a steady-state (presumably equilibrium) after 72 h of heating (up to 120 h). The hydrogen isotopic ratios of alkanes in the C2–C5 mixture shifted significantly in the presence of Rh/Al2O3 at 100 °C and C3-C5 compounds approached steady state after 70 h, whereas C2 had not yet reached a steady state D/H ratio after 216 h of heating. We evaluated the reaction progress of isotope exchange for each compound as a function of time with a reaction-network model informed by constraints on chemical kinetics. The model indicates that C3, iC4, nC4, iC5 and nC5 in the C2-C5 mixture reached internal isotope equilibrium after 70 h, hence we used their values to calculate experimental equilibrium results. We also calculated equilibrium isotope fractionations with the Bigeleisen-Mayer theorem using vibrational frequencies computed from the density functional theory (B3LYP/aug-cc-pVTZ). We included torsional conformers and explicit H positions in the calculations. We found that the experimental results from both the 200 ˚C C1–C3 experiment and 100 ˚C C2–C5 experiment are consistent with theoretical equilibrium results within experimental uncertainty. Additionally, our reaction-network model for catalyzed hydrogen isotope exchange between alkanes succeeded in fitting our experimental results. This modeling framework can be adjusted to simulate exchange in natural settings for constraining temperature histories and sources of natural gases.