Reaction Rate Rules of Intramolecular H-Migration Reaction Class for RIORIIOO·Radicals in Ether Combustion

Sun, Xiaohui; Li, Zerong

doi:10.3390/molecules29184387

Open AccessArticle

Reaction Rate Rules of Intramolecular H-Migration Reaction Class for R^IOR^IIOO·Radicals in Ether Combustion

by

Xiaohui Sun

^1,* and

Zerong Li

^2,*

¹

School of Energy Engineering, Shanxi College of Technology, Shuozhou 036000, China

²

College of Chemistry, Sichuan University, Chengdu 610064, China

^*

Authors to whom correspondence should be addressed.

Molecules 2024, 29(18), 4387; https://doi.org/10.3390/molecules29184387

Submission received: 19 August 2024 / Revised: 8 September 2024 / Accepted: 13 September 2024 / Published: 15 September 2024

(This article belongs to the Section Physical Chemistry)

Download

Browse Figures

Versions Notes

Abstract

:

The intramolecular H-migration reaction of R^IOR^IIOO· radicals constitute a key class of reactions in the low-temperature combustion mechanism of ethers. Despite this, there is a dearth of direct computations regarding the potential energy surface and rate constants specific to ethers, especially when considering large molecular systems and intricate branched-chain structures. Furthermore, combustion kinetic models for large molecular ethers generally utilize rate constants derived from those of structurally similar alcohols or alkane fuels. Consequently, chemical kinetic studies involve the calculation of energy barriers and rate rules for the intramolecular H-migration reaction class of R^IOR^IIOO· radicals, which are systematically conducted using the isodesmic reaction method (IRM). The geometries of the species participating in these reactions are optimized, and frequency calculations are executed using the M06–X method in tandem with the 6–31+G(d,p) basis set by the Gaussian 16 program. Moreover, the M06–2X/6–31+G(d,p) method acts as the low-level ab initio method, while the CBS–QB3 method is utilized as the high-level ab initio method for calculating single-point energies. Rate constants at the high-pressure-limit are computed based on the reaction class transition state theory (RC-TST) by ChemRate program, incorporating asymmetric Eckart tunneling corrections for intramolecular H-migration reactions across a temperature range of 500 to 2000 K. It was found that the isodesmic reaction method gives accurate energy barriers and rate constants, and the rate constants of the H-migration reaction for R^IOR^IIOO· radicals diverge from those of comparable reactions in alkanes and alcohol fuels. There are significant disparities in energy barriers and rate constants across the entire reaction classes of the H-migration reaction for R^IOR^IIOO· radicals, necessitating the subdivision of the H-migration reaction into subclasses. Rate rules are established by averaging the rate constants of representative reactions for each subclass, which is pivotal for the advancement of accurate low-temperature combustion reaction mechanisms for ethers.

Keywords:

intramolecular H-migration reaction; reaction rate rules; isodesmic reaction method (IRM); reaction class transition state theory (RC-TST)

1. Introduction

The burgeoning demand for energy, coupled with the limited supply of fossil fuels and their detrimental environmental consequences, has directed scientific inquiry towards the pursuit of alternative fuels as a means to gradually supplant conventional fossil-based energy sources. Biofuels, as an emerging category of alternative fuels, have increasingly captured the spotlight due to their renewable characteristics and their capacity to mitigate net CO₂ emissions [1,2,3,4]. A thorough comprehension of the combustion dynamics of biofuels is essential for the optimal utilization of these advanced renewable energy sources. Typically, investigating the combustion process of biofuels (and their blends) through experimental methods is costly and time-intensive. However, theoretical research provides a cost-effective and efficient means to simulate the combustion process of fuels across a broad spectrum of pressures, temperatures, and equivalence ratios, based on a detailed combustion reaction mechanism [2]. The accuracy of these simulations hinges on the completeness of the reaction pathways and the precision of the thermodynamic parameters within the combustion reaction mechanism. Nevertheless, a comprehensive reaction mechanism encompasses hundreds of species and thousands of elementary reactions, which are typically constructed using automated mechanism generation, tailored to various reaction classes [5]. Additionally, the reaction mechanism within mechanism programs is commonly segregated into two components: the C0–C4 core mechanism and the C5 and above extension mechanism (which includes reactions involving larger species), both generated in accordance with reaction rate rules [6,7,8,9]. When there exists a substantial disparity in energy barriers within a reaction class, it is further subdivided into distinct reaction subclasses, with specific rate rules assigned to each subclass [10].

Numerous studies have been conducted on the reaction mechanisms of hydrocarbon fuels to date [11,12,13,14]. In comparison, research into the combustion chemistry of biofuels is relatively sparse compared to that of hydrocarbon fuels. Ethers, a type of biofuel, have a general formula of R^IOR^II, wherein R^I and R^II represent alkyl groups that may be of identical or distinct structures. The low-temperature oxidation process of ether fuels, as shown in Scheme 1, is similar to the low-temperature reaction mechanism of alkanes. To differentiate the radicals generated from the oxidation of alkanes, the expression method for radicals produced during the oxidation of ethers in this study employs the same formula as that used by Tommaso et al. [3]. Here, alkyl radicals (R·), alkyl-peroxy radicals (ROO·), and hydroperoxyl-alkyl radicals (QOOH) are denoted as R^IOR^II·, R^IOR^IIOO·, and ·R^IOR^IIOOH, respectively. The detailed one-step oxygenation reaction process is illustrated in Scheme 1. Initially, ethers are consumed via abstraction reactions by hydrogen atoms or hydroxyl radicals, resulting in the generation of R^IOR^II· radicals. The activation energy for H-atom abstraction is low, and the low-temperature oxidation activity is higher than that of alkane fuels with the same chain length, owing to the presence of the oxygen atom in R^IOR^II· radicals, which weakens the α-C-H bond (near the oxygen atom). Analogous to alkyl radicals, R^IOR^II· radicals react with oxygen molecules to form R^IOR^IIOO· radicals, which can undergo isomerization reactions leading to the formation of ·R^IOR^IIOOH radicals. This class of reactions is not only crucial in the oxidation of ethers but also in the analogous process for hydrocarbons, and it is also one of the factors contributing to the negative temperature coefficient (NTC) behavior of ethers in low-temperature combustion.

