Next Article in Journal
Effective In Vitro Control of Two Phytopathogens of Agricultural Interest Using Cell-Free Extracts of Pseudomonas fluorescens and Chitosan
Next Article in Special Issue
Gas Chromatographic Fingerprint Analysis for the Comparison of Seized Cannabis Samples
Previous Article in Journal
Design and Synthesis of Novel Betulin Derivatives Containing Thio-/Semicarbazone Moieties as Apoptotic Inducers through Mitochindria-Related Pathways
Previous Article in Special Issue
Pure Ion Chromatograms Combined with Advanced Machine Learning Methods Improve Accuracy of Discriminant Models in LC–MS-Based Untargeted Metabolomics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Chromatographic Applications in the Multi-Way Calibration Field

by
Fabricio A. Chiappini
1,2,
Mirta R. Alcaraz
1,2,
Graciela M. Escandar
2,3,
Héctor C. Goicoechea
1,2 and
Alejandro C. Olivieri
2,3,*
1
Laboratorio de Desarrollo Analítico y Quimiometría (LADAQ), Cátedra de Química Analítica I, Facultad de Bioquímica y Ciencias Biológicas, Universidad Nacional del Litoral, Ciudad Universitaria, Santa Fe S3000ZAA, Argentina
2
Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Godoy Cruz C1425FQB, Argentina
3
Departamento de Química Analítica, Facultad de Ciencias Bioquímicas y Farmacéuticas, Universidad Nacional de Rosario, Instituto de Química de Rosario (IQUIR-CONICET), Suipacha 531, Rosario S2002LRK, Argentina
*
Author to whom correspondence should be addressed.
Molecules 2021, 26(21), 6357; https://doi.org/10.3390/molecules26216357
Submission received: 27 September 2021 / Revised: 11 October 2021 / Accepted: 18 October 2021 / Published: 21 October 2021
(This article belongs to the Special Issue Chromatography and Chemometrics 2021)

Abstract

:
In this review, recent advances and applications using multi-way calibration protocols based on the processing of multi-dimensional chromatographic data are discussed. We first describe the various modes in which multi-way chromatographic data sets can be generated, including some important characteristics that should be taken into account for the selection of an adequate data processing model. We then discuss the different manners in which the collected instrumental data can be arranged, and the most usually applied models and algorithms for the decomposition of the data arrays. The latter activity leads to the estimation of surrogate variables (scores), useful for analyte quantitation in the presence of uncalibrated interferences, achieving the second-order advantage. Recent experimental reports based on multi-way liquid and gas chromatographic data are then reviewed. Finally, analytical figures of merit that should always accompany quantitative calibration reports are described.

1. Introduction

In the last decades, the use of chromatographic techniques coupled to multidimensional detection systems has gained the attention of analytical chemists. This is due to the fact that its combination with multi-way calibration tools allows one to obtain a large amount of chemical information and resolve highly complex systems outperforming the advantages of classical univariate chromatographic methods [1,2,3]. It has been demonstrated that the use of multi-way chromatographic data analysis yields better analytical performance than those based on univariate calibration, achieving additional analytical benefits, such as reducing the time of analysis, decreasing the solvent consumption and avoiding sample pre-processing steps. All these characteristics are consequences of the fact that if properly processed, multi-way data allow one to achieve the well-known second- and third-order advantages [4].
In the context of higher-order calibration methods (second-, third- and higher-order), the concept of the second-order advantage relies on the fact that, in certain circumstances, the contribution of individual sample constituents can be accurately obtained, even in the presence of unmodeled or unexpected interferences [5]. This property has been extensively demonstrated in a wide variety of applications, as can be verified in a vast number of publications. Even though this property was first proposed for second-order/three-way calibration methods, it has been extended to higher-order calibration methodologies. However, in the case of third- or higher-order data (four-way calibration and beyond), the additional advantages are still in discussion, and there is no general consensus about the real nature of the third-order advantage. Nevertheless, experimental and theoretical work are continuously growing in order to investigate the benefits associated to the increment in the number of instrumental data modes [6].
Second-order data consists in a collection of a bidimensional data array for a given sample. The bidimensional signals acquired for a set of samples can then be joined to obtain a three-way data array. This tensorial object is characterized by three experimental modes, one corresponding to the number of samples and the additional ones representing the measured analytical information. In particular, second-order data are generated when the analytical signals are recorded in two independent instrumental modes. In this regard, chromatography with spectral detection represents the most reported analytical application in the second-order calibration field. For instance, liquid chromatography with spectroscopic diode array detection (LC-DAD), fast scanning fluorescence detection (LC-FSFD), spectral mass spectrometric detection (LC–MS) and gas chromatography with spectral mass spectrometric detection (GC–MS) allow one to obtain second-order data.
On the other hand, third-order data are obtained when an additional experimental or instrumental mode is incorporated. In principle, there is no theoretical limitation to the benefits that may be brought about by additional instrumental modes in higher-order calibration [7]. Nevertheless, there might be a limitation for the generation of multidimensional data, since most multi-way chemometric modelling approaches imply strong assumptions about the mathematical structure of the data arrays. This restraint may be apparently overcome by virtue of the development of novel instrumental technology.
Even though countless alternatives would be possible for implementing third-order data acquisition, third-order chromatographic methodologies are one of the most challenging approaches. In the first attempt of acquiring third-order chromatographic data, in 1981, Apellof and Davidson reported a chromatographic method with excitation-emission matrices (EEM) detection for qualitative analysis [8]. Since then, interest in this type of data has been growing, accompanied by the development of novel and robust chemometric models capable of exploiting the potentiality of the multi-way data arrays [6], and also by the progress and expansion of analytical instrumentation.
Notwithstanding the wide variety of chemometric models covering a broad range of possibilities, it should be noted that there is no single chemometric model able to fit all the types of data that can be experimentally measured in a laboratory. Among all the reported models, the most used ones in second- and third-order calibration are parallel factor analysis (PARAFAC) [9], multivariate curve resolution-alternating least-squares (MCR-ALS) [10], and partial least-squares-based techniques, followed by residual multilinearization, e.g., unfolded- or multiway- partial least-squares, followed by residual multilinearization (U-PLS/RML and N-PLS/RML) [11]. In addition, variants of the above-mentioned models have been developed, aiming to overcome some limitations and to improve its performance. PARAFAC2 [12], augmented PARAFAC (APARAFAC) [13] and the family of alternating multilinear decomposition models (AMLD) [7] are examples of PARAFAC variants.
Together with data modelling, another essential part in the development and validation of calibration methodologies concerns the estimation of analytical figures of merit (AFOMs). These figures are numerical parameters used to characterize the performance of a developed protocol and to compare the relative success among different methodologies [14]. In multi-way calibration, the theory of error propagation is the core of the AFOM estimators. This topic has been a focus of recent discussions, and important advances have been reported in the literature [14] regarding the estimation of crucial AFOMs, such as sensitivity, analytical sensitivity, limit of detection and limit of quantitation.
The present review is intended to provide a comprehensive coverage of the most relevant alternatives reported for the generation and analysis of second- and third-order chromatographic data and their application in the multi-way calibration field.

2. LC Multi-Way Data Generation

2.1. Second-Order Data

Second-order data are characterized by the presence of two different instrumental modes in the data array acquired for each experimental sample. From the experimental point of view, there are two conceptually different methods to generate second-order data: (1) using a single instrument and (2) connecting two instruments in tandem. In the former case, measurements are directly performed on a single equipment, either because two components of the instrument itself are able to provide each of the data modes, or because first-order (vectorial) measurements are made as a function of time, and these vectors are joined to produce a data matrix per sample. In the latter, each of the connected instruments provides a data mode or matrix direction to the measured matrix data.
Although not directly related to chromatography, EEM are, probably, the most explored second-order data measured in a single instrument [15], and are good examples to introduce the notion of trilinearity. EEM data are, under certain circumstances, classified as trilinear [16] and, thus, robust and often unique trilinear decomposition models can be applied to them. This is because the excitation and emission component profiles are independent phenomena, and (properly normalized) do not depend on the sample. However, second-order fluorescence data may not always be trilinear [16]. For example, if inner filter effects occur in one of the data modes, lack of trilinearity will be observed, because the fluorescence profile of a given component will be different across samples in the mode where the inner filter takes place [17]. This is due to the fact that the magnitude of the inner filter depends on the concentration of the constituent producing the effect. These data will be classified as non-trilinear type 1 if the inner filter affects a single data mode, and type 2 if it affects both instrumental modes [16].
Briefly, non-trilinear (and non-multilinear) data types are those where the multi-linearity is lost because either (1) component profiles are not constant along one or more modes or (2) there is mutual dependence between the phenomena taking place in the instrumental modes. Specifically, three-wat data are non-trilinear type 1 if there is a single trilinearity breaking mode (non-constant profiles along this mode), non-trilinear type 2 if there are two breaking modes and non-trilinear type 3 if there is mutual dependence between the two instrumental modes.
The connection of two instruments in tandem provides another convenient way of generating second-order data. In fact, chromatography with multivariate spectral detection shares with matrix fluorescence data the priority in the publication record regarding second-order calibration protocols. Detection using a diode array detector (DAD), a fast-scanning fluorescence detector (FLD) or a spectral mass spectrometer (MS) provides the spectral mode to a liquid chromatograph, which itself is responsible for the elution time mode of the measured data. In the case of gas chromatography, spectral MS detection is the method of choice for generating second-order data [1]. In a few cases, electrophoretic measurements have replaced the chromatographic separation mode [18].
All data stemming from chromatographic measurements should be considered, in principle, non-multilinear type 1 [2]. This is because the reproducibility across sample injections, both in the position and shape of chromatographic peaks, is never perfect. However, if the elution time mode is indeed reproducible or quasi-reproducible, either because the total experimental time of each run is short enough, or because the instrument itself provides reproducible data, e.g., a gas chromatograph, then the data could be considered trilinear [19].

2.2. Third-Order Data

