Detailed Kinetic Modeling of CO 2 - Based Fischer–Tropsch Synthesis

: The direct hydrogenation of CO 2 to long-chain hydrocarbons, so called CO 2 -based Fischer– Tropsch synthesis (FTS), is a viable future production route for various hydrocarbons used in the chemical industry or fuel applications. The detailed modeling of the reactant consumption and product distribution is very important for further process improvements but has gained only limited attention so far. We adapted proven modeling approaches from the traditional FTS and developed a detailed kinetic model for the CO 2 -FTS based on experiments with an Fe based catalyst in a lab-scale tubular reactor. The model is based on a direct CO 2 dissociation mechanism for the reverse water gas shift and the alkyl mechanism with an H-assisted CO dissociation step for the FTS. The model is able to predict the reactant consumption, as well as the hydrocarbon distribution, reliably within the experimental range studied (10bar, 280–320°C, 900–120,000 mL N h − 1 g − 1 and H 2 /CO 2 molar inlet ratios of 2–4) and demonstrates the applicability of traditional FTS models for the CO 2 -based synthesis. Peculiarities of the fractions of individual hydrocarbon classes (1-alkenes, n-alkanes, and iso-alkenes) are accounted for with chain-length-dependent kinetic parameters for branching and dissociative desorption. However, the reliable modeling of class fractions for high carbon number products (>C 12 ) remains a challenge not only from a modeling perspective but also from product collection and analysis.


