Anal. Chem. 2004, 76, 5436-5440
Development of Ion Drift-Chemical Ionization Mass Spectrometry Edward C. Fortner, Jun Zhao, and Renyi Zhang*
Department of Atmospheric Sciences, Texas A&M University, College Station, Texas 77843
An ion drift-chemical ionization mass spectrometry (IDCIMS) technique has been developed to detect and quantify trace gases, including volatile organic compounds and inorganic species. The trace species are chemically ionized into positive or negative product ions with a wellcontrolled ion-molecule reaction time. The ID-CIMS method allows for quantification of the trace gases without the necessity of performing calibrations with authentic standards for the trace gases. Demonstrations of the ability of ID-CIMS to accurately quantify isoprene and HNO3 in a laboratory setting are presented. The results illustrate that the ID-CIMS technique facilitates detection and quantification of organic and inorganic species in laboratory kinetic investigations and field measurements. Chemical ionization mass spectrometry (CIMS) has been widely employed in atmospheric chemistry research, including laboratory kinetic investigations and field trace gas measurements.1-6 In the CIMS method, a neutral species is ionized by a reagent ion to yield a product ion, which is analyzed by a quadrupole mass analyzer for species identification and abundance determination. One type of CIMS applications involves proton-transfer reaction mass spectrometry (PTR-MS), which was developed originally by Lindinger and co-workers for on-line quantification of volatile organic compounds (VOCs) in air.5 There is a wide variety of sources of VOCs to the Earth’s atmosphere ranging from natural to anthropogenic in origin. Many questions exist concerning the ultimate fate of these compounds in the atmosphere. Understanding the VOC loss processes is critical to resolving the atmospheric budget of organic carbon.7 The photooxidation of VOCs also plays a key role in the processes of tropospheric ozone and secondary organic aerosol generation.8,9 These two major air pollutants have * Corresponding author. E-mail:
[email protected]. (1) Zhang, R.; Lei, W.; Molina, L. T.; Molina, M. J. Int. J. Mass Spectrom. 2000, 194, 41-48. (2) Zhang, R.; Leu, M. T.; Keyser, L. F. Geophys. Res. Lett. 1995, 22, 14931496. (3) Zhang, D.; Zhang, R.; North, S. W. J. Phys. Chem. 2003, 107, 11013-11019. (4) Huey, L. G.; Dunlea, E. J.; Lovejoy, E. R.; Hanson, D. R.; Norton, R. B.; Fehsenfeld, F. C.; Howard, C. J. J. Geophys. Res. 1998, 103, 3355-3360. (5) Lindinger, W.; Hansel, A.; Jordan, A. Int. J. Mass Spectrom. 1998, 173, 191241. (6) de Gouw, J.; Warneke, C.; Karl, T.; Eerdekens, G.; van der Veen, C.; Fall, R. Int. J. Mass Spectrom. 2003, 223, 365-382. (7) Hewitt, C. N., Ed. Reactive Hydrocarbons in the Atmosphere; Academic Press: San Diego, 1999. (8) Zhao, J.; Zhang, R.; Fortner, E. C.; North, S. W. J. Am. Chem. Soc. 2004, 126, 2686-2687.
5436 Analytical Chemistry, Vol. 76, No. 18, September 15, 2004
detrimental health effects and have adverse consequences to the Earth’s climate. In the PTR-MS method, proton-transfer reactions with hydronium ions (H3O+) are employed to ionize VOCs, and the resulting protonated product ion is then detected by mass spectrometry.5,6 VOCs directly emitted from natural and anthropogenic sources as well as their oxidation products can be quantified simultaneously by PTR-MS. Some desirable features offered by PTR-MS that differentiate it from the traditional GC/ MS include a fast-time response on the order of seconds or less, a lack of air sample preconcentration and chromatography, and the ability to perform unmanned on-line measurements over extended time periods.5,6 A recent experimental study has also demonstrated the application of the PTR-MS method in quantifying product yields of hydrocarbon oxidation reactions.8 This method allows product yield measurements without the necessity of calibration and hence is advantageous because of the general difficulty to obtain the authentic standards for many products of hydrocarbon reactions. In particular, recent advances have been made in quantum-chemical calculations to investigate the oxidation mechanism of VOCs.10-15 The PTR-MS application of measuring product yields without calibration using authentic standards offers a great potential for validation of the theoretical predictions, which could have a profound impact on elucidation the VOC chemistry in the atmosphere. The CIMS method has also been widely used in atmospheric monitoring of inorganic species.4,16,17 For example, nitric acid is a critical member of the NOy (NO, NO2, NO3, HNO2, HNO3, HNO4, peroxyacetyl nitrate, and organic and inorganic nitrates) family that is created in the troposphere predominantly by the association of OH with NO2. The transformation of NOx (NO + NO2) to HNO3 is the primary scavenging method of NOx from the atmosphere. Due to the fundamental role that NOx plays in the production and destruction of O3 in the atmosphere, it is critical that the magnitude of this transformation is understood in order to resolve the global ozone budget. Huey and co-workers originally developed a selective and high-sensitivity CIMS method for in situ (9) Zhang, R.; Suh, I.; Zhao, J.; Zhang, D.; Fortner, E. C.; Tie, X.; Molina, L. T.; Molina, M. J. Science 2004, 304, 1487-1490. (10) Lei, W.; Zhang, R. J. Chem. Phys. 2000, 113, 153-157. (11) Lei, W.; Derecskei-Kovacs, A.; Zhang, R. J. Chem. Phys. 2000, 113, 53545360. (12) Suh, I.; Zhang, D.; Zhang, R.; Molina, L. T.; Molina, M. J. Chem. Phys. Lett. 2002, 364, 454-462. (13) Zhang, D.; Zhang, R. J. Chem. Phys. 2002, 116, 9721-9728. (14) Zhang, D.; Lei, W.; Zhang, R. Chem. Phys. Lett. 2002, 358, 171-179. (15) Zhang, D.; Zhang, R.; Allen, D. T. J. Chem. Phys. 2003, 118, 1794-1801. (16) Furutani, H.; Akimoto, H. J. Geophys. Res. 2002, 107, 4016. (17) Huey, G. L.; Lovejoy, E. R. Int. J. Mass Spectrom. 1996, 155, 133-140. 10.1021/ac0493222 CCC: $27.50
© 2004 American Chemical Society Published on Web 08/14/2004
Figure 1. Schematic depiction of the ID-CIMS apparatus.
measurements of HNO3.4,17 It has been shown that this method allows for detection of HNO3 with a time resolution on the order of a few seconds and a detection sensitivity of the order of partsper-trillion (ppt) levels.4,16,17 In this paper, we describe the development of an ion driftchemical ionization mass spectrometry (ID-CIMS) instrument, which is capable of providing on-line measurements of both VOCs and inorganic species. The method inherits the CIMS scheme yet with a controlled ion-molecule reaction time, hence extending the operation of PTR-MS for positive and negative reagent ions. A description of the instrument configuration follows. The methodology used to quantify VOCs and inorganic species will be presented, and the differences between the two detection schemes will be detailed. The results of laboratory calibrations of isoprene and HNO3 will then be presented.
1.4 cm through which the flow passed. The chemical ionization reactions between the reagent ion and the neutral compound to be detected and quantified occurred in this region. A voltage was applied to the rings to enhance the ion flow. At the downstream end of the drift tube, the majority of the flow was diverted into a pump while a small portion of the flow passed through an aperture (0.2 mm) into the MS analyzer (Extrel QC-150). The aperture was slightly biased (0-10 V) for ion-tuning purposes, and the polarity was dependent on the polarity of the reagent ions produced. The pressure in the drift tube was regulated between 1 and 4 Torr. The sampled ions were focused with an ion optics, analyzed with a quadrupole mass filter, and detected by an electron multiplier. The multiplier was connected with a preamplifier, which converted the raw signal into TTL pulses allowing them to be processed by a personal computer.
INSTRUMENTATION DESCRIPTION The ID-CIMS system consists of an ion source to produce positive or negative reagent ions, a drift tube where the ionmolecule reaction takes place, and a quadrupole mass spectrometer where the reagent ions are analyzed, as depicted schematically in Figure 1. In the current configuration, a flow tube is connected to the ID-CIMS for laboratory measurements and calibrations, and an inlet will be needed for its use in the field. For the purpose of the calibrations conducted in our laboratory, the flow tube was a Pyrex tubing 70 cm in length and 2.25-cm i.d. A coating of halocarbon wax was applied to the inner walls of the flow tube to reduce wall loss. The majority of the gas in the flow tube was diverted into a mechanical pump at the downstream end of the flow tube while only a small amount (∼1%) of the total flow passed through a 1-mm opening into the drift tube region. The pressure in the flow tube was maintained at ∼100 Torr. The ion source was a custom-made glow discharge device. For VOC detection, ambient air flowed through the discharge, and a positive voltage between 1 and 2 kV was applied to the electrode in the discharge producing H3O+ ions. For HNO3 detection, nitrogen flowed through the discharge mixed with SF6 at the ppm level, and a negative voltage (-5 kV) was applied to the electrode producing SF6- ions. The drift tube was constructed of a 9.5-cm-length Pyrex tube that contained a set of 10 stainless steel rings connected in series with 1-MΩ resistors between rings. The rings had an i.d. of
METHODOLOGY The principle that ID-CIMS utilizes is the detection of a neutral molecule X according to the ion-molecule reaction
A(+ or -) + X f X′(+ or -) + others,
k1
(1)
where A is the reagent ion, X′ is the product ion, and k1 is the ion-molecule reaction rate constant. The abundance of the neutral species X in the drift tube can then be quantified from the equation
[X′]∆t ) k1[A][X]∆t
(2)
∆t ) l/Ut
(3)
where
and l is the length of the drift tube and Ut is the total velocity of the reagent ion, determined from the ion drift velocity (Ui) and flow velocity (Uf) in the drift tube. The concentration of the neutral species in the flow tube, [Xft], can be expressed as
[Xft] ) P/∆P[X]
(4)
where P is pressure in the flow tube and ∆P is the pressure change in the drift tube due to addition of the gases from the Analytical Chemistry, Vol. 76, No. 18, September 15, 2004
5437
flow tube. ∆P was determined by monitoring the pressure rise in the drift tube upon addition of the gases from the flow tube. The ion drift velocity Ui is determined by the equation
Ui ) uE
(5)
where u is the ionic mobility of the ion and E is the voltage gradient in the drift tube. The ionic mobility is determined by
u ) u0(760/P)(T/273)
(6)
where P and T are respectively the pressure and temperature within the drift tube and u0 is the reduced ionic mobility taken from the literature.18 If there is no mass-dependent discrimination in the MS detection efficiency, the concentrations of product ions and the reagent ions are proportional to their respective ion count signals measured by the mass spectrometer. Hence, the concentration of X can be determined from the signal ratio of the product ion to the reagent ion according to eq 2. It has been suggested that no change in detection efficiency was measured in the mass range of 19-250 Da.5 For the ∆t term, the work by de Gouw et al. has shown that the reaction time can be calculated according to the ion drift velocity, which is in agreement with the measurement.6 Hence, eqs 2 and 4 allow for the quantification of the gas-phase concentrations without the necessity of calibration, if the ionmolecule reaction rate constant for eq 1 is accurately obtained. The ion-molecule reaction rate constants can be determined experimentally or theoretically.5,19-22 Recently, we have determined the proton-transfer rate constants for a suite of atmospherically important VOCs using the average dipole orientation (ADO) theory according to their permanent dipole moment and polarizability calculated from ab initio calculations.22 For VOC detection, H3O+ is generally used as the reagent ion and the proton-transfer reaction
H3O+ + X f XH+ + H2O,
(7)
k7
occurs in the drift tube, an approach commonly referred to as PTR-MS.5 The proton-transfer reaction is very efficient for hydrocarbons and oxygenated VOCs with a proton affinity greater than that of water (165.5 kcal mol-1).22 If water clusters are present in the drift tube the reaction
H3O+‚(H2O)n + X f XH+ + H2O(n + 1),
k8
(8)
H2O with X, respectively. The fractional percentage of the entire reagent ion signal is represented by a for H3O+ and b for H3O+‚ H2O. For PTR-MS, E was kept sufficiently large so that Ui . Uf in the drift tube and Uf was neglected. Within ∆t, Ui must be modified to account for the fractional percentages of the different reagent ions. The ionic mobility of the two reagent ions is different and this affects Ui. The reduced ionic mobilities are 2.76 and 2.28 cm2 V-1 s-1 for H3O+ and H3O+‚H2O, respectively.18 It has been experimentally determined that some hydrocarbons (such as aromatic hydrocarbons toluene and benzene) do not react with water clusters so eq 7 can be used to quantify these species regardless of the percentage of water cluster where H3O+ is used as the reagent ions.21 The method of employing PTR-MS to quantify atmospheric VOCs has been thoroughly discussed by Lindinger et al.5 and de Gouw et al.,6 while the PTR-MS method to quantify products from hydrocarbon oxidation reactions has been described by Zhao et al.8 For HNO3 detection, the ion-molecule reaction sequence
SF6- + SiF4 f SiF5- + SF5, k10 ) (6.7 ( 1.5) × 10-10 cm3 molecule-1 s-1 (10) SiF5- + HNO3 f SiF5-‚HNO3, k11 ) (3.8 ( 1) × 10-10 cm3 molecule-1 s-1 (11) is used to produce the reagent ion.17 It has been suggested that a potential product of the reaction of SiF5- with HNO3 is SiF4‚NO3-.17 However, we did not encounter a signal at this ion mass, and on the basis of calibrations, we believed that the impact of this potential ion-molecule reaction was minimal. It is likely though that this product ion formation is suppressed in the drift tube. HNO3 is thus quantified using
[HNO3] ) [SiF5-‚HNO3]/(k11[SiF5-]∆t)
(12)
The voltage gradient E was kept small in this case so Uf and Ui must be added to obtain ∆t. The reduced ionic mobility of SiF5in N2 has not been experimentally determined. We estimated its value on the basis of the literature value for SF6- in N2, assuming an exponential mass-dependent expression between the masses of the two respective ions. The literature value for the ionic mobility of SF6- in N2 is 1.9 cm2 V-1 s-1,23 and we estimated a value of 2.0 cm2 V-1 s-1 for SiF5- in N2.
where k7 and k8 are the reaction rate constants for H3O+ and H3O+‚
RESULTS AND DISCUSSION To demonstrate the application of the ID-CIMS technique for detection and quantification of organic and inorganic compounds, we performed experimental evaluations of the ID-CIMS method for isoprene and HNO3. The proton-transfer reactions between H3O+ and the majority of VOC molecules generally proceed nondissociately to produce the protonated molecules.5 However, for the reaction between H3O+ and isoprene, our observation showed that other ionic products were formed, dependent on the E/N ratio, where E is
(18) Dotan, I.; Albritton, D. L.; Lindinger, W.; Pahl, M. J. Chem. Phys. 1976, 65, 5028-5030. (19) Su, T.; Bowers, M. T. J. Chem. Phys. 1973, 58, 3027-3037. (20) Zhang, R.; Lei, W. J. Chem. Phys. 2000, 113, 8574-8579.
(21) Warneke, C.; van der Veen, C.; Luxembourg, S.; de Gouw, J. A.; Kok, A. Int. J. Mass Spectrom. 2001, 207, 167-182. (22) Zhao, J.; Zhang, R. Atmos. Environ. 2004, 38, 2177-2185. (23) Hinojosa, G.; de Urquijo, J. J. Phys. D: Appl. Phys. 2003, 36, 2510-2514.
will occur in the drift tube when energetically allowed. If water clusters can be primarily limited to n ) 1, the neutral species X can still be reliably quantified as
[XH′] ) (k7a[H3O+] + k8b[H3O+‚H2O])[X]∆t
5438
(9)
Analytical Chemistry, Vol. 76, No. 18, September 15, 2004
Figure 2. Spectral scan showing the H3O+ reagent ion, along a small H2O‚H3O+ peak, and protonated isoprene C5H9+. The isoprene peak is multiplied by a factor of 20 for clarity.
the electric gradient (V cm-1) and N is the number concentration of the buffer gas (molecule cm-3). A considerable peak at mass 39 was observed, and a significant peak at mass 67 was also obtained at an E/N ratio larger than 140 Td (1 Td ) 1017 V cm2 molecule-1) when using air as the carrier gas. These two mass peaks increased when a higher E/N ratio was applied. This implies that there are three possible channels for the reaction between H3O+ and isoprene at high E/N ratios,
C5H8 + H3O+ f C5H9+ + H2O
(13a)
f C5H7+ + H2 + H2O
(13b)
f C3H3+ + C2H6 + H2O
(13c)
Hence, a proper E/N ratio must be chosen in order to inhibit the dissociation of the protonated isoprene and ensure the major channel for the desired proton-transfer reaction. On the other hand, if the E/N ratio is too low, a significant amount of water ion cluster will be formed, which would complicate the interpretation of the mass spectra. In this experiment, we kept the E/N ratio at ∼100 Td, which greatly prevented the protoned isoprene from being dissociated while the water clusters were effectively minimized. Figure 2 shows the spectral scan from 10 to 100 amu. The three peaks at mass 19, 37, and 69 correspond to H3O+, H3O‚ H2O+, and protonated isoprene C5H9+, respectively. The H3O‚H2O+ signal was ∼10% of the H3O+ signal, while the water cluster H3O+‚ (H2O)2 had a much lower signal and was neglected along with all other higher order water clusters. We conducted calibrations using volumetrically prepared samples of helium doped with 0.1% isoprene (Aldrich 99.5%).24 For the ion-molecule reaction between H3O+ and isoprene, we used a value of 1.94 × 10-9 cm3 molecule-1 s-1.21 Using the ADO theory, we calculated a rate constant of 1.5 × 10-9 cm3 molecule-1 s-1 for the reaction between H3O‚H2O+ and isoprene. For isoprene calibration, the pressures were about 100 and 3 Torr in the flow tube and drift tube, respectively. A carrier gas flow of ∼26 standard liters per minute (SLPM) entered the flow tube, and a turbulent flow condition was effectively maintained.24 A small isoprene sample flow (0-2 standard cubic centimeters per minute, sccm) was introduced into the flow tube. Air passed through the corona discharge and into the drift tube at 1 SLPM. As discussed above, (24) Zhang, R.; Suh, I.; Clinkenbeard, A. D.; Lei, W.; North, S. W. J. Geophys. Res. 2000, 105, 24627-24635.
Figure 3. Comparison of isoprene concentrations calibrated from the known volumetric mixing ratio of the gas standard (horizontal axis) and calculated using eq 9 (vertical axis).
the E/N ratio potentially impacts the isoprene detection: a higher ratio produces less water clusters but significant fragment as it results in more ion kinetic energy, which may be able to drive endothermic processes. For reagent ions that are supposed to form clusters, a higher E/N ratio may actually decrease detection sensitivity. For isoprene detection, the aperture plate was held at ground. Figure 3 shows an example of isoprene calibration. The horizontal axis corresponds to isoprene concentrations determined from the known volumetric mixing ratio of the gas standard in the flow reactor, and the vertical axis represents those measured by the PTR-MS method according to eq 9. Using an E/N ratio of ∼100 Td, an excellent agreement was achieved between the isoprene concentrations estimated from the known volumetric mixing ratio of the gas standard in the flow reactor and measured by the PTR-MS method. The correlation coefficient (R2) in the isoprene calibration is within 99.9%. This comparison indicates that isoprene can be accurately quantified using eq 9. A detection limit of ∼10 ppt (1 ppt corresponds to 2.45 × 107 molecule cm-3 at 298 K and 1 atm) was calculated with respect to isoprene for a unity ratio of signal-to-noise and 1-s integration time. For HNO3 detection, the drift tube pressure was at ∼3 Torr, and the flow tube was maintained under the turbulent flow condition. Nitrogen doped with SF6 at the ppm level was passed into the drift tube at a rate of 1.3 SLPM after passing through the discharge. A N2 flow of 50 sccm doped with SiF4 at the ppm level proceeded directly into the drift tube. A voltage gradient of 4.5 V cm-1 was applied in the drift tube, and the aperture had a small voltage (-8 V) applied to it to assist in tuning. A spectrum of the reagent ion SiF5- and product ion SiF5-‚HNO3 is shown in Figure 4 when a small (∼1 sccm) flow of a 0.1% HNO3 in helium sample was detected by the ID-CIMS apparatus. In our experiments, the residual signal due to SF6- was always much (by at least 2 orders of magnitude) smaller than that the SiF5- reagent ions, and its effect on HNO3 detection was negligible. Figure 5 shows the comparison between the HNO3 concentration calculated using eq 12 (vertical axis) and that determined by the standard flow consideration knowing the mixing ratio of the HNO3 sample (horizontal axis). The correlation coefficient in the HNO3 calibration was within 99.8%, indicating that HNO3 can be reliably Analytical Chemistry, Vol. 76, No. 18, September 15, 2004
5439
detect the desired product ions in our laboratory experiments several parameters will likely change in a field environment. It should be expected that variables such as the percentage of water clusters and the reagent ion signal intensity would likely be affected by the changes, potentially impacting the performance of the instrument under ambient conditions.
Figure 4. Spectral scan showing the reagent ion SiF5- and the product ion SiF5-‚HNO3. The HNO3 peak is multiplied by a factor of 100 for clarity. A small peak at 146 amu is due to the presence of SF6-.
Figure 5. Comparison of HNO3 concentrations calibrated from the known volumetric mixing ratio of the gas standard (horizontal axis) and calculated using eq 12 (vertical axis).
quantified using eq 12. A detection limit of ∼10 ppt was estimated with respect to HNO3 for a unity ratio of signal-to-noise and 1-s integration time. We have also conducted calibrations of several other organic (such as toluene) and inorganic (such as NO2) compounds that were not featured in this paper. For example, NO2 was detected as NO2- using the SF6- reagent ion. The calculated concentrations consistently corresponded to the calibrated values with the correlation coefficient within 99%. It is worth noting that while the calibrations show the ability of the ID-CIMS instrument to (25) Holzinger, R.; Sandoval-Soto, L.; Rottenberger, S.; Crutzen, P. J.; Kesselmeier, J. J. Geophys. Res. 2000, 105, 20573-20579.
5440 Analytical Chemistry, Vol. 76, No. 18, September 15, 2004
CONCLUSION The PTR-MS method, which allows for a controlled ionmolecule reaction with the hydronium ion within the drift tube, has been previously demonstrated for allowing on-line detection of atmospherically important VOCs with high detection sensitivity and fast-responding time.5,6,21,25 This method, however, is limited to detection of organic compounds whose proton affinity is larger than of water. On the other hand, the conventional CIMS approach does not employ a well-controlled ion-molecule reaction time.17 We have demonstrated in this paper the application of ID-CIMS, which expands the operation of PTR-MS to allow for the use of both negative and positive reagent ions. The trace species are chemically ionized into positive or negative product ions with a well-controlled ion-molecule reaction time; quantifying the trace species is possible, provided that the rate constants for the ionmolecule reaction are known. Our recent work showed that ionmolecule reaction rate constants can be reliably determined using theoretical approaches.22 The laboratory experiments presented in this work show the capability of ID-CIMS to execute on-line measurements of organic and inorganic compounds, as exemplified for isoprene and HNO3. Our results demonstrate that the ID-CIMS method allows for accurate quantification of trace gases without the necessity of performing calibrations with authentic standards. Hence, the IDCIMS technique facilitates detection and quantification of organic and inorganic species in laboratory kinetic studies and field measurements. Future work will focus on designing such a system to allow for implementation of ID-CIMS in a field environment. ACKNOWLEDGMENT This work was supported by the Robert A. Welch Foundation (A-1417) and Houston Advanced Research Center HARC (H-252004C2), and U.S. Environmental Protection Agency EPA (R030132). We thank Professor Mario J. Molina of MIT for helpful discussions. Three reviewers provided valuable comments for improving the manuscript.
Received for review May 7, 2004. Accepted June 28, 2004. AC0493222