Synthesis of Propylene Glycol Methyl Ether Acetate: Reaction Kinetics and Process Simulation Using Heterogeneous Catalyst

: Propylene glycol methyl ether acetate (PGMEA) serves as a crucial solvent in semiconductor and display material processes, demanding high purity and low acidity. Despite its significance, its conventional synthesis method using homogeneous catalysts requires extensive purification. Our study explores the use of Amberlyst-15, a stable solid catalyst, to streamline this process. Through batch reactions with a 1:1 reactant ratio at various temperatures and modeling using an integrated reaction rate equation, we obtained kinetic parameters. These parameters were used to predict the kinetics under different reactant ratios and different catalyst amounts, and the predictions match well with experimental results, especially when we used the catalyst amount scaled by the amount of the limiting reactant (PGME) rather than the total amount of the reactants. This highlights the importance of reporting kinetic parameters with proper scaling for catalyst used. Furthermore, we integrated these parameters into process simulations to determine the length of a plug flow reactor (PFR), constructed a PFR system, and confirmed that the simulation results matched well with experimental data obtained from the PFR system. Our findings suggest Amberlyst-15’s potential in simplifying PGMEA synthesis, promising advancements in industrial applications.


Introduction
Propylene glycol methyl ether acetate (PGMEA), also known as PMA or methoxy propyl acetate, serves as a vital solvent with diverse applications across several industries.Its use ranges from electronic chemicals in semiconductor manufacturing to ultra-high purity applications [1], ultra-low acidity [2,3], and reduced toxicity compared to ethylene glycol-based products [4].PGMEA is also used as a solvent in paints, inks, dyes, cleaning agents, and photoresists, making it a preferred choice due to its low toxicity profile and minimal particle formation.Notably, PGMEA is a key component in various coating formulations [5].PGMEA synthesis primarily involves the esterification reaction of propylene glycol methyl ether (PGME or methoxy propanol) with acetic acid (AA).However, this process encounters challenges due to chemical equilibrium limitations, resulting in low PGME conversion rates.The synthesis reaction is represented by the following equilibrium equation: 2 of 14 Processes 2024, 12, x FOR PEER REVIEW 2 of 14 (1) where  and  denote the forward and reverse reaction rate constants, respectively, and A through D are PGME, AA, PGMEA, and water, respectively.As an equilibrium reaction, it necessitates the use of an acidic catalyst to facilitate the reaction kinetics.
Homogeneous catalysts, such as sulfuric acid and hydrofluoric acid, have long been favored for their high catalytic activity and selectivity in various chemical reactions.However, their use comes with drawbacks, challenges in catalyst separation, process complexity, and corrosion issues.In contrast, heterogeneous acidic catalysts have emerged as promising alternatives due to their high catalytic activity, low corrosivity, and ease of separation [6,7].Solid acids, in particular, have gained traction across a spectrum of chemical reactions, including esterification [8,9], alkanes isomerization [10], aldol condensation [11] and ketal reaction [12].Their utility lies in their ability to enhance kinetics inherent in esterification reactions.Recent studies have introduced innovative methods, such as reactive distillation (RD) and pressure swing techniques, aimed at improving energy efficiency and reducing production costs in PGMEA synthesis [13,14].The RD reaction process, comprising a transesterification reaction coupled with a distillation column for azeotrope separation and energy optimization, offers a viable solution to enhance the efficiency of esterification processes [4,15].
Furthermore, the availability of heterogeneous catalysts enables the production of less-acidic electronic-grade PGMEA, contributing to the enhancement of product quality [2][3][4].Amberlyst ® 15 (Amberlyst-15) is a widely used catalyst in various chemical processes, renowned for its excellent physical, thermal, and chemical stability [16].Available in bead form, Amberlyst-15 is insoluble in water, making it suitable for non-aqueous catalysis-a crucial aspect, particularly in reactions where water is a byproduct and catalyst dissolution is undesirable.On the other hand, this catalyst can dissolve or swell in other types of liquid.For example, Amberlryst-15 will swell in acetone, so this catalyst should be used in reactions where the environment is favorable for the catalyst.Its proven stability in harsh environments, including acidic conditions, makes it an attractive option for catalytic reactions.Despite its established stability and compatibility, only limited studies have explored the use of Amberlyst-15 in synthesizing PGMEA.Further investigation into its catalytic performance in PGMEA synthesis could unveil its potential benefits.
Oh et al. [17] conducted PGMEA synthesis with Amberlyst-15 catalyst using a batch system, exploring various [PGME]:[AA] ratios and temperatures to determine reaction kinetics.Subsequently, they applied these findings to design and operate reactive distillation processes.Despite observing two reaction sites inside the column with larger injection volumes, their study achieved nearly 100% conversion at lower flow rates and smaller injection volumes.However, the batch system model fitting was not sufficiently robust, warranting further clarification of parameter determination methods.Gadekar-Shinde et al. [4] proposed an RD process for PGMEA production, exploring reaction rates at different temperatures and reactant mole ratios.They introduced toluene as an entrainer to enhance product purity without compromising catalyst stability, presenting a more compact and cost-effective process.Agrawal et al. [2] optimized a simulated moving bed reactor (SMBR) for PGMEA synthesis using Amberlyst-15.Their model revealed that increased AA conversion negatively impacted PGMEA production rates, emphasizing the importance of process operation for high PGMEA recovery and reduced downstream separation costs.Utilizing kinetic model parameters from previous studies, they optimized the SMBR, observing similar phenomena of two reaction sites and enhanced conversion beyond equilibrium in chromatograph reactors. (1) where k 1 and k 2 denote the forward and reverse reaction rate constants, respectively, and A through D are PGME, AA, PGMEA, and water, respectively.As an equilibrium reaction, it necessitates the use of an acidic catalyst to facilitate the reaction kinetics.
Homogeneous catalysts, such as sulfuric acid and hydrofluoric acid, have long been favored for their high catalytic activity and selectivity in various chemical reactions.However, their use comes with drawbacks, challenges in catalyst separation, process complexity, and corrosion issues.In contrast, heterogeneous acidic catalysts have emerged as promising alternatives due to their high catalytic activity, low corrosivity, and ease of separation [6,7].Solid acids, in particular, have gained traction across a spectrum of chemical reactions, including esterification [8,9], alkanes isomerization [10], aldol condensation [11] and ketal reaction [12].Their utility lies in their ability to enhance kinetics inherent in esterification reactions.Recent studies have introduced innovative methods, such as reactive distillation (RD) and pressure swing techniques, aimed at improving energy efficiency and reducing production costs in PGMEA synthesis [13,14].The RD reaction process, comprising a transesterification reaction coupled with a distillation column for azeotrope separation and energy optimization, offers a viable solution to enhance the efficiency of esterification processes [4,15].
Furthermore, the availability of heterogeneous catalysts enables the production of lessacidic electronic-grade PGMEA, contributing to the enhancement of product quality [2][3][4].Amberlyst ® 15 (Amberlyst-15) is a widely used catalyst in various chemical processes, renowned for its excellent physical, thermal, and chemical stability [16].Available in bead form, Amberlyst-15 is insoluble in water, making it suitable for non-aqueous catalysis-a crucial aspect, particularly in reactions where water is a byproduct and catalyst dissolution is undesirable.On the other hand, this catalyst can dissolve or swell in other types of liquid.For example, Amberlryst-15 will swell in acetone, so this catalyst should be used in reactions where the environment is favorable for the catalyst.Its proven stability in harsh environments, including acidic conditions, makes it an attractive option for catalytic reactions.Despite its established stability and compatibility, only limited studies have explored the use of Amberlyst-15 in synthesizing PGMEA.Further investigation into its catalytic performance in PGMEA synthesis could unveil its potential benefits.
Oh et al. [17] conducted PGMEA synthesis with Amberlyst-15 catalyst using a batch system, exploring various [PGME]:[AA] ratios and temperatures to determine reaction kinetics.Subsequently, they applied these findings to design and operate reactive distillation processes.Despite observing two reaction sites inside the column with larger injection volumes, their study achieved nearly 100% conversion at lower flow rates and smaller injection volumes.However, the batch system model fitting was not sufficiently robust, warranting further clarification of parameter determination methods.Gadekar-Shinde et al. [4] proposed an RD process for PGMEA production, exploring reaction rates at different temperatures and reactant mole ratios.They introduced toluene as an entrainer to enhance product purity without compromising catalyst stability, presenting a more compact and cost-effective process.Agrawal et al. [2] optimized a simulated moving bed reactor (SMBR) for PGMEA synthesis using Amberlyst-15.Their model revealed that increased AA conversion negatively impacted PGMEA production rates, emphasizing the importance of process operation for high PGMEA recovery and reduced downstream separation costs.Utilizing kinetic model parameters from previous studies, they optimized the SMBR, observing similar phenomena of two reaction sites and enhanced conversion beyond equilibrium in chromatograph reactors.
Huang et al. [18] investigated the synthesis of PGMEA using a solid acidic catalyst (SO 4 2− /TiO 2 ) and optimized its preparation.Their study demonstrated the catalyst's high stability and consistent performance.Oh et al. [3] and Wang et al. [15] explored PGMEA synthesis via transesterification reactions using different catalysts, observing varying conversion rates and catalyst deactivation.Wang et al. [15] utilized kinetic data to design and operate RD systems while Oh et al. [3] investigated the synthesis of PGMEA using ethyl acetate instead of acetic acid as a precursor.Fan et al. [13] investigated PGMEA synthesis via transesterification of methyl acetate and PGME using an RD process.They provided transesterification kinetics and thermodynamic parameters, validating the reaction feasibility experimentally in a batch RD column.
In this study, we utilized Amberlyst-15, with the aim of determining its reaction kinetics and investigating its applicability in continuous processes based on the kinetics.The catalyst was used in its commercially available form without any additional chemical treatment, mirroring common industrial practices.Our experimental approach involved studying the reaction kinetics of the process described in Equation (1) across a range of temperatures (60 to 100 • C) using a 1:1 ratio of PGME to AA and a fixed amount of Amberlyst-15 in a batch system.We developed an integrated reaction rate model, which was not available in previous studies, and extracted kinetic model parameters, including Arrhenius parameters, from the experimental data.Subsequently, we conducted additional experiments with different PGME to AA ratios and varied amounts of Amberlyst-15, comparing the resulting data with model predictions based on the earlier determined kinetic parameters.To account for the influence of catalyst amount on the reaction rate constant, we scaled the kinetic parameters using two distinct methods of expressing catalyst amount and showed which method is better.Finally, we employed the determined kinetic parameters to simulate the behavior of a plug flow reactor, informing the design process for reactor systems in potential industrial applications.Our study provides insights into the reaction kinetics of the PGME and AA synthesis process and highlights the potential utility of Amberlyst-15 as a catalyst in industrial settings.The novel features of this study are providing kinetic parameters with an explicit reaction equation, which are readily usable for designing reaction systems for PGMEA and suggesting the method to scale the amount of catalyst for better model predictions.Future research directions may involve further validation of the proposed kinetic model and exploration of alternative catalyst formulations for enhanced reaction efficiency.
Figure 1 illustrates the batch reaction system utilized to determine the reaction kinetics.The round-bottom flask has a 2 L volume with a 3-neck apparatus including a reflux condenser, heater, and magnetic stirrer.It also features a temperature controller and a sampling port, through which sampling took place regularly using a syringe with a needle.To minimize the loss of reactants and products during experiments, chilled water at 5 • C was circulated through the reflux condenser.We took a specified mass of each reactant, AA and PGME, to achieve a desired reactant ratio, along with approximately a 1 kg total mass.Each component was then introduced into the batch reactor accordingly.The ratio of the reactants varied for different tests.The reactor's temperature was controlled to one of five predetermined temperatures: 60, 70, 80, 90, or 100 • C. Upon reaching the designated reaction temperature, a desired amount of catalyst was added to initiate the reaction.Stirring of the reaction mixture at 300 rpm commenced immediately after catalyst addition, with the stirring speed chosen based on previous studies indicating its minimal impact on reaction kinetics, especially the concentration of PGMEA [17].Samples were extracted from the reactor at specified time intervals to monitor the reaction progress.It was assumed that the reaction ceased in the sample due to the absence of catalyst.Gas chromatography (GC) analysis using a Varian CP3800 (Palo Alto, CA, USA) equipped with a flame ionization detector (FID) was employed to determine the composition of each sample.The GC conditions, summarized in Table 1, included parameters such as column temperature, injection volume, and carrier gas flow rate.Before conducting the experiment, an external calibration was conducted using the same chemicals (AA, PGME, and PGMEA) with known concentrations.The resulting calibration curve was utilized to determine the concentrations of these chemicals during the experiments.No internal standard was employed.An example of chromatogram is shown in the Supplementary Materials.
Processes 2024, 12, x FOR PEER REVIEW 4 of 14 lyst addition, with the stirring speed chosen based on previous studies indicating its minimal impact on reaction kinetics, especially the concentration of PGMEA [17].Samples were extracted from the reactor at specified time intervals to monitor the reaction progress.It was assumed that the reaction ceased in the sample due to the absence of catalyst.Gas chromatography (GC) analysis using a Varian CP3800 (Palo Alto, CA, USA) equipped with a flame ionization detector (FID) was employed to determine the composition of each sample.The GC conditions, summarized in Table 1, included parameters such as column temperature, injection volume, and carrier gas flow rate.Before conducting the experiment, an external calibration was conducted using the same chemicals (AA, PGME, and PGMEA) with known concentrations.The resulting calibration curve was utilized to determine the concentrations of these chemicals during the experiments.No internal standard was employed.An example of chromatogram is shown in the Supplementary Materials.Table 1.Operating conditions for GC.We conducted simulations using Aspen Plus 2006.5, employing a reactor model of RPLUG.Our aim was to design a plug flow reactor (PFR) by integrating kinetic parameters obtained from the batch system experiments performed in this investigation.The simulations were carried out at a temperature of 90 °C, with the bed porosity and density of the catalyst particles set at 0.3 and 0.5 g/cc, respectively.The reactor dimensions were specified with a diameter of 0.1 m and a length of 1.5 m.The flowsheet and conditions are detailed in the Supplementary Materials.Accordingly, a PFR system (Figure 2) was constructed, comprising a PFR reactor with dimensions of 0.1 m in diameter and 1.5 m in length.We replicated identical operational conditions to those used in the simulation to obtain the mole fractions of PGME and PGMEA at two sampling locations along the PFR: 0.5 m and 1.5 m (i.e., outlet) to compare the simulation and experimental results.We conducted simulations using Aspen Plus 2006.5, employing a reactor model of RPLUG.Our aim was to design a plug flow reactor (PFR) by integrating kinetic parameters obtained from the batch system experiments performed in this investigation.The simulations were carried out at a temperature of 90 • C, with the bed porosity and density of the catalyst particles set at 0.3 and 0.5 g/cc, respectively.The reactor dimensions were specified with a diameter of 0.1 m and a length of 1.5 m.The flowsheet and conditions are detailed in the Supplementary Materials.Accordingly, a PFR system (Figure 2) was constructed, comprising a PFR reactor with dimensions of 0.1 m in diameter and 1.5 m in length.We replicated identical operational conditions to those used in the simulation to obtain the mole fractions of PGME and PGMEA at two sampling locations along the PFR: 0.5 m and 1.5 m (i.e., outlet) to compare the simulation and experimental results.

