Dimers of Nineteen‐Electron Sandwich Compounds: Crystal and Electronic Structures,and Comparison of Reducing Strengths

The dimers of some Group 8 metal cyclopentadienyl/arene complexes and Group 9 metallocenes can be handled in air, yet are strongly reducing, making them useful n‐dopants in organic electronics. In this work, the X‐ray molecular structures are shown to resemble those of Group 8 metal cyclopentadienyl/pentadienyl or Group 9 metal cyclopentadienyl/diene model compounds. Compared to those of the model compounds, the DFT HOMOs of the dimers are significantly destabilized by interactions between the metal and the central C?C σ‐bonding orbital, accounting for the facile oxidation of the dimers. The lengths of these C?C bonds (X‐ray or DFT) do not correlate with DFT dissociation energies, the latter depending strongly on the monomer stabilities. Ru and Ir monomers are more reducing than their Fe and Rh analogues, but the corresponding dimers also exhibit much higher dissociation energies, so the estimated monomer cation/neutral dimer potentials are, with the exception of that of [RhCp₂]₂, rather similar (?1.97 to ?2.15 V vs. FeCp₂^+/0 in THF). The consequences of the variations in bond strength and redox potentials for the reactivity of the dimers are discussed.     

n‐Doping of Organic Electronic Materials Using Air‐Stable Organometallics: A Mechanistic Study of Reduction by Dimeric Sandwich Compounds

Several 19‐electron sandwich compounds are known to exist as “2×18‐electron” dimers. Recently it has been shown that, despite their air stability in the solid state, some of these dimers act as powerful reductants when co‐deposited from either the gas phase or from solution and that this behavior can be useful in n‐doping materials for organic electronics, including compounds with moderate electron affinities, such as 6,13‐bis[tri(isopropyl)silylethynyl]pentacene ( 3 ). This paper addresses the mechanisms by which the dimers of 1,2,3,4,5‐pentamethylrhodocene ( 1 b₂ ), (pentamethylcyclopentadienyl)(1,3,5‐trialkylbenzene)ruthenium (alkyl=Me, 2 a₂ ; alkyl=Et, 2 b₂ ), and (pentamethylcyclopentadienyl)(benzene)iron ( 2 c₂ ) react with 3 in solution. Vis/NIR and NMR spectroscopy, and X‐ray crystallography indicate that the products of these solution reactions are 3 ^.? salts of the monomeric sandwich cations. Vis/NIR kinetic studies for the Group 8 dimers are consistent with a mechanism whereby an endergonic electron transfer from the dimer to 3 is followed by rapid cleavage of the dimer cation. NMR crossover experiments with partially deuterated derivatives suggest that the C? C bond in the 1 b₂ dimer is much more readily broken than that in 2 a₂ ; consistent with this observation, Vis/NIR kinetic measurements suggest that the solution reduction of 3 by 1 b₂ can occur by both the mechanism established for the Group 8 species and by a mechanism in which an endergonic dissociation of the dimer is followed by rapid electron transfer from monomeric 1 b to 3 .     

Why Does the Intramolecular Trimethylene‐Bridged Frustrated Lewis Pair Mes2PCH2CH2CH2B(C6F5)2 Not Activate Dihydrogen?

The methyl labelled C₃‐bridged frustrated phosphane borane Lewis pair (P/B FLP) 2 b was prepared by treatment of Mes₂PCl with a methallyl Grignard reagent followed by anti‐Markovnikov hydroboration with Piers’ borane [HB(C₆F₅)₂)]. The FLP 2 b is inactive toward dihydrogen under typical ambient conditions, in contrast to the C₂‐ and C₄‐bridged FLP analogues. Dynamic NMR spectroscopy showed that this was not due to kinetically hindered P???B dissociation of 2 b . DFT calculations showed that the hydrogen‐splitting reaction of the parent compound 2 a is markedly endergonic. The PH⁺/BH^? H₂‐splitting product of 2 b was indirectly synthesized by a sequence of H⁺/H^? addition. It lost H₂ at ambient conditions and confirmed the result of the DFT analysis.     

Synthesis and structural study of organic two‐photon‐absorbing cycloalkanone chromophores

