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

Atmospheric Bases-Enhanced Iodic Acid Nucleation: Altitude-Dependent Characteristics and Molecular Mechanisms

Iodic acid (IA), the key driver of marine aerosols, is widely detected within the gas and particle phases in the marine boundary layer (MBL) and even the free troposphere (FT). Although atmospheric bases like dimethylamine (DMA) and ammonia (NH3) can enhance IA particles formation, their different efficiencies and spatial distributions make the dominant base-stabilization mechanisms of forming IA particles unclear. Herein, we investigated the IA-DMA-NH3 nucleation system through quantum chemical calculations at the DLPNO-CCSD(T)/aug-cc-pVTZ(-PP)//ωB97X-D/6-311++G(3df,3pd) + aug-cc-pVTZ-PP level of theory and cluster dynamics simulations. We provide molecular-level evidence that DMA and NH3 can jointly stabilize the IA clusters. The formation rates of IA clusters initially decline before rising from the MBL to the FT, owing to variations in mechanism. In the MBL, IA-DMA nucleation predominates, while the contribution of IA-DMA-NH3 synergistic nucleation cannot be overlooked in polar and NH3-polluted regions. In the lower FT, IA-DMA-NH3 nucleation prevails, whereas in the upper FT, IA-NH3 nucleation dominates. The efficiency of IA-DMA-NH3 nucleation is comparable to that of IA-iodous acid nucleation in the MBL and sulfuric acid-NH3 nucleation in the FT. Hence, the IA-DMA-NH3 mechanism holds promise for revealing the missing sources of tropospheric IA particles.

Overlooked significance of iodic acid in new particle formation in the continental atmosphere

New particle formation (NPF) substantially affects the global radiation balance and climate. Iodic acid (IA) is a key marine NPF driver that recently has also been detected inland. However, its impact on continental particle nucleation remains unclear. Here, we provide molecular-level evidence that IA greatly facilitates clustering of two typical land-based nucleating precursors: dimethylamine (DMA) and sulfuric acid (SA), thereby enhancing particle nucleation. Incorporating this mechanism into an atmospheric chemical transport model, we show that IA-induced enhancement could realize an increase of over 20% in the SA-DMA nucleation rate in iodine-rich regions of China. With declining anthropogenic pollution driven by carbon neutrality and clean air policies in China, IA could enhance nucleation rates by 1.5 to 50 times by 2060. Our results demonstrate the overlooked key role of IA in continental NPF nucleation and highlight the necessity for considering synergistic SA-IA-DMA nucleation in atmospheric modeling for correct representation of the climatic impacts of aerosols.

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

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

Shapeshifting Nucleophiles HO-(NH3)n React with Methyl Chloride

The microsolvated anions HO-(NH3)n were found to induce new nucleophile NH2-(H2O)(NH3)n-1 via intramolecular proton transfer. Hence, the ion-molecule nucleophilic substitution (SN2) reaction between CH3Cl and these shapeshifting nucleophiles lead to both the HO- path and NH2- path, meaning that the respective attacking nucleophile is HO- or NH2-. The CCSD(T) level of calculation was performed to characterize the potential energy surfaces. Calculations indicate that the HO- species are lower in energy than the NH2- species, and the SN2 reaction barriers are lower for the HO- path than the NH2--path. Incremental solvation increases the barrier for both paths. Comparison between HO-(NH3)n and HOO-(NH3)n confirmed the existence of an α-effect under microsolvated conditions. Comparison between HO-(NH3)n and HO-(H2O)n indicated that the more polarized H2O stabilizes the nucleophiles more than NH3, and thus, the hydrated systems have higher SN2 reaction barriers. The aforementioned barrier changes can be explained by the differential stabilization of the nucleophile and HOMO levels upon solvation, thus affecting the HOMO-LUMO interaction between the nucleophile and substrate. For the same kind of nucleophilic attacking atom, O or N, the reaction barrier has a good linear correlation with the HOMO level of the nucleophiles. Hence, the HOMO level or the binding energy of microsolvated nucleophiles is a good indicator to evaluate the order of barrier heights. This work expands our understanding of the microsolvation effect on prototype SN2 reactions beyond the water solvent.

Atmospheric fate of typical liquid crystal monomers in the tropospheric gas, liquid, and granular phases

