Effect of Biodegradable Hydrophilic and Hydrophobic Emulsifiers on the Oleogels Containing Sunflower Wax and Sunflower Oil

The use of an appropriate oleogelator in the structuring of vegetable oil is a crucial point of consideration. Sunflower wax (SFW) is used as an oleogelator and displays an excellent potential to bind vegetable oils. The current study aimed to look for the effects of hydrophobic (SPAN-80) and hydrophilic (TWEEN-80) emulsifiers on the oleogels prepared using SFW and sunflower oil (SO). The biodegradability and all formulations showed globular crystals on their surface that varied in size and number. Wax ester, being the most abundant component of SFW, was found to produce fibrous and needle-like entanglements capable of binding more than 99% of SO. The formulations containing 3 mg of liquid emulsifiers in 20 g of oleogels showed better mechanical properties such as spreadability and lower firmness than the other tested concentrations. Although the FTIR spectra of all the formulations were similar, which indicated not much variation in the molecular interactions, XRD diffractograms confirmed the presence of β′ form of fat crystals. Further, the mentioned formulations also showed larger average crystallite sizes, which was supported by slow gelation kinetics. A characteristic melting point (Tm~60 °C) of triglyceride was visualized through DSC thermograms. However, a higher melting point in the case of few formulations suggests the possibility of even a stable β polymorph. The formed oleogels indicated the significant contribution of diffusion for curcumin release. Altogether, the use of SFW and SO oleogels with modified properties using biodegradable emulsifiers can be beneficial in replacing saturated fats and fat-derived products.


Introduction
Solid fats/triglycerides are common in the bakery industry due to their unique capacity to provide texture, flavor, and aroma to baked products. Saturated and trans-fats are usually rich in solid fat. As the apprehension regarding the use of trans and saturated fats on human health is rising worldwide, the World Health Organization (WHO) and Food and Agriculture Organization (FAO) are endorsing to increase the amount of mono and polyunsaturated fatty acids in food products [1]. United States Food and Drug Administration (US-FDA) in 2015 recognized trans-fats as unsafe for human consumption. Accordingly, the use of trans-fats were wholly banned for use in food products [2]. The search for solid fat replacers bought scientists closer to the concept of oleogels. Oleogels have received remarkable achievements in medicine, cosmetics, food, and biomaterials [3].

Visual Appearance and Oil Binding Capacity
The SO, SFW, and liquid emulsifier mixture were heated at 80 °C in the water bat to form a clear homogenous solution (Figure 1a). Oleogels were formed after placing them for 30 min in the thermal cabinet, as discussed above (Figure 1b). This was followed b inverting the sample bottle and maintaining it at room temperature for ~30 min to observ the visual appearance and confirm the gel formation. All the formulations used in ou study could form compact structures that were self-standing, which is similar to prev ously reported work [19]. The visual appearance did not significantly differ among all th formulations. Though SPAN-80 formulation appeared identical to the control, the inclu sion of TWEEN-80 in the oleogels may slightly improve the formulations' whiteness. Th formed oleogels were smooth with easier spreadability when touched with hands. [S1, S3, S5 and S10: 1, 3, 5 and 10 mg of SPAN-80 in oleogels of SO containing 5% (w/w) of SFW; T1, T3, T5 and T10: 1, 3, 5 and 10 mg of TWEEN-80 in oleogels of SO containing 5% (w/w) of SFW]. The oleogels were also evaluated for their capacity to bind oil. The SFW was proven to be efficient in binding the SO, evident from the %OBC values. The %OBC was greater than 99% in all the cases. Previous work has reported a stronger mechanical strength of oleogels because of the higher OBC [20]. The addition of both the emulsifiers at the selected amount did not significantly affect the oil binding capacity of the oleoegels. Interestingly, the crystalline phase formed from a 5% wax concentration was sufficient to develop a network to hold sunflower oil into the oleogel. A similar wax concentration has been reported to structure soyabean oil, which is commonly used to replace saturated fats in food products [10]. Waxes containing a long chain of alkanes and esters are much more efficient in vegetable oil gel than waxes made of short-chain alkanes [21]. Many natural waxes such as carnauba wax, candelilla wax, and beeswax, etc., have been reported as an efficient oleogelator because of their ability to crystallize the vegetable oil at lower concentrations [22]. The lower wax concentration for formulating oleogels is attributed to the presence of long-chain esters in the SFW. The emulsifiers in the proposed concentrations did not significantly affect the oil binding capacity of the oleogels.

Colorimetric Analysis
The L*, a*, and b* or CIE Lab is defined by Commission Internationale de l'Eclairag and is widely accepted for color measurement in food products [23]. According to this, all the colors are a combination of red, green, and blue, whose receptors are present in the human eye. The luminance (L*) is referred to as the lightness component, ranging from 0 to 100. Since all our formulations had shown an L* value close to 100, it was inferred that oleogels were substantially luminous. This can be possible due to the presence of smaller fat structures on the surface of oleogels that could reflect most of the light [24].
Additionally, a* and b* ( Figure 2) are the chromatic components that range from green to red and blue to yellow, respectively [25]. The values of these chromatic components range from −120 to +120 [25]. For the values of a* the range follows −ve (red) and +ve (green). For all the formulations, a* value was found negative. This gives an idea about the presence of a better fraction of green hue. Similar results have been reported with oleogels formed from SFW in the previous work [26]. The addition of emulsifiers has shown more shift towards green as compared to the control. The addition of SPAN-80 showed a subsequent increase from the control in the a* value till S3. However, the extent of greenness was similar in S1 and S3. On increasing the emulsifier content in S5, a significantly higher a* value from the control was observed. This value was considerably lower than the a* value of S1 and S3. At the maximum SPAN-80 content, the extent of the greenness of the formulation S10 was similar to the control. Further, the addition of TWEEN-80 showed a subsequent increase from control in the a* value till T5. However, till T5, the rise in the green hue of T5 was only significantly different from T1. Formulation T10 showed a similar a* value as the control, T1, and T5. However, the reduction in the a* of T10 from T3 was significant. Among the formulations, T3 showed the highest share of green, followed by S3.
Further, the b* value in all the formulations except T1 appeared positive. The b* values display −ve (blue) to +ve (yellow). Thus, a positive b* value indicates a more significant proportion of yellow. Among the SPAN-80 formulations, the b* value followed the same trend as that of a*. However, the addition of TWEEN-80 in T1 showed a similar b* value as the control. A significant increase in the b* value of T3 from control and T1 was observed. A further rise in the TWEEN-80 content in T5 displayed a substantial surge in the yellowness from the control and T1. On comparing T5 with T3, there was no difference in the b* values. At the highest emulsifier content in T10, the yellowness of the formulation appeared similar to the control, T1 and T5. However, the reduced b* value in T10 from T3 was significant. Similar to a* value, T3 has shown the highest share of yellow, followed by S3. The absolute colour difference (∆E) is a numerical value obtained from L*, a*, and b* and is generally used to compare the samples with a specific standard (in our case, the control). The calculated values of ∆E are represented in Table S1. Gels 2021, 7, x FOR PEER REVIEW 5 of 29 was significant. Similar to a* value, T3 has shown the highest share of yellow, followed by S3. The absolute colour difference (ΔE) is a numerical value obtained from L*, a*, and b* and is generally used to compare the samples with a specific standard (in our case, the control). The calculated values of ΔE are represented in Table S1.  (Table 5).

Surface Topography
The surface topology images display a uniform distribution of fat crystals in all the formed oleogels ( Figure 3). Obtained topographs clearly show the semi-crystalline structure of oleogels, which has both crystalline and amorphous regions. The fat crystals appeared as globular structures. The fat crystals varied in number and size in all the formulations. For example, the inclusion of SPAN-80 increased the number and size of the crystals in S1 and S3 compared to control. The crystals appeared bright and were most prominent in S3. In S5, the size of fat crystals was reduced; however, S10 showed similar globular size and density to the control. Similarly, an increase in TWEEN-80 concentration in oleogels displayed a rise in crystal structure and distribution size. The increase in the size of the crystals can be due to the co-crystallization of these emulsifiers with triacylglycerol. The phenomenon of co-crystallization may have enhanced the crystal growth [27]. The overall surface topology appears as the dispersion of many small fat crystals, which is the desired property for food applications [28].  (Table 5).

