Detection of explosives on the surface of banknotes by Raman hyperspectral imaging and independent component analysis

The aim of this study was to develop a methodology using Raman hyperspectral imaging and chemometric methods for identification of pre- and post-blast explosive residues on banknote surfaces. The explosives studied were of military, commercial and propellant uses. After the acquisition of the hyperspectral imaging, independent component analysis (ICA) was applied to extract the pure spectra and the distribution of the corresponding image constituents. The performance of the methodology was evaluated by the explained variance and the lack of fit of the models, by comparing the ICA recovered spectra with the reference spectra using correlation coefficients and by the presence of rotational ambiguity in the ICA solutions. The methodology was applied to forensic samples to solve an automated teller machine explosion case. Independent component analysis proved to be a suitable method of resolving curves, achieving equivalent performance with the multivariate curve resolution with alternating least squares (MCR-ALS) method. At low concentrations, MCR-ALS presents some limitations, as it did not provide the correct solution. The detection limit of the methodology presented in this study was 50 μg cm⁻².     

Application of Multivariate curve resolution-alternating least square methods on the resolution of overlapping CE peaks from different separation conditions

Discussed in this paper is the development of a new strategy to improve resolution of overlapping CE peaks by using second-order multivariate curve resolution with alternating least square (second-order MCR-ALS) methods. Several kinds of organic reagents are added, respectively, in buffers and sets of overlapping peaks with different separations are obtained. Augmented matrix is formed by the corresponding matrices of the overlapping peaks and is then analyzed by the second-order MCR-ALS method in order to use all data information to improve the precision of the resolution. Similarity between the resolved unit spectrum and the true one is used to assess the quality of the solutions provided by the above method. 3,4-Dihydropyrimidin-2-one derivatives (DHPOs) are used as model components and mixed artificially in order to obtain overlapping peaks. Three different impurity levels, 100, 20, and 10% relative to the main component, are used. With this strategy, the concentration profiles and spectra of impurities, which are no more than 10% of the main component, can be resolved from the overlapping peaks without pure standards participant in the analysis. The effects of the changes in the components spectra in the buffer with different organic reagents on the resolution are also evaluated, which are slight and can thus be ignored in the analysis. Individual data matrices (two-way data) are also analyzed by using MCR-ALS and heuristic evolving latent projections (HELP) methods and their results are compared with those when MCR-ALS is applied to augmented data matrix (three-way data) analysis.     

Use of local rank‐based spatial information for resolution of spectroscopic images

Spectroscopic images are singular chemical measurements that enclose chemical and spatial information about samples. Resolution of spectroscopic images is focused on the recovery of the pure spectra and distribution maps of the image constituents from the sole raw spectroscopic measurement. In image resolution, constraints are generally limited to non‐negativity and the spatial information is generally not used. Local rank analysis methods have been adapted to describe the local spatial complexity of an image, providing specific pixel information. This local rank information combined with reference spectral information allows the identification of absent compounds in pixels with low compound overlap. The introduction of this information in the resolution process under the form of constraints helps to increase the performance of the resolution method and to decrease the ambiguity linked to the final solutions. Copyright © 2008 John Wiley & Sons, Ltd.     

Multivariate curve resolution-particle swarm optimization: A high-throughput approach to exploit pure information from multi-component hyphenated chromatographic signals

Multivariate curve resolution-particle swarm optimization (MCR-PSO) algorithm is proposed to exploit pure chromatographic and spectroscopic information from multi-component hyphenated chromatographic signals. This new MCR method is based on rotation of mathematically unique PCA solutions into the chemically meaningful MCR solutions. To obtain a proper rotation matrix, an objective function based on non-fulfillment of constraints is defined and is optimized using particle swarm optimization (PSO) algorithm. Initial values of rotation matrix are calculated using local rank analysis and heuristic evolving latent projection (HELP) method. The ability of MCR-PSO in resolving the chromatographic data is evaluated using simulated gas chromatography–mass spectrometry (GC–MS) and high-performance liquid chromatography–diode array detection (HPLC–DAD) data. To present a comprehensive study, different number of components and various levels of noise under proper constraints of non-negativity, unimodality and spectral normalization are considered. Calculation of the extent of rotational ambiguity in MCR solutions for different chromatographic systems using MCR-BANDS method showed that MCR-PSO solutions are always in the range of feasible solutions like true solutions. In addition, the performance of MCR-PSO is compared with other popular MCR methods of multivariate curve resolution-objective function minimization (MCR-FMIN) and multivariate curve resolution-alternating least squares (MCR-ALS). The results showed that MCR-PSO solutions are rather similar or better (in some cases) than other MCR methods in terms of statistical parameters. Finally MCR-PSO is successfully applied in the resolution of real GC–MS data. It should be pointed out that in addition to multivariate resolution of hyphenated chromatographic signals, MCR-PSO algorithm can be straightforwardly applied to other types of separation, spectroscopic and electrochemical data.     

