Homopolar- and Heteropolar Bond Dissociation Energies and Heats of Formation of Radicals and Ions in the Gas Phase. I. Data on organic molecules

The literature data on heteropolar and homopolar 2-center bond dissociation energies in organic molecules in the gas phase and the corresponding heats of formation of radicals and ions have been critically evaluated. Data for more than 500 bonds are represented in tabular form together with the pertinent literature references. Selected electron affinities and π-bond dissociation energies have also been incorporated. The follow-up paper will discuss some empirical general aspects of these data particularly regarding the effect of structure on the bond dissociation energies.     

Multi-structural thermodynamics of C-H bond dissociation in hexane and isohexane yielding seven isomeric hexyl radicals

The C-H bond dissociation processes of n-hexane and isohexane involve 23 and 13 conformational structures, respectively in the parent molecules and 14-45 conformational structures in each of the seven isomeric products that we studied. Here we use the recently developed multi-structural (MS) thermodynamics method and CCSD(T)-F12a/jul-cc-pVTZ//M06-2X/6-311+G(2df,2p) potential energy surfaces to calculate the enthalpy, entropy, and heat capacity of n-hexane, isohexane, and seven of the possible radical products of dissociation of C-H bonds. We compare our calculations with the limited experimental data and with values obtained by group additivity fits used to extend the experimental data. This work shows that using the MS method involving a full set of structural isomers with density functional geometries, scaled density functional frequencies, and coupled cluster single-point energies can predict thermodynamic functions of complex molecules and bond dissociation reactions with chemical accuracy. The method should be useful to obtain thermodynamic data for complex molecules for which such data has not been measured and to obtain thermodynamic data at temperatures outside the temperature range where measurements are available.     

Thermodynamic Properties of the Isomers of [HNOS], [HNO2S], and [HNOS2] and the Role of the Central Sulfur

Simple hydrides of compounds containing N, S, and O are of significant interest due to the role that they play in atmospheric chemistry and in biological pathways. There is a lack of quantitative thermodynamic data on these compounds. We have used a reliable computational chemistry approach based on valence CCSD(T) calculations extrapolated to the complete basis set limit with additional corrections to predict the heats of formation and bond dissociation energies of such compounds. The results show that compounds with the ability of the central S atom to effectively expand its valency leads to more stable isomers and, as a consequence, that those with the NSO structural motif are thermochemically more stable than those with the SNO motif. In addition, S?O bonds are preferred over N?O bonds.     

Calculation of bond dissociation energies for oxygen containing molecules by ab initio and density functional theory methods

Two ab initio (ROHF and MP2), one local (SVWN), four hybrid (BHandH, BHandHLYP, Becke3LYP, and Becke3P86), and two nonlocal (BLYP and BP86) density functional theory (DFT) methods are used for calculating the dissociation energies of molecules that contain H(SINGLE BOND)O, O(SINGLE BOND)O and O(SINGLE BOND)C bonds. The sensitivity to the basis set of the prediction of bond dissociation energies with DFT methods was tested with Becke3LYP on the H(SINGLE BOND)O dissociation energy of water. The 6–31 + G(d) methods are chosen as the smallest basis set which produces reasonable results. The calculated values for all other ab initio and DFT methods were performed with these basis sets and then compared with the experimental data. The suitability of DFT methods for computing reliable bond dissociation energies of oxygen containing molecules is discussed. © 1996 John Wiley & Sons, Inc.     

Some INDO calculations of 1J(PC) and 1J(PF) values

INDO parameterized calculations, employing phosphorus s, p and d valence orbitals, are reported for values of ¹J(PC) and ¹J(PF) relating to phosphorus in formal tri- and pentavalent states. The ¹J(PC), interactions are mainly controlled by the contact term. Thus, trivalent phosphorus compounds have negative values for ¹J(PC), whereas those for pentavalent phosphorus are positive due to the s lone-pair effect. The inclusion of phosphorus 3d orbitals is shown to be important for an understanding of the processes contributing to ¹J(PF) interactions. ¹J(PF) values are shown to be negative for both tri- and pentavalent phosphorus compounds. The contact and orbital interactions are significant for the trivalent phosphorus molecules, whereas in the pentavalent phosphorus case ¹J(PF) is dominated by the orbital term.     

