Variational Forward–Backward Charge Transfer Analysis Based on Absolutely Localized Molecular Orbitals: Energetics and Molecular Properties

Participation of transition metal atoms in noncovalent bonds

The existence of halogen, chalcogen, pnicogen, and tetrel bonds as variants of noncovalent σ and π-hole bonds is now widely accepted, and many of their properties have been elucidated. The ability of the d-block transition metals to potentially act as Lewis acids in a similar capacity is examined systematically by DFT calculations. Metals examined span the entire range of the d-block from Group 3 to 12, and are selected from several rows of the periodic table. These atoms are placed in a variety of neutral MXn molecules, with X = Cl and O, and paired with a NH3 nucleophile. The resulting M⋯N bonds tend to be stronger than their p-block analogues, many of them with a substantial degree of covalency. The way in which the properties of these bonds is affected by the row and column of the periodic table from which the M atom is drawn, and the number and nature of ligands, is elucidated.

Molecular enhancement of Cu-based catalysts for CO2 electroreduction

The electrochemical carbon dioxide reduction reaction (eCO2RR) represents an effective means of achieving renewable energy storage and a supply of carbon-based raw materials. However, there are still great challenges in selectively producing specific hydrocarbon compounds. The unique ability of the copper (Cu) catalyst to promote proton-coupled electron transfer processes offers clear advantages in generating value-added products. This review presents molecular enhancement strategies for Cu-based catalysts for CO2 electroreduction. We also elucidate the principles of each strategy for enhancing eCO2RR performance, discuss the structure-activity relationships, and propose some promising molecular enhancement strategies. This review will provide guidance for the development of organic-inorganic hybrid Cu-based catalysts as high-performance CO2 electroreduction catalysts.

Highly Activated Terminal Carbon Monoxide Ligand in an Iron-Sulfur Cluster Model of FeMco with Intermediate Local Spin State at Fe

Nitrogenases, the enzymes that convert N2 to NH3, also catalyze the reductive coupling of CO to yield hydrocarbons. CO-coordinated species of nitrogenase clusters have been isolated and used to infer mechanistic information. However, synthetic FeS clusters displaying CO ligands remain rare, which limits benchmarking. Starting from a synthetic cluster that models a cubane portion of the FeMo cofactor (FeMoco), including a bridging carbyne ligand, we report a heterometallic tungsten-iron-sulfur cluster with a single terminal CO coordination in two oxidation states with a high level of CO activation (νCO = 1851 and 1751 cm-1). The local Fe coordination environment (2S, 1C, 1CO) is identical to that in the protein making this system a suitable benchmark. Computational studies find an unusual intermediate spin electronic configuration at the Fe sites promoted by the presence the carbyne ligand. This electronic feature is partly responsible for the high degree of CO activation in the reduced cluster.

Force Decomposition Analysis: A Method to Decompose Intermolecular Forces into Physically Relevant Component Contributions

Computational quantum chemistry can be more than just numerical experiments when methods are specifically adapted to investigate chemical concepts. One important example is the development of energy decomposition analysis (EDA) to reveal the physical driving forces behind intermolecular interactions. In EDA, typically the interaction energy from a good-quality density functional theory (DFT) calculation is decomposed into multiple additive components that unveil permanent and induced electrostatics, Pauli repulsion, dispersion, and charge-transfer contributions to noncovalent interactions. Herein, we formulate, implement, and investigate decomposing the forces associated with intermolecular interactions into the same components. The resulting force decomposition analysis (FDA) is potentially useful as a complement to the EDA to understand chemistry, while also providing far more information than an EDA for data analysis purposes such as training physics-based force fields. We apply the FDA based on absolutely localized molecular orbitals (ALMOs) to analyze interactions of water with sodium and chloride ions as well as in the water dimer. We also analyze the forces responsible for geometric changes in carbon dioxide upon adsorption onto (and activation by) gold and silver anions. We also investigate how the force components of an EDA-based force field for water clusters, namely MB-UCB, compare to those from force decomposition analysis.

Generalization of ETS-NOCV and ALMO-COVP Energy Decomposition Analysis to Connect Any Two Electronic States and Comparative Assessment

