Ab Initio Prediction of the Gas- and Solution-Phase Acidities of Strong Brønsted Acids: The Calculation of pKa Values Less Than −10

Prediction of Redox Potentials for U, Np, Pu, and Am in Aqueous Solution

The redox properties of the actinides in aqueous solution are important for fuel production/reprocessing and understanding the environmental impact of nuclear waste. The redox potentials for U, Np, Pu, and Am in oxidation states from 0 up to VII (as appropriate) in aqueous solutions have been predicted at the density functional theory level with the B3LYP functional, Stuttgart small core pseudopotential basis sets for the actinides, and explicit (30H2O molecules)/implicit treatment of the aqueous solvent using the self-consistent reaction field COSMO and SMD approaches for the implicit solvation. The predictions of the structural parameters of clusters incorporating first and second solvation shells are consistent with the available experimental data. Our results are typically within 0.2 V of the available experimental data using two explicit solvation shells with an implicit solvent model. The use of the PW91 functional substantially improved the prediction of the Pu(VI/V) redox couple. The redox couples for An(VI/IV) and An(V/IV) which involve the addition of protons and removal of the actinyl oxygens led to slightly larger differences from an experiment. The An(IV/0) and An(III/0) couples were reliably predicted with our approach. Predictions of the unknown An(II/I) redox potentials were negative, consistent with expectations, and predictions for unknown An(VII/VI), An(III/II), and An(II/0) redox couples improve prior estimates.

Prediction of the pKa's of Hydrated Metal Carbonates and Bicarbonates for Mg, Ca, Mn, Fe, Co, Ni, Cu, and Zn Dications

The gas- and aqueous-phase acidities of hydrated metal dication carbonates, bicarbonates, and hydroxide complexes M(CO3)(H2O)n for n = 1 to 3, M(HCO3)2, M(HCO3)2(H2O)2, M(HCO3)(OH), and M(HCO3)(H2O)2(OH) for M = Mg, Ca, Mn, Fe, Co, Ni, Cu, and Zn were calculated at the CCSD(T)/aug-cc-pwCVDZ/cc-pwCVDZ level in the gas phase and at the B3LYP/aug-cc-pVTZ/cc-pVTZ(-PP) level with the COSMO self-consistent reaction field (SCRF) method in the aqueous phase. The composite correlated molecular orbital theory G3(MP2) and G3(MP2)B3 methods were used to predict the pKa's of the Mg structures and cis-cis carbonic acid to provide additional benchmarks. Using values scaled to experiment for H2CO3, the pKa's of bicarbonate ligands in group 2 and transition-metal complexes were compared to carbonic acid to gauge the effect of the metal complex on the bicarbonate. The group 2 metal complexes M(HCO3)2 and M(HCO3)(OH) decreased the acidity of the bicarbonate ligands, whereas their dihydrates were even less acidic. The transition-metal di-bicarbonate and bicarbonate hydroxide complexes generally made the bicarbonate more acidic especially when reduction of the metal occurs consistent with electron donation from the ligands; this is accompanied by spin transfer which typically increases in the order Mn < Fe < Co < Ni < Cu. The transition-metal dihydrates were less acidic than carbonic acid. Using values scaled to experiment for hydrated metal dications, the pKa's of water coordinated to group 2 and transition-metal complexes were generally more acidic than the hydrated metal dications, with the exception of Ca bicarbonate dihydrate, Co carbonate, Ni di-bicarbonate dihydrate, and Cu bicarbonate hydroxide di-bicarbonate.

Volatilization Potential of Per- and Poly-fluoroalkyl Substances from Airfield Pavements and during Recycling of Asphalt

Per- and poly-fluoroalkyl substances (PFAS) in water are typically present in their ionic (nonvolatile) forms; however, these can transition to their nonionic (volatile) forms when in contact with organic solvents and organic matrices. In particular, when PFAS are dissolved in organic solvents such as residues left from firefighting foams, fuels, and bitumen present in asphalt, the equilibrium between ionic and nonionic forms can trend toward more volatile nonionic forms of PFAS. We assessed the volatility of common PFAS based on calculated and available experimental data across ambient temperature ranges experienced by airfield pavements and at elevated temperatures associated with reworking asphalts for reuse. Volatilities are shown to be comparable to hydrocarbons in the semivolatile range, suggesting that volatilization is a viable loss mechanism for some PFAS that are nonvolatile in water. The present study points to future investigative needs for this unexplored mass loss mechanism and potential exposure pathway. Environ Toxicol Chem 2022;41:2202-2208. © 2022 Commonwealth of Australia. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Predicting the Mechanism and Products of CO2 Capture by Amines in the Presence of H2O

