Evaluation of the Miscibility of Novel Cocoa Butter Equivalents by Raman Mapping and Multivariate Curve Resolution–Alternating Least Squares

Cocoa butter (CB) is an ingredient traditionally used in the manufacturing of chocolates, but its availability is decreasing due to its scarcity and high cost. For this reason, other vegetable oils, known as cocoa butter equivalents (CBE), are used to replace CB partially or wholly. In the present work, two Peruvian vegetable oils, coconut oil (CNO) and sacha inchi oil (SIO), are proposed as novel CBEs. Confocal Raman microscopy (CRM) was used for the chemical differentiation and polymorphism of these oils with CB based on their Raman spectra. To analyze their miscibility, two types of blends were prepared: CB with CNO, and CB with SIO. Both were prepared at 5 different concentrations (5%, 15%, 25%, 35%, and 45%). Raman mapping was used to obtain the chemical maps of the blends and analyze their miscibility through distribution maps, histograms and relative standard deviation (RSD). These values were obtained with multivariate curve resolution–alternating least squares. The results show that both vegetable oils are miscible with CB at high concentrations: 45% for CNO and 35% for SIO. At low concentrations, their miscibility decreases. This shows that it is possible to consider these vegetable oils as novel CBEs in the manufacturing of chocolates.


Introduction
In the manufacture of chocolates, one of the main ingredients is cocoa butter (CB), which is the main contributor to the high fat content of chocolate. In total, 30 to 40% of the weight of chocolate is fat [1]. CB remains as a solid at 20 • C (with hard texture and snap) and melts rapidly at 33 • C, leading to the release of flavor and soft mouth-feel texture [2,3]. Cost and technical limitations have increased CB demand, although it is the main lipid phase in chocolate manufacturing. In addition, it is also used in combination with other vegetable oils, such as hydrogenated or partially hydrogenated soybean oil and palm oil [2][3][4]. Therefore, researchers have been searching for cheaper alternatives with similar characteristics to CB [2]. These alternatives should improve the physical properties and bloom resistance and reduce the health risks of the final product [3]. CB can be substituted with vegetable fats by blending to produce chocolate or replace CB either partially or wholly [2]. These are the so-called cocoa butter equivalents (CBE), which

Raman Mapping
Following the process of Mitsutake et al. [22] with some modifications, the samples were mapped using a Raman confocal microscope system (Horiba Scientific, XploRA plus, Montpellier, France). Chemical maps were obtained by a 532 nm laser as an excitation light with a 50% filter. The experimental conditions were as follows: 100 nm slit width, pinhole 100 µm, x50/0.90 NA Vis-LWD air objective, and 1 s acquisition time with 2 accumulations. The Raman signal was obtained using a 600 lines/mm grating centered between 800 and 3100 cm −1 . The acquired spectra were corrected in a range from 1000 to 1800 cm −1 , smoothed, and baseline corrected using LabSpec 6 Suite software. Each sample generated a cube of data with dimensions of 25 × 25 × 761, where 25 was the number of pixels at the x and y axes and 761 was the number of spectral variables.

Data Analysis of Chemical Maps
According to Vajna et al. [29], before chemometric evaluation, all spectra were baseline corrected (this was done by using the same baseline points for all maps and pure component spectra). Then, the spectral range from 1000 to 1800 cm −1 was used for the corresponding evaluation. Raman chemical map data were analyzed by using Solo+MIA software (Eigenvector, Research, Inc. Wenatchee, WA, USA). The raw 3-dimensional data were unfolded into a 2-dimensional matrix. The estimation of pure component spectra from the Raman chemical maps was carried out by MCR-ALS. This technique is based on the following bilinear model (Equation (1)): where X (p * λ) is the matrix containing the mapping spectra, S T (k * λ) is the set of pure component spectra, and C (p * k) contains the vectors of spectral concentrations (each row in C contains the concentrations of the k ingredients). The matrix E represents the residual noise. MCR-ALS generated both the concentration matrix C (scores) and recovered spectrum matrix S T (loadings) from the dataset X in an iterative manner, using an initial estimation for either C or S T and appropriate constraints.

