Next Article in Journal
Improving Environmental Efficiency of Reverse Filling Cementitious Materials through Packing Optimization and Fiber Incorporation
Previous Article in Journal
Design and Synthesis of π-Extended Resveratrol Analogues and In Vitro Antioxidant and Anti-Inflammatory Activity Evaluation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimisation of Epoxide Ring-Opening Reaction for the Synthesis of Bio-Polyol from Palm Oil Derivative Using Response Surface Methodology

by
Norsuhaili Kamairudin
1,
Seng Soi Hoong
2,*,
Luqman Chuah Abdullah
1,3,*,
Hidayah Ariffin
1,4 and
Dayang Radiah Awang Biak
3,5
1
Higher Education Centre of Excellence (HiCoE), Institute of Tropical Forestry and Forest Product, University Putra Malaysia, Serdang 43400, Selangor, Malaysia
2
Malaysian Palm Oil Board, No. 6, Persiaran Institusi, Bandar Baru Bangi, Kajang 43000, Selangor, Malaysia
3
Department of Chemical and Environmental Engineering, Faculty of Engineering, University Putra Malaysia, Serdang 43400, Selangor, Malaysia
4
Faculty of Biotechnology and Biomolecular Sciences, Universiti Putra Malaysia, Serdang 43400, Selangor, Malaysia
5
Institute of Advanced Technology, Universiti Putra Malaysia, Serdang 43400, Selangor, Malaysia
*
Authors to whom correspondence should be addressed.
Molecules 2021, 26(3), 648; https://doi.org/10.3390/molecules26030648
Submission received: 16 October 2020 / Revised: 2 December 2020 / Accepted: 3 December 2020 / Published: 27 January 2021

Abstract

:
The development of bio-polyol from vegetable oil and its derivatives is gaining much interest from polyurethane industries and academia. In view of this, the availability of methyl oleate derived from palm oil, which is aimed at biodiesel production, provides an excellent feedstock to produce bio-polyol for polyurethane applications. In this recent study, response surface methodology (RSM) with a combination of central composite rotatable design (CCRD) was used to optimise the reaction parameters in order to obtain a maximised hydroxyl value (OHV). Three reaction parameters were selected, namely the mole ratio of epoxidised methyl oleate (EMO) to glycerol (1:5–1:10), the amount of catalyst loading (0.15–0.55%) and reaction temperature (90–150 °C) on a response variable as the hydroxyl value (OHV). The analysis of variance (ANOVA) indicated that the quadratic model was significant at 98% confidence level with (p-value > 0.0001) with an insignificant lack of fit and the regression coefficient (R2) was 0.9897. The optimum reaction conditions established by the predicted model were: 1:10 mole ratio of EMO to glycerol, 0.18% of catalyst and 120 °C reaction temperature, giving a hydroxyl value (OHV) of 306.190 mg KOH/g for the experimental value and 301.248 mg KOH/g for the predicted value. This result proves that the RSM model is capable of forecasting the relevant response. FTIR analysis was employed to monitor the changes of functional group for each synthesis and the confirmation of this finding was analysed by NMR analysis. The viscosity and average molecular weight (MW) were 513.48 mPa and 491 Da, respectively.

1. Introduction

Polyol is known to be a very important and beneficial material in pharmaceutical, food science and polymer chemistry [1,2,3]. Polyol is commonly used as raw materials in polymer chemistry in order to manufacture polyurethane products (coatings, elastomers and foam) by reacting with isocyanates [4,5,6]. Therefore, due to its versatility, polyurethane is extensively used in daily life and it has been predicted that the demand will expand from 14.2 million tons in 2011 to 22.2 million tons by 2020 [7]. However, the recent polyol industry has been heavily dependent on petroleum resources, which aggravates the depletion of inadequate natural resources. In view of the current emphasis on sustainability and global warming, the production of bio-based polymers from renewable resources has become of great importance [8].
Historically, vegetable oils such as soy bean oil [9], palm oil [10,11], canola oil [12] and jatropha oil [13] are the major renewable resources used for synthesising bio-polyol. For the best commercial plantation, palm oil has a yield potential of usually about 4 to 6 tons of oil per hectare per year. Approximately, 17.86 million tons of crude palm oil were produced in 2019 in Malaysia alone [14]. It has been predicted to increase to over 20 million tons in 2020. Therefore, palm oil and its derivatives are a great eco-friendly renewable resource for bio-based polyurethane materials due to their low cost and abundance of supply [10,15]. Palm oil consists of 53.4% saturated and 46.6% unsaturated fatty acids. The major fatty acids are palmitic acid, oleic acid and linoleic acid [10]. However, due to the absence of hydroxyl group in the structure of palm oil and its derivatives, they are unreactive under the condition of polyurethane chemistry. Therefore, they need to be functionalised with the introduction of the hydroxyl group in the molecular structure, so that they can be employed to produce polyurethane products. The introduction of the hydroxyl group can be undertaken at the position of alkene by various methods [16]. Example of methods employed are ozonolysis and hydrogenation [17], epoxidation followed by ring opening [18], thiol-ene coupling [19], and microbial conversion [20].
Methyl oleate (MO) is one type of fatty acid methyl ester (FAME), also recognised as biodiesel, derived from the transesterification of fats or vegetable oils such as palm oil, soybean oil and rapeseed oil. Due to the abundant availability of palm oil, this situation drives Malaysia towards the development of biodiesel technology and production. The current biodiesel production capacity in Malaysia is about 10.2 million tons [21]. However, the Malaysian biodiesel industry faces some challenges due to insignificant public support and escalating nontariff barriers, which indirectly influence the demand on for palm oil biodiesel [22]. Diversification of methyl oleate (MO) into other chemical products would assist to decrease the excess supply and capacity of the biodiesel industry. Thus, functionalisation of methyl oleate by introducing the hydroxyl functional group may develop new and useful bio-based products such as bio-polyol for use in the production of polyurethane products.
As reported by Hazmi et al., [13], the epoxidation of jatropha oil followed by a hydroxylation process using methanol, distilled water and formic acid as a catalyst generates polyol with a hydroxyl value of 179 mg KOH/g. The same chemical modification can be seen through the study by Tuan Ismail et al., [21], in which the hydroxyl value for the FAME polyol synthesised via epoxidation and ring opening process with ethylene glycol and boron trifluoride diethyl ether complex was about 189 mg KOH/g, which limits its use in some polyurethane applications such as rigid polyurethane foam. Based on extensive studies in the literature concerning the synthesis of bio-polyol through epoxidation followed by ring opening reaction, there has been no report on the optimisation of the reaction condition for this process as a function of the hydroxyl value (OHV). Therefore, the optimisation process for the hydroxyl value can be carried out using a competent method to evaluate the effect of the individual factors of each reaction condition on the overall effectiveness. One of the methods is Response Surface Methodology (RSM) that was first introduced by Box and Wilson in the 1950s. RSM is a mathematical technique that has been widely used in various areas of industry and academic research in order to optimise experimental conditions [23]. In the normal technique, one factor at a time is employed to investigate the effect of the operational parameters on the respective response. However, this RSM tool can be used to discover the interaction between independent variables (factors) and dependent variables (responses) involving statistical analysis and mathematical modelling with a minimal number of experiments. Thus, it gives an indication of the optimal condition of the process [24,25].
The aim of this present study is to develop the RSM model to optimise the hydroxyl value of bio-polyol from methyl oleate. To the best of the knowledge of the authors, there is no reported study regarding the optimisation of the reaction condition for the epoxide ring opening reaction using RSM, even though it is used as a statistical tool for the optimisation technique in other chemistry areas such as extraction of flavonoids [23] and waste water treatment [26].

2. Results and Discussions

2.1. Epoxidation of Methyl Oleate (EMO)

Reaction Mechanism of Epoxidation of Methyl Oleate (EMO)

