Next Article in Journal
Development of a Tellurium Speciation Study Using IC-ICP-MS on Soil Samples Taken from an Area Associated with the Storage, Processing, and Recovery of Electrowaste
Next Article in Special Issue
Translating Analytical Techniques in Geochemistry to Environmental Health
Previous Article in Journal
Preclinical Assessment Addressing Intravenous Administration of a [68Ga]Ga-PSMA-617 Microemulsion: Acute In Vivo Toxicity, Tolerability, PET Imaging, and Biodistribution
Previous Article in Special Issue
Study of Long-Term Determination Accuracy for REEs in Geological Samples by Inductively Coupled Plasma Quadrupole Mass Spectrometry
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Combining Isotope Dilution and Standard Addition—Elemental Analysis in Complex Samples

1
Physikalisch-Technische Bundesanstalt (PTB), Bundesallee 100, 38116 Braunschweig, Germany
2
Bundesanstalt für Materialforschung und-Prüfung (BAM), Richard-Willstätter-Str. 11, 12489 Berlin, Germany
3
National Institute of Metrology (Thailand) (NIMT), 3/4-5 Moo 3, Klong 5, Klong Luang, Pathumthani 12120, Thailand
*
Author to whom correspondence should be addressed.
Molecules 2021, 26(9), 2649; https://doi.org/10.3390/molecules26092649
Submission received: 20 April 2021 / Revised: 23 April 2021 / Accepted: 28 April 2021 / Published: 30 April 2021

Abstract

:
A new method combining isotope dilution mass spectrometry (IDMS) and standard addition has been developed to determine the mass fractions w of different elements in complex matrices: (a) silicon in aqueous tetramethylammonium hydroxide (TMAH), (b) sulfur in biodiesel fuel, and (c) iron bound to transferrin in human serum. All measurements were carried out using inductively coupled plasma mass spectrometry (ICP–MS). The method requires the gravimetric preparation of several blends (bi)—each consisting of roughly the same masses (mx,i) of the sample solution (x) and my,i of a spike solution (y) plus different masses (mz,i) of a reference solution (z). Only these masses and the isotope ratios (Rb,i) in the blends and reference and spike solutions have to be measured. The derivation of the underlying equations based on linear regression is presented and compared to a related concept reported by Pagliano and Meija. The uncertainties achievable, e.g., in the case of the Si blank in extremely pure TMAH of urel (w(Si)) = 90% (linear regression method, this work) and urel (w(Si)) = 150% (the method reported by Pagliano and Meija) seem to suggest better applicability of the new method in practical use due to the higher robustness of regression analysis.

1. Introduction

Combined developments in inductively coupled plasma mass spectrometry (ICP–MS) techniques and instrumentation, with the continuous improvement of the calibration and evaluation of ICP–MS measurements, are currently considered the most accurate and precise tools in metrology in chemistry. This combination enables a remarkable decrease in measurement uncertainties in elemental analysis [1,2]. Related areas of analytical chemistry are geochemistry, environmental sciences, clinical chemistry, and forensics as well as fundamental metrology, to mention just a few [2,3,4,5,6,7]. Depending on sample composition (e.g., number of isotopes of the element of interest), concentration range, availability, experience of the analyst, availability of certified reference materials (CRMs), intended uncertainty, and other reasons, the determination of element content (in the following: mass fraction w) of a sample (x) can be conducted by a variety of calibration strategies [2]. Beneath simple external one-point calibration options with external standards [8], additional sophisticated methods have been developed. Among them, isotope dilution mass spectrometry (IDMS) concepts are usually the methods of choice, provided that the analyte element consists of at least two stable isotopes. The IDMS approach is considered a primary method in metrology in chemistry: it is fully understood and operated on the highest technological level. Moreover, it yields the smallest measurement uncertainties possible (most accurate and precise) according to the Guide to the Expression of Uncertainty in Measurement (GUM) [9]. Additionally, loss of the analyte after addition and equilibration of the spike does not affect the results. This method, based on the measurement of isotope ratios (intensity ratios), has been successfully applied in analytical chemistry for several decades and is still subject to further improvement [2,10,11,12,13]. Several variations of the IDMS method are known, based on the classical “single” IDMS approach. This needs a fully characterized spike material (same element as the analyte, with preferably inverted isotopic composition and known molar mass and purity) for the preparation of a blend consisting of the sample and the spike. However, often, the characterization of the spike is difficult or impossible, leading to “double” and higher-order (“triple”) IDMS approaches [2,14,15]. In “double” IDMS, a second blend is prepared from the spike and a well-characterized reference material, rendering the knowledge of purity and the molar mass of the spike obsolete. In “triple” IDMS, three blends are applied, with the additional advantage that the isotope ratio in the initial spike does not need to be measured. These higher-order approaches yield lower uncertainties at the cost of considerably more effort in sample preparation.
In the case of a monoisotopic system, the standard addition method is usually applied as an internal calibration method, being more sensitive and accurate than any external method [16,17,18]. This approach is chosen when the analyte is present in a complex matrix and/or a low concentration range. The regression analysis of the linear curve of the measured signal in several aliquots versus the ratio of the masses of the standard and sample aliquot yields in the x-axis intercept the desired mass fraction of the analyte element. A further-developed method is the standard addition method with an internal standard, useful in the case of long-term drifts and matrix effects [8,16]. This method requires the absence of the element of the internal standard in the sample, blank, and calibration standard. In most cases, the sample exists in a more or less different—and sometimes complex—matrix compared to the spike material.
Recently, Pagliano and Meija pioneered this problem in a sophisticated theoretical way by merging isotope dilution and standard addition methods [1]. They aimed at reducing or even circumventing matrix effects due to the matched matrix in all the measured solutions.
In this study, a similar and related method combining isotope dilution and standard addition methods has been developed independently. The intention of this work is to (1) reduce the impact of matrix effects in isotope dilution by combining it with a robust standard addition procedure, (2) simplify and linearize the working equations for the measurements of analyte mass fractions via linear regression to enable an uncertainty analysis according to GUM, (3) enable the determination of even trace amounts in, e.g., blank solutions, and (4) provide isotope ratio measurements with no need to correct for mass bias, which is a real benefit to the practical analyst. The method presented in the current work is demonstrated and validated experimentally using three completely different matrices in elemental analysis using ICP–MS: the determination of silicon with natural isotopic composition in an aqueous TMAH blank solution, the measurement of sulfur in biodiesel fuel (BDF), and the quantification of transferrin (TRF) in human serum.
A special need for the exact determination of the content of silicon with natural abundance—in aqueous TMAH used as a blank solution—is given in the Avogadro Project [19]. This was a multidisciplinary approach, initially set up for the realization of the Avogadro constant NA with the lowest associated measurement uncertainty and, after the SI revision in 2019, aimed at the dissemination of the SI units of kilogram and mole [19,20,21,22]. In that context, chemically highly pure silicon, extensively enriched in 28Si, with x(28Si) > 0.9999 mol/mol, was used. For the determination of the respective molar mass M of the silicon, urel(M) must be smaller than 1 × 10−8. This can only be achieved by reducing both the contamination and the quantification of the remaining natural silicon, which is described as an application of the new method presented.
The second example of the new combined method is the measurement of the mass fraction (wx) of sulfur in a biodiesel fuel matrix. These measurements were performed within an interlaboratory key comparison conducted by the IAWG of the CCQM (Inorganic Analysis Working Group of the Consultative Committee for Amount of Substance: Metrology in Chemistry and Biology): CCQM-K123 [23]. The mass fractions of several trace elements (impurities)—among them, sulfur—were determined. Originally, wx was determined using an established IDMS technique. In addition, one of the biodiesel fuel samples was measured using the new combined IDMS-Standard Addition approach as a complementary method (this work). This enables the comparability of the new approach with an established and validated method. The biodiesel fuel matrix is an extremely complex matrix because of its volatility and impurities. Sulfur occurs in different compounds in biodiesel fuel, and, for obtaining accurate results, the conversion to sulfate must be ensured in the digestion step. Additionally, sulfur contamination during the whole analytical procedure is an issue and must be controlled.
In the third part, the quantification of TRF in human serum is presented as an example for speciation analysis with a clinical background. The origin of this work was a doctoral thesis followed by the EMRP HLT05 project, which was a European metrology research program [15]. TRF is a biomarker for congenital disorders of glycosylation as well as certain cancers and alcohol abuse, and, therefore, determining its concentration in blood aids in effective clinical diagnosis [24]. Due to this, TRF is listed as an important clinical analyte in the Guidelines of the German Medical Association of 2019 (Rili-BÄK) and is routinely measured in clinical laboratories [25] even though the target values are method-dependent since no reference method is available [25]. In clinical laboratories, the common way of quantifying TRF in serum is immune-based methods like immunoturbidimetry [26].
In order to quantify the protein TRF via its iron content, high-performance liquid chromatography (HPLC) was coupled to an ICP–MS. This approach enables us to separate TRF from the complex blood matrix and identify TRF by its characteristic retention time. To this day, the quantification of TRF is performed by applying double or triple IDMS and considering the iron background and complex matrix of the samples [15,27]. The new approach presented in this publication, combining IDMS with standard addition, is used for the first time in a complex biological sample as proof of principle. Since the sample is a certified serum material, the results of the new approach are compared with the certified target value and its reported range.
A central part of this work is also the comparison of the new method presented here with the one reported in [1]. Both methods have the same or a similar intention: the improvement of measurement design and the reduction of uncertainty for the determination of mass fractions wx of elements in complex matrices. For better understanding and adaption to our notation, we have additionally derived the method reported in [1] from scratch for the three blends case: we ended up with exactly the same equation (which is Equation (9) in [1]). For the comparison, we evaluated all measurements using our new method (this work) and the one reported in [1].

