Gas-phase thermochemistry of polycyclic aromatic hydrocarbons: an approach integrating the quantum chemistry composite scheme and reaction generator

We introduce a protocol aimed at predicting the accurate gas-phase enthalpies of formation of polycyclic aromatic hydrocarbons (PAHs). Automatic generation of a dataset of equilibrated chemical reactions preserving the number of carbon atoms in each hybridization state on each side of equations is at the core of our scheme. The performed tests suggest the recommended enthalpy of formation to be derived via a two-step scheme. First, we consider the reactions with a minimal sum of the total number of particles involved, N, and the absolute difference between the total number of products and reactants, |ΔN|. Second, among these reactions, we identify the one with the smallest absolute reaction enthalpy change, . This approach has been applied to predict the gas-phase enthalpies of formation of 113 PAHs via the Feller-Peterson-Dixon approach. Our calculated values provide the mean absolute deviations of 1.7, 1.9, 4.2, 8.1, and 18.5 kJ mol^-1 with respect to the literature group-based error corrected (GBEC) G3MP2B3, ATOMIC (HC), group equivalent M06-2X, GBEC B3LYP, and G4MP2 values. Our predicted values give the mean signed and mean absolute errors of -7.5 and 12.9 kJ mol^-1 with respect to the experimental enthalpies of formation. The combination of our predicted and the experimental values provide the solid-state enthalpies of formation, , which are not available for a few species. Approaching these values as well as , producing large discrepancies from the experimental side, would be indispensable for testing and further tuning of computational chemistry approaches.

Learning to fly: thermochemistry of energetic materials by modified thermogravimetric analysis and highly accurate quantum chemical calculations

The standard state enthalpy of formation and the enthalpy of sublimation are essential thermochemical parameters determining the performance and application prospects of energetic materials (EM). Direct experimental measurements of these properties are complicated by low volatility and high heat release in bomb calorimetry experiments. As a result, the uncertainties in the reported enthalpies of formation for a number of even well-known CHNO-containing compounds might amount up to tens kJ mol-1, while for some novel high-nitrogen molecules they reach even hundreds of kJ mol-1. The present study reports a facile approach to determining the solid-state formation enthalpies comprised of complementary high-level quantum chemical calculations of the gas-phase thermochemistry and advanced thermal analysis techniques yielding sublimation enthalpies. The thermogravimetric procedure for the measurement of sublimation enthalpy was modified by using low external pressures (down to 0.2 Pa). This allows for observing sublimation/vaporization instead of thermal decomposition of the compounds studied. Extensive benchmarking on nonenergetic and energetic compounds reveals the average and maximal absolute errors of the sublimation enthalpies of 3.3 and 11.0 kJ mol-1, respectively. The comparison of the results with those obtained from the widely used Trouton-Williams empirical equation shows that the latter underestimates the sublimation enthalpy up to 140 kJ mol-1. Therefore, we performed a reparametrization of the latter equation with simple chemical descriptors that reduces the mean error down to 30 kJ mol-1. Highly accurate multi-level procedures W2-F12 and/or W1-F12 in conjunction with the atomization energy approach were used to calculate theoretically the gas-phase formation enthalpies. In several cases, the DLPNO-CCSD(T) enthalpies of isodesmic reactions were also employed to obtain the gas-phase thermochemistry for medium-sized important EMs. Combining the obtained thermochemical properties, we determined the solid-state enthalpies of formation for nearly 60 species containing various important explosophoric groups, from common nitroaromatics, nitroethers, and nitramines to novel nitrogen-rich heterocyclic species (e.g., the derivatives of pyrazole, tetrazole, furoxan, etc.). The large-scale benchmarking against the available experimental solid-state enthalpies of formation yielded the maximal inaccuracy of the proposed method of 25 kJ mol-1.

New Biosourced Flame Retardant Agents Based on Gallic and Ellagic Acids for Epoxy Resins

