Reactions of o-benzyne with propargyl and benzyl radicals: potential sources of polycyclic aromatic hydrocarbons in combustion

The kinetics and mechanisms of the reactions of o-benzyne with propargyl and benzyl radicals have been investigated computationally. The possible reaction pathways have been explored by quantum chemical calculations at the M06-2X/6-311+G(3df,2p)//B3LYP/6-311G(d,p) level and the mechanisms have been investigated by the Rice-Ramsperger-Kassel-Marcus theory/master-equation calculations. It was found that the o-benzyne associates with the propargyl and benzyl radicals without pronounced barriers and the activated adducts easily isomerize to five-membered ring species. Indenyl radical and fluorene + H were predicted to be dominantly produced by the reactions of o-benzyne with propargyl and benzyl radicals, respectively, with the rate constants close to the high-pressure limits at temperatures below 2000 K. The related reactions on the two potential energy surfaces, namely, the reaction between fulvenallenyl radical and acetylene and the decomposition reactions of indenyl and α-phenylbenzyl radicals were also investigated. The high reactivity of o-benzyne toward the resonance stabilized radicals suggested a potential role of o-benzyne as a precursor of polycyclic aromatic hydrocarbons in combustion.     

ENDOR AND ESR STUDIES OF RADICALS FORMED IN FRIEDEL-CRAFTS REACTIONS

The present study deals with the study of some Friedel-Crafts alkylating system and the accurate measurement of hyperfine coupling constants by ESR and ENDOR techniques. The results indicate that the observed ESR spectra are due to polycyclic aromatic radical cations formed from their parent hydrocarbons.It is suggested that benzyl halides produced in the Friedel-Crafts alkylating reaction undergo Scholl condensation reaction to give polycyclic aromatic hydrocarbons, which were converted into the corresponding polycyclic aromatic radical cations in the presence of aluminum chloride.     

戊烯自由基阳离子的密度泛函理论研究

使用密度泛函理论B3LYP方法和6-31G(d,p)、6-31＋G(d,p)、6-311G(d,p)及6-311＋G(d,p)基组，分别对2-C5H10＋和1-C5H10＋的各种构象进行了几何构型优化，并用B3LYP/6-311G(d,p)进行了频率分析计算.计算预言1-C5H10＋具有非平面构型，与以往报导的从头算计算结论相反.在两个自由基阳离子的各种构象的B3LYP几何构型上，进行了B3LYP和UMP2(full)方法的超精细偶合常数计算，得到了比以往更好的结果.     

Gas-phase ion/molecule reactions between dimethoxyphosphonium ions and aromatic hydrocarbons

Ion/molecule reactions between O=P(OCH(3))(2)(+) phosphonium ions and six aromatic hydrocarbons (benzene, toluene, 1,2,4-trimethylbenzene, naphthalene, acenaphthylene and fluorene) were performed in a quadrupole ion trap mass spectrometer. The O=P(OCH(3))(2)(+) phosphonium ions, formed by electron impact from neutral trimethyl phosphite, were found to react with aromatic hydrocarbons (ArHs) to give (i) an adduct [ArH, O=P(OCH(3))(2)](+) and (ii) for ArHs which have an ionization energy below or equal to 8.14 eV, a radical cation ArH(+ *) by charge transfer reaction. Collision-induced dissociation experiments, which produce fragment ions other than the O=P(OCH(3))(2)(+) ions, indicate that the adduct ions are covalent species. Isotope-labeled ArHs were used to elucidate fragmentation mechanisms. The charge transfer reactions were investigated using density functional theory at the B3LYP/6-311 + G(3df,2p)//B3LYP/6-31G(d,p) level of theory. The potential energy surface obtained from B3LYP/6-31G(d,p) calculations for the reaction between O=P(OCH(3))(2)(+) and benzene is described.     

Ab initio calculation of inner-sphere reorganization energies of arenediazonium ion couples