Surface Topography
The surface topology images display a uniform distribution of fat crystals in all the formed oleogels ( Figure 3). Obtained topographs clearly show the semi-crystalline structure of oleogels, which has both crystalline and amorphous regions. The fat crystals appeared as globular structures. The fat crystals varied in number and size in all the formulations. For example, the inclusion of SPAN-80 increased the number and size of the crystals in S1 and S3 compared to control. The crystals appeared bright and were most prominent in S3. In S5, the size of fat crystals was reduced; however, S10 showed similar globular size and density to the control. Similarly, an increase in TWEEN-80 concentration in oleogels displayed a rise in crystal structure and distribution size. The increase in the size of the crystals can be due to the co-crystallization of these emulsifiers with triacylglycerol. The phenomenon of co-crystallization may have enhanced the crystal growth [27]. The overall surface topology appears as the dispersion of many small fat crystals, which is the desired property for food applications [28].

Microstructure Analysis
The micrographs obtained using bright field and polarized light microscopes were used to visualize the network formed by fats in the oleogel. Since the formation of the wax-based oleogel depends upon the entrapment of the oil phase through wax crystals, looking at the morphology and polymorphism of these crystals is an interesting way to comprehend the physical properties of the oleogel. The bright-field micrographs showed a fibrous network ( Figure 4) in oleogels following previously reported work [26]. The reason here can be the high content of wax esters that are present in the SFW. Waxes or any low molecular weight oleogelators are known to form a three-dimensional-fibrous network. This is usually possible by molecular self-assembly of gelator molecules, which is governed by hydrogen bonding, π-π stacking, van der Waals forces, and hydrophobic forces. The molecular units of waxes are primarily linear. Hence, the crystals either grows in 1-dimension or 2-dimensions, thereby forming needle or plate-like structures, respectively [29]. The appearance of the fibre-like crystals results from the unidirectional growth of crystals from a single nucleation point [30]. A continuous branching was observed from various points of the fibres, which is typical for physical gels [31]. The selfassembled fibrillar network often forms these branches through crystallographic mismatch branching (CMB). The mechanism involves nucleation on the surface of the parent fibre of a daughter fibre which appears as a fork [32]. The fibre architecture of S1 seemed to be similar Gels 2021, 7, 133 6 of 28 to control in terms of fibre length, thickness, and branching. However, a slight increase in the emulsifier's amount in S3 increased the fibre length and hyper-branching. A high aspect ratio (length to diameter ratio) of fibres and their network is efficient in entrapping most of the oil [33]. A further rise in the amount SPAN-80 in S5 and S10 reduced the fibre length and branching. Similar to S1, the inclusion of TWEEN-80 in T1 did not significantly alter the fat crystal network's architecture compared to the control. However, the number of fibres, in this case, has increased slightly. The fibres appeared thinner in T3 and T5, which again became thicker in T10. Another interesting observation was regarding the hyper-branching of the observed fibres, which was prominent in T1. After that, the hyper-branching of the fat crystals was reduced with the increase in the TWEEN-80 till T5. Nevertheless, in T10, the fibres suddenly started appearing longer as well as being hyperbranched.  (Table 5).

Microstructure Analysis
The micrographs obtained using bright field and polarized light microscopes were used to visualize the network formed by fats in the oleogel. Since the formation of the wax-based oleogel depends upon the entrapment of the oil phase through wax crystals, looking at the morphology and polymorphism of these crystals is an interesting way to comprehend the physical properties of the oleogel. The bright-field micrographs showed a fibrous network ( Figure 4) in oleogels following previously reported work [26]. The reason here can be the high content of wax esters that are present in the SFW. Waxes or any low molecular weight oleogelators are known to form a three-dimensional-fibrous network. This is usually possible by molecular self-assembly of gelator molecules, which is governed by hydrogen bonding, π-π stacking, van der Waals forces, and hydrophobic forces. The molecular units of waxes are primarily linear. Hence, the crystals either grows in 1-dimension or 2-dimensions, thereby forming needle or plate-like structures, respectively [29]. The appearance of the fibre-like crystals results from the unidirectional growth of crystals from a single nucleation point [30]. A continuous branching was observed from various points of the fibres, which is typical for physical gels [31]. The self-assembled fibrillar network often forms these branches through crystallographic mismatch branching (CMB). The mechanism involves nucleation on the surface of the parent fibre of a daughter fibre which appears as a fork [32]. The fibre architecture of S1 seemed to be similar to control in terms of fibre length, thickness, and branching. However, a slight increase in the emulsifier's amount in S3 increased the fibre length and hyper-branching. A high aspect ratio (length to diameter ratio) of fibres and their network is efficient in entrapping most of the oil [33]. A further rise in the amount SPAN-80 in S5 and S10 reduced the fibre length and branching. Similar to S1, the inclusion of TWEEN-80 in T1 did not significantly alter the fat crystal network's architecture compared to the control. However, the number  (Table 5).   (Table 5).
Among the two types of liquid emulsifiers selected for this study, the SPAN-80 formulation displayed denser hyper-branching compared to TWEEN-80 formulations. This supports better fractal characteristics of lipid crystals in the SPAN-80 included oleogels [34]. The term fractal is associated with multiple crystal fractions, where each fraction holds a similar function to that of the parent crystal. Fractal dimension correlates with gels  (Table 5). Among the two types of liquid emulsifiers selected for this study, the SPAN-80 formulation displayed denser hyper-branching compared to TWEEN-80 formulations. This supports better fractal characteristics of lipid crystals in the SPAN-80 included oleogels [34]. The term fractal is associated with multiple crystal fractions, where each fraction holds a similar function to that of the parent crystal. Fractal dimension correlates with gels displaying an abundant number of crystals that are homogeneously distributed. These are usually calculated through computational methods such as the box-counting method [35]. The fractal dimension is associated with better oil binding capacity, an overall increase in the area for oil adsorption, and a consequent decrease in the void size of the 3-dimensional network of fat crystals [36].
The polarized light micrographs can be regarded as a better way to visualize the crystalline fat structures. The formed crystals appear bright, and the amorphous region appears as a black background. The obtained micrographs from polarized light microscopy further confirmed a mesh-like frame made up of needle/fibre morphology in all the formulations ( Figure 5). The formation of the mesh or needle network can occur through the entanglement of the fat crystals. This morphology results from arrangements made by n-alkanes and wax esters that form a needle network in the structured oil [37]. The van der Waal interaction between the hydrocarbon chain and ester is responsible for the molecular organization and lateral crystal growth to form needle morphology in the oleogels [36]. Specifically, monoglycerides present in SFW can be accountable for crystallizing the fatty acid chain through supercooling [38]. The existence of differences in melting and crystallization temperatures is referred to as the supercooling effect. Some researchers have reported the needle-like structure to be an artifact of platelet morphology [39]. In some formulations, it seems that the fat crystals appeared as a bundle of several fibres (white arrow). In these structures, the needles were arranged together in several sized arrays that created bright visible zones. The needle length is distinctly long in the control sample, which appears to exist on the top of solid architecture. This information was not evident through the bright field images. The possible explanation for such an appearance can be the crystalline region made up of many fibres. On adding SPAN-80, the crystalline regions were dominantly observed in S1, S3, and S10, except for S5, where the amorphous area (or less crystalline area) is distinctly visible (yellow arrow). This can be the role of critical concentration in S5. Another observation in S3 was made for the presence of large globular shaped amorphous regions (pink arrow), which was not majorly seen in any formulation. Similarly, the addition of TWEEN-80 expressed a few amorphous regions in T3, which were rare in T1, T5, and T10, where the crystalline areas were dominant and spread throughout the matrix.

Stress Relaxation Study
The viscoelastic properties of the prepared oleogel were studied using the stress relaxation (SR) profiles ( Figure 6a). Once the oleogel is externally deformed, the wax network and the fluid pressure of the entrapped oil together exert the stress on the probe, which is displayed in the SR profile [40]. The firmness of the formulation is predicted using the maximum attained force (F 0 ). Since the strained condition was maintained for 60 s, after achieving the maximum force, there occurs a decrease in the force values with an increase in time. On inclusion of SPAN-80, the F 0 values of S1 showed similarity to the control. However, a further rise in emulsifier content in S3 caused a marked reduction of the F 0 value of S1. The fragile nature of S3 may explain this observation. A further increase of SPAN-80 did not affect the firmness of oleogels and appeared similar to the control. Among the SPAN-80 formulations, S3, S5, and S10 displayed equal firmness to S1. Additionally, the F 0 values of S3 and S5 were also identical to one another. However, the increased firmness in S10 was statistically significant from both S3 and S5. The increased firmness in S10 can be due to the increased linkage points within the fibrous network. These linkage points could be modulated by the SPAN-80 molecules [41]. were rare in T1, T5, and T10, where the crystalline areas were dominant and spread throughout the matrix.  (Table 5).

