Abstract

Magnetic materials based on Keggin heteropolyacids immobilized on mesoporous silica-coated magnetite particles with a core-shell structure were synthesized. The activity of the catalyst was studied in the esterification reaction of levulinic acid with n-butanol, and its kinetics was studied by the systematic variation of several reaction parameters, such as stirring speed, catalyst loading, molar ratio of reactants, and temperature. It was also seen that the reaction was free from any external mass transfer as well as intraparticle diffusion limitations and was intrinsically kinetically controlled. A second-order kinetic equation was found to be consistent with the experimental data. Also, an experimental activation energy of 17 kcal/mol was found. A solvent-free condition for this reaction has also added the green chemistry perception to the reaction. In addition, the resulting catalyst can be used repeatedly without significant decrease in activity.

1. Introduction

In recent years, great interest has been aroused in biomass waste and the carbohydrates derived from it, as they are abundant and relatively inexpensive and are a promising renewable source of carbon. The biomass is a useful chemical platform for synthesizing a broad range of products such as fuels, additives, and fine chemicals in general, thus representing an alternative to petroleum-derived products [1, 2]. Biomass wastes include any organic material other than the raw material for which they were planted, e.g., waste from leaves and bagasse from the sugar cane during sugar production or paper or fruit peels and seeds in the food industry [3]. Thus, biorefineries are born, which are becoming increasingly important in the development of new sustainable technologies and are defined as an industry that, through different technological processes (mechanical, physical, and biochemical), separates renewable raw materials into their main components such as triglycerides, carbohydrates, and proteins to produce energy and value-added chemicals analogously to how oil refineries process crude oil to produce specialty chemicals [4, 5]. In 2004, researchers from the NREL (National Renewable Energy Laboratory) and PNNL (Pacific Northwest National Laboratory) conducted an extensive study to identify the most valuable building blocks based on sugars from lignocellulosic biomass [6]. Out of 300 candidates initially selected, the number was reduced to twelve, one of these top twelve promising building blocks being 4-oxopentanoic acid, known commercially as levulinic acid (LEA). In 2010, it was reconfirmed in the top 10 list, after an updated assessment of possible building blocks in which the advances of technology that led to the inclusion of a compound in the first listing were evaluated [7] (Scheme 1 shows some useful compounds derived from LEA). Hence, various investigations on how to synthesize LEA and different by-products from biomass derivatives are being carried out. Among them, the alkyl levulinates (ALE), especially n-butyl levulinate (BLE), are of great importance, because they currently have extensive applications, either as a fragrance [8] or as an additive in gasoline and biodiesel, among others [9].

Most of the current industrial processes for the production of LEA and ALE use sulfuric and phosphoric acids as catalysts in each of the steps [10], which involves disadvantages harmful to health and the environment due to its toxicity. For these reasons, this study focuses on providing a heterogeneous catalyst that facilitates the process and the performance of the industrial scale reaction. In this sense, the Fischer esterification of LEA with n-butanol (BOH) to synthesize BLE was chosen as the study reaction using a magnetic solid with tungstophosphoric acid (TPA) as active phase as catalyst, denominated Fe3O4@SiO2@SiO2/TPA400N2. To date, there are few reports on the kinetics of this reaction. In general, Bart et al. [11] gave kinetic information about the esterification of LEA with BOH with sulfuric acid as a homogeneous catalyst. Subsequently, Nandiwale and Bokade [12] studied the kinetics of the same reaction using modified H-ZSM-5 catalysts. Now, we intend to extend these kinetic studies and report data on diffusion effects, kinetic equation, and activation energy.

On the other hand, in the study of magnetic materials as catalysts, the core-shell structure has attracted increasing attention due to two main reasons: (i) the magnetic center allows controlling its dispersion in a reversible way by applying a magnetic field. Namely, without the presence of an external magnetic field, the nanoparticles are dispersed in the solution, redistributing their surface area facilitating the contact with the substrate molecules. After completing the reactions, the nanoparticles can be easily isolated from the system by approaching a magnet [13]. (ii) The coating serves as an anchor point for different materials and functional groups and thus, for example, immobilizes homogeneous catalysts [14].