The geometries of a series of substituted arenediazonium cations (p-NO2, p-CN, p-Cl, p-F, p-H, m-CH3, p-CH3, p-OH, p-OCH3, p-NH2) and the corresponding diazenyl radicals were optimized at the HF/6-31G, MP2/6-31G, B3LYP/6-31G, B3LYP/TZP, B3PW91/TZP, and CASSCF/6-31G levels of theory. Inner-sphere reorganization energies for the single electron-transfer reaction between the species were computed from the optimized geometries according to the NCG method and compared to experimental values determined by Doyle et al. All levels of theory predicted a CNN bond angle of 180 degrees in the cation. A bent neutral diazenyl radical was predicted at all levels of theory excepting B3LYP/TZP and B3PW91/TZP for the p-Cl-substituted compound. Inner-sphere reorganization energies determined at the HF, MP2, and CASSCF levels of theory correlated poorly with both experimental results and calculated geometries. Density functional methods correlated best with the experimental values, with B3LYP/6-31G yielding the most promising results, although the ROHF/6-31G survey also showed some promise. B3LYP/6-31G calculations correctly predicted the order of the inner-sphere reorganization energies for the series, excluding the halogen-substituted compounds, with values ranging from 42.8 kcal x mol(-1) for the p-NO2-substituted species to 55.1 kcal x mol(-1) for NH2. The magnitudes of these energies were lower than the experimental by a factor of 2. For the specific cases examined, the closed-shell cation geometries showed the expected geometry about the CNN bond, with variations in the CN and NN bond lengths correlating with the electron-donating/withdrawing capacity of the substituent. As predicted by Doyle et al., a large geometry change was observed upon reduction. The neutral diazenyl radicals showed a nominal CNN bond angle of 120 degrees and variations in the CN and NN bond lengths also correlated with the electron-donating/withdrawing capacity of the substituent. Changes in theta(CNN) and r(CN) both correlated well with calculated lambda(inner). The key parameters influencing inner-sphere reorganization energy were the CN and NN bond lengths and the CNN bond angle. This influence is explained qualitatively via resonance models produced from NRT analysis and is related to the amount of CN double bond character. Based on these observations, B3LYP/6-31G calculations are clearly the most amenable for calculating inner-sphere reorganization energies for the single electron-transfer reaction between cation/neutral arenediazonium ion couples.     

已烯自由基阳离子的密度泛函理论研究

使用密度泛函理论B3LYP方法和6-31G(d,p),6-31+G(d,p),6-311G(d,p)及6- 311+G(d,p)基组,分别对1-C_6H_(12)~+,2-C_6H_(12)~+和3-C_6H_(12)~+的各种构 象进行了几何构型优化,并在B3LYP/6-311G(d,p)水平上进行了频率分析计算,在 各优化构型上,使用B3LYP和MP2(full)方法进行了超精细结构的计算。计算的3- C_6H_(12)~+的超精细偶合常数比以往的计算结果更好;1-C_6H_(12)~+和2-C_6H_ (12)~+的超精细偶合常数目前尚无实验数据报道,本计算预言了它们的超精细偶合 常数和最稳定构型。     

Disclosure of Ground-State Zimmerman-Möbius Aromaticity in the Radical Anion of [6]Helicene and Evidence for 4π Periodic Aromatic Ring Currents in a Molecular “Metallic” Möbius Strip

In 1966, Zimmerman proposed a type of Möbius aromaticity that involves through-space electron delocalization; it has since been widely applied to explain reactivity in pericyclic reactions, but is considered to be limited to transition-state structures. Although the easily accessible hexahelicene radical anion has been known for more than half a century, it was overlooked that it exhibits a ground-state minimum and robust Zimmerman-Möbius aromaticity in its central noose-like opening, becoming, hence, the oldest existing Möbius aromatic system and with smallest Möbius cycle known. Despite its overall aromatic stabilization energy of 13.6 kcal mol⁻¹ (at B3LYP/6-311+G**), the radical also features a strong, globally induced paramagnetic ring current along its outer edge. Exclusive global paramagnetic currents can also be found in other fully delocalized radical anions of 4N+2 π-electron aromatic polycyclic benzenoid hydrocarbons (PAH), thus questioning the established magnetic criterion of antiaromaticity. As an example of a PAH with nontrivial topology, we studied a novel Möbius[16]cyclacene that has a non-orientable surface manifold and a stable closed-shell singlet ground state at several density functional theory levels. Its metallic monoanion radical (0.0095 eV band gap at HSE06/6-31G* level) is also wave-function stable and displays an unusual 4π-periodic, magnetically induced ring current (reminiscent of the transformation behaviour of spinors under spatial rotation), thus indicating the existence of a new, Hückel-rule-evading type of aromaticity.     

