MN15-L: A New Local Exchange-Correlation Functional for Kohn-Sham Density Functional Theory with Broad Accuracy for Atoms, Molecules, and Solids

Prediction of pKa s of Late Transition-Metal Hydrides via a QM/QM Approach

Three implicit solvation models, the conductor-like polarizable continuum model (C-PCM), the conductor-like screening model (COSMO), and universal implicit solvent model (SMD), combined with a hybrid two layer QM/QM approach (ONIOM), were utilized to calculate the pK_a values, using a direct thermodynamic scheme, of a set of Group 10 transition metal (TM) hydrides in acetonitrile. To obtain the optimal combination of quantum methods for ONIOM calculations with implicit solvation models, the influence of factors, such as the choice of density functional and basis set, the atomic radii used to build a cavity in the solvent, and the size of the model system in an ONIOM scheme, was examined. Additionally, the impact of Grimme's empirical dispersion correction and exact exchange was also investigated. The results were calibrated by experimental data. This investigation provides insight about effective models for the prediction of thermodynamic properties of TM-containing complexes with bulky ligands. © 2019 Wiley Periodicals, Inc.

Effect of ammonia and formic acid on the CH3O˙ + O2 reaction: a quantum chemical investigation

In the present work, the catalytic effect of ammonia and formic acid on the CH3O˙ + O2 reaction has been investigated employing the MN15L density functional. The investigations suggest that, in the presence of ammonia, the reaction can proceed through two different pathways, namely a single hydrogen atom transfer and a double hydrogen atom transfer path, but due to the high energy barrier associated with the double hydrogen atom transfer channel, it prefers the single hydrogen atom transfer channel. On the other hand, in the case of formic acid, only the single hydrogen atom transfer path is found to be feasible. Interestingly, it has been found that, in the presence of ammonia and formic acid, the reaction becomes a barrierless reaction. The calculated rate constant values at various temperatures indicate an anti-Arrhenius behavior for both the ammonia and formic acid catalyzed channels.

Anharmonic kinetics of the cyclopentane reaction with hydroxyl radical

Cyclopentane is one of the major constituents of transportation fuels, especially jet fuel and diesel, and also is a volatile organic compound with a significant presence in the atmosphere. Hydrogen abstraction from cyclopentane by hydroxyl radical plays a significant role in combustion and atmospheric chemistry. In this work we study the kinetics of this reaction at 200-2000 K using direct dynamics calculations in which the potential energy surface is obtained by quantum mechanical electronic structure calculations. The forward and reverse barrier heights and reaction energies obtained by the CCSD(T)-F12a/jun-cc-pVTZ coupled cluster calculations are used as a benchmark to select an accurate electronic structure method among 36 combinations of exchange-correlation functional and basis set. The selected M06-2X/MG3S method shows the best performance with a mean unsigned deviation from the benchmark of only 0.22 kcal mol-1 for reaction energies and barrier heights. A quadratic-quartic function is adopted to describe the ring bending potential of cyclopentane, and the quartic anharmonicity in the bending mode is treated by a one-dimensional Schrödinger equation using both Wentzel-Kramers-Brillouin (WKB) and Fourier Grid Hamiltonian (FGH) methods. The torsional anharmonicity in the transition state is treated in turn by the free rotor approximation, the one-dimensional hindered rotor approximation, and the multi-structural torsional anharmonicity method. Rate constants of the title reaction are computed by canonical variational transition state theory including tunneling by the multi-dimensional small-curvature tunneling approximation (CVT/SCT). The final rate constants include the quasiharmonic, quadratic-quartic, and torsional anharmonicity. Our calculations are in excellent agreement with all the experimental data available at both combustion and atmospheric temperatures with a deviation of less than 30%.

A molecular electron density theory (MEDT) study of the role of halogens (X2 = F2, Cl2, Br2 and I2) on the aza-Michael-addition reactions

The energy profile for gas phase F₂-catalyzed nucleophilic attack of PYR R1 on the β-conjugated carbon of MA R2 marked.

The one-electron self-interaction error in 74 density functional approximations: a case study on hydrogenic mono- and dinuclear systems

