The use of a new MultiNeb® nebulizer for the on-line isotope dilution inductively coupled plasma mass spectrometry (OID-ICP-MS)
1. Introduction
The development of analytical methodologies for the accurate and reliable determination of trace elements in samples of different kinds (e.g., sediment, biota organisms, or seawater) is in continual demand. Conventionally, testing laboratories employ analytical classical calibration techniques for this purpose which require the establishment of timeconsuming calibration graphs and further recalibration procedures. Whenever many samples need to be analyzed daily, the implementation of methods that do not require the establishment of a calibration graph may be advantageous as it involves significant savings, in time and cost.
The isotope dilution (ID) methodology is based on modifying the isotopic composition of the element to be determined in the sample by adding a known amount of an isotopically enriched standard of the same element (i.e., spike solution) [1]. The total amount of the element in the sample can be obtained from the isotope ratio measurement in the resulting blend. In addition, isotope dilution inductively coupled plasma mass spectrometry (ID-ICP-MS) has been identified to have the potential to be a primary method of measurement (PMM).
Hence, the ID-ICP-MS method presents several advantages over other analytical methodologies conventionally employed in laboratories for trace elements determination, including:
(i) once complete isotope equilibration between the sample and the spike is achieved, there is no need to know the pre-concentration or dilution factor of the sample or to consider any non-quantitative separation or evaporation process;
(ii) any variation of the instrumental sensitivity such as signal drift has not influence on the final results; and
(iii) the results are obtained with superior accuracy, and with a small combined uncertainty. In short, its basic advantage is that ID-ICP-MS measurements must be accepted without the use of standards (i.e., the calibration graph is avoided).
Traditionally, ID-ICP-MS is practiced only in highly sophisticated laboratories for reference measurements and for any cases where the high accuracy criterion are of great importance to the analytical results. Batch ID-ICP-MS is also regarded as being expensive and time-consuming since each sample is spiked individually. Sample spiking is tedious, and moreover, a main source of error as it occurs prior to isotope equilibration. However, it can be simplified by the on-line ID-ICP-MS (OID-ICP-MS), which consists of the on-line addition of the spike solution to the sample.
The new MultiNeb® allows for on-line sample spiking by simultaneously introducing the spike (i.e., isotopically enriched standard) and sample (i.e., natural isotopic composition) solutions and ensuring effective mixing between them. This is achieved by creating high-pressure turbulent conditions at the inner part of the nebulizer tip.
A new OID-ICP-MS methodology based on the MultiNeb® is presented. The objective is to enhance the productivity of routine laboratories in analyzing samples with complex matrices by implementing the MultiNeb® in on-line ID-ICP-MS analysis. This results in high-quality results, as well as a reduction in in human intervention required for individual spiking,
shorter total analysis time, and cost savings.

