Systematic kinetic study of magnesium production using magnesium oxide and carbonic materials at different temperatures

The main goal of this study was to determine the industrially best reductant for reduction of magnesium oxide to magnesium with wood charcoal and petroleum coke (petcoke) each in molar ratio 1:1 and 1:2 (oxidant:reductant) at high temperatures. In this study, a new and reliable combination of mathematical modeling and discrete numerical optimization theory by presenting 18 “mathematical filters” not relying only on statistical quantities of fitting (contrary to many similar researches) was introduced. The purpose of these filters was the determination of correct kinetic equation and therefore, the corresponding rate coefficient from among 18 equations most used at present in the challenging field of solid state chemical kinetics. With assistance of a new and fundamental mathematical function and the obtained values of rate coefficients, the function of rate coefficient in temperature was attained. The activation energy was then calculated as a function of temperature using the general definition of activation energy and the determined function for rate coefficient. The comparison between different reducing agents in the different conditions and with relevant previous study was accomplished to determine the best reducing agent from industry standpoint. Also, the areas under experimental data were calculated numerically and utilized for method validation and comparison. It turned out finally that relying only on fitting quantities in the solid state chemical kinetics can readily lead to wrong conclusions about the correct kinetic equation and about the most suitable reducing agent. It is obvious that the erroneous calculations and wrong decisions in the laboratory scale become significant and paramount in industry and this reveals the significance of rigorous mathematical analysis.


Introduction
The magnesium (Mg) element is one of metals in our universe with some interesting properties [1]. One significant feature of magnesium is the ability of it to produce much energy in combustion reaction. The considerable point in burning reaction of magnesium is that it does not produce carbon dioxide, a greenhouse gas. Therefore, it is clear that more attention to this metal can solve some fundamental problems of a society. If magnesium is produced through a clean route, such as using renewable energies, the process will be even "greener". In fact, because of many problems such as greenhouse effect resulting from fossil fuels and also, the decrease of these fuels in the future, a striking endeavor has started for producing clean energies from renewable resources [2][3][4][5][6][7]. In general, renewable energies are of different types and some of them are solar energy, wind energy, tidal energy, geothermal energy, and biomass energy, for example. Among various types of renewable energies, solar energy is an attractive subject for many researchers. It can be said that the human being noticed the importance of solar energy a little late but nonetheless, the subject of solar energy has seen a substantial amount of research projects in various fields [8][9][10][11][12][13][14][15][16][17]. It is obvious that the combination of magnesium production and utilizing renewable energies can be very beneficial for human being in this era.
In parallel, the clean and cheap fuels have been always like dreams for all people and especially for scientists, engineers and even, governors. One of the useful and suitable fuels is the hydrogen gas [18]. Therefore, the combination of solar energy and hydrogen production is of fundamental importance for improving industrial processes in order to deal with many difficulties and problems in a society. One of the substances suitable for production of hydrogen is metallic magnesium. Hence, the production of this metal from magnesium oxide (MgO) can be an interesting goal for research and development. Because of high temperatures needed for cleavage of stiff bonds in the crystal of MgO, a carbon-based reducing agent like carbon itself or methane is usually necessary. Therefore, the high energy stored in the concentrated solar radiation and the chemical energy in the carbon-based material can cooperate to overcome the powerful structure of MgO and to release magnesium in the pure metallic form. The metal produced can create hydrogen gas readily through reaction with water vapor and the hydrogen gas released can be used in the fuel cells for production of electric power. It must be noted that the reaction of magnesium with water vapor is very fast and therefore, this reaction must be performed in a controlled manner and with considering all safety issues.
Additionally, the usage of magnesium itself as a direct fuel has been suggested [19][20][21]. This novel idea is supported from two viewpoints: first, the heat of oxidation of magnesium is about ten times that of hydrogen and second, magnesium available in the seawater is enough for 300,000 years of humankind [19]. The direct usage of magnesium can be made in two modes, reaction with oxygen or reaction with water vapor. When metallic magnesium is mixed with oxygen or water vapor, it produces heat that can be used as a source of energy. The reaction of magnesium with water vapor has the benefit of producing hydrogen gas which can be used to provide energy and the magnesium oxide produced can then be converted back into magnesium using, for example, the solar laser and this forms a useful cycle [19]; the water vapor produced from reaction of hydrogen with oxygen in a fuel cell can be used in this cycle, too. All these studies as well as other fascinating properties of magnesium such as low weight (interesting for automobile and aerospace industries) and many others [22] prophesy about the important role that this metallic element will play in the future of our world. These facts reveal the high potentiality of such interesting fields and the importance and usefulness of research in them. Since chemical thermodynamics and chemical kinetics are the fundamental sciences in exploring chemical processes, it is a logical task to study the chemical reactions occurred in these fields thermodynamically and kinetically.
One of the useful studies carried out on magnesium oxide is the work of Gálvez et al. [23] in which the role of MgO in a solar thermochemical cycle was investigated both thermodynamically and kinetically. Also, Rongti et al. [24] explored the reaction of MgO with graphite powder kinetically.
The purpose of present paper is to perform a mathematical treatment of some significant but untreated points in the study of Gálvez et al. [23] and perhaps many others, more rigorously from standpoints of chemical kinetics (Chemistry) and optimization theory (Mathematics) coupled to each other in a large extent; this coupling, for first time, approaches the problem thoroughly with not relying only on curve-fitting quantities and therefore, it requires more consideration for understanding it.

