Heats of Formation of C6H5•, C6H5+, and C6H5NO by Threshold Photoelectron Photoion Coincidence and Active Thermochemical Tables Analysis

Dialumene as a Dimeric or Monomeric Al Synthon for C-F Activation in Monofluorobenzene

The activation of C-F bonds has long been regarded as the subject of research in organometallic chemistry, given their synthetic relevance and the fact that fluorine is the most abundant halogen in the Earth's crust. However, C-F bond activation remains a largely unsolved challenge due to the high bond dissociation energies, which was historically dominated by transition metal complexes. Main group elements that can cleave unactivated monofluorobenzene are still quite rare and restricted to s-block complexes with a biphilic nature. Herein, we demonstrate an Al-mediated activation of monofluorobenzene using a neutral dialumene, allowing for the synthesis of the formal oxidative addition products at either double or single aluminum centers. This neutral dialumene system introduces a novel methodology for C-F bond activation based on formal oxidative addition and reductive elimination processes around the two aluminum centers, as demonstrated by combined experimental and computational studies. A "masked" alumylene was unprecedentedly synthesized to prove the proposed reductive elimination pathway. Furthermore, the synthetic utility is highlighted by the functionalization of the resulting aryl-aluminum compounds.

The elusive phenylethynyl radical and its cation: synthesis, electronic structure, and reactivity

Alkynyl radicals and cations are crucial reactive intermediates in chemistry, but often evade direct detection. Herein, we report the direct observation of the phenylethynyl radical (C6H5CC˙) and its cation (C6H5CC+), which are two of the most reactive intermediates in organic chemistry. The radical is generated via pyrolysis of (bromoethynyl)benzene at temperatures above 1500 K and is characterized by photoion mass-selected threshold photoelectron spectroscopy (ms-TPES). Photoionization of the phenylethynyl radical yields the phenylethynyl cation, which has never been synthesized due to its extreme electrophilicity. Vibrationally-resolved ms-TPES assisted by ab initio calculations unveiled the complex electronic structure of the phenylethynyl cation, which appears at an adiabatic ionization energy (AIE) of 8.90 ± 0.05 eV and exhibits an uncommon triplet (3B1) ground state, while the closed-shell singlet (1A1) state lies just 2.8 kcal mol-1 (0.12 eV) higher in energy. The reactive phenylethynyl radical abstracts hydrogen to form ethynylbenzene (C6H5CCH) but also isomerizes via H-shift to the o-, m-, and p-ethynylphenyl isomers (C6H4CCH). These radicals are very reactive and undergo ring-opening followed by H-loss to form a mixture of C8H4 triynes, along with low yields of cyclic 3- and 4-ethynylbenzynes (C6H3CCH). At higher temperatures, dehydrogenation from the unbranched C8H4 triynes forms the linear tetraacetylene (C8H2), an astrochemically relevant polyyne.

Constructing Ultrastable Metallo-Cages via In Situ Deprotonation/Oxidation of Dynamic Supramolecular Assemblies

Coordination-driven self-assembly enables the spontaneous construction of metallo-supramolecules with high precision, facilitated by dynamic and reversible metal-ligand interactions. The dynamic nature of coordination, however, results in structural lability in many metallo-supramolecular assembly systems. Consequently, it remains a formidable challenge to achieve self-assembly reversibility and structural stability simultaneously in metallo-supramolecular systems. To tackle this issue, herein, we incorporate an acid-/base-responsive tridentate ligand into multitopic building blocks to precisely construct a series of metallo-supramolecular cages through coordination-driven self-assembly. These dynamic cagelike assemblies can be transformed to their static states through mild in situ deprotonation/oxidation, leading to ultrastable skeletons that can withstand high temperatures, metal ion chelators, and strong acid/base conditions. This in situ transformation provides a reliable and powerful approach to manipulate the kinetic features and stability of metallo-supramolecules and allows for modulation of encapsulation and release behaviors of metallo-cages when utilizing nanoscale quantum dots (QDs) as guest molecules.

Active Thermochemical Tables: Enthalpies of Formation of Bromo- and Iodo-Methanes, Ethenes and Ethynes

