High-Pressure Limit and Pressure-Dependent Rate Rules for β-Scission Reaction Class of Hydroperoxyl Alkyl Hydroperoxyl Radicals (•P(OOH)2) in Normal-Alkyl Cyclohexanes Combustion

Chemical kinetic studies of the β-scission reaction class of hydroperoxyl alkyl hydroperoxyl radicals (•P(OOH)2) from normal-alkyl cyclohexanes are carried out systematically through high-level ab initio calculations. Geometry optimizations and frequency calculations for all species involved in the reactions are performed at the B3LYP/CBSB7 level of theory. Electronic single-point energy calculations are calculated at the CBS-QB3 level of theory. Rate constants for the reactions of β-scission, in the temperature range of 500–1500 K and the pressure range of 0.01–100 atm, are calculated using transition state theory (TST) and Rice-Ramsberger-Kassel-Marcus/Master-Equation (RRKM/ME) theory taking asymmetric Eckart tunneling corrections and the one-dimensional hindered rotor approximation into consideration. The rate rules are obtained by averaging the rate constants of the representative reactions of this class. These rate rules can greatly assist in constructing more accurate low-temperature combustion mechanisms for normal-alkyl cyclohexanes.


Introduction
Hydrocarbons can be classified into several major groups, including alkanes, cycloalkanes (also known as naphthenes), aromatics, and alkenes [1].Cycloalkanes, especially those with C5 and C6 rings, are significant constituents of traditional transportation fuels like diesel and gasoline, accounting for approximately 10 to 30% of their composition [2,3].They also play a crucial role in jet fuels such as Jet-A/Jet-A1/JP-8 and RP-1, making up a considerable portion of their composition [4][5][6][7].Based on a comprehensive global study, it has been estimated that mono-, di-, and tri-cycloalkanes collectively constitute approximately 21% of the composition of jet fuel [8].Recently, there has been growing interest in cyclohexane and related compounds due to the impending depletion of conventional petroleum sources and the ever-increasing prices of fuels.Compared with non-cyclic counterparts such as alkanes, cycloalkanes are more prone to dehydrogenation due to their ring structure, resulting in the formation of unsaturated compounds such as cycloalkenes.These cycloalkenes can be further aromatized to form aromatics.Compared with non-aromatic compounds, aromatic compounds have a higher tendency to form soot during combustion.In addition, the existence of cyclic structures in cycloalkanes is conducive to the formation of polycyclic aromatic hydrocarbons (PAHs).PAHs are large, complex molecules made up of multiple carbon rings fused together.They are highly thermally stable and are considered to be significant contributors to soot emissions [9][10][11].To mitigate the production of unwanted pollutants and enhance the practical performance of engines, it is crucial to develop detailed combustion models that better represent low-temperature combustion techniques.These models play a crucial role in identifying critical reaction pathways and comprehending the low-temperature oxidation characteristics of large hydrocarbon fuels under different temperature and pressure conditions.This understanding is of immense significance in illuminating ignition phenomena within practical combustion systems like HCCI engines, jet engines, and spark ignition systems [12][13][14][15][16][17].
In the aspects of theoretical research, cycloalkanes [18][19][20][21] have received more attention in the past decade because they have the simplest molecular structure among all of the cyclic alkanes with a six-membered ring.The study of the oxidation reactivity of cyclohexanes with substituted alkyl-chains continues to be neglected, and very little information is available.There are only a few studies of the oxidation of methyl cyclohexane, ethyl cyclohexane, n-propyl cyclohexane, and n-butyl cyclohexane.In their study, Yang et al. [22] employed quantum chemical calculations, specifically CBS-QB3, to examine the uniqueness of the low-temperature oxidation of cycloalkanes.Their research focused on the influence of methyl substitution and analyzed the impacts of molecular structure, particularly in relation to potential energy surfaces (PES) involved in the isomerization reactions of ROO•, specifically 1,4 and 1,5 H-migration of cyclic alkylperoxyl radicals.Weber et al. [23] conducted calculations to determine the isomerization rate constants for the 2-methylcyclohexyl-1-peroxy radical, where the -OO group is attached at the beta position in relation to the methyl group.The researchers employed the CBS-QB3 composite method to calculate high-pressure limit rate constants, which were then utilized to update the kinetic model.Xing et al. [24,25] studied three different types of methyl cyclohexyl radicals, focusing on the low-temperature combustion reactions of these radicals.Furthermore, their work delved into the similarities and differences in reaction kinetics between acyclic and cyclic alkanes.Ning et al. [26] combined electron structure theory, transition state theory (TST), and Rice-Ramsberger-Kassel-Marcus/Master-Equation (RRKM/ME) theory to study the oxidation mechanism and kinetics of ethyl cyclohexane at low and medium temperatures.The high-pressure limit rate constants of tight transition state reactions are obtained from TST, and the high-pressure limit rate constants of barrierless reactions are obtained from the variational transition state theory (VTST).The rate constants for pressure-dependent reactions within the fall-off range were obtained using RRKM/ME theory.Through accurate computational analysis, valuable insights into the behavior of these reactions under varying temperature and pressure conditions were gained, thereby enhancing the low-temperature combustion model for cyclic alkanes.Liu et al. [27] conducted comprehensive chemical kinetic studies using high-level ab initio calculations to investigate hydrogen atom abstraction reactions by hydroperoxyl radicals from six alkyl cyclohexanes.Their systematic approach involves performing geometry optimizations and frequency calculations for all species involved in the reactions at the M06-2X/6-311++G(d,p) level of theory.Their approach also includes electronic single-point energy calculations, incorporating zero-point energy corrections, conducted at the UCCSD(T)-F12a/cc-pVDZ-F12 level of theory.They used the obtained kinetic and thermochemical data to update the mechanism of alkyl cyclohexane, and carried out further analysis to determine the impact of ignition delay time (in shock tube and fast compressor experiments) on the prediction.In our previous studies, we conducted systematic investigations into the isomerization reaction class of peroxyl alkyl radicals, the concerted elimination reaction class of peroxyl alkyl radicals, the β-scission reaction class of hydroperoxyl alkyl radicals, and the concerted elimination reaction class of peroxyl hydroperoxyl alkyl radicals derived from methyl, ethyl, n-propyl, and n-butyl cyclohexane.Through these studies, we developed high-pressure limit rate rules and pressure-dependent rate rules for these reactions.Ali et al. [28] examined the reactivity of n-butyl cyclohexane by employing the high-level quantum composite G3B3 method in conjunction with RRKM/ME simulations to elucidate reaction pathways.They obtained precise fuel ring-opening rate constants under high temperature conditions and extensively investigated seven distinct unimolecular decomposition channels for ring-opening.However, their study lacked theoretical kinetics and quantum chemical calculations at lower temperatures.Yang et al. [3] conducted a theoretical study on the H-migration reactions of cyclic alkylperoxy radicals and calculated the high-pressure limit rate constant and pressure-dependent rate constant of different reactions.All reactions are classified into seven reaction types, and these rules may be used in the development of low-temperature kinetic mechanisms for cycloalkanes.Generally speaking, as the simplest normal-alkyl cyclohexane, methyl cyclohexane has been widely studied in recent years.In contrast with methyl cyclohexane, the study on the combustion behavior of cyclohexanes with larger substituted alkyl-chains (the alkylic side chains from C2 to C4) remains too limited, and very little information is available.Therefore, increasing the credibility of a kinetic mechanism in jet fuels deserves much more effort, not only in experimental work but also in theoretical calculations.
The simplified kinetic scheme for the low-temperature oxidation of cycloalkanes, which is similar to alkanes, is illustrated in Scheme 1 [29].
Molecules 2024, 29, x FOR PEER REVIEW 3 of 25 opening rate constants under high temperature conditions and extensively investigated seven distinct unimolecular decomposition channels for ring-opening.However, their study lacked theoretical kinetics and quantum chemical calculations at lower temperatures.Yang et al. [3] conducted a theoretical study on the H-migration reactions of cyclic alkylperoxy radicals and calculated the high-pressure limit rate constant and pressuredependent rate constant of different reactions.All reactions are classified into seven reaction types, and these rules may be used in the development of low-temperature kinetic mechanisms for cycloalkanes.Generally speaking, as the simplest normal-alkyl cyclohexane, methyl cyclohexane has been widely studied in recent years.In contrast with methyl cyclohexane, the study on the combustion behavior of cyclohexanes with larger substituted alkyl-chains (the alkylic side chains from C2 to C4) remains too limited, and very li le information is available.Therefore, increasing the credibility of a kinetic mechanism in jet fuels deserves much more effort, not only in experimental work but also in theoretical calculations.
The simplified kinetic scheme for the low-temperature oxidation of cycloalkanes, which is similar to alkanes, is illustrated in Scheme 1 [29].
Scheme 1.The kinetic scheme for the low-temperature oxidation of cycloalkanes, with the specific class of reaction studied in this research highlighted in red [30].
Initially, alkane consumption takes place through H-atom abstraction, resulting in the formation of alkyl radicals (R•).These R• can react with the first molecular oxygen to generate alkyl peroxy radicals (ROO•), with the possibility of intermolecular hydrogen migration leading to the formation of hydroperoxyl alkyl radicals (•QOOH).Subsequently, •QOOH can react with the second molecular oxygen to produce hydroperoxylalkyl peroxyl radicals (•OOQOOH).Through intramolecular hydrogen migration, •OOQOOH can further transform into hydroperoxyl alkyl hydroperoxyl radicals (•P(OOH)2).In the low temperature regions, the decomposition of a •P(OOH)2 with a radical site beta to the -OOH moiety can result in the production of high molecular weight hydroperoxy olefin and HO2 radicals.This type of reaction is thought to be one of the reasons for the negative temperature coefficient (NTC) behavior of the fuel, as the reaction generates a relatively inactive HO2 radical, which has the potential to inhibit the ignition behavior of normal-alkyl cyclohexanes.
With the continuous advancement of engine technology and the introduction of new fuels, it remains crucial to effectively anticipate combustion behavior under various conditions of temperature, pressure, and equivalence ratio.This requires reliance on dependable combustion mechanisms and models that can accurately predict such behavior.The reliability of these combustion mechanisms depends on the comprehensiveness of the combustion reaction network as well as the precision of thermodynamic and kinetic parameters.However, precise kinetic parameters for various important reaction classes are difficult to obtain through experimental methods during the combustion process.Therefore, it becomes necessary to rely on advanced theoretical calculation methods to obtain rate constants for these crucial reactions.Internationally, numerous studies have been conducted to explore the formation of cyclohexenes and HO2 radicals through β-scission reactions of •QOOH derived from specific cycloalkanes.For instance, Xing et al. [24,25] focused on methyl cyclohexane and investigated potential energy surfaces, calculating pressure-dependent rate constants for the β-scission of •QOOH using the theoretical Scheme 1.The kinetic scheme for the low-temperature oxidation of cycloalkanes, with the specific class of reaction studied in this research highlighted in red [30].
Initially, alkane consumption takes place through H-atom abstraction, resulting in the formation of alkyl radicals (R•).These R• can react with the first molecular oxygen to generate alkyl peroxy radicals (ROO•), with the possibility of intermolecular hydrogen migration leading to the formation of hydroperoxyl alkyl radicals (•QOOH).Subsequently, •QOOH can react with the second molecular oxygen to produce hydroperoxylalkyl peroxyl radicals (•OOQOOH).Through intramolecular hydrogen migration, •OOQOOH can further transform into hydroperoxyl alkyl hydroperoxyl radicals (•P(OOH) 2 ).In the low temperature regions, the decomposition of a •P(OOH) 2 with a radical site beta to the -OOH moiety can result in the production of high molecular weight hydroperoxy olefin and HO 2 radicals.This type of reaction is thought to be one of the reasons for the negative temperature coefficient (NTC) behavior of the fuel, as the reaction generates a relatively inactive HO 2 radical, which has the potential to inhibit the ignition behavior of normal-alkyl cyclohexanes.
With the continuous advancement of engine technology and the introduction of new fuels, it remains crucial to effectively anticipate combustion behavior under various conditions of temperature, pressure, and equivalence ratio.This requires reliance on dependable combustion mechanisms and models that can accurately predict such behavior.The reliability of these combustion mechanisms depends on the comprehensiveness of the combustion reaction network as well as the precision of thermodynamic and kinetic parameters.However, precise kinetic parameters for various important reaction classes are difficult to obtain through experimental methods during the combustion process.Therefore, it becomes necessary to rely on advanced theoretical calculation methods to obtain rate constants for these crucial reactions.Internationally, numerous studies have been conducted to explore the formation of cyclohexenes and HO 2 radicals through β-scission reactions of •QOOH derived from specific cycloalkanes.For instance, Xing et al. [24,25] focused on methyl cyclohexane and investigated potential energy surfaces, calculating pressuredependent rate constants for the β-scission of •QOOH using the theoretical QCISD(T)/CBS level.Similarly, Ning et al. [26] examined the oxidation mechanisms and kinetics of ethyl cyclohexane at low and intermediate temperatures, and computed pressure-dependent rate constants for the β-scission of •QOOH using the CBS-QB3 level.In our previous work, we extensively investigated the potential energy surface of the oxidation pathway for the β-scission reactions of •QOOH in cyclohexanes during low-temperature combustion.Moreover, we accurately calculated the rate constants for this reaction class using the highly reliable CBS-QB3 theoretical level.Our findings strongly indicate that the β-scission reactions of •QOOH play a significant role in the oxidation pathways of cyclohexanes at low temperatures.However, research on the kinetics of the β-scission reactions of •P(OOH) 2 in cycloalkanes has been relatively limited.Consequently, in modeling studies of cycloalkane combustion, certain kinetic parameters are estimated by assuming similarities with reaction classes in alkanes.These approximations can introduce substantial errors in the modeling process.
Constructing a comprehensive combustion mechanism involves thousands of elementary reactions and hundreds of species.As hydrocarbon molecules increase in carbon atom count, the diversity and quantity of elementary reactions escalate exponentially, rendering manual construction of a detailed combustion mechanism impossible.A detailed hydrocarbon fuel mechanism typically comprises two primary components.First, a core mechanism encompasses smaller species (C0-C4) exhibiting high generality.Second, an extended mechanism caters to larger molecules (C5 and beyond) and is automatically generated by mechanism generation programs employing reaction class rules [31].This approach ensures an efficient and accurate representation of complex combustion processes.To the best of our knowledge, there is currently no existing literature that discusses rate rules for the essential β-scission process from •P(OOH) 2 of cyclohexanes.This process plays a critical role in constructing accurate low-and intermediate-temperature models for cycloalkanes.
Hence, the objective of this study is to comprehensively investigate the β-scission reactions of cyclohexanes to bridge the research gap in this field.Through theoretical calculations, we determine the rate constants at the high-pressure limit and a number of other pressures.Ultimately, this research aims to establish accurate high-pressure limit rate rules and the pressure-dependent rate rules for β-scission reactions of hydroperoxyl alkyl hydroperoxyl radicals in normal-alkyl cyclohexanes.These findings will serve as a theoretical foundation for constructing more comprehensive low-temperature combustion mechanisms for normal-alkyl cyclohexanes.