Kinetic Modeling
Activity coefficients, often estimated using UNIQUAC or UNIFAC methods, are utilized to correct for deviations from ideal behavior [19][20][21].For instance, Pöpken et al. [22] demonstrated that utilizing activities obtained from UNIQUAC calculations, rather than mole fractions, provided a better representation of experimental values in the synthesis of methyl acetate using Amberlyst-15.Despite the prevalence of activity-based models, kinetic models frequently rely on concentrations rather than activities [23][24][25][26].In our study, we adopted liquid concentrations (i.e., mole fractions) in the development of our kinetic model for simplicity, with a discussion of its predictive capabilities to follow.
We employed a pseudo-homogeneous model [22], assuming the ion-exchange resin serves as a source of solvated protons.The differential molar balance governing the formation of PGMEA can be expressed as a second-order kinetic model, representing a firstorder reaction with respect to each component.For our reversible reaction system, the reaction rate is modeled as where [i] (mol/L) is the concentration of a component i, (A, B, C, or D), as a function of time, and subscripts 1 and 2 denote forward and reverse reactions, respectively.The forward and reverse reaction rate constants,  and  , are expressed in units of L/(mol min).The concentration of each component in the reaction system was determined at 20 °C within the total volume of all components, which amounted to 1034 mL.The equilibrium constant,  , is defined as the ratio of the forward and reverse reaction rate constants at equilibrium and can be expressed as where the subscript  denotes equilibrium.

