Amines in the Earth’s Atmosphere: A Density Functional Theory Study of the Thermochemistry of Pre-Nucleation Clusters

Improved Configurational Sampling Protocol for Large Atmospheric Molecular Clusters

The nucleation process leading to the formation of new atmospheric particles plays a crucial role in aerosol research. Quantum chemical (QC) calculations can be used to model the early stages of aerosol formation, where atmospheric vapor molecules interact and form stable molecular clusters. However, QC calculations heavily depend on the chosen computational method, and when dealing with large systems, striking a balance between accuracy and computational cost becomes essential. We benchmarked the binding energies and structures and found the B97-3c method to be a good compromise between the accuracy and computational cost for studying large cluster systems. Further, we carefully assessed configurational sampling procedures for targeting large atmospheric molecular clusters containing up to 30 molecules (approximately 2 nm in diameter) and proposed a funneling approach with highly improved accuracy. We find that several parallel ABCluster explorations lead to better guesses for the cluster global energy minimum structures than one long exploration. This methodology allows us to bridge computational studies of molecular clusters, which typically reach only around 1 nm, with experimental studies that often measure particles larger than 2 nm. By employing this workflow, we searched for low-energy configurations of large sulfuric acid-ammonia and sulfuric acid-dimethylamine clusters. We find that the binding free energies of clusters containing dimethylamine are unequivocally more stable than those of the ammonia-containing clusters. Our improved configurational sampling protocol can in the future be applied to study the growth and dynamics of large clusters of arbitrary compositions.

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.

Clusterome: A Comprehensive Data Set of Atmospheric Molecular Clusters for Machine Learning Applications

Formation and growth of atmospheric molecular clusters into aerosol particles impact the global climate and contribute to the high uncertainty in modern climate models. Cluster formation is usually studied using quantum chemical methods, which quickly becomes computationally expensive when system sizes grow. In this work, we present a large database of ∼250k atmospheric relevant cluster structures, which can be applied for developing machine learning (ML) models. The database is used to train the ML model kernel ridge regression (KRR) with the FCHL19 representation. We test the ability of the model to extrapolate from smaller clusters to larger clusters, between different molecules, between equilibrium structures and out-of-equilibrium structures, and the transferability onto systems with new interactions. We show that KRR models can extrapolate to larger sizes and transfer acid and base interactions with mean absolute errors below 1 kcal/mol. We suggest introducing an iterative ML step in configurational sampling processes, which can reduce the computational expense. Such an approach would allow us to study significantly more cluster systems at higher accuracy than previously possible and thereby allow us to cover a much larger part of relevant atmospheric compounds.

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.

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

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

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.

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

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

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

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

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.

N-nitration of secondary aliphatic amines in the particle phase

Amines are frequently detected in atmospheric particles and are internally mixed with other particle-phase components. However, research on the further reactions of amine with reactive species after entering the particle phase is still limited. This study investigated the nitration reaction process of particulate dimethylamine (DMA), formed via a substitution reaction between DMA and (NH₄)₂SO₄, with NO_x. In situ attenuated total reflectance-infrared Fourier transform spectroscopy (in situ ATR-FTIR) and proton transfer reaction mass spectroscopy (PTR-MS), as well as DFT methods at the B3LYP level using the 6-311++G (d, p) basis set, were mainly used to confirm the formation of nitramine and nitrosamine in the nitration/nitrosation process of DMA. A hydrogen-bonding intermediate ([(CH₃)₂N⋯HONO]) is initially formed when particulate DMA reacts with NO₂ followed by aminyl radical formation, and then nitr- and nitros-amine form through addition reactions with NO₂ and NO, respectively. The dimer of NO₂ (i.e., N₂O₄) and the product of NO and NO₂ (i.e., N₂O₃) can also react with DMA to attack the lone pair electrons on the central N atom of DMA to finally form nitr- and nitros-amine. This study helps reveal the nitration reaction mechanism of organic amines in the particle phase. It also aids in understanding the process of nitrogen cycling in the atmosphere.

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.

Rapid sulfuric acid-dimethylamine nucleation enhanced by nitric acid in polluted regions

