Retro-[3 + 2]-cycloaddition reactions in the decomposition of five-membered nitrogen-containing heterocycles

Nitrogen-containing heterocycles form the basis for a new generation of high-energy density materials, and they serve as model compounds for nitrogen-containing fuels, such as coal and biomass, and they form the backbone of ionic liquids. A novel retro-[3 + 2]-cycloaddition to a three-membered diene and a two-membered dienophile, analogous to a retro-Diels-Alder reaction, may constitute an important initial reaction step in the thermal decomposition of these heterocyclic compounds. We investigate the kinetics and thermodynamics of these reactions for the heterocycles pyrrole, pyrazole, imidazole, 1,2,3-triazole, 1,2,4-triazole, 1,2,5-triazole, 1,3,4,-triazole, 1H-tetrazole, and 2H-tetrazole, using theoretical computational chemistry. The retro-cycloadditions are shown to form one of the three-membered products: hydrazoic acid (NH=N=N), nitrilimine (NH=N=CH), carbodiimide (NH=C=NH), or ketenimine (NH=C=CH2) plus one of the two-membered products acetylene, hydrogen cyanide, or N2. Accurate enthalpies of formation are calculated for the reaction products using the high-level W1 computational protocol, providing the previously undetermined enthalpy values of 70.09, 88.75, 35.03, and 44.28 kcal mol(-1) for hydrazoic acid, nitrilimine, carbodiimide, and ketenimine, respectively. We apply a variable-order form of the Marcus equation to the dissociation reactions in correlating the enthalpy of reaction with the activation enthalpy. Typical molecular elimination reactions from the heterocycles proceed with an intrinsic activation enthalpy of 36.8 kcal mol(-1) and intrinsic activation free energy of 42.1 kcal mol(-1). However, dissociation reactions resulting in the formation of either NH=C=NH or NH=C=CH2 demonstrate intrinsic barriers ca. 30 kcal mol(-1) higher, as a result of a concerted intramolecular hydrogen shift. Rate constants calculated between 300 and 3000 K indicate that the proposed dissociation reactions should be important in the decomposition of tetrazole and 1,2,3-triazole. This is confirmed by comparison with available experimental data. Decomposition of 1,2,4-triazole to HCN + nitrilimine may also be important at high temperatures. From extrapolation of our Marcus equation relationship, we predict pentazole to decompose to N2 + NHNN with an activation enthalpy of 19.5 kcal mol-1 and a half-life of only 14 s at 298 K.     

Enthalpies of formation, bond dissociation energies, and molecular structures of the n-aldehydes (acetaldehyde, propanal, butanal, pentanal, hexanal, and heptanal) and their radicals

