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.

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.

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×106 cm-3) DEA yields comparable nucleation rates to that of the IA-iodous acid system.

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.

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(T0)/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.

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)1(MSA)1 cluster, and subsequently by cluster collisions with either a (Gud)1(MSA)1 or (Gud)2(MSA)2 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.

Atmospheric Autoxidation of Organophosphate Esters.

Organophosphate esters (OPEs), widely used as flame retardants and plasticizers, have frequently been identified in the atmosphere. However, their atmospheric fate and toxicity associated with atmospheric transformations are unclear. Here, we performed quantum chemical calculations and computational toxicology to investigate the reaction mechanism of peroxy radicals of OPEs (OPEs-RO2•), key intermediates in determining the atmospheric chemistry of OPEs, and the toxicity of the reaction products. TMP-RO2• (R1) and TCPP-RO2• (R2) derived from trimethyl phosphate and tris(2-chloroisopropyl) phosphate, respectively, are selected as model systems. The results indicate that R1 and R2 can follow an H-shift-driven autoxidation mechanism under low NO concentration ([NO]) conditions, clarifying that RO2• from esters can follow an autoxidation mechanism. The unexpected autoxidation mechanism can be attributed to the distinct role of the ─(O)3P(═O) phosphate-ester group in facilitating the H-shift of OPEs-RO2• from commonly encountered ─OC(═O)─ and ─ONO2 ester groups in the atmosphere. Under high [NO] conditions, NO can mediate the autoxidation mechanism to form organonitrates and alkoxy radical-related products. The products from the autoxidation mechanism have low volatility and aquatic toxicity compared to their corresponding parent compounds. The proposed autoxidation mechanism advances our current understanding of the atmospheric RO2• chemistry and the environmental risk of OPEs.

Ozonolysis of α-Pinene and Δ3-Carene Mixtures: Formation of Dimers with Two Precursors.

The formation of secondary organic aerosol (SOA) from the structurally similar monoterpenes, α-pinene and Δ3-carene, differs substantially. The aerosol phase is already complex for a single precursor, and when mixtures are oxidized, products, e.g., dimers, may form between different volatile organic compounds (VOCs). This work investigates whether differences in SOA formation and properties from the oxidation of individual monoterpenes persist when a mixture of the monoterpenes is oxidized. Ozonolysis of α-pinene, Δ3-carene, and a 1:1 mixture of them was performed in the Aarhus University Research on Aerosol (AURA) atmospheric simulation chamber. Here, ∼100 ppb of monoterpene was oxidized by 200 ppb O3 under dark conditions at 20 °C. The particle number concentration and particle mass concentration for ozonolysis of α-pinene exceed those from ozonolysis of Δ3-carene alone, while their mixture results in concentrations similar to α-pinene ozonolysis. Detailed offline analysis reveals evidence of VOC-cross-product dimers in SOA from ozonolysis of the monoterpene mixture: a VOC-cross-product dimer likely composed of the monomeric units cis-caric acid and 10-hydroxy-pinonic acid and a VOC-cross-product dimer ester likely from the monomeric units caronaldehyde and terpenylic acid were tentatively identified by liquid chromatography-mass spectrometry. To improve the understanding of chemical mechanisms determining SOA, it is relevant to identify VOC-cross-products.

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(T0)-calculations.

Contribution of Methanesulfonic Acid to the Formation of Molecular Clusters in the Marine Atmosphere.

Because of the lack of long-term measurements, new particle formation (NPF) in the marine atmosphere remains puzzling. Using quantum chemical methods, this study elucidates the cluster formation and further growth of sulfuric acid-methanesulfonic acid-dimethylamine (SA-MSA-DMA) clusters, relevant to NPF in the marine atmosphere. The cluster structures and thermochemical parameters of (SA)n(MSA)m(DMA)l (n + m ≤ 4 and l ≤ 4) systems are calculated using density functional theory at the ωB97X-D/6-31++G(d,p) level of theory, and the single-point energies are calculated using high-level DLPNO-CCSD(T0)/aug-cc-pVTZ calculations. The calculated thermochemistry is used as input to the Atmospheric Cluster Dynamics Code (ACDC) to gain insight into the cluster dynamics. At ambient conditions (298.15 K, 1 atm), we find that the distribution of outgrowing clusters primarily consists of SA and DMA, with a minor contribution from the mixed SA-MSA-DMA clusters. At lower temperature (278.15 K, 1 atm) the distribution broadens, and clusters containing one or more MSA molecules emerge. These findings show that in the cold marine atmosphere MSA likely participates in atmospheric NPF.