For the general case of [A] ≤ [B]
, we can derive the integrated rate equation for Equation (1), i.e., the concentration of a component at a certain time during the reaction where z, a, and b are defined as

Kinetic Modeling
Activity coefficients, often estimated using UNIQUAC or UNIFAC methods, are utilized to correct for deviations from ideal behavior [19][20][21].For instance, Pöpken et al. [22] demonstrated that utilizing activities obtained from UNIQUAC calculations, rather than mole fractions, provided a better representation of experimental values in the synthesis of methyl acetate using Amberlyst-15.Despite the prevalence of activity-based models, kinetic models frequently rely on concentrations rather than activities [23][24][25][26].In our study, we adopted liquid concentrations (i.e., mole fractions) in the development of our kinetic model for simplicity, with a discussion of its predictive capabilities to follow.
We employed a pseudo-homogeneous model [22], assuming the ion-exchange resin serves as a source of solvated protons.The differential molar balance governing the formation of PGMEA can be expressed as a second-order kinetic model, representing a first-order reaction with respect to each component.For our reversible reaction system, the reaction rate is modeled as where [i] (mol/L) is the concentration of a component i, (A, B, C, or D), as a function of time, and subscripts 1 and 2 denote forward and reverse reactions, respectively.The forward and reverse reaction rate constants, k 1 and k 2 , are expressed in units of L/(mol min).The concentration of each component in the reaction system was determined at 20 • C within the total volume of all components, which amounted to 1034 mL.The equilibrium constant, K e , is defined as the ratio of the forward and reverse reaction rate constants at equilibrium and can be expressed as where the subscript e denotes equilibrium.
For the general case of [A] 0 ≤ [B] 0 , we can derive the integrated rate equation for Equation (1), i.e., the concentration of a component at a certain time during the reaction Processes 2024, 12, 865 6 of 14 where z, a, and b are defined as where [A] 0 and [B] 0 are the initial concentrations of components A and B, respectively.Then, the concentration of the other components at a certain time can be modeled as follows: [ We used [A] 0 and [B] 0 and measured [A] e and [B] e to calculate a, b, and z(t) with an initial guess of k 1 and eventually to obtain model [A](t) (Equation (4)).Using the least square method to minimize the sum of squared error, which is the difference between model [A](t) and data [A](t), k 1 was obtained (and naturally k 2 as well).
For the case of [A] 0 = [B] 0 , while Equations ( 4) and ( 5) can be used as those are, Equations ( 6) and ( 7) simplify to Equations ( 3) through ( 5) and ( 8) through (11) were utilized to determine the reaction rate constants by fitting experimental data obtained with a reactants' mole ratio of 1:1 at various temperatures to the model represented by Equation (4).This fitting assumed a pseudo-homogeneous reaction.Subsequently, the same reaction rate constants and Equations (3) through (9) were employed to forecast the reaction kinetics with a reactants' mole ratio other than 1:1, thereby validating the previously determined reaction rate constants.
The scaled forward and reverse reaction rate constants are expressed in units of L mol-PGME/(mol min g-catalyst).These constants are calculated using the following equations: where m cat in g-catalyst/mol-reactant is the catalyst amount relative to the number of moles of reactant (either total or limiting reactant).The scaled reaction rate constants k f and k r should remain unchanged regardless of the specific amount of catalyst utilized in the reaction system.As long as the reactant amount is used consistently, whether it is divided by the total amount of all the reactants or the amount of the limiting reactant, these constants can be universally applied across different catalyst amounts without alteration.
For reactions involving a different amount of catalyst, the prediction can be made using Equation (4) with the following updated reaction rate constants: where m ′ cat in g-catalyst/mol-reactant is the new amount of catalyst different from that used in this study, k ′ 1 in L/(mol min) represents the new forward reaction rate constant for the catalyst amount, and k f in L mol-PGME/(mol min g-catalyst) is the scaled forward reaction rate constant determined in this study.The new reverse reaction rate constant can be obtained by where K ′ e is the equilibrium constant obtained in the reaction with the catalyst amount, m ′ cat .Thus, Equation (4) should be reduced to the following universal equation: Again, Equations ( 5) through ( 9) and ( 16) were employed to forecast the reaction kinetics with catalyst amounts different from what used to determine k f and k r , thereby validating the previously determined reaction rate constants.
The effect of temperature on the reaction rate constants can be described by the Arrhenius equation: where E a,1 and E a,2 (J/mol) represent the activation energies for the forward and reverse reactions, respectively, R is the gas constant, equal to 8.31446 J/(mol K), A 1 and A 2 in L/(mol min) are the pre-exponential factors for the forward and reverse reaction rate constants, respectively, and A f and A r in L mol-PGME/(mol min g-catalyst) are the preexponential factors for the scaled forward and reverse reaction rate constants, respectively, as follows:

Effect of Temperature and Kinetic Parameters
First, we conducted PGMEA synthesis reactions with a 1:1 reactants mole ratio (i.e., [PGME] 0 = [AA] 0 ) at various temperatures.Figure 3 illustrates the effect of temperature on the concentration of PGME and the product PGMEA with time at five different temperatures ranging from 60 to 100 • C with a catalyst amount of 14.5 g-catalyst/mol-PGME (equivalent to 7.3 g-catalyst/mol-reactants or 0.1 g-catalyst/mL-reactants).The experiments were conducted over a significantly longer duration than depicted in Figure 3.However, for brevity and clarity, we present data only up to concentrations equivalent to those observed at later times when there was no more change in the concentrations.Consequently, the final concentrations depicted represent the equilibrium compositions.As expected, higher temperatures lead to shorter times to reach equilibrium, consistent with the well-known effect of temperature on reaction rates.The equilibrium concentrations of each reactant are approximately equal to that of the product, i.e., K e ≈ 1.
stants within the range of operating temperatures, except at 60 °C, where both rates are about the same.However, the activation energy for the reverse reaction is higher than that for the forward reaction, suggesting that the reaction rates of the reverse reaction are more temperature-sensitive compared to the forward reaction.Thus, at lower temperatures, such as 50 °C or below, we anticipate that the forward reaction rate will exceed that of the reverse reaction, leading to higher conversion rates.The lines are Arrhenius model predictions, Equations ( 17) and ( 18), with R 2 = 0.9981 and R 2 = 0.9984 for forward and reverse reactions, respectively.4) and ( 9), for PGMEA.Data and predictions after reaching equilibrium are not shown to avoid confusion.
Equations ( 3) through ( 5) and ( 8) through (11) were utilized to determine the reaction rate constants by fitting the data obtained with a reactants' mole ratio of 1:1 to the model, Equation (4), assuming a pseudo-homogeneous reaction.Figure 3 also presents the model predictions for both PGME and PGMEA.The former was obtained using Equation (4), while the latter was obtained using Equations (4) and (9).Generally, the data align well with the model predictions.The reaction rate constants as a function of temperature were naturally obtained and are depicted in Figure 4.These data were fitted to Equations ( 17) and (18) to determine the activation energy and pre-exponential factors.The resulting values are tabulated in Table 2.
while the latter was obtained using Equations ( 4) and (9).Generally, the data align well with the model predictions.The reaction rate constants as a function of temperature were naturally obtained and are depicted in Figure 4.These data were fitted to Equations ( 17) and ( 18) to determine the activation energy and pre-exponential factors.The resulting values are tabulated in Table 2.
The reverse reaction rate constants tend to be higher than the forward reaction constants within the range of operating temperatures, except at 60 °C, where both rates are about the same.However, the activation energy for the reverse reaction is higher than that for the forward reaction, suggesting that the reaction rates of the reverse reaction are more temperature-sensitive compared to the forward reaction.Thus, at lower temperatures, such as 50 °C or below, we anticipate that the forward reaction rate will exceed that of the reverse reaction, leading to higher conversion rates.The lines are Arrhenius model predictions, Equations ( 17) and ( 18), with R 2 = 0.9981 and R 2 = 0.9984 for forward and reverse reactions, respectively.17) and ( 18), with R 2 = 0.9981 and R 2 = 0.9984 for forward and reverse reactions, respectively.The reverse reaction rate constants tend to be higher than the forward reaction constants within the range of operating temperatures, except at 60 • C, where both rates are about the same.However, the activation energy for the reverse reaction is higher than that for the forward reaction, suggesting that the reaction rates of the reverse reaction are more temperature-sensitive compared to the forward reaction.Thus, at lower temperatures, such as 50 • C or below, we anticipate that the forward reaction rate will exceed that of the reverse reaction, leading to higher conversion rates.
It is intriguing to juxtapose our findings with those of other researchers using a similar system.Huang et al. [18] obtained kinetic parameters at higher temperatures than those investigated in this study.Notably, their reported reaction rate constants are not scaled with the catalyst amount.Thus, Figure 5 presents a comparison for the forward reaction rate constant without scaling, expressed in units of L/(min mol).Generally, the rate constants obtained in our study are higher than those reported by Huang et al. [18] even though those have the same order of magnitude.While those studies show different reaction rate constants, both datasets exhibit similar slopes, indicating comparable activation energies (Table 2).

