Treatment of ion peak areas in mass spectrometry.

Abstract

To provide a method of analyzing a signal generated by a mass analyzer.SOLUTION: A method disclosed herein comprises receiving a signal generated by a mass analyzer, determining area of a first ion peak of one or more ion peaks in the signal, and estimating the number of ions that contributed to the first ion peak. The number of ions that contributed to the first ion peak is estimated by determining a correction to be applied to the area of the first ion peak from a correction function, and applying the correction to the area of the first ion peak. The correction function describes a relationship between average single ion area and ion mass, mass-to-charge ratio and / or charge for the mass analyzer.SELECTED DRAWING: Figure 13

Description

[Technical Field]

[0001] The present invention relates to methods for analyzing ions, and in particular to mass spectrometry (MS) and mass analyzers. [Background technology]

[0002] Time-of-flight (ToF) mass analyzers exploit the property that the transit time of an ion in an electrostatic field is proportional to the square root of the ion's mass-to-charge ratio (m / z). Ions are emitted from an ion source, accelerated to a desired energy, and travel a specific distance before striking an ion detector. The signal produced by this detector is recorded, typically resulting in time-resolved peaks, each produced by one or more ions with the same mass-to-charge ratio (m / z).

[0003] Because the travel distance is essentially the same for all ions, the arrival time of the ions is used to determine the mass-to-charge ratio (m / z) of the ions, which can then be used for identification. In addition, the area of ​​a peak can be used to estimate the number of ions that contributed to that peak, which can in turn be used for quantification.

[0004] It is believed that there remains room for improvement in devices and methods for mass spectrometry. Summary of the Invention

[0005] A first aspect provides a method of analysing a signal produced by a mass analyser, the method comprising: receiving a signal produced by a mass analyzer, the signal including one or more ion peaks; determining an area of ​​a first ion peak of the one or more ion peaks; Estimating the number of ions that contributed to the first ion peak by (i) determining a correction to be applied to the area of ​​the first ion peak from a correction function, the correction function describing the relationship between the average single ion area and the ion mass (m), mass-to-charge ratio (m / z), and / or charge (z) for the mass analyzer, and (ii) applying the correction to the area of ​​the first ion peak.

[0006] An embodiment provides a method for analyzing a signal generated by a mass analyzer.

[0007] The mass analyzer may be a time-of-flight (ToF) mass analyzer. Accordingly, the analyzer may include an ion detector disposed at the end of the ion path. Packets of ions may be injected into the ion path, and once injected, the ions may travel along the ion path to the detector for detection. The detector may be configured to produce a signal indicative of the ion intensity received at the detector as a function of time (arrival).

[0008] The mass analyser may alternatively be an ion trap mass analyser or a quadrupole mass analyser. Thus the analyser may comprise either an RF ion trap configured to eject ions in order of mass to charge ratio (m / z), whereupon the ions pass to a detector for detection, or a quadrupole mass filter configured to filter ions according to their m / z, whereby ions transmitted by the mass filter pass to a detector for detection. The detector may be configured to produce a signal indicative of the ion intensity received at the detector as a function of time.

[0009] In this method, one or more ion peaks are identified in the signal, and then the area of ​​one or more, or each, of the identified ion peaks is determined, for example, by summing or integrating the area of ​​the signal under the peak. Then, for one or more, or each, ion peak, the number of ions that arrive at the detector and create the peak is estimated based on the area of ​​the ion peak by applying an appropriate correction to the ion peak area, for example, by dividing (or multiplying) the ion peak area by a correction factor.

[0010] In this method, the correction to be applied to the area of ​​a particular ion peak is determined from a correction function for the mass analyzer. This correction function describes the relationship between the average single ion area (SIA) and the ion mass, mass-to-charge ratio (m / z), and / or charge for the analyzer. In some embodiments, the correction function describes only the relationship between the average single ion area (SIA) and the ion mass or mass-to-charge ratio (m / z) for the analyzer. In other embodiments, the correction function also describes (e.g., separately) the relationship between the single ion area (SIA) and charge for the analyzer. For example, the correction function can provide an approximate average area of ​​an ion peak produced by the analyzer in response to detecting a single ion as a function of the ion's mass or m / z (or equivalently, as a function of the ion's velocity or arrival time), and optionally, as a function of the ion's charge.

[0011] A correction for a particular ion peak can be determined from the correction function by determining the mass or m / z (or equivalently, velocity or arrival time), and optionally charge, associated with the ion peak, and then determining the value of the correction function at that determined mass or m / z (or velocity or arrival time), and optionally charge. The value so determined can be applied to the ion peak area (e.g., by dividing (or multiplying) the ion peak area by that value) to estimate the number of ions that contributed to the ion peak. Thus, the correction can correct for the mass and / or charge dependence of the ion area generated by the analyzer.

[0012] In this regard, although the area of ​​an ion peak produced by a mass analyzer, such as a ToF analyzer, typically has a pronounced dependence on the mass and / or charge of the ion, known mass analyzers do not systematically correct for this mass and / or charge dependence. Furthermore, there are no known global correction functions that provide appropriate correction values ​​across the entire range of mass and / or charge associated with mass analyzers, such as time-of-flight mass analyzers, that can facilitate such correction.

[0013] As described in more detail below, embodiments provide a correction function that describes the relationship between the average single ion area and the ion mass m / z and / or charge over the entire operable mass m / z and / or charge range of the mass analyzer, thereby allowing systematic correction of the dependence of ion area on mass and / or charge in a particularly accurate and straightforward manner.

[0014] This method is particularly advantageous for time-of-flight mass analyzers in which ion packets are accumulated in the ion trap before being injected from the ion trap into the drift region of the analyzer. In such analyzers, it is often necessary to precisely control the total number of ions in each ion packet accumulated in the ion trap, optimizing the number of ions to be below, but as close as possible to, a limit of the ion trap, such as the space charge limit of the ion trap. This is commonly done using so-called automatic gain control (AGC) methods, which can control the number of ions in subsequent ion packets based on an accurate measurement or estimate of the number of ions in the ion packet. Each ion packet can contain ions spread over a wide range of masses and / or charges (and therefore have a wide range of corresponding SIAs), and thus correcting for the dependence of ion area on mass and / or charge can provide a significant improvement to the accuracy of the estimation of the number of ions in the ion packet and to the AGC method.

[0015] Embodiments also provide improved quantification of MS data, for example label-free quantification of data between ToF-MS.

[0016] It will therefore be appreciated that embodiments provide improved methods of analyzing data generated by an ion analyzer.

[0017] The mass analyzer can form part of an analytical instrument, such as a mass spectrometer. The instrument can include an ion source. Ions can be generated from a sample in the ion source. The ions can be passed from the ion source to the analyzer through one or more ion optical devices disposed between the ion source and the analyzer.

[0018] The one or more ion optical devices may comprise any suitable arrangement of one or more ion guides, one or more lenses, one or more gates, etc. The one or more ion optical devices may include one or more transmit ion guides for transmitting ions, and / or one or more mass selectors or filters for mass selecting ions, and / or one or more ion-cooled ion guides for cooling ions, and / or one or more collision or reaction cells for segmenting or reacting ions, etc. The one or more or each ion guide may include an RF ion guide, such as a multipole ion guide (e.g., a quadrupole ion guide, a hexapole ion guide, etc.), a segmented multipole ion guide, a stacked ring ion guide, etc.

[0019] The mass analyzer may be any suitable mass analyzer, such as a time-of-flight (ToF) mass analyzer, an ion trap mass analyzer, or a quadrupole mass analyzer.

[0020] In embodiments, the mass analyzer is a ToF analyzer configured to determine the mass-to-charge ratio (m / z) of ions from their arrival times. Typically, a ToF analyzer may include an ion injector located at the beginning of an ion path and an ion detector located at the end of the ion path. The analyzer may be configured to analyze ions by determining their arrival times at the detector (i.e., the time it takes for an ion to travel from the injector, through the ion path, and arrive at the detector).

[0021] The ion path may have any suitable form, such as being linear in the case of a linear ToF analyzer, or including one or more reflections in the case of a ToF analyzer comprising a reflectron or a multi-reflection time-of-flight (MR-ToF) analyzer.

