Benchmark Calculations for Bond Dissociation Enthalpies of Unsaturated Methyl Esters and the Bond Dissociation Enthalpies of Methyl Linolenate

Theoretical Study on the P-C Bond Dissociation Enthalpy in 9,10-Dihydro-9-oxa-10-phosphaphenanthrene-10-oxide Flame Retardants

The 9,10-dihydro-9-oxa-10-phosphaphenanthrene-10-oxide (DOPO) flame retardants (DOPO-FRs) have attracted more and more attention in the flame-retardant industry due to their high efficiency, environmental protection, and good molecular design. During the flame-retardant process, the breakage of P-C bonds is very important to the flame-retardant effect. Through the comparison of different density functional theories (DFTs) on P-C BDEs, it was found that MN12-L has the highest calculation accuracy, and the root-mean-square error is the smallest with 1.85 kcal/mol. Therefore, MN12-L was selected to investigate P-C BDEs of different DOPO-FRs including thiophen-amine, benzo[d]thiazol-amine, triazol-amine, and aniline DOPO-FRs. By comparing the theoretical calculation of BDE with the experimental parameters of high limiting oxygen index (LOI) and vertical combustion test (UL-94 test), it was found that the P-C BDEs have a certain correlation with the flame-retardant effect. Finally, based on P-C BDEs, substituent effects, and effective flame-retardant fragments, a series of new DOPO-FRs were designed. The results showed that when only one DOPO fragment was contained, the effective fragments of flame retardants were ranked as furan > thiophene > triazole > imidazole. When bis-DOPO fragment was contained, the flame-retardant effect of diamino-triazole fragments was better than that of benzyldimethylamine fragments. In addition, when the substituents on the effective fragment have two EDGs, the flame-retardant effect is superimposed, which makes the flame-retardant performance more excellent.

Computing accurate bond dissociation energies of emerging per- and polyfluoroalkyl substances: Achieving chemical accuracy using connectivity-based hierarchy schemes

Understanding the bond dissociation energies (BDEs) of per- and polyfluoroalkyl substances (PFAS) helps in devising their efficient degradation pathways. However, there is only limited experimental data on the PFAS BDEs, and there are uncertainties associated with the BDEs computed using density functional theory. Although quantum chemical methods like the G4 composite method can provide highly accurate BDEs (< 1 kcal mol-1), they are limited to small system sizes. To address DFT's accuracy limitations and G4's system size constraints, we examined the connectivity-based hierarchy (CBH) scheme and found that it can provide BDEs that are reasonably close to the G4 accuracy while retaining the computational efficiency of DFT. To further improve the accuracy, we modified the CBH scheme and demonstrated that BDEs calculated using it have a mean-absolute deviation of 0.7 kcal mol-1 from G4 BDEs. To validate the reliability of this new scheme, we computed the ground state free energies of seven PFAS compounds and BDEs for 44 C-C and C-F bonds at the G4 level of theory. Our results suggest that the modified CBH scheme can accurately compute the BDEs of both small and large PFAS at near G4 level accuracy, offering promise for more effective PFAS degradation strategies.

Mechanism of formation of p-benzylenephenol peroxide radical (p-PhC(O2•)HPhOH)

CONTEXT

The reactions of radicals with O2 play the important role in the biological, medicinal, and industrial processes. The mechanism of this reaction is studied previously for the alkane, alkene, alkyne, phenols, and close-related radicals. According to these studies, the formation of intermediates in these reactions is predicted only for the aromatic radicals. Thus, the Van der Waals complexes of O2 with phenyl or benzyl radicals are predicted, as well as the π-π cluster for benzene. However, the possibility of the formation of such intermediate π-π clusters in the case of bisphenol radicals and the thermochemistry of its formation is not studied. Bisphenols are one of the main components of bio-oil, produced during pyrolysis of lignin-contained biomass. Synthetic bisphenols are used in polycarbonate plastics, epoxy resins, and thermal papers. Their mechanism of oxidation is important for the determination of fire safety of these materials, the possibility of using them as additives for fuels for the decreasing and the description of the ignition delays, as well as for the determination of its health risk assessment in medicine.

METHODS