Energy decomposition analysis (EDA) is a useful tool for obtaining chemically meaningful insights into molecular interactions. The extended transition-state method with natural orbitals for chemical valence (ETS-NOCV) and the absolutely localized molecular orbital-based method with complementary occupied-virtual pairs (ALMO-COVP) are two successful EDA schemes. Working within ground-state generalized Kohn-Sham density functional theory (DFT), we extend these methods to perform EDA between any two electronic states that can be connected by a unitary transformation of density matrices. A direct proof that the NOCV eigenvalues are symmetric pairs is given, and we also prove that the charge and energy difference defined by ALMO are invariant under certain orbital rotations, allowing us to define COVPs. We point out that ETS is actually a 1-point quadrature to obtain the effective Fock matrix, and though it is reasonably accurate, it can be systematically further improved by adding more quadrature points. We explain why the calculated amount of transferred charge measured by ALMO-COVP is typically much smaller than that of ETS-NOCV and explain why the ALMO-COVP values should be preferred. While the two schemes are independent, ETS-NOCV and ALMO-COVP in fact give a very similar chemical picture for a variety of chemical interactions, including H-H⁺, the transition structure for the Diels-Alder reaction between ethene and butadiene, and two hydrogen-bonded complexes, H₂O···F^- and H₂O···HF.

Modulation of Metal Carbonyl Stretching Frequencies in the Second Coordination Sphere through the Internal Stark Effect

Spectroscopic and computational examination of a homologous series of rhodium(I) pybox carbonyl complexes has revealed a correlation between the conformation of the flanking aryl-substituted oxazoline donors and the carbonyl stretching frequency. This relationship is also observed experimentally for octahedral rhodium(III) and ruthenium(II) variants and cannot be explained through the classical, Dewar-Chatt-Duncanson, interpretation of metal-carbonyl bonding. Instead, these findings are reconciled by local changes in the magnitude of the electric field that is projected along the metal-carbonyl vector: the internal Stark effect.

Revisiting the Bonding Model for Gold(I) Species: The Importance of Pauli Repulsion Revealed in a Gold(I)-Cyclobutadiene Complex

Understanding the bonding of gold(I) species has been central to the development of gold(I) catalysis. Herein, we present the synthesis and characterization of the first gold(I)-cyclobutadiene complex, accompanied with bonding analysis by state-of-the-art energy decomposition analysis methods. Analysis of possible coordination modes for the new species not only confirms established characteristics of gold(I) bonding, but also suggests that Pauli repulsion is a key yet hitherto overlooked element. Additionally, we obtain a new perspective on gold(I)-bonding by comparison of the gold(I)-cyclobutadiene to congeners stabilized by p-, d-, and f-block metals. Consequently, we refine the gold(I) bonding model, with a delicate interplay of Pauli repulsion and charge transfer as the key driving force for various coordination motifs. Pauli repulsion is similarly determined as a significant interaction in AuI -alkyne species, corroborating this revised understanding of AuI bonding.

Dissociation of HCl in water nanoclusters: an energy decomposition analysis perspective

As known, small HCl-water nanoclusters display a particular dissociation behaviour, whereby at least four water molecules are required for the ionic dissociation of HCl. In this work, we examine how intermolecular interactions promote the ionic dissociation of such nanoclusters. To this end, a set of 45 HCl-water nanoclusters with up to four water molecules is introduced. Energy decomposition analysis based on absolutely localized molecular orbitals (ALMO-EDA) is employed in order to study the importance of frozen interaction, dispersion, polarization, and charge-transfer for the dissociation. The vertical ALMO-EDA scheme is applied to HCl-water clusters along a proton-transfer coordinate varying the amount of spectator water molecules. The corresponding ALMO-EDA results show a clear preference for the dissociated cluster only in the presence of four water molecules. Our analysis of adiabatic ALMO-EDA results reveals a push-pull mechanism for the destabilization of the HCl bond based on the synergy between forward and backward charge-transfer.

Variational Energy Decomposition Analysis of Charge-Transfer Interactions between Metals and Ligands in Carbonyl Complexes

Variational energy decomposition analyses have been presented to quantify the σ-dative, ligand-to-metal forward charge transfer (CT) and the π-conjugative, metal-to-ligand backward charge delocalization on a series of isolelectronic transition-metal carbonyl complexes M(CO)₆, including M = Ti^2-, V^-, Cr, Mn⁺, and Fe²⁺. Although the qualitative features of these energy terms are understood, well-defined quantitative studies have been scarce. Consistent with early findings, electrostatic and Pauli exchange effects play a key role in σ-donation, resulting in blue shifts in ligand vibrational frequency in the complex geometries. Excluding chemical bonding interactions between the CO ligand and the metal fragments in the energy decomposition analysis, we found that loosely bound electrostatic complexes can be formed at a longer metal-to-ligand distance due to the exponential decay of Pauli exchange. In all complexes, the overall binding stabilization can be attributed to CT effects, with opposing trends between σ-donation and π-back bonding that follows an order of Ti^2- (4.4) > V^1- (2.6) > Cr (1.5) > Mn¹⁺ (1.1) > Fe²⁺ (0.5) in π-to-σ CT ratio. These electronic and energetic features are mirrored in the vibrational frequency shifts induced by different factors. The present investigation may help stimulate the use of energy decomposition techniques to understand the structure and activity of metallocatalysts using density functional theory.

Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.