The magnetic nanoparticles (MNPs) were synthesized by the coprecipitation method, by means of iron salts (Fe2+/Fe3+) in ammonia solution under a nitrogen atmosphere. This method is widely used because it is simple and economical and has a high reaction yield and a narrow size distribution [15]. One of the disadvantages of the magnetite particles is that they are susceptible to oxidation by environmental factors or by acid reaction conditions. This makes it necessary to encapsulate them with an organic/inorganic polymer shell or stabilize them with a surfactant, further conferring an active surface to anchor the active phase of the catalyst [16, 17]. The most commonly used coating material is silicon oxide because it is chemically inert, nontoxic, and relatively easy to synthesize by the sol-gel method [18]. The TPA was included during the coating of the magnetic particles with the silicon oxide. There are a few studies on this type of catalyst, in which the heteropoly acid is generally added by impregnation and has never been evaluated in this esterification reaction [1922].

The characterization techniques used were IR/FT, XRD, BET, SEM, VSM, and potentiometric titration. The best reaction conditions, in terms of reactant molar ratio, catalyst loading, temperature, and reaction time, were established, as was the reuse of the catalyst. A possible catalytic mechanism for the esterification reaction is also proposed.

2. Experimental

2.1. General

All chemical reagents and solvents were of commercial analytical grade and were used without further purification. The progress and conversion of the reactions were quantified by GC/FID in a Shimadzu chromatograph model 2014, using a FID detector and a capillary column (SPB-1, length 30 m, I.D. 32 mm, and film thickness 1.00 μm). The carrier gas was high-purity nitrogen (≧99.999%), and the optimum conditions were as follows: splitless injection mode was used (1 μL), the injector temperature was set at 320°C, the initial column temperature was 150°C (held for 2 min), then ramped at 20°C/min to 200°C, and the FID temperature was 300°C. Conversions were determined by the relative peak areas of the substrates and products using a normalization method. The retention times of LEA and BLE were 3.42 and 6.15 min, respectively. For more details on the quantification method, see Supplementary Material (available here).

2.2. Catalyst Preparation
2.2.1. Synthesis of Fe3O4 Nanoparticles

Fe3O4 MNPs were synthesized through the chemical coprecipitation method of iron salts in an ammonia solution. The salts FeCl2·4H2O (1.0) and FeCl3·6H2O (2.7 g) (2 : 1 molar ratio) were dissolved in water (25 mL), mixed under vigorous stirring (1000 rpm) with nitrogen bubbling, and heated to 353 K. When they reached the indicated temperature, 12 mL of concentrated NH4OH solution (25% ) was added dropwise, controlling the pH between 11 and 13. Immediately, the solution turned black and was stirred for another 1 h. Then, the Fe3O4 nanoparticles were precipitated using a magnet. Finally, MNPs were washed repeatedly with water and ethanol to eliminate excess ammonia. Nitrogen was used in the whole process to prevent particle oxidation.

2.2.2. Synthesis of Fe3O4@SiO2

Silica coatings on the Fe3O4 MNPs were synthesized by the Stölber method with some modifications and basic hydrolysis, utilizing tetraethyl orthosilicate (TEOS) as SiO2 source. One g of Fe3O4 MNPs was sonicated for 30 min in 70 mL of ethanol and 24 mL of ammonia solution (25% ), at the molar ratio of 1 : 6 with respect to TEOS, heating to 318 K and with constant stirring. When the temperature was reached, TEOS was added dropwise to the mixture and reacted for 4 h. Then, the ethanol was evaporated at room temperature to dryness while stirring. The particles were separated using a magnet and washed four times with water and ethanol until neutral pH. The MNPs showed no detectable dissolution in a hydrochloric/nitric acid mixture for 24 h.

2.2.3. Synthesis of TPA/SiO2 Composite Coatings over Fe2O4@SiO2

A second coating was made with SiO2, in which the TPA was simultaneously added, so that it was immobilized in this layer. In this case, the coating was realized by acid hydrolysis of TEOS catalyzed by the same TPA. In a typical experiment, 1 g of the Fe3O4@SiO2 MNPs was redispersed in 70 mL of ethanol together with 1 g of TPA. The obtained mixture was heated to 318 K; when this temperature was reached, 4 mL of TEOS was slowly added and the heating was stopped. The sample was left under stirring at room temperature until dry. The product was washed with water and ethanol until neutral pH. Finally, the obtained solid was calcined at 673 K under nitrogen atmosphere for 2 h. The catalyst was called Fe3O4@SiO2@SiO2/TPA400N2. The amount of TPA was fixed in order to obtain a concentration of 40% TPA () in the catalyst.

