Thermodynamics of Metal Carbonates and Bicarbonates and Their Hydrates for Mg, Ca, Fe, and Cd Relevant to Mineral Energetics.

The heats of formation of the carbonate, bicarbonate, and bicarbonate/hydroxide metal complexes, including hydrates of Mg2+, Ca2+, Fe2+, and Cd2+, and the oxides, dichlorides, and dihydroxides are predicted from atomization energies using correlated molecular orbital theory at the CCSD(T) level extrapolated to the complete basis set limit following the Feller-Peterson-Dixon (FPD) approach. Using the calculated gas phase values and the available experimental solid-state values, we predicted the cohesive energies of selective minerals. The gas phase decomposition energies of MO, CO2, and H2O follow the order Mg ≈ Ca > Cd ≈ Fe and correlate with the hardness of the metal +2 ions. Gas phase hydration energies show that the order is Mg > Fe > Ca ≈ Cd. There are a number of bulk hydrated Mg and Ca complexes that occur as minerals but there are few if any for Fe and Cd, suggesting that a number of factors are important in determining the stability of the bulk mineral hydrates. The FPD heats of formation were used to benchmark a range of density functional theory exchange-correlation functionals, including those commonly used in solid-state mineral calculations. None of the functionals provided chemical accuracy agreement (±1 kcal/mol) with the FPD results. The best agreement to the FPD results is predicted for ωB97X and ωB97X-D functionals with an average unsigned error of 10 kcal/mol. The worst functionals are PW91, BP86, and PBE with average unsigned errors of 32-36 kcal/mol.

In Situ Natural Abundance 17O and 25Mg NMR Investigation of Aqueous Mg(OH)2 Dissolution in the Presence of Supercritical CO2.

We report an in situ high-pressure NMR capability that permits natural abundance 17O and 25Mg NMR characterization of dissolved species in aqueous solution and in the presence of supercritical CO2 fluid (scCO2). The dissolution of Mg(OH)2 (brucite) in a multiphase water/scCO2 fluid at 90 atm pressure and 50 °C was studied in situ, with relevance to geological carbon sequestration. 17O NMR spectra allowed identification and distinction of various fluid species including dissolved CO2 in the H2O-rich phase, scCO2, aqueous H2O, and HCO3-. The widely separated spectral peaks for various species can all be observed both dynamically and quantitatively at concentrations as low as 20 mM. Measurement of the concentrations of these individual species also allows an in situ estimate of the hydrogen ion concentration, or pCH+ values, of the reacting solutions. The concentration of Mg2+ can be observed by natural abundance 25Mg NMR at a concentration as low as 10 mM. Quantum chemistry calculations of the NMR chemical shifts on cluster models aided in the interpretation of the experimental results. Evidence for the formation of polymeric Mg2+ clusters at high concentrations in the H2O-rich phase, a possible critical step needed for magnesium carbonate formation, was found.

Prediction of the pKa's of aqueous metal ion +2 complexes.

Aqueous metal ions play an important role in many areas of chemistry. The acidities of [Be(H2O)4](2+), [M(H2O)6](2+), M = Mg(2+), Mn(2+), Fe(2+), Co(2+), Ni(2+), Cu(2+), Zn(2+), Cd(2+), and Hg(2+), and [M(H2O)n](2+), M = Ca(2+) and Sr(2+), n = 7 and 8, complexes have been predicted using density functional theory, second-order Møller-Plesset perturbation theory (MP2), and coupled cluster CCSD(T) theory in the gas phase. pKa's in aqueous solution were predicted by using self-consistent reaction field (SCRF) calculations with different solvation models. The most common binding motif of the majority of the metal +2 complexes is coordination number (CN) 6, with each hexaaquo cluster having reasonably high symmetry for the best arrangement of the water molecules in the first solvation shell. Be(2+) is tetracoordinated, but a second solvation shell of 8 waters is needed to predict the pKa. The Ca(2+) and Sr(2+) aquo clusters have a coordination number of 7 or 8 as found in terms of the energy of the reaction M(H2O)7(2+) + H2O → M(H2O)8(2+) and the pKa values. The calculated geometries are in reasonable agreement with experiment. The SCRF calculations with the conductor-like screening model (COSMO), and the conductor polarized continuum model (CPCM) using COSMO-RS radii, consistently agree best with experiment at the MP2/aug-cc-pVDZ and CCSD(T)/aug-cc-pVDZ levels of theory. The CCSD(T) level provides the most accurate pKa's, and the MP2 level also provides reliable predictions. Our predictions were used to elucidate the properties of metal +2 ion complexes. The pKa predictions provide confirmation of the size of the first solvation shell sizes. The calculations show that it is still difficult to predict pKa's using this cluster/implicit solvent approach to better than 1 pKa unit.

