Next Article in Journal
Incorporation of an Intermediate Polyelectrolyte Layer for Improved Interfacial Polymerization on PAI Hollow Fiber Membranes
Previous Article in Journal
Electrochemical Characterization of Charged Membranes from Different Materials and Structures via Membrane Potential Analysis
Previous Article in Special Issue
On the Use of Polymer Inclusion Membranes for the Selective Separation of Pb(II), Cd(II), and Zn(II) from Seawater
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

MID-FTIR-PLS Chemometric Analysis of Cr(VI) from Aqueous Solutions Using a Polymer Inclusion Membrane-Based Sensor

by
Armando Martínez de la Peña
,
Eduardo Rodríguez de San Miguel
* and
Josefina de Gyves
Departamento de Química Analítica, Facultad de Química, Universidad Nacional Autónoma de México (UNAM), Ciudad Universitaria, Mexico City 04510, Mexico
*
Author to whom correspondence should be addressed.
Membranes 2023, 13(8), 740; https://doi.org/10.3390/membranes13080740
Submission received: 30 July 2023 / Revised: 11 August 2023 / Accepted: 16 August 2023 / Published: 18 August 2023
(This article belongs to the Special Issue New Trends in Polymer Inclusion Membranes 2.0)

Abstract

:
A partial least squares (PLS) quantitative chemometric method based on the analysis of the mid-Fourier transform infrared spectroscopy (MID-FTIR) spectrum of polymer inclusion membranes (PIMs) used for the extraction of Cr(VI) from aqueous media is developed. The system previously optimized considering the variables membrane composition, extraction time, and pH, is characterized in terms of its adsorption isotherm, distribution coefficient, extraction percent, and enrichment factor. A Langmuir-type adsorption behavior with KL = 2199 cm3/mmol, qmax = 0.188 mmol/g, and 0 < RL < 1 indicates that metal adsorption is favorable. The characterization of the extraction reaction is performed as well, showing a 1:1 Cr(VI):Aliquat 336 ratio, in agreement with solvent extraction data. The principal component analysis (PCA) of the PIMs reveals a complex pattern, which is satisfactorily simplified and related to Cr(VI) concentrations through the use of a variable selection method (iPLS) in which the bands in the ranges 3451–3500 cm−1 and 3751–3800 cm−1 are chosen. The final PLS model, including the 100 wavelengths selected by iPLS and 10 latent variables, shows excellent parameter values with root mean square error of calibration (RMSEC) of 3.73115, root mean square error of cross-validation (RMSECV) of 6.82685, bias of −1.91847 × 10−13, cross-validation (CV) bias of 0.185947, R2 Cal of 0.98145, R2 CV of 0.940902, recovery% of 104.02 ± 4.12 (α = 0.05), sensitivity% of 0.001547 ppb, analytical sensitivity (γ) of 3.8 ppb, γ−1: 0.6 ppb−1, selectivity of 0.0155, linear range of 5.8–100 ppb, limit of detection (LD) of 1.9 ppb, and limit of quantitation (LQ) of 5.8 ppb. The developed PIM sensor is easy to implement as it requires few manipulations and a reduced number of chemical compounds in comparison to other similar reported systems.

1. Introduction

Due to the extensive and rapid development, environmental pollution has increased to alarming levels in nearby industrialized areas, which are generally contaminated with heavy metals. Among many others, hexavalent chromium, Cr(VI), is one of the most serious concerns, as there is credible evidence that through the water medium, it is a major contributor to the global burden of cancer in humans [1]. Chromium compounds are mainly used in industrial activities, such as corrosion control, oxidation process, leather tanning, electroplating, metallurgy, cement, textile dyeing, papermaking, inks, paints and pigments, and photographic industry [2]. Due to this wide use, there are locations where chromium compounds have been released to the environment via leakage, poor storage, or improper industrial disposal practices, so chromium pollution of waters and groundwaters represents a serious environmental concern. The USEPA regulates total chromium in drinking water and has set a Maximum Contaminant Level (MCL) of 0.1 mg/dm3 [3]. The World Health Organization (WHO) guideline is 0.05 mg/dm3 for total chromium [4]. Although this reference value has been questioned, as a practical measure, it has been retained as a provisional guideline value until additional information becomes available and chromium can be re-evaluated [4]. Chromium discharge limits in water are regulated on a national scale and often vary depending on the different types of industry or receiving water bodies (marine water, lake, river, or sewer system) [5]. The maximum discharge limit to the aquatic environment in the EU is 1 and 5 mg/dm3 for Cr(VI) and Crtotal, respectively [6].
There are currently several technologies for the remediation [5] and detection [7] of Cr(VI) from aqueous samples. As this species is usually present at low concentration levels, it is measured by atomic absorption and emission methods, providing high precision and sensitivity measurements. The selection of these techniques will depend largely on the sensitivity, reproducibility, detection and quantification limits, and simplicity of the method. However, environmental scientists are increasingly in need of measurement closer to the sample, in situations outside the laboratory, which restricts the use of such analytical techniques. Alternatively, UV/VIS, FTIR, Raman, and NIR portable spectrometers [8,9] are available to meet the growing need driven by the general trend in analytical instrumentation toward smaller size, improved reliability, and greater ease of operation. In such a context, the use of ion-selective optodes (ISOs), ionophore-based optical chemical sensors, represents an excellent alternative for the portable determination of ions [10]. In this regard, several systems for Cr(VI) monitoring using optical sensors have been developed in recent years [11,12,13,14].
On the other hand, besides being used for the removal and transport of metals [15,16], polymer inclusion membranes (PIMs) have also been employed for the detection of many chemical species [17]. These membranes are composed of a polymeric matrix, which serves as support of an encapsulated extractant (carrier) that is responsible for binding with the target analyte at the source solution/membrane interface and transporting it across the membrane, which can contain a plasticizer to improve mobility. In the case of Cr(IV), Aliquat 336 (methyltrioctylammonium chloride) has proved to be an effective and selective extractant in solvent extraction (SX) [18,19], supported liquid membranes (SLMs) [20], solvent-impregnated resins (SIR) [21,22], and PIM [20,23,24,25,26,27,28,29] systems.
In previous work, the capabilities of a PIM sensor to perform cadmium(II) determination in aqueous solutions by in-situ visible (VIS) and mid-Fourier transform infrared spectroscopy (MID-FTIR) analyses of the polymeric films using a partial least squares (PLS) chemometric approach were demonstrated [30]. One major advantage of the developed MID-FTIR-PLS PIM-based method was that it did not require the presence in the membrane of a chemical reagent with special properties, either a chromophore species that can complex the metal ion, i.e., acting as an ionophore [14,31], or a mixture of an ionophore and a chromophore in the same PIM [32], or a fluorescent reagent [33]. However, to extend the potentiality of the methodology to low analyte concentration ranges, a careful selection of the dielectric nature of the medium and the dipole moment of the bond associated with IR vibrations of the extracted complex was suggested [30], so that an increase in the band intensities could be achieved, i.e., an augment in the magnitude of the analytical signal, improving the detection capabilities of these sensors.
Taking advantage of the high dielectric constant values observed in PIMs with Aliquat 336 plasticized with 2-nitrophenyl octyl ether (NPOE) [34] and the excellent transport properties of Aliquat 336 for Cr(VI) in PIMs [20,23,24,25,26,27,28,29], in this work, such membranes were employed for MID-FTIR-PLS analysis of Cr(VI) from aqueous media. Once the system had been characterized and optimized, the calibration performance and the figures of merit (FOM) of the method were determined to show the potentiality of the proposed methodology.

2. Materials and Methods

2.1. Reagents and Apparatus