Scheme 1 illustrates the reaction mechanism of epoxidation of methyl oleate reaction. Step 1 shows that, after the addition of hydrogen peroxide, the performic acid was formed immediately. However, the performic acid was unstable and reactive. Thus, to prevent an extremely vigorous reaction, hydrogen peroxide was added dropwise. The oxirane ring (Step 2) was formed when the alkene group of methyl oleate reacted with performic acid, which yielded epoxidised methyl oleate. In order to get the highest oxirane conversion in the shortest time, excess hydrogen peroxide was needed to be used [27]. Therefore, in this present study, the mole ratio of EMO/formic acid/hydrogen peroxide was set at 1:1:2.5.
Figure 1 shows the Oxirane Oxygen Content (OOC) profile for epoxidation reaction of methyl oleate (MO). This epoxidation reaction was completed within 3 h. The maximum OOC recorded in this reaction was 5.30%. However, after the washing and neutralisation processes, the final OOC recorded was 5.10%. This is attributable to the ring opening of some of the epoxide groups [21]. The final percentage of conversion was 93.38%.
Figure 2 illustrates the comparative FTIR spectra for the MO and EMO. The chemical changes from MO to EMO can be clearly seen in the spectra. In the original state, the double bond (C=C) of the methyl oleate (MO) was detected as a small peak at 3007 cm−1. The peaks of the CH2 stretching vibration were presented at 2923 cm−1 and 2854 cm−1. The appearance of the characteristic ester carbonyl (C=O) was at 1742 cm−1, which confirmed the existence of the ester group in MO [21,28]. The peak of the ester functional group (C-O) appeared in the range of 1170 to 1240 cm−1 [21,29]. The peak at 722 cm−1 corresponded to the existence of the aliphatic chain. Upon the epoxidation reaction, the double bond of MO disappeared, and a new peak could be observed at 827 cm−1, which confirmed the formation of the epoxide group (C-O-C). According to Rayung et al., [30], the epoxide ring is usually detected at a bandwidth from 820–830 cm−1.
Figure 3 illustrates the proton (1H) NMR spectra for MO and EMO. The molecular structure of MO and EMO were revealed using 1H-NMR analysis. In Figure 3A, the signal at 5.3 ppm (peak g) corresponds to the unsaturated double bond (C=C) that is present in MO. According to Sathai et al., [31], an unsaturated double bond signal can be seen at 5.2–5.41 ppm whereas signals for the long chain aliphatic (CH2) and terminal CH3 would appear at 1.2–1.4 (peak b) and 0.8–0.9 ppm (peak a), respectively. During the epoxidation process, the double bond (C=C) in MO was converted to an epoxide ring (C-O-C) in EMO. As shown in Figure 3B, a new peak of the epoxide ring (C-O-C) was observed at 2.88 ppm (peak h), which confirmed the formation of the epoxide group. This fact was also supported by the FTIR spectrum when the epoxide peak appeared at 827 cm−1.

2.2. Ring Opening of Epoxidised Methyl Oleate with Glycerol

2.2.1. Reaction Mechanism of Epoxide Ring-Opening by Glycerol

Scheme 2 illustrates the reaction mechanism of epoxide ring-opening by glycerol. In the first step, the oxygen atom of epoxide was protonated by a Lewis acid (BF3·Et2O) catalyst and generated a good leaving group (Step 1). After that, the carbon-oxygen bond began to break (Step 2) and positive charge started to build up on the more substituted carbon. The nucleophile (from glycerol) attacked the electrophile carbon (Step 3). Then, the oxygen was deprotonated (Step 4) to form the MOG-polyol product.

2.2.2. Single Factor Experiment

Figure 4 illustrates the results of the hydroxyl value (OHV) obtained at different mole ratios of (EMO:Glycerol), percentages of catalyst and reaction temperature.

Effect of Mole Ratio (EMO:Glycerol) on Hydroxyl Value (OHV)

In this experiment, the effect of mole ratio between EMO and glycerol on the hydroxyl value (OHV) of MOG-polyol was examined. Other reaction parameters were kept constant: catalyst loading 0.5%, reaction temperature at 90 °C. The mole ratio EMO:Glycerol was evaluated in the range between 1:5–1:10. The result obtained is shown in Figure 4A. The analysis result showed that the hydroxyl values of polyols increased with increasing mole ratio of EMO:Glycerol. The reason is that with a higher concentration of nucleophilic reactant (glycerol), the epoxide ring-opening reaction was favoured compared to the oligomerisation which leads an increase of OHV [32]. However, in this study, the mole ratio of EMO:Glycerol was not increased to more than 1:10 due to the high usage of glycerol as a reactant.

Effect of Catalyst on Hydroxyl Value (OHV)

In this part, the effect of the catalyst on the hydroxyl value (OHV) of MOG-polyol was studied. Other reaction parameters were kept constant, in which the mole ratio of EMO:Glycerol was 1:10, the reaction temperature was 90 °C and the catalyst loading was varied in the range of 0.125–1.0%. The result obtained is shown in Figure 4B. The data showed that the hydroxyl values of the polyols increased with decreasing amount of catalyst. Specifically, the hydroxyl value (OHV) increased from 208 mg KOH/g to 287 mg KOH/g as the catalyst decreased from 1.0–0.125%.

Effect of Temperature on Hydroxyl Value (OHV)

In this experiment, the effect of reaction temperature on the hydroxyl value (OHV) was explored, while the other reaction parameters were kept constant (mole ratio of EMO:Glycerol at 1:10, catalyst loading = 0.5%). The temperature was varied in the range of 60–150 °C. The result obtained is shown in Figure 4C. The obtained result showed that the hydroxyl values (OHV) of the polyols increased with increasing reaction temperature until 120 °C, beyond which the hydroxyl value (OHV) showed a slight reduction.

2.2.3. Optimisation Reaction Conditions by Response Surface Methodology

Fitting the Model and Analysis of Variance (ANOVA)

Response surface modelling was applied to find the optimum reaction conditions for independent variables in order to obtain the maximum OHV for MOG-polyol. The experiments were programmed as per the CCRD layout using (Design Expert Software version 7.0, Stat-Ease, Minneapolis, MN, USA) containing 20 experiments with different combinations of independent variables. The design matrix of the actual and predicted OHV is shown in Table 1. Each experimental run was conducted in triplicate in order to obtain consistent results. The coefficients regressions were calculated based on CCRD and the results of the experiments. The second-order polynomial equation was established to obtain the hydroxyl value (OHV) as in Equation (1):
Y (Hydroxyl value) = + 280.83 + 16.43A − 18.15B + 2.05C + 7.08AB − 4.0AC − 0.81BC − 0.11A2 − 6.53B2 − 2.18C2
where A is the symbol of the mole ratio of EMO:Glycerol, B is the symbol of the catalyst, C is the symbol of the temperature and Y is the response (Hydroxyl value). A positive sign in the equation showed a synergistic effect. A negative sign indicated as an antagonistic effect [33,34].
Sometimes, a mathematical model developed after a fitting function to the data, may not satisfactorily express the experimental domain studied [33]. Therefore, a more reliable way to investigate the adequacy of a model and goodness fit for the response is by employing the analysis of variance (ANOVA) [35]. ANOVA analysis is also needed in order to discover whether the second-order polynomial model is significant or not.
The quadratic polynomial model is illustrated in Table 2. The p-value (<0.0001) indicated that the developed model was significant and adequate to demonstrate the actual interaction between the significant variables and the response (hydroxyl value) [36]. The lack of fit (F-value = 4.18) showed insignificant data, which was comparative to the pure error of the experiments, which indicated that the quadratic model data were convenient to depict the whole experimental data. Noise in the experiment contributed to the large value of lack of fit [37]. The p-value was also used to distinguish the significance of every coefficient. A small p-value and large F-value would mean that the independent variables have a large impact and would be significant for represented response [38]. In addition, when p-value has a small value, the corresponding coefficient will be more significant. Table 2 illustrates that the coefficients of A, B, C, AB, AC, B2 and C2 were identified as significant in model terms due to the p-value being below 0.05. From the model, catalyst (B) (F-value = 453.02; p-value < 0.0001) was the most significant variable followed by the mole ratio (A) and temperature (C). However, from this model, BC and A2 did not show any significant effect, thus indicating that there is no mutual interaction between catalyst and temperature (BC) and mole ratio (EMO:Glycerol), (A) multiplied by itself. Although these model terms were insignificant, they were kept in the model equation in order to support the hierarchy. As the model design could accurately predict up to a 2% residual standard error (RSE) according to the model validation that meant the insignificant terms still contributed to the model. Ordinarily, insignificant terms are not eliminated from the equation, but may be retained depending on their suitability in the model.
The coefficient of determination (R2) can be a useful tool to make sure of the accuracy of the quadratic polynomial model as shown in Table 3. A high coefficient determination (R2 = 0.9897) was presented by the polynomial model, which suggested that 98% of the variation in the hydroxyl value (OHV) response could be described by the independent variables and only 2% of the variation could not be described by the model. Additionally, the predicted R2 (0.9296) was also in reasonable agreement with the adjusted R2 (0.9804), with a difference of less than 0.2 [39]. Therefore, the relationship between the predicted and experimental values of hydroxyl value (OHV) response was proven to have a high R2-value, which suggested that the experimental values fitted well with the predicted value [40].
In addition, the adequate precision demonstrated that a signal-to-noise ratio of greater than 4 is desirable [41]. The model showed that the result of an adequate precision ratio of hydroxyl value (OHV) response was 37.588, which indicated that the model could be employed to navigate the design space. For the model to be adequate, the value of the coefficient of variation (CV) result needed to be less than 10%. Therefore, as the study indicated a very low value (1.15%) for the CV, this showed that the experimental value had very good reliability and high precision [42].
The residual plot is a crucial detection method to discover the systematic departures from the assumption employed in the regression equations [43]. As shown in Figure 5A, it confirmed the independence of residuals and the assumption of a normal distribution, as the data points on this plot were close to a straight line. Figure 5B illustrates the correlation between the experimental and predicted values for the hydroxyl value (OHV) response. The actual values represented the measured experimental response data for a specific experimental run, while the predicted values were the prediction results generated by the model. As shown in Figure 5B, the actual values were correlated to the predicted values with an R2 of 0.9897. Therefore, this model predicted accurately to determine the correlation between the independent variables and the hydroxyl value (OHV).

