Quantum Mechanical Study of Sulfuric Acid Hydration: Atmospheric Implications

Energetics and spectroscopic studies of CNO (-) (H 2 O) n clusters and the temperature dependencies of the isomers: An approach based on a combined recipe of parallel tempering and quantum chemical methods

A system associated with several number of weak interactions supports numerous number of stable structures within a narrow range of energy. Often, a deterministic search method fails to locate the global minimum geometry as well as important local minimum isomers for such systems. Therefore, in this work, the stochastic search technique, namely parallel tempering, has been executed on the quantum chemical surface of theCNO (-) (H 2 O) n system for n = 1 -8 to generate global minimum as well as several number of local minimum isomers. IR spectrum can act as the fingerprint property for such system to be identified. Thus, IR spectroscopic features have also been included in this work. Vertical detachment energy has also been calculated to obtain clear information about number of water molecules in several spheres around the central anion. In addition, in a real experimental scenario, not only the global but also the local minimum isomers play an important role in determining the average value of a particular physically observable property. Therefore, the relative conformational populations have been determined for all the evaluated structures for the temperature range between 20K and 400K. Further to understand the phase change behavior, the configurational heat capacities have also been calculated for different sizes.

Prebiotic dimer and trimer peptide formation in gas-phase atmospheric nanoclusters of water

Insight into the origin of prebiotic molecules is key to our understanding of how living systems evolved into the complex network of biological processes on Earth. By modelling diglycine and triglycine peptide formation in the prebiotic atmosphere, we provide a plausible pathway for peptide growth. By examining different transition states (TSs), we conclude that the formation of diglycine and triglycine in atmospheric nanoclusters of water in the prebiotic atmosphere kinetically favors peptide growth by an N-to-C synthesis of glycines through a trans conformation. Addition of water stabilizes the TS structures and lowers the Gibbs free activation energies. At temperatures that model the prebiotic atmosphere, the free energies of activation with a six water nanocluster as part of the TS are predicted to be 16 kcal mol-1 relative to the prereactive complex. Examination of the trans vs. cis six water transition states reveals that a homodromic water network that maximizes the acceptor/donor nature of the six waters is responsible for enhanced kinetic favorability of the trans N-to-C pathway. Compared to the non-hydrated trans TS, the trans six-water TS accelerates the reaction of diglycine and glycine to form triglycine by 13 orders of magnitude at 217 K. Nature uses the trans N-to-C pathway to synthesize proteins in the ribosome, and we note the similarities in hydrogen bond stabilization between the transition state for peptide synthesis in the ribosome and the transition states formed in nanoclusters of water in the same pathway. These results support the hypothesis that small oligomers formed in the prebiotic atmosphere and rained onto earth's surface.

Theoretical Study of the Gas-Phase Hydrolysis of Formaldehyde to Produce Methanediol and Its Implication to New Particle Formation

Aldehydes were speculated to be important precursor species in new particle formation (NPF). The direct involvement of formaldehyde (CH₂O) in sulfuric acid and water nucleation is negligible; however, whether its atmospheric hydrolysate, methanediol (CH₂(OH)₂), which contains two hydroxyl groups, participates in NPF is not known. This work investigates both CH₂O hydrolysis and NPF from sulfuric acid and CH₂(OH)₂ with quantum chemistry calculations and atmospheric cluster dynamics modeling. Kinetic calculation shows that reaction rates of the gas-phase hydrolysis of CH₂O catalyzed by sulfuric acid are 11-15 orders of magnitude faster than those of the naked path at 253-298 K. Based on structures and the calculated formation Gibbs free energies, the interaction between sulfuric acid/its dimer/its trimer and CH₂(OH)₂ is thermodynamically favorable, and CH₂(OH)₂ forms hydrogen bonds with sulfuric acid/its dimer/its trimer via two hydroxyl groups to stabilize clusters. Our further cluster kinetic calculations suggested that the particle formation rates of the system are higher than those of the binary system of sulfuric acid and water at ambient low sulfuric acid concentrations and low relative humidity. In addition, the formation rate is found to present a negative temperature dependence because evaporation rate constants contribute significantly to it. However, cluster growth is essentially limited by the weak formation of the largest clusters, which implies that other stabilizing vapors are required for stable cluster formation and growth.

Massive Assessment of the Binding Energies of Atmospheric Molecular Clusters

Quantum chemical studies of the formation and growth of atmospheric molecular clusters are important for understanding aerosol particle formation. However, the search for the lowest free-energy cluster configuration is extremely time consuming. This makes high-level benchmark data sets extremely valuable in the quest for the global minimum as it allows the identification of cost-efficient computational methodologies, as well as the development of high-level machine learning (ML) models. Herein, we present a highly versatile quantum chemical data set comprising a total of 11 749 (acid)_1-2(base)_1-2 cluster configurations, containing up to 44 atoms. Utilizing the LUMI supercomputer, we calculated highly accurate PNO-CCSD(F12*)(T)/cc-pVDZ-F12 binding energies of the full set of cluster configurations leading to an unprecedented data set both in regard to sheer size and with respect to the level of theory. We employ the constructed benchmark set to assess the performance of various semiempirical and density functional theory methods. In particular, we find that the r²-SCAN-3c method shows excellent performance across the data set related to both accuracy and CPU time, making it a promising method to employ during cluster configurational sampling. Furthermore, applying the data sets, we construct ML models based on Δ-learning and provide recommendations for future application of ML in cluster configurational sampling.