Mineral aerosol particles significantly impact environmental risk prediction of liquid crystal monomers (LCMs). In this work, we investigated the reaction mechanisms and kinetics of three typical LCMs (4-cyano-3,5-difluorophenyl 4-ethylbenzoate (CEB-2F), 4-cyano-3-fluorophenyl 4-ethylbenzoate (CEB-F), and 4-cyanophenyl 4-ethylbenzoate (CEB)) with ozone (O3) in the atmospheric gas, liquid, and particle phases employing density functional theory (DFT). Here, O3 is prone to add to the benzene ring without F atom(s) in the selected LCMs. The ozonolysis products are aldehydes, carboxylic acids, epoxides, and unsaturated hydrocarbons containing aromatic rings. Those products undergo secondary ozonolysis to generate small molecular compounds such as glyoxal, which is beneficial for generating secondary organic aerosol (SOA). Titanium dioxide (TiO2), an essential component of mineral aerosol particles, has good adsorption properties for LCMs; however, it slightly reduces the reactivity with O3. At 298 K, the reaction rate constant of the selected LCMs reacting with O3 in the gas and atmospheric liquid phases is (2.74‒5.53) × 10-24 cm3/(mol·sec) and 5.58 × 10-3‒39.1 L/(mol·sec), while CEB-2F reacting with O3 on (TiO2)6 cluster is 1.84 × 10-24 cm3/(mol·sec). The existence of TiO2 clusters increases the persistence and long-distance transportability of LCMs, which enlarges the contaminated area of LCMs.

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.

The aomogeneous and heterogeneous oxidation of organophosphate esters (OPEs) in the atmosphere: Take diphenyl phosphate (DPhP) as an example

Organophosphate esters (OPEs) are widely detected in the atmosphere. However, the atmospheric oxidative degradation mechanism of OPEs has not been closely examined. This work took density functional theory (DFT) to investigate the tropospheric ozonolysis of organophosphates, represented by diphenyl phosphate (DPhP), including adsorption mechanisms on the surface of titanium dioxide (TiO₂) mineral aerosols and oxidation reaction of hydroxyl groups (·OH) after photolysis. Besides, the reaction mechanism, reaction kinetics, adsorption mechanism, and ecotoxicity evaluation of the transformation products were also studied. At 298 K, the total reaction rate constants k_O3, k_OH, k_TiO2-O3, and k_TiO2-OH are 5.72 × 10^-15 cm³ molecule^-1 s^-1, 1.68 × 10^-13 cm³ molecule^-1 s^-1, 1.91 × 10^-23 cm³ molecule^-1 s^-1, and 2.30 × 10^-10 cm³ molecule^-1 s^-1. The atmospheric lifetime of DPhP ozonolysis in the near-surface troposphere is 4 min, much lower than that of hydroxyl radicals (·OH). Besides, the lower the altitude is, the stronger the oxidation is. The TiO₂ clusters carry DPhP promoting ·OH oxidation but inhibiting ozonolysis of DPhP. Finally, the main transformation products of this process are glyoxal, malealdehyde, aromatic aldehydes, etc., which are still ecotoxic. The findings shed new light on the atmospheric governance of OPEs.

Nucleophilic substitution reactions of microsolvated hydroperoxide anion HOO-(NH3)n with methyl chloride and comparison between ammonia and water as the solvent

Similar to microhydrated hydroperoxide anion HOO^-(H₂O)_n, the HOO^-(NH₃)_n=1-3 anion can induce alternative nucleophiles by proton transfer (PT) from the solvent molecule NH₃. The PT-induced species NH₂^-(H₂O₂)(NH₃)_n-1 is higher in energy than HOO^-(NH₃)_n, obeying the proton affinity (PA) prediction that HOO^- has a higher PA than NH₂^-. The potential energy profile of HOO^-(NH₃)_n reacting with CH₃Cl shows that the transition states of the traditional HOO^--S_N2 pathway are ∼10 kcal mol^-1 lower in energy than those of the PT-induced NH₂^--S_N2 pathway, indicating the latter path is unlikely to compete. The differential solvation energy for reactants and transition states with incremental solvation increases the barrier height of both HOO^--/NH₂^--S_N2 pathways and makes the transition structures more product-like. For HOO^-(sol)_n + CH₃Cl → CH₃OOH + Cl^-(sol)_n reactions, the barrier heights for sol = H₂O are higher than those for sol = NH₃, because H₂O is more polar than NH₃, and the electrostatic interaction is strengthened, hence H₂O molecules stabilize the microsolvated nucleophiles more. In addition, because the H₂O molecule is a better proton donor than the NH₃ molecule, the PT-induced HO^-S_N2 pathway is more likely to compete with the HOO^-S_N2 pathway. The HOMO level of nucleophiles, which negatively correlates with the S_N2 barrier heights, is found to be a good descriptor to predict the S_N2 barrier height of a microsolvated system with the same attacking nucleophile. This work adds to our understanding of the differential solvent effect on the prototype ion-molecule S_N2 reactions.

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

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