Structures and energetics of (MgCO3)n clusters (n ≤ 16).

There is significant interest in the role of carbonate minerals for the storage of CO2 and the role of prenucleation clusters in their formation. Global minima for (MgCO3)n (n ≤ 16) structures were optimized using a tree growth-hybrid genetic algorithm in conjunction with MNDO/MNDO/d semiempirical molecular orbital calculations followed by density functional theory geometry optimizations with the B3LYP functional. The most stable isomers for (MgCO3)n (n < 5) are approximately 2-dimensional. Mg can be bonded to one or two O atoms of a CO3(2-), and the 1-O bonding scheme is more favored as the cluster becomes larger. The average C-Mg coordination number increases as the cluster size increases, and at n = 16, the average C-Mg coordination number was calculated to be 5.2. The normalized dissociation energy to form monomers increases as n increases. At n = 16, the normalized dissociation energy is calculated to be 116.2 kcal/mol, as compared to the bulk value of 153.9 kcal/mol. The adiabatic reaction energies for the recombination reactions of (MgO)n clusters and CO2 to form (MgCO3)n were calculated. The exothermicity of the normalized recombination energy ⟨RE⟩CO2 decreases as n increases and converged to the experimental bulk limit rapidly. The normalized recombination energy ⟨RE⟩CO2 was calculated to be -52.2 kcal/mol for the monomer and -30.7 kcal/mol for n = 16, as compared to the experimental value of -27.9 kcal/mol for the solid phase reaction. Infrared spectra for the lowest energy isomers were calculated, and absorption bands in the previous experimental infrared studies were assigned with our density functional theory predictions. The (13)C, (17)O, and (25) Mg NMR chemical shifts for the clusters were predicted. The results provide insights into the structural and energetic transitions from nanoclusters of (MgCO3)n to the bulk and the spectroscopic properties of clusters for their experimental identification.

Density functional theory study of the complexation of the uranyl dication with anionic phosphate ligands with and without water molecules.

The structures, vibrational frequencies and energetics of anhydrous and hydrated complexes of UO2(2+) with the phosphate anions H2PO4(-), HPO4(2-), and PO4(3-) were predicted at the density functional theory (DFT) and MP2 molecular orbital theory levels as isolated gas phase species and in aqueous solution by using self-consistent reaction field (SCRF) calculations with different solvation models. The geometries and vibrational frequencies of the major binding modes for these complexes are compared to experiment where possible and good agreement is found. The uranyl moiety is nonlinear in many of the complexes, and the coordination number (CN) 5 in the equatorial plane is the predominant binding motif. The phosphates are found to bind in both monodentate and bidentate binding modes depending on the charge and the number of water molecules. The SCRF calculations were done with a variety of approaches, and different SCRF approaches were found to be optimal for different reaction types. The acidities of HxPO4(3-x) in HxPO4(3-x)(H2O)4, x = 0-3 complexes were calculated with different SCRF models and compared to experiment. Phosphate anions can displace water molecules from the first solvation shell at the uranyl exothermically. The addition of water molecules can cause the bonding of H2PO4(-) and HPO4(2-) to change from bidentate to monodentate exothermically while maintaining CN 5. The addition of water can generate monodentate structures capable of cross-linking to other uranyl phosphates to form the types of structures found in the solid state. [UO2(HPO4)(H2O)3] is predicted to be a strong base in the gas phase and in aqueous solution. It is predicted to be a much weaker acid than H3PO4 in the gas phase and in solution.

