Modeling study of the anti-knock tendency of substituted phenols as additives : an application of the reaction mechanism generator ( RMG ) †

This work presents kinetic modeling efforts to evaluate the anti-knock tendency of several substituted phenols if used as gasoline additives. They are p-cresol, m-cresol, o-cresol, 2,4-xylenol, 2-ethylphenol, and guaiacol. A detailed kinetic model was constructed to predict the ignition of blends of the phenols in n-butane with the help of reaction mechanism generator (RMG), an open-source software package. The resulting model, which has 1465 species and 27 428 reactions, was validated against literature n-butane ignition data in the low-to-intermediate temperature range. To rank the anti-knock tendency of the additives, engine-like simulations were performed in a closed adiabatic homogenous batch reactor with a volume history derived from the pressure profile of a real research octane number (RON) engine test. The ignition timings of the additive blends were compared to that of primary reference fuels (PRFs) to quantitatively predict the anti-knock ability. The model predictions agree well with experimental determinations of the changes in RON induced by the additives. This study explains the chemical mechanism by which methyl-substituted phenols increase RON, and demonstrates how fundamental chemical kinetics can be used to evaluate practical fuel additive performance.


Introduction
Additives are frequently incorporated into commercial formulations because they empirically improve the performance of the product. It is thought that many fuel additives work by altering the reaction kinetics, but for a newly developed additive, the exact chemical mechanism is usually unknown. Better understanding of the mechanism of action of additives is expected to accelerate the development of better fuels and other products, and to facilitate the introduction of new additives with environmental advantages such as those derived from waste biomass. Historically, the complexity of this kinetics has been a major challenge, but recent advances in computational methods for chemical kinetics make it possible to tackle these real-world problems.
Engine knock, an undesirable combustion phenomena caused by the auto-ignition of the unburned mixture (end gas), significantly limits the efficiency of spark-ignition engines. 1 Therefore, the anti-knock tendency, which is described by octane number (ON), is a good criterion to screen the candidates in the development of additives and drop-in fuels. However the standard octane rating procedures, including the rating of research octane number (RON) 2 and motor octane number (MON), 3 are capital and time intensive. First, the new proposed additives must be synthesized and purified, and then it requires a trained technician to operate a specially made cooperative fuel research (CFR) engine following a complicated empirical rating protocol. 2,3 After all of those efforts, the only result is a single number without revealing any mechanistic information. Although there are some studies trying to correlate the ON with other observables such as infrared spectroscopy, 4 distillation curve, 5 and dielectric spectroscopy, 6 these methods are developed by fitting the existing ON database; no studies have reported a method for accurately predicting the ON of fuel mixtures containing additives not yet synthesized. Even a perfect method for predicting ON would not provide perfect ranking of fuels, since fuels which have the same ON can have very different ignition behavior in modern engines. 7,8 An ideal model would be able to predict the ignition delay of any fuel mixture at any engine-relevant condition, not just at the conditions of the RON or MON test.
In an engine, the end gas is compressed by the piston motion and by the flame propagation around top-dead-center (TDC).
Knocking will occur if the end-gas auto-ignites before being consumed by the flame. [9][10][11] Therefore, the anti-knock tendency is mainly determined by the ignition delay time and the flame speed, which are both controlled by the combustion kinetics of the fuel. Consequently, the anti-knock tendency of an additive can be studied by simulation if (1) its combustion kinetics is known and (2) one can accurately estimate the conditions that the end gas experienced in the CFR test.
The former problem and the need for a model accurate at many reaction conditions can be solved by the reaction mechanism generator (RMG), 12 which is an open-source software package designed to automatically construct kinetic models using a flux-based algorithm for model expansion. As demonstrated in previously published papers, 13-15 a reasonable model can be obtained with minimal manual work, making it feasible to study a series of molecules in a short time. Moreover, with the help of the kinetic model, the chemical origin of the antiknock tendency of specific species can be analyzed, which will benefit the design of fuel additives in the future.
To estimate conditions that the end gas experienced in the CFR test, some studies [16][17][18] simply correlated the ON with the ignition delay in a constant volume batch reactor. However, more rigorous treatment also considers the actual pressure and temperature history experienced by the end gas. Westbrook et al. 19 suggested incorporating a pressure profile from a real engine test into a perfect stirred reactor (PSR) to study the relation between chemical structure and octane sensitivity. Similarly, Badra et al. 18 proposed a variable volume simulation, which incorporates a pressure-derived-volume profile into a homogenous batch reactor, to correlate the ignition timing with the RON and MON of the fuel. A variable volume similar to that suggested by Badra et al. 18 is used in this work to study the antiknock tendency of phenolics.
In the following sections, the above method will be applied to predict the anti-knock tendency of several substituted phenols, which are important products derived from the pyrolysis oil of lignocellulosic-biomass and have the potential to be used as drop-in fuels or fuel additives. [20][21][22] The methodology for model generation is introduced first. The base fuel model is validated with the literature data, and the effect variable volume condition has on base fuel chemistry is discussed. Next, the predicted changes in ignition delay at constant volume adiabatic conditions induced by adding each of the six phenols are presented. Then, the anti-knock performance for each phenol is predicted under variable-volume RON conditions, and the predictions are compared with experimental RON data on these additives blended in gasoline. Finally, the chemical mechanism causing the anti-knock property for each substituted phenol is discussed.
ESI † includes derivations for equations used to calculate volume profiles used in the engine-like simulation, a table of said volume profile, and sensitivity analysis for phenolics not discussed in detail in the main text. It also includes sample RMG input files used to generate the mechanisms, CHEMKIN files of the mechanism used in this work, and the species dictionary of the mechanism.