Computational chemistry of cluster: Understanding the mechanism of atmospheric new particle formation at the molecular level

New particle formation (NPF), which exerts significant influence over human health and global climate, has been a hot topic and rapidly expands field of research in the environmental and atmospheric chemistry recent years. Generally, NPF contains two processes: formation of critical nucleus and further growth of the nucleus. However, due to the complexity of the atmospheric nucleation, which is a multicomponent process, formation of critical clusters as well as their growth is still connected to large uncertainties. Detection limits of instruments in measuring specific gaseous aerosol precursors and chemical compositions at the molecular level call for computational studies. Computational chemistry could effectively compensate the deficiency of laboratory experiments as well as observations and predict the nucleation mechanisms. We review the present theoretical literatures that discuss nucleation mechanism of atmospheric clusters. Focus of this review is on different nucleation systems involving sulfur-containing species, nitrogen-containing species and iodine-containing species. We hope this review will provide a deep insight for the molecular interaction of nucleation precursors and reveal nucleation mechanism at the molecular level.

The driving effects of common atmospheric molecules for formation of prenucleation clusters: the case of sulfuric acid, formic acid, nitric acid, ammonia, and dimethyl amine

How secondary aerosols form is critical as aerosols' impact on Earth's climate is one of the main sources of uncertainty for understanding global warming. The beginning stages for formation of prenucleation complexes, that lead to larger aerosols, are difficult to decipher experimentally. We present a computational chemistry study of the interactions between three different acid molecules and two different bases. By combining a comprehensive search routine covering many thousands of configurations at the semiempirical level with high level quantum chemical calculations of approximately 1000 clusters for every possible combination of clusters containing a sulfuric acid molecule, a formic acid molecule, a nitric acid molecule, an ammonia molecule, a dimethylamine molecule, and 0-5 water molecules, we have completed an exhaustive search of the DLPNO-CCSD(T)/CBS//ωB97X-D/6-31++G** Gibbs free energy surface for this system. We find that the detailed geometries of each minimum free energy cluster are often more important than traditional acid or base strength. Addition of a water molecule to a dry cluster can enhance stabilization, and we find that the (SA)(NA)(A)(DMA)(W) cluster has special stability. Equilibrium calculations of SA, FA, NA, A, DMA, and water using our quantum chemical ΔG° values for cluster formation and realistic estimates of the concentrations of these monomers in the atmosphere reveals that nitric acid can drive early stages of particle formation just as efficiently as sulfuric acid. Our results lead us to believe that particle formation in the atmosphere results from the combination of many different molecules that are able to form highly stable complexes with acid molecules such as SA, NA, and FA.

Hydration effect of selected atmospheric gases with finite water clusters: A quantum chemical investigation towards atmospheric implications

Water vapor in atmosphere is ubiquitous, and it varies according to geographical locations. Various toxic and non-toxic gases co-exist with water vapor/moisture in the atmosphere. This computational study addresses the fact that how those gases interact with water vapor. We have done quantum chemical density functional theory calculations to probe the interaction of certain gases with a finite number of water molecules in gas phase with various functionals/basis sets. An ensemble of 14 gas molecules comprising various diatomic, triatomic, and polyatomic gases have been chosen for the investigations. The intermolecular interactions are understood from the interaction energy, electrostatic potential, frontier molecular orbitals, energy gap, and natural bond orbital analyses. Furthermore, quantum molecular descriptors such as electronegativity, chemical potential, chemical hardness and electrophilicity index are calculated to have deep insight on chemical nature of the gas molecules. Additionally, we have done implicit solvent modelling using PCM, and the corresponding solvation energies have been calculated. Interestingly, all the calculations and analyses have projected the similar results that Cl₂, SO_2, and NH₃ have very high interaction with the water clusters. To mimic various altitudes (0 km, 5 km and 10 km) in the atmosphere, thermochemistry calculations have been carried out at different temperature and pressure values. The Gibbs free energies of formation suggest that the hydration of Cl₂ is higher followed by O₂, SO₂ and NH₃ at all altitudes. Remarkably, it is found that the formation of hydrated clusters of Cl₂ and O₂ with 4H₂O are thermodynamically favourable. On the other hand, SO₂ and NH₃ requires 5H₂O and 3H₂O to form thermodynamically favourable clusters. In summary, it is anticipated that this kind of extensive computational studies facilitate to understand the structural, electronic, chemical and thermochemical properties of hydrated atmospheric gases that leads to the formation of prenucleation clusters followed by atmospheric aerosols.

Interaction of Glutamic Acid/Protonated Glutamic Acid with Amide and Water Molecules: A Theoretical Study

