Confocal Raman Spectroscopic Imaging for Evaluation of Distribution of Nano-Formulated Hydrophobic Active Cosmetic Ingredients in Hydrophilic Films

Film-forming systems are highly relevant to the topical administration of active ingredients (AI) to the body. Enhanced contact with the skin can increase the efficacy of delivery and penetration during prolonged exposure. However, after the evaporation of volatile solvents to form a thin film, the distribution of the ingredient should remain homogenous in order to ensure the effectiveness of the formula. This is especially critical for the use of hydrophobic molecules that have poor solubility in hydrophilic films. In order to address this concern, hydroxyphenethyl esters (PHE) of Punica granatum seed oil were prepared as a nanosuspension stabilised by poloxamers (NanoPHE). NanoPHE was then added to a formulation containing polyvinyl alcohol (PVA) as a film forming agent, Glycerol as a plasticiser and an antimicrobial agent, SepicideTM HB. Despite their reliability, reference methods such as high-performance liquid chromatography are increasingly challenged due to the need for consumables and solvents, which is contrary to current concerns about green industry in the cosmetics field. Moreover, such methods fail to provide spatially resolved chemical information. In order to investigate the distribution of ingredients in the dried film, Confocal Raman imaging (CRI) coupled to Non-negatively Constrained Least Squares (NCLS) analysis was used. The reconstructed heat maps from a range of films containing systematically varying PHE concentrations highlighted the changes in spectral contribution from each of the ingredients. First, using NCLS scores it was demonstrated that the distributions of PVA, Glycerol, SepicideTM HB and PHE were homogenous, with respective relative standard deviations (RSD) of 3.33%, 2.48%, 2.72% and 6.27%. Second, the respective relationships between ingredient concentrations in the films and their Raman responses, and the spectral abundance were established. Finally, a model for absolute quantification for PHE was be constructed using the percentage of spectral abundance. The prepared %w/w concentrations regressed against predicted %w/w concentrations, displaying high correlation (R2 = 0.995), while the Root Mean Squared Error (0.0869% w/w PHE) confirmed the precision of the analysis. The mean percent relative error of 3.75% indicates the accuracy to which the concentration in dried films could be determined, further supporting the suitability of CRI for analysis of composite solid film matrix. Ultimately, it was demonstrated that nanoformulation of hydrophobic PHE provides homogenous distribution in PVA based film-forming systems independent of the concentration of NanoPHE used in the formula.