Thermochemical properties of selenium fluorides, oxides, and oxofluorides.

The bond dissociation energies (BDEs), fluoride and fluorocation affinities, and electron affinities of SeF(n) (n = 1-6), SeOF(n) (n = 0-4), and SeO(2)F(n) (n = 0-2) have been predicted with coupled cluster CCSD(T) theory extrapolated to the complete basis set limit. To achieve near chemical accuracy, additional corrections were added to the complete basis set binding energies based on frozen core coupled cluster theory energies. These included corrections for core-valence effects, scalar relativistic effects, for first-order atomic spin-orbit effects, and vibrational zero point energies. The adiabatic BDEs contain contributions from product reorganization energies and, therefore, can be much smaller than the diabatic BDEs and can vary over a wide range. For thermochemical calculations, the adiabatic values must be used, whereas for bond strength and kinetic considerations, the diabatic values should be used when only small displacements of the atoms without change of the geometry of the molecule are involved. The adiabatic Se-F BDEs of SeF(n) (n = 1-6) are SeF(6) = 90, SeF(5) = 27, SeF(4) = 93, SeF(3) = 61, SeF(2) = 86, and SeF = 76 kcal/mol, and the corresponding diabatic values are SeF(6) = 90, SeF(5) = 88, SeF(4) = 93, SeF(3) = 74, SeF(2) = 86, and SeF = 76 kcal/mol. The adiabatic Se-O BDEs of SeO(n) (n = 1-3), SeOF(n) (n = 1-4), and SeO(2)F(n) (n = 1,2) range from 23 to 107 kcal/mol, whereas the diabatic ones range from 62 to 154 kcal/mol. The adiabatic Se-F BDEs of SeOF(n) (n = 1-4) and SeO(2)F(n) (n = 1,2) range from 20 to 88 kcal/mol, whereas the diabatic ones range from 73 to 112 kcal/mol. The fluoride affinities of SeF(n), (n = 1-6), SeO(n), (n = 1-3), SeOF(n), (n = 1-4), and SeO(2)F(n) (n = 1,2) range from 15 to 121 kcal/mol, demonstrating that the Lewis acidity of these species covers the spectrum from very weak (SeF(6)) to very strong (SeO(3)) acids. The electron affinities which are a measure of the oxidizing power of a species, span a wide range from 1.56 eV in SeF(4) to 5.16 eV in SeF(5) and for the free radicals are much higher than for the neutral molecules. Another interesting feature of these molecules and ions stems from the fact that many of them possess both a Se free valence electron pair and a free unpaired valence electron, raising the questions of their preferred location and their influence on the Se-F and Se═O bond strengths.

Methane to methanol conversion induced by thorium oxide through the CH3Th(O)H intermediate in solid argon.

Reactions of ThO molecules and CH(4) have been investigated in solid argon near 4 K. The CH(3)Th(O)H molecule is produced when the sample is exposed to UV irradiation. Identification of this new intermediate is substantiated by observation of the Th═O and Th-H stretching vibrational modes with isotopic substitution via matrix infrared spectroscopy, and the assignments are supported by electronic structure frequency calculations. Methanol absorptions increase together with formation of the CH(3)Th(O)H molecule, suggesting a methane to methanol conversion induced by thorium oxide proceeding through the CH(3)Th(O)H intermediate. The formation of CH(3)Th(O)H from ThO + CH(4) is exothermic (ΔH(rxn) = -11 kcal/mol) with an energy barrier of 30 kcal/mol at the CCSD(T)//B3LYP level. Decomposition of this intermediate to form methanol involves spin crossing, and the overall reaction from the intermediate is endothermic by 127 kcal/mol. There is no activation energy for the reaction of thorium atoms with methanol to give CH(3)Th(O)H, as observed in separate experiments with Th and CH(3)OH.

High-level ab initio predictions of the energetics of mCO2·(H2O)n (n = 1-3, m = 1-12) clusters.