Amino acids are important nitrogen-containing compounds and organic carbon components that exist widely in the atmosphere. The formation of atmospheric aerosols is affected by their interactions with amides. The dimers formed by glutamic acid (Glu) or protonated glutamic acid (Glu⁺) with three kinds of amide molecules (formamide FA, acetamide AA, urea U) and the hydrated clusters formed by Glu or Glu⁺, U molecules along with one to six water molecules were systematically studied at the M06-2X/6-311++G(3df,3pd) level. U is predicted to form a more stable structure with Glu/Glu⁺ than FA and AA by thermodynamics. If the concentration ratio of FA to U is less than 10⁴, U will play a critical role in NPF. The degree of hydration in Glu⁺-mU-nW is higher than that of Glu-mU-nW (m = 0, 1; n = 0-6) clusters. Notably, Glu contributes more to the Rayleigh scattering properties than glutaric acid and sulfuric acid, and thus may lead to the destruction of atmospheric visibility. This study is helpful to better understand the properties of organic aerosols containing amino acids or amides.

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.

Hydrogen-Bond Topology Is More Important Than Acid/Base Strength in Atmospheric Prenucleation Clusters

We explored the hypothesis that on the nanoscale level, acids and bases might exhibit different behavior than in bulk solution. Our study system consisted of sulfuric acid, formic acid, ammonia, and water. We calculated highly accurate Domain-based Local pair-Natural Orbital- Coupled-Cluster/Complete Basis Set (DLPNO-CCSD(T)/CBS) energies on DFT geometries and used the resulting Gibbs free energies for cluster formation to compute the overall equilibrium constants for every possible cluster. The equilibrium constants combined with the initial monomer concentrations were used to predict the formation of clusters at the top and the bottom of the troposphere. Our results show that formic acid is as effective as ammonia at forming clusters with sulfuric acid and water. The structure of formic acid is uniquely suited to form hydrogen bonds with sulfuric acid. Additionally, it can partner with water to form bridges from one side of sulfuric acid to the other, hence demonstrating that hydrogen bonding topology is more important than acid/base strength in these atmospheric prenucleation clusters.

Vacuum ultraviolet photochemistry of sulfuric acid vapor: a combined experimental and theoretical study

We present a vacuum ultraviolet (VUV) photoionization study of the gas-phase sulfuric acid (H₂SO₄) molecule in the 11-14 eV energy range by using the method of synchrotron radiation-based double imaging photoelectron photoion coincidence (i²PEPICO) spectroscopy complemented with accurate theoretical calculations. The slow photoelectron spectrum (SPES) of H₂SO₄ has been acquired and the three electronic states of H₂SO₄⁺, X²A, A²A and B²A have been populated and assigned. The adiabatic ionization energy of the H₂SO₄ molecule towards the X²A cationic ground state is measured at 11.684 ± 0.006 eV, in accordance with high-level calculated findings. With increasing photon energies, the H₂SO₄⁺ cation dissociates into HSO₃⁺ and OH fragments and their adiabatic appearance energy is measured at 13.498 ± 0.007 eV. Then, the enthalpies of formation for the species involved in the photoionization and dissociative photoionization have been determined through a thermochemical cycle.

Monomers of Glycine and Serine Have a Limited Ability to Hydrate in the Atmosphere

The role of atmospheric aerosols on climate change is one of the biggest uncertainties in most global climate models. Organic aerosols have been identified as potential cloud condensation nuclei (CCN), and amino acids are organic molecules that could serve as CCN. Amino acids make up a significant portion of the total organic material in the atmosphere, and herein we present a systematic study of hydration for two of the most common atmospheric amino acids, glycine and serine. We compute DLPNO/CCSD(T)//M08-HX/MG3S thermodynamic properties and atmospheric concentrations of Gly(H₂O)_n and Ser(H₂O)_n, where n = 1-5. We predict that serine-water clusters have higher concentrations at n = 1 and 5, while glycine-water clusters have higher concentrations at n = 2-4. However, both glycine and serine are inferred to exist primarily in their nonhydrated monomer forms in the absence of other species such as sulfuric acid.

Microwave Study of Triflic Acid Hydrates: Evidence for the Transition from Hydrogen-Bonded Clusters to a Microsolvated Ion Pair

Rotational spectra of the mono-, di-, and trihydrates of triflic acid, CF₃SO₃H···(H₂O)_n=1-3, have been recorded by pulsed nozzle Fourier transform microwave spectroscopy and spectroscopic constants obtained have been compared with values calculated at several levels of theory. The experimental results are consistent with the theoretical predictions presented here and elsewhere, indicating that with only one or two water molecules, triflic acid remains un-ionized in a cold molecular complex. The experiments further concur with theoretical predictions that the addition of a third water molecule transforms the system into what is best regarded as a hydrated hydronium triflate ion pair. Thus, only three water molecules are needed to induce ionization of triflic acid in a cold molecular cluster. This number is somewhat low compared with that for other simple protic acids and likely reflects the superacidity of triflic acid. Simple energetic arguments can be used to rationalize this result.

Calculating Reliable Gibbs Free Energies for Formation of Gas-Phase Clusters that Are Critical for Atmospheric Chemistry: (H2SO4)3

