Simultaneous Quantification of Multiple Polycyclic Aromatic Hydrocarbons in Aqueous Media using Micelle Assisted White Light Excitation Fluorescence

Qualitative and quantitative display of multiple fluorescent analytes is made simple and reliable in this micelle assisted methodology. The adopted method involves micelle assisted evincing of ppb level of PAHs in water; measurement of total fluorescence (white light excitation fluorescence, WLEF) and data deciphering using multivariate analysis. This protocol yields sensitive and accurate quantification of the cancerous pollutants (PAHs) in aqueous media with Limit of Quantification of the order 1–10 μg/L and accuracy of >98%. The use of WLEF enables the simultaneous acquisition of fluorescence signatures of all the PAHs. It has the additional advantage of being portable, layman-friendly and cost-effective. The optimized amount of surfactants for the simultaneous extraction of PAHs from real samples was estimated as 27.8 mg (19.3 mM) of SDS and 9.1 mg (5 mM) of CTAB. Also, the analytical fidelity of the quantification such as percentage recovery (98 ± 2%), linear dynamic range (2–250 μg/L), RMSEP (<0.5), etc. explains the veracity of methodology.

Synchronous Fluorescence Scan (SFS) and Total Synchronous Fluorescence Scan (TSFS) 19,20 . Though these techniques are fairly robust and accurate in predicting analyte composition, they are cumbersome, time-consuming and require scientific expertise for data collection as well as for analysis. Real-time samples often demand in situ analysis. For that optical spectrometer should reach the sample, which is almost impossible with conventional spectrometer design. i.e., the desktop design lacks the features of portability and flexibility in optimizing sample-specific geometry. White Light Excitation Fluorescence (WLEF) is evolved recently as a technique to measure total fluorescence in a two-dimensional (2-D) plot. It could investigate multiple fluorophores simultaneously without reconfiguring instrument parameters 21,22 . Also, the WLEF technique requires only a low power light source and hence minimizes the chance of photo-bleaching. These added advantages will help us to integrate a portable device for online/real-time quantification of multiple PAHs in water.

Objective
In this work, simultaneous real-time analysis of multiple PAHs (multi-probing) in aqueous media is proposed using a Dip Probe Fiber Optic Spectrometer and a highly sensitive WLEF technique with the assistance of the micellar system. The analytical complexity of multifluorophores due to self-quenching, energy transfer, etc. could be nullified under the condition that there is only one fluorophore per micelle on an average. This can be achieved by creating a micelle system using a suitable concentration of surfactant, which in turn restores spectral additive nature. This spectral additive nature of PAHs in the micelle driven system can be used for their simultaneous determination and precise quantification in aqueous media by employing multi-component regression analysis.