The thermochemistry of halocarbon species containing iodine and bromine is examined through an extensive interplay between new Feller-Peterson-Dixon (FPD) style composite methods and a detailed analysis of all available experimental and theoretical determinations using the thermochemical network that underlies the Active Thermochemical Tables (ATcT). From the computational viewpoint, a slower convergence of the components of composite thermochemistry methods is observed relative to species that solely contain first row elements, leading to a higher computational expense for achieving comparable levels of accuracy. Potential systematic sources of computational uncertainty are investigated, and, not surprisingly, spin-orbit coupling is found to be a critical component, particularly for iodine containing molecular species. The ATcT analysis of available experimental and theoretical determinations indicates that prior theoretical determinations have significantly larger uncertainties than originally reported, particularly in cases where molecular spin-orbit effects were ignored. Accurate and reliable heats of formation are reported for 38 halogen containing systems, based on combining the current computations with previous experimental and theoretical work via the ATcT approach.

Resolving the F2 bond energy discrepancy using coincidence ion pair production (cipp) spectroscopy

Coincidence ion pair production (cipp) spectra of F₂ were recorded on the DELICIOUS III coincidence spectrometer in the one-photon excitation region of 125 975-126 210 cm^-1. The F⁺ + F^- signal shows a rotational band head structure, corresponding to F₂ Rydberg states crossing over to the ion pair production surface. Spectral simulation and quantum defect analysis allowed the characterization of five new molecular Rydberg states (F₂**): one Π and four Σ states. The lowest-energy Rydberg state spectrum observed (T₀ = 125 999 cm^-1) lacked some of the predicted rotational structure, which allowed an accurate determination of the ion pair production threshold of 15.6229₄± 0.0004₃ eV. Using the well-known atomic fluorine ionization energy and electron affinity, this number leads to a ground state F-F dissociation energy of 1.6012₉± 0.0004₄ eV. Photoelectron photoion coincidence (PEPICO) experiments were also carried out on F₂ and the dissociative photoionization threshold to F⁺ + F was determined as 19.0242 ± 0.0006 eV. Using the atomic fluorine ionization energy, this can be converted to an F₂ dissociation energy of 1.6013₂± 0.0006₂ eV, further confirming the cipp-derived value above. Because the two experiments were independently energy-calibrated, they can be averaged to 1.6013₀± 0.0003₆ eV and this value can be used to derive the fluorine atom's 0 K heat of formation as 77.25₁± 0.01₇ kJ mol^-1. This latter is in excellent agreement with the latest Active Thermochemical Table (ATcT) value but improves its accuracy by almost a factor of three.

Modeling Gas-Phase Unimolecular Dissociation for Bond Dissociation Energies: Comparison of Statistical Rate Models within RRKM Theory

The Rice-Ramsperger-Kassel-Marcus (RRKM) theory provides a simple yet powerful rate theory for calculating microcanonical rate constants. In particular, it has found widespread use in combination with gas-phase kinetic experiments of unimolecular dissociations to extract experimental bond dissociation energies (BDEs). We have previously found several discrepancies between the computed BDE values and the respective experimental ones, obtained with our empirical rate model, named L-CID. To investigate the reliability of our rate model, we conducted a theoretical analysis and comparison of the performance of conventional rate models and L-CID within the RRKM framework. Using the previously published microcanonical rate data as well as reaction cross-section data, we show that the BDE values obtained with the L-CID model agree with the ones from the other rate models within the expected uncertainty bounds. Based on this agreement, we discuss the possible rationalization of the good performance of the L-CID model.

Bond Dissociation Energies in the Gas Phase for Large Molecular Ions by Threshold Collision-Induced Dissociation Experiments: Stretching the Limits