The five DFT (M06-2X (i = 1), B3LYP (i = 2), wB97XD (i = 3), M08HX (i = 4), MN15 (i = 5)) approaches with 6-311 +  + G(d,p) basis set are used for the determination of standard enthalpies of atomization (ΔraHo(Yi)) of considered compounds (molecules, radicals, and transition states). These values of ΔraHo(Yi) are corrected using the empirical linear calibration dependencies, reported previously. The different calibration dependencies are used for the hydrocarbons (including the aromatics and simple oxygenated derivatives) and for the peroxides. The corrected values of ΔraHo(Yi, CORR) are used according to Hess's law for the determination of ΔfHo(Yi, CORR). The most consistent values of ΔfHo(Y, MEAN) are derived from the coordination of the values of ΔfHo(Yi, CORR) using the intersection of their values of standard deviations (3SDi). These values of ΔfHo(Y, MEAN), as well as the B3LYP values of So(Y), which are accounting the frequency correction and internal rotations, as well as their temperature dependencies, are used for the determination of thermochemistry of considered reactions and of the calculation, within transition state theory (TST), of the values of high-pressure limits of the rate constant. The values of Ho(Yi), So(Yi), and Go(Yi) are calculated using the Gaussian 16w program. The temperature dependencies of thermochemical properties and the values of rate constants are determined using the ChemRate program (v.1.5). The optimized structures are visualized using the Chemcraft.

A Bond-Energy/Bond-Order and Populations Relationship

We report an analytical bond energy from bond orders and populations (BEBOP) model that provides intramolecular bond energy decompositions for chemical insight into the thermochemistry of molecules. The implementation reported here employs a minimum basis set Mulliken population analysis on well-conditioned Hartree-Fock orbitals to decompose total electronic energies into physically interpretable contributions. The model's parametrization scheme is based on atom-specific parameters for hybridization and atom pair-specific parameters for short-range repulsion and extended Hückel-type bond energy term fitted to reproduce CBS-QB3 thermochemistry data. The current implementation is suitable for molecules involving H, Li, Be, B, C, N, O, and F atoms, and it can be used to analyze intramolecular bond energies of molecular structures at optimized stationary points found from other computational methods. This first-generation model brings the computational cost of a Hartree-Fock calculation using a large triple-ζ basis set, and its atomization energies are comparable to those from widely used hybrid Kohn-Sham density functional theory (DFT, as benchmarked to 109 species from the G2/97 test set and an additional 83 reference species). This model should be useful for the community by interpreting overall ab initio molecular energies in terms of physically insightful bond energy contributions, e.g., bond dissociation energies, resonance energies, molecular strain energies, and qualitative energetic contributions to the activation barrier in chemical reaction mechanisms. This work reports a critical benchmarking of this method as well as discussions of its strengths and weaknesses compared to hybrid DFT (i.e., B3LYP, M062X, PBE0, and APF methods), and other cost-effective approximate Hamiltonian semiempirical quantum methods (i.e., AM1, PM6, PM7, and DFTB3).

Master Equation Study of Hydrogen Abstraction from HCHO by OH Via a Chemically Activated Intermediate

The abstraction reaction of hydrogen from formaldehyde by OH radical plays an important role in formaldehyde oxidation. The reaction involves a bimolecular association to form a chemically activated hydrogen-bonded reaction...

Multiple free radical scavenging reactions of aurones

A series of naturally occurring 3',4'-dihydroxy aurones have been studied with regard to multiple free radical scavenging reactions in the gas and two liquid phases using density functional theory (DFT). All of the aurones prefer to perform (2 + n)-HAT mechanism to trap 2 + n free radicals, where n is the sum of the numbers of catechol and guaiacyl units in the gas and benzene phases. The second HAT reaction favours occurring in the same catechol moiety of the first HAT mechanism occurring OH group due to the formation of a stable quinone and the highly exothermic step of the final stable product formation. The catechol and guaiacyl moieties show increased potency for the second and fourth H⁺/e^‒ reactions. In the water phase, aurones can perform multiple H⁺/e^‒ reactions through n₁PL-ET-n₂HAT-(n+1-n₂)ET mechanism, where n₁ is the number of OH groups and n₂ is the number of guaiacyl moieties.

Multistructural Variational Reaction Kinetics of the Simplest Unsaturated Methyl Ester: H-Abstraction from Methyl Acrylate by H, OH, CH3, and HO2 Radicals

