Multiple Bonding in AeN- (Ae = Ca, Sr, Ba)

Quantum chemical calculations using ab initio methods at the MRCI+Q(8,9)/def2-QZVPPD and CCSD(T)/def2-QZVPPD levels as well as density functional theory are reported for the diatomic molecules AeN- (Ae = Ca, Sr, Ba). The nature of the bonds is analyzed with a variety of methods. The anions CaN- and SrN- have electronic triplet (3Π) ground states with nearly identical bond dissociation energies De ~57 kcal/mol calculated at the MRCI+Q(8,9)/def2-QZVPPD level of theory. In contrast, the heavier homologue BaN- has a singlet (1Σ+) ground state, which is only 1.1 kcal/mol below the triplet (3Σ-) state. The computed bond dissociation energy of (1Σ+) BaN- is 68.4 kcal/mol. The calculations at the CCSD(T)-full/def2-QZVPPD and BP86-D3(BJ)/def2-QZVPPD levels of theory are in reasonable agreement with the MRCI+Q(8,9)/def2-QZVPPD data except for the singlet (1Σ+) state, which has a large multireference character. The calculated atomic partial charges given by the CM5, Voronoi and Hirshfeld methods suggest small to medium-sized charge donation toward Ae atom Ae←N- for most electronic states. In contrast, the NBO method predicts for all species medium to large electronic charge donation toward nitrogen Ae→N-, which is due to the neglect of the (n)p AOs of Ae atoms as genuine valence orbitals.

Coordination of a Phosphine-Tethered Aminoborane to Group 10 Metals

While π-complexes of C=C bonds are ubiquitous in organometallic chemistry, analogous complexes of the isoelectronic but strongly polarized B=N double bond of aminoboranes are extremely scarce. To address this gap, a diphosphine-aminoborane ligand (PhDPBAiPr) is introduced and its coordination with group 10 metals is investigated. The B=N bond does not coordinate to the metal in Pt(0) and Pd(II) complexes. In contrast, side-on coordination of the B=N bond is observed in the Ni(0) complex [PhDPBAiPr]Ni(NCPh), and the X-ray crystal structure reveals B-N bond elongation compared to the free ligand. The choice of co-ligand strongly influences the presence or absence of side-on coordination at Ni(0) as evidenced by NMR spectroscopy. While the B=N π-complex is geometrically similar to C=C analogues, a bonding analysis reveals that the interaction of the B=N motif with the electron-rich Ni(0) center is best described as 3c4e hyperbond, in which Ni and N are competing for the empty orbital on B.

Understanding chemistry with the symmetry-decomposed Voronoi deformation density charge analysis

The symmetry-decomposed Voronoi deformation density (VDD) charge analysis is an insightful and robust computational tool to aid the understanding of chemical bonding throughout all fields of chemistry. This method quantifies the atomic charge flow associated with chemical-bond formation and enables decomposition of this charge flow into contributions of (1) orbital interaction types, that is, Pauli repulsive or bonding orbital interactions; (2) per irreducible representation (irrep) of any point-group symmetry of interacting closed-shell molecular fragments; and now also (3) interacting open-shell (i.e., radical) molecular fragments. The symmetry-decomposed VDD charge analysis augments the symmetry-decomposed energy decomposition analysis (EDA) so that the charge flow associated with Pauli repulsion and orbital interactions can be quantified both per atom and per irrep, for example, for σ, π, and δ electrons. This provides detailed insights into fundamental aspects of chemical bonding that are not accessible from EDA.

Revisiting sp2 Dilithio Methandiides: From Geometric Curiosity to Simple Bonding Description

The reported tetracoordinate dilithio methandiide complex from Liddle and co-workers (1) is investigated from a coordination chemistry perspective, to probe the origin of its intriguing geometry. Through the application of a variety of computational techniques, non-covalent (steric, electrostatic) interactions are found to be dominant. Further, we arrive at a bonding description which emphasizes the tricoordinate sp2 -hybridized nature of the central methandiide carbon, differing somewhat from the original proposal. Thus, 1 is distinct from other dilithio methandiides since it contains only one C-Li σ-bond, and is found to be comparable to a simple aryllithium compound, phenyllithium.

Bismuth Infrared Star: Being at a Glance

Bismuth polycations have garnered significant attention from researchers due to their extraordinary and counter-intuitive structures and stoichiometries. Despite extensive experimental and theoretical investigations, understanding of the bonding in such clusters remains insufficient. An AdNDP bonding analysis was conducted to elucidate the bonding characteristics using both homoatomic and heteroatomic bismuth clusters with various stoichiometries. Analysis of the calculated nucleus-independent chemical shift data confirmed the aromatic nature of these species. Universal bonding patterns were identified that can be applied to a range of homoatomic and heteroatomic bismuth clusters. Additionally, calculations of absorbance and fluorescence spectra were performed to gain insights into the near-infrared emission and establish a potential correlation between absorbance and the identified bonding patterns.

Cage-size effects on the encapsulation of P2 by fullerenes

The classic pnictogen dichotomy stands for the great contrast between triply bonding very stable N₂ molecules and its heavier congeners, which appear as dimers or oligomers. A banner example involves phosphorus as it occurs in nature as P₄ instead of P₂ , given its weak π-bonds or strong σ-bonds. The P₂ synthetic value has brought Lewis bases and metal coordination stabilization strategies. Herein, we discuss the unrealized encapsulation alternative using the well-known fullerenes' capability to form endohedral and stabilize otherwise unstable molecules. We chose the most stable fullerene structures from C_n (n = 50, 60, 70, 80) and experimentally relevant from C_n (n = 90 and 100) to computationally study the thermodynamics and the geometrical consequences of encapsulating P₂ inside the fullerene cages. Given the size differences between P₂ and P₄ , we show that the fullerenes C₇₀ -C₁₀₀ are suitable cages to side exclude P₄ and host only one molecule of P₂ with an intact triple bond. The thermodynamic analysis indicates that the process is favorable, overcoming the dimerization energy. Additionally, we have evaluated the host-guest interaction to explain the origins of their stability using energy decomposition analysis.

Bonding situations in tricoordinated beryllium phenyl complexes

The bonding situation in the tricoordinated beryllium phenyl complexes [BePh₃ ]^- , [(pyridine)BePh₂ ] and [(trimethylsilyl-N-heterocyclic imine)BePh₂ ] is investigated experimentally and computationally. Comparison of the NMR spectroscopic properties of these complexes and of their structural parameters, which were determined by single crystal X-ray diffraction experiments, indicates the presence of π-interactions. Topology analysis of the electron density reveals elliptical electron density distributions at the bond critical points and the double bond character of the beryllium-element bonds is verified by energy decomposition analysis with the combination of natural orbital for chemical valence. The present beryllium-element bonds are highly polarized and the ligands around the central atom have a strong influence on the degree of π-delocalization. These results are compared to related triarylboranes.