Accurate bond dissociation energies for large molecules are difficult to obtain by either experimental or computational methods. The former methods are hampered by a range of physical and practical limitations in gas-phase measurement techniques, while the latter require incorporation of multiple approximations whose impact on accuracy may not always be clear. When internal benchmarks are not available, one hopes that experiment and theory can mutually support each other. A recent report found, however, a large discrepancy between gas-phase bond dissociation energies, measured mass spectrometrically, and the corresponding quantities computed using density functional theory (DFT)-D3 and DLPNO-CCSD(T) methods. With the widespread application of these computational methods to large molecular systems, the discrepancy needs to be resolved. We report a series of experimental studies that validate the mass spectrometric methods from small to large ions and find that bond dissociation energies extracted from threshold collision-induced dissociation experiments on large ions do indeed behave correctly. The implications for the computational studies are discussed.

Thermochemistry of Intermediates and Products in the Oxidation Reaction of 1,1,2-Trifluoroethene via OH Radical

Density functional theory (DFT) and composite ab initio based calculations are performed on trifluoroethane along with intermediate radicals, parent molecules of the radicals, and products related to the reaction of hydroxyl radical with 1,1,2-trifluoroethene, as a reference for hydrofluoroolefins (HFO). Potential energy barriers for internal rotations have been computed. Calculated torsional potentials are incorporated into the determination of entropy, S°₂₉₈, and heat capacities as a function of temperature, C_p(T), for each target molecule. Six isodesmic or isogyric reactions and five calculation methods are used to determine heats of formation at 298 K (Δ_fH₂₉₈) in kcal mol^-1 of each target species. The CBS-APNO method shows the best agreement with experimental data in comparisons from 16 reference reactions on Δ_rxnH of each method. The lowest configuration structures of each target species are reported. Intramolecular hydrogen bonds between the hydroxyl hydrogen atom and the fluorine atom on the adjacent carbon can stabilize molecules by up to 3 kcal mol^-1. R-OH bond dissociation energies are observed to increase with the number of fluorine atoms on the carbon connected to hydroxy group. Recommended Δ_fH₂₉₈ values in kcal mol^-1 derived from the most stable conformers are CF₂(OH)CH₂F (-213.0), CF₂(O^•)CH₂F (-148.6), CF₂(OH)C^•FH (-162.4), CHF₂CHFOH (-207.5), CHF₂C^•FOH (-158.3), C^•F₂CHFOH (-155.5), CHF₂CHFO^• (-150.4), CF₃CH₂OH (-212.5), and CF₃C^•HOH (-167.9).

Thermal Decomposition of Potential Ester Biofuels. Part I: Methyl Acetate and Methyl Butanoate

Two methyl esters were examined as models for the pyrolysis of biofuels. Dilute samples (0.06-0.13%) of methyl acetate (CH₃COOCH₃) and methyl butanoate (CH₃CH₂CH₂COOCH₃) were entrained in (He, Ar) carrier gas and decomposed in a set of flash-pyrolysis microreactors. The pyrolysis products resulting from the methyl esters were detected and identified by vacuum ultraviolet photoionization mass spectrometry. Complementary product identification was provided by matrix infrared absorption spectroscopy. Pyrolysis pressures in the pulsed microreactor were about 20 Torr and residence times through the reactors were roughly 25-150 μs. Reactor temperatures of 300-1600 K were explored. Decomposition of CH₃COOCH₃ commences at 1000 K, and the initial products are (CH₂═C═O and CH₃OH). As the microreactor is heated to 1300 K, a mixture of CH₂═C═O and CH₃OH, CH₃, CH₂═O, H, CO, and CO₂ appears. The thermal cracking of CH₃CH₂CH₂COOCH₃ begins at 800 K with the formation of CH₃CH₂CH═C═O and CH₃OH. By 1300 K, the pyrolysis of methyl butanoate yields a complex mixture of CH₃CH₂CH═C═O, CH₃OH, CH₃, CH₂═O, CO, CO₂, CH₃CH═CH₂, CH₂CHCH₂, CH₂═C═CH₂, HCCCH₂, CH₂═C═C═O, CH₂═CH₂, HC≡CH, and CH₂═C═O. On the basis of the results from the thermal cracking of methyl acetate and methyl butanoate, we predict several important decomposition channels for the pyrolysis of fatty acid methyl esters, R-CH₂-COOCH₃. The lowest-energy fragmentation will be a 4-center elimination of methanol to form the ketene RCH═C═O. At higher temperatures, concerted fragmentation to radicals will ensue to produce a mixture of species: (RCH₂ + CO₂ + CH₃) and (RCH₂ + CO + CH₂═O + H). Thermal cracking of the β C-C bond of the methyl ester will generate the radicals (R and H) as well as CH₂═C═O + CH₂═O. The thermochemistry of methyl acetate and its fragmentation products were obtained via the Active Thermochemical Tables (ATcT) approach, resulting in Δ_fH₂₉₈(CH₃COOCH₃) = -98.7 ± 0.2 kcal mol^-1, Δ_fH₂₉₈(CH₃CO₂) = -45.7 ± 0.3 kcal mol^-1, and Δ_fH₂₉₈(COOCH₃) = -38.3 ± 0.4 kcal mol^-1.