Mutual Effects of Process Parameters

The correlation between the independent variables provided significant impact on the hydroxyl value (OHV) response as shown in Equation (1). Instead of investigating one variable individually, a study about the correlation among variables could be done to determine their significance for more sophisticated optimisation. According to Mangili et al., [44], 2D (contour plot) and 3D (three dimensional graphs) of RSM model surfaces were illustrated as graphical presentation of the regression equation in order to have a better understanding of the interaction effect between two variables at a time within the range considered. Each 2D and 3D curve in Figure 6 depicts an infinite number of combinations of two significant variables while the other variables were kept at a fixed position at central level [41].
Figure 6A represents the interaction between the mole ratio of EMO to glycerol (A) and catalyst (B) on the hydroxyl value (OHV) response at a fixed temperature (C) of 120 °C. The OHV was positively correlated to the mole ratio (A), in which the OHV increased with a higher mole ratio. This was due to the high concentration of hydroxyl group of glycerol at a higher mole ratio. However, when the reaction was conducted with a higher loading of catalyst (B), the results showed lower OHV values. This may be due to the higher catalyst loading that promoted oligomerisation between the fatty ester chains to take place, which contributed to the lower hydroxyl values observed [45].
Figure 6B illustrates the relationship between the mole ratio (A) and temperature (C) as a function of the hydroxyl value (OHV) response at a fixed catalyst (B) of 0.35%. When the mole ratio of EMO:Glycerol was increased from 1:5 to 1:10, it caused the hydroxyl value to increase correspondingly. The same effect was observed in Figure 6A. The OHV increased with higher reaction temperature up to 120 °C, which was believed to be the optimum temperature. This might be due to deactivation of the catalyst (BF3·Et2O) because its boiling point is at 130 °C. Therefore, the catalyst was inefficient at the temperature above 130 °C.
Figure 6C illustrates the 2D and 3D plots of the correlation between catalyst (B) and the reaction temperature (C) on the OHV response at a fixed mole ratio (A) of (1:7.5). It was observed that the hydroxyl value (OHV) increased with decreasing catalyst loading such as when seen in the results shown in Figure 6A. For reactions conducted at temperatures in the range of 120–125 °C, the hydroxyl value (OHV) obtained ranged from 259–286 mg KOH/g.

2.3. Verification of the Model

A verification test was carried out to evaluate the validity of the model. In this part, three sets of experiment were produced from the Design-Expert Software (version 7.0, State Ease Inc., Minneapolis, MN, USA), in order to make a comparison between the predicted values and the experimental values. The maximum hydroxyl value (OHV) response was prescribed as a main target, while the other variables remained in the range studied (Table 4). Table 5 illustrates the actual and predicted responses from the optimum combination of the independent variables. For all experiments, the results showed a low residual standard error (RSE), which was less than 2%, thus indicating that the model was valid and could predict the OHV response up to 98%. RSE was calculated based on Equation (2). From Table 5, the highest OHV was obtained from Experiment 3. Thus, for this reason, Experiment 3 was chosen as the optimum reaction condition [46]. Therefore, based on CCRD-RSM, an optimised reaction condition with (1:10) mole ratio of (EMO:Glycerol), (0.18%) of catalyst and temperature of 120 °C was suggested, which would produce a polyol with hydroxyl value (OHV) of 306 mg KOH/g. Table 6 shows the optimised reaction conditions for the OHV response of the MOG-polyol with desirability about 0.976. A better-generated solution is given when the desirability is close to 1. Therefore, the optimum reaction condition for the maximum hydroxyl value (OHV) was successfully generated by the CCRD-RSM model.
R e s i d u a l   S t a n d a r d   E r r o r   R S E   % =   E x p e r i m e n t a l   v a l u e P r e d i c t e d   v a l u e P r e d i c t e d   v a l u e   × 100 %

2.4. Evaluation of the Importance of the Variables on the Hydroxyl Value (OHV) of Polyol

The important influence of efficient variables on the hydroxyl value was demonstrated by Pareto analysis as shown in Figure 7. The Pareto analysis expresses the influence of each variable, in percentages, on the response in relation to the coefficient of the coded variable in Equation (3) [33,47]. The pareto analysis showed that variable B (catalyst) gave the most important influence on the response (hydroxyl value) at 45.89% followed by variable A (mole ratio of EMO:Glycerol) at 37.61% and the interaction between A and B at 6.98%.
P i = b i 2 b i 2 × 100   i 0

2.5. Characterisation of Optimised MOG-Polyol

As shown in Figure 8, the disappearance of the epoxy group (C-O-C) vibration stretching at 827 cm−1 [47] for the MOG-polyol spectrum supported the fact that epoxide ring-opening reaction with glycerol had occurred. The reaction generated the hydroxyl group that was detected as a broad O-H peak at 3391 cm−1 that can be seen clearly in the FTIR spectrum [48]. The presence of this broad peak in this spectrum confirmed that the OH functional group had been added onto the epoxidised methyl oleate structure that originated from the glycerol molecule and the epoxide group.
Figure 9 presents the 1H-NMR spectrum of the MOG-polyol, which confirmed that the peak associated with the proton of the epoxide group (C-O-C) at 2.88 ppm was not observed. This observation indicated that the epoxide ring-opening reaction with glycerol had occurred. Upon occurrence of the ring-opening reaction, new proton peaks were detected at 3.03–3.06 ppm, 3.41–3.48 ppm and 3.64–3.81 ppm, showing the presence of the OH groups in the MOG-polyol. This is supported by the appearance of a broad peak at 3319 cm−1 in the FTIR spectrum, which proved the presence of OH groups in the polyol.
The iodine value of the EMO (Table 7), which quantifies the alkene group, was significantly low in comparison with the MO due to epoxidation of alkene to the epoxide (C-O-C) group. As shown in Scheme 2, MOG-polyol contains both the primary and secondary OH group in its molecular structure, leading to a higher molecular weight as compared to EMO and MO (Table 7). Additionally, the hydroxyl groups of MOG-polyol could be quantified by wet chemistry analysis and were about 300 mg KOH/g. This result is in agreement with the finding by Prociak et al., [49]. The viscosity of the MOG-polyol was higher compared to EMO and MO due to the higher hydroxyl value (OHV) and molecular weight (MW) of the polyol [50].