Preprocessing and Constraints
The preprocessing technique was normalized (1-norm, area = 1). The normalization of concentration profiles or resolved spectra or use of reference concentration values within the optimization helps to suppress the rotational ambiguity [30]. The applied constraints were non-negativity, which forces the profiles to be formed by positive values and can be implemented replacing negative values by zeros [30], and equality, which makes the concentration profile and/or spectra of a component equal to a certain known predefined shape [30]. Pure CB, CNO, and SIO spectra were used as equality constraints.

Miscibility
To determinate the miscibility of the vegetable oils proposed as novel CBEs with CB, the homogeneity of the samples was quantitatively and qualitatively determined. The relative standard deviation (RSD) is a commonly used tool in the pharmaceutical industry to estimate the homogeneity of a component within a blend [31] and describe the distribution of the components quantitatively. The RSD was calculated from the ratio of the standard deviation (σ) and the mean (µ) of each measured image score. Using RMarkdown software, the RSD of the scores within a chemical image was calculated. A lower RSD of the chemical image corresponded to a more homogeneous distribution of the respective ingredient [29,32,33] and, therefore, its miscibility. Qualitatively, the homogeneity of the samples was analyzed using histograms. According to Gendrin et al. [33], a histogram showing a symmetric distribution with a narrow base and sharp peak is representative of a low-contrast image and therefore of a homogeneous sample. Conversely, an asymmetric histogram with a wide base and flatter peak or several modes is representative of a more contrasted image, i.e., a heterogeneous sample. Figure 1 shows the characteristics of the Raman spectra of the pure components (CB, CNO, and SIO) in the full range (1000-3100 cm −1 ) at room temperature (20 • C). In the CB spectra, the C-H stretching region shows peaks at 2885.7 and 2886.1 cm −1 which are assigned to alkyl-chain methylene symmetric (ν s ) and antisymmetric (ν as ) stretching. A peak at 2936.5 cm −1 associated with the terminal methyl ν s (CH 3 ) stretch was observed. In the C=O stretching region, we can see 2 peaks at 1745.4 and 1733.8 cm −1 , which are representative of forms III and IV at room temperature ( Figure 2). The full width at half maximum (FWHM) of the peak at 1745.4 cm −1 was 17.38 cm −1 . In the C=C stretching region [ν s (C=C)] of the olefinic band, we can see a peak at 1662.7 cm −1 (Table 2), representing the solid state of CB form IV at room temperature. This peak is of greater intensity in SIO due to its liquid state and is also related to the proportion of oleic acid. A total of 2 peaks at 1445.9 and 1462.9 cm −1 in the CH 2 and CH 3 deformation regions can also be seen. In the CH 2 twisting region, we can also see a peak at 1301.3 cm −1 , which is related to the degree of coupling of the alkyl chains in the lipids. The CB spectra also reveal 3 characteristic peaks of the CB polyforms, located in the C-C stretching range (1030-1183 cm −1 ). The peak at 1102.21 cm −1 indicates the existence of CB in its solid state at room temperature. This peak is not present in CNO and SIO. The pure spectra of CB, CNO, and SIO look similar ( Figure 1); however, we can find some differences in some peaks. CNO has a peak at 1662.7 cm −1 whose intensity is lower. SIO has high intensity peaks at 1276.7 cm −1 and 3020.2 cm −1 that the others do not have. These peaks are assigned to plane =CH deforma- tion in an unconjugated cis C=C and asymmetric C=H stretch group. From Figure 2, we can see that the peak at 1745.3 cm −1 in CB is shown at 1747.93 cm −1 in CNO and at 1746.8 in SIO, and the peak at 1733.84 cm −1 in CB is shown at 1734.83 cm −1 in SIO. We found differences between the area ratios of the peaks at 1733.84 and 1745.43 cm −1 , which can be used as differentiation patterns between CB, CNO, and SIO.  [34], and Jiménez-Sanchidrián et al. [35]. the alkyl chains in the lipids. The CB spectra also reveal 3 characteristic peaks of the CB polyforms, located in the C-C stretching range (1030-1183 cm −1 ). The peak at 1102.21 cm −1 indicates the existence of CB in its solid state at room temperature. This peak is not present in CNO and SIO. The pure spectra of CB, CNO, and SIO look similar ( Figure 1); however, we can find some differences in some peaks. CNO has a peak at 1662.7 cm −1 whose intensity is lower. SIO has high intensity peaks at 1276.7 cm −1 and 3020.2 cm −1 that the others do not have. These peaks are assigned to plane =CH deformation in an unconjugated cis C=C and asymmetric C=H stretch group. From Figure 2, we can see that the peak at 1745.3 cm −1 in CB is shown at 1747.93 cm −1 in CNO and at 1746.8 in SIO, and the peak at 1733.84 cm −1 in CB is shown at 1734.83 cm −1 in SIO. We found differences between the area ratios of the peaks at 1733.84 and 1745.43 cm −1 , which can be used as differentiation patterns between CB, CNO, and SIO.     [34], and Jiménez-Sanchidrián et al. [35].  Table 3 shows the area ratios and FWHM of CB, CNO, and SIO in the mode of vibrations ν(C=O). The lower FWHM values indicate a better arrangement of the crystals and their solid state; therefore, at room temperature, the crystals of CB (17.38 cm −1 ) would have a better arrangement than those of CNO (27.52 cm −1 ), demonstrating its solid state.  Table 3 shows the area ratios and FWHM of CB, CNO, and SIO in the mode of vibrations ν(C=O). The lower FWHM values indicate a better arrangement of the crystals and their solid state; therefore, at room temperature, the crystals of CB (17.38 cm −1 ) would have a better arrangement than those of CNO (27.52 cm −1 ), demonstrating its solid state.  Figure 3 shows the spectral range used for the miscibility analysis between CB, CNO, and SIO. The Raman spectra of the pure components are similar due to the similarity in their chemical composition. However, there are important differences that are considered for the analysis by MCR-ALS. These differences are mainly found in the C-C stretching region (1000-1200 cm −1 ), the C=C stretching [ν s (C=C)], and the carbonyl C=O stretching region (1700-1800 cm −1 ). The peak located at 1662.77 cm −1 has a greater intensity in the SIO than in CB and CNO, and is characteristic of its liquid state.   Figure 3 shows the spectral range used for the miscibility analysis between CB, CNO, and SIO. The Raman spectra of the pure components are similar due to the similarity in their chemical composition. However, there are important differences that are considered for the analysis by MCR-ALS. These differences are mainly found in the C-C stretching region (1000-1200 cm −1 ), the C=C stretching [νs(C=C)], and the carbonyl C=O stretching region (1700-1800 cm −1 ). The peak located at 1662.77 cm −1 has a greater intensity in the SIO than in CB and CNO, and is characteristic of its liquid state.     Table 4 shows the MCR-ALS quality parameters and correlation coefficients between the original spectra and the spectra recovered (S T ) by MCR-ALS. The quality parameter related to the fit of the model was the percentage of explained variance, whose values were between 94.59 and 98.46%, acceptable for our work. The use of the spectra of the pure compounds as equality constraints produced correlation coefficients between 0.9999 and 0.9993. With these values, it can be verified that component 1 was related to CB and component 2 was related to vegetable oils (CNO or SIO) according to the analyzed samples.   Table 4 shows the MCR-ALS quality parameters and correlation coefficients between the original spectra and the spectra recovered ( ) by MCR-ALS. The quality parameter related to the fit of the model was the percentage of explained variance, whose values were between 94.59 and 98.46%, acceptable for our work. The use of the spectra of the pure compounds as equality constraints produced correlation coefficients between 0.9999 and 0.9993. With these values, it can be verified that component 1 was related to CB and component 2 was related to vegetable oils (CNO or SIO) according to the analyzed samples. Table 4. MCR-ALS quality results and correlation coefficients between recovered spectra by the model and real spectra.    Figure 5 shows a comparison between the spectrum recovered (S T ) by the MCR-ALS model and the original spectrum of the pure components. The restrictions used allowed for almost identical spectra with a good correlation. We can note that the spectrum of component 1 recovered by MCR-ALS is identical to the real spectra of CB ( Figure 4a) and component 2 is identical to the real spectra of CNO (Figure 4b).