Comparison between 6-31 + bond functions and 6-31 G* basis sets in UMP calculations of dissociation energies of the first row element compounds. I

A basis set with bond functions (6-31G + BF) has been tested for its applicability to calculation of dissociation energies of single and multiple bonds by Moeller-Plesset perturbation theory at the second and third orders. Results have been compared with those calculated in the 6-31G* basis set. The 6-31G + BF basis at the MP 2 and MP 3 levels yields better results than 6-31G* basis and the time consumption is less as well. Consideration of the bond functions on the bonds neighboring the bond being broken has no significant influence on the dissociation energies either at the SCF or at the MP 2 levels. If both reactants and products can be characterized by two-center bonds, the 6-31G + BF basis and UMP2 variant of perturbation theory can be recommended for practical calculation of D_e values, especially for the systems where the use of more exact bases is rather difficult.     

Unusual structural and energetic features of homolytic bond dissociation: from tetrakis(disyl)diphosphine to tetrakis(di-tert-butylsilyl)hydrazine

Quantum chemical calculations of the structures and thermodynamics of homolytic dissociation of the central P-P and N-N bonds in tetrakis(disyl)diphosphine and tetrakis(di-tert-butylsilyl)hydrazine have been performed. The theory predicted negative standard enthalpies for homolytic bond dissociation in both cases, -71.0 and -108.4 kJ mol(-1) for the diphosphine and hydrazine, respectively, using the ONIOM (MP2/6-31+G*:B3LYP/3-21G*) level. The dissociation is accompanied by considerable structural changes in the radicals as compared to the corresponding fragments of the parent molecules, resulting in low dissociation enthalpies. The most pronounced changes in both radicals are the relaxation of bond angles in the substituents and a conformational change in the orientation of the substituent groups. In addition, the bis(di-tert-butylsilyl)aminyl radical displays a considerable increase in Si-N-Si angle and shortening of the Si-N bonds upon dissociation. These changes are not associated with any appreciable delocalisation of the lone electron, as the spin density is found from the B3LYP/3-21G* calculations to be largely concentrated on the nitrogen atom. It has been also shown that although the dissociation energies are low for both compounds, the intrinsic energies of the central bonds are still high, 140.6 kJ mol(-1) for the P-P bond in tetrakis(disyl)diphosphine and 490.6 kJ mol(-1) for the N-N bond in tetrakis(di-tert-butylsilyl)hydrazine, using the ONIOM method. The calculations predict that the dissociation of tetrakis(disyl)diphosphine would have negative free energy even without taking relaxation of the fragments into account, while the full potential of releasing about 306 kJ mol(-1) of energy stored in the ligands of tetrakis(di-tert-butylsilyl)hydrazine is only fully realised upon a considerable separation of the fragments.     

Damage to model DNA fragments from very low-energy (<1 eV) electrons

Although electrons having enough energy to ionize or electronically excite DNA have long been known to cause strand breaks (i.e., bond cleavages), only recently has it been suggested that even lower-energy electrons (most recently 1 eV and below) can also damage DNA. The findings of the present work suggest that, while DNA bases can attach electrons having kinetic energies in the 1 eV range and subsequently undergo phosphate-sugar O-C sigma bond cleavage, it is highly unlikely (in contrast to recent suggestions) that electrons having kinetic energies near 0 eV can attach to the phosphate unit's P=O bonds. Electron kinetic energies in the 2-3 eV range are required to attach directly to DNA's phosphate group's P=O pi orbital and induce phosphate-sugar O-C sigma bond cleavages if the phosphate groups are rendered neutral (e.g., by nearby counterions). Moreover, significant activation barriers to C-O bond breakage render the rates of both such damage mechanisms (i.e., P=O-attached and base-attached) slow as compared to electron autodetachment and to other damage processes.     

Strong intramolecular hydrogen bonds and steric effects involving C═S groups: An NMR and computational study

