Molecular-level study on the role of methanesulfonic acid in iodine oxoacid nucleation

Iodic acid (HIO 3) and iodous acid (HIO 2) have been identified as nucleating effectively by the Cosmics Leaving Outdoor Droplets (CLOUD) experiment at CERN, yet it may be hard to explain all HIO 3 -induced nucleation. Given the complexity of marine atmosphere, other precursors may be involved. Methanesulfonic acid (MSA), as a widespread precursor over oceans, has been proven to play a vital role in facilitating nucleation. However, its kinetic impacts on the synergistic nucleation of iodine oxoacids remain unclear. Hence, we investigated the MSA-involved HIO 3 –HIO 2 nucleation process at the molecular level using density functional theory (DFT) and the Atmospheric Cluster Dynamics Code (ACDC). Our results show that MSA can form stable molecular clusters with HIO 3 and HIO 2 jointed via hydrogen bonds, halogen bonds, and electrostatic attraction after proton transfer to HIO 2. Thermodynamically, the MSA-involved clustering can occur nearly without a free-energy barrier, following the HIO 2 –MSA binary and HIO 3 –HIO 2 –MSA ternary pathway. Furthermore, our results show that considering MSA will significantly enhance the calculated rate of HIO 3 –HIO 2 -based cluster formation, by up to 104 -fold in cold marine regions containing rich MSA and scarce iodine, such as the polar regions Ny-Ålesund and Marambio. Thus, the proposed more efficient HIO 3 –HIO 2 –MSA nucleation mechanism may provide theoretical evidence for explaining the frequent and intensive bursts of marine iodine particles.

