Benchmarking of DFT functionals for the kinetics and mechanisms of atmospheric addition reactions of OH radicals with phenyl and substituted phenyl‐based organic pollutants

^•OH addition reactions play a pivotal role in the atmospheric transformation of a number of phenyl and substituted phenyl‐based persistent and toxic organic pollutants. Here, we screened appropriate DFT functionals to predict reaction mechanisms and rate constants (k_OH) of the ^•OH additions by taking benzene and substituted benzenes (C₆H₅F, C₆H₅Cl, C₆H₅Br, C₆H₅CH₃, C₆H₅OH) as model compounds. By comparing the k_OH values calculated with DFT methods to experimental values, we found that the ωB97 functional is the best among the 18 functionals considered (using the basis sets 6‐31 + G(d,p) for optimizations and 6‐311++G(3df,2pd) for single point energy calculations) in the temperature range of 230‐330 K. In addition, we found that some other functionals performed well in specific conditions, e.g., BMKD3 is good for benzene, halogenated benzenes and C₆H₅CH₃, and CAM‐B3LYP is good for the reaction of C₆H₅OH at room temperature. Based on the diversity of the electronic structures of the selected model compounds and the frequent occurrence of certain substituents ( CH₃,  OH,  F,  Cl, and  Br) in the target compounds, the functionals recommended here can be used for future study of the reaction mechanisms and k_OH values for ^•OH addition to phenyl and substituted phenyl‐based persistent and toxic organic pollutants.     

Charge inversion as a structural probe for C6H 5 + and C6H 6 +· cations

Charge reversal (⁺CR⁻) of cations to anions can be used to structurally differentiate isomeric C₆H₅⁺ and C₆H₆^+· hydrocarbon ions by means of tandem mass spectrometry. In view of the manifold of possible isomers, only a few prototype precursors are examined. Thus, charge inversion demonstrates that electron ionization of 2,4-hexadiyne yields an intact molecular ion, whereas the charge inversion spectra of C₆H₆^+· obtained from benzene, 1,5-hexadiyne, and fulvene are identical within experimental error. Similarly, the ⁺CR⁻ spectrum of the C₆H₅⁺ cation generated by dissociative ionization of 2,4-hexadiyne is significantly different from the ⁺CR⁻ spectrum of C₆H₅⁺ obtained from iodobenzene, suggesting the formation of a 2,4-hexadiynyl cation from the former precursor. Although charge inversion of cations to anions has a low efficiency and requires large precursor ion fluxes, the particular value of this method is that the spectra may not just differ in fragment ion intensities, but these differences can directly be related to the underlying ion structures.     

Energy and substituent effects on the mass spectra of substituted benzophenones

For two competing decompositions of the same molecular ion to give products [A₁⁺] and [A₂⁺], the ratio [A₁⁺]/[A₂⁺], is equal to the ratio of rate constants for the formation of the stable ions. Thr ratios, [Y C₇H₄O⁺]/[C₇H₅O⁺], were determined for several benzophenones for electron energies from 15 to 70 eV. Plots of log [Y C₇H₄O⁺]/[C₇H₅O⁺] vs.[ω⁺] gave good straight lines at all energies. Similar correlations have been reported for log [Y C₇H₆⁺]/[C₇H₇⁺] from substituted diphenyl ethanes and are also true for substituted acetophenones, log [YøCO⁺]/[CH₃CO⁺]. A few charge exchange data were obtained which showed the same general trend as the electron-impact data and emphasize the contribution of low energy ions in the 70 eV mass spectra. Relatively poor correlations were obtained for the [Y C₆H₄⁺] and [C₆H₅⁺] ions that are formed by both one-step and two-step decompositions.     

A structural investigation of [C6H6]2+ and [C6H6]+˙ ions involved in charge exchange reactions

Mass-analysed ion kinetic energy spectra for collisional activation (CA) of [C₆H₆]⁺˙ formed via electron capture by [C₆H₆]²⁺ ions in collision with neutral benzene molecules have been compared for the C₆H₆ isomers benzene, 1,5-hexadiyne and 2,4-hexadiyne. Comparisons of fragment abundance and total CA fragment yields were also made for [C₆H₆]⁺˙ ions generated by electron ionization (EI). CA conditions of ion velocity and collision gas pressure were identical in these comparisons. In general the fragment abundance patterns for the ions formed by charge exchange were very similar to those for singly charged benzene ions generated by EI. However, significant variations in CA fragment yield (the ratio of the total CA fragment ion abundance to the abundance of the incident unfragmented ions) were observed. It is not clear from the results whether these variations reflect structurally different ions or ions of different internal energies. The CA spectra of [C₆H₆]⁺˙ ions derived from charge exchange reactions between the benzene dication and the target gases He, Ne, Ar, Kr and Xe have also been recorded and, once again, very similar fragment abundance patterns were observed along with large variations in total CA fragment yields. Charge exchange efficiency measurements are reported for reactions between the benzene dication and the targets He, Ne, Ar, Kr, Xe and C₆H₆ (benzene) and also for the doubly charged ions derived from the linear C₆H₆ isomers. In the latter case Xe and benzene targets were used. The energetics and efficiency measurements for the former reactions suggest that for targets such as He and Ne the processes probably involve excited states of the doubly charged ions. The efficiencies measured for the latter reactions were distinctly different for the three C₆H₆ isomers and may indicate a strong dependence of charge exchange cross-section on doubly charged ion structure.     

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.     

Experimental and theoretical study on the photoionization of styrene

This study employed a vacuum ultraviolet synchrotron radiation source and reflectron time-of-flight mass spectrometry (TOF-MS) to investigate the photoionization and dissociation of styrene. By analyzing the photoionization mass spectrum and efficiency curve alongside G3B3 theoretical calculations, we determined the ionization energy of the molecular ion, appearance energy of fragment ions, and relevant dissociation pathways. The major ion peaks observed in the photoionization mass spectra of styrene correspond to C₈H₈⁺, C₈H₇⁺ and C₆H₆⁺. The ionization energy of styrene is measured as 8.46 ± 0.03 eV, whereas the appearance energies of C₈H₇⁺ and C₆H₆⁺ are found to be 12.42 ± 0.03 and 12.22 ± 0.03 eV, respectively, in agreement with theoretical values. The main channel for the photodissociation of styrene molecular ions is the formation of benzene ions, whereas the dissociation channel that loses hydrogen atoms is the secondary channel. Based on the experimental results and empirical formulas, the required dissociation energies (E_d) of C₈H₇⁺, C₈H₆⁺ and C₆H₆⁺ are calculated to be (3.96 ± 0.06), (4.00 ± 0.06) and (3.76 ± 0.06) eV, respectively. Combined with related thermochemical parameters, the standard enthalpies of formations of C₈H₈⁺, C₈H₇⁺, C₈H₆⁺ and C₆H₆⁺ are determined to be 964.2, 1346.3, 1350.2 and 1327.0 kJ/mol, respectively. Based on the theoretical study, the kinetic factors controlling the styrene dissociation reaction process are determined by using the Rice–Ramsperger–Kassel–Marcus (RRKM) theory. This provides a reference for further research on the atmospheric photooxidation reaction mechanism of styrene in atmospheric and interstellar environments.     

Charge separation reactions of doubly charged benzene ions

Photoion-photoion coincidence spectra of benzene and benzene-d₆ photoionized by He(II) light and synchrotron radiation show the existence of six major and eight minor charge-separation reactions of the [C₆H₆]²⁺ ion. Three main groups of ion pairs are related to [C₃H₃]⁺ + [C₃H₃]⁺, [C₂H₃]⁺ + [C₄H₃]⁺ and [CH₃]⁺ + [C₅H₃]⁺, with appearance energies of 32.2 ± 0.5 eV, 31.3 ± 0.5 eV and 28.4 ± 0.3 eV. The kinetic energy release is the same for all pairs within a group, irrespective of hydrogen number, but differs from group to group. Results are interpreted in terms of fast, direct charge separation of [C₆H₆]²⁺, and subsequent hydrogen loss by the singly charged fragments.     

The study of C1?C3 monosubstituted alkyl benzenes by the inverse convolution of first differential ionization efficiency curves

The electron impact ionization efficiency curves for the parent ions and the [C₇H₇]⁺ fragment ion formed from monosubstituted alkyl benzenes (R?CH₃? n-C₃H₇) have been studied by applying the inverse convolution technique of Vogt and Pascual to the first derivative ionization efficiency curves of the ions. Ionization and appearance energies measured for the ions at threshold are in good agreement with recently published photoionization values. Structures in the ionization efficiency curves (higher energy processes) are also reported for about 4 e V above threshold. The heats of formation calculated for [C₇H₇]⁺ fragment ions obtained from toluene and ethyl benzene at threshold are equal to 864 and 865 kJ mol^?1 respectively, and are consistent with the tropylium structure. However, for the [C₇H₇]⁺ fragment ion obtained from n-propyl benzene at threshold the calculated heat of formation is equal to 923 kJ mol^?1 and probably corresponds to a benzyl structure.     

Dissociative photoionization of isoprene: experiments and calculations

Vacuum ultraviolet (VUV) dissociative photoionization of isoprene in the energy region 8.5–18 eV was investigated with photoionization mass spectroscopy (PIMS) using synchrotron radiation (SR). The ionization energy (IE) of isoprene as well as the appearance energies (AEs) of its fragment ions C₅H₇⁺, C₅H₅⁺, C₄H₅⁺, C₃H₆⁺, C₃H₅⁺, C₃H₄⁺, C₃H₃⁺ and C₂H₃⁺ were determined with photoionization efficiency (PIE) curves. The dissociation energies of some possible dissociation channels to produce those fragment ions were also determined experimentally. The total energies of C₅H₈ and its main fragments were calculated using the Gaussian 03 program and the Gaussian‐2 method. The IE of C₅H₈, the AEs for its fragment ions, and the dissociation energies to produce them were predicted using the high‐accuracy energy model. According to our results, the experimental dissociation energies were in reasonable agreement with the calculated values of the proposed photodissociation channels of C₅H₈. Copyright © 2009 John Wiley & Sons, Ltd.     

Laser-induced dissociation of phenetole ions prepared by resonance-enhanced two-photon,two-color ionization

Laser-induced dissociation of phenetole ions has been carried out. The ions have been prepared in a well-defined excited state by resonance-enhanced two-photon, two-color ionization (R2PI/2C). Appearance energies of 2.431 and 2.82 eV have been determined for the first dissociation pathways leading to C₆H₆O⁺ and C₅H⁶⁺ ionic fragment. The dissociation spectrum as well as detailed dissociation rate constants versus internal energy have been obtained. The dissociation rates are compared with those expected from RRKM theory. A complex reaction scheme has been assumed in order to explain the observed discrepancy.     

Vacuum ultraviolet photoionization mass spectrometric study of cyclohexene

In this work, photoionization and dissociation of cyclohexene have been studied by means of coupling a reflectron time‐of‐flight mass spectrometer with the tunable vacuum ultraviolet (VUV) synchrotron radiation. The adiabatic ionization energy of cyclohexene as well as the appearance energies of its fragment ions C₆H₉⁺, C₆H₇⁺, C₅H₇⁺, C₅H₅⁺, C₄H₆⁺, C₄H₅⁺, C₃H₅⁺ and C₃H₃⁺ were derived from the onset of the photoionization efficiency (PIE) curves. The optimized structures for the transition states and intermediates on the ground state potential energy surfaces related to photodissociation of cyclohexene were characterized at the ωB97X‐D/6‐31+g(d,p) level. The coupled cluster method, CCSD(T)/cc‐pVTZ, was employed to calculate the corresponding energies with the zero‐point energy corrections by the ωB97X‐D/6‐31+g(d,p) approach. Combining experimental and theoretical results, possible formation pathways of the fragment ions were proposed and discussed in detail. The retro‐Cope rearrangement was found to play a crucial role in the formation of C₄H₆⁺, C₄H₅⁺ and C₃H₅⁺. Intramolecular hydrogen migrations were observed as dominant processes in most of the fragmentation pathways of cyclohexene. The present research provides a clear picture of the photoionization and dissociation processes of cyclohexene in the 8‐ to 15.5‐eV photon energy region. Copyright © 2016 John Wiley & Sons, Ltd.     

Reactions of Au+ and Au− with benzene and fluorine-substituted benzenes

Reactions of gold anions and cations generated by laser desorption/ionization were studied in the FTICR spectrometer. Au⁻ associated with C₆F₆ to give the novel Au⁻(C₆F₆) complex, whose binding energy was estimated to be 24 ± 4 kcal mol⁻¹ from analysis of the radiative association (RA) kinetics. Au⁺ associated with C₆F₅H to give Au⁺(C₆F₅H), with binding energy estimated to be 31 kcal mol⁻¹. Au⁺ reacted with C₆H₆ to form the well known Au⁺(C₆H₆) and Au⁺(C₆H₆)₂ complexes. The observation of rapid charge transfer from Au⁺(C₆H₆) to C₆H₆ was interpreted as showing that benzene binds more strongly to neutral Au than to Au⁺. The neutral Au–C₆H₆ bond is accordingly concluded to be stronger than about 70 kcal mol⁻¹.     

Quantum chemical metabolism‐based simulation of carcinogenic potency of benzene derivatives

Using an oxenoid model, we investigated dependences of carcinogenic potency of the benzenes C₆H₅‐X on a nature of substituents X. According to the model, a P450 enzyme breaks a dioxygen molecule and generates the oxens, which readily react with substrates. We suggest that a stability of the intermediate OC₆H₆‐X with tetrahedrally coordinated C atom relative to the molecule C₆H₅‐X determines a rate of substrate biotransformation. Using MO LCAO MNDO approach, we calculated the total energies of molecules C₆H₆‐X and arene oxides OC₆H₆‐X. A difference ΔE_min of these values determines activation energy of oxidation reaction. The compounds with the low ΔE_min values are noncarcinogenic. Benzene derivatives with high ΔE_min values belong to carcinogenic compounds series. The carcinogenicity of amino‐ and nitro‐substituted benzenes is also determined by N‐oxidation of amino and reduction of the nitro group. As the phenylhydroxylamines XC₆H₄NHOH and nitrenium ions XC₆OH₄NH⁺ are the common metabolites of the nitro‐ and amino‐substituted benzenes and nitrenium ions XC₆H₄NH⁺ are the ultimate carcinogens, we use the differences ΔE_N = E(XC₆H₄NH⁺) ? E(XC₆H₄NHOH) as the second parameter characterizing the carcinogenic activity of amino‐ and nitro‐substituted benzenes. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Charge neutralization of ions from benzene

Cross-sections have been measured for the charge neutrilization if ions from benzene in kiloelectron-volt collisions with benezene target molecules. Measured values range from 65 Ų for the symmetric [C₆H₆]⁺? C₆H₆ resonant reactions to 8 Ų for [C₃H₃]⁺? C₆H₆ reactions. Cross-sections computed using a simple resonance charge transfer model compare favourably with experimental data for the symmetric reactions. The cross-sections for asymmetric reactions are smaller that those for they symmetric system and magnitudes of the asymmetric cross-sections are correlated with recombination energies of the respective ions.     

Dissociative Photoionization of 1,4-Dioxane with Tunable VUV Synchrotron Radiation

The photoionization and photodissociation of 1,4-dioxane have been investigated with a reflectron time-of-flight photoionization mass spectrometry and a tunable vacuum ultraviolet synchrotron radiation in the energy region of 8.0-15.5 eV. Parent ion and fragment ions at m/z 88, 87, 58, 57, 45, 44, 43, 41, 31, 30, 29, 28 and 15 are detected under supersonic conditions. The ionization energy of DX as well as the appearance energies of its fragment ions C₄H₇O₂⁺, C₃H₆O⁺, C₃H₅O⁺, C₂H₅O⁺, C₂H₄O⁺, C₂H₃O⁺, C₃H₅⁺, CH₃O⁺, C₂H₆⁺, C₂H₅⁺/CHO⁺, C₂H₄⁺ and CH₃⁺ was determined from their photoionization efficiency curves. The optimized structures for the neutrals, cations, transition states and intermediates related to photodissociation of DX are characterized at the B3LYP/6-31+G(d,p) level and their energies are obtained by G3B3 method. Possible dissociative channels of the DX are proposed based on comparison of experimental AE values and theoretical predicted ones. Intramolecular hydrogen migrations are found to be the dominant processes in most of the fragmentation pathways of 1,4-dioxane.     

Reactions of NO2+ and solvated NO2+ ions with aromatic compounds and alkanes

The rate constants and modes of reaction of NO₂⁺ and C₂H₅ONO₂NO₂⁺ with aromatic compounds and alkanes have been determined in a pulsed ion cyclotron resonance mass spectrometer. Both ions undergo competing charge transfer and substitution reactions (NO₂⁺ + M → MO⁺ + NO; C₂H₅ONO₂NO₂⁺ + M → MNO₂⁺ + C₂H₅ONO₂) with aromatic molecules. In both cases, the probability that a collision results in charge transfer increases with increasing exothermicity of that process. The C₂H₅ONO₂NO₂⁺ ion does not undergo charge transfer with molecules having an ionization potential greater than about 212 kcal/mol (9.2 eV); this observation leads to an estimate of 13 kcal/mol for the binding energy between NO₂⁺ and C₂H₅ONO₂. The importance of the substitution reaction depends on the number of substituents on the aromatic ring and the molecular structure, and, in the case of C₂H₅ONO₂NO₂⁺ ions, on the energetics of the competing charge transfer process. Both NO₂⁺ and C₂H₅ONO₂NO₂⁺ undergo hydride transfer reactions with alkanes. For both these ions, k(hydride transfer)/k (collision) increases with increasing exothermicity of reaction, but in both cases the rate constants of reaction are unusually low when compared with other hydride transfer reactions of comparable exothermicity which have been reported in the literature. This is interpreted as evidence that the attack on the alkane preferentially involves the nitrogen atom (where the charge is localized) rather than one of the oxygen atoms of NO₂⁺.     

Multiple ionization,charge separation and charge stripping reactions involving polycyclic aromatic compounds

The 70 eV electron ionization mass spectra of polycyclic aromatic compounds are characterized by the presence of relatively stable multiply charged molecular ions [M]ⁿ⁺ (n=2–4). When generated from the compounds benzene, napthalene, anthracene, phenanthrene, 2,3-benzanthracene, 1,2-benzanthracene, chrysene, 9,10-benzophenanthrene and pyrene, the relative abundances of the multiply charged ions increase dramatically with the number of rings. These compounds form multiply charged molecular ions (n=2, 3) which undergo unimolecular decompositions indicative of considerable ionic rearrangement. The main charge separation processes observed here [M]²⁺→m₁⁺+m₂⁺, [M]³⁺˙→m₃⁺+m→+m₄²⁺) involve, in almost every case, one or more of the products [CH₃]⁺, [C₂H₃]⁺˙ and [C₃H₃]⁺. This suggests the existence of preferred structures amongst the metastable parent ions. Information on the relative importance of the various fragmentation pathways is presented here along with translational energy release data. Some tentative structural information about the metastable ions has been inferred from the translational energy release on the assumption that the released energy is due primarily to coulombic repulsion within the transition state structure. For the triply charged ions these interpretations have necessitated the use of a coulombic repulsion model which takes account of an extra charge. Vertical ionization energies for the process [M]ⁿ⁺+G→[M](ⁿ⁺¹)+G+e^? (charge stripping) have also been determined where possible for n=1 and 2 and the results from these experiments allow the derivation of simple empirical equations which relate successive ionization energies for the formation of [M]²⁺ and [M]³⁺˙ to the appearance energy of [M]⁺˙.     

Benzene as a selective chemical ionization reagent gas

Dilute mixtures of C₆H₆ or C₆D₆ in He provide abundant [C₆H₆]^+· or [C₆D₆]^+· ions and small amounts of [C₆H₇]⁺ or [C₆D₇]⁺ ions as chemical ionization (CI) reagent ions. The C₆H₆ or C₆D₆ CI spectra of alkylbenzenes and alkylanilines contain predominantly M^+· ions from reactions of [C₆H₆]^+· or [C₆D₆]^+· and small amounts of MH⁺ or MD⁺ ions from reactions of [C₆H₇]⁺ or [C₆D₇]⁺. Benzene CI spectra of aliphatic amines contain M^+·, fragment ions and sample-size-dependent MH⁺ ions from sample ion-sample molecules reactions. The C₆D₆ CI spectra of substituted pyridines contain M^+· and MD⁺ ions in different ratios depending on the substituent (which alters the ionization energy of the substituted pyridine), as well as sample-size-dependent MH⁺ ions from sample ion-sample molecule reactions. Two mechanisms are observed for the formation of MD⁺ ions: proton transfer from [C₆D₆]^+· or charge transfer from [C₆D₆]^+· to give M^+·, followed by deuteron transfer from C₆D₆ to M^+·. The mechanisms of reactions were established by ion cyclotron resonance (ICR) experiments. Proton transfer from [C₆H₆]^+· or [C₆D₆]^+· is rapid only for compounds for which proton transfer is exothermic and charge transfer is endothermic. For compounds for which both charge transfer and proton transfer are exothermic, charge transfer is the almost exclusive reaction.     

Gas-phase electrophilic addition reactions of chlorinated alkyl ions

Chemical ionization was used to study gas-phase electrophilic addition reactions of chloromethyl ions ([CH_xCl_3-x]⁺, x = 0, 1, 2) with a number of substituted benzenes (C₆H₅Y, Y = NH₂, OH, CHO, CN, NO₂). Mass-analyzed ion kinetic energy spectrometry was used to characterize the reaction products with respect to the site of electrophilic addition (ring v. substituent). In some cases, examination of secondary reaction products (ion–molecule adduct which has undergone an elimination reaction in the ion source) aided in establishing the original site of electrophilic addition. Aniline, benzonitrile and nitrobenzene exhibited preferential substituent interaction, while phenol and benzaldehyde gave a mixture of ring and substituent reaction products. These gas-phase results differ considerably from solution-phase Friedel–Crafts alkylation; however, they are consistent with the notion of preferential σ-bond formation at polarizable centers of negative charge.     

Molecular fragmentations: Part VII. Structures and energies of mass spectral fragments derived from benzene

Molecular structures and energies have been calculated, using MINDO/3, of the mass spectral ions arising from benzene: (C₆H₆)⁺ (three non-valence isomers); (C₆H₅⁺); (C₅H₃⁺) (four isomers); (C₄H₄)⁺ (three isomers); (C₄H₃)⁺ (two isomers); (C₄H₂)⁺ (four isomers); (C₃H₃)⁺; and (C₂H₂)⁺. Calculations have been made for the conjugate neutral fragments, allowing calculation of appearance potentials, and also for the ion (C₆H₇)⁺.