A number of 5-acyl rhodanines and thiorhodanines with bulky acyl groups (pivaloyl and adamantoyl), not previously available, have been synthesized. The compounds are shown to exist in the enol form. Structures have been calculated using both the MP2 approach and the B3LYP-GD3BJ functional and the 6-311++G(d,p) basis set. Hydrogen bond energies are estimated by subtracting energies of a structure with the OH group turned 180° from those of the intramolecularly hydrogen-bonded one. Properties such as OH chemical shifts, two-bond isotope effects on ¹³C chemical shifts, electron densities at the bond critical point from atoms in molecules analysis, and the hydrogen bond energies show that the sterically hindered compounds have stronger hydrogen bonds than methyl or isopropyl derivatives. The combination of oxygen and sulfur derivatives enables a detailed analysis of hydrogen bond energies.     

Interaction of α,α-difluoroazides with trivalent phosphorus compounds and triphenylantimony

α,α-Difluoroazides react with triphenylantimony and various compounds of trivalent phosphorus according to the oxidative fluorination scheme. In the case of trivalent phosphorus compound the primary products are phosphazenes, phosphazides or difluorophosphoranes that may undergo further transformations to the corresponding fluorine derivatives of pentavalent phosphorus.     

Can mass dissociation patterns of transition‐metal complexes be predicted from electrochemical data?

The Cooks kinetic method has been very convenient to correlate the relative dissociation rates obtained by collision-induced fragmentation experiments with the energies of two related bonds in molecules and complexes in the gas phase. Reliable bond energy data are, however, not always available, particularly for polynuclear transition-metal complexes, such as the triruthenium acetate clusters of the general formula [Ru(3) (micro(3)-O)(micro-CH(3)COO)(6)(py)(2)(L)](+), where L = ring substituted N-heterocyclic ligands. Accordingly, their gas-phase collision-induced tandem mass spectrometry (CID MS/MS) dissociation patterns have been analyzed pursuing a relationship with the more easily accessible redox potentials (E(1/2)) and Lever's E(L) parameters. In fact, excellent linear correlations of ln(1/2A(L)/A(py)), where A(py) and A(L) are the abundance of the fragments retaining the pyridine (py) and L ligand, respectively, with E(1/2) and E(L) were found. This result shows that those electrochemical parameters are correlated with bond energies and can be used in the analysis of the dissociation data. Such modified Cooks method can be used, for example, to determine the electronic effects of substituents on the metal-ligand bonds for a series of transition-metal complexes.     

Adducts of dicoordinated trivalent phosphorus with carbenes

1,4-Dichloro-3a,6a-diaza-1,4-diphosphapentalene reacts with sodium acetylacetonate to form an adduct of diacylcarbene with trivalent dicoordinated phosphorus. In solution this adduct dimerizes according to the head-to-head type. The phosphorus atoms undergo transformation from the hypervalent state (in the adduct), in which the lone electron pair on the phosphorus atom is not involved in the formation of additional bonds with this phosphorus(III) atom, to the pentavalent pentacoordinated and trivalent tricoordinated state (in the dimer).     

分子内结构环境对解离能的影响

为了理解化学键的这一结构效应, 本文对具有相同化学键而分子内结构环境不同的系列分子进行了计算研究, 讨论了化学键结构环境对解离能的影响.     

第一过渡系金属硅烯配合物的理论研究

以HF/6-311+G^*基组研究了硅烯SiH₂同第一过渡系金属的配合物MSiH₂的分子轨道特征及键解离能.MSiH₂为共平面构型.其中基态的3TiSiH₂和4CoSiH₂带有明显的双键特征.M-Si键具有共价性质.M-Si的键解离能,从Sc到Cu呈现周期性变化,这种变化趋势同M的金属离子激发能之间存在近似的线性关系.     

Density Functional Theory Study of the Interaction between Guanine and Catechin

The interacting patterns and mechanism of the catechin and guanine have been investigated with the density functional theory B3LYP method by 6‐31G* basis set. Fourteen stable structures for the catechin‐guanine complexes have been found which form two hydrogen bonds at least. The results indicate that the complexes are mainly stabilized by the hydrogen bonding interactions. At the same time, the number and strength of hydrogen bond play a co‐determinant parts in the stability of the complexes which can form two or more hydrogen bonds. Theories of atoms in molecules (AIM) and natural bond orbital (NBO) have been adopted to investigate the hydrogen bonds involved in all systems. The interaction energies of all complexes have been corrected for basis set superposition error (BSSE), ranging from ?38.86 to ?14.56 kJ/mol. The results showed that the hydrogen bonding contributes to the interaction energies dominantly. The corresponding bonds stretching motions in all complexes are red‐shifted relative to that of the monomer, which is in agreement with experimental results.     