Detailed kinetic modeling
RMG 12 is an automated reaction mechanism generator using a flux-based model enlargement algorithm. First, it reacts inputted core species at specified conditions to propose possible reactions and products in the model edge. The importance of an edge species is judged by calculating the chemical flux of the core towards this edge species in homogenous isothermal reactors. When a certain chemical flux is exceeded, an edge species (and all associated reactions) are moved to the core to subsequently generate additional edge products and reactions. Then, the simulation is restarted with the revised core. This cycle of adding core species continues until some convergence criterion, usually conversion of a key reactant, is reached. The converged core model is then exported as an RMG generated detailed kinetic model. In the above process, the thermochemistry of the generated species are mainly estimated by the Benson group additivity scheme, 23 which has been extensively validated and used in the thermochemistry estimation 23,25,26 and was stated to have an error of less than 1 kcal mol À1 for most species. 23 The kinetics of the generated reactions are estimated by mapping to rate rule templates or training reactions collected in the database. More details of RMG can be found in the literature. 12,24 The model used in this study was generated by RMG-Py (version 2.0.0) and RMG-database (version 2.0.0) with some additions. Special attention was paid to the thermochemistry estimation of substituted aromatics and the reactivity of the cyclic structures. Non-nearest-neighbor interactions for substituted aromatic compounds 25 and substituted aromatic radicals 26 were implemented into the Benson group additivity scheme to improve its accuracy. Samples of the RMG input files are included in the ESI. † Other than the libraries listed below, no other special treatments of RMG or quantum calculations were employed in this study. This ensures the method proposed in this study can be easily applied to the study of the anti-knock tendency of other compounds.
In this study, n-butane is selected as the base fuel because of the relatively small size and moderate RON of 94, 27 which is similar to that of practical gasolines. 27 The n-butane model is generated under conditions of 1 bar and 650-2000 K. Most of the reaction rate coefficients are estimated using high-P-limit Arrhenius forms in the RMG database, many derived from CBS-QB3 TST calculations. Foundational Fuel Chemistry Model Version 1 28 and the H 2 /O 2 mechanism from Burke et al. 29 were used as a seed mechanism to include pressure dependence for small species. Additionally, a few sub-mechanisms were used as libraries to provide better rate coefficient estimates for selected reactions: the methyl formate mechanism of Dooley et al. 30 provides some of the low temperature chemistry for butane oxidation, rates for phenol decomposition by Brezinsky et al., 31 and calculated rates for cyclopentadienyl pyrolysis by Long et al. 32 The models for the blends of substituted phenols were generated using the aforementioned n-butane inputs as a base. Six individual sub-models were respectively generated for the additive/butane blends, each with 2% mole fraction of the additive. The reactor conditions were 20 bar and 650-2000 K, with equivalence ratio unity. The six sub-models were merged together to obtain a unified model, which will be used in the following analysis. Finally, some species (and corresponding reactions including these species) were manually removed from the final merged mechanism. The removed species included polycyclic species and bi-radicals whose thermochemistry the current version of RMG poorly estimates. There is little flux towards these species, consequently, the removal of these species barely changes the simulation results. The model sizes of the sub-models and the merged model are shown in Table 1. The chemical structures of the important chemicals in this study are shown in Table 2. The final merged model is provided in the ESI. †

