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

Classical versus Collective Interactions in Asymmetric Trigonal Bipyramidal Alkaline Metal-Boron Halide Complexes

Collective interactions are a novel type of chemical bond formed between metals and electron-rich substituents around an electron-poor central atom. So far only a limited number of candidates for having collective interactions are reported. In this work, we extend the newly introduced concept of collective bonding to a series of neutral boron complexes with the general formula M2BX3 (M=Li, Na, and K; X=F, Cl, and Br). Our state-of-the-art ab initio computations suggest that these complexes form trigonal bipyramidal structures with a D3h to C3v distortion along the C3 axis of symmetry. The BX3 unit in the complexes distorts from planar to pyramidal akin to a sp3 hybridized atom. Interestingly, the interaction of the metals with the pyramidal side of BX3, where the lone pair in a hypothetical [BX3]2- should be located, is weaker than the interactions of metals with the inverted side, i. e., the middle of three halogen atoms. The origin of this stronger interaction can be explained by the formation of collective interactions between metals and halogen atoms as we explored via energy decomposition within the context of the theory of interacting quantum atoms, IQA.

BeM(CO)3- (M = Co, Rh, Ir) and BeM(CO)3 (M = Ni, Pd, Pt): Triply bonded terminal beryllium in zero oxidation state

We perform detailed potential energy surface explorations of BeM(CO)3- (M = Co, Rh, Ir) and BeM(CO)3 (M = Ni, Pd, Pt) using both single-reference and multireference-based methods. The present results at the CASPT2(12,12)/def2-QZVPD//M06-D3/def2-TZVPPD level reveal that the global minimum of BeM(CO)3- (M = Co, Rh, Ir) and BePt(CO)3 is a C3v symmetric structure with an 1A1 electronic state, where Be is located in a terminal position bonded to M along the center axis. For other cases, the C3v symmetric structure is a low-lying local minimum. Although the present complexes are isoelectronic with the recently reported BFe(CO)3- complex having a B-Fe quadruple bond, radial orbital-energy slope (ROS) analysis reveals that the highest occupied molecular orbital (HOMO) in the title complexes is slightly antibonding in nature, which bars a quadruple bonding assignment. Similar weak antibonding nature of HOMO in the previously reported BeM(CO)4 (M = Ru, Os) complexes is also noted in ROS analysis. The bonding analysis through energy decomposition analysis in combination with the natural orbital for chemical valence shows that the bonding between Be and M(CO)3q (q = -1 for M = Co, Rh, Ir and q = 0 for M = Ni, Pd, Pt) can be best described as Be in the ground state (1S) interacting with M(CO)30/- via dative bonds. The Be(spσ) → M(CO)3q σ-donation and the complementary Be(spσ) ← M(CO)3q σ-back donation make the overall σ bond, which is accompanied by two weak Be(pπ) ← M(CO)3q π-bonds. These complexes represent triply bonded terminal beryllium in an unusual zero oxidation state.

In Silico Design and Characterization of a New Molecular Electride: Li@Calix[3]Pyrrole

Electrides, in which anionic electrons are localized independently of the atoms in the compound, have shown promise, especially as catalysts and optoelectronic materials. Here, we present a new computationally designed molecular electride, Li@calix[3]pyrrole (Li@C3P). Electron density and electron localization function analyses unequivocally confirm the existence of localized electride electron density, outside the system, independent of any specific atoms. Non-covalent interaction plots further validate the character of the isolated localized electron, suggesting that the system can be accurately represented by Li+@calix[3]pyrrole ⋅ e-, denoting its distinct charge separation. The remarkable non-linear optical properties of Li@C3P, including average polarizability,α ‾ ${\bar{\alpha }}$ =412.4 au, first hyperpolarizability, β=4.46×104 au, and second hyperpolarizability,γ ∥ ${{\gamma }_{\parallel }}$ =18.40×106 au, are unparalleled in the previously reported and similar Li@C4P molecular electride. Furthermore, energy decomposition analysis in combination with natural orbital for chemical valence theory sheds light on the mechanism of electron density transfer from Li to the C3P cage, yielding the charge-separated Li@C3P complex. In addition to the electron transfer, a key factor to its electride nature is the electronic structure of the CnP cage, which has its lowest unoccupied molecular orbital located in the void adjacent to the N-H groups at the back of the bowl-shaped CnP cage.

MB16 - (M=Sc, Y, La): Perfect Bowl-Like Boron Clusters

The introduction of transition-metal doping has engendered a remarkable array of unprecedented boron motifs characterized by distinctive geometries and bonding, particularly those heretofore unobserved in pure boron clusters. In this study, we present a perfect (no defects) boron framework manifesting an inherently high-symmetry, bowl-like architecture, denoted as MB16 - (M=Sc, Y, La). In MB16 -, the B16 is coordinated to M atoms along the C5v-symmetry axis. The bowl-shaped MB16 - structure is predicted to be the lowest-energy structure with superior stability, owing to its concentric (2 π+10 π) dual π aromaticity. Notably, the C5v-symmetry bowl-like B16 - is profoundly stabilized through the doping of an M atom, facilitated by strong d-pπ interactions between M and boron motifs, in conjunction with additional electrostatic stabilization by an electron transfer from M to the boron motifs. This concerted interplay of covalent and electrostatic interactions between M and bowl-like B16 renders MB16 - a species of exceptional thermodynamic stability, thus making it a viable candidate for gas-phase experimental detection.