An extensive correlated molecular orbital theory study of the reactions of CO₂ with a range of substituted amines and H₂O in the gas phase and aqueous solution was performed at the G3(MP2) level with a self-consistent reaction field approach. The G3(MP2) calculations were benchmarked at the CCSD(T)/CBS level for NH₃ reactions. A catalytic NH₃ reduces the energy barrier more than a catalytic H₂O for the formation of H₂NCOOH and H₂CO₃. In aqueous solution, the barriers to form both H₂NCOOH and H₂CO₃ are reduced, with HCO₃^- formation possible with one amine present and H₂NCOO^- formation possible only with two amines. Further reactions of H₂NCOOH to form HNCO and urea via the Bazarov reaction have high barriers and are unlikely in both the gas phase and aqueous solution. Reaction coordinates for CH₃NH₂, CH₃CH₂NH₂, (CH₃)₂NH, CH₃CH₂CH₂NH₂, (CH₃)₃N, and DMAP were also calculated. The barrier for proton transfer correlates with amine basicity for alkylammonium carbamate (ΔG^‡_aq < 15 kcal/mol) and alkylammonium bicarbonate (ΔG^‡_aq < 30 kcal/mol) formation. In aqueous solution, carbamic acids, carbamates, and bicarbonates can all form in small amounts with ammonium carbamates dominating for primary and secondary alkylamines. These results have implications for CO₂ capture by amines in both the gas phase and aqueous solution as well as in the solid state, if enough water is present.

Bond Dissociation Energies in Heavy Element Chalcogen and Halogen Small Molecules

Thermodynamic properties including bond dissociation energies (BDEs), heats of formation, and gas-phase acidities for the hydrides and dimers of chalcogens and halogens, H₂Y, HX, Y₂, and X₂ for Y = Se, Te, and At and X = Br, I, and At, have been predicted using the Feller-Peterson-Dixon composite-correlated molecular orbital theory approach. A full four-component CCSD(T) approach was used to calculate the spin-orbit effects on thermodynamic properties, except for Se₂, where the AoC-DHF value was used due to strong multireference effects in Se₂ for the SO calculations. The calculated results show that the At₂ BDE is quite small, 19.5 kcal/mol, with much of the low bond energy due to spin-orbit effects. H₂Po is not predicted to be stable to dehydrogenation to Po + H₂ in terms of the free energy at 298 K. In the gas phase, HAt is predicted to be a stronger acid than H₂SO₄. The current results provide insights into potential difficulties in the actual experimental observation of such species for heavy elements.

The Effect of Intramolecular Hydrogen Bond Type on the Gas-Phase Deprotonation of ortho-Substituted Benzenesulfonic Acids. A Density Functional Theory Study

Structural factors have been identified that determine the gas-phase acidity of ortho-substituted benzenesulfonic acid, 2-XC6H4-SO3H, (X = -SO3H, -COOH, -NO2, -SO2F, -C≡N, -NH2, -CH3, -OCH3, -N(CH3)2, -OH). The DFT/B3LYP/cc-pVTZ method was used to perform conformational analysis and study the structural features of the molecular and deprotonated forms of these compounds. It has been shown that many of the conformers may contain anintramolecular hydrogen bond (IHB) between the sulfonic group and the substituent, and the sulfonic group can be an IHB donor or an acceptor. The Gibbs energies of gas-phase deprotonation ΔrG0298 (kJ mol-1) were calculated for all compounds. It has been set that in ortho-substituted benzenesulfonic acids, the formation of various types of IHB is possible, having a significant effect on the ΔrG0298 values of gas-phase deprotonation. If the -SO3H group is the IHB donor, then an ion without an IHB is formed upon deprotonation, and the deprotonation energy increases. If this group is an IHB acceptor, then a significant decrease in ΔrG0298 of gas-phase deprotonation is observed due to an increase in IHB strength and the A- anion additional stabilization. A proton donor ability comparative characteristic of the -SO3H group in the studied ortho-substituted benzenesulfonic acids is given, and the ΔrG0298 energies are compared with the corresponding values of ortho-substituted benzoic acids.