Numbers of Factor
Foods 2021, 10 Figure 5 shows a comparison between the spectrum recovered ( ) by the MCR-ALS model and the original spectrum of the pure components. The restrictions used allowed for almost identical spectra with a good correlation. We can note that the spectrum of component 1 recovered by MCR-ALS is identical to the real spectra of CB ( Figure 4a) and component 2 is identical to the real spectra of CNO (Figure 4b).    Figure 6 shows the distribution maps and their corresponding histograms. These maps are constructed by the MCR-ALS model from the concentration matrix C and show the distribution of the compounds in the blend. The reddest areas correspond to higher concentrations and the greenest areas correspond to the lowest concentrations. We can note that, in general, the histograms corresponding to the distribution maps are peak-shaped and symmetric, which indicates that both the CNO (Figure 5a-e) and the SIO (Figure 5f-j) form a homogeneous blend with the CB. However, there are differences in the shapes between each histogram caused by the different concentrations of vegetable oil used in each sample. Table 5 shows the quantitative analysis performed on the samples. The lowest RSD value of each component indicates its most homogeneous distribution in the sample. The RSD results show that CNO is more homogeneously distributed in the CB when its concentration is 45%, while its distribution is less homogeneous at 15% or 5%. The distribution of SIO in the sample is more homogeneous at 35% and less homogeneous at 15% or 5%. With these results, we can affirm that the CNO is more miscible with CB than SIO and that the miscibility of both oils improves by increasing their concentrations in the sample.   Table 5 shows the quantitative analysis performed on the samples. The lowest RSD value of each component indicates its most homogeneous distribution in the sample. The RSD results show that CNO is more homogeneously distributed in the CB when its concentration is 45%, while its distribution is less homogeneous at 15% or 5%. The distribution of SIO in the sample is more homogeneous at 35% and less homogeneous at 15% or 5%. With these results, we can affirm that the CNO is more miscible with CB than SIO and that the miscibility of both oils improves by increasing their concentrations in the sample.