Recent research [Wang et al., Nature 581, 184-189 (2020)] indicates nitric acid (NA) can participate in sulfuric acid (SA)-ammonia (NH₃) nucleation in the clean and cold upper free troposphere, whereas NA exhibits no obvious effects at the boundary layer with relatively high temperatures. Herein, considering that an SA-dimethylamine (DMA) nucleation mechanism was detected in megacities [Yao et al., Science 361, 278-281 (2018)], the roles of NA in SA-DMA nucleation are investigated. Different from SA-NH₃ nucleation, we found that NA can enhance SA-DMA-based particle formation rates in the polluted atmospheric boundary layer, such as Beijing in winter, with the enhancement up to 80-fold. Moreover, we found that NA can promote the number concentrations of nucleation clusters (up to 27-fold) and contribute 76% of cluster formation pathways at 280 K. The enhancements on particle formation by NA are critical for particulate pollution in the polluted boundary layer with relatively high NA and DMA concentrations.

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.

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.

A theoretical study of hydrogen-bonded molecular clusters of sulfuric acid and organic acids with amides

Amides, a series of significant atmospheric nitrogen-containing volatile organic compounds (VOCs), can participate in new particle formation (NPF) throught interacting with sulfuric acid (SA) and organic acids. In this study, we investigated the molecular interactions of formamide (FA), acetamide (AA), N-methylformamide (MF), propanamide (PA), N-methylacetamide (MA), and N,N-dimethylformamide (DMF) with SA, acetic acid (HAC), propanoic acid (PAC), oxalic acid (OA), and malonic acid (MOA). Global minimum of clusters were obtained through the association of the artificial bee colony (ABC) algorithm and density functional theory (DFT) calculations. The conformational analysis, thermochemical analysis, frequency analysis, and topological analysis were conducted to determine the interactions of hydrogen-bonded molecular clusters. The heterodimers formed a hepta or octa membered ring through four different types of hydrogen bonds, and the strength of the bonds are ranked in the following order: SOH•••O > COH•••O > NH•••O > CH•••O. We also evaluated the stability of the clusters and found that the stabilization effect of amides with SA is weaker than that of amines with SA but stronger than that of ammonia (NH₃) with SA in the dimer formation of nucleation process. Additionally, the nucleation capacity of SA with amides is greater than that of organic acids with amides.

Seed-Adsorbate Interactions as the Key of Heterogeneous Butanol and Diethylene Glycol Nucleation on Ammonium Bisulfate and Tetramethylammonium Bromide

Condensation particle counter (CPC) instruments are commonly used to detect atmospheric nanoparticles. They operate on the basis of condensing an organic working fluid on the nanoparticle seeds to grow the particles to a detectable size, and at the size of few nanometers, their efficiency depends on how well the working fluid interacts with the seeds under the measurement conditions. This study models the first steps of heterogeneous nucleation of two working fluids commonly used in CPCs (diethylene glycol (DEG) and n-butanol) onto two positively charged seeds, ammonium bisulfate and tetramethylammonium bromide. The nucleation process is modeled on a molecular level using a combination of systematic configurational sampling and density functional theory (DFT). We take into account the conformational flexibility of DEG and n-butanol and determine the key factors that can improve the efficiency of nanoparticle measurements by CPCs. The results show that hydrogen bonding between the seed and the working fluid molecules is central to the adsorption of the first DEG/n-butanol molecules onto the seeds. However, intermolecular hydrogen bonding between the adsorbed molecules can also enhance the nucleation process for the weakly adsorbing vapor molecules. Accordingly, the heterogeneous nucleation probability is higher for working fluid-nanoparticle combinations with a higher potential for hydrogen bonding; in this case, DEG and ammonium bisulfate. Moreover, conformational analysis and methodology evaluations indicate that the consideration of adsorbate conformers and step-wise addition of the vapor molecules to the seeds is not essential for qualitative modeling of heterogeneous nucleation systems, at least for systems where the adsorbate and seed chemical properties are clearly different. This is the first molecular-level modeling study reporting detailed chemical reasons for experimentally observed seed and working fluid preferences in CPCs and reproducing the experimental observations. Our presented approach can be likely used for predicting preferences in similar nucleating systems.

Determinant Factor for Thermodynamic Stability of Sulfuric Acid-Amine Complexes