Results and Discussion
In our previous study, it was revealed that •QOOH generated via the 1.5 H-migration of ROO• from normal-alkyl cyclohexane is the primary intermediate in the H-migration reaction pathway, and this •QOOH species further facilitates O 2 addition.Consequently, the initial reactant configuration selected for the β-scission reaction of •P(OOH) 2 corresponding to the •QOOH configuration resulting from the 1.5 H-migration of ROO•.For our present investigation, a comprehensive set of 44 representative reactions involving methyl, ethyl, n-propyl, and n-butyl cyclohexane was chosen.These selected reactions are displayed in Scheme 2.

Conformational Analysis
In this study system, the significance of intramolecular hydrogen bonds between -OOH and -OOH groups in the reactants and transition states emphasizes the need to investigate their impact on conformational variations.Specifically, this study focuses on how intramolecular hydrogen bonds formed between the hydrogen atoms in one -OOH group and the oxygen atoms in another -OOH group affect the conformational behavior of the species.The influence of hydrogen bonding on the formation of ring structures plays a crucial role in diverse conformations and ring-to-ring conversions.This study focuses on three specific scenarios regarding the placement of -OOH and -OOH groups within a cycloalkane: (1) both groups located on the ring (Scheme 3a); (2) one group on the ring and the other on the alkyl side chain (Scheme 3b); (3) both groups positioned on the alkyl side chain (Scheme 3c).
of the species.The influence of hydrogen bonding on the formation of ring structures plays a crucial role in diverse conformations and ring-to-ring conversions.This study focuses on three specific scenarios regarding the placement of -OOH and -OOH groups within a cycloalkane: (1) both groups located on the ring (Scheme 3a); (2) one group on the ring and the other on the alkyl side chain (Scheme 3b); (3) both groups positioned on the alkyl side chain (Scheme 3c).In order to demonstrate the significance of intramolecular hydrogen bonding, we have chosen the reaction R17 as a representative example.Our aim is to study and analyze the effects of the presence or absence of intramolecular hydrogen bonding on the configurations of the reactants and the transition states, and further on the energy barriers associated with reactions, and the results are shown in Table 1.In Case (1), where neither the reactants nor the transition states contain intramolecular hydrogen bonds, the energy barrier is measured at 16.4 kcal mol −1 .Contrarily, in Case (2) of Table 1, the presence of intramolecular hydrogen bonds in the reactant leads to a reduction in energy, resulting in a higher energy barrier of 19.9 kcal mol −1 .However, in Case (3), we find that the presence of intramolecular hydrogen bonds in the transition state leads to a decrease in energy, resulting in a lower energy barrier of 12.7 kcal mol −1 .Finally, in Case (4), both the reactant and transition state contain intramolecular hydrogen bonds, and the energies of both the reactant and transition state decrease, resulting in a lower energy barrier of 16.2 kcal mol −1 compared to Case (1) in Table 1.However, an analysis of the intrinsic reaction coordinates (IRCs) for the four reactions listed in Table 1 reveals that only the reactant structure containing an intramolecular hydrogen bond can connect with the transition state that also In order to demonstrate the significance of intramolecular hydrogen bonding, we have chosen the reaction R17 as a representative example.Our aim is to study and analyze the effects of the presence or absence of intramolecular hydrogen bonding on the configurations of the reactants and the transition states, and further on the energy barriers associated with reactions, and the results are shown in Table 1.In Case (1), where neither the reactants nor the transition states contain intramolecular hydrogen bonds, the energy barrier is measured at 16.4 kcal mol −1 .Contrarily, in Case (2) of Table 1, the presence of intramolecular hydrogen bonds in the reactant leads to a reduction in energy, resulting in a higher energy barrier of 19.9 kcal mol −1 .However, in Case (3), we find that the presence of intramolecular hydrogen bonds in the transition state leads to a decrease in energy, resulting in a lower energy barrier of 12.7 kcal mol −1 .Finally, in Case (4), both the reactant and transition state contain intramolecular hydrogen bonds, and the energies of both the reactant and transition state decrease, resulting in a lower energy barrier of 16.2 kcal mol −1 compared to Case (1) in Table 1.However, an analysis of the intrinsic reaction coordinates (IRCs) for the four reactions listed in Table 1 reveals that only the reactant structure containing an intramolecular hydrogen bond can connect with the transition state that also contains an intramolecular hydrogen bond.Conversely, the reactant structure without an intramolecular hydrogen bond can only connect with a transition state that lacks an intramolecular hydrogen bond.As a result, only reactions (1) and (4) from Table 1 are observed in practical scenarios.Furthermore, due to the significant reduction in energy barriers offered by intramolecular hydrogen bonds, this study considers the presence of intramolecular hydrogen bonds in both the reactant and transition state structures.