The three organic two‐photon‐absorbing cycloalkanone chromophores 2,4‐bis[4‐(diethylamino)benzylidene]cyclobutanone, C₂₆H₃₂N₂O ( I ), 2,5‐bis[4‐(diethylamino)benzylidene]cyclopentanone, C₂₇H₃₄N₂O ( II ), and 2,6‐bis[4‐(diethylamino)benzylidene]cyclohexanone, C₂₈H₃₆N₂O ( III ), were obtained by a reaction between 4‐(diethylamino)benzaldehyde and the corresponding cycloalkanone and were characterized by single‐crystal X‐ray diffraction studies, as well as density functional theory (DFT) quantum‐chemical calculations. Molecules of this series have three main fragments, i.e. central acceptor (A) and two terminal donors (D₁ and D₂) and represent examples of the D₁–π–A–π–D₂ molecular design. All three compounds crystallize with two crystallographically independent molecules in the asymmetric unit ( A and B ) and are distinguished by the conformations of both the molecular Et₂N—C₆H₄—C=C—C(=O)—C=C—C₆H₄—NEt₂ backbone (arcuate or linear) and the terminal diethylamino substituents (syn‐ or antiperiplanar to the plane of the molecule). The central four‐ and five‐membered rings in I and II are almost planar, and the six‐membered ring in III adopts a sofa conformation. In the crystals of I – III , the two independent molecules A and B form hydrogen‐bonded [ A … B ] dimers via intermolecular C—H…O hydrogen bonds. Furthermore, the [ A … B ] dimers in I are bound by intermolecular C—H…O hydrogen bonds into two‐tier puckered layers, whereas in the crystals of II and III , the [ A … B ] dimers are stacked along the c and a axes, respectively. Taking into account the decreasing steric strain upon expanding the central ring, compound I might be more efficient as a two‐photon absorption chromophore than compounds II and III , which corresponds to the results of spectroscopic studies.     

Theoretical Study of Weakly Bound Dimers between Hydrogen Fluoride and Some Polar Molecules

Weakly bound linear and bent dimers, FH—X (where X = CO, OC, CNH, NCH, N₂O and ON₂), are investigated using the DFT B3LYP and ab initio MP2 methods with the same basis sets (6–311++G(3df,2pd)). The strengths of the H—C or H—N H‐bonds in dimers FH—CO, FH—CNH, and FH—N₂O are compared with those of the H—O or H—N H‐bonds in dimers FH—OC, FH—NCH, and FH—ON₂. The results obtained for the H‐bond distances, the elongation effect of the HF bond, the red shift of the HF stretching frequency, and the energy difference between the dimer and the charge transfer reveal that the H‐bonds of the first group of dimers are stronger than those of the second. The Gibbs energies calculated for the six dimer formations indicate that the weakly bound dimers are unstable at room temperature (T = 298 K) (FH—X's → FH + X's, ΔG < 0).     

An Adaptable N‐Heterocyclic Carbene Macrocycle Hosting Copper in Three Oxidation States

A neutral hybrid macrocycle with two trans‐positioned N‐heterocyclic carbenes (NHCs) and two pyridine donors hosts copper in three oxidation states (+I–+III) in a series of structurally characterized complexes ( 1 – 3 ). Redox interconversion of [LCu]^+/2+/3+ is electrochemically (quasi)reversible and occurs at moderate potentials (E_1/2=?0.45 V and +0.82 V (vs. Fc/Fc⁺)). A linear C^NHC‐Cu‐C^NHC arrangement and hemilability of the two pyridine donors allows the ligand to adapt to the different stereoelectronic and coordination requirements of Cu^I versus Cu^II/Cu^III. Analytical methods such as NMR, UV/Vis, IR, electron paramagnetic resonance, and Cu Kβ high‐energy‐resolution fluorescence detection X‐ray absorption spectroscopies, as well as DFT calculations, give insight into the geometric and electronic structures of the complexes. The XAS signatures of 1 – 3 are textbook examples for Cu^I, Cu^II, and Cu^III species. Facile 2‐electron interconversion combined with the exposure of two basic pyridine N sites in the reduced Cu^I form suggest that [LCu]^+/2+/3+ may operate in catalysis via coupled 2 e^?/2 H⁺ transfer.     