As an extension of second-order, third-order data are characterized by the existence of three instrumental modes in the data that are collected for a given sample. In the chromatography field, third-order data are usually obtained through the hyphenation of two different instrument, i.e., a chromatograph coupled to a second-order detector (for example, EEM), or by means of a two-dimensional chromatographic instrument with vectorial signal detector, e.g., DAD or FLD, among others [20].
Literature reports regarding third-order data analysis are still scarce in comparison with second-order calibration methods for chromatographic applications. This fact could be a consequence of the complexity associated to the required instrumental arrangements and the intrinsic difficulties of monitoring in-flow analysis.
In time-dependent experiments, the sample composition changes with the time, which is monitored when registering the analytical signal. In these kinds of systems, measurements can be conducted by performing different experiments, either in steady-state conditions or in continuous-flow conditions. In steady-state flow applications, the system is monitored at different periods of time in a condition where it no longer evolves, and any property associated with the flow or time remains constant. Stopped-flow systems and fraction collection-based methodologies are examples of steady-state flow applications. Under these circumstances, the instrumental modes of the third-order data are fully independent between them. On the contrary, in continuous-flow systems, the first issue to be considered is the synchronization between the evolving rate of the system and the scan rate of the detector. If full synchronization exists, the instrumental modes are mutually independent and the data will fulfil one of the multilinearity conditions [16]. For instance, DAD systems coupled to LC allow acquiring an entire spectrum at every chromatographic time and then, bilinear second-order data (elution time × spectra) are obtained. However, for detectors based on spectral scanning, the synchronization issue is not trivial.
The most reported strategy for generating third-order data consists in the acquisition of EEM as a function of the chromatographic time (LC-EEM). To generate these kinds of data, several instrumental arrangements have been proposed. In this regard, it should be reminded that a fluorescence spectrometer is a second-order instrument that enables the acquisition of bidimensional arrays as a function of elution time. Nevertheless, the commercially available spectrofluorometers operate through mechanical motion of the gratings to generate complete spectra or EEMs. These motions demand a finite time, which is considerably larger than the one required to benchmark the evolving rate of the chromatographic system. Despite the fact that modern analytical instrumentation has simplified the generation of multi-way data, third-order chromatographic acquisition is still a challenging task from the instrumental standpoint, and also constitutes a challenge from the chemometric perspective.
A number of analytical methodologies based on the generation of third-order LC data using fluorescence detection have been reported. One of the strategies follows the path pioneered by Apellof and Davidson [8] who proposed the generation of third-order LC data with qualitative aims by acquiring EEM at discrete chromatographic times. This strategy is based on the collection of discrete fractions eluting from the chromatograph for which an EEM is then obtained using a conventional spectrofluorometer. In 1997, R. Bro [21] proposed, for the first time, an approach for the generation of third-order LC-EEM data with quantitative aims based on Apellof and Davidson’s idea. Samples were analysed under identical conditions and a four-way array was then built and subjected to chemometric decomposition. In 2014, Alcaraz et al. [22] implemented the same strategy for the quantitation of three fluoroquinolones in water samples using a custom-made fraction collector, which was connected to the end of the chromatographic column and enabled the collection of fractions in a 96 wells-ELISA plate. At the end of the sampling, the plate was placed into a conventional spectrofluorometer for registering the EEM of each well. Even though the instrumental modes are mutually independent and synchronization between rates is not demanded by this strategy, several issues are undesired from the chemometric standpoint. Despite the fact that the third-order data obtained for each sample are trilinear, the four-way array built for a set of samples does not behave as quadrilinear, but as non-quadrilinear type 1, because of the lack of time and peak shape reproducibility between runs. Recall that non-quadrilinear four-way data are type 1 for a single breaking mode, type 2 for two breaking modes, type 3 for three breaking modes and type 4 if there is dependence among pairs of instrumental modes.
A different chromatographic approach was first introduced by Muñoz de la Peña’s group [23] and was then implemented in other applications with quantitative purposes [24,25]. In this case, to avoid stopping the flow, the FLD capabilities were exploited and the third-order LC-EEM data were measured by performing several chromatographic runs of aliquots of the same sample. The corresponding second-order LC-FLD data matrix were then registered, changing the excitation wavelength for every injected aliquot. Hence, the three-way array was built by joining the data matrices acquired at each excitation wavelength. This approach is characterized by the fact that a single chromatograph is used for the acquisition of the third-order data, including an autosampler and an FLD. However, deviations of multilinearity are present when this approach is used. The first aspect to be considered is the fact that time shifts, and peak distortions may appear among runs, breaking the trilinearity of the third-order LC-EEM data (the same phenomenon occurring in second-order calibration). Hence, the four-way data array will not be quadrilinear, but non-quadrilinear type 4, which is a complex scenario for chemometric models. To the best of our knowledge, there are no efficient pre-processing tool for recovering the trilinearity of LC-EEM third-order data, and no adequate models dealing with non-quadrilinearity type 4. However, this data can be properly modelled by implementing a bilinear decomposition of a super-augmented bilinear data matrix (see Section 3). It should be noticed that only a small number of runs are performed per sample to reduce time, sample and reagent consumption, leading to imbalanced data arrays with many data points in the emission and the chromatographic directions and only a few points in the excitation direction. These issues may hinder its application for the analysis of complex systems with a large number of constituents.
Finally, the most explored and promising approach is the one consisting in the hyphenation of a chromatograph and a fast-scanning spectrofluorometer connected through a fluorescence flow-cell. The first work that report the implementation of an online EEM registering system was done by Goicoechea’s group [25], who described and analysed the advantages and disadvantages of this alternative in comparison with the two aforementioned ones. This approach presents the great advantage of simultaneously recording the EEM on-line with the LC procedure, achieving a drastic reduction of time analysis, reagents and sample. Notwithstanding, the main drawback is related to the strong dependence of the elution time mode with both spectral modes, which leads to a loss of trilinearity in the third-order data, and to non-quadrilinear data of type 4 for multi-sample analysis. This phenomenon occurs as a consequence of the lag between the elution and the fluorescence scanning rates. To cope with these limitations, researchers have implemented novel instrumental configurations [26], developed new spectrometers [27] and introduced new chemometric alternatives [28]. Escandar’s group implemented a chromatographic setup that decreased the time dependence effect by reducing the linear flow rate of the mobile phase, incorporating a large inner-diameter tube between the column and the flow-cell [26]. In this way, the authors reported that the time-dependence effect is negligible, and the third-order data of individual samples are indeed trilinear. In addition, due to the large chromatographic times and the slow chromatographic rate, reproducibility in the elution time among samples was observed leading to quadrilinear data. More recently, Alcaraz et al. have presented an ultra-fast multi-way detector that enables measuring a complete EEM by bidimensional excitation and emission spatial dispersion [27]. This device, based on the use of a CCD camera, allows acquiring fluorescence images in the order of milliseconds. The first advantage of this setup is the feasibility of acquiring trilinear third-order data with several data points in all the three modes. For instance, three-dimensional arrays of size 1450 × 240 × 320 for elution time, excitation and emission modes, respectively, were reported. In this case, non-quadrilinear type 1 four-way data arrays are generated, which can be easily decomposed by known chemometric models such as APARAFAC, MCR-ALS and U-PLS with residual quadrilinearization (RQL), among others. This device represents a step-forward in the field of third-order data acquisition for dynamic systems.
Totally different third-order chromatographic data can be generated by bidimensional (2D) chromatography in combination with vectorial detection or by three-dimensional (3D) chromatography. In the case of 2D chromatography, a sample is driven through two independent columns; the effluent from the first-dimension column is sequentially injected into the second-dimension column. At the end of the latter, a spectral detector registers a vectorial signal at each chromatographic time. In this way, the generated third-order data involve the first-dimension chromatographic time, the second-dimension injections and the corresponding spectra. In this regard, LC2-DAD [29] and GC2-MSTOF (time of flight) methodologies have been reported and four-way calibration method have been successfully implemented [30]. More recently, a new methodology based on 3D GC with univariate detection (FID or single quadrupole MS) was proposed as an alternative to yield third-order data with quantitative aims [31]. In this case, the three-dimensional data array is built by the combination of the three chromatographic modes. All these strategies share the same undesired particularity in chemometric terms, i.e., the elution time shifts and peak distortions that occur among samples in both the first- and the second-dimension (and the third one in the case of 3D chromatography), which, in principle, preclude the application of a multilinear model.

3. LC Multi-Way Data Analysis: Chemometric Models and Algorithms

As shown in the previous section, the generation of three- and four-way data implies the construction of a variety of data arrays, characterized by various mathematical properties. Hence, it is not surprising that a large diversity of models and algorithms have proliferated in recent decades.
The concept of data multilinearity, i.e., bilinearity, trilinearity and quadrilinearity is the common thread in method taxonomy and characterization [16]. These concepts should orientate the analyst in the selection of the most suitable method for a given calibration scenario. In this regard, it is important to consider that despite the structure of the raw experimental data, they can be subjected to mathematical operations prior to chemometric processing, in order to fulfil the conditions required for a successful chemometric decomposition. For instance, second-order matrices can be arranged as follows: (1) they can be stacked in a third mode to give rise to a three-way data array, (2) joined in-plane to produce an augmented data matrix in either of the instrumental directions, or (3) unfolded into vectors and then join the vectors to produce a single matrix. Two key factors determine the selection of any of the latter structures: the existence of phenomena producing constituent profiles which vary from sample to sample, and/or mutually dependent phenomena occurring in the two instrumental modes. The different data arrays that can be obtained from a second-order dataset, prior to chemometric modelling, are illustrated in Figure 1.
Following the same criteria described before, the third-order arrays in four-way calibration can be organized as follows: (1) they can be disposed in a fourth mode yielding a four-way data array; (2) organized in an augmented three-dimensional data array in either of the instrumental directions; (3) joined in-plane to generate a bidimensional data matrix with one or two augmented modes; and (4) completely unfolded into vectors and then stacked in a single matrix. These possibilities are graphically summarized in Figure 2.
In general terms, three families of chemometric modelling approaches can be distinguished for three- and four-way data analysis. Each group of methods is characterized by assuming different hypotheses about the data structure and exhibit different degrees of flexibility. This leads to distinctive model interpretabilities and capabilities of exploiting the benefits of multi-way data arrays in terms of analytical performance and solution uniqueness. Although a vast list of mathematical models and algorithms exists, only the most relevant ones are summarized and detailed in this report.
In accordance with the different types of required mathematical operations, stacked multi-way data that conform the property of multilinearity can be subjected to multilinear decomposition methods, mainly represented by PARAFAC and AMLD [20], which define the first family of models here described. They are based on a multilinear decomposition procedure, and present the following advantages: (1) no need of special initialization methods, since the solutions are often unique, (2) constraints may not be necessary, in general, to drive the optimization phase to the final solution, (3) figures of merit are well-known and (4) the decomposition of multi-dimensional data in their original structure is known to be more efficient than unfolding the data into arrays of lower dimensions [32]. The uniqueness property which is often achieved in this kind of decomposition is directly tied to the possibility of exploiting the relevant second-order advantage in analytical protocols.
Given a set of I bilinear second-order data with J and K data points in each of the instrumental modes and N responsive constituents, a three-way data arrangement X _ of size I × J × K can be obtained. If the tensorial object X _ obeys the property of low-rank trilinearity, then each element can be expressed as:
x i j k = n = 1 N a i n b j n c k n + e i j k
where a i n , b j n and c k n are the i th, j th and k th elements of the profile matrices A I × N , B J × N and C K × N in the n th column vectors, respectively. The scalar e i j k denotes a generic element of the three-way residual array E _ I × J × K . Under this assumption, X _ arrangement can be submitted to trilinear decomposition, according to the following model formulation:
X _ = n = 1 N a n b n c n + E _
where indicates the Kronecker product [9] and an, bn and cn are the nth. columns of A, B and C, respectively.
The same model formulation expressed by Equation (2). can be directly extended to a multilinear multi-way data arrangement of any order. For instance, in the case of third-order data, i.e., four-way calibration, and additional instrumental mode is included. If L data points are registered, then, quadrilinear decomposition follows the expression:
X _ = n = 1 N a n b n c n d n + E _
where D L × N contains the profiles of N constituents in the third instrumental mode and dn is the nth. column of D.
For three-way calibration, various optimization models for the estimation of A, B and C matrices have been proposed. In this sense, direct trilinear decomposition method (DTLD) [33] is based on eigenvalue-eigenvector decomposition, whereas PARAFAC [9] and AMLD [7] are based on alternating least-squares philosophy. In particular, AMLD includes a variety of strategies, such as alternating trilinear decomposition (ATLD) [34] or self-weighted ATLD (SWATLD) [35] for second-order data, and alternating penalty quadrilinear decomposition (APQLD) [36] or regularized self-weighted alternating quadrilinear decomposition (RSWAQLD) [37] for third-order data. The main advantage of these models, in general, relies in the fact that the multilinear decomposition is unique. From the analytical point of view, this property implies that, even when initial estimators are unknown, the optimization phase would retrieve the true constituent profiles along each of the instrumental modes, as well as the relative contribution of each component along the sample mode. The latter are then used to build the so-called pseudo-univariate calibration curve for quantitative purposes.
Notwithstanding, all these benefits can be achieved only if the modes are mutually independent from each other. This is not generally the common rule for chromatographic multi-way data. As stated above, multi-way chromatographic data are normally classified as non-trilinear type 1 because the reproducibility across sample injections is never perfect. In general, multilinear decomposition is the least flexible option, since it is sensitive to the lack of multilinearity [21]. This issue can be overcome by implementing different strategies. In certain cases, the lack of reproducibility in the sample mode can be solved by applying pre-processing procedures for chromatographic peak alignment [38]. Besides, different MLD variants derived from the previously mentioned approaches have also emerged aiming at dealing with the lack of multilinearity in chromatographic data. For instance, PARAFAC2 is a variant of the classical PARAFAC, which is based on an alternative form of Equation (2) [12]:
x i j k = n = 1 N a i n b j n ( i ) c k n + e i j k
where b j n ( i ) is the jth. element of the bn profile along the sample dependent mode and depends on the sample index i . In addition, to follow Equation (4), the PARAFAC2 model requires constant cross product between all pairs of bn profiles. This model formulation allows the chromatographic profiles not to be identical from run to run. For more details the reader must be referred to Ref [12]. Despite that this method proved to have some success in the field of chromatography, it is still a non-flexible model since it assumes that the degree of overlap between chromatographic peaks is constant for each pair of constituents in all samples [12,39].
The second group of multi-way calibration methodologies is represented by bidimensional decomposition and curve resolution methods, where MCR-ALS [40] emerges as the most important modelling approach for chromatographic multi-way data. Since its original publication, MCR-ALS has been the subject of extensive research in both fundamental and applied chemometrics. To put it succinctly, in contrast to multilinear decomposition methods, the MCR-ALS approach consists in performing a bilinear decomposition of a bilinear data arrangement X of size J × K , according to:
X = C S T + E
where C J × N and S K × N are matrices that capture the pure instrumental responses of N components in each of the instrumental modes, respectively, and E J × K collects the model residuals. The bilinear decomposition is achieved by virtue of the ALS algorithm. Besides, in contrast to the previously mentioned methods, ALS initialization with random values is not a convenient strategy. On the contrary, MCR-ALS is a more flexible model but suffers from ambiguity phenomena, which can have dramatic effects on the analytical performance of a calibration protocol [41,42]. Hence, in the usual application of this model, initial estimates of C or S are commonly obtained through the so-called purest variables methodology [43]. On the other hand, the study of MCR-ALS ambiguity is still a matter of extensive research in the chemometric field [41]. In general terms, the extent of RA can be mitigated through the incorporation of mathematical and chemically sensible constraints during ALS optimization, such as non-negativity, unimodality, selectivity, correspondence between species, among others [44,45,46,47,48,49,50,51].
For the particular case of three-way calibration, the most common approach for chromatographic data modelling is based on generating an augmented data array, i.e., a set of I bilinear matrices X are disposed in a column-wise augmented arrangement X aug of size J × K N , by placing each individual matrix one below the other. Then, X aug is submitted to bilinear decomposition according to the extended model formulation of MCR-ALS [52]:
X aug = C aug S T + E aug
After convergence, S captures the profiles of N species in the non-augmented mode (generally, the spectral mode) which is common to all samples, whereas C aug comprises the profiles of N species in the augmented mode, in each of the submatrices of X aug (generally, the elution time mode). Additionally, the area under the profiles captured within C aug is tied to the relative contributions of the individual components in each submatrix, which are then coupled to a regression model for quantitative purposes.
The possibility of performing a bilinear decomposition of augmented data arrays constitutes the fundamental aspect that gives this model both the optimal flexibility and versatility to deal with multilinearity deviation problems that characterize chromatographic data. This fact is key to understand why MCR-ALS has originated a myriad of applications in analytical calibration and has been extended to different kind of data, including third-order data.
For the specific case of four-way calibration, an alternative strategy known as APARAFAC was proposed by Bortolato et al. [13] combining the benefits of the PARAFAC uniqueness and the flexibility of MCR-ALS. This latter model can be implemented through a typical ALS process for performing a trilinear decomposition of an augmented three-way array. In this way, it allows to deal with non-quadrilinear type 1 data. The corresponding model for the X _ aug array of size IJ × K × L a can be represented by:
X _ aug = A B 3 W ( D 3 W C 3 W ) T + E _ aug
where the results of the decomposition are collected into three loading matrices AB3W (IJ × N), C3W (K × N) and D3W (L × N). The model residuals are retained in the E _ aug (IJ × K × L) array. Here, AB3W collects the unfolded profiles along the augmented chromatographic elution time modes and brings the relative contribution of the individual constituent present in every sample. The remaining decomposition matrices contain the profiles that enable a qualitative interpretation; for example, for LC-EEM data C and D will contain the excitation and emission spectral profiles of each responsive component.
Finally, the third category of multi-way chemometric approaches is constituted by latent variable-based models, essentially, the regression variants U-PLS and N-PLS. The PLS model was originally conceived for first-order calibration, i.e., senso strictum is only able to exploit the first-order advantage [53]. However, the PLS philosophy was extended to higher-order calibration, where the second-order advantage can be achieved by coupling the model to the RML methodology [54].
In U-PLS/RML, during the calibration step, a classical PLS model is built with a set of unfolded multi-way data, disposed in a matrix X PLS , which is then subjected to the following decomposition [53]:
X PLS = P T T + E PLS
In Equation (8), P and T are the PLS loading and score matrices, respectively, which aim to maximize both the explained variance in the X PLS data and the covariance with the nominal analyte concentration in the calibration samples. E PLS captures the PLS model residuals [55]. The PLS model returns a vector of latent regression coefficients v , which is usually employed to make predictions in test samples according to:
y = v T t
where t is the score vector of a given test sample, obtained by projection of an unfolded test signal onto the space spanned by the PLS calibration loadings. If a given test sample contains unmodelled constituents, the PLS prediction residual might be abnormally large compared to the expected noise level [53]. Under these circumstances, the RML methodology intends to decompose the part of the test data unexplained by PLS, assuming that the residuals can be rearranged into a multilinear array. For second-order calibration, RML is usually referred as RBL (bilinear) and its mathematical expression for a given test sample X test can be formulated as [54]:
X test = reshape ( P t RBL ) + B RBL T RBL T + E RBL
where B RBL T RBL T derives from the PCA model for the residual matrix X test reshape ( P t RBL ) with n RBL principal components. During the RBL procedure, a new sample score vector t RBL is calculated which only represents the analyte information. The t RBL vector is then used to make analyte predictions through Equation (9). RTL is a natural extension of Equation (10) for third-order data. On the other hand, N-PLS/RML represents a multi-way variant of PLS and is based on an analogous fundamental to those of U-PLS/RML. The key difference is that the original multi-way data structure is preserved during the PLS calibration/prediction stages.
From the qualitative point of view, in contrast to other model families, PLS loadings have no direct chemical interpretability. However, both latent models have shown to be advantageous in specific calibration scenarios. In particular, U-PLS/RML is the most flexible model and can be appropriate to model certain types of non-multilinear data where the lack of multilinearity occurs in more than one experimental mode. In principle, this is partially true if the second- or third-order advantage is to be achieved. Although PLS can satisfactory model a non-multilinear calibration dataset, RML procedures may fail during sample prediction if the residual matrix is not multilinear. In addition, for the particular case of chromatographic data, these methods are only applicable if the shifts in chromatographic profiles among runs are small. This means that U-PL/RML and N-PLS/RML are sensitive to the lack of multilinearity along the sample mode.
In order to graphically summarize all the information described in Section 2 and Section 3, the flow-chart shown in Figure 3 describes the most important modelling approaches for calibration purposes with chromatographic multi-way data, which have been here considered. If the generated data are multilinear, any of the presented models can in principle be implemented. However, multilinear models such as PARAFAC and AMLD variants should in this case be the first choice. When the data are not multilinear, or multilinearity cannot be restored by implementing convenient pre-processing procedures, MCR-ALS, PARAFAC2 or APARAFAC (for four-way arrays) are the most appropriate models to be applied. In any case, minor multilinear deviations can be tackled by PLS models. In all cases, the concentration of one or more analytes can be simultaneously obtained, even in the presence of unexpected sample constituents. In the case of PLS models, the concentration of one or more analytes are directly calculated by means of RML procedures. On the other, in multilinear decomposition and curve-resolution models, the relative contribution of the calibrated analytes is represented by the so-called component scores. They can be coupled to a pseudo-univariate regression model to estimate the analyte concentration in unknown samples.