Energy Barrier
In this study, the reaction barriers for all β-scission reactions of •P(OOH) 2 are calculated using the CBS-QB3 level and are documented in Table 2.At the same time, we calculate the average energy barrier and the maximum energy barrier deviation of the 44 reactions in this study, where the maximum energy barrier deviation represents the difference between the highest and lowest barriers, which are also listed in Table 2.It can be seen from Table 2 that the maximum energy barrier deviation is 2.5 kcal mol −1 , which indicates that for the β-scission reactions of •P(OOH) 2 , the difference in energy barriers of the 44 representative reactions selected is not particularly large.Therefore, a general reaction rate rule can be constructed for this reaction class.Table 1.A comparison of the energy barriers for the reaction R17 with and without hydrogen bond conformation (unit: kcal mol −1 ).

Reactants Transition States Energy Barriers
(1) R-no HB Molecules 2024, 29, x FOR PEER REVIEW 7 of 25 contains an intramolecular hydrogen bond.Conversely, the reactant structure without an intramolecular hydrogen bond can only connect with a transition state that lacks an intramolecular hydrogen bond.As a result, only reactions (1) and (4) from Table 1 are observed in practical scenarios.Furthermore, due to the significant reduction in energy barriers offered by intramolecular hydrogen bonds, this study considers the presence of intramolecular hydrogen bonds in both the reactant and transition state structures.
Table 1.A comparison of the energy barriers for the reaction R17 with and without hydrogen bond conformation (unit: kcal mol −1 ).

Energy Barrier
In this study, the reaction barriers for all β-scission reactions of •P(OOH)2 are calculated using the CBS-QB3 level and are documented in Table 2.At the same time, we calculate the average energy barrier and the maximum energy barrier deviation of the 44 reactions in this study, where the maximum energy barrier deviation represents the difference between the highest and lowest barriers, which are also listed in Table 2.It can be seen from Table 2 that the maximum energy barrier deviation is 2.5 kcal mol −1 , which indicates that for the β-scission reactions of •P(OOH)2, the difference in energy barriers of the 44 representative reactions selected is not particularly large.Therefore, a general reaction rate rule can be constructed for this reaction class.

TS-no HB
Molecules 2024, 29, x FOR PEER REVIEW 7 of 25 contains an intramolecular hydrogen bond.Conversely, the reactant structure without an intramolecular hydrogen bond can only connect with a transition state that lacks an intramolecular hydrogen bond.As a result, only reactions (1) and (4) from Table 1 are observed in practical scenarios.Furthermore, due to the significant reduction in energy barriers offered by intramolecular hydrogen bonds, this study considers the presence of intramolecular hydrogen bonds in both the reactant and transition state structures.
Table 1.A comparison of the energy barriers for the reaction R17 with and without hydrogen bond conformation (unit: kcal mol −1 ).

Energy Barrier
In this study, the reaction barriers for all β-scission reactions of •P(OOH)2 are calculated using the CBS-QB3 level and are documented in Table 2.At the same time, we calculate the average energy barrier and the maximum energy barrier deviation of the 44 reactions in this study, where the maximum energy barrier deviation represents the difference between the highest and lowest barriers, which are also listed in Table 2.It can be seen from Table 2 that the maximum energy barrier deviation is 2.5 kcal mol −1 , which indicates that for the β-scission reactions of •P(OOH)2, the difference in energy barriers of the 44 representative reactions selected is not particularly large.Therefore, a general reaction rate rule can be constructed for this reaction class.contains an intramolecular hydrogen bond.Conversely, the reactant structure without an intramolecular hydrogen bond can only connect with a transition state that lacks an intramolecular hydrogen bond.As a result, only reactions (1) and (4) from Table 1 are observed in practical scenarios.Furthermore, due to the significant reduction in energy barriers offered by intramolecular hydrogen bonds, this study considers the presence of intramolecular hydrogen bonds in both the reactant and transition state structures.
Table 1.A comparison of the energy barriers for the reaction R17 with and without hydrogen bond conformation (unit: kcal mol −1 ).

Energy Barrier
In this study, the reaction barriers for all β-scission reactions of •P(OOH)2 are calculated using the CBS-QB3 level and are documented in Table 2.At the same time, we calculate the average energy barrier and the maximum energy barrier deviation of the 44 reactions in this study, where the maximum energy barrier deviation represents the difference between the highest and lowest barriers, which are also listed in Table 2.It can be seen from Table 2 that the maximum energy barrier deviation is 2.5 kcal mol −1 , which indicates that for the β-scission reactions of •P(OOH)2, the difference in energy barriers of the 44 representative reactions selected is not particularly large.Therefore, a general reaction rate rule can be constructed for this reaction class.

TS-no HB
Molecules 2024, 29, x FOR PEER REVIEW 7 of 25 contains an intramolecular hydrogen bond.Conversely, the reactant structure without an intramolecular hydrogen bond can only connect with a transition state that lacks an intramolecular hydrogen bond.As a result, only reactions (1) and (4) from Table 1 are observed in practical scenarios.Furthermore, due to the significant reduction in energy barriers offered by intramolecular hydrogen bonds, this study considers the presence of intramolecular hydrogen bonds in both the reactant and transition state structures.
Table 1.A comparison of the energy barriers for the reaction R17 with and without hydrogen bond conformation (unit: kcal mol −1 ).

Energy Barrier
In this study, the reaction barriers for all β-scission reactions of •P(OOH)2 are calculated using the CBS-QB3 level and are documented in Table 2.At the same time, we calculate the average energy barrier and the maximum energy barrier deviation of the 44 reactions in this study, where the maximum energy barrier deviation represents the difference between the highest and lowest barriers, which are also listed in Table 2.It can be seen from Table 2 that the maximum energy barrier deviation is 2.5 kcal mol −1 , which indicates that for the β-scission reactions of •P(OOH)2, the difference in energy barriers of the 44 representative reactions selected is not particularly large.Therefore, a general reaction rate rule can be constructed for this reaction class.contains an intramolecular hydrogen bond.Conversely, the reactant structure without an intramolecular hydrogen bond can only connect with a transition state that lacks an intramolecular hydrogen bond.As a result, only reactions (1) and (4) from Table 1 are observed in practical scenarios.Furthermore, due to the significant reduction in energy barriers offered by intramolecular hydrogen bonds, this study considers the presence of intramolecular hydrogen bonds in both the reactant and transition state structures.
Table 1.A comparison of the energy barriers for the reaction R17 with and without hydrogen bond conformation (unit: kcal mol −1 ).

Energy Barrier
In this study, the reaction barriers for all β-scission reactions of •P(OOH)2 are calculated using the CBS-QB3 level and are documented in Table 2.At the same time, we calculate the average energy barrier and the maximum energy barrier deviation of the 44 reactions in this study, where the maximum energy barrier deviation represents the difference between the highest and lowest barriers, which are also listed in Table 2.It can be seen from Table 2 that the maximum energy barrier deviation is 2.5 kcal mol −1 , which indicates that for the β-scission reactions of •P(OOH)2, the difference in energy barriers of the 44 representative reactions selected is not particularly large.Therefore, a general reaction rate rule can be constructed for this reaction class.

TS-HB
Molecules 2024, 29, x FOR PEER REVIEW 7 of 25 contains an intramolecular hydrogen bond.Conversely, the reactant structure without an intramolecular hydrogen bond can only connect with a transition state that lacks an intramolecular hydrogen bond.As a result, only reactions (1) and (4) from Table 1 are observed in practical scenarios.Furthermore, due to the significant reduction in energy barriers offered by intramolecular hydrogen bonds, this study considers the presence of intramolecular hydrogen bonds in both the reactant and transition state structures.
Table 1.A comparison of the energy barriers for the reaction R17 with and without hydrogen bond conformation (unit: kcal mol −1 ).

