Minimizing false positives in ion mobility spectrometry by changing the field or pressure Larry A. Viehland
International Journal for Ion Mobility Spectrometry ISSN 1435-6163 Int. J. Ion Mobil. Spec. DOI 10.1007/s12127-015-0181-0
1 23
Your article is protected by copyright and all rights are held exclusively by SpringerVerlag Berlin Heidelberg. This e-offprint is for personal use only and shall not be selfarchived in electronic repositories. If you wish to self-archive your article, please use the accepted manuscript version for posting on your own website. You may further deposit the accepted manuscript version in any repository, provided it is only made publicly available 12 months after official publication or later and provided acknowledgement is given to the original source of publication and a link is inserted to the published article on Springer's website. The link must be accompanied by the following text: "The final publication is available at link.springer.com”.
1 23
Author's personal copy Int. J. Ion Mobil. Spec. DOI 10.1007/s12127-015-0181-0
ORIGINAL RESEARCH
Minimizing false positives in ion mobility spectrometry by changing the field or pressure Larry A. Viehland 1
Received: 28 April 2015 / Revised: 4 June 2015 / Accepted: 5 June 2015 # Springer-Verlag Berlin Heidelberg 2015
Abstract A method is proposed for minimizing false positives in ion mobility spectrometry by changing the value of E=N , the ratio of the electrostatic field strength to the gas number density. The feasibility of this method is demonstrated by simulation of the mobilities of 2,4,6-trinitrotolune and 2methyl-3,5-dinitrophenol in helium. The zero-field mobilities of these ions are nearly the same, while their dependences upon E/N are so different that the ions can be distinguished by modest changes in E/N at 473 K and by somewhat larger changes at 300 K. Keywords False positives . Ion mobility spectrometry . Gaseous ion transport . (n,4)-core model potential
overlapping IMS peaks. Using a gas chromatograph as a preIMS separation device or adding chemical dopants can decrease the incidence of false positives [8], as can combining IMS with mass spectrometry [6]. The use of secondary electrospray ionization has been suggested [9] as a further improvement, one that avoids the need for sample separation before analysis or the addition of chemical dopants. The problems with these techniques include capital costs and the need for trained operators. The purposes of this paper are to propose a method by which false positives can be minimized by making rather simple changes in the instrument used, and to demonstrate the feasibility of this method by simulation of the mobilities of two ions that have similar mobilities.
Introduction Background Among the trace detection systems used at airports and in military and policing activities is ion mobility spectrometry (IMS) [10]. Since the consequence of not detecting a security threat might be disastrous, the operating parameters of IMS are chosen to make false negatives unlikely. Unfortunately, this can cause the number of false positives to increase, resulting in substantial costs in time and money [11]. A primary reason for false positives is that the number of chemical species present with similar mobilities can lead to This article is submitted for the special issue on IMS modeling and simulation. * Larry A. Viehland
[email protected] 1
Department of Science, Chatham University, Pittsburgh, PA 15232, USA
It has been known for many years [13] that the average energy of the ion-neutral collisions in an IMS is governed by the gas temperature, T , and the ratio, E=N , of the electric field strength to the gas number density. A rather recent formula [26] for when E/N becomes unimportant is E
σ it is attractive. The ability to consider both long-range attraction and short-range repulsion is an important advantage of using the (12,4)-core potential, in contrast to the hard-sphere
projection method [17]. The parameter n is a measure of the steepness of the repulsive (short-range) potential; n ¼ ∞ for a rigid sphere, it is usually about 8 for the soft potentials that characterize the interactions of small ions and neutrals, but here we use n ¼ 12 as originally suggested by Mason et al. [14]. To determine the other parameters in Eq. (2), we used ChemBio3D [7] to generate a model of the ion bonded to a N2 molecule and then another of the ion along with (but not bonded to) N2. For each, we then used the MOPAC [19] interface from ChemBio3D to find the minimum energy and corresponding structure. We set ε0 equal to the difference in these minimum energies. The value of a was obtained from ChemBio3D by rotating the energy-minimized structure for the non-bonded ion and N2, and setting it equal to the smallest dimension found. Finally, rm was obtained by noting that at large r Eq. (2) becomes V ðrÞ→
c4 ; r4
ð4Þ
where c4 ¼
3 ε0 ðrm aÞ4 : 2
ð5Þ
Since the value of c4 is the same for all ions in the same neutral gas, and since it is known [13] to have the value 5.938 in atomic units when the gas is nitrogen, we can use this value along with ε0 and a to determine rm from Eq. (5). The values for the parameters for each system are given in Table 1. We followed the path outlined above to determine the standard (or reduced) mobilities. Since we are using a simulation and a rather crude estimate of the interaction potential, it is not surprising that the zero-field mobilities are only approximately equal to the experimental values at 473 K. It is possible to systematically vary the parameters (including n ) in the (n,4)-core potential until a set is found for each ion-air system for which the calculated zero-field mobility at 473 K matches the experimental value. This was not done here because adjusting four parameters subject to only two constraints (Eq. 5 and the measured K 0 ) does not yield a unique solution. Instead, we scaled the values with a constant factor chosen so that K 0 ¼ 1:50 cmÇ=Vs for TNT+ in He at E=N ¼ 0, and 1.55 cm²/Vs for DNOC+. The results are shown in Fig. 1, where typical experimental error bars on each Table 1
Parameter values for the (12,4) core potentials
Parameter
TNT+
DNOC+
ε0 (millihartree)
1.219
1.937
a (bohr)
2.143
1.070
rm (bohr)
9.692
7.794
σ (bohr)
8.723
6.931
Author's personal copy Int. J. Ion Mobil. Spec.
1.70
particular that now the DNOC+ mobilities increase more with E=N , since the thermal energy is much smaller than before, while the TNT+ mobilities vary little with E=N , since the thermal energy is slightly below ε0, a small increase in E=N is needed to surmount the mobility maximum, and only at still higher E=N values are the mobilities probing the repulsive part of the potential.
K0
1.60 1.50 1.40
1.30 1.20
Discussion
1.10
E/N
1.00 0
20
40
60
80
100
Fig. 1 Simulated standard mobilities, K 0 in cm2/Vs, of TNT+ (red curve) and DNOC+ (blue points) in He at 473 K, as functions of E=N in Td
set of results have been replaced by a single set of 5 % error bars on the DNOC+ results. What is immediately obvious is that it does not take a very large value of E=N for the TNT+ values to fall outside the error bars and hence for false positives for DNOC+ to no longer occur. The different E=N dependences shown in Fig. 1 can be understood in terms of the well depths given in Table 1. For TNT+ the value of ε0 corresponds to a thermal energy of 384 K, which means that, when the gas temperature is 473 K, the mobility will decrease monotonically from its zero-field value as E=N increases; this is clearly shown in Fig. 1. For DNOC+, however, the well depth corresponds to a thermal energy of 612 K, so the mobility will increase from its zero-field value until the combination of the thermal and field energies is sufficient to reach the mobility maximum; Fig. 1 shows this behavior and that the mobility maximum occurs for E=N ¼ 85 Td. As a confirmation of these arguments, the calculations were repeated for T ¼ 300 K. The results in Fig. 2 are generally the same as in Fig. 1, but clearly one would have to get to much higher values of E=N to minimize false positives for DNOC+ when the air is at 300 K than when it is at 473 K. Note in
This paper proposes that the incidence of false positives in IMS could be reduced by modifying the instrument so that rapid changes in E=N by a factor of 10–20 could be made. Numerical simulations of the mobility of TNT+ ions in air and of DNOC+ ions in air showed that, when the zero-field mobilities were nearly identical, it was possible to distinguish the ions by such changes in E=N . The main reason this works for these ions is the difference in well depths between the two ion-air interactions. When the well depth is larger than or approximately equal to the thermal energy corresponding to the gas temperature in the experiments, the zero-field mobility is most sensitive to the longrange, attractive portion of the ion-neutral interaction, so the mobility will increase as E=N increases. On the other hand, when the well depth is smaller than the thermal energy, the zero-field mobility reflects the short-range, repulsive portion of the interaction, so the mobility decreases with E=N . It should be noted that it may still be possible to distinguish ions with the same zerofield mobility by making changes in E=N even when both mobilities increase or decrease; all that is needed is sufficiently different slopes on graphs of K 0 vs. E=N . It is recommended that experimenters make the changes in their IMS apparatus needed to investigate the extent to which false positives can be minimized by rapid changes in E=N . It is also recommended that both theorists and experimenters pay attention to the magnitude, rather than simply the position, of the potential minima for the ion-neutral interactions probed by IMS.
1.80 K0
1.75
Acknowledgments This work completes and extends work done by Jacqueline L. Wilmot (nee Regan) under the supervision of the author, as part of an undergraduate thesis at Chatham University in 2009. The author is grateful for her interest in the project and her help in getting it started.
1.70 1.65 1.60
Conflict of interest The author declares that he has no conflicts of interest.
1.55 1.50
References
1.45 E/N
1.40 0
20
40
Fig. 2 Same as Fig. 1, at 300 K
60
80
100
1.
Abedi A, Sattar L, Gharibi M, Viehland LA (2014) Investigation of temperature, electric field and drift-gas composition effects on the mobility of NH4+ ions in He, Ar, N2 and CO2. Int J Mass Spectrom 370:101–106
Author's personal copy Int. J. Ion Mobil. Spec. 2.
3.
4.
5.
6. 7.
8. 9.
10. 11. 12.
13. 14.
Asbury GR, Hill HH Jr (2000) Evaluation of ultrahigh resolution ion mobility spectrometry as an analytical separation device in chromatographic terms. J Microcolumn Sep 12:172–178 Balla G, Koutselos AD (2003) Molecular dynamics simulation of ion transport in moderately dense gases in an electrostatic field. J Chem Phys 119:11374–11380 Benhenni M, Yousfi M, Bekstein A, Eichwald O, Merbahi N (2006) Analysis of ion mobility and diffusion in atmospheric gaseous mixtures from Monte Carlo simulation and macroscopic laws. J Phys D Appl Phys 39:4886–4893 Buchachenko AA, Viehland LA (2014) Mobility of singly-charged lanthanide cations in rare gases: theoretical assessment of the state specificity. J. Chem. Phys. 140:114309 Busch KL (2002) Chemical noise in mass spectrometry. Spectroscopy 17:32–36 ChemBio3D (2008) Perkin Elmer. Waltham, Massachusetts. http:// www.cambridgesoft.com/Ensemble_for_Chemistry/ChemBio3D/. Accessed April 24, 2015 Cottingham K (2003) Ion mobility rediscovered. Anal Chem 75: 435A–439A Crawford CL, Hill HH Jr (2013) Evaluation of false positive responses by mass spectrometry and ion mobility spectrometry for the detection of trace explosives in complex samples. Anal Chim Acta 795:36–43 Eiceman GA, Karpas A (2005) Ion mobility spectrometry, 2nd edn. CRC Press, Boca Raton Eiceman GA, Stone JA (2004) Ion mobility spectrometers in national defense. Anal Chem 76:390A–397A Gharibi M, Viehland LA, Abedi A, Jalili AH, Afsahi G, Behnejad H (2013) Interaction potential and gaseous ion mobility of CO+ ions in He. Mol Phys 111:909–921 Mason EA, McDaniel EW (1988) Transport properties of ions in gases. Wiley, New York Mason EA, O'Hara H, Smith FJ (1972) Mobilities of polyatomic ions in gases: core model. J Phys B 5:169–176
15.
Matz LM, Tornatore PS, Hill HH Jr (2001) Evaluation of suspected interferents for TNT detection by ion mobility spectrometry. Talanta 54:171–179 16. Mesleh MF, Hunter JM, Shvartsburg AA, Schatz GC, Jarrold MF (1996) Structural information from ion mobility measurements: effects of the long-range potential. J Phys Chem 100:16082–16086 17. Shvartsburg AA, Jarrold MF (1996) An exact hard spheres scattering model for the mobilities of polyatomic ions. Chem Phys Lett 261:86–91 18. SIMION® 8.1 Ion and Electron Optics Simulation Software (2015) Scientific Instrument Services. Ringoes, New Jersey. http://simion. com/info/simion81.html. Accessed June 3, 2015 19. Stewart JJP (2008) Stewart Computational Chemistry. Colorado Springs, Colorado. http://openmopac.net/MOPAC2007.html. Accessed April 24, 2015 20. Tabrizchi M, Rouholahnejad F (2005) Comparing the effect of pressure and temperature on ion mobilities. J Phys D Appl Phys 38:857–862 21. Viehland LA, Buchachenko AA (2014) Test of the interaction potential energy for Na+-H2 by gaseous ion transport data. J Chem Phys 141:114305 22. Viehland LA, Chang Y (2010) Transport cross sections for collisions between particles. Comput Phys Commun 181:1687–1696 23. Viehland LA, Chang Y (2012) Beyond the Monchick-Mason approximation: the mobility of Li+ ions in H2. Mol Phys 110:259–266 24. Viehland LA, Dickinson AS (1995) Transport of diatomic ions in atomic gases. Chem Phys 193:255–286 25. Viehland LA, Dickinson AS, Maclagan RGAR (1996) Transport coefficients for NO+ ions in helium gas: a test of the NO+ − He interaction potential. Chem Phys 211:1–15 26. Yousef A, Shrestha S, Viehland LA, Lee EPF, Gray BR, Ayles VL, Wright TG, Breckenridge WH (2007) Interaction potentials and transport properties of coinage metal cations in rare gases. J Chem Phys 127:154309