Materials and Methods
Chemicals. PAHs  to 250 μgL −1 in the test solution. The solutions were prepared by transferring the required amount of PAH solution in a volatile solvent (acetone) into a sample container followed by solvent evaporation. It is then dissolved in the required amount of surfactant solution. The concentration of CTAB and SDS in the sample test solution was optimized and maintained as 5 mM and 20 mM respectively for PAHs quantification.
Synthetic mixtures for multi-probing. A series of synthetic mixtures (5 components, 3 sets {A, B and C}), were made from the standard stock solutions (1000 ng mL −1 ) for each of the PAHs. The sample matrix is made such that there is enough spectral overlap between their emission profiles. The synthetic mixtures are solubilized randomly using anionic and cationic surfactants as follows, Set A and B are extracted with CTAB (5 mM) and the Set C is with SDS (20 mM). The WLEF intensities of the mixtures are measured using custom-fabricated fiber optic compatible portable spectrometer subsequent to micelle induced disaggregation. WLEF instrumentation and data acquisition. In our previous publication, WLEF has been introduced as the 2-D analogue of the 3-D total fluorescence spectrum, EEMF, which collects total fluorescence response from the sample and display a plot of emission intensity against wavelength 21 . WLEF intensity at a specific wavelength is proportional to the excitation efficiency, integral product of absorption profile of molecule and the lamp profile of an exciting source.
WLEF has two main advantages, viz., • It features simultaneous excitation (WLE) of all the fluorophores • WLE improves emission intensity by populating more molecules in the excited state.
The reduction of data dimension (3D to 2D), fast data acquisition, spectral veracity, sim Dip Probe Fiber Optic spectrometerplicity in analysis, cost-effectiveness and portability are additional features of the WLEF method. The perks of WLEF are already exploited in monitoring quality of fuel blends 23 , biofluids 24 , etc.
A fiber optic compatible and portable Dip Probe Fiber Optic Spectrometer 25 was fabricated for the WLEF measurement using a white light excitation source (270-900 nm, W-X lamp); sensitive CCD-Diode Array Detector (DAD, QE9000-Pro) and UV-VIS XRS Solarization Resistant fibers (1 m, transmission range:180-900 nm, diameter: 600 μm and NA: 0.22). The components were purchased from Ocean Optics Inc. The WLEF spectra were collected using a personal computer connected to a fiber-optic spectrometer. The Ocean View software procured from Ocean Optics Inc. eases data acquisition.
The schematic optical scheme of the Dip Probe Fiber Optic Spectrometer is shown in Fig. 1 (A photograph of the Dip Probe design is given as ESI-1). The dip probe is fabricated such that the fluorescence emission was collected at the right-angled geometry to the excitation beam, which minimizes spectral contamination by reflected and transmitted light. The volume of the test solution required for the measurement is 5 mL and the path length maintained as 10 mm. Optimized instrument parameters for WLEF measurement are 10 seconds of spectral integration, averaged over 5 spectra and boxcar smoothening across 5 consecutive wavelengths. These optimized parameters gave reasonable highest WLEF intensity for each species, which was chosen for the entire analysis.
Multivariate analysis of WLEF data. The WLEF intensity at each wavelength depends on two variables, viz., concentrations of species and their respective WLEF intensities, i.e., the data matrix has a bilinear structure. The uniqueness of Multivariate Curve Resolution-Alternating Least Square (MCR-ALS) analysis is that it provides analyte-specific pure spectral (instrument) responses along with their abundance as concentration profiles in the output. Though the exploratory methods like independent component analysis (ICA), principal component analysis (PCA), etc. are good for resolving bilinear data, the MCR-ALS algorithm was preferred over them to analyze the WELF data matrix. Also, ICA or PCA methods lack a genuine chemical or scientific model. MCR-ALS algorithm alternatively optimizes both the spectral and concentration profiles in each iterative cycle. Iterative methods do not require a structured concentration direction. MCR always furnishes a simple bilinear description based on the variation in the data especially the natural properties of the problem of interest. MCR-ALS algorithm has given the freedom of the incorporating chemical and mathematical constraints to make the data analysis task efficient. These characteristics make MCR outputs legible to non-chemometricians [26][27][28] .
The PLS-Toolbox 5.0.3 written in MATLAB language was used for MCR-ALS analyses as follows. The WLEF intensity of the mixture of PAHs at a particular emission wavelength, W i is given by where w i and c i are respectively the WLEF intensity and the fractional concentration of the i th component. The value of 'i' is ranging between 1 and n, the total number of the component present in the mixture. Hence the quantification needs to estimate concentration c i 's along with the unique molecular feature w i . A set of 'n' linear equations with 'n' fractional coefficients can be resolved as follows.
The relationship between the WLEF data matrix W i is bilinear in nature, i.e., it is the sum of the product of the concentration of each PAHs and the respective spectra as expected.
where W is the WLEF spectra of I samples against J wavelengths; C is the concentration of K number of PAHs in the analytes; S T are the pure WLEF spectra, its rows contain the pure spectra of respective species; E is the calibration error 26 .
Prior to MCR analysis, the number of components and the spectral window are iteratively optimized by means of PCA followed Evolving Factor Analysis (EFA) procedures 29 . The initial estimation for the spectral profile (determination of the purest variables) is done using EFA followed by MCR in order to initiate the iterative

