Synthesis of Biodiesel from Tall Oil Fatty Acids by Homogeneous and Heterogeneous Catalysis

: This study compared the yield of biodiesel produced from tall oil fatty acids (TOFA) via (i) homogeneous catalyst (sulfuric acid) and (ii) a heterogeneous catalyst (Amberlyst ® BD20, together with Ambersep BD 19 (Midcontinental Chemical Co., Olathe, KS, USA) ® using a batch reactor. The effect of operation conditions including temperature, catalyst concentration, methanol: oil ratio and reaction time on esteriﬁcation yield were investigated. Gas chromatographic data showed that the major fatty acids present in the TOFA are oleic acid (C18:1n9) and linoleic acid (C18:2n6). Homogenous catalysis yielded 96.76% biodiesel compared to 90.24% for heterogeneous catalysis. Optimized conditions for homogenous catalysis were at a catalyst concentration of 0.5 w / w %, 15:1 methanol: oil mass ratio at 55 ◦ C for 60 min. FTIR results also showed that the homogeneous catalyst yielded a more complete reaction toward biodiesel production in a shorter time (60 min) compared to the heterogeneous catalyst (4.7 h). For heterogeneous catalysis, the highest yield and the lowest acid value were achieved after a second recycling because the reactants were not fully in contact with the catalyst during the ﬁrst recycling. The catalyst did not show a reduction in catalytic activity even after the fourth recycling. However, the acid value was higher than that for ASTM standards for biodiesel.


Introduction
The continuous depletion of fossil fuels coupled with their negative impact on the environment has necessitated the search for alternate energy sources with lower environmental footprint. Biodiesel, mostly the methyl and ethyl esters of fatty acids, is currently the subject of much attention because it is a clean, renewable and biodegradable source of energy [1]. Crude tall oil (CTO), the third largest chemical by-product of the pulp and paper industry, is an ideal feedstock for biodiesel production [2]. CTO can be distilled to make TOFA, a mixture of nearly pure free fatty acids (FFA). The global market for TOFA is expected to reach USD 1 billion in 2020 [3]. As the demand for naturally derived fatty acid-based biofuel increases, there needs to be sustained research focused on newer sources of biodiesel as well as improvements in their conversion processes.
The majority of biodiesel is obtained from transterification of triglycerides from plant oils or animal fats [4][5][6][7]. Methanol or ethanol are commonly used for producing biodiesel in the presence of a base catalyst (e.g., sodium hydroxide), acid catalyst (e.g., sulfuric acid), with glycerol produced as a by-product [8]. The selection of an acid or base catalyst mainly depends on the free fatty acid (FFA) composition of oil. Additionally, FFA present in oils will affect the quality, yield, and fuel properties of biodiesel significantly [9]. For instance, when FFAs are high (typically > 1%), heterogeneous catalysis cannot be applied directly due to saponification (soap) problems. An acid-base catalyzed process can be used instead, where firstly FFAs are reduced to <0.5% by acid esterification followed by base-catalyzed conversion to biodiesel.

Synthesis of Biodiesel via Homogeneous Catalysis
Approximately 20 g of TOFA was weighed into 125 mL conical flasks. Based on the methanol: oil ratio, equivalent amounts of methanol and sulfuric acid were measured into a 250 mL conical flask. The catalyst-methanol mixture was added to the TOFA and stirred (100 rpm) at target temperature under reflux. Reaction time was calculated when the mixture attained the target temperature. At the end of the reaction period, the mixture was poured into a separating funnel and kept overnight for the separation of the water and methyl ester phases. The heavier bottom water phase was collected first by gravity, and the upper biodiesel phase (methyl ester or biodiesel) was washed with warm water several times to remove suspended impurities, unreacted methanol, catalyst, and any soap formed. The biodiesel phase was filtered through a bed of anhydrous sodium sulfate to remove all water and then stored in a refrigerator at approximately 6 • C until use.

Synthesis of Biodiesel via Heterogeneous Catalysis
Prior to esterification, TOFA was purified with Ambersep ® BD19. Heterogeneous catalysis was carried out in a water bath at a constant temperature range of 75-80 • C. Briefly, 40 g of TOFA was weighed into a conical flask and dry Amberlyst ® BD20 ion-exchange resin catalyst added. Methanol was added to the reaction mixture and then refluxed at 800 rpm for the required reaction time. The catalyst was filtered from the reaction mixture after completion. The mixture was heated in a Buchi Rotavapor R-210 under vacuum to separate unreacted methanol, residual water and any impurities from the biodiesel. The resultant biodiesel solution was poured into a separating funnel and the heavier water portion separated from the biodiesel. Finally, the biodiesel fraction was passed through a bed of anhydrous sodium sulfate (Na 2 SO 4 ) to remove residual water. The yield of biodiesel was determined gravimetrically.
TOFA samples were analyzed for elemental metal impurities using an inductively coupled plasma optical emission spectroscopy (ICP-OES). Samples were digested with concentrated HNO 3 before analysis. Acid value of the biodiesel produced was determined according to the AOCS Official Method Cd 3a-63 (Revised 1987, Reapproved 1989).
Fatty acid methyl ester (FAME) content was determined by capillary gas chromatography (GC) using a Varian CP-3800 GC equipped with a flame ionization detector (FID) and an Agilent Select Biodiesel column (100 m, 0.25 mm, 0.20 µm). GC inlet temperature was set to 250 • C with a 10:1 split, nitrogen carrier gas at 4 mL/min, and the initial oven temperature was set to 140 • C and held for 6 min, ramped up to 225 • C at 4 • C per min, and held for 15 min. Injector temperature was 270 • C and the detector temperature reached 300 • C. Identification of FAME based on comparison with retention times of FAME standards (Supleco quantitative standard FAME 37).
Fourier transform infrared spectroscopy (FTIR) was used to acquire the transmission spectra of the samples in the region between 4000-400 cm −1 . FTIR analysis was performed with a PerkinElmer Spectrum 100 FTIR. Analysis was based on the characteristic functional groups in the ester after removal of water and other impurities. Esterification of TOFA at any time under various experimental conditions was defined by the shape of the peaks in the carbonyl region of the spectrum.

Design of Experiments
In homogenous phase catalysis, a full factorial (2 4 ) experimental design-four factors each at two levels (Table 1)-was used to investigate their effect on biodiesel yield (response selected). The factors included methanol: oil molar ratio, reaction temperature, catalyst For heterogeneous catalysis, a central composite design (CCD) with temperature held constant was used. The treatment conditions were methanol: oil ratio, catalyst concentration, and reaction time. Each condition called for a low and a very high point, with intermediate points (Table 1). Each treatment condition generated was duplicated, except one central point for treatment 6, which was triplicated. In the optimization process, the influence of methanol: oil molar ratio, catalyst concentration, and reaction time on yield of TOFA biodiesel was investigated by applying response surface methodology (RSM).

Results and Discussion
The present study describes a comparison between two different catalysts used for TOFA esterification. The chosen catalysts for this study were a homogeneous phase catalyst, H 2 SO 4 , which is commonly used, and a commercial ion-exchange resin heterogeneous phase catalyst-Amberlyst BD20. FAME products from both catalysts were compared.

Elemental Analysis of TOFA
The results of the elemental analysis of TOFA are shown in Table 2. The results show that sodium and calcium account for the highest amounts in the oil sample. However, their amounts are lower than the minimum amount of 10 mg/L value for impurities in the oil sample that would cause the catalyst to lose its activity [26]. Impurities such as salts, phosphorus compounds and proteins should be removed from oil by a pretreatment process before biodiesel production because they negatively affect the production process (www.amberlyst.com) (accessed on 25 February 2021). Impurities such as proteins, phospholipids, metal ions (sodium, potassium, etc.) cause the catalyst to lose its activity. If any metal ions are present in the feed at levels greater than 10 ppm, full analysis of the oil feedstock for such impurities must be determined [26]. The oil was, however, purified using Ambersep BD 19 to forestall any deactivation of the catalyst. The purification process ensured the removal of the metal ions and other impurities that may impede the conversion process.

Physicochemical Properties of TOFA
TOFA is a mixture of acidic compounds namely resin acids, fatty acids and unsaponifiables (neutrals). Results from GC characterization of TOFA are shown in Table 3. The major fatty acids present were oleic acid (C18:1 cis-9; 57.72%) and linoleic acid (C18:2 cis-9, 12; 35.41%). Trace amounts of pentadecanoic acid (C15:1 cis-10, 1.4%), palmitoleic acid (C16:1 cis-9, 0.87%) and arachidic acid (C20:0, 2.61%) were also present. Trace amounts of minor fatty acids with a total composition of 1.98% were also recorded. These results are in good agreement with previously reported literature. For instance, Keskin et al. [28] reported TOFA composition of 52.7% oleic acid, 38.3% linoleic acid, and 6.9% linolenic acid. Variations in chemical composition of TOFA is attributed to the age, pine species, geographical location of the coniferous trees and pulping process used [29].  Table 3 also shows the fuel properties of TOFA and their methyl esters (TOME). There was a drop in kinematic viscosity of TOFA from 20.614 cSt to 5.216 cSt for biodiesel. This result is within the ASTM standard level of 1.9-6.0 cSt (ASTM D445).

FTIR Characterization of TOFA and Biodiesel
The FTIR spectra of TOFA and the biodiesel (TOME) produced via homogeneous catalysis are shown in Figure 1. The strong, sharp peaks at 2923 cm −1 and 2854 cm −1 are representative of aliphatic anti-symmetric CH 2 stretching and aliphatic symmetrical CH 2 stretching vibrations, respectively. The high intensity of these peaks is indicative of saturated oil; however, the small peak at the shoulder of the intense antisymmetric CH 2 peak at 3007 cm −1 represents unsaturation and is due to olefinic double bond [30]. Oils with a high proportion of linolenic or linoleic acyl groups show a higher frequency than those with oleic acyl groups. The absorption at 1741 cm −1 is due to the ester carbonyl (C=O), because carbonyl absorption of aliphatic esters usually appears in the 1750-1735 cm −1 region [31]. For α, β-unsaturated esters, the peaks appear in the region 1730-1715 cm −1 . The absence of a peak in this region indicates the absence of α, β-unsaturated ester, and hence there is no conjugation of cyclic compound present. This observation agrees with the GC data showing that the fatty acids present are long chain carboxylic acids.  Figure 1. FTIR spectra of TOFA and acid-catalyzed biodiesel. TOME = tall oil methyl ester (biodiesel); TOME 9.1 = biodiesel from treatment 9, replication1; TOME 1.1 = biodiesel from treatment 1, replication1; TOFA 2.1 = tall oil fatty acid.

Biodiesel Yield
The experimental matrix for the factorial design and the results of biodiesel %yield after homogeneous phase (H2SO4) catalysis of TOFA are shown in Table 4. The reaction parameters were (6:1, 15:1), (55, 75 °C), (0.5, 2% w/w), (30, 60 min) for methanol: oil ratio, temperature, catalyst concentration and time, respectively ( Table 1). The mathematical model representing biodiesel %yield as a function of the test variables in the experiment is expressed by the following equations: Statistical analysis of biodiesel %yield using JMP statistical software (JMP Pro 11) gave a p-value of 0.0001 at the 5% significance level. To test the fit of the model, the regression equation (Equation (1)) and coefficient of determination (R 2 ) were evaluated. A determination coefficient, R 2 , of 0.79 indicated that the variation in 79% TOFA biodiesel %yield obtained by homogeneous phase catalysis can be attributed to the independent variables selected (X1, X2, X3 and X4) and 21% of the total variation is not explained by the model. Other variables within the process may have affected biodiesel %yield. Statistical analysis of the experimental variables shows that methanol: oil molar ratio is the most important factor in the response. Figure 1. FTIR spectra of TOFA and acid-catalyzed biodiesel. TOME = tall oil methyl ester (biodiesel); TOME 9.1 = biodiesel from treatment 9, replication1; TOME 1.1 = biodiesel from treatment 1, replica-tion1; TOFA 2.1 = tall oil fatty acid.

Biodiesel Yield
The experimental matrix for the factorial design and the results of biodiesel %yield after homogeneous phase (H 2 SO 4 ) catalysis of TOFA are shown in Table 4. The reaction parameters were (6:1, 15:1), (55, 75 • C), (0.5, 2% w/w), (30, 60 min) for methanol: oil ratio, temperature, catalyst concentration and time, respectively ( Table 1). The mathematical model representing biodiesel %yield as a function of the test variables in the experiment is expressed by the following equations: (1) Statistical analysis of biodiesel %yield using JMP statistical software (JMP Pro 11) gave a p-value of 0.0001 at the 5% significance level. To test the fit of the model, the regression equation (Equation (1)) and coefficient of determination (R 2 ) were evaluated. A determination coefficient, R 2 , of 0.79 indicated that the variation in 79% TOFA biodiesel %yield obtained by homogeneous phase catalysis can be attributed to the independent variables selected (X 1 , X 2 , X 3 and X 4 ) and 21% of the total variation is not explained by the model. Other variables within the process may have affected biodiesel %yield. Statistical analysis of the experimental variables shows that methanol: oil molar ratio is the most important factor in the response.
For interaction variables, ratio and catalyst concentration had the most positive effect, and a slight influence from ratio and time. Temperature, catalyst concentration and time had a negative influence on the response (biodiesel yield). From Table 4, the optimum biodiesel %yield was determined with a methanol: oil ratio of 15:1, temperature of 55 • C, and catalyst concentration of 0.5% at a reaction time of 60 min. The optimum biodiesel %yield as predicted from the model was 96.76.
The main factor effects, including methanol: oil ratio, temperature, catalyst concentration, time and all two-way interaction terms, were significant at the 5% level except temperature and catalyst concentration (p-value of 0.68). Thus, it was concluded that all five factors investigated had significant effects, either as individual or as two-factor interactions. The R 2 value was 0.79, with an adjusted R 2 value of 0.73. There is therefore a high correlation between the factors. The regression explained 79% of the model.

Characterization of Homogenous Catalyzed TOFA Biodiesel Quality
FTIR spectroscopy was used to monitor the extent of conversion of TOFA to methyl ester (biodiesel) via the disappearance and appearance of certain absorption peaks. The homogeneous catalyzed conversion of TOFA to their methyl esters is a reversible reaction, and certain peaks or regions in the FTIR spectra can reveal the extent of biodiesel production. The wavenumber 3300-2500 cm −1 region gives a broad spectrum for TOFA, and this is attributed to -OH stretching absorption [31,32]. As the reaction proceeds to the formation of the ester, the broadness of this region begins to flatten and disappear completely when total conversion is achieved. The peaks observed at 3007 cm −1 were attributed to the C-H stretch of Cis-C=C bonds. The peaks at 2923 and 2854 cm −1 are related to the asymmetric and symmetric vibration of the CH 2 stretch, respectively [33].
The band at 1741 cm −1 is related to the C=O stretch of triglyceride ester linkage and the peak at 1710 cm −1 is the characteristic carbonyl stretching band of FFAs (C=OH). A single sharp peak at 1741 cm −1 (ester carbonyl C=O peak, Tome 9.1) is observed. This single peak is indicative of the higher purity of this sample compared with the spectra of the sample in Figure 1 (Tome 1.1), which has an acid C=O peak at 1709 cm −1 . The bands at 1434 and 1459 cm −1 are related to the C-H absorption of bending vibrations of CH 2 and CH 3 bands [31].
The bands of CH 2 bending vibrations can be observed at 1376 and 1236 cm −1 . The bands at 1171, 1197 cm −1 are attributed to C-O-C and O-CH3 stretching vibration of a terminal methoxy group of the ester (biodiesel) respectively [34,35]. The peaks at 1284 cm −1 and 1217 cm −1 observed for TOFA are due to the C-O stretching band (1320-1210 cm −1 ) of carboxylic acids from the interaction between C-O stretching and in-plane C-O-H bending vibrations. These are absent in the ester spectra of the biodiesel. Additionally, a peak at 1241-1244 cm −1 (in the ester band) occurred, but was absent in the carboxylic acid band of TOFA spectra. For TOFA, a peak at 1412 cm −1 due to the C-O-H bending vibration can be observed.
The band at 968 cm −1 is assigned to isolated trans-double bonds vibrational [30]. Furthermore, a peak around 930-940 cm −1 due to O-H out of plane bending vibration [31] can be observed. For TOFA (Figure 1), peaks due to acid -OH and C=O as well as C-OH bonds are more prominent. The CH 2 rocking vibration of long-chain alkanes can be clearly seen at 723 cm −1 [31].
The intensity (I) ratio of the band at 3007 cm −1 and 2854 cm −1 is a measure of the degree of unsaturation [36]. The higher the intensity ratio is, the higher the unsaturation. Acid carbonyl intensity ratios for the various treatment groups decreased, while the ester carbonyl intensity ratios increased. The intensity ratios show that treatments 8 exhibited the highest average intensity ratio of 0.9929, with 84.4% yield (due to the ester carbonyl), followed by treatment 7, with an average intensity ratio of 0.9816 and 85.3% yield (Supplementary Materials). This was followed by treatment 16 (average intensity ratio 0.9185; 93.3% yield), treatment 14 (average intensity ratio 0.9163; 94.7% yield), and treatment 3 (average intensity ratio 0.8929; 91.9% yield).
Although FTIR spectra monitoring demonstrated complete conversion of TOFA to biodiesel, statistical analysis showed that optimum biodiesel yield was achieved with treatment 3 (methanol: oil ratio of 15:1, temperature 55 • C, catalyst concentration 0.5% and reaction time of 60 min). Thus, a biodiesel yield of 96.76% was obtained under these optimum conditions, even though treatment 8 exhibited the highest intensity ratio.

Biodiesel Yield
The results of the heterogeneous catalysis of TOFA using the solid catalyst Amberlyst BD 20 are shown in Table 5. Statistical analysis of the heterogeneous catalysis of TOFA to produce biodiesel was done using the central composite design (CCD) and the response surface methodology (RSM). The results of the analysis gave an R 2 value of 0.53, with an adjusted R 2 value of 0.33. This indicates that there are other factors not in the model that affect the yield. ANOVA yielded a p-value of 0.03, which is significant at p = 0.05. The only factor that was significant at p = 0.01 was methanol: oil ratio. Catalyst concentration and time of reaction were both not significant at the p = 0.05 significance level. Likewise, all the square and interaction terms were insignificant. The regression equation generated for %Yield is given by: %Yield = 66.16 + 0.5 × Ratio.
(2) Optimum biodiesel production was obtained with an oil-to-methanol ratio of 15.5:1, a catalyst concentration with 6.4 g of Amberlyst BD20 (representing 16% of oil weight), and a time of 3.4 h. The optimum solution was therefore a Saddle Point with a %yield of 76.79 ( Figure 2). The model was maximized to predict a higher %yield of biodiesel of 90.24% Optimum biodiesel production was obtained with an oil-to-methanol ratio of 15.5:1, a catalyst concentration with 6.4 g of Amberlyst BD20 (representing 16% of oil weight), and a time of 3.4 h. The optimum solution was therefore a Saddle Point with a %yield of 76.79 (Figure 2). The model was maximized to predict a higher %yield of biodiesel of 90.24% using a methanol: oil ratio of 20.8, a catalyst concentration of 9.36 g (representing 23.4% of the oil weight), and a time of 4.7 h. For acid value, R 2 and adjusted R 2 values of 0.58 and 0.41, respectively, were obtained. ANOVA gave a p-value of 0.01, which was significant at the p = 0.05 significance level. For acid value, R 2 and adjusted R 2 values of 0.58 and 0.41, respectively, were obtained. ANOVA gave a p-value of 0.01, which was significant at the p = 0.05 significance level. Ratio was significant at a p-value of 0.03 and catalyst concentration was significant at p = 0.01 (p = 0.05 significance level). All the squared terms were insignificant. However, the interaction between catalyst concentration and time gave a p-value of 0.002 (at p = 0.05 significance level). Time was not significant, and hence the interaction term was left out of the model. The acid value regression equation is given by: The Acid Value at optimum production was predicted with a methanol: oil ratio of 20.2:1, catalyst concentration of 4.89 g (representing 12.23% of the oil weight), and time of reaction of 2.12 h. This optimal solution was a saddle point, with a predicted acid value of 27.21 (Figure 3). Minimizing the model in order to predict the smallest value of acid number gave a methanol: oil ratio of 21.4:1, and a catalyst concentration of 8.22 g (representing 20.6% of oil weight), and a time of 4.7 h. The minimization produced an acid value of 14.1 mg KOH/g oleic acid.
A graph of %yield against the acid value ( Figure 4) was linear, with an R 2 value of 0.39 and an adjusted R 2 value of 0.37. The model was significant, with a p-value of 0.0002. The linear equation is given by: %Yield = 84.01 − 0.17 × Acid Value (4) reaction of 2.12 h. This optimal solution was a saddle point, with a predicted acid value of 27.21 (Figure 3). Minimizing the model in order to predict the smallest value of acid number gave a methanol: oil ratio of 21.4:1, and a catalyst concentration of 8.22 g (representing 20.6% of oil weight), and a time of 4.7 h. The minimization produced an acid value of 14.1 mg KOH/g oleic acid.

FTIR Characterization of Heterogeneous Catalysis of TOFA
In the heterogeneous synthesis of biodiesel from TOFA ( Figure 5), FTIR spectra m itoring showed that there was a gradual shift from acid carbonyl peaks to ester carbo peaks. Ester carbonyl peaks formed became prominent attended by decrease in inten of peaks corresponding to acid carbonyl. The acid carbonyl peak finally disappeared. elimination of the acid peak is marked by the flatness of the curve beyond the 3000 c region, an indication that the -OH of the carboxylic acid is completely removed after reaction.

FTIR Characterization of Heterogeneous Catalysis of TOFA
In the heterogeneous synthesis of biodiesel from TOFA ( Figure 5), FTIR spectra monitoring showed that there was a gradual shift from acid carbonyl peaks to ester carbonyl peaks. Ester carbonyl peaks formed became prominent attended by decrease in intensity of peaks corresponding to acid carbonyl. The acid carbonyl peak finally disappeared. The elimination of the acid peak is marked by the flatness of the curve beyond the 3000 cm −1 region, an indication that the -OH of the carboxylic acid is completely removed after the reaction.
peaks. Ester carbonyl peaks formed became prominent attended by decrease in intensity of peaks corresponding to acid carbonyl. The acid carbonyl peak finally disappeared. The elimination of the acid peak is marked by the flatness of the curve beyond the 3000 cm −1 region, an indication that the -OH of the carboxylic acid is completely removed after the reaction. FTIR spectra of the various reaction conditions are shown in Figure 6. There is a gradual elimination of the acid functional groups and replacement by ester functional groups. This demonstrates the reversibility of the reaction process. The FTIR spectra in the 3375 and 2425 cm −1 regions showed a lower intensity due to complete removal of -OH acid bond resulting in flattening of the spectra curve in this region (HETBDH3). FTIR spectra monitoring the heterogeneous catalysis process showed that TOFA exhibited the highest FTIR spectra of the various reaction conditions are shown in Figure 6. There is a gradual elimination of the acid functional groups and replacement by ester functional groups. This demonstrates the reversibility of the reaction process. The FTIR spectra in the 3375 and 2425 cm −1 regions showed a lower intensity due to complete removal of -OH acid bond resulting in flattening of the spectra curve in this region (HETBDH3). FTIR spectra monitoring the heterogeneous catalysis process showed that TOFA exhibited the highest intensity of the acid carbonyl peak (44%T) followed by treatment HetBDH18 (61%T), treatment HetBDH16 (69%T), treatment HetBDH13 (76%T) and treatment HetBDH3 (88%T) in that order. This trend was reversed in the ester carbonyl peak intensities in the 1741-1742 cm −1 region. Here, treatment HetBDH3 exhibited the highest intensity of ester carbonyl peak (61%T) while treatment HetBDH18 gave the lowest peak intensity (75%T).
Sustain. Chem. 2021, 2, FOR PEER REVIEW 12 intensity of the acid carbonyl peak (44%T) followed by treatment HetBDH18 (61%T), treatment HetBDH16 (69%T), treatment HetBDH13 (76%T) and treatment HetBDH3 (88%T) in that order. This trend was reversed in the ester carbonyl peak intensities in the 1741-1742 cm −1 region. Here, treatment HetBDH3 exhibited the highest intensity of ester carbonyl peak (61%T) while treatment HetBDH18 gave the lowest peak intensity (75%T). The heterogeneously catalyzed treatments HetBDH3 and HetBDH7 are replicates for reaction conditions of Methanol:Oil mole ratio of 20:1, catalyst concentration of 20%, the time of reaction of 4 h within a temperature range of 75-80 °C (Figures 6 and 7, Supple mentary data). Treatments HetBDH3 and HetBDH7 exhibited the lowest Acid Value of 12.65 and 13.59 respectively with corresponding percent yields of 86.88% and 85.49% (average, 86.18%). However, statistical optimization of the process gave a molar ratio of The heterogeneously catalyzed treatments HetBDH3 and HetBDH7 are replicates for reaction conditions of Methanol:Oil mole ratio of 20:1, catalyst concentration of 20%, the time of reaction of 4 h within a temperature range of 75-80 • C (Figures 6 and 7, Supplementary Materials). Treatments HetBDH3 and HetBDH7 exhibited the lowest Acid Value of 12.65 and 13.59 respectively with corresponding percent yields of 86.88% and 85.49% (average, 86.18%). However, statistical optimization of the process gave a molar ratio of Methanol:Oil as 20.8:1, a catalyst concentration of 23.4% (of the weight of oil), and a reaction time of 4.7 h to give a yield of 90.24%. Importantly, the acid value did show that HetBDH3 and HetBDH7 had the lowest values consistent with the carbonyl peaks from the FTIR spectra. FTIR spectra monitoring of the reaction showed that treatments HetBDH3 and HetBDH7 were nearest to the optimum conditions.  However, the optimum conditions for the heterogeneous catalyst, Amberlyst BD 20, to produce a 90.24% yield was given by Methanol:Oil Molar ratio of 20.8:1, Catalyst concentration of 23.4% and Time of Reaction of 4.7 h. The acid value attained in this process was too high for use as a biodiesel. There is the need to lower the acid value of the biodiesel produced.

Recycling of Catalyst
The effect of the recycling of the catalyst on the biodiesel production process is presented in Table 6 and Figures 8 and 9. The catalyst recycling step was evaluated because of its impact on the overall cost of the production process. The %yield showed an increase after the first reaction, an indication that the catalyst However, the optimum conditions for the heterogeneous catalyst, Amberlyst BD 20, to produce a 90.24% yield was given by Methanol:Oil Molar ratio of 20.8:1, Catalyst concentration of 23.4% and Time of Reaction of 4.7 h. The acid value attained in this process was too high for use as a biodiesel. There is the need to lower the acid value of the biodiesel produced.

Recycling of Catalyst
The effect of the recycling of the catalyst on the biodiesel production process is presented in Table 6 and Figures 8 and 9. The catalyst recycling step was evaluated because of its impact on the overall cost of the production process. catalyst in the wavelength region of 1200-1100 cm −1 (Figure 9) shows that the first cycle had the highest acid carbonyl peak, and that the recycling reduces that peak while increasing the ester carbonyl peak. The highest catalytic activity is shown for the second cycle, and while it gradually decreases towards the fourth cycle, it remains still potent in its activity.

% Yield and Acid Value of recycled Treatment HetBD_H3
Acid Value % Yield  The %yield showed an increase after the first reaction, an indication that the catalyst was more potent after the first cycle. Similarly, the acid value gave the highest value of 14.53 after the first cycle, and the value continued to fall to a minimum of 9.72 after the fourth cycle ( Figure 8). This phenomenon whereby the catalyst increased the process yield to the maximum before becoming stable may be due to swelling, as reported by Lopez et al. [37]. Park et al. [26] also reported that Amberlyst BD 20 increases its catalytic activity from the second cycle and ascribed it to time required by the catalyst to contact the reactants. Consequently, the catalytic activity of the first cycle is retarded. The results also demonstrated that the catalyst could be recycled even to the fourth cycle without decrease in activity of the catalyst (Figure 8). This is consistent with the work of Telebian-Kiakalaieh et al. [38], Park et al. [26]. The overlaid FTIR spectra of the biodiesel from the recycled catalyst in the wavelength region of 1200-1100 cm −1 (Figure 9) shows that the first cycle had the highest acid carbonyl peak, and that the recycling reduces that peak while increasing the ester carbonyl peak. The highest catalytic activity is shown for the second cycle, and while it gradually decreases towards the fourth cycle, it remains still potent in its activity.

Conclusions
FTIR spectroscopy was used to monitor the progress of biodiesel synthesis from tall oil fatty acids using homogeneous catalyst (sulfuric acid), and heterogeneous catalyst (Amberlyst BD 20, from Dow Chemicals). The optimum reaction conditions for the homogeneous catalysis was a methanol: oil molar ratio of 15:1 at a Temperature of 55 • C, catalyst concentration of 0.5%, and a reaction time of 60 min. The yield obtained was 96.76%. The FTIR data and spectra agreed with the optimum treatment conditions showing a total conversion of the fatty acid into biodiesel. The acid value was reduced from 113.18 ± 0.01 mg KOH/g oleic for the TOFA to 4.24 ± 0.12 mg KOH/g oleic acid for the biodiesel. This acid value is too high above the ASTM standard value of 0.8 mg KOH/g oleic acid. Kinematic viscosity was reduced from a value of 20.614 ± 0.059 mm 2 /s for the TOFA to 5.216 ± 0.016 mm 2 /s for the biodiesel which is within the ASTM standard value of 1.9-6.0 mm 2 /s. For the heterogeneous catalysis process, the optimum conditions used for the biodiesel production were a methanol: oil molar ratio of 20.8:1, a catalyst concentration of 23.4% of the oil, and a reaction time of 4.7 h at a constant temperature range of 75-80 • C. The predicted optimum yield obtained was 90.24%. The acid value for the heterogeneous catalysis was higher, at a value of 14.1 mg KOH/g of oleic acid. There is a need to further reduce the acid value in order to meet the ASTM standards for biodiesel. The major fatty acids present in the TOFA are oleic and linoleic acids. The biodiesel percent yield from the homogeneous catalysis was higher than that of the heterogeneous catalysis. However, the heterogeneous catalyst could be recycled about four times without losing its catalytic activity. The FTIR spectra agree with the yield data, showing a better conversion of biodiesel for the homogeneous catalysis than the heterogeneous catalysis.