2. Theoretical Methods

The main advantage of the combination of IDMS and standard addition is its unique simplicity of using gravimetrically prepared blends, ensuring traceability to the International System of Units (SI). It is based on the measurement of virtually only one isotope ratio and three masses per blend. The new approach yields, finally, Equation (1) to describe linear regression. The regression parameters, slope a1 and y-intercept a0, are then used to calculate the aimed-at mass fraction wx of the analyte element in the original sample (x). The detailed derivation of Equations (1) and (2) is given in Appendix A.
m y , i m x , i × R y R b , i R b , i R x = y i = 1 w y M y M x R y R x × w z = a 1 × m z , i m x , i = x i + 1 w y M y M x R y R x × w x = a 0
w x = a 0 a 1 × w z
where wx is directly proportional to the mass fraction wz of a reference material (z; same element as the sample with the same isotopic composition as the analyte or, at least, very close to it); my,i, mx,i, and mz,i are the masses of a spike solution (y; same element as x with a preferably inverse isotopic composition), sample x, and reference material solution z, respectively. A number of blends (bi; e.g., five) containing almost the same masses of x and y, respectively, and different masses of z (mz,1 < mz,2 < mz,3 < mz,4 < mz,5) are required for the experiment and evaluation (Figure 1 and Figure 2). Ry, Rx, and Rb,i, denote isotope ratios measured in the respective spike, sample, and ith blend solution related to the isotope of major abundance (reference isotope; in the case of enriched silicon: 28Si; Rx = I(30Si)/I(28Si), with the intensities I of the measured isotope signals); Mx and My are the molar masses of the respective analyte and spike material, which need not be known. The only requirement is the knowledge of mass fraction wz of the element of interest in reference material z. This is given in a calibration certificate or accessible via a separate measurement (gravimetric preparation).
One advantage of this new method is based on the fact that the measured intensity ratios are entered into Equation (1) directly without the need for mass bias correction. Usually, in ICP–MS, measured isotope ratios suffer from mass discrimination effects (mass bias) and have to be corrected by the use of calibration K factors. The latter are accessible via reference materials, applications of empirical relations like the exponential law, or gravimetrical methods [3,5,28]. In the set of equations used in this work, however, the K factors applied to measured isotope ratios cancel out, enabling the direct use of the measured ratios and, thus, simplifying this method for convenience. The only exception is spike ratio Ry, for which no value measured in the matrix is available. Reasonable spike materials feature Ry > 100. In this case, mixing ratio Ry, measured in a matrix different from that of the sample matrix, with Rx and Rb,i measured in the sample matrix, changes the result (wx) well within its uncertainty and, therefore, is insignificant.
As stated above, Equation (9) in [1] has been analytically derived again in this work from scratch. This derivation is given in the Supplementary Materials. Equation (3) displays the analytical solution for mass fraction wx in the notation used in this work. This equation is exactly the same as Equation (9) in [1]. Additionally, mass fraction wz of the reference material used has to be known. Masses (mj,i) of components j (j = sample x, reference z, or spike y) in blends i (i = 1, 2, 3) have to be known (gravimetric preparation) as well as the respective isotope ratios, namely, the Rb3 reads:isotope ratio (monitor vs. reference isotope) in Blend 3, the Rz,2 reads:isotope ratio (monitor vs. reference isotope) in reference z, and xz,1 reads:amount-of-substance fraction x of Isotope 1 (reference isotope) in reference material z. Note that the “reference” isotope is the isotope of the highest abundance in the sample, whereas “reference” material is a characterized material with (almost) the same isotopic composition as the sample. In the case of different isotopic compositions of x and z, the molar masses, Mx (sample) and Mz (reference), and the amount-of-substance fractions, xz,1 and xx,1, have to be known. When applying Equation (3), three gravimetrically prepared blends are required.
w x = w z m y 1 m z 2 m y 3 R b 2 R z , 2 R b 3 R b 1 m y 1 m y 2 m z 3 R b 3 R z , 2 R b 2 R b 1 m z 1 m y 2 m y 3 R b 1 R z , 2 R b 3 R b 2 m y 1 m x 2 m y 3 R b 2 R x , 2 R b 3 R b 1 + m y 1 m y 2 m x 3 R b 3 R x , 2 R b 2 R b 1 + m x 1 m y 2 m y 3 R b 1 R x , 2 R b 3 R b 2 x z , 1 M x x x , 1 M z

3. Materials and Experimental Methods