The TPA content in the catalyst was estimated as the difference between the W amount contained in the heteropoly acid solution originally used for the impregnation and the amount of W in the water solutions obtained after the solid washing. The W content was determined by atomic absorption spectrometry using a Varian AA Model 240 spectrophotometer. The calibration curve method was used with standards prepared in the laboratory. The analyses were carried out at a wavelength of 254.9 nm, band width 0.3 nm, lamp current 15 mA, phototube amplification 800 V, burner height 4 mm, and acetylene–nitrous oxide flame 11 : 14.

2.3. Catalyst Characterization

SEM images were obtained on Philips SEM 505 equipment, using a voltage of 25 eV and increasing from 2500x to 10000x. The samples were premetallized with gold. The diameter of the material particles was calculated by direct observation of SEM images using an EDAX Apex 2 analyzer system, Apollo X Model, which has software for obtaining and analyzing digital images.

XRD diagrams were obtained in order to determine the presence of crystalline phases associated with Fe3O4 and TPA. They were recorded in PANalytical X’Pert Pro equipment, using Cu kα radiation and Ni filter. The scanning angle was performed in the range 10–90° of 2θ using a scan rate of 1 s and a step size of 0.04° s-1.

Textural characteristics of materials were obtained from N2 adsorption–desorption isotherms at liquid-nitrogen temperature, using Micromeritics ASAP 2020 equipment. Previously, the solids were degassed at 373 K for 2 h.

Magnetization measurements were obtained on a vibrating magnetometer VersaLab. For this purpose, 0.01 g of sample was weighed and a scan of -10000 to +10000 Oe was performed at 300 K.

The chemical groups present during the synthesis of the materials were analyzed on a Nicolet iS50 FTIR spectrophotometer coupled with attenuated total reflectance (ATR). For the analysis, 0.1 g of each solid was weighed and placed in the ZnSe cell at room temperature, with a scan of 32 scans and a resolution of 6 cm-1. The collected spectra were manipulated with the program OMNIC 5.1®.

The acid strength was determined by potentiometric titration with n-butylamine (for more information about the method, see Supplementary Materials).

2.4. Catalytic Studies

The esterification reaction of LEA with BOH was selected as a model reaction to evaluate the catalytic activity of Fe3O4@SiO2@SiO2/TPA400N2 material. For all the optimization tests, the same procedure was followed and the same equipment was used. A reaction tube with a water reflux condenser was used as a reactor. When the chosen temperature was reached, the predetermined amounts of reagents in the tube were added and agitated at 500 rpm. An experimental design with three variables was proposed; these include the reaction temperature (323, 333, 343, and 353 K), acid : alcohol molar ratio (1 : 2, 1 : 4, 1 : 6, and 1 : 8 mmol/mmol), and catalyst amounts (25, 50, 100, and 150 mg). Samples were taken at regular intervals by momentarily reducing the stirring rate and bringing a magnet closer to precipitate the catalyst. In all reactions, the amount of LEA (1 mmol) remained constant, BOH was in excess with respect to LEA, and the reaction conversion was calculated based on the limiting reagent, LEA. Each parameter was analyzed in triplicate with standard deviations less than ±2%.

2.5. Catalyst Reuse

To evaluate the reuse of the catalyst, four consecutive experiments were carried out under the same reaction conditions. After each reaction cycle, the catalyst was separated from the reaction system with the aid of a magnet, washed with toluene (3 × 2 mL), dried under vacuum, and reused by adjusting the amounts of reactants to the amount of catalyst recovered.

3. Results and Discussions

3.1. Catalyst Characterization

The morphological characteristics of each catalytic material obtained in the different stages of catalyst synthesis were analyzed by SEM. Figure 1(a) shows the SEM image of the MNPs prepared through the chemical mixture of ferric and ferrous salts, precipitating as Fe3O4 by the addition of ammonia, which have a uniform shape and a size around 10 nm. The chemical equation of Fe3O4 synthesis is shown below:

Then, the MNPs were sonicated in an alkaline ethanol-ammonia solution and the magnetic nucleus was coated with silica shells by hydrolysis of TEOS (Figure 1(b)), increasing the size of the nanoparticles to an average diameter of 300 nm. Finally, by means of the second coating with TEOS, TPA was immobilized (Figure 1(c)), increasing the average diameter to 620 nm. In all cases, the particles were spherical and of uniform size. However, some aggregates could be observed after the second SiO2 coating and remained after calcination at 673 K (Figure 1(d)).

