Rheological properties and physical stability of ecological emulsions stabilized by a surfactant derived from cocoa oil and high pressure homogenization

The goal of this work was to investigate the influence of the emulsification method on the rheological properties, droplet size distribution and physical stability of O/W green emulsions formulated with an eco-friendly surfactant derived from cocoa oil. The methodology used can be applied to other emulsions. Polyoxyethylene glycerol esters are non-ionic surfactants obtained from a renewable source which fulfill the environmental and toxicological requirements to be used as eco-friendly emulsifying agents. In the same way, N,NDimethyloctanamide and α-Pinene (solvents used as oil phase) could be considered green solvents. Emulsions with submicron mean diameters and slight shear thinning behavior were obtained regardless of the homogenizer, pressure or number of passes used. All emulsions exhibited destabilization by creaming and a further coalescence process which was applied to the coarse emulsion prepared with a rotor-stator homogenizer. The emulsion obtained with high pressure at 15000 psi and 1-pass was the most stable.


INTRODUCTION
Emulsions are dispersions of two or more immiscible liquid phases, usually water and oil.If the disperse phase is a lipophilic liquid and the continuous phase is a hydrophilic liquid, then it is an oil-inwater (O/W) emulsion.
Emulsions often contain a surface active agent which has two main functions: a) to decrease the interfacial tension between phases enabling easier formation of the emulsion, and b) to stabilize the dispersed-phase against coalescence once it is formed (Muñoz et al., 2007).Surfactants obtained directly from renewable natural materials are gaining more and more attention since they are more eco-friendly than traditional ones (Gómez, 2009;Ruiz-Márquez et al., 2010).Polyoxyethylene glycerol esters derived from cocoa oil are non-ionic surfactants obtained from a renewable source which fulfill the environmental and toxicological requirements to be used as emulsifying agents (Bermejo et al., 2002).These surfactants are completely innocuous for human skin and hair, and their properties are adequate for designing eco-friendly products (Castán and González, 2003;Lutz, 2006;Yu et al., 2008;Denolle et al., 2011).The safety data sheet provided by the supplier reports a value for oral toxicity (LD50) higher than 5000 mg•kg −1 of animals in tests carried out with rats.It is interesting to note that this value would be 3000 mg•kg −1 for salt (Hollinger, 2005).In addition, the interfacial properties at the air/water and α-pinene/water interface of these surfactants, namely equilibrium adsorption, dynamic surface tension and interfacial rheology, have been recently reported (Trujillo-Cayado et al., 2014a;Trujillo-Cayado et al., 2014b).
Regarding the oil phase, it is also important to highlight the trend towards using eco-friendly solvents such as fatty acid dimethylamides and α-Pinene.α-Pinene is a terpenic solvent that can be obtained from the essential oils of juniper, rosemary and eucalyptus (García et al., 2014).This renewable biosolvent is completely miscible with oils and insoluble in water (Bertouche et al., 2013).Fatty acid dimethylamides are among the green solvents found in agrochemical applications (Höfer and Bigorra, 2007).
In order to produce emulsions, the disperse phase must be distributed in finely divided form throughout the continuous phase.Droplet size is an important property because it determines the shelf life stability, consistency and rheological properties of the emulsion (Israelachvili, 1994).In many cases, the aim of emulsification is to produce droplets which are as fine as possible (d<1 μm).Fine emulsions can be produced in many different ways.The most frequent processes for emulsification are the rotorstator systems and the high pressure homogenizers (Schultz et al., 2004).High pressure homogenization based on micro-channel technology (Microfluidics) is the methodology of choice if fluid-like emulsions with submicron mean diameters and narrow droplet size distributions (DSD) are the target, since they can reach extremely high shear rates (Vladisavljevic ´ et al., 2004).Microfluidization processing implies that a coarse emulsion is passed through an interaction chamber using a high pressure pumping device (Jafari et al., 2007a).The interaction chamber consists of two flow channels which are designed to make two coarse emulsion streams collide with each other at a high velocity, thus creating a very high shearing action that provides an exceptionally fine emulsion (Salvia-Trujillo et al., 2014).
The main objective of this work was to study the influence of the emulsification method on the rheological properties, droplet size distribution and physical stability of green O/W emulsions formulated with an eco-friendly surfactant derived from cocoa oil.The methodology resulting from this work can be applied to any kind of emulsion.
A polyoxyethylene glycerol ester derived from cocoa oil, namely Levenol C-201 TM (Glycereth-17 cocoate, HLB=13), received as a gift from KAO, was used as the emulsifier.Deionized water was used for the preparation of all emulsions.

