Effect of Mixing Ammonia and Alkylamines on Sulfate Aerosol Formation

Photochemical Mechanisms in Atmospherically Relevant Iodine Oxide Clusters

Atmospheric new particle formation events can be driven by iodine oxides or oxoacids via both neutral and ionic mechanisms. Photolysis of new particles likely plays a significant role in their growth mechanisms, but their spectra and photolysis mechanisms remain difficult to characterize. We recorded ultraviolet (UV) photodissociation spectra of (I2O5)0-3(IO3-) clusters, observing loss of an O atom, I2O4, and (I2O5)1,2 in the atmospherically relevant range of 300-340 nm. With increasing cluster size, the intensity of absorption red shifts and generally increases, suggesting particles photolyze more frequently as they grow. Estimates of the rates indicate that even relatively small clusters are likely to undergo photolysis under high-UV conditions. Vibrational spectra identify the covalent moiety I3O8- as the likely chromophore, not IO3-. The I2O5 loss pathway competes with particle growth, while the slower O loss pathway likely produces 3O + 3(cluster) products that could drive subsequent intraparticle chemistry, particularly with co-adsorbed organic or amine species.

Stability and Structure of Potentially Atmospherically Relevant Glycine Ammonium Bisulfate Clusters

New particle formation (NPF) is the process by which trace atmospheric acids and bases cluster and grow into particles that ultimately impact climate. Sulfuric acid concentration drives NPF, but nitrogen-containing bases promote the formation of more stable clusters via salt bridge formation. Recent computational efforts have suggested that amino acids can enhance NPF, predicting that they can stabilize new particles via multiple protonation sites, but there has yet to be experimental validation of these predictions. We used mass spectrometry and infrared spectroscopy to study the structure and stability of cationic clusters composed of glycine, sulfuric acid, and ammonia. When collisionally activated, clusters were significantly more likely to eliminate ammonia or sulfuric acid than glycine, while quantum chemical calculations predicted lower binding free energies for ammonia but similar binding free energies for glycine and sulfuric acid. These calculations predicted several low-energy structures, so we compared experimental and computed vibrational spectra to attempt to validate the computationally predicted minimum energy structure. Unambiguous identification of the experimental structure by comparison to these calculations was made difficult by the complexity of the experimental spectra and the fact that the identity of the computed lowest-energy structure depended strongly on temperature. If their vapors are present, amino acids are likely to be enriched in new particles by displacing more weakly bound ammonia, similar to the behavior of other atmospheric amines. The carboxylic acid groups were found to preferentially interact with other carboxylic acids, suggesting incipient organic/inorganic phase separation even at these small sizes.

Computational Tools for Handling Molecular Clusters: Configurational Sampling, Storage, Analysis, and Machine Learning

Computational modeling of atmospheric molecular clusters requires a comprehensive understanding of their complex configurational spaces, interaction patterns, stabilities against fragmentation, and even dynamic behaviors. To address these needs, we introduce the Jammy Key framework, a collection of automated scripts that facilitate and streamline molecular cluster modeling workflows. Jammy Key handles file manipulations between varieties of integrated third-party programs. The framework is divided into three main functionalities: (1) Jammy Key for configurational sampling (JKCS) to perform systematic configurational sampling of molecular clusters, (2) Jammy Key for quantum chemistry (JKQC) to analyze commonly used quantum chemistry output files and facilitate database construction, handling, and analysis, and (3) Jammy Key for machine learning (JKML) to manage machine learning methods in optimizing molecular cluster modeling. This automation and machine learning utilization significantly reduces manual labor, greatly speeds up the search for molecular cluster configurations, and thus increases the number of systems that can be studied. Following the example of the Atmospheric Cluster Database (ACDB) of Elm (ACS Omega, 4, 10965-10984, 2019), the molecular clusters modeled in our group using the Jammy Key framework have been stored in an improved online GitHub repository named ACDB 2.0. In this work, we present the Jammy Key package alongside its assorted applications, which underline its versatility. Using several illustrative examples, we discuss how to choose appropriate combinations of methodologies for treating particular cluster types, including reactive, multicomponent, charged, or radical clusters, as well as clusters containing flexible or multiconformer monomers or heavy atoms. Finally, we present a detailed example of using the tools for atmospheric acid-base clusters.

Current and future machine learning approaches for modeling atmospheric cluster formation