Study of ε-caprolactone polymerization by NIR spectroscopy

Near-infrared (NIR) spectroscopy is proposed for the in-line quantitative and kinetic study of the polymerization of ε-caprolactone and eventually to facilitate real-time control of the manufacturing process. Spectra were acquired with a fibre-optic probe operating in transflectance mode immersed in the reactor. The NIR data acquired were processed using a multivariate curve resolution alternating least squares (MCR-ALS) algorithm. The proposed method allows calculation of the concentration and spectral profiles of the species involved in the reaction. The key point of this method is the lack of reference concentrations needed to perform the MCR-ALS method. The use of an extended spectral matrix using both process and pure analyte spectra solves the rank deficiency. The concentration profiles obtained were used to calculate a kinetic fitting of the reaction, but the method was improved by applying kinetic constraints (hard modelling). The rate constants of batches at different temperatures and the energy of activation for this reaction were calculated. Whenever possible, the hard modelling combined with the MCR-ALS method improves the fit of the experimental data: the results show good correlation between the NIR and reference data and allow the collection of high-quality kinetic information on the reaction (rate constants and energy of activation).     

A New Alternative Tool to Analyse Glycosylation in Monoclonal Antibodies Based on Drop-Coating Deposition Raman imaging: A Proof of Concept

Glycosylation is considered a critical quality attribute of therapeutic proteins as it affects their stability, bioactivity, and safety. Hence, the development of analytical methods able to characterize the composition and structure of glycoproteins is crucial. Existing methods are time consuming, expensive, and require significant sample preparation, which can alter the robustness of the analyses. In this context, we developed a fast, direct, and simple drop-coating deposition Raman imaging (DCDR) method combined with multivariate curve resolution alternating least square (MCR-ALS) to analyze glycosylation in monoclonal antibodies (mAbs). A database of hyperspectral Raman imaging data of glycoproteins was built, and the glycoproteins were characterized by LC-FLR-MS as a reference method to determine the composition in glycans and monosaccharides. The DCDR method was used and allowed the separation of excipient and protein by forming a “coffee ring”. MCR-ALS analysis was performed to visualize the distribution of the compounds in the drop and to extract the pure spectral components. Further, the strategy of SVD-truncation was used to select the number of components to resolve by MCR-ALS. Raman spectra were processed by support vector regression (SVR). SVR models showed good predictive performance in terms of RMSECV, R²_CV.     

Resolution and segmentation of hyperspectral biomedical images by multivariate curve resolution-alternating least squares

MCR-ALS is a resolution method that has been applied in many different fields, such as process analysis, environmental data and, recently, hyperspectral image analysis. In this context, the algorithm provides the distribution maps and the pure spectra of the image constituents from the sole information in the raw image measurement. Based on the distribution maps and spectra obtained, additional information can be easily derived, such as identification of constituents when libraries are available or quantitation within the image, expressed as constituent signal contribution. This work summarizes first the protocol followed for the resolution on two examples of kidney calculi, taken as representations of images with major and minor compounds, respectively. Image segmentation allows separating regions of images according to their pixel similarity and is also relevant in the biomedical field to differentiate healthy from non-healthy regions in tissues or to identify sample regions with distinct properties. Information on pixel similarity is enclosed not only in pixel spectra, but also in other smaller pixel representations, such as PCA scores. In this paper, we propose the use of MCR scores (concentration profiles) for segmentation purposes. K-means results obtained from different pixel representations of the data set are compared. The main advantages of the use of MCR scores are the interpretability of the class centroids and the compound-wise selection and preprocessing of the input information in the segmentation scheme.     

Flexible Implementation of the Trilinearity Constraint in Multivariate Curve Resolution Alternating Least Squares (MCR-ALS) of Chromatographic and Other Type of Data