Validation of Kinetic Parameters with Various Initial Concentrations and Catalyst Amount
It is crucial to assess whether the parameters for reaction kinetics obtained in the previous section for the 1:1 reactant mole ratio are still applicable for other reactant mole ratios and different catalyst amounts.Figure 6a displays the concentration profiles of two reactants and PGMEA over time for a [PGME]:[AA] ratio of 1:4 with 37 g-catalyst/mol-PGME at 80 °C.The model predictions are based on Equations (3) through (9) using the same scaled reaction rate constants ( ) derived in the previous section for the 1:1 reactant mole ratio with 14.5 gcatalyst/mol-PGME at the same temperature.The new reaction rate constants were obtained using Equations ( 14) and (15).Similarly, Figure 6b depicts the concentration profiles for a [PGME]:[AA] ratio of 1:2 with 20 g-catalyst/mol-PGME at 100 °C.The predictions were obtained using the same approach.Generally, the model predicts the behavior quite accurately 1 in this study with a similar study [18].
On the other hand, Oh et al. [17] and Gadekar-Shinde et al.
[4] did not explicitly detail how their reaction rate constants were scaled, making a direct comparison challenging.However, it is important to note that the method of scaling should primarily affect the pre-exponential factor rather than the activation energy.As shown in Table 2, the activation energy for the forward reaction remains consistent across the four studies, while that for the reverse reaction displays considerable variation.Particularly, Gadekar-Shinde et al. [4] reported a notably low activation energy for the reverse reaction, achieving high conversions.Their methodology involved stirring the reaction mixture at a significantly higher rate than in our study to minimize mass transfer resistance.Additionally, they applied chemical treatment to the catalyst, which could account for the observed differences.The challenge of achieving high conversion in esterification is closely tied to the removal of water from the reaction mixture, as water is one of the products in this reversible reaction.Le Chatelier's principle predicts that removing water will drive the equilibrium towards the formation of products, thus enhancing conversion.Additionally, modifying the catalyst through chemical treatments can also induce changes in its chemistry, further influencing the reaction dynamics.However, the implementation of an additional chemical treatment facility can be costly.Moreover, stirring the reaction mixture at high speeds in an industrial setting can pose risks.Considering these factors, along with the potential for uncontrolled mass transfer resistance in a plug flow reactor (PFR), we chose a lower stirring speed that does not influence the reaction [17], and furthermore, we used the catalyst in its original state without any chemical treatment.Our aim was to provide reaction kinetic parameters applicable to similar industrial setups.