Possible molecular hydrogen formation mediated by the radical cations of anthracene and pyrene

Hydrogen molecules cannot be formed readily by the association of gaseous hydrogen atoms. Possible H(2) formation mediated by the radical cations of typical polycyclic aromatic hydrocarbons (PAHs), anthracene and pyrene, was studied at the B3LYP/6-31G** level of theory. We presumed that H(2) is formed by way of two elementary reactions: the addition of an H atom to a PAH molecular cation, and the H abstraction from the resulting monohydro-PAH cation (i.e., arenium ion) by a second H atom to yield H(2). The first reaction takes place without any activation energy. The second reaction is also predicted to proceed along almost barrierless pathways, although it is far from being a typical ion-molecule reaction. There is a possibility that these reactions might constitute one of the mechanisms for H(2) formation in extremely cold interstellar space. Deuterium enrichment in PAH cations is possibly accompanied by such H(2) formation because deuteration lowers the energies of polyatomic PAH cations appreciably.     

HO3、HO+3、HO-3光谱和结构的MP2和DFT理论研究

用B3LYP和MP2量子化学理论,详细地计算了HO3、HO+3、HO-3可能的平衡构型、能量和光谱参数,分析了它们可能的解离通道和稳定性.还分析了HO3的电离势(IE)及电子亲和势(Ea),讨论了HO3、HO3-自由基的光电子能谱的特征,及用光电子能谱去探测HO3的可能性.     

氟代吡啶阳离子的理论研究

Cations of fluorinated pyridines（pentafluoropyridine,2,6-difluoropyridine,and 2-fluoropyridine）have been studied by using density functional B3LYP method in conjunction with 6-31G（d,p）,6-311G（d,p）,6-31+G（d,p）,and 6-311+G（d,p）basis sets. B3LYP geometry optimization and frequency analysis calculations indicate that the pentafluoropyridine cation,2,6-difluoropyridine cation,and 2-fluoropyridine cation have C2v,C2v,and Cs structures in the 2A2,2A2,and 2A" ground states,respectively. The calculated geometries of the cations and the parent molecules were compared. The natural population analysis calculations at the B3LYP level with different basis sets were performed on the three cations and the three parent molecules. The isotropic hyperfine coupling constants in the three cations（radicals）were calculated. The vertical and adiabatic ionization potential（VIP and AIP）values of the pentafluoropyridine,2,6-difluoropyridine,and 2-fluoropyridine molecules were calculated by using the B3LYP method,and the calculated VIP values are in excellent agreement with experiment.     

Functionalized nanodiamonds part I. An experimental assessment of diamantane and computational predictions for higher diamondoids

The structures, strain energies, and enthalpies of formation of diamantane 1, triamantane 2, isomeric tetramantanes 3-5, T(d)-pentamantane 6, and D(3d)-hexamantane 7, and the structures of their respective radicals, cations, as well as radical cations, were computed at the B3LYP/6-31G* level of theory. For the most symmetrical hydrocarbons, the relative strain (per carbon atom) decreases from the lower to the higher diamondoids. The relative stabilities of isomeric diamondoidyl radicals vary only within small limits, while the stabilities of the diamondoidyl cations increase with cage size and depend strongly on the geometric position of the charge. Positive charge located close to the geometrical center of the molecule is stabilized by 2-5 kcal mol(-1). In contrast, diamondoid radical cations preferentially form highly delocalized structures with elongated peripheral C-H bonds. The effective spin/charge delocalization lowers the ionization potentials of diamondoids significantly (down to 176.9 kcal mol(-1) for 7). The reactivity of 1 was extensively studied experimentally. Whereas reactions with carbon-centered radicals (Hal)(3)C(*) (Hal=halogen) lead to mixtures of all possible tertiary and secondary halodiamantanes, uncharged electrophiles (dimethyldioxirane, m-chloroperbenzoic acid, and CrO(2)Cl(2)) give much higher tertiary versus secondary selectivities. Medial bridgehead substitution dominates in the reactions with strong electrophiles (Br(2), 100 % HNO(3)), whereas with strong single-electron transfer (SET) acceptors (photoexcited 1,2,4,5-tetracyanobenzene) apical C(4)-H bridgehead substitution is preferred. For diamondoids that form well-defined radical cations (such as 1 and 4-7), exceptionally high selectivities are expected upon oxidation with outer-sphere SET reagents.     