Preparation of a Compound with a SiII -SiIV -SiII Bonding Arrangement

Herein, we report the synthesis of a rare bis-silylene, 1, in which two SiII atoms are bridged by a SiIV atom. Compound 1 contains an unusual SiII -SiIV -SiII bonding arrangement with SiII -SiIV bond distances of 2.4212(8) and 2.4157(7) Å. Treatment of 1 with Fe(CO)5 afforded a dinuclear Fe0 complex 2 with two unusually long Si-Si bonds (2.4515(8) and 2.4488(10) Å). We have also carried out a detailed computational study to understand the nature of the Si-Si bonds in these compounds. Natural bond orbital (NBO) and energy decomposition analysis-natural orbital for chemical valence (EDA-NOCV) analyses reveal that the Si-Si bonds in 1 and 2 are of an electron-sharing nature.

Unveiling the Anomaly of Reduction of Carborane-bis-silylene-Stabilised Silylone/Germylone Leading to Unusual Oxidation of Si0 /Ge0 to SiI /GeI with EDA-NOCV Analyses

Researchers have successfully isolated Si0 /Ge0 species, termed silylone and germylone, with two lone pairs of electrons on them. These elusive compounds have been stabilised in singlet ground states by using different donor base ligands. Driess et al. in particular have made strides in this area, isolating carborane-bis-silylene-stabilised silylone/germylone and their N+ /Pb analogues. Carborane (C2 B10 H10 ) plays a pivotal role as a redox-active ligand, converting from closo-carborane to nido-carborane with the addition of two electrons. Notably, anomalous oxidation of Si0 /Ge0 centres in carborane-bis-silylene-stabilised species to SiI /GeI has been reported, resulting in the formation of dimeric SiI -SiI /GeI -GeI di-cationic units. The energy decomposition analysis coupled with natural orbital for chemical valence (EDA-NOCV) study focuses on the carborane-bis-silylene ligand in the free state, and its three other species, including silylone/germylone species. Interestingly, it reveals that the carborane unit in an anionic doublet state tends to form one electron-sharing bond and one dative bond with the counter fragment in its cationic doublet state. This helps us to rationalise why the carborane unit undergoes intramolecular electronic rearrangements leading to the formation of a di-anionic carborane unit with a significantly elongated C-C bond (2.38-2.68 Å) and undergoes unusual oxidation of Si0 /Ge0 to SiI /GeI .

Anti-electrostatic cation⋯π-hole and cation⋯lp-hole interactions are stabilized via collective interactions

Collective interactions are a novel type of bond between metals and AX3 fragments with an electropositive central atom, A, and electronegative X substituents. Here, using electrostatic potential maps and state-of-the-art bonding analysis tools we have shown that collective interactions are anti-electrostatic cation⋯π-Hole or cation⋯lp-Hole interactions.

Open shell versus closed shell bonding interaction in cyclopropane derivatives: EDA-NOCV analyses

Cyclopropane ring is a very common motif in organic/bio-organic compounds. The chemical bonding of this strained ring is taught to all chemistry students. This three-membered cyclic, C3 ring is quite reactive which has attracted both, synthetic and theoretical chemists to rationalize/correlate its stability and bonding with its reactivity and physical properties over a century. There are a few bonding models (mainly the Bent-Bond model and Walsh model) of this C3 ring that are debated to date. Herein, we have carried out energy decomposition analysis coupled with natural orbital for chemical valence (EDA-NOCV) to study the two most reactive bonds of cyclopropane rings of 49 different organic compounds containing different functional groups to obtain a much deeper bonding insight toward a more general bonding model of this class of compounds. The EDA-NOCV analyses of fragment orbitals and susequent bond formation revealed that the nature of the CC bond of the cyclopropane (splitting two bonds at a time out of three CC bonds) ring is preferred to form two dative covalent CC bonds (between a singlet olefin-fragment and an excited singlet carbene-fragment with a vacant sp2 orbital and a filled p-orbital) for the majority (37/49) of compounds over two covalent electron sharing bonds in some (7/49) compounds (between an excited triplet olefin and triplet carbene), while a few (5/49) compounds show flexibility to adopt either the electron sharing or dative covalent bond as both are equally possible. The effects of functional groups on the nature of chemical bond in cyclopropane rings have been studied in detail. Our bonding analyses are in line with the QTAIM analyses which produce small negative values of the Laplacian, significantly positive values of bond ellipticity, and accumulation of electron densities around the ring critical point of C3 -rings. These corresponding QTAIM parameters of C3 -rings are quite different for CC single bonds of normal hydrocarbons as expected. The chemical bonding in the majority of cyclopropane rings can be very similar to those of metal-olefin systems.

Inter-Alkali-Metal Dative Bond in the MMN3 - (M=Alkali Metal) Cluster