Results and Discussion
Selection of PAHs and its sensitivity in WLEF measurement. The PAHs were chosen based on their insoluble nature and cancerous properties (Table 1). Since they last in the environment for quite a long time (half-life ~ 6 yrs), these living carcinogens are monitored in aqueous media. The acceptable limit of PAHs in the environment, especially water bodies, is less than 1 to 5 mg L −1 .
The WLEF spectra of selected PAHs measured in 5 mM CTAB are shown in Fig. 2.
The merits of this analytical procedure for each PAH are independently established by making a linear plot of WLEF intensity against the amount of PAH. WLEF intensity in the concentration ranges from 10 to 100 nM. The WLEF spectra and their respective calibration plots of selected PAHs are presented as ESI-2. The Linear Quantification Range (LQR), Pearson's regression coefficient (R 2 ) and Limit of Detection (LOD) were determined for each PAH and gathered in ESI-3.
The LDR for the analysis of PAHs by WLEF ranged from 2.5 to 250 μgL −1 with very good correlation (R 2 ~0.9989 to 0.9999). The LOD of various PAHs was found to be between 0.1 and 2 μgL −1 . Figure 2 clearly reveals that there is a robust spectral overlap among the emission spectrum of PAHs. The linearity plot of WLEF intensity against concentration in the analytical range (ESI-2) verifies that the possibility of energy transfers, inner filter effects, etc., is negligible in micellar media. In this regard, though the simultaneous quantification of multiple components has a good impact, a suitable data collection methodology followed by logical analysis is also required.
Micelle assisted probing of multiple PAHs in aqueous media. The micelle assisted multi-probing aims at instantaneous quantitative visualization of multiple PAHs in real samples. It involves 3 steps, viz., (a) micellar dissolution of PAHs to enhance the fluorescence as well as remove spectral interferences, (b) measurement of WLEF to obtain total fluorescence features in a very short time and (c) spectral deconvolution using chemometrics to provide qualitative and quantitative information.
Micelle assisted dissolution of PAHs. Micelles have two roles in this method, (a) minimize the problems of non-linearity and regain spectral additivity and (b) provide a hydrophobic environment to the analyte and make it visible. Micelle scavenges each PAH from aqueous media to its micellar cavity and increases the average distance between them which leads to enhanced sensitivity, reduced non-radiative energy transfer, and negligible collisional quenching 3 .
Micellar solubilization of PAHs (BaP as a representative example, which is insoluble in water) was studied by fixing the concentration of it and varying the amount of surfactant, CTAB (cmc = 0.96 mM; CAN = 62), in the test solution (Fig. 3). Enhanced solubility or fluorescence is observed for the sample above cmc of surfactant, due to the micellar scavenging of fluorescent and hydrophobic PAH. Saturation in fluorescence intensity was observed above 5 mM of CTAB (micelle concentration, [M] =65 μM). So the optimum concentration of the micellar medium, CTAB, was selected as 5 mM to scavenge PAH from aqueous media with minimal intensity quenching effects. Similar study has been done with an anionic surfactant, SDS (cmc = 8.3 mM; CAN = 63) and optimized concentration was found to be 19  The variables 'x' and 'n' are the PAH to micelle ratio and the number of PAH respectively. For the given micellar (CTAB, 5 mM) and PAH (0.75 μM) concentrations, the chance of one PAH being inside a micelle is about 1%. Likewise, the probability of two PAHs in a single micelle is 4.36 × 10 −5 . The table, ESI-5, clearly indicates that the probability of finding two or more PAHs in the same micelle is rarest or can be omitted.
The fluorescence negating phenomena of resonance energy transfer and the dynamic collisional quenching etc. are thus nullified by pushing each PAHs into separate micelles and hence by increasing the average spacing between any two PAHs. In addition, micellization enhances solubilization of hydrophobic PAHs and hence fluorescence intensity. The effect of the second fluorophore on the calibration curve of a target fluorophore was established by gradually varying the concentration of second (interfering) fluorophore. The interference between different fluorophore and spectral additivity of PAH mixtures were also investigated and were found to be negligible in micellar media. The optimized concentration of the surfactant regains the spectral additivity of target as well as interfering fluorophores with an analytically relevant range (2 to 250 μgL −1 ). The amount of surfactants required for the complete evincing of multifluorophoric mixture (sample volume -5 mL) was 9.1 mg (5 mM) and 27.8 mg (19.3 mM) respectively for CTAB and SDS. Since the pH and presence of salts (ionic strength) can significantly alter the cmc, the optimized surfactant concentration may vary as a function of the concentration of other ionic species in the real system. This must be factored in while taking real-life samples. Non-ionic surfactants with a zeta potential of the micelles close to zero may be the most ideal choice. MCR-ALS analyses can decompose the WLEF spectral data matrix into concentration and spectral data. To get meaningful information a series of logical and mathematical constraints are employed. Initially, non-negativity constraint is applied to both spectral and concentration profiles because of their real nature. Since the spectra of PAHs are having vibrational structures, unimodality constraints cannot be a logical choice, so unimodality constraint is applied only to the concentration profile 29,32 . The correlation constraint is smartly introduced to predict the exact concentration of analytes in an unknown environment.
A soft MCR-ALS procedure performed on calibration data with non-negativity constraints on both concentrations and spectra profiles, and unimodality on concentration profile. PCA initially estimates the number of components as 5 which captures the maximum variation of about 99.9% for all the three sets. This is matching with the composition of the spiked samples as well. EFA yields the spectral window (variable) in which the analytes www.nature.com/scientificreports www.nature.com/scientificreports/ are present; only two components have less overlap with the rest of the spectra, the other three are completely overlapped. The MCR-ALS resolved 5 species for all the three sample sets (A, B and C). Spectral deciphering and quantification of multiple analytes in complex mixture, Set A (mixtures Set B and C are shown in ESI-6a & b) using MCR-ALS multivariate analysis is shown in Fig. 4. The exactitude of resolution is evaluated by considering 3 parameters, viz., Standard Deviance of residuals, the fitting error and the variance and is tabulated in ESI-7.
The correlation constraint creates an MCR-ALS calibration model, which can be adopted for the quantitative visualization of the PAHs in real samples; even unknown spectral interferences have existed. Prior knowledge of the spectral nature helps us to identify the PAH along with its amount in the mixture accurately. Percentage (%) recovery of each PAH present in the 3 sets of calibration samples was determined (ESI-8) using the following relations.
where x is the spiked amount and x' is the recovered amount of PAH. The root mean square error in prediction (RMSEP) of the analyte combinations for each sample determined using the following relation.
Relative percentage recoveries of sample mixtures were found to be 98 ± 2% and lack-of-fit in all the analytes were well below 1%. So, this micelle assisted methodology for probing multiple PAHs (multi-probing) in aqueous media using WLEF could be a low cost, portable option. Also, this methodology can be easily adapted for the real-time monitoring of PAHs in different water bodies.
Real sample analysis (% recovery of PAHs). The real analytical system carries complexity with respect to matrix fluorescence (mainly, humic substance) and as well as multiple fluorophoric components. The reliability of the methodology was tested by a river water sample obtained from Vettar-river spiked with 3 synthetic mixtures of PAHs. The amount of SDS used to extract the total fluorescence was 30 mg (optimized value 27.8 mg). SDS is adopted for real sample analysis because of its environment-friendly nature. The WLEF spectra of the spiked samples are resolved using the developed model (Fig. 5). The river water sample was collected and filtered through disposable syringe filters (MILLEX, 33 mm, and 0.22 μm) to remove particulate matters, followed by PAHs spiking in order to minimize the scattering and particulate matters. The humic substance present in the fairly clear river water sample was found to be 142 ± 5 ppm 33 . Correlation constraint eases the adoption of developed MCR-ALS model to predict the quality of water even in the presence of background fluorescence interference from humic substances and other fluorescing contaminants. The calibration model designed for a test sample has been used to predict the quantity of the 'unknown' mixture in river water samples.
The percentage recovery of each PAHs present in the 3 sets of samples are tabulated with its figures of merit ( Table 2). The RMSEP values show that the spectral subtraction will predict the sample composition with more precision. The spectral subtraction using a suitable baseline (TDW is used here) removes unwanted and inherent scattering components such as Rayleigh and Raman scatterings. The percentage of recovery studies show that baseline-corrected measurement has a feeble deviation from the original values (~2%) and the RMSEP is well below 0.3 units. These parameters explain the reliability and fidelity of the methodology in real samples. The uniqueness of spectral features and the availability of fluorescence spectral libraries of PAHs widen the scope of our methodology. The simplicity of instrumentation, data collection in no time, ease of data processing, portability of equipment, etc. are the additional features of this micelle assisted multi-probing. It allows real-time monitoring and simultaneous quantification of a mixture of PAHs in water bodies. This proof-of-concept work well with  www.nature.com/scientificreports www.nature.com/scientificreports/ low scattering samples (even PAH spiked river samples) and merely pure with respect to other ionic contaminants and low salinity samples.