Noble Gas-Silicon Cations: Theoretical Insights into the Nature of the Bond

The structure, stability, and bonding situation of some exemplary noble gas-silicon cations were investigated at the MP2/aVTZ level of theory. The explored species include the mono-coordinated NgSiX3+ (Ng = He-Rn; X = H, F, Cl) and NgSiF22+ (Ng = He-Rn), the di-coordinated Ar2SiX3+ (X = H, F, Cl), and the “inserted” FNgSiF2+ (Ng = Kr, Xe, Rn). The bonding analysis was accomplished by the method that we recently proposed to assay the bonding situation of noblegas compounds. The Ng-Si bonds are generally tight and feature a partial contribution of covalency. In the NgSiX3+, the degree of the Ng-Si interaction mirrors the trends of two factors, namely the polarizability of Ng that increases when going from Ng = He to Ng = Rn, and the Lewis acidity of SiX3+ that decreases in the order SiF3+ > SiH3+ > SiCl3+. For the HeSiX3+, it was also possible to catch peculiar effects referable to the small size of He. When going from the NgSiF3+ to the NgSiF22+, the increased charge on Si promotes an appreciable increase inthe Ng-Si interaction, which becomes truly covalent for the heaviest Ng. The strength of the bond also increases when going from the NgSiF3+ to the “inserted” FNgSiF2+, likely due to the cooperative effect of the adjacent F atom. On the other hand, the ligation of a second Ar atom to ArSiX3+ (X = H, F, Cl), as to form Ar2(SiX3+), produces a weakening of the bond. Our obtained data were compared with previous findings already available in the literature.

Automated Bonding Analysis with Crystal Orbital Hamilton Populations

Understanding crystalline structures based on their chemical bonding is growing in importance. In this context, chemical bonding can be studied with the Crystal Orbital Hamilton Population (COHP), allowing for quantifying interatomic bond strength. Here we present a new set of tools to automate the calculation of COHP and analyze the results. We use the program packages VASP and LOBSTER, and the Python packages atomate and pymatgen. The analysis produced by our tools includes plots, a textual description, and key data in a machine-readable format. To illustrate those capabilities, we have selected simple test compounds (NaCl, GaN), the oxynitrides BaTaO2 N, CaTaO2 N, and SrTaO2 N, and the thermoelectric material Yb14 Mn1 Sb11 . We show correlations between bond strengths and stabilities in the oxynitrides and the influence of the Mn-Sb bonds on the magnetism in Yb14 Mn1 Sb11 . Our contribution enables high-throughput bonding analysis and will facilitate the use of bonding information for machine learning studies.

How 5 f Electron Polarisability Drives Covalency and Selectivity in Actinide N-Donor Complexes

We report a series of isostructural tetravalent actinide (Th, U−Pu) complexes with the N‐donor ligand N,N’‐ethylene‐bis((pyrrole‐2‐yl)methanimine) (H₂L, H₂pyren). Structural data from SC‐XRD analysis reveal [An(pyren)₂] complexes with different An−N_imine versus An−N_pyrrolide bond lengths. Quantum chemical calculations elucidated the bonding situation, including differences in the covalent character of the coordinative bonds. A comparison to the intensely studied analogous N,N′‐ethylene‐bis(salicylideneimine) (H₂salen)‐based complexes [An(salen)₂] displays, on average, almost equal electron sharing of pyren or salen with the An^IV, pointing to a potential ligand‐cage‐driven complex stabilisation. This is shown in the fixed ligand arrangement of pyren and salen in the respective An^IV complexes. The overall bond strength of the pure N‐donor ligand pyren to An^IV (An=Th, U, Np, Pu) is slightly weaker than to salen, with the exception of the Pa^IV complex, which exhibits extraordinarily high electron sharing of pyren with Pa^IV. Such an altered ligand preference within the early An^IV series points to a specificity of the 5f¹ configuration, which can be explained by polarisation effects of the 5 f electrons, allowing the strongest f electron backbonding from Pa^IV (5f¹) to the N donors of pyren.

Concerning the Role of σ-Hole in Non-Covalent Interactions: Insights from the Study of the Complexes of ArBeO with Simple Ligands

The structure, stability, and bonding character of some exemplary LAr and L-ArBeO (L = He, Ne, Ar, N₂, CO, F₂, Cl₂, ClF, HF, HCl, NH₃) were investigated by MP2 and coupled-cluster calculations, and by symmetry-adapted perturbation theory. The nature of the stabilizing interactions was also assayed by the method recently proposed by the authors to classify the chemical bonds in noble-gas compounds. The comparative analysis of the LAr and L-ArBeO unraveled geometric and bonding effects peculiarly related to the σ-hole at the Ar atom of ArBeO, including the major stabilizing/destabilizing role of the electrostatic interactionensuing from the negative/positive molecular electrostatic potential of L at the contact zone with ArBeO. The role of the inductive and dispersive components was also assayed, making it possible to discern the factors governing the transition from the (mainly) dispersive domain of the LAr, to the σ-hole domain of the L-ArBeO. Our conclusions could be valid for various types of non-covalent interactions, especially those involving σ-holes of respectable strength such as those occurring in ArBeO.

Ni Oxidation State and Ligand Saturation Impact on the Capability of Octaazamacrocyclic Complexes to Bind and Reduce CO2

Two 15-membered octaazamacrocyclic nickel(II) complexes are investigated by theoretical methods to shed light on their affinity forwards binding and reducing CO2. In the first complex 1[NiIIL]0, the octaazamacrocyclic ligand is grossly unsaturated (π-conjugated), while in the second 1[NiIILH]2+ one, the macrocycle is saturated with hydrogens. One and two-electron reductions are described using Mulliken population analysis, quantum theory of atoms in molecules, localized orbitals, and domain averaged fermi holes, including the characterization of the Ni-CCO2 bond and the oxidation state of the central Ni atom. It was found that in the [NiLH] complex, the central atom is reduced to Ni0 and/or NiI and is thus able to bind CO2 via a single σ bond. In addition, the two-electron reduced 3[NiL]2- species also shows an affinity forwards CO2.

Bonding in Barium Boryloxides, Siloxides, Phenoxides and Silazides: A Comparison with the Lighter Alkaline Earths