The oxidation mechanism of ether fuels plays a key role under atmospheric conditions, and it has also been reported in some studies in the literature [15,16,17,18]. For instance, Burke et al. [17] conducted measurements on the ignition delay and sensitivity analysis of dimethyl ether at high pressures. They discovered that intramolecular H-migration reactions involving R^IOR^IIOO· radicals are significant in low-temperature combustion chemistry. Nevertheless, high-pressure-limit rate constants for this class of reactions are scarce, with only a few available for small molecular systems of R^IOR^IIOO· radicals. Dimethyl ether (DME) and diethyl ether (DEE) are the most extensively studied members of this class, both experimentally and theoretically. Adamo et al. [18] proposed an autoxidation reaction mechanism for diethyl ether (DEE) using density functional theory (DFT) and the B3LYP/6–31+G(d,p) method to calculate the rate constants for reactions. However, there is a lack of direct calculation of the potential energy surface and rate constants for ether fuels, particularly for large molecular systems and complex branched-chain ethers, which possess a higher calorific value and are more suitable as direct fuels [19]. Due to the limited availability of accurate ab initio methods for ethers, the rate constants in combustion kinetic models of ethers typically employ rate constants derived from alcohols or alkane fuels with similar structures. Several theoretical studies have investigated the role that the different groups play in the H-migration reaction of peroxy radicals, such as Vereecken et al. [20] and Otkjær et al. [21]. The results show that when the carbon atoms connected by the extracted hydrogen atoms are replaced by different substitutes, the rate constants are very different, so it is necessary to calculate the rate constant accurately for ether fuels.

It is widely acknowledged that high-level quantum chemistry methods, including the Gaussian-n method [22,23,24,25], High accuracy Extrapolated Ab initio Thermochemistry (HEAT) method [26], Complete Basis Set (CBS) method [27,28], and the Coupled Cluster theory with Single and Double excitations and a quasi-perturbative treatment of connected Triple excitations [CCSD(T)] method, are extensively employed for calculating electronic energies with chemical accuracy of roughly 1 kcal/mol [29]. Nonetheless, these high-level ab initio methods are best suited for small- to medium-sized molecules and are computationally intensive. Consequently, accurately determining single-point energies for large molecular systems remains a formidable challenge. Building upon our prior research [30], we have applied the isodesmic reaction method (IRM) in conjunction with the reaction class transition state theory (RC-TST) to compute rate constants for concerted elimination reaction classes of hydroperoxyl-alkyl-peroxyl radicals, achieving results closely matching those from high-level methods. In the present study, we utilize the isodesmic reaction method in combination with reaction class transition state theory to determine accurate energy barriers and rate constants at high-pressure-limits for intramolecular H–migration reaction of R^IOR^IIOO· radicals. Within reaction class transition state theory, a representative reaction of smaller size within the reaction class is chosen as the reference reaction, while other reactions within the class are considered as target reactions. Then, the high-pressure-limit rate constants for the target reactions are computed using a low-level ab initio method, and subsequently, these results are refined through the isodesmic reaction correction scheme.

The objectives of this study are as follows: (1) To supply high-precision reaction barriers and rate constants at high-pressure-limits for the intramolecular H–migration reaction of R^IOR^IIOO· radicals. (2) To furnish reaction rate rules at high-pressure-limits for the intramolecular H–migration reaction of R^IOR^IIOO· radicals. These findings will lay a theoretical groundwork for the development of more reliable combustion mechanisms for ether fuels at low temperatures. (3) To systematically investigate and compare the rate constants for the intramolecular H–migration reaction of ROO·, R^IOR^IIOO·, and HOROO· radicals, thereby offering theoretical insights to enhance the application of biofuels.

2. Results and Discussion

In this study, we have selected a total of 41 reactions involving R^IOR^IIOO· radicals of C2–C8, as detailed in Table 1. These reactions are categorized into four classes based on the distance between the carbon atom that the hydrogen atom is shifted from, which determines the ring size of transition states. To reduce the uncertainty of the rate rules for these reaction classes, they are further subdivided into three subclasses, based on the type of carbon atoms, which connected the migrating hydrogen atoms. Specifically, the carbon atom that is connected with three hydrogen atoms is defined as the “p” site, whereas a carbon atom that is connected with two hydrogen atoms and one hydrogen atom is defined as the “s” and “t” site, respectively. Consequently, the 41 intramolecular H-migration reactions are classified into the following subclasses: 1,n-H(p), 1,n-H(s), and 1,n-H(t), where n equals 3, 5, 6, or 7. In the text, tables, and figures, a radical site is denoted by a bullet symbol.

The reaction centers for the transition states of various reaction classes are depicted in Scheme 2. It is evident that the reaction center for these transition states forms a ring structure, and the 1,n-H migration reactions proceed through transition states with (n + 1) member-ring sizes.

2.1. Energy Barriers

2.1.1. Validation of the Energy Barriers by the CBS–QB3 Method

The precision of energy barriers is crucial for ensuring the reliability of rate constants. High-precision computational methods, such as the CCSD(T) method, are capable of calculating these energy barriers and reaction enthalpies with an accuracy of 1 to 2.0 kcal/mol [29]. However, this method is time-consuming when calculating electronic energies and is not practical for large molecular systems. Consequently, the CBS–QB3 method has been employed as the high-level ab initio method in this study. To confirm the accuracy of the energy barriers computed by the CBS–QB3 method, reactions R1, R2, R5, and R21 from Table 1 are randomly selected for comparison with the benchmark CCSD(T)/cc-PVTZ method. The data presented in Figure 1 illustrate that the energy barrier discrepancies between the CCSD(T)/cc-PVTZ and CBS–QB3 methods for reactions R1, R2, R5, and R21 are 1.43 kcal/mol, 1.47 kcal/mol, 1.36 kcal/mol, and 1.51 kcal/mol, respectively. These values fall within the range of chemical accuracy, which is typically defined as 1 to 2 kcal/mol [29,31]. These findings suggest that the energy barriers obtained using the CBS–QB3 method are acceptable. Furthermore, the energy barrier difference of 1.51 kcal/mol results in a reaction rate constant ratio of 3.12 at the standard ignition temperature of 500 K, which increases by a factor of 3.68 at the standard combustion temperature of 1500 K. Both values fall within the order of magnitude and are within the margin of error.

2.1.2. Energy Barriers for R^IOR^IIOO· Radicals

In this study, the first reaction within each subclass has been chosen as the reference reaction, with the remaining reactions in each subclass considered as target reactions. The disparities in energy barriers and reaction enthalpies between the low-level M06–2X method and the high-level CBS–QB3 method for the reference reactions of each subclass are detailed in Table S1 of the Supporting Information. The energy barriers for all reactions presented in Table 1 have been computed using the M06–2X method, incorporating isodesmic reaction corrections, and are listed in Table S2 of the Supporting Information. Utilizing the data from Table S2, the average value and maximum deviation of energy barriers for each subclass are compiled in Table 2. Upon Table 2, with the exception of the 1,3-H migration reaction class, the maximum deviations in energy barriers for 1,5-H, 1,6-H, and 1,7-H are 7.19 kcal/mol, 6.85 kcal/mol, and 6.70 kcal/mol, respectively. When the reaction class is further subdivided into reaction subclasses, the range of maximum deviations in energy barriers spans from 0.23 to 2.69 kcal/mol. Consequently, it is both logical and essential to categorize the reaction class into distinct reaction subclasses.