The formation of strongly bound atmospheric molecular clusters is the first step towards forming new aerosol particles. Recent advances in the application of machine learning models open an enormous opportunity for complementing expensive quantum chemical calculations with efficient machine learning predictions. In this Perspective, we present how data-driven approaches can be applied to accelerate cluster configurational sampling, thereby greatly increasing the number of chemically relevant systems that can be covered.

Atmospheric Sulfuric Acid-Multi-Base New Particle Formation Revealed through Quantum Chemistry Enhanced by Machine Learning

The formation of molecular clusters and secondary aerosols in the atmosphere has a significant impact on the climate. Studies typically focus on the new particle formation (NPF) of sulfuric acid (SA) with a single base molecule (e.g., dimethylamine or ammonia). In this work, we examine the combinations and synergy of several bases. Specifically, we used computational quantum chemistry to perform configurational sampling (CS) of (SA)_0-4(base)_0-4 clusters with five different types of bases: ammonia (AM), methylamine (MA), dimethylamine (DMA), trimethylamine (TMA), and ethylenediamine (EDA). Overall, we studied 316 different clusters. We used a traditional multilevel funnelling sampling approach augmented by a machine-learning (ML) step. The ML made the CS of these clusters possible by significantly enhancing the speed and quality of the search for the lowest free energy configurations. Subsequently, the cluster thermodynamics properties were evaluated at the DLPNO-CCSD(T₀)/aug-cc-pVTZ//ωB97X-D/6-31++G(d,p) level of theory. The calculated binding free energies were used to evaluate the cluster stabilities for population dynamics simulations. The resultant SA-driven NPF rates and synergies of the studied bases are presented to show that DMA and EDA act as nucleators (although EDA becomes weak in large clusters), TMA acts as a catalyzer, and AM/MA is often overshadowed by strong bases.

The Competition between Hydrogen, Halogen, and Covalent Bonding in Atmospherically Relevant Ammonium Iodate Clusters

Iodine-containing clusters are expected to be central to new particle formation (NPF) events in polar and midlatitude coastal regions. Iodine oxoacids and iodine oxides are observed in newly formed clusters, and in more polluted midlatitude settings, theoretical studies suggest ammonia may increase growth rates. Structural information was obtained via infrared (IR) spectroscopy and quantum chemical calculations for a series of clusters containing ammonia, iodic acid, and iodine pentoxide. Structures for five of the smallest cationic clusters present in the mass spectrum were identified, and four of the structures were found to preferentially form halogen and/or covalent bonds over hydrogen bonds. Ammonia is important in proton transfer from iodic acid components and also provides a scaffold to template the formation of a halogen and covalent bonded backbone. The calculations executed for the two largest clusters studied suggested the formation of a covalent I₃O8- anion within the clusters.

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.

Clusteromics IV: The Role of Nitric Acid in Atmospheric Cluster Formation

Nitric acid (NA) has previously been shown to affect atmospheric new particle formation; however, its role still remains highly uncertain. Through the employment of state-of-the-art quantum chemical methods, we study the (acid)_1-2(base)_1-2 and (acid)₃(base)₂ clusters containing at least one nitric acid (NA) and sulfuric acid (SA) or methanesulfonic acid (MSA) with bases ammonia (A), methylamine (MA), dimethylamine (DMA), trimethylamine (TMA), and ethylenediamine (EDA). The initial cluster configurations are generated using the ABCluster program. PM7 and ωB97X-D/6-31++G(d,p) calculations are used to reduce the number of relevant configurations. The thermochemical parameters are calculated at the ωB97X-D/6-31++G(d,p) level of theory with the quasi-harmonic approximation, and the final single-point energies are calculated with high-level DLPNO-CCSD(T₀)/aug-cc-pVTZ calculations. The enhancing effect from the presence of nitric acid on cluster formation is studied using the calculated thermochemical data and cluster dynamics simulations. We find that when NA is in excess compared with the other acids, it has a substantial enhancing effect on the cluster formation potential.

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.

Clusteromics III: Acid Synergy in Sulfuric Acid-Methanesulfonic Acid-Base Cluster Formation