4. Applications of LC Multi-Way Data

4.1. Second-Order/Three-Way Chromatographic Calibration

The continuous publication of scientific reports devoted to obtaining second-order/three-way chromatographic data and their valuable analytical applications makes it relevant to update the subject. As it is known, the analysis time demanded for both sample pre-treatment step and the chromatographic run itself is drastically reduced when multi-way calibration methods are implemented, while obtaining accurate quantitative information. Overcoming the problems inherent to chromatographic analysis of multi-component systems, such as co-elution of analytes and/or non-calibrated components and elution time shifts among runs, mostly depends on the proper selection of the chemometric model to be applied.
Table 1 summarizes some examples reported from 2018 to date developed for the analyte determinations based on second-order/three-way chromatographic calibration. The investigated analytes and matrices, the model/s selected for data processing, some relevant characteristics of the system and the attained limits of detection are indicated. As can be seen, the most widely applied approach for second-order/three-way chromatographic data generation is the popular LC-DAD. This may be due to a number of reasons, e.g., (1) most organic compounds absorb electromagnetic radiation in some region of the UV-visible spectrum, (2) modern DAD in chromatographic instruments enjoy a great sensitivity and versatility, and (3) it requires accessible equipment and low-cost consumables.
On the other hand, since not all compounds show fluorescent properties, it is not surprising that the number of publications describing chromatographic second-order calibration with fluorescence detection is much smaller. As will be seen below, in multivariate calibration, fluorescent detection has been more popular for obtaining third-order chromatographic data, where the EEM provides two instrumental modes while the third mode is given by the elution time. This detection mode takes advantage of the maximum selectivity and sensitivity which is gained by scanning both the emission and the excitation spectra, in comparison with fluorescence detection at a single excitation wavelength. In the latter case, a compromise is required among the excitation maxima of all analytes.
Finally, a smaller number of studies applying second-order LC data with MS detection have been published in the evaluated time period in comparison with DAD, despite the great potential of the former methodology. Most GC-MS second-order calibration publications shown in Table 1 involved quantitative studies of the migration of selected analytes from packaging material to foods, food simulants and cosmetic creams. LC-MS was applied to a variety of multi-components systems and complex samples.
Regardless of the detection method, the most used models for data processing were PARAFAC, ATLD and MCR-ALS. The advantages obtained when carrying out the chromatographic quantitation of analytes in complex samples using second-order data, as well as the problems that arise in each analysis, can be found in Table 1 for each report.

4.2. Third-Order/Four-Way Chromatographic Calibration

The quantitation of analytes through third-order/four-way data is not extensively applied by the analytical community. Although some works on this topic have recently been developed and their remarkable advantages as powerful analytical tools have been steadily highlighted, the number of bibliographic citations is still significantly lower than those corresponding to second-order calibration.
Applications of this type of high-order quantitation in the evaluate period are shown in Table 2. As can be appreciated, the most frequent way to generate third-order/four-way chromatographic data consists in measuring LC-EEM data.
While fluorescence matrices are easily obtained in fluorescent systems, their coupling to the chromatographic elution time mode must be carefully evaluated in order to generate a four-way array that complies with the properties required by the selected model for successful data processing [1,28,93,94,95]. The generation of photoinduced fluorescence in systems with analytes that do not present native fluorescence has been achieved through the incorporation of a UV reactor in the chromatographic equipment [93].
The research group led by Synovec has carried out an important contribution in the gas chromatography area for obtaining high-quality third-order data involving GC2-TOFMS. In recent years, the following works may be mentioned: (1) the study of the influence of the column selection and modulation period in the trilinear deviation ratio (TDR) range and, consequently, in the correct PARAFAC deconvolution [30]; (2) the use of a pulse flow valve for an ultra-fast modulation period and an innovative data processing method coupled with MCR-ALS deconvolution [96]; and (3) the advantages of partial modulation in the negative pulse mode [97]. In addition, valuable improvements to comprehensive GC3-FID detection have been achieved by the same research group [98,99].

5. Analytical Figures of Merit

Ideally, every development of a new calibration protocol should be accompanied by a report providing the AFOMs. This is usually the case in univariate calibration, but unfortunately the practice is not universally extended to multi-way calibration procedures. The subject has been thoroughly reviewed in 2014 [14], but new developments have taken place in recent years, which deserve to be commented in the present review.
Since we advocate for the use of MCR-ALS as the model of choice for multi-way chromatographic data, the AFOMs will be discussed in the context of this chemometric tool [100]. Specifically, the component-wise classical parameters sensitivity (SENn), analytical sensitivity (gn), selectivity (SELn), limit of detection (LODn) and limit of quantitation (LOQn) will be discussed, together with a recently proposed figure of merit which is characteristic of MCR-ALS: the rotational ambiguity derived uncertainty (RMSERA) [101].
The SENn can be described in qualitative terms as the net analyte signal at unit concentration, which can be estimated for a specific sample component in MCR-ALS as [14]:
SEN n = m n J 1 / 2 {   δ n T   [   S exp T (   I S unx S unx + )   S exp ] 1 δ n } 1 / 2
where J is the number of sensors of each sub-matrix in the augmented mode, mn is the slope of the pseudo-univariate calibration line, Sexp is the matrix of profiles in the non-augmented mode for expected constituents in the calibration set, Sunx is the matrix of profiles in the non-augmented mode for the unexpected constituents (interferents) and dn is an analyte selector vector, having a value of 1 at the analyte index and 0s otherwise.
The Equation (10) is written in such a way that it resembles the SENn definition in terms of net analyte signal. An equivalent expression can be written in a more compact manner as [100]:
SEN n = m n J 1 / 2 ( S T S ) n n 1 / 2
where the subscript ‘nn’ indicates the (n,n) element of a matrix.
More useful than the plain SENn is the analytical sensitivity γn, defined as the ration between SENn and the instrumental noise level σx [102]:
γ n = SEN n σ x
because it has inverse concentration units and is independent on the type of measured signal.
The SELn, in turn, can be estimated as the ratio between the SENn and the slope of the analyte calibration graph:
SEL n = J 1 / 2 SEN n m n
The degree by which the product (J1/2 SENn) departs from mn in the latter equation depends on the level of overlapping among the component profiles. The value of SELn varies between 0 (null selectivity) and 1 (100% or full selectivity). In the chromatographic context, multi-way selectivity has shown to be directly related to the effective peak capacity of a chromatogram [103].
The definition of the LODn has been changing with time [104]. The old concept of LODn as three times the standard deviation for a blank sample has been abandoned by IUPAC, in favour of the modern view of the LODn as a function of Type I and Type II errors (also called a and b errors, or false positive and false negative errors) [105,106]. In addition, today it is widely accepted that the uncertainties in the measurement of the test sample should be added to those stemming from the uncertainties affecting the calibration procedure (including both the measurement of the instrumental signal and the preparation of the calibration standards).
Since MCR-ALS quantitation is usually performed through a pseudo-univariate procedure in which the test sample score is interpolated in a regression line of calibration scores against nominal analyte concentrations, the LODn can be estimated from the latter line by extending the procedure, which is usual in univariate calibration, i.e.,
LOD n = 3 . 3 ( SEN n 2 σ x 2 + h 0   SEN n 2 σ x 2 + h 0 σ y cal 2 ) 1 / 2
where h0 is the leverage for the blank sample, and the factor 3.3 is equal to (ta,n + tb,n) for a = 0.05 and b = 0.05 and a large value of n (a and b are the probabilities for Type I and Type II errors, respectively and n is the number of degrees of freedom), sx is a measure of the uncertainty in the instrumental signal, which is considered to stem from identically and independently distributed noise in the signal measurements, and sycal is the concentration uncertainty when preparing the calibration samples. The factor 3.3 in the latter equation may be changed, if needed, for other values of a and b. Notice the that the assumptions underlying the expression for the limit of detection are: (1) the LODn is close enough to the blank so that the leverage can be considered equal to h0, and (2) the distance from the blank to the LODn is a sum of two confidence intervals; a more rigorous treatment suggests the use of a non-centrality parameter of a non-central t distribution instead of a sum of classical t-coefficients [107].
Under similar assumptions, the LODn is defined as the concentration for which the relative prediction uncertainty is 10%, leading to:
LOQ n = 10 ( SEN n 2 σ x 2 + h 0   SEN n 2 σ x 2 + h 0 σ y cal 2 ) 1 / 2
Specifically, for MCR-ALS, a recently proposed figure of merit is required for analytical reports based on the use of this model: the rotational ambiguity uncertainty. The latter stems from the fact that all bilinear decomposition solutions, such as those provided by MCR-ALS, are not necessarily unique, even when all possible constraints are applied during the decomposition [41]. When this is the case, i.e., for non-unique solutions, an uncertainty remains in the estimation of the analyte concentration, because a range of feasible solutions are possible [108]. The RMSERA derived from the existence of the feasible solutions has been shown to be given by the following range [101]:
R M S E RA = [ δ RA 12 ; δ RA 3 ]
where δRA is the width of the range of feasible concentration values for the analyte, which, in turn, is equal to the difference between maximum and minimum areas of the analyte concentration profile, converted to concentration units through the calibration slope:
δ RA = m a x ( a test ) m i n ( a test ) m n
It is important to notice that, even when a significant RMSERA value can be computed by the latter equation, the MCR-ALS solutions may be driven to the correct bilinear solution by an adequate initialization procedure, as recently discussed [109]. In these cases, the value of RMSERA may be large, although this would not be reflected in the average prediction error for a set of test samples, which may ultimately be reasonably low. However, it is important to notice that for future test samples, having a varying chemical composition in terms of interferents, the chosen initialization may not lead to similarly good results. The subject has been explained in detail in a recent tutorial [41].
Finally, if PARAFAC or other AMLD models are employed to process multi-way data of chromatographic origin, the decomposition is often unique, with AFOMs that have already been discussed in detail [14,110]. These models may successfully process data where the position and shapes of the chromatographic peaks do not change across samples, but as explained above it cannot be considered as the tool of choice for the presently discussed multi-way data.