2. Experimental
Instrumentation
Acid digestion of the sample was carried out in a closed vessel device using a Milestone Start D microwave acid digestion system (Sorisole, Italy), equipped with polytetrafluoroethylene vessels. The samples were digested at conditions recommended by the manufacturer (Application note HPR-FO-17 for dried fish).
All isotope ratio measurements were performed with an inductively coupled plasma mass spectrometer (model 7700x, Agilent Technologies, Santa Clara, CA, United States), working in helium collision mode with the 3rd generation Octopole Reaction System (ORS3) for better polyatomic interferences removal. A standard quartz torch (i.e., 2.5 mm internal diameter) and standard nickel cones (i.e., sampler and skimmer) were used, while no instrument’s autosampler was used throughout the work. The ICPMS operating conditions, which were auto-tuned using the 7700 MassHunter software, are listed in Table 1. Total analysis time per sample, including wash-in and wash-out, was 3 minutes.
The selection of the isotopes to be measured in the ICP-MS was done with respect to the availability in the spike material, abundance of the isotopes, and possibility of spectral interferences. The signal intensities used for isotope ratio measurements were corrected for instrumental background, dead time, and possible spectral interferences.
Reagents and materials
High-quality Milli-Q water (18.2MΩ·cm resistivity) obtained from PureLab Flex 3 system (Lane End, United Kingdom) was used throughout this work. Ultra-pure HNO (Merck, Darmstadt, Germany) and HO (Merck) were used for sample digestion.
Quality control (QC) solutions, QC1 and QC2, were prepared by appropriate gravimetric dilution of single-element 1000 mg L-1 standard solutions with natural isotopic composition from High Purity Standards (Charleston, SC, USA) with a solution containing 2% (w w-1) nitric acid. The QC1 concentration is 20 ng g-1 for Cu, Ni, Pb and Zn and the QC2 concentration is 10 ng g-1 for Cu, Ni, Pb and Zn. These solutions were used for mass discrimination correction during ICP-MS measurements.
The natural isotopic abundances of the elements of interest, except for lead, were taken from IUPAC tables [2] and are shown in Table 2. Lead natural isotopic composition instead, was analyzed prior to isotope ratios measurements in the ICP-MS [3] and is also shown in Table 2. Finally, a fish homogenate Certified Reference Material IAEA-476, that was obtained from the International Atomic Energy Agency (IAEA, Principality of Monaco).

Table 1. ICP-MS operating conditions for the on-line isotope ratio determinations.

a Expanded uncertainty for each value was calculated as U = 2·u, where k=2 is the coverage factor for a 95% confidence interval and u is the combined standard uncertainty.
b From reference 2.
c Determined experimentally by ICP-MS measurements.
Sample pretreatment
A portion of ~0.2 g of the biota sample was weighed directly into a polytetrafluoroethylene vessel and 7 mL of HNO3 65% (w w-1) and 1 mL of H2O2 30% (w w-1) were subsequently added. The heating program was applied in two steps: (1) 15 min to reach 200 °C and (2) 15 min at 200 °C, and an additional 15-min cooling step. A maximum 1.5 kW of microwave power was applied. After completing the digestion and cooling down steps, the final digests were diluted to 50 g of final solution with Milli-Q water, filtered through syringe with 0.45 μm pore size polyvinylidene fluoride filters (Millipore, Madrid, Spain) and then stored in polyethylene tubes at 4ºC until the analysis.
In order to avoid memory effects from previous experiments, the set of polytetrafluoroethylene vessels employed for sample digestion were first cleansed using a microwave procedure, consisting of the addition of 8 g of concentrated HNO3 to each vessel, followed by a microwave treatment at 350 W during 10 min to reach 100 ºC and 10 min at 100 ºC (i.e., 20 min in total). This procedure was performed twice before sample digestion.
Four procedural blanks were also prepared and subjected to the entire analytical methodology together with the fish homogenate replicates in order to evaluate the contribution of the blank contamination. Six replicates were prepared for the analysis.
On-line isotope dilution
The nebulizer used was a multi nebulizer MultiNeb® (Ingeniatrics, Seville, Spain) (Figure 1) which incorporates two independent liquid inlets into a single nebulization body with a common nebulization gas inlet and a single outlet orifice. The liquid streams are mixed at the tip inside the nebulizer in the aerosol phase at high pressure and the resulting aerosol from the mixture of the liquids exits through the orifice. As a result, there is an increase in the mixing efficiency and in chemical reaction speed.