The effects of atmospheric aerosols on our climate are one of the biggest uncertainties in global climate models. Calculating the pathway for the formation of pre-nucleation clusters that become aerosols is challenging, requiring a comprehensive analysis of configurational space and highly accurate Gibbs free energy calculations. We identified a large set of minimum energy configurations of (H₂SO₄)₃ using a sampling technique based on a genetic algorithm and a stepwise density functional theory (DFT) approach and computed the thermodynamics of formation of these configurations with more accurate wavefunction-based electronic energies computed on the DFT geometries. The DLPNO-CCSD(T) methods always return more positive energies compared to the DFT energies. Within the DLPNO-CCSD(T) methods, extrapolating to the complete basis set limit gives more positive free energies compared to explicitly correlated single-point energies. The CBS extrapolation was shown to be robust as both the 4-5 inverse polynomial and Riemann zeta function schemes were within chemical accuracy of one another.

Atmospheric Chemistry of CH2OO: The Hydrolysis of CH2OO in Small Clusters of Sulfuric Acid

The hydrolysis of CH₂OO is not only a dominant sink for the CH₂OO intermediate in the atmosphere but also a key process in the formation of aerosols. Herein, the reaction mechanism and kinetics for the hydrolysis of CH₂OO catalyzed by the precursors of atmospheric aerosols, including H₂SO₄, H₂SO₄···H₂O, and (H₂SO₄)₂, have been studied theoretically at the CCSD(T)-F12a/cc-pVDZ-F12//B3LYP/6-311+G(2df,2pd) level. The calculated results show that the three catalysts decrease the energy barrier by over 10.3 kcal·mol^-1; at the same time, the product formation of HOCH₂OOH is more strongly bonded to the three catalysts than to the reactants CH₂OO and H₂O, revealing that small clusters of sulfuric acid promote the hydrolysis of CH₂OO both kinetically and thermodynamically. Kinetic simulations show that the H₂SO₄-assisted reaction is more favorable than the H₂SO₄···H₂O- (the pseudo-first-order rate constant being 27.9-11.5 times larger) and (H₂SO₄)₂- (between 2.8 × 10⁴ and 3.4 × 10⁵ times larger) catalyzed reactions. Additionally, due to relatively lower concentration of H₂SO₄, the hydrolysis of CH₂OO with H₂SO₄ cannot compete with the CH₂OO + H₂O or (H₂O)₂ reaction within the temperature range of 280-320 K, since its pseudo-first-order rate ratio is smaller by 4-7 or 6-8 orders of magnitude, respectively. However, the present results provide a good example of how small clusters of sulfuric acid catalyze the hydrolysis of an important atmospheric species.

Hydration motifs of ammonium bisulfate clusters show complex temperature dependence

The role of water in the formation of particles from atmospheric trace gases is not well understood, in large part due to difficulties in detecting its presence under atmospheric conditions and the variety of possible structures that must be screened computationally. Here, we use infrared spectroscopy and variable-temperature ion trap mass spectrometry to investigate the structural motifs adopted by water bound to ammonium bisulfate clusters and their temperature dependence. For clusters featuring only acid-base linkages, water adopts a bridging arrangement spanning an adjacent ammonium and bisulfate. For larger clusters, water can also insert into a bisulfate-bisulfate hydrogen bond, yielding hydration isomers with very similar binding energies. The population of these isomers shows a complex temperature evolution, as an apparent third isomer appears with a temperature dependence that is difficult to explain using simple thermodynamic arguments. These observations suggest that the thermodynamics of water binding to atmospheric clusters such as these may not be straightforward.

Theoretical study of the hydration effects on alkylamine and alkanolamine clusters and the atmospheric implication

Amines are important atmospheric nucleation precursors in polluted areas. However, specific roles of various amines in enhancing the stability of pre-nucleation clusters are poorly understood. Herein, different roles of trimethylamine (TMA) and monoethanolamine (MEA) in the formation of sulfuric acid (SA)-based pre-nucleation clusters were investigated. The hydration effects of up to four water (W) molecules on the interaction of the acid-base pairs of (TMA)(SA) and (MEA)(SA) were computationally investigated at the M06-2X/6-311++G (3df, 3pd) level of theory. Results show that the formation thermodynamics of key intermediate clusters are different with amines. Besides, MEA-enhanced formation of pre-nucleation clusters plays an important role in environments with high humidity while TMA may be the dominant nucleation precursors in dry conditions. The concentration of dominant MEA-containing pre-nucleation cluster is at least three orders of magnitude higher than that of TMA and dimethylamine near emission sources at conditions of T = 298.15 K and RH = 60%, indicating that alkanolamine may play an important role in atmospheric nucleation. Furthermore, the hydration of MEA is easier than that of alkylamines. Our results put forward the need to distinguish the performance of different types of amine under specific conditions to better model the new particle formation events in highly polluted areas. Besides, this study indicates that alkanolamines such as MEA are important participator in new particle formation.

Understanding Hygroscopic Nucleation of Sulfate Aerosols: Combination of Molecular Dynamics Simulation with Classical Nucleation Theory