Barium complexes ligated by bulky boryloxides [OBR₂ ]^- (where R=CH(SiMe₃ )₂ , 2,4,6-ⁱ Pr₃ -C₆ H₂ or 2,4,6-(CF₃ )₃ -C₆ H₂ ), siloxide [OSi(SiMe₃ )₃ ]^- , and/or phenoxide [O-2,6-Ph₂ -C₆ H₃ ]^- , have been prepared. A diversity of coordination patterns is observed in the solid state for both homoleptic and heteroleptic complexes, with coordination numbers ranging between 2 and 4. The identity of the bridging ligand in heteroleptic dimers [Ba(μ₂ -X₁ )(X₂ )]₂ depends largely on the given pair of ligands X₁ and X₂ . Experimentally, the propensity to fill the bridging position increases according to [OB{CH(SiMe₃ )₂ }₂ )]^- <[N(SiMe₃ )₂ ]^- <[OSi(SiMe₃ )₃ ]^- <[O(2,6-Ph₂ -C₆ H₃ )]^- <[OB(2,4,6-ⁱ Pr₃ -C₆ H₂ )₂ ]^- . This trend is the overall expression of 3 properties: steric constraints, electronic density and σ- and π-donating capability of the negatively charged atom, and ability to generate Ba ⋅ ⋅ ⋅ F, Ba ⋅ ⋅ ⋅ C(π) or Ba ⋅ ⋅ ⋅ H-C secondary interactions. The comparison of the structural motifs in the complexes [Ae{μ₂ -N(SiMe₃ )₂ }(OB{CH(SiMe₃ )₂ }₂ )]₂ (Ae = Mg, Ca, Sr and Ba) suggest that these observations may be extended to all alkaline earths. DFT calculations highlight the largely prevailing ionic character of ligand-Ae bonding in all compounds. The ionic character of the Ae-ligand bond encourages bridging coordination, whereas the number of bridging ligands is controlled by steric factors. DFT computations also indicate that in [Ba(μ₂ -X₁ )(X₂ )]₂ heteroleptic dimers, ligand predilection for bridging vs. terminal positions is dictated by the ability to establish secondary interactions between the metals and the ligands.

Highly Coordinated Heteronuclear Calcium-Iron Carbonyl Cation Complexes [CaFe(CO)n ]+ (n=5-12) with d-d Bonding

Heteronuclear calcium-iron carbonyl cation complexes in the form of [CaFe(CO)n ]+ (n=5-12) are produced in the gas phase. Infrared photodissociation spectroscopy in conjunction with quantum chemical calculations confirm that the n=10 complex is the coordination saturated ion where a Fe(CO)4 fragment is bonded with a Ca(CO)6 fragment through two side-on bridging carbonyl ligands. Bonding analysis indicates that it is best described by the bonding interactions between a [Ca(CO)6 ]2+ dication and an [Fe(CO)4 ]- anion forming a Fe→Ca d-d dative bond in the [(CO)6 Ca-Fe(CO)4 ]+ structure, which enriches the pool of experimentally observed complexes of calcium that mimic transition metal compounds. The molecule is the first example of a heteronuclear carbonyl complex featuring a d-d bond between calcium and a transition metal.

Not Carbon s-p Hybridization, but Coordination Number Determines C-H and C-C Bond Length

A fundamental and ubiquitous phenomenon in chemistry is the contraction of both C−H and C−C bonds as the carbon atoms involved vary, in s–p hybridization, along sp³ to sp² to sp. Our quantum chemical bonding analyses based on Kohn–Sham molecular orbital theory show that the generally accepted rationale behind this trend is incorrect. Inspection of the molecular orbitals and their corresponding orbital overlaps reveals that the above‐mentioned shortening in C−H and C−C bonds is not determined by an increasing amount of s‐character at the carbon atom in these bonds. Instead, we establish that this structural trend is caused by a diminishing steric (Pauli) repulsion between substituents around the pertinent carbon atom, as the coordination number decreases along sp³ to sp² to sp.

Chemical Bonding in Homoleptic Carbonyl Cations [M{Fe(CO)5 }2 ]+ (M=Cu, Ag, Au)

Syntheses of the copper and gold complexes [Cu{Fe(CO)5 }2 ][SbF6 ] and [Au{Fe(CO)5 }2 ][HOB{3,5-(CF3 )2 C6 H3 }3 ] containing the homoleptic carbonyl cations [M{Fe(CO)5 }2 ]+ (M=Cu, Au) are reported. Structural data of the rare, trimetallic Cu2 Fe, Ag2 Fe and Au2 Fe complexes [Cu{Fe(CO)5 }2 ][SbF6 ], [Ag{Fe(CO)5 }2 ][SbF6 ] and [Au{Fe(CO)5 }2 ][HOB{3,5-(CF3 )2 C6 H3 }3 ] are also given. The silver and gold cations [M{Fe(CO)5 }2 ]+ (M=Ag, Au) possess a nearly linear Fe-M-Fe' moiety but the Fe-Cu-Fe' in [Cu{Fe(CO)5 }2 ][SbF6 ] exhibits a significant bending angle of 147° due to the strong interaction with the [SbF6 ]- anion. The Fe(CO)5 ligands adopt a distorted square-pyramidal geometry in the cations [M{Fe(CO)5 }2 ]+ , with the basal CO groups inclined towards M. The geometry optimization with DFT methods of the cations [M{Fe(CO)5 }2 ]+ (M=Cu, Ag, Au) gives equilibrium structures with linear Fe-M-Fe' fragments and D2 symmetry for the copper and silver cations and D4d symmetry for the gold cation. There is nearly free rotation of the Fe(CO)5 ligands around the Fe-M-Fe' axis. The calculated bond dissociation energies for the loss of both Fe(CO)5 ligands from the cations [M{Fe(CO)5 }2 ]+ show the order M=Au (De =137.2 kcal mol-1 )>Cu (De =109.0 kcal mol-1 )>Ag (De =92.4 kcal mol-1 ). The QTAIM analysis shows bond paths and bond critical points for the M-Fe linkage but not between M and the CO ligands. The EDA-NOCV calculations suggest that the [Fe(CO)5 ]→M+ ←[Fe(CO)5 ] donation is significantly stronger than the [Fe(CO)5 ]←M+ →[Fe(CO)5 ] backdonation. Inspection of the pairwise orbital interactions identifies four contributions for the charge donation of the Fe(CO)5 ligands into the vacant (n)s and (n)p AOs of M+ and five components for the backdonation from the occupied (n-1)d AOs of M+ into vacant ligand orbitals.

