1. Introduction
Platinum-group metals, comprising palladium, platinum, rhodium, ruthenium, osmium and iridium, are extensively used in several industrial processes since their peculiar physical-chemical properties [
1]. For example, these metals are often employed as efficient catalysts to synthesize drug compounds [
2,
3,
4] Platinum-group metals are also an essential part of vehicle catalyst systems to reduce the emission of gaseous pollutants, such as nitrogen oxides, carbon monoxide and hydrocarbons [
5], but consequently due to the increased numbers of cars the released concentration of Pd, Rh and Pt have rapidly increased in the environment [
5,
6,
7,
8,
9,
10].
Toxicological studies have indicated that Pd, in particular, being transported to biological systems, accumulates in the food chain resulting in potential health risks [
11,
12].
For these reasons, it is fundamental to develop analytical methods for determining Pd in environmental and biological matrices.
The most employed instrumental analytical techniques for Pd determination are atomic absorption spectrometry (AAS), inductively coupled plasma–mass spectrometry (ICP-MS), inductively coupled plasma-atomic emission spectrometry (ICP-AES) and neutron activation analysis (NAA) [
13]. Most of these methods present high sensitivity but are bulky and expensive; besides, they require time-consuming sample preparation, sometimes including separation and preconcentration steps before the instrumental analysis [
13,
14,
15].
Few works reported the application of electrochemical methods, which use mercury electrodes for palladium determination. Their application in stripping voltammetry ensures high sensitivity, low detection limits, and excellent reproducibility. However, these methods are ineffective for routine analyses due to the problems related to mercury's toxicity and disposal [
16,
17,
18].
Significant advances have been made in the development of sensors devoted to Pd(II) detection. Despite a couple of examples where electrochemical sensors were proposed [
19,
20], in the majority of cases, these are fluorimetric and colorimetric [
21,
22,
23,
24].
The development of solid-substrate-based colorimetric sensors for both qualitative and quantitative determination of metal ions in water samples is one of the most interesting and promising strategies for metal ions sensing since these devices have several advantages: they are sensitive, practical, and adapt to inexperienced users.
Among colorimetric sensors, Paper-based Analytical Devices (PADs) were widely applied for environmental analyses, food controls and clinical trials due to their easy preparation and quick response advantages [
25,
26,
27].
The capillary properties of paper provided by the cellulosic fiber network eliminate the necessity of pumping methods instead required in conventional microfluidic devices. Besides, the paper's ample availability, biodegradability, low price, and ease of surface functionalization make it an attractive substrate in sensors' development [
27,
28,
29].
Using colored complexes is a helpful strategy for colorimetric sensing on PADs and has been applied for detecting metal ions in different matrixes [
30,
31,
32,
33].
Differently from the univariate approach, the capability of chemometrics tools to extract a vast quantity of information can improve the performances of the PADs to discriminate and classify samples or for multiple-analytes detection [
27,
34,
35,
36,
37,
38]. Moreover, the proper application of chemometric tools to paper-based devices increases their accuracy, reliability, and robustness, thus strengthening their potential [
38].
The present work fits into this scenario; chemometric-assisted colorimetric PADs were developed for Pd(II) detection.
TazoC (disodium 2-[(1H-5-tetrazolyl)azo]-1,8-dihydroxynaphthalene-3,6 disulphonate, C
11H
6O
8N
6S
2Na
3·3 H
2O), an azo dye, was chosen as the chromophore receptor. This ligand makes the solutions in which it is dissolved and the solid phases impregnated by it orange-red colored. TazoC is not commercially available, but it was already studied as a complexing ligand for several metal ions [
23,
39,
40,
41,
42]; in particular, it forms a very stable complex with Pd(II), even in strongly acidic solutions. The complex Pd(II)/TazoC is purple-blue, very different from the free ligand's red-orange: this color change is the propriety used in the development of optical sensors [
23]. The structure of the Pd(II)/TazoC complex is shown in
Figure 1.
PADs' ability to detect Pd(II) was first studied at low pH (pH = 2) based on previous knowledge of the complex formation in aqueous solution [
40]. The dye does not complex other metal ions at such an acidic media, making TazoC-based PADs highly selective to Pd(II) detection. However, at higher pH values, other cations, such as Cu(II) and Ni(II), can be complexed by TazoC; in this case, the Pd(II) quantification is still possible by applying a multivariate approach. In particular, Partial Least Square regression (PLS) modeling [
43,
44,
45] was applied here to relate the Pd(II) concentration to the PADs' visible spectrum.
2. Materials and Methods
2.1. Reagents, materials and instruments
Cellulose filter paper Whatman grade 1 (180 μm thickness, 11 μm - particle retention, 150 sec/100 mL speed (Herzberg), and 0.25 psi wet burst, was obtained from Laboindustria S.p.A. (Arzergrande, Padova, Italy). Hydrochloric acid, sodium hydroxide, acetic acid, and sodium acetate were obtained from Merk Life Science S.r.l. (Milan, Italy); they were used to prepare the buffer solutions. Palladium standard, TraceCERT®, 1 g/L Pd in hydrochloric acid, Copper standard, TraceCERT®, 1 g/L Cu in nitric acid and Nickel standard, TraceCERT®, 1 g/L Ni in nitric acid (Merk Life Science S.r.l., Milan, Italy) were used to prepare diluted standard solutions.
TazoC was synthesized and purified according to previously described procedures [
39]; its empirical formula and molecular weight (536.34 g/mol) were confirmed by elemental analysis. A 1 mM ligand stock solution was prepared by dissolving the weighted solid in ultrapure water.
Tap water samples were obtained from the drinking water supply of Pavia, Italy. Samples were collected from the lab's sink (Department of Chemistry, University of Pavia, Italy) after flushing cold water for 15 minutes; they were subsequently acidified to pH 2 with HCl.
UV/vis spectra of the PADs were recorded with a Jasco V-750 spectrophotometer equipped with an FLH-740 film holder and a homemade clip designed to make the spectrum acquisition quick and easy (see
Figure 2).
The pH of the solutions was checked by a pH meter SevenMulti with an InLab Pro combined glass electrode (Mettler Toledo S.p.A.- Milan, Italy).
2.2. PADs preparation and analysis
The cellulose filter paper was cut into 2 cm-side squares using scissors. Each square, placed on a flat and clean surface, was drop-coated with 0.05 mL of 1 mM TazoC solution and left to air dry. The so-obtained PADs were immersed in 2.5 mL of aqueous metal-ion solution and kept under gentle stirring on a reciprocating shaker for 30 minutes, i.e., the time required to obtain a homogeneous and stable coloration of the PAD. The PADs were removed from the solution using plastic tweezers and left to air dry at room temperature for 2-3 min on a flat, clean surface. Subsequently, the PADs were inserted into the sample holder using plastic tweezers for the UV-vis analysis. The spectrum was registered in the wavelength range from 300 to 800 nm (bandwidth 0.2 nm, scan rate 200 nm/min) against a blank PAD wetted only with the buffer solution.
2.3. Chemometric data treatment: Partial Least Square regression (PLS)
The PLS bases have been widely described [
43,
44,
45,
46,
47,
48,
49], so they are only brief details here. PLS is a multivariate calibration algorithm aimed to extract most of the information contained in the UV-vis spectra (X variables), which is correlated to the Pd(II) concentration (Y variable) through the so-called latent variables (LVs). Each LV is a linear combination of the spectrum's absorbances. Unlike PCA, PLS finds directions of the latent variables explaining the maximum variance and maximizing the correlation with the response. The new coordinates of each object (the UV-vis spectrum) in the space of these latent variables are called scores, whereas the contribution of each spectral variable on the latent variables is called weight or loading.
The PLS model can be written as Y = X· b and allows the concentration (Y) prediction from the measured spectrum (X). Specifically, X is the n × m matrix of input data (absorbance values) obtained from the UV-vis spectrum: n is the number of samples, and m is the number of wavelengths of the whole spectrum. Y is the n × c concentration matrix, whereas c is equal to 1 in PLS1 models (which means when one analyte at the time is determined), and b is the n × c column vector of regression coefficients, solved when the model is calibrated. Knowing b, the PLS model can predict new Y concentration values from the UV-vis spectra of samples at unknown Pd(II) content.
The PLS models were built by selecting a suitable data set (training set) obtained with standard solutions at different Pd(II) content to cover the entire experimental domain homogeneously; three PLS models were developed at three different pH values (pH 2, pH 4, and pH 5.5). The training set for each model comprised the data of three replicates of 8-point calibrations (24 rows and 103 columns, i.e., 8 concentrations per 3 replicates and the absorbance values per 103 wavelengths). Each PLS model was tested by a cross-validated procedure on the training set and then on an external test set. Independent Pd(II) solutions were prepared for this purpose; the test set matrix comprised 9 rows (3 concentrations per 3 replicates) and 103 columns (absorbance values per 103 wavelengths).
Moreover, other datasets, comprising mixtures of Pd(II) and Cu(II) at pH 4 and Pd(II), Cu(II) and Ni(II) at pH 5.5 at different metal ions content, were prepared to evaluate the robustness of the method submitted to interferent cations.
A new training set was also prepared to build the PLS model for matrix-matched Pd(II) solutions (tap water). Then three independent tap water samples were fortified with different Pd(II) concentrations and used as the test set. From the predicted concentrations, the recovery percentages were calculated.
For all PLS models, only the centering pretreatment of data was applied. The chemometric data treatment was performed by the open-source R-based software CAT (Chemometric Agile Tool) [
50].
The tables with the metal ions concentrations for each dataset, the experimental and fitted Pd(II) concentrations and the graph of PLS models' performances were reported in the
Supplementary Materials.
3. Results
3.1. PADs preparation and analysis
The PADs proposed here were prepared simply with filter paper. In the recent literature, similar devices were realized using different procedures of paper pattering, such as screen-printing, wax printing or laser printing, inkjet etching, photolithography and plasma treatment [
51,
52,
53,
54,
55,
56,
57]. These techniques are expensive and sometimes require complicated procedures, making them unaffordable for limited resources laboratories. Moreover, PADs' wax-patterning needs rigorous flow control of the hydrophobic barrier into the porous matrix of the paper, causing poor homogeneity and, consequently, irreproducibility [
54]. An alternative method is a paper cutting into the desired shapes and dimensions, bypassing the need to create the liquid's containment barriers [
28]. Since the present study aimed to develop low-cost PADs using green materials and simple techniques, this last approach was adopted by cutting a sheet of filter paper into 2 cm-side squares.
As stated in the introduction, TazoC was chosen as the chromophore receptor since it forms a purple-blue-colored stable complex with Pd(II) also in very acidic media; this property has been exploited to develop the present devices. Indeed, until pH 2, only Pd(II) forms a strong complex with TazoC, as evident in
Figure 3, where the spectra of the free TazoC and the equimolar aqueous solutions of the ligand with Pd(II), Cu(II) and Ni(II) at pH 2 are reported.
The typical peak of the free ligand at about 460 nm (ε= 17333 cm
-1 M
-1) is in agreement with the data previously reported [
40], and the spectrum does not change if other cations than Pd(II) are added at pH 2. Besides, it can be observed that the Pd(II)/TazoC spectrum presents a broad and not well-defined peak, justifying the need for a multivariate calibration tool to quantify Pd(II).
Figure 4a,b show the spectra obtained in aqueous solution at higher pH values, i.e., pH 4 and 5.5, respectively. These pH values have been chosen because they correspond to the quantitative formation of the complexes with Cu(II) and Ni(II) without starting the precipitation of the respective hydroxides.
From the spectra, it is clear that Cu(II) at pH 4 and both Cu(II) and NI(II) at pH 5.5 can cause interference with the Pd(II) determination since their peaks are partially overlapped with that of the target analyte. All the more reason, also, in this case, the use of the multivariate regression PLS is necessary to perform the quantitative analysis, especially for data treatment of the PADs' spectra.
0.05 mL of 1 mM TazoC was used to load the filter paper since this volume is enough to cover the PAD's surface entirely without obtaining overflow of the ligand. As an alternative method, the immersion of the PAD in 2.5 mL of 1mM TazoC solution was tested, but low sorption kinetic and inhomogeneity of the PAD's color were verified.
It was then necessary to decide the modality to contact the sample solution with the PADs since it could affect the color homogeneity and intensity. Two approaches were tested: drop-coating with 0.5 mL of the solution or immersion of the PAD in 2.5 mL of aqueous metal-ion solution. The better uniformity of the color was achieved with the second strategy, which was thus adopted for all experiments.
The time required to obtain uniform PAD's color and avoid the leaching of the ligand was about 30 minutes with a gentle stirring on a reciprocating shaker.
Other kinds of filter papers were tested, hoping to reduce the immersion time, but due to the analogous porosity, their use does not confer advantages.
Regarding PADs' stability, it should be highlighted that these devices are disposable; moreover, the leaching of the chromophore receptor or its complexes with the metal ions can occur if they are stored in aqueous solutions for more than 1 h, so the analysis with the PADs must be immediately performed after their preparation.
Figure 5 shows the TazoC-PADs after contact with the metal ions solutions at different concentrations and pHs, and
Figure 6 shows the corresponding UV-vis spectra.
3.2. Chemometric-assisted Pd(II) determination by TazoC-PADs
As stated in the introduction and as seen in the picture of
Figure 5, a progressive color change of the PADs can be observed by the naked eye at increasing Pd(II) concentration. However, the simple correlation of the color (for instance, RGB indexes) with the metal ion concentration would not be sufficiently sensitive and accurate for quantitative analysis. For this reason, the UV-vis spectra were registered, and as can be observed in
Figure 6 (a-c), it cannot be possible to identify a well-defined peak related to the Pd(II)/TazoC-PAD complex; consequently, a classical univariate calibration shall not be applicable.
The Partial Least Squares regression (PLS) was selected; indeed, this multivariate technique makes it possible to correlate the entire spectrum to the analyte concentration.
Without any other data treatment, conversely applied for very disturbed spectrum (i.e., baseline correction, lines' smoothing), the PLS tool makes it possible to empathize differences among entire spectra rather than differences within the spectra (noise).
In the present research, tailored PLS models were built for the different studied media.
The first PLS model is the one related to the simplest system, i.e., Pd(II)/TazoC-PADs in aqueous solutions at pH 2 (HCl 0.01M).
This model was built with a training set comprising the data of three replicates of an 8-point calibration with the Pd(II) concentration ranging from 2.6 to 50 μM. The test set to validate and prove the model's robustness comprised three replicates of different TazoC-PADs immersed in three Pd(II) solutions (7.5, 25.2, 42.5 μM).
A 5-latent variables model was developed: it ensures 98.61% of explained variance in Cross-Validation (% E.V. in CV), a global Root Mean Square Error in CV (RMSECV) of 2.05 μM and, as regards the test set samples, 1.86 μM of Root Mean Square Error in Prediction (RMSEP). The model performance graph is reported in the Supplementary Material (
Figure S1).
Figure 7a shows the Experimental vs. Fitted values plot for the training set (burgundy-colored points) and test set (light blue-colored points); samples are similarly distributed alongside the y=x straight line, and no significant difference in the fitting error appears. The residuals plot is reported in
Figure 7b; a random distribution around 0 was observed, and a point cloud without a particular structure (the points are placed between the values -3 and 3) justifies the model's linearity and homoscedasticity assumptions.
Operational values of LOD and LOQ were also obtained by projecting in the PLS model the spectra of 10 blank samples (i.e., TazoC-PADs contacted with solutions at pH 2, without Pd(II), see the yellow-colored points in
Figure 7a); the LOD and the LOQ were then computed as 3.3 times and 10 times, respectively, the standard deviation of the predicted blank concentrations. The values are reported in
Table 1. This procedure is not precisely rigorous in the case of PLS, but it seems a reasonable estimate here since the predicted values are around zero at low concentrations with a very low standard deviation [
23]. The authors are well aware that the correct procedure for LOD and LOQ estimates in PLS analysis is still an open question.
Figure 8 shows similar plots obtained for the systems Pd(II)/TazoC-PADs at pH 4 and Pd(II)/TazoC-PADs at pH 5.5.
The figures of merit and the LOD and LOQ values for all three PLS models were summarized in
Table 1.
The three models' robustness and predictive ability were demonstrated; indeed, the training set and test set samples are analogously distributed alongside the y=x straight line, and no significant differences in the fitting errors arise. As expected, the operational LOD and LOQ values are similar independently of the solution pH since the quantitative formation of the Pd(II)/TazoC-PADs complex also occurs in very acidic media.
Moreover, the lowest quantifiable concentration values (LOQ) achieved were similar to or, in most cases, lower than those obtained with previously proposed test strips based on different colorimetric or fluorescent probes (see
Table 2).
3.3. TazoC-PADs applications: interference tests and spiked tap waters analysis
As for any analytical method subject to interference or in the presence of complex real matrices, calibration with external standards is often ineffective for analyte quantification [
63,
64,
65]. As a demonstrative example,
Figure 9 shows the PLS plot of the model Pd(II)/TazoC-PADs at pH 4, where mixtures of Pd(II) and Cu(II) at different concentrations were used as test set samples. The model failed: the Pd(II) concentrations were incorrectly predicted.
The usual strategy to overcome interferences and complex matrix problems is the application of the standard additions method or the matrix-matched calibrations [
63,
64]. The second approach was adopted here, developing tailored PLS models both in the case of interferents and for spiked tap water samples analysis.
The following
Figure 10 and
Table 3 report the PLS plots and the figures of merits for the models named Pd(II)+Cu(II)/TazoC-PADs pH 4 and Pd(II)+Cu(II)+Ni(II)/TazoC-PADs pH 5.5. In these cases, mixtures of variable different and independent concentrations of the metal cations were considered for both training and test set samples (see
Supplementary Materials).
Both PLS models proved to be adequate in predicting the Pd(II) concentrations even in the presence of interfering cations since the pretty good agreement between experimental and fitted values. The residual plots' analysis demonstrates, also in these cases, the linearity and homoscedasticity of the models as the random distribution of the points.
The applicability of the TazoC-PADs to tap water samples was verified by adopting the same approach of matrix-matched calibration used for the interference tests.
As recommended from the preservation of water samples for metals' analysis rules [
66], acidification at pH < 2 must be done to avoid the precipitation of hydroxides. Accordingly, the PLS model was built using tap water samples, acidified at pH 2 and fortified with Pd(II) (model name: Pd(II)/TazoC-PADs TW). Three tap water samples were acidified and spiked with different Pd(II) concentrations (7.5, 25.2 and 44.6 µM) and used as the external dataset (3 replicates for each sample).
The PLS model's figures of merit are reported in
Table 4; in
Figure 11a, the Experimental vs. predicted values plot for the training set (light blue-colored points) and test set (yellow-colored points) is shown. In
Figure 11b, the residual plot is reported. The recovery percentage was computed from the test set samples' predicted concentration values and reported in
Table 5, together with the Pd(II) concentration values obtained by ICP-OES analysis for comparison.
The percentage of recoveries (Rc%) for each Pd(II) concentration is between 80% and 110%, i.e., the acceptable recovery range [
67,
68], demonstrating the suitability of the proposed method for Pd(II) determination in drinking waters. Moreover, the concentration values obtained are not significantly different from those achieved by the classical ICP-OES technique.
4. Conclusions
A cheap, selective colorimetric Paper-based Analytical Device (PADs) for Pd(II) sensing and quantification is proposed. The chromogenic receptor immobilized on filter paper is a non-commercial azo dye named TazoC (disodium 2-[(1H-5-tetrazolyl)azo]-1,8-dihydroxynaphthalene-3,6 disulphonate); it is selected because it forms a strong and stable purple-blue-colored complex with Pd(II) also in very acidic media (below pH 2).
A proper application of the chemometric tool, PLS, permitted to build spectrum/Pd(II) concentration correlation models, taking advantage of using the whole spectrum as the signal accounting for changes in shape and height of the spectrum peaks. Different tailored PLS models were developed and validated, highlighting the need to perform calibrations in the media of interest.
The lowest quantifiable concentration values achieved, of about 2.5 μM, were similar to or, in most cases, lower than those obtained with previously proposed test strips based on different colorimetric or fluorescent probes. Moreover, according to the WHO threshold limit for Pd content in drug chemicals (from 47.0 μM to 94.0 μM, i.e., from 5 mg/L to 10 mg/L) [
69], and the United Nations Food and Agriculture Organization (FAO) recommended maximum level for irrigation waters of 47.0 μM (5 mg/L) [
70], the lowest quantifiable Pd(II) concentration obtainable with the proposed PADs meet the requirements of WHO and FAO for palladium(II) detection in that matrixes.