To verify the precision of the correction scheme for energy barriers, a selection of 11 representative target reactions from Table 1 is made by contrasting the energy barriers computed using the isodesmic reaction method against those calculated by the CBS–QB3 method. The results are detailed in Table 3. Inspection of Table 3 reveals that the discrepancies in energy barriers between the isodesmic reaction method and the CBS–QB3 method range from −0.70 to 0.56 kcal/mol, figures that fall within the bounds of chemical accuracy [29]. This suggests that the energy barriers determined by the isodesmic reaction method are indeed acceptable.

2.2. Reaction Enthalpies

In this study, the enthalpies associated with intramolecular H–migration reaction for R^IOR^IIOO· radicals are calculated employing the isodesmic reaction method. To ensure the reliability of the correction scheme embedded within the isodesmic reaction method, 10 reactions are randomly selected from Table 1 for verification. Both the CBS–QB3 method and the isodesmic reaction method are utilized to compute the reaction enthalpies for these selected reactions. The results are detailed in Table 4. As can be observed from Table 4, the deviations for reaction enthalpies calculated by the isodesmic reaction method in conjunction with the CBS–QB3 method range from −0.18 to 0.32 kcal/mol, a range that falls within the bounds of chemical accuracy (1~2 kcal/mol) [29]. Furthermore, for reaction R10, the reaction enthalpy computed by the isodesmic reaction method in this study amounts to 9.77 kcal/mol, which is in close proximity to the 9.30 kcal/mol reported by Curran et al. [32]. Therefore, it is concluded that the correction scheme of the isodesmic reaction method is viable for calculating precise energy barriers and reaction enthalpies for large molecular systems, achieved by adjusting values from low-level ab initio methods. The reaction enthalpies for all reactions are listed in Table S2 of the Supporting Information.

2.3. Rate Constants and Rate Rules at High-Pressure-Limit

In this study, the high-pressure-limit rate constants are performed with the ChemRate program [33]. Subsequently, these rate constants are fitted to the modified Arrhenius equation [11]:

k = A T^{n} \exp (- E / R T)

in the form of (A, n, E) over a temperature range of 500 to 2000 K, in increments of 100 K. Here, A represents the pre-exponential factor, E denotes the activation energy, and R is the gas constant.

The rate rules for each subclass are derived using expression (1) and are presented in the format (A, n, E), as listed in Table S3 of the Supporting Information.

k_{r u l e}^{T} = \frac{1}{N} \sum_{i = 1}^{i = N} k_{i}^{T}

(1)

Here, k_rule and k_i represent the average rate constant and the rate constant of any reaction i within a subclass at temperature T, respectively. The total number of reactions within the reaction subclass is denoted by N.

2.3.1. Comparison of Rate Constants by Isodesmic Reaction Method with CBS–QB3 Method

In order to evaluate the reliability of the rate constants derived from the modified scheme using the isodesmic reaction method, the rate constants computed by both the isodesmic reaction method and the CBS–QB3 method for reactions R24, R28, R32, R36, and R40 from Table 1 are compared. A ratio factor μ (μ = k_IRM/k_CBS) is employed to assess the deviation of rate constants across the temperature range of 500 to 2000 K. To provide a clearer insight, Figure 2 presents a comparative analysis of rate constants derived using the isodesmic reaction method and the CBS–QB3 method at temperatures of 500 K, 1000 K, and 1500 K. The data reveal that the ratios of the rate constants are consistently within the order of magnitude, suggesting that the rate constants derived from the isodesmic reaction method closely align with those yielded by the CBS–QB3 method. Consequently, the correction scheme at the M06–2X/6–31+G(d,p) level of theory can be reliably employed for calculating accurate rate constants for the reactions of interest.

2.3.2. Comparison of the Rate Constants with Values in the Literature for Alkanes

In the preceding introduction, we noted that the current reaction mechanism for ether fuels frequently utilizes the rate constants of analogous alkane reactions that are devoid of specific rate constants. This approach can lead to considerable inaccuracies in the reaction mechanism of ether fuels. Consequently, we have chosen four reactions from Table 1 to assess the discrepancies in the rate constants for H–migration reactions of R^IOR^IIOO· radicals in ethers and comparable reactions involving ROO· radicals in alkanes, as detailed by Villano et al. [7] and Bugler et al. [34]. The results are presented in Table 5. Inspection of Table 5 reveals that the rate constants at the high-pressure-limit, as calculated in our study, deviate from those reported by Villano et al. [7] and Bugler et al. [34], with the ratio varying from 1.06 to 2.11 × 10⁴ at 500 K. Therefore, the impact of the R^IOR^II group cannot be ignored. It is imperative to establish rate rules at the high-pressure-limit for R^IOR^IIOO· radicals, which are also vital for the advancement of reaction mechanisms in automated reaction mechanism generators.

2.3.3. Comparison of the Rate Constants in Analogous Reaction Subclass for R^IOR^IIOO· in Ethers with HOROO· in Alcohols and ROO· in Alkanes

As outlined in the introduction, the prevalent reaction mechanism for ether fuels typically employs kinetic parameters derived from analogous reactions involving alcohols or alkanes. Kerschgens et al. [35] observed that the ignition properties of ether fuels diverge from those of alcohol fuels and alkanes. Furthermore, the comparison of different substitutes on the carbon atom that the hydrogen atom abstracted is conducted by Otkjær et al. [21]. Consequently, a comparison has been made between the kinetic parameters for 1,6-H migration of R^IOR^IIOO· radicals in ethers, the HOROO· radicals in alcohols [36], and the ROO· radicals in alkanes [34], as depicted in Figure 3. Figure 3 illustrates that kinetic parameters for the same reaction subclass vary when catalyzed by different peroxy radicals. Notably, the kinetic parameters for HOROO· radicals in alcohols are the highest among the three types of peroxy radicals. In the case of the 1,6-H(p) reaction subclass, the rate constant of the hydrogen migration reaction of R^IOR^IIOO· radicals are the smallest. For the 1,6-H(s) and 1,6-H(t) subclasses, the reactivity of the HOROO· and R^IOR^IIOO· radicals surpass that of the ROO· radical, suggesting that the presence of oxygen-bearing groups reduces the bond energy of the C-H bond. This conclusion is consistent with the findings of Otkjær et al. [21]. Furthermore, we also contrasted the rate constants for the hydrogen migration reaction within the same peroxy group, considering the varying types of carbon atoms, as depicted in Figure 4. From Figure 4, the rate constant for the 1,6-H migration reaction of peroxy radicals with varying types of carbon atoms exhibits the following trend: tertiary > secondary > primary. This outcome aligns with the observations made by Otkjær et al. [21]. The reactions forming secondary radicals are about seven orders of magnitude faster than those forming primary radicals, and reactions leading to tertiary radicals are around 27 times faster than those leading to secondary radicals at 1200 K.