6. An Example Comparing Second- and Third-Order Data

It is known, from theoretical considerations, that the SENn and SELn should increase with the increasing number of data modes when decomposing a multi-way array [14]. However, only in a few cases, this has been experimentally confirmed, by registering, for the same analytical system, second- and third-order data, and then comparing the resulting figures of merit [25].
Recently, such a comparison was made in connection with the simultaneous quantitation of various analytes in complex food samples [94]. The applied protocols were based on LC with fluorescence detection, employing two different detection setups for the same chromatograph. In one of them, second-order data were measured with fluorescence emission detection at a fixed excitation wavelength. In the second one, third-order data were collected by EEM detection. The analytes were ten different quinolone antibiotics, which were simultaneously analysed in edible animal tissues (chicken liver, bovine liver and bovine kidney). MCR-ALS was applied to augmented matrices for each test sample along the elution time direction. The use of MCR-ALS as data processing model provided excellent results with second-order data; the relative prediction errors spanned the range 4–12% in the case of the real samples at analyte concentrations, which are compatible with the maximum residue levels accepted by international regulations.
Third-order data were expected to provide better results, on account of the known sensitivity and selectivity increase by increasing the data order. In this case, the raw three-way arrays for each sample were unfolded into a matrix, by concatenating the emission and excitation modes into a single one prior to building the augmented data matrix along the elution time direction. The authors noticed that the third-order MCR-ALS analytical results were, in general, worse than for second-order data, with relative errors in the range 9–23%, and significantly larger detection limits for some analytes (Table 3). Furthermore, one of the analytes could not be resolved. Even when the theoretically estimated SENn was larger for third-order data (Table 3), as expected from the increase in the number of data modes, the relevant statistical indicators for analyte estimation (average prediction error and relative error of prediction) and the detection capabilities were indeed worse (Table 3).
Subsequently, the latent structured models U-PLS/RBL and U-PLS/RTL were applied to the second- and third-order data, respectively [94]. In this case, relative errors in the range 7–18% were obtained for second-order calibration and 5–27% for third-order calibration. The latter errors were significantly larger than those for MCR-ALS/second-order data, although the U-PLS/RTL model permitted the detection of all the studied analytes.
It was concluded from the comparison of these two different experimental methodologies, carried out by using different detectors, that matrix fluorescence detection in LC would still require some advances in instrumental setup, in order to provide maximum sensitivity and signal/noise ratio at very high scanning rates, something which is not currently possible with commercially available equipment. In the above commented report, the second-order data were measured using a commercial instrument which is fully optimized for maximum sensitivity, whereas third-order data collection required the in-house connection of a LC to a conventional spectrofluorometer, which is not optimized for emission measurements in a small flowing cell.
This example also illustrates the fact that the plain sensitivity parameter SENn may not be used for a proper comparison of analytical methods based on data measured with different instruments. As explained above, a better AFOM is the γn. The effect of a larger noise level could be the ultimate reason why third-order calibration failed to provide better analytical results in the example of ref. [94].

7. Conclusions

Chromatography is an ever-growing discipline. Advances in experimental activities have demanded the development of new chemometric models, whose ultimate aim is to extract useful information for developing powerful analytical protocols. This is particularly true in the area of multi-way chromatography, which is able to produce large data sets that can be arranged into mathematical objects of increasing complexity. Selecting the appropriate multi-way chemometric model to process these data sets is crucial for achieving one of the most important advantages offered by this field, i.e., the possibility of quantitating analytes in complex samples containing uncalibrated interferents (the second-order advantage). The present review provides a current view of the state-of-the-art, both from experimental and theoretical perspectives.

Author Contributions

