Clarifying notes on the bonding analysis adopted by the energy decomposition analysis

Bismuth as a Z-Type Ligand: an Unsupported Pt-Bi Donor-Acceptor Interaction and its Umpolung by Reaction with H2

Establishing unprecedented types of bonding interactions is one of the fundamental challenges in synthetic chemistry, paving the way to new (electronic) structures, physicochemical properties, and reactivity. In this context, unsupported element-element interactions are particularly noteworthy since they offer pristine scientific information about the newly identified structural motif. Here we report the synthesis, isolation, and full characterization of the heterobimetallic Bi/Pt compound [Pt(PCy3)2(BiMe2)(SbF6)] (1), bearing the first unsupported transition metal→bismuth donor/acceptor interaction as its key structural motif. 1 is surprisingly robust, its electronic spectra are interpreted in a fully relativistic approach, and it reveals an unprecedented reactivity towards H2.

The strongest dative bond in main-group compounds. Theoretical study of OAeF- (Ae = Be-Ba)

Quantum chemical calculations of the anions OAeF- (Ae = Be-Ba) have been carried out using ab initio methods at the CCSD(T)/def2-TZVPP level and density functional theory employing BP86 with various basis sets. The equilibrium structures have linear geometries for Ae = Be and Mg but they are strongly bent for Ae = Sr and Ba while the calcium species has a quasi-linear structure with a very low bending potential. The calculated bond dissociation energies suggest a record-high BDE of De = 144.08 kcal mol-1 for OBeF- at the CCSD(T)/def2-TZVPP level, which is the strongest BDE for a dative bond that has been found so far. The BDE of the heavier homologues have a continuously decreasing order for Ae with Be > Mg (113.01 kcal mol-1) > Ca (84.06 kcal mol-1) > Sr (72.06 kcal mol-1) > Ba (60.00 kcal mol-1). The calculation of the charge distribution reveals a significant charge donation OAe ← F- with a declining sequence for the heavier atoms Ae. The oxygen atom in OAeF- carries always a higher partial charge than the fluorine atom, which contradicts the standard electronegativities of the atoms. The surprising partial charges are explained with the bonding situation of the atoms in the actual electronic structure. The bonding analysis of the OAe-F- bonds using the EDA-NOCV method shows that the bonds have much more electrostatic character than the Ae-F- bonds in the diatomic anions. This finding is supported by the results of the LED partitioning approach. The dative interactions have three major and one minor component. The assignment of a quadruple bond for the heavier species with Ae = Ca, Sr, Ba is not reasonable. The driving force for the bent geometries is the accumulation of electronic charge in the lone-pair region at the Ae atoms, which enhances the electrostatic attraction with the other atoms. An adequate description of the bonding situation is given by the formula O--Ae+ ← F-.

InnTl4-nH+ (n = 0∼4): Tetracoordinate Hydrogen in a Planar Fashion?

The recent report of planar tetracoordinate hydrogen (ptH) in In4H+ is very intriguing in planar hypercoordinate chemistry. Our high-level CCSD(T) calculations revealed that the proposed D4h-symmetric ptH In4H+ is a first-order saddle point with an imaginary frequency in the out-of-plane mode of the hydrogen atom. In fact, at the CCSD(T)/aug-cc-pV5Z/aug-cc-pV5Z-PP level, the C4v isomer, with the H atom located 0.70 Å above the In4 plane, is 0.5 kcal/mol more stable than the D4h isomer. However, given the small perturbation from planarity and essentially barrierless C4v ↔ D4h ↔ C4v transition, the vibrationally averaged structure can still be considered as a planar. Extending our exploration to the InnTl4-nH+ (n = 0-3) systems, we found all these ptH structures, except for In2Tl2H+, to be the putative global minimum. The single σ-delocalized interaction between the central hydrogen atom and InnTl4-n ligand rings proves pivotal in establishing planarity and aromaticity and conferring substantial stability upon these rule-breaking ptH species.

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.

Occupied-Virtual Orbitals for Chemical Valence with Applications to Charge Transfer in Energy Decomposition Analysis

