Iodine Clusters in the Atmosphere I: Computational Benchmark and Dimer Formation of Oxyacids and Oxides

The contribution of iodine-containing compounds to atmospheric new particle formation is still not fully understood, but iodic acid and iodous acid are thought to be significant contributors. While several quantum chemical studies have been carried out on clusters containing iodine, there is no comprehensive benchmark study quantifying the accuracy of the applied methods. Here, we present the first study in a series that investigate the role of iodine species in atmospheric cluster formation. In this work, we have studied the iodic acid, iodous acid, iodine tetroxide, and iodine pentoxide monomers and their dimers formed with common atmospheric precursors. We have tested the accuracy of commonly applied methods for calculating the geometry of the monomers, thermal corrections of monomers and dimers, the contribution of spin-orbit coupling to monomers and dimers, and finally, the accuracy of the electronic energy correction calculated at different levels of theory. We find that optimizing the structures either at the ωB97X-D3BJ/aug-cc-pVTZ-PP or the M06-2X/aug-cc-pVTZ-PP level achieves the best thermal contribution to the binding free energy. The electronic energy correction can then be calculated at the ZORA-DLPNO-CCSD(T0) level with the SARC-ZORA-TZVPP basis for iodine and ma-ZORA-def2-TZVPP for non-iodine atoms. We applied this methodology to calculate the binding free energies of iodine-containing dimer clusters, where we confirm the qualitative trends observed in previous studies. However, we identify that previous studies overestimate the stability of the clusters by several kcal/mol due to the neglect of relativistic effects. This means that their contributions to the currently studied nucleation pathways of new particle formation are likely overestimated.

Natural Marine Precursors Boost Continental New Particle Formation and Production of Cloud Condensation Nuclei

Marine dimethyl sulfide (DMS) emissions are the dominant source of natural sulfur in the atmosphere. DMS oxidizes to produce low-volatility acids that potentially nucleate to form particles that may grow into climatically important cloud condensation nuclei (CCN). In this work, we utilize the chemistry transport model ADCHEM to demonstrate that DMS emissions are likely to contribute to the majority of CCN during the biological active period (May-August) at three different forest stations in the Nordic countries. DMS increases CCN concentrations by forming nucleation and Aitken mode particles over the ocean and land, which eventually grow into the accumulation mode by condensation of low-volatility organic compounds from continental vegetation. Our findings provide a new understanding of the exchange of marine precursors between the ocean and land, highlighting their influence as one of the dominant sources of CCN particles over the boreal forest.

Effect of an Electric Field on the Structure and Stability of Atmospheric Clusters

We study the influence of an applied electric field on the structure and stability of some common bimolecular clusters that are found in the atmosphere. These clusters play an important role in new particle formation (NPF). For low values of the electric field (i.e., |E| ≤ 0.01 V Å-1), we demonstrate that the field response of the clusters can be predicted from simply calculating the dipole moment of the cluster and the dipole moments of the constituent molecules and that the influence on the association energy of the cluster is minimal (i.e., <0.5 kcal mol-1). For higher field strengths |E| > 0.2 V Å-1, there can be more dramatic effects on both structure and energetics, as the induced dipole, charge transfer, and geometric distortion play a larger role. Although such large fields are not very relevant in the atmosphere, they do exist in some situations of experimental interest, such as near interfaces and in intense laser fields.

Basis Set Superposition Errors Are Partly Basis Set Imbalances

When calculating fragment interaction energies by electronic structure methods employing medium-sized atom-centered basis sets, it is often observed that the effect is systematically overestimated. The common interpretation is that the systematic error arises because the basis set for the complex is more complete than for the isolated fragments, and this is denoted basis set superposition errors. It has been observed, however, that the interaction energy in some cases is underestimated, which defies the interpretation in terms of basis set completeness, and instead suggests that the effect partly is due to basis set imbalance. The imbalance can be removed by explicit optimization of the basis sets for each structure, and it is shown that this to a significant extent reduces the systematic overestimation attributed to basis set superposition error.

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.

Significant influence of water molecules on the SO3 + HCl reaction in the gas phase and at the air-water interface

