Initial evaluation of the functionalized polymer Raman data

The Raman spectra measured with the i-Raman EX spectrometer of the low, medium and highly functionalized polymer exhibits excellent specificity in the 900-1300 cm^-1 region (see Figure 2). An additional feature observed at 981 cm^-1 derives from an inorganic species. A second derivative spectral pretreatment provides a well-defined baseline and resolution of the desired aromatic peaks for the initial and functionalized polymer.

A simple peak height ratio analysis is used to evaluate the quality of the Raman data to determine polymer functionalization. Twelve samples that varied in the extent of polymer functionalization as determined by ¹H NMR spectroscopy were selected. Next the single point peak height intensities of the second derivative spectra at 1126 and 1005 cm^-1 are determined. The peak height analysis of the two components (i = initial, f = final) utilized the ratio of these two unique bands for component band f divided by the sum of the two component bands (i + f). The expected range of this ratio will be 0 to 1, where 1 corresponds to a 100% functionalized polymer. This simple univariate band ratio analysis resulted in very good linearity (see Figure 2) confirming the suitability of the Raman spectral data to develop a robust quantitative PLS1 model.

Standard selection

A well-designed calibration set includes the expected possible variation of future samples. This includes concentration ranges of constituents, as well as process-induced and environmental variability. For the PCA-MD classification of the organic and aqueous phase samples, this included multiple production samples from both the organic and aqueous phases. For the final PLS1 quantitative analysis this included over seventy samples from initial research lab design-of-experiments, pilot plant batches, and production lots that were fully characterized by the ¹H NMR spectroscopy.

PCA-MD Discrimination analysis of organic and aqueous phases

Principal component analysis (PCA) is a chemometric method in which the original spectral data is transformed into a new coordinate system, which condenses the information found in the original inter-correlated variables into a few uncorrelated variables, called principal components. The PCA can be used to reduce the dimensionality by zeroing out one or more of the smallest principal components. This results in a lower-dimensional projection of the data that preserves the maximal data variance.

This multivariate analysis reveals the internal structure of the data in a way that enables a simple explanation of the variance in the data. As such, it is a useful tool for visualization of datasets. In any principal component rotation, the first component accounts for the maximum proportion of the variance and subsequent components account for the remaining variance.

PCA clearly enables grouping or clustering the Raman spectra into different classes. The cumulative explained ratio plot quantifies how much variance is explained in the first N components. In this application, three (N) variables explain 98% of the variance. Cluster analysis methods aggregate the Raman spectra stepwise, based upon the similarity of their spectral features (variables). A simple distance measurement of the individual spectrum relative to the various clusters provides the pattern recognition. A standard distance metric includes the Mahalanobis distance. Here, the Mahalanobis distance is a statistical tool that is coupled to PCA and can be considered as the distance of a point from the center of a principal component cluster. The smaller the Mahalanobis distance, the higher probability that the point belongs to the cluster.

A simple PCA-MD classification model was successfully developed using 37 calibration and 34 validation Raman spectra. The model was 100% accurate in classification of the training set Raman spectra for both the organic and aqueous phases.

PLS1-based quantitative determination of % functionalized polymer

The PLS regression is one of the most popular multivariate techniques presently used in vibrational spectroscopic methods for quantitative analyses. It is used to develop a linear model between two matrices, the Raman spectral data and the reference values. The PLS model finds the empirical relationship between the two matrices (X and Y) and enables prediction of a continuous property, such as the % polymer functionalization, from the training set. In the case of the Raman spectra of the standards, the variation in the 1002 cm^-1 and the 1127 cm^-1 spectral features are highly correlated with the % polymer functionalization.

The developed quantitative model for the % functionalized polymer employs a PLS1-based (partial least squares for a single variable) algorithm. The Raman spectrum uses the 995-1200 cm^-1 spectral region which includes diagnostic bands from both the starting and reacted polymer. Data pretreatment included a second derivative. The overall performance of the final quantitative model is summarized in the plot of the predicted (PLS) versus known (¹H NMR) % functionalized polymer (see Figure 3). 29 calibration (blue symbols) and 31 validation (red symbols) standards where used. Excellent linearity was observed with an R² of 0.95 for the calibration standards and 0.87 for the validation standards. The RMSE (root mean square error) is based upon the residuals of the calibration data and provides a simple measure of the predicted values. This was calculated for both the calibration (RMSEC) and validation (RMSEP) standards, with an RMSEC of 1.22 and an RMSEP of 1.43. The RME less than two (<2) demonstrates the model exhibits good predictive performance.

PLS1 quantitative Raman model accuracy and precision

The final developed model provides good linearity between the predicted (PLS) and assayed (¹H NMR) values for both the calibration and validation standards. Table 3 summarizes the precision and accuracy performance of the quantitative PLS Raman model. This includes the RMSEP from the validation standards (see Figure 4) a standard single-day and lot precision measurement and finally the average RMSE calculated for 16 subsequent manufactured lots (see Figure 4). The observed precision of 0.49% and the RMSE of 1.31 to 1.43 demonstrates the good statistical performance of this quantitative Raman method.

High-quality Raman spectra of an aqueous based water-soluble polymer (10-20% solids) are rapidly and easily measured using the i-Raman EX spectrometer. The Raman spectrum provides strong, unique bands for both the initial and fully reacted polymer. This enables development of a simple, robust quantitative analysis of the percent polymer functionalization using a PLS1 chemometric model with the vendor-supplied software (BWIQ). The final Raman method exhibits excellent linearity, accuracy and precision. This Raman method is now used routinely in a multiuser plant quality control (QC) laboratory.

We gratefully acknowledge Ian Anderson, Michael Lynch, Michele Astrove, Boyu Li, and Peter J. Larkin of Solvay Technology Solutions for providing the resources needed to carry out this study and for permission to publish this work. We also acknowledge Katherine Bakeev (B&W Tek) for useful discussions about classification methodologies.