The H-abstraction reaction kinetics of methyl acrylate (MA) + H/OH/CH₃/HO₂ radicals have been investigated theoretically in the present work. For these reactions, the reaction energies and barrier heights are first computed using several density functionals and compared to the coupled cluster CCSD(T)-F12/jun-cc-pVTZ benchmark calculations. The M062X/maug-cc-pVTZ method shows the best performance with the smallest mean unsigned deviation (MUD) of 0.42 kcal mol^-1. Combined with the electronic structure calculations using the M062X/maug-cc-pVTZ method, the multistructural canonical variational transition-state theory (MS-CVT) with small-curvature tunneling (SCT) is employed to calculate the reaction rate constants at 500-2000 K. The variational effect is between 0.56 and 1.0, the multistructural torsional anharmonicity factor ranges from 0.004 to 4.57, and the tunneling coefficient is in the range of 0.5-4.70. Notably, given the existence of reactant complexes (RCs) between reactants and transition states for the reaction systems MA + OH/HO₂, we further compare the rate constants under the low-pressure limit (LPL) kinetic model, which treats the reaction as a single-step process and neglects RCs, and the pre-equilibrium model, which takes RCs into account in the reaction and treats the reaction as a two-step process. The rate constants calculated by these two models are similar within the combustion temperature range, and apparent differences occur at lower temperatures. In addition, we determine the branching ratios as a function of temperature and find that the methyl site (S3) abstractions by OH and H radicals are dominant in the low- and high-temperature ranges, respectively. Moreover, we update the kinetic model with the calculated H-abstraction rate constants to simulate the ignition delay times of MA. The simulations of the updated model are in good agreement with experimental results. The accurate reaction kinetics determined in this work are useful for the understanding and prediction of consumption branching fractions and ignition properties of the unsaturated methyl esters.

Hydrogenation of CO2 Promoted by Silicon-Activated H2S: Origin and Implications

Unlike the usual method of CO_x (x = 1, 2) hydrogenation using H₂ directly, H₂S and HSiSH (silicon-activated H₂S) were selected as alternative hydrogen sources in this study for the CO_x hydrogenation reactions. Our results suggest that it is kinetically infeasible for hydrogen in the form of H₂S to transfer to CO_x at low temperatures. However, when HSiSH is employed instead, the title reaction can be achieved. For this approach, the activation of CO₂ is initiated by its interaction with the HSiSH molecule, a reactive species with both a hydridic H^δ− and protonic H^δ+. These active hydrogens are responsible for the successive C-end and O-end activations of CO₂ and hence the final product (HCOOH). This finding represents a good example of an indirect hydrogen source used in CO₂ hydrogenation through reactivity tuned by silicon incorporation, and thus the underlying mechanism will be valuable for the design of similar reactions.

BonDNet: a graph neural network for the prediction of bond dissociation energies for charged molecules

A broad collection of technologies, including e.g. drug metabolism, biofuel combustion, photochemical decontamination of water, and interfacial passivation in energy production/storage systems rely on chemical processes that involve bond-breaking molecular reactions. In this context, a fundamental thermodynamic property of interest is the bond dissociation energy (BDE) which measures the strength of a chemical bond. Fast and accurate prediction of BDEs for arbitrary molecules would lay the groundwork for data-driven projections of complex reaction cascades and hence a deeper understanding of these critical chemical processes and, ultimately, how to reverse design them. In this paper, we propose a chemically inspired graph neural network machine learning model, BonDNet, for the rapid and accurate prediction of BDEs. BonDNet maps the difference between the molecular representations of the reactants and products to the reaction BDE. Because of the use of this difference representation and the introduction of global features, including molecular charge, it is the first machine learning model capable of predicting both homolytic and heterolytic BDEs for molecules of any charge. To test the model, we have constructed a dataset of both homolytic and heterolytic BDEs for neutral and charged (-1 and +1) molecules. BonDNet achieves a mean absolute error (MAE) of 0.022 eV for unseen test data, significantly below chemical accuracy (0.043 eV). Besides the ability to handle complex bond dissociation reactions that no previous model could consider, BonDNet distinguishes itself even in only predicting homolytic BDEs for neutral molecules; it achieves an MAE of 0.020 eV on the PubChem BDE dataset, a 20% improvement over the previous best performing model. We gain additional insight into the model's predictions by analyzing the patterns in the features representing the molecules and the bond dissociation reactions, which are qualitatively consistent with chemical rules and intuition. BonDNet is just one application of our general approach to representing and learning chemical reactivity, and it could be easily extended to the prediction of other reaction properties in the future.

Role of C‒H bond in the antioxidant activities of rooperol and its derivatives: A DFT study

Rooperol and its derivatives, derived from the Hypoxis rooperi plant, are polyphenolic and norlignan compounds with excellent antioxidant activities. The reaction enthalpies for the free-radical scavenging by rooperol and its six derivatives were studied using density functional theory. We found that the C-H groups played a significant role in the antioxidant activities in non-polar phases. In the gas and benzene phases, rooperol and its derivatives preferentially underwent the free-radical scavenging process via the 3‒CH group by following the hydrogen atom transfer (HAT) mechanism. In polar phases, the sequential proton loss electron transfer (SPLET) was the most preferred mechanism, and the phenolic O‒H groups played a significant role. Additionally, we found that when the hydrogen atom in the OH group was replaced by a glucose moiety, the antioxidant activity of the adjacent OH group was reduced. ROP, DHROP-I, DHROP-II, ROP-4″-G and ROP-4'G have catechol moiety, they may proceed double step-wise mechanisms to trap free radicals. In the gas and benzene phases, the preferable mechanism is dHAT. In water phase, it is SPLHAT.

