The Metal Hydride Problem of Computational Chemistry: Origins and Consequences

Computed versus experimental energy barriers in solution: Influence of the type of the density functional approximation

Mechanistic investigations at the density functional theory level of organic and organometallic reactions in solution are now broadly accessible and routinely implemented to complement experimental investigations. The selection of an appropriate functional among the plethora of developed ones is the first challenge on the way to reliable energy barrier calculations. To provide guidelines for the choice of an initial and reliable computational level, the performances of commonly used non-empirical (PBE, PBE0, PBE0-DH) and empirical density functionals (BLYP, B3LYP, B2PLYP) were evaluated relative to experimental activation enthalpies. Most reactivity databases to assess density functional performances are primarily based on high level calculations, here a set of experimental activation enthalpies of organic and organometallic reactions performed in solution were selected from the literature. As a general trend, the non-empirical functionals outperform the empirical ones. The most accurate energy barriers are obtained with hybrid PBE0 and double-hybrid PBE0-DH density functionals, both providing similar performance. Regardless of the functional under consideration, the addition of the GD3-BJ empirical dispersion correction does not enhance the accuracy of computed energy barriers.

Automated, Consistent, and Even-Handed Selection of Active Orbital Spaces for Quantum Embedding

A widely used strategy to reduce the computational cost of quantum-chemical calculations is to partition the system into an active subsystem, which is the focus of the computational efforts, and an environment that is treated at a lower computational level. The system partitioning is mostly based on localized molecular orbitals. When reaction paths or energy differences are to be calculated, it is crucial to keep the orbital space consistent for all structures. Inconsistencies in orbital space can lead to unpredictable errors on the potential energy surface. While successful strategies to ensure this consistency have been established for organic and even metal-organic systems, these methods often fail for metal clusters or nanoparticles with a high density of near-degenerate and delocalized molecular orbitals. However, such systems are highly relevant for catalysis. Accurate yet feasible quantum-mechanical ab initio calculations are therefore highly desired. In this work, we present an approach based on the subsystem projected atomic orbital decomposition algorithm that allows us to ensure automated and consistent partitioning even for systems with delocalized and near-degenerate molecular orbitals and demonstrate the validity of this method for the binding energies of small molecules on transition-metal clusters.

del Campo J. Charge delocalization error in Piris Natural Orbital Functionals

Piris Natural Orbital Functionals (PNOF) have been recognized as a low-scaling alternative to study strong correlated systems. In this work, we address the performance of the fifth functional (PNOF5) and the seventh functional (PNOF7) to deal with another common problem, the charge delocalization error. The effects of this problem can be observed in charged systems of repeated well-separated fragments, where the energy should be the sum of the charged and neutral fragments, regardless of how the charge is distributed. In practice, an energetic overstabilization of fractional charged fragments leads to a preference for having the charge delocalized throughout the system. To establish the performance of PNOF functionals regarding charge delocalization error, charged chains of helium atoms and the W4-17-MR set molecules were used as base fragments and their energy, charge distribution and correlation regime were studied. It was found that PNOF5 prefers localized charge distributions, while PNOF7 improves the treatment of interpair static correlation and tends to the correct energetic limit for several cases, although a preference for delocalized charge distributions may arise in highly strong correlation regimes. Overall, it is concluded that PNOF functionals can simultaneously deal with static correlation and charge delocalization errors, resulting in a promising choice to study charge-related problems.

Heterodinuclear Mg(II)M(II) (M=Cr, Mn, Fe, Co, Ni, Cu and Zn) Complexes for the Ring Opening Copolymerization of Carbon Dioxide/Epoxide and Anhydride/Epoxide

