Chemometrics and Intelligent Laboratory Systems 126 (2013) 91–99 Contents lists available at SciVerse ScienceDirect Chemometrics and Intelligent Laboratory Systems journal homepage: www.elsevier.com/locate/chemolab How to correct inner filter effects altering 3D fluorescence spectra by using a mirrored cell X. Luciani ⁎, R. Redon, S. Mounier Université de Toulon, PROTEE, EA 3819, 83957 La Garde, France a r t i c l e i n f o Article history: Received 7 February 2013 Received in revised form 17 April 2013 Accepted 22 April 2013 Available online 2 May 2013 Keywords: EEM Inner filter effects PARAFAC Fluorescence Mirrored cell Multi-way a b s t r a c t In this paper we present a new correction method of inner filter effects that occurs when measuring fluorescence Excitation–Emission Matrices (EEM) of concentrated solutions. While traditional method requires absorption measurement or sample dilution(s), the Mirrored Cell Approach (MCA) only requires two different EEM of the considered sample: a first one using a traditional cell and a second one using a mirrored cell. The mathematical relationship between both models is originally exploited to obtain a simple numerical correction. Method is validated using a set of known mixtures. In addition we show that advanced multilinear analysis can be efficiently applied on to the corrected EEM. © 2013 Elsevier B.V. All rights reserved. 1. Introduction Standard fluorometers allow to measure the fluorescence intensity emitted by a solution for a given couple of excitation and emission wavelengths. A fluorescence Excitation Emission Matrix (EEM) (also called 3D spectrum) is then obtained by scanning both wavelength domains [1] and allows to characterize the solution. EEMs are now widely used in various scientific domains such as medicine [2], analytical chemistry [3] or environmental sciences, in particular for Dissolved Organic Matter (DOM) tracing and characterization purpose [4,5]. Measured signal relies on the spectroscopic properties of each fluorescent component of the solution (fluorophore) and their concentration in the solution. EEM analysis consists of deducing these underlying parameters from the measured spectrum and thus characterizing mixture components and possibly estimating their relative contribution to the sample signal. A significant panel of chemometric tools has been adapted or specifically developed to solve this inverse problem, from basic peak picking [6] to most advanced approaches based on multidimensional algebra [7]. All of these methods resort to a model of the measured fluorescence signal. Classical multilinear model assumes that the EEM of a single fluorophore, has a magnitude proportional to the fluorophore concentration in the solution and that its pattern is given by the outer product between the excitation spectrum and the emission ⁎ Corresponding author at: Laboratoire PROTEE, Batiment R, Université du Sud ToulonVar, BP 20132, 83957 La Garde Cedex, France. Tel.: +33 4 94 14 23 18. E-mail addresses: [email protected] (X. Luciani), [email protected] (R. Redon), [email protected] (S. Mounier). 0169-7439/$ – see front matter © 2013 Elsevier B.V. All rights reserved. http://dx.doi.org/10.1016/j.chemolab.2013.04.014 spectrum of the fluorophore. Hence according to this model, EEM of a mixture of fluorophores is described by a linear combination of the individual EEM of each fluorophore. This model is very useful in practice, notably when considering several mixtures (samples) of the same fluorophores. Indeed, the corresponding excitation–emission matrices set defines a three-way array and the linear model of fluorescence is nothing else than the Canonical Polyadic decomposition (CP) of this array [8], also known in the community as CANonical DECOMPosition (CANDECOMP) [9] or PARAllel FACtor analysis (PARAFAC) [10,7]. We keep here the acronym CP that honors the original work and also stands for CANDECOMP/PARAFAC. This multilinear decomposition usually admits a unique solution [11] and many algorithms allow to perform the decomposition [12–14]. Thereby it has rapidly become one of the most popular EEM analysis method. However it is well known that the pertinence of the linear model decreases with the solution absorbancy [1] meaning that the gradual absorption by the solution of both exciting and fluorescent lights cannot be neglected. These effects are known as inner filter effects (IFE) and affect both EEMs magnitude and pattern. They are observed and studied for a long time now [15,16] and it had been shown that IFE are still perceptible at low absorption and occur in practical applications such as DOM analysis [17,18]. The result is that traditional linear excitation–emission matrices analysis methods cannot be directly applied on the raw measured EEMs. Two main methods are currently used to restore EEM multilinear properties (we then speak of corrected EEM). The first one is to strongly dilute the solution so that absorption of the diluted solution is lower than a certain threshold [19]. This threshold is hard to define precisely so that it is usually preferable to realize dilution series and compare the corresponding 92 X. Luciani et al. / Chemometrics and Intelligent Laboratory Systems 126 (2013) 91–99 measured EEM to effectively unsure the linearity. Another drawback of this method is that contamination or physico-chemical changes can occur during the dilution, thus modifying the fluorescence properties of the sample. The most common alternative [20,21], namely the Absorption Correction Approach (ACA) resorts to mathematical model of inner filter effects derived from the Beer–Lambert law [22–24]: −ðAðλex ÞþAðλem ÞÞ 2 F ðλex ; λem Þ ¼ Lðλex ; λem Þ10 ; ð1Þ where F is the measured EEM, L the corrected one and A is the solution absorption. Rigorous physical justification of this mathematical model can be found in [25] and [18]. ACA uses the measured absorption spectrum as an estimation of A and then deduces L from F using Eq. (1). Similar methods were proposed in [25–27]. However ACA requires another experimental device. In addition the measured absorption spectrum is sometimes a poor estimate of A because in case of highly absorbing solutions it is obtained from the ratio of two weak signals. In order to avoid these limitations a mathematical correction method has been proposed in [28]. More recently we have introduced the Controlled Dilution Approach (CDA) [18]. CDA does not require absorption measurement but uses a second EEM of the sample to correct the first one. This second EEM is measured from a controlled (and reduced) dilution of the considered sample. The dilution is not intended to suppress IFE by itself so that the dilution factor can be chosen arbitrarily small in order to avoid the drawback of the dilution approach. The key point of CDA is that this second EEM can be modelized by a second equation involving F and L. Then it is shown that the linear term L can be directly deduced from the combination of both equations. Hence CDA still involves additional experimental modifications of the sample. Therefore we introduce in this paper a Mirrored Cell Approach (MCA) which is based on a similar idea but suppresses the dilution step. More precisely, instead of measuring the second EEM from a dilute sample, this one is obtained from the same sample but put into a mirrored cell. MCA principle and implementation are described in the next section. We first derive some new mathematical models of IFE that take into account main multiple reflections into the mirrored cell. Then we describe MCA as an EEM correction method. In the third section we show how to set MCA parameters and correct practical issues. Third and fourth sections are dedicated to the experimental validation. MCA ability for IFE correction is studied onto a set of various mixtures of three fluorophores, along with its effectiveness as a CP pretreatment of non-linear sets of EEMs. We then conclude about the reliability of this new approach. 2. Theory 2.1. Physical models of inner filter effects We derive here from physical considerations a new mathematical model of inner filter effects. The proposed model takes into account multiple reflections into the cell. Although several simplifications (described below) are still assumed, it makes this model suitable for mirrored cell. We extend here the reasoning used in [18]. We consider a mixture of N fluorophores. For each fluorophore n(n = 1 ⋯ N), cn, εn(λex), εn(λem), Φn, γn(λem) denotes the concentration in the solution, the molar extinction coefficient at the excitation wavelength λex, the molar extinction coefficient at the emission wavelength λem, the quantum yield and the emission probability at wavelength λem respectively. Absorption coefficient is defined by αn(λ) = cnεn(λ). In the following model, we also take into consideration possible presence of chromophores in the solution, i.e.: species that are absorbing light but not fluorescing. Hence we denote α0 the total absorption coefficient of chromophores and thus the total absorption coefficient α of the solution is given by: α = ∑nN = 0 αn. Fig. 1 recalls basically the experimental device of right angle standard fluorometers, with multiple reflection on the facets. 0 0 Fig. 1. Scheme of the mirrored cell and considered optical paths, view from above. In order to simplify model equations, one classically considers perfectly collimated light beams in both excitation and emission. Although, this is usually not the case in practice (notably with our system) this simplification is mathematically justified by small angle approximation. The excitation light I0(λex) is absorbed through the sample cell (length l) by the solution (primary inner filter effect). This induces the fluorescent light which is also partly absorbed by the solution (secondary inner filter effect). A fraction Fmir(λem) of this fluorescence signal is then collected perpendicularly to the exciting beam. We make the following main approximations. First of all, fraction of the exciting light which do not reach the “influence zone” Z in Fig. 1 is neglected as well as the fluorescence light issued from the region outside Z. In other words two main optical paths are considered which represent the excitation beam and emission beam in Fig. 1 scheme. Furthermore, diffusion and re-emission effects are neglected. A model of inner filter effects for mirrored cells has been proposed in [29] that take into account reflections due to the mirrored facet of the cell. In the present study we assume that the reflection coefficients of both the mirrored and non-mirrored facets are non-null. These are denoted Rmir and R respectively. Hence we consider the multiple reflexions depicted in Fig. 1 scheme. However variations of Rmir and R according to the wavelength are neglected. Finally, wavelength dependence of I0(λex) is corrected by the apparatus so that I0(λex) = I0. The cell is then divided into horizontal and vertical elementary strips of respective dimensions dy × l and l × dx. At the entry of the cell (x = 0) each horizontal strip is supposed to receive an equal elementary fraction of the exciting light from the rectilinear exciting beam: dyI Δ . Beer–Lambert law then quantifies the intensity transmitted to a position x in the cell, where a fraction αndx of this intensity is absorbed by the fluorophore n. Hence, taking into account the primary inner filter effect (α(λex)) and the multiple reflections that depend on R and Rmir values, the total intensity absorbed by fluorophore n in an elementary horizontal strip of the “influence zone” (An) is given by: 0 y An ðλex Þ ¼ ∞ X l þ Δx dyI0 p p −α ðλex Þð2plþxÞ α n ðλex Þ∫l−Δx R Rmir e þ Δy p¼0 ! ∞ X p−1 p −α ðλex Þð2pl−xÞ R Rmir e dx: 2 2 p¼1 ð2Þ
© Copyright 2024