A vacuum ultraviolet laser pulsed field ionization-photoion study of methane (CH4): determination of the appearance energy of methylium from methane with unprecedented precision and the resulting impact on the bond dissociation energies of CH4and CH4+

High-resolution VUV laser PFI-PI detection method for the study of quantum-state-selected unimolecular ion dissociation.

Advances in threshold photoelectron spectroscopy (TPES) and threshold photoelectron photoion coincidence (TPEPICO)

The history and evolution of molecular threshold photoelectron spectroscopy and threshold photoelectron photoion coincidence spectroscopy (TPEPICO) over the last fifty years are reviewed.

Thermochemistry of icosahedral closo-dicarboranes: a composite ab initio quantum-chemical perspective

We obtained accurate thermochemical properties for the ortho-, meta-, and para-dicarborane isomers (C2B10H12) by means of explicitly correlated high-level thermochemical procedures. The thermochemical properties include heats of formation, isomerization energies, C–H and B–H bond dissociation energies (BDEs), and ionization potentials. Of these only the ionization potentials are known experimentally. Our best theoretical ionization potentials, obtained by means of the ab initio W1–F12 thermochemical protocol, was 241.50 kcal mol–1 (para-dicarborane), 238.45 kcal mol–1 (meta-dicarborane), and 236.54 kcal mol–1 (ortho-dicarborane). These values agree with the experimental values adopted by the National Institute of Standards and Technology (NIST) thermochemical tables to within overlapping uncertainties. However, they suggest that the experimental values may represent significant underestimations. For all isomers, the C–H BDEs are systematically higher than the B–H BDEs because of the relative stability of the boron-centred radicals. The C–H BDEs for the three isomers cluster within a narrow energetic interval, namely between 110.8 kcal mol–1 (para-dicarborane) and 111.7 kcal mol–1 (meta-dicarborane). The B–H BDEs cluster within a larger interval ranging between 105.8 and 108.1 kcal mol–1 (both obtained for ortho-dicarborane). We used our benchmark W1–F12 data to assess the performance of a number of lower cost composite ab initio methods. We found that the Gaussian-3 procedures (G3(MP2)B3 and G3B3) result in excellent performance with overall root-mean-square deviations (RMSDs) of 0.3–0.4 kcal mol–1 for the isomerization, ionization, and bond dissociation energies. However, the Gaussian-4 procedures (G4, G4(MP2), and G4(MP2)-6X) showed relatively poor performance with overall RMSDs of 1.3–3.7 kcal mol–1.

First-Principles Prediction of Enthalpies of Formation for Polycyclic Aromatic Hydrocarbons and Derivatives

In this article, the first-principles prediction of enthalpies of formation is demonstrated for 669 polycyclic aromatic hydrocarbon (PAH) compounds and a number of related functionalized molecules. It is shown that by extrapolating density functional theory calculations to a large basis set limit and then applying a group based correction scheme that good results may be obtained. Specifically, a mean unsigned deviation and root mean squared deviation from the experimental enthalpies of formation data of 5.0 and 6.4 kJ/mol, respectively, are obtained using this scheme. This computational scheme is economical to compute and straightforward to apply, while yielding results of reasonable reliability. The results are also compared for a smaller set of molecules to the predictions given by the G3B3 and G3MP2B3 variants of the Gaussian-3 model chemistry with a mean unsigned deviation and root mean squared deviation from the experimental enthalpies of formation of 4.5 and 4.8 kJ/mol, respectively.