Theoretical treatment in solid state chemical kinetics (SSCK)
The reactants in the solid state do not have the ability of moving through matter and therefore, the concentration concept is undefined for them. Even if a diffusion process controls the mechanism of reaction, this diffusion is slow in the solid state. For example, the surface of reaction can goes through body of matter while the reactants are in front and products in behind. As a result, it can be possible to substitute the concentration concept with fraction-of-reaction quantity [25]. The suggested mechanisms in SSCK are very diverse. The nucleation is a typical mechanism while another important one is the diffusion in one or several directions. Totally, there are more than 29 important equations in SSCK which come from different theories [26].
Irrespective of large number of equations presented, it is sometimes possible to fit satisfactorily the experimental data with more than one kinetic equation. Generally, finding correct kinetic equation and in other words the real mechanism had been one of the most difficult tasks in SSCK. From equations suggested, 18 mathematical equations derived from previous scientific researches and represented in Table A1 in Appendix A in Online Resource (Supplementary Material) are used more than others [25][26][27]. For simplicity, the abbreviations suggested by Sharp et al. [28] were used in Table A1 and throughout this study.
Because of the nature of reactions in the solid state, the considerable amounts for fraction of reaction can be attained only at high temperatures and hence, the temperature can be considered as a powerful catalyst; this is rooted in the very high activation energies involved in SSCK, originated in turn from lack of collisions in the solid state.