Stress Relaxation Study
The viscoelastic properties of the prepared oleogel were studied using the stress relaxation (SR) profiles ( Figure 6a). Once the oleogel is externally deformed, the wax network and the fluid pressure of the entrapped oil together exert the stress on the probe, which is displayed in the SR profile [40]. The firmness of the formulation is predicted using the maximum attained force (F0). Since the strained condition was maintained for 60 s, after achieving the maximum force, there occurs a decrease in the force values with an increase in time. On inclusion of SPAN-80, the F0 values of S1 showed similarity to the control. However, a further rise in emulsifier content in S3 caused a marked reduction of the F0 value of S1. The fragile nature of S3 may explain this observation. A further increase of SPAN-80 did not affect the firmness of oleogels and appeared similar to the control. Among the SPAN-80 formulations, S3, S5, and S10 displayed equal firmness to S1. Additionally, the F0 values of S3 and S5 were also identical to one another. However, the increased firmness in S10 was statistically significant from both S3 and S5. The increased firmness in S10 can be due to the increased linkage points within the fibrous network. These linkage points could be modulated by the SPAN-80 molecules [41].
Similar to the SPAN-80 formulation, only T3 among TWEEN-80 formulations showed a significant reduction in the F0 value from the control. The rest all the formulations hold comparable firmness to control. Among TWEEN-80 formulations T5 and T10 have similar firmness to that of the control. However, the decrease in firmness value of T3 from both T5 and T10 was notable, wherein T5 and T10 had identical firmness.
It was observed that the reduction in the F0 values of S3 and T3 was significant compared to the control. This observation advocates that adding 3 mg of SPAN-80 and TWEEN-80 in the said oleogel acts as a critical concentration point where the oleogels become more fragile. The polarized light micrographs can explain the reduced firmness, which displayed more amorphous regions ( Figure 5) in these formulations. The micrographs can expound the poor mechanical strength in the said formulations [42]. A denser crystalline network of oleogelator provides better mechanical strength and improves the  (Table 5).
ability to hold the oil phase. A lower firmness of S3 and T3 oleogels also suggests their better ability to spread [43]. A detailed comparison of F0 values in all the prepared oleogels is represented in Figure 6b.
The force value at the end of the relaxation profile depicts the residual elastic force (F60), which is described in Figure 6c. This force decay occurs by the significant rearrangements in the gelator molecules, disturbance in network structure, and the break in the fibre network [44]. No significant difference was observed in the F60 values among the formulations. The possible insignificant rise of F60 value from control in both types of formulations is justified by the ring structure at the carbonyl group in sorbitan, which offers rigidity to the 3-dimensional network of fat [30]. Further, the %SR was calculated using the values of F0 and F60 (Equation (7)). The %SR value represents a sample's ability to absorb the energy during a strained condition. The values of %SR follow a different trend in SPAN-80 and TWEEN-80 formulation. The inclusion of SPAN-80 showed no effect on %SR of S1 from the control. However, a significant decrease from the control was observed in the %SR of S3 and S5. S3 showed the lowest %SR value among the SPAN-80 formulation. A lower %SR of the SPAN-80 containing oleogels suggests an increased rigidity to that of the control. Again, the %SR in S10 was observed to be similar to that of the control. Among the SPAN-80 formulations, not many alterations were observed in the %SR values except for a noteworthy increase from S3 to S10. Interestingly, T1 and T3 showed a significant decrease of %SR from control. However, SR behaviour was similar to the control at the higher TWEEN-80 content, i.e., T5 and T10. A trend of T5 ≈ T10 ≈ T1 > T3 was observed for %SR of TWEEN-80 formulations. The inclusion of both emulsifiers at lower content may reduce the reorganizational capability of gelator molecules. This significant decrease also means a high elastic component, as confirmed through FR values and low viscous components of the oleogels [45]. It is quite evident that even a lower amount of emulsifier affects the mechanical properties of the formulation [41]. The possible reason behind this can be the alterations observed through the micrographs in the size, morphology, and number of fat crystals.  (Table 5).
Similar to the SPAN-80 formulation, only T3 among TWEEN-80 formulations showed a significant reduction in the F 0 value from the control. The rest all the formulations hold Gels 2021, 7, 133 9 of 28 comparable firmness to control. Among TWEEN-80 formulations T5 and T10 have similar firmness to that of the control. However, the decrease in firmness value of T3 from both T5 and T10 was notable, wherein T5 and T10 had identical firmness.
It was observed that the reduction in the F 0 values of S3 and T3 was significant compared to the control. This observation advocates that adding 3 mg of SPAN-80 and TWEEN-80 in the said oleogel acts as a critical concentration point where the oleogels become more fragile. The polarized light micrographs can explain the reduced firmness, which displayed more amorphous regions ( Figure 5) in these formulations. The micrographs can expound the poor mechanical strength in the said formulations [42]. A denser crystalline network of oleogelator provides better mechanical strength and improves the ability to hold the oil phase. A lower firmness of S3 and T3 oleogels also suggests their better ability to spread [43]. A detailed comparison of F 0 values in all the prepared oleogels is represented in Figure 6b.
The force value at the end of the relaxation profile depicts the residual elastic force (F 60 ), which is described in Figure 6c. This force decay occurs by the significant rearrangements in the gelator molecules, disturbance in network structure, and the break in the fibre network [44]. No significant difference was observed in the F 60 values among the formulations. The possible insignificant rise of F 60 value from control in both types of formulations is justified by the ring structure at the carbonyl group in sorbitan, which offers rigidity to the 3-dimensional network of fat [30]. Further, the %SR was calculated using the values of F 0 and F 60 (Equation (7)). The %SR value represents a sample's ability to absorb the energy during a strained condition. The values of %SR follow a different trend in SPAN-80 and TWEEN-80 formulation. The inclusion of SPAN-80 showed no effect on %SR of S1 from the control. However, a significant decrease from the control was observed in the %SR of S3 and S5. S3 showed the lowest %SR value among the SPAN-80 formulation. A lower %SR of the SPAN-80 containing oleogels suggests an increased rigidity to that of the control. Again, the %SR in S10 was observed to be similar to that of the control. Among the SPAN-80 formulations, not many alterations were observed in the %SR values except for a noteworthy increase from S3 to S10. Interestingly, T1 and T3 showed a significant decrease of %SR from control. However, SR behaviour was similar to the control at the higher TWEEN-80 content, i.e., T5 and T10. A trend of T5 ≈ T10 ≈ T1 > T3 was observed for %SR of TWEEN-80 formulations. The inclusion of both emulsifiers at lower content may reduce the reorganizational capability of gelator molecules. This significant decrease also means a high elastic component, as confirmed through F R values and low viscous components of the oleogels [45]. It is quite evident that even a lower amount of emulsifier affects the mechanical properties of the formulation [41]. The possible reason behind this can be the alterations observed through the micrographs in the size, morphology, and number of fat crystals.
Further modelling of the SR profile of oleogels was done using Weichert's model. The model includes many Maxwell units (spring and dashpot) connected in parallel [46]. Weichert's model considers relaxation to occur in consequent series. Evaluating the model's usefulness, in our study, we have used this model to analyse the viscoelastic property of the formed oleogels with two characteristic times represented in Equation (1). The correlation coefficient (R 2 ) value among the experimental and model data in all the oleogels formulation was >0.99.
where P, P 0 , P 1 , and P 2 are spring elements, and τ 1 and τ 2 are time constants of dashpots. The P 0 value helps interpret residual force, mechanical stability, and inherent elastic properties of the oleogel at the end of the stress relaxation process (Table 1) [44]. In SPAN-80 formulations, the P 0 value of all formulations appeared similar to the control. However, the rise of P 0 in S5 from S1 and S3 was found significant. The addition of TWEEN-80 caused a substantial and subsequent increment in the P 0 value from control till T3. The increment of P 0 from T1 to T3 was found to be significant. A further rise in the emulsifier content did not affect the P 0 value. The P 0 value of T3 was the highest among all the formulations, which suggests a better inherent mechanical property in this case ( Table 2). The time constants, τ 1 and τ 2 , provide information regarding instantaneous relaxation time and delayed relaxation time. When stress is applied to the gel, the immediate relaxation time corresponds to the molecular rearrangement of the constituents of stressed oleogels. On adding SPAN-80 to oleogels, the value of τ 1 was significantly reduced from the control in S1 and S3. A further rise in emulsifier content in S5 caused a significant increase in the τ 1 from control, S1, and S3. This observation suggests slower molecular rearrangement inside S5 oleogels when the force is applied. The possible reason for this is the heterogeneous networking, as shown in polarized light micrographs ( Figure 5). At the maximum SPAN-80 content in S10, there was no significant difference in the value of τ 1 from control and other formulations. On adding TWEEN-80 to the formulations, the value of τ 1 was significantly reduced in T1 and T3, from the control, although τ 1 in T1 and T3 appeared similar. In T5 and T10, a considerably higher value of τ 1 from control was observed. The rise in τ 1 of T10 from T5 was found to be significant. The delayed relaxation time scripts the disruption in the oleogel network if the stress is maintained for a longer time [47]. On adding SPAN-80 to oleogels, the value of τ 2 was significantly reduced from the control in S1. Further addition in the emulsifier content in S3 displayed a similar τ 2 value to that of the control and S1. A significant reduction in the τ 2 value from the control and S3 was observed in S5. At the highest emulsifier content, the τ 2 value of S10 appeared similar to the control and other SPAN-80 formulations. However, among TWEEN-80, all the formulations showed a significant reduction from control in the τ 2 values. The longer relaxation time on the addition of TWEEN-80 to oleogels indicates the lesser likelihood of oleogels undergoing network architecture breakage [45]. In other words, the gelator network formed by the TWEEN-80 oleogels can maintain their structural integrity for an extended time under stressful conditions.  Figure 7 represents the stages of the spreadability cycle. These stages help in the visualization of the semi-solid nature of the formulated oleogels. A trigger force of 5 g was set for the analysis. Afterward, the male cone penetrated oleogels at a rate of 0.5 mm/s to a depth of 13 mm. At 26 s, after the penetration of the male cone, a clear annulus is visible (green arrow), whose diameter slowly decreases with increasing depth. The male cone compresses and extrudes the oleogel from the female cone (Figure 7, 43 s), spreading evenly between the male and female cones. Figure S1 represents the spreadability profiles of the formulations. The peak value in the positive curve represents the firmness (F 0 ), and the area under the positive curve depicts the amount of energy (C 0 ) required to deform the oleogels. The probe, while returning, incorporates the information of cohesiveness and resistance of the oleogels to separate from the cone. The peak value in the negative curve denotes adhesive force/stickiness (S 0 ). The area under the negative peak is the measure of adhesiveness (A 0 ). Unlike the SR study, we did not observe any differences in the spreadability properties as the parameters (Table S3) discussed above appeared similar in all the formulations. This observation suggests that the bulk properties of oleogels remained the same even after the addition of the emulsifiers.