Reactivity of superoxide radical anion with cyclic nitrones: role of intramolecular H-bond and electrostatic effects

Limitations exist among the commonly used cyclic nitrone spin traps for biological free radical detection using electron paramagnetic resonance (EPR) spectroscopy. The design of new spin traps for biological free radical detection and identification using EPR spectroscopy has been a major challenge due to the lack of systematic and rational approaches to their design. In this work, density functional theory (DFT) calculations and stopped-flow kinetics were employed to predict the reactivity of functionalized spin traps with superoxide radical anion (O2*-). Functional groups provide versatility and can potentially improve spin-trap reactivity, adduct stability, and target specificity. The effect of functional group substitution at the C-5 position of pyrroline N-oxides on spin-trap reactivity toward O2*- was computationally rationalized at the PCM/B3LYP/6-31+G(d,p)//B3LYP/6-31G(d) and PCM/mPW1K/6-31+G(d,p) levels of theory. Calculated free energies and rate constants for the reactivity of O2*- with model nitrones were found to correlate with the experimentally obtained rate constants using stopped-flow and EPR spectroscopic methods. New insights into the nucleophilic nature of O2*- addition to nitrones as well as the role of intramolecular hydrogen bonding of O2*- in facilitating this reaction are discussed. This study shows that using an N-monoalkylsubstituted amide or an ester as attached groups on the nitrone can be ideal in molecular tethering for improved spin-trapping properties and could pave the way for improved in vivo radical detection at the site of superoxide formation.     

Magnetic properties of single crystal alpha-benzoin oxime: An EPR study

The electron paramagnetic resonance (EPR) spectra of gamma irradiated single crystals of alpha-benzoinoxime (ABO) have been examined between 120 and 440 K. Considering the dependence on temperature and the orientation of the spectra of single crystals in the magnetic field, we identified two different radicals formed in irradiated ABO single crystals. To theoretically determine the types of radicals, the most stable structure of ABO was obtained by molecular mechanic and B3LYP/6-31G(d,p) calculations. Four possible radicals were modeled and EPR parameters were calculated for the modeled radicals using the B3LYP method and the TZVP basis set. Calculated values of two modeled radicals were in strong agreement with experimental EPR parameters determined from the spectra. Additional simulated spectra of the modeled radicals, where calculated hyperfine coupling constants were used as starting points for simulations, were well matched with experimental spectra.     

A DFT‐Based Investigation of Hydrogen Abstraction Reactions from Methylated Polycyclic Aromatic Hydrocarbons

The growth of polycyclic aromatic hydrocarbons (PAHs) is in many areas of combustion and pyrolysis of hydrocarbons an inconvenient side effect that warrants an extensive investigation of the underlying reaction mechanism, which is known to be a cascade of radical reactions. Herein, the focus lies on one of the key reaction classes within the coke formation process: hydrogen abstraction reactions induced by a methyl radical from methylated benzenoid species. It has been shown previously that hydrogen abstractions determine the global PAH formation rate. In particular, the influence of the polyaromatic environment on the thermodynamic and kinetic properties is the subject of a thorough exploration. Reaction enthalpies at 298 K, reaction barriers at 0 K, rate constants, and kinetic parameters (within the temperature interval 700–1100 K) are calculated by using B3LYP/6‐31+G(d,p) geometries and BMK/6‐311+G(3df,2p) single‐point energies. This level of theory has been validated with available experimental data for the abstraction at toluene. The enhanced stability of the product benzylic radicals and its influence on the reaction enthalpies is highlighted. Corrections for tunneling effects and hindered (or free) rotations of the methyl group are taken into account. The largest spreading in thermochemical and kinetic data is observed in the series of linear acenes, and a normal reactivity–enthalpy relationship is obtained. The abstraction of a methyl hydrogen atom at one of the center rings of large methylated acenes is largely preferred. Geometrical and electronic aspects lie at the basis of this striking feature. Comparison with hydrogen abstractions leading to arylic radicals is also made.     

