Ailin
Li
,
Sandilya V. B.
Garimella
and
Yehia M.
Ibrahim
*
Biological Sciences Division, Pacific Northwest National Laboratory, Richland, WA 99352, USA. E-mail: Yehia.ibrahim@pnnl.gov
First published on 30th October 2019
Probing molecular properties in the gas phase requires the integration of complementary ion manipulation approaches such as ion mobility spectrometry. Structures for lossless ion manipulations (SLIM) have recently been developed to perform ultra-high resolution ion mobility separations using traveling waves as well as providing other advanced capabilities. Despite its success, the design aspects of SLIM have not been fully explored and remained largely unchanged. Here, we report on a computational study using SIMION simulations of a number of traveling wave (TW) patterns that can be used in SLIM. The TW pattern used in the current SLIM device is a set of 8 electrodes where, at any time, 4 electrodes are held at high voltage (i.e., 1111), while the other 4 electrodes are held at low voltage (i.e., 0000), forming one micro-trapping region of 11110000 pattern. Ion trajectory simulations demonstrated the feasibility to simplify the 8-electrode set to a shorter pattern (e.g., 6-electrode or 4-electrode set) while maintaining or improving the performance. The RF and TW amplitudes, guard voltage, and TW speed were optimized subsequently on the symmetric patterns of the 4-, 6-, and 8-electrode sets to further improve the performance. The resolution, peak broadening, peak capacity, and peak generation rate of each pattern were evaluated, showing that the 111000 pattern of the 6-electrode set has comparable performance to the current 11110000 pattern and is always better than the 1100 pattern. This work provides insight into the feasibility for simplification and modification of the TW configuration in SLIM and other traveling wave devices.
, where q is the ion charge, E is the applied constant electric field, L is the length of the drift tube, kB is the Boltzmann constant, and T is the temperature. The performance of DTIMS is limited by the maximum electric field that can be applied (typically a few kV) and geometry constraints imposed by the linear drift tube length (more than 2 m length becomes cumbersome and impractical). This effectively limits the maximal achievable resolving power from DTIMS instruments. Unlike DTIMS, which requires a fixed high voltage, TWIMS uses a low amplitude traveling wave (TW) to create a dynamic field to drive and separate ions. Ions with higher mobility can keep up with the TW than lower mobility ions that tend to fall behind by the passing wave. This form of ion separation enables the use of a fixed voltage without the debilitating limitation of high voltage drops associated with DTIMS. When the TW approach is coupled to a flexible ion path design (e.g., serpentine path), enabled by Structures for Lossless Ion Manipulations (SLIM), there is potential for ultra-long path length separations. SLIM15 uses two planar surfaces patterned with arrays of electrodes that are supplied with RF, DC, and TW potentials to confine and drive ions in gases at a pressure of a few Torr. The electrode arrays can be patterned on printed circuit boards (PCBs) at a low cost and are easily fabricated. The ability to pattern such electrodes on surfaces enables complex ion manipulations such as turning around a corner,16 switching ions to the orthogonal path,17 and moving ions between levels using elevators and escalators.18 These manipulations enable extremely long path length and ultra-high resolution ion mobility separations through multi-pass operation.15,19 The serpentine path design and multi-pass operation have been demonstrated for ultra-high resolution with more than 500 m ion path length, indicating SLIM's ability to achieve high resolution IM separation.15
SLIM instruments built thus far have involved high resolution instrumentation targeted for biological analysis. The nominal path length was 13 m with the ability for multi-pass. The TW patterns applied used 8 electrodes with wave profiles that were 45° phase-shifted to create a TW based IM separation. However, the performance of TW created with different numbers of electrode sets to create a TW was not investigated. One of the most effective ways to influence the performance of SLIM is to optimize the TW configuration. The TW configuration is defined as the number of TW electrodes in a set, while the TW pattern is the number of electrodes held at high and low voltages. The current SLIM design has 8 electrodes forming a single TW configuration which is repeated throughout the entire path length. In the 8-electrode set, 4 consecutive electrodes are held at high potential, while the remaining 4 electrodes are held at low potential (i.e., a pattern of 11110000, where 1 represents an electrode with high potential and 0 an electrode with low potential). The TW is created when the pattern is stepped through the 8 electrodes one electrode at a time. However, this configuration requires a power supply with 8 separate channels. If the TW configuration can be reduced to less than 8 electrodes, power supply requirements can be reduced, SLIM can be simplified, and tighter corners can be designed. These features enable the construction of a longer path length SLIM in a smaller footprint and provide the basis for a small and portable SLIM device. Previously published work from our group compared the mobility resolution for different patterns of the 8-electrode set and indicated that the current pattern of 11110000 is optimum.20 Other work has also demonstrated that a 4 repeat pattern of the 8 TW configuration has 20% increased field strength in the Synapt G2 IM-MS instrument compared to the 6 repeat pattern of the 8 TW configuration in the Synapt G1 IM-MS.3,21
The investigation of different traveling wave configurations and patterns would enable designing simpler waveforms with similar or better peak characteristics. For instance, the nature of the serpentine path separation in SLIM requires ions to turn around corners. While straight motion of ions entails some diffusional broadening, turning around corners has to be accomplished with minimization of “race track” effects which could lead to a reduced separation performance. This is typically overcome in SLIM by making sharp 90° turns. This is accomplished by intersecting ions traversing a path with an orthogonally moving wave. Precise alignment of the TW electrodes is required to maintain wave continuity at the turn, which is simpler if there is flexibility in the TW configuration.
The efforts in SLIM development so-far have not examined a more diverse range of TW electrodes, simpler patterns, or smaller number of electrodes. Here, we perform an ion simulation study to explore the feasibility of simplifying the TW configuration and its effect on the peak properties. We demonstrate that other TW configurations can be used without the loss of resolution via optimizing the operating parameters. The resolution trends on different TW patterns are calculated to investigate the impact of the TW patterns. Other operating parameters including RF and TW amplitude, guard voltage, and TW speed of 6- and 4-electrode sets were optimized. The resolution of different TW patterns has been demonstrated to be enhanced and comparable to the currently used 11110000 pattern. This work provides insight into the feasibility of simplifying and modifying the TW configuration in SLIM and other traveling wave devices while maintaining or improving the performance.
The mean arrival time and its standard deviation for the ions that transverse the full path were recorded. The resolution of m/z 622 and 922 ions was calculated using eqn (1):
![]() | (1) |
![]() | ||
| Fig. 3 Resolution of each pattern of different TW electrode sets (3 replicate runs for each). (a) 4-Electrode set; (b) 6-electrode set; (c) 8-electrode set; (d) 12-electrode set. | ||
As shown in Fig. 3, the ion mobility resolutions of the symmetric of 6-, 8-, and 12-electrode sets are close, which is promising for using a lower number of TW electrodes. The resolution of the 1100 pattern is ∼20% lower than that of the other symmetric patterns which can be further optimized and improved. However, the 11100000 and 11111000 patterns show similar (within error bars) resolutions compared to the 11110000 pattern. The error bars were calculated from the results of three replicates. The variations in the simulation arise from the computation of collisional statistics by the SDS model, wherein the variance in the arrival time distribution decreases with increasing ion populations at the cost of increasing computational time.
The RF potentials are utilized in SLIM to confine ions between the two surfaces. When the RF amplitude is not sufficiently high, ions are not effectively confined, especially ions with high m/z. The TW is used to separate and move ions along the axial direction. When the TW amplitude is too high, microtraps are formed where ions are pushed forward by the TW voltage “wall”, causing ions to move at the same speed as the TW, termed “surfing”. However, if the TW amplitude is too low, the ion drift velocity is low and the moving wave frequently overtakes the ions, resulting in more rollover events. This causes the effective flight time of ions within the SLIM device to increase, leading to increased peak broadening.
The RF and TW amplitudes (for TW patterns of 1100, 111000 and 1111000) were optimized by calculating the resolution for a range of RF and TW amplitudes. The RF amplitude was ramped from 250 Vp–p (in 50 V steps) to reach a practical voltage limit of 400 Vp–p (onset of electrical breakdown in our current setup), while the TW amplitudes of 10, 20, 25, 30, 35, 40, and 50 V were used. A total of 100 ions of each of m/z 622 and 922 were utilized in a path length of 152 mm. A longer path length of 1320 mm (as used in other simulations here) and a larger number of ions were unrealistic in this case owing to an extremely long computational time and the large number of simulations conducted. Therefore, while the trends of computed resolutions should be similar, the absolute values are lower compared to the results in Fig. 3. As shown in Fig. 4a–c, the optimal RF amplitude is ∼400 Vp–p for all the TW patterns, and the optimal TW amplitudes for 1100, 111000, and 11110000 are 40 V, 30 V, and 25 V, respectively. The resolution increased slightly for all the TW patterns with the increase in the RF amplitude. A higher RF amplitude results in a better confinement and reduced electric field variation causing a slightly higher resolution. This was demonstrated by comparing the electric field of ion plumes at two different RF amplitudes of 250 Vp–p and 400 Vp–p with the same TW amplitude (30 V) using the 111000 TW pattern. Fig. 4d shows a cross-sectional view of SLIM where a tighter ion plume is observed at a higher RF amplitude. Fig. 4e shows the TW electric field distribution that an ion experiences moving in the axial-x dimension (i.e., the direction of ion motion). The results indicate that at higher RF amplitudes, ions experience a tighter TW electric field distribution, which provides less deviation in ion speed and less peak broadening, allowing higher resolution to be achieved than with a lower RF amplitude.
After the critical parameters have been optimized, guard voltage was optimized to investigate its impact on the resolution. Guard voltage is a low DC voltage that is used to constrain the ions along the two sides of the ion path in the lateral direction (i.e., perpendicular to ion motion). If the guard voltage is too low, ions are not efficiently constrained, leading to ion leaking from the two sides. If the guard voltage is too high, there may be a field penetration from the guard electrode sides to the electric field at the ion trapping center, causing a field distortion in the trapping region.
An RF amplitude of 400 Vp–p was used for all of the TW patterns, and the optimized TW amplitudes of 40 V, 30 V, and 25 V were applied for 1100, 111000, and 11110000 TW patterns, respectively. Fig. 5a shows the resolution trends of m/z 622 and 922 ions at the different guard voltages. The optimal guard voltages are 15 V, 17.5 V and 20 V for the 1100, 111000 and 11110000 TW patterns, respectively. The resolution was not calculated at lower voltages since the guard voltage was too low to prevent ion losses. The resolution trends indicate a slight impact of guard voltage on the resolution change, which is somewhat negligible in the resolution optimization.
To better understand the small effect of guard voltage, the ion plume distribution was plotted at 4 different guard voltages (12.5 V, 17.5 V, 22.5 V, and 27.5 V) using the same RF amplitude (400 Vp–p), TW amplitude (30 V), and 111000 pattern. Fig. 5b shows the ion plume distribution in the yz plane, presenting a more laterally constrained ion plume at higher guard voltages. Fig. 5c shows similar TW electric field profiles in the axial-x dimension at different guard voltages, which confirms the negligible effect of guard voltage on the resolution. Similar results were obtained at an RF amplitude of 250 Vp–p (ESI Fig. 1†) but with a wider ion plume distribution due to the weaker RF field. Overall, these results indicate that the guard voltage does not significantly affect the RF or TW electric field but slightly affects the ion confinement in the lateral direction.
TW speed was optimized as another parameter to improve the IM resolution for each TW pattern at their optimal RF and TW amplitudes. As shown in Fig. 6, the observed trend of the TW speed on the resolution is similar to previous experimental findings and theoretical predictions.15,20,23 When the TW speed is too low (or ions’ mobility is sufficiently high), ions will surf and no separation will occur. When the TW speed is increased, m/z 622 ions with higher mobility move at the same velocity as the traveling wave without rolling over, while m/z 922 ions with lower mobility roll over and begin to separate. The optimal resolution for a specific m/z is obtained at the TW speed where ion motion transitions from the surfing to separation mode. When the TW speed is increased further, the resolution shows a decreasing trend to a plateau. In this high TW speed range, rolling over commonly occurs, which does not significantly change the resolution as shown in Fig. 6. The optimal TW speeds of 1100, 111000 and 11110000 TW patterns are 82 m s−1, 88 m s−1 and 82 m s−1, respectively.
![]() | ||
| Fig. 7 Resolution comparison of 1100, 111000 and 11110000 TW patterns before and after the optimization. | ||
The improvement of the resolution can also be inferred from understanding the spatial distribution of ions across the TW bins. Ions experiencing a large number of rolling over events tend to spread across a large number of TW bins and vice versa. To only understand the impact of rolling over events, we fixed the number of TW bins to 14 for 1100, 111000, and 11110000 patterns. 1000 ions were used in this simulation. Shown in Fig. 8a–f is the spatial distribution of ions after traversing the 14 TW bins before and after the optimization. It is clear that the spatial distribution of ions for all three patterns studied is much narrower after optimization.
![]() | ||
| Fig. 8 Ion distribution histograms recorded for unoptimized (top) and optimized (bottom) wave patterns of (a) 1100, (b) 111000 and (c) 11110000. | ||
Device performance is commonly evaluated by the resolution, which gives information about how well two charged species can be separated according to their mobility. Other metrics such as peak capacity are also useful in the performance evaluation for TWIMS. Peak capacity is a dimensionless measurement representing the maximum number of peaks that can be resolved in a single separation. Peak capacity can be calculated using eqn (2):24,25
![]() | (2) |
| TW pattern | Peak capacity (unoptimized) | Peak capacity (optimized) | Improvement |
|---|---|---|---|
| 1100 | 30 | 30 | 0 |
| 111000 | 46 | 43 | −6% |
| 11110000 | 37 | 45 | 22% |
The peak generation rate is another figure of merit calculated as the ratio of the peak capacity to the range of the separation time, shown in eqn (3):20,26
![]() | (3) |
We estimated the peak generation rate of each TW pattern before and after the optimization (Table 2). Before the optimization, the measured peak generation rates of 1100, 111000, and 11110000 TW patterns were 258 s−1, 402 s−1, and 527 s−1, respectively. The 11110000 pattern shows a higher peak generation rate than the other patterns, while after optimizing the parameters, the peak generation rate of 1100 and 111000 has been significantly improved by 120–200%, showing a much higher improvement than the 11% of the 11110000 pattern. The highest peak generation rate among the three TW patterns is 910 s−1 for the 111000 pattern, which is 55% higher than that of the 11110000 pattern and approximately 3 orders of magnitude greater than those in condensed phases (e.g., the average peak generation rate of liquid chromatography is ∼0.1–1.2 s−1) due to the fast speeds of gas-phase ions.25,27 However, ion mobility separation is a post-ionization technique while liquid chromatography is a pre-ionization technique and their integration is beneficial for maximum separation power. This result shows that a higher peak generation rate can be achieved from the 111000 pattern than the 1100 and 11110000 patterns after the optimization. This increase in the peak generation rate will maximize the separation power of SLIM when integrated with other separation approaches such as liquid chromatography.
| TW pattern | Peak generation rate (before) | Peak generation rate (after) | Improvement |
|---|---|---|---|
| 1100 | 258 s−1 | 775 s−1 | 200% |
| 111000 | 402 s−1 | 910 s−1 | 126% |
| 11110000 | 527 s−1 | 587 s−1 | 11% |
Footnote |
| † Electronic supplementary information (ESI) available. See DOI: 10.1039/c9an01509d |
| This journal is © The Royal Society of Chemistry 2020 |