Nature of the Hydrogen Bond Enhanced Halogen Bond

The factors responsible for the enhancement of the halogen bond by an adjacent hydrogen bond have been quantitatively explored by means of state-of-the-art computational methods. It is found that the strength of a halogen bond is enhanced by ca. 3 kcal/mol when the halogen donor simultaneously operates as a halogen bond donor and a hydrogen bond acceptor. This enhancement is the result of both stronger electrostatic and orbital interactions between the XB donor and the XB acceptor, which indicates a significant degree of covalency in these halogen bonds. In addition, the halogen bond strength can be easily tuned by modifying the electron density of the aryl group of the XB donor as well as the acidity of the hydrogen atoms responsible for the hydrogen bond.

On the Proton-Bound Noble Gas Dimers (Ng-H-Ng)+ and (Ng-H-Ng')+ (Ng, Ng'= He-Xe): Relationships betweenStructure, Stability, and Bonding Character

The structure, stability, and bonding character of fifteen (Ng-H-Ng)⁺ and (Ng-H-Ng’)⁺ (Ng, Ng’ = He-Xe) compounds were explored by theoretical calculations performed at the coupled cluster level of theory. The nature of the stabilizing interactions was, in particular, assayed using a method recently proposed by the authors to classify the chemical bonds involving the noble-gas atoms. The bond distances and dissociation energies of the investigated ions fall in rather large intervals, and follow regular periodic trends, clearly referable to the difference between the proton affinity (PA) of the various Ng and Ng’. These variations are nicely correlated with the bonding situation of the (Ng-H-Ng)⁺ and (Ng-H-Ng’)⁺. The Ng-H and Ng’-H contacts range, in fact, between strong covalent bonds to weak, non-covalent interactions, and their regular variability clearly illustrates the peculiar capability of the noble gases to undergo interactions covering the entire spectrum of the chemical bond.

Generation and Characterization of the C3 O2 - Anion with an Unexpected Unsymmetrical Structure

The carbon suboxide anion C₃ O₂ ^- is generated in solid neon matrix. It is characterized by infrared absorption spectroscopy as well as quantum chemical calculations to have a planar C_s structure where two CO groups with significantly different bond lengths and angles are attached in a zigzag fashion to the central carbon atom. Bonding analysis indicates that it is best described by the bonding interactions between a neutral CO in a triplet excited state and a doublet excited state of CCO^- .

CO-Induced Dinitrogen Fixation and Cleavage Mediated by Boron

The boron atoms react with carbon monoxide and dinitrogen forming the end-on bonded NNBCO complex in solid neon or in nitrogen matrices. The NNBCO complex rearranges to the (η2 -N2 )BCO isomer with a more activated side-on bonded dinitrogen ligand upon visible light excitation. (η2 -N2 )BCO and its weakly CO-coordinated complexes further isomerize to the NBNCO and B(NCO)2 molecules with N-N bond being completely cleaved under UV light irradiation. The geometries, energies and vibrational spectra of the molecules are calculated with quantum chemical methods and the electronic structures are analyzed with charge- and energy-partitioning methods.

Generation and Identification of the Linear OCBNO and OBNCO Molecules with 24 Valence Electrons

Two structural isomers containing five second-row element atoms with 24 valence electrons were generated and identified by matrix-isolation IR spectroscopy and quantum chemical calculations. The OCBNO complex, which is produced by the reaction of boron atoms with mixtures of carbon monoxide and nitric oxide in solid neon, rearranges to the more stable OBNCO isomer on UV excitation. Bonding analysis indicates that the OCBNO complex is best described by the bonding interactions between a triplet-state boron cation with an electron configuration of (2s)0 (2pσ )0 (2pπ )2 and the CO/NO- ligands in the triplet state forming two degenerate electron-sharing π bonds and two ligand-to-boron dative σ bonds.

Group 6 Hexacarbonyls as Ligands for the Silver Cation: Syntheses, Characterization, and Analysis of the Bonding Compared with the Isoelectronic Group 5 Hexacarbonylates

The syntheses of the two novel complexes [Ag{Mo/W(CO)6 }2 ]+ [F-{Al(ORF )3 }2 ]- (RF =C(CF3 )3 ) are reported along with their structural and spectroscopic characterization. The X-ray structure shows that three carbonyl ligands from each M(CO)6 fragment bend towards the silver atom within binding Ag-C distance range. DFT calculations of the free cations [Ag{M(CO)6 }2 ]+ (M=Cr, Mo, W) in the electronic singlet state give equilibrium structures with C2 symmetry with two bridging carbonyl groups from each hexacarbonyl ligand. Similar structures with C2 symmetry (M=Nb) and D2 symmetry (M=V, Ta) are calculated for the isoelectronic group 5 anions [Ag{M(CO)6 }2 ]- (M=V, Nb, Ta). The electronic structure of the cations is analyzed with the QTAIM and EDA-NOCV methods, which provide detailed information about the nature of the chemical bonds between Ag+ and the {M(CO)6 }2 q (q = -2, M = V, Nb, Ta; q = 0, M = Cr, Mo, W) ligands.

Bonding Trends in Tetravalent Th-Pu Monosalen Complexes