Acid-base molecular clusters are an important stage in atmospheric new particle formation. While such clusters are most likely multicomponent in nature, there are very few reports on clusters consisting of multiple acid molecules and multiple base molecules. By applying state-of-the-art quantum chemical methods, we herein study electrically neutral (SA)₁(MSA)₁(base)_0-2 clusters with base = ammonia (A), methylamine (MA), dimethylamine (DMA), trimethylamine (TMA) and ethylenediamine (EDA). The cluster structures are obtained using a funneling approach employing the ABCluster program, semiempirical PM7 calculations and ωB97X-D/6-31++G(d,p) calculations. The final binding free energies are calculated at the DLPNO-CCSD(T₀)/aug-cc-pVTZ//ωB97X-D/6-31++G(d,p) level of theory using the quasi-harmonic approximation. Based on the calculated cluster geometries and thermochemistry (at 298.15 K and 1 atm), we find that the mixed (SA)₁(MSA)₁(base)_1-2 clusters more resemble the (SA)₂(base)_1-2 clusters compared to the (MSA)₂(base)_1-2 clusters. Hence, some of the steric hindrance and lack of hydrogen bond capacity previously observed in the (MSA)₂(base)_1-2 clusters is diminished in the corresponding (SA)₁(MSA)₁(base)_1-2 clusters. Cluster kinetics simulations reveal that the presence of an MSA molecule in the clusters enhances the cluster formation potential by up to a factor of 20. We find that the SA-MSA-DMA clusters have the highest cluster formation potential, and thus, this system should be further extended to larger sizes in future studies.

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.

Tri-Base Synergy in Sulfuric Acid-Base Clusters

Synergistic effects between different bases can greatly enhance atmospheric sulfuric acid (SA)-base cluster formation. However, only the synergy between two base components has previously been investigated. Here, we extend this concept to three bases by studying large atmospherically relevant (SA)3(base)3 clusters, with the bases ammonia (A), methylamine (MA), dimethylamine (DMA), trimethylamine (TMA) and ethylenediamine (EDA). Using density functional theory—ωB97X-D/6-31++G(d,p)—we calculate the cluster structures and vibrational frequencies. The thermochemical parameters are calculated at 29,815 K and 1 atm, using the quasi-harmonic approximation. The binding energies of the clusters are calculated using high level DLPNO-CCSD(T0)/aug-cc-pVTZ. We find that the cluster stability in general depends on the basicity of the constituent bases, with some noteworthy additional guidelines: DMA enhances the cluster stability, TMA decreases the cluster stability and there is high synergy between DMA and EDA. Based on our calculations, we find it highly likely that three, or potentially more, different bases, are involved in the growth pathways of sulfuric acid-base clusters.

Clusteromics II: Methanesulfonic Acid-Base Cluster Formation

The role of methanesulfonic acid (MSA) in atmospheric new particle formation remains highly uncertain. Using state-of-the-art computational methods, we study the electrically neutral (MSA)0-2(base)0-2 clusters, with base = ammonia (A), methylamine (MA), dimethylamine (DMA), trimethylamine (TMA), and ethylenediamine (EDA). The cluster configurations are obtained using the ABCluster program and the number of initial cluster configurations is reduced based on PM7 calculations. Thermochemical parameters are calculated using the quasi-harmonic approximation based on the ωB97X-D/6-31++G(d,p) cluster structures and vibrational frequencies. The single point energies are calculated at the DLPNO-CCSD(T0)/aug-cc-pVTZ level of theory. We find that MSA shows a different interaction pattern with the bases compared to sulfuric acid and does not simply follow the basicity of the bases for these small clusters. In all cases, we find that the MSA-base clusters show very low cluster formation potential, indicating that electrically neutral clusters consisting solely of MSA as the clustering acid are most likely not capable of forming and growing under realistic atmospheric conditions.

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.

Clusteromics I: Principles, Protocols, and Applications to Sulfuric Acid-Base Cluster Formation

We recently coined the term clusteromics as a holistic approach for obtaining insight into the chemical complexity of atmospheric molecular cluster formation and at the same time providing the foundation for thermochemical databases that can be utilized for developing machine learning models. Here, we present the first paper in the series that applies state-of-the-art computational methods to study multicomponent (SA)0-2(base)0-2 clusters, with SA = sulfuric acid and base = [ammonia (A), methylamine (MA), dimethylamine (DMA), trimethylamine (TMA), and ethylenediamine (EDA)] with all combinations of the five bases. The initial cluster configurations are obtained using the ABCluster program and the number of relevant configurations are reduced based on PM7 and ωB97X-D/6-31++G(d,p) calculations. Thermochemical parameters are calculated based on the ωB97X-D/6-31++G(d,p) cluster structures and vibrational frequencies using the quasi-harmonic approximation. The single-point energies are refined with a high-level DLPNO-CCSD(T0)/aug-cc-pVTZ calculation. Using the calculated thermochemical data, we perform kinetics simulations to evaluate the potential of these small (SA)0-2(base)0-2 clusters to grow into larger cluster sizes. In all cases we find that having more than one type of base molecule present in the cluster will increase the potential for forming larger clusters primarily due to the increased available vapor concentration.