Aldehydes are important intermediates and products in a variety of combustion and gas-phase oxidation processes, such as in low-temperature combustion, in the atmosphere, and in interstellar media. Despite their importance, the enthalpies of formation and bond dissociation energies (BDEs) for the aldehydes are not accurately known. We have determined enthalpies of formation for acetaldehyde, propanal, and butanal from thermodynamic cycles, using experimentally measured reaction and formation enthalpies. All enthalpy values used for reference molecules and reactions were first verified to be accurate to within around 1 kcal mol-1 using high-level ab initio calculations. Enthalpies of formation were found to be -39.72 +/- 0.16 kcal mol-1 for acetaldehyde, -45.18 +/- 1.1 kcal mol-1 for propanal, and -49.27 +/- 0.16 kcal mol-1 for butanal. Enthalpies of formation for these three aldehydes, as well as for pentanal, hexanal, and heptanal, were calculated using the G3, G3B3, and CBS-APNO theoretical methods, in conjunction with bond-isodesmic work reactions. On the basis of the results of our thermodynamic cycles, theoretical calculations using isodesmic work reactions, and existing experimental measurements, we suggest that the best available formation enthalpies for the aldehydes acetaldehyde, propanal, butanal, pentanal, hexanal, and heptanal are -39.72, -45.18, -50.0, -54.61, -59.37, and -64.2 kcal mol-1, respectively. Our calculations also identify that the literature enthalpy of formation of crotonaldehyde is in error by as much as 1 kcal mol-1, and we suggest a value of -25.1 kcal mol-1, which we calculate using isodesmic work reactions. Bond energies for each of the bonds in the aldehydes up to pentanal were calculated at the CBS-APNO level. Analysis of the BDEs reveals the R-CH(2)CH=O to be the weakest bond in all aldehydes larger than acetaldehyde, due to formation of the resonantly stabilized vinoxy radical (vinyloxy radical/formyl methyl radical). It is proposed that the vinoxy radical as well as the more commonly considered formyl and acetyl radicals are important products of aldehyde combustion and oxidation, and the reaction pathways of the vinoxy, formyl, and acetyl radicals are discussed. Group additivity values for the carbon-oxygen-hydrogen groups common to the aldehydes are also determined. Internal rotor profiles and electrostatic potential surfaces are used to study the dipole induced dipole-dipole interaction in the synperiplanar conformation of propanal. It is proposed that the loss of this dipole-dipole interaction in RC(.-)HCH(2)CH=O radicals causes a ca. 1-2 kcal mol-1 decrease in the aldehyde C-H and C-C bond energies corresponding to RC(.-)HCH(2)CH=O radical formation.     

Development of an ONIOM-G3B3 method to accurately predict C-H and N-H bond dissociation enthalpies of ribonucleosides and deoxyribonucleosides

The roles of ribonucleoside and deoxyribonucleoside radicals in DNA and RNA damage cannot be properly understood in the absence of knowledge of the C-H and N-H bond dissociation enthalpies (BDEs) depicting the energy cost to generate each of these radicals. However, because the nucleoside radicals tend to be extremely short-lived and it is very difficult to separate and identify different nucleoside radicals, experimental BDEs for nucleosides have remained elusive. Herein, we developed an ONIOM-G3B3 method in order to reliably predict the BDEs of nucleosides and we carefully benchmarked this new method against over 60 experimental BDEs of diverse sizable molecules. It was found that the accuracy of the ONIOM-G3B3 method was about 1.4 kcal/mol for BDE calculations. Using the ONIOM-G3B3 method, a full scale of C-H and N-H BDEs were obtained for the first time for ribonucleosides and deoxyribonucleosides with an estimated error bar of +/-1.4 kcal/mol. Discussions were then made about the interesting connections between these BDE values and previously reported experimental observations concerning radical-mediated DNA and RNA lesions. The significance of the work is twofold: (i) Nucleosides represent one of the most important groups of compounds in science. A full scale of reliable bond dissociation enthalpies for nucleosides is of fundamental importance. (ii) This work demonstrates the feasibility to accurately predict the bond strength of various sizable molecules ranging from nanosize molecular devices to biologically significant compounds.     

Thermochemistry of methyl and ethyl nitro, RNO2, and nitrite, RONO, organic compounds

Computational quantum theory is employed to determine the thermochemical properties of n-alkyl nitro and nitrite compounds: methyl and ethyl nitrites, CH3ONO and C2H5ONO, plus nitromethane and nitroethane, CH3NO2 and C2H5NO2, at 298.15 K using multilevel G3, CBS-QB3, and CBS-APNO composite methods employing both atomization and isodesmic reaction analysis. Structures and enthalpies of the corresponding aci-tautomers are also determined. The enthalpies of formation for the most stable conformers of methyl and ethyl nitrites at 298 K are determined to be -15.64 +/- 0.10 kcal mol-1 (-65.44 +/- 0.42 kJ mol-1) and -23.58 +/- 0.12 kcal mol-1 (-98.32 +/- 0.58 kJ mol-1), respectively. DeltafHo(298 K) of nitroalkanes are correspondingly evaluated at -17.67 +/- 0.27 kcal mol-1 (-74.1 +/- 1.12 kJ mol-1) and -25.06 +/- 0.07 kcal mol-1 (-121.2 +/- 0.29 kJ mol-1) for CH3NO2 and C2H5NO2. Enthalpies of formation for the aci-tautomers are calculated as -3.45 +/- 0.44 kcal mol-1 (-14.43 +/- 0.11 kJ mol-1) for aci-nitromethane and -14.25 +/- 0.44 kcal mol-1 (-59.95 +/- 1.84 kJ mol-1) for the aci-nitroethane isomers, respectively. Data are evaluated against experimental and computational values in the literature with recommendations. A set of thermal correction parameters to atomic (H, C, N, O) enthalpies at 0 K is developed, to enable a direct calculation of species enthalpy of formation at 298.15 K, using atomization reaction and computation outputs.     

