Response surface methodology for optimization of bio‐lubricant basestock synthesis from high free fatty acids castor oil

In this paper, an eco‐friendly single‐step process for the synthesis of biolubricant basestock from high free fatty acid (FFA) castor oil (CO) via epoxidation reaction was investigated. Influence of various process parameters on the structural modification of CO and their interaction with the maximum oxirane oxygen content (OOC) was optimized. Central composite design (CCD) as one of the tools in response surface methodology (RSM) was used to evaluate the effects of process variables on maximum OOC. Iodine value (IV) and OOC was used to monitor the progress of epoxidation. From the RSM study, the optimal condition inferred was H2O2, 1.65 mol; catalyst loading, 15.14 wt%; temperature, 52.81°C; and reaction time, 2.81 h. At this optimum condition, OOC was found to be 3.85 mass%. Further, the epoxide product was confirmed by 1H, 13C NMR spectral technique and OOC was determined by the standard HBr method. Finally, the significant physico‐chemical properties for the prepared epoxide were determined and compared with the castor oil.


Introduction
Currently, plant seed oil/vegetable oil-derived renewable products have been gaining much importance to replace the conventional source of energy due to their depletion at a faster rate [1,2]. In addition, the use of fossils has stimulated the search for eco-friendly alternatives to conventional resources. Nevertheless, inexhaustible renewable resources can supply a raw material basis for day-to-day life products, and this can avoid contribution to green house effects due to CO 2 discharge minimization [3]. In addition to that, proper utilization of renewable raw materials such as plant seed oils can bridge the gap between the fossil reserves' demand and consumption in the future. Above all, renewable materials could meet the principles of green chemistry in terms of easy degradation and lower toxicity [4]. Moreover, the plant oils offer a wide number of advantages discussed by various technocrats [1,2,4]. In this regard, lubricants are one of the areas, which demand an alternative to the conventional lubricant basestocks due to diverse environmental issues reported by many researchers [5,6]. Therefore, an effective utilization of the bio-based feedstocks for the lubricant basestock synthesis could bring down the dependency on imported petroleum as well as promote the sustainable agricultural initiative [6].
In general, liquid lubricants are the most common form of lubricants; their composition consists of 70-99% basestock and 30-1% of additives to improve the performance properties. However, the ultimate performance of lubricant depends on the basestock which can be synthesized by using plant seed oils [7]. Replacement of lubricant basestocks with the plant-based resources offers a wide range of advantages [5][6][7][8]. Even though the plant seed oils offer many advantages, they restrict their direct use due to the inadequate thermal and oxidative stability, hydrolytic stability, and poor cold flow properties [9]. These negative impacts are due to the presence of bisallylic protons in plant seed oils structure which are highly susceptible to free radical attack, and thereby it undergoes oxidative degradation to form polar oxy compounds [10][11][12][13]. Plant seed oil oxidizes similarly to the hydrocarbon mineral oil by following the same free radical oxidation mechanism, but the oxidation rate of plant seed oils is faster than the hydrocarbon mineral oils [14,15]. This fast rate of oxidation can be attributed to the presence of unsaturated fatty acids in its composition [16]. On the other hand, thermo-oxidative stability, hydrolytic stability, and low temperature performance of nonedible plant oils are much low and poor [17]. However, chemical modification of triglycerides can eliminate the poly unsaturation (for better thermo-oxidative stability), and an optimal extent of chemical change can ameliorate low temperature behavior [18].
Till date, many studies are available in the literature on structural modification of different plant seed oils and their methyl esters to prepare lubricant basestock by various methods, such as epoxidation (chemical or structural modification) [19], genetic modification [20], blending with additives [21], and hydrogenation of double bonds [22]. Among these methods, one of the most significant methods is hydrogenation of unsaturated double bonds in vegetable oils and its methyl esters [23,24]. Hydrogenation yields total saturation of double bonds, resulting in miserable cold flow properties [23]. Genetic modification approach is tried either to reduce the saturated fatty acid content of the plant seed oils or to reduce the polyunsaturated fatty acid content, as these constituents have a negative impact on the thermo-oxidative stability [23,25]. Genetically modified oils consistently exhibited improved oleic acid content ranging from 84 to 88% during multiple environmental conditions [26,27]. However, among all these methods, structural alteration of unsaturated bonds via epoxidation has gained much attention due to high reactivity of three-member oxirane rings. Chemical modification of plant oils at the double bond sites' results in improved thermo-oxidative stability of the modified product [28].
The present study is focused on epoxidation of castor oil (nonedible oil) which is abundantly available all over India. Castor plant has the botanical name of Ricinus Communis of the family Eurphorbiacae [29]. Castor plant is primitively a tree or shrub that can grow in most tropical and subtropical countries above 10 m high, reaching an age up to 4 years [30]. Especially, castor plant needed a temperature between 15 and 38°C with lower humidity throughout the growing season in order to receive maximum oil yields [30,31]. Castor oil is a colorless to pale yellow liquid with mild/no odor or no taste, on an average CO seed contains about 46-55% oil by weight [29]. Oil fraction of castor seeds contains a higher amount of the ricinoleic acid as a hydroxylated fatty acid, and this unique structure has given a unique identity to this biological source for industrial synthesis of the variety of compounds [32]. Around the globe, India is the world's largest exporter of castor oil, and the other major producers are China and Brazil [30]. The total world production of castor seeds is estimated around one million tons, and the oil extracted is about 500,000 tons with productivity of 470 kg of oil per hectare [30,33]. Irrespective of castor seed origin, season in which it has grown, its fatty acid composition remains unique [34].
Goud et al. [35] reported the epoxidation of castor oil with acetic acid and formic acid, using Amberlite IR-120 as heterogeneous acid catalysts [36,37]. Further, Salimon et al. [38] also discussed the synthesis of ricinoleic acid epoxide and its characterization. Likewise, Salih et al. [14] described the synthesis of biolubricant basestocks from chemically modified ricinoleic acid-based tetra-esters via epoxidation, ring opening, and esterification reactions [39]. Recently, Hajar et al. [32] reported the production of biolubricant from castor oil substrate using novozyme 435. But very scanty information is available in the literature on the optimization of high FFA castor oil epoxidation process. Therefore, in the current study, an attempt has been made to examine the behavior of the epoxidation reaction and the physicochemical properties of epoxide prepared from high FFA raw material. However, very scanty information is available in the literature on optimization of high FFA castor oil epoxidation to understand the effect of reaction variables on oxirane oxygen content. Therefore, the present communication is aimed to bridge this gap by conducting the study on the modification of high FFA castor oil structure via epoxidation. The reaction was aimed at higher OOC and the main focus was to understand the effect of process parameters and interaction among them. RSM was adopted to optimize the epoxidation parameters and CCD was applied to understand the effect of process variables, using analysis of variance (ANOVA).