In the first test system, the sample solutions, x, consisted of aqueous TMAH (w(TMAH) = 0.0006 g/g) as the matrix with traces of the analyte: silicon with natural isotopic composition. The sample solutions were taken from a study described in [29]. There, the aqueous TMAH served as a blank solution for the measurement of isotope ratios of enriched silicon samples. Briefly, the TMAH sample solutions were prepared by dissolving highly concentrated (w(TMAH) = 0.25 g/g), electronic grade (99.9999%) TMAH from Alfa Aesar™ (Thermo Fisher (Kandel) GmbH, Germany) in highly purified water (resistivity = 18 MΩ cm), the latter prepared using a commercial water purification system (Merck Millipore™, Burlington, MA, USA). All solutions were prepared and stored in precleaned labware (perfluoroalkoxy copolymer PFA) according to a protocol given in [29]. The spike solutions, y, originated from chemically pure silicon single crystals highly enriched in 30Si dissolved in aqueous TMAH, whereas reference material z was prepared from well-characterized silicon crystals (material code: WASO04) with known natural-like isotopic composition [30]. All stock solutions, dilutions, and blends in this study were prepared gravimetrically, applying air buoyancy correction [31]. The silicon crystals used for the stock solutions of y (spike) and z (reference) were cleaned and etched prior to dissolution with the final mass fractions prior to blend preparation: wz = 4 µg/g; wy = 2 µg/g. The isotopic composition of y as well as the respective “true” isotope ratios, Ry, were taken from [20]. For the determination of wx according to Equations (1) and (2), a series of five blends, bi, was prepared, each consisting of approximately mx,i = 10 g, my,i = 22.5 g, and stepwise increasing amounts of solution z with mz,1 = 0 g, mz,2 = 7.5 g, mz,3 = 10 g, mz,4 = 14.5 g, and mz,5 = 23 g (compare Figure 1).
The sample preparation for sulfur determination in the biodiesel fuel matrix is described in [23]. Briefly, an ampoule (15 mL) from the CCQM-K123 study containing commercial biodiesel fuel doped with the target analyte sulfur using dibutyl sulfide was used as the sample. Defined aliquots of reference z (NIST SRM 3154) and spike y (BAM S-34) material were added prior to the digestion. For the blend solutions, acid digestion was applied (HNO3/H2O2) using a high-pressure asher system (Anton Paar, Graz, Austria), with ϑmax ≈ 300 °C and pmax ≈ 130 bar. After digestion, the solution was evaporated to dryness, redissolved with 2 mL HNO3 (0.028 mol/L), and loaded on chromatographic columns filled with 1 mL AG 1X8 resin. The sample matrix was eluted with water, and, subsequently, the sulfur was eluted with 8 mL HNO3 (0.25 mol/L). Samples were evaporated to dryness and redissolved in HNO3 (0.02 kg/kg) to adjust a sulfur mass fraction of 2 mg/kg. Prior to the MC–ICP–MS measurements, sodium was added so that a final Na mass fraction of 4 mg/kg was achieved. The addition of sodium significantly increases sensitivity and prevents losses of sulfur through the membrane of the desolvating nebulizer system, presumably as sulfuric acid [32]. The complete analytical procedure was validated by applying the double IDMS calibration approach with reference material NIST SRM 2723a. The determined sulfur mass fraction was (10.94 ± 0.10) mg/kg, while the certified sulfur mass fraction was (10.90 ± 0.31) mg/kg, with k = 2 in both cases.
For the measurement of the third system, TRF in human serum, reference material ERM®-DA470k/IFCC “Human Serum” (Merck KGaA, Darmstadt, Germany) was applied as the reference, and SeronormTM Immunoprotein Lyo L-1 (Invicon GmbH, Munich, Germany), which has certified values and ranges in accordance with the Guidelines of the German Medical Association (Rili-BÄK) of 2014, was used as the sample [33]. Both materials were reconstituted according to the manufacturers’ protocols and thereafter used for sample preparation [33]. The iron saturation procedure was based on the method described by del Castillo Busto et al. and C. Frank et al. except for some optimizations [15,34]. For the individual blends, a solution containing 7.5 to 22.5 mg of an iron solution (250 µg/g in 2.5% HNO3, dilution from BAM A-primary-Fe-2; BAM, Berlin, Germany), 5 to 10 mg sodium carbonate (Na2CO3) solution (500 mmol/kg; BioUltra, Merck KGaA, Darmstadt, Germany), and 50 to 277 mg tris(hydroxymethyl)methylamine (Tris) buffer solution (12.5 mmol/kg, adjusted to pH 6.4 with acetic acid; BioUltra, Merck KGaA, Darmstadt, Germany) was mixed, and 0 to 200 mg of the reconstituted ERM®-DA470k/IFCC serum and 50 mg of the reconstituted SeronormTM Immunoprotein Lyo L-1 serum were added. The amount of iron and Na2CO3 solution was adjusted to the amount of added human serum, which increased from Blend 1 to Blend 4. All solutions were incubated at room temperature for 1.5 h. Then, 160 mg of the TRF spike was added. In Section 4.3, the added amounts of reference, sample, and spike solution for the different blends are given.
The TRF spike was prepared in advance from human apo-TRF purchased from Sigma-Aldrich (St Louis, MO, USA) and 54Fe spike from Trace Science International Corp. (Ontario, ON, Canada), with a certified isotopic abundance of 54Fe 99.86%, 56Fe 0.11%, 57Fe 0.01%, and 58Fe 0.02%. For this, 20 mg of apo TRF was dissolved in a solution containing 0.435 g Na2CO3 (500 mmol/kg), 4.3 g Tris (12.5 mmol/kg), and 0.4 g 54Fe (250 µg/g), adjusted to pH = 8 with 0.15 mmol/kg HNO3. After an incubation time of 3 h at room temperature, the solution was purified with a PD-10 desalting column (Cytiva GmbH, Freiburg im Breisgau, Germany), according to the manufacturer’s protocol. The spike solution was freeze-dried, and the solid was used for the preparation of a TRF spike solution of 2.3 mg/g, which was used to prepare the blends.
After the gravimetric sample and blend preparation, in all blends, lipoprotein precipitation was performed by adding 5 µL of a solution containing 200 mg magnesium chloride (Merck KGaA, Darmstadt, Germany) and 100 mg dextran sulfate sodium salt (Merck KGaA, Darmstadt, Germany) in 1 mL ultrapure water, incubating for 30 min at 4 °C and centrifuging the samples at 12,000× g for 10 min. Each supernatant was purified with PD MidiTrap columns (Cytiva GmbH, Freiburg im Breisgau, Germany), according to the manufacturer’s protocol. The solutions obtained were used for the HPLC/ICP–MS analysis.
All three sample/matrix systems were investigated using ICP–MS instrumentation. The respective operation parameters of the isotope ratio measurements are summarized in Table 1. As an example, the Si in TMAH measurements were run in six sequences measuring the intensities of all three stable Si isotopes (28Si, 29Si, and 30Si) in each blend, bi (i = 1…5).

4. Results and Discussion

4.1. Silicon in Aqueous TMAH

The first example of the application of the new combined IDMS-Standard Addition approach is the determination of silicon with natural isotopic composition in the matrix of aqueous TMAH. As an example, the essential input and output quantities of one experimental run (Sequence 1) are given in Table 2. The corresponding linear regression curve is displayed in Figure 3.
As a direct consequence of linear regression analysis, mass fraction wx of the analyte in the respective sample x (TMAHaq) was derived according to Equation (2). A set of 5 blends is sufficient to obtain proper regression statistics.
The individual results of wx(Si), obtained from six sequences (runs), are displayed in Figure 4 (left) and Table 3. They yield an average of wx(Si) = 0.081(73) µg/g. As proof of consistency of the individual measurement results, the concept of degrees of equivalence (di) was applied. The individual di of the respective measurements are plotted in Figure 4 (right). All data encompass the zero line with their associated expanded uncertainties, which is an approval of the complete consistency of the data.

4.2. Sulfur in Biodiesel Fuel (BDF)