The catalysed ring opening copolymerizations (ROCOP) of carbon dioxide/epoxide or anhydride/epoxide are controlled polymerizations that access useful polycarbonates and polyesters. Here, a systematic investigation of a series of heterodinuclear Mg(II)M(II) complexes reveals which metal combinations are most effective. The complexes combine different first row transition metals (M(II)) from Cr(II) to Zn(II), with Mg(II); all complexes are coordinated by the same macrocyclic ancillary ligand and by two acetate co-ligands. The complex syntheses and characterization data, as well as the polymerization data, for both carbon dioxide/cyclohexene oxide (CHO) and endo-norbornene anhydride (NA)/cyclohexene oxide, are reported. The fastest catalyst for both polymerizations is Mg(II)Co(II) which shows propagation rate constants (k_p ) of 34.7 mM^-1  s^-1 (CO₂ ) and 75.3 mM^-1  s^-1 (NA) (100 °C). The Mg(II)Fe(II) catalyst also shows excellent performances with equivalent rates for CO₂ /CHO ROCOP (k_p =34.7 mM^-1  s^-1 ) and may be preferable in terms of metallic abundance, low cost and low toxicity. Polymerization kinetics analyses reveal that the two lead catalysts show overall second order rate laws, with zeroth order dependencies in CO₂ or anhydride concentrations and first order dependencies in both catalyst and epoxide concentrations. Compared to the homodinuclear Mg(II)Mg(II) complex, nearly all the transition metal heterodinuclear complexes show synergic rate enhancements whilst maintaining high selectivity and polymerization control. These findings are relevant to the future design and optimization of copolymerization catalysts and should stimulate broader investigations of synergic heterodinuclear main group/transition metal catalysts.

On the Use of Normalized Metrics for Density Sensitivity Analysis in DFT

In the past years, there has been a discussion about how the errors in density functional theory might be related to errors in the self-consistent densities obtained from different density functional approximations. This, in turn, brings up the discussion about the different ways in which we can measure such errors and develop metrics that assess the sensitivity of calculated energies to changes in the density. It is important to realize that there cannot be a unique metric in order to look at this density sensitivity, simultaneously needing size-extensive and size-intensive metrics. In this study, we report two metrics that are widely applicable to any density functional approximation. We also show how they can be used to classify different chemical systems of interest with respect to their sensitivity to small variations in the density.

Mechanistic Insights into Co and Fe Quaterpyridine-Based CO2 Reduction Catalysts: Metal-Ligand Orbital Interaction as the Key Driving Force for Distinct Pathways

Both [Co^II(qpy)(H₂O)₂]²⁺ and [Fe^II(qpy)(H₂O)₂]²⁺ (with qpy = 2,2':6',2″:6'',2‴-quaterpyridine) are efficient homogeneous electrocatalysts and photoelectrocatalysts for the reduction of CO₂ to CO. The Co catalyst is more efficient in the electrochemical reduction, while the Fe catalyst is an excellent photoelectrocatalyst ( ACS Catal. 2018, 8, 3411-3417). This work uses density functional theory to shed light on the contrasting catalytic pathways. While both catalysts experience primarily ligand-based reductions, the second reduction in the Co catalyst is delocalized onto the metal via a metal-ligand bonding interaction, causing a spin transition and a distorted ligand framework. This orbital interaction explains the experimentally observed mild reduction potential and slow kinetics of the second reduction. The decreased hardness and doubly occupied d_z²-orbital facilitate a σ-bond with the CO₂-π* in an η¹-κC binding mode. CO₂ binding is only possible after two reductions resulting in an EEC mechanism (E = electron transfer, C = chemical reaction), and the second protonation is rate-limiting. In contrast, the Fe catalyst maintains a Lewis acidic metal center throughout the reduction process because the metal orbitals do not strongly mix with the qpy-π* orbitals. This allows binding of the activated CO₂ in an η²-binding mode. This interaction stabilizes the activated CO₂ via a π-type interaction of a Fe-t_2g orbital and the CO₂-π* and a dative bond of the oxygen lone pair. This facilitates CO₂ binding to a singly reduced catalyst resulting in an ECE mechanism. The barrier for CO₂ addition and the second protonation are higher than those for the Co catalyst and rate-limiting.