Engine-like simulation
The aim of this simulation is evaluating the effect on the reactivity when 2% of the fuel is replaced by the additives. A full engine CFD simulation is infeasible due to the size of the mechanism, but a constant volume adiabatic batch reactor would be too simple considering real engines have highly dynamic pressure and temperature profiles. Therefore, an engine-like simulation is needed to keep the simplicity of the adiabatic batch reactor but include the pressure and temperature profile from the real engine test.
Such engine-like simulation is based on the understanding of a simplified in-cylinder scenario of the RON test. After the intake valve is closed in the RON test, the mixture is compressed to a high pressure, high temperature condition by the piston. Then the spark plug ignites the mixture, initiating a flame from the spark plug that propagates to the cylinder wall. Due to the flame propagation, the in-cylinder pressure continues increasing and compressing the end gas. Knocking occurs if the end gas auto-ignites before being consumed by the flame. Therefore, the critical condition would be the auto-ignition and the completion of fuel consumption by the flame occurring at the same time. It is assumed the laminar flame speed and the specific heat ratio of the tested mixtures are similar because of the similar fuel compositions. Consequently, the environment of the end gas Table 1 The size of the sub-models and the merged model  can be characterized by a pressure profile from a RON test with the following assumptions: (1) Homogenous assumption: the end gas is spatially homogenous.
(2) Adiabatic assumption: the end gas is adiabatically compressed by the piston motion and by the expansion of the burned gases in other parts of the cylinder.
Pressure profiles of RON tests of PRF100 was obtained from a CFR engine at UC Berkeley. The compression ratio is 9.21. The spark timing is 15 crank angle degrees (CAD) before top-deadcenter (BTDC). The intake valve close timing is 34 CAD after bottom-dead-center (ABDC). Therefore, the pressure history from 146 CAD BTDC (i.e. 34 CAD ABDC) to 50 CAD after top-dead-center (ATDC) is used in the following simulation. Fig. 1 shows the pressure profile averaged from 500 cycles.
Based on the adiabatic assumption, the pressure profile was converted to a volume profile 33 using the equations: ð Tc In the above equations, g was calculated for a stoichiometric mixture of iso-octane and varied with temperature, but not composition, over time. This profile was imported into an adiabatic batch reactor in Chemkin-PRO 34 to simulate the behavior of the end gas. In the simulation, the initial pressure was 1.09 bar, which was the experimental pressure at valve closing (À146 CAD BTDC). The initial temperature is not available from the experiment because the standard RON test protocol only ensures the intake air temperature is 325 K in the manifold; the temperature in the cylinder is not measured. Heat transfer from the residual gas and cylinder walls is expected to heat the intake air to a higher temperature around valve closing. Therefore, we tuned the initial temperature at À146 CAD BTDC in the PRF100 simulation to approximately match the simulated end gas's ignition time with the time of the measured maximum pressure in the experiment. The mechanism used in the tuning process is the PRF100 (i.e. iso-octane) mechanism from the Lawrence Livermore National Laboratory (LLNL) 35 and resulted in a tuned initial temperature of 400 K, which is used in all of the enginelike simulations in this study.
The timing of the ignition in the engine-like simulation is related to the RON of the fuel. While the actual compression ratios and pressure traces vary slightly with fuel composition, it is still relatively accurate to use fixed reaction conditions for the study of the additives because the additive only accounts for 2% in the total fuel, and flame speed is known to be relatively insensitive to fuel composition. 36 Later ignition under these fixed conditions corresponds to better anti-knock behavior. Fig. 2 shows the computed ignition timing of PRFs and n-butane versus RON in this variable volume adiabatic batch reactor, using the LLNL PRF model 35 and the merged butane model described in Section 2.1. Ignition delay in these simulations is defined as time when T reaches 1500 K. As expected, the PRFs with lower RONs ignite faster than those with higher RONs. In the range of RON 90-100, the timing of the ignition varies nearly linearly with RON. By interpolating the predicted ignition timing of n-butane with those of PRFs, the RON of n-butane is predicted to be 97.0, which is reasonably close to its real experimental RON of 94, 27 indicating the good consistency of the merged kinetic model constructed here and the LLNL PRF model.  37 with the prediction of the merged model at j = 1 and diluent to oxygen ratio of 3.76. The experiments with ignition delay longer than 5 ms were conducted in a rapid compression machine (RCM). Data points with ignition delay shorter than 5 ms were conducted in a shock tube (ST). The simulations were performed using Cantera 38 in a constant volume adiabatic homogenous batch reactor with nitrogen as a diluent without accounting for heat losses. In this simulation, ignition time was taken to be the time when dP/dt reached its maximum.