Icosahedral Carborane Superacids and their Conjugate Bases Comprising H, F, Cl, and CN Substituents: A Theoretical Investigation of Monomeric and Dimeric Cages

Theoretical investigation of the H(CHB₁₁ X₁₁ ) (X=H, F, Cl, CN), H(CHB₁₁ X_n Y_11-n ) (X,Y=F, Cl; n=1,5), and dimeric (H(CHB₁₁ X₁₁ ))₂ (X=F, Cl) carborane superacids performed at the B3LYP/6-311++G(d,p) theory level revealed the similarity of their equilibrium structures and the possibility of nearly barrierless hydrogen atom migration among the substituents attached to one side of the icosahedral CB₁₁ cage. The vertical electron detachment energies predicted at the OVGF/6-311++G(3df,2pd) theory level for the conjugate bases (CHB₁₁ X₁₁ )^- were found to span the 5.82-9.00 ev range. The acid strengths (manifested by the Gibbs free deprotonation energies spanning the 213-266 kcal/mol range) predicted for the icosahedral H(CHB₁₁ X₁₁ ) carborane systems confirm their superacidic properties which might be increased even further by the attachment of the second carborane H(CHB₁₁ X₁₁ ) unit that leads to a dimeric structure mimicking a part of an experimentally observed H-bridged polymeric chain. The Gibbs free deprotonation energy of the dimeric (H(CHB₁₁ Cl₁₁ ))₂ acid was predicted to be smaller by 17 kcal/mol than that of the corresponding monomeric H(CHB₁₁ Cl₁₁ ) acid.

Mechanistic Study of Enhanced Protonation by Chromium(III) in Electrospray Ionization: A Superacid Bound to a Peptide

Addition of trivalent chromium, Cr(III), to solutions undergoing electrospray ionization (ESI) enhances protonation and leads to formation of [M + 2H]²⁺ for peptides that normally produce [M + H]⁺. This effect is explored using electronic structure calculations at the density functional theory (DFT) level to predict the energetics of various species that are potentially important to the mechanism. Gas- and solution-phase reaction free energies for glycine and its anion reacting with [Cr(III)(H₂O)₆]³⁺ and for dehydration of these species have been predicted, where glycine is used as a simple model for a peptide. For comparison, calculations were also performed with Fe(III), Al(III), Sc(III), Y(III), and La(III). Removal of water from these complexes, as would occur during the ESI desolvation process, results in species that are highly acidic. The calculated pK_a of Cr(III) with a single solvation shell is -10.8, making [Cr(III)(H₂O)₆]³⁺ a superacid that is more acidic than sulfuric acid (pK_a = -8.8). Binding to glycine requires removal of two aqua ligands, which gives [Cr(III)(H₂O)₄]³⁺ that has an extremely acidic pK_a of -28.8. Removal of additional water further enhances acidity, reaching a pK_a of -84.7 for [Cr(III)(H₂O)]³⁺. A mechanism for enhanced protonation is proposed that incorporates computational and experiment results, as well as information on the known chemistry of Cr(III), which is substitutionally inert. The initial step involves binding of [Cr(III)(H₂O)₄]³⁺ to the deprotonated C-terminus of a peptide. As the drying process during ESI strips water from the complex, the resulting superacid transfers protons to the bound peptide, eventually leading to formation of [M + 2H]²⁺.

Designing aromatic heterocyclic superacids in terms of Brønsted and Lewis perspectives

The unexplored area of organic superacids was investigated in terms of both Brønsted and Lewis concepts of acids and bases. The primary requirement of a superacid-high affinity for electron/fluoride ions was fulfilled using two strategies: (i) using the superhalogen-type heterocyclic framework and (ii) selecting systems that have an electron count one short of attaining (4n + 2) Hückel aromaticity. With these in mind, eleven systems were considered throughout the study, expected to cross the target of 100% H₂SO₄ acidity and/or the fluoride affinity of SbF₅. To enhance the pK_a and F^- affinity values of the considered systems, electron-withdrawing ligands F and CN were employed. The superhalogen and aromaticity properties were verified by vertical detachment energy (VDE) and nucleus independent chemical shift (NICS) calculations, respectively. Finally, the collective effect of the potential super Lewis acids was looked into using a BL₃ skeleton with them acting as ligands.

Validation and extension of the gas-phase superacidity scale