Density Functional Theory for Molecule-Metal Surface Reactions: When Does the Generalized Gradient Approximation Get It Right, and What to Do If It Does Not

While density functional theory (DFT) is perhaps the most used electronic structure theory in chemistry, many of its practical aspects remain poorly understood. For instance, DFT at the generalized gradient approximation (GGA) tends to fail miserably at describing gas-phase reaction barriers, while it performs surprisingly well for many molecule-metal surface reactions. GGA-DFT also fails for many systems in the latter category, and up to now it has not been clear when one may expect it to work. We show that GGA-DFT tends to work if the difference between the work function of the metal and the molecule's electron affinity is greater than ∼7 eV and to fail if this difference is smaller, with sticking of O₂ on Al(111) being a spectacular example. Using dynamics calculations we show that, for this system, the DFT problem may be solved as done for gas-phase reactions, i.e., by resorting to hybrid functionals, but using screening at long-range to obtain a correct description of the metal. Our results suggest the GGA error in the O₂ + Al(111) barrier height to be functional driven. Our results also suggest the possibility to compute potential energy surfaces for the difficult-to-treat systems with computationally cheap nonself-consistent calculations in which a hybrid functional is applied to a GGA density.

Dioxygen Binding to all 3d, 4d, and 5d Transition Metals from Coupled-Cluster Theory

Understanding how transition metals bind and activate dioxygen (O₂ ) is limited by experimental and theoretical uncertainties, making accurate quantum mechanical descriptors of interest. Here we report coupled-cluster CCSD(T) energies with large basis sets and vibrational and relativistic corrections for 160 3d, 4d, and 5d metal-O₂ systems. We define four reaction energies (120 in total for the 30 metals) that quantify O-O activation and reveal linear relationships between metal-oxygen and O-O binding energies. The CCSD(T) data can be combined with thermochemical cycles to estimate chemisorption and physisorption energies for each metal from metal oxide embedding energies, in good correlation with atomization enthalpies (R² =0.75). Spin-geometry variations can break the linearities, of interest to circumventing the Sabatier principle. Pt, Pd, Co, and Fe form a distinct group with the weakest O₂ binding. R² up to 0.84 between surface adsorption energies and our energies for MO₂ systems indicate relevance also to real catalytic systems.

Using electronegativity and hardness to test density functionals

Density functional theory (DFT) is used in thousands of papers each year, yet lack of universality reduces DFT's predictive capacity, and functionals may produce energy-density imbalances. The absolute electronegativity (χ) and hardness (η) directly reflect the energy-density relationship via the chemical potential ∂E/∂N and we thus hypothesized that they probe universality. We studied χ and η for atoms Z = 1-36 using 50 diverse functionals covering all major classes. Very few functionals describe both χ and η well. η benefits from error cancellation, whereas χ is marred by error propagation from IP and EA; thus, almost all standard GGA and hybrid functionals display a plateau in the MAE at ∼0.2 eV-0.3 eV for η. In contrast, variable performance for χ indicates problems in describing the chemical potential by DFT. The accuracy and precision of a functional is far from linearly related, yet for a universal functional, we expect linearity. Popular functionals such as B3LYP, PBE, and revPBE perform poorly for both properties. Density sensitivity calculations indicate large density-derived errors as occupation of degenerate p- and d-orbitals causes "non-universality" and large dependency on exact exchange. Thus, we argue that performance for χ for the same systems is a hallmark of an important aspect of universality by probing ∂E/∂N. With this metric, B98, B97-1, PW6B95D3, MN-15, rev-TPSS, HSE06, and APFD are the most "universal" among the tested functionals. B98 and B97-1 are accurate for very diverse metal-ligand bonds, supporting that a balanced description of ∂E/∂N and ∂E²/∂N², via χ and η, is probably a first simple probe of universality.

Large-scale comparison of 3d and 4d transition metal complexes illuminates the reduced effect of exchange on second-row spin-state energetics

The origin of distinct 3d vs. 4d transition metal complex sensitivity to exchange is explored over a large data set.