Implications of Nitrogen Doping on Geometrical and Electronic Structure of the Fullerene Dimers

Because of its unsaturated bonds, C₆₀ is susceptible to polymerize into dimers. The implications of nitrogen doping on the geometrical and electronic structure of C₆₀ dimers have been ambiguous for years. A quarter‐century after the discovery of azafullerene dimer (C₅₉N)₂, we reported its single crystallographic structure in 2019. Herein, the unambiguous crystal structure information of (C₅₉N)₂ is elucidated specifically, revealing that the inter‐cage C—C single bond length of (C₅₉N)₂ is comparable with that of an ordinary C(sp³)‐C(sp³) single bond, and that the most stable conformer of (C₅₉N)₂ is gauche‐conformer with a dihedral angle of 66°. To amend the structural deviations, geometrical structure of (C₅₉N)₂ is optimized by a B3LYP‐D3BJ function, which is proved to be more consistent with its single crystal structure than those by the commonly used B3LYP function. Moreover, the calculation method is also suitable for other representative fullerene dimers, such as (C₆₀)₂ and its divalent anion. Additionally, the dissociation of (C₅₉N)₂ at 473 K under mass spectrometric conditions suggests the inter‐cage C—C bond is relatively weaker than an ordinary C—C single bond, which can be explained by the interaction energies of inter‐cages.     

CO Blocking of D2 Dissociative Adsorption on Ru(0001)

The influence of pre‐adsorbed CO on the dissociative adsorption of D₂ on Ru(0001) is studied by molecular‐beam techniques. We determine the initial dissociation probability of D₂ as a function of its kinetic energy for various CO pre‐coverages between 0.00 and 0.67 monolayers (ML) at a surface temperature of 180 K. The results indicate that CO blocks D₂ dissociation and perturbs the local surface reactivity up to the nearest‐neighbour Ru atoms. Non‐activated sticking and dissociation become less important with increasing CO coverage, and vanish at θ_CO≈0.33 ML. In addition, at high D₂ kinetic energy (>35 kJ mol^?1) the site‐blocking capability of CO decreases rapidly. These observations are attributed to a CO‐induced activation barrier for D₂ dissociation in the vicinity of CO molecules.     

Structural,electronic, and magnetoresponsive properties of triangular lanthanide clusters and their free‐standing nitrides

The molecular and electronic structures, stabilities, bonding features, and magnetoresponsive properties of three‐membered [c‐Ln₃]^+/0/? (Ln = La, Ce, Pr, Nd, Gd, Lu) and heterocyclic six‐membered [c‐Ln₃E₃]^q (Ln = La, Ce, Pr, Nd, Gd, Lu; E = C, N; q = 0 or 1) rings have been investigated by means of electronic structure calculation methods at the DFT level. The [c‐Ln₃]^+/0/? clusters are predicted to be bound with respect to dissociation to their constituent atoms, the estimated binding energies ranging from 45.8 to 2056.4 kJ/mol. The [c‐Ln₃] rings capture easily a planar three‐coordinated nitrogen atom at the center or above the center of the ring yielding the lanthanide nitride clusters [c‐Ln₃(μ₃‐N)] adopting a planar geometry, except [c‐La₃(μ₃‐N)] which exhibits pyramidal geometry. The [c‐Ln₃(μ₃‐N)] clusters are predicted to be bound, with respect to dissociation to N (⁴S) atom and [c‐Ln₃] clusters in their ground states, the binding energies ranging from 53.9 to 257.9 kcal/mol. The six‐membered [c‐Ln₃E₃]^q rings are predicted to be bound with respect to dissociation to LnE^q monomers in their ground states with dissociation energies in the range of 173.8 to 318.0 kcal/mol. Calculation of the NICS_zz‐scan curves of the clusters predicted a “hermaphrodic” magnetic response of the [c‐Ln₃]^+/0/? and heterocyclic six‐membered [c‐Ln₃E₃]^q rings, manifested by the coexistence of successive diatropic (aromatic) and paratropic (antiaromatic) zones. The [c‐La₃]^+/0/? and [c‐Lu₃]^? are predicted to be weakly antiaromatic, the [c‐Lu₃]^0/+, [c‐Lu₃C₃]⁺, and [c‐Lu₃N₃] double (σ+π) aromatic, and the [c‐Gd₃C₃] and [c‐Gd₃N₃]⁺ rings (σ+δ)‐aromatic systems. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2011     