Kinetic modeling of methyl pentanoate pyrolysis based on ab initio calculations

Recently, methyl pentanoate (MP) was proposed as a viable biodiesel surrogate to petroleum-based fuels. To better understand the pyrolysis chemistry of MP, the unimolecular decomposition kinetics of MP is theoretically investigated on the basis of ab initio calculations; ten primary channels, including four intramolecular H-shifts and six C-C and C-O bond fissions, are identified. The geometries are optimized at the M06-2X/cc-pVTZ level of theory, and accurate barrier heights are determined using the DLPNO-CCSD(T)/CBS(T-Q) method, which shows a good performance against the CCSD(T)/CBS(T-Q) method with an uncertainty of 0.5 kcal mol-1 for small methyl esters. The atomization enthalpy method is adopted to obtain the thermodynamics of involved species. The Rice-Ramsperger-Kassel-Marcus/master equation theory coupled with one-dimensional hindered rotor approximation is employed to calculate the phenomenological rate constants at 500-2000 K and 0.01-100 atm. The branching ratio analysis indicates that two reactions, MP ↔ CH3OC([double bond, length as m-dash]O)CH3 + CH2CHCH3 and MP ↔ CH3OC([double bond, length as m-dash]O)CH2 + CH2CH2CH3, are the dominant channels at low and high temperatures, respectively. The model from Diévart et al. [Proc. Combust. Inst., 2013, 34(1), 821-829] is updated with our calculations, and the modified model can yield a better prediction in reproducing the ignition delay times of MP at high temperatures. This work provides a comprehensive investigation of MP unimolecular decomposition, and can serve as a prototype for understanding the pyrolysis of larger alkyl esters.

Prediction of organic homolytic bond dissociation enthalpies at near chemical accuracy with sub-second computational cost

Bond dissociation enthalpies (BDEs) of organic molecules play a fundamental role in determining chemical reactivity and selectivity. However, BDE computations at sufficiently high levels of quantum mechanical theory require substantial computing resources. In this paper, we develop a machine learning model capable of accurately predicting BDEs for organic molecules in a fraction of a second. We perform automated density functional theory (DFT) calculations at the M06-2X/def2-TZVP level of theory for 42,577 small organic molecules, resulting in 290,664 BDEs. A graph neural network trained on a subset of these results achieves a mean absolute error of 0.58 kcal mol^-1 (vs DFT) for BDEs of unseen molecules. We further demonstrate the model on two applications: first, we rapidly and accurately predict major sites of hydrogen abstraction in the metabolism of drug-like molecules, and second, we determine the dominant molecular fragmentation pathways during soot formation.

Reaction kinetics of hydrogen addition reactions to methyl butenoate

We study the chemical kinetics of hydrogen addition reactions of unsaturated methyl esters, methyl 2-butenoate and methyl 3-butenoate, and compare the rate constants with those of hydrogen abstraction reactions by H atom.

On the Prediction of Standard Enthalpy of Formation of C2-C4 Oxygenated Species

In this study, to determine an efficient and accurate method for predicting standard enthalpy of formation (Δ_fH^o) of oxygenated species, we calculated Δ_fH^o for several typical C2-C4 oxygenated species using atomization and isodesmic reactions in combination with various quantum chemical methods, including six density functional theory methods, three compound methods, and CCSD(T)/CBS. Compared with experimental values, at the same quantum chemical level, Δ_fH^o values predicted by using isodesmic reactions are more accurate than those using atomization reactions. Comparing various quantum chemical methods when isodesmic reactions are used, the performance of G4 is the best with a mean unsigned deviation (MUE) of 0.3 kcal/mol and a standard deviation (SD) of 0.3 kcal/mol, while M06-2X can predict Δ_fH^o efficiently and accurately with an MUE of 0.6 kcal/mol and SD of 0.5 kcal/mol. Using the best methods we have found, we calculated the enthalpies of formation and other thermodynamic properties for dimethyl carbonate (DMC) and its associated species and then applied them in a DMC combustion model for predicting ignition delay times. Better agreement with the experiments is achieved when the newly computed thermodynamic properties are adopted.

A Benchmark of Density Functional Approximations For Thermochemistry and Kinetics of Hydride Reductions of Cyclohexanones

