Formation and growth of molecular clusters containing sulfuric acid, water, ammonia, and dimethylamine

The structures and thermochemistry of molecular clusters containing sulfuric acid, water, ammonia, and/or dimethylamine ((CH3)2NH or DMA) are explored using a combination of Monte Carlo configuration sampling, semiempirical calculations, and density functional theory (DFT) calculations. Clusters are of the general form [(BH(+))n(HSO4(-))n(H2O)y], where B = NH3 or DMA, 2 ≤ n ≤ 8, and 0 ≤ y ≤ 10. Cluster formulas are written based on the computed structures, which uniformly show proton transfer from each sulfuric acid molecule to a base molecule while the water molecules remain un-ionized. Cluster formation is energetically favorable, owing to strong electrostatic attraction among the ions. Water has a minor effect on the energetics of cluster formation, lowering the free energy of formation by ∼ 10% depending on the cluster size and number of water molecules. Cluster growth (addition of one base molecule and one sulfuric acid molecule to a pre-existing cluster) and base substitution (substituting DMA for ammonia) are also energetically favorable processes for both anhydrous and hydrated clusters. However, the effect of water is different for different bases. Hydrated ammonium bisulfate clusters have a more favorable free energy for growth (i.e., incrementing n with fixed y) than anhydrous clusters, while the reverse is observed for dimethylammonium bisulfate clusters, where the free energy for growth is more favorable for anhydrous clusters. The substitution of DMA for ammonia in bisulfate clusters is favorable but exhibits a complex water dependence. Base substitution in smaller bisulfate clusters is enhanced by the presence of water, while base substitution in larger bisulfate clusters is less favorable for hydrated clusters than that for anhydrous clusters. While DMA substitution can stabilize small clusters containing one or a few sulfuric acid molecules, the free energy advantage of forming amine clusters relative to ammonia clusters becomes less pronounced at larger sizes, especially when the effect of water is considered.