Prediction of Bond Dissociation Energy for Organic Molecules Based on a Machine-Learning Approach

Bond dissociation energy(BDE),which refers to the enthalpy change for the homolysis of a specific covalent bond,is one of the basic thermodynamic properties of molecules.It is very important for understanding chemical reactivities,chemical properties and chemical transformations.Here,a machine learning-based comprehensive BDE prediction model was established based on the iBonD experimental BDE dataset and the calculated BDE dataset by St.John et al.Differential Structural and PhysicOChemical(D-SPOC)descriptors that reflected changes in molecules'structural and physicochemical features in the process of bond homolysis were designed as input features.

Diffusion Monte Carlo Study of Bond Dissociation Energies for BH2, B（OH）2, BCl2, and BCI

On basis of bond dissociation energies （BDEs） for BH2, B（OH）2, BCl2, and BCl, the diffusion Monte Carlo （DMC） method is applied to explore the BDEs of HB-H, HOB-OH, ClB-Cl, and B-Cl. The effect of the choice of orbitals, as well as the backflow transformation, is studied. The Slater-Jastrow DMC algorithm gives BDEs of 359.1±0.12 kJ/mol for HB-H, 410.5±0.50 kJ/mol for HOB-OH, 357.8±1.46 kJ/mol for ClB-Cl, and 504.5±0.96 kJ/mol for B-Cl using B3PW91 orbitals and similar BDEs when B3LYP orbitals are used. DMC with backflow corrections （BF-DMC） gives a HB-H BDE of 369.9±0.12 kJ/mol which is close to one of the available experimental value （375.8 kJ/mol）. In the case of HOB-OH BDE, the BF-DMC calculation is 446.04-1.84 k J/mol that is closer to the experimental BDE. The BF-DMC BDE for ClB-Cl is 343.2±2.34 kJ/mol and the BF-DMC B-Cl BDE is 523.3±0.33 kJ/mol, which are close to the experimental BDEs, 341.9 and 530.0 kJ/mol, respectively.

The evaluation of bond dissociation energies for NO2 scission in nitro compounds using density functional and complete basis set methods

By using the density functional theory （B3LYP） and four highly accurate complete basis set （CBS-Q, CBS-QB3, CBS-Lq, and CBS-4M）ab initio methods, the X（C, N, O）-NO2 bond dissociation energies （BDEs） for CH3NO2, C2H3NO2, C2H5NO2, HONO2, CH3ONO2, C2H5ONO2, NH2NO2 （CH3）2NNO2 are computed. By comparing the computed BDEs and experimental results, it is found that the B3LYP method is unable to predict satisfactorily the results of bond dissociation energy （BDE）; however, all four CBS models are generally able to give reliable predication of the X（C, N, O）-NO2 BDEs for these nitro compounds. Moreover, the CBS-4M calculation is the least computationally demanding among the four CBS methods considered, Therefore, we recommend CBS-4M method as a reliable method of computing the BDEs for this nitro compound system.

Quantum chemical calculations of bond dissociation energies for COOH scission and electronic structure in some acids

Quantum chemical calculations are performed to investigate the equilibrium C-COOH bond distances and the bond dissociation energies（BDEs） for 15 acids.These compounds are studied by utilizing the hybrid density functional theory（DFT）（B3LYP,B3PW91,B3P86,PBE1PBE） and the complete basis set（CBS-Q） method in conjunction with the 6311G^＊＊ basis as DFT methods have been found to have low basis sets sensitivity for small and medium molecules in our previous work.Comparisons between the computational results and the experimental values reveal that CBS-Q method,which can produce reasonable BDEs for some systems in our previous work,seems unable to predict accurate BDEs here.However,the B3P86 calculated results accord very well with the experimental values,within an average absolute error of 2.3 kcal/mol.Thus,B3P86 method is suitable for computing the reliable BDEs of C-COOH bond for carboxylic acid compounds.In addition,the energy gaps between the highest occupied molecular orbital（HOMO） and the lowest unoccupied molecular orbital（LUMO） of studied compounds are estimated,based on which the relative thermal stabilities of the studied acids are also discussed.

Density Functional Calculations of C–NO_2 Bond Dissociation Energies for Nitroalkanes Molecules

