Ring Conserved Isodesmic Reactions: A New Method for Estimating the Heats of Formation of Aromatics and PAHs

High-Throughput Predictions of Accurate Enthalpies of Formation for Larger Molecules Utilizing the Bond Difference Correction Method

The prediction of standard enthalpies of formation (EOFs) for larger molecules involves a trade-off between accuracy and cost, often resulting in non-negligible errors. The connectivity-based hierarchy (CBH) and simple bond additivity correction (BAC) are two promising means for evaluating EOFs, although they cannot achieve strict chemical accuracy. Calculated errors in the CBH are confirmed from accumulated systematic errors associated with bond differences in chemical environments. On the basis of a new set of bond descriptors, our developed bond difference correction (BDC) method effectively solves incremental errors with molecular size and inability applications for aromatic molecules. To balance the accuracy between non-aromatic and aromatic molecules, a more accurate BAC-based method with unpaired electrons and p hybrid orbitals (BAC-EP) is developed. With the incorporation of the two methods above, strict chemical accuracy by the largest deviation is achieved at low costs. These universal, ultrafast, and high-throughput methods greatly contribute to self-consistent thermodynamic parameters in combustion mechanisms.

Benchmark calculations for bond dissociation energies and enthalpy of formation of chlorinated and brominated polycyclic aromatic hydrocarbons

Thermodynamic properties, i.e., bond dissociation energies and enthalpy of formation, of chlorinated and brominated polycyclic aromatic hydrocarbons play a fundamental role in understanding their formation mechanisms and reactivity. Computational electronic structure calculations routinely used to predict thermodynamic properties of various species are limited for these compounds due to large computational cost to obtain accurate results by employing high-level wave function theory methods. In this work, a number of composite model chemistry methods (CBS-QB3, G3MP2, G3, and G4) are used to compute bond dissociation energies and enthalpies of formation of small to medium-size chlorinated and brominated polycyclic aromatic hydrocarbon compounds. The enthalpy of formation is derived via the atomization method and compared against the recommended values. Statistical analysis indicates that G4 is the best method. For comparison, three commonly used density functional theory (DFT) methods (M06-2X, ωB97X-D and B2PLYP-D3) with various basis sets including 6-311++G(d, p), cc-pVTZ, and cc-pVQZ in the prediction of bond dissociation energies and enthalpies of formation have been tested using the optimized geometries at the same M06-2X/6-311++G(d, p) level of theory. It is found that ωB97X-D/6-311++G(d, p) shows the best performance in computing the bond dissociation energies, while ωB97X-D/cc-pVTZ exhibits the best prediction in enthalpy of formation of the studied reaction systems. The structural effect on the bond dissociation energies and enthalpy of formation of chlorinated and brominated polycyclic aromatic hydrocarbons are then systematically analyzed. Based on comparisons of the various methods, reliable DFT methods are recommended for future theoretical studies on large chlorinated and brominated polycyclic aromatic hydrocarbons considering both accuracy and computational cost. This work, to the authors' knowledge, is the first to systematically benchmark theoretical methods for the accurate prediction of thermodynamic properties for chlorinated and brominated polycyclic aromatic hydrocarbons.

Benchmark calculations using state-of-the-art DFT functionals and composite methods for bond dissociation energy and enthalpy of formation of halogenated polycyclic aromatic hydrocarbons are performed.

Polycyclic aromatic hydrocarbons: from small molecules through nano-sized species towards bulk graphene

We have examined the use of systematic bond-separation reactions and purposely constructed chemistry-preserving isodesmic reactions for the thermochemical calculation of aromatic hydrocarbon species. The bond-separation approach yields somewhat disappointing accuracy even when the reaction energies are obtained with generally robust composite and double-hybrid (DH) density functional theory (DFT) methods. In contrast, for the purposely constructed reactions, we find a dramatic improvement in the accuracy for energies calculated with all methods examined. Notably, for medium-sized aromatic hydrocarbons, we find that an effective approach for formulating a well-balanced reaction is to split the target species into two halves with an aromatic overlapping region. Overall, the G4(MP2)-XK, MPW2PLYP, MN15, PBE, and DC-DFTB3 methods are reasonable within their respective classes of methods for the calculation of bond-separation as well as chemistry-preserving isodesmic reactions. We have further computed per-carbon atomization energy (AE) for a series of D_6h benzene-type molecules, and thus obtained a formula for extrapolation to the graphene limit [AE_n = 711.5 × (1 - 1/n^0.640) kJ mol^-1, where n = number of carbons]. It suggests that nano-graphene with a length larger than 10 nm would resemble properties of bulk graphene, and conversely, downsizing a nano-graphene beyond this point may lead to considerably altered properties from the bulk.