Energy Barrier
In this study, the reaction barriers for all β-scission reactions of •P(OOH)2 are calculated using the CBS-QB3 level and are documented in Table 2.At the same time, we calculate the average energy barrier and the maximum energy barrier deviation of the 44 reactions in this study, where the maximum energy barrier deviation represents the difference between the highest and lowest barriers, which are also listed in Table 2.It can be seen from Table 2 that the maximum energy barrier deviation is 2.5 kcal mol −1 , which indicates that for the β-scission reactions of •P(OOH)2, the difference in energy barriers of the 44 representative reactions selected is not particularly large.Therefore, a general reaction rate rule can be constructed for this reaction class.contains an intramolecular hydrogen bond.Conversely, the reactant structure without an intramolecular hydrogen bond can only connect with a transition state that lacks an intramolecular hydrogen bond.As a result, only reactions (1) and (4) from Table 1 are observed in practical scenarios.Furthermore, due to the significant reduction in energy barriers offered by intramolecular hydrogen bonds, this study considers the presence of intramolecular hydrogen bonds in both the reactant and transition state structures.
Table 1.A comparison of the energy barriers for the reaction R17 with and without hydrogen bond conformation (unit: kcal mol −1 ).

Energy Barrier
In this study, the reaction barriers for all β-scission reactions of •P(OOH)2 are calculated using the CBS-QB3 level and are documented in Table 2.At the same time, we calculate the average energy barrier and the maximum energy barrier deviation of the 44 reactions in this study, where the maximum energy barrier deviation represents the difference between the highest and lowest barriers, which are also listed in Table 2.It can be seen from Table 2 that the maximum energy barrier deviation is 2.5 kcal mol −1 , which indicates that for the β-scission reactions of •P(OOH)2, the difference in energy barriers of the 44 representative reactions selected is not particularly large.Therefore, a general reaction rate rule can be constructed for this reaction class.

TS-HB
Molecules 2024, 29, x FOR PEER REVIEW 7 of 25 contains an intramolecular hydrogen bond.Conversely, the reactant structure without an intramolecular hydrogen bond can only connect with a transition state that lacks an intramolecular hydrogen bond.As a result, only reactions ( 1) and ( 4) from Table 1 are observed in practical scenarios.Furthermore, due to the significant reduction in energy barriers offered by intramolecular hydrogen bonds, this study considers the presence of intramolecular hydrogen bonds in both the reactant and transition state structures.
Table 1.A comparison of the energy barriers for the reaction R17 with and without hydrogen bond conformation (unit: kcal mol −1 ).

Energy Barrier
In this study, the reaction barriers for all β-scission reactions of •P(OOH)2 are calculated using the CBS-QB3 level and are documented in Table 2.At the same time, we calculate the average energy barrier and the maximum energy barrier deviation of the 44 reactions in this study, where the maximum energy barrier deviation represents the difference between the highest and lowest barriers, which are also listed in Table 2.It can be seen from Table 2 that the maximum energy barrier deviation is 2.5 kcal mol −1 , which indicates that for the β-scission reactions of •P(OOH)2, the difference in energy barriers of the 44 representative reactions selected is not particularly large.Therefore, a general reaction rate rule can be constructed for this reaction class.

Reaction
Energy Barrier

Comparison of Energy Barriers for β-Scission Class
In our previous studies, we have conducted systematic investigations into t sion reaction class of •QOOH derived from methyl, ethyl, n-propyl, and n-butyl ane, and we have computed the energy barriers for the β-scission reactions of •Q lizing the CBS-QB3 method.In order to verify whether intramolecular hydrogen sulting from the presence of two -OOH groups in •P(OOH)2 affect the energy barr β-scission reactions of •P(OOH)2, we compared those with the energy barriers of t sion reactions of •QOOH (reactive species containing one -OOH group).The r shown in Table 3, and the difference in the energy barriers between two classes is d

Comparison of Energy Barriers for β-Scission Class
In our previous studies, we have conducted systematic investigations into t sion reaction class of •QOOH derived from methyl, ethyl, n-propyl, and n-butyl ane, and we have computed the energy barriers for the β-scission reactions of •Q lizing the CBS-QB3 method.In order to verify whether intramolecular hydrogen sulting from the presence of two -OOH groups in •P(OOH)2 affect the energy barr β-scission reactions of •P(OOH)2, we compared those with the energy barriers of sion reactions of •QOOH (reactive species containing one -OOH group).The r shown in Table 3, and the difference in the energy barriers between two classes is d

Comparison of Energy Barriers for β-Scission Class
In our previous studies, we have conducted systematic investigations into t sion reaction class of •QOOH derived from methyl, ethyl, n-propyl, and n-butyl ane, and we have computed the energy barriers for the β-scission reactions of •Q lizing the CBS-QB3 method.In order to verify whether intramolecular hydrogen sulting from the presence of two -OOH groups in •P(OOH)2 affect the energy barri β-scission reactions of •P(OOH)2, we compared those with the energy barriers of t sion reactions of •QOOH (reactive species containing one -OOH group).The re shown in Table 3, and the difference in the energy barriers between two classes is d

Comparison of Energy Barriers for β-Scission Class
In our previous studies, we have conducted systematic investigations into t sion reaction class of •QOOH derived from methyl, ethyl, n-propyl, and n-butyl ane, and we have computed the energy barriers for the β-scission reactions of •Q lizing the CBS-QB3 method.In order to verify whether intramolecular hydrogen sulting from the presence of two -OOH groups in •P(OOH)2 affect the energy barr β-scission reactions of •P(OOH)2, we compared those with the energy barriers of t sion reactions of •QOOH (reactive species containing one -OOH group).The r shown in Table 3, and the difference in the energy barriers between two classes is d

Comparison of Energy Barriers for β-Scission Class
In our previous studies, we have conducted systematic investigations into the βscission reaction class of •QOOH derived from methyl, ethyl, n-propyl, and n-butyl cyclohexane, and we have computed the energy barriers for the β-scission reactions of •QOOH utilizing the CBS-QB3 method.In order to verify whether intramolecular hydrogen bonds resulting from the presence of two -OOH groups in •P(OOH) 2 affect the energy barriers of the β-scission reactions of •P(OOH) 2 , we compared those with the energy barriers of the β-scission reactions of •QOOH (reactive species containing one -OOH group).The results are shown in Table 3, and the difference in the energy barriers between two classes is discussed.

Comparison of Energy Barriers for β-Scission Class
In our previous studies, we have conducted systematic investigations into the β-scission reaction class of •QOOH derived from methyl, ethyl, n-propyl, and n-butyl cyclohexane, and we have computed the energy barriers for the β-scission reactions of •QOOH utilizing the CBS-QB3 method.In order to verify whether intramolecular hydrogen bonds resulting from the presence of two -OOH groups in •P(OOH)2 affect the energy barriers of the β-scission reactions of •P(OOH)2, we compared those with the energy barriers of the β-scission reactions of •QOOH (reactive species containing one -OOH group).The results are shown in Table 3, and the difference in the energy barriers between two classes is discussed.

Comparison of Energy Barriers for β-Scission Class
In our previous studies, we have conducted systematic investigations into the β-scission reaction class of •QOOH derived from methyl, ethyl, n-propyl, and n-butyl cyclohexane, and we have computed the energy barriers for the β-scission reactions of •QOOH utilizing the CBS-QB3 method.In order to verify whether intramolecular hydrogen bonds resulting from the presence of two -OOH groups in •P(OOH)2 affect the energy barriers of the β-scission reactions of •P(OOH)2, we compared those with the energy barriers of the β-scission reactions of •QOOH (reactive species containing one -OOH group).The results are shown in Table 3, and the difference in the energy barriers between two classes is discussed.

Comparison of Energy Barriers for β-Scission Class
In our previous studies, we have conducted systematic investigations into the β-scission reaction class of •QOOH derived from methyl, ethyl, n-propyl, and n-butyl cyclohexane, and we have computed the energy barriers for the β-scission reactions of •QOOH utilizing the CBS-QB3 method.In order to verify whether intramolecular hydrogen bonds resulting from the presence of two -OOH groups in •P(OOH)2 affect the energy barriers of the β-scission reactions of •P(OOH)2, we compared those with the energy barriers of the β-scission reactions of •QOOH (reactive species containing one -OOH group).The results are shown in Table 3, and the difference in the energy barriers between two classes is discussed.

Comparison of Energy Barriers for β-Scission Class
In our previous studies, we have conducted systematic investigations into the β-scission reaction class of •QOOH derived from methyl, ethyl, n-propyl, and n-butyl cyclohexane, and we have computed the energy barriers for the β-scission reactions of •QOOH utilizing the CBS-QB3 method.In order to verify whether intramolecular hydrogen bonds resulting from the presence of two -OOH groups in •P(OOH)2 affect the energy barriers of the β-scission reactions of •P(OOH)2, we compared those with the energy barriers of the β-scission reactions of •QOOH (reactive species containing one -OOH group).The results are shown in Table 3, and the difference in the energy barriers between two classes is discussed.

