Search

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. [ABSTRACT FROM AUTHOR]

Planar hexacoordination (ph) is only rarely reported in the literature. So far, only a few neutral and cationic molecules possessing phE (E = C, Si, B, Al, Ga) in the most stable isomer are predicted theoretically. Present electronic structure calculations report hitherto unknown anionic planar hexcoordinate beryllium and magnesium, phBe/Mg, as the most stable isomer. Global minimum searches show that the lowest energy structure of BeC6M3− (M = Al, Ga) and MgC6M3− (M = Ga, In, Tl) is the D3h symmetric phBe/Mg clusters, where beryllium/magnesium is covalently bonded with six carbon centers and M is located in a bridging position between two carbon centers. These global minimum phBe/Mg clusters are highly kinetically stable against isomerization, facilitating the experimental confirmation by photoelectron spectroscopy. Noteworthy is the fact that the phBe/Mg center is linked with carbon centers through three 7c–2e delocalized σ bonds and three 7c–2e π bonds, making the cluster double aromatic (σ + π) in nature. The bonding between the Be/Mg and outer ring moiety can be best expressed as an electron-sharing σ-bond between the s orbital of Be+/Mg+ and C6M32− followed by three dative interactions involving empty pπ and two in-plane p orbitals of Be/Mg. Furthermore, Lewis basic M centers of the title clusters can be passivated through the complexation with bulky Lewis acid, 9-boratriptycene, lowering the overall reactivity of the cluster, which can eventually open up the possibility of their large-scale syntheses. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Quantum chemical calculations of the anions AeF− (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 detailed bonding analyses using different charge‐ and energy partitioning methods show that the molecules possess three distinctively different dative bonds in the lighter species with Ae=Be, Mg and four dative bonds when Ae=Ca, Sr, Ba. The occupied 2p atomic orbitals (AOs) and to a lesser degree the occupied 2s AO of F− donate electronic charge into the vacant spx(σ) and p(π) orbitals of Be and Mg which leads to a triple bond Ae F−. The heavier Ae atoms Ca, Sr, Ba use their vacant (n‐1)d AOs as acceptor orbitals which enables them to form a second σ donor bond with F− that leads to quadruply bonded Ae F− (Ae=Ca−Ba). The presentation of molecular orbitals or charge distribution using only one isodensity value may give misleading information about the overall nature of the orbital or charge distribution. Better insights are given by contour line diagrams. The ELF calculations provide monosynaptic and disynaptic basins of AeF− which nicely agree with the analysis of the occupied molecular orbitals and with the charge density difference maps. A particular feature of the covalent bonds in AeF− concerns the inductive interaction of F− with the soft valence electrons in the (n)s valence orbitals of Ae. The polarization of the (n)s2 electrons induces a (n)spx hybridized lone‐pair orbital at atom Ae, which yields a large dipole moment with the negative end at Ae. The concomitant formation of a vacant (n)spx AO of atom Ae, which overlaps with the occupied 2p(σ) AO of F−, leads to a strong covalent σ bond. [ABSTRACT FROM AUTHOR]

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 Energy Decomposition Analysis–Natural Orbitals of Chemical Valence (EDA-NOCV) method, and the results are critically discussed and compared with data from Quantum Theory of Atoms in Molecules, Natural Bond Orbital, 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 valence bond 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 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. The bond formation in LiF is dominated by orbital interactions (90%) between Li and F yielding a single bond, but the eventually formed bond comes mainly from the electrostatic attraction between Li+ and F−, where the minor orbital interactions (10%) have equally strong σ and π components. The symmetry allowed bond formation of BeO between Be in the 1S ground state and O in the excited 1D state is dominated (90%) by a strong dative Be → O σ bond with negligible π interactions. The final bond situation in BeO is best described by the interaction between Be+ and O−, where the Coulombic forces provide 60% of the attraction and the orbital interactions give equally strong σ and π bonds. The chemical bond in BN is analyzed in the X3Π ground state and the a1Σ+ excited state. Both states have triple bonds with strong π bonds, which are in the a1Σ+ state even stronger than the σ bond. [ABSTRACT FROM AUTHOR]

When (4n +2) π‐electrons are located in single planar ring, it conventionally qualifies as aromatic. According Hückel's rule, systems possessing ten π‐electrons should be aromatic. Herein we report a series of D5h Li6E5Li6 sandwich structures, representing the first global minima featuring ten π‐electrons E510− ring (E=Si−Pb). However, these π‐electrons localize as five π‐lone‐pairs rather than delocalized orbitals. The high symmetry structure achieved is a direct consequence of σ‐aromaticity, particularly favored in elements from Si to Pb, resulting in a pronounced diatropic ring current flow that contributes to the enhanced stability of these systems. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Quantum chemical calculations have been carried out at the BP86/def2‐SVP level on Ge8(Sit‐butyl2methyl)6 (1) and the bonding situation has been analyzed with a variety of methods. The calculated equilibrium geometry of 1 is in good agreement with the reported x‐ray structure analysis. The D3 correction for dispersion interactions as a sum of pairwise attractions leads to an overestimate of the effect of dispersion forces. Calculations at BP86‐D3(BJ)/def2‐SVP give shorter bonds for Ge(I)−Ge(I) than for Ge(0)−Ge(I), which is in contrast to the experimental values and the BP86/def2‐SVP results. The NBO analysis suggests that the best Lewis structure of 1 has lone‐pair orbitals at the Ge(0) atoms with occupation numbers of 1.70 e. A lone‐pair character at Ge(0) albeit with less weight is also suggested by the shape of the HOMO, which is an antibonding orbital between the Ge(0) atoms with small contributions from the Ge(I) atoms. The LUMO of 1 is the corresponding bonding combination of the Ge(0) AOs, which can be explained with the reluctance of the heavier main‐group atoms to s/p hybridization of the valence orbitals. The calculated bond order values suggest significant direct Ge(0)−Ge(0) interactions. This is supported by the shape of the HOMO and by the results of EDA‐NOCV calculations. The deformation densities and the orbitals associated with the pairwise orbital interaction show that there is a direct charge flow between the Ge(0) atoms of the two fragments, but it is not completely separated from the Ge(0)−Ge(I) and Ge(I)−Ge(I) bond formation. The QTAIM calculations suggest that 1 has a cubic structure with a cage critical point but not a bond critical point for the Ge(0)−Ge(0) interactions. The dispersion interactions of the large substituents in 1 have a significant influence on the stability of the compound. [ABSTRACT FROM AUTHOR]

Planar boron clusters have often been regarded as "π‐analogous" to aromatic arenes because of their similar delocalized π‐bonding. However, unlike arenes such as C5H5− and C6H6, boron clusters have not previously shown the ability to form sandwich complexes. In this study, we present the first sandwich complex involving beryllium and boron, B7Be6B7. The global minimum of this combination adopts a unique architecture having a D6h geometry, featuring an unprecedented monocyclic Be6 ring sandwiched between two quasi‐planar B7 motifs. The thermochemical and kinetic stability of B7Be6B7 can be attributed to strong electrostatic and covalent interactions between the fragments. Chemical bonding analysis shows that B7Be6B7 can be considered as a [B7]3−[Be6]6+[B7]3− complex. Moreover, there is a significant electron delocalization within this cluster, supported by the local diatropic contributions of the B7 and Be6 fragments. [ABSTRACT FROM AUTHOR]

Quantum chemical calculations have been carried out at the BP86/def2‐SVP level on Ge8(Sit‐butyl2methyl)6 (1) and the bonding situation has been analyzed with a variety of methods. The calculated equilibrium geometry of 1 is in good agreement with the reported x‐ray structure analysis. The D3 correction for dispersion interactions as a sum of pairwise attractions leads to an overestimate of the effect of dispersion forces. Calculations at BP86‐D3(BJ)/def2‐SVP give shorter bonds for Ge(I)−Ge(I) than for Ge(0)−Ge(I), which is in contrast to the experimental values and the BP86/def2‐SVP results. The NBO analysis suggests that the best Lewis structure of 1 has lone‐pair orbitals at the Ge(0) atoms with occupation numbers of 1.70 e. A lone‐pair character at Ge(0) albeit with less weight is also suggested by the shape of the HOMO, which is an antibonding orbital between the Ge(0) atoms with small contributions from the Ge(I) atoms. The LUMO of 1 is the corresponding bonding combination of the Ge(0) AOs, which can be explained with the reluctance of the heavier main‐group atoms to s/p hybridization of the valence orbitals. The calculated bond order values suggest significant direct Ge(0)−Ge(0) interactions. This is supported by the shape of the HOMO and by the results of EDA‐NOCV calculations. The deformation densities and the orbitals associated with the pairwise orbital interaction show that there is a direct charge flow between the Ge(0) atoms of the two fragments, but it is not completely separated from the Ge(0)−Ge(I) and Ge(I)−Ge(I) bond formation. The QTAIM calculations suggest that 1 has a cubic structure with a cage critical point but not a bond critical point for the Ge(0)−Ge(0) interactions. The dispersion interactions of the large substituents in 1 have a significant influence on the stability of the compound. [ABSTRACT FROM AUTHOR]