Electronic structure calculations at the correlated molecular orbital theory and density functional theory levels have been used to generate a reliable set of clustering energies for up to three water molecules in carbon dioxide clusters up to n = 12. The structures and energetics are dominated by Lewis acid-base interactions with hydrogen-bonding interactions playing a lesser energetic role. The actual binding energies are somewhat larger than might be expected. The correlated molecular orbital MP2 method and density functional theory with the ωB97X exchange-correlation functional provide good results for the energetics of the clusters, but the B3LYP and ωB97X-D functionals do not. Seven CO(2) molecules form the first solvent shell about a single H(2)O with four CO(2) molecules interacting with the H(2)O via Lewis acid-base interactions, two CO(2) interacting with the H(2)O by hydrogen bonds, and the seventh CO(2) completing the shell. The Lewis acid-base and weak hydrogen bond interactions between the water molecules and the CO(2) molecules are strong enough to disrupt the trimer ring configuration for as few as seven CO(2) molecules. Calculated (13)C NMR chemical shifts for mCO(2)·(H(2)O)(n) show little change with respect to the number of H(2)O or CO(2) molecules in the cluster. The O-H stretching frequencies do exhibit shifts that can provide information about the interactions between water and CO(2) molecules.

Matrix infrared spectra and theoretical studies of thorium oxide species: ThOx and Th2Oy.

Infrared spectra of three new thorium oxide species have been obtained in argon and neon matrixes. All of the products are experimentally characterized using isotopic oxygen samples with the aid of electronic structure calculations. Ground state thorium atoms react with O(2) to form the ThO(2) molecules, which can dimerize to give Th(2)O(4) products. Th(2)O(4) is predicted to have nonplanar C(2h) symmetry for its closed shell singlet ground state. The rhombus-shaped Th(2)O(2) molecule in the (1)A(g) (D(2h)) ground state is also observed and its formation is proposed via the reaction of Th(2) with O(2). In addition, electron capture of neutral thorium dioxide results in the formation of the ThO(2)(-) anion. It is predicted to have a doublet ground state with a geometry similar to that of the neutral ThO(2) molecule. Electronic structure calculations on the unobserved Th(2)O and Th(2)O(3) molecules are also provided.

Prediction of vibrational frequencies of UO2(2+) at the CCSD(T) level.

Electronic structure calculations at the coupled cluster (CCSD(T)) and density functional theory levels with relativistic effective core potentials and large basis sets were used to predict the isolated uranyl ion frequencies. The effects of anharmonicity and spin-orbit corrections on the harmonic frequencies were calculated. The anharmonic effects are larger than the spin-orbit corrections, but both are small. The anharmonic effects decreased all the frequencies, whereas the spin-orbit corrections increased the stretches and decreased the bend. Overall, these two corrections decreased the harmonic asymmetric stretch frequency by 6 cm-1, the symmetric stretch by 3 cm-1, and the bend by 3 cm-1. The best calculated values for UO22+ for the asymmetric stretch, symmetric stretch, and bend were 1113, 1032, and 174 cm-1, respectively. The separation between the asymmetric and the symmetric stretch band origins was predicted to be 81 cm-1, which is consistent with experimental trends for substituted uranyls in solution and in the solid state. The anharmonic vibrational frequencies of the isoelectronic ThO2 molecule also were calculated and compared to experiment to calibrate the UO22+ results.

Mechanism of the hydration of carbon dioxide: direct participation of H2O versus microsolvation.