Characterization of the Spectra of Cocoa Butter and Vegetable Oils
According to Carmona et al. [36], the spectra were examined separately in the wavenumber region from 1000 to 3100 cm −1 to find differences between them. Then, although the spectra for edible vegetable oils were similar (Figure 1), it could be seen that they exhibit some differences which are small but enable their discrimination [35]. Wang et al. [37] reported a peak at 3016 cm −1 in the Chinese-specific peony seed oil spectrum. This peak is located in the region =C-H stretching vibration of the methyl linoleate group (cis, cis diene) of RCH=CHR, and it is used in the evaluation of oils with different unsaturation degrees. The SIO spectrum has a peak at 3020.2 cm −1 (Figure 1, Table 2) that is not present in CB and CNO. This peak demonstrates the degree of unsaturation of the SIO. In the characterization of the Raman spectra of CB at 22 • C carried out by Bresson et al. [34], 2 peaks were reported at 1744 and 1732 cm −1 , which are representative of forms III and IV. In the present work, these peaks were identified at 1745.4 and 1733.8 cm −1 (Figure 2, Table 2) and show that the existing conformational differences depend on the CB polymorphism. Bresson et al. [38] observed 3 components between 1750 and 1725 cm −1 , 1730, 1735, and 1744 cm −1 , which were attributed to the peak at 1735 cm −1 in CB or the peak at 1736 cm −1 in CBE for form V or VI. This peak was not observed in the CB, CNO, and SIO spectra (Figure 2, Table 3), so we can deduce that the V form was not present in the CB; therefore, the previous statement is corroborated.
In the olefinic band in the C=C stretching range (1200-1800 cm −1 ), Bresson et al. [34] attributed the liquid form of CB to the intensity of the peak located at~1661 cm −1 , as well as to the functional group present in oleic acid. That is, the higher intensity of this peak characterizes the liquid state of CB, and the lower intensity characterizes the solid state (form IV). This coincides with what is observed in Figure 2a,c, in which this peak is located at 1662.7 cm −1 in CB, CNO, and SIO. There is a noticeable difference in the intensity of this peak, since it is higher in the SIO, which determines its liquid state at room temperature and its proportion of oleic acid. Likewise, this peak was observed at 1658 cm −1 by De Géa Neves et al. [18] in the Raman spectrum of CNO and was used as a differentiating pattern between CNO and other vegetable oils. The peak at 1445 cm −1 is related to the C-H deformation vibration, and the peak at 1658 cm −1 is assigned to cis C=C bonds, both of which provide the degree of unsaturation value [35,37]. Therefore, these peaks can be useful to determine the degree of unsaturation of CB, CNO, and SIO.
In the CB spectra, the stretching C-C region (1030-1183 cm −1 ) allows its different polyforms to be identified. According to Bresson et al. [34], the existence of three peaks (1066.26, 1102.21, and 1133.10 cm −1 ) (Figure 3) at room temperature allows us to recognize that it is the IV or V form and the solid state of CB. The physical state of these ingredients must also be taken into account; that is, at room temperature, the semisolid CNO only presented a peak at 1133.10 cm −1 , and the liquid SIO did not present any peak in this region. This characteristic of the SIO spectrum agrees with the spectra of peony seed oil, soybean oil, and extra virgin olive oil reported by Wang et al. [37]. Bresson et al. [34] affirms that the peak at 1100 cm −1 is representative of the solid state of CB and is not shown in the liquid state. This statement agrees with our results, since these peaks are not seen in CNO and SIO, which are not solid at room temperature.
To find differentiation patterns between CB and CBE, Bresson et al. [38] identified peaks at 1744, 1735, and 1730 cm −1 and found notable differences between the area ratios of these peaks. The results in Table 3 indicate that only two peaks were identified at 1733.84 and 1745.43 cm −1 , which correspond to the peaks mentioned above. The area ratios of CB and SIO are very different, making their differentiation possible. It was not possible to calculate the area ratio for the CNO because the peak at 1733.84 cm −1 was not found. FWHM is also an indicator of polymorphism in the sense that a decrease in this value indicates the transition of CB from a liquid to a solid due to the better arrangement of the crystals [34]. This statement agrees with the results of Table 3, since the semisolid state of CNO will produce a higher FWHM than CB, which is solid at room temperature. The same behavior does not occur with SIO, since this vegetable oil is completely liquid at room temperature.