Mathematical modeling and raw kinetic data
The inviolable science of Mathematics is the inseparable part of natural sciences and specially, the fundamental and strategic science of Chemistry named as Central Science  [29]. Like many cases in chemical kinetics, the starting point is the respective differential equation: Before determination of quantities in Eq.
(1), it should be said that according to the study of Gálvez et al. [23], in first step of a solar thermochemical cycle, MgO is reduced to Mg with assistance of carbonic reducing agent (charcoal or petcoke) and thermal energy from solar radiation and in second step, Mg is converted back to MgO through reaction with water vapor and the hydrogen gas is produced as the desired product in Ref. [23]; the focus of our study was on the first step of this solar thermochemical cycle (Eq. (10)). In Eq. (1), a(t) is the fraction of reduction reaction of MgO to Mg proceeded in time t and is dimensionless, k is the rate coefficient with dimension of time −1 , f(a) is a mathematical function taken from Table A1 and without dimension, m(t) is the remaining mass of magnesium oxide in time t and n(t), the number of moles of remaining magnesium oxide in time t. The numerical values of a have been obtained by TGA (thermogravimetric analysis) instrument (TG, Netzsch STA 409 CD) [23]; in fact, instead of solar radiation, Gálvez et al. [23] used the thermal energy produced by thermogravimetric instrument as the heat required for reaction. The basic principle of a thermogravimeter is the measurement of mass change over time as temperature changes (dynamic runs) or is fixed (isothermal runs). For kinetic studies, the isothermal runs are more favored and it is possible to guess kinetic parameters and mechanism from performance of a chemical reaction at different temperatures with temperature treated as a parameter (that is, a fixed temperature for an individual run but different relative to temperatures of other runs). The raw kinetic data used in our study were taken from isothermal runs of Ref. [23]. Usually, other instruments and techniques such as Fourier-transform infrared spectroscopy (FTIR), mass spectrometry, and gas chromatography (GC) are coupled to thermogravimetric operations for analytical purposes; in Ref. [23], GC method was used in combination with thermogravimetric operations. In general, Ref. [23] is a very comprehensive study with various experimental techniques used but only some parts of it is pertinent to the subject of our study.
After rearranging the differential equation in Eq. (1) and taking integral from both sides, the function g(a) in Table A1 can be created: In Ref. [23], the D3 model was used for fitting of charcoal curves and the R3 model was utilized for petcoke curves; our paper, however, illustrates the determination of correct kinetic equations more systematically and more mathematically and with assistance of the behavior of experimental kinetic data.

Calculation of rate coefficient (k)
The rate coefficient k is one of the most important quantities in chemical kinetics and in most cases, including those in this study, has a fixed value by fixing temperature. What should be done is to obtain the rate coefficients at certain temperatures from correct kinetic equations in order for calculating the activation energies (another important quantity) and to compare these energies for selecting the most powerful reducing agent for reduction of metal oxide. The correct kinetic equation is unknown initially and the strategy is to calculate the rate coefficients for 18 equations and the detection of correct kinetic equation by some "mathematical filters" (introduced shortly) which are exerted on the calculated values of k. For calculation of k for any 1 of 18 equations in a fixed temperature, 4 methods can be considered in initial sight but only one optimum method must be selected. After selection of best method, 3 other methods are eliminated from discussions. Then, various mathematical tests, all relevant to selected method and introduced shortly, are presented for selection of best kinetic equation from among 18 mathematical equations; the selected rate coefficient comes from this best equation. It is important not to confuse things; in short, at first, only 1 method is selected from among 4 methods and then, all 18 equations are examined by various mathematical tests based on the selected method to give one selected equation and therefore, one selected rate coefficient. The 4 methods are:

Method 1
To extract the mathematical functions a(t; k) from functions in Table A1 and to fit the experimental data of a versus t to the functions a(t; k) for finding the rate coefficient k.

Method 2
The same as method 1 except that in this method, t is in the role of dependent variable and a in the role of independent one. Since the purpose is to compare 18 equations together and the shift of variables is done in all 18 equations, this shift is of no problem from standpoint of mathematical logics.
In this method, the curve of experimental data of a versus t is plotted and then, the quantity da/dt is calculated in every point of the curve. The calculation of da/dt is possible from two ways. First is to use the numerical algorithms for calculation of da/dt but this method can have appreciable error. Second is to fit as exact as possible the experimental data to a mathematical function compatible to the boundary conditions of experiment. It is not necessary that this function to be from those in Table A1 because this function is only an auxiliary one; in fact, the fundamental properties of this function must be the exactness in fitting and the compatibility to the boundary conditions of experiment. Considering a = 0 in t = 0 and a = 1 in t = ∞ and the particular shape of the experimental data of a versus t, the correct selection and the best choice is a function of the general form given in Eq. (3): where s.t. denotes "subject to" or "such that". After determination of function a(t), or more precise a(t; b, c), the derivative of this function is calculated and by inserting the experimental data of t in it, the quantity da/dt is calculated in every point of the experimental curve. The quantity f(a) in these points can be calculated by means of experimental data of a and Table A1. Considering Eq. (1), the value of k can be calculated through two routes. In the first route, the curve of da/dt versus f(a) is plotted and the slope of linear fitting function gives the value of k. In the second route, the value of k is determined from Eq. (1) for each datum of a and the average of obtained values of k gives the final value for k. It turned out for our data that two routes give essentially identical value for k; the second route was adopted in this study.

