Sulfuric acid (SA), together with base molecules such as ammonia (AM), methylamine (MA), and dimethylamine (DMA), is known to play a central role in atmospheric new particle formation (NPF). NPF occurs through gas-to-particle conversion via the formation and growth of molecular clusters. While previous studies have demonstrated that mixtures of bases can strongly enhance nucleation rates, the influence of multiple acidic species for larger cluster stability and growth remains less explored.
In this work, we investigate the role of mixed-acid systems in atmospheric cluster formation using quantum chemical calculations, assisted by machine learning. Cluster structures containing SA, methanesulfonic acid (MSA), and atmospherically relevant bases (AM, MA, and DMA), with compositions up to 10 acid–base pairs, were generated through extensive configurational sampling using ABCluster and metadynamics simulations with CREST. The resulting structures were subsequently optimized at the B97-3c level of theory, while a PaiNN machine-learning model was used to accelerate the calculation.
Our results show that, in contrast to previously reported base synergy, the acid synergy between SA and MSA is weak and highly system dependent. SA consistently dominates the thermodynamic stability of the smallest clusters, and MSA-only acid–base interactions are insufficient to explain efficient initial particle formation. In particular, for systems involving DMA, the strong SA–DMA interaction governs the cluster energetics, with little contribution of MSA to the stability. However, MSA can influence cluster stability at larger sizes in systems involving weaker bases such as AM and MA, especially under conditions where the relative abundance of MSA is high. These findings indicate that MSA acts not as a primary nucleating acid but rather as a secondary species that participates in the early growth of clusters.
Overall, this work highlights that the roles of different atmospheric acids in NPF are fundamentally distinct: SA controls the initial nucleation step, whereas MSA may enhance subsequent cluster growth under specific atmospheric conditions.
European Research Council101040353Introduction
Atmospheric aerosol particles play a central role in air quality and climate. Aerosols originate from both primary emissions and secondary formation processes. While primary aerosols are emitted directly into the atmosphere, secondary aerosols are formed through gas-to-particle conversion processes that can produce freshly nucleated particles (FNPs). The initial step of this process involves the formation of molecular clusters stabilized by strong hydrogen-bonding interactions between atmospheric vapors . Clusters that are sufficiently stable against evaporation can grow through condensation and coagulation processes, eventually forming particles with diameters exceeding roughly 2 nm.
Fine particles (<1µm) and ultrafine particles (<100 nm) are particularly important due to their impacts on both human health and the Earth’s radiative balance. Fine particulate matter is one of the leading contributors to air-pollution-related mortality worldwide . In addition to their health effects, aerosols influence climate directly by scattering and absorbing solar radiation and indirectly by acting as cloud condensation nuclei (CCN), thereby affecting cloud formation and properties . Despite their importance, aerosol–cloud interactions remain one of the largest uncertainties in current climate models .
Inorganic acids and bases are widely recognized as key contributors in the initial steps of cluster formation. Sulfuric acid (SA), together with bases such as ammonia (AM) and amines, constitutes one of the most important nucleation pathways in many atmospheric environments . Other atmospheric components may further influence these processes, including ions produced by galactic cosmic rays and highly oxygenated organic molecules (HOMs) . Understanding new particle formation (NPF) therefore requires detailed knowledge of both the chemical composition of clusters and the concentrations of their precursor vapors.
Experimental characterization of clusters smaller than 2 nm remains extremely challenging. Conventional condensation particle counters (CPCs) have detection limits of approximately 2–3 nm, which prevents direct observation of the smallest clusters . Particle size magnifiers (PSMs) can go down to ∼ 1.5 nm, but they provide limited chemical information about cluster chemical composition . Techniques such as chemical ionization atmospheric pressure interface time-of-flight (CI-APi-TOF) mass spectrometer have enabled detailed investigation of cluster composition . However, they can easily alter cluster structures during ionization due to fragmentation effects . Moreover, these measurements are typically limited to clusters containing up to roughly 10 acid–base pairs (20 molecules) . As a result, the chemical composition and stability of clusters in the size range between approximately 1.0 and 2.0 nm remain poorly understood. This size regime is particularly important because it corresponds to the transition from molecular clusters to stable aerosol particles . Computational approaches have therefore become essential for bridging this gap. For example, introduced a framework that enables systematic exploration of atmospheric cluster formation from individual molecules to clusters approaching 2 nm in size. Such methods allow detailed investigation of cluster thermodynamics, intermolecular interactions, and growth pathways, providing insights into the mechanisms governing the formation of atmospheric freshly nucleated particles (FNPs).
AM and amines are emitted from a wide range of anthropogenic and natural sources, including agriculture, industrial activities, combustion processes, and marine environments . Among these, AM, methylamine (MA), and dimethylamine (DMA) play a particularly important role in enhancing SA-driven nucleation, significantly increasing particle formation rates compared to pure SA or SA–water systems . This enhancement arises from proton transfer reactions between acids and bases, leading to salt formation and increased cluster stability. Numerous quantum chemical studies have confirmed the stabilizing effect of AM , MA , and DMA in stabilizing initial SA clusters.
Methanesulfonic acid (MSA), which is frequently observed in atmospheric aerosol particles, has been proposed as an important contributor to NPF alongside SA . MSA is primarily produced through the oxidation of organosulfur compounds, often concurrently with SO2 (the main precursor of SA), from sources such as marine biological activity, biomass burning, and industrial emissions . Typical atmospheric gas-phase concentrations of MSA are on the order of ∼105–107 molec. cm−3, corresponding to approximately 10 %–100 % of SA concentrations .
Both experimental and theoretical studies have demonstrated that MSA can participate in cluster formation with bases such as AM and amines , particularly under conditions where SA concentrations are limited. Recent computational studies have further highlighted the role of MSA in stabilizing acid–base clusters and influencing their growth pathways, suggesting that its contribution to NPF may be more significant than previously assumed .
Despite significant recent progress, several important questions regarding mixed SA–MSA systems remain unresolved. Previous computational work has demonstrated that incorporating MSA into SA clusters can substantially enhance cluster formation potential by up to approximately 20-fold and identified the SA–MSA–DMA system as a particularly promising candidate for extending investigations to larger cluster sizes . Experimental studies have likewise demonstrated co-condensation of SA, MSA, and atmospheric bases, providing direct evidence for the formation of mixed-acid clusters under laboratory conditions . More recently, reported that MSA can substantially enhance both nucleation and particle growth through synergistic interactions with sulfuric acid under cold marine conditions. Consistently, showed through quantum-chemical calculations, cluster-dynamics simulations, and global modeling that MSA can enhance SA–NH3 nucleation rates by 1–3 orders of magnitude and may dominate marine upper-tropospheric new particle formation. However, the molecular-level thermodynamics, preferred compositions, and structural evolution of large mixed-acid clusters beyond the smallest cluster sizes remain largely unexplored. In particular, it is still unclear how the incorporation of MSA influences the stability and composition of particle-like clusters during the transition from initial molecular clusters toward freshly nucleated particles. The present work addresses this gap by extending quantum chemical calculations to mixed-acid clusters containing up to 10 acid–base pairs (20 molecules).
In this work, we address this challenge by extending computational investigations to chemically complex mixed-acid–base clusters. Specifically, we perform quantum chemical calculations on clusters of the form (MSA)n(base)n and (SA)n−1(MSA)1(base)n considering systems up to n=10 acid-base pairs (20 molecules).
MethodsComputational details
Density functional theory calculations used during the configurational sampling workflow – including single-point energies, geometry optimizations, and vibrational frequency calculations – were carried out with the empirically corrected B97-3c method as implemented in the ORCA 5.0.4 quantum chemistry package . Semi-empirical calculations were performed with the xTB 6.4.0 program using both the original GFN1-xTB method and a reparameterized variant, GFN1-xTBre-par, developed by on the basis of previously reported FNP structures . Initial cluster geometries were generated with ABCluster version 3.2 employing the CHARMM force field . Additional conformational exploration was performed using CREST in non-covalent interaction mode through metadynamics simulations .
Configurational sampling of (MSA)<sub>1</sub>(SA)<sub><italic>n</italic>−1</sub>(AM/MA/DMA)<inline-formula><mml:math id="M19" display="inline"><mml:mi>n</mml:mi></mml:math></inline-formula>, with <inline-formula><mml:math id="M20" display="inline"><mml:mrow><mml:mi>n</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:math></inline-formula>–5 clusters
We initially studied (MSA)1(SA)n−1(AM/MA/DMA)n clusters with n=1 to 5 acid molecules. The configurational sampling of the clusters was performed using our recently developed improved configuration sampling workflow . The configurational sampling procedure can be outlined as follows:
ABC⟶N=10000xTBOPT⟶N=10000B97-3cSP⟶filterN=10000B97-3cPART OPT⟶filterN=10000B97-3cFULL OPT.
Initial cluster geometries were generated using 10 parallel ABCluster runs, generating 10 000 local minimum configurations. Previous works have shown that performing several parallel ABCluster searches improves the likelihood of finding the global minimum of large molecular clusters compared with a single extended search. Ionic monomers were employed while maintaining overall charge neutrality of the clusters. This enforces proton transfer processes observed in the lowest free-energy structures. All generated structures were subsequently optimized using the semi-empirical GFN1-xTB method . DFT single-point energies were then evaluated at the B97-3c level of theory on the GFN1-xTB optimized conformers.
To reduce the number of total conformers, high-energy configurations were removed by retaining only the 1000 lowest-energy conformers for further refinement. These structures were subjected to partial geometry optimization with n×2 iterations (n being the number of acid–base pairs) at the B97-3c level to efficiently eliminate energetically unfavorable configurations. From this subset, the 100 most stable structures were selected for full geometry optimization followed by vibrational frequency calculations.
The lowest free-energy conformer was then selected for additional CREST exploration as suggested by , using the following workflow:
CREST⟶N=100B97-3cFULLOPT.
The configurational exploration was then carried out using CREST at the GFN1-xTB/GFN1-xTBre-par level in the non-covalent interaction mode. Then, the 100 lowest-energy geometries were selected for full optimization and vibrational frequency calculations.
Machine-learning-enhanced configurational sampling of (MSA)<sub>1</sub>(SA)<sub><italic>n</italic>−1</sub>(AM/MA/DMA)<inline-formula><mml:math id="M33" display="inline"><mml:mi>n</mml:mi></mml:math></inline-formula>, with <inline-formula><mml:math id="M34" display="inline"><mml:mrow><mml:mi>n</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="normal">6</mml:mn></mml:mrow></mml:math></inline-formula>–10 clusters
The PaiNN architecture, as implemented in SchNetPack , was used to learn the energy and forces at the B97-3c level of theory. We used the hyperparameters suggested by , who found the PaiNN model to perform well in reproducing the geometries at the B97-3c level of theory for atmospheric molecular clusters. The (MSA)1(SA)n−1(AM/MA/DMA)n cluster structures with n=1–5 identified in the previous section were used as a training/validation set. A separate model was trained for each base system. The training data consisted of out-of-equilibrium structures (the initial GFN1-xTB structures with B97-3c single-point energies and forces) and equilibrium structures calculated at the B97-3c level. The number of structures was filtered based on gyration radius, electronic energy, and dipole moment using arbitrary thresholds until a manageable amount of data were left. This resulted in a training/validation set of 9645, 6746, and 8681 structures for the SA–MSA–AM, SA–MSA–MA, and SA–MSA–DMA models, respectively. The models were trained for 600 epochs using the loss function described by , yielding the validation metrics reported in Table .
Mean absolute errors (MAEs) and root mean squared errors (RMSEs) for the validation set at the final epoch. Here, ΔE denotes the per-cluster binding electronic energy (kcal mol−1), and F denotes the force components (kcal mol−1 Å−1).
Using the trained neural networks, the configurational sampling of the n=6–10 clusters was performed as follows:
ABC⟶N=10000xTBOPT⟶N=10000NNOPT⟶N=10000B97-3cSP⟶filterN=100B97-3cFULL OPT,
where the optimization using the neural network was performed using ASE , with the BFGS optimizer, for 600 steps. Although we found the geometries to be reliable and closer to the true B97-3c minimum, a single-point calculation at the B97-3c was needed to get the correct energy order. The 100 lowest electronic energy geometries were selected for full optimization, and vibrational frequency calculations were performed. This was followed by the CREST track workflow.
Configurational sampling of (MSA)<sub><italic>n</italic></sub>(AM/MA/DMA)<sub><italic>n</italic></sub>, with <inline-formula><mml:math id="M52" display="inline"><mml:mrow><mml:mi>n</mml:mi><mml:mo>=</mml:mo><mml:mn mathvariant="normal">1</mml:mn></mml:mrow></mml:math></inline-formula>–10 clusters
The pure (MSA)n(AM/MA/DMA)n clusters were sampled without machine learning using the full workflow in Sect. for the set of clusters with n up to 10. The SA–base clusters were taken from . Nevertheless, all three cluster series underwent extensive configurational sampling protocol, and the final selected candidate structures were optimized consistently at the reference quantum–chemical level. Moreover, the machine-learning (ML)-assisted optimization is expected to improve the candidate geometries relative to those obtained from the initial xTB sampling, bringing them closer to low-energy regions of the reference potential-energy surface.
Cluster-binding free energies
The thermodynamic stability of the clusters was calculated through their standard binding free energies defined as the difference between the Gibbs free energy of the cluster and the sum of the Gibbs free energies of the isolated monomers. The binding free energy is therefore expressed as
ΔGbind=Gcluster-∑iGmonomer,i.
Equation () above describes the intrinsic thermochemistry of the clusters under standard conditions. The quasi-harmonic approximation, with a 100 cm−1 threshold, was employed to calculate the vibrational entropy contribution .
To evaluate effective binding free energies under specific atmospheric conditions, we employ the self-consistent distribution formalism ,
ΔGbind(p)=ΔGbind-RT⋅1-1n⋅∑ilnpipref.
Here, n denotes the total number of molecules in the cluster, pref denotes the reference pressure (1 atm), and pi represents the partial pressure of monomer i. This self-consistent formulation ensures that the reference free energies of the isolated monomers are properly defined as zero under the chosen thermodynamic reference state.
Binding Gibbs free energies (ΔGbind) at the B97-3c level of theory under standard conditions (STD) for clusters containing SA, MSA, and bases (AM, MA, DMA). The shaded region represents the range between minimum and maximum values for each cluster size, while annotations indicate the acid composition (SA, MSA).
Results and discussionBinding free energies under standard conditions
Using the configurational sampling workflow described above, we studied clusters containing SA, MSA, and atmospherically relevant bases, focusing on the (acid)n(base)n stoichiometries up to n=10. We studied SA–MSA–AM, SA–MSA–MA, and SA–MSA–DMA clusters, as well as the corresponding limiting case containing only MSA as the acidic component. The mixed-acid clusters investigated in this work represent a systematic substitution pathway designed to isolate the effect of incorporating a single MSA molecule into growing sulfuric acid clusters. In our previous work , we found that clusters containing a single MSA molecule represented the most likely composition under most conditions. Consequently, the present calculations provide molecular-level insight into the thermodynamic role of MSA incorporation rather than predicting the most abundant mixed-acid cluster compositions under atmospheric conditions. The calculated standard binding free energies at 298.15 K and 1 atm are shown as a function of the number of molecules in the cluster m, and each point is labeled according to the acid composition, e.g., (SA, MSA).
Based on the standard binding free energies (Fig. ), we observe a clear hierarchy: clusters containing SA are consistently more stable than the corresponding clusters containing only MSA. This trend is observed for all three bases studied here. This is expected chemically since SA is a considerably stronger acid than MSA and therefore leads to stronger acid–base interactions. As a result, SA-containing clusters are thermodynamically stabilized against evaporation. This trend also persists as cluster size increases. Replacing one SA molecule with MSA does not provide any substantial energetic benefit and, in many cases, leads to less favorable binding free energies. This behavior is particularly pronounced for the DMA-containing systems, where the strong SA–DMA interaction already defines the thermodynamic landscape, leaving little room for MSA to contribute meaningfully to stabilization. Hence, for DMA-containing systems, MSA does not appear to participate in any notable acid synergy at the earliest stages of cluster formation.
A particularly clear example is the DMA system at the largest size studied (20 monomers, i.e., 10 acids + 10 bases): the pure SA–DMA cluster has a binding free energy of roughly -264 kcal mol−1, while the pure MSA–DMA analogue is about -183 kcal mol−1, a difference of ∼ 81 kcal mol−1 in favor of SA. The same qualitative pattern holds for AM and MA (e.g., at 20 monomers, MSA-only clusters are ∼ 55–62 kcal mol−1 less stable than SA-only clusters for AM and MA, respectively).
We note that this mixed-acid behavior differs markedly from the mixed-base synergy previously observed in SA–AM–DMA systems . In those mixed-base systems, the introduction of one or two AM molecules into SA–DMA clusters produced a clear stabilization effect for larger clusters despite DMA being the intrinsically stronger base. This behavior was attributed to the higher hydrogen bond coordination capacity of AM, which enabled the formation of more favorable and increasingly particle-like bonding networks as the clusters grew.
In contrast, the present mixed-acid systems do not exhibit a comparable synergistic stabilization upon replacing SA with MSA. Although mixed-acid clusters remain less stable than their pure-SA counterparts, the energetic penalty associated with MSA substitution decreases systematically with increasing cluster size. This trend indicates that the stronger intrinsic acidity of SA dominates the thermodynamics of the smallest clusters, whereas incorporation of MSA becomes progressively more favorable during growth. The effect is particularly evident for the SA–MSA–AM and SA–MSA–MA systems, where the free-energy differences relative to the pure-SA clusters diminish at larger sizes, suggesting a greater role for MSA in later growth stages than in the initial nucleation process.
This suggests that, although MSA does not enhance stability relative to pure-SA clusters, its incorporation becomes increasingly favorable as cluster size increases. A plausible explanation is that larger clusters possess more extended hydrogen-bonding networks (Fig. ), allowing the methyl-group perturbation to be distributed across the cluster. In contrast, the more compact bonding motifs of smaller clusters are less able to accommodate this perturbation, resulting in a larger energetic penalty for MSA substitution.
Comparison of pure and mixed-acid clusters for the MA system. (a) (SA)8(MA)8, (b) (MSA)1(SA)7(MA)8, (c) (SA)9(MA)9, (d) (MSA)1(SA)8(MA)9.
Structurally, both the pure-SA and mixed-acid clusters exhibit a gradual transition from relatively compact molecular complexes toward increasingly particle-like three-dimensional hydrogen-bonding networks as cluster size increases. The incorporation of one MSA molecule does not qualitatively alter this cluster-to-particle transition behavior. Instead, the main effect of MSA appears to be a modest perturbation of the hydrogen-bonding arrangement due to the presence of the methyl group, which slightly reduces the hydrogen-bond coordination efficiency compared to SA. This effect may be particularly important in DMA-containing clusters, where the additional methyl groups of DMA already impose steric constraints on the cluster geometry. Consequently, simultaneous incorporation of both DMA and MSA may further reduce the ability of the cluster to form compact and efficiently coordinated hydrogen-bonding motifs.
Overall, the standard free-energy data indicate that MSA is not a strong nucleating acid on its own, and SA remains the dominant acid controlling the thermodynamics of initial cluster formation. However, MSA could potentially incorporate into the cluster at larger sizes, thus contributing to the cluster growth. In general, the contribution of MSA to cluster stabilization is highly dependent on both cluster size and base type.
Binding free energies under given conditions
Based on the calculated standard binding free energies described above, we recalculated the binding free energies under atmospherically relevant monomer concentrations and temperatures using Eq. (). Figures – present the binding Gibbs free energies of the mixed-acid clusters at 278.15 and 298.15 K under four representative concentration regimes. These include high-concentration conditions at 278.15 K ([SA]=[MSA]=108 molec. cm−3; denoted 278H), high-MSA conditions at 278.15 K ([SA]=106 and [MSA]=108 molec. cm−3; denoted 278Hmsa), low-concentration conditions at 298.15 K ([SA]=[MSA]=106 molec. cm−3; denoted 298L), and high-MSA conditions at 298.15 K ([SA]=106 and [MSA]=108 molec. cm−3; denoted 298Lmsa).
In all cases, the base concentrations were fixed at AM =10 ppb and amines (MA, DMA) =5 ppt. These conditions were chosen to mimic typical chamber and atmospheric environments and to assess how the relative abundance of SA and MSA influences the stability and composition of the earliest clusters. In the following, we discuss each scenario separately.
Low SA–low MSA concentrations at 298.15 K
Under low-concentration conditions at 298.15 K, the calculated free energies (Fig. ) show that clusters containing SA remain clearly favored over those dominated by MSA. This trend is particularly strong for the DMA-containing system, where the SA–MSA–DMA clusters essentially show no energetic advantage from incorporating MSA. The lowest-free-energy clusters remain those with the highest SA content, demonstrating that MSA does not compete effectively with SA when paired with a strong base such as DMA.
Binding Gibbs free energies (ΔGbind) at the B97-3c level of theory at 298.15 K under low-concentration conditions (298L), with [SA]=[MSA]=106 molec. cm−3. Base concentrations are fixed at AM =10 ppb and amines (MA, DMA) =5 ppt.
The AM- and MA-containing systems are somewhat more sensitive to acid composition; however, even in these cases, the pure or SA-rich clusters remain lowest in free energy at the smallest sizes. Thus, under warm and dilute conditions, MSA does not appear to be capable of driving the initial formation of stable clusters. Instead, it behaves as a weaker auxiliary acid that may only enter the cluster after SA has already established the primary acid–base scaffold. Hence, under such conditions, it is expected that MSA only affects the growth of FNPs. Overall, under ambient conditions and at low concentrations of the acids, the thermodynamic landscape is uphill, and, thereby, it is unlikely that FNP formation will occur.
This result has important implications for previous experimental observations reporting efficient particle formation in MSA–MA systems. The present calculations suggest that binary MSA–base interactions alone are unlikely to fully explain the experimentally observed particle formation efficiencies. In particular, our simulations do not include additional stabilizing species such as water molecules or ions, which are commonly present under experimental and atmospheric conditions and are known to substantially enhance cluster stability and suppress evaporation. Indeed, previous experimental studies have also highlighted the important role of water in promoting particle formation in MSA-containing systems . Therefore, the discrepancy between our calculations and experimentally observed particle formation rates is likely to be related to the simplified dry-cluster framework employed here rather than an inconsistency with the experimental observations themselves.
Low-SA–high-MSA concentrations at 298.15 K
At 298.15 K and higher MSA concentrations, the overall cluster free energies (Fig. ) become less positive. Under these conditions, the differences between cluster compositions narrow somewhat, and mixed-acid clusters become more competitive in the AM- and MA-containing systems.
Binding Gibbs free energies (ΔGbind) at the B97-3c level of theory at 298.15 K under low-SA- and high-MSA-concentration conditions (298Lmsa), with [SA]=106 and [MSA]=108 molec. cm−3. Base concentrations are fixed at AM =10 ppb and amines (MA, DMA) =5 ppt.
In particular, the SA–MSA–AM and SA–MSA–MA clusters show that MSA can begin to contribute to cluster stabilization when its relative abundance is sufficiently high. For instance, for SA–MSA–AM we see clusters containing one MSA being most stable across almost all size ranges. This does not change the overall conclusion that SA remains the dominant nucleating acid, but it suggests that MSA can influence the growth trajectory once the cluster has reached a sufficient size. Although this effect remains clearly weaker than that of SA, it indicates that MSA may become incorporated into growing clusters under favorable concentration conditions.
For the SA–MSA–DMA system, however, the conclusions remain unchanged. Even at elevated MSA concentrations, the strong SA–DMA interaction dominates the cluster energetics, and MSA contributes little to the stabilization. Thus, the presence of a strong base suppresses any meaningful acid synergy between SA and MSA. Overall, we see the role of MSA in cluster formation at 298.15 K despite having a 100-times-higher concentration than SA. However, the overall conclusion is similar to the low-SA–low-MSA concentrations at 298.15 K, i.e., it is unlikely that FNP formation will occur at these conditions.
High-MSA–low-SA concentrations at 278.15 K
At 278.15 K and high concentrations of MSA, the free-energy surfaces (Fig. ) become increasingly downhill for the larger clusters containing AM and MA, indicating more favorable growth conditions under these low temperatures and high concentrations of MSA. Nevertheless, both systems still exhibit an initial nucleation barrier at the smallest cluster sizes before growth becomes thermodynamically favorable. In contrast, the DMA-containing clusters remain strongly stabilized, primarily by SA, while the pure MSA–DMA clusters remain purely uphill despite the enhanced MSA concentration.
Binding Gibbs free energies (ΔGbind) at the B97-3c level of theory at 278.15 K under high-MSA-concentration conditions (278Hmsa), with [SA]=106 and [MSA]=108 molec. cm−3. Base concentrations are fixed at AM =10 ppb and amines (MA, DMA) =5 ppt.
Under these conditions, the AM- and MA-containing mixed-acid systems again show the greatest sensitivity to MSA incorporation. In these systems, MSA can contribute to cluster stabilization at intermediate and larger sizes, especially when SA is not overwhelmingly dominant in the gas phase. However, even under the most favorable conditions for MSA participation, the calculations do not support a scenario in which MSA replaces SA as the primary nucleating acid. Instead, MSA acts as a secondary growth-modifying species, entering the cluster once the main SA-driven acid–base clusters are established. In the SA–MSA–DMA system, the strong SA–DMA pair continues to dominate the thermodynamics, and MSA again contributes negligibly. Thus, even when both low temperatures and high concentrations favor clustering, MSA cannot compete with sulfuric acid in the presence of the strongest base.
High-MSA–high-SA concentrations at 278.15 K
At the lower temperature of 278.15 K and high concentrations (Fig. ), all clusters become more stable compared to the corresponding 298.15 K case. Nevertheless, SA-containing clusters are still significantly favored over MSA-only clusters, confirming that the intrinsic weakness of MSA relative to SA is not compensated for simply by lowering the temperature.
Binding Gibbs free energies (ΔGbind) at the B97-3c level of theory at 278.15 K under high-concentration conditions (278H), with [SA]=[MSA]=108 molec. cm−3. Base concentrations are fixed at AM =10 ppb and amines (MA, DMA) =5 ppt.
The SA–MSA–AM and SA–MSA–MA systems again exhibit compositions in which MSA becomes more competitive at larger cluster sizes, indicating a possible role during the early growth stages. Under these conditions, the free-energy barriers for the AM- and MA-containing systems also go downhill after the initial dimer cluster is formed compared to the warmer conditions, suggesting that particle formation may become thermodynamically more favorable. Nevertheless, the smallest thermodynamically viable clusters remain predominantly SA-dominated.
For the SA–MSA–DMA systems, the strong SA–DMA interaction continues to dominate the thermodynamics, and the corresponding free-energy surfaces become nearly fully downhill under these conditions. In contrast, the pure MSA-containing DMA clusters remain significantly less favorable, indicating that MSA still contributes only weakly to the stabilization pattern in the presence of DMA.
Overall, these results suggest that, although lowering the temperature enhances clustering and reduces nucleation barriers in general, the relative importance of MSA is governed more strongly by composition and concentration than by temperature alone. Thus, even under favorable low-temperature and high-MSA conditions, SA remains the key acid controlling the earliest nucleation steps.
Low-temperature and MSA-rich conditions at 258.15 K
Finally, we tested an additional condition at a relatively low temperature of 258.15 K with an SA : MSA ratio of 1:1000. However, the highest MSA concentration considered in this work (109 molec. cm−3) was intentionally chosen as a sensitivity test rather than to represent a typical atmospheric condition. This allows us to identify the concentration regime at which MSA begins to noticeably influence cluster stability and to better understand its potential thermodynamic role. Under more representative atmospheric MSA concentrations, the effect of MSA is expected to be correspondingly smaller. Under low-temperature conditions at 258.15 K with a strong excess of MSA relative to SA (Fig. ), the stability of MSA-containing clusters increases substantially compared to warmer conditions. In particular, the pure MSA-containing systems become significantly more competitive at larger cluster sizes, indicating that low temperature combined with high MSA abundance can enhance the thermodynamic contribution of MSA to cluster growth.
Binding Gibbs free energies (ΔGbind) at the B97-3c level of theory at 258.15 K under MSA-rich conditions (258Lmsa), with [SA]=106 and [MSA]=109 molec. cm−3. Base concentrations are fixed at AM =10 ppb and amines (MA, DMA) =5 ppt.
The effect is most pronounced for the DMA-containing clusters, where the free energies of the pure MSA systems decrease considerably relative to the corresponding higher-temperature cases. Nevertheless, even under these highly MSA-favored conditions, mixed SA-containing clusters generally remain among the most stable configurations, particularly at smaller cluster sizes. The SA–MSA–AM and SA–MSA–MA systems also exhibit enhanced stabilization from MSA incorporation at larger sizes, suggesting that MSA may become increasingly relevant during the cluster growth regime.
Overall, while SA-containing clusters still dominate the earliest nucleation steps, these results indicate that the thermodynamic importance of MSA increases strongly under cold and MSA-rich atmospheric conditions. This is consistent with the previous work by . The formation of MSA from the oxidation of DMS by OH radicals has been shown to be temperature dependent . Hence, under colder conditions, the MSA:SA ratio will be high due to a shift from the hydrogen abstraction pathway to the OH radical addition pathway. This indicates that MSA might contribute to the early growth of FNPs but is still unlikely to participate in the very initial particle formation steps when SA and DMA are present.
Conclusions
Here, we find that the acid synergy between SA and MSA is qualitatively different and significantly weaker than the mixed-base synergy previously observed for SA–AM–DMA clusters. In the mixed-base case, combining strong and weak bases showed clear synergistic stabilization effects in larger clusters, whereas the present mixed-acid systems remain dominated by the intrinsically stronger acidity of SA relative to MSA.
Overall, the calculations show that, within the neutral AM-, MA-, and DMA-containing systems and cluster compositions investigated here, pure MSA-containing acid–base clusters are generally thermodynamically less stable than the corresponding SA-containing clusters under most of the considered conditions. Instead, SA-containing clusters consistently dominate the earliest cluster formation steps, particularly in systems involving DMA, where the strong SA–DMA interaction leads to nearly barrierless cluster growth under favorable conditions. In contrast, the AM- and MA-containing systems exhibit larger nucleation barriers, although these barriers decrease substantially at lower temperatures and under elevated MSA concentrations.
Within the specific neutral amine-containing systems and single-MSA compositions investigated here, the calculations indicate that MSA incorporation becomes thermodynamically more favorable as the cluster size increases. This trend should not, however, be interpreted as evidence that MSA contributes exclusively or predominantly to growth rather than initial nucleation under all atmospheric conditions. An important limitation of the present calculations is the omission of water. Recent CLOUD experiments demonstrated that, between +10 and -10 °C, MSA drives rapid particle growth at high relative humidity (>50 %), whereas its contribution is negligible under low-humidity conditions (<16 %) . The comparatively weak contribution of MSA predicted by our unhydrated calculations is therefore not necessarily inconsistent with these experimental observations. Rather, the comparison highlights the important role of water in stabilizing MSA-containing clusters and indicates that hydration should be included in future thermodynamic and kinetic investigations. In particular, MSA incorporation becomes increasingly favorable at larger cluster sizes, especially in SA–MSA–AM and SA–MSA–MA systems under conditions of low temperature and elevated MSA:SA ratios. Under the most MSA-rich and coldest conditions considered here (258.15 K and SA:MSA=1:1000), the thermodynamic stability of MSA-containing clusters increases substantially, and the free-energy penalties associated with MSA incorporation decrease remarkably. Nevertheless, even under these strongly MSA-favored conditions, SA-containing clusters remain among the most stable structures at the smallest cluster sizes.
While these thermodynamic properties provide important insight into the potential role of MSA during cluster formation and early growth, quantitative predictions of atmospheric formation rates and particle growth kinetics require comprehensive kinetic modeling, which is beyond the scope of the present work.
From an atmospheric perspective, these findings suggest that the role of MSA in particle formation depends strongly on environmental conditions. Over marine environments and colder atmospheric regions, where dimethyl sulfide oxidation may favor MSA formation relative to SA, MSA could become increasingly important for stabilizing growing clusters and enhancing early particle growth. For the specific neutral amine-containing systems investigated here, the thermodynamic results indicate that lower temperatures and high MSA:SA ratios increase the favorability of MSA incorporation. However, the present restricted thermodynamic analysis cannot establish the atmospheric conditions under which MSA-containing pathways become kinetically competitive with SA-driven nucleation.
Finally, it is important to emphasize that this study does not include additional stabilizing species such as water molecules or ions. Previous experimental studies have highlighted the important role of water in promoting particle formation in MSA-containing systems, and inclusion of water may substantially reduce evaporation rates and enhance cluster stability. Thus, while the present work isolates the intrinsic thermodynamic role of MSA in mixed-acid clustering, future studies including hydration and ionic effects will be necessary for a more complete description of atmospheric particle formation under realistic marine conditions. Such studies should also consider NH3, clusters containing multiple MSA molecules, and cluster-dynamics simulations.
Data availability
All of the calculated structures and all of the thermochemistry are available in the Atmospheric Cluster Database (ACDB) at 10.1021/acsomega.9b00860.
Author contributions
Conceptualization: J. E.; methodology: G. H., Y. K., H. W., J. E.; formal analysis: G. H., Y. K., H. W.; investigation: G. H., Y. K., H. W.; resources: J. E.; writing (original draft): G. H.; writing (review and editing): G. H., Y. K., H. W., J. E.; visualization: G. H., H. W.; project administration: J. E.; funding acquisition: J. E.; supervision: J. E.
Competing interests
At least one of the (co-)authors is a member of the editorial board of Aerosol Research. The peer-review process was guided by an independent editor. The remaining authors have no conflict of interests to declare.
Disclaimer
This study was funded by the European Union. The views and opinions expressed are, however, those of the authors only and do not necessarily reflect those of the European Union, the European Research Executive Agency, or the European Research Council Executive Agency. Neither the European Union nor the granting authority can be held responsible for them.Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. The authors bear the ultimate responsibility for providing appropriate place names. Views expressed in the text are those of the authors and do not necessarily reflect the views of the publisher.
Acknowledgements
The numerical results presented in this work were obtained at the Centre for Scientific Computing, Aarhus – https://phys.au.dk/forskning/faciliteter/cscaa/ (last access: 13 October 2026).
Financial support
This work was funded by the European Union (ERC, ExploreFNP, project no. 101040353) and the Danish National Research Foundation (grant no. DNRF172) through the Center of Excellence for Chemistry of Clouds.
Review statement
This paper was edited by Bingbing Wang and reviewed by two anonymous referees.
ReferencesAlfaouri, D., Passananti, M., Zanca, T., Ahonen, L., Kangasluoma, J., Kubečka, J., Myllys, N., and Vehkamäki, H.: A study on the fragmentation of sulfuric acid and dimethylamine clusters inside an atmospheric pressure interface time-of-flight mass spectrometer, Atmos. Meas. Tech., 15, 11–19, 10.5194/amt-15-11-2022, 2022. Almeida, J., Schobesberger, S., Kürten, A., Ortega, I. K., Kupiainen-Määttä, O., Praplan, A. P., Adamov, A., Amorim, A., Bianchi, F., Breitenlechner, M., David, A., Dommen, J., Donahue, N. M., Downard, A., Dunne, E., Duplissy, J., Ehrhart, S., Flagan, R. C., Franchin, A., Guida, R., Hakala, J., Hansel, A., Heinritzi, M., Henschel, H., Jokinen, T., Junninen, H., Kajos, M., Kangasluoma, J., Keskinen, H., Kupc, A., Kurtén, T., Kvashin, A. N., Laaksonen, A., Lehtipalo, K., Leiminger, M., Leppä, J., Loukonen, V., Makhmutov, V., Mathot, S., McGrath, M. J., Nieminen, T., Olenius, T., Onnela, A., Petäjä, T., Riccobono, F., Riipinen, I., Rissanen, M., Rondo, L., Ruuskanen, T., Santos, F. D., Sarnela, N., Schallhart, S., Schnitzhofer, R., Seinfeld, J. H., Simon, M., Sipilä, M., Stozhkov, Y., Stratmann, F., Tomé, A., Tröstl, J., Tsagkogeorgas, G., Vaattovaara, P., Viisanen, Y., Virtanen, A., Vrtala, A., Wagner, P. E., Weingartner, E., Wex, H., Williamson, C., Wimmer, D., Ye, P., Yli-Juuti, T., Carslaw, K. S., Kulmala, M., Curtius, J., Baltensperger, U., Worsnop, D. R., Vehkamäki, H., and Kirkby, J.: Molecular Understanding of Sulphuric Acid–Amine Particle Nucleation in the Atmosphere, Nature, 502, 359–363, 2013.Baalbaki, R., Shen, J., Simon, M., Klebach, H., Ruhl, S., DeVivo, J., Wang, M., Scholz, W., Dada, L., Rörup, B., Stolzenburg, D., Manninen, H. E., Sommer, E., Caudillo-Plath, L., Marie, G., Friedrich, M., Yu, W., Leiminger, M., Alfaouri, D., Amorim, A., Arnoldi-Meadows, T., Beckmann, H., Berntheusel, M., Bräkling, S., Brasseur, Z., Chiu, R., Duplissy, J., Finkenzeller, H., Heinritzi, M., Kunkler, F., Lamkaddam, H., Lopez, B., Mahfouz, N., Makhmutov, V., Martinez, M., Marten, R., Massabo, D., Mauldin, R., Mentler, B., Müller, M., Philippov, M., Piedehierro, A. A., Rato, P., Reinecke, T., Richter, S., Russell, D. M., Schulze, B., Surdu, M., Thakur, R. C., Tham, Y. J., Tian, P., Tomé, A., Tong, Y., Top, J., Wagner, A. C., Wang, D. S., Wang, Y., Ward, R. X., Weber, S. K., Welti, A., Wu, Y., Zauner-Wieczorek, M., Zhang, J., Curtius, J., Donahue, N. M., El Haddad, I., Flagan, R. C., Hansel, A., Harder, H., Kürten, A., Petäjä, T., Schobesberger, S., Sipilä, M., Volkamer, R., Winkler, P. M., Worsnop, D. R., Christoudias, T., Pozzer, A., Kulmala, M., Kirkby, J., Lehtipalo, K., and He, X.-C.: Role of methanesulfonic acid in atmospheric particle nucleation and growth, Nature, pp. 1–4, 10.1038/s41586-026-10810-2, 2026.Bannwarth, C., Caldeweyher, E., Ehlert, S., Hansen, A., Pracht, P., Seibert, J., Spicher, S., and Grimme, S.: Extended Tight-Binding Quantum Chemistry Methods, WIREs Comput. Mol. Sci., 11, e1493, 10.1002/wcms.1493, 2021. Barnes, I., Hjorth, J., and Mihalopoulos, N.: Dimethyl Sulfide and Dimethyl Sulfoxide and Their Oxidation in the Atmosphere, Chem. Rev., 106, 940–975, 2006. Bates, T., Lamb, B., Guenther, A., Dignon, J., and Stoiber, R.: Sulfur Emissions to the Atmosphere from Natural Sources, J. Atmos. Chem., 14, 315–337, 1992. Bianchi, F., Tröstl, J., Junninen, H., Frege, C., Henne, S., Hoyle, C. R., Molteni, U., Herrmann, E., Adamov, A., Bukowiecki, N., Chen, X., Duplissy, J., Gysel, M., Hutterli, M., Kangasluoma, J., Kontkanen, J., Kürten, A., Manninen, H. E., Münch, S., Peräkylä, O., Petäjä, T., Rondo, L., Williamson, C., Weingartner, E., Curtius, J., Worsnop, D. R., Kulmala, M., Dommen, J., and Baltensperger, U.: New Particle Formation in the Free Troposphere: A Question of Chemistry and Timing, Science, 352, 1109–1112, 2016.Brandenburg, J. G., Bannwarth, C., Hansen, A., and Grimme, S.: B97-3c: A Revised Low-Cost Variant of the B97-D Density Functional Method, J. Chem. Phys., 148, 064104, 10.1063/1.5012601, 2018. Chen, H., Ezell, M. J., Arquero, K. D., Varner, M. E., Dawson, M. L., Gerber, R. B., and Finlayson-Pitts, B. J.: New Particle Formation and Growth from Methanesulfonic Acid, Trimethylamine and Water, Phys. Chem. Chem. Phys., 17, 13699–13709, 2015.Cooley, S. R., Schoeman, D. S., Bopp, L., Boyd, P., Donner, S., Ghebrehiwet, D. Y., Ito, S.-I., Kiessling, W., Martinetto, P., Ojea, E., Racault, M.-F., Rost, B., and Skern-Mauritzen, M.: Oceans and Coastal Ecosystems and Their Services, in: Climate Change 2022: Impacts, Adaptation and Vulnerability, pp. 379–550, Cambridge University Press, Cambridge, UK and New York, NY, USA, 10.1017/9781009325844.005, 2023.Cromar, K. and Lazrak, N.: Risk Communication of Ambient Air Pollution in the WHO European Region: Review of Air Quality Indexes and Lessons Learned, license: CC BY-NC-SA 3.0 IGO, https://iris.who.int/handle/10665/365787 (last access: 14 October 2026), 2023. DePalma, J. W., Bzdek, B. R., Doren, D. J., and Johnston, M. V.: Structure and Energetics of Nanometer-Size Clusters of Sulfuric Acid with Ammonia and Dimethylamine, J. Phys. Chem. A, 116, 1030–1040, 2012. DePalma, J. W., Doren, D. J., and Johnston, M. V.: Formation and Growth of Molecular Clusters Containing Sulfuric Acid, Water, Ammonia, and Dimethylamine, J. Phys. Chem. A, 118, 5464–5473, 2014.Eisele, F. and Tanner, D.: Measurement of the Gas-Phase Concentration of H2SO4 and Methane Sulfonic Acid and Estimates of H2SO4 Production and Loss in the Atmosphere, J. Geophys. Res.-Atmos., 98, 9001–9010, 1993.Elm, J.: An Atmospheric Cluster Database Consisting of Sulfuric Acid, Bases, Organics, and Water, ACS Omega, 4, 10965–10974, 10.1021/acsomega.9b00860, 2019. Elm, J.: Clusteromics III: Acid Synergy in Sulfuric Acid–Methanesulfonic Acid–Base Cluster Formation, ACS Omega, 7, 15206–15214, 2022. Facchini, M. C., Decesari, S., Rinaldi, M., Carbone, C., Finessi, E., Mircea, M., Fuzzi, S., Moretti, F., Tagliavini, E., Ceburnis, D., et al. : Important Source of Marine Secondary Organic Aerosol from Biogenic Amines, Environ. Sci. Technol., 42, 9116–9121, 2008. Ge, X., Wexler, A. S., and Clegg, S. L.: Atmospheric Amines–Part I. A Review, Atmos. Environ., 45, 524–546, 2011. Grimme, S.: Supramolecular Binding Thermodynamics by Dispersion-Corrected Density Functional Theory, Chem. Eur. J, 18, 9955–9964, 2012. Grimme, S.: Exploration of Chemical Compound, Conformer, and Reaction Space with Meta-Dynamics Simulations Based on Tight-Binding Quantum Chemical Calculations, J. Chem. Theory Comput., 15, 2847–2862, 2019.Grimme, S., Bannwarth, C., and Shushkov, P.: A Robust and Accurate Tight-Binding Quantum Chemical Method for Structures, Vibrational Frequencies, and Noncovalent Interactions of Large Molecular Systems Parametrized for All spd-Block Elements (Z=1–86), J. Chem. Theory Comput., 13, 1989–2009, 2017.Halonen, R.: A Consistent Formation Free Energy Definition for Multicomponent Clusters in Quantum Thermochemistry, J. Aerosol Sci., 162, 105974, 10.1016/j.jaerosci.2022.105974, 2022.Hasan, G., Wu, H., Knattrup, Y., and Elm, J.: Base synergy in freshly nucleated particles, Aerosol Research, 3, 101–111, 10.5194/ar-3-101-2025, 2025. Henschel, H., Navarro, J. C. A., Yli-Juuti, T., Kupiainen-Määttä, O., Olenius, T., Ortega, I. K., Clegg, S. L., Kurtén, T., Riipinen, I., and Vehkamäki, H.: Hydration of Atmospherically Relevant Molecular Clusters: Computational Chemistry and Classical Thermodynamics, J. Phys. Chem. A, 118, 2599–2611, 2014. Henschel, H., Kurtén, T., and Vehkamäki, H.: Computational Study on the Effect of Hydration on New Particle Formation in the Sulfuric Acid/Ammonia and Sulfuric Acid/Dimethylamine Systems, J. Phys. Chem. A, 120, 1886–1896, 2016. Herb, J., Nadykto, A. B., and Yu, F.: Large Ternary Hydrogen-Bonded Pre-Nucleation Clusters in the Earth's Atmosphere, Chem. Phys. Lett., 518, 7–14, 2011. Huang, J. and MacKerell Jr, A. D.: CHARMM36 All-Atom Additive Protein Force Field: Validation Based on Comparison to NMR Data, J. Comput. Chem., 34, 2135–2145, 2013.Ianni, J. C. and Bandy, A. R.: A Density Functional Theory Study of the Hydrates of NH3⋅H2SO4 and Its Implications for the Formation of New Atmospheric Particles, J. Phys. Chem. A, 103, 2801–2811, 1999.Johnson, J. S. and Jen, C. N.: Role of methanesulfonic acid in sulfuric acid–amine and ammonia new particle formation, ACS Earth & Space Chemistry, 7, 653, 10.1021/acsearthspacechem.3c00017, 2023.Jokinen, T., Sipilä, M., Junninen, H., Ehn, M., Lönn, G., Hakala, J., Petäjä, T., Mauldin III, R. L., Kulmala, M., and Worsnop, D. R.: Atmospheric sulphuric acid and neutral cluster measurements using CI-APi-TOF, Atmos. Chem. Phys., 12, 4117–4125, 10.5194/acp-12-4117-2012, 2012. Kirkby, J., Curtius, J., Almeida, J., Dunne, E., Duplissy, J., Ehrhart, S., Franchin, A., Gagné, S., Ickes, L., Kürten, A., Kupc, A., Metzger, A., Riccobono, F., Rondo, L., Schobesberger, S., Tsagkogeorgas, G., Wimmer, D., Amorim, A., Bianchi, F., Breitenlechner, M., David, A., Dommen, J., Downard, A., Ehn, M., Flagan, R. C., Haider, S., Hansel, A., Hauser, D., Jud, W., Junninen, H., Kreissl, F., Kvashin, A., Laaksonen, A., Lehtipalo, K., Lima, J., Lovejoy, E. R., Makhmutov, V., Mathot, S., Mikkilä, J., Minginette, P., Mogo, S., Nieminen, T., Onnela, A., Pereira, P., Petäjä, T., Schnitzhofer, R., Seinfeld, J. H., Sipilä, M., Stozhkov, Y., Stratmann, F., Tomé, A., Vanhanen, J., Viisanen, Y., Vrtala, A., Wagner, P. E., Walther, H., Weingartner, E., Wex, H., Winkler, P. M., Carslaw, K. S., Worsnop, D. R., Baltensperger, U., and Kulmala, M.: Role of Sulphuric Acid, Ammonia and Galactic Cosmic Rays in Atmospheric Aerosol Nucleation, Nature, 476, 429–433, 2011. Kirkby, J., Duplissy, J., Sengupta, K., Frege, C., Gordon, H., Williamson, C., Heinritzi, M., Simon, M., Yan, C., Almeida, J., Tröstl, J., Nieminen, T., Ortega, I. K., Wagner, R., Adamov, A., Amorim, A., Bernhammer, A.-K., Bianchi, F., Breitenlechner, M., Brilke, S., Chen, X., Craven, J., Dias, A., Ehrhart, S., Flagan, R. C., Franchin, A., Fuchs, C., Guida, R., Hakala, J., Hoyle, C. R., Jokinen, T., Junninen, H., Kangasluoma, J., Kim, J., Krapf, M., Kürten, A., Laaksonen, A., Lehtipalo, K., Makhmutov, V., Mathot, S., Molteni, U., Onnela, A., Peräkylä, O., Piel, F., Petäjä, T., Praplan, A. P., Pringle, K., Rap, A., Richards, N. A. D., Riipinen, I., Rissanen, M. P., Rondo, L., Sarnela, N., Schobesberger, S., Scott, C. E., Seinfeld, J. H., Sipilä, M., Steiner, G., Stozhkov, Y., Stratmann, F., Tomé, A., Virtanen, A., Vogel, A. L., Wagner, A. C., Wagner, P. E., Weingartner, E., Wimmer, D., Winkler, P. M., Ye, P., Zhang, X., Hansel, A., Dommen, J., Donahue, N. M., Worsnop, D. R., Baltensperger, U., Kulmala, M., Carslaw, K. S., and Curtius, J.: Ion-Induced Nucleation of Pure Biogenic Particles, Nature, 533, 521–526, 2016. Kirkby, J., Amorim, A., Baltensperger, U., Carslaw, K. S., Christoudias, T., Curtius, J., Donahue, N. M., El Haddad, I., Flagan, R. C., Gordon, H., Hansel, A., Harder, H., Junninen, H., Kulmala, M., Kürten, A., Laaksonen, A., Lehtipalo, K., Lelieveld, J., Möhler, O., Riipinen, I., Stratmann, F., Tomé, A., Virtanen, A., Volkamer, R., Winkler, P. M., and Worsnop, D. R.: Atmospheric New Particle Formation from the CERN CLOUD Experiment, Nat. Geosci., 16, 948–957, 2023.Knattrup, Y., Kubečka, J., Wu, H., Jensen, F., and Elm, J.: Reparameterization of GFN1-xTB for Atmospheric Molecular Clusters: Applications to Multi-Acid–Multi-Base Systems, RSC Adv., 14, 20048–20055, 10.1039/D4RA03021D, 2024.Kubečka, J., Ayoubi, D., Tang, Z., Knattrup, Y., Engsvang, M., Wu, H., and Elm, J.: Accurate Modeling of the Potential Energy Surface of Atmospheric Molecular Clusters Boosted by Neural Networks, Environ. Sci. Adv., 3, 1438–1451, 10.1039/D4VA00255E, 2024. Kulmala, M., Kontkanen, J., Junninen, H., Lehtipalo, K., Manninen, H. E., Nieminen, T., Petäjä, T., Sipilä, M., Schobesberger, S., Rantala, P., Franchin, A., Jokinen, T., Järvinen, E., Äijälä, M., Kangasluoma, J., Hakala, J., Aalto, P. P., Paasonen, P., Mikkilä, J., Vanhanen, J., Aalto, J., Hakola, H., Makkonen, U., Ruuskanen, T., Mauldin, R. L., III, Duplissy, J., Vehkamäki, H., Bäck, J., Kortelainen, A., Riipinen, I., Kurtén, T., Johnston, M. V., Smith, J. N., Ehn, M., Mentel, T. F., Lehtinen, K. E. J., Laaksonen, A., Kerminen, V.-M., and Worsnop, D. R.: Direct Observations of Atmospheric Aerosol Nucleation, Science, 339, 943–946, 2013.Kupiainen, O., Ortega, I. K., Kurtén, T., and Vehkamäki, H.: Amine substitution into sulfuric acid – ammonia clusters, Atmos. Chem. Phys., 12, 3591–3599, 10.5194/acp-12-3591-2012, 2012.Kurtén, T., Torpo, L., Sundberg, M. R., Kerminen, V.-M., Vehkamäki, H., and Kulmala, M.: Estimating the NH3:H2SO4 ratio of nucleating clusters in atmospheric conditions using quantum chemical methods, Atmos. Chem. Phys., 7, 2765–2773, 10.5194/acp-7-2765-2007, 2007.Kurtén, T., Loukonen, V., Vehkamäki, H., and Kulmala, M.: Amines are likely to enhance neutral and ion-induced sulfuric acid-water nucleation in the atmosphere more effectively than ammonia, Atmos. Chem. Phys., 8, 4095–4103, 10.5194/acp-8-4095-2008, 2008.Larsen, A. H., Mortensen, J. J., Blomqvist, J., Castelli, I. E., Christensen, R., Dułak, M., Friis, J., Groves, M. N., Hammer, B., Hargus, C., Hermes, E. D., Jennings, P. C., Jensen, P. B., Kermode, J., Kitchin, J. R., Kolsbjerg, E. L., Kubal, J., Kaasbjerg, K., Lysgaard, S., Maronsson, J. B., Maxson, T., Olsen, T., Pastewka, L., Peterson, A., Rostgaard, C., Schiøtz, J., Schütt, O., Strange, M., Thygesen, K. S., Vegge, T., Vilhelmsen, L., Walter, M., Zeng, Z., and Jacobsen, K. W.: The atomic simulation environment—a Python library for working with atoms, J. Phys.: Condensed Matter, 29, 273002, 10.1088/1361-648X/aa680e, 2017. Larson, L. J., Largent, A., and Tao, M.: Structure of the Sulfuric Acid–Ammonia System and the Effect of Water Molecules in the Gas Phase, J. Phys. Chem. A, 103, 6786–6792, 1999. Liu, Y., Xie, H.-B., Ma, F., Chen, J., and Elm, J.: Amine-Enhanced Methanesulfonic Acid-Driven Nucleation: Predictive Model and Cluster Formation Mechanism, Environ. Sci. Technol., 56, 7751–7760, 2022. Loeb, N. G. and Kato, S.: Top-of-Atmosphere Direct Radiative Effect of Aerosols over the Tropical Oceans from the Clouds and the Earth's Radiant Energy System (CERES) Satellite Instrument, J. Climate, 15, 1474–1484, 2002.Loukonen, V., Kurtén, T., Ortega, I. K., Vehkamäki, H., Pádua, A. A. H., Sellegri, K., and Kulmala, M.: Enhancing effect of dimethylamine in sulfuric acid nucleation in the presence of water – a computational study, Atmos. Chem. Phys., 10, 4961–4974, 10.5194/acp-10-4961-2010, 2010.Ma, Y., Chen, J., Jiang, S., Liu, Y., Huang, T., Miao, S., Wang, C., and Huang, W.: Characterization of the Nucleation Precursor (H2SO4-(CH3)2NH) Complex: Intra-Cluster Interactions and Atmospheric Relevance, RSC Adv., 6, 5824–5836, 2016. McMurry, P. H.: The History of Condensation Nucleus Counters, Aerosol Sci. Technol., 33, 297–322, 2000.Meinardi, S., Simpson, I. J., Blake, N. J., Blake, D. R., and Rowland, F. S.: Dimethyl Disulfide (DMDS) and Dimethyl Sulfide (DMS) Emissions from Biomass Burning in Australia, Geophys. Res. Lett., 30, 10.1029/2003GL016967, 2003. Nadykto, A. B. and Yu, F.: Strong Hydrogen Bonding Between Atmospheric Nucleation Precursors and Common Organics, Chem. Phys. Lett., 435, 14–18, 2007. Nadykto, A. B., Yu, F., Jakovleva, M. V., Herb, J., and Xu, Y.: Amines in the Earth's Atmosphere: A Density Functional Theory Study of the Thermochemistry of Pre-Nucleation Clusters, Entropy, 13, 554–569, 2011. Nadykto, A. B., Herb, J., Yu, F., and Xu, Y.: Enhancement in the Production of Nucleating Clusters Due to Dimethylamine and Large Uncertainties in the Thermochemistry of Amine-Enhanced Nucleation, Chem. Phys. Lett., 609, 42–49, 2014. Nadykto, A. B., Herb, J., Yu, F., Xu, Y., and Nazarenko, E. S.: Estimating the Lower Limit of the Impact of Amines on Nucleation in the Earth's Atmosphere, Entropy, 17, 2764–2780, 2015. Napari, I., Kulmala, M., and Vehkamäki, H.: Ternary Nucleation of Inorganic Acids, Ammonia, and Water, J. Chem. Phys., 117, 8418–8425, 2002.Neese, F.: Software update: The ORCA program system—Version 5.0, WIREs Comput. Mol. Sci., 12, e1606, 10.1002/wcms.1606, 2022. Ning, A., Mao, L., Zhao, B., Zu, H., Shen, J., Zhao, Y., Li, J., Deng, X., Liu, L., Zhang, H., et al.: Marine upper-tropospheric rapid particle formation dominated by methanesulfonic acid, P. Natl. Acad. Sci. USA, 123, e2606521123, 2026.Olenius, T., Kupiainen-Määttä, O., Ortega, I. K., Kurtén, T., and Vehkamäki, H.: Free Energy Barrier in the Growth of Sulfuric Acid–Ammonia and Sulfuric Acid–Dimethylamine Clusters, J. Chem. Phys., 139, 084312, 10.1063/1.4819024, 2013.Ortega, I. K., Kupiainen, O., Kurtén, T., Olenius, T., Wilkman, O., McGrath, M. J., Loukonen, V., and Vehkamäki, H.: From quantum chemical formation free energies to evaporation rates, Atmos. Chem. Phys., 12, 225–235, 10.5194/acp-12-225-2012, 2012. Passananti, M., Zapadinsky, E., Zanca, T., Kangasluoma, J., Myllys, N., Rissanen, M. P., Kurtén, T., Ehn, M., Attoui, M., and Vehkamäki, H.: How Well Can We Predict Cluster Fragmentation Inside a Mass Spectrometer?, Chem. Commun., 55, 5946–5949, 2019.Pelucchi, C., Negri, E., Gallus, S., Boffetta, P., Tramacere, I., and La Vecchia, C.: Long-Term Particulate Matter Exposure and Mortality: A Review of European Epidemiological Studies, BMC Public Health, 9, 453, 10.1186/1471-2458-9-453, 2009. Perraud, V., Xu, J., Gerber, R. B., and Finlayson-Pitts, B.: Integrated Experimental and Theoretical Approach to Probe the Synergistic Effect of Ammonia in Methanesulfonic Acid Reactions with Small Alkylamines, Environ. Sci.: Processes Impacts, 22, 305–328, 2020. Pracht, P. and Grimme, S.: Calculation of Absolute Molecular Entropies and Heat Capacities Made Simple, Chem. Sci., 12, 6551–6568, 2021. Pracht, P., Bauer, C. A., and Grimme, S.: Automated and Efficient Quantum Chemical Determination and Energetic Ranking of Molecular Protonation Sites, J. Comput. Chem., 38, 2618–2631, 2017. Pracht, P., Bohle, F., and Grimme, S.: Automated Exploration of the Low-Energy Chemical Space with Fast Quantum Chemical Methods, Phys. Chem. Chem. Phys., 22, 7169–7192, 2020. Rasmussen, F. R., Kubečka, J., and Elm, J.: Contribution of Methanesulfonic Acid to the Formation of Molecular Clusters in the Marine Atmosphere, J. Phys. Chem. A, 126, 7127–7136, 2022. Riccobono, F., Schobesberger, S., Scott, C. E., Dommen, J., Ortega, I. K., Rondo, L., Almeida, J., Amorim, A., Bianchi, F., Breitenlechner, M., David, A., Downard, A., Dunne, E. M., Duplissy, J., Ehrhart, S., Flagan, R. C., Franchin, A., Hansel, A., Junninen, H., Kajos, M., Keskinen, H., Kupc, A., Kürten, A., Kvashin, A. N., Laaksonen, A., Lehtipalo, K., Makhmutov, V., Mathot, S., Nieminen, T., Onnela, A., Petäjä, T., Praplan, A. P., Santos, F. D., Schallhart, S., Seinfeld, J. H., Sipilä, M., Spracklen, D. V., Stozhkov, Y., Stratmann, F., Tomé, A., Tsagkogeorgas, G., Vaattovaara, P., Viisanen, Y., Vrtala, A., Wagner, P. E., Weingartner, E., Wex, H., Wimmer, D., Carslaw, K. S., Curtius, J., Donahue, N. M., Kirkby, J., Kulmala, M., Worsnop, D. R., and Baltensperger, U.: Oxidation Products of Biogenic Emissions Contribute to Nucleation of Atmospheric Particles, Science, 344, 717–721, 2014. Rosenfeld, D., Andreae, M. O., Asmi, A., Chin, M., de Leeuw, G., Donovan, D. P., Kahn, R., Kinne, S., Kivekäs, N., Kulmala, M., Lau, W., Schmidt, K. S., Suni, T., Wagner, T., Wild, M., and Quaas, J.: Global Observations of Aerosol-Cloud-Precipitation-Climate Interactions, Rev. Geophys., 52, 750–808, 2014. Rosenfeld, P. E., Henry, C. L., Dills, R. L., and Harrison, R. B.: Comparison of Odor Emissions from Three Different Biosolids Applied to Forest Soil, Water Air Soil Pollut., 127, 173–191, 2001. Schobesberger, S., Junninen, H., Bianchi, F., Lönn, G., Ehn, M., Lehtipalo, K., Dommen, J., Ehrhart, S., Ortega, I. K., Franchin, A., Nieminen, T., Riccobono, F., Hutterli, M., Duplissy, J., Almeida, J., Amorim, A., Breitenlechner, M., Downard, A. J., Dunne, E. M., Flagan, R. C., Kajos, M., Keskinen, H., Kirkby, J., Kupc, A., Kürten, A., Kurtén, T., Laaksonen, A., Mathot, S., Onnela, A., Praplan, A. P., Rondo, L., Santos, F. D., Schallhart, S., Schnitzhofer, R., Sipilä, M., Tomé, A., Tsagkogeorgas, G., Vehkamäki, H., Wimmer, D., Baltensperger, U., Carslaw, K. S., Curtius, J., Hansel, A., Petäjä, T., Kulmala, M., Donahue, N. M., and Worsnop, D. R.: Molecular Understanding of Atmospheric Particle Formation from Sulfuric Acid and Large Oxidized Organic Molecules, Proc. Natl. Acad. Sci. USA, 110, 17223–17228, 2013.Schütt, K. T., Kessel, P., Gastegger, M., Nicoli, K. A., Tkatchenko, A., and Müller, K.-R.: SchNetPack: A Deep Learning Toolbox For Atomistic Systems, J. Chem. Theory Comput., 15, 448–455, 10.1021/acs.jctc.8b00908, 2019.Schütt, K. T., Unke, O. T., and Gastegger, M.: Equivariant Message Passing for the Prediction of Tensorial Properties and Molecular Spectra, arXiv, 10.48550/arXiv.2102.03150, 2021.Schütt, K. T., Hessmann, S. S. P., Gebauer, N. W. A., Lederer, J., and Gastegger, M.: SchNetPack 2.0: A Neural Network Toolbox for Atomistic Machine Learning, J. Chem. Phys., 158, 10.1063/5.0138367, 2023. Shen, J., Xie, H.-B., Elm, J., Ma, F., Chen, J., and Vehkamaki, H.: Methanesulfonic Acid-Driven New Particle Formation Enhanced by Monoethanolamine: A Computational Study, Environ. Sci. Technol., 53, 14387–14397, 2019. Shen, J., Elm, J., Xie, H.-B., Chen, J., Niu, J., and Vehkamaki, H.: Structural Effects of Amines in Enhancing Methanesulfonic Acid-Driven New Particle Formation, Environ. Sci. Technol., 54, 13498–13508, 2020. Shen, J., Scholz, W., He, X.-C., Zhou, P., Marie, G., Wang, M., Marten, R., Surdu, M., Rörup, B., Baalbaki, R., Amorim, A., Ataei, F., Bell, D. M., Bertozzi, B., Brasseur, Z., Caudillo, L., Chen, D., Chu, B., Dada, L., Duplissy, J., Finkenzeller, H., Granzin, M., Guida, R., Heinritzi, M., Hofbauer, V., Iyer, S., Kemppainen, D., Kong, W., Krechmer, J. E., Kürten, A., Lamkaddam, H., Lee, C. P., Lopez, B., Mahfouz, N. G. A., Manninen, H. E., Massabò, D., Mauldin, R. L., Mentler, B., Müller, T., Pfeifer, J., Philippov, M., Piedehierro, A. A., Roldin, P., Schobesberger, S., Simon, M., Stolzenburg, D., Tham, Y. J., Tomé, A., Umo, N. S., Wang, D., Wang, Y., Weber, S. K., Welti, A., Wollesen de Jonge, R., Wu, Y., Zauner-Wieczorek, M., Zust, F., Baltensperger, U., Curtius, J., Flagan, R. C., Hansel, A., Möhler, O., Petäjä, T., Volkamer, R., Kulmala, M., Lehtipalo, K., Rissanen, M., Kirkby, J., El-Haddad, I., Bianchi, F., Sipilä, M., Donahue, N. M., and Worsnop, D. R.: High Gas-Phase Methanesulfonic Acid Production in the OH-Initiated Oxidation of Dimethyl Sulfide at Low Temperatures, Environ. Sci. Technol., 56, 13931–13944, 2022. Sipilä, M., Berndt, T., Petäjä, T., Brus, D., Vanhanen, J., Stratmann, F., Patokoski, J., Mauldin, R. L., III, Hyvärinen, A.-P., Lihavainen, H., and Kulmala, M.: The Role of Sulfuric Acid in Atmospheric Nucleation, Science, 327, 1243–1246, 2010. Sorooshian, A., Padró, L. T., Nenes, A., Feingold, G., McComiskey, A., Hersey, S. P., Gates, H., Jonsson, H. H., Miller, S. D., Stephens, G. L., Flagan, R. C., and Seinfeld, J. H.: On the Link Between Ocean Biota Emissions, Aerosol, and Maritime Clouds: Airborne, Ground, and Satellite Measurements off the Coast of California, Glob. Biogeochem. Cy., 23, 2009. Spicher, S., Plett, C., Pracht, P., Hansen, A., and Grimme, S.: Automated Molecular Cluster Growing for Explicit Solvation by Efficient Force Field and Tight-Binding Methods, J. Chem. Theory Comput., 18, 3174–3189, 2022.Spracklen, D. V., Carslaw, K. S., Kulmala, M., Kerminen, V.-M., Mann, G. W., and Sihto, S.-L.: The contribution of boundary layer nucleation events to total particle concentrations on regional and global scales, Atmos. Chem. Phys., 6, 5631–5648, 10.5194/acp-6-5631-2006, 2006. VanderGheynst, J. S., Cogan, D. J., DeFelice, P. J., Gossett, J. M., and Walker, L. P.: Effect of Process Management on the Emission of Organosulfur Compounds and Gaseous Antecedents from Composting Processes, Environ. Sci. Technol., 32, 3713–3718, 1998. Vanhanen, J., Mikkilä, J., Lehtipalo, K., Sipilä, M., Manninen, H., Siivola, E., Petäjä, T., and Kulmala, M.: Particle Size Magnifier for Nano-CN Detection, Aerosol Sci. Technol., 45, 533–542, 2011. Wilemski, G. and Wyslouzil, B. E.: Binary Nucleation Kinetics. I. Self-Consistent Size Distribution, J. Chem. Phys., 103, 1127–1136, 1995. Wu, H., Engsvang, M., Knattrup, Y., Kubecka, J., and Elm, J.: Improved Configurational Sampling Protocol for Large Atmospheric Molecular Clusters, ACS Omega, 8, 45065–45077, 2023.Wu, H., Knattrup, Y., Jensen, A. B., and Elm, J.: Cluster-to-particle transition in atmospheric nanoclusters, Aerosol Research, 2, 303–314, 10.5194/ar-2-303-2024, 2024.Yu, W., Baalbaki, R., Shen, J., DeVivo, J., Caudillo-Plath, L., Sommer, E., Heitto, A., Klebach, H., Russell, D. M., Amorim, A., Beckmann, H., Bhattacharyya, N., Blankenship, V., Chassaing, A., Cruz-Simbron, R., Dada, L., Jacobshagen, A., Judmaier, B., Kaniyodical Sebastian, M., Konrat, R., Krüger, T., Kunkler, F., Lietzke, C. J., Liu, L., Mauldin, R., Mentler, B., Morawiec, A., Rato, P., Rörup, B., Ruhl, S., Scholz, W., Simon, M., Stinchfield, A., Tomé, A., Tong, Y., Top, J., Umo, N. S., Unfer, G. R., Vettikkat, L., Weissbacher, J., Xenofontos, C., Yang, B., Zauner-Wieczorek, M., Zhang, J., Zheng, Z., Christoudias, T., Curtius, J., El Haddad, I., Flagan, R., Hansel, A., Harder, H., Kürten, A., Möhler, O., Petäjä, T., Volkamer, R., Winkler, P. M., Worsnop, D. R., Kulmala, M., Donahue, N. M., Kirkby, J., Yli-Juuti, T., Riipinen, I., He, X.-C., and Lehtipalo, K.: Impact of humidity on aerosol growth from methanesulfonic acid, Environ. Sci.-Atmos., 10.1039/d5ea00123d, 2026. Zapadinsky, E., Passananti, M., Myllys, N., Kurtén, T., and Vehkamäki, H.: Modeling on Fragmentation of Clusters Inside a Mass Spectrometer, J. Phys. Chem. A, 123, 611–624, 2018.Zhang, J. and Dolg, M.: ABCluster: The Artificial Bee Colony Algorithm for Cluster Global Optimization, Phys. Chem. Chem. Phys., 17, 24173–24181, 10.1039/C5CP04060D, 2015.Zhang, J. and Dolg, M.: Global Optimization of Clusters of Rigid Molecules Using the Artificial Bee Colony Algorithm, Phys. Chem. Chem. Phys., 18, 3003–3010, 10.1039/C5CP06313B, 2016. Zhao, B., Donahue, N. M., Zhang, K., Mao, L., Shrivastava, M., Ma, P.-L., Shen, J., Wang, S., Sun, J., Gordon, H., Tang, S., Fast, J., Wang, M., Gao, Y., Yan, C., Singh, B., Li, Z., Huang, L., Lou, S., Lin, G., Wang, H., Jiang, J., Ding, A., Nie, W., Qi, X., Chi, X., and Wang, L.: Global Variability in Atmospheric New Particle Formation Mechanisms, Nature, 631, 98–105, 2024.