N-H and alpha(C-H) bond dissociation enthalpies of aliphatic amines

Bond dissociation enthalpies (BDEs) of a large series of aliphatic amines (21) were measured by means of photoacoustic calorimetry. Despite the different structures studied in the primary, secondary, and tertiary amine series, the alpha(C-H) BDEs were found to be very similar for unconstrained amines with values very close to 91 kcal/mol. alphaC- and N-alkylation or introduction of an hydroxy group only slightly affect the BDEs, a fact in perfect agreement with calculations performed at different CBS levels. This demonstrates the predominance of the two-orbital-three-electron interaction involving the N and alphaC(*) orbitals. On the other hand, the N-H BDE decreases when going from primary to secondary amines. This result is interpreted in term of a hyperconjugation in sigmaC-C bonds, which leads to a stabilization of the aminyl radical. For cyclized amines, the BDEs depend on the relative geometry of the singly occupied alphaC(*) orbital with respect to that of the N atom, disfavoring the two-orbital-three-electron interaction. However, such structures can exhibit through-bond interaction. For a crowded structure such as triisopropylamine, for which the alphaC(*) orbital is not coplanar with the nitrogen one, the relaxation of a strain energy allows the BDE to be comparable to flexible structures.     

Substituent effects on the bond dissociation enthalpies of aromatic amines

Bond dissociation enthalpy differences, Z-X DeltaBDE = BDE(4-YC(6)H(4)Z-X) - BDE(C(6)H(5)Z-X), for Z = CH(2) and O are largely independent of X and are determined mainly by the stabilization/destabilization effect of Y on the 4-YC(6)H(4)Z(*) radicals. The effects of Y are small (< or =2 kcal/mol for all Y) for Z = CH(2), but they are large for Z = O, where good correlations with sigma(p)(+)(Y) yield rho(+) = 6.5 kcal/mol. For Z = NH, two sets of electrochemically measured N-H DeltaBDEs correlate with sigma(p)(+)(Y), yielding rho(+) = 3.9 and 3.0 kcal/mol. However, in contrast to the situation with phenols, these data indicate that the strengthening effect on N-H BDEs of electron-withdrawing (EW) Y's is greater than the weakening effect of electron-donating (ED) Y's. Attempts to measure N-H DeltaBDEs in anilines using two nonelectrochemical techniques were unsuccessful; therefore, we turned to density functional theory. Calculations on 15 4-YC(6)H(4)NH(2) gave N-H DeltaBDEs correlating with sigma(p)(+) (rho(+) = 4.6 kcal/mol) and indicated that EW and ED Y's had comparable strengthening and weakening effects, respectively, on the N-H bonds. To validate theory by connecting it to experiment, the N-H DeltaBDEs of four 4,4'-disubstituted diphenylamines and five 3,7-disubstituted phenothiazines were both calculated and measured by the radical equilibration EPR technique. For all compounds, theory and experiment agreed to better than 1 kcal/mol. Dissection of N-H DeltaBDEs in 4-substituted anilines and O-H DeltaBDEs in 4-substituted phenols into interaction enthalpies between Y and NH(2)/OH (molecule stabilization/destabilization enthalpy, MSE) and NH*/O* (radical stabilization/destabilization enthalpy, RSE) reveals that for both groups of compounds, ED Y's destabilize the molecule and stabilize the radical, while the opposite holds true for EW Y's. However, in the phenols the effects of substituents on the radical are roughly 3 times as great as those in the molecule, whereas in the anilines the two effects are of comparable magnitudes. These differences arise from the stronger ED character of NH(2) vs OH and the weaker EW character of NH* vs O*. The relatively large contributions to N-H BDEs in anilines arising from interactions in the molecules suggested that N-X DeltaBDEs in 4-YC(6)H(4)NH-X would depend on X, in contrast to the lack of effect of X on O-X and CH(2)-X DeltaBDEs in 4-YC(6)H(4)O-X and 4-YC(6)H(4)CH(2)-X. This suggestion was confirmed for X = CH(3), H, OH, and F, for which the calculated NH-X DeltaBDEs yielded rho(+) = 5.0, 4.6, 4.0, and 3.0 kcal/mol, respectively.     