PIMs were prepared using cellulose triacetate (CTA, Honeywell Fluka, Charlotte, NC, USA) as polymer support, 2-nitrophenyl octyl ether (NPOE, ≥99.0% Honeywell Fluka) as a plasticizer, Aliquat 336 (≥97%, methyltrioctylammonium chloride, Sigma-Aldrich, Chem. Co., St. Louis, MO, USA) as an extracting agent, and dichloromethane (Merck, Kenilworth, NJ, USA) as a casting solvent. Working Cr(VI) solutions were prepared by dissolving the corresponding amounts of (NH4)2Cr2O7 (≥99.5%, Sigma-Aldrich) in deionized water (18 MΩ cm, MilliQ, Merck Millipore, Burlington, MA, USA). A 1000 mg/dm3 Sigma-Aldrich AAS standard solution (≥99.5%, Sigma-Aldrich, 1 g/dm3 Cr in 2% nitric acid, prepared with high purity (NH4)2Cr2O7, HNO3 and water) was diluted using deionized water for the preparation of the standards for flame atomic absorption spectrometry (FAAS) determinations. Tris (hydroxymethyl)aminomethane (TRIS, 99.8% Sigma-Aldrich), sodium acetate (99% Sigma-Aldrich), acetic acid (99.7% Sigma-Aldrich), hydrochloric acid (37% Sigma-Aldrich), and 2-(N-morpholino) ethanesulfonic acid (MES, 99.5% Sigma-Aldrich) were employed to assess the effect of pH on chromium extraction.
A model 75Wrist ActionTM shaker (Burrell Scientific Inc., Pittsburgh, PA, USA) with 10 speeds and a Metrohm 620 pH-meter (Herisau, Switzerland) were used. A 3100-flame atomic absorption spectrometer (Perkin Elmer Waltham, MA, USA) and a Spectrum GX IR spectrometer (Perkin Elmer) were employed to measure metal ion contents and obtain the IR spectra, respectively. A Fowler IP54 micrometer (Fowler High Precision, Newton, MA, USA) was used for measuring PIM thicknesses.

2.2. Preparation of PIMs

Membranes were prepared dissolving 30 mg of CTA, 75 mg of NPOE, and 30 mg of Aliquat 336 in 10 cm3 of CH2Cl2. This mixture was stirred for 1 h in a 50 cm3 beaker and then the solvent was allowed to evaporate for 24 h. After this time, the membranes were detached from the beakers and their diameter and thickness were measured with a Fowler IP54 micrometer. An average weight of (135.86 ± 2.05) mg, an average thickness of (95 ± 15) µm, and an average diameter of (3.00 ± 0.02) cm were determined (n = 54). The PIM thickness is almost double that previously reported in a similar sensor [30] to enhance the detectability of the analyte.

2.3. Liquid-Solid Extraction

The obtained membranes were placed in a 50 cm3 polypropylene falcon tube together with 30 cm3 of Cr(VI) solutions at different concentrations. Stirring was performed for 20 min, taking a 1 cm3 aliquot every 2 min. At the end, the membranes were removed from the aqueous phase and the aliquots were brought to a volume of 2 cm3 to be able to measure their concentration by FAAS. Experiments were performed on a duplicate basis with an average RSD of 5%.

2.4. Metal Quantification

The absorbances corresponding to the samples and standards were measured by FAAS according to the conditions established by the manufacturer (354.7 nm wavelength (λ), 7 nm slit, with an air-acetylene flame (99.6%, Praxair, Danbury, CT, USA and an oil-free GS-003—Air compressor, PG Instruments, Loughborough, Leicestershire, UK), from 0 to 7 mg/dm3, using a hollow cathode lamp, and sensitivity check of 5 mg/dm3). The calibration curves were performed at each of the experimental conditions employed, mismatching the standards and experimental samples.

2.5. Measuring Infrared Spectra

The IR spectra were obtained by direct analysis of the PIMs in the transmission mode for the quantitative analysis of chromium in the mid-infrared region. The PIM was mounted on the transmission accessory of the equipment after sandwiching the membrane between two Petri dishes to avoid wrinkles and movement. The spectrum of each sample was acquired in the 400–5000 cm−1 region in triplicate with 30 scans per measurement to reduce the within-sample variation, which was found to be, on average, 1% RSD, and the mean value of the spectra was used in further data treatment. This procedure ensures the reduction in sample differences due to random variations in experimental conditions (e.g., aqueous volume, PIM contact area, PIM composition, etc.).

2.6. Development of the Chemometric Model

The multivariate calibration method was built from the information provided by the IR spectra of the PIMs with a set of 54 samples covering 27 different concentrations each by duplicate with an interval ranging from 1.92 × 10−7 to 1.92 × 10−6 mol/dm3 (10–100 ppb). The spectral range used for quantitative analysis initially was 400–5000 cm−1, so the experimental matrix consisted of 54 rows corresponding to the samples and 4601 wavelengths. Further cross-validation was performed in the calibration set using Venetian blinds with 3 data splits and 2 samples per blind (thickness). A test set consisting of 10 samples acquired as the calibration test, but independently, had dimensions of 10 × 4601. Metal concentrations given along the work refer to the initial concentrations of the metal in the solution. Data processing consisted of baseline correction (automatic weighted least squares, order one) and mean-centering. The PLS-Toolbox 9.0 software (Eigenvector Research, Inc. Wenatchee, WA, USA) was employed for all chemometric analyses. An in-house-made MATHLAB 9.7 R2019b (Natick, Apple Hill Campus, MA, USA) program was used for the evaluation of the FOM.

3. Results and Discussion

3.1. PIM Composition

The composition of the PIMs used in this work was that reported by Kozlowski [20], substituting the plasticizer o-nitrophenyl pentyl ether for NPOE. With this membrane composition, the extraction experiments carried out at pH 6 in a 0.01 mol/dm3 MES buffer solution for 2 h with an initial metal concentration of 1.35 × 10−4 mol/dm3 (7 ppm) showed 90.2% extraction of the metal. This result agrees with the expectations from SX experiments [19]. Consequently, due to the high extraction percentage, it was decided to keep such proportions of the polymer (22.2% ± 0.7) w/w%, plasticizer (54.4 ± 1.3) w/w%, and extractant (23.4 ± 1.5) w/w%.

3.2. Optimization of the Extraction Time

The extraction time necessary to reach the equilibrium was determined with a solution of 1.35 × 10−4 mol/dm3 Cr(VI), taking aliquots each 10 min for a total of 50 min. As observed in Figure 1, the concentration of the metal in the medium decreases until it reaches equilibrium, while the concentration of chromium in the membrane increases; thus, it can be concluded that the time required for the extraction reaction to reach equilibrium was 20 min. Consequently, further extraction experiments were carried out for 20 min of agitation at pH 6 with 0.01 mol/dm3 MES buffer to obtain the maximum extraction percentages.

3.3. Influence of pH

Knowing the time required to bring the extraction reaction to its equilibrium state, the influence of pH variations on metal extraction was studied. A concentration of 6.73 × 10−4 mol/dm3 (35 ppm) was used in the pH range 4 to 9, as this is usually the range found in natural waters [35]. In the experiments carried out at pH 4 and 5, a 0.01 mol/dm3 acetate/acetic acid buffer solution was used. For pH 6 and 7, a 0.01 mol/dm3 MES buffer was employed, and finally, for pH 8 and 9, a 0.01 mol/dm3 TRIS buffer was used. From Figure 2, it can be observed that chromium extraction decreases as pH increases in such a form that at pH 9, the extraction is less than 20%; however, at pH 4, the extraction exceeds 80%. This behavior is mainly due to the change in metal speciation (from HCrO4 to CrO42−) with the increase in pH of the solution, leading to a modification in the stoichiometry of the extraction reaction to neutralize the charge in the extracted metal anion, consistent with what is observed in liquid–liquid extraction [36]. This means that the form in which Cr(VI) is present in the system determines its extraction behavior. It is also clear from the graph that variations in pH values strongly affect E%, with the 4–7 interval being where changes of approximately 20% in E% are noticed. Considering that the PIM system will not be employed as a recovery method for the analyte, but for quantitative analysis, this implies that a constant pH value within this range can also be used with allowable results because E% remains constant and high enough to permit reproducible adsorption of the metal in the PIM. Therefore, even though the highest percentage of extraction was obtained at low pH, further experiments were performed at pH 6 since natural waters mostly register a pH around this value.

3.4. Influence of the Initial Metal Concentration

3.4.1. Adsorption Isotherm