Exact-two-component block-localized wave function: A simple scheme for the automatic computation of relativistic ΔSCF

Block-localized wave function is a useful method for optimizing constrained determinants. In this article, we extend the generalized block-localized wave function technique to a relativistic two-component framework. Optimization of excited state determinants for two-component wave functions presents a unique challenge because the excited state manifold is often quite dense with degenerate states. Furthermore, we test the degree to which certain symmetries result naturally from the ΔSCF optimization such as time-reversal symmetry and symmetry with respect to the total angular momentum operator on a series of atomic systems. Variational optimizations may often break the symmetry in order to lower the overall energy, just as unrestricted Hartree-Fock breaks spin symmetry. Overall, we demonstrate that time-reversal symmetry is roughly maintained when using Hartree-Fock, but less so when using Kohn-Sham density functional theory. Additionally, maintaining total angular momentum symmetry appears to be system dependent and not guaranteed. Finally, we were able to trace the breaking of total angular momentum symmetry to the relaxation of core electrons.

Revealing the nature of electron correlation in transition metal complexes with symmetry breaking and chemical intuition

In this work, we provide a nuanced view of electron correlation in the context of transition metal complexes, reconciling computational characterization via spin and spatial symmetry breaking in single-reference methods with qualitative concepts from ligand-field and molecular orbital theories. These insights provide the tools to reliably diagnose the multi-reference character, and our analysis reveals that while strong (i.e., static) correlation can be found in linear molecules (e.g., diatomics) and weakly bound and antiferromagnetically coupled (monometal-noninnocent ligand or multi-metal) complexes, it is rarely found in the ground-states of mono-transition-metal complexes. This leads to a picture of static correlation that is no more complex for transition metals than it is, e.g., for organic biradicaloids. In contrast, the ability of organometallic species to form more complex interactions, involving both ligand-to-metal σ-donation and metal-to-ligand π-backdonation, places a larger burden on a theory's treatment of dynamic correlation. We hypothesize that chemical bonds in which inter-electron pair correlation is non-negligible cannot be adequately described by theories using MP2 correlation energies and indeed find large errors vs experiment for carbonyl-dissociation energies from double-hybrid density functionals. A theory's description of dynamic correlation (and to a less important extent, delocalization error), which affects relative spin-state energetics and thus spin symmetry breaking, is found to govern the efficacy of its use to diagnose static correlation.

From Intermolecular Interaction Energies and Observable Shifts to Component Contributions and Back Again: A Tale of Variational Energy Decomposition Analysis

Quantum chemistry in the form of density functional theory (DFT) calculations is a powerful numerical experiment for predicting intermolecular interaction energies. However, no chemical insight is gained in this way beyond predictions of observables. Energy decomposition analysis (EDA) can quantitatively bridge this gap by providing values for the chemical drivers of the interactions, such as permanent electrostatics, Pauli repulsion, dispersion, and charge transfer. These energetic contributions are identified by performing DFT calculations with constraints that disable components of the interaction. This review describes the second-generation version of the absolutely localized molecular orbital EDA (ALMO-EDA-II). The effects of different physical contributions on changes in observables such as structure and vibrational frequencies upon complex formation are characterized via the adiabatic EDA. Example applications include red- versus blue-shifting hydrogen bonds; the bonding and frequency shifts of CO, N₂, and BF bound to a [Ru(II)(NH₃)₅]^{2 +} moiety; and the nature of the strongly bound complexes between pyridine and the benzene and naphthalene radical cations. Additionally, the use of ALMO-EDA-II to benchmark and guide the development of advanced force fields for molecular simulation is illustrated with the recent, very promising, MB-UCB potential.

A non-perturbative pairwise-additive analysis of charge transfer contributions to intermolecular interaction energies

A non-perturbative scheme for complete decomposition of energy and charge associated with charge transfer interaction into pairwise additive components.