The aim of this work was an investigation of the ability of gallic (GA) and ellagic (EA) acids, which are phenolic compounds encountered in various plants, to act as flame retardants (FRs) for epoxy resins. In order to improve their fireproofing properties, GA and EA were treated with boric acid (to obtain gallic acid derivatives (GAD) and ellagic acid derivatives (EAD)) to introduce borate ester moieties. Thermogravimetric analysis (TGA) highlighted the good charring ability of GA and EA, which was enhanced by boration. The grafting of borate groups was also shown to increase the thermal stability of GA and EA that goes up respectively from 269 to 528 °C and from 496 to 628 °C. The phenolic-based components were then incorporated into an epoxy resin formulated from diglycidyl ether of bisphenol A (DGEBA) and isophorone diamine (IPDA) (72, 18, and 10 wt.% of DGEBA, IPDA, and GA or EA, respectively). According to differential scanning calorimetry (DSC), the glass transition temperature (T_g) of the thermosets was decreased. Its values ranged from 137 up to 108 °C after adding the phenolic-based components. A cone calorimeter was used to evaluate the burning behavior of the formulated thermosets. A significant reduction of the peak of heat release rate (pHRR) for combustion was detected. Indeed, with 10 wt.% of GA and EA, pHRR was reduced by 12 and 44%, respectively, compared to that for neat epoxy resin. GAD and EAD also induced the decrease of pHRR values by 65 and 33%, respectively. In addition, a barrier effect was observed for the resin containing GAD. These results show the important influence of the biobased phenolic compounds and their boron derivatives on the fire behavior of a partially biobased epoxy resin.

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.

Efficient DLPNO–CCSD(T)-Based Estimation of Formation Enthalpies for C-, H-, O-, and N-Containing Closed-Shell Compounds Validated Against Critically Evaluated Experimental Data

An accurate and cost-efficient methodology for the estimation of the enthalpies of formation for closed-shell compounds composed of C, H, O, and N atoms is presented and validated against critically evaluated experimental data. The computational efficiency is achieved through the use of the resolution-of-identity (RI) and domain-based local pair-natural orbital coupled cluster (DLPNO-CCSD(T)) approximations, which results in a drastic reduction in both the computational cost and the number of necessary steps for a composite quantum chemical method. The expanded uncertainty for the proposed methodology evaluated using a data set of 45 thoroughly vetted experimental values for molecules containing up to 12 heavy atoms is about 3 kJ·mol-1, competitive with those of typical calorimetric measurements. For the compounds within the stated scope, the methodology is shown to be superior to a representative, more general, and widely used composite quantum chemical method, G4.

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.

Heat of adsorption of naphthalene on Pt(111) measured by adsorption calorimetry

The heat of adsorption of naphthalene on Pt(111) at 300 K was measured with single-crystal adsorption calorimetry. The heat of adsorption on the ideal, defect-free surface is estimated to be (300 - 34 - 199(2)) kJ/mol. From this, a C-Pt bond energy for aromatic hydrocarbons on Pt(111) of approximately 30 kJ/mol is estimated, consistent with earlier results for benzene on Pt(111). There is higher heat of adsorption at very low coverage, attributed to step sites where the adsorption heat is >/=330 kJ/mol. Saturation coverage, = 1 ML, corresponds to 1.55 x 10(14) molecules/cm(2). Sticking probability measurements of naphthalene on Pt(111) give a high initial value of 1.0 and a Kisliuk-type coverage dependence that implies precursor-mediated sticking. The ratio of the hopping rate to the desorption rate of this precursor is approximately 51. Naphthalene adsorbs transiently on top of chemisorbed naphthalene molecules with a heat of adsorption of 83-87 kJ/mol.

Polycyclic benzenoids: why kinked is more stable than straight

The enhanced stability of bent or kinked polycyclic benzenoids over linear ones is well established, phenanthrene and anthracene being archetypal representatives. The question why kinked is more stable than linear is, however, still a matter of discussion. Recently, it has been proposed that H-H bonding interactions between the two hydrogen atoms in the bay region of phenanthrene are responsible for the larger stability of this molecule as compared to anthracene. This conclusion conflicts with the vast body of evidence for nonbonded steric repulsion between these hydrogen atoms. In this work, we provide new, complementary evidence for the repulsive character of the H-H interactions in phenanthrene's bay region. We have traced the origin of phenanthrene's enhanced stability to the more efficient bonding in the pi-electron system using, among others, a quantitative energy decomposition analysis of the bonding between the two constituting 2-methtriyl-phenyl fragments in both phenanthrene and anthracene (i.e., C14H10 = C6H4*-CH** + C6H4*-CH**). The scope of our study is extended to polycyclic benzenoids by analyzing also hexacene and various bent isomers of the latter. Our results once more falsify one of the core concepts of the theory of atoms-in-molecules (AIM), namely, that the presence of bond paths and the presence of bond critical points (they exist indeed between the two bay H atoms in phenanthrene) are sufficient indicators for a stabilizing interaction. Instead, our results confirm that these AIM parameters merely diagnose the proximity or contact between charge distributions, be this contact stabilizing or destabilizing.