Method 4
The same as method 3 except that in this method, Eq. (2) plays the role. Now, the advantage(s) and disadvantage(s) of these 4 methods are expressed in order to select only one optimal method and to throw away other three methods. About method 1, it should be said that it is not possible to extract a(t; k) for all functions in Table A1. Since the goal is to compare all 18 equations together, the method 1 is not an interesting one. Also, this method relies only on agreement between experimental data and fitting curve and uses only statistical quantities like RSQ (coefficient of determination or R-square or R 2 ), RMSE (root mean squared error), and SSE (sum of squares of error) [30] for finding k. One of the permanent problems in chemical kinetics and especially in SSCK has been the satisfactory agreement between data and fitting curve while this fitting function in fact does not represent the real mechanism of reaction but the apparent agreement is reflected well in statistical quantities of fitting [30]. Additionally, the values of a are small and application of them in statistical quantities, which use them several times, introduces some error. For method 2, it should be said that this way does not have the fault of method 1 for extracting the function t(a; k) from Table A1. However, the other faults of method 1 is existent for method 2. About methods 3 and 4, it should be said that the method 4 is more exact because the mathematical form of function g(a) is definite and analytically obtained and extra steps of method 3 and unwanted resulted errors do not exist. As a result, the method 4 can be taken as the fundamental and principal method and various mathematical tests are derived from it (method 4) for examination of 18 equations; these tests determine the correct kinetic equation and therefore, the correct rate coefficient in an especial fixed temperature. Appendix B in Online Resource describes the usage of method 4 more mathematically and rigorously, but for those uninterested in mathematical details, Appendix B * in Online Resource is essentially the same as Appendix B but without mathematical equations; only the results are mentioned in paper.
It is clear that the more powerful and flexible the mathematical programming used, the more precise the huge amounts of calculations and because of this key point, the calculations in this study were performed by several mathematical programs written in the powerful and fast languages of MATLAB and C and combination of them known as MEX Functions.