Marine aerosol, which is the primary natural aerosol (O'Dowd and de Leeuw, 2007), has a significant impact on global climate, radiation balance, and even human health (Wang et al., 2010; Pöschl, 2005). New particle formation (NPF) is a primary source of marine aerosols, which proceeds via nucleation and subsequent growth (Lee et al., 2019; Zhang, 2010; Kulmala et al., 2013; Zhang et al., 2012). Moreover, nucleation, forming critical clusters at 1–2 nm from gaseous precursors, is the pivotal step affecting NPF (Zhang, 2010; Kulmala et al., 2013). However, the chemicals involved in nucleation and the underlying mechanisms remain poorly understood due to technological limitations in molecular-level analysis. Additionally, a lack of comprehensive and long-term ocean observations further hinders our knowledge of marine NPF.

Recent field studies have suggested that marine NPF events are closely related to atmospheric iodine-bearing molecules emitted by algae (Yu et al., 2019; Baccarini et al., 2020; Beck et al., 2021). During NPF events in coastal areas (e.g., Mace Head and Zhejiang; Sipilä et al., 2016; Yu et al., 2019) and polar oceans (e.g., Arctic Ocean) (Baccarini et al., 2020), the nucleation processes are mainly driven by iodic acids (HIO ${}_{\mathrm{3}}$ ). However, the self-nucleation of HIO ${}_{\mathrm{3}}$ alone cannot fully explain the observed NPF rates (Rong et al., 2020). More recently, the Cosmics Leaving Outdoor Droplets (CLOUD) experiment at CERN has found that iodous acid (HIO ${}_{\mathrm{2}}$ ) plays a key role in stabilizing HIO ${}_{\mathrm{3}}$ , enabling effective nucleation by the sequential addition of HIO ${}_{\mathrm{3}}$ followed by HIO ${}_{\mathrm{2}}$ (He et al., 2021). Further theoretical studies have uncovered the fact that the stabilizing effect of HIO ${}_{\mathrm{2}}$ on HIO ${}_{\mathrm{3}}$ stems from its role as a base in clustering (S. Zhang et al., 2022; Liu et al., 2023). Although the efficient nucleation of HIO ${}_{\mathrm{3}}$ and HIO ${}_{\mathrm{2}}$ is overall consistent with the CLOUD measurement (R. Zhang et al., 2022), this mechanism does not account for all HIO ${}_{\mathrm{3}}$ -induced nucleation in the authentic atmosphere (Ma et al., 2023). Thus, other essential precursors in the marine atmosphere might potentially affect HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ nucleation, but which ones have an effect and how remain largely unexamined.

Methanesulfonic acid (MSA), as a typical marine sulfur precursor, is widespread over oceans (Saltzman et al., 1983; Read et al., 2008; Chen et al., 2012; Yan et al., 2019) with considerable atmospheric concentrations ( ${\mathrm{10}}^{\mathrm{5}}$ – ${\mathrm{10}}^{\mathrm{8}}$  molec. cm ${}^{-\mathrm{3}}$ ) (Eisele and Tanner, 1993; Dal Maso et al., 2002; Chen et al., 2018; Yan et al., 2019). Moreover, MSA has been shown to initiate nucleation with vital atmospheric precursors, such as ammonia and amines, enhancing cluster formation (O'Dowd et al., 2002; Bork et al., 2014; Shen et al., 2019, 2020; Brean et al., 2021; Liu et al., 2022). Importantly, current evidence suggests that MSA can also form stable clusters with HIO ${}_{\mathrm{3}}$ or HIO ${}_{\mathrm{2}}$ individually, but none of the resulting binary nucleation can explain field measurements well (Ning et al., 2022; Wu et al., 2023). Despite the stabilizing effect of MSA on iodine oxoacids, whether MSA can synergistically nucleate with HIO ${}_{\mathrm{3}}$ and HIO ${}_{\mathrm{2}}$ , as well as the induced kinetic impacts on clustering, remains unknown. Furthermore, given the coexistence of MSA and HIO ${}_{\mathrm{3}}$ in different marine regions (Quéléver et al., 2022; Beck et al., 2021), along with the consistent presence of HIO ${}_{\mathrm{3}}$ and HIO ${}_{\mathrm{2}}$ as homologous substances (Sipilä et al., 2016), the importance of the HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ –MSA nucleation mechanism may differ under distinct ambient conditions, but this remains unclear.

Herein, we have systematically investigated the HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ -based nucleation involved in MSA, including $\left({\mathrm{HIO}}_{\mathrm{3}}{)}_x\right({\mathrm{HIO}}_{\mathrm{2}}{)}_y$ (MSA) ${}_z$ ( $\mathrm{1}\le x+y+z\le \mathrm{5}$ , $\mathrm{0}\le z\le \mathrm{3}$ ) clusters, by combining a quantum chemical (QC) approach with the Atmospheric Cluster Dynamics Code (ACDC) (McGrath et al., 2012). To probe the nature of cluster formation and investigate the intermolecular interactions, wave-function analysis was performed. Additionally, the Gibbs free energies of cluster formation were calculated to evaluate cluster stability. Moreover, a series of ACDC simulations were executed to delve into the influence of MSA on nucleation rates and mechanisms under varying atmospheric conditions, such as regarding precursor concentration and temperature.

To locate the low-lying isomers of $\left({\mathrm{HIO}}_{\mathrm{3}}{)}_x\right({\mathrm{HIO}}_{\mathrm{2}}{)}_y$ (MSA) ${}_z$ ( $\mathrm{1}\le x+y+z\le \mathrm{5}$ , $\mathrm{0}\le z\le \mathrm{3}$ ) clusters, the multi-step conformer search was adopted here (details in the Supplement). The resulting stable clusters with the lowest energies were identified at the $\mathit{\omega }$ B97X-D/6-311 $++$ G(3df,3pd) (for C, H, O, and S atoms) $+$ aug-cc-pVTZ-PP with ECP28MDF (for the I atom) level of theory (Francl et al., 1982; Peterson et al., 2003), and the corresponding Cartesian coordinates are collected in Table S9 in the Supplement. In addition, the structures of pure HIO ${}_{\mathrm{3}}$ , pure HIO ${}_{\mathrm{2}}$ , pure MSA, HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ , HIO ${}_{\mathrm{3}}$ –MSA, and HIO ${}_{\mathrm{2}}$ –MSA clusters in the present study have been adopted from previous studies (Rong et al., 2020; S. Zhang et al., 2022; Liu et al., 2023; Ning et al., 2022; Wu et al., 2023). All density functional theory (DFT) calculations were carried out using the Gaussian 09 package (Frisch et al., 2009), where FineGrid and tight convergence were employed. The single-point energy was calculated at the RI-CC2/aug-cc-pVTZ (for C, H, and O atoms) $+$ aug-cc-pV(T $+$ d)Z (for the S atom) $+$ aug-cc-pVTZ-PP with ECP28MDF (for the I atom) level of theory (Hättig and Weigend, 2000) by the TURBOMOLE program (Ahlrichs et al., 1989) because of its success in fitting with the experiments (Lu et al., 2020; Kürten et al., 2018; Rong et al., 2020; Almeida et al., 2013). In the present study, the Gibbs formation free energy ( $\mathrm{\Delta }{G}_{\mathrm{ref}}$ , kcal mol ${}^{-\mathrm{1}}$ ) of the HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ –MSA clusters at the reference pressure (1 atm) was calculated as

1 $$\mathrm{\Delta }{\mathrm{G}}_{\mathrm{ref}}=\mathrm{\Delta }{E}_{\mathrm{RI}-\mathrm{CC}\mathrm{2}}+\mathrm{\Delta }{G}_{\mathrm{thermal}}^{\mathit{\omega }\mathrm{B}\mathrm{97}\mathrm{X}-\mathrm{D}},$$

where $\mathrm{\Delta }{E}_{\mathrm{RI}-\mathrm{CC}\mathrm{2}}$ is the electronic contribution and $\mathrm{\Delta }{\mathrm{G}}_{\mathrm{thermal}}^{\mathit{\omega }\mathrm{B}\mathrm{97}\mathrm{X}-\mathrm{D}}$ is the thermal contribution to free energy. The $\mathrm{\Delta }{G}_{\mathrm{ref}}$ values at different temperatures ( $T=$  258–298 K) were calculated using the Shermo 2.0 code (Lu and Chen, 2021) and are collected in Table S1. Further, given the effect of the precursor's vapor pressures, the $\mathrm{\Delta }{G}_{\mathrm{ref}}$ was converted into $\mathrm{\Delta }G\left(P_{\mathrm{1}},P_{\mathrm{2}},\mathrm{...},P_n\right)$ (Vehkamäki, 2006) by Eq. (2):

2 $$\mathrm{\Delta }G\left(P_{\mathrm{1}},P_{\mathrm{2}},\mathrm{...},P_n\right)=\mathrm{\Delta }{G}_{\mathrm{ref}}-k_{\mathrm{B}}T\sum _{i=\mathrm{1}}^n{N}_i\ln \left({\displaystyle \frac{P_i}{P_{\mathrm{ref}}}}\right),$$

where $n$ is the number of components within the cluster, $k_{\mathrm{B}}$ denotes the Boltzmann constant, $T$ signifies the temperature, $N_i$ refers to the number of molecules of type $i$ in the number of components in the cluster, and $P_i$ is the partial pressure of component $i$ in the vapor phase.

Wave-function analysis was carried out using Multiwfn 3.7 (Lu and Chen, 2012) to investigate the binding nature within molecular clusters. The electrostatic potential (ESP) on the van der Waals (vdW) surface was calculated to identify active interaction sites. Specifically, the negative ESP region is electron-rich, while the positive ESP region is electron-deficient, potentially leading to mutual non-covalent interactions, such as hydrogen bonds (HBs) and halogen bonds (XBs). To further quantify the bond strength, the electron density $\mathit{\rho }\left(r\right)$ , Laplacian electron density ${\mathrm{\nabla }}^{\mathrm{2}}\mathit{\rho }\left(r\right)$ , and energy density $H\left(r\right)$ at bond critical points (BCPs) were calculated based on the atoms-in-molecules (AIM) theory (Lane et al., 2013).

To explore nucleation kinetics, the Atmospheric Cluster Dynamics Code (ACDC) (McGrath et al., 2012) was adopted here to compute the cluster formation rates, steady-state concentrations, and formation pathways by explicit solution of the birth–death equations (Eq. 3).

3 $$\begin{array}{rl}{\displaystyle \frac{\mathrm{d}{C}_i}{\mathrm{d}t}}& ={\displaystyle \frac{\mathrm{1}}{\mathrm{2}}}\sum _{j<i}{\mathit{\beta }}_{j,(i-j)}{C}_j{C}_{(i-j)}+\sum _j{\mathit{\gamma }}_{(i+j)\to i}{C}_{i+j}\\ & -\sum _j{\mathit{\beta }}_{i,j}{C}_i{C}_j-{\displaystyle \frac{\mathrm{1}}{\mathrm{2}}}\sum _{j<i}{\mathit{\gamma }}_{i\to j}{C}_i+Q_i-S_i,\end{array}$$

where the subscripts $i$ , $j$ , $i-j$ , and $i+j$ denote different clusters or monomers; $C_i$ is the concentration of cluster $i$ ; and ${\mathit{\beta }}_{i,j}$ and ${\mathit{\gamma }}_{(i+j)\to i}$ represent the cluster collision and evaporation rate coefficient, respectively. $Q_i$ and $S_i$ denote the external source and sink terms, respectively. ${\mathit{\beta }}_{i,j}$ is calculated as follows:

4 $${\mathit{\beta }}_{i,j}={\left({\displaystyle \frac{\mathrm{3}}{\mathrm{4}\mathit{\pi }}}\right)}^{\mathrm{1}/\mathrm{6}}{\left({\displaystyle \frac{\mathrm{6}{k}_{\mathrm{B}}T}{m_i}}+{\displaystyle \frac{\mathrm{6}{k}_{\mathrm{B}}T}{m_j}}\right)}^{\mathrm{1}/\mathrm{2}}{\left(V_i^{\mathrm{1}/\mathrm{3}}+V_j^{\mathrm{1}/\mathrm{3}}\right)}^{\mathrm{2}},$$

where $m_i$ and $V_i$ represent the mass and volume of cluster $i$ , respectively. $V_i=\mathrm{3}/\mathrm{4}\times \mathit{\pi }\times {\left(d_i/\mathrm{2}\right)}^{\mathrm{3}}$ , where the diameter $d_i$ of cluster $i$ is derived from the cluster volume $V_i$ calculated by Multiwfn 3.7 (Lu and Chen, 2012). ${\mathit{\gamma }}_{(i+j)\to i}$ is calculated by Eq. (5):

5 $${\mathit{\gamma }}_{(i+j)\to i}={\mathit{\beta }}_{i,j}{\displaystyle \frac{P_{\mathrm{ref}}}{k_{\mathrm{B}}T}}\exp \left({\displaystyle \frac{\mathrm{\Delta }{G}_{i+j}-\mathrm{\Delta }{G}_i-\mathrm{\Delta }{G}_j}{k_{\mathrm{B}}T}}\right),$$

where $P_{\mathrm{ref}}$ is the reference pressure at 1 atm and $\mathrm{\Delta }G$ is the formation free energy of the cluster.

In the ACDC simulations performed, we considered all possible collision and evaporation processes, including monomer–monomer, monomer–cluster, and cluster–cluster collisions, as well as decomposition of parent clusters into monomers and clusters or into two smaller clusters. Additionally, whether the clusters in the simulated system are stable depends on whether the rate of collision frequencies exceeds the total evaporation rate coefficients ( $\mathit{\beta }C/\mathrm{\Sigma }\mathit{\gamma }>\mathrm{1}$ ) (Table S4). The setting of the boundary conditions of ACDC simulations is summarized in Table S3. An uncertainty analysis was performed in this study (Table S8 and Figs. S15, S16, and S17), with details provided in the Supplement.

Here, conformational analysis was first carried out to study how MSA affects intermolecular interactions in the HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ –MSA clusters. The thermodynamic analysis was employed to assess stability of the formed clusters. To gain insights into nucleation mechanisms, a series of ACDC simulations were executed under varying atmospheric conditions.

Strong interactions among nucleation precursors are pivotal for forming stable clusters. To evaluate the binding potential of MSA with HIO ${}_{\mathrm{3}}$ and HIO ${}_{\mathrm{2}}$ , we calculated the ESP-mapped molecular vdW surface to identify interaction sites. As illustrated in Fig. 1, MSA has a positive ESP maximum ( $+$ 63.95 kcal mol ${}^{-\mathrm{1}}$ ) at the H atom of its $-$ OH group, serving as an HB donor. The iodine atoms of HIO ${}_{\mathrm{3}}$ and HIO ${}_{\mathrm{2}}$ with positive ESP maximums ( $+$ 51.90 and $+$ 45.26 kcal mol ${}^{-\mathrm{1}}$ ) can act as XB donors. Additionally, the oxygen atoms in the S $=$ O group (from MSA) and I $=$ O group (from HIO ${}_{\mathrm{3}}$ and HIO ${}_{\mathrm{2}}$ ) with strong electronegativity can act as HB or XB acceptor sites due to the lone-pair electrons. Therefore, as shown in Fig. 1d, MSA, HIO ${}_{\mathrm{3}}$ , and HIO ${}_{\mathrm{2}}$ have the potential to form clusters via intermolecular HBs and XBs.

Graph: Figure 1 The ESP-mapped molecular vdW surface of (a) MSA, (b) HIO2, (c) HIO3, and (d) (HIO3)1(HIO2)1(MSA)1. The golden and cyan dots represent the positions of maximums and minimums of ESP (units: kcal mol-1), respectively. The dashed gray arrows signify the site-to-site interaction tendencies.

As presented in Fig. 2, all the identified HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ –MSA clusters are structurally stabilized by the network of HBs (dashed blue lines) and XBs (dashed red lines). Within these clusters, the inward-facing oxygen atom and hydroxyl ( $-$ OH) group in MSA facilitates it being involved in forming more HBs and XBs compared to the HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ clusters (Fig. S1 in the Supplement). Statistically, within HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ –MSA clusters, the percentage of XBs (61 %) is higher than that of HBs (39 %). Notably, during the HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ –MSA cluster formation, HIO ${}_{\mathrm{2}}$ behaves like a base and is protonated by MSA instead of HIO ${}_{\mathrm{3}}$ , likely due to stronger acidity of MSA than HIO ${}_{\mathrm{3}}$ . After the MSA-driven proton transfer to HIO ${}_{\mathrm{2}}$ , the resulting electrostatic interactions between the formed ion pairs (CH ${}_{\mathrm{3}}$ SO ${}_{\mathrm{3}}^{-}$ –H ${}_{\mathrm{2}}$ IO ${}_{\mathrm{2}}^{+}$ ) further stabilize the clusters. Taken together, MSA can form clusters with HIO ${}_{\mathrm{3}}$ and HIO ${}_{\mathrm{2}}$ via HBs, XBs, and electrostatic attraction between ion pairs after proton transfer. Additionally, taking the $\left({\mathrm{HIO}}_{\mathrm{3}}{)}_{\mathrm{1}}\right({\mathrm{HIO}}_{\mathrm{2}}{)}_{\mathrm{3}}$ (MSA) ${}_{\mathrm{1}}$ cluster for example, there are still some potential remaining unoccupied binding sites as shown in Fig. S2. This suggests that the large-size clusters studied still have unoccupied HB and XB sites that can potentially facilitate the condensation of precursors in the atmosphere, enhancing further growth of marine aerosols.

Graph: Figure 2 The most stable configurations of the HIO3–HIO2–MSA ternary clusters identified at the ωB97XD/6-311++G(3df, 3pd) (for C, H, O, and S atoms) + aug-cc-pVTZ-PP with ECP28MDF (for I atom) level of theory. The lengths of bonds are given in ångströms (Å).

To further quantify bond strength within HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ –MSA clusters, topological analysis was performed based on the atoms-in-molecules (AIM) theory. The electron density $\mathit{\rho }\left(r\right)$ , Laplacian electron density ${\mathrm{\nabla }}^{\mathrm{2}}\mathit{\rho }\left(r\right)$ , and energy density $H\left(r\right)$ at the corresponding bond critical points (BCPs) in the studied HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ –MSA clusters were calculated and are collected in Table S2. $\mathit{\rho }\left(r\right)$ is generally positively associated with the bond strength. For the HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ –MSA clusters, $\mathit{\rho }\left(r\right)$ values at the BCPs of the HBs range from 0.0090 to 0.0869 a.u. (atomic units), exceeding the reported threshold of HB (0.002–0.040 a.u.) (Koch and Popelier, 1995; Grabowski, 2004). Moreover, the associated values of ${\mathrm{\nabla }}^{\mathrm{2}}\mathit{\rho }\left(r\right)$ at these BCPs range from 0.0310 to 0.1180 a.u., falling within the suggested range for HB (0.014–0.139 a.u.) (Koch and Popelier, 1995; Grabowski, 2004). Meanwhile, for O–I $\mathrm{\cdots }$ O XBs, the $\mathit{\rho }\left(r\right)$ and ${\mathrm{\nabla }}^{\mathrm{2}}\mathit{\rho }\left(r\right)$ values within HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ –MSA clusters range from 0.0068 to 0.1999 and 0.0288 to 0.1744 a.u., respectively. Collectively, MSA can stabilize HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ clusters via more relatively strong HBs and XBs while also protonating HIO ${}_{\mathrm{2}}$ to form ion pairs.

To explore how MSA affects HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ -based nucleation kinetics, the ACDC simulations were employed to reveal the nucleation mechanism under varying atmospheric conditions. Based on field measurements (Berresheim et al., 2002; Chen et al., 2018; Sipilä et al., 2016; Beck et al., 2021), the ranges of $\left[\mathrm{MSA}\right]$ , $\left[{\mathrm{HIO}}_{\mathrm{3}}\right]$ , and $\left[{\mathrm{HIO}}_{\mathrm{2}}\right]$ are set to be ${\mathrm{10}}^{\mathrm{6}}$ – ${\mathrm{10}}^{\mathrm{8}}$ , ${\mathrm{10}}^{\mathrm{6}}$ – ${\mathrm{10}}^{\mathrm{8}}$ , and $\mathrm{2.0}\times {\mathrm{10}}^{\mathrm{4}}$ – $\mathrm{2.0}\times {\mathrm{10}}^{\mathrm{6}}$  molec. cm ${}^{-\mathrm{3}}$ , respectively, where $\left[{\mathrm{HIO}}_{\mathrm{3}}\right]/\left[{\mathrm{HIO}}_{\mathrm{2}}\right]$ is a constant. Here, the condensation sink (CS) coefficient is set to be $\mathrm{2.0}\times {\mathrm{10}}^{-\mathrm{3}}$  s ${}^{-\mathrm{1}}$ (Dal Maso et al., 2002) and the temperature ( $T$ ) is 278 K. Under such conditions, the molecular-level nucleation pathways and the corresponding branching ratios are depicted in Fig. 3a. The detailed branch ratio is also shown at 278 K (Fig. S3) and 268 K (Fig. S4). Furthermore, to comprehend how the growth occurs thermodynamically, we herein calculated the Gibbs free energies ( $\mathrm{\Delta }G$ , Eq. 2) along the main clustering pathway at the conditions of $T=$  268–278 K, $\left[{\mathrm{HIO}}_{\mathrm{3}}\right]$   $=\mathrm{1.0}\times {\mathrm{10}}^{\mathrm{6}}$ , $\left[{\mathrm{HIO}}_{\mathrm{2}}\right]$   $=\mathrm{2.0}\times {\mathrm{10}}^{\mathrm{4}}$ , and $\left[\mathrm{MSA}\right]$   $=\mathrm{5.0}\times {\mathrm{10}}^{\mathrm{6}}$  molec. cm ${}^{-\mathrm{3}}$ (Figs. 3b and S5).

Graph: Figure 3 (a) Left: the main cluster growth pathway of the HIO3–HIO2–MSA nucleation system at T= 278 K, CS =2.0×10-3 s-1, [HIO3] =1.0×106, [HIO2] =2.0×104, and [MSA] =5.0×106 molec. cm-3. Right: the branch ratio of outward flux under varying [MSA] (106–108 molec. cm-3) and [HIO3] (106–108 molec. cm-3) values. The yellow and purple arrows (or bars) denote MSA-involved and non-MSA outward flux (or branch ratio), respectively. (b) The Gibbs free energies of cluster formation (ΔG, kcal mol-1) based on the main clustering pathway in the HIO3–HIO2–MSA nucleation system. [HIO3]/[HIO2] is a constant.

As shown in Fig. 3a, the clustering pathways, at $T=\mathrm{278}$  K, CS  $=\mathrm{2.0}\times {\mathrm{10}}^{-\mathrm{3}}$  s ${}^{-\mathrm{1}}$ , $\left[{\mathrm{HIO}}_{\mathrm{3}}\right]$   $=\mathrm{1.0}\times {\mathrm{10}}^{\mathrm{6}}$ , $[{\mathrm{HIO}}_{\mathrm{3}}$ ]  $=\mathrm{2.0}\times {\mathrm{10}}^{\mathrm{4}}$ , and [MSA]  $=\mathrm{5.0}\times {\mathrm{10}}^{\mathrm{6}}$  molec. cm ${}^{-\mathrm{3}}$ , can be categorized into two main types: (i) MSA-involved pathways, including HIO ${}_{\mathrm{2}}$ –MSA and HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ –MSA nucleation, and (ii) non-MSA pathways, primarily involving HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ nucleation. For the HIO ${}_{\mathrm{2}}$ –MSA pathway, the initial formation of $({\mathrm{HIO}}_{\mathrm{2}}{)}_{\mathrm{1}}$ (MSA) ${}_{\mathrm{1}}$ heterodimers occurs without any energy barrier (Fig. 3b). The subsequent cluster growth mainly proceeds via sequential addition of HIO ${}_{\mathrm{2}}$ or MSA monomers, partly coupled with cluster collisions. Specifically, 63 % of $({\mathrm{HIO}}_{\mathrm{2}}{)}_{\mathrm{2}}$ (MSA) ${}_{\mathrm{2}}$ results from $({\mathrm{HIO}}_{\mathrm{2}}{)}_{\mathrm{2}}$ (MSA) ${}_{\mathrm{1}}$ colliding with MSA monomers with an energy barrier of 1.00 kcal mol ${}^{-\mathrm{1}}$ , while 36 % results from a barrierless combination of two $({\mathrm{HIO}}_{\mathrm{2}}{)}_{\mathrm{1}}$ (MSA) ${}_{\mathrm{1}}$ clusters. At this point, kinetics drives growth by colliding $({\mathrm{HIO}}_{\mathrm{2}}{)}_{\mathrm{2}}$ (MSA) ${}_{\mathrm{1}}$ with MSA monomers, instead of following the lowest-energy pathway. This is because the collision frequency of $({\mathrm{HIO}}_{\mathrm{2}}{)}_{\mathrm{1}}$ (MSA) ${}_{\mathrm{1}}$ and a HIO ${}_{\mathrm{2}}$ monomer is higher, stemming from the higher $\left[{\mathrm{HIO}}_{\mathrm{2}}\right]$ . Then, the formed $({\mathrm{HIO}}_{\mathrm{2}}{)}_{\mathrm{2}}$ (MSA) ${}_{\mathrm{2}}$ further collides with a HIO ${}_{\mathrm{2}}$ monomer, yielding the stable $({\mathrm{HIO}}_{\mathrm{2}}{)}_{\mathrm{3}}$ (MSA) ${}_{\mathrm{2}}$ cluster against evaporation. As for HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ –MSA nucleation, the formation of an $\left({\mathrm{HIO}}_{\mathrm{3}}{)}_{\mathrm{1}}\right({\mathrm{HIO}}_{\mathrm{2}}{)}_{\mathrm{3}}$ (MSA) ${}_{\mathrm{1}}$ cluster arises from the collision of (i)  $({\mathrm{HIO}}_{\mathrm{2}}{)}_{\mathrm{2}}$ (MSA) ${}_{\mathrm{1}}$ with $\left({\mathrm{HIO}}_{\mathrm{3}}{)}_{\mathrm{1}}\right({\mathrm{HIO}}_{\mathrm{2}}{)}_{\mathrm{1}}$ (18 %, energy barrier: 1.00 kcal mol ${}^{-\mathrm{1}}$ ) and (ii)  $\left({\mathrm{HIO}}_{\mathrm{3}}{)}_{\mathrm{1}}\right({\mathrm{HIO}}_{\mathrm{2}}{)}_{\mathrm{2}}$ with $({\mathrm{HIO}}_{\mathrm{2}}{)}_{\mathrm{1}}$ (MSA) ${}_{\mathrm{1}}$ (62 %, energy barrier: 0.17 kcal mol ${}^{-\mathrm{1}}$ ), as well as (iii)  $\left({\mathrm{HIO}}_{\mathrm{3}}{)}_{\mathrm{1}}\right({\mathrm{HIO}}_{\mathrm{2}}{)}_{\mathrm{3}}$ binding with a MSA monomer (20 %, energy barrier: 3.80 kcal mol ${}^{-\mathrm{1}}$ ). In addition, for the non-MSA pathway marked by purple arrows, 76 % of $\left({\mathrm{HIO}}_{\mathrm{3}}{)}_{\mathrm{2}}\right({\mathrm{HIO}}_{\mathrm{2}}{)}_{\mathrm{2}}$ cluster formation arises from the collision between two $\left({\mathrm{HIO}}_{\mathrm{3}}{)}_{\mathrm{1}}\right({\mathrm{HIO}}_{\mathrm{2}}{)}_{\mathrm{1}}$ clusters, which accords closely with the barrierless pathway shown in Fig. 3b.

Overall, the MSA-involved pathways contribute to 74 % of cluster formation, while the non-MSA path accounts for only 26 %. Although the HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ –MSA growth pathway is less favorable than the HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ and HIO ${}_{\mathrm{2}}$ –MSA pathways at $T=$  278 K, it can become barrierless at a lower temperature of 268 K (Fig. S5). This result may be explained by the fact that the lower temperature results in a decrease in the evaporation rates of the HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ –MSA clusters. The detailed cluster evaporation paths and corresponding $\mathit{\gamma }$ at 268 and 278 K are collected in Tables S5 and S6. Generally, stable clusters have lower evaporation rates. According to the calculated cluster evaporation rates ( $\mathrm{\Sigma }\mathit{\gamma }$ , s ${}^{-\mathrm{1}}$ ) at 278 K (Table S7), more than 40 % of the clusters have $\mathrm{\Sigma }\mathit{\gamma }$ values less than ${\mathrm{10}}^{-\mathrm{3}}$  s ${}^{-\mathrm{1}}$ , indicating relatively high stability ( $\mathit{\beta }C/\mathrm{\Sigma }\mathit{\gamma }>\mathrm{1}$ ). Among these resulting stable clusters (see Fig. S6), the majority (85 %) contain HIO ${}_{\mathrm{2}}$ . Moreover, the concentration of these stable clusters increases gradually with time, even reaching a maximum of ${\mathrm{10}}^{\mathrm{4}}$  molec. cm ${}^{-\mathrm{3}}$ (Fig. S6). Of these stable clusters, initial $\left({\mathrm{HIO}}_{\mathrm{3}}{)}_{\mathrm{1}}\right({\mathrm{HIO}}_{\mathrm{2}}{)}_{\mathrm{1}}$ , $({\mathrm{HIO}}_{\mathrm{2}}{)}_{\mathrm{2}}$ , and $\left(\mathrm{MSA}{)}_{\mathrm{1}}\right({\mathrm{HIO}}_{\mathrm{2}}{)}_{\mathrm{1}}$ dimers form rapidly, and at ​​​​​​​ $t\approx \mathrm{1}$  s, heterotrimers $\left({\mathrm{HIO}}_{\mathrm{3}}{)}_{\mathrm{1}}\right({\mathrm{HIO}}_{\mathrm{2}}{)}_{\mathrm{2}}$ and $\left(\mathrm{MSA}{)}_{\mathrm{1}}\right({\mathrm{HIO}}_{\mathrm{2}}{)}_{\mathrm{2}}$ begin to form, after which the larger-sized clusters also form. This time-dependent evidence suggests that MSA is involved in the whole clustering process, from the initial formation of smaller clusters to the large-sized nucleated clusters that potentially further grow. Taken together, these findings highlight the direct and significant involvement of MSA in HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ nucleation, facilitating cluster formation.

As presented on the right of Fig. 3a, the contribution of different clustering pathways to the outward flux varies with precursor concentrations. With increasing [MSA] from ${\mathrm{10}}^{\mathrm{6}}$ to ${\mathrm{10}}^{\mathrm{8}}$  molec. cm ${}^{-\mathrm{3}}$ , the contribution of MSA-involved pathways rises from 1 % to 99 % during nucleation. At the median [MSA] of ${\mathrm{10}}^{\mathrm{7}}$  molec. cm ${}^{-\mathrm{3}}$ , the MSA-involved pathway contributes 86 %, whereas the non-MSA pathway accounts for just 14 %. In contrast, the ratio of MSA-involved pathways decreases (from 86 % to 0 %) with increasing concentrations of iodine oxoacids. At higher $\left[{\mathrm{HIO}}_{\mathrm{3}}\right]$ values of ${\mathrm{10}}^{\mathrm{8}}$  molec. cm ${}^{-\mathrm{3}}$ , the HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ pathway dominates nucleation. Predictably, the kinetic impact of MSA on HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ nucleation is more pronounced in marine regions with richer MSA away from iodine sources.

Guided by the clustering pathway analysis, MSA has shown its potential to participate in the HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ -based nucleation, but its detailed impacts on cluster formation rates ( $J$ , cm ${}^{-\mathrm{3}}$  s ${}^{-\mathrm{1}}$ ) remain uncertain. The influence of MSA on $J$ under different atmospheric conditions is systematically analyzed below.

Figure 4 presents the simulated $J$ of the HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ –MSA (red bar) and HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ system (gray bar) against the varying temperatures ( $T=$  258–298 K) at CS  $=\mathrm{2.0}\times {\mathrm{10}}^{-\mathrm{3}}$  s ${}^{-\mathrm{1}}$ , $\left[{\mathrm{HIO}}_{\mathrm{3}}\right]$   $=\mathrm{1.0}\times {\mathrm{10}}^{\mathrm{7}}$ , $\left[{\mathrm{HIO}}_{\mathrm{2}}\right]$   $=\mathrm{2.0}\times {\mathrm{10}}^{\mathrm{5}}$ , and $\left[\mathrm{MSA}\right]$   $=\mathrm{1.0}\times {\mathrm{10}}^{\mathrm{7}}$  molec. cm ${}^{-\mathrm{3}}$ . Clearly, $J$ (HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ –MSA) is consistently higher than $J$ (HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ ), highlighting the enhancement effect of MSA on HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ -based clustering under common atmospheric temperatures. Specifically, both $J$ (HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ –MSA) and $J$ (HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ ) are negatively dependent on $T$ due to reduced cluster evaporation caused by low $T$ . As a result, MSA could promote nucleation with higher $J$ , especially in the colder regions, such as polar oceans.

Graph: Figure 4 Simulated cluster formation rates J (cm-3 s-1) against varying atmospheric temperatures: T= 258–298 K, CS =2.0×10-3 s-1, [HIO3] =1.0×107, [HIO2] =2.0×105, and [MSA] =1.0×107 molec. cm-3.

It is worth noting that apart from atmospheric temperature, precursor concentrations might also vary regionally or seasonally, which can further affect nucleation. So, to comprehensively reveal the effect of MSA, here we defined and calculated MSA-driven enhancement factor $R$ (Eq. 6) under varying concentrations of MSA, HIO ${}_{\mathrm{3}}$ , and HIO ${}_{\mathrm{2}}$ (i.e., [MSA], [HIO ${}_{\mathrm{3}}$ ], and [HIO ${}_{\mathrm{2}}$ ], units: molec. cm ${}^{-\mathrm{3}}$ ).

6 $$\begin{array}{rl}R& ={\displaystyle \frac{J({\mathrm{HIO}}_{\mathrm{3}}-{\mathrm{HIO}}_{\mathrm{2}}-\mathrm{MSA})}{J({\mathrm{HIO}}_{\mathrm{3}}-{\mathrm{HIO}}_{\mathrm{2}})}}\\ & ={\displaystyle \frac{J\left(\left[{\mathrm{HIO}}_{\mathrm{3}}\right]=x,\left[{\mathrm{HIO}}_{\mathrm{2}}\right]=y,\left[\mathrm{MSA}\right]=z\right)}{J\left(\left[{\mathrm{HIO}}_{\mathrm{3}}\right]=x,\left[{\mathrm{HIO}}_{\mathrm{2}}\right]=y\right)}},\end{array}$$

where $J$ (HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ –MSA) and $J$ (HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ ) represent the cluster formation rate of the HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ –MSA and HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ system, respectively. $x$ , $y$ , and $z$ are the adopted [HIO ${}_{\mathrm{3}}$ ], [HIO ${}_{\mathrm{2}}$ ], and [MSA], respectively.

Graph: Figure 5 Enhancement strength R of MSA acting on cluster formation rates at varying precursor concentrations: [HIO3] =106–108 and [HIO2] =2.0×104–2.0×106 molec. cm-3 for (a) [MSA] =1.0×106 molec. cm-3, (b) [MSA] =1.0×107 molec. cm-3, and (c) [MSA] =1.0×108 molec. cm-3; T= 278 K and CS =2.0×10-3 s-1.

As seen from Fig. 5, with [MSA] ranging from ${\mathrm{10}}^{\mathrm{6}}$ to ${\mathrm{10}}^{\mathrm{8}}$  molec. cm ${}^{-\mathrm{3}}$ (Fig. 5a–c), the maximum $R$ increases 5- to ${\mathrm{10}}^{\mathrm{4}}$ -fold, which is primarily due to MSA-mediated synergistic nucleation with HIO ${}_{\mathrm{3}}$ and HIO ${}_{\mathrm{2}}$ (recalling Sect. 3.2). Even at a median [MSA] of $\mathrm{1.0}\times {\mathrm{10}}^{\mathrm{7}}$  molec. cm ${}^{-\mathrm{3}}$ , the resulting $R$ can increase approximately ${\mathrm{10}}^{\mathrm{2}}$ -fold. In contrast, $R$ is decayed under conditions of higher [HIO ${}_{\mathrm{3}}$ ] and [HIO ${}_{\mathrm{2}}$ ]. Furthermore, under the conditions with lower $\left[{\mathrm{HIO}}_{\mathrm{3}}\right]/\left[{\mathrm{HIO}}_{\mathrm{2}}\right]$ , where $R$ is higher, the contribution of MSA nucleating with HIO ${}_{\mathrm{2}}$ increases due to the relative scarcity of HIO ${}_{\mathrm{3}}$ . Conversely, $R$ decreases at higher $\left[{\mathrm{HIO}}_{\mathrm{3}}\right]/\left[{\mathrm{HIO}}_{\mathrm{2}}\right]$ ; i.e., the impacts of MSA decrease. That is, the enhancing effect of MSA on $J$ is limited in near-iodine source regions. Naturally, in regions with sparser iodine, the promoting effect of MSA is significant. However, the atmospheric [HIO ${}_{\mathrm{3}}$ ] ranges widely from ${\mathrm{10}}^{\mathrm{6}}$ to ${\mathrm{10}}^{\mathrm{8}}$  molec. cm ${}^{-\mathrm{3}}$ . When [HIO ${}_{\mathrm{3}}$ ] is comparable to or higher than [MSA], the HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ pathway contributes more, and the $R$ of MSA decreases with the rising [HIO ${}_{\mathrm{3}}$ ]. It is worth noting that when [HIO ${}_{\mathrm{3}}$ ] is comparable to [MSA], the $R$ of MSA is greater than 2, as the contribution of MSA to clustering includes not only the direct formation of HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ –MSA clusters ( $\sim$  20 %), but also its "catalysis" role in facilitating the formation of initial HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ clusters (Fig. S7). To sum up, MSA can promote nucleation, particularly in marine regions characterized by lower $T$ and lower [HIO ${}_{\mathrm{3}}$ ] and [HIO ${}_{\mathrm{2}}$ ]. In addition, we also considered the conditions in relatively polluted (CS  $=\mathrm{1.0}\times {\mathrm{10}}^{-\mathrm{2}}$  s ${}^{-\mathrm{1}}$ ) and clean (CS  $=\mathrm{1.0}\times {\mathrm{10}}^{-\mathrm{4}}$  s ${}^{-\mathrm{1}}$ ) environments and found that, similarly to the environment with a CS value of $\mathrm{2.0}\times {\mathrm{10}}^{-\mathrm{3}}$  s ${}^{-\mathrm{1}}$ , MSA exhibits a significant promoting effect on iodine particle formation (Figs. S8–S11). Furthermore, the effect of HIO ${}_{\mathrm{2}}$ addition on the whole nucleation system was considered, as it is not only a rate-limiting step for cluster formation, leading to a significant increase in $J$ (HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ –MSA) compared to $J$ (HIO ${}_{\mathrm{3}}$ –MSA) (Fig. S12), but also thermodynamically favorable due to the HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ –MSA path being almost barrierless (1.24 kcal mol ${}^{-\mathrm{1}}$ ) compared to the HIO ${}_{\mathrm{3}}$ –MSA pathway (Fig. S13).

To further assess the atmospheric implications of the proposed HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ –MSA nucleation, we simulated $J$ , shown in Fig. 6, based on the ambient conditions of the typical polar regions (e.g., Ny-Ålesund and Marambio) and the mid-latitude marine regions (e.g., Mace Head and Réunion). Subsequently, we compared these simulation results with observed nucleation rates, and the definition of the cluster formation rate is detailed in the Supplement. As shown in Fig. 6a, the $J$ (HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ –MSA) simulated at $T=$  268 K, CS  $=\mathrm{4.0}\times {\mathrm{10}}^{-\mathrm{4}}$  s ${}^{-\mathrm{1}}$ , $\left[{\mathrm{HIO}}_{\mathrm{3}}\right]$   $={\mathrm{10}}^{\mathrm{5}}$ – ${\mathrm{10}}^{\mathrm{6}}$ , $\left[{\mathrm{HIO}}_{\mathrm{2}}\right]$   $=\mathrm{2.0}\times {\mathrm{10}}^{\mathrm{3}}$ – $\mathrm{2.0}\times {\mathrm{10}}^{\mathrm{4}}$ , and [MSA]  $={\mathrm{10}}^{\mathrm{6}}$ – ${\mathrm{10}}^{\mathrm{8}}$  molec. cm ${}^{-\mathrm{3}}$ was compared with field observations in coastal Ny-Ålesund (Beck et al., 2021; He et al., 2021). Both $J$ (HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ –MSA) and $J$ (HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ ) increase with the rising [HIO ${}_{\mathrm{3}}$ ] and [HIO ${}_{\mathrm{2}}$ ]. Importantly, the addition of MSA effectively promotes $J$ to a higher level (orange area), aligning with most field measurements ( $\mathrm{1.2}\times {\mathrm{10}}^{-\mathrm{3}}$ – $\mathrm{1.2}\times {\mathrm{10}}^{-\mathrm{1}}$  cm ${}^{-\mathrm{3}}$  s ${}^{-\mathrm{1}}$ , gray lines) (Beck et al., 2021). Even when [MSA] is as low as $\mathrm{1.0}\times {\mathrm{10}}^{\mathrm{6}}$  molec. cm ${}^{-\mathrm{3}}$ (the orange line below), the $J$ (HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ –MSA) can be 1 order of magnitude higher than the observed $J$ of $\mathrm{1.2}\times {\mathrm{10}}^{-\mathrm{3}}$  cm ${}^{-\mathrm{3}}$  s ${}^{-\mathrm{1}}$ (the gray line below). Moreover, the simulated $J$ in Fig. 6b was obtained under the conditions of coastal Marambio, Antarctica: $T=$  273 K, CS  $=\mathrm{1.0}\times {\mathrm{10}}^{-\mathrm{4}}$  s ${}^{-\mathrm{1}}$ , $\left[{\mathrm{HIO}}_{\mathrm{3}}\right]$   $={\mathrm{10}}^{\mathrm{5}}$ – ${\mathrm{10}}^{\mathrm{6}}$ , $\left[{\mathrm{HIO}}_{\mathrm{2}}\right]$   $=\mathrm{2.0}\times {\mathrm{10}}^{\mathrm{3}}$ – $\mathrm{2.0}\times {\mathrm{10}}^{\mathrm{4}}$ , and [MSA]  $={\mathrm{10}}^{\mathrm{6}}$ – ${\mathrm{10}}^{\mathrm{7}}$  molec. cm ${}^{-\mathrm{3}}$ (Quéléver et al., 2022; He et al., 2021). Compared to $J$ (HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ ), the MSA-enhanced $J$ (HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ –MSA) is overall higher, better fitting with the field observations ( $\mathrm{1.2}\times {\mathrm{10}}^{-\mathrm{1}}$ – $\mathrm{2.4}\times {\mathrm{10}}^{\mathrm{1}}$  cm ${}^{-\mathrm{3}}$  s ${}^{-\mathrm{1}}$ , gray line) (Quéléver et al., 2022). These findings imply that MSA potentially plays a vital role in cold polar oceanic regions, particularly with higher [MSA] during NPF events.

Graph: Figure 6 Comparison with the simulated cluster formation rates (J, cm-3 s-1) and field observations under the ambient conditions of (a) Ny-Ålesund (T= 268 K; CS =4.0×10-4 s-1; and [HIO3] =105–106, [HIO2] =2.0×103–2.0×104, and [MSA] =106–108 molec. cm-3) and (b) Marambio (T= 273 K; CS =1.0×10-4 s-1; and [HIO3] =105–106, [HIO2] =2.0×103–2.0×104, and [MSA] =106–107 molec. cm-3). The orange area, blue line, and gray line represent J(HIO3–HIO2–MSA), J(HIO3–HIO2), and J(field observation), respectively. [HIO3]/[HIO2] is a constant.

In addition, the influence of HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ –MSA nucleation over the relatively warmer mid-latitude marine areas, such as Mace Head and Réunion, was investigated here (Fig. S14). We found that $J$ (HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ –MSA) is slightly higher than $J$ (HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ ), especially in regions with high concentrations of iodine oxoacids (e.g., Mace Head), showing a relatively limited enhancement effect of MSA on nucleation. Based on the simulated $J$ (HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ ) ( $\sim$   $\mathrm{1.2}\times {\mathrm{10}}^{\mathrm{4}}$  cm ${}^{-\mathrm{3}}$  s ${}^{-\mathrm{1}}$ ), iodine nucleation can almost completely explain the NPF events of Mace Head (Fig. S14a), which provides potential theoretical evidence for explaining the previous findings (Sipilä et al., 2016).

Overall, in the mid-latitude oceans, especially near iodine sources like Mace Head, MSA may have a limited enhancement effect on nucleation. In this case, the clustering process is dominated by abundant iodine oxoacids, while in the colder polar regions, particularly those with higher [MSA] like Marambio, MSA indeed significantly facilitates HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ nucleation, suggesting a vital role in polar NPF.

The present study systematically investigates the HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ -based nucleation process enhanced by MSA at the molecular level by QC calculations and ACDC simulations. The results indicate that MSA can stabilize HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ -based clusters by building the intricate networks with more HBs and XBs. During clustering, MSA replaces HIO ${}_{\mathrm{3}}$ in protonating HIO ${}_{\mathrm{2}}$ to form ion pairs, resulting in relatively strong electrostatic attractions. In addition, thermodynamic analyses suggest that MSA-involved clustering is nearly barrierless. Compared to the HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ system reported previously, the MSA-involved synergistic nucleation with HIO ${}_{\mathrm{3}}$ and HIO ${}_{\mathrm{2}}$ proceeds more efficiently, through two additional clustering pathways: (i) HIO ${}_{\mathrm{2}}$ –MSA binary and (ii) HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ –MSA ternary pathways. Moreover, the resulting enhancement effect of MSA on nucleation is more substantial in colder regions, especially with richer MSA, but weaker in the environments near iodine sources. Further comparison with field observations indicates that HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ –MSA synergistic nucleation plays a limited role in mid-latitude ocean regions, particularly in regions with abundant iodine (e.g., Mace Head), but a potential role in colder polar regions (e.g., Ny-Ålesund and Marambio).

This study highlights the essential enhancing role of MSA in iodine oxoacids nucleation, and the HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ –MSA synergistic nucleation may help to explain the observation of abundant iodine particles during marine NPF events. Given the complex oceanic atmosphere, other potential nucleation precursors beyond MSA, such as sulfuric acid and amines, may also affect the HIO ${}_{\mathrm{3}}$ –HIO ${}_{\mathrm{2}}$ nucleation process and further contribute to the formation of marine iodine particles, deserving future investigations.

The supplement related to this article is available online at: https://doi.org/10.5194/acp-24-3989-2024-supplement.

We acknowledge the National Supercomputing Center in Shenzhen for providing the computational resources and the TURBOMOLE program.