View Article Online
The negative temperature coefficient (NTC) phenomena is successfully captured by the simulation, although it occurs at a lower temperature in the predictions than in experiments, likely because the prediction does not account for heat loss. The simulated ignition delays mostly agree within a factor of three with the experimental data, which is reasonably good considering the simplified model, experimental uncertainties and many imperfectly known rate and thermochemistry parameters. The low temperature oxidation of butane is only briefly discussed here, but a more detailed discussion for analogous pathways in propane is given in Merchant et al. 39 The main reaction pathways are outlined in Fig. 4 to reveal the importance of the OH radical in driving the low temperature combustion chemistry. At temperatures below 800 K, the fuel consumption relies heavily on reactions with OH radical and to a much smaller degree HO 2 radical. In the low temperature regime, one OH radical is consumed to create the initial radical of the fuel, and then secondary chemistry can result in formation of up to three more OH radicals, as shown by the auto-catalytic cycle shown with the black arrows in Fig. 4. Above B690 K, this reaction pathway creates less than one OH for each OH inputted because two key reactions shift away from this cycle: QOOH begins to favor HO 2 + alkene formation, and O 2 QOOH begins to decompose back to QOOH + O 2 . From 690 K o T o 800 K, the system loses more OH than it produces which is the main cause of the NTC effect and two-stage ignition. At these temperatures, the reactions shown with blue arrows in Fig. 4 become dominant. The net effect of these reactions is the build-up of H 2 O 2 and the general build-up of heat from many exothermic reactions. After temperatures reach about 800 K, the decomposition of H 2 O 2 to OH becomes favorable enough that the net flux of OH becomes strongly positive. This influx of OH radicals propels the system to final ignition.