Bond dissociation energies for the removal of nitrogen dioxide group in some nitroalkane energetic materials have been calculated by using the three hybrid density functional theory （DFT） methods B3LYP, B3PW91 and B3P86 with 6-31g^＊＊ and 6-311g^＊＊ basis sets. The computed BDEs have been compared with the available experimental results. It is found that the B3P86 method with 6-31g^＊＊ and 6-311g^＊＊ basis sets can obtain satisfactory bond dissociation energies （BDEs）, which are in extraordinary agreement with the experimental data. Considering the smaller mean absolute deviation and maximum difference, the reliable B3P86/6-311g^＊＊ method was recommended to compute the BDEs for the removal of nitrogen dioxide group in the nitroalkane energetic materials. Using the method, the BDEs of 8 other nitroalkane energetic materials have been calculated and the maximum difference from experimental value is 1.76 kcal.mo1^-1 （for the BDE of tC4Hg-NOz）, which further proves the reliability of B3P86/6-311g^＊＊ method. In addition, it is noted that the BDEs of C-NO2 bond change slightly for main chain nitroalkane compounds with the maximum difference of only 3.43 kcal mo1^-1.

Theoretical Study of the N-NO_2 Bond Dissociation Energies for Energetic Materials with Density Functional Theory

The N-NO2 bond dissociation energies （BDEs） for 7 energetic materials were computed by means of accurate density functional theory （B3LYP, B3PW91 and B3P86） with 6-31G＊＊ and 6-311G＊＊ basis sets. By comparing the computed energies and experimental results, we find that the B3P86/6-311G＊＊ method can give good results of BDE, which has the mean absolute deviation of 1.30kcal/mol. In addition, substituent effects were also taken into account. It is noted that the Hammett constants of substituent groups are related to the BDEs of the N-NO2 bond and the bond dissociation energies of the energetic materials studied decrease when increasing the number of NO2 group.

Is HOMO Energy Level a Good Parameter to Characterize Antioxidant Activity

Semiempirical quantum chemical method AM1 was employed to calculate the highest occupied molecular orbital (HOMO) energy levels (E-HOMO) for various types of antioxidants. It was verified that the correlation between logarithm of free radical scavenging rate constants (1gks) and E-HOMO substantially arises from the correlation between E-HOMO and O-H bond dissociation energies (BDE) of antioxidants. Furthermore, E-HOMO were poorly correlated with the logarithm of relative free radical scavenging rate constants (1gk(3)/k(1)) for various types of antioxidants that possess complex structures (r = 0.5602). So in a broad sense, E-HOMO was not an appropriate parameter to characterize the free radical scavenging activity of antioxidants.

Research on C—C Bond Length Distribution in Hydrocarbon Molecules

The C--C bond dissociation energy （BDE） is a very important data in research of hydrocarbon cracking reactions, because it reflects the difficulty level of chemical reactions. But it is very difficult to obtain the C--C bond dissociation energy （BDE） by experiments, so using quantum chemistry calculation such as density functional theory （DFT） to study the C--C bond dissociation energy is a very useful means. The impact of acceptor substituents and donor substituents on the C--C bond length distribution was studied.

Density Functional Theory Study of Marine Polybrominated Diphenyl Ethers in Anaerobic Degradation

Polybrominated diphenyl ethers(PBDEs)are a kind of serious pollutants in the ocean.Biodegradation is considered as an economical and safe way for PBDEs removal and reductive debromination dominates the initial pathway of anaerobic degradation.On the basis of experimental study,Octa-BDE 197,Hepta-BDE 183,Hexa-BDE 153,Penta-BDE 99 and Tetra-BDE 47 were selected as the initial degradation objects,and their debromination degradation were studied using density functional theory.The structures were optimized by Gaussian 09 program.Furthermore,the molecular orbitals and charge distribution were analyzed.All C-Br bond dissociation energies at different positions including ortho,meta and para bromine atoms were calculated and the sequence of debromination was obtained.There is a close relationship between molecular structure,charge,molecular orbital and C-Br bond.All PBDEs exhibited similar debromination pathways with preferential removal of meta and para bromines.

Theoretical Study of Remote Substituent Effects on X-H （X=CH2, NH, O） Bond Dissociation Energies of Azulene