The determination of sulfur in BDF, according to Equations (1) and (2), was performed in a similar way using the five blends, bi. The essential input and output quantities of one experimental run are given in Table 4. The corresponding linear regression curve of the respective dataset is displayed in Figure 5.
Linear regression analysis yields the mass fraction (wx) of the analyte (sulfur) in the respective sample x (BDF), according to Equation (2).
The determination of sulfur in BDF, according to Equations (1) and (2), yields an average mass fraction of wx,corr(S) = 7.36(11) µg/g. The individual results of wx,corr(S), measured in three runs, are shown in Figure 6 (left) and Table 5. Additionally, as proof of consistency, the concept of degrees of equivalence, di, was applied. The di of the respective measurements are plotted in Figure 6 (right). The single results encompass the zero line with their associated expanded uncertainties, which proves the complete consistency of the data. In the case of sulfur determination in this work, the respective sample was taken from a stock solution already used in the interlaboratory comparison CCQM-K123, as described in [23]. There, the sulfur determination was carried out using double IDMS, yielding wx = 7.394(46) µg/g. The corresponding result is shown in Figure 6 (left, red solid line) for comparison. Both the IDMS and the new combined IDMS/standard addition approach show excellent agreement within the limits of uncertainty.

4.3. Transferrin in Human Serum

The determination of TRF in human serum, according to Equations (1) and (2), was performed in a similar way using four blends, bi. The essential input and output quantities of one experimental run are given in Table 6. The corresponding linear regression curve of the respective dataset is displayed in Figure 7.
The linear regression analysis yields mass fraction wx of the analyte TRF in the respective sample x (human serum) according to Equation (2).
The individual results of wx(TRF), obtained from two independent measurement sequences, are displayed in Figure 8 (left). Both values, wx(TRF) = 2.128(45) mg/g and wx(TRF) = 2.340(76) mg/g for k = 1, are consistent with the certified target value and its allowed range [33]. As proof of consistency of the individual measurement results, the concept of degrees of equivalence, di, was applied. The individual di of the respective measurements are plotted in Figure 8 (right). All data encompass the zero line with their associated expanded uncertainties.

4.4. Comparison of Linear Regression (This Work) and Analytical Solution (Pagliano and Meija)

The comparison of the two approaches, combining IDMS and standard addition presented in this work and in [1], was performed using the dataset of silicon determination in the TMAH matrix (five available blends: b1, b2, b3, b4, b5). Figure 9 shows the results of wx(Si), determined in six sequences (runs), using five blends (this work) and three blends (b2, b3, b4; approach of [1] using Equation (3)). At first glance, the associated uncertainties of the results determined in this work (black circles) are significantly smaller than the respective uncertainties of the results determined using the approach of [1] (red circles).
Prior to this comparison, both approaches (this work and that of [1]) were checked for consistency using simulated ideal data: both methods yielded exactly the same result. However, when using real experimental data (e.g., the results of the five blends for silicon determination), the results of the approach of [1] depend strongly on the input data (the three blends that were used for analysis). Table 7 shows the different numerical results between the two approaches of a representative run. In the case of this work, all five blends were used (left column), whereas, for the approach of [1], triple blends b1, b2, b3 or b2, b3, b4, or b3, b4, b5 each yielded extremely different results. The analytical approach (Equation (9) of [1], which is equal to Equation (3) in this publication) is based on a set of three blends only, without the need for the measurement of Ry, which appears initially as a benefit. However, the analysis shows that the three-blend application leads, in some cases, to unrealistic (not accurate) results (results of triples b1, b2, b3 and b3, b4, b5). Equation (9) in [1] reacts extremely sensitively towards that kind of input data. In our case, the three inner blends yield the most reasonable (“accurate”) result, e.g., the first blend, b1, includes no reference material z (mz = 0 g). This suggests a rather unstable and sensitive applicability of Equation (9) of [1] compared to the linear regression approach of this work. Moreover, for practical applications, Equations (1) and (2) are much simpler and user-friendly, and the respective evaluations and measurements are much easier than in the case of the approach used in [1], although the latter is correct from the mathematical point of view.
For the comparison of the uncertainties associated with wx, which is the second-most important outcome of the respective approach, we chose the results of the blend triple b2, b3, b4 of the six sequences. As can be seen in Figure 9, the uncertainties obtained using the approach of this work are significantly smaller than those obtained using Equation (9) from [1].
An uncertainty analysis (Table 8) of Equation (3) (same as Equation (9) in [1]) according to GUM [9], using GUM Workbench Pro™ software (version 2.4.1 392, Metrodata GmbH, Germany), shows that the absolute values of the sensitivity coefficients of the main uncertainty contributions, Rb2, Rb1, and Rb3, range in the 101-102 region. This seems rather high and is mainly responsible for the elevated uncertainty compared to the smaller uncertainty obtained using linear regression in the approach of this work.

5. Conclusions

The new method presented—a combination of IDMS and standard addition—for absolute elemental quantification is recommended for sample analytes with two or more isotopes, especially when appearing in traces (ng/g) and/or in a complex matrix that is difficult or (almost) impossible to separate. The three different examples of this study demonstrate the convenient applicability of this approach successfully. Moreover, the second example (determination of sulfur in biodiesel fuel) acts as a validation system since the same sample material has been measured previously and independently using a double IDMS method in the context of an international key comparison (CCQM-K123), yielding excellently matching results within the limits of uncertainty. Mainly, since the blends are measured, all the “measurement” solutions are matrix-matched; thus, potential matrix effects are also canceled out. The comparison of the new method (linear regression) with the analytical approach reported in [1] shows the formal (mathematical) agreement of both methods. However, the use of experimental data suggests a more robust behavior of the approach presented in this work, yielding more “accurate” results with significantly lower associated uncertainties (see calculation examples in the Supplement). Moreover, the new approach is more practical in daily laboratory work due to a simpler evaluation algorithm that might be less fault-prone.

Supplementary Materials

The following are available online, Supplement_a_ID_STDAdd.docx: with Table S1: List of quantities and symbols in the general final Equation (3) and Supplement-b-StdAdd-IDMS.xlsm.

Author Contributions

Conceptualization, A.P. and O.R.; methodology, O.R.; validation, A.P. and J.V.; formal analysis, C.B., A.P., O.R., A.S., P.P., and J.V.; investigation, C.B., A.P., O.R., and J.V.; data curation, C.B., A.P., O.R., and J.V.; writing—original draft preparation, A.P.; writing—review and editing, C.B., A.P., O.R., and J.V. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data that support the findings of this study are available from the authors upon request.

Acknowledgments