In summary, when considering the discussion on the variance of energy barriers within reaction classes as detailed in Section 2.1.2, further stratification of these classes into distinct reaction subclasses can diminish the uncertainty of the rate rules derived from the rate rule method. Concurrently, in the context of constructing the reaction mechanisms for large molecular ether fuels, it is inappropriate to employ rate constants from analogous reactions involving alkanes and alcohol fuels. Consequently, this paper proposes to formulate specific reaction rate rules for each subclass using the rate rule method, as elaborated upon in Section 2.3.4 below.

2.3.4. The Rate Rules at High-Pressure-Limit for R^IOR^IIOO· Radicals

Rate rules are crucial for constructing reaction mechanisms in automated mechanism generators. Typically, identical rate constants are assigned to reactions within the same class. However, this approach can occasionally yield significant errors, necessitating a further subdivision of reaction classes into distinct subclasses based on the varying energy barriers within each class. The high-pressure-limit rate rules for these different subclasses are detailed in Table 6. Furthermore, a ratio factor f, which compares the maximum rate constants to the minimum rate constants, is employed to assess the uncertainty of the rate rules within a subclass. According to Table 6, the uncertainty for the 1,n-H(p) subclass is substantial, with f ranging from 2.46 × 10³ to 4.71 × 10³ at 500 K and from 1.59 × 10³ to 2.63 × 10³ at 1000 K. Additionally, the uncertainty for the rate rules of the 1,n-H(s) and 1,n-H(t) subclasses falls between 1.27 to 25.2 and 1.98 to 14.20, respectively. This indicates that the rate constant derived using the reaction rate rule method exhibits a certain degree of deviation. It is essential to calculate the rate constant for each subclass precisely while considering the effect of the conformation for reactants and the tunneling effect.

3. Methods

3.1. Electronic Structure Calculation

All quantum chemical computations are performed using the Gaussian 16 program [37]. The geometries of the species involved in the reactions are optimized using the M06–2X method in conjunction with the 6–31+G(d,p) basis set. Vibrational frequencies are computed using the same method, with a scaling factor of 0.979 applied. Intrinsic Reaction Coordinate (IRC) analysis [38] is employed to ensure that the transition states are correctly connecting the minima associated with the products and reactants. The CBS–QB3 method is chosen as the high-level ab initio method for calculating the single-point energies of reference reactions in the isodesmic reaction method. This method demonstrates an accuracy exceeding 1 kcal/mol in predicting the heats of formation for a large set of molecules during the tests [39,40].

3.2. Calculation of the Rate Constant

Regarding the calculation of rate constants for large molecular systems, the reaction class transition state theory, developed by Truong [41,42] and their colleagues, is designed for elementary reactions within a class using a relatively low-level ab initio method. The fundamental principle of RC-TST is that both the reference reaction and any target reaction within the same class share the same reaction center. The rate coefficient k_r for the reference reaction is known, having been determined experimentally or through high-level ab initio methods, whereas the rate constants k_t for the target reactions are unknown and require calculation. In our previous work [30], we provided a detailed explanation of the reaction class transition state theory and defined it as the reaction class isodesmic reaction method. Here, we offer a brief description as follows:

The reference reaction and target reaction within the reaction class are established as

A_{r} + B_{r} = {TS}_{r} \cdot Δ V_{r}^{\neq}

(2)

A_{t} {+ B}_{t} {= TS}_{t} \cdot Δ V_{t}^{\neq}

(3)

where

Δ V_{r}^{\neq}

and

Δ V_{t}^{\neq}

represent the energy barriers for the reference reaction (2) and the target reaction (3), respectively. When subtracting reaction (2) from reaction (3), the resulting reaction (4) can also be regarded as an isodesmic reaction, where ΔΔV^≠ is the energy barrier for reaction (4).

A_{t} {+ B}_{t} {+ TS}_{r} {= A}_{r} + B_{r} + {TS}_{t} \cdot Δ Δ V^{\neq} = {Δ V}_{t}^{\neq} - {Δ V}_{r}^{\neq}

(4)

The isodesmic reaction method in the reaction class allows for the calculation of accurate energy barriers for target reactions R_t using a straightforward correction scheme based on reference reactions R_r. Within this scheme, reaction (4) is treated as an isodesmic reaction, from which the following expression can be derived:

Δ Δ V^{\neq} = {Δ V}_{t}^{\neq ’} - {Δ V}_{r}^{\neq ’} = {Δ V}_{t}^{\neq} - {Δ V}_{r}^{\neq}

(5)

where

Δ V_{t}^{\neq ’}

and

Δ V_{r}^{\neq ’}

represent the accurate energy barriers obtained from high-level ab initio methods or experimental data, and

Δ V_{t}^{\neq}

and

Δ V_{r}^{\neq}

denote the energy barriers calculated using low-level ab initio methods. Consequently, the following expressions can be derived:

Δ V^{\neq} = {Δ V}_{r}^{\neq ’} - {Δ V}_{r}^{\neq}

(6)

{Δ V}_{t}^{\neq ’} = {Δ V}_{t}^{\neq} + Δ Δ V^{\neq}

(7)

where ΔΔV^≠ represents the energy barrier value for reference reactions between the high-level and low-level ab initio methods. This value is utilized to adjust the energy barriers computed using the low-level method within the same reaction class for any target reactions, as per Equation (7). Following the correction with Equation (7), high-precision energy barriers can be derived from a low-level ab initio method, such as the DFT method. Likewise, accurate reaction enthalpies for the target reactions can be acquired using Equation (8).

{Δ H}_{t}^{\neq ’} = {Δ H}_{t}^{\neq} + Δ Δ H^{\neq}

(8)

Here,

Δ H_{t}^{’}

represents the high-precision reaction enthalpies for the target reaction, and

Δ H_{t}

are low-level reaction enthalpies for the target reactions. Furthermore, by adjusting the approximate rate constant k at a low-level ab initio method in accordance with expression (9), the exact rate constant k’ for any target reaction within the same reaction subclass can be derived.

k^{’} = k e x p (\frac{- Δ Δ V^{\neq}}{R T})

(9)

Here, ΔΔV^≠ represents the correction scheme for energy barriers, as derived from the reference reaction in expression (6).

4. Conclusions

In this study, we employed the isodesmic reaction method, grounded in reaction class transition state theory, to compute precise energy barriers and rate constants for H–migration reaction involving R^IOR^IIOO· radicals. We then constructed high-pressure-limit rate rules by averaging the rate constants of representative reactions within a subclass. Our findings lead to the following conclusions:

(1): After adjusting the energy barriers and reaction enthalpies for the target reaction at the M06–2X level using the isodesmic reaction method, the calculated values align closely with CBS–QB3 results and data in the literature.
(2): The overall reaction class exhibits a maximum deviation in energy barriers and rate constants that exceeds chemical accuracy, necessitating the subdivision of H–migration reaction into distinct subclasses.
(3): The high-pressure-limit rate constants derived in this work diverge from values in the literature for alkanes, with the ratio ranging from 1.06 to 2.11 × 10⁴ at 500 K. The presence of C-O-C bonds results in a weakening of the C-H bond adjacent to the ether group.
(4): By comparing the rate constants for ROO·, R^IOR^IIOO·, and HOROO· radicals, we determined that the rate constants for the H–migration reaction of HOROO· radicals are greater than those for ROO· and R^IOR^IIOO· radicals. The maximum deviation in rate constants between R^IOR^IIOO· and HOROO· radicals span two orders of magnitude. Consequently, it is imperative to establish separate rate rules for the H–migration reaction of R^IOR^IIOO· radicals, rather than adopting kinetic parameters from alcohols or alkane fuels with analogous reactions. Additionally, the rate constant derived using the reaction rate rule method exhibits a certain degree of deviation. It is essential to calculate the rate constant for the elementary reaction of ether fuels precisely.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/molecules29184387/s1: The differences in energy barriers and reaction enthalpies between the M06–2X method and CBS–QB3 method for reference reactions are listed in Table S1. The energy barriers and reaction enthalpies for all reactions are listed in Table S2. The high-pressure-limit kinetic parameters (A, n, E) for all reactions are listed in Table S3. The full names of abbreviations related to this article are listed in Table S4.

Author Contributions

Investigation, X.S.; Data curation, X.S.; Writing—original draft, X.S.; Writing—review & editing, Z.L.; Supervision, Z.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the Research Initiation Fund Programs of Shanxi College of Technology (No. 200101).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article and Supplementary Materials.

Conflicts of Interest

The authors declare no competing financial interests.

References

