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

Implications for new particle formation in air of the use of monoethanolamine in carbon capture and storage

Alkanolamines are currently being deployed in carbon capture and storage (CCS) technology worldwide, and atmospheric emissions have been found to coincide with locations exhibiting elevated concentrations of methanesulfonic acid (MSA). It is thus critical to understand the fate and potential atmospheric reactions of these chemicals. This study reports the characterization of sub-10 nm nanoparticles produced through the acid-base reaction between gas phase monoethanolamine (MEA) and MSA, a product of organosulfur compound oxidation in air, using a flow reactor under dry and humid (up to ∼60% RH) conditions. Number size distribution measurements show that MEA is even more efficient than methylamine in forming nanoparticles on reaction with MSA. This is attributed to the fact that the MEA structure contains both an -NH2 and an -OH group that facilitate hydrogen bonding within the clusters, in addition to the electrostatic interactions. Due to this already strong H-bond network, water has a relatively small influence on new particle formation (NPF) and growth in this system, in contrast to MSA reactions with alkylamines. Acid/base molar ratios of unity for 4-12 nm particles were measured using thermal desorption chemical ionization mass spectrometry. The data indicate that reaction of MEA with MSA may dominate NPF under some atmospheric conditions. Thus, the unique characteristics of alkanolamines in NPF must be taken into account for accurate predictions of impacts of CCS on visibility, health and climate.

HIO3-HIO2-Driven Three-Component Nucleation: Screening Model and Cluster Formation Mechanism

Iodine oxoacids (HIO3 and HIO2)-driven nucleation has been suggested to efficiently contribute to new particle formation (NPF) in marine atmospheres. Abundant atmospheric nucleation precursors may further enhance HIO3-HIO2-driven nucleation through various multicomponent nucleation mechanisms. However, the specific enhancing potential (EP) of different precursors remains largely unknown. Herein, the EP-based screening model of precursors and enhancing mechanism of the precursor with the highest EP on HIO3-HIO2 nucleation were investigated. The formation free energies (ΔG), as critical parameters for evaluating EP, were calculated for the dimers of 63 selected precursors with HIO2. Based on the ΔG values, (1) a quantitative structure-activity relationship model was developed for evaluating ΔG of other precursors and (2) atmospheric concentrations of 63 (precursor)1(HIO2)1 dimer clusters were assessed to identify the precursors with the highest EP for HIO3-HIO2-driven nucleation by combining with earlier results for the nucleation with HIO3 as the partner. Methanesulfonic acid (MSA) was found to be one of the precursors with the highest EP. Finally, we found that MSA can effectively enhance HIO3-HIO2 nucleation at atmospheric conditions by studying larger MSA-HIO3-HIO2 clusters. These results augment our current understanding of HIO3-HIO2 and MSA-driven nucleation and may suggest a larger impact of HIO2 in atmospheric aerosol nucleation.

Enhancement of Atmospheric Nucleation Precursors on Iodic Acid-Induced Nucleation: Predictive Model and Mechanism

Iodic acid (IA) has recently been recognized as a key driver for new particle formation (NPF) in marine atmospheres. However, the knowledge of which atmospheric vapors can enhance IA-induced NPF remains limited. The unique halogen bond (XB)-forming capacity of IA makes it difficult to evaluate the enhancing potential (EP) of target compounds on IA-induced NPF based on widely studied sulfuric acid systems. Herein, we employed a three-step procedure to evaluate the EP of potential atmospheric nucleation precursors on IA-induced NPF. First, we evaluated the EP of 63 precursors by simulating the formation free energies (ΔG) of the IA-containing dimer clusters. Among all dimer clusters, 44 contained XBs, demonstrating that XBs are frequently formed. Based on the calculated ΔG values, a quantitative structure-activity relationship model was developed for evaluating the EP of other precursors. Second, amines and O/S-atom-containing acids were found to have high EP, with diethylamine (DEA) yielding the highest potential to enhance IA-induced nucleation by combining both the calculated ΔG and atmospheric concentration of considered 63 precursors. Finally, by studying larger (IA)_1-3(DEA)_1-3 clusters, we found that the IA-DEA system with merely 0.1 ppt (2.5×10⁶ cm^-3) DEA yields comparable nucleation rates to that of the IA-iodous acid system.

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.

Formation mechanism of typical aromatic sulfuric anhydrides and their potential role in atmospheric nucleation process