Are the three hydroxyphenyl radical isomers created equal?--The role of the phenoxy radical

We have investigated the thermal decomposition of the three hydroxyphenyl radicals (˙C6H4OH) in a heated microtubular reactor. Intermediates and products were identified isomer-selectively applying photoion mass-selected threshold photoelectron spectroscopy with vacuum ultraviolet synchrotron radiation. Similarly to the phenoxy radical (C6H5-O˙), hydroxyphenyl decomposition yields cyclopentadienyl (c-C5H5) radicals in a decarbonylation reaction at elevated temperatures. This finding suggests that all hydroxyphenyl isomers first rearrange to form phenoxy species, which subsequently decarbonylate, a mechanism which we also investigate computationally. Meta- and para-radicals were selectively produced and spectroscopically detectable, whereas the ortho isomer could not be traced due to its fast rethermalization and rapid decomposition in the reactor. A smaller barrier to isomerization to phenoxy was found to be the reason for this observation. Since hydroxyphenyl species may be present under typical sooting conditions in flames, the resonantly stabilized cyclopentadienyl radical adds to the hydrocarbon pool and can contribute to the formation of polycyclic aromatic hydrocarbons, which are precursors in soot formation.

A Database of Formation Enthalpies of Nitrogen Species by Compound Methods (CBS-QB3, CBS-APNO, G3, G4)

Accurate thermochemical data for compounds containing C/H/N/O are required to underpin kinetics simulation and modeling of the reactions of these species in different environments. There is a dearth of experimental data so computational quantum chemistry has stepped in to fill this breach and to verify whether particular experiments are in need of revision. A number of composite model chemistries (CBS-QB3, CBS-APNO, G3, and G4) are used to compute theoretical atomization energies and hence enthalpies of formation at 0 and 298.15 K, and these are benchmarked against the best available compendium of values, the Active Thermochemical Tables or ATcT. In general the agreement is very good for some 28 species with the only discrepancy being for hydrazine. It is shown that, although individually the methods do not perform that well, collectively the mean unsigned error is <1.7 kJ mol(-1); hence, this approach provides a useful tool to screen published values and validate new experimental results. Using multiple model chemistries does have some drawbacks but can produce good results even for challenging molecules like HOON and CN2O2. The results for these smaller validated molecules are then used as anchors for determining the formation enthalpies of larger species such as methylated hydrazines and diazenes, five- and six-membered heterocyclics via carefully chosen isodesmic working reactions with the aim of resolving some discrepancies in the literature and establishing a properly validated database. This expanded database could be useful in testing the performance of computationally less-demanding density function methods with newer functionals that have the capacity to treat much larger systems than those tested here.

Pathways for the OH + Cl2 → HOCl + Cl and HOCl + Cl → HCl + ClO Reactions

High level coupled-cluster theory, with spin-orbit coupling evaluated via the Breit-Pauli operator in the interacting-states approach, is used to investigate the OH radical reaction with Cl2 and the subsequent reaction HOCl + Cl. The entrance complex, transition state, and exit complex for both reactions have been determined using the CCSD(T) method with correlation consistent basis sets up to cc-pV6Z. Also reported are CCSDT computations. The OH + Cl2 reaction is predicted to be endothermic by 2.2 kcal/mol, compared to the best experiments, 2.0 kcal/mol. The above theoretical results include zero-point vibrational energy corrections and spin-orbit contributions. The activation energy (Ea) of the OH + Cl2 reaction predicted here, 2.3 kcal/mol, could be as much as 1 kcal/mol too high, but it falls among the four experimental Ea values, which span the range 1.1-2.5 kcal/mol. The exothermicity of the second reaction HOCl + Cl → HCl + ClO is 8.4 kcal/mol, compared to experiment 8.7 kcal/mol. The activation energy for latter reaction is unknown experimentally, but predicted here to be large, 11.5 kcal/mol. There are currently no experiments relevant to the theoretical entrance and exit complexes predicted here.

