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.

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.

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.

Mono-Ortho-Beryllated Carbodiphosphoranes: Synthesis, Structure, Bonding and Reactivity.

The reaction of organoberyllium compounds with hexaphenylcarbodiphosphorane yields mono-ortho-beryllated complexes, which feature a double dative Be=C bond. The bonding situation in these compounds together with a simple carbodiphosphorane and an N-heterocyclic carbene adduct was analysed with energy decomposition analysis in combination with natural orbital for chemical valence as well as with quantum theory of atoms-in-molecules. Furthermore, the driving forces accountable for mono-ortho-beryllation were elucidated along with the reactivity of the Be=C bond.

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.

Analysis of the Unusual Chemical Bonds and Dipole Moments of AeF- (Ae=Be-Ba): A Lesson in Covalent Bonding.

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.

Li6 E5 Li6 : Tetrel Sandwich Complexes with 10-π-Electrons.

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  Li6 E5 Li6 sandwich structures, representing the first global minima featuring ten π-electrons E5 10- 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.

Transition Metal Behavior of Heavier Alkaline Earth Elements in Doped Monocyclic and Tubular Boron Clusters.

Quantum chemical calculations are carried out to design highly symmetric-doped boron clusters by employing the transition metal behavior of heavier alkaline earth (Ae = Ca, Sr, and Ba) metals. Following an electron counting rule, a set of monocyclic and tubular boron clusters capped by two heavier Ae metals were tested, which leads to the highly symmetric Ae2B8, Ae2B18, and Ae2B30 clusters as true minima on the potential energy surface having a monocyclic ring, two-ring tubular, and three-ring tubular boron motifs, respectively. Then, a thorough global minimum (GM) structural search reveals that a monocyclic B8 ring capped with two Ae atoms is indeed a GM for Ca2B8 and Ba2B8, while for Sr2B8 it is a low-lying isomer. Similarly, the present search also unambiguously shows the most stable isomers of Ae2B18 and Ae2B30 to be highly symmetric two- and three-ring tubular boron motifs, respectively, capped with two Ae atoms on each side of the tube. In these Ae-doped boron clusters, in addition to the electrostatic interactions, a substantial covalent interaction, specifically the bonding occurring between (n - 1)d orbitals of Ae and delocalized orbitals of boron motifs, provides the essential driving force behind their highly symmetrical structures and overall stability.

Stabilization of Cyclic C4 by Four Donor Ligands: A Theoretical Study of (L)4C4 (L = Carbene).

Quantum chemical studies using density functional theory were carried out for the (L)4C4 complexes with L = cAAC, DAC, NHC, SNHC, MIC1, and MIC2. The results show that the title complexes are highly stable with respect to dissociation, (L)4C4 → C4 + 4L. However, their stability with respect to (L)4C4 → 2(L)2C2 is crucial for the assessment of their experimental viability. The (L)4C4 complexes with L = cAAC and DAC dissociate exergonically at room temperature into two (L)2C2 units. In contrast, the other (L)4C4 complexes with L = NHC, SNHC, MIC1, and MIC2 are thermochemically stable with respect to dissociation, (L)4C4 → 2(L)2C2. The computed adiabatic ionization potentials of (L)4C4 complexes with L = NHC, MIC1, and MIC2 are lower than those for the cesium atom. Particularly, (MIC1)4C4 and (MIC2)4C4 will very easily lose electrons to form cationic complexes. The SNHC ligand is the best for the experimental realization of (L)4C4 complexes, followed by NHC. The bonding analysis using charge and energy decomposition methods suggests that the (L)3C4-CL bond can be best described as a typical electron-sharing double bond with a strong σ-bond and a weaker π-bond. Therefore, the core bonding pictures in the title complexes resemble a [4]radialene. Larger substituents at the carbene ligands enhance the stability of the complexes (L)4C4 against dissociation.

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.

Clusters and bulky Lewis acid protected complexes with planar hexacoordinate beryllium and magnesium.

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.

Planar pentacoordinate s-block metals.

