Computational Study of the Hydration of Sulfuric Acid Dimers: Implications for Acid Dissociation and Aerosol Formation

We have investigated the thermodynamics of sulfuric acid dimer hydration using ab initio quantum mechanical methods. For (H2SO4)2(H2O)n where n = 0−6, we employed high-level ab initio calculations to locate the most stable minima for each cluster size. The results presented herein yield a detailed understanding of the first deprotonation of sulfuric acid as a function of temperature for a system consisting of two sulfuric acid molecules and up to six waters. At 0 K, a cluster of two sulfuric acid molecules and one water remains undissociated. Addition of a second water begins the deprotonation of the first sulfuric acid leading to the di-ionic species (the bisulfate anion HSO4−, the hydronium cation H3O+, an undissociated sulfuric acid molecule, and a water). Upon the addition of a third water molecule, the second sulfuric acid molecule begins to dissociate. For the (H2SO4)2(H2O)3 cluster, the di-ionic cluster is a few kcal mol−1 more stable than the neutral cluster, which is just slightly more stable than the tetra-ionic cluster (two bisulfate anions, two hydronium cations, and one water). With four water molecules, the tetra-ionic cluster, (HSO4−)2(H3O+)2(H2O)2, becomes as favorable as the di-ionic cluster H2SO4(HSO4−)(H3O+)(H2O)3 at 0 K. Increasing the temperature favors the undissociated clusters, and at room temperature we predict that the di-ionic species is slightly more favorable than the neutral cluster once three waters have been added to the cluster. The tetra-ionic species competes with the di-ionic species once five waters have been added to the cluster. The thermodynamics of stepwise hydration of sulfuric acid dimer is similar to that of the monomer; it is favorable up to n = 4−5 at 298 K. A much more thermodynamically favorable pathway forming sulfuric acid dimer hydrates is through the combination of sulfuric acid monomer hydrates, but the low concentration of sulfuric acid relative to water vapor at ambient conditions limits that process

Hydration of the Sulfuric Acid−Methylamine Complex and Implications for Aerosol Formation

The binary H2SO4−H2O nucleation is one of the most important pathways by which aerosols form in the atmosphere, and the presence of ternary species like amines increases aerosol formation rates. In this study, we focus on the hydration of a ternary system of sulfuric acid (H2SO4), methylamine (NH2CH3), and up to six waters to evaluate its implications for aerosol formation. By combining molecular dynamics (MD) sampling with high-level ab initio calculations, we determine the thermodynamics of forming H2SO4(NH2CH3)(H2O)n, where n = 0−6. Because it is a strong acid−base system, H2SO4−NH2CH3 quickly forms a tightly bound HSO4−−NH3CH3+ complex that condenses water more readily than H2SO4 alone. The electronic binding energy of H2SO4−NH2CH3 is −21.8 kcal mol−1 compared with −16.8 kcal mol−1 for H2SO4−NH3 and −12.8 kcal mol−1 for H2SO4−H2O. Adding one to two water molecules to the H2SO4−NH2CH3 complex is more favorable than adding to H2SO4 alone, yet there is no systematic difference for n ≥ 3. However, the average number of water molecules around H2SO4−NH2CH3 is consistently higher than that of H2SO4, and it is fairly independent of temperature and relative humidity

Computation of Atmospheric Concentrations of Molecular Clusters from \u3cem\u3eab initio\u3c/em\u3e Thermochemistry

The computational study of the formation and growth of atmospheric aerosols requires an accurate Gibbs free energy surface, which can be obtained from gas phase electronic structure and vibrational frequency calculations. These quantities are valid for those atmospheric clusters whose geometries correspond to a minimum on their potential energy surfaces. The Gibbs free energy of the minimum energy structure can be used to predict atmospheric concentrations of the cluster under a variety of conditions such as temperature and pressure. We present a computationally inexpensive procedure built on a genetic algorithm-based configurational sampling followed by a series of increasingly accurate screening calculations. The procedure starts by generating and evolving the geometries of a large set of configurations using semi-empirical models then refines the resulting unique structures at a series of high-level ab initio levels of theory. Finally, thermodynamic corrections are computed for the resulting set of minimum-energy structures and used to compute the Gibbs free energies of formation, equilibrium constants, and atmospheric concentrations. We present the application of this procedure to the study of hydrated glycine clusters under ambient conditions

Computation of Atmospheric Concentrations of Molecular Clusters from ab initio Thermochemistry

The computational study of the formation and growth of atmospheric aerosols requires an accurate Gibbs free energy surface, which can be obtained from gas phase electronic structure and vibrational frequency calculations. These quantities are valid for those atmospheric clusters whose geometries correspond to a minimum on their potential energy surfaces. The Gibbs free energy of the minimum energy structure can be used to predict atmospheric concentrations of the cluster under a variety of conditions such as temperature and pressure. We present a computationally inexpensive procedure built on a genetic algorithm-based configurational sampling followed by a series of increasingly accurate screening calculations. The procedure starts by generating and evolving the geometries of a large set of configurations using semi-empirical models then refines the resulting unique structures at a series of high-level ab initio levels of theory. Finally, thermodynamic corrections are computed for the resulting set of minimum-energy structures and used to compute the Gibbs free energies of formation, equilibrium constants, and atmospheric concentrations. We present the application of this procedure to the study of hydrated glycine clusters under ambient conditions

Broadband Fourier transform rotational spectroscopy for structure determination: The water heptamer