Multivariate Curve Resolution Alternating Least Squares (MCR-ALS) can analyze three-way data under the assumption of a trilinear model using the trilinearity constraint. However, the rigid application of this constraint can produce unrealistic solutions in practice due to the inadequacy of the analyzed data to the characteristics and requirements of the trilinear model. Different methods for the relaxation of the trilinear model data requirements have been proposed, like in the PARAFAC2 and in the direct non-trilinear decomposition (DNTD) methods. In this work, the trilinearity constraint of MCR-ALS is adapted to different data scenarios where the profiles of all or some of the components of the system are shifted (not equally synchronized) or even change their shape among different slices in one of their data modes. This adaptation is especially useful in gas and liquid chromatography (GC and LC) and in Flow Injection Analysis (FIA) with multivariate spectroscopic detection. In a first data example, a synthetic LC-DAD dataset is built to investigate the possibilities of the proposed method to handle systematic changes (shifts) in the retention times of the elution profiles and the results are compared with those obtained using alternative methods like ATLD, PARAFAC, PARAFAC2 and DNTD. In a second data example, multiple wine samples were simultaneously analyzed by GC-MS where elution profiles presented large deviations (shifts) in their peak retention times, although they still preserve the same peak shape. Different modelling scenarios are tested and the results are also compared. Finally, in the third example, sample mixtures of acid compounds were analyzed by FIA under a pH gradient and monitored by UV spectroscopy and also examined by different chemometric methods using a different number of components. In this case, however, the departure of the trilinear model comes from the acid base speciation of the system depending on the pH more than from the shifting of the FIA diffusion profiles.     

Multivariate curve resolution applied to the simultaneous analysis of electrochemical and spectroscopic data: study of the Cd(II)/glutathione-fragment system by voltammetry and circular dichroism spectroscopy

Multivariate curve resolution with alternating least squares (MCR-ALS) is applied for the first time to the simultaneous analysis of electrochemical and spectroscopic data. Then, a data analysis is done with augmented matrices constituted by Differential Pulse Polarography and Circular Dichroism data submatrices. The use of proper, and different for each submatrix, constrains in the iterative ALS optimization allows to obtain chemically meaningful results constituted by a common matrix containing the concentration profiles, and two matrices with the pure electrochemical and spectroscopic signals. MCR-ALS is applied to the study of the complexation of Cd by Cys-Gly, a glutathione-fragment of great interest for understanding metal-phytochelatins complexation.     

Characterisation of hydrogen bond perturbations in aqueous systems using aquaphotomics and multivariate curve resolution-alternating least squares

Aquaphotomics is a new discipline that provides a framework for understanding changes in the structure of water caused by various perturbations, such as variations in temperature or the addition of solutes, using near infrared spectroscopy (NIRS). One of the main purposes of aquaphotomics is to identify water bands as main coordinates of future absorbance patterns to be used as biomarkers. These bands appear as consequence of perturbations in the NIR spectra. Curve resolution techniques may help to resolve and find new water bands or confirm already known bands. The aim of this study is to investigate the application of multivariate curve resolution-alternating least squares (MCR-ALS) to characterise the effects of various perturbations on the NIR spectra of water in terms of hydrogen bonding. For this purpose, the perturbations created by temperature change and the addition of four solutions of different ionic strength and Lewis acidity were studied (NaCl, KCl, MgCl₂ and AlCl₃, with concentrations ranging from 0.2 to 1 mol L⁻¹ in steps of 0.2 mol L⁻¹). Transmission spectra of all salt solutions and pure water were obtained at temperatures ranging from 28 to 45 °C. We have found that three distinct components with varying temperature dependence are present in water perturbed by temperature. The salt solutions studied exhibited similar trends with respect to the temperature perturbation, while the peak locations of their MCR-ALS pure components varied according to the ionic strength of the salt used.     

Infrared study of aging of edible oils by oxidative spectroscopic index and MCR-ALS chemometric method

One of the most suitable analytical techniques used for edible oil quality control is Fourier transform mid infrared spectroscopy (FT-MIR). FT-MIR spectroscopy was used to continuously characterize the aging of various edible oils thanks to a specific aging cell. There were differences in the spectra of fresh and aged oils from different vegetable sources, which provide the basis of a method to classify them according to the oxidative spectroscopic index value. The use of chemometric treatment such as multivariate curve resolution-alternative least square (MCR-ALS) made it possible to extract the spectra of main formed and degraded species. The concentration profiles gave interesting information about the ability of the various oils to support the oxidative treatment and showed that all oils present the same aging process. Both methods led to concordant results in terms of induction times determined by the oxidative spectroscopic index and the appearance of oxidation products revealed by MCR-ALS.     

Novel combination of non-aqueous capillary electrophoresis and multivariate curve resolution-alternating least squares to determine phenolic acids in virgin olive oil