Pitsch, H.; Goeb, D.; Cai, L.; Willems, W. Potential of oxymethylene ethers as renewable diesel substitute. Prog. Energy Combust. Sci. 2024, 104, 101173. [Google Scholar] [CrossRef]
Michelbach, C.A.; Tomlin, A.S. Automatic mechanism generation for the combustion of advanced biofuels: A case study for diethyl ether. Int. J. Chem. Kinet. 2024, 56, 233–262. [Google Scholar] [CrossRef]
Tommaso, S.D.; Rotureau, P.; Adamo, C. Oxidation mechanism of aliphatic ethers: Theoretical insights on the main reaction channels. J. Phys. Chem A 2012, 116, 9010–9019. [Google Scholar] [CrossRef] [PubMed]
Li, T.; Li, J.Q.; Chen, S.Y.; Zhu, Q.; Li, Z.R. Investigating the kinetics of the intramolecular H-migration reaction class of methyl-ester peroxy radicals in low-temperature oxidation mechanisms of biodiesel. Phys. Chem. Chem. Phys. 2023, 25, 32078–32092. [Google Scholar] [CrossRef] [PubMed]
Van de Vijver, R.; Vandewiele, N.M.; Bhoorasingh, P.L.; Slakman, B.L.; Khanshan, F.S.; Carstensen, H.H.; Reyniers, M.F.; Marin, G.B.; West, R.H.; Geem, K.M.V. Automatic mechanism and kinetic model generation for gas and solution-phase processes: A perspective on best practices, recent advances, and future challenges. Int. J. Chem. Kinet. 2015, 47, 199–231. [Google Scholar] [CrossRef]
Miyoshi, A. Systematic computational study on the unimolecular reactions of alkylperoxy (RO₂), hydroperoxyalkyl (QOOH), and hydroperoxyalkylperoxy (O₂QOOH) radicals. J. Phys. Chem. A 2011, 115, 3301–3325. [Google Scholar] [CrossRef]
Villano, S.M.; Huynh, L.K.; Carstensen, H.H.; Dean, A.M. High-pressure rate rules for alkyl + O₂ reactions. 1. The dissociation, concerted elimination, and isomerization channels of the alkyl peroxy radical. J. Phys. Chem. A 2011, 115, 13425–13442. [Google Scholar] [CrossRef]
Villano, S.M.; Huynh, L.K.; Carstensen, H.H.; Dean, A.M. High-pressure rate rules for alkyl + O₂ reactions. 2. the isomerization, cyclic ether formation, and β-scission reactions of hydroperoxy alkyl radicals. J. Phys. Chem. A 2012, 116, 5068–5089. [Google Scholar] [CrossRef]
Villano, S.M.; Carstensen, H.H.; Dean, A.M. Rate rules, branching ratios, and pressure dependence of the HO₂ + olefin addition channels. J. Phys. Chem. A 2013, 117, 6458–6473. [Google Scholar] [CrossRef]
Yao, X.X.; Pang, W.Q.; Li, T.; Shentu, J.T.; Zhu, Q.; Li, X.X. High-Pressure-Limit and Pressure-Dependent Rate Rules for Unimolecular Reactions Related to Hydroperoxy Alkyl Radicals in Normal-Alkyl Cyclohexane Combustion. 2. Cyclization Reaction Class. J. Phys. Chem. A 2021, 125, 8959–8977. [Google Scholar] [CrossRef]
Yao, X.X.; Sun, X.L.; Zhu, Y.F. High-pressure limit and pressure-dependent rate rules for β-scission reaction class of hydroperoxyl alkyl hydroperoxyl radicals(•P(OOH)₂) in normal-alkyl cyclohexanes combustion. Molecules 2024, 29, 544. [Google Scholar] [CrossRef] [PubMed]
Battin-Leclerc, F. Detailed chemical kinetic models for the low-temperature combustion of hydrocarbons with application to gasoline and diesel fuel surrogates. Prog. Energy Combust. Sci. 2008, 34, 440–498. [Google Scholar] [CrossRef]
Ning, H.B.; Gong, C.M.; Tan, N.X.; Li, Z.R.; Li, X.Y. Low- and intermediate-temperature oxidation of ethylcyclohexane: A theoretical study. Combust. Flame 2015, 162, 4167–4182. [Google Scholar] [CrossRef]
Simmie, J.M. Detailed chemical kinetic models for the combustion of hydrocarbon fuels. Prog. Energy Combust. Sci. 2003, 29, 599–634. [Google Scholar] [CrossRef]
Orlando, J.J.; Tyndall, G.S.; Wallington, T.J. The atmospheric chemistry of alkoxy radicals. Chem. Rev. 2003, 103, 4657–4689. [Google Scholar] [CrossRef]
Atkinson, R. Rate constants for the atmospheric reactions of alkoxy radicals: An updated estimation method. Atmos Environ. 2007, 41, 8468–8485. [Google Scholar] [CrossRef]
Burke, U.; Somers, K.P.; O’Toole, P.; Zinner, C.M.; Marquet, N.; Bourque, G.; Petersen, E.L.; Metcalfe, W.K.; Serinyel, Z.; Curran, H.J. An ignition delay and kinetic modeling study of methane, dimethyl ether, and their mixtures at high pressures. Combust. Flame 2015, 162, 315–330. [Google Scholar] [CrossRef]
Di Tommaso, S.; Rotureau, P.; Crescenzi, O.; Adamo, C. Oxidation mechanism of diethyl ether: A complex process for a simple molecule. Phys. Chem. Chem. Phys. 2011, 13, 14636–14645. [Google Scholar] [CrossRef]
Yang, B.; Zhang, Y.J.; Li, Y.Y.; Chen, Z.; Huang, Z.H.; Qi, F. Research progress and prospect of combustion reaction kinetics for carbon neutrality. J. Eng. Therm. 2022, 43, 1993–2008. [Google Scholar]
Vereecken, L.; Nozière, B. H migration in peroxy radicals under atmospheric conditions. Atmos. Chem. Phys. 2020, 20, 7429–7458. [Google Scholar] [CrossRef]
Otkjær, R.V.; Jakobsen, H.H.; Tram, C.M.; Kjaergaard, H.G. Calculated hydrogen shift rate constants in substituted alkyl peroxy radicals. J. Phys. Chem. A 2018, 122, 8665–8673. [Google Scholar] [CrossRef] [PubMed]
Curtiss, L.A.; Raghavachari, K.; Trucks, G.W.; Pople, J.A. Gaussian-2 theory for molecular energies of first- and second-row compounds. J. Chem. Phys. 1991, 94, 7221–7230. [Google Scholar] [CrossRef]
Curtiss, L.A.; Redfern, P.C.; Raghavachari, K. Gn theory. Comput. Mol. Sci. 2011, 1, 810–825. [Google Scholar] [CrossRef]
Curtiss, L.A.; Raghavachari, K.; Redfern, P.C.; Rassolov, V.; Pople, J.A. Gaussian-3 (G3) theory for molecules containing first and second-row atoms. J. Chem. Phys. 1998, 109, 7764–7776. [Google Scholar] [CrossRef]
Curtiss, L.A.; Jones, C.; Trucks, G.W.; Raghavachari, K.; Pople, J.A. Gaussian-1 theory of molecular energies for second-row compounds. J. Chem. Phys. 1990, 93, 2537–2545. [Google Scholar] [CrossRef]
Tajti, A.; Szalay, P.G.; Császár, A.G.; Kállay, M.; Gauss, J.; Valeev, E.F.; Flowers, B.A.; Vázquez, J.; Stanton, J.F. HEAT: High accuracy extrapolated ab initio thermochemistry. J. Chem. Phys. 2004, 121, 11599–11613. [Google Scholar] [CrossRef]
Ochterski, J.W.; Petersson, G.A.; Montgomery Jr, J.A. A complete basis set model chemistry. V. Extensions to six or more heavy atoms. J. Chem. Phys. 1996, 104, 2598–2619. [Google Scholar] [CrossRef]
Montgomery Jr, J.A.; Ochterski, J.W.; Petersson, G.A. A complete basis set model chemistry. IV. An improved atomic pair natural orbital method. J. Chem. Phys. 1994, 101, 5900–5909. [Google Scholar] [CrossRef]
Helgaker, T.; Ruden, T.A.; Jørgensen, P.; Olsen, J.; Klopper, W. A priori calculation of molecular properties to chemical accuracy. J. Phys. Org. Chem. 2004, 17, 913–933. [Google Scholar] [CrossRef]
Sun, X.H.; Yao, Q.; Li, Z.R.; Wang, J.B.; Li, X.Y. Calculation of the rate constants for concerted elimination reaction class of hydroperoxyl-alkyl-peroxyl radicals. Theor. Chem. Acc. 2017, 136, 64. [Google Scholar] [CrossRef]
Sana, M.; Nguyen, M.T. Comment on the accurate theoretical determination of heats of formation. Chem. Phys. Lett. 1992, 196, 390–396. [Google Scholar] [CrossRef]
Curran, H.J.; Pitz, W.J.; Westbrook, C.K.; Dagaut, P.; Boettner, J.C.; Cathonnet, M. A wide range modeling study of dimethyl ether oxidation. Int. J. Chem. Kinet. 1998, 30, 229–241. [Google Scholar] [CrossRef]
Mokrushin, V.; Tsang, W. Chemrate v.1.5.8. National Institute of Standards and Technology: Gaithersburg, MD, USA, 2009. [Google Scholar]
Bugler, J.; Somers, K.P.; Silke, E.J.; Curran, H.J. Revisiting the kinetics and thermodynamics of the low-tempersture oxidation pathways of alkanes: A case study of the three pentane isomers. J. Phys. Chem. A 2015, 119, 7510–7527. [Google Scholar] [CrossRef] [PubMed]
Kerschgens, B.; Cai, L.; Pitsch, H.; Heuser, B. Di-n-buthylether, n-octanol, and n-octane as fuel candidates for diesel engine combustion. Combust. Flame 2016, 163, 66–78. [Google Scholar] [CrossRef]
Pu, J.F.; Yao, X.X.; Li, Z.R.; Li, X.Y. High-pressure limit rate rules for intramolecular H-migration reactions of α, β-hydroxyalkylperoxy radicals. Theor. Chem. Acc. 2021, 140, 147. [Google Scholar] [CrossRef]
Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, H.; et al. Gaussian 16, Revision, C. 01; Gaussian, Inc.: Wallingford, CT, USA, 2016. [Google Scholar]
Gonzalez, C.; Schlegel, H.B. Reaction path following in mass-weighted internal coordinates. J. Phys. Chem. 1990, 94, 5523–5527. [Google Scholar] [CrossRef]
Vandeputte, A.G.; Sabbe, M.K.; Reyniers, M.F.; Van Speybroeck, V.; Waroquier, M.; Marin, G.B. Theoretical Study of the Thermodynamics and Kinetics of Hydrogen Abstractions from Hydrocarbons. J. Phys. Chem. A 2007, 111, 11771–11786. [Google Scholar] [CrossRef]
Saeys, M.; Reyniers, M.F.; Marin, G.B.; Van Speybroeck, V.; Waroquier, M. Ab Initio Calculations for Hydrocarbons: Enthalpy of Formation, Transition State Geometry, and Activation Energy for Radical Reactions. J. Phys. Chem. A 2003, 107, 9147–9159. [Google Scholar] [CrossRef]
Truong, T.N. Reaction class transition state theory: Hydrogen abstraction reactions by hydrogen atoms as test cases. J. Chem. Phys. 2000, 113, 4957–4964. [Google Scholar] [CrossRef]
Ratkiewicz, A.; Huynh, L.K.; Truong, T.N. Performance of first-principles-based reaction class transition state theory. J. Phys. Chem. B 2016, 120, 1871–1884. [Google Scholar] [CrossRef]