Relativistic and Non‐Relativistic Electronic Molecular‐Structure Calculations for Dimers of 4p‐, 5p‐, and 6p‐Block Elements

We report results of non‐relativistic and two‐component relativistic single‐reference coupled‐cluster with single and double and perturbative triple excitations [CCSD(T)] treatments for the 4p‐block dimers Ga₂ to Br₂, the 5p‐block dimers In₂ to I₂, and their atoms. Extended basis sets up to pentuple zeta are employed and energies extrapolated to the complete basis‐set limit. Relativistic and non‐relativistic results for the dissociation energy D_e are in close agreement with each other and previously published data, provided non‐relativistic or scalar‐relativistic results are corrected for spin–orbit contributions taken from the literature. An exception is Te₂ where theoretical results scatter by 0.085 eV. By virtue of this agreement it is unexpected that comparison with the experimental D₀ or D_e dissociation energies (zero‐point vibrational effects are negligible in this context) reveal errors larger than 0.1 eV for Ga₂, Ge₂, and Sb₂. Only relativistic treatments are presented for the 6p‐block cases Tl₂ to At₂. Sufficient agreement with experimental data is found only for Pb₂ and Bi₂, the deviation of the computed and experimental D₀ values for Po₂ is again larger than 0.1 eV. Deviations of 0.1 eV between the computed and experimental D₀ values are a major reason for concern and call for additional investigations in both fields to clarify the situation.     

The confirmation of accurate combination of functional and basis set for transition‐metal dimers: Fe2, Co2, Ni2, Ru2, Rh2, Pd2, Os2, Ir2, and Pt2

Eight kinds of density functionals named B3LYP, PBE1PBE, B1B95, BLYP, BP86, G96PW91, mPWPW91, and SVWN along with two different valence basis sets (LANL2DZ and CEP‐121g) are employed to study the transition‐metal dimers for the elements of group VIII. By comparing the equilibrium bond distances, vibrational frequencies, and dissociation energies of the ground state of these dimers with the available experimental values and theoretical data, we show that the “pure” DFT methods (G96PW91, BLYP, and BP86) with great‐gradient approximation always give better results relative to the hybrid HF/DFT schemes (B3LYP, PBE1PBE, and B1B95). The striking case found by us is that the G96PW91 functional, which is not tested in previous systemic studies, always predicts the dissociation energy to be well. The Ru₂ and Os₂ dimers are sensitive to not only the functionals employed but also the valence basis sets adopted. The natural bond orbital population is analyzed, and the molecular orbitals of the unpaired electrons are determined. Furthermore, our results indicate that the s and d orbitals of these dimers always hybridize with each other except for Rh₂ and Pt₂ molecules. And by analyzing the electron configuration of the bonding atom, the dissociation limit of the ground state is obtained. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Reinterpretation of Dynamic Vibrational Spectroscopy to Determine the Molecular Structure and Dynamics of Ferrocene

Molecular distortion of dynamic molecules gives a clear signature in the vibrational spectra, which can be modeled to give estimates of the energy barrier and the sensitivity of the frequencies of the vibrational modes to the reaction coordinate. The reaction coordinate method (RCM) utilizes ab initio‐calculated spectra of the molecule in its ground and transition states together with their relative energies to predict the temperature dependence of the vibrational spectra. DFT‐calculated spectra of the eclipsed (D_5h) and staggered (D_5d) forms of ferrocene (Fc), and its deuterated analogue, within RCM explain the IR spectra of Fc in gas (350 K), solution (300 K), solid solution (7–300 K), and solid (7–300 K) states. In each case the D_5h rotamer is lowest in energy but with the barrier to interconversion between rotamers higher for solution‐phase samples (ca. 6 kJ mol^?1) than for the gas‐phase species (1–3 kJ mol^?1). The generality of the approach is demonstrated with application to tricarbonyl(η⁴‐norbornadiene)iron(0), Fe(NBD)(CO)₃. The temperature‐dependent coalescence of the ν(CO) bands of Fe(NBD)(CO)₃ is well explained by the RCM without recourse to NMR‐like rapid exchange. The RCM establishes a clear link between the calculated ground and transition states of dynamic molecules and the temperature‐dependence of their vibrational spectra.     