3. Materials and Methods

3.1. Materials

Methyl oleate from palm oil with purity 98% and iodine number 87g I2/100g was purchased from Carotino Sdn Bhd (Pasir Gudang, Johor, Malaysia). Hydrogen peroxide (H2O2) (50% concentration in H2O), formic acid (94% concentration) and glycerol (99.5%) were purchased from (R&M Chemicals, Cennai, India). Boron trifluoride diethyl ether complex (BF3·Et2O) (50% purity) was purchased from Merck Sdn Bhd (Petaling Jaya, Selangor, Malaysia). Sodium chloride (NaCl), sodium carbonate (Na2CO3) and sodium hydrogen carbonate (NaHCO3) were purchased from Chem AR (Systerm Chemicals, Shah Alam, Malaysia), respectively. All purchased chemicals were used as received without further purification unless otherwise noted.

FTIR and NMR Data of Methyl Oleate

υmax/cm−1 3007 (C=C), 2924, 2854 (CH2), 1742 (C=O), 1459 (CH2), 1436 (CH3), 1170 (C-O), 723 (CH2); 1H-NMR (600 MHz, CDCl3): δH = 5.30–5.24 (2H, m, CH2CH=CHCH2), 3.58 (3H, s, OCH3), 2.27 (2H, t, O=CCH2CH2), 2.00–1.94 (4H, m, CH2CH=CHCH2), 1.58–1.52 (2H, m, O=CCH2CH2CH2), 1.23–1.19 (20H, m, CH2CH2CH2), 0.81 (3H, t, CH2CH3); 13C-NMR (150 MHz CDCl3): δC = 174.1 (CH3OC=O), 130.1–129.6 (CH2CH=CHCH2), 51.3 (CH3O), 34.1 (O=CCH2CH2), 31.9–29.1 (CH2CH2CH2), 27.2–27.1 (CH2CH=CHCH2), 24.9 (O=CCH2CH2CH2), 22.6 (CH2CH3), 13.9 (CH2CH3).

3.2. Epoxidation of Methyl Oleate

The epoxidation of methyl oleate was performed in a 5-L water-jacketed glass reactor, equipped with a mechanical agitator, a thermometer couple, a dropping funnel and reflux condenser. Methyl oleate (1000 g) was heated to 40 °C. Performic acid, a premixed reagent of 50 wt% hydrogen peroxide and formic acid of the desired ratio was slowly added dropwise to the preheated methyl oleate. The mole ratio of the methyl oleate/formic acid/hydrogen peroxide was 1:1:2.5. Oxirane oxygen content (OOC) analysis was used to monitor the progress of the epoxidation reaction. Sampling was done every 30 min until the OOC of the reaction mixture reached more than 90% of the theoretical value, i.e., 5.47. In this study, the time taken for the epoxidation reaction to reach this point was 3 h. This method was adopted from Arniza et al., and Tuan Ismail et al., [10,21]. The epoxidised methyl oleate (EMO) was then separated from the spent acid and neutralised by washing with 1% sodium chloride solution and 0.5% sodium hydrogen carbonate solution until the pH of the EMO was about 6–7. Then, it was dried under vacuum for several hours at 80–90 °C. The white liquid product with a yield of 90 wt% was denoted as EMO. υmax/cm−1 2924, 2855 (CH2), 1739 (C=O), 1463 (CH2), 1436 (CH3), 1170 (C-O), 827 (C-O-C), 723 (CH2); 1H-NMR (600 MHz, CDCl3): δH = 3.62 (3H, s, OCH3), 2.87–2.84 (2H, m, CH2CHOCHCH2), 2.26 (2H, t, O=CCH2CH2), 1.61–1.56 (2H, m, O=CCH2CH2CH2), 1.48–1.44 (4H, m, CHOCH2CH2), 1.31–1.21 (20H, m, CH2CH2CH2), 0.84 (3H, t, CH2CH3); 13C-NMR (150 MHz CDCl3): δC = 174.3 (CH3OC=O), 57.2–56.7 (CH2CHOCHCH2), 51.4 (CH3O), 34.1 (O=CCH2CH2), 31.9–29.0 (CH2CH2CH2), 27.8 (CH2CHOCHCH2), 24.9 (O=CCH2CH2CH2), 22.7 (CH2CH3), 14.1 (CH2CH3).

3.3. Ring-Opening of Epoxidised Methyl Oleate with Glycerol

The reaction synthesis was carried out in a 1-L three-neck flask equipped with a magnetic stirrer and dropping funnel. An amount of 100 g of EMO was placed in the flask and preheated to a temperature of (90–120 °C) at a stirring speed of 200 rpm. A mixture of catalyst, Boron trifluoride diethyl ether complex (BF3·Et2O) (0.15–0.55%) and glycerol was added slowly to the preheated EMO. The mole ratio of the EMO to glycerol was (1:5–1:10). The reaction was started at the time when the mixture of glycerol and catalyst was added to the preheated EMO. The reaction mixture was sampled every hour to monitor the progress of reaction through the oxirane oxygen content (OOC) analysis. The reaction was considered completed when the OOC of the EMO was less than 0.1%. Then, the mixture component was poured into the separating funnel and left to separate into two layers. The bottom layer that contained excess glycerol was removed. The component was cooled to room temperature and ethyl acetate (100 mL) was added to dilute the component. Sodium chloride (200 mL) and sodium hydrogen carbonate (200 mL) was used in sequences to wash and neutralize the component until the pH was about 6–7. The component was dried over MgSO4. Ethyl acetate was removed by evaporation under vacuum to yield about 90 wt% as a yellowish liquid denoted as MOG-polyol. υmax/cm−1 3391 (O-H), 2923, 2854 (CH2), 1740 (C=O), 1460 (CH2), 1436 (CH3), 1171 (C-O), 723 (CH2); 1H-NMR (600 MHz, CDCl3): δH = 3.81–3.78 (1H, m, CH2OCHOHCH2OH), 3.66–3.64 (1H, m, CHOHCH2OH), 3.59 (3H, s, OCH3), 3.48–3.44 (1H, m, CHOCH2CHOH), 3.43–3.40 (1H, m, CH2CHOHCHO), 3.06–3.03 (1H, m, CHOHCHOCH2), 2.23 (2H, t, O=CCH2CH2), 1.55–1.52 (2H, m, O=CCH2CH2CH2), 1.47–1.35 (2H, m, CHOCH2CH2), 1.31–1.18 (20H, m, CH2CH2CH2), 0.80 (3H, t, CH2CH3); 13C-NMR (150 MHz CDCl3): δC = 174.5 (CH3OC=O), 84.2 (CHOHCHOCH2), 73.5 (CH2CHOHCHO), 72.7 (CHOHCH2OH), 71.8 (CHOCH2CHOH), 71.1 (CH2OCHOHCH2OH), 51.5 (CH3O), 34.1 (O=CCH2CH2), 33.2 (CH2CH2CHOH), 31.9–29.0 (CH2CH2CH2), 24.8 (O=CCH2CH2CH2), 22.6 (CH2CH3), 14.1 (CH2CH3).

3.3.1. Single Factor Experiment

Single factor experiments were performed to explore the effects of the following parameters on the hydroxyl value (OHV) response, namely mole ratio (EMO:Glycerol) (1:5, 1:7.5 and 1:10), catalyst (0.125%, 0.25%, 0.50%, 0.75% and 1.0%) and temperature (60 °C, 90 °C, 120 °C and 150 °C). A series of studies was then performed afterwards in order to find out the best reaction condition on hydroxyl value (OHV) using Response Surface Methodology (RSM).

3.3.2. Experimental Design and Optimisation via Response Surface Methodology (RSM)

