Comparison of Experimental vs Theoretical Abundances of 13CH3D and 12CH2D2 for Isotopically Equilibrated Systems from 1 to 500 °C
[1];
[4];
[2];
- Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States); Univ. of California, Berkeley, CA (United States)
- California Inst. of Technology (CalTech), Pasadena, CA (United States)
- Univ. of California, Berkeley, CA (United States)
- Princeton Univ., Princeton, NJ (United States)
Methane is produced and consumed via numerous microbial and chemical reactions in atmospheric, hydrothermal, and magmatic reactions. The stable isotopic composition of methane has been used extensively for decades to constrain the source of methane in the environment. A recently introduced isotopic parameter used to study the formation temperature and formational conditions of methane is the measurement of molecules of methane with multiple rare, heavy isotopes (“clumped”) such as 13CH3D and 12CH2D2. In order to place methaneclumped isotope measurements into a thermodynamic reference frame that allows calculations of clumped isotope-based temperatures (geothermometry) and comparison between laboratories, all past studies have calibrated their measurements using a combination of experiment and theory based on the temperature dependence of clumped isotopologue distributions for isotopically equilibrated systems. These have previously been performed at relatively high temperatures (>150 °C). Given that many natural occurrences of methane form below these temperatures, previous calibrations require extrapolation when calculating clumped isotope-based temperatures outside of this calibration range. We provide a new experimental calibration of the relative equilibrium abundances of 13CH3D and 12CH2D2 from 1 to 500 °C using a combination of γ-Al2O3- and Ni-based catalysts and compare them to new theoretical computations using Path Integral Monte Carlo (PIMC) methods and find 1:1 agreement (within ±1 standard error) for the observed temperature dependence of clumping between experiment and theory over this range. This demonstrates that measurements, experiments, and theory agree from 1 to 500 °C, providing confidence in the overall approaches. Polynomial fits to PIMC computations, which are considered the most rigorous theoretical approach available, are given as follows (valid T ≥ 270 K): Δ13CH3D ≅ 1000 × ln(K13CH3D) = (1.47348 × 1019)/T7 – (2.08648 × 1017)/T6 + (1.19810 × 1015)/T5 – (3.54757 × 1012)/T4 + (5.54476 × 109)/T3 – (3.49294 × 106)/T2 + (8.89370 × 102)/T and Δ12CH2D2 ≅ 1000 × ln(8/3K12CH2D2) = –(9.67634 × 1015)/T6 + (1.71917 × 1014)/T5 – (1.24819 × 1012)/T4 + (4.30283 × 109)/T3 – (4.48660 × 106)/T2 + (1.86258 × 103)/T. We additionally compare PIMC computations to those performed utilizing traditional approaches that are the basis of most previous calibrations (Bigeleisen, Mayer, and Urey model, BMU) and discuss the potential sources of error in the BMU model relative to PIMC computations.
- Research Organization:
- Lawrence Berkeley National Laboratory (LBNL), Berkeley, CA (United States)
- Sponsoring Organization:
- USDOE Office of Science (SC), Basic Energy Sciences (BES); National Science Foundation (NSF)
- Grant/Contract Number:
- AC02-05CH11231; EAR-1911296; CHE-1611581
- OSTI ID:
- 1575296
- Alternate ID(s):
- OSTI ID: 1580033
- Journal Information:
- ACS Earth and Space Chemistry, Vol. 3, Issue 12; ISSN 2472-3452
- Publisher:
- American Chemical SocietyCopyright Statement
- Country of Publication:
- United States
- Language:
- English
Web of Science
Similar Records
Clumped 13CH2D and 12CHD2 compositions of methyl groups from wood and synthetic monomers: Methods, experimental and theoretical calibrations, and initial results
Measurements of stable isotope ratios ({sup 13}CH{sub 4}/{sup 12}CH{sub 4}; {sup 12}CH{sub 3}/{sup 12}CH{sub 3}) in landfill methane using a tunable diode laser absorption spectrometer