We present a combined molecular dynamics (MD) and classical nucleation theory (CNT) approach to address many issues regarding the nucleation of inorganic aerosols. By taking parameters from MD simulations, we find the CNT predicts fairly reasonable free-energy profiles for the hygroscopic nucleation of aerosols. Moreover, we find that the ionization of sulfates can play a key role in stabilizing aqueous clusters and that both the size of the critical nucleus and the nucleation barrier can be significantly lowered by the H₂SO₄ and NH₄HSO₄, whereas the effect of NH₃ on nucleation is negligible. NH₄HSO₄ provides stronger enhancement effect to aerosol formation than H₂SO₄. In view of the consistency between the theoretical prediction and experimental observation, the combination of MD simulation and CNT appears to be a valuable approach to gain deeper understanding of how aerosol nucleation is affected by different chemical species.

A molecular-scale study on the hydration of sulfuric acid-amide complexes and the atmospheric implication

Amides are ubiquitous in atmosphere. However, the role of amides in new particle formation (NPF) is poorly understood. Herein, the interaction of urea and formamide with sulfuric acid (SA) and up to four water (W) molecules has been studied at the M06-2X/6-311++G(3df,3pd) level of theory. The structures and properties of (Formamide)(SA)(W)_n (n = 0-4) and (Urea)(SA)(W)_n (n = 0-4) clusters were investigated. Results show that the interaction of SA with the CO group of amides plays a more important role in amide clusters compared with the NH₂ group. Proton transfer to water molecule become dominant in highly hydrated amide clusters at lower temperatures. There is no proton transfer to CO group in formamide clusters. The Rayleigh light scattering intensities of amide clusters are comparable to that of amine and oxalic acid clusters reported previously. Moreover, unhydrated (Amide)(SA) clusters have similar or even higher ability than hydrated SA clusters to participate in ion-induced nucleation. In comparison with formamide, urea has more interacting sites and its clusters have higher Rayleigh light scattering intensities, larger dipole moment, stronger interaction with SA and lower water affinity. The intermolecular interaction in (Formamide)(SA) is slightly weaker than that of SA dimer, which may be compensated by the high concentration of formamide, thus enabling formamide to participate in initial steps of NPF. This study may bring new insight into the role of amides in initial steps of NPF from molecular scale and could help better understand the properties of amide-containing organic aerosol.

Molecular understanding of the interaction of amino acids with sulfuric acid in the presence of water and the atmospheric implication

Amino acids are important components of atmospheric aerosols. Despite the diversity of amino acids structures, however, the role of amino acids with additional non-characteristic functional groups in new particle formation (NPF) has almost remained unexplored. Herein, the interaction of serine (Ser) and threonine (Thr), which feature a hydroxyl group and differ by a methyl-substitution, with sulfuric acid (SA) and up to three water (W) molecules has been investigated at the M06-2X/6-311++G (3df, 3pd) level of theory. The effects of structural differences of amino acids on the structure and properties of clusters were also pointed out. Results show that serine may play more important role in stabilizing sulfuric acid to promote NPF in initial steps compared with threonine, glycine and alanine. Meanwhile, threonine may participate in ion-induced nucleation due to the high dipole moment of (Thr) (SA) isomers. Moreover, the effects of structure differences of amino acids can be seen in several aspects. Firstly, methyl substitution and hydroxyl group of amino acids have great influence on the structure of clusters. Secondly, hydrated (Ser) (SA) and (Tur) (SA) clusters could retain water even at low relative humidity, which may due to the hydroxyl group in serine and threonine. In addition, the Rayleigh light scattering intensities of amino acid-containing clusters are higher than trimethylamine, monoethanolamine and oxalic acid-involved counterparts. The effect of carboxyl group and methyl substitution on optical properties of clusters is also discussed. This study may bring new insight into the role of amino acids with additional non-characteristic functional groups in initial steps of NPF.

Mass spectrometry of aerosol particle analogues in molecular beam experiments

Nanometer-size particles such as ultrafine aerosol particles, ice nanoparticles, water nanodroplets, etc, play an important, however, not yet fully understood role in the atmospheric chemistry and physics. These species are often composed of water with admixture of other atmospherically relevant molecules. To mimic and investigate such particles in laboratory experiments, mixed water clusters with atmospherically relevant molecules can be generated in molecular beams and studied by various mass spectrometric methods. The present review demonstrates that such experiments can provide unprecedented details of reaction mechanisms, and detailed insight into the photon-, electron-, and ion-induced processes relevant to the atmospheric chemistry. After a brief outline of the molecular beam preparation, cluster properties, and ionization methods, we focus on the mixed clusters with various atmospheric molecules, such as hydrated sulfuric acid and nitric acid clusters, N_x O_y and halogen-containing molecules with water. A special attention is paid to their reactivity and solvent effects of water molecules on the observed processes.

Hydration of Atmospheric Molecular Clusters: A New Method for Systematic Configurational Sampling