Formation and Characterisation of the Silver Hydride Nanocluster Cation [Ag3H2((Ph2P)2CH2)]+ and Its Release of Hydrogen

Multistage mass spectrometry and density functional theory (DFT) were used to characterise the small silver hydride nanocluster, [Ag₃H₂L]⁺ (where L=(Ph₂P)₂CH₂) and its gas‐phase unimolecular chemistry. Collision‐induced dissociation (CID) yields [Ag₂HL]⁺ as the major product while laser‐induced dissociation (LID) proceeds via H₂ formation and subsequent release from [Ag₃H₂L]⁺, giving rise to [Ag₃L]⁺ as the major product. Deuterium labelling studies on [Ag₃D₂L]⁺ prove that the source of H₂ is from the hydrides and not from the ligand. Comparison of TD‐DFT absorption patterns obtained for the optimised structures with action spectroscopy results, allows assignment of the measured features to structures of precursors and products. Molecular dynamics “on the fly” reveal that AgH loss is favoured in the ground state, but H₂ formation and loss is preferred in the first excited state S₁, in agreement with CID and LID experimental findings. This indicates favourable photo‐induced formation of H₂ and subsequent release from [Ag₃H₂L]⁺, an important finding in context of metal hydrides as a hydrogen storage medium, which can subsequently be released by heating or irradiation with light.     

A CAS study on S‐loss and O‐loss dissociation mechanisms of the SO2+ ion in the C,D, and E states

In the present work, we mainly study dissociation of the C ²B₁, D²A₁, and E²B₂ states of the SO₂⁺ ion using the complete active‐space self‐consistent field (CASSCF) and multiconfiguration second‐order perturbation theory (CASPT2) methods. We first performed CASPT2 potential energy curve (PEC) calculations for S‐ and O‐loss dissociation from the X, A, B, C, D, and E primarily ionization states and many quartet states. For studying S‐loss predissociation of the C, D, and E states by the quartet states to the first, second, and third S‐loss dissociation limits, the CASSCF minimum energy crossing point (MECP) calculations for the doublet/quartet state pairs were performed, and then the CASPT2 energies and CASSCF spin‐orbit couplings were calculated at the MECPs. Our calculations predict eight S‐loss predissociation processes (via MECPs and transition states) for the C, D, and E states and the energetics for these processes are reported. This study indicates that the C and D states can adiabatically dissociate to the first O‐loss dissociation limit. Our calculations (PEC and MECP) predict a predissociation process for the E state to the first O‐loss limit. Our calculations also predict that the E²B₂ state could dissociate to the first S‐ and O‐loss limits via the A²B₂ ← E²B₂ transition. On the basis of the 13 predicted processes, we discussed the S‐ and O‐loss dissociation mechanisms of the C, D, and E states proposed in the previous experimental studies. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Synthesis,Structure, and Reactivity of Anionic sp2–sp3 Diboron Compounds: Readily Accessible Boryl Nucleophiles

Lewis base adducts of tetra‐alkoxy diboron compounds, in particular bis(pinacolato)diboron (B₂pin₂), have been proposed as the active source of nucleophilic boryl species in metal‐free borylation reactions. We report the isolation and detailed structural characterization (by solid‐state and solution NMR spectroscopy and X‐ray crystallography) of a series of anionic adducts of B₂pin₂ with hard Lewis bases, such as alkoxides and fluoride. The study was extended to alternative Lewis bases, such as acetate, and other diboron reagents. The B(sp²)–B(sp³) adducts exhibit two distinct boron environments in the solid‐state and solution NMR spectra, except for [(4‐tBuC₆H₄O)B₂pin₂]^?, which shows rapid site exchange in solution. DFT calculations were performed to analyze the stability of the adducts with respect to dissociation. Stoichiometric reaction of the isolated adducts with two representative series of organic electrophiles—namely, aryl halides and diazonium salts—demonstrate the relative reactivities of the anionic diboron compounds as nucleophilic boryl anion sources.     

