Trends in Strong Chemical Bonding in C2, CN, CN-, CO, N2, NO, NO+, and O2

Overlap properties of chemical bonds in generic systems including unusual bonding situations

Chemical bond is a ubiquitous and fundamental concept in chemistry, in which the overlap plays a defining role. By using a new approach based on localized molecular orbitals, the overlap properties, e.g., polarizability [Formula: see text], population p_OP, intra [Formula: see text], and inter [Formula: see text] repulsions, and density ρ_OP, of polyatomic systems were calculated, analyzed, and correlated. Several trends are shown for these properties, which are rationalized by the balance of some well-known effects, such as, electron donor/withdrawing character and electronegativity. The overlap properties of unusual bonds are also analyzed, revealing an OZn₄(OOCH)₆ structure with four equivalent Zn-O chemical bonds with overlap properties like the O-O bond in H₂O₂, while in protonated methane [Formula: see text], it is observed that a CH₃⋯[Formula: see text] bond pattern at the equilibrium structure changes to a [Formula: see text]⋯H₂ pattern upon dissociation. Charge-shift resonance energies, atom-in-molecule properties, and the lone-pair-bond-weakening effects are related to the overlap properties, which can provide alternative views and insights into chemical bonds. Graphical abstract A chemical bond analysis approach based on its overlap properties is presented for the first time. The model was applied directly to 25 diatomics and for 28 bonds in polytomics employing localized molecular orbitals. Correlations of the overlap properties with the charge-shift resonance energies and with atom-in-molecule (AIM) properties were uncovered. In addition, it provided insights into the Zn-O bonds in the unusual OZn₄(OOCH)₆ system as well as in the bonding patterns of [Formula: see text] at equilibrium and upon dissociation.

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.

The dicarbon bonding puzzle viewed with photoelectron imaging

Bonding in the ground state of C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}_{2}$$\end{document}2 is still a matter of controversy, as reasonable arguments may be made for a dicarbon bond order of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$2$$\end{document}2, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$3$$\end{document}3, or \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$4$$\end{document}4. Here we report on photoelectron spectra of the C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}_{2}^{-}$$\end{document}2− anion, measured at a range of wavelengths using a high-resolution photoelectron imaging spectrometer, which reveal both the ground \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${X}^{1}{\Sigma}_{\mathrm{g}}^{+}$$\end{document}X1Σg+ and first-excited \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${a}^{3}{\Pi}_{{\mathrm{u}}}$$\end{document}a3Πu electronic states. These measurements yield electron angular anisotropies that identify the character of two orbitals: the diffuse detachment orbital of the anion and the highest occupied molecular orbital of the neutral. This work indicates that electron detachment occurs from predominantly \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$s$$\end{document}s-like (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$3{\sigma}_{\mathrm{g}}$$\end{document}3σg) and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p$$\end{document}p-like (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1{\pi }_{{\mathrm{u}}}$$\end{document}1πu) orbitals, respectively, which is inconsistent with the predictions required for the high bond-order models of strongly \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$sp$$\end{document}sp-mixed orbitals. This result suggests that the dominant contribution to the dicarbon bonding involves a double-bonded configuration, with 2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pi$$\end{document}π bonds and no accompanying \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma$$\end{document}σ bond.

In spite of its apparent simplicity, the dicarbon molecule has a bonding structure which is matter of debate. Here the authors measure high-resolution spectra of the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\mathrm{C}}}_{2}$$\end{document}C2 anion by photoelectron imaging, revealing a bonding configuration dominated by a double \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pi$$\end{document}π bond, with no accompanying \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma$$\end{document}σ bond.

λ-Density Functional Valence Bond: A Valence Bond-Based Multiconfigurational Density Functional Theory With a Single Variable Hybrid Parameter

A new valence bond (VB)-based multireference density functional theory (MRDFT) method, named λ-DFVB, is presented in this paper. The method follows the idea of the hybrid multireference density functional method theory proposed by Sharkas et al. (2012). λ-DFVB combines the valence bond self-consistent field (VBSCF) method with Kohn–Sham density functional theory (KS-DFT) by decomposing the electron–electron interactions with a hybrid parameter λ. Different from the Toulouse's scheme, the hybrid parameter λ in λ-DFVB is variable, defined as a function of a multireference character of a molecular system. Furthermore, the E_C correlation energy of a leading determinant is introduced to ensure size consistency at the dissociation limit. Satisfactory results of test calculations, including potential energy surfaces, bond dissociation energies, reaction barriers, and singlet–triplet energy gaps, show the potential capability of λ-DFVB for molecular systems with strong correlation.