From the adsorption-desorption nitrogen isotherms at 77 K, the surface areas were estimated applying the Brunauer-Emmett-Teller (BET) equation. Table 1 shows the results obtained from the measurement of the surface area and pore size of the materials. The values determined for the Fe3O4 magnetic particles show a low corresponding to 6 m2/g (Table 1, entry 1). This area increases to 20 m2/g with the first SiO2 coating (Table 1, entry 2) and increases again to 65 m2/g when the second coating is made with SiO2 and with the addition of TPA (Table 1, entry 3). The contribution to the total area, due to the presence of micropores, is around 40% in the sample Fe3O4@SiO2 and decreases to 15% for Fe3O4@SiO2@SiO2/TPASC. When subjecting the catalyst to thermal treatment at 673 K under a nitrogen atmosphere, decreases to 45 m2/g, also decreasing the contribution of micropores (Table 1, entry 4), possibly due to the collapse of the smaller pores, evidenced by an increase in the pore diameter from 7 to 9 nm.

By analyzing the shape and hysteresis loop of the adsorption-desorption isotherms shown in Figure 2, it can be seen that the isotherm for the Fe3O4 sample is typically type II, characteristic of nonporous, slightly porous, or macroporous materials. Therefore, a pore size for this sample was not determined (Table 1, entry 1). The samples coated with silica and TPA (Table 1, entries 2 to 4) present isotherms type IV, characteristic of mesoporous solids, which is in agreement with the obtained pore diameter. However, the different cycles of hysteresis suggest that the incorporation of silica and TPA into Fe3O4 particles modifies the textural properties, which is why it significantly increases its surface area and pore size (Table 1). Thus, the N2-sorption isotherm of Fe3O4 and Fe3O4@SiO2 magnetic particles exhibits a barely perceptible hysteresis, which increases when the second layer of silica and TPA is added. Hysteresis loops in both uncalcined and calcined catalysts are type H2.

The FTIR-ATR spectra are shown in Figure 3. In the Fe3O4 spectrum, the main absorption band is centered at around 540 cm-1, corresponding to stretching vibration of Fe-O. In the case of sample Fe3O4@SiO2, the characteristic absorption band for the silica network is assigned as follows. The high-intensity broad band at 1080 cm-1 is due to Si-O-Si asymmetric stretching bonds in the SiO4 tetrahedron associated with the movement of oxygen in the Si-O-Si antisymmetric stretch. The signal at 793 cm-1 is assigned to the Si-O-Si symmetric stretch. The absorption band at 960 cm-1 corresponds to stretches of the Si-OH bond. And the absorption band at 564 cm-1 is an indication of the presence of Si-O-Fe; this band overlaps the Fe-O band of the magnetite [18]. The structure of the Keggin ions of [PW12O40]3- is well known, and their bands have been reported in previous works [23, 24] (they are also detailed in Supplementary Material). In the sample Fe3O4@SiO2@SiO2/TPA400N2, they are overlapped by the wide bands originated from the silica network [20, 22].

The XRD of the Fe3O4@SiO2 and Fe3O4@SiO2@SiO2/TPA400N2 samples are shown in Figure 4(a). In them, the diffraction peaks located at , 30.2°, 35.6°, 43.3°, 53.6°, and 62.6° are assigned to planes (111), (220), (311), (400), (422), and (440), respectively. They could be assigned to magnetite (Fe3O4) or maghemite (γ-Fe2O3) phases, because the diffraction pattern of the cubic form of γ-Fe2O3 is identical to that of Fe3O4 with some line shift towards higher angles [25]. This makes it difficult to define what iron oxide is synthesized by analyzing the XRD. However, the black color in our samples was maintained, indirectly indicating the preferential formation of magnetite, because magnetite is usually black and pure maghemite is usually brown [26]. In addition, it is evident that the coating with SiO2 and calcination in an inert atmosphere protect the magnetite from oxidizing, since in the sample Fe3O4@SiO2@SiO2/TPA400N2, the diffraction peaks are exactly the same, although a little less defined and wider because of a greater amount of SiO2, indicating that the silica is in an amorphous form and possibly its small and wide peaks are dominated by the sharp peaks of the magnetite. Even so, a broad peak at 10° and a shoulder between 25° and 30° characteristic of the amorphous silica are observed. Finally, the characteristic diffraction peaks of H3PW12O40·23H2O and H3PW12O40·6H2O (shown in Figure 4(b)) are not observed in Figure 4(a) due to a large dispersion of the active phase on the support surface [27, 28]. The small peaks at 30.9° and 38.8° (marked with an asterisk in Figure 4(a)) could be assigned to the presence of a small amount of H3PW12O40 crystals.