FTIR Spectroscopy
FTIR spectroscopy was used to estimate the chemical nature of raw components (i.e., SO, SFW, SPAN-80, and TWEEN-80) used to prepare the oleogels ( Figure S1). The spectra observed in SFO and SFW resemble triglycerides, which are crucial components in edible oils [48]. The characteristic spectra of SO included a band at 3007 cm −1 due to the C-H stretching vibration in =C-H (cis). Since SO is composed of linoleic and oleic acids consisting of a significant amount of unsaturated fats, they are the possible explanation behind the appearance of this peak. The bands representing the C-H stretching vibration in methylene and methyl groups are present at 2924 cm −1 and 2852 cm −1 , respectively [48]. The peaks further confirmed the presence of these two groups in spectra of SO at 1459 cm −1 and 1377 cm −1 , which can be ascribed to the C-H bending vibration in methylene and methyl groups [48]. A comparatively larger peak at 1742 cm −1 was observed due to the -C=O double bond stretching vibration from the ester groups. These groups are present in abundance in SFO. The observed bands at 1236 cm −1 , 1159 cm −1, and 1097 cm −1 represent a -C-O stretching vibration corresponding to the ester group [49]. Further, the bending vibration of trans -CH=CH-is noted at 967 cm −1 and the rocking vibration of -(CH2)n at 722 cm −1 [49]. The FTIR spectra of SFW showed asymmetric stretching vibration of methyl group 2.5. Molecular Analysis 2.5.1. FTIR Spectroscopy FTIR spectroscopy was used to estimate the chemical nature of raw components (i.e., SO, SFW, SPAN-80, and TWEEN-80) used to prepare the oleogels ( Figure S1). The spectra observed in SFO and SFW resemble triglycerides, which are crucial components in edible oils [48]. The characteristic spectra of SO included a band at 3007 cm −1 due to the C-H stretching vibration in =C-H (cis). Since SO is composed of linoleic and oleic acids consisting of a significant amount of unsaturated fats, they are the possible explanation behind the appearance of this peak. The bands representing the C-H stretching vibration in methylene and methyl groups are present at 2924 cm −1 and 2852 cm −1 , respectively [48]. The peaks further confirmed the presence of these two groups in spectra of SO at 1459 cm −1 and 1377 cm −1 , which can be ascribed to the C-H bending vibration in methylene and methyl groups [48]. A comparatively larger peak at 1742 cm −1 was observed due to the -C=O double bond stretching vibration from the ester groups. These groups are present in abundance in SFO. The observed bands at 1236 cm −1 , 1159 cm −1, and 1097 cm −1 represent a -C-O stretching vibration corresponding to the ester group [49]. Further, the bending vibration of trans -CH=CH-is noted at 967 cm −1 and the rocking vibration of -(CH 2 ) n at 722 cm −1 [49]. The FTIR spectra of SFW showed asymmetric stretching vibration of methyl group through a band positioned at 2954 cm −1 . The band located at 1375 cm −1 corresponds to the symmetrical bending vibration of CH 3 [50] . Further, a distinct peak at 2916 cm −1 and 2846 cm −1 resembles the symmetric axial deformation of CH 2 . The spectral band at 1471 cm −1 and 1463 cm −1 stands for symmetric angular stretching or bending vibration of CH 2 [50]. Two continuous bands at 728 cm −1 and 718 cm −1 correspond to the in-plane CH 2 deformation [51]. A distinct peak at 1732 cm −1 corresponds to the stretching vibration of -C=O present in ester.
Structurally, SPAN-80 and TWEEN-80 consist of hydrophilic groups such as sorbitan and polyethylene glycol, which comprise several polar groups. The core lipophilic groups in them are long hydrocarbon chains of fatty acids, fatty alcohols, and esters. The band at 2922 cm −1 , and 2852 cm −1 in SPAN-80 and 2926 cm −1 , and 2863 cm −1 in TWEEN-80 spectra corresponds to the C-H stretching vibrations. However, these two bands occur as distinct peaks in the case of SPAN-80 but merge for TWEEN-80. The possible reason can be less branching in the structure of SPAN-80 as compared to TWEEN-80 [52]. Additionally, peaks in the range 720-946 cm −1 implicate the vibration of C-H deformation in the two emulsifiers [53]. The spectral band at 1738 cm −1 of SPAN-80 and 1734 cm −1 of TWEEN-80 is for stretching vibration of C=O of ester groups present in these emulsifiers. The peak 1463 cm −1 and 1377 cm −1 along with 1457 cm −1 and 1350 cm −1 corresponds to C-H scissoring vibration in SPAN-80 and TWEEN-80, correspondingly [54]. Stretching vibration of C-O-C present in the ester group was depicted by the peak present at 1171 cm −1 in SPAN-80. However, in TWEEN-80, the C-O-C stretching vibration was identified with the band at 1095 cm −1 .
FTIR spectroscopy was used to obtain information regarding the interactions among the components in the control oleogels ( Figure S2). We have noticed that all our formulations ( Figure 8) showed an overall IR spectra pattern similar to that of the control. The recorded IR spectra showed single peaks at 722 cm −1 and 1461 cm −1 . This accounts for CH 2 rocking vibration and CH 2 and CH 3 bending, respectively. FTIR spectra are capable of providing information regarding the acyl chain packaging in the lipid systems. The group splitting of CH 2 rocking (722 cm −1 ) and CH 2 bending (1461 cm −1 ) mentioned before is common to orthorhombic subcell [55]. This splitting is specific for orthorhombic packing and is absent in subcell hexagonal and triclinic packing [56]. Another set of significant peaks occurred at 2852 cm −1 and 2922 cm −1 , resulting from symmetric and antisymmetric CH 2 stretching [57]. A distinct peak at 1163 cm −1 represents symmetric-stretching. The prominent peak at 1742 cm −1 corresponds to the C=O aliphatic ester groups of triglycerides present in SFO and SFW [58]. A minor peak at 3010 cm −1 was recorded in all the oleogel formulations, also present in SFO. This peak is contributed to the bond vibration of the alkene (=C-H) from SFO. The presence of this band in the IR spectra indicates a high degree of unsaturation. There was no marked fluctuation in the peaks mentioned above in the IR spectra of SPAN-80 and TWEEN-80 formulations, possibly due to the meagre amounts of emulsifiers. The absence of any band shift was confirmed by taking the instrument's spectral resolution (4 cm −1 ) into consideration. Since no change in the peak position was observed, the obtained FTIR results confirm no significant difference in the chemical interaction in the emulsifier-containing oleogels compared to the pristine oleogel. The gel formation is solely based on non-covalent interactions such as hydrogen bonding and van der Waals attraction.  (Table 5).