Bond Energy Scheme for Estimating Heats of Formation of Monomers and Polymers. VI. Sulfur Compounds

The bond energy scheme is extended to sulfur compounds and heats of formation and atomization energy terms derived from thermochemical data reviewed to 1977, for bonds of sulfur with carbon, hydrogen, halogens, and oxygen atoms. A precision of ± 1 kcal/mole was attainable for the covalent bonds of divalent sulfur in the lowest oxidation state S(± II). The higher valency states: S(IV) and S(VI) involve polar contributions depending upon the electrouegativity of the combining atom as well as (d^π -p^π) orbital promotion energies which are specific to the compound and transferable to other molecules only with a limited precision, no better than about ± 3 kcal/mole. The atomization energy terms (E_a ^25°C) of various bonds of sulfur a are found consistent with the experimental bond dissociation energies and bear a relationship with bond lengths and force constants as observed in the previous work. Heats of polymer-forming reactions and heats of formation of sulfur-containing monomers and polymers are estimated from the newly derived bond energy terms.     

The thermochemistry of group 15 tetrachloride anions

Bond strengths for a series of Group 15 tetrachloride anions ACl4 (A = P, As, Sb, and Bi) have been determined by measuring thresholds for collision-induced dissociation of the anions in a flowing afterglow-tandem mass spectrometer. The central atoms in these systems have ten electrons, which violates the octet rule: the bond dissociation energies for ACl4- help to clarify the effect of the central atom on hypervalent bond strengths. The 0 K bond energies in kJ mol(-1) are D(Cl3A-CL-) = 90 +/- 7,115 +/- 7,161 +/- 8, and 154 +/- 15, respectively. Computational results using the B3LYP/LANL2DZpd level of theory are higher than the experimental bond energies. Calculations give a geometry for BiCl4 that is essentially tetrahedral rather than the see-saw observed for the other tetrachlorides. NBO calculations predict that the phosphorus and arsenic systems have 3C-4E bonds, while the antimony and bismuth systems are more ionic.     

On a definition of bond energies based on localized molecular orbitals

A theoretical model is presented for defining bond energies based on localized molecular Orbitals. These bond energies are obtained by rearranging the total SCF energy including the nuclear repulsion term to a sum over orbital and orbital interaction terms and then to total orbital terms, which can be interpreted as the energies of localized orbitals in a molecule. A scaling procedure is used to obtain a direct connection with experimental bond dissociation energies. Two scale parameters are employed, the C-C and the C-H bond dissociation energy in C₂H₆ for A-B and C-H type bonds, respectively. The implications of this scaling procedure are discussed. Numerical applications to a number of organic molecules containing no conjugated bonds gives in general a very satisfactory agreement between experimental and theoretical bond energies.     

Strain Energies in Homoatomic Nitrogen Clusters N(4), N(6), and N(8)

Strain energies and resonance energies can be obtained as the energy changes for appropriate homodesmotic reactions using ab initio calculated total energies as the energies of the reactants and products involved. Homodesmotic reactions conserve bond types and preserve valence environments at all atoms, requirements that favor the cancellation of basis set and electron correlation errors in the ab initio energies. In this paper we calculate strain energies and resonance energies for N(4), N(6), and N(8) clusters in a number of chemically significant but, for nitrogen, hypothetical structural forms. The nitrogen cluster strain energies are generally of the same order of magnitude as those of isostructural hydrocarbon clusters, and individual differences can be explained by using the ring strain additivity rule and recognizing the effect of the presence of lone pairs of electrons on nitrogen clusters but not on the hydrocarbons. Resonance energies of the nitrogen clusters are much smaller than those of the comparable aromatic hydrocarbons. The differences can be rationalized by considering the relative strengths of CC and NN single and double bonds. Strain and resonance energies of nitrogen clusters are compared with those previously reported for homoatomic clusters of phosphorus and arsenic. Trends through the series are remarkably similar, but strain energies for clusters from lower periods are progressively smaller. Strain and resonance have been important organizing concepts in organic chemistry for many years. Estimates of corresponding parameters for inorganic analogs are only now becoming available.