The magnetization curves for the Fe3O4, Fe3O4@SiO2, and Fe3O4@SiO2@SiO2/TPA400N2 materials are shown in Figure 5. The applied magnetic field for these systems confirms the ferromagnetic response of these materials. It can be observed that the Fe3O4 and Fe3O4@SiO2 curves pass through the origin (Figure 5(a), blue and red lines), giving rise to solids that do not present coercivity or remanence, whish evidences their superparamagnetic behavior [29, 30]. The magnetization saturation () values for these two materials were 87 and 52 emu/g, respectively. The decrease in the value for Fe3O4@SiO2 can be attributed to the deposition of nonmagnetic SiO2 on the Fe3O4 surface, with the consequent lower amount of magnetite in the same sample weight [28]. However, in the Fe3O4@SiO2@SiO2/TPA400N2 solid (Figure 5(a), green line), the value drops considerably to 15 emu/g, evidencing that the SiO2 layer is much thicker, with the consequent lower amount of magnetite per unit area. In addition, a slight hysteresis loop can be observed, indicating that the catalyst has a small remanence, probably due to a significant increase in size of the catalyst particles, to a lower crystallinity, or to changes in the surface structure, as evidenced in the XRD analysis. The content of impurities or deficient crystallization of SiO2/TPA can affect the magnetic behavior of the powder synthesized, reducing the maximum magnetization and slightly increasing the coercive force [26, 31]. Even so, its magnetization is sufficient to easily separate the solid from the reaction medium with a conventional magnet. Finally, the magnetic behavior of the uncalcined and calcined materials, after two reaction cycles, was analyzed and compared to see whether the calcination temperature and the reaction medium had any effect on their magnetic properties. The resulting curves are shown in Figure 5(b). It can be seen that the value remains at 15 emu/g in the 4 samples and maintains the same hysteresis cycle, indicating that the SiO2 layers protect the magnetite from abrupt changes of temperature and reaction medium.

The acidity of the solids was determined by potentiometric titration with n-butylamine; the curves obtained for the Fe2O4@SiO2 support, the Fe3O4@SiO2@SiO2/TPA400N2 catalyst, and a catalyst sample after the first reaction cycle, Fe3O4@SiO2@SiO2/TPA400N2C1, are shown in Figure 6. The Fe2O4@SiO2 support showed very strong sites, with an Ei value of 110 mV. This value increased significantly to 410 mV after the TPA was included in a second SiO2 layer. However, as expected, by including TPA in Fe2O4@SiO2, the strength of the acid sites decreased compared to the value recorded for the bulk TPA, which has an Ei value of 620 mV. This is an effect that depends on the degree of interaction of TPA with the support, in this case indicating a strong interaction between the protons of H3PW12O40·6H2O and the oxygen of the Si-OH groups. Furthermore, a good dispersion of TPA on the support is demonstrated by the absence of leaching during the reaction [32]. It can be seen that the potentiometric titration curve of Fe3O4@SiO2@SiO2/TPA400N2C1 gives an Ei value of 395 mV, which is slightly smaller compared to the unused catalyst Ei value (410 mV). The form and the potentials of the two curves are very similar, demonstrating that the amount leached during the reaction is very small, which is in agreement with the results obtained in the reuse experiment (4 cycles), exhibiting an excellent activity in esterification reactions.

3.2. Kinetic Study
3.2.1. External Mass Transfer Resistance

One way of experimentally determining the external resistance to mass transfer is the stirring test, since, in the case of gas-solid and liquid-solid systems, when the conversion does not depend on the stirring rate, it can be said that the system is free of interfacial limitations [33]. The test consists of carrying out different reactions by keeping the molar ratio of reactants, catalyst weight, and temperature constant and varying the stirring rate for a given time. In the range of stirring rates in which the conversion is kept constant, it can be said that there are no diffusion effects. In this case, the experimental conditions were 313 K, an initial acid/alcohol molar ratio of 1 : 2, catalyst of 100 mg, and stirring rates of 100, 250, 500, 750, and 1000 rpm for 3 h. It was found that in the interval between 100 and 1000 rpm, there was no effect of the stirring speed on the conversion of LEA, since in all cases the conversion was . Thus, all the experiments were performed at a stirring rate of 500 rpm to ensure the absence of external resistance.