Calculation of activation energy (E a )
After detecting the correct kinetic equation and finding the corresponding rate coefficient, it should be noted that all actions performed (mentioned more rigorously in Appendix B in Online Resource) are for reaction at one certain temperature. Therefore, all those long steps must be accomplished also for other temperatures. After obtaining the values of k for all temperatures of reaction, the activation energy of reaction can be calculated. In some studies, the well-known equation of Arrhenius has been used [23,[31][32][33]: T, A, and R are thermodynamic temperature, pre-exponential factor and gases constant, respectively. In general, the calculation of activation energy by non-linear regression from non-linear Eq. (4) is more exact relative to linear regression of logarithmic form of Eq. (4); the reason has to do with the statistical weights and residuals of data and is a detailed statistical discussion. The Arrhenius equation is a form of Hood equation [34]: where B is a constant. The questionable point in Eq. (4) is the independency of activation energy from temperature. Since on one hand, the reactions in SSCK are particular reactions and occurred at high temperatures and on the other hand, it is completely possible to obtain the values of k at different temperatures from different kinetic equations (because of different mechanisms), it is a logical action to take into account the possibility of change of activation energy with temperature. Thus, the original form of Arrhenius equation can be a more exact definition of activation energy [35]: Using Eq. (6) analytically is more exact relative to numerical usage of it and this requires the function k(T) to be introduced. The function k(T) can be found by curve fitting of obtained values of k at different temperatures but the general mathematical form of function k(T) should be selected carefully. One of the proposed forms is as [34][35][36].
in which, s is usually a small positive integer related to statistical thermodynamic discussions [35]. Equation (7) is a generalized form of Arrhenius equation and acts more exact in some cases. It will be more exact to take E a /R in Eq. (7) as a constant like C and to utilize Eq. (6) to obtain activation energy as E a = CR + sRT. Considering these points, the perfect mathematical form of Eq. (8) can be a suitable choice and was utilized in this study: Since Eq. (8) contains 4 parameters, the number of data points (temperatures) needed for determination of 4 parameters from fitting should be at least 4 and therefore, for cases like molar ratio 1:2 with only 3 data points (introduced shortly), D was taken as 1. Combining Eq. (6) and Eq. (8), the corresponding activation energy can be calculated as Reactions of MgO with carbonic reductants The general form of reaction of MgO with carbonic reductants at high temperatures can be shown as Molar ratio of MgO to reductant in all reactions performed was 1:1 or 1:2 [23]. The trend of calculations is described in details for one reaction in a certain molar ratio and only the results are presented for other cases.  Table 1 is an illustration and the quantities with number 1 are for curves of g(a) versus t, the quantities with number 2 for curves of k versus t and the quantities with number 3 for curves of t versus a.
As can be seen in Table 1, from first 15 tests, only 3 tests confirm a mechanism other than D4. Considering the power of tests, it is obvious that the most logical selection for correct kinetic equation is D4. Therefore, the rate coefficient is taken from equation D4 and two other equations (D3 and F3) are deleted for this case. The values of k for Table 1 The results of applying different tests to the different kinetic equations for determination of average k for reaction of MgO with charcoal in molar ratio 1:1 and at temperature 1723. 15

15)
Range ✓  Table 2. It is useful to mention this point that F mechanisms have been borrowed from other branches of chemical kinetics in which, contrary to SSCK, the concepts of concentration and collision are meaningful. Therefore, F mechanisms are in a lower grade compared to other mechanisms of SSCK. The defeat in 13 tests and the better fitting by respective t(a; k) reflected in RSQ (3), RMSE (3), and SSE (3) in Table 1 for F3 equation is in fact a clear example of wrong conclusions based on only the statistical quantities of fitting like RSQ; as can be seen in Table 2, the value of k for F3 model is about two orders of magnitude larger than that of D4 model and this affects the calculation of activation energy significantly. In fact, the erroneous calculation of activation energy can be very detrimental with respect to some critical decisions in industry such as identification of the bottleneck of process. Figure 1a illustrates the experimental data as well as final selected equations for charcoal in molar ratio 1:1 at different temperatures; except for 1823.15 K, all other temperatures follow D4 equation. Fig. 1b shows the experimental data fitted by Eq. (3) for these cases; this equation lacks the fundamental discussions about mechanism and emphasizes only on fitting and hence, cannot be used for calculation of activation energy; Appendix C in Online Resource shows an interesting usage of Eq. (3) regarding the estimation of 90% of reaction time. Table 3 reports the abbreviations of selected kinetic equations and also the corresponding optimized values of k for reactions of MgO with carbonic reductants in different molar ratios at different temperatures. In the study of Gálvez et al. [23], the D3 equation was used for all charcoal cases and the R3 equation for all petcoke cases; RSQ and RMSE were the important criteria for selection of these models in Ref. [23]. In general, the unsatisfactory fitting by some selected kinetic equations results from the lack of enough evolution of theories in SSCK and this problem is considered as the largest challenge in SSCK. In fact, it can be said that among various fields of chemical kinetics, SSCK has the most amount of ambiguity and obscurity and the elucidation of its mysteries requires the more attention, time and cost, relative to what is in action at present.