The authors gratefully acknowledge animated discussions with Detlef Schiel (PTB). Special thanks to Dorit Becker and Maren Koenig for their solid support with the sulfur measurements.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Derivation of Combined IDMS and Standard Addition Method
Isotope ratios (Rb,i) in the blends (i) containing the analyte can be expressed as the ratio of the sums of the amount-of-substances (n) of the respective components of x (sample), y (spike), and z (reference). In our notation, Subscript 1 denotes the reference and Subscript 2 the monitor (spike) isotope.
R b , i = n x 2 , i + n z 2 , i + n y 2 , i n x 1 , i + n z 1 , i + n y 1 , i
R b , i = x x 2 n x , i + x z 2 n z , i + x y 2 n y , i x x 1 n x , i + x z 1 n z , i + x y 1 n y , i
R b , i = x x 2 m x , i w x M x + x z 2 m z , i w z M z + x y 2 m y , i w y M y x x 1 m x , i w x M x + x z 1 m z , i w z M z + x y 1 m y , i w y M y
In case of equal isotopic patterns in the sample and reference (the equation, valid in the case of different isotopic patterns, including its derivation, is presented in the Supplement):
x x 1 = x z 1 , x x 2 = x z 2 , M x = M z
Therefore, Rb,i can be rewritten accordingly
R b , i = x x 2 M x m x , i w x + m z , i w z + x y 2 M y m y , i w y x x 1 M x m x , i w x + m z , i w z + x y 1 M y m y , i w y
After rearranging and expanding
R b , i x x 1 M x m x , i w x + m z , i w z + x y 1 M y m y , i w y = x x 2 M x m x , i w x + m z , i w z + x y 2 M y m y , i w y
R b , i x x 1 M x m x , i w x + m z , i w z + R b , i x y 1 M y m y , i w y = x x 2 M x m x , i w x + m z , i w z + x y 2 M y m y , i w y
m x , i w x + m z , i w z 1 M x R b , i x x 1 x x 2 = m y , i w y 1 M y x y 2 R b , i x y 1
m y , i w y 1 M y x y 2 R b , i x y 1 = m x , i w x + m z , i w z 1 M x R b , i x x 1 x x 2
m y , i x y 2 R b , i x y 1 R b , i x x 1 x x 2 = 1 w y M y M x m x , i w x + m z , i w z
m y , i x y 2 R b , i x y 1 R b , i x x 1 x x 2 = 1 w y M y M x m x , i w x + 1 w y M y M x m z , i w z
m y , i m x , i x y 2 R b , i x y 1 R b , i x x 1 x x 2 = 1 w y M y M x m z , i w z m x , i + 1 w y M y M x w x
m y , i m x , i R y 2 R y R b , i R y 1 R y R b , i R x 1 R x R x 2 R x = 1 w y M y M x m z , i w z m x , i + 1 w y M y M x w x with   R y 1 = R x 1 = 1
m y , i m x , i 1 R y 1 R x R y 2 R b , i × 1 R b , i × 1 R x 2 = 1 w y M y M x m z , i w z m x , i + 1 w y M y M x w x
m y , i m x , i R x R y R y 2 R b , i × 1 R b , i × 1 R x 2 = 1 w y M y M x m z , i w z m x , i + 1 w y M y M x w x
Equation (1) already represents the linear equation required to describe the regression line.
m y , i m x , i × R y R b , i R b , i R x = y i = 1 w y M y M x R y R x × w z = a 1 × m z , i m x , i = x i + 1 w y M y M x R y R x × w x = a 0
The y-axis intercept a0 and slope a1 yield the mass fraction (wx) of the analyte element in the sample
w x = a 0 a 1 × w z

