Determinant Factor for Thermodynamic Stability of Sulfuric Acid-Amine Complexes

Proton transfer between sulfonic acids and various propylamines by density functional theory calculations

CONTEXT

Proton transfer in acid-base systems is not well understood. Some acid-base reactions do not proceed to the extent that is expected from the difference in the pK_a values between the base and acid in aqueous solutions, yet some do. In that regard, we have computationally studied the process of proton transfer from the acids of varying strength (benzenesulfonic acid (BSu), methansulfonic acid (MsO), and sulfuric acid (SA)) to the amines with different numbers of propyl substituents on the nitrogen atom (propylamine (PrA), dipropylamine (DPrA), and tripropylamine (TPrA)) upon complexation. Density functional theory calculations were used to thoroughly examine the energetic and structural aspects of the molecular complexes and/or ionic pairs resulting from the acid-base interaction. The potential energy curves along the proton transfer coordinate in these acid-amine systems were analyzed. The change in free energies accompanying the molecular complexes and ionic pair formations was calculated, and the relationship between the energy values and the ΔРА parameter (difference in proton affinity of the acid anion and amine) was established. The larger ΔРА values were found to be unfavorable for the formation of ionic pairs. Using structural, energy, QTAIM, and NBO analyses, we determined that the hydrogen bonds in the molecular complexes PrA-MsO and PrA-BSu are stronger than those in their corresponding ionic pairs. The ionic pairs with the TPrA cation possess the strongest hydrogen bonds of all the ionic pairs being studied, regardless of the anion. The results showed that hydrogen bonding interactions in the molecular complexes contribute significantly to the energies of the acid-base interaction, while in the ionic pairs, the most important energy contribution comes from Coulomb interactions, followed by hydrogen bonding and dispersion forces. The ionic pairs with propylammonium, dipropylammonium, and tripropylammonium cations have stronger ion-ion interactions than tetrapropylammonium (TetPrA)-containing ionic pairs with the same anions. This effect rises with the order of the cation: TetPrA → TPrA → DPrA → PrA, and the sequence of anions is SA → BSu → MsO. The results obtained here expand the concept of acid-base interaction and provide an alternative to experimental searches for suitable acids and bases to obtain new types of protic ionic liquids.

METHODS

All quantum-chemical calculations were carried out by using the DFT/B3LYP-GD3/6-31++G(d,p) level as implemented in the Gaussian 09 software package. For the resulting structures, the electron density distribution was analyzed by the "atoms in molecules" (QTAIM) and the natural bond orbital (NBO) methods on the wave functions obtained at the same level of theory by AIMAll Version 10.05.04 and Gaussian NBO Version 3.1 programs, respectively.

Amino Acids Compete with Ammonia in Sulfuric Acid-Based Atmospheric Aerosol Prenucleation: The Case of Glycine and Serine

We present a computational investigation of the sulfuric acid, glycine, serine, ammonia, and water system to understand if this system can form prenucleation clusters, which are precursors to larger aerosols in the atmosphere. We have performed a comprehensive configurational search of all possible clusters in this system, starting with the four different monomers and zero to five waters. Accurate Gibbs free energies of formation have been calculated with the DLPNO-CCSD(T)/complete basis set (CBS) method on ωb97xd/6-31++G** geometries. For the dry dimers of sulfuric acid, the weakest base, serine, is found to form the most stable complex, which is a consequence of the strong di-ionic complex formed between the bisulfate ion and the protonated serine cation. For the dry dimers without sulfuric acid, the glycine-serine complex is more stable than the glycine-ammonia or serine-ammonia complexes, stemming from the detailed structure and not related to base strength. For the larger complexes, sulfuric acid deprotonates and the proton is shifted to glycine, serine, or ammonia. The two amino acids and ammonia are almost interchangeable and there is no easy way to predict which molecule will be protonated without the calculated results. Assuming reasonable starting concentrations and a closed system of sulfuric acid, glycine, serine, ammonia, and five waters, we predict the concentrations of all possible complexes at two temperatures spanning the troposphere. The most negative ΔG° values are a function of the detailed molecular interactions of these clusters. These details are more important than the base strength of ammonia, glycine, and serine.