Validation of Kinetic Parameters with Various Initial Concentrations and Catalyst Amount
It is crucial to assess whether the parameters for reaction kinetics obtained in the previous section for the 1:1 reactant mole ratio are still applicable for other reactant mole ratios and different catalyst amounts.Figure 6a displays the concentration profiles of two reactants and PGMEA over time for a [PGME]:[AA] ratio of 1:4 with 37 g-catalyst/mol-PGME at 80 • C. The model predictions are based on Equations (3) through (9) using the same scaled reaction rate constants (k f ) derived in the previous section for the 1:1 reactant mole ratio with 14.5 g-catalyst/mol-PGME at the same temperature.The new reaction rate constants were obtained using Equations ( 14) and (15).Similarly, Figure 6b depicts the concentration profiles for a [PGME]:[AA] ratio of 1:2 with 20 g-catalyst/mol-PGME at 100 • C. The predictions were obtained using the same approach.Generally, the model predicts the behavior quite accurately despite changes in both the reactant ratio and catalyst amount.These results suggest that the kinetic parameters (k f and thus A and E a ) can be universally applied to design reactors or predict reaction kinetics without using the activity coefficient.
The catalyst amount should impact kinetics rather than equilibrium, influencing the time required to reach equilibrium [18].Therefore, the method of scaling the reaction rate constant will affect predictions on the time axis.A pertinent question arises regarding the optimal way to scale with the catalyst amount: should it be relative to the total amount of the reactants or the amount of the limiting reactant (PGME in this study)?Figure 6a,b, based on the former case, indicate predictions slower than the data in terms of time to reach equilibrium.Conversely, Figure 6c,d, based on the latter case, show predictions faster than the former case.When employing a specific catalyst mass for two different reactions-one with a limiting reactant and the other with a 1:1 reactant ratio-while maintaining the same total number of moles of reactants, the catalyst amount per the limiting reactant in the former case will appear higher than that in the latter case even though the catalyst amount per the total amount of reactants is the same for both cases.Since the reaction is controlled by the limiting reactant, the reaction should proceed faster from the limiting reactant's perspective in the former case compared to the reaction with a 1:1 reactant ratio if the absolute amount of catalyst remains constant.Thus, while kinetic data in Figure 6a,b outpace the model predictions using the catalyst amount calculated based on the total amount of reactants, the predictions using the catalyst amount calculated based on the amount of the limiting reactant as shown in Figure 6c,d better align with the data.Therefore, when reporting or using parameters for reaction kinetics, the pre-exponential factor or reaction rate constants scaled with the catalyst amount per the amount of the limiting reactant need to be reported or used.Thus, m cat and m ′ cat in Equations ( 12) through ( 14), ( 16), (21), and ( 22) should be as follows: mass of catalyst moles of limiting reactant (23) and m ′ cat = mass of catalyst moles of limiting reactant (24) This naturally leads to using A * f and A * r instead of A + f and A + r in Table 2 to predict kinetics or design a reactor.This naturally leads to using  * and  * instead of  and  in Table 2 to predict kinetics or design a reactor.