Scheme 1. Reaction path for ethers at low temperature [3].

Scheme 2. The geometry structure of transition states (TSs) for different reaction classes. (a)1,3–H migration (b) 1,5–H migration (c)1,6-H migration (d)1,7–H migration.

Figure 1. Comparison of the energy barriers between the CCSD(T)/cc-PVTZ and CBS–QB3 method.

Figure 2. Comparison of the rate constants between the CBS–QB3 method and isodesmic reaction method at 500 K, 1000 K, and 1500 K.

Figure 3. Comparison of average rate constants for 1,6 H–migration reaction of different peroxy radicals [30,34]. (a) 1,6-H(p) (b) 1,6-H(s) (c) 1,6-H(t).

Figure 4. Rate constants of 1,6 H–migration reaction (1200 K) as a function of the type of carbon atoms of different peroxy radicals.

Table 1. Reaction equation lists.

Reaction Class	Reaction Equation
1,3-H(s)
R1	CH₃OCH₂OO·→ CH₃OC·HOOH
R2	CH₃CH₂OCH₂OO· → CH₃CH₂OC·HOOH
R3	CH₃(CH₂)₂OCH₂OO· → CH₃(CH₂)₂OC·HOOH
R4	CH₃(CH₂)₃OCH₂OO· → CH3(CH₂)₃OC·HOOH
1,3-H(t)
R5	CH₃OCH(CH₃)OO· → CH₃OC·(CH₃)OOH
R6	CH₃OCH(C₂H₅)OO· → CH₃OC·(C₂H₅)OOH
R7	CH₃CH₂OCH(CH₃)OO· → CH₃CH₂OC·(CH₃)OOH
R8	CH₃(CH₂)₂OCH(CH₃)OO· → CH₃(CH₂)₂OC·(CH₃)OOH
R9	CH₃(CH₂)₃OCH(CH₃)OO· → CH3(CH₂)₃OC·(CH₃)OOH
1,5-H(p)
R10	CH₃OCH₂OO· → ·CH₂OCH₂OOH
R11	CH₃OCH(CH₃)OO· → ·CH₂OCHCH(CH₃)OOH
R12	CH₃OCH(C₂H₅)OO· → ·CH₂OCH(C₂H₅)OOH
1,5-H(s)
R13	CH₃CH₂OCH₂OO· → CH₃C·HOCH₂OOH
R14	CH₃(CH₂)₂OCH₂OO·→ CH₃CH₂C·HOCH₂OOH
R15	CH₃CH₂OCH(CH₃)OO· → CH₃C·HOCH(CH₃)OOH
R16	CH₃CH₂OCH(C₂H₅)OO· → CH₃C·HOCH(C₂H₅)OOH
R17	CH₃CH(CH₃)CH₂OCH(CH₃)OO· → CH₃CH(CH₃)C·HOCH(CH₃)OOH
1,5-H(t)
R18	CH₃CH(CH₃)OCH₂OO·→ CH₃C·(CH₃)OCH₂OOH
R19	CH₃CH(CH₃)OCH(CH₃)OO·→ CH₃C·(CH₃)OCH(CH₃)OOH
1,6-H(p)
R20	CH₃CH₂OCH₂OO· → ·CH₂CH₂OCH₂OOH
R21	CH₃CH₂OCH(CH₃)OO· → ·CH₂CH₂OCH(CH₃)OOH
R22	CH₃CH₂OCH(C₂H₅)OO· → ·CH₂CH₂OCH(C₂H₅)OOH
1,6-H(s)
R23	CH₃(CH₂)₂OCH₂OO· → CH₃C·HCH₂OCH₂OOH
R24	CH₃(CH₂)₂OCH(CH₃)OO·· → CH₃C·HCH₂OCH(CH₃)OOH
R25	CH₃(CH₂)₃OCH₂OO· → CH₃CH₂C·HCH₂OCH₂OOH
R26	CH₃CH(CH₃)CH₂CH(CH₃)OCH(CH₃)OO· → CH₃CH(CH₃)C·HCH(CH₃)OCH(CH₃)OOH
1,6-H(t)
R27	CH₃CH(CH₃)CH₂OCH₂OO· → CH₃C·(CH₃)CH₂OCH₂OOH
R28	CH₃CH₂CH(CH₃)CH₂OCH₂OO· → CH₃CH₂C·(CH₃)CH₂OCH₂OOH
R29	CH₃CH₂CH(CH₃)CH(CH₃)OCH₂OO· → CH₃CH₂C·(CH₃)CH(CH₃)OCH₂OOH
R30	CH₃CH₂CH(CH₃)CH(CH₃)OCH(CH₃)OO· → CH₃CH₂C·(CH₃)CH(CH₃)OCH(CH₃)OOH
1,7-H(p)
R31	CH₃CH₂CH₂OCH₂OO· → ·CH₂CH₂CH₂OCH₂OOH
R32	CH₃CH(CH₃)CH₂OCH₂OO· → ·CH₂CH(CH₃)CH₂OCH₂OOH
R33	CH₃CH(CH₃)CH(CH₃)OCH₂OO· → ·CH₂CH(CH₃)CH(CH₃)OCH₂OOH
R34	CH₃CH(CH₃)CH(CH₃)OCH(CH₃)OO· → ·CH₂CH(CH₃)CH(CH₃)OCH(CH₃)OOH
1,7-H(s)
R35	CH₃(CH₂)₃OCH₂OO· → CH₃C·H(CH₂)₂OCH₂OOH
R36	CH₃CH₂CH(CH₃)CH₂OCH₂OO· → CH₃C·HCH(CH₃)CH₂OCH₂OOH
R37	CH₃CH₂CH(CH₃)CH(CH₃)OCH₂OO· → CH₃C·HCH(CH₃)CH(CH₃)OCH₂OOH
R38	CH₃CH₂CH(CH₃)CH(CH₃)OCH(CH₃)OO· → CH₃C·HCH(CH₃)CH(CH₃)OCH(CH₃)OOH
1,7-H(t)
R39	CH₃CH(CH₃)CH₂CH₂OCH₂OO· → CH₃C·(CH₃)(CH₂)₂OCH₂OOH
R40	CH₃CH(CH₃)CH(CH₃)CH₂OCH₂OO· → CH₃C·(CH₃)CH(CH₃)CH₂OCH₂OOH
R41	CH₃CH(CH₃)CH(CH₃)CH(CH₃)OCH₂OO· → CH₃C·(CH₃)CH(CH₃)CH(CH₃)OCH₂OOH