Materials and analytical techniques
Castor seeds were collected from Cherukupalli (Andhra Pradesh, India). Hydrogen peroxide (purity = 50% v/v) was purchased from Rankem, ion-exchange resin (Amberlite IR-120, strong acid), glacial acetic acid (purity = 99-100%) was obtained from Merck India Ltd. All other reagents used for analysis were of analytical grade and used as received.
Extraction of CO was carried out according to the standard AOAC method. Similarly, physico-chemical characterization of all the samples was carried out as per the standard method reported earlier in some of our studies [16,40].

Thin layer chromatography
The silica coated aluminum thin layer chromatography (TLC) plates were used for the analysis. The solvent system used was a mixture of hexane and dichloromethane (CH 2 Cl 2 ) at 2:7 ratio to which few drops of glacial acetic acid were added, and the sample coated plates were visualized in an iodine chamber.

Epoxidation experimental design
Interaction between the process variables on maximum OOC was judged by RSM. A full factorial CCD technique was used for the optimization of the CO epoxidation process variables. Four reaction variables such as substrate molar ratio, that is C=C bonds to H 2 O 2 mole ratio, catalyst loading (wt%), temperature (°C), and reaction time (h) were chosen to understand the effect of these variables on CO epoxidation reaction, as these variables are highly responsible for maximum OOC. The range of epoxidaton reaction variables involved in this study are shown in Table 1 along with lower (−1), medium (0), and higher (+1) levels of the variables. Value of α (alpha) is fixed at level 2 (α = 2 4/4 ). A 2 4 full factorial CCD for four independent variables was used by giving a total number of 30 (=2 n + 2n + 6) experiments, where 'n' is the number of independent variables. During the optimization study, eight axial experimentations and 16 factorial runs were carried out with six extra replications at the center of design to estimate the pure error.

