Original Research ARTICLE
Association and Proton Transfer Reactions With a Series of Organic Molecules
- 1Institute of Ion Physics and Applied Physics, University of Innsbruck, Innsbruck, Austria
- 2Department of Chemistry and Institute for Atmospheric and Earth System Research (INAR), University of Helsinki, Helsinki, Finland
In this study, we present reactions of with a series of analytes (A): acetone (C3H6O), methyl vinyl ketone (C4H6O), methyl ethyl ketone (C4H8O), and eight monoterpene isomers (C10H16) using a Selective Reagent Ionization Time-of-Flight Mass Spectrometer (SRI-ToF-MS). We studied the ion-molecule reactions at collision energies of 55 and 80 meV. The ketones, having a substantially lower proton affinity than NH3, produce only cluster ions (A) in detectable amounts at 55 meV. At 80 meV, no cluster ions were detected meaning that these adduct ions are formed by strongly temperature dependent association reactions. Bond energies of cluster ions and proton affinities for most monoterpenes are not known and were estimated by high level quantum chemical calculations. The calculations reveal monoterpene proton affinities, which range from slightly smaller to substantially higher than the proton affinity of NH3. Proton affinities and cluster bond energies allow to group the monoterpenes as a function of the enthalpy for the dissociation reaction . We find that this enthalpy can be used to predict the (A) cluster ion yield. The present study explains product ion formation involving ion chemistry. This is of importance for chemical ionization mass spectrometry (CIMS) utilizing as well as (H2O) as reagent ions to quantitatively detect atmospherically important organic compounds in real-time.
In the 1990's proton transfer reaction mass spectrometry (Hansel et al., 1995; Lindinger et al., 1998b) using H3O+ reagent ions became a widely used analytical instrument with applications in environmental science, medical applications, and food technology due to the large amount of volatile organic molecules, which can be quantitatively ionized. H3O+ undergoes proton transfer reactions with every analyte having a higher proton affinity (PA) than water [PA(H2O) = 165.0 kcal/mol (Hunter and Lias, 1998)]. In contrast, ionization is more specific. Due to the higher proton affinity of ammonia PA(NH3) = 204.0 kcal/mol (Hunter and Lias, 1998), exothermic, thus fast, proton transfer reactions between and analyte are limited to a much smaller number of molecules. If the analyte (A) possesses a proton affinity sufficiently larger than NH3, then reaction (1):
is exothermic and will occur on every collisions, which means that the reaction rate is close to the collisional limit value (Lindinger et al., 1998b). Lindinger et al. (1998a) used proton transfer reactions of to separate the isomeric molecules α-pinene [PA = 204–209 kcal/mol (Lindinger et al., 1998a; Solouki and Szulejko, 2007)] and 2-ethyl-3,5-dimethylpyrazine [PA > 204 kcal/mol (Lindinger et al., 1998a)], both having a molecular mass of 136 Th, and found that α-pinene was not ionized by . Keough and Destefano (1981) discussed several factors affecting reactivity in ammonia chemical ionization (CI) mass spectrometry. At that time analyte molecules were introduced directly into the ion source where ammonia is present in large excess. The presence of large amounts of NH3 and analyte at a typical ion source pressure of one Torr complicates the interpretation of product ion formation due to secondary reactions. Since the introduction of PTR-MS in the early 1990's, a strict separation of the ion source from the reaction region (drift tube) was achieved. This is one key factor why chemical ionization using the PTR-MS design became a quantitative analytical instrument. Keough and Destefano (1981) investigated a series of organic compounds with - chemical ionization. They concluded that analytes having a PA < 188 kcal/mol do not yield useful intensities of (A) adduct ions. Very recently, ammonia chemical ionization was found to be an extremely sensitive method detecting quantitatively first generation oxidized molecules as well as highly oxidized organic molecules with adduct ion chemistry (Berndt et al., 2018a,b; Hansel et al., 2018). Zhou et al. (2018) applied chemical ionization in an atmospheric pressure chemical ionization tandem mass spectrometer and investigated the ionization mechanism of molecules with a hydroperoxide moiety.
Here we present detailed results on the mechanism of chemical ionization of eight monoterpenes (C10H16) by chemical ionization. Additionally, we investigated the reaction of acetone (C3H6O), methyl vinyl ketone (MVK, C4H6O) and methyl ethyl ketone (MEK, C4H8O) in reactions with . We have chosen the following atmospherically most relevant monoterpenes according to Sindelarova et al. (2014) and Smolander et al. (2014): α-pinene, β-pinene, limonene, ocimene, myrcene, sabinene, 3-carene, and camphene. Ion chemistry was performed at twice the thermal energy (KEcm = 55 meV) and at a somewhat elevated collision energy (KEcm = 80 meV). Additionally, we investigated the effect of absolute humidity on the outcome of the individual reactions. To confirm the experimental results, quantum mechanical calculations on the proton affinities, adduct ion geometries, cluster-bond energies, and reaction enthalpies were performed.
Materials and Methods
Experimental Setup and Measurement Procedure
The experimental setup is illustrated in Figure 1. For the investigations of the monoterpenes, a temperature stabilized diffusion source was built. Humidified synthetic air was used to dilute the calibration compound diffusing from the diffusion tube. Humidification of the synthetic air was achieved with a Liquid Calibration Unit (LCU, Ionicon Analytik), which also controlled the carrier gas flow set to 3 slm and the variable liquid water flow. For each calibration compound we measured seven different absolute humidities ranging from 4 ± 3 to 25 ± 3 ppth (parts per thousand).
Figure 1. Illustration of the experimental setup. Calibration compounds were either taken from the calibration gas standards or from the diffusion source, which was flushed with a constant synthetic air flow controlled by a mass flow controller (MFC). Different humidity steps were adjusted with the Liquid Calibration Unit (LCU) and monitored by an infrared gas analyzer (IRGA). The humidified synthetic air containing the respective calibration compound(s) was introduced to a selective reagent ion time of flight mass spectrometer (SRI-ToF-MS).
Liquid Calibration Unit (LCU)
The Liquid Calibration Unit (LCU, Ionicon Analytik, Austria) was used to quantitatively evaporate certain amounts of water into the synthetic air stream resulting in absolute humidities in the range of 3–30 ppth. The calibration of the ketones was performed by dynamic dilution of calibration gas standards (Apel Riemer Environmental Inc., Broomfield (CO), USA) in a humidified carrier gas generated by the LCU. Highly water-soluble compounds can be calibrated precisely with the LCU. As monoterpenes are non-polar compounds they are not quite soluble in water. Therefore, we decided to build a temperature-controlled diffusion device, which was combined with the humidified synthetic air stream from the LCU to generate known amounts of monoterpenes in the parts per billion range.
Temperature Stabilized Diffusion Tube
Each diffusion tube containing several milliliters of the respective liquid calibration compound consisted of a 1/8 inch PFA tubing plug (Parker-Hannifin Corporation, Tucson, USA) connected to a PEEK capillary (Vici Valco, Switzerland) of defined length and inner diameter. The capillary was connected through a gas tight 1/8–1/4 inch tee reducer with the carrier gas stream from the LCU. As described in Fuller et al. (1966), the diffusion rate of a diffusion source is strongly temperature dependent. For this purpose, the diffusion tube was placed in a water bath filled with 50 ml water and wrapped with a heating wire and isolation foam. The water temperature was kept constant at 303 ± 3 K and controlled by a temperature controller (Cal3300, CAL Controls Ltd., Libertyville, USA) connected to a type K thermocouple. The heat capacity of the thermally isolated water bath helped to achieve constant volume mixing ratios (VMR) of individual calibration compounds in the air stream. VMR calculations of the individual analyte in the carrier gas have been performed according to McKelvey and Hoelscher (1957). Assuming that saturation of the head space in the diffusion source has occurred, the diffusion rate r [g s−1] of the analyte can be calculated according to Equation (2):
where D [cm2 s−1] is the diffusion coefficient, M [g mol−1] the molecular weight of the analyte, P [Pa] the total pressure, A [cm2] the cross-sectional area of the capillary, L [cm] the length of the capillary, R = 8.314·106 cm3 Pa mol−1 K−1 the gas constant, and Ps [Pa] the saturation vapor pressure of the liquid analyte. Saturation vapor pressures Ps were calculated with MPBWIN v1.43 (© 2000 U.S. Environmental Protection Agency) for the temperature of the diffusion source. Gaussian error analysis estimates an accuracy of the diffusion rate r of ±1% and an accuracy of the obtained volume mixing ratios of ±10–15%. An overview of all values and estimated errors is given in Supplementary Table 2.
The Selective Reagent Ion Time-of-Flight Mass Spectrometer (SRI-ToF-MS) used in this study is based on the PTR-ToF-MS described by Graus et al. (2010) and is adapted for the use with different reagent ions. In principle, it inhibits several advantageous features. The ion source is flushed with ~100 sccm helium, producing only He+ and metastable He*. The chemical ionization gas NH3 (Linde AG, Pullach, Germany) is added later to the source where reactions with He+ and He* lead to the formation of reagent ions. Compared to a standard PTR-ToF-MS, the SRI-ToF-MS is equipped with an ion-funnel located between the ion source and the drift tube having a length of 9 cm. To prevent photochemical reactions in the drift tube, the shape of the ion funnel is constructed in such a way that photons created in the glow discharge ion source don't reach the drift tube. Moreover, the direction of the gas flow through the drift tube and through the ion funnel is opposite to the ion drift direction. This pumping architecture prevents He and CI gases, as well as radicals created in the discharge, from entering the drift region. The metal drift rings used in common PTR-MS instruments are replaced by conductive PEEK rings with a thickness of 6 mm and an inner diameter of 12 mm allowing the measurement of compounds, which catalytically react on metal surfaces. The drift rings are separated by Teflon© spacers of 6 mm thickness. In the present study, the SRI-ToF-MS was operated at 2.3 mbar drift pressure and 35°C drift temperature. The drift voltage was varied between 250 V (E = 27.8 V cm−1) and 400 V (E = 44.4 V cm−1) resulting in an E/N value of 51 and 81 Td, respectively. E is the electric field strength and N the gas number density (1 Td equals 10−17 V cm2). Data processing was performed with an adapted version of the data processing routine described in Breitenlechner et al. (2017) and further data analysis was done with Matlab2018a©. Subsequently, to compensate variations in the reagent ion signal and the mass-depended ion transmission in TOF mass spectrometers, product ion signals (e.g., compound i being detected at mass mi) measured in counts per second (cps) were duty-cycle corrected (dcps; ) and normalized to 106 cps of (normalized counts per second, ncps). To study collision induced dissociation (CID) of adduct ions, we ramped the extraction voltage applied to the lenses in the ion transfer region located between the drift tube and the mass spectrometer. For an E/N of 51 Td in the drift tube, we ramped the extraction voltage in 5 V steps between 15 and 25 V. For an E/N of 81 Td the extraction voltages were changed between 20 and 30 V. The errors of the averaged ion signals in each voltage and humidity step are the standard deviations.
Calculation of Reaction Thermodynamics
In the drift tube, the ions travel as a result of the applied electric field strength E with an increased drift velocity vd through the buffer gas (Lindinger et al., 1998a):
We used ion mobility values μ0 of in N2 from Abedi et al. (2014) and adapted them to the gas number density N in our instrument. N0 is the gas number density at standard temperature and pressure. Thus, the drift time t of the ions traveling through the drift tube of length l can be calculated by:
The reduction of reagent ions  in reactions with the analyte in the drift region is small, thus the reaction can be treated as a pseudo-first order reaction. The density of protonated analytes [AH+] is then given according to Lindinger et al. (1998a) by:
Where [A] is the density of analyte A, k the reaction rate coefficient and t the drift time. The ratio is proportional to the detected ion signal ratio . Commonly the sensitivity (ε) is defined as the detected analyte signal at a volume mixing ratio of 1 ppbv (parts per billion per volume, 1 ppbv = 10−9) normalized to a reagent ion signal i() of 106 cps (Lindinger et al., 1998b). Combining Equations (4) and (5) and using that the volume mixing ratio (VMR) of A is related to [A] by VMR(A)·N = [A], the theoretically maximum sensitivity εcalc can be calculated as:
We compare the calculated sensitivity with the experimentally observed sensitivity εmeas. By dynamically diluting either the calibration gas standard or the diffused analyte from the diffusion source, we obtain a known volume mixing ratio VMR(A) of the analyte A. To determine the sensitivity of an analyte A, the contributions of all product ions of the reaction with must be considered (Cappellin et al., 2012). The measured sensitivity (εmeas) of a substance A is then given by the slope of a linear fit through the scatter plot of the normalized and duty cycle corrected ion signals of all product ions vs. the volume mixing ratio. The efficiency (eff) of a reaction is then given by
The center-of-mass kinetic energy KEcm is calculated according to Lindinger et al. (1998a):
With KEion being the mean kinetic energy of ion drifting in the buffer gas:
where mb is the mass of the buffer gas, Mion the mass of the ions and kb the Boltzmann constant.
Collisional limiting rate coefficients (kc) of ion molecule reactions are calculated according to Su and Chesnavich (1982) and Su (1994) using dipole moments and polarizabilities of the respective analyte (see Supplementary Table 1). Table 1 gives an overview of all reactions. It is worth to mention that we don't have to consider the back reaction of reaction (1) even if the proton affinity of A is only a few kcal/mol higher than the one of NH3 in analogy to the H3O+–formaldehyde reaction system (Hansel et al., 1997). In our case the efficiency of the back reaction is negligible, as in the SRI-ToF the CI gas NH3 does not enter the drift region.
Table 1. Overview of collisional rate coefficients kc, collision energies KEcm, calculated sensitivities εcalc, measured sensitivities εmeas and reaction efficiencies (eff) at dry and humid conditions.
Quantum Chemical Calculations
Proton affinities are not available for most monoterpenes, thus quantum chemical calculations were performed. Known proton affinities for ammonia, acetone, methyl vinyl ketone and methyl ethyl ketone were compared with our calculations. Additionally, the change of standard enthalpies for the ion-molecule reactions, possible protonation sites and probable adduct ion structures were calculated. The conformers were sampled using the systematic conformer search algorithm and the MMFF94 force field on Spartan'16 (Wavefunction, 2016). All conformers were optimized at the B3LYP/6-31+G* level of theory using the Gaussian 09 program (Frisch et al., 2009). The geometries of all conformers within 2 kcal/mol of the lowest-energy conformer were optimized and the harmonic frequencies were calculated at the ωB97X-D/aug-cc-pVTZ level of theory using the ultrafine integration grid. The B3LYP/6-31+G* geometry optimization was omitted for some of the protonation products to avoid the breaking and forming of C-C bonds in the cation product. We found the lowest enthalpy protonation products of the monoterpenes by placing the proton to each of the double bond carbons of the compounds at a time. Generally, the most stable protonation product is found when the positively charged carbon of the protonated structure is tertiary. For 3-carene, we found an energetically more favorable structure where the proton was not added to a double bond carbon. Final single-point energies were calculated at the CCSD(T)-F12/VDZ-F12 level of theory for the lowest enthalpy conformers using the MOLPRO program version 2015.1 (Werner et al., 2012, 2015).
(R)-(-)-limonene (analytical standard), ocimene (mixture of isomers, >90%), camphene (>95%), sabinene (75%), (+)-α-pinene (>99%), and myrcene (analytical standard) were obtained from Sigma Aldrich (Vienna, Austria). (+)-3-carene (>98.5%) and (-)-β-pinene (>99.0%) were purchased from Fluka. Pressured synthetic air grade 5.8 was obtained from Messer (Gumpoldskirchen, Austria), the calibration standard gases were fabricated by Apel Riemer Environmental Inc. (Broomfield, United States). Bottled NH3 grade 3.8 was purchased from Linde AG (Pullach, Germany).
Results and Discussion
Characteristics of the Diffusion Source
The monoterpenes β-pinene and limonene are also present in our calibration gas standards from Apel Riemer Environmental Inc., Broomfield (CO), USA. They certify an accuracy of typically ±10%. We compared the gas calibration results with our home build diffusion source. The agreement between the estimated sensitivities of the diffusion source and the gas standard differed not more than 25% for these two compounds. As many physical and chemical properties of the investigated monoterpenes are not experimentally determined, we had to rely on calculated values for saturation vapor pressures Ps and diffusion coefficients D. Thus, the error of the diffusion source seems very reasonable. To understand the principles of the investigated ion-molecule reactions it is of greater importance that the diffusion rate of the analyte remains constant over the entire measurement period. In our experiment we detected volume mixing ratio drifts of the diffusion source in the range of ± 3% only. Overall we estimate a calibration error of less than ± 30% taking into account also dilution errors from calibrated flow controllers.
Reagent Ion Distribution
Typical reagent ion distributions are shown in Figure 2 as a function of humidity and CID voltage settings. The reagent ion distribution is dominated by . At all voltage settings and humidity steps the signal dominates with typically 106 dcps. The next prominent reagent ion is the hydrated ammonium ion (H2O) (m/z = 36.04 Th), which is typically two orders of magnitude lower in intensity even at humid conditions and low E/N settings. The most abundant other “impurity” ions are (NH3) (m/z = 35.04 Th), H3O+ (m/z = 19.02 Th), NO+ (m/z = 29.99 Th) and (m/z = 31.99 Th). Negligible amounts of H3O+(H2O)2 cluster ions are also present. The Gibb's free energies ΔG to dissociate the respective cluster ions at 310 K drift tube temperature are: (H2O) ΔG = 13.0 kcal/mol (Meot-Ner and Speller, 1986), H3O+H2O(H2O) ΔG = 13.2 kcal/mol (Kebarle et al., 1967), (NH3) ΔG = 17.9 kcal/mol (Payzant et al., 1973). In the drift tube we increased the ion energy above thermal applying E/N values of 51 and 81 Td, respectively. Thus, ΔG dissociation energies of the cluster ions become even smaller for higher E/N settings. This is one reason why cluster ion intensities at 81 Td are substantially lower than at 51 Td. The other reason is that the formation of cluster ions is suppressed as a function of ion collision energy in the drift tube (Hansel et al., 1997; Lindinger et al., 1998a). With increasing extraction voltages, the amount of weakly bound cluster ions decrease due to collision induced dissociation (CID) in the ion transfer region. Overall we can conclude that the prevailing reagent ion reacting with the analytes in the drift tube is .
Figure 2. Typical reagent ion distributions are shown at dry (left) and humid (right) conditions as a function of extraction voltage settings at an E/N value of 51 (top) and 81 Td (bottom), respectively. (H2O) (m/z 36.04), (NH3) (m/z 35.04), H3O+(H2O)2 (m/z 55.04), and H3O+ (m/z 19.02) account for <1% of the (m/z 18.04) signal.
Reactions of With Small Ketones
First, we investigated the reaction of with acetone, methyl vinyl ketone (MVK) and methyl ethyl ketone (MEK). In Figure 3, measured sensitivities of acetone, MVK and MEK at dry and humid conditions and at two E/N values, 51 and 81 Td, are shown. As illustrated in Figure 3, all ketones are detected as (A) cluster ions: acetone: m/z = 76.08 Th, MVK: m/z = 88.08 Th, MEK: m/z = 90.09 Th. At an elevated E/N value of 81 Td, no cluster ions were observed. This is in agreement with first studies of ammonia chemical ionization of ketones with mass analyzed ion kinetic energy (MIKE) spectrometry (Maquestiaut et al., 1980). No protonated ketones have been observed. For all three ketones, proton transfer reactions with are energetically unfavorable due to their lower proton affinities compared to ammonia of PA(NH3) = 204 kcal/mol [PA(acetone) = 194 kcal/mol, PA(MVK) = 199.5 kcal/mol, PA(MEK) = 197.7 kcal/mol (Hunter and Lias, 1998)]. In SRI-ToF-MS, (A) cluster formation proceeds prevailingly as ternary association reactions under dry conditions. The standard reaction mechanism for ternary (three body or collisionally stabilized) ion-neutral reactions proceeds as follows:
Figure 3. Sensitivities (εmeas) of three ketones (acetone: m/z 76.08, MKV: m/z 88.08, MEK: m/z 90.09) for dry (left) and humid (right) conditions as a function of extraction voltages. Cluster ion formation (A), which is the only energetically possible reaction, has been observed at low E/N of 51 Td, only.
where M is a third body, and k+ and k− the reactions rates in the respective direction (Ikezoe et al., 1987). Fast ligand switching reactions of type (11), where the water ligand is exchanged by a more strongly bound ketone, could in principle also produce (A) cluster ions but don't contribute much at 2.3 mbar and at elevated collision energies under dry conditions.
Figure 2 shows that (H2O) reagent ions are less abundant <1% even at humid conditions and at low E/N values.
Table 1 gives an overview of the calculated collisional rate coefficient (kc) at the respective collision energy (KEcm), the calculated sensitivity (εcalc) using kc, the measured sensitivity (εmeas) and the reaction efficiency (eff) for all compounds studied under dry and humid conditions, and at 51 and 81 Td, respectively. The ketones show reaction efficiencies ranging from 4.2% (acetone), and 5.9% (MVK) to 9.4% (MEK) at dry conditions and a KEcm of 0.055 eV. These efficiencies indicate a rather high effective binary rate coefficient, which is most likely due to the long lifetime of the intermediate (NH4A)+* against unimolecular decomposition (i.e., k− is small relative to the stabilization rate in collisions with M) for these polyatomic intermediates, which vary with complexity (Johnston, 1966). It is also indicative that the -ketone cluster is strongly bound, as confirmed by quantum chemical calculations, see Table 2 and (Frege et al., 2018). In the lowest energy geometry of the -ketone cluster, one of the hydrogen atoms of is hydrogen bonded to the oxygen of the carbonyl group (see Figure 4) resulting in typical bond energies of ~ 26 kcal/mol. Adams et al. (2003) studied the association reactions of with a series of organic molecules using a Selected Ion Flow Tube (SIFT). In this study, Adams et al. (2003) reported that the reaction efficiency for acetone is 22% at 300 K and at 0.5 Torr He pressure. When the efficiency is high as for acetone, this indicates that the ternary reaction is close to pressure saturation, i.e., independent of the He pressure, and in this case only a lower limit to the ternary rate coefficient can be obtained. Earlier studies have shown that the ternary rate coefficient k decreases dramatically as a function of temperature (collision energy) with k α T−n, where n ranges from 2 to 3 (Adams and Smith, 1981). An E/N of 51 Td corresponds to 0.055 eV (≈ 600 K). According to Breitenlechner (2011) the fluctuation of the electric field strength along the central axis of the drift tube of the SRI-ToF-MS is within ±10%. Using k α T−n with n = 2–3 reduces the reaction efficiency for the -acetone cluster ion formation at 51 Td to 22% × 2−(2to3) = 5.5–2.8% in excellent agreement with our measured reaction efficiency of 4.2%. At an E/N of 81 Td we could not observe any -ketone cluster ions reliably. This is most likely due to the short lifetime of the more excited (NH4A)+* intermediates at this enhanced collision energy. Under humid conditions the reaction efficiencies of all ketones are increased by about 5% at 51 Td. This increase is somewhat unexpected and might be due to a larger amount of (H2O) cluster ions in the drift tube than measured even at the lowest extraction voltage setting of 15 V. Figure 2 shows that about 1% of all reagent ions comprise the hydrated ammonium cluster ions. We would need about 5% of hydrated ammonium cluster ions, which undergo exothermic thus fast ligand switching reactions according reaction (11). Reaction (11) is exothermic for all three ketones. The bond energy of (H2O) is 20.6 kcal/mol (Meot-Ner and Speller, 1986) (see Table 1), which is smaller compared to the -ketone bond energies ranging from 25.9 kcal/mol (MEK), 26.4 kcal/mol (acetone) to 27.3 kcal/mol (MVK) (see Table 2). Increasing the extraction voltage from 15 to 25 V decreases the -ketone adduct ions due to CID in the extraction region by 60–70% (Figure 3). In comparison, the (H2O) reagent ions show an even more pronounced decrease as a function of extraction voltage (Figure 2). At 25 V, more than 95% of hydrated ammonium cluster ions are lost compared to 15 V extraction voltage. This demonstrates again that cluster ions with a lower bond energy are lost more efficiently in the ion transfer region, supporting our assumption that (H2O) cluster ions could be lost even at an extraction voltage of 15 V. In any case our measured amount of (H2O) cluster ions (1% compared to ) is a lower limit and the enhanced ketone reactivity at humid conditions indicate a higher amount of 5%.
Table 2. Overview of the bond energies (BE), reaction enthalpies (ΔHr), proton affinities (PA), and protonated structures, calculated at the CCSD(T)-F12/VDZ-F12//ωB97X-D/aug-cc-pVTZ level of theory at 298 K.
Figure 4. Geometries of the -A cluster ions for all investigated analytes. Color coding: gray = C, white = H, red = O, blue = N. The green bonds indicate C = C double bonds in the analytes.
Reactions of With Monoterpenes
Product Ion Formation
We investigated the reaction of with eight atmospherically most common monoterpenes at dry and humid conditions, and at two E/N values 51 and 81 Td. In contrast to the measured ketones, for which proton transfer reactions are energetically unfavorable, the proton affinities of the monoterpenes range from 201.2 kcal/mol (limonene) to 215.9 kcal/mol (sabinene), according to our quantum chemical calculations (Table 2). Experimental proton affinities of monoterpenes are rare and exist only for limonene (PA = 209.1 ± 1.2 kcal/mol; Tereza Fernandez et al., 1998), camphene (PA = 205.7 ± 3.2 kcal/mol; Solouki and Szulejko, 2007). For α-pinene (Lindinger et al., 1998a) suggested an upper limit PA < 204 kcal/mol whereas (Solouki and Szulejko, 2007) estimate ~209 kcal/mol. The few experimentally measured proton affinities are in reasonable agreement with our calculations. A list of calculated proton affinities is given in Table 2. Our calculated proton affinity of NH3 is 203.8 kcal/mol, which is in excellent agreement with experimental values of 204 kcal/mol (Hunter and Lias, 1998). All investigated monoterpenes, except limonene, have higher proton affinities than ammonia making the proton transfer reaction (1) exothermic. Our results show significant differences in the product ion distribution of the eight studied monoterpenes. As expected we identified the protonated terpene ion C10H16-H+ (m/z = 137.13 Th), but also fragment ion C6 (m/z = 81.07 Th) and for some monoterpenes additionally small amounts of fragments C7 (m/z = 95.08 Th) and C7 (m/z = 93.07 Th). Here we record only ion signals, which are detected with relative intensities >1%. Other fragments reported in the literature for H3O+ chemical ionization using PTR-MS instruments (Wang et al., 2003; Schoon et al., 2004; Tani et al., 2004; Materić et al., 2017) have not been observed. This could be explained by the smaller amount of transferred internal energy using instead of H3O+ as reagent ion. The difference in proton affinities between the precursor ion and the respective monoterpene is transferred to the product ion causing fragmentation. Besides proton transfer product ions and corresponding fragment ions, we observe also cluster ions attached to monoterpenes for all eight monoterpenes. Product ion distributions for the eight monoterpenes are shown in Figures 5–9 and Supplementary Figures 1–3 at dry and humid conditions, at two E/N values 51 and 81 Td and as a function of extraction voltages.
Figure 5. Camphene product ion distributions are shown at dry (3 ± 1 ppth; left) and humid (18 ± 1 ppth; right) conditions as a function of extraction voltage settings at an E/N value of 51 Td (top) and 81 Td (bottom), respectively.
First, we will discuss the results for dry conditions at an E/N of 51 Td and an extraction voltage of 20 V only. Camphene (Figure 5) and sabinene (Supplementary Figure 1) are the monoterpenes producing only 5% cluster ions -C10H16 (m/z 154.16). The main product ions are protonated monoterpenes C10H16-H+ (m/z = 137.13 Th) and the corresponding fragment ion C6 (m/z = 81.07 Th). In the case of sabinene (PA = 215.9 kcal/mol) 85% of product ions are protonated sabinene and 10% are found as C6 (m/z = 81.07 Th) fragment, while camphene (PA = 207.2 kcal/mol) produces 93% protonated camphene and only 2% C6 fragment ions. Sabinene, having a higher proton affinity compared to camphene, shows a higher amount of fragment ions. Hence more internal energy is generated in the proton transfer channel of sabinene explaining the higher amount of fragmentation. Ocimene (PA = 210.7 kcal/mol) (Supplementary Figure 2) and the two bicyclic monoterpenes β-pinene (PA = 208.7 kcal/mol) (Supplementary Figure 3) and α-pinene (PA = 206.3 kcal/mol) (Figure 6) produce 30–40% cluster ions -C10H16. The remaining fraction is found prevailingly as protonated monoterpene C10H16-H+ and to a lesser amount at the C6 fragment ion. For the acyclic monoterpene myrcene (PA = 204.2 kcal/mol) (Figure 7) and bicyclic 3-carene (PA = 205.4 kcal/mol) (Figure 8) we detect approximately 50% of the product ion signal at the protonated mass. 3-carene exhibited a slightly dominant cluster ion yield (~54%). Limonene (PA = 201.2 kcal/mol) (Figure 9) shows the highest yield of cluster ions -C10H16, namely ~85% and the rest is the unfragmented C10H16-H+ ion.
Figure 6. α-pinene product ion distributions are shown at dry (7 ± 1 ppth; left) and humid (26 ± 1 ppth; right) conditions as a function of extraction voltage settings at an E/N value of 51 Td (top) and 81 Td (bottom), respectively.
Figure 7. Myrcene product ion distributions are shown at dry (7 ± 1 ppth; left) and humid (26 ± 1 ppth; right) conditions as a function of extraction voltage settings at an E/N value of 51 Td (top) and 81 Td (bottom), respectively.
Figure 8. 3-Carene product ion distributions are shown at dry (5 ± 1 ppth; left) and humid (26 ± 1 ppth; right) conditions as a function of extraction voltage settings at an E/N value of 51 Td (top) and 81 Td (bottom), respectively.
Figure 9. Limonene product ion distributions are shown at dry (7 ± 1 ppth; left) and humid (26 ± 1 ppth; right) conditions as a function of extraction voltage settings at an E/N value of 51 Td (top) and 81 Td (bottom), respectively.
We performed detailed quantum chemical calculations to better understand the reaction mechanism with monoterpenes. Table 2 gives an overview of calculated proton affinities compared to literature values, the bond energies (BE, in enthalpy) of cluster ions -A, and reaction enthalpies ΔHr of reaction (12).
The structures of protonated compounds are also shown in Table 2. In the cases of 3-carene and limonene we calculated two structures with similar proton affinities, respectively. The lowest enthalpy geometries of all (A) cluster ions are shown in Figure 4. The calculations of the lowest enthalpy conformers of -monoterpene cluster ions reveal that at least one hydrogen of the forms a hydrogen bond to the C = C double bond with a typical bond energy of ~19 kcal/mol. Some monoterpenes have more than one C = C double bond offering the possibility to form a second hydrogen bond. This is the case for myrcene (BE = 20.9 kcal/mol), limonene (BE = 22.3 kcal/mol), and ocimene (BE = 26 kcal/mol). Bond energies of these monoterpenes are only slightly higher (3–7 kcal/mol) than singly bonding monoterpenes. This is in contrast to calculated bond energies of compounds containing several carbonyl groups. Introducing a second C = O group increases the stability of the cluster ion considerably (almost 2-fold). Additionally, the position of the second functional group to form an optimal hydrogen bond (with a 180° angle of N-H-O) strongly influences the stability of -carbonyl adduct ions (Frege et al., 2018). To predict the yield of the cluster ion formation channel for the monoterpene reactions we correlated the fraction of measured cluster ions as a function of monoterpene proton affinities resulting in a correlation coefficient of R2 = 0.5 (not shown). The assumption is that monoterpenes with highest proton affinities could perform a quite exothermic direct proton transfer or form an energetically excited (NH4-A)+* intermediate that quickly dissociates forming AH+ + NH3. We therefore expected to find no cluster ion signal for sabinene (PA = 215.9 kcal/mol). As already discussed and shown in Supplementary Figure 1 the cluster ion yield for sabinene is solely 5% under dry conditions.
Limonene (PA = 201.2 kcal/mol) has the smallest PA, which is even smaller than NH3, and should form exclusively adduct ions only. We observe a cluster ion yield of 85%, which is the highest one of all monoterpenes investigated. But still, there exists a 15% channel at 0.057 eV collision energy producing protonated limonene. At 0.084 eV and dry conditions the yield of the proton transfer channel (including fragment ions) increases to 30% at the lowest extraction voltage (20 V). Increasing the extraction voltage from 20 to 30 V increases this channel even further reaching 50% (Figure 9). This means that the collision energy in the drift tube also leads to additional excitation of the (NH4-A)+*- intermediate, which is needed for the AH+ + NH3 channel to become thermodynamically accessible. Increasing the extraction voltage at 0.084 eV collision energy shows a further decrease of the (A) cluster ions, and a gain of AH+ and corresponding fragment ions. This is seen not only for sabinene but also for all investigated monoterpenes (Figures 5–9 and Supplementary Figures 1–3). This means that even stabilized (A) cluster ions (and AH+ ions) keep more internal energy at higher collision energies in the drift tube and it is then easier to form additional AH+ (and corresponding fragment ions) in the ion transfer region through collision induced dissociation (CID). We found the best correlation coefficient (R2 = 0.79) to predict the yield of the cluster ion formation channel for monoterpenes (see Figure 10) when the fraction of measured adduct ions was correlated as a function of the reaction enthalpy ΔHr of reaction (12), which is the difference in proton affinities plus the bond energy BE: ΔHr = PA(NH3) – PA(A) + BE(-A). This is shown in Figure 10, meaning that highest adduct ion yields are formed when compound A is strongly bond to and the difference in proton affinity is small resulting in smallest internal energy of the (NH4A)+* intermediate.
Figure 10. The reaction enthalpy ΔHr of the reaction , where A is a monoterpene, shows a good correlation with the -A cluster ion yield.
Table 1 gives an overview of the calculated collisional rate coefficient (kc) at the respective collision energy (KEcm), the calculated sensitivity (εcalc) using kc, the measured sensitivity (εmeas) and the reaction efficiency (eff) for all compounds studied under dry and humid conditions and at 51 and 81 Td. The monoterpenes show reaction efficiencies ranging from 18.4% (ocimene) to 34.6% (camphene) at dry conditions and a KEcm of 0.057 eV (51 Td). These efficiencies take into account the proton transfer channel, as well as the adduct ion formation channel. Limonene (adduct channel 85%) has a reaction efficiency (eff.) of 20.2%, which indicates a very high effective binary rate coefficient of the cluster channel, which is most likely due to the long lifetime of the intermediate (NH4A)+* against unimolecular decomposition. Also 3-carene (eff. 32.6%) and myrcene (eff. 21.3%) both having an adduct channel yield of 50% have very high effective binary rate coefficients. Compared to the three ketones (acetone, MVK and MEK), the monoterpenes in general have larger effective binary rate coefficients (cluster ion channel) at 0.057 eV. While in the case of the ketones, no reaction products at all are observed at elevated collision energies (KEcm = 0.084 eV; 81 Td), the adduct ion channel for the monoterpenes is somewhat smaller at this energy compared to 0.057 eV (51 Td) but still quite prominent. The intrinsic difference between the ketone and the monoterpene reaction system is that in the case of the ketones, no rearrangement (exothermic proton transfer) is thermodynamically accessible in the (NH4A)+* intermediate. The proton in the (NH4A)+* intermediate, when A is a monoterpene, is not strictly localized at the NH3. This can be anticipated comparing the calculated structures of the protonated monoterpenes (Table 2) and the geometries of (A) cluster ions. The delocalization of the charge in the intermediate could enhance the lifetime of the intermediate. Another explanation for the enhanced lifetime of the (NH4A)+* intermediate is the higher molecular complexity in the case of A being a monoterpene compared to small ketones.
Effect of Humidity
The product ion distribution was slightly shifted when changed from dry to humid conditions. More cluster ions (~ 3%) were found under humid conditions (Figures 5–9, Supplementary Figures 1–3). As mentioned before, more water molecules in the sample air lead to an increased formation of (H2O) reagent ions. The hydrated ammonium might be able to undergo ligand switching with the monoterpenes, which could explain the slight increase of cluster ions.
Selective Detection of Monoterpene Isomers
Since more than a decade different methods have been tested to differentiate the monoterpene isomers using H3O+-CIMS technology. Müller et al. (2009) coupled a PTR front part to a Triple Quadrupole Tandem MS and a Linear Ion Trap and performed MS/MS studies. The CID spectra of mass selected protonated monoterpene isomers were too similar to specify individual monoterpenes in complex mixtures. Misztal et al. (2012) operated a PTR-MS in an alternating drift voltage mode using 9 ascending and 9 descending voltage steps. The different “time points” of fragmentation of the resulting monoterpene ion signals were used to differentiate between the isomers. Their conclusion was that this method extends the selectivity of the PTR-MS method but cannot compete with the gold standard of GC-MS identification. Here we presented another possibility to separate monoterpene isomers due to their different cluster ion formation in reactions with reagent ions. We used eight different monoterpenes and found that limonene produces with 85% yield (A) cluster ions while sabinene and camphene have a corresponding cluster ion yield of <5%. In order to predict the yield of the cluster ion formation channel of monoterpenes we propose using the reaction enthalpy of the reaction as shown in Figure 10. SRI-TOF-MS allows to switch between different reagent ions such as H3O+ and which could be used to extend the selectivity of this CIMS method. But the challenge in separating all monoterpene isomers remains. The fundamental problem arises from the large number of possible isomers that are present at the same time in the real atmosphere.
In this laboratory study, we investigated the reactions of with a series of organic analytes (A): acetone (C3H6O), methyl vinyl ketone (C4H6O), methyl ethyl ketone (C4H8O) and eight monoterpene isomers (C10H16). The reactivity and product ion distribution were studied at two different collision energies and as a function of absolute humidity. Compounds having a lower proton affinity than NH3 produced only cluster ions (A). This is the case for the ketones acetone, MVK and MEK, which were observed only at a low collision energy of 0.055 eV. At an elevated collision energy of 0.080 eV no cluster ions of the carbonyls could be detected, meaning that these product ions are formed by association reactions, which are strongly temperature dependent in agreement with earlier observations of Adams et al. (2003). Collision induced dissociation of cluster ions has been studied by varying the extraction voltage applied between the drift tube and the TOF mass analyzer giving a first indication of cluster ion bond energies. Bond energies of cluster ions and proton affinities for most of the compounds used here are not known and have been estimated in the present study by high level quantum chemical calculations. In addition to cluster ion formation, also proton transfer reactions were observed for compounds having a higher proton affinity than that of NH3. The monoterpenes have proton affinities ranging from slightly lower to substantially higher than NH3. Calculated proton affinities and cluster bond energies allow to group these compounds as a function of the enthalpy for the dissociation reaction . We find that this enthalpy can be used for the monoterpenes to predict the (A) clusters ion yield. The present study explains product ion formation involving ion chemistry. This is of importance for chemical ionization mass spectrometry (CIMS) utilizing as well as (H2O) as reagent ions to detect pure hydrocarbon precursor having at least one C = C double bond as well as oxygenated organic compounds in real-time (Berndt et al., 2018a). Here we demonstrated that not only carbonyl compounds, but also hydrocarbons having a higher proton affinity than NH3, such as the monoterpenes, can be quantitatively detected with reagent ions.
The raw data supporting the conclusions of this manuscript will be made available by the authors, without undue reservation, to any qualified researcher.
BS and EC ran the experiments and analyzed the data. NH performed the quantum chemical calculations. EC, NH, BS, LF, and AH took part in the data discussion. LF implemented the raw data analysis software. AH and EC wrote the manuscript. All authors contributed to improvements of the manuscript.
This study has received financial funding from the Austrian Ministry for Science and Research (BMBWF Austria) within the program Sparkling Science (SPA 06/222 – CHAMPIONS). NH thanks the European Research Council (Grant No. 692891-DAMOCLES) for funding and CSC-IT Center for Science, Finland, for computational resources.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fchem.2019.00191/full#supplementary-material
Abedi, A., Sattar, L., Gharibi, M., and Viehland, L. A. (2014). Investigation of temperature, electric field and drift-gas composition effects on the mobility of ions in He, Ar, N2, and CO2. Int. J. Mass Spectrom. 370, 101–106. doi: 10.1016/j.ijms.2014.06.014
Adams, N. G., Babcock, L. M., Mostefaoui, T. M., and Kerns, M. S. (2003). Selected ion flow tube study of association and of product switching reactions with a series of organic molecules. Int. J. Mass Spectrom. 223–224, 459–471. doi: 10.1016/S1387-3806(02)00932-6
Berndt, T., Mentler, B., Scholz, W., Fischer, L., Herrmann, H., Kulmala, M., et al. (2018a). Accretion product formation from ozonolysis and OH radical reaction of α-pinene: mechanistic insight and the influence of isoprene and ethylene. Environ. Sci. Technol. 52, 11069–11077. doi: 10.1021/acs.est.8b02210
Berndt, T., Scholz, W., Mentler, B., Fischer, L., Herrmann, H., Kulmala, M., et al. (2018b). Accretion product formation from self- and cross-reactions of RO2 radicals in the atmosphere. Angew. Chem. Int. Ed. 57, 3820–3824. doi: 10.1002/anie.201710989
Breitenlechner, M., Fischer, L., Hainer, M., Heinritzi, M., Curtius, J., and Hansel, A. (2017). PTR3: an instrument for studying the lifecycle of reactive organic carbon in the atmosphere. Anal. Chem. 89, 5824–5831. doi: 10.1021/acs.analchem.6b05110
Cappellin, L., Karl, T., Probst, M., Ismailova, O., Winkler, P. M., Soukoulis, C., et al. (2012). On quantitative determination of volatile organic compound concentrations using proton transfer reaction time-of-flight mass spectrometry. Environ. Sci. Technol. 46, 2283–2290. doi: 10.1021/es203985t
Frege, C., Ortega, I. K., Rissanen, M. P., Praplan, A. P., Steiner, G., Heinritzi, M., et al. (2018). Influence of temperature on the molecular composition of ions and charged clusters during pure biogenic nucleation. Atmos. Chem. Phys. 18, 65–79. doi: 10.5194/acp-18-65-2018
Graus, M., Müller, M., and Hansel, A. (2010). High resolution PTR-TOF: quantification and formula confirmation of VOC in real time. J. Am. Soc. Mass Spectrom. 21, 1037–1044. doi: 10.1016/j.jasms.2010.02.006
Hansel, A., Jordan, A., Holzinger, R., Prazeller, P., Vogel, W., and Lindinger, W. (1995). Proton transfer reaction mass spectrometry: on-line trace gas analysis at the ppb level. Int. J. Mass Spectrom. Ion Process. 149–150, 609–619. doi: 10.1016/0168-1176(95)04294-U
Hansel, A., Scholz, W., Mentler, B., Fischer, L., and Berndt, T. (2018). Detection of RO2 radicals and other products from cyclohexene ozonolysis with and acetate chemical ionization mass spectrometry. Atmos. Environ. 186, 248–255. doi: 10.1016/j.atmosenv.2018.04.023
Hansel, A., Singer, W., Wisthaler, A., Schwarzmann, M., and Lindinger, W. (1997). Energy dependencies of the proton transfer reactions H3O+ + CH2O < = > CH2OH+ + H2O*. Int. J. Mass Spectrom. Ion Process. 167–168, 697–703. doi: 10.1016/S0168-1176(97)00128-6
Kebarle, P., Searles, S. K., Zolla, A., Scarborough, J., and Arshadi, M. (1967). Solvation of the hydrogen ion by water molecules in the gas phase. Heats and entropies of solvation of individual reactions. H+(H2O)n−1 + H2O -> H+(H2O)n. J. Am. Chem. Soc. 89, 6393–6399. doi: 10.1021/ja01001a001
Lindinger, W., Hansel, A., and Jordan, A. (1998a). On-line monitoring of volatile organic compounds at pptv levels by means of proton-transfer-reaction mass spectrometry (PTR-MS) medical applications, food control and environmental research. Int. J. Mass Spectrom. Ion Process. 173, 191–241. doi: 10.1016/S0168-1176(97)00281-4
Lindinger, W., Hansel, A., and Jordan, A. (1998b). Proton-transfer-reaction mass spectrometry (PTR–MS): on-line monitoring of volatile organic compounds at pptv levels. Chem. Soc. Rev. 27:347. doi: 10.1039/a827347z
Maquestiaut, A., Flammang, R., Nielsens, L., Maquestiau, A., Flammang, R., and Nielsen, L. (1980). A study of some cations formed in the ammonia chemical ionization of ketones using mass analysed ion kinetic energy spectrometry. Org. Mass Spectrom. 15, 376–379. doi: 10.1002/oms.1210150712
Materić, D., Lanza, M., Sulzer, P., Herbig, J., Bruhn, D., Gauci, V., et al. (2017). Selective reagent ion-time of flight-mass spectrometry study of six common monoterpenes. Int. J. Mass Spectrom. 421, 40–50. doi: 10.1016/j.ijms.2017.06.003
Misztal, P. K., Heal, M. R., Nemitz, E., and Cape, J. N. (2012). Development of PTR-MS selectivity for structural isomers: monoterpenes as a case study. Int. J. Mass Spectrom. 310, 10–19. doi: 10.1016/j.ijms.2011.11.001
Müller, M., Mielke, L. H., Breitenlechner, M., McLuckey, S. A., Shepson, P. B., Wisthaler, A., et al. (2009). MS/MS studies for the selective detection of isomeric biogenic VOCs using a townsend discharge triple quadrupole tandem MS and a PTR-linear ion trap MS. Atmos. Meas. Tech. 2, 703–712. doi: 10.5194/amt-2-703-2009
Payzant, J. D., Cunningham, A. J., and Kebarle, P. (1973). Gas phase solvation of the ammonium ion by NH3 and H2O and stabilities of mixed clusters (NH3)n(H2O)w. Can. J. Chem. 51, 3242–3249. doi: 10.1139/v73-485
Schoon, N., Amelynck, C., Vereecken, L., Coeckelberghs, H., and Arijs, E. (2004). A selected ion flow tube study of the reactions of H3O+, NO+ and with some monoterpene oxidation products. Int. J. Mass Spectrom. 239, 7–16. doi: 10.1016/j.ijms.2004.09.003
Sindelarova, K., Granier, C., Bouarar, I., Guenther, A., Tilmes, S., Stavrakou, T., et al. (2014). Global data set of biogenic VOC emissions calculated by the MEGAN model over the last 30 years. Atmos. Chem. Phys. 14, 9317–9341. doi: 10.5194/acp-14-9317-2014
Smolander, S., He, Q., Mogensen, D., Zhou, L., Bäck, J., Ruuskanen, T., et al. (2014). Comparing three vegetation monoterpene emission models to measured gas concentrations with a model of meteorology, air chemistry and chemical transport. Biogeosciences. 11, 5425–5443. doi: 10.5194/bg-11-5425-2014
Solouki, T., and Szulejko, J. E. (2007). Bimolecular and unimolecular contributions to the disparate self-chemical ionizations of α-pinene and camphene isomers. J. Am. Soc. Mass Spectrom. 18, 2026–2039. doi: 10.1016/j.jasms.2007.08.016
Tani, A., Hayward, S., Hansel, A., and Hewitt, C. N. (2004). Effect of water vapour pressure on monoterpene measurements using proton transfer reaction-mass spectrometry (PTR-MS). Int. J. Mass Spectrom. 239, 161–169. doi: 10.1016/j.ijms.2004.07.020
Tereza Fernandez, M., Williams, C., Mason, R. S., and Costa Cabral, B. J. (1998). Experimental and theoretical proton affinity of limonene. J. Chem. Soc. Faraday Trans. 94, 1427–1430. doi: 10.1039/a800049b
Wang, T., Španěl, P., and Smith, D. (2003). Selected ion flow tube, SIFT, studies of the reactions of H3O+, NO+ and with eleven C10H16 monoterpenes. Int. J. Mass Spectrom. 228, 117–126. doi: 10.1016/S1387-3806(03)00271-9
Werner, H.-J., Knowles, P. J., Knizia, G., Manby, F. R., Schütz, M., Celani, P., et al. (2015). MOLPRO, Version 2015.1, A Package of ab Initio Programs. Available at: http://www.molpro.net.
Zhou, S., Rivera-rios, J. C., Keutsch, F. N., and Abbatt, J. P. D. (2018). Identification of organic hydroperoxides and peroxy acids using atmospheric pressure chemical ionization-tandem mass spectrometry (APCI-MS/MS): application to secondary organic aerosol. Atmos. Meas. Tech. 11, 3081–3089. doi: 10.5194/amt-11-3081-2018
Keywords: NH4+, chemical ionization, PTR-ToF-MS, association reactions, monoterpenes, acetone, methyl vinyl ketone (MVK), methyl ethyl ketone (MEK)
Citation: Canaval E, Hyttinen N, Schmidbauer B, Fischer L and Hansel A (2019) NH4+ Association and Proton Transfer Reactions With a Series of Organic Molecules. Front. Chem. 7:191. doi: 10.3389/fchem.2019.00191
Received: 07 January 2019; Accepted: 13 March 2019;
Published: 03 April 2019.
Edited by:Paolo Tosi, University of Trento, Italy
Reviewed by:Oh-Hoon Kwon, Ulsan National Institute of Science and Technology, South Korea
Fuminori Misaizu, Tohoku University, Japan
Pawel K. Misztal, University of California, Berkeley, United States
Copyright © 2019 Canaval, Hyttinen, Schmidbauer, Fischer and Hansel. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Armin Hansel, firstname.lastname@example.org
†Present Address: Noora Hyttinen, Nano and Molecular Systems Research Unit, University of Oulu, Oulu, Finland