In the study, the X-H （X=CH2, NH, O） bond dissociation energies （BDE） of para-substituted azulene （Y-C10H8X-H） were predicted theoretically for the first time using Density Functronal Theory （DFT） methods at UB3LYP/6-311 ＋ ＋g（2df,2p）//UB3LYP/6-31 ＋g（d） level. It was found that the substituents exerted similar effects on the X-H BDE of azulene as those on benzene, except for 6-substituted 2-methylazulene. Owing to the substituent-dipole interaction, the reaction constants （ρ^＋） of 2- and 6-Y-CIoHsX-H （X=NH and O only） varied violently. The origin of the substituent effects on the X-H BDE of azulene was found, by both GE/RE and SIE theory, to be directly associated with variation of the radical effects, although the ground effects also played a modest role in determining the net. substituent effects.

Relativistic density functional investigation of UX_6(X=F,Cl,Br and I)

The molecular structures and the vibrational frequencies of uranium hexahalides UX6 （X=F, Cl, Br and I） molecules are investigated by using local density approximation （LDA） and generalised gradient approximation （GGA） functions （BP, BLYP and RPBE） in combination with two different relativistic methods （scalar and scalar＋spin-orbit relativistic effects）. The calculated results show that the differences are trivial between scalar and scalar＋spin-orbit relativistic methods. The vibrational frequencies are also compared with existing experimental values, and overall, the RPBE approach gives the smallest error. The bond dissociation energies （BDEs） of UX6 are computed by using the RPBE function, thereby obtaining exact vibrational frequencies. In addition, the calculated magnitudes of the spin orbit effect on the BDE of UX6 （X=F, Cl, Br, and I） are found to be approximately -0.3198, 0.3218, -0.3609 and -0.4415 eV, respectively.

DFT Study on the Antioxidant Activity of a Modeled p-Terphenyl Derivative

Relationships between the structure characteristics of natural p-terphenyl com- pounds isolated from three edible mushrooms （Thelephora ganbajun, Thelephora aeronautical, and Boletopsis grisea） indigenous to China and their mechanism of antioxidant activity were studied. Geometry structures of terphenyl molecule and four corresponding radicals, bond dissociation energy （BDE）, frontier orbitals （HOMO and LUMO） and single electron density were calculated using DFT methods （B3LYP/6-311G＊＊）. The computational results which are consistent with the experimental data well show that terphenyl molecule scavenges DPPH radical by hydrogen abstract mechanism and the high antioxidant activity depends on the substitution position of hydroxyls. Two active 7-, 8-hydroxyls facilitate the hydrogen abstraction due to the intramolecular hydrogen bond and the resonance effect makes 4-hydroxyl radical more stable.

Stability,detonation properties and pyrolysis mechanisms of polynitrotriprismanes C_6H_(6-n)(NO_2)_n (n=1-6)

To further test whether polynitriprismanes are capable of being potential high energy density materials （HEDMs）, extensive theoretical calculations were carried out to investigate on a series of polynitrotriprismanes （PNNPs）： C6H6-.（NO2）. （n=1-6） Heats of formation （HOFs）, strain energies （SE）, and disproportionation energy （DE） were obtained using B3LYP/6-311＋G（2df, 2p）//B3LYP/6-31G＊ method by designing different isodesmic reactions, respectively. Detonation properties of PNNPs were obtained by the well-known KAMLET-JACOBS equations, using the predicted densities （p） obtained by Monte Carlo method and HOFs. It is found that they increase as the number of nitro groups n varies from 1 to 6, and PNNPs with n〉4 have excellent detonation properties The relative stability and the pyrolysis mechanism of PNNPs were evaluated by the calculated bond dissociation energy （BDE）. The comparison of BDE suggests that rupturing the C--C bond is the trigger for thermolysis of PNNPs. The computed BDE for cleavage of C--C bond （88.5 kJ/mol） further demonstrates that only the hexa-nitrotriprismane can be considered to be the target of HEDMs.

α-Substituent Effects on Si--H, P- H and S--H Bond Dissociation Energies