References

  1. Pagliano, E.; Meija, J. Reducing the matrix effects in chemical analysis: Fusion of isotope dilution and standard addition methods. Metrologia 2016, 53, 829–834. [Google Scholar] [CrossRef] [Green Version]
  2. Sargent, M.; Goenaga-Infante, H.; Inagaki, K.; Ma, L.; Meija, J.; Pramann, A.; Rienitz, O.; Sturgeon, R.; Vogl, J.; Wang, J.; et al. The role of ICP-MS in inorganic chemical metrology. Metrologia 2019, 56, 1–32. [Google Scholar] [CrossRef]
  3. Vanhaecke, F.; Degryse, P. MC-ICP-MS in Isotopic Analysis: Fundamentals and Applications Using ICP-MS; Wiley-VCH: Weinheim, Germany, 2012. [Google Scholar]
  4. Yang, L. Accurate and precise determination of isotopic ratios by MC-ICP-MS: A review. Mass Spectrom. Rev. 2009, 28, 990–1011. [Google Scholar] [CrossRef] [PubMed]
  5. Yang, L.; Tong, S.; Zhou, Z.; Hu, Z.; Mester, Z.; Meija, J. A critical review on isotopic fractionation correction methods for accurate isotope amount ratio measurements by MC-ICP-MS. J. Anal. At. Spectrom. 2018, 33, 1849–1861. [Google Scholar] [CrossRef]
  6. Andrén, H.; Rodushkin, I.; Stenberg, A.; Malinovskiy, D.; Baxter, D.C. Sources of mass bias and isotope ratio variation in multi-collector ICP-MS: Optimization of instrumental parameters based on experimental observations. J. Anal. At. Spectrom. 2004, 19, 1217–1224. [Google Scholar] [CrossRef]
  7. Pramann, A.; Rienitz, O.; Schiel, D.; Güttler, B.; Valkiers, S. Novel concept for the mass spectrometric determination of absolute isotopic abundances with improved measurement uncertainty: Part 3—Molar mass of silicon highly enriched in 28Si. Int. J. Mass Spectrom. 2011, 305, 58–68. [Google Scholar] [CrossRef]
  8. Vogl, J. Calibration strategies and quality assurance. In Inductively Coupled Plasma Mass Spectrometry Handbook; Nelms, S.M., Ed.; Blackwell Publishing Ltd: Oxford, UK, 2005; pp. 147–181. [Google Scholar]
  9. Working Group 1 of the Joint Committee for Guides in Metrology (JCGM/WG 1). JCGM 100:2008. Evaluation of Measurement Data—Guide to the Expression of Uncertainty in Measurement; JCGM: Sevres, France, 2008. [Google Scholar]
  10. Heumann, K.G. Isotope-dilution mass spectrometry of inorganic and organic substances. Fresenius Z. Anal. Chem. 1986, 325, 661–666. [Google Scholar] [CrossRef]
  11. De Bièvre, P. Isotope dilution mass spectrometry: What can it contribute to accuracy in trace analysis? Fresenius Z. Anal. Chem. 1990, 327, 766–771. [Google Scholar] [CrossRef]
  12. Heumann, K.G. Isotope dilution mass spectrometry (IDMS) of the elements. Mass Spectrom. Rev. 1992, 11, 41–67. [Google Scholar] [CrossRef]
  13. Sargent, M.; Harrington, C.; Harte, R. Guidelines for Achieving High Accuracy in Isotope Dilution Mass Spectrometry (IDMS); LGC Ltd: Cambridge, MA, USA, 2002. [Google Scholar]
  14. Vogl, J. Measurement uncertainty in single, double, and triple isotope dilution mass spectrometry. Rapid Comm. Mass Spectrom. 2012, 26, 275–281. [Google Scholar] [CrossRef]
  15. Frank, C.; Rienitz, O.; Swart, C.; Schiel, D. Improving species-specific IDMS: The advantages of triple IDMS. Anal. Bioanal. Chem. 2013, 405, 1913–1919. [Google Scholar] [CrossRef]
  16. Hauswaldt, A.-L.; Rienitz, O.; Jährling, R.; Fischer, N.; Schiel, D.; Labarraque, G.; Magnusson, B. Uncertainty of standard addition experiments: A novel approach to include the uncertainty associated with the standard in the model equation. Accred. Qual. Assur. 2012, 17, 129–138. [Google Scholar] [CrossRef]
  17. Rienitz, O.; Hauswaldt, A.-L.; Jährling, R. Standard addition challenge. Anal. Bioanal. Chem. 2012, 403, 2461–2462. [Google Scholar] [CrossRef] [PubMed]
  18. Rienitz, O.; Hauswaldt, A.-L.; Jährling, R. Solution to standard addition challenge. Anal. Bioanal. Chem. 2012, 404, 2117–2118. [Google Scholar] [CrossRef] [PubMed]
  19. Fujii, K.; Bettin, H.; Becker, P.; Massa, E.; Rienitz, O.; Pramann, A.; Nicolaus, A.; Kuramoto, N.; Busch, I.; Borys, M. Realization of the kilogram by the XRCD method. Metrologia 2016, 53, A19–A45. [Google Scholar] [CrossRef]
  20. Pramann, A.; Lee, K.S.; Noordmann, J.; Rienitz, O. Probing the homogeneity of the isotopic composition and molar mass of the ‘Avogadro’-crystal. Metrologia 2015, 52, 800–810. [Google Scholar] [CrossRef]
  21. Bartl, G.; Becker, P.; Beckhoff, B.; Bettin, H.; Beyer, E.; Borys, M.; Busch, I.; Cibik, L.; D’Agostino, G.; Darlatt, E.; et al. A new 28Si single crystal: Counting the atoms for the new kilogram definition. Metrologia 2017, 54, 693–715. [Google Scholar] [CrossRef]
  22. Güttler, B.; Rienitz, O.; Pramann, A. The avogadro constant for the definition and realization of the mole. Ann. Phys. 2018, 1800292. [Google Scholar] [CrossRef] [Green Version]
  23. Kuroiwa, T.; Zhu, Y.; Inagaki, K.; Long, S.E.; Christopher, S.J.; Puelles, M.; Borinsky, M.; Hatamleh, N.; Murby, J.; Merrick, J.; et al. Report of the CCQM-K123: Trace elements in biodiesel fuel. Metrologia 2017, 54, 08008. [Google Scholar] [CrossRef]
  24. Jou, Y.J.; Lin, C.D.; Lai, C.H.; Chen, C.H.; Kao, J.Y.; Chen, S.Y.; Tsai, M.H.; Huang, S.H.; Lin, C.W. Proteomic identification of salivary transferrin as a biomarker for early detection of oral cancer. Anal. Chim. Acta 2010, 681, 41–48. [Google Scholar] [CrossRef]
  25. Deutsches Ärzteblatt. Neufassung der Richtlinie der Bundesärztekammer zur Qualitätssicherung laboratoriumsmedizinischer Untersuchungen—Rili-BÄK. Dtsch. Arztebl. 2019, 116, A-2422 / B-1990 / C-1930. [Google Scholar]
  26. Feng, L.; Zhang, D.; Wang, J.; Li, H. Simultaneous quantification of proteins in human serum via sulfur and iron using HPLC coupled to post-column isotope dilution mass spectrometry. Anal. Methods 2014, 6, 7655–7662. [Google Scholar] [CrossRef]
  27. Larios, R.; Del Castillo Busto, M.E.; Garcia-Sar, D.; Ward-Deitrich, C.; Goenaga-Infante, H. Accurate quantification of carboplatin adducts with serum proteins by monolithic chromatography coupled to ICPMS with isotope dilution analysis. J. Anal. At. Spectrom. 2019, 34, 729–740. [Google Scholar] [CrossRef]
  28. Mana, G.; Rienitz, O. The calibration of Si-isotope ratio measurements. Int. J. Mass Spectrom. 2010, 291, 55–60. [Google Scholar] [CrossRef]
  29. Rienitz, O.; Pramann, A.; Vogl, J.; Lee, K.-S.; Yim, Y.-H.; Malinovskiy, D.; Hill, S.; Dunn, P.; Goenaga-Infante, H.; Ren, T.; et al. The comparability of the determination of the molar mass of silicon highly enriched in 28Si: Results of the CCQM-P160 interlaboratory comparison and additional external measurements. Metrologia 2020, 57, 1–13. [Google Scholar] [CrossRef]
  30. D’Agostino, G.; Di Luzio, M.; Mana, G.; Oddone, M.; Pramann, A.; Prata, M. 30Si mole fraction of a silicon material highly enriched in 28si determined by instrumental neutron activation analysis. Anal. Chem. 2015, 87, 5716–5722. [Google Scholar] [CrossRef] [Green Version]
  31. Kaltenbach, A.; Noordmann, J.; Görlitz, V.; Pape, C.; Richter, S.; Kipphardt, H.; Kopp, G.; Jährling, R.; Rienitz, O.; Güttler, B. Gravimetric preparation and characterization of primary reference solutions of molybdenum and rhodium. Anal. Bioanal. Chem. 2015, 407, 3093–3102. [Google Scholar] [CrossRef]
  32. Paris, G.; Sessions, A.L.; Subhas, A.V.; Adkins, J.F. MC-ICP-MS measurement of δ34S and ∆33S in small amounts of dissolved sulfate. Chem. Geol. 2013, 345, 50–61. [Google Scholar] [CrossRef]
  33. Deutsches Ärzteblatt. Neufassung der Richtlinie der Bundesärztekammer zur Qualitätssicherung laboratoriumsmedizinischer Untersuchungen—Rili-BÄK. Dtsch. Arztebl. 2014, 111, A-1583/B-1363/C-1295. [Google Scholar]
  34. Del Castillo Busto, M.E.; Montes-Bayón, M.; Sans-Mendel, A. Accurate determination of human serum transferrin isoforms: Exploring metal-specific isotope dilution analysis as a quantitative proteomic tool. Anal. Chem. 2006, 78, 8218–8226. [Google Scholar] [CrossRef]
