Panayot
Petrov
*,
Simon
Cowen
and
Heidi
Goenaga-Infante
LGC Ltd, National Measurement Laboratory, Queens Road, Teddington, TW11 0LY, UK. E-mail: Panayot.Petrov@lgcgroup.com
First published on 7th July 2021
Species-specific isotope dilution mass spectrometry (SS-IDMS) has been the calibration method of choice for high accuracy speciation analysis because it can correct for detector sensitivity drifts, matrix effects, and analyte loss during sample preparation and analysis. However, in many cases SS-IDMS calibration is either not applicable (e.g. for monoisotopic elements) or not feasible (e.g. limited by the cost and availability of like-for-like isotopically enriched species). The work presented here demonstrates the potential of a novel on-column species-specific internal calibration approach, which is based on the chromatographic injection of the same species of the analyte as the internal standard (IS), after the sample injection. It can compensate for on-column analyte losses and signal drift and can be applied with any detector capable of recording time-resolved data, provided that enough species resolution can be achieved. The feasibility of this novel calibration strategy for accurate quantitative elemental speciation in complex matrices is demonstrated here through the analysis of inorganic arsenic in rice. An expanded uncertainty (k = 2) of <10% was obtained for a mass fraction range of 60 to 300 μg kg−1 inorganic-As (i-As) in dry rice products. The method is currently used for the certification of i-As in baby food matrices to support Commission Regulation (EU) 2015/1006 in regard to the maximum levels of i-As in foodstuffs.
The two approaches used for internal standardisation in speciation analysis are continuous IS and discrete IS.7 A schematic diagram summarising the type of IS by point of introduction is presented in the ESI† (Diagram 1). The former (continuous IS) relies on constantly supplying a suitable and fixed amount and flow of element(s), which is merged with the analyte flow before reaching the detector. The disadvantage of this approach is the same as discussed above for total element determination, i.e. the difference in the behaviour of elements/species used as IS and the species of interest.8–10 The second approach (discrete IS) has several modifications, depending on how the IS has been introduced (e.g. within each sample as a surrogate standard or post-column). Details on the different types of internal standardisation approach are readily available in the literature.11–14 In summary, the post-column set-up includes a switching valve, positioned on the sample flow line, after the separation module (e.g. chromatographic column). An injection loop is filled with IS, which is then injected into the effluent. The two main types of discrete post column internal standardisation then are (1) IS injection immediately after the sample injection within a couple of seconds delay, needed for baseline estimation,7 and (2) IS injection after the elution of the strongest retained species. In both cases, the advantage of this approach compared to continuous IS is that the same element and even the same species can used as IS. It is possible to have more than one species as an internal standard by using an array of switching valves or even performing IS with MSA.15 For example, one of the IS species could be injected at the beginning of the separation and one at the end, although IS injections can be done so to elute at other retention times where no analytes are present. However, many researchers prefer to use only one species with ICP-MS detection, assuming equal species behaviour irrespective of their different volatility and therefore transport efficiency to the plasma. With more and more reports on species dependent ICP-MS sensitivity being published,16 the preferred option has shifted towards the use of an “average” internal standard consisting of species with sensitivity in-between the sensitivities of the other species of interest and if feasible, with similar chemical behaviour. Alternatively, the most abundant species or the species of particular scientific interest can be used as IS to minimise the bias from the inadequate IS correction. By employing these strategies, alone or in combination with MSA calibration, many of the known sources of bias can be efficiently mitigated. These include moderately sized low frequency noise (e.g. small sensitivity drift) and some matrix effects.4,9,17 Yet, all of the above being post-column corrections, they do not correct for any adverse effects accompanying the separation process. When samples with complex matrices are analysed, some of the target species can be either significantly more strongly retained on the separation column or eluted earlier compared to the respective standards. Additionally, peak splits have been reported. These have led to large ambiguity in species identification. Moreover, quality controls based on mass balance calculations will be blind to such adverse effects. As a consequence, the species-specific on-column losses are not compensated for by the common IS approaches and may negatively bias the results for the species of interest.16
In order to minimise the impact of the discussed potential sources of bias on the accuracy and precision of the speciation data, we propose a novel on-column species-specific internal standardisation approach as an alternative to SS-IDMS, the former being particularly attractive for mono-isotopic elemental species. With this approach, the drift effect is monitored and corrected for on a sample-to-sample basis. Both the analyte and internal standard are the same species, thus minimising the issues related with different species behaviours. The feasibility of this novel calibration strategy for accurate quantitative elemental speciation in complex matrices is investigated here through the analysis of i-As in rice samples in support of EU Commission Regulation 2015/1006.18 The method accuracy was assessed through analysis of certified reference materials and by recovery experiments on spiked samples. A full uncertainty budget indicating the main contributing factors has been estimated and reported. Moreover, the capability of the method to quantify more than one species under non-isocratic conditions was investigated through simultaneous determination of i-As and Methyl Arsonate (MA) in a standard mixture submitted to the same preparation procedure as the rice sample.
With this setup, the analyst can set the injection time and therefore precisely position the IS peak within the chromatogram for each particular speciation analysis. Since the valve is connected before the separation column, the operator should ensure that its operational pressure exceeds the separation column backpressure. This is a clear difference with the post-column IS approach, for which low-pressure valves can be used. Such a high pressure valve can be integrated inside the liquid chromatography (HPLC) compartment or be positioned just before the analytical column. Finally, the mass fraction of the analyte in the sample is calculated based on eqn (1) below:
![]() | (1) |
Eqn (1) is used for the calculation of the analyte mass fraction via on-column internal standardisation with a linear calibration curve, where CS is the mass fraction of the analyte, a and b are the slope and intercept of the built linear calibration curve and AS and AIS are the peak areas of the analyte and internal standard, respectively. The calibration curve is built once, in the beginning of the analysis, to calculate parameters a and b using linear regression (Y vs. X). For this purpose, the ratios of the peak areas of calibration standards to the injected (see Fig. 1) internal standards peak areas (ratio plotted on axis Y) for each calibration standard level are plotted against the respective calibration standard concentrations (plotted on axis X). By using MS Excel© array function Linest for Y vs. X, a and b are calculated. More details are presented in the Results and discussion section.
RF power | 1600 W |
Carrier Ar flow rate | 850 mL min−1 |
Makeup/dilution gas | 280 mL min−1 |
Sample/skimmer cones | Ni/Ni |
Spray chamber temp. | 2 °C |
Reaction gas (flow) | O2 (30% setting; ∼3 mL min−1) |
Data acquisition mode | TRA |
Points per spectral peak | 1 |
Acquisition mode | Tandem MS |
Monitored mass-shift | 75As–75As16O (m/z 75–>91) |
Integration time | 0.3 to 1.0 s (HPLC method dependent) |
The acquired data (as peak areas) were processed using laboratory developed MS Excel™ spreadsheets. Two different types of nebuliser have been used depending on the required flow; Micromist (Crawford Scientific, Strathaven, Scotland) performed better with lower signal RSDs for flows below 900 μL min−1 and Conikal (Crawford Scientific) was used for flows of 1 mL min−1 and above.
AEC column | PRP-X100, 250 mm × 4.1 mm id × 10 μm (PRP-X100, 50 mm × 4.1 mm id × 5 μm) |
Separation mode | Isocratic |
Eluent | 37 mM (55 mM) NH4HCO3 in 1% v/v CH3OH, pH(9) |
Eluent flow | 0.8 mL min−1 (1 mL min−1) |
Column temperature | 22 ± 2 °C |
Autosampler temperature | 22 ± 2 °C |
Sample throughput rate | 1.5 (7.5) per hour |
A bio-inert valve (600 bar, 6-ports, two positions) was installed on automatic valve driver 1290 Infinity model G1170A (both obtained from Agilent Technologies, Cheadle, UK) and used for the on-column IS injection. The valve was equipped with sampling loops made of PEEK tubing with nominal volumes between 20 and 1000 μL. The internal standard (stored in a 250 mL PTFE container, see Fig. 1) was delivered through a secondary medium pressure pump (model P 4.1S, Knauer, Northampton, UK) with a nominal flow rate of 20 μL min−1. The IS container lid was pierced at two positions to (1) allow fitting of the 1/16′′ IS pump inlet tubing, along with (2) another small (1/16′′) hole for pressure equilibration. The plumbing schematics are presented in Fig. 1.
Concentrated nitric acid UpA grade (mass fraction w(HNO3) ≥ 63%, density of ρ(HNO3) = 1.4 g mL−1) and hydrogen peroxide UpA grade (w(H2O2) ≥ 30%) were both purchased from Romil (Cambridge, UK). A mixture of 2% v/v HNO3 and 1% v/v H2O2 in water was then used for the extraction of i-As from 0.5 g rice sample (90 °C for 4 hours). The extracts were filtered through 0.22 μm syringe filters (33 mm, PES, Starlab Ltd, Milton Keynes, UK) and neutralised with 1:
1 v/v diluted ammonia (32%, HiPerSolv, VWR, Leighton Buzzard, UK) to a pH of 9.0 ± 0.5. Ammonium bicarbonate (Fluka, Sigma-Aldrich, Gillingham, UK) was dissolved in water (concentrations detailed in Table 2) and the pH was adjusted to 9.0 ± 0.2 with 1
:
1 (v/v) ammonia. The pH was measured using a FisherBrand Hydrus 300 pH meter (Fisher Scientific, Leicestershire, UK). Methanol (Romil low metal, Spa, Romil, Cambridge, UK) was finally added to the solution to a final concentration of 1% (v/v). The solution was used both as a mobile phase and diluent, and for internal standard preparation for the on-column IS calibration. The total dilution factor for both the samples and the matrix CRMs (see below) was 40 (1 + 39 w/w ±2%). Sodium Arsenate CRM 7912-a (NMIJ, Tsukuba, Japan) and sodium arsenate dibasic heptahydrate (99.995% trace metal basis, Sigma-Aldrich, Dorset, UK) were used for the preparation of the calibration standard and IS, respectively. Certified for i-As reference materials of white rice (ERM BC211, 124 ± 11 μg kg−1, k = 2) and brown rice (NMIJ 7532-a, 298 ± 8 μg kg−1, k = 2) were both purchased from LGC Ltd (LGC, Teddington, UK) and used as quality control samples for the demonstration of the method trueness. Baby rice food products (ground rice >99%) were purchased from local grocery stores in the UK. The moisture content of these samples was determined by drying a 0.5 g subsample in an oven at 103 ± 2 °C until constant weight, following the procedure described in the ERM BC211 certificate.
![]() | ||
Fig. 2 Overlaid chromatograms of a series of As(V) calibration standards measured with on-column IS. |
The ratio between the As(V)-analyte peak area (obtained from the calibration standard) and As(V) peak area (obtained from the internal standard; As(V)-IS) is then plotted against the calibration standard concentration and a typical calibration curve is built. The analyte concentrations are calculated from the obtained slope and intercept. Using this calibration strategy, the calibrations were linear in the range 25 ng kg−1 to 25 μg kg−1, at the instrument, with the squared Pearson coefficient (r2) exceeding 0.9998. If raw peak areas were used instead of on-column IS, the measurement accuracy would be significantly impaired due to the presence of drift. As an example, in one particular case when the sensitivity drift reached approximately 15% over the course of analysis, the linearity of the external calibration was significantly impaired and this led to lower r2 (0.998) and further biased results up to 15%.
Fig. 3 presents the signal variation with time from 50 injection pairs of 20 μL 2.5 μg kg−1 As(V) standard and 2.5 μg kg−1 As(V) internal standard represented with blue and red markers, respectively. The ratio between the analyte and IS (green triangles) was plotted against time and no statistically significant time dependence was found, although as evident from the figure the signal intensity drifted by more than 15% from its initial value, similarly to the discussed example above. The analyte to IS peak area ratio distribution had no significant deviation from the normal Gaussian profile. The slope of this ratio with time was in the order of 10−4 per injection when a linear fitting was used, which demonstrates excellent drift correction and no linear time trends. The RSD ratio was below 1% in this case. Discussion on the significance of this value can be found in the uncertainty estimation section below.
![]() | ||
Fig. 3 Efficient correction of significant sensitivity drift by the on-column internal standardisation approach. |
Similarly to SS-IDMS calibration, the value of the analyte-to-internal standard ratio can affect the accuracy of the quantitative speciation data. In this work we have investigated this effect within the range from 1:
4 to 4
:
1. The relative standard recoveries of reagent blanks spiked with As(V) standards were between 99.5% and 100.9% (ESI Fig. 1†). The latter results confirm that an exact analyte-to-IS signal matching is not essential to achieve calibration uncertainties in this range (1%), which also agrees with the findings by Yuzuru and Matsuda14 who obtained similar uncertainty using conventional IS calibration. The described order of injections of sample and IS in this work can be swapped and the IS could be injected before the sample. This would allow for mixing of the IS with the non-retained components of the sample matrix providing similar (or potentially even better) correction of the matrix effects. Additionally, this approach would be beneficial if there is a co-eluting peak at the analyte's tail but not a fronting one (resolution R > 2 requirement, see section Feasibility of application to more than one As species, below) expanding the applicability of the approach to more complex samples. Unfortunately, most of the commercial LC control software do not provide an easy option for valve control prior to the sample injection (see Fig. 1) and therefore this approach has not been further pursued.
u2(V) = u2(T) − u2(S) − u2(I) | (2) |
Eqn (2) is used for the calculation of the sample and internal standard injection volume uncertainties, where u(S) and u(I) are the standard uncertainties of the ICP-MS signal and integration, respectively, and u(T) is the combined, total uncertainty. Since the u(Dr) component is not accounted for in the equation above, it should be noted that it is still present in the calculated u(V). As for the relative u(S) and u(I), values of 0.45% and 0.09% have been obtained. These uncertainty sources are uncorrelated; therefore, covariance terms do not need to be included. Based on this, u(V) was calculated to be 0.625% and 0.627% for the sample and for the internal standard injections, respectively. These values are sufficiently low and the u(Dr) contribution can therefore be ignored.
i-As fraction [μg kg−1] | Sample/IS area ratios | Standard recoveries | Relative calibration uncertainty |
---|---|---|---|
0.1 | 0.0612 | 101.2% | 22.2% |
0.5 | 0.2886 | 99.9% | 4.8% |
0.7 | 0.4300 | 99.2% | 3.2% |
1.0 | 0.5814 | 100.0% | 2.4% |
5.0 | 2.9213 | 100.3% | 0.5% |
7.7 | 4.4796 | 99.9% | 0.3% |
10.1 | 5.8993 | 100.0% | 0.3% |
For verification purposes, 150 injections (within 3 batches with 50 injections per batch) of the 2.5 μg kg−1 As(V) standard were performed and the signal ratio for the IS injections was calculated. The relative standard deviation of the ratios (and therefore of the final concentrations) varied between 0.5% and 1.5% per batch, which is in agreement with the expected relative uncertainty between 2.4% and 0.5% for 2.5 μg kg−1 As (Table 3, above). As expected, getting closer to the 0.025 μg kg−1 limit of detection results in increased measurement uncertainty. Although the calibration approach proposed here provides accurate results where other calibration strategies might fail, it is worth noting that on-column or any other type of IS calibration is not always a necessity for performing accurate measurements. The calibration uncertainty with IS would always be higher, or at least equal to the linear regression one, when there are no drift and matrix effects to compensate for. With IS calibration two uncertainties from peak area integrations are combined for each measurement vs. only one if simple external calibration is used. With the method described in this work this effect was below 1% increase in the calibration uncertainty in case of on-column IS calibration vs. typical external calibration. However, if there is a drift present the contribution of the uncertainty of the IS peak area estimation becomes negligible compared to the improved accuracy (0.5% compared to up to 15% bias, as observed with our experiments).
Material code | A | B | C | D | E |
---|---|---|---|---|---|
Natural i-As content [μg kg−1] | 61 | 116 | 97 | 111 | 269 |
i-As spike recovery [%] | 104.6–105.4 | 99.9–101.4 | 100.2–101.8 | 103–105 | 97.8–98.7 |
A typical chromatogram of a baby food rice extract is presented in Fig. 4. To confirm the quantitative conversion of As(III) to As(V), and therefore representativeness of the As(V) to the i-As in the sample, within each batch of rice samples a random sample was spiked with As(III) before extraction. In all cases the species was quantitatively recovered as As(V), confirming the complete species oxidation within the matrix. Additionally, two matrix reference materials (white rice ERM-BC211 and brown rice NMIJ 7532a, certified for i-As) were extracted and analysed along with the rice samples. The relative i-As recoveries from both CRMs were quantitative (between 97% and 107%) and within the certified measurement uncertainties, demonstrating the accuracy of both the extraction procedure and the newly proposed calibration approach (t-test used in both cases). The relative expanded measurement uncertainty for the whole method (k = 2) was about 10% or lower in the range 60 to 300 μg kg−1 i-As in dry rice products. This includes the contribution of the materials' homogeneity and within-run and between-run uncertainty components. The on-column calibration expanded relative uncertainty was between 2% and 3% (k = 2) with ≤1.5% sample/IS peak area ratio relative uncertainty (k = 2, see the section on calibration performance above).
The applicability of the method to gradient elution chromatography is demonstrated in Fig. 6 where the same two species have been quantified with different amounts of carbon in the eluent (and therefore in ICP-MS plasma) at the time of their elution – 37 mM NH4HCO3 (eluent A) and 55 mM NH4HCO3 (eluent B), both pH adjusted to 9.0 (±0.2). The separation method involved 5 min isocratic 100% A, followed by a linear increase to 100% B in 20 min and 5 min column equilibration (100% A) before the next sample injection. The analytes [MA and As(V)] and the respective internal standards [MA IS and As(V) IS] were of the same concentration (2.5 μg kg−1 each) but in this case, as expected the species signals showed different sensitivities due to the different carbon loads in the plasma. The MA peak has smaller peak area than the MA IS peak and the As(V) peak has smaller peak area than that of As(V) IS, in contrast to observations under isocratic conditions (see Fig. 2) where sensitivities are equal. However, the analyte/IS ratio remained constant with time, e.g. throughout the batch, and the instrument drift was efficiently compensated for. Standards recoveries close to 100% (±2%) were obtained.
An additional “stress-test” was performed under the same conditions but introducing 5% MeOH in water in between samples to induce stronger drift. Although under such conditions the plasma became unstable and the peak areas increased for both species (MA and As(V)), their ratio to the respective IS species remained constant (Fig. 7a and b in the ESI†).
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d1ay00699a |
This journal is © The Royal Society of Chemistry 2021 |