Comparison of Energy Barriers for β-Scission Class
In our previous studies, we have conducted systematic investigations into the β-scission reaction class of •QOOH derived from methyl, ethyl, n-propyl, and n-butyl cyclohexane, and we have computed the energy barriers for the β-scission reactions of •QOOH utilizing the CBS-QB3 method.In order to verify whether intramolecular hydrogen bonds resulting from the presence of two -OOH groups in •P(OOH)2 affect the energy barriers of the β-scission reactions of •P(OOH)2, we compared those with the energy barriers of the β-scission reactions of •QOOH (reactive species containing one -OOH group).The results are shown in Table 3, and the difference in the energy barriers between two classes is discussed.

Comparison of Energy Barriers for β-Scission Class
In our previous studies, we have conducted systematic investigations into the β-scission reaction class of •QOOH derived from methyl, ethyl, n-propyl, and n-butyl cyclohexane, and we have computed the energy barriers for the β-scission reactions of •QOOH utilizing the CBS-QB3 method.In order to verify whether intramolecular hydrogen bonds resulting from the presence of two -OOH groups in •P(OOH)2 affect the energy barriers of the β-scission reactions of •P(OOH)2, we compared those with the energy barriers of the β-scission reactions of •QOOH (reactive species containing one -OOH group).The results are shown in Table 3, and the difference in the energy barriers between two classes is discussed.

Comparison of Energy Barriers for β-Scission Class
In our previous studies, we have conducted systematic investigations into the β-scission reaction class of •QOOH derived from methyl, ethyl, n-propyl, and n-butyl cyclohexane, and we have computed the energy barriers for the β-scission reactions of •QOOH utilizing the CBS-QB3 method.In order to verify whether intramolecular hydrogen bonds resulting from the presence of two -OOH groups in •P(OOH)2 affect the energy barriers of the β-scission reactions of •P(OOH)2, we compared those with the energy barriers of the β-scission reactions of •QOOH (reactive species containing one -OOH group).The results are shown in Table 3, and the difference in the energy barriers between two classes is discussed.

Comparison of Energy Barriers for β-Scission Class
In our previous studies, we have conducted systematic investigations into the β-scission reaction class of •QOOH derived from methyl, ethyl, n-propyl, and n-butyl cyclohexane, and we have computed the energy barriers for the β-scission reactions of •QOOH utilizing the CBS-QB3 method.In order to verify whether intramolecular hydrogen bonds resulting from the presence of two -OOH groups in •P(OOH)2 affect the energy barriers of the β-scission reactions of •P(OOH)2, we compared those with the energy barriers of the β-scission reactions of •QOOH (reactive species containing one -OOH group).The results are shown in Table 3, and the difference in the energy barriers between two classes is discussed.

Comparison of Energy Barriers for β-Scission Class
In our previous studies, we have conducted systematic investigations into the β-scission reaction class of •QOOH derived from methyl, ethyl, n-propyl, and n-butyl cyclohexane, and we have computed the energy barriers for the β-scission reactions of •QOOH utilizing the CBS-QB3 method.In order to verify whether intramolecular hydrogen bonds resulting from the presence of two -OOH groups in •P(OOH)2 affect the energy barriers of the β-scission reactions of •P(OOH)2, we compared those with the energy barriers of the β-scission reactions of •QOOH (reactive species containing one -OOH group).The results are shown in Table 3, and the difference in the energy barriers between two classes is discussed.

Comparison of Energy Barriers for β-Scission Class
In our previous studies, we have conducted systematic investigations into the β-scission reaction class of •QOOH derived from methyl, ethyl, n-propyl, and n-butyl cyclohexane, and we have computed the energy barriers for the β-scission reactions of •QOOH utilizing the CBS-QB3 method.In order to verify whether intramolecular hydrogen bonds resulting from the presence of two -OOH groups in •P(OOH)2 affect the energy barriers of the β-scission reactions of •P(OOH)2, we compared those with the energy barriers of the β-scission reactions of •QOOH (reactive species containing one -OOH group).The results are shown in Table 3, and the difference in the energy barriers between two classes is discussed.

Comparison of Energy Barriers for β-Scission Class
In our previous studies, we have conducted systematic investigations into the β-scission reaction class of •QOOH derived from methyl, ethyl, n-propyl, and n-butyl cyclohexane, and we have computed the energy barriers for the β-scission reactions of •QOOH utilizing the CBS-QB3 method.In order to verify whether intramolecular hydrogen bonds resulting from the presence of two -OOH groups in •P(OOH)2 affect the energy barriers of the β-scission reactions of •P(OOH)2, we compared those with the energy barriers of the β-scission reactions of •QOOH (reactive species containing one -OOH group).The results are shown in Table 3, and the difference in the energy barriers between two classes is discussed.

Comparison of Energy Barriers for β-Scission Class
In our previous studies, we have conducted systematic investigations into the β-scission reaction class of •QOOH derived from methyl, ethyl, n-propyl, and n-butyl cyclohexane, and we have computed the energy barriers for the β-scission reactions of •QOOH utilizing the CBS-QB3 method.In order to verify whether intramolecular hydrogen bonds resulting from the presence of two -OOH groups in •P(OOH)2 affect the energy barriers of the β-scission reactions of •P(OOH)2, we compared those with the energy barriers of the β-scission reactions of •QOOH (reactive species containing one -OOH group).The results are shown in Table 3, and the difference in the energy barriers between two classes is discussed.

Comparison of Energy Barriers for β-Scission Class
In our previous studies, we have conducted systematic investigations into the β-scission reaction class of •QOOH derived from methyl, ethyl, n-propyl, and n-butyl cyclohexane, and we have computed the energy barriers for the β-scission reactions of •QOOH utilizing the CBS-QB3 method.In order to verify whether intramolecular hydrogen bonds resulting from the presence of two -OOH groups in •P(OOH)2 affect the energy barriers of the β-scission reactions of •P(OOH)2, we compared those with the energy barriers of the β-scission reactions of •QOOH (reactive species containing one -OOH group).The results are shown in Table 3, and the difference in the energy barriers between two classes is discussed.

Comparison of Energy Barriers for β-Scission Class
In our previous studies, we have conducted systematic investigations into the β-scission reaction class of •QOOH derived from methyl, ethyl, n-propyl, and n-butyl cyclohexane, and we have computed the energy barriers for the β-scission reactions of •QOOH utilizing the CBS-QB3 method.In order to verify whether intramolecular hydrogen bonds resulting from the presence of two -OOH groups in •P(OOH)2 affect the energy barriers of the β-scission reactions of •P(OOH)2, we compared those with the energy barriers of the β-scission reactions of •QOOH (reactive species containing one -OOH group).The results are shown in Table 3, and the difference in the energy barriers between two classes is discussed.As can be seen from Table 3, when the radical sites are located on the primary and secondary carbon, the energy barriers of the β-scission reactions of •QOOH is basically similar to or lower than those of •P(OOH)2.When the radical sites are located on the tertiary carbon, the energy barriers of the β-scission reactions of •QOOH are higher than those of •P(OOH)2.However, when the reaction centers of the β-scission reactions occur As can be seen from Table 3, when the radical sites are located on the primary and secondary carbon, the energy barriers of the β-scission reactions of •QOOH is basically similar to or lower than those of •P(OOH)2.When the radical sites are located on the tertiary carbon, the energy barriers of the β-scission reactions of •QOOH are higher than those of •P(OOH)2.However, when the reaction centers of the β-scission reactions occur As can be seen from Table 3, when the radical sites are located on the primary and secondary carbon, the energy barriers of the β-scission reactions of •QOOH is basically similar to or lower than those of •P(OOH)2.When the radical sites are located on the tertiary carbon, the energy barriers of the β-scission reactions of •QOOH are higher than those of •P(OOH)2.However, when the reaction centers of the β-scission reactions occur As can be seen from Table 3, when the radical sites are located on the primary and secondary carbon, the energy barriers of the β-scission reactions of •QOOH is basically similar to or lower than those of •P(OOH)2.When the radical sites are located on the tertiary carbon, the energy barriers of the β-scission reactions of •QOOH are higher than those of •P(OOH)2.However, when the reaction centers of the β-scission reactions occur As can be seen from Table 3, when the radical sites are located on the primary and secondary carbon, the energy barriers of the β-scission reactions of •QOOH is basically similar to or lower than those of •P(OOH)2.When the radical sites are located on the tertiary carbon, the energy barriers of the β-scission reactions of •QOOH are higher than those of •P(OOH)2.However, when the reaction centers of the β-scission reactions occur As can be seen from Table 3, when the radical sites are located on the primary and secondary carbon, the energy barriers of the β-scission reactions of •QOOH is basically similar to or lower than those of •P(OOH)2.When the radical sites are located on the tertiary carbon, the energy barriers of the β-scission reactions of •QOOH are higher than those of •P(OOH)2.However, when the reaction centers of the β-scission reactions occur As can be seen from Table 3, when the radical sites are located on the primary and secondary carbon, the energy barriers of the β-scission reactions of •QOOH is basically similar to or lower than those of •P(OOH)2.When the radical sites are located on the tertiary carbon, the energy barriers of the β-scission reactions of •QOOH are higher than those of •P(OOH)2.However, when the reaction centers of the β-scission reactions occur As can be seen from Table 3, when the radical sites are located on the primary and secondary carbon, the energy barriers of the β-scission reactions of •QOOH is basically similar to or lower than those of •P(OOH)2.When the radical sites are located on the tertiary carbon, the energy barriers of the β-scission reactions of •QOOH are higher than those of •P(OOH)2.However, when the reaction centers of the β-scission reactions occur As can be seen from Table 3, when the radical sites are located on the primary and secondary carbon, the energy barriers of the β-scission reactions of •QOOH is basically similar to or lower than those of •P(OOH) 2 .When the radical sites are located on the tertiary carbon, the energy barriers of the β-scission reactions of •QOOH are higher than those of •P(OOH) 2 .However, when the reaction centers of the β-scission reactions occur on the ring, there are some β-scission reactions of •QOOH that have lower energy barriers than those of •P(OOH) 2 , while others have higher energy barriers than those of •P(OOH) 2 .These comparison results show that the presence of intramolecular hydrogen bonds in the reactants and transition states will affect the energy barrier of the reactions with similar reaction classes.