Conclusions
The micelle assisted fluorescence enhancement followed by multivariate deciphering methodology may be devised for the real-time analysis of complex, but transparent analytes. The micellar scavenging made the 'invisible' carcinogens (PAHs) visible by providing a more hydrophobic environment and increasing solubility of PAHs. The optimized amount of surfactants required for the micellar extraction of PAHs is determined as 19.3 mM for SDS and 5 mM for CTAB respectively. WLEF of micellar extracted multifluorophoric systems regains spectral addition feature by optimizing the maximum number of PAHs in a micelle to unity and hence minimize self-quenching, energy transfer, etc. The WLEF spectra of mixtures are deciphered into quality and quantity information using multivariate MCR-ALS. Also, the LOD of various PAHs in water was found to be of the order of micrograms per liter (subnanomolar). The reliability and fidelity of the methodology are tested with real samples which yield more than 97% of accuracy with ~0.3 of RMSEP. The issues of Rayleigh and Raman scattering in the data set are easily removed in this methodology by baseline subtraction. In summary, Micelle assisted extraction followed by MCR-ALS analysis of bilinear structured WLEF data can easily probe multiple fluorescent components in a facile manner. Also its application potential under various matrix conditions like scattering/intrinsic fluorescence/pH/salinity/hardness etc. are yet to be explored.