Epoxidation reaction procedure
Epoxidation of CO was carried out in 250 mL three necked glass reactor equipped with five blade glass stirrer and condenser; the entire setup was immersed in a heating oil bath. During the experiment CO, hydrogen peroxide and acetic acid were measured in a molar ratio; ion-exchange resin (Amberlite, IR-120) was added in weight% based on the organic phase. Initially, 20 g of CO (0.07 mol) was transferred into the reactor and heated to the desired reaction temperature (60°C). Then glacial acetic acid 2.1 g (0.035 mol) and other reactants, that is, hydrogen peroxide and catalysts were added to the reaction mixture. Processing time and temperature conditions used are mentioned in an experimental design matrix (Table S1). Addition of hydrogen peroxide was carried out drop wise for the first-half hours, when the temperature was 5°C below the reaction temperature to avoid the explosion. During epoxidation 14 g stirring speed was maintained to ensure consistent mixing. After complete addition of hydrogen peroxide, the reaction was continued for the desired time duration as mentioned in Table S1. Upon completion of the reaction and prior to the analysis the samples were washed repeatedly with warm Millipore water (40°C) to make it neutral. The sample was concentrated by rotary evaporator. The OOC was determined for the final epoxide product after each run. All the measurements were carried out in duplicate, and the average values are reported.

Preliminary study of process variables
Single parameter optimization process was followed to decide the optimum range of variables. Initially, the reaction time was optimized by continuing the reaction up to 9 h and samples were drawn at regular 1 h intervals to estimate the OOC. Linear increases in the epoxide content were observed up to the reaction time of 3 h, beyond which a gradual decrease was noticed in OOC ( Fig. 1). A similar trend was noticed for all other reaction variables. Hence, depending on the preliminary studies, ranges of all the process variables were chosen and coded using the reported expression [41].

Statistical analysis
Based on experimental data shown in Table S1, the regression coefficient was determined by design expert software 8.0.7.1 trial version to predict the process response as a function of independent variables and their interactions were used to understand the system behavior. The mathematical relationship between the process variables and response was calculated by the following quadratic polynomial expression: where Y is the response, that is, the epoxide content, X i and X j represent the independent variables, β 0 is constant, β i is linear term coefficient, β ii is the quadratic term coefficient, β ij is cross-term coefficient and 'n' is the number of process variables studied and optimized during the study. ANOVA was carried out to estimate the effects of process variables and their possible interaction effects on the maximum OOC in the response surface regression procedure. The goodness and best fit of the model was evaluated by a regression coefficient R 2 . The response surface and counter plots are obtained using the fitted quadratic polynomial equation generated from regression analysis by keeping two of the independent variables at central value (0) and varying the other two.

Results and Discussion
Model fitting and ANOVA analysis In order to optimize the epoxidation process variables for maximum epoxide content, a three level, four-factorial CCD was favored. ANOVA as a multivariate technique was studied to determine optimum reaction conditions. All the 30 designed experimental runs (Table S1) were performed and the results were analyzed by multiple regression analysis ( Table 2). A quadratic polynomial equation was obtained from the experimental data to predict the epoxide content as shown below in terms of coded variables. Figure 1. Preliminary study data for castor oil epoxidation to find the effects of individual reaction variables (time, catalyst loading, substrate ratio and temperature) on response.
In order to ensure a thorough model fit, to measure the analysis of variance on individual model coefficients test for lack-of-fit need to be estimated. The lack-of-fit is an assessment of failure of a model to represent the data that cannot be reported by random error [42]. Generally, the significant process variables are decided based on the F-value or P-value (also known as a probability of error value or "prob > F" value) [42]. Greater the magnitude of F-value and correspondingly smaller the "prob > F" value, more important is the corresponding coefficient [42]. The results of second-order response surface model in the form of ANOVA for the maximum epoxide content are summarized in Table 2. From the table, it can be seen that F-value of the model is 2261.27, and the corresponding P-value (prob > F) is very small, that is <0.0001 implying that the model is highly significant. The P-values are adopted as a tool to ensure the importance of each co-efficient. In this study, the main linear effects of time (A), temperature (B), catalyst loading (C), and hydrogen peroxide molar ratio (D), cross and quadratic effects of all four process variables (i.e., AB, AC, BC, BD, A 2 , B 2 , C 2 , and D 2 ) are highly significant, as the P-value is very less (<0.0001). The other model term (cross variable) time and hydrogen peroxide molar ratio (AD) are also significant variables, since the P-value is 0.0309 (<0.05). In the present model, the absence of insignificant parameters intending that all the liner, cross, and quadratic terms are highly important for maximum epoxide content. The insignificant lack-of-fit F-value of 4.26 indicates that lack-of-fit is considerably significant relative to the pure error, which signifies that the model is extremely accurate without any noise, and the results are reproducible. Furthermore, the actual values are very much close to the predicted as shown in Figure 2.
The precision of a model is judged by the regression coefficient (R 2 ). The R 2 value is always in between 0 and 1, and its order of magnitude suggests the aptness of the model [8]. For a good statistical model, the R 2 value should be close to one and the regression value for higher epoxide content is presented as 0.9996, which is close to 1, and it signifies that the 99.96% model behavior can be interpreted for higher epoxide content while only about 0.04% of the full variance cannot be explained by the model. Regression coefficient R 2 represents that the accuracy and general ability of the polynomial model is good. The predicted R 2 value (0.9982) is in reasonable agreement with the adjusted R 2 (0.9993) which recommends prominent corelational statistics between the remarked values and the predicted data. Thus, the regression model provides an excellent explanation of the relationship between the independent process variables and the response variable [42].