Toward a Holistic Understanding of the Formation and Growth of Atmospheric Molecular Clusters: A Quantum Machine Learning Perspective

The formation of atmospheric molecular clusters is an important stage in forming new particles in the atmosphere. Despite being a highly focused research area, the exact chemical species involved in the initial steps in new particle formation remain elusive. In this Perspective the main challenges and recent progression in the field are outlined with a special emphasis on the chemical complexity of the puzzle and prospect of modeling larger clusters. In general, there is a high demand for accurate and more complete quantum chemical data sets that can be applied in cluster distribution dynamics models and coupled to atmospheric chemical transport models. A view on how the community could reach this goal by applying data-driven machine learning approaches for more efficient exploration of cluster configurations is presented. A path toward larger clusters and direct molecular dynamics simulations of cluster formation and growth using machine learning models is discussed.

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.

Establishing the structural motifs present in small ammonium and aminium bisulfate clusters of relevance to atmospheric new particle formation

Atmospheric new particle formation is the process by which atmospheric trace gases, typically acids and bases, cluster and grow into potentially climatically relevant particles. Here, we evaluate the structures and structural motifs present in small cationic ammonium and aminium bisulfate clusters that have been studied both experimentally and computationally as seeds for new particles. For several previously studied clusters, multiple different minimum-energy structures have been predicted. Vibrational spectra of mass-selected clusters and quantum chemical calculations allow us to assign the minimum-energy structure for the smallest cationic cluster of two ammonium ions and one bisulfate ion to a C_S-symmetry structure that is persistent under amine substitution. We derive phenomenological vibrational frequency scaling factors for key bisulfate vibrations to aid in the comparison of experimental and computed spectra of larger clusters. Finally, we identify a previously unassigned spectral marker for intermolecular bisulfate-bisulfate hydrogen bonds and show that it is present in a class of structures that are all lower in energy than any previously reported structure. Tracking this marker suggests that this motif is prominent in larger clusters as well as ∼180 nm ammonium bisulfate particles. Taken together, these results establish a set of structural motifs responsible for binding of gases at the surface of growing clusters that fully explain the spectrum of large particles and provide benchmarks for efforts to improve structure predictions, which are critical for the accurate theoretical treatment of this process.

Influence of atmospheric conditions on sulfuric acid-dimethylamine-ammonia-based new particle formation

A recent quantitative measurement of rates of new particle formation (NPF) in urban Shanghai showed that the high rates of NPF can be largely attributed to the sulfuric acid (SA)-dimethylamine (DMA) nucleation due to relatively high DMA concentration in urban atmosphere (Yao et al., Science. 2018, 361, 278). In certain atmospheric conditions, the release of DMA is accompanied with the emission of high concentration of ammonia. As a result, the ammonia (A) may participate in SA-DMA-based NPF. However, the main sources of DMA and A can be different, thereby leading to different mechanism for the SA-DMA-A-based nucleation under different atmospheric conditions. Near industrial sources with relatively high DMA concentration of 10⁸ molecules cm^-3, the contribution of binary SA-DMA nucleation to cluster formation is 61% at 278 K, representing a dominant pathway for NPF. However, in the region not too close to major source of DMA emission, e.g., near agriculture farmland, the routes involving ternary SA-DMA-A nucleation make a 64% contribution at 278 K with DMA concentration of 10⁷ molecules cm^-3, showing that A has marked impact on the cluster formation. Under such a condition, we predict that coexisting DMA and A could be detected in the process of NPF. Moreover, at winter temperatures or at higher altitudes, our calculations suggest that the clustering of initial clusters likely involve ternary SA-DMA-A clusters rather than binary SA-DMA clusters. These new insights may be helpful to analyze and predict atmospheric-condition-dependent NFP in either urban or rural regions and/or in different season of the year.

Integrated experimental and theoretical approach to probe the synergistic effect of ammonia in methanesulfonic acid reactions with small alkylamines