Emulsifi cation process
Oil-in-water emulsions were produced by homogenizing 30 wt% oil (a mixture of oils) with a 70 wt% aqueous phase containing a 3 wt% surfactant.
Firstly, coarse emulsions were prepared with an Ultraturrax T-50/G45F rotor-stator device by applying a rotational speed of 4000 rpm for 120 s at room temperature.Levenol C-201 TM was dispersed in the corresponding amount of water.Then the oils were mixed.Afterwards, the oil phase was added slowly to the aqueous phase while the system was agitated using the aforementioned homogenizer.
The homogenizing chamber was cooled using a cooling jacket containing water at 20 °C in order to minimize the temperature rise.
The nomenclature used to denote the emulsions is as follows: X/Y stands for 30 wt% oil emulsions containing 3 wt% Levenol C-201 TM which have been prepared at X psi of pressure and Y number of passes through the high pressure homogenizer.

Droplet size distribution (DSD)
Droplet size distributions of emulsions aged for 1, 3, 10 and 30 days were measured using a static light scattering instrument (Mastersizer X; Malvern Instruments, UK) after samples had been diluted with water.The refraction and absorption indexes used for the continuous medium were 0.1 and 1, respectively, whereas the refraction index for the dispersed phase was 1.53.The full size distribution was obtained using a poly-disperse analysis.
The mean droplet diameters were expressed as Sauter diameter (D 3,2 ) and volume mean diameter (D 4,3 ): Eq. ( 1) Eq. ( 2) where d i is the droplet diameter, N is the total number of droplets and n i is the number of droplets having a diameter d i .
The "span" was used to determine the distribution width of droplet sizes, which was calculated as follows: where D[v,0.9] and D[v, 0.1] stand for the 90th and 10th percentiles and D[v,0.5] for the median.Determinations were conducted in triplicate.

Rheological measurements
Dynamic viscosities were calculated from flow curves determined at 20±0.1 °C with a controlled stress rheometer (Haake Mars, Thermo-Scientific, Germany), connected to a C5P Phoenix circulator (Thermo-Scientific, Germany).Viscous flow tests were carried out within a stress range of 0.05-2 Pa, using a sandblasted Z20 coaxial cylinder geometry (Ri=1cm, Re/Ri=1.085).This special surface sensor system was used to avoid slip effects.
The samples were loaded into the rheometer measurement cell and allowed to equilibrate for 600s before beginning the test to allow for stress relaxation.
The results represent the mean of at least two measurements taken on emulsions aged for 1, 3, 13, 21 and 40 days.

Multiple light scattering
The destabilization of emulsions was monitored by multiple light scattering using an optical scanning instrument (Turbiscan Lab Expert).This instrument measured the back scattering of monochromatic light (λ=880 nm) from an emulsion as a function of its height.Emulsions were placed into cylindrical glass tubes and stored at room temperature.The back scattering of light from emulsions was then measured with height at different times.The results may be presented either as a backscattering percentage (% BS) or in reference mode (deltabackscattering), i.e. by subtracting the first scan from all the subsequent scans made.