The products resulting from the reactions between atmospheric acids and SO3 have a catalytic effect on the formation of new particles in aerosols. However, the SO3 + HCl reaction in the gas-phase and at the air-water interface has not been considered. Herein, this reaction was explored exhaustively by using high-level quantum chemical calculations and Born Oppenheimer molecular dynamics (BOMD) simulations. The quantum calculations show that the gas-phase reaction of SO3 + HCl is highly unlikely to occur under atmospheric conditions with a high energy barrier of 22.6 kcal mol-1. H2O and (H2O)2 play obvious catalytic roles in reducing the energy barrier of the SO3 + HCl reaction by over 18.2 kcal mol-1. The atmospheric lifetimes of SO3 show that the (H2O)2-assisted reaction dominates over the H2O-assisted reaction within the altitude range of 0-5 km, whereas the H2O-assisted reaction is more favorable within an altitude range of 10-50 km. BOMD simulations show that H2O-induced formation of the ClSO3-⋯H3O+ ion pair and HCl-assisted formation of the HSO4-⋯H3O+ ion pair were identified at the air-water interface. These routes followed a stepwise reaction mechanism and proceeded at a picosecond time scale. Interestingly, the formed ClSO3H in the gas phase has a tendency to aggregate with sulfuric acids, ammonias, and water molecules to form stable clusters within 40 ns simulation time, while the interfacial ClSO3- and H3O+ can attract H2SO4, NH3, and HNO3 for particle formation from the gas phase to the water surface. Thus, this work will not only help in understanding the SO3 + HCl reaction driven by water molecules in the gas-phase and at the air-water interface, but it will also provide some potential routes of aerosol formation from the reaction between SO3 and inorganic acids.

Toward Modeling the Growth of Large Atmospheric Sulfuric Acid-Ammonia Clusters

Studying large atmospheric molecular clusters is needed to understand the transition between clusters and aerosol particles. In this work, we studied the (SA)n(AM)n clusters with n up to 30 and the (SA)m(AM)m±2 clusters, with m = 6-20. The cluster configurations are sampled using the ABCluster program, and the cluster geometries and thermochemical parameters are calculated using GFN1-xTB. The cluster binding energies are calculated using B97-3c. We find that the addition of sulfuric acid is preferred to the addition of ammonia. The addition free energies were found to have large uncertainties, which could potentially be attributed to errors in the applied level of theory. Based on DLPNO-CCSD(T0)/aug-cc-pVTZ benchmarks of the binding energies of the large (SA)8-9(AM)10 and (SA)10(AM)10-11 clusters, we find that ωB97X-D3BJ with a large basis set is required to yield accurate binding and addition energies. However, based on recalculations of the single-point energy at r2SCAN-3c and ωB97X-D3BJ/6-311++G(3df,3pd), we show that the single-point energy contribution is not the primary source of error. We hypothesize that a larger source of error might be present in the form of insufficient configurational sampling. Finally, we train Δ machine learning model on (SA)n(AM)n clusters with n up to 5 and show that we can predict the binding energies of clusters up to sizes of (SA)30(AM)30 with a binding energy error below 0.6 %. This is an encouraging approach for accurately modeling the binding energies of large acid-base clusters in the future.

Nitric Acid and Organic Acids Suppress the Role of Methanesulfonic Acid in Atmospheric New Particle Formation

Multicomponent atmospheric molecular clusters, typically comprising a combination of acids and bases, play a pivotal role in our climate system and contribute to the perplexing uncertainties embedded in modern climate models. Our understanding of cluster formation is limited by the lack of studies on complex mixed-acid-mixed-base systems. Here, we investigate multicomponent clusters consisting of mixtures of several acid and base molecules: sulfuric acid (SA), methanesulfonic acid (MSA), nitric acid (NA), formic acid (FA), along with methylamine (MA), dimethylamine (DMA), and trimethylamine (TMA). We calculated the binding free energies of a comprehensive set of 252 mixed-acid-mixed-base clusters at the DLPNO-CCSD(T0)/aug-cc-pVTZ//ωB97X-D/6-31++G(d,p) level of theory. Combined with the existing datasets, we simulated the new particle formation (NPF) rates using the Atmospheric Cluster Dynamics Code (ACDC). We find that the presence of NA and FA had a substantial impact, increasing the NPF rate by 60% at realistic conditions. Intriguingly, we find that NA and FA suppress the role of MSA in NPF. These findings suggest that even high concentration of MSA has a limited impact on NPF in polluted regions with high FA and NA. We outline a method for generating a lookup table that could potentially be used in climate models that sufficiently incorporates all the required chemistry. By unraveling the molecular mechanisms of mixed-acid-mixed-base clusters, we get one step closer to comprehending their implications for our global climate system.

Quantifying Intramolecular Basis Set Superposition Errors

We show that medium-sized Gaussian basis sets lead to significant intramolecular basis set superposition errors at Hartree-Fock and density functional levels of theory, with artificial stabilization of compact over extended conformations for a 186 atom deca-peptide. Errors of ∼80 and ∼10 kJ/mol are observed, with polarized double zeta and polarized triple zeta quality basis sets, respectively. Two different procedures for taking the basis set superposition error into account are tested. While both reduce the error, it appears that polarized quadruple zeta basis sets are required to reduce the error below a few kJ/mol. Alternatively, the basis set superposition error can be eliminated using multiresolution methods based on Multiwavelets.