Introduction
The Fischer-Tropsch synthesis (FTS) is an established industrial process for the production of a broad range of hydrocarbons.Its main application is currently the conversion of synthesis gas (H 2 and CO), derived from coal or natural gas, to high-value transportation fuels [1].Due to extensive international efforts to substitute energy carriers derived from fossil resources, FTS gained a lot of attention for the production of sustainable hydrocarbons as CO 2 may also be used as a carbon source for the required synthesis gas [2][3][4][5].When using CO 2 , one might follow a traditional two-step approach by first generating the synthesis gas, e.g., via reverse water gas shift (RWGS), with a subsequent CO-FTS.However, by using Fe based catalysts, CO 2 also can be directly hydrogenated to long-chain hydrocarbons in one reactor as they catalyze both the RWGS and FTS [2][3][4][5][6].In this configuration, CO is formed in situ in the FTS reactor via the RWGS reaction.This will be referred to as CO 2 -FTS and is the focus of this work.
There are numerous studies dealing with the modeling of kinetics and product distribution of traditional FTS with a greatly varying degree of detail [7,8].For CO 2 -FTS, however, only a very limited number of studies has been conducted so far.Several macrokinetic models have been reported [9][10][11][12][13], including the model we developed for a supported iron catalyst [14].Very recently, Panzone et al. published the first detailed kinetic model [15].They derived explicit Langmuir-Hinshelwood-Hougen-Watson (LHHW) type kinetic expressions based on a redox mechanism for the RWGS, an alkyl-mechanism for the formation of hydrocarbons, and a CO insertion mechanism for the formation of oxygenates.To the best of our knowledge, additional detailed models have not been reported so far.
It appears reasonable to adapt proven modeling approaches from traditional FTS for the CO 2 -based pathway.However, care has to be applied regarding the mechanistic implications of the used catalyst.The product spectrum of Co based catalysts is strongly determined by secondary reactions, while, for alkalized Fe based catalysts at elevated temperatures, it is mainly determined by primary reactions [16].Modeling approaches for Co based catalysts might thus be misleading and we will only discuss studies dealing with Fe based catalysts here.
In general, FTS can be considered as a polymerization reaction with three main steps: chain initiation, chain growth and chain termination.When assuming a constant chain growth probability α for all growing species one ends up with Equation (1) for the mole fraction of hydrocarbons with n carbon atoms.This is the well known Anderson-Schulz-Flory (ASF) distribution, which yields a straight line when plotting the logarithmic mole fraction vs. the carbon number.In reality, however, deviations from this ideal model are usually observed, i.e., a higher methane selectivity, lower C 2 selectivity, and an overprediction of long-chain hydrocarbons [17,18]: Despite almost 100 years of research, the mechanism of chain growth is still under debate and several mechanisms haven been proposed throughout the decades [1,[17][18][19].Presumably, the alkyl mechanism (based on the carbide theory originally proposed by Fischer and Tropsch [20,21]) and the CO insertion mechanism are currently the most favored pathways [1,17,18].Within these mechanisms, CO dissociation is considered to be a key step and may proceed unassisted or H-assisted.Density functional theory (DFT) calculations of different Fe phases seem to support the fact that H-assisted and unassisted CO dissociation proceed in parallel with varying ratios depending on the surface configuration [22][23][24][25].
Detailed FTS kinetic models seek to achieve the correct modeling of reactant consumption, as well as product distribution, by the mathematical adaption of the theoretically derived mechanisms.Such models have been reported by several authors [26][27][28][29][30]. Special emphasis is usually focused towards the experimentally observed deviations from an ideal ASF distribution.Some authors solely focused on the product distribution [31][32][33][34].These are sometimes referred to as "hydrocarbon selectivity" models.
Secondary growth of 1-alkenes may be an important reaction path for Co based catalysts.For alkalized Fe based catalysts, however, secondary reactions of 1-alkenes are reported to be mainly hydrogenation and isomerization to 2-alkenes ( [39][40][41]).The consideration of two different FTS sites or growth mechanisms may be a well working modeling approach.Currently, it does appear to be a rather empirical approach though.A completely fundamental approach would be physically most correct but is far beyond this semi-empirical study.Hence, we applied the remaining chain-length-dependent parameter approach proposed by Botes [38].It is reported to resemble the chemically intrinsic dependence of kinetic parameters on chain length to at least some degree [30].
Summing up, we propose a detailed kinetic model for the CO 2 -based FTS on a supported Fe catalyst which is based on a direct CO 2 dissociation mechanism for the RWGS and the alkyl mechanism with chain-length dependent parameters and H-assisted CO dissociation for FTS.The model allows a reliable prediction of the reactant consumption, the overall hydrocarbon distribution, and the fractions of hydrocarbon groups (1-alkenes, n-alkanes, and iso-alkanes) up to C 12 .Oxygenates are not considered due to experimental hurdles to determine them.

Kinetic Model
The kinetic model was developed from elementary reaction steps adapted from recent DFT studies of the RWGS and FTS under the conditions of CO 2 -FTS [42][43][44].To reduce the complexity of the model, we did not assume RWGS and FTS to proceed on separate catalytic sites.Otherwise, this would have led to a significant increase of the model parameters.It has to be pointed out that modeling the RWGS and FTS sites separately is not synonymous with a 2-α approach, which would assume two separate FTS sites/growth mechanisms.
The adsorption of the two reactants, H 2 and CO 2 , was assumed to be in equilibrium (steps 1 and 2, Table 1).This assumption was not possible for CO and H 2 O.They are not present in the feed gas and one can thus not solve for the initial surface coverage with an equilibrium constant.Therefore, two rate expressions were assigned for the adsorption and desorption, respectively (step 3 and 4).H 2 and H 2 O adsorb dissociatively while CO 2 and CO adsorb associatively.

No. Description Reaction Step Equilibrium / Kinetic Expression Notes
1 Direct CO 2 dissociation mechanism, see Table 2; 2 H-assisted CO dissociation mechanism, see Table 3.
For the RWGS reaction, we adopted a direct CO 2 dissociation mechanism that was proposed in several DFT studies [42][43][44] to be the dominant pathway for the reaction conditions examined in this work.The corresponding elementary steps are given in Table 2. Due to the model design, it is not reasonable to calculate the surface coverage of multiple oxygen species.They would be completely correlated and thereby meaningless.Thus, we considered only one kind of oxygen containing species (here: OH * ), which is basically a pseudocomponent for all oxygen containing species that might be present on the catalyst surface.Thus, step 5b was assumed to be in quasi-equilibrium to express Θ O in terms of Θ OH .This led to the overall expressions for step 5 as given in Tables 1 and 2.
In our previously published macrokinetic model, we ended up with an H-assisted CO 2 dissociation pathway leading to the overall best result [14].However, it was pointed out that the model results were hardly affected by the assumed RWGS mechanism within the experimental range studied.
Table 2. Elementary reaction steps and assigned equilibrium/kinetic expressions for the direct CO 2 dissociation mechanism as proposed by Han et al. [42].

No. Elementary Step
Equilibrium / Kinetic Expression For the FTS, we applied the alkyl mechanism with H-assisted CO dissociation.Hassisted CO dissociation is reported by Nie et al. [44] to be dominant under CO 2 -FTS conditions.The elementary steps as proposed by Ojeda et al. [22] for traditional FTS are given in Table 3.In agreement with the literature, we considered step 6b as the ratedetermining step (RDS) and 6a to be in quasi-equilibrium.This yields the overall expression for step 6 as given in Tables 1 and 3.
Table 3. Elementary reaction steps and assigned equilibrium/kinetic expressions for the H-assisted CO dissociation mechanism as proposed by Ojeda et al. [22].

No. Elementary Step
Equilibrium/Kinetic Expression Notes The alkyl mechanism was extended by the possibility of chain branching as proposed by Schulz et al. [45].A schematic of the mechanism is given in Figure 1.Adsorbed alkylspecies (R n ) can either grow via the addition of a methylene unit (step 8), terminate to an n-alkane (P n , n-paraffin) via associative desorption (step 9), terminate to an 1-alkene (Ol n , α-olefin) via dissociative desorption (β-hydride elimination, step 10), or undergo branching to an iso-alkyl species (step 11, IR n ).It is assumed that branched and unbranched species grow with the same rate.The branching reaction is reported by Schulz et al. [45] to be highly dependent on chain length due to spacious restrictions.In line with their results, we assumed an exponentially decaying rate for chain branching and assigned an extra value for the C 3 -species.Additionally, compliant with the chain-length-dependent desorption model of Botes [38], we assumed an exponentially decaying rate of dissociative desorption with increasing chain length.It was necessary to include an empirical damping factor for all termination steps (ΘΘ −1 OH ), which will be discussed in Section 4. The temperature dependency of the rate constants was implemented with a reparameterized Arrhenius equation (Equation ( 2)).For dissociative desorption and branching, it was extended by a chain-length-dependence parameter Γ.It turned out that considering every reaction as temperature dependent did not improve the model and was therefore discarded.It was sufficient to describe the temperature dependency of the reactant consumption with the activation energy of the water adsorption and desorption steps.The product distribution did not display a strong temperature dependency and we only considered an activation energy for the 1-alkene termination.It is important to note that, in this case, the obtained values for the activation energies are not related to intrinsic values and have very limited physical meaning since they are lumps of several reactions: Figure 1.Schematic of the assumed chain growth mechanism with the possibility of branching, adapted from Schulz et al. [45]; R: adsorbed alkyl species, IR: adsorbed iso-alkyl species, P: n-alkane (n-paraffin), Ol: 1-alkene (α-olefin), I: iso-alkene, n: number of carbon atoms.
In total, the model consisted of 23 adaptable parameters.Secondary reactions of 1-alkenes were not considered.The possible influence is discussed in Section 3.1.

Reactor Model
The lab-scale tubular reactor was modeled as an isothermal, isobaric and pseudohomogeneous plug flow reactor (PFR) at steady state.The change of molar flow rate ṅi of each considered component i for an incremental mass of catalyst m cat is given by the sum of the reactions r j , multiplied by its respective stoichiometric coefficient ν i,j .Hydrocarbons with up to 40 carbon atoms are considered here.This leads to 120 ordinary differential equations according to Equation (3).For the calculations, Equation (3) was scaled with the inlet volumetric flow rate at normal conditions to make the calculation independent from absolute flow rates.At steady state, the catalyst surface coverage is constant at each position in the reactor over time.Thus, for all adsorbed species k, 82 algebraic equations according to Equation (4) have to be fulfilled.As the closing condition, the catalytic site balance (Equation ( 5)) has to be satisfied: The resulting semi-explicit differential-algebraic system of equations was solved with ode15s in MATLAB.Consistent initial conditions for the catalyst surface coverage were computed with ode15s as well.The modeling approach was adapted from Visconti et al. [46].

Parameter Regression
The parameters of the kinetic model were regressed with a global optimization approach using weighted least-squares.The objective function was the sum of squared relative deviations of the molar flow rates of defined components i for all experiments l (see Equation ( 6)).The relative deviations were weighted with empirical factors ω i to ensure that certain properties of the model that are of special interest (e.g., the reactant consumption) are given a higher priority.In practice, this meant a stronger weighting of the reactant flow rates to ensure a good reproduction of the conversions.In the regression process, we considered the flow rates of H 2 , CO 2 , CO, n-alkanes, 1-alkenes and iso-alkenes individually up to C 15 , and the sum of hydrocarbon flow rates from C 16 to C 30 : For the global optimization, quasirandom parameters were generated in a predefined space using sobolset in MATLAB and used as the initial guesses for a local optimization solver (lsqnonlin).

General Results and Model Details
Results of the presented kinetic model in comparison with experimental data are shown in Figure 2. The product spectrum of the applied Fe/K@γ-Al 2 O 3 catalyst (see Section 4) is similar to a high temperature FTS process with short-chain 1-alkenes as the main product.More general information regarding the product distribution is provided in our previous publication [14].The hydrocarbon distribution (a) displays the typically experimentally observed deviations from an ideal ASF distribution: a higher methane selectivity, a lower C 2 selectivity, and a larger fraction of long-chain hydrocarbons (>C 10 ).These features can be well reproduced by our model up to ≈C 30 for all experiments.Additionally, we observed a rapid decline for >C 30 in the experiments that is not accounted for by the model.This is likely due to problems in obtaining a reliable wax sample from the experiments.The fraction of hydrocarbons with more than 30 carbon atoms is extremely small (<3% on weight basis), and we also observed some accumulation of wax at edges in front of the hot trap in our setup.It is likely that the drained samples were therefore biased.
Figure 2b,c show the fractions of 1-alkenes and n-alkanes within one carbon number.The 1-alkene fraction exhibits a maximum for propene and declines rapidly after that with an inflection point around C 12 .The n-alkane fraction is almost constant up to C 10 and starts to increase for longer chains.These experimental results can be well reproduced by the model up to C 12 .The initial rapid decline of the 1-alkene fraction in the model is caused by chain-length-dependent branching while the following decrease is due to the chain-length-dependence of the termination step.When analyzing the absolute flow rates of 1-alkenes and n-alkanes (not shown), it becomes clear that the deviations for >C 12 are caused by the insufficient description of the n-alkane formation rates.We have to note that the chromatographic resolution decreased significantly in this carbon number range, which introduces additional uncertainty for the interpretation.However, we suspect that secondary hydrogenation is the main reason for this trend, which is not accounted for in the model.
We tried implementing a secondary hydrogenation reaction while also accounting for a possible vapor-liquid-equilibrium (VLE), which would cause the predominant hydrogenation of long chains because they have an increased residence time in the reactor.Upon inclusion, a slight improvement for long-chain hydrocarbons could be achieved.However, the inclusion of a sufficiently reliable VLE model (here: Predictive Soave-Redlich-Kwong) increased the computational demand drastically.From our perspective, the additional computational demand in relation to the slight improvement did not justify the extension.The main body of the product spectrum consists of hydrocarbons with up to 10 carbon atoms that can be well described and appear to be mainly controlled by primary formation mechanisms.Additionally, the integral character of the applied tubular reactor makes the analysis very challenging.The experimental results may be the product of several superimposed effects that cannot be easily deconvolved.It appears more reasonable to us to apply a CSTR-type reactor (e.g., a spinning basket reactor) for detailed mechanistic investigations.
The simulated coverage of the catalyst surface with increasing CO 2 conversion is shown in Figure 3.In our model, the catalyst surface is mainly covered with H 2 , CO 2 , hydroxyl groups, and the chain growth monomer CH 2 .The overall reaction rate is primarily controlled by the hydroxyl surface coverage, which is reflected by the fact that we only had to consider the activation energies for the water adsorption/desorption steps to reproduce the temperature effects on reactant consumption.The surface coverage of adsorbed alkyl species is <0.0001 and was thus omitted in the graph.In Section 3, we mentioned the necessity to include an empirical damping factor (ΘΘ −1 OH ) for all termination reactions.The reason for this may be inferred from Figure 3.The monomer surface coverage (CH 2 * ) decreases significantly at high conversions.This leads to a decreasing chain growth rate and would, without the damping factor, subsequently lead to a significantly shorter product.However, this effect was less pronounced in the experiments at high conversions.Such an empirical modification is highly undesirable and may not be necessary when the surface coverage of the FTS sites is separately modeled from the RWGS sites.

Parameter Estimates
The parameter estimates that were obtained from the applied regression approach are given in Table 4.We have to point out that we cannot provide a meaningful statistical analysis of the parameter uncertainty here as several necessary assumptions are not fulfilled [47] (e.g., uncorrelated responses or random and normally distributed residuals).Thus, we do not report individual confidence intervals for the parameters here (as they would be meaningless).

Parity Plots
The parity plots for CO 2 and H 2 conversion are given in Figure 4. Overall, a good agreement could be achieved between the experiments and simulation results.The average relative error of the conversions were 5% and 7% for CO 2 and H 2 , respectively.This is a good result compared to other detailed FTS models [27,28,48].In our previously published macrokinetic model [14] (which was developed upon the same experimental data), however, we could achieve a maximum relative deviation of 5% for the CO 2 conversion.Thus, when a precise reproduction of the reactant consumption is necessary, a simple macrokinetic model may be preferred.
The parity plots for CO and CH 4 selectivity are given in Figure 5.The model can reproduce the CO selectivity with a relative error of 9% and CH 4 with a relative error of 14%.For CH 4 , there seems to be a systematic error in the model for high values.This error might be caused by an insufficient description of the methane formation reaction.From our point of view, though, it is more likely that this problem is caused by the assumption of identical catalyst sites for all reactions, which also urged us to implement the damping factor.
The parity plots for the selectivities to ethene, propene, 1-butene and 1-pentene are given in Figure 6.Short-chain 1-alkenes are the main products of the reaction and need to be well reproduced by the model.The average relative deviations were 14%, 9%, 7% and 8% for ethene, propene, 1-butene and 1-pentene, respectively.For ethene, there seems to be a similar problem as for CH 4 , which indicates a systematic problem within the model for high values.The experimental data suggest some secondary hydrogenation of ethene which is not implemented in the model.This may explain the deviation partly.For propene, 1-butene and 1-pentene, no systematic error can be observed from the parity plots and most data points are within the 10% parity lines.
It is noteworthy that the model does not predict a CO selectivity of 1 for X CO 2 → 0. Riedel et al. [9] also observed this in their experiments.To account for this, they introduced a direct formation pathway for hydrocarbons from CO 2 in their kinetic model.In the detailed model, however, the catalyst surface is already covered with alkyl-species at the reactor entrance at steady state.Thus, hydrocarbons can be directly formed for X CO 2 → 0. Additionally, the model does not contain any information of the RWGS equilibrium, which strongly influences the reaction [14].The detailed model implicitly contains this information through the ratios of the equilibrium/kinetic parameters of the adsorption steps of H 2 , CO 2 , CO and H 2 O.

Influence of Process Conditions on Product Distribution
The influence of process conditions (GHSV, T, H 2 /CO 2 ) on the experimental product distribution in comparison to simulation results is shown in Figures 7-10.
The space velocity (or rather conversion) has a contradictory effect.At high space velocities, the average chain length increases with an increasing residence time (see Figure 7, note that we could not experimentally extract a wax sample for high space velocities).At low space velocities, though, the average chain length decreases with an increasing residence time (see Figure 8).This effect is likely the result of strongly varying CO partial pressures (and thus chain growth monomer concentrations, see Figure 3).CO 2 -FTS is a consecutive reaction of RWGS and FTS.At high space velocities, the CO partial pressure is low as the formation is kinetically limited by the RWGS.The CO partial pressure exhibits a maximum and decreases afterwards due to equilibrium constraints of the RWGS [14].This also leads to a maximum for the average chain length.However, for relevant conversion levels (X CO 2 > 25%), the effect is not very pronounced.For a tubular reactor, the output is always an integral composition of the whole length, which smooths local differences.In a CSTR type reactor, this effect could be likely larger.The 1-alkene and n-alkane fractions seem to be unaffected by the space velocity.Both observations can be well reproduced by the model.
For the temperature and inlet reactant ratio, slight trends could be observed.To avoid a bias due to residence time (or rather conversion) effects, the influence of these parameters are compared at similar conversion levels of the experiments.The temperature influence is exemplary shown at a CO 2 conversion of ≈36% and for the inlet reactant ratio it is shown at a H 2 conversion level of ≈38% (the CO 2 conversion differs significantly upon varying the H 2 /CO 2 ratio).(c) n-alkane fraction vs. carbon number; data points correspond to experimental values and lines to simulation results; conditions: T = 300 °C, p = 10 bar, (H 2 /CO 2 ) in = 3, (I) GHSV = 900 mL N h −1 g −1 (X CO 2 = 41%), (II) GHSV = 7200 mL N h −1 g −1 (X CO 2 = 27%), catalyst: Fe/K@γ-Al 2 O 3 .
A lower temperature led to a slightly lower 1-alkene fraction and a slight increase of the average carbon number (see Figure 9).The model is not able to reproduce these effects properly.It has to be emphasized, though, that the temperature influence is little more than a slight trend.Within the experimental range studied, the model will still provide a reasonable prediction of the product distribution.However, we are aware of the unfavorable numeric structure of the applied Arrhenius expression (Equation ( 2)) and the awkward coupling of the termination reactions with the hydroxyl surface coverage, which is highly temperature dependent (see Section 3.1).Future model developments explicitly need to address these points.
For a higher H 2 /CO 2 inlet ratio, a lighter product and a lower 1-alkene fraction were observed (see Figure 10).A higher H 2 partial pressure leads to a more hydrogenated and shorter product as chain termination via hydrogenation is favored.In this case, the model correctly reproduces the experimental trends.However, considering the relatively drastic increase of H 2 /CO 2 from 2 to 4, the changes in the product distribution were rather marginal.
In general, the product distribution was only moderately affected within the range of conditions studied and the model can provide a reasonable estimate of the distribution within this range.The deviations for >C 12 of the individual hydrocarbon classes as well as the implementation of the temperature dependencies leave room for further improvements.The observed trends are in agreement with literature for the classic CO-based FTS [49,50].It is known that the product spectrum of Fe based FTS catalysts is only moderately sensitive to changes in process conditions [49,51].It is thus not surprising that the same is also observed for the conditions of CO 2 -FTS.Analogue results have been reported by Panzone et al. [13] for a similar catalyst.

Kinetic Experiments and Product Analysis
Detailed descriptions of the kinetic experiments, the catalyst preparation, and the product analysis are provided in the Supporting Information and have also been published in a previous work [14].Here, only a short summary of the applied procedures is provided.Kinetic experiments were conducted in a lab-scale tubular reactor with 0.5 g or 2 g of an in-house prepared potassium doped alumina supported iron catalyst (Fe/K@γ-Al 2 O 3 , mass ratio: 15/5.25/100),diluted with silicon carbide.Fe/K@γ-Al 2 O 3 was first proposed by Choi et al. [52] and was chosen as a model catalyst because of its proven applicability for CO 2 -FTS.The catalyst was activated in-situ with a combined reduction/carburization method, adapted from Landau et al. [53].In total, 31 different experimental conditions were investigated in the following ranges: p = 10 bar, T = 280-320 °C, GHSV = 900-120,000 mL N h −1 g −1 and H 2 /CO 2 molar inlet ratios of 2-4 (see Table S3).The product distribution was determined via GC-FID analysis of the obtained gas and wax phases.For experiments with a GHSV ≥ 12,000 mL N h −1 g −1 (12 out of 31), we did not obtain a wax sample and only considered the gas analysis up to C 8 .Due to the complexity of the obtained product spectrum with several hundred different species, we only allocated 1-alkenes and n-alkanes.All other hydrocarbons with the same carbon number were lumped into the group other, which mainly consisted of iso-alkenes (additionally linear alkenes with internal double bond, iso-alkanes, aromatic components, cyclo-alkanes and oxygenates).These species were modeled as single methyl-branched iso-alkenes.It has to be pronounced that short-chain oxygenates could not be measured and are thus not included in the product spectrum.

Data Analysis and Definitions
Space velocities are given with respect to catalyst mass m cat as the gas hourly space velocity (GHSV) of the inlet volumetric flow rate Vin,N at normal conditions (0 °C and 1 atm).
The conversion X of species i is given as the relative amount of converted reactants, where ṅi denotes the molar flow rate of reactant i: The selectivitiy S of species i is given on a carbon basis (molar flow rate of species multiplied with its carbon number n i ) with respect to the converted amount of carbon dioxide: For the ASF plots (usually logarithmic mole fraction vs. carbon number), we adopted a modified calculation approach.To facilitate the comparison of different experiments, as well as experimental and simulation results, we normalized the molar hydrocarbon flow rates with the amount of carbon that has been converted to hydrocarbons instead of the sum of all hydrocarbons.This measure is referred to as the modified mole fraction x i,mod and is calculated as follows: x i,mod = ṅi,out ṅCO 2 ,in − ṅCO 2 ,out − ṅCO,out (10)

Conclusions and Outlook
A detailed kinetic model for CO 2 -FTS on a supported iron catalyst was developed.It is based on a direct CO 2 dissociation mechanism for the RWGS step and the alkyl mechanism with H-assisted CO dissociation and chain-length-dependent kinetic parameters for the FTS step.The model allows a reliable prediction of the reactant consumption, the overall hydrocarbon distribution, and the fractions of hydrocarbons groups (1-alkenes, n-alkanes, and iso-alkenes) up to C 12 within the experimental range studied.However, the model suffers from the assumption of identical catalyst sites for the RWGS and FTS.This led to the inclusion of an empirical damping factor, which should be avoided in further development steps.
Further model improvements may be the inclusion of oxygenates and secondary hydrogenation reactions.Our currently available experimental data does not allow this though.The application of a PFR type reactor for the kinetic experiments may not have been an optimal choice.The integral character leads to the superposition of different reaction conditions that can hardly be deconvolved for mechanistic conclusions.It is desirable to conduct more detailed mechanistic investigations in a CSTR-type reactor (e.g., a spinning basket reactor), especially regarding the secondary hydrogenation of 1-alkenes.
The complexity of the model makes it probably difficult to implement it into a commercial flowsheet simulation software.To simplify the usage, it might be possible to omit the modeling of the catalyst surface coverage and convert our implicit model into an explicit LHHW-type model.Considering other detailed kinetic models [26,30], this should be possible with a few assumptions.
Nevertheless, we are convinced that it is important to model the reactant consumption and product distribution simultaneously for a realistic reactor under recycle conditions.Recycle operation would be necessary in a technical application to achieve sufficiently high conversions.It might additionally affect the product distribution as large fractions of shortchain hydrocarbons would also be recycled into the reactor.We are currently assessing these issues in a bench scale system and will report on the results in a future publication.
Supplementary Materials: Supplementary Materials are available online at https://www.mdpi.com/article/10.3390/catal12060630/s1providing detailed information about the applied catalyst (preparation and characterization) and kinetic experiments (setup, procedures, product analysis, and results).Table S1: Detailed preparation procedure for Fe/K@γ-Al 2 O 3 .Table S2: Key properties of the investigated catalyst determined via N 2 -physisorption and inductively coupled plasma optical emission spectrometry.Table S3: Conditions and key results of the experiments used in the kinetic study.Figure S1: Electron probe microanalysis map of Fe/K@γ-Al 2 O 3 with K and Fe distribution for an exemplary particle.
Figure S2: (a) Simplified flow scheme of the lab-scale setup (reprinted with permission from Brübach et al.Copyright 2022 American Chemical Society); PI: pressure indicator, TI: temperature indicator, TIC: temperature indicator controller, BPR: back pressure regulator; (b) Photograph of the lab-scale setup (date: 18 October 2021).

Table 4 .
Parameter estimates of the proposed kinetic model determined via weighted least squares optimization.