The presence of a delocalized π-bond is often considered an essential criterion for achieving planar hypercoordination. Herein, we show that σ-delocalization could be sufficient to make the planar configuration the most stable isomer in a series of planar pentacoordinate s-block metals. High-level ab initio computations reveal that the global minimum of a series of interalkali and interalkali-alkaline earth clusters (LiNa5, Li5Mg+, Na5Mg+, K5Ca+, CaRb5+, Rb5Sr+, and SrCs5+) adopts a singlet D5h structure with a planar pentacoordinate lithium or alkaline earth metal (AE = Mg, Ca, Sr). These clusters are unusual combinations to stabilize a planar pentacoordinate atom, as all their constituents are electropositive. Despite the absence of π-electrons, Hückel's rule is fulfilled by the six σ-electrons. Furthermore, the systems exhibit a diatropic ring current in response to an external magnetic field and a strong magnetic shielding, so they might be classified as σ-aromatic. Therefore, multicenter σ-bonds and the resulting σ-delocalization stabilize these clusters, even though they lack π-aromaticity.

Structural Characterization and Bonding Analysis of [Hg{Fe(CO)5}2]2+ [SbF6]-2.

The non-classical carbonyl complex [Hg{Fe(CO)5}2]2+ [SbF6]-2 is prepared by reaction of Hg(SbF6)2 and excess Fe(CO)5 in anhydrous HF. The single-crystal X-ray structure reveals a linear Fe-Hg-Fe moiety as well as an eclipsed conformation of the eight basal CO ligands. Interestingly, the Hg-Fe bond length of 2.5745(7) Å is relatively similar to the corresponding Hg-Fe bonds in literature-known [Hg{Fe(CO)4}2]2- dianions (2.52-2.55 Å), which intrigued us to analyze the bonding situation in both the dications and dianions with the energy decomposition analysis with natural orbitals for chemical valence (EDA-NOCV) method. Both species are best described as Hg(0) compounds, which are also confirmed by the shape of the HOMO-4 and HOMO-5 of the dication and dianion, respectively, in which the electron pair is located mainly at the Hg. Furthermore, for the dication and the dianion, the σ back-donation from Hg into the [Fe(CO)5]22+ or the [Fe(CO)4]22- fragment is the most dominant orbital interaction and surprisingly these interaction energies are also very similar even in absolute values. The fact that both iron-based fragments are missing two electrons explains their prominent σ-acceptor properties.

B7 Be6 B7 : A Boron-Beryllium Sandwich Complex.

Planar boron clusters have often been regarded as "π-analogous" to aromatic arenes because of their similar delocalized π-bonding. However, unlike arenes such as C5 H5 - and C6 H6 , 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, B7 Be6 B7 . 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 B7 Be6 B7 can be attributed to strong electrostatic and covalent interactions between the fragments. Chemical bonding analysis shows that B7 Be6 B7 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.

Global Planar Tetra-, Penta- and Hexa-coordinate Silicon Clusters Constructed by Decorating SiO3 with Alkali Metals.

The achievement of the rule-breaking planar hypercoordinate motifs (carbon and other elements) is mainly attributed to a practical electronic stabilization mechanism, where the bonding of the central atom pz π electrons is a crucial issue. We have demonstrated that strong multiple bonds between the central atom and partial ligands can be an effective approach to explore stable planar hypercoordinate species. A set of planar tetra-, penta- and hexa-coordinate silicon clusters were herein found to be the lowest-energy structure, which can be viewed as decorating SiO3 by alkali metals in the MSiO3 - , M2 SiO3 and M3 SiO3 + (M=Li, Na) clusters. The strong charge transfer from M atoms to SiO3 effectively results in [M]+ SiO3 2- , [M2 ]2+ SiO3 2- and [M3 ]3+ SiO3 2- salt complexes, where the Si-O multiple bonding and structural integrity of the Benz-like SiO3 framework is maintained better than the corresponding SiO3 2- motifs. The bonding between M atoms and SiO3 motif is best described as M+ forming a few dative interactions by employing its vacant s, p, and high-lying d orbitals. These considerable M←SiO3 interactions and Si-O multiple bonding give rise to the highly stable planar hypercoordinate silicon clusters.

