The Nature of the Chemical Bond Revisited: An Energy-Partitioning Analysis of Nonpolar Bonds

Revisiting the bonding evolution theory: a fresh perspective on the ammonia pyramidal inversion and bond dissociations in ethane and borazane

This work offers a comprehensive and fresh perspective on the bonding evolution theory (BET) framework, originally proposed by Silvi and collaborators [X. Krokidis, S. Noury and B. Silvi, Characterization of elementary chemical processes by catastrophe theory, J. Phys. Chem. A, 1997, 101, 7277-7282]. By underscoring Thom's foundational work, we identify the parametric function characterizing bonding events along a reaction pathway through a three-step sequence to establish such association rigorously, namely: (a) computing the determinant of the Hessian matrix at all potentially degenerate critical points, (b) computing the relative distance between these points, and (c) assigning the unfolding based on these computations and considering the maximum number of critical points for each unfolding. In-depth examination of the ammonia inversion and the dissociation of ethane and ammonia borane molecules yields a striking discovery: no elliptic umbilic flag is detected along the reactive coordinate for any of the systems, contradicting previous reports. Our findings indicate that the core mechanisms of these chemical reactions can be understood using only two folds, the simplest polynomial of Thom's theory, leading to considerable simplification. In contrast to previous reports, no signatures of the elliptic umbilic unfolding were detected in any of the systems examined. This finding dramatically simplifies the topological rationalization of electron rearrangements within the BET framework, opening new approaches for investigating complex reactions.

Frustrated Radical Pairs in Organic Synthesis

Frustrated radical pairs (FRPs) describe the phenomenon that two distinct radicals─which would otherwise annihilate each other to form a closed-shell covalent adduct─can coexist in solution, owing to steric repulsion or weak bonding association. FRPs are typically formed via spontaneous single-electron transfer between two sterically encumbered precursors─an oxidant and a reductant─under ambient conditions. The two components of a FRP exhibit orthogonal chemical properties and can often act in cooperativity to achieve interesting radical reactivities. Initially observed in the study of traditional frustrated Lewis pairs, FRPs have recently been shown to be capable of homolytically activating various chemical bonds. In this Perspective, we will discuss the discovery of FRPs, their fundamental reactivity in chemical bond activation, and recent developments of their use in synthetic organic chemistry, including in C-H bond functionalization. We anticipate that FRPs will provide new reaction strategies for solving challenging problems in modern organic synthesis.

Periodic Trends in the Stabilization of Actinyls in Their Higher Oxidation States Using Pyrrophen Ligands

Owing to the prominent existence and unique chemistry of actinyls, their complexation with suitable ligands is of significant interest. The complexation of high-valent actinyl moieties (An = U, Np, Pu and Am) with the acyclic sal-porphyrin analogue called "pyrrophen" (L₍₁₎) and its dimethyl derivative (L₍₂₎) with four nitrogen and two oxygen donor atoms was studied using relativistic density functional theory. Based on the periodic trends, the [U^VO₂-L₍₁₎/L₍₂₎]^1- complexes show shorter bond lengths and higher bond orders that increase across the series of pentavalent actinyl complexes mainly due to the localization of the 5f orbitals. Among the hexavalent complexes, the [U^VIO₂-L₍₁₎/L₍₂₎] complexes have the shortest bonds. Following the uranyl complex, due to the plutonium turn, the [Am^VIO₂-L₍₁₎/L₍₂₎] complexes exhibit comparable properties with those of the former. Charge analysis suggests the complexation to be facilitated through ligand-to-metal charge transfer (LMCT) mainly through σ donation. Thermodynamic feasibility of complexation was modeled using hydrated actinyl moieties in aqueous medium and was found to be spontaneous. The dimethylated pyrrophen (L₍₂₎) shows higher magnitudes of thermodynamic parameters indicating increased feasibility compared to the unsubstituted ligand (L₍₁₎). Energy decomposition analysis (EDA) along with extended transition-state-natural orbitals for chemical valence theory (ETS-NOCV) analysis shows that the dominant electrostatic contributions decrease across the series and are counteracted by Pauli repulsion. Slight but considerable covalency is provided to hexavalent actinyl complexes by orbital contributions; this was confirmed by molecular orbital (MO) analysis that suggests strong covalency in americyl (VI) complexes. In addition to the pentavalent and hexavalent actinyl moieties, heptavalent actinyl species of neptunyl, plutonyl, and americyl were studied. Beyond the influence of the charges, the geometric and electronic properties point to the stabilization of neptunyl (VII) in the pyrrophen ligand environment, while the others shift to a lower (+VI) and relatively stable OS on complexation.