F.A.C., M.R.A., G.M.E., H.C.G. and A.C.O. have equally contributed to writing the manuscript. A.C.O. has revised and approved it for submission. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors are grateful to CONICET (Consejo Nacional de Investigaciones Científicas y Técnicas, Project PIP 2015-0111) and ANPCyT (Agencia Nacional de Promoción Científica y Tecnológica, Projects PICT 2017-0340, PICT 2018-02533 and PICT 2016-1122) for financial support. F.A.C. thanks CONICET for his fellowship.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Escandar, G.M.; Olivieri, A.C. Multi-way chromatographic calibration—A review. J. Chromatogr. A 2019, 1587, 2–13. [Google Scholar] [CrossRef]
  2. Santana, I.M.; Breitkreitz, M.C.; Pinto, L. Multivariate Curve Resolution Alternating Least Squares Applied to Chromatographic Data: From the Basics to the Recent Advances. Braz. J. Anal. Chem. BrJAC 2021, 8, 22–44. [Google Scholar]
  3. Anzardi, M.B.; Arancibia, J.A.; Olivieri, A.C. Processing multi-way chromatographic data for analytical calibration, classification and discrimination: A successful marriage between separation science and chemometrics. TrAC Trends Anal. Chem. 2021, 134, 116128. [Google Scholar] [CrossRef]
  4. Booksh, K.S.; Kowalski, B.R. Theory of Analytical Chemistry. Anal. Chem. 1994, 66, 782A–791A. [Google Scholar] [CrossRef]
  5. Escandar, G.M.; Olivieri, A.C.; Faber, N.M.; Goicoechea, H.C.; Muñoz de la Peña, A.; Poppi, R.J. Second- and third-order multivariate calibration: Data, algorithms and applications. TrAC Trends Anal. Chem. 2007, 26, 752–765. [Google Scholar] [CrossRef]
  6. Alcaraz, M.R.; Monago-Maraña, O.; Goicoechea, H.C.; Muñoz de la Peña, A. Four- and five-way excitation-emission luminescence-based data acquisition and modeling for analytical applications. A review. Anal. Chim. Acta 2019, 1083, 41–57. [Google Scholar] [CrossRef]
  7. Wu, H.-L.; Wang, T.; Yu, R.-Q. Recent advances in chemical multi-way calibration with second-order or higher-order advantages: Multilinear models, algorithms, related issues and applications. TrAC Trends Anal. Chem. 2020, 130, 115954. [Google Scholar] [CrossRef]
  8. Appellof, C.J.; Davidson, E.R. Strategies for analyzing data from video fluorometric monitoring of liquid chromatographic effluents. Anal. Chem. 1981, 53, 2053–2056. [Google Scholar] [CrossRef]
  9. Bro, R. PARAFAC. Tutorial and applications. Chemom. Intell. Lab. Syst. 1997, 38, 149–171. [Google Scholar] [CrossRef]
  10. Tauler, R. Multivariate curve resolution applied to second order data. Chemom. Intell. Lab. Syst. 1995, 30, 133–146. [Google Scholar] [CrossRef]
  11. Muñoz de la Peña, A.; Goicoechea, H.C.; Escandar, G.M.; Olivieri, A.C. Fundamentals and Analytical Applications of Multiway Calibration; Elservier: Amsterdam, The Netherlands, 2015. [Google Scholar]
  12. Kiers, H.A.L.; ten Berge, J.M.F.; Bro, R. PARAFAC2—Part I. A direct fitting algorithm for the PARAFAC2 model. J. Chemom. 1999, 13, 275–294. [Google Scholar] [CrossRef]
  13. Bortolato, S.A.; Lozano, V.A.; Muñoz de la Peña, A.; Olivieri, A.C. Novel augmented parallel factor model for four-way calibration of high-performance liquid chromatography–fluorescence excitation–emission data. Chemom. Intell. Lab. Syst. 2015, 141, 1–11. [Google Scholar] [CrossRef]
  14. Olivieri, A.C. Analytical Figures of Merit: From Univariate to Multiway Calibration. Chem. Rev. 2014, 114, 5358–5378. [Google Scholar] [CrossRef] [PubMed]
  15. Murphy, K.R.; Stedmon, C.A.; Graeber, D.; Bro, R. Fluorescence spectroscopy and multi-way techniques. PARAFAC. Anal. Methods 2013, 5, 6557–6566. [Google Scholar] [CrossRef] [Green Version]
  16. Olivieri, A.C.; Escandar, G.M. Parallel Factor Analysis: Trilinear Data. In Practical Three-Way Calibration; Elsevier: Waltham, MA, USA, 2014. [Google Scholar]
  17. Piccirilli, G.N.; Escandar, G.M. Partial least-squares with residual bilinearization for the spectrofluorimetric determination of pesticides. A solution of the problems of inner-filter effects and matrix interferents. Analyst 2006, 131, 1012–1020. [Google Scholar] [CrossRef]
  18. Benavente, F.; Andón, B.; Giménez, E.; Olivieri, A.C.; Barbosa, J.; Sanz-Nebot, V. A multiway approach for classification and characterization of rabbit liver apothioneins by CE-ESI-MS. Electrophoresis 2008, 29, 4355–4367. [Google Scholar] [CrossRef]
  19. Prebihalo, S.E.; Berrier, K.L.; Freye, C.E.; Bahaghighat, H.D.; Moore, N.R.; Pinkerton, D.K.; Synovec, R.E. Multidimensional Gas Chromatography: Advances in Instrumentation, Chemometrics, and Applications. Anal. Chem. 2018, 90, 505–532. [Google Scholar] [CrossRef]
  20. Chen, Y.; Zou, C.M.; Bin, J.; Yang, M.; Kang, C. Multilinear Mathematical Separation in Chromatography. Separations 2021, 8, 31. [Google Scholar] [CrossRef]
  21. Bro, R. Multi-Way Analysis in the Food Industry; Royal Veterinary and Agricultural University: Compenhagen, Denmark, 1997. [Google Scholar]
  22. Alcaraz, M.R.; Siano, G.G.; Culzoni, M.J.; Muñoz de la Peña, A.; Goicoechea, H.C. Modeling four and three-way fast high-performance liquid chromatography with fluorescence detection data for quantitation of fluoroquinolones in water samples. Anal. Chim. Acta 2014, 809, 37–46. [Google Scholar] [CrossRef]
  23. Lozano, V.A.; Muñoz de la Peña, A.; Durán-Merás, I.; Espinosa Mansilla, A.; Escandar, G.M. Four-way multivariate calibration using ultra-fast high-performance liquid chromatography with fluorescence excitation–emission detection. Application to the direct analysis of chlorophylls a and b and pheophytins a and b in olive oils. Chemom. Intell. Lab. Syst. 2013, 125, 121–131. [Google Scholar] [CrossRef]
  24. Montemurro, M.; Pinto, L.; Véras, G.; de Araújo Gomes, A.; Culzoni, M.J.; Ugulino de Araújo, M.C.; Goicoechea, H.C. Highly sensitive quantitation of pesticides in fruit juice samples by modeling four-way data gathered with high-performance liquid chromatography with fluorescence excitation-emission detection. Talanta 2016, 154, 208–218. [Google Scholar] [CrossRef]
  25. Montemurro, M.; Siano, G.G.; Alcaraz, M.R.; Goicoechea, H.C. Third order chromatographic-excitation–emission fluorescence data: Advances, challenges and prospects in analytical applications. TrAC Trends Anal. Chem. 2017, 93, 119–133. [Google Scholar] [CrossRef]
  26. Carabajal, M.D.; Arancibia, J.A.; Escandar, G.M. On-line generation of third-order liquid chromatography–excitation-emission fluorescence matrix data. Quantitation of heavy-polycyclic aromatic hydrocarbons. J. Chromatogr. A 2017, 1527, 61–69. [Google Scholar] [CrossRef] [PubMed]
  27. Alcaraz, M.R.; Morzán, E.; Sorbello, C.; Goicoechea, H.C.; Etchenique, R. Multiway analysis through direct excitation-emission matrix imaging. Anal. Chim. Acta 2018, 1032, 32–39. [Google Scholar] [CrossRef]
  28. Carabajal, M.D.; Arancibia, J.A.; Escandar, G.M. Multivariate curve resolution strategy for non-quadrilinear type 4 third-order/four way liquid chromatography–excitation-emission fluorescence matrix data. Talanta 2018, 189, 509–516. [Google Scholar] [CrossRef]
  29. Bailey, H.P.; Rutan, S.C. Chemometric Resolution and Quantification of Four-Way Data Arising from Comprehensive 2D-LC-DAD Analysis of Human Urine. Chemom. Intell. Lab. Syst. 2011, 106, 131–141. [Google Scholar] [CrossRef] [Green Version]
  30. Prebihalo, S.E.; Pinkerton, D.K.; Synovec, R.E. Impact of comprehensive two-dimensional gas chromatography time-of-flight mass spectrometry experimental design on data trilinearity and parallel factor analysis deconvolution. J. Chromatogr. A 2019, 1605, 460368. [Google Scholar] [CrossRef]
  31. Ferreira, V.H.C.; Hantao, L.W.; Poppi, R.J. Consumable-free Comprehensive Three-Dimensional Gas Chromatography and PARAFAC for Determination of Allergens in Perfumes. Chromatographia 2020, 83, 581–592. [Google Scholar] [CrossRef]
  32. Xiangqian, L.; Sidiropoulos, N.D. Cramer-Rao lower bounds for low-rank decomposition of multidimensional arrays. IEEE Trans. Signal Process. 2001, 49, 2074–2086. [Google Scholar] [CrossRef]
  33. Sanchez, E.; Kowalski, B.R. Tensorial resolution: A direct trilinear decomposition. J. Chemom. 1990, 4, 29–45. [Google Scholar] [CrossRef]
  34. Wu, H.-L.; Shibukawa, M.; Oguma, K. An alternating trilinear decomposition algorithm with application to calibration of HPLC–DAD for simultaneous determination of overlapped chlorinated aromatic hydrocarbons. J. Chemom. 1998, 12, 1–26. [Google Scholar] [CrossRef]
  35. Chen, Z.-P.; Wu, H.-L.; Jiang, J.-H.; Li, Y.; Yu, R.-Q. A novel trilinear decomposition algorithm for second-order linear calibration. Chemom. Intell. Lab. Syst. 2000, 52, 75–86. [Google Scholar] [CrossRef]
  36. Xia, A.L.; Wu, H.-L.; Li, S.-F.; Zhu, S.-H.; Hu, L.-Q.; Yu, R.-Q. Alternating penalty quadrilinear decomposition algorithm for an analysis of four-way data arrays. J. Chemom. 2007, 21, 133–144. [Google Scholar] [CrossRef]
  37. Kang, C.; Wu, H.-L.; Yu, Y.-J.; Liu, Y.-J.; Zhang, S.-R.; Zhang, X.-H.; Yu, R.-Q. An alternative quadrilinear decomposition algorithm for four-way calibration with application to analysis of four-way fluorescence excitation–emission–pH data array. Anal. Chim. Acta 2013, 758, 45–57. [Google Scholar] [CrossRef] [PubMed]
  38. Yu, Y.-J.; Wu, H.-L.; Niu, J.-F.; Zhao, J.; Li, Y.-N.; Kang, C.; Yu, R.-Q. A novel chromatographic peak alignment method coupled with trilinear decomposition for three dimensional chromatographic data analysis to obtain the second-order advantage. Analyst 2013, 138, 627–634. [Google Scholar] [CrossRef] [PubMed]
  39. Olivieri, A.C.; Escandar, G.M. Parallel Factor Analysis: Nontrilinear Data of Type 1. In Practical Three-Way Calibration; Elsevier: Boston, MA, USA, 2014; Chapter 7. [Google Scholar]
  40. Tauler, R.; Maeder, M.; de Juan, A. 2.15-Multiset Data Analysis: Extended Multivariate Curve Resolution. In Comprehensive Chemometrics, 2nd ed.; Brown, S., Tauler, R., Walczak, B., Eds.; Elsevier: Oxford, UK, 2020. [Google Scholar]
  41. Olivieri, A.C. A down-to-earth analyst view of rotational ambiguity in second-order calibration with multivariate curve resolution—A tutorial. Anal. Chim. Acta 2021, 1156, 338206. [Google Scholar] [CrossRef] [PubMed]
  42. Golshan, A.; Abdollahi, H.; Beyramysoltan, S.; Maeder, M.; Neymeyr, K.; Rajkó, R.; Sawall, M.; Tauler, R. A review of recent methods for the determination of ranges of feasible solutions resulting from soft modelling analyses of multivariate data. Anal. Chim. Acta 2016, 911, 1–13. [Google Scholar] [CrossRef]
  43. Windig, W.; Guilment, J. Interactive self-modeling mixture analysis. Anal. Chem. 1991, 63, 1425–1432. [Google Scholar] [CrossRef]
  44. Tauler, R.; de Juan, A. Multivariate Curve Resolution for Quantitative Analysis. In Data Handling in Science and Technology; Muñoz de la Peña, A., Goicoechea, H.C., Escandar, G.M., Olivieri, A.C., Eds.; Elsevier: Amsterdam, The Netherlands, 2015; Chapter 5. [Google Scholar]
  45. Ahmadi, G.; Tauler, R.; Abdollahi, H. Multivariate calibration of first-order data with the correlation constrained MCR-ALS method. Chemom. Intell. Lab. Syst. 2015, 142, 143–150. [Google Scholar] [CrossRef]
  46. Akbari Lakeh, M.; Abdollahi, H. Known-value constraint in multivariate curve resolution. Anal. Chim. Acta 2018, 1030, 42–51. [Google Scholar] [CrossRef]
  47. Akbari Lakeh, M.; Abdollahi, H.; Gemperline, P.J. Soft known-value constraints for improved quantitation in multivariate curve resolution. Anal. Chim. Acta 2020, 1105, 64–73. [Google Scholar] [CrossRef] [PubMed]
  48. Tauler, R.; Marqués, I.; Casassas, E. Multivariate curve resolution applied to three-way trilinear data: Study of a spectrofluorimetric acid–base titration of salicylic acid at three excitation wavelengths. J. Chemom. 1998, 12, 55–75. [Google Scholar] [CrossRef]
  49. Alier, M.; Felipe-Sotelo, M.; Hernández, I.; Tauler, R. Trilinearity and component interaction constraints in the multivariate curve resolution investigation of NO and O-3 pollution in Barcelona. Anal. Bioanal. Chem. 2011, 399, 2015–2029. [Google Scholar] [CrossRef] [PubMed]
  50. Tavakkoli, E.; Abdollahi, H.; Gemperline, P.J. Soft-trilinear constraints for improved quantitation in multivariate curve resolution. Analyst 2020, 145, 223–232. [Google Scholar] [CrossRef] [PubMed]
  51. de Juan, A.; Tauler, R. Multivariate Curve Resolution: 50 years addressing the mixture analysis problem—A review. Anal. Chim. Acta 2021, 1145, 59–78. [Google Scholar] [CrossRef]
  52. Tauler, R.; Maeder, M.; Juan, A. Multiset Data Analysis: Extended Multivariate Curve Resolution; Elsevier: Amsterdam, The Netherlands, 2009. [Google Scholar]
  53. Olivieri, A.C.; Escandar, G.M. Partial Least-Squares with Residual Bilinearization. In Practical Three-Way Calibration; Elsevier: Boston, MA, USA, 2014; Chapter 9. [Google Scholar]
  54. Olivieri, A.C. On a versatile second-order multivariate calibration method based on partial least-squares and residual bilinearization: Second-order advantage and precision properties. J. Chemom. 2005, 19, 253–265. [Google Scholar] [CrossRef]
  55. Haaland, D.M.; Thomas, E.V. Partial least-squares methods for spectral analyses. 1. Relation to other quantitative calibration methods and the extraction of qualitative information. Anal. Chem. 1988, 60, 1193–1202. [Google Scholar] [CrossRef]
  56. Mortera, P.; Zuljan, F.A.; Magni, C.; Bortolato, S.A.; Alarcón, S.H. Multivariate analysis of organic acids in fermented food from reversed-phase high-performance liquid chromatography data. Talanta 2018, 178, 15–23. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Wang, T.; Wu, H.-L.; Xie, L.-X.; Liu, Z.; Long, W.-J.; Cheng, L.; Ding, Y.-J.; Yu, R.-Q. Simultaneous and interference-free determination of eleven non-steroidal anti-inflammatory drugs illegally added into Chinese patent drugs using chemometrics-assisted HPLC-DAD strategy. Sci. China Chem. 2018, 61, 739–749. [Google Scholar] [CrossRef]
  58. Liu, Q.; Wu, H.-L.; Liu, Z.; Xiao, R.; Wang, T.; Hu, Y.; Ding, Y.-J.; Yu, R.-Q. Chemometrics-assisted HPLC-DAD as a rapid and interference-free strategy for simultaneous determination of 17 polyphenols in raw propolis. Anal. Methods 2018, 10, 5577–5588. [Google Scholar] [CrossRef]
  59. Monzón, C.M.; Teglia, C.M.; Delfino, M.R.; Goicoechea, H.C. Multiway calibration strategy with chromatographic data exploiting the second-order advantage for quantitation of three antidiabetic and three antihypertensive drugs in serum samples. Microchem. J. 2018, 136, 185–192. [Google Scholar] [CrossRef]
  60. Vosough, M.; Tehrani, S.M. Development of a fast HPLC-DAD method for simultaneous quantitation of three immunosuppressant drugs in whole blood samples using intelligent chemometrics resolving of coeluting peaks in the presence of blood interferences. J. Chromatogr. B 2018, 1073, 69–79. [Google Scholar] [CrossRef]
  61. Yin, X.-L.; Gu, H.-W.; Jalalvand, A.R.; Liu, Y.-J.; Chen, Y.; Peng, T.-Q. Dealing with overlapped and unaligned chromatographic peaks by second-order multivariate calibration for complex sample analysis: Fast and green quantification of eight selected preservatives in facial masks. J. Chromatogr. A 2018, 1573, 18–27. [Google Scholar] [CrossRef]
  62. Rezaei, F.; Sheikholeslami, M.; Vosough, M.; Maeder, M. Handling of highly coeluted chromatographic peaks by multivariate curve resolution for a complex bioanalytical problem: Quantitation of selected corticosteroids and mycophenolic acid in human plasma. Talanta 2018, 187, 1–12. [Google Scholar] [CrossRef]
  63. Yin, X.-L.; Wu, H.-L.; Gu, H.-W.; Zhang, X.-H.; Sun, Y.-M.; Hu, Y.; Liu, L.; Rong, Q.-M.; Yu, R.-Q. Chemometrics-enhanced high performance liquid chromatography-diode array detection strategy for simultaneous determination of eight co-eluted compounds in ten kinds of Chinese teas using second-order calibration method based on alternating trilinear decomposition algorithm. J. Chromatogr. A 2014, 1364, 151–162. [Google Scholar]
  64. Sun, X.-D.; Wu, H.-L.; Liu, Z.; Chen, Y.; Chen, J.-C.; Cheng, L.; Ding, Y.-J.; Yu, R.-Q. Target-based metabolomics for fast and sensitive quantification of eight small molecules in human urine using HPLC-DAD and chemometrics tools resolving of highly overlapping peaks. Talanta 2019, 201, 174–184. [Google Scholar] [CrossRef]
  65. Qing, X.-D.; Zhou, H.-B.; Zhang, X.-H.; Wu, H.-L.; Chen, C.-Y.; Xu, S.-W.; Li, S.-S. Alternating trilinear decomposition of highly overlapped chromatograms for simultaneously targeted quantification of 15 PAHs in samples of pollution source. Microchem. J. 2019, 146, 742–752. [Google Scholar] [CrossRef]
  66. Wang, T.; Wu, H.-L.; Yu, Y.-J.; Long, W.-J.; Cheng, L.; Chen, A.-Q.; Yu, R.-Q. A simple method for direct modeling of second-order liquid chromatographic data with retention time shifts and holding the second-order advantage. J. Chromatogr. A 2019, 1605, 360360. [Google Scholar] [CrossRef] [PubMed]
  67. Teglia, C.M.; Santamaría, C.G.; Rodriguez, H.A.; Culzoni, M.J.; Goicoechea, H.C. Determination of 2-hydroxy-4-methoxybenzophenone in mice serum and human plasma by ultra-high-performance liquid chromatography enhanced by chemometrics. Microchem. J. 2019, 148, 35–41. [Google Scholar] [CrossRef]
  68. Akvan, N.; Azimi, G.; Parastar, H. Chemometric assisted determination of 16 PAHs in water samples by ultrasonic assisted emulsification microextraction followed by fast high-performance liquid chromatography with diode array detector. Microchem. J. 2019, 150, 104056. [Google Scholar] [CrossRef]
  69. Dinç, E.; Büker, E. Spectrochromatographic determination of dorzolamide hydrochloride and timolol maleate in an ophthalmic solution using three-way analysis methods. Talanta 2019, 191, 248–256. [Google Scholar] [CrossRef] [PubMed]
  70. Sun, X.-D.; Wu, H.-L.; Liu, Z.; Chen, Y.; Liu, Q.; Ding, Y.-J.; Yu, R.-Q. Rapid and Sensitive Detection of Multi-Class Food Additives in Beverages for Quality Control by Using HPLC-DAD and Chemometrics Methods. Food Anal. Methods 2019, 12, 381–393. [Google Scholar] [CrossRef]
  71. Zhang, X.-H.; Zhou, Q.; Liu, Z.; Qing, X.-D.; Zheng, J.-J.; Mu, S.-T.; Liu, P.-H. Comparison of three second-order multivariate calibration methods for the rapid identification and quantitative analysis of tea polyphenols in Chinese teas using high-performance liquid chromatography. J. Chromatogr. A 2020, 1618, 460905. [Google Scholar] [CrossRef]
  72. Sousa, E.S.; Schneider, M.P.; Pinto, L.; de Araujo, M.C.U.; de Araújo Gomes, A. Chromatographic quantification of seven pesticide residues in vegetable: Univariate and multiway calibration comparison. Microchem. J. 2020, 152, 104301. [Google Scholar] [CrossRef]
  73. Zhang, X.-H.; Ma, Y.-X.; Yi, C.; Qing, X.-D.; Liu, Z.; Zheng, J.-J.; Lin, F.; Lv, T.-F. Chemometrics-enhanced HPLC–DAD as a rapid and interference-free strategy for simultaneous quantitative analysis of flavonoids in Chinese propolis. Eur. Food Res. Technol. 2020, 246, 1909–1918. [Google Scholar] [CrossRef]
  74. Ortiz, M.C.; Sanllorente, S.; Herrero, A.; Reguera, C.; Rubio, L.; Oca, M.L.; Valverde-Som, L.; Arce, M.M.; Sánchez, M.S.; Sarabia, L.A. Three-way PARAFAC decomposition of chromatographic data for the unequivocal identification and quantification of compounds in a regulatory framework. Chemom. Intell. Lab. Syst. 2020, 200, 104003. [Google Scholar] [CrossRef]
  75. Arce, M.M.; Ruiz, S.; Sanllorente, S.; Ortiz, M.C.; Sarabia, L.A.; Sánchez, M.S. A new approach based on inversion of a partial least squares model searching for a preset analytical target profile. Application to the determination of five bisphenols by liquid chromatography with diode array detector. Anal. Chim. Acta 2021, 1149, 338217. [Google Scholar] [CrossRef]
  76. Arce, M.M.; Ortiz, M.C.; Sanllorente, S. HPLC-DAD and PARAFAC for the determination of bisphenol-A and another four bisphenols migrating from BPA-free polycarbonate glasses. Microchem. J. 2021, 168, 106413. [Google Scholar] [CrossRef]
  77. Long, W.-J.; Wu, H.-L.; Wang, T.; Dong, M.-Y.; Yu, R.-Q. Interference-free analysis of multi-class preservatives in cosmetic products using alternating trilinear decomposition modeling of liquid chromatography diode array detection data. Microchem. J. 2021, 162, 105847. [Google Scholar] [CrossRef]
  78. Chang, Y.-Y.; Wu, H.-L.; Fang, H.; Wang, T.; Ouyang, Y.-Z.; Sun, X.-D.; Tong, G.-Y.; Ding, Y.-J.; Yu, R.-Q. Comparison of three chemometric methods for processing HPLC-DAD data with time shifts: Simultaneous determination of ten molecular targeted anti-tumor drugs in different biological samples. Talanta 2021, 224, 121798. [Google Scholar] [CrossRef]
  79. Carabajal, M.; Teglia, C.M.; Maine, M.A.; Goicoechea, H.C. Multivariate optimization of a dispersive liquid-liquid microextraction method for the determination of six antiparasite drugs in kennel effluent waters by using second-order chromatographic data. Talanta 2021, 224, 121929. [Google Scholar] [CrossRef]
  80. Wang, Z.-Y.; Wu, H.-L.; Chang, Y.-Y.; Wang, T.; Chen, W.; Tong, G.-Y.; Yu, R.-Q. Simultaneous determination of nine tyrosine kinase inhibitors in three complex biological matrices by using high-performance liquid chromatography–diode array detection combined with a second-order calibration method. J. Sep. Sci. 2021, 1–10. [Google Scholar]
  81. Anzardi, M.B.; Arancibia, J.A. Chemometrics-assisted liquid chromatographic determination of quinolones in edible animal tissues. Microchem. J. 2020, 158, 105138. [Google Scholar] [CrossRef]
  82. Rubio, L.; Sarabia, L.A.; Ortiz, M.C. Effect of the cleaning procedure of Tenax on its reuse in the determination of plasticizers after migration by gas chromatography/mass spectrometry. Talanta 2018, 182, 505–522. [Google Scholar] [CrossRef]
  83. Rubio, L.; Valverde-Som, L.; Sarabia, L.A.; Ortiz, M.C. Improvement in the identification and quantification of UV filters and additives in sunscreen cosmetic creams by gas chromatography/mass spectrometry through three-way calibration techniques. Talanta 2019, 205, 120156. [Google Scholar] [CrossRef]
  84. Rubio, L.; Valverde-Som, L.; Sarabia, L.A.; Ortiz, M.C. The behaviour of Tenax as food simulant in the migration of polymer additives from food contact materials by means of gas chromatography/mass spectrometry and PARAFAC. J. Chromatogr. A 2019, 1589, 18–29. [Google Scholar] [CrossRef]
  85. Valverde-Som, L.; Reguera, C.; Herrero, A.; Sarabia, L.A.; Ortiz, M.C. Determination of polymer additive residues that migrate from coffee capsules by means of stir bar sorptive extraction-gas chromatography-mass spectrometry and PARAFAC decomposition. Food Packag. Shelf Life 2021, 28, 100664. [Google Scholar] [CrossRef]
  86. Qing, X.; Zhou, X.; Xu, L.; Zhang, J.; Huang, Y.; Lin, L.; Liu, Z.; Zhang, X. Application of alternating trilinear decomposition-assisted multivariate curve resolution to gas chromatography-mass spectrometric data for the quantification of polycyclic aromatic hydrocarbons in aerosols. R. Soc. Open Sci. 2021, 8, 210458. [Google Scholar] [CrossRef] [PubMed]
  87. Long, W.-J.; Wu, H.-L.; Wang, T.; Xie, L.-X.; Hu, Y.; Fang, H.; Cheng, L.; Ding, Y.-J.; Yu, R.-Q. Chemometrics-assisted liquid chromatography with full scan mass spectrometry for the interference-free determination of glucocorticoids illegally added to face masks. J. Sep. Sci. 2018, 41, 3527–3537. [Google Scholar] [CrossRef] [PubMed]
  88. Hu, Y.; Wu, H.-L.; Yin, X.-L.; Gu, H.-W.; Xiao, R.; Xie, L.-X.; Liu, Z.; Fang, H.; Wang, L.; Yu, R.-Q. Rapid and interference-free analysis of nine B-group vitamins in energy drinks using trilinear component modeling of liquid chromatography-mass spectrometry data. Talanta 2018, 180, 108–119. [Google Scholar] [CrossRef] [PubMed]
  89. Sun, X.-D.; Wu, H.-L.; Liu, Z.; Xie, L.-X.; Hu, Y.; Fang, H.; Wang, T.; Xiao, R.; Ding, Y.-J.; Yu, R.-Q. Chemometrics-assisted liquid chromatography-full scan mass spectrometry for simultaneous determination of multi-class estrogens in infant milk powder. Anal. Methods 2018, 10, 1459–1471. [Google Scholar] [CrossRef]
  90. Long, W.-J.; Wu, H.-L.; Wang, T.; Dong, M.-Y.; Yu, R.-Q. Exploiting second-order advantage from mathematically modeled liquid chromatography–mass spectrometry data for simultaneous determination of polyphenols in Chinese propolis. Microchem. J. 2020, 157, 105003. [Google Scholar] [CrossRef]
  91. Sheikholeslami, M.N.; Vosough, M.; Esfahani, H.M. On the performance of multivariate curve resolution to resolve highly complex liquid chromatography–full scan mass spectrometry data for quantification of selected immunosuppressants in blood and water samples. Microchem. J. 2020, 152, 104298. [Google Scholar] [CrossRef]
  92. Sun, X.-D.; Wu, H.-L.; Chen, J.-C.; Chen, A.-Q.; Chen, Y.; Ouyang, Y.-Z.; Ding, Y.-J.; Yu, R.-Q. Exploration advantages of data combination and partition: First chemometric analysis of liquid chromatography–mass spectrometry data in full scan mode with quadruple fragmentor voltages. Anal. Chim. Acta 2020, 1110, 158–168. [Google Scholar] [CrossRef] [PubMed]
  93. Pellegrino Vidal, R.B.; Olivieri, A.C.; Ibañez, G.A.; Escandar, G.M. Online Third-Order Liquid Chromatographic Data with Native and Photoinduced Fluorescence Detection for the Quantitation of Organic Pollutants in Environmental Water. ACS Omega 2018, 3, 15771–15779. [Google Scholar] [CrossRef] [PubMed]
  94. Anzardi, M.B.; Arancibia, J.A.; Olivieri, A.C. Using chemometric tools to investigate the quality of three- and four-way liquid chromatographic data obtained with two different fluorescence detectors and applied to the determination of quinolone antibiotics in animal tissues. Chemom. Intell. Lab. Syst. 2020, 199, 103972. [Google Scholar] [CrossRef]
  95. Siano, G.G.; Vera Candioti, L.; Giovanini, L.L. Chemometric handling of spectral-temporal dependencies for liquid chromatography data with online registering of excitation-emission fluorescence matrices. Chemom. Intell. Lab. Syst. 2020, 199, 103961. [Google Scholar] [CrossRef]
  96. Bahaghighat, H.D.; Freye, C.E.; Gough, D.V.; Synovec, R.E. Comprehensive two-dimensional gas chromatography and time-of-flight mass spectrometry detection with a 50 ms modulation period. J. Chromatogr. A 2019, 1583, 117–123. [Google Scholar] [CrossRef]
  97. Gough, D.V.; Schöneich, S.; Synovec, R.E. Chemometric decomposition of comprehensive two-dimensional gas chromatography time-of-flight mass spectrometry data employing partial modulation in the negative pulse mode. Talanta 2020, 210, 120670. [Google Scholar] [CrossRef]
  98. Gough, D.V.; Bahaghighat, H.D.; Synovec, R.E. Column selection approach to achieve a high peak capacity in comprehensive three-dimensional gas chromatography. Talanta 2019, 195, 822–829. [Google Scholar] [CrossRef]
  99. Trinklein, T.J.; Prebihalo, S.E.; Warren, C.G.; Ochoa, G.S.; Synovec, R.E. Discovery-based analysis and quantification for comprehensive three-dimensional gas chromatography flame ionization detection data. J. Chromatogr. A 2020, 1623, 461190. [Google Scholar] [CrossRef]
  100. Bauza, M.C.; Ibañez, G.A.; Tauler, R.; Olivieri, A.C. Sensitivity Equation for Quantitative Analysis with Multivariate Curve Resolution-Alternating Least-Squares: Theoretical and Experimental Approach. Anal. Chem. 2012, 84, 8697–8706. [Google Scholar] [CrossRef] [PubMed]
  101. Pellegrino Vidal, R.B.; Olivieri, A.C. A New Parameter for Measuring the Prediction Uncertainty Produced by Rotational Ambiguity in Second-Order Calibration with Multivariate Curve Resolution. Anal. Chem. 2020, 92, 9118–9123. [Google Scholar] [CrossRef] [PubMed]
  102. Skogerboe, R.K.; Grant, C.L. Comments OH the Definitions of the Terms Sensitivity and Detection Limit. Spectrosc. Lett. 1970, 3, 215–220. [Google Scholar] [CrossRef]
  103. Cantwell, M.T.; Porter, S.E.G.; Rutan, S.C. Evaluation of the multivariate selectivity of multi-way liquid chromatography methods. J. Chemom. 2007, 21, 335–345. [Google Scholar] [CrossRef]
  104. Currie, L.A. Limits for qualitative detection and quantitative determination. Application to radiochemistry. Anal. Chem. 1968, 40, 586–593. [Google Scholar] [CrossRef]
  105. Hubaux, A.; Vos, G. Decision and detection limits for calibration curves. Anal. Chem. 1970, 42, 849–855. [Google Scholar] [CrossRef]
  106. Currie, L.A. Detection and quantification limits: Origins and historical overview. Anal. Chim. Acta 1999, 391, 127–134. [Google Scholar] [CrossRef]
  107. Clayton, C.A.; Hines, J.W.; Elkins, P.D. Detection limits with specified assurance probabilities. Anal. Chem. 1987, 59, 2506–2514. [Google Scholar] [CrossRef]
  108. Pellegrino Vidal, R.B.; Olivieri, A.C.; Tauler, R. Quantifying the prediction error in analytical multivariate curve resolution studies of multicomponent systems. Anal. Chem. 2018, 90, 7040–7047. [Google Scholar] [CrossRef]
  109. Olivieri, A.C. Second-order multivariate calibration with the extended bilinear model: Effect of initialization, constraints, and composition of the calibration set on the extent of rotational ambiguity. J. Chemom. 2020, 34, e3130. [Google Scholar] [CrossRef]
  110. Olivieri, A.C.; Faber, N.M. New developments for the sensitivity estimation in four-way calibration with the quadrilinear parallel factor moder. Anal. Chem. 2012, 84, 186–193. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Possible data arrays that can be built with second-order data prior to chemometric modelling.