Simulation and Data for Plug Flow Reactor System
A plug flow reactor (PFR) presents a practical solution for employing solid catalysts and effectively controlling reaction temperature in chemical processes.In this study, we utilized Aspen Plus 2006.5 to conduct simulations with a reactor model of RPLUG to design a PFR system, incorporating the kinetic parameters obtained from the batch system experiments conducted in this investigation.The simulations were conducted at a temperature of 90 °C, with the bed porosity and density of the catalyst particles set at 0.3 and 0.5 g/cc, respectively.The reactor dimensions were set as a diameter of 0.1 m and a length of 1.5 m.The feed composition was characterized by the mole fractions of

Simulation and Data for Plug Flow Reactor System
A plug flow reactor (PFR) presents a practical solution for employing solid catalysts and effectively controlling reaction temperature in chemical processes.In this study, we utilized Aspen Plus 2006.5 to conduct simulations with a reactor model of RPLUG to design a PFR system, incorporating the kinetic parameters obtained from the batch system experiments conducted in this investigation.The simulations were conducted at a temperature of 90 • C, with the bed porosity and density of the catalyst particles set at 0.3 and 0.5 g/cc, respectively.The reactor dimensions were set as a diameter of 0.1 m and a length of 1.5 m.The feed composition was characterized by the mole fractions of [AA]:[PGME]:[PGMEA] = 0.475:0.475:0.05,with a feed rate of 1.43 kmol/h.Figure 7 illustrates the simulated mole fraction profile for PGME and PGMEA resulting from the PFR design.Notably, equilibrium conditions were achieved within a PFR length of 1.5 m, suggesting that this length represents the minimum requirement for designing a PFR under the specified conditions.This insight allows for an efficient design approach in implementing PFR systems tailored to specific reaction requirements and for showcasing its potential to optimize reactor design while considering variables such as energy cost and manufacturing expenses.
Processes 2024, 12, x FOR PEER REVIEW 12 of 14 [AA]:[PGME]:[PGMEA] = 0.475:0.475:0.05,with a feed rate of 1.43 kmol/h.Figure 7 illustrates the simulated mole fraction profile for PGME and PGMEA resulting from the PFR design.Notably, equilibrium conditions were achieved within a PFR length of 1.5 m, suggesting that this length represents the minimum requirement for designing a PFR under the specified conditions.This insight allows for an efficient design approach in implementing PFR systems tailored to specific reaction requirements and for showcasing its potential to optimize reactor design while considering variables such as energy cost and manufacturing expenses.Consequently, a PFR system (Figure 2) was constructed with identical dimensions to those employed in the simulation.To account for the swelling of the catalyst, we swelled the catalyst in deionized water and then filled the PFR reactor with the swollen catalyst up to 1.5 m.By applying the same operational conditions, we determined the mole fractions of PGME and PGMEA at two sampling points within the PFR reactor.These results closely correlate with the simulated outcomes illustrated in Figure 7. Hence, it is apparent that kinetic parameters obtained from batch reactor experiments can significantly contribute to simulation endeavors aimed at designing and constructing alternative reactor configurations.This alignment between simulation results and real-world operation underscores the effectiveness of utilizing such parameters in optimizing the performance of novel reaction systems.