Sulfuric anhydrides, generated from the cycloaddition reaction of SO₃ with carboxylic acids, have been revealed to be potential participants in the nucleation process of new particle formation (NPF). Hence the reaction mechanisms of typical aromatic acids (benzoic acid (BA), phenylacetic acid (PAA), phthalic acid (PA), isophthalic acid (mPA), and terephthalic acid (PTA)) with SO₃ to generate the corresponding aromatic sulfuric anhydrides were investigated by density functional theory calculations at the level of M06-2X/6-311++G(3df,3pd). As a result, these reactions were found to be feasible in the gas phase with barriers of 0.34, 0.30, 0.18, 0.08 and 0.12 kcal/mol to generate corresponding aromatic sulfuric anhydrides, respectively. The thermodynamic stabilities of clusters containing aromatic sulfuric anhydrides and atmospheric nucleation precursors (sulfuric acid, ammonia and dimethylamine) were further analyzed to identify the potential role of aromatic sulfuric anhydrides in NPF. As the thermodynamic stability of a cluster depends on both the number and strength of hydrogen bonds, the greater stability of the interactions between atmospheric nucleation precursors and aromatic sulfuric anhydrides than with aromatic acids make aromatic sulfuric anhydrides potential participators in the nucleation process of NPF. Moreover, compared with BA, the addition of a -CH₂- functional group in PAA has little influence on the reaction barrier with SO₃ but an inhibitive effect on the thermodynamic stability of clusters. The position of the two -COOH functional groups in PA, mPA and PTA does not have a consistent impact on the reaction barrier with SO₃ or the thermodynamic stability.

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.

Critical Role of Iodous Acid in Neutral Iodine Oxoacid Nucleation

Nucleation of neutral iodine particles has recently been found to involve both iodic acid (HIO₃) and iodous acid (HIO₂). However, the precise role of HIO₂ in iodine oxoacid nucleation remains unclear. Herein, we probe such a role by investigating the cluster formation mechanisms and kinetics of (HIO₃)_m(HIO₂)_n (m = 0-4, n = 0-4) clusters with quantum chemical calculations and atmospheric cluster dynamics modeling. When compared with HIO₃, we find that HIO₂ binds more strongly with HIO₃ and also more strongly with HIO₂. After accounting for ambient vapor concentrations, the fastest nucleation rate is predicted for mixed HIO₃-HIO₂ clusters rather than for pure HIO₃ or HIO₂ ones. Our calculations reveal that the strong binding results from HIO₂ exhibiting a base behavior (accepting a proton from HIO₃) and forming stronger halogen bonds. Moreover, the binding energies of (HIO₃)_m(HIO₂)_n clusters show a far more tolerant choice of growth paths when compared with the strict stoichiometry required for sulfuric acid-base nucleation. Our predicted cluster formation rates and dimer concentrations are acceptably consistent with those measured by the Cosmic Leaving Outdoor Droplets (CLOUD) experiment. This study suggests that HIO₂ could facilitate the nucleation of other acids beyond HIO₃ in regions where base vapors such as ammonia or amines are scarce.

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.

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.

A theoretical study on the formation mechanism of carboxylic sulfuric anhydride and its potential role in new particle formation

New particle formation (NPF) is the major source of atmospheric aerosol particles. However, the chemical species involved and the exact mechanism are still unclear. Cycloaddition reaction of SO₃ to carboxylic acids bas been identified as a possible formation mechanism of carboxylic sulfuric anhydrides which may be involved in NPF. Herein, energy profiles for forming diaterpenylic acetate sulfuric anhydride (DTASA) through cycloaddition of SO₃ to diaterpenylic acid acetate (DTAA) and the potential role of DTASA in NPF were studied through computational methods combined with atmospheric cluster dynamics code (ACDC). Gas phase reaction barriers for the two carboxyl groups of DTAA are 0.4 and 0.6 kcal mol⁻¹, respectively, illustrating a feasible formation mechanism for DTASA. According to thermodynamical analysis and dynamical simulations, atmospheric clusters containing DTASA and atmospheric nucleation precursors sulfuric acid (SA), ammonia (NH₃) and dimethylamine (DMA) possess both thermodynamically and dynamically higher stabilities than those of DTAA-contained clusters. Furthermore, DTASA–NH₃ and DTASA–DMA are more stable than SA–NH₃ and SA–DMA, enabling DTASA, even carboxylic sulfuric anhydrides, to become potential participants in the atmospheric NPF process which may hence promote a better understanding of NPF.

Organic acids could improve their nucleation ability through the cycloaddition reaction of SO₃ to generate corresponding carboxylic sulfuric anhydrides which may play a potential role in the atmospheric new particle formation.

Organic acid-ammonia ion-induced nucleation pathways unveiled by quantum chemical calculation and kinetics modeling: A case study of 3-methyl-1,2,3-butanetricarboxylic acid