This paper presents the development of a non-aqueous capillary electrophoresis method coupled to UV detection combined with multivariate curve resolution-alternating least-squares (MCR-ALS) to carry out the resolution and quantitation of a mixture of six phenolic acids in virgin olive oil samples. p-Coumaric, caffeic, ferulic, 3,4-dihydroxyphenylacetic, vanillic and 4-hydroxyphenilacetic acids have been the analytes under study. All of them present different absorption spectra and overlapped time profiles with the olive oil matrix interferences and between them. The modeling strategy involves the building of a single MCR-ALS model composed of matrices augmented in the temporal mode, namely spectra remain invariant while time profiles may change from sample to sample. So MCR-ALS was used to cope with the coeluting interferences, on accounting the second order advantage inherent to this algorithm which, in addition, is able to handle data sets deviating from trilinearity, like the data herein analyzed. The method was firstly applied to resolve standard mixtures of the analytes randomly prepared in 1-propanol and, secondly, in real virgin olive oil samples, getting recovery values near to 100% in all cases. The importance and novelty of this methodology relies on the combination of non-aqueous capillary electrophoresis second-order data and MCR-ALS algorithm which allows performing the resolution of these compounds simplifying the previous sample pretreatment stages.     

多元曲线分辨-交替最小二乘方法在联用数据分辨中的应用

通过对模拟数据和高效毛细管电泳实验数据的分析,讨论了多元曲线分辨-交替最小二乘方法(MCR-ALS)在毛细管电泳-二极管阵列检测(CE-DAD)联用数据分辨中的应用.讨论了几种因素对MCR-ALS单个数据矩阵分辨结果的影响,包括待分析物光谱间的相似程度、浓度曲线的重叠程度以及由渐进因子分析(EFA)所得到的浓度初始值等.MCR-ALS还可用于多个数据矩阵的同时分析,即二阶MCR-ALS.结果表明,与一阶MCR-ALS相比,二阶MCR-ALS方法能够更好地解决各种分辨问题,得到合理和满意的分辨结果.     

多元曲线分辨结合在线红外用于3-氨基-4-氨基肟基呋咱的合成反应机理

利用红外光谱在线监测丙二睛、亚硝酸钠和盐酸羟胺合成3-氨基-4-氨基肟基呋咱的反应过程,采用多元曲线分辨-交替最小二乘法（MCR-ALS）、直观推导式演进特征投影法(HELP)等化学计量学方法对反应过程所获得的实时红外光谱数据矩阵进行解析,得到了各组分纯物质的浓度变化曲线和对应的红外光谱,并将多元曲线分辨 交替最小二乘法与直观推导式演进特征投影法的分析结果进行比较,得出可相互验证的一致结论,据此推出该反应合理的反应机理。 2种方法得到的反应物与生成物的光谱与原光谱的相似度近似于1,说明该解析方法具有准确性和可靠性。 结果表明,化学计量学结合红外光谱可有效的应用于3-氨基-4-氨基肟基呋咱合成过程的机理推断。     

Multivariate curve resolution of pH gradient flow injection mixture analysis with correction of the Schlieren effect

Multivariate curve resolution using alternating least squares (MCR-ALS) was used to quantify ascorbic (AA) and acetylsalicylic (ASA) acids in four pharmaceutical samples using a flow injection analysis (FIA) system with pH gradient and a diode array (DAD) spectrometer as a detector. Four different pharmaceutical drugs were analyzed, giving a data array of dimensions 51 x 291 x 61, corresponding respectively to number of samples, FIA times and spectral wavelengths. MCR-ALS was applied to these large data sets using different constraints to have optimal resolution and optimal quantitative estimations of the two analytes (AA and ASA). Since both analytes give an acid-basic pair of species contributing to the UV recorded signal, at least four components sholuld be proposed to model AA and ASA in synthetic mixture samples. Moreover, one additional component was needed to resolve accurately the Schlieren effect and another additional component was also needed to model the presence of possible interferences (like caffeine) in the commercial drugs tablets, giving therefore a total number of 6 independent components needed. The best quantification relative errors were around 2% compared to the reference values obtained by HPLC and by the oxidation-reduction titrimetric method, for ASA and AA respectively. In this work, the application of MCR-ALS allowed for the first time the full resolution of the FIA diffusion profile due to the Schlieren effect as an independent signal contribution, suggesting that the proposed MCR-ALS method allows for its accurate correction in FIA-DAD systems.     

Multivariate resolution of rank-deficient near-infrared spectroscopy data from the reaction of curing epoxy resins using the rank augmentation strategy and multivariate curve resolution alternating least squares approach