Comparison of RCM and engine-like auto-ignition
While all three chemistry regimes discussed in Section 3.1 are controlling at various times in a constant-volume reactor, it is unclear which, if any, regimes are rate-limiting for the Fig. 3 Comparison of the simulated ignition delays (lines) of n-butane with the experimental data (points) by Healy et al. 37 Equivalence ratio is 1. Diluent to oxygen ratio is 3.76 corresponding to that of air. The diluent is nitrogen in the simulation. engine-like simulation. Comparing a simulation of the analogous RCM to the engine-like simulation gives a good basis for discussion. The RCM simulation uses the same volume profile derived for the engine-like simulation but stops compression at À55 CAD (t = 32 ms), at which time the temperature is approximately T = 650 K. The engine-like simulation completes the full engine stroke (with additional compression from the burning gas after À15 CAD) through 50 CAD. Fig. 5(a) shows temperature profile of an RCM and engine-like simulation for pure butane. The RCM simulation takes approximately 10 times longer to ignite than the engine-like, and spends significantly more time below 800 K. In an RCM, the ignition is successively limited by the three cycles shown in Fig. 4: chain-branching through ketohydroperoxide intermediates, heat generation by secondary chemistry, and finally OH dissociation from H 2 O 2 . On the other hand, the engine-like simulation seems to push past the first two regimes, as even pure nitrogen reaches a max temperature of B1100 K in these conditions. This contrast can be further seen in Fig. 5(b), which presents the OH mole fraction as a function of temperature. For the RCM simulation, the OH mole fraction peaks around 750 K (after the ketohydroperoxide cycle has become unfavorable), and does not start increasing again until after 800 K. Meanwhile in the enginelike simulation, the mole fraction of OH increases monotonically, but is much lower for 650 K o T o 800 K compared to RCM. Above 800 K, the OH mole fraction experiences similar increases in the two simulations. This comparison shows that the enginelike simulation rushes through the first two regimes, but experiences similar radical growth during the regime controlled by H 2 O 2 dissociation. As a result, the engine-like simulation will build up less H 2 O 2 prior to the third regime. Fig. 6 shows the simulated ignition delays at constant volume adiabatic conditions of stoichiometric f = 1 mixtures in air of the different additive/n-butane blends with 2 mol% additive in fuel. The initial conditions are P 0 = 20 bar and T 0 = 650-1000 K, and we take the ignition time to be when dP/dt achieves its maximum. Ethylphenol and m-cresol are predicted to have very little effect on the ignition throughout all temperature ranges. Guaiacol has ignition enhancing effects throughout all temperature ranges. The other three additives, p-cresol, o-cresol, and xylenol have a more complex effect. They accelerate ignition for T 0 o 780 K. However, above 780 K these three additives slow ignition. Based solely on these simulations, it is unclear how to rank the anti-knock behavior of the additives. A real engine will cross the entire temperature range during a piston cycle, but analysis from Section 3.1 suggests that the temperature range above 800 K is most relevant to the CFR engine.

Prediction of the anti-knock tendency
To address the shortcomings of the constant-volume simulation, engine-like simulations following the methodology outlined in Section 2.2 were performed to investigate the ignition timing of the mixture when being exposed to the end gas condition. All simulations were performed at f = 1, T 0 = 400 K, P 0 = 1.09 bar (taken from experimental CFR data) and a blending ratio of 2 mol% additive in n-butane. Fig. 7 shows the temperature history of the simulations. Guaiacol ignites fastest of all of the additive blends, which is consistent with the constant volume simulation results. Both m-cresol and ethylphenol have small effects, but do provide some anti-knock behavior. o-Cresol and p-cresol have about the same anti-knock behavior which agrees well with the constant-volume simulation. Among all of the additive blends, 2,4-xylenol has the longest ignition delay, representing the best anti-knock capability. Based on Fig. 7, the anti-knock tendency of the additives can be arranged as: 2,4-xylenol 4 p-cresol = o-cresol 4 m-cresol = 2-ethylphenol 4 guaiacol. This ranking suggests that constant-volume simulations starting above 800 K are more relevant to the RON when RON 4 90.
Following the RON test protocol, engine experiments were performed in a CFR engine to test the anti-knock tendency of the additives. Experimental RON's measured from the CFR test are summarized in Table 3. The blends were tested over several days explaining the slight difference in measured base fuel RON. However, the base fuel was tested just before each blend to give accurate DRON. Based on these experiments, the additives from most to least anti-knock performance are: 2,4-xylenol 4 p-cresol 4 o-cresol 4 m-cresol 4 2-ethylphenol 4 guaiacol. We note that experimental determinations of DRON between fuel blends have uncertainty of about AE0.5.
Quantitative predictions for the RON of additives blend can be made by interpolating the ignition timings of the additives into those of PRFs in Fig. 2. There is some error in this prediction because the merged model does not exactly match  Table 3. Qualitatively, these predictions compare favorably with the experimental data, predicting DRON for all six additives within 0.8 RON units and ranking the anti-knock behavior of the additives correctly except for switching m-cresol and ethylphenol. Free radicals, mainly OH radical and HO 2 radical, attack and abstract the hydrogen atom in the substituted methyl group and hydroxyl group, producing 2-methylphenoxy radical and 2-hydroxybenzyl radical. These two radicals then go through disproportionation reactions, consuming another radical and forming a conjugated ketone. These two steps are denoted as pathway P0 in Fig. 8. The reaction sequence in P0 net destroys two radicals. This conjugated ketone is relatively stable, so the flux toward it is much higher than the secondary chemistry afterwards.

