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.

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.

Metal-CO Bonding in Mononuclear Transition Metal Carbonyl Complexes

DFT calculations have been carried out for coordinatively saturated neutral and charged carbonyl complexes [M(CO) _n ] ^q where M is a metal atom of groups 2-10. The model compounds M(CO)₂ (M = Ca, Sr, Ba) and the experimentally observed [Ba(CO)]⁺ were also studied. The bonding situation has been analyzed with a variety of charge and energy partitioning approaches. It is shown that the Dewar-Chatt-Duncanson model in terms of M ← CO σ-donation and M → CO π-backdonation is a valid approach to explain the M-CO bonds and the trend of the CO stretching frequencies. The carbonyl ligands of the neutral complexes carry a negative charge, and the polarity of the M-CO bonds increases for the less electronegative metals, which is particularly strong for the group 4 and group 2 atoms. The NBO method delivers an unrealistic charge distribution in the carbonyl complexes, while the AIM approach gives physically reasonable partial charges that are consistent with the EDA-NOCV calculations and with the trend of the C-O stretching frequencies. The AdNDP method provides delocalized MOs which are very useful models for the carbonyl complexes. Deep insight into the nature of the metal-CO bonds and quantitative information about the strength of the [M] ← (CO)₈ σ-donation and [M(d)] → (CO)₈ π-backdonation visualized by the deformation densities are provided by the EDA-NOCV method. The large polarity of the M-CO π orbitals toward the CO end in the alkaline earth octacarbonyls M(CO)₈ (M = Ca, Sr, Ba) leads to small values for the delocalization indices δ(M-C) and δ(M···O) and significant overlap between adjacent CO groups, but the origin of the charge migration and the associated red-shift of the C-O stretching frequencies is the [M(d)] → (CO)₈ π-backdonation. The heavier alkaline earth metals calcium, strontium and barium use their s/d valence orbitals for covalent bonding. They are therefore to be assigned to the transition metals.

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.

Carbon monoxide insertion at a heavy p-block element: unprecedented formation of a cationic bismuth carbamoyl

Major advances in the chemistry of 5th and 6th row heavy p-block element compounds have recently uncovered intriguing reactivity patterns towards small molecules such as H2, CO2, and ethylene. However, well-defined, homogeneous insertion reactions with carbon monoxide, one of the benchmark substrates in this field, have not been reported to date. We demonstrate here, that a cationic bismuth amide undergoes facile insertion of CO into the Bi-N bond under mild conditions. This approach grants direct access to the first cationic bismuth carbamoyl species. Its characterization by NMR, IR, and UV/vis spectroscopy, elemental analysis, single-crystal X-ray analysis, cyclic voltammetry, and DFT calculations revealed intriguing properties, such as a reversible electron transfer at the bismuth center and an absorption feature at 353 nm ascribed to a transition involving σ- and π-type orbitals of the bismuth-carbamoyl functionality. A combined experimental and theoretical approach provided insight into the mechanism of CO insertion. The substrate scope could be extended to isonitriles.

Ligand Effect on Bonding in Gold(III) Carbonyl Complexes

We quantitatively assess the Dewar-Chatt-Duncanson (DCD) components of the Au(III)-CO bond and the charge density polarization at the CO, in a series of neutral, cationic, and dicationic bis- and monocyclometalated gold(III) complexes via charge-displacement (CD) analysis. A striking feature concerns the very small net electron charge flux from CO to the metal fragment which is unexpectedly stable toward both the charge of the complex and the oxidation state of gold (I, III). All systems exhibit a similar trend for the σ charge rearrangement in the region of the carbonyl bond, where, by contrast, the π back-donation trend variation is large, which is strictly correlated to the change in CO bond distance and the shift in CO stretching frequencies, in close analogy with the gold(I) carbonyl complexes. In the whole series of gold(III) compounds, a large Au(III) ← CO σ donation is measured (from 0.19 to 0.31 electrons), as well as a significant Au(III) → CO π back-donation (from -0.09 to -0.22 electrons), which however is not generally able to completely balance the polarization of the CO π electrons in the direction from oxygen to carbon (C ← O) induced by the presence of the metal fragment [LAu(III)]^0/+1/+2. Surprisingly, all the gold(III) complexes in the series are characterized by a very small anisotropy in the Au(III) → CO in-plane and out-of-plane π back-donation components, in sharp contrast with the marked anisotropy found before for the experimentally characterized [(C^N^C)Au(III)CO]⁺ complex. A first attempt to figure out a rationale on the bonding/reactivity relationship for Au(III)-CO is made by performing a comparative study with an isostructural [(N^N^C)Pt(II)CO]⁺ complex in a model water-gas shift (WGS) reaction.

How π back-donation quantitatively controls the CO stretching response in classical and non-classical metal carbonyl complexes

The CO stretching response upon coordination to a metal M to form [(L) n M(CO)] m complexes (L is an auxiliary ligand) is investigated in relation to the σ donation and π back-donation components of the M-CO bond and to the electrostatic effect exerted by the ligand-metal fragment. Our analysis encompasses over 30 carbonyls, in which the relative importance of donation, back-donation and electrostatics are varied either through the ligand in a series of [(L)Au(CO)]0/+ gold(i) complexes, or through the metal in a series of anionic, neutral and cationic homoleptic carbonyls. Charge-displacement analysis is used to obtain well-defined, consistent measures of σ donation and π back-donation charges, as well as to quantify the σ and π components of CO polarization. It is found that all complexes feature a comparable charge flow of σ symmetry (both in the M-CO bonding region and in the CO fragment itself), which is therefore largely uncorrelated to CO response. By contrast, π back-donation is exceptionally variable and is found to correlate tightly with the change in CO bond distance, with the shift in CO stretching frequency, and with the extent and direction (C → O or C ← O) of the CO π polarization. As a result, we conclusively show that π back-donation can be an important bond component also in non-classical carbonyls and we provide the framework in which the spectroscopic data on coordinated CO can be used to extract quantitative information on the π donor properties of metal-ligand moieties.