Figure 1. Schematic of a set of 5 blends (bi) gravimetrically prepared from approximately the same masses (mx,i) of analyte sample x, same masses (my,i) of spike solution y, and different masses (mz,i) of reference solution z (mz,1 < mz,2 < mz,3 < mz,4 < mz,5; in this work, mz,1 = 0 g).
Figure 1. Schematic of a set of 5 blends (bi) gravimetrically prepared from approximately the same masses (mx,i) of analyte sample x, same masses (my,i) of spike solution y, and different masses (mz,i) of reference solution z (mz,1 < mz,2 < mz,3 < mz,4 < mz,5; in this work, mz,1 = 0 g).
Molecules 26 02649 g001
Figure 2. Relationship of the initial components (sample x, spike y, and reference z) and the respective blends (bi), indicating (a) the abundances in which the components were blended and (b) the quantities, which have to be measured to be able to calculate analyte mass fraction wx in the sample.
Figure 2. Relationship of the initial components (sample x, spike y, and reference z) and the respective blends (bi), indicating (a) the abundances in which the components were blended and (b) the quantities, which have to be measured to be able to calculate analyte mass fraction wx in the sample.
Molecules 26 02649 g002
Figure 3. Linear regression evaluation of the Si mass fraction in TMAHaq according to Equation (2). Dataset in Table 2.
Figure 3. Linear regression evaluation of the Si mass fraction in TMAHaq according to Equation (2). Dataset in Table 2.
Molecules 26 02649 g003
Figure 4. (Left): Mass fraction wx (Si) in TMAHaq. Error bars denote combined uncertainties (k = 1). The red dashed line indicates the average value. (Right): Degrees of equivalence di of the respective measurement results. Error bars indicate expanded uncertainties (k = 2) associated with di. All single results are consistent with the average value since the respective uncertainties encompass the red dashed zero line.
Figure 4. (Left): Mass fraction wx (Si) in TMAHaq. Error bars denote combined uncertainties (k = 1). The red dashed line indicates the average value. (Right): Degrees of equivalence di of the respective measurement results. Error bars indicate expanded uncertainties (k = 2) associated with di. All single results are consistent with the average value since the respective uncertainties encompass the red dashed zero line.
Molecules 26 02649 g004
Figure 5. Linear regression of the S mass fraction in biodiesel fuel according to Equation (2). Dataset in Table 4.
Figure 5. Linear regression of the S mass fraction in biodiesel fuel according to Equation (2). Dataset in Table 4.
Molecules 26 02649 g005
Figure 6. (Left) Black circles: corrected mass fractions wx,corr(S) of sulfur in BDF and the associated uncertainties (k = 1); black dashed line: average of three single runs. Data from the new combined IDMS/standard addition approach (this work). Red solid line: average value of w(S) from BAM, applying IDMS and using the same solution [23]. Upper and lower associated uncertainties: red dotted lines. (Right) Degrees of equivalence (di) of the respective measurement results. Error bars indicate expanded uncertainties (k = 2) associated with di. All single results are consistent with the average value since the respective uncertainties encompass the red dashed zero line.
Figure 6. (Left) Black circles: corrected mass fractions wx,corr(S) of sulfur in BDF and the associated uncertainties (k = 1); black dashed line: average of three single runs. Data from the new combined IDMS/standard addition approach (this work). Red solid line: average value of w(S) from BAM, applying IDMS and using the same solution [23]. Upper and lower associated uncertainties: red dotted lines. (Right) Degrees of equivalence (di) of the respective measurement results. Error bars indicate expanded uncertainties (k = 2) associated with di. All single results are consistent with the average value since the respective uncertainties encompass the red dashed zero line.
Molecules 26 02649 g006
Figure 7. Linear regression evaluation of the TRF mass fraction in human serum according to Equation (2). Dataset in Table 6.
Figure 7. Linear regression evaluation of the TRF mass fraction in human serum according to Equation (2). Dataset in Table 6.
Molecules 26 02649 g007
Figure 8. (Left) Black circles: mass fractions wx(TRF) of transferrin in human serum and the associated uncertainties (k = 1); dashed black line: certified value of the Seronorm sample (for instrument, see Beckmann, AU). Data from the new combined IDMS/standard addition approach (this work). Red dashed lines: allowed target range according to Rili-BÄK [33]. (Right) Degrees of equivalence, di, of the respective measurement results. Error bars indicate expanded uncertainties (k = 2) associated with di. All single results are consistent with the certified value since the respective uncertainties encompass the red dashed zero line.
Figure 8. (Left) Black circles: mass fractions wx(TRF) of transferrin in human serum and the associated uncertainties (k = 1); dashed black line: certified value of the Seronorm sample (for instrument, see Beckmann, AU). Data from the new combined IDMS/standard addition approach (this work). Red dashed lines: allowed target range according to Rili-BÄK [33]. (Right) Degrees of equivalence, di, of the respective measurement results. Error bars indicate expanded uncertainties (k = 2) associated with di. All single results are consistent with the certified value since the respective uncertainties encompass the red dashed zero line.
Molecules 26 02649 g008
Figure 9. Black circles (error bars): mass fractions wx(Si) and the associated uncertainties (k = 1) of six runs determined using the approach of this work. Red circles (error bars): results of the same input data (blends b2, b3, b4 only) using the approach of [1].
Figure 9. Black circles (error bars): mass fractions wx(Si) and the associated uncertainties (k = 1) of six runs determined using the approach of this work. Red circles (error bars): results of the same input data (blends b2, b3, b4 only) using the approach of [1].
Molecules 26 02649 g009
Table 1. Operation parameters of the mass spectrometric isotope ratio measurements applied for the three sample/matrix systems. (SP = spike; K = K factor).
Table 1. Operation parameters of the mass spectrometric isotope ratio measurements applied for the three sample/matrix systems. (SP = spike; K = K factor).
Sample/MatrixSilicon/TMAHSulfur/Biodiesel FuelTRF/Human Serum
LaboratoryPTBBAMPTB
InstrumentThermo
MC–ICP–MS Neptune
Thermo
MC-ICP-MS Neptune Plus
Agilent
8900 ICP-QQQ-MS
Sample
Introduction
PFA nebulizer 100 µL/min
PEEK/PFA cyclonic
+ Scott chamber
sapphire torch + injector
BN shield
Ni sampler + Ni X-skimmer
Aridus II desolvating system
PFA nebulizer 100 µL/min
Aridus PFA spray chamber
standard torch and injector
quartz shieldNi sampler + Ni H-skimmer
PFA MicroFlow nebulizer 700 µL, Scott chamber at 3 °C
torch with 1 mm injector
Pt shield
Pt sampler and skimmer
Gas Flow Rates (Ar)cooling: 16 L min−1
auxiliary: 0.8 L min−1
sample: 1.0 L min−1
cooling: 16 L min−1
auxiliary: 0.9 L min−1
sample: 0.85 L min−1
cooling: 15 L min−1
auxiliary: 0.9 L min−1
nebulizer gas: 0.8 L min1
reaction gas (H2): 6.1 mL min1
Machine
Parameters
high resolution (M/∆M = 8000)
RF power 1180 W
integration time 4.194 s
idle time 3 s
number of blocks 6
cycles/block 3
rotating amplifiers: yes
Faraday cups: L3(28Si), C(29Si), H3(30Si)
high resolution (M/∆M = 8000)
RF power 1200 W
integration time 4.194 s
idle time 3 s
number of blocks 1
cycles/block 40
rotating amplifiers: no
Faraday cups: L3(32S), C(33S), H3(34S)
MS/MS mode
RF power 1550 W
Sample depth 8.0 mm
x-lens configuration
integration time 0.1 s
m/z 53, 54, 56, 57, 58, 60
Sequence Settingsrinse time 120 s
take-up time 60 s
measured samples/sequence
b1, b2, b3, b4, b5 (4 times each)
rinse time 30 s
take-up time 80 s
measured samples/sequence
b1, b2, b3, b4, b5 (3 times each, separated by a block of 5 standards)
rinse time + take-up not applicable: HPLC separationmeasured samples/sequence
b1, b2, b3, b4, blank, SP, K, blank (4 times)
Separation
Settings
Agilent
Bioinert 1260 HPLC system
Column: MonoQ® GL 5/50 from GE Healthcare (Uppsala, Sweden)
Mobile phase A: 12.5 mmol/L Tris at pH = 6.4
Mobile phase B: 12.5 mmol/L Tris + 125 mmol/L NH4Ac at pH = 6.4
Flow: 0.5 mL min−1
Gradient:
0 min → 0% B,
20 min → 100% B,
27 min → 100% B
Column oven 30 °C
Injection volume: 10 µL
MWD 254 nm, 280 nm
Table 2. Determination of wx(Si) in TMAHaq. The relevant input and output data of the linear regression analysis are given for a single representative dataset (Sequence 1), with wz = 4.0069 µg/g.
Table 2. Determination of wx(Si) in TMAHaq. The relevant input and output data of the linear regression analysis are given for a single representative dataset (Sequence 1), with wz = 4.0069 µg/g.
xzy
biTMAHaqWASO04“Si30”Rb,i Rx,iRy,i
imx,imz,imy,iI(30Si)/I(28Si)xiyiI(30Si)/I(28Si)I(30Si)/I(28Si)
gggV/Vg/gg/gV/Vmol/mol
110.08630.000022.8557113.777320.00001.800.03353204.19578
29.78367.885822.75771.596550.8060301.510.03353204.19578
39.670010.525522.31981.188471.0885405.710.03353204.19578
411.319215.300022.48640.835311.3517503.860.03353204.19578
510.044022.406122.76510.614052.2308794.840.03353204.19578
a1a0wx
(g/g)/(g/g)(g/g)µg/g
356.1006211.4740.13
Table 3. Mass fractions wx (Si) of silicon in TMAHaq and the associated uncertainties (k = 1).
Table 3. Mass fractions wx (Si) of silicon in TMAHaq and the associated uncertainties (k = 1).
Runwx (Si)u(wx (Si))
µg/gµg/g
10.130.11
20.120.10
30.120.10
40.0440.012
50.0400.011
60.0280.007
average0.0810.073
Table 4. Determination of wx in BDF. The relevant input and output data of linear regression analysis are given for a single representative dataset (1st run), with wz = 16.145 µg/g. A procedural blank was subtracted from the result. The procedural blank was determined during an external measurement, as described in [23].
Table 4. Determination of wx in BDF. The relevant input and output data of linear regression analysis are given for a single representative dataset (1st run), with wz = 16.145 µg/g. A procedural blank was subtracted from the result. The procedural blank was determined during an external measurement, as described in [23].
xzy
biBDFNIST SRM 3154BAM S-34Rb,i Rx,iRy,i
imx,imz,imy,iI(32S)/I(34S)xiyiI(32S)/I(34S)I(32S)/I(34S)
gggV/Vg/gg/gV/Vmol/mol
10.237480.000000.096700.200300.000000.0038721.166430.00099
20.241490.101250.098430.367310.419270.0071821.166430.00099
30.238190.208280.108990.484530.874430.0107021.166430.00099
40.251260.303750.106310.649941.208910.0133921.166430.00099
50.243110.406790.099660.838871.673280.0169021.166430.00099
a1a0wxwx,corr
(g/g)/(g/g)g/gµg/gµg/g
0.0077970.0038948.0637.36
Table 5. Blank corrected mass fractions wx,corr(S) of sulfur in biodiesel fuel (uncertainties with k = 1).
Table 5. Blank corrected mass fractions wx,corr(S) of sulfur in biodiesel fuel (uncertainties with k = 1).
Runwx,corr(S)u(wx,corr(S))
µg/gµg/g
17.360.13
27.360.13
37.3580.079
average7.360.11
Table 6. Determination of wx(TRF) in human serum. The relevant input and output data of linear regression analysis are given for a single representative dataset (M21-3), with wz = 2.296 mg/g.
Table 6. Determination of wx(TRF) in human serum. The relevant input and output data of linear regression analysis are given for a single representative dataset (M21-3), with wz = 2.296 mg/g.
xzy
biSeronormTM Immuno-Protein Lyo L-1ERM®-DA470k/IFCCIn-House Prepared TRF SpikeRb,i Rx,iRy,i
imx,imz,imy,iR(54Fe/56Fe)xiyiR(54Fe/56Fe)R(54Fe/56Fe)
gggmol/molg/gg/gmol/molmol/mol
10.048480.000000.151043.011040.00000262.40.063703251.22
22.999580.00000263.4
32.991540.00000264.1
42.988870.00000264.4
50.049230.058340.151231.312821.18510614.6
61.321261.18510610.5
71.316561.18510612.8
81.313651.18510614.2
90.049100.074240.149091.152791.51220697.3
101.154261.51220696.3
111.157331.51220694.4
121.148921.51220699.8
130.049230.145610.149510.749592.957881109.1
140.747222.957881113.0
150.744402.957881117.6
160.744942.957881116.7
a1a0wx
(g/g)/(g/g)g/gmg/g
287.078266.042.128
Table 7. A representative result of mass fractions wx(Si) determined using the approach of this work (using all five blends) and the approach of [1] (using the three triple blends: b1, b2, b3; b2, b3, b4; and b3, b4, b5). Only the combination of b2, b3, b4 (the inner blends) yielded a reasonable numerical result.
Table 7. A representative result of mass fractions wx(Si) determined using the approach of this work (using all five blends) and the approach of [1] (using the three triple blends: b1, b2, b3; b2, b3, b4; and b3, b4, b5). Only the combination of b2, b3, b4 (the inner blends) yielded a reasonable numerical result.
This WorkApproach of [1]
blendsblends
b1, b2, b3, b4, b5b1, b2, b3
wx(Si)wx(Si)
µg/gµg/g
0.1292−0.5004
Table 8. A representative uncertainty budget of mass fraction wx(Si) determined using the approach of [1] (using the blends b2, b3, b4). The main contributions originate from Rb2, Rb1, and Rb3, with the largest absolute values of respective sensitivity coefficients.
Table 8. A representative uncertainty budget of mass fraction wx(Si) determined using the approach of [1] (using the blends b2, b3, b4). The main contributions originate from Rb2, Rb1, and Rb3, with the largest absolute values of respective sensitivity coefficients.
QuantityUnitBest Estimate
(Value)
Standard UncertaintySensitivity CoefficientIndex
Xi[Xi]xiu(xi)ci
wzµg/g4.006946.01 × 10−30.0310.0%
my1g22.496601.00 × 10−32.20.0%
mz2g10.312701.00 × 10−3110.0%
my3g22.268501.00 × 10−32.80.0%
Rb2V/V1.239333.62 × 10−39274.4%
Rz2V/V0.033527335 × 10−6
Rb3V/V0.884721.50 × 10−3−738.0%
Rb1V/V1.688465.46 × 10−3−3017.5%
my2g22.434701.00 × 10−3−5.00.0%
mz3g14.571701.00 × 10−3−4.20.0%
mz1g7.468401.00 × 10−3−6.40.0%
mx2g9.227301.00 × 10−30.330.0%
Rx2V/V0.033527335 × 10−6100.0%
mx3g9.790001.00 × 10−3−0.130.0%
mx1g9.464901.00 × 10−3−0.200.0%
wxg0.1250.388
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Brauckmann, C.; Pramann, A.; Rienitz, O.; Schulze, A.; Phukphatthanachai, P.; Vogl, J. Combining Isotope Dilution and Standard Addition—Elemental Analysis in Complex Samples. Molecules 2021, 26, 2649. https://doi.org/10.3390/molecules26092649

AMA Style

Brauckmann C, Pramann A, Rienitz O, Schulze A, Phukphatthanachai P, Vogl J. Combining Isotope Dilution and Standard Addition—Elemental Analysis in Complex Samples. Molecules. 2021; 26(9):2649. https://doi.org/10.3390/molecules26092649

Chicago/Turabian Style

Brauckmann, Christine, Axel Pramann, Olaf Rienitz, Alexander Schulze, Pranee Phukphatthanachai, and Jochen Vogl. 2021. "Combining Isotope Dilution and Standard Addition—Elemental Analysis in Complex Samples" Molecules 26, no. 9: 2649. https://doi.org/10.3390/molecules26092649

APA Style

Brauckmann, C., Pramann, A., Rienitz, O., Schulze, A., Phukphatthanachai, P., & Vogl, J. (2021). Combining Isotope Dilution and Standard Addition—Elemental Analysis in Complex Samples. Molecules, 26(9), 2649. https://doi.org/10.3390/molecules26092649

Article Metrics

Back to TopTop