Coupled cluster,B2PLYP and M06-2X investigation of the thermochemistry of five-membered nitrogen containing heterocycles,furan, and thiophene

Herein, the thermochemical properties of five-membered rings heterocycles were studied employing the CCSD(T) methodology coupled with the correlation consistent basis sets and including corrections for relativistic and core-valence effects as well as anharmonicities of the potentials. For pyrrole, furan, imidazole, pyrazole, 1H-1,2,4-triazole, and 1H-tetrazole, the mean absolute deviation (MAD) of the \Updelta \textH_{\textf, 2 9 8}^\texto \Updelta {\text{H}}_{{{\text{f}}, 2 9 8}}^{\text{o}} , computed at the CCSD(T) level, is 0.5 kcal/mol with respect to the experimental values. In the case of 1H-1,2,3-triazole, 2H-1,2,3-triazole, 4H-1,2,3-triazole, 4H-1,2,4-triazole, 2H-tetrazole, and pentazole, we propose the following \Updelta \textH_{\textf, 2 9 8}^\texto \Updelta {\text{H}}_{{{\text{f}}, 2 9 8}}^{\text{o}} : 62.6, 59.2, 85.0, 54.2, 77.7, and 107.5 kcal/mol, respectively. For thiophene, we revisit our previous result and propose a value of 26.0 kcal/mol. The theoretical estimations were used to study the performance of the M06-2X and B2PLYP functionals. Also, the convergence toward the complete basis set limit (CBS) was analyzed. M06-2X did not show a smooth convergence toward the CBS limit. Particularly, for the cc-pVTZ and cc-pVQZ basis sets, some problems were detected. Yet, along the cc-pVQZ, cc-pV5Z, and cc-pV6Z basis sets, the TAE smoothly decreased. The diminution of the TAE upon increase in basis set was not expected because the opposite behavior is more frequently observed. The MAD of the total atomization energies determined at the M06-2X level was 0.42 kcal/mol, with respect to the CCSD(T) results. In the case of the double hybrid B2PLYP functional, a smooth convergence toward the CBS limit was detected, even though the performance seriously degradated when the basis set was increased. At the CBS limit, the MAD with respect to the CCSD(T) TAEs was 8.26 kcal/mol.     

Is the 1,3,5-tridehydrobenzene triradical a cyclopropenyl radical analogue?

