Abstract:Advances in position-specific deuterium (2H) analyses provide powerful and accurate means for evaluating organic H isotopic signatures and fractionation factors. Density functional theory can aid in determining 2H–H equilibrium fractionation factors at specific positions on organic compounds. This work evaluates a DFT method that provides improved accuracy with respect to experimental data over previous methods while allowing for calculations on higher molecular weight compounds that are useful organic geochemical markers (biomarkers). We used computational quantum chemistry and applied the density functional theory method B3LYP and multiple basis sets. Based on accuracy criteria, we selected the basis set, 6-311++G(d,p), to calculate 2H–H equilibrium fractionation factors (α) on primary, secondary, tertiary, sp2, and aromatic C atoms for a model set that included methane, ethane, propane, 2-methylbutane, 2,3-dimethylpentane, 2,4-dimethylpentane, 2,6-dimethyloctane, 2,6,10,14-tetramethylpentadecane (pristane), 2,6,10,14-tetramethylhexadecane (phytane), twist-boat and chair cyclohexane, axial and equatorial methylcyclohexane, 2-heptanone, (E) and (Z) 2-pentene, and benzene with water. Two conformers each of cyclohexane, 2-pentene, and pristane were used to explore how thermodynamic weighting using the Boltzmann partition function could improve the α results for models of organic compounds that occur in more than one isomeric form. The B3LYP method coupled with the 6-311++G(d,p) basis set provided the most accurate α results, when compared with results from B3LYP coupled with 17 other basis sets. The calculated equilibrium α values for 2H on primary, secondary, and tertiary C atoms follow the observed trend in the deuterium preference for carbon–hydrogen bonds relative to water: α3° > α2° > α1°. Further, ln(β) values predicted from the DFT calculations also agree with fractionation estimated within 4%, as calculated using observed FTIR vibrational frequencies. Application of B3LYP/6-311++G(d,p) will aid in the assessments of position-specific 2H–H equilibrium fractionation factors (α) for comparison with experimental studies. Estimated position-specific alpha values can be scaled with temperature and used to predict molecular-averaged delta values, such as determined using compound-specific isotope analyses (CSIA). The ability to adjust position and molecular-averaged calculations for dominant isomeric forms makes this a useful approach to study equilibrium H isotope distributions in larger compounds of organic geochemical, biogeochemical, and environmental interest.

Performance of ab initio and density functional methods for conformational equilibria of CnH2n+2 alkane isomers (n=2-8)

Benchmark thermochemistry of the C_nH_{2n+2} alkane isomers (n=2--8) and performance of DFT and composite ab initio methods for dispersion-driven isomeric equilibria

Theoretical Studies on Stabilities and Spectroscopy of C80CH2

Barrier to asymmetric internal rotation, conformational stability, vibrational spectra and assignments, andAb Initio calculations ofn-butane-d0,d5 andd10

Assessment of Density Functionals on Intramolecular Dispersion Interaction in Large Normal Alkanes

What dictates alkane isomerization? A combined density functional theory and information-theoretic approach study

At what chain length do unbranched alkanes prefer folded conformations?

Study of Sterically Crowded Alkanes: Assessment of Non-Empirical Density Functionals Including Double-Hybrid (Cost-Effective) Methods

Simulating Vibrationally-Resolved X-ray Photoelectron Spectra of Flexible Molecules: Linear Alkanes C$_{n}$H$_{2n+2}$ ($n$=1-8)

Comparison of some dispersion-corrected and traditional functionals with CCSD(T) and MP2 ab initio methods: dispersion, induction, and basis set superposition error

Enthalpy Difference between Conformers of n-Butane and the Potential Function Governing Conformational Interchange

How Well Can B3LYP Heats of Formation Be Improved by Dispersion Correction Models?

Position-Specific 2 H/H Equilibrium Isotopic Fractionation Factors in Alkane, Alkene, and Aromatic Molecules: A Density Functional Theory Approach

Big data benchmarking: how do DFT methods across the rungs of Jacob's ladder perform for a dataset of 122k CCSD(T) total atomization energies?

Can G4-like Composite Ab Initio Methods Accurately Predict Vibrational Harmonic Frequencies?

Kinetics and Thermodynamics of Reactions Involving Criegee Intermediates: An Assessment of Density Functional Theory and Ab Initio Methods Through Comparison with CCSDT(Q)/CBS Data

Accurate Entropy Calculation for Large Flexible Hydrocarbons Using a Multi-Structural 2-Dimensional Torsion Method

Effect of molecular constitution and conformation on positron binding and annihilation in alkanes

A Theoretical Study on the Isomerization of Cyclopropane to Propene with Ab Initio and DFT Methods

Effect of Hindered Internal Rotation Treatments on Predicting the Thermodynamic Properties of Alkanes

Conformational stability, barriers to internal rotation, structural parameters, ab initio calculations, and vibrational assignment of cyclopropanecarboxaldehyde