Generalization of Block-Localized Wave Function for Constrained Optimization of Excited Determinants

The block-localized wave function method is useful to provide insights on chemical bonding and intermolecular interactions through energy decomposition analysis. The method relies on block localization of molecular orbitals (MOs) by constraining the orbitals to basis functions within given blocks. Here, a generalized block-localized orbital (GBLO) method is described to allow both physically localized and delocalized MOs to be constrained in orbital-block definitions. Consequently, GBLO optimization can be conveniently tailored by imposing specific constraints. The GBLO method is illustrated by three examples: (1) constrained polarization response orbitals through dipole and quadrupole perturbation in a water dimer complex, (2) the ground and first excited-state potential energy curves of ethene about its C-C bond rotation, and (3) excitation energies of double electron excited states. Multistate density functional theory is used to determine the energies of the adiabatic ground and excited states using a minimal active space (MAS) comprising specifically charge-constrained and excited determinant configurations that are variationally optimized by the GBLO method. We find that the GBLO expansion that includes delocalized MOs in configurational blocks significantly reduces computational errors in comparison with physical block localization, and the computed ground- and excited-state energies are in good accordance with experiments and results obtained from multireference configuration interaction calculations.

The Nature of Nonclassical Carbonyl Ligands Explained by Kohn-Sham Molecular Orbital Theory

When carbonyl ligands coordinate to transition metals, their bond distance either increases (classical) or decreases (nonclassical) with respect to the bond length in the isolated CO molecule. C−O expansion can easily be understood by π‐back‐donation, which results in a population of the CO's π*‐antibonding orbital and hence a weakening of its bond. Nonclassical carbonyl ligands are less straightforward to explain, and their nature is still subject of an ongoing debate. In this work, we studied five isoelectronic octahedral complexes, namely Fe(CO)₆²⁺, Mn(CO)₆⁺, Cr(CO)₆, V(CO)₆⁻ and Ti(CO)₆²⁻, at the ZORA‐BLYP/TZ2P level of theory to explain this nonclassical behavior in the framework of Kohn–Sham molecular orbital theory. We show that there are two competing forces that affect the C−O bond length, namely electrostatic interactions (favoring C−O contraction) and π‐back‐donation (favoring C−O expansion). It is a balance between those two terms that determines whether the carbonyl is classical or nonclassical. By further decomposing the electrostatic interaction ΔV_elstat into four fundamental terms, we are able to rationalize why ΔV_elstat gives rise to the nonclassical behavior, leading to new insights into the driving forces behind C−O contraction.

Probing radical-molecule interactions with a second generation energy decomposition analysis of DFT calculations using absolutely localized molecular orbitals

Intermolecular interactions between radicals and closed-shell molecules are ubiquitous in chemical processes, ranging from the benchtop to the atmosphere and extraterrestrial space. While energy decomposition analysis (EDA) schemes for closed-shell molecules can be generalized for studying radical-molecule interactions, they face challenges arising from the unique characteristics of the electronic structure of open-shell species. In this work, we introduce additional steps that are necessary for the proper treatment of radical-molecule interactions to our previously developed unrestricted Absolutely Localized Molecular Orbital (uALMO)-EDA based on density functional theory calculations. A "polarize-then-depolarize" (PtD) scheme is used to remove arbitrariness in the definition of the frozen wavefunction, rendering the ALMO-EDA results independent of the orientation of the unpaired electron obtained from isolated fragment calculations. The contribution of radical rehybridization to polarization energies is evaluated. It is also valuable to monitor the wavefunction stability of each intermediate state, as well as their associated spin density profiles, to ensure the EDA results correspond to a desired electronic state. These radical extensions are incorporated into the "vertical" and "adiabatic" variants of uALMO-EDA for studies of energy changes and property shifts upon complexation. The EDA is validated on two model complexes, H₂O˙F and FH˙OH. It is then applied to several chemically interesting radical-molecule complexes, including the sandwiched and T-shaped benzene dimer radical cation, complexes of pyridine with benzene and naphthalene radical cations, binary and ternary complexes of the hydroxyl radical with water (˙OH(H₂O) and ˙OH(H₂O)₂), and the pre-reactive complexes and transition states in the ˙OH + HCHO and ˙OH + CH₃CHO reactions. These examples suggest that this second generation uALMO-EDA is a useful tool for furthering one's understanding of both energetic and property changes associated with radical-molecule interactions.