How well can B3LYP heats of formation be improved by dispersion correction models?

Thermal stabilities and conformational behaviors of isocyanurates and cyclotrimerization energies of isocyanates: a computational study

Isocyanurates are cyclic trimers of isocyanate molecules. They are generally known as highly thermostable compounds. However, it is interesting how the thermal stabilities of the isocyanurate molecules will be altered depending on the substituents of their three nitrogen atoms. We performed computational investigations on the thermochemical behaviors of isocyanurate molecules with various alkyl and phenyl substituents. The cyclotrimerization processes of isocyanates are highly exothermic. Our best estimate of the enthalpy change for the cyclotrimerization of methyl isocyanate into trimethyl isocyanurate was -66.4 kcal mol-1. Additional negative cyclotrimerization enthalpy changes were observed for n-alkyl-substituted isocyanates. This trend was enhanced with an extension of n-alkyl chains. Conversely, low negative cyclotrimerization enthalpy changes were shown for secondary and tertiary alkyl-substituted isocyanates. The n-alkyl-substituted isocyanurates were shown to be stabilized due to attractive dispersion interactions between the substituents. Meanwhile, the branched alkyl-substituted isocyanurates were destabilized due to the deformation of their isocyanurate rings. For various alkyl-substituted isocyanates, the sum of the deformation energy of the isocyanurate ring and the intramolecular inter-substituent nonbonding interaction energies was found to be linearly correlated with their cyclotrimerization energies. The cyclotrimerization energy for phenyl isocyanate was shown to have significantly deviated from the linear relationship observed for the alkyl-substituted isocyanurates. This is probably attributable to a remarkable change in the orbital resonance interactions during the cyclotrimerization of phenyl isocyanate to triphenyl isocyanurate.

Accurate heats of formation of polycyclic saturated hydrocarbons predicted by using the XYG3 type of doubly hybrid functionals

Polycyclic saturated hydrocarbons (PSHs) are attractive candidates as hydrocarbon propellants. To assess their potential values, one of the key factors is to determine their energy contents, such as to calculate their heats of formation (HOF). In this work, we have calculated HOFs for a set of 36 PSHs including exo-Tricyclo[5.2.1.0^(2,6) ] decane, the principal component of the high-energy density hydrocarbon fuel commonly identified as JP-10. The results from B3LYP, B3LYP-D3BJ, M06-2X, B2PLYP, B2PLYP-D3BJ, and the XYG3 type of doubly hybrid (xDH) functionals are presented. It is demonstrated here that the xDH functionals yield accurate HOFs in good agreement with those from experiments or the G4 theory. In particular, XYGJ-OS, a low scaling xDH functional, is shown to hold the promise for accurate prediction of HOFs for PSHs of larger sizes. © 2018 Wiley Periodicals, Inc.

Benchmarking the Effective Fragment Potential Dispersion Correction on the S22 Test Set

The usual modeling of dispersion interactions in density functional theory (DFT) is often limited by the use of empirically fitted parameters. In this study, the accuracies of the popular empirical dispersion corrections and the first-principles derived effective fragment potential (EFP) dispersion correction are compared by computing the DFT-D and HF-D equilibria interaction energies and intermolecular distances of the S22 test set dimers. Functionals based on the local density approximation (LDA) and generalized gradient approximation (GGA), as well as hybrid functionals, are compared for the DFT-D calculations using coupled cluster CCSD(T) at the complete basis set (CBS) limit as the reference method. In general, the HF-D(EFP) method provides accurate equilibrium dimerization energies and intermolecular distances for hydrogen-bonded systems compared to the CCSD(T)/CBS reference data without using any empirical parameters. For dispersion-dominant and mixed systems, the structures and interaction energies obtained with the B3LYP-D(EFP) method are similar to or better than those obtained with the other DFT-D and HF-D methods. Overall, the first-principles derived -D(EFP) correction presents a robust alternative to the empirical -D corrections when used with the B3LYP functional for dispersion-dominant and mixed systems or with Hartree-Fock for hydrogen-bonded systems.