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.

In search of CS2(H2O)n = 1–4 clusters

Gaussian-3 and MP2/aug-cc-pVnZ methods have been used to calculate geometries and thermochemistry of CS2(H2O)n, where n = 1–4. An extensive molecular dynamics search followed by optimization using these two methods located two dimers, six trimers, six tetramers, and two pentamers. The MP2/aug-cc-pVDZ structure matched best with the experimental result for the CS2(H2O) dimer, showing that diffuse functions are necessary to model the interactions found in this complex. For larger CS2(H2O)n clusters, the MP2/aug-cc-pVDZ minima are significantly different from the MP2(full)/6-31G* structures, revealing that the G3 model chemistry is not suitable for investigation of sulfur containing van der Waals complexes. Based on the MP2/aug-cc-pVTZ free energies, the concentration of saturated water in the atmosphere and the average amount of CS2 in the atmosphere, the concentrations of these clusters are predicted to be on the order of 105CS2(H2O) clusters∙cm−3 and 102 CS2(H2O)2 clusters∙cm−3 at 298.15 K. The MP2/aug-cc-pVDZ scaled harmonic and anharmonic frequencies of the most abundant dimer cluster at 298 K are presented, along with the MP2/aug-cc-pVDZ scaled harmonic frequencies for the CS2(H2O)n structures predicted to be present in a low-temperature molecular beam experiment.

Hydration of OCS with One to Four Water Molecules in Atmospheric and Laboratory Conditions

Carbonyl sulfide is the most abundant sulfur gas in the atmosphere. We have used MP2 and CCSD(T) theory to study the structures and thermochemistries of carbonyl sulfide interacting with one to four water molecules. We have completed an extensive search for clusters of OCS(H2O)n, where n = 1−4. We located three dimers, two trimers, five tetramers, and four pentamers with the MP2/aug-cc-pVDZ method. In each of the complexes with two or more waters, OCS preferentially interacts with low-energy water clusters. Our results match current theoretical and experimental literature, showing correlation with available geometries and frequencies for the OCS(H2O) species. The CCSD(T)/aug-cc-pVTZ thermochemical values combined with the average amount of OCS and the saturated concentration of H2O in the troposphere, lead to the prediction of 106 OCS(H2O) clusters·cm−3 and 102 OCS(H2O)2 clusters·cm−3 at 298 K. We predict the structures of OCS(H2O)n, n = 1−4 that should predominate in a low-temperature molecular beam and identify specific infrared vibrations that can be used to identify these different clusters.

Synthesis of Mid-Chain Functionalized Polymers

Polymers with mid-chain alkoxyamine functionality were synthesized by activating monohalogenated polymers in the presence of nitroso or nitrone radical traps. The resulting polymers were either polystyrene (PSt) homopolymers with a mid-chain alkoxyamine or PSt-poly(methyl acrylate) (PMA) diblock copolymers with an alkoxyamine unit at the junction between the segments. Monohalogenated polymers where synthesized by atom transfer radical polymerization (ATRP) and were then reacted to form polymer radicals in the presence of a radical trap, nitrone or nitroso. When only polystyrene radicals were reacted with the radical trap a dimer was formed with an alkoxyamine functionality in the center of the polymer chain. This functionality allowed the polymer chain to be cleaved in order to visualize the extent of the alkoxyamine functionality incorporation into the polymer chains. It was found that near quantitative alkoxyamine mid-chain functionality could be achieved by activating the PStBr in the presence of 10 equivalents of nitrone, 5 equivalents of copper bromide, and 2 equivalents of copper metal. Further reducing the amount of copper metal led to incomplete coupling, while increasing the equivalents beyond 2 generated polymer dimers with less than quantitative mid-chain functionality. Monochlorinated polystyrene (PStCl) precursors gave much poorer coupling results compared to reactions with PStBr, which is consistent with the stronger C-Cl bond resisting activation and the formation of the polystyryl radicals. When poly (methyl acrylate) (PMABr) is reacted with PStBr in the presence of a nitroso group at reduced temperatures (30 oC) block copolymers were selectively formed with an alkoxyamine functionality in the center. This was done by first activating the PSt-Br to form a polymer radical that would react with the radical trap to form a persistent radical on the oxygen. The PMA-Br, once activated, reacted with the radical on the oxygen to form the block copolymer. To test the amount of functionality incorporated, a coupling reaction was performed with no nitroso present, and found that no reaction occurred. This showed that the radical trap is essential for the coupling to occur, and cleavage of the diblock indicated that the alkoxyamine functionality was indeed incorporated into the diblock.