Miscibility of Cocoa Butter and Vegetable Oils
Following Mitsutake et al. [22], to start the analysis by MCR-ALS, the Raman range from 1000 to 1800 cm −1 was chosen, because it contains those peaks that allow differences to be found between the three materials studied. Therefore, Figure 3 shows the range of analysis and the main peaks that differentiate CB from CNO and SIO. The most striking difference is the high intensity of the peak of the SIO spectrum located at 1662.77 cm, which is related to its liquid state at room temperature. On the other hand, De Géa Neves et al. [18] reported the existence of peaks at 1264 and 1658 cm −1 in the CNO spectrum that differentiated it from other vegetable oils. In the present work, these peaks were shown at 1268.8 and 1662.2 cm −1 , and their intensity was very low with respect to CB and SIO. According to Castro et al. [39], to remove noise signals and optimize the results of the MCR-ALS model, the data were preprocessed using the Savitzky-Golay filter from LabSpec 6 and the 1-norm normalization method from Solo + MIA. Figure 4b,d shows the preprocessed data showing low scattering contributions, with which they were corrected, obtaining acceptable results according to Zhang et al. [39,40].
The analysis of mixtures has been a constant concern in any scientific domain [41]. MCR-ALS solves this problem by providing a chemically (scientifically) significant additive bi-linear model of pure contributions from an original data matrix [30]. This bilinear model could produce several solutions, known as rotational ambiguity, which is the primary source of uncertainty. Thus, selecting the appropriate constraints is essential to obtain optimal solutions [30,41,42]. Once the optimization process has been finished, the MCR-ALS results are the set of concentration profiles, spectra and quality parameters (explained variance) related to the model [30]. Therefore, nonnegativity and equality restrictions were applied to our data. With these considerations, the explained variance of the MCR-ALS model for all samples was between 94.59 and 98.46% (Table 4), which means that this model is capable of representing the original data with high precision. Mitsusake et al. [22] reported explained variance percentages between 98.9 and 99.6% in MCR-ALS when it was applied to blends formulated using natural excipients, and Zhang et al. [40] reported values between 99.43 and 99.56% in the analysis of the constituents of commercial chocolate samples. The authors conclude that their results are well adjusted, and therefore, the MCR-ALS model is capable of constructing chemical maps of the samples. Although these values are higher than the results found in our work, we can say that our data fit the MCR-ALS model in such percentages.
Zhang et al. [40] worked with white chocolates, making a comparison between the spectrum of the pure components such as sucrose, lactose, butter and whey. The correlation coefficients between the pure spectra and those recovered by MCR-ALS with data preprocessing were between 0.6701 and 0.9910, which were considered satisfactory. The equality constraint fixes the recovered spectra or concentrations to specific known values [43]. This is the reason why the values of the correlation coefficient between the recovered spectra and the pure compounds are high (Table 4). Additionally, Figure 5 shows that there is no rotational ambiguity, because the recovered spectra (S T ) are identical to the original spectra. The same results were obtained for the other samples. Based on these results, we can affirm that the first spectrum recovered by MCR-ALS (component 1) corresponds to CB, while the second spectrum (component 2) corresponds to vegetable oil, according to the sample analyzed.
The mathematical analysis of each image allows for the extraction of parameters that are helpful in the interpretation of the images and in understanding of the blending process studied [31]. The data provided by Raman mapping contain spectral and spatial information; then, MCR-ALS can be used to visualize the concentration distribution maps of the different components present in a sample based on their individual spectral signals [44]. However, the quality of the Raman mapping is limited by the spectra of the individual compounds and their concentration, so the analysis becomes complex if the spectra of the components have common peaks [21], as is the case with CB, CNO and SIO (Figure 3). This was another reason why we decided to use the spectra of the pure compounds as equality constraints. It is important to mention that the values in matrix C are related to the concentrations, but they are not the real concentrations of the components of the blends, so the mean value should not be compared with the real concentration of each component [22]. Figure 6 shows the distribution maps of CB, CNO, and SIO in all the samples analyzed constructed from the matrix of concentrations C obtained by MCR-ALS. Figure 6a shows a better distribution of CB and CNO at concentrations of 55 and 45%, respectively. Similarly, Figure 6g shows a better distribution of CB and SIO at concentrations of 65 and 35%, respectively. Both figures show better distribution than the others. Similar results were obtained by Scoutaris et al. [21] when analyzing mixtures of paracetamol (PMOL) and compritol 888 (C-888).
We consider that the miscibility of two components can be determined by their distribution in a chemical map, and homogeneous distribution is an indicator of good miscibility. The homogeneity of the samples is also analyzed using the histograms of the distribution map. According to Gendrin et al. [45], a histogram that exhibits a symmetric distribution with a narrow base and a sharp peak is representative of an image with a low contrast, and therefore a homogeneous sample. The histograms corresponding to each distribution map show a symmetric shape in all cases ( Figure 6), with differences according to the actual concentration of each sample. Lyon et al. [46] prepared tablets composed of furosemide and excipients at five different degrees of mixing, reporting that the most homogeneous distribution was obtained in those samples whose histograms were symmetric.
The homogeneity of the samples can be quantitatively analyzed using the RSD. According to Scoutaris et al. [21], RSD has been used to compare the homogeneity of a sample; a low RSD value is interpreted as signifying higher homogeneity. Mitsusake et al. [22] showed that the standard deviation of histograms (STD) is used to evaluate the miscibility for the preformulation stage of pharmaceutical tablets. Therefore, we consider that both parameters are comparable. Lyon et al. [46] used the RSD of the histograms generated by the image scores to determine the homogeneity of the distribution of furosemide in tablets, observing a progressive increase in RSD as the degree of homogeneity decreased. Mitsutake et al. [22] observed an intermediate miscibility (STD = 6.9) between CB and CNO at real concentrations of 75 and 25%, respectively. Similar results were obtained in the present study (Table 5), in which it is shown that blends containing 45% CNO and 35% SIO have the lowest RSD; therefore, they form a more homogeneous blend with CB and are more miscible at those concentrations. In the elaboration of CBEs, the candidate vegetable oil must have an SOS triglyceride concentration similar to that of CB. To achieve this, it undergoes a fractionation process [8]. In the case of CNO and SIO, they were used in their natural state, showing good miscibility with CB; therefore, they would be good candidates to be used as CBE.
Food products are complex mixtures of heterogeneous nature; obtaining chemical and spatial information from them is crucial for food safety and quality control [47]. For this, food detection technologies play a fundamental role [25]; therefore, it is urgent to develop rapid methods of nondestructive analysis to control the quality and safety of food and thus control its circulation in the market [47]. Chemical Raman imaging (CRI), in combination with chemometrics, can provide spectral information and spatial distributions of specific chemicals, analyzing them non-destructively [42,[47][48][49][50]. However, in the food field, only a few investigations on the application of CRI have been reported [50]. CRM allows the application of Raman mapping with MCR-ALS to obtain the chemical characteristics of CB, CNO, and SIO through their Raman spectral fingerprint ( Figure 1) and to identify the miscibility of these three components ( Figure 6, Table 5), which demonstrates the usefulness of this methodology in initiating the development of new products in the chocolate industry. Some authors have also used this methodology, such as Liu et al. [51], who used chemical Raman mapping to study the compatibility between hydroxypropyl methylcellulose (HPMC) and gelatin. It was found that HPMC was easily adapted to form continuous and intermediate phases from the molecular interactions between both components. Mitsutake et al. [52] studied the miscibility and structural changes (polymorphism) in mixtures of natural and synthetic beeswax (BW) with copaiba oil using Raman mapping and MCR-ALS. Structural changes were found in the synthetic BWs, and the miscibility between both BWs with copaiba oil was not significantly different. It was also observed that the differences between the freshly prepared mixtures and those with three months of storage were more significant when the amount of oil was increased. On the other hand, Rodríguez et al. [12] studied the compatibility of shea butter and CB mixtures using the isosolid diagram but did not use Raman mapping.
Lauric acid and partially hydrogenated fats are not recommended in the chocolate industry because they can increase LDL cholesterol levels [53]. Norazlina et al. [2] reported the following CBE candidates: mean fraction of palm oil, mango seed fat, shea stearin, kokum fat, illipe butter, high oleic and stearic sunflower oil, palm stearin, and bambangan kernel fat. However, this author does not report CNO and SIO. The results of the present work show CNO and SIO as possible candidates for novel CBE, as they demonstrate some advantages, such as their high degree of unsaturation, which makes them healthy fats, as well as their excellent molecular compatibility with CB demonstrated by their miscibility. However, it is necessary to carry out some additional studies, such as the investigation of their thermal and rheological properties at different concentrations of solid and liquid lipids and the investigation of their behaviors over time.

Conclusions
In the present work, the usefulness of the confocal Raman microscopy (CRM) technique to identify the chemical properties of cocoa butter, coconut oil and sacha inchi oil is demonstrated. These latter vegetable oils are proposed as candidates to be cocoa butter equivalents in the manufacture of chocolates. The main differences are in the physical state and the degree of unsaturation, which are differentiated by the intensity of the peaks in the Raman spectra. Likewise, the usefulness of the chemometric technique known as multivariate curve resolution-alternating least squares to analyze the miscibility of these vegetable oils with cocoa butter is demonstrated. We conclude that coconut oil is more miscible with cocoa butter at a 45% concentration, and sacha inchi oil is more miscible at a 35% concentration. Between the two vegetable oils, coconut oil is more miscible than sacha inchi oil. We consider that this work is the first step in finding novel CBEs for developing new chocolates. Further work is necessary to evaluate their thermal, rheological, and sensorial properties.