Atmospheric Chemistry of Allylic Radicals from Isoprene: A Successive Cyclization-Driven Autoxidation Mechanism.

The atmospheric chemistry of isoprene has broad implications for regional air quality and the global climate. Allylic radicals, taking 13-17% yield in the isoprene oxidation by •Cl, can contribute as much as 3.6-4.9% to all possible formed intermediates in local regions at daytime. Considering the large quantity of isoprene emission, the chemistry of the allylic radicals is therefore highly desirable. Here, we investigated the atmospheric oxidation mechanism of the allylic radicals using quantum chemical calculations and kinetics modeling. The results indicate that the allylic radicals can barrierlessly combine with O2 to form peroxy radicals (RO2•). Under ≤100 ppt NO and ≤50 ppt HO2• conditions, the formed RO2• mainly undergo two times "successive cyclization and O2 addition" to finally form the product fragments 2-alkoxy-acetaldehyde (C2H3O2•) and 3-hydroperoxy-2-oxopropanal (C3H4O4). The presented reaction illustrates a novel successive cyclization-driven autoxidation mechanism. The formed 3-hydroperoxy-2-oxopropanal product is a new isomer of the atmospheric C3H4O4 family and a potential aqueous-phase secondary organic aerosol precursor. Under >100 ppt NO condition, NO can mediate the cyclization-driven autoxidation process to form C5H7NO3, C5H7NO7, and alkoxy radical-related products. The proposed novel autoxidation mechanism advances our current understanding of the atmospheric chemistry of both isoprene and RO2•.

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.

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)2(MSA)2, 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)2(MSA)2 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 of Low-Volatile Products and Unexpected High Formaldehyde Yield from the Atmospheric Oxidation of Methylsiloxanes.

With stricter regulation of atmospheric volatile organic compounds (VOCs) originating from fossil fuel-based vehicles and industries, the use of volatile chemical products (VCPs) and the transformation mechanism of VCPs have become increasingly important to quantify air quality. Volatile methylsiloxanes (VMS) are an important class of VCPs and high-production chemicals. Using quantum chemical calculations and kinetics modeling, we investigated the reaction mechanism of peroxy radicals of VMS, which are key intermediates in determining the atmospheric chemistry of VMS. L2-RSiCH2O2• and D3-RSiCH2O2• derived from hexamethyldisiloxane and hexamethylcyclotrisiloxane, respectively, were selected as representative model systems. The results indicated that L2-RSiCH2O2• and D3-RSiCH2O2• follow a novel Si-C-O rearrangement-driven autoxidation mechanism, leading to the formation of low volatile silanols and high yield of formaldehyde at low NO/HO2• conditions. At high NO/HO2• conditions, L2-RSiCH2O2• and D3-RSiCH2O2• react with NO/HO2• to form organic nitrate, hydroperoxide, and active alkoxy radicals. The alkoxy radicals further follow a Si-C-O rearrangement step to finally form formate esters. The novel Si-C-O rearrangement mechanism of both peroxy and alkoxy radicals are supported by available experimental studies on the oxidation of VMS. Notably, the high yield of formaldehyde is estimated to significantly contribute to formaldehyde pollution in the indoor environment, especially during indoor cleaning.

The reaction of isotope-substituted hydrated iodide I(HO)- with ozone: the reactive influence of the solvent water molecule.