Under the extraction conditions established, as mentioned, the effect of the initial metal concentration was studied. This study was performed using initial metal concentrations from 1.35 × 10−4 mol/dm3 (7 ppm) to 1.93 × 10−3 mol/dm3 (100 ppm). From Figure 3, it is seen that the extraction with high initial metal concentrations is less effective than the extraction with low metal concentrations, i.e., the lower the initial metal concentration the more effective the extraction. It is also observed that the time to rich equilibrium depends on the initial concentration of the metal; however, above 15 min, all systems reach such condition.
It is known that there is a quantity of metal that, when extracted, may saturate the active sites of the membrane, preventing it from continuing to extract the metal ion, meaning that the adsorption capacities depend on the different initial concentrations. In general, an adsorption isotherm is a curve that describes the phenomenon that governs the mobility of a substance from an aqueous medium to a solid phase at a constant temperature. Commonly, the mathematical correlation of this phenomenon is expressed graphically using the loading of the analyte in the solid phase (qe) as a function of the residual equilibrium concentration in the aqueous phase (Ce) [37]. Over the years, a wide variety of adsorption isotherm models have been developed to explain the kinetics, thermodynamics, and potential of this phenomenon [38]. The adsorption isotherm for the Cr(VI)/PIM is shown in Figure 4A.
From this, a Langmuir-type behavior is inferred. This model assumes monolayer analyte adsorption that occurs on localized sites that are equivalent without considering steric hindrance or adjacent interactions between adsorbed molecules, and it is described by the following equation [39]:
q e = q m a x K L C e 1 + K L   C e
where Ce is the solute concentration at equilibrium, qe is the amount of solute adsorbed at equilibrium, qmax is the maximum loading capacity, and KL is Langmuir’s constant. The previous equation can be linearized as:
C e q e = C e q m a x + 1 q m a x K L
So when plotting Ce/qe as a function of Ce, qmax and KL can be obtained from the abscissa and the slope of the line, respectively. Another parameter associated with the Langmuir adsorption model is the separation factor (RL) defined by [40]:
R L = 1 1 + K L   C 0
which indicates whether the nature of the adsorption is favorable or not (favorable (0 < RL < 1), linear (RL = 1), unfavorable (RL > 1), or irreversible (RL = 0)).
The linearized form of the adsorption isotherm is shown in Figure 4B. A determination coefficient of 0.999 indicates a good fit of the data to the Langmuir-type model. The parameters of the model are KL = 2199 cm3/mmol and qmax = 0.188 mmol/g. In Figure 4A, the points represent the data and the line the fitting to the model showing good agreement. Table 1 provides the parameters and proof of adequacy of the model. Furthermore, the RL values at the different initial concentrations are shown in Table 2. The 0 < RL < 1 values indicate that metal adsorption is favorable for this system [40].

3.4.2. Distribution Quotient

The distribution quotient (D) is defined as the ratio of the sum of the concentrations of all the chemical forms of the compound in the aqueous phase and in the organic phase according to:
D = [ Cr ( VI ) ] ¯ [ Cr ( VI ) ]
where the bar stands for the membrane phase. Figure 5A shows D as a function of the initial concentration of chromium from 1.35 × 10−4 to 1.93 × 10−3 mol/dm3. At low metal concentrations, D is ca10,000 units, which indicates a great affinity of the metal towards the membrane; on the other hand, a decay in its value at high concentrations is indicative of lower affinity, likely due to the saturation of the active sites, i.e., the extraction will be more effective at low metal concentrations. In general, the distribution into the organic phase of a metal complex increases with increasing temperature for complexes with significant hydrophobic character. The introduction of a complex into the organic phase involves several processes that can be associated with important changes in enthalpy (solvation processes) and entropy (solvent orientation and restructuring), leading to considerable temperature effects. In addition, to maintain electroneutrality and solute uphill pumping, the extraction systems require a coupling ion to be counter-transported along with the solute ion. Because the coupling ion must also cross the organic phase, it is bound to influence extraction efficiency [41].

3.4.3. Extraction Percent

In Figure 5B, the extraction percentages are presented for D values. However, if the volumes of the phases are considered, E%, may also be defined as:
E % = m m o l C r ¯ m m o l C r   0 × 100 %
At concentrations lower than 5 × 10−4 mol/dm3, E% is quantitatively high, with a value of 97.31%, and remains practically constant for the first 4 points. On the other hand, at higher concentrations, it decays, because of the occupation of the available sites on the membrane surface. From the point of view of quantitative analysis, initial concentrations lower than 5.5 × 10−4 mol/dm3 will be preferred to have constant extraction percentages, independent of the initial Cr(VI) content.

3.4.4. Enrichment Factor

The enrichment factor is a parameter that indicates how many times the metal is more concentrated in the membrane than in the solution, and is defined by [30]:
E = [ Cr ( VI ) ] ¯ [ Cr ( VI ) ] 0    
E was evaluated within the range 1.35 × 10−4–1.95 × 10−4 mol/dm3 (Figure 6). From the slope of the graph, E ≈ 18. However, the graph clearly shows two trends according to the concentration range below and above 6 × 10−4 mol/dm3 (Figure 6). Dividing the values according to this point, the highest enrichment factor is obtained at low concentrations, with E ≈ 33 where the extraction percentage is constant and practically independent of the initial chromium concentration. On the other hand, the enrichment factor at high concentrations corresponds to E ≈ 15, with the extraction percentage depending on the initial chromium concentration. This observation highlights the need to work in the lower range of concentrations to obtain better results, since the higher the metal is loaded in the PIM the more sensitive the analytical method. The breakpoint noted in Figure 6 is related to the saturation of the active sites of the PIM, as previously discussed; where they are limited, their fulfillment is dependent on the initial metal concentration due to the competition from them established by the ions in solution. In contrast, when they are not limited, this competition is reduced due to the excess of sites available.

3.5. Determination of the Extraction Equilibrium

For this determination, the concentration of Aliquat 336 in the membrane was varied, keeping the amount of CTA and NPOE constant. The mass values of the Aliquat 336 used were 1.5, 5.1, 9, 11.1, 12, 15, 22.5, and 30 mg. According to Kebiche-Senhadji [23], extraction occurs by the following reaction:
H i C r O 4 n + n L + C l ¯ H i C r O 4 n L n + ¯ + n C l
where L + C l ¯ stands for the extractant and H i C r O 4 n L n + ¯ for the extracted species in the PIM phase. The extraction constant, Kext, is defined by:
K e x t = [ H i C r O 4 n L n + ] ¯   [ C l ] n [ H i C r O 4 n ]   [ L + C l ¯ ] n
Considering the definition of D, Kext can be rewritten as follows:
K e x t = D   [ C l ] n [ L + C l ¯ ] n
Taking logarithms to both sides of the equation and rearranging:
l o g D = l o g K e x t + n l o g [ L + C l ¯ ] n l o g [ C l ]
Figure 7 shows logD as a function of the unreacted log[Aliquat 336]; a linear relationship with a determination coefficient of 0.9902 was obtained. The slope value of 1.16 indicates that the Cr(VI): Aliquat 336 ratio is 1:1, so the final extraction reaction becomes:
H C r O 4 + L + C l ¯ H C r O 4 L + ¯ + C l
This result is consistent with that reported in the literature, where the same stoichiometry was observed in SX [42].

3.6. Chemometric Analyses

To meet the chromium determination requirements in aqueous media [3,4,6] further experiments were performed at a lower concentration range (10–100 ppb) where, according to the previous discussions, the PIM system will perform the best.

3.6.1. PCA

The sample score plots of the PCA analysis are shown in Figure 8A,B. A five-component model accounted for 95.38% of the variance in the X-block with 82.5% accounting exclusively for PC1. The RMSEC and RMSECV values were 0.10036 and 0.126592, respectively. A distinctive V-shape pattern observed in structured dependent data was identified in the plot PC1 vs. PC2 (Figure 8A). In gene population studies, PCA plots often appear triangular due to the underlying genetic population structure and the way genetic variation is distributed across individuals, as individuals with similar genetic ancestry tend to cluster together, at the three corners of the triangle, which represent different groups of individuals with unique genetic backgrounds [43,44]. This pattern also occurs in chemistry mixture problems with the three vertices being samples that contain a single component each, the samples falling on a line between two vertices being binary mixtures of the three analytes, and the remaining points being ternary mixtures of the three analytes [45]. Similar behavior has also been observed in optodes with 1-(2-pyridylazo)-2-naphthol as a chromophore where the proportions of the analytes (Cu(II), Zn(II), and Pb(II)) were systematically varied giving rise to mixtures of colors in the PIM [46]. However, in the present case, the pattern seems to be more complex as observed when analyzing the contribution of the first three PCs (Figure 8B), where a progressive change in the spectra apparently not related to the concentration of the analyte is observed. Although a complete understanding of factors giving rise to the observed shape is outside the present work, it can be related to dependent structured modifications in the PIM medium leading to non-linearities of the system, as the polar nature of water molecules adsorbed can induce local electric dipoles resulting in an increase in the effective polarizability of the material, affecting the penetration depth of the evanescent wave into the sample.

3.6.2. Selection of the Spectral Wavelength Range

