Theoretical constraints on the effects of pH, salinity, and temperature on clumped isotope signatures of dissolved inorganic carbon species and precipitating carbonate minerals

In order to determine the most reliable and efficient modeling methods to represent aqueous dissolved inorganic carbon (DIC) species and carbonate minerals, we performed convergence and sensitivity testing on several different levels of theory. We used 4 different techniques for modeling the hydration of DIC: gas phase, implicit solvation (PCM and SMD), explicit solvation (ion with 3 water molecules) and supermolecular clusters (ion plus 21 to 32 water molecules with geometries generated by molecular dynamics). For each solvation technique, we performed sensitivity testing by combining different levels of theory (up to 8 ab initio/hybrid methods, each with up to 5 different sizes of basis sets) to understand the limits of each technique. We looked at the degree of convergence with the most complex (and accurate) models in order to select the most reliable and efficient modeling methods. The B3LYP method combined with the 6-311++G(2d,2p) basis set with supermolecular clusters worked well.