The synthesis of three complex series of the form [AnCl2 (salen)(Pyx)2 ] (H2 salen=N,N'-bis(salicylidene)ethylenediamine; Pyx=pyridine, 4-methylpyridine, 3,5-dimethylpyridine) with tetravalent early actinides (An=Th, U, Np, Pu) is reported with the goal to elucidate the affinity of these heavy elements for small neutral N-donor molecules. Structure determination by single-crystal XRD and characterization of bulk powders with infrared spectroscopy reveals isostructurality within each respective series and the same complex conformation in all reported structures. Although the trend of interatomic distances for An-Cl and An-N (imine nitrogen of salen or pyridyl nitrogen of Pyx) was found to reflect an ionic behavior, the trend of the An-O distances can only be described with additional covalent interactions for all elements heavier than thorium. All experimental results are supported by quantum chemical calculations, which confirm the mostly ionic character in the An-N and An-Cl bonds, as well as the highest degree of covalency of the An-O bonds. Structurally, the calculations indicate just minor electronic or steric effects of the additional Pyx substituents on the complex properties.

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.

Beryllium Atom Mediated Dinitrogen Activation via Coupling with Carbon Monoxide

The reactions of laser-ablated beryllium atoms with dinitrogen and carbon monoxide mixtures form the end-on bonded NNBeCO and side-on bonded (η2 -N2 )BeCO isomers in solid argon, which are predicted by quantum chemical calculations to be almost isoenergetic. The end-on bonded complex has a triplet ground state while the side-on bonded isomer has a singlet electronic ground state. The complexes rearrange to the energetically lowest lying NBeNCO isomer upon visible light excitation, which is characterized to be an isocyanate complex of a nitrene derivative with a triplet electronic ground state. A bonding analysis using a charge- and energy decomposition procedure reveals that the electronic reference state of Be in the NNBeCO isomers has an 2s0 2p2 excited configuration and that the metal-ligand bonds can be described in terms of N2 →Be←CO σ donation and concomitant N2 ←Be→CO π backdonation. The results demonstrate that the activation of N2 with the N-N bond being completely cleaved can be achieved via coupling with carbon monoxide mediated by a main group atom.

NBO/NRT Two-State Theory of Bond-Shift Spectral Excitation

We show that natural bond orbital (NBO) and natural resonance theory (NRT) analysis methods provide both optimized Lewis-structural bonding descriptors for ground-state electronic properties as well as suitable building blocks for idealized "diabatic" two-state models of the associated spectroscopic excitations. Specifically, in the framework of single-determinant Hartree-Fock or density functional methods for a resonance-stabilized molecule or supramolecular complex, we employ NBO/NRT descriptors of the ground-state determinant to develop a qualitative picture of the associated charge-transfer excitation that dominates the valence region of the electronic spectrum. We illustrate the procedure for the elementary bond shifts of SN2-type halide exchange reaction as well as the more complex bond shifts in a series of conjugated cyanine dyes. In each case, we show how NBO-based descriptors of resonance-type 3-center, 4-electron (3c/4e) interactions provide simple estimates of spectroscopic excitation energy, bond orders, and other vibronic details of the excited-state PES that anticipate important features of the full multi-configuration description. The deep 3c/4e connections to measurable spectral properties also provide evidence for NBO-based estimates of ground-state donor-acceptor stabilization energies (sometimes criticized as "too large" compared to alternative analysis methods) that are also found to be of proper magnitude to provide useful estimates of excitation energies and structure-dependent spectral shifts.

Interacting Quantum Atoms-A Review

The aim of this review is threefold. On the one hand, we intend it to serve as a gentle introduction to the Interacting Quantum Atoms (IQA) methodology for those unfamiliar with it. Second, we expect it to act as an up-to-date reference of recent developments related to IQA. Finally, we want it to highlight a non-exhaustive, yet representative set of showcase examples about how to use IQA to shed light in different chemical problems. To accomplish this, we start by providing a brief context to justify the development of IQA as a real space alternative to other existent energy partition schemes of the non-relativistic energy of molecules. We then introduce a self-contained algebraic derivation of the methodological IQA ecosystem as well as an overview of how these formulations vary with the level of theory employed to obtain the molecular wavefunction upon which the IQA procedure relies. Finally, we review the several applications of IQA as examined by different research groups worldwide to investigate a wide variety of chemical problems.

Filling a Gap: The Coordinatively Saturated Group 4 Carbonyl Complexes TM(CO)8 (TM=Zr, Hf) and Ti(CO)7

Homoleptic Group 4 metal carbonyl cation and neutral complexes were prepared in the gas phase and/or in solid neon matrix. Infrared spectroscopy studies reveal that both zirconium and hafnium form eight-coordinate carbonyl neutral and cation complexes. In contrast, titanium forms only the six-coordinate Ti(CO)6 + and seven-coordinate Ti(CO)7 . Titanium octacarbonyl Ti(CO)8 is unstable as a result of steric repulsion between the CO ligands. The 20-electron Zr(CO)8 and Hf(CO)8 complexes represent the first experimentally observed homoleptic octacarbonyl neutral complexes of transition metals. The molecules still fulfill the 18-electron rule, because one doubly occupied valence orbital does not mix with any of the metal valence atomic orbitals. Zr(CO)8 and Hf(CO)8 are stable against the loss of one CO because the CO ligands encounter less steric repulsion than Zr(CO)7 and Hf(CO)7 . The heptacarbonyl complexes have shorter metal-CO bonds than that of the octacarbonyl complexes due to stronger electrostatic and covalent bonding, but the significantly smaller repulsive Pauli term makes the octacarbonyl complexes stable.

Molecular Bismuth Cations: Assessment of Soft Lewis Acidity

Three-coordinate cationic bismuth compounds [Bi(diaryl)(EPMe3 )][SbF6 ] have been isolated and fully characterized (diaryl=[(C6 H4 )2 C2 H2 ]2- , E=S, Se). They represent rare examples of molecular complexes with Bi⋅⋅⋅EPR3 interactions (R=monoanionic substituent). The 31 P NMR chemical shift of EPMe3 has been found to be sensitive to the formation of LA⋅⋅⋅EPMe3 Lewis acid/base interactions (LA=Lewis acid). This corresponds to a modification of the Gutmann-Beckett method and reveals information about the hardness/softness of the Lewis acid under investigation. A series of organobismuth compounds, bismuth halides, and cationic bismuth species have been investigated with this approach and compared to traditional group 13 and cationic group 14 Lewis acids. Especially cationic bismuth species have been shown to be potent soft Lewis acids that may prefer Lewis pair formation with a soft (S/Se-based) rather than a hard (O/N-based) donor. Analytical techniques applied in this work include (heteronuclear) NMR spectroscopy, single-crystal X-ray diffraction analysis, and DFT calculations.

LaTe1.82(1): modulated crystal structure and chemical bonding of a chalcogen-deficient rare earth metal polytelluride

Crystals of the rare earth metal polytelluride LaTe_1.82(1), namely, lanthanum telluride (1/1.8), have been grown by molten alkali halide flux reactions and vapour-assisted crystallization with iodine. The two-dimensionally incommensurately modulated crystal structure has been investigated by X-ray diffraction experiments. In contrast to the tetragonal average structure with unit-cell dimensions of a = 4.4996 (5) and c = 9.179 (1) Å at 296 (1) K, which was solved and refined in the space group P4/nmm (No. 129), the satellite reflections are not compatible with a tetragonal symmetry but enforce a symmetry reduction. Possible space groups have been derived by group-subgroup relationships and by consideration of previous reports on similar rare earth metal polychalcogenide structures. Two structural models in the orthorhombic superspace group, i.e. Pmmn(α,β,1/2)000(-α,β,1/2)000 (No. 59.2.51.39) and Pm2₁n(α,β,1/2)000(-α,β,1/2)000 (No. 31.2.51.35), with modulation wave vectors q₁ = αa* + βb* + 1/2c* and q₂ = -αa* + βb* + 1/2c* [α = 0.272 (1) and β = 0.314 (1)], have been established and evaluated against each other. The modulation describes the distribution of defects in the planar [Te] layer, coupled to a displacive modulation due to the formation of different Te anions. The bonding situation in the planar [Te] layer and the different Te anion species have been investigated by density functional theory (DFT) methods and an electron localizability indicator (ELI-D)-based bonding analysis on three different approximants. The temperature-dependent electrical resistance revealed a semiconducting behaviour with an estimated band gap of 0.17 eV.

Reply to the Comment on "Realization of Lewis Basic Sodium Anion in the NaBH3 - Cluster"

We reply to the comment by S. Pan and G. Frenking who challenged our interpretation of the Na^- :→BH₃ dative bond in the recently synthesized NaBH₃ ^- cluster. Our conclusion remains the same as that in our original paper (https://doi.org/10.1002/anie.201907089 and https://doi.org/10.1002/ange.201907089). This conclusion is additionally supported by the energetic pathways and NBO charges calculated at UCCSD and CASMP2(4,4) levels of theory. We also discussed the suitability of the Laplacian of electron density (QTAIM) and Adaptive Natural Density Partitioning (AdNDP) method for bond type assignment. It seems that AdNDP yields more sensible results. This discussion reveals that the complex realm of bonding is full of semantic inconsistencies, and we invite experimentalists and theoreticians to elaborate this topic and find solutions incorporating different views on the dative bond.

Comment on "Realization of Lewis Basic Sodium Anion in the NaBH3 - Cluster"

We challenge the interpretation of the chemical bond in NaBH₃ ^- proposed by Liu et al. We argue that NaBH₃ ^- has an electron-sharing Na-BH₃ ^- covalent bond rather than a dative bond Na^- →BH₃ .

Complexes of helium with neutral molecules: Progress toward a quantitative scale of bonding character

The complexes of helium with nearly 30 neutral molecules (M) were investigated by various techniques of bonding analysis and symmetry-adapted perturbation theory (SAPT). The main investigated function was the local electron energy density H(r), analyzed, in particular, so to estimate the degree of polarization (DoP) of He in the various He(M). As we showed recently (Borocci et al., J. Comput. Chem., 2019, 40, 2318-2328), the DoP is a quantitative index that is generally informative about the role of polarization (induction plus charge transfer [CT]) and dispersion in noncovalent noble gas complexes. As further evidence in this regard, we presently ascertained quantitative correlations between the DoP(He) of the He(M) and indices based on the electron density ρ(r), including the molecular electrostatic potential at the HeM bond critical point, as well as the percentage contributions of induction and dispersion to the SAPT binding energies. Based also on the explicit evaluation of the CT, accomplished through the study of the charge-displacement function, we derived a quantitative scale that ranks the He(M) according to their dispersive, inductive, and CT bonding character. Our taken approach could be conceivably extended to other types of noncovalent complexes.

Embedding a Planar Hypercoordinate Carbon Atom into a [4n+2] π-System

Through delicate tuning of the electronic structure, we report herein a rational design of seventeen new putative global minimum energy structures containing a planar tetra- or pentacoordinate carbon atom embedded in an aromatic hydrocarbon. These structures are the result of replacing three consecutive hydrogen atoms of an aromatic hydrocarbon by less electronegative groups, forming a multicenter σ-bond with the planar hypercoordinate carbon atom and participating in the π-electron delocalization. This strategy that maximizes both mechanical and electronic effects through aromatic architectures can be extended to several molecular combinations to achieve new and diverse compounds containing planar hypercoordinate carbon centers.

Transition-Metal Chemistry of Alkaline-Earth Elements: The Trisbenzene Complexes M(Bz)3 (M=Sr, Ba)

We report the synthesis and spectroscopic identification of the trisbenzene complexes of strontium and barium M(Bz)3 (M=Sr, Ba) in low-temperature Ne matrix. Both complexes are characterized by a D3 symmetric structure involving three equivalent η6 -bound benzene ligands and a closed-shell singlet electronic ground state. The analysis of the electronic structure shows that the complexes exhibit metal-ligand bonds that are typical for transition metal compounds. The chemical bonds can be explained in terms of weak donation from the π MOs of benzene ligands into the vacant (n-1)d AOs of M and strong backdonation from the occupied (n-1)d AO of M into vacant π* MOs of benzene ligands. The metals in these 20-electron complexes have 18 effective valence electrons, and, thus, fulfill the 18-electron rule if only the metal-ligand bonding electrons are counted. The results suggest that the heavier alkaline earth atoms exhibit the full bonding scenario of transition metals.

Octacarbonyl Ion Complexes of Actinides [An(CO)8 ]+/- (An=Th, U) and the Role of f Orbitals in Metal-Ligand Bonding

The octacarbonyl cation and anion complexes of actinide metals [An(CO)₈ ]^+/- (An=Th, U) are prepared in the gas phase and are studied by mass-selected infrared photodissociation spectroscopy. Both the octacarbonyl cations and anions have been characterized to be saturated coordinated complexes. Quantum chemical calculations by using density functional theory show that the [Th(CO)₈ ]⁺ and [Th(CO)₈ ]^- complexes have a distorted octahedral (D_4h ) equilibrium geometry and a doublet electronic ground state. Both the [U(CO)₈ ]⁺ cation and the [U(CO)₈ ]^- anion exhibit cubic structures (O_h ) with a ⁶ A_1g ground state for the cation and a ⁴ A_1g ground state for the anion. The neutral species [Th(CO)₈ ] (O_h ; ¹ A_1g ) and [U(CO)₈ ] (D_4h ; ⁵ B_1u ) have also been calculated. Analysis of their electronic structures with the help on an energy decomposition method reveals that, along with the dominating 6d valence orbitals, there are significant 5f orbital participation in both the [An]←CO σ donation and [An]→CO π back donation interactions in the cations and anions, for which the electronic reference state of An has both occupied and vacant 5f AOs. The trend of the valence orbital contribution to the metal-CO bonds has the order of 6d≫5f>7s≈7p, with the 5f orbitals of uranium being more important than the 5f orbitals of thorium.

Realization of Lewis Basic Sodium Anion in the NaBH3 - Cluster

We report a Na:^- →B dative bond in the NaBH₃ ^- cluster, which was designed on the principle of minimum-energy rupture, prepared by laser vaporization, and characterized by a synergy of anion photoelectron spectroscopy and electronic structure calculations. The global minimum of NaBH₃ ^- features a Na-B bond. Its preferred heterolytic dissociation conforms with the IUPAC definition of dative bond. The lone electron pair revealed on Na and the negative Laplacian of electron density at the bond critical point further confirm the dative nature of the Na-B bond. This study represents the first example of a Lewis adduct with an alkalide as the Lewis base.

Noncovalent Complexes of the Noble-Gas Atoms: Analyzing the Transition from Physical to Chemical Interactions

The bonding character of the noncovalent complexes of the noble-gas (Ng) atoms ranges from nearly purely dispersive contacts to interactions featuring appreciable contributions of induction and charge transfer. In this study, we discuss a new quantitative index that seems peculiarly informative about these diverse bonding situations. This index was termed as the degree of polarization (DoP) of Ng, as it measures, in essence, the Ng polarization promoted by the binding partner. The definition of the DoP(Ng) relies on the analysis of the local electron energy density H(r), and its physical meaning was best appreciated by studying also the charge-displacement function and the molecular electrostatic potential of the investigated benchmark species, that include nearly 60 Ngs complexes of different bonding character. The DoP(Ng) appears of general applicability, and is also positively correlated with other bonding character indices. © 2019 Wiley Periodicals, Inc.

Bent Phosphaallenes With "Hidden" Lone Pairs as Ligands

Phosphaheteroallenes R-P=C=L, with L = N-heterocyclic carbenes (NHCs), can be viewed to a certain extent as phosphaisonitriles stabilized with NHCs, R-P=C:←L. The suitability of these molecules as ligands for coinage-metal ions was investigated and coordination through the central carbon center was observed in most cases. A combination of experiments, spectroscopic methods, and DFT calculations indicates the presence of a hidden electron pair at the carbon center of R-P=C:←L. Remarkably, this lone pair also inserts intramolecularly in C-H bonds showing the carbene-type reactivity which is expected for phosphaisonitriles.

Pressure-Induced Polymerization and Electrical Conductivity of a Polyiodide

We report the high-pressure structural characterization of an organic polyiodide salt in which a progressive addition of iodine to triiodide groups occurs. Compression leads to the initial formation of discrete heptaiodide units, followed by polymerization to a 3D anionic network. Although the structural changes appear to be continuous, the insulating salt becomes a semiconducting polymer above 10 GPa. The features of the pre-reactive state and the polymerized state are revealed by analysis of the computed electron and energy densities. The unusually high electrical conductivity can be explained with the formation of new bonds.

Octacarbonyl Anion Complexes of the Late Lanthanides Ln(CO)8 - (Ln=Tm, Yb, Lu) and the 32-Electron Rule

The lanthanide octacarbonyl anion complexes Ln(CO)₈ ^- (Ln=Tm, Yb, Lu) were produced in the gas phase and detected by mass-selected infrared photodissociation spectroscopy in the carbonyl stretching-frequency region. By comparison of the experimental CO-stretching frequencies with calculated data, which are strongly red-shifted with respect to free CO, the Yb(CO)₈ ^- and Lu(CO)₈ ^- complexes were determined to possess octahedral (O_h ) symmetry and a doublet X² A_2u (Yb) and singlet X¹ A_1g (Lu) electronic ground state, whereas Tm(CO)₈ ^- exhibits a D_4h equilibrium geometry and a triplet X³ B_1g ground state. The analysis of the electronic structures revealed that the metal-CO attractive forces come mainly from covalent orbital interactions, which are dominated by [Ln(d)]→(CO)₈ π backdonation and [Ln(d)]←(CO)₈ σ donation (contributes ≈77 and 16 % to covalent bonding, respectively). The metal f orbitals play a very minor role in the bonding. The electronic structure of all three lanthanide complexes obeys the 32-electron rule if only those electrons that occupy the valence orbitals of the metal are considered.

Dispersion-mediated steering of organic adsorbates on a precovered silicon surface

The chemistry of organic adsorbates on surfaces is often discussed in terms of Pauli repulsion as limiting factor regarding the packing of molecules. Here we show that the attractive part of the van der Waals potential can be similarly decisive. For the semiconductor surface Si(001), an already covalently bonded molecule of cyclooctyne steers a second incoming molecule via dispersion interactions onto the neighbouring adsorption site. This helps in understanding the nonstatistical pattern formation for this surface–adsorbate system and hints toward an inclusion of dispersion attraction as another determining factor for surface adsorption.

RNA versus DNA G-Quadruplex: The Origin of Increased Stability

DNA quadruplexes have been the subject of investigation because of their biological relevance and because of their potential application in supramolecular chemistry. Similarly, RNA quadruplexes are now gaining increasing attention. Although DNA and RNA quadruplexes are structurally very similar, the latter show higher stability. In this study we report dispersion‐corrected density functional theory (DFT‐D) quantum chemical calculations that were undertaken to understand the difference in stabilities of RNA and DNA quadruplexes. The smallest meaningful model of a stack of quartets, interacting with alkali metal cations, was simulated in an aqueous environment. The energy decomposition analysis allows for in‐depth examination of the interaction energies, emphasising the role of noncovalent interactions and better electrostatics in determining RNA‐GQs higher stabilities, particularly pinpointing the role of the extra 2′‐OH groups. Furthermore, our computations present new insights on why the cation is required for self‐assembly: unexpectedly the cation is not necessary to relieve the repulsion between the oxygen atoms in the central cavity, but it is needed to overcome the entropic penalty.

Covalency and Ionicity Do Not Oppose Each Other-Relationship Between Si-O Bond Character and Basicity of Siloxanes

Covalency and ionicity are orthogonal rather than antipodal concepts. We demonstrate for the case of siloxane systems [R₃ Si-(O-SiR₂ )_n -O-SiR₃ ] that both covalency and ionicity of the Si-O bonds impact on the basicity of the Si-O-Si linkage. The relationship between the siloxane basicity and the Si-O bond character has been under debate since previous studies have presented conflicting explanations. It has been shown with natural bond orbital methods that increased hyperconjugative interactions of LP(O)→σ*(Si-R) type, that is, increased orbital overlap and hence covalency, are responsible for the low siloxane basicity at large Si-O-Si angles. On the other hand, increased ionicity towards larger Si-O-Si angles has been revealed with real-space bonding indicators. To resolve this ostensible contradiction, we perform a complementary bonding analysis, which combines orbital-space, real-space, and bond-index considerations. We analyze the isolated disiloxane molecule H₃ SiOSiH₃ with varying Si-O-Si angles, and n-membered cyclic siloxane systems Si₂ H₄ O(CH₂ )_n-3 . All methods from quite different realms show that both covalent and ionic interactions increase simultaneously towards larger Si-O-Si angles. In addition, we present highly accurate absolute hydrogen-bond interaction energies of the investigated siloxane molecules with water and silanol as donors. It is found that intermolecular hydrogen bonding is significant at small Si-O-Si angles and weakens as the Si-O-Si angle increases until no stable hydrogen-bond complexes are obtained beyond φ_SiOSi =168°, angles typically displayed by minerals or polymers. The maximum hydrogen-bond interaction energy, which is obtained at an angle of 105°, is 11.05 kJ mol^-1 for the siloxane-water complex and 18.40 kJ mol^-1 for the siloxane-silanol complex.

Response to "Covalent Bonding and Charge Shift Bonds: Comment on 'The Carbon-Nitrogen Bonds in Ammonium Compounds Are Charge Shift Bonds'"

A response to the comment by Gernot Frenking, outlining common ground, as well as differences, with regard to a recent paper on charge-shift bonding in quaternary ammonium cations.

Covalent Bonding and Charge Shift Bonds: Comment on "The Carbon-Nitrogen Bonds in Ammonium Compounds Are Charge Shift Bonds"

The paper by Gershoni-Poranne and Chen (R. Gershoni-Poranne, P. Chen, Chem. Eur. J. 2017, 23, 4659) gives an incorrect definition of covalent bonding. Furthermore, the assignment of so-called charge shift bonds in ammonium compounds has no physical foundation and is conceptually redundant.

Hydrogen-Bond Strength of CC and GG Pairs Determined by Steric Repulsion: Electrostatics and Charge Transfer Overruled

Theoretical and experimental studies have elucidated the bonding mechanism in hydrogen bonds as an electrostatic interaction, which also exhibits considerable stabilization by charge transfer, polarization, and dispersion interactions. Therefore, these components have been used to rationalize the differences in strength of hydrogen‐bonded systems. A completely new viewpoint is presented, in which the Pauli (steric) repulsion controls the mechanism of hydrogen bonding. Quantum chemical computations on the mismatched DNA base pairs CC and GG (C=cytosine, G=guanine) show that the enhanced stabilization and shorter distance of GG is determined entirely by the difference in the Pauli repulsion, which is significantly less repulsive for GG than for CC. This is the first time that evidence is presented for the Pauli repulsion as decisive factor in relative hydrogen‐bond strengths and lengths.

The Bonding Situation in Metalated Ylides

Quantum chemical calculations have been carried out to study the electronic structure of metalated ylides particularly in comparison to their neutral analogues, the bisylides. A series of compounds of the general composition Ph3 P-C-L with L being either a neutral or an anionic ligand were analyzed and the impact of the nature of the substituent L and the total charge on the electronics and bonding situation was studied. The charge at the carbon atom as well as the dissociation energies, bond lengths, and Wiberg bond indices strongly depend on the nature of L. Here, not only the charge of the ligand but also the position of the charge within the ligand backbone plays an important role. Independent of the substitution pattern, the NBO analysis reveals the preference of unsymmetrical bonding situations (P=C-L or P-C=L) for almost all compounds. However, Lewis structures with two lone-pair orbitals at the central carbon atom are equally valid for the description of the bonding situation. This is confirmed by the pronounced lone-pair character of the frontier orbitals. Energy decomposition analysis mostly reveals the preference of several bonding situations, mostly with dative and ylidic electron-sharing bonds (e.g., P→C- -L). In general, the anionic systems show a higher preference of the ylidic bonding situations compared to the neutral analogues. However, in most of the cases different resonance structures have to be considered for the description of the "real" bonding situation.

Carbones as Ligands in Novel Main-Group Compounds E[C(NHC)2 ]2 (E=Be, B+ , C2+ , N3+ , Mg, Al+ , Si2+ , P3+ ): A Theoretical Study

Quantum chemical calculations of the main-group compounds E[C(NHC^Me )₂ ]₂ (E=Be, B⁺ , C²⁺ , N³⁺ , Mg, Al⁺ , Si²⁺ , P³⁺ ) have been carried out using density functional theory at the BP86/def2-TZVPP and BP86-D3(BJ)/def2-TZVPP levels of theory. The geometry optimization at BP86/def2-TZVPP gives equilibrium structures with two-coordinated species E and bending angles C-E-C between 152.5° (E=Be) and 110.5° (E=Al). Inclusion of dispersion forces at BP86-D3(BJ)/def2-TZVPP yields a three-coordinated beryllium compound Be[C(NHC^Me )₂ ]₂ as the only energy minimum form. Three-coordinated isomers are found besides the two-coordinated energy minima for the boron and carbon cations B[C(NHC^Me )₂ ]₂⁺ and C[C(NHC^Me )₂ ]₂²⁺ . The three-coordinated form of the boron compound is energetically lower lying than the two-coordinated form, while the opposite trend is calculated for the carbon species. The theoretically predicted bond dissociation energies suggest that all compounds are viable species for experimental studies. The X-ray structure of the benzoannelated homologue of P[C(NHC^Me )₂ ]₂³⁺ that was recently reported by Dordevic et al. agrees quite well with the calculated geometry of the molecule. A detailed bonding analysis using charge and energy decomposition methods shows that the two-coordinated neutral compounds Be[C(NHC^Me )₂ ]₂ and Mg[C(NHC^Me )₂ ]₂ possess strongly positively charged atoms Be and Mg. The carbodicarbene groups C(NHC^Me )₂ serve as acceptor ligands in the compounds and may be sketched with dative bonds (NHC^Me )₂ C←E→C(NHC^Me )₂ (E=Be, Mg). Dative bonds in which the carbones C(NHC^Me )₂ are donor ligands are suggested for the cations (NHC^Me )₂ C→E←C(NHC^Me )₂ (E=B⁺ , Al⁺ ). The dications and trications possess electron-sharing bonds in which the bonding situation is best described with the formula [(NHC^Me )₂ C]⁺ -E-[C(NHC^Me )₂ ]⁺ (E=C, Si, N⁺ , P⁺ ).

NHC-Stabilised Acetylene-How Far Can the Analogy Be Pushed?

Experimental studies suggest that the compound (NHC^bz )₂ C₂ H₂ can be considered as a complex of a distorted acetylene fragment which is stabilised by benzoannelated N-heterocyclic carbene ligands (NHC^bz )→(C₂ H₂ )←(NHC^bz ). A quantum chemical analysis of the electronic structures shows that the description with dative bonds is more favourable than with electron-sharing double bonds (NHC^bz )=(C₂ H₂ )=(NHC^bz ).