We report an investigation of the reaction of isotope-substituted hydrated iodide I(HO)- with ozone 16O3 to examine the involvement of the water molecules in the oxidation reactions that terminate with the formation of IO3-. Experimentally, we studied the reaction in the gas phase as elementary reactions using a radio-frequency (RF) ion-trap combined with a quadrupole mass spectrometer (QMS). In approximately 1.2% of the reactions of I(HO)- and 16O3, the 18O atom is found to appear in iodine oxide anions, thus giving evidence for a close involvement of the water molecule in a non-negligible number of the reactions towards IO3-. As a part of the experimental investigation, the reaction rate constant for the exchange reaction I(HO)- + HO → I(HO)- + HO at 300 K was found to be (1.3 ± 0.1) × 10-8 cm3 s-1. Quantum chemical calculations are exploited to establish the energetic difference between I(HO)- and I(HO)-.

Hydration of Atmospheric Molecular Clusters III: Procedure for Efficient Free Energy Surface Exploration of Large Hydrated Clusters.

Sampling the shallow free energy surface of hydrated atmospheric molecular clusters is a significant challenge. Using computational methods, we present an efficient approach to obtain minimum free energy structures for large hydrated clusters of atmospheric relevance. We study clusters consisting of two to four sulfuric acid (sa) molecules and hydrate them with up to five water (w) molecules. The structures of the "dry" clusters are obtained using the ABCluster program to yield a large pool of low-lying conformer minima with respect to free energy. The conformers (up to ten) lowest in free energy are then hydrated using our recently developed systematic hydrate sampling technique. Using this approach, we identify a total of 1145 unique (sa)2-4(w)1-5 cluster structures. The cluster geometries and thermochemical parameters are calculated at the ωB97X-D/6-31++G(d,p) level of theory, at 298.15 K and 1 atm. The single-point energy of the most stable clusters is calculated using a high-level DLPNO-CCSD(T0)/aug-cc-pVTZ method. Using the thermochemical data, we calculate the equilibrium hydrate distribution of the clusters under atmospheric conditions and find that the larger (sa)3 and (sa)4 clusters are significantly more hydrated than the smaller (sa)2 cluster or the sulfuric acid (sa)1 molecule. These findings indicate that more than five water molecules might be required to fully saturate the sulfuric acid clusters with water under atmospheric conditions. The presented methodology gives modelers a tool to take the effect of water explicitly into account in atmospheric particle formation models based on quantum chemistry.

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 CO2 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)1(MSA)1 cluster, followed by alternately adding one MSA and one MEA molecule, differing from the case of MA-MSA clusters.

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 CO2 capture. Participating in H2SO4 (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)1(SA)1 cluster, the cluster mainly grows by gradual addition of SA or PZ monomer, followed by addition of (PZ)1(SA)1 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.

Reply to the 'Comment on "Atmospheric chemistry of iodine anions: elementary reactions of I-, IO-, and IO with ozone studied in the gas-phase at 300 K using an ion trap"' by D. Britz, Phys. Chem. Chem. Phys., 2019, 21, C9CP03851E.

We reply to the comment by Dieter Britz on two recent papers in Physical Chemistry Chemical Physics. The comment presents a valuable, however, less flexible, alternative to the analysis performed in these papers and as such has no impact on any of the scientific results reported in the two publications.

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(H2O)- with ozone O3 at room temperature performed using a radio-frequency ion trap combined with a quadrupole mass spectrometer. Contrary to the oxidation reaction of the bare I- ion, the hydrated iodide I(H2O)- primarily reacts to form I- and IO2- with significant absolute reaction rate constants of 2.0 ± 0.3 × 10-10 cm3 molecule-1 s-1 and 2.5 ± 0.3 × 10-10 cm3 molecule-1 s-1 while direct pathways to IO- and IO3- are much weaker. Quantum chemical calculations indicate that in aqueous phase and for atmospherically relevant temperatures, the presence of hydrated iodides are favored over bare I- ions, thus suggesting that the chemistry of the hydrated ions is relevant for understanding and modeling atmospheric processes at the air-water interface.