Droplet size distribution
Figure 1 shows the droplet size distributions (DSDs) of the emulsions aged for 24 h as a function of the different emulsification processes.
The emulsions prepared with a high pressure homogenizer showed two populations of droplets regardless of the pressure or the number of passes.The main peak is below 1 micron and the second population is centered at about three microns.The occurrence of the second population of droplets is due to a re-coalescence phenomenon of new droplets as a consequence of some over-processing (Jafari et al., 2007b).During homogenization, the residence time of the surfactant in the chamber of the homogenizer might not be sufficient to allow for its adsorption on the droplet surfaces before droplet-droplet collisions occurred.The high turbulent intensity that is acting on the sample could lead a higher collision rate between drops, thus resulting in a higher coalescence rate (Floury et al., 2004).This means that the extremely high energy per unit volume provided by the high pressure homogenizer during the emulsification process partially breaks the O/W interface formed by glycereth-17 cocoate surfactant, inducing a re-coalescence phenomenon.This is consistent with the fact that the appearance of this second peak occurred only for emulsions processed with the high pressure homogenizer.
The emulsions processed with the Ultraturrax T-50 (coarse emulsion) showed a mono-modal DSD.The mean diameters and span are listed in Table 1 as a function of both the emulsification method and the processing variables.They reveal that the Microfluidizer M110P yielded emulsions with slightly lower mean diameters than the emulsion obtained with the Ultraturrax T50.Also shown is how they were not significantly affected by the homogenization pressure and that two passes through the high pressure homogenizer resulted in worse results, i.e. greater D 3,2 , D 4,3 and span values.This last result may be due to a higher incidence of the re-coalescence phenomenon mentioned above.Interestingly, both devices (a rotor-stator and a high pressure homogenizer based on micro-channel technology) provided emulsions with similar submicron sizes.These results may be interpreted by taking into account the nature of both solvents.AMD-8 TM is a slightly polar solvent and partially soluble in water.This partial solubility is a disadvantage for the development of oil-in-water emulsions because it may lead to a destabilization process by Ostwald ripening (Tadros, 2009).Conversely, α-Pinene is a strongly non-polar solvent with a high interfacial tension (γ=35.5±0.9 mN/m).This may be a disadvantage during the emulsification process since lower interfacial tension results in a higher ability to break into droplets (Pérez-Mosqueda et al., 2014).Thus, by using these solvents alone we obtain emulsions with droplet sizes above 1 μm.However, the use of a mixture of these solvents allows us to prepare very small size emulsions with low energy.According to Santos et al. (2014) an increase in the AMD-10 TM content of an AMD-10/D-Limonene mixture caused a progressive decrease in the interfacial tension.
It should be noted that the addition of just 1 wt% of AMD-10 TM to D-Limonene reduced interfacial tension by 33%.This could be due to the fact that AMD-10 TM was able to migrate to the interface.There, both solvents were distributed according to   a concentration gradient which could be similar to a "core-shell" model (Lee and Lim, 2003) where the AMD-10 TM would tend to stay nearer the interface and D-Limonene in the core of droplets.This distribution makes it easier to create a higher interfacial area by applying the same amount of energy during the emulsification process.Taking into account the similarities between the AMD-10 TM /D-Limonene system and the AMD-8 TM /α-Pinene mixture, a similar situation may be taking place in these emulsions.
The evolution of both the droplet size distribution and the mean diameters with aging time allow for the detection of destabilization phenomena such as the coalescence or Ostwald ripening.A flocculation process cannot be considered since the flocks are disrupted under the action of stirring and pumping during a laser diffraction measurement.By way of example, Figures 2a and 2b show the evolution of DSD and the mean diameters of the emulsion prepared at 15000 psi with 1 pass.All emulsions showed an increase in the D 4,3 with aging time.However, substantial changes in D 3,2 were not observed.Furthermore, both an increase in the second peak and a decrease in the population with smaller size were detected.This usually points to the occurrence of a destabilization process by coalescence instead of an Ostwald ripening phenomenon.For Ostwald ripening the particle size distribution should be shifted towards larger sizes without changing its shape, whereas with coalescence a bi-modal distribution is usually observed (McClements, 2007).If we compare the coarse emulsion with the emulsion obtained at high pressure, the worst behavior in terms of destabilization by coalescence was shown by the coarse emulsion.

Rheological measurements
All the emulsions exhibited shear-thinning behavior, which fitted the power-law equation (R 2 >0.999): where τ is the shear stress, k is the consistency index, γ˙ is the shear rate and n is the power law index.The consistency index has Pa•s n units, which depend on the power law index values.This prevents a direct comparison of this parameter for samples with different values of the power law index.To avoid this, we use a modified power law equation: Eq. ( 5) where τ 1 is the shear stress value at 1s −1 .For shear thinning materials, 0<n<1.A solid material would show n=0, while a Newtonian liquid would show n=1.
Table 2 shows the values of fitting parameters for the emulsions aged for 24 hours as a function of processing variables.Furthermore, the apparent viscosity at a constant stress of 2 Pa (η 2Pa ) is exhibited, which is calculated as follows: All emulsions showed almost Newtonian behavior with a similar flow index, although the emulsion prepared with the Ultra Turrax T50 is closest to this behavior.Two passes through the Microfluidizer M110P gave a lower flow index regardless of the pressure applied.The fluid-like behavior of these emulsions was supported by their low apparent viscosity values (η 2Pa ranged from 6.53 to 7.52 mPa•s at 20 °C).τ 1 and η 2Pa were similar for all emulsions prepared with high pressure but slightly higher than the coarse emulsion.These results are consistent with the occurrence of lower mean diameters (Barnes, 1994).The particle-size effect in the rheology is particularly important for fine dispersions with droplets considerably smaller than 1 μm (Pal, 1996).It should be noted that the fitting parameters of emulsions processed with a high pressure homogenizer showed no significant differences regardless of the homogenization pressure or the number of passes.
Figure 3 shows the shear stress at 1s −1 and flow index values as a function of aging time for the emulsion processed by microfluidization at 15000 psi and one pass, by way of example.All the emulsions exhibited a progressive increase of τ 1 and a decrease of n with the aging time.These results point to the development of a more viscous and shear-thinning structure at the top of the sample container (region where the sample is taken before being placed on the sensor system prior to measurement), which is likely due to a destabilization process by incipient creaming (Santos et al., 2014).As such, the incipient creaming of these emulsions led to a coalescence phenomenon.