While new particle formation events have been observed worldwide, our fundamental understanding of the precursors remains uncertain. It has been previously shown that small alkylamines and ammonia (NH3) are key actors in sub-3 nm particle formation through reactions with acids such as sulfuric acid (H2SO4) and methanesulfonic acid (CH3S(O)(O)OH, MSA), and that water also plays a role. Because NH3 and amines co-exist in air, we carried out combined experimental and theoretical studies examining the influence of the addition of NH3 on particle formation from the reactions of MSA with methylamine (MA) and trimethylamine (TMA). Experiments were performed in a 1 m flow reactor at 1 atm and 296 K. Measurements using an ultrafine condensation particle counter (CPC) and a scanning mobility particle sizer (SMPS) show that new particle formation was systematically enhanced upon simultaneous addition of NH3 to the MSA + amine binary system, with the magnitude depending on the amine investigated. For the MSA + TMA reaction system, the addition of NH3 at ppb concentrations produced a much greater effect (i.e. order of magnitude more particles) than the addition of ∼12 000 ppm water (corresponding to ∼45-50% relative humidity). The effect of NH3 on the MSA + MA system, which is already very efficient in forming particles on its own, was present but modest. Calculations of energies, partial charges and structures of small cluster models of the multi-component particles likewise suggest synergistic effects due to NH3 in the presence of MSA and amine. The local minimum structures and the interactions involved suggest mechanisms for this effect.

Enhancing Potential of Trimethylamine Oxide on Atmospheric Particle Formation

The role of an oxidation product of trimethylamine, trimethylamine oxide, in atmospheric particle formation is studied using quantum chemical methods and cluster formation simulations. Molecular-level cluster formation mechanisms are resolved, and theoretical results on particle formation are confirmed with mass spectrometer measurements. Trimethylamine oxide is capable of forming only one hydrogen bond with sulfuric acid, but unlike amines, trimethylamine oxide can form stable clusters via ion–dipole interactions. That is because of its zwitterionic structure, which causes a high dipole moment. Cluster growth occurs close to the acid:base ratio of 1:1, which is the same as for other monoprotic bases. Enhancement potential of trimethylamine oxide in particle formation is much higher than that of dimethylamine, but lower compared to guanidine. Therefore, at relatively low concentrations and high temperatures, guanidine and trimethylamine oxide may dominate particle formation events over amines.

Piperazine Enhancing Sulfuric Acid-Based New Particle Formation: Implications for the Atmospheric Fate of Piperazine

Piperazine (PZ), a cyclic diamine, is one of 160 detected atmospheric amines and an alternative solvent to the widely used monoethanolamine in post-combustion CO₂ capture. Participating in H₂SO₄ (sulfuric acid, SA)-based new particle formation (NPF) could be an important removal pathway for PZ. Here, we employed quantum chemical calculations and kinetics modeling to evaluate the enhancing potential of PZ on SA-based NPF by examining the formation of PZ-SA clusters. The results indicate that PZ behaves more like a monoamine in stabilizing SA and can enhance SA-based NPF at the parts per trillion (ppt) level. The enhancing potential of PZ is less than that of the chainlike diamine putrescine and greater than that of dimethylamine, which is one of the strongest enhancing agents confirmed by ambient observations and experiments. After the initial formation of the (PZ)₁(SA)₁ cluster, the cluster mainly grows by gradual addition of SA or PZ monomer, followed by addition of (PZ)₁(SA)₁ cluster. We find that the ratio of PZ removal by NPF to that by the combination of NPF and oxidations is 0.5-0.97 at 278.15 K. As a result, we conclude that participation in the NPF pathway could significantly alter the environmental impact of PZ compared to only considering oxidation pathways.

Configurational Sampling of Noncovalent (Atmospheric) Molecular Clusters: Sulfuric Acid and Guanidine

We studied the configurational sampling of noncovalently bonded molecular clusters relevant to the atmosphere. In this article, we discuss possible approaches to searching for optimal configurations and present one alternative based on systematic configurational sampling, which seems able to overcome the typical problems associated with searching for global minima on multidimensional potential energy surfaces. Since atmospheric molecular clusters are usually held together by intermolecular bonds, we also present a cost-effective strategy for treating hydrogen bonding and proton transferring by using rigid molecules and ions in different protonation states and illustrate its performance on clusters containing guanidine and sulfuric acid.

Hydration motifs of ammonium bisulfate clusters of relevance to atmospheric new particle formation

We have analyzed the binding motifs of water bound to a prototypical cluster containing three ammonium cations and two bisulfate anions using mass-selective vibrational spectroscopy and quantum chemical calculations.