In this article, we introduce the occupied-virtual orbitals for chemical valence (OVOCV). The OVOCVs can replace or complement the closely related idea of the natural orbitals for chemical valence (NOCV). The input is a difference density matrix connecting any initial single determinant to any final determinant, at a given molecular geometry, and a given one-particle basis. This arises in problems such as orbital rearrangement or charge transfer (CT) in energy decomposition analysis (EDA). The OVOCVs block-diagonalize the density difference operator into 2 × 2 blocks, which are spanned by one level that is filled in the initial state (the occupied OVOCV) and one that is empty (the virtual OVOCV). By contrast, the NOCVs fully diagonalize the density difference matrix and therefore are orbitals with mixed occupied-virtual character. Use of the OVOCVs makes it much easier to identify the donor and acceptor orbitals. We also introduce two different types of EDA methods with the OVOCVs and, most importantly, a charge decomposition analysis method that fixes the unreasonably large CT amount obtained directly from NOCV analysis. The square of the CT amount associated with each NOCV pair emerges as the appropriate value from the OVOCV analysis. When connecting the same initial and final states, this value is identical to the CT amount obtained from the independent absolutely localized molecular orbital (ALMO) complementary occupied-virtual orbital pair (COVP) analysis. The total, summed over all pairs, is also exactly the same as the independently suggested excitation number, as proved herein. Several examples are presented to compare NOCVs and OVOCVs: stretched H2+, a strong halogen bond between tetramethylthiourea and iodine, coordination of ethene in Zeise's salt, and binding in the Cp3La···C≡NCy complex.

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.

Revisiting the Structure and Bonding in Li5H6- and the Exploration of Reactivity: Planar Pentacoordinate Hydrogen

Recently, Guha and co-workers (Sarmah, K.; Kalita, A.; Purkayastha, S.; Guha, A. K. Pushing The Extreme of Multicentre Bonding: Planar Pentacoordinate Hydride. Angew. Chem. Int. Ed. 2024, e202318741) reported a highly intriguing bonding motif: planar pentacoordinate hydrogen (ppH) in Li5H6-, featuring C2v symmetry in the singlet state with two distinct H-Li (center-ring) bond distances. We herein revisited the potential energy surface of Li5H6- by using a target-oriented genetic algorithm. Our investigation revealed that the lowest-energy structure of Li5H6- exhibits a ppH configuration with very high D5h symmetry and a 1A1' electronic state. We did not find any electronic effect like Jahn-Teller distortion that could be responsible for lowering its symmetry. Moreover, our calculations demonstrated significant differences in the relative energies of other low-lying isomers. An energetically very competitive planar tetracoordinate hydrogen (ptH) isomer is also located, but it corresponds to a very shallow minimum on the potential energy surface depending on the used level of theory. Chemical bonding analyses, including AdNDP and EDA-NOCV, uncover that the optimal Lewis structure for Li5H6- involves H- ions stabilized by the Li5H5 crown. Surprisingly, despite the dominance of electrostatic interactions, the contribution from covalent bonding is also significant between ppH and the Li5H5 moiety, derived from H-(1s) → Li5H5 σ donation. Magnetically induced current density analysis revealed that due to minimal orbital overlap and the highly polar nature of the H-Li covalent interaction, the ppH exhibits local diatropic ring currents around the H centers, which fails to result in a global aromatic ring current. The coordination of Li5H6- with Lewis acids, BH3 and BMe3, instantly converts the ppH configuration to (quasi) ptH. These Lewis acid-bound ptH complexes show high electronic stability and high thermochemical stability against dissociation and, therefore, will be ideal candidates for the experimental realization.

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 using density functional theory are reported for the diatomic molecules AeN- (Ae=Ca, Sr, Ba). 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. 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 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 Ae←N- charge donation for most electronic states. In contrast, the NBO method predicts for all species medium to large Ae→N- electronic charge donation, which is due to the neglect of the (n)p AOs of Ae atoms as genuine valence orbitals. Neither the bond orders nor the bond lengths correlate with the bond dissociation energies. The EDA-NOCV calculations show that the heavier alkaline earth atoms Ca, Sr, Ba use their (n)s and (n-1)d orbitals for covalent bonding.