CBS-Q and G3 methods were used to generate a large number of reliable Si--H, P---H and S--H bond dissociation energies （BDEs） for the first time. It was found that the Si--H BDE displayed dramatically different substituent effects compared with the C--H BDE. On the other hand, the P---H and S--H BDE exhibited patterns of substituent effects similar to those of the N--H and O--H BDE. Further analysis indicated that increasing the positive charge on Si of XSiH3 would strengthen the Si--H bond whereas increasing the positive charge on P and S of XPH2 and XSH would weaken the P---H and S--H bonds. Meanwhile, increasing the positive charge on Si of XSiH2^＋ stabilized the silyl radical whereas increasing the positive charge on P and S in XPH＂ and XS＊ destabilized P- and S-centered radicals. These behaviors could be reasonalized by the fact that Si is less electronegative than H while P and S are not. Finally, it was demonstrated that the spin-delocalization effect was valid for the Si-, P- and S-centered radicals.

Theoretical studies on a series of nitroaliphatic energetic compounds

Density functional theory calculations at the B3LYP/6-311G＊＊ level are performed to study the geometric and elec- tronic structures of a series of nitroaliphatic compounds. The heats of formation （HOF） are predicted through the designed isodesmic reactions. Thermal stabilities are evaluated via the homolytic bond dissociation energies （BDEs）. Further, the correlation is developed between impact sensitivity h50% and the ratio （BDE/E） of the weakest BDE to the total energy E containing zero point energy correction. In addition, the relative stability of the title compounds is evaluated based on the analysis of calculated Mulliken population and the energy gaps between the frontier orbitals. The calculated BDEs, HOFs, and energy gaps consistently indicate that compound 1,1,1,6,6,6-hexanitro-3-hexyne is the most unstable and the compound 3,3,4,4,-tetranitro-hexane is the most stable. These results provide basic information for the molecular design of novel high energetic density materials.

Computational studies on the structure and detonation properties of nitro-substituted triazole-furazan derivatives

Some nitro-substituted triazole-furazan derivatives are considered as potential candidates for high energy density compounds through quantum chemical treatment. Their geometric and electronic structures,band gap,thermodynamic properties and detonation properties were studied using the density functional theory at the B3 LYP /6- 311 ＋ G＊＊level. The calculated energy of explosion,density,and detonation properties of model compounds were comparable to 1,3,5-trinitro-1,3,5-triazinane（ RDX） and 1,3,5,7-tetranitro-1,3,5,7-tetrazocane（ HMX）. The heats of formation and bond dissociation energy were also analysed to understand the nature of thermal stabilities and the trigger bond in the pyrolysis process.

Theoretical Investigations on the Structure,Density,Thermodynamic and Performance Properties of Bis(2,2-dinitropropyl) formal

Density functional method was used to investigate the IR spectrum, heat of forma- tion and thermal stability of a new energetic material bis（2,2-dinitropropyl） formal （BDNPF）. The detonation velocity and pressure were evaluated by using the Kamlet-Jacobs equations based on the theoretical density and heat of formation. The bond dissociation energies for the weakest bonds were analyzed to investigate the thermal stability of the title compound. The results show that the C（I ）-N（I ） bond is predicted to be the trigger bond during pyrolysis. The crystal structure obtained by molecular mechanics belongs to the P21 space group, with the lattice parameters to be Z = 2, a = 11.5254, b = 6.2168, c = 9.5000 A andp= 1.66 g/cm3.

Heterolytic and homolytic C-D bond dissociation energies of NADH models in acetonitrile and primary isotopic effects on hydride versus hydrogen atom transfer reactions

Heterolytic and homolytic C D bond dissociation energies of three NADH models: BNAH-4,4-d 2 , HEH-4,4-d 2 and AcrD 2 in acetonitrile were first estimated by using an efficient method. The results showed that the heterolytic C D bond dissociation energies are 65.2, 70.2, and 81.9 kcal/mol and the homolytic C D bond dissociation energies are 72.66, 70.69, and 74.95 kcal/mol for BNAH-4,4-d 2 , HEH-4,4-d 2 , and AcrD 2 , respectively. According to the bond dissociation energy differences of isotope isomers, an interesting conclusion can be made that the primary kinetic isotope effects are dependent not only on the zero-point energy difference of the isotope isomers, but also on the types of C D bond dissociations, and the C D bond homolytic dissociations should have much larger primary kinetic isotope effects (26.9 28.8) than the corresponding C D bond heterolytic dissociations (3.9-5.4).