Preliminary full-spectrum PLS modeling did not give satisfactory results, as a 10 latent variable model gave RMSEC of 8.41851, RMSECV of 36.9153, Bias of 2.84217 × 10−14, CV Bias of 0.208479, R2 Cal of 0.905564, and R2 CV of 0.000624484, clearly showing model overfitting. This result is somehow expected from the previous PCA analysis where no evident trend with Cr(VI) concentrations was observed. Furthermore, no improvement was observed by changing the preprocessing method or the number of latent variables. However, high improvement was attained when some regions of the spectra were discarded from the calculations, so the forward interval PLS algorithm (iPLS), a variable selection method, was implemented. iPLS selects a subset of variables, which will give superior prediction compared to using all the variables in the dataset. It performs a sequential, exhaustive search for the best variable or combination of variables [47]. Figure 9 shows the obtained results. It was observed that when wavelengths in the 3050–3890 cm−1 region were considered (in red color in Figure 9), the best predictions were attained, with a minimum RMSECV value when bands in the ranges of 3451–3500 cm−1 and 3751–3800 cm−1 were chosen (in green in Figure 9).
The extraction mechanism of Cr(VI) with Aliquat 336 involves ion-pair formation. Aliquat 336, being a cationic extractant, can form an ion pair with the hydrogen chromate ion through electrostatic interactions. The broad band in the FTIR spectra with the maximum located near 3500 cm−1 can be attributed both to hydroxyl stretching in CTA [47] and to stretching vibrations of adsorbed water in Aliquat 336 [22,28,48,49]. Fontàs et al. have found that modifications in the surface composition of water-equilibrated Aliquat 336-PIMs may be associated with the solubilization of the extractant in the water solution, which, therefore, may affect the reactivity of the membrane’s surface, but not the bulk properties of the membrane itself [50]. Such Aliquat 336 solubilized molecules may form micelles at the PIM [34], containing hydrated nonpolar cavities, such that the Me3N+ headgroup in Aliquat 336 and the HCrO4 anion likely form solvent-separated ion pairs at the micelle surface [51]. This explains why this IR region performs the best in the developed method and can also be the reason for similar behavior, in which changes in such a band range can be noticed after the adsorption of Cr(VI) in an Aliquat 336 Dowex 1 × 8 impregnated resin [52] and the band shifting in the -OH stretching region reported for a CTA/1,5-diphenylcarbazide/Aliquat 336 optode [14].

3.6.3. PLS Modeling

Hence, the final PLS model included the 100 wavelengths selected by iPLS and 10 latent variables according to the RMSEC and RMSECV values obtained using the calibration and cross-validation results. In Table 3, the percent of variance accounting for the different latent variables is shown. With the 10 latent variables, the variance captured by the model in the X-block reaches 99.99% while that in the y-block (Cr(VI) concentrations) is 97.17%.
A summary of this final model is graphically presented in Figure 10A–D. No important outliers were detected in the data (Figure 10A,B) and the first two latent variables accounted for 99.24% of the variance in the X-block with almost all data being inside the Hotelling T2 ellipse (Figure 10D). Interestingly, the score plot did not show the complex pattern previously observed in PCA analysis this time, indicating efficient removal of variance not related to Cr(VI) concentrations modification, which this time moves sequentially from the third to the first quadrant in the plot as concentration increases (Figure 10D). The parity plot, i.e., measured vs. predicted concentrations (Figure 10C), showed excellent parameter values with RMSEC of 3.73115, RMSECV of 6.82685, Bias of −1.91847 × 10−13, CV Bias of 0.185947, R2 Cal of 0.98145, and R2 CV of 0.940902.

3.7. Model Validation and Application

3.7.1. Figures of Merit (FOM) [53,54,55]

When an analytical method is proposed, it is necessary to investigate whether it will perform adequately under the conditions where it will be further applied. This procedure is referred to as the validation of the method and it is performed through the determination of the FOM. These are very important parameters in characterizing, comparing, and developing new multivariate methods. Many of these figures are closely related to the net analyte signal (NAS) concept, as presented by Lorber [56] and defined by:
N A S i = ( x i · b ) · ( b T · b ) 1 · b T
where xi is a sample spectrum after preprocessing and b is a column vector of the PLS regression coefficients. The concept of NAS arises naturally in multivariate calibration from the fact that a predicted sample spectrum can have various contributions from all the sample components. Therefore, it is logical to decompose the spectrum into two orthogonal parts: A part that can be exclusively assigned to the analyte of interest (NAS) and the remaining part that contains the contribution, possibly variable, of other components [57]. This decomposition is carried out by the regression algorithms in such a way that the NAS is proportional to the concentration of the analyte of interest. Since the NAS is the only part of the spectrum that is used for prediction, no information is lost by transforming the NAS vector into a scalar. The natural choice is to take the Euclidean norm, i.e., its length so that the scalar NAS is obtained as r* = ||r*||, where r is an arbitrary vector. Using the NAS, a multivariate calibration model can be represented on a pseudo-univariate plot, with this representation being exact, not approximate. In other words, the calibration curve is obtained by replacing the measured instrumental signal used in univariate calibration (e.g., absorbance at a single wavelength) with the NAS (r*) allowing a simpler interpretation of the response signal [53].

Accuracy (RMSE)

This parameter indicates the closeness of the relationship between the reference value and that found by the model. In multivariate calibration, it is usually expressed as the root mean square error of calibration (RMSEC), root mean square error of cross-validation (RMSECV), and root mean square error of prediction (RMSEP), according to the following equation:
R M S E = i = 1 n ( y i y i ^ ) 2 n
where y i and y i ^ are the estimated value of the model and the reference value of sample i, respectively, and n is the number of samples. The estimated value varied according to what is measured, the calibration (RMSEC), cross-validation (RMSECV), or the test (RMSEP) results.

Selectivity (sel)

This indicates the part of the total signal that is not lost due to spectral overlapping and can be defined in the multivariate context through the NAS calculation:
s e l = | | s k * | | | | s k | |
where | | s k | | stands for the norm of the sensitivity coefficients of the spectra containing the analyte k at unit concentration, and | | s k * | | for that corresponding to its NAS.

Sensitivity (sen)

Sensitivity measures the changes in the response, as a function of the concentration of a particular analyte, and is given by the following equation:
s e n = | | s k * | |
However, it can also be calculated without the use of NAS theory according to the equation:
s e n = 1 | | b | |
A more useful FOM is the analytical sensitivity (γ), which is defined, by analogy with univariate calibration, as the ratio between sen and the instrumental noise (δx). The inverse of γ (γ−1) provides an estimation of the minimum concentration difference that is discernible by the analytical method considering the random experimental noise as the only source of error, regardless of the specific technique employed [57].

Limit of Detection (LD)

Following the IUPAC recommendations, LD can be defined as the minimum detectable value of the net signal for which the false negative (β) and false positive (α) probabilities are 0.05. LD can be calculated analogously, as for univariate calibration, according to the equation:
L D = 3.3 δ x 1 s e n

Limit of Quantitation (LQ)

Quantizability is generally expressed in terms of the signal or analyte concentration value that will produce estimates with a specified relative standard deviation, usually 10% RSD. Following the same assumptions described above, the quantitation limit in the multivariate calibration can be calculated by:
L Q = 10 δ x 1 s e n
The calculated FOM of the PLS model is given in Table 4. The absence of systematic errors in accuracy determination was verified with an F-test over the parity plot during cross-validation (Figure 10C). In the absence of errors, it is expected that the slope obtained from this graph is equal to 1 (β1 = 1) and the abscissa is 0 (β0 = 0) (null hypothesis). The alternative hypothesis states that at least one of these parameters is different from the expected value.
This simultaneous null hypothesis can be tested from [30,56]:
F = ( β 0 b 0 ) 2 + 2 x ¯ ( β 0 b 0 ) ( β 1 b 1 ) + ( x i 2 / n ) ( β 1 b 1 ) 2 2 ( S e 2 / n )
where:
  • β0 = 0,
  • β1 = 1,
  • b0 = observed abscissa,
  • b1 = observed slope,
  • and S e = ( y i y ^ i ) 2 n 2 .