Influence of various process variables on maximum OOC
In order to estimate the best reaction condition for maximum epoxide content, the effects of linear, cross and quadratic reaction variable on the epoxidation was studied. The three-dimensional response surface plots and the two-dimensional counter plots, which is the graphical histrionics of the regression equation, are obtained by employing design expert software [26,43]. These graphical representations of the plots are shown by varying two process variables at a time while keeping the other two variables at a central level (0). During each experiment, samples were withdrawn at regular intervals, washed, neutralized, and analyzed for OOC and α-Glycol content (i.e., oxirane cleavage).
Effect of time and temperature on OOC Figure S1a and b describes the two-dimensional counter plot and corresponding three-dimensional response surface plot expressing the cross effect of time (A) and temperature (B) on the maximum OOC with respect to other process variables. Figure S1a and b shows that OOC increased with an increase in reaction time from 2 to 4 h, the increasing trend was noticed up to a certain time period of 3 h, beyond which gradual depletion in OOC was noticed (Table S1). Similarly, temperature was varied from 40 to 80°C in order to study the effect of temperature on OOC with respect to time. Results in Figure S1a and b show that OOC content increased with an increase in temperature up to 60°C, further increase in temperature, say 70°C OOC decreases (Table S1). The behavior of 2D and 3D figures revealed that longer reaction time and elevated temperature (beyond 60°C) lowers the OOC which might be due to oxirane cleavage (Fig. S2) [44]. The maximum OOC of 3.85 mass% was attained at the medium temperatures of 60°C and 3 h reaction time. Similar behavior was noticed by Salimon et al. [38] during their study on epoxidation of ricinoleic acid. In the current work, less OOC was obtained, which may be due to higher acid value of castor oil. Okieimen et al. [45] observed similar results during their study on epoxidation of high acid value rubber seed oil. Hence, 60°C temperature was considered as suitable for epoxidation of CO, when molar ratio of unsaturation: acetic acid: hydrogen peroxide was 1:0.5:1.5 and catalysts loading 15 wt%.