RATIONALE

In recent years it has become increasingly evident that the previously reported experimental gas-phase acidity (GA) values of several strong acids differ markedly from the corresponding high-level computational values. In this work, the superacidic part of the current gas-phase acidity scale was validated and extended.

METHODS

For that, the strongly acidic section of the gas-phase acidity scale was remeasured using the equilibrium Fourier transform ion cyclotron resonance (FTICR-MS) method, adding new compounds and introducing methodological changes. In particular, a novel approach for anchoring the scale was used - the results were anchored to the computational (W1BD) GA values of trifluoromethanesulfonic acid and bis(fluorosulfonyl)imide (291.3 and 286.2 kcal mol-1 , respectively).

RESULTS

The newly measured section consists of 20 gas-phase superacids and its consistency standard deviation is 0.2 kcal mol-1 , indicating good consistency. In contrast to the previously reported experimental gas-phase acidities for a number of important superacids, the current results are consistent with high-level theoretical GA values. Structure-acidity relationships based on the current results as well as available MeCN and DCE acidity data were described and explained.

CONCLUSIONS

The introduced methodological innovations were found to be adequate and strong evidence is presented in support of the current GA values of the strong acids.

Quantum Chemical Modeling of the Structure and H Bonding in Triethanolammonium-Based Protic Ionic Liquids with Sulfonic Acids

The results of electronic structure calculations based on density functional theory (DFT) for protic ionic liquids (PILs) consisting of triethanolammonium cation paired with anion of different sulfonic acids are reported. The influence of the anion nature on the structure and interactions in the ion pairs that are formed in these PILs is discussed in detail. Multiple H-bonding interactions exist between the protons in the NH/OH groups of the cation and different oxygen atoms of the acid anion in the ion pairs. The quantum theory of "atoms in molecules" has been used to estimate the individual contributions of each hydrogen bond to the stability of the ion pair. The hydrogen-bonding interactions in the ion pairs vary in their strength ranging from weak to moderately strong. In addition to these hydrogen bonds, there are other dispersion and electrostatic-dominant interactions that play an important role in the overall stability of PILs and their physicochemical properties. Aided by results from our previous DFT studies of triethanolammonium class of PILs with inorganic anions, these new data allow us to gain an improved understanding of the structure-property relationships in the studied ionic liquids. Close to linear correlation, in particular, has been found between the melting points and the binding energies of the cation and anion in the ion pairs.

The Acid Strength of the Lewis-Brønsted Superacids - A QSPR Study

The acidity of Lewis-Brønsted superacids can be derived from the theoretical calculations as the Gibbs free energy of the deprotonation reaction (ΔG_acid ), which describes the tendency of a studied compound to donate a proton. This paper presents the first Quantitative Structure - Property Relationship (QSPR) model that correlates the ΔG_acid of superacid (HF/MeX₃ formula (X=F, Cl, Br)) with their structure. Developed model is well fitted, roubustness, has good predictive abilities, fulfills all OECD recommendation for good model. Obtained results provide the insight into the relation of structural features of superacids, which are responsible for their acid strength - the structures characterized by strong F-Me dative bond (with relatively large vibrational frequency), small positive partial atomic charge on Me central atom, possibly large polarity exhibit large acid strength. Such assumption can be used in the future as valuable information in the process of the designing new, stronger, more effective superacids.

Superhalogens as building blocks of a new series of superacids

A new series of superacids by the protonation of B_nH_3n+1⁻ superhalogen anions has been proposed. The resulting B_nH_3n+2 species behave as superacids for n ≥ 2 due to their smaller free energy of deprotonation than that of H₂SO₄. These B_nH_3n+2 superacids do not require Brønsted/Lewis acid components.

Structural and Theoretical Study of Salts of the [B9 H14 ]- Ion: Isolation of Multiple Isomers and Implications for Energy Storage

Boranes and boron hydrides are well known for their novel molecular structures and useful chemical reactivity, with [B₉ H₁₄ ]^- notable in particular for its ease of isolation, unusual structure, and tautomerization. We report an experimental and theoretical investigation of the structure of [B₉ H₁₄ ]^- and the energetics of some of its reactions. Salts of [B₉ H₁₄ ]^- with 1-ethyl-3-methylimidazolium and N-butyl-N-methylpyrrolidinium were characterized by single-crystal X-ray diffraction and demonstrate the stabilization of an isomer of [B₉ H₁₄ ]^- not previously observed in the solid state. Heats of formation and acid dissociation constants of [B₉ H₁₄ ]^- and closely related structures were calculated. The results suggest a mechanism for particularly energetic hypergolic ignition induced by protonation and suggest potential for reversible H₂ storage. These results encourage further investigation of [B₉ H₁₄ ]^- as an energy-storage medium in ionic systems.