Exploring the Use of "Honorary Transition Metals" To Push the Boundaries of Planar Hypercoordinate Alkaline-Earth Metals

The quest for planar hypercoordinate atoms (phA) beyond six has predominantly focused on transition metals, with dodecacoordination being the highest reported thus far. Extending this bonding scenario to main-group elements, which typically lack d orbitals despite their larger atomic radius, has posed significant challenges. Intrigued by the potentiality of covalent bonding formation using the d orbitals of the heavier alkaline-earth metals (Ae = Ca, Sr, Ba), the so-called "honorary transition metals", we aim to push the boundaries of planar hypercoordination. By including rings formed by 12-15 atoms of boron-carbon and Ae centers, we propose a design scheme of 180 candidates with a phA. Further systematic screening, structural examination, and stability assessments identified 10 potential clusters with a planar hypercoordinate alkaline-earth metal (phAe) as the lowest-energy form. These unconventional structures embody planar dodeca-, trideca-, tetradeca-, and pentadecacoordinate atoms. Chemical bonding analyses reveal the important role of Ae d orbitals in facilitating covalent interactions between the central Ae atom and the surrounding boron-carbon rings, thereby establishing a new record for coordination numbers in the two-dimensional realm.

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.

Unraveling Reactivity Pathways: Dihydrogen Activation and Hydrogenation of Multiple Bonds by Pyramidalized Boron-Based Frustrated Lewis Pairs

The activation of H2 by pyramidalized boron-based frustrated Lewis Pairs (FLPs) (B/E-FLP systems where "E" refers to N, P, As, Sb, and Bi) have been explored using density functional theory (DFT) based computational study. The activation pathway for the entire process is accurately characterized through the utilization of the activation strain model (ASM) of reactivity, shedding light on the underlying physical factors governing the process. The study also explores the hydrogenation process of multiple bonds with the help of B/N-FLP. The research findings demonstrate that the liberation of activated dihydrogen occurs in a synchronized, albeit noticeably asynchronous, fashion. The transformation is extensively elucidated using the activation strain model and the energy decomposition analysis. This approach suggests a co-operative double hydrogen-transfer mechanism, where the B-H hydride triggers a nucleophilic attack on the carbon atom of the multiple bonds, succeeded by the migration of the protic N-H.

Stabilizing Monoatomic Two-Coordinate Bismuth(I) and Bismuth(II) Using a Redox Noninnocent Bis(germylene) Ligand

The formation of isolable monatomic BiI complexes and BiII radical species is challenging due to the pronounced reducing nature of metallic bismuth. Here, we report a convenient strategy to tame BiI and BiII atoms by taking advantage of the redox noninnocent character of a new chelating bis(germylene) ligand. The remarkably stable novel BiI cation complex 4, supported by the new bis(iminophosphonamido-germylene)xanthene ligand [(P)GeII(Xant)GeII(P)] 1, [(P)GeII(Xant)GeII(P) = Ph2P(NtBu)2GeII(Xant)GeII(NtBu)2PPh2, Xant = 9,9-dimethyl-xanthene-4,5-diyl], was synthesized by a two-electron reduction of the cationic BiIIII2 precursor complex 3 with cobaltocene (Cp2Co) in a molar ratio of 1:2. Notably, owing to the redox noninnocent character of the germylene moieties, the positive charge of BiI cation 4 migrates to one of the Ge atoms in the bis(germylene) ligand, giving rise to a germylium(germylene) BiI complex as suggested by DFT calculations and X-ray photoelectron spectroscopy (XPS). Likewise, migration of the positive charge of the BiIIII2 cation of 3 results in a bis(germylium)BiIIII2 complex. The delocalization of the positive charge in the ligand engenders a much higher stability of the BiI cation 4 in comparison to an isoelectronic two-coordinate Pb0 analogue (plumbylone; decomposition below -30 °C). Interestingly, 4[BArF] undergoes a reversible single-electron transfer (SET) reaction (oxidation) to afford the isolable BiII radical complex 5 in 5[BArF]2. According to electron paramagnetic resonance (EPR) spectroscopy, the unpaired electron predominantly resides at the BiII atom. Extending the redox reactivity of 4[OTf] employing AgOTf and MeOTf affords BiIII(OTf)2 complex 7 and BiIIIMe complex 8, respectively, demonstrating the high nucleophilic character of BiI cation 4.

