Highly oxygenated organic molecule cluster decomposition in atmospheric pressure interface time-of-flight mass spectrometers
- 1Institute for Atmospheric and Earth System Research/Physics, Faculty of Science, University of Helsinki, P.O. Box 64, 00014 Helsinki, Finland
- 2Dipartimento di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
- 3Institute for Atmospheric and Earth System Research/Chemistry, Faculty of Science, University of Helsinki, P.O. Box 55, 00014 Helsinki, Finland
Correspondence: Tommaso Zanca (firstname.lastname@example.org)
Identification of atmospheric molecular clusters and measurement of their concentrations by atmospheric pressure interface time-of-flight (APi-TOF) mass spectrometers may be affected by systematic error due to possible decomposition of clusters inside the instrument. Here, we perform numerical simulations of decomposition in an APi-TOF mass spectrometers and formation in the atmosphere of a set of clusters which involve a representative kind of highly oxygenated organic molecule (HOM), with the molecular formula C10H16O8. This elemental composition corresponds to one of the most common mass peaks observed in experiments on ozone-initiated autoxidation of α-pinene. Our results show that decomposition is highly unlikely for the considered clusters, provided their bonding energy is large enough to allow formation in the atmosphere in the first place.
Recent developments in mass spectrometry have brought huge advancements to the field of atmospheric science. For example, mass spectrometers are now able to detect parts-per-quadrillion-level concentrations of both clusters and precursor vapors in atmospheric samples (Junninen et al., 2010; Jokinen et al., 2012), as well as directly explore the chemistry of new-particle formation (NPF) in the atmosphere (Kulmala et al., 2014; Almeida et al., 2013; Bianchi et al., 2016; Ehn et al., 2014; Thomas et al., 2016; Hogan and de la Mora, 2010).
One of the most common mass spectrometers used to measure online cluster composition and concentration in the atmosphere is the atmospheric pressure interface time-of-flight mass spectrometer (APi-TOF MS). It can detect naturally charged molecules and clusters, or it can be used in combination with a chemical ionization chamber (CI-APi-TOF) to detect neutral molecules and clusters. Unfortunately, the detection process in the APi-TOF MS involves energetic interactions between the carrier gas and the clusters, possibly leading to their decomposition and thus altering the measurement results. Specifically, the ions are guided and focused by an electric field inside the atmospheric pressure interface (APi) through a series of three vacuum chambers before arriving to the time-of-flight mass spectrometer itself. The pressure decreases between successive chambers, until the final value of 10−6 mbar is reached in the time of flight (TOF). Inside the APi, the clusters can collide with carrier gas molecules and can possibly be decomposed in the process.
Cluster decomposition in the APi-TOF MS is one of the main sources of uncertainty in the measurements. Lack of accuracy for a quantitative estimate of decomposition for example makes it difficult to draw definitive conclusions on the presence (or absence) of certain molecular clusters in the atmosphere. Often the absence of observations of specific clusters by the APi-TOF MS has led to speculation about decomposition inside the mass spectrometer (Olenius et al., 2013). Recently, a numerical model to study decomposition in the APi-TOF MS has been developed by Zapadinsky et al. (2018). The decomposition model has been tested on simple clusters involving one bisulfate anion and two sulfuric acid molecules, giving very good agreement with experimental results (Passananti et al., 2019).
The uncertainties on cluster concentration measurements by the APi-TOF MS and the lack of a comprehensive understanding of decomposition inside the instrument are the motivations of the present work. Here we use a theoretical model to study in detail the decomposition of clusters involving so-called highly oxygenated organic molecules (HOMs), which have recently been identified as a key contributor to NPF (Bianchi et al., 2019). HOMs are molecules formed in the atmosphere from volatile organic compounds (VOCs). Some VOCs with suitable functional groups can undergo an autoxidation process involving peroxy radicals, generating polyfunctional low-volatility vapors (i.e., HOMs) that subsequently condense onto pre-existing particles. HOMs thus contribute to secondary organic aerosol (SOA), which constitutes a significant fraction of the submicron organic aerosol mass and is known to affect the Earth's radiation balance (Jimenez et al., 2009; Donahue et al., 2009; Hallquist et al., 2009). Recent chamber experiments have shown that NPF induced by multicomponent systems, such as sulfuric acid, ammonia and HOMs, could correctly reproduce the NPF events observed in boreal forests (Lehtipalo et al., 2018).
Our study involves a specific kind of representative HOM (C10H16O8) in the APi. This elemental composition corresponds to one of the most common mass peaks observed in experiments on ozone-initiated autoxidation of α-pinene, which also fulfills the “HOM” definition of Bianchi et al. (2019). The precise molecular structure was adopted from Kurtén et al. (2016) and corresponds to the lowest-volatility structural isomer of the three C10H16O8 compounds investigated in that study. The structure of the molecule is shown in Fig. 1.
The main scope of this work is to determine to what extent we are able to perform measurements of atmospheric cluster concentrations using APi-TOF mass spectrometers. More specifically, we want to determine whether decomposition can possibly be responsible for the lack of observations of some HOM-containing clusters in an APi-TOF. Basically, the formation of clusters in the atmosphere (often referred to as “nucleation”) is driven by the bigger stability of the cluster with respect to its separated molecular (and sometimes ionic) components (fragments). The degree of stability is given by the energy difference between the fragments and the cluster in their (electronic and rotational–vibrational) ground states, which is called either reaction, binding or decomposition energy, and it is equal to the amount of energy necessary to decompose the cluster (or, vice versa, the amount of energy released by the clustering). Higher decomposition energy implies lower decomposition and evaporation rates, and thus higher net formation rates in the atmosphere. On the contrary, decomposition in the APi-TOF MS is enhanced when the decomposition energy decreases, since less energy is required to break the cluster. It is clear at this point that decomposition and new-particle formation (“nucleation”) have opposite dependences on the cluster decomposition energy. For a given temperature and vapor concentration, two situations are possible (Fig. 2):
There is a limited range of values for the decomposition energy which allows both decomposition in the APi-TOF and particle formation in the atmosphere.
The smallest decomposition energy that allows particle formation in the atmosphere is still too large to permit decomposition in the APi-TOF.
Here, we predict both an upper bound for decomposition energy necessary for decomposition in the APi-TOF, and a lower bound for new-particle formation in the atmosphere given realistic vapor concentrations.
The former can be evaluated using the numerical model
developed by Zapadinsky et al. (2018), while the latter is computed using the Atmospheric Cluster Dynamics Code (ACDC) (McGrath et al., 2012).
In the present work the simulations on decomposition have been performed by a new
C++ version of the decomposition code, which keeps the same basic algorithm but performs faster.
Another purpose of this work is to analyze the dynamics of complex clusters in the APi, showing where and why the decomposition events take place.
In Sect. 2 we present the decomposition model. In Sect. 3 we show the main results on HOM cluster decomposition and we compare them with ACDC simulations. Finally, in Sect. 4, we present the conclusions.
The model computes the probability of decomposition of a single negatively charged cluster from a large number (≈103) of independent stochastic realizations of its dynamics in the APi. The algorithm of the code is presented in the simplified flowchart shown in Fig. 3. We give only a brief description of the algorithm here; for full details see Zapadinsky et al. (2018). The dynamics start with the ionized cluster accelerating in the mass spectrometer under the effect of an electric field generated by the electrodes inside the APi. While moving, a random time interval to the next collision is computed from the cumulative distribution function 𝒫coll(t), which expresses the probability of encountering a collision after a time t:
Here, Υ is the collision frequency which depends on the cluster velocity v, while t=0 is the moment when the previous collision occurred. Simultaneously, another random time interval is computed for the decomposition event from a Poisson distribution, which corresponds to the time-dependent survival probability 𝒫surv after collision:
where 𝒫dec is the decomposition probability, and k, the decomposition rate constant, is the inverse of the statistical average of decomposition time, which depends on the cluster excess energy ΔE beyond its decomposition energy threshold. The decomposition time is interpreted as the time the cluster spends intact before decomposing. The decomposition rate constant is derived from phase space theory (PST) of chemical reactions. Basically, it is determined by the ratio of the densities of states of products to the densities of states of reactant, exploiting the detailed balance condition (Zapadinsky et al., 2018).
At this point, if the time required by the cluster to escape from the simulated region of the mass spectrometer is less than both the collision and decomposition times, the single realization is completed and the count of intact clusters is increased by one. A second possibility is that the decomposition time is smaller than collision and escape times, in which case the count of decomposed clusters is increased by 1. Finally, as a third possibility, the collision time may be the smallest, in which case a collision between the cluster and a carrier gas molecule takes place.
The dynamics of the collision are described by a stochastic process, where random velocities for the carrier gas molecules are computed from the Maxwell–Boltzmann distribution. In this case, the motion of the gas molecules is considered as solely translational. During the collision, the kinetic energy of the two colliding objects is partially transferred to the internal degrees of freedom of the cluster (vibrational and rotational modes), according to the microcanonical ensemble approach: any configuration of the system with the same energy is equally probable. The fundamental relation governing these statistics is given by the proportionality between the probability density function of energies for two different interacting degrees of freedom and their densities of states:
where ϵi is the energy of mode i, ρi is the density of states of mode i and the Dirac delta function δ ensures conservation of total energy E. After integrating Eq. (3) over all the possible values that ϵb can assume, we get the probability density function for ϵa:
Notice that this formula holds at equilibrium. This assumption is justified from the fact that the energy transfer takes place at timescales much shorter than the timescale between two consecutive collisions. Specifically, vibrational–vibrational energy exchange takes place in about 10−13–10−12 s, rotational–vibrational exchange in 10−11 s and collisions every 10−9–10−5 s.
After collision, the energy is then redistributed between rotational and vibrational modes following the same principle, and the dynamics continue with the acceleration of the cluster in the electric field.
In region IV, where the quadrupole is located, the cluster is subjected to additional acceleration in the transverse directions, because of the alternating electric field generated by the quadrupole. The dynamics is determined by the renowned Mathieu equation (Miller and Denton, 1986):
with y one of the two transverse displacements from the quadrupole axis, t the time, e the elementary charge, U the DC component, V the AC amplitude of the electric potential, m the cluster mass, r0 the half-distance between quadrupole rods and ω the AC angular frequency.
The computation of the density of vibrational and rotational states requires the knowledge of the corresponding energy levels, which are computed by the quantum chemistry program Gaussian, using the PM7 semi-empirical method, which is the newest and generally best semi-empirical method available in Gaussian (Frisch et al., 2016). We note that while semi-empirical methods are unable to accurately model the energetics of molecular clustering, this is not a problem in the present study, as the binding/decomposition energy is treated here as a freely variable parameter (i.e., the unreliable PM7 binding energy is not actually used). The purpose of the PM7 optimizations and frequency calculations is simply to provide a qualitatively correct distribution of rotational and vibrational energy levels. As vibrational modes with wavenumbers above 2500 cm−1 are strongly underestimated by the PM7 method, we have further rescaled the frequencies following the relation suggested in Rozanska et al. (2014):
The vibrational density of states is then computed using the harmonic approximation, which is most likely the biggest source of error in the evaluation of the cluster survival probability. (See Sect. 3 for a sensitivity analysis on the effect of varying the vibrational frequencies.)
3.1 Decomposition in the APi
The simulation of decomposition involves only a portion of the total length of the APi-TOF mass spectrometer. Specifically, the simulations take place in the most critical region, ranging from the end of the first to the end of the second chamber of the APi, where the pressure values make decomposition possible: at the beginning of the first chamber the collisions are not energetic enough, while in the third chamber the carrier gas is so sparse that no collision happens (Zapadinsky et al., 2018).
In Fig. 4 we have sketched the simulated section of the instrument. We subdivide the simulated section into five regions with different longitudinal electric field and pressure values. Region I is at the end of the first chamber; region II defines the interface between the first and the second chambers; and regions III, IV and V are located in the second chamber of the APi. The pressure in the first region is P1, while in regions III, IV and V it is P2. In region II, where the skimmer is located, the pressure changes continuously from P1 to P2. The electric field takes different values in regions I, III, IV and V, while in II it is set to zero. The voltage and pressure configuration inside the APi chambers, used in the CLOUD10 experiments (Lehtipalo et al., 2018), are shown in Table 1, where VDC and VAC are the direct and alternating components of the quadrupole electric potential in region IV, while Ω is its radio frequency. We will adopt the same parameter configuration in our simulations. It is important to notice here that the electric field takes very different values in different regions, varying over 3 orders of magnitude, which will greatly diversify the probability of decomposition at different locations.
In this study the cluster decomposition energy Ef is treated as a free parameter. This allows us to explore the behavior of cluster survival probability as a function of Ef over a large energy range. Moreover, keeping Ef as a variable is also useful since its computation by quantum chemistry calculations (especially low-level methods such as PM7 used here) is affected by large errors.
The clusters studied here are formed by one bisulfate anion, one to two sulfuric acid molecules and one to three HOM molecules. The clusters are assumed to decompose by losing one HOM, as follows:
with , and HOM10 is the structure shown in Fig. 1. Our specific choice for HOM is not unique: this is one among many potential HOM produced in the α-pinene +O3 reaction. The compound chosen is broadly representative of autoxidation products as it contains both hydroperoxide, ketone and carboxylic acid groups. The clusters have been constructed by first maximizing H-bonds between the core ion and other molecules, and then maximizing other H-bonds without creating too much strain. We note that the cluster conformers obtained in this fashion are unlikely to correspond precisely to the global energy (or free energy) minima. However, this mainly affects the computed binding energy – which (as pointed out earlier) is not actually used in this study. For the purposes of generating a representative ensemble of vibrational and rotational energy levels, the conformer generation approach used here is adequate.
For each kind of cluster we performed several simulations at different Ef. The simulated carrier gas is air, modeled to consist of 80 % nitrogen and 20 % oxygen.
The final survival probabilities Psurv are shown in Fig. 5. As expected, Psurv is monotonic with respect to Ef. From these results we can identify a boundary at Ef≈25 kcal mol−1 beyond which decomposition is highly unlikely. The weakest bonds between non-hydrogen atoms in HOM molecules are likely to be the O−O bonds of peroxide or hydroperoxide groups. These typically have bond dissociation energies around 35–50 kcal mol−1 (Bach et al., 1996; Schweitzer-Chaput et al., 2015, 2017). Vinyl peroxide systems have much lower dissociation energies, but (with the exception of short-lived vinyl hydroperoxides generated in ozonolysis) these are unlikely to form in gas-phase oxidation processes. The threshold cluster decomposition energies reported here are thus much lower than the energies required to dissociate even the weakest covalent bonds in the HOM molecules.
Larger clusters are found to decompose more easily. This can be understood looking at the probability density function (PDF) of vibrational energy at equilibrium at temperature T=300 K (Fig. 6a). The distributions have been computed multiplying the density of states of vibrational modes by the Boltzmann factor, and then renormalized:
Larger clusters contain a higher number of bonds, which, for fixed decomposition energy, will each increase the internal energy of the cluster, because of the equipartition theorem. This leads to a higher probability of exceeding the decomposition energy, even at equilibrium conditions, i.e., in the absence of an electric field. Suppose, for example, Ef = 20 kcal mol−1. Figure 6a shows that the smallest cluster has roughly 50 % probability of having an internal energy higher than Ef, while the biggest cluster is certainly at higher energy than Ef, with an average near 50 kcal mol−1. This effect is in part counterbalanced by the behavior of the decomposition rate constant as a function of cluster size. As we can see in Fig. 6b, the decomposition rate at fixed internal energy tends to be smaller for larger clusters. Indeed, a larger number of bonds reduces the probability of finding a large quantity of energy in a single bond.
The model requires as input parameters the vibrational and rotational frequencies of both the reactant and products from an external quantum chemistry program (in our case Gaussian 16). These input data are affected by errors that depend on the level of accuracy of the computational method. It is therefore interesting to study how errors in the inputs affect the survival probability. For this reason, we generated new frequencies, adding random noise to the original ones. The new frequencies are generated from normal distributions centered at the original frequencies, with a standard deviation equal to 20 % of their values. Figure 7 shows the deviations of the survival probability for 10 different sets of random input frequencies. Here we can see that, even with large errors in the input parameters (20 %), the final results change by no more than 10 %. These results demonstrate the low sensitivity of the model on deviations in the input frequencies, thus validating our use of rather crude approaches for both the quantum chemical calculations, and the conformer generation.
Let us now analyze the locations where collisions take place in the APi. In particular, we are interested in the highest-energy collisions, which lead to decomposition. We dub them as fatal collisions. For this study we used the voltage configuration and pressures used in CLOUD10 and we collected data from 5000 realizations of cluster dynamics at Ef=17.88 kcal mol−1, corresponding to a survival probability Psurv≈50 % (see Fig. 5).
Figure 8a shows the average number of total and fatal collisions per single cluster evolution in different regions of the APi. The average number of fatal collisions is computed by dividing the total number of fatal collisions by the number of realizations (dark bars in Fig. 8a), and it can be interpreted as the probability of the cluster to decompose in the correspondent region in a complete realization. Summing up these values of all five regions gives the total probability of decomposition. Figure 8b instead shows the ratio between fatal and total collisions, which can be interpreted as the probability of the cluster to decompose after experiencing a collision in the correspondent region. Therefore, Fig. 8b identifies the regions where the collisions are most energetic (regions III and V).
The cluster dynamics starts in the first region (end of first APi chamber), where the high pressure (201 Pa) generates a large number of collisions (≈200 collisions per realization). Here the ratio between fatal collisions and total collisions is about 0.07 % (Fig. 8b), a sign of low-energy collisions despite the very strong electric field (7189 V m−1). When the cluster moves to region II, in the skimmer, we see a drop in collisions, reflecting the negative pressure gradient. Moreover, the absence of electric field in this region makes fatal collisions very unlikely to happen (≈0.017 % probability). In region III, at the beginning of the second APi chamber, the pressure is low (2.96 Pa), and the cluster experiences few collisions. However, the strong electric field (804.9 V m−1) leads to cluster acceleration and an increase in fatal collisions, about 0.4 % of total collisions. The passage to region IV is characterized by a high speed of the cluster acquired in the previous region. This leads to a large number of decompositions, as shown in Fig. 8a. Subsequently, the cluster decelerates because of the very weak electric field in this region (18.04 V m−1). Since the cluster moves for a long distance at low velocities in region IV, the average probability of decomposition per collision drops to ≈0.05 %. Finally, in region V, the cluster accelerates again due to a strong electric field (1104 V m−1), which is visible from the large increase in the probability of fatal collisions (≈2 %). The cumulative survival probability of clusters as they travel through regions I to V is illustrated by the red dots in Fig. 9. The decomposition probabilities in each respective region are indicated by the histogram bars.
Although the decomposition is mainly induced by the voltages applied between regions II and III as shown in Passananti et al. (2019), these results demonstrate that the decomposition takes place mainly in region IV (Fig. 8a). This is due to the large acceleration acquired by the cluster in region III, with subsequent possible collision at high energy in the beginning of the next region. These results are complementary to previous data reported in the literature (Passananti et al., 2019; Lopez-Hilfiker et al., 2016) and help to define the terms “declustering voltages” as the voltages that induce decomposition (between regions II and III) and the “declustering region” as the region of the APi where the decomposition takes place (region IV).
In region IV a quadrupole is used as ion guide in order to focus the cluster trajectories, through an alternating transverse electric field (Miller and Denton, 1986). For this reason, the cluster is also subjected to a transverse acceleration in addition to the longitudinal one, increasing the velocity and, consequently, the probability of decomposition. In order to separate the contributions to decomposition from the longitudinal and the transverse electric fields, we performed additional simulations in which the quadrupole electric potential was set to VAC=0. In this situation the motion in the transverse directions is negligible (as in the other regions of the APi). The increase in the final survival probabilities can be seen in Fig. 10a, which shows a minor deviation from the case with VAC=200 V. More specifically, Fig. 10b shows the slight reduction in decomposition probability within the quadrupole region when VAC=0. Hence, with the present voltage configuration, the quadrupole voltage does not affect the probability of decomposition significantly.
3.2 Formation in the atmosphere
In parallel to the simulations on decomposition in the APi, we used the ACDC code (McGrath et al., 2012) to compute the concentrations of HOM10 needed to provide a significant enhancement of new-particle formation rate in the atmosphere (details on ACDC simulations provided in Appendix A). Specifically, at given H2SO4 monomer concentration and HOM cluster decomposition energy Ef, we computed the concentration of HOM10 at which the new-particle formation rate J is increased by 1 cm−3 s−1 with respect to the pure sulfuric acid system (i.e., , with ΔJ=1 cm−3 s−1).
Figure 11a is useful to understand the framework of the simulations: the scheme shows the clusters involved in the ACDC simulations: along the x axis we have clusters with an increasing number of sulfuric acid molecules, while along the y axis we have clusters with an increasing number of HOM10. Clusters grow through collisions with bisulfate anions, sulfuric acid and HOM10 molecules, until they reach a size bigger than three HOM10 or three sulfuric acid (SA) molecules. At this point the clusters exit the simulation region and they are accounted as outgrowing (stable) clusters. ACDC code computes the rate of formation and evaporation of all the clusters in the simulation region. We are interested in the outgoing fluxes (blue and red arrows in Fig. 11a), which are the ones that contribute to new-particle formation. In particular, we are interested in the enhancement of the total outgoing flux (formation rate) when HOM10 molecules are added to a system of sulfuric acid particles. The increment of the formation rate (ΔJ) is represented by the red arrows, and we want to find the concentration of HOM10 that corresponds to ΔJ=1 cm−3 s−1. This value is indicative, and it serves as a reference for a reasonable NPF process in the atmosphere (Lehtipalo et al., 2018). The concentration of is kept constant at 700 cm−3, which corresponds to the steady-state concentration at the ion formation rate Jions=4 cm−3 s−1, reproducing the galactic cosmic ray (GCR) conditions (Kirkby et al., 2016).
Different simulations have been performed for different sulfuric acid monomer concentrations and different decomposition energies Ef. Figure 11b shows the dependence of HOM10 concentration on the SA monomer concentration, in order to maintain ΔJ=1 cm−3 s−1, for different decomposition energies.
As we can see, the H2SO4 and HOM10 concentrations decrease when the decomposition energy increases, reaching experimental conditions (e.g., non-nitrate HOM10 and H2SO4 concentrations measured in the CLOUD experiment; Lehtipalo et al., 2018) at Ef>30 kcal mol−1, where HOM10 and H2SO4 concentrations do not exceed 108 cm−3. At this decomposition energy, decomposition in APi is negligible (Psurv≈1 in Fig. 5). Thus we can conclude that for the case of SA–HOM clusters, rapid formation in the atmosphere (given typical vapor concentrations) and significant decomposition in the APi are mutually incompatible situations (case b in Fig. 2).
In this work, we have presented the numerical results on decomposition inside an APi-TOF MS instrument of a specific class of atmospheric clusters that involve sulfuric acid and HOM molecules. A previously reported low-volatility α-pinene ozonolysis product, with the molecular formula C10H16O8, was used as a representative HOM.
There are three main results from our simulations. First, decomposition of HOM clusters in the APi requires a range of cluster decomposition energies which is incompatible with efficient cluster formation in the atmosphere given sub-parts-per-billion vapor concentrations. This result has been obtained by computing the cluster survival probability in the APi as a function of its decomposition energy, and then comparing the range of energies leading to decomposition with those needed to obtain a reasonably high new-particle formation rate in atmospheric conditions. The two ranges have no overlap – the highest energy allowing decomposition differs from the lowest energy allowing atmospheric new-particle formation by roughly 10 kcal mol−1. Observations of SA–HOM clusters in CLOUD experiments (Lehtipalo et al., 2018) validate the results of our model.
Our second main result is the identification of the locations of the highest-energy collisions that lead to decomposition of SA–HOM clusters in the APi. The simulation shows that they are mainly localized in the quadrupole region (IV). This is due to the large velocity of the cluster acquired in the previous region (III), caused by the strong electric field. Moreover, our results show that the alternating field of the quadrupole in region IV has only a minor effect on decomposition probability.
As a third main result, we have shown that the model displays low sensitivity to changes in cluster vibrational and rotational frequencies. This feature allows us to also obtain reliable results with low-level input data (e.g., semi-empirical methods and limited configurational sampling).
This study was performed for a small set of HOM clusters with particular functional groups (e.g., hydroperoxides, ketones), but in future similar simulations could be performed for other types of clusters, for example to assess whether they should be detectable in an APi-TOF MS or not, or whether they could possibly affect new-particle formation, despite not being directly detected.
The configurational sampling of the H2SO4 molecule and molecular clusters has been performed as described by Kubečka et al. (2019) to obtain their thermodynamic (T=300 K) properties at a high level of theory (DLPNO-CCSD(T)/aug-cc-pVTZ//ωB97X-D/6-31++G(d,p) (Kubečka et al., 2019; Myllys et al., 2016; Riplinger and Neese, 2013; Riplinger et al., 2013, 2016; Chai and Head-Gordon, 2008; Dunning, 1989; Kendall et al., 1992; Liakos et al., 2015)).
The formation Gibbs free energies were used as input for the Atmospheric Cluster Dynamic Code (McGrath et al., 2012), which uses them to obtain cluster stabilities (i.e., evaporation rates). Bisulfate ions in the atmosphere are formed by collisions of “generic” anions and sulfuric acid H2SO4 (i.e., as H2SO4 is practically the strongest acid in the atmosphere, collision of H2SO4 with almost any anion will tend to produce ).
In the simulation, we assume a constant concentration of bisulfate ions equal to 700 cm−3, which approximately corresponds to the steady-state anion concentration at a representative total ion pair formation rate due to galactic cosmic rays (GCRs) of Jions=4 cm−3 s−1 (Kirkby et al., 2016). We assume that formation of (outgrowing cluster) or larger clusters leads to irreversible new-particle formation. Thus, at a specific H2SO4 monomer concentration, we simulated new-particle formation of sulfuric acid aerosols until a steady state had been reached. For each monomer concentration, the new-particle formation rate Jsa were noted down. (For simplicity, we neglect losses due to ion–ion recombination, wall loss, dilution or coagulation, as our focus is on comparing the relative effect of HOM.)
As described in Kubečka et al. (2019), we performed “weak” configurational sampling of molecules HOM10, H2SO4 and molecular clusters to obtain representative thermodynamic properties of these clusters at a low level of theory (PM7, Stewart, 2013). By “weak” configurational sampling we mean that just one conformation of HOM10 is included for generating of HOM-containing clusters. However, as only the rotational and vibrational energy levels from the PM7 calculations are actually used, and not the electronic energies, this approach is sufficient for our purposes. We further assume that the evaporation rate of H2SO4 from clusters is identical to the evaporation rate of H2SO4 from the previously calculated corresponding clusters. The evaporation rate γ of HOM10 from is calculated as
where kB is Boltzmann's constant, T is temperature, β is the collision rate of HOM10 and , pref is the reference pressure, and ΔG is the free energy change of the reaction
ΔG is calculated as
where Ef is the parametrized value of the decomposition energy (varying between 15 and 30 kcal mol−1 as shown), ΔGPM7 and are the reaction Gibbs free energy and electronic energy (respectively) computed at the PM7 level of theory, and corresponds to the thermal (and zero-point vibrational) correction to the reaction Gibbs free energy at the PM7 level.
Again, in the same way as in the decomposition model, we used the evaporation rates as input for the ACDC program assuming or to be the outgrowing clusters. Next, at given H2SO4 monomer concentration, and for a given value of Ef, we searched for the HOM10 concentration at which the new-particle formation rate is increased by 1 cm−3 s−1 (i.e., cm−3 s−1). The resulting graph is shown in the main article in Fig. 11b.
The above computations involved the following programs: the ABCluster program (Zhang and Dolg, 2015, 2016), the GFN2-xTB program (Grimme et al., 2017; Bannwarth et al., 2019), the Gaussian 16 program Revision A.03 (Frisch et al., 2016), the GoodVibes program (Funes-Ardois and Paton, 2016) and the Orca program version 18.104.22.168 (Neese, 2012).
TK and HV conceived the project. TZ and EZ wrote the code of the decomposition model. TZ performed simulations of the decomposition model, carried out the data analysis and wrote the first draft of the paper. JK performed simulations of the ACDC code, produced the relative figures and wrote the supplementary information on ACDC simulations. MP helped in planning the project and interpreting the results in relation to experimental data. All authors contributed to writing the paper.
The authors declare that they have no conflict of interest.
We acknowledge the ERC, the ATMATH Project and the Academy of Finland for funding, and the CSC-IT Center for Science in Espoo, Finland, for computational resources. We would like to thank Federico Bianchi and Siegfried Schobesberger for valuable discussions, and Chao Yan for providing CLOUD campaign data.
This research has been supported by the ERC (DAMOCLES (grant no. 692891)), the ATMATH Project and the Academy of Finland.
Open access funding provided by Helsinki University Library.
This paper was edited by Hartmut Herrmann and reviewed by two anonymous referees.
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, https://doi.org/10.1038/nature12663, 2013. a
Bach, R. D., Ayala, P. Y., and Schlegel, H. B.: A Reassessment of the Bond Dissociation Energies of Peroxides. Anab InitioStudy, J. Am. Chem. Soc., 118, 12758–12765, https://doi.org/10.1021/ja961838i, 1996. a
Bannwarth, C., Ehlert, S., and Grimme, S.: GFN2-xTB–An accurate and broadly parametrized self-consistent tight-binding quantum chemical method with multipole electrostatics and density-dependent dispersion contributions, J. Chem. Theory Comput., 15, 1652–1671, 2019. a
Bianchi, F., Trostl, 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., Kurten, A., Manninen, H. E., Munch, S., Perakyla, O., Petaja, 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, https://doi.org/10.1126/science.aad5456, 2016. a
Bianchi, F., Kurtén, T., Riva, M., Mohr, C., Rissanen, M. P., Roldin, P., Berndt, T., Crounse, J. D., Wennberg, P. O., Mentel, T. F., Wildt, J., Junninen, H., Jokinen, T., Kulmala, M., Worsnop, D. R., Thornton, J. A., Donahue, N., Kjaergaard, H. G., and Ehn, M.: Highly Oxygenated Organic Molecules (HOM) from Gas-Phase Autoxidation Involving Peroxy Radicals: A Key Contributor to Atmospheric Aerosol, Chem. Rev., 119, 3472–3509, https://doi.org/10.1021/acs.chemrev.8b00395, 2019. a, b
Chai, J.-D. and Head-Gordon, M.: Long-range corrected hybrid density functionals with damped atom-atom dispersion corrections, Phys. Chem. Chem. Phys., 10, 6615–6620, 2008. a
Donahue, N. M., Robinson, A. L., and Pandis, S. N.: Atmospheric organic particulate matter: From smoke to secondary organic aerosol, Atmos. Environ., 43, 94–106, https://doi.org/10.1016/j.atmosenv.2008.09.055, 2009. a
Dunning, T. H.: Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen, J. Chem. Phys., 90, 1007–1023, 1989. a
Ehn, M., Thornton, J. A., Kleist, E., Sipilä, M., Junninen, H., Pullinen, I., Springer, M., Rubach, F., Tillmann, R., Lee, B., Lopez-Hilfiker, F., Andres, S., Acir, I.-H., Rissanen, M., Jokinen, T., Schobesberger, S., Kangasluoma, J., Kontkanen, J., Nieminen, T., Kurtén, T., Nielsen, L. B., Jørgensen, S., Kjaergaard, H. G., Canagaratna, M., Maso, M. D., Berndt, T., Petäjä, T., Wahner, A., Kerminen, V.-M., Kulmala, M., Worsnop, D. R., Wildt, J., and Mentel, T. F.: A large source of low-volatility secondary organic aerosol, Nature, 506, 476–479, https://doi.org/10.1038/nature13032, 2014. a
Frisch, M. J., Trucks, G. W., Schlegel, H. B., Scuseria, G. E., Robb, M. A., Cheeseman, J. R., Scalmani, G., Barone, V., Petersson, G. A., Nakatsuji, H., Li, X., Caricato, M., Marenich, A. V., Bloino, J., Janesko, B. G., Gomperts, R., Mennucci, B., Hratchian, H. P., Ortiz, J. V., Izmaylov, A. F., Sonnenberg, J. L., Williams-Young, D., Ding, F., Lipparini, F., Egidi, F., Goings, J., Peng, B., Petrone, A., Henderson, T., Ranasinghe, D., Zakrzewski, V. G., Gao, J., Rega, N., Zheng, G., Liang, W., Hada, M., Ehara, M., Toyota, K., Fukuda, R., Hasegawa, J., Ishida, M., Nakajima, T., Honda, Y., Kitao, O., Nakai, H., Vreven, T., Throssell, K., Montgomery, Jr., J. A., Peralta, J. E., Ogliaro, F., Bearpark, M. J., Heyd, J. J., Brothers, E. N., Kudin, K. N., Staroverov, V. N., Keith, T. A., Kobayashi, R., Normand, J., Raghavachari, K., Rendell, A. P., Burant, J. C., Iyengar, S. S., Tomasi, J., Cossi, M., Millam, J. M., Klene, M., Adamo, C., Cammi, R., Ochterski, J. W., Martin, R. L., Morokuma, K., Farkas, O., Foresman, J. B., and Fox, D. J.: Gaussian 16, Revision A.03, gaussian Inc., Wallingford, CT, USA, 2016. a, b
Grimme, S., Bannwarth, C., and Shuskov, 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 (), J. Chem. Theory Comput., 13, 1989–2009, 2017. a
Hallquist, M., Wenger, J. C., Baltensperger, U., Rudich, Y., Simpson, D., Claeys, M., Dommen, J., Donahue, N. M., George, C., Goldstein, A. H., Hamilton, J. F., Herrmann, H., Hoffmann, T., Iinuma, Y., Jang, M., Jenkin, M. E., Jimenez, J. L., Kiendler-Scharr, A., Maenhaut, W., McFiggans, G., Mentel, Th. F., Monod, A., Prévôt, A. S. H., Seinfeld, J. H., Surratt, J. D., Szmigielski, R., and Wildt, J.: The formation, properties and impact of secondary organic aerosol: current and emerging issues, Atmos. Chem. Phys., 9, 5155–5236, https://doi.org/10.5194/acp-9-5155-2009, 2009. a
Jimenez, J. L., Canagaratna, M. R., Donahue, N. M., Prevot, A. S. H., Zhang, Q., Kroll, J. H., DeCarlo, P. F., Allan, J. D., Coe, H., Ng, N. L., Aiken, A. C., Docherty, K. S., Ulbrich, I. M., Grieshop, A. P., Robinson, A. L., Duplissy, J., Smith, J. D., Wilson, K. R., Lanz, V. A., Hueglin, C., Sun, Y. L., Tian, J., Laaksonen, A., Raatikainen, T., Rautiainen, J., Vaattovaara, P., Ehn, M., Kulmala, M., Tomlinson, J. M., Collins, D. R., Cubison, M. J., Dunlea, J., Huffman, J. A., Onasch, T. B., Alfarra, M. R., Williams, P. I., Bower, K., Kondo, Y., Schneider, J., Drewnick, F., Borrmann, S., Weimer, S., Demerjian, K., Salcedo, D., Cottrell, L., Griffin, R., Takami, A., Miyoshi, T., Hatakeyama, S., Shimono, A., Sun, J. Y., Zhang, Y. M., Dzepina, K., Kimmel, J. R., Sueper, D., Jayne, J. T., Herndon, S. C., Trimborn, A. M., Williams, L. R., Wood, E. C., Middlebrook, A. M., Kolb, C. E., Baltensperger, U., and Worsnop, D. R.: Evolution of Organic Aerosols in the Atmosphere, Science, 326, 1525–1529, https://doi.org/10.1126/science.1180353, 2009. a
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, https://doi.org/10.5194/acp-12-4117-2012, 2012. a
Junninen, H., Ehn, M., Petäjä, T., Luosujärvi, L., Kotiaho, T., Kostiainen, R., Rohner, U., Gonin, M., Fuhrer, K., Kulmala, M., and Worsnop, D. R.: A high-resolution mass spectrometer to measure atmospheric ion composition, Atmos. Meas. Tech., 3, 1039–1053, https://doi.org/10.5194/amt-3-1039-2010, 2010. a
Kendall, A. R., Dunning, T. H., and Harrison, R. J.: Electron affinities of the first-row atoms revisited. Systematic basis sets and wave functions, J. Chem. Phys., 96, 6796–6806, 1992. a
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, https://doi.org/10.1038/nature17953, 2016. a, b
Kubečka, J., Besel, V., Kurtén, T., Myllys, N., and Vehkamäki, H.: Configurational sampling of noncovalent (atmospheric) molecular clusters: sulfuric acid and guanidine, J. Phys. Chem. A, 123, 6022–6033, 2019. a, b, c
Kulmala, M., Petäjä, T., Ehn, M., Thornton, J., Sipilä, M., Worsnop, D., and Kerminen, V.-M.: Chemistry of Atmospheric Nucleation: On the Recent Advances on Precursor Characterization and Atmospheric Cluster Composition in Connection with Atmospheric New Particle Formation, Annu. Rev. Phys. Chem., 65, 21–37, https://doi.org/10.1146/annurev-physchem-040412-110014, 2014. a
Kurtén, T., Tiusanen, K., Roldin, P., Rissanen, M., Luy, J.-N., Boy, M., Ehn, M., and Donahue, N.: α-Pinene Autoxidation Products May Not Have Extremely Low Saturation Vapor Pressures Despite High O:C Ratios, J. Phys. Chem. A, 120, 2569–2582, https://doi.org/10.1021/acs.jpca.6b02196, 2016. a
Lehtipalo, K., Yan, C., Dada, L., Bianchi, F., Xiao, M., Wagner, R., Stolzenburg, D., Ahonen, L. R., Amorim, A., Baccarini, A., Bauer, P. S., Baumgartner, B., Bergen, A., Bernhammer, A.-K., Breitenlechner, M., Brilke, S., Buchholz, A., Mazon, S. B., Chen, D., Chen, X., Dias, A., Dommen, J., Draper, D. C., Duplissy, J., Ehn, M., Finkenzeller, H., Fischer, L., Frege, C., Fuchs, C., Garmash, O., Gordon, H., Hakala, J., He, X., Heikkinen, L., Heinritzi, M., Helm, J. C., Hofbauer, V., Hoyle, C. R., Jokinen, T., Kangasluoma, J., Kerminen, V.-M., Kim, C., Kirkby, J., Kontkanen, J., Kürten, A., Lawler, M. J., Mai, H., Mathot, S., Mauldin, R. L., Molteni, U., Nichman, L., Nie, W., Nieminen, T., Ojdanic, A., Onnela, A., Passananti, M., Petäjä, T., Piel, F., Pospisilova, V., Quéléver, L. L. J., Rissanen, M. P., Rose, C., Sarnela, N., Schallhart, S., Schuchmann, S., Sengupta, K., Simon, M., Sipilä, M., Tauber, C., Tomé, A., Tröstl, J., Väisänen, O., Vogel, A. L., Volkamer, R., Wagner, A. C., Wang, M., Weitz, L., Wimmer, D., Ye, P., Ylisirniö, A., Zha, Q., Carslaw, K. S., Curtius, J., Donahue, N. M., Flagan, R. C., Hansel, A., Riipinen, I., Virtanen, A., Winkler, P. M., Baltensperger, U., Kulmala, M., and Worsnop, D. R.: Multicomponent new particle formation from sulfuric acid, ammonia, and biogenic vapors, Science Advances, 4, eaau5363, https://doi.org/10.1126/sciadv.aau5363, 2018. a, b, c, d, e
Liakos, D. G., Sparta, M., Kesharwani, M. K., Martin, J. M. L., and Neese, F.: Exploring the accuracy limits of local pair natural orbital coupled-cluster theory, J. Chem. Theory Comput., 11, 1525–1539, 2015. a
Lopez-Hilfiker, F. D., Iyer, S., Mohr, C., Lee, B. H., D'Ambro, E. L., Kurtén, T., and Thornton, J. A.: Constraining the sensitivity of iodide adduct chemical ionization mass spectrometry to multifunctional organic molecules using the collision limit and thermodynamic stability of iodide ion adducts, Atmos. Meas. Tech., 9, 1505–1512, https://doi.org/10.5194/amt-9-1505-2016, 2016. a
McGrath, M. J., Olenius, T., Ortega, I. K., Loukonen, V., Paasonen, P., Kurtén, T., Kulmala, M., and Vehkamäki, H.: Atmospheric Cluster Dynamics Code: a flexible method for solution of the birth-death equations, Atmos. Chem. Phys., 12, 2345–2355, https://doi.org/10.5194/acp-12-2345-2012, 2012. a, b, c
Myllys, N., Elm, J., Halonen, R., Kurtén, T., and Vehkamäki, H.: Coupled cluster evaluation of the stability of atmospheric acid–base clusters with up to 10 molecules, J. Phys. Chem. A, 120, 621–630, 2016. a
Neese, F.: The ORCA program system, WIREs Comput. Mol. Sci., 2, 73–78, 2012. a
Olenius, T., Schobesberger, S., Kupiainen-Määttä, O., Franchin, A., Junninen, H., Ortega, I. K., Kurtén, T., Loukonen, V., Worsnop, D. R., Kulmala, M., and Vehkamäki, H.: Comparing simulated and experimental molecular cluster distributions, Faraday Discuss., 165, 75–89, https://doi.org/10.1039/c3fd00031a, 2013. a
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, https://doi.org/10.1039/c9cc02896j, 2019. a, b, c
Riplinger, C., Sandhoefer, B., Hansen, A., and Neese, F.: Natural triple excitations in local coupled cluster calculations with pair natural orbitals, J. Chem. Phys., 139, 134101, https://doi.org/10.1063/1.4821834, 2013. a
Riplinger, C., Pinski, P., Becker, U., Valeev, E. F., and Neese, F.: Sparse maps – A systematic infrastructure for reduced-scaling electronic structure methods. II. Linear scaling domain based pair natural orbital coupled cluster theory, J. Chem. Phys., 144, 024109, https://doi.org/10.1063/1.4939030, 2016. a
Rozanska, X., Stewart, J. J. P., Ungerer, P., Leblanc, B., Freeman, C., Saxe, P., and Wimmer, E.: High-Throughput Calculations of Molecular Properties in the MedeA Environment: Accuracy of PM7 in Predicting Vibrational Frequencies, Ideal Gas Entropies, Heat Capacities, and Gibbs Free Energies of Organic Molecules, J. Chem. Eng. Data, 59, 3136–3143, https://doi.org/10.1021/je500201y, 2014. a
Schweitzer-Chaput, B., Kurtén, T., and Klussmann, M.: Acid-Mediated Formation of Radicals or Baeyer-Villiger Oxidation from Criegee Adducts, Angew. Chem. Int. Edit., 54, 11848–11851, https://doi.org/10.1002/anie.201505648, 2015. a
Schweitzer-Chaput, B., Kurtén, T., and Klussmann, M.: Corrigendum: Acid-Mediated Formation of Radicals or Baeyer-Villiger Oxidation from Criegee Adducts, Angew. Chem. Int. Edit., 56, 7024–7024, https://doi.org/10.1002/anie.201704761, 2017. a
Stewart, J. J. P.: Optimization of parameters for semiempirical methods VI: More modifications to the NDDO approximations and re-optimization of parameters, J. Mol. Model., 120, 1–32, 2013. a
Thomas, J. M., He, S., Larriba-Andaluz, C., DePalma, J. W., Johnston, M. V., and Hogen, Jr., C. J.: Ion mobility spectrometry-mass spectrometry examination of the structures, stabilities, and extents of hydration of dimethylamine–sulfuric acid clusters, Phys. Chem. Chem. Phys., 18, 22962–22972, https://doi.org/10.1039/c6cp03432b, 2016. a
Zanca, T.: HOM-cluster-decomposition-in-APi-TOF-MS_Data, GitHub, available at: https://github.com/tommyzeta/HOM-cluster-decomposition-in-APi-TOF-MS_Data, last access: 30 June 2020. a
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, https://doi.org/10.1021/acs.jpca.8b10744, 2018. a, b, c, d, e
Zhang, J. and Dolg, M.: ABCluster: the artificial bee colony algorithm for cluster global optimization, Phys. Chem. Chem. Phys., 17, 24173–24181, 2015. a
Zhang, J. and Dolg, M.: Global optimization of rigid molecular clusters by the artificial bee colony algorithm, Phys. Chem. Chem. Phys., 18, 3003–3010, 2016. a