3.2.2. Internal Mass Transfer Resistance

The Weisz-Prater (WP) criterion () [34] was used to evaluate the influence of intraparticle diffusion resistance at each temperature. Such a calculation makes it possible to establish the absence of any pore diffusion limitation or to warn of its presence. The criterion states that an value below 0.3 for a reaction order in the substrate of 2 or less is indicative of the absence of significant diffusion limitations. To use the WP criterion successfully, the effective diffusivity coefficient () of a reagent in a given catalytic system must be determined accurately. In this case, it would involve determining the effective diffusivity of LEA (A) in BOH (B) (calculated from the Wilke-Chang criterion, porosity, and tortuosity), reaction rate, particle radius, and concentration of the reactants on the outer surface of the particle. For the latter, it was assumed that when the external mass transfer limitations can be discarded, the concentration of substances on the catalyst surface is taken equal to the mass concentration of the substances (). Thus, the mass concentration of A at the beginning of the reaction is known for all experiments. Since the values determined for each temperature used during the catalytic test in all cases were much lower than 0.3, it can be established that there is no resistance to intraparticle diffusion, which is important, since, in a kinetic study, it is necessary to generate data that are free of mass and energy transport limitations (the equations and calculations performed are shown in the Supplementary Material; similarly, the determined data are provided in Table S1).

3.3. Optimization of the Reaction Conditions
3.3.1. Influence of the Catalyst Amount

The catalyst amount was initially studied in the range from 25 to 150 mg under the same reaction conditions, namely, an acid/alcohol molar ratio equal to 1 : 2 and a temperature of 353 K for 5 h. The reaction profiles are shown in Figure 7. Since there are no external diffusion limitations, the reaction rate is proportional to the amount of catalyst, because there are a greater number of acid sites. In the absence of catalyst, a conversion of 25% in 5 h was only obtained, as the acidity of LEA catalyzed the reaction, which was further favored by an increase in temperature, since at 323 K it gave a conversion of 15% in 24 h. However, this is a lower value compared to the conversion obtained in the presence of the catalyst Fe3O4@SiO2@SiO2/TPA400N2, in which case with 25 mg of catalyst, a conversion of 82% was reached and with amounts of 50, 100, and 150 mg the conversion were 95%, 98%, and 100%, respectively. The values were significantly close, so a load of 100 mg per mole of LEA was chosen as the optimum catalyst amount.

3.3.2. Influence of Temperature

The influence of the reaction temperature on the conversion of BLE was studied in a range between 323 and 353 K (Figure 8). Not only did a significant increase in the conversion of the reaction was observed when the temperature was raised, but the reaction rate also increased, achieving a conversion of 100% in 5 h at 353 K compared to the 42% conversion in 7 h that was reached at 323 K. This suggests a kinetically controlled mechanism. An increase in the reaction temperature leads to drastic increases in collisions between reactant species and the catalyst. In consequence, the reaction rate increases, which favors the formation of BLE. In addition, the esterification reaction is an equilibrium reaction that generates water as one of the products. As the temperature increases, the water is removed from the system, displacing the equilibrium towards ester formation. It can be seen in Figure 8 that at 323 K, the reaction is slow and the reaction rate significantly increases at 343 and 353 K.

3.3.3. Influence of the Acid/Alcohol Molar Ratio

Another way of shifting the equilibrium towards the formation of BLE is by using an excess of alcohol, but in this case (Figure 9), an opposite effect was obtained, which agrees with the one published by Patil et al. [35], who found that a low molar ratio of alcohol to LEA, in the range of 2 to 5, favors the reaction. The said behavior is explained because an overly diluted system could involve limitations of the mass transfer of reactant to the catalyst and stimulate the intermolecular dehydration of the alcohol as a by-product, which occurs at low temperatures and a higher reactant molar ratio. In our case, the best molar ratio was 1 : 2, which favors the reaction rate and selectivity and allows it to be run in solvent-free conditions, with the additional advantage of improving the process of eco-compatibility.

3.4. Kinetic Model