Natural resonance-theoretic conceptions of extreme electronic delocalization in soft materials

In the broad context of Dalton's atomic hypothesis and subsequent classical vs. quantum understanding of macroscopic materials, we show how Pauling's resonance-type conceptions, as quantified in natural resonance theory (NRT) analysis of modern wavefunctions, can be modified to unify description of interatomic interactions from the Lewis-like limit of localized e-pair covalency in molecules to the extreme delocalized limit of supramolecular "soft matter" aggregation. Such "NRT-centric" integration of NRT bond orders for hard- and soft-matter interactions is illustrated with application to a long-predicted and recently synthesized organometallic sandwich-type complex ("diberyllocene") that exhibits bond orders ranging from the soft limit (bBeC ≈ 0.01) to the typical values (bCC ≈ 1.35) of molecular resonance-covalency in the organic domain, with intermediate value (bBeBe ≈ 0.86) for intermetallic Be⋯Be interaction.

The effect of {O,N}=X⋯M={Ti,Zr,Hf} interactions on the sensitivity of CNO2 trigger bonds in FOX-7: Approach based on the QTAIM/EDA-NOCV analysis

The local chemical reactivity of FOX-7 (1,1-diamino-2,2-nitroethylene, also known as DADNE from DiAminoDiNitroEthylene) was elucidated through a quantitative study of the electrostatic potential on the molecular surface, topological analysis based on Bader's theory, and the EDA-NOCV method. Unlike (O2N)2CC(NH2)H2N⋯Cp2MCH3+ complexes, which exhibit both σ-donor and π-acceptor features, the situation is different concerning the (H2N)2CC(NO2)(O)NO⋯Cp2MCH3+ complexes, where both charge transfers correspond to the σ-donation. The two charge transfers reinforce each other, resulting in increased stability for (H2N)2CC(NO2)(O)NO⋯Cp2MCH3+. This seems to strengthen the (H2N)2CC(NO2)(O)NO⋯M={Ti,Zr,Hf} bond, which may explain the high stability of (H2N)2CC(NO2)(O)NO⋯Cp2MCH3+ compared to (O2N)2CC(NH2)-H2N⋯Cp2MCH3+. Results from topological analysis revealed that the decreased sensitivity to decomposition of CNO2 bonds depends on the chemical nature of the interacting metal, and the best achievements are obtained for the Hf-based complex. Our results demonstrate that the interaction of M={Ti,Zr,Hf} with CNO2 is more favourable than that with CNH2, this specific action on the trigger bond may support the use of Metallocene Methyl Cations (MMC) as possible neutralisers.

Element effects in endohedral metal-metal-bonding fullerenes M2@C82 (M = Sc, Y, La, Lu)

Endohedral metal-metal-bonding fullerenes have recently emerged, in which encapsulated metals form a metal-metal bond. However, the physical reasons why some metal elements prefer to form metal-metal bonds inside fullerene are still unclear. Herein, we reported first-principles calculations on electronic structures, bonding properties, dynamics, and thermodynamic stabilities of endohedral metallofullerenes M2@C82 (M = Sc, Y, La, Lu). Multiple bonding analysis approaches unambiguously reveal the existence of one two-center two-electron σ covalent metal-metal bond in M2@C82 (M = Sc, Y, Lu); however, the La-La bonding interaction in La2@C82 is weaker and could not be categorized as one metal-metal covalent bond. The energy decomposition analysis on bonding interactions between an encapsulated metal dimer and fullerene cages suggested that there exist two electron-sharing bonds between a metal dimer and fullerene cages. The reasons why La2 prefers to donate electrons to fullerene cages rather than form a standard σ covalent metal-metal bond are mainly attributed to two following facts: La2 has a lower ionization potential, while the hybridization of ns, (n - 1)d, and np atomic orbitals in La2 is higher. Ab initio molecular dynamic simulations reveal that the M-M bond length at room temperature follows the trend of Sc < Lu < Y. The statistical thermodynamics calculations at different temperatures reveal that the experimentally observed endohedral metal-metal-bonding fullerenes M2@C82 have high concentrations in the endohedral fullerene formation temperature range.