Effect of time and catalyst loading on OOC
Keeping time (A) and catalyst loading (C) at central levels, that is, 3 h and 15 wt%, respectively, combined the effect of temperature (B) and hydrogen peroxide (D) was studied aiming at the higher epoxide content. From Figure S3a and b, it can be seen that OOC content increased almost linearly with a certain time period along with catalyst loading, thereafter OOC decreases (Table S1). The obtained results attributed to the hypothesis that availability of enough active surface area of catalyst for longer reaction time during the reaction results in cleavage of oxirane oxygen [44]. Hwang and Erhan [46] reported that presence of higher FFA content led to hydrolysis reaction in acid media thereby decreases the epoxide content due to oxirane cleavage. Hence, within the experimental conditions used in this study, the most favorable catalyst loading and reaction duration appeared to be 15 wt% and 3 h.
Effect of time and substrate molar ratio on OOC Figure S4a and b shows the effect of varying substrate molar ratio and reaction time on the maximum epoxide content. Increase in the substrate molar ratio from 0.5 to 1.5 mol (Table S1) leads to linear increased in OOC content due to the formation of more peracetic acid and thereby OOC. The maximum OOC was attained with the mole ratio, 1.5 mol and the stability of oxirane observed with this ratio was similar to that observed with 1 molar ratio. However, an increase in the final oxirane value was relatively less when the ratio was further increased to 2.5 mol [47]. The use of a higher substrate molar ratio raised additional problem of agitation and decreases the mass transfer rate thereby decreases the OOC. The maximum substrate ratio and higher reaction time provided an opportunity to react oxirane rings with excess hydrogen peroxide, acetic acid, and by-product water [48]. Therefore, for subsequent experimentation, optimum reaction time, H 2 O 2 moles considered are 3 h and 1.5 mol.

Effect of temperature and catalyst loading on OOC
Effect of reaction temperature (B) and catalyst loading (C) on the maximum OOC was investigated by varying catalyst loading from 5 to 25 wt% (Table S1) and temperature from 40 to 80°C with an interval of 10°C. Increase in the temperature ( Fig. S5a and b) showed a favorable effect on epoxidation reaction and maximum OOC was obtained (3.85 mass%) at 60°C and catalyst loading 15 wt%. From Figure S5a and b, it can also be seen that beyond 60°C decrease in OOC was observed. Epoxidation at higher temperatures (>60°C) acts as a medium to oxirane cleavage thereby decreasing the OOC value [49]. Likewise, altering the catalyst loading from 5 to 25% roaring in the epoxide content was observed up to 15 wt%. Beyond which inadequate OOC was noticed due to excess catalyst loading, which leads to oxirane cleavage [39]. However, under a given reaction condition and 15 wt% catalyst loading, it is assumed that reaction is free from mass transfer resistance and the maximum OOC can be obtained at a moderate reaction temperature of 60°C [39]. Hence, 15 wt% catalyst loading and 60°C temperature was considered as the optimum parameters for subsequent experimentation.

Effect of temperature and substrate molar ratio on OOC
The effect of the substrate molar ratio (D) on reaction temperature (B) and their combined interaction during the epoxidation process at constant catalyst loading of 15 wt% and 3 h of reaction time is represented by 3D response surface plot and corresponding 2D contour plot in Figure S6a and b. From the figures, it can be depicted with an increase in both temperature and substrate molar ratio, there was a progressive increase in the OOC which indicates that both the variables have a significant interaction between each other. However, the decrease in OOC was noticed beyond 60°C [39] and substrate molar ratio 1.5. Derawi et al. [50] reported that the rate of epoxidation increased as the concentration of hydrogen peroxide increased in the system, but the stability of oxirane rings was poor at higher substrate molar ratios [51]. On the other hand, Dinda et al. [39] have shown that oxirane ring was quite stable at lower hydrogen peroxide concentrations. Finally from the interaction of these variables, it can be seen that at a medium reaction temperature, the OOC increases with an increase in substrate molar ratio. However, higher reaction temperature and molar ratio showed a negative impact on OOC. Hence, a temperature 60°C and 1.5 substrate molar ratio was considered as optimum condition for further experimentations.

Effect of catalyst loading and substrate molar ratio on OOC
The influence of catalyst loading on varying hydrogen peroxide molar ratio was investigated at various catalyst loadings and various H 2 O 2 molar ratios. Enhanced rate of per acid formation was anticipated with an increase catalyst loading and hydrogen peroxide concentration followed by an increase in the OOC [49,52]. Figure S7a and b describes the effect of catalyst loading and hydrogen peroxide on the course of epoxide formation. Higher selectivity was achieved at 1.5 moles of hydrogen peroxide per mol of ethylenic unsaturation and catalyst loading 15 wt%. Concentrations of hydrogen peroxide higher than 1.5 mol leads to a higher rate of epoxy ring decomposition as resolved by most of the researchers [53] and also discussed in the previous section. This observation was in agreement with the results shown in Figure S7a and b. Above all, lower epoxide content for high FFA CO was due to the hydrolysis reaction in the aqueous phase and its mechanism is well explained by Blee et al. [54].

