1. Technical Held
Embodiments of the subject matter disclosed herein relate generally to apparatuses, methods and systems and, more particularly, to devices, processes, mechanisms and techniques for measuring size distribution and/or concentration of particles.
2. Description of Related Art
The ability to measure and to quantify the type, behavior, and/or characteristics of particles (e.g., measurement of individual size and local number and/or mass concentrations) is of utmost importance in a large number of applications of interest, including, for example, medical diagnostics, food preparation, environmental sciences in general, and the control of industrial and vehicular emissions, to name just a few. More specifically, in the medical diagnostics area, for example, it is well established that a risk for a particular disease (e.g., Coronary Heart Disease, or CHD) can be reasonably assessed by measuring mobility spectra (i.e., the mass distribution as a function of particle size) of lipoprotein particles in the blood of a patient. There are also many other medical diagnostics that are based on the ability to measure protein sizes and concentrations, including the ability to identify a genetic disease based on the variation of size and concentration of hemoglobin. In addition, in the area of aerosol and nanotechnology sciences, as applied to environmental concerns, for example, there exist several requirements to measure the size and/or type of a particle that can be related to emissions controls, environmental measurements of how much soot or particulate matter a person may be exposed to in ambient airflows, how much dust may be exhausting into the atmosphere from industrial source emissions, or the quantification of size and concentration of particulate matter in the atmosphere and their effect on the earth's climate. The need to perform particle-related measurement is further evidenced by the large number, and ever more stringent, governmental regulations dealing with the subject.
A variety of measurement techniques have been used to determine the size of particles smaller than 1 μm in diameter. Optical particle detectors determine particle size from the amount of light scattered into a photodetector when a particle passes through a single laser beam. This method is limited to particles larger than about 0.05 μm in diameter. Aerodynamic particle-acceleration lens systems have been designed to accelerate particles according to their size—smaller particles attaining higher velocity than larger particles when a gas containing such particles is accelerated through a nozzle. When the accelerated particles are directed to pass through two precisely spaced laser beams, light scattered by the transiting particle produces two pulses on a photodetector, thus revealing the time it takes a particles to pass through the two laser beams, from which particle velocity can be deduced. This technique is also limited to particles larger than about 0.05 μm in diameter. Aerosol impactors can be used to deposit selected sizes of particles onto a collection surface or into a collection fluid. Chemical analysis and the weight of the collected particles can then be used to construct particle size distributions.
In applications involving macromolecules (e.g., DMA, RNA, and proteins, including their fragments as well as small particles less than 100 nm in diameter), several conventional sizing techniques and/or devices are known, namely, gel electrophoresis, Differential Mobility Analyzers (or DMA), drift tubes, and mass spectrometers, although the latter is used to determine the size of molecules, it actually provides only molecular mass.
Gel electrophoresis is used in clinical chemistry to separate proteins by charge and or size and in biochemistry and molecular biology to separate a mixed population of DNA and RNA fragments by length, to estimate the size of DNA and RNA fragments, or to separate proteins by charge. Nucleic acid molecules are separated by applying an electric field to move the negatively charged molecules through a gel matrix. Shorter molecules move faster and migrate farther than longer ones because shorter molecules migrate more easily through the pores of the gel. Gel electrophoresis can also be used for separation of nanoparticles. Those of ordinary skill will however recognize that gel electrophoresis is a technique that requires a substantial amount of time for completion of any given measurement. Additionally, the position of a band in a gel electrophoresis lane needs to be compared to size standards, typically molecules of known molecular weight, in order to estimate the molecular weight of the material in the band. The need to calibrate gel lanes adds to the effort involved and makes the method a relative measurement technique.
Differential electrical mobility analyzers may be used to determine the size distribution of particles smaller than a micrometer in diameter, in this method, a cloud of charged aerosol particles is drawn between two electrodes, such as the annular space between two concentric cylinders. Voltage applied to the cylinders deflects particles of a predictable size into a particle detector. By scanning the voltage applied to the cylinders, the size distribution of the particles is obtained. U.S. Pat. No. 6,230,572 (which is incorporated herein by reference in its entirety) discloses an example of such an apparatus.
It has been shown (see, for example, U.S. Pat. Nos. 7,259,018, 7,851,224, and 7,713,744, the entire contents of which are incorporated herein by reference) that electrical mobility spectra of lipoprotein particles isolated from human serum reveal simultaneously the size distributions of HDL, LDL, IDL and VLDL particles in a plasma sample, thus revealing a useful technique for assessing risk of CHD. The capability of mobility measurements to span the HDL and LDL lipoprotein size range in one spectrum may be one of the advantages of the use of mobility technology. Conventional gel electrophoresis systems require two different types of gels to reveal the size distributions of HDL particles separately (Agar gel) from LDL, IDL and VLDL particles (gradient density polyacrylamide gel). An example of a conventional measurement of the mobility spectra of lipoproteins obtained from five patients participating in a cholesterol study is illustrated in
Another conventional device used to make measurements of ions is a drift tube.
Therefore, based at least on the above-noted challenges with conventional devices to measure the concentration and size of particles, it would be advantageous to have improved devices to accomplish the summarized tasks, among others, with increased measurement accuracy (particularly in embodiments operating on first principles without the need for calibration), lower cost of manufacturing and operation, reduction on the time required for measurements, and minimization or elimination of the effect of purge gas velocity on the velocity of the particles being measured, while, in some embodiments of the subject matter disclosed herein, increasing the resolution of such measurements by mathematically deconvolving from the measurements the effect of a spread in arrival times due to diffusion and non-ideal background flow velocity variations.
One or more of the above-summarized needs or others known in the art are addressed by apparatuses, methods, and processes to measure average size and concentration of particles as disclosed herein. These apparatuses include a body defining a flow chamber, an exhaust port disposed on the distal end portion of the body; a particle detector disposed inside of the chamber; a focusing electrode disposed inside of the chamber; an electronic gate disposed downstream of the focusing electrode, the focusing electrode and the electronic gate further dividing the chamber into a jet relaxation region a drift region; a plurality of electrodes disposed on the body in the drift region downstream of the electronic gate; and a timing device configured to measure a time taken by the charged particles to travel from the electronic gate to the detector, wherein a measurement of the size and concentration of the charged particles is based on an output signal from the detector and the time measured by the timing device.
The subject matter disclosed herein also includes methods and processes to measure size and concentration of particles. These methods and processes include generating charged particles with a particle source; introducing these particles into a flow chamber of a drift tube; reducing the velocity of the charged particles while the particles are maintained substantially along a centerline of the chamber; controlling the flow of the particles with an electronic gate; controlling the drift velocity of the particles allowed to pass through the electronic gate; measuring the time taken by the charged particles to travel from the exit of the electronic gate to a particle detector disposed at an end portion of the flow chamber; and measuring a current generated by the charged particles impacting the particle detector, wherein a measurement of size and concentration of the charged particles is based on the current generated by the charged particles and the time taken the charged particles to travel from the electronic gate to the detector.
The accompanying drawings (not drawn to scale), which are incorporated in and constitute a part of the specification, illustrate one or more embodiments and, together with the description, explain these embodiments. In the drawings:
The following description of the exemplary embodiments refers to the accompanying drawings. The same reference numbers in different drawings identify the same or similar elements. The following detailed description does not limit the invention. Instead, the scope of the invention is defined by the appended claims. The following embodiments are discussed, for simplicity, with regard to the terminology and structure of apparatuses, systems, or methods for measuring size and concentration of charged particles. However, the embodiments to be discussed next are not limited to these sets, but may be applied to other apparatuses, systems, or methods, including, but not limited to, the measurement of aerosol particle diameter, the measurement of the size distribution of clouds of aerosol particles, the determination of the concentration of particles suspended in a gas, the determination of size and concentration of particles in a liquid by transferring them out of the liquid and into a gas and then analyzing them, and the determination of additional properties of the particles by first separating particles according to their size before they are presented to ancillary analytical methods, such as mass spectrometry, laser induced fluorescence spectrometry and light scatter. In other embodiments the subject matter disclosed extends ion mobility spectrometry and ion drift tubes to the analysis of particles, which cannot be accomplished with current ion mobility spectrometers or ion drift tubes.
As used herein throughout, the term “particle” includes all forms of minute quantities of materials smaller than one micrometer in diameter and suspended in a gas or liquid. The materials may include pure substances, such as any molecule or clusters of identical molecules characterized with a molecular weight. The materials may also include mixtures of materials, such as lipoprotein particles, soot particles, paint droplets, or smoke particles. The term charged-particle refers to any minute form of material that carries an electric charge such as a lipoprotein particle with excess positive or negative charges. Examples of particles include, but are not limited to, particles isolated from a human or an animal serum, particles released from a combustion process, smoke particles, dust particles, nanoparticles, assemblies of molecules, clusters of protein and/or lipid molecules, atmospheric aerosols, nanolipid particles, nanolipid disks, and dumps of such particles. Particles also include charged ions having a molecular weight less than 3×1011 Daltons, equivalent to a 1 micrometer diameter particle.
Reference throughout the specification to “one embodiment” or “an embodiment” means that a particular feature, structure, or characteristic described in connection with an embodiment is included in at least one embodiment of the subject matter disclosed. Thus, the appearance of the phrases “in one embodiment” or “in an embodiment” in various places throughout the specification is not necessarily referring to the same embodiment. Further, the particular features, structures, or characteristics may be combined in any suitable manner in one or more embodiments.
An exemplary embodiment of a drift tube 30 according to the subject matter disclosed is illustrated in
After the jet relaxation region 38, the particles flow through a focusing electrode 42 configured to minimize and/or eliminate the tendency of the jet flow to cause the particles to diverge from the centerline of the drift region 48 of the tube 30, thereby guiding more particles onto a target area, as further explained below, for the purpose of improving detection efficiency. Although the embodiment illustrated in
By substantially reducing the pressure from the atmospheric pressure region 32 to the jet relaxation region 38, the challenge of measuring particles having low drift velocities (for example, particles in the lipoprotein size range) is eliminated and/or substantially ameliorated because aerodynamic drag is reduced, thus allowing particles to travel faster for a given electric field strength than the background gas velocity, it is known that drift velocity scales inversely with pressure, and to a first approximation, gas flow velocity in the drift region also scales inversely with pressure. One of the advantageous aspects of operation at reduced pressure is that the jet relaxation region 38 provides a method of transition and transport of charged particles into a uniform low velocity gas, along with acceptable drift times for the particles (for example, drift times on the order of 400 ms in some embodiments). Embodiments similar to the one illustrated in
In one embodiment, an electrospray source is used for generating charged particles suspended in a gas. Particles flow through the nozzle 36 (which, in one embodiment, comprises a small orifice having a diameter in the range of about 10 to about 1000 microns or preferably from about 200 to about 600 microns) at sonic velocity, where the downstream pressure in the jet relaxation region 38 is any value less than 0.5 bar so as to assure choked flow conditions at the nozzle 36, i.e., the gas jet and particles have sonic velocities at the point of minimum area of the nozzle 36. Downstream of the nozzle 36, the particle-laden jet 40 expands, with axial velocities decreasing to values of about 1 m/s within a short distance (e.g., less than 20 cm axial distance), which is approximately the location of the focusing electrode 42 illustrated in
Placement of the electrodes for the electronic gate 44 (or gate electrodes) downstream of the focusing electrode 42 in a region where the axial jet velocity is relaxed to lower values, as explained, allows the charged particles flowing near the centerline of the particle-laden jet 40 to be stopped by the voltage applied to the electronic gate 44. In some embodiments, the electronic gate 44 is held at this voltage for a period of time ranging from about 0.01 to about 10 seconds, so that the flow of charged particles can be collected at this gate position, increasing the ultimate signal levels of charged particles arriving at the detector 48, following their release from the region near the electronic gate 44.
Those of ordinary skill in the applicable arts will appreciate that, during the time particles are accumulated in the electronic gate 44, particle diffusion will tend to disperse the charged particles radially, although the focusing electrode 42 around the region where the electronic gate 44 is located will help to confine incoming charged particles around the centerline of the jet. After sufficient charged particles have been collected, the voltage of the gate electrode 44 and the focusing electrode 42 is raised above the voltage needed to stop charged particles to initiate the drift of all particles towards the detector 48. By proper choice of the various physical characteristic of the drift tube 30 (such as overall length and diameter), all particles in the range from about 5 to about 50 nm will have drift velocities greater than or equal to the background gas flow velocity, for example, about 30 cm/s, and travel nearly straight towards the detector
Those of ordinary skill in the applicable arts will also appreciate, after consideration of the subject matter disclosed herein, that the use of a nozzle and the subsequent jet expansion into an area of low pressure provides several benefits as compared to conventional designs. First, the particle laden jet transitions to the low pressure drift tube environment and the gas-particle interactions confine the particles near to the tube centerline, allowing for the particles to drift in a near-uniform electric field in the central core of the drift tube 30. Secondly, the expanding jet provides a way to reduce gas velocity, so that the gas velocity is small compared to the ultimate particle drift velocity in the drift region of the instrument. Thirdly, the diameter of the orifice may be chosen (e.g., values in a range from about 10 to about 1000 microns, or, preferably, from a range from about 200 to about 600 microns) to achieve sufficient flow of charged particles to give adequate signal at the current sensing detector 48, yet not too large a flow that increases the background gas flow velocity in the drift region 46. Indeed the final diameter of the drift tube (e.g., a value around 7 cm for some embodiments) may be selected so as to maintain a low background flow velocity). Fourthly, the exit pressure of the drift tube may be defined by a vacuum pump, which can readily achieve pressures below 0.01 atm in some embodiments, giving higher drift velocities and thus shorter analysis times. Vacuum pump operating pressures in the range from about 0.1 to about 0.01 atm may be used in some applications. Fifthly, the overall length of the instrument (jet relaxation region 38 plus the drift region 46) may be chosen to be about 35 cm and the flow tube diameter to be about 7 cm, thus resulting in transition of particles across the jet relaxation region in about 25 ms, while drift times may range in some embodiments from about 5 ms for 5 nm singly-charged particles to about 200 ms for 50 nm singly-charged particles at an operating pressure of about 0.03 bar, a drift tube length of about 9 cm, and an electric field of about 1 kV/cm.
As charged particles arrive at the detector, they are drawn into the porous metal material where they are collected by filtration and deposit their charge which passes to a current measuring electronic circuit. The particles have an arrival time distribution caused by their size distribution. Smaller particles arrive at the detector before larger ones. The charge delivered to the detector by the particles generates a varying current in the electronic circuit and the amplitude of the current as a function of arrival time can be used to determine particle size and particle concentration. The drift tube is designed to measure predominantly singly charged particles so that the magnitude of the defector current at any specific drift time is proportional to the number of particles of a specific size. Those familiar with the art of time-of-flight measurements will appreciate the mathematical relationship between particle charge, particle diameter, electric field intensity, particle velocity, and gas pressure, along with background gas velocity. The resulting detector signal is a time-varying current from which particle size distributions can be mathematically extracted.
Therefore, the jet relaxation section 38, the gate region 54, and the drift region 46 make up the drift tube 50. In operation, gas, a small number of neutral particles, and an appropriate number of charged particles are drawn into the drift tube via an inlet orifice by virtue of the reduced pressure inside of the drift tube 50. The in-flowing particle-laden gas results in the formation of an expanding jet. Gas conditions in the jet, i.e., gas velocity and pressure, are such as to carry the particle-laden gas towards the gate region 54. Initially, the particle-laden gas jet is confined closely to the centerline of the drift tube 50. After the gas travels a certain distance (e.g., about 30 cm for some embodiments), the jet relaxes and expands across the bore of the drift tube 50. Placement of the gate region 54 at a distance from the inlet orifice is selected so as to allow the charged particles to be stopped by the voltage applied to one of the grids 56, 58, 60 mounted perpendicularly to the longitudinal axis of the drift tube 50. If the grids 56, 58, 60 are mounted too close to the inlet orifice, the gas velocity will be higher and will push particles through the grids even when high voltage is applied to them. If the grids are mounted farther downstream, the jet may have expanded to the wall of the drift tube 50 and a fraction of the particles will become lost to the wall or fail to reach the detector. In some embodiments, the grids 56, 58, 60 are screens extending across the entire cross-sectional area of the tube at the axial location where each is disposed. In some embodiments, only one grid is used. In another, two grids are used. Also, although three grids are illustrated in
In embodiments using three grids at the gate region 54, the drift tube 50 may be operated in three different modes, i.e.: (1) a blocking mode configured to prevent charged particles from entering the drift region 46; (2) a trapping mode, configured to accumulate charged particles in the region around the gate region 54; and (3) an injection mode configured to inject a burst of charged particles into the drift region 46.
The blocking operating mode is configured to prevent most particles from entering the drift region 46 for a period of time during which the background gas sweeps particles already in the drift region 46 towards the exit. This action improves the response of the detector 72 by eliminating or minimizing background signals. In this mode, the grids 58 and 60 are powered at blocking potentials while the grid 56 shields incoming particles from the electric field produced by the grid 58.
In the trapping operating mode, a stopping potential is applied to grid 60 and a lower trapping potential is applied to grid 58, while grid 56 shields the incoming charged particles from the potential applied to grid 58. The trapping condition creates an electrostatic potential valley perpendicular to the bore of the drift tube 50 in which charged particles are trapped. This condition provides a way to accumulate a substantial number of charged particles in the gate region 54 so that when they are eventually injected into the drift section 46, a larger detector signal is obtained, thus increasing the resolution and accuracy of the final measurement while reducing the time required for each measurement.
The injection mode starts by switching the potential applied to grid 58 and focusing electrode 62 during the accumulation mode to values higher than that applied to the grid 60. This raises the trapped particles to a potential that allows them to roll down a “potential hill” in the electric field inside the drift region 46 in the direction of the detector. It also blocks charged particles from entering the drift region 46 while the injected particles are rolling down the “potential hill.”
In some embodiments, electric fields are established in the drift tube by means of ring electrodes and ring electrodes equipped with a fine mesh grid. The ring electrodes in the drift region are fabricated from stainless steel. A fine mesh grid having an open area of 85 percent is attached to a supporting ring electrode. The combination of a ring electrode and a fine mesh grid is a grid electrode. Practitioners in the art of fabricating mass spectrometers will appreciate the care needed to prevent electrical discharge from high voltage electrodes that are operated in a reduced pressure environment. The edges of the ring electrodes were polished to remove burrs and thus minimize the possibility of electrostatic discharge from burrs. The grid electrodes and ring electrodes are separated from each other in the drift region and trapping region by means of dielectric spacers. In some embodiments, electrical contact between a high voltage power supply and each ring electrode or grid electrode was accomplished by pressing a 1/16″ diameter metal rod through a plastic tube fitting (⅛″ NPT pipe thread× 1/16″ compression fitting) that was mounted externally on a radius of the insulating cylinder that surrounded the drift region. The bore of the tube fitting was aligned with the outer edge of a ring electrode. Electrical contact between the rod and a ring electrode was secured by friction and the opposite end of the rod was attached to a high voltage power supply by means of a corona ball. The corona ball was machined with internally threaded passages that accommodated set screws and served as a mechanical fixture for connecting the rod to a wire leading to the power supply and furthermore provided a way to minimize corona discharge by covering the sharp edges of the wire and the rod.
A representation of the axial variation of the relative electric potentials in a drift tube according to one embodiment of the subject matter disclosed for the three modes of operation of the electronic gate disposed in the electrode gate region 54, as just explained, is presented in
In the drift region 46 charged particles are exposed to a fixed electric field generated by maintaining an approximately constant voltage drop between each electrode, thus the electric field is constant along the bore of the drift region. This means that regardless of where particles are located radially or longitudinally, they will experience a well-defined electric force that moves positively charge particles towards the detector. The profile of the gas velocity in this section is nearly flat, but not perfectly flat, and this imperfection contributes to a variation in particle velocity across the diameter of the drift region. Knowledge of the background gas flow field, as obtained from CFD calculations, allows this variation to be remedied by means of mathematical algorithms, i.e., deconvolution of the arrival time distribution signal, as further explained below.
The detector is preceded by a grounded grid (the detector grid 68 shown in
In addition, as previously explained, when voltage is applied to the focusing electrode 62, charged particles are focused towards the center of the bore of the drift tube 50 as shown by the particle trajectory 64b in
An up-potential hill is depicted in
Another exemplary embodiment of the subject matter disclosed is illustrated in the drift tube 80 of
In the embodiment shown in
In the embodiment illustrated, one or more valves 89 may also be provided so as to control the background gas flow during operation of the drift tube 80. In operation, after charged particles are accumulated in the gate region, the valve 89 is dosed so as to momentarly stop the flow of background gas in order to minimize the effect of that variable on the measurement. Once the measurement is finished, the valve 89 is opened and the gas flow starts again. In another embodiment, a second valve 91 is installed at the inlet of the drift tube 80, as shown in
After evaluation of the subject matter disclosed herein, those of ordinary skill in the applicable arts will appreciate the existence of many variations and equivalent embodiments and associated methods and systems to be discussed further later. For example, the background gas could be introduced ‘backwards’ through the porous detector plate and/or through the secondary porous wall surrounding the detector. Such embodiments would be advantageous at least in part due to optimized operating conditions for small molecules when solvent vapor might degrade performance. In other embodiments, the electrode gate could be designed to focus charged particles more severely before they are released by a 2-stage voltage switching protocol, thus such a focusing could then be designed to minimize the influence of counter-current gas flow and prevent charged particles from being blown backwards up the trapping region. In other embodiments, the counter-flow gas could be a reactive gas that is used to identify specific charged particles due to a shift in drift time that corresponds to an adducted species of charged particles. In yet other embodiments, elimination of the reactive charged particle pulse due to the co-transport of negatively charged particles being electrostatically dragged along by positively charged particles would be advantageous. Also, operation in negative polarity modes would be desirable in some embodiments because some materials such as DNA are more readily charged with a negative charge. Finally, the embodiments explained here to prepare and deliver charged particles to the drift region so as to allow the measurement of their size and concentration could also be used for other purposes. For example, the first chamber containing the relaxation region and the electrode gate may be used to hold charged particles for secondary purposes, such as to generate an e-beam to alter charge states; to illuminate the accumulated charged particles with a free electron laser (FEL) for coherent x-ray diffraction; to measure the fluorescence properties of the trapped particles; or to hold them for release into a mass spectrometer.
It has been reported in the literature that some molecules change shape (conformation) in concordance with the level of charge they carry. It is useful in some studies to alter the charge state of trapped molecules for the purpose of measuring their cross section as a function of the number of charges they carry. This information may be used to investigate the stability of individual molecules or the stability of molecular clusters. Fluorescence measurements of trapped molecules could be used for example as a way to measure the efficiency of chemical derivatization reactions that are commonly used to label protein molecules. Coherent x-ray diffraction provides 3-D information about a molecule's structure. The present way to perform coherent x-ray diffraction studies is to introduce a beam of nanocrystals or molecular complexes into a rapidly switched x-ray beam and collect diffracted photons before the x-rays explode the target molecule. Current technology consumes inordinate amounts of carefully isolated molecules. By trapping the molecules in the drift tube and releasing them synchronously with the pulsed x-ray laser it will be feasible to use samples more efficiently. Some forms of ion mobility have already been operated in combination with a mass spectrometer. The development of a charged particle mobility spectrometer, as described here, will further enable the study of molecular conformation with molecular identification by means of mass spectrometry.
Those of ordinary skill in the applicable arts will appreciate that the subject matter disclosed can be further broadened and/or optimized by the use of a mathematical model of the applicable physical concepts, such as fluid flow, particle interactions with the flow of a gas, selection of operational voltage levels for the various electrodes disclosed and explained, etc. In the paragraphs that follow, a summary is presented of the modeling analysis performed by the inventors to further investigate the many benefits of the subject matter disclosed by use of a combination of simple analytical relationships that guided the general concepts, followed by detailed flow calculations using a CFD commercially available software. The applicable analytical equations provided general functional dependencies on the important parameters, while the CFD software identified the detailed influence of turbulence, flow recirculation, gas velocity and pressure within a drift tube. Furthermore, by use of user-defined functions, the CFD software was further modified to include the simulation of an electrical field and particle dispersion and/or behavior in order to simulate charged-particle motion, gas and particle trajectories, and travel time caused by the gas flow and electric fields.
Based on the values of the Reynolds numbers, or Red, applicable to the drift tubes disclosed herein, the CFD software was configured to simulate turbulent flow conditions based on values of tube diameter, jet exit velocities, and kinematic gas viscosities. Two different turbulence models were tested. Velocity contour plots generated using k-ε and S-A turbulence models gave similar results, indicating that the particular choice of turbulence model was not critical. Because the k-ε model is the most general viscous model used in fluid dynamics and simulates viscous effects better than the S-A model, the calculations for particle trajectory carried out in the CFD code were all based on the k-ε model.
Initial model calculations performed direct flow calculations without consideration for particle effects. Following confirmation of a feasible jet flow, particles were then seeded into the gas flow calculations with an electric field and the resulting motion of singly charged particles was modeled. Given the required number of charged particles for feasible measurements, particle concentrations are low, and thus do not affect the gas-phase flow. In addition, the particles are small and follow the primary flow field because equilibration times with the flow for 20 nm particles at 0.1 atm are less than a microsecond.
Small particle slip was modeled considering values of the Knudsen numbers, or Kn, which characterizes the relationship between gas phase and particle flow based on particle diameter, mean free path of gas phase molecules, and pressure. For small particles and low pressures, Kn is large, on the order of 100 or greater, which suggests the use of a correction factor (Cunningham, Cc), when particle trajectories are calculated. For the small particles modeled in the drift tube, they generally follow the primary gas jet, although they have sufficient momentum to avoid recirculation as discussed below. This can be shown to be true by consideration of a characteristic stopping distance, i.e., the distance that the particles would traverse if injected into a static environment with an initial velocity Vo. Conversely, the stopping distance can also be viewed as a “starting” distance to accelerate a particle to a final drift velocity. It has been shown that about five stopping distances are required for a particle to equilibrate within 1% of the gas velocity. For small particles at low pressures, characteristic stopping times are small (<1 μs), and thus the particles (in the absence of an electric field) tend to follow the flow very closely. For example, the corresponding stopping distance for Vo equal to 10 m/s is only about 1.4 μm. Thus particles equilibrate with the gas phase flow very quickly. There is a small velocity lag at the initial highest velocities, but micro-second lag times are small compared to total residence times in the jet of 25 ms. Particle drift velocities that may be present when a particle is accelerated in an electric field can also be modeled by taking into consideration a terminal steady state velocity of the particles in a given electric field, resulting in an equation for a total drift time that depends only on the background gas velocity, drift velocity, and the length of the drift tube. Finally, the effects of diffusion during the drift period can lead to a wider packet of particles, reducing time resolution. In addition, during the gate collection period, diffusion can contribute to radial spread if electric field confinement is not used. This effect has been taken into account by considering a RMS diffusion distance over the drift tube length, resulting in a conclusion that the irreducible diffusion broadening of a drift tube is independent of particle size, and only dependent on the electric field and drift tube length; variables that are relevant in defining the resolution limits of the instrument. In practice, the ratio of drift velocities to background flow velocity is more important in determining the resolution for the largest particles. Singly charged smaller molecules travel faster in the drift tube and therefore have a shorter time to diffuse, while larger ones travel more slowly, but since their diffusion velocity is also smaller, the two end up diffusing about the same distance as they fly through a fixed-length drift tube.
In order to evaluate the extent to which the gas flow influences particle trajectories the dispersion of one hundred 20-nm-diameter particles was simulated by computationally releasing them across the jet orifice and using the same expanding jet conditions shown in
The simulated particle trajectories show that it is feasible to install an electrostatic gate, for stopping and releasing particles, for the purpose of making time-of-flight particle electrical mobility measurements. The particle trajectory calculations show that there is an optimal gate location, where the velocity across the jet radius is nearly uniform, yet the particle jet radius has not over-expanded into lower particle velocity regions.
Following a gate accumulation period from about 0.5 to about 10 seconds (to increase signal levels), but not limited to this time range, particles are released into the uniform electric field drift region, where particles become separated in time and space according to their diameter. This behavior is simulated in the CFD model, again by gating one hundred 20-nm-diameter singly charged particles into the expanding jet at a location where the local flow velocity, in absence of an electric field, is approximately 1 m/s. A uniform electric field in the range of 1 kV/cm is also introduced into the model. Because the particles are confined closely to the centerline, it is assumed that a practical drift-tube electric field will be nearly uniform along the longitudinal-axis. The results of this calculation now represent the behavior of charged particles in the flow field combined with the electric field drift velocities, i.e., the performance of a classical drift tube.
Using the results of the CFD computations, one can compute the variation in TOF to the detector by two different methods. The first method is to compute directly, using the CFD code, the variation in TOF for each particle size. A second simpler and more flexible method is to use the variation in axial velocity profiles (as a function of radius), to compute the variation in TOF by computing the sum of background flow velocity and drift velocity as a function of size. The latter method allows decoupling of the TOF computation from the CFD calculations, once an optimal geometry for the jet flow has been determined, using two axial velocity profiles at the centerline and maximum radius to determine the maximum variation in TOF.
One of the conventional methods of characterizing drift tube performance is to optionally specify a minimum size resolution based on the maximum observed variation in TOF. In practice, the resolution based on standard deviation of the time-of-flight signal may be a factor of two better, so this minimum size resolution is a conservative estimate of performance.
One parameter of choice of drift tube voltages, discussed in more detail below in the section on electrostatics, is the Paschen limit for corona discharge, which must be avoided. The Paschen limit predicts an upper limit for establishing an electric field without electrical discharge for a range of operating pressures. The selection of voltage parameters for the gate and drift region are chosen so as not to exceed the Paschen limit. In general, there appears to be improvement in resolution by operating at the highest possible voltage for a given pressure.
As understood by those of ordinary skill, there may exist additional fundamental and practical sources of resolution degradation. One fundamental limit is based on diffusion during the drift time. Simulations carried out for an electric field of 1000V/cm, and a drift tube length of 10 cm indicated that diffusion appears to be negligible for particles in the lipoprotein size range. In addition, there is a limitation given by the physical width of the potential valley in which particles are trapped, relative to the total length of the drift tube, which provides a significantly better estimate of resolution compared to the more deleterious effects of variation in background flow velocity. An overall resolution of the system may be given by incorporating the effects of diffusion, gate spread (width of trapped particle cloud), and jet background velocity variations. In this case, comparisons with current differential mobility methods give a maximum resolution value (1 sigma) value in the range of 25, decreasing with particle size and pressure, in sharp contrast to the subject matter disclosed.
In order to simulate the electric charges at the drift tube gate and electrodes commercially available software, Simion, along with a modified Statistical Diffusion Simulation (SDS) Model was used to evaluate a number of different designs for the gate electrodes and the electrodes in the drift region. Simion calculates the shape of the electric field around electrodes and the SDS Model calculates the motion of charged particles as they diffuse, or are dragged, through a reduced pressure gas under the influence of an electric field superimposed by Simion. The results of these calculations are particle trajectories and particle velocity at any location along a particle trajectory.
One way to present the result of the Simion and SDS calculations is to plot particle trajectories superimposed on a potential energy surface. A simple way to think of the potential energy surface is to consider the path of a golf ball on a putting green. A charged particle responds to the shape of the potential energy surface in the drift tube similarly to the way a golf ball moves on a putting green—it goes downhill and moves away from hills.
In
As shown, the particle trajectory 120 deviates from the centerline of the jet due to the influence of the expanding jet and arrives at the first portion of the gate electrode 116 in the form of a fine mesh metal grid in some embodiments. After passing through this first portion of the gate electrode 116, the particle experiences the effect of a focusing electrode 117 that introduces a small potential ridge that turns the particle towards the center of the drift tube. The focusing electrode 117 serves to confine particles (not shown) near the center of the drift tube, as already explained. The particle then rolls down a potential hill where it is stopped by an abruptly rising potential upslope. The upslope defines a trapping condition that catches incoming particles and prevents them from travelling axially. Here particles are confined in a potential groove in which they can move only radially (as a result of diffusion), yet are also constrained radially by the influence of the focusing electrode's potential. The shape of the valley in the trapping region is another advantageous design feature in the operation of the drift tube. The trapping valley provides a way to accumulate particles exiting the nozzle and holds them in a groove that serves to line them up along a starting line. The sharpness of the potential groove, i.e., the narrowness of the starting line, enhances resolution of a TOF measurement because after release all particles begin their drift towards the detector at nearly the same position. The trapping valley can be disposed in several embodiments, including, for example (1) as described; (2) with two trapping grids for further focusing the particles on the starting line; (3) with one or two curved grids; or (4) a gridless design wherein the particles are trapped in an up-sloping potential valley.
After review of the subject matter disclosed herein, it will become apparent to one skilled in the applicable arts that many other embodiments with different shapes and numbers of grids and screens are within the scope of this disclosure.
In
In
In
In
Once an optimal drift tube configuration is designed, it remains possible to improve resolution by incorporating knowledge of the physical mechanisms that degrade performance, such as multiple charging, diffusion and gate spread, and flow velocity variations. If the spread mechanisms are approximately constant and known, deconvolution methods can be used to correct for non-ideal effects. For example, if by measurement calibration or calculation the shape of the actual arrival time distribution for a discrete particle size is determined, it will be possible to post-process full length TOF spectra and correct for the arrival time spread, as explained next.
The ideal current at a given time is specified by the number of particles of size and charge related to the TOF. However, because of the resolution factors described above, there is a spread in arrival times due to diffusion and non-ideal background flow velocity variations. However, once a flow geometry has been defined, the resulting TOF spread function for each size can be measured and/or modeled. Thus the current measured at each time is a convolution of different particle sizes and charge levels.
For specified instrument operating conditions, the relationship between particle drift velocity, Vd and particle size, dp is given by:
where / is the drift tube length, and td is the drift time. The drift time in turn is related to the measured particle velocity Vp by Vd=Vp−Vg, where Vg(r,x) is the gas velocity experienced by each particle size. For matrix deconvolution, the canonical form of Eq. (1) is used to characterize the correct form of size binning so that
To give a size resolution of 100 or greater, the time bins should be scaled to give Δtd/td<0.02. This measurement accuracy can readily be achieved.
Using the above time and size scaling conditions, one can then write a recursive set of equations that relates the charged particle number count per unit time, Cn (or measured charged particle current) at a measured time interval following charged particle release from the electrode gate to the flow of various charged particle sizes in the drift tube. The charged particle flow Cn is given by the convolution of charged particles of different sizes traversing different trajectories at varying gas velocities, yet arriving in the same time bin, Cn(1/tn) with time width defined in Eq. (2). In the following, the number concentration density, Nn(1/dn2), also has a bin width defined by Eq. (2). In the following, the charged particle concentration Nn for all sizes is assumed to be uniform. If the charged particle are not uniformly distributed in the grid trap prior to release, the following equations can readily be modified by a weighting function (of radius) that accounts for the radial distribution of charged particles and trajectories.
The primary equation for arbitrary Cn is written algebraically, following the general method described by Holve and Self (Holve, D., and Self, S. “Optical particle sizing for in situ measurements: Parts 1 and 2, Applied Optics, 18(10), 1632-1652, 1979, the contents of which are incorporated herein in their entirety by reference):
Cn(1/tn)=Qn,nNn+Qn,n+1Nn+1+ . . . +Qn,n+jNn+j (3)
Here, the volumetric flow of charged particles follows a specific trajectory of specified gas and drift velocity in a symmetric area annulus around the axis of the drift tube, i.e.,
where Vnk (cm/s) is the particle velocity corresponding to the sum of the average gas velocity, Vgl(rn) (integrated over the drift tube length) at a defined radius, plus the charged particle drift velocity Vd(dk). The cross-section, ΔSk [cm2] is dependent only on the radius, and not on the charged particle size (=2πrn Δrn).
Thus smaller charged particles (with faster drift velocities) traveling on slower gas velocity trajectories at larger radius can have the same total particle velocity and arrive at time interval tn. The number of terms, j, in the sum depends on the overall range of possible gas velocity trajectories that hit the detector ranging from r=0 to rd. This also determines the size range of faster charged particles that can arrive at the time interval tn.
For a specified set of drift tube instrument parameters (pressure, voltage, and geometry), the values of Qn,k can be determined through CFD calculations for each particle size. Once these values are specified, the system of equations represented in Eq. (3) can be readily solved by back substitution. These equations range from C1 and d1, which is the slowest and largest particle of interest (longest transit time), to the fastest charged particles of interest, Cm and dm. Back substitution begins with the fastest charged particles that have little or no overlap in transit time with other charged particles so that the initial conditions for starting the solution is just 1 term, i.e.:
Cn(1/tm)=Qm,mNm. (5)
As mentioned above, if the charged particles are not uniformly distributed in the trap, which is likely to be the case for a particular embodiment, a weighting function wn(rn) can be incorporated to multiply each Qn,k by wn (Σwn≦1) so as to account for the non-uniform distribution of charged particles with radius. This weighting function can be determined by using the results of CFD calculations that determine the radial gas velocities, coupled with a charged particle flow simulation program (e.g., Simion) that calculates the radial distribution of charged particle sizes at the electrode gate prior to release into the drift tube. Thus if the population of particles is confined to a smaller radial annulus, the variation of gas velocity and charged particle sizes that can contribute to a measured Cn is reduced, thereby increasing the initial size resolution prior to deconvolution. Note that the charged particle distribution sum is less than or equal to 1, accounting for possible radial loss of charged particles.
Furthermore, the above-summarized equations can also be modified to include multiple charging effects if they are significant. The drift velocities are proportional to the square root of the charge, which in the above analysis is assumed to be 1. However, if there are multiple charged particles that potentially overlap slower drift particles, these terms can be added to the set of equation represented by Eq. (3), correcting the final result for faster particle overlap on larger particle sizes. The probabilities of multiple charging have been calculated by Fuchs (see, for example, Fuchs, N. A. (1963), On the Stationary Charge Distribution on Aerosol Particles in a Bipolar Ionic Atmosphere, Geofisica Pura e Applicata, 56, 185, the contents of which are incorporated herein by reference in their entirety). The number concentrations in each size bin (Ni) may be readily converted to standard number concentration frequencies, dN/dlogd.
Eq. (3) can be simplified to varying degrees, the simplest being the assumption that Cn is related to one average particle size traversing the drift tube at one average drift and gas velocity. This effectively corresponds to a rolling average solution of all the terms up to index j, with one overall term in each row equation. This is the method currently used to show the results given by current experimental measurements further described below, but those of ordinary skill in the art will understand that this assumption is not intended to limit the scope of the subject matter disclosed herein in any way. As will be explained further below, the effects of particle overlap increase with increasing size, thereby reducing resolution. However, once the instrument geometry is defined and optimal configuration determined, it is straightforward to incorporate all the known parameters, and to solve the set of equations represented by Eq. (3) iteratively and rapidly in real time to obtain significantly improved high resolution results for all particle sizes of interest.
In order to illustrate the effect of the above-summarized advantageous deconvolution, experiments were performed using particles of a known size. It should be understood that these results are not presented here to limit the subject matter disclosed, but simply to illustrate the contents of this disclosure.
Given the variation in electrospray conditions, some variation in results is expected. Repeat measurements for fixed instrument parameters are very precise, within 1% or so of the indicated size. Peak Raw signal levels for the experimental measurement results shown in
The results above are based on 512 repeat scans, each of which is 0.6 seconds long, or a total of 5 minutes. The electrode gate accumulates and stops particle flow for 0.6 seconds (a variable) and then releases the particles into a stable electric and background gas flow field that transports the small particles faster than the large ones, giving a signal (proportional to the number of charged particles at each size) that is a function of time. This repeat scan process improves the S/N by sqrt(512), or by a factor of 22 over a single 0.6 second scan. With an optimized A/D design and processing, the computational overhead will be negligible, and spectra will display in real time. The fact that the signal/noise keeps improving with increased number of scans shows how stable the supersonic jet and gas flow is in the system, along with the switching high voltage system of the electrode gate. For these exemplary measurements, the output of the detector integrating amplifier signal is cabled to an oscilloscope, where the end result is downloaded into an Excel file.
For the exemplary measurements just shown, the computation comprises the following steps. First, rolling average smoothing of integrator output data following a numerical differentiation of the smoothed data. Secondly, summing the integrator output and the differentiated signal (times feedback constant=0.5 sec), leading to the input charged particle voltage as a function of time (10,000 data points from the particular data acquisition system used in these exemplary measurements). Subsequently, converting the time scale to size, dp=C* Sqrt (E/PVd), where C is a constant (including viscosity), E is electric field (volts/cm), P is the DT pressure, and Vd is the drift given the difference in particle and gas velocities, i.e., Vp−Vgavg. Vp, is given by l/tm, is the measured particle velocity detected as time in the data acquisition system, tm, and l is the drift tube length. Vgavg(r) is the average gas velocity at each radial trajectory from the electrode gate to the detector (9 cm in this embodiment), and varies with axial position, x, and radius r. If Vd>>Vgavg, then using a constant value is a good approximation. Currently, all the calculations use a constant value, so the size is not as accurate for larger particles above 20-30 nm, but those of ordinary skill will know readily how to relax this assumption based on the subject matter disclosed.
In assessing accuracy and error analysis for the exemplary experimental results just described, the general formula above can be re-arranged to give explicit results in terms of time, i.e.:
In Eq. (6), if the drift velocity is fast relative to the background gas velocity, (i.e., tg>>tm or tg/(tg−tm)≈1), sizing accuracy is just dependent on accuracy of time measurements, tm, E, P, and l. These values can be specified to better than 2%. For large particles the greatest uncertainty is the variation in tg, giving an absolute uncertainty around 8% for 50 nm particles.
In order to obtain maximum resolution a more detailed evaluation of the variation in Vg(x,r) is desired, which is given by the CFD calculations already explained. The radial velocity values in the jet expansion region explain why charged particles (and signal) are lost radially in a conventional electrode gating system. These radial velocities in the region of the electrode gate are on the order of 10 cm/s, and show that the maximum charged particle accumulation times are on the order of 50-100 ms. At longer times they flow to the wall outside the 1.5 cm detector radius, and are not measured. A modified focusing electrode and trapping grid has been designed (using Simion) to constrain all particles within a specified detector diameter, (currently 35 mm), eliminating or substantially reducing radial charged particle loss, and boost signal levels by an estimated factor of ten or more.
The subject matter disclosed herein also includes methods and processes to measure size and concentration of particles.
Based on the subject matter disclosed herein and the exemplary method disclosed in
One or more of the steps of the methods comprising the subject matter disclosed may be implemented in a computing system specifically configured to calculate the size and concentration of charged particles as explained hereinabove. An example of a representative computing system capable of carrying out operations in accordance with the exemplary embodiments is illustrated in
The exemplary computing system 1200 suitable for performing the activities described in the exemplary embodiments may include server 1201. Such a server 1201 may include a central processor (CPU) 1202 coupled to a random access memory (RAM) 1204 and to a read-only memory (ROM) 1206. The ROM 1206 may also be other types of storage media to store programs, such as programmable ROM (PROM), erasable PROM (EPROM), etc. The processor 1202 may communicate with other internal and external components through input/output (I/O) circuitry 1208 and bussing 1210, to provide control signals and the like. The processor 1202 carries out a variety of functions as is known in the art, as dictated by software and/or firmware instructions.
The server 1201 may also include one or more data storage devices, including a hard drive 1212, CD-ROM drives 1214, and other hardware capable of reading and/or storing information such as DVD, etc. In one embodiment, software for carrying out the above discussed steps may be stored and distributed on a CD-ROM 1216, removable memory device 1218 or other form of media capable of portably storing information. These storage media may be inserted into, and read by, devices such as the CD-ROM drive 1214, the disk drive 1212, etc. The server 1201 may be coupled to a display 1220, which may be any type of known display or presentation screen, such as LCD displays, LED displays, plasma display, cathode ray tubes (CRT), etc. A user input interface 1222 is provided, including one or more user interface mechanisms such as a mouse, keyboard, microphone, touch pad, touch screen, voice-recognition system, etc.
The server 1201 may be coupled to other computing devices, such as the landline and/or wireless terminals via a network. The server may be part of a larger network configuration as in a global area network (GAN) such as the Internet 1228, which allows ultimate connection to the various landline and/or mobile client devices. The computing device may be implemented on a vehicle that performs a land seismic survey.
Exemplary embodiments of software comprising the subject matter disclosed include those to analyze charged particle samples automatically. The software may be configured control a pipetting robot to dispense samples into 96 well plates, from which aliquots are removed and treated with standard chemical methodologies to isolate lipoproteins from serum, or whatever might be analyzed by this method. Further automated pipetting steps, controlled by software, may deliver purified lipoproteins to a fresh 96 well plate that will be stored for later analysis by drift tube mobility spectrometry. The stored samples may be further processed to deliver them to an electrospray source that pumps the lipoprotein sample at a controlled flow rate for the purpose of generating gas-borne lipoproteins that are drawn into the drift tube.
Furthermore, these software products may record the output of the drift tube detector, apply various algorithms to clean the data, such as smoothing, and may then mathematically process the output signal to generate a secondary signal that can be further processed mathematically to produce charged particle mobility spectra. The spectra may be treated with deconvolution algorithms to generate information useful to classify lipoproteins according to the way cardiologists use particle size analysis to diagnose coronary heart disease and establish medical treatment regimens for these or other applicable health conditions.
Software products are also within the scope of the subject matter disclosed that will generate a plurality of reportable summaries. A first example is a software to interface into the laboratory information management system of analytical labs. A second example is a software appropriate for communicating test results to the prescribing physician and to the patient.
The disclosed exemplary embodiments provide apparatuses, methods, and systems for measuring size and concentration of charged particles as well the other uses hereinabove summarized. It should be understood that this description is not intended to limit the invention. On the contrary, the exemplary embodiments are intended to cover alternatives, modifications and equivalents, which are included in the spirit and scope of the invention as defined by the appended claims. Further, in the detailed description of the exemplary embodiments, numerous specific details are set forth in order to provide a comprehensive understanding of the claimed invention. However, one skilled in the art would understand that various embodiments may be practiced without such specific details.
Although the features and elements of the present exemplary embodiments are described in the embodiments in particular combinations, each feature or element can be used alone without the other features and elements of the embodiments or in various combinations with or without other features and elements disclosed herein.
This written description uses examples of the subject matter disclosed to enable any person skilled in the art to practice the same, including making and using any devices or systems and performing any incorporated methods. The patentable scope of the subject matter is defined by the claims, and may include other examples that occur to those skilled in the art. Such other examples are intended to be within the scope of the claims.
While the disclosed embodiments of the subject matter described herein have been shown in the drawings and fully described above with particularity and detail in connection with several exemplary embodiments, it will be apparent to those of ordinary skill in the art that many modifications, changes, and omissions are possible without materially departing from the novel teachings, the principles and concepts set forth herein, and advantages of the subject matter recited in the appended claims. Hence, the proper scope of the disclosed innovations should be determined only by the broadest interpretation of the appended claims so as to encompass all such modifications, changes, and omissions. In addition, the order or sequence of any process or method steps may be varied or re-sequenced according to alternative embodiments. Finally, in the claims, any means-plus-function clause is intended to cover the structures described herein as performing the recited function and not only structural equivalents, but also equivalent structures.
This application claims priority to U.S. patent application of Ser. No. 61/781,734, filed on Mar. 14, 2013, entitled “Apparatuses, Processes, and Systems for Measuring Particle Size Distribution and Concentration,” the contents of which are herein incorporated by reference in their entirety.
At least part of the materials disclosed herein was made with Government support under Contract No. 1 R43 HL080933-01. As such, the Government may have certain rights therein as provided for by the terms of these contracts.
Number | Name | Date | Kind |
---|---|---|---|
3845301 | Wernlund | Oct 1974 | A |
4777363 | Eiceman | Oct 1988 | A |
5719392 | Franzen | Feb 1998 | A |
5747799 | Franzen | May 1998 | A |
5905258 | Clemmer | May 1999 | A |
5936242 | De La Mora | Aug 1999 | A |
5949001 | Willeke | Sep 1999 | A |
5969349 | Budovich | Oct 1999 | A |
6012343 | Boulaud | Jan 2000 | A |
6051832 | Bradshaw | Apr 2000 | A |
6229143 | Wernlund | May 2001 | B1 |
6230572 | Pui et al. | May 2001 | B1 |
6263744 | Russell | Jul 2001 | B1 |
6323482 | Clemmer | Nov 2001 | B1 |
6389912 | Wood | May 2002 | B1 |
6408704 | Willeke | Jun 2002 | B1 |
6412359 | Zborowski | Jul 2002 | B1 |
6418802 | Wood | Jul 2002 | B1 |
6491872 | Wick | Dec 2002 | B1 |
6498342 | Clemmer | Dec 2002 | B1 |
6583407 | Fischer | Jun 2003 | B1 |
6586732 | Lee | Jul 2003 | B2 |
6593567 | Abdel-Rahman | Jul 2003 | B1 |
6601464 | Downing | Aug 2003 | B1 |
6748815 | Povey | Jun 2004 | B2 |
6796195 | Povey | Sep 2004 | B2 |
6809314 | Yamada | Oct 2004 | B2 |
6829955 | Mahgerefteh | Dec 2004 | B1 |
6862536 | Rosin | Mar 2005 | B2 |
6919367 | Bacon | Jul 2005 | B2 |
6931950 | Malachowski | Aug 2005 | B2 |
6959618 | Larsen | Nov 2005 | B1 |
7066037 | Keskinen | Jun 2006 | B2 |
7140266 | Marjamaki | Nov 2006 | B2 |
7164122 | Fuhrer | Jan 2007 | B2 |
7213475 | Coghill | May 2007 | B2 |
7213476 | Cheng | May 2007 | B2 |
7244931 | Zimmermann | Jul 2007 | B2 |
7259018 | Benner et al. | Aug 2007 | B2 |
7268347 | Blanchard | Sep 2007 | B1 |
7298127 | Golder | Nov 2007 | B2 |
7713744 | Benner et al. | May 2010 | B2 |
7851224 | Benner et al. | Dec 2010 | B2 |
20010032930 | Gillig | Oct 2001 | A1 |
20020011560 | Sheehan | Jan 2002 | A1 |
Number | Date | Country | |
---|---|---|---|
20140260702 A1 | Sep 2014 | US |
Number | Date | Country | |
---|---|---|---|
61781734 | Mar 2013 | US |