We present a new systematic configurational sampling algorithm for investigating the potential energy surface of hydrated atmospheric molecular clusters. The algorithm is based on creating a Fibonacci sphere around each atom in the cluster and adding water molecules to each point in nine different orientations. For the sampling of water molecules to existing hydrogen bonds, the cluster is displaced along the hydrogen bond, and a water molecule is placed in between in three different orientations. Generated redundant structures are eliminated based on minimizing the root-mean-square distance of different conformers. Initially, the clusters are sampled using the semiempirical PM6 method and subsequently using density functional theory (M06-2X and ωB97X-D) with the 6-31++G(d,p) basis set. Applying the developed algorithm, we study the hydration of sulfuric acid with up to 15 water molecules. We find that the addition of the first four water molecules "saturate" the sulfuric acid molecule and that they are more thermodynamically favorable than the addition of water molecules 5-15. Using the large generated set of conformers, we assess the performance of approximate methods (ωB97X-D, M06-2X, PW91, and PW6B95-D3) in calculating the binding energies and assigning the global minimum conformation compared to high level CCSD(T)-F12a/VDZ-F12 reference calculations. The tested DFT functionals systematically overestimate the binding energies compared to coupled cluster calculations, and we find that this deficiency can be corrected by a simple scaling factor.

Closed-Shell Organic Compounds Might Form Dimers at the Surface of Molecular Clusters

The role of covalently bound dimer formation is studied using high-level quantum chemical methods. Reaction free energy profiles for dimer formation between common oxygen-containing functional groups are calculated, and based on the Gibbs free energy differences between transition states and reactants, we show that none of the studied two-component gas-phase reactions are kinetically feasible at 298.15 K and 1 atm. Therefore, the catalyzing effect of water, base, or acid molecules is calculated, and sulfuric acid is identified to lower the activation free energies significantly. We find that the reactions yielding hemiacetal, peroxyhemiacetal, α-hydroxyester, and geminal diol products occur with activation free energies of less than 10 kcal/mol with sulfuric acid as a catalyst, indicating that these reactions could potentially take place on the surface of sulfuric acid clusters. Additionally, the formed dimer products bind stronger onto the pre-existing cluster than the corresponding reagent monomers do. This implies that covalent dimerization reactions stabilize the existing cluster thermodynamically and make it less likely to evaporate. However, the studied small organic compounds, which contain only one functional group, are not able to form dimer products that are stable against evaporation at atmospheric conditions. Calculations of dimer formation onto a cluster surface and the clustering ability of dimer products should be extended to large terpene oxidation products in order to estimate the real atmospheric significance.

Mechanistic and kinetic study on the catalytic hydrolysis of COS in small clusters of sulfuric acid

The catalytic hydrolysis of carbonyl sulfide (COS) and the effect of small clusters of H₂O and H₂SO₄ have been studied by theoretical calculations. The addition of H₂SO₄ could increase the enthalpy change (ΔH<0) and decrease relative energy of products (relative energy<0), resulting in hydrolysis reaction changed from an endothermic reaction to an exothermic reaction. Further, H₂SO₄ decreases the energy barrier by 5.25 kcal/mol, and it enhances the catalytic hydrolysis through the hydrogen transfer effect. The (COS + H₂SO₄-H₂O) reaction has the lowest energy barrier of 29.97 kcal/mol. Although an excess addition of H₂O and H₂SO₄ increases the energy barrier, decreases the catalytic hydrolysis, which is consistent with experimental observations. The order of the energy barriers for the three reactions from low to high are as follows: COS + H₂SO₄-H₂O < COS + H₂O + H₂SO₄-H₂O < COS + H₂O+(H₂SO₄)₂. Kinetic simulations show that the addition of H₂SO₄ can increase the reaction rate constants. Consequently, adding an appropriate amount of sulfuric acid promotes the catalytic hydrolysis of COS both kinetically and thermodynamically.

Hydration of the methanesulfonate–ammonia/amine complex and its atmospheric implications

Methanesulfonate (MSA⁻), found in substantial concentrations in the atmosphere, is expected to enhance aerosol nucleation and the growth of nanoparticles, but the details of methanesulfonate clusters are poorly understood. In this study, MSA⁻ was chosen along with ammonia (NH₃) or three common amines and water (H₂O) to discuss the roles of ternary homogeneous nucleation and ion-induced nucleation in aerosol formation. We studied the structural characteristics and thermodynamics of the clusters using density functional theory at the PW91PW91/6-311++G(3df,3pd) level. The analysis of noncovalent interactions predicts that the amines can form more stable clusters with MSA⁻ than NH₃, in agreement with the results from structures and thermodynamics; however, the enhancement in stability for amines is not large enough to overcome the difference in the concentrations of NH₃ and amines under typical atmospheric conditions. In addition, the favorable free energies of formation for the (MSA⁻)(NH₃/amines)(H₂O)_n (n = 0–3) clusters at 298.15 K show that MSA⁻ could contribute to the aerosol nucleation process with binding NH₃/amines and H₂O up to n = 3. There are strong temperature and humidity dependences for the formation of complexes; higher humidity and temperature promote the formation of larger hydrates. Finally, for the (MSA⁻)(NH₃/amines)(H₂O)_n clusters, the evaporation rates were determined to further investigate the atmospheric implications.

Methanesulfonate (MSA⁻), found in substantial concentrations in the atmosphere, is expected to enhance aerosol nucleation and the growth of nanoparticles, but the details of methanesulfonate clusters are poorly understood.

Nanoparticles grown from methanesulfonic acid and methylamine: microscopic structures and formation mechanism