Important reaction pathways of the additives
There are three main consumption pathways of the conjugated ketone, denoted as P1, P2, and P3 respectively in Fig. 8. In the pathway P1, the ketone adds a methyl radical to its terminal carbon and disproportionates. Next, HO 2 abstracts to form H 2 O 2 , which at this temperature quickly dissociates into OH radicals. The P2 pathway adds methyl radical to the benzene ring instead. The resulting alkoxy radical performs intramolecular hydrogen migration moving the radical to carbon at the ortho position. Finally, oxygen disproportionates the alkyl radical to form another conjugated ketone and the lower reactivity radical HO 2 . The P3 pathway adds H radical to the benzene ring and is followed by addition of oxygen to form peroxy radicals. Including the two radicals terminated by P0, the P1 pathway is radical neutral, P2 consumes 2 radicals, and P3 consumes 3 radicals.
Sensitivity analysis for OH radical was performed at the same conditions and time point. Because the phenolic compounds are in low concentration, the most sensitive reactions in the blend coincide with those of pure butane. To emphasize the contribution of the additives, the 50 most sensitive reactions from butane simulations were removed from consideration. The remaining sensitive reactions are shown in Fig. 9. The three most sensitive Fig. 6 Simulated ignition delays of the additive/n-butane blends in a constant volume adiabatic homogenous batch reactor as a function of initial temperature T 0 . The blend ratio is 2%, f = 1, and initial pressure is 20 bar. Fig. 7 The predicted temperature history of the engine-like simulations using different fuel blends. PRF90: 90% iso-octane/n-heptane; PRF100: pure iso-octane; pure n-butane; 2% additive/n-butane. Conditions: f = 1, T 0 = 400 K, P 0 = 1.09 bar, adiabatic, effective V(t) defined by RON100 P(t).
reactions occur on the P0 pathway. The remaining reactions only show secondary reactions of butane. It follows that the secondary chemistry of o-cresol has a much smaller effect on the ignition than the initial P0 pathway.
Local first order uncertainty analysis was performed at the same condition and time point with respect to OH radical using the uncertainty analysis module newly developed in RMG. 40 The uncertainty of the OH radical concentration predicted by the RMG generated model originates from both the uncertainty of thermochemistry and the uncertainty of the forward rate coefficients; here all the uncertainties were treated as uncorrelated.
As shown in Fig. 10a, the most important uncertain thermochemistry is that of the methylperoxy radical. However, this has been fairly well-studied e.g. in the active thermochemical tables, 43 so it would not be easy to significantly reduce its error bar. The uncertainty contribution of kinetics, as shown in Fig. 10b, is dominated by reaction R1, disproportionation of 2-methylphenoxy radical by hydroperoxyl radical producing o-quinone methide and hydrogen peroxide, and reaction R2, which is the decomposition of hydrogen peroxide. The kinetics of the latter reaction is well studied and considered to have a small uncertainty. However, the kinetics of the former reaction Table 3 Comparison of predictions by engine-like simulation and experimentally determined DRON for fuel blends with additives. The top set of data corresponds to prediction of standard RON engine-like simulations with 2% mole fraction phenol additives in butane. RON was computed by linear interpolation of calculations for different PRF fuels in our engine simulations, using the LLNL PRF mechanism, see Fig. 2  has not been studied in the literature; here it is estimated by averaging two rate rules in the disproportionation reaction family. Therefore we suggest work to improve the accuracy of the rate coefficient of reaction R1 to reduce the uncertainty of the current mechanism. Global uncertainty analysis for the o-cresol case was also performed on selected highly uncertain parameters, using the methodology and uncertainty estimation procedure detailed by Gao. 40 The uncertainty of the kinetics of reaction R1 and the thermochemistry of o-quinone methide were estimated based on their sources, and subsequently propagated to compute the uncertainty of OH radical mole fraction, which was calculated to be 26%. In this case, the uncertainty estimates from global uncertainty analysis are similar to the results from local first-order analysis.
The main reaction pathways for p-cresol, 2,4-xylenol, and ethylphenol are analogous to o-cresol. As shown in Fig. 11, their consumption also starts with the hydrogen abstraction reaction and is followed by disproportionation, resulting in conjugated ketones; their own version of the P0 pathway. The secondary chemistry is not shown as, similarly to o-cresol, it has a much smaller effect than the P0 pathways featured. Because 2,4-xylenol has two methyl groups, either of which can form conjugated ketones, its P0 sequence has a higher reaction rate and stronger anti-knock effect than either p-cresol or o-cresol. Sensitivity analyses for these other additives, at engine-like conditions, are given in ESI; † these also show that the P0 pathways have the largest influence on ignition.
The main reaction pathway for guaiacol at 925 K is direct dissociation to methyl radical and alkoxy radical as shown in Fig. 12. This one reaction is the largest contributor to its cetane boosting effect at the conditions shown. Some of the o-hydroxy phenoxy radicals formed abstract H atoms, further accelerating ignition, but some recombine or disproportionate with O 2 or  other radicals reducing the reactivity. Two secondary pathways of guaiacol are also significant. H-Abstraction from the phenol group induces O-C cleavage, so it is radical neutral. H-Abstraction from the methyl group accelerates ignition by adding O 2 at low T, or by releasing formaldehyde and a reactive hydroxy phenyl radical as shown on the right side of Fig. 12. Sensitivity analysis shown in Fig. 13 shows the bond dissociation pathway having the largest influence on OH concentration. At this time point a different HO 2 disproportionation to form conjugated ketone is important. This gives a large negative sensitivity coefficient as HO 2 forms O 2 instead of H 2 O 2 in this sequence. The third most sensitive reaction shows methyl peroxy radical disproportionating with o-hydroxy phenoxy to form conjugated ketone. This gives positive sensitivity for OH sensitivity because the resulting methyl hydroperoxide will dissociate rapidly at these temperatures to methoxy and OH.
The reaction pathway of m-cresol is different from those of the other five additives as it cannot follow the P0 pathway. As shown in Fig. 14, after the hydrogen abstraction from the methyl site, 3-hydroxybenzyl radical cannot produce the conjugated ketone through disproportionation reaction because of the positioning of the methyl group. Instead, the 3-hydroxybenzyl radical will combine with HO 2 . The resulting hydroperoxide breaks down to OH radical at these temperatures. The competing pathway forming 3-methylphenoxy radical can still form conjugated ketones but only by adding a radical, and there are only three conjugated bonds as opposed to the four seen for all the other additives. The reduced conjugation of its intermediates and the competing pathway forming OH, gives meta-cresol a lower net octane-boosting effect than the previous additives.
OH radical is the most important intermediate species in the ignition process. The anti-knock ability of additives mainly originates from the fact that they can quench OH radicals or precursors effectively. The reaction path analysis shown in Fig. 8 reveals the importance of conversion of the reactive phenol to a stable conjugated ketone in the anti-knock behavior of substituted phenols. Two active radicals, one during initial H-abstraction and a second during subsequent disproportionation, are consumed in the formation of the conjugated ketone through the first two steps (P0). Most of the produced conjugated ketone will stay relatively inert and not be consumed until hot ignition. Therefore, the formation of the conjugated ketone is the key to understand the anti-knock ability of these additives. p-Cresol and o-cresol have similar P0 pathways, thus they have similar anti-knock tendency. 2,4-Xylenol has three channels for the hydrogen abstraction reaction and two conjugated ketones as the products of the pathway P0, which means it consumes active radicals faster than p-cresol and o-cresol. Therefore, it has a higher anti-knock tendency. Guaiacol has a decomposition channel which creates methyl radicals on the way to its conjugated ketone, net producing radicals rather than consuming them. m-Cresol cannot do the P0 pathway at all, giving it the smallest effect of any of the additives.
The above understanding, which is that the formation of conjugated ketone radical and radical consumption via the P0 pathway leads to the anti-knock ability of substituted phenols, can help in the search for better additives. For example, we speculate that 2,4,6-trimethylphenol will have better anti-knock ability than 2,4-xylenol according to the above analysis.
The main pathways leading to the anti-knock ability were explored by RMG without any prior knowledge for specific additives. Although there is no direct study focusing on the oxidation kinetics of cresols, a few earlier studies 41,42 investigated the oxidation of xylenes, which have structures similar to those of cresols. Prior researchers proposed that hydrogen abstraction followed by disproportionation is one of the main reaction pathways in the oxidation of xylenes, 41,42 which agrees Fig. 11 The P0 reaction pathways of p-cresol, 2,4-xylenol, and xylenol.

