Microbial production and consumption of methane are widespread in natural and artificial environments, with important economic and climatic implications. Attempts to use the isotopic composition of methane to constrain its sources are complicated by incomplete understanding of the mechanisms of variation in methane's isotopic composition. Knowledge of the equilibrium isotope fractionations among the large organic intracellular intermediates in the microbial pathways of methane production and consumption must form the basis of any exploration of the mechanisms of isotopic variation, but estimates of these equilibrium isotope fractionations are currently unavailable. To address this gap, we calculated the equilibrium isotopic fractionation of carbon (13 C/ 12 C) and hydrogen (D/H) isotopes among compounds in anaerobic methane metabolisms, as well as the abundance of multiple isotope substitutions ("clumping," e.g., 13 C-D) in these compounds. The Density Functional Theory calculations employed the M06-L/def2-TZVP level of theory and the SMD implicit solvation model, which we have recently optimized for large organic molecules and tested against measured equilibrium isotope fractionations. The computed 13 β and 2 β values decrease with decreasing average oxidation state of the carbon atom in the molecules, resulting in a preference for enrichment of the molecules with more oxidized carbon in 13 C and D. Using the computed β values, we calculated the equilibrium isotope fractionation factors in the prominent methanogenesis pathways (hydrogenotrophic, methylotrophic and acetoclastic) and in the pathway for anaerobic oxidation of methane (AOM) over a temperature range of 0-700 • C. Our calculated equilibrium fractionation factors compare favorably with experimental constrains, where available, and we used them to investigate the relation between the apparent isotope fractionation during