The distribution of stable isotopes in organic compounds within geochemical systems reflects the reactions they undergo, which are often kinetically governed and can be elucidated through quantification based on first principles: employing quantum chemical simulations to derive energy barriers and kinetic isotope effects (KIEs), followed by kinetic calculations using these parameters. In this paper, the distribution of 13C in hydrocarbon gases during the thermal generation and decomposition of ethane and propane has been addressed. The work is based on abundant and systematic isotopic data, including the compound- and position-specific 13C compositions, recent improvements in kinetic networks in hydrocarbon generation, and state-of-the-art quantum chemical methods that can process organic geochemical systems involving solid phases. Density Functional Theory computations validate the carbonium path of hydrocarbon generation and the free-radical path of their decomposition under geological conditions. The first step of each path (alkyl sidechain protonation for generation and hydrogen abstraction for decomposition) is the rate-limiting step and determines the 13C fractionations of gaseous hydrocarbons, with the reaction rate constrained by the concentration or accessibility of the active species (hydrated protons and surface free radicals in the two paths, respectively). Alkyl protonation has normal 13C KIEs, but hydrogen abstraction on the surfaces of solid organic matter has inverse ones (substituting 12C by 13C accelerates the reaction). These inverse KIEs explain three isotopic phenomena in gaseous hydrocarbons in their late thermal evolution: 1) the depletion of 13C in ethane and propane accompanied by their decreased fraction in hydrocarbon gases (isotope “rollover”); 2) the reversed 13C sequence with carbon number (δ13Cmethane > δ13Cethane > δ13Cpropane); and 3) the depletion of 13C in the center carbon atoms of residual propane. Isotopic fractionation between gaseous hydrocarbons upon generation is evaluated using the Successive Random Chain-Scission model. The covariation of δ13C in alkanes shows a depletion in CH4 compared with the values expected from the reciprocal of the carbon numbers, indicating that this deviation may be intrinsic instead of mixing with microbial methane. The comparison of intramolecular 13C fractionation of propane between computational and observed data suggests the contribution of isoprenoids at the beginning of the generation, in addition to isotopic homogenization between the center and end carbon atoms through the ring closure of intermediate cations. The 13C–13C clumping exhibits a weak KIE in both the generation and decomposition of ethane; the clumping is essentially stochastically governed. This work demonstrates how first-principle calculations can provide significant insights into the details of geochemical reactions governing the isotopic distributions.
Read full abstract