Central composite rotatable design (CCRD) of RSM was employed to optimise the process parameter conditions of the three independent variables, namely mole ratio (EMO:Glycerol) (A), catalyst (B) and temperature (C). A summary of the independent variables and their coded level is illustrated in Table 8. This design was used to evaluate the combined effect of these variables towards the hydroxyl value (OHV) response (Y) of the MOG-polyol. A total number of 20 experiment runs (Table 9) based on CCRD was developed to identify the optimised levels of the significant independent variables and their relationship between independent variables generated by the Design-Expert version 7.0 Software (State Ease Inc., Minneapolis, MN, USA). These three independent variables were studied at three different levels. Therefore, the total 20 experiment runs involved 6 axial points, 6 central points and 8 factorial points, gained using Equation (3):
N = 2n + 2n +nc = 23 + 2(3) + 6 = 20
where N is the total number of experiments run, n is number of factors and nc is the number of replicates in the central points. The experiments were run randomly to minimise the influence of the extraneous factors [38].

3.3.3. Statistical Analysis

Analysis of variance (ANOVA) was utilised to analyse the output data. By using ANOVA, it is a more reliable option to evaluate the quality of the fitted model. ANOVA is an important analysis in order to study the significance of the second order polynomial in the model. The second order polynomial regression model was used to fit the experimental data, and can be expressed using Equation (5):
Y = β 0 +   i = 1 k β i x i +   i = 1 k β i i x i 2 +   i = 1 k 1 j > i k β i j x i x j
where Y is the dependent variable (response); β0 is the constant coefficient; βi, βii and βij are coefficients of linear, quadratic and interaction terms, respectively; k is the number of tested variables (k = 4); xi and xj are the independent variables.

3.4. Verification of the Model

The verification of the final model was performed in order to inquire the adequacy of the predicted values. Three sets of the experiment data were suggested from the design software. The actual and predicted values were employed to calculate the percentage of residual standard error (RSE%), at less than 2%.

3.5. Characterisation

Wet chemistry analyses were performed following the American Oil Chemists’ Society Official (AOCS) methods [51]. The oxirane oxygen content (OOC) analysis of epoxidised methyl oleate (EMO) and methyl oleate-polyol (MOG-polyol) were determined according to AOCS Official method Cd 9-57. This analysis was conducted to measure the epoxide content of both samples. The iodine value (IV) was analysed following the AOCS Official method Cd 1d-92, in order to quantify the alkene group present in a sample and it was conducted on the MO and EMO. The acid value is defined as the amount of potassium hydroxide in milligrams needed to neutralise the free fatty acid in 1 g of sample. It was conducted according to the AOCS Official method Te 2a-64. The hydroxyl value (OHV) was determined following the AOCS Official method Cd 13-60, in order to quantify the hydroxyl group present in a sample. This analysis was performed on all the samples.
The infrared (FTIR) spectra of all samples, namely MO, EMO and MOG-polyol were investigated using a Perkin-Elmer FTIR Spectrum100 (Walthman, MA, USA) with an attenuated total reflectance (ATR) accessory. The FTIR was performed in the frequency range of 4000–650 cm−1 for 16 repeated scans with a resolution of 4 cm−1 in transmittance mode at room temperature. FTIR Spectroscopy was employed in order to identify the functional group that existed in the samples.
Nuclear magnetic resonance spectroscopy (1H-NMR) and (13C-NMR) was conducted using a JOEL JNM-ECZ600R (Peabody, MA, USA) at 600 MHz and 150 MHz, respectively. Deuterated chloroform was used as a solvent to dissolve the sample with approximately 10% w/w of sample. NMR Spectroscopy was utilized to identify the functional groups and the molecular structure of the samples. NMR analysis was performed on all samples, MO, EMO and MOG-polyol.
Gel permeation chromatography (GPC) was carried out using a Varian PL-GPC 50 Plus, (Polymer Laboratories Ltd., Church Stretton, Shropshire, UK) equipped with a differential refractive index (DRI)/viscometer and an autosampler combined with a detector. The sample was dissolved in THF with a concentration of 2 mg/mL. A PL gel mixed D column with THF as an eluent at a flow rate of 1 mL/minute was used to analyse the samples. The molecular weight was measured based on the universal calibration curve created using polystyrene standards with range from 162 to 1 × 105 Da. Cirrus GPC/SEC system PL-GPC 50 plus (Shropshire, UK) was used to analyse the data. GPC was utilised to determine the molecular weight distribution and average molecular weight of a sample. This analysis was performed for all samples.
The viscosity of a sample was measured by using a Brookfield Digital Rheometer, Model DV-111+ (Brookfield Engineering, Middleborough, MA, USA) at 25.0 °C ± 1. Approximately, 1 mL of sample was put on the sample platform for testing. A flat plate spindle (PP41) was moved onto the sample in rotation mode with a gap of 1 mm. The viscosity was measured using plate-plate geometry mode and analysed using Rheo Min Software (version 5.0, Middleborough, MA, USA). Viscosity analysis was performed on all samples.

4. Conclusions

The optimum reaction conditions for the synthesis of bio-polyol (MOG-polyol) by applying Response surface methodology (RSM) with a combination of central composite rotatable design (CCRD) with three variables at five levels (−2, −1, 0, +1, +2) was successfully carried out. The effect of the three main variables was investigated using an empirical model and it showed a linear effect on the response (hydroxyl value). Based on the ANOVA analysis, the predicted model was fitted, with the p-value > 0.0001 and F-value (106.47), and an insignificant lack of fit and a high coefficient of determination (R2) = 0.9897. The verification of this model was achieved by testing three experiments suggested by RSM, and the RSE obtained was less than 2%. The optimised reaction condition for synthesis MOG-polyol was 1:10 (mole ratio of EMO to glycerol), 0.18% of catalyst and 120 °C (temperature) which resulted a 306.19 mg KOH/g of hydroxyl value (OHV) and a desirability of 0.976. The molecular structure of the MOG-polyol was confirmed by FTIR and NMR analysis. Therefore, from this study, the result obtained showed that RSM is an excellent tool to optimise the reaction conditions for synthesising the bio-polyol (MOG-polyol) from methyl oleate as the model provided a desirable prediction of the hydroxyl value (OHV). Hence, a scaled-up process for MOG-polyol may be carried out using these optimum reaction conditions.

Author Contributions

Conceptualization, N.K., S.S.H. and L.C.A.; Methodology, N.K. and S.S.H.; Software, N.K.; Validation, N.K., S.S.H. and L.C.A.; Formal analysis, N.K. and S.S.H.; Investigation, N.K. and S.S.H.; Writing—original draft preparation, N.K.; Writing—review and editing, N.K., S.S.H., L.C.A., H.A. and D.R.A.B.; Supervision, S.S.H. and L.C.A.; Funding acquisition, S.S.H. and L.C.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministry of Higher Education Funding FRGS Grant (FRGS/1/2019/TK05/UPM/01/1) and the Malaysian Palm Oil Board (MPOB). The APC was funded by the Research Management Centre of Universiti Putra Malaysia.

Data Availability Statement

Not applicable for this study.

Acknowledgments

