Andriy
Pysanenko
a,
Karolína
Fárníková
a,
Jozef
Lengyel
*b,
Eva
Pluhařová
*a and
Michal
Fárník
*a
aJ. Heyrovský Institute of Physical Chemistry, v.v.i., The Czech Academy of Sciences, Dolejškova 2155/3, Prague 8 182 23, Czech Republic. E-mail: michal.farnik@jh-inst.cas.cz; eva.pluharova@jh-inst.cas.cz
bChair of Physical Chemistry, TUM School of Natural Sciences, Technical University of Munich, Garching 85748, Germany. E-mail: jozef.lengyel@tum.de
First published on 26th August 2022
Mixed nitric acid/water clusters with dimethylamine (DMA) represent a suitable model system for understanding acid–base chemistry in atmospherically relevant clusters. We investigate these clusters in a detailed molecular-beam experiment accompanied by ab initio calculations. The (HNO3)M(H2O)N clusters are produced by supersonic expansion into vacuum and doped by DMA molecules in a pickup process. Two complementary mass spectrometry approaches are employed to analyze the resulting (DMA)K(HNO3)M(H2O)N clusters: (i) electron impact ionization at 70 eV to form positive cluster ions and (ii) low-energy electron attachment at 0–10 eV to form negative clusters. The positive ion spectra contain mainly protonated (DMA)k(HNO3)mH+ clusters with k = m + 1, whereas the negative ones are dominated by (DMA)k(HNO3)mNO−3 with m ≥ k, followed by (DMA)k(HNO3)mNO2− (m > k) ions with low abundances. These observations are rationalized by our calculations, which exhibit the tendency of the mixed clusters to maximize the number DMA·H+ and NO3− ions in the clusters. In the neutral clusters, this is fulfilled for 1:1 ratio of DMA and HNO3, while the positively charged (DMA)k(HNO3)mH+ clusters satisfy this condition for k = m + 1. The protonated clusters always contain the DMA·H+ moiety. For the negatively charged cluster ions, thermochemistry supports the prevailing formation of NO−3 and m ≥ k ion composition. Furthermore, the NO−3-containing cluster ions can form when an electron attaches to the protonated moiety of the DMA·H+⋯NO−3 ion pair in the cluster, which leads to H atom evaporation. From the gas phase HNO3 molecule, where NO−2 is formed exclusively upon an electron attachment, the tendency to form NO−3 increases to hydrated HNO3 clusters, where both NO−2 and NO−3 ions are generated in approximately equal abundances, to the DMA doped clusters, where NO−3 strongly prevails NO−2.
Environmental significanceNew particle formation (NPF) in atmospheric aerosol generation is still insufficiently understood. It is well-established that acid-base chemistry plays a key role in NPF, yet, it is difficult to obtain a truly molecular-level understanding from field and aerosol chamber investigations. Studies of isolated clusters in molecular beams can provide such picture. Nitric acid and dimethylamine (DMA) are ubiquitous throughout the atmosphere and are involved in NPF. We dope hydrated nitric acid clusters with DMA in a molecular beam experiment complemented by theoretical calculations, and provide a molecular-level insight into these species and their chemistry. Discussed relationship between neutral cluster structure and positively and negatively charged cluster ions is relevant for implementation of mass spectrometry in aerosol research. |
Significant fraction of atmospheric aerosol particles originates from molecular clusters formed from single molecules in the gas phase via gas-to-particle conversion known as new particle formation (NPF). They grow to larger sizes and eventually become stable against evaporation. However, the detection of the smallest particles in the nanometer size range directly in the atmosphere is difficult.6 Nevertheless, these sizes, where the transition from the gas phase clusters to stable particles occurs, are important.7 The laboratory alternative to investigate the clusters in the NPF process is offered by various aerosol chamber experiments such as CESAM in Paris8 or CLOUD in CERN9 and others,10 or various aerosol flow-tube experiments, see, e.g., ref. 11–13. A different perspective is offered by experiments with clusters in molecular beams in vacuum. They can provide a bottom-up approach and truly molecular-level insight.
It has been well established that the key components in the NPF process are acids, especially sulfuric acid, which have to be stabilized in the clusters by bases such as ammonia and the simple alkyl amines abundant in the atmosphere.2,14–17 Other species, such as highly oxygenated organic molecules (HOMs) might stabilize the clusters as well,18–21 and water vapour is indispensable in the NPF process too.2,14,15 However, dehydrated clusters are often detected in laboratory experiments,17,22–24 thus the exact role of water still has to be better understood.15
Though nitric acid (HNO3) is mainly considered in polar stratospheric cloud particles,3,25 it is ubiquitous throughout the entire atmosphere and its involvement in the particle growth in the atmosphere has been studied.26–29 Recent work has indicated a rapid NPF and growth by nitric acid and ammonia condensation.9 These species are more abundant in polluted areas and thus nitric acid can be responsible for urban smog. Although, there are no indications that clusters of nitric acid HNO3 with dimethylamine (CH3)2NH (DMA) are directly involved in NPF in atmosphere, they represent an excellent model system for understanding detailed acid–base chemistry in the NPF process. Nanoparticle formation and growth from the nitric acid HNO3 and dimethylamine (CH3)2NH (DMA) was investigated before in a flow tube experiment accompanied by theoretical modelling.13 The authors deduced acid:base ratio close to 1:1 in all their relative humidity conditions, which was different from the previous work on sulfuric acid/base particles, namely H2SO4–DMA.12,30,31 Here we approach similar system, DMA–HNO3–H2O clusters, from a complementary experimental and theoretical perspective.
Molecular beam studies of clusters with atmospheric relevance started from the early cluster experiments in 1970s (e.g., see recent reviews 32,33). The early experiments with the hydrated HNO3/H2O clusters performed in Castleman's group,34 revealed the acid dissociation and H3O+⋯NO−3 ion pair generation in these clusters. These early studies were extended by our experiments with Na-doping of the (HNO3)M(H2O)N clusters,35 and electron attachment to these clusters.36,37 Later on, we have investigated also uptake of various atmospherically relevant molecules by these clusters.33,38,39
Here, we dope the preexisting hydrated nitric acid clusters with the dimethylamine molecules: the small well defined (HNO3)M(H2O)N clusters fly isolated in a molecular beam in vacuum and adsorb individual DMA molecules in a pickup chamber.33 Due to the finite size of the clusters, this method allows us observing acid–base reactions and water evaporation directly in the small (DMA)K(HNO3)M(H2O)N clusters. The mass spectrometry of the positively as well as negatively charged clusters accompanied by ab initio quantum chemical calculations of these clusters offers an unprecedented detailed insight into the reaction and ionization processes in these clusters. This truly molecular-level insight provides a basis for detailed understanding of the NPF processes involving nitric acid and dimethylamine.
The cluster beam was skimmed and passed to a pickup chamber, where dimethylamine (99.8% Sigma-Aldrich) was introduced via a needle valve with a controlled pressure ranging from pc = 2 × 10−4 mbar to 2 × 10−3 mbar. Mass spectra at different DMA pressures were recorded. It ought to be stressed that the individual clusters fly isolated without any interaction in high vacuum, and even at the highest pickup pressure pc = 2 × 10−3 mbar they undergo only several individual collisions with isolated molecules (note that this pickup pressure is still almost 6 orders of magnitude below the atmospheric pressure).
After flying another ≈ 1.5 m through two differentially pumped chambers in ultrahigh vacuum (≤10−7 mbar) the clusters entered the vacuum chamber with a perpendicularly mounted reflectron time-of-flight mass spectrometer (RTOF), which was operated either in the positive or negative ion mode.
In the positive mode, the clusters were ionized by a pulsed electron gun at a frequency of 10 kHz at 70 eV electron energy. A typical electron current was 10 μA. The ionization pulse width was 5 μs. After 0.5 μs delay, to avoid any free electron effects, the ions were extracted by a 4 kV pulse and further accelerated to 8 keV into the time-of-flight region. After ≈ 95 cm flight path, the ions were detected by an MCP detector.
In the negative ion mode, the voltages on the extraction plates were inverted. The clusters were ionized at low electron energies Ee 0–14 eV, and the spectra were recorded in the electron energy steps of 0.2 eV and integrated. The electron pulse was 5 μs long, and after 0.5 μs delay a 3 μs extraction pulse was applied to extract the negative ions. A typical electron current at the present experiment was 0.2 μA at 5 eV. The individual mass spectra recorded at different electron energies allow for the evaluation of the electron energy dependence of the negative ion yield for individual mass peaks.37 The energy spectra were normalized to the electron current. Since the current is very low at electron energies close to zero, below Ee ≈ 1.5 eV the spectra can be subject to a large error.37
To sample the configurational space of the clusters, we used Born–Oppenheimer molecular dynamics40 using the CP2K program package.41 For treatment of the electronic structure we applied revPBE42 functional with the D3 dispersion correction.43 The core electrons were described by the Goedecker–Tetter–Hutter norm-conserving pseudopotentials,44 with Kohn–Sham orbitals being expanded in TZV2P MOLOPT basis set.45 For the auxiliary plane-basis set a cutoff of 400 Ry was used. The equations of motion were integrated with a time step of 0.5 fs. To keep the temperature constant a CSVR thermostat46 was used with a coupling constant of 50 fs. We generated at least one 20 ps trajectory at 300 K for each of the systems, specifically one for clusters containing one HNO3 unit, two for dehydrated clusters with two HNO3 units and three for larger clusters. Periodic boundary conditions were applied. Box length of the cubic simulation cell was 14 Å.
To calculate the most stable structure of each system, we optimized structures at first at the PM6 (ref. 47) level of theory which was then followed up by the M06-2X/aug-cc-pVDZ48,49 level of theory. The DFT functional has been chosen in order to be consistent with previous studies of systems containing hydrated nitric acid36,39 where it was benchmarked against CCSD(T). M06-2X has been shown to provide reliable results for atmospherically relevant clusters.50 We gathered these structures from snapshots of the previously calculated trajectories. Snapshots were spaced by 1 ps. We performed frequency analysis to confirm the local minimum. The energy span of the resulting structures was 20–90 kJ mol−1, the structure with the lowest energy was then used to evaluate the binding energy of either dimethylamine or water to the cluster at M06-2X/aug-cc-pVDZ level. The binding energy includes zero-point vibrational energy. We performed all quantum chemical calculations in the Gaussian 16.51
Panels below show the spectra evolution with increasing pressure of DMA in the pickup chamber (b) 8.6 × 10−4 mbar and (c) 1.1 × 10−3 mbar. Additional mass spectra measured at different pc pressures are shown in ESI.† The spectra were measured under the same conditions except for the pickup gas pressure pc. Thus, the differences between panels in the relative signal intensities reflect the decreasing cluster beam intensity with increasing DMA pressure due to the cluster scattering out of the beam. The relative increase of peak intensities at the higher mass side of the spectrum with respect to the lower mass side is simply the consequence of the easier scattering of smaller (lighter) clusters out of the beam.
The important difference between the spectra are the pronounced peaks that clearly stand out as the DMA pressure increases. They are labeled by their mass-to-charge ratio (m/z) in panel (b). Although there are mass coincidences with the original (HNO3)m(H2O)nH+ masses, these peaks must be due to the picked up DMA molecules. Therefore, they are assigned to (DMA)k(HNO3)mH+ ions with k = m + 1, as indicated in the last panel (c) where these peaks are labeled by k:m. Upon careful inspection, further mass peaks corresponding to the (DMA)k(HNO3)mH+ series with k = m could be unambiguously identified as indicated also in panel (c). However, this series is very weak compared to the k = m + 1 series. No hydrated species (DMA)k(HNO3)m(H2O)nH+ could be identified unambiguously. They would partly coincide with the original (HNO3)m(H2O)nH+ series, however, we have not observed any intensity anomalies which would point to the generation of the (DMA)k(HNO3)m(H2O)nH+ mass peaks upon the DMA uptake.
The negative ion spectrum after the pickup is shown in the bottom panel (Fig. 2b). The change in the spectrum with pickup pressures was relatively abrupt and appeared just at higher pc. Therefore only the spectrum at pc = 8.4 × 10−4 mbar is shown in (Fig. 2b), which corresponds to the positive ion spectrum in (Fig. 1b). Due to the scattering of the clusters from the beam, the ion abundances are very low at higher pc pressures.
Clearly, new major series appear in the spectrum, which could be unambiguously assigned to (DMA)k(HNO3)mNO−3 ions (triangles up). The series distinguished by different colors and connected by dashed lines always start with k = m ion followed by ions with more HNO3 than DMA molecules m > k. Series with k = 1–3 and m up to k + 4 and k + 3 are indicated in (Fig. 2b). New series with DMA molecules and NO−2 ions, (DMA)k(HNO3)mNO−2, could be identified in the spectrum too, however, at very low intensities. The series indicated by the orange diamonds corresponds to the k = 1 case. The first member of this series corresponds to m = k + 1. A very weak k = 2 series could be found upon a closer inspection (not indicated to avoid congestion).
Traces of the series from the top spectrum can be found in (Fig. 2b) too, and are labeled by the same symbols as in the top panel (circles and triangles down). We note that the relative intensities of the original NO−2 and NO−3 series change upon the DMA pickup: while in the pure cluster spectrum (a), the NO−2 containing peaks (circles) are higher than the corresponding NO−3 containing peaks (triangles down), in the bottom spectrum with DMA pickup (b) the ratio is reversed. The total integrated intensities of the NO−3 ion containing series (including the ones with DMA) constitute over 80% of the mass spectrum, i.e., the ratio of the cluster ion series containing NO−3:NO−2 is about 4.3 ± 0.5 (the larger error bars are due to the lower intensity and ambiguity in the assignment of some less abundant peaks due to mass coincidences).
The positive ion mass spectra indicated no clear evidence for any mass peaks containing DMA, HNO3 and water. Since we start with strongly hydrated clusters that pick up DMA molecules, the question arises, if the water evaporates upon the DMA pickup and subsequent acid–base reaction, or if it is the result of the cluster ionization process eventually followed by ion–molecule reactions. Therefore, we seek for evidence of hydrated mixed clusters upon the electron attachment. The attachment might be a softer ionization leading less to water evaporation than the 70 eV electron impact ionization. There is a series of 3 peaks following the first (DMA)k(HNO3)mNO−3, k = m = 1, which are displaced by m/z = 18. They could correspond to the hydrated cluster ions (DMA)k(HNO3)m(H2O)nNO−3, n = 1–3. Due to mass coincidences we discuss these peaks in more detail in ESI.† Nevertheless, we can safely conclude that the abundance of hydrated cluster ions is strongly suppressed with respect to the non-hydrated ones.
For comparison, the grey line shows the Ee-dependence of (HNO3)NO−3, m/z = 125 mass peak in the pure spectrum without pickup (note, that its intensity was scaled down by more that two orders of magnitude). This spectrum is identical with the Ee-dependencies of electron attachment to the pure (HNO3)M(H2O)N clusters measured and analyzed in detail previously.37 Thus, Fig. 3 demonstrates that the energy spectra of the DMA-containing clusters exhibit essentially the same character as the spectra of the pure clusters without the DMA pickup, and the electron attachment proceeds via the same mechanism to the pure (HNO3)M(H2O)N as well as to the (DMA)K(HNO3)M(H2O)N clusters with the adsorbed DMA molecules. Due to the relatively low signal-to-noise ratio of the negative ion signal and mass coincidences, no further unambiguous conclusions could be drawn from the energy dependencies, and therefore we present somewhat extended discussion only in ESI† for interested reader.
The left panel (Fig. 4a) illustrates that binding of a water molecule is about 45 kJ mol−1 nearly independently of the clusters size. Therefore, the pickup of DMA is sufficient to evaporate two water molecules. The pickup of the second DMA molecule is possible, but it releases less energy, because the molecule remains neutral and has a binding strength comparable to water.
DMA is a stronger base than NO−3 or water, thus clusters with n ≥ 1 contain the DMA·H+⋯NO−3 ion pair. This feature is preserved in larger structures (DMA:HNO3 2:2 or 3:3) which possess acidically dissociated structures even without any water (see Fig. 4 bottom).
The experimental mass spectra contain predominantly the series (DMA)k(HNO3)mH+ ions with k = m + 1 and small signals corresponding to k = m. This can be explained by the binding energy of each of the molecules (Fig. 5 top). The most likely molecule to evaporate upon the ionization is water with binding energy of around 55 kJ mol−1 in protonated clusters of different composition. The binding of water to the protonated clusters containing DMA (red and blue circles, full lines) is slightly weaker than to those without DMA (black circles). The overall positive charge of the cluster strengthens water binding by about 10 kJ mol−1 in comparison with the neutral structure. The second weakest bound moiety is nitric acid with about 85 kJ mol−1 cost for evaporation from DMA2(HNO3)2(H2O)3H+ (green line). The most tightly bound species is DMA (about 150 kJ mol−1), because it is protonated in all clusters. The binding energy of the first DMA is by about 25 kJ mol−1 higher than of the second one for n = 3. This difference decreases with increasing number of water molecules. Thus, evaporation of HNO3 from DMA2(HNO3)2(H2O)nH+ giving DMA2(HNO3)(H2O)nH+ is plausible and corroborates with the dominant k = m + 1 series in the mass spectra.
The tight binding of both DMA molecules in the positively charged clusters corresponds to the fact that in order to evaporate them, the extra proton has to be accepted by a weaker base at first. From structural point of view, the clusters maximizing the number of DMA·H+ and NO−3 ions are the optimal ones. This is fulfilled for K = M in the neutral clusters and for k = m + 1 in the protonated ones (Fig. 5 left).
DMAK(HNO3)M(H2O)N + e− → DMAK(HNO3)M−1(H2O)NNO3− + H | (1) |
DMAK(HNO3)M(H2O)N + e− → DMAK(HNO3)M−1(H2O)NNO2− + OH | (2) |
DMAK(HNO3)M(H2O)N + e− → DMAK(HNO3)M−2(H2O)N+1NO3− + NO2 | (3) |
Based on the reaction enthalpies (Fig. 6), the predominant pathway is the release of NO2 radical (∼−130 kJ mol−1) and formation of NO−3-containing ions. The reaction enthalpy does not depend significantly on the amount of DMA, HNO3 and H2O in the clusters. The second, but less likely possibility, is the formation of OH radical and generation of NO−2-containing ions which is slightly exothermic (∼−20 kJ mol−1) in our clusters. The hydrogen radical is unlikely to form, because the reactions are endothermic (∼40 kJ mol−1). Thus, the energetics suggests that the clusters containing NO−3 will be more populated than those with NO−2.
If just the reactions triggered by electron attachment to the DMAK(HNO3)M(H2O)N clusters containing the same number of DMA and HNO3 molecules K = M were considered, the resulting structures should contain more DMA than HNO3. This is not the case, since we observe DMAk(HNO3)mNO−3 cluster ions with m ≥ k. Therefore, evaporation of molecules from the clusters generated immediately after the electron attachment has to be quantified too. The binding strengths of DMA and H2O to the negatively charged clusters are similar to those in neutral clusters. The binding energy of a water molecule is about 50 kJ mol−1 (see ESI).† Binding of DMA in the intermediate clusters is on average by a few kJ mol−1 weaker (Fig. 7 black, blue and magenta lines). Clusters containing NO−3 or NO−2 do not differ in this aspect. It is likely that water molecules and DMA molecules which cannot be protonated evaporate leading to the cluster composition with m = k (see labels of the black, blue and magenta lines).
In contrast to that, HNO3 is tightly bound (at least 150 kJ mol−1) in the DMAk(HNO3)m(H2O)nNO−3 or DMAk(HNO3)m(H2O)nNO−2m = k clusters (Fig. 7 red). DMA binding to the m = k cluster (Fig. 7 green) is also stronger than to the above mentioned m ≤ k clusters, but the interaction is weaker than in case of HNO3. The corresponding values increase from 80 kJ mol−1 for less hydrated clusters containing NO−3 to 120 kJ mol−1 more hydrated clusters. The lower values for analogous process in DMAk(HNO3)m(H2O)nNO−2 can be rationalized by the presence of HONO (not NO−2) (see ESI).† Thus, the evaporation of DMA from the m = k clusters leading to the experimentally observed m ≥ k series is more likely than the hypothetical HNO3 evaporation. Even structures with m ≥ k bind HNO3 quite strongly, even though not as much as in the m = k case. The evaporation costs more than 90 kJ mol−1 (Fig. 7 orange).
Fig. 7 (right) shows some representative structures. The intermediate DMA3(HNO3)NO−3 formed after electron attachment to the neutral DMA3(HNO3)2 clusters has one DMA·H+ symmetrically surrounded by NO−3 and two neutral DMA molecules. The facile evaporation of the later leads to DMA(HNO3)NO−3 which is in fact DMA·H+(NO−3)2. The tendency to maximize the number of protonated DMA·H+ and anionic units is preserved in negatively charged clusters too. It occurs for k = m in DMAk(HNO3)mNO−3,2. In addition, the evaporation of DMA from the m = k clusters is more facile than evaporation of HNO3. The resulting m ≥ k structures observed in the experimental spectrum are still relatively stable with respect to the HNO3 evaporation.
The calculations suggest that the optimal structures are those maximizing the number of DMA·H+ and NO−3 ions. This occurs for K = M in neutral clusters, for k = m + 1 in protonated ones and for k = m in DMAk(HNO3)mNO−3,2. Fig. 7 rationalizes the experimental observation that the major peaks in the negative spectra correspond to m > k (specifically, m = k + 1 and m = k + 2), since additional DMA molecules can evaporate from the negatively charged clusters, while the evaporation of HNO3 is less likely, as discussed above. In addition, we have to explain why the clusters with m ≥ k + 2 are still observed experimentally. The experimental spectrum reflects the stability of the ions as well as the population of the original neutral species. If we start from neutral clusters with more HNO3 than DMA molecules M > K, it is easier to end up with the m > k negatively charged ionic clusters. It ought to be mentioned that the neutral (DMA)K(HNO3)M(H2O)N clusters are generated by the pickup of K molecules of DMA on (HNO3)M(H2O)N, and thus the initial K:M composition can be different from the optimal K = M neutral cluster structure, and the clusters with M > K can be more abundant. Upon the positive ionization, there is enough energy delivered in the ionization process to reach the optimal k = m + 1 structure of the protonated (DMA)k(HNO3)mH+ cluster ions, while the clusters with m ≥ k + 2 are still observed upon the electron attachment.
Second point to discuss is the dehydration of clusters. We start the experiment with highly hydrated clusters before the DMA pickup. In the positive (HNO3)m(H2O)nH+ spectrum of pure clusters without DMA pickup, the average cluster ion contains six to three H2O molecules per one HNO3 molecule. In addition, at least OH and probably also some water molecules are evaporated upon the electron ionization at 70 eV. Therefore, the number of H2O molecules in the neutral cluster is even larger than in the observed cluster ion fragment, N > n. Upon the DMA pickup mostly cluster ions without any water are observed in the mass spectrum. Does the water evaporate from the neutral cluster due to the acid–base reaction, or upon the subsequent ionization? The comparison of the binding energies suggest that the uptake of one DMA and subsequent acid–base reaction is sufficient to evaporate about two water molecules. In addition, the average collision energy of DMA molecules with clusters in our experiments corresponds to about 58 kJ mol−1, which could lead to another single water molecule evaporation. However, this is not quite sufficient for the complete dehydration of the clusters, thus, the remaining water molecules are evaporated by both, positive and negative ionization processes.
Finally, we concentrate on the effect of cluster medium observed in the negative ion spectra. Let us remind that the DEA of an isolated HNO3 molecule yielded essentially only NO−2 ions with the relative abundance of 96.5% with almost no traces of NO−3.36 Going from the gas phase HNO3 to the hydrated HNO3 clusters, our previous investigations36,37 showed that the relative population of NO−3 ions upon the electron attachment substantially increased. In the present experiment, the ratio of integrated intensities of the ion series containing NO−3 and NO−2 is 1.0 ± 0.1. Upon further addition of the DMA molecule, the propensity for the NO−3 ion based clusters is even stronger with the ratio of the cluster ion series containing NO−3:NO−2 of about 4.3 ± 0.5.
The change of the NO−3:NO−2 ratio upon the uptake of DMA can be discussed in terms of a simple mechanism of electron attachment to an ion pair formed in cluster. If a neutral cluster contains an ion pair – either H3O+⋯NO−3 or DMA·H+⋯NO−3 – the electron is likely to attach to the protonated species leading to the hydrogen atom ejection from the cluster and the NO−3-containing cluster ion formation. For HNO3(H2O)N clusters with single HNO3 molecule, the H3O+⋯NO−3 ion pair formation occurs at n ≥ 4.34–36,52,53 Our present calculations indicate that the DMA·H+⋯NO−3 ion pair is already formed in clusters with a single water molecule because of DMA higher basicity. With addition of second HNO3 molecule, the number of water molecules required to form an ion pair is reduced even further. For example, with two HNO3 molecules only three water molecules are needed for the ion pair to be energetically preferred in (HNO3)2(H2O)3 cluster, and the presence of DMA in the cluster induces spontaneous proton transfer even without any water. Thus, the relative abundance of ion-pair containing clusters increases upon the DMA addition. Since they lead to the NO−3-containing cluster ions upon the electron attachment, the relative abundance of these ions increases upon the DMA addition, as observed experimentally.
However, this simple mechanism cannot explain the observation of the NO−2 containing cluster ions, which are present in the spectra too, albeit at low intensities. Therefore, we present the argumentation based on the reaction energetics summarized in Fig. 6: the formation of NO−3 is more exothermic than NO−2, but NO−2-containing clusters can be formed as well.
• Optimal neutral cluster composition corresponds to DMA and HNO3 in 1:1 ratio. This is rationalized by the tendency to maximize the number DMA·H+ and NO−3 ions in the clusters and by binding energies of cluster constituents.
• The dominant series in the positive mass spectra is k = m + 1 which is nicely supported by calculated structures where the additional H+ resides on the extra DMA and prevents its evaporation.
• Upon an electron attachment, a formation of NO−3 is more likely than NO−2 which is due to the presence of the DMA·H+⋯NO−3 ion pair structures and reaction energetics (shown in Fig. 6). In addition, calculations give k = m as the optimal composition. The observation m > k is due to the fact that the negatively charged cluster evaporate DMA more easily than HNO3 according to the calculated evaporation enthalpies, and also due to the experimental neutral cluster composition.
An interesting point is the effect of cluster environment on the product ions observed in the electron attachment experiments. The population of NO−3 and NO−2 ions depends strongly on the medium. The dissociative electron attachment (DEA) to an isolated HNO3 molecule yields essentially exclusively NO−2 ions, while the NO−3 ion abundance increases with hydration: in the present pure HNO3/H2O clusters the ratio of integrated intensities of the ion series containing NO−3 and NO−2 ions is 1.0 ± 0.1. With addition of the strong base DMA, this ratio increased to about 4.3 ± 0.5. The propensity for NO−3-containing cluster ions can be explained by a simple model of electron attachment to the protonated species in the neutral clusters containing H3O+⋯NO−3 or DMA·H+⋯NO−3 ion pairs, which are generated due to the proton transfer in the acid–base reaction. Nevertheless, to justify the generation of both NO−3 and NO−2-containing cluster ions, we calculate the reaction energetics, which explains the observed tendencies qualitatively.
We also demonstrate that the observed ionic species strongly depend on the ionization procedure. While the positive ion spectrum suggests prevalence of DMA molecules in the clusters, the negative spectrum clearly points to more HNO3 molecules in clusters. In addition, the observed major cluster ions do not contain any water, while the spectra of pure clusters without DMA pickup clearly exhibit a high level of hydration in the neutral clusters. Therefore, extreme caution has to be exercised when conclusions about neutral cluster composition are made from mass spectra, which should be considered in any implementation of mass spectrometry in aerosol research.
Footnote |
† Electronic supplementary information (ESI) available: Experimental details, additional mass spectra and their discussion, further computational results, cluster structures and coordinates. See https://doi.org/10.1039/d2ea00094f |
This journal is © The Royal Society of Chemistry 2022 |