Process optimization for epoxidation of high FFA CO
In the current study, this model is aspired to find the epoxidation process for the best fit of variables that gives maximum OOC. The predicted values of epoxidation reaction obtained from the model equation for a maximum amount of epoxide content are substrate ratio 1.65 mol, catalyst loading 15.14 wt%, and 2.81 h of reaction time at 52.8°C reaction temperature (Table S2). These values were chosen from the optimum solutions proposed by RSM optimization tool. The model expects the maximum OOC that can be obtained at these optimum conditions is 4.09 mass%. To confirm the model prediction, the best response variables are tested at the theoretical condition by conducting the laboratory experiment. At the optimized process condition epoxide content was found to be 3.85 mass% which agrees well with the model predicted value suggests that the formulated model was believed to be accurate and reliable (Table S2). In the present study, all the 30 experiments (Table S1) including confirmatory experiments (Table S2) were performed in duplicate at laboratory conditions and the average values are reported.
Physico-chemical characterization of prepared high FFA CO epoxide product

NMR spectroscopy
Epoxidation of high FFA CO was confirmed by 1 H NMR spectroscopy; it is one of the significant techniques to confirm the product formation and monitor the progress of reaction. 1 H NMR spectrum of the starting material (CO) and final product (CO epoxide) are shown in Figure  S8a and b. Significant signals in the 1 H NMR spectra of CO shows the presence of unsaturation bands (-CH=CH-) from 5.25 to 5.58 ppm range (Fig. S8a). However, the partial disappearance of unsaturation bands have been noticed (Fig. S8b) after completion of the reaction, which signifies the presence of unsaturation in the end product at optimum condition. Similarly, the peaks (-CH-O-CH-) at 2.8 to 3.1 ppm range indicate the formation of the epoxide product, which are absent in the oil spectra. Furthermore, characterization of the epoxide product by 1 H NMR revealed the presence of the ring open product, which was confirmed by the additional peak at 3.41 ppm. Although 1 H NMR confirmed the complete conversion of unsaturation, but minimal oxirane cleavage was noticed due to the high FFA content of oil. Padmasiri et al. [47]; Derawi and Salimon [50]; Blee et al. [54] also noticed similar results during their study on mee and rubber oil. Similarly, Figure 3A and B show 13 C NMR spectra of CO and its epoxide. Appearance of peaks at 125-133 ppm represents the olefinic carbons in CO (Fig. 3A). Comparison of Figure 3A and B, revealed the formation of epoxy protons at 53-59 ppm (Fig. 3B) and disappearance of olefinic carbons in the epoxide product at 125-133 ppm (Fig. 3A). The appearance of additional peaks at 53 and 59 ppm in the CO confirms the epoxidation of high FFA CO. Similar results are noticed by Madankar et al. [55] during their study on canola oil.

TLC analysis
Thin layer chromatography analysis of CO and its epoxide was carried out in order to confirm the product and mixture of compounds present in the final product (Fig. 4).
The TLC spectra of the epoxide product indicate the presence of polar compounds, which was found to be absent in the CO confirms the formation of epoxide. From the CO epoxide TLC spectra, it was clear that some starting material is still present in the end product, which indicates an incomplete reaction. The same can be noticed from IV of the epoxide product (Table 3) and 1 H NMR spectrum ( Fig. S8a and b) [47].

Pour point determination by DSC
Estimation of pour point is important in order to determine the low-temperature flow behavior of the modified (i.e. epoxide) and unmodified CO. It is a rough indication of the lowest temperature at which the epoxide is promptly pumpable [14]. In the present communication, cold flow properties of CO and its epoxide are determined using the method followed by Borugadda et al. [16]. DSC thermograms of CO and its epoxide are shown in Figure 5A and B. From the thermograms pour point of CO and CO epoxide was found to be −20 and −15 o C, respectively, which revealed that after epoxidation (i.e., the structural modification of CO), the pour point reduced significantly (5°C). This behavior is attributed to the fact that conversion of unsaturation content into epoxide  altered the cold flow properties (PP). However, Soriano et al. [56] also observed similar kind of results during PP analysis of oils and esters. Elaborate literature survey on vegetable oils cold flow properties revealed that the cold flow property of plant oils is extremely inadequate, and this restricts their use at lower operating temperatures [10]. Plant oils have a tendency to form macro-crystalline structures at low temperatures through uniform stacking of the triglyceride backbone [10]. Formation of the macrocrystals restricts a free flow of the fluid due to loss of kinetic energy of individual molecules during the selfstacking [10]. In general, the pour point should be low enough to ensure that the epoxide is pumpable at lower temperatures in an application point of view [8]. Finally from this study, it can be seen that the modification in the structure of CO can alter the PP which was regarded as one of the major concerns for the plant oils to be used as an alternative to fossil resources. However, further improvement in the low-temperature properties can be done by additivation.