XRD Diffraction Patterns
The polymorphism of the oleogels was studied through the XRD diffractograms. The diffraction patterns of oleogels prepared by adding a varied amount of emulsifiers are shown in Figure 9. The lateral packing of the fatty acid chains in the 3-dimensional networks is obtained by the short d-spacings, which can be calculated from the diffractograms at wide angles. The lateral packing of triacylglycerol has been reported through three classic organizations. These organizations include the least stable α form, metastable β' form, and most stable β polymorphs [59]. These polymorphs exist in subcell packing as hexagonal, orthorhombic perpendicular, and triclinic parallel [60]. Further, palmitic acid in sunflower oil promotes the β' polymorph in the oleogels [61]. As previously mentioned, wax comprises various molecules, i.e., hydrocarbons, esters, fatty alcohols, and fatty acids. The diffraction pattern of control showed a broad peak at 22.6° 2θ, corresponding to Bragg's distance (dspacing; interplanar spacing) of 4.56 Å (Equation (1)). A sharp peak was observed at 25.06° 2θ with a d-spacing of 4.12 Å. The third sharp peak, positioned at 27.78° 2θ, has an intensity lower than the previous one and a d-spacing of 3.72 Å. The three mentioned peaks are called "short spacing peaks," associated with hydrocarbon chains' lateral packing. The d-spacing value of 4.12 Å and 3.72 Å is an indication of the presence of β' polymorph of triacylglycerols [62]. These d-spacing values correspond to the orthorhombic subcell packing [36]. The β' polymorph of fats in food products such as margarine and shortenings is responsible for the even texture, spreadability, and mouthfeel.
Incorporating the emulsifiers increased the intensity of these peaks, which gives an idea that the emulsifiers, at the content used, have assisted in the lateral packing of chain in the fat network. Inclusion of SPAN-80 has displayed these peaks roughly at 22.94° 2θ, 25.145° 2θ, and 27.92° 2θ. The interplanar spacing displayed the following values respective to the mentioned peaks, i.e., 4.49 Å, 4.10 Å, and 3.707 Å. Similarly, in TWEEN-80 oleogels, the broad peaks occur at ~22.83° 2θ, having an interplanar spacing of 4.51 Å. The other two peaks were positioned at 25.12° 2θ and 27.92° 2θ with the corresponding interplanar spacings of 4.11 Å and 3.71 Å. These peak positions are related to the presence of β' polymorphs of wax crystals whose stability is somewhere between polymorphs α and  (Table 5).