The saturation of the gas phase acidity of nHF/AlF3 and nHF/GeF4 (n = 1–6) superacids caused by increasing the number of surrounding HF molecules

The acidic strength of selected Brønsted/Lewis superacids is evaluated on the basis of theoretical calculations carried out at the QCISD/6-311++G(d,p) level.

Chemistry of Fluorinated Carbon Acids: Synthesis, Physicochemical Properties, and Catalysis

The bis[(trifluoromethyl)sulfonyl]methyl (Tf2CH; Tf=SO2CF3) group is known to be one of the strongest carbon acid functionalities. The acidity of such carbon acids in the gas phase is stronger than that of sulfuric acid. Our recent investigations have demonstrated that this type of carbon acids work as novel acid catalysts. In this paper, recent achievements in carbon acid chemistry by our research group, including synthesis, physicochemical properties, and catalysis, are summarized.

An Experimental and Computational Study of the Gas-Phase Acidities of the Common Amino Acid Amides

Using proton-transfer reactions in a Fourier transform ion cyclotron resonance mass spectrometer and correlated molecular orbital theory at the G3(MP2) level, gas-phase acidities (GAs) and the associated structures for amides corresponding to the common amino acids have been determined for the first time. These values are important because amino acid amides are models for residues in peptides and proteins. For compounds whose most acidic site is the C-terminal amide nitrogen, two ions populations were observed experimentally with GAs that differ by 4-7 kcal/mol. The lower energy, more acidic structure accounts for the majority of the ions formed by electrospray ionization. G3(MP2) calculations predict that the lowest energy anionic conformer has a cis-like orientation of the [-C(═O)NH](-) group whereas the higher energy, less acidic conformer has a trans-like orientation of this group. These two distinct conformers were predicted for compounds with aliphatic, amide, basic, hydroxyl, and thioether side chains. For the most acidic amino acid amides (tyrosine, cysteine, tryptophan, histidine, aspartic acid, and glutamic acid amides) only one conformer was observed experimentally, and its experimental GA correlates with the theoretical GA related to side chain deprotonation.

Nanoscale metal fluorides: a new class of heterogeneous catalysts

Nanoscale metal fluorides and hydroxide fluorides prepared according the fluorolytic sol–gel synthesis represent a powerful class of bi-acidic heterogeneous catalysts.

Development of Monopole Interaction Models for Ionic Compounds. Part I: Estimation of Aqueous Henry's Law Constants for Ions and Gas Phase pKa Values for Acidic Compounds

The SPARC (SPARC Performs Automated Reasoning in Chemistry) physicochemical mechanistic models for neutral compounds have been extended to estimate Henry's Law Constant (HLC) for charged species by incorporating ionic electrostatic interaction models. Combinations of absolute aqueous pKa values, relative pKa values in the gas phase, and aqueous HLC for neutral compounds have been used to develop monopole interaction models that quantify the energy differences upon moving an ionic solute molecule from the gas phase to the liquid phase. Inter-molecular interaction energies were factored into mechanistic contributions of monopoles with polarizability, dipole, H-bonding, and resonance. The monopole ionic models were validated by a wide range of measured gas phase pKa data for 450 acidic compounds. The RMS deviation error and R(2) for the OH, SH, CO2 H, CH3 and NR2 acidic reaction centers (C) were 16.9 kcal/mol and 0.87, respectively. The calculated HLCs of ions were compared to the HLCs of 142 ions calculated by quantum mechanics. Effects of inter-molecular interaction of the monopoles with polarizability, dipole, H-bonding, and resonance on acidity of the solutes in the gas phase are discussed.

Dehydration, dehydrogenation, and condensation of alcohols on supported oxide catalysts based on cyclic (WO3)3 and (MoO3)3 clusters

The review summarizes recent synthesis and reactivity studies of model oxide catalysts prepared by the deposition of gas phase cyclic (WO₃)₃ and (MoO₃)₃ clusters.