The self-interaction error (SIE), i.e. unphysical interactions of electrons with themselves, has plagued developers and users of Density Functional Approximations (DFAs) since the inception of Density Functional Theory (DFT). Formally, it can be separated into the one-electron and many-electron SIE; herein we present one of the most comprehensive studies of the first. While we focus mostly on the total SIE, we also make use of two different decompositions. The first is a separation into functional and density-driven errors as championed by Sim, Burke and co-workers [J. Phys. Chem. Lett., 2018, 9, 6385-6392]; the second separates the error into exchange, correlation, and one-electron components, with the latter being a density error that has not been discussed in this form before. After investigating the familiar hydrogen atom and dihydrogen cation, we establish a relationship between the SIE and the nuclear charge with the help of a series of heavier hydrogenic analogues. For the mononuclear systems and the diatomics at the dissociation limit, this relationship is linear in nature with prominent exceptions, mostly belonging to the Minnesota and range-separated (double-)hybrid DFAs. For the first time, we also show how the magnitude of the SIE depends on the underlying atomic-orbital basis set and how DFAs that rely on a popular van-der-Waals DFT type London-dispersion term exhibit "self-dispersion". We find that range separation is not a panacea for solving the one-electron SIE. DFAs that have been developed to be one-electron SIE free for one system, such as the hydrogen atom, show larger errors for heavier hydrogenic systems. Often, one-electron SIE-free DFAs rely on fortuitous error cancellation between their exchange and correlation components. An analysis of the most robust methods for general applications to date reveals that they suffer moderately from the one-electron SIE, while DFAs that are nearly SIE-free do not perform well in applications. Implicit in the continued existence of the one-electron SIE is that well-performing DFAs continue to suffer insufficiencies at their fundamental levels that are being compensated for by the SIE. Our analysis includes more than 250 000 datapoints, resulting in multiple insights that may drive future developments of new DFAs or SIE corrections.

Polymorphism in a π stacked Blatter radical: structures and magnetic properties of 3-(phenyl)-1-(pyrid-2-yl)-1,4-dihydrobenzo[e][1,2,4]triazin-4-yl

The first polymorphism example in Blatter radicals, is reported.

Fundamental roles of ZnO and ZrO2 in the conversion of ethanol to 1,3-butadiene over ZnO–ZrO2/SiO2

ZnO and ZrO₂ sites in the ZnO–ZrO₂/SiO₂ catalyst act cooperatively to catalyze different steps of the conversion of ethanol to 1,3-butadiene.

Sulfur hydrogen bonding and internal dynamics in the monohydrates of thenyl mercaptan and thenyl alcohol

Water forms weak H-bonds with thenyl compounds, simultaneously retaining internal mobility in the dimer.

Unraveling the kinetics and molecular mechanism of gas phase pyrolysis of cubane to [8]annulene

The kinetic and electron density flows are studied theoretically for the gas phase pyrolysis of cubane via its cage opening to reach bicyclooctatriene and then thermal rearrangement of bicyclooctatriene to produce [8]annulene which is the experimentally observed major product. The observed kinetic data at the MN15-L/maug-cc-pVTZ level of theory were in good agreement with the experimental results as compared to the CBS-QB3 method. The cage opening and the thermal rearrangement steps at the experimentally employed temperature of 520 K were exergonic and exothermic. The atmospheric rate constants calculated by means of the RRKM theory show that the cage opening is the rate-determining step. The temperature dependence of the rate constant for the cage opening step at the MN15-L level can be expressed as log(k/s⁻¹)^1bar_MN15-L = (15.63) − (48.99 kcal mol⁻¹)/RT ln 10. The molecular mechanism of the reactions has been investigated by means of the bonding evolution theory (BET) at the B3LYP/6-311G (d,p) level of theory. The cage opening course is described topologically by cleaving of C1–C2, C4–C8, and C5–C6 single bonds and electron saturation of the C1–C4, C2–C6, and C5–C8 bonds, while the rearrangement of bicyclooctatriene is described by C3–C7 bond rupture, depopulation of C1–C4 and C5–C8 double bonds, and electron saturation of C1–C5, C3–C4, and C7–C8 bonds. Electron density rearrangement along the two successive steps are asynchronous and the sequence of catastrophes can be represented as: η-1-13-C^†C^†FFFC^†C^†FFFC^†C^†-2-6-[C]₂C^†[F]₂[C^†]₂C^†-0.

Kinetic and electron density flows are studied theoretically for the gas phase pyrolysis of cubane via its cage opening to bicyclooctatriene and then thermal rearrangement of bicyclooctatriene to produce [8]annulene which is the experimentally observed major product.