Over the recent years chirped-pulse, Fourier-transform microwave (CP-FTMW) spectrometers have chan- ged the scope of rotational spectroscopy. The broad frequency and large dynamic range make possible structural determinations in molecular systems of increasingly larger size from measurements of heavy atom (13C, 15N, 18O) isotopes recorded in natural abundance in the same spectrum as that of the parent isotopic species. The design of a broadband spectrometer operating in the 2–8 GHz frequency range with further improvements in sensitivity is presented. The current CP-FTMW spectrometer performance is benchmarked in the analyses of the rotational spectrum of the water heptamer, (H2O)7, in both 2– 8 GHz and 6–18 GHz frequency ranges. Two isomers of the water heptamer have been observed in a pulsed supersonic molecular expansion. High level ab initio structural searches were performed to pro- vide plausible low-energy candidates which were directly compared with accurate structures provided from broadband rotational spectra. The full substitution structure of the most stable species has been obtained through the analysis of all possible singly-substituted isotopologues (H218O and HDO), and a least-squares rm(1) geometry of the oxygen framework determined from 16 different isotopic species compares with the calculated O–O equilibrium distances at the 0.01 Å level

The Role of Non-Covalent Interactions on Cluster Formation: Pentamer, Hexamers and Heptamer of Difluoromethane

The role of non-covalent interactions (NCIs) has broadened with the inclusion of new types of interactions and a plethora of weak donor/acceptor partners. This work illustrates the potential of chirped-pulse Fourier transform microwave technique, which has revolutionized the field of rotational spectroscopy. In particular, it has been exploited to reveal the role of NCIs' in the molecular self-aggregation of difluoromethane where a pentamer, two hexamers and a heptamer were detected. The development of a new automated assignment program and a sophisticated computational screening protocol was essential for identifying the homoclusters in conditions of spectral congestion. The major role of dispersion forces leads to less directional interactions and more distorted structures than those found in polar clusters, although a detailed analysis demonstrates that the dominant interaction energy is the pairwise interaction. The tetramer cluster is identified as a structural unit in larger clusters, representing the maximum expression of bond between dimers.We thank MINECO (CTQ2017-89150-R), Basque Government (IT1162-19 and PIBA2018-11), the UPV/EHU (PPG17/10, GIU18/207), CSIC (2018FR0036, LINKA20249), University of Bologna (RFO), Fondazione CARISBO (2018/0353) and NSF (CHE-1903871 and CHE-2018427) for the financial support. C.C thanks MINECO for a Juan de la Cierva contract. L.E. was supported by Marie Curie fellowship PIOF-GA-2012-32840

Computation of Molecular Properties at the Ab Initio Limit

The accuracy of a quantum chemical calculation inherently depends on the ability to account for the completeness of the one- and n-particle spaces. The size of the basis set used can be systematically increased until it reaches the complete one-particle basis set limit (CBS) while the n-particle space approaches its exact full configuration interaction (FCI) limit by following a hierarchy of electron correlation methods developed over the last seventy years. If extremely high accuracy is desired, properly correcting for very small effects such as those resulting the Born-Oppenheimer approximation and the neglect of relativistic effects becomes indispensable. For a series of chemically interesting and challenging systems, we identify the limits of conventional approaches and use state-of-the-art quantum chemical methods along with large basis sets to get the “right answer for the right reasons.” First, we quantify the importance of small effects that are ignored in conventional quantum chemical calculations and manage to achieve spectroscopic accuracy (agreement of 1 cm−1 or less with experimental harmonic vibrational frequencies) for BH, CH+ and NH. We then definitively resolve the global minimum structure for Li₆ , Li₆⁺ , and Li₆- using high accuracy calculations of the binding energies, ionization potentials, electron affinities and vertical excitation spectra for the competing isomers. The same rigorous approach is used to study a series of hydrogen transfer reactions and validate the necessary parameters for the hydrogen abstraction and donation steps in the mechanosynthesis of diamondoids. Finally, in an effort to overcome the steep computational scaling of most high-level methods, a new hybrid methodology which scales as O(N⁵) but performs comparably to O(N⁶) methods is benchmarked for its performance in the equilibrium and dissociation regimes.Ph.D.Committee Chair: Sherrill, C. David; Committee Member: Brédas, Jean-Luc; Committee Member: Hernandez, Rigoberto; Committee Member: Merkle, Ralph C.; Committee Member: Whetten, Robert L

Hybrid correlation models based on active-space partitioning: Seeking accurate O(N ⁵) ab initio methods for bond breaking

© 2006 American Institute of Physics. The electronic version of this article is the complete one and can be found at: http://dx.doi.org/10.1063/1.2222350DOI: 10.1063/1.2222350Møller-Plesset second-order (MP2) perturbation theory remains the least expensive standard ab initio method that includes electron correlation, scaling as O(N ⁵) with the number of molecular orbitals N. Unfortunately, when restricted Hartree-Fock orbitals are employed, the potential energy curves calculated with this method are of little use at large interatomic separations because of the divergent behavior of MP2 in these regions. In our previous study [J. Chem. Phys. 122, 234110 (2005)] we combined the MP2 method with the singles and doubles coupled cluster (CCSD) method to produce a hybrid method that retains the computational scaling of MP2 and improves dramatically the shape of the MP2 curves. In this work we expand the hybrid methodology to several other schemes. We investigate a new, improved MP2-CCSD method as well as a few other O(N ⁵) methods related to the Epstein-Nesbet pair correlation theory. Nonparallelity errors across the dissociation curve as well as several spectroscopic constants are computed for BH, HF, H₂O, CH+, CH₄, and Li₂ molecules with the 6-31G* basis set and compared with the corresponding full configuration interaction results. We show that among the O(N ⁵) methods considered, our new hybrid MP2-CCSD method is the most accurate and significantly outperforms MP2 not only at large interatomic separations, but also near equilibrium geometries