Dispersion-corrected energy decomposition analysis for intermolecular interactions based on the BLW and dDXDM methods

As the simplest variant of the valence bond (VB) theory, the block-localized wave function (BLW) method defines the intermediate electron-localized state self-consistently at the DFT level and can be used to explore the nature of intermolecular interactions in terms of several physically intuitive energy components. Yet, it is unclear how the dispersion interaction affects such a kind of energy decomposition analysis (EDA) as standard density functional approximations neglect the long-range dispersion attractive interactions. Three electron densities corresponding to the initial electron-localized state, optimal electron-localized state, and final electron-delocalized state are involved in the BLW-ED approach; a density-dependent dispersion correction, such as the recently proposed dDXDM approach, can thus uniquely probe the impact of the long-range dispersion effect on EDA results computed at the DFT level. In this paper, we incorporate the dDXDM dispersion corrections into the BLW-ED approach and investigate a range of representative systems such as hydrogen-bonding systems, acid-base pairs, and van der Waals complexes. Results show that both the polarization and charge-transfer energies are little affected by the inclusion of the long-range dispersion effect, which thus can be regarded as an independent energy component in EDA.

Energy decomposition analysis based on a block-localized wavefunction and multistate density functional theory

An interaction energy decomposition analysis method based on the block-localized wavefunction (BLW-ED) approach is described. The first main feature of the BLW-ED method is that it combines concepts of valence bond and molecular orbital theories such that the intermediate and physically intuitive electron-localized states are variationally optimized by self-consistent field calculations. Furthermore, the block-localization scheme can be used both in wave function theory and in density functional theory, providing a useful tool to gain insights on intermolecular interactions that would otherwise be difficult to obtain using the delocalized Kohn-Sham DFT. These features allow broad applications of the BLW method to energy decomposition (BLW-ED) analysis for intermolecular interactions. In this perspective, we outline theoretical aspects of the BLW-ED method, and illustrate its applications in hydrogen-bonding and π-cation intermolecular interactions as well as metal-carbonyl complexes. Future prospects on the development of a multistate density functional theory (MSDFT) are presented, making use of block-localized electronic states as the basis configurations.

The B-H...H-P dihydrogen bonding in ion pair complexes [(CF(3))(3)BH(-)][HPH(3-n)(Me)(n)(+)] (n = 0-3) and its implication in H(2) elimination and activation reactions

The B-H(delta-)...(delta+)H-P dihydrogen bonding (DHB) in ion pair complexes [(CF(3))(3)BH(-)][HPH(3-n)(Me)(n)(+)] (n = 0-3) and its role in the combination of proton and hydride with the release of H(2) or, reversibly, the heterolytic activation of H(2) by Lewis pairs (CF(3))(3)BPH(3-n)(Me)(n) have been theoretically investigated at the MP2 and DFT levels. It is found that the B-H...H-P bonds behave similarly to those in neutral pairs and ion-molecule complexes in most respects, such as the linearity of the H...H-P moiety, the characteristics of the electron transfer and rearrangement, and the topological properties of the DHB critical point, except that in certain cases, a blue-shifting of the H-bond vibrational frequency is observed. In [(CF(3))(3)BH(-)][HPH(3-n)(Me)(n)(+)], the proton shifting within the complexes leads to the formation of the dihydrogen complex B(CF(3))(3)(eta(2)-H(2)), which is followed by a subsequent H(2) release. The stability of B(CF(3))(3)(eta(2)-H(2)) (D(e)/D(0) = 10.8/6.0 kcal/mol) makes the proton-hydride combination proceed in a fashion similar to the protonation reactions in transition-metal hydrides rather than those in group 13 hydrides EH(4)(-) (E = B, Al, Ga). As for the H(2)-splitting reaction R(3)BPR'(3) + H(2) --> [R(3)BH(-)][HPR'(3)(+)], classical Lewis pair (CLP) (CF(3))(3)BPH(3) exhibits a high barrier and results in an unstable ion pair product [(CF(3))(3)BH(-)][HPH(3)(+)] compared with the "frustrated Lewis pair" (FLP) (C(6)F(5))(3)BP(tBu)(3). A detailed analysis of the mechanistic aspects of H(2) activation by (CF(3))(3)BPH(3) and (C(6)F(5))(3)BP(tBu)(3), supported by another CLP (CF(3))(3)BP(tBu)(3) which has a binding energy comparable to (CF(3))(3)BPH(3) but a reaction exothermicity comparable to (C(6)F(5))(3)BP(tBu)(3), allows us to suggest that the low stability of FLP (C(6)F(5))(3)BP(tBu)(3) is the determining factor for the low reaction barrier. The relative stability and other properties of the ion pair products [R(3)BH(-)][HPR'(3)(+)] have also been analyzed. Results strongly support the view proposed by Rokob et al. [ Rokob , T. A. ; Hamza , A. ; Stirling , A. ; Soos , T. ; Papai , I. Angew. Chem., Int. Ed. 2008 , 47 , 2435 ] that the frustration energy lowers the energy barrier and increases the exothermicity of the reaction.