Thermo-oxidative stability
Thermo-oxidative stability of CO and its epoxide was estimated by onset (under nitrogen atmosphere) and oxidative onset (under oxygen atmosphere) temperatures. TGA thermo grams of ( Fig. S9a and b) CO and its epoxide furnish the detailed information related to its degradation behavior in the presence of the inert gas (N 2 ) and reactive gas (O 2 ). Onset temperature is reported as the first perceptible temperature at which the degradation of CO and its epoxides starts. From Figure S9a, it can be seen that CO and its epoxide are stable up to 310 and 308°C (onset temperatures), respectively, in an inert atmosphere. Similarly, the maximum decomposition temperature represents the temperature at which utmost weight loss of CO and epoxide samples occurred, that is, 371.5 and 350.5°C, respectively. Lower thermal stability of epoxide signifies the degradation of epoxide at lower temperature compared to the unmodified CO. This may be due to the presence of high acid media in the mixture which is responsible for the lower thermal stability of CO epoxide [54].
Park et al. [57] reported that oil with saturated fatty acids, and mono unsaturation content has a positive influence and thermally more stable than polyunsaturation. Since the double bonds are converted to epoxide, hence higher thermal stability of epoxide was anticipated compared to unmodified CO, but due to the aforementioned reason, lower thermal stability was noticed. Study on thermal stability of epoxide is of great importance for various applications, and it depends on the chemical composition and structure of the epoxide.  Similarly, oxidative stability was found as a quality indicative parameter under oxygen atmosphere. In the present communication, oxidative stability was specified as the resistance of the epoxide against oxidation in the air atmosphere. Figure S9b describes the TGA plots of CO and its epoxide, from the oxidative onset temperature of the thermogram, it was found to be 320°C for both modified and unmodified CO (Fig. S9b). Similarly, the maximum decomposition temperature was found to be 382 and 354°C, respectively, for CO and its epoxide. In addition to the aforementioned elucidation, this attitude can be justified by the presence of unsaturation in the epoxide sample. Finally, it can be concluded that the CO showed more thermo-oxidative stability than the structurally modified CO epoxide.

Viscosity and viscosity index
In the current study, kinematic viscosities of CO and its epoxide were found to be 193.13 and 249.84 cSt, respectively at 40°C (Table 3). Increase in the epoxide molecular weight can be attributed to the addition of an oxygen molecule at the unsaturation sites thereby forming the oxirane ring [45]. The results of the study depict that CO epoxide viscosity (249.84 cSt) has been improved over the unmodified CO, and this signifies that the synthesized CO epoxide can serve as an acceptable lubricant base stock with enhanced viscosity to reduce friction. The ability of substance to resist free flow is one of the highly significant attributes for many heavy duty and industrial materials such as fuels, lubricants, and surfactants [58]. Therefore, it is desirable that the viscosity must be high enough all the time to keep good oil film between the moving parts to reduce friction [58]. Otherwise, due to loss in the lubricant, there is a tendency to increase the friction thereby resulting in power loss and rapid wear on the machine parts [58]. Viscosity is one of the crucial parameters while selecting a lubricant for a specific application, and failure to use the right lubricant with the required viscosity results in extreme temperature that may result in poor lubrication, equipment failure, and damage [59]. As discussed above, viscosity of acceptable range cannot be covered by conventional vegetable oils, hence in this study, structural modification was attempted to enhance the viscosity and molecular weight. However, viscosity can also be improved by addition of a viscosity enhancer [60].
Similarly, it is also desirable to check the viscosity in terms of the viscosity index (VI) which indicates the lower sensitivity of viscosity variations at higher temperatures. During this study, CO and its epoxide VI was found to be 99.52 and 188.92, respectively. Higher VI of CO epoxide signifies that the addition of an oxygen molecule at the unsaturation sites leads to an overall increase in the molecular weight of the final epoxide product and the viscosity index [14]. The higher VI of epoxide (188.92) almost doubles compared with the viscosity index of castor oil (99.52) indicating that the prepared epoxide can act as a high-temperature lubricant basestock.