[0022] In certain embodiments, the analyzer is a multi-reflection time-of-flight (MR-ToF) analyzer, such as a tilted mirror multi-reflection time-of-flight mass analyzer of the type described in U.S. Pat. No. 9,136,101, or a single focus lens multi-reflection time-of-flight mass analyzer of the type described in British Patent No. 2,580,089.

[0023] The analyzer may thus include two ion mirrors spaced apart and facing each other in a first direction X, each mirror generally elongated along a drift direction Y between a first end and a second end, the drift direction Y being orthogonal to the first direction X. The ion injector may be disposed near the first end of the ion mirrors and configured to inject ions into the space between the ion mirrors. The detector may be disposed near the first end of the ion mirrors and configured to detect the ions after they have completed multiple reflections between the ion mirrors. The analyzer may be configured to analyze ions by injecting ions from the ion injector into the space between the ion mirrors, whereupon the ions (a) drift from the deflector along the drift direction Y toward the second end of the ion mirror, (b) reverse their drift direction velocity near the second end of the ion mirror, and (c) take a zigzag ion path with multiple reflections in direction X between the ion mirrors while drifting along the drift direction Y back to the deflector for detection.

[0024] The ion implanter can have any suitable form, such as, for example, one or more (e.g., orthogonal) accelerating electrodes. However, in certain embodiments, the ion implanter comprises an ion trap. The ion trap can be configured to receive ions (from an ion source via one or more ion optical devices) and accumulate ion packets, for example, by accumulating the ions during an accumulation period. The ion trap can be configured to inject each accumulated packet ion into an ion path (e.g., by accelerating the ion packet along an ion path), whereupon the packet ions can travel along the ion path to a detector. The ion trap can have any suitable form, such as being an extraction trap, for example, as described in GB Patent No. 2,580,089.

[0025] The mass analyser may alternatively be an ion trap mass analyser or a quadrupole mass analyser. Thus the analyser may comprise either an RF ion trap configured to eject ions in order of mass to charge ratio (m / z), whereupon the ions pass to a detector for detection, or a quadrupole mass filter configured to filter ions according to their m / z, whereby ions transmitted by the mass filter pass to a detector for detection. The detector may be configured to produce a signal indicative of the ion intensity received at the detector as a function of time.

[0026] In various embodiments, the detector may be any suitable ion detector, such as one or more conversion dynodes, optionally followed by one or more electron multipliers, one or more scintillators, and / or one or more photon multipliers, etc. The detector may be configured to detect ions received at the detector and may be configured to produce a signal indicative of the intensity of the ions received at the detector as a function of time (arrival).

[0027] The detector may include a digitizer, such as a time-to-digital converter (TDC) or an analogue-to-digital converter (ADC), configured to digitise each signal to produce a digitised signal. Thus, each signal produced by the analyzer may be a digitised signal. The digitiser may have a single channel or multiple channels, e.g., the signal from the detector may be split between a first high-gain channel and a second low-gain channel (to increase the dynamic range).

[0028] Each signal may include one or more ion peaks, such as multiple ion peaks, each peak being produced by an ion (or multiple ions having substantially the same mass-to-charge ratio (m / z)) detected by the detector.

[0029] In some embodiments, the ToF analyzer can be configured to generate a signal for each ion packet. For example, multiple such ion packets can be sequentially injected into the ToF ion path and detected by a detector. Thus, the method can include repeatedly (i) injecting ion packets into the ToF ion path, (ii) detecting the ion packets at the detector, and (iii) producing respective signals for the ion packets. The ion packets can be injected into the ToF ion path (and signals produced) at any desired rate, such as a rate greater than about 100 Hz and less than about 10 kHz, e.g., about 200 Hz.

[0030] Each signal received and processed in an embodiment may be generated from signals relating to a single ion packet, or each signal received and processed may be generated from signals relating to multiple ion packets (i.e., multiple ion injections), whereby multiple (digitized) signals are combined (e.g., averaged).

[0031] In the method, one or more ion peaks may be identified within the signal. This may be done in any suitable manner, for example, using any suitable peak detection algorithm. The method may proceed by determining one or more characteristics of the one or more identified ion peaks or of each identified ion peak. Suitable characteristics may include, for example, the centroid and / or intensity of the ion peak. This may be done in any suitable manner, for example, by fitting a suitable peak model to the identified ion peaks. Any suitable peak model may be used, for example, a Gaussian or an asymmetric Gaussian.

[0032] One or more characteristics of each ion peak can be used as needed. For example, physicochemical properties associated with each ion peak, such as its mass, charge, and / or mass-to-charge ratio (m / z), can be determined using one or more characteristics. The physicochemical properties associated with an ion peak can be determined using only one or more characteristics associated with that peak, or by using both one or more characteristics associated with that peak and one or more characteristics associated with one or more other peaks. For example, the m / z associated with an ion peak can be determined from the arrival time (i.e., center of gravity) associated with that peak. The charge associated with an ion peak can be determined from prior knowledge of the sample and / or instrument settings, or by taking into account the associated isotope pattern and / or neighboring charge state ion peaks. The mass associated with an ion peak can be determined from the m / z and charge associated with the peak.

[0033] In addition to these characteristics, the method determines the area of ​​one or more or each ion peak identified in the signal. This may be done in any suitable manner, for example, by summing or integrating the area of ​​the signal under the peak (e.g., relative to a noise threshold) and / or by fitting a suitable ion peak model to the peak and determining the area of ​​the ion peak model.

[0034] The ion peak area so determined is then used to estimate the number of ions that contributed to the ion peak, i.e., how many ions impinged on the ion detector dynode to produce the ion peak. This is done by applying a correction to the ion peak area. This correction may be applied in any suitable manner, for example, by dividing or multiplying the ion peak area by a correction factor, or by performing any equivalent operation.

[0035] If multiple different ion peaks are identified in the signal (corresponding to ions having different masses or m / z), the method can include estimating the number of ions that contributed to each of the ion peaks by applying a different correction to each of the ion peaks. The estimated number of ions that contributed to each ion peak in the signal can then be summed to estimate the total number of ions that contributed to the signal, e.g., the total number of ions in the ion packet that produced the signal.

[0036] In this method, the correction to be applied to the area of ​​a particular ion peak is determined from a correction function for the analyzer. The correction can be determined from the correction function by determining the mass or m / z (or equivalently, velocity or arrival time) associated with the ion peak (i.e., determining the mass or m / z of the ion giving rise to the ion peak) and then determining the value of the correction function at the determined mass or m / z (or velocity or arrival time). The value so determined can be used as the correction factor, for example, by dividing (or multiplying) the ion peak area by the correction factor.

[0037] The correction function can take any suitable form that describes (at least) the relationship between the average single ion area (SIA) and ion mass or m / z for the analyzer. As used herein, "single ion area" or "SIA" is the area of ​​an ion peak generated by the analyzer in response to detecting a single ion. For example, the correction function may be or include a curve of the average SIA as a function of mass or m / z (or equivalently, as a function of ion velocity or arrival time). In certain embodiments, the correction function provides an (approximate) average SIA as a function of ion mass or m / z (or equivalently, as a function of ion velocity or arrival time) for the analyzer. Thus, the number of ions that contributed to a particular ion peak can be estimated by obtaining the average SIA from the correction function at the mass or m / z associated with the ion peak and dividing the measured area of ​​the ion peak by this value.

[0038] In embodiments, the correction function describes the relationship between the average SIA and ion mass or m / z for an analyzer over the m / z range of interest for the analyzer, such as over the entire operable m / z range of the analyzer. As used herein, the "operable m / z range" of an analyzer is the range of ion m / z that the analyzer is designed to analyze. In certain embodiments, the mass analyzer is designed to analyze ions in the life sciences, e.g., small and large organic molecules, biomolecules, DNA, RNA, proteins, peptides, and fragments thereof. Thus, the correction function can describe the relationship between the average SIA and ion m / z over a range from about 25, 50, 75, or 100 to about 6000, 8000, 10,000, 15,000, or more. In some embodiments, the correction function describes the relationship between the average SIA and ion m / z over an m / z range of (at least) about 50 to about 8000.