Nucleation of organic acids (OAs) and H₂SO₄ is an important source for new particle formation in the atmosphere. However, it is still unclear whether organic acids can produce nanoparticles independent of H₂SO₄. In this study, 3-methyl-1,2,3-butanetricarboxylic acid (MBTCA) was adopted as a model of OAs. Pathways of clustering from MBTCA, ammonia and ions (NH₄⁺ and NO₃^-) to form a 1.9 nm nucleus were investigated by quantum chemical calculation and kinetic modeling. Results show recombination of charged clusters/ions plays an essential role in the nucleation processes. Cluster formation rates increase by a factor of 10³ when NH₃ increases from 2.6 × 10⁸ molecules·cm^-3 (under clean conditions) to 2.6 × 10¹¹ molecules·cm^-3 (under polluted conditions), as NH₃ can stabilize MBTCA clusters and change ion compositions from H₃O⁺ to NH₄⁺. Although the proposed new mechanism cannot compete with H₂SO₄-NH₃-H₂O or H₂SO₄-OA nucleation currently, it may be important in the future with the decline of SO₂ concentration.

Vibrational Spectra of Atmospherically Relevant Hydrogen-Bonded MSA···(H2SO4)n (n = 1-3) Clusters

Methanesulfonic acid (CH₃SO₃H), also known as MSA, has been found to be capable of forming a strong hydrogen-bonded interaction with sulfuric acid (H₂SO₄) under ambient conditions. The energetic stability of the MSA···H₂SO₄ clusters increases with decreasing temperature at higher altitudes in the troposphere, which is relevant in the context of atmospheric aerosol formation. We have performed, in the present work, a detailed and systematic quantum-chemical calculation with high-level density functional theory to characterize the hydrogen bond formation in the binary MSA···H₂SO₄, ternary MSA···(H₂SO₄)₂, and quaternary MSA···(H₂SO₄)₃ clusters. The five different conformations of MSA···(H₂SO₄)₂ and six conformations of MSA···(H₂SO₄)₃, considered in the present work for the spectroscopic analysis, have been taken from our previous work [J. Phys. Chem. A. 2020, 124, 11072-11085]. The hydrogen bonds were analyzed on the basis of infrared vibrational frequencies of different O-H stretching modes and quantum theory of atoms in molecules (QTAIM). A strong positive correlation has been observed between the red shift of the OH groups in MSA and H₂SO₄ and the corresponding O-H elongation as a result of hydrogen bond formation. Topological analysis employing QTAIM shows that most of the charge density and the Laplacian values at bond critical points (BCPs) of the hydrogen bonds of the MSA···(H₂SO₄)_n (n = 1-3) complexes fall within the standard hydrogen-bond criteria. However, those outside these criteria fall in the category of a very strong hydrogen bond with a hydrogen bond length as low as 1.41 Å and an O-H bond elongation as high as 0.096 Å. In general, the charge densities of the BCPs located on hydrogen bonds increase as the hydrogen-bond lengths decrease. Proportionately, a larger number of hydrogen bonds in ternary MSA···(H₂SO₄)₂ demonstrate a partial covalent character when compared with the quaternary clusters.

TiO2 Coating Strategy for Robust Catalysis of the Metal-Organic Framework toward Energy-Efficient CO2 Capture

High energy duty restricts the application of amine-based absorption in CO2 capture and limits the achievement of carbon neutrality. Although regenerating the amine solvent with solid acid catalysts can increase energy efficiency, inactivation of the catalyst must be addressed. Here, we report a robust metal-organic framework (MOF)-derived hybrid solid acid catalyst (SO42-/ZIF-67-C@TiO2) with improved acidity for promoting amine regeneration. The TiO2 coating effectively prevented the active components stripping from the surface of the catalyst, thus prolonging its lifespan. The well-protected Co-Nx sites and protonated groups introduced onto the TiO2 surface increased the amount and rate of CO2 desorption by more than 64.5 and 153%, respectively. Consequently, the energy consumption decreased by approximately 36%. The catalyzed N-C bond rupture and proton transfer mechanisms are proposed. This work provides an effective protection strategy for robust acid catalysts, thus advancing the CO2 capture with less energy duty.

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.

Autoxidation mechanism for atmospheric oxidation of tertiary amines: Implications for secondary organic aerosol formation