Multistructural Anharmonicity Controls the Radical Generation Process in Biofuel Combustion

The OH radical plays an important role in combustion, and isopentanol (3-methylbutan-1-ol) is a promising sustainable fuel additive and second-generation biofuel. The abstractions of H atoms from fuel molecules are key initiation steps for chain branching in combustion chemistry. In comparison with the more frequently studied ethanol, isopentanol has a longer carbon chain that allows a greater number of products, and experimental work is unavailable for the branching fractions to the various products. However, the site-dependent kinetics of isopentanol with OH radicals are usually experimentally unavailable. Alcohol oxidation by OH is also important in the atmosphere, and in the present study we calculate the rate constants and branching fractions of the hydrogen abstraction reaction of isopentanol by OH radical in a broad temperature range of 298-2400 K, covering temperatures important for atmospheric chemistry and those important for combustion. The calculations are done by multipath variational transition state theory (MP-VTST). With a combination of electronic structure calculations, we determine previously missing thermochemical data. With MP-VTST, a multidimensional tunneling approximation, multiple-structure anharmonicity, and torsional potential anharmonicity, we carried out more realistic rate constant calculations than can be computed by conventional single-structure harmonic transition state theory or by the empirical relations that are currently used in atmospheric and combustion modeling. The roles of various factors in determining the rates are elucidated, and we show that recrossing, tunneling, and multiple structures are all essential for accurate work. We conclude that the multiple structure anharmonicity is the most important correction to conventional transition state theory for this reaction, although recrossing effects and tunneling are by no means insignificant and the tunneling depends significantly on the path. The thermodynamic and kinetics data determined in this work are indispensable for the gas-phase degradation of alcohols in the atmosphere and for the detailed understanding and prediction of ignition mechanisms of biofuels in combustion.

Machine Learning the Physical Nonlocal Exchange-Correlation Functional of Density-Functional Theory

We train a neural network as the universal exchange-correlation functional of density-functional theory that simultaneously reproduces both the exact exchange-correlation energy and the potential. This functional is extremely nonlocal but retains the computational scaling of traditional local or semilocal approximations. It therefore holds the promise of solving some of the delocalization problems that plague density-functional theory, while maintaining the computational efficiency that characterizes the Kohn-Sham equations. Furthermore, by using automatic differentiation, a capability present in modern machine-learning frameworks, we impose the exact mathematical relation between the exchange-correlation energy and the potential, leading to a fully consistent method. We demonstrate the feasibility of our approach by looking at one-dimensional systems with two strongly correlated electrons, where density-functional methods are known to fail, and investigate the behavior and performance of our functional by varying the degree of nonlocality.

Kohn-Sham Density Functional Theory with Complex, Spin-Restricted Orbitals: Accessing a New Class of Densities without the Symmetry Dilemma

We show that using complex, spin-restricted orbitals in Kohn-Sham (KS) density functional theory allows one to access a new class of densities that is not accessible by either spin-restricted (RKS) or spin-unrestricted (UKS) orbitals. We further show that the real part of a complex RKS (CRKS) density matrix can be nonidempotent when the imaginary part of the density matrix is not zero. Using CRKS orbitals shows significant improvements in the triplet-singlet gaps of a benchmark set, called TS12, for well-established, widely used density functionals. Moreover, it was shown that RKS and UKS yield qualitatively wrong charge densities and spin densities, respectively, leading to worse energetics. We demonstrate that representative modern density functionals show surprisingly no improvement even with a qualitatively more accurate density from CRKS orbitals. To this end, our work not only provides a way to escape the symmetry dilemma whenever there exists a CRKS solution, but also suggests a new route to design better approximate density functionals.

Cost-effective density functional theory (DFT) calculations of equilibrium isotopic fractionation in large organic molecules

The application of stable isotopes to address a wide range of biochemical, microbiological and environmental problems is hindered by the experimental difficulty and the computational cost of determining equilibrium isotopic fractionations (EIF) of large organic molecules. Here, we evaluate the factors that impact the accuracy of computed EIFs and develop a framework for cost-effective and accurate computation of EIFs by density functional theory (DFT). We generated two benchmark databases of experimentally determined EIFs, one for H isotopes and another for the isotopes of the heavy atoms C, N and O. The accuracy of several DFT exchange-correlation functionals in calculating EIFs was then evaluated by comparing the computational results to these experimental datasets. We find that with the def2-TZVP basis set, O3LYP had the lowest mean absolute deviation (21‰ and 3.9‰ for the isotopic fractionation of H and the heavier atoms, respectively), but the GGA/meta-GGA functionals τHCTH_D3BJ, τHCTH and HCTH have similar performances (22‰ and 4.1‰, respectively, for τHCTH_D3BJ). Leveraging the good performance of computationally efficient functionals, we provide a robust, practical, experimentally validated framework for using DFT to accurately predict EIFs of large organic molecules, including uncertainty estimates.