The gas-phase formation mechanism of iodic acid as an atmospheric aerosol source

Iodine is a reactive trace element in atmospheric chemistry that destroys ozone and nucleates particles. Iodine emissions have tripled since 1950 and are projected to keep increasing with rising O₃ surface concentrations. Although iodic acid (HIO₃) is widespread and forms particles more efficiently than sulfuric acid, its gas-phase formation mechanism remains unresolved. Here, in CLOUD atmospheric simulation chamber experiments that generate iodine radicals at atmospherically relevant rates, we show that iodooxy hypoiodite, IOIO, is efficiently converted into HIO₃ via reactions (R1) IOIO + O₃ → IOIO₄ and (R2) IOIO₄ + H₂O → HIO₃ + HOI + ⁽¹⁾O₂. The laboratory-derived reaction rate coefficients are corroborated by theory and shown to explain field observations of daytime HIO₃ in the remote lower free troposphere. The mechanism provides a missing link between iodine sources and particle formation. Because particulate iodate is readily reduced, recycling iodine back into the gas phase, our results suggest a catalytic role of iodine in aerosol formation.

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.

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.

Iodous acid - a more efficient nucleation precursor than iodic acid

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

Insights into the Chemistry of Iodine New Particle Formation: The Role of Iodine Oxides and the Source of Iodic Acid

Iodine chemistry is an important driver of new particle formation in the marine and polar boundary layers. There are, however, conflicting views about how iodine gas-to-particle conversion proceeds. Laboratory studies indicate that the photooxidation of iodine produces iodine oxides (I_xO_y), which are well-known particle precursors. By contrast, nitrate anion chemical ionization mass spectrometry (CIMS) observations in field and environmental chamber studies have been interpreted as evidence of a dominant role of iodic acid (HIO₃) in iodine-driven particle formation. Here, we report flow tube laboratory experiments that solve these discrepancies by showing that both I_xO_y and HIO₃ are involved in atmospheric new particle formation. I₂O_y molecules (y = 2, 3, and 4) react with nitrate core ions to generate mass spectra similar to those obtained by CIMS, including the iodate anion. Iodine pentoxide (I₂O₅) produced by photolysis of higher-order I_xO_y is hydrolyzed, likely by the water dimer, to yield HIO₃, which also contributes to the iodate anion signal. We estimate that ∼50% of the iodate anion signals observed by nitrate CIMS under atmospheric water vapor concentrations originate from I₂O_y. Under such conditions, iodine-containing clusters and particles are formed by aggregation of I₂O_y and HIO₃, while under dry laboratory conditions, particle formation is driven exclusively by I₂O_y. An updated mechanism for iodine gas-to-particle conversion is provided. Furthermore, we propose that a key iodine reservoir species such as iodine nitrate, which we observe as a product of the reaction between iodine oxides and the nitrate anion, can also be detected by CIMS in the atmosphere.

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.

Potential Application of Machine-Learning-Based Quantum Chemical Methods in Environmental Chemistry

It is an important topic in environmental sciences to understand the behavior and toxicology of chemical pollutants. Quantum chemical methodologies have served as useful tools for probing behavior and toxicology of chemical pollutants in recent decades. In recent years, machine learning (ML) techniques have brought revolutionary developments to the field of quantum chemistry, which may be beneficial for investigating environmental behavior and toxicology of chemical pollutants. However, the ML-based quantum chemical methods (ML-QCMs) have only scarcely been used in environmental chemical studies so far. To promote applications of the promising methods, this Perspective summarizes recent progress in the ML-QCMs and focuses on their potential applications in environmental chemical studies that could hardly be achieved by the conventional quantum chemical methods. Potential applications and challenges of the ML-QCMs in predicting degradation networks of chemical pollutants, searching global minima for atmospheric nanoclusters, discovering heterogeneous or photochemical transformation pathways of pollutants, as well as predicting environmentally relevant end points with wave functions as descriptors are introduced and discussed.

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.