The performance of density functionals and wavefunction methods for describing the thermodynamics and kinetics of hydride reductions of 2-substituted cyclohexanones has been evaluated for the first time. A variety of exchange correlation functionals ranging from generalized gradient approximations to double hybrids have been tested and their performance to describe the facial selectivity of hydride reductions of cyclohexanones has been carefully assessed relative to the CCSD(T) method. Among the tested methods, an approach in which single-point energy calculations using the double hybrid B2PLYP-D3 functional on ωB97X-D optimized geometries provides the most accurate transition state energies for these kinetically-controlled reactions. Moreover, the role of torsional strain, temperature, solvation, noncovalent interactions on the stereoselectivity of these reductions was elucidated. Our results indicate a prominent role of the substituent on the cis/trans ratios driven by the delicate interplay between torsional strain and dispersion interactions.

Effect of hindered internal rotation treatments on predicting the thermodynamic properties of alkanes

We compare different vibrational analysis methods and examine the effect of hindered internal rotation treatments on predicting thermodynamic properties.

Computational Study on N-N Homolytic Bond Dissociation Enthalpies of Hydrazine Derivatives

The hydrazine derivatives have been regarded as the important building blocks in organic chemistry for the synthesis of organic N-containing compounds. It is important to understand the structure-activity relationship of the thermodynamics of N-N bonds, in particular, their strength as measured by using the homolytic bond dissociation enthalpies (BDEs). We calculated the N-N BDEs of 13 organonitrogen compounds by eight composite high-level ab initio methods including G3, G3B3, G4, G4MP2, CBS-QB3, ROCBS-QB3, CBS-Q, and CBS-APNO. Then 25 density functional theory (DFT) methods were selected for calculating the N-N BDEs of 58 organonitrogen compounds. The M05-2X method can provide the most accurate results with the smallest root-mean-square error (RMSE) of 8.9 kJ/mol. Subsequently, the N-N BDE predictions of different hydrazine derivatives including cycloalkylhydrazines, N-heterocyclic hydrazines, arylhydrazines, and hydrazides as well as the substituent effects were investigated in detail by using the M05-2X method. In addition, the analysis including the natural bond orbital (NBO) as well as the energies of frontier orbitals were performed in order to further understand the essence of the N-N BDE change patterns.

Kinetics and branching fractions of the hydrogen abstraction reaction from methyl butenoates by H atoms

We studied hydrogen abstraction reactions at various sites of unsaturated methyl esters by H atoms, including variational effects, tunneling and multi-structural torsional anharmonicity.

Reaction kinetics of hydrogen abstraction from polycyclic aromatic hydrocarbons by H atoms

We studied how the PAH structure, site, and size affect the rate constants of hydrogen abstraction reactions of PAH systematically.

MN15: A Kohn-Sham global-hybrid exchange-correlation density functional with broad accuracy for multi-reference and single-reference systems and noncovalent interactions

Kohn-Sham density functionals are widely used; however, no currently available exchange-correlation functional can predict all chemical properties with chemical accuracy. Here we report a new functional, called MN15, that has broader accuracy than any previously available one. The properties considered in the parameterization include bond energies, atomization energies, ionization potentials, electron affinities, proton affinities, reaction barrier heights, noncovalent interactions, hydrocarbon thermochemistry, isomerization energies, electronic excitation energies, absolute atomic energies, and molecular structures. When compared with 82 other density functionals that have been defined in the literature, MN15 gives the second smallest mean unsigned error (MUE) for 54 data on inherently multiconfigurational systems, the smallest MUE for 313 single-reference chemical data, and the smallest MUE on 87 noncovalent data, with MUEs for these three categories of 4.75, 1.85, and 0.25 kcal mol-1, respectively, as compared to the average MUEs of the other 82 functionals of 14.0, 4.63, and 1.98 kcal mol-1. The MUE for 17 absolute atomic energies is 7.4 kcal mol-1 as compared to an average MUE of the other 82 functionals of 34.6 kcal mol-1. We further tested MN15 for 10 transition-metal coordination energies, the entire S66x8 database of noncovalent interactions, 21 transition-metal reaction barrier heights, 69 electronic excitation energies of organic molecules, 31 semiconductor band gaps, seven transition-metal dimer bond lengths, and 193 bond lengths of 47 organic molecules. The MN15 functional not only performs very well for our training set, which has 481 pieces of data, but also performs very well for our test set, which has 823 data that are not in our training set. The test set includes both ground-state properties and molecular excitation energies. For the latter MN15 achieves simultaneous accuracy for both valence and Rydberg electronic excitations when used with linear-response time-dependent density functional theory, with an MUE of less than 0.3 eV for both types of excitations.