Organometallic Allene [(μ-C)(Fe(CO)4 )2 ]: Bridging Carbon Showing Transformation from Classical Electron-Sharing Bonding to Double σ-Donor and Double π-Acceptor Ligation

Allenes (R2 C=C=CR2 ) have been traditionally perceived to feature localized orthogonal π-bonds between the carbon centres. We have carried out quantum-mechanical studies of the organometallic allenes envisioned by the isolobal replacement of the terminal CH2 groups by the d8 Fe(CO)4 fragment. Our studies have identified two organometallic allenes viz. D2d symmetric [(μ-C)(Fe(CO)4 )2 ] (2) and D3 symmetric [(μ-C)(Fe(CO)4 )2 ] (3) with trigonal bipyramidal coordination at the Fe atoms. Compound 2 features the bridging carbon atom in an equatorial position with respect to the ligands on the TM centre, while 3 features the central carbon atom in an axial position. The bis-pseudoallylic anionic delocalisation proposed in the C2-C1-C3 spine of organic allene is retained in the organometallic allene 2, and is transformed to a typical three-centre bis-allylic anionic delocalisation in the organometallic allene 3. The topological analysis of electron density also indicates a bis-allylic anionic type delocalisation in the organometallic allenes. The quantitative bonding analysis using the EDA-NOCV method suggests a transition from classical electron-sharing bonding between the central carbon atom and the terminal groups in 1 to donor-acceptor bonding in 3. Meanwhile, both electron-sharing and donor-acceptor bonding models are found to be probable heuristic bonding representations in the organometallic allene 2.

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.

Mimicking the C2 molecule: M2B2 and M3B2+ clusters (M = Li, Na) and the reactivity of the N-heterocyclic carbene bound Li2B2 complex

C2 has attracted considerable attention from the scientific community for its debatable bonding situation. Herein, we show that the global minima of M2B2 and M3B2+ (M = Li, Na) possess similar covalent bonding patterns to C2. Because of strong charge transfer from M2/M3 to B2 dimer, they can be better described as [M2]2+[B2]2- and [M3]3+[B2]2- salt complexes with the B22- core surrounded perpendicularly by two and three M+ atoms, respectively. The energy decomposition analyses in combination with the natural orbital for chemical valence theory give four bonding components in C2, M2B2, and M3B2+ clusters. However, the fourth component does not arise from a bonding interaction but from polarization/hybridization. Considering the effect of Pauli repulsion in σ-space, the attractive covalent interaction in these molecules mainly comes from the two π-bonds. We further presented stable N-heterocyclic carbene (NHC) and triphenylphosphine (PPh3) ligands bound Li2B2(NHC)2 and Li2B2(PPh3)2 complexes. A comparative study of reactivity towards L = CO2, CO, and N2 between Li2B2(NHC)2 and B2(NHC)2 is also performed. L-Li2B2(NHC)2 is highly stable against L dissociation at room temperature for L = CO2 and CO, and the stability is markedly higher than that in L-B2(NHC)2. The larger B2→L π-backdonation in L-Li2B2(NHC)2 also makes L more activated than in L-B2(NHC)2.

Comment on "The oxidation state in low-valent beryllium and magnesium compounds" by M. Gimferrer, S. Danés, E. Vos, C. B. Yildiz, I. Corral, A. Jana, P. Salvador and D. M. Andrada, Chem. Sci. 2022, , 6583