Table 2. The maximum deviation of the energy barriers for intramolecular H–migration reaction of R^IOR^IIOO· radicals (unit: kcal/mol).

Subclass	1,3-H	1,5-H	1,6-H	1,7-H
p	NA ^a	0.33	0.54	1.13
s	0.23	1.34	0.80	1.25
t	0.64	2.69	0.94	0.90
Avg ^b	1.16 *	7.19 *	6.85 *	6.70 *

The value with an asterisk in the last line is the maximum deviation of energy barriers for the reaction class. ^a NA represents there is no value for this reaction class. ^b Avg in the last line are the average energy barriers for reaction class.

Table 3. Comparison of the energy barriers by isodesmic reaction method with CBS–QB3 method (unit: kcal/mol).

Reaction Class	Reaction	CBS–QB3	IRM	ΔV
1,3-H(s)	R2	44.37	44.34	0.03
1,3-H(t)	R6	44.67	44.77	−0.10
1,5-H(p)	R11	23.08	23.02	0.06
1,5-H(s)	R14	21.98	22.21	−0.23
1,5-H(t)	R19	15.46	16.27	−0.81
1,6-H(p)	R21	27.39	27.40	−0.01
1,6-H(s)	R24	24.11	24.43	−0.32
1,6-H(t)	R28	20.04	20.74	−0.70
1,7-H(p)	R32	25.19	24.63	0.56
1,7-H(s)	R36	21.35	21.32	0.03
1,7-H(t)	R40	18.13	18.18	−0.05

ΔV represents the difference in energy barriers as determined by the isodesmic reaction method and the CBS–QB3 method.

Table 4. Comparison of the reaction enthalpies by isodesmic reaction method with CBS–QB3 method (unit: kcal/mol).

Reaction Class	Reaction	CBS–QB3	IRM	Curran et al. ^a [32]	ΔH ^b
1,3-H(s)	R2	−46.33	−46.26	NA	−0.07
1,3-H(t)	R6	−46.33	−46.60	NA	0.27
1,5-H(p)	R10	NA	9.77	9.30	NA
1,5-H(p)	R11	12.16	12.34	NA	−0.18
1,5-H(t)	R19	10.53	10.21	NA	0.32
1,6-H(p)	R21	17.08	17.08	NA	0.00
1,6-H(s)	R24	14.29	14.08	NA	0.21
1,6-H(t)	R28	11.29	11.43	NA	−0.14
1,7-H(p)	R32	15.38	15.42	NA	−0.04
1,7-H(s)	R36	12.54	12.28	NA	0.26
1,7-H(t)	R40	9.76	9.46	NA	0.30

^a The value of the reaction enthalpy from the work of Curran et al. [32] ^b ΔH represents the difference in reaction enthalpies as determined by the isodesmic reaction method and the CBS–QB3 method. NA represents that there is no value for this reaction class.

Table 5. Comparison of the rate constants at 500 K in this work with the analogous reactions in alkanes.

Reaction		R^IOR^IIOO·	ROO·
Reaction		This Work	Villano et al. ^a [7]	Bugler et al. ^b [34]
R10	COCOO· → ·COCOO	7.22	7.67 (1.06)	2.10 × 10¹ (2.91)
R13	CCOCOO· → CC·OCOO	3.67 × 10³	2.84 × 10² (12.92)	5.68 × 10² (6.46)
R18	C₂COCOO· → C₂C·OCOO	8.61 × 10⁵	1.61 × 10³ (5.35 × 10²)	8.22 × 10³ (1.05 × 10²)
R20	CCOCOO· → ·CCOCOO	6.76 × 10⁴	3.21 (2.11 × 10⁴)	6.36 (1.06 × 10⁴)

^a The ratios of average rate constants to the values of the corresponding average rate constants proposed by Villano et al. [7]. ^b The average rate constants to the values of the corresponding average rate constants proposed by Bugler et al. [34].

Table 6. Rate rules at high-pressure-limit for intramolecular H–migration reaction of R^IOR^IIOO· radicals.

Reaction Class	Modified Arrhenius Parameters			500 K	1000 K
Reaction Class	A (s⁻¹)	n	E (kcal/mol)	f	f
1,3-H(s)	1.36 × 10⁹	1.26	37.74	1.27	1.38
1,3-H(t)	9.58 × 10⁹	0.86	39.05	1.98	2.34
1,5-H(p)	2.03 × 10¹³	0.31	16.96	3.54 × 10³	2.43 × 10³
1,5-H(s)	4.47 × 10⁹	0.56	16.26	25.20	14.90
1,5-H(t)	1.22 × 10¹⁰	0.58	11.06	14.20	4.17
1,6-H(p)	1.19 × 10¹³	0.26	20.27	4.71 × 10³	2.63 × 10³
1,6-H(s)	2.15 × 10⁹	0.55	17.78	2.87	3.83
1,6-H(t)	2.62 × 10⁹	0.65	14.79	9.01	4.81
1,7-H(p)	7.37 × 10¹²	0.16	17.84	2.46 × 10³	1.59 × 10³
1,7-H(s)	2.32 × 10⁹	0.56	14.87	20.00	18.20
1,7-H(t)	2.19 × 10¹⁰	0.47	12.97	2.78	4.84

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Sun, X.; Li, Z. Reaction Rate Rules of Intramolecular H-Migration Reaction Class for R^IOR^IIOO·Radicals in Ether Combustion. Molecules 2024, 29, 4387. https://doi.org/10.3390/molecules29184387

AMA Style

Sun X, Li Z. Reaction Rate Rules of Intramolecular H-Migration Reaction Class for R^IOR^IIOO·Radicals in Ether Combustion. Molecules. 2024; 29(18):4387. https://doi.org/10.3390/molecules29184387

Chicago/Turabian Style

Sun, Xiaohui, and Zerong Li. 2024. "Reaction Rate Rules of Intramolecular H-Migration Reaction Class for R^IOR^IIOO·Radicals in Ether Combustion" Molecules 29, no. 18: 4387. https://doi.org/10.3390/molecules29184387

APA Style

Sun, X., & Li, Z. (2024). Reaction Rate Rules of Intramolecular H-Migration Reaction Class for R^IOR^IIOO·Radicals in Ether Combustion. Molecules, 29(18), 4387. https://doi.org/10.3390/molecules29184387

Article Menu

Reaction Rate Rules of Intramolecular H-Migration Reaction Class for R^IOR^IIOO·Radicals in Ether Combustion

Abstract

1. Introduction