Physico-chemical characterization
Significant physico-chemical properties of CO and its epoxide are estimated and given in Table 3. Acid value (AV) of CO and its epoxide was found to be 45.6 and 1.46 mg KOH/g, respectively. Lower acid value of CO epoxide signifies the smooth operation and functioning of the equipment during its usage. FFA is always considered as half of the AV, which signifies the formation of soap when mixed with water. The obtained FFA value was very less (0.73 mg KOH/g) and signifies trouble-free performance of the epoxide. The density of CO was found to be 790.74 kg/m 3 and corresponding epoxide density was 837.24 kg/m 3 . Improved density of CO epoxide was assigned to increase molecular weight of epoxide by addition of an oxygen molecule in the midst of unsaturation sites. Determination of the IV (i.e., quantity of the double bonds) after epoxidation reaction is one of the best ways to confirm and support the completion of the epoxidation reaction. Initial IV of CO was found to be 89.69 (g I 2 /100 g) whereas after epoxidation the value was 51.86 (g I 2 /100 g). Higher IV of epoxide signifies incomplete epoxidation reaction and the same has been observed from 1 H NMR and TLC spectral analysis ( Table 3). The presence of hydroxyl groups was confirmed by α-Glycol content analysis. The theoretical and experimental values of αglycol content along with the relative percentage conversion of αglycol is calculated and reported in Table 3.
Another significant property of epoxide is moisture content indicating the presence of water in the epoxide. Presence of moisture in the epoxide supports the bacterial growth which leads to an undesired performance upon usage thereby increasing AV, viscosity, and formation of free radical compounds via oxidation [21]. In this study, the CO epoxide moisture content was found to be 0.21 wt% (Table 3), which indicates safe performance of epoxide during usage. Refractive index for both the samples was found to be almost similar 1.477 and 1.479 (Table 3), which conveys that a very small quantity of heat energy can pass through the CO and epoxide samples, which helps to avoid the thermal degradation of end products during its usage and storage.

Conclusion
This is the first comprehensive report on the synthesis of castor oil epoxide from high FFA castor oil. The optimum condition for CO epoxidation occurred at temperature, 52.81°C; hydrogen peroxide to ethylenic unsaturation molar ratio, 1.65; catalyst loading, 15.14 wt %; and reaction time, 2.81 h. At this optimal condition, maximum epoxide content and relative conversion to oxirane was found to be 3.85mass% and 0.71%. Further, the product was confirmed by 1 H-NMR, 13 C-NMR, IV, and oxirane analysis. Quadratic polynomial model and ANOVA has well explained the interaction between the process variables. Further the model was examined and validated for the best fit. Comparative evaluation of physico-chemical properties of castor oil and its epoxide revealed surprising results. Among all the evaluated properties, only a few properties showed significant improvement, whereas thermo-oxidative stability and cold flow property results were unsatisfactory which may be due to high FFA content in CO or higher oxirane cleavage. However, the obtained viscosity index of epoxide indicates that the product can be used as hightemperature lubricant basestocks. However, from the present study, it can be concluded that low FFA CO could act as a suitable feedstock for lubricant synthesis with improved physico-chemical properties and anticipated to give improved thermo-oxidative stability and cold flow properties.

Supporting Information
Additional Supporting Information may be found in the online version of this article: Table S1. Experimental design matrix and results of oxygen oxirane content (OOC) of CO as affected by time, temperature, hydrogen peroxide molar ratio and catalyst loading. Table S2. Optimum reaction conditions for maximum OOC of castor oil. Figure S1. (a) Counter and (b) Response surface plots for the effect of the time (A) and temperature (B) on the maximum oxirane oxygen content. Figure S2. Mechanism of excess reagent's reaction with the epoxy rings to cause unwanted side reactions. Figure S3.      Figure S8. 1 H NMR spectrum of CO (a) and its epoxide (b). Figure S9. TGA thermogram for CO and its epoxide at 10°C/min, under inert (N 2 ) atmosphere (thermal stability) (a), (reactive) oxygen atmosphere (oxidative stability) (b).