Conference paper

ALS Scheme using Extent-based Constraints for the analysis of Chemical Reaction Systems

Multivariate curve resolution via alternating least squares (ALS) is used to resolve the concentration profiles C and the pure component spectra E of S species from the multivariate absorbance data A, assuming the bilinear model A = C E. Due to the possible permutations of profiles and the presence of intensity and rotational ambiguities, soft constraints such as nonnegativity of C and C as well as unimodality, monotonicity, closure, and local rank selectivity of C are typically used to obtain tighter solution bounds for C and E [1].

In addition, hard constraints in the form of kinetic models are also often used. Unfortunately, these models are subject to structural plant-model mismatch and parametric uncertainty, which weakens their impact. As an alternative, this paper proposes to use constraints based on variant states called extents. The computation of these extents does not require any information on the rate processes, that is, x(t) = T n(t), with x(t) and n(t) the vectors of extents and numbers of moles at time t, respectively, and T a matrix known from the reaction stoichiometry, the inlet composition and the initial conditions [2]. Expressing the S concentrations in terms of d extents and q = S - d invariants reduces the dimensionality of the problem from S to d. Each column of the extent matrix X describes the extent of a single rate process, for example of a reaction or an inlet flow. It turns out that the unknown concentration matrix C can be expressed in terms of the lower-dimensional matrix X as C = V -1 X T -T, where V is a diagonal matrix containing the volume profile. Since each column of X describes a single rate process, additional constraints can be enforced on X such as monotonicity and convexity/concavity [3]. Furthermore, the q invariant relationships can be used as constraints in the least-squares problem. As a consequence, the use of extents in ALS reduces the ambiguity between C and E, yielding faster convergence and tighter solutions.

The use of extent-based constraints also opens up new perspectives for hard-soft ALS methods, since hard kinetic models can be identified individually (that is independently of the other rates) for some selected processes, while soft extent-based constraints are used for the unknown processes. Another feature involves the possibility of initializing or constraining the ALS scheme with a concentration submatrix (of dimension at least S × S) estimated from multiple experiments performed under well-designed conditions and local-rank information. Together with the corresponding absorbance data, this submatrix can efficiently replace the traditional initialization via factor analysis and be used to compute a better initial estimate of E.

After a brief review of the mathematical properties of extents for batch and open reactors, this talk will present the modified ALS scheme and illustrate it via simulated examples.

[1] A. de Juan, J. Jaumot, R. Tauler, Anal. Methods 6, 4964 (2014).
[2] D. Rodrigues, S. Srinivasan, J. Billeter, D. Bonvin, Comp. Chem. Eng. 73, 23 (2015)
[3] S. Srinivasan, D.M.D. Kumar, J. Billeter, S. Narasimhan, D. Bonvin, IFAC Symposium DYCOPS (2016)

Related material