The Metal Hydride Problem of Computational Chemistry: Origins and Consequences

Formation and breaking of metal-hydrogen bonds are central to many important catalytic processes such as transition-metal catalyzed ammonia synthesis, hydrogenation reactions, and water splitting, and thus, they require an adequate theoretical description. We studied a data set of all 30 M-H and 30 M⁺-H bonds of the 3d, 4d, and 5d transition series; 50 of these systems have experimentally known bond dissociation enthalpies (BDE). To probe both the limit of low and high coordination number, we also studied a data set of 19 ML _nH complexes. The BDEs were computed using Hartree-Fock (HF), MP2, CCSD, CCSD(T), and 10 diverse density functionals including local, GGA, hybrid GGA, meta hybrid, range-separated, and double hybrids. Our ten most important findings are as follows: (1) HF fails completely to describe the metal hydrogen bond due to its lack of static correlation; (2) this makes post-HF methods such as MP2 and even CCSD(T) perform worse than many density functionals; (3) DFT requires much more HF exchange (∼35% on average) to describe the pure M-H bonds than to describe other metal ligand bonds (0-20%); (4) we design a test to determine if self-interaction error (SIE) is important by correlating DFT errors against a one-electron SIE metric; (5) we show that SIE correlates directly with the DFT errors and thus causes most of the problem; (6) HF-DFT cannot handle these systems because the HF method is too pathological already at the density level; (7) instead, we define and apply a simple metric of electronic abnormality as the difference in PBE energy computed at the self-consistent PBE0 and SVWN densities, and this metric gives appropriate spread and effectively captures density-derived errors; (8) the low electronegativity of the metal enforces a diffuse hydride-like electron density, which make the metal hydrides primary examples of many-electron systems exhibiting SIE already at equilibrium geometries; (9) in the coordinatively saturated ML _nH systems, much less HF exchange is required; i.e., the HF exchange requirements vary drastically with coordination number. Accordingly, DFT is unbalanced for any catalytic process involving both M-H and M-L bonds and changing coordination numbers; (10) importantly, the range-separated and double-hybrid functionals CAM-B3LYP and B2PLYP alone perform well for both M-H and M-L systems and in both limits of low and high coordination number, and at least as well as CCSD(T). This lends hope to a balanced treatment of computational chemistry for all types of M-L bonds at variable coordination number, as required for real catalytic reactions.

Chemical Bond Energies of 3d Transition Metals Studied by Density Functional Theory

Despite their vast importance to inorganic chemistry, materials science, and catalysis, the accuracy of modeling the formation or cleavage of metal-ligand (M-L) bonds depends greatly on the chosen functional and the type of bond in a way that is not systematically understood. In order to approach a state of high-accuracy DFT for rational prediction of chemistry and catalysis, such system-dependencies need to be resolved. We studied 30 different density functionals applied to a "balanced data set" of 60 experimental diatomic M-L bond energies; this data set has no bias toward any d^q configuration, metal, bond type, or ligand as all of these occur to the same extent, and we can therefore identify accuracy bottlenecks. We show that the performance of a functional is very dependent on data set choice, and we dissect these effects into system type. In addition to the use of balanced data sets, we also argue that the precision (rather than just accuracy) of a functional is of interest, measured by standard deviations of the errors. There are distinct system dependencies both in the ligand and metal series: Hydrides are best described by a very large HF exchange percentage, possibly due to self-interaction error, whereas halides are best described by very small (0-10%) HF exchange fractions, and double-bond enforcing oxides and sulfides favor 10-25% HF exchange, as is also average for the full data set. Thus, average HF requirements hide major system-dependent requirements. For late transition metals Co-Zn, HF percentage of 0-10% is favored, whereas for the early transition metals Sc-Fe hybrid functionals with 20% HF exchange or higher are commonly favored. Accordingly, B3LYP is an excellent choice for early d-block but a poor choice for late transition metals. We conclude that DFT intrinsically underestimates the bond strengths of late vs early transition metals, correlating with increased effective nuclear charge. Thus, the revised RPBE, which reduces the overbinding tendency of PBE, is mainly an advantage for the early and mid transition metals and not very much for the late transition metals, i.e. there is a metal-dependent effect of the relative performance of RPBE vs PBE, which are widely used to study adsorption energetics on metal surfaces. Overall, the best performing functionals are PW6B95, the MN15 and MN15-L functionals, and the double hybrid B2PLYP.