Alkali metals are generally Lewis acids. On the contrary, Lewis basic character of alkali metals forming donor - acceptor complexes is a very rare phenomenon. In this contribution, I have theoretically designed an anionic cluster MMN₃ ^- (M=alkali metals) on the basis of experimentally known reagent, alkali salt of azide ion MN₃ , which shows unprecedented M:^- →M donor-acceptor interaction. To the best of author's knowledge, the characterization of such donor-acceptor interaction among alkali metals is unprecedented. Formation of the 2c-2e donor-acceptor bonds have been confirmed by quantum theory of atoms in molecules and electron localization function analyses. The calculated bond dissociation energies are significant suggesting their possible spectroscopic identification.

Stability and bonding of carbon(0)-iron-N2 complexes relevant to nitrogenase co-factor: EDA-NOCV analyses

The factors/structural features which are responsible for the binding, activation and reduction of N₂ to NH₃ by FeMoco of nitrogenase have not been completely understood well. Several relevant model complexes by Holland et al. and Peters et al. have been synthesized, characterized and studied by theoretical calculations. For a matter of fact, those complexes are much different than real active N₂ -binding Fe-sites of FeMoco, which possesses a central C(4-) ion having an eight valence electrons as an μ₆ -bridge. Here, a series of [(S₃ C(0))Fe(II/I/0)-N₂ ]^n- complexes in different charged/spin states containing a coordinated σ- and π-donor C(0)-atom which possesses eight outer shell electrons [carbone, (Ph₃ P)₂ C(0); Ph₃ P→C(0)←PPh₃ ] and three S-donor sites (i.e. ^- S-Ar), have been studied by DFT, QTAIM, and EDA-NOCV calculations. The effect of the weak field ligand on Fe-centres and the subsequent N₂ -binding has been studied by EDA-NOCV analysis. The role of the oxidation state of Fe and N₂ -binding in different charged and spin states of the complex have been investigated by EDA-NOCV analyses. The intrinsic interaction energies of the Fe-N₂ bond are in the range from -42/-35 to -67 kcal/mol in their corresponding ground states. The S₃ C(0) donor set is argued here to be closer to the actual coordination environment of one of the six Fe-centres of nitrogenase. In comparison, the captivating model complexes reported by Holland et al. and Peter et al. possess a stronger π-acceptor C-ring (S₂ C_ring donor, π-C donor) and stronger donor set like CP₃ (σ-C donor) ligands, respectively.

Designing a Four-Ring Tubular Boron Motif through Metal Doping

Tubular boron clusters represent a class of extremely unusual geometries that can be regarded as a key indicator for the 2D-to-3D boron structural evolution as well as the embryos for boron nanotubes. While a good number of pure boron or metal-doped boron tubular clusters have been reported so far, most of them are two-ring tubular structures, and their higher-ring analogues are very scarce. We report herein the first example of a four-ring tubular boron motif in the cagelike global minimum of Be₂B₂₄⁺. Global-minimum searches of MB₂₄^q and M₂B₂₄^q (M = alkali/alkaline-earth metals; q = 1+, 0, 1-) reveal that the most stable structure of Be₂B₂₄⁺ is a C_2v-symmetric cage having a four-ring tubular boron moiety, whereas it is a high-lying isomer for those having a two/three-ring tubular boron motif for all other systems. The B₂₄ framework in Be₂B₂₄⁺ can be viewed as consisting of two two-ring B₁₂ tubular structures linked together at one side of the B₆ rings along the high-symmetry axis and two offside B₆ rings capped by two Be atoms. The Be₂-B₂₄ bonding is best described as Be₂²⁺ in an excited triplet state, forming two highly polarized covalent bonds with B₂₄^- in a quartet spin state. The unique ability of beryllium to make strong covalent and electrostatic interactions makes the Be₂B₂₄⁺ cluster stable in such an unusual geometry.

The Labile Nature of Air Stable Ni(II)/Ni(0)-phosphine/Olefin Catalysts/Intermediates: EDA-NOCV Analysis

Metal ions-based inorganic-organic hybrid composites are often reported acting as good to excellent catalysts with various substrate scopes under milder reaction conditions. The active catalyst of a catalytic cycle is sometimes proposed to be a short-lived reactive intermediate species. A three coordinate (L-Me)Ni(II) intermediate species [L-Me = O 2 N donor dianionic ligand] can bind with short-lived carbene-ester ligands to produce four coordinate Ni(II) species which can act as carbene transfer intermediates under suitable reaction conditions for C-H functionalization or cyclopropanation reactions. The dissociation of phosphine (PPh 3 ) from the Ni(II) centre of (L-Me)Ni(II)(PPh 3 ) ( 1a ) and binding of short lived carbene esters (:CR 1 -CO 2 R 2 ; R 1 = H, Ph; R 2 = aliphatic group; 2-4 and other carbenes; 5-10 ) to Ni(II) rationalize the phenomenon in solution. Air stable Ni(0)-olefin complexes/intermediates ( 12-18 ) have recently been shown to mediate a variety of organic transformations. This analysis will further help organic/organometallic chemists to rationalize the design and synthesis of future catalysts for organic transformation. EDA-NOCV calculations have been performed to shed light on the stability and bonding of those species. Additionally, our analysis provides a proper reason why the analogous (L-Me)Pd-PPh 3 complex ( 1b ) does not dissociate in solution and hence, a similar catalytic product has not been isolated from identical reaction conditions. The stability and the labile nature of Ni(II/0) complexes has been investigated by state-of-the-art EDA-NOCV analyses.