Conclusions
Through our study, we conducted the synthesis of PGMEA using AA and PGME alongside a commercially available catalyst, Amberlyst-15.This endeavor aimed to elucidate the critical parameters governing reaction kinetics, alongside the development of an integrated kinetic model.These aspects are both novel and indispensable in the design of other types of reactors.Successfully modeling the reaction kinetics across various temperatures with a 1:1 reactant ratio, we employed an integrated second-order reversible reaction equation along with the Arrhenius model.This approach allowed for the natural determination of kinetic and Arrhenius parameters, pivotal for understanding reaction dynamics.Utilizing these parameters, we extended our analysis to predict reaction kinetics for scenarios involving different reactant ratios and catalyst amounts.Notably, our predictions aligned well with experimental data, particularly when the catalyst amount was scaled by the amount of PGME, the limiting reactant, rather than that of all the reactants.Consequently, a PFR system (Figure 2) was constructed with identical dimensions to those employed in the simulation.To account for the swelling of the catalyst, we swelled the catalyst in deionized water and then filled the PFR reactor with the swollen catalyst up to 1.5 m.By applying the same operational conditions, we determined the mole fractions of PGME and PGMEA at two sampling points within the PFR reactor.These results closely correlate with the simulated outcomes illustrated in Figure 7. Hence, it is apparent that kinetic parameters obtained from batch reactor experiments can significantly contribute to simulation endeavors aimed at designing and constructing alternative reactor configurations.This alignment between simulation results and real-world operation underscores the effectiveness of utilizing such parameters in optimizing the performance of novel reaction systems.