Figure 1. Possible data arrays that can be built with second-order data prior to chemometric modelling.
Molecules 26 06357 g001
Figure 2. Possible data arrays that can be built with third-order data prior to chemometric modelling.
Figure 2. Possible data arrays that can be built with third-order data prior to chemometric modelling.
Molecules 26 06357 g002
Figure 3. Flow-chart showing different modelling approaches to analyse LC multi-way data with quantitative aims. All the models are indicated as examples. APARAFAC* is implemented only for third-order data analysis.
Figure 3. Flow-chart showing different modelling approaches to analyse LC multi-way data with quantitative aims. All the models are indicated as examples. APARAFAC* is implemented only for third-order data analysis.
Molecules 26 06357 g003
Table 1. Reports from 2018 to date based on chromatographic second-order/three-way data for quantitation purposes.
Table 1. Reports from 2018 to date based on chromatographic second-order/three-way data for quantitation purposes.
Analytes and SamplesModelRemarksRef.
LC-DAD
Malic, oxalic, formic, lactic, acetic, citric, pyruvic, succinic, tartaric, propionic and α-cetoglutaric acids in yoghurt, cultured milk, cheese and winePARAFAC
U-PLS/RBL
ET 10 min (isocratic mode).
LODs: 0.15–10.0 mmol L−1 in validation samples
[56]
Meloxicam, flurbiprofen, phenylbutazone, ibuprofen, diclofenac, mefenamic acid, celecoxib, naproxen, ketoprofen, diflunisal (non-steroidal anti-inflammatories) in Chinese patent drugs and health productsATLDET 14.5 min (gradient elution).
To simplify data processing, the retention time mode was subdivided in four regions.
LODs: 0.01–0.12 μg mL−1
[57]
Chlorogenic acid, (−)-epicatechin, caffeic acid, taxifolin, p-coumaric acid, hesperetin, naringenin, chrysin, apigenin, kaempferol, luteolin, quercetin, myricetin, rutin, (+)-catechin, ferulic acid, isorhamnetin (polyphenols) in raw propolisATLDET 16.5 min (gradient elution).
To simplify data processing, the elution time mode was subdivided in eight regions.
LODs: 0.01–0.38 μg mL−1
[58]
Gliclazide, glibenclamide, glimepiride (antidiabetics), atenolol, enalapril, amlodipine (antihypertensives) in serumMCR-ALS
U-PLS/RBL
ET 3 min (isocratic mode).
Elution time mode was subdivided in two regions.
LODs: <30 ng mL−1; better for U-PLS/RBL
[59]
Tacrolimus, everolimus, cyclosporine A (immunosuppressants) in whole bloodMCR-ALSET 2.7 min (isocratic mode).
Minimum sample preparation steps. The time mode was subdivided in three regions. Sample-added calibration strategy for matrix effect.
LODs: 0.56 μg L−1 (tracolimus), 0.08 µg L−1 (everolimus), 7.6 µg L−1 (cyclosporine A)
[60]
Methylparaben, ethylparaben, propyl-paraben, butylparaben, phenoxyethanol salicylic acid, methylisothiazolinone, 3-iodo-2-propynyl-n-butylcarbamate (preservatives) in facial masksATLD
MCR-ALS
ET 8.2 min (gradient elution).
Elution time mode was subdivided in four regions. Satisfactory and statistically comparable results were obtained with both models, except for phenoxyethanol in one studied sample (this fact was attributed to matrix interferences).
LODs: 1.2 ng mL−1 (butylparaben), 1466 ng mL−1 (3-iodo-2-propynyl-n-butylcarbamate)
[61]
Prednisolone, methylprednisolone (corticosteroids), mycophenolic acid (immunosuppressant) in human plasmaMCR-ALSET < 4 min.
Two isocratic elution methods with two different mobile phases were used.
Data array was divided into two regions (matrix augmentation in spectra and retention time direction were implemented for the first and second regions, respectively).
LODs: 0.9 µg L−1 (prednisolone), 1.3 µg L−1 (methylprednisolone), 300 µg L−1 (mycophenolic acid)
[62]
1,2-Dinitrobenzene, 1,3-dinitrobenzene, 2,4,6-trinitrotoluene, 2,4-dinitrotoluene, 2-nitrotoluene, 3-nitrotoluene and 4-nitrotoluene (explosives, agrochemical,
textile dyes and chemical intermediates)
in river and pond waters
MCR-ALSET 10 min (isocratic mode).
Very similar analyte structures.
LODs: 0.05–0.12 µg mL−1
[63]
Uric acid, creatinine, tyrosine, homovanillic acid, hippuric acid, indole-3-acetic acid, tryptophan, 2-methylhippuric acid (small molecules related to early diseases diagnosis) in human urineATLD
MCR-ALS
ET 6 min (isocratic mode).
Both models rendered comparable recoveries and root mean square error of predictions.
LODs: 29.9–464.2 ng mL−1 (ATLD); 11.7–127.1 ng mL−1 (MCR-ALS)
[64]
Chrysene, naphtalene, acenaphthylene, fluorene, phenanthrene, acenaphthene, anthracene, pyrene, benzo[a]anthracene, guaiazulene, benzo[e]pyrene, fluoranthene, benzo[a]pyrene, benzo[b]fluoranthene, benzo[k]fluoranthene (PAHs) in flue-dust and greasy dirt samplesATLDET 18 min (isocratic mode).
InertSustain®-C18 (5.0 μm, 4.6 mm × 250 mm) reversed phase column. Elution region was divided into four sub-segments.
LODs: 0.94–48.86 ng mL−1
[65]
Naproxen, ketoprofen, meloxicam
(non-steroidal anti-inflammatories) in Chinese patent drugs
ATLD-MCRATLD-MCR model was compared with PARAFAC, ATLD and MCR-ALS. Two simulated HPLC-DAD data sets, one simulated LC-MS data set, and a semi-simulated LC-MS data set were evaluated. ATLD-MCR proved to be able for handling chromatographic data with time shifts and signal overlapping. b[66]
2-Hydroxy-4-methoxybenzophenone (UV filter) in mice serum and human plasmaMCR-ALSET 2 min (isocratic mode).
Ultra-HPLC using two different experimental methods was applied. One of these methods rendered better results.
LOD: 0.66 ng mL−1
[67]
Sixteen PAHs (US EPA priority pollutants) in river waterMCR-ALSUSAEME as extraction method. Region-based analysis had improvement regarding to the whole data analysis.
LODs: 4.77–16.44 ng mL−1
[68]
Dorzolamide hydrochloride (carbonic
anhydrase-II inhibitor) and timolol maleate (non-specific adrenergic blocker) in an ophthalmic solution
PARAFAC
3W-PLS
U-PLS
ET 0.5 min (isocratic mode).
LODs: 0.57 µg mL−1 (dorzolamide hydrochloride), 0.66 µg mL−1 (timolol maleate)
[69]
Tartrazine, sunset yellow, carmine, amaranth, brilliant blue, aspartame, acesulfame potassium, sodium saccharin, caffeine, benzoic acid, sorbic acid, glycyrrhizin acid (food additives) in beverages (cola, grape, lemon, and orange sodas, green and black teas, orange and apple juices, milk drinks and grape wine)ATLDET 8.1 min (gradient elution).
ET mode was subdivided in four regions.
LODs: 1.40–165.1 ng mL−1
[70]
Gallic acid, epigallocatechin, epicatechin, epigallocatechin gallate, epicatechin gallate (polyphenols) in red, green, black and clinacanthus nutans Chinese teasATLD
MCR-ALS
ATLD-MCR
ET 5.4 min (isocratic mode).
Data array was divided in two sub-regions.
MCR-ALS and ATLD-MCR were better than ATLD in the case of larger time shifts
[71]
Carbendazim, thiabendazole, fuberidazole, carbofuran, carbaryl, flutriafol, 1-naphthol (pesticides) in vegetables (lettuce, cabbage leaf, carrot, beet, tomato, green bell pepper)MCR-ALSET 25 min (gradient elution).
Standard addition method due to matrix effect. Advantages of multi-way calibration in comparison with the univariate one when interferents are present (vegetable samples) were demonstrated
[72]
Epicatechin, myricetin, fisetin, quercetin, hesperidin, kaempferol, rutin (flavonoids) in raw and purified Chinese propolisATLDET 3.5 min (isocratic mode).
LODs: 0.01–0.20 μg mL−1
[73]
Melamine (plastic) in a food simulant migrated from kitchenwarePARAFAC
PARAFAC2
ET 2.2 min (isocratic mode).
Both models avoided the overestimation of migrated melamine amount despite coelution of the interferent with the analyte.
LOD: 0.58 mg L−1
[74]
Bisphenol A, bis(4-hydroxyphenyl) methane, 4,4′-cyclohexylidenebisphenol,
4,4′-(hexafluoroisopropylidene)diphenol, bis(4-hydroxyphenyl) sulfone [bisphenols] in methanol synthetic solutions
U-PLSET 4 min (isocratic mode). After quantitation, optimization of experimental conditions was made by inversion of PLS. LODs: 334–1156 μg L−1[75]
4,4′-isopropylidenediphenol, 4,4′-methylenediphenol, 4,4′-cyclohexylidene bisphenol, 4,4′-(hexafluoroisopropylidene) diphenol, 4,4′-sulfonyldiphenol [bisphenols] in a food simulantPARAFAC ET 4 min (isocratic mode).
Migration from BPA-free PC glasses is studied.
The maximum amount of BPA migrated from PC glasses (5.60 μg L−1) was lower than the established migration limit for a non-authorized substance.
LODs: 44.0–138.9 μg L−1.
[76]
Benzoic acid, methylisothiazolinone, sorbic acid, phenoxyethanol, methylparaben, ethylparaben, propylparaben, butylparaben, 3-iodo-2-propynyl-n-butylcarbamate, clorophene, triclocarban, triclosan (preservatives) in emulsion, cream, powder, gel and facial masksATLDET 10.5 min (gradient elution).
LC-DAD data were divided into five temporal regions for data processing.
LODs: 9.57 × 10−4–0.33 μg mL−1
[77]
Lenalidomide, gefitinib, crizotinib, chidamide, dasatinib, axitinib, lapatinib, erlotinib, nilotinib, idelalisib (anti-tumor drugs) in plasma, urine, cell culture mediumATLD, MCR-ALS ATLD-MCRET 6.5 min (gradient elution).
MCR-ALS and ATLD-MCR rendered better recoveries than ATLD.
LODs: 5.4–398 ng mL−1 (ATLD-MCR), 0.1–536 ng mL−1 (MCR-ALS)
[78]
Imidacloprid, albendazole, fenbendazole, praziquantel, fipronil, permethrin (veterinary active ingredients) in water from a wetland system used for the treatment of waste from a dog breeding plantMCR-ALSET 12 min (gradient elution).
Optimized DLLME. Data array was divided into six regions.
LODs: 1.3–8.5 ng mL−1
[79]
Gefitinib, crizotinib, dasatinib, axitinib, lapatinib, erlotinib, pexidartinib, nilotinib, LDK378 (tyrosine kinase inhibitors) in plasma, urine and a cell culture mediumATLD-MCRET 7 min (gradient elution).
Data were divided into three regions.
LODs: 0.02–0.24 μg mL−1 (plasma), 0.01–0.14 μg mL−1 (urine), 0.02–2.44 μg mL−1 (cell culture)
[80]
Sudan I, sudan II, sudan III, sudan IV, sudan red B, sudan red G, sudan red 7B, para red, diethyl yellow, methyl red, butter yellow, toluidine red (edible azo dyes) in chili sauces, saffron, ketchup, chili powderATLD
MCR-ALS
ET 6.5 min (isocratic mode).
Three-way data array was divided into four regions on the basis of elution time.
LODs: 0.01–2.56 mg kg−1 (ATLD), 0.01–2.95 mg kg−1 (MCR-ALS).
[80]
LC-FLD
Pipemidic acid, marbofloxacin, enoxacin,
ofloxacin, norfloxacin, ciprofloxacin, lomefloxacin, danofloxacin, enrofloxacin, sarafloxacin (quinolone antibiotics) in chicken liver, bovine liver and kidney
MCR-ALSET 4.7 min (isocratic mode).
LC-FLD data were divided into three temporal regions for data processing.
LODs: 7–125 μg kg−1
[81]
GC-MS
Butylated hydroxytoluene (BHT) [antioxidant], diisobutyl phthalate (DiBP), bis(2-ethylhexyl) adipate (DEHA), diisononyl phthalate (DiNP) [plasticizers], benzophenone (BP) [UV stabilizer] in Tenax aPARAFACET 19.1 min.
Data were acquired in SIM mode using five acquisition windows.
LODs 2.28 µg L−1 (BHT), 7.87 µg L−1 (DiBP), 3.04 µg L−1 (DEHA), 124.8 µg L−1 (DiNP), 10.57 µg L−1 (BP). Tenax could not be reused in this multiresidue determination
[82]
Butylated hydroxytoluene (BHT), benzophenone (BP), benzophenone-3 (BP3), diisobutyl phthalate (DiBP) [filters and additives] in sunscreen cosmetic creamsPARAFAC
PARAFAC2
ET 15.1 min.
Data were acquired in SIM mode using four acquisition windows.
LODs 7.93 µg L−1 (BHT), 12.40 µg L−1 (BP), 11.65 µg L−1 (DiBP), 279.8 µg L−1 (BP3). Analyte identification using the univariate standard method was incorrect. Multivariate calibration avoided false negative results
[83]
Butylated hydroxytoluene (BHT) [antioxidant], diisobutyl phthalate (DiBP), bis(2-ethylhexyl) adipate (DEHA), diisononyl phthalate (DiNP) [plasticizers] and benzophenone (BP) [UV stabilizer] in Tenax aPARAFAC
PARAFAC2
ET 19.1 min.
Migration from PE, PVC, and PP is studied. Data were acquired in SIM mode. Five acquisition windows were considered.
Presence of BHT, DiBP and DEHA was confirmed in Tenax blanks in some of the analysis. BP, DiBP migrated from both PVC film and PP coffee capsules, whereas DEHA migrated from PVC film
LODs: 3.48–360.2 µg L−1.
[84]
Bisphenol A (plasticizer) in a food simulant migrated from polycarbonate tableware, dichlobenil (pesticide) in onion, and oxybenzone (aromatic ketone) in sunscreen cosmetic creamsPARAFAC
PARAFAC2
Both models allowed analyte quantitation. The analyzed cases were: presence of interferents with overlapping peaks to the IS and analyte, coeluting compounds which share ions with the IS, retention time shifts from sample to sample, and coelution of interferents[74]
Butylated hydroxytoluene (BHT) [antioxidant], diisobutyl phthalate (DiBP) [plasticizer], benzophenone (BP) [UV stabilizer] in coffeePARAFACET 19.1 min.
Migration from plastic capsules is studied. SBSE for analyte extraction and concentration. Standard addition method due to matrix effect. Data acquired in SIM mode using three acquisition windows. Traces of the analytes found in the Milli-Q water samples were taken into account in the analysis. Found levels in coffee were below or around (DiBP case) than the migration established limits
[85]
Fluoranthene, benzo[b]fluoranthene,
chrysene, benzo[a]anthracene, pyrene (PAHs) in aerosol samples collected from Loudi City (China) in functional zones
ATLD-MCR
MCR-ALS
ET < 8 min.
Filters sample were extracted with the Soxhlet method. Scan mode was used for mass spectrum detection.
In real samples, ATLD-MCR
provided results which were better than or similar to MCR-ALS.
LODs: 0.003–0.087 µg mL−1.
[86]
LC-MS
Betamethasone, dexamethasone, triamcinolone acetonide, cortisone 21-acetate, dexamethasone 21-acetate, budesonide, triamcinolone acetonide acetate, fluocinonide, clobetasol 17-propionate, betamethasone dipropionate, beclomethasone dipropionate, beclomethasone, fluoromethalone,
fluticasone propionate, betamethasone 17-valerate (glucocorticoids) in face masks
ATLDET 11 min (gradient elution).
ESI interface operating in positive mode. LC–MS analysis in full scan mode. The three-way data array was divided into six sub-regions.
Betamethasone and dexamethasone (epimers) were simultaneously quantified under a simple elution program.
LODs: 0.56–13.55 ng mL−1.
[87]
Thiamine, riboflavin, nicotinic acid, biotin, nicotinamide, D-pantothenic acid, pyridoxine, folic acid, cyanocobalamin (B-group vitamins) in energy drinksATLD
APTLD
ET < 4.5 min (gradient elution).
ESI interface operating in positive mode. LC–MS analysis in full scan mode. Data array was divided into three sub-regions.
Both models rendered similar recovery and statistical results
LODs: 2 × 10−3–2.5 × 10−2 µg mL−1 (ATLD), 1 × 10−3–2.5 × 10−2 µg mL−1 (APTLD).
[88]
Estriol, 17α-estradiol, 17β-estradiol, estrone, ethinyl estradiol, diethylstilbestrol (estrogens), bisphenol A (xenoestrogen) in infant milk powderATLDET < 7 min (gradient elution).
ESI interface operating in negative mode.
LC–MS analysis in full scan mode. Data array was subdivided into four sub-regions on the basis of the elution ranges of estrogen.
LODs: 0.07–2.49 ng mL−1
[89]
Gallic acid, chlorogenic acid, caffeic acid,
(+)-catechin, p-coumaric acid, taxifolin, (−)-epicatechin, ferulic acid, myricetin, luteolin, quercetin (polyphenols) in Chinese propolis
ATLD
MCR-ALS
ET < 7.0 min (gradient elution).
ESI interface operating in negative mode. LC–MS analysis in full scan mode. Data array was subdivided into six sub-regions on the basis of the retention time.
LODs: 2.8–80.0 ng mL−1 (ATLD), 0.9–54.5 ng mL−1 (MCR-ALS)
[90]
Cyclosporine-A and tacrolimus (immunosuppressants) in blood and surface waterMCR-ALSESI interface operating in positive mode.
LC–MS analysis in full scan mode. The regions of interest method of the LC–MS data was employed for data compression. Matrix-matched calibration strategy was employed due to matrix effect.
LODs: (blood) 5.8 ng mL−1 (cyclosporine-A), 4.8 ng mL−1 (tacrolimus).
LODs (water) 2.3 × 10−2 ng mL−1 (cyclosporine-A), 9.0 × 10−2 ng mL−1 (tacrolimus)
[91]
17-β-estradiol, estrone, diethylstilbestrol (estrogens), bisphenol A (xenoestrogen) in river water [system I], L-glutamic acid, L-tyrosine, L-tryptophan, L-phenylalanine (amino acids), xanthine, hypoxanthine (purines), kynurenic acid, L-kynurenine (metabolites) in human urine [system II]ATLDCombination and partition of the MS1 full scan ion peaks recorded at different fragmentor voltages. Combined data and partitioned in two ways were compared using two systems.
System I: ET 4.4 min (gradient elution)
LODs: 0.18–2.72 ng mL−1(combined data), 0.25–2.35 ng mL−1 (partitioned data, ofv), 0.04–0.54 ng mL−1 (partitioned data, hfv).
System II: ET 4.4 min (gradient elution).
LODs: 0.99–5.43 ng mL−1 (combined data), 0.04–7.69 ng mL−1 (partitioned data, ofv), 2.96–7.38 ng mL−1 (partitioned data, hfv).
In most cases, data combination rendered higher sensitivity and more reliable results. Data partition provided higher selectivity in some cases but in others was unable to quantify analytes
[92]
a Tenax: food simulant for testing migration from plastics into dry foodstuffs. b For sample pre-processing and elution conditions in real samples see ref. [66]. Abbreviations: APTLD, alternating penalty trilinear decomposition; DLLME, dispersive liquid-liquid microextraction; ESI, electrospray ionization; ET, elution time; hfv, highest fragmentor voltage; IS, internal standard; LOD, limit of detection; ofv, optimum fragmentor voltage; PAH, polycyclic aromatic hydrocarbon, PE, polyethylene; PP, polypropylene; PVC, polyvinyl chloride; SBSE, stir bar sorptive extraction; SIM, single ion monitoring; USAEME, ultrasonic assisted emulsification micro-extraction; US EPA, United States Environmental Protection Agency; 3W-PLS, three-way partial least-squares.
Table 2. Reports from 2018 to date based on chromatographic third-order/four-way data for quantitation purposes.
Table 2. Reports from 2018 to date based on chromatographic third-order/four-way data for quantitation purposes.
Analytes and SamplesModelRemarksRef.
LC-EEFM
Rimsulfuron (herbicide), fuberidazole (fungicide), carbaryl (insecticide), naproxen (non-steroidal anti-inflammatory), albendazole (antihelminthic agent), tamoxifen (anticancer agent) in well and river watersMCR-ALSET < 12 min (gradient elution mode).
Native and photoinduced fluorescence (using a post-column UV reactor) were measured. Quadrilinearity was broken due to temporal shifts.
LODs: 0.02–0.27 ng mL−1 (spiked water samples) after a preconcentration SPE step
[93]
Benz[a]anthracene, chrysene, benzo[b]fluoranthene, benzo[a]pyrene (PAHs) in tea leavesMCR-ALSET ~ 9 min (isocratic mode).
Non-quadrilinear LC-EEM data type 4 were successfully processed with an MCR-ALS strategy. LODs: 1.0–1.4 ng mL−1 (validation samples).
LODs: 1.3–2.9 ng mL−1 (samples with interferents)
[28]
Pipemidic acid, marbofloxacin, enoxacin,
ofloxacin, norfloxacin, ciprofloxacin, lomefloxacin, danofloxacin, enrofloxacin, sarafloxacin (quinolone antibiotics) in animal tissues (chicken liver, bovine liver and kidney)
MCR-ALS
U-PLS/RTL
Third-order/four-way data results were compared with second-order/three-way data for the same system employing two different fluorescence detectors. MCR-ALS gave suitable results with second-order data but could not resolve all the analytes in the third-order/four-way system. U-PLS (with RBL or RTL) rendered good results in both cases, but the statistical indicators were not better than MCR-ALS second-order data. For more discussion, see section on Figures of Merit[94]
Pyridoxine (vitamin B6) in the presence of
L-tyrosine, L-tryptophan, 4-aminophenol
in synthetic aqueous samples
PARAFAC
APARAFAC
ET ~ 200 sec (isocratic mode).
Multilinearity was restored by chemometric processing.
LOD: 7 mg L−1
[95]
GC3/univariate detection
Citronellol, eugenol, farnesol, geraniol, menthol, trans-anethole, carvone, β-pinene (allergens) in perfumesPARAFACTwo detectors (FID and a mass analyzer) were used for data acquisition. GC3 system involved a first modulator (thermal desorption modulator, mp 6 s) that interfaced the first two columns, and a second modulator (differential flow modulator, mp 300 ms) which connected the last two columns. For data processing, smaller subsections of the chromatogram were used.
LODs (GC3-FID): 2.1–6.8 μL L−1,
LODs (GC3-MS): 4.8–8.5 μL L−1
[31]
Abbreviations: GC3, three-dimensional gas chromatography, ET, elution time; FID, flame ionization detector; LOD, limit of detection; mp, modulation period; PAH, polycyclic aromatic hydrocarbon; SPE solid phase extraction.
Table 3. AFOMs from second- and third-order calibration of a number of fluoroquinolone antibiotics using MCR-ALS. Data from ref. [94].
Table 3. AFOMs from second- and third-order calibration of a number of fluoroquinolone antibiotics using MCR-ALS. Data from ref. [94].
Analyte aSecond-OrderThird-Order
SEN bLOD cLOQ cSEN bLOD cLOQ c
PIPE1450.010.031500.10.3
MARBO370.070.2450- d- d
ENO90.10.3N.F. eN.F. eN.F. e
OFLO1800.040.121000.040.1
NOR2260.010.03340- d- d
CIPRO260.030.186- d- d
LOME480.030.16500.10.3
DANO2300.010.0340000.010.03
ENRO380.010.037200.080.2
SARA1200.020.069100.10.3
a PIPE, pipemidic acid, (MARBO), marbofloxacin, (ENO), enoxacin, (OFLO), ofloxacin, (NOR), norfloxacin, (CIPRO), ciprofloxacin, (LOME), lomefloxacin, (DANO), danofloxacin (ENRO), enrofloxacin, (SARA), sarafloxacin. b SEN = sensitivity in arbitrary fluorescence units × L mg–1. c LOD = limit of detection, LOQ = limit of quantitation, both in mg L–1. d LOD values are close to the maximum calibration concentrations. e N.F. = not found.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chiappini, F.A.; Alcaraz, M.R.; Escandar, G.M.; Goicoechea, H.C.; Olivieri, A.C. Chromatographic Applications in the Multi-Way Calibration Field. Molecules 2021, 26, 6357. https://doi.org/10.3390/molecules26216357

AMA Style

Chiappini FA, Alcaraz MR, Escandar GM, Goicoechea HC, Olivieri AC. Chromatographic Applications in the Multi-Way Calibration Field. Molecules. 2021; 26(21):6357. https://doi.org/10.3390/molecules26216357

Chicago/Turabian Style

Chiappini, Fabricio A., Mirta R. Alcaraz, Graciela M. Escandar, Héctor C. Goicoechea, and Alejandro C. Olivieri. 2021. "Chromatographic Applications in the Multi-Way Calibration Field" Molecules 26, no. 21: 6357. https://doi.org/10.3390/molecules26216357

APA Style

Chiappini, F. A., Alcaraz, M. R., Escandar, G. M., Goicoechea, H. C., & Olivieri, A. C. (2021). Chromatographic Applications in the Multi-Way Calibration Field. Molecules, 26(21), 6357. https://doi.org/10.3390/molecules26216357

Article Metrics

Back to TopTop