Ab initio inter-radical potential for p-phenylenediamine radical cation dimer

Ab initio molecular orbital calculations were carried out to investigate the inter-radical interaction of the paired p-phenylenediamine radical cations in the singlet state. After initial optimization of the dimer in the parallel sandwich (D2h) and parallel displaced (Cs) configurations at the B3LYP/6-31G* theoretical level, the MP2/6-31G* and B3LYP/6-31G* single energies of the dimer were calculated as a function of the inter-radical distance R. The depths of the potential minima near R = 3.2 A were estimated to be in the order of the hydrogen bonding energy, assuming that the electrostatic contribution between the cations is canceled out by the attractive contributions due to the counter anions on the aspect of a simple electrostatic model. This can be related to the indications of the cation dimer formation in solution in the presence of counter anions at a low temperature reported previously in the literature by resonance Raman and electronic absorption spectra. The inter-radical (Raman active) frequencies of the dimer were calculated, one of which corresponds to the reported value at 161 cm(-1) observed in the resonance Raman spectrum in ethanol at 200 K by Yokoyama and Maeda (Chem. Phys. Lett. 48 (1977) 59).     

(CH2)2N和(CH3)2NH+的密度泛函理论计算

用密度泛函理论UB3LYP／6-31G（d,p）方法研究了二甲胺自由基（CH3）2N及其质子化离子（CH3）2NH 的构象和超精细结构．在由构象研究得到的两种自由基的最稳定结构上,用密度泛函的UB3LYP和UB3PW91方法及从头计算的UHF、UMP2（FULL）方法计算了α-质子、β-质子和N核上的超精细偶合常数A（Hα）、A（Hβ）和A（N）结果表明：两种自由基中甲基内旋转的位垒均很低,分别为0．46kJ·mol-1（（CH3）2NH ）和2．6kJ·mol-1（（CH3）2N）．UB3LYP／6-31G（d,p）和UB3PW91／6-31G（d,p）计算的A（Hα）、A（Hβ）和A（N）与ESR实验结果符合得很好,UMP2／6-31G（d,p）方法的计算值与实验值符合得也较好．     

Superoxide radical anion adduct of 5,5-dimethyl-1-pyrroline n-oxide (DMPO). 1. The thermodynamics of formation and its acidity

The nitrone 5,5-dimethyl-1-pyrroline N-oxide (DMPO) has been the most widely used spin trap for the detection of transient free radicals in chemical, biological, and biomedical research using electron paramagnetic resonance (EPR) spectroscopy. A density functional theory (DFT) approach was used to predict the thermodynamics of formation of the superoxide anion/hydroperoxyl radical (O2*-/*O2H) adduct of DMPO as well as its pK(a) in aqueous systems. At the B3LYP/6-31+G(d,p)//B3LYP/6-31G(d) level, we predicted (in the gas phase and with a polarizable continuum model (PCM) for water) three conformational minima for both the DMPO-O2- and DMPO-O2H adducts. Using DFT and the PCM solvation method, the pK(a) of DMPO-O2H was predicted to be 14.9 +/- 0.5. On the basis of free energy considerations, the formation of DMPO-O2H at neutral pH proceeds via initial addition of O2*- to DMPO to form the DMPO-O2- adduct and then subsequent protonation by water (or other acidic sources) to form DMPO-O2H. Under acidic conditions, the addition of *O2H to DMPO is predicted to be more exoergic than the addition of O2*- and is consistent with available experimental kinetic data.     

Thermochemical properties and bond dissociation energies of C3-C5 cycloalkyl hydroperoxides and peroxy radicals: cycloalkyl radical + (3)O2 reaction thermochemistry

Cyclic aliphatic hydrocarbons are major components in modern fuels; they can be present in the reactants, and they can be formed during the gas-phase oxidation processes. In combustion and thermal oxidation processes, these cyclics will form radicals that react with (3)O(2) to form peroxy radicals. In this study, density functional theory and higher level ab initio calculations are used to calculate thermochemical properties and bond dissociation energies of 3-5-membered cycloalkanes, corresponding hydroperoxides, hydroperoxycycloalkyl radicals, and cycloalkyl radicals that occur in these reaction systems. Geometries, vibration frequencies, and thermochemical properties, ΔH(f 298)°, are calculated with the B3LYP/6-31 g(d,p), B3LYP/6-31 g(2d,2p), composite CBS-QB3, and G3MP2B3 methods. Standard enthalpies of formation at 298 K are evaluated using isodesmic reaction schemes with several work reactions for each species. Group additivity contributions are developed, and application of group additivity with comparison to calculated values is illustrated. Entropy and heat capacities, S°(T) and C(p)°(T) (5 K ≤ T ≤ 5000), are determined using geometric parameters and frequencies from the B3LYP/6-31 g(d,p) calculations.     

Reaction of ketenyl radical with acetylene: a promising route for cyclopropenyl radical

The reaction of the ketenyl radical (HCCO) with acetylene (C(2)H(2)) is very relevant to the oxygen-acetylene flames and fuel-rich combustion process for nitrogen-containing compounds. Unfortunately, except for several rate constant measurements, the mechanism is completely unknown for this reaction. In this paper, detailed theoretical investigations are performed for the HCCO + C(2)H(2) reaction at the G3B3 level using the B3LYP/6-31G(d), B3LYP/6-311++G(d,p), and QCISD/6-31G(d) geometries. The exclusive fragmentation channel is the formation of the cyclopropenyl radical (c-C(3)H(3)) and carbon monoxide (CO) via the chainlike OCCHCHCH and three-membered ring OC-cCHCHCH intermediates. Thus, the mass spectroscopic peak of C(3)H(3)(+) in a previous experiment can be explained. The calculated overall reaction barrier is 4.4, 4.4, and 5.3 kcal/mol at the G3B3//B3LYP/6-31G(d), G3B3//B3LYP/6-311++G(d,p), and G3B3//QCISD/6-31G(d) levels, respectively. The title reaction may provide an effective route for generating the long-sought cyclopropenyl radical in the laboratory, which has been the long-standing subject of numerous theoretical studies as the simplest cyclic conjugate radical, and its bulky derivatives were already known. Future experimental investigations for the HCCO + C(2)H(2) reaction are greatly desired to test the predicted fragmentation channel. The implication of the present study in combustion and interstellar processes is discussed.     

Spin trapping by 5-carbamoyl-5-methyl-1-pyrroline N-oxide (AMPO): theoretical and experimental studies

The nitrone 5-carbamoyl-5-methyl-1-pyrroline N-oxide (AMPO) was synthesized and characterized. Spin trapping of various radicals by AMPO was demonstrated for the first time by electron paramagnetic resonance (EPR) spectroscopy. The resulting spin adducts for each of these radicals gave unique spectral profiles. The hyperfine splitting constants for the superoxide adduct are as follows: isomer I (80%), a(nitronyl)(-)(N) = 13.0 G and a(beta)(-)(H) = 10.8 G; isomer II (20%), a(nitronyl)(-)(N) = 13.1 G, a(beta)(-)(H) = 12.5 G, and a(gamma)(-)(H) = 1.75 G. The half-life of the AMPO-O(2)H was about 8 min, similar to that observed for EMPO but significantly shorter than that of the DEPMPO-O(2)H with t(1/2) approximately 16 min. However, the spectral profile of AMPO-O(2)H at high S/N ratio is distinguishable from the spectrum of the (*)OH adduct. Theoretical analyses using density functional theory calculations at the B3LYP/6-31+G//B3LYP/6-31G level were performed on AMPO and its corresponding superoxide adduct. Calculations predicted the presence of intramolecular H-bonding in both AMPO and its superoxide adduct. The H-bonding interaction was further confirmed by an X-ray structure of AMPO, and of the novel and analogous amido nitrone 2-amino-5-carbamoyl-5-methyl-1-pyrroline N-oxide (NH(2)-AMPO). The thermodynamic quantities for superoxide radical trapping by various nitrones have been found to predict favorable formation of certain isomers. The measured partition coefficient in an n-octanol/buffer system of AMPO was similar to those of DMPO and DEPMPO. This study demonstrates the suitability of the AMPO nitrone for use as a spin trap to study radical production in aqueous systems.