Active Thermochemical Tables: Sequential Bond Dissociation Enthalpies of Methane, Ethane, and Methanol and the Related Thermochemistry

Active Thermochemical Tables (ATcT) thermochemistry for the sequential bond dissociations of methane, ethane, and methanol systems were obtained by analyzing and solving a very large thermochemical network (TN). Values for all possible C-H, C-C, C-O, and O-H bond dissociation enthalpies at 298.15 K (BDE298) and bond dissociation energies at 0 K (D0) are presented. The corresponding ATcT standard gas-phase enthalpies of formation of the resulting CHn, n = 4-0 species (methane, methyl, methylene, methylidyne, and carbon atom), C2Hn, n = 6-0 species (ethane, ethyl, ethylene, ethylidene, vinyl, ethylidyne, acetylene, vinylidene, ethynyl, and ethynylene), and COHn, n = 4-0 species (methanol, hydroxymethyl, methoxy, formaldehyde, hydroxymethylene, formyl, isoformyl, and carbon monoxide) are also presented. The ATcT thermochemistry of carbon dioxide, water, hydroxyl, and carbon, oxygen, and hydrogen atoms is also included, together with the sequential BDEs of CO2 and H2O. The provenances of the ATcT enthalpies of formation, which are quite distributed and involve a large number of relevant determinations, are analyzed by variance decomposition and discussed in terms of principal contributions. The underlying reasons for periodic appearances of remarkably low and/or unusually high BDEs, alternating along the dissociation sequences, are analyzed and quantitatively rationalized. The present ATcT results are the most accurate thermochemical values currently available for these species.

Thermochemical factors affecting the dehalogenation of aromatics

Halogenated aromatics are one of the largest chemical classes of environmental contaminants, and dehalogenation remains one of the most important processes by which these compounds are degraded and detoxified. The thermodynamic constraints of aromatic dehalogenation reactions are thus important for understanding the feasibility of such reactions and the redox conditions necessary for promoting them. Accordingly, the thermochemical properties of the (poly)fluoro-, (poly)chloro-, and (poly)bromobenzenes, including standard enthalpies of formation, bond dissociation enthalpies, free energies of reaction, and the redox potentials of Ar-X/Ar-H couples, were investigated using a validated density functional protocol combined with continuum solvation calculations when appropriate. The results highlight the fact that fluorinated aromatics stand distinct from their chloro- and bromo- counterparts in terms of both their relative thermodynamic stability toward dehalogenation and how different substitution patterns give rise to relevant properties, such as bond strengths and reduction potentials.

A classical trajectory study of the dissociation and isomerization of C2H5

Motivated by photodissociation experiments in which non-RRKM nanosecond lifetimes of the ethyl radical were reported, we have performed a classical trajectory study of the dissociation and isomerization of C2H5 over the energy range 100-150 kcal/mol. We used a customized version of the AIREBO semiempirical potential (Stuart, S. J.; et al. J. Chem. Phys. 2000, 112, 6472-6486) to more accurately describe the gas-phase decomposition of C2H5. This study constitutes one of the first gas-phase applications of this potential form. At each energy, 10,000 trajectories were run and all underwent dissociation in less than 100 ps. The calculated dissociation rate constants are consistent with RRKM models; no evidence was found for nanosecond lifetimes. An analytic kinetics model of isomerization/dissociation competition was developed that incorporated incomplete mode mixing through a postulated divided phase space. The fits of the model to the trajectory data are good and represent the trajectory results in detail through repeated isomerizations at all energies. The model correctly displays single exponential decay at lower energies, but at higher energies, multiexponential decay due to incomplete mode mixing becomes more apparent. At both ends of the energy range, we carried out similar trajectory studies on CD2CH3 to examine isotopic scrambling. The results largely support the assumption that a H or a D atom is equally likely to dissociate from the mixed-isotope methyl end of the molecule. The calculated fraction of products that have the D atom dissociation is ∼20%, twice the experimental value available at one energy within our range. The calculated degree of isotopic scrambling is non-monotonic with respect to energy due to a non-monotonic ratio of the isomerization to dissociation rate constants.