XRD Diffraction Patterns
The polymorphism of the oleogels was studied through the XRD diffractograms. The diffraction patterns of oleogels prepared by adding a varied amount of emulsifiers are shown in Figure 9. The lateral packing of the fatty acid chains in the 3-dimensional networks is obtained by the short d-spacings, which can be calculated from the diffractograms at wide angles. The lateral packing of triacylglycerol has been reported through three classic organizations. These organizations include the least stable α form, metastable β form, and most stable β polymorphs [59]. These polymorphs exist in subcell packing as hexagonal, orthorhombic perpendicular, and triclinic parallel [60]. Further, palmitic acid in sunflower oil promotes the β polymorph in the oleogels [61]. As previously mentioned, wax comprises various molecules, i.e., hydrocarbons, esters, fatty alcohols, and fatty acids. The diffraction pattern of control showed a broad peak at 22.6 • 2θ, corresponding to Bragg's distance (d-spacing; interplanar spacing) of 4.56 Å (Equation (1)). A sharp peak was observed at 25.06 • 2θ with a d-spacing of 4.12 Å. The third sharp peak, positioned at 27.78 • 2θ, has an intensity lower than the previous one and a d-spacing of 3.72 Å. The three mentioned peaks are called "short spacing peaks," associated with hydrocarbon chains' lateral packing. The d-spacing value of 4.12 Å and 3.72 Å is an indication of the presence of β polymorph of triacylglycerols [62]. These d-spacing values correspond to the orthorhombic subcell packing [36]. The β polymorph of fats in food products such as margarine and shortenings is responsible for the even texture, spreadability, and mouthfeel.
Incorporating the emulsifiers increased the intensity of these peaks, which gives an idea that the emulsifiers, at the content used, have assisted in the lateral packing of chain in the fat network. Inclusion of SPAN-80 has displayed these peaks roughly at 22.94 • 2θ, 25.145 • 2θ, and 27.92 • 2θ. The interplanar spacing displayed the following values respective to the mentioned peaks, i.e., 4.49 Å, 4.10 Å, and 3.707 Å. Similarly, in TWEEN-80 oleogels, the broad peaks occur at~22.83 • 2θ, having an interplanar spacing of 4.51 Å. The other two peaks were positioned at 25.12 • 2θ and 27.92 • 2θ with the corresponding interplanar spacings of 4.11 Å and 3.71 Å. These peak positions are related to the presence of β polymorphs of wax crystals whose stability is somewhere between polymorphs α and β. Additionally, the β form of wax crystal displays needle and fine grain-like microstructures, which was evident in micrographs from our study [63]. A better understanding of d-spacing and their trend with the usage of emulsifiers is discussed in the next section. The observed shifts in the peak position of oleogels added with emulsifier compared to control can be due to the defect caused in the fat network. β. Additionally, the β' form of wax crystal displays needle and fine grain-like microstructures, which was evident in micrographs from our study [63]. A better understanding of d-spacing and their trend with the usage of emulsifiers is discussed in the next section. The observed shifts in the peak position of oleogels added with emulsifier compared to control can be due to the defect caused in the fat network.  Table 5).
The profiles were deconvoluted in Origin Pro software using the Gauss peak fitting function to understand the crystal properties better. The obtained deconvoluted data was capable of displaying five characteristic peaks in all the formulations. Further, this data helped to calculate the crystallite size (D), lattice strain, and dislocation density ( Table 2). On adding SPAN-80 in oleogels, the average d-spacing values showed an increment in S3, and the rest of the formulations showed values similar to control. However, the addition of TWEEN-80 did not show any significant change in the average d-spacing values. The addition of SPAN-80 and TWEEN-80 to the oleogels has improvised the crystallite size compared to control. Additionally, the large crystal size is responsible for high crystallinity, allowing wax molecules to form stable polymorphs [64]. This suggests the possibility of liquid emulsifiers used in our study to act as crystal modifiers. A deep analysis of the parameters clear that although the S3 micrographs have displayed more amorphous regions, the crystallite size in them is most significant compared to the control and other SPAN-80 formulations. The large crystallite size may be due to the slow crystallization rate of the fat molecules. The inclusion of 3 mg of SPAN-80 in the oleogel supported the fat crystal growth. Similarly, the reduced value of lattice strain and dislocation density backs fewer crystal defects in S3, which promoted the formation of the larger fat crystals. Among the TWEEN-80 formulation, T3 has shown the largest crystallite size and lower lattice strain and dislocation density. Previous work has demonstrated the potential of  Table 5).
The profiles were deconvoluted in Origin Pro software using the Gauss peak fitting function to understand the crystal properties better. The obtained deconvoluted data was capable of displaying five characteristic peaks in all the formulations. Further, this data helped to calculate the crystallite size (D), lattice strain, and dislocation density ( Table 2). On adding SPAN-80 in oleogels, the average d-spacing values showed an increment in S3, and the rest of the formulations showed values similar to control. However, the addition of TWEEN-80 did not show any significant change in the average d-spacing values. The addition of SPAN-80 and TWEEN-80 to the oleogels has improvised the crystallite size compared to control. Additionally, the large crystal size is responsible for high crystallinity, allowing wax molecules to form stable polymorphs [64]. This suggests the possibility of liquid emulsifiers used in our study to act as crystal modifiers. A deep analysis of the parameters clear that although the S3 micrographs have displayed more amorphous regions, the crystallite size in them is most significant compared to the control and other SPAN-80 formulations. The large crystallite size may be due to the slow crystallization rate of the fat molecules. The inclusion of 3 mg of SPAN-80 in the oleogel supported the fat crystal growth. Similarly, the reduced value of lattice strain and dislocation density backs fewer crystal defects in S3, which promoted the formation of the larger fat crystals. Among the TWEEN-80 formulation, T3 has shown the largest crystallite size and lower lattice strain and dislocation density. Previous work has demonstrated the potential of polysorbates to co-crystalize with fats, thus improving crystal growth [27]. In a nutshell, it can be inferred that the inclusion of both the liquid emulsifiers has introduced fluctuations in the structural architecture of the wax network. However, formulations S3 and T3 have been found to support fat crystal growth. Gelation/crystallization is a process of arranging the triacylglycerols into a compact structure. This arrangement is possible due to the physical and chemical bonds between the triacylglycerol, thus restricting its movement [65]. The gelation behaviour of fats and its understanding is essential to safeguard required industrial applications of fat-based products. The conventional mechanism of lipid structuring consists of three stages, i.e., nucleation, crystal growth, and maturation. These three stages assist in the formation of the fat crystal lattice. The graph of gelation kinetics ( Figure 10) depicts the three stages of fat crystallization. The stages include initial, intermediate, and final saturation phases [57]. The initial step is marked by a sharp decline in the temperature of formulations and is marked in blue arrows in the representative temperature versus time graphs. This can be correlated to the nucleation phase of lipid structuring and is a sign of secondary crystallization. The secondary crystallization of fats in the control is marked at a time point of 621 s, and this point in the graph differs among all the formulations. An essential aspect of our study is understanding the impact of the chosen liquid emulsifiers on fat crystallization. The effect of emulsifiers on the gelation kinetics can be due to the different organization of crystals in the arrangement or by forming imperfections. Emulsifiers with dissimilar hydrophobic properties cause substantial effects on the crystallization of fats [66]. This is usually done by altering the kinetics of crystallization or by varying transitions in the polymorphs. The addition of SPAN-80 in the oleogel lengthened the onset of secondary crystallization in the case of S1, S3, and S5 upon increasing the emulsifier's amount (Table 3). However, S10 showed a rapid onset (908 s) of secondary crystallization in the oleogels. This agrees with the previously reported data where a higher proportion of lipophilic emulsifiers have displayed reduced onset of crystallization induction time [67]. In TWEEN-80 formulations, the onset of secondary crystallization was slow in T1, T3, and T10 and was comparatively rapid in the case of T5. lization, it took less time than other formulations to reach thermal equilibrium (~2092 s). A similar trend in the initiation stage and thermal equilibrium was observed in TWEEN-80 formulations. Most of the TWEEN-80 oleogels showed delayed onset of secondary crystallization and attainment of thermal equilibrium between the two formulations. The bulky hydrophilic head and kinked carbon chain of TWEEN-80 are possible reasons that can affect the process of fat crystallization [27].  Table 5).
The crystallization curve in the early portion (0-200 s) was then analysed in-depth. The initial portion of the curve was then fitted to an exponential decay function, which is represented as; where 'a' is the initial temperature (°C), 'k' is the crystallization rate, and 't' is time (s). As a general observation, the inclusion of both the emulsifiers reduced the crystallization rate compared to the control. This suggests the possibility of alteration in the kinetics of gelation. This observation can be reasoned to the presence of sorbitan esters. The sorbitan esters not only slow down or inhibit the polymorphic transitions and stabilize the β' polymorphs of triacylglycerols [15]. A low crystallization rate of S3 further supports  (Table 5). Beyond the blue arrow, a continuous transformation is visualized in the crystallization process of all the formulations. Finally, the green arrow scripts the accomplishment of thermal equilibrium in the crystallization kinetics curve, which, like the onset point, differs among the formulations. In our study, control and S1 have reached equilibrium at the earliest~1900 s. Further, an increase in the amount of SPAN-80 delayed the attainment time of thermal equilibrium till S5. Since S10 showed a rapid initiation of secondary crystallization, it took less time than other formulations to reach thermal equilibrium (~2092 s). A similar trend in the initiation stage and thermal equilibrium was observed in TWEEN-80 formulations. Most of the TWEEN-80 oleogels showed delayed onset of secondary crystallization and attainment of thermal equilibrium between the two formulations. The bulky hydrophilic head and kinked carbon chain of TWEEN-80 are possible reasons that can affect the process of fat crystallization [27].
The crystallization curve in the early portion (0-200 s) was then analysed in-depth. The initial portion of the curve was then fitted to an exponential decay function, which is represented as; y = ae −kt (2) where 'a' is the initial temperature ( • C), 'k' is the crystallization rate, and 't' is time (s). As a general observation, the inclusion of both the emulsifiers reduced the crystallization rate compared to the control. This suggests the possibility of alteration in the kinetics of gelation. This observation can be reasoned to the presence of sorbitan esters. The sorbitan esters not only slow down or inhibit the polymorphic transitions and stabilize the β polymorphs of triacylglycerols [15]. A low crystallization rate of S3 further supports the fewer crystal defects, as confirmed from XRD. Similarly, an even lower value of rate in S5 can be the reason for the delay (~2300 s) in thermal equilibrium. Again, the addition of TWEEN-80 showed reduced crystallization rates to a greater extent in all the formulations, notably in T3 and T5. The possible reason can be the observed larger crystallite size in these formulations. In addition, the giant hydrophilic head and the unsaturated carbon chain of TWEEN-80 can cause a hindrance and thus delay the crystallization rate [27]. A research group has also attributed the delayed crystallization to the interference of emulsifiers in the packing of crystals [27]. The reduction in the rate values of crystallization is a combined effect of nucleation rate, crystal type, morphology, size, crystal arrangement, and growth [68].