Influence of water on the CH3O˙ + O2 → CH2O + HO2˙ reaction

Electronic structure calculations employing density functional theory have been used to study the effect of a single water molecule on the CH3O˙ + O2 → CH2O + HO2˙ reaction. The investigation suggests that in the presence of water the reaction barrier reduces from 3.01 kcal mol-1 to -1.86 kcal mol-1. Consequently, when we consider the bimolecular rate constants for the water catalyzed channel, they were found to be 104 to 105 times higher than that of the uncatalyzed reaction. Interestingly, the Arrhenius plot indicates a negative temperature dependency of the catalyzed channel (anti-Arrhenius behavior); as a result of this the domination of the catalyzed channel over the bare reaction increases with the lowering of the temperature. But the effective bimolecular rate constant values for the catalyzed channel were found to be approximately four orders of magnitude lower than that of the uncatalyzed one, which implies that the contribution of the catalyzed channels to the overall rate of the reaction is very small.

The CUAGAU Set of Coupled-Cluster Reference Data for Small Copper, Silver, and Gold Compounds and Assessment of DFT Methods

We have obtained benchmark data for a set of small molecular systems of Cu, Ag, and Au using coupled-cluster methods. Using this collection of reference data (that we termed the CUAGAU set) for assessing DFT-type methods, we find the MN15-L nonhybrid DFT to be cost-effective for geometry optimization [mean absolute deviation (MAD) in bond length = 0.20 Å], with an accuracy that is comparable to that for the double-hybrid (DH) DFT method DSD-PBEP86 (MAD = 0.19 Å). For the computation of thermochemical properties, among "conventional" (non-MP2-based) DFT methods, the best performance is found for the global-hybrid meta-GGA functional MN15, with an MAD of 11.4 kJ mol^-1. We also find the nonhybrid method B97M-rV to have a reasonable performance (MAD = 14.4 kJ mol^-1), and it may serve as a cost-effective means for qualitative study. If we look beyond conventional functionals, we find DSD-PBEP86 (MAD = 7.3 kJ mol^-1) to be more accurate than even MN15. Nonetheless, this level of accuracy is still not sufficient for quantitative studies. In this regard, high-level wave function methods such as composite procedures that are based on coupled cluster are still indispensable for obtaining reliable reference data for transition-metal species.

Prediction of Excited-State Properties of Oligoacene Crystals Using Fragment-Based Quantum Mechanical Method

A fundamental understanding of the excited-state properties of molecular crystals is of central importance for their optoelectronics applications. In this study, we developed the electrostatically embedded generalized molecular fractionation (EE-GMF) method for the quantitative characterization of the excited-state properties of locally excited molecular clusters. The accuracy of the EE-GMF method is systematically assessed for oligoacene crystals. Our result demonstrates that the EE-GMF method is capable of providing the lowest vertical singlet (S₁) and triplet excitation energies (T₁), in excellent agreement with the full-system quantum mechanical calculations. Using this method, we also investigated the performance of different density functionals in predicting the excited-state properties of the oligoacene crystals. Among the 13 tested functionals, B3LYP and MN15 give the two lowest overall mean unsigned errors with reference to the experimental S₁ and T₁ excitation energies. The EE-GMF approach can be readily utilized for studying the excited-state properties of large-scale organic solids at diverse ab initio levels.

Rapid unimolecular reaction of stabilized Criegee intermediates and implications for atmospheric chemistry

Elucidating atmospheric oxidation mechanisms is necessary for estimating the lifetimes of atmospheric species and understanding secondary organic aerosol formation and atmospheric oxidation capacity. We report an unexpectedly fast mechanistic pathway for the unimolecular reactions of large stabilized Criegee intermediates, which involves the formation of bicyclic structures from large Criegee intermediates containing an aldehyde group. The barrier heights of the mechanistic pathways are unexpectedly low – about 2–3 kcal/mol – and are at least 10 kcal/mol lower than those of hydrogen shift processes in large syn Criegee intermediates; and the calculated rate constants show that the mechanistic pathways are 10⁵-10⁹ times faster than those of the corresponding hydrogen shift processes. The present findings indicate that analogous low-energy pathways can now also be expected in other large Criegee intermediates and that oxidative capacity of some Criegee intermediates is smaller than would be predicted by existing models.