Thermochemical parameters of carbonic acid and the stationary points on the neutral hydration pathways of carbon dioxide, CO 2 + nH 2O --> H 2CO 3 + ( n - 1)H 2O, with n = 1, 2, 3, and 4, were calculated using geometries optimized at the MP2/aug-cc-pVTZ level. Coupled-cluster theory (CCSD(T)) energies were extrapolated to the complete basis set limit in most cases and then used to evaluate heats of formation. A high energy barrier of approximately 50 kcal/mol was predicted for the addition of one water molecule to CO 2 ( n = 1). This barrier is lowered in cyclic H-bonded systems of CO 2 with water dimer and water trimer in which preassociation complexes are formed with binding energies of approximately 7 and 15 kcal/mol, respectively. For n = 2, a trimeric six-member cyclic transition state has an energy barrier of approximately 33 (gas phase) and a free energy barrier of approximately 31 (in a continuum solvent model of water at 298 K) kcal/mol, relative to the precomplex. For n = 3, two reactive pathways are possible with the first having all three water molecules involved in hydrogen transfer via an eight-member cycle, and in the second, the third water molecule is not directly involved in the hydrogen transfer but solvates the n = 2 transition state. In the gas phase, the two transition states have comparable energies of approximately 15 kcal/mol relative to separated reactants. The first path is favored over in aqueous solution by approximately 5 kcal/mol in free energy due to the formation of a structure resembling a (HCO 3 (-)/H 3OH 2O (+)) ion pair. Bulk solvation reduces the free energy barrier of the first path by approximately 10 kcal/mol for a free energy barrier of approximately 22 kcal/mol for the (CO 2 + 3H 2O) aq reaction. For n = 4, the transition state, in which a three-water chain takes part in the hydrogen transfer while the fourth water microsolvates the cluster, is energetically more favored than transition states incorporating two or four active water molecules. An energy barrier of approximately 20 (gas phase) and a free energy barrier of approximately 19 (in water) kcal/mol were derived for the CO 2 + 4H 2O reaction, and again formation of an ion pair is important. The calculated results confirm the crucial role of direct participation of three water molecules ( n = 3) in the eight-member cyclic TS for the CO 2 hydration reaction. Carbonic acid and its water complexes are consistently higher in energy (by approximately 6-7 kcal/mol) than the corresponding CO 2 complexes and can undergo more facile water-assisted dehydration processes.

Fundamental thermochemical properties of amino acids: gas-phase and aqueous acidities and gas-phase heats of formation.

The gas-phase acidities of the 20 L-amino acids have been predicted at the composite G3(MP2) level. A broad range of structures of the neutral and anion were studied to determine the lowest energy conformer. Excellent agreement is found with the available experimental gas-phase deprotonation enthalpies, and the calculated values are within experimental error. We predict that tyrosine is deprotonated at the CO(2)H site. Cysteine is predicted to be deprotonated at the SH but the proton on the CO(2)H is shared with the S(-) site. Self-consistent reaction field (SCRF) calculations with the COSMO parametrization were used to predict the pK(a)'s of the non-zwitterion form in aqueous solution. The differences in the non-zwitterion pK(a) values were used to estimate the free energy difference between the zwitterion and nonzwitterion forms in solution. The heats of formation of the neutral compounds were calculated from atomization energies and isodesmic reactions to provide the first reliable set of these values in the gas phase. Further calculations were performed on five rare amino acids to predict their heats of formation, acidities, and pK(a) values.

Quantum-chemical calculations of carbon-isotope fractionation in CO2(g), aqueous carbonate species, and carbonate minerals.

Quantum chemical calculations on large supermolecular carbonate-water and carbonate mineral clusters are used to predict equilibrium constants for 13,12C-isotope-exchange reactions between CO2(g), aqueous carbonate species, and the common carbonate minerals. For the aqueous species, we evaluate the influence of the size and conformational variability of the solvation shell, the exchange-correlation functional, and the basis set. The choice of exchange-correlation functional (PBE vs B3LYP), the basis set (6-31G* vs aug-cc-pVDZ), and solvation shell size (first shell only vs first shell and a partial second shell) each produce changes of approximately 5-10 per mil in the reduced partition function ratio. Conformational variability gives rise to a standard error of approximately 0.5 per mil using approximately 10 solute-solvent conformations. The best results are obtained with the B3LYP/aug-cc-pVDZ combination, but because the improvements in the basis set and exchange correlation functional drive the reduced partition function ratios in opposite directions, reasonably good results are also obtained with the PBE/6-31G* combination. To construct molecular clusters representative of mineral environments, a new method is introduced on the basis of conservation of Pauling bond strength. Using these clusters as models for minerals, calculations of mineral-gas and mineral-aqueous carbon-isotope fractionation factors, are in good agreement with experimental measurements. Carbon-isotope fractionation factors for gas, aqueous, and mineral phases are thus integrated into a single theoretical/computational framework.