[0039] The correction function may be continuous over most or all of this m / z range, and the correction function may vary continuously over most or all of this m / z range, which allows for systematic correction of the dependence of ion area on mass or m / z in a particularly accurate and straightforward manner.

[0040] The correction functions can have any suitable form and can be derived in any suitable manner. For example, the correction functions may be derived based on a "first principles" theoretical analysis of the interaction of ions with the detector. However, the inventors have recognized that due to the large mass and / or charge parameter space and the inherent complexity of molecular ions typically analyzed in the life science MS field (e.g., including a wide variety of different higher-order structures), such first principles analysis is not well suited to the present context.

[0041] Thus, in certain embodiments, the correction function is derived by fitting (parameters of) the model to experimentally acquired SIA data for that analyzer. This experimentally acquired SIA data can include experimental measurements of SIA for the analyzer for ions across most or all of the mass m / z and / or charge range of interest. The SIA data can include experimental measurements of SIA for ions of the same (e.g., single) charge. The SIA data can include multiple measurements of SIA at each of multiple different masses or m / z.

[0042] The SIA data to which the model is fitted can be obtained in any suitable manner. For example, it would be possible to obtain SIA data for each and every individual analyzer and derive a correction function for each respective analyzer by fitting (parameters of) a model to the SIA data for that analyzer. However, in certain embodiments, SIA data is obtained for only one (or a few) specific representative analyzers of a class of analyzers, and a global correction function is derived by fitting (parameters of) a model to that data. The global correction function can then be used as, or to derive, a correction function for one or more different analyzers of the same class (as described further below).

[0043] The SIA data to which the model is fitted can be the raw SIA data recorded for a (representative) analyzer, or the raw SIA data can be appropriately processed and / or corrected before the model is fitted to the processed and / or corrected data.

[0044] In embodiments, the data is first cleaned to remove effects such as noise and / or peak ringing. This may be done by truncating the raw SIA data points for a given mass using a lower and / or upper threshold. In these embodiments, the lower and / or upper threshold may be determined in any suitable manner. For example, the lower threshold may be a fixed threshold, and the upper threshold may be dynamically calculated, for example, by rejecting data points that are more than a specified number of standard deviations from the median. Other data cleansing procedures are also possible.

[0045] In embodiments, the mean value of the (remaining) SIA values ​​for each given mass is determined, and the model is fitted to these mean SIA values. The mean SIA values ​​may be corrected to account for data excluded by the cleaning step, for example, by fitting an expected distribution to the cleaned data and taking the mean of the fitted distribution rather than the mean of the cleaned data.

[0046] In a further embodiment, the SIA data may be corrected for detection efficiency before a model is fitted to the corrected data. As described above, the detector may include one or more electron multiplication stages followed by a conversion dynode. The conversion dynode may be configured to produce one or more secondary electrons when an ion strikes the conversion dynode. However, it is recognized that there is a non-zero probability that zero secondary electrons will be produced by the conversion dynode when an ion strikes it. Such zero secondary electron events do not produce any measurable signal, and therefore, the average SIA value of experimentally measured SIA data is likely to be an overestimate of the actual average SIA. Furthermore, the likelihood of a zero secondary electron event is higher for masses that provide lower SIA because ions at these masses are converted to electrons less efficiently (i.e., produce fewer secondary electrons per ion).

[0047] Thus, in an embodiment, the SIA data (or average SIA values) are corrected to account for the possibility of zero secondary electron events, and the model is fitted to the average of the corrected SIA values ​​(or the corrected average SIA values). It has been found that by taking into account zero secondary electron events, the accuracy of the correction function can be significantly improved.

[0048] Corrections to SIA data may be made in any suitable manner and may take any suitable form. In embodiments, the magnitude of the correction applied to a particular SIA value at a particular mass (or to a particular average SIA value) depends on the average SIA value at the particular mass, thus taking into account, for example, the variation in the probability of a zero secondary electron event as a function of the average SIA. As explained further below, the correction may be obtained by assuming a quasi-Poisson probability distribution for the ion-to-electron conversion process and using an estimate of the secondary electron yield (SEY) at a particular mass (i.e., the average number of secondary electrons produced by a dynode in response to an ion striking the dynode). An estimate of the SEY at a particular mass may then be obtained by comparing experimentally obtained SIA data with Monte Carlo simulations.

[0049] The model fitted to the (optionally processed and / or corrected) experimental data can take any suitable form. In particular, the model may be a semi-empirically derived model. In this regard, as explained in more detail below, theoretical considerations generally suggest that the mass dependence of SIA has a maximum value at a transition mass that depends on the characteristics of the dynode, that SIA increases with increasing mass below the transition mass, and that SIA decreases with increasing mass above the transition mass. According to embodiments, a model has been created that takes these (and other) theoretical considerations into account, but is otherwise experimentally derived. Thus, in embodiments, the mass dependence of the model (and correction function) has a maximum value at the transition mass, increases with increasing mass below the transition mass, and decreases with increasing mass above the transition mass.

[0050] The transition mass can have a suitable value depending on the characteristics of the dynode, for example, the material from which the dynode is formed. A suitable transition mass is in the range of 100-200 Da. For example, if the dynode is a stainless steel dynode, the transition mass may have a value of about 140 Da.

[0051] The inventors have found that a particularly suitable model (and therefore correction function) has the form:

number

number

[0052] The inventors have found this model to be particularly well-suited for the present embodiment when it is desired to establish a global correction function that can be used across multiple different instruments of the same class and in the context of analyzers used to analyze a wide range of different analytes. This global correction function can be defined for multiple analyzers that have at least the same instrument geometry, use the same dynode materials, and are configured so that ions arrive at the detector at the same angle of incidence, etc. In this regard, the inventors have recognized that the global correction functions of various embodiments can be calibrated to each individual analyzer in a particularly straightforward manner, for example, by simple scaling. Furthermore, the use of this relatively simple power law model prevents over-fitting to experimental data, which allows the subsequently obtained curve to be used for multiple different instruments of the same class and for a wide range of different analytes as appropriate.

[0053] As described above, the global correction function can be used as a correction function for one or more different analyzers of the same class (i.e., analyzers having at least the same instrument geometry, using the same dynode material, and configured so that ions arrive at the detector at the same angle of incidence, etc.). In certain embodiments, the global correction function is used to derive correction functions for different analyzers of the same class. Other classes of analyzers, or similar analyzers using other materials, for example, for the dynodes, may require recalibration of all free parameters of the correction function. In certain embodiments, the correction function for a particular analyzer is derived by calibrating the global correction function to the particular analyzer, e.g., by vertically scaling it (i.e., by scaling its SIA axis). This calibration can include determining the average SIA value of the analyzer at a particular mass or m / z, e.g., when the particular mass or m / z is at or near a transition mass. The correction function for the analyzer can then be derived by scaling (the SIA axis of) the global correction function according to the average SIA value. This can include scaling the SIA axis of the global correction function so that the value of the scaled correction function at a particular mass or m / z is equivalent to the determined average SIA value for the analyzer at that particular mass or m / z.