With the experimental data presented previously, a kinetic model was constructed according to the reaction under study. The reaction can be described by either a pseudohomogeneous (pH), Langmuir–Hinshelwood (LH), or Eley–Rideal (ER) model. An overall rate of expression assuming that the slowest step is the surface reaction and not the adsorption and desorption steps is expressed as where for the pH model, for the ER model, and for the LH model. This equation has been described by various authors in the esterification of organic acids [3639]. Assuming that for the esterification of LEA, there is a weak sorption of reaction components on the catalyst surface (equilibrium constants of adsorption are low) so that the term ∑kiCi in the denominator is almost negligible [40]; the expression of the pH model is reduced to where  = catalyst loading.

This expression is commonly used in the esterification of LEA [12, 41]. In this sense, equation (3) can be rewritten in terms of fractional conversion as follows: where , and and are the concentrations of LEA at any reaction time and initial time, respectively. Solving the integral form of (4), we have

Plots of against reaction time at different temperatures are shown in Figure 10. In all cases, the curves are linear indicating that the esterification of LEA with BOH is of the second order.

To improve the adjustment of the kinetic model, the adsorption steps of the reactants or the products should be considered. It is well established that the esterification reactions can be ascribed to either a Langmuir–Hinshelwood or Eley–Rideal pathway [42]. To avoid the contribution of products and simplify the overall reaction scheme, the initial reaction rates were only considered. The initial reaction rate was calculated using a differential method at time zero and the following equation:

The plots of the initial reaction rate as a function of or (not shown) describe a linear dependence with increasing LEA, but decreases linearly with increasing BOH. This behavior is explained due to saturation of the catalyst surface with the alcohol, thereby stopping the reaction by blocking the adsorption of LEA. The reaction occurs between the two reagents, levulinic acid (A) and n-butanol (B), in only one active acid site (S), but as was previously mentioned, the alcohol has an inhibiting effect. Considering an LH expression and assuming that the adsorption of BOH is stronger than the adsorption of LEA, then , equation (2) could be simplified to

The linearization of equation (7) allows a plot of against (Figure 11) yielding a and an . The fit was found to be good at different temperatures, and the reaction rate and adsorption constants thus obtained are summarized in Table 2.

An apparent activation energy of 17 kcal/mol was determined; besides, the constants of adsorption decrease with the temperatures, which could indicate that the enthalpies and entropies of adsorption are negative, according to their thermodynamic consistency. Considering the saturation of the catalyst surface with the alcohol and a minor adsorption of LEA on one active site and based on this LH model, we propose a plausible reaction mechanism for the esterification of LEA with BOH, which is outlined in Figure 12. The mechanism involves adsorption of BOH on a Brønsted acid site (B), and subsequently, LEA is activated and attacked by BOH through an ethoxide mechanism. The stronger adsorption of alcohol in this figure is represented as two molecules of alcohol by the active site, which inhibits the adsorption of LEA and consequently decreases the reaction rate.

4. Conclusion

It was possible to include the TPA in materials composed of magnetite encapsulated in silicon oxide (Fe3O4@SiO2). The synthesis was monitored and standardized by analyzing the size and morphology of the NPMs by SEM images, determining that they are uniform spheres and they are not affected by factors such as the acidic medium or temperature if they are calcined in an inert atmosphere. The catalyst named Fe3O4@SiO2@SiO2/TPA400N2 was shown to be active in the synthesis of BLE from LEA and BOH, a product of interest in biomass valorization, with a conversion of 95%, a selectivity of 100% at 353 K, and up to 4 consecutive reaction cycles without appreciable loss of activity due to poisoning or leaching of the active phase. The influence of various reaction parameters on LEA conversion was studied, determining that the reaction is of second order and was validated successfully by the proposed kinetic model. An activation energy of 17 kcal/mol was determined. A plausible reaction mechanism based on the ethoxide route is outlined.

Data Availability

The datasets used to support the finding of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.

Acknowledgments

This work was financed by the projects of CONICET (PIP 003), ANPCyT, Universidad Nacional de La Plata, and COLCIENCIAS 110965843004, contract: Colciencias-UPTC 047-2015.

Supplementary Materials

The supplementary material file contains the graphical abstract, the estimated values of the effective diffusion coefficient of levulinic acid (LEA) in n-butanol, additional experimental details (acidity measurements, FT-IR TPA analysis, and reaction product quantification), and related bibliography. (Supplementary Materials)