18) SSE (3) ✓
The mechanisms of chemical reactions have been always a very challenging subject in Chemistry, and this difficulty is even more severe in the solid state. As can be seen in Table 3, except for charcoal in molar ratio 1:1 and at temperature 1823.15 K that obeys D3 mechanism, other cases follow D4 mechanism. In all, it is inferred that the diffusion is the rate-determining step in the processes involved in this study. Since the boiling temperature of Mg metal is about 1363 K and considering the higher temperatures studied in this paper, the produced magnesium in Eq. (10) will be as gaseous and hence, it can be said that the products of reaction exit from reaction zone and thus, the rate-determining step of diffusion concerns probably only the reactants. The significant point is that the melting point of a solid has inverse relation with diffusion [37] and hence, the high melting points of reactants in Eq. (10) are compatible with the conclusion that diffusion is the rate-determining step. Since D4 is an extension to D3 model, the numerical values of rate coefficients for D3 and D4 models in the identical conditions are not very different and this fact was observed in this study frequently; Table 2 contains one instance. Therefore, it can be said that the D4 equation can be taken as the fundamental kinetic model for describing the kinetic behavior of reaction of MgO with charcoal and with petcoke, each in molar ratios 1:1 and 1:2 and at different temperatures. The derivation of mathematical equation for D4 mechanism using first Fick law is discussed in detail in Ref. [27].
The subtle and significant point in Table 3 is the very small values of rate coefficients for reactions in the solid state and especially, for reaction of MgO with carbonic materials according to Eq. (10), in spite of very large temperatures involved. In fact, it can be said that if Eq. (10) is a part of a chemical process in a facility, it will be a limiting step in entire process. This limiting property is highlighted when comparison with the reactions, for example, in solution; for instance, the reaction of C 2 H 5 I with OH¯in  ethanol has a rate coefficient of 0.119 L mol −1 s −1 at 363.8 K [35], where L denotes liter. This limiting nature of reactions in the solid state will be reflected at high values for activation energy, discussed shortly.

Activation energies at different conditions
The calculation of activation energies requires a good fitting of obtained values of k. Figure 2 shows the final values of k from Table 3 as well as curves of k(T) obtained from fitting by Eq. (8); in figures of this study, C11 is applied to charcoal 1:1, C12 to charcoal 1:2, P11 to petcoke 1:1 and P12 to petcoke 1:2. The parameters obtained from fitting and the relevant statistical quantities of fitting [30] are reported in Table 4. The interesting point in Table 4 is that the values of parameter D for charcoal 1:1 and petcoke 1:1 are practically 1; as pointed out below Eq. (8), because of 5 values of k for 5 temperatures for molar ratio 1:1 for both charcoal and petcoke, D can be set as fourth parameter in 1:1 cases (unlike molar ratio 1:2 with D equal to 1 from beginning). Because of low number of points in the case of molar ratio 1:2, the corresponding fittings in Fig. 2 and the resulted activation energies are not very reliable; in fact, because of equality of number of parameters and number of experimental data points in molar ratio 1:2, the quantity of RMSE is infinity, as can be seen in Table 4 and verifiable from Eq. (B7) in Appendix B in Online Resource. In general, the more the number of experimental data, the more precise the calculations. The values of E a at different temperatures, the values of SD, RSD and range for E a , the mean value of E a , E a obtained from non-linear regression of Arrhenius equation, E a obtained from linear regression of Arrhenius equation (using logarithmic linear form as in Ref. [23]) and E a from Ref. [23] are all holus-bolus collected in Table 5. The reasons of equal E a for each of molar ratio 1:2 are the essentially zero value of B and exact value of 1 for D that according to Eq. (9), cause this property. In Ref. [23], E a was obtained only for 1:1 molar ratio for both reducing agents and 1: 2 molar ratio was not treated in this regard; according to Table 5, contrary to Ref. [23], E a for 1:1 molar ratio for petcoke case is larger than that of charcoal case; introduced  shortly, the quantitative validation of method using experimental data will clarify that our treatment gives more real values for E a , relative to treatment of Re [23].. It is also noteworthy that in Table 5, the values of E a for an individual carbonic reductant but in different molar ratios are different; more has been said about this difference in Appendix D in Online Resource. In the study of Rongti et al. [24], only partly similar to our work, they examined the reaction between MgO and graphite powders in a dynamic (non-isothermal) manner in the temperature range 293 to 1973 K and obtained the average activation energy equal to 208.29 kJ mol −1 for 0.1 < a < 0.25 and 374.13 kJ mol −1 for 0.25 < a < 0.5; these values were found for molar ratio 1:2 (MgO:C); they also calculated the activation energy using transition state theory (TST) as 470 kJ mol −1 and introduced a chemical mechanism but did not discuss the kinetic equations. In a study by Shu et al. [38], the isothermal reduction of Scheelite (CaWO 4 ) by silicon (Si) from 1423 to 1523 K was accomplished from a kinetic viewpoint and it turned out that both D3 and D4 models can describe the kinetics of reaction well; Arrhenius equation (in logarithmic linear form) was used for calculation of E a and the values of 379.93 and 387.16 kJ mol −1 were obtained for D3 and D4 models, respectively. One point clearly observable in Table 5 is the high values of activation energy for reactions in the solid state, relative to those in other phases; for example, the reaction of C 2 H 5 I with OH¯in ethanol has E a equal to 90.3744 kJ mol −1 , obtained from Arrhenius equation in logarithmic linear form [35]. In fact, the highness for values of E a in Table  5 is rooted in the lowness of k values in Table 3. Consequently, it can be inferred that the reaction of MgO with carbonic materials (Eq. (10)) can be considered as the limiting step in a process.