Uncovering the hidden reactivity of benzyne/aryne precursors utilized under milder condition: Bonding and stability studies by EDA-NOCV analyses

Arenes [C₆ H₃ R(TMS)(OTf); also called benzyne/aryne precursors] containing inter-related leaving groups Me₃ Si (TMS) and CF₃ SO₃ (OTf) on the adjacent positions (1,2-position) are generally converted to their corresponding aryne-intermediates via the addition of fluoride anion (F^- ) and subsequent elimination of TMS and OTf groups. This reaction is believed to proceed via the formation of an anionic intermediate [C₆ H₄ (TMS-F)(OTf)]^- . The EDA-NOCV analysis (EDA-NOCV = energy decomposition analysis-natural orbital for chemical valence) of over 35 such precursors of varied types have been reported to reveal bonding and stability of C_Ar Si and COTf bonds. EDA-NOCV showed that the nature of the C_Ar Si bond of C₆ H₃ R(TMS)(OTf) can be expressed as both dative and electron sharing [C_Ar Si, C_Ar →Si]. The C_Ar OTf bond, on the other hand, can be described explicitly as dative [C_Ar ←OTf]. The nature of C_Ar Si bond of [C₆ H₄ (TMS-F)(OTf)]^- exclusively changes to covalent dative σ-bond C_Ar →S(Me)3F on the attachment of F^- to the TMS group of C₆ H₄ (TMS)(OTf). Introduction of σ-electron withdrawing group (like OMe, NMe₂ , and NO₂ ) to the ortho-position of the TMS group of functionalized arynes C₆ H₃ R(TMS)(OTf) prefer to have a covalent dative σ-bond (C_Ar →Si) over an electron-sharing covalent σ-bond (C_Ar Si). If this σ-electron withdrawing group is shifted from ortho-position to meta- and para-positions, then the preference for a dative bond decreases significantly, implying that the electronic effect on the nature of chemical bonds affects through bond paths. This effect dies with distance, similar to the well-known inductive effect.

Can an alkalide act as a perfect Lewis base?

The Lewis basic character of alkali metals forming donor-acceptor complexes is a very rare phenomenon. No Lewis adduct with an alkalide as the Lewis basic centre has ever been reported. Herein, we theoretically designed EXH₂^- (E = Li, Na, K; X = Be, Mg, Ca) clusters which represent the first true example of Lewis adducts with alkalides as the two-electron donor basic sites. Our high level ab initio calculations reveal the formation of an unprecedented E:^- → XH₂ donor-acceptor interaction. Topological analysis within the realm of the electron localization function confirms this bonding scenario. The bonding scenario is exactly replicated in all the clusters, rendering support to our proposal. The calculated bond dissociation energies are significant, suggesting their possible spectroscopic identification.

Collective interactions among organometallics are exotic bonds hidden on lab shelves

Recent discovery of an unusual bond between Na and B in NaBH₃^- motivated us to look for potentially similar bonds, which remained unnoticed among systems isoelectronic with NaBH₃^-. Here, we report a novel family of collective interactions and a measure called exchange-correlation interaction collectivity index (ICI_XC; [Formula: see text]) to characterize the extent of collective versus pairwise bonding. Unlike conventional bonds in which ICI_XC remains close to one, in collective interactions ICI_XC may approach zero. We show that collective interactions are commonplace among widely used organometallics, as well as among boron and aluminum complexes with the general formula [M^a+AR₃]^b- (A: C, B or Al). In these species, the metal atom interacts more efficiently with the substituents (R) on the central atoms than the central atoms (A) upon forming efficient collective interactions. Furthermore, collective interactions were also found among fluorine atoms of XF_n systems (X: B or C). Some of organolithium and organomagnesium species have the lowest ICI_XC among the more than 100 studied systems revealing the fact that collective interactions are rather a rule than an exception among organometallic species.

Bonding analysis of the C2 precursor Me3E–C2–I(Ph)FBF3 (E = C, Si, Ge)

A series of possible precursors for generating C2 with the general formula Me3E–C2–I(Ph)FBF3 [E = C (1), Si (2), and Ge (3)] has been theoretically investigated using quantum chemical calculations. The equilibrium geometries of all species show a linear E–C2–I+ backbone. The inspection of the electronic structure of the Me3E–C2 bond by energy decomposition analysis coupled with the natural orbital for chemical valence (EDA-NOCV) method suggests a combination of electron sharing C–C σ-bond and v weak π-dative bond between Me3C and C2 fragments in the doublet state for species 1 (E = C). For species 2 (Si) and 3 (Ge), the analysis reveals σ-dative Me3E–C2 bonds (E = Si, Ge; Me3E←C2) resulting from the interaction of singly charged (Me3E)+ and (C2–IPh(BF4))− fragments in their singlet states. The C2–I bond is diagnosed as an electron sharing σ-bond in all three species, 1, 2 and 3.

Bonding and stability of dinitrogen-bonded donor base-stabilized Si(0)/Ge(0) species [(cAACMe-Si/Ge)2(N2)]: EDA-NOCV analysis