Atmospheric amines are thought to play significant roles in the nucleation of sulfuric acid-mediated aerosol particles. Their enhancing effects on the stabilization of the related complexes have formerly been correlated with the amine base strength, but there are a few exceptions reported. In this work, the influence of seven alkylamines on the thermodynamic stability of sulfuric acid-amine complexes has been theoretically investigated, e.g., ethylamine, propylamine, isopropylamine, tert-butylamine, dimethylamine, ethylmethylamine, and trimethylamine. For all primary and secondary amine-mediated complexes, a dual hydrogen bond configuration is generally suggested in the most stable isomer. The stabilization of this special structure predicted by the electrostatic potential distribution on the molecular surface of amines exactly agrees with the base strength sequence, providing crucial evidence for the previous deduction of correlation between the base strength and the enhancing effect. Meanwhile, the considerable van der Waals interactions are found between the free hydroxyl of sulfuric acid and the β-methyl group of amine, resulting in the extra stability for sulfuric acid-dimethylamine and sulfuric acid-ethylmethylamine complexes. Therefore, the electrostatic potential distribution of amines is the essential determinant factor for the thermodynamic stability of the relevant complexes. Our conclusions provide new insight into a way to evaluate the enhancing abilities of amines in aerosol particle nucleation.

Interactions of sulfuric acid with common atmospheric bases and organic acids: Thermodynamics and implications to new particle formation

Interactions of the three common atmospheric bases, dimethylamine ((CH₃)₂NH), methylamine (CH₃NH₂), ammonia (NH₃), all considered to be efficient stabilizers of binary clusters in the Earth's atmosphere, with H₂SO₄, the key atmospheric precursor, and 14 common atmospheric organic acids (COAs) (formic, acetic, oxalic, malonic, succinic, glutaric acid, adipic, benzoic, phenylacetic, pyruvic, maleic acid, malic, tartaric and pinonic acids) have been studied using the density functional theory (DFT) and composite high-accuracy G3MP2 method. The thermodynamic stability of mixed (COA)(H₂SO₄), (COA)(B1), (COA)(B2) and (COA)(B3) dimers and (COA)(H₂SO₄)(B1), (COA)(H₂SO₄)(B2) and (COA)(H₂SO₄)(B3) trimers, where B1, B2 and B3 refer to (CH₃)₂NH, CH₃NH₂ and NH₃, respectively, have been investigated and their impacts on the thermodynamic stability of clusters containing H₂SO₄ have been studied. Our investigation shows that interactions of H₂SO₄ with COA, (CH₃)₂NH, CH₃NH₂ and NH₃ lead to the formation of more stable mixed dimers and trimers than (H₂SO₄)₂ and (H₂SO₄)₂(base), respectively, and emphasize the importance of common organic species for early stages of atmospheric nucleation. We also show that although amines are generally confirmed to be more active than NH₃ as stabilizers of binary clusters, in some cases mixed trimers containing NH₃ are more stable thermodynamically than those containing CH₃NH₂. This study indicates an important role of COA, which coexist and interact with that H₂SO₄ and common atmospheric bases in the Earth atmosphere, in formation of stable pre-nucleation clusters and suggests that the impacts of COA on new particle formation (NPF) should be studied in further details.

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.

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

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

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

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

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.

A molecular understanding of the interaction of typical aromatic acids with common aerosol nucleation precursors and their atmospheric implications

Aromatic acids, which are generated from numerous anthropogenic emissions and secondary transformations, have been considered to play a crucial role in new particle formation. In this study, we performed theoretical calculations at the PW91PW91/6-311++G(3df,3pd) level to investigate the interaction between typical aromatic acids namely benzoic acid (BA), phenylacetic acid (PAA), phthalic acid (PA), isophthalic acid (mPA), and terephthalic acid (PTA) and common atmospheric nucleation precursors namely sulfuric acid (SA), water (H2O), ammonia (NH3), methylamine (MA), dimethylamine (DMA), and trimethylamine (TMA). The geometric analysis, Gibbs free energy analysis, OH/NH-stretching vibrational frequency calculation, and atoms in molecules (AIM) analysis were conducted to determine the interactions in the complexes. The heterodimers formed a six to eight membered ring through four types of hydrogen bond, and the bond strength could be ranked in descending order: SO-H⋯O > O-H⋯O/N > N-H⋯O. The BA/PAA/mPA/PTA-SA complexes had the lowest Gibbs free energy values. PA was more likely to interact with NH3 or amines rather than SA due to an intra-molecular hydrogen bond. Additionally, the aromatic acids have similar ability to interact with SA and NH3 as monocarboxylic/dicarboxylic acid. The formation potential of the heterodimers from aromatic acids with common nucleation precursors in ambient atmosphere was investigated.

A density functional theory study of the molecular interactions between a series of amides and sulfuric acid