Methanesulfonic acid and iodous acid nucleation: a novel mechanism for marine aerosols

By seeding clouds, new particle formation (NPF) has a substantial impact on radiation balance, bio-geochemical cycles and global climate. Over oceans, both methanesulfonic acid (CH₃S(O)₂OH, MSA) and iodous acid (HIO₂) have been reported to be closely associated with NPF events; however, much less is known about whether they can jointly nucleate to form nanoclusters. Hence, quantum chemical calculations and Atmospheric Cluster Dynamics Code (ACDC) simulations were performed to investigate the novel mechanism of MSA-HIO₂ binary nucleation. The results indicate that MSA and HIO₂ can form stable clusters via multiple interactions including hydrogen bonds, halogen bonds, and electrostatic forces between ion pairs after proton transfer, which are more diverse than those in MSA-iodic acid (HIO₃) and MSA-dimethylamine (DMA) clusters. Interestingly, HIO₂ can be protonated by MSA exhibiting base-like behavior, but it differs from base nucleation precursors by self-nucleation rather than solely binding to MSA. Due to the greater stability of MSA-HIO₂ clusters, the formation rate of MSA-HIO₂ clusters can be even higher than that of MSA-DMA clusters, suggesting that MSA-HIO₂ nucleation is a non-negligible source of marine NPF. This work proposes a novel mechanism of MSA-HIO₂ binary nucleation for marine aerosols and provides deeper insights into the distinctive nucleation characteristics of HIO₂, which can help in constructing a more comprehensive sulfur- and iodine-bearing nucleation model for marine NPF.

Chemical Implications of Rapid Reactive Absorption of I2O4 at the Air-Water Interface

Marine aerosol formation involving iodine-bearing species significantly affects the global climate and radiation balance. Although recent studies outline the critical role of iodine oxide in nucleation, much less is known about its contribution to aerosol growth. This paper presents molecular-level evidence that the air-water interfacial reaction of I₂O₄ mediated by potent atmospheric chemicals, such as sulfuric acid (H₂SO₄) and amines [e.g., dimethylamine (DMA) and trimethylamine (TMA)], can occur rapidly on a picosecond time scale by Born-Oppenheimer molecular dynamics simulations. The interfacial water bridges the reactants while facilitating the DMA-mediated proton transfer and stabilizing the ionic products of H₂SO₄-involved reactions. The identified heterogeneous mechanisms exhibit the dual contribution to aerosol growth: (i) the ionic products (e.g., IO₃^-, DMAH⁺, TMAH⁺, and HSO₄^-) formed by reactive adsorption possess less volatility than the reactants and (ii) these ions, such as alkylammonium salts (e.g., DMAH⁺), are also highly hydrophilic, further facilitating hygroscopic growth. This investigation enhances not only our understanding of heterogeneous iodine chemistry but also the impact of iodine oxide on aerosol growth. Also, these findings can bridge the gap between the abundance of I₂O₄ in the laboratory and its absence in field-collected aerosols and provide an explanation for the missing source of IO₃^-, HSO₄^-, and DMAH⁺ in marine aerosols.

Heterogeneous Ion-Induced Nucleation of Water and Butanol Vapors Studied via Computational Quantum Chemistry beyond Prenucleation and Critical Cluster Sizes

Water and butanol are used as working fluids in condensation particle counters, and condensation of a single vapor onto an ion can be used as a simple model system for the study of ion-induced nucleation in the atmosphere. Motivated by this, we examine heterogeneous nucleation of water (H₂O) and n-butanol (BuOH) vapors onto three positively (Li⁺, Na⁺, K⁺) and three negatively charged (F^-, Cl^-, Br^-) ions using classical nucleation theory and computational quantum chemistry methods. We study phenomena that cannot be captured by Kelvin-Thomson equation for small nucleation ion cores. Our quantum chemistry calculations reveal the molecular mechanism behind ion-induced nucleation for each studied system. Typically, ions become solvated from all sides after several vapor molecules condense onto the ion. However, we show that the clusters of water and large negatively charged ions (Cl^- and Br^-) thermodynamically prefer the ion being migrated to the cluster surface. Although our methods generally do not show clear sign-preference for ion-water nucleation, we identified positive sign-preference for ion-butanol nucleation caused by the possibility to form stabilizing hydrogen bonds between butanol molecules condensed onto a positively charged ion. These bonds cannot form when butanol condenses onto a negatively charged ion. Therefore, we show that ion charge, its sign, as well as vapor properties have effects on the prenucleation and critical cluster/droplet sizes and also on the molecular mechanism of ion-induced nucleation.