Introduction
The effectiveness of a topical product depends on the physicochemical properties of the active substances and on the capacity of the formula to maintain close contact with the skin for prolonged exposure, enabling the diffusion of the molecules in or through the skin. Topical film-forming systems are increasingly studied for drug delivery because they adhere to the skin, forming a thin transparent film in situ and providing delivery of the active ingredients (AI) to the body tissue [1][2][3]. A cosmetic film-forming system contains three types of components: the active ingredient(s), the excipients (film-forming polymer, plasticizers or additives), and volatile solvent(s) [4] that evaporate upon contact with the skin to form the film.
The resulting film should form a homogenous reservoir of the active ingredient on the skin. It has been reported in the literature that only drugs in molecular (as opposed to solid) form in the residual film can diffuse [4]. During the film-forming process, the evaporation of volatile components can change the characteristics of the formulation and even produce drug precipitation on the skin surface (so called vehicle metamorphosis [12]). This phenomenon is particularly observed for hydrophobic AI dispersed in a hydrophilic film. To prevent the precipitation of AI, effective precipitation inhibitors [13] can be added to the mixture, or the hydrophobic AI can be encapsulated [14].
Physicochemical characterisation of newly developed products can be a pivotal step in optimising and refine formulae to reach the desired texture and sensorial properties without loss of efficacy. Commonly, topical films are characterised in terms of their morphology and mechanical properties. For example, surface texture (smooth or rough) of the film can be observed using light microscopy, scanning electron microscopy (SEM) [15][16][17][18] and atomic force microscopy (AFM) [19][20][21]. Analytical control to confirm the concentration of ingredients in films can be performed with gold standard techniques such as highperformance liquid chromatography (HPLC) [17]. However, to determine the content of an AI, fastidious extraction protocols need to be applied prior to analysis. Firstly, such approaches fail to deliver spatially resolved information about AI distribution in the films. Secondly, there is an increasingly high demand from the cosmetic industry for greener methods that are environmentally responsible, and hence alternatives with lower requirements of consumables and solvents are sought.
Vibrational spectroscopy, i.e., infrared absorption and Raman scattering spectroscopy, is a non-destructive, label-free and reagent (solvent)-free technique that provides molecular fingerprints, enabling the chemical composition of samples to be probed. Fourier transform infrared spectroscopy (FTIR) has also been used for the characterisation of films [1,15] and or to identify the molecular interactions between active ingredients and the other constituent components of films [16,19]. Use of Raman spectroscopy has been reported for study the composition of various thin polymer films [14,20,21]. Most studies have relied on point measurements [1,15,[19][20][21], whereas coupling the molecular specificity of the analysis with the micrometric spatial resolution of Raman confocal microscope offers a powerful non-destructive tool to study the local distribution of AI in cosmetic films. Confocal Raman Imaging (CRI) has been extensively studied in, for example, the sub-cellular analysis of single cells to investigate physiological processes [22] or the internalization of anticancer drugs [23] and nanoparticles [24] as well as the cartography of tissue sections from different organs for diagnostic purposes [25], monitoring AI penetration into the skin [26] and the analysis of solid forms such as pharmaceutical tablets for quality control [27] and identifying counterfeits [28].
More recently, a study using Raman cartography highlighted that hydrophobic compounds from Punica granatum seed oil hydroxyphenethyl esters (PHE) stabilised in a poloxamer-based nanosystem showed very homogeneous distribution in hydrophilic PVA films compared to the free form introduced as an ethanol solution [14]. While quantitative analysis of AI concentrations in cosmetic and pharmaceutical products using Raman spectroscopy is increasingly documented for solutions and semi-solids forms (gels, cream) [29][30][31], application of the technique to solid complex mixtures remains largely limited to semi-quantitative analysis [32][33][34].
In the current study, Raman spectroscopy coupled with Non-negatively Constrained Least Squares (NCLS) analysis is used to deliver reconstructed chemical maps to evaluate not only the distribution but also the constituent chemical content of nanoformulated Punica granatum seed oil hydroxyphenethyl esters (NanoPHE) in thin hydrophilic PVA films containing Glycerol as a plasticiser along with an antimicrobial agent. Increasing concentrations of NanoPHE were incorporated in PVA-based films to study the homogeneity of the AI distribution and evaluate the compatibility of the formula with a broad range of concentrations. Moreover, achieving quantification of the ingredients using Confocal Raman Imaging (CRI) is addressed, while the technique is proposed as a non-destructive and reagent-free approach that can support the shift towards green chemistry and environmentally responsible industrial processes in the cosmetic field.

Sample Preparation
Nanoformulation of PHE (NanoPHE: NanoPHE were prepared by a green emulsionsolvent evaporation process assisted by ultrasound, as described in a previous study [14]. Briefly, the aqueous phase was prepared by dissolving 600 mg of Pluronic ® F127 in 10.5 g of ultrapure water. The organic phase was prepared by dissolving 1225 mg of DL and 80 mg of Vitamin E in 2.70 g of dimethyl carbonate. The nanoemulsion was obtained by sonicating the mixture during 2 min in a glass rosette cell (20 kHz, amplitude 30%, Vibracell 75041-Sonics). The organic solvent was evaporated (60 mbar 1 h, 25 • C, LABOROTA 40001-Heidolph) to obtain nanodispersion (NanoPHE). The volume of the preparation was adjusted to 20 mL with ultrapure water to ensure the reproducibility of the concentration.
Characterisation of NanoPHE: The hydrodynamic diameter (D H ) and polydispersity index (PdI) of NanoPHE were measured at 25 • C using a dynamic light scattering (DLS) instrument (NanoZS, Malvern Instruments, Orsay, France) piloted by Zetasizer Software (zetasizer software.ink, version 7.10, Malvern panalytical, Worcester, UK). Each sample was diluted 1:10 in ultrapure water, then measured using a 633 nm laser with a detection angle of 173 • . NanoPHE zeta potential was measured at the same dilution and with the same equipment, with a detection angle of 13 • . All measurements were carried out in triplicate.
PVA solution was prepared under heating (90 °C) and vigorous stirring until complete dissolution in ultrapure water (approximately 4 h), then cooled to room temperature. In order to avoid potential degradation of other components, PVA solutions were prepared upstream, then Glycerol, AMA, and NanoPHE were added. Glycerol and AMA were first added to the PVA solution and stirred for 20 min; finally, the NanoPHE suspension was added and stirred for 30 min before film preparation. For each sample, 10 g of the final film-forming solution was prepared ( Table 1). The initial PVA solution was adjusted for C0-C5 to keep the final concentration at ~10% w/w. Glycerol and AMA had respective final concentrations of ~10% w/w and 1% w/w. Nano-PHE concentration ranged from 1.15% w/w (C1) to 10.218% w/w (C5). The final nominal concentration (% w/w) for each component in the dried films is provided in Table 2.  PVA solution was prepared under heating (90 • C) and vigorous stirring until complete dissolution in ultrapure water (approximately 4 h), then cooled to room temperature. In order to avoid potential degradation of other components, PVA solutions were prepared upstream, then Glycerol, AMA, and NanoPHE were added. Glycerol and AMA were first added to the PVA solution and stirred for 20 min; finally, the NanoPHE suspension was added and stirred for 30 min before film preparation.
For each sample, 10 g of the final film-forming solution was prepared ( Table 1). The initial PVA solution was adjusted for C0-C5 to keep the final concentration at~10% w/w. Glycerol and AMA had respective final concentrations of~10% w/w and 1% w/w. NanoPHE concentration ranged from 1.15% w/w (C1) to 10.218% w/w (C5). The final nominal concentration (%w/w) for each component in the dried films is provided in Table 2. Films were prepared on Superfrost ® histological glass slides (ThermoScientific, Waltham, MA, USA; Menzel-Gläser, Braunschweig, Germany) by scraping the solution between four adhesive tapes used to form a rectangular shaped mould. Another clean glass slide was used to spread the sample over the entire length of the substrate, leaving a liquid film with a constant thickness. Samples were left to air-dry overnight in darkness. Films for each concentration were prepared in triplicate and identified as R1, R2 and R3.

Confocal Raman Imaging
Acquisition: Raman images were acquired using an Alpha 300 R microspectrometer (WITec, Ulm, Germany). Spectral data were recorded using a 532 nm excitation laser source with a power of 20 mW at the sample. The instrument is calibrated daily using a two-step procedure. First, the TrueCal function of Project 5 (WITec, Ulm, Germany) was used; this is an automatic multipoint calibration routine performed with a Mercury-Argon (HgAr) light source integrated in the Raman microscope. Second, prior to data acquisition, verification was done using the peak at 520.7 cm −1 from a silicon substrate. The laser was focused onto the sample through a 50× objective (Zeiss, LD Epiplan HD, N.A. 0.5). Backscattered light was detected by dispersing Raman-shifted radiation onto a Deep Depletion CCD detector (1024 × 128 pixels) using a 600 lines/mm grating centred at 2500 cm −1 , resulting in a spectral resolution of ≈4 cm −1 . The acquisition time was set to 10 s, with two accumulations. 100 × 100 µm 2 Raman maps were recorded with an X and Y lateral step size of 4 µm (corresponding to 625 spectra). For each concentration and condition, three separate films were prepared (R1, R2, R3), and the measurements were performed once for each sample. Additionally, reference Raman spectra were recorded for all constituent compounds of the films.
Data handling: Raman maps were pre-processed and analysed using Matlab ® (Mathworks, Portola Valley, CA, USA). Spectra were cut in the range 300 to 1800 cm −1 and then subjected to offset correction followed by vector normalisation.
Principal Components Analysis (PCA): PCA is a multivariate analysis technique that is widely used to simplify a complex data set of multiple dimensions. It permits a reduction in the number of variables in a multidimensional data set while retaining most of the variation within the data set. The other advantage of this method is the derivation of component (PC) loadings, which represent the variance of each variable (wavenumber) for a given PC. Analysis of the loading vector of a PC can provide information about the source of the variability inside a data set derived from variations in the chemical components contributing to the spectra [26].

Non-negatively Constrained Least Squares analysis (NCLS):
NCLS is a method for the unmixing of spectral data described by Kwan et al. [35]. The NCLS algorithm calculates the signal contribution of the molecule of interest in each experimental Raman spectrum collected from the sample based on a set of reference spectra. A simplified description of the NCLS outcome can be defined as Equation (1): where S S -simulated spectrum for a given pixel of the Raman map, S Ri -reference spectra, F i -abundance fraction estimated, i-number of reference spectra included in the model and R-the residual. The algorithm aims to calculate an S S which is as close as possible to the experiential spectrum, with minimal residual. The NCLS algorithm requires reference spectra in order to calculate the relative spectral contributions of each compound for every pixel of the Raman map. Calculation cannot be performed with a single reference spectrum [26], and it has been demonstrated by Miloudi et al. that increasing the number of replicates for a given reference compound significantly increases the reliability of the outcome [26]. Therefore, for each component (free or encapsulated PHE, Vitamin E, Pluronic ® F127, Glycerol, and AMA), ten spectra were used in the NCLS. Concerning the PVA, a control film containing only the polymer was prepared and a total of 1875 reference spectra were collected and used in the NCLS. For clarity, the results are presented as the mean contribution of each compound, calculated from the n replicates of reference spectra.
First, the results were presented as reconstructed heat maps, with a scale bar ranging from blue (low concentration) to red (high concentration) reflecting the spectral contribution calculated from reference spectra, as determined by NCLS fitting. Individual maps were therefore obtained for each individual reference compound used for the NCLS calculations. The obtained NCLS score corresponds to the spectral abundance fraction of a given ingredient, and it is between 0 and 1. A value of 1 is equal to 100% of the reference spectrum used in the calculation.
Second, NCLS results were expressed as the percentage of spectral abundance. The contribution F of a given compound was normalised to 1 by dividing its score obtained by NCLS by the sum of the scores of all references used, then multiplied by 100, as follows: where X-PVA, Glycerol, AMA or PHE, Cn-the concentration of NanoPHE in the film, n-the number of concentrations tested (from 0 to 5), F-the spectral abundance fraction estimated from NCLS, i-the number of reference compounds included in the model. The spectral abundance was then used to construct linear regression plots. It should be noted, however, that while the terms spectral abundance and percent spectral abundance are used in the description above, each constituent compound has a different Raman cross-section, and therefore can contribute more or less strongly to the Raman spectrum of the film. The percentage-based spectral maps thus provide a semiquantitative representation of the relative contributions, but do not directly correspond to the actual percentages of the concentrations.
Thirdly, in an effort to take the final step towards full quantification, a calibration model for PHE is constructed using the percentage of spectral abundance as a function of %w/w PHE in dried films. The results for maps R1 of each concentration (C0 to C5) were used as data points for the calibration curve, then the percentages of spectral abundances calculated from the maps R2 and R3 were projected in the model as unknown samples to be determined. A linear regression of the predicted concentration against the prepared concentrations resulted, with the linearity (coefficient of determination-R 2 ) and Root Mean Square Error (RMSE) used as criteria of the model's reliability, with RMSE defined as follows: where n-number of observations, y i -observed values, andŷ i -predicted values.

Preparation and Characterisation of NanoPHE
Pluronic ® F127 is a non-ionic amphiphilic triblock polymer, poly(ethylene oxide) 101poly(propylene oxide) 56 -poly(ethylene oxide) 101 . Composed of two hydrophilic blocks (polyethylene oxide, PEO) and one hydrophobic (polypropylene oxide, PPO), the micellar organisation of Pluronic F127 in aqueous solution induces a core-shell nanostructure. In this study, the hydrophobic active ingredient PHE has been housed in the polypropylene oxide hydrophobic core, producing NanoPHE structures. The nanoformulation size was characterized by DLS. NanoPHE aqueous suspensions typically displayed hydrodynamic diameter values of 126.3 ± 1.0 nm and exhibited low PdI values (0.107 ± 0.002), indicative of the monodispersity of the suspension. The zeta potential was slightly negative. Electrostatic repulsion of the particles was not an issue, as the colloidal stability was expected to be ensured by the steric hindrance of the PEO chains of Pluronic ® F127, stabilising the system [36].
The properties of NanoPHE were similar to those previously reported, highlighting the reproducibility of the protocols [14]. These nanostructures, with a D H close to 120 nm, can be considered as suitable for skin administration [37].

Spectral Characterisation of NanoPHE
The Raman spectrum of PHE ( Figure 2A) exhibits two strong features, at 1636 cm −1 , corresponding to C=C stretching, and at 1169 cm −1 , corresponding to C-C and C-O stretching. Other bands in the range of 800 cm −1 -1450 cm −1 display lower intensities, but can still be easily observed in the enlarged inset of Figure 2. The corresponding vibrational mode assignments are listed in Table 3.
poly(propylene oxide)56-poly(ethylene oxide)101. Composed of two hydrophilic blocks (polyethylene oxide, PEO) and one hydrophobic (polypropylene oxide, PPO), the micellar organisation of Pluronic F127 in aqueous solution induces a core-shell nanostructure. In this study, the hydrophobic active ingredient PHE has been housed in the polypropylene oxide hydrophobic core, producing NanoPHE structures.
The nanoformulation size was characterized by DLS. NanoPHE aqueous suspensions typically displayed hydrodynamic diameter values of 126.3 ± 1.0 nm and exhibited low PdI values (0.107 ± 0.002), indicative of the monodispersity of the suspension. The zeta potential was slightly negative. Electrostatic repulsion of the particles was not an issue, as the colloidal stability was expected to be ensured by the steric hindrance of the PEO chains of Pluronic ® F127, stabilising the system [36].
The properties of NanoPHE were similar to those previously reported, highlighting the reproducibility of the protocols [14]. These nanostructures, with a DH close to 120 nm, can be considered as suitable for skin administration [37].

Spectral Characterisation of NanoPHE
The Raman spectrum of PHE ( Figure 2A) exhibits two strong features, at 1636 cm −1 , corresponding to C=C stretching, and at 1169 cm −1 , corresponding to C-C and C-O stretching. Other bands in the range of 800 cm −1 -1450 cm −1 display lower intensities, but can still be easily observed in the enlarged inset of Figure 2. The corresponding vibrational mode assignments are listed in Table 3. The Raman spectrum of NanoPHE ( Figure 2B) is similar to that of pure PHE, with a significant band shift in the C=C stretching of PHE from 1636 cm −1 to 1639 cm −1 , explained by molecular interaction with other components from the shell or the core of the nanosystem. The other features do not appear to be modified by the nanoformulation.
The PCA scatter plot in Figure 3A shows a clear separation along Principal Component 1 (PC-1) (94% of the total explained variance) between the data corresponding to PHE  The Raman spectrum of NanoPHE ( Figure 2B) is similar to that of pure PHE, with a significant band shift in the C=C stretching of PHE from 1636 cm −1 to 1639 cm −1 , explained by molecular interaction with other components from the shell or the core of the nanosystem. The other features do not appear to be modified by the nanoformulation.
The PCA scatter plot in Figure 3A shows a clear separation along Principal Component 1 (PC-1) (94% of the total explained variance) between the data corresponding to PHE spectra (black dots) and NanoPHE spectra (red dots). The loading plot corresponding to PC-1 (Figure 3(Ba)) shows that the main features driving the discrimination between PHE and NanoPHE are the 1636-1639 cm −1 peaks, corresponding to C=C stretching of the aromatic ring of PHE.

Reference Compounds
Reference Raman spectra of NanoPHE, PVA, Glycerol and AMA were recorded to identify the characteristic vibrational bands of the molecules (Figure 4).

Reference Compounds
Reference Raman spectra of NanoPHE, PVA, Glycerol and AMA were recorded to identify the characteristic vibrational bands of the molecules (Figure 4).  [41]. The Raman spectrum of Glycerol ( Figure 4B), in the 1000-1500 cm −1 spectral region, contains two broad bands that encompass modes due to CCO stretching, CH2 twisting, and COH stretching vibrations of the primary -H groups. The bands at 1060 and 1111 cm −1 are assigned to the symmetric stretching mode of the primary and secondary -OH groups, respectively. The bands in the 1207-1310 cm −1 spectral region and the band at 1468 cm −1 are assigned to CH2 twisting [42].

PVA-GS-NanoPHE Films
Each pixel from a Raman map contains a full spectrum resulting from the combined contribution of all constituents. Figure 5 displays the mean spectra obtained from PVA-GS-NanoPHE films with increasing concentrations of NanoPHE (C0 to C5).
The mean spectrum of C0 exhibits the spectral features of the films formed from the mixture without NanoPHE. The signature is mainly due to the combined contribution of PVA and Glycerol, most of whose peaks overlap, for example in the spectral ranges 825-855 cm −1 , 920-928 cm −1 , 1060-1128 cm −1 , and 1440-1468 cm −1 . The contribution from AMA is limited to the 1001 cm −1 and 1612 cm −1 features.  [41]. The Raman spectrum of Glycerol ( Figure 4B), in the 1000-1500 cm −1 spectral region, contains two broad bands that encompass modes due to CCO stretching, CH 2 twisting, and COH stretching vibrations of the primary -H groups. The bands at 1060 and 1111 cm −1 are assigned to the symmetric stretching mode of the primary and secondary -OH groups, respectively. The bands in the 1207-1310 cm −1 spectral region and the band at 1468 cm −1 are assigned to CH 2 twisting [42].

PVA-GS-NanoPHE Films
Each pixel from a Raman map contains a full spectrum resulting from the combined contribution of all constituents. Figure 5 displays the mean spectra obtained from PVA-GS-NanoPHE films with increasing concentrations of NanoPHE (C0 to C5).
The mean spectrum of C0 exhibits the spectral features of the films formed from the mixture without NanoPHE. The signature is mainly due to the combined contribution of PVA and Glycerol, most of whose peaks overlap, for example in the spectral ranges 825-855 cm −1 , 920-928 cm −1 , 1060-1128 cm −1 , and 1440-1468 cm −1 . The contribution from AMA is limited to the 1001 cm −1 and 1612 cm −1 features.
The emergence of the strong PHE feature at 1639 cm −1 as a function of NanoPHE concentrations added to the films is particularly noticeable (Figure 5, C1-C5). At first, the PHE band appears as a separate band beside features from AMA ( Figure 5, C1 and C2) but for higher concentrations, the dominance of the PHE band results in a merging with the weaker AMA features ( Figure 5, C3, C4 and C5). Over the range of 300-1500 cm −1 , the spectral variations linked to PHE are not as obvious. However, PCA performed in the 700 and 1500 cm −1 region leads to a clear separation of data clustered as a function of PHE concentrations along the principal component 1 axis (PC1, 68.5% of the total explained variance) ( Figure 6A). The emergence of the strong PHE feature at 1639 cm −1 as a function of NanoPHE concentrations added to the films is particularly noticeable ( Figure 5, C1-C5). At first, the PHE band appears as a separate band beside features from AMA ( Figure 5, C1 and C2) but for higher concentrations, the dominance of the PHE band results in a merging with the weaker AMA features ( Figure 5, C3, C4 and C5). Over the range of 300-1500 cm −1 , the spectral variations linked to PHE are not as obvious. However, PCA performed in the 700 and 1500 cm −1 region leads to a clear separation of data clustered as a function of PHE concentrations along the principal component 1 axis (PC1, 68.5% of the total explained variance) ( Figure 6A).
Interestingly, the loading plot ( Figure 6Bb) highlights spectral differences matching features specific to PHE (Figure 6Ba), which are negative, consistent with the increasing PHE concentration ( Figure 6A). The negative peak at 1001 cm −1 cannot be assigned to PHE but indicates the contributions of AMA (See Figure 5) to the composite spectrum. Positive peaks between 800-950 cm −1 are PVA-and Glycerol-related. It can be observed on the scatter plot that PC 2 accounts for 20% of the explained variance. The PC 2 loading vector is presented in Figure 6Bc. There is no evidence of any positive or negative features linked to PHE. However, the positive peak around 800 cm −1 and the sharp positive band at ~1000 cm −1 can be assigned to AMA ( Figure 4C). The other bands are mixed contributions from PVA and glycerol ( Figure 4A,B).
PCA provides clear evidence that, in the data collected from films, the PHE contribution is not limited to the most intense feature at 1639 cm −1 , but that the presence of the AI at different concentrations is identifiable across the full spectrum. Spectral unmixing approaches like NCLS can be performed using the full range, and the fit can be improved when multiple features from reference spectra of each compound can be identified in the film spectra. Interestingly, the loading plot ( Figure 6(Bb)) highlights spectral differences matching features specific to PHE (Figure 6(Ba)), which are negative, consistent with the increasing PHE concentration ( Figure 6A). The negative peak at 1001 cm −1 cannot be assigned to PHE but indicates the contributions of AMA (See Figure 5) to the composite spectrum. Positive peaks between 800-950 cm −1 are PVA-and Glycerol-related. It can be observed on the scatter plot that PC 2 accounts for 20% of the explained variance. The PC 2 loading vector is presented in Figure 6(Bc). There is no evidence of any positive or negative features linked to PHE. However, the positive peak around 800 cm −1 and the sharp positive band at 1000 cm −1 can be assigned to AMA ( Figure 4C). The other bands are mixed contributions from PVA and glycerol ( Figure 4A,B).
PCA provides clear evidence that, in the data collected from films, the PHE contribution is not limited to the most intense feature at 1639 cm −1 , but that the presence of the AI at different concentrations is identifiable across the full spectrum. Spectral unmixing approaches like NCLS can be performed using the full range, and the fit can be improved when multiple features from reference spectra of each compound can be identified in the film spectra. Figure 7 presents reconstructed false colour maps obtained from NCLS analysis applied to the first replicate (R1) of each PVA-GS film prepared for each concentration (C0 to C5). Rows represent the NCLS scores (spectral abundance fractions) calculated for PHE, PVA, Glycerol and AMA. A score equal to zero indicates no contribution for a given compound, while 1 corresponds to 100%, i.e., a spectrum from the map matching 100% with a pure spectrum from a reference. None of the images displayed NCLS scores higher than~0.5-0.6, and therefore intensities are displayed between 0 and 0.7 for clarity and to help direct comparison of the reconstructed heat maps.  Figure 7 presents reconstructed false colour maps obtained from NCLS analysis applied to the first replicate (R1) of each PVA-GS film prepared for each concentration (C0 to C5). Rows represent the NCLS scores (spectral abundance fractions) calculated for PHE, PVA, Glycerol and AMA. A score equal to zero indicates no contribution for a given compound, while 1 corresponds to 100%, i.e., a spectrum from the map matching 100% with a pure spectrum from a reference. None of the images displayed NCLS scores higher than ~0.5-0.6, and therefore intensities are displayed between 0 and 0.7 for clarity and to help direct comparison of the reconstructed heat maps.

Distribution and Content of PHE in PVA-GS-Based Films
For C0, the NCLS analysis of Figure 7 correctly returns a 0% contribution for PHE, while the PVA, Glycerol, and AMA contributions are ~0.42, 0.53 and 0.10, respectively. The initial NCLS score for PHE for C0 (0.000394) is of great importance to verify that the model constructed using NCLS is reliable. It demonstrates there is no false positive, hence no overestimation in the determination of PHE in the films containing 0%. As previously discussed by Miloudi et al., the number of reference spectra used plays a key role in the specificity of the NCLS [26]. Notably, it is important that the reference spectra account for the matrix (i.e., PVA) in order to ensure that the algorithm builds the model on a robust For C0, the NCLS analysis of Figure 7 correctly returns a 0% contribution for PHE, while the PVA, Glycerol, and AMA contributions are~0.42, 0.53 and 0.10, respectively. The initial NCLS score for PHE for C0 (0.000394) is of great importance to verify that the model constructed using NCLS is reliable. It demonstrates there is no false positive, hence no overestimation in the determination of PHE in the films containing 0%. As previously discussed by Miloudi et al., the number of reference spectra used plays a key role in the specificity of the NCLS [26]. Notably, it is important that the reference spectra account for the matrix (i.e., PVA) in order to ensure that the algorithm builds the model on a robust basis to then estimate other constituents. For the purpose of this study, 1875 spectra collected from three PVA films were used in the NCLS.
Ratios between compounds are clearly seen to evolve systematically across the formulations C1-C5. While PHE maps display a change of colour from 0% contribution (deep blue) in C0 to~45% in C5, reflecting the increase of concentrations in films, PVA and Glycerol have the opposite behaviour, with decreasing NCLS scores from C0 to C5. For AMA, the heat maps do not exhibit significant changes in intensity. basis to then estimate other constituents. For the purpose of this study, 1875 spectra collected from three PVA films were used in the NCLS. Ratios between compounds are clearly seen to evolve systematically across the formulations C1-C5. While PHE maps display a change of colour from 0% contribution (deep blue) in C0 to ~45% in C5, reflecting the increase of concentrations in films, PVA and Glycerol have the opposite behaviour, with decreasing NCLS scores from C0 to C5. For AMA, the heat maps do not exhibit significant changes in intensity.
Irrespective of the PHE concentration tested, the NCLS reconstructed color-coded maps demonstrate the homogeneous distribution of the constituent components of the films. For concentration C5, the mean and standard deviations for NCLS scores are, respectively, 0.433 ± 0.014, 0.351 ± 0.009, 0.419 ± 0.011 and 0.112 ± 0.007 for PHE, PVA, Glycerol and AMA, corresponding to respective relative standard deviations (RSD) of 3.33%, 2.48%, 2.72% and 6.27%. These values confirm that the variations in concentration across the Raman maps of the dried films are minimal. The value, however, is slightly higher for AMA, suggesting a more heterogeneous distribution in the maps. Taking the mean from the three replicates, R1, R2 and R3, the RSD gives ranges from 0.58% (C1) to 1.74% (C5) for PVA and from 0.53% (C0) to 1.83% (C5) for Glycerol. For PHE, the mean RSD values are 4.90%, (C3), 4.43% (C4) and 3.70% (C5). For C0, C1 and C2, RSD values above 10% are reached, which can be explained by the low NCLS scores that progressively decrease towards 0, leading to increased weight of heterogeneity compared the concentration in the films. Nevertheless, it can be concluded that nanoformulation is an efficient way to ensure even distribution of a hydrophobic AI such as PHE in the hydrophilic matrix formed by polymers. Irrespective of the PHE concentration tested, the NCLS reconstructed color-coded maps demonstrate the homogeneous distribution of the constituent components of the films. For concentration C5, the mean and standard deviations for NCLS scores are, respectively, 0.433 ± 0.014, 0.351 ± 0.009, 0.419 ± 0.011 and 0.112 ± 0.007 for PHE, PVA, Glycerol and AMA, corresponding to respective relative standard deviations (RSD) of 3.33%, 2.48%, 2.72% and 6.27%. These values confirm that the variations in concentration across the Raman maps of the dried films are minimal. The value, however, is slightly higher for AMA, suggesting a more heterogeneous distribution in the maps. Taking the mean from the three replicates, R1, R2 and R3, the RSD gives ranges from 0.58% (C1) to 1.74% (C5) for PVA and from 0.53% (C0) to 1.83% (C5) for Glycerol. For PHE, the mean RSD values are 4.90%, (C3), 4.43% (C4) and 3.70% (C5). For C0, C1 and C2, RSD values above 10% are reached, which can be explained by the low NCLS scores that progressively decrease towards 0, leading to increased weight of heterogeneity compared the concentration in the films. Nevertheless, it can be concluded that nanoformulation is an efficient way to ensure even distribution of a hydrophobic AI such as PHE in the hydrophilic matrix formed by polymers.
As regards estimating the respective content of ingredients in the films, NCLS scores cannot be used directly, as the relative spectral abundances are weighted by the respective spectral responses of each compound. Sums of spectral abundance exceed a value of 1, especially when the PHE concentration increases; hence, the correlation of NCLS results with actual constituent concentrations is not evident at first. For instance, in Figure 7 the sum of the NCLS scores corresponding to PHE, PVA, Glycerol and AMA for, respectively, C0, C1, C2, C3, C4 and C5 are 1.11, 1.14, 1.17, 1.26, 1.31 and 1.35. In this study, Raman analysis was performed on dry films; therefore, changing the amount of PHE inevitably modifies the concentration, expressed as %w/w, of all constituents. It is necessary to convert this notion to NCLS scores in order to have a better appreciation of the systematic evolution of the spectral abundance fractions calculated by NCLS as a function of PHE concentrations. Once the percentages of spectral abundance were calculated (see Material and Methods), it was possible to construct linear regression plots (Figure 8).
sum of the NCLS scores corresponding to PHE, PVA, Glycerol and AMA for, respectively, C0, C1, C2, C3, C4 and C5 are 1.11, 1.14, 1.17, 1.26, 1.31 and 1.35. In this study, Raman analysis was performed on dry films; therefore, changing the amount of PHE inevitably modifies the concentration, expressed as %w/w, of all constituents. It is necessary to convert this notion to NCLS scores in order to have a better appreciation of the systematic evolution of the spectral abundance fractions calculated by NCLS as a function of PHE concentrations. Once the percentages of spectral abundance were calculated (see Material and Methods), it was possible to construct linear regression plots (Figure 8). Raman response or Raman cross section dictates the relative strength of the signal originating from the compounds found in the analysed mixture. Depending on the chemical structures involved, the contribution of molecules can significantly differ from one to another. In the spectral map of C0, the relative spectral abundances of PVA and Glycerol are 0.42 and 0.53 respectively (PVA/Glycerol = 0.792). In Table 2, however, their nominal, as prepared %w/w contents are 46.30 and 48.96% (PVA/Glycerol = 0.946), indicating that the spectral contribution of PVA is relatively weaker than that of Glycerol by a factor of ~0.84. As can be seen from the linear regressions in Figure 8, this ratio of relative percentage of spectral abundance is maintained across the films as the concentration of PHE is increased. PVA and Glycerol have relative concentrations that are inversely correlated with PHE content, as also clearly observed in spectral abundance fractions in the NCLS reconstructed images (Figure 7, rows 2 and 3). In practice, Raman spectra are collected from a voxel defining the volume of analysis where the laser is focused. When analysis is performed on dry films, adding PHE inevitably leads to a decrease in PVA, Glycerol and AMA within this voxel. Linear regressions result in good fitting, with R² = 0.9971 (Y = −6.247x + 49.567) and R² = 0.9916 (Y = −4.8209x + 40.314), respectively, for Glycerol and Raman response or Raman cross section dictates the relative strength of the signal originating from the compounds found in the analysed mixture. Depending on the chemical structures involved, the contribution of molecules can significantly differ from one to another. In the spectral map of C0, the relative spectral abundances of PVA and Glycerol are 0.42 and 0.53 respectively (PVA/Glycerol = 0.792). In Table 2, however, their nominal, as prepared %w/w contents are 46.30 and 48.96% (PVA/Glycerol = 0.946), indicating that the spectral contribution of PVA is relatively weaker than that of Glycerol by a factor of~0.84. As can be seen from the linear regressions in Figure 8, this ratio of relative percentage of spectral abundance is maintained across the films as the concentration of PHE is increased. PVA and Glycerol have relative concentrations that are inversely correlated with PHE content, as also clearly observed in spectral abundance fractions in the NCLS reconstructed images (Figure 7, rows 2 and 3). In practice, Raman spectra are collected from a voxel defining the volume of analysis where the laser is focused. When analysis is performed on dry films, adding PHE inevitably leads to a decrease in PVA, Glycerol and AMA within this voxel. Linear regressions result in good fitting, with R 2 = 0.9971 (Y = −6.247x + 49.567) and R 2 = 0.9916 (Y = −4.8209x + 40.314), respectively, for Glycerol and PVA (Figure 8, blue and green). Similarly for AMA, taking into account the relative spectral abundance of 0.102 observed in image C0, it can be determined that its spectral contribution is~2.4 stronger compared to PVA and~2 stronger than Glycerol. AMA, however, has a lower R 2 with 0.6096 (Y = −1.2187x + 11.055). The consistency of the NCLS analysis is evidenced by the fact that the regression slope for PHE, 12.287, is equal to the sum of the slopes for PVA, Glycerol and AMA (Figure 8).
The PHE percentage of spectral abundance of 0.0373% ± 0.008 for C0 increases to a maximum of 33.90 ± 1.21% for C5 (Figure 8, red). The linear regression obtained from the mean values of each concentration delivers a relationship of Y = 12.287x − 0.9352 with R 2 = 0.9952, highlighting the proportionality of the Raman response according to the AI concentration in the dried films. However, while the concentration of PHE ranges from 0% w/w (C0) to 2.86% w/w (C5) ( Table 2), the observed percent spectral abundance is~12 times higher, suggesting the AI has a particularly strong Raman response compared to the other ingredients. This is confirmed by the spectral map of C1 (Figure 7), in which the relative spectral abundances of PHE and PVA are equal to 0.0334 and 0.427 (PHE/PVA = 0.0782), respectively. In Table 2, however, their nominal, as prepared %w/w contents are 0.33% and 47.2% (PHE/PVA = 0.00699), respectively, indicating that the spectral contribution of PHE is significantly stronger than that of PVA, by a factor of~11.2. Similarly, in the spectral map of C1, the relative spectral abundance of Glycerol is 0.514, and hence PHE/Glycerol = 0.0650. In Table 3, however, its nominal, as prepared %w/w content is 47.5%, a ratio of PHE/Glycerol = 0.00694, indicating that the spectral contribution of PHE is also significantly stronger than that of Glycerol, by a factor of~9.37.
The relative spectral abundance of AMA in the spectral map of C1 is 0.139, and hence PHE/AMA = 0.240 for a prepared concentration of 4.8% w/w ( Table 2). The ratio PHE/AMA = 0.06875 suggest that the spectral contribution of PHE is also significantly stronger than that of AMA, by a factor of~3.49.

Absolute Quantification of PHE in PVA-GS Films
NCLS analysis is an unmixing method that delivers results as scores, i.e., abundance fractions, for each compound introduced into the calculation for references. The expression of results as percentages of spectral abundance enables construction of a linear regression for PHE as a function of concentration in the film (Figure 8). Thanks to the proportionality between %w/w PHE in films and percent spectral abundance in the Raman data, a predictive model can be constructed. To demonstrate this, the map of R1 was used for calibration to define the equation of the linear regression (Y = 12.711x − 1.2796, R 2 = 0.9939), and then the percent spectral abundance for all data points from maps R2 and R3 were determined as unknown samples. Figure 9 presents the prepared concentrations plotted against the predicted concentrations. This representation, widely used for spectral data subjected to quantitative data analysis [45][46][47][48], provides a visualisation of the quality of the fitting achieved. With Y = 0.9499x + 0.0406 and R 2 = 0.995, the high correlation between the concentration of the AI in the dried films and the concentrations estimated using the NCLS results is confirmed. The RMSE is a standard way to measure the error of a model in predicting quantitative data. For instance, an RMSE = 0.0868% w/w PHE, that is, 6.07% of the median concentration of the range analysed, confirms the reliability of the analysis. The mean percent relative errors in the predicted concentration is equal to 3.75%, highlighting the accuracy of the analysis. Considering the concentrations individually, it is observed that the percentages of relative errors are, respectively, 1.61 ± 0.79, 13.6 ± 0.15, 3.42 ± 1.81 and 1.92 ± 0.70 for C1, C2, C3, C4 and C5. Except for C2, the concentrations are predicted with relative errors below 5%.  NCLS analysis of Raman microspectroscopic profiles can therefore be employed to perform quantitative analysis of the percentage of AI content within a composite solid film matrix with a high degree of accuracy. The methodology of mapping provides a qualitative analysis of the film in terms of, for example, the degree of homogeneity of the composition, and also provides relatively large datasets to construct and validate calibration models for absolute quantitation purposes.

Conclusions
First, chemical maps reconstructed from Raman maps using NCLS confirm the homogeneous distribution of PHE, the active ingredient, after drying. The results support stabilisation of hydrophobic active cosmetic ingredients such as Pluronic ® F127 stabilised nanoformulation as a promising approach to prepare film-forming cosmetics for topical application. In addition, NCLS analysis enables the distribution of other ingredients to be visualised, specifically thanks to a collection of reference Raman spectra collected from pure compounds. Linear regressions for PVA, Glycerol and PHE displayed R 2 > 0.99, highlighting their high correlation in Raman features with varying amounts of the AI within the films. Ultimately, it is further demonstrated that NCLS scores, i.e., calculated spectral abundance for each ingredient, can be used to perform absolute quantification. The R 2 = 0.995 between the %w/w prepared and %w/w predicted concentrations strengthens the potential of Raman spectroscopy coupled to multivariate analysis as a tool to support the development of new film-forming topical products. Second, this study demonstrates that confocal Raman Imaging is a valuable non-destructive, label-free and reagent-free tool that can support the development of green analytical protocols in the cosmetics industry. In addition to delivering spatially resolved chemical maps with micrometric resolution, this technique can contribute greatly to reducing the consumption of consumables and solvents from the early stages of development of new formulae to final quality control of products. Coupling Raman confocal imaging with advanced multivariate analysis such as NCLS opens perspectives for standardisation and automation of data mining protocols, and thus for its implementation for daily use in the industry.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.