Rate Constants and Rate Rules
In this study, rate constants for the investigated reactions are calculated within the temperature range of 500 K to 1500 K.These rate constants are then fitted to the modified Arrhenius Equation: Which is given in the form (A, n, E).To derive the rate rule for a class, the rate constant of each reaction at each temperature is fitted to Equation (2): where k T avg denotes the average rate constant at temperature T, the rate constant for the i reaction at temperature T is denoted by k T i , and N represents the total number of reactions within a particular class.By utilizing this method, a three-parameter (A, n, E) representation for the rate rule of the class is obtained through fitting to the modified Arrhenius Equation.
To measure the uncertainty and the applicability of the rate rule for a class, an uncertainty factor f of a reaction class is defined as: f = k max /k min , where k max and k min are the maximum and the minimum rate constant in a class, respectively.k/k ave , which is defined as the ratio of the rate constant of a reaction to the average rate constant of the reactions in a class, are also calculated at a given temperature to measure the suitability of the rate rule.

High-Pressure Limit Rate Constants and Rate Rules
The calculated high-pressure limit rate constants at temperatures from 500-1500 K for all studied reactions are given in the form of (A, n, E) in Table 4.The f = k max /k min and k/k ave in each reaction subclass are also given in these Table 4. Observing Table 4, it is apparent that the ratios between reaction rate constants and the average rate constant for the class, along with the uncertainty factor for the class, fall within a range of one order of magnitude.This signifies that the high-pressure limit rate rules formulated for the β-scission class are both reasonable and valid.Zou et al. [32] successfully developed a comprehensive combustion mechanism for ethyl cyclohexane at low temperature.Within this mechanism, the high-pressure limit rate constants for the β-scission reactions of •P(OOH) 2 were determined using kinetic data from similar reaction classes.Specifically, these rate constants were derived from the computational results of Fernandes et al. [21] for cyclohexane and Xing et al. [24,25] for methylcyclohexane.Recently, Liu et al. [27] developed an intricate chemical kinetic model that encompasses the entire combustion process of n-propyl cyclohexane, spanning from low to high temperatures.Within this comprehensive mechanism, the rate constants for the β-scission reactions of •P(OOH) 2 were determined based on the kinetic data obtained from the β-scission reactions of •QOOH with similar structural characteristics.When the reaction centers reside on the alkyl side chain of n-propyl cyclohexane, the rate constants of the reactions in the mechanism are derived from the kinetic data obtained from βscission reactions occurring within analogous alkane structures [33].Conversely, when the reaction center is situated within the cyclic structure, the rate constants of the reactions in the mechanism primarily rely on theoretical kinetic studies of cyclohexane with similar structural characteristics [21].Mao et al. [34,35] have devised mechanisms to replicate the ignition and oxidation of n-butyl cyclohexane across diverse operational scenarios.Nevertheless, the detailed kinetic models they proposed omitted the β-scission reactions of •P(OOH) 2 within n-butyl cyclohexane.Figure 1 presents a comparative analysis of rate constants between our calculated values for high-pressure limit reactions and those from the models developed by Zou et al. for ethyl cyclohexane and Liu et al. for n-propyl cyclohexane.These comparisons were conducted using a temperature range of 500 K to 1500 K. Figure 1 clearly illustrates that the rate constants derived from the β-scission of •P(OOH) 2 rate rule in our study consistently surpass the values of similar reactions found in the mechanisms.This discrepancy is evident with ratio ranges of 14.9 to 44.7, 6.3 to 14.4, and 15.5 to 21.5, respectively.These findings highlight a substantial deviation between our calculated rate constants and those reported in the mechanisms by Zou et  •P(OOH)2 rate rule in our study consistently surpass the values of similar reactions found in the mechanisms.This discrepancy is evident with ratio ranges of 14.9 to 44.7, 6.3 to 14.4, and 15.5 to 21.5, respectively.These findings highlight a substantial deviation between our calculated rate constants and those reported in the mechanisms by Zou et al., and Liu et al.Villano et al. [33] systematically investigated the β-scission reactions of C2-C6 hydroperoxy alkyl radicals in non-cyclic systems using the CBS-QB3 level of theory.In their study, they assigned a single rate rule for the β-scission reaction within this reaction class.Miyoshi [36] reported high-pressure limit rate rules for the β-scission class in non-cyclic systems at the CBS-QB3 level of theory.In his study, the β-scission class was further categorized into nine distinct subclasses based on the specific carbon atom bonded to the hydroperoxy group and the characteristics of the radical site.For each subclass, Miyoshi developed a unique rate rule.Figure 2 provides comparisons with the literature.In Figure Villano et al. [33] systematically investigated the β-scission reactions of C2-C6 hydroperoxy alkyl radicals in non-cyclic systems using the CBS-QB3 level of theory.In their study, they assigned a single rate rule for the β-scission reaction within this reac-tion class.Miyoshi [36] reported high-pressure limit rate rules for the β-scission class in non-cyclic systems at the CBS-QB3 level of theory.In his study, the β-scission class was further categorized into nine distinct subclasses based on the specific carbon atom bonded to the hydroperoxy group and the characteristics of the radical site.For each subclass, Miyoshi developed a unique rate rule.Figure 2 provides comparisons with the literature.In Figure 2a, the rate constants of the β-scission of •P(OOH) 2 rate rule derived in this study are consistently higher than those obtained from the β-scission of QOOH rate rule in the study conducted by Villano et al., with a ratio ranging from 7.8 to 12.2.Similarly, Figure 2b illustrates that the rate constants of the β-scission of •P(OOH) 2 rate rule derived in this study are all greater than the values obtained from the β-scission of QOOH subclasses rate rules in the study conducted by Miyoshi et   (a) the reported rules by Villano et al. [33]; (b) the reported rules by Miyoshi [36].

Pressure-Dependent Rate Constants and Rate Rules
The β-scission of •P(OOH)2 is classified as a pressure-dependent channel.To determine the rate constants for this reaction, the RRKM/ME theory is employed.The pressuredependent rate constants are calculated over a range of pressures (0.01-100 atm) and temperatures (500-1500 K).The calculated pressure-dependent rate constants are presented in Table 5, denoted in the form of (A, n, E).The simple rate rule method is utilized to determine the pressure-dependent rate rules.This method entails averaging the pressuredependent rate constants for all reactions within each subclass.The resulting pressuredependent rate rules are fi ed into the form of (A, n, E) and presented in Table 5.To enable easier comparisons between different reactions or rate rules, Table 4 also includes the pressure-dependent rate constants at 800 K.The ratios of the rate constants for each reaction to the average rate constant within this class are presented in Table 5, providing a visual representation of the deviation of the rate constants within the class.To evaluate the uncertainty of the rate rule, we define the uncertainty factor as the ratio of the maximum and minimum pressure-dependent rate constants within the class.The uncertainty factors for this class are also listed in Table 5.Based on the findings presented in Table 5, it is evident that the ratios between the rate constants for reactions and the average rate constants within the class, along with the uncertainty factors for the class, fall within the range of 1 to 2 orders of magnitude.This result implies that the pressure-dependent rate rules formulated for the β-scission of •P(OOH)2 are both reasonable and satisfactory.(a) the reported rules by Villano et al. [33]; (b) the reported rules by Miyoshi [36].