Clusteromics V: Organic Enhanced Atmospheric Cluster Formation

Formic acid (FA) is a prominent candidate for organic enhanced nucleation due to its high abundance and stabilizing effect on smaller clusters. Its role in new particle formation is studied through the use of state-of-the-art quantum chemical methods on the cluster systems (acid)_1-2(FA)₁(base)_1-2 with the acids being sulfuric acid (SA)/methanesulfonic acid (MSA) and the bases consisting of ammonia (A), methylamine (MA), dimethylamine (DMA), trimethylamine (TMA), and ethylenediamine (EDA). A funneling approach is used to determine the cluster structures with initial configurations generated through the ABCluster program, followed by semiempirical PM7 and ωB97X-D/6-31++G(d,p) calculations. The final binding free energy is calculated at the DLPNO-CCSD(T₀)/aug-cc-pVTZ//ωB97X-D/6-31++G(d,p) level of theory using the quasi-harmonic approximation. Cluster dynamics simulations show that FA has a minuscule or negligible effect on the MSA-FA-base systems as well as most of the SA-FA-base systems. The SA-FA-DMA cluster system shows the highest influence from FA with an enhancement of 21%, compared to its non-FA counterpart.

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.

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.

Amine-Enhanced Methanesulfonic Acid-Driven Nucleation: Predictive Model and Cluster Formation Mechanism

Atmospheric amines are considered to be an effective enhancer for methanesulfonic acid (MSA)-driven nucleation. However, out of the 195 detected atmospheric amines, the enhancing potential (EP) has so far only been studied for five amines. This severely hinders the understanding of the contribution of amines to MSA-driven nucleation. Herein, a two-step procedure was employed to probe the EP of various amines on MSA-driven nucleation. Initially, the formation free energies (ΔG) of 50 MSA-amine dimer clusters were calculated. Based on the calculated ΔG values, a robust quantitative structure-activity relationship (QSAR) model was built and utilized to predict the ΔG values of the remaining 145 amines. The QSAR model identified two guanidino-containing compounds as the potentially strongest enhancer for MSA-driven nucleation. Second, the EP of guanidino-containing compounds was studied by employing larger clusters and selecting guanidine (Gud) as a representative. The results indicate that Gud indeed has the strongest EP. The Gud-MSA system presents a unique clustering mechanism, proceeding via the initial formation of the (Gud)₁(MSA)₁ cluster, and subsequently by cluster collisions with either a (Gud)₁(MSA)₁ or (Gud)₂(MSA)₂ cluster. The developed QSAR model and the identification of amines with the strongest EP provide a foundation for comprehensively evaluating the contribution of atmospheric amines to MSA-driven nucleation.

Iodous acid - a more efficient nucleation precursor than iodic acid

Iodous acid (HIO₂), a vital iodine oxyacid, potentially plays an important role in the formation of new particles in marine areas (He et al., Science, 2021, 371, 589-595). However, the nucleation mechanism of HIO₂ is still poorly understood. Herein, the self-nucleation of HIO₂ under different atmospheric conditions is investigated by a combination of quantum chemical calculations and the Atmospheric Cluster Dynamics Code (ACDC) simulations. The results indicate that HIO₂ can form relatively stable molecular clusters through hydrogen bonds and halogen bonds, and the self-nucleation of HIO₂ proceeds by sequential addition of HIO₂ or HIO₂-based small clusters. Besides, in order to better illustrate the role of HIO₂ in new particle formation (NPF) in marine areas, we compare its nucleation properties with those of iodic acid (HIO₃), a significant iodine-containing nucleation precursor in marine regions. We find that the cluster formation rate of the self-nucleation of HIO₂ is higher than that of the self-nucleation of HIO₃ although [HIO₂] is lower than [HIO₃], which indicates that the HIO₂ molecule is a more efficient nucleation precursor than the HIO₃ molecule. Therefore, the self-nucleation of HIO₂ could become one of the most important sources for NPF in marine areas, which could provide potential theoretical evidence for explaining the intensive NPF events observed in these areas.

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.

Gas-phase catalytic hydration of I2O5 in the polluted coastal regions: Reaction mechanisms and atmospheric implications