[0054] In embodiments, each individual analyzer can also be calibrated for its operable range of detector voltages (i.e., the range of voltages that can be applied across the detector's electron multiplication stages, thereby providing different levels of amplification). This calibration can be used to derive an SIA vs. detector voltage calibration curve for each individual analyzer, and can be repeated, for example, periodically to account for detector aging. This calibration function can be combined with a correction function derived from the global correction function described above, i.e., whereby the ion-to-electron conversion efficiency can be assumed to be the same for all analyzers of the same class. Thus, the above calibration describes the mass or m / z dependence of the ion-to-electron calibration, while a second calibration describes the detector voltage dependence of electron multiplication (or, in some cases, electron-to-photon conversion, followed by PMT). Because ion-to-electron conversion and electron multiplication are entirely independent processes, these two calibrations can be combined to model the entire detector signal produced by a single ion.

[0055] For example, if an analyzer can be used to analyze only ions of the same charge, the correction function could describe the SIA as a function of mass or m / z only. However, in further embodiments, the correction function also describes the relationship between the average single-ion area and charge for the analyzer. For example, the correction function can be a surface that describes the average SIA as a function of mass or m / z (or equivalently, as a function of ion velocity or arrival time) and as a function of charge. In certain embodiments, the correction function provides an (approximate) average SIA for the analyzer as a function of ion mass or m / z (or equivalently, as a function of ion velocity or arrival time) and as a function of charge. Thus, the number of ions that contributed to a particular ion peak can be estimated by obtaining the average SIA from the correction function at the mass or m / z and charge associated with the ion peak and dividing the measured area of ​​the ion peak by this value.

[0056] In embodiments, the correction function describes the relationship between the average SIA and ion charge for the analyzer over the charge range of interest for the analyzer, such as from about 1 elementary charge to about 10 elementary charges, 15 elementary charges, 20 elementary charges or more, etc. In certain embodiments, the correction function describes the relationship between the average SIA and ion charge over the range of from about 1 elementary charge to about 15 elementary charges.

[0057] In these embodiments, the relationship between average single ion area and charge can be re-derived by fitting the model (parameters of the model) to the analyzer's experimentally acquired SIA data. The experimentally acquired SIA data can include experimental measurements of the analyzer's SIA for ions across most or all of the charge range of interest. In certain embodiments, the model (and correction function) increases monotonically linearly with charge.

[0058] Various embodiments also extend to generating correction functions.

[0059] Accordingly, a further aspect provides a method for determining a correction function for a mass analyzer, the method comprising: analyzing a plurality of single ions using a mass analyzer, the plurality of single ions including ions having mass m / z and / or charge spanning most or all of the mass m / z and / or charge range of interest for the analyzer; generating single ion area (SIA) data by determining an area of ​​each ion peak of the plurality of ion peaks produced by the analyzer in response to detecting the plurality of single ions; and deriving a correction function by fitting the model to the SIA data.

[0060] This aspect may include, and in embodiments includes, any one or more or each of the optional features described herein.

[0061] Thus, for example, the SIA data may be processed, corrected, and / or averaged SIA data, e.g., as described above. The mass or m / z dependence of the model (and correction function) may have the form described above, e.g.,

number

number

[0062] The SIA data can be representative analyzer SIA data, the correction function can be a global correction function, and the correction function for another analyzer in the same class can be derived from the global correction function by scaling the global correction function, for example, as described above. The correction function can be used to estimate the number of ions that contributed to an ion peak produced by the analyzer, for example, as described above.

[0063] A further aspect provides a method of operating an analytical instrument, such as a mass spectrometer, comprising an ion source and a mass analyzer, the method comprising: generating ions in an ion source; analyzing the ions using a mass analyzer to generate a signal; and analyzing the signal using the methods described above.

[0064] This aspect may include, and in embodiments includes, any one or more or each of the optional features described herein.

[0065] Thus, for example, the analyzer may be a ToF analyzer including an ion trap disposed at the beginning of an ion path and an ion detector disposed at the end of the ion path. The ion trap may be configured to receive ions and accumulate packets of ions. The ion trap may be configured to inject each accumulated packet ion into the ion path, whereupon the packet ions can travel along the ion path to the detector.

[0066] The method can include estimating the total number of ions in the ion packet that generated the signal. The estimated total number of ions can then be used to control the total number of ions in the ion packet that are subsequently stored in the ion trap. The number of ions in the ion packet that are subsequently stored in the ion trap can be controlled by controlling the storage time of ions in the ion trap (e.g., fill time) (i.e., with an automatic gain control (AGC) method). Thus, the method can include using the estimated total number of ions in the first ion packet to control the storage time of ions in the ion trap (e.g., fill time) for a subsequent ion packet. Generally, the fill time can be determined based on the estimated total number of ions in one or more previous ion packets, the known fill time of one or more previous ion packets into the ion trap, and a desired or target maximum number of ions in the ion trap.

[0067] A further aspect provides a non-transitory computer readable storage medium storing computer software code that, when executed on a processor, performs the above-described method.

[0068] A further aspect provides a control system for an analytical instrument, such as a mass spectrometer, the control system being configured to cause the analytical instrument to perform the above-described method.

[0069] A further aspect provides an analytical instrument, such as a mass spectrometer, comprising an ion analyzer and a control system as described above.

[0070] A further aspect provides an analytical instrument, the analytical instrument comprising: a mass analyzer; 1. A control system comprising: receiving a signal produced by a mass analyzer, the signal including one or more ion peaks; determining an area of ​​a first ion peak of the one or more ion peaks; and a control system configured to estimate the number of ions contributed to the first ion peak by (i) determining a correction to be applied to the area of ​​the first ion peak from a correction function, the correction function describing the relationship between an average single ion area and ion mass, mass-to-charge ratio, and / or charge for the mass analyzer, and (ii) applying the correction to the area of ​​the first ion peak.

[0071] These aspects and embodiments may, and in embodiments do, include any one or more or each of the optional features described herein.

[0072] Thus, for example, the analytical instrument may be a mass spectrometer, and the mass analyzer may be a time-of-flight (ToF) mass analyzer, such as a multi-reflecting time-of-flight (MR-ToF) analyzer. The time-of-flight (ToF) mass analyzer may include an ion trap configured to inject ions into a drift region of the analyzer. The control system may be configured to cause the instrument to perform any one or more, or each, of the methods described above.

[0073] Various embodiments will now be described in more detail with reference to the accompanying drawings. [Brief explanation of the drawings]

[0074] [Figure 1] 1 illustrates a schematic representation of an analytical instrument according to an embodiment; [Figure 2] 1 illustrates a schematic representation of a multi-reflecting time-of-flight mass analyzer according to an embodiment. [Figure 3] 1 illustrates a schematic diagram of a process for detecting ions in an ion analyzer, according to an embodiment. [Figure 4A] Figure 4B shows a histogram containing the core density estimates of SIAs recorded for selected species, and Figure 4B shows the mean values ​​of their corresponding SIA distributions. [Figure 4B] Figure 4B shows a histogram containing the core density estimates of SIAs recorded for selected species, and Figure 4B shows the mean values ​​of their corresponding SIA distributions. [Figure 5] Box plots of raw SIA data are shown along with the mean values ​​of the washed data. [Figure 6] The best fit parameters and regression lines for SIA as a function of incident velocity are shown. [Figure 7] The best-fit average SIA as a function of mass is shown. [Figure 8A] The pulse area distribution data is shown along with Monte Carlo simulated outputs with the same mean SIA and different ion-to-electron conversion values, and FIG. 8B shows the corresponding best fit to the pulse area distribution data. [Figure 8B] The pulse area distribution data is shown along with Monte Carlo simulated outputs with the same mean SIA and different ion-to-electron conversion values, and FIG. 8B shows the corresponding best fit to the pulse area distribution data. [Figure 9] Detection efficiency as a function of SEY

number

[0075] Detailed Description 1 is a schematic illustration of an analytical instrument, such as a mass spectrometer, that may be operated in accordance with embodiments. As shown in FIG. 1, the analytical instrument includes an ion source 10, one or more ion transfer stages 20, and an analyzer 30.

[0076] The ion source 10 is configured to generate ions from a sample. The ion source 10 can be any suitable continuous or pulsed ion source, such as an electrospray ionization (ESI) ion source, a MALDI ion source, an atmospheric pressure ionization (API) ion source, a plasma ion source, an electron ionization ion source, or a chemical ionization ion source. In some embodiments, two or more ion sources can be provided and used. The ions can be of any suitable type to be analyzed, such as small and large organic molecules, biomolecules, DNA, RNA, proteins, peptides, fragments thereof, etc.

[0077] The ion source 10 may optionally be coupled to a separation device, such as a liquid chromatography separation device or a capillary electrophoresis separation device (not shown), whereby the sample that is ionized in the ion source 10 comes from the separation device.

[0078] The ion transport stage 20 is positioned downstream of the ion source 10 and may include an atmospheric pressure interface and one or more ion guides, lenses and / or other ion optical devices configured to enable the transport of some or most of the ions generated by the ion source 10 from the ion source 10 to the analyzer 30. The ion transport stage 20 may include any suitable number and configuration of ion optical devices, for example, optionally including any one or more of one or more RF and / or multipole ion guides, one or more ion guides for cooling ions, one or more mass selective ion guides, etc.

[0079] The analyzer 30 is positioned downstream of the ion transmission stage 20 and is configured to receive ions from the ion transmission stage 20. The analyzer is configured to analyze the ions to determine physicochemical properties of the ions, such as their mass or mass-to-charge ratio. To do this, the analyzer 30 is configured to pass the ions through a detector. The instrument may be configured to determine the physicochemical properties of the ions from a signal measured by the detector. The instrument may be configured to produce a spectrum of the analyzed ions, such as a mass spectrum.

[0080] The analyzer 30 may be any suitable mass analyzer, such as a time-of-flight (ToF) mass analyzer, an ion trap mass analyzer, or a quadrupole mass analyzer.

[0081] In certain embodiments, analyzer 30 is a time-of-flight (ToF) mass analyzer, e.g., configured to determine the mass-to-charge ratio (m / z) of ions by passing the ions along an ion path in a drift region of the analyzer, which is maintained in a high vacuum (e.g., <1×10 -5 The pressure in the drift region is maintained at a pressure of 1000 psi (mbar). Ions may be accelerated by the electric field into the drift region and detected by an ion detector located at the end of the ion path. The acceleration may cause ions with relatively low mass-to-charge ratios to achieve relatively high velocities and arrive at the ion detector before ions with relatively high mass-to-charge ratios. As such, ions arrive at the ion detector after a time determined by their velocity and the length of the ion path, allowing the ion's mass-to-charge ratio to be determined. Each ion or group of ions arriving at the detector may be sampled by the detector, and the signal from the detector may be digitized. A processor may then determine a value indicative of the time-of-flight and / or mass-to-charge ratio ("m / z") of the ion or group of ions. Data for multiple ions may be collected and combined to generate a time-of-flight ("ToF") spectrum and / or a mass spectrum.

[0082] 1 is merely schematic, and the analytical instrument can, and in embodiments does, include any number of one or more additional components. For example, in some embodiments, the analytical instrument includes a collision or reaction cell for fragmenting or reacting ions, and the ions analyzed by the analyzer 30 can be fragment or product ions created by fragmenting or reacting parent ions generated by the ion source 10.

[0083] 1, the instrument is under the control of a control unit 40, such as a suitably programmed computer, which controls the operation of the various components of the instrument, including the analyzer 30. The control unit 40 is also capable of receiving and processing data from the various components, including the detector according to embodiments described herein.

[0084] 2 illustrates schematically the details of one exemplary embodiment of analyzer 30. In this embodiment, analyzer 30 is a multi-reflecting time-of-flight (MR-ToF) mass analyzer.

[0085] 2, the multi-reflecting time-of-flight analyzer 30 includes a pair of ion mirrors 31, 32 spaced apart and facing each other in a first direction X. The ion mirrors 31, 32 are elongated along an orthogonal drift direction Y between first and second ends.

[0086] An ion source (injector) 33, which may be in the form of an ion trap, is located at one end (first end) of the analyzer. The ion source 33 may be positioned and configured to receive ions from the ion transport stage 20. Ions may be accumulated within the ion source 33 before being injected into the space between the ion mirrors 31, 32. As shown in Figure 2, ions may be injected from the ion source 33 at a relatively small injection angle or drift direction velocity to produce zigzag ion trajectories, whereby different oscillations between the mirrors 31, 32 are separated in space.

[0087] One or more lenses and / or deflectors may be positioned along the ion path between the ion source 33 and the ion mirror 32 that the ions first encounter. For example, as shown in Figure 2, a first out-of-plane lens 34, an injection deflector 35, and a second out-of-plane lens 36 may be positioned along the ion path between the ion source 33 and the ion mirror 32 that the ions first encounter. Other configurations may be possible. In general, the one or more lenses and / or deflectors may be configured to suitably condition, focus, and / or deflect the ion beam, i.e., to cause the ion beam to follow a desired trajectory through the analyzer.

[0088] The analyzer also includes another deflector 37 positioned along the ion path between ion mirrors 31, 32. As shown in Figure 2, deflector 37 may be positioned approximately equidistant along the ion path between ion mirrors 31, 32 after its first ion mirror reflection (at ion mirror 32) and before its second ion mirror reflection (at the other ion mirror 31).

[0089] The analyzer also includes a detector 38. The detector 38 may be any suitable ion detector configured to detect ions and, for example, record the intensity and arrival time associated with the arrival of the ions at the detector. Suitable detectors include, for example, one or more conversion dynodes, optionally followed by one or more electron multipliers, or the like.

[0090] To analyze ions, ions may be injected from the ion source 33 into the space between the ion mirrors 31, 32 such that the ions follow a zigzag ion path that makes multiple reflections in the X direction between the ion mirrors 31, 32 while (a) drifting along drift direction Y toward the opposite (second) ends of the ion mirrors 31, 32, (b) reversing their drift direction velocity near the second ends of the ion mirrors 31, 32, and then (c) drifting back along drift direction Y toward the deflector 37. The ions may then be transferred from the deflector 37 to the detector 38 for detection.

[0091] In the analyzer of Figure 2, both ion mirrors 31, 32 are tilted with respect to the X direction and / or the drift Y direction. Alternatively, it would be possible to tilt only one of the ion mirrors 31, 32, for example, with the other ion mirror 31, 32 positioned parallel to the drift Y direction. In general, the ion mirrors are at a non-constant distance from each other in the X direction along their length in the drift direction Y. The drift direction velocity of ions towards the second end of the ion mirror is opposed by an electric field resulting from the non-constant distance of the two mirrors from each other, which causes the ions to reverse their drift direction velocity near the second end of the ion mirror and drift back along the drift direction towards the deflector 37.

[0092] The analyzer shown in Figure 2 further comprises a pair of correction stripe electrodes 39. Ions traveling down the drift length are slightly deflected as they pass each mirror 31, 32, and the additional stripe electrodes 39 are used to correct for time-of-flight errors caused by changes in the distance between the mirrors. For example, the stripe electrodes 39 are electrically biased so that the period of ion oscillation between the mirrors is substantially constant along the entire drift length (despite the fact that the distance between the two mirrors is not constant). The ions are eventually reflected back down the drift space and focused onto a detector 38.

[0093] Further details of the tilted mirror multi-reflecting time-of-flight mass analyzer of FIG. 2 are described in US Pat. No. 9,136,101.

[0094] It should be noted that in general the analyser 30 may be any suitable type of mass analyser or time-of-flight (ToF) mass analyser, for example the analyser may be a single lens type multi-reflecting time-of-flight mass analyser, such as that described in GB Patent No. 2,580,089.

[0095] Generally, time-of-flight (ToF) mass analyzers with ion impact detectors exploit the property that the transit time of ions in an electrostatic field is proportional to the square root of the ion's mass-to-charge ratio (m / z). Ions are simultaneously ejected from an ion source (e.g., an orthogonal accelerator or radio frequency ion trap), accelerated to a desired energy, and travel a specific distance before striking an ion detector. Similarly, in ion trap and quadrupole mass analyzers, ions ultimately strike an ion detector.

[0096] Figure 3 shows the ion 3 illustrates schematically the process by which a packet of ions 50 is detected by detector 38. In particular, Figure 3 illustrates the ion-to-electron conversion process and the subsequent signal processing chain.

[0097] As shown in Figure 3, the detector 38 comprises a conversion dynode 51 followed by one or more electron multiplier stages 53. ion When a packet of ions 50 strikes the conversion dynode 51, n e Secondary electrons 52 are produced. ion A packet of ions is expressed as n with a conversion factor γ e are converted into secondary electrons.

[0098] The secondary electrons 52 are then amplified by one or more electron multiplier stages 53 to produce a signal indicative of the intensity of the ions 50 received at the conversion dynode 51 as a function of time. The one or more stages of electron multiplication 53 are amplified by a gain factor g em provides a signal increase with

[0099] The generated signal is recorded by data acquisition electronics 54, such as a digitizer, which may be either a time-to-digital converter (TDC) or an analog-to-digital converter (ADC). The analog-to-digital conversion stage 54 may be further enhanced by another gain factor g sw Introduce.

[0100] In the embodiment shown in Figure 3, detector 38 comprises a conversion dynode 51 followed by one or more electron multiplier stages 53. However, it should be noted that in general, detector 38 may include one or more intermediate conversion processes in addition to those shown in Figure 3, such as, for example, one or more scintillator crystals, photomultiplier tubes (PMTs), etc.

[0101] As shown in Figure 3, the detection process results in time-resolved peaks 55, where each peak is produced by ions of substantially the same mass-to-charge ratio. In a ToF analyzer, the travel distance is substantially the same for all of the ions 50, so the ion arrival time can be used to determine the mass-to-charge ratio (m / z) of the ion, which can then be used for ion identification.

[0102] Furthermore, the area S of a time-resolved peak can be used to determine the number of ions that contributed to that peak, which can then be used for quantification. The final signal S is obtained by peak-by-peak integration of the digitized signal counts over the arrival time axis, and is herein referred to as the "ion area," "ion peak area" (shaded area in FIG. 3), or, more specifically, the limit n ion = 1 is called the "single ion area" (SIA).

[0103] The integrated digitized signal S and the initial number of incident ions n ion The effective gain factor G between ion It has been recognized that ρ is not constant, but has a clear dependence on the statistical properties of the incident ions, the different conversion and amplification stages, and on the mass and charge state of the ions.

[0104] Figure 4A shows a histogram containing core density estimates of SIAs recorded for a selected species (singly charged MRFA peptides) using an analyzer of the type shown in Figure 2. To clean the raw input data, a lower SIA threshold (near the estimated noise level) and a consistently determined upper SIA threshold were applied. Figure 4B shows the corresponding distribution of SIA means, which demonstrate a clear mass dependence for singly charged ions. More sophisticated data analysis and cleaning procedures (not shown) allowed the identification of m tr A unimodal shape of the SIA is suggested across the mass curve with a maximum near a transition mass of ≈150 Da.

[0105] In the following references, the secondary electron yield (SEY) during ion impact on a (metallic or semiconducting) conversion dynode as a function of mass is derived from first principles or determined semi-heuristically for a narrow range of masses and incident ion species: (i) Beuhler, RJ and L. Friedman, "A Model of Secondary Electron Yields from Atomic and Polyatomic Ion Impacts on Copper and Tungsten Surfaces Based upon Stopping-power Calculations," Journal of Applied Physics 48, no. 9 (September 1, 1977): 3928-3936. (ii) Liu, Ranran, Qiyao Li, and Lloyd M. Smith, “Detection of Large Ions in Time-of-Flight Mass Spectrometry: Effects of Ion Mass and Acceleration Voltage on Microchannel Plate Detector Response,” Journal of the American Society for Mass Spectrometry 25, no. 8 (August 2014): 1374-83. (iii) Moshammer, R. and R. Matthaus, "Secondary Electron Emission by Low Energy Ion Impact" Journal de Physique Colloques 50, no. C2 (1989): C2-111-C2-120. (iv) Sternglass, EJ "Theory of Secondary Electron Emission by High-Speed ​​Ions" Physical Review 108, no. 1 (October 1, 1957): 1-12.(v) Westmacott, G., W. Ens, and KG Standing. "Secondary Ion and Electron Yield Measurements for Surfaces Bombarded with Large Molecular Ions" Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms 108, no. 3 (February 1, 1996): 282-89.

[0106] These estimates are typically based on the pioneering work of Bloch, Bethe, and Bohr on the stopping of material power on charged particles, which is the average energy loss per transmission distance,

number

number

[0107] These various studies each refer to very different mass ranges, starting from elementary particles up to atomic number Z≈100, and extending to the detection of large molecular ions (m≈1000-100,000 Da). Also, the linear velocity dependence of SEY has long been assumed and used in the art, with masses above the transition mass (m > m tr ) (See, for example, Ingram, Hayden, and Hess, "Mass Spectrometry in Physics Research," National Bureau of Standards (US), Circ. 522, 1953, p. 257).

[0108] Despite literature and experimental evidence showing that the SEY in ToF analyzers has a clear dependence on ion mass and charge, known mass spectrometers do not systematically correct for the mass and / or charge dependence of the ion-electron conversion process, which effectively amounts to approximating the data points in Figure 4B by a constant horizontal line.

[0109] Furthermore, there have been no successful attempts to find a global SEY(m,z) curve that best fits experimental data over the entire range of mass-to-charge ratios (e.g., m / z ≈ 50 to over 15,000) relevant to a ToF-MS analytical instrument, such as an MR-ToF instrument of the type illustrated in Figure 2, using an accelerating voltage of 10 kV.

[0110] Additionally, in known analytical instruments, it can be seen from the statistical nature of the conversion process and the subsequent detection tool chain that the measured SIA does not include true negative and / or zero secondary electron events that would occur with a probability that would not actually disappear. Corresponding corrections based on scaling of detection efficiency < 1 are not made in known analytical instruments.

[0111] Therefore, embodiments provide correction functions that can be used to correct for these mass and charge dependencies, as well as detection efficiency. In particular, embodiments provide correction functions that describe the relationship between SIA, ion mass, and charge across the entire operating parameter space of a mass analyzer. This allows for systematic correction of the dependence of ion area on mass and / or charge in a particularly accurate and straightforward manner.

[0112] To provide an accurate estimation of the number of incident ions from a measured peak area (e.g., for quantification), several distinct but related problems are addressed by embodiments. First, a model of SEY, or equivalently SIA, as a function of mass (or equivalently velocity) and charge is provided. Second, the measured SIA is corrected for zero-electron events by estimating the detection efficiency. This can be obtained by leveraging known statistical properties of the conversion process, but requires as input an initial estimate of the SEY itself. Thus, an estimate of the conversion of initial ions to (secondary) electrons is provided. By addressing these individual problems, embodiments facilitate an accurate estimation of the total number of incident ions from a measured ion peak area.

[0113] In this regard, the inventors have recognized that the large mass and charge parameter space and inherent complexity (including, for example, their different conformations) of molecular ions typically analyzed using ToF analyzers and other analyzers in the life sciences make it highly unlikely that first-principles analysis will be feasible for the problem at hand. Accordingly, embodiments provide correction functions derived by fitting semi-empirical globally valid models to masses from experimentally acquired SIA data for ToF-MS (or other MS) instruments.

[0114] In general, embodiments include: (i) fitting a semi-empirical, globally valid model to the mass / velocity and charge dependence of SIA and SEY based on experimental mass analyzer data; (ii) estimating SEY by fitting measured, fully corrected, and algorithmically cleaned SIA data to model-based predictions from Monte Carlo (MC) simulations; (iii) using the estimated SEY to correct the initially acquired SIA data with a detection efficiency factor obtained based on statistical properties of the electron multiplication process; and (iv) thereby facilitating accurate prediction of the number of incident ions, e.g., for label-free quantification of ToF-MS data.

[0115] Thus, embodiments provide improved quantification of incident ion abundance for mass spectrometry, particularly time-of-flight mass spectrometry (ToF-MS), based on improved (heuristic and statistical) models for the mass and charge dependence of the ion-to-electron conversion process and the subsequent secondary electron yield of the detection tool chain.

[0116] The following description provides a representative example of a workflow based on raw SIA data to illustrate an embodiment and the algorithm steps used therein. For clarity, the mass-dependent correction for singly charged ions and the charge-dependent additional correction are presented separately. As used herein, all references to mass m (which, when used in the context of mass spectra without describing the charge state, means charge number z=1) should be understood to be provided in units of Daltons (Da) unless otherwise stated, and the dimensionless m / z is the quotient of mass and charge number in units of Daltons.

[0117] 1. Mass-dependent correction for singly charged species First, for a given acceleration and detector voltage, SIA data is recorded over a sufficiently large mass range for singly charged species.

[0118] Figure 5 shows an excerpt of a boxplot of raw SIA data for an MR-ToF instrument of the type shown in Figure 2, along with the average values ​​of the cleaned data, and refers to the exact same data set as shown in Figure 4. The data shown in Figure 5 covers the mass range m / z ≈ 50-8,000, but SIA data over any desired mass range can be acquired in the same way. Charge-dependent data can also be acquired, for example, for a selected set of ion masses (as described in Section 2 below).

[0119] In an MR-ToF instrument of the type shown in FIG. 2 , a slight shift, e.g., <5%, in the voltage applied to the ion mirrors 31, 32 is sufficient to completely defocus the ion packet to a single ion at the detector 38. For example, an ion with m / z of 500, which would otherwise be focused to a peak of <4 ns, can spread over a period of microseconds. This allows the signal to be split into approximately hundreds of individual peaks, each with a width of approximately 1 ns. The raw data shown in FIG. 5 was acquired using this technique, which allows single-ion data to be acquired in a relatively fast and straightforward manner, e.g., without the need to attenuate the ion beam. This technique is also compatible with other types of ToF analyzers, e.g., simpler single-reflection ToF analyzers, although the level of defocusing must be greater or the number of ions must be smaller.

[0120] 1.1 Obtaining the best-fit curve for the initial mass-dependent SIA Due to the large variability in the input data, and to account for additional effects such as noise thresholds and peak ringing, the set of raw SIA data points for a given mass is first clipped using a fixed lower threshold and a dynamically calculated upper threshold obtained from iterative sigma clipping around the data median (successively rejecting values ​​beyond a specified number of standard deviations from the median).

[0121] The mean of the remaining SIAs is then used for a given mass, and (for simplicity) we assume that the washed distribution is of a shape that can be adequately approximated by a Poisson distribution. This approximation is motivated by the nature of the underlying governing physical processes: the Poisson distribution is completely determined by a single parameter corresponding to the mean.

[0122] These average SIAs are then calculated using the kinetic energy equation

number

[0123] To provide a simple semi-empirical average v dependence with the desired properties (i.e., a power-law dependence with a maximum in the dynode material-dependent transition rate),

number

[0124] The inventors have found this model to be particularly suitable in the present situation, where it is desirable to establish a relatively general model that can be used across multiple different instruments of the same class and in the context of analyzing a wide range of different analytes. The use of this relatively simple power law model prevents overfitting to the experimental data, thereby allowing the resulting curves to be used for multiple different instruments of the same design and for a wide range of different analytes.

[0125] While this form is particularly advantageous, it should be noted that the particular form of the model may be selected as desired, and embodiments are not limited to the particular form described herein. Other correction functions may provide less critical but acceptable performance, or may provide good performance over a limited mass and / or charge range, or may provide good, but not necessarily generalizable, performance for particular types of molecular species.

[0126] A nonlinear regression of the model is performed on the experimentally obtained data (using an appropriate regression library function). The regression line and best fit parameters obtained for the SIA as a function of incident velocity are shown in Figure 6. As shown in Figure 6, for this example, a = 7.22e-4, b = 1.11, and v = 124469 (in units corresponding to the depiction in Figure 6). The transition velocity is v tr =138706 m / s.

[0127] The result is then transformed back to the mass domain to obtain the transition mass m Tr and overall behavior. In the mass domain,

number

[0128] Without additional correction for non-uniform detection efficiency (described below), the number of ions can now be obtained approximately using this curve, e.g., as the measured peak area divided by the best-fit SIA(m) for any measured mass of interest.

[0129] 1.2 Correction for detection efficiency Due to the statistical nature of the ion-to-electron conversion process and the subsequent signal processing chain, there is a finite probability of a zero-electron event.

[0130] In particular, once one or more secondary electrons are produced by the dynode 51, they are multiplied in a positive (Poissonian) manner by one or more stages of electron multiplication 53. However, there is a chance that the dynode 51 will respond to an ion but not produce any secondary electrons. This can be seen in Figure 8A, which shows a simulated SIA distribution. As shown in Figure 8A, most of the distribution is approximately Gaussian, but there is a triangular peak at zero due to the zero-electron event.

[0131] This trend is not accounted for in the curves of Figure 7. Zero secondary electron events do not produce any measurable signal, and therefore the average SIA values ​​in Figure 7 are likely to be overestimates of the actual average SIA. Furthermore, the likelihood of zero secondary electron events is higher for masses that provide lower SIA because ions at these masses are converted to electrons less efficiently (i.e., produce fewer secondary electrons per ion). Therefore, further embodiments employ the following correction scheme:

[0132] Fixed mass m≒m Tr (Ideally, high

number

number

[0133] Figure 8 shows MC modeling of the initial ion-to-electron conversion, and hence the SEY, for a given SIA. Figure 8A shows pulse area distribution data along with simulated outputs with the same average SIA and different ion-to-electron conversion values ​​(1 to 7 electrons). The various curves represent MC simulations of the initial electron conversion process, assuming different numbers of secondary electrons are produced in response to a single ion at the dynode 51.

[0134] Figure 8B shows the corresponding best fit, which can be used to obtain an estimate for SEY. In the example shown in Figure 8, a particular stainless steel dynode is shown to produce an average of four electrons in response to each single ion.

[0135] A first estimate of the uncorrected mean abundance at other masses c0 is

number

[0136] If the probability distribution of the ion-electron conversion process is assumed to be (to a good approximation) quasi-Poisson, its efficiency can generally be obtained as follows:

number

[0137] this is,

number

number

[0138] FIG. 9 shows the resulting detection efficiency curve ε as a function of the (initially estimated) SEY.

[0139] Starting value mismatch

number

number

number

number

[0140] The corrected SIA forecast is finally

number

number

[0141] Since these statistical properties are defined on a microscopic scale, n is formally discretized to positive integers,

number

[0142] When applied to the macroscopically averaged SIA / SEY, any value with an initial guess c<1 will project to zero. This is the case in the non-physical boundary (the macroscopically averaged SEY cannot drop below 1 without a zero-electron contribution, so by definition, the microscopically defined

number

number

[0143] Figure 10 shows two real-valued branches of the efficiency solution near c0 = 1 (macroscopic case) and n = 1 (microscopic case).

[0144] Applying the full correction scheme to the initially obtained SIA values ​​resulted in a significantly modified velocity / mass dependence of the SIA data, highlighting the effectiveness of the full correction scheme as presented herein. This can also be understood from the initial mass dependence of the SIA, as higher masses have lower SIA, resulting in lower SEY and lower conversion efficiency.

[0145] Figure 11 shows the initial SIA data (black dots) and the efficiency-corrected SIA data (red dots), along with the new best-fit parameters and new regression line. Correction for zero-electron events results in the mass-dependent power-law exponent b / 2 (in this example,

number

[0146] 2. Charge dependence of SIA The following describes an additional correction for the dependence of SIA on charge state.

[0147] Figure 12A shows the charge dependence of SIA measured without mass-dependent correction over a mass range of approximately m≈600-1300 Da. This charge-dependent SIA data was obtained from two consecutive measurements on an instrument of the type shown in Figure 2 (the mass degrees of freedom are not shown for clarity). A highly accentuated dependence of SIA on charge state can be seen, with an apparent linear average trend.

[0148] Figure 12B shows the same data, but with the mass-dependent correction from the previous section applied. Additionally, the ion-electron conversion estimate based on the MC simulation is used to obtain an estimate that includes the number of secondary electrons (using the mass correction curve corresponds to normalizing the SEY to a common z=1 curve). Clearly, the linear trend is further accentuated.

[0149] Therefore, in an embodiment, the linear dependence of SIA on charge is included in the correction function.

[0150] 3. Using the correction function The correction function derived in the manner described above can be used as a global correction function for all individual analyzers of the same class, e.g., all analytical instruments of the type illustrated in FIG. 2 that have at least the same instrument geometry, use the same dynode material, and are configured so that ions arrive at the detector at the same angle of incidence. The global correction function is used to derive a correction function for each individual analyzer by scaling the global correction function to each individual analyzer. This can include, for example, measuring the average single ion area (SIA) for ions of a particular mass when the mass is at or near the transition mass, and then scaling the global correction function based on the measured average SIA. This can include scaling the SIA axis of the global correction function so that the value of the scaled correction function at the particular mass is equivalent to the determined average SIA value for the analyzer at the particular mass.

[0151] FIG. 13 illustrates a schematic diagram of a method for determining the number of ions that contributed to an ion peak, according to an embodiment.

[0152] As shown in FIG. 13, a particular (“ith”) ion peak is identified in the digitized signal produced by detector 38 and its peak area S i is initially determined (step 60), for example, by integrating the area of ​​the signal under the peak. The mass m i also optionally has its charge z i (step 62). This can be done as needed. For example, the mass-to-charge ratio (m / z) of an ion peak can be determined from its arrival time, its charge z can be determined from the experimental context and / or by taking into account related isotope patterns and / or neighboring charge state ion peaks, and its mass m can be determined from its m / z and charge z.

[0153] Next, m of the i-th ion peak i and / or charge zi Correction factor SIA for ion peaks using i (m i , z i ) is found from the (scaled) correction function SIA(m,z) (produced in the manner described above) (step 64). Finally, the number of ions that contributed to the ith ion peak

number

number

[0154] In this regard, it is noted that the peak area is known to increase linearly with the number of ions in the peak. Saturation effects can be corrected, for example, using any suitable known technique.

[0155] This process of determining the number of ions that contributed to an ion peak can be repeated for one or more ion peaks, or for each other ion peak, in the signal. The estimated numbers of ions that contributed to each ion peak in the signal can then be summed to estimate the total number of ions that contributed to the signal, e.g., to estimate the total number of ions in the ion packet that generated the signal.

[0156] This information can be used, for example, to quantify specific analytes within a sample.

[0157] Additionally or alternatively, this information can be used for so-called automatic gain control (AGC) methods. As mentioned above, in time-of-flight mass analyzers of the type illustrated in Figure 2, ion packets are accumulated in the ion trap 33 before being injected from the ion trap 33 into the drift region between the mirrors 31, 32. In these analyzers, it is necessary to accurately control the total number of ions in each ion packet accumulated in the ion trap 33, for example, to optimize the number of ions below, but as close as possible to, a limit of the ion trap 33, such as the space charge limit of the ion trap 33. This is typically done using so-called automatic gain control (AGC) methods, which control the number of ions in subsequent ion packets based on an accurate measurement or estimate of the number of ions in the ion packet.

[0158] Embodiments allow the total number of ions in an ion packet that produced a signal to be accurately estimated by considering only that signal. Each ion packet can contain ions spread over a wide range of masses and charges (with a wide range of corresponding SIAs), and therefore correcting for the dependence of ion area on mass and / or charge provides a significant improvement over the accuracy of the estimation of the number of ions in an ion packet and over AGC methods.

[0159] Figures 14A and 14B compare mass spectra of Pierce™ Flexmix™ calibration solutions with a mass range of 150 to 3000. The spectrum in Figure 14B has m / z-dependent single ion area correction applied, whereas the spectrum in Figure 14A does not. It can be observed that higher m / z species are relatively attenuated without correction, which has the effect of substantially hindering any attempt to quantify those species.

[0160] From the above, it will be appreciated that embodiments provide a comprehensive and self-contained scheme that can be used to obtain improved estimates of the number of initial ions contributing to an ion peak. Monte Carlo simulation and statistical correction are used to obtain a heuristically motivated correction scheme that spans a large mass (and charge state) range. This avoids the need to rely on first-principles models, whose generalizability to the broad range of life science fields is not feasible.

[0161] While the invention has been described with reference to various embodiments, it will be understood that various modifications can be made without departing from the scope of the invention as set forth in the appended claims.

Claims

1. 1. A method for analyzing a signal generated by a mass analyzer, the method comprising: receiving a signal produced by the mass analyzer, the signal including one or more ion peaks; determining an area of ​​a first ion peak of the one or more ion peaks; (i) correcting the number of ions that contributed to the first ion peak; determining a correction to be applied to the area of ​​the first ion peak from a function, the correction function describing a relationship between an average single ion area and an ion mass (m), mass to charge ratio (m / z), and / or charge (z) for the mass analyzer; and (ii) applying the correction to the area of ​​the first ion peak. A method wherein the correction function is determined by fitting a model to measured single ion area (SIA) data.

2. 2. The method of claim 1, wherein the correction function describes the relationship between average single ion area and ion mass or m / z for the analyzer over a m / z range, the m / z range being from about 25, 50, 75, or 100 to about 6000, 8000, 10,000, or 15,000.

3. The method of claim 2 , wherein the mass or m / z dependence of the correction function is continuous across the m / z range.

4. 4. The method of claim 1, wherein the mass or m / z dependence of the correction function has a maximum value at a transition mass, increases with increasing mass at masses below the transition mass, and decreases with increasing mass at masses above the transition mass.

5. The method of claim 4 , wherein the analyzer includes an ion detector comprising at least a dynode, and the transition mass depends on one or more properties of the dynode.

6. The mass or m / z dependence of the correction function is of the form [0010] where T is the kinetic energy of the ion, m is the mass of the ion, v is the velocity of the ion, and a, b, and v 0 The method according to any one of claims 1 to 3, wherein: ? is a best-fit parameter.

7. 4. The method of claim 1, wherein the correction function describes the relationship between average single ion area and charge for the analyzer over a charge range from 1 elementary charge to 10 elementary charges or more.

8. The method of claim 7 , wherein the charge dependence of the correction function varies linearly with charge.

9. The method according to any one of claims 1 to 3, wherein the correction function for the analyzer is obtained by scaling a global correction function.

10. determining an area of ​​one or more further ion peaks of the one or more ion peaks; estimating, for each of the one or more further ion peaks, a number of ions contributed to the one or more further ion peaks by: (i) determining a correction to be applied to the area of ​​the ion peak from the correction function; and (ii) applying the correction to the area of ​​the ion peak; The method of any one of claims 1 to 3, further comprising:

11. The method of claim 9 , wherein the global correction function is determined by fitting a model to measured single ion area (SIA) data.

12. 1. A method for determining a correction function for a mass analyzer, the method comprising: analyzing a plurality of single ions using a mass analyzer, the plurality of single ions including ions having mass m / z and / or charge spanning most or all of a mass m / z and / or charge range of interest for the analyzer; generating single ion area (SIA) data by determining an area of ​​each ion peak of a plurality of ion peaks produced by the analyzer in response to analyzing the plurality of single ions; and determining a correction function by fitting a model to the SIA data.

13. 13. The method of claim 1 or 12, wherein the single ion area (SIA) data is corrected SIA data derived from the raw SIA data by correcting the raw SIA data for detector efficiency.

14. 1. A method of operating an analytical instrument comprising an ion source and a mass analyzer, the method comprising: generating ions in the ion source; analyzing the ions using the mass analyzer to generate a signal; Analysing the signal using the method of any one of claims 1, 2 and 3.

15. 13. The method of any one of claims 1, 2, 3 and 12, wherein the mass analyzer is a time-of-flight (ToF) mass analyzer.

16. the time-of-flight mass analyzer comprises an ion trap, the method comprising: accumulating a first packet of ions in the ion trap; analyzing the first ion packet to generate a first signal; analyzing the first signal to estimate a total number of ions in the first ion packet; and storing a second packet of ions in the ion trap. The method of claim 15 , wherein the estimated total number of ions in the first ion packet is used to control the total number of ions in the second ion packet.

17. A non-transitory computer readable storage medium storing computer software code which, when executed on a processor, performs the method of any one of claims 1, 2, 3 and 12.

18. 13. A control system for an analytical instrument, said control system being configured to cause said analytical instrument to perform the method of any one of claims 1, 2, 3 and 12.

19. 20. An analytical instrument comprising an ion analyzer and a control system according to claim 18.

20. An analytical instrument comprising: A mass analyzer; 1. A control system comprising: receiving a signal produced by the mass analyzer, the signal including one or more ion peaks; determining an area of ​​a first ion peak of the one or more ion peaks; (i) correcting the number of ions that contributed to the first ion peak; a control system configured to estimate by (i) determining a correction to be applied to the area of ​​the first ion peak from a function, the correction function describing a relationship between an average single ion area and an ion mass, mass to charge ratio, and / or charge for the mass analyzer, and (ii) applying the correction to the area of ​​the first ion peak; The analytical instrument, wherein the correction function is determined by fitting a model to measured single ion area (SIA) data.

21. the mass analyzer is a time-of-flight mass analyzer with an ion trap, the instrument comprising: accumulating a first packet of ions in the ion trap; analyzing the first ion packet to generate a first signal; estimating a total number of ions in the first ion packet; and using the estimated total number of ions in the first ion packet to control the number of ions in a second ion packet stored in the ion trap.

Images