Pressure-Dependent Rate Constants and Rate Rules
The β-scission of •P(OOH) 2 is classified as a pressure-dependent channel.To determine the rate constants for this reaction, the RRKM/ME theory is employed.The pressure-dependent rate constants are calculated over a range of pressures (0.01-100 atm) and temperatures (500-1500 K).The calculated pressure-dependent rate constants are presented in Table 5, denoted in the form of (A, n, E).The simple rate rule method is utilized to determine the pressure-dependent rate rules.This method entails averaging the pressure-dependent rate constants for all reactions within each subclass.The resulting pressure-dependent rate rules are fitted into the form of (A, n, E) and presented in Table 5.
To enable easier comparisons between different reactions or rate rules, Table 4 also includes the pressure-dependent rate constants at 800 K.The ratios of the rate constants for each reaction to the average rate constant within this class are presented in Table 5, providing a visual representation of the deviation of the rate constants within the class.To evaluate the uncertainty of the rate rule, we define the uncertainty factor as the ratio of the maximum and minimum pressure-dependent rate constants within the class.The uncertainty factors for this class are also listed in Table 5.Based on the findings presented in Table 5, it is evident that the ratios between the rate constants for reactions and the average rate constants within the class, along with the uncertainty factors for the class, fall within the range of 1 to 2 orders of magnitude.This result implies that the pressure-dependent rate rules formulated for the β-scission of •P(OOH) 2 are both reasonable and satisfactory.In this study, we performed a comparative analysis between the pressure-dependent rate rule for the β-scission reactions of •P(OOH) 2 obtained in our research, and the values derived from the combustion mechanisms proposed by Zou et al. [32] for ethyl cyclohexane and Liu et al. [27] for n-propyl cyclohexane at atmospheric pressure.The results are illustrated in Figure 3a and b, respectively.The results presented in Figure 3a clearly demonstrate that the pressure-dependent rate constants calculated in this study are consistently higher than the corresponding values obtained from the combustion mechanism, with ratios ranging from 12.4 to 131.2, 2.3 to 61.4, and 11.8 to 61.3.In Liu's propyl cyclohexane mechanism, all pressure-dependent rate constants associated with reactions occurring on the rings are based on the same set of kinetic data, derived from the rate constants of β-scission reactions observed in similar reaction types in cyclohexane.Therefore, in Figure 3b, we solely focus on comparing the pressure-dependent rate rule obtained in this study at 1 atm with the rate constants of one of these reactions.Based on Figure 3b, it is evident that the pressure-dependent rate rule calculated in our study for the β-scission reaction surpasses the corresponding values present in the mechanism, particularly in the lower temperature range.In fact, the ratios can reach almost 60 times at 500 K.
with ratios ranging from 12.4 to 131.2, 2.3 to 61.4, and 11.8 to 61.3.In Liu's propyl cyclohexane mechanism, all pressure-dependent rate constants associated with reactions occurring on the rings are based on the same set of kinetic data, derived from the rate constants of β-scission reactions observed in similar reaction types in cyclohexane.Therefore, in Figure 3b, we solely focus on comparing the pressure-dependent rate rule obtained in this study at 1 atm with the rate constants of one of these reactions.Based on Figure 3b, it is evident that the pressure-dependent rate rule calculated in our study for the β-scission reaction surpasses the corresponding values present in the mechanism, particularly in the lower temperature range.In fact, the ratios can reach almost 60 times at 500 K.

Comparison of the Pressure-Dependent Rate Rules at Different Pressures
The rate constants for the β-scission reactions of •P(OOH)2 are depicted in Figure 4a, showcasing their dependence on temperatures across the range of 500-1500 K. Figure 4b exhibits the rate constants for this reaction class as a function of pressures, spanning from 0.01 to 100 atm.Notably, Figure 4a reveals that the β-scission reactions are highly sensitive to pressure, with the impact of pressure on the rate constants intensifying with increasing temperature.Particularly, this influence becomes significant at temperatures surpassing 700 K.By examining Figure 4b, it becomes evident that reactions belonging to the class of β-scission of •P(OOH)2 reach their high-pressure limit rate constants at varying pressures at different temperatures.When the temperature is relatively low, for instance, 500 K, the reaction rate constants swiftly reach the high-pressure limit at around 10 atm.However, as the temperature escalates to 800 K, the reaction rate constants exhibit continuous growth across a wide pressure range, spanning from 0.01 atm to 100 atm.Astonishingly, even at the highest pressure of 100 atm, the rate constants have not yet achieved the highpressure limit.Based on the aforementioned comparison, it is evident that the rate constants associated with the β-scission of •P(OOH)2 exhibit significant dependence on pressure.This implies that the pressure effect of the β-scission reactions within this class cannot be disregarded, particularly at elevated temperatures.

Comparison of the Pressure-Dependent Rate Rules at Different Pressures
The rate constants for the β-scission reactions of •P(OOH) 2 are depicted in Figure 4a, showcasing their dependence on temperatures across the range of 500-1500 K. Figure 4b exhibits the rate constants for this reaction class as a function of pressures, spanning from 0.01 to 100 atm.Notably, Figure 4a reveals that the β-scission reactions are highly sensitive to pressure, with the impact of pressure on the rate constants intensifying with increasing temperature.Particularly, this influence becomes significant at temperatures surpassing 700 K.By examining Figure 4b, it becomes evident that reactions belonging to the class of β-scission of •P(OOH) 2 reach their high-pressure limit rate constants at varying pressures at different temperatures.When the temperature is relatively low, for instance, 500 K, the reaction rate constants swiftly reach the high-pressure limit at around 10 atm.However, as the temperature escalates to 800 K, the reaction rate constants exhibit continuous growth across a wide pressure range, spanning from 0.01 atm to 100 atm.Astonishingly, even at the highest pressure of 100 atm, the rate constants have not yet achieved the high-pressure limit.Based on the aforementioned comparison, it is evident that the rate constants associated with the β-scission of •P(OOH) 2 exhibit significant dependence on pressure.This implies that the pressure effect of the β-scission reactions within this class cannot be disregarded, particularly at elevated temperatures.