The obtained F-value is then compared to an F distribution reference value with 2 and n − 2 degrees of freedom at the chosen significance level.
From Figure 10C, b0 = 4.22147, b1 = 0.924373, x ¯ = 5 4.4814, n = 54, and Se = 5.916 giving F = 3.32. As this value is smaller than the tabulated F0.05 2, 52 = 3.96, the hypothesis of β1 = 1, β0 = 0 is accepted, i.e., no systematic bias is present.
The model presents a sensitivity value of 0.0015%, indicating that it is capable of distinguishing samples with concentration differences of 0.6 ppb. The selectivity value indicates that 1.5% of the information on the analyte contained in the sample is orthogonal to the interference space. In other words, the selectivity indicates that approximately 98.5% of the interfering analytical signals were removed during the NAS calculation. This result is expected due to the low concentration range of the analyte in relation to the other PIM components, meaning that only slight changes in the PIM spectra were related to Cr(VI) concentration variation. However, even such low variability was effectively accounted by the PLS algorithm. On the other hand, the well-probed chemical selectivity of Aliquat 336 for Cr(VI) over Co(II), Cd(II), Pb(II), Ni(II), Cu(II), Zn(II), and Fe(II) [23,27,28,58] favors the testing of the PIM on real samples.
Table 5 compares the results of the developed sensor to those reported in the literature for Cr(VI) analysis with comparable systems. Similar FOM with other Cr(VI) optical quantitation methods were found. However, as most of the works reported in Table 5 rely on the use 1,5-diphenylcarbazide as a chromophore, they are limited to work at low pH values, where the highest absorbance values are attained [13,14], with those methods then being not suitable for in-situ analysis of water samples. In comparison to some of them, the developed PIM sensor is, in many cases, much easier to implement as it requires few manipulations and a reduced number of chemical compounds.

3.7.2. Application

The developed PLS method was applied to each of the test samples, and the t-value for the comparison between the measured (reference) and predicted (found) concentrations was calculated. The results are shown in Table 6. As the obtained p-value > 0.05, there were no significant differences between both quantities at the 95% confidence level. The RMSEP value given in the same table is close to the RMSEC and RMSECV values, meaning that model overfitting was prevented. An average recovery of 104.02 ± 4.12 (α = 0.05) was obtained.

4. Conclusions

Using infrared spectroscopy and multivariate statistical analysis, a MID-FTIR-PLS PIM-based sensor for the quantitative determination of Cr(VI) from an aqueous medium was developed with satisfactory results. Optimized conditions for the system consisted of (22.2 ± 0.7) w/w%, CTA, (54.4 ± 1.3) w/w% NPOE and (23.4 ± 1.5) w/w% Aliquat 336, pH 6, and 20 min of extraction time, although pH values in the range of 4–7 could also be used. It was observed that chromium extraction decreases as the pH increases, in such a form that at pH 9, the extraction is less than 20%. However, at pH 4, the extraction exceeds 80%; this behavior is mainly due to the change in metal speciation (from HCrO4 to CrO42−) with increasing pH of the solution. Cr(VI) adsorption followed a Langmuir-type isotherm with KL = 2199 cm3/mmol, qmax = 0.188 mmol/g and 0 < RL < 1. At low metal concentrations, D is ca.10,000 units, which indicates a great affinity of the metal towards the membrane; on the other hand, a decay in its value at high concentrations is indicative of lower affinity, as also observed in the E% profile. The highest enrichment factor was obtained at low analyte concentrations, with E ≈ 33, where the extraction percentage is practically constant. On the other hand, the enrichment factor at high concentrations corresponds to E ≈ 15. Characterization of the extraction reaction indicated a 1:1 Cr(VI): Aliquat 336 ratio. The PCA analysis of the PIMs revealed a complex pattern, which was satisfactorily simplified and related to Cr(VI) concentrations through the use of a variable selection method (iPLS) in which bands in the ranges of 3451–3500 cm−1 and 3751–3800 cm−1 were chosen. The final PLS model including the 100 wavelengths selected by iPLS and 10 latent variables shows excellent FOM values. The developed PIM sensor is suitable for in situ analysis of aqueous samples, and it is easy to implement as it requires few manipulations and a reduced number of chemical compounds in comparison to other similar reported systems.

Author Contributions

Conceptualization, A.M.d.l.P. and E.R.d.S.M.; methodology, A.M.d.l.P.; formal analysis, A.M.d.l.P. and E.R.d.S.M.; investigation, A.M.d.l.P.; resources, E.R.d.S.M.; writing—original draft preparation, A.M.d.l.P.; writing—review and editing, E.R.d.S.M. and J.d.G.; supervision, E.R.d.S.M. and J.d.G.; project administration, E.R.d.S.M.; funding acquisition, E.R.d.S.M. and J.d.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by CONACYT (project 2011CB167148) and DGAPA-UNAM (project PAPIIT IN229219).

Institutional Review Board Statement

Not applicable.

Data Availability Statement

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

Acknowledgments

Q. Nadia Marcela Munguía Acevedo and Q.F.B. María Guadalupe Espejel Maya are thanked for their help in technical services. A.M.d.l.P. acknowledges CONACYT for his MsSc scholarship and the Consejo Mexiquense de Ciencia y Tecnología for the financial support provided through its scholarship-thesis program.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

Aliquat 336Methyltrioctylammonium chloride
CalCalibration
CTACellulose triacetate
CVCross-validation
DDistribution coefficient
EEnrichment factor
E%Extraction percentage
FAASFlame atomic absorption spectroscopy
FOMFigures of merit
FTIRFourier-transform infrared spectroscopy
iPLSInterval PLS
IRInfrared spectroscopy
ISOsIon-selective optodes
KextExtraction constant
LDLimit of detection
LQLimit of quantitation
MES2-(N-morpholino) ethanesulfonic acid
MID-FTIRMid Fourier transform infrared spectroscopy
NASNet analyte signal
NPOE2-nitrophenyl octyl ether
PCAPrincipal component analysis
PIM(s)Polymer inclusion membrane(s)
PLSPartial least squares
R2Determination coefficient
RMSECRoot mean square error of calibration
RMSECVRoot mean square error of cross-validation
RMSEPRoot mean square error of prediction
RSDRelative standard deviation
selSelectivity
senSensitivity
SIRSolvent-impregnated resins
SLMSupported liquid membranes
SXSolvent extraction
TRISTris(hydroxymethyl)aminomethane
VISVisible spectroscopy