DSC Analysis
Plant waxes are of natural origin, which attributes to their compositional differences and physical behaviour. The melting and crystallization behaviour is one of the physical parameters that were observed for the prepared oleogels. The colligative properties of wax help them by solubilizing into the vegetable oil. This results in the decrement of melting and crystallization temperature of wax oleogels. These temperatures of wax-based oleogels are usually lower in comparison to neat waxes [36]. Any change in the physical state of oleogel or polymorphic transitions is associated with the difference in the absorption (endothermic) or release (exothermic) of heat. This heat flow can be measured through DSC as a function of time and temperature. Figure 11 represents the DSC heating and cooling profiles of all the formulations. The obtained heating profiles of all the formulations displayed a broad peak and one shoulder peak in the temperature range 44-70 • C. The melting onset in all our formulations is from~40 • C. Deconvolution of the broad peak using Origin Pro software was performed to reveal the positions and area of the peaks. It was found out that the broad peak comprised of two endothermic signals (Table S3) in the temperature range 60-70 • C. When present in different polymorphic forms, fats usually appear as more than one endotherm [69]. The control sample's melting temperature/peak temperature (T m ) was observed through a significant peak at~62 • C along with a shoulder peak at 65 • C. Previous work reported a T m of 58.4 • C for oleogel containing 3% SFW along with the possibility of increased values at higher wax concentration [70]. The peak melting temperature of SFW and olive oil oleogels has been reported to be between 58 and 63 • C in literature [71]. The range of observed melting points in our study corresponds with that of triglycerides melting behaviour [72]. β form, the most stable form of fats, has the highest melting point [73]. SFW is a high melting wax due to the presence of highly saturated or long-chain fatty acids. Since our XRD diffractograms confirmed the presence of β polymorph, there is a possibility of melt-mediated transformation in the formulations. This transformation suggests the appearance of a stable (β) fat polymorph when a metastable (β ) form is heated over its melting point [73]. However, incorporation of a emulsifier at a certain amount slightly changed the peak positions. On adding SPAN-80 to oleogels, the melting temperature reduced compared to the control in all t formulations except for S3, where T m was observed at 64 • C. This is interesting because a high melting temperature suggests better formulation stability and correlates with slower gelation kinetics and fewer crystal defects in the S3 discussed above. Further, the conversion of β to β is associated with the formation of a crystal with a large size [74]. The presence of a larger crystallite size displayed through XRD analysis in S3 further supports its better stability among SPAN-80 formulations. However, the addition of TWEEN-80 has increased the T m value in the case of T3 and T10. The result is yet again in line with the results obtained from the XRD analysis. Crystallization onset and peak temperature have also been tabulated (Table S4). The crystallization of waxes in liquid oil starts with its supersaturation in the oil and, thus, can be affected by the amount of wax used as a gelator [75]. For every formulation, the peak in the cooling profile occurs at~60 • C irrespective of the addition of a liquid emulsifier. The appearance of a similar peak indicates the presence of similar polymorphs i.e., β form in all the formulations.
Gels 2021, 7, x FOR PEER REVIEW 20 of 29 (β') form is heated over its melting point [73]. However, incorporation of a emulsifier at a certain amount slightly changed the peak positions. On adding SPAN-80 to oleogels, the melting temperature reduced compared to the control in all t formulations except for S3, where Tm was observed at 64 °C. This is interesting because a high melting temperature suggests better formulation stability and correlates with slower gelation kinetics and fewer crystal defects in the S3 discussed above. Further, the conversion of β' to β is associated with the formation of a crystal with a large size [74]. The presence of a larger crystallite size displayed through XRD analysis in S3 further supports its better stability among SPAN-80 formulations. However, the addition of TWEEN-80 has increased the Tm value in the case of T3 and T10. The result is yet again in line with the results obtained from the XRD analysis. Crystallization onset and peak temperature have also been tabulated (Table S4). The crystallization of waxes in liquid oil starts with its supersaturation in the oil and, thus, can be affected by the amount of wax used as a gelator [75]. For every formulation, the peak in the cooling profile occurs at ~60 °C irrespective of the addition of a liquid emulsifier. The appearance of a similar peak indicates the presence of similar polymorphs i.e., β form in all the formulations.  (Table 5).

Drug Release Analysis
Curcumin is a lipophilic compound with enormous health benefits as a herbal medicine [76]. It is crucial to develop curcumin formulations that have improved bioavailability. The bioavailability at the site of action will depend on the release of curcumin from the oleogels. This release was calculated as cumulative percentage drug release (CPDR). At the end of the experiment (after 180 min), the emulsifier-containing oleogels showed a marked reduction in the %CPDR compared to the control ( Figure 12). The inclusion of hydrophobic surfactants in the oleogel might have increased the lipophilicity of the wax oleogels. Among the SPAN-80 formulations, S1 and S3 showed similar %CPDR. S5 showed a substantial increase in the %CPDR from S1 and S3. At the highest SPAN-80 content, the release was found to be identical to S5; however, it was greater than S1 and S3. The addition of TWEEN-80 in T1 and T3 showed similarity in the drug release after  (Table 5).

Drug Release Analysis
Curcumin is a lipophilic compound with enormous health benefits as a herbal medicine [76]. It is crucial to develop curcumin formulations that have improved bioavailability. The bioavailability at the site of action will depend on the release of curcumin from the oleogels. This release was calculated as cumulative percentage drug release (CPDR). At the end of the experiment (after 180 min), the emulsifier-containing oleogels showed a marked reduction in the %CPDR compared to the control ( Figure 12). The inclusion of hydrophobic surfactants in the oleogel might have increased the lipophilicity of the wax oleogels. Among the SPAN-80 formulations, S1 and S3 showed similar %CPDR. S5 showed Gels 2021, 7, 133 20 of 28 a substantial increase in the %CPDR from S1 and S3. At the highest SPAN-80 content, the release was found to be identical to S5; however, it was greater than S1 and S3. The addition of TWEEN-80 in T1 and T3 showed similarity in the drug release after 180 min. The formulation T5 led to a similar %CPDR to T3, though this value was higher than T1. At the highest TWEEN-80 content, there was a significant rise in the %CPDR from T1, T3, and T5.   (Table 5).
A further understanding of the drug release mechanism was developed using Peppas Sahlin (PS) model fitting with correlation coefficient >0.99 (Table 4), which is represented in Equation (2) [77]. The obtained correlation coefficient was an indication of the good fit of the experimental data with the modelled data. The PS model combines the effect of Fickian diffusional release and relaxation release of drugs from the polymers. The PS model was initially developed to evaluate the drug release from the hydrophilic matrix; however, attainment of good fit in this study suggested its implementation to hydrophobic polymers [78].
where F indicates the portion of the solute released; Kd and Kr are constants related to Fickian kinetics and relaxation kinetics, respectively; m is the Fickian diffusional exponent, which prevents the controlled release of the drug, and t is the time at which sampling took place. The Fickian diffusion constant (Kd) for every SPAN-80 formulation was significantly reduced in comparison to the control. The Kd value subsequently decreased till S3, however, this value showed a rise in S5 from S1 and S3. An increase of SPAN-80 content in S10 showed the least Kd among all the formulations. Slightly higher Kd values of S1 and S5 among SPAN-80 formulations are attributed to the void spaces observed in these formulations' micrographs. The voids may have aided in a comparatively faster release of curcumin. The addition of TWEEN-80 to the formulations displayed a subsequent decrease in the Kd value till T3. Formulation T3 had the lowest Kd value among all the for-  (Table 5).
A further understanding of the drug release mechanism was developed using Peppas Sahlin (PS) model fitting with correlation coefficient >0.99 (Table 4), which is represented in Equation (2) [77]. The obtained correlation coefficient was an indication of the good fit of the experimental data with the modelled data. The PS model combines the effect of Fickian diffusional release and relaxation release of drugs from the polymers. The PS model was initially developed to evaluate the drug release from the hydrophilic matrix; however, attainment of good fit in this study suggested its implementation to hydrophobic polymers [78].
where F indicates the portion of the solute released; K d and K r are constants related to Fickian kinetics and relaxation kinetics, respectively; m is the Fickian diffusional exponent, which prevents the controlled release of the drug, and t is the time at which sampling took place. The Fickian diffusion constant (K d ) for every SPAN-80 formulation was significantly reduced in comparison to the control. The K d value subsequently decreased till S3, however, this value showed a rise in S5 from S1 and S3. An increase of SPAN-80 content in S10 showed the least K d among all the formulations. Slightly higher K d values of S1 and S5 among SPAN-80 formulations are attributed to the void spaces observed in these formulations' micrographs. The voids may have aided in a comparatively faster release of curcumin. The addition of TWEEN-80 to the formulations displayed a subsequent decrease Gels 2021, 7, 133 21 of 28 in the K d value till T3. Formulation T3 had the lowest K d value among all the formulations, suggesting an equal contribution of diffusion and relaxation in the release process. Further rise in the emulsifier content in T5 displayed a significantly higher K d value than T3, but this value was found to be similar to the control. Formulation T10 showed the highest K d value, which was even higher than the control, suggesting a dominance of diffusional mechanism in this case. The higher values of K d than K r (or K d /K r > 1) in all the formulations indicated the importance of diffusional contribution towards the curcumin release from the SFW-based oleogels [79].

Conclusions
A continuous change in the global scenario of dietary patterns and guidelines has prompted scientists to discover the alternative of trans fats. An opposite substitute of trans-fats should not compromise the functional integration of the food products, which otherwise will not be accepted by the consumers. The development of oleogels in lipid research has emerged with an exciting potential to replace saturated fats. The present study showed the potential of SFW to bind the SFO, which was formulated using the heating-and-cooling method. Further, the effects of liquid hydrophilic and hydrophobic emulsifiers were studied through various techniques such as colorimetry, microscopy, texture analyser, FTIR spectroscopy, XRD, crystallization, and DSC. Irrespective of the emulsifier treatments, all the 5% SFW oleogels have shown more than 99% OBC. The formulated oleogels had an ivory appearance and smooth touch, ideal for many food applications, including shortenings and margarine [71]. Our study revealed an amount of 3 mg of both SPAN-80 and TWEEN-80 in 20 g of oleogels is optimum for structural, molecular, and thermal parameters when included in 5% SFW oleogels. The surface topography images displayed the uniform distribution of globular structures. The addition of both the emulsifiers has shown an alteration in globular size and distribution. Topographs of S3 showed crystals that were most noticeable and bright. Comparable globular structures were also observed in T3. The wax esters of SFW were able to form a fibrous network of varied length, thickness, and branching, which were clearly visible under bright field micrographs. The SPAN-80 formulations, notably S3, were more hyperbranched and dense than TWEEN-80 formulations, which specify a good mechanical strength and oil binding capacity. The examined viscoelastic properties of the oleogels through SR profiles revealed alteration in the firmness, residual elastic force, and %SR. These properties were not entirely dependent on the emulsifier's content. The formulations, S3 and T3, displayed lower spreadability and higher inherent elastic properties. All the formulations showed an overall IR spectra pattern similar to that of the control, thus suggesting a not significant difference in molecular aspects. The XRD diffractograms revealed the presence of metastable β form in all the formulations. Further, S3 and T3 showed large crystallite sizes compared to the control sample, suggesting that the emulsifiers were able to induce optimum conditions that supported the crystal growth of the fats. The crystallization kinetics study under isothermal conditions showed that both the emulsifiers had reduced the crystallization rate in all the formulations as compared to the control sample. A slower crystallization rate in S3 and T3 supported the formation of larger crystal sizes with the possibility of fewer crystal defects. An interesting case of melt-mediated transition of β to β was revealed through the heating profile of DSC thermograms. This was prominently observed in S3 and T3, further suggesting the presence of a stable polymorph in these formulations. Based on the observation mentioned above, it can be said that the minimum acceptable amount of emulsifiers in food products still have the potential to act as crystal modifiers. SFW oleogel with improvised properties in this way can act as a fat replacer in many bakeries and other food products.

Oleogel Preparations
Five percent of SFW was found to be the critical gelation concentration (CGC) for SFO. Hence, an oleogel of SFO containing 5% (w/w) of SFW was prepared, which served as the control formulation. In gist, accurately weighed SO and SFW were taken in a beaker placed in a water bath (80 • C) for~20 min. The mixture was heated until all the wax was dissolved in the oil, and a homogeneous solution was formed. This clear solution was then kept at room temperature (25 • C) in the thermal cabinet for 30 min to encourage oleogel formation. Stock solutions of 0.1% (w/w) of SPAN-80 and TWEEN-80 were prepared in SO. A series of 1 mg, 3 mg, 5 mg, and 10 mg of SPAN-80 and TWEEN-80 added formulations were prepared following the composition as per Table 5. The rest of the process for the development of the formulations was the same.  The study of oil released was performed initially to find out the OBC of formed oleogels. For this, an empty 2 mL previously weighed (a) Eppendorf was taken, and 1 mL of the molten sample was placed in the Eppendorf. The sample was then allowed to solidify, followed by weighing the Eppendorf (b). The Eppedorfs were kept at 4 • C for 24 h, after which they were centrifuged with Remi C-24 BL refrigerated centrifuge at 10,000 rpm (18 • C) for 15 min. The Eppendorfs were then inverted to drain the released oil. The remaining oil over the surface of the oleogels was removed using filter paper. The %oil released and OBC was calculated using Equations (3) and (4), respectively [78].
] × 100 (4) OBC = 100 − %Oil released (5) where, a is the 2 mL previously weighed Eppendorf; b is the solidified sample containing Eppendorf; and c is the final weight of Eppendorf after centrifugation and draining.