Direct Link between Structure and Hydration in Ammonium and Aminium Bisulfate Clusters Implicated in Atmospheric New Particle Formation

The acid-base chemistry of amines and sulfuric acid promotes growth in the early stages of atmospheric new particle formation, with more basic amines enhancing growth rates. Hydration of these particles has been proposed to depend on acidity or basicity but is difficult to quantify; therefore, the role of water in this process is not well understood. Using tandem mass spectrometry coupled to a temperature-controlled ion trap, we show that water uptake by aminium bisulfate clusters depends on the total number of free hydrogen bond donors in the cluster and is unaffected by the interchange of amines featuring the same number of substituents but differing gas-phase basicity. Analyzing this trend reveals site-specific propensities for hydration. These results indicate that hydration is determined by structural factors and that reported dependences on acidity or basicity arise from the weaker correlation between the number of hydrogen bond donors of amines and their gas-phase basicity.

Guanidine: A Highly Efficient Stabilizer in Atmospheric New-Particle Formation

The role of a strong organobase, guanidine, in sulfuric acid-driven new-particle formation is studied using state-of-the-art quantum chemical methods and molecular cluster formation simulations. Cluster formation mechanisms at the molecular level are resolved, and theoretical results on cluster stability are confirmed with mass spectrometer measurements. New-particle formation from guanidine and sulfuric acid molecules occurs without thermodynamic barriers under studied conditions, and clusters are growing close to a 1:1 composition of acid and base. Evaporation rates of the most stable clusters are extremely low, which can be explained by the proton transfers and symmetrical cluster structures. We compare the ability of guanidine and dimethylamine to enhance sulfuric acid-driven particle formation and show that more than 2000-fold concentration of dimethylamine is needed to yield as efficient particle formation as in the case of guanidine. At similar conditions, guanidine yields 8 orders of magnitude higher particle formation rates compared to dimethylamine. Highly basic compounds such as guanidine may explain experimentally observed particle formation events at low precursor vapor concentrations, whereas less basic and more abundant bases such as ammonia and amines are likely to explain measurements at high concentrations.

The Interplay Between Hydrogen Bonding and Coulombic Forces in Determining the Structure of Sulfuric Acid-Amine Clusters

Acid-base cluster chemistry drives atmospheric new particle formation (NPF), but the details of the growth mechanisms are difficult to experimentally probe. Clusters of ammonia, alkylamines, and sulfuric acid, species fundamental to NPF, are probed by infrared spectroscopy. These spectra show that substitution of amines for ammonia, which is linked to accelerated growth, induces profound structural rearrangement in clusters with initial compositions (NH₄⁺) _n+1(HSO₄^-) _n (1 ≤ n ≤ 3). This rearrangement is driven by the loss of N-H hydrogen bond donors, yielding direct bisulfate-bisulfate hydrogen bonds, and its onset with respect to cluster composition indicates that more substituted amines induce rearrangement at smaller sizes. A simple model counting hydrogen bond donors and acceptors explains these observations. The presence of direct hydrogen bonds between formal anions shows that hydrogen bonding can compete with Coulombic forces in determining cluster structure. These results suggest that NPF mechanisms may be highly dependent on amine identity.

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.

Interaction of oxalic acid with methylamine and its atmospheric implications

Oxalic acid, which is one of the most common dicarboxylic acids, is expected to be an important component of atmospheric aerosols. However, the contribution of oxalic acid to the generation of new particles is still poorly understood. In this study, the structural characteristics and thermodynamics of (C2H2O4)(CH3NH2) n (n = 1-4) were investigated at the PW91PW91/6-311++G(3df,3pd) level of theory. We found that clusters formed by oxalic acid and methylamine are relatively stable, and the more the atoms participating in the formation of a ring-like structure, the more stable is the cluster. In addition, via the analysis of atmospheric relevance, it can be revealed that clusters of (C2H2O4)(CH3NH2) n (n = 1-4) have a noteworthy concentration in the atmosphere, which indicates that these clusters could be participating in new particle formation. Moreover, by comparison with (H2C2O4)(NH3) n (n = 1-6) species, it can be seen that oxalic acid is more readily bound to methylamine than to ammonia, which promotes nucleation or new particle formation. Finally, the Rayleigh scattering properties of clusters of (C2H2O4)(CH3NH2) n (n = 1-4) were investigated for the first time to determine their atmospheric implications.