References

  1. Georgaki, M.-N.; Charalambous, M.; Kazakis, N.; Talias, M.A.; Georgakis, C.; Papamitsou, T.; Mytiglaki, C. Chromium in Water and Carcinogenic Human Health Risk. Environments 2023, 10, 33. [Google Scholar] [CrossRef]
  2. Bhattacharya, M.; Dutta, S.K.; Sikder, J.; Mandal, M.K. Computational and experimental study of chromium (VI) removal in direct contact membrane distillation. J. Membr. Sci. 2014, 450, 447–456. [Google Scholar] [CrossRef]
  3. Chromium in Drinking Water. United States Environmental Protection Agency (EPA). Available online: https://www.epa.gov/sdwa/chromium-drinking-water (accessed on 13 June 2023).
  4. Chromium in Drinking-Water. World Health Organization. Available online: https://cdn.who.int/media/docs/default-source/wash-documents/wash-chemicals/chromium.pdf?sfvrsn=37abd598_6 (accessed on 13 June 2023).
  5. Tumolo, M.; Ancona, V.; De Paola, D.; Losacco, D.; Campanale, C.; Massarelli, C.; Uricchio, V.T. Chromium Pollution in European Water, Sources, Health Risk, and Remediation Strategies: An Overview. Int. J. Environ. Res. Public Health 2020, 17, 5438. [Google Scholar] [CrossRef] [PubMed]
  6. Vaiopoulou, E.; Gikas, P. Regulations for chromium emissions to the aquatic environment in Europe and elsewhere. Chemosphere 2020, 254, 126876. [Google Scholar] [CrossRef] [PubMed]
  7. Parks, J.L.; McNeill, L.; Frey, M.; Eaton, A.D.; Haghani, A.; Ramirez, L.; Edwards, M. Determination of total chromium in environmental water samples. Water Res. 2004, 38, 2827–2838. [Google Scholar] [CrossRef]
  8. Capitán-Vallvey, L.F.; Palma, A.J. Recent developments in handheld and portable optosensing—A review. Anal. Chim. Acta 2011, 696, 27–46. [Google Scholar] [CrossRef]
  9. Vandenabeele, P.; Castro, K.; Hargreaves, M.; Moens, L.; Madariaga, J.M.; Edwards, H.G.M. Comparative study of mobile Raman instrumentation for art analysis. Anal. Chim. Acta 2007, 588, 108–116. [Google Scholar] [CrossRef]
  10. Du, X.; Xie, X. Ion-Selective optodes: Alternative approaches for simplified fabrication and signaling. Sens. Actuators B Chem. 2021, 335, 129368. [Google Scholar] [CrossRef]
  11. Scindia, Y.M.; Pandey, A.K.; Reddy, A.V.R.; Manohar, S.B. Chemically selective membrane optode for Cr(VI) determination in aqueous samples. Anal. Chim. Acta 2004, 515, 311–321. [Google Scholar] [CrossRef]
  12. Carrington, N.A.; Thomas, G.H.; Rodman, D.L.; Beach, D.B.; Xue, Z.-L. Optical determination of Cr(VI) using regenerable, functionalized sol-gel monoliths. Anal. Chim. Acta 2007, 581, 232–240. [Google Scholar] [CrossRef]
  13. Castilleja-Rivera, W.L.; Hinojosa-Reyes, L.; Guzmán-Mar, J.L.; Hernández-Ramírez, A.; Ruíz-Ruíz, E.; Cerda, V. Sensitive determination of chromium(VI)in paint samples using a membrane optode coupled to a multisyringe flow injection system. Talanta 2012, 99, 730–736. [Google Scholar] [CrossRef]
  14. Munandar, R. A Cellulose Triacetate-Based Optode Membrane for Cr(VI) Detection in Water: Fabrication and Characterization. Bachelor’s Thesis, IPB University, Bogor, Indonesia, 2022. [Google Scholar]
  15. Keskin, B.; Zeytuncu-Gökoğlu, B.; Koyuncu, I. Polymer inclusion membrane applications for transport of metal ions: A critical review. Chemosphere 2021, 279, 130604. [Google Scholar] [CrossRef] [PubMed]
  16. Ulewicz, M. Polymer Membranes in Separation Process. Membranes 2021, 11, 831. [Google Scholar] [CrossRef] [PubMed]
  17. Almeida, M.I.G.S.; Cattrall, R.W.; Kolev, S.D. Polymer inclusion membranes (PIMs) in chemical analysis—A review. Anal. Chim. Acta 2017, 987, 1–14. [Google Scholar] [CrossRef] [PubMed]
  18. Salazar, E.; Ortiz, M.I.; Urtiaga, A.M.; Irabien, J.A. Equilibrium and Kinetics of Cr(V1) Extraction with Aliquat 336. Ind. Eng. Chem. Res. 1992, 31, 1516–1522. [Google Scholar] [CrossRef]
  19. Senol, A. Amine extraction of chromium(VI) from aqueous acidic solutions. Sep. Purif. Technol. 2004, 36, 63–75. [Google Scholar] [CrossRef]
  20. Kozlowski, C.; Walkowiak, W. Applicability of liquid membranes in chromium(VI) transport with amines as ion carriers. J. Membr. Sci. 2005, 266, 143–150. [Google Scholar] [CrossRef]
  21. Kabaya, N.; Arda, M.; Saha, B.; Streat, M. Removal of Cr(VI) by solvent impregnated resins (SIR) containing aliquat 336. React. Funct. Polym. 2003, 54, 103–115. [Google Scholar] [CrossRef]
  22. Wojcik, G.; Hubicki, Z.; Ryczkowski, J. Investigations of Chromium(III) and (VI) Ions Sorption on SIR by using Photoacoustic and DRS Methods. Acta Phys. Pol. A 2009, 116, 432–434. [Google Scholar] [CrossRef]
  23. Kebiche-Senhadji, O.; Tingry, S.; Seta, P.; Benamor, M. Selective extraction of Cr(VI) over metallic species by polymer inclusion membrane (PIM) using anion (Aliquat 336) as carrier. Desalination 2010, 258, 59–65. [Google Scholar] [CrossRef]
  24. Gherasim, C.V.; Bourceanu, G. Removal of chromium(VI) from aqueous solutions using a polyvinyl-chloride inclusion membrane: Experimental study and modelling. Chem. Eng. J. 2013, 220, 24–34. [Google Scholar] [CrossRef]
  25. Kagaya, S.; Maeno, T.; Ito, K.; Gemmei-Ide, M.; Cattrall, R.W.; Kolev, S.D. Improvement of chromium(VI) extraction from acidic solutions using a poly(vinyl chloride)-based polymer inclusion membrane with Aliquat 336 as the carrier. Anal. Sci. 2017, 33, 643–646. [Google Scholar] [CrossRef]
  26. Semghouni, H.; Bey, S.; Figoli, A.; Criscuoli, A.; Russo, F.; Mohamed, B.; Drioli, E. Chromium(VI) Removal by Polyvinyl Chloride(PVC)/Aliquat-336 Polymeric Inclusion Membranes in a Multiframe Flat Sheet Membrane Module. Appl. Sci. 2019, 9, 2994. [Google Scholar] [CrossRef]
  27. Sellami, F.; Kebiche-Senhadji, O.; Marais, S.; Couvrat, N.; Fatyeyeva, K. Polymer inclusion membranes based on CTA/PBAT blend containing Aliquat 336 as extractant for removal of Cr(VI): Efficiency, stability and selectivity. React. Funct. Polym. 2019, 139, 120–132. [Google Scholar] [CrossRef]
  28. Sellami, F.; Kebiche-Senhadji, O.; Marais, S.; Colasse, L.; Fatyeyeva, K. Enhanced removal of Cr(VI) by polymer inclusion membrane based on poly(vinylidene fluoride) and Aliquat 336. Sep. Purif. Technol. 2020, 248, 117038. [Google Scholar] [CrossRef]
  29. Bahrami, S.; Reza Yaftian, M.; Najvak, P.; Dolatyari, L.; Shayani-Jam, H.; Kolev, S.D. PVDF-HFP based polymer inclusion membranes containing Cyphos® IL 101 and Aliquat® 336 for the removal of Cr(VI) from sulfate solutions. Sep. Purif. Technol. 2020, 250, 117251. [Google Scholar] [CrossRef]
  30. González-Albarrán, R.; de Gyves, J.; Rodríguez de San Miguel, E. Determination of Cadmium (II) in Aqueous Solutions by In Situ MID-FTIR-PLS Analysis Using a Polymer Inclusion Membrane-Based Sensor: First Considerations. Molecules 2020, 25, 3436. [Google Scholar] [CrossRef]
  31. Mohd Suah, F.B.; Ahmad, M.; Heng, L.Y. A novel polymer inclusion membrane based optode for sensitive determination of Al3+. Spectrochim. Acta A 2015, 144, 81–87. [Google Scholar] [CrossRef]
  32. Ngarisan, N.I.; Zanariah, C.W.; Ngah, C.W.; Ahmad, M.; Kuswandi, B. Optimization of polymer inclusion membranes (PIMs) preparation for immobilization of Chrome Azurol S for optical sensing of aluminum (III). Sens. Actuator B Chem. 2014, 203, 465–470. [Google Scholar] [CrossRef]
  33. Mohd Suah, F.B.; Ahmad, M.; Heng, L.Y. Highly sensitive fluorescence optode for aluminium(III) based on non-plasticized polymer inclusion membrane. Sens. Actuators B Chem. 2014, 201, 490–495. [Google Scholar] [CrossRef]
  34. O’Rourke, M.; Duffy, N.; De Marco, R.; Potter, I. Electrochemical Impedance Spectroscopy—A Simple Method for the Characterization of Polymer Inclusion Membranes Containing Aliquat 336. Membranes 2011, 1, 132–148. [Google Scholar] [CrossRef] [PubMed]
  35. Health Canada. Guidelines for Canadian Drinking Water Quality: Guideline Technical Document—pH. Water and Air Quality Bureau, Healthy Environments and Consumer Safety Branch, Health Canada, Ottawa, Ontario, 2015. (Catalogue No H144-28/2016E-PDF). Available online: https://www.canada.ca/content/dam/hc-sc/documents/services/publications/healthy-living/guidelines-canadian-drinking-water-quality-guideline-technical-document-ph-eng.pdf (accessed on 15 July 2023).
  36. Oh, J.-K.; Lee, H.-Y.; Kim, S.-G. Solvent Extraction of Chromium from Dilute Aqueous Solution. Geosyst. Eng. 2000, 3, 49–52. [Google Scholar] [CrossRef]
  37. Ncibi, M.C. Applicability of some statistical tools to predict optimum adsorption isotherm after linear and non-linear regression analysis. J. Hazard. Mater. 2008, 153, 207–212. [Google Scholar] [CrossRef]
  38. Malek, A.; Farooq, S. Comparison of isotherms models for hydrocarbon adsorption on activated carbon. AIChE J. 1996, 42, 3191–3201. [Google Scholar] [CrossRef]
  39. Foo, K.Y.; Hameed, B.H. Insights into the modeling of adsorption isotherm systems. Chem. Eng. J. 2010, 156, 2–10. [Google Scholar] [CrossRef]
  40. Zheng, H.; Liu, D.; Zheng, Y.; Liang, S.; Liu, Z. Sorption isotherm and kinetic modeling of aniline on Cr-bentonite. J. Hazard. Mater. 2009, 167, 141–147. [Google Scholar] [CrossRef] [PubMed]
  41. Kislik, V.S. Solvent Extraction: Classical and Novel Approaches, 1st ed.; Elsevier: Amsterdam, The Netherlands, 2012; pp. 44–45. [Google Scholar]
  42. Senol, A. Extractive removal of Cr(VI) from aqueous acidic media by aliquat 336/xylene system: Optimization and modelling of equilibrium. Indian J. Chem. Technol. 2017, 24, 269–283. [Google Scholar]
  43. Bulbul, O.; Filoglu, G.; Zorlu, T.; Altuncul, H.; Freire-Aradas, A.; Söchtig, J.; Ruiz, Y.; Klintschar, M.; Triki-Fendri, S.; Rebai, A.; et al. Inference of biogeographical ancestry across central regions of Eurasia. Int. J. Leg. Med. 2015, 130, 73–79. [Google Scholar] [CrossRef]
  44. Tay, G.K.; Henschel, A.; Daw Elbait, G.; Al Safar, H.S. Genetic Diversity and Low Stratification of the Population of the United Arab Emirates. Front. Genet. 2020, 11, 608. [Google Scholar] [CrossRef]
  45. El Haddad, J.; de Miollis, F.; Bou Sleiman, J.; Canioni, L.; Mounaix, P.; Bousquet, B. Chemometrics applied to quantitative analysis of ternary mixtures by Terahertz spectroscopy. Anal. Chem. 2014, 86, 4927–4933. [Google Scholar] [CrossRef]
  46. Mancilla-Rico, A.J.; Rodríguez de San Miguel, E. Simultaneous determination of Cu(II), Zn(II), and Pb(II) from aqueous solutions using a polymer inclusion membrane (PIM) based-sensor with 1-(2-pyridylazo)-2-naphthol (PAN) as chromophore and chemometric methods. Front. Anal. Sci. 2022, 2, 971352. [Google Scholar] [CrossRef]
  47. Nørgaard, L.; Saudland, A.; Wagner, J.; Nielsen, J.P.; Munck, L.; Engelsen, S.B. Interval Partial Least-Squares Regression (iPLS): A Comparative Chemometric Study with an Example from Near-Infrared Spectroscopy. App. Spectrosc. 2000, 54, 413–419. [Google Scholar] [CrossRef]
  48. Nunes, S.; Ramacciotti, F.; Neves, A.; Angelin, E.M.; Ramos, A.M.; Roldão, E.; Wallaszkovits, N.; Armijo, A.A.; Melo, M.J. A diagnostic tool for assessing the conservation condition of cellulose nitrate and acetate in heritage collections: Quantifying the degree of substitution by infrared spectroscopy. Herit. Sci. 2020, 8, 33. [Google Scholar] [CrossRef]
  49. Naz, G.; Othaman, Z.; Shamsuddin, M.; Ghoshal, S.K. Aliquat 336 stabilized multi-faceted gold nanoparticles with minimal ligand density. App. Surf. Sci. 2016, 363, 74–82. [Google Scholar] [CrossRef]
  50. Fontàs, C.; Vera, R.; Anticó, E.; Martínez de Yuso, M.d.V.; Rodríguez-Castellón, E.; Benavente, J. New Insights on the Effects of Water on Polymer Inclusion Membranes Containing Aliquat 336 Derivatives as Carriers. Membranes 2022, 12, 192. [Google Scholar] [CrossRef] [PubMed]
  51. van der Vegt, N.F.A.; Haldrup, K.; Roke, S.; Zheng, J.; Lund, M.; Bakker, H.J. Water-Mediated Ion Pairing: Occurrence and Relevance. Chem. Rev. 2016, 116, 7626–7641. [Google Scholar] [CrossRef]
  52. Kalidhasan, S.; Santhana, A.; Kumar, K.; Rajesh, V.; Rajesha, N. An efficient ultrasound assisted approach for the impregnation of room temperature ionic liquid onto Dowex 1×8 resin matrix and its application toward the enhanced adsorption of chromium (VI). J. Hazard. Mater. 2012, 213–214, 249–257. [Google Scholar] [CrossRef]
  53. Olivieri, A.C.; Faber, N.M.; Ferré, J.; Boqué, R.; Kalivas, J.H.; Mark, H. Uncertainty estimation and figures of merit for multivariate calibration (IUPAC Technical Report). Pure Appl. Chem. 2006, 78, 633–661. [Google Scholar] [CrossRef]
  54. Braga, J.W.B.; Trevizan, L.C.; Nunes, L.C.; Rufini, I.A.; Santos, D., Jr.; Krug, F.J. Comparison of univariate and multivariate calibration for the determination of micronutrients in pellets of plant materials by laser induced breakdown spectrometry. Spectrochim. Acta B 2010, 65, 66–74. [Google Scholar] [CrossRef]
  55. de Carvalho Rocha, W.F.; Nogueira, R.; Vaz, B.G. Validation of model of multivariate calibration: An application to the determination of biodiesel blend levels in diesel by near-infrared spectroscopy. J. Chemometr. 2012, 26, 456–461. [Google Scholar] [CrossRef]
  56. Lorber, A. Error propagation and figures of merit for quantification by solving matrix equations. Anal. Chem. 1986, 58, 1167–1172. [Google Scholar] [CrossRef]
  57. Ferreira, M.H.; Braga, J.W.B.; Sena, M.M. Development and validation of a chemometric method for direct determination of hydrochlorothiazide in pharmaceutical samples by diffuse reflectance near infrared spectroscopy. Microchem. J. 2013, 109, 158–164. [Google Scholar] [CrossRef]
  58. Gherasim, C.V.; Bourceanu, G.; Olariu, R.J.; Arsene, C. A novel polymer inclusion membrane applied in chromium (VI) separation from aqueous solutions. J. Hazard. Mater. 2011, 197, 244–253. [Google Scholar] [CrossRef] [PubMed]
  59. Sreevani, I.; Reddy, P.R.; Reddy, V.K. Rapid and Simple Spectrophotometric Determination of Traces of Chromium (VI) in Wastewater Samples and in Soil samples by using 2-Hydroxy, 3-Methoxy Benzaldehyde Thiosemicarbazone (HMBATSC). IOSR J. App. Phys. 2013, 3, 40–45. [Google Scholar] [CrossRef]
  60. Lace, A.; Ryan, D.; Bowkett, M.; Cleary, J. Chromium Monitoring in Water by Colorimetry Using Optimised 1,5-Diphenylcarbazide Method. Int. J. Environ. Res. Public Health 2019, 16, 1803. [Google Scholar] [CrossRef]
  61. Revanasiddappa, H.D.; Kiran Kumar, T.N. Spectrophotometric Determination of Trace Amounts of Chromium with Citrazinic Acid. J. Anal. Chem. 2001, 56, 1084–1088. [Google Scholar] [CrossRef]
  62. Sun, X.; Li, B.; Qi, A.; Tian, C.; Han, J.; Shi, Y.; Lin, B.; Chen, L. Improved assessment of accuracy and performance using a rotational paper-based device for multiplexed detection of heavy metals. Talanta 2018, 178, 426–431. [Google Scholar] [CrossRef]