Quantitative validation of method using experimental data
Because of significance of activation energy in the design of chemical processes in industry and for suitable treatment of slow reactions (chemical bottlenecks), the comparison of E a values is of critical importance. One of the important results obtained from Fig. 2 and Table 5 is that, contrary to Ref. [23], the values of average E a in both molar ratios are larger for petcoke than those for charcoal. In this section, the quantitative validation of adopted method in present paper is introduced to show the validity of conclusions about activation energy; before that, it will be appropriate to mention a related subtle point regarding change of rate coefficient with temperature.
In general, for finding the amount of sensitivity of a reaction to temperature change, the quantity of dk/dT is a good choice. Rewriting Eq. (6) as Eq. (11), it is obvious that both k and E a appear as effective factors but usually, it is the quantity of dlnk/dT that is used in this regard. The point is that dlnk = dk/k is, in fact, the limit of Δk/k when ΔT → 0 and the fraction of Δk/k is a more suitable quantity for investigation of sensitivity compared to Δk. Therefore, Eq. (6) is more appropriate compared to Eq. (11) for evaluation of sensitivity to temperature. Rearranging Eq. (6) as Eq. (12), and it is obvious that k does not appear as a contributing factor in right hand side of Eq. (12). Therefore, according to Eq. (12), the larger the activation energy for a reaction, the higher the sensitivity of that reaction to temperature; this sensitivity approaches zero when T → ∞ and approaches infinity when T → 0. Considering Eq. (12) and Table 5, petcoke with molar ratio 1:1 will have the most sensitivity to temperature. Figure 3a, b, practically identical but drawn in the different styles, confirms this conclusion. If the number of data would be more and equal for all cases, the sensitivity of petcoke to temperature would flaunt even more. In Fig. 2, it is also obvious that the numerical values of k and therefore, the corresponding fitting curve for petcoke 1:1 are lower than those of charcoal 1:1 at not very high temperatures but are higher at very high temperatures; the same situation is observed for numerical values of k for petcoke 1:2 relative to charcoal 1:2; the reason of not observing this trend for corresponding fitting curve of petcoke 1:2 relative to charcoal 1:2 is the low number of data points and high curvature of fitting function required for petcoke 1:2, leading to disability of fitting method in giving a good fit.
More quantitatively and in principle, the areas under experimental data in Fig. 3a, b can be calculated numerically and compared; this comparison would be more accurate if initial and final times would be the same for all experimental runs. Notwithstanding, for an approximate exploration, these areas were calculated for charcoal 1:1 and petcoke 1:1 at temperature interval T = 1773.15 K to T = 1823.15 K with the initial point taken in 40 min and the final point in 348 min; the reason of such selections was the relatively good coverage of experimental data, as can be seen in Fig. 3a, b. Figure 4a shows this quantitative comparison of AUC, the area under curve, obtained by numerical integration using a program written in MATLAB; in accordance with Eq. (12) and following the "overtaking" behavior observed in Figs. 2 and 3a, b, petcoke 1:1 is more sensitive to temperature change and this sensitivity is especially large in lower temperatures. Also, by noting Fig. 3a, b, it can be said that the temperature T = 1773.15 K is a "jumping temperature" for both charcoal 1:1 and petcoke 1:1 that "excites" the reactants for occurrence of chemical reaction; this behavior does not mean that the reaction must be performed at lower temperatures and the point is that the lower the temperature, the more the sensitivity to temperature change. Figure 4b shows this behavior better by plotting where AUC 1 and AUC 2 are the areas corresponding to the first and second temperatures in each subinterval in Fig. 4b, respectively; for example, for first temperature interval in Fig. 4b, ΔAUC = AUC(T = 1798.15 K) -AUC(T = 1773.15K). Also, the elapse of enough time for a reaction is an important factor in achieving the desired result from that reaction. Therefore, totally, from industry standpoint and in a reasonable amount of reaction time, charcoal with excess molar ratio is the recommended selection for reducing agent at not very high temperatures. At very high temperatures, petcoke with extra molar ratio is a suitable choice and even perhaps, more suitable than charcoal with extra molar ratio, especially if it is possible to provide long time for reaction. It must be noted that the difference in the behavior of carbonic reducing agents reported in this study is itself a detailed scientific discussion and is relevant to the different microstructural properties of different forms of carbon [39,40].
In end, it can be said that, like many other reactions in the solid state, diffusion is the dominant process controlling the rate of reactions in this study, especially that the reactants have high melting points. Also, it turned out that the D4 kinetic equation is a powerful mathematical model for kinetic treatment of reactions in the solid state including reactions of this study; Appendices E and F in Online Resource discuss D4 model more thoroughly and Refs [41][42][43][44][45][46][47][48][49][50][51][52][53][54]. are related to Appendix F which addresses the usefulness of D4 equation in analyzing various reactions in Chemistry.