Conclusions
Through our study, we conducted the synthesis of PGMEA using AA and PGME alongside a commercially available catalyst, Amberlyst-15.This endeavor aimed to elucidate the critical parameters governing reaction kinetics, alongside the development of an integrated kinetic model.These aspects are both novel and indispensable in the design of other types of reactors.Successfully modeling the reaction kinetics across various temperatures with a 1:1 reactant ratio, we employed an integrated second-order reversible reaction equation along with the Arrhenius model.This approach allowed for the natural determination of kinetic and Arrhenius parameters, pivotal for understanding reaction dynamics.Utilizing these parameters, we extended our analysis to predict reaction kinetics for scenarios involving different reactant ratios and catalyst amounts.Notably, our predictions aligned well with experimental data, particularly when the catalyst amount was scaled by the amount

Figure 2 .
Figure 2. PFR system with the same dimensions used in simulation.

Figure 2 .
Figure 2. PFR system with the same dimensions used in simulation.

Figure 3 .
Figure 3. Effects of temperature and time on the concentrations () of PGME and PGMEA during reaction with 14.5 g catalyst/mol PGME, which is equivalent to 7.3 g catalyst/mol reactants or 0.1 g catalyst/mL reactants.Dashed lines are model predictions, Equation (4), for PGME, while solid lines are model predictions, Equations (4) and (9), for PGMEA.Data and predictions after reaching equilibrium are not shown to avoid confusion.

Figure 4 .
Figure 4. Effects of temperature on the forward and reverse reaction rate constants in L/(mol min).The lines are Arrhenius model predictions, Equations (17) and (18), with R 2 = 0.9981 and R 2 = 0.9984 for forward and reverse reactions, respectively.

Figure 3 .
Figure 3. Effects of temperature and time on the concentrations () of PGME and PGMEA during reaction with 14.5 g catalyst/mol PGME, which is equivalent to 7.3 g catalyst/mol reactants or 0.1 g catalyst/mL reactants.Dashed lines are model predictions, Equation (4), for PGME, while solid lines are model predictions, Equations (4) and (9), for PGMEA.Data and predictions after reaching equilibrium are not shown to avoid confusion.

Figure 4 .
Figure 4. Effects of temperature on the forward and reverse reaction rate constants in L/(mol min).The lines are Arrhenius model predictions, Equations (17) and (18), with R 2 = 0.9981 and R 2 = 0.9984 for forward and reverse reactions, respectively.

Figure 5 .
Figure 5.Comparison of the forward reaction rate constant for the case of [PGME]:[AA] = 1:1 in this study with a similar study [18].

Figure 5 .
Figure 5.Comparison of the forward reaction rate constant for the case of [PGME]:[AA] = 1:1 in this study with a similar study [18].

Figure 7 .
Figure 7. Simulated and experimental mole fraction profiles along with the length of the PFR at 90 °C.

Figure 7 .
Figure 7. Simulated and experimental mole fraction profiles along with the length of the PFR at 90 • C.

Table 1 .
Operating conditions for GC.

Table 2 .
Parameters of reaction kinetics.