Criegee intermediates have received much attention in atmospheric chemistry because of their importance in ozonolysis mechanisms. Here, using quantum mechanical kinetics, the authors reveal an unexpectedly fast mechanistic pathway for unimolecular reactions of large stabilized Criegee intermediates.

Synthesis, characterization and DNA binding studies of a new ibuprofen-platinum(II) complex

The study is focused on the synthesis of a novel complex of ibuprofen and Platinum(II). The formation of the product was characterized through analytical tools including Fourier-transform infrared spectroscopy, proton nuclear magnetic resonance, ultraviolet-visible spectroscopy, mass spectrometry as well as density functional theory. Using the continuous variation method, the stoichiometry of Pt(HIb)₂(Cl)₂ binding on DNA (ct-DNA) determines as a single class of binding. Based on the results of Stern-Volmer analysis on the fluorescence quenching data, the quenching mechanism was determined to be static in nature. The studies indicated that the complex could bind DNA molecules via groove binding for four major reasons. Initially, the complex-DNA binding constant determined based on spectrophotometric data were found to be comparable to those of groove-binding drugs. In addition, the competitive fluorimetric based on the applications of Hoechst 33258 proved the ability of Pt(HIb)₂(Cl)₂ to displace with Hoechst in its DNA-bounded form, reflecting the competition between Pt(HIb)₂(Cl)₂ and Hoechst for groove binding. Further, no considerable changes were observed in the intensity of the methylene blue (MB)-DNA system after adding the Pt(HIb)₂(Cl)₂ complex, reflecting the stability of MB molecules in the DNA helix and a non-intercalative bonds of Pt(HIb)₂(Cl)₂ interaction on DNA. Finally, minor changes in the viscosity of DNA in the presence of Pt(HIb)₂(Cl)₂, indicated that changes in the length of DNA in the presence of the complex are negligible, supporting the assumption of DNA groove-binding. Also induced CD spectral changes and docking simulations were in favor of the groove mechanism for the Pt(HIb)₂(Cl)₂-DNA binding.Communicated by Ramaswamy H. Sarma.

Evaluating Transition Metal Barrier Heights with the Latest Density Functional Theory Exchange-Correlation Functionals: The MOBH35 Benchmark Database

A new database of transition metal reaction barrier heights (MOBH35) is presented. Benchmark energies (forward and reverse barriers and reaction energy) are calculated using DLPNO-CCSD(T) extrapolated to the complete basis set limit using a Weizmann-1-like scheme. Using these benchmark energies, the performance of a wide selection of density functional theory (DFT) exchange-correlation functionals, including the latest from the Martin, Truhlar, and Head-Gordon groups, is evaluated. It was found, using the def2-TZVPP basis set, that the ωB97M-V (MAD 1.7 kcal/mol), ωB97M-D3BJ (MAD 1.9 kcal/mol), ωB97X-V (MAD 2.0 kcal/mol), and revTPSS0-D4 (MAD 2.2 kcal/mol) hybrid functionals are recommended. The double-hybrid functionals B2K-PLYP (MAD 1.7 kcal/mol) and revDOD-PBEP86-D4 (MAD 1.8 kcal/mol) also performed well, but this has to be balanced by their increased computational cost.

Benchmarking quantum chemistry methods for spin-state energetics of iron complexes against quantitative experimental data