Conclusion
A detailed kinetic model was automatically built using the reaction mechanism generator (RMG) for the study of antiknock tendency of n-butane blended with six additives: p-cresol, o-cresol, m-cresol, 2,4-xylenol, 2-ethylphenol, and guaiacol. Ignition delays were simulated in a constant volume batch reactor.
It was found that some of the additives have opposite performance in different temperature ranges. Therefore, engine-like simulations were performed to investigate the ignition timing of the blends at end gas condition. These simulations were conducted by specifying a volume history, which is derived from a pressure history obtained in a RON test of PRF100, to a batch reactor. Based on the predicted ignition timing in the engine-like simulation, the anti-knock tendency of the additives can be arranged as: 2,4-xylenol 4 p-cresol = o-cresol 4 m-cresol 4 2-ethylphenol 4 guaiacol. Quantitative predictions of RON change induced by the additives were made by interpolating the ignition timings of the additives into those of PRFs. The procedure for predicting the increase in RON agrees with experimental measurements on blends of these additives in real gasoline in a standard RON test almost within experimental uncertainty.
The key mechanism determining the anti-knock behavior was found to be the consumption of radicals as the phenol is converted into a conjugated ketone. This consists of a two-step pathway: hydrogen abstraction reactions of additives followed by disproportionation, consuming two net radicals. The differences in the formation of conjugated ketone and subsequent secondary chemistry were used to explain the anti-knock ranking of p-cresol, o-cresol, 2,4-xylenol, 2-ethylphenol, and guaiacol.
This study is a first attempt to apply the automated mechanism generation technique to study of anti-knock additives at realistic engine-like conditions. The good consistency of the  predictions and the experimental results shows this method is a promising way to rank proposed anti-knock additives even before performing any experiments. It is hoped that this example will encourage future efforts to bring advanced chemical kinetics techniques to bear on real-world problems for the benefit of society.

Conflicts of interest
There are no conflicts to declare.