Recently, dinitrogen (N₂) binding and its activation have been achieved by non-metal compounds like intermediate cAAC-borylene as (cAAC)₂(B-Dur)₂(N₂) [cAAC = cyclic alkyl(amino) carbene; Dur = aryl group, 2,3,5,6-tetramethylphenyl; B-Dur = borylene]. It has attracted a lot of scientific attention from different research areas because of its future prospects as a potent species towards the metal free reduction of N₂ into ammonia (NH₃) under mild conditions. Two (cAAC)(B-Dur) units, each of which possesses six valence electrons around the B-centre, are shown to accept σ-donations from the N₂ ligand (B ← N₂). Two B-Dur further provide π-backdonations (B → N₂) to a central N₂ ligand to strengthen the B–N₂–B bond, providing maximum stability to the compound (cAAC)₂(B-Dur)₂(N₂) since the summation of each pair wise interaction accounted for the total stabilization energy of the molecule. (cAAC)(B-Dur) unit is isolobal to cAAC–E (E = Si, Ge) fragment. Herein, we report on the stability and bonding of cAAC–E bonded N₂-complex (cAAC–E)₂(N₂) (1–2; Si, Ge) by NBO, QTAIM and EDA-NOCV analyses (EDA-NOCV = energy decomposition analysis coupled with natural orbital for chemical valence; QTAIM = quantum theory of atoms in molecule). Our calculation suggested that syntheses of elusive (cAAC–E)₂(N₂) (1–2; Si, Ge) species may be possible with cAAC ligands having bulky substitutions adjacent to the C_cAAC atom by preventing the homo-dimerization of two (cAAC)(E) units which can lead to the formation of (cAAC–E)₂. The formation of E Created by potrace 1.16, written by Peter Selinger 2001-2019 ]]> E bond is thermodynamically more favourable (E = Si, Ge) over binding energy of N₂ inbetween two cAAC–E units.

Dinitrogen (N₂) binding and its activation have been achieved by non-metal compounds like intermediate cAACborylene with the general formula of (cAAC)₂(B-Dur)₂(N₂) [cAAC = cyclic alkyl(amino)carbene; Dur = aryl group, 2,3,5,6-tetramethylphenyl; B-Dur = aryl-borylene].

Path-dependency of energy decomposition analysis & the elusive nature of bonding

Here, we provide evidence of the path-dependency of the energy components of the energy decomposition analysis scheme, EDA, by studying a set of thirty-one closed-shell model systems with the D_2h symmetry point group. For each system, we computed EDA components from nine different pathways and numerically showed that the relative magnitudes of the components differ substantially from one path to the other. Not surprisingly, yet unfortunately, the most significant variations in the relative magnitudes of the EDA components appear in the case of species with bonds within the grey zone of covalency and ionicity. We further discussed that the role of anions and their effect on arbitrary Pauli repulsion energy components affects the nature of bonding defined by EDA. The outcome variation by the selected partitioning scheme of EDA might bring arbitrariness when a careful comparison is overlooked.

Here, we provide evidence of the path-dependency of the energy components of the energy decomposition analysis scheme, EDA, by studying a set of thirty-one closed-shell model systems with the D_2h symmetry point group.

Estimations of Fe0/−1–N2 interaction energies of iron(0)-dicarbene and its reduced analogue by EDA-NOCV analyses: crucial steps in dinitrogen activation under mild conditions

Metal complexes containing low valence iron atoms are often experimentally observed to bind with the dinitrogen (N₂) molecule. This phenomenon has attracted the attention of industrialists, chemists and bio-chemists since these N₂-bonded iron complexes can produce ammonia under suitable chemical or electrochemical conditions. The higher binding affinity of the Fe-atom towards N₂ is a bit ‘mysterious’ compared to that of the other first row transition metal atoms. Fine powders of α-Fe⁰ are even part of industrial ammonia production (Haber–Bosch process) which operates at high temperature and high pressure. Herein, we report the EDA-NOCV analyses of the previously reported dinitrogen-bonded neutral molecular complex (cAAC^R)₂Fe⁰–N₂ (1) and mono-anionic complex (cAAC^R)₂Fe⁻¹–N₂ (2) to give deeper insight of the Fe–N₂ interacting orbitals and corresponding pairwise intrinsic interaction energies (cAAC^R = cyclic alkyl(amino) carbene; R = Dipp or Me). The Fe⁰ atom of 1 prefers to accept electron densities from N₂via σ-donation while the comparatively electron rich Fe⁻¹ centre of 2 donates electron densities to N₂via π-backdonation. However, major stability due to the formation of an Fe–N₂ bond arises due to Fe → N₂ π-backdonation in both 1 and 2. The cAAC^R ligands act as a charge reservoir around the Fe centre. The electron densities drift away from cAAC ligands during the binding of N₂ molecules mostly via π-backdonation. EDA-NOCV analysis suggests that N₂ is a stronger π-acceptor rather than a σ-donor.

The stable Fe–N₂ bond of stable complex should have a sufficiently high interaction energy.

Na⋯B bond in NaBH 3 - : An induced spin-polarized bond

The nature of the Na⋯B bond, in the recently synthesized NaBH 3 - adduct, is analyzed on the light of the Na^- propensity to polarize along the bond axis as a consequence of the electric field produced by the BH₃ fragment. The observed induced polarization has two consequences: (i) the energetic stabilization of the Na^- , and (ii) the split of its valence electrons into two opposite lobes along the bond axis. Additionally, an analysis of the electron localization is presented using the information content of the correlated conditional pair density that reveals a significant delocalization between one lobe of the polarized Na^- anion and the BH₃ fragment at the equilibrium distance. Our findings reported here complement previous works on this system.