Unveiling the role of intra and interatomic interactions in the energetics of reaction schemes: a quantum chemical topology analysis

In this work we present a detailed analysis of selected reaction schemes in terms of the atomic components of the electronic energy defined by the quantum theory of atoms in molecules and the interacting quantum atoms method. The aim is to provide an interpretation tool for the energy change involved in a chemical reaction by means of the atomic and interaction contributions to the energies of the molecules involved. Ring strain in cyclic alkanes, the resonance energy of aromatic and antiaromatic molecules, local aromaticity in polycyclic aromatic hydrocarbons, intermolecular bonding in hydrogen fluoride clusters, and hydration of d-block metal dications were selected for the study. It was found that in addition to the changes in the strong C-C interactions in the carbon skeleton of the organic molecular rings, other contributions not usually considered to be important such as those between C and H atoms (either bonded or not) need to be considered in order to account for the net energy changes. The analysis unveils the role of the ionic and covalent contributions to the hydrogen bonding in HF clusters and the energetic origin and extent of cooperative effects involved. Moreover, the "double-hump" behavior observed for the hydration energy trend of [M(H2O)6]2+ complexes is explained in terms of the deformation energy of the metal cation and the increasingly covalent metal-water interactions. In addition, proper comparisons with the description provided by other methodologies are briefly discussed. The topological approach proposed in this contribution proves to be useful for the description of energy changes of apposite reaction schemes in chemically meaningful terms.

Gas-Phase Heat of Formation Values for Buckminsterfullerene (C60), C70 Fullerene (C70), Corannulene, Coronene, Sumanene, and Other Polycyclic Aromatic Hydrocarbons Calculated Using Density Functional Theory (M06 2X) Coupled with a Versatile Inexpensive Group-Equivalent Approach

A straightforward procedure using density functional theory (M06 2X) coupled with a group-equivalent approach is described that was used to calculate gas-phase heat of formation (Δ_f H°_g,298) values for buckminsterfullerene (C₆₀), C70 fullerene (C₇₀), corannulene, coronene, and sumanene. This procedure was also used to calculate exceptionally accurate Δ_f H°_g,298 values for a variety of single-ring aromatic and 2-7 ring polycyclic aromatic hydrocarbons (PAHs) as well as a large selection of other hydrocarbons and phenols. The approach described herein is internally consistent, and results for C₆₀, C₇₀, corannulene, coronene, and sumanene are in very close agreement with results reported by others who used higher-level computational theory. Statistical analysis of a test set containing benzene and 18 two to seven ring PAHs demonstrated that by using this approach a mean absolute deviation (MAD) and a root-mean-square deviation (RMSD) of 0.8 and 1.3 kJ/mol, respectively, were achieved for reference/experimental Δ_f H°_g,298 values versus calculated/predicted Δ_f H°_g,298 values. For statistical analysis of a larger test set containing 235 aromatic and aliphatic hydrocarbons and phenols, a MAD and a RMSD of 1.2 and 1.9 kJ/mol, respectively, were achieved for reference/experimental Δ_f H°_g,298 values versus calculated/predicted Δ_f H°_g,298 values.

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.

Chemistry of polycyclic aromatic hydrocarbons formation from phenyl radical pyrolysis and reaction of phenyl and acetylene

An experimental investigation of phenyl radical pyrolysis and the phenyl radical + acetylene reaction has been performed to clarify the role of different reaction mechanisms involved in the formation and growth of polycyclic aromatic hydrocarbons (PAHs) serving as precursors for soot formation. Experiments were conducted using GC/GC-MS diagnostics coupled to the high-pressure single-pulse shock tube present at the University of Illinois at Chicago. For the first time, comprehensive speciation of the major stable products, including small hydrocarbons and large PAH intermediates, was obtained over a wide range of pressures (25-60 atm) and temperatures (900-1800 K) which encompass the typical conditions in modern combustion devices. The experimental results were used to validate a comprehensive chemical kinetic model which provides relevant information on the chemistry associated with the formation of PAH compounds. In particular, the modeling results indicate that the o-benzyne chemistry is a key factor in the formation of multi-ring intermediates in phenyl radical pyrolysis. On the other hand, the PAHs from the phenyl + acetylene reaction are formed mainly through recombination between single-ring aromatics and through the hydrogen abstraction/acetylene addition mechanism. Polymerization is the common dominant process at high temperature conditions.