Figure 1. Side view (left) and front view (right) of the MultiNeb® multiple nebulizer prototype.
The liquid sample introduction system was composed by a MultiNeb® nebulizer coupled to a quartz double pass Peltier-cooled (2ºC) spray chamber without any additional modification required, as the MultiNeb® is built on the right dimensions to allow the easy connection to any commercial spray chamber conventionally used in ICP-based techniques. Two Tygon® peristaltic tubes for aqueous solutions (R-3607, white-white code color, i.d. 1.02 mm, Ismatec Cole-Parmer GmbH, Wertheim, Germany) were used for the sample/QC solutions and spike solutions.
A scheme of the sample introduction system for OIDICP-MS is shown in Figure 2. The on-line sample spiking was achieved by using both liquid inlets of the MultiNeb®: one inlet (Channel 1 in Figure 2, flow rate: 0.3 mL min-1) for the continuous multispike solution addition and the other one (Channel 2 in Figure 2, flow rate: 0.3 mL min- 1) for the sequential introduction of the procedural blanks, the liquid samples and the QC solutions. Thus, the two solution streams were mixed inside the tip of the nebulizer, leading to a primary aerosol of the resulting blend.

Figure 2. Scheme of the on-line isotope dilution system using the MultiNeb®.
Analyte concentrations in the sample digests, Cs, were calculated using the following simplified equation, obtained from the isotope dilution theory [1]:

where Cs is the concentration of the analyte in the sample digest solution (s); Csp, the concentration of the analyte in the spike solution (sp); Fst, molar flow of the spike solution (Channel 1); Fs , molar flow of the sample solution (Channel 2); Rm, the measured isotope ratio of the mixed sample and spike solution (s + sp); Ab sp, abundance of isotope b (Table 2) of the spike solution (sp); Aa sp, abundance of isotope a (Table 2) of the spike solution (sp); Aa s , abundance of isotope a (Table 2) of the sample solution (s); Ab s , abundance of isotope b (Table 2) of the sample solution (s).
From the equation above, it seems clear that the only parameters that needed to be experimentally determined for the calculation of Cs were Rm and Rst. Both isotope ratios were determined by OID-ICP-MS measurements and then incorporated into the equation.
Quality control solutions were assumed to contain natural isotopic composition of Cu, Ni and Zn and their associated uncertainties were taken from IUPAC [2]. Lead isotopic composition instead, was determined experimentally by ICP-MS measurements within the same analytical run (Table 2).
Analysis Procedure
A scheme of the analysis procedure is shown in Figure 3. The standard-sample bracketing method [4] was employed for calibration and an effective mass discrimination correction. Quality control solution (QC1: 20 ng g-1 for trace elements) with known isotopic composition was analyzed at the beginning of the analysis. Then, procedural blank solutions were analyzed. The sample was later analyzed.
Finally, QC2 solution with natural isotopic composition but lower concentration values (QC2: 10 ng g-1 for trace elements) was also analyzed every three sample replicates in order to correct for mass bias and to verify that the obtained recovery values were close to 100%. Each sample replicate was measured three times and average final concentration as well as standard deviation was calculated for each element.

Figure 3. Analytical procedure for the on-line isotope dilution analysis using the MultiNeb® nebulizer with two independent liquid inlets (Channel 1 and 2).
3. Results and discussion
CRM analysis and uncertainty estimation
In traditional isotope dilution analysis, the optimal amount of spike to be added to a sample is determined by calculating the ideal ratio using the «error magnification factor“. This factor depends on the isotopic abundances of the enriched element spike, as well as the natural isotopic abundances of that element in the samples. However, for on-line isotope dilution analysis, the spike amount is typically constant, but it is preferable to “overspike” the samples to yield better counting statistics and therefore less uncertainty in the isotope ratio measurements [5].
The analytical performance of the proposed OID-ICPMS method was evaluated. As it can be seen in Table 3, the concentrations found in the fish homogenate reference material were in good agreement with the certified and information values for all the analytes. Accuracy results obtained with the proposed method was verified by evaluating recovery values of each analyte. Recovery values were found within the range of 94-109 %.