Amides, a class of nitrogen-containing organic pollutants in the atmosphere, may affect the formation of atmospheric aerosols by the interactions with sulfuric acid. Here, the molecular interactions of sulfuric acid with formamide, methylformamide, dimethylformamide, acetamide, methylacetamide and dimethylacetamide was investigated by density functional theory. Geometry optimization and Gibbs free energy calculation were carried out at M06-2X/6-311++G(3df,3pd) level. The results indicate that the addition of amides to H₂SO₄ might have a promoting effect on atmospheric new particle formation at 298.15 K and 1 atm. In the initial stage of new particle formation, the binding capacity of amides and sulfuric acid is stronger than ammonia, but weaker than methylamine. It is worth noting that the trans-methylacetamide could have similar capabilities of stabilizing sulfuric acid as dimethylamine. In the presence of water, amides are found to only have a weak enhancement capability on new particle formation. In addition, we can infer from evaporation rate that the small molecule clusters of formamide and sulfuric acid may be more energetically favorable than macromolecule clusters.

The role of nitric acid in atmospheric new particle formation

Nitric acid, an air pollutant with strong acidity and oxidizability, can be found in considerable quantities in the gas and aerosol phase. Understanding the role of nitric acid in atmospheric new particle formation is essential to study the complicated nucleation mechanism. Using density functional theory combined with the Atmospheric Clusters Dynamic Code (ACDC), the role of nitric acid in the formation of new particles has been investigated under different atmospheric conditions (different precursor concentrations and temperatures). The results show that nitric acid can form clusters with sulfuric acid and ammonia by hydrogen bond or even proton-transfer interactions. The concentrations of clusters involving nitric acid can be comparable with those of sulfuric acid-ammonia-based clusters, considering the thermodynamic stability combined with the realistic atmospheric concentrations of precursors. Within the atmospheric concentration range, nitric acid can enhance the formation rates of sulfuric acid-ammonia clusters, especially at low temperature, low sulfuric acid concentration and high ammonia concentration. In addition, the new particle formation mechanism indicates that nitric acid can contribute to the cluster formation and the role of nitric acid in the cluster formation pathway is as a "bridge" connecting the smaller and larger clusters.

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

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

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

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

Hydration of 3-hydroxy-4,4-dimethylglutaric acid with dimethylamine complex and its atmospheric implications

Changes in temperature affects the distribution of isomers, which facilitates the understanding of new particle formation in the atmosphere.

Uptake of water by an acid–base nanoparticle: theoretical and experimental studies of the methanesulfonic acid–methylamine system

Uptake of water by nanoparticles composed by methanesulfonic acid and methylamine using a combination of theoretical calculations and laboratory experiments.

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

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

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

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

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

Atmospheric Fate of Monoethanolamine: Enhancing New Particle Formation of Sulfuric Acid as an Important Removal Process

Monoethanolamine (MEA), a potential atmospheric pollutant from the capture unit of a leading CO₂ capture technology, could be removed by participating H₂SO₄-based new particle formation (NPF) as simple amines. Here we evaluated the enhancing potential of MEA on H₂SO₄-based NPF by examining the formation of molecular clusters of MEA and H₂SO₄ using combined quantum chemistry calculations and kinetics modeling. The results indicate that MEA at the parts per trillion (ppt) level can enhance H₂SO₄-based NPF. The enhancing potential of MEA is less than that of dimethylamine (DMA), one of the strongest enhancing agents, and much greater than methylamine (MA), in contrast to the order suggested solely by their basicity (MEA < MA < DMA). The unexpectedly high enhancing potential is attributed to the role of -OH of MEA in increasing cluster binding free energies by acting as both a hydrogen bond donor and acceptor. After the initial formation of one H₂SO₄ and one MEA cluster, the cluster growth mainly proceeds by first adding one H₂SO₄, and then one MEA, which differs from growth pathways in H₂SO₄-DMA and H₂SO₄-MA systems. Importantly, the effective removal rate of MEA due to participation in NPF is comparable to that of oxidation by hydroxyl radicals at 278.15 K, indicating NPF as an important sink for MEA.

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

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

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

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

Interaction of oxalic acid with dimethylamine and its atmospheric implications

Oxalic acid and dimethylamine are the most common organic acid and base in the atmosphere, and are recognized as significant precursor species in atmospheric new particle formation.

Hydrogen Bonding Interaction between Atmospheric Gaseous Amides and Methanol

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