Conclusion
In this study, a new and detailed combination of mathematical modeling and discrete numerical optimization theory was used by presenting 18 mathematical filters not relying only on statistical quantities of fitting. The purpose of these filters was the determination of correct kinetic equation and therefore, the corresponding rate coefficient from among 18 equations most used at present in solid state chemical kinetics. Then, the obtained values of rate coefficient were fitted by a new fundamental mathematical function. Using the general definition of activation energy and the determined function for rate coefficients, the activation energy was introduced as a function of temperature. The comparison with previous study was performed and it became clear that, contrary to it, the values of average E a for petcoke are larger than those for charcoal in both molar ratios and besides, the values of E a for an individual carbonic reductant but in different molar ratios are different. Also, the areas under experimental data were calculated and utilized for method validation and comparison. In all, with assistance of experimental data and mathematical discussions, it became clear that for reduction of magnesium oxide to metallic magnesium in a reasonable amount of reaction time, at not very high temperatures, charcoal with excess molar ratio is the most beneficial choice industrially but at very high temperatures, petcoke with extra molar ratio is also suitable and even perhaps, more suitable than charcoal with extra molar ratio, especially if it is possible to provide long time for reaction.
Also, it can be said that, like many other reactions in solid state, diffusion is the dominant process controlling the rate of reactions in this study, especially that the reactants have high melting points. Besides, it turned out that the D4 kinetic equation is a powerful mathematical model for kinetic treatment of reactions in the solid state including reactions of this study.