Quest of Quadruple Bonding Between Two Main-Group Atoms in AeB- and AeC (Ae=Ca, Sr, Ba) and the Role of d Orbitals of Heavier Alkaline-Earth Atoms in Covalent Interactions.

Quantum chemical calculations using ab initio methods at the MRCI+Q(6,8)/def2-QZVPP and CCSD(T)/def2-QZVPP levels as well as density functional theory are reported for the diatomic molecules AeB- and isoelectronic AeC (Ae=Ca, Sr, Ba). The boride anions AeB- have an electronic triplet (3 Σ- ) ground state. The quintet (5 Σ- ) state is 5.8-12.3 kcal/mol higher in energy and the singlet (1 Δ) state is 13.1-15.3 kcal/mol above the triplet. The isoelectronic AeC molecules are also predicted to have a low-lying triplet (3 Σ- ) state but the quintet (5 Σ- ) state is only 2.2 kcal/mol (SrC) and 2.9 kcal/mol (CaC) above the triplet state. The triplet (3 Σ- ) and quintet (5 Σ- ) states of BaC are nearly isoenergetic. All systems have rather strong bonds. The calculated bond dissociation energies of the triplet (3 Σ- ) state are between De =38.3-41.7 kcal/mol for AeB- and De =49.4-57.5 kcal/mol for AeC. The barium species have always the strongest bonds whereas the calcium and strontium compounds have similar BDEs. The bonding analysis indicates that there is little charge migration in AeB- in the direction Ae→B- where the alkaline earth atoms carry positive charges between 0.09 e-0.22 e. The positive charges at the Ae atoms are much larger in AeC where the charge migration Ae→C is between 0.90 e-0.91 e. A detailed analysis of the interatomic interactions with the EDA-NOCV method shows that all diatomic species AeB- and AeC are built from dative interactions between Ae (1 S, ns2 ) and B- or C (3 P, 2 s2 2pπ 1 2pπ' 1 ). The eventually formed bonds in AeC are better described in terms of interactions between the ions Ae+ (2 S, ns1 )+C- (4 S, 2 s2 2pπ 1 2pπ' 1 2pσ 1 ). Inspection of the orbital interactions suggests that the alkaline earth atoms Ca, Sr, Ba use mainly their (n-1)d AOs besides the (n)s AOs for the covalent bonds. This creates a second energetically low-lying σ-bonding MO in the molecules, which feature valence orbitals with the order ϕ1 (σ-bonding)<ϕ2 (σ-bonding)<ϕ3 (degenerate π-bonding). All four occupied valence MOs of AeB- and AeC are bonding orbitals. Since the degenerate π orbitals ϕ3 are only singly occupied, the formal bond order is three.

Synthesis of a rhodium(III) dinitrogen complex using a calix[4]arene-based diphosphine ligand.

The synthesis and characterisation of the rhodium(III) dinitrogen complex [Rh(2,2'-biphenyl)(CxP2)(N2)]+ are described, where CxP2 is a trans-spanning calix[4]arene-based diphosphine and the dinitrogen ligand is projected into the cavity of the macrocycle.

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, 13, 6583.

We challenge the assignment of the oxidation state +2 for beryllium and magnesium in the complexes Be(cAACDip)2 and Mg(cAACDip)2 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 1D (np2) state. We also correct some misleading statements about the EDA method.

BH4 Ng+ (Ar-Rn): Viable Compounds with a B-Ng Covalent Bond.

In this work, we explore, using high-level calculations, the ability of BH4 + to interact with noble gases. The He system is energetically unstable, while the Ne system could only be observed at cryogenic temperatures. In the case of the Ar, Kr and Xe systems, all are energetically stable, even at room temperature. The different chemical bond descriptors reveal a covalent character between B and the noble gas from Ar to Rn. However, this interaction gradually weakens the multicentric bond between the boron atom and the H2 fragment. Thus, although BH4 Rn+ exhibits a strong covalent bond, it tends to dissociate at room temperature into BH2 Rn+ +H2 .

Bonding situations in tricoordinated beryllium phenyl complexes.

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