We challenge the assignment of the oxidation state +2 for beryllium and magnesium in the complexes Be(cAAC^Dip)₂ and Mg(cAAC^Dip)₂ as suggested by Gimferrer et al., Chem. Sci. 2022, 13, 6583 in a recent study. A careful review of the data in the ESI contradicts their own statement and shows that the results support the earlier suggestion that the metals are in the zero oxidation state. The authors reported wrong data for the excitation energies of Be and Mg to the ¹D (np²) state. We also correct some misleading statements about the EDA method.

Hitting the Bull's Eye: Stable HeBeOH+ Complex

It is now known that the heavier noble gases (Ng=Ar-Rn) show some varying degrees of reactivity with a gradual increase in reactivity along Ar-Rn. However, because of their very small size and very high ionization potential, helium and neon are the hardest targets to crack. Although few neon complexes are isolated at very low temperatures, helium needs very extreme situations like very high pressure. Here, we find that protonated BeO, BeOH⁺ can bind helium and neon spontaneously at room temperature. Therefore, extreme conditions like very low temperature and/or high pressure will not be required for their experimental isolation. The Ng-Be bond strength is very high for their heavier homologs and the bond strength shows a gradual increase from He to Rn. Moreover, the Ng-Be attractive energy is almost exclusively originated from the orbital interaction which is composed of one Ng(s/p_σ )→BeOH⁺ σ-donation and two weaker Ng(p_π )→BeOH⁺ π-donations, except for helium. Helium uses its low-lying vacant 2p orbitals to accept π-electron density from BeOH⁺ . Previously, such electron-accepting ability of helium was used to explain a somewhat stronger helium bond than neon for neutral complexes. However, the present results indicate that such π-back donations are too weak in nature to decide any energetic trend between helium and neon.

A Combined DFT, Energy Decomposition, and Data Analysis Approach to Investigate the Relationship Between Noncovalent Interactions and Selectivity in a Flexible DABCOnium/Chiral Anion Catalyst System

Developing strategies to study reactivity and selectivity in flexible catalyst systems has become an important topic of research. Herein, we report a combined experimental and computational study aimed at understanding the mechanistic role of an achiral DABCOnium cofactor in a regio- and enantiodivergent bromocyclization reaction. It was found that electron-deficient aryl substituents enable rigidified transition states via an anion-π interaction with the catalyst, which drives the selectivity of the reaction. In contrast, electron-rich aryl groups on the DABCOnium result in significantly more flexible transition states, where interactions between the catalyst and substrate are more important. An analysis of not only the lowest-energy transition state structures but also an ensemble of low-energy transition state conformers via energy decomposition analysis and machine learning was crucial to revealing the dominant noncovalent interactions responsible for observed changes in selectivity in this flexible system.

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.

Bare and ligand protected planar hexacoordinate silicon in SiSb3M3 + (M = Ca, Sr, Ba) clusters

The occurrence of planar hexacoordination is very rare in main group elements. We report here a class of clusters containing a planar hexacoordinate silicon (phSi) atom with the formula SiSb3M3 + (M = Ca, Sr, Ba), which have D 3h (1A1') symmetry in their global minimum structure. The unique ability of heavier alkaline-earth atoms to use their vacant d atomic orbitals in bonding effectively stabilizes the peripheral ring and is responsible for covalent interaction with the Si center. Although the interaction between Si and Sb is significantly stronger than the Si-M one, sizable stabilization energies (-27.4 to -35.4 kcal mol-1) also originated from the combined electrostatic and covalent attraction between Si and M centers. The lighter homologues, SiE3M3 + (E = N, P, As; M = Ca, Sr, Ba) clusters, also possess similar D 3h symmetric structures as the global minima. However, the repulsive electrostatic interaction between Si and M dominates over covalent attraction making the Si-M contacts repulsive in nature. Most interestingly, the planarity of the phSi core and the attractive nature of all the six contacts of phSi are maintained in N-heterocyclic carbene (NHC) and benzene (Bz) bound SiSb3M3(NHC)6 + and SiSb3M3(Bz)6 + (M = Ca, Sr, Ba) complexes. Therefore, bare and ligand-protected SiSb3M3 + clusters are suitable candidates for gas-phase detection and large-scale synthesis, respectively.