The gas-phase heat of formation (DeltaH(f,298)) of the 1,3,5-tridehydrobenzene triradical has been determined by using a negative ion thermochemical cycle. The first three measurements carried out were of the gas-phase acidity of 3,5-dichlorobenzoic acid, the enthalpy for decarboxylation of 3,5-dichlorobenzoate, and the enthalpy for chloride loss from 3,5,-dichlorophenide and constitute the measurement of the heat of formation for 5-chloro-m-benzyne. The last two measurements, the electron affinity of 5-chloro-m-benzyne, and the threshold for chloride loss from 5-chloro-m-benzyne, when combined with DeltaH(f,298) of 5-chloro-m-benzyne, give the heat of formation of the triradical. The 5-chloro-m-benzyne heat of formation is 116.2 +/- 3.7 kcal/mol. The heat of formation of the 1,3,5-tridehydrobenzene triradical measured in this work is 179.1 +/- 4.6 kcal/mol. This heat of formation was used to derive the bond dissociation energy (BDE) at the 5-position of m-benzyne, a third BDE in benzene. The BDE, at 109.2 +/- 5.6 kcal/mol, is ca. 4 kcal/mol lower than the first BDE in benzene (112.9 kcal/mol) and significantly higher than the BDE of phenyl radical at the meta position. The agreement between the first and third BDEs implies that the triradical is best described as a phenyl radical that interacts little with a m-benzyne moiety. The experimentally measured BDE is in good agreement with multireference configuration interaction calculations, which predict a (2)A(1) ground state for the Jahn-Teller distorted triradical. The trends in the first, second, and third BDEs of benzene are similar to those found for cyclopropane, suggesting a cyclopropenyl-like electronic structure within the six-membered ring of the 1,3,5-benzene triradical.     

C-H bond dissociation enthalpies in norbornane. An experimental and computational study

Gas-phase C-H bond dissociation enthalpies (BDEs) in norbornane were determined by quantum chemistry calculations and the C2-H BDE was experimentally obtained for the first time by time-resolved photoacoustic calorimetry. CBS-Q and CBS-QB3 methods were used to derive the values DH degrees (C1-H) = 449 kJ mol-1, DH degrees (C7-H) = 439 kJ mol-1, and DH degrees (C2-H) = 413 kJ mol-1. The experimental result DH degrees (C2-H) = 414.6 +/- 5.4 kJ mol-1 is in excellent agreement with the theoretical value. The trend DH degrees (C1-H) > DH degrees (C7-H) > DH degrees (C2-H) is discussed.     

Reaction of crown ethers with five-membered heterocycles: 1H-imidazole-4,5-dicarbonitrile, 3-nitro-1,2,4-triazole, and 1H-tetrazole

By the reaction of 1H-imidazole-4,5-dicarbonitrile, 3-nitro-1,2,4-triazole, and 1H-tetrazole, containing one proton at each of the two proton-donating atoms separated by one covalent bond, it was shown that uncharged crystalline supramolecular compounds can be obtained on account of the proton-donating characteristics of the aromatic heterocycles. __________ Translated from Khimiya Geterotsiklicheskikh Soedinenii, No. 9, pp. 1302–1305, September, 2006.     

Direct sp3 C-H bond activation adjacent to nitrogen in heterocycles

Activation of sp(3) C-H bonds adjacent to nitrogen in heterocycles is an attractive transformation that is emerging as a practical method in organic synthesis. This tutorial review aims to summarize the key examples of direct functionalization of nitrogen-containing heterocycles via metal-mediated and metal-catalyzed processes, which is meant to serve as a foundation for future investigations into this rapidly developing area of research. The review covers functionalization of N-heterocycles via alpha-lithiation with alkyllithium/diamine complexes, alpha-amino radical formation, metal-catalyzed direct C-H activation, C-H oxidations and oxidative couplings, and metal-catalyzed carbene insertions.     

Molecular energetics of cytosine revisited: a joint computational and experimental study