The authors wish to thank Universiti Putra Malaysia and the Advanced Oleochemical Technology Division, Malaysian Palm Oil Board (MPOB) for providing the research facilities.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cui, S.; Liu, Z.; Li, Y. Bio-polyols synthesized from crude glycerol and applications on polyurethane wood adhesives. Ind. Crops Prod. 2017, 108, 798–805. [Google Scholar] [CrossRef]
  2. Vanlede, K.; Kluijtmans, L.A.; Monnens, L.; Levtchenko, E. Urinary excretion of polyols and sugars in children with chronic kidney disease. Pediatr. Nephrol. 2015, 30, 1537–1540. [Google Scholar] [CrossRef] [PubMed]
  3. Chatterjee, A.; Maity, B.; Ahmed, S.A.; Seth, D. Photophysics and rotational diffusion of hydrophilic molecule in polymer and polyols. J. Phys. Chem. B 2014, 118, 12680–12691. [Google Scholar] [CrossRef] [PubMed]
  4. Saalah, S.; Abdullah, L.C.; Aung, M.M.; Salleh, M.Z.; Awang Biak, D.R.; Basri, M.; Jusoh, E.R. Waterborne polyurethane dispersions synthesized from jatropha oil. Ind. Crops Prod. 2015, 64, 194–200. [Google Scholar] [CrossRef]
  5. Prociak, A.; Malewska, E.; Kurańska, M.; Bąk, S.; Budny, P. Flexible polyurethane foams synthesized with palm oil-based bio-polyols obtained with the use of different oxirane ring opener. Ind. Crops Prod. 2018, 115, 69–77. [Google Scholar] [CrossRef]
  6. Moser, P.; Cornelio, M.L.; Nicoletti Telis, V.R. Influence of the concentration of polyols on the rheological and spectral characteristics of guar gum. LWT 2013, 53, 29–36. [Google Scholar] [CrossRef] [Green Version]
  7. Prociak, A.; Szczepkowski, L.; Ryszkowska, J.; Kurańska, M.; Auguścik, M.; Malewska, E.; Gloc, M.; Michałowski, S. Influence of chemical structure of petrochemical polyol on properties of bio-polyurethane foams. J. Polym. Environ. 2019, 27, 2360–2368. [Google Scholar] [CrossRef] [Green Version]
  8. Desroches, M.; Escouvois, M.; Auvergne, R.; Caillol, S.; Boutevin, B. From vegetable oils to polyurethanes: Synthetic routes to polyols and main industrial products. Polym. Rev. 2012, 52, 38–79. [Google Scholar] [CrossRef] [Green Version]
  9. Bakhshi, H.; Yeganeh, H.; Mehdipour-Ataei, S.; Shokrgozar, M.A.; Yari, A.; Saeedi-Eslami, S.N. Synthesis and characterization of antibacterial polyurethane coatings from quaternary ammonium salts functionalized soybean oil based polyols. Mater. Sci. Eng. C Mater. Biol. Appl. 2013, 33, 153–164. [Google Scholar] [CrossRef]
  10. Arniza, M.Z.; Hoong, S.S.; Idris, Z.; Yeong, S.K.; Hassan, H.A.; Din, A.K.; Choo, Y.M. Synthesis of transesterified palm olein-based polyol and rigid polyurethanes from this polyol. J. Am. Oil Chem. Soc. 2015, 92, 243–255. [Google Scholar] [CrossRef] [Green Version]
  11. Sabrina, S.; Rohah, A.M.; Rahmat, A.R. Novel palm oil based polyols with amine functionality synthesis via ring opening reaction of epoxidized palm oil. J. Teknol. 2018, 80, 1–7. [Google Scholar]
  12. Kong, X.; Liu, G.; Curtis, J.M. Characterization of canola oil based polyurethane wood adhesives. Int. J. Adhes. Adhes. 2011, 31, 559–564. [Google Scholar] [CrossRef]
  13. Hazmi, A.S.A.; Aung, M.M.; Abdullah, L.C.; Salleh, M.Z.; Mahmood, M.H. Producing Jatropha oil-based polyol via epoxidation and ring opening. Ind. Crops Prod. 2013, 50, 563–567. [Google Scholar] [CrossRef]
  14. Malaysian Palm Oil Board (MPOB). Pocketbook Palm Oil Statistics; Malaysian Palm Oil Board (MPOB): Kajang, Selangor, Malaysia, 2019; pp. 1–29.
  15. Abril-Milan, D.; Valdes, O.; Mirabal-Gallardo, Y.; Alexander, F.d.l.T.; Bustamante, C.; Contreras, J. Preparation of renewable bio-polyols from two species of Colliguaja for rigid polyurethane foams. Materials 2018, 11, 2244. [Google Scholar] [CrossRef] [Green Version]
  16. Pillai, P.K.S.; Li, S.; Bouzidi, L.; Narine, S.S. Metathesized palm oil & novel polyol derivatives: Structure, chemical composition and physical properties. Ind. Crops Prod. 2016, 84, 205–223. [Google Scholar] [CrossRef]
  17. Narine, S.S.; Yue, J.; Kong, X. Production of polyols from canola oil and their chemical identification and physical properties. J. Am. Oil Chem. Soc. 2006, 84, 173–179. [Google Scholar] [CrossRef]
  18. Garrison, T.F.; Kessler, M.R.; Larock, R.C. Effects of unsaturation and different ring-opening methods on the properties of vegetable oil-based polyurethane coatings. Polymer 2014, 55, 1004–1011. [Google Scholar] [CrossRef] [Green Version]
  19. Desroches, M.; Caillol, S.; Lapinte, V.; Auvergne, R.m.; Boutevin, B. Synthesis of biobased polyols by thiol-ene coupling from vegetable oils. Macromolecules 2011, 44, 2489–2500. [Google Scholar] [CrossRef] [Green Version]
  20. Hou, C.T.; Lin, J.-T. Methods for microbial screening and production of polyol oils from soybean oil through bioprocess. Biocatal. Agric. Biotechnol. 2013, 2, 1–6. [Google Scholar] [CrossRef]
  21. Tuan Ismail, T.N.M.; Ibrahim, N.A.; Mohd Noor, M.A.; Hoong, S.S.; Poo Palam, K.D.; Yeong, S.K.; Idris, Z.; Schiffman, C.M.; Sendijarevic, I.; Abd Malek, E. Oligomeric composition of polyols from fatty acid methyl ester: The effect of ring-opening reactants of epoxide groups. J. Am. Oil Chem. Soc. 2018, 95, 509–523. [Google Scholar] [CrossRef] [Green Version]
  22. Lim, S.; Teong, L.K. Recent trends, opportunities and challenges of biodiesel in Malaysia: An overview. Renew. Sust. Energ. Rev. 2010, 14, 938–954. [Google Scholar] [CrossRef]
  23. Azahar, N.F.; Gani, S.S.A.; Mohd Mokhtar, N.F. Optimization of phenolics and flavonoids extraction conditions of Curcuma Zedoaria leaves using response surface methodology. Chem. Cent. J. 2017, 11, 96. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Izadiyan, Z.; Basri, M.; Fard Masoumi, H.R.; Abedi Karjiban, R.; Salim, N.; Shameli, K. Modeling and optimization of nanoemulsion containing Sorafenib for cancer treatment by response surface methodology. Chem. Cent. J. 2017, 11, 21. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Pilkington, J.L.; Preston, C.; Gomes, R.L. Comparison of response surface methodology (RSM) and artificial neural networks (ANN) towards efficient extraction of artemisinin from Artemisia annua. Ind. Crops Prod. 2014, 58, 15–24. [Google Scholar] [CrossRef]
  26. Naresh Yadav, D.; Anand Kishore, K.; Saroj, D. A study on removal of methylene blue dye by photo catalysis integrated with nanofiltration using statistical and experimental approaches. Environ. Technol. 2020, 1–14. [Google Scholar] [CrossRef] [Green Version]
  27. Lee, P.L.; Wan Yunus, W.M.Z.; Yeong, S.K.; Abdullah, D.K.; Lim, W.H. Optimization of the epoxidation of methyl ester of palm fatty acid distillate. J. Oil Palm Res. 2009, 21, 675–682. [Google Scholar]
  28. Tajau, R.; Rohani, R.; Wan Isahak, W.N.R.; Salleh, M.Z.; Ghazali, Z. Development of new bio-based polyol ester from palm oil for potential polymeric drug carrier. Adv. Polym. Technol. 2018, 37, 3552–3560. [Google Scholar] [CrossRef] [Green Version]
  29. Kai Ling, C.; Aung, M.M.; Rayung, M.; Chuah Abdullah, L.; Lim, H.N.; Mohd Noor, I.S. Performance of ionic transport properties in vegetable oil-based polyurethane acrylate gel polymer electrolyte. ACS Omega 2019, 4, 2554–2564. [Google Scholar] [CrossRef]
  30. Rayung, M.; Aung, M.M.; Ahmad, A.; Su’ait, M.S.; Abdullah, L.C.; Ain Md Jamil, S.N. Characteristics of ionically conducting jatropha oil-based polyurethane acrylate gel electrolyte doped with potassium iodide. Mater. Chem. Phys. 2019, 222, 110–117. [Google Scholar] [CrossRef]
  31. Saithai, P.; Tanrattanakul, V. Epoxidized soybean oil-PMMA copolymer and nano-titanium dioxide. Adv. Sci. Lett. 2013, 19, 862–865. [Google Scholar] [CrossRef]
  32. Hoang, T.T.K.; Kim, I. Epoxidation and ring-opening of palm oil to produce high-functionality polyols. Aust. J. Basic Appl. Sci. 2015, 9, 89–93. [Google Scholar]
  33. Khairudin, N.; Basri, M.; Fard Masoumi, H.R.; Samson, S.; Ashari, S.E. Enhancing theBBioconversion of azelaic acid to its derivatives by response surface methodology. Molecules 2018, 23, 397. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Kamairudin, N.; Gani, S.S.; Masoumi, H.R.; Hashim, P. Optimization of natural lipstick formulation based on pitaya (Hylocereus polyrhizus) seed oil using D-optimal mixture experimental design. Molecules 2014, 19, 16672–16683. [Google Scholar] [CrossRef] [Green Version]
  35. Yang, B.; Zhang, M.; Weng, H.; Xu, Y.; Zeng, L. Optimization of ultrasound assisted extraction (UAE) of Kinsenoside compound from Anoectochilus roxburghii (Wall.) Lindl by response surface methodology (RSM). Molecules 2020, 25, 193. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Samiun, W.S.; Ashari, S.E.; Salim, N.; Ahmad, S. Optimization of processing parameters of nanoemulsion containing Aripiprazole using response surface methodology. Int. J. Nanomed. 2020, 15, 1585–1594. [Google Scholar] [CrossRef] [Green Version]
  37. Azman, N.Z.M.; Zainal, P.N.S.; Alang Ahmad, S.A. Enhancement the electrochemical conductivity of a modified reduced graphene oxide/calixarene screen-printed electrode using response surface methodology. PLoS ONE 2020, 15, e0234148. [Google Scholar] [CrossRef]
  38. Roselan, M.A.; Ashari, S.E.; Faujan, N.H.; Mohd Faudzi, S.M.; Mohamad, R. An improved nanoemulsion formulation containing Kojic Monooleate: Optimization, characterization and in vitro studies. Molecules 2020, 25, 2616. [Google Scholar] [CrossRef]
  39. Che Sulaiman, I.S.; Basri, M.; Fard Masoumi, H.R.; Chee, W.J.; Ashari, S.E.; Ismail, M. Effects of temperature, time, and solvent ratio on the extraction of phenolic compounds and the anti-radical activity of Clinacanthus nutans Lindau leaves by response surface methodology. Chem. Cent. J. 2017, 11, 54. [Google Scholar] [CrossRef]
  40. Bezerra, M.A.; Santelli, R.E.; Oliveira, E.P.; Villar, L.S.; Escaleira, L.A. Response surface methodology (RSM) as a tool for optimization in analytical chemistry. Talanta 2008, 76, 965–977. [Google Scholar] [CrossRef]
  41. Jibril, Z.I.; Ramli, A.; Jumbri, K.; Mohamad Yunus, N. Phoenix dactylifera L. seed pretreatment for oil extraction and optimization studies for biodiesel production using Ce-Zr/Al-MCM-41 catalyst. Catalysts 2020, 10, 764. [Google Scholar] [CrossRef]
  42. Zhang, L.; Cheng, Z.; Guo, X.; Jiang, X.; Liu, R. Process optimization, kinetics and equilibrium of orange G and acid orange 7 adsorptions onto chitosan/surfactant. J. Mol. Liq. 2014, 197, 353–367. [Google Scholar] [CrossRef]
  43. Yetilmezsoy, K.; Demirel, S.; Vanderbei, R.J. Response surface modeling of Pb(II) removal from aqueous solution by Pistacia vera L.: Box-Behnken experimental design. J. Hazard. Mater. 2009, 171, 551–562. [Google Scholar] [CrossRef] [PubMed]
  44. Mangili, I.; Lasagni, M.; Huang, K.; Isayev, A.I. Modeling and optimization of ultrasonic devulcanization using the response surface methodology based on central composite face-centered design. Chemometr. Intell. Lab. Syst. 2015, 144, 1–10. [Google Scholar] [CrossRef]
  45. Soi, H.S.; Kian, Y.S.; Idris, Z. Synthesis of novel oligomeric polyols from epoxidised methyl oleate and epoxidised palm olein for elastic polyurethane products. J. Oil Palm Res. 2017, 29, 401–413. [Google Scholar] [CrossRef]
  46. Abdullah, A.; Abd Gani, S.S.; Yun Hin, T.Y.; Haiyee, Z.A.; Zaidan, U.H.; Kassim, M.A.; Effendi Halmi, M.I. Lipase-catalyzed synthesis of red pitaya (Hylocereus polyrhizus) seed oil esters for cosmeceutical applications: Process optimization using response surface methodology. RSC Adv. 2019, 9, 5599–5609. [Google Scholar] [CrossRef] [Green Version]
  47. Sohrabi, M.R.; Amiri, S.; Masoumi, H.R.F.; Moghri, M. Optimization of direct yellow 12 dye removal by nanoscale zero-valent iron using response surface methodology. J. Ind. Eng. Chem. 2014, 20, 2535–2542. [Google Scholar] [CrossRef]
  48. Ang, K.P.; Lee, C.S.; Cheng, S.F.; Chuah, C.H. Synthesis of palm oil-based polyester polyol for polyurethane adhesive production. J. Appl. Polym. Sci. 2014, 131, 39967. [Google Scholar] [CrossRef]
  49. Prociak, A.; Kurańska, M.; Cabulis, U.; Ryszkowska, J.; Leszczyńska, M.; Uram, K.; Kirpluks, M. Effect of bio-polyols with different chemical structures on foaming of polyurethane systems and foam properties. Ind. Crops Prod. 2018, 120, 262–270. [Google Scholar] [CrossRef]
  50. Hoong, S.S.; Arniza, M.Z.; Mariam, N.M.D.N.S.; Armylisas, A.H.N.; Yeong, S.K. Synthesis and physicochemical properties of novel lauric acid capped estolide esters and amides made from oleic acid and their evaluations for biolubricant basestock. Ind. Crops Prod. 2019, 140. [Google Scholar] [CrossRef]
  51. American Oil Chemists’ Society (AOCS). Acid Value Te 1a-64; Hydroxyl Value Cd 13-60; Iodine Value Cd 1d-92; Oxirane Oxygen Content Cd 9-57. Official Methods and Recommended Practices of the American oil Chemists’ Society, 6th ed.; AOCS Press: Champaign, IL, USA, 2009. [Google Scholar]