A model of the curing reaction between phenyl glycidyl ether (PGE) and aniline as the curing agent was studied isothermally at 95 °C and monitored in situ by near-infrared spectroscopy (NIR). The spectra were recorded every 5 min. The ubiquitous problem of rank deficiency in reaction network systems was solved by assembling an augmented column-wise matrix containing five process runs from different initial conditions. The data were analyzed using a two-way multivariate curve resolution alternating least squares method (MCR-ALS). Initial estimates of spectra required by MCR-ALS were given by a SIMPLe-to-use Interactive Self-modeling Mixture Analysis (SIMPLISMA) approach. The reactants, product and intermediate spectra were successfully resolved and the concentration profiles properly represented the system studied. The performance of the model was evaluated by two parameters: ALS lack of fit (lof=0.88%) and explained variance (R²=99.99%). To validate the MCR-ALS results, the similarity coefficients (r) between the recovered spectra and the pure species spectra were calculated. These were: PGE (r=0.998), aniline (r=0.994) and tertiary amine (r=0.999).     

直观推导式演进特征投影法对酶催化反应的过程分析

以丁香酸(syringic acid)为模型化合物, 研究了漆酶(laccase)催化降解木质素(lignin)的复杂生化反应过程, 设计了一种反应过程动态量测系统, 该系统可以在5 s的时间间隔测定反应体系在190～800 nm波长范围内的实时光谱信息. 利用固定窗口因子分析-直观推导式演进特征投影法(FSMWEFA-HELP), 解析反应过程中测得的实时动力学光谱数据矩阵, 得到反应物和中间体的数目及其浓度的变化和纯物质的光谱曲线, 并推导出合理的反应机理. 将所得结果与多元曲线分辨-交替最小二乘(MCR-ALS)方法进行了比较.     

Spectroscopic and quantitative analysis of spiroorthoester synthesis by two-dimensional correlation and multivariate curve resolution methods of NIR data

The spiroorthoester synthesis includes several competitives reactions. A way of determining the reactions that are taking place and their sequential order, is presented. The reaction between the phenylglycidylether and gamma-butyrolactone to obtain a spiroorthoester has been monitored by near-infrared spectroscopy (NIR). In addition to the formation of the corresponding spiroorthoester, some parallel processes can occur. By means of two-dimensional correlation analysis, only one reaction is postulated, the one corresponding to the spiroorthoester formation. This was confirmed by recording the NMR spectra of the final product. Applying multivariate curve resolution-alternating least squares (MCR-ALS) to the NIR spectra obtained during the reaction, it has been possible to obtain the concentration values of the species involved in the reaction. The recovered spectra were compared with the experimentally recorded spectra for the reagents (phenylglycidylether, gamma-butyrolactone) and the final product (spiroorthoester) and the correlation coefficients were, in all cases, higher than 0.990. The maximum and minimum limits associated with the ALS solutions were calculated, making it possible to limit to a considerable extent the ambiguity that is characteristic of these curve resolution methods.     

Multivariate curve resolution applied to Fourier transform infrared spectra of macromolecules: structural characterisation of the acid form and the salt form of humic acids in interaction with lead

The significance of evolving mixtures structural spectroscopic studies might appear limited when the experimental spectra do not present a sufficient quality for a precise interpretation. It is the case when the chemical behaviour of macromolecules is studied on the basis of infrared spectra. If the effective resolution is low, the spectral profiles appear similar despite the applied chemical conditions change. This makes impossible the interpretation of the raw spectra and mathematical treatments are required to separate the different contributions that overlap.To determine the behaviour of the reactive sites of humic acids in the binding with heavy metals, infrared spectra are recorded under various chemical conditions. The cation to be considered is Pb²⁺ and the two chemical variables to be studied are pH and initial lead concentration. Four series of FTIR spectra are recorded, but no visible difference can be directly assigned to the different chemical states of the macromolecules. Multivariate self-modelling curve resolution is thus proposed as a tool for resolving these complex and strong overlapping datasets. First, initial estimates are obtained from pure variable detection methods: it comes out that two spectra are enough to reconstruct the experimental matrices. In a further step, the application of the multivariate curve resolution-alternating least squares (MCR-ALS) algorithm with additional constraints on each individual dataset, as well as on column-wise augmented matrices, allows to optimise the profiles and spectra that appear to be highly characterising the acid and the salt form of the molecule. Moreover, the concentrations profiles associated to these two limit spectral forms allow interpreting the analytical measurements made during the reactions between humic acids and H⁺ or Pb²⁺. Consequently, depending on the initial state of the humic acid, two distinct reactional mechanisms are proposed.