A static bomb calorimeter has been used to measure the standard molar energy of combustion, in oxygen, at T = 298.15 K, of a commercial sample of cytosine. From this energy, the standard (p degrees = 0.1 MPa) molar enthalpy of formation in the crystalline state was derived as -(221.9 +/- 1.7) kJ.mol(-1). This value confirms one experimental value already published in the literature but differs from another literature value by 13.5 kJ.mol(-1). Using the present standard molar enthalpy of formation in the condensed phase and the enthalpy of sublimation due to Burkinshaw and Mortimer [J. Chem. Soc., Dalton Trans. 1984, 75], (155.0 +/- 3.0) kJ.mol(-1), results in a value for the gas-phase standard molar enthalpy of formation for cytosine of -66.9 kJ.mol(-1). A similar value, -65.1 kJ.mol(-1), has been estimated after G3MP2B3 calculations combined with the reaction of atomization on three different tautomers of cytosine. In agreement with experimental evidence, the hydroxy-amino tautomer is the most stable form of cytosine in the gas phase. The enthalpies of formation of the other two tautomers were also estimated as -60.7 kJ.mol(-1) and -57.2 kJ.mol(-1) for the oxo-amino and oxo-imino tautomers, respectively. The same composite approach was also used to compute other thermochemical data, which is difficult to be measured experimentally, such as C-H, N-H, and O-H bond dissociation enthalpies, gas-phase acidities, and ionization enthalpies.     

Sulfur and selenium ylide bond enthalpies

The bond dissociation enthalpies (BDEs) of sulfur and selenium ylides have been estimated by applying MP2/6-311++G(3df,2p)//MP2/6-31G(d,p), G3, and other computational methods. Computed sulfoxide bond enthalpies were compared to experimental results to ensure the reliability of the computational methods before extending to related compounds. The examined ylides include the following: sulfoxides, sulfilimines, S,C-sulfonium ylides, and selenoxides. Selenoxides have BDEs about 10 kcal/mol smaller than the corresponding sulfoxides. N-H sulfilimines and CH2-S,C-sulfonium ylides have low BDEs, unless the sulfilimine or S,C-sulfonium ylide is stabilized by an electronegative substituent on N or C, respectively. Incorporation of the S or Se into a thiophene or selenophene-type ring lowers the BDE for the ylide.     

Homolytic C-H and N-H bond dissociation energies of strained organic compounds

High-level computations at G3, CBS-Q, and G3B3 levels were conducted, and good-quality C-H and N-H bond dissociation energies (BDEs) were obtained for a variety of saturated and unsaturated strained hydrocarbons and amines for the first time. From detailed NBO analyses, we found that the C-H BDEs of hydrocarbons are determined mainly by the hybridization of the parent compound, the hybridization of the radical, and the extent of spin delocalization of the radical. The ring strain has a significant effect on the C-H BDE because it forces the parent compound and radical to adopt certain undesirable hybridization. A structure-activity relationship equation (i.e., BDE (C-H) = 61.1-227.8 (p(parent)% - 0.75)(2) + 152.9 (p(radical)% - 1.00)(2) + 40.4 spin) was established, and it can predict the C-H BDEs of a variety of saturated and unsaturated strained hydrocarbons fairly well. For the C-H BDEs associated with the bridgehead carbons of the highly rigid strained compounds, we found that the strength of the C-H bond can also be predicted from the H-C-C bond angles of the bridgehead carbon. Finally, we found that N-H BDEs show less dependence on the ring strain than C-H BDEs.     

Modeling the competitive dissociation of protonated 2,3-butanedione. The enthalpy of formation of methylhydroxycarbene

The enthalpy of formation of methylhydroxycarbene, CH(3)COH, has been determined from measurements of the threshold energy for collision-induced dissociation of protonated 2,3-butanedione in a flowing afterglow-triple quadrupole mass spectrometer and found to be 16 +/- 4 kcal/mol, 57 +/- 4 kcal/mol higher than that of acetaldehyde. From the measured enthalpy of formation, the difference between the first and second C-H BDEs in ethanol is found to be 17 kcal/mol, which implies a singlet-triplet splitting of 28 kcal/mol in the carbene. The activation energies for loss of ketene and carbon monoxide from protonated butanedione are found to be 60 +/- 4 and 50 +/- 4 kcal/mol, respectively. On the basis of experimental and computational results, the loss of carbon monoxide is proposed to proceed through a tight transition state. Although calculations also suggest a tight transition state for loss of ketene, the experimental data indicate that it occurs via a loose transition state, possibly forming by proton transfer along the direct dissociation pathway.     