The accuracy of relative spin-state energetics predicted by selected quantum chemistry methods: coupled cluster theory at the CCSD(T) level, multiconfigurational perturbation theory (CASPT2, NEVPT2), multireference configuration interaction at the MRCISD+Q level, and a number of DFT methods, is quantitatively evaluated by comparison with the experimental data of four octahedral iron complexes. The available experimental data, either spin-forbidden transition energies or spin crossover enthalpies, are corrected for relevant environmental effects in order to derive the quantitative benchmark set of iron spin-state energetics. Comparison of theory predictions with the resulting reference data: (1) validates the high accuracy of the CCSD(T) method, particularly when based on Kohn-Sham orbitals, giving the maximum error below 2 kcal mol-1 and the mean absolute error (MAE) below 1 kcal mol-1; (2) corroborates the tendency of CASPT2 to systematically overstabilize higher-spin states by up to 5.5 kcal mol-1; (3) confirms that the latter problem is partly remedied by the recently proposed CASPT2/CC approach [Phung et al., J. Chem. Theory Comput., 2018, 14, 2446-2455]; (4) demonstrates that NEVPT2 performs worse than CASPT2, by giving errors up to 7 kcal mol-1; (5) shows that the accuracy of MRCISD+Q spin-state energetics strongly depends on the size-consistency correction: the Davidson-Silver and Pople corrections perform best (MAE < 3 kcal mol-1), whereas the standard Davidson correction is not recommended (MAE of 7 kcal mol-1). Only a few DFT methods (including the best performing ones identified in this study: B2PLYP-D3 and OPBE) are able to provide a balanced description of the spin-state energetics for all four studied iron complexes simultaneously, corroborating the non-universality problem of approximate density functionals.

Benchmark study of popular density functionals for calculating binding energies of three-center two-electron bonds

Density functional theory (DFT) can be used to study the three-center two-electron (3c2e) bonding mode, which is universal in catalysts containing alkaline-earth (Ae) and boron-group (Bg) elements. However, because of the delocalization pattern of the 3c2e bond, the wavefunction cannot be accurately described by DFT methods. The calculated energies of Ae and Bg catalysts therefore fluctuate greatly when different functionals are used, largely because of inconsistent DFT-calculated binding energies of 3c2e bonds. Nevertheless, with the development of supercomputers and theoretical calculation software, the DFT method is becoming increasingly popular for studying Ae and Bg catalysts. In this study, we compared the performances of 21 functionals with the high-level composite G3B3 method in calculations for the binding energies of 3c2e bonds. Several frequently used post-Hartree-Fock methods were also tested. The calculation results indicate that the M06-2X, MN12-L, and MN15 functionals give consistent and reliable binding energies for common 3c2e bonds. © 2018 Wiley Periodicals, Inc.

Ferromagnetic interactions in a 1D Heisenberg linear chain of 1-phenyl-3,7-bis(trifluoromethyl)-1,4-dihydro-1,2,4-benzotriazin-4-yls

1-Phenyl-3,7-bis(trifluoromethyl)-1,4-dihydro-1,2,4-benzotriazin-4-yl (2), was characterized by single crystal X-ray diffractometry and variable temperature SQUID magnetometry to investigate its structure-magnetism correlation.

Statistically representative databases for density functional theory via data science

Cluster analysis applied to quantum chemistry: a new broad database of chemical properties with a reasonable computational cost.

The effects of active site and support on hydrogen elimination over transition-metal-functionalized yttria-decorated metal–organic frameworks

Transition-metal catalysts supported on a metal–organic framework have been screened computationally to reveal the best catalytic candidates for hydrogen elimination reactions, which are critical in many catalytic cycles.

Screened hybrid meta-GGA exchange–correlation functionals for extended systems

Screened Hartree–Fock exchange integrated with semilocal exchange–correlation functionals often proficiently predict several solid-state properties.

Computational kinetics of the hydrogen abstraction reactions of n-propanol and iso-propanol by OH radical

The total reaction rate constants show a significant negative dependence on temperature in the low temperature regime and approach the capture rate for the formation of the pre-reactive complex when temperature is down to the ultracold regime.

A Trip to the Density Functional Theory Zoo: Warnings and Recommendations for the User