Mechanisms of particle formation and growth in the atmosphere are of great interest due to their impacts on climate, health and visibility. However, the microscopic structures and related properties of the smallest nanoparticles are not known. In this paper we pursue computationally a microscopic description for the formation and growth of methanesulfonic acid (MSA) and methylamine (MA) particles under dry conditions. Energetic and dynamics simulations were used to assess the stabilities of proposed model structures for these particles. Density functional theory (DFT) and semi-empirical (PM3) calculations suggest that (MSA-MA)₄ is a major intermediate in the growth process, with the dissociation energies, enthalpies and free energies indicating considerable stability for this cluster. Dynamics simulations show that this species is stable for at least 100 ps at temperatures up to 500 K, well above atmospheric temperatures. In order to reach experimentally detectable sizes (>1.4 nm), continuing growth is suggested to occur via clustering of (MSA-MA)₄. The dimer (MSA-MA)₄(MSA-MA)₄ may be one of the smaller experimentally measured particles. Step by step addition of MSA to (MSA-MA)₄, is also a likely potential growth mechanism when MSA is excess. In addition, an MSA-MA crystal is predicted to exist. These studies demonstrate that computations of particle structure and dynamics in the nano-size range can be useful for molecular level understanding of processes that grow clusters into detectable particles.

Does H4SO5 exist?

The possible existence of H₄SO₅ in aqueous sulfuric acid is analyzed in detail. For bare H₄SO₅, the computed free energy barrier for the exergonic transformation of H₄SO₅ into the H₂SO₄H₂O complex is only 3.8 kcal mol^-1. The presence of water or sulfuric acid catalyzes the dehydration to such an extent that it becomes almost a barrierless process. In the gas phase, dehydration of H₄SO₅ is an autocatalytic reaction as the water molecule produced by the decomposition of one H₄SO₅ molecule induces further dissociation. Thus, in solution, the surrounding water molecules make the para-sulfuric acid a very vulnerable species to exist. The simulated Raman spectra also corroborate the absence of H₄SO₅ in solution.

Contribution of methane sulfonic acid to new particle formation in the atmosphere

Methane sulfonic acid (MSA) is present in substantial concentrations in the gas phase over oceans and coastal regions. We present an investigation into the contribution of MSA to new particle formation with the common atmospheric aerosol nucleation precursors including MSA, methanol, formic acid, acetone, dimethylether, formaldehyde, methyl formate, by making use a quantum chemical approach. Density functional theory calculations indicate that these bimolecular complexes are characterized by the presence of strong inter-molecular hydrogen bonds (SOH⋯O) with large binding energies and thermodynamic equilibrium constants. Topological analysis employing quantum theory of atoms in molecules shows that the charge density of the SOH⋯O hydrogen bonds of the MSA complexes falls in the range of hydrogen bonding criteria, but the Laplacian at bond critical points exceeds the range, which is due to the strong hydrogen bonding interactions. In all the studied complexes, the electrostatic interactions are found to be the main attractive force by localized molecular orbital energy decomposition analysis. All these indicate the environmental fate of MSA could play the role of nucleation centers in new particle formation. The effect of the atmospheric heights (0-12 km) was also considered. The Gibbs free energy of formation decreases with the increase of the atmospheric height owing to the decrease of the atmospheric temperature and pressure. The calculated Gibbs free energies of formation within the atmospheric temperature and pressure range could help to understand the atmospheric pollution.

Formation of atmospheric molecular clusters consisting of sulfuric acid and C8H12O6 tricarboxylic acid

We present the structures and thermochemical properties of (MBTCA)₁₋₃(H₂SO₄)₁₋₄ atmospheric molecular clusters.

Basis set convergence of the binding energies of strongly hydrogen-bonded atmospheric clusters

We present the first binding energy benchmark set at the CBS limit of strongly hydrogen bonded atmospheric molecular clusters.

Particle formation and growth from oxalic acid, methanesulfonic acid, trimethylamine and water: a combined experimental and theoretical study

A combined experimental-theoretical study on the effect of oxalic acid on particle formation and growth from the reaction of MSA with trimethylamine in the absence and presence of water.

Hydrogen Bonding Interaction between Atmospheric Gaseous Amides and Methanol

Amides are important atmospheric organic–nitrogen compounds. Hydrogen bonded complexes of methanol (MeOH) with amides (formamide, N-methylformamide, N,N-dimethylformamide, acetamide, N-methylacetamide and N,N-dimethylacetamide) have been investigated. The carbonyl oxygen of the amides behaves as a hydrogen bond acceptor and the NH group of the amides acts as a hydrogen bond donor. The dominant hydrogen bonding interaction occurs between the carbonyl oxygen and the OH group of methanol as well as the interaction between the NH group of amides and the oxygen of methanol. However, the hydrogen bonds between the CH group and the carbonyl oxygen or the oxygen of methanol are also important for the overall stability of the complexes. Comparable red shifts of the C=O, NH- and OH-stretching transitions were found in these MeOH–amide complexes with considerable intensity enhancement. Topological analysis shows that the electron density at the bond critical points of the complexes fall in the range of hydrogen bonding criteria, and the Laplacian of charge density of the O–H∙∙∙O hydrogen bond slightly exceeds the upper value of the Laplacian criteria. The energy decomposition analysis further suggests that the hydrogen bonding interaction energies can be mainly attributed to the electrostatic, exchange and dispersion components.