Energetic differences between the five- and six-membered ring hydrocarbons: strain energies in the parent and radical molecules

The C-H bond dissociation enthalpies (BDEs) for the five- and six-membered ring alkanes, alkenes, and dienes were investigated and discussed in terms of conventional strain energies (SEs). New determinations are reported for cyclopentane and cyclohexane by time-resolved photoacoustic calorimetry and quantum chemistry methods. The C-H BDEs for the alkenes yielding the alkyl radicals cyclopenten-4-yl and cyclohexen-4-yl and the alpha-C-H BDE in cyclopentene were also calculated. The s-homodesmotic model was used to determine SEs for both the parent molecules and the radicals. When the appropriate s-homodesmotic model is chosen, the obtained SEs are in good agreement with the ones derived from group additivity schemes. The different BDEs in the title molecules are explained by the calculated SEs in the parent molecules and their radicals: (1) BDEs leading to alkyl radicals are ca. 10 kJ mol (-1) lower in cyclopentane and cyclopentene than in cyclohexane and cyclohexene, due to a smaller eclipsing strain in the five-membered radicals relative to the parent molecules (six-membered hydrocarbons and their radicals are essentially strain free). (2) C-H BDEs in cyclopentene and cyclohexene leading to the allyl radicals are similar because cyclopenten-3-yl has almost as much strain as its parent molecule, due to a synperiplanar configuration. (3) The C-H BDE in 1,3-cyclopentadiene is 27 kJ mol (-1) higher than in 1,4-cyclohexadiene due to the stabilizing effect of the conjugated double bond in 1,3-cyclopentadiene and not to a destabilization of the cyclopentadienyl radical. The chemical insight afforded by group additivity methods in choosing the correct model for SE estimation is highlighted.     

Density Functional Calculations of C-NO2 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.     

The applicability of AM1 and PM3 semi-empirical methods for the study of N-H bond dissociation enthalpies and ionisation potentials of amine type antioxidants

A study of 40 para-substituted anilines is presented. These serve as model structures of amine type antioxidants. Molecules and their radical structures were studied using the AM1 and PM3 quantum chemical methods in order to calculate the N-H bond dissociation enthalpies (BDEs) and ionisation potentials (IPs) which are among the most important characteristics of antioxidants. Calculated BDEs were compared with available experimental values and the results of DFT calculations to ascertain the suitability of AM1 and PM3 methods for amine BDEs calculation. The results show that both methods reproduce experimental BDEs and DFT data satisfactorily. Comparison with experimental data shows that AM1 and PM3 methods overestimate the IP values. The results also indicate that dependences of N-H bond BDEs and IPs on Hammett constants of the substituents are linear.     

C-H bond strengths and acidities in aromatic systems: effects of nitrogen incorporation in mono-, di-, and triazines

The negative ion chemistry of five azine molecules has been investigated using the combined experimental techniques of negative ion photoelectron spectroscopy to obtain electron affinities (EA) and tandem flowing afterglow-selected ion tube (FA-SIFT) mass spectrometry to obtain deprotonation enthalpies (Δ(acid)H(298)). The measured Δ(acid)H(298) for the most acidic site of each azine species is combined with the EA of the corresponding radical in a thermochemical cycle to determine the corresponding C-H bond dissociation energy (BDE). The site-specific C-H BDE values of pyridine, 1,2-diazine, 1,3-diazine, 1,4-diazine, and 1,3,5-triazine are 110.4 ± 2.0, 111.3 ± 0.7, 113.4 ± 0.7, 107.5 ± 0.4, and 107.8 ± 0.7 kcal mol(-1), respectively. The application of complementary experimental methods, along with quantum chemical calculations, to a series of nitrogen-substituted azines sheds light on the influence of nitrogen atom substitution on the strength of C-H bonds in six-membered rings.