Colorimetric Analysis
Colorimetric analysis of the oleogels was done with a lab-developed colorimeter. The hardware of the colorimeter device consists of a light-emitting diode (LED) as a light source and Picam for imaging purposes [79]. The instrument was calibrated before the experiment with the standard white and black tiles. The analysis was done using oleogels that were placed in 35 mm Petri-dishes. Thereafter, the colour coordinates L*, a*, and b* values were measured.

Surface Topography
The surface topography of the formulations was visualized using a Stereo Zoom Microscope (Model: SM-2TZ; Make: AMscope, Irvine, CA, USA). The microscope was equipped with an external eyepiece lens camera (AMscope MD500, Irvine, CA, USA). The samples for the study were prepared by transferring 5 g of the molten mixture in 35 mm petri-dish and solidifying the same as per the method reported in Section 2.1.

Microstructure Visualization
The molten oleogels drop was placed on a glass slide followed by covering with a coverslip. The microstructure arrangement of the prepared oleogel was visualized using an upright bright-field compound microscope (Leica Microsystems, model: DM750, GmbH, Wetzlar, Germany), coupled with an in-house built polarizer. The micrographs were visualized both under bright-field and polarizing modes.

Stress Relaxation
Mechanical properties of the oleogel were examined by stress relaxation (SR) studies. The study was conducted using the texture analyser HD plus instrument (Godalming, Surrey, UK). For the SR study, the oleogels were prepared in 100 mL polypropylene beakers. The initial trigger force was set as 5 g. The solidified oleogel was penetrated with a constant strain of 1 mm at a rate of 0.5 mm/s using an acrylic male conical probe (angle 45 • ). The strained condition was maintained for 60 s, and the values for change in forces were recorded. Then, the probe was brought back to its usual height, i.e., 30 mm from the surface. Using the force values, %SR was calculated for each oleogel sample using Equation (5) mentioned below: The spreadability of the oleogels was performed texture analyser HD plus instrument (Stable Microsystems, Godalming, UK). The molten formulations were kept in a female cone and allowed to solidify at 25 • C for 90 min. The initial trigger force was set as 5 g. The 45 • male perspex cone was allowed to penetrate the female cone at 0.5 mm/s, keeping a clearance of 5 mm from the base. After that, the male cone was brought back to its usual height, i.e., 30 mm from the surface. 4.7. Molecular Characterization 4.7.1. Fourier Transform Infrared (FTIR) Spectroscopy IR spectra of all the oleogels were studied using an FTIR spectrophotometer (Alpha-E; Bruker, Billerica, MA, USA) that was working in the attenuated total reflectance (ATR) mode. The oleogels were scanned within the wavenumber range 600-4000 cm −1 , each with 25 scans. The spectral resolution of the instrument was 4 cm −1 . The IR spectra of pure SFO, SFW, SPAN-80, and TWEEN-80 were also collected similarly.

X-ray Diffraction Study
XRD of the oleogels was performed using an X-ray diffractometer (Model: Bruker D8 Advance, Austin, TX, USA) that was fitted with a Co-Kα radiation source (wavelength = 1.79 Å). The analysis was done at a voltage of 35 kV and a current of 25 mA. The X-ray scanning was carried out in the 2θ range from 5 • to 50 • and a scan rate of 5 • 2θ/min. The XRD parameters of d-spacing (d) were calculated using Brag's law (Equation (6)), and the crystallite size (D) was calculated using Debye-Scherrer equation (Equation (7)). Further lattice strain ( ) and dislocation density (δ) were calculated using Equations (8) and (9).
where λ = 1.79Å, i.e., the wavelength of the X-ray, n is an integer value, and θ is the diffraction angle.
where k is Scherrer constant, β is in radian and represents the full width at half maxima (FWHM) at a scattering angle 2θ.

Gelation Kinetics
The gelation kinetics of oleogels were studied with a lab-built temperature sensor system. For the experiment, 10 g of molten oleogels were transferred in a glass bottle (volume: 15 mL). The oleogel containing glass bottles were heated at 80 • C in a water bath. The bottles were then connected to the temperature sensor system and placed in the refrigerated water bath at 5 • C. Gelation of the formed oleogel was recorded for 90 min after the molten solution achieved a temperature of 50 • C through a time v/s temperature plot.

Differential Scanning Calorimetry (DSC) Analysis
The thermal properties of the oleogels were inspected using a differential scanning calorimeter (200 F3 DSC, Maia, Netzsch, Burlington, MA, USA) at a thermal scan of 5 • C/min for their melting and crystallization capabilities. For this study,~10 mg of oleogel was weighed in an aluminium pan and sealed using a pierced lid. A single program was run for each oleogel where the melting profiles from 0 • C to 100 • C and a cooling profile from 100 • C to 0 • C were recorded. In between each run, the oleogels were maintained at 100 • C for 5 min.