This account is written for general users of density functional theory (DFT) methods as well as experimental researchers who are new to the field and would like to conduct such calculations. Its main emphasis lies on how to find a way through the confusing ‘zoo’ of DFT by addressing common misconceptions and highlighting those modern methods that should ideally be used in calculations of energetic properties and geometries. A particular focus is on highly popular methods and the important fact that popularity does not imply accuracy. In this context, we present a new analysis of the openly available data published in Swart and co-workers’ famous annual ‘DFT poll’ (http://www.marcelswart.eu/dft-poll/) to demonstrate the existing communication gap between the DFT user and developer communities. We show that despite considerable methodological advances in the field, the perception of some parts of the user community regarding their favourite approaches has changed little. It is hoped that this account makes a contribution towards changing this status and that users are inspired to adjust their current computational protocols to accommodate strategies that are based on proven robustness, accuracy, and efficiency rather than popularity.

Performance of density functional theory for describing hetero-metallic active-site motifs for methane-to-methanol conversion in metal-exchanged zeolites

Methane-to-methanol conversion (MMC) can be facilitated with high methanol selectivities by copper-exchanged zeolites. There are however two open questions regarding the use of these zeolites to facilitate the MMC process. The first concerns the possibility of operating the three cycles in the stepwise MMC process by these zeolites in an isothermal fashion. The second concerns the possibility of improving the methanol yields by systematic substitution of some copper centers in these active sites with other earth-abundant transition metals. Quantum-mechanical computations can be used to compare methane activation by copper oxide species and analogous mixed-metal systems. To carry out such screening, it is important that we use theoretical methods that are accurate and computationally affordable for describing the properties of the hetero-metallic catalytic species. We have examined the performance of 47 exchange-correlation density functionals for predicting the relative spin-state energies and chemical reactivities of six hetero-metallic [M-O-Cu]²⁺ and [M-O₂ -Cu]²⁺ , (where MCo, Fe, and Ni), species by comparison with coupled cluster theory including iterative single, double excitations as well as perturbative treatment of triple excitations, CCSD(T). We also performed multireference calculations on some of these systems. We considered two types of reactions (hydrogen addition and oxygen addition) that are relevant to MMC. We recommend the use of τ-HCTH and OLYP to determine the spin-state energy splittings in the hetero-metallic motifs. ωB97, ωB97X, ωB97X-D3, and MN15 performed best for predicting the energies of the hydrogen and oxygen addition reactions. In contrast, local, and semilocal functionals do poorly for chemical reactivity. Using [Fe-O-Cu]²⁺ as a test, we see that the nonlocal functionals perform well for the methane CH activation barrier. In contrast, the semilocal functionals perform rather poorly. © 2018 Wiley Periodicals, Inc.

Bond dissociation energies of FeSi, RuSi, OsSi, CoSi, RhSi, IrSi, NiSi, and PtSi

Resonant two-photon ionization spectroscopy has been used to investigate the spectra of the diatomic late transition metal silicides, MSi, M = Fe, Ru, Os, Co, Rh, Ir, Ni, and Pt, in the vicinity of the bond dissociation energy. In these molecules, the density of vibronic states is so large that the spectra appear quasicontinuous in this energy range. When the excitation energy exceeds the ground separated atom limit, however, a new decay process becomes available-molecular dissociation. This occurs so rapidly that the molecule falls apart before it can absorb another photon and be ionized. The result is a sharp drop to the baseline in the ion signal, which we identify as occurring at the thermochemical 0 K bond dissociation energy, D₀. On this basis, the measured predissociation thresholds provide D₀ = 2.402(3), 4.132(3), 4.516(3), 2.862(3), 4.169(3), 4.952(3), 3.324(3), and 5.325(9) eV for FeSi, RuSi, OsSi, CoSi, RhSi, IrSi, NiSi, and PtSi, respectively. Using thermochemical cycles, the enthalpies of formation of the gaseous MSi molecules are derived as 627(8), 700(10), 799(10), 595(8), 599(8), 636(10), 553(12), and 497(8) kJ/mol for FeSi, RuSi, OsSi, CoSi, RhSi, IrSi, NiSi, and PtSi, respectively. Likewise, combining these results with other data provides the ionization energies of CoSi and NiSi as 7.49(7) and 7.62(7) eV, respectively. Chemical bonding trends among the diatomic transition metal silicides are discussed.

Delocalization Errors in Density Functional Theory Are Essentially Quadratic in Fractional Occupation Number

Approximate functionals used in practical density functional theory (DFT) deviate from the piecewise linear behavior of the exact functional for fractional charges. This deviation causes excess charge delocalization, which leads to incorrect densities, molecular properties, barrier heights, band gaps, and excitation energies. We present a simple delocalization function for characterizing this error and find it to be almost perfectly linear vs the fractional electron number for systems spanning in size from the H atom to the C₁₂H₁₄ polyene. This causes the delocalization energy error to be a quadratic polynomial in the fractional electron number, which permits us to assess the comparative performance of 47 popular and recent functionals through the curvature. The quadratic form further suggests that information about a single fractional charge is sufficient to eliminate the principal source of delocalization error. Generalizing traditional two-point information like ionization potentials or electron affinities to account for a third, fractional charge-based data point could therefore permit fitting/tuning of functionals with lower delocalization error.