Atoms-In-Molecules' Faces of Chemical Hardness by Conceptual Density Functional Theory

The chemical hardness concept and its realization within the conceptual density functional theory is approached with innovative perspectives, such as the electronegativity and hardness equalization of atoms in molecules connected with the softness kernel, in order to examine the structure-reactivity equalization ansatz between the electronic sharing index and the charge transfer either in the additive or geometrical mean picture of bonding. On the other hand, the maximum hardness principle presents a relation with the chemical stability of the hardness concept. In light of the inverse relation between hardness and polarizability, the minimum polarizability principle has been proposed. Additionally, this review includes important applications of the chemical hardness concept to solid-state chemistry. The mentioned applications support the validity of the electronic structure principles regarding chemical hardness and polarizability in solid-state chemistry.

Demarcating Noncovalent and Covalent Bond Territories: Imine-CO2 Complexes and Cooperative CO2 Capture

Chemical bond territory is rich with covalently bonded molecules wherein a strong bond is formed by equal or unequal sharing of a quantum of electrons. The noncovalent version of the bonding scenarios expands the chemical bonding territory to a weak domain wherein the interplay of electrostatic and π-effects, dipole-dipole, dipole-induced dipole, and induced dipole-induced dipole interactions, and hydrophobic effects occur. Here we study both the covalent and noncovalent interactive behavior of cyclic and acyclic imine-based functional molecules (XN) with CO₂. All parent XN systems preferred the formation of noncovalent (nc) complex XN···CO₂, while more saturated such systems (XN') produced both nc and covalent (c) complexes XN'⁺-(CO₂)^-. In all such cases, crossover from an nc to c complex is clearly demarcated with the identification of a transition state (ts). The complexes XN'···CO₂ and XN'⁺-(CO₂)^- are bond stretch isomers, and they define the weak and strong bonding territories, respectively, while the ts appears as the demarcation point of the two territories. Cluster formation of XN with CO₂ reinforces the interaction between them, and all become covalent clusters of general formula (XN⁺-(CO₂)^-)_n. The positive cooperativity associated with the NH···OC hydrogen bond formation between any two XN'⁺-(CO₂)^- units strengthened the N-C coordinate covalent bond and led to massive stabilization of the cluster. For instance, the stabilizing interaction between the XN unit with CO₂ is increased from 2-7 kcal/mol range in a monomer complex to 14-31 kcal/mol range for the octamer cluster (XN'⁺-(CO₂)^-)₈. The cooperativity effect compensates for the large reduction in the entropy of cluster formation. Several imine systems showed the exergonic formation of the cluster and are predicted as potential candidates for CO₂ capture and conversion.

The Nature of the Polar Covalent Bond

Quantum chemical calculations using density functional theory are reported for the diatomic molecules LiF, BeO, and BN. The nature of the interatomic interactions is analyzed with the EDA-NOCV method, and the results are critically discussed and compared with data from QTAIM, NBO and Mayer approaches. Polar bonds, like nonpolar bonds, are caused by the interference of wave functions, which lead to an accumulation of electronic charge in the bonding region. Polar bonds generally have a larger percentage of electrostatic bonding to the total attraction, but nonpolar bonds may also possess large contributions from Coulombic interaction. The term "ionic contribution" refers to VB structures and is misleading because it refers to separate fragments with negligible overlap that occur only in the solid state and in solution, not in a molecule. The EDA-NOCV method gives detailed information about the individual orbital contributions, which can nicely be identified by visual inspection of the associated deformation densities. It is very important, particularly for polar bonds to distinguish between the interatomic interactions of the final dissociation products after bond rupture and the interactions between the fragments in the eventually formed bond.

Dipolar repulsion in α-halocarbonyl compounds revisited

The concept of dipolar repulsion has been widely used to explain several phenomena in organic chemistry, including the conformational preferences of carbonyl compounds. This model, in which atoms and bonds are viewed as point charges and dipole moment vectors, respectively, is however oversimplified. To provide a causal model rooted in quantitative molecular orbital theory, we have analyzed the rotational isomerism of haloacetaldehydes OHC–CH₂X (X = F, Cl, Br, I), using relativistic density functional theory. We have found that the overall trend in the rotational energy profiles is set by the combined effects of Pauli repulsion (introducing a barrier around gauche that separates minima at syn and anti), orbital interactions (which can pull the anti minimum towards anticlinal to maximize hyperconjugation), and electrostatic interactions. Only for X = F, not for X = Cl–I, electrostatic interactions push the preference from syn to anti. Our bonding analyses show how this trend is related to the compact nature of F versus the more diffuse nature of the heavier halogens.

Beyond point charges! The point charge concept within dipolar repulsion model is valid for compact atoms like fluorine. This model breaks down for larger halogens, for which the electrostatic attraction between nuclei and charge densities dominates.

TWOE Code: An Efficient Tool for Explicit Partition of Coupled Cluster and Configuration Interaction Energies into Atomic and Diatomic Contributions

The efficient implementation of the TWOE program for evaluating the atomic and interatomic energy components at post-HF level was developed. The systematic convergence of these terms up to a near full-CI limit was performed for the first time for a series of coupled cluster methods: CCSD → CCSDT → CCSDTQ → CCSDTQP. A comparison with corresponding CI approaches (up to fifth excitation level) is additionally discussed. For a set of diatomic systems, it was demonstrated that, along with a full molecular energy convergence, all its components are also converged but with different patterns. It was found that not all components are decreased in their values at increasing computational rank. For instance, atomic energy parts are decreased while interatomic (interaction) energies are increased as the limiting level is approached. Two schemes were employed for atomic partition of molecules: the Baders approach and planes dissection. Influence of dynamical correlation effects on atomic energy components was analyzed in detail. Current TWOE implementation allows one, in principle, to work with any ab initio method providing the two-particle density matrix. It is believed that the developed program will be a useful tool for a real space energy decomposition that helps to reveal the most peculiar points in the structure of the total and correlation energies of a molecule.

The Valence Orbitals of the Alkaline-Earth Atoms

Quantum chemical calculations of the alkaline‐earth oxides, imides and dihydrides of the alkaline‐earth atoms (Ae=Be, Mg, Ca, Sr, Ba) and the calcium cluster Ca₆H₉[N(SiMe₃)₂]₃(pmdta)₃ (pmdta=N,N,N′,N′′,N′′‐pentamethyldiethylenetriamine) have been carried out by using density functional theory. Analysis of the electronic structures by charge and energy partitioning methods suggests that the valence orbitals of the lighter atoms Be and Mg are the (n)s and (n)p orbitals. In contrast, the valence orbitals of the heavier atoms Ca, Sr and Ba comprise the (n)s and (n−1)d orbitals. The alkaline‐earth metals Be and Mg build covalent bonds like typical main‐group elements, whereas Ca, Sr and Ba covalently bind like transition metals. The results not only shed new light on the covalent bonds of the heavier alkaline‐earth metals, but are also very important for understanding and designing experimental studies.

Distortion-Controlled Redshift of Organic Dye Molecules

It is shown, quantum chemically, how structural distortion of an aromatic dye molecule can be leveraged to rationally tune its optoelectronic properties. By using a quantitative Kohn-Sham molecular orbital (KS-MO) approach, in combination with time-dependent DFT (TD-DFT), the influence of various structural and electronic tuning parameters on the HOMO-LUMO gap of a benzenoid model dye have been investigated. These parameters include 1) out-of-plane bending of the aromatic core, 2) bending of the bridge with respect to the core, 3) the nature of the bridge itself, and 4) π-π stacking. The study reveals the coupling of multiple structural distortions as a function of bridge length and number of bridges in benzene to be chiefly responsible for a decreased HOMO-LUMO gap, and consequently, red-shifting of the absorption wavelength associated with the lowest singlet excitation (λ≈560 nm) in the model cyclophane systems. These physical insights together with a rational approach for tuning the oscillator strength were leveraged for the proof-of-concept design of an intense near-infrared (NIR) absorbing cyclophane dye at λ=785 nm. This design may contribute to a new class of distortion-controlled NIR absorbing organic dye molecules.

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.

Curly arrows, electron flow, and reaction mechanisms from the perspective of the bonding evolution theory

Despite the usefulness of curly arrows in chemistry, their relationship with real electron density flows is still imprecise, and even their direct connection to quantum chemistry is still controversial. The paradigmatic description - from first principles - of the mechanistic aspects of a given chemical process is based mainly on the relative energies and geometrical changes at the stationary points of the potential energy surface along the reaction pathway; however, it is not sufficient to describe chemical systems in terms of bonding aspects. Probing the electron density distribution during a chemical reaction can provide important insights, enabling us to understand and control chemical reactions. This aim has required an extension of the relationships between the concepts of traditional chemistry and those of quantum mechanics. Bonding evolution theory (BET), which combines the topological analysis of the electron localization function (ELF) and Thom's catastrophe theory (CT), provides a powerful method that offers insight into the molecular mechanism of chemical rearrangements. In agreement with the laws of physical and aspects of quantum theory, BET can be considered an appropriate tool to tackle chemical reactivity with a wide range of possible applications. In this work, BET is applied to address a long-standing problem: the ability to monitor the flow of electron density. BET analysis shows a connection between quantum mechanics and bond making/forming processes. Likewise, the present approach retrieves the classical curly arrows used to describe the rearrangements of chemical bonds and provides detailed physical grounds for this type of representation. We demonstrate this procedure using the test set of prototypical examples of thermal ring apertures, and the degenerated Cope rearrangement of semibullvalene.

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.

Theoretical study of azido gauche effect and its origin

The strength and origin of the azido gauche effect were studied by ab initio calculations and compared with the well-known fluorine gauche effect.

The Chemical Bond in C2

Quantum chemical calculations using the complete active space of the valence orbitals have been carried out for Hn CCHn (n=0-3) and N2. The quadratic force constants and the stretching potentials of Hn CCHn have been calculated at the CASSCF/cc-pVTZ level. The bond dissociation energies of the C-C bonds of C2 and HC≡CH were computed using explicitly correlated CASPT2-F12/cc-pVTZ-F12 wave functions. The bond dissociation energies and the force constants suggest that C2 has a weaker C-C bond than acetylene. The analysis of the CASSCF wavefunctions in conjunction with the effective bond orders of the multiple bonds shows that there are four bonding components in C2, while there are only three in acetylene and in N2. The bonding components in C2 consist of two weakly bonding σ bonds and two electron-sharing π bonds. The bonding situation in C2 can be described with the σ bonds in Be2 that are enforced by two π bonds. There is no single Lewis structure that adequately depicts the bonding situation in C2. The assignment of quadruple bonding in C2 is misleading, because the bond is weaker than the triple bond in HC≡CH.

Direct estimate of the internal π-donation to the carbene centre within N-heterocyclic carbenes and related molecules

Fifteen cyclic and acylic carbenes have been calculated with density functional theory at the BP86/def2-TZVPP level. The strength of the internal X→p(π) π-donation of heteroatoms and carbon which are bonded to the C(II) atom is estimated with the help of NBO calculations and with an energy decomposition analysis. The investigated molecules include N-heterocyclic carbenes (NHCs), the cyclic alkyl(amino)carbene (cAAC), mesoionic carbenes and ylide-stabilized carbenes. The bonding analysis suggests that the carbene centre in cAAC and in diamidocarbene have the weakest X→p(π) π-donation while mesoionic carbenes possess the strongest π-donation.

Description of Polar Chemical Bonds from the Quantum Mechanical Interference Perspective

The Generalized Product Function Energy Partitioning (GPF-EP) method has been applied to a set of molecules, AH (A = Li, Be, B, C, N, O, F), CO and LiF with quite different dipole moments, in order to investigate the role played by the quantum interference effect in the formation of polar chemical bonds. The calculations were carried out with GPF wave functions treating all the core electrons as a single Hartree-Fock group and the bonding electrons at the Generalized Valence Bond Perfect-Pairing (GVB-PP) level, with the cc-pVTZ basis set. The results of the energy partitioning into interference and quasi-classical contributions along the respective Potential Energy Surfaces (PES) show that the main contribution to the depth of the potential wells comes from the interference term, which is an indication that all the molecules mentioned above form typical covalent bonds. In all cases, the stabilization promoted by the interference term comes from the kinetic contribution, in agreement with previous results. The analysis of the effect of quantum interference on the electron density reveals that while polarization effects (quasi-classical) tend to displace electronic density from the most polarizable atom toward the less polarizable one, interference (quantum effects) counteracts by displacing electronic density to the bond region, giving rise to the right electronic density and dipole moment.

The boron-boron triple bond in NHC→B 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 1111111111111111111111111111111111 1111111111111111111111111111111111 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 1111111111111111111111111111111111 1111111111111111111111111111111111 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 1111111111111111111111111111111111 1111111111111111111111111111111111 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 0000000000000000000000000000000000 B←NHC

Thorough examination of the electronic structure of the compound B₂(NHC^Me)₂ provides convincing evidence for a B Created by potrace 1.16, written by Peter Selinger 2001-2019 ]]> B triple bond.

Quantum chemical calculations of the compound B₂(NHC^Me)₂ and a thorough examination of the electronic structure with an energy decomposition analysis provide strong evidence for the appearance of boron–boron triple bond character. This holds for the model compound and for the isolated diboryne B₂(NHC^R)₂ of Braunschweig which has an even slightly shorter B–B bond. The bonding situation in the molecule is best described in terms of NHC^Me→B₂←NHC^Me donor–acceptor interactions and concomitant π-backdonation NHC^Me←B₂→NHC^Me which weakens the B–B bond, but the essential features of a triple bond are preserved. An appropriate formula which depicts both interactions is the sketch NHC^Me⇄B Created by potrace 1.16, written by Peter Selinger 2001-2019 ]]> B⇄NHC^Me. Calculations of the stretching force constants F_BB which take molecules that have genuine single, double and triple bonds as references suggest that the effective bond order of B₂(NHC^Me)₂ has the value of 2.34. The suggestion by Köppe and Schnöckel that the strength of the boron–boron bond in B₂(NHC^H)₂ is only between a single and a double bond is repudiated. It misleadingly takes the force constant F_BB of OBBO as the reference value for a B–B single bond which ignores π bonding contributions. The alleged similarity between the B–O bonds in OBBO and the B–C bonds in B₂(NHC^Me)₂ is a mistaken application of the principle of isolable relationship.

Combined activation strain model and energy decomposition analysis methods: a new way to understand pericyclic reactions

The recently introduced activation strain model (ASM) has allowed us to gain more insight into the intimacies of different fundamental processes in chemistry. In combination with the energy decomposition analysis (EDA) method, we have nowadays a very useful tool to quantitatively understand the physical factors that govern the activation barriers of reactions within organic and organometallic chemistry. In this Perspective article, we present selected illustrative examples of the application of this method to pericyclic reactions (Diels-Alder and double group transfer reactions) to show that this methodology nicely complements other more traditional, widely used theoretical methods.

Ferrocene/fullerene hybrids showing large second-order nonlinear optical activities: impact of the cage unit size

Ferrocene/fullerene complexes through face-to-face fusion enjoy the merits of both ferrocene and fullerene due to their strong donor–acceptor interactions.

2-Adamantylidene and its heavier analogues: hyperconjugation versus lone pair stability and electrophilicity versus nucleophilicity

Too heavy to bend: The Cieplak-type hyperconjugative interaction mainly contributes to the stability of 2C and 2Si whereas the inertness of lone pairs is the major factor for 2Ge and 2Sn. These molecules (2X; X = C–Sn) can be classified as a class of ambiphilic compounds.

Energy decomposition analysis of gauche preference in 2-haloethanol, 2-haloethylamine (halogen = F, Cl), their protonated forms and anti preference in 1-chloro-2-fluoroethane

Small, electronegative elements contribute more electrostatic and orbital stabilization to the anti → gauche isomerization, and greater steric repulsion. The first and the latter actually oppose our traditional view of conformational equilibria.

Neutral noble gas compounds exhibiting a Xe-Xe bond: structure, stability and bonding situation

The structure and stability towards decomposition of eight novel noble gas compounds having a Xe-Xe bond, which have not been experimentally observed so far, have been studied computationally. In addition, the nature of the Xe-Xe interaction has been analysed by a combination of the most popular methods to study the bonding situation of molecules, i.e. Natural Bond Orbital, Atom in Molecules and Energy Decomposition Analysis methods. Two related series of compounds have been considered: HXeXeX (X = F to I) and RXeXeR' (R = halogen atom). Our calculations indicate that the replacement of the fluorine atom by a heavier group 17 congener in the HXeXeX series leads to a less stable compound, thus making more difficult its experimental observation. The same effect occurs in the RXeXeR' series, but these species are more kinetically protected against the decomposition reaction and therefore, their experimental detection is more likely.

Borylene complexes (BH)L2 and nitrogen cation complexes (N+)L2: isoelectronic homologues of carbones CL2

Quantum chemical calculations using DFT (BP86, M05-2X) and ab initio methods (CCSD(T), SCS-MP2) have been carried out on the borylene complexes (BH)L(2) and nitrogen cation complexes (N(+))L(2) with the ligands L=CO, N(2), PPh(3), NHC(Me), CAAC, and CAAC(model). The results are compared with those obtained for the isoelectronic carbones CL(2). The geometries and bond dissociation energies of the ligands, the proton affinities, and adducts with the Lewis acids BH(3) and AuCl were calculated. The nature of the bonding has been analyzed with charge and energy partitioning methods. The calculated borylene complexes (BH)L(2) have trigonal planar coordinated boron atoms which possess rather short B-L bonds. The calculated bond dissociation energies (BDEs) of the ligands for complexes where L is a carbene (NHC or CAAC) are very large (D(e) =141.6-177.3 kcal mol(-1)) which suggest that such species might become isolated in a condensed phase. The borylene complexes (BH)(PPh(3))(2) and (BH)(CO)(2) have intermediate bond strengths (D(e) =90.1 and 92.6 kcal mol(-1)). Substituted homologues with bulky groups at boron which protect the boron atom from electrophilic attack might also be stable enough to become isolated. The BDE of (BH)(N(2))(2) is much smaller (D(e) =31.9 kcal mol(-1)), but could become observable in a low-temperature matrix. The proton affinities of the borylene complexes are very large, particularly for the bulky adducts with L=PPh(3), NHC(Me), CAAC(model) and CAAC and thus, they are superbases. All (BH)L(2) molecules bind strongly AuCl either η(1) (L=N(2), PPh(3), NHC(Me), CAAC) or η(2) (L=CO, CAAC(model)). The BDEs of H(3)B-(BH)L(2) adducts which possess a hitherto unknown boron→boron donor-acceptor bond are smaller than for the AuCl complexes. The strongest bonded BH(3) adduct that might be isolable is (BH)(PPh(3))(2)-BH(3) (D(e) =36.2 kcal mol(-1)). The analysis of the bonding situation reveals that (BH)-L(2) bonding comes mainly from the orbital interactions which has three major contributions, that is, the donation from the symmetric (σ) and antisymmetric (π(||)) combination of the ligand lone-pair orbitals into the vacant MOs of BH L→(BH)←L and the L←(BH)→L π backdonation from the boron lone-pair orbital. The nitrogen cation complexes (N(+))L(2) have strongly bent L-N-L geometries, in which the calculated bending angle varies between 113.9° (L=N(2)) and 146.9° (L=CAAC). The BDEs for (N(+))L(2) are much larger than those of the borylene complexes. The carbene ligands NHC and CAAC but also the phosphane ligands PPh(3) bind very strongly between D(e) =358.4 kcal mol(-1) (L=PPh(3)) and D(e) =412.5 kcal mol(-1) (L=CAAC(model)). The proton affinities (PA) of (N(+))L(2) are much smaller and they bind AuCl and BH(3) less strongly compared with (BH)L(2). However, the PAs (N(+))L(2) for complexes with bulky ligands L are still between 139.9 kcal mol(-1) (L=CAAC(model)) and 168.5 kcal mol(-1) (L=CAAC). The analysis of the (N(+))-L(2) bonding situation reveals that the binding interactions come mainly from the L→(N(+))←L donation while L←(N(+) )→L π backdonation is rather weak.