Multiple light scattering
Figures 4a and 4b show back scattering (BS) profiles obtained as a function of time for the emulsions prepared with a) Ultraturrax T50 and b) Microfluidizer M110P at 15000 psi and 1pass, by way of example.These figures also include an inset where BS is plotted at the bottom part of the tube with time in the reference mode (delta-back scattering) to obtain a higher resolution of backs cattering changes.Both emulsions exhibited an increase in the size of the serum layer that displays a decrease in BS at the bottom part of the measurement cell.This result indicates the occurrence of a migration of droplets to the upper part of the cell due to a difference in density between the organic and aqueous phases, i.e., a slight clarification due to an incipient creaming in this zone of the measurement cell (García et al., 2014;Mengual et al., 1999;Santos et al., 2013).Furthermore, the clear decrease in back scattering at the top of the tube (higher decrease for coarse emulsion than for the one prepared by high pressure) indicates a further destabilization process by coalescence, which is less significant than the creaming phenomenon.
The extent of creaming can be characterized by using the creaming index (CI) (McClements, 2007): Eq. ( 7) where H E is the total height of the emulsion and H S is the height of the serum layer.
In Figure 5, the creaming index (CI) is plotted as a function of aging time for all emulsions.Obviously, the value of the creaming index will depend on the time that the measurement is made.Normally, the creaming index should start at zero and increase   Eq. ( 8) It was found that the creaming velocity was very similar for all systems prepared with the high-pressure homogenizer.However, an important change in this parameter was found for the coarse emulsion, which exhibited the highest value.This is consistent with the previous results since the creaming rate is influenced by the droplet size and the viscosity of emulsions.Thus, an emulsion with low viscosity and high mean diameters will be more unstable against creaming (McClements, 2007).Figure 6 illustrates the aspect of the emulsions shown in Figure 4 (coarse emulsion (0/0) and the emulsion prepared by high-pressure homogenization at 15000 psi and 1pass (15000/1)) after an aging time of 40 days.It is important to note that destabilization by coalescence could only be observed by the naked eye for the emulsion prepared with the rotorstator device (upper part of the measuring cell).By contrast, destabilization by creaming could be visually observed for all emulsions.

CONCLUSIONS
Regardless of the emulsification method used, the use of a mixture of green solvents and a polyoxyethylene glycerol fatty acid ester derived from cocoa oil as emulsifier led to the creation of emulsions with submicron mean droplet diameters.The high pressure homogenizer achieved a higher fraction of smaller droplets.However, all emulsions prepared by microfluidization showed bimodal distributions.The occurrence of a second peak has been related with a re-coalescence phenomenon induced by an excess of mechanical energy during the emulsifying process.After a short time, the influence of the homogenization pressure on the droplet size distribution of emulsions was not significant, but a second pass through the Microfluidizer led to a slight increase in the mean diameters and span.The ageing time also slightly increases the mean diameters, particularly the volume mean diameter, related to a coalescence process.This process was more significant for the coarse emulsion.
Flow curves of all emulsions exhibited almost Newtonian behavior.The emulsions obtained with the high pressure homogenizer showed higher apparent viscosity and lower flow index values than the coarse emulsion.This result is in agreement with lower mean diameters.
The Multiple light scattering technique allowed for the detection of creaming destabilization and a further coalescence process, which was more significant for the emulsion obtained with the rotor stator homogenizer.The processing variables (pressure and number of passes) did not influence the physical stability of emulsions.
The most stable emulsion was the one prepared with the high-pressure homogenizer at 15000 psi and 1 pass.

FIGURE 1 .
FIGURE 1. Droplet size distributions for emulsions aged 24 hours as a function of the emulsification variables.

FIGURE 2 .
FIGURE 2. (A) Droplet size distributions and (B) Sauter mean diameter and volume mean diameter for the emulsion prepared with Microfluidizer M110P at 15000 psi and 1 pass as a function of aging time.

FIGURE 3 .
FIGURE 3. Shear stress values at 1 s −1 and flow index for the emulsion prepared with Microfluidizer M110P at 15000 psi and 1 pass as a function of aging time.

FIGURE 4 .
FIGURE 4. Changes in Back scattering profiles as a function of the tube height with the storage time in quiescent conditions for (A) the coarse emulsion and (B) the emulsion prepared with Microfluidizer M110P at 15000 psi and 1 pass.Insets illustrate delta back scattering values at the bottom of the measuring cell.Tube length: 55 mm, temperature: 20 °C.

FIGURE 5 .
FIGURE 5. Creaming index (CI) as a function of aging time for selected emulsions.Samples kept under storage at 20 °C.

TABLE 2 .
Flow curve fitting parameters for the power law model for the studied emulsions as a function of emulsification method.Aging time of 24 h