Effects of Global and Local Anharmonicities on the Thermodynamic Properties of Sulfuric Acid Monohydrate

We use state-of-the-art electronic structure calculation methods and large basis sets to obtain reliable values for the thermodynamic properties of sulfuric acid monohydrate and study the effects of vibrational anharmonicity on these properties. We distinguish between two forms of vibrational anharmonicity: local anharmonicity, which refers to the anharmonicity of the vibrational modes of a given cluster conformer, and global anharmonicity, which originates from accounting for the presence of different conformers in the first place. In our most accurate approach, we solve the nuclear Schrödinger equation variationally for the intermolecular large-amplitude motions, thus quantum-mechanically accounting for the presence of higher-energy conformers for both reactants and products, while using the standard vibrational perturbational approach for the other vibrational modes. This results in a value of -11.0 kJ/mol for the reaction Gibbs free energy at 298.15 K. When standard vibrational perturbational approaches are employed, the effects of local anharmonicity depend heavily on the choice of the electronic structure calculation basis set. In fact, better results can often be achieved by combining a simple harmonic treatment for the vibrational partition function with a statistical mechanical accounting of global anharmonicity. Thus, we recommend that future studies that intend to include anharmonicity start by accounting for the presence of higher-energy conformers and only then consider whether local anharmonicity calculations are feasible and necessary.

Sulphur Kβ emission spectra reveal protonation states of aqueous sulfuric acid

In this paper we report an X-ray emission study of bulk aqueous sulfuric acid. Throughout the range of molarities from 1 M to 18 M the sulfur Kβ emission spectra from H2SO4 (aq) depend on the molar fractions and related deprotonation of H2SO4. We compare the experimental results with results from emission spectrum calculations based on atomic structures of single molecules and structures from ab initio molecular dynamics simulations. We show that the S Kβ emission spectrum is a sensitive probe of the protonation state of the acid molecules. Using non-negative matrix factorization we are able to extract the fractions of different protonation states in the spectra, and the results are in good agreement with the simulation for the higher part of the concentration range.

Characterization of the nucleation precursor (H2SO4–(CH3)2NH) complex: intra-cluster interactions and atmospheric relevance

Amines have been proposed to participate in the nucleation process, but the electron density analysis and the determination of a temperature dependence of the clusters are still lacking.

Disulfuric acid dissociated by two water molecules: ab initio and density functional theory calculations

We have studied geometries, energies and vibrational spectra of disulfuric acid (H2S2O7) and its anion (HS2O7(-)) hydrated by a few water molecules, using density functional theory (M062X) and ab initio theory (SCS-MP2 and CCSD(T)). The most noteworthy result is found in H2S2O7(H2O)2 in which the lowest energy conformer shows deprotonated H2S2O7. Thus, H2S2O7 requires only two water molecules, the fewest number of water molecules for deprotonation among various hydrated monomeric acids reported so far. Even the second deprotonation of the first deprotonated species HS2O7(-) needs only four water molecules. The deprotonation is supported by vibration spectra, in which acid O-H stretching peaks disappear and specific three O-H stretching peaks for H3O(+) (eigen structure) appear. We have also kept track of variations in several geometrical parameters, atomic charges, and hybrid orbital characters upon addition of water. As the number of water molecules added increases, the S-O bond weakens in the case of H2S2O7, but strengthens in the case of HS2O7(-). It implies that the decomposition leading to H2SO4 and SO3 hardly occurs prior to the 2nd deprotonation at low temperatures.

Hydrolysis of Sulfur Dioxide in Small Clusters of Sulfuric Acid: Mechanistic and Kinetic Study

The deposition and hydrolysis reaction of SO2 + H2O in small clusters of sulfuric acid and water are studied by theoretical calculations of the molecular clusters SO2-(H2SO4)n-(H2O)m (m = 1,2; n = 1,2). Sulfuric acid exhibits a dramatic catalytic effect on the hydrolysis reaction of SO2 as it lowers the energy barrier by over 20 kcal/mol. The reaction with monohydrated sulfuric acid (SO2 + H2O + H2SO4 - H2O) has the lowest energy barrier of 3.83 kcal/mol, in which the cluster H2SO4-(H2O)2 forms initially at the entrance channel. The energy barriers for the three hydrolysis reactions are in the order SO2 + (H2SO4)-H2O > SO2 + (H2SO4)2-H2O > SO2 + H2SO4-H2O. Furthermore, sulfurous acid is more strongly bonded to the hydrated sulfuric acid (or dimer) clusters than the corresponding reactant (monohydrated SO2). Consequently, sulfuric acid promotes the hydrolysis of SO2 both kinetically and thermodynamically. Kinetics simulations have been performed to study the importance of these reactions in the reduction of atmospheric SO2. The results will give a new insight on how the pre-existing aerosols catalyze the hydrolysis of SO2, leading to the formation and growth of new particles.