Electronic Structure Calculation
The energies of reactants, products, intermediates, and transition states (TS) involved in the reactions are determined using the CBS-QB3 composite method [37] implemented in the Gaussian 16 quantum chemistry software package [38].In the CBS-QB3 composite method, the geometry optimization and frequency analysis are carried out at the B3LYP/CBSB7 level, and the single-point energy calculations are performed at the CCSD(T)/6−31+G(d') and MP4SDQ/CBSB4 levels.The total energy is extrapolated to the infinite-basis-set limit using pair natural-orbital energies at the MP2/CBSB3 level and an additive correction to the CCSD(T) level.A correction for spin contamination in open-shell species is added to the total energy.In Supplemental Material-SI Table S1, the comparison between the spin eigenvalues with the S 2 computed by the Gaussian 16 package for the species involved in this study is given.It is worth noting that the maximum relative error observed is 4.3%.As a result, the spin contamination involved in this study can be safely ignored.To validate the reaction pathways and the connection between reactants and products, we performed an intrinsic reaction coordinate (IRC) analysis [39] at the B3LYP/CBSB7 level.The representative IRCs involving transition states for β-scission of •P(OOH)2 reactions are provided in the Figure S1 of the Supplemental Material-SI.The optimized structures of all species in this study can be accessed in Supplemental Material-SI, including their corresponding Cartesian coordinates.

Calculation of Rate Constant
In this study, we utilized the Chemrate program [40] to compute the high-pressure limit rate constants and pressure-dependent rate constants for all reactions.The computation of these rate constants is based on the principles of traditional TST [41,42].The expression for these rate constants is as follows: The expression provided includes several important parameters.(), ,  , and ℎ represent the Eckart tunneling correction factor, reaction path degeneracy, Bol mann constant, and Planck constant, respectively.The parameter T signifies the temperature, and  indicates the magnitude of the energy barrier.Furthermore,  () and  () represent the partition functions of the transition state and the reactants, respectively.

Electronic Structure Calculation
The energies of reactants, products, intermediates, and transition states (TS) involved in the reactions are determined using the CBS-QB3 composite method [37] implemented in the Gaussian 16 quantum chemistry software package [38].In the CBS-QB3 composite method, the geometry optimization and frequency analysis are carried out at the B3LYP/CBSB7 level, and the single-point energy calculations are performed at the CCSD(T)/6−31+G(d') and MP4SDQ/CBSB4 levels.The total energy is extrapolated to the infinite-basis-set limit using pair natural-orbital energies at the MP2/CBSB3 level and an additive correction to the CCSD(T) level.A correction for spin contamination in open-shell species is added to the total energy.In Supplemental Material-SI Table S1, the comparison between the spin eigenvalues with the S 2 computed by the Gaussian 16 package for the species involved in this study is given.It is worth noting that the maximum relative error observed is 4.3%.As a result, the spin contamination involved in this study can be safely ignored.To validate the reaction pathways and the connection between reactants and products, we performed an intrinsic reaction coordinate (IRC) analysis [39] at the B3LYP/CBSB7 level.The representative IRCs involving transition states for β-scission of •P(OOH) 2 reactions are provided in the Figure S1 of the Supplemental Material-SI.The optimized structures of all species in this study can be accessed in Supplemental Material-SI, including their corresponding Cartesian coordinates.

Calculation of Rate Constant
In this study, we utilized the Chemrate program [40] to compute the high-pressure limit rate constants and pressure-dependent rate constants for all reactions.The computation of these rate constants is based on the principles of traditional TST [41,42].The expression for these rate constants is as follows: The expression provided includes several important parameters.κ(T), rpd, κ B , and h represent the Eckart tunneling correction factor, reaction path degeneracy, Boltzmann constant, and Planck constant, respectively.The parameter T signifies the temperature, and V ̸ = indicates the magnitude of the energy barrier.Furthermore, Q TS (T) and Q R (T) represent the partition functions of the transition state and the reactants, respectively.
In this study, we demonstrate that the β-scission reactions arising from •P(OOH) 2 exhibit pressure-dependent characteristics.The rate constants' sensitivity to pressure is thoroughly investigated using RRKM/ME theory [43].To simulate the reaction environment, argon (Ar) is chosen as the bath gas, and the collision frequency between the reactant and the bath gas is estimated utilizing the Lennard-Jones (L-J) parameters σ (in Å) and ε (in K).To determine these L-J parameters, we employ the method proposed by Wang and Frenklach [44].The σ and ε values employed in this study are outlined in Table S2 of Supplemental Material-SI.
In this study, we apply the one-dimensional asymmetric Eckart tunneling correction method [45,46] to account for tunneling effects.To rectify the substantial errors associated with approximating low-frequency vibrations of single bond torsions harmonically [47], we utilize a one-dimensional (1-D) hindered rotor model [48,49] to accurately represent the low-frequency vibrations of single bond torsions in the reaction.To determine the torsional potential energy of the reactants, we employed the B3LYP/CBSB7 theoretical level and conducted a dihedral angle scan at 10-degree intervals.To ensure the integrity of the reaction center, the atoms involved are initially fixed during the scanning process for the transition state.The potential energy profiles for internal rotations about C-C, C-OO, and CO-OH single bonds of the reactant, the transition state, and the product in a representative reaction are provided in the Figure S2 of the Supplemental Material-SI.In addition, in order to facilitate combustion modeling using widely-used software such as Chemkin-PRO (15092) [50], the thermodynamic parameters, namely enthalpy, entropy, and heat capacity, were fitted using the NASA format's 14 parameters.The outcomes of this fitting process can be found in Supplemental Material-SII.
In conclusion, the rate constants within the temperature range of 500 K to 1500 K are determined by fitting them using the three-parameter form of the Arrhenius Equation:

Conclusions
In this research study, systematic calculations were performed to determine the energy barriers, high-pressure limit rate constants, and pressure-dependent rate constants of β-scission reactions of •P(OOH) 2 compounds in normal-alkyl cyclohexanes.The calculations were carried out over a temperature range of 500-1500 K and a pressure range of 0.01-100 atm.The following conclusions can be drawn from our study: (1) The energy barriers and rate constants for the β-scission of •P(OOH) 2 reactions in our calculations did not exhibit significant variations.As a result, there is no need to further categorize these reactions into different subclasses.Instead, a comprehensive high-pressure limit rate rule and pressure-dependent rate rule can be established for this reaction class.(2) The maximum deviation in energy barriers and the uncertainty factor associated with the rate rule within this class were found to be within the bounds of chemical accuracy.This demonstrates the reliability and acceptability of our classification schemes.(3) The rate constants for the β-scission reactions of •P(OOH) 2 in normal-alkyl cyclohexanes exhibit significant discrepancies compared to the kinetic data provided in certain published mechanisms.This highlights the critical need to establish accurate rate rules for the β-scission reactions of •P(OOH) 2 of normal-alkyl cyclohexanes.(4) Pressure exerts a substantial influence on the rate constants of the β-scission reactions of •P(OOH) 2 of normal-alkyl cyclohexanes.Hence, our precisely-computed high-pressure limit rate rule and pressure-dependent rate rules are of immense importance, as they contribute to advancing the automatic generation of combustion mechanisms for normal-alkyl cyclohexanes.

Scheme 3 .
Scheme 3. The presence of hydrogen bonds in reactants (left) and transition state (right) for three different scenarios: (a) where both -OOH groups are positioned on the ring; (b) where one -OOH group is situated on the ring and one on the alkyl side chain; (c) where both -OOH groups are placed on the alkyl side chain.

Scheme 3 .
Scheme 3. The presence of hydrogen bonds in reactants (left) and transition state (right) for three different scenarios: (a) where both -OOH groups are positioned on the ring; (b) where one -OOH group is situated on the ring and one on the alkyl side chain; (c) where both -OOH groups are placed on the alkyl side chain.

a
The average energy barrier for reactions within this class.bThe maximum deviation for reactions within this class.

Figure 1 .
Figure 1.Comparison of the high-pressure limit rate constants: β-scission reaction rate rule of •P(OOH) 2 correspond to ECH1Q3Q4J = ECH3N1Q + HO2 in the mechanism from Zou et al.[32], PCHQ14-2 = HO2 + PCH1N4Q and PCHQ14-3 = HO2 + PCH3N1Q in the mechanism from Liu et al.[27], respectively.2.4.3.Comparison of the High-Pressure Limit Rate Rule with the Kinetic Data for the β-Scission Class in Non-Cyclic SystemsVillano et al.[33] systematically investigated the β-scission reactions of C2-C6 hydroperoxy alkyl radicals in non-cyclic systems using the CBS-QB3 level of theory.In their study, they assigned a single rate rule for the β-scission reaction within this reac-

Figure 2 .
Figure 2. Comparisons of the rate rule derived in this work (black lines) to previously reported rules:(a) the reported rules by Villano et al.[33];(b) the reported rules by Miyoshi[36].

Figure 3 .
Figure 3.Comparison of the pressure−dependent rate rule derived in this work with the values from the published mechanisms: (a) the combustion mechanism constructed by Zou et al. [32]; (b) the combustion mechanism constructed by Liu et al. [27].

Figure 3 .
Figure 3.Comparison of the pressure−dependent rate rule derived in this work with the values from the published mechanisms: (a) the combustion mechanism constructed by Zou et al. [32]; (b) the combustion mechanism constructed by Liu et al. [27].

Figure 4 .
Figure 4. Plot of the calculated average pressure-dependent rate constants of the class β−scission of •P(OOH)2:(a) the dependence on temperatures across the range of 500-1500 K; (b) the dependence on pressures across the range of 0.01-100 atm.

Figure 4 .
Figure 4. Plot of the calculated average pressure-dependent rate constants of the class β-scission of •P(OOH) 2 :(a) the dependence on temperatures across the range of 500-1500 K; (b) the dependence on pressures across the range of 0.01-100 atm.

Table 2 .
Cont.The average energy barrier for reactions within this class.b The maximum deviation fo within this class. a

Table 3 .
Comparison of the energy barriers between the β-scission reactions of •QOOH of •P(OOH)2 (kcal mol −1 The average energy barrier for reactions within this class.bThemaximum deviation fo within this class. ). β-Scission of •P(OOH)2 a

Table 3 .
Comparison of the energy barriers between the β-scission reactions of •QOOH and those of •P(OOH) 2 (kcal mol −1

2 This Work Energy Barrier β-Scission of •QOOH Previous Work Energy Barrier
).β-Scission of •P(OOH) a The average energy barrier for reactions within this class.b The maximum deviation for reactions within this class.

Table 3 .
Comparison of the energy barriers between the β-scission reactions of •QOOH and those of •P(OOH)2 (kcal mol −1

Table 3 .
Comparison of the energy barriers between the β-scission reactions of •QOOH and those of •P(OOH)2 (kcal mol −1

Scission of •P(OOH)2 This Work Energy Barrier β-Scission of •QOOH Previous Work Energy Barrier
).β-aThe average energy barrier for reactions within this class.b The maximum deviation for reactions within this class.

Table 3 .
Comparison of the energy barriers between the β-scission reactions of •QOOH and those of •P(OOH)2 (kcal mol −1
β-scission class (15.5 a , 2.8 b ) a The average energy barrier for reactions within this class.b The maximum deviation for reactions within this class.
β-scission class (15.5 a , 2.8 b ) a The average energy barrier for reactions within this class.b The maximum deviation for reactions within this class.
β-scission class (15.5 a , 2.8 b ) a The average energy barrier for reactions within this class.b The maximum deviation for reactions within this class.

Table 3 .
Comparison of the energy barriers between the β-scission reactions of •QOOH and those of •P(OOH)2 (kcal mol −1

Table 3 .
Comparison of the energy barriers between the β-scission reactions of •QOOH and those of •P(OOH)2 (kcal mol −1

Table 3 .
Comparison of the energy barriers between the β-scission reactions of •QOOH and those of •P(OOH)2 (kcal mol −1

Table 3 .
Comparison of the energy barriers between the β-scission reactions of •QOOH and those of •P(OOH)2 (kcal mol −1

Table 3 .
Comparison of the energy barriers between the β-scission reactions of •QOOH and those of •P(OOH)2 (kcal mol −1
a , 2.8 b ) a The average energy barrier for reactions within this class.b The maximum deviation for reactions within this class.
a , 2.8 b ) a The average energy barrier for reactions within this class.b The maximum deviation for reactions within this class.

Table 4 .
Analysis of rate constants, rate ratios, and uncertainty factor in β-scission class at high pressure.Comparison of the High-Pressure Limit Rate Constants with the Kinetic Data for the Corresponding Reactions in Published Mechanisms al., and Liu et al.

Table 5 .
Calculated rate constants, ratios, and uncertainty factor in β-scission class at different pressures.Comparisons of the rate rule derived in this work (black lines) to previously reported rules:

Table 5 .
Calculated rate constants, ratios, and uncertainty factor in β-scission class at different pressures.

-scission class rate rule 0.01 3.02 × 10 −10 5.70 −579.1 1.4 × 20.0 # 0.1 4.78 × 10 9 0.005 7720.47 3.9 × 10.9 # 1 7.93 × 10 3 1.98 5569.14 1.4 ×
k ave is the average rate constants at 800 K for class.#Theuncertainty factor values of the pressure-dependent rate rule at 800 K for class in our work (The ratios of the largest rate constant to the smallest rate constant at 800 K for class). e