Synthesis,crystal structure and magnetic properties of (acetato‐κ2O,O′)bis(5,5′‐dimethyl‐2,2′‐bipyridine‐κ2N,N′)nickel(II) perchlorate monohydrate

The title hydrated ionic complex, [Ni(CH₃COO)(C₁₂H₁₂N₂)₂]ClO₄·H₂O or [Ni(ac)(5,5′‐dmbpy)₂]ClO₄·H₂O (where 5,5′‐dmbpy is 5,5′‐dimethyl‐2,2′‐bipyridine and ac is acetate), (1), was isolated as violet crystals from the aqueous ethanolic nickel acetate–5,5′‐dmbpy–KClO₄ system. Within the complex cation, the Ni^II atom is hexacoordinated by two chelating 5,5′‐dmbpy ligands and one chelating ac ligand. The mean Ni—N and Ni—O bond lengths are 2.0628 (17) and 2.1341 (15) Å, respectively. The water solvent molecule is disordered over two partially occupied positions and links two complex cations and two perchlorate anions into hydrogen‐bonded centrosymmetric dimers, which are further connected by π–π interactions. The magnetic properties of (1) at low temperatures are governed by the action of single‐ion anisotropy, D, which arises from the reduced local symmetry of the cis‐NiO₂N₄ chromophore. The fitting of the variable‐temperature magnetic data (2–300 K) gives g_iso = 2.134 and D/hc = 3.13 cm⁻¹.     

Nature of closed‐ and open‐shell interactions between noble metals and rare gas atoms

Interactions between noble metals and rare gases have become an interesting topic over the last few years. In this work, a computational study of the open‐shell (d¹⁰s¹) and closed‐shell (d¹⁰s and d¹⁰s²) noble metals (M = Cu, Ag, and Au) with three heaviest rare gas atoms (Rg = Kr, Xe, and Rn) has been performed. Potential energy curves based on ab initio [MP2, MP4, QCISD, and CCSD(T)] and DFT functionals (M06‐2X and CAM‐B3LYP) were obtained for ionic and neutral AuXe complexes. Dissociation energies indicate that neutral metals have the lowest and cationic metals have the highest affinities for interaction with rare gas atoms. For the same metals, there is a continuous increase in dissociation energies (D_e) from Kr to Rn. The nature of bonding and the trend of D_e and equilibrium bond lengths (R_e) have been interpreted by means of quantum theory of atoms in molecules, natural bond orbital, and energy decomposition analysis. © 2013 Wiley Periodicals, Inc.     

EHF theory of chemical reactions V. Nature of manganese–oxygen bonds by hybrid density functional theory (DFT) and coupled‐cluster (CC) methods

Hybrid density functional theory (DFT) and post‐Hartree–Fock methods are compared by depicting potential energy curves of the O–O dissociation of hydroperoxide and the M–O dissociation of transition‐metal oxides. The former approach includes BLYP, B2LYP, B3LYP, and more general hybrid DFT methods, while the unrestricted Hartree–Fock (UHF) coulpled‐cluster (UCC) SD(T) method is considered as the latter approach. The hybrid DFT methods can reproduce the potential curve of the O–O dissociation process and the dissociation energy of HOOH by UCCSD(T). The methods are also useful for depicting potential curves of copper oxide (CuO) and manganese oxide (MnO), and reproduce the experimental M–O binding energies. The nature of Mn–O bonds in the naked Mn–O, Mn–O porphyrine system and model complexes, XH₃Mn(IV)O₂Mn(IV)H₃Y (X,Y=O,H), are examined in relation to the possible mechanisms of oxygenation reactions. It is found that the radical character of Mn–O bonds increases with the increase of the oxidation number of the Mn ion in these systems. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

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.