Scheme 1. Reaction mechanism of epoxidation reaction for methyl oleate.
Scheme 1. Reaction mechanism of epoxidation reaction for methyl oleate.
Molecules 26 00648 sch001
Figure 1. Oxirane oxygen profile for epoxidation reaction of methyl oleate.
Figure 1. Oxirane oxygen profile for epoxidation reaction of methyl oleate.
Molecules 26 00648 g001
Figure 2. FTIR spectra of Methyl oleate (MO) and Epoxidised methyl oleate (EMO).
Figure 2. FTIR spectra of Methyl oleate (MO) and Epoxidised methyl oleate (EMO).
Molecules 26 00648 g002
Figure 3. Proton NMR Spectra of (A) Methyl oleate and (B) Epoxidised methyl oleate (EMO).
Figure 3. Proton NMR Spectra of (A) Methyl oleate and (B) Epoxidised methyl oleate (EMO).
Molecules 26 00648 g003
Scheme 2. Reaction mechanism of epoxide ring-opening reaction with glycerol.
Scheme 2. Reaction mechanism of epoxide ring-opening reaction with glycerol.
Molecules 26 00648 sch002
Figure 4. Effects of three factors on hydroxyl value. (A) Effect of mole ratio (EMO:Glycerol), (B) Effect of catalyst and (C) Effect of Temperature.
Figure 4. Effects of three factors on hydroxyl value. (A) Effect of mole ratio (EMO:Glycerol), (B) Effect of catalyst and (C) Effect of Temperature.
Molecules 26 00648 g004
Figure 5. (A) The residual plot of experimental runs from central composite rotatable design (CCRD) and (B) scatter plot of predicted hydroxyl value (mg KOH/g) versus actual hydroxyl value (mg KOH/g).
Figure 5. (A) The residual plot of experimental runs from central composite rotatable design (CCRD) and (B) scatter plot of predicted hydroxyl value (mg KOH/g) versus actual hydroxyl value (mg KOH/g).
Molecules 26 00648 g005
Figure 6. Response surface plots: (A) mole ratio (EMO:Glycerol) versus catalyst (%), (B) temperature (°C) versus Mole ratio (EMO:Glycerol) and (C) catalyst (%) versus temperature (°C).
Figure 6. Response surface plots: (A) mole ratio (EMO:Glycerol) versus catalyst (%), (B) temperature (°C) versus Mole ratio (EMO:Glycerol) and (C) catalyst (%) versus temperature (°C).
Molecules 26 00648 g006
Figure 7. Pareto plot analysis describes the effects of each variables in percentage (%).
Figure 7. Pareto plot analysis describes the effects of each variables in percentage (%).
Molecules 26 00648 g007
Figure 8. FTIR spectra for epoxidised methyl oleate (EMO) and MOG-polyol.
Figure 8. FTIR spectra for epoxidised methyl oleate (EMO) and MOG-polyol.
Molecules 26 00648 g008
Figure 9. Proton NMR Spectra of MOG-polyol.
Figure 9. Proton NMR Spectra of MOG-polyol.
Molecules 26 00648 g009
Table 1. Experimental data of actual and predicted values of hydroxyl value of MOG-polyol.
Table 1. Experimental data of actual and predicted values of hydroxyl value of MOG-polyol.
Independent VariableHydroxyl Value (mg KOH/g)
Run OrderMole Ratio (EMO:Glycerol)Catalyst (%)Temperature (°C)ActualPredicted
11:3.300.35120.00254.95252.90
21:7.500.35120.00282.40280.83
31:5.000.55150.00232.02234.87
41:7.500.35120.00281.82280.83
51:5.000.1590.00272.51274.68
61:7.500.3569.50272.27271.21
71:10.000.15150.00299.00297.66
81:10.000.1590.00303.53299.92
91:7.500.35120.00282.40280.83
101:7.500.35120.00278.35280.83
111:7.500.69120.00233.72231.84
121:7.500.35120.00281.82280.83
131:7.500.35120.00278.35280.83
141:11.700.35120.00305.04308.15
151:10.000.5590.00278.83280.85
161:7.500.01120.00289.95292.89
171:10.000.55150.00278.26275.34
181:5.000.5590.00223.80224.39
191:5.000.15150.00291.19288.41
201:7.500.35170.50275.99278.12
Table 2. ANOVA results for quadratic model developed for Hydroxyl value of MOG-polyol.
Table 2. ANOVA results for quadratic model developed for Hydroxyl value of MOG-polyol.
SourceSum of SquaresdfMean SquareF-Valuep-ValueSignificant
Model9517.6291057.51106.47<0.0001*
A3685.2513685.25371.04<0.0001*
B4499.4314499.43453.02<0.0001*
C57.64157.645.80.0367**
AB487.341487.3449.07<0.0001*
AC128112812.890.0049*
BC5.2815.280.530.4826***
A20.1610.160.0160.9014***
B2613.761613.7661.79<0.0001*
C268.43168.436.890.0254**
Residual99.32109.93
Lack of Fit80.14516.034.180.0714***
Pure Error19.1953.84
Corrected Total9616.9419
Note: *** in the table indicates insignificant difference at p-value > 0.05; ** indicates significant difference at 0.05 level; * indicates significant difference at 0.01 level.
Table 3. Regression coefficients of final reduced model.
Table 3. Regression coefficients of final reduced model.
SD3.15R20.9897
Mean274.81Adjusted R20.9804
CV%1.15Predicted R20.9296
PRESS677.12Adequate Precision37.588
Table 4. Constraints applied for optimisation.
Table 4. Constraints applied for optimisation.
Independent VariablesGoalLower LimitUpper Limit
Mole ratio (A)In range5.0010.00
Catalyst (B)In range0.150.55
Temperature (C)In range90.00150.00
Hydroxyl value Maximise223.80305.04
Table 5. Actual and predicted responses at optimum combinations.
Table 5. Actual and predicted responses at optimum combinations.
ExpIndependent VariableHydroxyl Value (mg KOH/g)RSE (%)
AB (%)C (°C)ActualPredicted
11:10.000.18120306.19301.2481.64
21:10.000.19120304.06301.2460.93
31:10.000.20120303.03301.2250.60
Note: In the table indicates A: Mole ratio (EMO:Glycerol); B: Catalyst; C: Temperature.
Table 6. Optimum condition for synthesis of MOG-polyol.
Table 6. Optimum condition for synthesis of MOG-polyol.
Independent VariableHydroxyl Value (mg KOH/g)Desirability
AB (%)C (°C)ActualPredicted
1:10.000.18120306.19301.2480.976
Note: In the table indicates A: Mole ratio (EMO:Glycerol); B: Catalyst; C: Temperature.
Table 7. Characteristics of MO, EMO and MOG-polyol.
Table 7. Characteristics of MO, EMO and MOG-polyol.
ParameterMOEMOMOG-Polyol
Iodine value (g I2 100/g)91.860.63nd
Oxirane oxygen content (%)nd5.100.05
Acid value (mg KOH/g)1.121.161.18
Hydroxyl value (mg KOH/g)16.5020.60306.19
Viscosity at 25 °C (m Pa)5.4913.59513.48
Average molecular weight, MW (Da) 315297488
Polydispersity Index1.00961.01751.1402
Note: nd in the table indicates ‘not done’.
Table 8. Summary of independent variables and their coded levels.
Table 8. Summary of independent variables and their coded levels.
Independent VariableUnitsCoded Level
−2−10+1+2
A-3.305.007.5010.0011.70
B%0.010.150.350.550.69
C°C69.5090.00120.00150.00170.50
Note: In the table indicates A: Mole ratio (EMO:Glycerol); B: Catalyst; C: Temperature.
Table 9. CCRD design matrix with 20 experiments.
Table 9. CCRD design matrix with 20 experiments.
Independent Variable
Run OrderMole Ratio (EMO:Glycerol)Catalyst (%)Temperature (°C)
11:3.300.35120.00
21:7.500.35120.00
31:5.000.55150.00
41:7.500.35120.00
51:5.000.1590.00
61:7.500.3569.50
71:10.000.15150.00
81:10.000.1590.00
91:7.500.35120.00
101:7.500.35120.00
111:7.500.69120.00
121:7.500.35120.00
131:7.500.35120.00
141:11.700.35120.00
151:10.000.5590.00
161:7.500.01120.00
171:10.000.55150.00
181:5.000.5590.00
191:5.000.15150.00
201:7.500.35170.50
Sample Availability: Samples of the bio-polyol are available from the authors.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kamairudin, N.; Hoong, S.S.; Abdullah, L.C.; Ariffin, H.; Biak, D.R.A. Optimisation of Epoxide Ring-Opening Reaction for the Synthesis of Bio-Polyol from Palm Oil Derivative Using Response Surface Methodology. Molecules 2021, 26, 648. https://doi.org/10.3390/molecules26030648

AMA Style

Kamairudin N, Hoong SS, Abdullah LC, Ariffin H, Biak DRA. Optimisation of Epoxide Ring-Opening Reaction for the Synthesis of Bio-Polyol from Palm Oil Derivative Using Response Surface Methodology. Molecules. 2021; 26(3):648. https://doi.org/10.3390/molecules26030648

Chicago/Turabian Style

Kamairudin, Norsuhaili, Seng Soi Hoong, Luqman Chuah Abdullah, Hidayah Ariffin, and Dayang Radiah Awang Biak. 2021. "Optimisation of Epoxide Ring-Opening Reaction for the Synthesis of Bio-Polyol from Palm Oil Derivative Using Response Surface Methodology" Molecules 26, no. 3: 648. https://doi.org/10.3390/molecules26030648

APA Style

Kamairudin, N., Hoong, S. S., Abdullah, L. C., Ariffin, H., & Biak, D. R. A. (2021). Optimisation of Epoxide Ring-Opening Reaction for the Synthesis of Bio-Polyol from Palm Oil Derivative Using Response Surface Methodology. Molecules, 26(3), 648. https://doi.org/10.3390/molecules26030648

Article Metrics

Back to TopTop