Figure 1. Cr(VI) extraction profile as a function of time (pH 6 in a 0.01 mol/dm3 MES buffer solution, [Cr(VI)]0 = 1.35 × 10−4 mol/dm3).
Figure 1. Cr(VI) extraction profile as a function of time (pH 6 in a 0.01 mol/dm3 MES buffer solution, [Cr(VI)]0 = 1.35 × 10−4 mol/dm3).
Membranes 13 00740 g001
Figure 2. Cr(VI) extraction profile as a function of pH of the aqueous sample. [Cr(VI)]0 = 6.73 × 10−4 mol/dm3.
Figure 2. Cr(VI) extraction profile as a function of pH of the aqueous sample. [Cr(VI)]0 = 6.73 × 10−4 mol/dm3.
Membranes 13 00740 g002
Figure 3. Cr(VI) extraction profile for different initial concentrations of metal at pH 6.
Figure 3. Cr(VI) extraction profile for different initial concentrations of metal at pH 6.
Membranes 13 00740 g003
Figure 4. (A) Cr(VI) adsorption isotherm at pH 6. (B) Linearized form of the Cr(VI) Langmuir adsorption isotherm.
Figure 4. (A) Cr(VI) adsorption isotherm at pH 6. (B) Linearized form of the Cr(VI) Langmuir adsorption isotherm.
Membranes 13 00740 g004
Figure 5. Variation of the distribution quotient of Cr(VI) (A) and the extraction percentage (B) as a function of the initial concentration of the metal.
Figure 5. Variation of the distribution quotient of Cr(VI) (A) and the extraction percentage (B) as a function of the initial concentration of the metal.
Membranes 13 00740 g005
Figure 6. Determination of the enrichment factor, E, at different initial Cr(VI) concentrations, considering all points (black dotted line, slope = 18.39), first 4 points at low concentrations (red line, slope = 33.30), and last four points at high concentrations (green line, slope = 15.09).
Figure 6. Determination of the enrichment factor, E, at different initial Cr(VI) concentrations, considering all points (black dotted line, slope = 18.39), first 4 points at low concentrations (red line, slope = 33.30), and last four points at high concentrations (green line, slope = 15.09).
Membranes 13 00740 g006
Figure 7. Determination of the stoichiometric ratio of the extracted Aliquat 336–Cr(VI) complex.
Figure 7. Determination of the stoichiometric ratio of the extracted Aliquat 336–Cr(VI) complex.
Membranes 13 00740 g007
Figure 8. Score plot of the PCA analysis. The points are colored according to the initial Cr(VI) concentrations in aqueous solution. (A) PC1 vs. PC2, (B) PC1 vs. PC2 vs. PC3.
Figure 8. Score plot of the PCA analysis. The points are colored according to the initial Cr(VI) concentrations in aqueous solution. (A) PC1 vs. PC2, (B) PC1 vs. PC2 vs. PC3.
Membranes 13 00740 g008
Figure 9. Forward iPLS results on the FTIR mean spectrum of the samples (black) showing the regions where minimum values of RMSECV were attained (red) and the optimal ones (green).
Figure 9. Forward iPLS results on the FTIR mean spectrum of the samples (black) showing the regions where minimum values of RMSECV were attained (red) and the optimal ones (green).
Membranes 13 00740 g009
Figure 10. Summary of the results of the final PLS regression model. (A,B) Outliers plot, (C) parity plot, (D) plot of the first two latent variables showing the Hotelling T2 ellipse.
Figure 10. Summary of the results of the final PLS regression model. (A,B) Outliers plot, (C) parity plot, (D) plot of the first two latent variables showing the Hotelling T2 ellipse.
Membranes 13 00740 g010
Table 1. Values of the parameters and the fitting results of the Langmuir model equation.
Table 1. Values of the parameters and the fitting results of the Langmuir model equation.
ParameterValue
qmax0.188 mmol/g
KL2199 cm3/mmol
Reduced Chi-Sqr0.02628
Adjusted R20.90717
ANOVA:
Regression sum of squares0.0247
Residual sum of squares0.0024
Regression mean square0.0123
Residual mean square0.00048
F-value25.6361
p-value0.0023
Table 2. Values of the RL parameter of the Langmuir model at different initial Cr(VI) concentrations.
Table 2. Values of the RL parameter of the Langmuir model at different initial Cr(VI) concentrations.
[ C r ( V I ) ] 0 ( mg / dm 3 ) RL
6.810.0064
14.530.0031
20.650.0021
27.100.0016
36.740.0012
40.730.0010
48.700.0009
Table 3. Percent of variance captured by the PLS Regression Model.
Table 3. Percent of variance captured by the PLS Regression Model.
X-Blocky-Block
ComponentThis ComponentTotalThis ComponentTotal
193.2393.2310.7110.71
26.0099.2414.7125.41
30.4899.7251.5576.96
40.2399.957.3584.31
50.0399.983.3087.62
60.0199.997.4995.10
70.0099.992.0797.17
Table 4. Analytical figures of merit for the final PLS regression model.
Table 4. Analytical figures of merit for the final PLS regression model.
FOM Results
AccuracyRMSEC3.73115
RMSECV6.82685
RMSEP3.3229
Bias−1.91847 × 10−13 (Cal)
0.185947 (CV)
R20.98145 (Cal)
0.940902 (CV)
Recovery%104.02 ± 4.12 * (Test)
sen 0.00001547 ppb
γ 3.8 ppb
γ−1 0.6 ppb−1
sel 0.0155
Linear range 5.8–100 ppb
LD 1.9 ppb
LQ 5.8 ppb
* 95% confidence level.
Table 5. Comparison of some representative works for Cr(VI) quantitation reported in the literature *.
Table 5. Comparison of some representative works for Cr(VI) quantitation reported in the literature *.
Detection MethodCarrier/ChromophoreLinear RangepHLDLQReference
Colorimetry2-hydroxy, 3-methoxy benzaldehyde thiosemicarbazone0.260–2.60 μg/cm360.014 μg/cm30.041 μg/cm3[59]
Colorimetry1,5-diphenylcarbazide0.03–3 μg/cm32.20.023 μg/cm30.076 μg/cm3[60]
Colorimetrydiazonium salt and citrazinic
acid
0.2–1.5 μg/cm3Alkaline medium0.04 μg/cm3 [61]
Rotational microfluidic paper-based device1,5-diphenylcarbazide0.5–10 μg/cm3Very acidic0.18 μg/cm3 [62]
Sol-gel monolithspyridine-functionalized sol-gel monoliths and diphenylcarbazide 1about 0.010 μg/cm3 [12]
Optodealiquat 336 and 1,5-diphenylcarbazide0.020–0.397 μg/cm330.011 μg/cm30.013 μg/cm3[11]
Optode1,5-diphenylcarbazide0.0024–1 μg/cm310.0007 μg/cm30.0024 μg/cm3[13]
Optodealiquat 336 and 1,5-diphenylcarbazide0.02–0.40 μg/cm330.0055 μg/cm30.0165 μg/cm3[14]
FTIR optodealiquat 3360.0058–0.1 μg/cm360.0019 μg/cm30.0058 μg/cm3This work
* Units of concentration have been standardized to facilitate the comparison of the different reported methods.
Table 6. Results of the analysis of the test samples by the final PLS regression model.
Table 6. Results of the analysis of the test samples by the final PLS regression model.
SampleMeasured (ppb)Predicted (ppb)Recovery%RMSEPp-Value
Test 11313.99107.66
Test 22023.57117.86
Test 33637.25103.47
Test 44042.32105.81
Test 56363.23100.37
Test 67375.85103.90
Test 78077.8997.48
Test 88689.07103.57
Test 99090.90101.01
Test 1010099.0099.00
Average 104.02
3.32290.0611
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Martínez de la Peña, A.; Rodríguez de San Miguel, E.; de Gyves, J. MID-FTIR-PLS Chemometric Analysis of Cr(VI) from Aqueous Solutions Using a Polymer Inclusion Membrane-Based Sensor. Membranes 2023, 13, 740. https://doi.org/10.3390/membranes13080740

AMA Style

Martínez de la Peña A, Rodríguez de San Miguel E, de Gyves J. MID-FTIR-PLS Chemometric Analysis of Cr(VI) from Aqueous Solutions Using a Polymer Inclusion Membrane-Based Sensor. Membranes. 2023; 13(8):740. https://doi.org/10.3390/membranes13080740

Chicago/Turabian Style

Martínez de la Peña, Armando, Eduardo Rodríguez de San Miguel, and Josefina de Gyves. 2023. "MID-FTIR-PLS Chemometric Analysis of Cr(VI) from Aqueous Solutions Using a Polymer Inclusion Membrane-Based Sensor" Membranes 13, no. 8: 740. https://doi.org/10.3390/membranes13080740

APA Style

Martínez de la Peña, A., Rodríguez de San Miguel, E., & de Gyves, J. (2023). MID-FTIR-PLS Chemometric Analysis of Cr(VI) from Aqueous Solutions Using a Polymer Inclusion Membrane-Based Sensor. Membranes, 13(8), 740. https://doi.org/10.3390/membranes13080740

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop