Practical Considerations on Indirect Calibration in Analytical Chemistry

Indirect or methodological calibration in chemical analysis is outlined. The establish- ment of calibration curves is introduced and discussed. Linear calibration is presented and considered in three scenarios commonly faced in chemical analysis: external cali- bration (EC) when there are no matrix effects in the sample analysis; standard addition calibration (SAC) when these effects are present and internal standard calibration (ISC) in cases of intrinsic variability of the analytical signal or possible losses of the analyte in stages prior to the measurement. In each kind of calibration, the uncertainty and confi- dence interval for the determined analyte concentration are given. evaluated from the calibration data.


Introduction
Direct absolute methods such as gravimetry, titrimetry or coulometry (among others) are directly traceable to SI units. Thus, traceability of contemporary instrumental methods is accomplished by applying indirect calibration procedures. In a direct calibration, the value of the standard (reference value) is expressed in the same quantity as the measurement of the equipment (for instance, the calibration of an analytical balance). In an indirect calibration, the value of the standard is expressed in a quantity different from the output one, that is, the measurement and the measurand are different. This is the most common kind of calibration in chemical analysis, for example, the calibration of a spectrophotometric method. Accordingly, the indirect calibration in analytical chemistry, also known as methodological calibration, is the operation that determines the functional relationship between measured values and analytical quantities, characterizing types of analytes, and their amounts. In this chapter, the establishment and validation of the mathematical model for the calibration function will be studied and discussed as well as the habitual scenarios concerning interferences coming from the chemical environment (matrix effects) and physical/instrumental lack of control leading to signal modification (standard additions and internal standard methodology). Confidence intervals for the calculated analyte concentration will be outlined and discussed.

The calibration in analytical chemistry
Calibration, as previously defined, can be assimilated to a mathematical function, where Y is the analytical signal or response corresponding to the analyte concentration x. The major analytical aim consists of finding this functionality. When applying absolute methods of analysis [1], where traceability is assured, such as gravimetry, titrimetry or coulometry, there is no need for indirect calibration. The analyte amount is evaluated from the analytical signal with the use of physicochemical constants (atomic mass, Faraday constant) and the concentration of the standardized titration solution in titrimetry, leading to a typical linear response model x ¼ KY: • Gravimetry: x ¼ Gðgravimetric factorÞ Â Yðmass of weighing formÞ • Titrimetry: x ¼ pðstoichiometryÞ Â Cðtitrant concentrationÞ Â Yðtitrant volumeÞ • Coulometry: x ¼ Yðtotal chargeÞ nðelectrons transferredÞFðFaraday constantÞ But in the field of relative methods (the majority of instrumental ones), traceability is reached just by performing an indirect calibration, that is by establishing the relationship between the analyte concentration and the analytical response. There are some theoretical relationships [2][3][4][5] verified for special analytical techniques as depicted in Table 1.
Nevertheless, in the common situations, the response function has to be empirically established by using standard analyte solutions. Many response functions exhibit linear zones, generally at low concentrations of analyte and other zones where a curvature appears, and in some cases,

Response function
Reference Analytical technique Beer's law Absorption spectroscopy Atomic emission spectrometry Wagenaar et al. [3] Atomic absorption spectrometry Andrews et al. [4] y ¼ AÀD there are regions where the response signal is independent of the analyte concentration [6].
Analysts are interested to the portion of the response function where the variation of the analytical signal with the analyte concentration contains useful analytical information. This portion of response function with analytical interest for calibration purposes is called the calibration curve. From the calibration curve, the amount of analyte in an unknown sample is evaluated from interpolation. The calibration step is of utmost importance within the realm of method validation.
In many situations, the calibration curve is linear, and a calibration straight line is obtained. From the mathematical models applied for establishing the response function, the most straightforward, studied, and easy to handle is the linear one. Accordingly, the linear calibration model will be considered throughout this chapter.
In case of non-linear response, there are several alternatives. The use of linearizing transformations is a common tool [7], but when this procedure does not work, curve-fitting methods are chosen. The best procedure is to try with polynomials of degrees successively larger until the F-test of residual variances indicates that the systematic error due to the lack of fit is negligible. If the plot has "N" points, the major degree polynomial to be used is of degree NÀ1. But the blind use of high-order polynomial may lead to overfitting. This kind of fitting is solved by multilinear regression [8]. This technique sometimes fails because the coefficient matrix is nearly singular. To avoid this, we can use orthogonal polynomials. The use of these polynomials leads to a diagonal coefficient matrix, overcoming singularities, and simplifying calculations. The orthogonal polynomials commonly used in curve fittings are the Chebischev's polynomials [9] and the Forsythe ones [10].
Aside from the advantages and applications of orthogonal polynomials, they are not at all the ultimate weapon. Rice proposed rational polynomial functions of the type FðxÞ ¼ that present a higher flexibility than orthogonal polynomials for adjusting purposes [11]. Another approach is to fit the points to a curve consisting of several linked sections of different geometrical shapes. This is the basis of the spline functions. Cubic spline [12] is the most used. They approximate the data to a series of cubic equations. These cubic links overlap in p interpolation points called "knots," and it is essential that splines show continuity at such points. This continuity applies to the spline function and its first derivatives. A total cubic spline has pÀ1 links, with four coefficients (S ¼ a þ bx þ cx 2 þ dx 3 ). Thus, 4(pÀ1) coefficients have to be calculated. This technique has been successfully applied in radioimmunoassay, gasliquid chromatography, and atomic absorption spectrometry [13].
The most usual technique for establishing a calibration straight line is the method of least squares. This consists of minimizing the function Q ¼ Three main requisites must be fulfilled before using this method [14], namely: • The x variable is free from error εðxÞ ¼ 0.
• The variance of response Y, σ 2 ðYÞ, remains uniform in the dynamic range of x (homoscedasticity).
In analytical calibrations, the analyte concentration is known with high accuracy and precision and, accordingly, the requirement (i) is accomplished. The condition (ii) is assumed by many researchers without a previous testing. There are several statistical assays for testing normality [13], and they should be performed before embarking in the fitting. Analysts have paid much more attention to the requirement (iii). In situations of heteroscedasticity (non-constant variance), the method of least squares can be applied but by using the so-called weighing factors [15], which are defined as w i ¼ 1 . Thus, the function to be minimized now is leading to expressions similar to the one obtained in simple linear regression. This is the weighted regression [13].
Let us assume that we deal with a situation often found in routine analysis where the three mentioned requirements are fulfilled. In the following, we consider the different scenarios we can face.

Metrological foundations on indirect calibration
Consider a new proposed analytical method which is applied to dissolved test portions of a given sample within the linear dynamic range of the linear analytical response (Y). This response may be expressed by the following linear relationship involving both analyte and matrix amounts [16]:Ŷ whereŶ is the estimated analytical response and A, B, C and D are constants.
A is a constant that does not change when the concentrations of the matrix, z, and/or the analyte, x, change. It is obviously related to the constant error blank correction. The blank must account for signals coming from reagents and solvents used in the assay as well as any bias resulting from interactions between the analyte and the sample's matrix. It is well known that the calibration blank and the reagent blank compensate for signals from reagents and solvents, but neither of them can correct for a bias resulting from an interaction between the analyte and the sample's matrix. The suitable blank must include both the sample's matrix and the analyte, and so it must be determined using the sample itself. The term A is called the true sample blank and can be estimated from the Youden sample plot, which is defined as the "sample response curve" [17]. Thus, by applying the selected analytical method to different test portions, namely m (a different mass taken from the test sample), different analytical responses Y are obtained. The plot of Y versus m is the Youden sample plot, and the intercept of the corresponding regression line is the so-called total Youden blank (TYB) which is the true sample blank [17][18][19]. However, when a "matrix without analyte" is available, the term A can be determined by evaluating the system blank (calibration and reagent blank).
Bx is the essential term that justifies the analytical method because it directly deals with the sensitivity to the presence of analyte.
Cz refers to the signal contribution from the matrix, depending only on its amount, z. When this term occurs, the matrix is called interferent. This contribution must be absent, because a validated analytical method should be selective enough with respect to the potential interferences appearing in the samples where the analyte is determined. Accordingly, the majority of validated methods do not suffer from such a direct matrix interference.
Dxz is an interaction analyte/matrix term. This matrix effect occurs when the sensitivity of the instrument to the analyte is dependent on the presence of the other species (the matrix) in the sample [20]. For the sake of determining analytes, this effect may be overcome by using the method of standard additions as we consider later.
Thus, the calibration function remains as: This function has to be established by using standards and could be applied to samples according to different methodologies. Calibration standards are prepared from primary standards containing the analyte or a surrogate, that is, a pure substance equivalent to analyte in chemical composition, separation and measuring that is taken as representative of the native analyte. It must be absent in the sample. Commonly, a surrogate is used in an internal methodology and in this case is termed as internal standard (IS) [21].
Three different scenarios can be considered for establishing the calibration function in order to determine the analyte in the sample: the external calibration (EC) (applicable when there is no matrix effect); the standard addition calibration (SAC) (used when matrix effect is present); and the internal calibration (IC) (applied for compensate uncontrolled analytical signal variations). These methodologies are outlined in the following section.

The external calibration
The external calibration (EC) is the most commonly used calibration methodology. It is named so because the calibration standards are not made up of the sample test portion. Instead, they are prepared and analysed separately from samples [21]. Accordingly, the signals recorded accounts for the analyte added as primary standard, reagents, solvents and other agents according to the analytical procedure, except the sample matrix. Accordingly, because EC is established in a free matrix environment, it can be applied for analyte determination only when sample matrix effects are absent. Thus, as a preliminary step within the method, validation to assess constant and proportional bias due to matrix effects has to be performed [22]. Being a matrix free calibration scenario, z ¼ 0, B is the slope of EC, b EC , and A can be taken as the system blank, a EC . In order to evaluate the goodness of the fit, the regression analysis of the analytical signal on the analyte concentrations established in the calibration set yields the calibration curve for the predicted responses. The simplest model is the linear one, very often found in analytical methodology, leading to predicted responses according tô Eq. (3) must be checked for goodness of fit.
The correlation coefficient r ¼ , although commonly used, especially in linear models, is not appropriate owing to the little value of this parameter for detecting curvature [23,24]. In statistical theory, correlation is the measure of the association between two random variables, but in our case, x and Y are strongly related. Thus, there is no correlation in its mathematical sense. Values of r near þ1 or À1 provide an environment of respectability but not much else. Some authors apply statistical tests for significance of the correlation coefficient, for instance, the student t-test [13] , but they cannot ward off danger because the null hypothesis is that the variables are uncorrelated (zero correlation), and accordingly, a small r value can be considered significantly different from r ¼ 0. As Thompson [23] pointed out, "certainly it is true that, if the calibration points are tightly clustered around a straight line, the experimental value of r will be close to unity. But the converse is not true". Thus, some more suitable criteria should be considered. A very simple way to prove that the linear model suitably fits the experimental data and is right for searching possible calibration pathologies is the analysis of residuals [8,13,25]. So, if the model is suitable, the residuals should be normally distributed. This can be assessed by plotting them on a normal probability graph. The presence of curvature reveals a lack of fit due to non-linear behaviour. A residual segmented pattern may indicate heteroscedasticity in the data, and a weighted linear regression could be used.
Another parameter measuring the goodness of the fit is called the on-line linearity [26] and is a parameter that measures the dispersion of the points around the calibration straight line and is evaluated as the relative standard deviation of its slope: on-line linearity Nevertheless, the best way to test the goodness of fit is by comparing the variance of the lack of fit against the pure error variance [27]. For an adequate assess of the lack of fit of the linear model, a suitable experimental design for performing the calibration is needed as indicated in the following [28]: i. At least six calibration points spaced over the concentration range of the method scope are required for establishing the calibration straight line. ii.
Calibration standards should be measured over 5 days for suitably covering the possible sources of uncertainty.
iii. Each calibration standard should be measured in triplicate to account for pure error variance.
From these data, we can test homoscedasticity. Accordingly, we have a triplicate of responses for each calibration standard and hence an estimation of the pure error variance of the response is available at each calibration point. We can apply the Cochran's assay because the number of observations is same for all concentration levels of analyte. Thus, if the number of calibration standards is N and they are replicated n times, the Cochran statistics is calculated as: where s 2 i is the response variance at the concentration level i and s 2 max is the maximum variance. This value is compared against the critical tabulated value C tab ðN, n, PÞ, P being the selected confidence level. If C ≤ C tab , then the response variances can be considered to be uniform across the range of analyte concentrations, and an estimated pooled sum of squares due to pure errors, SS PE , can be obtained: The residual sum of squares of the model SS R is given by where Y ij is the recorded analytical signal of the calibration point i at the replication j andŶ ij . This value can be split into two terms: the sum of squares corresponding to pure error (SS PE ) and the sum of squares corresponding to the lack of fit (SS LOF ): The pure error variance is SS PE /(nÀ1), and the variance of the lack of fit, by considering NÀ2 degrees of freedom for SS R , is SS LOF /(NÀnÀ1). So, for assessing the adequacy of the model, the Fisher F-test is applied: The calibration model is considered suitable if less than the one-tailed tabulated value F tab ðN À n À 1, n À 1, PÞ exists at a P given confidence level.
Once the model is adequate for application, analyte determination is carried out by interpolating the analytical signal of the sample in the calibration model. Typical statistical calculations for evaluating the variances of slope, intercept, its covariance as well as the uncertainty associated to the estimated analyte concentration can be found in several texts, for instance, Miller and Miller [13]. Thus, if Y 0 is the response signal recorded by applying the analytical method on the sample, the concentration of native analyte, x 0 , is given byx In order to evaluate its standard deviation, and the corresponding expanded uncertainty, the theorem of variance propagation is applied. The propagation of variance is the common approach for evaluating the uncertainty of indirect measurements according to the current edition of the guide for the expression of uncertainty measurement (GUM). However, an essential limitation has to be taken into account. The non-linearity of the function (here the calibration function) must be negligible. This is fundamental because the function is expanded in a Taylor series, and then, it is truncated by neglecting second-and higherorder terms. To avoid this drawback, the propagation of distributions instead of the propagation of variance is a very suitable way for estimating the measurement uncertainty. The application of Monte-Carlo method to carry out the propagation of distributions is very effective [29].
Saying that brute-force Monte-Carlo (MC) methods are "very effective" may seem strange to some readers, as one major problem of MC is their methodological in-efficiency. It is due to large sampling variance of the relatively small samples acceptable in computationally demanding applications. In other words, any acceptable sample of 100 values may have a large random unknown error, generally different from any other sample of comparable size. To overcome this inefficiency, approximate simplified surrogate models are often used to allow for sampling a much as 10 6 times, just to reduce sampling variability. I would thus rather call MC methods 'general', 'useful', 'simple' and 'powerful' etc., as they apply to any parametric model and any distribution (if a random generator can be found), and can be utilized by anybody with very little statistical training.
But in our case, where the calibration function has been considered linear, the use of theorem of variance propagation can be applied without risks: Considering the following equivalences (see Ref. [13] for instance): After some algebraical manipulations, we get If the signal Y 0 is obtained as the average of m measurements, we have And the corresponding expanded uncertainty can be evaluated from the tabulated Student tstatistics or by assuming a Gaussian distribution and using the z score at a given confidence level (generally P ¼ 95%) and so: EC is adequate for analytical procedures that could be considered as methods free from matrix effects, but it has the main limitation coming from the assumption that the different environments (matrices) of the calibration standards (solvent, buffer,…) and of the samples are equivalent, and they have no effect on the calibration function [21]. If this assumption is incorrect, additive and/or proportional systematic errors may appear. Accordingly, in a preliminary stage within the method validation, constant and proportional bias due to matrix effects must be investigated with the help of standard addition calibration and Youden plot [22].

Standard addition calibration
The standard addition calibration (SAC) or standard addition method was originally proposed in 1937 by Hans Hohn in polarographic studies [30]. He used this strategy in order to avoid the matrix effects on the intensity of emission signal, and nowadays, it is widely used in chemical analysis. SAC can be applied with three fundamental goals [31]: • To determine analytes in samples where the analyte-matrix interactions lead to inaccurate results when the EC is used.
• To determine analytes where the content in the sample is smaller than the quantitation limit but within the range of analytical sensitivity.
• To check the accuracy of an analytical result when no reference materials or reference method is available (recovery assay).
In essence, the calibration for the two first purposes comprises three steps [32]: a. Measure the analytical response produced by the test solution.
b. Spike the test solution with one or more amounts of analyte to get corresponding solutions and measure the new responses.
c. From the responses, calculate a straight-line fit of the experimental data and from that evaluate the concentration that produced the response obtained from the untreated test solution.
The SAC can be performed either at a final fixed volume or at a variable volume [19]. In this discussion, we only consider the first case by working at constant final volume.
Consider now the application of the analytical procedure to a dissolved test portion of an unknown sample within the linear working range. The analyte concentration x is the sum of the fixed native concentration coming from the sample (volume of test portion V 0 ) and the variable spiked concentration (spiked volume, V spike ) and keeping a final constant volume V.
The amount of matrix in the test portion (z) is constant. Accordingly, the analytical response can be now modelled as: where C native is the actual concentration of the analyte in the unspiked sample, C spike the actual concentration of the spiked analyte and a SAC and b SAC are the intercept and the slope of the SAC calibration straight line. If we try to estimate the analyte concentration of a spiked sample by using the external calibration line, we obtain an estimation of the total observed analyte concentration: For the unspiked sample, C spike ¼ 0, we obtain According to Eqs. (16) and (17), the spiked concentration of the analyte is estimated from the external calibration as: From Eq. (18), an overall estimation of the overall consensus recovery is calculated as: When proportional bias is absent, we have b SAC ¼ b EC , and that implies Rec ¼ 1. This must be tested for statistical significance by using the student t-test [22]: with the recovery standard deviation given by: Thus, if the degrees of freedom ν corresponding to the uncertainty of consensus recovery are known, student t-statistic is compared with the critical two-tailed tabulated value, t tab (ν,P), at P% confidence. If t ⩽ t tab , the consensus recovery is not significantly different from 1. Alternatively, instead of t tab , a coverage factor k taken as z score may be used for the comparison. Typical values are k ¼ 2 or k ¼ 3 for 95 or 99% confidence, respectively [22], so • if jRec À1j sRec ≤ k, the recovery is not significantly different from 1.
• if jRec À1j sRec > k, the recovery is significantly different from 1, and the results have to be corrected by Rec.
Although recovery is sometimes considered a separate validation parameter, it should be established as a part of method validation because it is directly related to the trueness assessment [33]. Aside from the statistical testing considered above, the Association of Official Analytical Chemists (AOAC) has published tables of acceptable recovery percentages as a function of the level of analyte in the sample (see Table 1 of [22]). The relative uncertainty for proportional bias owing to matrix effects is taken as s Rec Rec according to SAC.
The relationships between the analytical signal and the analyte concentration when a matrix effect is present are given by Eq. (15). The independent term "A" is the total Youden blank, which is included in the intercept of the SAC calibration (a SAC ¼ A þ b SAC C native ). The Youden's plot [17][18][19] consists of plotting the analytical response (Y) against the amount of the test portion taken for analysis: The intercept of this plot is an evaluation of the TYB, which is the sum of the system blank (SB) corresponding to the intercept of the EC (a EC ) and the YB associated with the constant bias in the method [13]. Thus, we can equate TYB ¼ A, SB ¼ a EC and YB ¼ Aa EC . We can define the method constant bias as: The uncertainty of the constant bias can be obtained by the law of variance propagation [22]: The variances s 2 a EC , s 2 b EC and the covariance are obtained from the statistical parameters of the EC straight line and s 2 A from the Youden's plot. Once s θ is calculated, the constant bias may be assessed for significance as in the case of recovery.
• If jθj sθ ≤ k, the constant bias is not significantly different from 0.
• If jθj sθ > k, the constant bias is significantly different from 0, and the results have to be corrected by θ.
Accordingly, if after performing the assessment of proportional and constant bias matrix effects are present, the uncorrected result x 0 , found by EC, must be suitably corrected as Another way of getting the correct result from the reading of analytical signal Y 0 is On the other hand, when using the SAC for evaluating the analyte concentration x 0 of a sample, its standard deviation can be obtained by applying the theorem of variance propagation to the function leading to After several algebraical manipulations, we obtain But many workers apply the SAC without considering the true blank, that is, by setting A ¼ 0 and s A ¼ 0 leading to This expression is presented in several standard analytical textbooks, for instance [13,19,32,34]. However, Ortiz et al. [35] pointed out that when extrapolating, the analyte concentration is obtained by setting Y 0 ¼ 0 and calculating x 0 ¼ Àa SAC =b SAC , but even in this case, the uncertainty of the signal must be included in calculations, leading to The SAC, as it has been outlined, is considered as an extrapolation method but an interpolation approach is available [32,36]. A plot of the data obtained from SAC and how the analyte concentration is predicted by extrapolation are depicted in Figure 1. Nevertheless, an interpolation alternative is also gathered there. The latter is discussed in the following.
What value of the analytical signal Y 0 will correspond to a spiked x value that is equal to the concentration of the native analyte? That is: And if we disregard the true blank, we get Thus, the native analyte concentration can be obtained by interpolation by setting the analytical signal for the sample as the double of the signal, corresponding to the unspiked sample minus the true blank: This leads to a variance for the native analyte According to Andrade et al. [32,36], the use of extrapolation in the SAC is a risky practice because it may lead to biased prediction and uncertainties substantially different from interpolation. Confidence interval from extrapolation is always higher than those obtained by interpolation.

The internal standard calibration
The method of internal standard calibration (ISC) was first applied in the 1950s in several analytical fields [37][38][39][40][41]. This method is especially useful when the analytical response varies slightly from run to run due to different causes, for instance: • Temperature fluctuations in atomic emission spectrometry.
• Changes in the capillary characteristics in polarography.
• Inhomogeneities in the effective magnetic field due to shielding effect in nuclear magnetic resonance (NMR).
• Variability in the injection volume in gas chromatography (manual injection).
• Irreproducibility of automatic injectors in capillary electrophoresis.
• Differences in the nature of particulate matter in the sample in X-ray fluorescence.
The use of an internal standard is also needed for analytical methods where there are multiple sample preparation steps, especially when volumetric recovery at each step may vary (extraction with separation cartridges) or when involving chemical derivatizations with low or variable yields of reaction.
An internal standard is a substance different from the analyte but that has physicochemical properties very similar to the analyte. Evidently, the internal standard cannot be a component of the sample.
It is added to the sample, and the patterns in known amounts and the signal produced by both the analyte and the internal standard are measured. If in repeated measurements, there is signal oscillation, it will occur both in the analyte and in the internal standard, and the ratio of the signals of both will not change.
Thus, instead of the response Y, the ratio of responses Y/Y IS is used in the calibration procedure. Assuming that in the instrumental method the signal is in direct proportion to the analyte and internal standard, we get: Here, Y is the analytical signal due to the analyte and Y IS is the analytical signal corresponding to the internal standard. The calibration straight line is performed as in EC by preparing standards at several analyte concentrations and with the same concentration fixed for internal standard x IS . Thus, the calibration constant F is evaluated. Whereas the dispersion of the calibration straight line Y ¼ kx may be significant, the one obtained with the ISC is negligible.
The sample is then treated in the same way by spiking the internal standard at the same concentration in the standards. Thus, if the reading of the sample is By applying the variance propagation law and considering negligible variance of x IS , we get The main advantage of ISC is that this quantification method does not need a previous calibration because it is implicit in the quantification [21]. Accordingly, the use of onepoint calibration method can be used. It only requires the addition of known and equal amounts of internal standards to the standard analyte solution and to sample solution and measures the analytical signals of analyte and internal standard in the standard and in the sample. Evidently, the signals of analyte and internal standard must be distinguishable without overlapping. Thus, Another exclusive feature of ISC is the possibility of performing the quantification of several analytes of the same chemical family in the same test portion and in a unique internal calibration with a single internal standard. Consequently, it could be possible to evaluate the mass fraction of each analyte according to [21].
Accordingly, ISC is a very powerful method for congener analysis (for instance in fat analysis, determination of waxes, sterols, aliphatic alcohol and so on) by using only a unique internal standard.

Synthesis
Indirect calibration is a key concept for method validation. Instrumental analysis involving indirect calibration is a common feature in routine analysis, and three typical scenarios can be found depending on the analyte-matrix interaction and the uncontrolled variation of the analytical signal owing to intrinsic characteristics of the analytical process. Thus, when the interaction of the matrix of sample is negligible, the external calibration is the normal choice. Otherwise, the Standard Addition calibration together with the Youden plot have to be applied. In cases where there are non-random signal variations run to run or possible analyte losses due to sample preparation procedures or derivatization reactions, Internal Standard calibration must be considered. These three approaches have been outlined and discussed. Uncertainty values for the analyte concentration coming from the calibration step are considered and evaluated from the calibration data.

Antonio Gustavo González
Address all correspondence to: agonzale@us.es Department of Analytical Chemistry, Faculty of Chemistry, University of Seville, Seville, Spain