A Critical Look at Linus Pauling's Influence on the Understanding of Chemical Bonding

The influence of Linus Pauling on the understanding of chemical bonding is critically examined. Pauling deserves credit for presenting a connection between the quantum theoretical description of chemical bonding and Gilbert Lewis's classical bonding model of localized electron pair bonds for a wide range of chemistry. Using the concept of resonance that he introduced, he was able to present a consistent description of chemical bonding for molecules, metals, and ionic crystals which was used by many chemists and subsequently found its way into chemistry textbooks. However, his one-sided restriction to the valence bond method and his rejection of the molecular orbital approach hindered further development of chemical bonding theory for a while and his close association of the heuristic Lewis binding model with the quantum chemical VB approach led to misleading ideas until today.

Energy Decomposition Analysis Coupled with Natural Orbitals for Chemical Valence and Nucleus-Independent Chemical Shift Analysis of Bonding, Stability, and Aromaticity of Functionalized Fulvenes: A Bonding Insight

The Donor base ligand-stabilized cyclopentadienyl-carbene compounds L-C5H4 (L = H2C, aAAC; (CO2Me)2C, Py; aNHC, NHC, PPh3; SNHC; aAAC = acyclic alkyl(amino) carbene, aNHC = acyclic N-hetero cyclic carbene, NHC = cyclic N-hetero cyclic carbene, SNHC = saturated N-hetero cyclic carbene, Py = pyridine) (1a-1d, 2a-2c, 3) have been theoretically investigated by energy decomposition analysis coupled with natural orbitals for chemical valence calculation. Among all these compounds, aNHC=C5H4 (2a) and Ph3P=C5H4 (2c) had been reported five decades ago. The bonding analysis of compounds with the general formula L=C5H4 (1a-1d) [L = (H2C, aAAC, (CO2Me)2C, Py] showed that they possess one electron-sharing σ bond and electron-sharing π bond between L and C5H4 neutral fragments in their triplet states as expected. Interestingly, the bonding scenarios have completely changed for L = aNHC, NHC, PPh3, SNHC. The aNHC analogue (2a) prefers to form one electron-sharing σ bond (CL-CC5H4) and dative π bond (CL ← CC5H4) between cationic (aNHC)+ and anionic C5H4 - fragments in their doublet states. Similar bonding scenarios have been observed for NHC (2b) and PPh3 (2c) (PL-CC5H4, PL ← CC5H4) analogues. In contrast, the SNHC and C5H4 neutral fragments of SNHC=C5H4 (3) prefer to form a dative σ bond (CSNHC → CC5H4) and a dative π bond (CSNHC ← CC5H4) in their singlet states. The pyridine analogue 1d is quite different from 2c from the bonding and aromaticity point of view. The nucleus-independent chemical shifts of all the abovementioned species (1-3) corresponding to aromaticity have been computed using the gauge-independent atomic orbital approach.

Stabilization of group 14 elements E = C, Si, Ge by hetero-bileptic ligands cAAC, MCOn with push-pull mechanism

The stability and bonding of a series of hetero-diatomic molecules with general formula (cAAC)EM(CO)_n , where cAAC = cyclic alkyl(amino) carbene; E = group 14 elements (C, Si, and Ge); M = transition metal (Ni, Fe, and Cr) have been studied by quantum chemical calculations using density functional theory (DFT) and energy decomposition analysis-natural orbital chemical valence (EDA-NOCV). The equilibrium geometries were calculated at the BP86/def2-TZVPP level of theory. The tri-coordinated group 14 complex (1a, 4a, and 7a) in which one of the CO groups is migrated to the central group 14 element from adjacent metal is theoretically found to be more stable when the central atom (E) is carbon. On the other hand, the two-coordinate group 14 element containing metal-complexes (2, 5, 8, 3, 6, and 9) are found to be more stable with their corresponding heavier analogues. The electronic structures of all the molecules have been analyzed by molecular orbital, topological analysis of electron density and natural bond orbital (NBO) analysis at the M06/def2-TZVPP//BP86/def2-TZVPP level of theory. The nature of the cAACE and EM bonds has been studied by EDA-NOCV calculations at BP86-D3(BJ)/TZ2P level of theory. The EDA analysis suggests that the bonding of cAACC(CO) can be best represented by electron sharing σ and π interactions, whereas, C(CO)M(CO)_n-1 by dative σ and π interactions. On the other hand, EDA-NOCV calculations suggests both dative σ and π interactions for cAACE and EM(CO)_n bonds of the corresponding Si and Ge analogues having stronger σ- and relatively weaker π-bonds. The topological analysis of electron density supports the closed-shell interaction for the Si and Ge complexes and open-shell interaction for the carbon complexes. The calculated proton affinity and hydride affinity values corroborated well with the present bonding description. This class of complexes might act as efficient future catalysts for different organic transformations due to the presence of electron rich group 14 element and metal carbonyl.

Diboron Bonds Between BX3 (X=H, F, CH3 ) and BYZ2 (Y=H, F; Z=CO, N2 , CNH)

The ability of B atoms on two different molecules to engage with one another in a noncovalent diboron bond is studied by ab initio calculations. Due to electron donation from its substituents, the trivalent B atom of BYZ₂ (Z=CO, N₂ , and CNH; Y=H and F) has the ability to in turn donate charge to the B of a BX₃ molecule (X=H, F, and CH₃ ), thus forming a B⋅⋅⋅B diboron bond. These bonds are of two different strengths and character. BH(CO)₂ and BH(CNH)₂ , and their fluorosubstituted analogues BF(CO)₂ and BF(CNH)₂ , engage in a typical noncovalent bond with B(CH₃ )₃ and BF₃ , with interaction energies in the 3-8 kcal/mol range. Certain other combinations result in a much stronger diboron bond, in the 26-44 kcal/mol range, and with a high degree of covalent character. Bonds of this type occur when BH₃ is added to BH(CO)₂ , BH(CNH)₂ , BH(N₂ )₂ , and BF(CO)₂ , or in the complexes of BH(N₂ )₂ with B(CH₃ )₃ and BF₃ . The weaker noncovalent bonds are held together by roughly equal electrostatic and dispersion components, complemented by smaller polarization energy, while polarization is primarily responsible for the stronger ones.

Bonding in M(NHBMe)2 and M[Mn(CO)5]2 complexes (M=Zn, Cd, Hg; NHBMe=(HCNMe)2B): divalent group 12 metals with zero oxidation state

Quantum chemical studies using density functional theory were carried out on M(NHBMe)2 and M[Mn(CO)5]2 (M=Zn, Cd, Hg) complexes. The calculations suggest that M(NHBMe)2 and M[Mn(CO)5]2 have D2d and D4d symmetry, respectively, with a 1A1 electronic ground state. The bond dissociation energies of the ligands have the order of Zn > Cd > Hg. A thorough bonding analysis using charge and energy decomposition methods suggests that the title complexes are best represented as NHBMe⇆M0⇄NHBMe and Mn(CO)5⇆M0⇄Mn(CO)5 where the metal atom M in the electronic ground state with an ns2 electron configuration is bonded to the (NHBMe)2 and [Mn(CO)5]2 ligands through donor–acceptor interaction. These experimentally known complexes are the first examples of mononuclear complexes with divalent group 12 metals with zero oxidation state that are stable at ambient condition. These complexes represent the rare situation where the ligands act as a strong acceptor and the metal center acts as strong donor. The relativistic effect of Hg leads to a weaker electron donating strength of the 6s orbital, which explains the trend of the bond dissociation energy.

Na⋅⋅⋅B Bond in NaBH3 - : Solving the Conundrum

Bonding in the recently synthesized NaBH₃ ^- cluster is investigated using the high level Valence Bond BOVB method. Contrary to earlier conclusions, the Na-B bond is found to be neither a genuine dative bond, nor a standard polar-covalent bond at equilibrium. It is rather revealed as a split and polarized weakly coupled electron-pair, which allows this cluster to be more effectively stabilized by a combination of (major) dipole-dipole electrostatic interaction and (secondary) resonant one-electron bonding mechanism. Our analysis of this unprecedented bonding situation extends to similar clusters, and the VB model unifies and articulates the previously published variegated views on this exotic "bond".

Neither too Classic nor too Exotic: One-Electron Na⋅B Bond in NaBH3 - Cluster

It is here reported that the NaBH₃ ^- cluster exhibits a Na⋅B one-electron bond, a well-established type of electron-deficient bonding in the literature. The topological analysis of the electron localization function, at the correlated level, reveals that Na^- , when approaching the bonding distance, fairly distributes its valence electron pair between two lobes. One of these electrons is used to bond with BH₃ , which participates through its boron empty p-orbital. Furthermore, the bonding situation of LiBH₃ ^- , KBH₃ ^- , MgBH₃ , and CaBH₃ global minima structures are similar to that of NaBH₃ ^- , extending the family of these new one-electron bond systems with biradicaloid character.

Metal-metal bonded alkaline-earth distannyls

The first families of alkaline-earth stannylides [Ae(SnPh3)2·(thf) x ] (Ae = Ca, x = 3, 1; Sr, x = 3, 2; Ba, x = 4, 3) and [Ae{Sn(SiMe3)3}2·(thf) x ] (Ae = Ca, x = 4, 4; Sr, x = 4, 5; Ba, x = 4, 6), where Ae is a large alkaline earth with direct Ae-Sn bonds, are presented. All complexes have been characterised by high-resolution solution NMR spectroscopy, including 119Sn NMR, and by X-ray diffraction crystallography. The molecular structures of [Ca(SnPh3)2·(thf)4] (1'), [Sr(SnPh3)2·(thf)4] (2'), [Ba(SnPh3)2·(thf)5] (3'), 4, 5 and [Ba{Sn(SiMe3)3}2·(thf)5] (6'), most of which crystallised as higher thf solvates than their parents 1-6, were established by XRD analysis; the experimentally determined Sn-Ae-Sn' angles lie in the range 158.10(3)-179.33(4)°. In a given series, the 119Sn NMR chemical shifts are slightly deshielded upon descending group 2 from Ca to Ba, while the silyl-substituted stannyls are much more shielded than the phenyl ones (δ 119Sn/ppm: 1', -133.4; 2', -123.6; 3', -95.5; 4, -856.8; 5, -848.2; 6', -792.7). The bonding and electronic properties of these complexes were also analysed by DFT calculations. The combined spectroscopic, crystallographic and computational analysis of these complexes provide some insight into the main features of these unique families of homoleptic complexes. A comprehensive DFT study (Wiberg bond index, QTAIM and energy decomposition analysis) points at a primarily ionic Ae-Sn bonding, with a small covalent contribution, in these series of complexes; the Sn-Ae-Sn' angle is associated with a flat energy potential surface around its minimum, consistent with the broad range of values determined by experimental and computational methods.

Bonding and stability of donor ligand-supported heavier analogues of cyanogen halides (L')PSi(X)(L)

Fluoro- and chloro-phosphasilynes [X-Si[triple bond, length as m-dash]P (X = F, Cl)] belong to a class of illusive chemical species which are expected to have Si[triple bond, length as m-dash]P multiple bonds. Theoretical investigations of the bonding and stability of the corresponding Lewis base-stabilized species (L')PSi(X)(L) [L' = cAAC^Me (cyclic alkyl(amino) carbene); L = cAAC^Me, NHC^Me (N-heterocyclic carbene), PMe₃, aAAC (acyclic alkyl(amino) carbene); X = Cl, F] have been studied using the energy decomposition analysis-natural orbitals for chemical valence (EDA-NOCV) method. The variation of the ligands (L) on the Si-atom leads to different bonding scenarios depending on their σ-donation and π-back acceptance properties. The ligands with higher lying HOMOs prefer profoundly different bonding scenarios than the ligands with lower lying HOMOs. The type of halogen (Cl or F) on the Si-atom was also found to have a significant influence on the overall bonding scenario. The reasonably higher value and endergonic nature of the dissociation energies along with the appreciable HOMO-LUMO energy gap may corroborate to the synthetic viability of the homo and heteroleptic ligand-stabilized elusive PSi(Cl/F) species in the laboratory.

Revisiting the Bonding Scenario of Two Donor Ligand Stabilized C2 Species

Quantum chemical calculations using density functional methods were performed for complexes of type L₂C₂ with L = NHC^Me (1), SNHC^Me (2) (S = saturated), cAAC^Me (3), and diamidocarbene (DAC^Me) (4). The equilibrium structures of 1-4 possess almost linear C₄ cores. A high thermochemical stability of the complexes with respect to dissociation, L₂C₂ → C₂ + 2L, is indicated by the large bond dissociation energy following the order 3 > 4 > 2 > 1. The results show that the use of SNHC^Me and DAC^Me as ligands is preferable over NHC^Me. The bonding analysis using charge and energy decomposition methods reveals that (cAAC^Me)₂C₂ and (DAC^Me)₂C₂ possess genuine cumulene C₄ moieties, which results from the electron-sharing bonding between quintet L₂ and quintet C₂ fragments. In contrast, the bonding in (NHC^Me)₂C₂ and (SNHC^Me)₂C₂ comes from a combination of dative and electron-sharing interactions between doublet L₂⁺ and doublet C₂^- fragments.

Beyond the Classical Electron-Sharing and Dative Bond Picture: Case of the Spin-Polarized Bond

Chemical bonds are traditionally assigned as electron-sharing or donor-acceptor/dative. External criteria such as the nature of the dissociation process, energy partitioning schemes, or quantum chemical topology are invoked to assess the bonding situation. However, for systems with marked multi-reference character, this binary categorization might not be precise enough to render the bonding properties. A third scenario can be foreseen: spin polarized bonds. To illustrate this, the case of a NaBH₃ ^- cluster is presented. According to the analysis NaBH₃ ^- exhibits a strong diradical character and cannot be classified as either electron-sharing or a dative bond. Elaborated upon are the common problems of popular bonding descriptions. Additionally, a simple model, based on the bond order and local spin indicators, which discriminates between all three bonding situations, is provided.

The Na⋅⋅⋅B Bond in NaBH3 - : A Different Type of Bond

A newly introduced Na-B bond in NaBH₃ ^- has been a challenge for the chemical bonding community. Here, a series of MBH₃ ^- (M=Li, Na, K) species and NaB(CN)₃ ^- are studied within the context of quantum chemical topology approaches. The analyses suggest that M-B interaction cannot be classified as an ordinary covalent, dative, or even simple ionic interaction. The interactions are controlled by coulombic forces between the metals and the substituents on boron, for example, H or CN, more than the direct M-B interaction. On the other hand, while the characteristics of the (3, -1) critical points of the bonds are comparable to weak hydrogen bonds, not covalent bonds, the metal and boron share a substantial sum of electrons. To the best of the author's knowledge, the characteristics of these bonds are unprecedented among known molecules. Considering all paradoxical properties of these bonds, they are herein described as ionic-enforced covalent bonds.

Intriguing structural, bonding and reactivity features in some beryllium containing complexes

We highlighted our contributions to Be chemistry which include bond-stretch isomerism in Be₃²⁻ species, Be complexes bound with noble gas, CO, and N₂, Be based nanorotors, and intriguing bonding situations in some Be complexes.