Neutral, ion gas-phase energetics and structural properties of hydroxybenzophenones

We have carried out a study of the energetics, structural, and physical properties of o-, m-, and p-hydroxybenzophenone neutral molecules, C(13)H(10)O(2), and their corresponding anions. In particular, the standard enthalpies of formation in the gas phase at 298.15 K for all of these species were determined. A reliable experimental estimation of the enthalpy associated with intramolecular hydrogen bonding in chelated species was experimentally obtained. The gas-phase acidities (GA) of benzophenones, substituted phenols, and several aliphatic alcohols are compared with the corresponding aqueous acidities (pK(a)), covering a range of 278 kJ.mol(-1) in GA and 11.4 in pK(a). A computational study of the various species shed light on structural effects and further confirmed the self-consistency of the experimental results.

Thermochemistry and gas-phase ion energetics of 2-hydroxy-4-methoxy-benzophenone (oxybenzone)

We have investigated the thermochemistry and ion energetics of the oxybenzone (2-hydroxy-4-methoxy-benzophenone, C14H12O3, 1H) molecule. The following parameters have been determined for this species: gas-phase enthalpy for the of neutral molecule at 298.15K, (Delta(f)H0(m)(g) = -303.5 +/- 5.1 kJ x mol-1), the intrinsic (gas-phase) acidity (GA(1H) = 1402.1 +/- 8.4 kJ x mol-1), enthalpy of formation for the oxybenzone anion (Delta(f)H0(m)(1-,g) = -402.3 +/- 9.8 kJ x mol-1). We also have obtained the enthalpy of formation of, 4-hydroxy-4'-methoxybenzophenone (Delta(f)H0(m)(g) = -275.4 +/- 10 kJ x mol-1) and 3-methoxyphenol anion (Delta(f)H0(m)(C7H7O2-,g) = -317.7 +/- 8.7 kJ x mol-1). A reliable experimental estimation of enthalpy related to intramolecular hydrogen bonding in oxybenzone has also been obtained (30.1 +/- 6.3 kJ x mol-1) and compared with our theoretical calculations at the B3LYP/6-311++G** level of theory, by means of an isodesmic reaction scheme. In addition, heat capacities, temperature, and enthalpy of fusion have been determined for this molecule by differential scanning calorimetry.

Integrating process safety with molecular modeling-based risk assessment of chemicals within the REACH regulatory framework: benefits and future challenges

Registration, evaluation and authorization of chemicals (REACH) represents a recent regulatory initiative by the European union commission to protect human health and the environment from potentially hazardous chemicals. Under REACH, all stakeholders must submit (thermo)physical, thermochemical, and toxicological data for certain chemicals. The commission's impact assessment studies estimate that the costs of REACH will be approximately 3-5 billion Euros. The present study advocates the systematic incorporation of computational chemistry and computer-assisted chemical risk assessment methods into REACH to reduce regulatory compliance costs. Currently powerful computer-aided ab initio techniques can be used to generate predictions of key properties of broad classes of chemicals, without resorting to costly experimentation and potentially hazardous testing. These data could be integrated into a centralized IT decision and compliance support system, and stored in a retrievable, easily communicable manner should new regulatory and/or production requirements necessitate the introduction of different uses of chemicals under different conditions. For illustration purposes, ab initio calculations are performed on heterocyclic nitrogen-containing compounds which currently serve as high energy density materials in the chemical industry. Since investigations of these compounds are still in their infancy, stability studies are imperative regarding their safe handling and storage, as well as registration under REACH.

High Pressure Pyrolysis of Toluene. 2. Modeling Benzyl Decomposition and Formation of Soot Precursors

The pyrolysis of toluene, the simplest methyl-substituted aromatic molecule, has been studied behind reflected shock waves using a single pulse shock tube. Part 1 in this two-part series focused on the high-pressure experimental results and the high-pressure limiting rate coefficients for the primary steps in toluene decomposition. The present work focuses on the modeling of benzyl decomposition and the growth of key soot precursors (C2H2, C4H2, C8H6, and indene) from toluene pyrolysis with 81 among the 262 reactions in the detailed toluene model representing the chemistry that describes the formation and decomposition of these species. Feasible pathways for benzyl decomposition as well as phenylacetylene and indene formation have been tested. The simulations also show very good agreement with the single pulse shock tube profiles for the growth of key soot precursors such as C2H2, C4H2, C8H6, and indene.