Marine aerosols play an important role in the global aerosol system. In polluted coastal regions, ultra-fine particles have been recognized to be related to iodine-containing species and is more serious due to the impact of atmospheric pollutants. Many previous studies have identified iodine pentoxide (I₂O₅, IP) to be the key species in new particles formation (NPF) in marine regions, but the role of IP in the polluted coastal atmosphere is far to be fully understood. Considering the high humidity and concentrations of pollutants in the polluted coastal regions, the gas-phase hydration of IP catalyzed by sulfuric acid (SA), nitric acid (NA), dimethylamine (DMA), and ammonia (A) have been investigated at DLPNO-CCSD(T)//ωB97X-D/aug-cc-pVTZ + aug-cc-pVTZ-PP with ECP28MDF (for iodine) level of theory. The results show that the hydration of IP involves a significant energy barrier of 22.33 kcal/mol, while the pollutants SA, NA, DMA, and A all could catalyze the hydration of IP. Especially, with SA and DMA as catalysts, the hydration reactions of IP present extremely low barriers and high rate constants. It is suggested that IP is unstable under the catalysis of SA and DMA to generate iodic acid, which is the key component in NPF in marine regions. Thus, the catalytic hydration of IP is very likely to trigger the formation of iodine-containing particles. Our research provides a clear picture of the catalytic hydration of IP as well as theoretical guidance for NPF in the polluted coastal atmosphere.

DFT Study of the Formation of Atmospheric Aerosol Precursors from the Interaction between Sulfuric Acid and Benzenedicarboxylic Acid Molecules

Dicarboxylic acids are ubiquitous products of the photooxidation of volatile organic compounds which are believed to play a significant role in the formation of secondary organic aerosols in the atmosphere. In this paper, we report high-level quantum investigations of the clustering properties of sulfuric acid and benzenedicarboxylic acid molecules. Up to four molecules have been considered in the calculations, and the behavior of the three isomers of the organic diacid species have been compared. The most stable geometries have been characterized together with the corresponding thermodynamic data. From an atmospheric point of view, the results of the DFT calculations show that the organic diacid molecules may significantly enhance the nucleation of small atmospheric clusters, at least from an energetic point of view. In this respect, the phthalic acid isomer seems more efficient than the two other isomers of the benzenedicarboxylic acid, in particular because the internal distance between the two carboxyl groups in the organic diacids appears to play an important role in the stabilization of the H-bond network inside the corresponding heterocluster formed with sulfuric acid molecules.

Real-time monitoring of aerosol particle formation from sulfuric acid vapor at elevated concentrations and temperatures

In the present study, time-resolved aerosol particle formation from sulfuric acid vapor is examined with special attention to the stabilization of molecular clusters in the early phase of unary nucleation. An important factor governing this process is the amount of condensable acid vapor. Here it is produced from fast gas-phase reactions in a batch-type reaction cell for which we introduce modifications enabling real-time monitoring. The key component for size- and time-resolved detection of ultrafine particles is a new 1 nm-SMPS. With this new tool at hand, the effect of varying the precursor concentration over two orders of magnitude is investigated. We demonstrate the ability to tune between different growth scenarios as indicated by the size-resolved particle traces which exhibit a transition from sigmoidal over quasi-stationary to peak-like shape. The second key parameter relevant for nucleation studies is the temperature-dependent cluster evaporation. Due to a temperature rise during the mixing stage of the experiment, evaporation is strongly promoted in the early phase. Therefore, the present study extends the T-range used in, e.g., smog chambers. We investigate this temperature effect in a kinetic simulation and can successfully combine simulated and measured data for validating theoretical evaporation rates obtained from DLPNO-CCSD(T₀)-calculations.

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.

New Particle Formation and Growth from Dimethyl Sulfide Oxidation by Hydroxyl Radicals

Dimethyl sulfide (DMS) is produced by plankton in oceans and constitutes the largest natural emission of sulfur to the atmosphere. In this work, we examine new particle formation from the primary pathway of oxidation of gas-phase DMS by OH radicals. We particularly focus on particle growth and mass yield as studied experimentally under dry conditions using the atmospheric simulation chamber AURA. Experimentally, we show that aerosol mass yields from oxidation of 50-200 ppb of DMS are low (2-7%) and that particle growth rates (8.2-24.4 nm/h) are comparable with ambient observations. An HR-ToF-AMS was calibrated using methanesulfonic acid (MSA) to account for fragments distributed across both the organic and sulfate fragmentation table. AMS-derived chemical compositions revealed that MSA was always more dominant than sulfate in the secondary aerosols formed. Modeling using the Aerosol Dynamics, gas- and particle-phase chemistry kinetic multilayer model for laboratory CHAMber studies (ADCHAM) indicates that the Master Chemical Mechanism gas-phase chemistry alone underestimates experimentally observed particle formation and that DMS multiphase and autoxidation chemistry is needed to explain observations. Based on quantum chemical calculations, we conclude that particle formation from DMS oxidation in the ambient atmosphere will most likely be driven by mixed sulfuric acid/MSA clusters clustering with both amines and ammonia.

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.

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.

Bond Dissociation Energies in the Gas Phase for Large Molecular Ions by Threshold Collision-Induced Dissociation Experiments: Stretching the Limits

Accurate bond dissociation energies for large molecules are difficult to obtain by either experimental or computational methods. The former methods are hampered by a range of physical and practical limitations in gas-phase measurement techniques, while the latter require incorporation of multiple approximations whose impact on accuracy may not always be clear. When internal benchmarks are not available, one hopes that experiment and theory can mutually support each other. A recent report found, however, a large discrepancy between gas-phase bond dissociation energies, measured mass spectrometrically, and the corresponding quantities computed using density functional theory (DFT)-D3 and DLPNO-CCSD(T) methods. With the widespread application of these computational methods to large molecular systems, the discrepancy needs to be resolved. We report a series of experimental studies that validate the mass spectrometric methods from small to large ions and find that bond dissociation energies extracted from threshold collision-induced dissociation experiments on large ions do indeed behave correctly. The implications for the computational studies are discussed.

Neutral Sulfuric Acid-Water Clustering Rates: Bridging the Gap between Molecular Simulation and Experiment

The role of sulfuric acid during atmospheric new particle formation is an ongoing topic of discussion. In this work, we provide quantitative experimental constraints for quantum chemically calculated evaporation rates for the smallest H₂SO₄-H₂O clusters, characterizing the mechanism governing nucleation on a kinetic, single-molecule level. We compare experimental particle size distributions resulting from a highly supersaturated homogeneous H₂SO₄ gas phase with the results from kinetic simulations employing quantum chemically derived decomposition rates of electrically neutral H₂SO₄ molecular clusters up to the pentamer at a large range of relative humidities. By using high H₂SO₄ concentrations, we circumvent the uncertainties concerning contaminants and competing reactions present in studies at atmospheric conditions. We show good agreement between molecular simulation and experimental measurements and provide the first evaluation of theoretical predictions of the stabilization provided by water molecules.

Assessment of the DLPNO Binding Energies of Strongly Noncovalent Bonded Atmospheric Molecular Clusters

This work assesses the performance of DLPNO-CCSD(T0), DLPNO-MP2, and density functional theory methods in calculating the binding energies of a representative test set of 45 atmospheric acid-acid, acid-base, and acid-water dimer clusters. The performance of the approximate methods is compared to high level explicitly correlated CCSD(F12*)(T)/complete basis set (CBS) reference calculations. Out of the tested density functionals, ωB97X-D3(BJ) shows the best performance with a mean deviation of 0.09 kcal/mol and a maximum deviation of 0.83 kcal/mol. The RI-CC2/aug-cc-pV(T+d)Z level of theory severely overpredicts the cluster binding energies with a mean deviation of -1.31 kcal/mol and a maximum deviation up to -3.00 kcal/mol. Hence, RI-CC2/aug-cc-pV(T+d)Z should not be utilized for studying atmospheric molecular clusters. The DLPNO variants are tested both with and without the inclusion of explicit correlation (F12) in the wavefunction, with different pair natural orbital (PNO) settings (loosePNO, normalPNO, and tightPNO) and using both double and triple zeta basis sets. The performance of the DLPNO-MP2 methods is found to be independent of PNO settings and yield low mean deviations of -0.84 kcal/mol or below. However, DLPNO-MP2 requires explicitly correlated wavefunctions to yield maximum deviations below 1.40 kcal/mol. For obtaining high accuracy, with maximum deviation below ∼1.0 kcal/mol, either DLPNO-CCSD(T0)/aug-cc-pVTZ (normalPNO) calculations or DLPNO-CCSD(T0)-F12/cc-pVTZ-F12 (normalPNO) calculations are required. The most accurate level of theory is found to be DLPNO-CCSD(T0)-F12/cc-pVTZ-F12 using a tightPNO criterion which yields a mean deviation of 0.10 kcal/mol, with a maximum deviation of 0.20 kcal/mol, compared to the CCSD(F12*)(T)/CBS reference.

Role of glycine on sulfuric acid-ammonia clusters formation: Transporter or participator

Glycine (Gly) is ubiquitous in the atmosphere and plays a vital role in new particle formation (NPF). However, the potential mechanism of its on sulfuric acid (SA) - ammonia (A) clusters formation under various atmospheric conditions is still ambiguous. Herein, a (Gly)_x·(SA)_y·(A)_z (z ≤ x + y ≤ 3) multicomponent system was investigated by using density functional theory (DFT) combined with Atmospheric Cluster Dynamics Code (ACDC) at different temperatures and precursor concentrations. The results show that Gly, with one carboxyl (-COOH) and one amine (-NH₂) group, can interact strongly with SA and A in two directions through hydrogen bonds or proton transfer. Within the relevant range of atmospheric concentrations, Gly can enhance the formation rate of SA-A-based clusters, especially at low temperature, low [SA], and median [A]. The enhancement (R) of Gly on NPF can be up to 340 at T = 218.15 K, [SA] = 10⁴, [A] = 10⁹, and [Gly] = 10⁷ molecules/cm³. In addition, the main growth paths of clusters show that Gly molecules participate into cluster formation in the initial stage and eventually leave the cluster by evaporation in subsequent cluster growth at low [Gly], it acts as an important "transporter" to connect the smaller and larger cluster. With the increase of [Gly], it acts as a "participator" directly participating in NPF.

Effect of NH3 and HCOOH on the H2O2 + HO → HO2 + H2O reaction in the troposphere: competition between the one-step and stepwise mechanisms

The H₂O₂ + HO → HO₂ + H₂O reaction is an important reservoir for both radicals of HO and HO₂ catalyzing the destruction of O₃. Here, this reaction assisted by NH₃ and HCOOH catalysts was explored using the CCSD(T)-F12a/cc-pVDZ-F12//M06-2X/aug-cc-pVTZ method and canonical variational transition state theory with small curvature tunneling. Two possible sets of mechanisms, (i) one-step routes and (ii) stepwise processes, are possible. Our results show that in the presence of both NH₃ and HCOOH catalysts under relevant atmospheric temperature, mechanism (i) is favored both energetically and kinetically than the corresponding mechanism (ii). At 298 K, the relative rate for mechanism (i) in the presence of NH₃ (10, 2900 ppbv) and HCOOH (10 ppbv) is respectively 3–5 and 2–4 orders of magnitude lower than that of the water-catalyzed reaction. This is due to a comparatively lower concentration of NH₃ and HCOOH than H₂O which indicates the positive water effect under atmospheric conditions. Although NH₃ and HCOOH catalysts play a negligible role in the reservoir for both radicals of HO and HO₂ catalyzing the destruction of O₃, the current study provides a comprehensive example of how acidic and basic catalysts assisted the gas-phase reactions.

The H₂O₂ + HO → HO₂ + H₂O reaction is an important reservoir for both radicals of HO and HO₂ catalyzing the destruction of O₃.

Methanesulfonic Acid-driven New Particle Formation Enhanced by Monoethanolamine: A Computational Study

Amines are recognized as significant enhancing species on methanesulfonic acid (MSA)-driven new particle formation (NPF). Monoethanolamine (MEA) has been detected in the atmosphere, and its concentration could be significantly increased once MEA-based postcombustion CO₂ capture technology is widely implemented. Here, we evaluated the enhancing potential of MEA on MSA-driven NPF by examining the formation of MEA-MSA clusters using a combination of quantum chemical calculations and kinetics modeling. The results indicate that the -OH group of MEA can form at least one hydrogen bond with MSA or MEA in all MEA-containing clusters. The enhancing potential of MEA is higher than that of the strongest enhancing agent known so far, methylamine (MA), for MSA-driven NPF. Such high enhancing potential can be ascribed to not only the higher gas-phase basicity but also the role of the additional -OH group of MEA in increasing the binding free energy by forming additional hydrogen bonds. This clarifies the importance of hydrogen-bonding capacity from the nonamino group of amines in enhancing MSA-driven NPF. The main growth pathway for MEA-MSA clusters proceeds via the initial formation of the (MEA)₁(MSA)₁ cluster, followed by alternately adding one MSA and one MEA molecule, differing from the case of MA-MSA clusters.

Unexpected quenching effect on new particle formation from the atmospheric reaction of methanol with SO3

Despite the high abundance in the atmosphere, alcohols in general and methanol in particular are believed to play a small role in atmospheric new particle formation (NPF) largely due to the weak binding abilities of alcohols with the major nucleation precursors, e.g., sulfuric acid (SA) and dimethylamine (DMA). Herein, we identify a catalytic reaction that was previously overlooked, namely, the reaction between methanol and SO₃, catalyzed by SA, DMA, or water. We found that alcohols can have unexpected quenching effects on the NPF process, particularly in dry and highly polluted regions with high concentrations of alcohols. Specifically, the catalytic reaction between methanol and SO₃ can convert methanol into a less-volatile species--methyl hydrogen sulfate (MHS). The latter was initially thought to be a good nucleation agent for NPF. However, our simulation results suggest that the formation of MHS consumes an appreciable amount of atmospheric SO₃, disfavoring further reactions of SO₃ with H₂O. Indeed, we found that MHS formation can cause a reduction of SA concentration up to 87%, whereas the nucleation ability of MHS toward new particles is not as good as that of SA. Hence, a high abundance of methanol in the atmosphere can lower the particle nucleation rate by as much as two orders of magnitude. Such a quenching effect suggests that the recently identified catalytic reactions between alcohols and SO₃ need to be considered in atmospheric modeling in order to predict SA concentration from SO₂, while also account for their potentially negative effect on NPF.

Strong Even/Odd Pattern in the Computed Gas-Phase Stability of Dicarboxylic Acid Dimers: Implications for Condensation Thermodynamics

The physical properties of small straight-chain dicarboxylic acids are well known to exhibit even/odd alternations with respect to the carbon chain length. For example, odd numbered diacids have lower melting points and higher saturation vapor pressures than adjacent even numbered diacids. This alternation has previously been explained in terms of solid-state properties, such as higher torsional strain of odd number diacids. Using quantum chemical methods, we demonstrate an additional contribution to this alternation in properties resulting from gas-phase dimer formation. Due to a combination of hydrogen bond strength and torsional strain, dimer formation in the gas phase occurs efficiently for glutaric acid (C5) and pimelic acid (C7) but is unfavorable for succinic acid (C4) and adipic acid (C6). Our results indicate that a significant fraction of the total atmospheric gas-phase concentration of glutaric and pimelic acid may consist of dimers.

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.

Enhancement of Atmospheric Nucleation by Highly Oxygenated Organic Molecules: A Density Functional Theory Study

New particle formation (NPF) by gas-particle conversion is the main source of atmospheric aerosols. Highly oxygenated organic molecules (HOMs) and sulfuric acid (SA) are important NPF participants. 2-Methylglyceric acid (MGA), a kind of HOMs, is a tracer of isoprene-derived secondary organic aerosols. The nucleation mechanisms of MGA with SA were studied using density functional theory and atmospheric cluster dynamics simulation in this study, along with that of MGA with methanesulfonic acid (MSA) as a comparison. Our theoretical works indicate that the (MGA)(SA) and (MGA)(MSA) clusters are the most stable ones in the (MGA) _i(SA) _j ( i = 1-2, j = 1-2) and (MGA) _i(MSA) _j ( i = 1-2, j = 1-2) clusters, respectively. Both the formation rates of (MGA)(SA) and (MGA)(MSA) clusters are quite large and could have significant contributions to NPF. The results imply that the homomolecular nucleation of MGA is unlikely to occur in the atmosphere, and MGA and SA can effectively contribute to heteromolecular nucleation mainly in the form of heterodimers. MSA exhibits properties similar to SA in its ability to form clusters with MGA but is slightly weaker than SA.

The reaction of hydrated iodide I(H2O)− with ozone: a new route to IO2− products

We report on an experimental characterization of the isolated reaction of hydrated iodide I(H₂O)⁻ with ozone O₃ at room temperature performed using a radio-frequency ion trap combined with a quadrupole mass spectrometer.

Effects of water, ammonia and formic acid on HO2 + Cl reactions under atmospheric conditions: competition between a stepwise route and one elementary step

Quantum chemical calculations at M06-2X and CCSD(T) levels of theory have been performed to investigate the effects of H₂O, NH₃, and HCOOH on the HO₂ + Cl → HCl + O₂ reaction. The results show that catalyzed reactions with three catalysts could proceed through two different mechanisms, namely a stepwise route and one elementary step, where the former reaction is more favorable than the latter. Meanwhile, for the stepwise route, a single hydrogen atom transfer pathway in the presence of all catalysts has more advantages than the respective double hydrogen atom transfer pathway. Then, the relative impacts of catalysts under tropospheric conditions were investigated by considering the temperature dependence of the rate constants and the altitude dependence of catalyst concentrations. The calculated results show that at 0 km altitude, the HO₂ + Cl → HCl + O₂ reaction with catalysts, such as H₂O, NH₃, or HCOOH, cannot compete with the reaction without a catalyst, as the effective rate constant with a catalyst is smaller by 2–6 orders of magnitude than the naked reaction within the temperature range 280–320 K. The calculated results also show that at altitudes of 5, 10 and 15 km, the effective rate constant of the HCOOH-catalyzed reaction increases obviously with an increase in altitude. At 15 km altitude, its value is up to 9.63 × 10⁻¹¹ cm³ per molecule per s, which is close to the corresponding value of the reaction without a catalyst, showing that the contribution of HCOOH to the HO₂ + Cl → HCl + O₂ reaction cannot be neglected at high altitudes. The new findings in this investigation are not only of great necessity and importance for elucidating the gas-phase reaction of HO₂ with Cl in the presence of acidic, neutral and basic catalysts, but are also of great interest for understanding the importance of other types of hydrogen abstraction in the atmosphere.

The effects of acidic (FA), neutral (WM) and basic (AM) catalysts on the energetic and kinetic aspects of the HO₂ + Cl reaction have been studied. At 298 K, the catalytic order of FA, WM and AM is WM > FA > AM.

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.

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.