a Mean value expanded uncertainty with a coverage factor k=2.
b n is the number of laboratories participating in the characterization of the IAEA-476 candidate CRM.
c Information value
The relative expanded uncertainties (k=2) for the obtained values of the investigated elements were as follow: on the Cu (U´ =3.5 %), Ni (U´=2.6 %), Pb (U´=2.6 %) and Zn (U´=2.4 %). These low uncertainty values offer a great advantage over other conventional analytical methodologies commonly employed for the same purpose, as the analysis of trace metals in environmental samples often generates data with large uncertainty values. It was demonstrated that OID-ICP-MS was less susceptible to matrix effect that other calibration methodologies (e.g., external calibration).
There are several practical benefits of implementing the OID-ICP-MS methodology using the MultiNeb®. Firstly, the individual spiking step is avoided, which is often a main error source of the expanded uncertainty estimation as it is prior to isotope equilibrium, thereby also reducing sample handling and analysis time per sample. This fact has a direct beneficial impact on the costs of the analysis. The only sample pretreatment required in this method was sample digestion, that could even be omitted in the case of liquid samples (e.g., wastewater). Moreover, it should be noted that the analysis of the fish sample using OID-ICP-MS just required the preparation of two QC solutions for all the analytes and that they could be used for the analysis of more samples. The method has been implemented in the IAEA-476 CRM simply for the convenience, but several samples with high matrix, which are commonly employed in toxicological and ecological risk monitoring (e.g., wastewater, sediment, seawater) could have been used instead.
4. Conclusions
It has been demonstrated that the new MultiNeb® is suitable for on-line isotope dilution analysis of trace elements in biota sample by ICP-MS.
Its implementation for the online isotope dilution provides significant advantages: (i) it simultaneously corrects matrix effects (i.e., matrix-based polyatomic interferences) and any variation of the instrumental sensitivity (i.e., signal drift); (ii) it significantly reduces the total time of the analysis as the timeconsuming individual spiking step of conventional IDA is avoided, no external calibration graph is needed (making unnecessary any recalibration and re-analysis procedures), and only one standard solution per sample is necessary (i.e., reference standard solution); (iii) it provides ease of operation, as the system is simple (i.e., the use of valves and specific components is not required), robust, easy to handle, the multi nebulizer perfectly couples to commercial spray chambers (without any additional modification), and multi-element isotope enriched spike solutions are commercially available; (iv) it provides high-quality reliable information with lower uncertainty values (i.e., within 1-5 % relative combined uncertainty for the evaluated elements) compared to conventional quantification methods and excellent agreement with IAEA-476 CRM certified values; and (v) it could be applied to the analysis of many analytes and with different samples.
All these features result in a promising solution for the trace element analysis of real-world samples which simplifies operation and significantly increases sample throughput and, thus, enhances productivity, with associated economic benefits. It yields savings in total cost of the analysis per sample, despite the relatively high cost of the enriched isotope standard.
5. References
[1] J. I. García Alonso and P. Rodríguez-González, Isotope Dilution Mass Spectrometry, Isotope Dilution Mass Spectrometry, Royal Society of Chemistry, Cambridge, 2013.
[2] J. Meija, T. B Coplen, M. Berglund, W. Brand, P. De Bièvre, M. Gröning, N. Holden, J. Irrgeher, R. D Loss, T. Walczyk and T. Prohaska, Isotopic compositions of the elements 2013 (IUPAC Technical Report), Pure Appl. Chem., 2016, 88, 293– 306.
[3] J. W. Valley and D. R. Cole, Stable isotope geochemistry, Mineralogical Society of America: Washington DC, 2019.
[4] E. Bolea-Fernández, L. Balcaen, M. Resano and F. Vanhaecke, Tandem ICP-Mass Spectrometry for Sr isotopic analysis without prior Rb/Sr separation, J. Anal. At. Spectrom., 2016, 31, 303–310.
[5] G. Centineo, J. A. Rodríguez-Castrillón and E. Muñoz Agudo, On-line isotope dilution analysis with the 7700 Series ICP-MS : Analysis of trace elements in high matrix samples, Appl. Note – Agil. Technol., 2011, 1–7.