Tertiary amines are one kind of identified amines in the atmosphere. Here, the atmospheric oxidation mechanism and kinetics of tertiary amines were investigated by using computational methods. As proxies of these amines, trimethylamine (TMA) and triethylamine (TEA) have been selected. Results indicate that N-containing peroxy radicals (NRO₂⋅), which are key intermediates in ⋅OH initiated oxidation of TMA and TEA, can follow a so-called autoxidation mechanism (a chain reaction of H-shift followed by O₂ addition) even on the condition of high NO/HO₂⋅ concentration. Such unique mechanism can be ascribed to the ability of N-atom in facilitating the unimolecular H-shift of NRO₂⋅ and the absence of H-atoms on N-atom. However, different from TMA reaction system, the pathway dissociating into fragmental products can compete with the autoxidation pathway for TEA system. More importantly, TEA reaction system cannot lead to the formation of products with high O/C ratio due to the autoxidation pathway terminated by the release of fragmental molecules. Such difference can be corroborated by previously observing lower secondary organic aerosol yield of TEA oxidation than that of TMA oxidation. The unveiled mechanism enhances current understanding on atmospheric fate of amines and autoxidation mechanism.

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.

Atmospherically Relevant Hydrogen-Bonded Interactions between Methanesulfonic Acid and H2SO4 Clusters: A Computational Study

A detailed and systematic quantum-chemical calculation has been performed with high-level density functional theory (DFT) to analyze the electrostatic interaction of methanesulfonic acid (CH₃SO₃H), also known as MSA, with pre-formed clusters of sulfuric acid (H₂SO₄) molecules in ambient conditions. Both MSA and H₂SO₄ are considered as atmospheric molecules that might play active roles in aerosol formation. The interactions between MSA and H₂SO₄ clusters lead to the formation of MSA···(H₂SO₄)_n (n = 2, 3) complexes stabilized by the formation of different types of intermolecular hydrogen bond networks. Analyses of cluster binding energies and free energy changes associated with their formation indicate that MSA could bring additional stability into the atmospheric molecular clusters responsible for aerosol formation. Variations of Gibbs free energy with temperature and pressure have been analyzed. The lower temperatures and pressures at the higher altitudes of the troposphere are found to play in favor of higher stability of the MSA···(H₂SO₄)_n clusters. Effects of hydrogen bond formation on dipole moment, mean polarizability, and anisotropy of polarizability of the clusters have been analyzed. Rayleigh scattering intensities are found to increase many-fold when light interacts with the MSA···(H₂SO₄)_n clusters.

Structural Effects of Amines in Enhancing Methanesulfonic Acid-Driven New Particle Formation

Atmospheric amines can enhance methanesulfonic acid (MSA)-driven new particle formation (NPF), but the mechanism is fundamentally different compared to that of the extensively studied sulfuric acid (SA)-driven process. Generally, the enhancing potentials of amines in SA-driven NPF follow the basicity, while this is not the case for MSA-driven NPF, where structural effects dominate, making MSA-driven NPF more prominent for methylamine (MA) compared to dimethylamine (DMA). Therefore, probing structural factors determining the enhancing potentials of amines on MSA-driven NPF is key to fully understanding the contribution of MSA to NPF. Here, we performed a comparative study on DMA and MA enhancing MSA-driven NPF by examining cluster formation using computational methods. The results indicate that DMA-MSA clusters are more stable than the corresponding MA-MSA clusters for cluster sizes up to (DMA)₂(MSA)₂, indicating that the basicity of amines dominates the initial cluster formation. The methyl groups of DMA were found to present significant steric hindrance beyond the (DMA)₂(MSA)₂ cluster and this adds to the lower hydrogen bonding capacity of DMA, making the cluster growth less favorable compared to MA. This study implies that several amines could synergistically enhance MSA-driven NPF by maximizing the advantage of different amines in different amine-MSA cluster growth stages.

Formation Mechanisms of Iodine-Ammonia Clusters in Polluted Coastal Areas Unveiled by Thermodynamics and Kinetic Simulations

It has been revealed that iodine species play important roles in atmospheric new particle formations (NPFs) in pristine coastal areas. However, it is unclear whether other atmospheric species, such as NH₃, for which the levels in coastal areas of China are >2.5 × 10¹⁰ molecules·cm^-3 are involved in the NPFs of iodine species, although NH₃ has been proved to promote particle formation of H₂SO₄. Via high-level quantum chemical calculations and atmospheric cluster dynamic code simulations, this study unveiled new mechanisms of nucleation, in which NH₃ mediates the formation of iodine particles by assisting hydrolysis of I₂O₅ or reacting with HIO₃. The simulated formation rates of iodine-ammonia clusters via the new mechanisms are much higher than those simulated via sequential addition of HIO₃ with subsequent release of H₂O, under the condition that NH₃ concentrations are higher than 10¹⁰ molecules·cm^-3. The new mechanisms can well explain the observed cluster formation rates at a coastal site in Zhejiang of China. The findings not only expand the current understandings of the role of NH₃ in NPFs but also highlight the importance of monitoring and evaluating NPFs via the iodine-ammonia cluster pathway in the coastal areas of China and other regions worldwide.

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.